void PlotPointingBias()
{
//=========Macro generated from canvas: c_m/c
//=========  (Mon May 23 10:28:35 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_Hyy_stack_1 = new TH1F("hmean_Hino_135_Hyy_stack_1","",22, xAxis7);
   hmean_Hino_135_Hyy_stack_1->SetBinContent(1,170.6705);
   hmean_Hino_135_Hyy_stack_1->SetBinContent(2,104.3045);
   hmean_Hino_135_Hyy_stack_1->SetBinContent(3,70.30111);
   hmean_Hino_135_Hyy_stack_1->SetBinContent(4,56.75356);
   hmean_Hino_135_Hyy_stack_1->SetBinContent(5,39.49113);
   hmean_Hino_135_Hyy_stack_1->SetBinContent(6,30.55589);
   hmean_Hino_135_Hyy_stack_1->SetBinContent(7,25.43349);
   hmean_Hino_135_Hyy_stack_1->SetBinContent(8,19.24012);
   hmean_Hino_135_Hyy_stack_1->SetBinContent(9,14.38252);
   hmean_Hino_135_Hyy_stack_1->SetBinContent(10,8.80328);
   hmean_Hino_135_Hyy_stack_1->SetBinContent(11,3.471184);
   hmean_Hino_135_Hyy_stack_1->SetBinContent(12,-1.501376);
   hmean_Hino_135_Hyy_stack_1->SetBinContent(13,-6.014585);
   hmean_Hino_135_Hyy_stack_1->SetBinContent(14,-12.25136);
   hmean_Hino_135_Hyy_stack_1->SetBinContent(15,-16.55743);
   hmean_Hino_135_Hyy_stack_1->SetBinContent(16,-22.00263);
   hmean_Hino_135_Hyy_stack_1->SetBinContent(17,-29.07327);
   hmean_Hino_135_Hyy_stack_1->SetBinContent(18,-35.29497);
   hmean_Hino_135_Hyy_stack_1->SetBinContent(19,-51.68882);
   hmean_Hino_135_Hyy_stack_1->SetBinContent(20,-71.81535);
   hmean_Hino_135_Hyy_stack_1->SetBinContent(21,-95.71429);
   hmean_Hino_135_Hyy_stack_1->SetBinContent(22,-160.7384);
   hmean_Hino_135_Hyy_stack_1->SetBinError(1,3.331125);
   hmean_Hino_135_Hyy_stack_1->SetBinError(2,2.182423);
   hmean_Hino_135_Hyy_stack_1->SetBinError(3,2.060911);
   hmean_Hino_135_Hyy_stack_1->SetBinError(4,1.637876);
   hmean_Hino_135_Hyy_stack_1->SetBinError(5,1.083342);
   hmean_Hino_135_Hyy_stack_1->SetBinError(6,1.243525);
   hmean_Hino_135_Hyy_stack_1->SetBinError(7,0.8342272);
   hmean_Hino_135_Hyy_stack_1->SetBinError(8,0.6630896);
   hmean_Hino_135_Hyy_stack_1->SetBinError(9,0.4466271);
   hmean_Hino_135_Hyy_stack_1->SetBinError(10,0.3443177);
   hmean_Hino_135_Hyy_stack_1->SetBinError(11,0.2915438);
   hmean_Hino_135_Hyy_stack_1->SetBinError(12,0.3098544);
   hmean_Hino_135_Hyy_stack_1->SetBinError(13,0.3714651);
   hmean_Hino_135_Hyy_stack_1->SetBinError(14,0.4172754);
   hmean_Hino_135_Hyy_stack_1->SetBinError(15,0.6201798);
   hmean_Hino_135_Hyy_stack_1->SetBinError(16,0.826638);
   hmean_Hino_135_Hyy_stack_1->SetBinError(17,1.072709);
   hmean_Hino_135_Hyy_stack_1->SetBinError(18,1.020681);
   hmean_Hino_135_Hyy_stack_1->SetBinError(19,1.604808);
   hmean_Hino_135_Hyy_stack_1->SetBinError(20,1.978211);
   hmean_Hino_135_Hyy_stack_1->SetBinError(21,2.278204);
   hmean_Hino_135_Hyy_stack_1->SetBinError(22,2.964532);
   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 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_Zee_stack_2 = new TH1F("hmean_Hino_135_Zee_stack_2","",22, xAxis8);
   hmean_Hino_135_Zee_stack_2->SetBinContent(1,156.0689);
   hmean_Hino_135_Zee_stack_2->SetBinContent(2,93.99088);
   hmean_Hino_135_Zee_stack_2->SetBinContent(3,68.90588);
   hmean_Hino_135_Zee_stack_2->SetBinContent(4,54.87856);
   hmean_Hino_135_Zee_stack_2->SetBinContent(5,38.96018);
   hmean_Hino_135_Zee_stack_2->SetBinContent(6,28.85543);
   hmean_Hino_135_Zee_stack_2->SetBinContent(7,24.13808);
   hmean_Hino_135_Zee_stack_2->SetBinContent(8,19.26167);
   hmean_Hino_135_Zee_stack_2->SetBinContent(9,15.15855);
   hmean_Hino_135_Zee_stack_2->SetBinContent(10,9.369331);
   hmean_Hino_135_Zee_stack_2->SetBinContent(11,4.314612);
   hmean_Hino_135_Zee_stack_2->SetBinContent(12,-1.162064);
   hmean_Hino_135_Zee_stack_2->SetBinContent(13,-7.166437);
   hmean_Hino_135_Zee_stack_2->SetBinContent(14,-11.50228);
   hmean_Hino_135_Zee_stack_2->SetBinContent(15,-17.25619);
   hmean_Hino_135_Zee_stack_2->SetBinContent(16,-22.30745);
   hmean_Hino_135_Zee_stack_2->SetBinContent(17,-28.57574);
   hmean_Hino_135_Zee_stack_2->SetBinContent(18,-37.84972);
   hmean_Hino_135_Zee_stack_2->SetBinContent(19,-49.94757);
   hmean_Hino_135_Zee_stack_2->SetBinContent(20,-65.83729);
   hmean_Hino_135_Zee_stack_2->SetBinContent(21,-87.73516);
   hmean_Hino_135_Zee_stack_2->SetBinContent(22,-151.8253);
   hmean_Hino_135_Zee_stack_2->SetBinError(1,2.596382);
   hmean_Hino_135_Zee_stack_2->SetBinError(2,2.001997);
   hmean_Hino_135_Zee_stack_2->SetBinError(3,2.103908);
   hmean_Hino_135_Zee_stack_2->SetBinError(4,1.641478);
   hmean_Hino_135_Zee_stack_2->SetBinError(5,1.085096);
   hmean_Hino_135_Zee_stack_2->SetBinError(6,1.131107);
   hmean_Hino_135_Zee_stack_2->SetBinError(7,0.8429134);
   hmean_Hino_135_Zee_stack_2->SetBinError(8,0.6238338);
   hmean_Hino_135_Zee_stack_2->SetBinError(9,0.4941693);
   hmean_Hino_135_Zee_stack_2->SetBinError(10,0.3908168);
   hmean_Hino_135_Zee_stack_2->SetBinError(11,0.3100371);
   hmean_Hino_135_Zee_stack_2->SetBinError(12,0.3550544);
   hmean_Hino_135_Zee_stack_2->SetBinError(13,0.3794574);
   hmean_Hino_135_Zee_stack_2->SetBinError(14,0.4931897);
   hmean_Hino_135_Zee_stack_2->SetBinError(15,0.6779155);
   hmean_Hino_135_Zee_stack_2->SetBinError(16,0.8058581);
   hmean_Hino_135_Zee_stack_2->SetBinError(17,1.202637);
   hmean_Hino_135_Zee_stack_2->SetBinError(18,1.038054);
   hmean_Hino_135_Zee_stack_2->SetBinError(19,1.58167);
   hmean_Hino_135_Zee_stack_2->SetBinError(20,2.025134);
   hmean_Hino_135_Zee_stack_2->SetBinError(21,2.040831);
   hmean_Hino_135_Zee_stack_2->SetBinError(22,2.822352);
   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 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_Zee_stack_3 = new TH1F("hmean_Zee_stack_3","",22, xAxis9);
   hmean_Zee_stack_3->SetBinContent(6,28.70661);
   hmean_Zee_stack_3->SetBinContent(7,22.97734);
   hmean_Zee_stack_3->SetBinContent(8,18.80374);
   hmean_Zee_stack_3->SetBinContent(9,14.32856);
   hmean_Zee_stack_3->SetBinContent(10,9.112808);
   hmean_Zee_stack_3->SetBinContent(11,4.147033);
   hmean_Zee_stack_3->SetBinContent(12,-1.190624);
   hmean_Zee_stack_3->SetBinContent(13,-6.012897);
   hmean_Zee_stack_3->SetBinContent(14,-10.84521);
   hmean_Zee_stack_3->SetBinContent(15,-15.54524);
   hmean_Zee_stack_3->SetBinContent(16,-20.21454);
   hmean_Zee_stack_3->SetBinContent(17,-26.7359);
   hmean_Zee_stack_3->SetBinError(6,2.724128);
   hmean_Zee_stack_3->SetBinError(7,0.6084954);
   hmean_Zee_stack_3->SetBinError(8,0.2534539);
   hmean_Zee_stack_3->SetBinError(9,0.1280938);
   hmean_Zee_stack_3->SetBinError(10,0.07906213);
   hmean_Zee_stack_3->SetBinError(11,0.06541807);
   hmean_Zee_stack_3->SetBinError(12,0.06353211);
   hmean_Zee_stack_3->SetBinError(13,0.07900152);
   hmean_Zee_stack_3->SetBinError(14,0.1317429);
   hmean_Zee_stack_3->SetBinError(15,0.2397296);
   hmean_Zee_stack_3->SetBinError(16,0.5658332);
   hmean_Zee_stack_3->SetBinError(17,1.894437);
   hmean_Zee_stack_3->SetEntries(12);

   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,"");
   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_data_all_stack_4 = new TH1F("hmean_data_all_stack_4","",22, xAxis10);
   hmean_data_all_stack_4->SetBinContent(6,28.36755);
   hmean_data_all_stack_4->SetBinContent(7,23.51038);
   hmean_data_all_stack_4->SetBinContent(8,18.85677);
   hmean_data_all_stack_4->SetBinContent(9,14.27549);
   hmean_data_all_stack_4->SetBinContent(10,9.400945);
   hmean_data_all_stack_4->SetBinContent(11,4.25146);
   hmean_data_all_stack_4->SetBinContent(12,-1.209663);
   hmean_data_all_stack_4->SetBinContent(13,-6.381004);
   hmean_data_all_stack_4->SetBinContent(14,-11.3866);
   hmean_data_all_stack_4->SetBinContent(15,-15.69635);
   hmean_data_all_stack_4->SetBinContent(16,-20.94263);
   hmean_data_all_stack_4->SetBinContent(17,-25.08592);
   hmean_data_all_stack_4->SetBinError(6,1.317145);
   hmean_data_all_stack_4->SetBinError(7,0.3240914);
   hmean_data_all_stack_4->SetBinError(8,0.110384);
   hmean_data_all_stack_4->SetBinError(9,0.04478233);
   hmean_data_all_stack_4->SetBinError(10,0.0256005);
   hmean_data_all_stack_4->SetBinError(11,0.01978529);
   hmean_data_all_stack_4->SetBinError(12,0.02042837);
   hmean_data_all_stack_4->SetBinError(13,0.02740584);
   hmean_data_all_stack_4->SetBinError(14,0.05156872);
   hmean_data_all_stack_4->SetBinError(15,0.1191345);
   hmean_data_all_stack_4->SetBinError(16,0.3638064);
   hmean_data_all_stack_4->SetBinError(17,1.461523);
   hmean_data_all_stack_4->SetEntries(12);
   hmean_data_all_stack_4->SetMarkerStyle(20);
   hmean_data_all_stack_4->SetMarkerSize(2);
   hmean_data_all_stack_4->GetXaxis()->SetLabelFont(42);
   hmean_data_all_stack_4->GetXaxis()->SetTitleOffset(1);
   hmean_data_all_stack_4->GetXaxis()->SetTitleFont(42);
   hmean_data_all_stack_4->GetYaxis()->SetLabelFont(42);
   hmean_data_all_stack_4->GetYaxis()->SetTitleFont(42);
   hmean_data_all_stack_4->GetZaxis()->SetLabelFont(42);
   hmean_data_all_stack_4->GetZaxis()->SetTitleOffset(1);
   hmean_data_all_stack_4->GetZaxis()->SetTitleFont(42);
   hs_m->Add(hmean_data_all_stack_4,"");
   hs_m->Draw("nostack p");
   
   TLegend *leg = new TLegend(0.65,0.6,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);
   entry=leg->AddEntry("hmean_data_all_stack_4","Zee data","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);
}