Changeset 7189 for trunk/MagicSoft/Mars/mcalib/MCalibrationChargeCalc.cc

Timestamp:

07/14/05 16:03:32 (19 years ago)

Author:

tbretz

Message:

*** empty log message ***

File:

• trunk/MagicSoft/Mars/mcalib/MCalibrationChargeCalc.cc (modified) (28 diffs)

trunk/MagicSoft/Mars/mcalib/MCalibrationChargeCalc.cc

@@ -184 +184 @@
 //   + Float_t fPheErrUpperLimit; // Upper limit acceptance nr. phe's w.r.t. area idx mean (in sigmas)
 //
+//
 //  ClassVersion 3:
 //   + Bool_t  fUseExtractorRes;  // Include extractor resolution in F-Factor method
+//
+//  ClassVersion 3:
+//   +  Float_t fUnsuitablesLimit;                   // Limit for relative number of unsuitable pixels
+//   +  Float_t fUnreliablesLimit;                   // Limit for relative number of unreliable pixels
+//   +  Float_t fExternalNumPhes;                    // External mean number of photo-electrons set from outside
+//   +  Float_t fExternalNumPhesRelVar;              // External rel. var. number of photo-electrons set from outside  
 //
 //////////////////////////////////////////////////////////////////////////////
@@ -198 +205 @@
 
 #include "MParList.h"
-
-#include "MStatusDisplay.h"
 
 #include "MCalibrationPattern.h"
@@ -247 +252 @@
 const Float_t MCalibrationChargeCalc::fgFFactorErrLimit        = 4.5;
 const Float_t MCalibrationChargeCalc::fgArrTimeRmsLimit        = 3.5;
+const Float_t MCalibrationChargeCalc::fgUnsuitablesLimit       = 0.1;
+const Float_t MCalibrationChargeCalc::fgUnreliablesLimit       = 0.3;
 const TString MCalibrationChargeCalc::fgNamePedestalCam = "MPedestalCam";
 
@@ -263 +270 @@
 //
 // Initializes:
-// - fArrTimeRmsLimit   to fgArrTimeRmsLimit
-// - fChargeLimit       to fgChargeLimit 
-// - fChargeErrLimit    to fgChargeErrLimit   
-// - fChargeRelErrLimit to fgChargeRelErrLimit
-// - fFFactorErrLimit   to fgFFactorErrLimit
-// - fLambdaCheckLimit  to fgLambdaCheckLimit 
-// - fLambdaErrLimit    to fgLambdaErrLimit 
-// - fNamePedestalCam   to fgNamePedestalCam
-// - fPheErrLowerLimit  to fgPheErrLowerLimit
-// - fPheErrUpperLimit  to fgPheErrUpperLimit
-// - fPulserColor       to MCalibrationCam::kCT1
-// - fOutputPath        to "."
-// - fOutputFile        to "ChargeCalibStat.txt"
-// - flag debug         to kFALSE
+// - fArrTimeRmsLimit     to fgArrTimeRmsLimit
+// - fChargeLimit         to fgChargeLimit
+// - fChargeErrLimit      to fgChargeErrLimit
+// - fChargeRelErrLimit   to fgChargeRelErrLimit
+// - fFFactorErrLimit     to fgFFactorErrLimit
+// - fLambdaCheckLimit    to fgLambdaCheckLimit
+// - fLambdaErrLimit      to fgLambdaErrLimit
+// - fNamePedestalCam     to fgNamePedestalCam
+// - fPheErrLowerLimit    to fgPheErrLowerLimit
+// - fPheErrUpperLimit    to fgPheErrUpperLimit
+// - fPulserColor         to MCalibrationCam::kCT1
+// - fOutputPath          to "."
+// - fOutputFile          to "ChargeCalibStat.txt"
+// - flag debug           to kFALSE
+// - IsUseExtractorRes    to kFALSE
+// - IsUseExternalNumPhes to kFALSE
+// - fExternalNumPhes     to 0.
+// - fExternalNumPhesRelVar to 0.
 //
 // Sets all checks
