Changeset 3268 for trunk/MagicSoft/Mars/mcalib/MExtractBlindPixel.cc

Timestamp:

02/23/04 22:17:10 (21 years ago)

Author:

gaug

Message:

*** empty log message ***

File:

• trunk/MagicSoft/Mars/mcalib/MExtractBlindPixel.cc (modified) (4 diffs)

trunk/MagicSoft/Mars/mcalib/MExtractBlindPixel.cc

@@ -40 +40 @@
 #include "MRawEvtPixelIter.h"
 
+#include "MPedestalCam.h"
+#include "MPedestalPix.h"
+
 #include "MExtractedSignalBlindPixel.h"
 
@@ -49 +52 @@
 const Byte_t MExtractBlindPixel::fgSaturationLimit = 254;
 const Byte_t MExtractBlindPixel::fgFirst =  3;
-const Byte_t MExtractBlindPixel::fgLast  = 14;
+const Byte_t MExtractBlindPixel::fgLast  = 16;
 
 // --------------------------------------------------------------------------
@@ -98 +101 @@
     }
 
+    fPedestals = (MPedestalCam*)pList->FindObject(AddSerialNumber("MPedestalCam"));
+    if (!fPedestals)
+    {
+        *fLog << err << AddSerialNumber("MPedestalCam") << " not found... aborting." << endl;
+        return kFALSE;
+    }
 
     fBlindPixel = (MExtractedSignalBlindPixel*)pList->FindCreateObj(AddSerialNumber("MExtractedSignalBlindPixel"));
@@ -177 +186 @@
       *fLog << warn << "WARNING - saturation occurred in the Blind Pixel " << endl;
 
+    //
+    // catch the pointer to the event, because MHCalibrationChargeBlindPix will need it later
+    //
+    fPixelIter = &pixel;
+
     fBlindPixel->SetReadyToSave();
 
     return kTRUE;
+}
+
+Int_t MExtractBlindPixel::PostProcess()
+{
+
+    MPedestalPix &pedpix  = (*fPedestals)[fBlindPixelIdx];    
+
+    if (&pedpix)
+    {
+         fPed       = pedpix.GetPedestal()   * fNumSamples;
+         fPedErr    = pedpix.GetPedestalRms()* fNumSamples / TMath::Sqrt((Float_t)fPedestals->GetTotalEntries());
+         fPedRms    = pedpix.GetPedestalRms()* TMath::Sqrt((Float_t)fNumSamples);
+         fPedRmsErr = fPedErr/2.;
+    }
+
+    return kTRUE;
+
 }