void PlotPointingBiasVrho()
{
//=========Macro generated from canvas: c_m/c
//=========  (Thu Jun  2 10:10:59 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 Resolution [mm];");
   hs_m->SetMinimum(-200);
   hs_m->SetMaximum(200);
   Double_t xAxis1[10] = {0, 20, 40, 60, 80, 100, 120, 300, 500, 2000}; 
   
   TH1F *hs_m_stack_1 = new TH1F("hs_m_stack_1","",9, xAxis1);
   hs_m_stack_1->SetMinimum(-200);
   hs_m_stack_1->SetMaximum(200);
   hs_m_stack_1->SetDirectory(0);
   hs_m_stack_1->SetStats(0);

   Int_t ci;      // for color index setting
   TColor *color; // for color definition with alpha
   ci = TColor::GetColor("#000099");
   hs_m_stack_1->SetLineColor(ci);
   hs_m_stack_1->GetXaxis()->SetTitle("|z_{origin}| [mm]");
   hs_m_stack_1->GetXaxis()->SetLabelFont(42);
   hs_m_stack_1->GetXaxis()->SetTitleOffset(1);
   hs_m_stack_1->GetXaxis()->SetTitleFont(42);
   hs_m_stack_1->GetYaxis()->SetTitle(" Pointing Resolution [mm]");
   hs_m_stack_1->GetYaxis()->SetLabelFont(42);
   hs_m_stack_1->GetYaxis()->SetTitleFont(42);
   hs_m_stack_1->GetZaxis()->SetLabelFont(42);
   hs_m_stack_1->GetZaxis()->SetTitleOffset(1);
   hs_m_stack_1->GetZaxis()->SetTitleFont(42);
   hs_m->SetHistogram(hs_m_stack_1);
   
   Double_t xAxis2[10] = {0, 20, 40, 60, 80, 100, 120, 300, 500, 2000}; 
   
   TH1F *hmean_Hino_2ns_stack_1 = new TH1F("hmean_Hino_2ns_stack_1","",9, xAxis2);
   hmean_Hino_2ns_stack_1->SetBinContent(1,-1.774653);
   hmean_Hino_2ns_stack_1->SetBinContent(2,-6.6545);
   hmean_Hino_2ns_stack_1->SetBinContent(3,-12.70475);
   hmean_Hino_2ns_stack_1->SetBinContent(4,-15.75124);
   hmean_Hino_2ns_stack_1->SetBinContent(5,-21.44527);
   hmean_Hino_2ns_stack_1->SetBinContent(6,-27.0001);
   hmean_Hino_2ns_stack_1->SetBinContent(7,-49.76247);
   hmean_Hino_2ns_stack_1->SetBinContent(8,-130.2184);
   hmean_Hino_2ns_stack_1->SetBinContent(9,-256.5955);
   hmean_Hino_2ns_stack_1->SetBinError(1,0.03553059);
   hmean_Hino_2ns_stack_1->SetBinError(2,0.04189272);
   hmean_Hino_2ns_stack_1->SetBinError(3,0.05024572);
   hmean_Hino_2ns_stack_1->SetBinError(4,0.06155321);
   hmean_Hino_2ns_stack_1->SetBinError(5,0.07036816);
   hmean_Hino_2ns_stack_1->SetBinError(6,0.08981939);
   hmean_Hino_2ns_stack_1->SetBinError(7,0.06379192);
   hmean_Hino_2ns_stack_1->SetBinError(8,0.1494207);
   hmean_Hino_2ns_stack_1->SetBinError(9,1.128038);
   hmean_Hino_2ns_stack_1->SetEntries(9);

   ci = TColor::GetColor("#00ff00");
   hmean_Hino_2ns_stack_1->SetLineColor(ci);

   ci = TColor::GetColor("#00ff00");
   hmean_Hino_2ns_stack_1->SetMarkerColor(ci);
   hmean_Hino_2ns_stack_1->SetMarkerStyle(26);
   hmean_Hino_2ns_stack_1->SetMarkerSize(2);
   hmean_Hino_2ns_stack_1->GetXaxis()->SetLabelFont(42);
   hmean_Hino_2ns_stack_1->GetXaxis()->SetTitleOffset(1);
   hmean_Hino_2ns_stack_1->GetXaxis()->SetTitleFont(42);
   hmean_Hino_2ns_stack_1->GetYaxis()->SetLabelFont(42);
   hmean_Hino_2ns_stack_1->GetYaxis()->SetTitleFont(42);
   hmean_Hino_2ns_stack_1->GetZaxis()->SetLabelFont(42);
   hmean_Hino_2ns_stack_1->GetZaxis()->SetTitleOffset(1);
   hmean_Hino_2ns_stack_1->GetZaxis()->SetTitleFont(42);
   hs_m->Add(hmean_Hino_2ns_stack_1,"");
   Double_t xAxis3[10] = {0, 20, 40, 60, 80, 100, 120, 300, 500, 2000}; 
   
   TH1F *hmean_Hino_325_Hyy_stack_2 = new TH1F("hmean_Hino_325_Hyy_stack_2","",9, xAxis3);
   hmean_Hino_325_Hyy_stack_2->SetBinContent(1,-2.80478);
   hmean_Hino_325_Hyy_stack_2->SetBinContent(2,-5.994893);
   hmean_Hino_325_Hyy_stack_2->SetBinContent(3,-10.82783);
   hmean_Hino_325_Hyy_stack_2->SetBinContent(4,-13.98789);
   hmean_Hino_325_Hyy_stack_2->SetBinContent(5,-17.43292);
   hmean_Hino_325_Hyy_stack_2->SetBinContent(6,-24.94314);
   hmean_Hino_325_Hyy_stack_2->SetBinContent(7,-45.65496);
   hmean_Hino_325_Hyy_stack_2->SetBinContent(8,-121.5579);
   hmean_Hino_325_Hyy_stack_2->SetBinContent(9,-255.0621);
   hmean_Hino_325_Hyy_stack_2->SetBinError(1,0.168328);
   hmean_Hino_325_Hyy_stack_2->SetBinError(2,0.1560754);
   hmean_Hino_325_Hyy_stack_2->SetBinError(3,0.2050919);
   hmean_Hino_325_Hyy_stack_2->SetBinError(4,0.2663996);
   hmean_Hino_325_Hyy_stack_2->SetBinError(5,0.3711365);
   hmean_Hino_325_Hyy_stack_2->SetBinError(6,0.4502929);
   hmean_Hino_325_Hyy_stack_2->SetBinError(7,0.3273074);
   hmean_Hino_325_Hyy_stack_2->SetBinError(8,0.9510366);
   hmean_Hino_325_Hyy_stack_2->SetBinError(9,0.7223145);
   hmean_Hino_325_Hyy_stack_2->SetEntries(9);

   ci = TColor::GetColor("#ffcc00");
   hmean_Hino_325_Hyy_stack_2->SetLineColor(ci);

   ci = TColor::GetColor("#ffcc00");
   hmean_Hino_325_Hyy_stack_2->SetMarkerColor(ci);
   hmean_Hino_325_Hyy_stack_2->SetMarkerStyle(26);
   hmean_Hino_325_Hyy_stack_2->SetMarkerSize(2);
   hmean_Hino_325_Hyy_stack_2->GetXaxis()->SetLabelFont(42);
   hmean_Hino_325_Hyy_stack_2->GetXaxis()->SetTitleOffset(1);
   hmean_Hino_325_Hyy_stack_2->GetXaxis()->SetTitleFont(42);
   hmean_Hino_325_Hyy_stack_2->GetYaxis()->SetLabelFont(42);
   hmean_Hino_325_Hyy_stack_2->GetYaxis()->SetTitleFont(42);
   hmean_Hino_325_Hyy_stack_2->GetZaxis()->SetLabelFont(42);
   hmean_Hino_325_Hyy_stack_2->GetZaxis()->SetTitleOffset(1);
   hmean_Hino_325_Hyy_stack_2->GetZaxis()->SetTitleFont(42);
   hs_m->Add(hmean_Hino_325_Hyy_stack_2,"");
   Double_t xAxis4[10] = {0, 20, 40, 60, 80, 100, 120, 300, 500, 2000}; 
   
   TH1F *hmean_Hino_135_Zee_stack_3 = new TH1F("hmean_Hino_135_Zee_stack_3","",9, xAxis4);
   hmean_Hino_135_Zee_stack_3->SetBinContent(1,-1.803193);
   hmean_Hino_135_Zee_stack_3->SetBinContent(2,-6.400372);
   hmean_Hino_135_Zee_stack_3->SetBinContent(3,-12.13272);
   hmean_Hino_135_Zee_stack_3->SetBinContent(4,-15.80955);
   hmean_Hino_135_Zee_stack_3->SetBinContent(5,-19.20107);
   hmean_Hino_135_Zee_stack_3->SetBinContent(6,-24.08948);
   hmean_Hino_135_Zee_stack_3->SetBinContent(7,-49.4388);
   hmean_Hino_135_Zee_stack_3->SetBinContent(8,-130.9747);
   hmean_Hino_135_Zee_stack_3->SetBinContent(9,-244.3493);
   hmean_Hino_135_Zee_stack_3->SetBinError(1,0.1121489);
   hmean_Hino_135_Zee_stack_3->SetBinError(2,0.1238994);
   hmean_Hino_135_Zee_stack_3->SetBinError(3,0.1319501);
   hmean_Hino_135_Zee_stack_3->SetBinError(4,0.1985226);
   hmean_Hino_135_Zee_stack_3->SetBinError(5,0.2034695);
   hmean_Hino_135_Zee_stack_3->SetBinError(6,0.2683608);
   hmean_Hino_135_Zee_stack_3->SetBinError(7,0.1850917);
   hmean_Hino_135_Zee_stack_3->SetBinError(8,0.4531363);
   hmean_Hino_135_Zee_stack_3->SetBinError(9,2.726679);
   hmean_Hino_135_Zee_stack_3->SetEntries(9);

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

   ci = TColor::GetColor("#ff00ff");
   hmean_Hino_135_Zee_stack_3->SetMarkerColor(ci);
   hmean_Hino_135_Zee_stack_3->SetMarkerStyle(26);
   hmean_Hino_135_Zee_stack_3->SetMarkerSize(2);
   hmean_Hino_135_Zee_stack_3->GetXaxis()->SetLabelFont(42);
   hmean_Hino_135_Zee_stack_3->GetXaxis()->SetTitleOffset(1);
   hmean_Hino_135_Zee_stack_3->GetXaxis()->SetTitleFont(42);
   hmean_Hino_135_Zee_stack_3->GetYaxis()->SetLabelFont(42);
   hmean_Hino_135_Zee_stack_3->GetYaxis()->SetTitleFont(42);
   hmean_Hino_135_Zee_stack_3->GetZaxis()->SetLabelFont(42);
   hmean_Hino_135_Zee_stack_3->GetZaxis()->SetTitleOffset(1);
   hmean_Hino_135_Zee_stack_3->GetZaxis()->SetTitleFont(42);
   hs_m->Add(hmean_Hino_135_Zee_stack_3,"");
   Double_t xAxis5[10] = {0, 20, 40, 60, 80, 100, 120, 300, 500, 2000}; 
   
   TH1F *hmean_Hino_135_2ns_stack_4 = new TH1F("hmean_Hino_135_2ns_stack_4","",9, xAxis5);
   hmean_Hino_135_2ns_stack_4->SetBinContent(1,-1.309783);
   hmean_Hino_135_2ns_stack_4->SetBinContent(2,-6.934951);
   hmean_Hino_135_2ns_stack_4->SetBinContent(3,-13.78039);
   hmean_Hino_135_2ns_stack_4->SetBinContent(4,-16.25996);
   hmean_Hino_135_2ns_stack_4->SetBinContent(5,-21.13138);
   hmean_Hino_135_2ns_stack_4->SetBinContent(6,-25.48223);
   hmean_Hino_135_2ns_stack_4->SetBinContent(7,-50.7233);
   hmean_Hino_135_2ns_stack_4->SetBinContent(8,-132.8711);
   hmean_Hino_135_2ns_stack_4->SetBinContent(9,-252.8059);
   hmean_Hino_135_2ns_stack_4->SetBinError(1,0.09502247);
   hmean_Hino_135_2ns_stack_4->SetBinError(2,0.1181133);
   hmean_Hino_135_2ns_stack_4->SetBinError(3,0.1357861);
   hmean_Hino_135_2ns_stack_4->SetBinError(4,0.1622826);
   hmean_Hino_135_2ns_stack_4->SetBinError(5,0.1828527);
   hmean_Hino_135_2ns_stack_4->SetBinError(6,0.2363051);
   hmean_Hino_135_2ns_stack_4->SetBinError(7,0.1598374);
   hmean_Hino_135_2ns_stack_4->SetBinError(8,0.369418);
   hmean_Hino_135_2ns_stack_4->SetBinError(9,3.374964);
   hmean_Hino_135_2ns_stack_4->SetEntries(9);

   ci = TColor::GetColor("#cc00ff");
   hmean_Hino_135_2ns_stack_4->SetLineColor(ci);

   ci = TColor::GetColor("#cc00ff");
   hmean_Hino_135_2ns_stack_4->SetMarkerColor(ci);
   hmean_Hino_135_2ns_stack_4->SetMarkerStyle(26);
   hmean_Hino_135_2ns_stack_4->SetMarkerSize(2);
   hmean_Hino_135_2ns_stack_4->GetXaxis()->SetLabelFont(42);
   hmean_Hino_135_2ns_stack_4->GetXaxis()->SetTitleOffset(1);
   hmean_Hino_135_2ns_stack_4->GetXaxis()->SetTitleFont(42);
   hmean_Hino_135_2ns_stack_4->GetYaxis()->SetLabelFont(42);
   hmean_Hino_135_2ns_stack_4->GetYaxis()->SetTitleFont(42);
   hmean_Hino_135_2ns_stack_4->GetZaxis()->SetLabelFont(42);
   hmean_Hino_135_2ns_stack_4->GetZaxis()->SetTitleOffset(1);
   hmean_Hino_135_2ns_stack_4->GetZaxis()->SetTitleFont(42);
   hs_m->Add(hmean_Hino_135_2ns_stack_4,"");
   Double_t xAxis6[10] = {0, 20, 40, 60, 80, 100, 120, 300, 500, 2000}; 
   
   TH1F *hmean_data_skimmed_stack_5 = new TH1F("hmean_data_skimmed_stack_5","",9, xAxis6);
   hmean_data_skimmed_stack_5->SetBinContent(1,-2.445387);
   hmean_data_skimmed_stack_5->SetBinContent(2,-7.094371);
   hmean_data_skimmed_stack_5->SetBinContent(3,-11.73411);
   hmean_data_skimmed_stack_5->SetBinContent(4,-16.34245);
   hmean_data_skimmed_stack_5->SetBinContent(5,-20.92836);
   hmean_data_skimmed_stack_5->SetBinContent(6,-25.78485);
   hmean_data_skimmed_stack_5->SetBinContent(7,-30.86056);
   hmean_data_skimmed_stack_5->SetBinError(1,0.005535387);
   hmean_data_skimmed_stack_5->SetBinError(2,0.00689502);
   hmean_data_skimmed_stack_5->SetBinError(3,0.01048189);
   hmean_data_skimmed_stack_5->SetBinError(4,0.01763284);
   hmean_data_skimmed_stack_5->SetBinError(5,0.03836241);
   hmean_data_skimmed_stack_5->SetBinError(6,0.08680015);
   hmean_data_skimmed_stack_5->SetBinError(7,0.2868017);
   hmean_data_skimmed_stack_5->SetEntries(9);
   hmean_data_skimmed_stack_5->SetMarkerStyle(20);
   hmean_data_skimmed_stack_5->SetMarkerSize(2);
   hmean_data_skimmed_stack_5->GetXaxis()->SetLabelFont(42);
   hmean_data_skimmed_stack_5->GetXaxis()->SetTitleOffset(1);
   hmean_data_skimmed_stack_5->GetXaxis()->SetTitleFont(42);
   hmean_data_skimmed_stack_5->GetYaxis()->SetLabelFont(42);
   hmean_data_skimmed_stack_5->GetYaxis()->SetTitleFont(42);
   hmean_data_skimmed_stack_5->GetZaxis()->SetLabelFont(42);
   hmean_data_skimmed_stack_5->GetZaxis()->SetTitleOffset(1);
   hmean_data_skimmed_stack_5->GetZaxis()->SetTitleFont(42);
   hs_m->Add(hmean_data_skimmed_stack_5,"");
   Double_t xAxis7[10] = {0, 20, 40, 60, 80, 100, 120, 300, 500, 2000}; 
   
   TH1F *hmean_Zee_truZ_stack_6 = new TH1F("hmean_Zee_truZ_stack_6","",9, xAxis7);
   hmean_Zee_truZ_stack_6->SetBinContent(1,-2.523327);
   hmean_Zee_truZ_stack_6->SetBinContent(2,-6.668474);
   hmean_Zee_truZ_stack_6->SetBinContent(3,-11.17236);
   hmean_Zee_truZ_stack_6->SetBinContent(4,-16.20149);
   hmean_Zee_truZ_stack_6->SetBinContent(5,-19.17834);
   hmean_Zee_truZ_stack_6->SetBinContent(6,-24.79618);
   hmean_Zee_truZ_stack_6->SetBinContent(7,-32.69033);
   hmean_Zee_truZ_stack_6->SetBinError(1,0.01883723);
   hmean_Zee_truZ_stack_6->SetBinError(2,0.02140598);
   hmean_Zee_truZ_stack_6->SetBinError(3,0.03042137);
   hmean_Zee_truZ_stack_6->SetBinError(4,0.04860204);
   hmean_Zee_truZ_stack_6->SetBinError(5,0.09033302);
   hmean_Zee_truZ_stack_6->SetBinError(6,0.1904913);
   hmean_Zee_truZ_stack_6->SetBinError(7,0.4199159);
   hmean_Zee_truZ_stack_6->SetEntries(9);

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

   ci = TColor::GetColor("#0000ff");
   hmean_Zee_truZ_stack_6->SetMarkerColor(ci);
   hmean_Zee_truZ_stack_6->SetMarkerStyle(20);
   hmean_Zee_truZ_stack_6->SetMarkerSize(2);
   hmean_Zee_truZ_stack_6->GetXaxis()->SetLabelFont(42);
   hmean_Zee_truZ_stack_6->GetXaxis()->SetTitleOffset(1);
   hmean_Zee_truZ_stack_6->GetXaxis()->SetTitleFont(42);
   hmean_Zee_truZ_stack_6->GetYaxis()->SetLabelFont(42);
   hmean_Zee_truZ_stack_6->GetYaxis()->SetTitleFont(42);
   hmean_Zee_truZ_stack_6->GetZaxis()->SetLabelFont(42);
   hmean_Zee_truZ_stack_6->GetZaxis()->SetTitleOffset(1);
   hmean_Zee_truZ_stack_6->GetZaxis()->SetTitleFont(42);
   hs_m->Add(hmean_Zee_truZ_stack_6,"");
   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_2ns_stack_1","Hino_2ns","lpf");
   entry->SetFillStyle(1001);

   ci = TColor::GetColor("#00ff00");
   entry->SetLineColor(ci);
   entry->SetLineStyle(1);
   entry->SetLineWidth(1);

   ci = TColor::GetColor("#00ff00");
   entry->SetMarkerColor(ci);
   entry->SetMarkerStyle(26);
   entry->SetMarkerSize(2);
   entry->SetTextFont(42);
   entry=leg->AddEntry("hmean_Hino_325_Hyy_stack_2","Hino_325_Hyy","lpf");
   entry->SetFillStyle(1001);

   ci = TColor::GetColor("#ffcc00");
   entry->SetLineColor(ci);
   entry->SetLineStyle(1);
   entry->SetLineWidth(1);

   ci = TColor::GetColor("#ffcc00");
   entry->SetMarkerColor(ci);
   entry->SetMarkerStyle(26);
   entry->SetMarkerSize(2);
   entry->SetTextFont(42);
   entry=leg->AddEntry("hmean_Hino_135_Zee_stack_3","Hino_135_Zee","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_4","Hino_135_2ns","lpf");
   entry->SetFillStyle(1001);

   ci = TColor::GetColor("#cc00ff");
   entry->SetLineColor(ci);
   entry->SetLineStyle(1);
   entry->SetLineWidth(1);

   ci = TColor::GetColor("#cc00ff");
   entry->SetMarkerColor(ci);
   entry->SetMarkerStyle(26);
   entry->SetMarkerSize(2);
   entry->SetTextFont(42);
   entry=leg->AddEntry("hmean_data_skimmed_stack_5","data_skimmed","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);
   entry=leg->AddEntry("hmean_Zee_truZ_stack_6","Zee_truZ","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(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);
}