1 // Steer TRD QA train for Reconstruction (Clusterizer, Tracking and PID).
4 // AddTrainPerformanceTRD.C(MC, friends, tasks)
5 // tasks : "ALL" or one/more of the following:
6 // "EFF" : TRD Tracking Efficiency
7 // "EFFC" : TRD Tracking Efficiency Combined (barrel + stand alone) - only in case of simulations
8 // "MULT" : TRD single track selection
9 // "RES" : TRD tracking Resolution
10 // "CLRES": clusters Resolution
11 // "CAL" : TRD calibration
12 // "ALGN" : TRD alignment
13 // "PID" : TRD PID - pion efficiency
14 // "PIDR" : TRD PID - reference data
15 // "DET" : Basic TRD Detector checks
16 // "NOFR" : Data set does not have AliESDfriends.root
17 // "NOMC" : Data set does not have Monte Carlo Informations (real data), so all tasks which rely
18 // on MC information are switched off
21 // Don't forget to load first the libraries
22 // gSystem->Load("libMemStat.so")
23 // gSystem->Load("libMemStatGui.so")
24 // gSystem->Load("libANALYSIS.so")
25 // gSystem->Load("libANALYSISalice.so")
26 // gSystem->Load("libPWG1.so");
29 // Alex Bercuci (A.Bercuci@gsi.de)
30 // Markus Fasel (m.Fasel@gsi.de)
32 #if ! defined (__CINT__) || defined (__MAKECINT__)
34 #include <Riostream.h>
36 #include "TStopwatch.h"
38 #include "TMemStatViewerGUI.h"
47 #include "TAlienCollection.h"
48 #include "TGridCollection.h"
49 #include "TGridResult.h"
50 #include "TGeoGlobalMagField.h"
53 #include "AliTracker.h"
55 #include "AliCDBManager.h"
56 #include "AliGRPManager.h"
57 #include "AliGeomManager.h"
58 #include "AliAnalysisManager.h"
59 #include "AliAnalysisDataContainer.h"
60 #include "AliMCEventHandler.h"
61 #include "AliESDInputHandler.h"
63 #include "TRD/AliTRDtrackerV1.h"
64 #include "TRD/AliTRDcalibDB.h"
66 #include "PWG1/TRD/macros/AliTRDperformanceTrain.h"
67 #include "PWG1/TRD/macros/AddTRDcheckESD.C"
68 #include "PWG1/TRD/macros/AddTRDinfoGen.C"
69 #include "PWG1/TRD/macros/AddTRDcheckDET.C"
70 #include "PWG1/TRD/macros/AddTRDefficiency.C"
71 #include "PWG1/TRD/macros/AddTRDresolution.C"
72 #include "PWG1/TRD/macros/AddTRDcheckPID.C"
73 #include "PWG1/TRD/macros/AddTRDv0Monitor.C"
76 #include "../TRD/macros/AliTRDperformanceTrain.h"
77 #include "../TRD/macros/helper.C"
80 const Char_t* Translate(Bool_t doCheckESD=kTRUE, Bool_t doCheckDET=kTRUE, Bool_t doEffic=kTRUE, Bool_t doResolution=kTRUE, Bool_t doCheckPID=kTRUE, Bool_t doV0Monitor=kTRUE);
81 Bool_t AddTrainPerformanceTRD(Char_t *trd="ALL", const Char_t *addMacroPath = "$ALICE_ROOT/PWG1/TRD/macros")
83 AliAnalysisManager *mgr = AliAnalysisManager::GetAnalysisManager();
85 Error("AddTrainPerformanceTRD", "AliAnalysisManager not set!");
89 // TRD data containers
90 AliAnalysisDataContainer *ci[kNOutSlots];
91 AliAnalysisDataContainer *ce[5];
93 Info("AddTrainPerformanceTRD", Form("Add Macros taken from %s", addMacroPath));
94 Info("AddTrainPerformanceTRD", Form("TRD wagons \"%s\"", trd));
95 Int_t bitmap = ParseOptions(trd);
96 for(Int_t it=0; it<NTRDQATASKS; it++){
97 if(!TSTBIT(bitmap, it)) continue;
98 if(gROOT->LoadMacro(Form("%s/Add%s.C+", addMacroPath, TString(fgkTRDtaskClassName[it])(3,20).Data()))) {
99 Error("AddTrainPerformanceTRD()", Form("Error loading %s task.", fgkTRDtaskClassName[it]));
105 AddTRDcheckESD(mgr); break;
107 AddTRDinfoGen(mgr, 0, NULL, ci); break;
110 ce[0]=ci[kEventInfo];
111 ce[1]=ci[kTracksBarrel];
113 ce[3]=ci[kTracksKink];
114 AddTRDcheckDET(mgr, bitmap, ce);
118 ce[0]=ci[kTracksBarrel];
120 ce[2]=ci[kTracksKink];
121 AddTRDefficiency(mgr, bitmap, ce);
125 ce[0]=ci[kTracksBarrel];
127 ce[2]=ci[kTracksKink];
128 AddTRDresolution(mgr, bitmap, ce);
132 ce[0]=ci[kTracksBarrel];
134 AddTRDcheckPID(mgr, bitmap, ce, &ce[2]);
137 // slots already mapped by checkPID
138 AddTRDv0Monitor(mgr, 0, ce);
141 Warning("AddTrainPerformanceTRD()", Form("No performance task registered at slot %d.", it));
147 const Char_t* Translate(Bool_t doCheckESD, Bool_t doCheckDET, Bool_t doEffic, Bool_t doResolution, Bool_t doCheckPID, Bool_t doCheckV0)
149 if( doCheckESD==kTRUE &&
152 doResolution==kTRUE &&
157 Bool_t kINDENT(kFALSE);
163 if(kINDENT) opt.Append(" ");
168 if(kINDENT) opt.Append(" ");
173 if(kINDENT) opt.Append(" ");
178 if(kINDENT) opt.Append(" ");
183 if(kINDENT) opt.Append(" ");
189 return (const Char_t*)opt.Data();