1 #if !defined(__CINT__) || defined(__MAKECINT__)
4 #include <TStopwatch.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"
16 #include "AliITSRawStreamSDDCompressed.h"
17 #include "TPaveStats.h"
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)
27 void AnalyzeSDDNoiseAllMod(Char_t *datafil,
34 const Int_t kTotDDL=24;
35 const Int_t kModPerDDL=12;
37 Bool_t writtenoutput=kFALSE;
39 AliITSOnlineSDDBase **base=new AliITSOnlineSDDBase*[kTotDDL*kModPerDDL*kSides];
40 TH2F **histo=new TH2F*[kTotDDL*kModPerDDL*kSides];
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;
47 base[index]=new AliITSOnlineSDDBase(iddl,imod,isid);
48 if(adcfreq==40) base[index]->SetLastGoodTB(254);
49 else base[index]->SetLastGoodTB(126);
50 sprintf(hisnam,"h%02dc%02ds%d",iddl,imod,isid);
51 histo[index]=new TH2F(hisnam,"",256,-0.5,255.5,256,-0.5,255.5);
57 TCanvas* c0 = new TCanvas("c0","Ev Display",900,900);
58 gStyle->SetPalette(1);
63 if(strstr(datafil,".root")!=0){
64 rd=new AliRawReaderRoot(datafil,iev);
66 rd=new AliRawReaderDate(datafil,iev);
68 TLatex *t0=new TLatex();
70 t0->SetTextSize(0.06);
75 c0->Divide(4,6,0.001,0.001);
76 printf("Event # %d ",iev);
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();
87 UChar_t cdhAttr=AliITSRawStreamSDD::ReadBlockAttributes(rd);
88 UInt_t amSamplFreq=AliITSRawStreamSDD::ReadAMSamplFreqFromCDH(cdhAttr);
89 AliITSRawStream* s=AliITSRawStreamSDD::CreateRawStreamSDD(rd,cdhAttr);
91 printf("Use %s raw stream, sampling frequency %d MHz\n",s->ClassName(),amSamplFreq);
96 Int_t iDDL=rd->GetDDLID();
97 Int_t iCarlos=s->GetCarlosId();
98 if(s->IsCompletedModule()) continue;
99 if(s->IsCompletedDDL()) continue;
100 if(iDDL>=0 && iDDL<kTotDDL){
101 Int_t index=kSides*(kModPerDDL*iDDL+iCarlos)+s->GetChannel();
102 histo[index]->Fill(s->GetCoord2(),s->GetCoord1(),s->GetSignal());
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]);
113 Int_t index2=kSides*imod+isid;
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);
124 }while(rd->NextEvent()&&iev<=lastEv);
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!!!!!!!!!
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];
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);
148 if(adcfreq==40) corr[index]->SetLastGoodTB(254);
149 else corr[index]->SetLastGoodTB(126);
157 if(strstr(datafil,".root")!=0){
158 rd2=new AliRawReaderRoot(datafil,iev);
160 rd2=new AliRawReaderDate(datafil,iev);
164 c0->Divide(4,6,0.001,0.001);
165 printf("Event # %d ",iev);
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();
176 UChar_t cdhAttr=AliITSRawStreamSDD::ReadBlockAttributes(rd2);
177 UInt_t amSamplFreq=AliITSRawStreamSDD::ReadAMSamplFreqFromCDH(cdhAttr);
178 AliITSRawStream* s=AliITSRawStreamSDD::CreateRawStreamSDD(rd2,cdhAttr);
180 printf("Use %s raw stream, sampling frequency %d MHz\n",s->ClassName(),amSamplFreq);
184 Int_t iDDL=rd2->GetDDLID();
185 Int_t iCarlos=s->GetCarlosId();
186 if(s->IsCompletedModule()) continue;
187 if(s->IsCompletedDDL()) continue;
188 if(iDDL>=0 && iDDL<kTotDDL){
189 Int_t index=kSides*(kModPerDDL*iDDL+iCarlos)+s->GetChannel();
190 histo[index]->Fill(s->GetCoord2(),s->GetCoord1(),s->GetSignal());
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]);
202 Int_t index2=kSides*imod+isid;
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);
213 }while(rd2->NextEvent()&&iev<=lastEv);
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);
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;
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);
238 htotbaseq->Fill(basleq);
239 htotnoise->Fill(noi);
240 htotnoisecorr->Fill(cornoi);
247 // Draw Statistics of baselines and noise
248 TCanvas *call=new TCanvas("call","General stats",700,700);
252 htotbas->GetXaxis()->SetTitle("Baselines");
253 htotbas->GetXaxis()->SetTitleSize(0.07);
254 htotbas->GetXaxis()->SetTitleOffset(0.6);
257 htotbaseq->GetXaxis()->SetTitle("Baselines after equalization");
258 htotbaseq->GetXaxis()->SetTitleSize(0.07);
259 htotbaseq->GetXaxis()->SetTitleOffset(0.6);
261 htotnoisecorr->SetLineColor(2);
262 htotnoisecorr->Draw();
264 TPaveStats *st1=(TPaveStats*)htotnoisecorr->GetListOfFunctions()->FindObject("stats");
267 htotnoisecorr->GetXaxis()->SetTitle("Noise");
268 htotnoisecorr->GetXaxis()->SetTitleSize(0.07);
269 htotnoisecorr->GetXaxis()->SetTitleOffset(0.6);
270 htotnoise->Draw("SAMES");
272 TPaveStats *st2=(TPaveStats*)htotnoise->GetListOfFunctions()->FindObject("stats");
278 hstatus->GetXaxis()->SetTitle("Anode Status (0=bad 1=good)");
279 hstatus->GetXaxis()->SetTitleSize(0.07);
280 hstatus->GetXaxis()->SetTitleOffset(0.6);
282 call->SaveAs("GenStatsPedestal.gif");
284 // Draw baselines and noise for all modules of the selected DDL
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");
306 t1->SetTextSize(0.05);
307 TLatex *t2=new TLatex(0.4,0.2,"Corrected Noise");
309 t2->SetTextSize(0.05);
311 TLatex *t3=new TLatex();
313 t3->SetTextSize(0.06);
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();
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);
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");
357 t3->DrawLatex(0.15,0.92,text);
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);
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");
377 t3->DrawLatex(0.15,0.92,text);
382 c1->SaveAs("Baselines.gif");
383 c2->SaveAs("Noise.gif");
384 c3->SaveAs("BaselinesDist.gif");
385 c4->SaveAs("NoiseDist.gif");
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);
399 void AnalyzeSDDNoiseAllMod(Int_t nrun, Int_t n2, Int_t year=2009, Char_t* dir="LHC09b_SDD",
405 TGrid::Connect("alien:",0,0,"t");
407 sprintf(filnam,"alien:///alice/data/%d/%s/%09d/raw/%02d%09d%03d.10.root",year,dir,nrun,year-2000,nrun,n2);
408 printf("Open file %s\n",filnam);
409 AnalyzeSDDNoiseAllMod(filnam,adcfreq,nDDL,firstEv,lastEv);