Index: trunk/MagicSoft/Mars/manalysis/MHPedestalCam.cc
===================================================================
--- trunk/MagicSoft/Mars/manalysis/MHPedestalCam.cc	(revision 3768)
+++ trunk/MagicSoft/Mars/manalysis/MHPedestalCam.cc	(revision 3770)
@@ -312,4 +312,8 @@
   }
 
+
+  fLambda.Set(npixels);
+  fLambdaErr.Set(npixels);
+
   return kTRUE;
 
@@ -442,6 +446,8 @@
                   MBadPixelsPix::kLoGainOscillating);
 
+
   return kTRUE;
   
+
 }
 
@@ -476,8 +482,17 @@
 // 12: Fourier spectrum not OK
 //
+// Poisson fit:
+// ============
+//
+// 13: Poisson fit lambda
+// 14: Poisson fit lambda error 
+//
 Bool_t MHPedestalCam::GetPixelContent(Double_t &val, Int_t idx, const MGeomCam &cam, Int_t type) const
 {
 
   if (fHiGainArray->GetSize() <= idx)
+    return kFALSE;
+
+  if ((*this)[idx].IsExcluded())
     return kFALSE;
 
@@ -540,4 +555,14 @@
         val = 1.;
       break;
+    case 13:
+      if (!(*this)[idx].IsPoissonFitOK())
+        return kFALSE;
+      val = fLambda[idx];
+      break;
+    case 14:
+      if (!(*this)[idx].IsPoissonFitOK())
+        return kFALSE;
+      val = fLambdaErr[idx];
+      break;
     default:
       return kFALSE;
Index: trunk/MagicSoft/Mars/manalysis/MHPedestalPix.cc
===================================================================
--- trunk/MagicSoft/Mars/manalysis/MHPedestalPix.cc	(revision 3768)
+++ trunk/MagicSoft/Mars/manalysis/MHPedestalPix.cc	(revision 3770)
@@ -42,4 +42,7 @@
 #include "MHPedestalPix.h"
 
+#include "MLog.h"
+#include "MLogManip.h"
+
 #include <TH1.h>
 
@@ -48,6 +51,6 @@
 using namespace std;
 //
-const Int_t   MHPedestalPix::fgChargeNbins    = 450 ;
-const Axis_t  MHPedestalPix::fgChargeFirst    = -0.5;
+const Int_t   MHPedestalPix::fgChargeNbins    = 500;
+const Axis_t  MHPedestalPix::fgChargeFirst    = -49.5;
 const Axis_t  MHPedestalPix::fgChargeLast     = 449.5;
 // --------------------------------------------------------------------------
@@ -103,7 +106,4 @@
 {
 
-  if (fMean == fMeanErr == fSigma == fSigmaErr == 0.)
-    return;
-
   //
   // One never knows...
@@ -113,5 +113,5 @@
   
   const Float_t sqslices = TMath::Sqrt(fNSlices);
-  
+
   SetMean     ( GetMean()    / fNSlices  );
   //