@@ -284 +295 @@
 //
 MCalibrationChargeCalc::MCalibrationChargeCalc(const char *name, const char *title)
-    : fUseExtractorRes(kFALSE),
-    fGeom(NULL), fSignal(NULL), fCalibPattern(NULL), fExtractor(NULL)
+    : fGeom(NULL), fSignal(NULL), fCalibPattern(NULL), fExtractor(NULL)
 {
         
@@ -309 +319 @@
   SetPheErrLowerLimit  ();
   SetPheErrUpperLimit  ();
+  SetUnsuitablesLimit  ();
+  SetUnreliablesLimit  ();
+  SetUseExtractorRes   ();
+  SetUseUnreliables    ();
+  SetUseExternalNumPhes(kFALSE);
+ 
+  SetExternalNumPhes      ();
+  SetExternalNumPhesRelVar();
 
   SetCheckArrivalTimes     ();
@@ -560 +578 @@
     }
 
+  fResultFlags.Set(fGeom->GetNumAreas());
+
   return kTRUE;
 }
@@ -777 +797 @@
     {
 
-      const MPedestalPix     &ped = fPedestals->GetAverageArea(aidx);
+      const MPedestalPix    &ped = fPedestals->GetAverageArea(aidx);
       MCalibrationChargePix &pix = (MCalibrationChargePix&)chargecam->GetAverageArea(aidx);
+      const MArrayI &arr         = fHCam->GetAverageAreaNum();
 
       FinalizePedestals(ped,pix,aidx);
+
+      //
+      // Correct for sqrt(number of valid pixels) in the pedestal RMS
+      // (already done for calibration sigma in MHCalibrationCam::CalcAverageSigma()
+      //
+      pix.SetPedRMS(pix.GetPedRms()*TMath::Sqrt((Float_t)arr[aidx]),
+                    pix.GetPedRmsErr()*TMath::Sqrt((Float_t)arr[aidx]));
+      pix.SetSigma (pix.GetSigma()/pix.GetFFactorFADC2Phe());
+
       FinalizeCharges(pix, chargecam->GetAverageBadArea(aidx),"area id");
       FinalizeArrivalTimes(pix, chargecam->GetAverageBadArea(aidx), "area id");
@@ -847 +877 @@
   }
 
-  //
-  // Finalize calibration statistics
-  //
-  FinalizeUnsuitablePixels();
-
-  chargecam->SetReadyToSave();
-  qecam    ->SetReadyToSave();
-  badcam   ->SetReadyToSave();
-  
-  if (blindcam)
-    blindcam->SetReadyToSave();
-  if (fPINDiode)
-    fPINDiode->SetReadyToSave();
-  
   *fLog << inf << endl;
   *fLog << GetDescriptor() << ": Fatal errors statistics:" << endl;  
@@ -884 +900 @@
   PrintUncalibrated(MBadPixelsPix::kDeviatingNumPhes,
                     "Deviating number of phes:                         ");
+  PrintUncalibrated(MBadPixelsPix::kLoGainBlackout,
+                    "Too many blackout events in low gain:             ");
   PrintUncalibrated(MBadPixelsPix::kPreviouslyExcluded,
                     "Previously excluded:                              ");
@@ -906 +924 @@
       SetLogStream(oldlog);
 
-  return kTRUE;
+  chargecam->SetReadyToSave();
+  qecam    ->SetReadyToSave();
+  badcam   ->SetReadyToSave();
+
+  if (blindcam)
+    blindcam->SetReadyToSave();
+  if (fPINDiode)
+    fPINDiode->SetReadyToSave();
+
+  //
+  // Finalize calibration statistics
+  //
+  return FinalizeUnsuitablePixels();
 }
 
@@ -1042 +1072 @@
     }
 
-  if (!fUseExtractorRes)
+  if (!IsUseExtractorRes())
     return kTRUE;
 
@@ -1059 +1089 @@
     : resinphes/cal.GetMeanConvFADC2Phe();
 
