correct the cluster rejection by energy range and number of cells cut

[u/mrichter/AliRoot.git] / PWG4 / macros / electrons / anaJete.C

diff --git a/PWG4/macros/electrons/anaJete.C b/PWG4/macros/electrons/anaJete.C
index 220ace96c5b098bb341735b758642a695587e589..02e2382515c480cb0ce979cc00fb09dfa4eb4eb8 100755 (executable)
--- a/PWG4/macros/electrons/anaJete.C
+++ b/PWG4/macros/electrons/anaJete.C
@@ -36,12 +36,12 @@ char * kmydataset = "/COMMON/COMMON/LHC09a4_run8101X";
 //Put name of file containing xsection \r
 //Put number of events per ESD file\r
 //This is an specific case for normalization of Pythia files.\r
-const Bool_t kGetXSectionFromFileAndScale = kTRUE;\r
+Bool_t kGetXSectionFromFileAndScale = kTRUE;\r
 const char * kXSFileName = "pyxsec.root";\r
 const Int_t kNumberOfEventsPerFile = 200; \r
 //---------------------------------------------------------------------------\r
 \r
-const Bool_t kMC = kTRUE; //With real data kMC = kFALSE\r
+Bool_t kMC = kTRUE; //With real data kMC = kFALSE\r
 TString kInputData = "ESD";//ESD, AOD, MC\r
 TString kTreeName = "esdTree";\r
 //const   Bool_t kMergeAODs = kTRUE; //uncomment for AOD merging\r
@@ -65,7 +65,11 @@ void anaJete()
 \r
   //Process environmental variables from command line:\r
   ProcessEnvironment();        \r
-  printf("Final    Variables: kInputData %s, mode %d, config2 %s, config3 %s, sevent %d\n",kInputData.Data(),mode,sconfig2,sconfig3,sevent);\r
+  printf("Final    Variables: kInputData %s, kMC %d, mode %d, config2 %s, config3 %s, sevent %d\n",kInputData.Data(),kMC,mode,sconfig2,sconfig3,sevent);\r
+\r
+  if(!kMC) {\r
+    kGetXSectionFromFileAndScale = kFALSE;\r
+  }\r
 \r
   //--------------------------------------------------------------------\r
   // Load analysis libraries\r
@@ -112,7 +116,7 @@ void anaJete()
       //plugin->AddRunNumber(30010); //dummy\r
       plugin->AddDataFile("mycollect.xml");\r
       plugin->SetGridWorkingDir("work3");\r
-      plugin->SetAdditionalLibs("anaJet.C ConfigJetAnalysisFastJet.C ConfigAnalysisElectron.C ANALYSIS.par ANALYSISalice.par AOD.par ESD.par STEERBase.par JETAN.par FASTJETAN.par");\r
+      plugin->SetAdditionalLibs("anaJet.C ConfigJetAnalysisFastJet.C ConfigAnalysisElectron.C ANALYSIS.par ANALYSISalice.par AOD.par EMCALUtils.par ESD.par PHOSUtils.par STEERBase.par JETAN.par FASTJETAN.par");\r
       plugin->SetJDLName("anaJet.jdl");\r
       plugin->SetExecutable("anaJet.sh");\r
       plugin->SetOutputFiles("histos.root");\r
