Changeset 5609 for trunk/MagicSoft/Mars/mhcalib/MHCalibrationRelTimeCam.cc

Timestamp:

12/16/04 17:32:32 (20 years ago)

Author:

gaug

Message:

*** empty log message ***

File:

• trunk/MagicSoft/Mars/mhcalib/MHCalibrationRelTimeCam.cc (modified) (6 diffs)

trunk/MagicSoft/Mars/mhcalib/MHCalibrationRelTimeCam.cc

@@ -413 +413 @@
   MBadPixelsCam   *badcam = fIntensBad ? fIntensBad->GetCam() : fBadPixels;
 
+  const Int_t nareas = fAverageHiGainAreas->GetSize();
+  const Int_t nsectors = fAverageHiGainSectors->GetSize();
+
+  TArrayI satarea(nareas);
+  TArrayI satsect(nsectors);
+  fNumareahi  .Reset(); 
+  fNumsectorhi.Reset(); 
+
   for (Int_t i=0; i<fHiGainArray->GetSize(); i++)
     {
       
       MHCalibrationPix &histhi = (*this)[i];
-      
+
       if (histhi.IsExcluded())
         continue;
       
+      const Int_t aidx   = (*fGeom)[i].GetAidx();
+      const Int_t sector = (*fGeom)[i].GetSector();
+
       MCalibrationRelTimePix &pix = (MCalibrationRelTimePix&)(*relcam)[i] ;
 
+      fNumareahi[aidx]++;
+      fNumsectorhi[sector]++;
+      //
+      // Check saturation
+      // 
       if (histhi.GetSaturated() > fNumHiGainSaturationLimit*histhi.GetHGausHist()->GetEntries())
+        {
+          pix.SetHiGainSaturation();
+          histhi.SetExcluded();
+          satarea[aidx]++;
+          satsect[sector]++;
+        }
+      else
+        if (IsLoGain())
+          (*this)(i).SetExcluded();
+
+      //
+      // Check histogram overflow
+      // 
+      CheckOverflow(histhi);
+      if (IsLoGain())
+        CheckOverflow((*this)(i));
+
+    }
+
+  for (Int_t j=0; j<nareas; j++)
+    {
+      
+      MHCalibrationPix     &histhi = GetAverageHiGainArea(j);      
+      MCalibrationRelTimePix &pix = (MCalibrationRelTimePix&)relcam->GetAverageArea(j);
+
+      if (satarea[j] > 0.5*fNumareahi[j])
         {
           pix.SetHiGainSaturation();
@@ -430 +472 @@
       else
         if (IsLoGain())
-          (*this)(i).SetExcluded();
-
-      Stat_t overflow = histhi.GetHGausHist()->GetBinContent(histhi.GetHGausHist()->GetNbinsX()+1);
-      if (overflow > 0.1)
-        {
-          *fLog << warn << "HiGain Hist-overflow occurred " << overflow
-                << " times in pix: " << i << " (w/o saturation!) " << endl;
-          //          bad.SetUncalibrated( MBadPixelsPix::kHiGainOverFlow ); 
-        }
-
-      overflow = histhi.GetHGausHist()->GetBinContent(0);
-      if (overflow > 0.1)
-        {
-          *fLog << warn << "HiGain Hist-underflow occurred " << overflow
-                << " times in pix: " << i << " (w/o saturation!) " << endl;
-          //          bad.SetUncalibrated( MBadPixelsPix::kHiGainOverFlow ); 
-        }
-    }
-
-  for (Int_t j=0; j<fAverageHiGainAreas->GetSize(); j++)
-    {
-      
-      MHCalibrationPix &histhi = GetAverageHiGainArea(j);      
-
-      if (histhi.GetSaturated() > fNumHiGainSaturationLimit*histhi.GetHGausHist()->GetEntries())
-        {
-          MCalibrationRelTimePix &pix = (MCalibrationRelTimePix&)relcam->GetAverageArea(j);
-          pix.SetHiGainSaturation();
-          histhi.SetExcluded();
-        }
-      else
-        if (IsLoGain())
           GetAverageLoGainArea(j).SetExcluded();
 
+      //
+      // Check histogram overflow
+      // 
+      CheckOverflow(histhi);
+      if (IsLoGain())
+        CheckOverflow(GetAverageLoGainArea(j));
    }
   
@@ -472 +488 @@
       MCalibrationRelTimePix  &pix = (MCalibrationRelTimePix&)relcam->GetAverageSector(j) ;
 
-      if (histhi.GetSaturated() > fNumHiGainSaturationLimit*histhi.GetHGausHist()->GetEntries())
+      if (satsect[j] > 0.5*fNumsectorhi[j])
         {
           pix.SetHiGainSaturation();
@@ -480 +496 @@
         if (IsLoGain())        
           GetAverageLoGainSector(j).SetExcluded();
+
+      //
+      // Check histogram overflow
+      // 
+      CheckOverflow(histhi);
+      if (IsLoGain())
+        CheckOverflow(GetAverageLoGainSector(j));
     }
 
@@ -632 +655 @@
        // Ask for Hi-Gain saturation
        //
-       if (hipix.GetSaturated() > fNumHiGainSaturationLimit*hipix.GetHGausHist()->GetEntries() && IsLoGain())
-         {
+       if (hipix.IsExcluded() && IsLoGain())
+        {
            MHCalibrationPix &lopix = GetAverageLoGainArea(i);
            DrawDataCheckPixel(lopix,0.);
@@ -641 +664 @@
     }      
 }
+
+void MHCalibrationRelTimeCam::CheckOverflow( MHCalibrationPix &pix )
+{
+
+  if (pix.IsExcluded())
+    return;
+
+  TH1F *hist = pix.GetHGausHist();
+  
+  Stat_t overflow = hist->GetBinContent(hist->GetNbinsX()+1);
+  if (overflow > 0.0005*hist->GetEntries())
+    {
+      *fLog << warn << "HiGain Hist-overflow occurred " << overflow
+            << " times in " << pix.GetName() << " (w/o saturation!) " << endl;
+    }
+  
+  overflow = hist->GetBinContent(0);
+  if (overflow > 0.0005*hist->GetEntries())
+    {
+      *fLog << warn << "HiGain Hist-underflow occurred " << overflow
+            << " times in " << pix.GetName() << " (w/o saturation!) " << endl;
+    }
+}
+
 
 // -----------------------------------------------------------------------------