Changeset 8978 for trunk/MagicSoft/Mars/mbase

Timestamp:

06/18/08 21:39:50 (17 years ago)

Author:

tbretz

Message:

*** empty log message ***

Location:

trunk/MagicSoft/Mars/mbase

Files:

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

trunk/MagicSoft/Mars/mbase/MMath.cc

@@ -1 +1 @@
 /* ======================================================================== *\
-! $Name: not supported by cvs2svn $:$Id: MMath.cc,v 1.38 2008-06-14 15:55:50 tbretz Exp $
+! $Name: not supported by cvs2svn $:$Id: MMath.cc,v 1.39 2008-06-18 20:39:48 tbretz Exp $
 ! --------------------------------------------------------------------------
 !
@@ -811 +811 @@
     return 3;
 }
+
+// --------------------------------------------------------------------------
+//
+// Format a value and its error corresponding to the rules (note
+// this won't work if the error is more then eight orders smaller than
+// the value)
+//
+void MMath::Format(Double_t &v, Double_t &e)
+{
+    // Valid digits
+    Int_t i = TMath::FloorNint(TMath::Log10(v))-TMath::FloorNint(TMath::Log10(e));
+
+    // Check if error starts with 1 or 2. In this case use one
+    // more valid digit
+    TString error = Form("%.0e", e);
+    if (error[0]=='1' || error[0]=='2')
+    {
+        i++;
+        error = Form("%.1e", e);
+    }
+
+    const char *fmt = Form("%%.%de", i);
+
+    v = atof(Form(fmt, v));
+    e = error.Atof();
+}

trunk/MagicSoft/Mars/mbase/MMath.h

@@ -1 +1 @@
 #ifndef MARS_MMath
 #define MARS_MMath
+
+#ifndef _MATH_H
+#include <math.h>  // Needed for darwin
+#endif
 
 #ifndef ROOT_TMath
@@ -69 +73 @@
     TArrayD LeastSqFitPowerLaw(Int_t n, Double_t *x, Double_t *y);
 
-    inline Int_t ModF(Double_t dbl, Double_t &frac) { Double_t rc; frac = modf(dbl, &rc); return TMath::Nint(rc); }
+    inline Int_t ModF(Double_t dbl, Double_t &frac) { Double_t rc; frac = ::modf(dbl, &rc); return TMath::Nint(rc); }
 
-    inline Double_t Sqrt3(Double_t x) { return cbrt(x); }
+    inline Double_t Sqrt3(Double_t x) { return ::cbrt(x); }
 
     inline Double_t Sgn(Double_t d) { return d<0 ? -1 : 1; }
+
+    void Format(Double_t &v, Double_t &e);
 }