4 TFile *f = new TFile("SSD_his.root");
6 //gStyle->SetOptStat(1111111);
7 //gStyle->SetOptLogy();
8 //TCanvas *c1 = new TCanvas("c1","SPD clusters",400,10,600,700);
9 TCanvas *c2 = new TCanvas("c2","SSD clusters",400,10,600,700);
13 ///////////////////////// Ntuple analysis ///////////////////////////////
15 // ntuple is created inside the hit loop for the hits corresponding to the
18 // ntuple1 is created after a finish of the hit loop if one or more hits
19 // correspond to the recpoint;
21 // ntuple2 is created befor the hit loop for all recpoints;
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 // -------------------------------------------------------------------------
46 gPad->SetFillColor(33);
47 ntuple->SetFillColor(42);
48 ntuple->Draw("dx","lay == 5 && hitprim == 1");
50 gPad->SetFillColor(33);
51 ntuple->SetFillColor(46);
52 ntuple->Draw("dz","lay == 5 && hitprim == 1");
54 gPad->SetFillColor(33);
55 ntuple->SetFillColor(42);
56 ntuple->Draw("dx","lay == 6 && hitprim == 1");
58 gPad->SetFillColor(33);
59 ntuple->SetFillColor(46);
60 ntuple->Draw("dz","lay == 6 && hitprim == 1");
65 gPad->SetFillColor(33);
66 ntuple->SetFillColor(42);
67 ntuple->Draw("dx","lay == 5 && hitprim == 1&&ntrover==1");
69 gPad->SetFillColor(33);
70 ntuple->SetFillColor(46);
71 ntuple->Draw("dz","lay == 5 && hitprim == 1&&ntrover==1");
73 gPad->SetFillColor(33);
74 ntuple->SetFillColor(42);
75 ntuple->Draw("dx","lay == 6 && hitprim == 1&&ntrover==1");
77 gPad->SetFillColor(33);
78 ntuple->SetFillColor(46);
79 ntuple->Draw("dz","lay == 6 && hitprim == 1&&ntrover==1");
85 gPad->SetFillColor(33);
86 ntuple->SetFillColor(42);
87 ntuple->Draw("dx","lay == 5 && hitprim == 1&&abs(dx)<200");
89 gPad->SetFillColor(33);
90 ntuple->SetFillColor(46);
91 ntuple->Draw("dz","lay == 5 && hitprim == 1&&abs(dz)<5000");
93 gPad->SetFillColor(33);
94 ntuple->SetFillColor(42);
95 ntuple->Draw("dx","lay == 6 && hitprim == 1&&abs(dx)<200");
97 gPad->SetFillColor(33);
98 ntuple->SetFillColor(46);
99 ntuple->Draw("dz","lay == 6 && hitprim == 1&&abs(dz)<5000");
104 gPad->SetFillColor(33);
105 ntuple->SetFillColor(42);
106 ntuple->Draw("dx","lay == 5 && hitprim == 1&&abs(dx)<200&&ntrover==1");
108 gPad->SetFillColor(33);
109 ntuple->SetFillColor(46);
110 ntuple->Draw("dz","lay == 5 && hitprim == 1&&abs(dz)<5000&&ntrover==1");
112 gPad->SetFillColor(33);
113 ntuple->SetFillColor(42);
114 ntuple->Draw("dx","lay == 6 && hitprim == 1&&abs(dx)<200&&ntrover==1");
116 gPad->SetFillColor(33);
117 ntuple->SetFillColor(46);
118 ntuple->Draw("dz","lay == 6 && hitprim == 1&&abs(dz)<5000&&ntrover==1");
123 gPad->SetFillColor(33);
124 ntuple1->SetFillColor(42);
125 ntuple1->Draw("nxP","lay == 5&&noverprim>=0");
127 gPad->SetFillColor(33);
128 ntuple1->SetFillColor(46);
129 ntuple1->Draw("nxN","lay == 5&&noverprim>=0");
131 gPad->SetFillColor(33);
132 ntuple1->SetFillColor(42);
133 ntuple1->Draw("nxP","lay == 6 && noverprim>=0");
135 gPad->SetFillColor(33);
136 ntuple1->SetFillColor(46);
137 ntuple1->Draw("nxN","lay == 6 && noverprim>=0");
142 gPad->SetFillColor(33);
143 ntuple1->SetFillColor(42);
144 ntuple1->Draw("dx","lay == 5&&noverprim>=0&&abs(dx)<50&&nxP==2&&nxN==2");
146 gPad->SetFillColor(33);
147 ntuple1->SetFillColor(46);
148 ntuple1->Draw("dz","lay == 5&&noverprim>=0&&abs(dz)<2000&&nxP==2&&nxN==2");
150 gPad->SetFillColor(33);
151 ntuple1->SetFillColor(42);
152 ntuple1->Draw("dx","lay == 6 && noverprim>=0&&abs(dx)<200");
154 gPad->SetFillColor(33);
155 ntuple1->SetFillColor(46);
156 ntuple1->Draw("dz","lay == 6 && noverprim>=0&&abs(dz)<5000");
161 gPad->SetFillColor(33);
162 ntuple1->SetFillColor(42);
163 ntuple1->Draw("dx","lay == 5&&noverprim>=0&&noverlaps==0&&qcut<0.18");
165 gPad->SetFillColor(33);
166 ntuple1->SetFillColor(46);
167 ntuple1->Draw("dz","lay == 5&&noverprim>=0&&noverlaps==0&&qcut<0.18");
169 gPad->SetFillColor(33);
170 ntuple1->SetFillColor(42);
171 ntuple1->Draw("dx","lay == 6 && noverprim>=0&&noverlaps==0&&qcut<0.18");
173 gPad->SetFillColor(33);
174 ntuple1->SetFillColor(46);
175 ntuple1->Draw("dz","lay == 6 && noverprim>=0&&noverlaps==0&&qcut<0.18");
180 gPad->SetFillColor(33);
181 ntuple1->SetFillColor(42);
182 ntuple1->Draw("dx","lay == 5&&noverprim>=0&&noverlaps==0&&abs(dx)<100&&qcut<0.18");
184 gPad->SetFillColor(33);
185 ntuple1->SetFillColor(46);
186 ntuple1->Draw("dz","lay == 5&&noverprim>=0&&noverlaps==0&&abs(dz)<5000&&qcut<0.18");
188 gPad->SetFillColor(33);
189 ntuple1->SetFillColor(42);
190 ntuple1->Draw("dx","lay == 6&&noverprim>=0&&noverlaps==0&&abs(dx)<100&&qcut<0.18");
192 gPad->SetFillColor(33);
193 ntuple1->SetFillColor(46);
194 ntuple1->Draw("dz","lay == 6&&noverprim>=0&&noverlaps==0&&abs(dz)<5000&&qcut<0.18");
199 gPad->SetFillColor(33);
200 ntuple1->SetFillColor(42);
201 ntuple1->Draw("noverprim","lay == 5&&noverprim>=-1");
203 gPad->SetFillColor(33);
204 ntuple1->SetFillColor(46);
205 ntuple1->Draw("noverlaps","lay == 5&&noverprim>=-1");
207 gPad->SetFillColor(33);
208 ntuple1->SetFillColor(42);
209 ntuple1->Draw("noverprim","lay == 6 && noverprim>=-1");
211 gPad->SetFillColor(33);
212 ntuple1->SetFillColor(46);
213 ntuple1->Draw("noverlaps","lay == 6 && noverprim>=-1");
218 gPad->SetFillColor(33);
219 ntuple1->SetFillColor(42);
220 ntuple1->Draw("ntrover","lay == 5&&noverprim>=0");
222 gPad->SetFillColor(33);
223 ntuple1->SetFillColor(46);
224 ntuple1->Draw("noverlaps","lay == 5&&noverprim>=0");
226 gPad->SetFillColor(33);
227 ntuple1->SetFillColor(42);
228 ntuple1->Draw("ntrover","lay == 6 && noverprim>=0");
230 gPad->SetFillColor(33);
231 ntuple1->SetFillColor(46);
232 ntuple1->Draw("noverlaps","lay == 6 && noverprim>=0");
237 gPad->SetFillColor(33);
238 ntuple1->SetFillColor(42);
239 ntuple1->Draw("qclP","lay == 5&&noverprim>=0");
241 gPad->SetFillColor(33);
242 ntuple1->SetFillColor(42);
243 ntuple1->Draw("qclN","lay == 5&&noverprim>=0");
245 gPad->SetFillColor(33);
246 ntuple1->SetFillColor(46);
247 ntuple1->Draw("qclP","lay == 6&&noverprim>=0");
249 gPad->SetFillColor(33);
250 ntuple1->SetFillColor(46);
251 ntuple1->Draw("qclN","lay == 6&&noverprim>=0");
255 ///////////////////// Histogramm/ntuple analysis ////////////////////////
259 gPad->SetFillColor(33);
260 adcPadcN5all->SetFillColor(42);
261 adcPadcN5all->Draw();
263 gPad->SetFillColor(33);
264 adcPadcN6all->SetFillColor(42);
265 adcPadcN6all->Draw();
267 gPad->SetFillColor(33);
268 adcPadcN5cut->SetFillColor(42);
269 adcPadcN5cut->Draw();
271 gPad->SetFillColor(33);
272 adcPadcN6cut->SetFillColor(46);
273 adcPadcN6cut->Draw();
278 gPad->SetFillColor(33);
279 ntuple1->SetFillColor(42);
280 ntuple1->Draw("qclP","noverprim>=0&&qclP<500");
282 gPad->SetFillColor(33);
283 ntuple1->SetFillColor(42);
284 ntuple1->Draw("qclN","noverprim>=0&&qclN<500");
286 gPad->SetFillColor(33);
287 adcPadcN5cut->SetFillColor(42);
288 adcPadcN5cut->Draw();
290 gPad->SetFillColor(33);
291 adcPadcN6cut->SetFillColor(46);
292 adcPadcN6cut->Draw();
298 c2->Print("ssd_res.ps");
303 c2->Print("spd_clsize.ps");