Changeset 3029 for trunk/MagicSoft/Mars/mcalib/MCalibrationCalc.cc

Timestamp:

02/05/04 17:28:52 (21 years ago)

Author:

gaug

Message:

*** empty log message ***

File:

• trunk/MagicSoft/Mars/mcalib/MCalibrationCalc.cc (modified) (10 diffs)

trunk/MagicSoft/Mars/mcalib/MCalibrationCalc.cc

@@ -68 +68 @@
 //                fits)
 //
+//                Hi-Gain vs. Lo-Gain Calibration (very memory-intensive)
+//                can be skipped with the command:
+//                MalibrationCam::SkipHiLoGainCalibration()
+//
 //  Input Containers:
 //   MRawEvtData
@@ -116 +120 @@
 const Int_t MCalibrationCalc::fBlindPixelId = 559;
 const Int_t MCalibrationCalc::fPINDiodeId   = 9999;
+const Byte_t MCalibrationCalc::fgSaturationLimit = 254;
+const Byte_t MCalibrationCalc::fgBlindPixelFirst = 3;
+const Byte_t MCalibrationCalc::fgBlindPixelLast  = 12;
 
 // --------------------------------------------------------------------------
@@ -144 +151 @@
     SETBIT(fFlags, kUseCosmicsRejection);
     SETBIT(fFlags, kUseQualityChecks);
+    SETBIT(fFlags, kHiLoGainCalibration);
 
     CLRBIT(fFlags, kBlindPixelOverFlow);
@@ -496 +504 @@
                 }
             }
-        }
+        }  /* if Use Times */
       
       switch(pixid)
@@ -502 +510 @@
           
         case fBlindPixelId:
-          
-          if (!blindpixel.FillCharge(sumhi)) 
-            *fLog << warn << 
-              "Overflow or Underflow occurred filling Blind Pixel sum = " << sumhi << endl;
-          
-          if (TESTBIT(fFlags,kUseTimes))
+
+          if (TESTBIT(fFlags,kUseBlindPixelFit))
             {
-              if (!blindpixel.FillTime(reltime))
+          
+              Float_t blindpixelsum = 0.;
+              //
+              // We need a dedicated signal extractor for the blind pixel
+              //
+              MPedestalPix &ped  = (*fPedestals)[pixid];
+              if (!CalcSignalBlindPixel(pixel.GetHiGainSamples(), blindpixelsum, ped.GetPedestal()))
+                return kFALSE;
+              
+              if (!blindpixel.FillCharge(blindpixelsum)) 
                 *fLog << warn << 
-                  "Overflow or Underflow occurred filling Blind Pixel time = " << reltime << endl;
-            }
-          
-          if (!TESTBIT(fFlags,kBlindPixelOverFlow))
-            if (!blindpixel.FillRChargevsTime(sumhi,fEvents))
-              {
+                  "Overflow or Underflow occurred filling Blind Pixel sum = " << blindpixelsum << endl;
+              
+              if (TESTBIT(fFlags,kUseTimes))
+                {
+                  if (!blindpixel.FillTime(reltime))
+                    *fLog << warn << 
+                      "Overflow or Underflow occurred filling Blind Pixel time = " << reltime << endl;
+                }
+              
+              if (!TESTBIT(fFlags,kBlindPixelOverFlow))
+                if (!blindpixel.FillRChargevsTime(blindpixelsum,fEvents))
+                  {
+                    *fLog << warn << 
+                      "Overflow or Underflow occurred filling Blind Pixel eventnr = " << fEvents << endl;
+                    SETBIT(fFlags,kBlindPixelOverFlow);
+                  }
+            } /* if use blind pixel */
+          
+          break;
+              
+        case fPINDiodeId:
+
+          if (TESTBIT(fFlags,kUsePinDiodeFit))
+            {
+          
+              if (!pindiode.FillCharge(sumhi)) 
                 *fLog << warn << 
-                  "Overflow or Underflow occurred filling Blind Pixel eventnr = " << fEvents << endl;
-                SETBIT(fFlags,kBlindPixelOverFlow);
-              }
+                  "Overflow or Underflow occurred filling PINDiode: sum = " << sumhi << endl;
+              
+              if (TESTBIT(fFlags,kUseTimes))
+                {
+                  if (!pindiode.FillAbsTime(abstime))
+                    *fLog << warn << 
+                      "Overflow or Underflow occurred filling PINDiode abs. time = " << abstime << endl;
+                  if (!pindiode.FillRelTime(reltime))
+                    *fLog << warn << 
+                      "Overflow or Underflow occurred filling PINDiode rel. time = " << reltime << endl;
+                }
+              
+              if (!TESTBIT(fFlags,kPINDiodeOverFlow))
+                if (!pindiode.FillRChargevsTime(sumhi,fEvents))
+                  {
+                    *fLog << warn 
+                          << "Overflow or Underflow occurred filling PINDiode: eventnr = " 
+                          << fEvents << endl;
+                    SETBIT(fFlags,kPINDiodeOverFlow);
+                  }
+              
+            } /* if use PIN Diode */
           
           break;
