Changeset 12625 for trunk/Mars/mpedestal/MPedestalSubtractedEvt.cc

Timestamp:

11/23/11 14:09:26 (13 years ago)

Author:

tbretz

Message:

Added possibility to return mean and rms in GetPixelContant

File:

• trunk/Mars/mpedestal/MPedestalSubtractedEvt.cc (modified) (3 diffs)

trunk/Mars/mpedestal/MPedestalSubtractedEvt.cc

@@ -32 +32 @@
 #include "MPedestalSubtractedEvt.h"
 
+#include <TMath.h>
+
 #include "MLogManip.h"
 
@@ -127 +129 @@
 }
 
+void MPedestalSubtractedEvt::GetStat(const Int_t idx, Float_t &mean, Float_t &rms) const
+{
+    if (fNumSamples<20)
+        return;
+
+    // Get pointer to first slice to be considered
+    Float_t const *sam = GetSamples(idx);
+    Float_t const *beg = sam;
+
+    Double_t sum  = 0;
+    Double_t sum2 = 0;
+
+    for (const Float_t *ptr=beg+10; ptr<sam+fNumSamples-10; ptr++)
+    {
+        sum  += *ptr;
+        sum2 += *ptr * *ptr;
+    }
+
+    sum  /= fNumSamples-20;
+    sum2 /= fNumSamples-20;
+
+    mean = sum;
+    rms  = TMath::Sqrt(sum2 - sum * sum);
+}
+
 // --------------------------------------------------------------------------
 //
@@ -206 +233 @@
     }
     *fLog << endl;
+}
+
+Bool_t MPedestalSubtractedEvt::GetPixelContent(Double_t &val, Int_t idx, const MGeomCam &cam, Int_t type) const
+{
+    switch (type)
+    {
+    case 0:
+    case 1:
+        {
+            Float_t mean, rms;
+            GetStat(idx, mean, rms);
+            val = type==0 ? mean : rms;
+        }
+        break;
+    case 2:
+        {
+            Float_t flt;
+            GetMaxPos(idx, flt);
+            val = flt;
+        }
+        break;
+    }
+
+    return kTRUE;
 }