Changeset 4334 for trunk/MagicSoft

Timestamp:

06/23/04 00:34:54 (20 years ago)

Author:

gaug

Message:

*** empty log message ***

Location:

trunk/MagicSoft/Mars

Files:

• mcalib/MCalibrationChargeCalc.cc (modified) (4 diffs)
• mcalib/MCalibrationChargeCalc.h (modified) (3 diffs)
• mcalib/MCalibrationChargePix.cc (modified) (12 diffs)
• mjobs/MJCalibration.cc (modified) (2 diffs)
• mjobs/MJCalibration.h (modified) (2 diffs)

trunk/MagicSoft/Mars/mcalib/MCalibrationChargeCalc.cc

@@ -482 +482 @@
   fNumHiGainSamples  =  fCam->GetNumHiGainFADCSlices();
   fNumLoGainSamples  =  fCam->GetNumLoGainFADCSlices();
+
   fSqrtHiGainSamples =  TMath::Sqrt(fNumHiGainSamples);
   fSqrtLoGainSamples =  TMath::Sqrt(fNumLoGainSamples);
@@ -504 +505 @@
         }
       
+      if (IsDebug())
+        pix.SetDebug();
     }
 
@@ -743 +746 @@
   const Float_t num    = TMath::Sqrt((Float_t)fPedestals->GetTotalEntries());
 
+
   //
   // set them in the calibration camera
@@ -815 +819 @@
       return kFALSE;
     }
