2 |
* FTrkCalibQLookExpert.cxx |
* FTrkCalibQLookExpert.cxx |
3 |
* |
* |
4 |
* autor: D.Fedele |
* autor: D.Fedele |
5 |
* version v1r12 |
* version v1r15 |
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==1){ |
537 |
maxhist=histoped[nn]->GetMaximum(); |
maxhist=histoped[nn]->GetMaximum(); |
568 |
b.SetFillStyle(3954); |
b.SetFillStyle(3954); |
569 |
b.DrawBox(512.,700.,768.,maxhist); |
b.DrawBox(512.,700.,768.,maxhist); |
570 |
} |
} |
571 |
|
else if(nn==9){ |
572 |
|
maxhist=histoped[nn]->GetMaximum(); |
573 |
|
b.SetFillColor(107); |
574 |
|
b.SetFillStyle(3954); |
575 |
|
b.DrawBox(256.,700.,384.,maxhist); |
576 |
|
} |
577 |
else if(nn==11){ |
else if(nn==11){ |
578 |
maxhist=histoped[nn]->GetMaximum(); |
maxhist=histoped[nn]->GetMaximum(); |
579 |
b.SetFillColor(6); |
b.SetFillColor(6); |
583 |
b.SetFillColor(107); |
b.SetFillColor(107); |
584 |
b.SetFillStyle(3954); |
b.SetFillStyle(3954); |
585 |
b.DrawBox(0.,700.,512.,maxhist); |
b.DrawBox(0.,700.,512.,maxhist); |
586 |
b.DrawBox(1920.,700.,2048.,maxhist); |
b.DrawBox(1920.,700.,2560.,maxhist); |
|
b.DrawBox(2176.,700.,2304.,maxhist); |
|
587 |
} |
} |
588 |
} |
// } |
589 |
else histoped[nn]->Draw("axis"); |
// else histoped[nn]->Draw("axis"); |
590 |
histoaped[nn]->Draw("same"); |
histoaped[nn]->Draw("same"); |
591 |
if((nn+1)%2==1) { |
if((nn+1)%2==1) { |
592 |
li.DrawLine(1024.5,2200,1024.5,3200); |
li.DrawLine(1024.5,2200,1024.5,3200); |
596 |
li.DrawLine(1024.5,700,1024.5,1700); |
li.DrawLine(1024.5,700,1024.5,1700); |
597 |
li.DrawLine(2048.5,700,2048.5,1700); |
li.DrawLine(2048.5,700,2048.5,1700); |
598 |
} |
} |
599 |
|
c1->Update();//draw pads in canvas |
600 |
|
|
601 |
|
|
602 |
|
|
619 |
histosig[nn]->GetYaxis()->CenterTitle(); |
histosig[nn]->GetYaxis()->CenterTitle(); |
620 |
histoasig[nn]->SetLineColor(5); |
histoasig[nn]->SetLineColor(5); |
621 |
histoasig[nn]->SetLineWidth(1); |
histoasig[nn]->SetLineWidth(1); |
622 |
if(ctrk.good0[0]==1 && ctrk.good0[1]==1){ |
// if(ctrk.good0[0]==1 && ctrk.good0[1]==1){ |
623 |
histosig[nn]->Draw("b"); |
histosig[nn]->Draw("b"); |
624 |
if(nn==1){ |
if(nn==1){ |
625 |
maxhist=histosig[nn]->GetMaximum(); |
maxhist=histosig[nn]->GetMaximum(); |
656 |
b.SetFillStyle(3954); |
b.SetFillStyle(3954); |
657 |
b.DrawBox(512.,0.,768.,maxhist); |
b.DrawBox(512.,0.,768.,maxhist); |
658 |
} |
} |
659 |
|
else if(nn==9){ |
660 |
|
maxhist=histosig[nn]->GetMaximum(); |
661 |
|
b.SetFillColor(107); |
662 |
|
b.SetFillStyle(3954); |
663 |
|
b.DrawBox(256.,0.,384.,maxhist); |
664 |
|
} |
665 |
else if(nn==11){ |
else if(nn==11){ |
666 |
maxhist=histosig[nn]->GetMaximum(); |
maxhist=histosig[nn]->GetMaximum(); |
667 |
b.SetFillColor(6); |
b.SetFillColor(6); |
671 |
b.SetFillColor(107); |
b.SetFillColor(107); |
672 |
b.SetFillStyle(3954); |
b.SetFillStyle(3954); |
673 |
b.DrawBox(0.,0.,512.,maxhist); |
b.DrawBox(0.,0.,512.,maxhist); |
674 |
b.DrawBox(1920.,0.,2048.,maxhist); |
b.DrawBox(1920.,0.,2560.,maxhist); |
|
b.DrawBox(2176.,0.,2304.,maxhist); |
|
675 |
} |
} |
676 |
} |
// } |
677 |
else histosig[nn]->Draw("axis"); |
// else histosig[nn]->Draw("axis"); |
678 |
histoasig[nn]->Draw("same"); |
histoasig[nn]->Draw("same"); |
679 |
li.DrawLine(1024.5,0,1024.5,max); |
li.DrawLine(1024.5,0,1024.5,max); |
680 |
li.DrawLine(2048.5,0,2048.5,max); |
li.DrawLine(2048.5,0,2048.5,max); |
692 |
histosiglad[nn][ii]->GetXaxis()->SetTitle("SIG (ADC channels)"); |
histosiglad[nn][ii]->GetXaxis()->SetTitle("SIG (ADC channels)"); |
693 |
histosiglad[nn][ii]->GetXaxis()->CenterTitle(); |
histosiglad[nn][ii]->GetXaxis()->CenterTitle(); |
694 |
histosiglad[nn][ii]->Draw(""); |
histosiglad[nn][ii]->Draw(""); |
695 |
} |
} |
696 |
|
|
697 |
};//end loop on views |
};//end loop on views |
698 |
c1->Update();//draw pads in canvas |
c1->Update();//draw pads in canvas |
699 |
c2->Update();//draw pads in canvas |
c2->Update();//draw pads in canvas |