Changeset 2669 for trunk/MagicSoft/Mars

Timestamp:

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

Author:

moralejo

Message:

*** empty log message ***

Location:

trunk/MagicSoft/Mars

Files:

• Changelog (modified) (1 diff)
• manalysis/MMcCalibrationCalc.cc (modified) (7 diffs)
• manalysis/MMcCalibrationCalc.h (modified) (2 diffs)

trunk/MagicSoft/Mars/Changelog

@@ -4 +4 @@
 
                                                  -*-*- END OF LINE -*-*-
+ 2003/12/12: Abelardo Moralejo
+
+   * manalysis/MMcCalibrationCalc.[h,cc]
+     - added fADC2PheInner, fADC2PheOuter and setter functions. Put 
+       these values into MCalibrationCam. Look for MExtractedSignalCam
+       to check which gain branch has been used (relevant for the
+       conversion factor) for each pixel.
+
 
  2003/12/12: Markus Gaug

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.);
+
     }
 

trunk/MagicSoft/Mars/manalysis/MMcCalibrationCalc.h

@@ -8 +8 @@
 class MCalibrationCam;
 class MMcFadcHeader;
+class MExtractedSignalCam;
 
 class MMcCalibrationCalc : public MTask
@@ -17 +18 @@
     Bool_t ReInit(MParList *pList);
 
-    MCalibrationCam *fCalCam;
-    MMcFadcHeader   *fHeaderFadc;
+    MCalibrationCam     *fCalCam;
+    MMcFadcHeader       *fHeaderFadc;
+    MExtractedSignalCam *fSignalCam;
+
+    Float_t fADC2PheInner;
+    Float_t fADC2PheOuter;
 
 public:
     MMcCalibrationCalc(const char *name=NULL, const char *title=NULL);
+
+    void SetADC2PheInner(Float_t x) {fADC2PheInner = x; }
+    void SetADC2PheOuter(Float_t x) {fADC2PheOuter = x; }
 
     ClassDef(MMcCalibrationCalc, 0)   // Task which obtains, for MC files, the pedestal mean and rms, and the calibration factor from ADC counts to photons.