]>
Commit | Line | Data |
---|---|---|
91282711 | 1 | #if !defined(__CINT__) || defined(__MAKECINT__) |
2 | #include <TH2F.h> | |
3 | #include <TCanvas.h> | |
4 | #include <TStopwatch.h> | |
5 | #include <TStyle.h> | |
6 | #include <TGrid.h> | |
7 | #include <TSystem.h> | |
8 | #include <TLatex.h> | |
9 | #include <TFile.h> | |
10 | #include "AliRawReader.h" | |
11 | #include "AliRawReaderDate.h" | |
12 | #include "AliRawReaderRoot.h" | |
13 | #include "AliITSOnlineSDDBase.h" | |
14 | #include "AliITSOnlineSDDCMN.h" | |
15 | #include "AliITSRawStreamSDD.h" | |
a8cd2ab2 | 16 | #include "AliITSRawStreamSDDCompressed.h" |
91282711 | 17 | #include "TPaveStats.h" |
18 | #endif | |
19 | ||
20 | // Macro for the analysis of PEDESTAL runs (equivalent to ITSSDDBASda.cxx) | |
21 | // Two functions named AnalyzeSDDNoiseAllModules: | |
22 | // The first is for analyzing a local raw data file and takes as agrument the file name. | |
23 | // The second is for running on ALIEN | |
24 | // All DDLs are analyzed, the argument nDDL selects the DDL to be plotted | |
25 | // Origin: F. Prino (prino@to.infn.it) | |
26 | ||
a8cd2ab2 | 27 | void AnalyzeSDDNoiseAllMod(Char_t *datafil, |
28 | Int_t adcfreq=20, | |
29 | Int_t nDDL=0, | |
e5b2f7f2 | 30 | Int_t firstEv=18, |
31 | Int_t lastEv=20){ | |
a8cd2ab2 | 32 | |
5a4e63d5 | 33 | |
91282711 | 34 | const Int_t kTotDDL=24; |
35 | const Int_t kModPerDDL=12; | |
36 | const Int_t kSides=2; | |
e5b2f7f2 | 37 | Bool_t writtenoutput=kFALSE; |
91282711 | 38 | |
39 | AliITSOnlineSDDBase **base=new AliITSOnlineSDDBase*[kTotDDL*kModPerDDL*kSides]; | |
40 | TH2F **histo=new TH2F*[kTotDDL*kModPerDDL*kSides]; | |
41 | ||
42 | Char_t hisnam[20]; | |
43 | for(Int_t iddl=0; iddl<kTotDDL;iddl++){ | |
44 | for(Int_t imod=0; imod<kModPerDDL;imod++){ | |
45 | for(Int_t isid=0;isid<kSides;isid++){ | |
46 | Int_t index=kSides*(kModPerDDL*iddl+imod)+isid; | |
b255f56d | 47 | base[index]=new AliITSOnlineSDDBase(iddl,imod,isid); |
48 | if(adcfreq==40) base[index]->SetLastGoodTB(254); | |
49 | else base[index]->SetLastGoodTB(126); | |
91282711 | 50 | sprintf(hisnam,"h%02dc%02ds%d",iddl,imod,isid); |
e8fc1111 | 51 | histo[index]=new TH2F(hisnam,"",256,-0.5,255.5,256,-0.5,255.5); |
91282711 | 52 | } |
53 | } | |
54 | } | |
55 | ||
56 | ||
57 | TCanvas* c0 = new TCanvas("c0","Ev Display",900,900); | |
58 | gStyle->SetPalette(1); | |
59 | Char_t text[50]; | |
60 | ||
61 | Int_t iev=firstEv; | |
62 | AliRawReader *rd; | |
63 | if(strstr(datafil,".root")!=0){ | |
64 | rd=new AliRawReaderRoot(datafil,iev); | |
65 | }else{ | |
66 | rd=new AliRawReaderDate(datafil,iev); | |
67 | } | |
68 | TLatex *t0=new TLatex(); | |
69 | t0->SetNDC(); | |
70 | t0->SetTextSize(0.06); | |
71 | t0->SetTextColor(4); | |
72 | ||
73 | do{ | |
74 | c0->Clear(); | |
75 | c0->Divide(4,6,0.001,0.001); | |
76 | printf("Event # %d ",iev); | |
91282711 | 77 | rd->Reset(); |
78 | for(Int_t iddl=0; iddl<kTotDDL;iddl++){ | |
79 | for(Int_t imod=0; imod<kModPerDDL;imod++){ | |
80 | for(Int_t isid=0;isid<kSides;isid++){ | |
81 | Int_t index=kSides*(kModPerDDL*iddl+imod)+isid; | |
82 | histo[index]->Reset(); | |
83 | } | |
84 | } | |
85 | } | |
e5b2f7f2 | 86 | |
87 | UChar_t cdhAttr=AliITSRawStreamSDD::ReadBlockAttributes(rd); | |
88 | UInt_t amSamplFreq=AliITSRawStreamSDD::ReadAMSamplFreqFromCDH(cdhAttr); | |
89 | AliITSRawStream* s=AliITSRawStreamSDD::CreateRawStreamSDD(rd,cdhAttr); | |
90 | if(!writtenoutput){ | |
91 | printf("Use %s raw stream, sampling frequency %d MHz\n",s->ClassName(),amSamplFreq); | |
92 | writtenoutput=kTRUE; | |
a8cd2ab2 | 93 | } |
e5b2f7f2 | 94 | |
a8cd2ab2 | 95 | while(s->Next()){ |
91282711 | 96 | Int_t iDDL=rd->GetDDLID(); |
a8cd2ab2 | 97 | Int_t iCarlos=s->GetCarlosId(); |
98 | if(s->IsCompletedModule()) continue; | |
99 | if(s->IsCompletedDDL()) continue; | |
de075dae | 100 | if(iDDL>=0 && iDDL<kTotDDL){ |
a8cd2ab2 | 101 | Int_t index=kSides*(kModPerDDL*iDDL+iCarlos)+s->GetChannel(); |
102 | histo[index]->Fill(s->GetCoord2(),s->GetCoord1(),s->GetSignal()); | |
91282711 | 103 | } |
104 | } | |
a8cd2ab2 | 105 | delete s; |
91282711 | 106 | iev++; |
107 | for(Int_t iddl=0; iddl<kTotDDL;iddl++){ | |
108 | for(Int_t imod=0; imod<kModPerDDL;imod++){ | |
109 | for(Int_t isid=0;isid<kSides;isid++){ | |
110 | Int_t index=kSides*(kModPerDDL*iddl+imod)+isid; | |
111 | base[index]->AddEvent(histo[index]); | |
112 | if(iddl==nDDL){ | |
113 | Int_t index2=kSides*imod+isid; | |
114 | c0->cd(index2+1); | |
115 | histo[index]->DrawCopy("colz"); | |
116 | sprintf(text,"DDL %d channel %d Side %d",nDDL,imod,isid); | |
117 | t0->DrawLatex(0.15,0.92,text); | |
118 | c0->Update(); | |
119 | } | |
120 | } | |
121 | } | |
122 | } | |
123 | printf(" --- OK\n"); | |
124 | }while(rd->NextEvent()&&iev<=lastEv); | |
125 | ||
126 | for(Int_t iddl=0; iddl<kTotDDL;iddl++){ | |
127 | for(Int_t imod=0; imod<kModPerDDL;imod++){ | |
128 | for(Int_t isid=0;isid<kSides;isid++){ | |
129 | Int_t index=kSides*(kModPerDDL*iddl+imod)+isid; | |
130 | base[index]->ValidateAnodes(); | |
131 | base[index]->WriteToASCII(); // fondamentale!!!!!!!!! | |
132 | delete base[index]; | |
133 | } | |
134 | } | |
135 | } | |
136 | delete rd; | |
137 | delete [] base; | |
138 | ||
139 | printf("Start second analysis for Common Mode correction\n"); | |
140 | AliITSOnlineSDDCMN **corr=new AliITSOnlineSDDCMN*[kTotDDL*kModPerDDL*kSides]; | |
141 | Bool_t isFilled[kTotDDL*kModPerDDL*kSides]; | |
142 | ||
143 | for(Int_t iddl=0; iddl<kTotDDL;iddl++){ | |
144 | for(Int_t imod=0; imod<kModPerDDL;imod++){ | |
145 | for(Int_t isid=0;isid<kSides;isid++){ | |
146 | Int_t index=kSides*(kModPerDDL*iddl+imod)+isid; | |
147 | corr[index]=new AliITSOnlineSDDCMN(iddl,imod,isid); | |
b255f56d | 148 | if(adcfreq==40) corr[index]->SetLastGoodTB(254); |
149 | else corr[index]->SetLastGoodTB(126); | |
91282711 | 150 | isFilled[index]=0; |
151 | } | |
152 | } | |
153 | } | |
154 | ||
155 | iev=firstEv; | |
156 | AliRawReader *rd2; | |
157 | if(strstr(datafil,".root")!=0){ | |
158 | rd2=new AliRawReaderRoot(datafil,iev); | |
159 | }else{ | |
160 | rd2=new AliRawReaderDate(datafil,iev); | |
161 | } | |
162 | do{ | |
163 | c0->Clear(); | |
164 | c0->Divide(4,6,0.001,0.001); | |
165 | printf("Event # %d ",iev); | |
91282711 | 166 | rd2->Reset(); |
167 | for(Int_t iddl=0; iddl<kTotDDL;iddl++){ | |
168 | for(Int_t imod=0; imod<kModPerDDL;imod++){ | |
169 | for(Int_t isid=0;isid<kSides;isid++){ | |
170 | Int_t index=kSides*(kModPerDDL*iddl+imod)+isid; | |
171 | histo[index]->Reset(); | |
172 | } | |
173 | } | |
174 | } | |
175 | ||
e5b2f7f2 | 176 | UChar_t cdhAttr=AliITSRawStreamSDD::ReadBlockAttributes(rd2); |
177 | UInt_t amSamplFreq=AliITSRawStreamSDD::ReadAMSamplFreqFromCDH(cdhAttr); | |
178 | AliITSRawStream* s=AliITSRawStreamSDD::CreateRawStreamSDD(rd2,cdhAttr); | |
179 | if(!writtenoutput){ | |
180 | printf("Use %s raw stream, sampling frequency %d MHz\n",s->ClassName(),amSamplFreq); | |
181 | writtenoutput=kTRUE; | |
a8cd2ab2 | 182 | } |
183 | while(s->Next()){ | |
91282711 | 184 | Int_t iDDL=rd2->GetDDLID(); |
a8cd2ab2 | 185 | Int_t iCarlos=s->GetCarlosId(); |
186 | if(s->IsCompletedModule()) continue; | |
187 | if(s->IsCompletedDDL()) continue; | |
de075dae | 188 | if(iDDL>=0 && iDDL<kTotDDL){ |
a8cd2ab2 | 189 | Int_t index=kSides*(kModPerDDL*iDDL+iCarlos)+s->GetChannel(); |
190 | histo[index]->Fill(s->GetCoord2(),s->GetCoord1(),s->GetSignal()); | |
91282711 | 191 | isFilled[index]=1; |
192 | } | |
193 | } | |
a8cd2ab2 | 194 | delete s; |
91282711 | 195 | iev++; |
196 | for(Int_t iddl=0; iddl<kTotDDL;iddl++){ | |
197 | for(Int_t imod=0; imod<kModPerDDL;imod++){ | |
198 | for(Int_t isid=0;isid<kSides;isid++){ | |
199 | Int_t index=kSides*(kModPerDDL*iddl+imod)+isid; | |
200 | if(isFilled[index]) corr[index]->AddEvent(histo[index]); | |
201 | if(iddl==nDDL){ | |
202 | Int_t index2=kSides*imod+isid; | |
203 | c0->cd(index2+1); | |
204 | histo[index]->DrawCopy("colz"); | |
205 | sprintf(text,"DDL %d channel %d Side %d",nDDL,imod,isid); | |
206 | t0->DrawLatex(0.15,0.92,text); | |
207 | c0->Update(); | |
208 | } | |
209 | } | |
210 | } | |
211 | } | |
212 | printf(" --- OK\n"); | |
213 | }while(rd2->NextEvent()&&iev<=lastEv); | |
214 | ||
215 | TH1F *htotbas=new TH1F("htotbas","",100,0.,150.); | |
216 | TH1F *htotbaseq=new TH1F("htotbaseq","",100,0.,150.); | |
217 | TH1F *htotnoise=new TH1F("htotnoise","",100,0.,10.); | |
218 | TH1F *htotnoisecorr=new TH1F("htotnoisecorr","",100,0.,10.); | |
219 | TH1F *hstatus=new TH1F("hstatus","",2,-0.5,1.5); | |
220 | ||
221 | TFile *outfil=new TFile("SDDbase-results.root","recreate"); | |
222 | for(Int_t iddl=0; iddl<kTotDDL;iddl++){ | |
223 | for(Int_t imod=0; imod<kModPerDDL;imod++){ | |
224 | for(Int_t isid=0;isid<kSides;isid++){ | |
225 | Int_t index=kSides*(kModPerDDL*iddl+imod)+isid; | |
226 | if(isFilled[index]){ | |
227 | corr[index]->ValidateAnodes(); | |
228 | corr[index]->WriteToASCII(); | |
229 | corr[index]->WriteToROOT(outfil); | |
230 | for(Int_t ian=0; ian<256;ian++){ | |
231 | Float_t basl=corr[index]->GetAnodeBaseline(ian); | |
232 | Float_t basleq=corr[index]->GetAnodeEqualizedBaseline(ian); | |
233 | Float_t noi=corr[index]->GetAnodeRawNoise(ian); | |
234 | Float_t cornoi=corr[index]->GetAnodeCorrNoise(ian); | |
235 | Int_t anstatus=corr[index]->IsAnodeGood(ian); | |
236 | hstatus->Fill(anstatus); | |
237 | htotbas->Fill(basl); | |
238 | htotbaseq->Fill(basleq); | |
239 | htotnoise->Fill(noi); | |
240 | htotnoisecorr->Fill(cornoi); | |
241 | } | |
242 | } | |
243 | } | |
244 | } | |
245 | } | |
246 | outfil->Close(); | |
247 | // Draw Statistics of baselines and noise | |
248 | TCanvas *call=new TCanvas("call","General stats",700,700); | |
249 | call->Divide(2,2); | |
250 | call->cd(1); | |
251 | htotbas->Draw(); | |
252 | htotbas->GetXaxis()->SetTitle("Baselines"); | |
253 | htotbas->GetXaxis()->SetTitleSize(0.07); | |
254 | htotbas->GetXaxis()->SetTitleOffset(0.6); | |
255 | call->cd(2); | |
256 | htotbaseq->Draw(); | |
257 | htotbaseq->GetXaxis()->SetTitle("Baselines after equalization"); | |
258 | htotbaseq->GetXaxis()->SetTitleSize(0.07); | |
259 | htotbaseq->GetXaxis()->SetTitleOffset(0.6); | |
260 | call->cd(3); | |
261 | htotnoisecorr->SetLineColor(2); | |
262 | htotnoisecorr->Draw(); | |
263 | call->Update(); | |
264 | TPaveStats *st1=(TPaveStats*)htotnoisecorr->GetListOfFunctions()->FindObject("stats"); | |
265 | st1->SetY1NDC(0.51); | |
266 | st1->SetY2NDC(0.7); | |
267 | htotnoisecorr->GetXaxis()->SetTitle("Noise"); | |
268 | htotnoisecorr->GetXaxis()->SetTitleSize(0.07); | |
269 | htotnoisecorr->GetXaxis()->SetTitleOffset(0.6); | |
270 | htotnoise->Draw("SAMES"); | |
271 | call->Update(); | |
272 | TPaveStats *st2=(TPaveStats*)htotnoise->GetListOfFunctions()->FindObject("stats"); | |
273 | st2->SetY1NDC(0.71); | |
274 | st2->SetY2NDC(0.9); | |
275 | ||
276 | call->cd(4); | |
277 | hstatus->Draw(); | |
278 | hstatus->GetXaxis()->SetTitle("Anode Status (0=bad 1=good)"); | |
279 | hstatus->GetXaxis()->SetTitleSize(0.07); | |
280 | hstatus->GetXaxis()->SetTitleOffset(0.6); | |
281 | call->Update(); | |
282 | call->SaveAs("GenStatsPedestal.gif"); | |
283 | ||
284 | // Draw baselines and noise for all modules of the selected DDL | |
285 | ||
286 | TH1F** hbas = new TH1F*[kSides*kModPerDDL]; | |
287 | TH1F** hrawn = new TH1F*[kSides*kModPerDDL]; | |
288 | TH1F** hcorrn = new TH1F*[kSides*kModPerDDL]; | |
289 | TH1F** hdbas = new TH1F*[kSides*kModPerDDL]; | |
290 | TH1F** hdrawn = new TH1F*[kSides*kModPerDDL]; | |
291 | TH1F** hdcorrn = new TH1F*[kSides*kModPerDDL]; | |
292 | TCanvas *c1=new TCanvas("c1","DDL: Baselines vs anode",900,900); | |
293 | c1->SetBottomMargin(0.14); | |
294 | c1->Divide(4,6,0.001,0.001); | |
295 | TCanvas *c2=new TCanvas("c2","DDL: Noise vs anode",900,900); | |
296 | c2->SetBottomMargin(0.14); | |
297 | c2->Divide(4,6,0.001,0.001); | |
298 | TCanvas *c3=new TCanvas("c3","DDL: Baselines distr",900,900); | |
299 | c3->SetBottomMargin(0.14); | |
300 | c3->Divide(4,6,0.001,0.001); | |
301 | TCanvas *c4=new TCanvas("c4","DDL: Noise Distr",900,900); | |
302 | c4->SetBottomMargin(0.14); | |
303 | c4->Divide(4,6,0.001,0.001); | |
304 | TLatex *t1=new TLatex(0.15,0.2,"Raw Noise"); | |
305 | t1->SetNDC(); | |
306 | t1->SetTextSize(0.05); | |
307 | TLatex *t2=new TLatex(0.4,0.2,"Corrected Noise"); | |
308 | t2->SetNDC(); | |
309 | t2->SetTextSize(0.05); | |
310 | t2->SetTextColor(2); | |
311 | TLatex *t3=new TLatex(); | |
312 | t3->SetNDC(); | |
313 | t3->SetTextSize(0.06); | |
314 | t3->SetTextColor(4); | |
315 | ||
316 | for(Int_t imod=0; imod<kModPerDDL;imod++){ | |
317 | for(Int_t isid=0;isid<kSides;isid++){ | |
318 | Int_t index1=kSides*(kModPerDDL*nDDL+imod)+isid; | |
319 | Int_t index2=kSides*imod+isid; | |
320 | sprintf(text,"DDL %d channel %d Side %d",nDDL,imod,isid); | |
321 | hbas[index2]=corr[index1]->GetBaselineAnodeHisto(); | |
322 | hrawn[index2]=corr[index1]->GetRawNoiseAnodeHisto(); | |
323 | hcorrn[index2]=corr[index1]->GetCorrNoiseAnodeHisto(); | |
324 | hdbas[index2]=corr[index1]->GetBaselineHisto(); | |
325 | hdrawn[index2]=corr[index1]->GetRawNoiseHisto(); | |
326 | hdcorrn[index2]=corr[index1]->GetCorrNoiseHisto(); | |
327 | c1->cd(index2+1); | |
328 | hbas[index2]->Draw(); | |
329 | hbas[index2]->SetMinimum(0); | |
330 | hbas[index2]->SetMaximum(75); | |
331 | hbas[index2]->GetXaxis()->SetTitle("Anode"); | |
332 | hbas[index2]->GetYaxis()->SetTitle("Baseline"); | |
333 | hbas[index2]->GetXaxis()->SetTitleSize(0.07); | |
334 | hbas[index2]->GetYaxis()->SetTitleSize(0.07); | |
335 | hbas[index2]->GetXaxis()->SetTitleOffset(0.6); | |
336 | hbas[index2]->GetYaxis()->SetTitleOffset(0.7); | |
337 | t3->DrawLatex(0.15,0.92,text); | |
338 | c1->Update(); | |
339 | ||
340 | ||
341 | c2->cd(index2+1); | |
342 | hrawn[index2]->SetMinimum(1.); | |
343 | hrawn[index2]->SetMaximum(6.); | |
344 | hrawn[index2]->Draw(); | |
345 | hrawn[index2]->GetXaxis()->SetTitle("Anode"); | |
346 | hrawn[index2]->GetYaxis()->SetTitle("Noise"); | |
347 | hrawn[index2]->GetXaxis()->SetTitleSize(0.07); | |
348 | hrawn[index2]->GetYaxis()->SetTitleSize(0.07); | |
349 | hrawn[index2]->GetXaxis()->SetTitleOffset(0.6); | |
350 | hrawn[index2]->GetYaxis()->SetTitleOffset(0.7); | |
351 | gStyle->SetOptStat(0); | |
352 | hrawn[index2]->SetStats(0); | |
353 | hcorrn[index2]->SetLineColor(2); | |
354 | hcorrn[index2]->Draw("SAME"); | |
355 | t1->Draw(); | |
356 | t2->Draw(); | |
357 | t3->DrawLatex(0.15,0.92,text); | |
358 | c2->Update(); | |
359 | ||
360 | c3->cd(index2+1); | |
361 | hdbas[index2]->Draw(); | |
362 | hdbas[index2]->GetXaxis()->SetTitle("Baseline"); | |
363 | hdbas[index2]->GetXaxis()->SetTitleSize(0.07); | |
364 | hdbas[index2]->GetXaxis()->SetTitleOffset(0.6); | |
365 | t3->DrawLatex(0.15,0.92,text); | |
366 | c3->Update(); | |
367 | ||
368 | c4->cd(index2+1); | |
369 | hdrawn[index2]->Draw(); | |
370 | hdrawn[index2]->GetXaxis()->SetTitle("Noise"); | |
371 | hdrawn[index2]->GetXaxis()->SetTitleSize(0.07); | |
372 | hdrawn[index2]->GetXaxis()->SetTitleOffset(0.6); | |
373 | hdcorrn[index2]->SetLineColor(2); | |
374 | hdcorrn[index2]->Draw("SAME"); | |
375 | t1->Draw(); | |
376 | t2->Draw(); | |
377 | t3->DrawLatex(0.15,0.92,text); | |
378 | c4->Update(); | |
379 | } | |
380 | } | |
381 | ||
382 | c1->SaveAs("Baselines.gif"); | |
383 | c2->SaveAs("Noise.gif"); | |
384 | c3->SaveAs("BaselinesDist.gif"); | |
385 | c4->SaveAs("NoiseDist.gif"); | |
386 | ||
387 | Char_t delfil[100]; | |
388 | for(Int_t iddl=0; iddl<kTotDDL;iddl++){ | |
389 | for(Int_t imod=0; imod<kModPerDDL;imod++){ | |
390 | for(Int_t isid=0;isid<kSides;isid++){ | |
391 | sprintf(delfil,"rm SDDbase_step1_ddl%02dc%02d_sid%d.data",iddl,imod,isid); | |
392 | gSystem->Exec(delfil); | |
393 | } | |
394 | } | |
395 | } | |
396 | ||
397 | } | |
398 | ||
e5b2f7f2 | 399 | void AnalyzeSDDNoiseAllMod(Int_t nrun, Int_t n2, Int_t year=2009, Char_t* dir="LHC09b_SDD", |
a8cd2ab2 | 400 | Int_t adcfreq=20, |
401 | Int_t nDDL=0, | |
e5b2f7f2 | 402 | Int_t firstEv=18, |
403 | Int_t lastEv=20){ | |
a8cd2ab2 | 404 | |
91282711 | 405 | TGrid::Connect("alien:",0,0,"t"); |
406 | Char_t filnam[200]; | |
e5b2f7f2 | 407 | sprintf(filnam,"alien:///alice/data/%d/%s/%09d/raw/%02d%09d%03d.10.root",year,dir,nrun,year-2000,nrun,n2); |
91282711 | 408 | printf("Open file %s\n",filnam); |
e5b2f7f2 | 409 | AnalyzeSDDNoiseAllMod(filnam,adcfreq,nDDL,firstEv,lastEv); |
91282711 | 410 | } |