-  if (resinfadc > 1.5*cal.GetPedRms() )
-    {
-      *fLog << warn << " Extractor Resolution: " << resinfadc << " bigger than Pedestal RMS " << cal.GetPedRms() << endl;
-      resinfadc = cal.GetPedRms();
+  if (resinfadc > 3.0*cal.GetPedRms() )
+    {
+      *fLog << warn << " Extractor Resolution: " << resinfadc << " bigger than 3 Pedestal RMS "
+            << cal.GetPedRms() << endl;
+      resinfadc = 3.0*cal.GetPedRms();
     }
 
@@ -1152 +1183 @@
           if (bad.IsUncalibrated( MBadPixelsPix::kChargeIsPedestal))
             bad.SetUnsuitable(   MBadPixelsPix::kUnsuitableRun   );
- 
-          if (bad.IsUncalibrated( MBadPixelsPix::kChargeErrNotValid ))
-            bad.SetUnsuitable(   MBadPixelsPix::kUnsuitableRun    );
-
-          if (bad.IsUncalibrated( MBadPixelsPix::kChargeRelErrNotValid ))
-            bad.SetUnsuitable(   MBadPixelsPix::kUnsuitableRun    );
         }
       
@@ -1268 +1293 @@
       if (bad.IsUncalibrated(MBadPixelsPix::kChargeSigmaNotValid))
         continue;
-      
+
+      if (!IsUseUnreliables())
+        if (bad.IsUnsuitable(MBadPixelsPix::kUnreliableRun))
+          continue;
+
       const Float_t nphe  = pix.GetPheFFactorMethod();
       const Int_t   aidx  = (*fGeom)[i].GetAidx();
@@ -1322 +1351 @@
       if (ndf < 5)
         {
-          *fLog << warn << GetDescriptor() << ": Cannot use a Gauss fit to the number of photo-electrons " 
-                << "in the camera with area index: " << i << endl;
-          *fLog << warn << GetDescriptor() << ": Number of dof.: " << ndf << " is smaller than 5 " << endl;
-          *fLog << warn << GetDescriptor() << ": Will use the simple mean and rms " << endl;
-          delete hist;
-          continue;
+            *fLog << warn << GetDescriptor() << ": Cannot use a Gauss fit to the number of photo-electrons " << endl;
+            *fLog << " in the camera with area index: " << i << endl;
+            *fLog << " Number of dof.: " << ndf << " is smaller than 5 " << endl;
+            *fLog << " Will use the simple mean and rms " << endl;
+            delete hist;
+            SetPheFitOK(i,kFALSE);
+            continue;
         }
       
@@ -1334 +1364 @@
       if (prob < 0.001)
         {
-          *fLog << warn << GetDescriptor() << ": Cannot use a Gauss fit to the number of photo-electrons " 
-                << "in the camera with area index: " << i << endl;
-          *fLog << warn << GetDescriptor() << ": Fit probability " << prob 
-                << " is smaller than 0.001 " << endl;
-          *fLog << warn << GetDescriptor() << ": Will use the simple mean and rms " << endl;
-          delete hist;
-          continue;
+            *fLog << warn << GetDescriptor() << ": Cannot use a Gauss fit to the number of photo-electrons " << endl;
+            *fLog << " in the camera with area index: " << i << endl;
+            *fLog << " Fit probability " << prob << " is smaller than 0.001 " << endl;
+            *fLog << " Will use the simple mean and rms " << endl;
+            delete hist;
+            SetPheFitOK(i,kFALSE);
+            continue;
         }
       
       if (mean < 0.)
         {
-          *fLog << inf << GetDescriptor() << ": Fitted mean number of photo-electrons " 
-                << "with area idx " << i << ": " << mean << " is smaller than 0. " << endl;
-          *fLog << warn << GetDescriptor() << ": Will use the simple mean and rms " << endl;          
-          delete hist;
-          continue;
+            *fLog << inf << GetDescriptor() << ": Fitted mean number of photo-electrons " << endl;
+            *fLog << " with area idx " << i << ": " << mean << " is smaller than 0. " << endl;
+            *fLog << warn << " Will use the simple mean and rms " << endl;
+            SetPheFitOK(i,kFALSE);
+            delete hist;
+            continue;
         }
       
