Index: trunk/Mars/mpedestal/MPedestalSubtractedEvt.cc
===================================================================
--- trunk/Mars/mpedestal/MPedestalSubtractedEvt.cc	(revision 11873)
+++ trunk/Mars/mpedestal/MPedestalSubtractedEvt.cc	(revision 12625)
@@ -32,4 +32,6 @@
 #include "MPedestalSubtractedEvt.h"
 
+#include <TMath.h>
+
 #include "MLogManip.h"
 
@@ -127,4 +129,29 @@
 }
 
+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,4 +233,28 @@
     }
     *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;
 }
 
Index: trunk/Mars/mpedestal/MPedestalSubtractedEvt.h
===================================================================
--- trunk/Mars/mpedestal/MPedestalSubtractedEvt.h	(revision 11873)
+++ trunk/Mars/mpedestal/MPedestalSubtractedEvt.h	(revision 12625)
@@ -93,4 +93,5 @@
     }
 
+    void GetStat(const Int_t idx, Float_t &mean, Float_t &rms) const;
 
     Int_t GetSaturation(const Int_t pixidx, Int_t limit) const
@@ -116,5 +117,5 @@
     void Print(Option_t *o="") const;
 
-    Bool_t GetPixelContent(Double_t &val, Int_t idx, const MGeomCam &cam, Int_t type=0) const { return kTRUE; }
+    Bool_t GetPixelContent(Double_t &val, Int_t idx, const MGeomCam &cam, Int_t type=0) const;
     void   DrawPixelContent(Int_t num) const { }