From c6a8e83808c0e8fd0b26ca68339ff84a3e03478a Mon Sep 17 00:00:00 2001 From: coppedis Date: Tue, 30 Jun 2009 14:16:11 +0000 Subject: [PATCH] Macro for a manual check of the automatic fit on pedestal data --- ZDC/CheckPedestal.C | 196 ++++++++++++++++++++++++++++++++++++++++++++ 1 file changed, 196 insertions(+) create mode 100644 ZDC/CheckPedestal.C diff --git a/ZDC/CheckPedestal.C b/ZDC/CheckPedestal.C new file mode 100644 index 00000000000..d1764ad35a6 --- /dev/null +++ b/ZDC/CheckPedestal.C @@ -0,0 +1,196 @@ +#if !defined(__CINT__) || defined(__MAKECINT__) + +#include +#include +#include +#include "TClassTable.h" +#include +#include +#include +#include +#include +#include +#include +#include +#include +#include + +#endif + +void CheckPedestal( Int_t optPlot = 1, const char * histoFName = "ZDCPedHisto.root") +{ + + int const kNChannels = 24; + // + TFile * file = new TFile(histoFName, "READ"); + file->cd(); + TH1F::AddDirectory(0); + // + TH1F *hPedhg[kNChannels], *hPedOutOfTimehg[kNChannels]; + TH2F *hPedCorrhg[kNChannels]; + TH1F *hPedlg[kNChannels], *hPedOutOfTimelg[kNChannels]; + TH2F *hPedCorrlg[kNChannels]; + // + char namhist1hg[50], namhist2hg[50], namhist3hg[50]; + char namhist1lg[50], namhist2lg[50], namhist3lg[50]; + for(Int_t j=0; j=5 && j<=9){ // ZPC + sprintf(namhist1hg,"PedZPChg_%d",j-5); + sprintf(namhist2hg,"PedZPChgOutOfTime_%d",j-5); + sprintf(namhist3hg,"PedCorrZPChg_%d",j-5); + // + sprintf(namhist1lg,"PedZPClg_%d",j-5); + sprintf(namhist2lg,"PedZPClgOutOfTime_%d",j-5); + sprintf(namhist3lg,"PedCorrZPClg_%d",j-5); + } + else if(j==10 || j==11){ // ZEM + sprintf(namhist1hg,"PedZEMhg_%d",j-10); + sprintf(namhist2hg,"PedZEMhgOutOfTime_%d",j-10); + sprintf(namhist3hg,"PedCorrZEMhg_%d",j-10); + // + sprintf(namhist1lg,"PedZEMlg_%d",j-10); + sprintf(namhist2lg,"PedZEMlgOutOfTime_%d",j-10); + sprintf(namhist3lg,"PedCorrZEMlg_%d",j-10); + } + else if(j>=12 && j<=16){ // ZNA + sprintf(namhist1hg,"PedZNAhg_%d",j-12); + sprintf(namhist2hg,"PedZNAhgOutOfTime_%d",j-12); + sprintf(namhist3hg,"PedCorrZNAhg_%d",j-12); + // + sprintf(namhist1lg,"PedZNAlg_%d",j-12); + sprintf(namhist2lg,"PedZNAlgOutOfTime_%d",j-12); + sprintf(namhist3lg,"PedCorrZNAlg_%d",j-12); + } + else if(j>=17 && j<=21){ // ZPA + sprintf(namhist1hg,"PedZPAhg_%d",j-17); + sprintf(namhist2hg,"PedZPAhgOutOfTime_%d",j-17); + sprintf(namhist3hg,"PedCorrZPAhg_%d",j-17); + // + sprintf(namhist1lg,"PedZPAlg_%d",j-17); + sprintf(namhist2lg,"PedZPAlgOutOfTime_%d",j-17); + sprintf(namhist3lg,"PedCorrZPAlg_%d",j-17); + } + else if(j==22 || j==24){ //Reference PMs + sprintf(namhist1hg,"PedRefhg_%d",j-22); + sprintf(namhist2hg,"PedRefhgOutOfTime_%d",j-22); + sprintf(namhist3hg,"PedCorrRefhg_%d",j-22); + // + sprintf(namhist1lg,"PedReflg_%d",j-22); + sprintf(namhist2lg,"PedReflgOutOfTime_%d",j-22); + sprintf(namhist3lg,"PedCorrReflg_%d",j-22); + } + // --- High gain chain histos + hPedhg[j] = (TH1F*) file->Get(namhist1hg); + hPedOutOfTimehg[j] = (TH1F*) file->Get(namhist2hg); + hPedCorrhg[j] = (TH2F*) file->Get(namhist3hg); + // --- Low gain chain histos + hPedlg[j] = (TH1F*) file->Get(namhist1lg); + hPedOutOfTimelg[j] =(TH1F*) file->Get(namhist2lg); + hPedCorrlg[j] = (TH2F*) file->Get(namhist3lg); + } + + if(optPlot==1){ + // Plot the retrieved histos + //*********************************************************** + // #### ROOT initialization + gROOT->Reset(); + gStyle->SetCanvasColor(10); + gStyle->SetFrameFillColor(10); + gStyle->SetOptTitle(0); + gStyle->SetOptStat(1111); + gStyle->SetOptFit(111); + gStyle->SetTitleTextColor(9); + gStyle->SetStatTextColor(4); + gStyle->SetLineColor(1); + gStyle->SetPalette(1); + //*********************************************************** + TCanvas *c6 = new TCanvas("c6","Side C correlations",0,200,1000,800); + c6->Divide(5,4); + for(Int_t t=0; t<10; t++){ + c6->cd(t+1); + hPedCorrhg[t]->Draw(); + c6->cd(t+11); + hPedCorrlg[t]->Draw(); + } + c6->Print("CorrSideC.ps"); + // + TCanvas *c7 = new TCanvas("c7","Side A correlations",300,200,1000,800); + c7->Divide(5,4); + for(Int_t t=0; t<10; t++){ + c7->cd(t+1); + hPedCorrhg[t+12]->Draw(); + c7->cd(t+11); + hPedCorrlg[t+12]->Draw(); + } + c7->Print("CorrSideA.ps"); + // + TCanvas *c8 = new TCanvas("c8","ZEM correlations",400,200,400,400); + c8->Divide(2,2); + for(Int_t t=0; t<2; t++){ + c8->cd(t+1); + hPedCorrhg[t+10]->Draw(); + c8->cd(t+3); + hPedCorrlg[t+10]->Draw(); + } + c8->Print("CorrZEM.ps"); + //*********************************************************** + TCanvas *c1 = new TCanvas("c1","ZNC pedestals",0,0,1000,400); + c1->Divide(5,2); + for(Int_t y=0; y<5; y++){ + c1->cd(y+1); + hPedhg[y]->Draw(); + c1->cd(y+6); + hPedlg[y]->Draw(); + } + c1->Print("ZNCPed.ps"); + // + TCanvas *c2 = new TCanvas("c2","ZPC pedestals",300,0,1000,400); + c2->Divide(5,2); + for(Int_t y=0; y<5; y++){ + c2->cd(y+1); + hPedhg[y+5]->Draw(); + c2->cd(y+6); + hPedlg[y+5]->Draw(); + } + c2->Print("ZPCPed.ps"); + // + TCanvas *c3 = new TCanvas("c3","ZEM pedestals",400,0,400,400); + c3->Divide(2,2); + for(Int_t y=0; y<2; y++){ + c3->cd(y+1); + hPedhg[y+10]->Draw(); + c3->cd(y+3); + hPedlg[y+10]->Draw(); + } + c3->Print("ZEMPed.ps"); + // + TCanvas *c4 = new TCanvas("c4","ZNA pedestals",0,400,1000,400); + c4->Divide(5,2); + for(Int_t y=0; y<5; y++){ + c4->cd(y+1); + hPedhg[y+12]->Draw(); + c4->cd(y+6); + hPedlg[y+12]->Draw(); + } + c4->Print("ZNAPed.ps"); + // + TCanvas *c5 = new TCanvas("c5","ZPA pedestals",300,400,1000,400); + c5->Divide(5,2); + for(Int_t y=0; y<5; y++){ + c5->cd(y+1); + hPedhg[y+17]->Draw(); + c5->cd(y+6); + hPedlg[y+17]->Draw(); + } + c5->Print("ZPAPed.ps"); + } +} -- 2.43.0