Changeset 2728 for trunk/MagicSoft/Mars/manalysis/MCalibrate.cc

Timestamp:

12/20/03 13:46:17 (21 years ago)

Author:

tbretz

Message:

*** empty log message ***

File:

• trunk/MagicSoft/Mars/manalysis/MCalibrate.cc (modified) (5 diffs)

trunk/MagicSoft/Mars/manalysis/MCalibrate.cc

@@ -24 +24 @@
 
 //////////////////////////////////////////////////////////////////////////////
-//                                                                          //
-//   MCalibrate                                                             //
-//                                                                          //
-//   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.                      //
-//                                                                          //
-//   Input Containers:                                                      //
-//                                                                          //     
-//   MGeomCam                                                               //
-//   MExtractedSingal                                                       //
-//   MCalibrationCam                                                        //
-//                                                                          //
-//   Output Containers:                                                     //
-//                                                                          //
-//   MCerPhotEvt                                                            //
-//                                                                          //
+//
+//   MCalibrate
+//
+//   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.
+//
+//   Input Containers:
+//    MExtractedSingal
+//    MCalibrationCam
+//
+//   Output Containers:
+//    MCerPhotEvt
+//
 //////////////////////////////////////////////////////////////////////////////
-
 #include "MCalibrate.h"
 #include "MCalibrationConfig.h"
@@ -74 +70 @@
 MCalibrate::MCalibrate(const char *name, const char *title)
 {
-
     fName  = name  ? name  : "MCalibrate";
     fTitle = title ? title : "Task to calculate the number of photons in one event";
-
 }
 
@@ -94 +88 @@
 Int_t MCalibrate::PreProcess(MParList *pList)
 {
-
-    fGeomCam = (MGeomCam*)pList->FindObject("MGeomCam");
-    if (!fGeomCam)
-      {
-        *fLog << err << dbginf << "MGeomCam not found ... aborting." << endl;        
-        return kFALSE;
-      }
-
     fSignals = (MExtractedSignalCam*)pList->FindObject("MExtractedSignalCam");
     if (!fSignals)
-      {
-        *fLog << err << dbginf << "MExtractedSignalCam not found ... aborting" << endl;
+    {
+        *fLog << err << "MExtractedSignalCam not found ... aborting" << endl;
         return kFALSE;
-      }
+    }
 
     fCalibrations = (MCalibrationCam*)pList->FindObject("MCalibrationCam");
     if (!fCalibrations)
-      {
-        *fLog << err << dbginf << "MCalibrationCam not found ... aborting." << endl;        
+    {
+        *fLog << err << "MCalibrationCam not found ... aborting." << endl;
         return kFALSE;
-      }
+    }
 
     fCerPhotEvt = (MCerPhotEvt*)pList->FindCreateObj("MCerPhotEvt");
     if (!fCerPhotEvt)
-      {
-        *fLog << err << dbginf << "Cannot create MCerPhotEvt ... aborting" << endl;
         return kFALSE;
-      }
     
     return kTRUE;
 }
-
 
 // --------------------------------------------------------------------------
@@ -132 +114 @@
 Int_t MCalibrate::Process()
 {
-    UInt_t imaxnumpix = fGeomCam->GetNumPixels();
+    if (fCalibrations->GetNumPixels() != (UInt_t)fSignals->GetSize())
+    {
+        // FIXME: MExtractedSignal must be of variable size -
+        //        like MCerPhotEvt - because we must be able
+        //        to reduce size by zero supression
+        //        For the moment this check could be done in ReInit...
+        *fLog << err << "MExtractedSignal and MCalibrationCam have different sizes... abort." << endl;
+        return kFALSE;
+    }
+
+    const UInt_t imaxnumpix = fSignals->GetSize();
     
-    for (UInt_t pixid = 0; pixid < imaxnumpix; pixid++)
+    for (UInt_t pixidx=0; pixidx<imaxnumpix; pixidx++)
     {
-        MCalibrationPix &pix = (*fCalibrations)[pixid];
+        const MCalibrationPix &pix = (*fCalibrations)[pixidx];
 
+        if (!pix.IsBlindPixelMethodValid())
+            continue;
 
-        if (pix.IsBlindPixelMethodValid())
-        {
-            MExtractedSignalPix &sig =  (*fSignals)[pixid];
+        MExtractedSignalPix &sig =  (*fSignals)[pixidx];
 
-            Float_t signal;
-            Float_t signalErr = 0.;
+        Float_t signal;
+        Float_t signalErr = 0.;
 
-            if (sig.IsLoGainUsed())
-              {
-                signal    = sig.GetExtractedSignalLoGain()*pix.GetConversionHiLo();
-                signalErr = signal*pix.GetConversionHiLoError();
-              }
-            else
-              {
-                signal = sig.GetExtractedSignalHiGain();
-              }
-            
-            //      Float_t calibrationConvertionFactor = pix.GetMeanConversionFFactorMethod();
-            Float_t calibrationConversionFactor      = pix.GetMeanConversionBlindPixelMethod();
-            Float_t calibrationConversionFactorError = pix.GetErrorConversionBlindPixelMethod();
+        if (sig.IsLoGainUsed())
+        {
+            signal    = sig.GetExtractedSignalLoGain()*pix.GetConversionHiLo();
+            signalErr = signal*pix.GetConversionHiLoError();
+        }
+        else
+        {
+            signal = sig.GetExtractedSignalHiGain();
+        }
 
-            Float_t nphot    = signal*calibrationConversionFactor;
-            Float_t nphotErr = signal*calibrationConversionFactorError
-                              *signal*calibrationConversionFactorError;
-            nphotErr        += signalErr*calibrationConversionFactor
-                              *signalErr*calibrationConversionFactor;
+        //      Float_t calibrationConvertionFactor = pix.GetMeanConversionFFactorMethod();
+        const Float_t calibrationConversionFactor      = pix.GetMeanConversionBlindPixelMethod();
+        const Float_t calibrationConversionFactorError = pix.GetErrorConversionBlindPixelMethod();
 
-            nphotErr  = TMath::Sqrt(nphotErr);
+        const Float_t nphot = signal*calibrationConversionFactor;
+        Float_t nphotErr    = signal*calibrationConversionFactorError
+                             *signal*calibrationConversionFactorError;
+        nphotErr += signalErr*calibrationConversionFactor
+                   *signalErr*calibrationConversionFactor;
 
-            fCerPhotEvt->AddPixel(pixid, nphot, nphotErr);
-        }
-        else
-          {
-            fCerPhotEvt->AddPixel(pixid, 0., 0.);
-            (*fCerPhotEvt)[pixid].SetPixelUnused();
-          }
+        nphotErr  = TMath::Sqrt(nphotErr);
+
+        fCerPhotEvt->AddPixel(pixidx, nphot, nphotErr);
     }
 
@@ -182 +168 @@
     return kTRUE;
 }
-