Changeset 3757 for trunk/MagicSoft/Mars/mcalib

Timestamp:

04/15/04 17:13:25 (21 years ago)

Author:

gaug

Message:

*** empty log message ***

File:

• trunk/MagicSoft/Mars/mcalib/MCalibrateRelTimes.cc (modified) (1 diff)

trunk/MagicSoft/Mars/mcalib/MCalibrateRelTimes.cc

@@ -174 +174 @@
       
       Float_t signal;
+      Float_t sigerr;
             
       if (sig.IsLoGainUsed())
-        signal = sig.GetArrivalTimeLoGain();
+        {
+          signal = sig.GetArrivalTimeLoGain();
+          sigerr = sig.GetArrivalTimeLoGainError();
+        }
       else
-        signal = sig.GetArrivalTimeHiGain();
+        {
+          signal = sig.GetArrivalTimeHiGain();
+          sigerr = sig.GetArrivalTimeHiGainError();
+        }
       
       const Float_t time = signal - offset;
-
-      //
-      // The following part is the outcommented first version of the error calculation
-      // Contact Markus Gaug for questions (or wait for the next documentation update...)
-      //
-      /*
-        nphotErr = signal    > 0 ? signalErr*signalErr / (signal * signal)  : 0.
-                 + calibConv > 0 ? calibConvVar  / (calibConv * calibConv ) : 0.
-                 + calibQE   > 0 ? calibQEVar    / (calibQE   * calibQE   ) : 0.;
-        nphotErr  = TMath::Sqrt(nphotErr) * nphot;
-      */
+      Float_t err        = sigerr*sigerr + precision*precision;
+      if (err > 0)
+        err = TMath::Sqrt(err);
 
       fArrivalTime->SetTime(pixidx,time);
+      fArrivalTime->SetTimeErr(pixidx,err);
 
     } /* for (UInt_t pixidx=0; pixidx<npix; pixidx++) */