4 TFile *f = new TFile("SPD_his.root");
5 //TH1F *Pix = (TH1F*)f->Get("Pix");
7 //gStyle->SetOptStat(1111111);
9 //TCanvas *c1 = new TCanvas("c1","SPD clusters",400,10,600,700);
10 TCanvas *c2 = new TCanvas("c2","SPD clusters",400,10,600,700);
15 ///////////////////////// Ntuple analysis ///////////////////////////////
17 // ntuple is created inside the hit loop for the hits in the cluster region;
19 // ntuple1 is created after a finish of the hit loop if one or more hits
20 // are in the cluster region;
22 // ntuple2 is created befor the hit loop for all clusters;
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 finded in this
34 // macros in the cluster region;
35 // ntrover - number of particles overlapping in one cluster finded in the
36 // AliITSClusterFinderSPD class;
37 // noverprim - number of primary particles overlapping in one cluster;
38 // qcl - cluster charge in electron number
39 // -------------------------------------------------------------------------
43 gPad->SetFillColor(33);
44 ntuple->SetFillColor(42);
45 ntuple->Draw("dx","lay == 1 && hitprim == 1");
47 gPad->SetFillColor(33);
48 ntuple->SetFillColor(46);
49 ntuple->Draw("dz","lay == 1 && hitprim == 1");
51 gPad->SetFillColor(33);
52 ntuple->SetFillColor(42);
53 ntuple->Draw("dx","lay == 2 && hitprim == 1");
55 gPad->SetFillColor(33);
56 ntuple->SetFillColor(46);
57 ntuple->Draw("dz","lay == 2 && hitprim == 1");
62 gPad->SetFillColor(33);
63 ntuple->SetFillColor(42);
64 ntuple->Draw("dx","lay == 1 && hitprim == 1");
66 gPad->SetFillColor(33);
67 ntuple->SetFillColor(46);
68 ntuple->Draw("dx","lay == 1 && hitprim == 1 && nx < 15");
70 gPad->SetFillColor(33);
71 ntuple->SetFillColor(42);
72 ntuple->Draw("dx","lay == 1 && hitprim == 1 && nx < 5");
74 gPad->SetFillColor(33);
75 ntuple->SetFillColor(46);
76 ntuple->Draw("dx","lay == 1 && hitprim == 1 && nx < 4");
82 gPad->SetFillColor(33);
83 ntuple->SetFillColor(42);
84 ntuple->Draw("dx","lay==1&&hitprim==1&&dx>-150&&dx<150");
86 gPad->SetFillColor(33);
87 ntuple->SetFillColor(46);
88 ntuple->Draw("dz","lay==1&&hitprim==1&&dz>-500&&dz<500");
90 gPad->SetFillColor(33);
91 ntuple->SetFillColor(42);
92 ntuple->Draw("dx","lay==2&&hitprim==1&&dx>-150&&dx<150");
94 gPad->SetFillColor(33);
95 ntuple->SetFillColor(46);
96 ntuple->Draw("dz","lay==2&&hitprim==1&&dz>-500&&dz<500");
101 gPad->SetFillColor(33);
102 ntuple->SetFillColor(42);
103 ntuple->Draw("dx","lay==1&&hitprim==1&&nx<4");
105 gPad->SetFillColor(33);
106 ntuple->SetFillColor(46);
107 ntuple->Draw("dz","lay==1&&hitprim==1&&nz<4");
109 gPad->SetFillColor(33);
110 ntuple->SetFillColor(42);
111 ntuple->Draw("dx","lay==2&&hitprim==1&&nx<4");
113 gPad->SetFillColor(33);
114 ntuple->SetFillColor(46);
115 ntuple->Draw("dz","lay==2&&hitprim==1&&nz<4");
120 gPad->SetFillColor(33);
121 ntuple1->SetFillColor(42);
122 ntuple1->Draw("noverlaps","lay == 1");
124 gPad->SetFillColor(33);
125 ntuple1->SetFillColor(46);
126 ntuple1->Draw("noverlaps","lay == 2");
128 gPad->SetFillColor(33);
129 ntuple1->SetFillColor(42);
130 ntuple1->Draw("ntrover","lay == 1");
132 gPad->SetFillColor(33);
133 ntuple1->SetFillColor(46);
134 ntuple1->Draw("ntrover","lay == 2");
138 gPad->SetFillColor(33);
139 ntuple1->SetFillColor(42);
140 ntuple1->Draw("noverprim","lay == 1");
142 gPad->SetFillColor(33);
143 ntuple1->SetFillColor(46);
144 ntuple1->Draw("noverprim","lay == 2");
146 gPad->SetFillColor(33);
147 ntuple1->SetFillColor(42);
148 ntuple1->Draw("noverprim","lay == 1 && noverprim > 1");
150 gPad->SetFillColor(33);
151 ntuple1->SetFillColor(46);
152 ntuple1->Draw("noverprim","lay == 2 && noverprim > 1");
157 gPad->SetFillColor(33);
158 ntuple1->SetFillColor(42);
159 ntuple1->Draw("qcl","lay == 1&&noverprim>0");
161 gPad->SetFillColor(33);
162 ntuple1->SetFillColor(46);
163 ntuple1->Draw("qcl","lay == 2&&noverprim>0");
168 gPad->SetFillColor(33);
169 ntuple1->SetFillColor(42);
170 ntuple1->Draw("nx","lay == 1&&noverprim>0");
172 gPad->SetFillColor(33);
173 ntuple1->SetFillColor(46);
174 ntuple1->Draw("nx","lay == 2&&noverprim>0");
176 gPad->SetFillColor(33);
177 ntuple1->SetFillColor(42);
178 ntuple1->Draw("nz","lay == 1 && noverprim>0");
180 gPad->SetFillColor(33);
181 ntuple1->SetFillColor(46);
182 ntuple1->Draw("nz","lay == 2 && noverprim>0");
187 gPad->SetFillColor(33);
188 ntuple1->SetFillColor(42);
189 ntuple1->Draw("dx","lay == 1&&noverprim>0");
191 gPad->SetFillColor(33);
192 ntuple1->SetFillColor(46);
193 ntuple1->Draw("dx","lay == 2&&noverprim>0");
195 gPad->SetFillColor(33);
196 ntuple1->SetFillColor(42);
197 ntuple1->Draw("dz","lay == 1 && noverprim>0");
199 gPad->SetFillColor(33);
200 ntuple1->SetFillColor(46);
201 ntuple1->Draw("dz","lay == 2 && noverprim>0");
206 gPad->SetFillColor(33);
207 ntuple1->SetFillColor(42);
208 ntuple1->Draw("dx","lay == 1");
210 gPad->SetFillColor(33);
211 ntuple1->SetFillColor(46);
212 ntuple1->Draw("dx","lay == 2");
214 gPad->SetFillColor(33);
215 ntuple1->SetFillColor(42);
216 ntuple1->Draw("dz","lay == 1");
218 gPad->SetFillColor(33);
219 ntuple1->SetFillColor(46);
220 ntuple1->Draw("dz","lay == 2");
225 gPad->SetFillColor(33);
226 ntuple2->SetFillColor(42);
227 ntuple2->Draw("nx","lay == 1");
229 gPad->SetFillColor(33);
230 ntuple2->SetFillColor(46);
231 ntuple2->Draw("nz","lay == 1");
233 gPad->SetFillColor(33);
234 ntuple2->SetFillColor(42);
235 ntuple2->Draw("nx","lay == 2");
237 gPad->SetFillColor(33);
238 ntuple2->SetFillColor(46);
239 ntuple2->Draw("nz","lay == 2");
242 ///////////////////// Histogramm analysis ////////////////////////
246 gPad->SetFillColor(33);
247 Xres1->SetFillColor(42);
250 gPad->SetFillColor(33);
251 Zres1->SetFillColor(42);
254 gPad->SetFillColor(33);
255 Xres2->SetFillColor(46);
258 gPad->SetFillColor(33);
259 Zres2->SetFillColor(46);
264 gPad->SetFillColor(33);
265 Nzpix1->SetFillColor(42);
268 gPad->SetFillColor(33);
269 Nxpix1->SetFillColor(42);
272 gPad->SetFillColor(33);
273 Nzpix2->SetFillColor(46);
276 gPad->SetFillColor(33);
277 Nxpix2->SetFillColor(46);
283 gPad->SetFillColor(33);
284 Zpix1->SetFillColor(42);
287 gPad->SetFillColor(33);
288 Xpix1->SetFillColor(46);
291 gPad->SetFillColor(33);
292 Zpix2->SetFillColor(42);
295 gPad->SetFillColor(33);
296 Xpix2->SetFillColor(46);
302 c2->Print("spd_res.ps");
306 c2->Print("spd_clsize.ps");