Changeset 4224

Timestamp:

05/27/04 18:02:15 (20 years ago)

Author:

gaug

Message:

*** empty log message ***

Location:

trunk/MagicSoft/Mars

Files:

• Changelog (modified) (1 diff)
• mcalib/MCalibrationChargeCalc.cc (modified) (2 diffs)

trunk/MagicSoft/Mars/Changelog

@@ -39 +39 @@
        error remains around 20%). 
 
+   * mcalib/MCalibrationChargeCalc.cc
+     - write the number of photons directy into MCalibrationChargeCam. 
+       These numbers would otherwise have to be assembled by different 
+       parts in separate code. 
 
  2004/05/27: Thomas Bretz

trunk/MagicSoft/Mars/mcalib/MCalibrationChargeCalc.cc

@@ -1239 +1239 @@
 
   const Float_t avphotons   = avpix.GetPheFFactorMethod()       
-                           / qepix.GetQEFFactor(fPulserColor)
+                           / qepix.GetDefaultQE(fPulserColor)
                            / fQECam->GetPlexiglassQE();
 
   const Float_t avphotrelvar = avpix.GetPheFFactorMethodRelVar() 
-                            + qepix.GetQEFFactorRelVar(fPulserColor)
+                            + qepix.GetDefaultQERelVar(fPulserColor)
                             + fQECam->GetPlexiglassQERelVar();
 
@@ -1365 +1365 @@
 
   const UInt_t npixels  = fGeom->GetNumPixels();
-  
+
+  //
+  // Set the results in the MCalibrationChargeCam
+  //
+  if (fBlindPixel)
+    {
+      if (fBlindPixel->IsFluxInsidePlexiglassAvailable())
+        {
+          const Float_t photons = fBlindPixel->GetFluxInsidePlexiglass() * (*fGeom)[0].GetA()
+                               / fQECam->GetPlexiglassQE();
+          fCam->SetNumPhotonsBlindPixelMethod   (photons);
+          const Float_t photrelvar = fBlindPixel->GetFluxInsidePlexiglassRelVar() 
+                                   + fQECam->GetPlexiglassQERelVar();
+          if (photrelvar > 0.)
+            fCam->SetNumPhotonsBlindPixelMethodErr(TMath::Sqrt( photrelvar * photons * photons));  
+        }
+    }
   //
   //  With the knowledge of the overall photon flux, calculate the