Changeset 4716 for trunk/MagicSoft/Mars/mbase

Timestamp:

08/23/04 17:49:53 (20 years ago)

Author:

tbretz

Message:

*** empty log message ***

Location:

trunk/MagicSoft/Mars/mbase

Files:

• MMath.cc (modified) (3 diffs)
• MMath.h (modified) (1 diff)

trunk/MagicSoft/Mars/mbase/MMath.cc

@@ -47 +47 @@
     const Double_t f = s+k*k*b;
 
-    return f==0 ? 0 : (s-b)/TMath::Sqrt(f);
+    return f==0 ? 0 : (s-b)/Sqrt(f);
 }
 
@@ -86 +86 @@
         return -1;
 
-    const Double_t l = s*TMath::Log(s/sum*(alpha+1)/alpha);
-    const Double_t m = b*TMath::Log(b/sum*(alpha+1)      );
+    const Double_t l = s*Log(s/sum*(alpha+1)/alpha);
+    const Double_t m = b*Log(b/sum*(alpha+1)      );
 
-    return l+m<0 ? -1 : TMath::Sqrt((l+m)*2);
+    return l+m<0 ? -1 : Sqrt((l+m)*2);
 }
 
@@ -104 +104 @@
         return 0;
 
-    return TMath::Sign(sig, s-b);
+    return Sign(sig, s-b);
 }
+
+// --------------------------------------------------------------------------
+//
+// Returns: 2/(sigma*sqrt(2))*integral[0,x](exp(-(x-mu)^2/(2*sigma^2)))
+//
+Double_t MMath::GaussProb(Double_t x, Double_t sigma, Double_t mean)
+{
+    static const Double_t sqrt2 = Sqrt(2.);
+    return Erf((x-mean)/(sigma*sqrt2));
+}
+

trunk/MagicSoft/Mars/mbase/MMath.h

@@ -9 +9 @@
 {
 public:
+    static Double_t GaussProb(Double_t x, Double_t sigma, Double_t mean=0);
+
     static Double_t Significance(Double_t s, Double_t b);
     static Double_t SignificanceSym(Double_t s, Double_t b);