File indexing completed on 2025-08-05 08:13:19
0001 void SetHist(TH2* h, string xt ="", string yt ="",int color = 4, int marker = 20,int width = 3, float size = 1.0);
0002 void SetHist(TH2* h, int color = 4);
0003 void SetLeg(TLegend* l,float txtsize=0.03);
0004
0005
0006 void DrawUECentralityPlots(string infile = "UEvsEtaCentrality.root"){
0007
0008 gStyle->SetOptStat(0);
0009
0010 TFile *_file0 = TFile::Open(infile.c_str());
0011
0012 TH2* hv2_cent = (TH2*)gROOT->FindObject("hv2_cent");
0013 TH2* hPsi2_cent = (TH2*)gROOT->FindObject("hPsi2_cent");
0014 TCanvas* can1 = new TCanvas("can1","",800,400);
0015 can1->Divide(2,1,0,0);
0016
0017 can1->cd(1);
0018 gPad->SetLeftMargin(0.15);
0019 gPad->SetRightMargin(.15);
0020 gPad->SetTopMargin(.07);
0021 gPad->SetBottomMargin(.10);
0022 SetHist(hv2_cent,"Centrality %","v2 value",1);
0023 hv2_cent->SetTitle("v2 vs Centrality");
0024 hv2_cent->Draw("ColZ");
0025
0026 TProfile* v2profile = hv2_cent->ProfileX("Mean Values of v2");
0027
0028 gStyle->SetTitleTextColor(1);
0029 v2profile->Draw("same");
0030
0031 can1->cd(2);
0032 gPad->SetLeftMargin(0.15);
0033 gPad->SetRightMargin(.15);
0034 gPad->SetTopMargin(.07);
0035 gPad->SetBottomMargin(.10);
0036 SetHist(hPsi2_cent,"Centrality %","Psi2 value",1);
0037 hPsi2_cent->SetTitle("Psi2 vs Centrality");
0038 hPsi2_cent->Draw("ColZ");
0039
0040 TProfile* Psi2profile = hPsi2_cent->ProfileX("Mean Values of Psi2");
0041
0042 gStyle->SetTitleTextColor(1);
0043 Psi2profile->Draw("same");
0044
0045 TH2* hUEiHcalEta_Cent0_20 = (TH2*)gROOT->FindObject("hUEiHcalEta_Cent0_20");
0046 TH2* hUEoHcalEta_Cent0_20 = (TH2*)gROOT->FindObject("hUEoHcalEta_Cent0_20");
0047 TH2* hUEemcalEta_Cent0_20 = (TH2*)gROOT->FindObject("hUEemcalEta_Cent0_20");
0048
0049 TH2* hUEiHcalEta_Cent20_50 = (TH2*)gROOT->FindObject("hUEiHcalEta_Cent20_50");
0050 TH2* hUEoHcalEta_Cent20_50 = (TH2*)gROOT->FindObject("hUEoHcalEta_Cent20_50");
0051 TH2* hUEemcalEta_Cent20_50 = (TH2*)gROOT->FindObject("hUEemcalEta_Cent20_50");
0052
0053 TH2* hUEiHcalEta_Cent50_100 = (TH2*)gROOT->FindObject("hUEiHcalEta_Cent50_100");
0054 TH2* hUEoHcalEta_Cent50_100 = (TH2*)gROOT->FindObject("hUEoHcalEta_Cent50_100");
0055 TH2* hUEemcalEta_Cent50_100 = (TH2*)gROOT->FindObject("hUEemcalEta_Cent50_100");
0056
0057 TCanvas* can2 = new TCanvas("can2","",1200,400);
0058 can2->Divide(3,1,0,0);
0059
0060 can2->cd(1);
0061 gPad->SetLeftMargin(0.15);
0062 gPad->SetRightMargin(.10);
0063 gPad->SetTopMargin(.07);
0064 gPad->SetBottomMargin(.10);
0065 SetHist(hUEiHcalEta_Cent0_20,"Eta","UE Size/Tower (GeV)",1);
0066 hUEiHcalEta_Cent0_20->SetTitle("UE Size vs Eta 0 - 20 % Centrality Range iHCal");
0067 hUEiHcalEta_Cent0_20->Draw("ColZ");
0068 hUEiHcalEta_Cent0_20->ProfileX()->Draw("same");
0069
0070 can2->cd(2);
0071 gPad->SetLeftMargin(0.15);
0072 gPad->SetRightMargin(.15);
0073 gPad->SetTopMargin(.07);
0074 gPad->SetBottomMargin(.10);
0075 SetHist(hUEiHcalEta_Cent20_50,"Eta","UE Size/Tower (GeV)",1);
0076 hUEiHcalEta_Cent20_50->SetTitle("UE Size vs Eta 20 - 50 % Centrality Range iHCal");
0077 hUEiHcalEta_Cent20_50->GetYaxis()->SetRangeUser(0.0,0.05);
0078 hUEiHcalEta_Cent20_50->Draw("ColZ");
0079 hUEiHcalEta_Cent20_50->ProfileX()->Draw("same");
0080
0081 can2->cd(3);
0082 gPad->SetLeftMargin(0.15);
0083 gPad->SetRightMargin(.10);
0084 gPad->SetTopMargin(.07);
0085 gPad->SetBottomMargin(.10);
0086 SetHist(hUEiHcalEta_Cent50_100,"Eta","UE Size/Tower (GeV)",1);
0087 hUEiHcalEta_Cent50_100->SetTitle("UE Size vs Eta 50 - 100 % Centrality Range iHCal");
0088 hUEiHcalEta_Cent50_100->GetYaxis()->SetRangeUser(0.0,0.03);
0089 hUEiHcalEta_Cent50_100->Draw("ColZ");
0090 hUEiHcalEta_Cent50_100->ProfileX()->Draw("same");
0091
0092 TCanvas* can3 = new TCanvas("can3","",1200,400);
0093 can3->Divide(3,1,0,0);
0094
0095 can3->cd(1);
0096 gPad->SetLeftMargin(0.15);
0097 gPad->SetRightMargin(.10);
0098 gPad->SetTopMargin(.07);
0099 gPad->SetBottomMargin(.10);
0100 SetHist(hUEoHcalEta_Cent0_20,"Eta","UE Size/Tower (GeV)",1);
0101 hUEoHcalEta_Cent0_20->SetTitle("UE Size vs Eta 0 - 20 % Centrality Range oHCal");
0102 hUEoHcalEta_Cent0_20->Draw("ColZ");
0103 hUEoHcalEta_Cent0_20->ProfileX()->Draw("same");
0104
0105 can3->cd(2);
0106 gPad->SetLeftMargin(0.15);
0107 gPad->SetRightMargin(.15);
0108 gPad->SetTopMargin(.07);
0109 gPad->SetBottomMargin(.10);
0110 SetHist(hUEoHcalEta_Cent20_50,"Eta","UE Size/Tower (GeV)",1);
0111 hUEoHcalEta_Cent20_50->SetTitle("UE Size vs Eta 20 - 50 % Centrality Range oHCal");
0112 hUEoHcalEta_Cent20_50->GetYaxis()->SetRangeUser(0.0,0.25);
0113 hUEoHcalEta_Cent20_50->Draw("ColZ");
0114 hUEoHcalEta_Cent20_50->ProfileX()->Draw("same");
0115
0116 can3->cd(3);
0117 gPad->SetLeftMargin(0.15);
0118 gPad->SetRightMargin(.10);
0119 gPad->SetTopMargin(.07);
0120 gPad->SetBottomMargin(.10);
0121 SetHist(hUEoHcalEta_Cent50_100,"Eta","UE Size/Tower (GeV)",1);
0122 hUEoHcalEta_Cent50_100->SetTitle("UE Size vs Eta 50 - 100 % Centrality Range oHCal");
0123 hUEoHcalEta_Cent50_100->GetYaxis()->SetRangeUser(0.0,0.15);
0124 hUEoHcalEta_Cent50_100->Draw("ColZ");
0125 hUEoHcalEta_Cent50_100->ProfileX()->Draw("same");
0126
0127 TCanvas* can4 = new TCanvas("can4","",1200,400);
0128 can4->Divide(3,1,0,0);
0129
0130
0131 can4->cd(1);
0132 gPad->SetLeftMargin(0.15);
0133 gPad->SetRightMargin(.10);
0134 gPad->SetTopMargin(.07);
0135 gPad->SetBottomMargin(.10);
0136 SetHist(hUEemcalEta_Cent0_20,"Eta","UE Size/Tower (GeV)",1);
0137 hUEemcalEta_Cent0_20->SetTitle("UE Size vs Eta 0 - 20 % Centrality Range EMCal");
0138 hUEemcalEta_Cent0_20->Draw("ColZ");
0139 hUEemcalEta_Cent0_20->ProfileX()->Draw("same");
0140
0141 can4->cd(2);
0142 gPad->SetLeftMargin(0.15);
0143 gPad->SetRightMargin(.15);
0144 gPad->SetTopMargin(.07);
0145 gPad->SetBottomMargin(.10);
0146 SetHist(hUEemcalEta_Cent20_50,"Eta","UE Size/Tower (GeV)",1);
0147 hUEemcalEta_Cent20_50->SetTitle("UE Size vs Eta 20 - 50 % Centrality Range EMCal");
0148 hUEemcalEta_Cent20_50->GetYaxis()->SetRangeUser(0.0,0.60);
0149 hUEemcalEta_Cent20_50->Draw("ColZ");
0150 hUEemcalEta_Cent20_50->ProfileX()->Draw("same");
0151
0152 can4->cd(3);
0153 gPad->SetLeftMargin(0.15);
0154 gPad->SetRightMargin(.10);
0155 gPad->SetTopMargin(.07);
0156 gPad->SetBottomMargin(.10);
0157 SetHist(hUEemcalEta_Cent50_100,"Eta","UE Size/Tower (GeV)",1);
0158 hUEemcalEta_Cent50_100->SetTitle("UE Size vs Eta 50 - 100 % Centrality Range EMCal");
0159 hUEemcalEta_Cent50_100->GetYaxis()->SetRangeUser(0.0,0.30);
0160 hUEemcalEta_Cent50_100->Draw("ColZ");
0161 hUEemcalEta_Cent50_100->ProfileX()->Draw("same");
0162
0163
0164 TProfile *IH0_20 = hUEiHcalEta_Cent0_20->ProfileX();
0165 TProfile *IH20_50 = hUEiHcalEta_Cent20_50->ProfileX();
0166 TProfile *IH50_100 = hUEiHcalEta_Cent50_100->ProfileX();
0167 TProfile *OH0_20 = hUEoHcalEta_Cent0_20->ProfileX();
0168 TProfile *OH20_50 = hUEoHcalEta_Cent20_50->ProfileX();
0169 TProfile *OH50_100 = hUEoHcalEta_Cent50_100->ProfileX();
0170 TProfile *EM0_20 = hUEemcalEta_Cent0_20->ProfileX();
0171 TProfile *EM20_50 = hUEemcalEta_Cent20_50->ProfileX();
0172 TProfile *EM50_100 = hUEemcalEta_Cent50_100->ProfileX();
0173
0174 TCanvas* can5 = new TCanvas("can5","",1200,400);
0175 can5->Divide(3,1,0,0);
0176 can5->cd(1);
0177 gPad->SetLeftMargin(0.15);
0178 gPad->SetRightMargin(.10);
0179 gPad->SetTopMargin(.07);
0180 gPad->SetBottomMargin(.10);
0181 IH0_20->SetTitle("Mean Values of UE Size vs Eta iHCal");
0182 IH0_20->GetYaxis()->SetTitle("UE Size/Tower (GeV)");
0183 IH0_20->GetYaxis()->SetRangeUser(0.0,0.08);
0184 IH0_20->GetXaxis()->SetTitle("Eta");
0185 IH0_20->Draw();
0186 IH20_50->SetMarkerColor(1);
0187 IH50_100->SetMarkerColor(2);
0188 IH20_50->Draw("same");
0189 IH50_100->Draw("same");
0190
0191 TLegend* leg1 = new TLegend(0.4,0.7,0.7,0.85);
0192 SetLeg(leg1);
0193 leg1->AddEntry(IH0_20,"0-20 % Centrality","P");
0194 leg1->AddEntry(IH20_50,"20-50 % Centrality","P");
0195 leg1->AddEntry(IH50_100,"50-100 % Centrality","P");
0196 leg1->Draw();
0197
0198
0199 can5->cd(2);
0200 gPad->SetLeftMargin(0.15);
0201 gPad->SetRightMargin(.10);
0202 gPad->SetTopMargin(.07);
0203 gPad->SetBottomMargin(.10);
0204 OH0_20->SetTitle("Mean Values of UE Size vs Eta oHCal");
0205 OH0_20->GetYaxis()->SetTitle("UE Size/Tower (GeV)");
0206 OH0_20->GetXaxis()->SetTitle("Eta");
0207 OH0_20->GetYaxis()->SetRangeUser(0.0,0.30);
0208 OH0_20->Draw();
0209 OH20_50->SetMarkerColor(1);
0210 OH50_100->SetMarkerColor(2);
0211 OH20_50->Draw("same");
0212 OH50_100->Draw("same");
0213
0214 TLegend* leg2 = new TLegend(0.4,0.7,0.7,0.85);
0215 SetLeg(leg2);
0216 leg2->AddEntry(OH0_20,"0-20 % Centrality","P"); leg2->AddEntry(OH20_50,"20-50 % Centrality","P"); leg2->AddEntry(OH50_100,"50-100 % Centrality","P"); leg2->Draw();
0217
0218 can5->cd(3);
0219 gPad->SetLeftMargin(0.15);
0220 gPad->SetRightMargin(.10);
0221 gPad->SetTopMargin(.07);
0222 gPad->SetBottomMargin(.10);
0223 EM0_20->SetTitle("Mean Values of UE Size vs Eta EMCal");
0224 EM0_20->GetYaxis()->SetTitle("UE Size/Tower (GeV)");
0225 EM0_20->GetXaxis()->SetTitle("Eta");
0226 EM0_20->GetYaxis()->SetRangeUser(0.0,0.70);
0227 EM0_20->Draw();
0228 EM20_50->SetMarkerColor(1);
0229 EM50_100->SetMarkerColor(2);
0230 EM20_50->Draw("same");
0231 EM50_100->Draw("same");
0232
0233
0234 TLegend* leg3 = new TLegend(0.4,0.7,0.7,0.85);
0235 SetLeg(leg3);
0236 leg3->AddEntry(EM0_20,"0-20 % Centrality","P"); leg3->AddEntry(EM20_50,"20-50 % Centrality","P"); leg3->AddEntry(EM50_100,"50-100 % Centrality","P"); leg3->Draw();
0237 }
0238
0239 void SetHist(TH2* h, string xt ="", string yt ="",int color = 4, int marker = 24,int width = 3, float size = 1.0)
0240 {
0241 h->SetLineWidth(width);
0242 h->SetLineColor(color);
0243 h->SetMarkerColor(3);
0244 h->SetMarkerSize(0.5);
0245 h->SetMarkerStyle(marker);
0246 h->GetYaxis()->SetTitle(yt.c_str());
0247 h->GetYaxis()->SetTitleOffset(1.6);
0248 h->GetXaxis()->SetTitle(xt.c_str());
0249 }
0250
0251 void SetHist(TH2* h, int color = 4)
0252 {
0253 h->SetLineWidth(3);
0254 h->SetLineColor(color);
0255 h->SetMarkerColor(color);
0256 h->SetMarkerSize(0.25);
0257 h->GetYaxis()->SetTitleOffset(1.6);
0258
0259 }
0260
0261
0262 void SetLeg(TLegend* l,float txtsize=0.04){
0263 l->SetBorderSize(1);
0264 l->SetFillColor(0);
0265 l->SetTextSize(txtsize);
0266 }
