PWG2/FEMTOSCOPY/macros/Local/runLocal.C

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