Changeset 3671 for trunk/MagicSoft/Mars/mcalib/MCalibrationPix.cc

Timestamp:

04/06/04 16:54:53 (21 years ago)

Author:

gaug

Message:

*** empty log message ***

File:

• trunk/MagicSoft/Mars/mcalib/MCalibrationPix.cc (modified) (5 diffs)

trunk/MagicSoft/Mars/mcalib/MCalibrationPix.cc

@@ -127 +127 @@
 // --------------------------------------------------------------------------
 //
+// Get the Relative Variance of either High Gain or Low Gain Mean 
+// depending on IsHighGainSaturation()
+//
+// If variance is smaller than 0. return -1.
+//
+Float_t MCalibrationPix::GetMeanRelVar()  const
+{
+
+  if (IsHiGainSaturation())
+    if (fLoGainMeanVar < 0.)
+      return -1.;
+    else
+      return fLoGainMeanVar / (fLoGainMean * fLoGainMean);
+  else
+    if (fHiGainMeanVar < 0.)
+      return -1.;
+    else
+      return fHiGainMeanVar / (fHiGainMean * fHiGainMean);
+}
+
+// --------------------------------------------------------------------------
+//
+// Get the Relative Variance of either High Gain or Low Gain Sigma 
+// depending on IsHighGainSaturation()
+//
+// If variance is smaller than 0. return -1.
+//
+Float_t MCalibrationPix::GetSigmaRelVar()  const
+{
+
+  if (IsHiGainSaturation())
+    if (fLoGainSigmaVar < 0.)
+      return -1.;
+    else
+      return fLoGainSigmaVar / (fLoGainSigma * fLoGainSigma);
+  else
+    if (fHiGainSigmaVar < 0.)
+      return -1.;
+    else
+      return fHiGainSigmaVar / (fHiGainSigma * fHiGainSigma);
+}
+
+// --------------------------------------------------------------------------
+//
 // Get the High Gain Mean Error: Takes square root of fHiGainMeanVar
 //
 Float_t MCalibrationPix::GetHiGainMeanErr()  const
 {
-  return TMath::Sqrt(fHiGainMeanVar);
-}
+  if (fLoGainMeanVar < 0.)
+    return -1.;
+
+  return TMath::Sqrt(fLoGainMeanVar);
+}
+
 
 // --------------------------------------------------------------------------
@@ -140 +188 @@
 Float_t MCalibrationPix::GetHiGainSigmaErr()  const
 {
+  if (fHiGainSigmaVar < 0.)
+    return -1.;
+
   return TMath::Sqrt(fHiGainSigmaVar);
 }
@@ -145 +196 @@
 // --------------------------------------------------------------------------
 //
-// Get the Low Gain Mean Error: Takes square root of fHiGainMeanVar
+// Get the Low Gain Mean Error: Takes square root of fLoGainMeanVar
 //
 Float_t MCalibrationPix::GetLoGainMeanErr()  const
 {
+  if (fLoGainMeanVar < 0.)
+    return -1.;
+
   return TMath::Sqrt(fLoGainMeanVar);
 }
 
-
 // --------------------------------------------------------------------------
 //
@@ -159 +212 @@
 Float_t MCalibrationPix::GetLoGainSigmaErr()  const
 {
+  if (fLoGainSigmaVar < 0.)
+    return -1.;
+
   return TMath::Sqrt(fLoGainSigmaVar);
 }
@@ -168 +224 @@
 Bool_t MCalibrationPix::IsHiGainSaturation()    const
 { 
-   return TESTBIT(fFlags,kHiGainSaturation);  
+  return TESTBIT(fFlags,kHiGainSaturation);  
 }