void PlotPointingBias()
{
//=========Macro generated from canvas: c_m/c
//=========  (Mon May 16 05:45:10 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(-875,-250,875,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 xAxis5[23] = {-700, -400, -300, -250, -200, -150, -125, -100, -75, -50, -25, 0, 25, 50, 75, 100, 125, 150, 200, 250, 300, 400, 700}; 
   
   TH1F *hs_m_stack_2 = new TH1F("hs_m_stack_2","",22, xAxis5);
   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 xAxis6[23] = {-700, -400, -300, -250, -200, -150, -125, -100, -75, -50, -25, 0, 25, 50, 75, 100, 125, 150, 200, 250, 300, 400, 700}; 
   
   TH1F *hmean_Hino_135_Hyy_stack_1 = new TH1F("hmean_Hino_135_Hyy_stack_1","",22, xAxis6);
   hmean_Hino_135_Hyy_stack_1->SetBinContent(1,141.0901);
   hmean_Hino_135_Hyy_stack_1->SetBinContent(2,92.46945);
   hmean_Hino_135_Hyy_stack_1->SetBinContent(3,62.98962);
   hmean_Hino_135_Hyy_stack_1->SetBinContent(4,49.10697);
   hmean_Hino_135_Hyy_stack_1->SetBinContent(5,34.60996);
   hmean_Hino_135_Hyy_stack_1->SetBinContent(6,28.63077);
   hmean_Hino_135_Hyy_stack_1->SetBinContent(7,23.16259);
   hmean_Hino_135_Hyy_stack_1->SetBinContent(8,17.92728);
   hmean_Hino_135_Hyy_stack_1->SetBinContent(9,13.292);
   hmean_Hino_135_Hyy_stack_1->SetBinContent(10,8.530674);
   hmean_Hino_135_Hyy_stack_1->SetBinContent(11,3.297803);
   hmean_Hino_135_Hyy_stack_1->SetBinContent(12,-1.423034);
   hmean_Hino_135_Hyy_stack_1->SetBinContent(13,-5.821613);
   hmean_Hino_135_Hyy_stack_1->SetBinContent(14,-11.82986);
   hmean_Hino_135_Hyy_stack_1->SetBinContent(15,-15.68517);
   hmean_Hino_135_Hyy_stack_1->SetBinContent(16,-20.87623);
   hmean_Hino_135_Hyy_stack_1->SetBinContent(17,-25.44413);
   hmean_Hino_135_Hyy_stack_1->SetBinContent(18,-31.8734);
   hmean_Hino_135_Hyy_stack_1->SetBinContent(19,-47.39612);
   hmean_Hino_135_Hyy_stack_1->SetBinContent(20,-62.77481);
   hmean_Hino_135_Hyy_stack_1->SetBinContent(21,-82.30357);
   hmean_Hino_135_Hyy_stack_1->SetBinContent(22,-129.5905);
   hmean_Hino_135_Hyy_stack_1->SetBinError(1,2.813724);
   hmean_Hino_135_Hyy_stack_1->SetBinError(2,2.093819);
   hmean_Hino_135_Hyy_stack_1->SetBinError(3,1.992733);
   hmean_Hino_135_Hyy_stack_1->SetBinError(4,1.745631);
   hmean_Hino_135_Hyy_stack_1->SetBinError(5,1.020037);
   hmean_Hino_135_Hyy_stack_1->SetBinError(6,1.160713);
   hmean_Hino_135_Hyy_stack_1->SetBinError(7,0.779781);
   hmean_Hino_135_Hyy_stack_1->SetBinError(8,0.6076049);
   hmean_Hino_135_Hyy_stack_1->SetBinError(9,0.4645832);
   hmean_Hino_135_Hyy_stack_1->SetBinError(10,0.3203643);
   hmean_Hino_135_Hyy_stack_1->SetBinError(11,0.2731525);
   hmean_Hino_135_Hyy_stack_1->SetBinError(12,0.2838026);
   hmean_Hino_135_Hyy_stack_1->SetBinError(13,0.3658999);
   hmean_Hino_135_Hyy_stack_1->SetBinError(14,0.4043272);
   hmean_Hino_135_Hyy_stack_1->SetBinError(15,0.6278621);
   hmean_Hino_135_Hyy_stack_1->SetBinError(16,0.8637497);
   hmean_Hino_135_Hyy_stack_1->SetBinError(17,0.9940882);
   hmean_Hino_135_Hyy_stack_1->SetBinError(18,1.052372);
   hmean_Hino_135_Hyy_stack_1->SetBinError(19,1.675);
   hmean_Hino_135_Hyy_stack_1->SetBinError(20,2.144946);
   hmean_Hino_135_Hyy_stack_1->SetBinError(21,2.155177);
   hmean_Hino_135_Hyy_stack_1->SetBinError(22,2.770614);
   hmean_Hino_135_Hyy_stack_1->SetEntries(22);

   ci = TColor::GetColor("#ff00ff");
   hmean_Hino_135_Hyy_stack_1->SetLineColor(ci);

   ci = TColor::GetColor("#ff00ff");
   hmean_Hino_135_Hyy_stack_1->SetMarkerColor(ci);
   hmean_Hino_135_Hyy_stack_1->SetMarkerStyle(26);
   hmean_Hino_135_Hyy_stack_1->SetMarkerSize(2);
   hmean_Hino_135_Hyy_stack_1->GetXaxis()->SetLabelFont(42);
   hmean_Hino_135_Hyy_stack_1->GetXaxis()->SetTitleOffset(1);
   hmean_Hino_135_Hyy_stack_1->GetXaxis()->SetTitleFont(42);
   hmean_Hino_135_Hyy_stack_1->GetYaxis()->SetLabelFont(42);
   hmean_Hino_135_Hyy_stack_1->GetYaxis()->SetTitleFont(42);
   hmean_Hino_135_Hyy_stack_1->GetZaxis()->SetLabelFont(42);
   hmean_Hino_135_Hyy_stack_1->GetZaxis()->SetTitleOffset(1);
   hmean_Hino_135_Hyy_stack_1->GetZaxis()->SetTitleFont(42);
   hs_m->Add(hmean_Hino_135_Hyy_stack_1,"");
   Double_t xAxis7[23] = {-700, -400, -300, -250, -200, -150, -125, -100, -75, -50, -25, 0, 25, 50, 75, 100, 125, 150, 200, 250, 300, 400, 700}; 
   
   TH1F *hmean_Hino_135_Zee_stack_2 = new TH1F("hmean_Hino_135_Zee_stack_2","",22, xAxis7);
   hmean_Hino_135_Zee_stack_2->SetBinContent(1,136.7319);
   hmean_Hino_135_Zee_stack_2->SetBinContent(2,86.74226);
   hmean_Hino_135_Zee_stack_2->SetBinContent(3,62.9591);
   hmean_Hino_135_Zee_stack_2->SetBinContent(4,51.51276);
   hmean_Hino_135_Zee_stack_2->SetBinContent(5,36.39552);
   hmean_Hino_135_Zee_stack_2->SetBinContent(6,26.92623);
   hmean_Hino_135_Zee_stack_2->SetBinContent(7,22.53745);
   hmean_Hino_135_Zee_stack_2->SetBinContent(8,18.39727);
   hmean_Hino_135_Zee_stack_2->SetBinContent(9,14.81699);
   hmean_Hino_135_Zee_stack_2->SetBinContent(10,8.983024);
   hmean_Hino_135_Zee_stack_2->SetBinContent(11,4.47338);
   hmean_Hino_135_Zee_stack_2->SetBinContent(12,-1.107453);
   hmean_Hino_135_Zee_stack_2->SetBinContent(13,-7.12547);
   hmean_Hino_135_Zee_stack_2->SetBinContent(14,-11.31614);
   hmean_Hino_135_Zee_stack_2->SetBinContent(15,-17.09657);
   hmean_Hino_135_Zee_stack_2->SetBinContent(16,-21.82865);
   hmean_Hino_135_Zee_stack_2->SetBinContent(17,-26.48529);
   hmean_Hino_135_Zee_stack_2->SetBinContent(18,-34.75954);
   hmean_Hino_135_Zee_stack_2->SetBinContent(19,-43.361);
   hmean_Hino_135_Zee_stack_2->SetBinContent(20,-60.72621);
   hmean_Hino_135_Zee_stack_2->SetBinContent(21,-77.22147);
   hmean_Hino_135_Zee_stack_2->SetBinContent(22,-124.4573);
   hmean_Hino_135_Zee_stack_2->SetBinError(1,2.610443);
   hmean_Hino_135_Zee_stack_2->SetBinError(2,1.8621);
   hmean_Hino_135_Zee_stack_2->SetBinError(3,2.156815);
   hmean_Hino_135_Zee_stack_2->SetBinError(4,1.415511);
   hmean_Hino_135_Zee_stack_2->SetBinError(5,1.035622);
   hmean_Hino_135_Zee_stack_2->SetBinError(6,1.022248);
   hmean_Hino_135_Zee_stack_2->SetBinError(7,0.8757718);
   hmean_Hino_135_Zee_stack_2->SetBinError(8,0.6048418);
   hmean_Hino_135_Zee_stack_2->SetBinError(9,0.4810902);
   hmean_Hino_135_Zee_stack_2->SetBinError(10,0.3826401);
   hmean_Hino_135_Zee_stack_2->SetBinError(11,0.3783161);
   hmean_Hino_135_Zee_stack_2->SetBinError(12,0.3389598);
   hmean_Hino_135_Zee_stack_2->SetBinError(13,0.3857252);
   hmean_Hino_135_Zee_stack_2->SetBinError(14,0.4786068);
   hmean_Hino_135_Zee_stack_2->SetBinError(15,0.7320178);
   hmean_Hino_135_Zee_stack_2->SetBinError(16,0.8056339);
   hmean_Hino_135_Zee_stack_2->SetBinError(17,1.147328);
   hmean_Hino_135_Zee_stack_2->SetBinError(18,1.017741);
   hmean_Hino_135_Zee_stack_2->SetBinError(19,1.478474);
   hmean_Hino_135_Zee_stack_2->SetBinError(20,2.060233);
   hmean_Hino_135_Zee_stack_2->SetBinError(21,1.822655);
   hmean_Hino_135_Zee_stack_2->SetBinError(22,2.688695);
   hmean_Hino_135_Zee_stack_2->SetEntries(22);

   ci = TColor::GetColor("#0000ff");
   hmean_Hino_135_Zee_stack_2->SetLineColor(ci);

   ci = TColor::GetColor("#0000ff");
   hmean_Hino_135_Zee_stack_2->SetMarkerColor(ci);
   hmean_Hino_135_Zee_stack_2->SetMarkerStyle(26);
   hmean_Hino_135_Zee_stack_2->SetMarkerSize(2);
   hmean_Hino_135_Zee_stack_2->GetXaxis()->SetLabelFont(42);
   hmean_Hino_135_Zee_stack_2->GetXaxis()->SetTitleOffset(1);
   hmean_Hino_135_Zee_stack_2->GetXaxis()->SetTitleFont(42);
   hmean_Hino_135_Zee_stack_2->GetYaxis()->SetLabelFont(42);
   hmean_Hino_135_Zee_stack_2->GetYaxis()->SetTitleFont(42);
   hmean_Hino_135_Zee_stack_2->GetZaxis()->SetLabelFont(42);
   hmean_Hino_135_Zee_stack_2->GetZaxis()->SetTitleOffset(1);
   hmean_Hino_135_Zee_stack_2->GetZaxis()->SetTitleFont(42);
   hs_m->Add(hmean_Hino_135_Zee_stack_2,"");
   Double_t xAxis8[23] = {-700, -400, -300, -250, -200, -150, -125, -100, -75, -50, -25, 0, 25, 50, 75, 100, 125, 150, 200, 250, 300, 400, 700}; 
   
   TH1F *hmean_Zee_stack_3 = new TH1F("hmean_Zee_stack_3","",22, xAxis8);
   hmean_Zee_stack_3->SetBinContent(5,26.2472);
   hmean_Zee_stack_3->SetBinContent(6,28.63855);
   hmean_Zee_stack_3->SetBinContent(7,22.41809);
   hmean_Zee_stack_3->SetBinContent(8,18.74179);
   hmean_Zee_stack_3->SetBinContent(9,14.1796);
   hmean_Zee_stack_3->SetBinContent(10,9.053682);
   hmean_Zee_stack_3->SetBinContent(11,4.319993);
   hmean_Zee_stack_3->SetBinContent(12,-0.9074668);
   hmean_Zee_stack_3->SetBinContent(13,-6.027654);
   hmean_Zee_stack_3->SetBinContent(14,-10.61783);
   hmean_Zee_stack_3->SetBinContent(15,-15.37624);
   hmean_Zee_stack_3->SetBinContent(16,-20.17956);
   hmean_Zee_stack_3->SetBinContent(17,-25.93547);
   hmean_Zee_stack_3->SetBinContent(18,-31.49138);
   hmean_Zee_stack_3->SetBinError(5,11.1769);
   hmean_Zee_stack_3->SetBinError(6,2.265167);
   hmean_Zee_stack_3->SetBinError(7,0.5977097);
   hmean_Zee_stack_3->SetBinError(8,0.2405699);
   hmean_Zee_stack_3->SetBinError(9,0.1288551);
   hmean_Zee_stack_3->SetBinError(10,0.07493018);
   hmean_Zee_stack_3->SetBinError(11,0.06582834);
   hmean_Zee_stack_3->SetBinError(12,0.06207859);
   hmean_Zee_stack_3->SetBinError(13,0.08085803);
   hmean_Zee_stack_3->SetBinError(14,0.1289176);
   hmean_Zee_stack_3->SetBinError(15,0.2325543);
   hmean_Zee_stack_3->SetBinError(16,0.5768661);
   hmean_Zee_stack_3->SetBinError(17,1.781097);
   hmean_Zee_stack_3->SetBinError(18,6.682741);
   hmean_Zee_stack_3->SetEntries(15);

   ci = TColor::GetColor("#ff0000");
   hmean_Zee_stack_3->SetLineColor(ci);

   ci = TColor::GetColor("#ff0000");
   hmean_Zee_stack_3->SetMarkerColor(ci);
   hmean_Zee_stack_3->SetMarkerStyle(24);
   hmean_Zee_stack_3->SetMarkerSize(2);
   hmean_Zee_stack_3->GetXaxis()->SetLabelFont(42);
   hmean_Zee_stack_3->GetXaxis()->SetTitleOffset(1);
   hmean_Zee_stack_3->GetXaxis()->SetTitleFont(42);
   hmean_Zee_stack_3->GetYaxis()->SetLabelFont(42);
   hmean_Zee_stack_3->GetYaxis()->SetTitleFont(42);
   hmean_Zee_stack_3->GetZaxis()->SetLabelFont(42);
   hmean_Zee_stack_3->GetZaxis()->SetTitleOffset(1);
   hmean_Zee_stack_3->GetZaxis()->SetTitleFont(42);
   hs_m->Add(hmean_Zee_stack_3,"");
   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_Hyy_stack_1","135GeV Hyy 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_Zee_stack_2","135GeV Zee 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_Zee_stack_3","Zee MC","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);
   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);
}