5 {
6 gROOT->SetStyle("Pub");
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
30Int_t
plateArray[13]={10,12,14,16,18,20,23,26,29,32,35,40,45};
32
33
34Double_t E[12] = {0.5,1.0,2.0, 4.0, 0.75, 1.5, 3.0, 6.0, 8.0, 16.0, 32.0 ,64.0};
35Int_t E_ASCEND[12] = {0,4,1,5,2,6,3,7,8,9,10,11};
36Int_t NPL[13]= {10,12,14,16,18,20,23,26,29,32,35,40,45};
37
38Int_t N_E=12; N_E=12;
39Int_t N_NPL=13; N_NPL=13;
40
41TGraph* graph_MeanEnergy__vs__Energy[13];
42TGraph* graph_MeanEnergy__vs__Npl[12];
43TGraph* graph_RMSEnergy__vs__Energy[13];
44TGraph* graph_RMSEnergy__vs__Npl[12];
45
46TGraph* graph_MeanEnergyOverE__vs__Energy[13];
47TGraph* graph_MeanEnergyOverE__vs__Npl[12];
48TGraph* graph_RMSEnergyOverE__vs__Energy[13];
49TGraph* graph_RMSEnergyOverE__vs__Npl[12];
50
51TGraph* graph_MeanEnergyRelativeOverE__vs__Energy[13];
52TGraph* graph_MeanEnergyRelativeOverE__vs__Npl[12];
53TGraph* graph_RMSEnergyRelativeOverE__vs__Energy[13];
54TGraph* graph_RMSEnergyRelativeOverE__vs__Npl[12];
55
56TGraph* graph_MeanEnergyMinusNominalE__vs__Energy[13];
57TGraph* graph_MeanEnergyMinusNominalE__vs__Npl[12];
58
59for (int loop_E=0; loop_E<N_E; loop_E++) {
60 graph_MeanEnergy__vs__Npl[loop_E] = new TGraph(N_NPL);
61 graph_RMSEnergy__vs__Npl[loop_E] = new TGraph(N_NPL);
62 graph_MeanEnergyOverE__vs__Npl[loop_E] = new TGraph(N_NPL);
63 graph_RMSEnergyOverE__vs__Npl[loop_E] = new TGraph(N_NPL);
64 graph_MeanEnergyRelativeOverE__vs__Npl[loop_E] = new TGraph(N_NPL);
65 graph_RMSEnergyRelativeOverE__vs__Npl[loop_E] = new TGraph(N_NPL);
66 graph_MeanEnergyMinusNominalE__vs__Npl[loop_E] = new TGraph(N_NPL);
67}
68for (int loop_NPL=0; loop_NPL<N_NPL; loop_NPL++) {
69 graph_MeanEnergy__vs__Energy[loop_NPL] = new TGraph(N_E);
70 graph_RMSEnergy__vs__Energy[loop_NPL] = new TGraph(N_E);
71 graph_MeanEnergyOverE__vs__Energy[loop_NPL] = new TGraph(N_E);
72 graph_RMSEnergyOverE__vs__Energy[loop_NPL] = new TGraph(N_E);
73 graph_MeanEnergyRelativeOverE__vs__Energy[loop_NPL] = new TGraph(N_E);
74 graph_RMSEnergyRelativeOverE__vs__Energy[loop_NPL] = new TGraph(N_E);
75 graph_MeanEnergyMinusNominalE__vs__Energy[loop_NPL] = new TGraph(N_E);
76}
77
78for (int loop_NPL=0; loop_NPL<N_NPL; loop_NPL++) {
79 graph_MeanEnergy__vs__Energy[loop_NPL]->SetTitle("MeanEnergy__vs__Energy");
80 graph_RMSEnergy__vs__Energy[loop_NPL]->SetTitle("RMSEnergy__vs__Energy");
81
82 graph_MeanEnergyOverE__vs__Energy[loop_NPL]->SetTitle("graph_MeanEnergyOverE__vs__Energy");
83 graph_RMSEnergyOverE__vs__Energy[loop_NPL]->SetTitle("graph_RMSEnergyOverE__vs__Energy");
84
85 graph_MeanEnergyRelativeOverE__vs__Energy[loop_NPL]->SetTitle("graph_MeanEnergyRelativeOverE__vs__Energy");
86 graph_RMSEnergyRelativeOverE__vs__Energy[loop_NPL]->SetTitle("graph_RMSEnergyRelativeOverE__vs__Energy");
87
88 graph_MeanEnergyMinusNominalE__vs__Energy[loop_NPL]->SetTitle("graph_MeanEnergyMinusNominalE__vs__Energy");
89}
90
91
92
93
95 const char*
name =
"libShower_Energy_Systematics_Electron.txt";
97 for (int i=0;i<13;i++) {
98 loop_NPL=i;
99 cout << loop_NPL << "---------------------------------" << endl;
101 for (int j=0;j<12;j++) {
103
104 graph_RMSEnergy__vs__Energy[loop_NPL]->SetPoint(j,E[E_ASCEND[j]],
dummyf);
105 graph_RMSEnergyRelativeOverE__vs__Energy[loop_NPL]->SetPoint(j,E[E_ASCEND[j]],
dummyf/E[E_ASCEND[j]]);
106
107 }
108 }
111
112
113
114
115
116
117
118
119TCanvas* Canv_RMSEnergyOverE__vs__Energy =
new TCanvas(
"Canv_RMSEnergyOverE__vs__Energy_SYSTEMATICS_ALLCOMBINED",
"Canv_RMSEnergyOverE__vs__Energy_SYSTEMATICS_ALLCOMBINED",1050,750);
120
121
122graph_RMSEnergyRelativeOverE__vs__Energy[5]->Draw("APLS");
123graph_RMSEnergyRelativeOverE__vs__Energy[5]->SetTitle("Electron");
124graph_RMSEnergyRelativeOverE__vs__Energy[5]->GetXaxis()->SetTitle("Energy [GeV]");
125graph_RMSEnergyRelativeOverE__vs__Energy[5]->GetYaxis()->SetTitle("#sigma (E) / E (\% /100) Syst. Result");
126graph_RMSEnergyRelativeOverE__vs__Energy[5]->GetYaxis()->SetTitleOffset(1.2);
127graph_RMSEnergyRelativeOverE__vs__Energy[5]->GetXaxis()->SetTitleColor(1);
128graph_RMSEnergyRelativeOverE__vs__Energy[5]->GetYaxis()->SetTitleColor(1);
129graph_RMSEnergyRelativeOverE__vs__Energy[5]->GetXaxis()->CenterTitle();
130graph_RMSEnergyRelativeOverE__vs__Energy[5]->GetYaxis()->CenterTitle();
131
132graph_RMSEnergyRelativeOverE__vs__Energy[5]->GetXaxis()->SetRangeUser(0,17);
133
134
135
136graph_RMSEnergyRelativeOverE__vs__Energy[0]->Draw("LSsame");
137graph_RMSEnergyRelativeOverE__vs__Energy[0]->SetLineStyle(2);
138
139
140graph_RMSEnergyRelativeOverE__vs__Energy[9]->Draw("LSsame");
141graph_RMSEnergyRelativeOverE__vs__Energy[9]->SetLineStyle(5);
142
143
144graph_RMSEnergyRelativeOverE__vs__Energy[11]->Draw("LSsame");
145graph_RMSEnergyRelativeOverE__vs__Energy[11]->SetLineStyle(3);
146
147
148TLegend* leg_selected=
new TLegend(0.65,0.55,0.95,0.85,
NULL,
"brNDC");
149leg_selected->AddEntry(graph_RMSEnergyRelativeOverE__vs__Energy[0],"NPl=10","L");
150leg_selected->AddEntry(graph_RMSEnergyRelativeOverE__vs__Energy[5],"NPl=20","L");
151leg_selected->AddEntry(graph_RMSEnergyRelativeOverE__vs__Energy[9],"NPl=32","L");
152leg_selected->AddEntry(graph_RMSEnergyRelativeOverE__vs__Energy[11],"NPl=40","L");
153leg_selected->SetFillStyle(0);
154leg_selected->SetBorderSize(0);
155leg_selected->Draw();
156
157TPaveText*
pt =
new TPaveText(0.8,0.93,0.98,0.98,
"brNDC");
164
165
166TString canvname;
167if (
PARTICLE_TYPE==0) canvname=TString(Form(
"Canv_SYSTEMATICS_ALLCOMBINED__RMSEnergyOverE__vs__Energy__ELECTRON",0));
168if (
PARTICLE_TYPE==1) canvname=TString(Form(
"Canv_SYSTEMATICS_ALLCOMBINED__RMSEnergyOverE__vs__Energy__GAMMA",0));
169Canv_RMSEnergyOverE__vs__Energy->Print(canvname+".C");
170Canv_RMSEnergyOverE__vs__Energy->Print(canvname+".pdf");
171Canv_RMSEnergyOverE__vs__Energy->Print(canvname+".eps");
172Canv_RMSEnergyOverE__vs__Energy->Print(canvname+".png");
173
174
175
176TCanvas* Canv_RMSEnergy__vs__Energy =
new TCanvas(
"Canv_RMSEnergy__vs__Energy",
"Canv_RMSEnergy__vs__Energy",1050,750);
177
178
179graph_RMSEnergy__vs__Energy[5]->Draw("APLS");
180graph_RMSEnergy__vs__Energy[5]->SetTitle("Electron");
181graph_RMSEnergy__vs__Energy[5]->GetXaxis()->SetTitle("Energy [GeV]");
182graph_RMSEnergy__vs__Energy[5]->GetYaxis()->SetTitle("#sigma (E) [GeV] Syst. Result");
183graph_RMSEnergy__vs__Energy[5]->GetYaxis()->SetTitleOffset(1.2);
184graph_RMSEnergy__vs__Energy[5]->GetXaxis()->SetTitleColor(1);
185graph_RMSEnergy__vs__Energy[5]->GetYaxis()->SetTitleColor(1);
186graph_RMSEnergy__vs__Energy[5]->GetXaxis()->CenterTitle();
187graph_RMSEnergy__vs__Energy[5]->GetYaxis()->CenterTitle();
188
189graph_RMSEnergy__vs__Energy[5]->GetXaxis()->SetRangeUser(0,17);
190
191
192graph_RMSEnergy__vs__Energy[0]->Draw("LSsame");
193graph_RMSEnergy__vs__Energy[0]->SetLineStyle(2);
194
195
196graph_RMSEnergy__vs__Energy[9]->Draw("LSsame");
197graph_RMSEnergy__vs__Energy[9]->SetLineStyle(5);
198
199
200graph_RMSEnergy__vs__Energy[11]->Draw("LSsame");
201graph_RMSEnergy__vs__Energy[11]->SetLineStyle(3);
202
203
204TLegend* leg_selected=
new TLegend(0.18,0.55,0.5,0.85,
NULL,
"brNDC");
205leg_selected->AddEntry(graph_RMSEnergy__vs__Energy[0],"NPl=10","L");
206leg_selected->AddEntry(graph_RMSEnergy__vs__Energy[5],"NPl=20","L");
207leg_selected->AddEntry(graph_RMSEnergy__vs__Energy[9],"NPl=32","L");
208leg_selected->AddEntry(graph_RMSEnergy__vs__Energy[11],"NPl=40","L");
209leg_selected->SetFillStyle(0);
210leg_selected->SetBorderSize(0);
211leg_selected->Draw();
212
214
215if (
PARTICLE_TYPE==0) canvname=TString(Form(
"Canv_SYSTEMATICS_ALLCOMBINED__RMSEnergy__vs__Energy__ELECTRON",0));
216if (
PARTICLE_TYPE==1) canvname=TString(Form(
"Canv_SYSTEMATICS_ALLCOMBINED__RMSEnergy__vs__Energy__GAMMA",0));
217Canv_RMSEnergy__vs__Energy->Print(canvname+".C");
218Canv_RMSEnergy__vs__Energy->Print(canvname+".pdf");
219Canv_RMSEnergy__vs__Energy->Print(canvname+".eps");
220Canv_RMSEnergy__vs__Energy->Print(canvname+".png");
221
222
223
224
225
226
227
228
229
230
231
232
233
234
235Double_t x,y;
236TString latexname;
238
239
240
241if (
PARTICLE_TYPE==0) latexname=TString(Form(
"table__SYSTEMATICS_ALLCOMBINED__RMSEnergyRelativeOverE__vs__Energy__ELECTRON_ALLPLATES.tex",
NULL));
242if (
PARTICLE_TYPE==1) latexname=TString(Form(
"table__SYSTEMATICS_ALLCOMBINED__RMSEnergyRelativeOverE__vs__Energy__GAMMA_ALLPLATES.tex",
NULL));
243cout << endl << endl << latexname << endl << endl;
244
247
248for (
int p=0;
p<80;
p++) cout <<
"-";
249cout << endl;
250for (
int p=0;
p<80;
p++) cout <<
"-";
251cout << endl;
252cout << setw(6) << setprecision(2) << " E\Npl ";
253for (
int p=0;
p<N_E;
p++) cout << setw(6) << setprecision(2) << E[E_ASCEND[
p]] <<
" ";
254cout << endl;
255for (
int o=0;
o<N_NPL;
o++) {
256cout << setw(6) << setprecision(2) << NPL[
o];
257for (
int p=0;
p<N_E;
p++) {
258 graph_RMSEnergyRelativeOverE__vs__Energy[
o]->GetPoint(
p,x,y);
259 cout << setw(6) << setprecision(2) << y << " ";
260}
261cout << endl;
262}
263for (
int p=0;
p<80;
p++) cout <<
"-";
264cout << endl;
265for (
int p=0;
p<80;
p++) cout <<
"-";
266cout << endl;
267
268cout << "\\begin{tabular}{c c c c c c c c c c c c c}" << endl;
269cout << "\\hline" << endl;
270cout << "\\hline" << endl;
271cout << " Energy/NPl & ";
272for (
int p=0;
p<N_E-1;
p++) cout << setw(6) << setprecision(2) << E[E_ASCEND[
p]] <<
" & ";
273cout << setw(6) << setprecision(2) << E[E_ASCEND[N_E-1]] << " ";
274cout << "\\\\ " << endl;
275cout << "\\hline" << endl;
276for (
int o=0;
o<N_NPL;
o++) {
277cout << setw(6) << setprecision(2) << NPL[
o] <<
" ";
278for (
int p=0;
p<N_E;
p++) {
279 cout << " & ";
280 graph_RMSEnergyRelativeOverE__vs__Energy[
o]->GetPoint(
p,x,y);
281 cout << setw(6) << setprecision(2) << y << " ";
282}
283cout << " \\\\" << endl;
284}
285cout << "\\hline" << endl;
286cout << "\\hline" << endl;
287cout << "\\end{tabular}" << endl;
288
289
290stream << "\\begin{tabular}{c c c c c c c c c c c c c}" << endl;
291stream << "\\hline" << endl;
292stream << "\\hline" << endl;
293stream << " Energy/NPl & ";
294for (
int p=0;
p<N_E-1;
p++) stream << setw(6) << setprecision(2) << E[E_ASCEND[
p]] <<
" & ";
295stream << setw(6) << setprecision(2) << E[E_ASCEND[N_E-1]] << " ";
296stream << "\\\\ " << endl;
297stream << "\\hline" << endl;
298for (
int o=0;
o<N_NPL;
o++) {
299stream << setw(6) << setprecision(2) << NPL[
o] <<
" ";
300for (
int p=0;
p<N_E;
p++) {
301 stream << " & ";
302 graph_RMSEnergyRelativeOverE__vs__Energy[
o]->GetPoint(
p,x,y);
303 stream << setw(6) << setprecision(2) << y << " ";
304}
305stream << " \\\\" << endl;
306}
307stream << "\\hline" << endl;
308stream << "\\hline" << endl;
309stream << "\\end{tabular}" << endl;
310
311
312
313graph_RMSEnergyRelativeOverE__vs__Energy[0]->Print();
314
315
316
317
318
319if (
PARTICLE_TYPE==0) latexname=TString(Form(
"table__SYSTEMATICS_ALLCOMBINED__RMSEnergy__vs__Energy__ELECTRON_ALLPLATES.tex",
NULL));
320if (
PARTICLE_TYPE==1) latexname=TString(Form(
"table__SYSTEMATICS_ALLCOMBINED__RMSEnergy__vs__Energy__GAMMA_ALLPLATES.tex",
NULL));
321cout << endl << endl << latexname << endl << endl;
324
325for (
int p=0;
p<80;
p++) cout <<
"-";
326cout << endl;
327for (
int p=0;
p<80;
p++) cout <<
"-";
328cout << endl;
329cout << setw(6) << setprecision(2) << " E\Npl ";
330for (
int p=0;
p<N_E;
p++) cout << setw(6) << setprecision(2) << E[E_ASCEND[
p]] <<
" ";
331cout << endl;
332for (
int o=0;
o<N_NPL;
o++) {
333cout << setw(6) << setprecision(2) << NPL[
o];
334for (
int p=0;
p<N_E;
p++) {
335 graph_RMSEnergy__vs__Energy[
o]->GetPoint(
p,x,y);
336 cout << setw(6) << setprecision(2) << y << " ";
337}
338cout << endl;
339}
340for (
int p=0;
p<80;
p++) cout <<
"-";
341cout << endl;
342for (
int p=0;
p<80;
p++) cout <<
"-";
343cout << endl;
344
345cout << "\\begin{tabular}{c c c c c c c c c c c c c}" << endl;
346cout << "\\hline" << endl;
347cout << "\\hline" << endl;
348cout << " Energy/NPl & ";
349for (
int p=0;
p<N_E-1;
p++) cout << setw(6) << setprecision(2) << E[E_ASCEND[
p]] <<
" & ";
350cout << setw(6) << setprecision(2) << E[E_ASCEND[N_E-1]] << " ";
351cout << "\\\\ " << endl;
352cout << "\\hline" << endl;
353for (
int o=0;
o<N_NPL;
o++) {
354cout << setw(6) << setprecision(2) << NPL[
o] <<
" ";
355for (
int p=0;
p<N_E;
p++) {
356 cout << " & ";
357 graph_RMSEnergy__vs__Energy[
o]->GetPoint(
p,x,y);
358 cout << setw(6) << setprecision(2) << y << " ";
359}
360cout << " \\\\" << endl;
361}
362cout << "\\hline" << endl;
363cout << "\\hline" << endl;
364cout << "\\end{tabular}" << endl;
365
366
367stream << "\\begin{tabular}{c c c c c c c c c c c c c}" << endl;
368stream << "\\hline" << endl;
369stream << "\\hline" << endl;
370stream << " Energy/NPl & ";
371for (
int p=0;
p<N_E-1;
p++) stream << setw(6) << setprecision(2) << E[E_ASCEND[
p]] <<
" & ";
372stream << setw(6) << setprecision(2) << E[E_ASCEND[N_E-1]] << " ";
373stream << "\\\\ " << endl;
374stream << "\\hline" << endl;
375for (
int o=0;
o<N_NPL;
o++) {
376stream << setw(6) << setprecision(2) << NPL[
o] <<
" ";
377for (
int p=0;
p<N_E;
p++) {
378 stream << " & ";
379 graph_RMSEnergy__vs__Energy[
o]->GetPoint(
p,x,y);
380 stream << setw(6) << setprecision(2) << y << " ";
381}
382stream << " \\\\" << endl;
383}
384stream << "\\hline" << endl;
385stream << "\\hline" << endl;
386stream << "\\end{tabular}" << endl;
387
388
389
390gSystem->Exec("cp -v table__SYSTEMATICS_ALLCOMBINED__RMSEnergyRelativeOverE__vs__Energy__ELECTRON_ALLPLATES.tex /home/meisel/opus/tables/.");
391gSystem->Exec("cp -v table__SYSTEMATICS_ALLCOMBINED__RMSEnergy__vs__Energy__ELECTRON_ALLPLATES.tex /home/meisel/opus/tables/.");
392
393
394gSystem->Exec("cp -v Canv_SYSTEMATICS_ALLCOMBINED__*.png /home/meisel/opus/plotsHIGHRES_png/.");
395gSystem->Exec("cp -v Canv_SYSTEMATICS_ALLCOMBINED__*.eps /home/meisel/opus/plotsHIGHRES_eps/.");
396}
TPaveText * pt
Definition: Canv_SYSTEMATICS_ALLCOMBINED__RMSEnergy__vs__Energy__ELECTRON.C:160
Int_t PARTICLE_TYPE
Definition: Draw__SYSTEMATICS_ALLCOMBINED_Electron.C:2
const char filename[256]
Definition: RecDispNU.C:83
AcqOdyssey * o
Definition: hwinit.C:2
Int_t dummy
Definition: merge_Energy_SytematicSources_Electron.C:7
pFile
Definition: merge_Energy_SytematicSources_Electron.C:25
Int_t plateArray[13]
Definition: merge_Energy_SytematicSources_Electron.C:8
Float_t dummyf
Definition: merge_Energy_SytematicSources_Electron.C:9
const char * name
Definition: merge_Energy_SytematicSources_Electron.C:24
#define NULL
Definition: nidaqmx.h:84
p
Definition: testBGReduction_AllMethods.C:8