Changeset 3583 for trunk/MagicSoft/Mars/mhist

Timestamp:

03/23/04 09:11:11 (21 years ago)

Author:

tbretz

Message:

*** empty log message ***

File:

• trunk/MagicSoft/Mars/mhist/MHFalseSource.cc (modified) (5 diffs)

trunk/MagicSoft/Mars/mhist/MHFalseSource.cc

@@ -119 +119 @@
 #include "MAstroSky2Local.h"
 #include "MStatusDisplay.h"
+#include "MMath.h"
 
 #include "MBinning.h"
@@ -248 +249 @@
 Double_t MHFalseSource::Significance(Double_t s, Double_t b)
 {
-    const Double_t k = b==0 ? 0 : s/b;
-    const Double_t f = s+k*k*b;
-
-    return f==0 ? 0 : (s-b)/TMath::Sqrt(f);
+    return MMath::Significance(s, b);
 }
 
@@ -265 +263 @@
 Double_t MHFalseSource::SignificanceLiMa(Double_t s, Double_t b, Double_t alpha)
 {
-    const Double_t sum = s+b;
-
-    if (sum<=0 || alpha<=0)
-        return 0;
-
-    const Double_t l = s*TMath::Log(s/sum*(alpha+1)/alpha);
-    const Double_t m = b*TMath::Log(b/sum*(alpha+1)      );
-
-    return TMath::Sqrt((l+m)*2);
+    return MMath::SignificanceLiMa(s, b);
 }
 
@@ -712 +702 @@
     const Int_t nx = fHist.GetXaxis()->GetNbins();
     const Int_t ny = fHist.GetYaxis()->GetNbins();
+    const Int_t nr = nx*nx+ny*ny;
 
     Double_t maxalpha0=0;
@@ -825 +816 @@
                 h6.Fill(sig);
 
-            if (sig>maxs)
+            if (sig>maxs && ix*ix+iy*iy<nr*nr/9)
             {
                 maxs = sig;