ANALYSIS/macros/QAtrainsim.C

   1 #include "Riostream.h"
   2 void LoadLibraries();
   3 void AddAnalysisTasks();
   4 void QAmerge(const char *, Int_t);
   5
   6 Int_t iCollisionType = 1; // 0=pp, 1=PbPb
   7 // Trigger mask.
   8
   9 UInt_t kTriggerInt = AliVEvent::kAnyINT;
  10 UInt_t kTriggerMuonAll = AliVEvent::kMUL7 | AliVEvent::kMUSH7 | AliVEvent::kMUU7 | AliVEvent::kMUS7
  11                        | AliVEvent::kMUSPB | AliVEvent::kMUSHPB | AliVEvent::kMuonLikePB | AliVEvent::kMuonUnlikePB;
  12 UInt_t kTriggerMuonBarell = AliVEvent::kMUU7;
  13 UInt_t kTriggerEMC   = AliVEvent::kEMC7 | AliVEvent::kEMCEJE | AliVEvent::kEMCEGA;
  14 UInt_t kTriggerHM   = AliVEvent::kHighMult;
  15 // Main trigger mask used:
  16 UInt_t kTriggerMask = kTriggerInt;
  17
  18 Int_t runNumbers[5] = {158626};
  19
  20 Bool_t doCDBconnect   = 1;
  21 Bool_t doEventStat    = 1;
  22 Bool_t doCentrality   = 1;
  23 Bool_t doQAsym        = 1;
  24 Bool_t doVZERO        = 1;   // there is a 2nd file
  25 Bool_t doVZEROPbPb    = 1;
  26 Bool_t doVertex       = 1;
  27 Bool_t doSPD          = 1;   // needs RP
  28 Bool_t doTPC          = 1;
  29 Bool_t doHLT          = 1;
  30 Bool_t doSDD          = 1;   // needs RP
  31 Bool_t doSSDdEdx      = 1;
  32
  33 Bool_t doTRD          = 1;
  34 Bool_t doITS          = 1;
  35 Bool_t doITSsaTracks  = 1;
  36 Bool_t doITSalign     = 1;
  37 Bool_t doCALO         = 1;
  38 Bool_t doMUONTrig     = 1;
  39 Bool_t doImpParRes    = 1;
  40 Bool_t doMUON         = 1;
  41 Bool_t doTOF          = 1;
  42 Bool_t doHMPID        = 1;
  43 Bool_t doT0           = 1;
  44 Bool_t doZDC          = 1;
  45 Bool_t doPIDResponse  = 1;
  46 Bool_t doPIDqa        = 1; //new
  47 Bool_t doFMD          = 1; // new
  48 Bool_t doPHOS         = 1; // new
  49 Bool_t doEMCAL        = 1;
  50 Bool_t doFBFqa        = 1; // new - not ported yet to revision
  51
  52                // Debug level
  53 Int_t       debug_level        = 1;        // Debugging
  54 Int_t       run_number = 0;
  55
  56 void QAtrain(Int_t run = 0,
  57              const char *xmlfile   = "wn.xml",
  58              Int_t  stage          = 0) /*0 = QA train, 1...n - merging stage*/
  59 {
  60   run_number = run;
  61
  62   TGrid::Connect("alien://");
  63   if (!gGrid || !gGrid->IsConnected()) {
  64     ::Error("QAtrain", "No grid connection");
  65     return;
  66   }
  67   // Set temporary merging directory to current one
  68   gSystem->Setenv("TMPDIR", gSystem->pwd());
  69   // Set temporary compilation directory to current one
  70   gSystem->SetBuildDir(gSystem->pwd(), kTRUE);
  71   // Load libraries
  72   LoadLibraries();
  73   printf("Include path: %s\n", gSystem->GetIncludePath());
  74   // Create manager
  75   AliAnalysisManager *mgr  = new AliAnalysisManager("PilotAnalysis_sim", "Production train");
  76   mgr->SetRunFromPath(run_number);
  77   // Input handler
  78   AliESDInputHandlerRP *esdHandler = new AliESDInputHandlerRP();
  79   esdHandler->SetReadFriends(kTRUE);
  80   esdHandler->SetActiveBranches("ESDfriend");
  81   mgr->SetInputEventHandler(esdHandler);
  82   mgr->SetDebugLevel(debug_level);
  83
  84   // Monte Carlo handler
  85   AliMCEventHandler* mcHandler = new AliMCEventHandler();
  86   mgr->SetMCtruthEventHandler(mcHandler);
  87   mcHandler->SetPreReadMode(1);
  88   mcHandler->SetReadTR(kTRUE);
  89
  90   // AnalysisTasks
  91   AddAnalysisTasks();
  92   if (stage>0) {
  93     QAmerge(xmlfile, stage);
  94     return;
  95   }
  96   // Input chain
  97   TChain *chain = new TChain("esdTree");
  98   chain->Add("AliESDs.root");
  99   TStopwatch timer;
 100   timer.Start();
 101   if (mgr->InitAnalysis()) {
 102     mgr->PrintStatus();
 103     mgr->SetSkipTerminate(kTRUE);
 104 //    mgr->SetNSysInfo(1);
 105     mgr->StartAnalysis("local", chain);
 106   }
 107   timer.Print();
 108 }
 109
 110 void LoadLibraries()
 111 {
 112   gSystem->SetIncludePath("-I. -I$ROOTSYS/include -I$ALICE_ROOT/include -I$ALICE_ROOT -I$ALICE_ROOT/ITS -I$ALICE_ROOT/TRD -I$ALICE_ROOT/PWGPP -I$ALICE_ROOT/PWGPP/TRD");
 113   gSystem->Load("libANALYSIS");
 114   gSystem->Load("libANALYSISalice");
 115   gSystem->Load("libCORRFW");
 116   gSystem->Load("libTENDER");
 117   gSystem->Load("libPWGPP.so");
 118   gSystem->Load("libAliHLTTrigger.so");
 119
 120   if (doEMCAL || doCALO || doPHOS) {
 121      gSystem->Load("libEMCALUtils");
 122      gSystem->Load("libPHOSUtils");
 123      gSystem->Load("libPWGCaloTrackCorrBase");
 124      gSystem->Load("libPWGGACaloTrackCorrelations");
 125      gSystem->Load("libPWGGACaloTasks");
 126      gSystem->Load("libPWGGAPHOSTasks");
 127      gSystem->Load("libPWGGAEMCALTasks");
 128   }
 129   if(doMUON || doMUONTrig) {
 130      gSystem->Load("libPWGmuon");
 131      gSystem->Load("libPWGPPMUONlite");
 132      gSystem->Load("libPWGmuondep");
 133   }
 134   if (doFMD) {
 135      gSystem->Load("libPWGLFforward2");
 136   }
 137 }
 138
 139 void AddAnalysisTasks()
 140 {
 141   AliAnalysisManager *mgr = AliAnalysisManager::GetAnalysisManager();
 142   mgr->SetCommonFileName("QAresults.root");
 143   // Statistics task
 144   mgr->AddStatisticsTask(kTriggerMask);
 145   //
 146   // CDB connection
 147   //
 148   if (doCDBconnect) {
 149     gROOT->LoadMacro("$ALICE_ROOT/PWGPP/PilotTrain/AddTaskCDBconnect.C");
 150     AliTaskCDBconnect *taskCDB = AddTaskCDBconnect();
 151     if (!taskCDB) return;
 152     AliCDBManager *cdb = AliCDBManager::Instance();
 153     cdb->SetDefaultStorage("raw://");
 154     taskCDB->SetRunNumber(run_number);
 155   }
 156
 157   //
 158   // Event Statistics (Jan Fiete)
 159   //
 160   if (doEventStat) {
 161       gROOT->LoadMacro("$ALICE_ROOT/ANALYSIS/macros/AddTaskPhysicsSelection.C");
 162       AliPhysicsSelectionTask* physSelTask = AddTaskPhysicsSelection(kTRUE /*MC*/);
 163   }
 164
 165   //
 166   // Centrality (A. Toia)
 167   //
 168   if (doCentrality) {
 169      if (!iCollisionType) {
 170         printf("Disabling centrality task for p-p\n");
 171         doCentrality = kFALSE;
 172      } else {
 173         gROOT->LoadMacro("$ALICE_ROOT/ANALYSIS/macros/AddTaskCentrality.C");
 174         AliCentralitySelectionTask *taskCentrality = AddTaskCentrality();
 175         taskCentrality->SetMCInput();
 176      }
 177   }
 178
 179   // Vertexing (A. Dainese)
 180   //
 181   if (doVertex) {
 182     gROOT->LoadMacro("$ALICE_ROOT/PWGPP/macros/AddTaskVertexESD.C");
 183     // Specific setting for MC
 184     AliAnalysisTaskVertexESD* taskvertexesd =  AddTaskVertexESD(kTRUE, kTriggerMask);
 185     taskvertexesd->SelectCollisionCandidates(kTriggerMask);
 186   }
 187
 188   // TPC QA (E. Sicking)
 189   //
 190   if (doQAsym) {
 191   // offline trigger in AddTask
 192     gROOT->LoadMacro("$ALICE_ROOT/PWGPP/PilotTrain/AddTaskQAsym.C");
 193     AliAnalysisTaskSE * taskqasim = AddTaskQAsym(0, kTriggerMask, kTriggerHM, kTriggerEMC, kTriggerMuonBarell);
 194   }
 195   //
 196   // VZERO QA  (C. Cheshkov)
 197   //
 198   if (doVZERO) {
 199     gROOT->LoadMacro("$ALICE_ROOT/PWGPP/PilotTrain/AddTaskVZEROQA.C");
 200     AliAnalysisTaskSE * taskv0qa = AddTaskVZEROQA(0);
 201 //  taskv0qa->SelectCollisionCandidates();
 202   }
 203   if (doVZEROPbPb && iCollisionType==1) {
 204     gROOT->LoadMacro("$ALICE_ROOT/PWGPP/VZERO/AddTaskVZEROPbPb.C");
 205     AliAnaVZEROPbPb* taskV0PbPb = (AliAnaVZEROPbPb*)AddTaskVZEROPbPb(0);
 206     taskV0PbPb->SetClassesNames("CTRUE-,C0HWU-,CPBI2WU-,CPBI2-,CPBI2WU_B1-,CPBI2_B1-,CPBI1WU-,CPBI1-,CVHNWU-,CVHN-,CVHN_R2-,CVHNWU_R2-,CVLNWU-,CVLN-,CVLN_R1-,CVLN_B2-,CVLNWU_R1-,CVLNWU_B2-,CSEMI_R1-,CSEMIWU_R1-,CCENT_R2-,CCENTWU_R2-");
 207   }
 208   //
 209   // TPC (Jacek Otwinowski & Michael Knichel)
 210   //
 211   //
 212   // Optionally MC information can be used by setting the 1st argument to true
 213   // Optionally friends information can be switched off by setting the 2st argument
 214   // to false
 215   // Optionally highMult axis can be used by setting the 3st argument to true (for PbPb)
 216   if (doTPC) {
 217     gROOT->LoadMacro("$ALICE_ROOT/PWGPP/TPC/macros/AddTaskPerformanceTPCdEdxQA.C");
 218     AliPerformanceTask *tpcQA = 0;
 219     if (iCollisionType) {
 220        // High multiplicity Pb-Pb
 221        tpcQA = AddTaskPerformanceTPCdEdxQA(kTRUE, kTRUE, kTRUE);
 222     } else {
 223       // Low multiplicity (pp)
 224        tpcQA = AddTaskPerformanceTPCdEdxQA(kTRUE, kTRUE, kFALSE);
 225     }
 226     tpcQA->SelectCollisionCandidates(kTriggerMask);
 227   }
 228
 229   // HLT (Alberica Toia)
 230   if (doHLT) {
 231     gROOT->LoadMacro("$ALICE_ROOT/PWGPP/TPC/macros/AddTaskPerformanceTPCdEdxQA.C");
 232     AliPerformanceTask *hltQA = AddTaskPerformanceTPCdEdxQA(kTRUE, kTRUE, kFALSE,0,kTRUE);
 233     hltQA->SelectCollisionCandidates(kTriggerMask);
 234   }
 235   //
 236   // SPD (A. Mastroserio)
 237   //
 238   if (doSPD) {
 239     gROOT->LoadMacro("$ALICE_ROOT/PWGPP/PilotTrain/AddTaskSPDQA.C");
 240     AliAnalysisTaskSPD* taskspdqa = (AliAnalysisTaskSPD*)AddTaskSPDQA();
 241     // Request from Annalisa
 242     if (iCollisionType) taskspdqa->SetHeavyIonMode();
 243     taskspdqa->SelectCollisionCandidates(kTriggerMask);
 244     taskspdqa->SetOCDBInfo(run_number, "raw://");
 245   }
 246   //
 247   // SDD (F. Prino)
 248   //
 249   if (doSDD) {
 250     gROOT->LoadMacro("$ALICE_ROOT/PWGPP/PilotTrain/AddSDDPoints.C");
 251     AliAnalysisTaskSE* tasksdd = AddSDDPoints();
 252     tasksdd->SelectCollisionCandidates(kTriggerMask);
 253   }
 254   //
 255   // SSD dEdx (Marek Chojnacki)
 256   //
 257   if (doSSDdEdx) {
 258     gROOT->LoadMacro("$ALICE_ROOT/PWGPP/PilotTrain/AddTaskdEdxSSDQA.C");
 259     AliAnalysisTaskSE* taskssddedx = AddTaskdEdxSSDQA();
 260     taskssddedx->SelectCollisionCandidates(kTriggerMask);
 261   }
 262
 263   //
 264   // ITS
 265   //
 266   if (doITS) {
 267   // hardcoded non-zero trigger mask
 268       gROOT->LoadMacro("$ALICE_ROOT/PWGPP/macros/AddTaskPerformanceITS.C");
 269       AliAnalysisTaskITSTrackingCheck *itsQA = 0;
 270       AliAnalysisTaskITSTrackingCheck *itsQACent0010 = 0;
 271       AliAnalysisTaskITSTrackingCheck *itsQACent3050 = 0;
 272       AliAnalysisTaskITSTrackingCheck *itsQACent6080 = 0;
 273       if(iCollisionType==0) {
 274         itsQA = AddTaskPerformanceITS(kTRUE);
 275       } else {
 276         itsQA = AddTaskPerformanceITS(kTRUE);
 277         itsQACent0010 = AddTaskPerformanceITS(kTRUE,kFALSE,kFALSE,3500,10000);
 278         itsQACent3050 = AddTaskPerformanceITS(kTRUE,kFALSE,kFALSE,590,1570);
 279         itsQACent6080 = AddTaskPerformanceITS(kTRUE,kFALSE,kFALSE,70,310);
 280       }
 281   }
 282   //
 283   // ITS saTracks, align (F.Prino)
 284   //
 285   if (doITSsaTracks) {
 286   // offline trigger in AddTask
 287      gROOT->LoadMacro("$ALICE_ROOT/PWGPP/macros/AddTaskITSsaTracks.C");
 288      AliAnalysisTaskITSsaTracks *itssaTracks = AddTaskITSsaTracks(kTRUE,kFALSE);
 289      itssaTracks->SelectCollisionCandidates(kTriggerMask);
 290   }
 291   if (doITSalign) {
 292   // no offline trigger selection
 293      gROOT->LoadMacro("$ALICE_ROOT/PWGPP/macros/AddTaskITSAlign.C");
 294      AliAnalysisTaskITSAlignQA *itsAlign = AddTaskITSAlign(0,2011);
 295   }
 296   //
 297   // TRD (Alex Bercuci, M. Fasel)
 298   //
 299   if(doTRD) {
 300   // no offline trigger selection
 301       gROOT->LoadMacro("$ALICE_ROOT/PWGPP/macros/AddTrainPerformanceTRD.C");
 302       // steer individual TRD tasks
 303       Bool_t
 304       doCheckESD(kTRUE),  // AliTRDcheckESD
 305       doCheckDET(kTRUE),  // AliTRDcheckDET
 306       doEffic(kTRUE),     // AliTRDefficiency
 307       doResolution(kTRUE),// AliTRDresolution
 308       doCheckPID(kTRUE),  // AliTRDcheckPID
 309       doV0Monitor(kFALSE);// AliTRDv0Monitor
 310       AddTrainPerformanceTRD(Translate(doCheckESD, doCheckDET, doEffic, doResolution, doCheckPID, doV0Monitor));
 311   }
 312
 313   //
 314   // ZDC (Chiara Oppedisano)
 315   //
 316   if(doZDC) {
 317   // hardcoded kMB trigger mask
 318      gROOT->LoadMacro("$ALICE_ROOT/PWGPP/ZDC/AddTaskZDCQA.C");
 319      AliAnalysisTaskSE *taskZDC = AddTaskZDCQA();
 320      taskZDC->SelectCollisionCandidates(kTriggerMask);
 321   }
 322   //
 323   // Calorimetry (Gustavo Conesa)
 324   //
 325
 326   if(doCALO) {
 327       gROOT->LoadMacro("$ALICE_ROOT/PWGGA/CaloTrackCorrelations/macros/QA/AddTaskCalorimeterQA.C");
 328       AliAnalysisTaskCaloTrackCorrelation *taskCaloQA = AddTaskCalorimeterQA("ESD", 2011, kFALSE, kTRUE);
 329       taskCaloQA->SetDebugLevel(0);
 330       // offline mask set in AddTask to kMB
 331       taskCaloQA->SelectCollisionCandidates(kTriggerMask);
 332       // Add a new calo task with EMC1 trigger only
 333       taskCaloQA = AddTaskCalorimeterQA("ESD", 2011, kFALSE, kTRUE, "", "EMC7");
 334       taskCaloQA->SetDebugLevel(0);
 335       taskCaloQA->SelectCollisionCandidates(kTriggerEMC);
 336   }
 337
 338   //
 339   // Muon Trigger
 340   //
 341
 342   if(doMUONTrig) {
 343   // no offline trigger selection
 344       gROOT->LoadMacro("$ALICE_ROOT/PWGPP/macros/AddTaskMTRchamberEfficiency.C");
 345       AliAnalysisTaskTrigChEff *taskMuonTrig = AddTaskMTRchamberEfficiency();
 346   }
 347
 348   //
 349   // Impact parameter resolution (xianbao.yuan@pd.infn.it, andrea.dainese@pd.infn.it)
 350   //
 351   if (doImpParRes) {
 352     gROOT->LoadMacro("$ALICE_ROOT/PWGPP/macros/AddTaskImpParRes.C");
 353     AliAnalysisTaskSE* taskimpparres=0;
 354     // Specific setting for MC
 355     if(iCollisionType==0) {
 356        taskimpparres= AddTaskImpParRes(kTRUE);
 357     } else {
 358        taskimpparres= AddTaskImpParRes(kTRUE,-1,kTRUE,kFALSE);
 359     }
 360     taskimpparres->SelectCollisionCandidates(kTriggerMask);
 361   }
 362   //
 363   // MUON QA (Philippe Pillot)
 364   //
 365   if (doMUON) {
 366   // trigger analysis internal
 367     gROOT->LoadMacro("$ALICE_ROOT/PWGPP/PilotTrain/AddTaskMuonQA.C");
 368     AliAnalysisTaskSE* taskmuonqa= AddTaskMuonQA();
 369   }
 370   //
 371   // TOF (Francesca Bellini)
 372   //
 373   if (doTOF) {
 374     gROOT->LoadMacro("$ALICE_ROOT/PWGPP/TOF/AddTaskTOFQA.C");
 375     AliAnalysisTaskTOFqa *tofQA = AddTaskTOFQA();
 376     tofQA->SelectCollisionCandidates(kTriggerMask);
 377   }
 378    //
 379   // PIDResponse(JENS)
 380   //
 381   if (doPIDResponse) {
 382     gROOT->LoadMacro("$ALICE_ROOT/ANALYSIS/macros/AddTaskPIDResponse.C");
 383     AliAnalysisTaskPIDResponse *PIDResponse = AddTaskPIDResponse(kTRUE);
 384     PIDResponse->SelectCollisionCandidates(kTriggerMask);
 385   }
 386
 387   //
 388   // PIDqa(JENS)
 389   //
 390   if (doPIDqa) {
 391     gROOT->LoadMacro("$ALICE_ROOT/ANALYSIS/macros/AddTaskPIDqa.C");
 392     AliAnalysisTaskPIDqa *PIDQA = AddTaskPIDqa();
 393     PIDQA->SelectCollisionCandidates(kTriggerMask);
 394   }
 395
 396   //
 397   // HMPID QA (Giacomo Volpe)
 398   //
 399   if (doHMPID) {
 400     gROOT->LoadMacro("$ALICE_ROOT/PWGPP/HMPID/AddTaskHmpidQA.C");
 401     AliAnalysisTaskSE* taskhmpidqa= AddTaskHmpidQA(kTRUE);
 402       // offline mask set in AddTask to kMB
 403     taskhmpidqa->SelectCollisionCandidates(kTriggerMask);
 404   }
 405   // T0 QA (Alla Mayevskaya)
 406   if (doT0) {
 407   // no offline trigger selection
 408     gROOT->LoadMacro("$ALICE_ROOT/PWGPP/T0/AddTaskT0QA.C");
 409     AliT0AnalysisTaskQA* taskt0qa= AddTaskT0QA();
 410     taskt0qa->SelectCollisionCandidates(kTriggerMask);
 411   }
 412   // FMD QA (Christian Holm Christiansen)
 413   if (doFMD) {
 414     gROOT->LoadMacro("$ALICE_ROOT/PWGLF/FORWARD/analysis2/AddTaskForwardQA.C");
 415     // Parameters: usemc, usecentrality
 416     AliAnalysisTaskSE *forwardQA = (AliAnalysisTaskSE *)AddTaskForwardQA(kTRUE, (Bool_t)doCentrality);
 417     // No offline trigger config. needed (see #84077)
 418   }
 419    //
 420   // PHOS QA (Boris Polishchuk)
 421   //
 422   if (doPHOS) {
 423     gROOT->LoadMacro("$ALICE_ROOT/PWGGA/PHOSTasks/CaloCellQA/macros/AddTaskCaloCellsQA.C");
 424     AliAnalysisTaskCaloCellsQA *taskPHOSCellQA1 = AddTaskCaloCellsQA(4, 1, NULL,"PHOSCellsQA_AnyInt");
 425     taskPHOSCellQA1->SelectCollisionCandidates(kTriggerMask);
 426     taskPHOSCellQA1->GetCaloCellsQA()->SetClusterEnergyCuts(0.3,0.3,1.0);
 427     AliAnalysisTaskCaloCellsQA *taskPHOSCellQA2 = AddTaskCaloCellsQA(4, 1, NULL,"PHOSCellsQA_PHI7");
 428     taskPHOSCellQA2->SelectCollisionCandidates(AliVEvent::kPHI7);
 429     taskPHOSCellQA2->GetCaloCellsQA()->SetClusterEnergyCuts(0.3,0.3,1.0);
 430     // Pi0 QA fo PbPb
 431     if (iCollisionType) {
 432       gROOT->LoadMacro("$ALICE_ROOT/PWGGA/PHOSTasks/PHOS_PbPbQA/macros/AddTaskPHOSPbPb.C");
 433       AliAnalysisTaskPHOSPbPbQA* phosPbPb = AddTaskPHOSPbPbQA(0);
 434     }
 435   }
 436   //
 437   // EMCAL QA (Gustavo Conesa)
 438   //
 439   if (doEMCAL) {
 440      gROOT->LoadMacro("$ALICE_ROOT/PWGGA/EMCALTasks/macros/AddTaskEMCALTriggerQA.C");
 441      AliAnalysisTaskEMCALTriggerQA *emctrig = AddTaskEMCALTriggerQA();
 442   }
 443   //
 444   // FLOW and BF QA (C.Perez && A.Rodriguez)
 445   //
 446   if (doFBFqa) {
 447     gROOT->LoadMacro("$ALICE_ROOT/PWGPP/macros/AddTaskFBFqa.C");
 448     AliAnalysisTaskSE *qaFBFMB = (AliAnalysisTaskSE*) AddTaskFBFqa("qaFBFmb",kFALSE);
 449     qaFBFMB->SelectCollisionCandidates(AliVEvent::kMB);
 450     AliAnalysisTaskSE *qaFBFSC = (AliAnalysisTaskSE*) AddTaskFBFqa("qaFBFsc",kFALSE);
 451     qaFBFSC->SelectCollisionCandidates(AliVEvent::kSemiCentral);
 452     AliAnalysisTaskSE *qaFBFCE = (AliAnalysisTaskSE*) AddTaskFBFqa("qaFBFce",kFALSE);
 453     qaFBFCE->SelectCollisionCandidates(AliVEvent::kCentral);
 454   }
 455 }
 456
 457 void QAmerge(const char *dir, Int_t stage)
 458 {
 459 // Merging method
 460   TStopwatch timer;
 461   timer.Start();
 462   TString outputDir = dir;
 463   TString outputFiles = "QAresults.root,EventStat_temp.root,RecoQAresults.root";
 464   TString mergeExcludes = "";
 465   TObjArray *list = outputFiles.Tokenize(",");
 466   TIter *iter = new TIter(list);
 467   TObjString *str;
 468   TString outputFile;
 469   Bool_t merged = kTRUE;
 470   while((str=(TObjString*)iter->Next())) {
 471     outputFile = str->GetString();
 472     // Skip already merged outputs
 473     if (!gSystem->AccessPathName(outputFile)) {
 474        printf("Output file <%s> found. Not merging again.",outputFile.Data());
 475        continue;
 476     }
 477     if (mergeExcludes.Contains(outputFile.Data())) continue;
 478     merged = AliAnalysisAlien::MergeOutput(outputFile, outputDir, 10, stage);
 479     if (!merged) {
 480        printf("ERROR: Cannot merge %s\n", outputFile.Data());
 481        return;
 482     }
 483   }
 484   // read the analysis manager from file
 485   if (!outputDir.Contains("Stage")) {
 486     ofstream out;
 487     out.open("outputs_valid", ios::out);
 488     out.close();
 489     return;
 490   }
 491   AliAnalysisManager *mgr = AliAnalysisManager::GetAnalysisManager();
 492   mgr->SetRunFromPath(mgr->GetRunFromAlienPath(dir));
 493   mgr->SetSkipTerminate(kFALSE);
 494   if (!mgr->InitAnalysis()) return;
 495   mgr->PrintStatus();
 496   AliLog::SetGlobalLogLevel(AliLog::kError);
 497   TTree *tree = NULL;
 498   mgr->StartAnalysis("gridterminate", tree);
 499   ofstream out;
 500   out.open("outputs_valid", ios::out);
 501   out.close();
 502   timer.Print();
 503 }