Changeset 7386 for trunk/MagicSoft/Mars/mhflux

Timestamp:

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

Author:

tbretz

Message:

*** empty log message ***

Location:

trunk/MagicSoft/Mars/mhflux

Files:

• MAlphaFitter.cc (modified) (1 diff)
• MAlphaFitter.h (modified) (1 diff)
• MHAlpha.cc (modified) (4 diffs)

trunk/MagicSoft/Mars/mhflux/MAlphaFitter.cc

@@ -34 +34 @@
 // container in the parlisrt
 //
-// PRELIMINARY!
+// Version 2:
+// ----------
+//  +Double_t fSignificanceExc;  // significance of excess
+//
 //
 //////////////////////////////////////////////////////////////////////////////

trunk/MagicSoft/Mars/mhflux/MAlphaFitter.h

@@ -220 +220 @@
     Int_t ReadEnv(const TEnv &env, TString prefix, Bool_t print=kFALSE);
 
-    ClassDef(MAlphaFitter, 1)
+    ClassDef(MAlphaFitter, 2)
 };
 

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();
         }