PWG2/FEMTOSCOPY/macros/Batch/runBatch.C

   1 void runBatch(const char *chainlistfile) {
   2   TStopwatch timer;
   3   timer.Start();
   4
   5   printf("*** Connect to AliEn ***\n");
   6   TGrid::Connect("alien://");
   7   gSystem->Load("libProofPlayer.so");
   8
   9   //____________________________________________________//
  10   //_____________Setting up STEERBase.par_______________//
  11   //____________________________________________________//
  12   setupPar("STEERBase");
  13   gSystem->Load("libSTEERBase.so");
  14
  15   //____________________________________________________//
  16   //_____________Setting up ESD.par_____________________//
  17   //____________________________________________________//
  18   setupPar("ESD");
  19   gSystem->Load("libVMC.so");
  20   gSystem->Load("libESD.so");
  21
  22   //____________________________________________________//
  23   //_____________Setting up AOD.par_____________________//
  24   //____________________________________________________//
  25   setupPar("AOD");
  26   gSystem->Load("libAOD.so");
  27
  28   //_________________________________________________________//
  29   //_____________Setting up ANALYSIS.par_____________________//
  30   //_________________________________________________________//
  31   setupPar("ANALYSIS");
  32   gSystem->Load("libANALYSIS.so");
  33
  34   //_________________________________________________________//
  35   //_____________Setting up ANALYSISalice.par________________//
  36   //_________________________________________________________//
  37   setupPar("ANALYSISalice");
  38   gSystem->Load("libANALYSISalice.so");
  39
  40   //____________________________________________________//
  41   //_____________Setting up PWG2AOD.par_________________//
  42   //____________________________________________________//
  43   setupPar("PWG2AOD");
  44   gSystem->Load("libPWG2AOD.so");
  45
  46   //____________________________________________________//
  47   //_____________Setting up PWG2femtoscopy.par__________//
  48   //____________________________________________________//
  49   setupPar("PWG2femtoscopy");
  50   gSystem->Load("libPWG2femtoscopy.so");
  51
  52   //____________________________________________________//
  53   //_____________Setting up PWG2femtoscopyUser.par______//
  54   //____________________________________________________//
  55   setupPar("PWG2femtoscopyUser");
  56   gSystem->Load("libPWG2femtoscopyUser.so");
  57
  58   //ANALYSIS PART
  59   gSystem->SetIncludePath("-I$ROOTSYS/include  -I\"/usr/local/CERN/root/include\" -I./PWG2femtoscopy/FEMTOSCOPY/AliFemto -I./PWG2femtoscopyUser/FEMTOSCOPY/AliFemtoUser -I./ESD -I./AOD -I./ANALYSIS -I./PWG2AOD/AOD");
  60   gROOT->LoadMacro("ConfigFemtoAnalysis.C++");
  61   gROOT->LoadMacro("AliAnalysisTaskFemto.cxx+");
  62   //  const char *collectionfile = "wn.xml";
  63
  64   //____________________________________________//
  65   //Usage of event tags
  66   AliTagAnalysis *analysis = new AliTagAnalysis();
  67   TChain *chain = 0x0;
  68   chain = analysis->GetChainFromCollection(chainlistfile,"esdTree");
  69
  70   //____________________________________________//
  71   // Make the analysis manager
  72   AliAnalysisManager *mgr = new AliAnalysisManager("TestManager");
  73   AliESDInputHandler* esdH = new AliESDInputHandler;
  74   esdH->SetInactiveBranches("FMD CaloCluster");
  75   mgr->SetInputEventHandler(esdH);
  76   //____________________________________________//
  77   // 1st Pt task
  78   AliAnalysisTaskFemto *task1 = new AliAnalysisTaskFemto("TaskFemto");
  79
  80   mgr->AddTask(task1);
  81
  82   // Create containers for input/output
  83   AliAnalysisDataContainer *cinput1 = mgr->GetCommonInputContainer();
  84   AliAnalysisDataContainer *coutput1 = mgr->CreateContainer("clist1", TList::Class(),AliAnalysisManager::kOutputContainer,"Femto.ESD.root");
  85
  86   //____________________________________________//
  87 //  cinput1->SetData(chain);
  88   mgr->ConnectInput(task1,0,cinput1);
  89   mgr->ConnectOutput(task1,0,coutput1);
  90
  91   if (!mgr->InitAnalysis()) return;
  92   mgr->PrintStatus();
  93   mgr->StartAnalysis("local",chain);
  94
  95   timer.Stop();
  96   timer.Print();
  97 }
  98
  99 Int_t setupPar(const char* pararchivename) {
 100   ///////////////////
 101   // Setup PAR File//
 102   ///////////////////
 103   if (pararchivename) {
 104     char processline[1024];
 105     sprintf(processline,".! tar xvzf %s.par",pararchivename);
 106     gROOT->ProcessLine(processline);
 107     const char* ocwd = gSystem->WorkingDirectory();
 108     gSystem->ChangeDirectory(pararchivename);
 109
 110     // check for BUILD.sh and execute
 111     if (!gSystem->AccessPathName("PROOF-INF/BUILD.sh")) {
 112       printf("*******************************\n");
 113       printf("*** Building PAR archive    ***\n");
 114       printf("*******************************\n");
 115
 116       if (gSystem->Exec("PROOF-INF/BUILD.sh")) {
 117         Error("runProcess","Cannot Build the PAR Archive! - Abort!");
 118         return -1;
 119       }
 120     }
 121     // check for SETUP.C and execute
 122     if (!gSystem->AccessPathName("PROOF-INF/SETUP.C")) {
 123       printf("*******************************\n");
 124       printf("*** Setup PAR archive       ***\n");
 125       printf("*******************************\n");
 126       gROOT->Macro("PROOF-INF/SETUP.C");
 127     }
 128
 129     gSystem->ChangeDirectory("../");
 130   }
 131
 132   return 1;
 133 }