Changeset 3683

Timestamp:

04/08/04 19:56:26 (20 years ago)

Author:

gaug

Message:

*** empty log message ***

Location:

trunk/MagicSoft/Mars/mcalib

Files:

• MCalibrationChargeCalc.cc (modified) (3 diffs)
• MCalibrationChargeCam.cc (modified) (1 diff)
• MCalibrationChargePix.cc (modified) (4 diffs)
• MCalibrationChargePix.h (modified) (1 diff)

trunk/MagicSoft/Mars/mcalib/MCalibrationChargeCalc.cc

@@ -1009 +1009 @@
   MCalibrationQEPix     &qepix = (MCalibrationQEPix&)  fQECam->GetAverageArea(0);  
 
-  const Float_t avphotonflux = avpix.GetPheFFactorMethod()       
-                             / qepix.GetQEFFactor(fPulserColor)
-                             / fQECam->GetPlexiglassQE();
-
-  const Float_t avfluxrelvar = avpix.GetPheFFactorMethodRelVar() 
-                             + qepix.GetQEFFactorRelVar(fPulserColor)
-                             + fQECam->GetPlexiglassQERelVar();
+  const Float_t avphotons   = avpix.GetPheFFactorMethod()       
+                           / qepix.GetQEFFactor(fPulserColor)
+                           / fQECam->GetPlexiglassQE();
+
+  const Float_t avphotrelvar = avpix.GetPheFFactorMethodRelVar() 
+                            + qepix.GetQEFFactorRelVar(fPulserColor)
+                            + fQECam->GetPlexiglassQERelVar();
 
   const UInt_t npixels  = fGeom->GetNumPixels();
@@ -1031 +1031 @@
         }
       
-      const Float_t photons = avphotonflux / fGeom->GetPixRatio(i);
+      const Float_t photons = avphotons / fGeom->GetPixRatio(i);
       const Float_t qe      = pix.GetPheFFactorMethod() / photons ;
 
-      if (!pix.CalcMeanFFactor( photons , avfluxrelvar ))
+      if (!pix.CalcMeanFFactor( photons , avphotrelvar ))
         {
           pix.SetFFactorMethodValid(kFALSE);
@@ -1041 +1041 @@
         }
 
-      const Float_t qerelvar = avfluxrelvar +  pix.GetPheFFactorMethodRelVar();
+      const Float_t qerelvar = avphotrelvar +  pix.GetPheFFactorMethodRelVar();
 
       qepix.SetQEFFactor    ( qe            , fPulserColor );

trunk/MagicSoft/Mars/mcalib/MCalibrationChargeCam.cc

@@ -497 +497 @@
       if (pix.IsExcluded() || !pix.IsFFactorMethodValid())
         return kFALSE;
-      val = pix.GetMeanFFactorFADC2Phot();
+      val = pix.GetMeanFFactorFADC2PhotErr();
       break;
     case 15:

trunk/MagicSoft/Mars/mcalib/MCalibrationChargePix.cc

@@ -659 +659 @@
 //
 // Calculate the total F-Factor with the formula:
-//   fMeanFFactorFADC2Phot = fRSigmaSquare / GetMeanSquare() * sqrt(photflux)
+//   fMeanFFactorFADC2Phot = Sqrt ( fRSigmaSquare / GetMeanSquare() ) * nphotons
 //
 // Calculate the error of the total F-Factor
 //
-Bool_t MCalibrationChargePix::CalcMeanFFactor( const Float_t photflux, const Float_t photfluxrelvar )
-{
-
-  if (photflux <= 0.)
+Bool_t MCalibrationChargePix::CalcMeanFFactor( const Float_t nphotons, const Float_t nphotonsrelvar )
+{
+
+  if (nphotons <= 0.)
     {
       *fLog << warn << GetDescriptor() << ": Assumed photon flux is smaller or equal 0." << endl;
@@ -672 +672 @@
     }
 
-  if (photfluxrelvar < 0.)
+  if (nphotonsrelvar < 0.)
     {
       *fLog << warn << GetDescriptor() << ": Assumed photon flux variance is smaller than 0." << endl;
@@ -678 +678 @@
     }
 
-  fMeanFFactorFADC2Phot =  fRSigmaSquare / GetMeanSquare() * TMath::Sqrt(photflux);
+  fMeanFFactorFADC2Phot =  TMath::Sqrt(fRSigmaSquare) / GetMean() * nphotons;
   
   if (fMeanFFactorFADC2Phot < 0.)
@@ -686 +686 @@
     }
   
-  const Float_t ffactorrelvar = fRSigmaSquareVar / ( fRSigmaSquare * fRSigmaSquare) 
-                              + 4.   * GetMeanRelVar() 
-                              + 0.25 * photfluxrelvar;
+  const Float_t ffactorrelvar = 0.25 * fRSigmaSquareVar / ( fRSigmaSquare * fRSigmaSquare) 
+                              + GetMeanRelVar() 
+                              + nphotonsrelvar;
   
   fMeanFFactorFADC2PhotVar    = ffactorrelvar * fMeanFFactorFADC2Phot * fMeanFFactorFADC2Phot;

trunk/MagicSoft/Mars/mcalib/MCalibrationChargePix.h

@@ -102 +102 @@
   Bool_t CalcReducedSigma  ();
   Bool_t CalcFFactorMethod ();
-  Bool_t CalcMeanFFactor   ( const Float_t photflux, const Float_t phtfluxrelvar );
+  Bool_t CalcMeanFFactor   ( const Float_t nphotons, const Float_t nphotonsrelvar );
   
   ClassDef(MCalibrationChargePix, 1)    // Container Charge Calibration Results Pixel