Changeset 861 for trunk/MagicSoft/Mars/mmontecarlo/MMcThresholdCalc.cc

Timestamp:

07/10/01 12:41:13 (23 years ago)

Author:

jlopez

Message:

*** empty log message ***

File:

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

trunk/MagicSoft/Mars/mmontecarlo/MMcThresholdCalc.cc

@@ -24 +24 @@
 \* ======================================================================== */
 
-/////////////////////////////////////////////////////////////////////////////
-//                                                                         //
-//  MMcThresholdCalc                                                       //
-//                                                                         //
-//  Input Containers:                                                      //
-//   MMcEvt                                                                //
-//                                                                         //
-//  Output Containers:                                                     //
-//   MMcTrig;*, MHMcEnergy;*                                               //
-//                                                                         //
+///////////////////////////////////////////////////////////////////////////
+//
+//  MMcThresholdCalc
+//                                                                         
+//  Input Containers:
+//   MMcEvt, MMcTrig;*
+//
+//  Output Containers:
+//   MHMcEnergies
+//
 /////////////////////////////////////////////////////////////////////////////
 
@@ -52 +52 @@
 ClassImp(MMcThresholdCalc)
 
+// --------------------------------------------------------------------------
+//
+// Default Constructor.
+//
 MMcThresholdCalc::MMcThresholdCalc(const UInt_t dim, const char* name,
                                    const char* title) : fDimension(dim)
@@ -58 +62 @@
     *fTitle = title ? title : "Task to calculate the energy threshold from Monte Carlo";
 
+    // Arrays of MMcTrig* and MHMcEnergy* are created in order to be 
+    // able to work with root files with several trigger conditions.
     fMcTrig     = new MMcTrig*[fDimension];
     fHMcEnergy  = new MHMcEnergy*[fDimension];
@@ -66 +72 @@
 }
 
+// -------------------------------------------------------------------------
+//
+// Default Destructor.
+//
 MMcThresholdCalc::~MMcThresholdCalc()
 {
@@ -81 +91 @@
     // connect Monte Carlo data with this task
 
+    // This task has into accout if the root file has one trigger
+    // condition (MMcTrig) or severl of them (MMcTrig;#.)
+
     fMcEvt = (MMcEvt*)pList->FindObject("MMcEvt");
     if (!fMcEvt)
@@ -93 +106 @@
     {
         if (fDimension>1)
-            sprintf(auxname+7, ";%i", i+1);
+            sprintf(auxname+7, ";%i.", i+1);
 
         fMcTrig[i] = (MMcTrig*)pList->FindObject(auxname);
@@ -110 +123 @@
     for (unsigned int i=0; i<fDimension; i++)
     {
-        if (fDimension>1)
-            sprintf(auxname+10, ";%i", i+1);
+        if (fDimension>1&&i!=0)
+            sprintf(auxname+10, ";%i", i);
 
         fHMcEnergy[i] = (MHMcEnergy*)pList->FindObject(auxname);
@@ -119 +132 @@
         *fLog << dbginf << "'" << auxname << "' not found in list... creating." << endl;
 
-        fHMcEnergy[i] = new MHMcEnergy(fDimension>1 ? i+1 : 0);
+        fHMcEnergy[i] = new MHMcEnergy(fDimension>1&&i!=0 ? i : 0);
         fMustDelete[i] = kTRUE;
         pList->AddToList(fHMcEnergy[i]);
@@ -129 +142 @@
 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);
@@ -147 +165 @@
 {
     // 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);
@@ -152 +171 @@
     for (unsigned int i=0; i<fDimension; i++)
     {
-        MHMcEnergy &h = *fHMcEnergy[i];
+      Float_t peak;
+      Float_t sigma;
 
-        const Float_t peak  = h.GetGaussPeak();
-        const Float_t sigma = h.GetGaussSigma();
-
-        h.Fit(1, 3);
-        h.Fit(peak - 2.   *sigma, peak + 2.   *sigma);
-        h.Fit(peak - sqrt2*sigma, peak + sqrt2*sigma);
+      fHMcEnergy[i]->Fit(1, 3);
+      
+      peak  = fHMcEnergy[i]->GetGaussPeak();
+      sigma = fHMcEnergy[i]->GetGaussSigma();
+      fHMcEnergy[i]->Fit(peak - 2.   *sigma, peak + 2.   *sigma);
+        
+      peak  = fHMcEnergy[i]->GetGaussPeak();
+      sigma = fHMcEnergy[i]->GetGaussSigma();
+      fHMcEnergy[i]->Fit(peak - sqrt2*sigma, peak + sqrt2*sigma);
     }
     return kTRUE;