Faster calculation of invariant mass (Yifei)

[u/mrichter/AliRoot.git] / PWG3 / vertexingHF / RunAnalysisAODVertexingHF.C

diff --git a/PWG3/vertexingHF/RunAnalysisAODVertexingHF.C b/PWG3/vertexingHF/RunAnalysisAODVertexingHF.C
index fab76b5bd9390ebd8feed45eadec50d640c10acc..547e831cc31a314d9f56d18026929f91efea3969 100644 (file)
--- a/PWG3/vertexingHF/RunAnalysisAODVertexingHF.C
+++ b/PWG3/vertexingHF/RunAnalysisAODVertexingHF.C
@@ -27,7 +27,7 @@ void RunAnalysisAODVertexingHF()
   TString pluginmode="test";\r
   Bool_t saveProofToAlien=kFALSE;\r
   TString proofOutdir = "";\r
-  TString loadMacroPath="$ALICE_ROOT/PWG3/vertexingHF/";\r
+  TString loadMacroPath="$ALICE_ROOT/PWG3/vertexingHF/macros/";\r
   //TString loadMacroPath="./"; // this is normally needed for CAF\r
   //\r
 \r
@@ -164,7 +164,7 @@ void RunAnalysisAODVertexingHF()
 \r
   // Input\r
   AliAODInputHandler *inputHandler = new AliAODInputHandler();\r
