Changeset 3923 for trunk/MagicSoft/Mars/mcalib

Timestamp:

05/01/04 09:57:15 (20 years ago)

Author:

gaug

Message:

*** empty log message ***

Location:

trunk/MagicSoft/Mars/mcalib

Files:

• MCalibrationChargeCalc.cc (modified) (6 diffs)
• MHCalibrationRelTimePix.cc (modified) (1 diff)

trunk/MagicSoft/Mars/mcalib/MCalibrationChargeCalc.cc

@@ -229 +229 @@
 const Float_t MCalibrationChargeCalc::fgLambdaCheckLimit       = 0.2;
 const Float_t MCalibrationChargeCalc::fgPheErrLimit            = 4.;
+const Float_t MCalibrationChargeCalc::fgFFactorErrLimit        = 3.;
 // --------------------------------------------------------------------------
 //
@@ -245 +246 @@
 // - fChargeErrLimit    to fgChargeErrLimit   
 // - fChargeRelErrLimit to fgChargeRelErrLimit
+// - fFFactorErrLimit   to fgFFactorErrLimit
 // - fLambdaCheckLimit  to fgLambdaCheckLimit 
 // - fLambdaErrLimit    to fgLambdaErrLimit 
@@ -256 +258 @@
     : fBadPixels(NULL), fCam(NULL), fBlindPixel(NULL), fPINDiode(NULL), 
       fQECam(NULL), fGeom(NULL), fPedestals(NULL), fEvtTime(NULL)
-      {
-
+{
+        
   fName  = name  ? name  : "MCalibrationChargeCalc";
   fTitle = title ? title : "Task to calculate the calibration constants and MCalibrationCam ";
@@ -269 +271 @@
   SetChargeErrLimit();  
   SetChargeRelErrLimit();
+  SetFFactorErrLimit();
   SetLambdaCheckLimit();
   SetLambdaErrLimit();
@@ -1207 +1210 @@
                             + fQECam->GetPlexiglassQERelVar();
 
+  const UInt_t nareas   = fGeom->GetNumAreas();
+
+  Float_t lowlim           [nareas];
+  Float_t upplim           [nareas];
+  Float_t avffactorphotons [nareas];
+  Float_t avffactorphotvar [nareas];
+  Int_t   numffactor       [nareas];
+
+  memset(lowlim          ,0, nareas   * sizeof(Float_t));
+  memset(upplim          ,0, nareas   * sizeof(Float_t));
+  memset(avffactorphotons,0, nareas   * sizeof(Float_t));
+  memset(avffactorphotvar,0, nareas   * sizeof(Float_t));
+  memset(numffactor      ,0, nareas   * sizeof(Int_t));
+
   const UInt_t npixels  = fGeom->GetNumPixels();
 
@@ -1240 +1257 @@
         *fLog << warn << GetDescriptor() 
               << ": Cannot update Quantum efficiencies with the F-Factor Method" << endl;
+
+      const Int_t aidx = (*fGeom)[i].GetAidx();
+
+      avffactorphotons[aidx] += pix.GetMeanFFactorFADC2Phot();
+      avffactorphotvar[aidx] += pix.GetMeanFFactorFADC2PhotVar();
+      numffactor[aidx]++;
+    }
+
+  for (UInt_t i=0; i<nareas; i++)
+    {
+      avffactorphotons[i] /= numffactor[i];
+      avffactorphotvar[i] /= numffactor[i];
+      lowlim  [i] = avffactorphotons[i] - fFFactorErrLimit*TMath::Sqrt(avffactorphotvar[i]);
+      upplim  [i] = avffactorphotons[i] + fFFactorErrLimit*TMath::Sqrt(avffactorphotvar[i]);
+    }
+  
+  for (UInt_t i=0; i<npixels; i++)
+    {
+      
+      MCalibrationChargePix &pix = (MCalibrationChargePix&)(*fCam)[i];
+
+      if (!pix.IsFFactorMethodValid())
+        continue;
+      
+      const Float_t ffactor = pix.GetMeanFFactorFADC2Phot();
+      MBadPixelsPix    &bad = (*fBadPixels)[i];
+
+      const Int_t   aidx   = (*fGeom)[i].GetAidx();
+
+      if ( ffactor < lowlim[aidx] || ffactor > upplim[aidx] )
+        {
+          *fLog << warn << GetDescriptor() << ": Deviating F-Factor: " 
+                << Form("%4.2f",ffactor) << " out of accepted limits: [" 
+                << Form("%4.2f%s%4.2f",lowlim[aidx],",",upplim[aidx]) << "] in pixel " << i << endl;
+          bad.SetUncalibrated( MBadPixelsPix::kDeviatingFFactor );
+          bad.SetUnsuitable  ( MBadPixelsPix::kUnsuitableRun    );
+          pix.SetExcluded();
+        }
     }
 }

trunk/MagicSoft/Mars/mcalib/MHCalibrationRelTimePix.cc

@@ -41 +41 @@
 //
 const Int_t   MHCalibrationRelTimePix::fgRelTimeNbins    = 900;
-const Axis_t  MHCalibrationRelTimePix::fgRelTimeFirst    = -13.;
-const Axis_t  MHCalibrationRelTimePix::fgRelTimeLast     =  13.;
+const Axis_t  MHCalibrationRelTimePix::fgRelTimeFirst    = -5.;
+const Axis_t  MHCalibrationRelTimePix::fgRelTimeLast     =  5.;
 const Float_t MHCalibrationRelTimePix::fgFADCSliceWidth  =  3.3333;
 // --------------------------------------------------------------------------