]> git.uio.no Git - u/mrichter/AliRoot.git/blame - ITS/SPD_ntuple.C
git://git.uio.no / u / mrichter / AliRoot.git / blame
? search:
summary | shortlog | log | commit | commitdiff | tree
blob | blame (incremental) | history | HEAD
Remove some shielding to accomodate compensator magnet.
[u/mrichter/AliRoot.git] / ITS / SPD_ntuple.C
CommitLineData
e8189707 1void SPD_ntuple()
2{
3
4TFile *f = new TFile("SPD_his.root");
5//TH1F *Pix = (TH1F*)f->Get("Pix");
6
7//gStyle->SetOptStat(1111111);
536aa5ec 8//gStyle->SetOptLogy();
e8189707 9//TCanvas *c1 = new TCanvas("c1","SPD clusters",400,10,600,700);
10TCanvas *c2 = new TCanvas("c2","SPD clusters",400,10,600,700);
11//c1->Divide(2,2);
12c2->Divide(2,2);
e8189707 13
14///////////////////////// Ntuple analysis ///////////////////////////////
15
16// ntuple is created inside the hit loop for the hits in the cluster region;
17
18// ntuple1 is created after a finish of the hit loop if one or more hits
19// are in the cluster region;
20
21// ntuple2 is created befor the hit loop for all clusters;
22
23// -----------------------------------------------------------------------
24// lay - number of ITS layer;
536aa5ec 25// x - coordinates in r*phi(x) direction (mm);
26// z - coordinates in z direction (mm);
e8189707 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;
536aa5ec 34// noverlaps - number of particles overlapping in one cluster found in this
e2ce6aca 35// macros in the cluster region;
536aa5ec 36// ntrover - number of particles overlapping in one cluster found in the
e2ce6aca 37// AliITSClusterFinderSPD class;
e8189707 38// noverprim - number of primary particles overlapping in one cluster;
39// qcl - cluster charge in electron number
40// -------------------------------------------------------------------------
41
536aa5ec 42/*
43c2->cd(1);
44gPad->SetFillColor(33);
45 ntuple->SetFillColor(42);
46 ntuple->Draw("dx","lay == 1");
47c2->cd(2);
48gPad->SetFillColor(33);
49 ntuple->SetFillColor(46);
50 ntuple->Draw("dz","lay == 1");
51c2->cd(3);
52gPad->SetFillColor(33);
53 ntuple->SetFillColor(42);
54 ntuple->Draw("dx","lay == 2");
55c2->cd(4);
56gPad->SetFillColor(33);
57 ntuple->SetFillColor(46);
58 ntuple->Draw("dz","lay == 2");
59*/
60
e8189707 61/*
62c2->cd(1);
63gPad->SetFillColor(33);
64 ntuple->SetFillColor(42);
65 ntuple->Draw("dx","lay == 1 && hitprim == 1");
66c2->cd(2);
67gPad->SetFillColor(33);
68 ntuple->SetFillColor(46);
69 ntuple->Draw("dz","lay == 1 && hitprim == 1");
70c2->cd(3);
71gPad->SetFillColor(33);
72 ntuple->SetFillColor(42);
73 ntuple->Draw("dx","lay == 2 && hitprim == 1");
74c2->cd(4);
75gPad->SetFillColor(33);
76 ntuple->SetFillColor(46);
77 ntuple->Draw("dz","lay == 2 && hitprim == 1");
78*/
79
536aa5ec 80/*
81c2->cd(1);
82gPad->SetFillColor(33);
83 ntuple->SetFillColor(42);
bf83d539 84 ntuple->Draw("dx","lay == 1 && hitprim == 1 && ntrover>=1&&nx==4");
536aa5ec 85c2->cd(2);
86gPad->SetFillColor(33);
87 ntuple->SetFillColor(46);
bf83d539 88 ntuple->Draw("dz","lay == 1 && hitprim == 1 && ntrover>=1&&nz==4");
536aa5ec 89c2->cd(3);
90gPad->SetFillColor(33);
91 ntuple->SetFillColor(42);
bf83d539 92 ntuple->Draw("dx","lay == 2 && hitprim == 1 && ntrover>=1&&nx==4");
536aa5ec 93c2->cd(4);
94gPad->SetFillColor(33);
95 ntuple->SetFillColor(46);
bf83d539 96 ntuple->Draw("dz","lay == 2 && hitprim == 1 && ntrover>=1&&nz==4");
536aa5ec 97*/
98
e2ce6aca 99
e8189707 100c2->cd(1);
101gPad->SetFillColor(33);
102 ntuple->SetFillColor(42);
103 ntuple->Draw("dx","lay==1&&hitprim==1&&dx>-150&&dx<150");
104c2->cd(2);
105gPad->SetFillColor(33);
106 ntuple->SetFillColor(46);
107 ntuple->Draw("dz","lay==1&&hitprim==1&&dz>-500&&dz<500");
108c2->cd(3);
109gPad->SetFillColor(33);
110 ntuple->SetFillColor(42);
111 ntuple->Draw("dx","lay==2&&hitprim==1&&dx>-150&&dx<150");
112c2->cd(4);
113gPad->SetFillColor(33);
114 ntuple->SetFillColor(46);
115 ntuple->Draw("dz","lay==2&&hitprim==1&&dz>-500&&dz<500");
bf83d539 116
536aa5ec 117
118/*
119c2->cd(1);
120gPad->SetFillColor(33);
121 ntuple->SetFillColor(42);
122 ntuple->Draw("dx","lay==1&&hitprim==1&&dx>-150&&dx<150&&ntrover==1");
123c2->cd(2);
124gPad->SetFillColor(33);
125 ntuple->SetFillColor(46);
126 ntuple->Draw("dz","lay==1&&hitprim==1&&dz>-500&&dz<500&&ntrover==1");
127c2->cd(3);
128gPad->SetFillColor(33);
129 ntuple->SetFillColor(42);
130 ntuple->Draw("dx","lay==2&&hitprim==1&&dx>-150&&dx<150&&ntrover==1");
131c2->cd(4);
132gPad->SetFillColor(33);
133 ntuple->SetFillColor(46);
134 ntuple->Draw("dz","lay==2&&hitprim==1&&dz>-500&&dz<500&&ntrover==1");
135*/
e2ce6aca 136
bf83d539 137 /*
e8189707 138c2->cd(1);
139gPad->SetFillColor(33);
140 ntuple->SetFillColor(42);
bf83d539 141 ntuple->Draw("dx","lay==1&&hitprim==1&&nx>4");
e8189707 142c2->cd(2);
143gPad->SetFillColor(33);
144 ntuple->SetFillColor(46);
bf83d539 145 ntuple->Draw("dz","lay==1&&hitprim==1&&nz>4");
e8189707 146c2->cd(3);
147gPad->SetFillColor(33);
148 ntuple->SetFillColor(42);
bf83d539 149 ntuple->Draw("dx","lay==2&&hitprim==1&&nx>4");
e8189707 150c2->cd(4);
151gPad->SetFillColor(33);
152 ntuple->SetFillColor(46);
bf83d539 153 ntuple->Draw("dz","lay==2&&hitprim==1&&nz>4");
154 */
e8189707 155
536aa5ec 156/*
e8189707 157c2->cd(1);
158gPad->SetFillColor(33);
159 ntuple1->SetFillColor(42);
bf83d539 160 ntuple1->Draw("noverlaps","lay == 1&&noverprim>0&&noverlaps>0");
e8189707 161c2->cd(2);
162gPad->SetFillColor(33);
163 ntuple1->SetFillColor(46);
bf83d539 164 ntuple1->Draw("noverlaps","lay == 2&&noverprim>0&&noverlaps>0");
e8189707 165c2->cd(3);
166gPad->SetFillColor(33);
167 ntuple1->SetFillColor(42);
bf83d539 168 ntuple1->Draw("ntrover","lay == 1&&noverprim>0&&ntrover>0");
e8189707 169c2->cd(4);
170gPad->SetFillColor(33);
171 ntuple1->SetFillColor(46);
bf83d539 172 ntuple1->Draw("ntrover","lay == 2&&noverprim>0&&ntrover>0");
536aa5ec 173*/
e2ce6aca 174
bf83d539 175 /*
e8189707 176c2->cd(1);
177gPad->SetFillColor(33);
178 ntuple1->SetFillColor(42);
179 ntuple1->Draw("noverprim","lay == 1");
180c2->cd(2);
181gPad->SetFillColor(33);
182 ntuple1->SetFillColor(46);
183 ntuple1->Draw("noverprim","lay == 2");
184c2->cd(3);
185gPad->SetFillColor(33);
186 ntuple1->SetFillColor(42);
536aa5ec 187 ntuple1->Draw("noverprim","lay == 1 && noverprim < 1");
e8189707 188c2->cd(4);
189gPad->SetFillColor(33);
190 ntuple1->SetFillColor(46);
536aa5ec 191 ntuple1->Draw("noverprim","lay == 2 && noverprim < 1");
bf83d539 192 */
e2ce6aca 193
e8189707 194/*
195c2->cd(1);
196gPad->SetFillColor(33);
197 ntuple1->SetFillColor(42);
536aa5ec 198 ntuple1->Draw("nx:nz","lay == 1");
e8189707 199c2->cd(2);
200gPad->SetFillColor(33);
201 ntuple1->SetFillColor(46);
536aa5ec 202 ntuple1->Draw("nx:nz","lay == 2");
e8189707 203*/
204
bf83d539 205 /*
536aa5ec 206c2->cd(1);
207gPad->SetFillColor(33);
208 ntuple1->SetFillColor(42);
209 ntuple1->Draw("qcl","lay == 1");
210c2->cd(2);
211gPad->SetFillColor(33);
212 ntuple1->SetFillColor(46);
213 ntuple1->Draw("qcl","lay == 2");
bf83d539 214 */
536aa5ec 215
216/*
217c2->cd(1);
218gPad->SetFillColor(33);
219 ntuple1->SetFillColor(42);
220 ntuple1->Draw("nx","lay == 1");
221c2->cd(2);
222gPad->SetFillColor(33);
223 ntuple1->SetFillColor(46);
224 ntuple1->Draw("nx","lay == 2");
225c2->cd(3);
226gPad->SetFillColor(33);
227 ntuple1->SetFillColor(42);
228 ntuple1->Draw("nz","lay == 1");
229c2->cd(4);
230gPad->SetFillColor(33);
231 ntuple1->SetFillColor(46);
232 ntuple1->Draw("nz","lay == 2");
233*/
234
235 /*
e8189707 236c2->cd(1);
237gPad->SetFillColor(33);
238 ntuple1->SetFillColor(42);
239 ntuple1->Draw("nx","lay == 1&&noverprim>0");
240c2->cd(2);
241gPad->SetFillColor(33);
242 ntuple1->SetFillColor(46);
243 ntuple1->Draw("nx","lay == 2&&noverprim>0");
244c2->cd(3);
245gPad->SetFillColor(33);
246 ntuple1->SetFillColor(42);
247 ntuple1->Draw("nz","lay == 1 && noverprim>0");
248c2->cd(4);
249gPad->SetFillColor(33);
250 ntuple1->SetFillColor(46);
251 ntuple1->Draw("nz","lay == 2 && noverprim>0");
536aa5ec 252 */
253
254/*
255c2->cd(1);
256gPad->SetFillColor(33);
257 ntuple1->SetFillColor(42);
258 ntuple1->Draw("x","lay == 1&&noverprim>=0");
259c2->cd(2);
260gPad->SetFillColor(33);
261 ntuple1->SetFillColor(46);
262 ntuple1->Draw("x","lay == 2&&noverprim>=0");
263c2->cd(3);
264gPad->SetFillColor(33);
265 ntuple1->SetFillColor(42);
266 ntuple1->Draw("z","lay == 1 && noverprim>=0");
267c2->cd(4);
268gPad->SetFillColor(33);
269 ntuple1->SetFillColor(46);
270 ntuple1->Draw("z","lay == 2 && noverprim>=0");
e2ce6aca 271*/
536aa5ec 272
273 /*
274c2->cd(1);
275gPad->SetFillColor(33);
276 ntuple1->SetFillColor(42);
277 ntuple1->Draw("nx","lay == 1&&noverprim>0 && ntrover == 1");
278c2->cd(2);
279gPad->SetFillColor(33);
280 ntuple1->SetFillColor(46);
281 ntuple1->Draw("nx","lay == 2&&noverprim>0 && ntrover == 1");
282c2->cd(3);
283gPad->SetFillColor(33);
284 ntuple1->SetFillColor(42);
285 ntuple1->Draw("nz","lay == 1 && noverprim>0 && ntrover == 1");
286c2->cd(4);
287gPad->SetFillColor(33);
288 ntuple1->SetFillColor(46);
289 ntuple1->Draw("nz","lay == 2 && noverprim>0 && ntrover ==1");
290 */
e8189707 291
292/*
293c2->cd(1);
294gPad->SetFillColor(33);
295 ntuple1->SetFillColor(42);
296 ntuple1->Draw("dx","lay == 1&&noverprim>0");
297c2->cd(2);
298gPad->SetFillColor(33);
299 ntuple1->SetFillColor(46);
300 ntuple1->Draw("dx","lay == 2&&noverprim>0");
301c2->cd(3);
302gPad->SetFillColor(33);
303 ntuple1->SetFillColor(42);
304 ntuple1->Draw("dz","lay == 1 && noverprim>0");
305c2->cd(4);
306gPad->SetFillColor(33);
307 ntuple1->SetFillColor(46);
308 ntuple1->Draw("dz","lay == 2 && noverprim>0");
309*/
310
311 /*
312c2->cd(1);
313gPad->SetFillColor(33);
314 ntuple1->SetFillColor(42);
315 ntuple1->Draw("dx","lay == 1");
316c2->cd(2);
317gPad->SetFillColor(33);
318 ntuple1->SetFillColor(46);
319 ntuple1->Draw("dx","lay == 2");
320c2->cd(3);
321gPad->SetFillColor(33);
322 ntuple1->SetFillColor(42);
323 ntuple1->Draw("dz","lay == 1");
324c2->cd(4);
325gPad->SetFillColor(33);
326 ntuple1->SetFillColor(46);
327 ntuple1->Draw("dz","lay == 2");
328 */
e2ce6aca 329
bf83d539 330 /*
e8189707 331c2->cd(1);
332gPad->SetFillColor(33);
333 ntuple2->SetFillColor(42);
536aa5ec 334 ntuple2->Draw("x","lay==1");
e8189707 335c2->cd(2);
336gPad->SetFillColor(33);
337 ntuple2->SetFillColor(46);
536aa5ec 338 ntuple2->Draw("z","lay==1");
e8189707 339c2->cd(3);
340gPad->SetFillColor(33);
341 ntuple2->SetFillColor(42);
536aa5ec 342 ntuple2->Draw("x","lay==2");
e8189707 343c2->cd(4);
344gPad->SetFillColor(33);
345 ntuple2->SetFillColor(46);
536aa5ec 346 ntuple2->Draw("z","lay==2");
bf83d539 347 */
536aa5ec 348
349 /*
350c2->cd(1);
351gPad->SetFillColor(33);
352 ntuple2->SetFillColor(42);
353 ntuple2->Draw("nx","lay==1");
354c2->cd(2);
355gPad->SetFillColor(33);
356 ntuple2->SetFillColor(46);
357 ntuple2->Draw("nz","lay==1");
358c2->cd(3);
359gPad->SetFillColor(33);
360 ntuple2->SetFillColor(42);
361 ntuple2->Draw("nx","lay==2");
362c2->cd(4);
363gPad->SetFillColor(33);
364 ntuple2->SetFillColor(46);
365 ntuple2->Draw("nz","lay==2");
e8189707 366 */
536aa5ec 367
368/*
369c2->cd(1);
370gPad->SetFillColor(33);
371 ntuple2->SetFillColor(42);
372 ntuple2->Draw("qcl","lay==1");
373c2->cd(2);
374gPad->SetFillColor(33);
375 ntuple2->SetFillColor(46);
376 ntuple2->Draw("nx","lay==1");
377c2->cd(3);
378gPad->SetFillColor(33);
379 ntuple2->SetFillColor(42);
380 ntuple2->Draw("qcl","lay==2");
381c2->cd(4);
382gPad->SetFillColor(33);
383 ntuple2->SetFillColor(46);
384 ntuple2->Draw("nz","lay==2");
385*/
e8189707 386
387///////////////////// Histogramm analysis ////////////////////////
388
536aa5ec 389/*
e8189707 390c2->cd(1);
391gPad->SetFillColor(33);
392 Xres1->SetFillColor(42);
393 Xres1->Draw();
394c2->cd(2);
395gPad->SetFillColor(33);
396 Zres1->SetFillColor(42);
397 Zres1->Draw();
398c2->cd(3);
399gPad->SetFillColor(33);
400 Xres2->SetFillColor(46);
401 Xres2->Draw();
402c2->cd(4);
403gPad->SetFillColor(33);
404 Zres2->SetFillColor(46);
405 Zres2->Draw();
536aa5ec 406*/
407
e8189707 408 /*
409c2->cd(1);
410gPad->SetFillColor(33);
411 Nzpix1->SetFillColor(42);
412 Nzpix1->Draw();
413c2->cd(2);
414gPad->SetFillColor(33);
415 Nxpix1->SetFillColor(42);
416 Nxpix1->Draw();
417c2->cd(3);
418gPad->SetFillColor(33);
419 Nzpix2->SetFillColor(46);
420 Nzpix2->Draw();
421c2->cd(4);
422gPad->SetFillColor(33);
423 Nxpix2->SetFillColor(46);
424 Nxpix2->Draw();
425 */
426
427 /*
428c2->cd(1);
429gPad->SetFillColor(33);
430 Zpix1->SetFillColor(42);
431 Zpix1->Draw();
432c2->cd(2);
433gPad->SetFillColor(33);
434 Xpix1->SetFillColor(46);
435 Xpix1->Draw();
436c2->cd(3);
437gPad->SetFillColor(33);
438 Zpix2->SetFillColor(42);
439 Zpix2->Draw();
440c2->cd(4);
441gPad->SetFillColor(33);
442 Xpix2->SetFillColor(46);
443 Xpix2->Draw();
444 */
536aa5ec 445
446
447/*
448c2->cd(1);
449gPad->SetFillColor(33);
450 Theta1->SetFillColor(42);
451 Theta1->Draw();
452c2->cd(2);
453gPad->SetFillColor(33);
454 Phi1->SetFillColor(46);
455 Phi1->Draw();
456c2->cd(3);
457gPad->SetFillColor(33);
458 Theta2->SetFillColor(42);
459 Theta2->Draw();
460c2->cd(4);
461gPad->SetFillColor(33);
462 Phi2->SetFillColor(46);
463 Phi2->Draw();
464*/
465
466/*
467c2->cd(1);
468gPad->SetFillColor(33);
469 Ptot1->SetFillColor(42);
470 Ptot1->Draw();
471c2->cd(2);
472gPad->SetFillColor(33);
473 Pz1->SetFillColor(46);
474 Pz1->Draw();
475c2->cd(3);
476gPad->SetFillColor(33);
477 Ptot2->SetFillColor(42);
478 Ptot2->Draw();
479c2->cd(4);
480gPad->SetFillColor(33);
481 Pz2->SetFillColor(46);
482 Pz2->Draw();
483*/
484
485/*
486c2->cd(1);
487gPad->SetFillColor(33);
488 Eta1->SetFillColor(42);
489 Eta1->Draw();
490c2->cd(2);
491gPad->SetFillColor(33);
492 Y1->SetFillColor(46);
493 Y1->Draw();
494c2->cd(3);
495gPad->SetFillColor(33);
496 Eta2->SetFillColor(42);
497 Eta2->Draw();
498c2->cd(4);
499gPad->SetFillColor(33);
500 Y2->SetFillColor(46);
501 Y2->Draw();
502*/
503
504 /*
505c2->cd(1);
506gPad->SetFillColor(33);
507 Eta1Den->SetFillColor(42);
508 Eta1Den->Draw();
509c2->cd(2);
510gPad->SetFillColor(33);
511 Y1Den->SetFillColor(46);
512 Y1Den->Draw();
513c2->cd(3);
514gPad->SetFillColor(33);
515 Eta2Den->SetFillColor(42);
516 Eta2Den->Draw();
517c2->cd(4);
518gPad->SetFillColor(33);
519 Y2Den->SetFillColor(46);
520 Y2Den->Draw();
521 */
522
523 /*
524c2->cd(1);
525gPad->SetFillColor(33);
526 Eta1DenA->SetFillColor(42);
527 Eta1DenA->Draw();
528c2->cd(2);
529gPad->SetFillColor(33);
530 Y1DenA->SetFillColor(46);
531 Y1DenA->Draw();
532c2->cd(3);
533gPad->SetFillColor(33);
534 Eta2DenA->SetFillColor(42);
535 Eta2DenA->Draw();
536c2->cd(4);
537gPad->SetFillColor(33);
538 Y2DenA->SetFillColor(46);
539 Y2DenA->Draw();
540 */
e8189707 541
542 /*
543c2->Draw();
544c2->Print("spd_res.ps");
545 */
546 /*
547c2->Draw();
548c2->Print("spd_clsize.ps");
549 */
550
551}
552
553
554
555
556
557
558
559