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