Index: /trunk/Mars/mjobs/MJSimulation.cc
===================================================================
--- /trunk/Mars/mjobs/MJSimulation.cc	(revision 17662)
+++ /trunk/Mars/mjobs/MJSimulation.cc	(revision 17663)
@@ -350,15 +350,4 @@
 Bool_t MJSimulation::Process(const MArgs &args, const MSequence &seq)
 {
-    /*
-    if (!fSequence.IsValid())
-    {
-        *fLog << err << "ERROR - Sequence invalid!" << endl;
-        return kFALSE;
-    }
-    */
-
-//    if (!HasWritePermission(CombinePath(fPathOut, GetOutFile(seq))))
-//        return kFALSE;
-
     *fLog << inf;
     fLog->Separator(GetDescriptor());
@@ -381,13 +370,10 @@
 
     // --------------------------------------------------------------------------------
-
     // Setup Parlist
     MParList plist;
     plist.AddToList(this); // take care of fDisplay!
-
     // setup TaskList
     MTaskList tasks;
     plist.AddToList(&tasks);
-
     // --------------------------------------------------------------------------------
 
@@ -444,15 +430,4 @@
     header.SetObservation("On", "MonteCarlo");
     plist.AddToList(&header);
-    // ++++++++ FIXME FIXME FIXME +++++++++++++
-
-    /*
-    MPedestalCam pedcam;
-    pedcam.Init(geomcam.GetNumPixels());
-    for (UInt_t i=0; i<geomcam.GetNumPixels(); i++)
-        pedcam[i].Set(128./header.GetScale(), 22.5/header.GetScale());
-    plist.AddToList(&pedcam);
-    */
-
-    // -------------------------------------------------------------------
 
     MCorsikaRead read;
@@ -677,4 +652,5 @@
     write3af.AddContainer("IntendedPulsePos", "RunHeaders");
     write3af.AddContainer("MRawEvtData",      "Events");
+    write3af.AddContainer("MTruePhotonsPerPixelCont", "Events");
     
     if (fOperationMode == kModeData){
Index: /trunk/Mars/msimcamera/MSimCamera.cc
===================================================================
--- /trunk/Mars/msimcamera/MSimCamera.cc	(revision 17662)
+++ /trunk/Mars/msimcamera/MSimCamera.cc	(revision 17663)
@@ -47,4 +47,6 @@
 #include "MLogManip.h"
 
+#include "MTruePhotonsPerPixelCont.h"
+
 #include "MSpline3.h"
 #include "MParSpline.h"
@@ -149,4 +151,11 @@
     {
         *fLog << err << "CrosstalkCoeffParam [MParameterD] not found... aborting." << endl;
+        return kFALSE;
+    }
+
+    fTruePhotons = (MTruePhotonsPerPixelCont*)pList->FindCreateObj("MTruePhotonsPerPixelCont");
+    if (!fTruePhotons)
+    {
+        *fLog << err << "MTruePhotonsPerPixelCont not found... aborting." << endl;
         return kFALSE;
     }
@@ -344,5 +353,33 @@
 
         if (ph.GetPrimary()!=MMcEvt::kNightSky && ph.GetPrimary()!=MMcEvt::kArtificial)
+        {
             tot += ph.GetWeight();
+            (*fTruePhotons->cherenkov_photons_weight)[idx] += ph.GetWeight();
+            (*fTruePhotons->cherenkov_photons_number)[idx] += 1;
+
+            (*fTruePhotons->cherenkov_arrival_time_mean)[idx] += t;
+            (*fTruePhotons->cherenkov_arrival_time_variance)[idx] += t*t;
+
+            if (ph.GetPrimary()==MMcEvt::kMUON)
+            {
+                (*fTruePhotons->muon_cherenkov_photons_weight)[idx] += ph.GetWeight();
+                (*fTruePhotons->muon_cherenkov_photons_number)[idx] += 1;
+            }
+
+            // find min
+            if (t < (*fTruePhotons->cherenkov_arrival_time_min)[idx] )
+            {
+                (*fTruePhotons->cherenkov_arrival_time_min)[idx] = t;
+            }
+            // find max
+            if (t > (*fTruePhotons->cherenkov_arrival_time_max)[idx] )
+            {
+               (*fTruePhotons->cherenkov_arrival_time_max)[idx] = t;
+            }
+        }
+        else
+        {
+            (*fTruePhotons->noise_photons_weight)[idx] += ph.GetWeight();
+        }
 
         // Sorry, the name "pedestal" is misleading here
@@ -352,4 +389,13 @@
         // === FIXME === FIXME === FIXME === Frequency!!!!
         (*fCamera)[idx].AddPulse(*fSpline, t, ph.GetWeight()*gain);
+    }
+
+    for (unsigned int i=0 ; i < 1440 ; i++)
+    {
+        float number = (*fTruePhotons->cherenkov_photons_number)[i];
+        (*fTruePhotons->cherenkov_arrival_time_mean)[i] /= number;
+        float mean = (*fTruePhotons->cherenkov_arrival_time_mean)[i];
+        float sum_tt = (*fTruePhotons->cherenkov_arrival_time_variance)[i];
+        (*fTruePhotons->cherenkov_arrival_time_variance)[i] = (sum_tt / number - mean*mean) /(number - 1);
     }
 
Index: /trunk/Mars/msimcamera/MSimCamera.h
===================================================================
--- /trunk/Mars/msimcamera/MSimCamera.h	(revision 17662)
+++ /trunk/Mars/msimcamera/MSimCamera.h	(revision 17663)
@@ -5,4 +5,6 @@
 #include "MTask.h"
 #endif
+
+#include "MArrayF.h"
 
 class MMcEvt;
@@ -13,4 +15,6 @@
 class MAnalogChannels;
 class MPedestalCam;
+class MArrayF;
+class MTruePhotonsPerPixelCont;
 
 class MSpline3;
@@ -31,4 +35,5 @@
 
     MParameterD       *fCrosstalkCoeffParam;
+    MTruePhotonsPerPixelCont    *fTruePhotons;  //! Container to store the number of photons per pixel
 
     const MSpline3    *fSpline;          // Pulse Shape