Changeset 1127 for trunk/MagicSoft

Timestamp:

12/18/01 17:51:08 (23 years ago)

Author:

blanch

Message:

Remove fNumLoGainFadcSlices and fNumHiGainFadcSlices members and it has been
fixed a bug in the convention of Pedestals and Pedestals fluctuations (it did
not affect previous result since the bug was done twice).

File:

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

trunk/MagicSoft/Mars/manalysis/MCerPhotCalc.cc

@@ -83 +83 @@
 Bool_t MCerPhotCalc::PreProcess(MParList *pList)
 {
-    MRawRunHeader *RunHeader = (MRawRunHeader*)pList->FindObject("MRawRunHeader");
-    if (!RunHeader)
+    fRunHeader = (MRawRunHeader*)pList->FindObject("MRawRunHeader");
+    if (!fRunHeader)
     {
         *fLog << dbginf << "MRawRunHeader not found... aborting." << endl;
         return kFALSE;
     }
-    fNumLoGainFadcSlices= RunHeader->GetNumSamplesLoGain();
-    fNumHiGainFadcSlices= RunHeader->GetNumSamplesHiGain();
 
     fRawEvt = (MRawEvtData*)pList->FindObject("MRawEvtData");
@@ -133 +131 @@
         Float_t nphot = (Float_t)pixel.GetSumHiGainSamples();
 
-        if(nphot>0) nphot -= fNumHiGainFadcSlices*ped.GetMean();
+        // We check that the pixel is not empty, if it is empty
+        // we won't substract the pedestal. Empty means that it has
+        // 0 signal in all the slices.
 
-        fCerPhotEvt->AddPixel(pixid, nphot, sqrt(fNumHiGainFadcSlices)*ped.GetMeanRms());
+        if(nphot!=0) nphot -= fRunHeader->GetNumSamplesHiGain()*ped.GetMean();
+
+        fCerPhotEvt->AddPixel(pixid, nphot, sqrt(fRunHeader->GetNumSamplesHiGain())*ped.GetSigma());
 
         // FIXME! Handling of Lo Gains is missing!
@@ -144 +146 @@
     return kTRUE;
 }
+