2 |
* FTrkCalibQLookExpert.cxx |
* FTrkCalibQLookExpert.cxx |
3 |
* |
* |
4 |
* autor: D.Fedele |
* autor: D.Fedele |
5 |
* version v1r08 |
* version v1r20 |
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 |
|
|
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 |
|
ULong64_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; |
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"); |
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); |
550 |
|
|
551 |
b.SetFillColor(107); |
b.SetFillColor(107); |
552 |
b.SetFillStyle(3954); |
b.SetFillStyle(3954); |
553 |
|
//b.DrawBox(384.,700.,512.,maxhist); |
554 |
b.DrawBox(2816.,700.,2944.,maxhist); |
b.DrawBox(2816.,700.,2944.,maxhist); |
555 |
b.DrawBox(2048.,700.,2176.,maxhist); |
b.DrawBox(2048.,700.,2176.,maxhist); |
556 |
} |
} |
557 |
|
else if(nn==4){ |
558 |
|
b.SetFillColor(107); |
559 |
|
b.SetFillStyle(3954); |
560 |
|
b.DrawBox(384.,2200.,512.,maxhist); |
561 |
|
} |
562 |
else if(nn==6){ |
else if(nn==6){ |
|
maxhist=histoped[nn]->GetMaximum(); |
|
563 |
b.SetFillColor(6); |
b.SetFillColor(6); |
564 |
b.SetFillStyle(3945); |
b.SetFillStyle(3945); |
565 |
b.DrawBox(2560.,2200.,2816.,maxhist); |
b.DrawBox(2560.,2200.,2816.,maxhist); |
566 |
b.DrawBox(1024.,2200.,1792.,maxhist); |
b.DrawBox(1024.,2200.,1535.,maxhist); |
567 |
|
|
568 |
b.SetFillColor(107); |
b.SetFillColor(107); |
569 |
b.SetFillStyle(3954); |
b.SetFillStyle(3954); |
570 |
b.DrawBox(512.,2200.,768.,maxhist); |
b.DrawBox(512.,2200.,768.,maxhist); |
571 |
|
b.DrawBox(1536.,2200.,1792.,maxhist); |
572 |
} |
} |
573 |
else if(nn==7){ |
else if(nn==7){ |
|
maxhist=histoped[nn]->GetMaximum(); |
|
574 |
b.SetFillColor(107); |
b.SetFillColor(107); |
575 |
b.SetFillStyle(3954); |
b.SetFillStyle(3954); |
576 |
b.DrawBox(512.,700.,768.,maxhist); |
b.DrawBox(512.,700.,768.,maxhist); |
577 |
} |
} |
578 |
|
else if(nn==8){ |
579 |
|
b.SetFillColor(107); |
580 |
|
b.SetFillStyle(3954); |
581 |
|
b.DrawBox(512.,2200.,768.,maxhist); |
582 |
|
} |
583 |
|
else if(nn==9){ |
584 |
|
b.SetFillColor(107); |
585 |
|
b.SetFillStyle(3954); |
586 |
|
b.DrawBox(256.,700.,384.,maxhist); |
587 |
|
//b.DrawBox(1280.,700.,1408.,maxhist); |
588 |
|
//b.DrawBox(1792.,700.,1920.,maxhist); |
589 |
|
} |
590 |
|
else if(nn==10){ |
591 |
|
b.SetFillColor(107); |
592 |
|
b.SetFillStyle(3954); |
593 |
|
b.DrawBox(2048.,2200.,3070.,maxhist); |
594 |
|
} |
595 |
else if(nn==11){ |
else if(nn==11){ |
|
maxhist=histoped[nn]->GetMaximum(); |
|
596 |
b.SetFillColor(6); |
b.SetFillColor(6); |
597 |
b.SetFillStyle(3945); |
b.SetFillStyle(3945); |
598 |
b.DrawBox(768.,700.,1024.,maxhist); |
b.DrawBox(768.,700.,1024.,maxhist); |
600 |
b.SetFillColor(107); |
b.SetFillColor(107); |
601 |
b.SetFillStyle(3954); |
b.SetFillStyle(3954); |
602 |
b.DrawBox(0.,700.,512.,maxhist); |
b.DrawBox(0.,700.,512.,maxhist); |
603 |
b.DrawBox(1920.,700.,2048.,maxhist); |
b.DrawBox(1920.,700.,2560.,maxhist); |
604 |
} |
} |
605 |
} |
// } |
606 |
else histoped[nn]->Draw("axis"); |
// else histoped[nn]->Draw("axis"); |
607 |
histoaped[nn]->Draw("same"); |
histoaped[nn]->Draw("same"); |
608 |
if((nn+1)%2==1) { |
if((nn+1)%2==1) { |
609 |
li.DrawLine(1024.5,2200,1024.5,3200); |
for(int va=1; va<24; va++){ |
610 |
li.DrawLine(2048.5,2200,2048.5,3200); |
va1x=128*va; |
611 |
|
liva1.DrawLine(va1x,2200.,va1x,3200.); |
612 |
|
} |
613 |
|
li.DrawLine(1024.5,2200.,1024.5,3200.); |
614 |
|
li.DrawLine(2048.5,2200.,2048.5,3200.); |
615 |
} |
} |
616 |
if((nn+1)%2==0) { |
if((nn+1)%2==0) { |
617 |
|
for(int va=1; va<24; va++){ |
618 |
|
va1x=128*va; |
619 |
|
liva1.DrawLine(va1x,700.,va1x,1700.); |
620 |
|
} |
621 |
li.DrawLine(1024.5,700,1024.5,1700); |
li.DrawLine(1024.5,700,1024.5,1700); |
622 |
li.DrawLine(2048.5,700,2048.5,1700); |
li.DrawLine(2048.5,700,2048.5,1700); |
623 |
} |
} |
624 |
|
c1->Update();//draw pads in canvas |
625 |
|
|
626 |
|
|
627 |
|
|
644 |
histosig[nn]->GetYaxis()->CenterTitle(); |
histosig[nn]->GetYaxis()->CenterTitle(); |
645 |
histoasig[nn]->SetLineColor(5); |
histoasig[nn]->SetLineColor(5); |
646 |
histoasig[nn]->SetLineWidth(1); |
histoasig[nn]->SetLineWidth(1); |
647 |
if(ctrk.good0[0]==1 && ctrk.good0[1]==1){ |
// if(ctrk.good0[0]==1 && ctrk.good0[1]==1){ |
648 |
histosig[nn]->Draw("b"); |
histosig[nn]->Draw("b"); |
649 |
if(nn==1){ |
maxhist=histosig[nn]->GetMaximum(); |
650 |
maxhist=histosig[nn]->GetMaximum(); |
if(nn==0){ |
651 |
|
b.SetFillColor(107); |
652 |
|
b.SetFillStyle(3945); |
653 |
|
b.DrawBox(768.,0.,2047.,maxhist); |
654 |
|
} |
655 |
|
else if(nn==1){ |
656 |
b.SetFillColor(6); |
b.SetFillColor(6); |
657 |
b.SetFillStyle(3945); |
b.SetFillStyle(3945); |
658 |
b.DrawBox(2944.,0.,3060.,maxhist); |
b.DrawBox(2944.,0.,3070.,maxhist); |
659 |
|
|
660 |
b.SetFillColor(107); |
b.SetFillColor(107); |
661 |
b.SetFillStyle(3954); |
b.SetFillStyle(3954); |
662 |
|
//b.DrawBox(384.,0.,512.,maxhist); |
663 |
b.DrawBox(2816.,0.,2944.,maxhist); |
b.DrawBox(2816.,0.,2944.,maxhist); |
664 |
b.DrawBox(2048.,0.,2176.,maxhist); |
b.DrawBox(2048.,0.,2176.,maxhist); |
665 |
} |
} |
666 |
|
else if(nn==4){ |
667 |
|
b.SetFillColor(107); |
668 |
|
b.SetFillStyle(3954); |
669 |
|
b.DrawBox(384.,0.,512.,maxhist); |
670 |
|
} |
671 |
else if(nn==6){ |
else if(nn==6){ |
|
maxhist=histosig[nn]->GetMaximum(); |
|
672 |
b.SetFillColor(6); |
b.SetFillColor(6); |
673 |
b.SetFillStyle(3945); |
b.SetFillStyle(3945); |
674 |
b.DrawBox(2560.,0.,2816.,maxhist); |
b.DrawBox(2560.,0.,2816.,maxhist); |
675 |
b.DrawBox(1024.,0.,1792.,maxhist); |
b.DrawBox(1024.,0.,1535.,maxhist); |
676 |
|
|
677 |
b.SetFillColor(107); |
b.SetFillColor(107); |
678 |
b.SetFillStyle(3954); |
b.SetFillStyle(3954); |
679 |
b.DrawBox(512.,0.,768.,maxhist); |
b.DrawBox(512.,0.,768.,maxhist); |
680 |
|
b.DrawBox(1536.,0.,1792.,maxhist); |
681 |
} |
} |
682 |
else if(nn==7){ |
else if(nn==7){ |
|
maxhist=histosig[nn]->GetMaximum(); |
|
683 |
b.SetFillColor(107); |
b.SetFillColor(107); |
684 |
b.SetFillStyle(3954); |
b.SetFillStyle(3954); |
685 |
b.DrawBox(512.,0.,768.,maxhist); |
b.DrawBox(512.,0.,768.,maxhist); |
686 |
} |
} |
687 |
|
else if(nn==8){ |
688 |
|
b.SetFillColor(107); |
689 |
|
b.SetFillStyle(3954); |
690 |
|
b.DrawBox(512.,0.,768.,maxhist); |
691 |
|
} |
692 |
|
else if(nn==9){ |
693 |
|
b.SetFillColor(107); |
694 |
|
b.SetFillStyle(3954); |
695 |
|
b.DrawBox(256.,0.,384.,maxhist); |
696 |
|
//b.DrawBox(1280.,0.,1408.,maxhist); |
697 |
|
//b.DrawBox(1792.,0.,1920.,maxhist); |
698 |
|
} |
699 |
|
else if(nn==10){ |
700 |
|
b.SetFillColor(107); |
701 |
|
b.SetFillStyle(3954); |
702 |
|
b.DrawBox(2048.,0.,3070.,maxhist); |
703 |
|
} |
704 |
else if(nn==11){ |
else if(nn==11){ |
|
maxhist=histosig[nn]->GetMaximum(); |
|
705 |
b.SetFillColor(6); |
b.SetFillColor(6); |
706 |
b.SetFillStyle(3945); |
b.SetFillStyle(3945); |
707 |
b.DrawBox(768.,0.,1024.,maxhist); |
b.DrawBox(768.,0.,1024.,maxhist); |
709 |
b.SetFillColor(107); |
b.SetFillColor(107); |
710 |
b.SetFillStyle(3954); |
b.SetFillStyle(3954); |
711 |
b.DrawBox(0.,0.,512.,maxhist); |
b.DrawBox(0.,0.,512.,maxhist); |
712 |
b.DrawBox(1920.,0.,2048.,maxhist); |
b.DrawBox(1920.,0.,2560.,maxhist); |
713 |
} |
} |
714 |
} |
// } |
715 |
else histosig[nn]->Draw("axis"); |
// else histosig[nn]->Draw("axis"); |
716 |
histoasig[nn]->Draw("same"); |
histoasig[nn]->Draw("same"); |
717 |
|
for(int va=1; va<24; va++){ |
718 |
|
va1x=128*va; |
719 |
|
liva1.DrawLine(va1x,0.,va1x,max); |
720 |
|
} |
721 |
li.DrawLine(1024.5,0,1024.5,max); |
li.DrawLine(1024.5,0,1024.5,max); |
722 |
li.DrawLine(2048.5,0,2048.5,max); |
li.DrawLine(2048.5,0,2048.5,max); |
723 |
|
c2->Update();//draw pads in canvas |
724 |
|
|
725 |
for(int ii=0;ii<3;ii++){ |
for(int ii=0;ii<3;ii++){ |
726 |
sig->cd(); |
sig->cd(); |
735 |
histosiglad[nn][ii]->GetXaxis()->SetTitle("SIG (ADC channels)"); |
histosiglad[nn][ii]->GetXaxis()->SetTitle("SIG (ADC channels)"); |
736 |
histosiglad[nn][ii]->GetXaxis()->CenterTitle(); |
histosiglad[nn][ii]->GetXaxis()->CenterTitle(); |
737 |
histosiglad[nn][ii]->Draw(""); |
histosiglad[nn][ii]->Draw(""); |
738 |
} |
} |
739 |
|
|
740 |
};//end loop on views |
};//end loop on views |
741 |
c1->Update();//draw pads in canvas |
c1->Update();//draw pads in canvas |
742 |
c2->Update();//draw pads in canvas |
c2->Update();//draw pads in canvas |