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

Timestamp:

01/20/04 11:34:33 (21 years ago)

Author:

gaug

Message:

*** empty log message ***

File:

• trunk/MagicSoft/Mars/mcalib/MCalibrate.cc (modified) (1 diff)

trunk/MagicSoft/Mars/mcalib/MCalibrate.cc

@@ -125 +125 @@
   */
 
-    const UInt_t imaxnumpix = fSignals->GetSize();
+  UInt_t numpix = fSignals->GetSize();
 
-    for (UInt_t pixidx=0; pixidx<imaxnumpix; pixidx++)
+  for (UInt_t pixid=0; pixid< numpix; pixid++)
     {
-        const MCalibrationPix &pix = (*fCalibrations)[pixidx];
+      const MCalibrationPix &pix = (*fCalibrations)[pixid];
+      
+      if (!pix.IsFitValid())
+        continue;
 
-        if (!pix.IsFitValid())
-            continue;
+      MExtractedSignalPix &sig =  (*fSignals)[pixid];
+      
+      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 calibrationConversionFactor = pix.GetMeanConversionFFactorMethod();
+      const Float_t calibrationConversionFactor      = pix.GetMeanConversionBlindPixelMethod();
+      const Float_t calibrationConversionFactorError = pix.GetErrorConversionBlindPixelMethod();
+      
+      const Float_t nphot = signal*calibrationConversionFactor;
 
-        MExtractedSignalPix &sig =  (*fSignals)[pixidx];
+      Float_t nphotErr    = signal*calibrationConversionFactorError
+                           *signal*calibrationConversionFactorError;
 
-        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 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;
-
-        nphotErr  = TMath::Sqrt(nphotErr);
-
-
-        fCerPhotEvt->AddPixel(pixidx, nphot, nphotErr);
-
-        if (sig.GetNumLoGainSaturated() > 0)
-          (*fCerPhotEvt)[pixidx].SetPixelSaturated();
+      nphotErr += signalErr*calibrationConversionFactor
+                 *signalErr*calibrationConversionFactor;
+      
+      nphotErr  = TMath::Sqrt(nphotErr);
+      
+      fCerPhotEvt->AddPixel(pixid, nphot, nphotErr);
+      
+      if (sig.GetNumLoGainSaturated() > 0)
+        (*fCerPhotEvt)[pixid].SetPixelSaturated();
     }
-
-    fCerPhotEvt->FixSize();
-    fCerPhotEvt->SetReadyToSave();
-
-    return kTRUE;
+  
+  fCerPhotEvt->FixSize();
+  fCerPhotEvt->SetReadyToSave();
+  
+  return kTRUE;
 }