]> git.uio.no Git - u/mrichter/AliRoot.git/blame - ITS/SPD_ntuple.C
Containers definition
[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;
e2ce6aca 33// noverlaps - number of particles overlapping in one cluster finded in this
34// macros in the cluster region;
35// ntrover - number of particles overlapping in one cluster finded in the
36// AliITSClusterFinderSPD class;
e8189707 37// noverprim - number of primary particles overlapping in one cluster;
38// qcl - cluster charge in electron number
39// -------------------------------------------------------------------------
40
41/*
42c2->cd(1);
43gPad->SetFillColor(33);
44 ntuple->SetFillColor(42);
45 ntuple->Draw("dx","lay == 1 && hitprim == 1");
46c2->cd(2);
47gPad->SetFillColor(33);
48 ntuple->SetFillColor(46);
49 ntuple->Draw("dz","lay == 1 && hitprim == 1");
50c2->cd(3);
51gPad->SetFillColor(33);
52 ntuple->SetFillColor(42);
53 ntuple->Draw("dx","lay == 2 && hitprim == 1");
54c2->cd(4);
55gPad->SetFillColor(33);
56 ntuple->SetFillColor(46);
57 ntuple->Draw("dz","lay == 2 && hitprim == 1");
58*/
59
60/*
61c2->cd(1);
62gPad->SetFillColor(33);
63 ntuple->SetFillColor(42);
64 ntuple->Draw("dx","lay == 1 && hitprim == 1");
65c2->cd(2);
66gPad->SetFillColor(33);
67 ntuple->SetFillColor(46);
68 ntuple->Draw("dx","lay == 1 && hitprim == 1 && nx < 15");
69c2->cd(3);
70gPad->SetFillColor(33);
71 ntuple->SetFillColor(42);
72 ntuple->Draw("dx","lay == 1 && hitprim == 1 && nx < 5");
73c2->cd(4);
74gPad->SetFillColor(33);
75 ntuple->SetFillColor(46);
76 ntuple->Draw("dx","lay == 1 && hitprim == 1 && nx < 4");
77*/
78
e2ce6aca 79
e8189707 80 /*
81c2->cd(1);
82gPad->SetFillColor(33);
83 ntuple->SetFillColor(42);
84 ntuple->Draw("dx","lay==1&&hitprim==1&&dx>-150&&dx<150");
85c2->cd(2);
86gPad->SetFillColor(33);
87 ntuple->SetFillColor(46);
88 ntuple->Draw("dz","lay==1&&hitprim==1&&dz>-500&&dz<500");
89c2->cd(3);
90gPad->SetFillColor(33);
91 ntuple->SetFillColor(42);
92 ntuple->Draw("dx","lay==2&&hitprim==1&&dx>-150&&dx<150");
93c2->cd(4);
94gPad->SetFillColor(33);
95 ntuple->SetFillColor(46);
96 ntuple->Draw("dz","lay==2&&hitprim==1&&dz>-500&&dz<500");
97 */
e2ce6aca 98
e8189707 99/*
100c2->cd(1);
101gPad->SetFillColor(33);
102 ntuple->SetFillColor(42);
103 ntuple->Draw("dx","lay==1&&hitprim==1&&nx<4");
104c2->cd(2);
105gPad->SetFillColor(33);
106 ntuple->SetFillColor(46);
107 ntuple->Draw("dz","lay==1&&hitprim==1&&nz<4");
108c2->cd(3);
109gPad->SetFillColor(33);
110 ntuple->SetFillColor(42);
111 ntuple->Draw("dx","lay==2&&hitprim==1&&nx<4");
112c2->cd(4);
113gPad->SetFillColor(33);
114 ntuple->SetFillColor(46);
115 ntuple->Draw("dz","lay==2&&hitprim==1&&nz<4");
116*/
117
e2ce6aca 118
e8189707 119c2->cd(1);
120gPad->SetFillColor(33);
121 ntuple1->SetFillColor(42);
122 ntuple1->Draw("noverlaps","lay == 1");
123c2->cd(2);
124gPad->SetFillColor(33);
125 ntuple1->SetFillColor(46);
126 ntuple1->Draw("noverlaps","lay == 2");
e8189707 127c2->cd(3);
128gPad->SetFillColor(33);
129 ntuple1->SetFillColor(42);
e2ce6aca 130 ntuple1->Draw("ntrover","lay == 1");
e8189707 131c2->cd(4);
132gPad->SetFillColor(33);
133 ntuple1->SetFillColor(46);
e2ce6aca 134 ntuple1->Draw("ntrover","lay == 2");
135
e8189707 136/*
137c2->cd(1);
138gPad->SetFillColor(33);
139 ntuple1->SetFillColor(42);
140 ntuple1->Draw("noverprim","lay == 1");
141c2->cd(2);
142gPad->SetFillColor(33);
143 ntuple1->SetFillColor(46);
144 ntuple1->Draw("noverprim","lay == 2");
145c2->cd(3);
146gPad->SetFillColor(33);
147 ntuple1->SetFillColor(42);
148 ntuple1->Draw("noverprim","lay == 1 && noverprim > 1");
149c2->cd(4);
150gPad->SetFillColor(33);
151 ntuple1->SetFillColor(46);
152 ntuple1->Draw("noverprim","lay == 2 && noverprim > 1");
153*/
e2ce6aca 154
e8189707 155/*
156c2->cd(1);
157gPad->SetFillColor(33);
158 ntuple1->SetFillColor(42);
159 ntuple1->Draw("qcl","lay == 1&&noverprim>0");
160c2->cd(2);
161gPad->SetFillColor(33);
162 ntuple1->SetFillColor(46);
163 ntuple1->Draw("qcl","lay == 2&&noverprim>0");
164*/
165
e2ce6aca 166/*
e8189707 167c2->cd(1);
168gPad->SetFillColor(33);
169 ntuple1->SetFillColor(42);
170 ntuple1->Draw("nx","lay == 1&&noverprim>0");
171c2->cd(2);
172gPad->SetFillColor(33);
173 ntuple1->SetFillColor(46);
174 ntuple1->Draw("nx","lay == 2&&noverprim>0");
175c2->cd(3);
176gPad->SetFillColor(33);
177 ntuple1->SetFillColor(42);
178 ntuple1->Draw("nz","lay == 1 && noverprim>0");
179c2->cd(4);
180gPad->SetFillColor(33);
181 ntuple1->SetFillColor(46);
182 ntuple1->Draw("nz","lay == 2 && noverprim>0");
e2ce6aca 183*/
e8189707 184
185/*
186c2->cd(1);
187gPad->SetFillColor(33);
188 ntuple1->SetFillColor(42);
189 ntuple1->Draw("dx","lay == 1&&noverprim>0");
190c2->cd(2);
191gPad->SetFillColor(33);
192 ntuple1->SetFillColor(46);
193 ntuple1->Draw("dx","lay == 2&&noverprim>0");
194c2->cd(3);
195gPad->SetFillColor(33);
196 ntuple1->SetFillColor(42);
197 ntuple1->Draw("dz","lay == 1 && noverprim>0");
198c2->cd(4);
199gPad->SetFillColor(33);
200 ntuple1->SetFillColor(46);
201 ntuple1->Draw("dz","lay == 2 && noverprim>0");
202*/
203
204 /*
205c2->cd(1);
206gPad->SetFillColor(33);
207 ntuple1->SetFillColor(42);
208 ntuple1->Draw("dx","lay == 1");
209c2->cd(2);
210gPad->SetFillColor(33);
211 ntuple1->SetFillColor(46);
212 ntuple1->Draw("dx","lay == 2");
213c2->cd(3);
214gPad->SetFillColor(33);
215 ntuple1->SetFillColor(42);
216 ntuple1->Draw("dz","lay == 1");
217c2->cd(4);
218gPad->SetFillColor(33);
219 ntuple1->SetFillColor(46);
220 ntuple1->Draw("dz","lay == 2");
221 */
e2ce6aca 222
e8189707 223 /*
224c2->cd(1);
225gPad->SetFillColor(33);
226 ntuple2->SetFillColor(42);
227 ntuple2->Draw("nx","lay == 1");
228c2->cd(2);
229gPad->SetFillColor(33);
230 ntuple2->SetFillColor(46);
231 ntuple2->Draw("nz","lay == 1");
232c2->cd(3);
233gPad->SetFillColor(33);
234 ntuple2->SetFillColor(42);
235 ntuple2->Draw("nx","lay == 2");
236c2->cd(4);
237gPad->SetFillColor(33);
238 ntuple2->SetFillColor(46);
239 ntuple2->Draw("nz","lay == 2");
240 */
241
242///////////////////// Histogramm analysis ////////////////////////
243
244 /*
245c2->cd(1);
246gPad->SetFillColor(33);
247 Xres1->SetFillColor(42);
248 Xres1->Draw();
249c2->cd(2);
250gPad->SetFillColor(33);
251 Zres1->SetFillColor(42);
252 Zres1->Draw();
253c2->cd(3);
254gPad->SetFillColor(33);
255 Xres2->SetFillColor(46);
256 Xres2->Draw();
257c2->cd(4);
258gPad->SetFillColor(33);
259 Zres2->SetFillColor(46);
260 Zres2->Draw();
261 */
262 /*
263c2->cd(1);
264gPad->SetFillColor(33);
265 Nzpix1->SetFillColor(42);
266 Nzpix1->Draw();
267c2->cd(2);
268gPad->SetFillColor(33);
269 Nxpix1->SetFillColor(42);
270 Nxpix1->Draw();
271c2->cd(3);
272gPad->SetFillColor(33);
273 Nzpix2->SetFillColor(46);
274 Nzpix2->Draw();
275c2->cd(4);
276gPad->SetFillColor(33);
277 Nxpix2->SetFillColor(46);
278 Nxpix2->Draw();
279 */
280
281 /*
282c2->cd(1);
283gPad->SetFillColor(33);
284 Zpix1->SetFillColor(42);
285 Zpix1->Draw();
286c2->cd(2);
287gPad->SetFillColor(33);
288 Xpix1->SetFillColor(46);
289 Xpix1->Draw();
290c2->cd(3);
291gPad->SetFillColor(33);
292 Zpix2->SetFillColor(42);
293 Zpix2->Draw();
294c2->cd(4);
295gPad->SetFillColor(33);
296 Xpix2->SetFillColor(46);
297 Xpix2->Draw();
298 */
299
300 /*
301c2->Draw();
302c2->Print("spd_res.ps");
303 */
304 /*
305c2->Draw();
306c2->Print("spd_clsize.ps");
307 */
308
309}
310
311
312
313
314
315
316
317