]> git.uio.no Git - u/mrichter/AliRoot.git/blobdiff - PWGJE/FlavourJetTasks/macros/AnalysisTrainCorrJetsLocal.C
TENDER becomes Tender
[u/mrichter/AliRoot.git] / PWGJE / FlavourJetTasks / macros / AnalysisTrainCorrJetsLocal.C
index 98540288cb10eb35740f2c77244213b0712a20a0..5e2a84baa5f78e409ffcd1b3787812128b3d22bc 100644 (file)
@@ -1,5 +1,3 @@
-// $Id$
-
 // runEMCalJetAnalysis.C
 // =====================
 // This macro can be used to run a jet analysis within the EMCal Jet Framework.
@@ -32,31 +30,34 @@ AliAnalysisGrid* CreateAlienHandler(const char* uniqueName, const char* gridDir,
 void AnalysisTrainCorrJetsLocal (
          const char*    dataType            = "AOD",                       // set the analysis type, AOD, ESD or sESD
          Bool_t         useGrid             = kTRUE,                      // local or grid
-        TString localfilename = "file_aodlhc10d.txt",
+        TString localfilename = "/data/Work/jets/testEMCalJetFramework/ptJdistrAug2nd/setMCtasks/map/inheritDev/files_LHC10f7a.txt",
          const char*    gridMode            = "test",                      // set the grid run mode (can be "full", "test", "offline", "submit" or "terminate")
-         const char*    pattern             = "*ESDs/pass2/AOD099/*AliAOD.root", //"*/*/AliAOD.root" "*ESDs/pass1/AOD106/*AliAOD.root",    // file pattern (here one can specify subdirs like passX etc.) (used on grid)
-         const char*    gridDir             = "/alice/data/2010/LHC10d",   // /alice/data/2011/LHC11d /alice/sim/2012/LHC12f2b   dir on alien, where the files live (used on grid)
-         const char*    runNumbers          = /*126437*/" 126432 126425 126424 126422 126409 126408 126407 126406 126405 126404 126403 126359 126352 126351 126350 126285 126284 126283 126168 126167 126160 126158 126097 126090 126088 126082 126081 126078 126073 126008 126007 126004 125855 125851 125850 125849 125848 125847 125844 125843 125842 125633 125632 125630 125296 125134 125101 125100 125097 125085 125023 124751 122375 122374",             // considered run numbers (used on grid) /*LHC12g 188359 188362, LHC11a 146860 146859*/ /*LHC12f2b 158285 159582 */
+         const char*    pattern             = "*AOD136/*AliAOD.root", //"*/*/AliAOD.root" "*ESDs/pass1/AOD106/*AliAOD.root",    // file pattern (here one can specify subdirs like passX etc.) (used on grid)
+         const char*    gridDir             = "/alice/sim/LHC10f7a",   // /alice/data/2011/LHC11d /alice/sim/2012/LHC12f2b   dir on alien, where the files live (used on grid)
+         const char*    runNumbers          = /*114931 */" 115186 115193 115393 115401 115414 116102 116288 116402 116403 116562 116571 116574 116643 116645 117048 117050 117052 117053 117054 117059 117060 117063 117065 117077 117086 117092 117099 117109 117112 117116 117220 117222 119159 119161 119163 119841 119842 119844 119845 119846 119849 119853 119856 119859 119862 120067 120069 120072 120073 120076 120079 120244 120503 120504 120505 120616 120617 120671 120741 120750 120758 120820 120821 120822 120823 120824 120825 120829 122374 122375 124187 124191 124355 124358 124362 124367 124378 124380 124381 124604 124605 124606 124607 124608 124702 124746 124750 124751 125023 125085 125097 125100 125101 125133 125134 125139 125140 125156 125186 125296 125628 125630 125632 125633 125842 125843 125844 125847 125848 125849 125850 125851 125855 126004 126007 126008 126073 126078 126081 126082 126088 126090 126097 126158 126160 126167 126168 126283 126284 126285 126350 126351 126352 126359 126403 126404 126405 126406 126407 126408 126409 126422 126424 126425 126432 126437 127719 127724 127729 127730 127814 127819 127930 127940 128263 128778 128913 129536 129599 129639 129641 129654 129659 129666 129667 129723 129725 129726 129729 129735 129736 129738 129742 129744 129959 129960 129961 129962 129966 129983 130149 130151 130157 130158 130172 130178 130179 130342 130343 130354 130356 130358 130360 130375 130479 130480 130481 130517 130519 130520 130524 130526 130601 130608 130620 130621 130623 130628 130696 130704 130793 130795 130798 130799 130834 130840 130842 130844 130847 130848",             // considered run numbers (used on grid) /*LHC12g 188359 188362, LHC11a 146860 146859*/ /*LHC12f2b 158285 159582 */
         const Int_t nrunspermaster= 100,
-         UInt_t         numLocalFiles       = 3,                          // number of files analyzed locally  
-         const char*    runPeriod           = "LHC10d",                    // set the run period (used on grid)
-         const char*    uniqueName          = "DJetNewCode",     // sets base string for the name of the task on the grid
+         UInt_t         numLocalFiles       = 10,                          // number of files analyzed locally  
+         const char*    runPeriod           = "LHC10f7a",                    // set the run period (used on grid)
+         const char*    uniqueName          = "DJetNewCodeMCSandBchJ",     // sets base string for the name of the task on the grid
          UInt_t         pSel                = AliVEvent::kAny,             // used event selection for every task except for the analysis tasks
          Bool_t         useTender           = kFALSE,                      // trigger, if tender task should be used
-         Bool_t         isMC                = kFALSE,                      // trigger, if MC handler should be used
-
+         Bool_t         isMC                = kTRUE,                      // trigger, if MC handler should be used
+        Bool_t         isReco                = kTRUE,
+        Bool_t         isMap                = kTRUE,
          // Here you have to specify additional code files you want to use but that are not in aliroot
-         const char*    addCXXs             = "AliAnalysisTaskRecoJetCorrelations.cxx AliAnalysisTaskSEDmesonsForJetCorrelations.cxx",
-         const char*    addHs               = "AliAnalysisTaskRecoJetCorrelations.h AliAnalysisTaskSEDmesonsForJetCorrelations.h",
+         const char*    addCXXs             = "",
+         const char*    addHs               = "",
 
          // These two settings depend on the dataset and your quotas on the AliEN services
          Int_t          maxFilesPerWorker   = 20,
          Int_t          workerTTL           = 86000,
-        Int_t          nfiletestmode       = 1
+        Int_t          nfiletestmode       = 3
          )
 {
 
   // Some pre-settings and constants
+  TStopwatch watch;
+  watch.Start();
 
   enum AlgoType {kKT, kANTIKT};
   enum JetType  {kFULLJETS, kCHARGEDJETS, kNEUTRALJETS};
@@ -113,21 +114,24 @@ void AnalysisTrainCorrJetsLocal (
   if(!useGrid)
     cout << "Using " << localFiles.Data() << " as input file list.\n";
 
-  gROOT->LoadMacro("AddTasksCorrJets.C");
+  gROOT->LoadMacro("$ALICE_ROOT/PWGJE/FlavourJetTasks/macros/AddTasksFlavourJet.C");
   //List of arguments:
 
-// const Int_t iCandtype=1 /*0 = D0, 1=Dstar...*/,
-// const TString sCutFile = "cutsHF/D0toKpiCutsppRecVtxNoPileupRejNoEMCAL.root",
-// const Double_t dJetPtCut   = 1.,
-// const Double_t dJetAreaCut = 0.,
-// const Int_t iAccCut = 1,
-// const TString sRunPeriod = "LHC10b",
-// const Int_t    uBeamType   = 0,
-// const UInt_t uTriggerMask = AliVEvent::kMB, /*for jets and phys sel; the D mesons trigger is defined in the cut object*/
-// const Bool_t bIsMC=kFALSE,
-// TString sText=""/*completes the name of the candidate task lists*/
-
-     AddTasksCorrJets(1,"/data/Work/jets/testEMCalJetFramework/CutFilesMB/DStartoKpipiCuts_new.root",10.,0.,1,runPeriod,0,pSel,isMC,"");
+  // const Int_t iCandType = 1 /*0 = D0, 1=Dstar...*/,
+  //  const TString sCutFile = "cutsHF/D0toKpiCutsppRecVtxNoPileupRejNoEMCAL.root",
+  //  const Double_t dJetPtCut   = 1.,
+  //  const Double_t dJetAreaCut = 0.,
+  //  const char *acctype = "TPC",
+  //  const TString sRunPeriod = "LHC10b",
+  //  const Int_t    uBeamType = 0,
+  //  const UInt_t uTriggerMask = AliVEvent::kMB, /*for jets; the D mesons trigger is defined in the cut object*/
+  //  const Bool_t bIsMC = kFALSE,
+  //  const Bool_t bIsReco = kFALSE,
+  //  const Bool_t bIsMap = kFALSE,
+  //  TString sText=""/*completes the name of the candidate task lists*/
+
+  //AddTasksFlavourJet(1,"/data/Work/jets/testEMCalJetFramework/CutFilesMB/DStartoKpipiCuts_new.root",10.,0.,"TPC",runPeriod,0,pSel,isMC,isReco,isMap,"");
+  AddTasksFlavourJet(1,"/data/Work/jets/testEMCalJetFramework/CutFilesMB/DStartoKpipiCuts.root",5.,0.,"TPC",runPeriod,0,pSel,isMC,isReco,isMap,"");
 
   // Set the physics selection for all given tasks
   TObjArray *toptasks = mgr->GetTasks();
@@ -184,6 +188,8 @@ void AnalysisTrainCorrJetsLocal (
     mgr->SetDebugLevel(10);
     mgr->StartAnalysis("local", chain);
   }
+  watch.Stop();
+  watch.Print();
 }
 
 //______________________________________________________________________________
@@ -210,56 +216,65 @@ void LoadLibs()
   gSystem->Load("libEVGEN");
   gSystem->Load("libANALYSISalice");
   gSystem->Load("libCORRFW");
-  gSystem->Load("libTOFbase");
+  //gSystem->Load("libTOFbase");
   //gSystem->Load("libTOFrec");
-  gSystem->Load("libRAWDatabase.so");
-  gSystem->Load("libRAWDatarec.so");
-  gSystem->Load("libTPCbase.so");
-  gSystem->Load("libTPCrec.so");
-  gSystem->Load("libITSbase.so");
-  gSystem->Load("libITSrec.so");
-  gSystem->Load("libTRDbase.so");
-  gSystem->Load("libTENDER.so");
-  gSystem->Load("libSTAT.so");
-  gSystem->Load("libTRDrec.so");
-  gSystem->Load("libHMPIDbase.so");
-  gSystem->Load("libPWGPP.so");
+  //gSystem->Load("libRAWDatabase.so");
+  //gSystem->Load("libRAWDatarec.so");
+  //gSystem->Load("libTPCbase.so");
+  //gSystem->Load("libTPCrec.so");
+  //gSystem->Load("libITSbase.so");
+  //gSystem->Load("libITSrec.so");
+  //gSystem->Load("libTRDbase.so");
+  //gSystem->Load("libTender.so");
+  //gSystem->Load("libSTAT.so");
+  //gSystem->Load("libTRDrec.so");
+  //gSystem->Load("libHMPIDbase.so");
+  gSystem->Load("libPWGTools.so");
+  //gSystem->Load("libPWGPP.so");
   gSystem->Load("libPWGHFbase");
-  gSystem->Load("libPWGDQdielectron");
-  gSystem->Load("libPWGHFhfe");
+  //gSystem->Load("libPWGDQdielectron");
+  //gSystem->Load("libPWGHFhfe");
   gSystem->Load("libPWGflowBase.so");
   gSystem->Load("libPWGflowTasks.so");
   gSystem->Load("libPWGHFvertexingHF");
   gSystem->Load("libEMCALUtils");
-  gSystem->Load("libPHOSUtils");
-  gSystem->Load("libPWGCaloTrackCorrBase");
+  //gSystem->Load("libPHOSUtils");
+  //gSystem->Load("libPWGCaloTrackCorrBase");
   gSystem->Load("libEMCALraw");
   gSystem->Load("libEMCALbase");
-  gSystem->Load("libEMCALrec");
-  gSystem->Load("libTRDbase");
-  gSystem->Load("libVZERObase");
-  gSystem->Load("libVZEROrec");
-  gSystem->Load("libTENDER");
-  gSystem->Load("libTENDERSupplies");
+  //gSystem->Load("libEMCALrec");
+  //gSystem->Load("libTRDbase");
+  //gSystem->Load("libVZERObase");
+  //gSystem->Load("libVZEROrec");
+  gSystem->Load("libTender");
+  gSystem->Load("libTenderSupplies");
+  gSystem->Load("libESDfilter");
   gSystem->Load("libPWGEMCAL");
   gSystem->Load("libPWGGAEMCALTasks");
-  gSystem->Load("libPWGTools");
-  gSystem->Load("libPWGCFCorrelationsBase");
-  gSystem->Load("libPWGCFCorrelationsDPhi");
+  //gSystem->Load("libPWGCFCorrelationsBase");
+  //gSystem->Load("libPWGCFCorrelationsDPhi");
 
-  //load CGAL, Fastjet and SISCone
+  //load CGAL, Fastjet
   //gSystem->Load("libCGAL");
   gSystem->Load("libfastjet");
-  gSystem->Load("libSISConePlugin");
-
+  gSystem->Load("libfastjet");
+  gSystem->Load("libsiscone");
+  gSystem->Load("libsiscone_spherical");
+  gSystem->Load("libfastjetplugins");
+  gSystem->Load("libfastjettools");
+  gSystem->Load("libfastjetcontribfragile");
   gSystem->Load("libJETAN");
   gSystem->Load("libFASTJETAN");
+  gSystem->Load("libPWGGAEMCALTasks");
+
   gSystem->Load("libPWGJEEMCALJetTasks");
+  gSystem->Load("libPWGJEFlavourJetTasks");
 
 
   // include paths
   gSystem->AddIncludePath("-Wno-deprecated");
-  gSystem->AddIncludePath("-I$ALICE_ROOT -I$ALICE_ROOT/include -I$ALICE_ROOT/EMCAL");
+  gSystem->AddIncludePath("-I$ALICE_ROOT -I$ALICE_ROOT/include -I$ALICE_ROOT/EMCAL -I$ALICE_ROOT/PWG/EMCAL -I$ALICE_ROOT/PWGJE/EMCALJetTasks");
+  gSystem->AddIncludePath(" -I$ALICE_ROOT/PWGHF/ -I$ALICE_ROOT/PWGHF/base -I$ALICE_ROOT/PWGHF/vertexingHF -I$ALICE_ROOT/PWGJE/FlavourJetTasks");
   gSystem->AddIncludePath("-I$ALICE_ROOT/PWGDQ/dielectron -I$ALICE_ROOT/PWGHF/hfe");
   gSystem->AddIncludePath("-I$ALICE_ROOT/JETAN -I$ALICE_ROOT/JETAN/fastjet");
 }
@@ -281,7 +296,7 @@ AliAnalysisGrid* CreateAlienHandler(const char* uniqueName, const char* gridDir,
   }else tmpName +="_20130412_122423";
   */
   TString tmpAdditionalLibs("");
-  tmpAdditionalLibs = Form("libTree.so libVMC.so libGeom.so libGui.so libXMLParser.so libMinuit.so libMinuit2.so libProof.so libPhysics.so libSTEERBase.so libESD.so libAOD.so libOADB.so libANALYSIS.so libCDB.so libRAWDatabase.so libSTEER.so libANALYSISalice.so libCORRFW.so libTOFbase.so libRAWDatabase.so libRAWDatarec.so libTPCbase.so libTPCrec.so libITSbase.so libITSrec.so libTRDbase.so libTENDER.so libSTAT.so libTRDrec.so libHMPIDbase.so libPWGPP.so libPWGHFbase.so libPWGDQdielectron.so libPWGHFhfe.so libPWGflowBase.so libPWGflowTasks.so libPWGHFvertexingHF.so libEMCALUtils.so libPHOSUtils.so libPWGCaloTrackCorrBase.so libEMCALraw.so libEMCALbase.so libEMCALrec.so libTRDbase.so libVZERObase.so libVZEROrec.so libTENDER.so libTENDERSupplies.so libPWGEMCAL.so libPWGGAEMCALTasks.so libPWGTools.so libPWGCFCorrelationsBase.so libPWGCFCorrelationsDPhi.so  libCGAL.so libJETAN.so libfastjet.so libSISConePlugin.so libFASTJETAN.so libPWGJE.so libPWGmuon.so libPWGJEEMCALJetTasks.so %s %s",additionalCode.Data(),additionalHeaders.Data());
+  tmpAdditionalLibs = Form("libTree.so libVMC.so libGeom.so libGui.so libXMLParser.so libMinuit.so libMinuit2.so libProof.so libPhysics.so libSTEERBase.so libESD.so libESDfilter.so libAOD.so libOADB.so libANALYSIS.so libCDB.so libRAWDatabase.so libSTEER.so libANALYSISalice.so libCORRFW.so libPWGTools.so libPWGHFbase.so libPWGflowBase.so libPWGflowTasks.so libPWGHFvertexingHF.so libEMCALUtils.so libEMCALraw.so libEMCALbase.so libTender.so libPWGEMCAL.so libPWGGAEMCALTasks.so libPWGTools.so libCGAL.so libfastjet.so libsiscone.so libsiscone_spherical.so libfastjetplugins.so libfastjettools.so libfastjetcontribfragile.so libJETAN.so libFASTJETAN.so  libPWGJE.so libPWGJEEMCALJetTasks.so libPWGJEFlavourJetTasks.so %s %s",additionalCode.Data(),additionalHeaders.Data());
 
 
   TString macroName("");
@@ -297,8 +312,8 @@ AliAnalysisGrid* CreateAlienHandler(const char* uniqueName, const char* gridDir,
 
   // Here you can set the (Ali)ROOT version you want to use
   plugin->SetAPIVersion("V1.1x");
-  plugin->SetROOTVersion("v5-34-05");
-  plugin->SetAliROOTVersion("v5-04-38-AN");
+  plugin->SetROOTVersion(":v5-34-08-6");
+  plugin->SetAliROOTVersion("vAN-20141005");
 
   plugin->SetGridDataDir(gridDir); // e.g. "/alice/sim/LHC10a6"
   plugin->SetDataPattern(pattern); //dir structure in run directory
@@ -314,10 +329,10 @@ AliAnalysisGrid* CreateAlienHandler(const char* uniqueName, const char* gridDir,
 
   plugin->SetAnalysisSource(additionalCode.Data());
   plugin->SetAdditionalLibs(tmpAdditionalLibs.Data());
-  plugin->AddIncludePath("-I. -I$ROOTSYS/include -I$ALICE_ROOT -I$ALICE_ROOT/include -I$ALICE_ROOT/STEER/STEER -I$ALICE_ROOT/STEER/STEERBase -I$ALICE_ROOT/STEER/ESD -I$ALICE_ROOT/STEER/AOD -I$ALICE_ROOT/ANALYSIS  -I$ALICE_ROOT/OADB -I$ALICE_ROOT/PWGHF -I$ALICE_ROOT/PWGHF/base -I$ALICE_ROOT/PWGHF/vertexingHF -I$ALICE_ROOT/PWG/FLOW/Base -I$ALICE_ROOT/PWG/FLOW/Tasks  -I$ALICE_ROOT/PWGJE  -I$ALICE_ROOT/JETAN -I$ALICE_ROOT/PWGJE/EMCALJetTasks -g");
-  plugin->AddExternalPackage("boost::v1_43_0");
-  plugin->AddExternalPackage("cgal::v3.6");
-  plugin->AddExternalPackage("fastjet::v2.4.2");
+  plugin->AddIncludePath("-I. -I$ROOTSYS/include -I$ALICE_ROOT -I$ALICE_ROOT/include -I$ALICE_ROOT/STEER/STEER -I$ALICE_ROOT/STEER/STEERBase -I$ALICE_ROOT/STEER/ESD -I$ALICE_ROOT/STEER/AOD -I$ALICE_ROOT/ANALYSIS  -I$ALICE_ROOT/OADB -I$ALICE_ROOT/PWGHF -I$ALICE_ROOT/PWGHF/base -I$ALICE_ROOT/PWGHF/vertexingHF -I$ALICE_ROOT/PWG/FLOW/Base -I$ALICE_ROOT/PWG/FLOW/Tasks  -I$ALICE_ROOT/PWGJE  -I$ALICE_ROOT/JETAN -I$ALICE_ROOT/PWGJE/EMCALJetTasks -I$ALICE_ROOT/PWGJE/FlavourJetTasks -g");
+  plugin->AddExternalPackage("boost::v1_53_0");
+  plugin->AddExternalPackage("cgal::v4.4");
+  plugin->AddExternalPackage("fastjet::v3.0.6_1.012");
 
   plugin->SetDefaultOutputs(kTRUE);
   // merging via jdl