| 2 |
* FTrkCalibQLookExpert.cxx |
* FTrkCalibQLookExpert.cxx |
| 3 |
* |
* |
| 4 |
* autor: D.Fedele |
* autor: D.Fedele |
| 5 |
* version v1r12 |
* version v1r17 |
| 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 |
| 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.,2200.,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); |
| 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 |
|
maxhist=histoped[nn]->GetMaximum(); |
| 579 |
|
b.SetFillColor(107); |
| 580 |
|
b.SetFillStyle(3954); |
| 581 |
|
b.DrawBox(512.,2200.,768.,maxhist); |
| 582 |
|
} |
| 583 |
|
else if(nn==9){ |
| 584 |
|
maxhist=histoped[nn]->GetMaximum(); |
| 585 |
|
b.SetFillColor(107); |
| 586 |
|
b.SetFillStyle(3954); |
| 587 |
|
b.DrawBox(256.,700.,384.,maxhist); |
| 588 |
|
} |
| 589 |
else if(nn==11){ |
else if(nn==11){ |
| 590 |
maxhist=histoped[nn]->GetMaximum(); |
maxhist=histoped[nn]->GetMaximum(); |
| 591 |
b.SetFillColor(6); |
b.SetFillColor(6); |
| 595 |
b.SetFillColor(107); |
b.SetFillColor(107); |
| 596 |
b.SetFillStyle(3954); |
b.SetFillStyle(3954); |
| 597 |
b.DrawBox(0.,700.,512.,maxhist); |
b.DrawBox(0.,700.,512.,maxhist); |
| 598 |
b.DrawBox(1920.,700.,2048.,maxhist); |
b.DrawBox(1920.,700.,2560.,maxhist); |
|
b.DrawBox(2176.,700.,2304.,maxhist); |
|
| 599 |
} |
} |
| 600 |
} |
// } |
| 601 |
else histoped[nn]->Draw("axis"); |
// else histoped[nn]->Draw("axis"); |
| 602 |
histoaped[nn]->Draw("same"); |
histoaped[nn]->Draw("same"); |
| 603 |
if((nn+1)%2==1) { |
if((nn+1)%2==1) { |
| 604 |
li.DrawLine(1024.5,2200,1024.5,3200); |
li.DrawLine(1024.5,2200,1024.5,3200); |
| 608 |
li.DrawLine(1024.5,700,1024.5,1700); |
li.DrawLine(1024.5,700,1024.5,1700); |
| 609 |
li.DrawLine(2048.5,700,2048.5,1700); |
li.DrawLine(2048.5,700,2048.5,1700); |
| 610 |
} |
} |
| 611 |
|
c1->Update();//draw pads in canvas |
| 612 |
|
|
| 613 |
|
|
| 614 |
|
|
| 631 |
histosig[nn]->GetYaxis()->CenterTitle(); |
histosig[nn]->GetYaxis()->CenterTitle(); |
| 632 |
histoasig[nn]->SetLineColor(5); |
histoasig[nn]->SetLineColor(5); |
| 633 |
histoasig[nn]->SetLineWidth(1); |
histoasig[nn]->SetLineWidth(1); |
| 634 |
if(ctrk.good0[0]==1 && ctrk.good0[1]==1){ |
// if(ctrk.good0[0]==1 && ctrk.good0[1]==1){ |
| 635 |
histosig[nn]->Draw("b"); |
histosig[nn]->Draw("b"); |
| 636 |
if(nn==1){ |
if(nn==0){ |
| 637 |
|
maxhist=histosig[nn]->GetMaximum(); |
| 638 |
|
b.SetFillColor(107); |
| 639 |
|
b.SetFillStyle(3945); |
| 640 |
|
b.DrawBox(768.,0.,2047.,maxhist); |
| 641 |
|
} |
| 642 |
|
else if(nn==1){ |
| 643 |
maxhist=histosig[nn]->GetMaximum(); |
maxhist=histosig[nn]->GetMaximum(); |
| 644 |
b.SetFillColor(6); |
b.SetFillColor(6); |
| 645 |
b.SetFillStyle(3945); |
b.SetFillStyle(3945); |
| 674 |
b.SetFillStyle(3954); |
b.SetFillStyle(3954); |
| 675 |
b.DrawBox(512.,0.,768.,maxhist); |
b.DrawBox(512.,0.,768.,maxhist); |
| 676 |
} |
} |
| 677 |
|
else if(nn==8){ |
| 678 |
|
maxhist=histosig[nn]->GetMaximum(); |
| 679 |
|
b.SetFillColor(107); |
| 680 |
|
b.SetFillStyle(3954); |
| 681 |
|
b.DrawBox(512.,0.,768.,maxhist); |
| 682 |
|
} |
| 683 |
|
else if(nn==9){ |
| 684 |
|
maxhist=histosig[nn]->GetMaximum(); |
| 685 |
|
b.SetFillColor(107); |
| 686 |
|
b.SetFillStyle(3954); |
| 687 |
|
b.DrawBox(256.,0.,384.,maxhist); |
| 688 |
|
} |
| 689 |
else if(nn==11){ |
else if(nn==11){ |
| 690 |
maxhist=histosig[nn]->GetMaximum(); |
maxhist=histosig[nn]->GetMaximum(); |
| 691 |
b.SetFillColor(6); |
b.SetFillColor(6); |
| 695 |
b.SetFillColor(107); |
b.SetFillColor(107); |
| 696 |
b.SetFillStyle(3954); |
b.SetFillStyle(3954); |
| 697 |
b.DrawBox(0.,0.,512.,maxhist); |
b.DrawBox(0.,0.,512.,maxhist); |
| 698 |
b.DrawBox(1920.,0.,2048.,maxhist); |
b.DrawBox(1920.,0.,2560.,maxhist); |
|
b.DrawBox(2176.,0.,2304.,maxhist); |
|
| 699 |
} |
} |
| 700 |
} |
// } |
| 701 |
else histosig[nn]->Draw("axis"); |
// else histosig[nn]->Draw("axis"); |
| 702 |
histoasig[nn]->Draw("same"); |
histoasig[nn]->Draw("same"); |
| 703 |
li.DrawLine(1024.5,0,1024.5,max); |
li.DrawLine(1024.5,0,1024.5,max); |
| 704 |
li.DrawLine(2048.5,0,2048.5,max); |
li.DrawLine(2048.5,0,2048.5,max); |
| 716 |
histosiglad[nn][ii]->GetXaxis()->SetTitle("SIG (ADC channels)"); |
histosiglad[nn][ii]->GetXaxis()->SetTitle("SIG (ADC channels)"); |
| 717 |
histosiglad[nn][ii]->GetXaxis()->CenterTitle(); |
histosiglad[nn][ii]->GetXaxis()->CenterTitle(); |
| 718 |
histosiglad[nn][ii]->Draw(""); |
histosiglad[nn][ii]->Draw(""); |
| 719 |
} |
} |
| 720 |
|
|
| 721 |
};//end loop on views |
};//end loop on views |
| 722 |
c1->Update();//draw pads in canvas |
c1->Update();//draw pads in canvas |
| 723 |
c2->Update();//draw pads in canvas |
c2->Update();//draw pads in canvas |