| 2 |
* FTrkCalibQLookExpert.cxx |
* FTrkCalibQLookExpert.cxx |
| 3 |
* |
* |
| 4 |
* autor: D.Fedele |
* autor: D.Fedele |
| 5 |
* version v1r08 |
* version v1r16 |
| 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 |
|
|
| 276 |
Int_t risposta=0; |
Int_t risposta=0; |
| 277 |
stringstream fromfile; |
stringstream fromfile; |
| 278 |
|
|
| 279 |
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))]<<" --"; |
| 280 |
|
|
| 281 |
gStyle->SetLabelSize(0.07,"x"); |
gStyle->SetLabelSize(0.07,"x"); |
| 282 |
gStyle->SetLabelSize(0.07,"y"); |
gStyle->SetLabelSize(0.07,"y"); |
| 418 |
for(Int_t n = 0; n<12; n++){ |
for(Int_t n = 0; n<12; n++){ |
| 419 |
if(ctrk.ncalev[n]==0 && ctrk.calfl[n]==0)calok = 1;//GOOD |
if(ctrk.ncalev[n]==0 && ctrk.calfl[n]==0)calok = 1;//GOOD |
| 420 |
|
|
| 421 |
|
if(ctrk.dspnum[n]==0) |
| 422 |
|
continue; |
| 423 |
|
|
| 424 |
nn=ctrk.dspnum[n]-1; |
nn=ctrk.dspnum[n]-1; |
| 425 |
/*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.* |
/*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.* |
| 426 |
* |
* |
| 493 |
/******************************************************/ |
/******************************************************/ |
| 494 |
/* fill histos */ |
/* fill histos */ |
| 495 |
for(Int_t j = 0; j < 3072; j++){ |
for(Int_t j = 0; j < 3072; j++){ |
| 496 |
histosig[nn]->Fill((Float_t)j,ctrk.dspsig[nn][j]); |
histosig[nn]->Fill(j,ctrk.dspsig[n][j]); |
| 497 |
histoped[nn]->Fill((Float_t)j,ctrk.dspped[nn][j]); |
histoped[nn]->Fill(j,ctrk.dspped[n][j]); |
| 498 |
if(j<1024) histosiglad[nn][0]->Fill(ctrk.dspsig[nn][j]); |
if(j<1024) histosiglad[nn][0]->Fill(ctrk.dspsig[n][j]); |
| 499 |
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]); |
| 500 |
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]); |
| 501 |
}; |
}; |
| 502 |
histoasig[nn]->Fill(1,ctrk.sig1[nn]); |
histoasig[nn]->Fill(1,ctrk.sig1[n]); |
| 503 |
histoasig[nn]->Fill(1025,ctrk.sig2[nn]); |
histoasig[nn]->Fill(1025,ctrk.sig2[n]); |
| 504 |
histoasig[nn]->Fill(2049,ctrk.sig3[nn]); |
histoasig[nn]->Fill(2049,ctrk.sig3[n]); |
| 505 |
histoaped[nn]->Fill(1,ctrk.ped1[nn]); |
histoaped[nn]->Fill(1,ctrk.ped1[nn]); |
| 506 |
histoaped[nn]->Fill(1025,ctrk.ped2[nn]); |
histoaped[nn]->Fill(1025,ctrk.ped2[n]); |
| 507 |
histoaped[nn]->Fill(2049,ctrk.ped3[nn]); |
histoaped[nn]->Fill(2049,ctrk.ped3[n]); |
| 508 |
/******************************************************/ |
/******************************************************/ |
| 509 |
|
|
| 510 |
TLine li; |
TLine li; |
| 531 |
if((nn+1)%2==0) histoped[nn]->GetYaxis()->SetRangeUser(700,1700); |
if((nn+1)%2==0) histoped[nn]->GetYaxis()->SetRangeUser(700,1700); |
| 532 |
histoaped[nn]->SetLineColor(5); |
histoaped[nn]->SetLineColor(5); |
| 533 |
histoaped[nn]->SetLineWidth(1); |
histoaped[nn]->SetLineWidth(1); |
| 534 |
if(ctrk.good0[0]==1 && ctrk.good0[1]==1){ |
// if(ctrk.good0[0]==1 && ctrk.good0[1]==1){ |
| 535 |
histoped[nn]->Draw("b"); |
histoped[nn]->Draw("b"); |
| 536 |
if(nn==1){ |
if(nn==0){ |
| 537 |
|
maxhist=histoped[nn]->GetMaximum(); |
| 538 |
|
b.SetFillColor(107); |
| 539 |
|
b.SetFillStyle(3945); |
| 540 |
|
b.DrawBox(768.,700.,2047.,maxhist); |
| 541 |
|
} |
| 542 |
|
else if(nn==1){ |
| 543 |
maxhist=histoped[nn]->GetMaximum(); |
maxhist=histoped[nn]->GetMaximum(); |
| 544 |
b.SetFillColor(6); |
b.SetFillColor(6); |
| 545 |
b.SetFillStyle(3945); |
b.SetFillStyle(3945); |
| 550 |
b.DrawBox(2816.,700.,2944.,maxhist); |
b.DrawBox(2816.,700.,2944.,maxhist); |
| 551 |
b.DrawBox(2048.,700.,2176.,maxhist); |
b.DrawBox(2048.,700.,2176.,maxhist); |
| 552 |
} |
} |
| 553 |
|
else if(nn==4){ |
| 554 |
|
maxhist=histoped[nn]->GetMaximum(); |
| 555 |
|
b.SetFillColor(107); |
| 556 |
|
b.SetFillStyle(3954); |
| 557 |
|
b.DrawBox(384.,2200.,512.,maxhist); |
| 558 |
|
} |
| 559 |
else if(nn==6){ |
else if(nn==6){ |
| 560 |
maxhist=histoped[nn]->GetMaximum(); |
maxhist=histoped[nn]->GetMaximum(); |
| 561 |
b.SetFillColor(6); |
b.SetFillColor(6); |
| 562 |
b.SetFillStyle(3945); |
b.SetFillStyle(3945); |
| 563 |
b.DrawBox(2560.,2200.,2816.,maxhist); |
b.DrawBox(2560.,2200.,2816.,maxhist); |
| 564 |
b.DrawBox(1024.,2200.,1792.,maxhist); |
b.DrawBox(1024.,2200.,1535.,maxhist); |
| 565 |
|
|
| 566 |
b.SetFillColor(107); |
b.SetFillColor(107); |
| 567 |
b.SetFillStyle(3954); |
b.SetFillStyle(3954); |
| 568 |
b.DrawBox(512.,2200.,768.,maxhist); |
b.DrawBox(512.,2200.,768.,maxhist); |
| 569 |
|
b.DrawBox(1536.,2200.,1792.,maxhist); |
| 570 |
} |
} |
| 571 |
else if(nn==7){ |
else if(nn==7){ |
| 572 |
maxhist=histoped[nn]->GetMaximum(); |
maxhist=histoped[nn]->GetMaximum(); |
| 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==9){ |
| 578 |
|
maxhist=histoped[nn]->GetMaximum(); |
| 579 |
|
b.SetFillColor(107); |
| 580 |
|
b.SetFillStyle(3954); |
| 581 |
|
b.DrawBox(256.,700.,384.,maxhist); |
| 582 |
|
} |
| 583 |
else if(nn==11){ |
else if(nn==11){ |
| 584 |
maxhist=histoped[nn]->GetMaximum(); |
maxhist=histoped[nn]->GetMaximum(); |
| 585 |
b.SetFillColor(6); |
b.SetFillColor(6); |
| 589 |
b.SetFillColor(107); |
b.SetFillColor(107); |
| 590 |
b.SetFillStyle(3954); |
b.SetFillStyle(3954); |
| 591 |
b.DrawBox(0.,700.,512.,maxhist); |
b.DrawBox(0.,700.,512.,maxhist); |
| 592 |
b.DrawBox(1920.,700.,2048.,maxhist); |
b.DrawBox(1920.,700.,2560.,maxhist); |
| 593 |
} |
} |
| 594 |
} |
// } |
| 595 |
else histoped[nn]->Draw("axis"); |
// else histoped[nn]->Draw("axis"); |
| 596 |
histoaped[nn]->Draw("same"); |
histoaped[nn]->Draw("same"); |
| 597 |
if((nn+1)%2==1) { |
if((nn+1)%2==1) { |
| 598 |
li.DrawLine(1024.5,2200,1024.5,3200); |
li.DrawLine(1024.5,2200,1024.5,3200); |
| 602 |
li.DrawLine(1024.5,700,1024.5,1700); |
li.DrawLine(1024.5,700,1024.5,1700); |
| 603 |
li.DrawLine(2048.5,700,2048.5,1700); |
li.DrawLine(2048.5,700,2048.5,1700); |
| 604 |
} |
} |
| 605 |
|
c1->Update();//draw pads in canvas |
| 606 |
|
|
| 607 |
|
|
| 608 |
|
|
| 625 |
histosig[nn]->GetYaxis()->CenterTitle(); |
histosig[nn]->GetYaxis()->CenterTitle(); |
| 626 |
histoasig[nn]->SetLineColor(5); |
histoasig[nn]->SetLineColor(5); |
| 627 |
histoasig[nn]->SetLineWidth(1); |
histoasig[nn]->SetLineWidth(1); |
| 628 |
if(ctrk.good0[0]==1 && ctrk.good0[1]==1){ |
// if(ctrk.good0[0]==1 && ctrk.good0[1]==1){ |
| 629 |
histosig[nn]->Draw("b"); |
histosig[nn]->Draw("b"); |
| 630 |
if(nn==1){ |
if(nn==0){ |
| 631 |
|
maxhist=histosig[nn]->GetMaximum(); |
| 632 |
|
b.SetFillColor(107); |
| 633 |
|
b.SetFillStyle(3945); |
| 634 |
|
b.DrawBox(768.,0.,2047.,maxhist); |
| 635 |
|
} |
| 636 |
|
else if(nn==1){ |
| 637 |
maxhist=histosig[nn]->GetMaximum(); |
maxhist=histosig[nn]->GetMaximum(); |
| 638 |
b.SetFillColor(6); |
b.SetFillColor(6); |
| 639 |
b.SetFillStyle(3945); |
b.SetFillStyle(3945); |
| 644 |
b.DrawBox(2816.,0.,2944.,maxhist); |
b.DrawBox(2816.,0.,2944.,maxhist); |
| 645 |
b.DrawBox(2048.,0.,2176.,maxhist); |
b.DrawBox(2048.,0.,2176.,maxhist); |
| 646 |
} |
} |
| 647 |
|
else if(nn==4){ |
| 648 |
|
maxhist=histosig[nn]->GetMaximum(); |
| 649 |
|
b.SetFillColor(107); |
| 650 |
|
b.SetFillStyle(3954); |
| 651 |
|
b.DrawBox(384.,0.,512.,maxhist); |
| 652 |
|
} |
| 653 |
else if(nn==6){ |
else if(nn==6){ |
| 654 |
maxhist=histosig[nn]->GetMaximum(); |
maxhist=histosig[nn]->GetMaximum(); |
| 655 |
b.SetFillColor(6); |
b.SetFillColor(6); |
| 656 |
b.SetFillStyle(3945); |
b.SetFillStyle(3945); |
| 657 |
b.DrawBox(2560.,0.,2816.,maxhist); |
b.DrawBox(2560.,0.,2816.,maxhist); |
| 658 |
b.DrawBox(1024.,0.,1792.,maxhist); |
b.DrawBox(1024.,0.,1535.,maxhist); |
| 659 |
|
|
| 660 |
b.SetFillColor(107); |
b.SetFillColor(107); |
| 661 |
b.SetFillStyle(3954); |
b.SetFillStyle(3954); |
| 662 |
b.DrawBox(512.,0.,768.,maxhist); |
b.DrawBox(512.,0.,768.,maxhist); |
| 663 |
|
b.DrawBox(1536.,0.,1792.,maxhist); |
| 664 |
} |
} |
| 665 |
else if(nn==7){ |
else if(nn==7){ |
| 666 |
maxhist=histosig[nn]->GetMaximum(); |
maxhist=histosig[nn]->GetMaximum(); |
| 668 |
b.SetFillStyle(3954); |
b.SetFillStyle(3954); |
| 669 |
b.DrawBox(512.,0.,768.,maxhist); |
b.DrawBox(512.,0.,768.,maxhist); |
| 670 |
} |
} |
| 671 |
|
else if(nn==9){ |
| 672 |
|
maxhist=histosig[nn]->GetMaximum(); |
| 673 |
|
b.SetFillColor(107); |
| 674 |
|
b.SetFillStyle(3954); |
| 675 |
|
b.DrawBox(256.,0.,384.,maxhist); |
| 676 |
|
} |
| 677 |
else if(nn==11){ |
else if(nn==11){ |
| 678 |
maxhist=histosig[nn]->GetMaximum(); |
maxhist=histosig[nn]->GetMaximum(); |
| 679 |
b.SetFillColor(6); |
b.SetFillColor(6); |
| 683 |
b.SetFillColor(107); |
b.SetFillColor(107); |
| 684 |
b.SetFillStyle(3954); |
b.SetFillStyle(3954); |
| 685 |
b.DrawBox(0.,0.,512.,maxhist); |
b.DrawBox(0.,0.,512.,maxhist); |
| 686 |
b.DrawBox(1920.,0.,2048.,maxhist); |
b.DrawBox(1920.,0.,2560.,maxhist); |
| 687 |
} |
} |
| 688 |
} |
// } |
| 689 |
else histosig[nn]->Draw("axis"); |
// else histosig[nn]->Draw("axis"); |
| 690 |
histoasig[nn]->Draw("same"); |
histoasig[nn]->Draw("same"); |
| 691 |
li.DrawLine(1024.5,0,1024.5,max); |
li.DrawLine(1024.5,0,1024.5,max); |
| 692 |
li.DrawLine(2048.5,0,2048.5,max); |
li.DrawLine(2048.5,0,2048.5,max); |
| 704 |
histosiglad[nn][ii]->GetXaxis()->SetTitle("SIG (ADC channels)"); |
histosiglad[nn][ii]->GetXaxis()->SetTitle("SIG (ADC channels)"); |
| 705 |
histosiglad[nn][ii]->GetXaxis()->CenterTitle(); |
histosiglad[nn][ii]->GetXaxis()->CenterTitle(); |
| 706 |
histosiglad[nn][ii]->Draw(""); |
histosiglad[nn][ii]->Draw(""); |
| 707 |
} |
} |
| 708 |
|
|
| 709 |
};//end loop on views |
};//end loop on views |
| 710 |
c1->Update();//draw pads in canvas |
c1->Update();//draw pads in canvas |
| 711 |
c2->Update();//draw pads in canvas |
c2->Update();//draw pads in canvas |