Changeset 11426 for trunk/Mars/mbase/MMath.cc

Timestamp:

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

Author:

tbretz

Message:

Added GaussProb2D

File:

• trunk/Mars/mbase/MMath.cc (modified) (2 diffs)

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