Changeset 4793 for trunk/MagicSoft/Mars/mcalib/MHCalibrationChargeCam.cc

Timestamp:

08/29/04 17:04:39 (20 years ago)

Author:

gaug

Message:

*** empty log message ***

File:

• trunk/MagicSoft/Mars/mcalib/MHCalibrationChargeCam.cc (modified) (13 diffs)

trunk/MagicSoft/Mars/mcalib/MHCalibrationChargeCam.cc

@@ -154 +154 @@
 using namespace std;
 
+const Int_t   MHCalibrationChargeCam::fgAverageNbinsHiGain       = 3300;
+const Int_t   MHCalibrationChargeCam::fgAverageNbinsLoGain       = 1950;
 const Float_t MHCalibrationChargeCam::fgNumHiGainSaturationLimit = 0.01;
 const Float_t MHCalibrationChargeCam::fgNumLoGainSaturationLimit = 0.005;
@@ -175 +177 @@
 // - fTimeLowerLimit           to fgTimeLowerLimit 
 // - fTimeUpperLimit           to fgTimeUpperLimit 
+// - fAverageNbins             to fgAverageNbinsHiGain
 //
 MHCalibrationChargeCam::MHCalibrationChargeCam(const char *name, const char *title)
     : fRawEvt(NULL)
 {
+
   fName  = name  ? name  : "MHCalibrationChargeCam";
   fTitle = title ? title : "Class to fill the calibration histograms ";
@@ -187 +191 @@
   SetTimeUpperLimit();
 
-
+  SetAverageNbins(fgAverageNbinsHiGain);
 }
 
@@ -220 +224 @@
 // Initializes, if empty to MGeomCam::GetNumPixels():
 // - MHCalibrationCam::fHiGainArray, MHCalibrationCam::fLoGainArray
+// - MHCalibrationCam::fHiGainOverflow, MHCalibrationCam::fLoGainOverflow
 //
 // Initializes, if empty to MGeomCam::GetNumAreas() for:
@@ -243 +248 @@
 {
 
-  MExtractedSignalCam *signal = (MExtractedSignalCam*)pList->FindObject(AddSerialNumber("MExtractedSignalCam"));
+  MExtractedSignalCam *signal = 
+    (MExtractedSignalCam*)pList->FindObject(AddSerialNumber("MExtractedSignalCam"));
   if (!signal)
   {
@@ -301 +307 @@
         InitHists((*this)[i],(*fBadPixels)[i],i);
       }
+      fHiGainOverFlow.Set(npixels);
   }
 
@@ -310 +317 @@
       {
           (*fLoGainArray)[i] = new MHCalibrationChargeLoGainPix;
-          InitHists((*this)(i),(*fBadPixels)[i],i);
+          MHGausEvents &lopix = (*this)(i);
+          //
+          // Adapt the range for the case, the intense blue is used:
+          // FIXME: this is a nasty workaround, but for the moment necessary 
+          // in order to avoid default memory space.
+          //
+          if (fGeom->InheritsFrom("MGeomCamMagic"))
+            {
+              if ( fColor == MCalibrationCam::kBLUE)
+                {
+                  lopix.SetLast(999.5);
+                  lopix.SetNbins(1150);
+                }
+            }
+          InitHists(lopix,(*fBadPixels)[i],i);
       }
-      
+      fLoGainOverFlow.Set(npixels);  
   }
 
@@ -357 +378 @@
 
         MHCalibrationChargePix &hist = (MHCalibrationChargePix&)GetAverageLoGainArea(j);
-
-        hist.SetNbins(fAverageNbins);
+        hist.SetNbins(fgAverageNbinsLoGain);
+
+        //
+        // Adapt the range for the case, the intense blue is used:
+        // FIXME: this is a nasty workaround, but for the moment necessary 
+        // in order to avoid default memory space.
+        //
+        if (fGeom->InheritsFrom("MGeomCamMagic"))
+            {
+              if ( fColor == MCalibrationCam::kBLUE)
+                {
+                  hist.SetFirst(-10.5);
+                  hist.SetLast(999.5);
+                  hist.SetNbins(3030);
+                }
+            }
+
         hist.GetHAbsTime()->SetTitle("Absolute Arrival Time average LoGain Area Idx ");
 
@@ -410 +446 @@
 
           hist.GetHGausHist()->SetTitle("Summed FADC slices average LoGain Sector ");
-          hist.SetNbins(fAverageNbins);
+          hist.SetNbins(fgAverageNbinsLoGain);
+          
+          //
+          // Adapt the range for the case, the intense blue is used:
+          // FIXME: this is a nasty workaround, but for the moment necessary 
+          // in order to avoid default memory space.
+          //
+          if (fGeom->InheritsFrom("MGeomCamMagic"))
+            {
+              if ( fColor == MCalibrationCam::kBLUE)
+                {
+                  hist.SetFirst(-10.5);
+                  hist.SetLast(999.5);
+                  hist.SetNbins(3030);
+                }
+            }
+          
           hist.GetHAbsTime()->SetTitle("Absolute Arrival Time average LoGain Sector ");
 
           InitHists(hist,fCam->GetAverageBadSector(j),j);
-          
       }
   }
@@ -494 +545 @@
 
       if (!histhi.FillHistAndArray(sumhi))
-        fHiGainOverFlow++;
+        fHiGainOverFlow[i]++;
+
       if (!histlo.FillHistAndArray(sumlo))
-        fLoGainOverFlow++;
-      
+        fLoGainOverFlow[i]++;
+
       const Int_t sathi = (Int_t)pix.GetNumHiGainSaturated();
       const Int_t satlo = (Int_t)pix.GetNumLoGainSaturated();
@@ -615 +667 @@
   *fLog << endl;
 
-  if (fHiGainOverFlow)
-    *fLog << warn << GetDescriptor()
-          << ": Histogram Overflow has occurred " << fHiGainOverFlow << " in the High-Gain! " << endl;
-  if (fLoGainOverFlow)
-    *fLog << warn << GetDescriptor()
-          << ": Histogram Overflow has occurred " << fLoGainOverFlow << " in the Low-Gain! " << endl;
-
   for (Int_t i=0; i<fHiGainArray->GetSize(); i++)
     {
@@ -639 +684 @@
         }
 
+      if (fHiGainOverFlow[i])
+        {
+          *fLog << warn << GetDescriptor()
+                << ": Hi-Gain Histogram Overflow occurred " << fHiGainOverFlow[i] 
+                << " times in pixel: " << i << " (without saturation!) " << endl;
+          bad.SetUncalibrated( MBadPixelsPix::kHiGainOverFlow ); 
+        }
+      
       FinalizeAbsTimes(histhi, pix, bad, fFirstHiGain, fLastHiGain);
     }
@@ -659 +712 @@
         }
   
+      if (fLoGainOverFlow[i])
+        {
+          *fLog << warn << GetDescriptor()
+                << ": Lo-Gain Histogram Overflow occurred " << fLoGainOverFlow[i] 
+                << " times in pixel: " << i << " (without saturation!) " << endl;
+          bad.SetUncalibrated( MBadPixelsPix::kLoGainOverFlow ); 
+        }
+
       MCalibrationChargePix &pix    = (MCalibrationChargePix&)(*fCam)[i];