Changeset 4781

Timestamp:

08/27/04 20:54:37 (20 years ago)

Author:

gaug

Message:

*** empty log message ***

Location:

trunk/MagicSoft/Mars

Files:

• Changelog (modified) (1 diff)
• mcalib/MCalibrationChargeCalc.cc (modified) (2 diffs)

trunk/MagicSoft/Mars/Changelog

@@ -34 +34 @@
    * mjobs/MJPedestal.h
      - added function SetUseData() and IsUseData()
+
+   * mcalib/MCalibrationQEPix.cc
+     - changed initialization of fNorm... variables from 1. to -1. 
+       Like this, an accidental not-filling of MCalibrationQECam will
+       be recognized easier.
 
    * mcalib/MCalibrationChargeCalc.cc

trunk/MagicSoft/Mars/mcalib/MCalibrationChargeCalc.cc

@@ -789 +789 @@
   const Float_t pedes  = ped.GetPedestal();
   const Float_t prms   = ped.GetPedestalRms();
-  const Float_t num    = TMath::Sqrt((Float_t)fPedestals->GetTotalEntries());
-
+  const Int_t   num    = fPedestals->GetTotalEntries();
+  
+  //
+  // RMS error set by PedCalcFromLoGain, 0 in case MPedCalcPedRun was used.
+  //
+  const Float_t prmserr = num>0 ? prms/TMath::Sqrt(2.*num) : ped.GetPedestalRmsError();
 
   //
@@ -797 +801 @@
   if (cal.IsHiGainSaturation())
     {
-      cal.SetPedestal(pedes* fNumLoGainSamples,
-                      prms * fSqrtLoGainSamples,
-                      prms * fNumLoGainSamples / num);
+      cal.SetPedestal(pedes   * fNumLoGainSamples,
+                      prms    * fSqrtLoGainSamples,
+                      prmserr * fSqrtLoGainSamples);
       cal.CalcLoGainPedestal((Float_t)fNumLoGainSamples, aidx);
     }
   else
     {
-      cal.SetPedestal(pedes* fNumHiGainSamples,
-                      prms * fSqrtHiGainSamples,
-                      prms * fNumHiGainSamples / num);
+      cal.SetPedestal(pedes   * fNumHiGainSamples,
+                      prms    * fSqrtHiGainSamples,
+                      prmserr * fSqrtHiGainSamples);
     }