Changeset 7672 for trunk/MagicSoft/Mars/mbase

Timestamp:

04/30/06 21:01:33 (19 years ago)

Author:

tbretz

Message:

*** empty log message ***

File:

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

trunk/MagicSoft/Mars/mbase/MMath.cc

@@ -78 +78 @@
 //  The significance has the same (positive!) value for s>b and b>s.
 //
-//  Returns -1 if sum<0 or alpha<0 or the argument of sqrt<0
-//  Returns  0 if s+b==0, s==0 or b==0
+//  Returns -1 if s<0 or b<0 or alpha<0 or the argument of sqrt<0
 //
 // Here is some eMail written by Daniel Mazin about the meaning of the arguments:
@@ -107 +106 @@
     const Double_t sum = s+b;
 
-    if (s==0 || b==0 || sum==0)
-        return 0;
-
-    if (sum<0 || alpha<=0)
+    if (s<0 || b<0 || alpha<=0)
         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==0 ? 0 : s*TMath::Log(s/sum*(alpha+1)/alpha);
+    const Double_t m = b==0 ? 0 : b*TMath::Log(b/sum*(alpha+1)      );
 
     return l+m<0 ? -1 : TMath::Sqrt((l+m)*2);