@@ -1359 +1390 @@
 
       delete hist;
+
+      SetPheFitOK(i,kTRUE);
     }
 
@@ -1383 +1416 @@
       if (bad.IsUncalibrated(MBadPixelsPix::kChargeSigmaNotValid))
         continue;
+
+      if (!IsUseUnreliables())
+        if (bad.IsUnsuitable(MBadPixelsPix::kUnreliableRun))
+          continue;
 
       const Float_t nvar  = pix.GetPheFFactorMethodVar();
@@ -1730 +1767 @@
   MCalibrationQEPix     &qepix = (MCalibrationQEPix&)    qecam->GetAverageArea(0);
 
-  const Float_t avphotons   = avpix.GetPheFFactorMethod()       
+  if (IsDebug())
+    *fLog << dbginf << "External Phes: " << fExternalNumPhes
+          << " Internal Phes: " << avpix.GetPheFFactorMethod() << endl;
+
+  const Float_t avphotons   = ( IsUseExternalNumPhes()
+                                ? fExternalNumPhes
+                                : avpix.GetPheFFactorMethod() )
                            / qepix.GetDefaultQE(fPulserColor) 
                            / qepix.GetPMTCollectionEff()
@@ -1736 +1779 @@
                            / qecam->GetPlexiglassQE();
 
-  const Float_t avphotrelvar = avpix.GetPheFFactorMethodRelVar() 
+  const Float_t avphotrelvar = ( IsUseExternalNumPhes()
+                                 ? fExternalNumPhesRelVar
+                                 : avpix.GetPheFFactorMethodRelVar() )
                             + qepix.GetDefaultQERelVar(fPulserColor) / fNumInnerFFactorMethodUsed
                             + qepix.GetPMTCollectionEffRelVar()
@@ -1846 +1891 @@
       if (ndf < 2)
         {
-          *fLog << warn << GetDescriptor() << ": Cannot use a Gauss fit to the F-Factor " 
-                << "in the camera with area index: " << i << endl;
-          *fLog << "Number of dof.: " << ndf << " is smaller than 2 " << endl;
-          *fLog << "Will use the simple mean and rms." << endl;
+          *fLog << warn << GetDescriptor() << ": Cannot use a Gauss fit to the F-Factor " << endl;
+          *fLog << " in the camera with area index: " << i << endl;
+          *fLog << " Number of dof.: " << ndf << " is smaller than 2 " << endl;
+          *fLog << " Will use the simple mean and rms." << endl;
           delete hist;
+          SetFFactorFitOK(i,kFALSE);
           continue;
         }
@@ -1858 +1904 @@
       if (prob < 0.001)
         {
-          *fLog << warn << GetDescriptor() << ": Cannot use a Gauss fit to the F-Factor " 
-                << "in the camera with area index: " << i << endl;
-          *fLog << "Fit probability " << prob
-                << " is smaller than 0.001 " << endl;
-          *fLog << "Will use the simple mean and rms." << endl;
+          *fLog << warn << GetDescriptor() << ": Cannot use a Gauss fit to the F-Factor " << endl;
+          *fLog << " in the camera with area index: " << i << endl;
+          *fLog << " Fit probability " << prob << " is smaller than 0.001 " << endl;
+          *fLog << " Will use the simple mean and rms." << endl;
           delete hist;
-          continue;
+          SetFFactorFitOK(i,kFALSE);
+         continue;
         }
 
       *fLog << inf << GetDescriptor() << ": Mean F-Factor " 
-          << "with area index #" << i << ": "
-            << Form("%4.2f+-%4.2f",mean,sigma) << endl;
+          << "with area index #" << i << ": " << Form("%4.2f+-%4.2f",mean,sigma) << endl;
 
       lowlim  [i] = 1.;
@@ -1875 +1920 @@
 
       delete hist;
+
+      SetFFactorFitOK(i,kTRUE);
     }
   
