Changeset 5031 for trunk/MagicSoft/Mars/mcalib/MCalibrationBlindCam.cc

Timestamp:

09/15/04 15:31:22 (20 years ago)

Author:

gaug

Message:

*** empty log message ***

File:

• trunk/MagicSoft/Mars/mcalib/MCalibrationBlindCam.cc (modified) (1 diff)

trunk/MagicSoft/Mars/mcalib/MCalibrationBlindCam.cc

@@ -106 +106 @@
 }
 
+// --------------------------------------------------------------------------
+//
+// Return true, if any of the blind pixels have an available photon flux
+//
+Bool_t MCalibrationBlindCam::IsFluxInsidePlexiglassAvailable() const
+{
+  for (Int_t i=0; i<GetSize(); i++)
+    {
+      MCalibrationBlindPix &bp = (MCalibrationBlindPix&)(*this)[i];      
+      if (bp.IsFluxInsidePlexiglassAvailable())
+        return kTRUE;
+    }
+  
+  return kFALSE;
+}
 
+// --------------------------------------------------------------------------
+//
+// Returns weighted average of the flux calculated by each blind pixel
+//
+Float_t MCalibrationBlindCam::GetFluxInsidePlexiglass() const
+{
+
+  Float_t flux    = 0.;
+  Float_t fluxvar = 0.;  
+
+  for (Int_t i=0; i<GetSize(); i++)
+    {
+      MCalibrationBlindPix &bp = (MCalibrationBlindPix&)(*this)[i];
+      if (bp.IsFluxInsidePlexiglassAvailable())
+        {
+          const Float_t weight = 1./ 
+                                bp.GetFluxInsidePlexiglassErr() 
+                                / bp.GetFluxInsidePlexiglassErr();
+          flux    += weight * bp.GetFluxInsidePlexiglass();
+          fluxvar += weight;
+        }
+    }
+  return fluxvar > 0.0001 ? flux / fluxvar : -1.;
+}
+
+// --------------------------------------------------------------------------
+//
+// Returns weighted variance of the flux calculated by each blind pixel
+//
+Float_t MCalibrationBlindCam::GetFluxInsidePlexiglassVar() const
+{
+
+  Float_t fluxvar = 0.;  
+
+  for (Int_t i=0; i<GetSize(); i++)
+    {
+      MCalibrationBlindPix &bp = (MCalibrationBlindPix&)(*this)[i];
+      if (bp.IsFluxInsidePlexiglassAvailable())
+        {
+          const Float_t weight = 1./ 
+                                bp.GetFluxInsidePlexiglassErr() 
+                                / bp.GetFluxInsidePlexiglassErr();
+          fluxvar += weight;
+        }
+    }
+  return fluxvar > 0.0001 ? 1. / fluxvar : -1.;
+}
+
+// --------------------------------------------------------------------------
+//
+// Returns weighted rel. variance of the flux calculated by each blind pixel
+//
+Float_t MCalibrationBlindCam::GetFluxInsidePlexiglassRelVar() const
+{
+
+  Float_t flux    = 0.;
+  Float_t fluxvar = 0.;  
+
+  for (Int_t i=0; i<GetSize(); i++)
+    {
+      MCalibrationBlindPix &bp = (MCalibrationBlindPix&)(*this)[i];
+      if (bp.IsFluxInsidePlexiglassAvailable())
+        {
+          const Float_t weight = 1./ 
+                                bp.GetFluxInsidePlexiglassErr() 
+                                / bp.GetFluxInsidePlexiglassErr();
+          flux    += weight * bp.GetFluxInsidePlexiglass();
+          fluxvar += weight;
+        }
+    }
+
+  return fluxvar > 0.0001 ? flux * fluxvar : -1.;
+}
+
+