TENDER becomes Tender, removing .so
[u/mrichter/AliRoot.git] / ANALYSIS / macros / QAtrain_duo.C
1 // QAtrain_duo - assumes having 2 sets of input data files using 2 suffixes:
2 //    AliESDs_barrel.root AliESDfriends_barrel.root ...  suffix = _barrel
3 //    AliESDs_outer.root AliESDfriends_outer.root ...    suffix = _outer
4 // HOW TO RUN (for run_no=158285):
5 // QA job: should call:
6 //    ln -fs AliESDs_barrel.root AliESDs.root
7 //    ln -fs AliESDfriends_barrel.root AliESDfriends.root
8 //    ... same for all specific reco files
9 //    aliroot -b -q QAtrain_duo.C\(\"_barrel\"\,158285\) producing:
10 //      -> EventStat_temp_barrel.root QAresults_barrel.root 400_400_0_350.stat
11 //
12 //    ln -fs AliESDs_outer.root AliESDs.root
13 //    ln -fs AliESDfriends_outer.root AliESDfriends.root
14 //    ... same for all specific reco files
15 //    aliroot QAtrain_duo.C\(\"_outer\"\,158285\) producing:
16 //      -> EventStat_temp_outer.root QAresults_outer.root 400_400_0_380.stat
17 //
18 // Each merging job:
19 // for stages i < 5
20 //    aliroot -b -q QAtrain_duo.C\(\"_barrel\"\,158285\,\"wn.xml\"\,i\)
21 //    aliroot -b -q QAtrain_duo.C\(\"_outer\"\,158285\,\"wn.xml\"\,i\)
22 //      -> same output files as the QA jobs (except .stat)
23 // for stage 5
24 //    aliroot -b -q QAtrain_duo.C\(\"_barrel\"\,158285\,\"Stage_5.xml\"\,5\)
25 //      -> event_stat_barrel.root trending_barrel.root 145230_145230_0_120345.stat
26 //    aliroot -b -q QAtrain_duo.C\(\"_outer\"\,158285\,\"Stage_5.xml\"\,5\)
27 //      -> event_stat_outer.root trending_outer.root 145230_145230_0_132897.stat
28 #include "Riostream.h"
29 void LoadLibraries();
30 void AddAnalysisTasks(const char *, const char *); 
31 void QAmerge(const char *,const char *, Int_t);
32 void ProcessEnvironmentVars();
33 void SetDoQA(Bool_t &var, Bool_t onoff, const char* name);
34 Bool_t CheckEnvS(const char* var,TString& envString);
35
36 Int_t iCollisionType = 0; // 0=pp, 1=PbPb
37 // Trigger mask.
38
39 UInt_t kTriggerInt = AliVEvent::kAnyINT;
40 UInt_t kTriggerMuonBarrel = AliVEvent::kMUU7 | AliVEvent::kMuonUnlikeLowPt8 | AliVEvent::kMuonUnlikeLowPt0;
41 UInt_t kTriggerEMC   = AliVEvent::kEMC7 | AliVEvent::kEMC8 | AliVEvent::kEMCEJE | AliVEvent::kEMCEGA;
42 UInt_t kTriggerHM   = AliVEvent::kHighMult;
43 // Main trigger mask used:
44 UInt_t kTriggerMask = kTriggerInt;
45
46 Int_t runNumbers[5] = {158626};
47
48 Bool_t doStatistics   = 1;
49 Bool_t doCDBconnect   = 1;
50 Bool_t doEventStat    = 1;
51 Bool_t doCentrality   = 0;
52 Bool_t doQAsym        = 0;
53 Bool_t doVZERO        = 1;   // there is a 2nd file
54 Bool_t doVZEROPbPb    = 1; 
55 Bool_t doVertex       = 1;
56 Bool_t doSPD          = 1;   // needs RP   
57 Bool_t doTPC          = 1;
58 Bool_t doHLT          = 1;
59 Bool_t doSDD          = 1;   // needs RP
60 Bool_t doSSDdEdx      = 1;
61
62 Bool_t doTRD          = 1;
63 Bool_t doITS          = 1;
64 Bool_t doITSsaTracks  = 1; 
65 Bool_t doITSalign     = 1;  
66 Bool_t doCALO         = 1;
67 Bool_t doMUONTrig     = 1;
68 Bool_t doImpParRes    = 1;
69 Bool_t doMUON         = 1;
70 Bool_t doTOF          = 1;
71 Bool_t doHMPID        = 1;
72 Bool_t doT0           = 1;
73 Bool_t doZDC          = 1;
74 Bool_t doPIDResponse  = 1;
75 Bool_t doPIDqa        = 1; //new
76 Bool_t doFMD          = 1; // new
77 Bool_t doPHOS         = 1; // new
78 Bool_t doPHOSTrig     = 1; // new
79 Bool_t doEMCAL        = 0;
80 Bool_t doFBFqa        = 1; // new - not ported yet to revision
81
82 Bool_t doTaskFilteredTree        = 0;      // high pt filter task
83
84                // Debug level
85 Int_t       debug_level        = 1;        // Debugging
86 Int_t       run_number = 0;
87
88
89 void QAtrain_duo(const char *suffix="", Int_t run = 0, 
90              const char *xmlfile   = "wn.xml",
91              Int_t  stage          = 0, /*0 = QA train, 1...n - merging stage*/
92              const char *cdb     = "raw://")
93 //             const char *cdb     = "local://$ALICE_ROOT/OCDB")
94 {
95   run_number = run;
96   TString ss(suffix);
97   ss.ToLower();
98   Bool_t ibarrel = (ss.Contains("barrel"))?kTRUE:kFALSE;
99
100   ProcessEnvironmentVars();
101
102   TString cdbString(cdb);
103   if (cdbString.Contains("raw://"))
104   {
105     TGrid::Connect("alien://");
106     if (!gGrid || !gGrid->IsConnected()) {
107       ::Error("QAtrain", "No grid connection");
108       return;
109     }  
110   }
111   // Set temporary merging directory to current one
112   gSystem->Setenv("TMPDIR", gSystem->pwd());
113   // Set temporary compilation directory to current one
114   gSystem->SetBuildDir(gSystem->pwd(), kTRUE);
115   // Load libraries
116   LoadLibraries();
117   printf("Include path: %s\n", gSystem->GetIncludePath());
118   // Create manager
119   AliAnalysisManager *mgr  = new AliAnalysisManager("PilotAnalysis", "Production train");
120   mgr->SetRunFromPath(run_number);
121   mgr->SetCacheSize(0);
122   // Input handler
123   AliESDInputHandlerRP *esdHandler = new AliESDInputHandlerRP();
124   if (ibarrel) {
125     esdHandler->SetReadFriends(kTRUE);
126     esdHandler->SetActiveBranches("ESDfriend");
127   }
128   mgr->SetInputEventHandler(esdHandler);
129   mgr->SetDebugLevel(debug_level);
130   mgr->SetFileInfoLog("fileinfo.log"); 
131   
132   // AnalysisTasks
133   AddAnalysisTasks(suffix, cdb);
134   if (stage>0) {
135     QAmerge(suffix, xmlfile, stage);
136     return;
137   }   
138   // Input chain
139   TChain *chain = new TChain("esdTree");
140   chain->Add("AliESDs.root");
141   TStopwatch timer;
142   timer.Start();
143   if (mgr->InitAnalysis()) {                                                                                                              
144     mgr->PrintStatus(); 
145     mgr->SetSkipTerminate(kTRUE);
146  //   mgr->SetNSysInfo(100);
147     mgr->StartAnalysis("local", chain);
148   }
149   timer.Print();
150 }
151
152 void LoadLibraries()
153 {
154   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 -I$ALICE_ROOT/ANALYSIS");
155   gSystem->Load("libANALYSIS");
156   gSystem->Load("libANALYSISalice");
157   gSystem->Load("libESDfilter");
158   gSystem->Load("libCORRFW");
159   gSystem->Load("libTender");
160   gSystem->Load("libPWGPP");
161   gSystem->Load("libAliHLTTrigger");
162
163   if (doEMCAL || doPHOS || doCALO) {
164      gSystem->Load("libEMCALUtils");
165      gSystem->Load("libPHOSUtils");
166      gSystem->Load("libPWGCaloTrackCorrBase");
167      gSystem->Load("libPWGGACaloTrackCorrelations");
168      gSystem->Load("libPWGGACaloTasks");
169      gSystem->Load("libPWGGAPHOSTasks");
170      gSystem->Load("libPWGTools");
171      gSystem->Load("libPWGEMCAL");
172      gSystem->Load("libPWGGAEMCALTasks");
173   }  
174   if(doMUON || doMUONTrig) {
175      gSystem->Load("libPWGmuon");
176      gSystem->Load("libPWGPPMUONlite");
177      gSystem->Load("libPWGmuondep");
178   }
179   if (doFMD) {
180      gSystem->Load("libPWGLFforward2");
181   }      
182 }
183
184 void AddAnalysisTasks(const char *suffix, const char *cdb_location)
185 {
186   TString ss(suffix);
187   ss.ToLower();
188   Bool_t ibarrel = (ss.Contains("barrel"))?kTRUE:kFALSE;
189   Bool_t iall = ss.IsNull();
190   AliAnalysisManager *mgr = AliAnalysisManager::GetAnalysisManager();
191   mgr->SetCommonFileName(Form("QAresults%s.root",suffix));
192   // Statistics task
193   if (doStatistics) mgr->AddStatisticsTask(kTriggerMask);
194
195   // Clean Geometry: Ruben
196   gROOT->LoadMacro("$ALICE_ROOT/PWGPP/CalibMacros/commonMacros/CleanGeom.C++");
197   CleanGeom* clgmTask = new CleanGeom("cleanGeom");
198   mgr->AddTask(clgmTask);
199   AliAnalysisDataContainer *dummyInp = mgr->GetCommonInputContainer();
200   if (dummyInp) mgr->ConnectInput(clgmTask,0,dummyInp);  
201
202   //
203   // CDB connection
204   //
205   if (doCDBconnect) {
206     gROOT->LoadMacro("$ALICE_ROOT/PWGPP/PilotTrain/AddTaskCDBconnect.C");
207     AliTaskCDBconnect *taskCDB = AddTaskCDBconnect(cdb_location, run_number);    
208 //    AliTaskCDBconnect *taskCDB = AddTaskCDBconnect();
209     if (!taskCDB) return;
210 //    AliCDBManager *cdb = AliCDBManager::Instance();
211 //    cdb->SetDefaultStorage(cdb_location);
212 //    taskCDB->SetRunNumber(run_number);
213   }    
214   
215   //
216   // Event Statistics (Jan Fiete)
217   //
218   if (doEventStat) {
219       gROOT->LoadMacro("$ALICE_ROOT/ANALYSIS/macros/AddTaskPhysicsSelection.C");
220       AliPhysicsSelectionTask* physSelTask = AddTaskPhysicsSelection(kFALSE /*MC*/);
221       AliAnalysisDataContainer *cstatsout = (AliAnalysisDataContainer*)mgr->GetOutputs()->FindObject("cstatsout");
222       cstatsout->SetFileName(Form("EventStat_temp%s.root", suffix));
223   }
224   
225    //
226   // PIDResponse(JENS)
227   //
228   if (doPIDResponse && (ibarrel || iall)) {
229     gROOT->LoadMacro("$ALICE_ROOT/ANALYSIS/macros/AddTaskPIDResponse.C"); 
230     AliAnalysisTaskPIDResponse *PIDResponse = AddTaskPIDResponse();
231     PIDResponse->SetUserDataRecoPass(1);
232     PIDResponse->SelectCollisionCandidates(kTriggerMask);
233   }  
234
235
236   //
237   // Centrality (A. Toia)
238   //
239   if (doCentrality) {
240 //     if (!iCollisionType) {
241 //        printf("Disabling centrality task for p-p\n");
242 //        doCentrality = kFALSE;
243 //     } else {           
244         gROOT->LoadMacro("$ALICE_ROOT/ANALYSIS/macros/AddTaskCentrality.C");
245         AliCentralitySelectionTask *taskCentrality = AddTaskCentrality();
246       AliAnalysisDataContainer *centralityout = (AliAnalysisDataContainer*)mgr->GetOutputs()->FindObject("CentralityStat");
247       centralityout->SetFileName(Form("EventStat_temp%s.root", suffix));
248 //     }   
249   }   
250   
251   // Vertexing (A. Dainese)
252   // 
253   if (doVertex && (ibarrel || iall)) {
254     gROOT->LoadMacro("$ALICE_ROOT/PWGPP/macros/AddTaskVertexESD.C");
255     AliAnalysisTaskVertexESD* taskvertexesd =  AddTaskVertexESD(kFALSE, kTriggerMask);
256     taskvertexesd->SelectCollisionCandidates(kTriggerMask);
257   }  
258
259   // TPC QA (E. Sicking)
260   //
261   if (doQAsym && (ibarrel || iall)) {
262   // offline trigger in AddTask
263     gROOT->LoadMacro("$ALICE_ROOT/PWGPP/PilotTrain/AddTaskQAsym.C");
264     AliAnalysisTaskSE * taskqasim = AddTaskQAsym(0, kTriggerMask, kTriggerHM, kTriggerEMC, kTriggerMuonBarrel);
265   }  
266   //
267   // VZERO QA  (C. Cheshkov)
268   //
269   if (doVZERO && (ibarrel || iall)) {
270     gROOT->LoadMacro("$ALICE_ROOT/PWGPP/PilotTrain/AddTaskVZEROQA.C");
271     AliAnalysisTaskSE * taskv0qa = AddTaskVZEROQA(0);
272 //  taskv0qa->SelectCollisionCandidates();
273     gROOT->LoadMacro("$ALICE_ROOT/PWGPP/PilotTrain/AddTaskVZEROQATrig.C");
274     AliAnaVZEROQA *taskv0qatrig = AddTaskVZEROQATrig(0);
275     taskv0qatrig->SelectCollisionCandidates(AliVEvent::kINT8);
276
277   }
278   if (doVZEROPbPb && iCollisionType==0) {
279     gROOT->LoadMacro("$ALICE_ROOT/PWGPP/VZERO/AddTaskVZEROPbPb.C");
280     AliAnaVZEROPbPb* taskV0PbPb = (AliAnaVZEROPbPb*)AddTaskVZEROPbPb(run_number);
281 //    taskV0PbPb->SetClassesNames("CTRUE-,C0HWU-,CPBI2WU_B1-,CPBI2_B1-,CPBI1WU-,CPBI1-,CVHNWU-,CVHN-,CVLNWU-,CVLN-");
282 //    taskV0PbPb->SetClassesNames("CTRUE-,C0HWU-,CPBI2WU,CPBI2,CPBI1WU-,CPBI1-,CVHNWU,CVHN,CVLNWU,CVLN");
283 //    taskV0PbPb->SetClassesNames("CTRUE-,C0HWU-,CPBI2WU-,CPBI2-,CPBI2WU_B1-,CPBI2_B1-,CPBI1WU-,CPBI1-,CVHNWU-,CVHN-,CVHN_R2-,CVHNWU_R2-,CVLNWU-,CVLN-,CVLN_B2-,CVLNWU_B2-");
284 //    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-");
285 //    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-");
286   }
287   //
288   // TPC (Jacek Otwinowski & Michael Knichel)
289   //
290   //
291   // Optionally MC information can be used by setting the 1st argument to true
292   // Optionally friends information can be switched off by setting the 2st argument 
293   // to false
294   // Optionally highMult axis can be used by setting the 3st argument to true (for PbPb)
295   if (doTPC && (ibarrel || iall)) {
296     gROOT->LoadMacro("$ALICE_ROOT/PWGPP/TPC/macros/AddTaskPerformanceTPCdEdxQA.C");
297     AliPerformanceTask *tpcQA = 0;
298     if (iCollisionType) {
299        // High multiplicity Pb-Pb
300        tpcQA = AddTaskPerformanceTPCdEdxQA(kFALSE, kTRUE, kTRUE);
301     } else {
302       // Low multiplicity (pp)
303        tpcQA = AddTaskPerformanceTPCdEdxQA(kFALSE, kTRUE, kFALSE);
304     }
305     tpcQA->SelectCollisionCandidates(kTriggerMask);
306   }  
307
308   // HLT (Alberica Toia)
309   if (doHLT && (ibarrel || iall)) {
310     gROOT->LoadMacro("$ALICE_ROOT/PWGPP/TPC/macros/AddTaskPerformanceTPCdEdxQA.C");
311     AliPerformanceTask *hltQA = AddTaskPerformanceTPCdEdxQA(kFALSE, kTRUE, kFALSE,0,kTRUE);
312     hltQA->SelectCollisionCandidates(kTriggerMask);
313   }  
314   //
315   // SPD (A. Mastroserio)
316   //
317   if (doSPD && (ibarrel || iall)) {
318     gROOT->LoadMacro("$ALICE_ROOT/PWGPP/PilotTrain/AddTaskSPDQA.C");
319     AliAnalysisTaskSPD* taskspdqa = (AliAnalysisTaskSPD*)AddTaskSPDQA();
320     // Request from Annalisa
321     if (iCollisionType) taskspdqa->SetHeavyIonMode();
322     taskspdqa->SelectCollisionCandidates(kTriggerMask);
323     taskspdqa->SetOCDBInfo(run_number, cdb_location);
324   }  
325   //
326   // SDD (F. Prino)
327   //
328   if (doSDD && (ibarrel || iall)) {
329     gROOT->LoadMacro("$ALICE_ROOT/PWGPP/PilotTrain/AddSDDPoints.C");
330     AliAnalysisTaskSE* tasksdd = AddSDDPoints();
331     tasksdd->SelectCollisionCandidates(kTriggerMask);
332   }
333   //
334   // SSD dEdx (Marek Chojnacki)
335   //
336   if (doSSDdEdx && (ibarrel || iall)) {
337     gROOT->LoadMacro("$ALICE_ROOT/PWGPP/PilotTrain/AddTaskdEdxSSDQA.C");
338     AliAnalysisTaskSE* taskssddedx = AddTaskdEdxSSDQA();
339     taskssddedx->SelectCollisionCandidates(kTriggerMask);
340   }
341
342   //
343   // ITS
344   //
345   if (doITS && (ibarrel || iall)) {
346   // hardcoded non-zero trigger mask
347       gROOT->LoadMacro("$ALICE_ROOT/PWGPP/macros/AddTaskPerformanceITS.C");
348       AliAnalysisTaskITSTrackingCheck *itsQA = 0;
349       AliAnalysisTaskITSTrackingCheck *itsQACent0010 = 0;
350       AliAnalysisTaskITSTrackingCheck *itsQACent3050 = 0;
351       AliAnalysisTaskITSTrackingCheck *itsQACent6080 = 0;
352       if(iCollisionType==0) {
353         itsQA = AddTaskPerformanceITS(kFALSE);
354       } else {
355         itsQA = AddTaskPerformanceITS(kFALSE);
356         itsQACent0010 = AddTaskPerformanceITS(kFALSE,kFALSE,kFALSE,3500,10000);
357         itsQACent3050 = AddTaskPerformanceITS(kFALSE,kFALSE,kFALSE,590,1570);
358         itsQACent6080 = AddTaskPerformanceITS(kFALSE,kFALSE,kFALSE,70,310);
359       }
360   }
361   //
362   // ITS saTracks, align (F.Prino)
363   //
364   if (doITSsaTracks && (ibarrel || iall)) {
365   // offline trigger in AddTask
366      gROOT->LoadMacro("$ALICE_ROOT/PWGPP/macros/AddTaskITSsaTracks.C");
367      AliAnalysisTaskITSsaTracks *itssaTracks = AddTaskITSsaTracks(kFALSE,kFALSE);
368      itssaTracks->SelectCollisionCandidates(kTriggerMask);
369   }   
370   if (doITSalign && (ibarrel || iall)) {
371   // no offline trigger selection
372      gROOT->LoadMacro("$ALICE_ROOT/PWGPP/macros/AddTaskITSAlign.C");
373      AliAnalysisTaskITSAlignQA *itsAlign = AddTaskITSAlign(0,2012);
374   }   
375   //
376   // TRD (Alex Bercuci, M. Fasel) 
377   //
378   if(doTRD && (ibarrel || iall)) {
379   // no offline trigger selection
380       gROOT->LoadMacro("$ALICE_ROOT/PWGPP/macros/AddTrainPerformanceTRD.C");
381       // steer individual TRD tasks
382       Bool_t 
383       doCheckESD(kTRUE),  // AliTRDcheckESD
384       doCheckDET(kTRUE),  // AliTRDcheckDET
385       doEffic(kTRUE),     // AliTRDefficiency
386       doResolution(kTRUE),// AliTRDresolution
387       doCheckPID(kTRUE),  // AliTRDcheckPID
388       doV0Monitor(kFALSE);// AliTRDv0Monitor
389       AddTrainPerformanceTRD(Translate(doCheckESD, doCheckDET, doEffic, doResolution, doCheckPID, doV0Monitor));
390   }
391
392   //
393   // ZDC (Chiara Oppedisano) 
394   //
395   if(doZDC && (ibarrel || iall)) {
396   // hardcoded kMB trigger mask
397      gROOT->LoadMacro("$ALICE_ROOT/PWGPP/ZDC/AddTaskZDCQA.C");
398      AliAnalysisTaskSE *taskZDC = AddTaskZDCQA();
399      taskZDC->SelectCollisionCandidates(kTriggerMask);
400   }   
401   //
402   // Calorimetry (Gustavo Conesa)
403   //
404
405   if(doCALO) {
406   // - In case of MC analysis, it skips the triggered events wagon (given that
407   //the normal wagon is always suffixed as "default")
408   //- No need to pass the type of data
409   //- No need to specify the year. In principle the year is only needed when
410   //setting the size of some histograms (change in number of SM along the years),
411   //I do not know how to access automatically the run number when the histograms
412   // are created in the UserCreateOutput. By default I set them to the maximum
413   //expected size, but I still let the possibility to check the year.
414
415   //So the way to initialize the task in trunk is now
416
417   //AddTaskCalorimeterQA(Bool_t kSimulation = kFALSE,
418                                    //  const char *suffix="default",
419                                    // TString outputFile = "", 
420                                    //  Int_t year = 2012, 
421                                    //  Bool_t kPrintSettings = kFALSE)
422
423   //So, in principle only the first 2 need to be specified.
424
425   
426       gROOT->LoadMacro("$ALICE_ROOT/PWGGA/CaloTrackCorrelations/macros/QA/AddTaskCalorimeterQA.C");
427       AliAnalysisTaskCaloTrackCorrelation *taskCaloQA = AddTaskCalorimeterQA("default");
428       taskCaloQA->SetDebugLevel(0);
429       // offline mask set in AddTask to kMB
430       taskCaloQA->SelectCollisionCandidates(kTriggerMask);
431       // Add a new calo task with EMC1 trigger only
432       // taskCaloQA = AddTaskCalorimeterQA("trigEMC");        // for 2011 data, not 2010
433       // taskCaloQA->SelectCollisionCandidates(kTriggerEMC);  // for 2011 data, not 2010
434       taskCaloQA->SetDebugLevel(0);
435   }
436
437   //
438   // Muon Trigger
439   //
440   
441   if(doMUONTrig) {
442   // no offline trigger selection
443       gROOT->LoadMacro("$ALICE_ROOT/PWGPP/macros/AddTaskMTRchamberEfficiency.C");
444       AliAnalysisTaskTrigChEff *taskMuonTrig = AddTaskMTRchamberEfficiency();
445 //      taskMuonTrig->GetMuonTrackCuts()->SetCustomParamFromRun(160000,2);       
446   }
447
448   //
449   // Impact parameter resolution (xianbao.yuan@pd.infn.it, andrea.dainese@pd.infn.it)
450   //
451   if (doImpParRes && (ibarrel || iall)) {
452     gROOT->LoadMacro("$ALICE_ROOT/PWGPP/macros/AddTaskImpParRes.C");
453     AliAnalysisTaskSE* taskimpparres=0;
454     if(iCollisionType==0) {
455        taskimpparres= AddTaskImpParRes();
456     } else {
457        taskimpparres= AddTaskImpParRes(kFALSE,-1,kFALSE,kFALSE);
458     }
459     taskimpparres->SelectCollisionCandidates(kTriggerMask);
460   }  
461   //
462   // MUON QA (Philippe Pillot)
463   //
464   if (doMUON) {
465   // trigger analysis internal
466     gROOT->LoadMacro("$ALICE_ROOT/PWGPP/PilotTrain/AddTaskMuonQA.C");
467     AliAnalysisTaskMuonQA* taskmuonqa= AddTaskMuonQA();
468 //    taskmuonqa->GetTrackCuts()->SetCustomParamFromRun(160000,2);
469   }  
470   //
471   // TOF (Francesca Bellini)
472   //
473   if (doTOF && (ibarrel || iall)) {
474     gROOT->LoadMacro("$ALICE_ROOT/PWGPP/TOF/AddTaskTOFQA.C");
475     AliAnalysisTaskTOFqa *tofQA = AddTaskTOFQA(kFALSE);
476     tofQA->SelectCollisionCandidates(kTriggerMask);
477   } 
478   //
479   // PIDqa(JENS)
480   //
481   if (doPIDqa && (ibarrel || iall)) {
482     gROOT->LoadMacro("$ALICE_ROOT/ANALYSIS/macros/AddTaskPIDqa.C");
483     AliAnalysisTaskPIDqa *PIDQA = AddTaskPIDqa();
484     PIDQA->SelectCollisionCandidates(kTriggerMask);
485   }  
486  
487   //
488   // HMPID QA (Giacomo Volpe)
489   //
490   if (doHMPID && (ibarrel || iall)) {
491     gROOT->LoadMacro("$ALICE_ROOT/PWGPP/HMPID/AddTaskHmpidQA.C");
492     AliAnalysisTaskSE* taskhmpidqa= AddTaskHmpidQA(kFALSE);
493       // offline mask set in AddTask to kMB
494     taskhmpidqa->SelectCollisionCandidates(kTriggerMask);
495   }      
496   
497   
498   // T0 QA (Alla Mayevskaya)
499   if (doT0 && (ibarrel || iall)) {
500   // no offline trigger selection
501     gROOT->LoadMacro("$ALICE_ROOT/PWGPP/T0/AddTaskT0QA.C");
502     AliT0AnalysisTaskQA* taskt0qa= AddTaskT0QA();
503     taskt0qa->SelectCollisionCandidates(kTriggerMask);
504   }      
505   // FMD QA (Christian Holm Christiansen)
506   if (doFMD && (ibarrel || iall)) {
507     gROOT->LoadMacro("$ALICE_ROOT/PWGLF/FORWARD/analysis2/AddTaskForwardQA.C");
508     // Parameters: usemc, usecentrality
509     AliForwardQATask *forwardQA = (AliAnalysisTaskSE *)AddTaskForwardQA(kFALSE, kFALSE);
510     forwardQA->SetDebug(0);
511     // No offline trigger config. needed (see #84077)
512   }
513    //     
514   // PHOS QA (Boris Polishchuk)
515   //
516   if (doPHOS) {
517     gROOT->LoadMacro("$ALICE_ROOT/PWGGA/PHOSTasks/CaloCellQA/macros/AddTaskCaloCellsQA.C");
518     AliAnalysisTaskCaloCellsQA *taskPHOSCellQA1 = AddTaskCaloCellsQA(4, 1, NULL,"PHOSCellsQA_AnyInt"); 
519     taskPHOSCellQA1->SelectCollisionCandidates(kTriggerMask);
520     taskPHOSCellQA1->GetCaloCellsQA()->SetClusterEnergyCuts(0.3,0.3,1.0);
521     AliAnalysisTaskCaloCellsQA *taskPHOSCellQA2 = AddTaskCaloCellsQA(4, 1, NULL,"PHOSCellsQA_PHI7"); 
522     taskPHOSCellQA2->SelectCollisionCandidates(AliVEvent::kPHI7);
523     taskPHOSCellQA2->GetCaloCellsQA()->SetClusterEnergyCuts(0.3,0.3,1.0);
524     // Pi0 QA fo PbPb
525     if (iCollisionType == 0) {
526       gROOT->LoadMacro("$ALICE_ROOT/PWGGA/PHOSTasks/PHOS_PbPbQA/macros/AddTaskPHOSPbPb.C");
527       AliAnalysisTaskPHOSPbPbQA* phosPbPb = AddTaskPHOSPbPbQA(0);
528     }
529   }    
530   if (doPHOSTrig) {
531      gROOT->LoadMacro("$ALICE_ROOT/PWGGA/PHOSTasks/PHOS_TriggerQA/macros/AddTaskPHOSTriggerQA.C");
532      AliAnalysisTaskPHOSTriggerQA *taskPHOSTrig = AddTaskPHOSTriggerQA(NULL);
533   }
534   //
535   // EMCAL QA (Gustavo Conesa)
536   //
537   if (doEMCAL) {
538      gROOT->LoadMacro("$ALICE_ROOT/PWGGA/EMCALTasks/macros/AddTaskEMCALTriggerQA.C");
539      AliAnalysisTaskEMCALTriggerQA *emctrig = AddTaskEMCALTriggerQA();
540 //     emctrig->GetRecoUtils()->SwitchOffBadChannelsRemoval();
541   }   
542   //     
543   // FLOW and BF QA (C.Perez && A.Rodriguez)
544   //
545   if (doFBFqa && (ibarrel || iall)) {
546     gROOT->LoadMacro("$ALICE_ROOT/PWGPP/macros/AddTaskFBFqa.C");
547     AliAnalysisTaskSE *qaFBFMB = (AliAnalysisTaskSE*) AddTaskFBFqa("qaFBFmb",kFALSE);
548     qaFBFMB->SelectCollisionCandidates(AliVEvent::kMB);
549     AliAnalysisTaskSE *qaFBFSC = (AliAnalysisTaskSE*) AddTaskFBFqa("qaFBFsc",kFALSE);
550     qaFBFSC->SelectCollisionCandidates(AliVEvent::kSemiCentral);
551     AliAnalysisTaskSE *qaFBFCE = (AliAnalysisTaskSE*) AddTaskFBFqa("qaFBFce",kFALSE);
552     qaFBFCE->SelectCollisionCandidates(AliVEvent::kCentral);
553   }
554   
555 //Jacek
556    if (gSystem->Getenv("QA_TaskFilteredTree")) doTaskFilteredTree=1;
557    if (doTaskFilteredTree) {
558       gROOT->LoadMacro("$ALICE_ROOT/PWGPP/macros/AddTaskFilteredTree.C");
559       AddTaskFilteredTree("FilterEvents_Trees.root");
560    }   
561    
562 }
563
564 void QAmerge(const char *suffix, const char *dir, Int_t stage)
565 {
566 // Merging method
567   TStopwatch timer;
568   timer.Start();
569   TString outputDir = dir;
570   TString outputFiles = Form("QAresults%s.root,EventStat_temp%s.root,RecoQAresults%s.root,FilterEvents_trees%s.root",suffix,suffix,suffix);
571   TString mergeExcludes = "";
572   TObjArray *list = outputFiles.Tokenize(",");
573   TIter *iter = new TIter(list);
574   TObjString *str;
575   TString outputFile;
576   Bool_t merged = kTRUE;
577   while((str=(TObjString*)iter->Next())) {
578     outputFile = str->GetString();
579     // Skip already merged outputs
580     if (!gSystem->AccessPathName(outputFile)) {
581        printf("Output file <%s> found. Not merging again.",outputFile.Data());
582        continue;
583     }
584     if (mergeExcludes.Contains(outputFile.Data())) continue;
585     merged = AliAnalysisAlien::MergeOutput(outputFile, outputDir, 10, stage);
586     if (!merged && !outputFile.Contains("RecoQAresults") && !outputFile.Contains("FilterEvents_trees")) {
587        printf("ERROR: Cannot merge %s\n", outputFile.Data());
588        return;
589     }
590   }
591   // read the analysis manager from file
592   TString infolog = "fileinfo.log";
593   AliAnalysisAlien::MergeInfo(infolog, dir); 
594   if (!outputDir.Contains("Stage")) {
595     ofstream out;
596     out.open("outputs_valid", ios::out);
597     out.close();    
598     return;
599   }
600   AliAnalysisManager *mgr = AliAnalysisManager::GetAnalysisManager();
601   mgr->SetRunFromPath(mgr->GetRunFromAlienPath(dir));
602   mgr->SetSkipTerminate(kFALSE);
603   if (!mgr->InitAnalysis()) return;
604   mgr->PrintStatus();
605   AliLog::SetGlobalLogLevel(AliLog::kError);
606   TTree *tree = NULL;
607   gROOT->cd();
608   mgr->StartAnalysis("gridterminate", tree);
609   if (strlen(suffix)) {
610      if (gSystem->Exec(Form("mv trending.root trending%s.root", suffix)))
611         ::Error("QAmerge", "File trending.root was not produced");
612      if (gSystem->Exec(Form("mv event_stat.root event_stat%s.root", suffix)))
613         ::Error("QAmerge", "File trending.root was not produced");
614   }   
615   ofstream out;
616   out.open("outputs_valid", ios::out);
617   out.close();
618   timer.Print();
619 }
620
621 void ProcessEnvironmentVars()
622 {
623   int var=0;
624   TString envS;
625   //
626   if ( !(envS=gSystem->Getenv("CDBconnect")).IsNull() && CheckEnvS("CDBconnect",envS) ) {
627     doCDBconnect = atoi(envS.Data());
628     printf("Set doCDBconnect=%d according to environment variables\n",doCDBconnect);
629   }
630   //
631   if ( !(envS=gSystem->Getenv("Statistics")).IsNull() && CheckEnvS("Statistics",envS) ) {
632     doStatistics = atoi(envS.Data());
633     printf("Set doStatistics=%d according to environment variables\n",doStatistics);
634   }
635   //
636   if ( !(envS=gSystem->Getenv("EventStat")).IsNull() && CheckEnvS("EventStat",envS) ) {
637     doEventStat = atoi(envS.Data());
638     printf("Set doEventStat=%d according to environment variables\n",doEventStat);
639   }
640   //
641   if ( !(envS=gSystem->Getenv("Centrality")).IsNull() && CheckEnvS("Centrality",envS) ) {
642     doCentrality = atoi(envS.Data());
643     printf("Set doCentrality=%d according to environment variables\n",doCentrality);
644   }
645   //
646   if ( !(envS=gSystem->Getenv("QAsym")).IsNull() && CheckEnvS("QAsym",envS) ) {
647     doQAsym = atoi(envS.Data());
648     printf("Set doQAsym=%d according to environment variables\n",doQAsym);
649   }
650   //
651   if ( !(envS=gSystem->Getenv("VZERO")).IsNull() && CheckEnvS("VZERO",envS) ) {
652     doVZERO = atoi(envS.Data());
653     printf("Set doVZERO=%d according to environment variables\n",doVZERO);
654   }
655   //
656   if ( !(envS=gSystem->Getenv("VZEROPbPb")).IsNull() && CheckEnvS("VZEROPbPb",envS) ) {
657     doVZEROPbPb = atoi(envS.Data());
658     printf("Set doVZEROPbPb=%d according to environment variables\n",doVZEROPbPb);
659   }
660   //
661   if ( !(envS=gSystem->Getenv("Vertex")).IsNull() && CheckEnvS("Vertex",envS) ) {
662     doVertex = atoi(envS.Data());
663     printf("Set doVertex=%d according to environment variables\n",doVertex);
664   }
665   //
666   if ( !(envS=gSystem->Getenv("SPD")).IsNull() && CheckEnvS("SPD",envS) ) {
667     doSPD = atoi(envS.Data());
668     printf("Set doSPD=%d according to environment variables\n",doSPD);
669   }
670   //
671   if ( !(envS=gSystem->Getenv("TPC")).IsNull() && CheckEnvS("TPC",envS) ) {
672     doTPC = atoi(envS.Data());
673     printf("Set doTPC=%d according to environment variables\n",doTPC);
674   }
675   //
676   if ( !(envS=gSystem->Getenv("HLT")).IsNull() && CheckEnvS("HLT",envS) ) {
677     doHLT = atoi(envS.Data());
678     printf("Set doHLT=%d according to environment variables\n",doHLT);
679   }
680   //
681   if ( !(envS=gSystem->Getenv("SDD")).IsNull() && CheckEnvS("SDD",envS) ) {
682     doSDD = atoi(envS.Data());
683     printf("Set doSDD=%d according to environment variables\n",doSDD);
684   }
685   //
686   if ( !(envS=gSystem->Getenv("SSDdEdx")).IsNull() && CheckEnvS("SSDdEdx",envS) ) {
687     doSSDdEdx = atoi(envS.Data());
688     printf("Set doSSDdEdx=%d according to environment variables\n",doSSDdEdx);
689   }
690   //
691   if ( !(envS=gSystem->Getenv("TRD")).IsNull() && CheckEnvS("TRD",envS) ) {
692     doTRD = atoi(envS.Data());
693     printf("Set doTRD=%d according to environment variables\n",doTRD);
694   }
695   //
696   if ( !(envS=gSystem->Getenv("ITS")).IsNull() && CheckEnvS("ITS",envS) ) {
697     doITS = atoi(envS.Data());
698     printf("Set doITS=%d according to environment variables\n",doITS);
699   }
700   //
701   if ( !(envS=gSystem->Getenv("ITSsaTracks")).IsNull() && CheckEnvS("ITSsaTracks",envS) ) {
702     doITSsaTracks = atoi(envS.Data());
703     printf("Set doITSsaTracks=%d according to environment variables\n",doITSsaTracks);
704   }
705   //
706   if ( !(envS=gSystem->Getenv("ITSalign")).IsNull() && CheckEnvS("ITSalign",envS) ) {
707     doITSalign = atoi(envS.Data());
708     printf("Set doITSalign=%d according to environment variables\n",doITSalign);
709   }
710   //
711   if ( !(envS=gSystem->Getenv("CALO")).IsNull() && CheckEnvS("CALO",envS) ) {
712     doCALO = atoi(envS.Data());
713     printf("Set doCALO=%d according to environment variables\n",doCALO);
714   }
715   //
716   if ( !(envS=gSystem->Getenv("MUONTrig")).IsNull() && CheckEnvS("MUONTrig",envS) ) {
717     doMUONTrig = atoi(envS.Data());
718     printf("Set doMUONTrig=%d according to environment variables\n",doMUONTrig);
719   }
720   //
721   if ( !(envS=gSystem->Getenv("ImpParRes")).IsNull() && CheckEnvS("ImpParRes",envS) ) {
722     doImpParRes = atoi(envS.Data());
723     printf("Set doImpParRes=%d according to environment variables\n",doImpParRes);
724   }
725   //
726   if ( !(envS=gSystem->Getenv("MUON")).IsNull() && CheckEnvS("MUON",envS) ) {
727     doMUON = atoi(envS.Data());
728     printf("Set doMUON=%d according to environment variables\n",doMUON);
729   }
730   //
731   if ( !(envS=gSystem->Getenv("TOF")).IsNull() && CheckEnvS("TOF",envS) ) {
732     doTOF = atoi(envS.Data());
733     printf("Set doTOF=%d according to environment variables\n",doTOF);
734   }
735   //
736   if ( !(envS=gSystem->Getenv("HMPID")).IsNull() && CheckEnvS("HMPID",envS) ) {
737     doHMPID = atoi(envS.Data());
738     printf("Set doHMPID=%d according to environment variables\n",doHMPID);
739   }
740   //
741   if ( !(envS=gSystem->Getenv("T0")).IsNull() && CheckEnvS("T0",envS) ) {
742     doT0 = atoi(envS.Data());
743     printf("Set doT0=%d according to environment variables\n",doT0);
744   }
745   //
746   if ( !(envS=gSystem->Getenv("ZDC")).IsNull() && CheckEnvS("ZDC",envS) ) {
747     doZDC = atoi(envS.Data());
748     printf("Set doZDC=%d according to environment variables\n",doZDC);
749   }
750   //
751   if ( !(envS=gSystem->Getenv("PIDResponse")).IsNull() && CheckEnvS("PIDResponse",envS) ) {
752     doPIDResponse = atoi(envS.Data());
753     printf("Set doPIDResponse=%d according to environment variables\n",doPIDResponse);
754   }
755   //
756   if ( !(envS=gSystem->Getenv("PIDqa")).IsNull() && CheckEnvS("PIDqa",envS) ) {
757     doPIDqa = atoi(envS.Data());
758     printf("Set doPIDqa=%d according to environment variables\n",doPIDqa);
759   }
760   //
761   if ( !(envS=gSystem->Getenv("FMD")).IsNull() && CheckEnvS("FMD",envS) ) {
762     doFMD = atoi(envS.Data());
763     printf("Set doFMD=%d according to environment variables\n",doFMD);
764   }
765   //
766   if ( !(envS=gSystem->Getenv("PHOS")).IsNull() && CheckEnvS("PHOS",envS) ) {
767     doPHOS = atoi(envS.Data());
768     printf("Set doPHOS=%d according to environment variables\n",doPHOS);
769   }
770   //
771   if ( !(envS=gSystem->Getenv("PHOSTrig")).IsNull() && CheckEnvS("PHOSTrig",envS) ) {
772     doPHOSTrig = atoi(envS.Data());
773     printf("Set doPHOSTrig=%d according to environment variables\n",doPHOSTrig);
774   }
775   //
776   if ( !(envS=gSystem->Getenv("EMCAL")).IsNull() && CheckEnvS("EMCAL",envS) ) {
777     doEMCAL = atoi(envS.Data());
778     printf("Set doEMCAL=%d according to environment variables\n",doEMCAL);
779   }
780   //
781   if ( !(envS=gSystem->Getenv("FBFqa")).IsNull() && CheckEnvS("FBFqa",envS) ) {
782     doFBFqa = atoi(envS.Data());
783     printf("Set doFBFqa=%d according to environment variables\n",doFBFqa);
784   }
785   //
786   if ( !(envS=gSystem->Getenv("TaskFilteredTree")).IsNull() && CheckEnvS("TaskFilteredTree",envS) ) {
787     doTaskFilteredTree = atoi(envS.Data());
788     printf("Set doTaskFilteredTree=%d according to environment variables\n",doTaskFilteredTree);
789   }
790   //-----------------------------------------------------------  
791   // combinations
792   if ( !(envS=gSystem->Getenv("QAGROUP")).IsNull() ) {
793     int qaGRP = atoi(envS.Data());
794     switch (qaGRP) {
795     case 0: 
796       SetDoQA(doStatistics,  1, "Statistics");
797       SetDoQA(doVZERO,       1, "VZERO");
798       SetDoQA(doVZEROPbPb,   1, "VZEROPbPb");
799       SetDoQA(doVertex,      1, "Vertex");
800       SetDoQA(doSPD,         1, "SPD");
801       SetDoQA(doTPC,         0, "TPC");
802       SetDoQA(doHLT,         0, "HLT");
803       SetDoQA(doSDD,         0, "SDD");
804       SetDoQA(doSSDdEdx,     0, "SSDdEdx");
805       SetDoQA(doTRD,         1, "TRD");
806       SetDoQA(doITS,         0, "ITS");
807       SetDoQA(doITSsaTracks, 0, "ITSsaTracks");
808       SetDoQA(doITSalign,    0, "ITSalign");
809       SetDoQA(doCALO ,       1, "CALO");
810       SetDoQA(doMUONTrig,    0, "MUONTrig");
811       SetDoQA(doImpParRes,   1, "ImpParRes");
812       SetDoQA(doMUON,        0, "MUON");
813       SetDoQA(doTOF,         0, "TOF");
814       SetDoQA(doHMPID,       1, "HMPID");
815       SetDoQA(doT0,          1, "T0");
816       SetDoQA(doZDC,         1, "ZDC");
817       //      SetDoQA(doPIDResponse, 1, "PIDResponse");
818       SetDoQA(doPIDqa,       1, "PIDqa");
819       SetDoQA(doFMD,         1, "FMD");
820       SetDoQA(doPHOS,        1, "PHOS");
821       SetDoQA(doPHOSTrig,    1, "PHOSTrig");
822       SetDoQA(doEMCAL,       0, "EMCAL");
823       SetDoQA(doFBFqa,       1, "FBFqa");
824       break;
825     case 1: 
826       SetDoQA(doStatistics,  0, "Statistics");
827       SetDoQA(doVZERO,       0, "VZERO");
828       SetDoQA(doVZEROPbPb,   0, "VZEROPbPb");
829       SetDoQA(doVertex,      0, "Vertex");
830       SetDoQA(doSPD,         0, "SPD");
831       SetDoQA(doTPC,         1, "TPC");
832       SetDoQA(doHLT,         1, "HLT");
833       SetDoQA(doSDD,         0, "SDD");
834       SetDoQA(doSSDdEdx,     0, "SSDdEdx");
835       SetDoQA(doTRD,         0, "TRD");
836       SetDoQA(doITS,         0, "ITS");
837       SetDoQA(doITSsaTracks, 0, "ITSsaTracks");
838       SetDoQA(doITSalign,    0, "ITSalign");
839       SetDoQA(doCALO ,       0, "CALO");
840       SetDoQA(doMUONTrig,    0, "MUONTrig");
841       SetDoQA(doImpParRes,   0, "ImpParRes");
842       SetDoQA(doMUON,        0, "MUON");
843       SetDoQA(doTOF,         0, "TOF");
844       SetDoQA(doHMPID,       0, "HMPID");
845       SetDoQA(doT0,          0, "T0");
846       SetDoQA(doZDC,         0, "ZDC");
847       //      SetDoQA(doPIDResponse, 0, "PIDResponse");
848       SetDoQA(doPIDqa,       0, "PIDqa");
849       SetDoQA(doFMD,         0, "FMD");
850       SetDoQA(doPHOS,        0, "PHOS");
851       SetDoQA(doPHOSTrig,    0, "PHOSTrig");
852       SetDoQA(doEMCAL,       0, "EMCAL");
853       SetDoQA(doFBFqa,       0, "FBFqa");
854       break;
855     case 2: 
856       SetDoQA(doStatistics,  0, "Statistics");
857       SetDoQA(doVZERO,       0, "VZERO");
858       SetDoQA(doVZEROPbPb,   0, "VZEROPbPb");
859       SetDoQA(doVertex,      0, "Vertex");
860       SetDoQA(doSPD,         0, "SPD");
861       SetDoQA(doTPC,         0, "TPC");
862       SetDoQA(doHLT,         0, "HLT");
863       SetDoQA(doSDD,         0, "SDD");
864       SetDoQA(doSSDdEdx,     0, "SSDdEdx");
865       SetDoQA(doTRD,         0, "TRD");
866       SetDoQA(doITS,         0, "ITS");
867       SetDoQA(doITSsaTracks, 1, "ITSsaTracks");
868       SetDoQA(doITSalign,    1, "ITSalign");
869       SetDoQA(doCALO ,       0, "CALO");
870       SetDoQA(doMUONTrig,    0, "MUONTrig");
871       SetDoQA(doImpParRes,   0, "ImpParRes");
872       SetDoQA(doMUON,        0, "MUON");
873       SetDoQA(doTOF,         0, "TOF");
874       SetDoQA(doHMPID,       0, "HMPID");
875       SetDoQA(doT0,          0, "T0");
876       SetDoQA(doZDC,         0, "ZDC");
877       //      SetDoQA(doPIDResponse, 0, "PIDResponse");
878       SetDoQA(doPIDqa,       0, "PIDqa");
879       SetDoQA(doFMD,         0, "FMD");
880       SetDoQA(doPHOS,        0, "PHOS");
881       SetDoQA(doPHOSTrig,    0, "PHOSTrig");
882       SetDoQA(doEMCAL,       0, "EMCAL");
883       SetDoQA(doFBFqa,       0, "FBFqa");
884       break;
885     case 3: 
886       SetDoQA(doStatistics,  0, "Statistics");
887       SetDoQA(doVZERO,       0, "VZERO");
888       SetDoQA(doVZEROPbPb,   0, "VZEROPbPb");
889       SetDoQA(doVertex,      0, "Vertex");
890       SetDoQA(doSPD,         0, "SPD");
891       SetDoQA(doTPC,         0, "TPC");
892       SetDoQA(doHLT,         0, "HLT");
893       SetDoQA(doSDD,         1, "SDD");
894       SetDoQA(doSSDdEdx,     1, "SSDdEdx");
895       SetDoQA(doTRD,         0, "TRD");
896       SetDoQA(doITS,         1, "ITS");
897       SetDoQA(doITSsaTracks, 0, "ITSsaTracks");
898       SetDoQA(doITSalign,    0, "ITSalign");
899       SetDoQA(doCALO ,       0, "CALO");
900       SetDoQA(doMUONTrig,    0, "MUONTrig");
901       SetDoQA(doImpParRes,   0, "ImpParRes");
902       SetDoQA(doMUON,        0, "MUON");
903       SetDoQA(doTOF,         1, "TOF");
904       SetDoQA(doHMPID,       0, "HMPID");
905       SetDoQA(doT0,          0, "T0");
906       SetDoQA(doZDC,         0, "ZDC");
907       //      SetDoQA(doPIDResponse, 0, "PIDResponse");
908       SetDoQA(doPIDqa,       0, "PIDqa");
909       SetDoQA(doFMD,         0, "FMD");
910       SetDoQA(doPHOS,        0, "PHOS");
911       SetDoQA(doPHOSTrig,    0, "PHOSTrig");
912       SetDoQA(doEMCAL,       0, "EMCAL");
913       SetDoQA(doFBFqa,       0, "FBFqa");
914       break;
915     case 4: 
916       SetDoQA(doStatistics,  0, "Statistics");
917       SetDoQA(doVZERO,       0, "VZERO");
918       SetDoQA(doVZEROPbPb,   0, "VZEROPbPb");
919       SetDoQA(doVertex,      0, "Vertex");
920       SetDoQA(doSPD,         0, "SPD");
921       SetDoQA(doTPC,         0, "TPC");
922       SetDoQA(doHLT,         0, "HLT");
923       SetDoQA(doSDD,         0, "SDD");
924       SetDoQA(doSSDdEdx,     0, "SSDdEdx");
925       SetDoQA(doTRD,         0, "TRD");
926       SetDoQA(doITS,         0, "ITS");
927       SetDoQA(doITSsaTracks, 0, "ITSsaTracks");
928       SetDoQA(doITSalign,    0, "ITSalign");
929       SetDoQA(doCALO ,       0, "CALO");
930       SetDoQA(doMUONTrig,    1, "MUONTrig");
931       SetDoQA(doImpParRes,   0, "ImpParRes");
932       SetDoQA(doMUON,        1, "MUON");
933       SetDoQA(doTOF,         0, "TOF");
934       SetDoQA(doHMPID,       0, "HMPID");
935       SetDoQA(doT0,          0, "T0");
936       SetDoQA(doZDC,         0, "ZDC");
937       //      SetDoQA(doPIDResponse, 0, "PIDResponse");
938       SetDoQA(doPIDqa,       0, "PIDqa");
939       SetDoQA(doFMD,         0, "FMD");
940       SetDoQA(doPHOS,        0, "PHOS");
941       SetDoQA(doPHOSTrig,    0, "PHOSTrig");
942       SetDoQA(doEMCAL,       1, "EMCAL");
943       SetDoQA(doFBFqa,       0, "FBFqa");
944       break;
945     };
946
947     printf("Set doTaskFilteredTree=%d according to environment variables\n",doTaskFilteredTree);
948   }
949
950
951 }
952
953 void SetDoQA(Bool_t& var, Bool_t onoff, const char* name)
954 {
955   // set to ON only if previous value is not off
956   if (var) var = onoff;
957   if (name) printf("Set %-15s to %s\n",name,var ? "ON":"OFF");
958 }
959
960 Bool_t CheckEnvS(const char* var,TString& envString)
961 {
962   if (envString=="0" || envString=="1") return kTRUE;
963   else printf("Ignoring wrong value %s for environment variable %s\n",envString.Data(),var);
964   return kFALSE;
965 }