1 #if !defined(__CINT__) || defined(__MAKECINT__)
6 #include "TClassTable.h"
7 #include <TStopwatch.h>
9 #include <TClassTable.h>
14 #include <TFunction.h>
20 void CheckLaser(const char * histoFName = "ZDCLaserHisto.root", Int_t optPlot = 1)
23 TFile * file = new TFile(histoFName, "READ");
25 TH1F::AddDirectory(0);
27 // --- Histos for reference PMTs (high gain chains)
28 TH1F *hPMRefChg = new TH1F("hPMRefChg","hPMRefChg", 100,0.,1000.);
29 TH1F *hPMRefAhg = new TH1F("hPMRefAhg","hPMRefAhg", 100,0.,1000.);
30 TH1F *hPMRefClg = new TH1F("hPMRefClg","hPMRefClg", 100,0.,4000.);
31 TH1F *hPMRefAlg = new TH1F("hPMRefAlg","hPMRefAlg", 100,0.,4000.);
33 hPMRefChg = (TH1F*) file->Get("hPMRefChg");
34 hPMRefAhg = (TH1F*) file->Get("hPMRefAhg");
35 hPMRefClg = (TH1F*) file->Get("hPMRefClg");
36 hPMRefAlg = (TH1F*) file->Get("hPMRefAlg");
37 // --- Histos for detector PMTs
38 TH1F *hZNChg[5], *hZPChg[5], *hZNAhg[5], *hZPAhg[5], *hZEMhg[2];
39 TH1F *hZNClg[5], *hZPClg[5], *hZNAlg[5], *hZPAlg[5], *hZEMlg[2];
40 char hnamZNChg[20], hnamZPChg[20], hnamZNAhg[20], hnamZPAhg[20];
41 char hnamZNClg[20], hnamZPClg[20], hnamZNAlg[20], hnamZPAlg[20];
42 char hnamZEMhg[20], hnamZEMlg[20];
43 for(Int_t j=0; j<5; j++){
44 sprintf(hnamZNChg,"ZNChg-tow%d",j);
45 sprintf(hnamZPChg,"ZPChg-tow%d",j);
46 sprintf(hnamZNAhg,"ZNAhg-tow%d",j);
47 sprintf(hnamZPAhg,"ZPAhg-tow%d",j);
49 hZNChg[j] = new TH1F(hnamZNChg, hnamZNChg, 100, 0., 1000.);
50 hZPChg[j] = new TH1F(hnamZPChg, hnamZPChg, 100, 0., 1000.);
51 hZNAhg[j] = new TH1F(hnamZNAhg, hnamZNAhg, 100, 0., 1000.);
52 hZPAhg[j] = new TH1F(hnamZPAhg, hnamZPAhg, 100, 0., 1000.);
54 sprintf(hnamZNClg,"ZNClg-tow%d",j);
55 sprintf(hnamZPClg,"ZPClg-tow%d",j);
56 sprintf(hnamZNAlg,"ZNAlg-tow%d",j);
57 sprintf(hnamZPAlg,"ZPAlg-tow%d",j);
59 hZNClg[j] = new TH1F(hnamZNClg, hnamZNClg, 100, 0., 4000.);
60 hZPClg[j] = new TH1F(hnamZPClg, hnamZPClg, 100, 0., 4000.);
61 hZNAlg[j] = new TH1F(hnamZNAlg, hnamZNAlg, 100, 0., 4000.);
62 hZPAlg[j] = new TH1F(hnamZPAlg, hnamZPAlg, 100, 0., 4000.);
65 sprintf(hnamZEMhg,"ZEM%dhg",j);
66 sprintf(hnamZEMlg,"ZEM%dlg",j);
68 hZEMhg[j] = new TH1F(hnamZEMhg, hnamZEMhg, 100, 0., 1000.);
69 hZEMlg[j] = new TH1F(hnamZEMlg, hnamZEMlg, 100, 0., 4000.);
72 hZNChg[j] = (TH1F*) file->Get(hnamZNChg);
73 hZPChg[j] = (TH1F*) file->Get(hnamZPChg);
74 hZNAhg[j] = (TH1F*) file->Get(hnamZNAhg);
75 hZPAhg[j] = (TH1F*) file->Get(hnamZPAhg);
77 hZNClg[j] = (TH1F*) file->Get(hnamZNClg);
78 hZPClg[j] = (TH1F*) file->Get(hnamZPClg);
79 hZNAlg[j] = (TH1F*) file->Get(hnamZNAlg);
80 hZPAlg[j] = (TH1F*) file->Get(hnamZPAlg);
86 // Plot the retrieved histos
87 //***********************************************************
88 // #### ROOT initialization
90 gStyle->SetCanvasColor(10);
91 gStyle->SetFrameFillColor(10);
92 gStyle->SetOptTitle(0);
93 gStyle->SetOptStat(1111);
94 gStyle->SetOptFit(111);
95 gStyle->SetTitleTextColor(9);
96 gStyle->SetStatTextColor(4);
97 gStyle->SetLineColor(1);
98 gStyle->SetPalette(1);
99 //***********************************************************
100 TCanvas *c1 = new TCanvas("c1","ZNC",0,0,1000,400);
102 for(Int_t y=0; y<5; y++){
108 c1->Print("ZNCLaser.ps");
110 TCanvas *c2 = new TCanvas("c2","ZPC",300,0,1000,400);
112 for(Int_t y=0; y<5; y++){
118 c2->Print("ZPCLaser.ps");
120 TCanvas *c3 = new TCanvas("c3","ZEM",400,0,400,400);
122 for(Int_t y=0; y<2; y++){
128 c3->Print("ZEMLaser.ps");
130 TCanvas *c4 = new TCanvas("c4","ZNA",0,400,1000,400);
132 for(Int_t y=0; y<5; y++){
138 c4->Print("ZNALaser.ps");
140 TCanvas *c5 = new TCanvas("c5","ZPA",300,400,1000,400);
142 for(Int_t y=0; y<5; y++){
148 c5->Print("ZPALaser.ps");
150 TCanvas *c6 = new TCanvas("c6","Ref",400,0,400,400);
161 c6->Print("RefLaser.ps");