Changeset 7867 for trunk/MagicSoft/Mars/mhist/MHCamera.cc

Timestamp:

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

Author:

tbretz

Message:

*** empty log message ***

File:

• trunk/MagicSoft/Mars/mhist/MHCamera.cc (modified) (2 diffs)

trunk/MagicSoft/Mars/mhist/MHCamera.cc

@@ -84 +84 @@
 
 #include "MArrayD.h"
+#include "MMath.h"    // MMath::GaussProb
 
 #define kItemsLegend 48 // see SetPalette(1,0)
@@ -312 +313 @@
     for (Int_t idx=0; idx<fNcells-2; idx++)
         arr[idx] ? SetUsed(idx) : ResetUsed(idx);
+}
+
+// ------------------------------------------------------------------------
+//
+// Return the median value on the y-axis (profile option is correctly taken
+// into account)
+//
+Stat_t MHCamera::GetMedian() const
+{
+    // Just for speed reasons
+    if (!TestBit(kProfile))
+        return TMath::Median(GetSize()-2, GetArray()+1);
+
+    // Copy profiled data into new array (FIXME: Should we take errors into account?)
+    TArrayD arr(fNcells-2);
+    for (int i=1; i<fNcells-1; i++)
+        arr[i-1] = GetBinContent(i);
+
+    // return Median of the profile data
+    return TMath::Median(arr.GetSize(), arr.GetArray());
+}
+
+// ------------------------------------------------------------------------
+//
+// Return the median value (divided by MMath::GausProb(1.0)) of the
+// distribution of abs(y[i]-Median). This is my Median equivalent of the RMS
+//
+Stat_t MHCamera::GetMedianDev() const
+{
+    // Just for speed reasons
+    if (!TestBit(kProfile))
+        return MMath::MedianDev(GetSize()-2, GetArray()+1);
+
+    // Copy profiled data into new array (FIXME: Should we take errors into account?)
+    TArrayD arr(fNcells-2);
+    for (int i=1; i<fNcells-1; i++)
+        arr[i-1] = GetBinContent(i);
+
+    // return MedianDev of the profile data
+    return MMath::MedianDev(arr.GetSize(), arr.GetArray())/MMath::GaussProb(1.0);
 }