@@ -337,6 +341,8 @@ void  LoadLibraries(const anaModes mode) {
       SetupPar("AOD");\r
       SetupPar("ANALYSIS");\r
       SetupPar("ANALYSISalice");\r
+      SetupPar("PHOSUtils");\r
+      SetupPar("EMCALUtils");\r
       if( kDoJetTask ){\r
         cerr<<"Now Loading JETAN"<<endl;\r
         SetupPar("JETAN");\r
@@ -358,6 +364,8 @@ void  LoadLibraries(const anaModes mode) {
       gSystem->Load("libAOD");\r
       gSystem->Load("libANALYSIS");\r
       gSystem->Load("libANALYSISalice");\r
+      gSystem->Load("libPHOSUtils");\r
+      gSystem->Load("libEMCALUtils");\r
       if( kDoJetTask ){\r
         gSystem->Load("libJETAN");\r
         gSystem->Load("libFASTJETAN");\r
@@ -384,7 +392,7 @@ void  LoadLibraries(const anaModes mode) {
     //TProof::Reset("myproofname",kTRUE);\r
     gEnv->SetValue("XSec.GSI.DelegProxy","2");   \r
     //TProof::Mgr(myproofname)->ShowROOTVersions();\r
-    //TProof::Mgr(myproofname)->SetROOTVersion("v5-23-04");\r
+    //TProof::Mgr(myproofname)->SetROOTVersion("v5-24-00");\r
     TProof::Open(myproofname);\r
 \r
     // gProof->ClearPackages();\r
@@ -394,10 +402,14 @@ void  LoadLibraries(const anaModes mode) {
     // gProof->ClearPackage("AOD");\r
     // gProof->ClearPackage("ANALYSIS");\r
     // gProof->ClearPackage("ANALYSISalice");\r
+    // gProof->ClearPackage("PHOSUtils");\r
+    // gProof->ClearPackage("EMCALUtils");\r
     // if( kDoJetTask ){\r
     //   gProof->ClearPackage("JETAN");\r
     //   gProof->ClearPackage("FASTJETAN");\r
     // }\r
+    // gProof->ClearPackage("PWG4PartCorrBase");\r
+    // gProof->ClearPackage("PWG4PartCorrDep");\r
     // gProof->ShowEnabledPackages();\r
 \r
     // Enable the STEERBase Package\r
@@ -415,6 +427,12 @@ void  LoadLibraries(const anaModes mode) {
     // Enable the Analysis Package\r
     gProof->UploadPackage("ANALYSISalice.par");\r
     gProof->EnablePackage("ANALYSISalice");\r
+    // Enable the PHOSUtils Package\r
+    gProof->UploadPackage("PHOSUtils.par");\r
+    gProof->EnablePackage("PHOSUtils");\r
+    // Enable the EMCALUtils Package\r
+    gProof->UploadPackage("EMCALUtils.par");\r
+    gProof->EnablePackage("EMCALUtils");\r
     if( kDoJetTask ){\r
       // Enable JETAN analysis\r
       gProof->UploadPackage("JETAN.par");\r
@@ -423,6 +441,12 @@ void  LoadLibraries(const anaModes mode) {
       gProof->UploadPackage("FASTJETAN.par");\r
       gProof->EnablePackage("FASTJETAN");\r
     }\r
+    // Enable the PWG4PartCorrBase Package\r
+    gProof->UploadPackage("PWG4PartCorrBase.par");\r
+    gProof->EnablePackage("PWG4PartCorrBase");\r
+    // Enable the PWG4PartCorrDep Package\r
+    gProof->UploadPackage("PWG4PartCorrDep.par");\r
+    gProof->EnablePackage("PWG4PartCorrDep");\r
 \r
     gProof->ShowEnabledPackages();\r
   }  \r
@@ -682,6 +706,9 @@ void ProcessEnvironment(){
   if (gSystem->Getenv("anaInputData"))\r
      kInputData = gSystem->Getenv("anaInputData");\r
 \r
+  if (gSystem->Getenv("anakMC"))\r
+     kMC = atoi(gSystem->Getenv("anakMC"));\r
+\r
   if (gSystem->Getenv("MODE"))\r
      mode = atoi(gSystem->Getenv("MODE"));\r
 \r
@@ -697,6 +724,6 @@ void ProcessEnvironment(){
   if (gSystem->Getenv("SEVENT"))\r
       sevent = atoi (gSystem->Getenv("SEVENT"));\r
        \r
-      printf("Process: Variables: kInputData %s, mode %d, config2 %s, config3 %s, sevent %d\n",kInputData.Data(),mode,sconfig2,sconfig3,sevent);\r
+      printf("Process: Variables: kInputData %s, kMC %d, mode %d, config2 %s, config3 %s, sevent %d\n",kInputData.Data(),kMC,mode,sconfig2,sconfig3,sevent);\r
 \r
 }\r