]> git.uio.no Git - u/mrichter/AliRoot.git/blame - ITS/SSD_ntuple.C
git://git.uio.no / u / mrichter / AliRoot.git / blame
? search:
summary | shortlog | log | commit | commitdiff | tree
blob | blame (incremental) | history | HEAD
Use TMath::Abs instead of fabs (Alpha)
[u/mrichter/AliRoot.git] / ITS / SSD_ntuple.C
CommitLineData
6b8f55ce 1void SSD_ntuple()
2{
3
4TFile *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 9TCanvas *c2 = new TCanvas("c2","SSD clusters",400,10,600,700);
6b8f55ce 10//c1->Divide(2,2);
11c2->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/*
45c2->cd(1);
46gPad->SetFillColor(33);
47 ntuple->SetFillColor(42);
48 ntuple->Draw("dx","lay == 5 && hitprim == 1");
49c2->cd(2);
50gPad->SetFillColor(33);
51 ntuple->SetFillColor(46);
52 ntuple->Draw("dz","lay == 5 && hitprim == 1");
53c2->cd(3);
54gPad->SetFillColor(33);
55 ntuple->SetFillColor(42);
56 ntuple->Draw("dx","lay == 6 && hitprim == 1");
57c2->cd(4);
58gPad->SetFillColor(33);
59 ntuple->SetFillColor(46);
60 ntuple->Draw("dz","lay == 6 && hitprim == 1");
61*/
6b8f55ce 62
6298ec22 63/*
64c2->cd(1);
65gPad->SetFillColor(33);
66 ntuple->SetFillColor(42);
67 ntuple->Draw("dx","lay == 5 && hitprim == 1&&ntrover==1");
68c2->cd(2);
69gPad->SetFillColor(33);
70 ntuple->SetFillColor(46);
71 ntuple->Draw("dz","lay == 5 && hitprim == 1&&ntrover==1");
72c2->cd(3);
73gPad->SetFillColor(33);
74 ntuple->SetFillColor(42);
75 ntuple->Draw("dx","lay == 6 && hitprim == 1&&ntrover==1");
76c2->cd(4);
77gPad->SetFillColor(33);
78 ntuple->SetFillColor(46);
79 ntuple->Draw("dz","lay == 6 && hitprim == 1&&ntrover==1");
80*/
81
82
83/*
6b8f55ce 84c2->cd(1);
85gPad->SetFillColor(33);
86 ntuple->SetFillColor(42);
87 ntuple->Draw("dx","lay == 5 && hitprim == 1&&abs(dx)<200");
88c2->cd(2);
89gPad->SetFillColor(33);
90 ntuple->SetFillColor(46);
91 ntuple->Draw("dz","lay == 5 && hitprim == 1&&abs(dz)<5000");
92c2->cd(3);
93gPad->SetFillColor(33);
94 ntuple->SetFillColor(42);
95 ntuple->Draw("dx","lay == 6 && hitprim == 1&&abs(dx)<200");
96c2->cd(4);
97gPad->SetFillColor(33);
98 ntuple->SetFillColor(46);
99 ntuple->Draw("dz","lay == 6 && hitprim == 1&&abs(dz)<5000");
6298ec22 100*/
6b8f55ce 101
6298ec22 102/*
103c2->cd(1);
104gPad->SetFillColor(33);
105 ntuple->SetFillColor(42);
106 ntuple->Draw("dx","lay == 5 && hitprim == 1&&abs(dx)<200&&ntrover==1");
107c2->cd(2);
108gPad->SetFillColor(33);
109 ntuple->SetFillColor(46);
110 ntuple->Draw("dz","lay == 5 && hitprim == 1&&abs(dz)<5000&&ntrover==1");
111c2->cd(3);
112gPad->SetFillColor(33);
113 ntuple->SetFillColor(42);
114 ntuple->Draw("dx","lay == 6 && hitprim == 1&&abs(dx)<200&&ntrover==1");
115c2->cd(4);
116gPad->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 122c2->cd(1);
123gPad->SetFillColor(33);
124 ntuple1->SetFillColor(42);
125 ntuple1->Draw("nxP","lay == 5&&noverprim>=0");
126c2->cd(2);
127gPad->SetFillColor(33);
128 ntuple1->SetFillColor(46);
129 ntuple1->Draw("nxN","lay == 5&&noverprim>=0");
130c2->cd(3);
131gPad->SetFillColor(33);
132 ntuple1->SetFillColor(42);
133 ntuple1->Draw("nxP","lay == 6 && noverprim>=0");
134c2->cd(4);
135gPad->SetFillColor(33);
136 ntuple1->SetFillColor(46);
137 ntuple1->Draw("nxN","lay == 6 && noverprim>=0");
0315d466 138*/
6b8f55ce 139
6298ec22 140/*
141c2->cd(1);
142gPad->SetFillColor(33);
143 ntuple1->SetFillColor(42);
144 ntuple1->Draw("dx","lay == 5&&noverprim>=0&&abs(dx)<50&&nxP==2&&nxN==2");
145c2->cd(2);
146gPad->SetFillColor(33);
147 ntuple1->SetFillColor(46);
148 ntuple1->Draw("dz","lay == 5&&noverprim>=0&&abs(dz)<2000&&nxP==2&&nxN==2");
149c2->cd(3);
150gPad->SetFillColor(33);
151 ntuple1->SetFillColor(42);
152 ntuple1->Draw("dx","lay == 6 && noverprim>=0&&abs(dx)<200");
153c2->cd(4);
154gPad->SetFillColor(33);
155 ntuple1->SetFillColor(46);
156 ntuple1->Draw("dz","lay == 6 && noverprim>=0&&abs(dz)<5000");
157*/
158
159/*
160c2->cd(1);
161gPad->SetFillColor(33);
162 ntuple1->SetFillColor(42);
163 ntuple1->Draw("dx","lay == 5&&noverprim>=0&&noverlaps==0&&qcut<0.18");
164c2->cd(2);
165gPad->SetFillColor(33);
166 ntuple1->SetFillColor(46);
167 ntuple1->Draw("dz","lay == 5&&noverprim>=0&&noverlaps==0&&qcut<0.18");
168c2->cd(3);
169gPad->SetFillColor(33);
170 ntuple1->SetFillColor(42);
171 ntuple1->Draw("dx","lay == 6 && noverprim>=0&&noverlaps==0&&qcut<0.18");
172c2->cd(4);
173gPad->SetFillColor(33);
174 ntuple1->SetFillColor(46);
175 ntuple1->Draw("dz","lay == 6 && noverprim>=0&&noverlaps==0&&qcut<0.18");
176*/
177
178/*
179c2->cd(1);
180gPad->SetFillColor(33);
181 ntuple1->SetFillColor(42);
182 ntuple1->Draw("dx","lay == 5&&noverprim>=0&&noverlaps==0&&abs(dx)<100&&qcut<0.18");
183c2->cd(2);
184gPad->SetFillColor(33);
185 ntuple1->SetFillColor(46);
186 ntuple1->Draw("dz","lay == 5&&noverprim>=0&&noverlaps==0&&abs(dz)<5000&&qcut<0.18");
187c2->cd(3);
188gPad->SetFillColor(33);
189 ntuple1->SetFillColor(42);
190 ntuple1->Draw("dx","lay == 6&&noverprim>=0&&noverlaps==0&&abs(dx)<100&&qcut<0.18");
191c2->cd(4);
192gPad->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/*
198c2->cd(1);
199gPad->SetFillColor(33);
200 ntuple1->SetFillColor(42);
201 ntuple1->Draw("noverprim","lay == 5&&noverprim>=-1");
202c2->cd(2);
203gPad->SetFillColor(33);
204 ntuple1->SetFillColor(46);
205 ntuple1->Draw("noverlaps","lay == 5&&noverprim>=-1");
206c2->cd(3);
207gPad->SetFillColor(33);
208 ntuple1->SetFillColor(42);
209 ntuple1->Draw("noverprim","lay == 6 && noverprim>=-1");
210c2->cd(4);
211gPad->SetFillColor(33);
212 ntuple1->SetFillColor(46);
213 ntuple1->Draw("noverlaps","lay == 6 && noverprim>=-1");
214*/
215
216 /*
217c2->cd(1);
218gPad->SetFillColor(33);
219 ntuple1->SetFillColor(42);
220 ntuple1->Draw("ntrover","lay == 5&&noverprim>=0");
221c2->cd(2);
222gPad->SetFillColor(33);
223 ntuple1->SetFillColor(46);
224 ntuple1->Draw("noverlaps","lay == 5&&noverprim>=0");
225c2->cd(3);
226gPad->SetFillColor(33);
227 ntuple1->SetFillColor(42);
228 ntuple1->Draw("ntrover","lay == 6 && noverprim>=0");
229c2->cd(4);
230gPad->SetFillColor(33);
231 ntuple1->SetFillColor(46);
232 ntuple1->Draw("noverlaps","lay == 6 && noverprim>=0");
233 */
234
6b8f55ce 235/*
236c2->cd(1);
237gPad->SetFillColor(33);
238 ntuple1->SetFillColor(42);
239 ntuple1->Draw("qclP","lay == 5&&noverprim>=0");
240c2->cd(2);
241gPad->SetFillColor(33);
242 ntuple1->SetFillColor(42);
243 ntuple1->Draw("qclN","lay == 5&&noverprim>=0");
244c2->cd(3);
245gPad->SetFillColor(33);
246 ntuple1->SetFillColor(46);
247 ntuple1->Draw("qclP","lay == 6&&noverprim>=0");
248c2->cd(4);
249gPad->SetFillColor(33);
250 ntuple1->SetFillColor(46);
251 ntuple1->Draw("qclN","lay == 6&&noverprim>=0");
252*/
253
254
255///////////////////// Histogramm/ntuple analysis ////////////////////////
256
6298ec22 257
258c2->cd(1);
259gPad->SetFillColor(33);
260 adcPadcN5all->SetFillColor(42);
261 adcPadcN5all->Draw();
262c2->cd(2);
263gPad->SetFillColor(33);
264 adcPadcN6all->SetFillColor(42);
265 adcPadcN6all->Draw();
266c2->cd(3);
267gPad->SetFillColor(33);
268 adcPadcN5cut->SetFillColor(42);
269 adcPadcN5cut->Draw();
270c2->cd(4);
271gPad->SetFillColor(33);
272 adcPadcN6cut->SetFillColor(46);
273 adcPadcN6cut->Draw();
274
275
6b8f55ce 276 /*
277c2->cd(1);
278gPad->SetFillColor(33);
0315d466 279 ntuple1->SetFillColor(42);
6298ec22 280 ntuple1->Draw("qclP","noverprim>=0&&qclP<500");
6b8f55ce 281c2->cd(2);
282gPad->SetFillColor(33);
0315d466 283 ntuple1->SetFillColor(42);
6298ec22 284 ntuple1->Draw("qclN","noverprim>=0&&qclN<500");
6b8f55ce 285c2->cd(3);
286gPad->SetFillColor(33);
287 adcPadcN5cut->SetFillColor(42);
288 adcPadcN5cut->Draw();
289c2->cd(4);
290gPad->SetFillColor(33);
291 adcPadcN6cut->SetFillColor(46);
292 adcPadcN6cut->Draw();
293 */
294
295
296 /*
297c2->Draw();
298c2->Print("ssd_res.ps");
299 */
300
301 /*
302c2->Draw();
303c2->Print("spd_clsize.ps");
304 */
305
306}
307
308
309
310
311
312
313
314
315