-  if(analysisMode=="proof") {\r
+  if(analysisMode=="proof" ) {\r
     inputHandler->AddFriend("AliAOD.VertexingHF.root");\r
     if(saveProofToAlien) mgr->SetSpecialOutputLocation(proofOutdir);\r
   }\r
@@ -176,6 +176,14 @@ void RunAnalysisAODVertexingHF()
   // Analysis tasks (wagons of the train)   \r
   //\r
   TString taskName;\r
+  ////// ADD THE FULL D2H TRAIN\r
+  taskName="AddD2HTrain.C"; taskName.Prepend(loadMacroPath.Data());\r
+  gROOT->LoadMacro(taskName.Data());\r
+  Bool_t readMC=kFALSE;\r
+  AddD2HTrain(readMC);//,1,0,0,0,0,0,0,0,0,0,0);\r
+\r
+  ////// OR ADD INDIVIDUAL TASKS\r
+  /*\r
   \r
   //taskName="AddTaskCompareHF.C"; taskName.Prepend(loadMacroPath.Data());\r
   //gROOT->LoadMacro(taskName.Data());\r
@@ -190,6 +198,10 @@ void RunAnalysisAODVertexingHF()
   gROOT->LoadMacro(taskName.Data());\r
   AliAnalysisTaskSEDplus *dplusTask = AddTaskDplus();\r
   \r
+  taskName="AddTaskDs.C"; taskName.Prepend(loadMacroPath.Data());\r
+  gROOT->LoadMacro(taskName.Data());\r
+  AliAnalysisTaskSEDs *dsTask = AddTaskDs();\r
+\r
   //taskName="AddTaskSelectHF.C"; taskName.Prepend(loadMacroPath.Data());\r
   //gROOT->LoadMacro(taskName.Data());\r
   //AliAnalysisTaskSESelectHF *seleTask = AddTaskSelectHF();\r
@@ -213,8 +225,8 @@ void RunAnalysisAODVertexingHF()
   taskName="AddTaskCharmFraction.C"; taskName.Prepend(loadMacroPath.Data());\r
   gROOT->LoadMacro(taskName.Data());\r
   Int_t switchMC[5]={1,1,1,1,1};\r
-  AliAnalysisTaskSECharmFraction *cFractTask = AddTaskCharmFraction("d0D0.root",switchMC);\r
-  \r
+  AliAnalysisTaskSECharmFraction *cFractTask = AddTaskCharmFraction("d0D0",switchMC);\r
+  */\r
   // attach a private task (not committed)\r
   // (the files MyTask.h MyTask.cxx AddMyTask.C have to be declared in plugin\r
   // configuration, see below)\r
@@ -267,7 +279,7 @@ AliAnalysisGrid* CreateAlienHandler(TString pluginmode="test",Bool_t useParFiles
    // Set the run mode (can be "full", "test", "offline", "submit" or "terminate")\r
    plugin->SetRunMode(pluginmode.Data());\r
    plugin->SetUser("dainesea");\r
-   plugin->SetNtestFiles(2);\r
+   plugin->SetNtestFiles(1);\r
    // Set versions of used packages\r
    plugin->SetAPIVersion("V2.4");\r
    plugin->SetROOTVersion("v5-24-00");\r
@@ -275,12 +287,12 @@ AliAnalysisGrid* CreateAlienHandler(TString pluginmode="test",Bool_t useParFiles
    // Declare input data to be processed.\r
    // Method 1: Create automatically XML collections using alien 'find' command.\r
    // Define production directory LFN\r
-   //plugin->SetGridDataDir("/alice/cern.ch/user/r/rbala/newtrain/out_lhc08x/");\r
-   //plugin->SetGridDataDir("/alice/cern.ch/user/m/mgheata/analysisESD/output_train_default_28May2009_09h33/");\r
-   plugin->SetGridDataDir("/alice/sim/PDC_09/LHC09a5/AOD1/");\r
+   plugin->SetGridDataDir("/alice/cern.ch/user/r/rbala/data_pass4_good_runCINT1B_8thfeb/");\r
+   //plugin->SetGridDataDir("/alice/sim/PDC_09/LHC09a4/AOD3/");\r
    // Set data search pattern\r
    plugin->SetDataPattern("AliAOD.root");\r
-   plugin->SetFriendChainName("AliAOD.VertexingHF.root");\r
+   plugin->SetFriendChainName("./AliAOD.VertexingHF.root");\r
+   //plugin->SetFriendChainName("deltas/AliAOD.VertexingHF.root");\r
    // ...then add run numbers to be considered\r
    //plugin->AddRunNumber(529007);\r
    //  or\r
@@ -300,7 +312,7 @@ AliAnalysisGrid* CreateAlienHandler(TString pluginmode="test",Bool_t useParFiles
    //plugin->SetAnalysisSource("AliDStarJets.cxx");\r
    // Declare all libraries (other than the default ones for the framework. These will be\r
    // loaded by the generated analysis macro. Add all extra files (task .cxx/.h) here.\r
-   plugin->SetAdditionalLibs("libPWG3vertexingHF.so libPWG3base.so libPWG3muon.so libPWG4PartCorrBase.so libPWG4PartCorrDep.so MakeAODInputChain.C"/* AliDStarJets.cxx AliDStarJets.h"*/);\r
+   plugin->SetAdditionalLibs("libPWG3vertexingHF.so libPWG3base.so libPWG3muon.so");\r
    // use par files\r
    if(useParFiles) {\r
      plugin->EnablePackage("STEERBase.par");\r
@@ -312,8 +324,6 @@ AliAnalysisGrid* CreateAlienHandler(TString pluginmode="test",Bool_t useParFiles
      plugin->EnablePackage("PWG3base.par");\r
      plugin->EnablePackage("PWG3vertexingHF.par");\r
      plugin->EnablePackage("PWG3muon.par");\r
-     plugin->EnablePackage("PWG4PartCorrBase.par");\r
-     plugin->EnablePackage("PWG4PartCorrDep.par");\r
    }\r
    plugin->AddIncludePath("-I. -I$ROOTSYS/include -I$ALICE_ROOT -I$ALICE_ROOT/include -I$ALICE_ROOT/ITS -I$ALICE_ROOT/TPC -I$ALICE_ROOT/CONTAINERS -I$ALICE_ROOT/STEER -I$ALICE_ROOT/TRD -I$ALICE_ROOT/macros -I$ALICE_ROOT/ANALYSIS -I$ALICE_ROOT/PWG3 -I$ALICE_ROOT/PWG3/vertexingHF -g");\r
    // Declare the output file names separated by blancs.\r