Back to home page

sPhenix code displayed by LXR
 
 

    

File indexing completed on 2025-08-05 08:13:19 

0001 #include "sPhenixStyle.h"
0002 #include "sPhenixStyle.C"
0003 #include <iostream>
0004 #include "utilities.h"
0005 #include <TGraphErrors.h>
0006 #include <TRandom3.h>
0007 
0008 void HIJING_SubJets() 
0009 {
0010 
0011 
0012   SetsPhenixStyle();
0013   TH1::SetDefaultSumw2();
0014   TH2::SetDefaultSumw2();
0015 
0016   using namespace std;{
0017 
0018     gStyle->SetOptStat(0);
0019 
0020   
0021     TFile * fin = TFile::Open("../macro/hist_jets.root");  
0022     TH1F * _h_R04_z_g_15_20 = (TH1F*)fin->Get("R04_z_g_15_20");
0023     _h_R04_z_g_15_20->Sumw2(); // Enable Sumw2 for this histogram
0024     TH1F * _h_R04_z_g_20_25 = (TH1F*)fin->Get("R04_z_g_20_25");
0025     _h_R04_z_g_20_25->Sumw2(); // Enable Sumw2 for this histogram
0026     TH1F * _h_R04_z_g_25_30 = (TH1F*)fin->Get("R04_z_g_25_30");
0027     _h_R04_z_g_25_30->Sumw2(); // Enable Sumw2 for this histogram
0028     TH1F * _h_R04_z_g_30_40 = (TH1F*)fin->Get("R04_z_g_30_40");
0029     _h_R04_z_g_30_40->Sumw2(); // Enable Sumw2 for this histogram
0030     TH1F * _h_R04_theta_g_15_20 = (TH1F*)fin->Get("R04_theta_g_15_20");
0031     _h_R04_theta_g_15_20->Sumw2(); // Enable Sumw2 for this histogram
0032     TH1F * _h_R04_theta_g_20_25 = (TH1F*)fin->Get("R04_theta_g_20_25");
0033     _h_R04_theta_g_20_25->Sumw2(); // Enable Sumw2 for this histogram
0034     TH1F * _h_R04_theta_g_25_30 = (TH1F*)fin->Get("R04_theta_g_25_30");
0035     _h_R04_theta_g_25_30->Sumw2(); // Enable Sumw2 for this histogram
0036     TH1F * _h_R04_theta_g_30_40 = (TH1F*)fin->Get("R04_theta_g_30_40");
0037     _h_R04_theta_g_30_40->Sumw2(); // Enable Sumw2 for this histogram
0038     TFile * fun = TFile::Open("../macro/hist_jets.root");  
0039     TH1F * _h_R04_z_sj_15_20 = (TH1F*)fun->Get("R04_z_sj_15_20");
0040     _h_R04_z_sj_15_20->Sumw2(); // Enable Sumw2 for this histogram
0041     TH1F * _h_R04_z_sj_20_25 = (TH1F*)fun->Get("R04_z_sj_20_25");
0042     _h_R04_z_sj_20_25->Sumw2(); // Enable Sumw2 for this histogram
0043     TH1F * _h_R04_z_sj_25_30 = (TH1F*)fun->Get("R04_z_sj_25_30");
0044     _h_R04_z_sj_25_30->Sumw2(); // Enable Sumw2 for this histogram
0045     TH1F * _h_R04_z_sj_30_40 = (TH1F*)fun->Get("R04_z_sj_30_40");
0046     _h_R04_z_sj_30_40->Sumw2(); // Enable Sumw2 for this histogram
0047     TH1F * _h_R04_theta_sj_15_20 = (TH1F*)fun->Get("R04_theta_sj_15_20");
0048     _h_R04_theta_sj_15_20->Sumw2(); // Enable Sumw2 for this histogram
0049     TH1F * _h_R04_theta_sj_20_25 = (TH1F*)fun->Get("R04_theta_sj_20_25");
0050     _h_R04_theta_sj_20_25->Sumw2(); // Enable Sumw2 for this histogram
0051     TH1F * _h_R04_theta_sj_25_30 = (TH1F*)fun->Get("R04_theta_sj_25_30");
0052     _h_R04_theta_sj_25_30->Sumw2(); // Enable Sumw2 for this histogram
0053     TH1F * _h_R04_theta_sj_30_40 = (TH1F*)fun->Get("R04_theta_sj_30_40");
0054     _h_R04_theta_sj_30_40->Sumw2(); // Enable Sumw2 for this histogram
0055     TH2D * correlation_theta_15_20 = (TH2D*)fun->Get("correlation_theta_15_20");
0056     TH2D * correlation_theta_20_25 = (TH2D*)fun->Get("correlation_theta_20_25");
0057     TH2D * correlation_theta_25_30 = (TH2D*)fun->Get("correlation_theta_25_30");
0058     TH2D * correlation_theta_30_40 = (TH2D*)fun->Get("correlation_theta_30_40");
0059 
0060     TH2D * correlation_z_15_20 = (TH2D*)fun->Get("correlation_z_15_20");
0061     TH2D * correlation_z_20_25 = (TH2D*)fun->Get("correlation_z_20_25");
0062     TH2D * correlation_z_25_30 = (TH2D*)fun->Get("correlation_z_25_30");
0063     TH2D * correlation_z_30_40 = (TH2D*)fun->Get("correlation_z_30_40");
0064 
0065     TCanvas * csubjetR01 = new TCanvas("csubjetR01", " ", 1500, 500);
0066   
0067     makeMultiPanelCanvasWithGap(csubjetR01, 4, 2, 0.2, 0.2, 0.2, 0.2, 0.08, 0.01);
0068  
0069 
0070     // TH1F * hdummy_jetpT = (TH1F*) _h_R04_z_g_15_20->Clone("hdummy_jetpT");
0071     // hdummy_jetpT->Reset();
0072     // hdummy_jetpT->SetAxisRange(0, 5, "Y");
0073   
0074     csubjetR01->cd(1);
0075     gStyle->SetErrorX(0);
0076     gPad->SetTickx();
0077     gPad->SetTicky();
0078     _h_R04_z_sj_15_20->SetLineStyle(1);
0079     _h_R04_z_sj_15_20->SetLineWidth(1);
0080     _h_R04_z_sj_15_20->SetTitle(" ");
0081     _h_R04_z_sj_15_20->SetYTitle("(1/N_{jet}) dN/dz_{sj}");
0082     _h_R04_z_sj_15_20->GetYaxis()->SetTitleOffset(0.8);
0083     _h_R04_z_sj_15_20->GetYaxis()->SetTitleSize(0.07);
0084     _h_R04_z_sj_15_20->GetYaxis()->SetLabelSize(0.04);
0085     _h_R04_z_sj_15_20->GetYaxis()->CenterTitle(true);
0086 
0087     _h_R04_z_sj_15_20->SetXTitle("z_{sj}");
0088     _h_R04_z_sj_15_20->GetXaxis()->SetTitleSize(0.07);
0089     _h_R04_z_sj_15_20->GetXaxis()->CenterTitle(true);
0090     _h_R04_z_sj_15_20->GetXaxis()->SetLabelSize(0.04);
0091     _h_R04_z_sj_15_20->SetMarkerSize(1.0);
0092     _h_R04_z_sj_15_20->SetAxisRange(0, 5.5, "Y");
0093     _h_R04_z_sj_15_20->SetAxisRange(0, 0.5, "X");
0094     //  _h_R04_z_sj_15_20 ->Draw("p E2 same");
0095     _h_R04_z_sj_15_20->Draw("hist C same");
0096     // Create a TGraphErrors for error bars
0097     int nbins15 = _h_R04_z_sj_15_20->GetNbinsX();
0098     double* x15 = new double[nbins15];
0099     double* y15 = new double[nbins15];
0100     double* ex15 = new double[nbins15];
0101     double* ey15 = new double[nbins15];
0102     for (int i = 1; i <= nbins15; ++i) {
0103       x15[i - 1] = _h_R04_z_sj_15_20 ->GetBinCenter(i);
0104       y15[i - 1] = _h_R04_z_sj_15_20 ->GetBinContent(i);
0105       ex15[i - 1] = 0.0; // Assumes uniform bin width
0106       ey15[i - 1] = _h_R04_z_sj_15_20 ->GetBinError(i);
0107     }
0108     TGraphErrors* errorGraph15 = new TGraphErrors(nbins15, x15, y15, ex15, ey15);
0109     errorGraph15->SetMarkerStyle(29);
0110     errorGraph15->SetMarkerColor(2);
0111     errorGraph15->Draw("p same");
0112     
0113     
0114     csubjetR01->cd(2);
0115     gStyle->SetErrorX(0);
0116     gPad->SetTickx();
0117     gPad->SetTicky();
0118     // cjetSpectra->cd(1)->SetLogy();
0119     // hdummy_jetpT->Draw();
0120     _h_R04_z_sj_20_25->SetTitle(" ");
0121     _h_R04_z_sj_20_25->SetXTitle("z_{sj}");
0122     _h_R04_z_sj_20_25->GetXaxis()->SetTitleSize(0.07);
0123     _h_R04_z_sj_20_25->GetXaxis()->CenterTitle(true);
0124     _h_R04_z_sj_20_25->GetXaxis()->SetLabelSize(0.04);
0125     _h_R04_z_sj_20_25->SetAxisRange(0, 5.5, "Y");
0126     _h_R04_z_sj_20_25->SetAxisRange(0, 0.5, "X");
0127     _h_R04_z_sj_20_25->SetMarkerSize(1.0);
0128     // _h_R04_z_sj_20_25->Draw("p E1 same");
0129     _h_R04_z_sj_20_25->Draw("hist C same");
0130 
0131     // Create a TGraphErrors for error bars
0132     int nbins20 = _h_R04_z_sj_20_25->GetNbinsX();
0133     double* x20 = new double[nbins20];
0134     double* y20 = new double[nbins20];
0135     double* ex20 = new double[nbins20];
0136     double* ey20 = new double[nbins20];
0137     for (int i = 1; i <= nbins20; ++i) {
0138       x20[i - 1] = _h_R04_z_sj_20_25->GetBinCenter(i);
0139       y20[i - 1] = _h_R04_z_sj_20_25->GetBinContent(i);
0140       ex20[i - 1] = 0.0; // Assumes uniform bin width
0141       ey20[i - 1] = _h_R04_z_sj_20_25->GetBinError(i);
0142     }
0143     TGraphErrors* errorGraph20 = new TGraphErrors(nbins20, x20, y20, ex20, ey20);
0144     errorGraph20->SetMarkerStyle(29);
0145     errorGraph20->SetMarkerColor(2);
0146     errorGraph20->Draw("p same");
0147 
0148     csubjetR01->cd(3);
0149     gStyle->SetErrorX(0);
0150     gPad->SetTickx();
0151     gPad->SetTicky();
0152     // cjetSpectra->cd(1)->SetLogy();
0153     // hdummy_jetpT->Draw();
0154     _h_R04_z_sj_25_30->SetMarkerStyle(29);
0155     _h_R04_z_sj_25_30->SetMarkerColor(2);
0156     _h_R04_z_sj_25_30->SetTitle(" ");
0157     _h_R04_z_sj_25_30->SetXTitle("z_{sj}");
0158     _h_R04_z_sj_25_30->GetXaxis()->SetTitleSize(0.07);
0159     _h_R04_z_sj_25_30->GetXaxis()->CenterTitle(true);
0160     _h_R04_z_sj_25_30->GetXaxis()->SetLabelSize(0.04);
0161     _h_R04_z_sj_25_30->SetAxisRange(0, 5.5, "Y");
0162     _h_R04_z_sj_25_30->SetAxisRange(0, 0.5, "X");
0163     _h_R04_z_sj_25_30->SetMarkerSize(1.0);
0164     //  _h_R04_z_sj_25_30->Draw("p E1 same");
0165     _h_R04_z_sj_25_30->Draw("hist C same");
0166     // Create a TGraphErrors for error bars
0167     int nbins25 = _h_R04_z_sj_25_30->GetNbinsX();
0168     double* x25 = new double[nbins25];
0169     double* y25 = new double[nbins25];
0170     double* ex25 = new double[nbins25];
0171     double* ey25 = new double[nbins25];
0172     for (int i = 1; i <= nbins25; ++i) {
0173       x25[i - 1] = _h_R04_z_sj_25_30->GetBinCenter(i);
0174       y25[i - 1] = _h_R04_z_sj_25_30->GetBinContent(i);
0175       ex25[i - 1] = 0.0; // Assumes uniform bin width
0176       ey25[i - 1] = _h_R04_z_sj_25_30->GetBinError(i);
0177     }
0178     TGraphErrors* errorGraph25 = new TGraphErrors(nbins25, x25, y25, ex25, ey25);
0179     errorGraph25->SetMarkerStyle(29);
0180     errorGraph25->SetMarkerColor(2);
0181     errorGraph25->Draw("p same");
0182 
0183     csubjetR01->cd(4);
0184     gStyle->SetErrorX(0);
0185     gPad->SetTickx();
0186     gPad->SetTicky();
0187     _h_R04_z_sj_30_40->SetMarkerStyle(29);
0188     _h_R04_z_sj_30_40->SetMarkerColor(2);
0189     _h_R04_z_sj_30_40->SetTitle(" ");
0190     _h_R04_z_sj_30_40->SetXTitle("z_{sj}");
0191     _h_R04_z_sj_30_40->GetXaxis()->SetTitleSize(0.07);
0192     _h_R04_z_sj_30_40->GetXaxis()->CenterTitle(true);
0193     _h_R04_z_sj_30_40->GetXaxis()->SetLabelSize(0.04);
0194     _h_R04_z_sj_30_40->GetXaxis()->SetTitleSize(0.07);
0195     _h_R04_z_sj_30_40->SetAxisRange(0, 5.5, "Y");
0196     _h_R04_z_sj_30_40->SetAxisRange(0, 0.5, "X");
0197     _h_R04_z_sj_30_40->SetMarkerSize(1.0);
0198     //  _h_R04_z_sj_30_40->Draw("p E1 same");
0199     _h_R04_z_sj_30_40->Draw("hist C same");
0200   
0201     drawText("p+p #sqrt{s} = 200 GeV", 0.3, 0.8, 14);
0202     drawText("anti-k_{T} R=0.4 jets, |#eta_{jet}| < 0.6", 0.3, 0.7, 14);
0203     drawText("R_{sj}=0.1", 0.3, 0.6, 14);
0204 
0205     // Create a TGraphErrors for error bars
0206     int nbins30 = _h_R04_z_sj_30_40->GetNbinsX();
0207     double* x30 = new double[nbins30];
0208     double* y30 = new double[nbins30];
0209     double* ex30 = new double[nbins30];
0210     double* ey30 = new double[nbins30];
0211     for (int i = 1; i <= nbins30; ++i) {
0212       x30[i - 1] = _h_R04_z_sj_30_40->GetBinCenter(i);
0213       y30[i - 1] = _h_R04_z_sj_30_40->GetBinContent(i);
0214       ex30[i - 1] = 0.0; // Assumes uniform bin width
0215       ey30[i - 1] = _h_R04_z_sj_30_40->GetBinError(i);
0216     }
0217     TGraphErrors* errorGraph30 = new TGraphErrors(nbins15, x30, y30, ex30, ey30);
0218     errorGraph30->SetMarkerStyle(29);
0219     errorGraph30->SetMarkerColor(2);
0220     errorGraph30->Draw("p same");
0221 
0222     csubjetR01->cd(5);
0223     gStyle->SetErrorX(0);
0224     gPad->SetTickx();
0225     gPad->SetTicky();
0226     _h_R04_theta_sj_15_20->SetMarkerStyle(29);
0227     _h_R04_theta_sj_15_20->SetMarkerColor(2);
0228     _h_R04_theta_sj_15_20->SetMarkerSize(1.0);
0229     _h_R04_theta_sj_15_20->SetTitle(" ");
0230     _h_R04_theta_sj_15_20->SetYTitle("(1/N_{jet}) dN/d#theta_{sj}");
0231     _h_R04_theta_sj_15_20->SetAxisRange(0, 9.0, "Y");
0232     _h_R04_theta_sj_15_20->SetAxisRange(0, 0.5, "X");
0233     _h_R04_theta_sj_15_20->GetYaxis()->SetTitleOffset(0.8);
0234     _h_R04_theta_sj_15_20->GetYaxis()->SetTitleSize(0.07);
0235     _h_R04_theta_sj_15_20->GetYaxis()->SetLabelSize(0.04);
0236     _h_R04_theta_sj_15_20->GetYaxis()->CenterTitle(true);
0237 
0238     _h_R04_theta_sj_15_20->SetXTitle("#theta_{sj}");
0239     _h_R04_theta_sj_15_20->GetXaxis()->SetTitleSize(0.07);
0240     _h_R04_theta_sj_15_20->GetXaxis()->CenterTitle(true);
0241     _h_R04_theta_sj_15_20->GetXaxis()->SetLabelSize(0.04);
0242 
0243     _h_R04_theta_sj_15_20->Draw("hist C same");
0244     // Create a TGraphErrors for error bars
0245     int nbins15t = _h_R04_theta_sj_15_20->GetNbinsX();
0246     double* x15t = new double[nbins15t];
0247     double* y15t = new double[nbins15t];
0248     double* ex15t = new double[nbins15t];
0249     double* ey15t = new double[nbins15t];
0250     for (int i = 1; i <= nbins15t; ++i) {
0251       x15t[i - 1] = _h_R04_theta_sj_15_20->GetBinCenter(i);
0252       y15t[i - 1] = _h_R04_theta_sj_15_20->GetBinContent(i);
0253       ex15t[i - 1] = 0.0; // Assumes uniform bin width
0254       ey15t[i - 1] = _h_R04_theta_sj_15_20->GetBinError(i);
0255     }
0256     TGraphErrors* errorGraph15t = new TGraphErrors(nbins15t, x15t, y15t, ex15t, ey15t);
0257     errorGraph15t->SetMarkerStyle(29);
0258     errorGraph15t->SetMarkerColor(2);
0259     errorGraph15t->Draw("p same");
0260 
0261     drawText("15 < p_{T} < 20 [GeV/c]", 0.6, 0.8, 12);
0262     
0263     csubjetR01->cd(6);
0264     gStyle->SetErrorX(0);
0265     gPad->SetTickx();
0266     gPad->SetTicky();
0267     // cjetSpectra->cd(5)->SetLogy();
0268     // hdummy_jetpT->Draw();
0269     _h_R04_theta_sj_20_25->SetMarkerStyle(29);
0270     _h_R04_theta_sj_20_25->SetMarkerColor(2);
0271     _h_R04_theta_sj_20_25->SetMarkerSize(1.0);
0272     _h_R04_theta_sj_20_25->SetTitle(" ");
0273     _h_R04_theta_sj_20_25->SetXTitle("#theta_{sj}");
0274     _h_R04_theta_sj_20_25->GetXaxis()->SetTitleSize(0.07);
0275     _h_R04_theta_sj_20_25->GetXaxis()->CenterTitle(true);
0276     _h_R04_theta_sj_20_25->GetXaxis()->SetLabelSize(0.04);
0277     _h_R04_theta_sj_20_25->SetAxisRange(0, 9.0, "Y");
0278     _h_R04_theta_sj_20_25->SetAxisRange(0, 0.5, "X");
0279     _h_R04_theta_sj_20_25->Draw("hist C same");
0280     //  _h_R04_theta_sj_20_25->Draw("p E1 same");
0281     // Create a TGraphErrors for error bars
0282     int nbins20t = _h_R04_theta_sj_20_25->GetNbinsX();
0283     double* x20t = new double[nbins20t];
0284     double* y20t = new double[nbins20t];
0285     double* ex20t = new double[nbins20t];
0286     double* ey20t = new double[nbins20t];
0287     for (int i = 1; i <= nbins15t; ++i) {
0288       x20t[i - 1] = _h_R04_theta_sj_20_25->GetBinCenter(i);
0289       y20t[i - 1] = _h_R04_theta_sj_20_25->GetBinContent(i);
0290       ex20t[i - 1] = 0.0; // Assumes uniform bin width
0291       ey20t[i - 1] = _h_R04_theta_sj_20_25->GetBinError(i);
0292     }
0293     TGraphErrors* errorGraph20t = new TGraphErrors(nbins20t, x20t, y20t, ex20t, ey20t);
0294     errorGraph20t->SetMarkerStyle(29);
0295     errorGraph20t->SetMarkerColor(2);
0296     errorGraph20t->Draw("p same");
0297     drawText("20 < p_{T} < 25 [GeV/c]", 0.5, 0.8, 12);
0298 
0299     csubjetR01->cd(7);
0300     gStyle->SetErrorX(0);
0301     gPad->SetTickx();
0302     gPad->SetTicky();
0303     // cjetSpectra->cd(5)->SetLogy();
0304     // hdummy_jetpT->Draw();
0305     _h_R04_theta_sj_25_30->SetMarkerStyle(29);
0306     _h_R04_theta_sj_25_30->SetMarkerColor(2);
0307     _h_R04_theta_sj_25_30->SetMarkerSize(1.0);
0308     _h_R04_theta_sj_25_30->SetTitle(" ");
0309     _h_R04_theta_sj_25_30->SetXTitle("#theta_{sj}");
0310     _h_R04_theta_sj_25_30->GetXaxis()->SetTitleSize(0.07);
0311     _h_R04_theta_sj_25_30->GetXaxis()->CenterTitle(true);
0312     _h_R04_theta_sj_25_30->GetXaxis()->SetLabelSize(0.04);
0313     _h_R04_theta_sj_25_30->SetAxisRange(0, 9.0, "Y");
0314     _h_R04_theta_sj_25_30->SetAxisRange(0, 0.5, "X");
0315     // _h_R04_theta_sj_25_30->Draw("p E1 same");
0316 
0317     _h_R04_theta_sj_25_30->Draw("hist C same");
0318     // Create a TGraphErrors for error bars
0319     int nbins25t = _h_R04_theta_sj_25_30->GetNbinsX();
0320     double* x25t = new double[nbins25t];
0321     double* y25t = new double[nbins25t];
0322     double* ex25t = new double[nbins25t];
0323     double* ey25t = new double[nbins25t];
0324     for (int i = 1; i <= nbins25t; ++i) {
0325       x25t[i - 1] = _h_R04_theta_sj_25_30->GetBinCenter(i);
0326       y25t[i - 1] = _h_R04_theta_sj_25_30->GetBinContent(i);
0327       ex25t[i - 1] = 0.0; // Assumes uniform bin width
0328       ey25t[i - 1] = _h_R04_theta_sj_25_30->GetBinError(i);
0329     }
0330     TGraphErrors* errorGraph25t = new TGraphErrors(nbins25t, x25t, y25t, ex25t, ey25t);
0331     errorGraph25t->SetMarkerStyle(29);
0332     errorGraph25t->SetMarkerColor(2);
0333     errorGraph25t->Draw("p same");
0334     drawText("25 < p_{T} < 30 [GeV/c]", 0.5, 0.8, 12);
0335 
0336     csubjetR01->cd(8);
0337     gStyle->SetErrorX(0);
0338     gPad->SetTickx();
0339     gPad->SetTicky();
0340     _h_R04_theta_sj_30_40->SetMarkerStyle(29);
0341     _h_R04_theta_sj_30_40->SetMarkerColor(2);
0342     _h_R04_theta_sj_30_40->SetMarkerSize(1.0);
0343     _h_R04_theta_sj_30_40->SetTitle(" ");
0344     _h_R04_theta_sj_30_40->SetXTitle("#theta_{sj}");
0345     _h_R04_theta_sj_30_40->GetXaxis()->SetTitleSize(0.07);
0346     _h_R04_theta_sj_30_40->GetXaxis()->CenterTitle(true);
0347     _h_R04_theta_sj_30_40->GetXaxis()->SetLabelSize(0.04);
0348     _h_R04_theta_sj_30_40->GetXaxis()->SetTitleSize(0.07);
0349     _h_R04_theta_sj_30_40->SetAxisRange(0, 9.0, "Y");
0350     _h_R04_theta_sj_30_40->SetAxisRange(0, 0.5, "X");
0351     //  _h_R04_theta_sj_30_40->Draw("p E1 same");
0352     _h_R04_theta_sj_30_40->Draw("hist C same");
0353     // Create a TGraphErrors for error bars
0354     int nbins30t = _h_R04_theta_sj_30_40->GetNbinsX();
0355     double* x30t = new double[nbins30t];
0356     double* y30t = new double[nbins30t];
0357     double* ex30t = new double[nbins30t];
0358     double* ey30t = new double[nbins30t];
0359     for (int i = 1; i <= nbins30t; ++i) {
0360       x30t[i - 1] = _h_R04_theta_sj_30_40->GetBinCenter(i);
0361       y30t[i - 1] = _h_R04_theta_sj_30_40->GetBinContent(i);
0362       ex30t[i - 1] = 0.0; // Assumes uniform bin width
0363       ey30t[i - 1] = _h_R04_theta_sj_30_40->GetBinError(i);
0364     }
0365     TGraphErrors* errorGraph30t = new TGraphErrors(nbins30t, x30t, y30t, ex30t, ey30t);
0366     errorGraph30t->SetMarkerStyle(29);
0367     errorGraph30t->SetMarkerColor(2);
0368     errorGraph30t->Draw("p same");
0369 
0370     drawText("30 < p_{T} < 40 [GeV/c]", 0.5, 0.8, 12);
0371   
0372     csubjetR01->SaveAs("RHIC_subjetR04_AKTB0z1.pdf", "RECREATE");
0373 
0374 
0375     TCanvas * csubjetR04 = new TCanvas("csubjetR04", " ", 1500, 500);
0376   
0377     makeMultiPanelCanvasWithGap(csubjetR04, 4, 2, 0.2, 0.2, 0.2, 0.2, 0.08, 0.01);
0378  
0379     csubjetR04->cd(1);
0380     gStyle->SetErrorX(0);
0381     gPad->SetTickx();
0382     gPad->SetTicky();
0383     _h_R04_z_g_15_20 ->SetLineStyle(1);
0384     _h_R04_z_g_15_20->SetLineWidth(1);
0385     _h_R04_z_g_15_20->SetTitle(" ");
0386     _h_R04_z_g_15_20->SetYTitle("(1/N_{jet}) dN/dz_{g}");
0387     _h_R04_z_g_15_20->GetYaxis()->SetTitleOffset(0.8);
0388     _h_R04_z_g_15_20->GetYaxis()->SetTitleSize(0.07);
0389     _h_R04_z_g_15_20->GetYaxis()->SetLabelSize(0.04);
0390     _h_R04_z_g_15_20->GetYaxis()->CenterTitle(true);
0391 
0392     _h_R04_z_g_15_20->SetXTitle("z_{g}");
0393     _h_R04_z_g_15_20->GetXaxis()->SetTitleSize(0.07);
0394     _h_R04_z_g_15_20->GetXaxis()->CenterTitle(true);
0395     _h_R04_z_g_15_20->GetXaxis()->SetLabelSize(0.04);
0396     _h_R04_z_g_15_20->SetMarkerSize(1.0);
0397     _h_R04_z_g_15_20->SetAxisRange(0, 6.5, "Y");
0398     _h_R04_z_g_15_20->SetAxisRange(0, 0.5, "X");
0399     //  _h_R04_z_g_15_20->Draw("p E2 same");
0400     _h_R04_z_g_15_20->Draw("hist C same");
0401     // Create a TGraphErrors for error bars
0402     int nbins15zg = _h_R04_z_g_15_20->GetNbinsX();
0403     double* x15zg = new double[nbins15zg];
0404     double* y15zg = new double[nbins15zg];
0405     double* ex15zg = new double[nbins15zg];
0406     double* ey15zg = new double[nbins15zg];
0407     for (int i = 1; i <= nbins15zg; ++i) {
0408       x15zg[i - 1] = _h_R04_z_g_15_20->GetBinCenter(i);
0409       y15zg[i - 1] = _h_R04_z_g_15_20->GetBinContent(i);
0410       ex15zg[i - 1] = 0.0; // Assumes uniform bin width
0411       ey15zg[i - 1] = _h_R04_z_g_15_20->GetBinError(i);
0412     }
0413     TGraphErrors* errorGraph15zg = new TGraphErrors(nbins15zg, x15zg, y15zg, ex15zg, ey15zg);
0414     errorGraph15zg->SetMarkerStyle(29);
0415     errorGraph15zg->SetMarkerColor(2);
0416     errorGraph15zg->Draw("p same");
0417   
0418 
0419     csubjetR04->cd(2);
0420     gStyle->SetErrorX(0);
0421     gPad->SetTickx();
0422     gPad->SetTicky();
0423     // cjetSpectra->cd(1)->SetLogy();
0424     // hdummy_jetpT->Draw();
0425     _h_R04_z_g_20_25->SetTitle(" ");
0426     _h_R04_z_g_20_25->SetXTitle("z_{g}");
0427     _h_R04_z_g_20_25->GetXaxis()->SetTitleSize(0.07);
0428     _h_R04_z_g_20_25->GetXaxis()->CenterTitle(true);
0429     _h_R04_z_g_20_25->GetXaxis()->SetLabelSize(0.04);
0430     _h_R04_z_g_20_25->SetAxisRange(0, 6.5, "Y");
0431     _h_R04_z_g_20_25->SetAxisRange(0, 0.5, "X");
0432     _h_R04_z_g_20_25->SetMarkerSize(1.0);
0433     // _h_R04_z_g_20_25->Draw("p E1 same");
0434     _h_R04_z_g_20_25->Draw("hist C same");
0435 
0436     // Create a TGraphErrors for error bars
0437     int nbins20zg = _h_R04_z_g_20_25->GetNbinsX();
0438     double* x20zg = new double[nbins20zg];
0439     double* y20zg = new double[nbins20zg];
0440     double* ex20zg = new double[nbins20zg];
0441     double* ey20zg = new double[nbins20zg];
0442     for (int i = 1; i <= nbins20; ++i) {
0443       x20zg[i - 1] = _h_R04_z_g_20_25->GetBinCenter(i);
0444       y20zg[i - 1] = _h_R04_z_g_20_25->GetBinContent(i);
0445       ex20zg[i - 1] = 0.0; // Assumes uniform bin width
0446       ey20zg[i - 1] = _h_R04_z_g_20_25->GetBinError(i);
0447     }
0448     TGraphErrors* errorGraph20zg = new TGraphErrors(nbins20, x20zg, y20zg, ex20zg, ey20zg);
0449     errorGraph20zg->SetMarkerStyle(29);
0450     errorGraph20zg->SetMarkerColor(2);
0451     errorGraph20zg->Draw("p same");
0452 
0453     csubjetR04->cd(3);
0454     gStyle->SetErrorX(0);
0455     gPad->SetTickx();
0456     gPad->SetTicky();
0457     // cjetSpectra->cd(1)->SetLogy();
0458     // hdummy_jetpT->Draw();
0459     _h_R04_z_g_25_30->SetMarkerStyle(29);
0460     _h_R04_z_g_25_30->SetMarkerColor(2);
0461     _h_R04_z_g_25_30->SetTitle(" ");
0462     _h_R04_z_g_25_30->SetXTitle("z_{g}");
0463     _h_R04_z_g_25_30->GetXaxis()->SetTitleSize(0.07);
0464     _h_R04_z_g_25_30->GetXaxis()->CenterTitle(true);
0465     _h_R04_z_g_25_30->GetXaxis()->SetLabelSize(0.04);
0466     _h_R04_z_g_25_30->SetAxisRange(0, 6.5, "Y");
0467     _h_R04_z_g_25_30->SetAxisRange(0, 0.5, "X");
0468     _h_R04_z_g_25_30->SetMarkerSize(1.0);
0469     //  _h_R04_z_g_25_30->Draw("p E1 same");
0470     _h_R04_z_g_25_30->Draw("hist C same");
0471     // Create a TGraphErrors for error bars
0472     int nbins25zg = _h_R04_z_g_25_30->GetNbinsX();
0473     double* x25zg = new double[nbins25zg];
0474     double* y25zg = new double[nbins25zg];
0475     double* ex25zg = new double[nbins25zg];
0476     double* ey25zg = new double[nbins25zg];
0477     for (int i = 1; i <= nbins25zg; ++i) {
0478       x25zg[i - 1] = _h_R04_z_g_25_30->GetBinCenter(i);
0479       y25zg[i - 1] = _h_R04_z_g_25_30->GetBinContent(i);
0480       ex25zg[i - 1] = 0.0; // Assumes uniform bin width
0481       ey25zg[i - 1] = _h_R04_z_g_25_30->GetBinError(i);
0482     }
0483     TGraphErrors* errorGraph25zg = new TGraphErrors(nbins25zg, x25zg, y25zg, ex25zg, ey25zg);
0484     errorGraph25zg->SetMarkerStyle(29);
0485     errorGraph25zg->SetMarkerColor(2);
0486     errorGraph25zg->Draw("p same");
0487 
0488     csubjetR04->cd(4);
0489     gStyle->SetErrorX(0);
0490     gPad->SetTickx();
0491     gPad->SetTicky();
0492     _h_R04_z_g_30_40->SetMarkerStyle(29);
0493     _h_R04_z_g_30_40->SetMarkerColor(2);
0494     _h_R04_z_g_30_40->SetTitle(" ");
0495     _h_R04_z_g_30_40->SetXTitle("z_{g}");
0496     _h_R04_z_g_30_40->GetXaxis()->SetTitleSize(0.07);
0497     _h_R04_z_g_30_40->GetXaxis()->CenterTitle(true);
0498     _h_R04_z_g_30_40->GetXaxis()->SetLabelSize(0.04);
0499     _h_R04_z_g_30_40->GetXaxis()->SetTitleSize(0.07);
0500     _h_R04_z_g_30_40->SetAxisRange(0, 6.5, "Y");
0501     _h_R04_z_g_30_40->SetAxisRange(0, 0.5, "X");
0502     _h_R04_z_g_30_40->SetMarkerSize(1.0);
0503     //  _h_R04_z_g_30_40->Draw("p E1 same");
0504     _h_R04_z_g_30_40->Draw("hist C same");
0505   
0506     drawText("p+p #sqrt{s} = 200 GeV", 0.4, 0.85, 12);
0507     drawText("anti-k_{T} R=0.4 jets, |#eta_{jet}| < 0.6", 0.4, 0.75, 12);
0508     drawText("#beta=0.0, z_{cut}=0.1", 0.4, 0.65, 12);
0509    
0510     // Create a TGraphErrors for error bars
0511     int nbins30zg = _h_R04_z_g_30_40->GetNbinsX();
0512     double* x30zg = new double[nbins30zg];
0513     double* y30zg = new double[nbins30zg];
0514     double* ex30zg = new double[nbins30zg];
0515     double* ey30zg = new double[nbins30zg];
0516     for (int i = 1; i <= nbins20zg; ++i) {
0517       x30zg[i - 1] = _h_R04_z_g_30_40->GetBinCenter(i);
0518       y30zg[i - 1] = _h_R04_z_g_30_40->GetBinContent(i);
0519       ex30zg[i - 1] = 0.0; // Assumes uniform bin width
0520       ey30zg[i - 1] = _h_R04_z_g_30_40->GetBinError(i);
0521     }
0522     TGraphErrors* errorGraph30zg = new TGraphErrors(nbins15zg, x30zg, y30zg, ex30zg, ey30zg);
0523     errorGraph30zg->SetMarkerStyle(29);
0524     errorGraph30zg->SetMarkerColor(2);
0525     errorGraph30zg->Draw("p same");
0526 
0527     csubjetR04->cd(5);
0528     gStyle->SetErrorX(0);
0529     gPad->SetTickx();
0530     gPad->SetTicky();
0531     _h_R04_theta_g_15_20->SetMarkerStyle(29);
0532     _h_R04_theta_g_15_20->SetMarkerColor(2);
0533     _h_R04_theta_g_15_20->SetMarkerSize(1.0);
0534     _h_R04_theta_g_15_20->SetTitle(" ");
0535     _h_R04_theta_g_15_20->SetYTitle("(1/N_{jet}) dN/dR_{g}");
0536     _h_R04_theta_g_15_20->SetAxisRange(0, 9.0, "Y");
0537     _h_R04_theta_g_15_20->SetAxisRange(0, 0.5, "X");
0538     _h_R04_theta_g_15_20->GetYaxis()->SetTitleOffset(0.8);
0539     _h_R04_theta_g_15_20->GetYaxis()->SetTitleSize(0.07);
0540     _h_R04_theta_g_15_20->GetYaxis()->SetLabelSize(0.04);
0541     _h_R04_theta_g_15_20->GetYaxis()->CenterTitle(true);
0542 
0543     _h_R04_theta_g_15_20->SetXTitle("R_{g}");
0544     _h_R04_theta_g_15_20->GetXaxis()->SetTitleSize(0.07);
0545     _h_R04_theta_g_15_20->GetXaxis()->CenterTitle(true);
0546     _h_R04_theta_g_15_20->GetXaxis()->SetLabelSize(0.04);
0547 
0548     _h_R04_theta_g_15_20->Draw("hist C same");
0549     // Create a TGraphErrors for error bars
0550     int nbins15tgt = _h_R04_theta_g_15_20->GetNbinsX();
0551     double* x15gt = new double[nbins15tgt];
0552     double* y15gt = new double[nbins15tgt];
0553     double* ex15gt = new double[nbins15tgt];
0554     double* ey15gt = new double[nbins15tgt];
0555     for (int i = 1; i <= nbins15tgt; ++i) {
0556       x15gt[i - 1] = _h_R04_theta_g_15_20->GetBinCenter(i);
0557       y15gt[i - 1] = _h_R04_theta_g_15_20->GetBinContent(i);
0558       ex15gt[i - 1] = 0.0; // Assumes uniform bin width
0559       ey15gt[i - 1] = _h_R04_theta_g_15_20->GetBinError(i);
0560     }
0561     TGraphErrors* errorGraph15gt = new TGraphErrors(nbins15tgt, x15gt, y15gt, ex15gt, ey15gt);
0562     errorGraph15gt->SetMarkerStyle(29);
0563     errorGraph15gt->SetMarkerColor(2);
0564     errorGraph15gt->Draw("p same");
0565 
0566     drawText("15 < p_{T} < 20 GeV/c", 0.7, 0.8, 12);
0567     
0568     csubjetR04->cd(6);
0569     gStyle->SetErrorX(0);
0570     gPad->SetTickx();
0571     gPad->SetTicky();
0572     // cjetSpectra->cd(5)->SetLogy();
0573     // hdummy_jetpT->Draw();
0574     _h_R04_theta_g_20_25->SetMarkerStyle(29);
0575     _h_R04_theta_g_20_25->SetMarkerColor(2);
0576     _h_R04_theta_g_20_25->SetMarkerSize(1.0);
0577     _h_R04_theta_g_20_25->SetTitle(" ");
0578     _h_R04_theta_g_20_25->SetXTitle("R_{g}");
0579     _h_R04_theta_g_20_25->GetXaxis()->SetTitleSize(0.07);
0580     _h_R04_theta_g_20_25->GetXaxis()->CenterTitle(true);
0581     _h_R04_theta_g_20_25->GetXaxis()->SetLabelSize(0.04);
0582     _h_R04_theta_g_20_25->SetAxisRange(0, 9.0, "Y");
0583     _h_R04_theta_g_20_25->SetAxisRange(0, 0.5, "X");
0584     _h_R04_theta_g_20_25->Draw("hist C same");
0585     //  _h_R04_theta_g_20_25->Draw("p E1 same");
0586     // Create a TGraphErrors for error bars
0587     int nbins20tgt = _h_R04_theta_g_20_25->GetNbinsX();
0588     double* x20gt = new double[nbins20tgt];
0589     double* y20gt = new double[nbins20tgt];
0590     double* ex20gt = new double[nbins20tgt];
0591     double* ey20gt = new double[nbins20tgt];
0592     for (int i = 1; i <= nbins15tgt; ++i) {
0593       x20gt[i - 1] = _h_R04_theta_g_20_25->GetBinCenter(i);
0594       y20gt[i - 1] = _h_R04_theta_g_20_25->GetBinContent(i);
0595       ex20gt[i - 1] =0.0; // Assumes uniform bin width
0596       ey20gt[i - 1] = _h_R04_theta_g_20_25->GetBinError(i);
0597     }
0598     TGraphErrors* errorGraph20gt = new TGraphErrors(nbins20tgt, x20gt, y20gt, ex20gt, ey20gt);
0599     errorGraph20gt->SetMarkerStyle(29);
0600     errorGraph20gt->SetMarkerColor(2);
0601     errorGraph20gt->Draw("p same");
0602     drawText("20 < p_{T} < 25 GeV/c", 0.6, 0.8, 12);
0603 
0604     csubjetR04->cd(7);
0605     gStyle->SetErrorX(0);
0606     gPad->SetTickx();
0607     gPad->SetTicky();
0608     // cjetSpectra->cd(5)->SetLogy();
0609     // hdummy_jetpT->Draw();
0610     _h_R04_theta_g_25_30->SetMarkerStyle(29);
0611     _h_R04_theta_g_25_30->SetMarkerColor(2);
0612     _h_R04_theta_g_25_30->SetMarkerSize(1.0);
0613     _h_R04_theta_g_25_30->SetTitle(" ");
0614     _h_R04_theta_g_25_30->SetXTitle("R_{g}");
0615     _h_R04_theta_g_25_30->GetXaxis()->SetTitleSize(0.07);
0616     _h_R04_theta_g_25_30->GetXaxis()->CenterTitle(true);
0617     _h_R04_theta_g_25_30->GetXaxis()->SetLabelSize(0.04);
0618     _h_R04_theta_g_25_30->SetAxisRange(0, 9.0, "Y");
0619     _h_R04_theta_g_25_30->SetAxisRange(0, 0.5, "X");
0620     // _h_R04_theta_g_25_30->Draw("p E1 same");
0621 
0622     _h_R04_theta_g_25_30->Draw("hist C same");
0623     // Create a TGraphErrors for error bars
0624     int nbins25tgt = _h_R04_theta_g_25_30->GetNbinsX();
0625     double* x25gt = new double[nbins25tgt];
0626     double* y25gt = new double[nbins25tgt];
0627     double* ex25gt = new double[nbins25tgt];
0628     double* ey25gt = new double[nbins25tgt];
0629     for (int i = 1; i <= nbins25tgt; ++i) {
0630       x25gt[i - 1] = _h_R04_theta_g_25_30->GetBinCenter(i);
0631       y25gt[i - 1] = _h_R04_theta_g_25_30->GetBinContent(i);
0632       ex25gt[i - 1] = 0.0;
0633       //ex25t[i - 1] = _h_R04_theta_g_25_30->GetBinWidth(i) / 2.0; // Assumes uniform bin width
0634       ey25gt[i - 1] = _h_R04_theta_g_25_30->GetBinError(i);
0635     }
0636     TGraphErrors* errorGraph25gt = new TGraphErrors(nbins25tgt, x25gt, y25gt, ex25gt, ey25gt);
0637     errorGraph25gt->SetMarkerStyle(29);
0638     errorGraph25gt->SetMarkerColor(2);
0639     errorGraph25gt->Draw("p same");
0640     drawText("25 < p_{T} < 30 GeV/c", 0.6, 0.8, 12);
0641 
0642     csubjetR04->cd(8);
0643     gStyle->SetErrorX(0);
0644     gPad->SetTickx();
0645     gPad->SetTicky();
0646     _h_R04_theta_g_30_40->SetMarkerStyle(29);
0647     _h_R04_theta_g_30_40->SetMarkerColor(2);
0648     _h_R04_theta_g_30_40->SetMarkerSize(1.0);
0649     _h_R04_theta_g_30_40->SetTitle(" ");
0650     _h_R04_theta_g_30_40->SetXTitle("R_{g}");
0651     _h_R04_theta_g_30_40->GetXaxis()->SetTitleSize(0.07);
0652     _h_R04_theta_g_30_40->GetXaxis()->CenterTitle(true);
0653     _h_R04_theta_g_30_40->GetXaxis()->SetLabelSize(0.04);
0654     _h_R04_theta_g_30_40->GetXaxis()->SetTitleSize(0.07);
0655     _h_R04_theta_g_30_40->SetAxisRange(0, 9.0, "Y");
0656     _h_R04_theta_g_30_40->SetAxisRange(0, 0.5, "X");
0657     //  _h_R04_theta_g_30_40->Draw("p E1 same");
0658     _h_R04_theta_g_30_40->Draw("hist C same");
0659     // Create a TGraphErrors for error bars
0660     int nbins30tgt = _h_R04_theta_g_30_40->GetNbinsX();
0661     double* x30gt = new double[nbins30tgt];
0662     double* y30gt = new double[nbins30tgt];
0663     double* ex30gt = new double[nbins30tgt];
0664     double* ey30gt = new double[nbins30tgt];
0665     for (int i = 1; i <= nbins20tgt; ++i) {
0666       x30gt[i - 1] = _h_R04_theta_g_30_40->GetBinCenter(i);
0667       y30gt[i - 1] = _h_R04_theta_g_30_40->GetBinContent(i);
0668       // ex30tt[i - 1] = _h_R04_theta_g_30_40->GetBinWidth(i) / 2.0; // Assumes uniform bin width
0669       ex30gt[i - 1] = 0.0;
0670       ey30gt[i - 1] = _h_R04_theta_g_30_40->GetBinError(i);
0671     }
0672     TGraphErrors* errorGraph30gt = new TGraphErrors(nbins20tgt, x30gt, y30gt, ex30gt, ey30gt);
0673     errorGraph30gt->SetMarkerStyle(29);
0674     errorGraph30gt->SetMarkerColor(2);
0675     errorGraph30gt->Draw("p same");
0676 
0677     drawText("30 < p_{T} < 40 GeV/c", 0.55, 0.8, 12);
0678   
0679     csubjetR04->SaveAs("RHIC_subjetR04_CAB0z1.pdf", "RECREATE");
0680 
0681 
0682     TCanvas * canvas = new TCanvas("canvas", " ", 1500, 500);
0683    
0684     makeMultiPanelCanvasWithGap(canvas, 4, 2, 0.2, 0.2, 0.2, 0.2, 0.09, 0.01);
0685 
0686     canvas->cd(1);
0687     gStyle->SetErrorX(0);
0688     gPad->SetTickx();
0689     gPad->SetTicky();
0690    
0691     correlation_z_15_20->SetMarkerStyle(29);
0692     correlation_z_15_20->SetMarkerColor(2);
0693     correlation_z_15_20->SetMarkerSize(1.0);
0694     correlation_z_15_20->GetXaxis()->SetTitleSize(0.07);
0695     correlation_z_15_20->GetYaxis()->SetTitleSize(0.07);
0696     correlation_z_15_20->GetXaxis()->CenterTitle(true);
0697     correlation_z_15_20->GetXaxis()->SetLabelSize(0.04);
0698     correlation_z_15_20->GetXaxis()->SetTitleSize(0.07);
0699     correlation_z_15_20->SetAxisRange(0, 0.6, "Y");
0700     correlation_z_15_20->SetAxisRange(0, 0.5, "X");
0701     correlation_z_15_20->SetXTitle("z_{g}");
0702     correlation_z_15_20->SetYTitle("z_{sj}");
0703     correlation_z_15_20->GetYaxis()->SetLabelSize(0.04);
0704     correlation_z_15_20->SetTitle(" ");
0705     correlation_z_15_20->Draw("colz"); // Use "colz" option to draw the color map
0706     drawText("15 < p_{T} < 20 GeV/c", 0.4, 0.75, 12);
0707 
0708    
0709     canvas->cd(2);
0710     gStyle->SetErrorX(0);
0711     gPad->SetTickx();
0712     gPad->SetTicky();
0713    
0714     correlation_z_20_25->SetMarkerStyle(29);
0715     correlation_z_20_25->SetMarkerColor(2);
0716     correlation_z_20_25->SetMarkerSize(1.0);
0717     correlation_z_20_25->GetXaxis()->SetTitleSize(0.07);
0718     correlation_z_20_25->GetXaxis()->CenterTitle(true);
0719     correlation_z_20_25->GetXaxis()->SetLabelSize(0.04);
0720     correlation_z_20_25->GetXaxis()->SetTitleSize(0.07);
0721     correlation_z_20_25->SetAxisRange(0, 0.6, "Y");
0722     correlation_z_20_25->SetAxisRange(0, 0.5, "X");
0723     correlation_z_20_25->SetXTitle("z_{g}");
0724     // correlation_z_20_25->SetYTitle("z_{sj}");
0725     correlation_z_20_25->SetTitle(" ");
0726     correlation_z_20_25->Draw("colz"); // Use "colz" option to draw the color map
0727     drawText("20 < p_{T} < 25 GeV/c", 0.3, 0.75, 12);
0728 
0729    
0730     canvas->cd(3);
0731     gStyle->SetErrorX(0);
0732     gPad->SetTickx();
0733     gPad->SetTicky();
0734    
0735     correlation_z_25_30->SetMarkerStyle(29);
0736     correlation_z_25_30->SetMarkerColor(2);
0737     correlation_z_25_30->SetMarkerSize(1.0);
0738     correlation_z_25_30->GetXaxis()->SetTitleSize(0.07);
0739     correlation_z_25_30->GetXaxis()->CenterTitle(true);
0740     correlation_z_25_30->GetXaxis()->SetLabelSize(0.04);
0741     correlation_z_25_30->GetXaxis()->SetTitleSize(0.07);
0742     correlation_z_25_30->SetAxisRange(0, 0.6, "Y");
0743     correlation_z_25_30->SetAxisRange(0, 0.5, "X");
0744     correlation_z_25_30->SetXTitle("z_{g}");
0745     // correlation_z_25_30->SetYTitle("z_{sj}");
0746     correlation_z_25_30->SetTitle(" ");
0747     correlation_z_25_30->Draw("colz"); // Use "colz" option to draw the color map
0748 
0749     drawText("25 < p_{T} < 30 GeV/c", 0.3, 0.75, 12);
0750     canvas->cd(4);
0751     gStyle->SetErrorX(0);
0752     gPad->SetTickx();
0753     gPad->SetTicky();
0754    
0755     correlation_z_30_40->SetMarkerStyle(29);
0756     correlation_z_30_40->SetMarkerColor(2);
0757     correlation_z_30_40->SetMarkerSize(1.0);
0758     correlation_z_30_40->GetXaxis()->SetTitleSize(0.07);
0759     correlation_z_30_40->GetXaxis()->CenterTitle(true);
0760     correlation_z_30_40->GetXaxis()->SetLabelSize(0.04);
0761     correlation_z_30_40->GetXaxis()->SetTitleSize(0.07);
0762     correlation_z_30_40->SetAxisRange(0, 0.6, "Y");
0763     correlation_z_30_40->SetAxisRange(0, 0.5, "X");
0764     correlation_z_30_40->SetXTitle("z_{g}");
0765     // correlation_z_30_40->SetYTitle("z_{sj}");
0766     correlation_z_30_40->SetTitle(" ");
0767     correlation_z_30_40->Draw("colz"); // Use "colz" option to draw the color map
0768    
0769     drawText("30 < p_{T} < 40 GeV/c", 0.3, 0.75, 12);
0770 
0771     canvas->cd(5);
0772     gStyle->SetErrorX(0);
0773     gPad->SetTickx();
0774     gPad->SetTicky();
0775    
0776    
0777     correlation_theta_15_20->SetMarkerStyle(29);
0778     correlation_theta_15_20->SetMarkerColor(2);
0779     correlation_theta_15_20->SetMarkerSize(1.0);
0780     correlation_theta_15_20->SetTitle(" ");
0781     correlation_theta_15_20->SetXTitle("#theta_{g}");
0782     correlation_theta_15_20->GetXaxis()->SetTitleSize(0.07);
0783     correlation_theta_15_20->GetXaxis()->CenterTitle(true);
0784     correlation_theta_15_20->GetXaxis()->SetLabelSize(0.04);
0785     correlation_theta_15_20->GetYaxis()->SetTitleSize(0.07);
0786     correlation_theta_15_20->SetAxisRange(0, 0.5, "Y");
0787     correlation_theta_15_20->SetAxisRange(0, 0.5, "X");
0788     correlation_theta_15_20->SetXTitle("R_{g}");
0789     correlation_theta_15_20->SetYTitle("#theta_{sj}");
0790     correlation_theta_15_20->GetYaxis()->SetLabelSize(0.04);
0791     correlation_theta_15_20->SetTitle(" ");
0792     correlation_theta_15_20->Draw("colz"); // Use "colz" option to draw the color map
0793    
0794     canvas->cd(6);
0795     gStyle->SetErrorX(0);
0796     gPad->SetTickx();
0797     gPad->SetTicky();
0798 
0799     correlation_theta_20_25->SetMarkerStyle(29);
0800     correlation_theta_20_25->SetMarkerColor(2);
0801     correlation_theta_20_25->SetMarkerSize(1.0);
0802     correlation_theta_20_25->SetTitle(" ");
0803     correlation_theta_20_25->SetXTitle("#theta_{g}");
0804     correlation_theta_20_25->GetXaxis()->SetTitleSize(0.07);
0805     correlation_theta_20_25->GetXaxis()->CenterTitle(true);
0806     correlation_theta_20_25->GetXaxis()->SetLabelSize(0.04);
0807     correlation_theta_20_25->GetXaxis()->SetTitleSize(0.07);
0808     correlation_theta_20_25->SetAxisRange(0, 0.5, "Y");
0809     correlation_theta_20_25->SetAxisRange(0, 0.5, "X");
0810     correlation_theta_20_25->SetXTitle("R_{g}");
0811     // correlation_theta_20_25->SetYTitle("#theta_{sj}");
0812     correlation_theta_20_25->SetTitle(" ");
0813     correlation_theta_20_25->Draw("colz"); // Use "colz" option to draw the color map
0814     drawText("p+p #sqrt{s} = 200 GeV", 0.05, 0.9, 12);
0815    
0816     canvas->cd(7);
0817     gStyle->SetErrorX(0);
0818     gPad->SetTickx();
0819     gPad->SetTicky();
0820    
0821     correlation_theta_25_30->SetMarkerStyle(29);
0822     correlation_theta_25_30->SetMarkerColor(2);
0823     correlation_theta_25_30->SetMarkerSize(1.0);
0824     correlation_theta_25_30->SetTitle(" ");
0825     correlation_theta_25_30->SetXTitle("#theta_{g}");
0826     correlation_theta_25_30->GetXaxis()->SetTitleSize(0.07);
0827     correlation_theta_25_30->GetXaxis()->CenterTitle(true);
0828     correlation_theta_25_30->GetXaxis()->SetLabelSize(0.04);
0829     correlation_theta_25_30->GetXaxis()->SetTitleSize(0.07);
0830     correlation_theta_25_30->SetAxisRange(0, 0.5, "Y");
0831     correlation_theta_25_30->SetAxisRange(0, 0.5, "X");
0832     correlation_theta_25_30->SetXTitle("R_{g}");
0833     // correlation_theta_25_30->SetYTitle("#theta_{sj}");
0834     correlation_theta_25_30->SetTitle(" ");
0835     correlation_theta_25_30->Draw("colz"); // Use "colz" option to draw the color map
0836     drawText("anti-k_{T} R=0.4 jets, |#eta_{jet}| < 0.6", 0.03, 0.9, 12);
0837 
0838     canvas->cd(8);
0839     gStyle->SetErrorX(0);
0840     gPad->SetTickx();
0841     gPad->SetTicky();
0842    
0843     correlation_theta_30_40->SetMarkerStyle(29);
0844     correlation_theta_30_40->SetMarkerColor(2);
0845     correlation_theta_30_40->SetMarkerSize(1.0);
0846     correlation_theta_30_40->SetTitle(" ");
0847     correlation_theta_30_40->GetXaxis()->SetTitleSize(0.07);
0848     correlation_theta_30_40->GetXaxis()->CenterTitle(true);
0849     correlation_theta_30_40->GetXaxis()->SetLabelSize(0.04);
0850     correlation_theta_30_40->SetAxisRange(0, 0.5, "Y");
0851     correlation_theta_30_40->SetAxisRange(0, 0.5, "X");
0852     correlation_theta_30_40->SetXTitle("R_{g}");
0853     // correlation_theta_30_40->SetYTitle("#theta_{sj}");
0854     correlation_theta_30_40->SetTitle(" ");
0855     correlation_theta_30_40->Draw("colz"); // Use "colz" option to draw the color map
0856     drawText("#beta=0.0, z_{cut}=0.1", 0.05, 0.9, 12);
0857    
0858    
0859     canvas->SaveAs("CorellBeta0z1.pdf", "RECREATE");
0860    
0861   }
0862 
0863 }
0864
[ Source navigation ] [ Diff markup ] [ Identifier search ] [ general search ]

This page was automatically generated by the 2.3.7 LXR engine.

The LXR team
Valid CSS 2.1! Valid HTML 4.01!