@@ -2134 +2181 @@
 // - store numbers of bad pixels of each type in fCam or fIntensCam
 //
-void MCalibrationChargeCalc::FinalizeUnsuitablePixels()
+Bool_t MCalibrationChargeCalc::FinalizeUnsuitablePixels()
 {
   
@@ -2143 +2190 @@
   const Int_t nareas = fGeom->GetNumAreas();
 
+  TArrayI suit(nareas);
+  TArrayI unsuit(nareas);
+  TArrayI unrel(nareas);
+
+  const MBadPixelsCam *badcam = fIntensBad ? (MBadPixelsCam*)fIntensBad->GetCam() : fBadPixels;
+  MCalibrationChargeCam *chargecam = fIntensCam ? (MCalibrationChargeCam*)fIntensCam->GetCam() : fCam;
+
+  Int_t unsuitcnt=0;
+  Int_t unrelcnt =0;
+
+  // Count number of succesfully calibrated pixels
+  for (Int_t aidx=0; aidx<nareas; aidx++)
+  {
+      suit[aidx]   = badcam->GetNumSuitable(MBadPixelsPix::kUnsuitableRun,   fGeom, aidx);
+      unsuit[aidx] = badcam->GetNumUnsuitable(MBadPixelsPix::kUnsuitableRun, fGeom, aidx);
+      unrel[aidx]  = badcam->GetNumUnsuitable(MBadPixelsPix::kUnreliableRun, fGeom, aidx);
+
+      unsuitcnt   += unsuit[aidx];
+      unrelcnt    += unrel[aidx];
+
+      chargecam->SetNumUnsuitable(unsuit[aidx], aidx);
+      chargecam->SetNumUnreliable(unrel[aidx],  aidx);
+  }
+
   TArrayI counts(nareas);
-
-  MBadPixelsCam         *badcam    = fIntensBad 
-    ? (MBadPixelsCam*)fIntensBad->GetCam()  : fBadPixels;
-  MCalibrationChargeCam *chargecam = fIntensCam 
-    ? (MCalibrationChargeCam*)fIntensCam->GetCam()  : fCam;
-
-  for (Int_t i=0; i<badcam->GetSize(); i++)
-    {
-      MBadPixelsPix &bad = (*badcam)[i];
-      if (!bad.IsUnsuitable(MBadPixelsPix::kUnsuitableRun))
-        {
+  for (Int_t i=0; i<chargecam->GetSize(); i++)
+  {
+      MCalibrationPix &pix = (*chargecam)[i];
+      if (pix.IsHiGainSaturation())
+      {
           const Int_t aidx = (*fGeom)[i].GetAidx();
           counts[aidx]++;
-        }
-    }
+      }
+  }
 
   if (fGeom->InheritsFrom("MGeomCamMagic"))
-    *fLog << " " << setw(7) << "Successfully calibrated Pixels: " 
-          << Form("%s%3i%s%3i","Inner: ",counts[0]," Outer: ",counts[1]) << endl;
-  
-  counts.Reset();
-
-  for (Int_t i=0; i<badcam->GetSize(); i++)
-    {
-      MBadPixelsPix &bad =  (*badcam)[i];
-
-      if (bad.IsUnsuitable(MBadPixelsPix::kUnsuitableRun))
-        {
-          const Int_t aidx = (*fGeom)[i].GetAidx();
-          counts[aidx]++;
-        }
-    }
-
-  for (Int_t aidx=0; aidx<nareas; aidx++)
-    chargecam->SetNumUnsuitable(counts[aidx], aidx);
-
-  if (fGeom->InheritsFrom("MGeomCamMagic"))
-    *fLog << " " << setw(7) << "Uncalibrated Pixels:            " 
-          << Form("%s%3i%s%3i","Inner: ",counts[0]," Outer: ",counts[1]) << endl;
-
-  counts.Reset();
-
-  for (Int_t i=0; i<badcam->GetSize(); i++)
-    {
-
-      MBadPixelsPix &bad = (*badcam)[i];
-
-      if (bad.IsUnsuitable(MBadPixelsPix::kUnreliableRun))
-        {
-          const Int_t aidx = (*fGeom)[i].GetAidx();
-          counts[aidx]++;
-        }
-    }
-
-  for (Int_t aidx=0; aidx<nareas; aidx++)
-    chargecam->SetNumUnreliable(counts[aidx], aidx);
-
-  *fLog << " " << setw(7) << "Unreliable Pixels:              "
-        << Form("%s%3i%s%3i","Inner: ",counts[0]," Outer: ",counts[1]) << endl;
-
+  {
+      *fLog << " Successfully calibrated Pixels: Inner: "
+          << Form("%3i",suit[0])   << " Outer: " << Form("%3i",suit[1])   << endl;
+      *fLog << " Uncalibrated Pixels:            Inner: "
+          << Form("%3i",unsuit[0]) << " Outer: " << Form("%3i",unsuit[1]) << endl;
+      *fLog << " Unreliable Pixels:              Inner: "
+          << Form("%3i",unrel[0])  << " Outer: " << Form("%3i",unrel[1])  << endl;
+      *fLog << " High-gain saturated Pixels:     Inner: "
+          << Form("%3i",counts[0]) << " Outer: " << Form("%3i",counts[1]) << endl;
+      *fLog << endl;
+  }
+
+  return kTRUE;
+
+  if (unsuitcnt > fUnsuitablesLimit*fGeom->GetNumPixels())
+  {
+      *fLog << err << "Number of unsuitable pixels: " << 100.*unsuitcnt/fGeom->GetNumPixels()
+          << "% exceeds limit of " << fUnsuitablesLimit*100 << "%" <<  endl;
+      return kFALSE;
+  }
+
+  if (unrelcnt > fUnreliablesLimit*fGeom->GetNumPixels())
+  {
+      *fLog << err << "Relative number of unreliable pixels: " << 100.*unrelcnt/fGeom->GetNumPixels()
+          << "% exceeds limit of " << fUnreliablesLimit*100 << "%" << endl;
+      return kFALSE;
+  }
+  return kTRUE;
 }
 
