2 |
* FTrkCalibQLookExpert.cxx |
* FTrkCalibQLookExpert.cxx |
3 |
* |
* |
4 |
* autor: D.Fedele |
* autor: D.Fedele |
5 |
* version v1r12 |
* 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 |
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); |
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); |
|
b.DrawBox(2176.,700.,2304.,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); |
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); |
|
b.DrawBox(2176.,0.,2304.,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 |