1 void SPD_ntuple_Dubna()
4 TFile *f = new TFile("SPD_his_dubna.root");
5 //TH1F *Pix = (TH1F*)f->Get("Pix");
7 //gStyle->SetOptStat(1111111);
8 //gStyle->SetOptLogy();
9 //TCanvas *c1 = new TCanvas("c1","SPD clusters",400,10,600,700);
10 TCanvas *c2 = new TCanvas("c2","SPD clusters",400,10,600,700);
14 ///////////////////////// Ntuple analysis ///////////////////////////////
16 // ntuple is created inside the hit loop for the hits in the cluster region;
18 // ntuple1 is created after a finish of the hit loop if one or more hits
19 // are in the cluster region;
21 // ntuple2 is created befor the hit loop for all clusters;
23 // -----------------------------------------------------------------------
24 // lay - number of ITS layer;
25 // x - coordinates in r*phi(x) direction (mm);
26 // z - coordinates in z direction (mm);
27 // nx - cluster size in the r*phi(x) direction;
28 // nz - cluster size in the z direction;
29 // hitprim - primary particle(hit) flag ( = 1 for primery particle);
30 // dx - difference of hit(mediate) and reconstructed (from cluster)
31 // coordinates in r*phi(x) direction;
32 // dz - difference of hit(mediate) and reconstructed (from cluster)
33 // coordinates in z direction;
34 // noverlaps - number of particles overlapping in one cluster found in this
35 // macros in the cluster region;
36 // ntrover - number of particles overlapping in one cluster found in the
37 // AliITSClusterFinderSPD class;
38 // noverprim - number of primary particles overlapping in one cluster;
39 // qcl - cluster charge in electron number
40 // -------------------------------------------------------------------------
44 gPad->SetFillColor(33);
45 ntuple->SetFillColor(42);
46 ntuple->Draw("dx","lay == 1");
48 gPad->SetFillColor(33);
49 ntuple->SetFillColor(46);
50 ntuple->Draw("dz","lay == 1");
52 gPad->SetFillColor(33);
53 ntuple->SetFillColor(42);
54 ntuple->Draw("dx","lay == 2");
56 gPad->SetFillColor(33);
57 ntuple->SetFillColor(46);
58 ntuple->Draw("dz","lay == 2");
63 gPad->SetFillColor(33);
64 ntuple->SetFillColor(42);
65 ntuple->Draw("dx","lay == 1 && hitprim == 1");
67 gPad->SetFillColor(33);
68 ntuple->SetFillColor(46);
69 ntuple->Draw("dz","lay == 1 && hitprim == 1");
71 gPad->SetFillColor(33);
72 ntuple->SetFillColor(42);
73 ntuple->Draw("dx","lay == 2 && hitprim == 1");
75 gPad->SetFillColor(33);
76 ntuple->SetFillColor(46);
77 ntuple->Draw("dz","lay == 2 && hitprim == 1");
82 gPad->SetFillColor(33);
83 ntuple->SetFillColor(42);
84 ntuple->Draw("dx","lay == 1 && hitprim == 1 && ntrover>=1&&nx==1");
86 gPad->SetFillColor(33);
87 ntuple->SetFillColor(46);
88 ntuple->Draw("dz","lay == 1 && hitprim == 1 && ntrover>=1&&nz==1");
90 gPad->SetFillColor(33);
91 ntuple->SetFillColor(42);
92 ntuple->Draw("dx","lay == 2 && hitprim == 1 && ntrover>=1&&nx==1");
94 gPad->SetFillColor(33);
95 ntuple->SetFillColor(46);
96 ntuple->Draw("dz","lay == 2 && hitprim == 1 && ntrover>=1&&nz==1");
102 gPad->SetFillColor(33);
103 ntuple->SetFillColor(42);
104 ntuple->Draw("dx","lay==1&&hitprim==1&&abs(dx)<150&&nx==1");
106 gPad->SetFillColor(33);
107 ntuple->SetFillColor(46);
108 ntuple->Draw("dz","lay==1&&hitprim==1&&abs(dz)<500&&nz==1");
110 gPad->SetFillColor(33);
111 ntuple->SetFillColor(42);
112 ntuple->Draw("dx","lay==2&&hitprim==1&&abs(dx)<150&&nx==1");
114 gPad->SetFillColor(33);
115 ntuple->SetFillColor(46);
116 ntuple->Draw("dz","lay==2&&hitprim==1&&abs(dz)<500&&nz==1");
121 gPad->SetFillColor(33);
122 ntuple->SetFillColor(42);
123 ntuple->Draw("dx","lay==1&&hitprim==1&&abs(dx)<150&&ntrover==1&&nx>0");
125 gPad->SetFillColor(33);
126 ntuple->SetFillColor(46);
127 ntuple->Draw("dz","lay==1&&hitprim==1&&abs(dz)<500&&ntrover==1&&nz>0");
129 gPad->SetFillColor(33);
130 ntuple->SetFillColor(42);
131 ntuple->Draw("dx","lay==2&&hitprim==1&&abs(dx)<150&&ntrover==1&&nx>0");
133 gPad->SetFillColor(33);
134 ntuple->SetFillColor(46);
135 ntuple->Draw("dz","lay==2&&hitprim==1&&abs(dz)<500&&ntrover==1&&nz>0");
140 gPad->SetFillColor(33);
141 ntuple->SetFillColor(42);
142 ntuple->Draw("dx","lay==1&&hitprim==1&&nx==1");
144 gPad->SetFillColor(33);
145 ntuple->SetFillColor(46);
146 ntuple->Draw("dz","lay==1&&hitprim==1&&nz==1");
148 gPad->SetFillColor(33);
149 ntuple->SetFillColor(42);
150 ntuple->Draw("dx","lay==2&&hitprim==1&&nx==1");
152 gPad->SetFillColor(33);
153 ntuple->SetFillColor(46);
154 ntuple->Draw("dz","lay==2&&hitprim==1&&nz==1");
159 gPad->SetFillColor(33);
160 ntuple1->SetFillColor(42);
161 ntuple1->Draw("noverlaps","lay == 1&&noverprim>0&&noverlaps>1");
163 gPad->SetFillColor(33);
164 ntuple1->SetFillColor(46);
165 ntuple1->Draw("noverlaps","lay == 2&&noverprim>0&&noverlaps>1");
167 gPad->SetFillColor(33);
168 ntuple1->SetFillColor(42);
169 ntuple1->Draw("ntrover","lay == 1&&noverprim>0&&ntrover>1");
171 gPad->SetFillColor(33);
172 ntuple1->SetFillColor(46);
173 ntuple1->Draw("ntrover","lay == 2&&noverprim>0&&ntrover>1");
178 gPad->SetFillColor(33);
179 ntuple1->SetFillColor(42);
180 ntuple1->Draw("noverprim","lay == 1");
182 gPad->SetFillColor(33);
183 ntuple1->SetFillColor(46);
184 ntuple1->Draw("noverprim","lay == 2");
186 gPad->SetFillColor(33);
187 ntuple1->SetFillColor(42);
188 ntuple1->Draw("noverprim","lay == 1 && noverprim < 1");
190 gPad->SetFillColor(33);
191 ntuple1->SetFillColor(46);
192 ntuple1->Draw("noverprim","lay == 2 && noverprim < 1");
197 gPad->SetFillColor(33);
198 ntuple1->SetFillColor(42);
199 ntuple1->Draw("nx:nz","lay == 1");
201 gPad->SetFillColor(33);
202 ntuple1->SetFillColor(46);
203 ntuple1->Draw("nx:nz","lay == 2");
208 gPad->SetFillColor(33);
209 ntuple1->SetFillColor(42);
210 ntuple1->Draw("qcl","lay == 1");
212 gPad->SetFillColor(33);
213 ntuple1->SetFillColor(46);
214 ntuple1->Draw("qcl","lay == 2");
219 gPad->SetFillColor(33);
220 ntuple1->SetFillColor(42);
221 ntuple1->Draw("nx","lay == 1");
223 gPad->SetFillColor(33);
224 ntuple1->SetFillColor(46);
225 ntuple1->Draw("nx","lay == 2");
227 gPad->SetFillColor(33);
228 ntuple1->SetFillColor(42);
229 ntuple1->Draw("nz","lay == 1");
231 gPad->SetFillColor(33);
232 ntuple1->SetFillColor(46);
233 ntuple1->Draw("nz","lay == 2");
238 gPad->SetFillColor(33);
239 ntuple1->SetFillColor(42);
240 ntuple1->Draw("nx","lay == 1&&noverprim>0");
242 gPad->SetFillColor(33);
243 ntuple1->SetFillColor(46);
244 ntuple1->Draw("nx","lay == 2&&noverprim>0");
246 gPad->SetFillColor(33);
247 ntuple1->SetFillColor(42);
248 ntuple1->Draw("nz","lay == 1 && noverprim>0");
250 gPad->SetFillColor(33);
251 ntuple1->SetFillColor(46);
252 ntuple1->Draw("nz","lay == 2 && noverprim>0");
257 gPad->SetFillColor(33);
258 ntuple1->SetFillColor(42);
259 ntuple1->Draw("x","lay == 1&&noverprim>=0");
261 gPad->SetFillColor(33);
262 ntuple1->SetFillColor(46);
263 ntuple1->Draw("x","lay == 2&&noverprim>=0");
265 gPad->SetFillColor(33);
266 ntuple1->SetFillColor(42);
267 ntuple1->Draw("z","lay == 1 && noverprim>=0");
269 gPad->SetFillColor(33);
270 ntuple1->SetFillColor(46);
271 ntuple1->Draw("z","lay == 2 && noverprim>=0");
276 gPad->SetFillColor(33);
277 ntuple1->SetFillColor(42);
278 ntuple1->Draw("nx","lay == 1&&noverprim>0 && ntrover == 1");
280 gPad->SetFillColor(33);
281 ntuple1->SetFillColor(46);
282 ntuple1->Draw("nx","lay == 2&&noverprim>0 && ntrover == 1");
284 gPad->SetFillColor(33);
285 ntuple1->SetFillColor(42);
286 ntuple1->Draw("nz","lay == 1 && noverprim>0 && ntrover == 1");
288 gPad->SetFillColor(33);
289 ntuple1->SetFillColor(46);
290 ntuple1->Draw("nz","lay == 2 && noverprim>0 && ntrover ==1");
295 gPad->SetFillColor(33);
296 ntuple1->SetFillColor(42);
297 ntuple1->Draw("dx","lay == 1&&noverprim>0");
299 gPad->SetFillColor(33);
300 ntuple1->SetFillColor(46);
301 ntuple1->Draw("dx","lay == 2&&noverprim>0");
303 gPad->SetFillColor(33);
304 ntuple1->SetFillColor(42);
305 ntuple1->Draw("dz","lay == 1 && noverprim>0");
307 gPad->SetFillColor(33);
308 ntuple1->SetFillColor(46);
309 ntuple1->Draw("dz","lay == 2 && noverprim>0");
314 gPad->SetFillColor(33);
315 ntuple1->SetFillColor(42);
316 ntuple1->Draw("dx","lay == 1");
318 gPad->SetFillColor(33);
319 ntuple1->SetFillColor(46);
320 ntuple1->Draw("dx","lay == 2");
322 gPad->SetFillColor(33);
323 ntuple1->SetFillColor(42);
324 ntuple1->Draw("dz","lay == 1");
326 gPad->SetFillColor(33);
327 ntuple1->SetFillColor(46);
328 ntuple1->Draw("dz","lay == 2");
333 gPad->SetFillColor(33);
334 ntuple2->SetFillColor(42);
335 ntuple2->Draw("x","lay==1");
337 gPad->SetFillColor(33);
338 ntuple2->SetFillColor(46);
339 ntuple2->Draw("z","lay==1");
341 gPad->SetFillColor(33);
342 ntuple2->SetFillColor(42);
343 ntuple2->Draw("x","lay==2");
345 gPad->SetFillColor(33);
346 ntuple2->SetFillColor(46);
347 ntuple2->Draw("z","lay==2");
352 gPad->SetFillColor(33);
353 ntuple2->SetFillColor(42);
354 ntuple2->Draw("nx","lay==1");
356 gPad->SetFillColor(33);
357 ntuple2->SetFillColor(46);
358 ntuple2->Draw("nz","lay==1");
360 gPad->SetFillColor(33);
361 ntuple2->SetFillColor(42);
362 ntuple2->Draw("nx","lay==2");
364 gPad->SetFillColor(33);
365 ntuple2->SetFillColor(46);
366 ntuple2->Draw("nz","lay==2");
371 gPad->SetFillColor(33);
372 ntuple2->SetFillColor(42);
373 ntuple2->Draw("qcl","lay==1");
375 gPad->SetFillColor(33);
376 ntuple2->SetFillColor(46);
377 ntuple2->Draw("nx","lay==1");
379 gPad->SetFillColor(33);
380 ntuple2->SetFillColor(42);
381 ntuple2->Draw("qcl","lay==2");
383 gPad->SetFillColor(33);
384 ntuple2->SetFillColor(46);
385 ntuple2->Draw("nz","lay==2");
388 ///////////////////// Histogramm analysis ////////////////////////
392 gPad->SetFillColor(33);
393 Xres1->SetFillColor(42);
396 gPad->SetFillColor(33);
397 Zres1->SetFillColor(42);
400 gPad->SetFillColor(33);
401 Xres2->SetFillColor(46);
404 gPad->SetFillColor(33);
405 Zres2->SetFillColor(46);
411 gPad->SetFillColor(33);
412 Nzpix1->SetFillColor(42);
415 gPad->SetFillColor(33);
416 Nxpix1->SetFillColor(42);
419 gPad->SetFillColor(33);
420 Nzpix2->SetFillColor(46);
423 gPad->SetFillColor(33);
424 Nxpix2->SetFillColor(46);
430 gPad->SetFillColor(33);
431 Zpix1->SetFillColor(42);
434 gPad->SetFillColor(33);
435 Xpix1->SetFillColor(46);
438 gPad->SetFillColor(33);
439 Zpix2->SetFillColor(42);
442 gPad->SetFillColor(33);
443 Xpix2->SetFillColor(46);
450 gPad->SetFillColor(33);
451 Theta1->SetFillColor(42);
454 gPad->SetFillColor(33);
455 Phi1->SetFillColor(46);
458 gPad->SetFillColor(33);
459 Theta2->SetFillColor(42);
462 gPad->SetFillColor(33);
463 Phi2->SetFillColor(46);
469 gPad->SetFillColor(33);
470 Ptot1->SetFillColor(42);
473 gPad->SetFillColor(33);
474 Pz1->SetFillColor(46);
477 gPad->SetFillColor(33);
478 Ptot2->SetFillColor(42);
481 gPad->SetFillColor(33);
482 Pz2->SetFillColor(46);
488 gPad->SetFillColor(33);
489 Eta1->SetFillColor(42);
492 gPad->SetFillColor(33);
493 Y1->SetFillColor(46);
496 gPad->SetFillColor(33);
497 Eta2->SetFillColor(42);
500 gPad->SetFillColor(33);
501 Y2->SetFillColor(46);
507 gPad->SetFillColor(33);
508 Eta1Den->SetFillColor(42);
511 gPad->SetFillColor(33);
512 Y1Den->SetFillColor(46);
515 gPad->SetFillColor(33);
516 Eta2Den->SetFillColor(42);
519 gPad->SetFillColor(33);
520 Y2Den->SetFillColor(46);
526 gPad->SetFillColor(33);
527 Eta1DenA->SetFillColor(42);
530 gPad->SetFillColor(33);
531 Y1DenA->SetFillColor(46);
534 gPad->SetFillColor(33);
535 Eta2DenA->SetFillColor(42);
538 gPad->SetFillColor(33);
539 Y2DenA->SetFillColor(46);
545 c2->Print("spd_res.ps");
549 c2->Print("spd_clsize.ps");