2 |
* FTrkCalibQLookExpert.cxx |
* FTrkCalibQLookExpert.cxx |
3 |
* |
* |
4 |
* autor: D.Fedele |
* autor: D.Fedele |
5 |
* version v1r08 |
* version v1r28 |
6 |
* Parameters: |
* Parameters: |
7 |
* file - the data file to analyze |
* file - the data file to analyze |
8 |
* step - select =1 in order to analyze one event at time |
* step - select =1 in order to analyze one event at time |
77 |
Int_t nevents=0; |
Int_t nevents=0; |
78 |
Int_t minevent = 0; |
Int_t minevent = 0; |
79 |
Int_t maxevent = 0; |
Int_t maxevent = 0; |
80 |
ULong64_t FOBT[2]; |
ULong_t FOBT[2]; |
81 |
|
|
82 |
FOBT[0]=0; |
FOBT[0]=0; |
83 |
FOBT[1]=0; |
FOBT[1]=0; |
176 |
} else if (toevent > nevents) { |
} else if (toevent > nevents) { |
177 |
maxevent = nevents; |
maxevent = nevents; |
178 |
} else { |
} else { |
179 |
maxevent = fromevent; |
maxevent = nevents; |
180 |
}; |
}; |
181 |
}; |
}; |
182 |
|
|
187 |
//********************************************************************** |
//********************************************************************** |
188 |
|
|
189 |
Int_t hcevent=hotr->GetEntries(); |
Int_t hcevent=hotr->GetEntries(); |
190 |
ULong64_t HOBT[hcevent], TOBT[hcevent]; |
Int_t tcevent=totr->GetEntries(); |
191 |
|
ULong_t HOBT[hcevent], TOBT[tcevent]; |
192 |
for (Int_t i = 0; i < hcevent; i++){ |
for (Int_t i = 0; i < hcevent; i++){ |
|
totr->GetEntry(i); |
|
193 |
hotr->GetEntry(i); |
hotr->GetEntry(i); |
194 |
ph = eh4->GetPscuHeader(); |
ph = eh4->GetPscuHeader(); |
195 |
HOBT[i]= ph->GetOrbitalTime(); |
HOBT[i]= ph->GetOrbitalTime(); |
196 |
|
} |
197 |
|
for (Int_t i = 0; i < tcevent; i++){ |
198 |
|
totr->GetEntry(i); |
199 |
ph = eh3->GetPscuHeader(); |
ph = eh3->GetPscuHeader(); |
200 |
TOBT[i]= ph->GetOrbitalTime(); |
TOBT[i]= ph->GetOrbitalTime(); |
201 |
} |
} |
212 |
ph = eh1->GetPscuHeader(); |
ph = eh1->GetPscuHeader(); |
213 |
cod = eh1->GetCounter(); |
cod = eh1->GetCounter(); |
214 |
FOBT[0]= ph->GetOrbitalTime(); |
FOBT[0]= ph->GetOrbitalTime(); |
215 |
ctrk.daqmode[trk1->DSPnumber[m]-1]=trk1->DAQmode[m]; |
if(trk1->DSPnumber[m]>0){ |
216 |
ctrk.dspnum[trk1->DSPnumber[m]-1]=trk1->DSPnumber[m]; |
if(trk1->DSPnumber[m]<13){ |
217 |
ctrk.calibnum[trk1->DSPnumber[m]-1]=trk1->calibnumber[m]; |
ctrk.daqmode[trk1->DSPnumber[m]-1]=trk1->DAQmode[m]; |
218 |
ctrk.ncalev[trk1->DSPnumber[m]-1]=trk1->ncalib_event[m]; |
ctrk.dspnum[trk1->DSPnumber[m]-1]=trk1->DSPnumber[m]; |
219 |
ctrk.ped1[trk1->DSPnumber[m]-1]=trk1->ped_l1[m]; |
ctrk.calibnum[trk1->DSPnumber[m]-1]=trk1->calibnumber[m]; |
220 |
ctrk.ped2[trk1->DSPnumber[m]-1]=trk1->ped_l2[m]; |
ctrk.ncalev[trk1->DSPnumber[m]-1]=trk1->ncalib_event[m]; |
221 |
ctrk.ped3[trk1->DSPnumber[m]-1]=trk1->ped_l3[m]; |
ctrk.ped1[trk1->DSPnumber[m]-1]=trk1->ped_l1[m]; |
222 |
ctrk.sig1[trk1->DSPnumber[m]-1]=trk1->sig_l1[m]; |
ctrk.ped2[trk1->DSPnumber[m]-1]=trk1->ped_l2[m]; |
223 |
ctrk.sig2[trk1->DSPnumber[m]-1]=trk1->sig_l2[m]; |
ctrk.ped3[trk1->DSPnumber[m]-1]=trk1->ped_l3[m]; |
224 |
ctrk.sig3[trk1->DSPnumber[m]-1]=trk1->sig_l3[m]; |
ctrk.sig1[trk1->DSPnumber[m]-1]=trk1->sig_l1[m]; |
225 |
ctrk.nbad1[trk1->DSPnumber[m]-1]=trk1->nbad_l1[m]; |
ctrk.sig2[trk1->DSPnumber[m]-1]=trk1->sig_l2[m]; |
226 |
ctrk.nbad2[trk1->DSPnumber[m]-1]=trk1->nbad_l2[m]; |
ctrk.sig3[trk1->DSPnumber[m]-1]=trk1->sig_l3[m]; |
227 |
ctrk.nbad3[trk1->DSPnumber[m]-1]=trk1->nbad_l3[m]; |
ctrk.nbad1[trk1->DSPnumber[m]-1]=trk1->nbad_l1[m]; |
228 |
ctrk.calfl[trk1->DSPnumber[m]-1]=trk1->cal_flag[m]; |
ctrk.nbad2[trk1->DSPnumber[m]-1]=trk1->nbad_l2[m]; |
229 |
ctrk.crc_c[trk1->DSPnumber[m]-1][0]=trk1->crc_cal[m][0]; |
ctrk.nbad3[trk1->DSPnumber[m]-1]=trk1->nbad_l3[m]; |
230 |
ctrk.crc_c[trk1->DSPnumber[m]-1][1]=trk1->crc_cal[m][1]; |
ctrk.calfl[trk1->DSPnumber[m]-1]=trk1->cal_flag[m]; |
231 |
ctrk.crc_c[trk1->DSPnumber[m]-1][2]=trk1->crc_cal[m][2]; |
ctrk.crc_c[trk1->DSPnumber[m]-1][0]=trk1->crc_cal[m][0]; |
232 |
ctrk.crc_hc[trk1->DSPnumber[m]-1]=trk1->crc_hcal[m]; |
ctrk.crc_c[trk1->DSPnumber[m]-1][1]=trk1->crc_cal[m][1]; |
233 |
for (Int_t j = 0; j < 3072; j++){ |
ctrk.crc_c[trk1->DSPnumber[m]-1][2]=trk1->crc_cal[m][2]; |
234 |
ctrk.dspped[trk1->DSPnumber[m]-1][j]=trk1->DSPped_par[m][j]; |
ctrk.crc_hc[trk1->DSPnumber[m]-1]=trk1->crc_hcal[m]; |
235 |
ctrk.dspsig[trk1->DSPnumber[m]-1][j]=trk1->DSPsig_par[m][j]; |
for (Int_t j = 0; j < 3072; j++){ |
236 |
ctrk.dspbad[trk1->DSPnumber[m]-1][j]=trk1->DSPbad_par[m][j]; |
ctrk.dspped[trk1->DSPnumber[m]-1][j]=trk1->DSPped_par[m][j]; |
237 |
}; |
ctrk.dspsig[trk1->DSPnumber[m]-1][j]=trk1->DSPsig_par[m][j]; |
238 |
|
ctrk.dspbad[trk1->DSPnumber[m]-1][j]=trk1->DSPbad_par[m][j]; |
239 |
|
}; |
240 |
|
} |
241 |
|
} |
242 |
ph = eh2->GetPscuHeader(); |
ph = eh2->GetPscuHeader(); |
243 |
FOBT[1]= ph->GetOrbitalTime(); |
FOBT[1]= ph->GetOrbitalTime(); |
244 |
ctrk.daqmode[trk2->DSPnumber[m]-1]=trk2->DAQmode[m]; |
if(trk2->DSPnumber[m]>0){ |
245 |
ctrk.dspnum[trk2->DSPnumber[m]-1]=trk2->DSPnumber[m]; |
if(trk2->DSPnumber[m]<13){ |
246 |
ctrk.calibnum[trk2->DSPnumber[m]-1]=trk2->calibnumber[m]; |
ctrk.daqmode[trk2->DSPnumber[m]-1]=trk2->DAQmode[m]; |
247 |
ctrk.ncalev[trk2->DSPnumber[m]-1]=trk2->ncalib_event[m]; |
ctrk.dspnum[trk2->DSPnumber[m]-1]=trk2->DSPnumber[m]; |
248 |
ctrk.ped1[trk2->DSPnumber[m]-1]=trk2->ped_l1[m]; |
ctrk.calibnum[trk2->DSPnumber[m]-1]=trk2->calibnumber[m]; |
249 |
ctrk.ped2[trk2->DSPnumber[m]-1]=trk2->ped_l2[m]; |
ctrk.ncalev[trk2->DSPnumber[m]-1]=trk2->ncalib_event[m]; |
250 |
ctrk.ped3[trk2->DSPnumber[m]-1]=trk2->ped_l3[m]; |
ctrk.ped1[trk2->DSPnumber[m]-1]=trk2->ped_l1[m]; |
251 |
ctrk.sig1[trk2->DSPnumber[m]-1]=trk2->sig_l1[m]; |
ctrk.ped2[trk2->DSPnumber[m]-1]=trk2->ped_l2[m]; |
252 |
ctrk.sig2[trk2->DSPnumber[m]-1]=trk2->sig_l2[m]; |
ctrk.ped3[trk2->DSPnumber[m]-1]=trk2->ped_l3[m]; |
253 |
ctrk.sig3[trk2->DSPnumber[m]-1]=trk2->sig_l3[m]; |
ctrk.sig1[trk2->DSPnumber[m]-1]=trk2->sig_l1[m]; |
254 |
ctrk.nbad1[trk2->DSPnumber[m]-1]=trk2->nbad_l1[m]; |
ctrk.sig2[trk2->DSPnumber[m]-1]=trk2->sig_l2[m]; |
255 |
ctrk.nbad2[trk2->DSPnumber[m]-1]=trk2->nbad_l2[m]; |
ctrk.sig3[trk2->DSPnumber[m]-1]=trk2->sig_l3[m]; |
256 |
ctrk.nbad3[trk2->DSPnumber[m]-1]=trk2->nbad_l3[m]; |
ctrk.nbad1[trk2->DSPnumber[m]-1]=trk2->nbad_l1[m]; |
257 |
ctrk.calfl[trk2->DSPnumber[m]-1]=trk2->cal_flag[m]; |
ctrk.nbad2[trk2->DSPnumber[m]-1]=trk2->nbad_l2[m]; |
258 |
ctrk.crc_c[trk1->DSPnumber[m]-1][0]=trk2->crc_cal[m][0]; |
ctrk.nbad3[trk2->DSPnumber[m]-1]=trk2->nbad_l3[m]; |
259 |
ctrk.crc_c[trk1->DSPnumber[m]-1][1]=trk2->crc_cal[m][1]; |
ctrk.calfl[trk2->DSPnumber[m]-1]=trk2->cal_flag[m]; |
260 |
ctrk.crc_c[trk1->DSPnumber[m]-1][2]=trk2->crc_cal[m][2]; |
ctrk.crc_c[trk1->DSPnumber[m]-1][0]=trk2->crc_cal[m][0]; |
261 |
ctrk.crc_hc[trk1->DSPnumber[m]-1]=trk2->crc_hcal[m]; |
ctrk.crc_c[trk1->DSPnumber[m]-1][1]=trk2->crc_cal[m][1]; |
262 |
for (Int_t j = 0; j < 3072; j++){ |
ctrk.crc_c[trk1->DSPnumber[m]-1][2]=trk2->crc_cal[m][2]; |
263 |
ctrk.dspped[trk2->DSPnumber[m]-1][j]=trk2->DSPped_par[m][j]; |
ctrk.crc_hc[trk1->DSPnumber[m]-1]=trk2->crc_hcal[m]; |
264 |
ctrk.dspsig[trk2->DSPnumber[m]-1][j]=trk2->DSPsig_par[m][j]; |
for (Int_t j = 0; j < 3072; j++){ |
265 |
ctrk.dspbad[trk2->DSPnumber[m]-1][j]=trk2->DSPbad_par[m][j]; |
ctrk.dspped[trk2->DSPnumber[m]-1][j]=trk2->DSPped_par[m][j]; |
266 |
|
ctrk.dspsig[trk2->DSPnumber[m]-1][j]=trk2->DSPsig_par[m][j]; |
267 |
|
ctrk.dspbad[trk2->DSPnumber[m]-1][j]=trk2->DSPbad_par[m][j]; |
268 |
|
} |
269 |
|
} |
270 |
} |
} |
271 |
} |
} |
272 |
|
|
273 |
|
|
274 |
// |
// |
275 |
// other variables definitions |
// other variables definitions |
276 |
|
|
277 |
Int_t risposta=0; |
Int_t risposta=0; |
278 |
stringstream fromfile; |
stringstream fromfile,isfile; |
279 |
|
|
280 |
fromfile<<"FTrkCalibQLook_EXPERT File: "<<ffile<<" -- CalibHeader OBT= "<<HOBT[(cod->Get(pctp->CalibHeader))-1]<<" -- Calib pkt OBT= "<<FOBT[0]<<" -- CalibTrailer OBT= "<<TOBT[(cod->Get(pctp->CalibTrailer))-1]<<" --"; |
fromfile<<"FTrkCalibQLook_EXPERT File: "<<ffile<<" -- CalibHeader OBT= "<<HOBT[(cod->Get(pctp->CalibHeader))-1]<<" -- Calib pkt OBT= "<<FOBT[0]<<" -- CalibTrailer OBT= "<<TOBT[(cod->Get(pctp->CalibTrailer))]<<" --"; |
281 |
|
|
282 |
gStyle->SetLabelSize(0.07,"x"); |
gStyle->SetLabelSize(0.07,"x"); |
283 |
gStyle->SetLabelSize(0.07,"y"); |
gStyle->SetLabelSize(0.07,"y"); |
305 |
Int_t canvasy=900; |
Int_t canvasy=900; |
306 |
TCanvas *c1 = new TCanvas("c1","FTrkCalibQLook_EXPERT_ped",canvasx,canvasy); |
TCanvas *c1 = new TCanvas("c1","FTrkCalibQLook_EXPERT_ped",canvasx,canvasy); |
307 |
c1->SetFillColor(10); |
c1->SetFillColor(10); |
308 |
tzz->DrawLatex(.01,0.98,fromfile.str().c_str()); |
tzz->DrawLatex(.01,.98,fromfile.str().c_str()); |
309 |
tzz->DrawLatex(.90,0.98,"PEDESTAL"); |
tzz->DrawLatex(.90,.98,"PEDESTAL"); |
310 |
|
|
311 |
|
isfile<<"Chips with white line at least once showed anomalous behaviour"; |
312 |
|
tzz->SetTextColor(17); |
313 |
|
tzz->SetTextSize(0.018); |
314 |
|
tzz->DrawLatex(.01,.96,isfile.str().c_str()); |
315 |
|
isfile.str(""); |
316 |
|
|
317 |
|
isfile<<"Boxes so colored point out anomalous chips (not necessarily broken)"; |
318 |
|
tzz->SetTextColor(107); |
319 |
|
// tzz->DrawLatex(.05,.96,isfile.str().c_str()); |
320 |
|
isfile.str(""); |
321 |
|
|
322 |
TCanvas *c2 = new TCanvas("c2","FTrkCalibQLook_EXPERT_sig",canvasx,canvasy); |
TCanvas *c2 = new TCanvas("c2","FTrkCalibQLook_EXPERT_sig",canvasx,canvasy); |
323 |
c2->SetFillColor(10); |
c2->SetFillColor(10); |
324 |
tzz->DrawLatex(.01,0.98,fromfile.str().c_str()); |
tzz->SetTextColor(1); |
325 |
tzz->DrawLatex(.90,0.98,"SIGMA"); |
tzz->DrawLatex(.01,.98,fromfile.str().c_str()); |
326 |
|
tzz->DrawLatex(.90,.98,"SIGMA"); |
327 |
|
|
328 |
|
isfile<<"Chips with white line at least once showed anomalous behaviour"; |
329 |
|
tzz->SetTextColor(17); |
330 |
|
tzz->SetTextSize(0.018); |
331 |
|
tzz->DrawLatex(.01,.96,isfile.str().c_str()); |
332 |
|
isfile.str(""); |
333 |
|
|
334 |
|
isfile<<"Boxes so colored point out anomalous chips (not necessarily broken)"; |
335 |
|
tzz->SetTextColor(107); |
336 |
|
// tzz->DrawLatex(.05,.96,isfile.str().c_str()); |
337 |
|
isfile.str(""); |
338 |
|
|
339 |
|
|
340 |
TCanvas *sig=new TCanvas("sig","FTrkCalibQLook_EXPERT_histosig",canvasx,canvasy); |
TCanvas *sig=new TCanvas("sig","FTrkCalibQLook_EXPERT_histosig",canvasx,canvasy); |
341 |
sig->SetFillColor(10); |
sig->SetFillColor(10); |
342 |
tzz->DrawLatex(.01,0.98,fromfile.str().c_str()); |
tzz->SetTextColor(1); |
343 |
tzz->DrawLatex(.85,0.97,"Histograms of the sigmas"); |
tzz->DrawLatex(.01,.98,fromfile.str().c_str()); |
344 |
|
tzz->DrawLatex(.85,.97,"Histograms of the sigmas"); |
345 |
|
|
346 |
|
|
347 |
|
|
444 |
for(Int_t n = 0; n<12; n++){ |
for(Int_t n = 0; n<12; n++){ |
445 |
if(ctrk.ncalev[n]==0 && ctrk.calfl[n]==0)calok = 1;//GOOD |
if(ctrk.ncalev[n]==0 && ctrk.calfl[n]==0)calok = 1;//GOOD |
446 |
|
|
447 |
|
if(ctrk.dspnum[n]==0) |
448 |
|
continue; |
449 |
|
|
450 |
nn=ctrk.dspnum[n]-1; |
nn=ctrk.dspnum[n]-1; |
451 |
/*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.* |
/*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.* |
452 |
* |
* |
519 |
/******************************************************/ |
/******************************************************/ |
520 |
/* fill histos */ |
/* fill histos */ |
521 |
for(Int_t j = 0; j < 3072; j++){ |
for(Int_t j = 0; j < 3072; j++){ |
522 |
histosig[nn]->Fill((Float_t)j,ctrk.dspsig[nn][j]); |
histosig[nn]->Fill(j,ctrk.dspsig[n][j]); |
523 |
histoped[nn]->Fill((Float_t)j,ctrk.dspped[nn][j]); |
histoped[nn]->Fill(j,ctrk.dspped[n][j]); |
524 |
if(j<1024) histosiglad[nn][0]->Fill(ctrk.dspsig[nn][j]); |
if(j<1024) histosiglad[nn][0]->Fill(ctrk.dspsig[n][j]); |
525 |
if(j>=1024 && j<2048) histosiglad[nn][1]->Fill(ctrk.dspsig[nn][j]); |
if(j>=1024 && j<2048) histosiglad[nn][1]->Fill(ctrk.dspsig[n][j]); |
526 |
if(j>=2048 && j<3072) histosiglad[nn][2]->Fill(ctrk.dspsig[nn][j]); |
if(j>=2048 && j<3072) histosiglad[nn][2]->Fill(ctrk.dspsig[n][j]); |
527 |
}; |
}; |
528 |
histoasig[nn]->Fill(1,ctrk.sig1[nn]); |
histoasig[nn]->Fill(1,ctrk.sig1[n]); |
529 |
histoasig[nn]->Fill(1025,ctrk.sig2[nn]); |
histoasig[nn]->Fill(1025,ctrk.sig2[n]); |
530 |
histoasig[nn]->Fill(2049,ctrk.sig3[nn]); |
histoasig[nn]->Fill(2049,ctrk.sig3[n]); |
531 |
histoaped[nn]->Fill(1,ctrk.ped1[nn]); |
histoaped[nn]->Fill(1,ctrk.ped1[nn]); |
532 |
histoaped[nn]->Fill(1025,ctrk.ped2[nn]); |
histoaped[nn]->Fill(1025,ctrk.ped2[n]); |
533 |
histoaped[nn]->Fill(2049,ctrk.ped3[nn]); |
histoaped[nn]->Fill(2049,ctrk.ped3[n]); |
534 |
/******************************************************/ |
/******************************************************/ |
535 |
|
|
536 |
TLine li; |
TLine li,liva1; |
537 |
li.SetLineColor(38); |
li.SetLineColor(38); |
538 |
li.SetLineStyle(3); |
li.SetLineStyle(4); |
539 |
li.SetLineWidth(2); |
li.SetLineWidth(2); |
540 |
|
liva1.SetLineColor(42); |
541 |
Float_t maxhist=0; |
liva1.SetLineStyle(3); |
542 |
|
liva1.SetLineWidth(1); |
543 |
|
|
544 |
|
Float_t maxhist=0,va1x=0,minhist=0.; |
545 |
TBox b; |
TBox b; |
546 |
/* plot PEDESTAL */ |
/* plot PEDESTAL */ |
547 |
c1->cd(); |
c1->cd(); |
560 |
if((nn+1)%2==0) histoped[nn]->GetYaxis()->SetRangeUser(700,1700); |
if((nn+1)%2==0) histoped[nn]->GetYaxis()->SetRangeUser(700,1700); |
561 |
histoaped[nn]->SetLineColor(5); |
histoaped[nn]->SetLineColor(5); |
562 |
histoaped[nn]->SetLineWidth(1); |
histoaped[nn]->SetLineWidth(1); |
563 |
if(ctrk.good0[0]==1 && ctrk.good0[1]==1){ |
// if(ctrk.good0[0]==1 && ctrk.good0[1]==1){ |
564 |
histoped[nn]->Draw("b"); |
histoped[nn]->Draw("b"); |
565 |
if(nn==1){ |
maxhist=histoped[nn]->GetMaximum(); |
566 |
maxhist=histoped[nn]->GetMaximum(); |
if(nn%2==0) minhist=2200; |
567 |
b.SetFillColor(6); |
else minhist=700; |
568 |
b.SetFillStyle(3945); |
b.SetFillColor(19); |
569 |
b.DrawBox(2944.,700.,3060.,maxhist); |
b.SetFillStyle(3954); |
570 |
|
if(nn==0){ |
571 |
b.SetFillColor(107); |
b.DrawBox(0.,minhist,2047.,maxhist); |
572 |
b.SetFillStyle(3954); |
} |
573 |
b.DrawBox(2816.,700.,2944.,maxhist); |
else if(nn==1){ |
574 |
b.DrawBox(2048.,700.,2176.,maxhist); |
b.DrawBox(128.,minhist,256.,maxhist); |
575 |
} |
b.DrawBox(384.,minhist,512.,maxhist); |
576 |
else if(nn==6){ |
b.DrawBox(896.,minhist,1024.,maxhist); |
577 |
maxhist=histoped[nn]->GetMaximum(); |
b.DrawBox(2048.,minhist,2432.,maxhist); |
578 |
b.SetFillColor(6); |
b.DrawBox(2816.,minhist,2944.,maxhist); |
579 |
b.SetFillStyle(3945); |
b.DrawBox(2944.,minhist,3070.,maxhist); |
580 |
b.DrawBox(2560.,2200.,2816.,maxhist); |
} |
581 |
b.DrawBox(1024.,2200.,1792.,maxhist); |
else if(nn==3){ |
582 |
|
b.DrawBox(0.,minhist,256.,maxhist); |
583 |
b.SetFillColor(107); |
b.DrawBox(2816.,minhist,3070.,maxhist); |
584 |
b.SetFillStyle(3954); |
} |
585 |
b.DrawBox(512.,2200.,768.,maxhist); |
else if(nn==4){ |
586 |
} |
b.DrawBox(256.,minhist,512.,maxhist); |
587 |
else if(nn==7){ |
b.DrawBox(1792.,minhist,1920.,maxhist); |
588 |
maxhist=histoped[nn]->GetMaximum(); |
b.DrawBox(2816.,minhist,3070.,maxhist); |
589 |
b.SetFillColor(107); |
} |
590 |
b.SetFillStyle(3954); |
else if(nn==5){ |
591 |
b.DrawBox(512.,700.,768.,maxhist); |
b.DrawBox(0.,minhist,256.,maxhist); |
592 |
} |
b.DrawBox(896.,minhist,1024.,maxhist); |
593 |
else if(nn==11){ |
b.DrawBox(1664.,minhist,1792.,maxhist); |
|
maxhist=histoped[nn]->GetMaximum(); |
|
|
b.SetFillColor(6); |
|
|
b.SetFillStyle(3945); |
|
|
b.DrawBox(768.,700.,1024.,maxhist); |
|
|
|
|
|
b.SetFillColor(107); |
|
|
b.SetFillStyle(3954); |
|
|
b.DrawBox(0.,700.,512.,maxhist); |
|
|
b.DrawBox(1920.,700.,2048.,maxhist); |
|
|
} |
|
594 |
} |
} |
595 |
else histoped[nn]->Draw("axis"); |
else if(nn==6){ |
596 |
|
b.DrawBox(512.,minhist,768.,maxhist); |
597 |
|
b.DrawBox(1024.,minhist,1280.,maxhist); |
598 |
|
b.DrawBox(1280.,minhist,1792.,maxhist); |
599 |
|
b.DrawBox(2560.,minhist,2816.,maxhist); |
600 |
|
} |
601 |
|
else if(nn==7){ |
602 |
|
b.DrawBox(0.,minhist,1535.,maxhist); |
603 |
|
b.DrawBox(2024.,minhist,2280.,maxhist); |
604 |
|
} |
605 |
|
else if(nn==8){ |
606 |
|
b.DrawBox(512.,minhist,768.,maxhist); |
607 |
|
} |
608 |
|
else if(nn==9){ |
609 |
|
b.DrawBox(0.,minhist,128.,maxhist); |
610 |
|
b.DrawBox(256.,minhist,384.,maxhist); |
611 |
|
b.DrawBox(512.,minhist,640.,maxhist); |
612 |
|
b.DrawBox(896.,minhist,1152.,maxhist); |
613 |
|
b.DrawBox(1280.,minhist,1535.,maxhist); |
614 |
|
b.DrawBox(1664.,minhist,1920.,maxhist); |
615 |
|
b.DrawBox(2048.,minhist,2304.,maxhist); |
616 |
|
} |
617 |
|
else if(nn==10){ |
618 |
|
b.DrawBox(0.,minhist,512.,maxhist); |
619 |
|
b.DrawBox(1024.,minhist,3070.,maxhist); |
620 |
|
} |
621 |
|
else if(nn==11){ |
622 |
|
b.DrawBox(0.,minhist,512.,maxhist); |
623 |
|
b.DrawBox(768.,minhist,1024.,maxhist); |
624 |
|
b.DrawBox(1536.,minhist,1664.,maxhist); |
625 |
|
b.DrawBox(1920.,minhist,2560.,maxhist); |
626 |
|
b.DrawBox(2816.,minhist,3070.,maxhist); |
627 |
|
} |
628 |
|
// } |
629 |
|
// else histoped[nn]->Draw("axis"); |
630 |
histoaped[nn]->Draw("same"); |
histoaped[nn]->Draw("same"); |
631 |
if((nn+1)%2==1) { |
if((nn+1)%2==1) { |
632 |
li.DrawLine(1024.5,2200,1024.5,3200); |
for(int va=1; va<24; va++){ |
633 |
li.DrawLine(2048.5,2200,2048.5,3200); |
va1x=128*va; |
634 |
|
liva1.DrawLine(va1x,2200.,va1x,3200.); |
635 |
|
} |
636 |
|
li.DrawLine(1024.5,2200.,1024.5,3200.); |
637 |
|
li.DrawLine(2048.5,2200.,2048.5,3200.); |
638 |
} |
} |
639 |
if((nn+1)%2==0) { |
if((nn+1)%2==0) { |
640 |
|
for(int va=1; va<24; va++){ |
641 |
|
va1x=128*va; |
642 |
|
liva1.DrawLine(va1x,700.,va1x,1700.); |
643 |
|
} |
644 |
li.DrawLine(1024.5,700,1024.5,1700); |
li.DrawLine(1024.5,700,1024.5,1700); |
645 |
li.DrawLine(2048.5,700,2048.5,1700); |
li.DrawLine(2048.5,700,2048.5,1700); |
646 |
} |
} |
647 |
|
c1->Update();//draw pads in canvas |
648 |
|
|
649 |
|
|
650 |
|
|
667 |
histosig[nn]->GetYaxis()->CenterTitle(); |
histosig[nn]->GetYaxis()->CenterTitle(); |
668 |
histoasig[nn]->SetLineColor(5); |
histoasig[nn]->SetLineColor(5); |
669 |
histoasig[nn]->SetLineWidth(1); |
histoasig[nn]->SetLineWidth(1); |
670 |
if(ctrk.good0[0]==1 && ctrk.good0[1]==1){ |
// if(ctrk.good0[0]==1 && ctrk.good0[1]==1){ |
671 |
histosig[nn]->Draw("b"); |
histosig[nn]->Draw("b"); |
672 |
if(nn==1){ |
maxhist=histosig[nn]->GetMaximum(); |
673 |
maxhist=histosig[nn]->GetMaximum(); |
minhist=0.; |
674 |
b.SetFillColor(6); |
b.SetFillColor(19); |
675 |
b.SetFillStyle(3945); |
b.SetFillStyle(3945); |
676 |
b.DrawBox(2944.,0.,3060.,maxhist); |
if(nn==0){ |
677 |
|
b.DrawBox(0.,minhist,2047.,maxhist); |
678 |
b.SetFillColor(107); |
} |
679 |
b.SetFillStyle(3954); |
else if(nn==1){ |
680 |
b.DrawBox(2816.,0.,2944.,maxhist); |
b.DrawBox(128.,minhist,256.,maxhist); |
681 |
b.DrawBox(2048.,0.,2176.,maxhist); |
b.DrawBox(384.,minhist,512.,maxhist); |
682 |
} |
b.DrawBox(896.,minhist,1024.,maxhist); |
683 |
else if(nn==6){ |
b.DrawBox(2048.,minhist,2432.,maxhist); |
684 |
maxhist=histosig[nn]->GetMaximum(); |
b.DrawBox(2816.,minhist,2944.,maxhist); |
685 |
b.SetFillColor(6); |
b.DrawBox(2944.,minhist,3070.,maxhist); |
|
b.SetFillStyle(3945); |
|
|
b.DrawBox(2560.,0.,2816.,maxhist); |
|
|
b.DrawBox(1024.,0.,1792.,maxhist); |
|
|
|
|
|
b.SetFillColor(107); |
|
|
b.SetFillStyle(3954); |
|
|
b.DrawBox(512.,0.,768.,maxhist); |
|
|
} |
|
|
else if(nn==7){ |
|
|
maxhist=histosig[nn]->GetMaximum(); |
|
|
b.SetFillColor(107); |
|
|
b.SetFillStyle(3954); |
|
|
b.DrawBox(512.,0.,768.,maxhist); |
|
|
} |
|
|
else if(nn==11){ |
|
|
maxhist=histosig[nn]->GetMaximum(); |
|
|
b.SetFillColor(6); |
|
|
b.SetFillStyle(3945); |
|
|
b.DrawBox(768.,0.,1024.,maxhist); |
|
|
|
|
|
b.SetFillColor(107); |
|
|
b.SetFillStyle(3954); |
|
|
b.DrawBox(0.,0.,512.,maxhist); |
|
|
b.DrawBox(1920.,0.,2048.,maxhist); |
|
|
} |
|
686 |
} |
} |
687 |
else histosig[nn]->Draw("axis"); |
else if(nn==3){ |
688 |
|
b.DrawBox(0.,minhist,256.,maxhist); |
689 |
|
b.DrawBox(2816.,minhist,3070.,maxhist); |
690 |
|
} |
691 |
|
else if(nn==4){ |
692 |
|
b.DrawBox(256.,minhist,512.,maxhist); |
693 |
|
b.DrawBox(1792.,minhist,1920.,maxhist); |
694 |
|
b.DrawBox(2816.,minhist,3070.,maxhist); |
695 |
|
} |
696 |
|
else if(nn==5){ |
697 |
|
b.DrawBox(0.,minhist,256.,maxhist); |
698 |
|
b.DrawBox(896.,minhist,1024.,maxhist); |
699 |
|
b.DrawBox(1664.,minhist,1792.,maxhist); |
700 |
|
} |
701 |
|
else if(nn==6){ |
702 |
|
b.DrawBox(512.,minhist,768.,maxhist); |
703 |
|
b.DrawBox(1024.,minhist,1280.,maxhist); |
704 |
|
b.DrawBox(1280.,minhist,1792.,maxhist); |
705 |
|
b.DrawBox(2560.,minhist,2816.,maxhist); |
706 |
|
} |
707 |
|
else if(nn==7){ |
708 |
|
b.DrawBox(0.,minhist,1535.,maxhist); |
709 |
|
b.DrawBox(2024.,minhist,2280.,maxhist); |
710 |
|
} |
711 |
|
else if(nn==8){ |
712 |
|
b.DrawBox(512.,minhist,768.,maxhist); |
713 |
|
} |
714 |
|
else if(nn==9){ |
715 |
|
b.DrawBox(0.,minhist,128.,maxhist); |
716 |
|
b.DrawBox(256.,minhist,384.,maxhist); |
717 |
|
b.DrawBox(512.,minhist,640.,maxhist); |
718 |
|
b.DrawBox(896.,minhist,1152.,maxhist); |
719 |
|
b.DrawBox(1280.,minhist,1535.,maxhist); |
720 |
|
b.DrawBox(1664.,minhist,1920.,maxhist); |
721 |
|
b.DrawBox(2048.,minhist,2304.,maxhist); |
722 |
|
} |
723 |
|
else if(nn==10){ |
724 |
|
b.DrawBox(0.,minhist,512.,maxhist); |
725 |
|
b.DrawBox(1024.,minhist,3070.,maxhist); |
726 |
|
} |
727 |
|
else if(nn==11){ |
728 |
|
b.DrawBox(0.,minhist,512.,maxhist); |
729 |
|
b.DrawBox(768.,minhist,1024.,maxhist); |
730 |
|
b.DrawBox(1536.,minhist,1664.,maxhist); |
731 |
|
b.DrawBox(1920.,minhist,2560.,maxhist); |
732 |
|
b.DrawBox(2816.,minhist,3070.,maxhist); |
733 |
|
} |
734 |
|
// } |
735 |
|
// else histosig[nn]->Draw("axis"); |
736 |
histoasig[nn]->Draw("same"); |
histoasig[nn]->Draw("same"); |
737 |
|
for(int va=1; va<24; va++){ |
738 |
|
va1x=128*va; |
739 |
|
liva1.DrawLine(va1x,0.,va1x,max); |
740 |
|
} |
741 |
li.DrawLine(1024.5,0,1024.5,max); |
li.DrawLine(1024.5,0,1024.5,max); |
742 |
li.DrawLine(2048.5,0,2048.5,max); |
li.DrawLine(2048.5,0,2048.5,max); |
743 |
|
c2->Update();//draw pads in canvas |
744 |
|
|
745 |
for(int ii=0;ii<3;ii++){ |
for(int ii=0;ii<3;ii++){ |
746 |
sig->cd(); |
sig->cd(); |
755 |
histosiglad[nn][ii]->GetXaxis()->SetTitle("SIG (ADC channels)"); |
histosiglad[nn][ii]->GetXaxis()->SetTitle("SIG (ADC channels)"); |
756 |
histosiglad[nn][ii]->GetXaxis()->CenterTitle(); |
histosiglad[nn][ii]->GetXaxis()->CenterTitle(); |
757 |
histosiglad[nn][ii]->Draw(""); |
histosiglad[nn][ii]->Draw(""); |
758 |
} |
} |
759 |
|
|
760 |
};//end loop on views |
};//end loop on views |
761 |
c1->Update();//draw pads in canvas |
c1->Update();//draw pads in canvas |
762 |
c2->Update();//draw pads in canvas |
c2->Update();//draw pads in canvas |