Changeset 3597 for trunk/MagicSoft/Mars/mbase

Timestamp:

03/24/04 15:08:19 (21 years ago)

Author:

tbretz

Message:

*** empty log message ***

File:

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

trunk/MagicSoft/Mars/mbase/MMath.cc

@@ -59 +59 @@
 //  alpha = t_on/t_off;  // t: observation time
 //
-//  Returns -1 if calculation failed...
+//  Returns -1 if sum<0 or alpha<0 or the argument of sqrt<0
+//  Returns  0 if s+b==0
 //
 Double_t MMath::SignificanceLiMa(Double_t s, Double_t b, Double_t alpha)
@@ -65 +66 @@
     const Double_t sum = s+b;
 
-    if (sum<=0 || alpha<=0)
+    if (sum==0)
+        return 0;
+
+    if (sum<0 || alpha<=0)
         return -1;
 
@@ -71 +75 @@
     const Double_t m = b*TMath::Log(b/sum*(alpha+1)      );
 
-    return TMath::Sqrt((l+m)*2);
+    return l+m<0 ? -1 : TMath::Sqrt((l+m)*2);
 }