]>
Commit | Line | Data |
---|---|---|
96b3109a | 1 | //------------------------------------------------------------------------------ |
2 | // Implementation of AliComparisonRes class. It keeps information from | |
3 | // comparison of reconstructed and MC particle tracks. In addtion, | |
4 | // it keeps selection cuts used during comparison. The comparison | |
5 | // information is stored in the ROOT histograms. Analysis of these | |
6 | // histograms can be done by using Analyse() class function. The result of | |
3baa4bfd | 7 | // the analysis (histograms/graphs) are stored in the folder which is |
8 | // a data member of AliComparisonRes. | |
9 | // | |
96b3109a | 10 | // Author: J.Otwinowski 04/02/2008 |
11 | //------------------------------------------------------------------------------ | |
12 | ||
13 | /* | |
3baa4bfd | 14 | |
15 | // after running comparison task, read the file, and get component | |
26e10df5 | 16 | gROOT->LoadMacro("$ALICE_ROOT/PWG1/Macros/LoadMyLibs.C"); |
17 | LoadMyLibs(); | |
18 | ||
96b3109a | 19 | TFile f("Output.root"); |
3baa4bfd | 20 | AliComparisonRes * compObj = (AliComparisonRes*)f.Get("AliComparisonRes"); |
21 | ||
22 | // analyse comparison data | |
23 | compObj->Analyse(); | |
24 | ||
25 | // the output histograms/graphs will be stored in the folder "folderRes" | |
26 | compObj->GetAnalysisFolder()->ls("*"); | |
27 | ||
28 | // user can save whole comparison object (or only folder with anlysed histograms) | |
29 | // in the seperate output file (e.g.) | |
30 | TFile fout("Analysed_Res.root","recreate"); | |
31 | compObj->Write(); // compObj->GetAnalysisFolder()->Write(); | |
32 | fout.Close(); | |
96b3109a | 33 | |
96b3109a | 34 | */ |
35 | ||
36 | #include <iostream> | |
37 | ||
38 | #include "TFile.h" | |
39 | #include "TCint.h" | |
40 | #include "TH3F.h" | |
41 | #include "TH2F.h" | |
42 | #include "TF1.h" | |
43 | #include "TProfile.h" | |
44 | #include "TProfile2D.h" | |
45 | #include "TGraph2D.h" | |
46 | #include "TCanvas.h" | |
47 | #include "TGraph.h" | |
48 | ||
49 | #include "AliESDEvent.h" | |
50 | #include "AliESD.h" | |
51 | #include "AliESDfriend.h" | |
52 | #include "AliESDfriendTrack.h" | |
53 | #include "AliESDVertex.h" | |
54 | #include "AliRecInfoCuts.h" | |
55 | #include "AliMCInfoCuts.h" | |
56 | #include "AliLog.h" | |
57 | #include "AliTracker.h" | |
58 | ||
59 | #include "AliMathBase.h" | |
60 | #include "AliTreeDraw.h" | |
61 | ||
62 | #include "AliMCInfo.h" | |
63 | #include "AliESDRecInfo.h" | |
64 | #include "AliComparisonRes.h" | |
65 | ||
66 | using namespace std; | |
67 | ||
68 | ClassImp(AliComparisonRes) | |
69 | ||
70 | //_____________________________________________________________________________ | |
71 | AliComparisonRes::AliComparisonRes(): | |
3baa4bfd | 72 | AliComparisonObject("AliComparisonRes"), |
96b3109a | 73 | |
74 | // Resolution | |
75 | fPtResolLPT(0), // pt resolution - low pt | |
76 | fPtResolHPT(0), // pt resolution - high pt | |
77 | fPtPullLPT(0), // pt resolution - low pt | |
78 | fPtPullHPT(0), // pt resolution - high pt | |
decf0997 | 79 | fPhiResolTan(0), //-> angular resolution |
80 | fTanResolTan(0), //-> angular resolution | |
81 | fPhiPullTan(0), //-> angular resolution | |
82 | fTanPullTan(0), //-> angular resolution | |
83 | ||
96b3109a | 84 | // |
85 | // Resolution constrained param | |
86 | // | |
87 | fCPhiResolTan(0), // angular resolution - constrained | |
88 | fCTanResolTan(0), // angular resolution - constrained | |
89 | fCPtResolTan(0), // pt resolution - constrained | |
90 | fCPhiPullTan(0), // angular resolution - constrained | |
91 | fCTanPullTan(0), // angular resolution - constrained | |
92 | fCPtPullTan(0), // pt resolution - constrained | |
93 | ||
94 | // | |
95 | // Parametrisation histograms | |
96 | // | |
97 | ||
3baa4bfd | 98 | f1Pt2ResolS1PtTPC(0), |
99 | f1Pt2ResolS1PtTPCITS(0), | |
100 | fYResolS1PtTPC(0), | |
101 | fYResolS1PtTPCITS(0), | |
102 | fZResolS1PtTPC(0), | |
103 | fZResolS1PtTPCITS(0), | |
104 | fPhiResolS1PtTPC(0), | |
105 | fPhiResolS1PtTPCITS(0), | |
106 | fThetaResolS1PtTPC(0), | |
107 | fThetaResolS1PtTPCITS(0), | |
96b3109a | 108 | |
109 | // constrained | |
3baa4bfd | 110 | fC1Pt2ResolS1PtTPC(0), |
111 | fC1Pt2ResolS1PtTPCITS(0), | |
112 | fCYResolS1PtTPC(0), | |
113 | fCYResolS1PtTPCITS(0), | |
114 | fCZResolS1PtTPC(0), | |
115 | fCZResolS1PtTPCITS(0), | |
116 | fCPhiResolS1PtTPC(0), | |
117 | fCPhiResolS1PtTPCITS(0), | |
118 | fCThetaResolS1PtTPC(0), | |
119 | fCThetaResolS1PtTPCITS(0), | |
96b3109a | 120 | |
121 | // vertex | |
122 | fVertex(0), | |
123 | ||
124 | // Cuts | |
125 | fCutsRC(0), | |
3baa4bfd | 126 | fCutsMC(0), |
127 | ||
128 | // histogram folder | |
129 | fAnalysisFolder(0) | |
96b3109a | 130 | { |
3baa4bfd | 131 | Init(); |
96b3109a | 132 | |
133 | // vertex (0,0,0) | |
134 | fVertex = new AliESDVertex(); | |
135 | fVertex->SetXv(0.0); | |
136 | fVertex->SetYv(0.0); | |
137 | fVertex->SetZv(0.0); | |
138 | } | |
139 | ||
140 | //_____________________________________________________________________________ | |
141 | AliComparisonRes::~AliComparisonRes(){ | |
142 | ||
143 | // Resolution histograms | |
144 | if(fPtResolLPT) delete fPtResolLPT; fPtResolLPT=0; | |
145 | if(fPtResolHPT) delete fPtResolHPT; fPtResolHPT=0; | |
146 | if(fPtPullLPT) delete fPtPullLPT; fPtPullLPT=0; | |
147 | if(fPtPullHPT) delete fPtPullHPT; fPtPullHPT=0; | |
148 | ||
decf0997 | 149 | if(fPhiResolTan) delete fPhiResolTan; fPhiResolTan=0; |
150 | if(fTanResolTan) delete fTanResolTan; fTanResolTan=0; | |
151 | if(fPhiPullTan) delete fPhiPullTan; fPhiPullTan=0; | |
152 | if(fTanPullTan) delete fTanPullTan; fTanPullTan=0; | |
153 | ||
96b3109a | 154 | // Resolution histograms (constrained param) |
155 | if(fCPhiResolTan) delete fCPhiResolTan; fCPhiResolTan=0; | |
156 | if(fCTanResolTan) delete fCTanResolTan; fCTanResolTan=0; | |
157 | if(fCPtResolTan) delete fCPtResolTan; fCPtResolTan=0; | |
158 | if(fCPhiPullTan) delete fCPhiPullTan; fCPhiPullTan=0; | |
159 | if(fCTanPullTan) delete fCTanPullTan; fCTanPullTan=0; | |
160 | if(fCPtPullTan) delete fCPtPullTan; fCPtPullTan=0; | |
161 | ||
162 | // Parametrisation histograms | |
163 | // | |
3baa4bfd | 164 | if(f1Pt2ResolS1PtTPC) delete f1Pt2ResolS1PtTPC; f1Pt2ResolS1PtTPC=0; |
165 | if(f1Pt2ResolS1PtTPCITS) delete f1Pt2ResolS1PtTPCITS; f1Pt2ResolS1PtTPCITS=0; | |
166 | if(fYResolS1PtTPC) delete fYResolS1PtTPC; fYResolS1PtTPC=0; | |
167 | if(fYResolS1PtTPCITS) delete fYResolS1PtTPCITS; fYResolS1PtTPCITS=0; | |
168 | if(fZResolS1PtTPC) delete fZResolS1PtTPC; fZResolS1PtTPC=0; | |
169 | if(fZResolS1PtTPCITS) delete fZResolS1PtTPCITS; fZResolS1PtTPCITS=0; | |
170 | if(fPhiResolS1PtTPC) delete fPhiResolS1PtTPC; fPhiResolS1PtTPC=0; | |
171 | if(fPhiResolS1PtTPCITS) delete fPhiResolS1PtTPCITS; fPhiResolS1PtTPCITS=0; | |
172 | if(fThetaResolS1PtTPC) delete fThetaResolS1PtTPC; fThetaResolS1PtTPC=0; | |
173 | if(fThetaResolS1PtTPCITS) delete fThetaResolS1PtTPCITS; fThetaResolS1PtTPCITS=0; | |
96b3109a | 174 | |
175 | // constrained | |
3baa4bfd | 176 | if(fC1Pt2ResolS1PtTPC) delete fC1Pt2ResolS1PtTPC; fC1Pt2ResolS1PtTPC=0; |
177 | if(fC1Pt2ResolS1PtTPCITS) delete fC1Pt2ResolS1PtTPCITS; fC1Pt2ResolS1PtTPCITS=0; | |
178 | if(fCYResolS1PtTPC) delete fCYResolS1PtTPC; fCYResolS1PtTPC=0; | |
179 | if(fCYResolS1PtTPCITS) delete fCYResolS1PtTPCITS; fCYResolS1PtTPCITS=0; | |
180 | if(fCZResolS1PtTPC) delete fCZResolS1PtTPC; fCZResolS1PtTPC=0; | |
181 | if(fCZResolS1PtTPCITS) delete fCZResolS1PtTPCITS; fCZResolS1PtTPCITS=0; | |
182 | if(fCPhiResolS1PtTPC) delete fCPhiResolS1PtTPC; fCPhiResolS1PtTPC=0; | |
183 | if(fCPhiResolS1PtTPCITS) delete fCPhiResolS1PtTPCITS; fCPhiResolS1PtTPCITS=0; | |
184 | if(fCThetaResolS1PtTPC) delete fCThetaResolS1PtTPC; fCThetaResolS1PtTPC=0; | |
185 | if(fCThetaResolS1PtTPCITS) delete fCThetaResolS1PtTPCITS; fCThetaResolS1PtTPCITS=0; | |
96b3109a | 186 | |
187 | if(fVertex) delete fVertex; fVertex=0; | |
188 | ||
3baa4bfd | 189 | if(fAnalysisFolder) delete fAnalysisFolder; fAnalysisFolder=0; |
96b3109a | 190 | } |
191 | ||
192 | //_____________________________________________________________________________ | |
3baa4bfd | 193 | void AliComparisonRes::Init(){ |
96b3109a | 194 | |
195 | // Init histograms | |
196 | fCPhiResolTan = new TH2F("CPhiResolTan","CPhiResolTan",50, -2,2,200,-0.025,0.025); | |
197 | fCPhiResolTan->SetXTitle("tan(#theta)"); | |
198 | fCPhiResolTan->SetYTitle("#Delta#phi"); | |
199 | ||
200 | fCTanResolTan = new TH2F("CTanResolTan","CTanResolTan",50, -2,2,200,-0.025,0.025); | |
201 | fCTanResolTan->SetXTitle("tan(#theta)"); | |
202 | fCTanResolTan->SetYTitle("#Delta#theta"); | |
203 | ||
204 | fCPtResolTan=new TH2F("CPtResol","CPtResol",50, -2,2,200,-0.2,0.2); | |
205 | fCPtResolTan->SetXTitle("Tan(#theta)"); | |
206 | fCPtResolTan->SetYTitle("#Deltap_{t}/p_{t}"); | |
207 | ||
208 | fCPhiPullTan = new TH2F("CPhiPullTan","CPhiPullTan",50, -2,2,200,-5,5); | |
209 | fCPhiPullTan->SetXTitle("Tan(#theta)"); | |
210 | fCPhiPullTan->SetYTitle("#Delta#phi/#Sigma"); | |
211 | ||
212 | fCTanPullTan = new TH2F("CTanPullTan","CTanPullTan",50, -2,2,200,-5,5); | |
213 | fCTanPullTan->SetXTitle("Tan(#theta)"); | |
214 | fCTanPullTan->SetYTitle("#Delta#theta/#Sigma"); | |
215 | ||
216 | fCPtPullTan=new TH2F("CPtPull","CPtPull",50, -2,2,200,-5,5); | |
217 | fCPtPullTan->SetXTitle("Tan(#theta)"); | |
218 | fCPtPullTan->SetYTitle("(1/mcp_{t}-1/p_{t})/#Sigma"); | |
219 | ||
220 | fPtResolLPT = new TH2F("Pt_resol_lpt","pt resol",10, 0.1,3,200,-0.2,0.2); | |
221 | fPtResolLPT->SetXTitle("p_{t}"); | |
222 | fPtResolLPT->SetYTitle("#Deltap_{t}/p_{t}"); | |
223 | ||
224 | fPtResolHPT = new TH2F("Pt_resol_hpt","pt resol",10, 2,100,200,-0.3,0.3); | |
225 | fPtResolHPT->SetXTitle("p_{t}"); | |
226 | fPtResolHPT->SetYTitle("#Deltap_{t}/p_{t}"); | |
227 | ||
228 | fPtPullLPT = new TH2F("Pt_pull_lpt","pt pull",10, 0.1,3,200,-6,6); | |
229 | fPtPullLPT->SetXTitle("p_{t}"); | |
230 | fPtPullLPT->SetYTitle("#Deltap_{t}/#Sigma"); | |
231 | ||
232 | fPtPullHPT = new TH2F("Pt_pull_hpt","pt pull",10,2,100,200,-6,6); | |
233 | fPtPullHPT->SetXTitle("p_{t}"); | |
234 | fPtPullHPT->SetYTitle("#Deltap_{t}/#Sigma"); | |
decf0997 | 235 | |
236 | fPhiResolTan = new TH2F("PhiResolTan","PhiResolTan",50, -2,2,200,-0.025,0.025); | |
237 | fPhiResolTan->SetXTitle("tan(#theta)"); | |
238 | fPhiResolTan->SetYTitle("#Delta#phi"); | |
239 | ||
240 | fTanResolTan = new TH2F("TanResolTan","TanResolTan",50, -2,2,200,-0.025,0.025); | |
241 | fTanResolTan->SetXTitle("tan(#theta)"); | |
242 | fTanResolTan->SetYTitle("#Delta#theta"); | |
243 | ||
244 | fPhiPullTan = new TH2F("PhiPullTan","PhiPullTan",50, -2,2,200,-5,5); | |
245 | fPhiPullTan->SetXTitle("Tan(#theta)"); | |
246 | fPhiPullTan->SetYTitle("#Delta#phi/#Sigma"); | |
247 | ||
248 | fTanPullTan = new TH2F("TanPullTan","TanPullTan",50, -2,2,200,-5,5); | |
249 | fTanPullTan->SetXTitle("Tan(#theta)"); | |
250 | fTanPullTan->SetYTitle("#Delta#theta/#Sigma"); | |
96b3109a | 251 | |
252 | // | |
253 | // Parametrisation histograms | |
254 | // | |
255 | ||
3baa4bfd | 256 | f1Pt2ResolS1PtTPC = new TH2F("f1Pt2ResolS1PtTPC","(1/mcpt-1/pt)/(1+1/mcpt)^2 vs sqrt(1/pt))",100,0,3,200,-0.010,0.010); |
257 | f1Pt2ResolS1PtTPC->SetXTitle("#sqrt{1/mcp_{t}}"); | |
258 | f1Pt2ResolS1PtTPC->SetYTitle("(1/mcp_{t}-1/p_{t})/(1+1/mcp_{t})^2)"); | |
96b3109a | 259 | |
3baa4bfd | 260 | f1Pt2ResolS1PtTPCITS = new TH2F("f1Pt2ResolS1PtTPCITS","(1/mcpt-1/pt)/(1+1/mcpt)^2 vs sqrt(1/pt))",100,0,3,200,-0.010,0.010); |
261 | f1Pt2ResolS1PtTPCITS->SetXTitle("#sqrt{1/mcp_{t}}"); | |
262 | f1Pt2ResolS1PtTPCITS->SetYTitle("(1/mcp_{t}-1/p_{t})/(1+1/mcp_{t})^2)"); | |
96b3109a | 263 | |
3baa4bfd | 264 | fYResolS1PtTPC = new TH2F("fYResolS1PtTPC","fYResolS1PtTPC",100, 0,3,200,-1.0,1.0); |
265 | fYResolS1PtTPC->SetXTitle("#sqrt{1/mcp_{t}}"); | |
266 | fYResolS1PtTPC->SetYTitle("#DeltaY"); | |
96b3109a | 267 | |
3baa4bfd | 268 | fYResolS1PtTPCITS = new TH2F("fYResolS1PtTPCITS","fYResolS1PtTPCITS",100, 0,3,200,-0.05,0.05); |
269 | fYResolS1PtTPCITS->SetXTitle("#sqrt{1/mcp_{t}}"); | |
270 | fYResolS1PtTPCITS->SetYTitle("#DeltaY"); | |
96b3109a | 271 | |
3baa4bfd | 272 | fZResolS1PtTPC = new TH2F("fZResolS1PtTPC","fZResolS1PtTPC",100, 0,3,200,-1.0,1.0); |
273 | fZResolS1PtTPC->SetXTitle("#sqrt{1/mcp_{t}}"); | |
274 | fZResolS1PtTPC->SetYTitle("#DeltaZ"); | |
96b3109a | 275 | |
3baa4bfd | 276 | fZResolS1PtTPCITS = new TH2F("fZResolS1PtTPCITS","fZResolS1PtTPCITS",100, 0,3,200,-0.05,0.05); |
277 | fZResolS1PtTPCITS->SetXTitle("#sqrt{1/mcp_{t}}"); | |
278 | fZResolS1PtTPCITS->SetYTitle("#DeltaZ"); | |
96b3109a | 279 | |
3baa4bfd | 280 | fPhiResolS1PtTPC = new TH2F("fPhiResolS1PtTPC","fPhiResolS1PtTPC",100, 0,3,200,-0.025,0.025); |
281 | fPhiResolS1PtTPC->SetXTitle("#sqrt{1/mcp_{t}}"); | |
282 | fPhiResolS1PtTPC->SetYTitle("#Delta#phi"); | |
96b3109a | 283 | |
3baa4bfd | 284 | fPhiResolS1PtTPCITS = new TH2F("fPhiResolS1PtTPCITS","fPhiResolS1PtTPCITS",100, 0,3,200,-0.01,0.01); |
285 | fPhiResolS1PtTPCITS->SetXTitle("#sqrt{1/mcp_{t}}"); | |
286 | fPhiResolS1PtTPCITS->SetYTitle("#Delta#phi"); | |
96b3109a | 287 | |
3baa4bfd | 288 | fThetaResolS1PtTPC = new TH2F("fThetaResolS1PtTPC","fThetaResolS1PtTPC",100, 0,3,200,-0.025,0.025); |
289 | fThetaResolS1PtTPC->SetXTitle("#sqrt{1/mcp_{t}}"); | |
290 | fThetaResolS1PtTPC->SetYTitle("#Delta#theta"); | |
96b3109a | 291 | |
3baa4bfd | 292 | fThetaResolS1PtTPCITS = new TH2F("fThetaResolS1PtTPCITS","fThetaResolS1PtTPCITS",100, 0,3,200,-0.01,0.01); |
293 | fThetaResolS1PtTPCITS->SetXTitle("#sqrt{1/mcp_{t}}"); | |
294 | fThetaResolS1PtTPCITS->SetYTitle("#Delta#theta"); | |
96b3109a | 295 | |
296 | // constrained | |
3baa4bfd | 297 | fC1Pt2ResolS1PtTPC = new TH2F("fC1Pt2ResolS1PtTPC","(1/mcpt-1/pt)/(1+1/mcpt)^2 vs 1/pt)",100,0,3,200,-0.010,0.010); |
298 | fC1Pt2ResolS1PtTPC->SetXTitle("#sqrt{1/mcp_{t}}"); | |
299 | fC1Pt2ResolS1PtTPC->SetYTitle("(1/mcp_{t}-1/p_{t})/(1+1/mcp_{t})^2)"); | |
96b3109a | 300 | |
3baa4bfd | 301 | fC1Pt2ResolS1PtTPCITS = new TH2F("fC1Pt2ResolS1PtTPCITS","(1/mcpt-1/pt)/(1+1/mcpt)^2 vs 1/pt)",100,0,3,200,-0.010,0.010); |
302 | fC1Pt2ResolS1PtTPCITS->SetXTitle("#sqrt{1/mcp_{t}}"); | |
303 | fC1Pt2ResolS1PtTPCITS->SetYTitle("(1/mcp_{t}-1/p_{t})/(1+1/mcp_{t})^2)"); | |
96b3109a | 304 | |
3baa4bfd | 305 | fCYResolS1PtTPC = new TH2F("fCYResolS1PtTPC","fCYResolS1PtTPC",100, 0,3,200,-1.0,1.0); |
306 | fCYResolS1PtTPC->SetXTitle("#sqrt{1/mcp_{t}}"); | |
307 | fCYResolS1PtTPC->SetYTitle("#DeltaY"); | |
96b3109a | 308 | |
3baa4bfd | 309 | fCYResolS1PtTPCITS = new TH2F("fCYResolS1PtTPCITS","fCYResolS1PtTPCITS",100, 0,3,200,-0.01,0.01); |
310 | fCYResolS1PtTPCITS->SetXTitle("#sqrt{1/mcp_{t}}"); | |
311 | fCYResolS1PtTPCITS->SetYTitle("#DeltaY"); | |
96b3109a | 312 | |
3baa4bfd | 313 | fCZResolS1PtTPC = new TH2F("fCZResolS1PtTPC","fCZResolS1PtTPC",100, 0,3,200,-1.0,1.0); |
314 | fCZResolS1PtTPC->SetXTitle("#sqrt{1/mcp_{t}}"); | |
315 | fCZResolS1PtTPC->SetYTitle("#DeltaZ"); | |
96b3109a | 316 | |
3baa4bfd | 317 | fCZResolS1PtTPCITS = new TH2F("fCZResolS1PtTPCITS","fCZResolS1PtTPCITS",100, 0,3,200,-0.025,0.025); |
318 | fCZResolS1PtTPCITS->SetXTitle("#sqrt{1/mcp_{t}}"); | |
319 | fCZResolS1PtTPCITS->SetYTitle("#DeltaZ"); | |
96b3109a | 320 | |
3baa4bfd | 321 | fCPhiResolS1PtTPC = new TH2F("fCPhiResolS1PtTPC","fCPhiResolS1PtTPC",100, 0,3,200,-0.025,0.025); |
322 | fCPhiResolS1PtTPC->SetXTitle("#sqrt{1/mcp_{t}}"); | |
323 | fCPhiResolS1PtTPC->SetYTitle("#Delta#phi"); | |
96b3109a | 324 | |
3baa4bfd | 325 | fCPhiResolS1PtTPCITS = new TH2F("fCPhiResolS1PtTPCITS","fCPhiResolS1PtTPCITS",100, 0,3,200,-0.003,0.003); |
326 | fCPhiResolS1PtTPCITS->SetXTitle("#sqrt{1/mcp_{t}}"); | |
327 | fCPhiResolS1PtTPCITS->SetYTitle("#Delta#phi"); | |
96b3109a | 328 | |
3baa4bfd | 329 | fCThetaResolS1PtTPC = new TH2F("fCThetaResolS1PtTPC","fCThetaResolS1PtTPC",100, 0,3,200,-0.025,0.025); |
330 | fCThetaResolS1PtTPC->SetXTitle("#sqrt{1/mcp_{t}}"); | |
331 | fCThetaResolS1PtTPC->SetYTitle("#Delta#theta"); | |
96b3109a | 332 | |
3baa4bfd | 333 | fCThetaResolS1PtTPCITS = new TH2F("fCThetaResolS1PtTPCITS","fCThetaResolS1PtTPCITS",100, 0,3,200,-0.005,0.005); |
334 | fCThetaResolS1PtTPCITS->SetXTitle("#sqrt{1/mcp_{t}}"); | |
335 | fCThetaResolS1PtTPCITS->SetYTitle("#Delta#theta"); | |
96b3109a | 336 | |
96b3109a | 337 | // Init cuts |
338 | if(!fCutsMC) | |
339 | AliDebug(AliLog::kError, "ERROR: Cannot find AliMCInfoCuts object"); | |
340 | if(!fCutsRC) | |
341 | AliDebug(AliLog::kError, "ERROR: Cannot find AliRecInfoCuts object"); | |
3baa4bfd | 342 | |
343 | // init folder | |
344 | fAnalysisFolder = CreateFolder("folderRes","Analysis Resolution Folder"); | |
96b3109a | 345 | } |
346 | ||
347 | //_____________________________________________________________________________ | |
348 | void AliComparisonRes::Process(AliMCInfo* infoMC, AliESDRecInfo *infoRC) | |
349 | { | |
350 | // Fill resolution comparison information | |
351 | AliExternalTrackParam *track = 0; | |
96b3109a | 352 | Double_t field = AliTracker::GetBz(); // nominal Bz field [kG] |
353 | Double_t kMaxD = 123456.0; // max distance | |
354 | ||
355 | Int_t clusterITS[200]; | |
356 | Double_t dca[2], cov[3]; // dca_xy, dca_z, sigma_xy, sigma_xy_z, sigma_z | |
357 | ||
decf0997 | 358 | Float_t deltaPt, pullPt, deltaPhi,pullPhi, deltaTan, pullTan, delta1Pt2, deltaY1Pt, deltaZ1Pt, deltaPhi1Pt, deltaTheta1Pt; |
96b3109a | 359 | Float_t deltaPtTPC, pullPtTPC, deltaPhiTPC, deltaTanTPC, delta1Pt2TPC, deltaY1PtTPC, deltaZ1PtTPC, deltaPhi1PtTPC, deltaTheta1PtTPC; |
360 | ||
361 | Float_t mcpt = infoMC->GetParticle().Pt(); | |
3baa4bfd | 362 | Float_t s1mcpt = TMath::Sqrt(1./infoMC->GetParticle().Pt()); |
decf0997 | 363 | Float_t tantheta = TMath::Tan(infoMC->GetParticle().Theta()-TMath::Pi()*0.5); |
96b3109a | 364 | |
365 | // distance to Prim. vertex | |
366 | const Double_t* dv = infoMC->GetVDist(); | |
367 | Bool_t isPrim = TMath::Sqrt(dv[0]*dv[0] + dv[1]*dv[1])<fCutsMC->GetMaxR() && TMath::Abs(dv[2])<fCutsMC->GetMaxVz(); | |
368 | ||
369 | // Check selection cuts | |
370 | if (fCutsMC->IsPdgParticle(TMath::Abs(infoMC->GetParticle().GetPdgCode())) == kFALSE) return; | |
371 | if (!isPrim) return; | |
96b3109a | 372 | if (infoRC->GetStatus(1)!=3) return; // TPC refit |
373 | if (!infoRC->GetESDtrack()) return; | |
374 | if (infoRC->GetESDtrack()->GetTPCNcls()<fCutsRC->GetMinNClustersTPC()) return; | |
375 | ||
376 | // calculate and set prim. vertex | |
377 | fVertex->SetXv( infoMC->GetParticle().Vx() - dv[0] ); | |
378 | fVertex->SetYv( infoMC->GetParticle().Vy() - dv[1] ); | |
379 | fVertex->SetZv( infoMC->GetParticle().Vz() - dv[2] ); | |
380 | ||
381 | deltaPt= (mcpt-infoRC->GetESDtrack()->Pt())/mcpt; | |
382 | pullPt= (1/mcpt-infoRC->GetESDtrack()->OneOverPt())/TMath::Sqrt(infoRC->GetESDtrack()->GetSigma1Pt2()); | |
383 | deltaPhi = TMath::ATan2(infoRC->GetESDtrack()->Py(),infoRC->GetESDtrack()->Px())- | |
384 | TMath::ATan2(infoMC->GetParticle().Py(),infoMC->GetParticle().Px()); | |
decf0997 | 385 | pullPhi = deltaPhi/TMath::Sqrt(infoRC->GetESDtrack()->GetSigmaSnp2()); |
96b3109a | 386 | |
387 | deltaTan = TMath::ATan2(infoRC->GetESDtrack()->Pz(),infoRC->GetESDtrack()->Pt())- | |
388 | TMath::ATan2(infoMC->GetParticle().Pz(),infoMC->GetParticle().Pt()); | |
decf0997 | 389 | pullTan = deltaPhi/TMath::Sqrt(infoRC->GetESDtrack()->GetSigmaSnp2()); |
96b3109a | 390 | |
391 | delta1Pt2 = (1/mcpt-infoRC->GetESDtrack()->OneOverPt())/TMath::Power(1+1/mcpt,2); | |
392 | deltaY1Pt = (infoMC->GetParticle().Vy()-infoRC->GetESDtrack()->GetY()) / (0.2+1/mcpt); | |
393 | deltaZ1Pt = (infoMC->GetParticle().Vz()-infoRC->GetESDtrack()->GetZ()) / (0.2+1/mcpt); | |
394 | deltaPhi1Pt = deltaPhi / (0.1+1/mcpt); | |
395 | deltaTheta1Pt = deltaTan / (0.1+1/mcpt); | |
396 | ||
96b3109a | 397 | // calculate track parameters at vertex |
398 | const AliExternalTrackParam *innerTPC = 0; | |
399 | if ((innerTPC = infoRC->GetESDtrack()->GetTPCInnerParam()) != 0) | |
400 | { | |
401 | if ((track = new AliExternalTrackParam(*infoRC->GetESDtrack()->GetTPCInnerParam())) != 0 ) | |
402 | { | |
96b3109a | 403 | Bool_t bDCAStatus = track->PropagateToDCA(fVertex,field,kMaxD,dca,cov); |
404 | ||
405 | // Fill parametrisation histograms (only TPC track) | |
3baa4bfd | 406 | if(bDCAStatus) |
96b3109a | 407 | { |
408 | deltaPtTPC= (mcpt-innerTPC->Pt())/mcpt; | |
409 | pullPtTPC= (1/mcpt-innerTPC->OneOverPt())/TMath::Sqrt(innerTPC->GetSigma1Pt2()); | |
410 | deltaPhiTPC = TMath::ATan2(innerTPC->Py(),innerTPC->Px())- | |
411 | TMath::ATan2(infoMC->GetParticle().Py(),infoMC->GetParticle().Px()); | |
412 | ||
413 | deltaTanTPC = TMath::ATan2(innerTPC->Pz(),innerTPC->Pt())- | |
414 | TMath::ATan2(infoMC->GetParticle().Pz(),infoMC->GetParticle().Pt()); | |
415 | ||
416 | delta1Pt2TPC = (1/mcpt-innerTPC->OneOverPt())/TMath::Power(1+1/mcpt,2); | |
417 | deltaY1PtTPC= (infoMC->GetParticle().Vy()-innerTPC->GetY()) / (0.2+1/mcpt); | |
418 | deltaZ1PtTPC = (infoMC->GetParticle().Vz()-innerTPC->GetZ()) / (0.2+1/mcpt); | |
419 | deltaPhi1PtTPC = deltaPhiTPC / (0.1+1/mcpt); | |
420 | deltaTheta1PtTPC = deltaTanTPC / (0.1+1/mcpt); | |
421 | ||
3baa4bfd | 422 | f1Pt2ResolS1PtTPC->Fill(s1mcpt,delta1Pt2TPC); |
423 | fYResolS1PtTPC->Fill(s1mcpt,deltaY1PtTPC); | |
424 | fZResolS1PtTPC->Fill(s1mcpt,deltaZ1PtTPC); | |
425 | fPhiResolS1PtTPC->Fill(s1mcpt,deltaPhi1PtTPC); | |
426 | fThetaResolS1PtTPC->Fill(s1mcpt,deltaTheta1PtTPC); | |
96b3109a | 427 | } |
428 | delete track; | |
429 | } | |
430 | } | |
431 | ||
432 | // TPC and ITS (nb. of clusters >2) in the system | |
433 | if(infoRC->GetESDtrack()->GetITSclusters(clusterITS)>2) | |
434 | { | |
3baa4bfd | 435 | f1Pt2ResolS1PtTPCITS->Fill(s1mcpt,delta1Pt2); |
436 | fYResolS1PtTPCITS->Fill(s1mcpt,deltaY1Pt); | |
437 | fZResolS1PtTPCITS->Fill(s1mcpt,deltaZ1Pt); | |
438 | fPhiResolS1PtTPCITS->Fill(s1mcpt,deltaPhi1Pt); | |
439 | fThetaResolS1PtTPCITS->Fill(s1mcpt,deltaTheta1Pt); | |
96b3109a | 440 | } |
441 | ||
442 | // Fill histograms | |
443 | fPtResolLPT->Fill(mcpt,deltaPt); | |
444 | fPtResolHPT->Fill(mcpt,deltaPt); | |
445 | fPtPullLPT->Fill(mcpt,pullPt); | |
446 | fPtPullHPT->Fill(mcpt,pullPt); | |
decf0997 | 447 | |
448 | fPhiResolTan->Fill(tantheta,deltaPhi); | |
449 | fPhiPullTan->Fill(tantheta,pullPhi); | |
450 | fTanResolTan->Fill(tantheta,deltaTan); | |
451 | fTanPullTan->Fill(tantheta,pullTan); | |
452 | ||
96b3109a | 453 | } |
454 | ||
455 | //_____________________________________________________________________________ | |
456 | void AliComparisonRes::ProcessConstrained(AliMCInfo* infoMC, AliESDRecInfo *infoRC) | |
457 | { | |
458 | // Fill resolution comparison information (constarained parameters) | |
459 | // | |
460 | AliExternalTrackParam *track = 0; | |
96b3109a | 461 | Double_t field = AliTracker::GetBz(); // nominal Bz field [kG] |
462 | Double_t kMaxD = 123456.0; // max distance | |
463 | ||
464 | Int_t clusterITS[200]; | |
465 | Double_t dca[2], cov[3]; // dca_xy, dca_z, sigma_xy, sigma_xy_z, sigma_z | |
466 | ||
467 | Float_t deltaPt, pullPt, deltaPhi, pullPhi, deltaTan, pullTan, delta1Pt2, deltaY1Pt, deltaZ1Pt, deltaPhi1Pt, deltaTheta1Pt; | |
468 | Float_t deltaPtTPC, pullPtTPC, deltaPhiTPC, deltaTanTPC, delta1Pt2TPC, deltaY1PtTPC, deltaZ1PtTPC, deltaPhi1PtTPC, deltaTheta1PtTPC; | |
469 | ||
470 | Float_t mcpt = infoMC->GetParticle().Pt(); | |
3baa4bfd | 471 | Float_t s1mcpt = TMath::Sqrt(1./infoMC->GetParticle().Pt()); |
96b3109a | 472 | Float_t tantheta = TMath::Tan(infoMC->GetParticle().Theta()-TMath::Pi()*0.5); |
473 | ||
474 | // distance to Prim. vertex | |
475 | const Double_t* dv = infoMC->GetVDist(); | |
476 | Bool_t isPrim = TMath::Sqrt(dv[0]*dv[0] + dv[1]*dv[1])<fCutsMC->GetMaxR() && TMath::Abs(dv[2])<fCutsMC->GetMaxVz(); | |
477 | ||
478 | // Check selection cuts | |
479 | if (fCutsMC->IsPdgParticle(TMath::Abs(infoMC->GetParticle().GetPdgCode())) == kFALSE) return; | |
480 | if (!isPrim) return; | |
481 | if (infoRC->GetStatus(1)!=3) return; | |
482 | if (!infoRC->GetESDtrack()) return; | |
483 | if (infoRC->GetESDtrack()->GetTPCNcls()<fCutsRC->GetMinNClustersTPC()) return; | |
484 | if (!infoRC->GetESDtrack()->GetConstrainedParam()) return; | |
485 | ||
486 | // calculate and set prim. vertex | |
487 | fVertex->SetXv( infoMC->GetParticle().Vx() - dv[0] ); | |
488 | fVertex->SetYv( infoMC->GetParticle().Vy() - dv[1] ); | |
489 | fVertex->SetZv( infoMC->GetParticle().Vz() - dv[2] ); | |
490 | ||
491 | // constrained parameters resolution | |
492 | const AliExternalTrackParam * cparam = infoRC->GetESDtrack()->GetConstrainedParam(); | |
493 | deltaPt= (mcpt-cparam->Pt())/mcpt; | |
494 | pullPt= (1/mcpt-cparam->OneOverPt())/TMath::Sqrt(cparam->GetSigma1Pt2()); | |
495 | deltaPhi = TMath::ATan2(cparam->Py(),cparam->Px())- | |
496 | TMath::ATan2(infoMC->GetParticle().Py(),infoMC->GetParticle().Px()); | |
497 | pullPhi = deltaPhi/TMath::Sqrt(cparam->GetSigmaSnp2()); | |
498 | deltaTan = TMath::ATan2(cparam->Pz(),cparam->Pt())-TMath::ATan2(infoMC->GetParticle().Pz(),infoMC->GetParticle().Pt()); | |
499 | pullTan = deltaPhi/TMath::Sqrt(cparam->GetSigmaSnp2()); | |
500 | ||
501 | ||
502 | delta1Pt2 = (1/mcpt-cparam->OneOverPt())/TMath::Power(1+1/mcpt,2); | |
503 | ||
504 | deltaY1Pt = (infoMC->GetParticle().Vy()-cparam->GetY()) / (0.2+1/mcpt); | |
505 | deltaZ1Pt = (infoMC->GetParticle().Vz()-cparam->GetZ()) / (0.2+1/mcpt); | |
506 | deltaPhi1Pt = deltaPhi / (0.1+1/mcpt); | |
507 | deltaTheta1Pt = deltaTan / (0.1+1/mcpt); | |
508 | ||
509 | // calculate track parameters at vertex | |
510 | const AliExternalTrackParam *innerTPC = 0; | |
511 | if ((innerTPC = infoRC->GetESDtrack()->GetTPCInnerParam()) != 0) | |
512 | { | |
513 | if ((track = new AliExternalTrackParam(*infoRC->GetESDtrack()->GetTPCInnerParam())) != 0 ) | |
514 | { | |
96b3109a | 515 | Bool_t bDCAStatus = track->PropagateToDCA(fVertex,field,kMaxD,dca,cov); |
516 | ||
517 | // Fill parametrisation histograms (only TPC track) | |
3baa4bfd | 518 | if(bDCAStatus) |
96b3109a | 519 | { |
520 | deltaPtTPC= (mcpt-innerTPC->Pt())/mcpt; | |
521 | pullPtTPC= (1/mcpt-innerTPC->OneOverPt())/TMath::Sqrt(innerTPC->GetSigma1Pt2()); | |
522 | deltaPhiTPC = TMath::ATan2(innerTPC->Py(),innerTPC->Px())- | |
523 | TMath::ATan2(infoMC->GetParticle().Py(),infoMC->GetParticle().Px()); | |
524 | ||
525 | deltaTanTPC = TMath::ATan2(innerTPC->Pz(),innerTPC->Pt())- | |
526 | TMath::ATan2(infoMC->GetParticle().Pz(),infoMC->GetParticle().Pt()); | |
527 | ||
528 | delta1Pt2TPC = (1/mcpt-innerTPC->OneOverPt())/TMath::Power(1+1/mcpt,2); | |
529 | deltaY1PtTPC= (infoMC->GetParticle().Vy()-innerTPC->GetY()) / (0.2+1/mcpt); | |
530 | deltaZ1PtTPC = (infoMC->GetParticle().Vz()-innerTPC->GetZ()) / (0.2+1/mcpt); | |
531 | deltaPhi1PtTPC = deltaPhiTPC / (0.1+1/mcpt); | |
532 | deltaTheta1PtTPC = deltaTanTPC / (0.1+1/mcpt); | |
533 | ||
3baa4bfd | 534 | fC1Pt2ResolS1PtTPC->Fill(s1mcpt,delta1Pt2TPC); |
535 | fCYResolS1PtTPC->Fill(s1mcpt,deltaY1PtTPC); | |
536 | fCZResolS1PtTPC->Fill(s1mcpt,deltaZ1PtTPC); | |
537 | fCPhiResolS1PtTPC->Fill(s1mcpt,deltaPhi1PtTPC); | |
538 | fCThetaResolS1PtTPC->Fill(s1mcpt,deltaTheta1PtTPC); | |
96b3109a | 539 | } |
540 | delete track; | |
541 | } | |
542 | } | |
543 | ||
544 | // TPC and ITS (nb. of clusters >2) in the system | |
545 | if(infoRC->GetESDtrack()->GetITSclusters(clusterITS)>2) | |
546 | { | |
3baa4bfd | 547 | fC1Pt2ResolS1PtTPCITS->Fill(s1mcpt,delta1Pt2); |
548 | fCYResolS1PtTPCITS->Fill(s1mcpt,deltaY1Pt); | |
549 | fCZResolS1PtTPCITS->Fill(s1mcpt,deltaZ1Pt); | |
550 | fCPhiResolS1PtTPCITS->Fill(s1mcpt,deltaPhi1Pt); | |
551 | fCThetaResolS1PtTPCITS->Fill(s1mcpt,deltaTheta1Pt); | |
96b3109a | 552 | } |
553 | ||
554 | // Fill histograms | |
555 | fCPtResolTan->Fill(tantheta,deltaPt); | |
556 | fCPtPullTan->Fill(tantheta,pullPt); | |
557 | fCPhiResolTan->Fill(tantheta,deltaPhi); | |
558 | fCPhiPullTan->Fill(tantheta,pullPhi); | |
559 | fCTanResolTan->Fill(tantheta,deltaTan); | |
560 | fCTanPullTan->Fill(tantheta,pullTan); | |
561 | } | |
562 | ||
563 | //_____________________________________________________________________________ | |
564 | void AliComparisonRes::Exec(AliMCInfo* infoMC, AliESDRecInfo *infoRC){ | |
565 | ||
566 | // Process comparison information | |
567 | Process(infoMC,infoRC); | |
568 | ProcessConstrained(infoMC,infoRC); | |
569 | } | |
570 | ||
571 | //_____________________________________________________________________________ | |
572 | TH1F* AliComparisonRes::MakeResol(TH2F * his, Int_t integ, Bool_t type){ | |
573 | // Create resolution histograms | |
574 | ||
575 | TH1F *hisr, *hism; | |
576 | if (!gPad) new TCanvas; | |
577 | hisr = AliTreeDraw::CreateResHistoI(his,&hism,integ); | |
578 | if (type) return hism; | |
579 | else | |
580 | return hisr; | |
581 | } | |
582 | ||
583 | //_____________________________________________________________________________ | |
584 | void AliComparisonRes::Analyse(){ | |
3baa4bfd | 585 | // Analyse comparison information and store output histograms |
586 | // in the folder "folderRes" | |
587 | // | |
588 | ||
589 | TH1::AddDirectory(kFALSE); | |
590 | ||
96b3109a | 591 | AliComparisonRes * comp=this; |
592 | TH1F *hiss=0; | |
b4126c69 | 593 | TObjArray *aFolderObj = new TObjArray; |
3baa4bfd | 594 | |
595 | // write results in the folder | |
96b3109a | 596 | |
597 | TCanvas * c = new TCanvas("Phi resol Tan","Phi resol Tan"); | |
598 | c->cd(); | |
599 | // | |
600 | hiss = comp->MakeResol(comp->fCPtResolTan,1,0); | |
601 | hiss->SetXTitle("Tan(#theta)"); | |
602 | hiss->SetYTitle("#sigmap_{t}/p_{t}"); | |
603 | hiss->Draw(); | |
3baa4bfd | 604 | hiss->SetName("CptResolTan"); |
605 | ||
b4126c69 | 606 | aFolderObj->Add(hiss); |
3baa4bfd | 607 | |
96b3109a | 608 | // |
decf0997 | 609 | hiss = comp->MakeResol(comp->fPtResolLPT,1,0); |
610 | hiss->SetXTitle("p_{t} (GeV)"); | |
611 | hiss->SetYTitle("#sigmap_{t}/p_{t}"); | |
612 | hiss->Draw(); | |
613 | hiss->SetName("PtResolPt"); | |
614 | ||
615 | aFolderObj->Add(hiss); | |
616 | ||
617 | // | |
618 | hiss = comp->MakeResol(comp->fPtResolLPT,1,1); | |
619 | hiss->SetXTitle("p_{t} (GeV)"); | |
620 | hiss->SetYTitle("mean #Deltap_{t}/p_{t}"); | |
621 | hiss->Draw(); | |
622 | hiss->SetName("PtMeanPt"); | |
623 | ||
624 | aFolderObj->Add(hiss); | |
625 | ||
626 | // | |
627 | hiss = comp->MakeResol(comp->fPhiResolTan,1,0); | |
96b3109a | 628 | hiss->SetXTitle("Tan(#theta)"); |
629 | hiss->SetYTitle("#sigma#phi (rad)"); | |
630 | hiss->Draw(); | |
3baa4bfd | 631 | hiss->SetName("PhiResolTan"); |
632 | ||
b4126c69 | 633 | aFolderObj->Add(hiss); |
96b3109a | 634 | // |
decf0997 | 635 | hiss = comp->MakeResol(comp->fTanResolTan,1,0); |
96b3109a | 636 | hiss->SetXTitle("Tan(#theta)"); |
637 | hiss->SetYTitle("#sigma#theta (rad)"); | |
638 | hiss->Draw(); | |
3baa4bfd | 639 | hiss->SetName("ThetaResolTan"); |
640 | ||
decf0997 | 641 | aFolderObj->Add(hiss); |
642 | ||
643 | // | |
644 | hiss = comp->MakeResol(comp->fPhiResolTan,1,1); | |
645 | hiss->SetXTitle("Tan(#theta)"); | |
646 | hiss->SetYTitle("mean #Delta#phi (rad)"); | |
647 | hiss->Draw(); | |
648 | hiss->SetName("PhiMeanTan"); | |
649 | ||
650 | aFolderObj->Add(hiss); | |
651 | // | |
652 | hiss = comp->MakeResol(comp->fTanResolTan,1,1); | |
653 | hiss->SetXTitle("Tan(#theta)"); | |
654 | hiss->SetYTitle("mean #Delta#theta (rad)"); | |
655 | hiss->Draw(); | |
656 | hiss->SetName("ThetaMeanTan"); | |
657 | ||
658 | aFolderObj->Add(hiss); | |
659 | ||
660 | // | |
661 | hiss = comp->MakeResol(comp->fCPhiResolTan,1,0); | |
662 | hiss->SetXTitle("Tan(#theta)"); | |
663 | hiss->SetYTitle("#sigma#phi (rad)"); | |
664 | hiss->Draw(); | |
665 | hiss->SetName("CPhiResolTan"); | |
666 | ||
667 | aFolderObj->Add(hiss); | |
668 | // | |
669 | hiss = comp->MakeResol(comp->fCTanResolTan,1,0); | |
670 | hiss->SetXTitle("Tan(#theta)"); | |
671 | hiss->SetYTitle("#sigma#theta (rad)"); | |
672 | hiss->Draw(); | |
673 | hiss->SetName("CThetaResolTan"); | |
674 | ||
b4126c69 | 675 | aFolderObj->Add(hiss); |
96b3109a | 676 | // |
677 | hiss = comp->MakeResol(comp->fCPtPullTan,1,0); | |
678 | hiss->SetXTitle("Tan(#theta)"); | |
679 | hiss->SetYTitle("1/mcp_{t}-1/p_{t}/#Sigma(1/p_{t})"); | |
680 | hiss->Draw(); | |
3baa4bfd | 681 | hiss->SetName("CptPullTan"); |
682 | ||
b4126c69 | 683 | aFolderObj->Add(hiss); |
96b3109a | 684 | // |
3baa4bfd | 685 | hiss = comp->MakeResol(comp->fC1Pt2ResolS1PtTPC,1,0); |
b4126c69 | 686 | hiss->SetXTitle("#sqrt(1/mcp_{t})"); |
96b3109a | 687 | hiss->SetYTitle("1/mcp_{t}-1/p_{t}/(1+1/p_{t})^2"); |
688 | hiss->Draw(); | |
3baa4bfd | 689 | hiss->SetName("C1Pt2ResolS1PtTPC"); |
690 | ||
b4126c69 | 691 | aFolderObj->Add(hiss); |
96b3109a | 692 | |
3baa4bfd | 693 | hiss = comp->MakeResol(comp->fC1Pt2ResolS1PtTPCITS,1,0); |
b4126c69 | 694 | hiss->SetXTitle("#sqrt(1/mcp_{t})"); |
96b3109a | 695 | hiss->SetYTitle("1/mcp_{t}-1/p_{t}/(1+1/p_{t})^2"); |
696 | hiss->Draw(); | |
3baa4bfd | 697 | hiss->SetName("C1Pt2ResolS1PtTPCITS"); |
698 | ||
b4126c69 | 699 | aFolderObj->Add(hiss); |
96b3109a | 700 | // |
3baa4bfd | 701 | hiss = comp->MakeResol(comp->fCYResolS1PtTPC,1,0); |
b4126c69 | 702 | hiss->SetXTitle("#sqrt(1/mcp_{t})"); |
96b3109a | 703 | hiss->SetYTitle("(mcy-y)/(0.2+1/mcp_{t})"); |
704 | hiss->Draw(); | |
3baa4bfd | 705 | hiss->SetName("CYResolS1PtTPC"); |
706 | ||
b4126c69 | 707 | aFolderObj->Add(hiss); |
96b3109a | 708 | |
3baa4bfd | 709 | hiss = comp->MakeResol(comp->fCYResolS1PtTPCITS,1,0); |
b4126c69 | 710 | hiss->SetXTitle("#sqrt(1/mcp_{t})"); |
96b3109a | 711 | hiss->SetYTitle("(mcy-y)/(0.2+1/mcp_{t})"); |
712 | hiss->Draw(); | |
3baa4bfd | 713 | hiss->SetName("CYResolS1PtTPCITS"); |
714 | ||
b4126c69 | 715 | aFolderObj->Add(hiss); |
96b3109a | 716 | // |
3baa4bfd | 717 | hiss = comp->MakeResol(comp->fCZResolS1PtTPC,1,0); |
b4126c69 | 718 | hiss->SetXTitle("#sqrt(1/mcp_{t})"); |
96b3109a | 719 | hiss->SetYTitle("(mcz-z)/(0.2+1/mcp_{t})"); |
720 | hiss->Draw(); | |
3baa4bfd | 721 | hiss->SetName("CZResolS1PtTPC"); |
722 | ||
b4126c69 | 723 | aFolderObj->Add(hiss); |
96b3109a | 724 | |
3baa4bfd | 725 | hiss = comp->MakeResol(comp->fCZResolS1PtTPCITS,1,0); |
b4126c69 | 726 | hiss->SetXTitle("#sqrt(1/mcp_{t})"); |
96b3109a | 727 | hiss->SetYTitle("(mcz-z)/(0.2+1/mcp_{t})"); |
728 | hiss->Draw(); | |
3baa4bfd | 729 | hiss->SetName("CZResolS1PtTPCITS"); |
730 | ||
b4126c69 | 731 | aFolderObj->Add(hiss); |
96b3109a | 732 | // |
3baa4bfd | 733 | hiss = comp->MakeResol(comp->fCPhiResolS1PtTPC,1,0); |
b4126c69 | 734 | hiss->SetXTitle("#sqrt(1/mcp_{t})"); |
96b3109a | 735 | hiss->SetYTitle("(mc#phi-#phi)/(0.1+1/mcp_{t})"); |
736 | hiss->Draw(); | |
3baa4bfd | 737 | hiss->SetName("CPhiResolS1PtTPC"); |
738 | ||
b4126c69 | 739 | aFolderObj->Add(hiss); |
96b3109a | 740 | |
3baa4bfd | 741 | hiss = comp->MakeResol(comp->fCPhiResolS1PtTPCITS,1,0); |
b4126c69 | 742 | hiss->SetXTitle("#sqrt(1/mcp_{t})"); |
96b3109a | 743 | hiss->SetYTitle("(mc#phi-#phi)/(0.1+1/mcp_{t})"); |
744 | hiss->Draw(); | |
3baa4bfd | 745 | hiss->SetName("CPhiResolS1PtTPCITS"); |
746 | ||
b4126c69 | 747 | aFolderObj->Add(hiss); |
96b3109a | 748 | // |
3baa4bfd | 749 | hiss = comp->MakeResol(comp->fCThetaResolS1PtTPC,1,0); |
b4126c69 | 750 | hiss->SetXTitle("#sqrt(1/mcp_{t})"); |
96b3109a | 751 | hiss->SetYTitle("(mc#theta-#theta)/(0.1+1/mcp_{t})"); |
752 | hiss->Draw(); | |
3baa4bfd | 753 | hiss->SetName("CThetaResolS1PtTPC"); |
754 | ||
b4126c69 | 755 | aFolderObj->Add(hiss); |
96b3109a | 756 | |
3baa4bfd | 757 | hiss = comp->MakeResol(comp->fCThetaResolS1PtTPCITS,1,0); |
b4126c69 | 758 | hiss->SetXTitle("#sqrt(1/mcp_{t})"); |
96b3109a | 759 | hiss->SetYTitle("(mc#theta-#theta)/(0.1+1/mcp_{t})"); |
760 | hiss->Draw(); | |
3baa4bfd | 761 | hiss->SetName("CThetaResolS1PtTPCITS"); |
762 | ||
b4126c69 | 763 | aFolderObj->Add(hiss); |
96b3109a | 764 | |
765 | // | |
3baa4bfd | 766 | hiss = comp->MakeResol(comp->f1Pt2ResolS1PtTPC,1,0); |
b4126c69 | 767 | hiss->SetXTitle("#sqrt(1/mcp_{t})"); |
96b3109a | 768 | hiss->SetYTitle("1/mcp_{t}-1/p_{t}/(1+1/p_{t})^2"); |
769 | hiss->Draw(); | |
3baa4bfd | 770 | hiss->SetName("OnePt2ResolS1PtTPC"); |
771 | ||
b4126c69 | 772 | aFolderObj->Add(hiss); |
96b3109a | 773 | |
3baa4bfd | 774 | hiss = comp->MakeResol(comp->f1Pt2ResolS1PtTPCITS,1,0); |
b4126c69 | 775 | hiss->SetXTitle("#sqrt(1/mcp_{t})"); |
96b3109a | 776 | hiss->SetYTitle("1/mcp_{t}-1/p_{t}/(1+1/p_{t})^2"); |
777 | hiss->Draw(); | |
3baa4bfd | 778 | hiss->SetName("OnePt2ResolS1PtTPCITS"); |
779 | ||
b4126c69 | 780 | aFolderObj->Add(hiss); |
96b3109a | 781 | // |
3baa4bfd | 782 | hiss = comp->MakeResol(comp->fYResolS1PtTPC,1,0); |
b4126c69 | 783 | hiss->SetXTitle("#sqrt(1/mcp_{t})"); |
96b3109a | 784 | hiss->SetYTitle("(mcy-y)/(0.2+1/mcp_{t})"); |
785 | hiss->Draw(); | |
3baa4bfd | 786 | hiss->SetName("YResolS1PtTPC"); |
787 | ||
b4126c69 | 788 | aFolderObj->Add(hiss); |
96b3109a | 789 | |
3baa4bfd | 790 | hiss = comp->MakeResol(comp->fYResolS1PtTPCITS,1,0); |
b4126c69 | 791 | hiss->SetXTitle("#sqrt(1/mcp_{t})"); |
96b3109a | 792 | hiss->SetYTitle("(mcy-y)/(0.2+1/mcp_{t})"); |
793 | hiss->Draw(); | |
3baa4bfd | 794 | hiss->SetName("YResolS1PtTPCITS"); |
795 | ||
b4126c69 | 796 | aFolderObj->Add(hiss); |
96b3109a | 797 | // |
3baa4bfd | 798 | hiss = comp->MakeResol(comp->fZResolS1PtTPC,1,0); |
b4126c69 | 799 | hiss->SetXTitle("#sqrt(1/mcp_{t})"); |
96b3109a | 800 | hiss->SetYTitle("(mcz-z)/(0.2+1/mcp_{t})"); |
801 | hiss->Draw(); | |
3baa4bfd | 802 | hiss->SetName("ZResolS1PtTPC"); |
803 | ||
b4126c69 | 804 | aFolderObj->Add(hiss); |
96b3109a | 805 | |
3baa4bfd | 806 | hiss = comp->MakeResol(comp->fZResolS1PtTPCITS,1,0); |
b4126c69 | 807 | hiss->SetXTitle("#sqrt(1/mcp_{t})"); |
96b3109a | 808 | hiss->SetYTitle("(mcz-z)/(0.2+1/mcp_{t})"); |
809 | hiss->Draw(); | |
3baa4bfd | 810 | hiss->SetName("ZResolS1PtTPCITS"); |
811 | ||
b4126c69 | 812 | aFolderObj->Add(hiss); |
96b3109a | 813 | // |
3baa4bfd | 814 | hiss = comp->MakeResol(comp->fPhiResolS1PtTPC,1,0); |
b4126c69 | 815 | hiss->SetXTitle("#sqrt(1/mcp_{t})"); |
96b3109a | 816 | hiss->SetYTitle("(mc#phi-#phi)/(0.1+1/mcp_{t})"); |
817 | hiss->Draw(); | |
3baa4bfd | 818 | hiss->SetName("PhiResolS1PtTPC"); |
819 | ||
b4126c69 | 820 | aFolderObj->Add(hiss); |
96b3109a | 821 | |
3baa4bfd | 822 | hiss = comp->MakeResol(comp->fPhiResolS1PtTPCITS,1,0); |
b4126c69 | 823 | hiss->SetXTitle("#sqrt(1/mcp_{t})"); |
96b3109a | 824 | hiss->SetYTitle("(mc#phi-#phi)/(0.1+1/mcp_{t})"); |
825 | hiss->Draw(); | |
3baa4bfd | 826 | hiss->SetName("PhiResolS1PtTPCITS"); |
827 | ||
b4126c69 | 828 | aFolderObj->Add(hiss); |
96b3109a | 829 | // |
3baa4bfd | 830 | hiss = comp->MakeResol(comp->fThetaResolS1PtTPC,1,0); |
b4126c69 | 831 | hiss->SetXTitle("#sqrt(1/mcp_{t})"); |
96b3109a | 832 | hiss->SetYTitle("(mc#theta-#theta)/(0.1+1/mcp_{t})"); |
833 | hiss->Draw(); | |
3baa4bfd | 834 | hiss->SetName("ThetaResolS1PtTPC"); |
835 | ||
b4126c69 | 836 | aFolderObj->Add(hiss); |
96b3109a | 837 | |
3baa4bfd | 838 | hiss = comp->MakeResol(comp->fThetaResolS1PtTPCITS,1,0); |
b4126c69 | 839 | hiss->SetXTitle("#sqrt(1/mcp_{t})"); |
96b3109a | 840 | hiss->SetYTitle("(mc#theta-#theta)/(0.1+1/mcp_{t})"); |
841 | hiss->Draw(); | |
3baa4bfd | 842 | hiss->SetName("ThetaResolS1PtTPCITS"); |
843 | ||
b4126c69 | 844 | aFolderObj->Add(hiss); |
845 | ||
846 | // export objects to analysis folder | |
847 | fAnalysisFolder = ExportToFolder(aFolderObj); | |
848 | ||
849 | // delete only TObjArray | |
850 | if(aFolderObj) delete aFolderObj; | |
851 | } | |
852 | ||
853 | //_____________________________________________________________________________ | |
854 | TFolder* AliComparisonRes::ExportToFolder(TObjArray * array) | |
855 | { | |
856 | // recreate folder avery time and export objects to new one | |
857 | // | |
858 | AliComparisonRes * comp=this; | |
859 | TFolder *folder = comp->GetAnalysisFolder(); | |
860 | ||
861 | TString name, title; | |
862 | TFolder *newFolder = 0; | |
863 | Int_t i = 0; | |
864 | Int_t size = array->GetSize(); | |
865 | ||
866 | if(folder) { | |
867 | // get name and title from old folder | |
868 | name = folder->GetName(); | |
869 | title = folder->GetTitle(); | |
870 | ||
871 | // delete old one | |
872 | delete folder; | |
873 | ||
874 | // create new one | |
875 | newFolder = CreateFolder(name.Data(),title.Data()); | |
876 | newFolder->SetOwner(); | |
877 | ||
878 | // add objects to folder | |
879 | while(i < size) { | |
880 | newFolder->Add(array->At(i)); | |
881 | i++; | |
882 | } | |
883 | } | |
96b3109a | 884 | |
b4126c69 | 885 | return newFolder; |
96b3109a | 886 | } |
887 | ||
888 | //_____________________________________________________________________________ | |
889 | Long64_t AliComparisonRes::Merge(TCollection* list) | |
890 | { | |
891 | // Merge list of objects (needed by PROOF) | |
892 | ||
893 | if (!list) | |
894 | return 0; | |
895 | ||
896 | if (list->IsEmpty()) | |
897 | return 1; | |
898 | ||
899 | TIterator* iter = list->MakeIterator(); | |
900 | TObject* obj = 0; | |
901 | ||
902 | // collection of generated histograms | |
903 | Int_t count=0; | |
904 | while((obj = iter->Next()) != 0) | |
905 | { | |
906 | AliComparisonRes* entry = dynamic_cast<AliComparisonRes*>(obj); | |
907 | if (entry == 0) continue; | |
908 | ||
909 | fPtResolLPT->Add(entry->fPtResolLPT); | |
910 | fPtResolHPT->Add(entry->fPtResolHPT); | |
911 | fPtPullLPT->Add(entry->fPtPullLPT); | |
912 | fPtPullHPT->Add(entry->fPtPullHPT); | |
decf0997 | 913 | fPhiResolTan->Add(entry->fPhiResolTan); |
914 | fTanResolTan->Add(entry->fTanResolTan); | |
915 | fPhiPullTan->Add(entry->fPhiPullTan); | |
916 | fTanPullTan->Add(entry->fTanPullTan); | |
96b3109a | 917 | |
918 | // Histograms for 1/pt parameterisation | |
3baa4bfd | 919 | f1Pt2ResolS1PtTPC->Add(entry->f1Pt2ResolS1PtTPC); |
920 | fYResolS1PtTPC->Add(entry->fYResolS1PtTPC); | |
921 | fZResolS1PtTPC->Add(entry->fZResolS1PtTPC); | |
922 | fPhiResolS1PtTPC->Add(entry->fPhiResolS1PtTPC); | |
923 | fThetaResolS1PtTPC->Add(entry->fThetaResolS1PtTPC); | |
924 | ||
925 | f1Pt2ResolS1PtTPCITS->Add(entry->f1Pt2ResolS1PtTPCITS); | |
926 | fYResolS1PtTPCITS->Add(entry->fYResolS1PtTPCITS); | |
927 | fZResolS1PtTPCITS->Add(entry->fZResolS1PtTPCITS); | |
928 | fPhiResolS1PtTPCITS->Add(entry->fPhiResolS1PtTPCITS); | |
929 | fThetaResolS1PtTPCITS->Add(entry->fThetaResolS1PtTPCITS); | |
96b3109a | 930 | |
931 | // Resolution histograms (constrained param) | |
932 | fCPhiResolTan->Add(entry->fCPhiResolTan); | |
933 | fCTanResolTan->Add(entry->fCTanResolTan); | |
934 | fCPtResolTan->Add(entry->fCPtResolTan); | |
935 | fCPhiPullTan->Add(entry->fCPhiPullTan); | |
936 | fCTanPullTan->Add(entry->fCTanPullTan); | |
937 | fCPtPullTan->Add(entry->fCPtPullTan); | |
938 | ||
939 | // Histograms for 1/pt parameterisation (constrained) | |
3baa4bfd | 940 | fC1Pt2ResolS1PtTPC->Add(entry->fC1Pt2ResolS1PtTPC); |
941 | fCYResolS1PtTPC->Add(entry->fCYResolS1PtTPC); | |
942 | fCZResolS1PtTPC->Add(entry->fCZResolS1PtTPC); | |
943 | fCPhiResolS1PtTPC->Add(entry->fCPhiResolS1PtTPC); | |
944 | fCThetaResolS1PtTPC->Add(entry->fCThetaResolS1PtTPC); | |
945 | ||
946 | fC1Pt2ResolS1PtTPCITS->Add(entry->fC1Pt2ResolS1PtTPCITS); | |
947 | fCYResolS1PtTPCITS->Add(entry->fCYResolS1PtTPCITS); | |
948 | fCZResolS1PtTPCITS->Add(entry->fCZResolS1PtTPCITS); | |
949 | fCPhiResolS1PtTPCITS->Add(entry->fCPhiResolS1PtTPCITS); | |
950 | fCThetaResolS1PtTPCITS->Add(entry->fCThetaResolS1PtTPCITS); | |
96b3109a | 951 | |
952 | count++; | |
953 | } | |
954 | ||
955 | return count; | |
956 | } | |
3baa4bfd | 957 | |
958 | //_____________________________________________________________________________ | |
959 | TFolder* AliComparisonRes::CreateFolder(TString name,TString title) { | |
960 | // create folder for analysed histograms | |
961 | // | |
962 | TFolder *folder = 0; | |
963 | folder = new TFolder(name.Data(),title.Data()); | |
964 | ||
965 | return folder; | |
966 | } |