Changeset 7386 for trunk/MagicSoft/Mars/mhflux/MHAlpha.cc

Timestamp:

10/11/05 14:34:04 (19 years ago)

Author:

tbretz

Message:

*** empty log message ***

File:

• trunk/MagicSoft/Mars/mhflux/MHAlpha.cc (modified) (4 diffs)

trunk/MagicSoft/Mars/mhflux/MHAlpha.cc

@@ -163 +163 @@
         {
             fHEnergy.SetBinContent(i, fit.GetEventsExcess());
-            if (fit.GetSignificance()>1)
-                fHEnergy.SetBinError(i, fit.GetEventsExcess()/fit.GetSignificance());
+            if (fit.GetSignificanceExc()>1)
+                fHEnergy.SetBinError(i, fit.GetEventsExcess()/fit.GetSignificanceExc());
             else
                 fHEnergy.SetBinError(i, fit.GetEventsExcess());
 
-            if (fit.GetSignificance()>1)
+            if (fit.GetSignificanceExc()>1)
             {
-                mean += fit.GetSignificance();
+                mean += fit.GetSignificanceExc()*fit.GetSignificanceExc();
                 num++;
             }
         }
     }
-    return mean/num;
+    return TMath::Sqrt(mean)/num;
 }
 
@@ -192 +192 @@
         {
             fHTheta.SetBinContent(i, fit.GetEventsExcess());
-            if (fit.GetSignificance()>1)
-                fHTheta.SetBinError(i, fit.GetEventsExcess()/fit.GetSignificance());
+            if (fit.GetSignificanceExc()>1)
+                fHTheta.SetBinError(i, fit.GetEventsExcess()/fit.GetSignificanceExc());
             else
                 fHTheta.SetBinError(i, fit.GetEventsExcess());
@@ -430 +430 @@
     //
     fHTime.SetBinContent(n+1, fit.GetEventsExcess());
-    if (fit.GetSignificance()>1)
-        fHTime.SetBinError(n+1, fit.GetEventsExcess()/fit.GetSignificance());
+    if (fit.GetSignificanceExc()>1)
+        fHTime.SetBinError(n+1, fit.GetEventsExcess()/fit.GetSignificanceExc());
     else
         fHTime.SetBinError(n+1, fit.GetEventsExcess());
@@ -810 +810 @@
         if (hof ? fit.Fit(*hon, *hof, alpha) : fit.Fit(*hon))
         {
-            *fLog << dbg << "Bin " << i << ": sigma=" << fit.GetSignificance() << " omega=" << fit.GetGausSigma() << " events=" << fit.GetEventsExcess() << " scale=" << fit.GetScaleFactor() << endl;
+            *fLog << dbg << "Bin " << i << ": sigmaexc=" << fit.GetSignificanceExc() << " omega=" << fit.GetGausSigma() << " events=" << fit.GetEventsExcess() << " scale=" << fit.GetScaleFactor() << endl;
             fit.PaintResult();
         }