Changeset 3195 for trunk/MagicSoft/Mars/mcalib/MExtractPINDiode.cc

Timestamp:

02/16/04 14:20:26 (21 years ago)

Author:

gaug

Message:

*** empty log message ***

File:

• trunk/MagicSoft/Mars/mcalib/MExtractPINDiode.cc (modified) (3 diffs)

trunk/MagicSoft/Mars/mcalib/MExtractPINDiode.cc

@@ -117 +117 @@
     }
 
+    const MPedestalPix &ped   = (*fPedestals)[fPINDiodeIdx]; 
+
+    if (&ped)
+      {
+        fPedestal = ped.GetPedestal();
+        fPedRms   = ped.GetPedestalRms();
+      }
+    else
+      {
+        *fLog << err << " Cannot find MPedestalPix of the PIN Diode (idx=" 
+              << fPINDiodeIdx << ")" << endl;
+        return kFALSE;
+      }
+
     return kTRUE;
 }
@@ -163 +177 @@
                sum, sum2, sat, max);
     
-    const MPedestalPix &ped   = (*fPedestals)[fPINDiodeIdx]; 
-
-    Float_t pedes  = -999.;
-    Float_t pedrms = -999.;
-
-    if (&ped)
-      {
-        pedes  = ped.GetPedestal();
-        pedrms = ped.GetPedestalRms();
-      }
-    else
-      {
-        *fLog << warn << " Cannot find MPedestalPix of the PIN Diode (idx=" 
-              << fPINDiodeIdx << ")" << endl;
-
-      }
 
     const Float_t var = ((Float_t)sum2 - (Float_t)sum*sum/fNumSamples)/(fNumSamples-1);
@@ -186 +184 @@
     // FIXME: The following formulae have to be revised!!
     //
-    fPINDiode->SetExtractedSignal(sum - pedes*fNumSamples, pedrms*fSqrtSamples);
+    fPINDiode->SetExtractedSignal(sum - fPedestal*fNumSamples, fPedRms*fSqrtSamples);
     fPINDiode->SetExtractedRms   (rms, rms/2./fSqrtSamples);
     fPINDiode->SetExtractedTime  (max, rms/fSqrtSamples);