Index: /trunk/MagicSoft/Mars/mmontecarlo/MMcThresholdCalc.cc
===================================================================
--- /trunk/MagicSoft/Mars/mmontecarlo/MMcThresholdCalc.cc	(revision 861)
+++ /trunk/MagicSoft/Mars/mmontecarlo/MMcThresholdCalc.cc	(revision 862)
@@ -87,9 +87,11 @@
 }
 
+// --------------------------------------------------------------------------
+//
+// connect Monte Carlo data with this task
+//
 Bool_t MMcThresholdCalc::PreProcess(MParList* pList)
 {
-    // connect Monte Carlo data with this task
-
-    // This task has into accout if the root file has one trigger
+    // This task takes into accout if the root file has one trigger
     // condition (MMcTrig) or severl of them (MMcTrig;#.)
 
@@ -140,11 +142,12 @@
 }
 
+// --------------------------------------------------------------------------
+//
+// The histograms are filled with log10 of the energy for triggered
+// events and weighted with 1/E because it is needed the dN/dE vs. logE
+// distribution to get the energy threshold.
+//
 Bool_t MMcThresholdCalc::Process()
 {
-
-    // The histograms are filled with log10 of the energy for triggered
-    // events and weighted with 1/E because it is needed the dN/dE vs. logE
-    // distribution to get the energy threshold.
-
     const Float_t energy   = fMcEvt->GetEnergy();
     const Float_t lg10     = log10(energy);
@@ -152,19 +155,17 @@
 
     for (unsigned int i=0; i<fDimension; i++)
-    {
-        if (fMcTrig[i]->GetFirstLevel()<=0)
-            continue;
-
-        fHMcEnergy[i]->Fill(lg10, reciproc);
-    }
+        if (fMcTrig[i]->GetFirstLevel()>0)
+            fHMcEnergy[i]->Fill(lg10, reciproc);
 
     return kTRUE;
 }
 
+// --------------------------------------------------------------------------
+//
+// fit the energy distribution to get the threshold
+// Some iterations are done to be sure the fit parameters converge.
+//
 Bool_t MMcThresholdCalc::PostProcess()
 {
-    // fit the energy distribution to get the threshold
-    // Some iterations are done to be sure the fit parameters converge.
-
     const Float_t sqrt2 = sqrt(2);