void PlotPointingBias() { //=========Macro generated from canvas: c_m/c //========= (Mon May 16 09:09:49 2022) by ROOT version 6.20/06 TCanvas *c_m = new TCanvas("c_m", "c",0,0,1600,1200); gStyle->SetOptStat(0); c_m->Range(-937.5,-250,937.5,250); c_m->SetFillColor(0); c_m->SetBorderMode(0); c_m->SetBorderSize(2); c_m->SetFrameBorderMode(0); c_m->SetFrameBorderMode(0); THStack *hs_m = new THStack(); hs_m->SetName("hs_m"); hs_m->SetTitle(";z_{origin} [mm]; Pointing Bias [mm];"); hs_m->SetMinimum(-200); hs_m->SetMaximum(200); Double_t xAxis6[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_m_stack_2 = new TH1F("hs_m_stack_2","",22, xAxis6); hs_m_stack_2->SetMinimum(-200); hs_m_stack_2->SetMaximum(200); hs_m_stack_2->SetDirectory(0); hs_m_stack_2->SetStats(0); Int_t ci; // for color index setting TColor *color; // for color definition with alpha ci = TColor::GetColor("#000099"); hs_m_stack_2->SetLineColor(ci); hs_m_stack_2->GetXaxis()->SetTitle("z_{origin} [mm]"); hs_m_stack_2->GetXaxis()->SetLabelFont(42); hs_m_stack_2->GetXaxis()->SetTitleOffset(1); hs_m_stack_2->GetXaxis()->SetTitleFont(42); hs_m_stack_2->GetYaxis()->SetTitle(" Pointing Bias [mm]"); hs_m_stack_2->GetYaxis()->SetLabelFont(42); hs_m_stack_2->GetYaxis()->SetTitleFont(42); hs_m_stack_2->GetZaxis()->SetLabelFont(42); hs_m_stack_2->GetZaxis()->SetTitleOffset(1); hs_m_stack_2->GetZaxis()->SetTitleFont(42); hs_m->SetHistogram(hs_m_stack_2); Double_t xAxis7[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 *hmean_Hino_135_10ns_stack_1 = new TH1F("hmean_Hino_135_10ns_stack_1","",22, xAxis7); hmean_Hino_135_10ns_stack_1->SetBinContent(1,155.7873); hmean_Hino_135_10ns_stack_1->SetBinContent(2,94.46177); hmean_Hino_135_10ns_stack_1->SetBinContent(3,67.31615); hmean_Hino_135_10ns_stack_1->SetBinContent(4,52.35643); hmean_Hino_135_10ns_stack_1->SetBinContent(5,39.36731); hmean_Hino_135_10ns_stack_1->SetBinContent(6,30.02062); hmean_Hino_135_10ns_stack_1->SetBinContent(7,25.63226); hmean_Hino_135_10ns_stack_1->SetBinContent(8,16.66452); hmean_Hino_135_10ns_stack_1->SetBinContent(9,16.32619); hmean_Hino_135_10ns_stack_1->SetBinContent(10,9.742337); hmean_Hino_135_10ns_stack_1->SetBinContent(11,4.542204); hmean_Hino_135_10ns_stack_1->SetBinContent(12,-1.892682); hmean_Hino_135_10ns_stack_1->SetBinContent(13,-5.369662); hmean_Hino_135_10ns_stack_1->SetBinContent(14,-13.1304); hmean_Hino_135_10ns_stack_1->SetBinContent(15,-17.48446); hmean_Hino_135_10ns_stack_1->SetBinContent(16,-25.75374); hmean_Hino_135_10ns_stack_1->SetBinContent(17,-28.36901); hmean_Hino_135_10ns_stack_1->SetBinContent(18,-35.96225); hmean_Hino_135_10ns_stack_1->SetBinContent(19,-46.42785); hmean_Hino_135_10ns_stack_1->SetBinContent(20,-63.39075); hmean_Hino_135_10ns_stack_1->SetBinContent(21,-86.72802); hmean_Hino_135_10ns_stack_1->SetBinContent(22,-148.7264); hmean_Hino_135_10ns_stack_1->SetBinError(1,1.855414); hmean_Hino_135_10ns_stack_1->SetBinError(2,1.147195); hmean_Hino_135_10ns_stack_1->SetBinError(3,0.9558603); hmean_Hino_135_10ns_stack_1->SetBinError(4,0.8669459); hmean_Hino_135_10ns_stack_1->SetBinError(5,0.6206155); hmean_Hino_135_10ns_stack_1->SetBinError(6,0.5733142); hmean_Hino_135_10ns_stack_1->SetBinError(7,0.4026507); hmean_Hino_135_10ns_stack_1->SetBinError(8,0.3785535); hmean_Hino_135_10ns_stack_1->SetBinError(9,0.2497781); hmean_Hino_135_10ns_stack_1->SetBinError(10,0.1951774); hmean_Hino_135_10ns_stack_1->SetBinError(11,0.15266); hmean_Hino_135_10ns_stack_1->SetBinError(12,0.1735797); hmean_Hino_135_10ns_stack_1->SetBinError(13,0.2385695); hmean_Hino_135_10ns_stack_1->SetBinError(14,0.2799898); hmean_Hino_135_10ns_stack_1->SetBinError(15,0.3325641); hmean_Hino_135_10ns_stack_1->SetBinError(16,0.4657469); hmean_Hino_135_10ns_stack_1->SetBinError(17,0.577792); hmean_Hino_135_10ns_stack_1->SetBinError(18,0.480839); hmean_Hino_135_10ns_stack_1->SetBinError(19,1.011702); hmean_Hino_135_10ns_stack_1->SetBinError(20,0.9685172); hmean_Hino_135_10ns_stack_1->SetBinError(21,1.151406); hmean_Hino_135_10ns_stack_1->SetBinError(22,1.563923); hmean_Hino_135_10ns_stack_1->SetEntries(22); ci = TColor::GetColor("#ff00ff"); hmean_Hino_135_10ns_stack_1->SetLineColor(ci); ci = TColor::GetColor("#ff00ff"); hmean_Hino_135_10ns_stack_1->SetMarkerColor(ci); hmean_Hino_135_10ns_stack_1->SetMarkerStyle(26); hmean_Hino_135_10ns_stack_1->SetMarkerSize(2); hmean_Hino_135_10ns_stack_1->GetXaxis()->SetLabelFont(42); hmean_Hino_135_10ns_stack_1->GetXaxis()->SetTitleOffset(1); hmean_Hino_135_10ns_stack_1->GetXaxis()->SetTitleFont(42); hmean_Hino_135_10ns_stack_1->GetYaxis()->SetLabelFont(42); hmean_Hino_135_10ns_stack_1->GetYaxis()->SetTitleFont(42); hmean_Hino_135_10ns_stack_1->GetZaxis()->SetLabelFont(42); hmean_Hino_135_10ns_stack_1->GetZaxis()->SetTitleOffset(1); hmean_Hino_135_10ns_stack_1->GetZaxis()->SetTitleFont(42); hs_m->Add(hmean_Hino_135_10ns_stack_1,""); Double_t xAxis8[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 *hmean_Hino_135_2ns_stack_2 = new TH1F("hmean_Hino_135_2ns_stack_2","",22, xAxis8); hmean_Hino_135_2ns_stack_2->SetBinContent(1,137.3021); hmean_Hino_135_2ns_stack_2->SetBinContent(2,85.45183); hmean_Hino_135_2ns_stack_2->SetBinContent(3,60.46997); hmean_Hino_135_2ns_stack_2->SetBinContent(4,52.5466); hmean_Hino_135_2ns_stack_2->SetBinContent(5,35.15818); hmean_Hino_135_2ns_stack_2->SetBinContent(6,27.21475); hmean_Hino_135_2ns_stack_2->SetBinContent(7,22.95615); hmean_Hino_135_2ns_stack_2->SetBinContent(8,17.98761); hmean_Hino_135_2ns_stack_2->SetBinContent(9,13.79463); hmean_Hino_135_2ns_stack_2->SetBinContent(10,8.701632); hmean_Hino_135_2ns_stack_2->SetBinContent(11,3.905642); hmean_Hino_135_2ns_stack_2->SetBinContent(12,-1.352747); hmean_Hino_135_2ns_stack_2->SetBinContent(13,-6.124126); hmean_Hino_135_2ns_stack_2->SetBinContent(14,-11.25723); hmean_Hino_135_2ns_stack_2->SetBinContent(15,-15.85189); hmean_Hino_135_2ns_stack_2->SetBinContent(16,-20.67181); hmean_Hino_135_2ns_stack_2->SetBinContent(17,-24.84672); hmean_Hino_135_2ns_stack_2->SetBinContent(18,-33.7615); hmean_Hino_135_2ns_stack_2->SetBinContent(19,-46.43964); hmean_Hino_135_2ns_stack_2->SetBinContent(20,-59.56991); hmean_Hino_135_2ns_stack_2->SetBinContent(21,-80.2435); hmean_Hino_135_2ns_stack_2->SetBinContent(22,-128.963); hmean_Hino_135_2ns_stack_2->SetBinError(1,1.04143); hmean_Hino_135_2ns_stack_2->SetBinError(2,0.6939977); hmean_Hino_135_2ns_stack_2->SetBinError(3,0.6844119); hmean_Hino_135_2ns_stack_2->SetBinError(4,0.4390868); hmean_Hino_135_2ns_stack_2->SetBinError(5,0.2877842); hmean_Hino_135_2ns_stack_2->SetBinError(6,0.3192736); hmean_Hino_135_2ns_stack_2->SetBinError(7,0.2411902); hmean_Hino_135_2ns_stack_2->SetBinError(8,0.1875558); hmean_Hino_135_2ns_stack_2->SetBinError(9,0.1323023); hmean_Hino_135_2ns_stack_2->SetBinError(10,0.1009348); hmean_Hino_135_2ns_stack_2->SetBinError(11,0.08853987); hmean_Hino_135_2ns_stack_2->SetBinError(12,0.09979649); hmean_Hino_135_2ns_stack_2->SetBinError(13,0.106231); hmean_Hino_135_2ns_stack_2->SetBinError(14,0.129392); hmean_Hino_135_2ns_stack_2->SetBinError(15,0.1920774); hmean_Hino_135_2ns_stack_2->SetBinError(16,0.2452257); hmean_Hino_135_2ns_stack_2->SetBinError(17,0.3078361); hmean_Hino_135_2ns_stack_2->SetBinError(18,0.3237766); hmean_Hino_135_2ns_stack_2->SetBinError(19,0.5093541); hmean_Hino_135_2ns_stack_2->SetBinError(20,0.6468439); hmean_Hino_135_2ns_stack_2->SetBinError(21,0.7057011); hmean_Hino_135_2ns_stack_2->SetBinError(22,1.01096); hmean_Hino_135_2ns_stack_2->SetEntries(22); ci = TColor::GetColor("#0000ff"); hmean_Hino_135_2ns_stack_2->SetLineColor(ci); ci = TColor::GetColor("#0000ff"); hmean_Hino_135_2ns_stack_2->SetMarkerColor(ci); hmean_Hino_135_2ns_stack_2->SetMarkerStyle(26); hmean_Hino_135_2ns_stack_2->SetMarkerSize(2); hmean_Hino_135_2ns_stack_2->GetXaxis()->SetLabelFont(42); hmean_Hino_135_2ns_stack_2->GetXaxis()->SetTitleOffset(1); hmean_Hino_135_2ns_stack_2->GetXaxis()->SetTitleFont(42); hmean_Hino_135_2ns_stack_2->GetYaxis()->SetLabelFont(42); hmean_Hino_135_2ns_stack_2->GetYaxis()->SetTitleFont(42); hmean_Hino_135_2ns_stack_2->GetZaxis()->SetLabelFont(42); hmean_Hino_135_2ns_stack_2->GetZaxis()->SetTitleOffset(1); hmean_Hino_135_2ns_stack_2->GetZaxis()->SetTitleFont(42); hs_m->Add(hmean_Hino_135_2ns_stack_2,""); Double_t xAxis9[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 *hmean_Hino_135_Hyy_stack_3 = new TH1F("hmean_Hino_135_Hyy_stack_3","",22, xAxis9); hmean_Hino_135_Hyy_stack_3->SetBinContent(1,144.2077); hmean_Hino_135_Hyy_stack_3->SetBinContent(2,91.05521); hmean_Hino_135_Hyy_stack_3->SetBinContent(3,60.76589); hmean_Hino_135_Hyy_stack_3->SetBinContent(4,51.38012); hmean_Hino_135_Hyy_stack_3->SetBinContent(5,36.20037); hmean_Hino_135_Hyy_stack_3->SetBinContent(6,28.72125); hmean_Hino_135_Hyy_stack_3->SetBinContent(7,24.11675); hmean_Hino_135_Hyy_stack_3->SetBinContent(8,17.97845); hmean_Hino_135_Hyy_stack_3->SetBinContent(9,13.74108); hmean_Hino_135_Hyy_stack_3->SetBinContent(10,8.668206); hmean_Hino_135_Hyy_stack_3->SetBinContent(11,3.594204); hmean_Hino_135_Hyy_stack_3->SetBinContent(12,-1.95629); hmean_Hino_135_Hyy_stack_3->SetBinContent(13,-5.464774); hmean_Hino_135_Hyy_stack_3->SetBinContent(14,-11.70699); hmean_Hino_135_Hyy_stack_3->SetBinContent(15,-16.17486); hmean_Hino_135_Hyy_stack_3->SetBinContent(16,-21.47552); hmean_Hino_135_Hyy_stack_3->SetBinContent(17,-26.77675); hmean_Hino_135_Hyy_stack_3->SetBinContent(18,-32.9808); hmean_Hino_135_Hyy_stack_3->SetBinContent(19,-47.40245); hmean_Hino_135_Hyy_stack_3->SetBinContent(20,-62.06454); hmean_Hino_135_Hyy_stack_3->SetBinContent(21,-84.28474); hmean_Hino_135_Hyy_stack_3->SetBinContent(22,-129.7012); hmean_Hino_135_Hyy_stack_3->SetBinError(1,1.001125); hmean_Hino_135_Hyy_stack_3->SetBinError(2,0.7029557); hmean_Hino_135_Hyy_stack_3->SetBinError(3,0.6775733); hmean_Hino_135_Hyy_stack_3->SetBinError(4,0.5681656); hmean_Hino_135_Hyy_stack_3->SetBinError(5,0.3288972); hmean_Hino_135_Hyy_stack_3->SetBinError(6,0.3347657); hmean_Hino_135_Hyy_stack_3->SetBinError(7,0.2347389); hmean_Hino_135_Hyy_stack_3->SetBinError(8,0.2170261); hmean_Hino_135_Hyy_stack_3->SetBinError(9,0.1414754); hmean_Hino_135_Hyy_stack_3->SetBinError(10,0.1109624); hmean_Hino_135_Hyy_stack_3->SetBinError(11,0.09075633); hmean_Hino_135_Hyy_stack_3->SetBinError(12,0.1053682); hmean_Hino_135_Hyy_stack_3->SetBinError(13,0.120986); hmean_Hino_135_Hyy_stack_3->SetBinError(14,0.1383942); hmean_Hino_135_Hyy_stack_3->SetBinError(15,0.2041889); hmean_Hino_135_Hyy_stack_3->SetBinError(16,0.2592609); hmean_Hino_135_Hyy_stack_3->SetBinError(17,0.3575097); hmean_Hino_135_Hyy_stack_3->SetBinError(18,0.399888); hmean_Hino_135_Hyy_stack_3->SetBinError(19,0.4749649); hmean_Hino_135_Hyy_stack_3->SetBinError(20,0.6811654); hmean_Hino_135_Hyy_stack_3->SetBinError(21,0.7420164); hmean_Hino_135_Hyy_stack_3->SetBinError(22,0.9416993); hmean_Hino_135_Hyy_stack_3->SetEntries(22); ci = TColor::GetColor("#ff0000"); hmean_Hino_135_Hyy_stack_3->SetLineColor(ci); ci = TColor::GetColor("#ff0000"); hmean_Hino_135_Hyy_stack_3->SetMarkerColor(ci); hmean_Hino_135_Hyy_stack_3->SetMarkerStyle(24); hmean_Hino_135_Hyy_stack_3->SetMarkerSize(2); hmean_Hino_135_Hyy_stack_3->GetXaxis()->SetLabelFont(42); hmean_Hino_135_Hyy_stack_3->GetXaxis()->SetTitleOffset(1); hmean_Hino_135_Hyy_stack_3->GetXaxis()->SetTitleFont(42); hmean_Hino_135_Hyy_stack_3->GetYaxis()->SetLabelFont(42); hmean_Hino_135_Hyy_stack_3->GetYaxis()->SetTitleFont(42); hmean_Hino_135_Hyy_stack_3->GetZaxis()->SetLabelFont(42); hmean_Hino_135_Hyy_stack_3->GetZaxis()->SetTitleOffset(1); hmean_Hino_135_Hyy_stack_3->GetZaxis()->SetTitleFont(42); hs_m->Add(hmean_Hino_135_Hyy_stack_3,""); Double_t xAxis10[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 *hmean_Zee_stack_4 = new TH1F("hmean_Zee_stack_4","",22, xAxis10); hmean_Zee_stack_4->SetBinContent(6,31.29541); hmean_Zee_stack_4->SetBinContent(7,23.09946); hmean_Zee_stack_4->SetBinContent(8,18.91688); hmean_Zee_stack_4->SetBinContent(9,14.40222); hmean_Zee_stack_4->SetBinContent(10,9.203378); hmean_Zee_stack_4->SetBinContent(11,4.232938); hmean_Zee_stack_4->SetBinContent(12,-1.207011); hmean_Zee_stack_4->SetBinContent(13,-6.176314); hmean_Zee_stack_4->SetBinContent(14,-10.9162); hmean_Zee_stack_4->SetBinContent(15,-15.50382); hmean_Zee_stack_4->SetBinContent(16,-20.80232); hmean_Zee_stack_4->SetBinContent(17,-24.81242); hmean_Zee_stack_4->SetBinError(6,1.121656); hmean_Zee_stack_4->SetBinError(7,0.2480441); hmean_Zee_stack_4->SetBinError(8,0.09694975); hmean_Zee_stack_4->SetBinError(9,0.04781965); hmean_Zee_stack_4->SetBinError(10,0.03018516); hmean_Zee_stack_4->SetBinError(11,0.02424172); hmean_Zee_stack_4->SetBinError(12,0.0237044); hmean_Zee_stack_4->SetBinError(13,0.03050967); hmean_Zee_stack_4->SetBinError(14,0.04744555); hmean_Zee_stack_4->SetBinError(15,0.09222348); hmean_Zee_stack_4->SetBinError(16,0.1987363); hmean_Zee_stack_4->SetBinError(17,0.6346914); hmean_Zee_stack_4->SetEntries(12); hmean_Zee_stack_4->SetMarkerStyle(20); hmean_Zee_stack_4->SetMarkerSize(2); hmean_Zee_stack_4->GetXaxis()->SetLabelFont(42); hmean_Zee_stack_4->GetXaxis()->SetTitleOffset(1); hmean_Zee_stack_4->GetXaxis()->SetTitleFont(42); hmean_Zee_stack_4->GetYaxis()->SetLabelFont(42); hmean_Zee_stack_4->GetYaxis()->SetTitleFont(42); hmean_Zee_stack_4->GetZaxis()->SetLabelFont(42); hmean_Zee_stack_4->GetZaxis()->SetTitleOffset(1); hmean_Zee_stack_4->GetZaxis()->SetTitleFont(42); hs_m->Add(hmean_Zee_stack_4,""); hs_m->Draw("nostack p"); TLegend *leg = new TLegend(0.65,0.7,0.82,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("hmean_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("hmean_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("hmean_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("hmean_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_m->Modified(); c_m->cd(); c_m->SetSelected(c_m); }