Changeset 3678 for trunk/MagicSoft/Mars/mcalib/MCalibrationChargeBlindPix.cc

Timestamp:

04/08/04 17:58:54 (21 years ago)

Author:

gaug

Message:

*** empty log message ***

File:

• trunk/MagicSoft/Mars/mcalib/MCalibrationChargeBlindPix.cc (modified) (4 diffs)

trunk/MagicSoft/Mars/mcalib/MCalibrationChargeBlindPix.cc

@@ -53 +53 @@
 /////////////////////////////////////////////////////////////////////////////
 #include "MCalibrationChargeBlindPix.h"
+#include "MCalibrationCam.h"
 
 #include <TH1.h>
@@ -184 +185 @@
 // --------------------------------------------------------------------------
 //
+// Return -1 if fFluxInsidePlexiglassVar is smaller than 0.
+// Return -1 if fFluxInsidePlexiglass    is 0.
+// Return fFluxInsidePlexiglassVar / fFluxInsidePlexiglass^2
+// 
+Float_t MCalibrationChargeBlindPix::GetFluxInsidePlexiglassRelVar() const
+{
+  if (fFluxInsidePlexiglassVar < 0.)
+    return -1.;
+
+  if (fFluxInsidePlexiglass == 0.)
+    return -1.;
+  
+  return fFluxInsidePlexiglassVar / (fFluxInsidePlexiglass * fFluxInsidePlexiglass) ;
+}
+
+// --------------------------------------------------------------------------
+//
 // Return -1 if fLambdaVar is smaller than 0.
 // Return square root of fLambdaVar
@@ -279 +297 @@
   return gkBlindPixelQECT1Err * gkBlindPixelQECT1Err / gkBlindPixelQECT1 / gkBlindPixelQECT1 ;
 }
-
 
 // --------------------------------------------------------------------------
@@ -355 +372 @@
   switch (fColor)
     {
-    case kGREEN:
+    case MCalibrationCam::kGREEN:
       fFluxInsidePlexiglass    = fLambda * gkBlindPixelQEGreen * TMath::Power(10,gkBlindPixelAttGreen);   
       // attenuation has negligible error
       fFluxInsidePlexiglassVar = GetLambdaRelVar() + GetBlindPixelQEGreenRelVar();
       break;
-    case kBLUE:
+    case MCalibrationCam::kBLUE:
       fFluxInsidePlexiglass    = fLambda * gkBlindPixelQEBlue * TMath::Power(10,gkBlindPixelAttBlue);   
       // attenuation has negligible error
       fFluxInsidePlexiglassVar = GetLambdaRelVar() + GetBlindPixelQEBlueRelVar();
       break;
-    case kUV:
+    case MCalibrationCam::kUV:
       fFluxInsidePlexiglass    = fLambda * gkBlindPixelQEUV * TMath::Power(10,gkBlindPixelAttUV);   
       // attenuation has negligible error
       fFluxInsidePlexiglassVar = GetLambdaRelVar() + GetBlindPixelQEUVRelVar();
       break;
-    case kCT1:
+    case MCalibrationCam::kCT1:
     default:
       fFluxInsidePlexiglass    = fLambda * gkBlindPixelQECT1 * TMath::Power(10,gkBlindPixelAttCT1);