Changeset 1291 for trunk/MagicSoft/Mars/manalysis

Timestamp:

04/24/02 09:23:07 (23 years ago)

Author:

tbretz

Message:

*** empty log message ***

Location:

trunk/MagicSoft/Mars/manalysis

Files:

• MCerPhotCalc.cc (modified) (3 diffs)
• MCerPhotCalc.h (modified) (2 diffs)

trunk/MagicSoft/Mars/manalysis/MCerPhotCalc.cc

@@ -44 +44 @@
 #include "MLog.h"
 #include "MLogManip.h"
+
+#include "MMcRunHeader.hxx"
 
 #include "MRawRunHeader.h"
@@ -113 +115 @@
 // --------------------------------------------------------------------------
 //
+// Check for the run type and camera version.
+// If the file is a MC file and the used camera version is <= 40
+// we enable a fix for truncated pedestal means in this version.
+//
+Bool_t MCerPhotCalc::ReInit(MParList *pList)
+{
+    const MRawRunHeader *runheader = (MRawRunHeader*)pList->FindObject("MRawRunHeader");
+    if (!runheader)
+    {
+        *fLog << warn << dbginf << "Warning - cannot check file type, MRawRunHeader not found." << endl;
+        return kTRUE;
+    }
+
+    if (runheader->GetRunType() != kRTMonteCarlo)
+        return kTRUE;
+
+    MMcRunHeader *mcrunheader = (MMcRunHeader*)pList->FindObject("MMcRunHeader");
+    if (!mcrunheader)
+    {
+        *fLog << warn << dbginf << "Warning - cannot check for camera version, MC run header not found." << endl;
+        return kTRUE;
+    }
+
+    if (mcrunheader->GetCamVersion() <= 40)
+        fEnableFix = kTRUE;
+
+    return kTRUE;
+}
+
+// --------------------------------------------------------------------------
+//
 // Calculate the integral of the FADC time slices and store them as a new
 // pixel in the MCerPhotEvt container.
@@ -144 +177 @@
         // 0 signal in all the slices.
         //
+        const Double_t mean = fEnableFix ? ped.GetMean()-0.5 : ped.GetMean();
         if (nphot!=0)
-            nphot -= fRunHeader->GetNumSamplesHiGain()*ped.GetMean();
+            nphot -= fRunHeader->GetNumSamplesHiGain()*mean;
 
         fCerPhotEvt->AddPixel(pixid, nphot, sqrt(fRunHeader->GetNumSamplesHiGain())*ped.GetSigma());

trunk/MagicSoft/Mars/manalysis/MCerPhotCalc.h

@@ -26 +26 @@
     MCerPhotEvt   *fCerPhotEvt; // Cerenkov Photon Event used for calculation
     MRawRunHeader *fRunHeader;  //  RunHeader information
- 
+
+    Bool_t         fEnableFix;  // fix for a bug in files from older camera versions (<=40)
+
 public:
     MCerPhotCalc(const char *name=NULL, const char *title=NULL);
@@ -33 +35 @@
     Bool_t Process();
 
+    Bool_t ReInit(MParList *pList);
+
     ClassDef(MCerPhotCalc, 0)   // Task to calculate cerenkov photons from raw data
 };