1 #if !defined(__CINT__) || defined(__MAKECINT__)
6 #include <TStopwatch.h>
8 #include <TClassTable.h>
13 #include <TFunction.h>
20 void CheckPedestalHistos(Int_t nRun=0, Int_t optPlot = 1)
23 printf("\n\n YOU MUST PROVIDE A RUN NUMBER!!! \n\n");
27 TGrid::Connect("alien:",0,0,"t");
30 sprintf(histoFName,"alien:///alice/data/2012/Reference/ZDC/%d_pedestalReference.root",nRun);
32 int const kNChannels = 24;
34 TFile *file = TFile::Open(histoFName);
36 TH1F::AddDirectory(0);
38 TH1F *hPedhg[kNChannels], *hPedOutOfTimehg[kNChannels];
39 TH2F *hPedCorrhg[kNChannels];
40 TH1F *hPedlg[kNChannels], *hPedOutOfTimelg[kNChannels];
41 TH2F *hPedCorrlg[kNChannels];
43 char namhist1hg[50], namhist2hg[50], namhist3hg[50];
44 char namhist1lg[50], namhist2lg[50], namhist3lg[50];
45 for(Int_t j=0; j<kNChannels; j++){
47 sprintf(namhist1hg,"PedZNChg_%d",j);
48 sprintf(namhist2hg,"PedZNChgOutOfTime_%d",j);
49 sprintf(namhist3hg,"PedCorrZNChg_%d",j);
51 sprintf(namhist1lg,"PedZNClg_%d",j);
52 sprintf(namhist2lg,"PedZNClgOutOfTime_%d",j);
53 sprintf(namhist3lg,"PedCorrZNClg_%d",j);
55 else if(j>=5 && j<=9){ // ZPC
56 sprintf(namhist1hg,"PedZPChg_%d",j-5);
57 sprintf(namhist2hg,"PedZPChgOutOfTime_%d",j-5);
58 sprintf(namhist3hg,"PedCorrZPChg_%d",j-5);
60 sprintf(namhist1lg,"PedZPClg_%d",j-5);
61 sprintf(namhist2lg,"PedZPClgOutOfTime_%d",j-5);
62 sprintf(namhist3lg,"PedCorrZPClg_%d",j-5);
64 else if(j==10 || j==11){ // ZEM
65 sprintf(namhist1hg,"PedZEMhg_%d",j-9);
66 sprintf(namhist2hg,"PedZEMhgOutOfTime_%d",j-9);
67 sprintf(namhist3hg,"PedCorrZEMhg_%d",j-9);
69 sprintf(namhist1lg,"PedZEMlg_%d",j-9);
70 sprintf(namhist2lg,"PedZEMlgOutOfTime_%d",j-9);
71 sprintf(namhist3lg,"PedCorrZEMlg_%d",j-9);
73 else if(j>=12 && j<=16){ // ZNA
74 sprintf(namhist1hg,"PedZNAhg_%d",j-12);
75 sprintf(namhist2hg,"PedZNAhgOutOfTime_%d",j-12);
76 sprintf(namhist3hg,"PedCorrZNAhg_%d",j-12);
78 sprintf(namhist1lg,"PedZNAlg_%d",j-12);
79 sprintf(namhist2lg,"PedZNAlgOutOfTime_%d",j-12);
80 sprintf(namhist3lg,"PedCorrZNAlg_%d",j-12);
82 else if(j>=17 && j<=21){ // ZPA
83 sprintf(namhist1hg,"PedZPAhg_%d",j-17);
84 sprintf(namhist2hg,"PedZPAhgOutOfTime_%d",j-17);
85 sprintf(namhist3hg,"PedCorrZPAhg_%d",j-17);
87 sprintf(namhist1lg,"PedZPAlg_%d",j-17);
88 sprintf(namhist2lg,"PedZPAlgOutOfTime_%d",j-17);
89 sprintf(namhist3lg,"PedCorrZPAlg_%d",j-17);
91 else if(j==22 || j==23){ //Reference PMs
92 sprintf(namhist1hg,"PedRefhg_%d",j-22);
93 sprintf(namhist2hg,"PedRefhgOutOfTime_%d",j-22);
94 sprintf(namhist3hg,"PedCorrRefhg_%d",j-22);
96 sprintf(namhist1lg,"PedReflg_%d",j-22);
97 sprintf(namhist2lg,"PedReflgOutOfTime_%d",j-22);
98 sprintf(namhist3lg,"PedCorrReflg_%d",j-22);
100 // --- High gain chain histos
101 hPedhg[j] = (TH1F*) file->Get(namhist1hg);
102 hPedOutOfTimehg[j] = (TH1F*) file->Get(namhist2hg);
103 hPedCorrhg[j] = (TH2F*) file->Get(namhist3hg);
104 // --- Low gain chain histos
105 hPedlg[j] = (TH1F*) file->Get(namhist1lg);
106 hPedOutOfTimelg[j] =(TH1F*) file->Get(namhist2lg);
107 hPedCorrlg[j] = (TH2F*) file->Get(namhist3lg);
111 // Plot the retrieved histos
112 //***********************************************************
113 // #### ROOT initialization
115 gStyle->SetCanvasColor(10);
116 gStyle->SetFrameFillColor(10);
117 gStyle->SetOptTitle(0);
118 gStyle->SetOptStat(1111);
119 gStyle->SetOptFit(111);
120 gStyle->SetTitleTextColor(9);
121 gStyle->SetStatTextColor(4);
122 gStyle->SetLineColor(1);
123 gStyle->SetPalette(1);
124 //***********************************************************
125 TCanvas *c6 = new TCanvas("c6","Side C correlations",0,200,1000,800);
127 for(Int_t t=0; t<10; t++){
129 hPedCorrhg[t]->Draw();
131 hPedCorrlg[t]->Draw();
133 c6->Print("CorrSideC.gif");
135 TCanvas *c7 = new TCanvas("c7","Side A correlations",300,200,1000,800);
137 for(Int_t t=0; t<10; t++){
139 hPedCorrhg[t+12]->Draw();
141 hPedCorrlg[t+12]->Draw();
143 c7->Print("CorrSideA.gif");
145 TCanvas *c8 = new TCanvas("c8","ZEM correlations",400,200,400,400);
147 for(Int_t t=0; t<2; t++){
149 hPedCorrhg[t+10]->Draw();
151 hPedCorrlg[t+10]->Draw();
153 c8->Print("CorrZEM.gif");
154 //***********************************************************
155 TCanvas *c1 = new TCanvas("c1","ZNC pedestals",0,0,1000,400);
157 for(Int_t y=0; y<5; y++){
163 c1->Print("ZNCPed.gif");
165 TCanvas *c2 = new TCanvas("c2","ZPC pedestals",300,0,1000,400);
167 for(Int_t y=0; y<5; y++){
173 c2->Print("ZPCPed.gif");
175 TCanvas *c3 = new TCanvas("c3","ZEM pedestals",400,0,400,400);
177 for(Int_t y=0; y<2; y++){
179 hPedhg[y+10]->Draw();
181 hPedlg[y+10]->Draw();
183 c3->Print("ZEMPed.gif");
185 TCanvas *c4 = new TCanvas("c4","ZNA pedestals",0,400,1000,400);
187 for(Int_t y=0; y<5; y++){
189 hPedhg[y+12]->Draw();
191 hPedlg[y+12]->Draw();
193 c4->Print("ZNAPed.gif");
195 TCanvas *c5 = new TCanvas("c5","ZPA pedestals",300,400,1000,400);
197 for(Int_t y=0; y<5; y++){
199 hPedhg[y+17]->Draw();
201 hPedlg[y+17]->Draw();
203 c5->Print("ZPAPed.gif");