PWGLF/STRANGENESS/Correlations/macros/runGridLambdaOverK0sJets.C

   1
   2 // "ESDs/pass2/AOD086/*/AliAOD.root"
   3 // "/alice/data/2010/LHC10h"
   4
   5 class AliAnalysisAlien;
   6
   7 void runGridLambdaOverK0sJets(TString  runMode    = "full",
   8                               TString  alirootVer = "vAN-20140708",
   9                               TString  rootVer    = "v5-34-08-6",
  10                               TString  dataPath   = /*"AOD124/* /AliAOD.root"*/"ESDs/pass2/AOD115/*/AliAOD.root",
  11                               TString  dataDir    = /*"/alice/sim/2012/LHC12a17d_fix"*/"/alice/data/2011/LHC11h_2",
  12                               TString  workDir    = "test",
  13                               TString  name       = "LambdaOverK0sRatio",
  14                               TString  data       = "PbPb2011",
  15                               Float_t  minCen     = 0.,
  16                               Float_t  maxCen     = 40.,
  17                               Bool_t   sepInjec   = kFALSE,
  18                               Bool_t   isMC       = kFALSE,
  19                               Bool_t   usePID     = kFALSE,
  20                               Bool_t   doQA       = kTRUE,
  21                               Bool_t   useEtaCut  = kFALSE,
  22                               Int_t    run        = 169838/*137530*//*138624*/){
  23
  24
  25
  26   Printf("   \nThe parameters of the programm are : \n ");
  27   Printf(" \t Analysis mode:\t %s\n \t Centrality:\t %.1lf - %.1lf\n \t Use MC Data?:\t %s\n \t Use PID?:\t %s\n",
  28          "Grid",minCen,maxCen,
  29          (isMC) ? "Yes" : "No",
  30          (usePID) ? "Yes" : "No");
  31
  32   // _____________________________________________________ //
  33
  34   InitAndLoadLibs();
  35
  36   AliAnalysisManager *mgr = new AliAnalysisManager("Manager");
  37
  38   AliAnalysisGrid *alienHandler = CreateAlienHandler(runMode,alirootVer,rootVer,dataPath,dataDir,workDir,isMC,run);
  39   if (!alienHandler) return;
  40   mgr->SetGridHandler(alienHandler);
  41
  42   AliAODInputHandler* aodH = new AliAODInputHandler;
  43   mgr->SetInputEventHandler(aodH);
  44
  45   //PID
  46   gROOT->LoadMacro("$ALICE_ROOT/ANALYSIS/macros/AddTaskPIDResponse.C");
  47   AliAnalysisTask *pidTask = AddTaskPIDResponse(isMC);
  48   //AliAnalysisTask *pidTask = AddTaskPIDResponse(isMC,kTRUE);
  49   if(!pidTask) { printf("no PIDtask\n"); return; }
  50
  51   //Float_t checkIDTrig= kTRUE;
  52
  53   // My task
  54   gROOT->LoadMacro("AliAnalysisTaskLambdaOverK0sJets.cxx++g");
  55   //gSystem->Load("libPWGLFSTRANGENESS");
  56   gROOT->LoadMacro("AddTaskLambdaOverK0sJets.C");
  57   AliAnalysisTaskLambdaOverK0sJets *task = AddTaskLambdaOverK0sJets(name,data,minCen,maxCen,sepInjec,isMC,doQA,useEtaCut);
  58    // _____________________________________________________ //
  59
  60    if (!mgr->InitAnalysis()) return;
  61    mgr->PrintStatus();
  62
  63    mgr->StartAnalysis("grid");
  64 }
  65
  66 // ______________________________________________________________
  67
  68 void InitAndLoadLibs() {
  69
  70     gSystem->Load("libCore");
  71     gSystem->Load("libTree");
  72     gSystem->Load("libGeom");
  73     gSystem->Load("libVMC");
  74     gSystem->Load("libPhysics");
  75     gSystem->Load("libMinuit");
  76     gSystem->Load("libProof");
  77     gSystem->Load("libGui");
  78     gSystem->Load("libXMLParser");
  79     gSystem->Load("libProofPlayer");
  80     gSystem->Load("libXMLIO");
  81
  82     gSystem->Load("libSTEERBase");
  83     gSystem->Load("libESD");
  84     gSystem->Load("libAOD");
  85     gSystem->Load("libCDB");
  86     gSystem->Load("libANALYSIS");
  87     gSystem->Load("libANALYSISalice");
  88     gSystem->Load("libCORRFW");
  89     gSystem->Load("libJETAN");
  90     gSystem->Load("libRAWDatabase");
  91     gSystem->Load("libSTEER");
  92     gSystem->Load("libCORRFW");
  93
  94
  95     gSystem->AddIncludePath("-I$ALICE_ROOT/include");
  96 }
  97
  98 // ___________________________________________________________________ //
  99
 100 // **** It is change 'AliAnalysisGrid' by 'AliAnalysisAlien'
 101 AliAnalysisAlien* CreateAlienHandler(TString runMode,TString alirootVer,
 102                                      TString rootVer,TString dataPath,
 103                                      TString dataDir,TString workDir,
 104                                      Bool_t isMC,Int_t kRun) {
 105
 106   AliAnalysisAlien *plugin = new AliAnalysisAlien();
 107   plugin->SetCheckCopy(kFALSE);
 108
 109   plugin->SetRunMode(runMode);
 110   // Set versions of used packages
 111   plugin->SetAPIVersion("V1.1x");
 112   plugin->SetROOTVersion(rootVer);
 113   plugin->SetAliROOTVersion(alirootVer);
 114
 115   /////////////////////////////////////////////////////////////////
 116
 117   // Declare input data to be processed.
 118   // Method 1: Create automatically XML collections using alien 'find' command.
 119   plugin->SetGridDataDir(dataDir);
 120
 121   if (!isMC)
 122     plugin->SetRunPrefix("000");
 123   plugin->SetDataPattern(dataPath);
 124
 125   plugin->AddRunNumber(kRun);
 126
 127   // Method 2: Declare existing data files (raw collections, xml collections, root file)
 128   const char working_dir[250];
 129
 130   sprintf(working_dir, "%s/%d",workDir.Data(),kRun);
 131   TString path = TString(working_dir);
 132   plugin->SetGridWorkingDir(path);
 133
 134   //plugin->SetGridWorkingDir(workDir);
 135   // Declare alien output directory. Relative to working directory.
 136   plugin->SetGridOutputDir("output"); // In this case will be $HOME/work/output
 137   // Declare the analysis source files names separated by blancs. To be compiled runtime
 138   plugin->SetAnalysisSource("AliAnalysisTaskLambdaOverK0sJets.cxx");
 139   //plugin->SetAdditionalLibs("AliAnalysisTaskMultiplicity.h AliAnalysisTaskMultiplicity.cxx");
 140   plugin->SetAdditionalLibs("AliAnalysisTaskLambdaOverK0sJets.cxx AliAnalysisTaskLambdaOverK0sJets.h");
 141   // Optionally set a name for the generated analysis macro (default MyAnalysis.C)
 142   plugin->SetAnalysisMacro("AnalysisProduction.C");
 143   //plugin->SetAnalysisMacro("mytask.C");
 144   // Optionally set maximum number of input files/subjob (default 100, put 0 to ignore)
 145   plugin->SetSplitMaxInputFileNumber(20);
 146   // Optionally set number of failed jobs that will trigger killing waiting sub-jobs.
 147   plugin->SetMaxInitFailed(12);
 148   // Optionally resubmit threshold.
 149   plugin->SetMasterResubmitThreshold(90);
 150   // Optionally set time to live (default 30000 sec)
 151   plugin->SetTTL(30000);
 152   // Optionally set input format (default xml-single)
 153   plugin->SetInputFormat("xml-single");
 154   // Optionally modify the name of the generated JDL (default analysis.jdl)
 155   plugin->SetJDLName("TaskProduction.jdl");
 156   //plugin->SetJDLName("mytask.jdl");
 157   plugin->SetMergeViaJDL(kTRUE);
 158   // Optionally modify job price (default 1)
 159   plugin->SetPrice(1);
 160   // Optionally modify split mode (default 'se')
 161   plugin->SetSplitMode("se");
 162   return plugin;
 163
 164 }