1 #if !defined(__CINT__) || defined(__MAKECINT__)
6 #include <TStopwatch.h>
8 #include <TClassTable.h>
13 #include <TFunction.h>
20 void CheckLaserHistos(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/2010/Reference/ZDC/%d_laserReference.root",nRun);
32 TFile *file = TFile::Open(histoFName);
34 TH1F::AddDirectory(0);
36 // --- Histos for reference PMTs (high gain chains)
37 TH1F *hPMRefChg = new TH1F("hPMRefChg","hPMRefChg", 100,0.,1000.);
38 TH1F *hPMRefAhg = new TH1F("hPMRefAhg","hPMRefAhg", 100,0.,1000.);
39 TH1F *hPMRefClg = new TH1F("hPMRefClg","hPMRefClg", 100,0.,4000.);
40 TH1F *hPMRefAlg = new TH1F("hPMRefAlg","hPMRefAlg", 100,0.,4000.);
42 hPMRefChg = (TH1F*) file->Get("hPMRefChg");
43 hPMRefAhg = (TH1F*) file->Get("hPMRefAhg");
44 hPMRefClg = (TH1F*) file->Get("hPMRefClg");
45 hPMRefAlg = (TH1F*) file->Get("hPMRefAlg");
46 // --- Histos for detector PMTs
47 TH1F *hZNChg[5], *hZPChg[5], *hZNAhg[5], *hZPAhg[5], *hZEMhg[2];
48 TH1F *hZNClg[5], *hZPClg[5], *hZNAlg[5], *hZPAlg[5], *hZEMlg[2];
49 char hnamZNChg[20], hnamZPChg[20], hnamZNAhg[20], hnamZPAhg[20];
50 char hnamZNClg[20], hnamZPClg[20], hnamZNAlg[20], hnamZPAlg[20];
51 char hnamZEMhg[20], hnamZEMlg[20];
52 for(Int_t j=0; j<5; j++){
53 sprintf(hnamZNChg,"ZNChg-tow%d",j);
54 sprintf(hnamZPChg,"ZPChg-tow%d",j);
55 sprintf(hnamZNAhg,"ZNAhg-tow%d",j);
56 sprintf(hnamZPAhg,"ZPAhg-tow%d",j);
58 hZNChg[j] = new TH1F(hnamZNChg, hnamZNChg, 100, 0., 1000.);
59 hZPChg[j] = new TH1F(hnamZPChg, hnamZPChg, 100, 0., 1000.);
60 hZNAhg[j] = new TH1F(hnamZNAhg, hnamZNAhg, 100, 0., 1000.);
61 hZPAhg[j] = new TH1F(hnamZPAhg, hnamZPAhg, 100, 0., 1000.);
63 sprintf(hnamZNClg,"ZNClg-tow%d",j);
64 sprintf(hnamZPClg,"ZPClg-tow%d",j);
65 sprintf(hnamZNAlg,"ZNAlg-tow%d",j);
66 sprintf(hnamZPAlg,"ZPAlg-tow%d",j);
68 hZNClg[j] = new TH1F(hnamZNClg, hnamZNClg, 100, 0., 4000.);
69 hZPClg[j] = new TH1F(hnamZPClg, hnamZPClg, 100, 0., 4000.);
70 hZNAlg[j] = new TH1F(hnamZNAlg, hnamZNAlg, 100, 0., 4000.);
71 hZPAlg[j] = new TH1F(hnamZPAlg, hnamZPAlg, 100, 0., 4000.);
74 sprintf(hnamZEMhg,"ZEM%dhg",j);
75 sprintf(hnamZEMlg,"ZEM%dlg",j);
77 hZEMhg[j] = new TH1F(hnamZEMhg, hnamZEMhg, 100, 0., 1000.);
78 hZEMlg[j] = new TH1F(hnamZEMlg, hnamZEMlg, 100, 0., 4000.);
81 hZNChg[j] = (TH1F*) file->Get(hnamZNChg);
82 hZPChg[j] = (TH1F*) file->Get(hnamZPChg);
83 hZNAhg[j] = (TH1F*) file->Get(hnamZNAhg);
84 hZPAhg[j] = (TH1F*) file->Get(hnamZPAhg);
86 hZNClg[j] = (TH1F*) file->Get(hnamZNClg);
87 hZPClg[j] = (TH1F*) file->Get(hnamZPClg);
88 hZNAlg[j] = (TH1F*) file->Get(hnamZNAlg);
89 hZPAlg[j] = (TH1F*) file->Get(hnamZPAlg);
95 // Plot the retrieved histos
96 //***********************************************************
97 // #### ROOT initialization
99 gStyle->SetCanvasColor(10);
100 gStyle->SetFrameFillColor(10);
101 gStyle->SetOptTitle(0);
102 gStyle->SetOptStat(1111);
103 gStyle->SetOptFit(111);
104 gStyle->SetTitleTextColor(9);
105 gStyle->SetStatTextColor(4);
106 gStyle->SetLineColor(1);
107 gStyle->SetPalette(1);
108 //***********************************************************
109 TCanvas *c1 = new TCanvas("c1","ZNC",0,0,1000,400);
111 for(Int_t y=0; y<5; y++){
117 c1->Print("ZNCLaser.ps");
119 TCanvas *c2 = new TCanvas("c2","ZPC",300,0,1000,400);
121 for(Int_t y=0; y<5; y++){
127 c2->Print("ZPCLaser.ps");
129 TCanvas *c3 = new TCanvas("c3","ZEM",400,0,400,400);
131 for(Int_t y=0; y<2; y++){
137 c3->Print("ZEMLaser.ps");
139 TCanvas *c4 = new TCanvas("c4","ZNA",0,400,1000,400);
141 for(Int_t y=0; y<5; y++){
147 c4->Print("ZNALaser.ps");
149 TCanvas *c5 = new TCanvas("c5","ZPA",300,400,1000,400);
151 for(Int_t y=0; y<5; y++){
157 c5->Print("ZPALaser.ps");
159 TCanvas *c6 = new TCanvas("c6","Ref",400,0,400,400);
170 c6->Print("RefLaser.ps");