void PlotPointingBiasVrho()
{
//=========Macro generated from canvas: c_m/c
//=========  (Mon May 16 10:01:42 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(-250,-250,2250,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[6] = {0, 80, 160, 300, 520, 2000}; 
   
   TH1F *hs_m_stack_2 = new TH1F("hs_m_stack_2","",5, 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[6] = {0, 80, 160, 300, 520, 2000}; 
   
   TH1F *hmean_Hino_135_10ns_stack_1 = new TH1F("hmean_Hino_135_10ns_stack_1","",5, xAxis7);
   hmean_Hino_135_10ns_stack_1->SetBinContent(1,0.8173263);
   hmean_Hino_135_10ns_stack_1->SetBinContent(2,1.850752);
   hmean_Hino_135_10ns_stack_1->SetBinContent(3,1.302153);
   hmean_Hino_135_10ns_stack_1->SetBinContent(4,-0.04461334);
   hmean_Hino_135_10ns_stack_1->SetBinContent(5,1.536736);
   hmean_Hino_135_10ns_stack_1->SetBinError(1,0.2042759);
   hmean_Hino_135_10ns_stack_1->SetBinError(2,0.199898);
   hmean_Hino_135_10ns_stack_1->SetBinError(3,0.2472084);
   hmean_Hino_135_10ns_stack_1->SetBinError(4,0.3263111);
   hmean_Hino_135_10ns_stack_1->SetBinError(5,0.5024515);
   hmean_Hino_135_10ns_stack_1->SetEntries(5);

   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[6] = {0, 80, 160, 300, 520, 2000}; 
   
   TH1F *hmean_Hino_135_2ns_stack_2 = new TH1F("hmean_Hino_135_2ns_stack_2","",5, xAxis8);
   hmean_Hino_135_2ns_stack_2->SetBinContent(1,0.9946722);
   hmean_Hino_135_2ns_stack_2->SetBinContent(2,1.806061);
   hmean_Hino_135_2ns_stack_2->SetBinContent(3,1.075706);
   hmean_Hino_135_2ns_stack_2->SetBinContent(4,1.086931);
   hmean_Hino_135_2ns_stack_2->SetBinContent(5,0.7431733);
   hmean_Hino_135_2ns_stack_2->SetBinError(1,0.0998381);
   hmean_Hino_135_2ns_stack_2->SetBinError(2,0.09722299);
   hmean_Hino_135_2ns_stack_2->SetBinError(3,0.1227678);
   hmean_Hino_135_2ns_stack_2->SetBinError(4,0.176088);
   hmean_Hino_135_2ns_stack_2->SetBinError(5,0.2846239);
   hmean_Hino_135_2ns_stack_2->SetEntries(5);

   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[6] = {0, 80, 160, 300, 520, 2000}; 
   
   TH1F *hmean_Hino_135_Hyy_stack_3 = new TH1F("hmean_Hino_135_Hyy_stack_3","",5, xAxis9);
   hmean_Hino_135_Hyy_stack_3->SetBinContent(1,0.9771867);
   hmean_Hino_135_Hyy_stack_3->SetBinContent(2,2.25213);
   hmean_Hino_135_Hyy_stack_3->SetBinContent(3,0.7182798);
   hmean_Hino_135_Hyy_stack_3->SetBinContent(4,1.32906);
   hmean_Hino_135_Hyy_stack_3->SetBinContent(5,0.3770981);
   hmean_Hino_135_Hyy_stack_3->SetBinError(1,0.09640732);
   hmean_Hino_135_Hyy_stack_3->SetBinError(2,0.1088653);
   hmean_Hino_135_Hyy_stack_3->SetBinError(3,0.14698);
   hmean_Hino_135_Hyy_stack_3->SetBinError(4,0.2146725);
   hmean_Hino_135_Hyy_stack_3->SetBinError(5,0.3025756);
   hmean_Hino_135_Hyy_stack_3->SetEntries(5);

   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[6] = {0, 80, 160, 300, 520, 2000}; 
   
   TH1F *hmean_Zee_stack_4 = new TH1F("hmean_Zee_stack_4","",5, xAxis10);
   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);
}