]>
Commit | Line | Data |
---|---|---|
6d1c79ca | 1 | //------------------------------------------------------------------------------\r |
2 | // Implementation of AliComparisonEff class. It keeps information from \r | |
3 | // comparison of reconstructed and MC particle tracks. In addtion, \r | |
4 | // it keeps selection cuts used during comparison. The comparison \r | |
5 | // information is stored in the ROOT histograms. Analysis of these \r | |
6 | // histograms can be done by using Analyse() class function. The result of \r | |
3baa4bfd | 7 | // the analysis (histograms/graphs) are stored in the folder which is \r |
8 | // a data member of AliComparisonEff.\r | |
9 | // \r | |
6d1c79ca | 10 | // Author: J.Otwinowski 04/02/2008 \r |
11 | //------------------------------------------------------------------------------\r | |
12 | \r | |
13 | /*\r | |
3baa4bfd | 14 | \r |
15 | // after running comparison task, read the file, and get component\r | |
b4126c69 | 16 | gROOT->LoadMacro("$ALICE_ROOT/PWG1/Macros/LoadMyLibs.C");\r |
17 | LoadMyLibs();\r | |
6d1c79ca | 18 | TFile f("Output.root");\r |
35771050 | 19 | //AliComparisonEff * compObj = (AliComparisonEff*)f.Get("AliComparisonEff");\r |
20 | AliComparisonEff * compObj = (AliComparisonEff*)cOutput->FindObject("AliComparisonEff");\r | |
6d1c79ca | 21 | \r |
b4126c69 | 22 | // Analyse comparison data\r |
3baa4bfd | 23 | compObj->Analyse();\r |
24 | \r | |
25 | // the output histograms/graphs will be stored in the folder "folderEff" \r | |
26 | compObj->GetAnalysisFolder()->ls("*");\r | |
27 | \r | |
28 | // user can save whole comparison object (or only folder with anlysed histograms) \r | |
29 | // in the seperate output file (e.g.)\r | |
30 | TFile fout("Analysed_Eff.root","recreate");\r | |
31 | compObj->Write(); // compObj->GetAnalysisFolder()->Write();\r | |
32 | fout.Close();\r | |
6d1c79ca | 33 | \r |
6d1c79ca | 34 | */\r |
35 | \r | |
36 | \r | |
37 | #include <iostream>\r | |
38 | \r | |
39 | #include "TFile.h"\r | |
40 | #include "TCint.h"\r | |
41 | #include "TH3F.h"\r | |
42 | #include "TH2F.h"\r | |
43 | #include "TF1.h"\r | |
44 | #include "TProfile.h"\r | |
45 | #include "TProfile2D.h"\r | |
46 | #include "TGraph2D.h"\r | |
47 | #include "TCanvas.h"\r | |
48 | #include "TGraph.h"\r | |
49 | // \r | |
50 | #include "AliESDEvent.h" \r | |
51 | #include "AliESD.h"\r | |
52 | #include "AliESDfriend.h"\r | |
53 | #include "AliESDfriendTrack.h"\r | |
54 | #include "AliRecInfoCuts.h" \r | |
55 | #include "AliMCInfoCuts.h" \r | |
56 | #include "AliLog.h" \r | |
57 | //\r | |
58 | #include "AliMathBase.h"\r | |
59 | #include "AliTreeDraw.h" \r | |
f7a1cc68 | 60 | #include "AliMagF.h" \r |
6d1c79ca | 61 | #include "AliESDVertex.h" \r |
62 | #include "AliExternalTrackParam.h" \r | |
63 | #include "AliTracker.h" \r | |
64 | \r | |
65 | #include "AliMCInfo.h" \r | |
66 | #include "AliESDRecInfo.h" \r | |
67 | #include "AliComparisonEff.h" \r | |
68 | \r | |
69 | using namespace std;\r | |
70 | \r | |
71 | \r | |
72 | ClassImp(AliComparisonEff)\r | |
73 | \r | |
74 | //_____________________________________________________________________________\r | |
75 | AliComparisonEff::AliComparisonEff():\r | |
3baa4bfd | 76 | AliComparisonObject("AliComparisonEff"),\r |
6d1c79ca | 77 | \r |
78 | // histograms\r | |
79 | \r | |
80 | fMCPt(0),\r | |
81 | fMCRecPt(0),\r | |
82 | fMCRecPrimPt(0),\r | |
83 | fMCRecSecPt(0),\r | |
84 | \r | |
85 | fEffTPCPt(0),\r | |
86 | fEffTPCPtMC(0),\r | |
87 | fEffTPCPtF(0),\r | |
88 | //\r | |
89 | fEffTPCPt_P(0),\r | |
90 | fEffTPCPtMC_P(0),\r | |
91 | fEffTPCPtF_P(0),\r | |
92 | //\r | |
93 | fEffTPCPt_Pi(0),\r | |
94 | fEffTPCPtMC_Pi(0),\r | |
95 | fEffTPCPtF_Pi(0),\r | |
96 | //\r | |
97 | fEffTPCPt_K(0),\r | |
98 | fEffTPCPtMC_K(0),\r | |
99 | fEffTPCPtF_K(0),\r | |
100 | \r | |
101 | fEffTPCTan(0),\r | |
102 | fEffTPCTanMC(0),\r | |
103 | fEffTPCTanF(0),\r | |
104 | //\r | |
105 | fEffTPCPtTan(0),\r | |
106 | fEffTPCPtTanMC(0),\r | |
107 | fEffTPCPtTanF(0),\r | |
35771050 | 108 | // TPC+ITS\r |
109 | fEffTPCITSPt(0),\r | |
110 | fEffTPCITSTan(0),\r | |
6d1c79ca | 111 | \r |
112 | // Cuts \r | |
113 | fCutsRC(0), \r | |
114 | fCutsMC(0),\r | |
115 | \r | |
3baa4bfd | 116 | fVertex(0),\r |
117 | \r | |
118 | // histogram folder \r | |
119 | fAnalysisFolder(0)\r | |
6d1c79ca | 120 | {\r |
121 | // init vertex\r | |
122 | fVertex = new AliESDVertex();\r | |
123 | fVertex->SetXv(0.0); fVertex->SetYv(0.0); fVertex->SetZv(0.0); \r | |
124 | \r | |
125 | for(Int_t i=0; i<4; ++i)\r | |
126 | {\r | |
127 | fTPCPtDCASigmaIdeal[i]=0;\r | |
128 | fTPCPtDCASigmaFull[i]=0;\r | |
129 | fTPCPtDCASigmaDay0[i]=0;\r | |
130 | \r | |
131 | fTPCPtDCAXY[i]=0;\r | |
132 | fTPCPtDCAZ[i]=0;\r | |
133 | \r | |
134 | fTPCPtDCASigmaIdealPid[i]=0;\r | |
135 | fTPCPtDCASigmaFullPid[i]=0;\r | |
136 | fTPCPtDCASigmaDay0Pid[i]=0;\r | |
137 | \r | |
138 | fTPCPtDCAXYPid[i]=0; \r | |
139 | fTPCPtDCAZPid[i]=0; \r | |
140 | }\r | |
b4126c69 | 141 | \r |
3baa4bfd | 142 | Init();\r |
6d1c79ca | 143 | }\r |
144 | \r | |
145 | //_____________________________________________________________________________\r | |
146 | AliComparisonEff::~AliComparisonEff(){\r | |
147 | \r | |
148 | // \r | |
3baa4bfd | 149 | if(fMCPt) delete fMCPt; fMCPt=0;\r |
6d1c79ca | 150 | if(fMCRecPt) delete fMCRecPt; fMCRecPt=0;\r |
151 | if(fMCRecPrimPt) delete fMCRecPrimPt; fMCRecPrimPt=0;\r | |
152 | if(fMCRecSecPt) delete fMCRecSecPt; fMCRecSecPt=0;\r | |
153 | \r | |
154 | // \r | |
155 | if(fEffTPCPt) delete fEffTPCPt; fEffTPCPt=0;\r | |
156 | if(fEffTPCPtMC) delete fEffTPCPtMC; fEffTPCPtMC=0;\r | |
157 | if(fEffTPCPtF) delete fEffTPCPtF; fEffTPCPtF=0;\r | |
158 | \r | |
159 | // \r | |
160 | if(fEffTPCPt_P) delete fEffTPCPt_P; fEffTPCPt_P=0;\r | |
161 | if(fEffTPCPtMC_P) delete fEffTPCPtMC_P; fEffTPCPtMC_P=0;\r | |
162 | if(fEffTPCPtF_P) delete fEffTPCPtF_P; fEffTPCPtF_P=0;\r | |
163 | \r | |
164 | // \r | |
165 | if(fEffTPCPt_Pi) delete fEffTPCPt_Pi; fEffTPCPt_Pi=0;\r | |
166 | if(fEffTPCPtMC_Pi) delete fEffTPCPtMC_Pi; fEffTPCPtMC_Pi=0;\r | |
167 | if(fEffTPCPtF_Pi) delete fEffTPCPtF_Pi; fEffTPCPtF_Pi=0;\r | |
168 | \r | |
169 | // \r | |
170 | if(fEffTPCPt_K) delete fEffTPCPt_K; fEffTPCPt_K=0;\r | |
171 | if(fEffTPCPtMC_K) delete fEffTPCPtMC_K; fEffTPCPtMC_K=0;\r | |
172 | if(fEffTPCPtF_K) delete fEffTPCPtF_K; fEffTPCPtF_K=0;\r | |
173 | \r | |
174 | // \r | |
175 | if(fEffTPCTan) delete fEffTPCTan; fEffTPCTan=0;\r | |
176 | if(fEffTPCTanMC) delete fEffTPCTanMC; fEffTPCTanMC=0;\r | |
177 | if(fEffTPCTanF) delete fEffTPCTanF; fEffTPCTanF=0;\r | |
178 | \r | |
179 | //\r | |
180 | if(fEffTPCPtTan) delete fEffTPCPtTan; fEffTPCPtTan=0;\r | |
181 | if(fEffTPCPtTanMC) delete fEffTPCPtTanMC; fEffTPCPtTanMC=0;\r | |
182 | if(fEffTPCPtTanF) delete fEffTPCPtTanF; fEffTPCPtTanF=0;\r | |
183 | \r | |
35771050 | 184 | if(fEffTPCITSPt) delete fEffTPCITSPt; fEffTPCITSPt=0;\r |
185 | if(fEffTPCITSTan) delete fEffTPCITSTan; fEffTPCITSTan=0;\r | |
186 | \r | |
6d1c79ca | 187 | for(Int_t i=0; i<4; ++i)\r |
188 | {\r | |
189 | if(fTPCPtDCASigmaIdeal[i]) delete fTPCPtDCASigmaIdeal[i]; fTPCPtDCASigmaIdeal[i]=0;\r | |
190 | if(fTPCPtDCASigmaFull[i]) delete fTPCPtDCASigmaFull[i]; fTPCPtDCASigmaFull[i]=0;\r | |
191 | if(fTPCPtDCASigmaDay0[i]) delete fTPCPtDCASigmaDay0[i]; fTPCPtDCASigmaDay0[i]=0;\r | |
192 | \r | |
193 | if(fTPCPtDCAXY[i]) delete fTPCPtDCAXY[i]; fTPCPtDCAXY[i]=0;\r | |
194 | if(fTPCPtDCAZ[i]) delete fTPCPtDCAZ[i]; fTPCPtDCAZ[i]=0;\r | |
195 | \r | |
196 | if(fTPCPtDCASigmaIdealPid[i]) delete fTPCPtDCASigmaIdealPid[i]; fTPCPtDCASigmaIdealPid[i]=0;\r | |
197 | if(fTPCPtDCASigmaFullPid[i]) delete fTPCPtDCASigmaFullPid[i]; fTPCPtDCASigmaFullPid[i]=0;\r | |
198 | if(fTPCPtDCASigmaDay0Pid[i]) delete fTPCPtDCASigmaDay0Pid[i]; fTPCPtDCASigmaDay0Pid[i]=0;\r | |
199 | \r | |
200 | if(fTPCPtDCAXYPid[i]) delete fTPCPtDCAXYPid[i]; fTPCPtDCAXYPid[i]=0;\r | |
201 | if(fTPCPtDCAZPid[i]) delete fTPCPtDCAZPid[i]; fTPCPtDCAZPid[i]=0;\r | |
202 | }\r | |
3baa4bfd | 203 | \r |
204 | if(fAnalysisFolder) delete fAnalysisFolder; fAnalysisFolder=0;\r | |
6d1c79ca | 205 | }\r |
206 | \r | |
207 | //_____________________________________________________________________________\r | |
3baa4bfd | 208 | void AliComparisonEff::Init(){\r |
6d1c79ca | 209 | \r |
210 | // Init histograms\r | |
211 | //\r | |
212 | fMCPt = new TH1F("fMCPt","fMCPt",50,0.1,3); \r | |
213 | fMCPt->SetXTitle("p_{t}");\r | |
214 | fMCPt->SetYTitle("yield");\r | |
215 | \r | |
216 | fMCRecPt = new TH1F("fMCRecPt","fMCRecPt",50,0.1,3); \r | |
217 | fMCRecPt->SetXTitle("p_{t}");\r | |
218 | fMCRecPt->SetYTitle("yield");\r | |
219 | \r | |
220 | fMCRecPrimPt = new TH1F("fMCRecPrimPt","fMCRecPrimPt",50,0.1,3); \r | |
221 | fMCRecPrimPt->SetXTitle("p_{t}");\r | |
222 | fMCRecPrimPt->SetYTitle("yield");\r | |
223 | \r | |
224 | fMCRecSecPt = new TH1F("fMCRecSecPt","fMCRecSecPt",50,0.1,3); \r | |
225 | fMCRecSecPt->SetXTitle("p_{t}");\r | |
226 | fMCRecSecPt->SetYTitle("yield");\r | |
227 | \r | |
228 | // Efficiency as function of pt\r | |
229 | fEffTPCPt = new TProfile("Eff_pt","Eff_Pt",50,0.1,3); \r | |
230 | fEffTPCPt->SetXTitle("p_{t}");\r | |
231 | fEffTPCPt->SetYTitle("TPC Efficiency");\r | |
232 | \r | |
35771050 | 233 | // Efficiency as function of pt (TPC+ITS)\r |
234 | fEffTPCITSPt = new TProfile("Eff_pt_TPCITS","Eff_Pt_TPCITS",50,0.1,3); \r | |
235 | fEffTPCITSPt->SetXTitle("p_{t}");\r | |
236 | fEffTPCITSPt->SetYTitle("TPCITS Efficiency");\r | |
237 | \r | |
238 | \r | |
6d1c79ca | 239 | fEffTPCPtMC = new TProfile("MC_Eff_pt","MC_Eff_Pt",50,0.1,3); \r |
240 | fEffTPCPtMC->SetXTitle("p_{t}");\r | |
241 | fEffTPCPtMC->SetYTitle("MC TPC Efficiency");\r | |
242 | \r | |
243 | fEffTPCPtF = new TProfile("F_Eff_pt","F_Eff_Pt",50,0.1,3); \r | |
244 | fEffTPCPtF->SetXTitle("p_{t}");\r | |
245 | fEffTPCPtF->SetYTitle("TPC Findable Efficiency");\r | |
246 | \r | |
247 | // Efficiency as function of pt protons\r | |
248 | fEffTPCPt_P = new TProfile("Eff_pt_P","Eff_Pt_P",50,0.1,3); \r | |
249 | fEffTPCPt_P->SetXTitle("p_{t}");\r | |
250 | fEffTPCPt_P->SetYTitle("TPC Efficiency");\r | |
251 | \r | |
252 | fEffTPCPtMC_P = new TProfile("MC_Eff_pt_P","MC_Eff_Pt_P",50,0.1,3); \r | |
253 | fEffTPCPtMC_P->SetXTitle("p_{t}");\r | |
254 | fEffTPCPtMC_P->SetYTitle("MC TPC Efficiency");\r | |
255 | \r | |
256 | fEffTPCPtF_P = new TProfile("F_Eff_pt_P","F_Eff_Pt_P",50,0.1,3); \r | |
257 | fEffTPCPtF_P->SetXTitle("p_{t}");\r | |
258 | fEffTPCPtF_P->SetYTitle("TPC Findable Efficiency");\r | |
259 | \r | |
260 | // Efficiency as function of pt pions\r | |
261 | fEffTPCPt_Pi = new TProfile("Eff_pt_Pi","Eff_Pit_Pi",50,0.1,3); \r | |
262 | fEffTPCPt_Pi->SetXTitle("p_{t}");\r | |
263 | fEffTPCPt_Pi->SetYTitle("TPC Efficiency");\r | |
264 | \r | |
265 | fEffTPCPtMC_Pi = new TProfile("MC_Eff_pt_Pi","MC_Eff_Pit_Pi",50,0.1,3); \r | |
266 | fEffTPCPtMC_Pi->SetXTitle("p_{t}");\r | |
267 | fEffTPCPtMC_Pi->SetYTitle("MC TPC Efficiency");\r | |
268 | \r | |
269 | fEffTPCPtF_Pi = new TProfile("F_Eff_pt_Pi","F_Eff_Pit_Pi",50,0.1,3); \r | |
270 | fEffTPCPtF_Pi->SetXTitle("p_{t}");\r | |
271 | fEffTPCPtF_Pi->SetYTitle("TPC Findable Efficiency");\r | |
272 | \r | |
273 | // Efficiency as function of pt kaons\r | |
274 | fEffTPCPt_K = new TProfile("Eff_pt_K","Eff_Kt_K",50,0.1,3); \r | |
275 | fEffTPCPt_K->SetXTitle("p_{t}");\r | |
276 | fEffTPCPt_K->SetYTitle("TPC Efficiency");\r | |
277 | \r | |
278 | fEffTPCPtMC_K = new TProfile("MC_Eff_pt_K","MC_Eff_Kt_K",50,0.1,3); \r | |
279 | fEffTPCPtMC_K->SetXTitle("p_{t}");\r | |
280 | fEffTPCPtMC_K->SetYTitle("MC TPC Efficiency");\r | |
281 | \r | |
282 | fEffTPCPtF_K = new TProfile("F_Eff_pt_K","F_Eff_Kt_K",50,0.1,3); \r | |
283 | fEffTPCPtF_K->SetXTitle("p_{t}");\r | |
284 | fEffTPCPtF_K->SetYTitle("TPC Findable Efficiency");\r | |
285 | \r | |
286 | // Efficiency as function of tan(theta) \r | |
287 | fEffTPCTan = new TProfile("Eff_tan","Eff_tan",50,-2.5,2.5); \r | |
288 | fEffTPCTan->SetXTitle("tan(#theta)");\r | |
289 | fEffTPCTan->SetYTitle("TPC Efficiency");\r | |
290 | \r | |
35771050 | 291 | // Efficiency as function of tan(theta) (TPC+ITS)\r |
292 | fEffTPCITSTan = new TProfile("Eff_tan_TPCITS","Eff_tan_TPCITS",50,-2.5,2.5); \r | |
293 | fEffTPCITSTan->SetXTitle("tan(#theta)");\r | |
294 | fEffTPCITSTan->SetYTitle("TPCITS Efficiency");\r | |
295 | \r | |
296 | \r | |
6d1c79ca | 297 | fEffTPCTanMC = new TProfile("MC_Eff_tan","MC_Eff_tan",50,-2.5,2.5); \r |
298 | fEffTPCTanMC->SetXTitle("tan(#theta)");\r | |
299 | fEffTPCTanMC->SetYTitle("MC TPC Efficiency");\r | |
300 | \r | |
301 | fEffTPCTanF = new TProfile("F_Eff_tan","F_Eff_tan",50,-2.5,2.5); \r | |
302 | fEffTPCPtF->SetXTitle("tan(#theta)");\r | |
303 | fEffTPCPtF->SetYTitle("TPC Findable Efficiency");\r | |
304 | \r | |
305 | // Efficiency as function of pt and tan(theta) \r | |
306 | fEffTPCPtTan = new TProfile2D("Eff_pt_tan","Eff_Pt_tan",10,0.1,3,20,-2.,2.);\r | |
307 | fEffTPCPtTan->SetXTitle("tan(#theta)");\r | |
308 | fEffTPCPtTan->SetYTitle("p_{t}");\r | |
309 | \r | |
310 | fEffTPCPtTanMC = new TProfile2D("MC_Eff_pt_tan_MC","MC Eff Pt",10,0.1,3,20,-2.,2.);\r | |
311 | fEffTPCPtTanMC->SetXTitle("tan(#theta)");\r | |
312 | fEffTPCPtTanMC->SetYTitle("p_{t}");\r | |
313 | \r | |
314 | fEffTPCPtTanF = new TProfile2D("MC_Eff_pt_tan_F","MC Eff Pt",10,0.1,3,20,-2.,2.);\r | |
315 | fEffTPCPtTanF->SetXTitle("tan(#theta)");\r | |
316 | fEffTPCPtTanF->SetYTitle("p_{t}");\r | |
317 | \r | |
318 | char name[256];\r | |
319 | for(Int_t i=0; i<4; ++i)\r | |
320 | {\r | |
321 | sprintf(name, "fTPCPtDCASigmaIdeal_%d",i);\r | |
322 | fTPCPtDCASigmaIdeal[i] = new TH2F(name,name,50,0.1,3,100,0,100);\r | |
323 | \r | |
324 | sprintf(name, "fTPCPtDCASigmaFull_%d",i);\r | |
325 | fTPCPtDCASigmaFull[i] = new TH2F(name,name,50,0.1,3,100,0,100);\r | |
326 | \r | |
327 | sprintf(name, "fTPCPtDCASigmaDay0_%d",i);\r | |
328 | fTPCPtDCASigmaDay0[i] = new TH2F(name,name,50,0.1,3,100,0,100);\r | |
329 | \r | |
330 | sprintf(name, "fTPCPtDCAXY_%d",i);\r | |
331 | fTPCPtDCAXY[i]= new TH2F(name,name,50,0.1,3,100,0,100);\r | |
332 | sprintf(name, "fTPCPtDCAZ_%d",i);\r | |
333 | fTPCPtDCAZ[i]= new TH2F(name,name,50,0.1,3,100,0,100);\r | |
334 | \r | |
335 | sprintf(name, "fTPCPtDCASigmaIdealPid_%d",i);\r | |
336 | fTPCPtDCASigmaIdealPid[i] = new TH3F(name,name,50,0.1,3,100,0,100,5,0,5);\r | |
337 | \r | |
338 | sprintf(name, "fTPCPtDCASigmaFullPid_%d",i);\r | |
339 | fTPCPtDCASigmaFullPid[i]= new TH3F(name,name,50,0.1,3,100,0,100,5,0,5);\r | |
340 | \r | |
341 | sprintf(name, "fTPCPtDCASigmaDay0Pid_%d",i);\r | |
342 | fTPCPtDCASigmaDay0Pid[i]= new TH3F(name,name,50,0.1,3,100,0,100,5,0,5);\r | |
343 | \r | |
344 | sprintf(name, "fTPCPtDCAXYPid_%d",i);\r | |
345 | fTPCPtDCAXYPid[i]= new TH3F(name,name,50,0.1,3,100,0,100,5,0,5);\r | |
346 | \r | |
347 | sprintf(name, "fTPCPtDCAZPid_%d",i);\r | |
348 | fTPCPtDCAZPid[i]= new TH3F(name,name,50,0.1,3,100,0,100,5,0,5);\r | |
349 | }\r | |
6d1c79ca | 350 | \r |
3baa4bfd | 351 | // init cuts\r |
6d1c79ca | 352 | if(!fCutsMC) \r |
353 | AliDebug(AliLog::kError, "ERROR: Cannot find AliMCInfoCuts object");\r | |
354 | if(!fCutsRC) \r | |
355 | AliDebug(AliLog::kError, "ERROR: Cannot find AliRecInfoCuts object");\r | |
3baa4bfd | 356 | \r |
357 | // init folder\r | |
358 | fAnalysisFolder = CreateFolder("folderEff","Analysis Efficiency Folder");\r | |
6d1c79ca | 359 | }\r |
360 | \r | |
361 | //_____________________________________________________________________________\r | |
362 | void AliComparisonEff::Process(AliMCInfo* infoMC, AliESDRecInfo *infoRC)\r | |
363 | {\r | |
364 | // Fill efficiency comparison information\r | |
365 | \r | |
366 | AliExternalTrackParam *track = 0;\r | |
6d1c79ca | 367 | Double_t field = AliTracker::GetBz(); // nominal Bz field [kG]\r |
368 | Double_t kMaxD = 123456.0; // max distance\r | |
35771050 | 369 | Int_t clusterITS[200];\r |
6d1c79ca | 370 | \r |
371 | // systematics\r | |
372 | const Double_t kSigma2Full_xy = 0.25; // ExB full systematics [cm]\r | |
b4126c69 | 373 | const Double_t kSigma2Full_z = 5.0; // drift velocity (goofie) [cm] \r |
6d1c79ca | 374 | \r |
375 | const Double_t kSigma2Day0_xy = 0.02; // ExB [cm]\r | |
35771050 | 376 | const Double_t kSigma2Day0_z = 0.1; // drift velocity (goofie) [cm] \r |
6d1c79ca | 377 | \r |
378 | // \r | |
35771050 | 379 | Double_t DCASigmaIdeal=0;\r |
6d1c79ca | 380 | Double_t DCASigmaFull=0;\r |
381 | Double_t DCASigmaDay0=0;\r | |
382 | \r | |
383 | Double_t dca[2], cov[3]; // dca_xy, dca_z, sigma_xy, sigma_xy_z, sigma_z\r | |
384 | \r | |
385 | // distance to Prim. vertex \r | |
386 | const Double_t* dv = infoMC->GetVDist(); \r | |
387 | \r | |
388 | Float_t mcpt = infoMC->GetParticle().Pt();\r | |
389 | Float_t tantheta = TMath::Tan(infoMC->GetParticle().Theta()-TMath::Pi()*0.5);\r | |
390 | Bool_t isPrim = TMath::Sqrt(dv[0]*dv[0] + dv[1]*dv[1])<fCutsMC->GetMaxR() && TMath::Abs(dv[2])<fCutsMC->GetMaxVz();\r | |
391 | \r | |
392 | // calculate and set prim. vertex\r | |
393 | fVertex->SetXv( infoMC->GetParticle().Vx() - dv[0] );\r | |
394 | fVertex->SetYv( infoMC->GetParticle().Vy() - dv[1] );\r | |
395 | fVertex->SetZv( infoMC->GetParticle().Vz() - dv[2] );\r | |
396 | \r | |
397 | // Check selection cuts\r | |
398 | if (fCutsMC->IsPdgParticle(TMath::Abs(infoMC->GetParticle().GetPdgCode())) == kFALSE) return; \r | |
399 | \r | |
400 | // transform Pdg to Pid\r | |
401 | // Pdg convension is different for hadrons and leptons \r | |
402 | // (e.g. K+/K- = 321/-321; e+/e- = -11/11 ) \r | |
403 | Double_t pid = -1;\r | |
404 | if( TMath::Abs(infoMC->GetParticle().GetPdgCode())==fCutsMC->GetEM() ) pid = 0; \r | |
405 | if( TMath::Abs(infoMC->GetParticle().GetPdgCode())==fCutsMC->GetMuM() ) pid = 1; \r | |
406 | if( TMath::Abs(infoMC->GetParticle().GetPdgCode())==fCutsMC->GetPiP() ) pid = 2; \r | |
407 | if( TMath::Abs(infoMC->GetParticle().GetPdgCode())==fCutsMC->GetKP() ) pid = 3; \r | |
408 | if( TMath::Abs(infoMC->GetParticle().GetPdgCode())==fCutsMC->GetProt() ) pid = 4; \r | |
409 | \r | |
410 | //cout << "dv[0] " << dv[0] << " dv[1] " << dv[1] << " dv[2] " << dv[2] << endl; \r | |
411 | //cout << "v[0] " << fVertex->GetXv() << " v[1] " << fVertex->GetYv() << " v[2] " << fVertex->GetZv()<< endl; \r | |
412 | if (TMath::Abs(tantheta)<fCutsRC->GetMaxAbsTanTheta())\r | |
413 | {\r | |
18e1354a | 414 | if (infoRC->GetESDtrack() && infoRC->GetESDtrack()->GetTPCInnerParam()) \r |
6d1c79ca | 415 | {\r |
416 | if ((track = new AliExternalTrackParam(*infoRC->GetESDtrack()->GetTPCInnerParam())) != 0 )\r | |
417 | {\r | |
6d1c79ca | 418 | Bool_t bDCAStatus = track->PropagateToDCA(fVertex,field,kMaxD,dca,cov);\r |
3baa4bfd | 419 | if(bDCAStatus) {\r |
6d1c79ca | 420 | //\r |
421 | cov[2] = track->GetCovariance()[2];\r | |
422 | \r | |
423 | // Eff = infoRC->GetStatus(1)==3 && isPrim / isPrim;\r | |
424 | // Pt vs ( dca[0]^2/cov[0]^2 + dca[1]^2/cov[2]^2 ) \r | |
425 | // Pt vs ( dca[0]^2/(cov[0]^2 + kSigma2Full_xy) + dca[1]^2/( cov[2]^2 + kSigma2Full_z ) \r | |
426 | // Pt vs ( dca[0]^2/(cov[0]^2 + kSigma2_xy) + dca[1]^2/( cov[2]^2 + kSigma2_z ) \r | |
427 | \r | |
428 | if(cov[0]>0.0 && cov[2]>0.0)\r | |
429 | {\r | |
430 | DCASigmaIdeal = TMath::Power(dca[0],2)/cov[0] \r | |
431 | + TMath::Power(dca[1],2)/cov[2]; \r | |
432 | \r | |
433 | DCASigmaFull = TMath::Power(dca[0],2)/(cov[0]+kSigma2Full_xy) \r | |
434 | + TMath::Power(dca[1],2)/(cov[2]+kSigma2Full_z); \r | |
435 | \r | |
436 | DCASigmaDay0 = TMath::Power(dca[0],2)/(cov[0]+kSigma2Day0_xy) \r | |
437 | + TMath::Power(dca[1],2)/(cov[2]+kSigma2Day0_z); \r | |
438 | \r | |
439 | //cout << "dca[0] " << dca[0] << " dca[1] " << dca[1] << endl; \r | |
440 | //cout << "cov[0] " << cov[0] << " cov[2] " << cov[2] << endl; \r | |
441 | //cout << "DCASigmaIdeal " << DCASigmaIdeal << " DCASigmaFull " << DCASigmaFull << " DCASigmaDay0 " <<DCASigmaDay0 << endl; \r | |
442 | //cout << " -------------------------------------------------------- "<< endl; \r | |
443 | }\r | |
444 | \r | |
445 | // MC pt\r | |
446 | fMCPt->Fill(mcpt); \r | |
447 | if(infoRC->GetStatus(1)==3) fMCRecPt->Fill(mcpt); \r | |
448 | if(infoRC->GetStatus(1)==3 && isPrim) fMCRecPrimPt->Fill(mcpt); \r | |
449 | if(infoRC->GetStatus(1)==3 && !isPrim) fMCRecSecPt->Fill(mcpt); \r | |
450 | \r | |
451 | if(isPrim)\r | |
452 | {\r | |
453 | fTPCPtDCASigmaIdeal[0]->Fill(mcpt,DCASigmaIdeal);\r | |
454 | fTPCPtDCASigmaFull[0]->Fill(mcpt,DCASigmaFull);\r | |
455 | fTPCPtDCASigmaDay0[0]->Fill(mcpt,DCASigmaDay0);\r | |
456 | \r | |
457 | fTPCPtDCAXY[0]->Fill(mcpt,dca[0]);\r | |
458 | fTPCPtDCAZ[0]->Fill(mcpt,dca[1]);\r | |
459 | \r | |
460 | fTPCPtDCASigmaIdealPid[0]->Fill(mcpt,DCASigmaIdeal,pid);\r | |
461 | fTPCPtDCASigmaFullPid[0]->Fill(mcpt,DCASigmaFull,pid);\r | |
462 | fTPCPtDCASigmaDay0Pid[0]->Fill(mcpt,DCASigmaDay0,pid);\r | |
463 | \r | |
464 | fTPCPtDCAXYPid[0]->Fill(mcpt,dca[0],pid);\r | |
465 | fTPCPtDCAZPid[0]->Fill(mcpt,dca[1],pid);\r | |
466 | \r | |
467 | if(infoRC->GetStatus(1)==3)\r | |
468 | {\r | |
469 | fTPCPtDCASigmaIdeal[1]->Fill(mcpt,DCASigmaIdeal);\r | |
470 | fTPCPtDCASigmaFull[1]->Fill(mcpt,DCASigmaFull);\r | |
471 | fTPCPtDCASigmaDay0[1]->Fill(mcpt,DCASigmaDay0);\r | |
472 | \r | |
473 | fTPCPtDCAXY[1]->Fill(mcpt,dca[0]);\r | |
474 | fTPCPtDCAZ[1]->Fill(mcpt,dca[1]);\r | |
475 | \r | |
476 | fTPCPtDCASigmaIdealPid[1]->Fill(mcpt,DCASigmaIdeal,pid);\r | |
477 | fTPCPtDCASigmaFullPid[1]->Fill(mcpt,DCASigmaFull,pid);\r | |
478 | fTPCPtDCASigmaDay0Pid[1]->Fill(mcpt,DCASigmaDay0,pid);\r | |
479 | \r | |
480 | fTPCPtDCAXYPid[1]->Fill(mcpt,dca[0],pid);\r | |
481 | fTPCPtDCAZPid[1]->Fill(mcpt,dca[1],pid);\r | |
482 | }\r | |
483 | }\r | |
484 | \r | |
485 | // Cont = infoRC->GetStatus(1)==3 && !isPrim / infoRC->GetStatus(1)==3 \r | |
486 | // Pt vs ( dz[0]^2/cov[0]^2 + dz[1]^2/cov[2]^2 ) \r | |
487 | // Pt vs ( dz[0]^2/(cov[0]^2 + kSigma2Full_xy) + dz[1]^2/( cov[2]^2 + kSigma2Full_z ) \r | |
488 | // Pt vs ( dz[0]^2/(cov[0]^2 + kSigma2_xy) + dz[1]^2/( cov[2]^2 + kSigma2_z ) \r | |
489 | \r | |
490 | if(infoRC->GetStatus(1)==3) \r | |
491 | {\r | |
492 | fTPCPtDCASigmaIdeal[2]->Fill(mcpt,DCASigmaIdeal);\r | |
493 | fTPCPtDCASigmaFull[2]->Fill(mcpt,DCASigmaFull);\r | |
494 | fTPCPtDCASigmaDay0[2]->Fill(mcpt,DCASigmaDay0);\r | |
495 | \r | |
496 | fTPCPtDCAXY[2]->Fill(mcpt,dca[0]);\r | |
497 | fTPCPtDCAZ[2]->Fill(mcpt,dca[1]);\r | |
498 | \r | |
499 | fTPCPtDCASigmaIdealPid[2]->Fill(mcpt,DCASigmaIdeal,pid);\r | |
500 | fTPCPtDCASigmaFullPid[2]->Fill(mcpt,DCASigmaFull,pid);\r | |
501 | fTPCPtDCASigmaDay0Pid[2]->Fill(mcpt,DCASigmaDay0,pid);\r | |
502 | \r | |
503 | fTPCPtDCAXYPid[2]->Fill(mcpt,dca[0],pid);\r | |
504 | fTPCPtDCAZPid[2]->Fill(mcpt,dca[1],pid);\r | |
505 | \r | |
506 | if(isPrim==0)\r | |
507 | {\r | |
508 | fTPCPtDCASigmaIdeal[3]->Fill(mcpt,DCASigmaIdeal);\r | |
509 | fTPCPtDCASigmaFull[3]->Fill(mcpt,DCASigmaFull);\r | |
510 | fTPCPtDCASigmaDay0[3]->Fill(mcpt,DCASigmaDay0);\r | |
511 | \r | |
512 | fTPCPtDCAXY[3]->Fill(mcpt,dca[0]);\r | |
513 | fTPCPtDCAZ[3]->Fill(mcpt,dca[1]);\r | |
514 | \r | |
515 | fTPCPtDCASigmaIdealPid[3]->Fill(mcpt,DCASigmaIdeal,pid);\r | |
516 | fTPCPtDCASigmaFullPid[3]->Fill(mcpt,DCASigmaFull,pid);\r | |
517 | fTPCPtDCASigmaDay0Pid[3]->Fill(mcpt,DCASigmaDay0,pid);\r | |
518 | \r | |
519 | fTPCPtDCAXYPid[3]->Fill(mcpt,dca[0],pid);\r | |
520 | fTPCPtDCAZPid[3]->Fill(mcpt,dca[1],pid);\r | |
521 | }\r | |
522 | }\r | |
523 | delete track;\r | |
524 | }\r | |
525 | }\r | |
526 | } \r | |
527 | else \r | |
528 | {\r | |
529 | if(isPrim)\r | |
530 | {\r | |
531 | fTPCPtDCASigmaIdeal[0]->Fill(mcpt,0.0);\r | |
532 | fTPCPtDCASigmaFull[0]->Fill(mcpt,0.0);\r | |
533 | fTPCPtDCASigmaDay0[0]->Fill(mcpt,0.0);\r | |
534 | \r | |
535 | fTPCPtDCAXY[0]->Fill(mcpt,0.0);\r | |
536 | fTPCPtDCAZ[0]->Fill(mcpt,0.0);\r | |
537 | \r | |
538 | fTPCPtDCASigmaIdealPid[0]->Fill(mcpt,0.0,pid);\r | |
539 | fTPCPtDCASigmaFullPid[0]->Fill(mcpt,0.0,pid);\r | |
540 | fTPCPtDCASigmaDay0Pid[0]->Fill(mcpt,0.0,pid);\r | |
541 | \r | |
542 | fTPCPtDCAXYPid[0]->Fill(mcpt,0.0,pid);\r | |
543 | fTPCPtDCAZPid[0]->Fill(mcpt,0.0,pid);\r | |
544 | }\r | |
545 | }\r | |
546 | }\r | |
547 | \r | |
548 | // only primary particles\r | |
549 | if (!isPrim) return;\r | |
550 | \r | |
551 | // pt\r | |
552 | if (TMath::Abs(tantheta)<fCutsRC->GetMaxAbsTanTheta()){\r | |
553 | \r | |
554 | fEffTPCPt->Fill(mcpt, infoRC->GetStatus(1)==3);\r | |
555 | fEffTPCPtMC->Fill(mcpt, infoMC->GetRowsWithDigits()>fCutsMC->GetMinRowsWithDigits());\r | |
556 | if (infoMC->GetRowsWithDigits()>fCutsMC->GetMinRowsWithDigits()){\r | |
557 | fEffTPCPtF->Fill(mcpt, infoRC->GetStatus(1)==3);\r | |
558 | }\r | |
35771050 | 559 | \r |
560 | fEffTPCITSPt->Fill(mcpt, infoRC->GetStatus(1)==3 && (infoRC->GetESDtrack()->GetITSclusters(clusterITS)>fCutsRC->GetMinNClustersITS()));\r | |
6d1c79ca | 561 | \r |
562 | // protons\r | |
563 | if(TMath::Abs(infoMC->GetParticle().GetPdgCode())==fCutsMC->GetProt()) { \r | |
564 | fEffTPCPt_P->Fill(mcpt, infoRC->GetStatus(1)==3);\r | |
565 | fEffTPCPtMC_P->Fill(mcpt, infoMC->GetRowsWithDigits()>fCutsMC->GetMinRowsWithDigits());\r | |
566 | \r | |
567 | if(infoMC->GetRowsWithDigits()>fCutsMC->GetMinRowsWithDigits()) {\r | |
568 | fEffTPCPtF_P->Fill(mcpt, infoRC->GetStatus(1)==3);\r | |
569 | }\r | |
570 | }\r | |
571 | \r | |
572 | // pions\r | |
573 | if(TMath::Abs(infoMC->GetParticle().GetPdgCode())==fCutsMC->GetPiP()) {\r | |
574 | fEffTPCPt_Pi->Fill(mcpt, infoRC->GetStatus(1)==3);\r | |
575 | fEffTPCPtMC_Pi->Fill(mcpt, infoMC->GetRowsWithDigits()>fCutsMC->GetMinRowsWithDigits());\r | |
576 | \r | |
577 | if(infoMC->GetRowsWithDigits()>fCutsMC->GetMinRowsWithDigits()) {\r | |
578 | fEffTPCPtF_Pi->Fill(mcpt, infoRC->GetStatus(1)==3);\r | |
579 | }\r | |
580 | }\r | |
581 | \r | |
582 | // kaons\r | |
583 | if(TMath::Abs(infoMC->GetParticle().GetPdgCode())==fCutsMC->GetKP()) {\r | |
584 | fEffTPCPt_K->Fill(mcpt, infoRC->GetStatus(1)==3);\r | |
585 | fEffTPCPtMC_K->Fill(mcpt, infoMC->GetRowsWithDigits()>fCutsMC->GetMinRowsWithDigits());\r | |
586 | \r | |
587 | if(infoMC->GetRowsWithDigits()>fCutsMC->GetMinRowsWithDigits()) {\r | |
588 | fEffTPCPtF_K->Fill(mcpt, infoRC->GetStatus(1)==3);\r | |
589 | }\r | |
35771050 | 590 | }\r |
6d1c79ca | 591 | }\r |
592 | \r | |
593 | // theta\r | |
594 | if (TMath::Abs(mcpt)>fCutsRC->GetPtMin()){\r | |
595 | fEffTPCTan->Fill(tantheta, infoRC->GetStatus(1)==3);\r | |
596 | fEffTPCTanMC->Fill(tantheta, infoMC->GetRowsWithDigits()>fCutsMC->GetMinRowsWithDigits());\r | |
597 | if (infoMC->GetRowsWithDigits()>fCutsMC->GetMinRowsWithDigits()){\r | |
598 | fEffTPCTanF->Fill(tantheta, infoRC->GetStatus(1)==3);\r | |
599 | }\r | |
35771050 | 600 | \r |
601 | fEffTPCITSTan->Fill(tantheta, infoRC->GetStatus(1)==3 && (infoRC->GetESDtrack()->GetITSclusters(clusterITS)>fCutsRC->GetMinNClustersITS()));\r | |
6d1c79ca | 602 | }\r |
603 | \r | |
604 | // pt-theta\r | |
605 | fEffTPCPtTan->Fill(mcpt,tantheta,infoRC->GetStatus(1)==3);\r | |
606 | fEffTPCPtTanMC->Fill(mcpt,tantheta,infoMC->GetRowsWithDigits()>fCutsMC->GetMinRowsWithDigits()); \r | |
607 | if (infoMC->GetRowsWithDigits()>fCutsMC->GetMinRowsWithDigits()){\r | |
608 | fEffTPCPtTanF->Fill(mcpt,tantheta,infoRC->GetStatus(1)==3); \r | |
609 | }\r | |
610 | }\r | |
611 | \r | |
612 | //_____________________________________________________________________________\r | |
613 | void AliComparisonEff::Exec(AliMCInfo* infoMC, AliESDRecInfo *infoRC)\r | |
614 | {\r | |
615 | // Process comparison information\r | |
616 | Process(infoMC,infoRC);\r | |
617 | }\r | |
618 | \r | |
619 | //_____________________________________________________________________________\r | |
620 | Long64_t AliComparisonEff::Merge(TCollection* list) \r | |
621 | {\r | |
622 | // Merge list of objects (needed by PROOF)\r | |
623 | \r | |
624 | if (!list)\r | |
625 | return 0;\r | |
626 | \r | |
627 | if (list->IsEmpty())\r | |
628 | return 1;\r | |
629 | \r | |
630 | TIterator* iter = list->MakeIterator();\r | |
631 | TObject* obj = 0;\r | |
632 | \r | |
633 | // collection of generated histograms\r | |
634 | \r | |
635 | Int_t count=0;\r | |
636 | while((obj = iter->Next()) != 0) \r | |
637 | {\r | |
638 | AliComparisonEff* entry = dynamic_cast<AliComparisonEff*>(obj);\r | |
3baa4bfd | 639 | if (entry == 0) continue; \r |
6d1c79ca | 640 | \r |
641 | fMCPt->Add(entry->fMCPt);\r | |
642 | fMCRecPt->Add(entry->fMCRecPt);\r | |
643 | fMCRecPrimPt->Add(entry->fMCRecPrimPt);\r | |
644 | fMCRecSecPt->Add(entry->fMCRecSecPt);\r | |
645 | \r | |
646 | fEffTPCPt->Add(entry->fEffTPCPt);\r | |
647 | fEffTPCPtMC->Add(entry->fEffTPCPtMC);\r | |
648 | fEffTPCPtF->Add(entry->fEffTPCPtF);\r | |
649 | \r | |
650 | fEffTPCPt_P->Add(entry->fEffTPCPt_P);\r | |
651 | fEffTPCPtMC_P->Add(entry->fEffTPCPtMC_P);\r | |
652 | fEffTPCPtF_P->Add(entry->fEffTPCPtF_P);\r | |
653 | \r | |
654 | fEffTPCPt_Pi->Add(entry->fEffTPCPt_Pi);\r | |
655 | fEffTPCPtMC_Pi->Add(entry->fEffTPCPtMC_Pi);\r | |
656 | fEffTPCPtF_Pi->Add(entry->fEffTPCPtF_Pi);\r | |
657 | \r | |
658 | fEffTPCPt_K->Add(entry->fEffTPCPt_K);\r | |
659 | fEffTPCPtMC_K->Add(entry->fEffTPCPtMC_K);\r | |
660 | fEffTPCPtF_K->Add(entry->fEffTPCPtF_K);\r | |
661 | \r | |
662 | fEffTPCTan->Add(entry->fEffTPCTan);\r | |
663 | fEffTPCTanMC->Add(entry->fEffTPCTanMC);\r | |
664 | fEffTPCTanF->Add(entry->fEffTPCTanF);\r | |
665 | \r | |
666 | fEffTPCPtTan->Add(entry->fEffTPCPtTan);\r | |
667 | fEffTPCPtTanMC->Add(entry->fEffTPCPtTanMC);\r | |
668 | fEffTPCPtTanF->Add(entry->fEffTPCPtTanF);\r | |
669 | \r | |
35771050 | 670 | fEffTPCITSPt->Add(entry->fEffTPCITSPt);\r |
671 | fEffTPCITSTan->Add(entry->fEffTPCITSTan);\r | |
672 | \r | |
6d1c79ca | 673 | for(Int_t i=0; i<4; ++i)\r |
674 | {\r | |
675 | fTPCPtDCASigmaIdeal[i]->Add(entry->fTPCPtDCASigmaIdeal[i]);\r | |
676 | fTPCPtDCASigmaFull[i]->Add(entry->fTPCPtDCASigmaFull[i]);\r | |
677 | fTPCPtDCASigmaDay0[i]->Add(entry->fTPCPtDCASigmaDay0[i]);\r | |
678 | \r | |
679 | fTPCPtDCAXY[i]->Add(entry->fTPCPtDCAXY[i]);\r | |
680 | fTPCPtDCAZ[i]->Add(entry->fTPCPtDCAXY[i]);\r | |
681 | \r | |
682 | fTPCPtDCASigmaIdealPid[i]->Add(entry->fTPCPtDCASigmaIdealPid[i]);\r | |
683 | fTPCPtDCASigmaFullPid[i]->Add(entry->fTPCPtDCASigmaFullPid[i]);\r | |
684 | fTPCPtDCASigmaDay0Pid[i]->Add(entry->fTPCPtDCASigmaDay0Pid[i]);\r | |
685 | \r | |
686 | fTPCPtDCAXYPid[i]->Add(entry->fTPCPtDCAXYPid[i]);\r | |
687 | fTPCPtDCAZPid[i]->Add(entry->fTPCPtDCAXYPid[i]);\r | |
688 | }\r | |
689 | \r | |
690 | count++;\r | |
691 | }\r | |
692 | \r | |
693 | return count;\r | |
694 | }\r | |
695 | \r | |
696 | //_____________________________________________________________________________\r | |
697 | void AliComparisonEff::Analyse() \r | |
698 | {\r | |
3baa4bfd | 699 | // Analyse comparison information and store output histograms\r |
700 | // in the folder "folderEff" \r | |
701 | //\r | |
6d1c79ca | 702 | \r |
3baa4bfd | 703 | TH1::AddDirectory(kFALSE);\r |
704 | \r | |
6d1c79ca | 705 | AliComparisonEff * comp=this;\r |
b4126c69 | 706 | TObjArray *aFolderObj = new TObjArray;\r |
6d1c79ca | 707 | \r |
3baa4bfd | 708 | // calculate efficiency and contamination (4 sigma) \r |
6d1c79ca | 709 | TH1 *h_sigmaidealpid[20];\r |
710 | TH1 *h_sigmafullpid[20];\r | |
711 | TH1 *h_sigmaday0pid[20];\r | |
712 | \r | |
6d1c79ca | 713 | TH1 *h_sigmaidealpidtot[4];\r |
714 | TH1 *h_sigmafullpidtot[4];\r | |
715 | TH1 *h_sigmaday0pidtot[4];\r | |
716 | \r | |
6d1c79ca | 717 | char name[256];\r |
718 | char name1[256];\r | |
719 | Int_t idx;\r | |
720 | \r | |
721 | for(Int_t i=0; i<4; ++i)\r | |
722 | {\r | |
723 | //total\r | |
724 | comp->fTPCPtDCASigmaIdealPid[i]->GetYaxis()->SetRange(1,4);\r | |
725 | sprintf(name,"h_sigmaidealpidtot_%d",i);\r | |
726 | h_sigmaidealpidtot[i] = comp->fTPCPtDCASigmaIdealPid[i]->Project3D();\r | |
727 | h_sigmaidealpidtot[i]->SetName(name);\r | |
728 | \r | |
729 | comp->fTPCPtDCASigmaFullPid[i]->GetYaxis()->SetRange(1,4);\r | |
730 | sprintf(name,"h_sigmafullpidtot_%d",i);\r | |
731 | h_sigmafullpidtot[i] = comp->fTPCPtDCASigmaFullPid[i]->Project3D();\r | |
732 | h_sigmafullpidtot[i]->SetName(name);\r | |
733 | \r | |
734 | comp->fTPCPtDCASigmaDay0Pid[i]->GetYaxis()->SetRange(1,4);\r | |
735 | sprintf(name,"h_sigmaday0pidtot_%d",i);\r | |
736 | h_sigmaday0pidtot[i] = comp->fTPCPtDCASigmaDay0Pid[i]->Project3D();\r | |
737 | h_sigmaday0pidtot[i]->SetName(name);\r | |
738 | \r | |
739 | // pid wise\r | |
740 | for(Int_t j=0; j<5; ++j)\r | |
741 | {\r | |
742 | idx = i*5 + j;\r | |
743 | \r | |
744 | comp->fTPCPtDCASigmaIdealPid[i]->GetYaxis()->SetRange(1,4);\r | |
745 | comp->fTPCPtDCASigmaIdealPid[i]->GetZaxis()->SetRange(j+1,j+1);\r | |
746 | \r | |
747 | sprintf(name,"h_sigmaidealpid_%d",idx);\r | |
748 | h_sigmaidealpid[idx] = comp->fTPCPtDCASigmaIdealPid[i]->Project3D();\r | |
749 | h_sigmaidealpid[idx]->SetName(name);\r | |
750 | \r | |
751 | \r | |
752 | comp->fTPCPtDCASigmaFullPid[i]->GetYaxis()->SetRange(1,4);\r | |
753 | comp->fTPCPtDCASigmaFullPid[i]->GetZaxis()->SetRange(j+1,j+1);\r | |
754 | \r | |
755 | sprintf(name,"h_sigmafullpid_%d",idx);\r | |
756 | h_sigmafullpid[idx] = comp->fTPCPtDCASigmaFullPid[i]->Project3D();\r | |
757 | h_sigmafullpid[idx]->SetName(name);\r | |
758 | \r | |
759 | comp->fTPCPtDCASigmaDay0Pid[i]->GetYaxis()->SetRange(1,4);\r | |
760 | comp->fTPCPtDCASigmaDay0Pid[i]->GetZaxis()->SetRange(j+1,j+1);\r | |
761 | \r | |
762 | sprintf(name,"h_sigmaday0pid_%d",idx);\r | |
763 | h_sigmaday0pid[idx] = comp->fTPCPtDCASigmaDay0Pid[i]->Project3D();\r | |
764 | h_sigmaday0pid[idx]->SetName(name);\r | |
6d1c79ca | 765 | } \r |
766 | }\r | |
767 | \r | |
768 | // calculate efficiency and contamination (all pids)\r | |
769 | h_sigmaidealpidtot[0]->Sumw2();\r | |
770 | h_sigmaidealpidtot[1]->Divide(h_sigmaidealpidtot[0]);\r | |
771 | h_sigmaidealpidtot[2]->Sumw2();\r | |
772 | h_sigmaidealpidtot[3]->Divide(h_sigmaidealpidtot[2]);\r | |
773 | \r | |
774 | h_sigmafullpidtot[0]->Sumw2();\r | |
775 | h_sigmafullpidtot[1]->Divide(h_sigmafullpidtot[0]);\r | |
776 | h_sigmafullpidtot[2]->Sumw2();\r | |
777 | h_sigmafullpidtot[3]->Divide(h_sigmafullpidtot[2]);\r | |
778 | \r | |
779 | h_sigmaday0pidtot[0]->Sumw2();\r | |
780 | h_sigmaday0pidtot[1]->Divide(h_sigmaday0pidtot[0]);\r | |
781 | h_sigmaday0pidtot[2]->Sumw2();\r | |
782 | h_sigmaday0pidtot[3]->Divide(h_sigmaday0pidtot[2]);\r | |
783 | \r | |
784 | // calculate efficiency pid wise\r | |
07933e2a | 785 | for(Int_t jdx = 0; jdx<5; jdx++)\r |
6d1c79ca | 786 | {\r |
07933e2a | 787 | h_sigmaidealpid[jdx]->Sumw2();\r |
788 | h_sigmaidealpid[jdx+5]->Divide(h_sigmaidealpid[jdx]);\r | |
6d1c79ca | 789 | \r |
07933e2a | 790 | h_sigmafullpid[jdx]->Sumw2();\r |
791 | h_sigmafullpid[jdx+5]->Divide(h_sigmafullpid[jdx]);\r | |
6d1c79ca | 792 | \r |
07933e2a | 793 | h_sigmaday0pid[jdx]->Sumw2();\r |
794 | h_sigmaday0pid[jdx+5]->Divide(h_sigmaday0pid[jdx]);\r | |
6d1c79ca | 795 | }\r |
796 | \r | |
797 | // calculate cont. pid wise\r | |
07933e2a | 798 | for(Int_t jdx = 0; jdx<5; jdx++)\r |
6d1c79ca | 799 | {\r |
07933e2a | 800 | h_sigmaidealpid[jdx+15]->Divide(h_sigmaidealpidtot[2]);\r |
801 | h_sigmafullpid[jdx+15]->Divide(h_sigmafullpidtot[2]);\r | |
802 | h_sigmaday0pid[jdx+15]->Divide(h_sigmaday0pidtot[2]);\r | |
6d1c79ca | 803 | }\r |
804 | \r | |
6d1c79ca | 805 | TCanvas * c = new TCanvas("Efficiency","Track efficiency");\r |
806 | c->cd();\r | |
b4126c69 | 807 | c->Divide(1,2);\r |
6d1c79ca | 808 | \r |
6d1c79ca | 809 | //\r |
b4126c69 | 810 | c->cd(1);\r |
6d1c79ca | 811 | comp->fEffTPCTanF->SetXTitle("Tan(#theta)");\r |
812 | comp->fEffTPCTanF->SetYTitle("eff_{findable}");\r | |
3baa4bfd | 813 | comp->fEffTPCTanF->SetName("EffTanFindable");\r |
b4126c69 | 814 | comp->fEffTPCTanF->Draw();\r |
6d1c79ca | 815 | //\r |
b4126c69 | 816 | c->cd(2);\r |
6d1c79ca | 817 | comp->fEffTPCTan->SetXTitle("Tan(#theta)");\r |
818 | comp->fEffTPCTan->SetYTitle("eff_{all}");\r | |
3baa4bfd | 819 | comp->fEffTPCTan->SetName("EffTanAll");\r |
b4126c69 | 820 | comp->fEffTPCTan->Draw();\r |
6d1c79ca | 821 | \r |
b4126c69 | 822 | aFolderObj->Add(comp->fEffTPCTanF);\r |
823 | aFolderObj->Add(comp->fEffTPCTan);\r | |
6d1c79ca | 824 | \r |
b4126c69 | 825 | h_sigmaidealpidtot[1]->SetXTitle("p_{t}");\r |
826 | h_sigmaidealpidtot[1]->SetYTitle("efficiency");\r | |
827 | h_sigmaidealpidtot[1]->SetTitle("Eff_SigmaIdeal");\r | |
3baa4bfd | 828 | h_sigmaidealpidtot[1]->SetName("Eff_SigmaIdeal");\r |
b4126c69 | 829 | \r |
830 | h_sigmaidealpidtot[3]->SetXTitle("p_{t}");\r | |
831 | h_sigmaidealpidtot[3]->SetYTitle("contamination");\r | |
832 | h_sigmaidealpidtot[3]->SetTitle("Cont_SigmaIdeal");\r | |
3baa4bfd | 833 | h_sigmaidealpidtot[3]->SetName("Cont_SigmaIdeal");\r |
6d1c79ca | 834 | \r |
b4126c69 | 835 | aFolderObj->Add(h_sigmaidealpidtot[1]);\r |
836 | aFolderObj->Add(h_sigmaidealpidtot[3]);\r | |
3baa4bfd | 837 | \r |
b4126c69 | 838 | h_sigmafullpidtot[1]->SetXTitle("p_{t}");\r |
839 | h_sigmafullpidtot[1]->SetYTitle("efficiency");\r | |
840 | h_sigmafullpidtot[1]->SetTitle("Eff_SigmaFull");\r | |
3baa4bfd | 841 | h_sigmafullpidtot[1]->SetName("Eff_SigmaFull");\r |
b4126c69 | 842 | \r |
843 | h_sigmafullpidtot[3]->SetXTitle("p_{t}");\r | |
844 | h_sigmafullpidtot[3]->SetYTitle("contamination");\r | |
845 | h_sigmafullpidtot[3]->SetTitle("Cont_SigmaFull");\r | |
3baa4bfd | 846 | h_sigmafullpidtot[3]->SetName("Cont_SigmaFull");\r |
847 | \r | |
b4126c69 | 848 | aFolderObj->Add(h_sigmafullpidtot[1]);\r |
849 | aFolderObj->Add(h_sigmafullpidtot[3]);\r | |
3baa4bfd | 850 | \r |
b4126c69 | 851 | h_sigmaday0pidtot[1]->SetXTitle("p_{t}");\r |
852 | h_sigmaday0pidtot[1]->SetYTitle("efficiency");\r | |
853 | h_sigmaday0pidtot[1]->SetTitle("Eff_SigmaDay0");\r | |
3baa4bfd | 854 | h_sigmaday0pidtot[1]->SetName("Eff_SigmaDay0");\r |
b4126c69 | 855 | \r |
856 | h_sigmaday0pidtot[3]->SetXTitle("p_{t}");\r | |
857 | h_sigmaday0pidtot[3]->SetYTitle("contamination");\r | |
858 | h_sigmaday0pidtot[3]->SetTitle("Cont_SigmaDay0");\r | |
3baa4bfd | 859 | h_sigmaday0pidtot[3]->SetName("Cont_SigmaDay0");\r |
860 | \r | |
b4126c69 | 861 | aFolderObj->Add(h_sigmaday0pidtot[1]);\r |
862 | aFolderObj->Add(h_sigmaday0pidtot[3]);\r | |
6d1c79ca | 863 | \r |
07933e2a | 864 | for(Int_t jdx = 0; jdx<5; jdx++)\r |
6d1c79ca | 865 | {\r |
07933e2a | 866 | sprintf(name,"Eff_SigmaIdeal_%d",jdx);\r |
867 | sprintf(name1,"Cont_SigmaIdeal_%d",jdx);\r | |
6d1c79ca | 868 | \r |
b4126c69 | 869 | \r |
07933e2a | 870 | h_sigmaidealpid[jdx+5]->SetXTitle("p_{t}");\r |
871 | h_sigmaidealpid[jdx+5]->SetYTitle("efficiency");\r | |
872 | h_sigmaidealpid[jdx+5]->SetTitle(name);\r | |
873 | h_sigmaidealpid[jdx+5]->SetName(name);\r | |
b4126c69 | 874 | \r |
07933e2a | 875 | h_sigmaidealpid[jdx+15]->SetXTitle("p_{t}");\r |
876 | h_sigmaidealpid[jdx+15]->SetYTitle("contamination");\r | |
877 | h_sigmaidealpid[jdx+15]->SetTitle(name1);\r | |
878 | h_sigmaidealpid[jdx+15]->SetName(name1);\r | |
3baa4bfd | 879 | \r |
07933e2a | 880 | aFolderObj->Add(h_sigmaidealpid[jdx+5]);\r |
881 | aFolderObj->Add(h_sigmaidealpid[jdx+15]);\r | |
6d1c79ca | 882 | \r |
07933e2a | 883 | sprintf(name,"Eff_SigmaFull_%d",jdx);\r |
884 | sprintf(name1,"Cont_SigmaFull_%d",jdx);\r | |
6d1c79ca | 885 | \r |
07933e2a | 886 | h_sigmafullpid[jdx+5]->SetXTitle("p_{t}");\r |
887 | h_sigmafullpid[jdx+5]->SetYTitle("efficiency");\r | |
888 | h_sigmafullpid[jdx+5]->SetTitle(name);\r | |
889 | h_sigmafullpid[jdx+5]->SetName(name);\r | |
b4126c69 | 890 | \r |
07933e2a | 891 | h_sigmafullpid[jdx+15]->SetXTitle("p_{t}");\r |
892 | h_sigmafullpid[jdx+15]->SetYTitle("contamination");\r | |
893 | h_sigmafullpid[jdx+15]->SetTitle(name1);\r | |
894 | h_sigmafullpid[jdx+15]->SetName(name1);\r | |
3baa4bfd | 895 | \r |
07933e2a | 896 | aFolderObj->Add(h_sigmafullpid[jdx+5]);\r |
897 | aFolderObj->Add(h_sigmafullpid[jdx+15]);\r | |
6d1c79ca | 898 | \r |
07933e2a | 899 | sprintf(name,"Eff_SigmaDay0_%d",jdx);\r |
900 | sprintf(name1,"Cont_SigmaDay0_%d",jdx);\r | |
6d1c79ca | 901 | \r |
07933e2a | 902 | h_sigmaday0pid[jdx+5]->SetXTitle("p_{t}");\r |
903 | h_sigmaday0pid[jdx+5]->SetYTitle("efficiency");\r | |
904 | h_sigmaday0pid[jdx+5]->SetTitle(name);\r | |
905 | h_sigmaday0pid[jdx+5]->SetName(name);\r | |
b4126c69 | 906 | \r |
07933e2a | 907 | h_sigmaday0pid[jdx+15]->SetXTitle("p_{t}");\r |
908 | h_sigmaday0pid[jdx+15]->SetYTitle("contamination");\r | |
909 | h_sigmaday0pid[jdx+15]->SetTitle(name1);\r | |
910 | h_sigmaday0pid[jdx+15]->SetName(name1);\r | |
3baa4bfd | 911 | \r |
07933e2a | 912 | aFolderObj->Add(h_sigmaday0pid[jdx+5]);\r |
913 | aFolderObj->Add(h_sigmaday0pid[jdx+15]);\r | |
6d1c79ca | 914 | }\r |
915 | \r | |
b4126c69 | 916 | // export objects to analysis folder\r |
917 | fAnalysisFolder = ExportToFolder(aFolderObj);\r | |
3baa4bfd | 918 | \r |
b4126c69 | 919 | // delete only TObjArray\r |
920 | if(aFolderObj) delete aFolderObj;\r | |
3baa4bfd | 921 | }\r |
922 | \r | |
b4126c69 | 923 | //_____________________________________________________________________________\r |
924 | TFolder* AliComparisonEff::ExportToFolder(TObjArray * array) \r | |
925 | {\r | |
926 | // recreate folder avery time and export objects to new one\r | |
927 | //\r | |
928 | AliComparisonEff * comp=this;\r | |
929 | TFolder *folder = comp->GetAnalysisFolder();\r | |
930 | \r | |
931 | TString name, title;\r | |
932 | TFolder *newFolder = 0;\r | |
933 | Int_t i = 0;\r | |
934 | Int_t size = array->GetSize();\r | |
935 | \r | |
936 | if(folder) { \r | |
937 | // get name and title from old folder\r | |
938 | name = folder->GetName(); \r | |
939 | title = folder->GetTitle(); \r | |
940 | \r | |
941 | // delete old one\r | |
942 | delete folder;\r | |
943 | \r | |
944 | // create new one\r | |
945 | newFolder = CreateFolder(name.Data(),title.Data());\r | |
946 | newFolder->SetOwner();\r | |
947 | \r | |
948 | // add objects to folder\r | |
949 | while(i < size) {\r | |
950 | newFolder->Add(array->At(i));\r | |
951 | i++;\r | |
952 | }\r | |
953 | }\r | |
954 | \r | |
955 | return newFolder;\r | |
956 | }\r | |
957 | \r | |
958 | \r | |
3baa4bfd | 959 | //_____________________________________________________________________________\r |
960 | TFolder* AliComparisonEff::CreateFolder(TString name,TString title) { \r | |
961 | // create folder for analysed histograms\r | |
962 | //\r | |
963 | TFolder *folder = 0;\r | |
964 | folder = new TFolder(name.Data(),title.Data());\r | |
965 | \r | |
966 | return folder;\r | |
6d1c79ca | 967 | }\r |