Changeset 2721 for trunk/MagicSoft/Mars/manalysis

Timestamp:

12/18/03 20:14:53 (21 years ago)

Author:

moralejo

Message:

*** empty log message ***

File:

• trunk/MagicSoft/Mars/manalysis/MMcCalibrationCalc.cc (modified) (8 diffs)

trunk/MagicSoft/Mars/manalysis/MMcCalibrationCalc.cc

@@ -54 +54 @@
 #include "MExtractedSignalCam.h"
 #include "MExtractedSignalPix.h"
+#include "MGeomCam.h"
 #include "MPedPhotCam.h"
 #include "MPedPhotPix.h"
@@ -69 +70 @@
     fTitle = title ? title : "Write MC pedestals and conversion factors into MCalibration Container";
 
-    fADC2PheInner = 1.;
-    fADC2PheOuter = 1.;
+    fADC2PhInner = 1.;
+    fADC2PhOuter = 1.;
 
 }
@@ -110 +111 @@
     }
 
+  if ( ! pList->FindCreateObj(AddSerialNumber("MPedPhotCam")))
+    {
+      *fLog << err << dbginf << "Cannot create MPedPhotCam... aborting." << endl;
+      return kFALSE;
+    }
+
+  fGeom = (MGeomCam*) pList->FindObject(AddSerialNumber("MGeomCam"));
+  if ( ! fGeom )
+    {
+      *fLog << err << dbginf << "Cannot find MGeomCam... aborting." << endl;
+      return kFALSE;
+    }
+
   fSignalCam = (MExtractedSignalCam*) pList->FindObject(AddSerialNumber("MExtractedSignalCam"));
   if ( ! fSignalCam)
@@ -124 +138 @@
 //
 // Check for the runtype.
-// Search for MCalibrationCam and MMcFadcHeader.
+// Search for MCalibrationCam, MPedPhotCam and MMcFadcHeader.
 //
 Bool_t MMcCalibrationCalc::ReInit(MParList *pList)
@@ -154 +168 @@
     }
 
+    fPedPhotCam = (MPedPhotCam*)pList->FindObject(AddSerialNumber("MPedPhotCam"));
+
+    if (!fPedPhotCam)
+    {
+        *fLog << err << dbginf << "Could not create MPedPhotCam... aborting. " << endl;
+        return kFALSE;
+    }
+
     return kTRUE;
 }
@@ -160 +182 @@
 // --------------------------------------------------------------------------
 //
-// Fill the MCalibrationCam object
+// Fill the MCalibrationCam and MCerPhotPed objects
 //
 Int_t MMcCalibrationCalc::Process()
@@ -171 +193 @@
         MExtractedSignalPix &sigpix = (*fSignalCam)[i];
 
+        //
+        // ped mean and rms per pixel, in ADC counts, according to signal 
+        // calculation (hi or low gain and number of integrated slices):
+        //
         const Float_t pedestmean = sigpix.IsLoGainUsed()? 
           fSignalCam->GetNumUsedLoGainFADCSlices()*fHeaderFadc->GetPedestal(i) : 
@@ -179 +205 @@
           sqrt((Double_t)(fSignalCam->GetNumUsedHiGainFADCSlices())) * fHeaderFadc->GetPedestalRmsHigh(i);
 
-        // ped mean and rms per pixel, according to signal calculation:
-
         calpix.SetPedestal(pedestmean, pedestrms);
-
-        if (i < 397) // FIXME!
-          calpix.SetConversionBlindPixelMethod(fADC2PheInner, 0., 0.);
+        calpix.SetBlindPixelMethodValid();
+
+        if ( fGeom->GetPixRatio(i) > fGeom->GetPixRatio(0) )
+          calpix.SetConversionBlindPixelMethod(fADC2PhOuter, 0., 0.);
         else
-          calpix.SetConversionBlindPixelMethod(fADC2PheOuter, 0., 0.);
+          calpix.SetConversionBlindPixelMethod(fADC2PhInner, 0., 0.);
 
     }