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   gROOT->LoadMacro("CreateESDChain.C");
  72   chain = CreateESDChain(chainlistfile,500);
  73
  74   //____________________________________________//
  75   // Make the analysis manager
  76   AliAnalysisManager *mgr = new AliAnalysisManager("TestManager");
  77   AliESDInputHandler* esdH = new AliESDInputHandler;
  78   esdH->SetInactiveBranches("FMD CaloCluster");
  79   mgr->SetInputEventHandler(esdH);
  80   //____________________________________________//
  81   // 1st Pt task
  82   AliAnalysisTaskFemto *task1 = new AliAnalysisTaskFemto("TaskFemto");
  83
  84   mgr->AddTask(task1);
  85
  86   // Create containers for input/output
  87   AliAnalysisDataContainer *cinput1 = mgr->CreateContainer("cchain1",TChain::Class(),AliAnalysisManager::kInputContainer);
  88   AliAnalysisDataContainer *coutput1 = mgr->CreateContainer("clist1", TList::Class(),AliAnalysisManager::kOutputContainer,"Femto.ESD.root");
  89
  90   //____________________________________________//
  91   cinput1->SetData(chain);
  92   mgr->ConnectInput(task1,0,cinput1);
  93   mgr->ConnectOutput(task1,0,coutput1);
  94
  95   if (!mgr->InitAnalysis()) return;
  96   mgr->PrintStatus();
  97   mgr->StartAnalysis("local",chain);
  98
  99   timer.Stop();
 100   timer.Print();
 101 }
 102
 103 Int_t setupPar(const char* pararchivename) {
 104   ///////////////////
 105   // Setup PAR File//
 106   ///////////////////
 107   if (pararchivename) {
 108     char processline[1024];
 109     sprintf(processline,".! tar xvzf %s.par",pararchivename);
 110     gROOT->ProcessLine(processline);
 111     const char* ocwd = gSystem->WorkingDirectory();
 112     gSystem->ChangeDirectory(pararchivename);
 113
 114     // check for BUILD.sh and execute
 115     if (!gSystem->AccessPathName("PROOF-INF/BUILD.sh")) {
 116       printf("*******************************\n");
 117       printf("*** Building PAR archive    ***\n");
 118       printf("*******************************\n");
 119
 120       if (gSystem->Exec("PROOF-INF/BUILD.sh")) {
 121         Error("runProcess","Cannot Build the PAR Archive! - Abort!");
 122         return -1;
 123       }
 124     }
 125     // check for SETUP.C and execute
 126     if (!gSystem->AccessPathName("PROOF-INF/SETUP.C")) {
 127       printf("*******************************\n");
 128       printf("*** Setup PAR archive       ***\n");
 129       printf("*******************************\n");
 130       gROOT->Macro("PROOF-INF/SETUP.C");
 131     }
 132
 133     gSystem->ChangeDirectory("../");
 134   }
 135
 136   return 1;
 137 }