]> git.uio.no Git - u/mrichter/AliRoot.git/blame - PWG1/macros/AddTaskImpParRes.C
Added d0 vs phi and cluster type
[u/mrichter/AliRoot.git] / PWG1 / macros / AddTaskImpParRes.C
CommitLineData
a979993d 1AliAnalysisTaskSEImpParRes *AddTaskImpParRes(Bool_t readMC=kFALSE,Int_t selPdg=-1,Bool_t diamond=kTRUE)
77e570bf 2{
3 //
4 // Configuration for the study of the impact parameter resolution
5 //
6 // xianbao.yuan@pd.infn.it, andrea.dainese@pd.infn.it
7 //
8
9
10 // Get the pointer to the existing analysis manager via the static access method.
11 //==============================================================================
12 AliAnalysisManager *mgr = AliAnalysisManager::GetAnalysisManager();
13 if (!mgr) {
14 ::Error("AddTaskImpParRes", "No analysis manager to connect to.");
15 return NULL;
16 }
77e570bf 17
18 // Aanalysis task
19 AliAnalysisTaskSEImpParRes *d0ResTask = new AliAnalysisTaskSEImpParRes("d0ResAnalysis");
20 d0ResTask->SetDebugLevel(2);
21 d0ResTask->SetReadMC(readMC);
22 d0ResTask->SetSelectedPdg(selPdg);
a979993d 23 d0ResTask->SetUseDiamond(diamond);
77e570bf 24 mgr->AddTask(d0ResTask);
cbddc2a0 25
26 TString fname="ImpParRes.Performance";
27 if(selPdg>0) {fname+=selPdg;}
28 fname.Append(".root");
29
77e570bf 30
31 //
32 // Create containers for input/output
33 AliAnalysisDataContainer *cinputd0distr = mgr->CreateContainer("cinputd0distr",TChain::Class(),
34 AliAnalysisManager::kInputContainer);
35
36 AliAnalysisDataContainer *coutputd0ITSpureSARec = mgr->CreateContainer("coutputd0ITSpureSARec",TList::Class(),
37 AliAnalysisManager::kOutputContainer,
cbddc2a0 38 fname.Data());
77e570bf 39
40 AliAnalysisDataContainer *coutputd0ITSpureSASkip = mgr->CreateContainer("coutputd0ITSpureSASkip",TList::Class(),
41 AliAnalysisManager::kOutputContainer,
cbddc2a0 42 fname.Data());
77e570bf 43
44 AliAnalysisDataContainer *coutputd0allPointRec = mgr->CreateContainer("coutputd0allPointRec",TList::Class(),
45 AliAnalysisManager::kOutputContainer,
cbddc2a0 46 fname.Data());
77e570bf 47
48 AliAnalysisDataContainer *coutputd0allPointSkip = mgr->CreateContainer("coutputd0allPointSkip",TList::Class(),
49 AliAnalysisManager::kOutputContainer,
cbddc2a0 50 fname.Data());
77e570bf 51
52 AliAnalysisDataContainer *coutputd0partPointRec = mgr->CreateContainer("coutputd0partPointRec",TList::Class(),
53 AliAnalysisManager::kOutputContainer,
cbddc2a0 54 fname.Data());
77e570bf 55
56 AliAnalysisDataContainer *coutputd0partPointSkip = mgr->CreateContainer("coutputd0partPointSkip",TList::Class(),
57 AliAnalysisManager::kOutputContainer,
cbddc2a0 58 fname.Data());
77e570bf 59
60 AliAnalysisDataContainer *coutputd0onepointSPDRec = mgr->CreateContainer("coutputd0onepointSPDRec",TList::Class(),
61 AliAnalysisManager::kOutputContainer,
cbddc2a0 62 fname.Data());
77e570bf 63
64 AliAnalysisDataContainer *coutputd0onepointSPDSkip = mgr->CreateContainer("coutputd0onepointSPDSkip",TList::Class(),
65 AliAnalysisManager::kOutputContainer,
cbddc2a0 66 fname.Data());
77e570bf 67
68 AliAnalysisDataContainer *coutputd0postvTracRec = mgr->CreateContainer("coutputd0postvTracRec",TList::Class(),
69 AliAnalysisManager::kOutputContainer,
cbddc2a0 70 fname.Data());
77e570bf 71
72
73 AliAnalysisDataContainer *coutputd0postvTracSkip = mgr->CreateContainer("coutputd0postvTracSkip",TList::Class(),
74 AliAnalysisManager::kOutputContainer,
cbddc2a0 75 fname.Data());
77e570bf 76
77 AliAnalysisDataContainer *coutputd0negtvTracRec = mgr->CreateContainer("coutputd0negtvTracRec",TList::Class(),
78 AliAnalysisManager::kOutputContainer,
cbddc2a0 79 fname.Data());
77e570bf 80
81 AliAnalysisDataContainer *coutputd0negtvTracSkip = mgr->CreateContainer("coutputd0negtvTracSkip",TList::Class(),
82 AliAnalysisManager::kOutputContainer,
cbddc2a0 83 fname.Data());
77e570bf 84
85 AliAnalysisDataContainer *coutputd0pullAllpointRec = mgr->CreateContainer("coutputd0pullAllpointRec",TList::Class(),
86 AliAnalysisManager::kOutputContainer,
cbddc2a0 87 fname.Data());
77e570bf 88
89
90 AliAnalysisDataContainer *coutputd0pullAllpointSkip = mgr->CreateContainer("coutputd0pullAllpointSkip",TList::Class(),
91 AliAnalysisManager::kOutputContainer,
cbddc2a0 92 fname.Data());
77e570bf 93
94 AliAnalysisDataContainer *coutputd0onlyRefitRec = mgr->CreateContainer("coutputd0onlyRefitRec",TList::Class(),
95 AliAnalysisManager::kOutputContainer,
cbddc2a0 96 fname.Data());
77e570bf 97
98 AliAnalysisDataContainer *coutputd0onlyRefitSkip = mgr->CreateContainer("coutputd0onlyRefitSkip",TList::Class(),
99 AliAnalysisManager::kOutputContainer,
cbddc2a0 100 fname.Data());
101
77e570bf 102
0013ea25 103 AliAnalysisDataContainer *coutputd0sinThetaRec = mgr->CreateContainer("coutputd0sinThetaRec",TList::Class(),
cbddc2a0 104 AliAnalysisManager::kOutputContainer,
105 fname.Data());
106
107
0013ea25 108 AliAnalysisDataContainer *coutputd0sinThetaSkip = mgr->CreateContainer("coutputd0sinThetaSkip",TList::Class(),
77e570bf 109 AliAnalysisManager::kOutputContainer,
cbddc2a0 110 fname.Data());
77e570bf 111
cbddc2a0 112
113 AliAnalysisDataContainer *coutputd0allPointTrue = mgr->CreateContainer("coutputd0allPointTrue",TList::Class(),
114 AliAnalysisManager::kOutputContainer,
115 fname.Data());
116
117 AliAnalysisDataContainer *coutputd0postvTracTrue = mgr->CreateContainer("coutputd0postvTracTrue",TList::Class(),
118 AliAnalysisManager::kOutputContainer,
119 fname.Data());
120
121 AliAnalysisDataContainer *coutputd0negtvTracTrue = mgr->CreateContainer("coutputd0negtvTracTrue",TList::Class(),
122 AliAnalysisManager::kOutputContainer,
123 fname.Data());
124
125 AliAnalysisDataContainer *coutputd0pullAllpointTrue = mgr->CreateContainer("coutputd0pullAllpointTrue",TList::Class(),
126 AliAnalysisManager::kOutputContainer,
127 fname.Data());
128
129
130 AliAnalysisDataContainer *coutputd0phiAllpointSkip = mgr->CreateContainer("coutputd0phiAllpointSkip",TList::Class(),
131 AliAnalysisManager::kOutputContainer,
132 fname.Data());
133
134 AliAnalysisDataContainer *coutputd0phiPostvtracSkip = mgr->CreateContainer("coutputd0phiPostvtracSkip",TList::Class(),
135 AliAnalysisManager::kOutputContainer,
136 fname.Data());
137
138 AliAnalysisDataContainer *coutputd0phiNegtvtracSkip = mgr->CreateContainer("coutputd0phiNegtvtracSkip",TList::Class(),
139 AliAnalysisManager::kOutputContainer,
140 fname.Data());
141
142
143 AliAnalysisDataContainer *coutputd0clusterTypeSPD01Skip = mgr->CreateContainer("coutputd0clusterTypeSPD01Skip",TList::Class(),
144 AliAnalysisManager::kOutputContainer,
145 fname.Data());
146
147 AliAnalysisDataContainer *coutputd0clusterTypeSPD02Skip = mgr->CreateContainer("coutputd0clusterTypeSPD02Skip",TList::Class(),
148 AliAnalysisManager::kOutputContainer,
149 fname.Data());
150
151 AliAnalysisDataContainer *coutputd0clusterTypeSPD03Skip = mgr->CreateContainer("coutputd0clusterTypeSPD03Skip",TList::Class(),
152 AliAnalysisManager::kOutputContainer,
153 fname.Data());
154
155 AliAnalysisDataContainer *coutputd0clusterTypeSPD11Skip = mgr->CreateContainer("coutputd0clusterTypeSPD11Skip",TList::Class(),
156 AliAnalysisManager::kOutputContainer,
157 fname.Data());
158
159 AliAnalysisDataContainer *coutputd0clusterTypeSPD12Skip = mgr->CreateContainer("coutputd0clusterTypeSPD12Skip",TList::Class(),
160 AliAnalysisManager::kOutputContainer,
161 fname.Data());
162
163 AliAnalysisDataContainer *coutputd0clusterTypeSPD13Skip = mgr->CreateContainer("coutputd0clusterTypeSPD13Skip",TList::Class(),
164 AliAnalysisManager::kOutputContainer,
165 fname.Data());
166
167 AliAnalysisDataContainer *coutputd0Pt = mgr->CreateContainer("coutputd0Pt",TList::Class(),
168 AliAnalysisManager::kOutputContainer,
169 fname.Data());
170
77e570bf 171 AliAnalysisDataContainer *coutputNentries = mgr->CreateContainer("coutputNentries",TH1F::Class(),
172 AliAnalysisManager::kOutputContainer,
cbddc2a0 173 fname.Data());
77e570bf 174
175 AliAnalysisDataContainer *coutputEstimVtx = mgr->CreateContainer("coutputEstimVtx",TH1F::Class(),
176 AliAnalysisManager::kOutputContainer,
cbddc2a0 177 fname.Data());
178
77e570bf 179 mgr->ConnectInput(d0ResTask,0,mgr->GetCommonInputContainer());
180 //mgr->ConnectOutput(d0ResTask,0,mgr->GetCommonOutputContainer());
181 mgr->ConnectOutput(d0ResTask,1,coutputd0ITSpureSARec);
182 mgr->ConnectOutput(d0ResTask,2,coutputd0ITSpureSASkip);
183 mgr->ConnectOutput(d0ResTask,3,coutputd0allPointRec);
184 mgr->ConnectOutput(d0ResTask,4,coutputd0allPointSkip);
185 mgr->ConnectOutput(d0ResTask,5,coutputd0partPointRec);
186 mgr->ConnectOutput(d0ResTask,6,coutputd0partPointSkip);
187 mgr->ConnectOutput(d0ResTask,7,coutputd0onepointSPDRec);
188 mgr->ConnectOutput(d0ResTask,8,coutputd0onepointSPDSkip);
189 mgr->ConnectOutput(d0ResTask,9,coutputd0postvTracRec);
190 mgr->ConnectOutput(d0ResTask,10,coutputd0postvTracSkip);
191 mgr->ConnectOutput(d0ResTask,11,coutputd0negtvTracRec);
192 mgr->ConnectOutput(d0ResTask,12,coutputd0negtvTracSkip);
193 mgr->ConnectOutput(d0ResTask,13,coutputd0pullAllpointRec);
194 mgr->ConnectOutput(d0ResTask,14,coutputd0pullAllpointSkip);
195 mgr->ConnectOutput(d0ResTask,15,coutputd0onlyRefitRec);
196 mgr->ConnectOutput(d0ResTask,16,coutputd0onlyRefitSkip);
0013ea25 197 mgr->ConnectOutput(d0ResTask,17,coutputd0sinThetaRec);
198 mgr->ConnectOutput(d0ResTask,18,coutputd0sinThetaSkip);
cbddc2a0 199 mgr->ConnectOutput(d0ResTask,19,coutputd0allPointTrue);
200 mgr->ConnectOutput(d0ResTask,20,coutputd0postvTracTrue);
201 mgr->ConnectOutput(d0ResTask,21,coutputd0negtvTracTrue);
202 mgr->ConnectOutput(d0ResTask,22,coutputd0pullAllpointTrue);
203 mgr->ConnectOutput(d0ResTask,23,coutputd0phiAllpointSkip);
204 mgr->ConnectOutput(d0ResTask,24,coutputd0phiPostvtracSkip);
205 mgr->ConnectOutput(d0ResTask,25,coutputd0phiNegtvtracSkip);
206 mgr->ConnectOutput(d0ResTask,26,coutputd0clusterTypeSPD01Skip);
207 mgr->ConnectOutput(d0ResTask,27,coutputd0clusterTypeSPD02Skip);
208 mgr->ConnectOutput(d0ResTask,28,coutputd0clusterTypeSPD03Skip);
209 mgr->ConnectOutput(d0ResTask,29,coutputd0clusterTypeSPD11Skip);
210 mgr->ConnectOutput(d0ResTask,30,coutputd0clusterTypeSPD12Skip);
211 mgr->ConnectOutput(d0ResTask,31,coutputd0clusterTypeSPD13Skip);
212 mgr->ConnectOutput(d0ResTask,32,coutputd0Pt);
213 mgr->ConnectOutput(d0ResTask,33,coutputNentries);
214 mgr->ConnectOutput(d0ResTask,34,coutputEstimVtx);
77e570bf 215 return d0ResTask;
216}