Changeset 5194 for trunk/MagicSoft/Mars

Timestamp:

10/06/04 11:44:00 (20 years ago)

Author:

gaug

Message:

*** empty log message ***

Location:

trunk/MagicSoft/Mars

Files:

• Changelog (modified) (1 diff)
• mcalib/MCalibrationChargeCalc.cc (modified) (1 diff)
• mcalib/MCalibrationChargePix.cc (modified) (4 diffs)
• mcalib/MCalibrationChargePix.h (modified) (3 diffs)

trunk/MagicSoft/Mars/Changelog

@@ -22 +22 @@
  2004/10/05: Markus Gaug
    
+  * mcalib/MCalibrationChargePix.[h,cc]
+  * mcalib/MCalibrationChargeCalc.cc
+    - after talking with Florian, we reduced the electronic noise for the 
+      low-gain pedestal calculation from 1.5 to 0.8, independent from
+      inner or outer pixels. 
+
   * mcalib/MCalibrationChargeCalc.cc
     - fixed small bug for the case that multiple colours are found, but 

trunk/MagicSoft/Mars/mcalib/MCalibrationChargeCalc.cc

@@ -910 +910 @@
                       prms    * fSqrtLoGainSamples,
                       prmserr * fSqrtLoGainSamples);
-      cal.CalcLoGainPedestal((Float_t)fNumLoGainSamples, aidx);
+      cal.CalcLoGainPedestal((Float_t)fNumLoGainSamples);
     }
   else

trunk/MagicSoft/Mars/mcalib/MCalibrationChargePix.cc

@@ -81 +81 @@
 using namespace std;
 
-const Float_t MCalibrationChargePix::gkElectronicPedRmsInner    = 1.5;
-const Float_t MCalibrationChargePix::gkElectronicPedRmsOuter    = 1.8;
-const Float_t MCalibrationChargePix::gkElectronicPedRmsErr      = 0.35;
-const Float_t MCalibrationChargePix::gkFFactor                  = 1.15;
-const Float_t MCalibrationChargePix::gkFFactorErr               = 0.02;
+const Float_t MCalibrationChargePix::gkElectronicPedRms    = 0.8;
+const Float_t MCalibrationChargePix::gkElectronicPedRmsErr = 0.1;
+const Float_t MCalibrationChargePix::gkFFactor             = 1.15;
+const Float_t MCalibrationChargePix::gkFFactorErr          = 0.02;
 
 const Float_t MCalibrationChargePix::fgConversionHiLo           = 10.;
@@ -244 +243 @@
 //
 // Get the relative variance of the electronics pedestal RMS
-// - If aidx is 0, return rel. variance of gkElectronicPedRmsInner
-// - If aidx is 1, return rel. variance of gkElectronicPedRmsOuter
-//
-const Float_t MCalibrationChargePix::GetElectronicPedRmsRelVar(const Int_t aidx) const 
-{
-
-  if (aidx == 0)
-    return gkElectronicPedRmsErr * gkElectronicPedRmsErr / gkElectronicPedRmsInner / gkElectronicPedRmsInner;
-
-  if (aidx == 1)
-    return gkElectronicPedRmsErr * gkElectronicPedRmsErr / gkElectronicPedRmsOuter / gkElectronicPedRmsOuter;
-
-  return -1.;
+// - returns rel. variance of gkElectronicPedRms
+//
+const Float_t MCalibrationChargePix::GetElectronicPedRmsRelVar() const 
+{
+  return gkElectronicPedRmsErr * gkElectronicPedRmsErr / gkElectronicPedRms / gkElectronicPedRms;
 }
  
@@ -1006 +997 @@
 //   Var (Low Gain Ped RMS Square) = Var(LoGain NSB square) + Var(elec.ped. square)
 //
-void MCalibrationChargePix::CalcLoGainPedestal(Float_t logainsamples, const Int_t aidx)
+void MCalibrationChargePix::CalcLoGainPedestal(Float_t logainsamples)
 {
 
@@ -1015 +1006 @@
     return;
 
-  const Float_t elecPedRms     = (aidx == 0 ? gkElectronicPedRmsInner : gkElectronicPedRmsOuter ) 
-                               * TMath::Sqrt(logainsamples) / fConversionHiLo;
-  const Float_t elecPedRmsVar  = ( GetElectronicPedRmsRelVar(aidx) + GetConversionHiLoRelVar() )
+  const Float_t elecPedRms     = gkElectronicPedRms * TMath::Sqrt(logainsamples) / fConversionHiLo;
+  const Float_t elecPedRmsVar  = ( GetElectronicPedRmsRelVar() + GetConversionHiLoRelVar() )
                                * elecPedRms * elecPedRms;
   

trunk/MagicSoft/Mars/mcalib/MCalibrationChargePix.h

@@ -10 +10 @@
 private:
 
-  static const Float_t gkElectronicPedRmsInner;  //! Electronic component of ped. RMS (now set to: 1.5)
-  static const Float_t gkElectronicPedRmsOuter;  //! Electronic component of ped. RMS (now set to: 1.8)
-  static const Float_t gkElectronicPedRmsErr;    //! Error Electr. component ped. RMS (now set to: 0.35)
-  static const Float_t gkFFactor;                //! Laboratory F-factor PMTs         (now set to: 1.15)
-  static const Float_t gkFFactorErr;             //! Laboratory F-factor Error PMTs   (now set to: 0.02)
+  static const Float_t gkElectronicPedRms;     //! Electronic component of ped. RMS from opt. receiver on (now set to: 0.8)
+  static const Float_t gkElectronicPedRmsErr;  //! Error Electr. component ped. RMS (now set to: 0.1)
+  static const Float_t gkFFactor;              //! Laboratory F-factor PMTs         (now set to: 1.15)
+  static const Float_t gkFFactorErr;           //! Laboratory F-factor Error PMTs   (now set to: 0.02)
 
   static const Float_t fgConversionHiLo;         //! Default fConversionHiLo          (now set to: 10.)
@@ -45 +44 @@
   enum  { kFFactorMethodValid   };          // Possible bits to be set
 
-  const Float_t GetElectronicPedRmsRelVar(const Int_t aidx) const;
+  const Float_t GetElectronicPedRmsRelVar() const;
   const Float_t GetConversionHiLoRelVar()  const;
   const Float_t GetFFactorRelVar()     const;
@@ -108 +107 @@
 
   // Calculations
-  void   CalcLoGainPedestal       ( const Float_t logainsamples, const Int_t aidx );
+  void   CalcLoGainPedestal       ( const Float_t logainsamples);
   Bool_t CalcReducedSigma  (  );
   Bool_t CalcFFactor       (  );