1 #if !defined(__CINT__) || defined(__MAKECINT__)
12 #include "AliHMPIDDigit.h"
13 #include "AliHMPIDRawStream.h"
19 TH2F hgPedMapMean[14][6];
20 TH2F hgPedMapSigma[14][6];
21 TH1F hgPedMapMean1D[14][6];
22 TH1F hgPedMapSigma1D[14][6];
23 TH1F hgPedMapSigma1Db[14][6];
34 //++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
35 void Convert(Int_t ddl,Int_t r,Int_t d,Int_t a,Int_t &ch, Int_t &pc, Int_t &px, Int_t &py)
39 Int_t a2y[6]={3,2,4,1,5,0};//pady for a given address (for single DILOGIC chip)
42 Int_t tmp=(24-r)/8; Int_t pc=(ddl%2)? 5-2*tmp:2*tmp;
44 tmp=(ddl%2)?r-1:(24-r); Int_t py=6*(tmp%8)+a2y[a%6];
46 //++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
47 void PlotPedHisto(Char_t *tag,Int_t runNum, Int_t ddl=0,Int_t row=1,Int_t dil=1,Int_t pad=0)
50 if( 0<=ddl && ddl<=13 && 1<=dil && dil<=10 && 1<=row && row<=24 && 0<=pad && pad<=47)
55 hgtmp=(TH1I*)(fgin[ddl]->Get(Form("hDDL_%d_Row_%d_Dil_%d_Pad_%d",ddl,row,dil,pad)))->Clone();
56 hgtmp2=(TH1I*)(fgin[ddl]->Get(Form("hDDL_%d_Row_%d_Dil_%d_Pad_%d",ddl,row,dil,pad)))->Clone();
58 TCanvas *c1=new TCanvas(Form("hDDL%s_%d_Row_%d_Dil_%d_Pad_%d",tag,ddl,row,dil,pad),Form("hDDL%s_%d_Row_%d_Dil_%d_Pad_%d",tag,ddl,row,dil,pad));
60 hgtmp->SetXTitle("ADC");
61 hgtmp->SetYTitle("Entries");
63 hgtmp->SetAxisRange(0,350);
64 hgtmp->Draw("hist same");
65 hgtmp2->SetFillColor(5);
66 hgtmp2->Draw("hist same");
69 c1->SaveAs(Form("hDDL%s_%d_Row_%d_Dil_%d_Pad_%d.eps",tag,ddl,row,dil,pad));
70 c1->SaveAs(Form("hDDL%s_%d_Row_%d_Dil_%d_Pad_%d.gif",tag,ddl,row,dil,pad));
84 if( ddl < 0 || ddl > 13) {Printf("Not a valid DDL, exiting %d ...",ddl); }
85 if( dil < 1 || dil > 10) {Printf("Not a valid DIL, exiting %d ...",dil); }
86 if( row < 1 || row > 24) {Printf("Not a valid ROW, exiting %d ...",row); }
87 if( pad < 0 || pad > 48) {Printf("Not a valid PAD, exiting %d ...",pad); }
90 //++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
91 void PlotHisto1D(TH1& hin,Char_t* name)
93 TCanvas *c1=new TCanvas("c1","c1");
95 c1->SaveAs(Form("%s.eps",name));
96 c1->SaveAs(Form("%s.gif",name));
98 //++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
99 void PlotHisto2D(TH2 &hin,Char_t* name,Char_t* opt)
101 TCanvas *c1=new TCanvas("c1","c1");
103 c1->SaveAs(Form("%s.eps",name));
104 c1->SaveAs(Form("%s.gif",name));
106 //++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
107 void ProcPed(Int_t nDDL,Int_t mode)
110 if(gSystem->IsFileInIncludePath(Form("./HmpidPedDdl%02i.txt",nDDL)))
111 ifstream infile(Form("./HmpidPedDdl%02i.txt",nDDL));
114 for(Int_t ni=0;ni<6;ni++)
116 hgPedMapMean[nDDL][ni] = new TH2F(Form("hPedMapMean_DDL%d_PC_%d" ,nDDL,ni), Form("Pedestal: Mean, DDL(0-13): %d PC(0-5): %d;padx;pady" ,nDDL,ni),80,0,80,48,0,48);
117 hgPedMapSigma[nDDL][ni] = new TH2F(Form("hPedMapSigma_DDL%d_PC_%d",nDDL,ni), Form("Pedestal: Sigma, DDL(0-13): %d PC(0-5): %d;padx;pady",nDDL,ni),80,0,80,48,0,48);
118 hgPedMapMean1D[nDDL][ni] = new TH1F(Form("hPedMapMean1D_DDL%d_PC_%d" ,nDDL,ni), Form("Pedestal: Mean, DDL(0-13): %d PC(0-5): %d;pad;Mean" ,nDDL,ni),3841,-0.5,3840.5);
119 hgPedMapSigma1D[nDDL][ni] = new TH1F(Form("hPedMapSigma1D_DDL%d_PC_%d" ,nDDL,ni), Form("Pedestal: Sigma, DDL(0-13): %d PC(0-5): %d;pad;Sigma" ,nDDL,ni),3841,-0.5,3840.5);
120 hgPedMapSigma1Db[nDDL][ni] = new TH1F(Form("hPedMapSigma1Db_DDL%d_PC_%d" ,nDDL,ni), Form("Pedestal: Sigma, DDL(0-13): %d PC(0-5): %d;pad;Sigma" ,nDDL,ni),3841,-0.5,3840.5);
121 hgPedMapSigma1Db[nDDL][ni]->SetMaximum(2.5);
124 Int_t nSigCut,r,d,a,hard; Float_t mean,sigma;
125 Int_t ch=0,pc=0,px=0,py=0;
127 Int_t runNumber=-99999;
132 Int_t nBadEv;Float_t nBadEvPer;
133 Printf("Start reading DDL: %d ...",nDDL);
134 infile>>tName>>runNumber;fgRunNum=runNumber;
135 infile>>tName>>ldcId;
136 infile>>tName>>timeStamp;
138 infile>>tName>>nDdlEv;
139 infile>>tName>>nBadEv;
140 infile>>tName>>nBadEvPer;
141 infile>>tName>>nSigCut;
143 Printf("RunNumber: %d, LdcId: %d TimeStamp: %d nEv: %d",fgRunNum,ldcId,timeStamp,nEv);
145 while(!infile.eof()){
146 infile>>dec>>r>>d>>a>>mean>>sigma>>hex>>hard;
150 Convert(nDDL,r,d,a,ch,pc,px,py);
151 hgPedMapMean[nDDL][pc].Fill(px,py,mean);
152 hgPedMapSigma[nDDL][pc].Fill(px,py,sigma);
153 hgPedMapMean1D[nDDL][pc].Fill(cnt%3840,mean);
154 hgPedMapSigma1D[nDDL][pc].Fill(cnt%3840,sigma);
155 hgPedMapSigma1Db[nDDL][pc].Fill(cnt%3840,sigma);
158 if(sigma>3 && sigma!=1000) PlotPedHisto("Sigma",runNumber,nDDL,r,d,a);
161 if(sigma==1000) PlotPedHisto("NoEntry",runNumber,nDDL,r,d,a);
169 Printf("Stop reading DDL: %d ...",nDDL);
171 /* fill the overall histos */
175 for(Int_t npc=0;npc<6;npc++)
177 /* plot pedestal mean values */
178 if(hgPedMapMean[nDDL][npc].GetEntries()!=0) {
179 hgPedMapMean[nDDL][npc].SetStats(0);
180 PlotHisto2D(hgPedMapMean[nDDL][npc],Form("Run%d_PedMeanA_DDL%d_PC%d",fgRunNum,nDDL,npc),"colz");
181 PlotHisto2D(hgPedMapMean[nDDL][npc],Form("Run%d_PedMeanB_DDL%d_PC%d",fgRunNum,nDDL,npc),"surf1");
183 /* plot pedestal sigma values */
184 if(hgPedMapSigma[nDDL][npc].GetEntries()!=0) {
185 hgPedMapSigma[nDDL][npc].SetStats(0);
186 PlotHisto2D(hgPedMapSigma[nDDL][npc],Form("Run%d_PedSigmaA_DDL%d_PC%d",fgRunNum,nDDL,npc),"colz");
187 PlotHisto2D(hgPedMapSigma[nDDL][npc],Form("Run%d_PedSigmaB_DDL%d_PC%d",fgRunNum,nDDL,npc),"surf1");
189 /* plot pedestal sigma values */
190 if(hgPedMapSigma[nDDL][npc].GetEntries()!=0) {
191 hgPedMapSigma[nDDL][npc].SetStats(0);
192 hgPedMapSigma[nDDL][npc].SetMaximum(2.5);
193 PlotHisto2D(hgPedMapSigma[nDDL][npc],Form("Run%d_PedSigmaC_DDL%d_PC%d",fgRunNum,nDDL,npc),"colz");
194 PlotHisto2D(hgPedMapSigma[nDDL][npc],Form("Run%d_PedSigmaD_DDL%d_PC%d",fgRunNum,nDDL,npc),"surf1");
196 /* plot pedestal mean values */
197 if(hgPedMapMean1D[nDDL][npc].GetEntries()!=0) {
198 hgPedMapMean1D[nDDL][npc].SetStats(0);
199 PlotHisto1D(hgPedMapMean1D[nDDL][npc],Form("Run%d_PedMean1DA_DDL%d_PC%d",fgRunNum,nDDL,npc));
201 /* plot pedestal sigma values */
202 if(hgPedMapSigma1D[nDDL][npc].GetEntries()!=0) {
203 hgPedMapSigma1D[nDDL][npc].SetStats(0);
204 PlotHisto1D(hgPedMapSigma1D[nDDL][npc],Form("Run%d_PedSigma1DA_DDL%d_PC%d",fgRunNum,nDDL,npc));
206 /* plot pedestal sigma values */
207 if(hgPedMapSigma1Db[nDDL][npc].GetEntries()!=0) {
208 hgPedMapSigma1Db[nDDL][npc].SetStats(0);
209 PlotHisto1D(hgPedMapSigma1D[nDDL][npc],Form("Run%d_PedSigma1DA_DDL%d_PC%d",fgRunNum,nDDL,npc));
216 //++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
217 void ProcErr(Int_t nDDL,Int_t mode)
219 if(gSystem->IsFileInIncludePath(Form("./HmpidErrorsDdl%02i.txt",nDDL)))
220 ifstream infile(Form("./HmpidErrorsDdl%02i.txt",nDDL));
224 hgDdlErr[nDDL] = new TH1I(Form("hPedErr_DDL%d",nDDL),Form("DDL(%d) Decoding Errors",nDDL),AliHMPIDRawStream::kSumErr+1,-0.5,AliHMPIDRawStream::kSumErr+0.5);
225 for(Int_t ilabel=0; ilabel< AliHMPIDRawStream::kSumErr; ilabel++) {
226 hgDdlErr[nDDL].SetStats(0);
227 hgDdlErr[nDDL].GetXaxis()->CenterLabels(kTRUE);
228 hgDdlErr[nDDL].GetXaxis()->SetBinLabel((ilabel+1),Form("%i %s",ilabel+1,AliHMPIDRawStream::GetErrName(ilabel)));
229 hgDdlErr[nDDL].SetYTitle("Error #");
231 hgDdlErr[nDDL].SetFillColor(5);
234 Int_t runNumber=-99999;
239 Int_t rerr;Int_t nBadEv;Float_t nBadEvPer;
240 Int_t row,dil,pad,nzero;
242 Printf("Start reading Error File DDL: %d ...",nDDL);
243 infile>>tName>>runNumber;fgRunNum=runNumber;
244 infile>>tName>>ldcId;
245 infile>>tName>>timeStamp;
247 infile>>tName>>nDdlEv;
248 infile>>tName>>nBadEv;
249 infile>>tName>>nBadEvPer;
250 Printf("RunNumber: %d, LdcId: %d TimeStamp: %d nEv: %d",runNumber,ldcId,timeStamp,nEv);
251 for(Int_t ierr=0;ierr<AliHMPIDRawStream::kSumErr;ierr++)
253 infile>>rerr;hgDdlErr[nDDL].SetBinContent(ierr+1,rerr);
259 TCanvas *cerr = new TCanvas("cped","cped");
261 hgDdlErr[nDDL]->Draw();
262 cerr->SaveAs(Form("PedError_DDL%d.eps",nDDL));
263 cerr->SaveAs(Form("PedError_DDL%d.gif",nDDL));
268 while(!infile.eof()){
269 infile>>dec>>row>>dil>>pad>>nzero;
270 if( 0<=nDDL && nDDL<=13 && 1<=dil && dil<=10 && 1<=row && row<=24 && 0<=pad && pad<=47)
272 PlotPedHisto("ZeroQ",runNumber,nDDL,row,dil,pad);
283 //++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
284 void HplotDA(Int_t runNumber)
286 gStyle->SetPalette(1);
289 Printf("********************************************************************");
290 Printf("********************** HplotDA *************************************");
292 Printf("*** Choose from the following options: ***");
294 Printf("*** 1.) Draw errors for all DDL ***");
295 Printf("*** 2.) Draw pedestals for all DDL ***");
296 Printf("*** 3.) Draw ADC histo for pads with sigma > 3 ***");
297 Printf("*** 4.) Draw ADC histo for pads with zero charge ***");
298 Printf("*** 5.) Draw ADC histo for pads with no valid charge ***");
299 Printf("*** 6.) Run all ***");
301 Printf("********************************************************************");
302 Printf("*** Please select: ");
308 for(Int_t i=0;i<14;i++)
319 fgin[i]=new TFile(Form("Run%d_DDL%d.root",runNumber,i),"read");
324 fgin[i]=new TFile(Form("Run%d_DDL%d.root",runNumber,i),"read");
329 fgin[i]=new TFile(Form("Run%d_DDL%d.root",runNumber,i),"read");
334 if(i==0) fgout=new TFile(Form("SummaryOfRun%d.root",fgRunNum),"recreate");
338 fgin[i]=new TFile(Form("Run%d_DDL%d.root",runNumber,i),"read");
346 for(Int_t ipc=0;ipc<6;ipc++)
348 hgPedMapMean[i][ipc].Write();
349 hgPedMapSigma[i][ipc].Write();
350 hgPedMapMean1D[i][ipc].Write();
351 hgPedMapSigma1D[i][ipc].Write();
352 hgPedMapSigma1Db[i][ipc].Write();
358 Printf("Not a valid selection bye-bye....");
363 if(fgout!=0x0)fgout->Close();
366 //++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++