Changeset 7868

Timestamp:

08/09/06 10:39:45 (18 years ago)

Author:

tbretz

Message:

*** empty log message ***

File:

• trunk/MagicSoft/Mars/mbase/MMath.cc (modified) (1 diff)

trunk/MagicSoft/Mars/mbase/MMath.cc

@@ -168 +168 @@
 Double_t MMath::MedianDevImp(Size n, const Element *a)
 {
+    // Sanity check
+    if (n <= 0 || !a)
+        return 0;
+
     // 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
+    Double_t arr[n];
     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 TMath::Median(n, arr);
+    /*
+    // Sort distribution
+    Long64_t idx[n];
+    TMath::SortImp(n, arr, idx, kTRUE);
+
+    // Define where to divide
+    const Int_t div = TMath::Nint(n*MMath::GaussProb(1.0));
+
+    // Calculate result
+    Double_t dev = TMath::KOrdStat(n, arr, div, idx);
+    if (n%2 == 0)
+    {
+        dev += TMath::KOrdStat(n, arr, div-1, idx);
+        dev /= 2;
+    }
+
+    return dev;
+*/
 }