]> git.uio.no Git - u/mrichter/AliRoot.git/blame_incremental - HBTAN/AliHBTQResolutionFctns.cxx
git://git.uio.no / u / mrichter / AliRoot.git / blame_incremental
? search:
summary | shortlog | log | commit | commitdiff | tree
blob | blame (non-incremental) | history | HEAD
Updated to the new way AOD stores particles
[u/mrichter/AliRoot.git] / HBTAN / AliHBTQResolutionFctns.cxx
... / ...
CommitLineData
1#include "AliHBTQResolutionFctns.h"
2
3//__________________________________________________________________
4////////////////////////////////////////////////////////////////////
5// //
6// General Remark: //
7// LCMS means //
8// Center of Mass System Longitudially Co-moving //
9// //
10// //
11// This class creates resolution function of Qout //
12// (difference of simulated pair Qout and recontructed pair) //
13// in function of QInv //
14// it inherits from AliHBTTwoPairFctn2D //
15// it needs two pairs to compare //
16// and is two dimentional: numerator and denominator are TH2D //
17// //
18////////////////////////////////////////////////////////////////////
19
20
21
22/******************************************************************/
23/******************************************************************/
24/******************************************************************/
25ClassImp( AliHBTKtResolVsQInvFctn )
26AliHBTKtResolVsQInvFctn::
27AliHBTKtResolVsQInvFctn(Int_t nXbins, Double_t maxXval, Double_t minXval,
28 Int_t nYbins, Double_t maxYval, Double_t minYval):
29 AliHBTTwoPairFctn2D(nXbins,maxXval,minXval,nYbins,maxYval,minYval)
30{
31//ctor
32 Rename("KtResolVsQInv","K_{t} Resolution vs. Q_{Inv}");
33}
34/******************************************************************/
35
36void AliHBTKtResolVsQInvFctn::GetValues(AliHBTPair* trackpair, AliHBTPair* partpair, Double_t& x, Double_t& y) const
37{
38//returns values of the functiion
39 y = partpair->GetKt() - trackpair->GetKt();
40 x = partpair->GetQInv();
41}
42/******************************************************************/
43/******************************************************************/
44/******************************************************************/
45ClassImp( AliHBTQInvResolVsQInvFctn )
46AliHBTQInvResolVsQInvFctn::
47AliHBTQInvResolVsQInvFctn(Int_t nXbins, Double_t maxXval, Double_t minXval,
48 Int_t nYbins, Double_t maxYval, Double_t minYval):
49 AliHBTTwoPairFctn2D(nXbins,maxXval,minXval,nYbins,maxYval,minYval)
50{
51//ctor
52 Rename("QInvResolVsQInv","Q_{Inv} Resolution vs. Q_{Inv}");
53}
54/******************************************************************/
55void AliHBTQInvResolVsQInvFctn::GetValues(AliHBTPair* trackpair, AliHBTPair* partpair, Double_t& x, Double_t& y) const
56{
57//returns values of the functiion
58 y = partpair->GetQInv() - trackpair->GetQInv();
59 x = partpair->GetQInv();
60}
61/******************************************************************/
62/******************************************************************/
63/******************************************************************/
64
65ClassImp( AliHBTQOutResolVsQInvFctn )
66AliHBTQOutResolVsQInvFctn::
67AliHBTQOutResolVsQInvFctn(Int_t nXbins, Double_t maxXval, Double_t minXval,
68 Int_t nYbins, Double_t maxYval, Double_t minYval):
69 AliHBTTwoPairFctn2D(nXbins,maxXval,minXval,nYbins,maxYval,minYval)
70{
71//ctor
72 Rename("QOutResolVsQInv","Q_{Out} Resolution vs. Q_{Inv}");
73}
74/******************************************************************/
75
76void AliHBTQOutResolVsQInvFctn::GetValues(AliHBTPair* trackpair, AliHBTPair* partpair, Double_t& x, Double_t& y) const
77{
78 //returns Qoutsim-Qoutrec for y
79 //returns Qinv for x
80 Double_t tqout = trackpair->GetQOutLCMS();
81 y = partpair->GetQOutLCMS() - tqout;
82 if (tqout < 0.0) y = -y;
83 x = partpair->GetQInv();
84}
85
86/******************************************************************/
87/******************************************************************/
88/******************************************************************/
89
90ClassImp( AliHBTQSideResolVsQInvFctn )
91
92AliHBTQSideResolVsQInvFctn::
93AliHBTQSideResolVsQInvFctn(Int_t nXbins, Double_t maxXval, Double_t minXval,
94 Int_t nYbins, Double_t maxYval, Double_t minYval):
95 AliHBTTwoPairFctn2D(nXbins,maxXval,minXval,nYbins,maxYval,minYval)
96{
97//ctor
98 Rename("QSideResolVsQInv","Q_{Side} Resolution vs. Q_{Inv}");
99}
100/******************************************************************/
101
102void AliHBTQSideResolVsQInvFctn::GetValues(AliHBTPair* trackpair, AliHBTPair* partpair, Double_t& x, Double_t& y) const
103{
104 //returns Qsidesim-Qsiderec for y
105 //returns Qinv for x
106 y = partpair->GetQSideLCMS() - trackpair->GetQSideLCMS();
107 if (trackpair->GetQSideLCMS() < 0.0) y = -y;
108 x = partpair->GetQInv();
109}
110
111/******************************************************************/
112/******************************************************************/
113/******************************************************************/
114
115ClassImp( AliHBTQLongResolVsQInvFctn )
116
117AliHBTQLongResolVsQInvFctn::
118AliHBTQLongResolVsQInvFctn(Int_t nXbins, Double_t maxXval, Double_t minXval,
119 Int_t nYbins, Double_t maxYval, Double_t minYval):
120 AliHBTTwoPairFctn2D(nXbins,maxXval,minXval,nYbins,maxYval,minYval)
121{
122//ctor
123 Rename("QLongResolVsQInv","Q_{Long} Resolution vs. Q_{Inv}");
124}
125/******************************************************************/
126
127void AliHBTQLongResolVsQInvFctn::GetValues(AliHBTPair* trackpair, AliHBTPair* partpair, Double_t& x, Double_t& y) const
128{
129 //returns Qlongsim-Qlongrec for y
130 //returns Qinv for x
131 y = partpair->GetQLongLCMS() - trackpair->GetQLongLCMS();
132 if (trackpair->GetQLongLCMS() < 0.0) y = -y;
133 x = partpair->GetQInv();
134}
135
136/******************************************************************/
137/******************************************************************/
138/******************************************************************/
139
140ClassImp( AliHBTQInvResolVsKtFctn )
141
142AliHBTQInvResolVsKtFctn::
143AliHBTQInvResolVsKtFctn(Int_t nXbins, Double_t maxXval, Double_t minXval,
144 Int_t nYbins, Double_t maxYval, Double_t minYval):
145 AliHBTTwoPairFctn2D(nXbins,maxXval,minXval,nYbins,maxYval,minYval)
146{
147//ctor
148 Rename("QInvResolVsKt","Q_{Inv} Resolution vs. K_{t}");
149}
150/******************************************************************/
151
152void AliHBTQInvResolVsKtFctn::GetValues(AliHBTPair* trackpair, AliHBTPair* partpair, Double_t& x, Double_t& y) const
153{
154 //returns values of the function
155 y = partpair->GetQInv() - trackpair->GetQInv();
156 x = partpair->GetKt();
157}
158
159/******************************************************************/
160/******************************************************************/
161/******************************************************************/
162ClassImp( AliHBTQOutResolVsKtFctn )
163
164AliHBTQOutResolVsKtFctn::
165AliHBTQOutResolVsKtFctn(Int_t nXbins, Double_t maxXval, Double_t minXval,
166 Int_t nYbins, Double_t maxYval, Double_t minYval):
167 AliHBTTwoPairFctn2D(nXbins,maxXval,minXval,nYbins,maxYval,minYval)
168{
169//ctor
170 Rename("QOutResolVsKt","Q_{Out} Resolution vs. K_{t} ");
171}
172/******************************************************************/
173
174void AliHBTQOutResolVsKtFctn::GetValues(AliHBTPair* trackpair, AliHBTPair* partpair, Double_t& x, Double_t& y) const
175{
176 //returns Qoutsim-Qoutrec for y
177 //returns Kt for x
178 y = partpair->GetQOutLCMS() - trackpair->GetQOutLCMS();
179 if (trackpair->GetQOutLCMS() < 0.0) y = -y;
180 x = partpair->GetKt();
181}
182
183/******************************************************************/
184/******************************************************************/
185/******************************************************************/
186ClassImp( AliHBTQSideResolVsKtFctn )
187
188AliHBTQSideResolVsKtFctn::
189AliHBTQSideResolVsKtFctn(Int_t nXbins, Double_t maxXval, Double_t minXval,
190 Int_t nYbins, Double_t maxYval, Double_t minYval):
191 AliHBTTwoPairFctn2D(nXbins,maxXval,minXval,nYbins,maxYval,minYval)
192{
193//ctor
194 Rename("QSideResolVsKt","Q_{Side} Resolution vs. K_{t} ");
195}
196/******************************************************************/
197
198void AliHBTQSideResolVsKtFctn::GetValues(AliHBTPair* trackpair, AliHBTPair* partpair, Double_t& x, Double_t& y) const
199{
200 //returns Qsidesim-Qsiderec for y
201 //returns Kt for x
202 y = partpair->GetQSideLCMS() - trackpair->GetQSideLCMS();
203 if (trackpair->GetQSideLCMS() < 0.0) y = -y;
204 x = partpair->GetKt();
205}
206
207/******************************************************************/
208/******************************************************************/
209/******************************************************************/
210ClassImp( AliHBTQLongResolVsKtFctn )
211
212AliHBTQLongResolVsKtFctn::
213AliHBTQLongResolVsKtFctn(Int_t nXbins, Double_t maxXval, Double_t minXval,
214 Int_t nYbins, Double_t maxYval, Double_t minYval):
215 AliHBTTwoPairFctn2D(nXbins,maxXval,minXval,nYbins,maxYval,minYval)
216{
217//ctor
218 Rename("QLongResolVsKt","Q_{Long} Resolution vs. K_{t} ");
219}
220/******************************************************************/
221
222void AliHBTQLongResolVsKtFctn::GetValues(AliHBTPair* trackpair, AliHBTPair* partpair, Double_t& x, Double_t& y) const
223{
224 //returns Qlongsim-Qlongrec for y
225 //returns Kt for x
226 y = partpair->GetQLongLCMS() - trackpair->GetQLongLCMS();
227 if (trackpair->GetQLongLCMS() < 0.0) y = -y;
228 x = partpair->GetKt();
229}
230/******************************************************************/
231/******************************************************************/
232/******************************************************************/
233
234ClassImp( AliHBTQOutResolVsQOutFctn)
235
236AliHBTQOutResolVsQOutFctn::
237AliHBTQOutResolVsQOutFctn(Int_t nXbins, Double_t maxXval, Double_t minXval,
238 Int_t nYbins, Double_t maxYval, Double_t minYval):
239 AliHBTTwoPairFctn2D(nXbins,maxXval,minXval,nYbins,maxYval,minYval)
240{
241//ctor
242 Rename("QOutResolVsQOut","Q_{Out} Resolution vs. Q_{Out} ");
243}
244/******************************************************************/
245
246void AliHBTQOutResolVsQOutFctn::GetValues(AliHBTPair* trackpair, AliHBTPair* partpair, Double_t& x, Double_t& y) const
247{
248//returns values of the function
249 x = partpair->GetQOutLCMS();
250 y = x - trackpair->GetQOutLCMS();
251}
252
253/******************************************************************/
254/******************************************************************/
255/******************************************************************/
256ClassImp( AliHBTQSideResolVsQSideFctn )
257
258AliHBTQSideResolVsQSideFctn::
259AliHBTQSideResolVsQSideFctn(Int_t nXbins, Double_t maxXval, Double_t minXval,
260 Int_t nYbins, Double_t maxYval, Double_t minYval):
261 AliHBTTwoPairFctn2D(nXbins,maxXval,minXval,nYbins,maxYval,minYval)
262{
263//ctor
264 Rename("QSideResolVsQSide","Q_{Side} Resolution vs. Q_{Side} ");
265}
266/******************************************************************/
267
268void AliHBTQSideResolVsQSideFctn::GetValues(AliHBTPair* trackpair, AliHBTPair* partpair, Double_t& x, Double_t& y) const
269{
270//returns values of the function
271 x = partpair->GetQSideLCMS();
272 y = x - trackpair->GetQSideLCMS();
273}
274
275/******************************************************************/
276/******************************************************************/
277/******************************************************************/
278ClassImp( AliHBTQLongResolVsQLongFctn )
279
280AliHBTQLongResolVsQLongFctn::
281AliHBTQLongResolVsQLongFctn(Int_t nXbins, Double_t maxXval, Double_t minXval,
282 Int_t nYbins, Double_t maxYval, Double_t minYval):
283 AliHBTTwoPairFctn2D(nXbins,maxXval,minXval,nYbins,maxYval,minYval)
284{
285//ctor
286 Rename("QLongResolVsQLong","Q_{Long} Resolution vs. Q_{Long} ");
287}
288/******************************************************************/
289
290void AliHBTQLongResolVsQLongFctn::GetValues(AliHBTPair* trackpair, AliHBTPair* partpair, Double_t& x, Double_t& y) const
291{
292//returns values of the function
293 x = partpair->GetQLongLCMS();
294 y = x - trackpair->GetQLongLCMS();
295}
296
297/******************************************************************/
298/******************************************************************/
299/******************************************************************/
300
301ClassImp( AliHBTPairThetaResolVsQInvFctn )
302
303AliHBTPairThetaResolVsQInvFctn::
304AliHBTPairThetaResolVsQInvFctn(Int_t nXbins, Double_t maxXval, Double_t minXval,
305 Int_t nYbins, Double_t maxYval, Double_t minYval):
306 AliHBTTwoPairFctn2D(nXbins,maxXval,minXval,nYbins,maxYval,minYval)
307{
308//ctor
309 Rename("PairThetaResolVsQInv","Pair Theta Angle Resolution vs. Q_{Inv} ");
310}
311/******************************************************************/
312
313void AliHBTPairThetaResolVsQInvFctn::GetValues(AliHBTPair* trackpair, AliHBTPair* partpair, Double_t& x, Double_t& y) const
314 {
315 //returns Pair Theta sim - Pair Theta rec for y
316 //returns Qinv for x
317 Double_t partTheta = partpair->Particle1()->Theta() - partpair->Particle2()->Theta();
318 Double_t trackTheta = trackpair->Particle1()->Theta() - trackpair->Particle2()->Theta();
319 y = partTheta - trackTheta;
320 x = partpair->GetQInv();
321 }
322/******************************************************************/
323/******************************************************************/
324/******************************************************************/
325
326ClassImp( AliHBTPairThetaResolVsPairThetaFctn )
327
328AliHBTPairThetaResolVsPairThetaFctn::
329AliHBTPairThetaResolVsPairThetaFctn(Int_t nXbins, Double_t maxXval, Double_t minXval,
330 Int_t nYbins, Double_t maxYval, Double_t minYval):
331 AliHBTTwoPairFctn2D(nXbins,maxXval,minXval,nYbins,maxYval,minYval)
332{
333//ctor
334 Rename("PairThetaResolVsPairTheta","Pair Theta Angle Resolution vs. Pair Theta ");
335}
336/******************************************************************/
337void AliHBTPairThetaResolVsPairThetaFctn::GetValues(AliHBTPair* trackpair, AliHBTPair* partpair, Double_t& x, Double_t& y) const
338{
339 //returns Pair Theta sim - Pair Theta rec for y
340 //returns Pair Theta sim for x
341 Double_t partTheta = partpair->Particle1()->Theta() - partpair->Particle2()->Theta();
342 Double_t trackTheta = trackpair->Particle1()->Theta() - trackpair->Particle2()->Theta();
343 y = partTheta - trackTheta;
344 if (trackTheta < 0.0) y = -y;
345 x = trackTheta;
346}
347/******************************************************************/
348/******************************************************************/
349/******************************************************************/
350
351ClassImp( AliHBTPairPhiResolVsQInvFctn )
352
353AliHBTPairPhiResolVsQInvFctn::
354AliHBTPairPhiResolVsQInvFctn(Int_t nXbins, Double_t maxXval, Double_t minXval,
355 Int_t nYbins, Double_t maxYval, Double_t minYval):
356 AliHBTTwoPairFctn2D(nXbins,maxXval,minXval,nYbins,maxYval,minYval)
357{
358//ctor
359 Rename("PairPhiResolVsQInv","Pair Phi Angle Resolution vs. Q_{Inv} ");
360}
361/******************************************************************/
362
363void AliHBTPairPhiResolVsQInvFctn::GetValues(AliHBTPair* trackpair, AliHBTPair* partpair, Double_t& x, Double_t& y) const
364 {
365 //returns Pair Phi sim - Pair Phi rec for y
366 //returns QInv sim for x
367 Double_t partPhi = partpair->Particle1()->Phi() - partpair->Particle2()->Phi();
368 Double_t trackPhi = trackpair->Particle1()->Phi() - trackpair->Particle2()->Phi();
369 y = partPhi - trackPhi;
370 x = partpair->GetQInv();
371 }
372
373/******************************************************************/
374/******************************************************************/
375/******************************************************************/
376
377
378ClassImp( AliHBTPairThetaResolVsKtFctn )
379
380AliHBTPairThetaResolVsKtFctn::
381AliHBTPairThetaResolVsKtFctn(Int_t nXbins, Double_t maxXval, Double_t minXval,
382 Int_t nYbins, Double_t maxYval, Double_t minYval):
383 AliHBTTwoPairFctn2D(nXbins,maxXval,minXval,nYbins,maxYval,minYval)
384{
385//ctor
386 Rename("PairThetaResolVsKt","Pair Theta Angle Resolution vs. K_{t} ");
387}
388/******************************************************************/
389
390void AliHBTPairThetaResolVsKtFctn::GetValues(AliHBTPair* trackpair, AliHBTPair* partpair, Double_t& x, Double_t& y) const
391{
392 //returns Pair Theta sim - Pair Theta rec for y
393 //returns Kt sim for x
394 Double_t partTheta = partpair->Particle1()->Theta() - partpair->Particle2()->Theta();
395 Double_t trackTheta = trackpair->Particle1()->Theta() - trackpair->Particle2()->Theta();
396 y = partTheta - trackTheta;
397 x = partpair->GetKt();
398}
399/******************************************************************/
400/******************************************************************/
401/******************************************************************/
402
403ClassImp( AliHBTPairPhiResolVsKtFctn )
404
405AliHBTPairPhiResolVsKtFctn::
406AliHBTPairPhiResolVsKtFctn(Int_t nXbins, Double_t maxXval, Double_t minXval,
407 Int_t nYbins, Double_t maxYval, Double_t minYval):
408 AliHBTTwoPairFctn2D(nXbins,maxXval,minXval,nYbins,maxYval,minYval)
409{
410//ctor
411 Rename("PairPhiResolVsKt","Pair Phi Angle Resolution vs. K_{t} ");
412}
413/******************************************************************/
414
415void AliHBTPairPhiResolVsKtFctn::GetValues(AliHBTPair* trackpair, AliHBTPair* partpair, Double_t& x, Double_t& y) const
416{
417 //returns Pair Phi sim - Pair Phi rec for y
418 //returns Kt sim for x
419 Double_t partPhi = partpair->Particle1()->Phi() - partpair->Particle2()->Phi();
420 Double_t trackPhi = trackpair->Particle1()->Phi() - trackpair->Particle2()->Phi();
421 y = partPhi - trackPhi;
422 x = partpair->GetKt();
423}
424
425/******************************************************************/
426/******************************************************************/
427/******************************************************************/
428
429
430ClassImp( AliHBTPairPhiResolVsPairPhiFctn )
431
432AliHBTPairPhiResolVsPairPhiFctn::
433AliHBTPairPhiResolVsPairPhiFctn(Int_t nXbins, Double_t maxXval, Double_t minXval,
434 Int_t nYbins, Double_t maxYval, Double_t minYval):
435 AliHBTTwoPairFctn2D(nXbins,maxXval,minXval,nYbins,maxYval,minYval)
436{
437//ctor
438 Rename("PairPhiResolVsPairPhi","Pair Phi Angle Resolution vs. Pair Phi ");
439}
440/******************************************************************/
441
442void AliHBTPairPhiResolVsPairPhiFctn::GetValues(AliHBTPair* trackpair, AliHBTPair* partpair, Double_t& x, Double_t& y) const
443{
444 //returns Pair Phi sim - Pair Phi rec for y
445 //returns Pair Phi sim for x
446 Double_t partPhi = partpair->Particle1()->Phi() - partpair->Particle2()->Phi();
447 Double_t trackPhi = trackpair->Particle1()->Phi() - trackpair->Particle2()->Phi();
448 y = partPhi - trackPhi;
449 x = trackPhi;
450}
451
452/******************************************************************/
453/******************************************************************/
454/******************************************************************/