1 #include "sPhenixStyle.h"
2 #include "sPhenixStyle.C"
5 #include <TGraphErrors.h>
13 TH1::SetDefaultSumw2();
14 TH2::SetDefaultSumw2();
18 gStyle->SetOptStat(0);
21 TFile * fin = TFile::Open(
"../macro/hist_jets.root");
22 TH1F * _h_R04_z_g_15_20 = (TH1F*)fin->Get(
"R04_z_g_15_20");
23 _h_R04_z_g_15_20->Sumw2();
24 TH1F * _h_R04_z_g_20_25 = (TH1F*)fin->Get(
"R04_z_g_20_25");
25 _h_R04_z_g_20_25->Sumw2();
26 TH1F * _h_R04_z_g_25_30 = (TH1F*)fin->Get(
"R04_z_g_25_30");
27 _h_R04_z_g_25_30->Sumw2();
28 TH1F * _h_R04_z_g_30_40 = (TH1F*)fin->Get(
"R04_z_g_30_40");
29 _h_R04_z_g_30_40->Sumw2();
30 TH1F * _h_R04_theta_g_15_20 = (TH1F*)fin->Get(
"R04_theta_g_15_20");
31 _h_R04_theta_g_15_20->Sumw2();
32 TH1F * _h_R04_theta_g_20_25 = (TH1F*)fin->Get(
"R04_theta_g_20_25");
33 _h_R04_theta_g_20_25->Sumw2();
34 TH1F * _h_R04_theta_g_25_30 = (TH1F*)fin->Get(
"R04_theta_g_25_30");
35 _h_R04_theta_g_25_30->Sumw2();
36 TH1F * _h_R04_theta_g_30_40 = (TH1F*)fin->Get(
"R04_theta_g_30_40");
37 _h_R04_theta_g_30_40->Sumw2();
38 TFile * fun = TFile::Open(
"../macro/hist_jets.root");
39 TH1F * _h_R04_z_sj_15_20 = (TH1F*)fun->Get(
"R04_z_sj_15_20");
40 _h_R04_z_sj_15_20->Sumw2();
41 TH1F * _h_R04_z_sj_20_25 = (TH1F*)fun->Get(
"R04_z_sj_20_25");
42 _h_R04_z_sj_20_25->Sumw2();
43 TH1F * _h_R04_z_sj_25_30 = (TH1F*)fun->Get(
"R04_z_sj_25_30");
44 _h_R04_z_sj_25_30->Sumw2();
45 TH1F * _h_R04_z_sj_30_40 = (TH1F*)fun->Get(
"R04_z_sj_30_40");
46 _h_R04_z_sj_30_40->Sumw2();
47 TH1F * _h_R04_theta_sj_15_20 = (TH1F*)fun->Get(
"R04_theta_sj_15_20");
48 _h_R04_theta_sj_15_20->Sumw2();
49 TH1F * _h_R04_theta_sj_20_25 = (TH1F*)fun->Get(
"R04_theta_sj_20_25");
50 _h_R04_theta_sj_20_25->Sumw2();
51 TH1F * _h_R04_theta_sj_25_30 = (TH1F*)fun->Get(
"R04_theta_sj_25_30");
52 _h_R04_theta_sj_25_30->Sumw2();
53 TH1F * _h_R04_theta_sj_30_40 = (TH1F*)fun->Get(
"R04_theta_sj_30_40");
54 _h_R04_theta_sj_30_40->Sumw2();
55 TH2D * correlation_theta_15_20 = (TH2D*)fun->Get(
"correlation_theta_15_20");
56 TH2D * correlation_theta_20_25 = (TH2D*)fun->Get(
"correlation_theta_20_25");
57 TH2D * correlation_theta_25_30 = (TH2D*)fun->Get(
"correlation_theta_25_30");
58 TH2D * correlation_theta_30_40 = (TH2D*)fun->Get(
"correlation_theta_30_40");
60 TH2D * correlation_z_15_20 = (TH2D*)fun->Get(
"correlation_z_15_20");
61 TH2D * correlation_z_20_25 = (TH2D*)fun->Get(
"correlation_z_20_25");
62 TH2D * correlation_z_25_30 = (TH2D*)fun->Get(
"correlation_z_25_30");
63 TH2D * correlation_z_30_40 = (TH2D*)fun->Get(
"correlation_z_30_40");
65 TCanvas * csubjetR01 =
new TCanvas(
"csubjetR01",
" ", 1500, 500);
67 makeMultiPanelCanvasWithGap(csubjetR01, 4, 2, 0.2, 0.2, 0.2, 0.2, 0.08, 0.01);
78 _h_R04_z_sj_15_20->SetLineStyle(1);
79 _h_R04_z_sj_15_20->SetLineWidth(1);
80 _h_R04_z_sj_15_20->SetTitle(
" ");
81 _h_R04_z_sj_15_20->SetYTitle(
"(1/N_{jet}) dN/dz_{sj}");
82 _h_R04_z_sj_15_20->GetYaxis()->SetTitleOffset(0.8);
83 _h_R04_z_sj_15_20->GetYaxis()->SetTitleSize(0.07);
84 _h_R04_z_sj_15_20->GetYaxis()->SetLabelSize(0.04);
85 _h_R04_z_sj_15_20->GetYaxis()->CenterTitle(
true);
87 _h_R04_z_sj_15_20->SetXTitle(
"z_{sj}");
88 _h_R04_z_sj_15_20->GetXaxis()->SetTitleSize(0.07);
89 _h_R04_z_sj_15_20->GetXaxis()->CenterTitle(
true);
90 _h_R04_z_sj_15_20->GetXaxis()->SetLabelSize(0.04);
91 _h_R04_z_sj_15_20->SetMarkerSize(1.0);
92 _h_R04_z_sj_15_20->SetAxisRange(0, 5.5,
"Y");
93 _h_R04_z_sj_15_20->SetAxisRange(0, 0.5,
"X");
95 _h_R04_z_sj_15_20->Draw(
"hist C same");
97 int nbins15 = _h_R04_z_sj_15_20->GetNbinsX();
98 double* x15 =
new double[nbins15];
99 double* y15 =
new double[nbins15];
100 double* ex15 =
new double[nbins15];
101 double* ey15 =
new double[nbins15];
102 for (
int i = 1;
i <= nbins15; ++
i) {
103 x15[
i - 1] = _h_R04_z_sj_15_20 ->GetBinCenter(
i);
104 y15[
i - 1] = _h_R04_z_sj_15_20 ->GetBinContent(
i);
106 ey15[
i - 1] = _h_R04_z_sj_15_20 ->GetBinError(
i);
108 TGraphErrors* errorGraph15 =
new TGraphErrors(nbins15, x15, y15, ex15, ey15);
109 errorGraph15->SetMarkerStyle(29);
110 errorGraph15->SetMarkerColor(2);
111 errorGraph15->Draw(
"p same");
115 gStyle->SetErrorX(0);
120 _h_R04_z_sj_20_25->SetTitle(
" ");
121 _h_R04_z_sj_20_25->SetXTitle(
"z_{sj}");
122 _h_R04_z_sj_20_25->GetXaxis()->SetTitleSize(0.07);
123 _h_R04_z_sj_20_25->GetXaxis()->CenterTitle(
true);
124 _h_R04_z_sj_20_25->GetXaxis()->SetLabelSize(0.04);
125 _h_R04_z_sj_20_25->SetAxisRange(0, 5.5,
"Y");
126 _h_R04_z_sj_20_25->SetAxisRange(0, 0.5,
"X");
127 _h_R04_z_sj_20_25->SetMarkerSize(1.0);
129 _h_R04_z_sj_20_25->Draw(
"hist C same");
132 int nbins20 = _h_R04_z_sj_20_25->GetNbinsX();
133 double* x20 =
new double[nbins20];
134 double* y20 =
new double[nbins20];
135 double* ex20 =
new double[nbins20];
136 double* ey20 =
new double[nbins20];
137 for (
int i = 1;
i <= nbins20; ++
i) {
138 x20[
i - 1] = _h_R04_z_sj_20_25->GetBinCenter(
i);
139 y20[
i - 1] = _h_R04_z_sj_20_25->GetBinContent(
i);
141 ey20[
i - 1] = _h_R04_z_sj_20_25->GetBinError(
i);
143 TGraphErrors* errorGraph20 =
new TGraphErrors(nbins20, x20, y20, ex20, ey20);
144 errorGraph20->SetMarkerStyle(29);
145 errorGraph20->SetMarkerColor(2);
146 errorGraph20->Draw(
"p same");
149 gStyle->SetErrorX(0);
154 _h_R04_z_sj_25_30->SetMarkerStyle(29);
155 _h_R04_z_sj_25_30->SetMarkerColor(2);
156 _h_R04_z_sj_25_30->SetTitle(
" ");
157 _h_R04_z_sj_25_30->SetXTitle(
"z_{sj}");
158 _h_R04_z_sj_25_30->GetXaxis()->SetTitleSize(0.07);
159 _h_R04_z_sj_25_30->GetXaxis()->CenterTitle(
true);
160 _h_R04_z_sj_25_30->GetXaxis()->SetLabelSize(0.04);
161 _h_R04_z_sj_25_30->SetAxisRange(0, 5.5,
"Y");
162 _h_R04_z_sj_25_30->SetAxisRange(0, 0.5,
"X");
163 _h_R04_z_sj_25_30->SetMarkerSize(1.0);
165 _h_R04_z_sj_25_30->Draw(
"hist C same");
167 int nbins25 = _h_R04_z_sj_25_30->GetNbinsX();
168 double* x25 =
new double[nbins25];
169 double* y25 =
new double[nbins25];
170 double* ex25 =
new double[nbins25];
171 double* ey25 =
new double[nbins25];
172 for (
int i = 1;
i <= nbins25; ++
i) {
173 x25[
i - 1] = _h_R04_z_sj_25_30->GetBinCenter(
i);
174 y25[
i - 1] = _h_R04_z_sj_25_30->GetBinContent(
i);
176 ey25[
i - 1] = _h_R04_z_sj_25_30->GetBinError(
i);
178 TGraphErrors* errorGraph25 =
new TGraphErrors(nbins25, x25, y25, ex25, ey25);
179 errorGraph25->SetMarkerStyle(29);
180 errorGraph25->SetMarkerColor(2);
181 errorGraph25->Draw(
"p same");
184 gStyle->SetErrorX(0);
187 _h_R04_z_sj_30_40->SetMarkerStyle(29);
188 _h_R04_z_sj_30_40->SetMarkerColor(2);
189 _h_R04_z_sj_30_40->SetTitle(
" ");
190 _h_R04_z_sj_30_40->SetXTitle(
"z_{sj}");
191 _h_R04_z_sj_30_40->GetXaxis()->SetTitleSize(0.07);
192 _h_R04_z_sj_30_40->GetXaxis()->CenterTitle(
true);
193 _h_R04_z_sj_30_40->GetXaxis()->SetLabelSize(0.04);
194 _h_R04_z_sj_30_40->GetXaxis()->SetTitleSize(0.07);
195 _h_R04_z_sj_30_40->SetAxisRange(0, 5.5,
"Y");
196 _h_R04_z_sj_30_40->SetAxisRange(0, 0.5,
"X");
197 _h_R04_z_sj_30_40->SetMarkerSize(1.0);
199 _h_R04_z_sj_30_40->Draw(
"hist C same");
201 drawText(
"p+p #sqrt{s} = 200 GeV", 0.3, 0.8, 14);
202 drawText(
"anti-k_{T} R=0.4 jets, |#eta_{jet}| < 0.6", 0.3, 0.7, 14);
203 drawText(
"R_{sj}=0.1", 0.3, 0.6, 14);
206 int nbins30 = _h_R04_z_sj_30_40->GetNbinsX();
207 double* x30 =
new double[nbins30];
208 double* y30 =
new double[nbins30];
209 double* ex30 =
new double[nbins30];
210 double* ey30 =
new double[nbins30];
211 for (
int i = 1;
i <= nbins30; ++
i) {
212 x30[
i - 1] = _h_R04_z_sj_30_40->GetBinCenter(
i);
213 y30[
i - 1] = _h_R04_z_sj_30_40->GetBinContent(
i);
215 ey30[
i - 1] = _h_R04_z_sj_30_40->GetBinError(
i);
217 TGraphErrors* errorGraph30 =
new TGraphErrors(nbins15, x30, y30, ex30, ey30);
218 errorGraph30->SetMarkerStyle(29);
219 errorGraph30->SetMarkerColor(2);
220 errorGraph30->Draw(
"p same");
223 gStyle->SetErrorX(0);
226 _h_R04_theta_sj_15_20->SetMarkerStyle(29);
227 _h_R04_theta_sj_15_20->SetMarkerColor(2);
228 _h_R04_theta_sj_15_20->SetMarkerSize(1.0);
229 _h_R04_theta_sj_15_20->SetTitle(
" ");
230 _h_R04_theta_sj_15_20->SetYTitle(
"(1/N_{jet}) dN/d#theta_{sj}");
231 _h_R04_theta_sj_15_20->SetAxisRange(0, 9.0,
"Y");
232 _h_R04_theta_sj_15_20->SetAxisRange(0, 0.5,
"X");
233 _h_R04_theta_sj_15_20->GetYaxis()->SetTitleOffset(0.8);
234 _h_R04_theta_sj_15_20->GetYaxis()->SetTitleSize(0.07);
235 _h_R04_theta_sj_15_20->GetYaxis()->SetLabelSize(0.04);
236 _h_R04_theta_sj_15_20->GetYaxis()->CenterTitle(
true);
238 _h_R04_theta_sj_15_20->SetXTitle(
"#theta_{sj}");
239 _h_R04_theta_sj_15_20->GetXaxis()->SetTitleSize(0.07);
240 _h_R04_theta_sj_15_20->GetXaxis()->CenterTitle(
true);
241 _h_R04_theta_sj_15_20->GetXaxis()->SetLabelSize(0.04);
243 _h_R04_theta_sj_15_20->Draw(
"hist C same");
245 int nbins15t = _h_R04_theta_sj_15_20->GetNbinsX();
246 double* x15t =
new double[nbins15t];
247 double* y15t =
new double[nbins15t];
248 double* ex15t =
new double[nbins15t];
249 double* ey15t =
new double[nbins15t];
250 for (
int i = 1;
i <= nbins15t; ++
i) {
251 x15t[
i - 1] = _h_R04_theta_sj_15_20->GetBinCenter(
i);
252 y15t[
i - 1] = _h_R04_theta_sj_15_20->GetBinContent(
i);
254 ey15t[
i - 1] = _h_R04_theta_sj_15_20->GetBinError(
i);
256 TGraphErrors* errorGraph15t =
new TGraphErrors(nbins15t, x15t, y15t, ex15t, ey15t);
257 errorGraph15t->SetMarkerStyle(29);
258 errorGraph15t->SetMarkerColor(2);
259 errorGraph15t->Draw(
"p same");
261 drawText(
"15 < p_{T} < 20 [GeV/c]", 0.6, 0.8, 12);
264 gStyle->SetErrorX(0);
269 _h_R04_theta_sj_20_25->SetMarkerStyle(29);
270 _h_R04_theta_sj_20_25->SetMarkerColor(2);
271 _h_R04_theta_sj_20_25->SetMarkerSize(1.0);
272 _h_R04_theta_sj_20_25->SetTitle(
" ");
273 _h_R04_theta_sj_20_25->SetXTitle(
"#theta_{sj}");
274 _h_R04_theta_sj_20_25->GetXaxis()->SetTitleSize(0.07);
275 _h_R04_theta_sj_20_25->GetXaxis()->CenterTitle(
true);
276 _h_R04_theta_sj_20_25->GetXaxis()->SetLabelSize(0.04);
277 _h_R04_theta_sj_20_25->SetAxisRange(0, 9.0,
"Y");
278 _h_R04_theta_sj_20_25->SetAxisRange(0, 0.5,
"X");
279 _h_R04_theta_sj_20_25->Draw(
"hist C same");
282 int nbins20t = _h_R04_theta_sj_20_25->GetNbinsX();
283 double* x20t =
new double[nbins20t];
284 double* y20t =
new double[nbins20t];
285 double* ex20t =
new double[nbins20t];
286 double* ey20t =
new double[nbins20t];
287 for (
int i = 1;
i <= nbins15t; ++
i) {
288 x20t[
i - 1] = _h_R04_theta_sj_20_25->GetBinCenter(
i);
289 y20t[
i - 1] = _h_R04_theta_sj_20_25->GetBinContent(
i);
291 ey20t[
i - 1] = _h_R04_theta_sj_20_25->GetBinError(
i);
293 TGraphErrors* errorGraph20t =
new TGraphErrors(nbins20t, x20t, y20t, ex20t, ey20t);
294 errorGraph20t->SetMarkerStyle(29);
295 errorGraph20t->SetMarkerColor(2);
296 errorGraph20t->Draw(
"p same");
297 drawText(
"20 < p_{T} < 25 [GeV/c]", 0.5, 0.8, 12);
300 gStyle->SetErrorX(0);
305 _h_R04_theta_sj_25_30->SetMarkerStyle(29);
306 _h_R04_theta_sj_25_30->SetMarkerColor(2);
307 _h_R04_theta_sj_25_30->SetMarkerSize(1.0);
308 _h_R04_theta_sj_25_30->SetTitle(
" ");
309 _h_R04_theta_sj_25_30->SetXTitle(
"#theta_{sj}");
310 _h_R04_theta_sj_25_30->GetXaxis()->SetTitleSize(0.07);
311 _h_R04_theta_sj_25_30->GetXaxis()->CenterTitle(
true);
312 _h_R04_theta_sj_25_30->GetXaxis()->SetLabelSize(0.04);
313 _h_R04_theta_sj_25_30->SetAxisRange(0, 9.0,
"Y");
314 _h_R04_theta_sj_25_30->SetAxisRange(0, 0.5,
"X");
317 _h_R04_theta_sj_25_30->Draw(
"hist C same");
319 int nbins25t = _h_R04_theta_sj_25_30->GetNbinsX();
320 double* x25t =
new double[nbins25t];
321 double* y25t =
new double[nbins25t];
322 double* ex25t =
new double[nbins25t];
323 double* ey25t =
new double[nbins25t];
324 for (
int i = 1;
i <= nbins25t; ++
i) {
325 x25t[
i - 1] = _h_R04_theta_sj_25_30->GetBinCenter(
i);
326 y25t[
i - 1] = _h_R04_theta_sj_25_30->GetBinContent(
i);
328 ey25t[
i - 1] = _h_R04_theta_sj_25_30->GetBinError(
i);
330 TGraphErrors* errorGraph25t =
new TGraphErrors(nbins25t, x25t, y25t, ex25t, ey25t);
331 errorGraph25t->SetMarkerStyle(29);
332 errorGraph25t->SetMarkerColor(2);
333 errorGraph25t->Draw(
"p same");
334 drawText(
"25 < p_{T} < 30 [GeV/c]", 0.5, 0.8, 12);
337 gStyle->SetErrorX(0);
340 _h_R04_theta_sj_30_40->SetMarkerStyle(29);
341 _h_R04_theta_sj_30_40->SetMarkerColor(2);
342 _h_R04_theta_sj_30_40->SetMarkerSize(1.0);
343 _h_R04_theta_sj_30_40->SetTitle(
" ");
344 _h_R04_theta_sj_30_40->SetXTitle(
"#theta_{sj}");
345 _h_R04_theta_sj_30_40->GetXaxis()->SetTitleSize(0.07);
346 _h_R04_theta_sj_30_40->GetXaxis()->CenterTitle(
true);
347 _h_R04_theta_sj_30_40->GetXaxis()->SetLabelSize(0.04);
348 _h_R04_theta_sj_30_40->GetXaxis()->SetTitleSize(0.07);
349 _h_R04_theta_sj_30_40->SetAxisRange(0, 9.0,
"Y");
350 _h_R04_theta_sj_30_40->SetAxisRange(0, 0.5,
"X");
352 _h_R04_theta_sj_30_40->Draw(
"hist C same");
354 int nbins30t = _h_R04_theta_sj_30_40->GetNbinsX();
355 double* x30t =
new double[nbins30t];
356 double* y30t =
new double[nbins30t];
357 double* ex30t =
new double[nbins30t];
358 double* ey30t =
new double[nbins30t];
359 for (
int i = 1;
i <= nbins30t; ++
i) {
360 x30t[
i - 1] = _h_R04_theta_sj_30_40->GetBinCenter(
i);
361 y30t[
i - 1] = _h_R04_theta_sj_30_40->GetBinContent(
i);
363 ey30t[
i - 1] = _h_R04_theta_sj_30_40->GetBinError(
i);
365 TGraphErrors* errorGraph30t =
new TGraphErrors(nbins30t, x30t, y30t, ex30t, ey30t);
366 errorGraph30t->SetMarkerStyle(29);
367 errorGraph30t->SetMarkerColor(2);
368 errorGraph30t->Draw(
"p same");
370 drawText(
"30 < p_{T} < 40 [GeV/c]", 0.5, 0.8, 12);
372 csubjetR01->SaveAs(
"RHIC_subjetR04_AKTB0z1.pdf",
"RECREATE");
375 TCanvas * csubjetR04 =
new TCanvas(
"csubjetR04",
" ", 1500, 500);
377 makeMultiPanelCanvasWithGap(csubjetR04, 4, 2, 0.2, 0.2, 0.2, 0.2, 0.08, 0.01);
380 gStyle->SetErrorX(0);
383 _h_R04_z_g_15_20 ->SetLineStyle(1);
384 _h_R04_z_g_15_20->SetLineWidth(1);
385 _h_R04_z_g_15_20->SetTitle(
" ");
386 _h_R04_z_g_15_20->SetYTitle(
"(1/N_{jet}) dN/dz_{g}");
387 _h_R04_z_g_15_20->GetYaxis()->SetTitleOffset(0.8);
388 _h_R04_z_g_15_20->GetYaxis()->SetTitleSize(0.07);
389 _h_R04_z_g_15_20->GetYaxis()->SetLabelSize(0.04);
390 _h_R04_z_g_15_20->GetYaxis()->CenterTitle(
true);
392 _h_R04_z_g_15_20->SetXTitle(
"z_{g}");
393 _h_R04_z_g_15_20->GetXaxis()->SetTitleSize(0.07);
394 _h_R04_z_g_15_20->GetXaxis()->CenterTitle(
true);
395 _h_R04_z_g_15_20->GetXaxis()->SetLabelSize(0.04);
396 _h_R04_z_g_15_20->SetMarkerSize(1.0);
397 _h_R04_z_g_15_20->SetAxisRange(0, 6.5,
"Y");
398 _h_R04_z_g_15_20->SetAxisRange(0, 0.5,
"X");
400 _h_R04_z_g_15_20->Draw(
"hist C same");
402 int nbins15zg = _h_R04_z_g_15_20->GetNbinsX();
403 double* x15zg =
new double[nbins15zg];
404 double* y15zg =
new double[nbins15zg];
405 double* ex15zg =
new double[nbins15zg];
406 double* ey15zg =
new double[nbins15zg];
407 for (
int i = 1;
i <= nbins15zg; ++
i) {
408 x15zg[
i - 1] = _h_R04_z_g_15_20->GetBinCenter(
i);
409 y15zg[
i - 1] = _h_R04_z_g_15_20->GetBinContent(
i);
411 ey15zg[
i - 1] = _h_R04_z_g_15_20->GetBinError(
i);
413 TGraphErrors* errorGraph15zg =
new TGraphErrors(nbins15zg, x15zg, y15zg, ex15zg, ey15zg);
414 errorGraph15zg->SetMarkerStyle(29);
415 errorGraph15zg->SetMarkerColor(2);
416 errorGraph15zg->Draw(
"p same");
420 gStyle->SetErrorX(0);
425 _h_R04_z_g_20_25->SetTitle(
" ");
426 _h_R04_z_g_20_25->SetXTitle(
"z_{g}");
427 _h_R04_z_g_20_25->GetXaxis()->SetTitleSize(0.07);
428 _h_R04_z_g_20_25->GetXaxis()->CenterTitle(
true);
429 _h_R04_z_g_20_25->GetXaxis()->SetLabelSize(0.04);
430 _h_R04_z_g_20_25->SetAxisRange(0, 6.5,
"Y");
431 _h_R04_z_g_20_25->SetAxisRange(0, 0.5,
"X");
432 _h_R04_z_g_20_25->SetMarkerSize(1.0);
434 _h_R04_z_g_20_25->Draw(
"hist C same");
437 int nbins20zg = _h_R04_z_g_20_25->GetNbinsX();
438 double* x20zg =
new double[nbins20zg];
439 double* y20zg =
new double[nbins20zg];
440 double* ex20zg =
new double[nbins20zg];
441 double* ey20zg =
new double[nbins20zg];
442 for (
int i = 1;
i <= nbins20; ++
i) {
443 x20zg[
i - 1] = _h_R04_z_g_20_25->GetBinCenter(
i);
444 y20zg[
i - 1] = _h_R04_z_g_20_25->GetBinContent(
i);
446 ey20zg[
i - 1] = _h_R04_z_g_20_25->GetBinError(
i);
448 TGraphErrors* errorGraph20zg =
new TGraphErrors(nbins20, x20zg, y20zg, ex20zg, ey20zg);
449 errorGraph20zg->SetMarkerStyle(29);
450 errorGraph20zg->SetMarkerColor(2);
451 errorGraph20zg->Draw(
"p same");
454 gStyle->SetErrorX(0);
459 _h_R04_z_g_25_30->SetMarkerStyle(29);
460 _h_R04_z_g_25_30->SetMarkerColor(2);
461 _h_R04_z_g_25_30->SetTitle(
" ");
462 _h_R04_z_g_25_30->SetXTitle(
"z_{g}");
463 _h_R04_z_g_25_30->GetXaxis()->SetTitleSize(0.07);
464 _h_R04_z_g_25_30->GetXaxis()->CenterTitle(
true);
465 _h_R04_z_g_25_30->GetXaxis()->SetLabelSize(0.04);
466 _h_R04_z_g_25_30->SetAxisRange(0, 6.5,
"Y");
467 _h_R04_z_g_25_30->SetAxisRange(0, 0.5,
"X");
468 _h_R04_z_g_25_30->SetMarkerSize(1.0);
470 _h_R04_z_g_25_30->Draw(
"hist C same");
472 int nbins25zg = _h_R04_z_g_25_30->GetNbinsX();
473 double* x25zg =
new double[nbins25zg];
474 double* y25zg =
new double[nbins25zg];
475 double* ex25zg =
new double[nbins25zg];
476 double* ey25zg =
new double[nbins25zg];
477 for (
int i = 1;
i <= nbins25zg; ++
i) {
478 x25zg[
i - 1] = _h_R04_z_g_25_30->GetBinCenter(
i);
479 y25zg[
i - 1] = _h_R04_z_g_25_30->GetBinContent(
i);
481 ey25zg[
i - 1] = _h_R04_z_g_25_30->GetBinError(
i);
483 TGraphErrors* errorGraph25zg =
new TGraphErrors(nbins25zg, x25zg, y25zg, ex25zg, ey25zg);
484 errorGraph25zg->SetMarkerStyle(29);
485 errorGraph25zg->SetMarkerColor(2);
486 errorGraph25zg->Draw(
"p same");
489 gStyle->SetErrorX(0);
492 _h_R04_z_g_30_40->SetMarkerStyle(29);
493 _h_R04_z_g_30_40->SetMarkerColor(2);
494 _h_R04_z_g_30_40->SetTitle(
" ");
495 _h_R04_z_g_30_40->SetXTitle(
"z_{g}");
496 _h_R04_z_g_30_40->GetXaxis()->SetTitleSize(0.07);
497 _h_R04_z_g_30_40->GetXaxis()->CenterTitle(
true);
498 _h_R04_z_g_30_40->GetXaxis()->SetLabelSize(0.04);
499 _h_R04_z_g_30_40->GetXaxis()->SetTitleSize(0.07);
500 _h_R04_z_g_30_40->SetAxisRange(0, 6.5,
"Y");
501 _h_R04_z_g_30_40->SetAxisRange(0, 0.5,
"X");
502 _h_R04_z_g_30_40->SetMarkerSize(1.0);
504 _h_R04_z_g_30_40->Draw(
"hist C same");
506 drawText(
"p+p #sqrt{s} = 200 GeV", 0.4, 0.85, 12);
507 drawText(
"anti-k_{T} R=0.4 jets, |#eta_{jet}| < 0.6", 0.4, 0.75, 12);
508 drawText(
"#beta=0.0, z_{cut}=0.1", 0.4, 0.65, 12);
511 int nbins30zg = _h_R04_z_g_30_40->GetNbinsX();
512 double* x30zg =
new double[nbins30zg];
513 double* y30zg =
new double[nbins30zg];
514 double* ex30zg =
new double[nbins30zg];
515 double* ey30zg =
new double[nbins30zg];
516 for (
int i = 1;
i <= nbins20zg; ++
i) {
517 x30zg[
i - 1] = _h_R04_z_g_30_40->GetBinCenter(
i);
518 y30zg[
i - 1] = _h_R04_z_g_30_40->GetBinContent(
i);
520 ey30zg[
i - 1] = _h_R04_z_g_30_40->GetBinError(
i);
522 TGraphErrors* errorGraph30zg =
new TGraphErrors(nbins15zg, x30zg, y30zg, ex30zg, ey30zg);
523 errorGraph30zg->SetMarkerStyle(29);
524 errorGraph30zg->SetMarkerColor(2);
525 errorGraph30zg->Draw(
"p same");
528 gStyle->SetErrorX(0);
531 _h_R04_theta_g_15_20->SetMarkerStyle(29);
532 _h_R04_theta_g_15_20->SetMarkerColor(2);
533 _h_R04_theta_g_15_20->SetMarkerSize(1.0);
534 _h_R04_theta_g_15_20->SetTitle(
" ");
535 _h_R04_theta_g_15_20->SetYTitle(
"(1/N_{jet}) dN/dR_{g}");
536 _h_R04_theta_g_15_20->SetAxisRange(0, 9.0,
"Y");
537 _h_R04_theta_g_15_20->SetAxisRange(0, 0.5,
"X");
538 _h_R04_theta_g_15_20->GetYaxis()->SetTitleOffset(0.8);
539 _h_R04_theta_g_15_20->GetYaxis()->SetTitleSize(0.07);
540 _h_R04_theta_g_15_20->GetYaxis()->SetLabelSize(0.04);
541 _h_R04_theta_g_15_20->GetYaxis()->CenterTitle(
true);
543 _h_R04_theta_g_15_20->SetXTitle(
"R_{g}");
544 _h_R04_theta_g_15_20->GetXaxis()->SetTitleSize(0.07);
545 _h_R04_theta_g_15_20->GetXaxis()->CenterTitle(
true);
546 _h_R04_theta_g_15_20->GetXaxis()->SetLabelSize(0.04);
548 _h_R04_theta_g_15_20->Draw(
"hist C same");
550 int nbins15tgt = _h_R04_theta_g_15_20->GetNbinsX();
551 double* x15gt =
new double[nbins15tgt];
552 double* y15gt =
new double[nbins15tgt];
553 double* ex15gt =
new double[nbins15tgt];
554 double* ey15gt =
new double[nbins15tgt];
555 for (
int i = 1;
i <= nbins15tgt; ++
i) {
556 x15gt[
i - 1] = _h_R04_theta_g_15_20->GetBinCenter(
i);
557 y15gt[
i - 1] = _h_R04_theta_g_15_20->GetBinContent(
i);
559 ey15gt[
i - 1] = _h_R04_theta_g_15_20->GetBinError(
i);
561 TGraphErrors* errorGraph15gt =
new TGraphErrors(nbins15tgt, x15gt, y15gt, ex15gt, ey15gt);
562 errorGraph15gt->SetMarkerStyle(29);
563 errorGraph15gt->SetMarkerColor(2);
564 errorGraph15gt->Draw(
"p same");
566 drawText(
"15 < p_{T} < 20 GeV/c", 0.7, 0.8, 12);
569 gStyle->SetErrorX(0);
574 _h_R04_theta_g_20_25->SetMarkerStyle(29);
575 _h_R04_theta_g_20_25->SetMarkerColor(2);
576 _h_R04_theta_g_20_25->SetMarkerSize(1.0);
577 _h_R04_theta_g_20_25->SetTitle(
" ");
578 _h_R04_theta_g_20_25->SetXTitle(
"R_{g}");
579 _h_R04_theta_g_20_25->GetXaxis()->SetTitleSize(0.07);
580 _h_R04_theta_g_20_25->GetXaxis()->CenterTitle(
true);
581 _h_R04_theta_g_20_25->GetXaxis()->SetLabelSize(0.04);
582 _h_R04_theta_g_20_25->SetAxisRange(0, 9.0,
"Y");
583 _h_R04_theta_g_20_25->SetAxisRange(0, 0.5,
"X");
584 _h_R04_theta_g_20_25->Draw(
"hist C same");
587 int nbins20tgt = _h_R04_theta_g_20_25->GetNbinsX();
588 double* x20gt =
new double[nbins20tgt];
589 double* y20gt =
new double[nbins20tgt];
590 double* ex20gt =
new double[nbins20tgt];
591 double* ey20gt =
new double[nbins20tgt];
592 for (
int i = 1;
i <= nbins15tgt; ++
i) {
593 x20gt[
i - 1] = _h_R04_theta_g_20_25->GetBinCenter(
i);
594 y20gt[
i - 1] = _h_R04_theta_g_20_25->GetBinContent(
i);
596 ey20gt[
i - 1] = _h_R04_theta_g_20_25->GetBinError(
i);
598 TGraphErrors* errorGraph20gt =
new TGraphErrors(nbins20tgt, x20gt, y20gt, ex20gt, ey20gt);
599 errorGraph20gt->SetMarkerStyle(29);
600 errorGraph20gt->SetMarkerColor(2);
601 errorGraph20gt->Draw(
"p same");
602 drawText(
"20 < p_{T} < 25 GeV/c", 0.6, 0.8, 12);
605 gStyle->SetErrorX(0);
610 _h_R04_theta_g_25_30->SetMarkerStyle(29);
611 _h_R04_theta_g_25_30->SetMarkerColor(2);
612 _h_R04_theta_g_25_30->SetMarkerSize(1.0);
613 _h_R04_theta_g_25_30->SetTitle(
" ");
614 _h_R04_theta_g_25_30->SetXTitle(
"R_{g}");
615 _h_R04_theta_g_25_30->GetXaxis()->SetTitleSize(0.07);
616 _h_R04_theta_g_25_30->GetXaxis()->CenterTitle(
true);
617 _h_R04_theta_g_25_30->GetXaxis()->SetLabelSize(0.04);
618 _h_R04_theta_g_25_30->SetAxisRange(0, 9.0,
"Y");
619 _h_R04_theta_g_25_30->SetAxisRange(0, 0.5,
"X");
622 _h_R04_theta_g_25_30->Draw(
"hist C same");
624 int nbins25tgt = _h_R04_theta_g_25_30->GetNbinsX();
625 double* x25gt =
new double[nbins25tgt];
626 double* y25gt =
new double[nbins25tgt];
627 double* ex25gt =
new double[nbins25tgt];
628 double* ey25gt =
new double[nbins25tgt];
629 for (
int i = 1;
i <= nbins25tgt; ++
i) {
630 x25gt[
i - 1] = _h_R04_theta_g_25_30->GetBinCenter(
i);
631 y25gt[
i - 1] = _h_R04_theta_g_25_30->GetBinContent(
i);
634 ey25gt[
i - 1] = _h_R04_theta_g_25_30->GetBinError(
i);
636 TGraphErrors* errorGraph25gt =
new TGraphErrors(nbins25tgt, x25gt, y25gt, ex25gt, ey25gt);
637 errorGraph25gt->SetMarkerStyle(29);
638 errorGraph25gt->SetMarkerColor(2);
639 errorGraph25gt->Draw(
"p same");
640 drawText(
"25 < p_{T} < 30 GeV/c", 0.6, 0.8, 12);
643 gStyle->SetErrorX(0);
646 _h_R04_theta_g_30_40->SetMarkerStyle(29);
647 _h_R04_theta_g_30_40->SetMarkerColor(2);
648 _h_R04_theta_g_30_40->SetMarkerSize(1.0);
649 _h_R04_theta_g_30_40->SetTitle(
" ");
650 _h_R04_theta_g_30_40->SetXTitle(
"R_{g}");
651 _h_R04_theta_g_30_40->GetXaxis()->SetTitleSize(0.07);
652 _h_R04_theta_g_30_40->GetXaxis()->CenterTitle(
true);
653 _h_R04_theta_g_30_40->GetXaxis()->SetLabelSize(0.04);
654 _h_R04_theta_g_30_40->GetXaxis()->SetTitleSize(0.07);
655 _h_R04_theta_g_30_40->SetAxisRange(0, 9.0,
"Y");
656 _h_R04_theta_g_30_40->SetAxisRange(0, 0.5,
"X");
658 _h_R04_theta_g_30_40->Draw(
"hist C same");
660 int nbins30tgt = _h_R04_theta_g_30_40->GetNbinsX();
661 double* x30gt =
new double[nbins30tgt];
662 double* y30gt =
new double[nbins30tgt];
663 double* ex30gt =
new double[nbins30tgt];
664 double* ey30gt =
new double[nbins30tgt];
665 for (
int i = 1;
i <= nbins20tgt; ++
i) {
666 x30gt[
i - 1] = _h_R04_theta_g_30_40->GetBinCenter(
i);
667 y30gt[
i - 1] = _h_R04_theta_g_30_40->GetBinContent(
i);
670 ey30gt[
i - 1] = _h_R04_theta_g_30_40->GetBinError(
i);
672 TGraphErrors* errorGraph30gt =
new TGraphErrors(nbins20tgt, x30gt, y30gt, ex30gt, ey30gt);
673 errorGraph30gt->SetMarkerStyle(29);
674 errorGraph30gt->SetMarkerColor(2);
675 errorGraph30gt->Draw(
"p same");
677 drawText(
"30 < p_{T} < 40 GeV/c", 0.55, 0.8, 12);
679 csubjetR04->SaveAs(
"RHIC_subjetR04_CAB0z1.pdf",
"RECREATE");
682 TCanvas * canvas =
new TCanvas(
"canvas",
" ", 1500, 500);
684 makeMultiPanelCanvasWithGap(canvas, 4, 2, 0.2, 0.2, 0.2, 0.2, 0.09, 0.01);
687 gStyle->SetErrorX(0);
691 correlation_z_15_20->SetMarkerStyle(29);
692 correlation_z_15_20->SetMarkerColor(2);
693 correlation_z_15_20->SetMarkerSize(1.0);
694 correlation_z_15_20->GetXaxis()->SetTitleSize(0.07);
695 correlation_z_15_20->GetYaxis()->SetTitleSize(0.07);
696 correlation_z_15_20->GetXaxis()->CenterTitle(
true);
697 correlation_z_15_20->GetXaxis()->SetLabelSize(0.04);
698 correlation_z_15_20->GetXaxis()->SetTitleSize(0.07);
699 correlation_z_15_20->SetAxisRange(0, 0.6,
"Y");
700 correlation_z_15_20->SetAxisRange(0, 0.5,
"X");
701 correlation_z_15_20->SetXTitle(
"z_{g}");
702 correlation_z_15_20->SetYTitle(
"z_{sj}");
703 correlation_z_15_20->GetYaxis()->SetLabelSize(0.04);
704 correlation_z_15_20->SetTitle(
" ");
705 correlation_z_15_20->Draw(
"colz");
706 drawText(
"15 < p_{T} < 20 GeV/c", 0.4, 0.75, 12);
710 gStyle->SetErrorX(0);
714 correlation_z_20_25->SetMarkerStyle(29);
715 correlation_z_20_25->SetMarkerColor(2);
716 correlation_z_20_25->SetMarkerSize(1.0);
717 correlation_z_20_25->GetXaxis()->SetTitleSize(0.07);
718 correlation_z_20_25->GetXaxis()->CenterTitle(
true);
719 correlation_z_20_25->GetXaxis()->SetLabelSize(0.04);
720 correlation_z_20_25->GetXaxis()->SetTitleSize(0.07);
721 correlation_z_20_25->SetAxisRange(0, 0.6,
"Y");
722 correlation_z_20_25->SetAxisRange(0, 0.5,
"X");
723 correlation_z_20_25->SetXTitle(
"z_{g}");
725 correlation_z_20_25->SetTitle(
" ");
726 correlation_z_20_25->Draw(
"colz");
727 drawText(
"20 < p_{T} < 25 GeV/c", 0.3, 0.75, 12);
731 gStyle->SetErrorX(0);
735 correlation_z_25_30->SetMarkerStyle(29);
736 correlation_z_25_30->SetMarkerColor(2);
737 correlation_z_25_30->SetMarkerSize(1.0);
738 correlation_z_25_30->GetXaxis()->SetTitleSize(0.07);
739 correlation_z_25_30->GetXaxis()->CenterTitle(
true);
740 correlation_z_25_30->GetXaxis()->SetLabelSize(0.04);
741 correlation_z_25_30->GetXaxis()->SetTitleSize(0.07);
742 correlation_z_25_30->SetAxisRange(0, 0.6,
"Y");
743 correlation_z_25_30->SetAxisRange(0, 0.5,
"X");
744 correlation_z_25_30->SetXTitle(
"z_{g}");
746 correlation_z_25_30->SetTitle(
" ");
747 correlation_z_25_30->Draw(
"colz");
749 drawText(
"25 < p_{T} < 30 GeV/c", 0.3, 0.75, 12);
751 gStyle->SetErrorX(0);
755 correlation_z_30_40->SetMarkerStyle(29);
756 correlation_z_30_40->SetMarkerColor(2);
757 correlation_z_30_40->SetMarkerSize(1.0);
758 correlation_z_30_40->GetXaxis()->SetTitleSize(0.07);
759 correlation_z_30_40->GetXaxis()->CenterTitle(
true);
760 correlation_z_30_40->GetXaxis()->SetLabelSize(0.04);
761 correlation_z_30_40->GetXaxis()->SetTitleSize(0.07);
762 correlation_z_30_40->SetAxisRange(0, 0.6,
"Y");
763 correlation_z_30_40->SetAxisRange(0, 0.5,
"X");
764 correlation_z_30_40->SetXTitle(
"z_{g}");
766 correlation_z_30_40->SetTitle(
" ");
767 correlation_z_30_40->Draw(
"colz");
769 drawText(
"30 < p_{T} < 40 GeV/c", 0.3, 0.75, 12);
772 gStyle->SetErrorX(0);
777 correlation_theta_15_20->SetMarkerStyle(29);
778 correlation_theta_15_20->SetMarkerColor(2);
779 correlation_theta_15_20->SetMarkerSize(1.0);
780 correlation_theta_15_20->SetTitle(
" ");
781 correlation_theta_15_20->SetXTitle(
"#theta_{g}");
782 correlation_theta_15_20->GetXaxis()->SetTitleSize(0.07);
783 correlation_theta_15_20->GetXaxis()->CenterTitle(
true);
784 correlation_theta_15_20->GetXaxis()->SetLabelSize(0.04);
785 correlation_theta_15_20->GetYaxis()->SetTitleSize(0.07);
786 correlation_theta_15_20->SetAxisRange(0, 0.5,
"Y");
787 correlation_theta_15_20->SetAxisRange(0, 0.5,
"X");
788 correlation_theta_15_20->SetXTitle(
"R_{g}");
789 correlation_theta_15_20->SetYTitle(
"#theta_{sj}");
790 correlation_theta_15_20->GetYaxis()->SetLabelSize(0.04);
791 correlation_theta_15_20->SetTitle(
" ");
792 correlation_theta_15_20->Draw(
"colz");
795 gStyle->SetErrorX(0);
799 correlation_theta_20_25->SetMarkerStyle(29);
800 correlation_theta_20_25->SetMarkerColor(2);
801 correlation_theta_20_25->SetMarkerSize(1.0);
802 correlation_theta_20_25->SetTitle(
" ");
803 correlation_theta_20_25->SetXTitle(
"#theta_{g}");
804 correlation_theta_20_25->GetXaxis()->SetTitleSize(0.07);
805 correlation_theta_20_25->GetXaxis()->CenterTitle(
true);
806 correlation_theta_20_25->GetXaxis()->SetLabelSize(0.04);
807 correlation_theta_20_25->GetXaxis()->SetTitleSize(0.07);
808 correlation_theta_20_25->SetAxisRange(0, 0.5,
"Y");
809 correlation_theta_20_25->SetAxisRange(0, 0.5,
"X");
810 correlation_theta_20_25->SetXTitle(
"R_{g}");
812 correlation_theta_20_25->SetTitle(
" ");
813 correlation_theta_20_25->Draw(
"colz");
814 drawText(
"p+p #sqrt{s} = 200 GeV", 0.05, 0.9, 12);
817 gStyle->SetErrorX(0);
821 correlation_theta_25_30->SetMarkerStyle(29);
822 correlation_theta_25_30->SetMarkerColor(2);
823 correlation_theta_25_30->SetMarkerSize(1.0);
824 correlation_theta_25_30->SetTitle(
" ");
825 correlation_theta_25_30->SetXTitle(
"#theta_{g}");
826 correlation_theta_25_30->GetXaxis()->SetTitleSize(0.07);
827 correlation_theta_25_30->GetXaxis()->CenterTitle(
true);
828 correlation_theta_25_30->GetXaxis()->SetLabelSize(0.04);
829 correlation_theta_25_30->GetXaxis()->SetTitleSize(0.07);
830 correlation_theta_25_30->SetAxisRange(0, 0.5,
"Y");
831 correlation_theta_25_30->SetAxisRange(0, 0.5,
"X");
832 correlation_theta_25_30->SetXTitle(
"R_{g}");
834 correlation_theta_25_30->SetTitle(
" ");
835 correlation_theta_25_30->Draw(
"colz");
836 drawText(
"anti-k_{T} R=0.4 jets, |#eta_{jet}| < 0.6", 0.03, 0.9, 12);
839 gStyle->SetErrorX(0);
843 correlation_theta_30_40->SetMarkerStyle(29);
844 correlation_theta_30_40->SetMarkerColor(2);
845 correlation_theta_30_40->SetMarkerSize(1.0);
846 correlation_theta_30_40->SetTitle(
" ");
847 correlation_theta_30_40->GetXaxis()->SetTitleSize(0.07);
848 correlation_theta_30_40->GetXaxis()->CenterTitle(
true);
849 correlation_theta_30_40->GetXaxis()->SetLabelSize(0.04);
850 correlation_theta_30_40->SetAxisRange(0, 0.5,
"Y");
851 correlation_theta_30_40->SetAxisRange(0, 0.5,
"X");
852 correlation_theta_30_40->SetXTitle(
"R_{g}");
854 correlation_theta_30_40->SetTitle(
" ");
855 correlation_theta_30_40->Draw(
"colz");
856 drawText(
"#beta=0.0, z_{cut}=0.1", 0.05, 0.9, 12);
859 canvas->SaveAs(
"CorellBeta0z1.pdf",
"RECREATE");