Changeset 7867 for trunk/MagicSoft/Mars/mbase

Timestamp:

08/08/06 17:51:35 (18 years ago)

Author:

tbretz

Message:

*** empty log message ***

Location:

trunk/MagicSoft/Mars/mbase

Files:

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

trunk/MagicSoft/Mars/mbase/MMath.cc

@@ -161 +161 @@
 }
 
+// ------------------------------------------------------------------------
+//
+// Return the median value of the distribution of abs(a[i]-Median)
+//
+template <class Size, class Element>
+Double_t MMath::MedianDevImp(Size n, const Element *a)
+{
+    // Get median of distribution
+    const Double_t med = TMath::Median(n, a);
+
+    // Allocate space for distribution
+    Double_t *arr = new Double_t[n];
+
+    // Create the abs(a[i]-med) distribution
+    for (int i=0; i<n; i++)
+        arr[i] = TMath::Abs(a[i]-med);
+
+    // FIXME: GausProb() is a workaround. It should be taken into account in Median!
+    const Double_t rc = TMath::Median(n, arr);
+
+    // delete space
+    delete arr;
+
+    // return result
+    return rc;
+}
+
+// ------------------------------------------------------------------------
+//
+// Return the median value of the distribution of abs(a[i]-Median)
+//
+Double_t MMath::MedianDev(Long64_t n, const Short_t *a)
+{
+    return MedianDevImp(n, a);
+}
+
+// ------------------------------------------------------------------------
+//
+// Return the median value of the distribution of abs(a[i]-Median)
+//
+Double_t MMath::MedianDev(Long64_t n, const Int_t *a)
+{
+    return MedianDevImp(n, a);
+}
+
+// ------------------------------------------------------------------------
+//
+// Return the median value of the distribution of abs(a[i]-Median)
+//
+Double_t MMath::MedianDev(Long64_t n, const Float_t *a)
+{
+    return MedianDevImp(n, a);
+}
+
+// ------------------------------------------------------------------------
+//
+// Return the median value of the distribution of abs(a[i]-Median)
+//
+Double_t MMath::MedianDev(Long64_t n, const Double_t *a)
+{
+    return MedianDevImp(n, a);
+}
+
+// ------------------------------------------------------------------------
+//
+// Return the median value of the distribution of abs(a[i]-Median)
+//
+Double_t MMath::MedianDev(Long64_t n, const Long_t *a)
+{
+    return MedianDevImp(n, a);
+}
+
+// ------------------------------------------------------------------------
+//
+// Return the median value of the distribution of abs(a[i]-Median)
+//
+Double_t MMath::MedianDev(Long64_t n, const Long64_t *a)
+{
+    return MedianDevImp(n, a);
+}
+
 // --------------------------------------------------------------------------
 //

trunk/MagicSoft/Mars/mbase/MMath.h

@@ -18 +18 @@
 
     Double_t GaussProb(Double_t x, Double_t sigma=1, Double_t mean=0);
+
+    template <class Size, class Element> Double_t MedianDevImp(Size n, const Element *a);
+    Double_t  MedianDev(Long64_t n, const Short_t  *a);
+    Double_t  MedianDev(Long64_t n, const Int_t    *a);
+    Double_t  MedianDev(Long64_t n, const Float_t  *a);
+    Double_t  MedianDev(Long64_t n, const Double_t *a);
+    Double_t  MedianDev(Long64_t n, const Long_t   *a);
+    Double_t  MedianDev(Long64_t n, const Long64_t *a);
 
     Double_t Significance(Double_t s, Double_t b);