2 |
* FTrkCalibQLookExpert.cxx |
* FTrkCalibQLookExpert.cxx |
3 |
* |
* |
4 |
* autor: D.Fedele |
* autor: D.Fedele |
5 |
* version v1r10 |
* version v1r18 |
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 |
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 |
|
|
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; |
279 |
|
|
419 |
for(Int_t n = 0; n<12; n++){ |
for(Int_t n = 0; n<12; n++){ |
420 |
if(ctrk.ncalev[n]==0 && ctrk.calfl[n]==0)calok = 1;//GOOD |
if(ctrk.ncalev[n]==0 && ctrk.calfl[n]==0)calok = 1;//GOOD |
421 |
|
|
422 |
|
if(ctrk.dspnum[n]==0) |
423 |
|
continue; |
424 |
|
|
425 |
nn=ctrk.dspnum[n]-1; |
nn=ctrk.dspnum[n]-1; |
426 |
/*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.* |
/*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.* |
427 |
* |
* |
494 |
/******************************************************/ |
/******************************************************/ |
495 |
/* fill histos */ |
/* fill histos */ |
496 |
for(Int_t j = 0; j < 3072; j++){ |
for(Int_t j = 0; j < 3072; j++){ |
497 |
histosig[nn]->Fill((Float_t)j,ctrk.dspsig[nn][j]); |
histosig[nn]->Fill(j,ctrk.dspsig[n][j]); |
498 |
histoped[nn]->Fill((Float_t)j,ctrk.dspped[nn][j]); |
histoped[nn]->Fill(j,ctrk.dspped[n][j]); |
499 |
if(j<1024) histosiglad[nn][0]->Fill(ctrk.dspsig[nn][j]); |
if(j<1024) histosiglad[nn][0]->Fill(ctrk.dspsig[n][j]); |
500 |
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]); |
501 |
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]); |
502 |
}; |
}; |
503 |
histoasig[nn]->Fill(1,ctrk.sig1[nn]); |
histoasig[nn]->Fill(1,ctrk.sig1[n]); |
504 |
histoasig[nn]->Fill(1025,ctrk.sig2[nn]); |
histoasig[nn]->Fill(1025,ctrk.sig2[n]); |
505 |
histoasig[nn]->Fill(2049,ctrk.sig3[nn]); |
histoasig[nn]->Fill(2049,ctrk.sig3[n]); |
506 |
histoaped[nn]->Fill(1,ctrk.ped1[nn]); |
histoaped[nn]->Fill(1,ctrk.ped1[nn]); |
507 |
histoaped[nn]->Fill(1025,ctrk.ped2[nn]); |
histoaped[nn]->Fill(1025,ctrk.ped2[n]); |
508 |
histoaped[nn]->Fill(2049,ctrk.ped3[nn]); |
histoaped[nn]->Fill(2049,ctrk.ped3[n]); |
509 |
/******************************************************/ |
/******************************************************/ |
510 |
|
|
511 |
TLine li; |
TLine li,liva1; |
512 |
li.SetLineColor(38); |
li.SetLineColor(38); |
513 |
li.SetLineStyle(3); |
li.SetLineStyle(4); |
514 |
li.SetLineWidth(2); |
li.SetLineWidth(2); |
515 |
|
liva1.SetLineColor(42); |
516 |
Float_t maxhist=0; |
liva1.SetLineStyle(3); |
517 |
|
liva1.SetLineWidth(1); |
518 |
|
|
519 |
|
Float_t maxhist=0,va1x=0; |
520 |
TBox b; |
TBox b; |
521 |
/* plot PEDESTAL */ |
/* plot PEDESTAL */ |
522 |
c1->cd(); |
c1->cd(); |
535 |
if((nn+1)%2==0) histoped[nn]->GetYaxis()->SetRangeUser(700,1700); |
if((nn+1)%2==0) histoped[nn]->GetYaxis()->SetRangeUser(700,1700); |
536 |
histoaped[nn]->SetLineColor(5); |
histoaped[nn]->SetLineColor(5); |
537 |
histoaped[nn]->SetLineWidth(1); |
histoaped[nn]->SetLineWidth(1); |
538 |
if(ctrk.good0[0]==1 && ctrk.good0[1]==1){ |
// if(ctrk.good0[0]==1 && ctrk.good0[1]==1){ |
539 |
histoped[nn]->Draw("b"); |
histoped[nn]->Draw("b"); |
540 |
if(nn==1){ |
maxhist=histoped[nn]->GetMaximum(); |
541 |
maxhist=histoped[nn]->GetMaximum(); |
if(nn==0){ |
542 |
|
b.SetFillColor(107); |
543 |
|
b.SetFillStyle(3945); |
544 |
|
b.DrawBox(768.,2200.,2047.,maxhist); |
545 |
|
} |
546 |
|
else if(nn==1){ |
547 |
b.SetFillColor(6); |
b.SetFillColor(6); |
548 |
b.SetFillStyle(3945); |
b.SetFillStyle(3945); |
549 |
b.DrawBox(2944.,700.,3060.,maxhist); |
b.DrawBox(2944.,700.,3060.,maxhist); |
554 |
b.DrawBox(2048.,700.,2176.,maxhist); |
b.DrawBox(2048.,700.,2176.,maxhist); |
555 |
} |
} |
556 |
else if(nn==4){ |
else if(nn==4){ |
|
maxhist=histoped[nn]->GetMaximum(); |
|
557 |
b.SetFillColor(107); |
b.SetFillColor(107); |
558 |
b.SetFillStyle(3954); |
b.SetFillStyle(3954); |
559 |
b.DrawBox(384.,2200.,512.,maxhist); |
b.DrawBox(384.,2200.,512.,maxhist); |
560 |
} |
} |
561 |
else if(nn==6){ |
else if(nn==6){ |
|
maxhist=histoped[nn]->GetMaximum(); |
|
562 |
b.SetFillColor(6); |
b.SetFillColor(6); |
563 |
b.SetFillStyle(3945); |
b.SetFillStyle(3945); |
564 |
b.DrawBox(2560.,2200.,2816.,maxhist); |
b.DrawBox(2560.,2200.,2816.,maxhist); |
565 |
b.DrawBox(1024.,2200.,1792.,maxhist); |
b.DrawBox(1024.,2200.,1535.,maxhist); |
566 |
|
|
567 |
b.SetFillColor(107); |
b.SetFillColor(107); |
568 |
b.SetFillStyle(3954); |
b.SetFillStyle(3954); |
569 |
b.DrawBox(512.,2200.,768.,maxhist); |
b.DrawBox(512.,2200.,768.,maxhist); |
570 |
|
b.DrawBox(1536.,2200.,1792.,maxhist); |
571 |
} |
} |
572 |
else if(nn==7){ |
else if(nn==7){ |
|
maxhist=histoped[nn]->GetMaximum(); |
|
573 |
b.SetFillColor(107); |
b.SetFillColor(107); |
574 |
b.SetFillStyle(3954); |
b.SetFillStyle(3954); |
575 |
b.DrawBox(512.,700.,768.,maxhist); |
b.DrawBox(512.,700.,768.,maxhist); |
576 |
} |
} |
577 |
|
else if(nn==8){ |
578 |
|
b.SetFillColor(107); |
579 |
|
b.SetFillStyle(3954); |
580 |
|
b.DrawBox(512.,2200.,768.,maxhist); |
581 |
|
} |
582 |
|
else if(nn==9){ |
583 |
|
b.SetFillColor(107); |
584 |
|
b.SetFillStyle(3954); |
585 |
|
b.DrawBox(256.,700.,384.,maxhist); |
586 |
|
} |
587 |
|
else if(nn==10){ |
588 |
|
b.SetFillColor(107); |
589 |
|
b.SetFillStyle(3954); |
590 |
|
b.DrawBox(2560.,2200.,3070.,maxhist); |
591 |
|
} |
592 |
else if(nn==11){ |
else if(nn==11){ |
|
maxhist=histoped[nn]->GetMaximum(); |
|
593 |
b.SetFillColor(6); |
b.SetFillColor(6); |
594 |
b.SetFillStyle(3945); |
b.SetFillStyle(3945); |
595 |
b.DrawBox(768.,700.,1024.,maxhist); |
b.DrawBox(768.,700.,1024.,maxhist); |
597 |
b.SetFillColor(107); |
b.SetFillColor(107); |
598 |
b.SetFillStyle(3954); |
b.SetFillStyle(3954); |
599 |
b.DrawBox(0.,700.,512.,maxhist); |
b.DrawBox(0.,700.,512.,maxhist); |
600 |
b.DrawBox(1920.,700.,2048.,maxhist); |
b.DrawBox(1920.,700.,2560.,maxhist); |
601 |
} |
} |
602 |
} |
// } |
603 |
else histoped[nn]->Draw("axis"); |
// else histoped[nn]->Draw("axis"); |
604 |
histoaped[nn]->Draw("same"); |
histoaped[nn]->Draw("same"); |
605 |
if((nn+1)%2==1) { |
if((nn+1)%2==1) { |
606 |
li.DrawLine(1024.5,2200,1024.5,3200); |
for(int va=1; va<24; va++){ |
607 |
li.DrawLine(2048.5,2200,2048.5,3200); |
va1x=128*va; |
608 |
|
liva1.DrawLine(va1x,2200.,va1x,3200.); |
609 |
|
} |
610 |
|
li.DrawLine(1024.5,2200.,1024.5,3200.); |
611 |
|
li.DrawLine(2048.5,2200.,2048.5,3200.); |
612 |
} |
} |
613 |
if((nn+1)%2==0) { |
if((nn+1)%2==0) { |
614 |
|
for(int va=1; va<24; va++){ |
615 |
|
va1x=128*va; |
616 |
|
liva1.DrawLine(va1x,700.,va1x,1700.); |
617 |
|
} |
618 |
li.DrawLine(1024.5,700,1024.5,1700); |
li.DrawLine(1024.5,700,1024.5,1700); |
619 |
li.DrawLine(2048.5,700,2048.5,1700); |
li.DrawLine(2048.5,700,2048.5,1700); |
620 |
} |
} |
621 |
|
c1->Update();//draw pads in canvas |
622 |
|
|
623 |
|
|
624 |
|
|
641 |
histosig[nn]->GetYaxis()->CenterTitle(); |
histosig[nn]->GetYaxis()->CenterTitle(); |
642 |
histoasig[nn]->SetLineColor(5); |
histoasig[nn]->SetLineColor(5); |
643 |
histoasig[nn]->SetLineWidth(1); |
histoasig[nn]->SetLineWidth(1); |
644 |
if(ctrk.good0[0]==1 && ctrk.good0[1]==1){ |
// if(ctrk.good0[0]==1 && ctrk.good0[1]==1){ |
645 |
histosig[nn]->Draw("b"); |
histosig[nn]->Draw("b"); |
646 |
if(nn==1){ |
maxhist=histosig[nn]->GetMaximum(); |
647 |
maxhist=histosig[nn]->GetMaximum(); |
if(nn==0){ |
648 |
|
b.SetFillColor(107); |
649 |
|
b.SetFillStyle(3945); |
650 |
|
b.DrawBox(768.,0.,2047.,maxhist); |
651 |
|
} |
652 |
|
else if(nn==1){ |
653 |
b.SetFillColor(6); |
b.SetFillColor(6); |
654 |
b.SetFillStyle(3945); |
b.SetFillStyle(3945); |
655 |
b.DrawBox(2944.,0.,3060.,maxhist); |
b.DrawBox(2944.,0.,3070.,maxhist); |
656 |
|
|
657 |
b.SetFillColor(107); |
b.SetFillColor(107); |
658 |
b.SetFillStyle(3954); |
b.SetFillStyle(3954); |
660 |
b.DrawBox(2048.,0.,2176.,maxhist); |
b.DrawBox(2048.,0.,2176.,maxhist); |
661 |
} |
} |
662 |
else if(nn==4){ |
else if(nn==4){ |
|
maxhist=histosig[nn]->GetMaximum(); |
|
663 |
b.SetFillColor(107); |
b.SetFillColor(107); |
664 |
b.SetFillStyle(3954); |
b.SetFillStyle(3954); |
665 |
b.DrawBox(384.,0.,512.,maxhist); |
b.DrawBox(384.,0.,512.,maxhist); |
666 |
} |
} |
667 |
else if(nn==6){ |
else if(nn==6){ |
|
maxhist=histosig[nn]->GetMaximum(); |
|
668 |
b.SetFillColor(6); |
b.SetFillColor(6); |
669 |
b.SetFillStyle(3945); |
b.SetFillStyle(3945); |
670 |
b.DrawBox(2560.,0.,2816.,maxhist); |
b.DrawBox(2560.,0.,2816.,maxhist); |
671 |
b.DrawBox(1024.,0.,1792.,maxhist); |
b.DrawBox(1024.,0.,1535.,maxhist); |
672 |
|
|
673 |
b.SetFillColor(107); |
b.SetFillColor(107); |
674 |
b.SetFillStyle(3954); |
b.SetFillStyle(3954); |
675 |
b.DrawBox(512.,0.,768.,maxhist); |
b.DrawBox(512.,0.,768.,maxhist); |
676 |
|
b.DrawBox(1536.,0.,1792.,maxhist); |
677 |
} |
} |
678 |
else if(nn==7){ |
else if(nn==7){ |
|
maxhist=histosig[nn]->GetMaximum(); |
|
679 |
b.SetFillColor(107); |
b.SetFillColor(107); |
680 |
b.SetFillStyle(3954); |
b.SetFillStyle(3954); |
681 |
b.DrawBox(512.,0.,768.,maxhist); |
b.DrawBox(512.,0.,768.,maxhist); |
682 |
} |
} |
683 |
|
else if(nn==8){ |
684 |
|
b.SetFillColor(107); |
685 |
|
b.SetFillStyle(3954); |
686 |
|
b.DrawBox(512.,0.,768.,maxhist); |
687 |
|
} |
688 |
|
else if(nn==9){ |
689 |
|
b.SetFillColor(107); |
690 |
|
b.SetFillStyle(3954); |
691 |
|
b.DrawBox(256.,0.,384.,maxhist); |
692 |
|
} |
693 |
|
else if(nn==10){ |
694 |
|
b.SetFillColor(107); |
695 |
|
b.SetFillStyle(3954); |
696 |
|
b.DrawBox(2560.,0.,3070.,maxhist); |
697 |
|
} |
698 |
else if(nn==11){ |
else if(nn==11){ |
|
maxhist=histosig[nn]->GetMaximum(); |
|
699 |
b.SetFillColor(6); |
b.SetFillColor(6); |
700 |
b.SetFillStyle(3945); |
b.SetFillStyle(3945); |
701 |
b.DrawBox(768.,0.,1024.,maxhist); |
b.DrawBox(768.,0.,1024.,maxhist); |
703 |
b.SetFillColor(107); |
b.SetFillColor(107); |
704 |
b.SetFillStyle(3954); |
b.SetFillStyle(3954); |
705 |
b.DrawBox(0.,0.,512.,maxhist); |
b.DrawBox(0.,0.,512.,maxhist); |
706 |
b.DrawBox(1920.,0.,2048.,maxhist); |
b.DrawBox(1920.,0.,2560.,maxhist); |
707 |
} |
} |
708 |
} |
// } |
709 |
else histosig[nn]->Draw("axis"); |
// else histosig[nn]->Draw("axis"); |
710 |
histoasig[nn]->Draw("same"); |
histoasig[nn]->Draw("same"); |
711 |
|
for(int va=1; va<24; va++){ |
712 |
|
va1x=128*va; |
713 |
|
liva1.DrawLine(va1x,0.,va1x,max); |
714 |
|
} |
715 |
li.DrawLine(1024.5,0,1024.5,max); |
li.DrawLine(1024.5,0,1024.5,max); |
716 |
li.DrawLine(2048.5,0,2048.5,max); |
li.DrawLine(2048.5,0,2048.5,max); |
717 |
|
c2->Update();//draw pads in canvas |
718 |
|
|
719 |
for(int ii=0;ii<3;ii++){ |
for(int ii=0;ii<3;ii++){ |
720 |
sig->cd(); |
sig->cd(); |
729 |
histosiglad[nn][ii]->GetXaxis()->SetTitle("SIG (ADC channels)"); |
histosiglad[nn][ii]->GetXaxis()->SetTitle("SIG (ADC channels)"); |
730 |
histosiglad[nn][ii]->GetXaxis()->CenterTitle(); |
histosiglad[nn][ii]->GetXaxis()->CenterTitle(); |
731 |
histosiglad[nn][ii]->Draw(""); |
histosiglad[nn][ii]->Draw(""); |
732 |
} |
} |
733 |
|
|
734 |
};//end loop on views |
};//end loop on views |
735 |
c1->Update();//draw pads in canvas |
c1->Update();//draw pads in canvas |
736 |
c2->Update();//draw pads in canvas |
c2->Update();//draw pads in canvas |