Changeset 11426 for trunk/Mars

Timestamp:

07/15/11 15:46:14 (13 years ago)

Author:

tbretz

Message:

Added GaussProb2D

Location:

trunk/Mars/mbase

Files:

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

trunk/Mars/mbase/MMath.cc

@@ -223 +223 @@
 }
 
+// --------------------------------------------------------------------------
+//
+// Returns: 1 - exp((x/sigma)^2 / 2)
+//
+Double_t MMath::GaussProb2D(Double_t x, Double_t sigma)
+{
+    const Double_t xs = x/sigma;
+    return 1 - TMath::Exp(-xs*xs/2);
+}
+
 // ------------------------------------------------------------------------
 //
@@ -720 +730 @@
 //     Q := (3*b - a^2)/9
 //     R := (9*b*a - 27*c - 2*a^3)/54
-//     D := Q^3 + R^2                          % polynomial discriminant
+//     D := R^2 -Q^3                           % polynomial discriminant
 // 
 //     if (D >= 0) then                        % complex or duplicate roots

trunk/Mars/mbase/MMath.h

@@ -17 +17 @@
 {
     Double_t GaussProb(Double_t x, Double_t sigma=1, Double_t mean=0);
+    Double_t GaussProb2D(Double_t x, Double_t sigma=1);
 
     template <class Size, class Element> void ReSortImp(Size n, Element *a, Bool_t down=kFALSE);