Changeset 862

Timestamp:

07/10/01 14:19:56 (23 years ago)

Author:

tbretz

Message:

*** empty log message ***

File:

• trunk/MagicSoft/Mars/mmontecarlo/MMcThresholdCalc.cc (modified) (3 diffs)

trunk/MagicSoft/Mars/mmontecarlo/MMcThresholdCalc.cc

@@ -87 +87 @@
 }
 
+// --------------------------------------------------------------------------
+//
+// 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 +142 @@
 }
 
+// --------------------------------------------------------------------------
+//
+// 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 +155 @@
 
     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);