Changeset 12920 for trunk/Mars

Timestamp:

02/22/12 19:02:50 (13 years ago)

Author:

tbretz

Message:

Added a function to return the median of the whole camera.

Location:

trunk/Mars/mgui

Files:

• MCamEvent.cc (modified) (4 diffs)
• MCamEvent.h (modified) (1 diff)

trunk/Mars/mgui/MCamEvent.cc

@@ -43 +43 @@
 
 #include "MGeomCam.h"
+#include "MArrayD.h"
 
 ClassImp(MCamEvent);
@@ -60 +61 @@
 // --------------------------------------------------------------------------
 //
-// Return the mean of all corresponding GetPixelContent
+// Return the median of all corresponding GetPixelContent
 //
-Double_t MCamEvent::GetCameraMean(const MGeomCam &cam, Int_t type) const
+Double_t MCamEvent::GetCameraMedian(const MGeomCam &cam, Int_t type) const
 {
-    return GetCameraStat(cam, type)[0];
-    /*
-    Int_t    num  = 0;
-    Double_t mean = 0;
+    Int_t    num = 0;
+    Double_t sum = 0;
+    Double_t sq  = 0;
+
+    MArrayD arr(cam.GetNumPixels());
     for (unsigned int i=0; i<cam.GetNumPixels(); i++)
     {
@@ -74 +76 @@
             continue;
 
-        mean += val;
+        arr[i] = val;
         num ++;
     }
-    return num == 0 ? 0 : mean/num;
-    */
+
+    return TMath::Median(num, arr.GetArray());
+}
+
+// --------------------------------------------------------------------------
+//
+// Return the mean of all corresponding GetPixelContent
+//
+Double_t MCamEvent::GetCameraMean(const MGeomCam &cam, Int_t type) const
+{
+    return GetCameraStat(cam, type)[0];
 }
 
@@ -88 +99 @@
 {
     return GetCameraStat(cam, type)[1];
-    /*
-    Int_t    num  = 0;
-    Double_t sum = 0;
-    Double_t sq  = 0;
-
-    for (unsigned int i=0; i<cam.GetNumPixels(); i++)
-    {
-        Double_t val;
-        if (!GetPixelContent(val, i, cam, type))
-            continue;
-
-        sum += val;
-        sq  += val;
-        num ++;
-    }
-
-    if (num==0)
-        return 0;
-
-    sum /= num;
-    sq  /= num;
-
-    return TMath::Sqrt(sq-sum*sum);
-    */
 }
 

trunk/Mars/mgui/MCamEvent.h

@@ -20 +20 @@
 
     virtual Double_t GetCameraMean(const MGeomCam &cam, Int_t type=0) const;
+    virtual Double_t GetCameraMedian(const MGeomCam &cam, Int_t type=0) const;
     virtual Double_t GetCameraRMS(const MGeomCam &cam, Int_t type=0) const;
     virtual TArrayD  GetCameraStat(const MGeomCam &cam, Int_t type=0) const;