]>
Commit | Line | Data |
---|---|---|
70f9c629 | 1 | void RICHgainvar (Int_t evNumber1=0,Int_t evNumber2=0) |
2 | { | |
3 | ||
4 | gClassTable->GetID("AliRun"); | |
5 | ||
6 | ||
7 | // Dynamically link some shared libs | |
8 | ||
9 | if (gClassTable->GetID("AliRun") < 0) { | |
10 | gROOT->LoadMacro("loadlibs.C"); | |
11 | loadlibs(); | |
12 | }else { | |
13 | delete gAlice; | |
14 | gAlice = 0; | |
15 | } | |
16 | ||
17 | gAlice=0; | |
18 | ||
19 | // Connect the Root Galice file containing Geometry, Kine and Hits | |
20 | ||
21 | TFile *file = (TFile*)gROOT->GetListOfFiles()->FindObject("galice.root"); | |
22 | if (!file) file = new TFile("galice.root","UPDATE"); | |
23 | ||
24 | // Get AliRun object from file or create it if not on file | |
25 | ||
26 | if (!gAlice) { | |
27 | gAlice = (AliRun*)file->Get("gAlice"); | |
28 | if (gAlice) printf("AliRun object found on file\n"); | |
29 | if (!gAlice) gAlice = new AliRun("gAlice","Alice test program"); | |
30 | } | |
31 | else { | |
32 | delete gAlice; | |
33 | gAlice = (AliRun*)file->Get("gAlice"); | |
34 | if (gAlice) printf("AliRun object found on file\n"); | |
35 | if (!gAlice) gAlice = new AliRun("gAlice","Alice test program"); | |
36 | } | |
37 | ||
38 | AliRICH *RICH = (AliRICH*)gAlice->GetDetector("RICH"); | |
39 | ||
40 | TH1F *gainX = new TH1F("gainX","Gain along the wires",40,0,40); | |
41 | TH1F *gainIn = new TH1F("gainIn","Gain near the center of photocathode",300,0,300); | |
42 | TH1F *gainOut = new TH1F("gainOut","Gain near the edge of photocathode",300,0,300); | |
7087b66c | 43 | /*TH1F *gainPhi1 = new TH1F("gainPhi1","Gain variation along wires, 0-30 degrees",300,0,300); |
70f9c629 | 44 | TH1F *gainPhi2 = new TH1F("gainPhi2","Gain variation along wires, 30-60 degrees",300,0,300); |
45 | TH1F *gainPhi3 = new TH1F("gainPhi3","Gain variation along wires, 60-90 degrees",300,0,300); | |
46 | TH1F *gainPhi4 = new TH1F("gainPhi4","Gain variation along wires, 90-120 degrees",300,0,300); | |
47 | TH1F *gainPhi5 = new TH1F("gainPhi5","Gain variation along wires, 120-150 degrees",300,0,300); | |
48 | TH1F *gainPhi6 = new TH1F("gainPhi6","Gain variation along wires, 150-180 degrees",300,0,300); | |
49 | TH1F *gainPhi7 = new TH1F("gainPhi7","Gain variation along wires, 180-210 degrees",300,0,300); | |
50 | TH1F *gainPhi8 = new TH1F("gainPhi8","Gain variation along wires, 210-240 degrees",300,0,300); | |
51 | TH1F *gainPhi9 = new TH1F("gainPhi9","Gain variation along wires, 240-270 degrees",300,0,300); | |
52 | TH1F *gainPhi10 = new TH1F("gainPhi10","Gain variation along wires, 270-300 degrees",300,0,300); | |
53 | TH1F *gainPhi11 = new TH1F("gainPhi11","Gain variation along wires, 300-330 degrees",300,0,300); | |
54 | TH1F *gainPhi12 = new TH1F("gainPhi12","Gain variation along wires, 330-360 degrees",300,0,300); | |
7087b66c | 55 | TH1F *gainPhi = new TH1F("gainPhi","Gain variation along phi",36,0,360);*/ |
70f9c629 | 56 | |
7087b66c | 57 | |
58 | TH1F *gainPhi1 = new TH1F("gainPhi1","Gain variation along wires, 315 to 45 degrees",300,0,300); | |
59 | TH1F *gainPhi2 = new TH1F("gainPhi2","Gain variation along wires, 45 to 135 degrees",300,0,300); | |
60 | TH1F *gainPhi3 = new TH1F("gainPhi3","Gain variation along wires, 135 to 225 degrees",300,0,300); | |
61 | TH1F *gainPhi4 = new TH1F("gainPhi4","Gain variation along wires, 225 to 315 degrees",300,0,300); | |
62 | TH1F *gainPhi = new TH1F("gainPhi","Gain variation along phi",8,0,360); | |
63 | ||
70f9c629 | 64 | // Start loop over events |
65 | ||
66 | for (int nev=0; nev<= evNumber2; nev++) { | |
67 | Int_t nparticles = gAlice->GetEvent(nev); | |
68 | ||
69 | ||
70 | printf ("\n**********************************\nProcessing Event: %d\n",nev); | |
71 | printf ("Particles : %d\n\n",nparticles); | |
72 | if (nev < evNumber1) continue; | |
73 | if (nparticles <= 0) return; | |
74 | ||
75 | // Get pointers to RICH detector and Hits containers | |
76 | ||
77 | ||
78 | TTree *TH = gAlice->TreeH(); | |
79 | Int_t ntracks = TH->GetEntries(); | |
80 | Int_t nent=(Int_t)gAlice->TreeR()->GetEntries(); | |
81 | gAlice->TreeR()->GetEvent(nent-1); | |
82 | TClonesArray *Rawclusters = RICH->RawClustAddress(2); // Raw clusters branch | |
83 | //printf ("Rawclusters:%p",Rawclusters); | |
84 | Int_t nrawclusters = Rawclusters->GetEntriesFast(); | |
85 | //printf (" nrawclusters:%d\n",nrawclusters); | |
86 | ||
87 | ||
88 | // Start loop on tracks in the hits containers | |
89 | Int_t Nc=0; | |
90 | for (Int_t track=0; track<ntracks;track++) { | |
91 | printf ("\nProcessing Track: %d\n",track); | |
92 | gAlice->ResetHits(); | |
93 | Int_t nbytes += TH->GetEvent(track); | |
94 | if (RICH) | |
95 | TClonesArray *Hits = RICH->Hits(); // Hits branch | |
96 | //see hits distribution | |
97 | ||
98 | Int_t nhits = Hits->GetEntriesFast(); | |
99 | printf("Hits : %d\n",nhits); | |
100 | for (Int_t hit=0;hit<nhits;hit++) { | |
101 | mHit = (AliRICHHit*) Hits->UncheckedAt(hit); | |
102 | Float_t mx = mHit->X(); // x-pos of hit | |
103 | Float_t my = mHit->Z(); // y-pos | |
104 | } | |
105 | ||
106 | ||
107 | if (nrawclusters) { | |
108 | printf("Raw Clusters : %d\n",nrawclusters); | |
109 | for (Int_t hit=0;hit<nrawclusters;hit++) { | |
110 | rcHit = (AliRICHRawCluster*) Rawclusters->UncheckedAt(hit); | |
111 | //Int_t nchamber = rcHit->fChamber; // chamber number | |
112 | //Int_t nhit = cHit->fHitNumber; // hit number | |
113 | Int_t qtot = rcHit->fQ; // charge | |
114 | Int_t fx = rcHit->fX; // x-position | |
115 | Int_t fy = rcHit->fY; // y-position | |
116 | Int_t mult = rcHit->fMultiplicity; // How many pads form the cluster | |
117 | ||
118 | Float_t radfid = TMath::Sqrt((TMath::Power(fx-mx,2)+(TMath::Power(fy-my,2)))); | |
119 | //printf("Radius:%f\n", radfid); | |
120 | ||
121 | ||
122 | if(qtot<200 && radfid>=9 && radfid<=15) | |
123 | { | |
124 | if (fy>=0) | |
125 | gainX->Fill(fy-40,(float) qtot); | |
126 | else | |
127 | gainX->Fill(fy+40,(float) qtot); | |
128 | ||
129 | if(fy<=0 && fy>=-4) | |
130 | gainOut->Fill(qtot,(float) 1); | |
131 | if(fy<=-11 && fy>=-15) | |
132 | gainIn->Fill(qtot,(float) 1); | |
133 | ||
134 | Float_t phi = TMath::ATan2(fx-mx, fy-my); | |
135 | Float_t phi_deg = phi/TMath::Pi()*180; | |
136 | ||
137 | if ((fx-mx)<0) | |
138 | { | |
139 | phi_deg = 360+phi_deg; | |
140 | //printf("X:%f, Y:%f, Phi:%f\n",-fy+my,fx-mx,phi_deg); | |
141 | } | |
142 | else | |
143 | { | |
144 | //printf("X:%f, Y:%f, Phi:%f\n",-fy+my,fx-mx,phi_deg); | |
145 | } | |
146 | ||
147 | ||
7087b66c | 148 | /*if (phi_deg>=0 && phi_deg<30) |
70f9c629 | 149 | gainPhi1->Fill(qtot,(float) 1); |
150 | if (phi_deg>=30 && phi_deg<60) | |
151 | gainPhi2->Fill(qtot,(float) 1); | |
152 | if (phi_deg>=60 && phi_deg<90) | |
153 | gainPhi3->Fill(qtot,(float) 1); | |
154 | if (phi_deg>=90 && phi_deg<120) | |
155 | gainPhi4->Fill(qtot,(float) 1); | |
156 | if (phi_deg>=120 && phi_deg<150) | |
157 | gainPhi5->Fill(qtot,(float) 1); | |
158 | if (phi_deg>=150 && phi_deg<180) | |
159 | gainPhi6->Fill(qtot,(float) 1); | |
160 | if (phi_deg>=180 && phi_deg<210) | |
161 | gainPhi7->Fill(qtot,(float) 1); | |
162 | if (phi_deg>=210 && phi_deg<240) | |
163 | gainPhi8->Fill(qtot,(float) 1); | |
164 | if (phi_deg>=240 && phi_deg<270) | |
165 | gainPhi9->Fill(qtot,(float) 1); | |
166 | if (phi_deg>=270 && phi_deg<300) | |
167 | gainPhi10->Fill(qtot,(float) 1); | |
168 | if (phi_deg>=300 && phi_deg<330) | |
169 | gainPhi11->Fill(qtot,(float) 1); | |
170 | if (phi_deg>=330 && phi_deg<360) | |
7087b66c | 171 | gainPhi12->Fill(qtot,(float) 1);*/ |
172 | ||
173 | if (phi_deg>=315 || phi_deg<45) | |
174 | gainPhi1->Fill(qtot,(float) 1); | |
175 | if (phi_deg>=45 && phi_deg<135) | |
176 | gainPhi2->Fill(qtot,(float) 1); | |
177 | if (phi_deg>=135 && phi_deg<225) | |
178 | gainPhi3->Fill(qtot,(float) 1); | |
179 | if (phi_deg>=225 && phi_deg<315) | |
180 | gainPhi4->Fill(qtot,(float) 1); | |
70f9c629 | 181 | |
182 | } | |
183 | ||
184 | //printf("Y:%d, Q:%d %d in %d\n",fy, qtot,hit,nrawclusters); | |
185 | } | |
186 | } | |
187 | } | |
188 | } | |
189 | ||
190 | ||
191 | TCanvas *c16 = new TCanvas("c16","Gain Variation over Phi",50,50,1200,700); | |
7087b66c | 192 | //c16->Divide(4,2); |
70f9c629 | 193 | //c16->SetFillColor(42); |
194 | ||
7087b66c | 195 | c16->Divide(2,2); |
196 | ||
70f9c629 | 197 | c16->cd(1); |
198 | c16_1->SetLogy(); | |
199 | gainPhi1->SetFillColor(5); | |
200 | gainPhi1->SetXTitle("ADC counts"); | |
201 | gainPhi1->Fit("expo","M"); | |
202 | expo->SetLineColor(2); | |
203 | expo->SetLineWidth(3); | |
204 | Float_t inv_slope1 = -1/expo->GetParameter(1); | |
7087b66c | 205 | //gainPhi->Fill(15, (float) inv_slope1); |
206 | gainPhi->Fill(0, (float) inv_slope1); | |
70f9c629 | 207 | gainPhi1->Draw(); |
208 | ||
209 | c16->cd(2); | |
210 | c16_2->SetLogy(); | |
211 | gainPhi2->SetFillColor(5); | |
212 | gainPhi2->SetXTitle("ADC counts"); | |
213 | gainPhi2->Fit("expo","M"); | |
214 | expo->SetLineColor(2); | |
215 | expo->SetLineWidth(3); | |
216 | Float_t inv_slope2 = -1/expo->GetParameter(1); | |
7087b66c | 217 | //gainPhi->Fill(45, (float) inv_slope2); |
218 | gainPhi->Fill(90, (float) inv_slope2); | |
70f9c629 | 219 | gainPhi2->Draw(); |
220 | ||
221 | c16->cd(3); | |
222 | c16_3->SetLogy(); | |
223 | gainPhi3->SetFillColor(5); | |
224 | gainPhi3->SetXTitle("ADC counts"); | |
225 | gainPhi3->Fit("expo","M"); | |
226 | expo->SetLineColor(2); | |
227 | expo->SetLineWidth(3); | |
228 | Float_t inv_slope3 = -1/expo->GetParameter(1); | |
7087b66c | 229 | //gainPhi->Fill(75, (float) inv_slope3); |
230 | gainPhi->Fill(180, (float) inv_slope3); | |
70f9c629 | 231 | gainPhi3->Draw(); |
232 | ||
233 | c16->cd(4); | |
234 | c16_4->SetLogy(); | |
235 | gainPhi4->SetFillColor(5); | |
236 | gainPhi4->SetXTitle("ADC counts"); | |
237 | gainPhi4->Fit("expo","M"); | |
238 | expo->SetLineColor(2); | |
239 | expo->SetLineWidth(3); | |
240 | Float_t inv_slope4 = -1/expo->GetParameter(1); | |
7087b66c | 241 | //gainPhi->Fill(105, (float) inv_slope4); |
242 | gainPhi->Fill(270, (float) inv_slope4); | |
70f9c629 | 243 | gainPhi4->Draw(); |
244 | ||
7087b66c | 245 | /*c16->cd(5); |
70f9c629 | 246 | c16_5->SetLogy(); |
247 | gainPhi5->SetFillColor(5); | |
248 | gainPhi5->SetXTitle("ADC counts"); | |
249 | gainPhi5->Fit("expo","M"); | |
250 | expo->SetLineColor(2); | |
251 | expo->SetLineWidth(3); | |
252 | Float_t inv_slope5 = -1/expo->GetParameter(1); | |
253 | gainPhi->Fill(135, (float) inv_slope5); | |
254 | gainPhi5->Draw(); | |
255 | ||
256 | c16->cd(6); | |
257 | c16_6->SetLogy(); | |
258 | gainPhi6->SetFillColor(5); | |
259 | gainPhi6->SetXTitle("ADC counts"); | |
260 | gainPhi6->Fit("expo","M"); | |
261 | expo->SetLineColor(2); | |
262 | expo->SetLineWidth(3); | |
263 | Float_t inv_slope6 = -1/expo->GetParameter(1); | |
264 | gainPhi->Fill(165, (float) inv_slope6); | |
265 | gainPhi6->Draw(); | |
266 | ||
267 | c16->cd(7); | |
268 | c16_7->SetLogy(); | |
269 | gainPhi7->SetFillColor(5); | |
270 | gainPhi7->SetXTitle("ADC counts"); | |
271 | gainPhi7->Fit("expo","M"); | |
272 | expo->SetLineColor(2); | |
273 | expo->SetLineWidth(3); | |
274 | Float_t inv_slope7 = -1/expo->GetParameter(1); | |
275 | gainPhi->Fill(195, (float) inv_slope7); | |
276 | gainPhi7->Draw(); | |
277 | ||
278 | c16->cd(8); | |
279 | c16_8->SetLogy(); | |
280 | gainPhi8->SetFillColor(5); | |
281 | gainPhi8->SetXTitle("ADC counts"); | |
282 | gainPhi8->Fit("expo","M"); | |
283 | expo->SetLineColor(2); | |
284 | expo->SetLineWidth(3); | |
285 | Float_t inv_slope8 = -1/expo->GetParameter(1); | |
286 | gainPhi->Fill(225, (float) inv_slope8); | |
287 | gainPhi8->Draw(); | |
288 | ||
289 | c16->cd(9); | |
290 | c16_9->SetLogy(); | |
291 | gainPhi9->SetFillColor(5); | |
292 | gainPhi9->SetXTitle("ADC counts"); | |
293 | gainPhi9->Fit("expo","M"); | |
294 | expo->SetLineColor(2); | |
295 | expo->SetLineWidth(3); | |
296 | Float_t inv_slope9 = -1/expo->GetParameter(1); | |
297 | gainPhi->Fill(255, (float) inv_slope9); | |
298 | c16->Update(); | |
299 | gainPhi9->Draw(); | |
300 | ||
301 | ||
302 | c16->cd(10); | |
303 | c16_10->SetLogy(); | |
304 | gainPhi10->SetFillColor(5); | |
305 | gainPhi10->SetXTitle("ADC counts"); | |
306 | gainPhi10->Fit("expo","M"); | |
307 | expo->SetLineColor(2); | |
308 | expo->SetLineWidth(3); | |
309 | Float_t inv_slope10 = -1/expo->GetParameter(1); | |
310 | gainPhi->Fill(285, (float) inv_slope10); | |
311 | gainPhi10->Draw(); | |
312 | ||
313 | c16->cd(11); | |
314 | c16_11->SetLogy(); | |
315 | gainPhi11->SetFillColor(5); | |
316 | gainPhi11->SetXTitle("ADC counts"); | |
317 | gainPhi11->Fit("expo","M"); | |
318 | expo->SetLineColor(2); | |
319 | expo->SetLineWidth(3); | |
320 | Float_t inv_slope11 = -1/expo->GetParameter(1); | |
321 | gainPhi->Fill(315, (float) inv_slope11); | |
322 | gainPhi11->Draw(); | |
323 | ||
324 | c16->cd(12); | |
325 | c16_12->SetLogy(); | |
326 | gainPhi12->SetFillColor(5); | |
327 | gainPhi12->SetXTitle("ADC counts"); | |
328 | gainPhi12->Fit("expo","M"); | |
329 | expo->SetLineColor(2); | |
330 | expo->SetLineWidth(3); | |
331 | Float_t inv_slope12 = -1/expo->GetParameter(1); | |
332 | gainPhi->Fill(335, (float) inv_slope12); | |
7087b66c | 333 | gainPhi12->Draw();*/ |
70f9c629 | 334 | |
335 | TCanvas *c17 = new TCanvas("c17","Gain Variation",50,50,600,700); | |
336 | c17->Divide(2,2); | |
337 | ||
338 | //c11->SetFillColor(42); | |
339 | ||
340 | c17->cd(1); | |
341 | c17_1->SetLogy(); | |
342 | gainIn->SetFillColor(5); | |
343 | gainIn->SetXTitle("ADC counts"); | |
344 | gainIn->Fit("expo","M"); | |
345 | expo->SetLineColor(2); | |
346 | expo->SetLineWidth(3); | |
347 | gainIn->Draw(); | |
348 | ||
349 | c17->cd(2); | |
350 | c17_2->SetLogy(); | |
351 | gainOut->SetFillColor(5); | |
352 | gainOut->SetXTitle("ADC counts"); | |
353 | gainOut->Fit("expo","M"); | |
354 | expo->SetLineColor(2); | |
355 | expo->SetLineWidth(3); | |
356 | gainOut->Draw(); | |
357 | ||
358 | c17->cd(3); | |
359 | gainX->SetFillColor(5); | |
360 | gainX->SetYTitle("Total ADC counts"); | |
361 | gainX->SetXTitle("(pads)"); | |
362 | gainX->Draw(); | |
363 | ||
364 | c17->cd(4); | |
365 | gainPhi->SetFillColor(5); | |
366 | gainPhi->SetYTitle("Gain (ADC)"); | |
367 | gainPhi->SetXTitle("phi (degrees)"); | |
368 | gainPhi->Draw(); | |
369 | ||
370 | printf("\nEnd of macro\n"); | |
371 | printf("**********************************\n"); | |
372 | ||
373 | } |