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