Changeset 2272 for trunk/MagicSoft/Mars/macros

Timestamp:

07/10/03 12:34:10 (21 years ago)

Author:

wittek

Message:

*** empty log message ***

File:

• trunk/MagicSoft/Mars/macros/CT1EgyEst.C (modified) (8 diffs)

trunk/MagicSoft/Mars/macros/CT1EgyEst.C

@@ -33 +33 @@
 #include "MHMatrix.h"
 #include "MEnergyEstParam.h"
+#include "MEnergyEstParamDanielMkn421.h"
 #include "MMatrixLoop.h"
 #include "MChisqEval.h"
@@ -138 +139 @@
   Optimize.SetOutFile(paramout);
 
-  /*
+  
   MFCT1SelFinal filterhadrons;
   filterhadrons.SetHadronnessName(hadronnessName);
@@ -144 +145 @@
   filterhadrons.SetInverted();
   Optimize.SetEventFilter(&filterhadrons);
-  */
+  
 
   Optimize.SetNevents(howMany);
@@ -205 +206 @@
   //
 
-  MEnergyEstParam eest2(hilName);
+  //MEnergyEstParam eest2(hilName);
+  //eest2.Add(hilSrcName);
+
+  //eest2.SetCoeffA(parA);
+  //eest2.SetCoeffB(parB);
+
+  // estimate energy using Daniel's parameters
+  MEnergyEstParamDanielMkn421 eest2(hilName);
   eest2.Add(hilSrcName);
-
-  eest2.SetCoeffA(parA);
-  eest2.SetCoeffB(parB);
 
 
@@ -238 +243 @@
   // Select events
   //
-
-  /*
+  
   cout << "Select events with hadronness < " << maxhadronness 
         << " and |alpha| < " << maxalpha << endl; 
@@ -247 +251 @@
 
   MContinue cont(&hcut2);
-  */
-
-  parlist2.AddToList(&tlist2);
-
-  //********************************
-  // Entries in MTaskList
-
-  tlist2.AddToList(&read2);
-  //  tlist2.AddToList(&cont);
-  tlist2.AddToList(&eest2);
-
-  //
+  
+  //-----------------------------------------------
   // Create Object MHMcEnergyMigration containing useful histograms,
   // and task MHMcEnergyMigration to fill them:
   //
 
-  MHMcEnergyMigration mighist;
-
-  parlist2.AddToList(&mighist);
-
-  MFillH migfill(&mighist, "MMcEvt");
-
-  tlist2.AddToList(&migfill);
-
-  parlist2.AddToList(binsE);
-  parlist2.AddToList(binsTheta);
+  //MHMcEnergyMigration mighist;
+
+  MFillH migfill("MHMcEnergyMigration", "MMcEvt");
 
   MBinning BinningDE("BinningDE");
@@ -279 +266 @@
   BinningDE.SetEdges(60, -1.2, 1.2);
   BinningImpact.SetEdges(50, 0., 400.);
+
+  //********************************
+  // Entries in MParList
+
+  parlist2.AddToList(&tlist2);
+  //parlist2.AddToList(&mighist);
+  parlist2.AddToList(binsE);
+  parlist2.AddToList(binsTheta);
   parlist2.AddToList(&BinningDE);
   parlist2.AddToList(&BinningImpact);
 
+  //********************************
+  // Entries in MTaskList
+
+  tlist2.AddToList(&read2);
+  tlist2.AddToList(&cont);
+  tlist2.AddToList(&eest2);
+  tlist2.AddToList(&migfill);
+
+  //-----------------------------------------------
   cout << "Event loop was setup" << endl; 
   MProgressBar bar;
@@ -291 +295 @@
   if (!evtloop2.Eventloop())
     return;
+
+  tlist2.PrintStatistics(0, kTRUE);
+
+  MHMcEnergyMigration &mighist = 
+    *(MHMcEnergyMigration*)parlist2->FindObject("MHMcEnergyMigration");
+  if (!&mighist)
+  {
+    cout << "CT1EgyEst.C : object 'MHMcEnergyMigration' not found ... aborting"
+         << endl;
+    return;
+  }
 
   TString paramout(outPath);