Changeset 2945 for trunk/MagicSoft/Mars/mcalib/MCalibrate.cc

Timestamp:

01/28/04 15:05:30 (21 years ago)

Author:

rico

Message:

*** empty log message ***

File:

• trunk/MagicSoft/Mars/mcalib/MCalibrate.cc (modified) (8 diffs)

trunk/MagicSoft/Mars/mcalib/MCalibrate.cc

@@ -17 +17 @@
 !
 !   Author(s): Javier Lopez 12/2003 <mailto:jlopez@ifae.es>
-!
-!   Copyright: MAGIC Software Development, 2000-2001
+!   Modified by: Javier Rico  01/2004 <mailto:jrico@ifae.es>
+!
+!   Copyright: MAGIC Software Development, 2000-2004
 !
 !
@@ -29 +30 @@
 //   This task takes the integrated charge from MExtractedSignal and apply
 //   the calibration constants from MCalibraitionCam to the charge. Then
-//   stores number of photons obtained in MCerPhotEvt.
+//   stores number of photons obtained in MCerPhotEvt. Selection of different
+//   calibration methods is allowed through SetCalibrationMode member function
 //
 //   Input Containers:
@@ -65 +67 @@
 // Default constructor. 
 //
-MCalibrate::MCalibrate(const char *name, const char *title)
+MCalibrate::MCalibrate(CalibrationMode_t calmode,const char *name, const char *title) : fCalibrationMode(calmode)
 {
     fName  = name  ? name  : "MCalibrate";
@@ -85 +87 @@
 Int_t MCalibrate::PreProcess(MParList *pList)
 {
+
     fSignals = (MExtractedSignalCam*)pList->FindObject(AddSerialNumber("MExtractedSignalCam"));
+
     if (!fSignals)
     {
@@ -92 +96 @@
     }
 
-    fCalibrations = (MCalibrationCam*)pList->FindObject(AddSerialNumber("MCalibrationCam"));
-    if (!fCalibrations)
-    {
-      *fLog << err << AddSerialNumber("MCalibrationCam") << " not found ... aborting." << endl;
-        return kFALSE;
-    }
+   if(fCalibrationMode>kNone)
+      {
+
+        fCalibrations = (MCalibrationCam*)pList->FindObject(AddSerialNumber("MCalibrationCam"));
+        if (!fCalibrations)
+          {
+            *fLog << err << AddSerialNumber("MCalibrationCam") << " not found ... aborting." << endl;
+            return kFALSE;
+          }
+
+      }
 
     fCerPhotEvt = (MCerPhotEvt*)pList->FindCreateObj(AddSerialNumber("MCerPhotEvt"));
@@ -108 +117 @@
 // --------------------------------------------------------------------------
 //
+// Check for validity of the selected calibration method, switch to a 
+// different one in case of need
+//
+Int_t MCalibrate::ReInit()
+{
+  if(fCalibrationMode == kBlindPixel && !fCalibrations->IsBlindPixelMethodValid())
+    {
+      *fLog << warn << GetDescriptor() << "Warning: Blind pixel calibration method not valid, switching to F-factor method" << endl;
+      fCalibrationMode = kFfactor;
+    }
+
+  if(fCalibrationMode == kPinDiode && !fCalibrations->IsPINDiodeMethodValid())
+    { 
+      *fLog << warn << GetDescriptor() << "Warning: PIN diode calibration method not valid, switching to F-factor method" << endl;
+      fCalibrationMode = kFfactor;
+    }
+  return kTRUE;
+}
+// --------------------------------------------------------------------------
+//
+// Apply the calibration factors to the extracted signal according to the 
+// selected calibration method
 //
 Int_t MCalibrate::Process()
@@ -123 +154 @@
   */
 
-  UInt_t numpix = fSignals->GetSize();
-
-  for (UInt_t pixid=0; pixid< numpix; pixid++)
-    {
-      const MCalibrationPix &pix = (*fCalibrations)[pixid];
-      
-      if (!pix.IsFitValid())
-        continue;
-
-      MExtractedSignalPix &sig =  (*fSignals)[pixid];
+  UInt_t npix = fSignals->GetSize();
+
+  Float_t hiloconv   = 1.;
+  Float_t hiloconverr= 0.;
+  Float_t calibrationConversionFactor = 1.;
+  Float_t calibrationConversionFactorError = 0.;
+  
+  for (UInt_t pixidx=0; pixidx<npix; pixidx++)
+    {
+      if(fCalibrationMode!=kNone)
+        {
+          const MCalibrationPix &pix = (*fCalibrations)[pixidx];       
+          
+          if (!pix.IsFitValid())
+            continue;
+          
+          hiloconv   = pix.GetConversionHiLo();
+          hiloconverr= pix.GetConversionHiLoError();
+          
+          switch(fCalibrationMode)
+            {
+            case kBlindPixel:
+              calibrationConversionFactor      = pix.GetMeanConversionBlindPixelMethod();
+              calibrationConversionFactorError = pix.GetErrorConversionBlindPixelMethod();
+              break;
+            case kFfactor:
+              calibrationConversionFactor      = pix.GetMeanConversionFFactorMethod();
+              calibrationConversionFactorError = pix.GetErrorConversionFFactorMethod();
+              break;
+            default:
+              *fLog << warn << "MCalibrate::Process Warning: Calibration mode value ("<<fCalibrationMode<<") not known" << endl;
+              break;
+            }
+        }
+      
+      MExtractedSignalPix &sig =  (*fSignals)[pixidx];
       
       Float_t signal;
       Float_t signalErr = 0.;
-      
+      Float_t nphot,nphotErr;
+            
       if (sig.IsLoGainUsed())
         {
-          signal    = sig.GetExtractedSignalLoGain()*pix.GetConversionHiLo();
-          signalErr = signal*pix.GetConversionHiLoError();
+          signal    = sig.GetExtractedSignalLoGain()*hiloconv;
+          signalErr = signal*hiloconverr;
         }
       else
@@ -152 +210 @@
             signal = sig.GetExtractedSignalHiGain();
         }
-
-      //      Float_t calibrationConversionFactor = pix.GetMeanConversionFFactorMethod();
-      const Float_t calibrationConversionFactor      = pix.GetMeanConversionBlindPixelMethod();
-      const Float_t calibrationConversionFactorError = pix.GetErrorConversionBlindPixelMethod();
-      
-      const Float_t nphot = signal*calibrationConversionFactor;
-
-      Float_t nphotErr    = signal*calibrationConversionFactorError
-                           *signal*calibrationConversionFactorError;
-
-      nphotErr += signalErr*calibrationConversionFactor
-                 *signalErr*calibrationConversionFactor;
+      
+      nphot    = signal*calibrationConversionFactor;
+      nphotErr = signal*calibrationConversionFactorError
+        *signal*calibrationConversionFactorError
+        +signalErr*calibrationConversionFactor
+        *signalErr*calibrationConversionFactor;
       
       nphotErr  = TMath::Sqrt(nphotErr);
       
-      MCerPhotPix *cpix = fCerPhotEvt->AddPixel(pixid, nphot, nphotErr);
+      MCerPhotPix *cpix = fCerPhotEvt->AddPixel(pixidx, nphot, nphotErr);
 
       if (sig.GetNumLoGainSaturated() > 0)
-          cpix->SetPixelSaturated();
+        cpix->SetPixelSaturated();
     }