Changeset 2669 for trunk/MagicSoft/Mars/manalysis/MMcCalibrationCalc.cc

Timestamp:

12/12/03 15:00:48 (21 years ago)

Author:

moralejo

Message:

*** empty log message ***

File:

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

trunk/MagicSoft/Mars/manalysis/MMcCalibrationCalc.cc

@@ -28 +28 @@
 //
 //  This task looks for the ìnformation about FADC pedestals in
-//  MMcFadcHeader and translates it to the pedestal mean and rms,
+//  MMcFadcHeader and translates it to the pedestal mean and rms (in adc counts),
 //  and in the conversion factor between photons and ADC counts in
 //  MCalibrationCam.
+//  Then we create and fill also the MPedPhotCam object containing the pedestal
+//  rms in units of photons.
 //
 //  Input Containers:
@@ -38 +40 @@
 //  Output Containers:
 //   MCalibrationCam
+//   MPedPhotCam
 //
 /////////////////////////////////////////////////////////////////////////////
@@ -49 +52 @@
 #include "MCalibrationPix.h"
 #include "MCalibrationCam.h"
+#include "MExtractedSignalCam.h"
+#include "MExtractedSignalPix.h"
+#include "MPedPhotCam.h"
+#include "MPedPhotPix.h"
 
 #include "MRawRunHeader.h"
@@ -61 +68 @@
     fName  = name  ? name  : "MMcCalibrationCalc";
     fTitle = title ? title : "Write MC pedestals and conversion factors into MCalibration Container";
+
+    fADC2PheInner = 1.;
+    fADC2PheOuter = 1.;
 
 }
@@ -94 +104 @@
     return kTRUE;
 
-  return pList->FindCreateObj(AddSerialNumber("MCalibrationCam")) ? kTRUE : kFALSE;
+  if ( ! pList->FindCreateObj(AddSerialNumber("MCalibrationCam")))
+    {
+      *fLog << err << dbginf << "Cannot create MCalibrationCam... aborting." << endl;
+      return kFALSE;
+    }
+
+  fSignalCam = (MExtractedSignalCam*) pList->FindObject(AddSerialNumber("MExtractedSignalCam"));
+  if ( ! fSignalCam)
+    {
+      *fLog << err << dbginf << "Cannot find MExtractedSignalCam... aborting." << endl;
+      return kFALSE;
+    }
+
+  return kTRUE;
+
 }
 
@@ -122 +146 @@
     }
 
-    fCalCam = (MCalibrationCam*)pList->FindCreateObj(AddSerialNumber("MCalibrationCam"));
+    fCalCam = (MCalibrationCam*)pList->FindObject(AddSerialNumber("MCalibrationCam"));
 
     if (!fCalCam)
@@ -144 +168 @@
     for (int i=0; i<num; i++)
     {
-        MCalibrationPix &pix = (*fCalCam)[i];
+        MCalibrationPix &calpix = (*fCalCam)[i];
+        MExtractedSignalPix &sigpix = (*fSignalCam)[i];
 
-        const Float_t pedestmean = fHeaderFadc->GetPedestal(i);
-        const Float_t pedestrms  = fHeaderFadc->GetPedestalRmsHigh(i); 
-        // per slice
+        const Float_t pedestmean = sigpix.IsLoGainUsed()? 
+          fSignalCam->GetNumUsedLoGainFADCSlices()*fHeaderFadc->GetPedestal(i) : 
+          fSignalCam->GetNumUsedHiGainFADCSlices()*fHeaderFadc->GetPedestal(i);
 
-        pix.SetPedestal(pedestmean, pedestrms);
+        const Float_t pedestrms  = sigpix.IsLoGainUsed()? 
+          sqrt((Double_t)(fSignalCam->GetNumUsedLoGainFADCSlices())) * fHeaderFadc->GetPedestalRmsLow(i) : 
+          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.);
+        else
+          calpix.SetConversionBlindPixelMethod(fADC2PheOuter, 0., 0.);
+
     }