In mapping:
[u/mrichter/AliRoot.git] / PWG1 / AliAnalysisTaskQASym.cxx
CommitLineData
9f9631f3 1#include "TChain.h"
2#include "TTree.h"
3#include "TH1F.h"
4#include "TH2F.h"
1d55ca60 5#include "TH3F.h"
9f9631f3 6#include "TCanvas.h"
7#include "TList.h"
8#include "TParticle.h"
9#include "TParticlePDG.h"
10#include "TProfile.h"
11#include "TNtuple.h"
12#include "TFile.h"
13
14#include "AliAnalysisTask.h"
15#include "AliAnalysisManager.h"
16
17#include "AliESDEvent.h"
e4ab7356 18#include "AliLog.h"
9f9631f3 19#include "AliESDVertex.h"
20#include "AliESDInputHandler.h"
21#include "AliESDtrackCuts.h"
22#include "AliMultiplicity.h"
23
9f9631f3 24#include "AliAnalysisTaskQASym.h"
25#include "AliExternalTrackParam.h"
26#include "AliTrackReference.h"
9f9631f3 27#include "AliHeader.h"
28#include "AliGenEventHeader.h"
29#include "AliGenDPMjetEventHeader.h"
30
1d55ca60 31// Analysis Task for basic QA on the ESD
32
33// Authors: Jan Fiete Grosse-Oetringhaus, Christian Klein-Boesing,
34// Andreas Morsch, Eva Sicking
35
9f9631f3 36ClassImp(AliAnalysisTaskQASym)
37
1d55ca60 38 //________________________________________________________________________
39 AliAnalysisTaskQASym::AliAnalysisTaskQASym(const char *name)
40 : AliAnalysisTaskSE(name)
41 ,fTrackType(0)
2efff183 42 ,fFieldOn(kTRUE)
43 ,fHists(0)
44 ,fHistRECpt(0)
45 ,fEta(0)
46 ,fEtaPhi(0)
47 ,fEtaPt(0)
48 ,fQPt(0)
49 ,fDca(0)
50 ,fqRec(0)
51 ,fsigmaPt(0)
cbac53a7 52
2efff183 53 ,fRecPtPos(0)
54 ,fRecPtNeg(0)
55 ,fRecPhiPos(0)
56 ,fRecPhiNeg(0)
57 ,fRecEtaPos(0)
58 ,fRecEtaNeg(0)
59 ,fRecEtaPtPos(0)
60 ,fRecEtaPtNeg(0)
61 ,fRecDcaPos(0)
62 ,fRecDcaNeg(0)
63 ,fRecDcaNegInv(0)
64 ,fRecDPos(0)
65 ,fRecDNeg(0)
cbac53a7 66
67
2efff183 68 ,fRecQPtPosEta(0)
69 ,fRecQPtNegEta(0)
70 ,fRecPtPosEta(0)
71 ,fRecPtNegEta(0)
72 ,fRecPhiPosEta(0)
73 ,fRecPhiNegEta(0)
74 ,fRecDcaPosEta(0)
75 ,fRecDcaNegEta(0)
76 ,fRecDPosEta(0)
77 ,fRecDNegEta(0)
cbac53a7 78
2efff183 79 ,fRecPtPosVz(0)
80 ,fRecPtNegVz(0)
81 ,fRecEtaPosVz(0)
82 ,fRecEtaNegVz(0)
83 ,fRecPhiPosVz(0)
84 ,fRecPhiNegVz(0)
85 ,fSignedDcaPosVz(0)
86 ,fSignedDcaNegVz(0)
87 ,fRecQPtPosEtaVz(0)
88 ,fRecQPtNegEtaVz(0)
89 ,fRecEtaPtPosVz(0)
90 ,fRecEtaPtNegVz(0)
cbac53a7 91
92
2efff183 93 ,fDeltaPhiAll(0)
94 ,fDeltaPhiLeading(0)
95 ,fDiffDcaD(0)
cbac53a7 96
2efff183 97 ,fPhiRec(0)
98 ,fThetaRec(0)
99 ,fNumber(0)
100 ,fVx(0)
101 ,fVy(0)
102 ,fVz(0)
cbac53a7 103 ,fVertexX(0)
104 ,fVertexY(0)
105 ,fVertexZ(0)
1d55ca60 106 ,test(0)
9f9631f3 107
1d55ca60 108 ,fRecDcaPosPhi(0)
109 ,fRecDcaNegPhi(0)
110 ,fRecPtPosPhi(0)
111 ,fRecPtNegPhi(0)
112 ,fRecEtaPosPhi(0)
113 ,fRecEtaNegPhi(0)
114 ,fRecQPtPhi(0)
115 ,fRecEtaPtPosPhi(0)
116 ,fRecEtaPtNegPhi(0)
117
118 ,fRecPtPosEtaPos(0)
119 ,fRecPtNegEtaPos(0)
120 ,fRecPtPosEtaNeg(0)
121 ,fRecPtNegEtaNeg(0)
122
123 ,fRec1PtPosEtaPos(0)
124 ,fRec1PtNegEtaPos(0)
125 ,fRec1PtPosEtaNeg(0)
126 ,fRec1PtNegEtaNeg(0)
127
128 ,fRecPhiPosEtaPos(0)
129 ,fRecPhiNegEtaPos(0)
130 ,fRecPhiPosEtaNeg(0)
131 ,fRecPhiNegEtaNeg(0)
132
133 ,fRecDcaPosPhiEtaPos(0)
134 ,fRecDcaNegPhiEtaPos(0)
135 ,fRecDcaPosPhiEtaNeg(0)
136 ,fRecDcaNegPhiEtaNeg(0)
98a5f582 137
138 ,fRecDcaPosPtEtaPos(0)
139 ,fRecDcaNegPtEtaPos(0)
140 ,fRecDcaPosPtEtaNeg(0)
141 ,fRecDcaNegPtEtaNeg(0)
142
1d55ca60 143 ,fRecPtPosPhiEtaPos(0)
144 ,fRecPtNegPhiEtaPos(0)
145 ,fRecPtPosPhiEtaNeg(0)
146 ,fRecPtNegPhiEtaNeg(0)
147
148
8ae6f5d0 149// ,fRecDcaPhiPtPosEtaPos(0)
150// ,fRecDcaPhiPtNegEtaPos(0)
151// ,fRecDcaPhiPtPosEtaNeg(0)
152// ,fRecDcaPhiPtNegEtaNeg(0)
153
1d55ca60 154 ,fEtavPt(0)
98a5f582 155 ,fCompareTPCparam(0)
9f9631f3 156
1d55ca60 157
158 ,sdca(0)
159 ,xy(0)
160 ,z(0)
161 ,xvertexcor(0)
162 ,yvertexcor(0)
163
9f9631f3 164 ,fCuts(0)
1d55ca60 165
9f9631f3 166{
1d55ca60 167 // Constructor
9f9631f3 168 for(Int_t i = 0;i<18;++i){
169 fRecPtTpcSector[i] = 0;
170 fRecEtaTpcSector[i] = 0;
171 fSignedDcaTpcSector[i] = 0;
172 fRecQPtTpcSector[i] = 0;
173 fRecEtaPtTpcSector[i] = 0;
174 }
175
176 for(Int_t i = 0;i< 7;++i){
177 fRecPtPosLadder[i] = 0;
178 fRecPtNegLadder[i] = 0;
179 fRecPhiPosLadder[i] = 0;
180 fRecPhiNegLadder[i] = 0;
181 fRecEtaPosLadder[i] = 0;
182 fRecEtaNegLadder[i] = 0;
183 fSignDcaPos[i] = 0;
184 fSignDcaNeg[i] = 0;
185 fSignDcaNegInv[i] = 0;
186 fPtSigmaPos[i] =0;
187 fPtSigmaNeg[i] =0;
188 fqPtRec[i] =0;
189 fDcaSigmaPos[i] =0;
190 fDcaSigmaNeg[i] =0;
191 }
192
193 DefineOutput(1, TList::Class());
194
195
196
197}
198
199
200//________________________________________________________________________
201void AliAnalysisTaskQASym::UserCreateOutputObjects()
202{
203 // Create histograms
204 // Called once
205
1d55ca60 206 Bool_t oldStatus = TH1::AddDirectoryStatus();
207 TH1::AddDirectory(kFALSE);
208
8ae6f5d0 209 Double_t range = 0.3;
9f9631f3 210 Double_t pt = 20.;
211
212 fHists = new TList();
1d55ca60 213 // test = new TNtuple("test","test",
214 // "pt:phi:theta:x:y:z:charge");
9f9631f3 215 fHistRECpt = new TH1F("fHistRECpt",
216 " p_{T}",
cbac53a7 217 200, 0., pt);
9f9631f3 218 fEta = new TH1F("fEta",
1d55ca60 219 " #eta",
220 200, -2., 2.);
98a5f582 221 fEtavPt = new TH2F("fEtavPt",
222 " #eta -p_{T}",
223 200, -2., 2.,
224 100, 0, 1.5);
225 fCompareTPCparam = new TH2F("fCompareTPCparam",
226 "fCompareTPCparam",
227 100, -1., 1.,100,-5, 5);
228
e4ab7356 229 fEtaPhi = new TH2F("fEtaPhi",
1d55ca60 230 " #eta - #phi",
231 200, -2., 2., 128, 0., 2. * TMath::Pi());
e4ab7356 232
9f9631f3 233 fThetaRec = new TH1F("fThetaRec",
1d55ca60 234 " #theta",
235 180, 0., TMath::Pi());
9f9631f3 236 fPhiRec = new TH1F("fPhiRec",
1d55ca60 237 " #phi",
238 180, 0., 2*TMath::Pi());
9f9631f3 239 fNumber = new TH1F("fNumber",
240 "number of tracks per event",
cbac53a7 241 200, -0.5, 199.5);
9f9631f3 242 fVx = new TH1F("fVx",
cbac53a7 243 "X of first track point",
1d55ca60 244 100, -1., 1.);
9f9631f3 245 fVy = new TH1F("fVy",
cbac53a7 246 "Y of first track point",
1d55ca60 247 100, -1., 1.);
9f9631f3 248 fVz = new TH1F("fVz",
cbac53a7 249 "Z of first track point",
1d55ca60 250 200, -50., 50.);
cbac53a7 251 fVertexX = new TH1F("fVerteX",
252 "X of vertex",
253 100, -1., 1.);
254 fVertexY = new TH1F("fVertexY",
255 "Y of vertex",
256 100, -1., 1.);
257 fVertexZ = new TH1F("fVertexZ",
258 "Z of vertex",
259 200, -50., 50.);
260
9f9631f3 261 fEtaPt = new TH1F("fEtaPt",
1d55ca60 262 " #eta/p_{T} ",
263 100, -1., 1.);
9f9631f3 264
265 fQPt = new TH1F("fQPt",
1d55ca60 266 " charge/p_{T} ",
267 100, -1., 1.);
9f9631f3 268
269 fDca = new TH1F("fDca",
1d55ca60 270 " dca ",
271 200, -range*(1+Int_t(fTrackType/2)*9), range*(1+Int_t(fTrackType/2)*9));
9f9631f3 272
273
274 fqRec = new TH1F("fqRec",
1d55ca60 275 " charge all reconstructed particle",
276 21, -9.5, 10.5);
9f9631f3 277
278 fsigmaPt = new TH1F("fsigmaPt",
1d55ca60 279 "Log_{10}(#sigma_{p_{T}})",
280 200, -4., 8.);
9f9631f3 281
282
283
284
285 //------------
286 for(Int_t ITSlayer_case=0;ITSlayer_case<7;ITSlayer_case++){
287
288 fSignDcaPos[ITSlayer_case] = new TH1F(Form("fSignDcaPos%d", ITSlayer_case),
1d55ca60 289 " Signed dca",
290 200, -range*(1+Int_t(fTrackType/2)*9), range*(1+Int_t(fTrackType/2)*9));
9f9631f3 291 fSignDcaPos[ITSlayer_case]->GetXaxis()->SetTitle("dca");
292 fSignDcaPos[ITSlayer_case]->GetYaxis()->SetTitle("");
293
294
295 fSignDcaNeg[ITSlayer_case] = new TH1F(Form("fSignDcaNeg%d", ITSlayer_case),
1d55ca60 296 " Signed dcas",
297 200, -range*(1+Int_t(fTrackType/2)*9), range*(1+Int_t(fTrackType/2)*9));
9f9631f3 298 fSignDcaNeg[ITSlayer_case]->GetXaxis()->SetTitle("dca");
299 fSignDcaNeg[ITSlayer_case]->GetYaxis()->SetTitle("");
300
301 fSignDcaNegInv[ITSlayer_case] = new TH1F(Form("fSignDcaNegInv%d", ITSlayer_case),
1d55ca60 302 " inverse Signed dca ",
303 200, -range*(1+Int_t(fTrackType/2)*9), range*(1+Int_t(fTrackType/2)*9));
9f9631f3 304 fSignDcaNegInv[ITSlayer_case]->GetXaxis()->SetTitle("-dca");
305 fSignDcaNegInv[ITSlayer_case]->GetYaxis()->SetTitle("");
306
307
308
309
310 fPtSigmaPos[ITSlayer_case] = new TH1F(Form("fPtSigmaPos%d", ITSlayer_case),
1d55ca60 311 " #sigma_{pT} ",
312 208, -4., 8.);
9f9631f3 313 fPtSigmaPos[ITSlayer_case]->GetXaxis()->SetTitle("Log_{10}(#sigma_{pT})");
314 fPtSigmaPos[ITSlayer_case]->GetYaxis()->SetTitle("");
315
316
317 fPtSigmaNeg[ITSlayer_case] = new TH1F(Form("fPtSigmaNeg%d",ITSlayer_case),
1d55ca60 318 " #sigma_{pT}",
319 208, -4., 8.);
9f9631f3 320 fPtSigmaNeg[ITSlayer_case]->GetXaxis()->SetTitle("Log_{10}(#sigma_{pT})");
321 fPtSigmaNeg[ITSlayer_case]->GetYaxis()->SetTitle("");
322
323
324
325
326
327 fqPtRec[ITSlayer_case] = new TH1F(Form("fqPtRec%d",ITSlayer_case),
1d55ca60 328 "q/ p_{T}",
329 200, -100., 100.);
9f9631f3 330 fqPtRec[ITSlayer_case]->GetXaxis()->SetTitle("q_{tr}/p_{T, tr} (GeV/c)");
331 fqPtRec[ITSlayer_case]->GetYaxis()->SetTitle("");
332
333
334
335
336
337 fDcaSigmaPos[ITSlayer_case] = new TH2F(Form("fDcaSigmaPos%d", ITSlayer_case),
1d55ca60 338 " p_{T} shift vs #sigma_{pT} ",
339 200, -range*(1+Int_t(fTrackType/2)*9), range*(1+Int_t(fTrackType/2)*9),200, -4., 4. );
9f9631f3 340 fDcaSigmaPos[ITSlayer_case]->GetXaxis()->SetTitle("signed DCA)");
341 fDcaSigmaPos[ITSlayer_case]->GetYaxis()->SetTitle("log_{10}(#sigma_{pT})");
342
343
344 fDcaSigmaNeg[ITSlayer_case] = new TH2F(Form("fDcaSigmaNeg%d", ITSlayer_case),
1d55ca60 345 " p_{T} shift vs #sigma_{pT} ",
346 200, -range*(1+Int_t(fTrackType/2)*9), range*(1+Int_t(fTrackType/2)*9),200, -4., 4. );
9f9631f3 347 fDcaSigmaNeg[ITSlayer_case]->GetXaxis()->SetTitle("signed DCA");
348 fDcaSigmaNeg[ITSlayer_case]->GetYaxis()->SetTitle("log_{10}(#sigma_{pT})");
1d55ca60 349
350
2efff183 351 }
9f9631f3 352
353
354
1d55ca60 355 // YIELDs---------- positive and negative particles
9f9631f3 356
1d55ca60 357 fRecPtPos = new TH1F("fRecPtPos",
358 " p_{T}",
359 100, 0.,pt);
360 fRecPtPos->GetXaxis()->SetTitle("p_{T} (GeV/c)");
361 fRecPtNeg = new TH1F("fRecPtNeg",
362 " p_{T} ",
363 100, 0., pt);
364 fRecPtNeg->GetXaxis()->SetTitle("p_{T} (GeV/c)");
9f9631f3 365
366
1d55ca60 367 fRecPhiPos = new TH1F("fRecPhiPos",
368 "#phi",
369 361, 0., 360.);
370 fRecPhiPos->GetXaxis()->SetTitle("#phi (deg)");
9f9631f3 371
1d55ca60 372 fRecPhiNeg = new TH1F("fRecPhiNeg",
373 "#phi ",
374 361, 0., 360.);
375 fRecPhiNeg->GetXaxis()->SetTitle("#phi (deg)");
9f9631f3 376
1d55ca60 377 fRecEtaPos = new TH1F("fRecEtaPos",
378 "#eta",
379 200, -2., 2.);
380 fRecEtaPos->GetXaxis()->SetTitle("#eta");
9f9631f3 381
1d55ca60 382 fRecEtaNeg = new TH1F("fRecEtaNeg",
383 "#eta",
384 200, -2., 2.);
385 fRecEtaNeg->GetXaxis()->SetTitle("#eta");
386
387 fRecEtaPtPos = new TH1F("fRecEtaPtPos",
388 "#eta/p_{T}",
389 200, -0.1, .1);
390 fRecEtaPtPos->GetXaxis()->SetTitle("#eta/p_{T}");
391
392 fRecEtaPtNeg = new TH1F("fRecEtaPtNeg",
393 "#eta/p_{T}",
394 200, -.1, .1);
395 fRecEtaPtNeg->GetXaxis()->SetTitle("#eta/p_{T}");
396
397 fRecDcaPos = new TH1F("fRecDcaPos",
398 " dca",
399 100, -range*(1+Int_t(fTrackType/2)*9), range*(1+Int_t(fTrackType/2)*9));
400 fRecDcaPos->GetXaxis()->SetTitle("dca (cm)");
401 fRecDcaNeg = new TH1F("fRecDcaNeg",
402 " dca",
403 100, -range*(1+Int_t(fTrackType/2)*9), range*(1+Int_t(fTrackType/2)*9));
404 fRecDcaNeg->GetXaxis()->SetTitle("dca (cm)");
405
406 fRecDcaNegInv = new TH1F("fRecDcaNegInv",
407 " dca",
408 100, -range*(1+Int_t(fTrackType/2)*9), range*(1+Int_t(fTrackType/2)*9));
409 fRecDcaNegInv->GetXaxis()->SetTitle("dca (cm)");
410
411
412 fRecDPos = new TH1F("fRecDPos",
413 " d",
414 100, -range*(1+Int_t(fTrackType/2)*9), range*(1+Int_t(fTrackType/2)*9));
415 fRecDPos->GetXaxis()->SetTitle("d (cm)");
416 fRecDNeg = new TH1F("fRecDNeg",
417 "d",
418 100, -range*(1+Int_t(fTrackType/2)*9), range*(1+Int_t(fTrackType/2)*9));
419 fRecDNeg->GetXaxis()->SetTitle("d (cm)");
420
421
422 // YIELDs ---------------- positive and negative eta
9f9631f3 423
424
1d55ca60 425 fRecQPtPosEta = new TH1F("fRecQPtPosEta",
426 "q/p_{T}",
427 200, -0.5, 0.5);
428 fRecQPtPosEta->GetXaxis()->SetTitle("q/p_{T} ");
429
430 fRecQPtNegEta = new TH1F("fRecQPtNegEta",
431 "q/p_{T}",
432 200, -0.5, 0.5);
433 fRecQPtNegEta->GetXaxis()->SetTitle("q/p_{T}");
9f9631f3 434
1d55ca60 435 fRecPtPosEta = new TH1F("fRecPtPosEta",
436 " p_{T} ",
437 100, 0., pt);
438 fRecPtPosEta->GetXaxis()->SetTitle("p_{T} (GeV/c)");
439
440 fRecPtNegEta = new TH1F("fRecPtNegEta",
441 " p_{T}",
442 100, 0., pt);
443 fRecPtNegEta->GetXaxis()->SetTitle("p_{T} (GeV/c)");
9f9631f3 444
1d55ca60 445 fRecPhiPosEta = new TH1F("fRecPhiPosEta",
446 "#phi",
447 361, 0., 360);
448 fRecPhiPosEta->GetXaxis()->SetTitle("#phi (deg)");
449
450 fRecPhiNegEta = new TH1F("fRecPhiNegEta",
451 "#phi ",
452 361, 0, 360);
453 fRecPhiNegEta->GetXaxis()->SetTitle("#phi (deg)");
454
455 fRecDcaPosEta = new TH1F("fRecDcaPosEta",
456 " dca ",
457 100, -range*(1+Int_t(fTrackType/2)*9), range*(1+Int_t(fTrackType/2)*9));
458 fRecDcaPosEta->GetXaxis()->SetTitle("dca (cm)");
459 fRecDcaNegEta = new TH1F("fRecDcaNegEta",
460 " dca",
461 100, -range*(1+Int_t(fTrackType/2)*9), range*(1+Int_t(fTrackType/2)*9));
462 fRecDcaNegEta->GetXaxis()->SetTitle("dca (cm)");
463
464 fRecDPosEta = new TH1F("fRecDPosEta",
465 " d",
466 100, -range*(1+Int_t(fTrackType/2)*9), range*(1+Int_t(fTrackType/2)*9));
467 fRecDPosEta->GetXaxis()->SetTitle("d (cm)");
468 fRecDNegEta = new TH1F("fRecDNegEta",
9f9631f3 469 "d",
470 100, -5., 5.);
1d55ca60 471 fRecDNegEta->GetXaxis()->SetTitle("d (cm)");
472
473 fRecDcaPosPhi = new TH2F("fRecDcaPosPhi",
474 " dca vs. phi",
475 100, -range*(1+Int_t(fTrackType/2)*9), range*(1+Int_t(fTrackType/2)*9), 180, 0, TMath::Pi()*2);
476 fRecDcaPosPhi->GetXaxis()->SetTitle("dca (cm)");
477 fRecDcaPosPhi->GetYaxis()->SetTitle("#phi (rad.)");
478 fRecDcaNegPhi = new TH2F("fRecDcaNegPhi",
479 " dca vs. phi",
480 100, -range*(1+Int_t(fTrackType/2)*9), range*(1+Int_t(fTrackType/2)*9), 180, 0, TMath::Pi()*2);
481 fRecDcaNegPhi->GetXaxis()->SetTitle("dca (cm)");
482 fRecDcaNegPhi->GetYaxis()->SetTitle("#phi (rad.)");
483
484 fRecPtPosPhi = new TH2F("fRecPtPosPhi",
485 " log(p_T) vs. phi",
486 100, -2.5, 2., 180, 0, TMath::Pi()*2);
487 fRecPtPosPhi->GetXaxis()->SetTitle("log_{10}(p_{T})");
488 fRecPtPosPhi->GetYaxis()->SetTitle("#phi (rad.)");
489 fRecPtNegPhi = new TH2F("fRecPtNegPhi",
490 " log(p_T) vs. phi",
491 100,-2.5 , 2., 180, 0, TMath::Pi()*2);
492 fRecPtNegPhi->GetXaxis()->SetTitle("log_{10}(p_{T})");
493 fRecPtNegPhi->GetYaxis()->SetTitle("#phi (rad.)");
494
495 fRecEtaPosPhi = new TH2F("fRecEtaPosPhi",
496 " eta vs. phi",
497 100, -1.5, 1.5, 180, 0, TMath::Pi()*2);
498 fRecEtaPosPhi->GetXaxis()->SetTitle("#eta");
499 fRecEtaPosPhi->GetYaxis()->SetTitle("#phi (rad.)");
500 fRecEtaNegPhi = new TH2F("fRecEtaNegPhi",
501 " eta vs. phi",
502 100, -1.5, 1.5, 180, 0, TMath::Pi()*2);
503 fRecEtaNegPhi->GetXaxis()->SetTitle("#eta");
504 fRecEtaNegPhi->GetYaxis()->SetTitle("#phi (rad.)");
505
506 fRecQPtPhi = new TH2F("fRecQPtPhi",
507 " charge/p_T vs. phi",
508 100,-1. , 1., 180, 0, TMath::Pi()*2);
509 fRecQPtPhi->GetXaxis()->SetTitle("charge/p_{T}");
510 fRecQPtPhi->GetYaxis()->SetTitle("#phi (rad.)");
511
512 fRecEtaPtPosPhi = new TH2F("fRecEtaPtPosPhi",
513 " eta/p_T vs. phi",
514 100, -5, 5., 180, 0, TMath::Pi()*2);
515 fRecEtaPtPosPhi->GetXaxis()->SetTitle("#eta/p_{T}");
516 fRecEtaPtPosPhi->GetYaxis()->SetTitle("#phi (rad.)");
517 fRecEtaPtNegPhi = new TH2F("fRecEtaPtNegPhi",
518 " eta/p_T vs. phi",
519 100,-5 , 5., 180, 0, TMath::Pi()*2);
520 fRecEtaPtNegPhi->GetXaxis()->SetTitle("#eta/p_{T}");
521 fRecEtaPtNegPhi->GetYaxis()->SetTitle("#phi (rad.)");
522
523
524
525
526
527 fRecDcaPosPhiEtaPos = new TH2F("fRecDcaPosPhiEtaPos",
528 " dca vs. phi",
529 100, -range*(1+Int_t(fTrackType/2)*9), range*(1+Int_t(fTrackType/2)*9), 180, 0, TMath::Pi()*2);
530 fRecDcaPosPhiEtaPos->GetXaxis()->SetTitle("dca (cm)");
531 fRecDcaPosPhiEtaPos->GetYaxis()->SetTitle("#phi (rad.)");
532 fRecDcaNegPhiEtaPos = new TH2F("fRecDcaNegPhiEtaPos",
533 " dca vs. phi",
534 100, -range*(1+Int_t(fTrackType/2)*9), range*(1+Int_t(fTrackType/2)*9), 180, 0, TMath::Pi()*2);
535 fRecDcaNegPhiEtaPos->GetXaxis()->SetTitle("dca (cm)");
536 fRecDcaNegPhiEtaPos->GetYaxis()->SetTitle("#phi (rad.)");
537
538 fRecPtPosPhiEtaPos = new TH2F("fRecPtPosPhiEtaPos",
539 " log(p_T) vs. phi",
540 100, -2.5, 2., 180, 0, TMath::Pi()*2);
541 fRecPtPosPhiEtaPos->GetXaxis()->SetTitle("log_{10}(p_{T})");
542 fRecPtPosPhiEtaPos->GetYaxis()->SetTitle("#phi (rad.)");
543 fRecPtNegPhiEtaPos = new TH2F("fRecPtNegPhiEtaPos",
544 " log(p_T) vs. phi",
545 100,-2.5 , 2., 180, 0, TMath::Pi()*2);
546 fRecPtNegPhiEtaPos->GetXaxis()->SetTitle("log_{10}(p_{T})");
547 fRecPtNegPhiEtaPos->GetYaxis()->SetTitle("#phi (rad.)");
548
549
550 fRecDcaPosPhiEtaNeg = new TH2F("fRecDcaPosPhiEtaNeg",
551 " dca vs. phi",
552 100, -range*(1+Int_t(fTrackType/2)*9), range*(1+Int_t(fTrackType/2)*9), 180, 0, TMath::Pi()*2);
553 fRecDcaPosPhiEtaNeg->GetXaxis()->SetTitle("dca (cm)");
554 fRecDcaPosPhiEtaNeg->GetYaxis()->SetTitle("#phi (rad.)");
555 fRecDcaNegPhiEtaNeg = new TH2F("fRecDcaNegPhiEtaNeg",
556 " dca vs. phi",
557 100, -range*(1+Int_t(fTrackType/2)*9), range*(1+Int_t(fTrackType/2)*9), 180, 0, TMath::Pi()*2);
558 fRecDcaNegPhiEtaNeg->GetXaxis()->SetTitle("dca (cm)");
559 fRecDcaNegPhiEtaNeg->GetYaxis()->SetTitle("#phi (rad.)");
560
561 fRecPtPosPhiEtaNeg = new TH2F("fRecPtPosPhiEtaNeg",
562 " log(p_T) vs. phi",
563 100, -2.5, 2., 180, 0, TMath::Pi()*2);
564 fRecPtPosPhiEtaNeg->GetXaxis()->SetTitle("log_{10}(p_{T})");
565 fRecPtPosPhiEtaNeg->GetYaxis()->SetTitle("#phi (rad.)");
566 fRecPtNegPhiEtaNeg = new TH2F("fRecPtNegPhiEtaNeg",
567 " log(p_T) vs. phi",
568 100,-2.5 , 2., 180, 0, TMath::Pi()*2);
569 fRecPtNegPhiEtaNeg->GetXaxis()->SetTitle("log_{10}(p_{T})");
570 fRecPtNegPhiEtaNeg->GetYaxis()->SetTitle("#phi (rad.)");
98a5f582 571
572 //new
573 fRecDcaPosPtEtaPos = new TH2F("fRecDcaPosPtEtaPos",
574 " dca vs. pt",
575 100, -range*(1+Int_t(fTrackType/2)*9), range*(1+Int_t(fTrackType/2)*9), 200, -2, 2);
576 fRecDcaPosPtEtaPos->GetXaxis()->SetTitle("dca (cm)");
577 fRecDcaPosPtEtaPos->GetYaxis()->SetTitle("log_{10}(p_{T})");
578
579 fRecDcaPosPtEtaNeg = new TH2F("fRecDcaPosPtEtaNeg",
580 " dca vs. pt",
581 100, -range*(1+Int_t(fTrackType/2)*9), range*(1+Int_t(fTrackType/2)*9), 200, -2, 2);
582 fRecDcaPosPtEtaNeg->GetXaxis()->SetTitle("dca (cm)");
583 fRecDcaPosPtEtaNeg->GetYaxis()->SetTitle("log_{10}(p_{T})");
584
585 fRecDcaNegPtEtaPos = new TH2F("fRecDcaNegPtEtaPos",
586 " dca vs. pt",
587 100, -range*(1+Int_t(fTrackType/2)*9), range*(1+Int_t(fTrackType/2)*9), 200, -2, 2);
588 fRecDcaNegPtEtaPos->GetXaxis()->SetTitle("dca (cm)");
589 fRecDcaNegPtEtaPos->GetYaxis()->SetTitle("log_{10}(p_{T})");
590
591 fRecDcaNegPtEtaNeg = new TH2F("fRecDcaNegPtEtaNeg",
592 " dca vs. pt",
593 100, -range*(1+Int_t(fTrackType/2)*9), range*(1+Int_t(fTrackType/2)*9), 200, -2, 2);
594 fRecDcaNegPtEtaNeg->GetXaxis()->SetTitle("dca (cm)");
595 fRecDcaNegPtEtaNeg->GetYaxis()->SetTitle("log_{10}(p_{T})");
596
1d55ca60 597
9f9631f3 598
1d55ca60 599 // YIELDs ---------------- for TPC sectors
600 for(Int_t sector=0; sector<18;sector++){
9f9631f3 601
602
1d55ca60 603 fRecPtTpcSector[sector] = new TH1F(Form("fRecPtTpcSector%02d",sector),
604 Form("p_{T} distribution: TPC sector %d",
605 sector),100, 0., pt);
606 fRecPtTpcSector[sector]->GetXaxis()->SetTitle("p_{T} (GeV/c)");
607
608 fRecEtaTpcSector[sector] = new TH1F(Form("fRecEtaTpcSector%02d",sector),
609 Form("#eta distribution: TPC sector %d",
610 sector),200, -2., 2.);
611 fRecEtaTpcSector[sector]->GetXaxis()->SetTitle("p_{T} (GeV/c)");
612
613
614 fSignedDcaTpcSector[sector] = new TH1F(Form("fSignedDcaTpcSector%02d",sector),
615 Form("dca distribution: TPC sector %d",
616 sector),200, -range*(1+Int_t(fTrackType/2)*9), range*(1+Int_t(fTrackType/2)*9) );
617 fSignedDcaTpcSector[sector]->GetXaxis()->SetTitle("dca");
618
619 fRecQPtTpcSector[sector] = new TH1F(Form("fRecQPtTpcSector%02d",sector),
620 Form("Q/ p_{T} distribution: TPC sector %d",
621 sector),100, -1., 1.);
622 fRecQPtTpcSector[sector]->GetXaxis()->SetTitle("Q/p_{T} (GeV/c)");
623
624 fRecEtaPtTpcSector[sector] = new TH1F(Form("fRecEtaPtTpcSector%02d",sector),
625 Form("#eta/ p_{T} distribution: TPC sector %d",
626 sector),100, -1., 1.);
627 fRecEtaPtTpcSector[sector]->GetXaxis()->SetTitle("#eta/p_{T} (GeV/c)");
628
629 }
630 // YIELDS ITS ladder
631 for(Int_t i=0;i<7;i++){
632 fRecPtPosLadder[i] = new TH1F(Form("fRecPtPosLadder%d", i),
633 " p_{T} distribution",
634 100, 0., pt);
635 fRecPtPosLadder[i]->GetXaxis()->SetTitle("p_{T} (GeV/c)");
636 fRecPtNegLadder[i] = new TH1F(Form("fRecPtNegLadder%d",i),
637 " p_{T} distribution ",
638 100, 0., pt);
639 fRecPtNegLadder[i]->GetXaxis()->SetTitle("p_{T} (GeV/c)");
640
641
642 fRecPhiPosLadder[i] = new TH1F(Form("fRecPhiPosLadder%d",i),
643 "#phi distribution: all pos eta",
644 361, 0., 360);
645 fRecPhiPosLadder[i]->GetXaxis()->SetTitle("#phi (deg)");
646
647 fRecPhiNegLadder[i] = new TH1F(Form("fRecPhiNegLadder%d", i),
648 "#phi distribution: all neg eta",
649 361, 0, 360);
650 fRecPhiNegLadder[i]->GetXaxis()->SetTitle("#phi (deg)");
651
652
653
654 fRecEtaPosLadder[i] = new TH1F(Form("fRecEtaPosLadder%d",i),
655 "#eta distribution",
656 200, -2., 2.);
657 fRecEtaPosLadder[i]->GetXaxis()->SetTitle("#eta)");
658
659 fRecEtaNegLadder[i] = new TH1F(Form("fRecEtaNegLadder%d", i),
660 "#eta distribution",
661 200, -2., 2.);
662 fRecEtaNegLadder[i]->GetXaxis()->SetTitle("#eta");
663 }
664
665 Double_t vzmax = 15;
666
667 fRecPtPosVz = new TH2F("fRecPtPosVz",
668 "p_{T} distribution vs Vz()",
669 100, -1., 2., 200,-vzmax,vzmax);
670 fRecPtPosVz->GetXaxis()->SetTitle("log_{10}(p_{T})");
9f9631f3 671
1d55ca60 672 fRecPtNegVz = new TH2F("fRecPtNegVz",
673 "p_{T} distribution vs Vz()",
674 100, -1., 2.,200,-vzmax,vzmax);
675 fRecPtNegVz->GetXaxis()->SetTitle("Log_{10}(p_{T})");
9f9631f3 676
677
1d55ca60 678 fRecEtaPosVz= new TH2F("fRecEtaPosVz",
679 "#eta distribution vs Vz()",
680 100, -2., 2., 200,-vzmax,vzmax);
681 fRecEtaPosVz->GetXaxis()->SetTitle("#eta");
682 fRecEtaNegVz = new TH2F("fRecEtaNegVz",
9f9631f3 683 "#eta distribution vs Vz()",
1d55ca60 684 100, -2., 2.,200,-vzmax,vzmax);
685 fRecEtaNegVz->GetXaxis()->SetTitle("#eta");
686
687 fRecPhiPosVz= new TH2F("fRecPhiPosVz",
688 "#eta distribution vs Vz()",
689 361, 0., 360., 200,-vzmax,vzmax);
690 fRecPhiPosVz->GetXaxis()->SetTitle("#phi (deg)");
691 fRecPhiNegVz = new TH2F("fRecPhiNegVz",
692 "dca vs Vz()",
693 361, 0., 360.,200,-vzmax,vzmax);
694 fRecPhiNegVz->GetXaxis()->SetTitle("#phi (deg)");
695
696 fSignedDcaPosVz= new TH2F("fSignedDcaPosVz",
697 "#eta distribution vs Vz()",
698 200, -range*(1+Int_t(fTrackType/2)*9), range*(1+Int_t(fTrackType/2)*9), 200,-vzmax,vzmax);
699 fSignedDcaPosVz->GetXaxis()->SetTitle("dca (cm)");
700 fSignedDcaNegVz = new TH2F("fSignedDcaNegVz",
701 "dca vs Vz()",
702 200, -range*(1+Int_t(fTrackType/2)*9), range*(1+Int_t(fTrackType/2)*9),200,-vzmax,vzmax);
703 fSignedDcaNegVz->GetXaxis()->SetTitle("dca (cm)");
704
705 fRecQPtPosEtaVz= new TH2F("fRecQPtPosEtaVz",
706 " Q/p_{T} distribution vs Vz()",
707 100, -1., 1., 200,-vzmax,vzmax);
708 fRecQPtPosEtaVz->GetXaxis()->SetTitle("Q/p_{T}");
709 fRecQPtNegEtaVz = new TH2F("fRecQPtNegEtaVz",
710 " Q/p_{T} distribution vs Vz()",
711 100, -1., 1.,200,-vzmax,vzmax);
712 fRecQPtNegEtaVz->GetXaxis()->SetTitle("Q/p_{T}");
9f9631f3 713
714
1d55ca60 715 fRecEtaPtPosVz= new TH2F("fRecEtaPtPosVz",
716 " #eta/p_{T} distribution vs Vz()",
717 100, -1., 1., 200,-vzmax,vzmax);
718 fRecEtaPtPosVz->GetXaxis()->SetTitle("#eta/p_{T");
719 fRecEtaPtNegVz = new TH2F("fRecEtaPtNegVz",
720 " #eta/p_{T} distribution vs Vz()",
721 100, -1., 1.,200,-vzmax,vzmax);
722 fRecEtaPtNegVz->GetXaxis()->SetTitle("#eta/p_{T}");
723
724 //new
725 fDeltaPhiAll = new TH1F("fDeltaPhiAll",
726 " #Delta #phi",200,-360,360);
727 fDeltaPhiAll->GetXaxis()->SetTitle("#Delta #phi");
728
729
730 fDeltaPhiLeading = new TH2F("fDeltaPhiLeading",
731 " #Delta #phi",361,-360,360, 361,0, 360);
732 fDeltaPhiLeading->GetXaxis()->SetTitle("#Delta #phi (deg.)");
733 fDeltaPhiLeading->GetYaxis()->SetTitle("#phi_{leading particle} (deg.)");
734
735 fDiffDcaD = new TH1F("fDiffDcaD",
736 "dca-d",
737 200, -1., 1.);
738
739
740 fRecPtPosEtaPos = new TH1F("fRecPtPosEtaPos",
741 "p_{T} distribution",100,0,pt);
742 fRecPtPosEtaPos->GetXaxis()->SetTitle("p_{T} (GeV/c)");
743
744 fRecPtNegEtaPos = new TH1F("fRecPtNegEtaPos",
745 "p_{T} distribution",100,0,pt);
746 fRecPtNegEtaPos->GetXaxis()->SetTitle("p_{T} (GeV/c)");
747
748 fRecPtPosEtaNeg = new TH1F("fRecPtPosEtaNeg",
749 "p_{T} distribution",100,0,pt);
750 fRecPtPosEtaNeg->GetXaxis()->SetTitle("p_{T} (GeV/c)");
751
752 fRecPtNegEtaNeg = new TH1F("fRecPtNegEtaNeg",
753 "p_{T} distribution",100,0,pt);
754 fRecPtNegEtaNeg->GetXaxis()->SetTitle("p_{T} (GeV/c)");
755
756
757
758 fRec1PtPosEtaPos = new TH1F("fRec1PtPosEtaPos",
759 "1/p_{T} distribution",100,0,0.5);
760 fRec1PtPosEtaPos->GetXaxis()->SetTitle("p_{T} (c/GeV)");
761
762 fRec1PtNegEtaPos = new TH1F("fRec1PtNegEtaPos",
763 "1/p_{T} distribution",100,0,0.5);
764 fRec1PtNegEtaPos->GetXaxis()->SetTitle("p_{T} (c/GeV)");
765
766 fRec1PtPosEtaNeg = new TH1F("fRec1PtPosEtaNeg",
767 "1/p_{T} distribution",100,0,0.5);
768 fRec1PtPosEtaNeg->GetXaxis()->SetTitle("p_{T} (c/GeV)");
769
770 fRec1PtNegEtaNeg = new TH1F("fRec1PtNegEtaNeg",
771 "1/p_{T} distribution",100,0,0.5);
772 fRec1PtNegEtaNeg->GetXaxis()->SetTitle("1/p_{T} (c/GeV)");
773
774
775
776 fRecPhiPosEtaPos = new TH1F("fRecPhiPosEtaPos",
777 "#phi",180,0,2*TMath::Pi());
778 fRecPhiPosEtaPos->GetXaxis()->SetTitle("#phi (rad.)");
779
780 fRecPhiNegEtaPos = new TH1F("fRecPhiNegEtaPos",
781 "#phi",180,0,2*TMath::Pi());
782 fRecPhiNegEtaPos->GetXaxis()->SetTitle("#phi (rad.)");
783
784 fRecPhiPosEtaNeg = new TH1F("fRecPhiPosEtaNeg",
785 "#phi",180,0,2*TMath::Pi());
786 fRecPhiPosEtaNeg->GetXaxis()->SetTitle("#phi (rad.)");
787
788 fRecPhiNegEtaNeg = new TH1F("fRecPhiNegEtaNeg",
789 "#phi",180,0,2*TMath::Pi());
790 fRecPhiNegEtaNeg->GetXaxis()->SetTitle("#phi (rad.)");
791
792
793
8ae6f5d0 794// fRecDcaPhiPtPosEtaPos = new TH3F("fRecDcaPhiPtPosEtaPos",
795// "#phi- p_{T} - DCA",
796// 180,0,2*TMath::Pi(),
797// 100,0,pt,
798// 100, -range*(1+Int_t(fTrackType/2)*9), range*(1+Int_t(fTrackType/2)*9));
799// fRecDcaPhiPtPosEtaPos->GetXaxis()->SetTitle("#phi (rad.)");
800// fRecDcaPhiPtPosEtaPos->GetYaxis()->SetTitle("p_{T} (GeV/c)");
801// fRecDcaPhiPtPosEtaPos->GetZaxis()->SetTitle("dca (cm)");
1d55ca60 802
8ae6f5d0 803// fRecDcaPhiPtPosEtaNeg = new TH3F("fRecDcaPhiPtPosEtaNeg",
804// "#phi- p_{T} - DCA",
805// 180,0,2*TMath::Pi(),
806// 100,0,pt,
807// 100, -range*(1+Int_t(fTrackType/2)*9), range*(1+Int_t(fTrackType/2)*9));
808// fRecDcaPhiPtPosEtaNeg->GetZaxis()->SetTitle("dca (cm)");
809// fRecDcaPhiPtPosEtaNeg->GetXaxis()->SetTitle("#phi (rad.)");
810// fRecDcaPhiPtPosEtaNeg->GetYaxis()->SetTitle("p_{T} (GeV/c)");
1d55ca60 811
8ae6f5d0 812// fRecDcaPhiPtNegEtaPos = new TH3F("fRecDcaPhiPtNegEtaPos",
813// "#phi- p_{T} - DCA",
814// 180,0,2*TMath::Pi(),
815// 100,0,pt,
816// 100, -range*(1+Int_t(fTrackType/2)*9), range*(1+Int_t(fTrackType/2)*9));
817// fRecDcaPhiPtNegEtaPos->GetZaxis()->SetTitle("dca (cm)");
818// fRecDcaPhiPtNegEtaPos->GetXaxis()->SetTitle("#phi (rad.)");
819// fRecDcaPhiPtNegEtaPos->GetYaxis()->SetTitle("p_{T} (GeV/c)");
1d55ca60 820
8ae6f5d0 821// fRecDcaPhiPtNegEtaNeg = new TH3F("fRecDcaPhiPtNegEtaNeg",
822// "#phi- p_{T} - DCA",
823// 180,0,2*TMath::Pi(),
824// 100,0,pt,
825// 100, -range*(1+Int_t(fTrackType/2)*9), range*(1+Int_t(fTrackType/2)*9));
826// fRecDcaPhiPtNegEtaNeg->GetZaxis()->SetTitle("dca (cm)");
827// fRecDcaPhiPtNegEtaNeg->GetYaxis()->SetTitle("#phi (rad.)");
828// fRecDcaPhiPtNegEtaNeg->GetXaxis()->SetTitle("p_{T} (GeV/c)");
1d55ca60 829
830
831
9f9631f3 832
833 fHists->SetOwner();
834
9f9631f3 835 fHists->Add(fHistRECpt);
836 fHists->Add(fEta);
1d55ca60 837 fHists->Add(fEtavPt);
98a5f582 838 fHists->Add(fCompareTPCparam);
e4ab7356 839 fHists->Add(fEtaPhi);
9f9631f3 840 fHists->Add(fThetaRec);
841 fHists->Add(fPhiRec);
842 fHists->Add(fNumber);
843 fHists->Add(fVx);
844 fHists->Add(fVy);
845 fHists->Add(fVz);
cbac53a7 846 fHists->Add(fVertexX);
847 fHists->Add(fVertexY);
848 fHists->Add(fVertexZ);
9f9631f3 849
850 fHists->Add(fEtaPt);
851 fHists->Add(fQPt);
852 fHists->Add(fDca);
853
854 fHists->Add(fDeltaPhiAll);
855 fHists->Add(fDeltaPhiLeading);
856 fHists->Add(fDiffDcaD);
857
858 fHists->Add(fqRec);
859 fHists->Add(fsigmaPt);
860
861 fHists->Add(fRecPtPos);
862 fHists->Add(fRecPtNeg);
863 fHists->Add(fRecPhiPos);
864 fHists->Add(fRecPhiNeg);
865 fHists->Add(fRecEtaPos);
866 fHists->Add(fRecEtaNeg);
867 fHists->Add(fRecEtaPtPos);
868 fHists->Add(fRecEtaPtNeg);
869 fHists->Add(fRecDcaPos);
870 fHists->Add(fRecDcaNeg);
871 fHists->Add(fRecDcaNegInv);
872 fHists->Add(fRecDPos);
873 fHists->Add(fRecDNeg);
874
875
876 fHists->Add(fRecQPtPosEta);
877 fHists->Add(fRecQPtNegEta);
878 fHists->Add(fRecPtPosEta);
879 fHists->Add(fRecPtNegEta);
880 fHists->Add(fRecPhiPosEta);
881 fHists->Add(fRecPhiNegEta);
882 fHists->Add(fRecDcaPosEta);
883 fHists->Add(fRecDcaNegEta);
884 fHists->Add(fRecDPosEta);
885 fHists->Add(fRecDNegEta);
886
887
888 for(Int_t i=0;i<18;i++){
889 fHists->Add(fRecPtTpcSector[i]);
890 fHists->Add(fRecEtaTpcSector[i]);
891 fHists->Add(fSignedDcaTpcSector[i]);
892 fHists->Add(fRecQPtTpcSector[i]);
893 fHists->Add(fRecEtaPtTpcSector[i]);
894 }
895
896 for(Int_t i=0;i<7;i++){
897 fHists->Add(fRecPtPosLadder[i]);
898 fHists->Add(fRecPtNegLadder[i]);
899 fHists->Add(fRecPhiPosLadder[i]);
900 fHists->Add(fRecPhiNegLadder[i]);
901 fHists->Add(fRecEtaPosLadder[i]);
902 fHists->Add(fRecEtaNegLadder[i]);
903 }
904
905 fHists->Add(fRecPtPosVz);
906 fHists->Add(fRecPtNegVz);
907 fHists->Add(fRecEtaPosVz);
908 fHists->Add(fRecEtaNegVz);
909 fHists->Add(fRecPhiPosVz);
910 fHists->Add(fRecPhiNegVz);
911 fHists->Add(fSignedDcaPosVz);
912 fHists->Add(fSignedDcaNegVz);
913 fHists->Add(fRecQPtPosEtaVz);
914 fHists->Add(fRecQPtNegEtaVz);
915 fHists->Add(fRecEtaPtPosVz);
916 fHists->Add(fRecEtaPtNegVz);
917
918
919 for(Int_t i=0;i<7;i++){
920 fHists->Add(fSignDcaPos[i]);
921 fHists->Add(fSignDcaNeg[i]);
922 fHists->Add(fSignDcaNegInv[i]);
923
924 fHists->Add(fPtSigmaPos[i]);
925 fHists->Add(fPtSigmaNeg[i]);
926 fHists->Add(fqPtRec[i]);
927
928 fHists->Add(fDcaSigmaPos[i]);
929 fHists->Add(fDcaSigmaNeg[i]);
930
931
932 }
933
1d55ca60 934 fHists->Add(fRecDcaPosPhi);
935 fHists->Add(fRecDcaNegPhi);
936 fHists->Add(fRecPtPosPhi);
937 fHists->Add(fRecPtNegPhi);
938 fHists->Add(fRecEtaPosPhi);
939 fHists->Add(fRecEtaNegPhi);
940 fHists->Add(fRecQPtPhi);
941 fHists->Add(fRecEtaPtPosPhi);
942 fHists->Add(fRecEtaPtNegPhi);
943
944 fHists->Add(fRecPtPosEtaPos);
945 fHists->Add(fRecPtNegEtaPos);
946 fHists->Add(fRecPtPosEtaNeg);
947 fHists->Add(fRecPtNegEtaNeg);
948
949 fHists->Add(fRec1PtPosEtaPos);
950 fHists->Add(fRec1PtNegEtaPos);
951 fHists->Add(fRec1PtPosEtaNeg);
952 fHists->Add(fRec1PtNegEtaNeg);
953
954
955 fHists->Add(fRecPhiPosEtaPos);
956 fHists->Add(fRecPhiNegEtaPos);
957 fHists->Add(fRecPhiPosEtaNeg);
958 fHists->Add(fRecPhiNegEtaNeg);
959
960 fHists->Add(fRecDcaPosPhiEtaPos);
961 fHists->Add(fRecDcaNegPhiEtaPos);
962 fHists->Add(fRecPtPosPhiEtaPos);
963 fHists->Add(fRecPtNegPhiEtaPos);
964 fHists->Add(fRecDcaPosPhiEtaNeg);
965 fHists->Add(fRecDcaNegPhiEtaNeg);
966 fHists->Add(fRecPtPosPhiEtaNeg);
967 fHists->Add(fRecPtNegPhiEtaNeg);
968
98a5f582 969 fHists->Add(fRecDcaPosPtEtaPos);
970 fHists->Add(fRecDcaNegPtEtaPos);
971 fHists->Add(fRecDcaPosPtEtaNeg);
972 fHists->Add(fRecDcaNegPtEtaNeg);
973
8ae6f5d0 974 // fHists->Add(fRecDcaPhiPtPosEtaPos);
975 // fHists->Add(fRecDcaPhiPtPosEtaNeg);
976 // fHists->Add(fRecDcaPhiPtNegEtaPos);
977 // fHists->Add(fRecDcaPhiPtNegEtaNeg);
9f9631f3 978
979
1d55ca60 980
981
982
983// for (Int_t i=0; i<fHists->GetEntries(); ++i) {
984// TH1 *h1 = dynamic_cast<TH1*>(fHists->At(i));
985// if (h1){
986// // Printf("%s ",h1->GetName());
987// h1->Sumw2();
988// }
989// }
990 // BKC
991
992 TH1::AddDirectory(oldStatus);
9f9631f3 993}
994
995//__________________________________________________________
996
997void AliAnalysisTaskQASym::UserExec(Option_t *)
998{
999 AliVEvent *event = InputEvent();
1000 if (!event) {
1d55ca60 1001 Printf("ERROR: Could not retrieve event");
1002 return;
9f9631f3 1003 }
1004
1005
1006 if(Entry()==0){
1007 AliESDEvent* esd = dynamic_cast<AliESDEvent*>(event);
1008 if(esd){
1009 Printf("We are reading from ESD");
1010 }
1011
1012 }
1013
1d55ca60 1014
1015
1016 if(fDebug>1)Printf("There are %d tracks in this event", event->GetNumberOfTracks());
1017
9f9631f3 1018
50dbca26 1019 Int_t leadingTrack = 0;
1020 Float_t leadingEnergy = -20.;
1021 Float_t leadingPhi = 0;//TMath::Pi();
9f9631f3 1022
1023
1024 if(event->GetNumberOfTracks()!=0) fNumber->Fill(event->GetNumberOfTracks());
1025
e4ab7356 1026 const AliVVertex* vertex = event->GetPrimaryVertex();
8ae6f5d0 1027 if(vertex->GetNContributors()==0) return;
cbac53a7 1028 Float_t vx = vertex->GetX();
1029 Float_t vy = vertex->GetY();
e4ab7356 1030 Float_t vz = vertex->GetZ();
cbac53a7 1031
1032 fVertexX->Fill(vx);
1033 fVertexY->Fill(vy);
1034 fVertexZ->Fill(vz);
1035
e4ab7356 1036 if (TMath::Abs(vz) > 10.) return;
1037
1fcb92ef 1038 AliESDtrack *tpcP = 0x0;
1039
9f9631f3 1040 for (Int_t iTrack = 0; iTrack < event->GetNumberOfTracks(); iTrack++) {
1041
1fcb92ef 1042 //prevent mem leak for TPConly track
1043 if(fTrackType==2&&tpcP){
1044 delete tpcP;
1045 tpcP = 0;
1046 }
9f9631f3 1047
1048 AliVParticle *track = event->GetTrack(iTrack);
1049 AliESDtrack *esdtrack = dynamic_cast<AliESDtrack*>(track);
1d55ca60 1050 esdtrack->PropagateToDCA(event->GetPrimaryVertex(),
1051 event->GetMagneticField(), 10000.);
1052
9f9631f3 1053 if (!track) {
1054 Printf("ERROR: Could not receive track %d", iTrack);
1055 continue;
1056 }
1d55ca60 1057 //__________
1058 // run Task for global tracks or ITS tracks or TPC tracks
1d55ca60 1059 const AliExternalTrackParam *tpcPin = 0x0;
1060 Double_t phiIn=0;
cbac53a7 1061
1d55ca60 1062 if(fTrackType==0){
1063 //Fill all histograms with global tracks
1064 tpcP = esdtrack;
1065 phiIn = tpcP->Phi();
1066 if (!tpcP) continue;
1067 if (!fCuts->AcceptTrack(tpcP)) continue;
1068 }
1069 else if(fTrackType==1){
1070 //Fill all histograms with ITS tracks
1071 tpcP = esdtrack;
1072 phiIn = tpcP->Phi();
1073 if (!tpcP) continue;
1074 if (!fCuts->AcceptTrack(tpcP)) continue;
1075 }
1076 else if(fTrackType==2){
1077 //Fill all histograms with TPC track information
1078 tpcPin = esdtrack->GetInnerParam();
1079 if (!tpcPin) continue;
1080 phiIn=tpcPin->Phi();
1081
1082 tpcP = AliESDtrackCuts::GetTPCOnlyTrack(dynamic_cast<AliESDEvent*>(event),esdtrack->GetID());
1083 if (!tpcP) continue;
1084 if (!fCuts->AcceptTrack(tpcP)) continue;
1085 if(tpcP->GetNcls(1)>160)continue;//jacek's track cut
1086 if(tpcP->GetConstrainedChi2TPC()<0)continue; // jacek's track cut
1087 }
1088 else{
1089 Printf("ERROR: wrong track type \n");
1090 continue;
1091 }
1092 //___________
1093 //
1094
1095
50dbca26 1096 if(tpcP->E()>leadingEnergy){
1097 leadingTrack=iTrack;
1098 leadingEnergy=tpcP->E();
1d55ca60 1099 leadingPhi=phiIn;
9f9631f3 1100 }
1101
98a5f582 1102
9f9631f3 1103 fqRec->Fill(tpcP->Charge());
1104
1105
1106 Double_t sigmapt = tpcP->GetSigma1Pt2();
1107 sigmapt= sqrt(sigmapt);
1108 sigmapt= sigmapt *(tpcP->Pt()*tpcP->Pt());
1109
1d55ca60 1110 if(sigmapt == 0.)continue;
9f9631f3 1111 fsigmaPt->Fill(TMath::Log10(sigmapt));
1d55ca60 1112
9f9631f3 1113
9f9631f3 1114 // hits in ITS layer
1115 Int_t cas=-1;
1d55ca60 1116 if(tpcP->HasPointOnITSLayer(0))
9f9631f3 1117 cas=0;
1d55ca60 1118 else if(!tpcP->HasPointOnITSLayer(0)
1119 && tpcP->HasPointOnITSLayer(1))
9f9631f3 1120 cas=1;
1d55ca60 1121 else if(!tpcP->HasPointOnITSLayer(0)
1122 && !tpcP->HasPointOnITSLayer(1)
1123 && tpcP->HasPointOnITSLayer(2))
9f9631f3 1124 cas=2;
1d55ca60 1125 else if(!tpcP->HasPointOnITSLayer(0)
1126 && !tpcP->HasPointOnITSLayer(1)
1127 && !tpcP->HasPointOnITSLayer(2)
1128 && tpcP->HasPointOnITSLayer(3))
9f9631f3 1129 cas=3;
1d55ca60 1130 else if(!tpcP->HasPointOnITSLayer(0)
1131 && !tpcP->HasPointOnITSLayer(1)
1132 && !tpcP->HasPointOnITSLayer(2)
1133 && !tpcP->HasPointOnITSLayer(3)
1134 && tpcP->HasPointOnITSLayer(4))
9f9631f3 1135 cas=4;
1d55ca60 1136 else if( !tpcP->HasPointOnITSLayer(0)
1137 && !tpcP->HasPointOnITSLayer(1)
1138 && !tpcP->HasPointOnITSLayer(2)
1139 && !tpcP->HasPointOnITSLayer(3)
1140 && !tpcP->HasPointOnITSLayer(4)
1141 && tpcP->HasPointOnITSLayer(5))
9f9631f3 1142 cas=5;
1143 else
1144 cas=6;
1145
1146
1147
1148 //-------------------
1149
1d55ca60 1150 xvertexcor = tpcP->Xv() - vertex->GetX(); // coordinate corrected for vertex position
1151 yvertexcor = tpcP->Yv() - vertex->GetY(); // "
1152 sdca = (tpcP->Py()*xvertexcor - tpcP->Px()*yvertexcor)/tpcP->Pt();
9f9631f3 1153
1154
1155 fqPtRec[cas]->Fill(tpcP->Charge()/tpcP->Pt());
1156
1d55ca60 1157
1158
9f9631f3 1159 fHistRECpt->Fill(tpcP->Pt());
1160 fEta->Fill(tpcP->Eta());
1d55ca60 1161 fEtavPt->Fill(tpcP->Eta(), tpcP->Pt());
1162 fEtaPhi->Fill(tpcP->Eta(), phiIn);
9f9631f3 1163 fThetaRec->Fill(tpcP->Theta());
1d55ca60 1164 fPhiRec->Fill(phiIn);
9f9631f3 1165 fVx->Fill(tpcP->Xv());
1166 fVy->Fill(tpcP->Yv());
1167 fVz->Fill(tpcP->Zv());
1168
1169
1170 fEtaPt->Fill(tpcP->Eta()/tpcP->Pt());
1171 fQPt->Fill(tpcP->Charge()/tpcP->Pt());
1d55ca60 1172 fDca->Fill(sdca);
1173 fRecQPtPhi->Fill(tpcP->Charge()/tpcP->Pt(), phiIn);
9f9631f3 1174
1175
1d55ca60 1176 tpcP->GetImpactParameters(xy,z);
1177 fDiffDcaD->Fill(sdca+xy);
9f9631f3 1178
98a5f582 1179 if(fTrackType==2) fCompareTPCparam->Fill(z,tpcPin->GetTgl());
9f9631f3 1180
1181 //for positive particles
e4ab7356 1182
9f9631f3 1183 if(tpcP->Charge()>0){
1184 fRecPtPos->Fill(tpcP->Pt());
1185 fRecPtPosLadder[cas]->Fill(tpcP->Pt());
1186 fRecPtPosVz->Fill(TMath::Log10(tpcP->Pt()),tpcP->Zv());
1d55ca60 1187 fRecPhiPos->Fill(TMath::RadToDeg()*phiIn);
9f9631f3 1188
1189
1d55ca60 1190 fRecPhiPosLadder[cas]->Fill(TMath::RadToDeg()*phiIn);
1191 fRecPhiPosVz->Fill(TMath::RadToDeg()*phiIn,tpcP->Zv());
1192 fSignedDcaPosVz->Fill(sdca,tpcP->Zv());
9f9631f3 1193
1194 fRecEtaPos->Fill(tpcP->Eta());
1195 fRecEtaPosLadder[cas]->Fill(tpcP->Eta());
1196 fRecEtaPtPos->Fill(tpcP->Eta()/tpcP->Pt());
1197 fRecEtaPosVz->Fill(tpcP->Eta(),tpcP->Zv());
1198 fRecEtaPtPosVz->Fill(tpcP->Eta()/tpcP->Pt(),tpcP->Zv());
1199
1d55ca60 1200 fRecDcaPos->Fill(sdca);
1201 fRecDcaPosPhi->Fill(sdca, phiIn);
1202 fRecPtPosPhi->Fill(TMath::Log10(tpcP->Pt()), phiIn);
1203 fRecEtaPtPosPhi->Fill(tpcP->Eta()/tpcP->Pt(), phiIn);
1204 fRecEtaPosPhi->Fill(tpcP->Eta(), phiIn);
9f9631f3 1205 fRecDPos->Fill(xy);
1d55ca60 1206 fSignDcaPos[cas]->Fill(sdca);
9f9631f3 1207
1208
1d55ca60 1209 fDcaSigmaPos[cas]->Fill(sdca, TMath::Log10(sigmapt));
9f9631f3 1210
1211 fPtSigmaPos[cas]->Fill(TMath::Log10(sigmapt));
1d55ca60 1212 //pos eta
1213 if(tpcP->Eta()>0){
1214 fRecPtPosEtaPos->Fill(tpcP->Pt());
1215 fRec1PtPosEtaPos->Fill(1/tpcP->Pt());
1216 fRecPhiPosEtaPos->Fill(phiIn);
1217 fRecDcaPosPhiEtaPos->Fill(sdca, phiIn);
98a5f582 1218 fRecDcaPosPtEtaPos->Fill(sdca, TMath::Log10(tpcP->Pt()));
1d55ca60 1219 fRecPtPosPhiEtaPos->Fill(TMath::Log10(tpcP->Pt()), phiIn);
8ae6f5d0 1220 //fRecDcaPhiPtPosEtaPos->Fill(phiIn, tpcP->Pt(), sdca);
1d55ca60 1221 }
1222 //neg eta
1223 else{
1224 fRecPtPosEtaNeg->Fill(tpcP->Pt());
1225 fRec1PtPosEtaNeg->Fill(1/tpcP->Pt());
1226 fRecPhiPosEtaNeg->Fill(phiIn);
1227 fRecDcaPosPhiEtaNeg->Fill(sdca, phiIn);
98a5f582 1228 fRecDcaPosPtEtaNeg->Fill(sdca, TMath::Log10(tpcP->Pt()));
1d55ca60 1229 fRecPtPosPhiEtaNeg->Fill(TMath::Log10(tpcP->Pt()), phiIn);
8ae6f5d0 1230 //fRecDcaPhiPtPosEtaNeg->Fill(phiIn, tpcP->Pt(), sdca);
1d55ca60 1231 }
1232
9f9631f3 1233 }
1234 //and negative particles
1235 else {
1236 fRecPtNeg->Fill(tpcP->Pt());
1237 fRecPtNegLadder[cas]->Fill(tpcP->Pt());
1238 fRecPtNegVz->Fill(TMath::Log10(tpcP->Pt()),tpcP->Zv());
1239
1d55ca60 1240 fRecPhiNeg->Fill(TMath::RadToDeg()*phiIn);
1241 fRecPhiNegLadder[cas]->Fill(TMath::RadToDeg()*phiIn);
1242 fRecPhiNegVz->Fill(TMath::RadToDeg()*phiIn,tpcP->Zv());
1243 fSignedDcaNegVz->Fill(sdca,tpcP->Zv());
9f9631f3 1244 fRecEtaPtNegVz->Fill(tpcP->Eta()/tpcP->Pt(),tpcP->Zv());
1245
1246 fRecEtaNeg->Fill(tpcP->Eta());
1247 fRecEtaNegLadder[cas]->Fill(tpcP->Eta());
1248 fRecEtaPtNeg->Fill(tpcP->Eta()/tpcP->Pt());
1249 fRecEtaNegVz->Fill(tpcP->Eta(),tpcP->Zv());
1250
1d55ca60 1251 fRecDcaNeg->Fill(sdca);
1252 fRecDcaNegInv->Fill(-sdca);
1253 fRecDcaNegPhi->Fill(sdca, phiIn);
1254 fRecPtNegPhi->Fill(TMath::Log10(tpcP->Pt()), phiIn);
1255 fRecEtaNegPhi->Fill(tpcP->Eta(), phiIn);
1256 fRecEtaPtNegPhi->Fill(tpcP->Eta()/tpcP->Pt(), phiIn);
9f9631f3 1257 fRecDNeg->Fill(xy);
1d55ca60 1258 fSignDcaNeg[cas]->Fill(sdca);
1259 fSignDcaNegInv[cas]->Fill(-sdca);
9f9631f3 1260
1261
1d55ca60 1262 fDcaSigmaNeg[cas]->Fill(sdca,TMath::Log10(sigmapt));
9f9631f3 1263
1264 fPtSigmaNeg[cas]->Fill(TMath::Log10(sigmapt));
1d55ca60 1265
1266 //pos eta
1267 if(tpcP->Eta()>0){
1268 fRecPtNegEtaPos->Fill(tpcP->Pt());
1269 fRec1PtNegEtaPos->Fill(1/tpcP->Pt());
1270 fRecPhiNegEtaPos->Fill(phiIn);
1271 fRecDcaNegPhiEtaPos->Fill(sdca, phiIn);
98a5f582 1272 fRecDcaNegPtEtaPos->Fill(sdca, TMath::Log10(tpcP->Pt()));
1d55ca60 1273 fRecPtNegPhiEtaPos->Fill(TMath::Log10(tpcP->Pt()), phiIn);
8ae6f5d0 1274 //fRecDcaPhiPtNegEtaPos->Fill(phiIn, tpcP->Pt(), sdca);
1d55ca60 1275 }
1276 //neg eta
1277 else{
1278 fRecPtNegEtaNeg->Fill(tpcP->Pt());
1279 fRec1PtNegEtaNeg->Fill(1/tpcP->Pt());
1280 fRecPhiNegEtaNeg->Fill(phiIn);
1281 fRecDcaNegPhiEtaNeg->Fill(sdca, phiIn);
98a5f582 1282 fRecDcaNegPtEtaNeg->Fill(sdca, TMath::Log10(tpcP->Pt()));
1d55ca60 1283 fRecPtNegPhiEtaNeg->Fill(TMath::Log10(tpcP->Pt()), phiIn);
8ae6f5d0 1284 //fRecDcaPhiPtNegEtaNeg->Fill(phiIn, tpcP->Pt(), sdca);
1d55ca60 1285 }
1286
9f9631f3 1287 }
1288
1289
1290
1291 //all particles with positive eta
1292 if(tpcP->Eta()>0){
1293 fRecQPtPosEta->Fill(tpcP->Charge()/tpcP->Pt());
1294 fRecPtPosEta->Fill(tpcP->Pt());
1d55ca60 1295 fRecPhiPosEta->Fill(TMath::RadToDeg()*phiIn);
9f9631f3 1296 fRecQPtPosEtaVz->Fill(tpcP->Charge()/tpcP->Pt(),tpcP->Zv());
1d55ca60 1297 fRecDcaPosEta->Fill(sdca);
9f9631f3 1298 fRecDPosEta->Fill(xy);
1299 }
1300 //all particles with negative eta (and eta==0)
1301 else{
1302 fRecQPtNegEta->Fill(tpcP->Charge()/tpcP->Pt());
1303 fRecPtNegEta->Fill(tpcP->Pt());
1d55ca60 1304 fRecPhiNegEta->Fill(TMath::RadToDeg()*phiIn);
9f9631f3 1305 fRecQPtNegEtaVz->Fill(tpcP->Charge()/tpcP->Pt(),tpcP->Zv());
1d55ca60 1306 fRecDcaNegEta->Fill(sdca);
9f9631f3 1307 fRecDNegEta->Fill(xy);
1308
1309 }
9f9631f3 1310
1d55ca60 1311
1312 fRecPtTpcSector[Int_t(phiIn*
9f9631f3 1313 TMath::RadToDeg()/20)]->Fill(tpcP->Pt());
1d55ca60 1314 fRecEtaTpcSector[Int_t(phiIn*
1315 TMath::RadToDeg()/20)]->Fill(tpcP->Eta());
1316 fSignedDcaTpcSector[Int_t(phiIn*
1317 TMath::RadToDeg()/20)]->Fill(sdca);
1318 fRecQPtTpcSector[Int_t(phiIn*
1319 TMath::RadToDeg()/20)]->Fill(tpcP->Charge()/tpcP->Pt());
1320 fRecEtaPtTpcSector[Int_t(phiIn*
1321 TMath::RadToDeg()/20)]->Fill(tpcP->Eta()/tpcP->Pt());
9f9631f3 1322
1323
1324
1d55ca60 1325// // another track loop
1326// for (Int_t iTrack2 = 0; iTrack2 < event->GetNumberOfTracks(); iTrack2++) {
9f9631f3 1327
1d55ca60 1328// if(LeadingTrack==iTrack2) continue;
1329
1330// AliVParticle *track2 = event->GetTrack(iTrack2);
1331// AliESDtrack* esdtrack2 = dynamic_cast<AliESDtrack*>(track2);
1332// if (!track2) {
1333// Printf("ERROR: Could not receive track %d", iTrack);
1334// continue;
1335// }
1336// if (!fCuts->AcceptTrack(esdtrack2)) continue;
1337// //propagate to dca
1338// esdtrack2->PropagateToDCA(event->GetPrimaryVertex(),
1339// event->GetMagneticField(), 10000.);
9f9631f3 1340
1d55ca60 1341// fDeltaPhiLeading->Fill((LeadingPhi-esdtrack2->Phi())*TMath::RadToDeg(),
1342// LeadingPhi*TMath::RadToDeg() );
9f9631f3 1343
1344
1345
1d55ca60 1346// }//second track loop
98a5f582 1347
1fcb92ef 1348 // if(fTrackType==2) delete tpcP; // delete in case of TPCOnlyTrack
98a5f582 1349
9f9631f3 1350 }//first track loop
1351
1fcb92ef 1352 //prevent mem leak for TPConly track
1353 if(fTrackType==2&&tpcP){
1354 delete tpcP;
1355 tpcP = 0;
1356 }
9f9631f3 1357
9f9631f3 1358
1359 // Post output data.
1360 // PostData(1, fHistPt);
1361 PostData(1, fHists);
1362}
1363
1364
1365
1366
1367
1368//________________________________________________________________________
1369void AliAnalysisTaskQASym::Terminate(Option_t *)
1370{
1d55ca60 1371
9f9631f3 1372
1373}
1374
1375
1376
1377
1378