Changeset 4781 for trunk/MagicSoft/Mars/mcalib

Timestamp:

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

Author:

gaug

Message:

*** empty log message ***

File:

• trunk/MagicSoft/Mars/mcalib/MCalibrationChargeCalc.cc (modified) (2 diffs)

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);
     }