New class for ITS coordiante transformations used by AliITSgeom nearly
[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);
8gStyle->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//c2->SetLogy();
14
15///////////////////////// Ntuple analysis ///////////////////////////////
16
17// ntuple is created inside the hit loop for the hits in the cluster region;
18
19// ntuple1 is created after a finish of the hit loop if one or more hits
20// are in the cluster region;
21
22// ntuple2 is created befor the hit loop for all clusters;
23
24// -----------------------------------------------------------------------
25// lay - number of ITS layer;
26// nx - cluster size in the r*phi(x) direction;
27// nz - cluster size in the z direction;
28// hitprim - primary particle(hit) flag ( = 1 for primery particle);
29// dx - difference of hit(mediate) and reconstructed (from cluster)
30// coordinates in r*phi(x) direction;
31// dz - difference of hit(mediate) and reconstructed (from cluster)
32// coordinates in z direction;
33// noverlaps - number of particles overlapping in one cluster;
34// noverprim - number of primary particles overlapping in one cluster;
35// qcl - cluster charge in electron number
36// -------------------------------------------------------------------------
37
38/*
39c2->cd(1);
40gPad->SetFillColor(33);
41 ntuple->SetFillColor(42);
42 ntuple->Draw("dx","lay == 1 && hitprim == 1");
43c2->cd(2);
44gPad->SetFillColor(33);
45 ntuple->SetFillColor(46);
46 ntuple->Draw("dz","lay == 1 && hitprim == 1");
47c2->cd(3);
48gPad->SetFillColor(33);
49 ntuple->SetFillColor(42);
50 ntuple->Draw("dx","lay == 2 && hitprim == 1");
51c2->cd(4);
52gPad->SetFillColor(33);
53 ntuple->SetFillColor(46);
54 ntuple->Draw("dz","lay == 2 && hitprim == 1");
55*/
56
57/*
58c2->cd(1);
59gPad->SetFillColor(33);
60 ntuple->SetFillColor(42);
61 ntuple->Draw("dx","lay == 1 && hitprim == 1");
62c2->cd(2);
63gPad->SetFillColor(33);
64 ntuple->SetFillColor(46);
65 ntuple->Draw("dx","lay == 1 && hitprim == 1 && nx < 15");
66c2->cd(3);
67gPad->SetFillColor(33);
68 ntuple->SetFillColor(42);
69 ntuple->Draw("dx","lay == 1 && hitprim == 1 && nx < 5");
70c2->cd(4);
71gPad->SetFillColor(33);
72 ntuple->SetFillColor(46);
73 ntuple->Draw("dx","lay == 1 && hitprim == 1 && nx < 4");
74*/
75
76 /*
77c2->cd(1);
78gPad->SetFillColor(33);
79 ntuple->SetFillColor(42);
80 ntuple->Draw("dx","lay==1&&hitprim==1&&dx>-150&&dx<150");
81c2->cd(2);
82gPad->SetFillColor(33);
83 ntuple->SetFillColor(46);
84 ntuple->Draw("dz","lay==1&&hitprim==1&&dz>-500&&dz<500");
85c2->cd(3);
86gPad->SetFillColor(33);
87 ntuple->SetFillColor(42);
88 ntuple->Draw("dx","lay==2&&hitprim==1&&dx>-150&&dx<150");
89c2->cd(4);
90gPad->SetFillColor(33);
91 ntuple->SetFillColor(46);
92 ntuple->Draw("dz","lay==2&&hitprim==1&&dz>-500&&dz<500");
93 */
94/*
95c2->cd(1);
96gPad->SetFillColor(33);
97 ntuple->SetFillColor(42);
98 ntuple->Draw("dx","lay==1&&hitprim==1&&dx>-150&&dx<150&&nx<5");
99c2->cd(2);
100gPad->SetFillColor(33);
101 ntuple->SetFillColor(46);
102 ntuple->Draw("dz","lay==1&&hitprim==1&&dz>-500&&dz<500&&nz<5");
103c2->cd(3);
104gPad->SetFillColor(33);
105 ntuple->SetFillColor(42);
106 ntuple->Draw("dx","lay==2&&hitprim==1&&dx>-150&&dx<150&&nx<5");
107c2->cd(4);
108gPad->SetFillColor(33);
109 ntuple->SetFillColor(46);
110 ntuple->Draw("dz","lay==2&&hitprim==1&&dz>-500&&dz<500&&nz<5");
111*/
112/*
113c2->cd(1);
114gPad->SetFillColor(33);
115 ntuple->SetFillColor(42);
116 ntuple->Draw("dx","lay==1&&hitprim==1&&nx<4");
117c2->cd(2);
118gPad->SetFillColor(33);
119 ntuple->SetFillColor(46);
120 ntuple->Draw("dz","lay==1&&hitprim==1&&nz<4");
121c2->cd(3);
122gPad->SetFillColor(33);
123 ntuple->SetFillColor(42);
124 ntuple->Draw("dx","lay==2&&hitprim==1&&nx<4");
125c2->cd(4);
126gPad->SetFillColor(33);
127 ntuple->SetFillColor(46);
128 ntuple->Draw("dz","lay==2&&hitprim==1&&nz<4");
129*/
130
131/*
132c2->cd(1);
133gPad->SetFillColor(33);
134 ntuple1->SetFillColor(42);
135 ntuple1->Draw("noverlaps","lay == 1");
136c2->cd(2);
137gPad->SetFillColor(33);
138 ntuple1->SetFillColor(46);
139 ntuple1->Draw("noverlaps","lay == 2");
140
141c2->cd(3);
142gPad->SetFillColor(33);
143 ntuple1->SetFillColor(42);
144 ntuple1->Draw("noverlaps","lay == 1 && noverprim > 0");
145c2->cd(4);
146gPad->SetFillColor(33);
147 ntuple1->SetFillColor(46);
148 ntuple1->Draw("noverlaps","lay == 2 && noverprim > 0");
149*/
150/*
151c2->cd(1);
152gPad->SetFillColor(33);
153 ntuple1->SetFillColor(42);
154 ntuple1->Draw("noverprim","lay == 1");
155c2->cd(2);
156gPad->SetFillColor(33);
157 ntuple1->SetFillColor(46);
158 ntuple1->Draw("noverprim","lay == 2");
159c2->cd(3);
160gPad->SetFillColor(33);
161 ntuple1->SetFillColor(42);
162 ntuple1->Draw("noverprim","lay == 1 && noverprim > 1");
163c2->cd(4);
164gPad->SetFillColor(33);
165 ntuple1->SetFillColor(46);
166 ntuple1->Draw("noverprim","lay == 2 && noverprim > 1");
167*/
168/*
169c2->cd(1);
170gPad->SetFillColor(33);
171 ntuple1->SetFillColor(42);
172 ntuple1->Draw("qcl","lay == 1&&noverprim>0");
173c2->cd(2);
174gPad->SetFillColor(33);
175 ntuple1->SetFillColor(46);
176 ntuple1->Draw("qcl","lay == 2&&noverprim>0");
177*/
178
179
180c2->cd(1);
181gPad->SetFillColor(33);
182 ntuple1->SetFillColor(42);
183 ntuple1->Draw("nx","lay == 1&&noverprim>0");
184c2->cd(2);
185gPad->SetFillColor(33);
186 ntuple1->SetFillColor(46);
187 ntuple1->Draw("nx","lay == 2&&noverprim>0");
188c2->cd(3);
189gPad->SetFillColor(33);
190 ntuple1->SetFillColor(42);
191 ntuple1->Draw("nz","lay == 1 && noverprim>0");
192c2->cd(4);
193gPad->SetFillColor(33);
194 ntuple1->SetFillColor(46);
195 ntuple1->Draw("nz","lay == 2 && noverprim>0");
196
197
198/*
199c2->cd(1);
200gPad->SetFillColor(33);
201 ntuple1->SetFillColor(42);
202 ntuple1->Draw("dx","lay == 1&&noverprim>0");
203c2->cd(2);
204gPad->SetFillColor(33);
205 ntuple1->SetFillColor(46);
206 ntuple1->Draw("dx","lay == 2&&noverprim>0");
207c2->cd(3);
208gPad->SetFillColor(33);
209 ntuple1->SetFillColor(42);
210 ntuple1->Draw("dz","lay == 1 && noverprim>0");
211c2->cd(4);
212gPad->SetFillColor(33);
213 ntuple1->SetFillColor(46);
214 ntuple1->Draw("dz","lay == 2 && noverprim>0");
215*/
216
217 /*
218c2->cd(1);
219gPad->SetFillColor(33);
220 ntuple1->SetFillColor(42);
221 ntuple1->Draw("dx","lay == 1");
222c2->cd(2);
223gPad->SetFillColor(33);
224 ntuple1->SetFillColor(46);
225 ntuple1->Draw("dx","lay == 2");
226c2->cd(3);
227gPad->SetFillColor(33);
228 ntuple1->SetFillColor(42);
229 ntuple1->Draw("dz","lay == 1");
230c2->cd(4);
231gPad->SetFillColor(33);
232 ntuple1->SetFillColor(46);
233 ntuple1->Draw("dz","lay == 2");
234 */
235 /*
236c2->cd(1);
237gPad->SetFillColor(33);
238 ntuple2->SetFillColor(42);
239 ntuple2->Draw("nx","lay == 1");
240c2->cd(2);
241gPad->SetFillColor(33);
242 ntuple2->SetFillColor(46);
243 ntuple2->Draw("nz","lay == 1");
244c2->cd(3);
245gPad->SetFillColor(33);
246 ntuple2->SetFillColor(42);
247 ntuple2->Draw("nx","lay == 2");
248c2->cd(4);
249gPad->SetFillColor(33);
250 ntuple2->SetFillColor(46);
251 ntuple2->Draw("nz","lay == 2");
252 */
253
254///////////////////// Histogramm analysis ////////////////////////
255
256 /*
257c2->cd(1);
258gPad->SetFillColor(33);
259 Xres1->SetFillColor(42);
260 Xres1->Draw();
261c2->cd(2);
262gPad->SetFillColor(33);
263 Zres1->SetFillColor(42);
264 Zres1->Draw();
265c2->cd(3);
266gPad->SetFillColor(33);
267 Xres2->SetFillColor(46);
268 Xres2->Draw();
269c2->cd(4);
270gPad->SetFillColor(33);
271 Zres2->SetFillColor(46);
272 Zres2->Draw();
273 */
274 /*
275c2->cd(1);
276gPad->SetFillColor(33);
277 Nzpix1->SetFillColor(42);
278 Nzpix1->Draw();
279c2->cd(2);
280gPad->SetFillColor(33);
281 Nxpix1->SetFillColor(42);
282 Nxpix1->Draw();
283c2->cd(3);
284gPad->SetFillColor(33);
285 Nzpix2->SetFillColor(46);
286 Nzpix2->Draw();
287c2->cd(4);
288gPad->SetFillColor(33);
289 Nxpix2->SetFillColor(46);
290 Nxpix2->Draw();
291 */
292
293 /*
294c2->cd(1);
295gPad->SetFillColor(33);
296 Zpix1->SetFillColor(42);
297 Zpix1->Draw();
298c2->cd(2);
299gPad->SetFillColor(33);
300 Xpix1->SetFillColor(46);
301 Xpix1->Draw();
302c2->cd(3);
303gPad->SetFillColor(33);
304 Zpix2->SetFillColor(42);
305 Zpix2->Draw();
306c2->cd(4);
307gPad->SetFillColor(33);
308 Xpix2->SetFillColor(46);
309 Xpix2->Draw();
310 */
311
312 /*
313c2->Draw();
314c2->Print("spd_res.ps");
315 */
316 /*
317c2->Draw();
318c2->Print("spd_clsize.ps");
319 */
320
321}
322
323
324
325
326
327
328
329