Changeset 5773 for trunk

Timestamp:

01/10/05 17:29:37 (20 years ago)

Author:

gaug

Message:

*** empty log message ***

Location:

trunk/MagicSoft/Mars

Files:

• Changelog (modified) (1 diff)
• mhcalib/MHCalibrationChargeCam.cc (modified) (4 diffs)
• mhcalib/MHCalibrationChargeCam.h (modified) (3 diffs)

trunk/MagicSoft/Mars/Changelog

@@ -41 +41 @@
    * msignal/MExtractor.h
      - added Getter for fOffsetLoGain (was still missing)
+
+   * mhcalib/MHCalibrationChargeCam.[h,cc]
+     - remove pixels with histogram overflow only if their percentage 
+       passes a certain limit (default: 0.5%). Before, every overflow 
+       lead to an unsuitable pixel.
 
 

trunk/MagicSoft/Mars/mhcalib/MHCalibrationChargeCam.cc

@@ -172 +172 @@
 const Float_t MHCalibrationChargeCam::fgNumHiGainSaturationLimit = 0.01;
 const Float_t MHCalibrationChargeCam::fgNumLoGainSaturationLimit = 0.005;
+const Float_t MHCalibrationChargeCam::fgNumOverflowLimit         = 0.005;
 const Float_t MHCalibrationChargeCam::fgTimeLowerLimit           = 1.;
 const Float_t MHCalibrationChargeCam::fgTimeUpperLimit           = 2.;
@@ -219 +220 @@
   SetNumHiGainSaturationLimit(fgNumHiGainSaturationLimit);
   SetNumLoGainSaturationLimit(fgNumLoGainSaturationLimit);
+  SetNumOverflowLimit(fgNumOverflowLimit);
+
   SetTimeLowerLimit();
   SetTimeUpperLimit();
@@ -925 +928 @@
 
       Stat_t overflow = h->GetBinContent(h->GetNbinsX()+1);
-      if (overflow > 0.1)
+      if (overflow > fNumOverflowLimit*histhi.GetHGausHist()->GetEntries())
         {
           *fLog << warn
                 << "HiGain Hist-overflow " << overflow
-                << " times in pix: " << i << " (w/o saturation!) " << endl;
+                << " times in: " << histhi.GetName() << " (w/o saturation!) " << endl;
           bad.SetUncalibrated( MBadPixelsPix::kHiGainOverFlow ); 
         }
 
       overflow = h->GetBinContent(0);
-      if (overflow > 0.1)
+      if (overflow >  fNumOverflowLimit*histhi.GetHGausHist()->GetEntries())
         {
           *fLog << warn
                 << "HiGain Hist-underflow " << overflow
-                << " times in pix: " << i << " (w/o saturation!) " << endl;
+                << " times in pix: " << histhi.GetName() << " (w/o saturation!) " << endl;
           bad.SetUncalibrated( MBadPixelsPix::kHiGainOverFlow ); 
         }
@@ -967 +970 @@
 
         Stat_t overflow = h->GetBinContent(h->GetNbinsX()+1);
-        if (overflow > 0.1)
+        if (overflow >  fNumOverflowLimit*histlo.GetHGausHist()->GetEntries())
           {
               *fLog << warn
                   << "LoGain Hist-overflow " << overflow
-                  << " times in pix: " << i << " (w/o saturation!) " << endl;
+                  << " times in: " << histlo.GetName() << " (w/o saturation!) " << endl;
             bad.SetUncalibrated( MBadPixelsPix::kLoGainOverFlow ); 
           }
 
         overflow = h->GetBinContent(0);
-        if (overflow > 0.1)
+        if (overflow > fNumOverflowLimit*histlo.GetHGausHist()->GetEntries())
           {
             *fLog << warn
                   << "LoGain Hist-underflow " << overflow
-                  << " times in pix: " << i << " (w/o saturation!) " << endl;
+                  << " times in: " << histlo.GetName() << " (w/o saturation!) " << endl;
             bad.SetUncalibrated( MBadPixelsPix::kLoGainOverFlow ); 
           }

trunk/MagicSoft/Mars/mhcalib/MHCalibrationChargeCam.h

@@ -44 +44 @@
   static const Float_t fgNumHiGainSaturationLimit;   //! The default for fNumHiGainSaturationLimit (now at: 0.01)
   static const Float_t fgNumLoGainSaturationLimit;   //! The default for fNumLoGainSaturationLimit (now at: 0.005)
+  static const Float_t fgNumOverflowLimit;           //! The default for fNumOverflowLimit         (now at: 0.005)
+
   static const Float_t fgTimeLowerLimit;             //! Default for fTimeLowerLimit    (now set to: 1.) 
   static const Float_t fgTimeUpperLimit;             //! Default for fTimeUpperLimit    (now set to: 2.)
@@ -56 +58 @@
   Axis_t  fLoGainLast;                               // Upper histogram limit low gain  
 
+  Float_t fNumOverflowLimit;                         // Part of maximum allowed overflow events
+  
   TString fAbsHistName;                              // Histogram names abs.times
   TString fAbsHistTitle;                             // Histogram titles abs. times
@@ -113 +117 @@
   void SetLoGainLast        ( const Axis_t f )       { fLoGainLast    = f; }  
 
+  void SetNumOverflowLimit ( const Float_t f=fgNumOverflowLimit ) { fNumOverflowLimit = f; }
+  
   void SetTimeLowerLimit    ( const Float_t f=fgTimeLowerLimit ) { fTimeLowerLimit = f; }
   void SetTimeUpperLimit    ( const Float_t f=fgTimeUpperLimit ) { fTimeUpperLimit = f; }