-          
-        case fPINDiodeId:
-
-          if (!pindiode.FillCharge(sumhi)) 
-            *fLog << warn << 
-              "Overflow or Underflow occurred filling PINDiode: sum = " << sumhi << endl;
-
-          if (TESTBIT(fFlags,kUseTimes))
-            {
-              if (!pindiode.FillAbsTime(abstime))
-                *fLog << warn << 
-                  "Overflow or Underflow occurred filling PINDiode abs. time = " << abstime << endl;
-              if (!pindiode.FillRelTime(reltime))
-                *fLog << warn << 
-                  "Overflow or Underflow occurred filling PINDiode rel. time = " << reltime << endl;
-            }
-          
-          if (!TESTBIT(fFlags,kPINDiodeOverFlow))
-            if (!pindiode.FillRChargevsTime(sumhi,fEvents))
-              {
-                *fLog << warn 
-                      << "Overflow or Underflow occurred filling PINDiode: eventnr = " 
-                      << fEvents << endl;
-                SETBIT(fFlags,kPINDiodeOverFlow);
-              }
-          break;
-          
+              
         default:
-          
-          pix.FillChargesInGraph(sumhi,sumlo);
-
+              
           if (!TESTBIT(fFlags,kLoGainOverFlow))
             if (!pix.FillRChargevsTimeLoGain(sumlo,fEvents))
@@ -561 +585 @@
                       << fEvents << endl;
                 SETBIT(fFlags,kLoGainOverFlow);
+                SkipHiLoGainCalibration();
               }
           
@@ -570 +595 @@
                       << fEvents << endl;
                 SETBIT(fFlags,kHiGainOverFlow);
+                SkipHiLoGainCalibration();
               }
 
-          
+
+          if (TESTBIT(fFlags,kHiLoGainCalibration))
+            pix.FillChargesInGraph(sumhi,sumlo);
+
           if (sig.IsLoGainUsed())
             {
@@ -670 +699 @@
           // Set the corresponding values
           //
-          const Float_t nslices     = (Float_t)fSignals->GetNumUsedFADCSlices();
+          const Float_t nslices     = (Float_t)(fgBlindPixelLast-fgBlindPixelFirst+1);
+          const Float_t sqrslice    = TMath::Sqrt(nslices);
           const ULong_t nentries    = fPedestals->GetTotalEntries();
           
           const Float_t peddiff     = (pedhist.GetChargeMean()-pedpix.GetPedestal())*nslices;
-
-          const Float_t sqrslice    = TMath::Sqrt(nslices);
 
           Float_t pederr  = pedhist.GetChargeMeanErr()*pedhist.GetChargeMeanErr();
@@ -721 +749 @@
 
       //
-      // get the pedestals
-      //
-      const Float_t ped    = (*fPedestals)[pixid].GetPedestal() * fNumHiGainSamples;
-      const Float_t prms   = (*fPedestals)[pixid].GetPedestalRms() * fSqrtHiGainSamples;
-
-      //
-      // set them in the calibration camera
-      //
-      pix.SetPedestal(ped,prms);
-
-
-      //
       // Check if the pixel has been excluded from the fits
       //
       if (pix.IsExcluded())
         continue;
+
+      //
+      // get the pedestals
+      //
+      const Float_t ped    = (*fPedestals)[pixid].GetPedestal();
+      const Float_t prms   = (*fPedestals)[pixid].GetPedestalRms();
+
+      //
+      // set them in the calibration camera
+      //
+      pix.SetPedestal(ped,prms,(Float_t)fNumHiGainSamples,(Float_t)fNumLoGainSamples);
 
       //
@@ -793 +820 @@
 }
 
+Bool_t MCalibrationCalc::CalcSignalBlindPixel(Byte_t *ptr, Float_t &signal, const Float_t ped) const
+{
+
+  Byte_t *start = ptr + fgBlindPixelFirst;
+  Byte_t *end   = ptr + fgBlindPixelLast;
+
+  Byte_t sum = 0;
+  Int_t  sat = 0;
+  
+  ptr = start;
+
+  while (ptr<end)
+    {
+      sum += *ptr;
+      
+      if (*ptr++ >= fgSaturationLimit)
+            sat++;
+    }
+
+  if (sat)
+    {
+      *fLog << err << "HI Gain Saturation occurred in the blind pixel! " 
+            << " Do not know yet how to treat this ... aborting " << endl;
+      return kFALSE;
+    }
+
+  signal = (Float_t)sum - ped*(Float_t)(fgBlindPixelLast-fgBlindPixelFirst+1);
+
+  return kTRUE;
+}