void PlotPointingRes() { //=========Macro generated from canvas: c/c //========= (Mon May 16 09:09:49 2022) by ROOT version 6.20/06 TCanvas *c = new TCanvas("c", "c",0,0,1600,1200); gStyle->SetOptStat(0); c->Range(-937.5,-18.75,937.5,168.75); c->SetFillColor(0); c->SetBorderMode(0); c->SetBorderSize(2); c->SetFrameBorderMode(0); c->SetFrameBorderMode(0); THStack *hs = new THStack(); hs->SetName("hs"); hs->SetTitle(";z_{origin} [mm]; Pointing Resolution [mm];"); hs->SetMinimum(0); hs->SetMaximum(150); Double_t xAxis1[23] = {-750, -400, -300, -250, -200, -150, -125, -100, -75, -50, -25, 0, 25, 50, 75, 100, 125, 150, 200, 250, 300, 400, 750}; TH1F *hs_stack_1 = new TH1F("hs_stack_1","",22, xAxis1); hs_stack_1->SetMinimum(0); hs_stack_1->SetMaximum(150); hs_stack_1->SetDirectory(0); hs_stack_1->SetStats(0); Int_t ci; // for color index setting TColor *color; // for color definition with alpha ci = TColor::GetColor("#000099"); hs_stack_1->SetLineColor(ci); hs_stack_1->GetXaxis()->SetTitle("z_{origin} [mm]"); hs_stack_1->GetXaxis()->SetLabelFont(42); hs_stack_1->GetXaxis()->SetTitleOffset(1); hs_stack_1->GetXaxis()->SetTitleFont(42); hs_stack_1->GetYaxis()->SetTitle(" Pointing Resolution [mm]"); hs_stack_1->GetYaxis()->SetLabelFont(42); hs_stack_1->GetYaxis()->SetTitleFont(42); hs_stack_1->GetZaxis()->SetLabelFont(42); hs_stack_1->GetZaxis()->SetTitleOffset(1); hs_stack_1->GetZaxis()->SetTitleFont(42); hs->SetHistogram(hs_stack_1); Double_t xAxis2[23] = {-750, -400, -300, -250, -200, -150, -125, -100, -75, -50, -25, 0, 25, 50, 75, 100, 125, 150, 200, 250, 300, 400, 750}; TH1F *hsig_Hino_135_10ns_stack_1 = new TH1F("hsig_Hino_135_10ns_stack_1","",22, xAxis2); hsig_Hino_135_10ns_stack_1->SetBinContent(1,120.9265); hsig_Hino_135_10ns_stack_1->SetBinContent(2,77.62166); hsig_Hino_135_10ns_stack_1->SetBinContent(3,52.98224); hsig_Hino_135_10ns_stack_1->SetBinContent(4,57.36145); hsig_Hino_135_10ns_stack_1->SetBinContent(5,42.40071); hsig_Hino_135_10ns_stack_1->SetBinContent(6,32.64183); hsig_Hino_135_10ns_stack_1->SetBinContent(7,20.42934); hsig_Hino_135_10ns_stack_1->SetBinContent(8,20.43621); hsig_Hino_135_10ns_stack_1->SetBinContent(9,14.33404); hsig_Hino_135_10ns_stack_1->SetBinContent(10,11.61202); hsig_Hino_135_10ns_stack_1->SetBinContent(11,8.859221); hsig_Hino_135_10ns_stack_1->SetBinContent(12,10.73035); hsig_Hino_135_10ns_stack_1->SetBinContent(13,14.32106); hsig_Hino_135_10ns_stack_1->SetBinContent(14,16.40014); hsig_Hino_135_10ns_stack_1->SetBinContent(15,17.67178); hsig_Hino_135_10ns_stack_1->SetBinContent(16,23.50431); hsig_Hino_135_10ns_stack_1->SetBinContent(17,28.89335); hsig_Hino_135_10ns_stack_1->SetBinContent(18,32.05732); hsig_Hino_135_10ns_stack_1->SetBinContent(19,56.29406); hsig_Hino_135_10ns_stack_1->SetBinContent(20,56.46368); hsig_Hino_135_10ns_stack_1->SetBinContent(21,74.47071); hsig_Hino_135_10ns_stack_1->SetBinContent(22,110.021); hsig_Hino_135_10ns_stack_1->SetBinError(1,2.586166); hsig_Hino_135_10ns_stack_1->SetBinError(2,1.464358); hsig_Hino_135_10ns_stack_1->SetBinError(3,1.160547); hsig_Hino_135_10ns_stack_1->SetBinError(4,1.00079); hsig_Hino_135_10ns_stack_1->SetBinError(5,0.7629911); hsig_Hino_135_10ns_stack_1->SetBinError(6,0.6334668); hsig_Hino_135_10ns_stack_1->SetBinError(7,0.4778916); hsig_Hino_135_10ns_stack_1->SetBinError(8,0.5031314); hsig_Hino_135_10ns_stack_1->SetBinError(9,0.3267311); hsig_Hino_135_10ns_stack_1->SetBinError(10,0.2414033); hsig_Hino_135_10ns_stack_1->SetBinError(11,0.1928448); hsig_Hino_135_10ns_stack_1->SetBinError(12,0.2373107); hsig_Hino_135_10ns_stack_1->SetBinError(13,0.322394); hsig_Hino_135_10ns_stack_1->SetBinError(14,0.3905554); hsig_Hino_135_10ns_stack_1->SetBinError(15,0.4166364); hsig_Hino_135_10ns_stack_1->SetBinError(16,0.5563527); hsig_Hino_135_10ns_stack_1->SetBinError(17,0.7585679); hsig_Hino_135_10ns_stack_1->SetBinError(18,0.6196121); hsig_Hino_135_10ns_stack_1->SetBinError(19,1.512515); hsig_Hino_135_10ns_stack_1->SetBinError(20,1.144912); hsig_Hino_135_10ns_stack_1->SetBinError(21,1.471564); hsig_Hino_135_10ns_stack_1->SetBinError(22,1.989396); hsig_Hino_135_10ns_stack_1->SetEntries(22); ci = TColor::GetColor("#ff00ff"); hsig_Hino_135_10ns_stack_1->SetLineColor(ci); ci = TColor::GetColor("#ff00ff"); hsig_Hino_135_10ns_stack_1->SetMarkerColor(ci); hsig_Hino_135_10ns_stack_1->SetMarkerStyle(26); hsig_Hino_135_10ns_stack_1->SetMarkerSize(2); hsig_Hino_135_10ns_stack_1->GetXaxis()->SetTitle("|z_{origin}| [mm]"); hsig_Hino_135_10ns_stack_1->GetXaxis()->SetLabelFont(42); hsig_Hino_135_10ns_stack_1->GetXaxis()->SetTitleOffset(1); hsig_Hino_135_10ns_stack_1->GetXaxis()->SetTitleFont(42); hsig_Hino_135_10ns_stack_1->GetYaxis()->SetTitle(" Pointing Resolution [mm]"); hsig_Hino_135_10ns_stack_1->GetYaxis()->SetLabelFont(42); hsig_Hino_135_10ns_stack_1->GetYaxis()->SetTitleFont(42); hsig_Hino_135_10ns_stack_1->GetZaxis()->SetLabelFont(42); hsig_Hino_135_10ns_stack_1->GetZaxis()->SetTitleOffset(1); hsig_Hino_135_10ns_stack_1->GetZaxis()->SetTitleFont(42); hs->Add(hsig_Hino_135_10ns_stack_1,""); Double_t xAxis3[23] = {-750, -400, -300, -250, -200, -150, -125, -100, -75, -50, -25, 0, 25, 50, 75, 100, 125, 150, 200, 250, 300, 400, 750}; TH1F *hsig_Hino_135_2ns_stack_2 = new TH1F("hsig_Hino_135_2ns_stack_2","",22, xAxis3); hsig_Hino_135_2ns_stack_2->SetBinContent(1,114.3606); hsig_Hino_135_2ns_stack_2->SetBinContent(2,77.69562); hsig_Hino_135_2ns_stack_2->SetBinContent(3,64.13433); hsig_Hino_135_2ns_stack_2->SetBinContent(4,49.57817); hsig_Hino_135_2ns_stack_2->SetBinContent(5,36.95729); hsig_Hino_135_2ns_stack_2->SetBinContent(6,30.781); hsig_Hino_135_2ns_stack_2->SetBinContent(7,23.70428); hsig_Hino_135_2ns_stack_2->SetBinContent(8,19.70467); hsig_Hino_135_2ns_stack_2->SetBinContent(9,14.61741); hsig_Hino_135_2ns_stack_2->SetBinContent(10,11.72085); hsig_Hino_135_2ns_stack_2->SetBinContent(11,10.2396); hsig_Hino_135_2ns_stack_2->SetBinContent(12,11.49671); hsig_Hino_135_2ns_stack_2->SetBinContent(13,12.21777); hsig_Hino_135_2ns_stack_2->SetBinContent(14,14.53042); hsig_Hino_135_2ns_stack_2->SetBinContent(15,19.98182); hsig_Hino_135_2ns_stack_2->SetBinContent(16,23.41485); hsig_Hino_135_2ns_stack_2->SetBinContent(17,30.14566); hsig_Hino_135_2ns_stack_2->SetBinContent(18,39.53489); hsig_Hino_135_2ns_stack_2->SetBinContent(19,53.26426); hsig_Hino_135_2ns_stack_2->SetBinContent(20,61.01122); hsig_Hino_135_2ns_stack_2->SetBinContent(21,82.08636); hsig_Hino_135_2ns_stack_2->SetBinContent(22,115.046); hsig_Hino_135_2ns_stack_2->SetBinError(1,1.306249); hsig_Hino_135_2ns_stack_2->SetBinError(2,0.8624832); hsig_Hino_135_2ns_stack_2->SetBinError(3,0.8713448); hsig_Hino_135_2ns_stack_2->SetBinError(4,0.5281473); hsig_Hino_135_2ns_stack_2->SetBinError(5,0.348496); hsig_Hino_135_2ns_stack_2->SetBinError(6,0.3711733); hsig_Hino_135_2ns_stack_2->SetBinError(7,0.2953853); hsig_Hino_135_2ns_stack_2->SetBinError(8,0.2501384); hsig_Hino_135_2ns_stack_2->SetBinError(9,0.1735626); hsig_Hino_135_2ns_stack_2->SetBinError(10,0.130568); hsig_Hino_135_2ns_stack_2->SetBinError(11,0.1168149); hsig_Hino_135_2ns_stack_2->SetBinError(12,0.1460116); hsig_Hino_135_2ns_stack_2->SetBinError(13,0.1407996); hsig_Hino_135_2ns_stack_2->SetBinError(14,0.1646287); hsig_Hino_135_2ns_stack_2->SetBinError(15,0.2536571); hsig_Hino_135_2ns_stack_2->SetBinError(16,0.302886); hsig_Hino_135_2ns_stack_2->SetBinError(17,0.3789356); hsig_Hino_135_2ns_stack_2->SetBinError(18,0.4028325); hsig_Hino_135_2ns_stack_2->SetBinError(19,0.6566924); hsig_Hino_135_2ns_stack_2->SetBinError(20,0.775235); hsig_Hino_135_2ns_stack_2->SetBinError(21,0.8291415); hsig_Hino_135_2ns_stack_2->SetBinError(22,1.282282); hsig_Hino_135_2ns_stack_2->SetEntries(22); ci = TColor::GetColor("#0000ff"); hsig_Hino_135_2ns_stack_2->SetLineColor(ci); ci = TColor::GetColor("#0000ff"); hsig_Hino_135_2ns_stack_2->SetMarkerColor(ci); hsig_Hino_135_2ns_stack_2->SetMarkerStyle(26); hsig_Hino_135_2ns_stack_2->SetMarkerSize(2); hsig_Hino_135_2ns_stack_2->GetXaxis()->SetTitle("|z_{origin}| [mm]"); hsig_Hino_135_2ns_stack_2->GetXaxis()->SetLabelFont(42); hsig_Hino_135_2ns_stack_2->GetXaxis()->SetTitleOffset(1); hsig_Hino_135_2ns_stack_2->GetXaxis()->SetTitleFont(42); hsig_Hino_135_2ns_stack_2->GetYaxis()->SetTitle(" Pointing Resolution [mm]"); hsig_Hino_135_2ns_stack_2->GetYaxis()->SetLabelFont(42); hsig_Hino_135_2ns_stack_2->GetYaxis()->SetTitleFont(42); hsig_Hino_135_2ns_stack_2->GetZaxis()->SetLabelFont(42); hsig_Hino_135_2ns_stack_2->GetZaxis()->SetTitleOffset(1); hsig_Hino_135_2ns_stack_2->GetZaxis()->SetTitleFont(42); hs->Add(hsig_Hino_135_2ns_stack_2,""); Double_t xAxis4[23] = {-750, -400, -300, -250, -200, -150, -125, -100, -75, -50, -25, 0, 25, 50, 75, 100, 125, 150, 200, 250, 300, 400, 750}; TH1F *hsig_Hino_135_Hyy_stack_3 = new TH1F("hsig_Hino_135_Hyy_stack_3","",22, xAxis4); hsig_Hino_135_Hyy_stack_3->SetBinContent(1,117.859); hsig_Hino_135_Hyy_stack_3->SetBinContent(2,81.39754); hsig_Hino_135_Hyy_stack_3->SetBinContent(3,63.04055); hsig_Hino_135_Hyy_stack_3->SetBinContent(4,57.8531); hsig_Hino_135_Hyy_stack_3->SetBinContent(5,38.07898); hsig_Hino_135_Hyy_stack_3->SetBinContent(6,32.16523); hsig_Hino_135_Hyy_stack_3->SetBinContent(7,20.9966); hsig_Hino_135_Hyy_stack_3->SetBinContent(8,20.50955); hsig_Hino_135_Hyy_stack_3->SetBinContent(9,14.39726); hsig_Hino_135_Hyy_stack_3->SetBinContent(10,11.96467); hsig_Hino_135_Hyy_stack_3->SetBinContent(11,9.993485); hsig_Hino_135_Hyy_stack_3->SetBinContent(12,11.45258); hsig_Hino_135_Hyy_stack_3->SetBinContent(13,13.21957); hsig_Hino_135_Hyy_stack_3->SetBinContent(14,13.88452); hsig_Hino_135_Hyy_stack_3->SetBinContent(15,19.08692); hsig_Hino_135_Hyy_stack_3->SetBinContent(16,23.00681); hsig_Hino_135_Hyy_stack_3->SetBinContent(17,28.3835); hsig_Hino_135_Hyy_stack_3->SetBinContent(18,40.12212); hsig_Hino_135_Hyy_stack_3->SetBinContent(19,51.52132); hsig_Hino_135_Hyy_stack_3->SetBinContent(20,62.56757); hsig_Hino_135_Hyy_stack_3->SetBinContent(21,81.73407); hsig_Hino_135_Hyy_stack_3->SetBinContent(22,117.2871); hsig_Hino_135_Hyy_stack_3->SetBinError(1,1.279771); hsig_Hino_135_Hyy_stack_3->SetBinError(2,0.8416983); hsig_Hino_135_Hyy_stack_3->SetBinError(3,0.8505508); hsig_Hino_135_Hyy_stack_3->SetBinError(4,0.732224); hsig_Hino_135_Hyy_stack_3->SetBinError(5,0.4145166); hsig_Hino_135_Hyy_stack_3->SetBinError(6,0.3663025); hsig_Hino_135_Hyy_stack_3->SetBinError(7,0.2705767); hsig_Hino_135_Hyy_stack_3->SetBinError(8,0.3029827); hsig_Hino_135_Hyy_stack_3->SetBinError(9,0.1975384); hsig_Hino_135_Hyy_stack_3->SetBinError(10,0.1446507); hsig_Hino_135_Hyy_stack_3->SetBinError(11,0.1179982); hsig_Hino_135_Hyy_stack_3->SetBinError(12,0.152962); hsig_Hino_135_Hyy_stack_3->SetBinError(13,0.1615089); hsig_Hino_135_Hyy_stack_3->SetBinError(14,0.186218); hsig_Hino_135_Hyy_stack_3->SetBinError(15,0.275355); hsig_Hino_135_Hyy_stack_3->SetBinError(16,0.3196969); hsig_Hino_135_Hyy_stack_3->SetBinError(17,0.4509352); hsig_Hino_135_Hyy_stack_3->SetBinError(18,0.5695532); hsig_Hino_135_Hyy_stack_3->SetBinError(19,0.5321193); hsig_Hino_135_Hyy_stack_3->SetBinError(20,0.864853); hsig_Hino_135_Hyy_stack_3->SetBinError(21,0.9977375); hsig_Hino_135_Hyy_stack_3->SetBinError(22,1.199031); hsig_Hino_135_Hyy_stack_3->SetEntries(22); ci = TColor::GetColor("#ff0000"); hsig_Hino_135_Hyy_stack_3->SetLineColor(ci); ci = TColor::GetColor("#ff0000"); hsig_Hino_135_Hyy_stack_3->SetMarkerColor(ci); hsig_Hino_135_Hyy_stack_3->SetMarkerStyle(24); hsig_Hino_135_Hyy_stack_3->SetMarkerSize(2); hsig_Hino_135_Hyy_stack_3->GetXaxis()->SetTitle("|z_{origin}| [mm]"); hsig_Hino_135_Hyy_stack_3->GetXaxis()->SetLabelFont(42); hsig_Hino_135_Hyy_stack_3->GetXaxis()->SetTitleOffset(1); hsig_Hino_135_Hyy_stack_3->GetXaxis()->SetTitleFont(42); hsig_Hino_135_Hyy_stack_3->GetYaxis()->SetTitle(" Pointing Resolution [mm]"); hsig_Hino_135_Hyy_stack_3->GetYaxis()->SetLabelFont(42); hsig_Hino_135_Hyy_stack_3->GetYaxis()->SetTitleFont(42); hsig_Hino_135_Hyy_stack_3->GetZaxis()->SetLabelFont(42); hsig_Hino_135_Hyy_stack_3->GetZaxis()->SetTitleOffset(1); hsig_Hino_135_Hyy_stack_3->GetZaxis()->SetTitleFont(42); hs->Add(hsig_Hino_135_Hyy_stack_3,""); Double_t xAxis5[23] = {-750, -400, -300, -250, -200, -150, -125, -100, -75, -50, -25, 0, 25, 50, 75, 100, 125, 150, 200, 250, 300, 400, 750}; TH1F *hsig_Zee_stack_4 = new TH1F("hsig_Zee_stack_4","",22, xAxis5); hsig_Zee_stack_4->SetBinContent(6,30.3194); hsig_Zee_stack_4->SetBinContent(7,21.41113); hsig_Zee_stack_4->SetBinContent(8,16.27529); hsig_Zee_stack_4->SetBinContent(9,13.76733); hsig_Zee_stack_4->SetBinContent(10,12.02766); hsig_Zee_stack_4->SetBinContent(11,11.13904); hsig_Zee_stack_4->SetBinContent(12,10.971); hsig_Zee_stack_4->SetBinContent(13,12.08243); hsig_Zee_stack_4->SetBinContent(14,13.25971); hsig_Zee_stack_4->SetBinContent(15,15.96522); hsig_Zee_stack_4->SetBinContent(16,17.53958); hsig_Zee_stack_4->SetBinContent(17,23.94714); hsig_Zee_stack_4->SetBinError(6,1.480913); hsig_Zee_stack_4->SetBinError(7,0.3359183); hsig_Zee_stack_4->SetBinError(8,0.1242792); hsig_Zee_stack_4->SetBinError(9,0.06090065); hsig_Zee_stack_4->SetBinError(10,0.0398131); hsig_Zee_stack_4->SetBinError(11,0.03227324); hsig_Zee_stack_4->SetBinError(12,0.03126813); hsig_Zee_stack_4->SetBinError(13,0.0405321); hsig_Zee_stack_4->SetBinError(14,0.06297045); hsig_Zee_stack_4->SetBinError(15,0.1195467); hsig_Zee_stack_4->SetBinError(16,0.2569084); hsig_Zee_stack_4->SetBinError(17,0.7873421); hsig_Zee_stack_4->SetEntries(12); hsig_Zee_stack_4->SetMarkerStyle(20); hsig_Zee_stack_4->SetMarkerSize(2); hsig_Zee_stack_4->GetXaxis()->SetTitle("|z_{origin}| [mm]"); hsig_Zee_stack_4->GetXaxis()->SetLabelFont(42); hsig_Zee_stack_4->GetXaxis()->SetTitleOffset(1); hsig_Zee_stack_4->GetXaxis()->SetTitleFont(42); hsig_Zee_stack_4->GetYaxis()->SetTitle(" Pointing Resolution [mm]"); hsig_Zee_stack_4->GetYaxis()->SetLabelFont(42); hsig_Zee_stack_4->GetYaxis()->SetTitleFont(42); hsig_Zee_stack_4->GetZaxis()->SetLabelFont(42); hsig_Zee_stack_4->GetZaxis()->SetTitleOffset(1); hsig_Zee_stack_4->GetZaxis()->SetTitleFont(42); hs->Add(hsig_Zee_stack_4,""); hs->Draw("nostack p"); TLegend *leg = new TLegend(0.55,0.7,0.72,0.88,NULL,"brNDC"); leg->SetBorderSize(0); leg->SetTextSize(0.03); leg->SetLineColor(1); leg->SetLineStyle(1); leg->SetLineWidth(1); leg->SetFillColor(0); leg->SetFillStyle(1001); TLegendEntry *entry=leg->AddEntry("hsig_Hino_135_10ns_stack_1","135GeV 2ns signal","lpf"); entry->SetFillStyle(1001); ci = TColor::GetColor("#ff00ff"); entry->SetLineColor(ci); entry->SetLineStyle(1); entry->SetLineWidth(1); ci = TColor::GetColor("#ff00ff"); entry->SetMarkerColor(ci); entry->SetMarkerStyle(26); entry->SetMarkerSize(2); entry->SetTextFont(42); entry=leg->AddEntry("hsig_Hino_135_2ns_stack_2","135GeV 10ns signal","lpf"); entry->SetFillStyle(1001); ci = TColor::GetColor("#0000ff"); entry->SetLineColor(ci); entry->SetLineStyle(1); entry->SetLineWidth(1); ci = TColor::GetColor("#0000ff"); entry->SetMarkerColor(ci); entry->SetMarkerStyle(26); entry->SetMarkerSize(2); entry->SetTextFont(42); entry=leg->AddEntry("hsig_Hino_135_Hyy_stack_3","135 GeV Hyy signal","lpf"); entry->SetFillStyle(1001); ci = TColor::GetColor("#ff0000"); entry->SetLineColor(ci); entry->SetLineStyle(1); entry->SetLineWidth(1); ci = TColor::GetColor("#ff0000"); entry->SetMarkerColor(ci); entry->SetMarkerStyle(24); entry->SetMarkerSize(2); entry->SetTextFont(42); entry=leg->AddEntry("hsig_Zee_stack_4","Zee MC","lpf"); entry->SetFillStyle(1001); entry->SetLineColor(1); entry->SetLineStyle(1); entry->SetLineWidth(1); entry->SetMarkerColor(1); entry->SetMarkerStyle(20); entry->SetMarkerSize(2); entry->SetTextFont(42); leg->Draw(); TLatex * tex = new TLatex(0.14,0.84,"ATLAS"); tex->SetNDC(); tex->SetTextFont(72); tex->SetTextSize(0.045); tex->SetLineWidth(2); tex->Draw(); tex = new TLatex(0.25,0.84,"Internal"); tex->SetNDC(); tex->SetTextFont(42); tex->SetTextSize(0.045); tex->SetLineWidth(2); tex->Draw(); tex = new TLatex(0.14,0.8,"#sqrt{s} = 13 TeV"); tex->SetNDC(); tex->SetTextFont(42); tex->SetTextSize(0.03375); tex->SetLineWidth(2); tex->Draw(); c->Modified(); c->cd(); c->SetSelected(c); }