Analysis Software
Documentation for sPHENIX simulation software
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Friends Macros Groups Pages
HIJING_SubJets.C
Go to the documentation of this file. Or view the newest version in sPHENIX GitHub for file HIJING_SubJets.C
1 #include "sPhenixStyle.h"
2 #include "sPhenixStyle.C"
3 #include <iostream>
4 #include "utilities.h"
5 #include <TGraphErrors.h>
6 #include <TRandom3.h>
7 
9 {
10 
11 
13  TH1::SetDefaultSumw2();
14  TH2::SetDefaultSumw2();
15 
16  using namespace std;{
17 
18  gStyle->SetOptStat(0);
19 
20 
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(); // Enable Sumw2 for this histogram
24  TH1F * _h_R04_z_g_20_25 = (TH1F*)fin->Get("R04_z_g_20_25");
25  _h_R04_z_g_20_25->Sumw2(); // Enable Sumw2 for this histogram
26  TH1F * _h_R04_z_g_25_30 = (TH1F*)fin->Get("R04_z_g_25_30");
27  _h_R04_z_g_25_30->Sumw2(); // Enable Sumw2 for this histogram
28  TH1F * _h_R04_z_g_30_40 = (TH1F*)fin->Get("R04_z_g_30_40");
29  _h_R04_z_g_30_40->Sumw2(); // Enable Sumw2 for this histogram
30  TH1F * _h_R04_theta_g_15_20 = (TH1F*)fin->Get("R04_theta_g_15_20");
31  _h_R04_theta_g_15_20->Sumw2(); // Enable Sumw2 for this histogram
32  TH1F * _h_R04_theta_g_20_25 = (TH1F*)fin->Get("R04_theta_g_20_25");
33  _h_R04_theta_g_20_25->Sumw2(); // Enable Sumw2 for this histogram
34  TH1F * _h_R04_theta_g_25_30 = (TH1F*)fin->Get("R04_theta_g_25_30");
35  _h_R04_theta_g_25_30->Sumw2(); // Enable Sumw2 for this histogram
36  TH1F * _h_R04_theta_g_30_40 = (TH1F*)fin->Get("R04_theta_g_30_40");
37  _h_R04_theta_g_30_40->Sumw2(); // Enable Sumw2 for this histogram
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(); // Enable Sumw2 for this histogram
41  TH1F * _h_R04_z_sj_20_25 = (TH1F*)fun->Get("R04_z_sj_20_25");
42  _h_R04_z_sj_20_25->Sumw2(); // Enable Sumw2 for this histogram
43  TH1F * _h_R04_z_sj_25_30 = (TH1F*)fun->Get("R04_z_sj_25_30");
44  _h_R04_z_sj_25_30->Sumw2(); // Enable Sumw2 for this histogram
45  TH1F * _h_R04_z_sj_30_40 = (TH1F*)fun->Get("R04_z_sj_30_40");
46  _h_R04_z_sj_30_40->Sumw2(); // Enable Sumw2 for this histogram
47  TH1F * _h_R04_theta_sj_15_20 = (TH1F*)fun->Get("R04_theta_sj_15_20");
48  _h_R04_theta_sj_15_20->Sumw2(); // Enable Sumw2 for this histogram
49  TH1F * _h_R04_theta_sj_20_25 = (TH1F*)fun->Get("R04_theta_sj_20_25");
50  _h_R04_theta_sj_20_25->Sumw2(); // Enable Sumw2 for this histogram
51  TH1F * _h_R04_theta_sj_25_30 = (TH1F*)fun->Get("R04_theta_sj_25_30");
52  _h_R04_theta_sj_25_30->Sumw2(); // Enable Sumw2 for this histogram
53  TH1F * _h_R04_theta_sj_30_40 = (TH1F*)fun->Get("R04_theta_sj_30_40");
54  _h_R04_theta_sj_30_40->Sumw2(); // Enable Sumw2 for this histogram
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");
59 
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");
64 
65  TCanvas * csubjetR01 = new TCanvas("csubjetR01", " ", 1500, 500);
66 
67  makeMultiPanelCanvasWithGap(csubjetR01, 4, 2, 0.2, 0.2, 0.2, 0.2, 0.08, 0.01);
68 
69 
70  // TH1F * hdummy_jetpT = (TH1F*) _h_R04_z_g_15_20->Clone("hdummy_jetpT");
71  // hdummy_jetpT->Reset();
72  // hdummy_jetpT->SetAxisRange(0, 5, "Y");
73 
74  csubjetR01->cd(1);
75  gStyle->SetErrorX(0);
76  gPad->SetTickx();
77  gPad->SetTicky();
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);
86 
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");
94  // _h_R04_z_sj_15_20 ->Draw("p E2 same");
95  _h_R04_z_sj_15_20->Draw("hist C same");
96  // Create a TGraphErrors for error bars
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);
105  ex15[i - 1] = 0.0; // Assumes uniform bin width
106  ey15[i - 1] = _h_R04_z_sj_15_20 ->GetBinError(i);
107  }
108  TGraphErrors* errorGraph15 = new TGraphErrors(nbins15, x15, y15, ex15, ey15);
109  errorGraph15->SetMarkerStyle(29);
110  errorGraph15->SetMarkerColor(2);
111  errorGraph15->Draw("p same");
112 
113 
114  csubjetR01->cd(2);
115  gStyle->SetErrorX(0);
116  gPad->SetTickx();
117  gPad->SetTicky();
118  // cjetSpectra->cd(1)->SetLogy();
119  // hdummy_jetpT->Draw();
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);
128  // _h_R04_z_sj_20_25->Draw("p E1 same");
129  _h_R04_z_sj_20_25->Draw("hist C same");
130 
131  // Create a TGraphErrors for error bars
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);
140  ex20[i - 1] = 0.0; // Assumes uniform bin width
141  ey20[i - 1] = _h_R04_z_sj_20_25->GetBinError(i);
142  }
143  TGraphErrors* errorGraph20 = new TGraphErrors(nbins20, x20, y20, ex20, ey20);
144  errorGraph20->SetMarkerStyle(29);
145  errorGraph20->SetMarkerColor(2);
146  errorGraph20->Draw("p same");
147 
148  csubjetR01->cd(3);
149  gStyle->SetErrorX(0);
150  gPad->SetTickx();
151  gPad->SetTicky();
152  // cjetSpectra->cd(1)->SetLogy();
153  // hdummy_jetpT->Draw();
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);
164  // _h_R04_z_sj_25_30->Draw("p E1 same");
165  _h_R04_z_sj_25_30->Draw("hist C same");
166  // Create a TGraphErrors for error bars
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);
175  ex25[i - 1] = 0.0; // Assumes uniform bin width
176  ey25[i - 1] = _h_R04_z_sj_25_30->GetBinError(i);
177  }
178  TGraphErrors* errorGraph25 = new TGraphErrors(nbins25, x25, y25, ex25, ey25);
179  errorGraph25->SetMarkerStyle(29);
180  errorGraph25->SetMarkerColor(2);
181  errorGraph25->Draw("p same");
182 
183  csubjetR01->cd(4);
184  gStyle->SetErrorX(0);
185  gPad->SetTickx();
186  gPad->SetTicky();
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);
198  // _h_R04_z_sj_30_40->Draw("p E1 same");
199  _h_R04_z_sj_30_40->Draw("hist C same");
200 
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);
204 
205  // Create a TGraphErrors for error bars
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);
214  ex30[i - 1] = 0.0; // Assumes uniform bin width
215  ey30[i - 1] = _h_R04_z_sj_30_40->GetBinError(i);
216  }
217  TGraphErrors* errorGraph30 = new TGraphErrors(nbins15, x30, y30, ex30, ey30);
218  errorGraph30->SetMarkerStyle(29);
219  errorGraph30->SetMarkerColor(2);
220  errorGraph30->Draw("p same");
221 
222  csubjetR01->cd(5);
223  gStyle->SetErrorX(0);
224  gPad->SetTickx();
225  gPad->SetTicky();
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);
237 
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);
242 
243  _h_R04_theta_sj_15_20->Draw("hist C same");
244  // Create a TGraphErrors for error bars
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);
253  ex15t[i - 1] = 0.0; // Assumes uniform bin width
254  ey15t[i - 1] = _h_R04_theta_sj_15_20->GetBinError(i);
255  }
256  TGraphErrors* errorGraph15t = new TGraphErrors(nbins15t, x15t, y15t, ex15t, ey15t);
257  errorGraph15t->SetMarkerStyle(29);
258  errorGraph15t->SetMarkerColor(2);
259  errorGraph15t->Draw("p same");
260 
261  drawText("15 < p_{T} < 20 [GeV/c]", 0.6, 0.8, 12);
262 
263  csubjetR01->cd(6);
264  gStyle->SetErrorX(0);
265  gPad->SetTickx();
266  gPad->SetTicky();
267  // cjetSpectra->cd(5)->SetLogy();
268  // hdummy_jetpT->Draw();
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");
280  // _h_R04_theta_sj_20_25->Draw("p E1 same");
281  // Create a TGraphErrors for error bars
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);
290  ex20t[i - 1] = 0.0; // Assumes uniform bin width
291  ey20t[i - 1] = _h_R04_theta_sj_20_25->GetBinError(i);
292  }
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);
298 
299  csubjetR01->cd(7);
300  gStyle->SetErrorX(0);
301  gPad->SetTickx();
302  gPad->SetTicky();
303  // cjetSpectra->cd(5)->SetLogy();
304  // hdummy_jetpT->Draw();
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");
315  // _h_R04_theta_sj_25_30->Draw("p E1 same");
316 
317  _h_R04_theta_sj_25_30->Draw("hist C same");
318  // Create a TGraphErrors for error bars
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);
327  ex25t[i - 1] = 0.0; // Assumes uniform bin width
328  ey25t[i - 1] = _h_R04_theta_sj_25_30->GetBinError(i);
329  }
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);
335 
336  csubjetR01->cd(8);
337  gStyle->SetErrorX(0);
338  gPad->SetTickx();
339  gPad->SetTicky();
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");
351  // _h_R04_theta_sj_30_40->Draw("p E1 same");
352  _h_R04_theta_sj_30_40->Draw("hist C same");
353  // Create a TGraphErrors for error bars
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);
362  ex30t[i - 1] = 0.0; // Assumes uniform bin width
363  ey30t[i - 1] = _h_R04_theta_sj_30_40->GetBinError(i);
364  }
365  TGraphErrors* errorGraph30t = new TGraphErrors(nbins30t, x30t, y30t, ex30t, ey30t);
366  errorGraph30t->SetMarkerStyle(29);
367  errorGraph30t->SetMarkerColor(2);
368  errorGraph30t->Draw("p same");
369 
370  drawText("30 < p_{T} < 40 [GeV/c]", 0.5, 0.8, 12);
371 
372  csubjetR01->SaveAs("RHIC_subjetR04_AKTB0z1.pdf", "RECREATE");
373 
374 
375  TCanvas * csubjetR04 = new TCanvas("csubjetR04", " ", 1500, 500);
376 
377  makeMultiPanelCanvasWithGap(csubjetR04, 4, 2, 0.2, 0.2, 0.2, 0.2, 0.08, 0.01);
378 
379  csubjetR04->cd(1);
380  gStyle->SetErrorX(0);
381  gPad->SetTickx();
382  gPad->SetTicky();
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);
391 
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");
399  // _h_R04_z_g_15_20->Draw("p E2 same");
400  _h_R04_z_g_15_20->Draw("hist C same");
401  // Create a TGraphErrors for error bars
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);
410  ex15zg[i - 1] = 0.0; // Assumes uniform bin width
411  ey15zg[i - 1] = _h_R04_z_g_15_20->GetBinError(i);
412  }
413  TGraphErrors* errorGraph15zg = new TGraphErrors(nbins15zg, x15zg, y15zg, ex15zg, ey15zg);
414  errorGraph15zg->SetMarkerStyle(29);
415  errorGraph15zg->SetMarkerColor(2);
416  errorGraph15zg->Draw("p same");
417 
418 
419  csubjetR04->cd(2);
420  gStyle->SetErrorX(0);
421  gPad->SetTickx();
422  gPad->SetTicky();
423  // cjetSpectra->cd(1)->SetLogy();
424  // hdummy_jetpT->Draw();
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);
433  // _h_R04_z_g_20_25->Draw("p E1 same");
434  _h_R04_z_g_20_25->Draw("hist C same");
435 
436  // Create a TGraphErrors for error bars
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);
445  ex20zg[i - 1] = 0.0; // Assumes uniform bin width
446  ey20zg[i - 1] = _h_R04_z_g_20_25->GetBinError(i);
447  }
448  TGraphErrors* errorGraph20zg = new TGraphErrors(nbins20, x20zg, y20zg, ex20zg, ey20zg);
449  errorGraph20zg->SetMarkerStyle(29);
450  errorGraph20zg->SetMarkerColor(2);
451  errorGraph20zg->Draw("p same");
452 
453  csubjetR04->cd(3);
454  gStyle->SetErrorX(0);
455  gPad->SetTickx();
456  gPad->SetTicky();
457  // cjetSpectra->cd(1)->SetLogy();
458  // hdummy_jetpT->Draw();
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);
469  // _h_R04_z_g_25_30->Draw("p E1 same");
470  _h_R04_z_g_25_30->Draw("hist C same");
471  // Create a TGraphErrors for error bars
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);
480  ex25zg[i - 1] = 0.0; // Assumes uniform bin width
481  ey25zg[i - 1] = _h_R04_z_g_25_30->GetBinError(i);
482  }
483  TGraphErrors* errorGraph25zg = new TGraphErrors(nbins25zg, x25zg, y25zg, ex25zg, ey25zg);
484  errorGraph25zg->SetMarkerStyle(29);
485  errorGraph25zg->SetMarkerColor(2);
486  errorGraph25zg->Draw("p same");
487 
488  csubjetR04->cd(4);
489  gStyle->SetErrorX(0);
490  gPad->SetTickx();
491  gPad->SetTicky();
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);
503  // _h_R04_z_g_30_40->Draw("p E1 same");
504  _h_R04_z_g_30_40->Draw("hist C same");
505 
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);
509 
510  // Create a TGraphErrors for error bars
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);
519  ex30zg[i - 1] = 0.0; // Assumes uniform bin width
520  ey30zg[i - 1] = _h_R04_z_g_30_40->GetBinError(i);
521  }
522  TGraphErrors* errorGraph30zg = new TGraphErrors(nbins15zg, x30zg, y30zg, ex30zg, ey30zg);
523  errorGraph30zg->SetMarkerStyle(29);
524  errorGraph30zg->SetMarkerColor(2);
525  errorGraph30zg->Draw("p same");
526 
527  csubjetR04->cd(5);
528  gStyle->SetErrorX(0);
529  gPad->SetTickx();
530  gPad->SetTicky();
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);
542 
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);
547 
548  _h_R04_theta_g_15_20->Draw("hist C same");
549  // Create a TGraphErrors for error bars
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);
558  ex15gt[i - 1] = 0.0; // Assumes uniform bin width
559  ey15gt[i - 1] = _h_R04_theta_g_15_20->GetBinError(i);
560  }
561  TGraphErrors* errorGraph15gt = new TGraphErrors(nbins15tgt, x15gt, y15gt, ex15gt, ey15gt);
562  errorGraph15gt->SetMarkerStyle(29);
563  errorGraph15gt->SetMarkerColor(2);
564  errorGraph15gt->Draw("p same");
565 
566  drawText("15 < p_{T} < 20 GeV/c", 0.7, 0.8, 12);
567 
568  csubjetR04->cd(6);
569  gStyle->SetErrorX(0);
570  gPad->SetTickx();
571  gPad->SetTicky();
572  // cjetSpectra->cd(5)->SetLogy();
573  // hdummy_jetpT->Draw();
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");
585  // _h_R04_theta_g_20_25->Draw("p E1 same");
586  // Create a TGraphErrors for error bars
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);
595  ex20gt[i - 1] =0.0; // Assumes uniform bin width
596  ey20gt[i - 1] = _h_R04_theta_g_20_25->GetBinError(i);
597  }
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);
603 
604  csubjetR04->cd(7);
605  gStyle->SetErrorX(0);
606  gPad->SetTickx();
607  gPad->SetTicky();
608  // cjetSpectra->cd(5)->SetLogy();
609  // hdummy_jetpT->Draw();
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");
620  // _h_R04_theta_g_25_30->Draw("p E1 same");
621 
622  _h_R04_theta_g_25_30->Draw("hist C same");
623  // Create a TGraphErrors for error bars
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);
632  ex25gt[i - 1] = 0.0;
633  //ex25t[i - 1] = _h_R04_theta_g_25_30->GetBinWidth(i) / 2.0; // Assumes uniform bin width
634  ey25gt[i - 1] = _h_R04_theta_g_25_30->GetBinError(i);
635  }
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);
641 
642  csubjetR04->cd(8);
643  gStyle->SetErrorX(0);
644  gPad->SetTickx();
645  gPad->SetTicky();
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");
657  // _h_R04_theta_g_30_40->Draw("p E1 same");
658  _h_R04_theta_g_30_40->Draw("hist C same");
659  // Create a TGraphErrors for error bars
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);
668  // ex30tt[i - 1] = _h_R04_theta_g_30_40->GetBinWidth(i) / 2.0; // Assumes uniform bin width
669  ex30gt[i - 1] = 0.0;
670  ey30gt[i - 1] = _h_R04_theta_g_30_40->GetBinError(i);
671  }
672  TGraphErrors* errorGraph30gt = new TGraphErrors(nbins20tgt, x30gt, y30gt, ex30gt, ey30gt);
673  errorGraph30gt->SetMarkerStyle(29);
674  errorGraph30gt->SetMarkerColor(2);
675  errorGraph30gt->Draw("p same");
676 
677  drawText("30 < p_{T} < 40 GeV/c", 0.55, 0.8, 12);
678 
679  csubjetR04->SaveAs("RHIC_subjetR04_CAB0z1.pdf", "RECREATE");
680 
681 
682  TCanvas * canvas = new TCanvas("canvas", " ", 1500, 500);
683 
684  makeMultiPanelCanvasWithGap(canvas, 4, 2, 0.2, 0.2, 0.2, 0.2, 0.09, 0.01);
685 
686  canvas->cd(1);
687  gStyle->SetErrorX(0);
688  gPad->SetTickx();
689  gPad->SetTicky();
690 
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"); // Use "colz" option to draw the color map
706  drawText("15 < p_{T} < 20 GeV/c", 0.4, 0.75, 12);
707 
708 
709  canvas->cd(2);
710  gStyle->SetErrorX(0);
711  gPad->SetTickx();
712  gPad->SetTicky();
713 
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}");
724  // correlation_z_20_25->SetYTitle("z_{sj}");
725  correlation_z_20_25->SetTitle(" ");
726  correlation_z_20_25->Draw("colz"); // Use "colz" option to draw the color map
727  drawText("20 < p_{T} < 25 GeV/c", 0.3, 0.75, 12);
728 
729 
730  canvas->cd(3);
731  gStyle->SetErrorX(0);
732  gPad->SetTickx();
733  gPad->SetTicky();
734 
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}");
745  // correlation_z_25_30->SetYTitle("z_{sj}");
746  correlation_z_25_30->SetTitle(" ");
747  correlation_z_25_30->Draw("colz"); // Use "colz" option to draw the color map
748 
749  drawText("25 < p_{T} < 30 GeV/c", 0.3, 0.75, 12);
750  canvas->cd(4);
751  gStyle->SetErrorX(0);
752  gPad->SetTickx();
753  gPad->SetTicky();
754 
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}");
765  // correlation_z_30_40->SetYTitle("z_{sj}");
766  correlation_z_30_40->SetTitle(" ");
767  correlation_z_30_40->Draw("colz"); // Use "colz" option to draw the color map
768 
769  drawText("30 < p_{T} < 40 GeV/c", 0.3, 0.75, 12);
770 
771  canvas->cd(5);
772  gStyle->SetErrorX(0);
773  gPad->SetTickx();
774  gPad->SetTicky();
775 
776 
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"); // Use "colz" option to draw the color map
793 
794  canvas->cd(6);
795  gStyle->SetErrorX(0);
796  gPad->SetTickx();
797  gPad->SetTicky();
798 
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}");
811  // correlation_theta_20_25->SetYTitle("#theta_{sj}");
812  correlation_theta_20_25->SetTitle(" ");
813  correlation_theta_20_25->Draw("colz"); // Use "colz" option to draw the color map
814  drawText("p+p #sqrt{s} = 200 GeV", 0.05, 0.9, 12);
815 
816  canvas->cd(7);
817  gStyle->SetErrorX(0);
818  gPad->SetTickx();
819  gPad->SetTicky();
820 
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}");
833  // correlation_theta_25_30->SetYTitle("#theta_{sj}");
834  correlation_theta_25_30->SetTitle(" ");
835  correlation_theta_25_30->Draw("colz"); // Use "colz" option to draw the color map
836  drawText("anti-k_{T} R=0.4 jets, |#eta_{jet}| < 0.6", 0.03, 0.9, 12);
837 
838  canvas->cd(8);
839  gStyle->SetErrorX(0);
840  gPad->SetTickx();
841  gPad->SetTicky();
842 
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}");
853  // correlation_theta_30_40->SetYTitle("#theta_{sj}");
854  correlation_theta_30_40->SetTitle(" ");
855  correlation_theta_30_40->Draw("colz"); // Use "colz" option to draw the color map
856  drawText("#beta=0.0, z_{cut}=0.1", 0.05, 0.9, 12);
857 
858 
859  canvas->SaveAs("CorellBeta0z1.pdf", "RECREATE");
860 
861  }
862 
863 }
864