Changeset 8304 for trunk/MagicSoft/Mars/mhcalib/MHCalibrationChargeCam.cc

Timestamp:

02/04/07 15:33:08 (18 years ago)

Author:

tbretz

Message:

*** empty log message ***

File:

• trunk/MagicSoft/Mars/mhcalib/MHCalibrationChargeCam.cc (modified) (14 diffs)

trunk/MagicSoft/Mars/mhcalib/MHCalibrationChargeCam.cc

@@ -1 +1 @@
 /* ======================================================================== *\
-! $Name: not supported by cvs2svn $:$Id: MHCalibrationChargeCam.cc,v 1.49 2007-01-30 14:16:45 tbretz Exp $
+! $Name: not supported by cvs2svn $:$Id: MHCalibrationChargeCam.cc,v 1.50 2007-02-04 15:29:09 tbretz Exp $
 ! --------------------------------------------------------------------------
 !
@@ -114 +114 @@
 //  + Float_t fNumLoGainBlackoutLimit; // Rel. amount blackout logain events until pixel is declared unsuitable
 //
+// Class Version 3:
+//  - Byte_t  fFirstHiGain;                              // First used slice High Gain
+//  - Byte_t  fLastHiGain;                               // Last used slice High Gain
+//  - Byte_t  fFirstLoGain;                              // First used slice Low Gain
+//  - Byte_t  fLastLoGain;                               // Last used slice Low Gain
+//
 /////////////////////////////////////////////////////////////////////////////
 #include "MHCalibrationChargeCam.h"
@@ -219 +225 @@
 //
 MHCalibrationChargeCam::MHCalibrationChargeCam(const char *name, const char *title)
-    : fRawEvt(NULL)
+    : fRawEvt(NULL), fSignal(NULL)
 {
 
@@ -367 +373 @@
 {
 
-  MExtractedSignalCam *signal = 
-    (MExtractedSignalCam*)pList->FindObject(AddSerialNumber("MExtractedSignalCam"));
-  if (!signal)
+  fSignal = (MExtractedSignalCam*)pList->FindObject(AddSerialNumber("MExtractedSignalCam"));
+  if (!fSignal)
   {
       *fLog << err << "MExtractedSignalCam not found... abort." << endl;
@@ -377 +382 @@
   if (!InitCams(pList,"Charge"))
     return kFALSE;
-
-  fFirstHiGain = signal->GetFirstUsedSliceHiGain();
-  fLastHiGain  = signal->GetLastUsedSliceHiGain();
-  fFirstLoGain = signal->GetFirstUsedSliceLoGain();
-  fLastLoGain  = signal->GetLastUsedSliceLoGain();
 
   const Int_t npixels  = fGeom->GetNumPixels();
@@ -900 +900 @@
   MBadPixelsCam   *badcam    = fIntensBad ? fIntensBad->GetCam() : fBadPixels;
 
+  const Int_t hifirst = fSignal->GetFirstUsedSliceHiGain();
+  const Int_t hilast  = fSignal->GetLastUsedSliceHiGain();
+  const Int_t lofirst = fSignal->GetFirstUsedSliceLoGain();
+  const Int_t lolast  = fSignal->GetLastUsedSliceLoGain();
+
 
   for (Int_t i=0; i<fHiGainArray->GetSize(); i++)
@@ -944 +949 @@
           bad.SetUncalibrated( MBadPixelsPix::kHiGainOverFlow ); 
         }
-      
-      FinalizeAbsTimes(histhi, pix, bad, fFirstHiGain, fLastHiGain);
+
+      FinalizeAbsTimes(histhi, pix, bad, hifirst, hilast);
     }
 
@@ -991 +996 @@
           }
         
-        FinalizeAbsTimes(histlo, pix, bad, fFirstLoGain, fLastLoGain);
+        FinalizeAbsTimes(histlo, pix, bad, lofirst, lolast);
       }
 
@@ -1009 +1014 @@
 
       MBadPixelsPix &bad = chargecam->GetAverageBadArea(j);
-      FinalizeAbsTimes(histhi, pix, bad, fFirstHiGain, fLastHiGain);
+      FinalizeAbsTimes(histhi, pix, bad, hifirst, hilast);
    }
 
@@ -1031 +1036 @@
           {
             MBadPixelsPix &bad = chargecam->GetAverageBadArea(j);
-            FinalizeAbsTimes(histlo, pix, bad, fFirstLoGain, fLastLoGain);
+            FinalizeAbsTimes(histlo, pix, bad, lofirst, lolast);
           }
         
@@ -1051 +1056 @@
 
       MBadPixelsPix &bad = chargecam->GetAverageBadSector(j);
-      FinalizeAbsTimes(histhi, pix, bad, fFirstHiGain, fLastHiGain);
+      FinalizeAbsTimes(histhi, pix, bad, hifirst, hilast);
     }
   
@@ -1073 +1078 @@
 
         if (pix.IsHiGainSaturation())
-          FinalizeAbsTimes(histlo, pix, bad, fFirstLoGain, fLastLoGain);
+          FinalizeAbsTimes(histlo, pix, bad, lofirst, lolast);
       }
   
@@ -1131 +1136 @@
 //
 void MHCalibrationChargeCam::FinalizeAbsTimes(MHCalibrationChargePix &hist, MCalibrationChargePix &pix, MBadPixelsPix &bad, 
-                                              Byte_t first, Byte_t last)
+                                              Int_t first, Int_t last)
 {
     const Float_t mean = hist.GetAbsTimeMean();
@@ -1139 +1144 @@
     pix.SetAbsTimeRms(rms);
 
-    const Float_t lowerlimit = (Float_t)first;// + fTimeLowerLimit;
-    const Float_t upperlimit = (Float_t)last ;// - fTimeUpperLimit;
+    const Float_t lowerlimit = (Float_t)first+1;// + fTimeLowerLimit;
+    const Float_t upperlimit = (Float_t)last -1;// - fTimeUpperLimit;
 
     // FIXME: instead of checking whether the maximum is in the first or