@@ -2232 +2282 @@
     }
 
-  *fLog << " " << setw(7) << text  
-        << Form("%s%3i%s%3i","Inner: ",countinner," Outer: ",countouter) << endl;
+  *fLog << " " << setw(7) << text << "Inner: " << Form("%3i",countinner)
+      << " Outer: " << Form("%3i", countouter) << endl;
 }
 
@@ -2361 +2411 @@
       rc = kTRUE;
     }
-  if (IsEnvDefined(env, prefix, "UseExtractorRes", print))
-    {
-      SetUseExtractorRes(GetEnvValue(env, prefix, "UseExtractorRes", fUseExtractorRes));
+   if (IsEnvDefined(env, prefix, "UseExtractorRes", print))
+     {
+       SetUseExtractorRes(GetEnvValue(env, prefix, "UseExtractorRes", IsUseExtractorRes()));
       rc = kTRUE;
     }
+  if (IsEnvDefined(env, prefix, "UseUnreliables", print))
+    {
+      SetUseUnreliables(GetEnvValue(env, prefix, "UseUnreliables", IsUseUnreliables()));
+      rc = kTRUE;
+    }
+
+  if (IsEnvDefined(env, prefix, "UseExternalNumPhes", print))
+    {
+      SetUseExternalNumPhes(GetEnvValue(env, prefix, "UseExternalNumPhes", IsUseExternalNumPhes()));
+      rc = kTRUE;
+    }
+
+  if (IsEnvDefined(env, prefix, "UnsuitablesLimit", print))
+    {
+      SetUnsuitablesLimit(GetEnvValue(env, prefix, "UnsuitablesLimit", fUnsuitablesLimit));
+      rc = kTRUE;
+    }
+
+  if (IsEnvDefined(env, prefix, "UnreliablesLimit", print))
+    {
+       SetUnreliablesLimit(GetEnvValue(env, prefix, "UnreliablesLimit", fUnreliablesLimit));
+       rc = kTRUE;
+     }
+ 
 
   return rc;