Changeset 2241

Timestamp:

06/26/03 20:57:18 (22 years ago)

Author:

moralejo

Message:

*** empty log message ***

Location:

trunk/MagicSoft/Mars

Files:

• Changelog (modified) (1 diff)
• macros/CT1EgyEst.C (modified) (9 diffs)

trunk/MagicSoft/Mars/Changelog

@@ -1 +1 @@
                                                  -*-*- END OF LINE -*-*-
+
+ 2003/06/26: Abelardo Moralejo
+  * macros/CT1EgyEst.C
+    - Added to the argument list of CT1EEst the binnings in theta
+      and energy. Commented out the MLog.h and MLogManip.h, headers,
+      otherwise it does not run (don't know why). Commented out filters 
+      (my test file does not contain hadronness), please uncomment 
+      them to run on files with hadronness info.
 
  2003/06/26: Abelardo Moralejo

trunk/MagicSoft/Mars/macros/CT1EgyEst.C

@@ -35 +35 @@
 #include "MMatrixLoop.h"
 #include "MChisqEval.h"
+
+/*
 #include "MLog.h"
 #include "MLogManip.h"
+*/
 
 void CT1EgyEst()
@@ -61 +64 @@
   Float_t maxdist       = 1.;
 
+  MBinning BinningE("BinningE");
+  MBinning BinningTheta("BinningTheta");
+  BinningE.SetEdgesLog(50, 500, 50000.);
+  BinningTheta.SetEdges(5, 0., 50.);
+
   CT1EEst(inPath,  fileNameIn, 
-            outPath, energyParName, 
-            hilName, hilSrcName,    hadronnessName, 
-            howMany, maxhadronness, maxalpha, maxdist);
+          outPath, energyParName, 
+          hilName, hilSrcName,    hadronnessName, 
+          howMany, maxhadronness, maxalpha, maxdist, 
+          &BinningE, &BinningTheta);
 }
 
@@ -84 +93 @@
 
 void CT1EEst(TString inPath,  TString fileNameIn, 
-               TString outPath, TString energyParName,
-               TString hilName, TString hilSrcName, TString hadronnessName,
-               Int_t howMany, 
-               Float_t maxhadronness, Float_t maxalpha, Float_t maxdist)
+             TString outPath, TString energyParName,
+             TString hilName, TString hilSrcName, TString hadronnessName,
+             Int_t howMany, 
+             Float_t maxhadronness, Float_t maxalpha, Float_t maxdist,
+             MBinning* binsE, MBinning* binsTheta)
 {
   cout << "================================================================"
@@ -125 +135 @@
   Optimize.SetOutFile(paramout);
 
+  /*
   MFCT1SelFinal filterhadrons;
   filterhadrons.SetHadronnessName(hadronnessName);
@@ -130 +141 @@
   filterhadrons.SetInverted();
   Optimize.SetEventFilter(&filterhadrons);
+  */
 
   Optimize.SetNevents(howMany);
@@ -223 +235 @@
   // Select events
   //
+
+  /*
   cout << "Select events with hadronness < " << maxhadronness 
         << " and |alpha| < " << maxalpha << endl; 
@@ -230 +244 @@
 
   MContinue cont(&hcut2);
+  */
 
   parlist2.AddToList(&tlist2);
@@ -237 +252 @@
 
   tlist2.AddToList(&read2);
-  tlist2.AddToList(&cont);
+  //  tlist2.AddToList(&cont);
   tlist2.AddToList(&eest2);
 
@@ -253 +268 @@
   tlist2.AddToList(&migfill);
 
-  MBinning BinningE("BinningE");
-  MBinning BinningTheta("BinningTheta");
-  BinningE.SetEdgesLog(50, 500, 50000.);
-  BinningTheta.SetEdges(5, 0., 50.);
-  
-  parlist2.AddToList(&BinningE);
-  parlist2.AddToList(&BinningTheta);
+  parlist2.AddToList(binsE);
+  parlist2.AddToList(binsTheta);
 
   MBinning BinningDE("BinningDE");