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