/
EMCAL_gamma.C
274 lines (267 loc) · 11.3 KB
/
EMCAL_gamma.C
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
168
169
170
171
172
173
174
175
176
177
178
179
180
181
182
183
184
185
186
187
188
189
190
191
192
193
194
195
196
197
198
199
200
201
202
203
204
205
206
207
208
209
210
211
212
213
214
215
216
217
218
219
220
221
222
223
224
225
226
227
228
229
230
231
232
233
234
235
236
237
238
239
240
241
242
243
244
245
246
247
248
249
250
251
252
253
254
255
256
257
258
259
260
261
262
263
264
265
266
267
268
269
270
271
272
273
274
{
//=========Macro generated from canvas: c1/c1
//========= (Fri Jul 22 13:12:33 2016) by ROOT version5.34/30
TCanvas *c1 = new TCanvas("c1", "c1",65,52,700,500);
gStyle->SetOptStat(0);
c1->Range(-1.75,-118.7813,10.75,1069.031);
c1->SetFillColor(0);
c1->SetBorderMode(0);
c1->SetBorderSize(2);
c1->SetFrameBorderMode(0);
c1->SetFrameBorderMode(0);
TH1D *gammaHist80EMC = new TH1D("gammaHist80EMC","gammaHistEMC #chi^{2} <= 80",100,-0.5,9.5);
gammaHist80EMC->SetBinContent(6,106);
gammaHist80EMC->SetBinContent(7,48);
gammaHist80EMC->SetBinContent(8,54);
gammaHist80EMC->SetBinContent(9,50);
gammaHist80EMC->SetBinContent(10,58);
gammaHist80EMC->SetBinContent(11,46);
gammaHist80EMC->SetBinContent(12,60);
gammaHist80EMC->SetBinContent(13,57);
gammaHist80EMC->SetBinContent(14,68);
gammaHist80EMC->SetBinContent(15,80);
gammaHist80EMC->SetBinContent(16,67);
gammaHist80EMC->SetBinContent(17,64);
gammaHist80EMC->SetBinContent(18,72);
gammaHist80EMC->SetBinContent(19,57);
gammaHist80EMC->SetBinContent(20,61);
gammaHist80EMC->SetBinContent(21,101);
gammaHist80EMC->SetBinContent(22,89);
gammaHist80EMC->SetBinContent(23,87);
gammaHist80EMC->SetBinContent(24,79);
gammaHist80EMC->SetBinContent(25,82);
gammaHist80EMC->SetBinContent(26,86);
gammaHist80EMC->SetBinContent(27,86);
gammaHist80EMC->SetBinContent(28,110);
gammaHist80EMC->SetBinContent(29,112);
gammaHist80EMC->SetBinContent(30,92);
gammaHist80EMC->SetBinContent(31,96);
gammaHist80EMC->SetBinContent(32,97);
gammaHist80EMC->SetBinContent(33,119);
gammaHist80EMC->SetBinContent(34,124);
gammaHist80EMC->SetBinContent(35,129);
gammaHist80EMC->SetBinContent(36,123);
gammaHist80EMC->SetBinContent(37,153);
gammaHist80EMC->SetBinContent(38,199);
gammaHist80EMC->SetBinContent(39,213);
gammaHist80EMC->SetBinContent(40,257);
gammaHist80EMC->SetBinContent(41,367);
gammaHist80EMC->SetBinContent(42,483);
gammaHist80EMC->SetBinContent(43,659);
gammaHist80EMC->SetBinContent(44,814);
gammaHist80EMC->SetBinContent(45,905);
gammaHist80EMC->SetBinContent(46,839);
gammaHist80EMC->SetBinContent(47,731);
gammaHist80EMC->SetBinContent(48,658);
gammaHist80EMC->SetBinContent(49,511);
gammaHist80EMC->SetBinContent(50,405);
gammaHist80EMC->SetBinContent(51,309);
gammaHist80EMC->SetBinContent(52,208);
gammaHist80EMC->SetBinContent(53,129);
gammaHist80EMC->SetBinContent(54,115);
gammaHist80EMC->SetBinContent(55,90);
gammaHist80EMC->SetBinContent(56,84);
gammaHist80EMC->SetBinContent(57,60);
gammaHist80EMC->SetBinContent(58,48);
gammaHist80EMC->SetBinContent(59,39);
gammaHist80EMC->SetBinContent(60,47);
gammaHist80EMC->SetBinContent(61,11);
gammaHist80EMC->SetBinContent(62,10);
gammaHist80EMC->SetBinContent(63,14);
gammaHist80EMC->SetBinContent(64,6);
gammaHist80EMC->SetBinContent(65,4);
gammaHist80EMC->SetBinContent(66,5);
gammaHist80EMC->SetBinContent(67,1);
gammaHist80EMC->SetBinContent(68,1);
gammaHist80EMC->SetBinContent(69,1);
gammaHist80EMC->SetBinContent(70,1);
gammaHist80EMC->SetBinContent(71,1);
gammaHist80EMC->SetBinContent(101,2);
gammaHist80EMC->SetEntries(10940);
Int_t ci; // for color index setting
TColor *color; // for color definition with alpha
ci = TColor::GetColor("#000099");
gammaHist80EMC->SetLineColor(ci);
gammaHist80EMC->GetXaxis()->SetLabelFont(42);
gammaHist80EMC->GetXaxis()->SetLabelSize(0.035);
gammaHist80EMC->GetXaxis()->SetTitleSize(0.035);
gammaHist80EMC->GetXaxis()->SetTitleFont(42);
gammaHist80EMC->GetYaxis()->SetLabelFont(42);
gammaHist80EMC->GetYaxis()->SetLabelSize(0.035);
gammaHist80EMC->GetYaxis()->SetTitleSize(0.035);
gammaHist80EMC->GetYaxis()->SetTitleFont(42);
gammaHist80EMC->GetZaxis()->SetLabelFont(42);
gammaHist80EMC->GetZaxis()->SetLabelSize(0.035);
gammaHist80EMC->GetZaxis()->SetTitleSize(0.035);
gammaHist80EMC->GetZaxis()->SetTitleFont(42);
gammaHist80EMC->Draw("");
TPaveText *pt = new TPaveText(0.2866954,0.9314407,0.7133046,0.995,"blNDC");
pt->SetName("title");
pt->SetBorderSize(0);
pt->SetFillColor(0);
pt->SetFillStyle(0);
pt->SetTextFont(42);
TText *text = pt->AddText("gammaHistEMC #chi^{2} <= 80");
pt->Draw();
TH1D *gammaHist200EMC = new TH1D("gammaHist200EMC","gammaHistEMC 80 < #chi^{2} <= 200",100,-0.5,9.5);
gammaHist200EMC->SetBinContent(6,1);
gammaHist200EMC->SetBinContent(8,1);
gammaHist200EMC->SetBinContent(10,1);
gammaHist200EMC->SetBinContent(11,3);
gammaHist200EMC->SetBinContent(12,2);
gammaHist200EMC->SetBinContent(13,1);
gammaHist200EMC->SetBinContent(14,2);
gammaHist200EMC->SetBinContent(16,3);
gammaHist200EMC->SetBinContent(17,4);
gammaHist200EMC->SetBinContent(18,3);
gammaHist200EMC->SetBinContent(19,2);
gammaHist200EMC->SetBinContent(20,1);
gammaHist200EMC->SetBinContent(21,2);
gammaHist200EMC->SetBinContent(22,5);
gammaHist200EMC->SetBinContent(23,4);
gammaHist200EMC->SetBinContent(24,1);
gammaHist200EMC->SetBinContent(25,5);
gammaHist200EMC->SetBinContent(26,2);
gammaHist200EMC->SetBinContent(27,3);
gammaHist200EMC->SetBinContent(28,2);
gammaHist200EMC->SetBinContent(29,4);
gammaHist200EMC->SetBinContent(30,3);
gammaHist200EMC->SetBinContent(31,4);
gammaHist200EMC->SetBinContent(32,3);
gammaHist200EMC->SetBinContent(33,2);
gammaHist200EMC->SetBinContent(34,1);
gammaHist200EMC->SetBinContent(35,1);
gammaHist200EMC->SetBinContent(36,4);
gammaHist200EMC->SetBinContent(37,2);
gammaHist200EMC->SetBinContent(38,4);
gammaHist200EMC->SetBinContent(39,7);
gammaHist200EMC->SetBinContent(40,6);
gammaHist200EMC->SetBinContent(41,8);
gammaHist200EMC->SetBinContent(42,13);
gammaHist200EMC->SetBinContent(43,21);
gammaHist200EMC->SetBinContent(44,17);
gammaHist200EMC->SetBinContent(45,23);
gammaHist200EMC->SetBinContent(46,17);
gammaHist200EMC->SetBinContent(47,23);
gammaHist200EMC->SetBinContent(48,18);
gammaHist200EMC->SetBinContent(49,13);
gammaHist200EMC->SetBinContent(50,9);
gammaHist200EMC->SetBinContent(51,7);
gammaHist200EMC->SetBinContent(52,3);
gammaHist200EMC->SetBinContent(53,5);
gammaHist200EMC->SetBinContent(54,3);
gammaHist200EMC->SetBinContent(55,7);
gammaHist200EMC->SetBinContent(56,2);
gammaHist200EMC->SetBinContent(57,1);
gammaHist200EMC->SetBinContent(58,2);
gammaHist200EMC->SetBinContent(59,4);
gammaHist200EMC->SetBinContent(60,1);
gammaHist200EMC->SetBinContent(61,1);
gammaHist200EMC->SetBinContent(63,2);
gammaHist200EMC->SetBinContent(64,2);
gammaHist200EMC->SetBinContent(66,2);
gammaHist200EMC->SetBinContent(101,2);
gammaHist200EMC->SetEntries(295);
ci = TColor::GetColor("#ff0000");
gammaHist200EMC->SetLineColor(ci);
gammaHist200EMC->GetXaxis()->SetLabelFont(42);
gammaHist200EMC->GetXaxis()->SetLabelSize(0.035);
gammaHist200EMC->GetXaxis()->SetTitleSize(0.035);
gammaHist200EMC->GetXaxis()->SetTitleFont(42);
gammaHist200EMC->GetYaxis()->SetLabelFont(42);
gammaHist200EMC->GetYaxis()->SetLabelSize(0.035);
gammaHist200EMC->GetYaxis()->SetTitleSize(0.035);
gammaHist200EMC->GetYaxis()->SetTitleFont(42);
gammaHist200EMC->GetZaxis()->SetLabelFont(42);
gammaHist200EMC->GetZaxis()->SetLabelSize(0.035);
gammaHist200EMC->GetZaxis()->SetTitleSize(0.035);
gammaHist200EMC->GetZaxis()->SetTitleFont(42);
gammaHist200EMC->Draw("same");
TH1D *gammaHistBeyondEMC = new TH1D("gammaHistBeyondEMC","gammaHistEMC #chi^{2} > 200",100,-0.5,9.5);
gammaHistBeyondEMC->SetBinContent(6,8);
gammaHistBeyondEMC->SetBinContent(7,3);
gammaHistBeyondEMC->SetBinContent(8,1);
gammaHistBeyondEMC->SetBinContent(9,3);
gammaHistBeyondEMC->SetBinContent(10,9);
gammaHistBeyondEMC->SetBinContent(12,6);
gammaHistBeyondEMC->SetBinContent(13,3);
gammaHistBeyondEMC->SetBinContent(14,4);
gammaHistBeyondEMC->SetBinContent(15,3);
gammaHistBeyondEMC->SetBinContent(16,6);
gammaHistBeyondEMC->SetBinContent(17,4);
gammaHistBeyondEMC->SetBinContent(18,1);
gammaHistBeyondEMC->SetBinContent(19,5);
gammaHistBeyondEMC->SetBinContent(20,12);
gammaHistBeyondEMC->SetBinContent(21,4);
gammaHistBeyondEMC->SetBinContent(22,6);
gammaHistBeyondEMC->SetBinContent(23,12);
gammaHistBeyondEMC->SetBinContent(24,8);
gammaHistBeyondEMC->SetBinContent(25,8);
gammaHistBeyondEMC->SetBinContent(26,10);
gammaHistBeyondEMC->SetBinContent(27,15);
gammaHistBeyondEMC->SetBinContent(28,8);
gammaHistBeyondEMC->SetBinContent(29,4);
gammaHistBeyondEMC->SetBinContent(30,13);
gammaHistBeyondEMC->SetBinContent(31,10);
gammaHistBeyondEMC->SetBinContent(32,9);
gammaHistBeyondEMC->SetBinContent(33,12);
gammaHistBeyondEMC->SetBinContent(34,15);
gammaHistBeyondEMC->SetBinContent(35,15);
gammaHistBeyondEMC->SetBinContent(36,12);
gammaHistBeyondEMC->SetBinContent(37,10);
gammaHistBeyondEMC->SetBinContent(38,24);
gammaHistBeyondEMC->SetBinContent(39,14);
gammaHistBeyondEMC->SetBinContent(40,26);
gammaHistBeyondEMC->SetBinContent(41,36);
gammaHistBeyondEMC->SetBinContent(42,43);
gammaHistBeyondEMC->SetBinContent(43,46);
gammaHistBeyondEMC->SetBinContent(44,42);
gammaHistBeyondEMC->SetBinContent(45,37);
gammaHistBeyondEMC->SetBinContent(46,32);
gammaHistBeyondEMC->SetBinContent(47,32);
gammaHistBeyondEMC->SetBinContent(48,18);
gammaHistBeyondEMC->SetBinContent(49,22);
gammaHistBeyondEMC->SetBinContent(50,14);
gammaHistBeyondEMC->SetBinContent(51,5);
gammaHistBeyondEMC->SetBinContent(52,9);
gammaHistBeyondEMC->SetBinContent(53,8);
gammaHistBeyondEMC->SetBinContent(54,2);
gammaHistBeyondEMC->SetBinContent(55,2);
gammaHistBeyondEMC->SetBinContent(56,2);
gammaHistBeyondEMC->SetBinContent(57,3);
gammaHistBeyondEMC->SetBinContent(58,3);
gammaHistBeyondEMC->SetBinContent(59,3);
gammaHistBeyondEMC->SetBinContent(60,5);
gammaHistBeyondEMC->SetBinContent(61,1);
gammaHistBeyondEMC->SetBinContent(63,1);
gammaHistBeyondEMC->SetBinContent(64,1);
gammaHistBeyondEMC->SetBinContent(65,1);
gammaHistBeyondEMC->SetBinContent(69,2);
gammaHistBeyondEMC->SetBinContent(70,3);
gammaHistBeyondEMC->SetBinContent(75,1);
gammaHistBeyondEMC->SetBinContent(76,1);
gammaHistBeyondEMC->SetBinContent(81,1);
gammaHistBeyondEMC->SetBinContent(85,1);
gammaHistBeyondEMC->SetBinContent(89,1);
gammaHistBeyondEMC->SetBinContent(93,1);
gammaHistBeyondEMC->SetBinContent(101,2);
gammaHistBeyondEMC->SetEntries(684);
ci = TColor::GetColor("#00ff00");
gammaHistBeyondEMC->SetLineColor(ci);
gammaHistBeyondEMC->GetXaxis()->SetLabelFont(42);
gammaHistBeyondEMC->GetXaxis()->SetLabelSize(0.035);
gammaHistBeyondEMC->GetXaxis()->SetTitleSize(0.035);
gammaHistBeyondEMC->GetXaxis()->SetTitleFont(42);
gammaHistBeyondEMC->GetYaxis()->SetLabelFont(42);
gammaHistBeyondEMC->GetYaxis()->SetLabelSize(0.035);
gammaHistBeyondEMC->GetYaxis()->SetTitleSize(0.035);
gammaHistBeyondEMC->GetYaxis()->SetTitleFont(42);
gammaHistBeyondEMC->GetZaxis()->SetLabelFont(42);
gammaHistBeyondEMC->GetZaxis()->SetLabelSize(0.035);
gammaHistBeyondEMC->GetZaxis()->SetTitleSize(0.035);
gammaHistBeyondEMC->GetZaxis()->SetTitleFont(42);
gammaHistBeyondEMC->Draw("same");
c1->Modified();
c1->cd();
c1->SetSelected(c1);
}