]>
Commit | Line | Data |
---|---|---|
6b8f55ce | 1 | void SSD_ntuple() |
2 | { | |
3 | ||
4 | TFile *f = new TFile("SSD_his.root"); | |
5 | ||
6 | //gStyle->SetOptStat(1111111); | |
7 | //gStyle->SetOptLogy(); | |
8 | //TCanvas *c1 = new TCanvas("c1","SPD clusters",400,10,600,700); | |
6298ec22 | 9 | TCanvas *c2 = new TCanvas("c2","SSD clusters",400,10,600,700); |
6b8f55ce | 10 | //c1->Divide(2,2); |
11 | c2->Divide(2,2); | |
12 | ||
13 | ///////////////////////// Ntuple analysis /////////////////////////////// | |
14 | ||
15 | // ntuple is created inside the hit loop for the hits corresponding to the | |
16 | // recpoint; | |
17 | ||
18 | // ntuple1 is created after a finish of the hit loop if one or more hits | |
19 | // correspond to the recpoint; | |
20 | ||
21 | // ntuple2 is created befor the hit loop for all recpoints; | |
22 | ||
23 | // ----------------------------------------------------------------------- | |
24 | // lay - number of ITS layer; | |
25 | // lad - number of ITS ladder; | |
26 | // det - number of ITS detector; | |
27 | // nxP/N - cluster size in the r*phi(x) direction for P/N sides; | |
28 | // hitprim - primary particle(hit) flag ( = 1 for primery particle); | |
29 | // x - x local coordinate in mm; | |
30 | // z - z local coordinate in mm; | |
31 | // dx - difference of hit(mediate) and reconstructed (from cluster) | |
32 | // coordinates in r*phi(x) direction in microns; | |
33 | // dz - difference of hit(mediate) and reconstructed (from cluster) | |
34 | // coordinates in z direction in microns; | |
35 | // noverlaps - number of particles overlapping in one cluster; | |
36 | // noverprim - number of primary particles overlapping in one cluster; | |
37 | // qclP/N - cluster signals in ADC normalized to the path length in Si | |
38 | // qrec - recpoint signal (maximum from qclP and qclN) | |
39 | // pmod - particle momentum at the vertex in MeV/c | |
40 | // partcode - particle code | |
41 | // ------------------------------------------------------------------------- | |
42 | ||
43 | ||
6298ec22 | 44 | /* |
45 | c2->cd(1); | |
46 | gPad->SetFillColor(33); | |
47 | ntuple->SetFillColor(42); | |
48 | ntuple->Draw("dx","lay == 5 && hitprim == 1"); | |
49 | c2->cd(2); | |
50 | gPad->SetFillColor(33); | |
51 | ntuple->SetFillColor(46); | |
52 | ntuple->Draw("dz","lay == 5 && hitprim == 1"); | |
53 | c2->cd(3); | |
54 | gPad->SetFillColor(33); | |
55 | ntuple->SetFillColor(42); | |
56 | ntuple->Draw("dx","lay == 6 && hitprim == 1"); | |
57 | c2->cd(4); | |
58 | gPad->SetFillColor(33); | |
59 | ntuple->SetFillColor(46); | |
60 | ntuple->Draw("dz","lay == 6 && hitprim == 1"); | |
61 | */ | |
6b8f55ce | 62 | |
6298ec22 | 63 | /* |
64 | c2->cd(1); | |
65 | gPad->SetFillColor(33); | |
66 | ntuple->SetFillColor(42); | |
67 | ntuple->Draw("dx","lay == 5 && hitprim == 1&&ntrover==1"); | |
68 | c2->cd(2); | |
69 | gPad->SetFillColor(33); | |
70 | ntuple->SetFillColor(46); | |
71 | ntuple->Draw("dz","lay == 5 && hitprim == 1&&ntrover==1"); | |
72 | c2->cd(3); | |
73 | gPad->SetFillColor(33); | |
74 | ntuple->SetFillColor(42); | |
75 | ntuple->Draw("dx","lay == 6 && hitprim == 1&&ntrover==1"); | |
76 | c2->cd(4); | |
77 | gPad->SetFillColor(33); | |
78 | ntuple->SetFillColor(46); | |
79 | ntuple->Draw("dz","lay == 6 && hitprim == 1&&ntrover==1"); | |
80 | */ | |
81 | ||
82 | ||
83 | /* | |
6b8f55ce | 84 | c2->cd(1); |
85 | gPad->SetFillColor(33); | |
86 | ntuple->SetFillColor(42); | |
87 | ntuple->Draw("dx","lay == 5 && hitprim == 1&&abs(dx)<200"); | |
88 | c2->cd(2); | |
89 | gPad->SetFillColor(33); | |
90 | ntuple->SetFillColor(46); | |
91 | ntuple->Draw("dz","lay == 5 && hitprim == 1&&abs(dz)<5000"); | |
92 | c2->cd(3); | |
93 | gPad->SetFillColor(33); | |
94 | ntuple->SetFillColor(42); | |
95 | ntuple->Draw("dx","lay == 6 && hitprim == 1&&abs(dx)<200"); | |
96 | c2->cd(4); | |
97 | gPad->SetFillColor(33); | |
98 | ntuple->SetFillColor(46); | |
99 | ntuple->Draw("dz","lay == 6 && hitprim == 1&&abs(dz)<5000"); | |
6298ec22 | 100 | */ |
6b8f55ce | 101 | |
6298ec22 | 102 | /* |
103 | c2->cd(1); | |
104 | gPad->SetFillColor(33); | |
105 | ntuple->SetFillColor(42); | |
106 | ntuple->Draw("dx","lay == 5 && hitprim == 1&&abs(dx)<200&&ntrover==1"); | |
107 | c2->cd(2); | |
108 | gPad->SetFillColor(33); | |
109 | ntuple->SetFillColor(46); | |
110 | ntuple->Draw("dz","lay == 5 && hitprim == 1&&abs(dz)<5000&&ntrover==1"); | |
111 | c2->cd(3); | |
112 | gPad->SetFillColor(33); | |
113 | ntuple->SetFillColor(42); | |
114 | ntuple->Draw("dx","lay == 6 && hitprim == 1&&abs(dx)<200&&ntrover==1"); | |
115 | c2->cd(4); | |
116 | gPad->SetFillColor(33); | |
117 | ntuple->SetFillColor(46); | |
118 | ntuple->Draw("dz","lay == 6 && hitprim == 1&&abs(dz)<5000&&ntrover==1"); | |
119 | */ | |
6b8f55ce | 120 | |
0315d466 | 121 | /* |
6b8f55ce | 122 | c2->cd(1); |
123 | gPad->SetFillColor(33); | |
124 | ntuple1->SetFillColor(42); | |
125 | ntuple1->Draw("nxP","lay == 5&&noverprim>=0"); | |
126 | c2->cd(2); | |
127 | gPad->SetFillColor(33); | |
128 | ntuple1->SetFillColor(46); | |
129 | ntuple1->Draw("nxN","lay == 5&&noverprim>=0"); | |
130 | c2->cd(3); | |
131 | gPad->SetFillColor(33); | |
132 | ntuple1->SetFillColor(42); | |
133 | ntuple1->Draw("nxP","lay == 6 && noverprim>=0"); | |
134 | c2->cd(4); | |
135 | gPad->SetFillColor(33); | |
136 | ntuple1->SetFillColor(46); | |
137 | ntuple1->Draw("nxN","lay == 6 && noverprim>=0"); | |
0315d466 | 138 | */ |
6b8f55ce | 139 | |
6298ec22 | 140 | /* |
141 | c2->cd(1); | |
142 | gPad->SetFillColor(33); | |
143 | ntuple1->SetFillColor(42); | |
144 | ntuple1->Draw("dx","lay == 5&&noverprim>=0&&abs(dx)<50&&nxP==2&&nxN==2"); | |
145 | c2->cd(2); | |
146 | gPad->SetFillColor(33); | |
147 | ntuple1->SetFillColor(46); | |
148 | ntuple1->Draw("dz","lay == 5&&noverprim>=0&&abs(dz)<2000&&nxP==2&&nxN==2"); | |
149 | c2->cd(3); | |
150 | gPad->SetFillColor(33); | |
151 | ntuple1->SetFillColor(42); | |
152 | ntuple1->Draw("dx","lay == 6 && noverprim>=0&&abs(dx)<200"); | |
153 | c2->cd(4); | |
154 | gPad->SetFillColor(33); | |
155 | ntuple1->SetFillColor(46); | |
156 | ntuple1->Draw("dz","lay == 6 && noverprim>=0&&abs(dz)<5000"); | |
157 | */ | |
158 | ||
159 | /* | |
160 | c2->cd(1); | |
161 | gPad->SetFillColor(33); | |
162 | ntuple1->SetFillColor(42); | |
163 | ntuple1->Draw("dx","lay == 5&&noverprim>=0&&noverlaps==0&&qcut<0.18"); | |
164 | c2->cd(2); | |
165 | gPad->SetFillColor(33); | |
166 | ntuple1->SetFillColor(46); | |
167 | ntuple1->Draw("dz","lay == 5&&noverprim>=0&&noverlaps==0&&qcut<0.18"); | |
168 | c2->cd(3); | |
169 | gPad->SetFillColor(33); | |
170 | ntuple1->SetFillColor(42); | |
171 | ntuple1->Draw("dx","lay == 6 && noverprim>=0&&noverlaps==0&&qcut<0.18"); | |
172 | c2->cd(4); | |
173 | gPad->SetFillColor(33); | |
174 | ntuple1->SetFillColor(46); | |
175 | ntuple1->Draw("dz","lay == 6 && noverprim>=0&&noverlaps==0&&qcut<0.18"); | |
176 | */ | |
177 | ||
178 | /* | |
179 | c2->cd(1); | |
180 | gPad->SetFillColor(33); | |
181 | ntuple1->SetFillColor(42); | |
182 | ntuple1->Draw("dx","lay == 5&&noverprim>=0&&noverlaps==0&&abs(dx)<100&&qcut<0.18"); | |
183 | c2->cd(2); | |
184 | gPad->SetFillColor(33); | |
185 | ntuple1->SetFillColor(46); | |
186 | ntuple1->Draw("dz","lay == 5&&noverprim>=0&&noverlaps==0&&abs(dz)<5000&&qcut<0.18"); | |
187 | c2->cd(3); | |
188 | gPad->SetFillColor(33); | |
189 | ntuple1->SetFillColor(42); | |
190 | ntuple1->Draw("dx","lay == 6&&noverprim>=0&&noverlaps==0&&abs(dx)<100&&qcut<0.18"); | |
191 | c2->cd(4); | |
192 | gPad->SetFillColor(33); | |
193 | ntuple1->SetFillColor(46); | |
194 | ntuple1->Draw("dz","lay == 6&&noverprim>=0&&noverlaps==0&&abs(dz)<5000&&qcut<0.18"); | |
195 | */ | |
196 | ||
197 | /* | |
198 | c2->cd(1); | |
199 | gPad->SetFillColor(33); | |
200 | ntuple1->SetFillColor(42); | |
201 | ntuple1->Draw("noverprim","lay == 5&&noverprim>=-1"); | |
202 | c2->cd(2); | |
203 | gPad->SetFillColor(33); | |
204 | ntuple1->SetFillColor(46); | |
205 | ntuple1->Draw("noverlaps","lay == 5&&noverprim>=-1"); | |
206 | c2->cd(3); | |
207 | gPad->SetFillColor(33); | |
208 | ntuple1->SetFillColor(42); | |
209 | ntuple1->Draw("noverprim","lay == 6 && noverprim>=-1"); | |
210 | c2->cd(4); | |
211 | gPad->SetFillColor(33); | |
212 | ntuple1->SetFillColor(46); | |
213 | ntuple1->Draw("noverlaps","lay == 6 && noverprim>=-1"); | |
214 | */ | |
215 | ||
216 | /* | |
217 | c2->cd(1); | |
218 | gPad->SetFillColor(33); | |
219 | ntuple1->SetFillColor(42); | |
220 | ntuple1->Draw("ntrover","lay == 5&&noverprim>=0"); | |
221 | c2->cd(2); | |
222 | gPad->SetFillColor(33); | |
223 | ntuple1->SetFillColor(46); | |
224 | ntuple1->Draw("noverlaps","lay == 5&&noverprim>=0"); | |
225 | c2->cd(3); | |
226 | gPad->SetFillColor(33); | |
227 | ntuple1->SetFillColor(42); | |
228 | ntuple1->Draw("ntrover","lay == 6 && noverprim>=0"); | |
229 | c2->cd(4); | |
230 | gPad->SetFillColor(33); | |
231 | ntuple1->SetFillColor(46); | |
232 | ntuple1->Draw("noverlaps","lay == 6 && noverprim>=0"); | |
233 | */ | |
234 | ||
6b8f55ce | 235 | /* |
236 | c2->cd(1); | |
237 | gPad->SetFillColor(33); | |
238 | ntuple1->SetFillColor(42); | |
239 | ntuple1->Draw("qclP","lay == 5&&noverprim>=0"); | |
240 | c2->cd(2); | |
241 | gPad->SetFillColor(33); | |
242 | ntuple1->SetFillColor(42); | |
243 | ntuple1->Draw("qclN","lay == 5&&noverprim>=0"); | |
244 | c2->cd(3); | |
245 | gPad->SetFillColor(33); | |
246 | ntuple1->SetFillColor(46); | |
247 | ntuple1->Draw("qclP","lay == 6&&noverprim>=0"); | |
248 | c2->cd(4); | |
249 | gPad->SetFillColor(33); | |
250 | ntuple1->SetFillColor(46); | |
251 | ntuple1->Draw("qclN","lay == 6&&noverprim>=0"); | |
252 | */ | |
253 | ||
254 | ||
255 | ///////////////////// Histogramm/ntuple analysis //////////////////////// | |
256 | ||
6298ec22 | 257 | |
258 | c2->cd(1); | |
259 | gPad->SetFillColor(33); | |
260 | adcPadcN5all->SetFillColor(42); | |
261 | adcPadcN5all->Draw(); | |
262 | c2->cd(2); | |
263 | gPad->SetFillColor(33); | |
264 | adcPadcN6all->SetFillColor(42); | |
265 | adcPadcN6all->Draw(); | |
266 | c2->cd(3); | |
267 | gPad->SetFillColor(33); | |
268 | adcPadcN5cut->SetFillColor(42); | |
269 | adcPadcN5cut->Draw(); | |
270 | c2->cd(4); | |
271 | gPad->SetFillColor(33); | |
272 | adcPadcN6cut->SetFillColor(46); | |
273 | adcPadcN6cut->Draw(); | |
274 | ||
275 | ||
6b8f55ce | 276 | /* |
277 | c2->cd(1); | |
278 | gPad->SetFillColor(33); | |
0315d466 | 279 | ntuple1->SetFillColor(42); |
6298ec22 | 280 | ntuple1->Draw("qclP","noverprim>=0&&qclP<500"); |
6b8f55ce | 281 | c2->cd(2); |
282 | gPad->SetFillColor(33); | |
0315d466 | 283 | ntuple1->SetFillColor(42); |
6298ec22 | 284 | ntuple1->Draw("qclN","noverprim>=0&&qclN<500"); |
6b8f55ce | 285 | c2->cd(3); |
286 | gPad->SetFillColor(33); | |
287 | adcPadcN5cut->SetFillColor(42); | |
288 | adcPadcN5cut->Draw(); | |
289 | c2->cd(4); | |
290 | gPad->SetFillColor(33); | |
291 | adcPadcN6cut->SetFillColor(46); | |
292 | adcPadcN6cut->Draw(); | |
293 | */ | |
294 | ||
295 | ||
296 | /* | |
297 | c2->Draw(); | |
298 | c2->Print("ssd_res.ps"); | |
299 | */ | |
300 | ||
301 | /* | |
302 | c2->Draw(); | |
303 | c2->Print("spd_clsize.ps"); | |
304 | */ | |
305 | ||
306 | } | |
307 | ||
308 | ||
309 | ||
310 | ||
311 | ||
312 | ||
313 | ||
314 | ||
315 |