-  
+
   if (!cal.CalcFFactorMethod())
     {

trunk/MagicSoft/Mars/mcalib/MCalibrationChargeCalc.h

@@ -81 +81 @@
 
   // enums
-  enum  { kHiLoGainCalibration };
+  enum  { kHiLoGainCalibration, kDebug };
 
   // functions
@@ -108 +108 @@
 
   void Clear(const Option_t *o="");
+
+  Bool_t IsDebug() const   {  return TESTBIT(fFlags,kDebug); }
   
   void SetChargeLimit       ( const Float_t f=fgChargeLimit            ) { fChargeLimit       = f;   }
@@ -122 +124 @@
   void SkipHiLoGainCalibration ( const Bool_t b=kTRUE )
       { b ? CLRBIT(fFlags, kHiLoGainCalibration) : SETBIT(fFlags, kHiLoGainCalibration); }
+  void SetDebug              ( const Bool_t b=kTRUE )
+      { b ? CLRBIT(fFlags, kDebug              ) : SETBIT(fFlags, kDebug              ); }
 
   ClassDef(MCalibrationChargeCalc, 1)   // Task calculating Calibration Containers and Quantum Efficiencies

trunk/MagicSoft/Mars/mcalib/MCalibrationChargePix.cc

@@ -695 +695 @@
   const Float_t pedRmsSquareVar = IsHiGainSaturation() ? fLoGainPedRmsSquareVar : 0.25*fPedVar*pedRmsSquare;
 
+  if (IsDebug())
+    {
+      *fLog << dbginf << "ID: " << GetPixId() 
+            << " HiGainSaturation: " << IsHiGainSaturation() 
+            << " Sigma: " << sigma 
+            << " Var.Sigma: " << sigmavar
+            << " PedRmsSquare: " << pedRmsSquare
+            << " pedRmsSquareVar: " << pedRmsSquareVar
+            << endl;
+    }
+  
   const Float_t sigmaSquare    =      sigma     * sigma;
   const Float_t sigmaSquareVar = 4. * sigmavar  * sigmaSquare;
@@ -703 +714 @@
   fRSigmaSquare = sigmaSquare - pedRmsSquare;
 
+  if (IsDebug())
+    {
+      *fLog << dbginf << "ID: " << GetPixId() 
+            << " Red.Sigma Square: " << fRSigmaSquare
+            << endl;
+    }
+  
   if (fRSigmaSquare <= 0.)
     {
@@ -711 +729 @@
     }
 
+  
   fRSigmaSquareVar = 4. * (sigmaSquareVar + pedRmsSquareVar);
+
+  if (IsDebug())
+    {
+      *fLog << dbginf << "ID: " << GetPixId() 
+            << " Var.Red.Sigma Square: " << fRSigmaSquareVar
+            << endl;
+    }
 
   return kTRUE;
@@ -746 +772 @@
 {
 
+
   if (fRSigmaSquare < 0.)
     return kFALSE;
@@ -765 +792 @@
   fPheFFactorMethod = ffactorsquare * meanSquare / fRSigmaSquare;
 
+  if (IsDebug())
+    {
+      *fLog << dbginf << "ID: " << GetPixId() 
+            << " F-Factor Square: " << ffactorsquare
+            << " Mean Square: " << meanSquare
+            << " Red.Sigma Square: " << fRSigmaSquare
+            << " Photo-electrons: " << fPheFFactorMethod
+            << endl;
+    }
+
   if (fPheFFactorMethod < fPheFFactorMethodLimit)
     return kFALSE;
@@ -774 +811 @@
   fPheFFactorMethodVar =  pheRelVar * fPheFFactorMethod * fPheFFactorMethod;
 
+  if (IsDebug())
+    {
+      *fLog << dbginf << "ID: " << GetPixId() 
+            << " Rel.Var.F-Factor Square: " << ffactorsquareRelVar
+            << " Rel.Var. Mean Square: " << meanSquareRelVar
+            << " Rel.Var. Red.Sigma Square: " << rsigmaSquareRelVar
+            << " Rel.Var. Photo-electrons: " << pheRelVar
+            << endl;
+    }
+
   if (fPheFFactorMethodVar < 0. )
     return kFALSE;
-
-  fMeanConvFADC2Phe    =  fPheFFactorMethod / GetConvertedMean();
-  
+  
+  const Float_t convmean = GetConvertedMean();
+  
+
+  if (convmean > 0.)
+    fMeanConvFADC2Phe    =  fPheFFactorMethod / GetConvertedMean();
+  else
+    fMeanConvFADC2Phe    =  -1.;
+
+  if (IsDebug())
+    {
+      *fLog << dbginf << "ID: " << GetPixId() 
+            << " Converted Mean: " << convmean 
+            << " Conversion FADC2Phe: " << fMeanConvFADC2Phe
+            << endl;
+    }
+
   if (fMeanConvFADC2Phe < 0. )
     return kFALSE;
@@ -790 +851 @@
   const Float_t limit      = IsHiGainSaturation() ? fConvFFactorRelVarLimit * 4. : fConvFFactorRelVarLimit;
 
+  if (IsDebug())
+    {
+      *fLog << dbginf << "ID: " << GetPixId() 
+            << " Rel.Var.Red.Sigma: " << rsigmaSquareRelVar
+            << " Rel.Var.Mean: " << GetMeanRelVar()
+            << " Rel.Var.F-Factor: " << ffactorsquareRelVar
+            << " Rel.Var.Conversion FADC2Phe: " << convrelvar
+            << endl;
+    }
+
   if (convrelvar > limit || convrelvar < 0.)
     {
@@ -816 +887 @@
 {
 
+
+  if (IsDebug())
+    {
+      *fLog << dbginf << "ID: " << GetPixId() 
+            << " Number photons: " << nphotons
+            << " Rel.Var.Number photons: " << nphotonsrelvar
+            << " Red.Sigma Square: " << fRSigmaSquare
+            << " Mean: " << GetMean()
+            << endl;
+    }
+
+
   if (nphotons <= 0.)
     {
@@ -830 +913 @@
   fMeanFFactorFADC2Phot =  TMath::Sqrt(fRSigmaSquare * nphotons) / GetMean() ;
   
+  if (IsDebug())
+    {
+      *fLog << dbginf << "ID: " << GetPixId() 
+            << " F-Factor FADC2Phot: " << fMeanFFactorFADC2Phot
+            << endl;
+    }
+
   if (fMeanFFactorFADC2Phot < 0.)
     {
@@ -841 +931 @@
   
   fMeanFFactorFADC2PhotVar    = ffactorrelvar * fMeanFFactorFADC2Phot * fMeanFFactorFADC2Phot;
+
+  if (IsDebug())
+    {
+      *fLog << dbginf << "ID: " << GetPixId() 
+            << " Rel.Var.Red.Sigma: " << 0.25 * fRSigmaSquareVar / ( fRSigmaSquare * fRSigmaSquare) 
+            << " Rel.Var.Mean: " << GetMeanRelVar()
+            << " Rel.Var.photons: " << 0.25 * nphotonsrelvar
+            << " Rel.Var.F-Factor FADC2Phot: " << ffactorrelvar
+            << endl;
+    }
 
   return kTRUE;
@@ -896 +996 @@
   const Float_t elecRmsSquareVar = 4.*elecPedRmsVar * elecRmsSquare;
   
+  if (IsDebug())
+    {
+      *fLog << dbginf << "ID: " << GetPixId() 
+            << " Ped.Rms Square:  " << pedRmsSquare 
+            << " Elec.Rms Square: " << elecRmsSquare 
+            << " Ped.Rms.Square Var.: " << pedRmsSquareVar 
+            << " Elec.Rms Square Var.: " << elecRmsSquareVar 
+            << endl;
+    }
+  
+
   Float_t higainNsbSquare        =  pedRmsSquare    - elecRmsSquare;
-  Float_t higainNsbSquareRelVar  = (pedRmsSquareVar + elecRmsSquareVar)
-                                 / (higainNsbSquare * higainNsbSquare) ;
-  
-  if (higainNsbSquare < 0.)
+  Float_t higainNsbSquareRelVar  = (pedRmsSquareVar + elecRmsSquareVar);
+
+
+  if (higainNsbSquare < 0.001)
     higainNsbSquare = 0.;
+  else
+    higainNsbSquareRelVar  /= (higainNsbSquare * higainNsbSquare) ;
+  
+  if (IsDebug())
+    {
+      *fLog << dbginf << "ID: " << GetPixId() 
+            << " HiGain NSB Square:  " << higainNsbSquare
+            << " Rel.Var.HiGain NSB Square:  " << higainNsbSquareRelVar
+            << endl;
+    }
   
   //
@@ -917 +1038 @@
   fLoGainPedRmsSquareVar = logainNsbSquareVar + elecRmsSquareVar;
 
+  if (IsDebug())
+    {
+      *fLog << dbginf << "ID: " << GetPixId() 
+            << " LoGain Ped Rms Square:  " << fLoGainPedRmsSquare
+            << " Var.Ped Rms Square:  " << fLoGainPedRmsSquareVar
+            << endl;
+    }
+  
+  
 }
  

trunk/MagicSoft/Mars/mjobs/MJCalibration.cc

@@ -155 +155 @@
     : fRuns(0), fExtractor(NULL), fTimeExtractor(NULL), 
       fColor(MCalibrationCam::kNONE), fDisplayType(kNormalDisplay),
-      fRelTimes(kFALSE), fDataCheck(kFALSE)
+      fRelTimes(kFALSE), fDataCheck(kFALSE), fDebug(kFALSE)
 {
 
@@ -1004 +1004 @@
   calcalc.SetOutputPath(fOutputPath);
   calcalc.SetOutputFile(Form("%s-ChargeCalibStat.txt",(const char*)fRuns->GetRunsAsFileName()));
+  if (fDebug)
+    calcalc.SetDebug();
 
   MCalibrationRelTimeCalc  timecalc;

trunk/MagicSoft/Mars/mjobs/MJCalibration.h

@@ -57 +57 @@
   Bool_t fRelTimes;                                        // Flag if relative times have to be calibrated
   Bool_t fDataCheck;                                       // Flag if the data check is run on raw data
-  
+  Bool_t fDebug;
+
   void   DisplayResult(MParList &plist);
   Bool_t WriteResult();
@@ -97 +98 @@
   void SetDataCheck         (const Bool_t b=kTRUE) { fDataCheck        = b; SetDataCheckDisplay(); }
 
+  // Debug
+  void SetDebug             (const Bool_t b=kTRUE) { fDebug           = b; }
+  
   // Devices
   void SetUseBlindPixel( const Bool_t b=kTRUE );