Changeset 1092 for trunk

Timestamp:

11/28/01 10:15:33 (23 years ago)

Author:

blanch

Message:

Change in the methode to substract the pedestal. The new
methode is valid either for MC data and real data but only for
the current integration (GetSumHiGainSamples()).

File:

• trunk/MagicSoft/Mars/manalysis/MCerPhotCalc.cc (modified) (7 diffs)

trunk/MagicSoft/Mars/manalysis/MCerPhotCalc.cc

@@ -45 +45 @@
 #include "MLogManip.h"
 
+#include "MRawRunHeader.h"
 #include "MRawEvtData.h"       // MRawEvtData::GetNumPixels
 #include "MCerPhotEvt.h"
@@ -66 +67 @@
     AddToBranchList("fHiGainFadcSamples");
     AddToBranchList("fLoGainFadcSamples");
+
 }
 
@@ -71 +73 @@
 //
 // The PreProcess searches for the following input containers:
+//  - MRawRunHeader
 //  - MRawEvtData
 //  - MPedestalCam
@@ -80 +83 @@
 Bool_t MCerPhotCalc::PreProcess(MParList *pList)
 {
+    MRawRunHeader *RunHeader = (MRawRunHeader*)pList->FindObject("MRawRunHeader");
+    if (!RunHeader)
+    {
+        *fLog << dbginf << "MRawRunHeader not found... aborting." << endl;
+        return kFALSE;
+    }
+    fNumLoGainFadcSlices= RunHeader->GetNumSamplesLoGain();
+    fNumHiGainFadcSlices= RunHeader->GetNumSamplesHiGain();
+
     fRawEvt = (MRawEvtData*)pList->FindObject("MRawEvtData");
     if (!fRawEvt)
@@ -103 +115 @@
 // --------------------------------------------------------------------------
 //
-// Calculate the integral of the FADC time slaices and store them as a new
+// Calculate the integral of the FADC time slices and store them as a new
 // pixel in the MCerPhotEvt container.
 //
@@ -111 +123 @@
 
     MRawEvtPixelIter pixel(fRawEvt);
+
     while (pixel.Next())
     {
+
         const UInt_t pixid = pixel.GetPixelId();
 
         const MPedestalPix &ped = (*fPedestals)[pixid];
 
-        const Float_t nphot = (Float_t)pixel.GetSumHiGainSamples() - ped.GetMean();
+        Float_t nphot = (Float_t)pixel.GetSumHiGainSamples();
 
-        fCerPhotEvt->AddPixel(pixid, nphot, ped.GetMeanRms());
+        if(nphot>0) nphot -= fNumHiGainFadcSlices*ped.GetMean();
+
+        fCerPhotEvt->AddPixel(pixid, nphot, sqrt(fNumHiGainFadcSlices)*ped.GetMeanRms());
 
         // FIXME! Handling of Lo Gains is missing!
@@ -128 +144 @@
     return kTRUE;
 }
-