void PlotPointingRes()
{
//=========Macro generated from canvas: c/c
//=========  (Tue May 10 10:04:49 2022) by ROOT version 6.20/06
   TCanvas *c = new TCanvas("c", "c",0,0,1600,1200);
   gStyle->SetOptStat(0);
   c->Range(-87.50001,-18.75,787.5,168.75);
   c->SetFillColor(0);
   c->SetBorderMode(0);
   c->SetBorderSize(2);
   c->SetFrameBorderMode(0);
   c->SetFrameBorderMode(0);
   
   THStack *hs = new THStack();
   hs->SetName("hs");
   hs->SetTitle(";|z_{origin}| [mm]; Pointing Resolution [mm];");
   hs->SetMinimum(0);
   hs->SetMaximum(150);
   Double_t xAxis1[10] = {0, 25, 50, 75, 100, 125, 150, 250, 400, 700}; 
   
   TH1F *hs_stack_1 = new TH1F("hs_stack_1","",9, xAxis1);
   hs_stack_1->SetMinimum(0);
   hs_stack_1->SetMaximum(150);
   hs_stack_1->SetDirectory(0);
   hs_stack_1->SetStats(0);

   Int_t ci;      // for color index setting
   TColor *color; // for color definition with alpha
   ci = TColor::GetColor("#000099");
   hs_stack_1->SetLineColor(ci);
   hs_stack_1->GetXaxis()->SetTitle("|z_{origin}| [mm]");
   hs_stack_1->GetXaxis()->SetLabelFont(42);
   hs_stack_1->GetXaxis()->SetTitleOffset(1);
   hs_stack_1->GetXaxis()->SetTitleFont(42);
   hs_stack_1->GetYaxis()->SetTitle(" Pointing Resolution [mm]");
   hs_stack_1->GetYaxis()->SetLabelFont(42);
   hs_stack_1->GetYaxis()->SetTitleFont(42);
   hs_stack_1->GetZaxis()->SetLabelFont(42);
   hs_stack_1->GetZaxis()->SetTitleOffset(1);
   hs_stack_1->GetZaxis()->SetTitleFont(42);
   hs->SetHistogram(hs_stack_1);
   
   Double_t xAxis2[10] = {0, 25, 50, 75, 100, 125, 150, 250, 400, 700}; 
   
   TH1F *hsig_Hino_135_Hyy_stack_1 = new TH1F("hsig_Hino_135_Hyy_stack_1","",9, xAxis2);
   hsig_Hino_135_Hyy_stack_1->SetBinContent(1,7.913764);
   hsig_Hino_135_Hyy_stack_1->SetBinContent(2,9.138239);
   hsig_Hino_135_Hyy_stack_1->SetBinContent(3,9.627999);
   hsig_Hino_135_Hyy_stack_1->SetBinContent(4,14.5734);
   hsig_Hino_135_Hyy_stack_1->SetBinContent(5,21.93221);
   hsig_Hino_135_Hyy_stack_1->SetBinContent(6,23.89254);
   hsig_Hino_135_Hyy_stack_1->SetBinContent(7,36.56425);
   hsig_Hino_135_Hyy_stack_1->SetBinContent(8,67.53034);
   hsig_Hino_135_Hyy_stack_1->SetBinContent(9,92.77773);
   hsig_Hino_135_Hyy_stack_1->SetBinError(1,0.571364);
   hsig_Hino_135_Hyy_stack_1->SetBinError(2,0.5884143);
   hsig_Hino_135_Hyy_stack_1->SetBinError(3,0.8755084);
   hsig_Hino_135_Hyy_stack_1->SetBinError(4,1.308099);
   hsig_Hino_135_Hyy_stack_1->SetBinError(5,2.61858);
   hsig_Hino_135_Hyy_stack_1->SetBinError(6,1.643088);
   hsig_Hino_135_Hyy_stack_1->SetBinError(7,1.4761);
   hsig_Hino_135_Hyy_stack_1->SetBinError(8,4.891775);
   hsig_Hino_135_Hyy_stack_1->SetBinError(9,7.406205);
   hsig_Hino_135_Hyy_stack_1->SetEntries(9);

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

   ci = TColor::GetColor("#ff00ff");
   hsig_Hino_135_Hyy_stack_1->SetMarkerColor(ci);
   hsig_Hino_135_Hyy_stack_1->SetMarkerStyle(26);
   hsig_Hino_135_Hyy_stack_1->SetMarkerSize(2);
   hsig_Hino_135_Hyy_stack_1->GetXaxis()->SetTitle("|z_{origin}| [mm]");
   hsig_Hino_135_Hyy_stack_1->GetXaxis()->SetLabelFont(42);
   hsig_Hino_135_Hyy_stack_1->GetXaxis()->SetTitleOffset(1);
   hsig_Hino_135_Hyy_stack_1->GetXaxis()->SetTitleFont(42);
   hsig_Hino_135_Hyy_stack_1->GetYaxis()->SetTitle(" Pointing Resolution [mm]");
   hsig_Hino_135_Hyy_stack_1->GetYaxis()->SetLabelFont(42);
   hsig_Hino_135_Hyy_stack_1->GetYaxis()->SetTitleFont(42);
   hsig_Hino_135_Hyy_stack_1->GetZaxis()->SetLabelFont(42);
   hsig_Hino_135_Hyy_stack_1->GetZaxis()->SetTitleOffset(1);
   hsig_Hino_135_Hyy_stack_1->GetZaxis()->SetTitleFont(42);
   hs->Add(hsig_Hino_135_Hyy_stack_1,"");
   Double_t xAxis3[10] = {0, 25, 50, 75, 100, 125, 150, 250, 400, 700}; 
   
   TH1F *hsig_Hino_135_Zee_stack_2 = new TH1F("hsig_Hino_135_Zee_stack_2","",9, xAxis3);
   hsig_Hino_135_Zee_stack_2->SetBinContent(1,7.438476);
   hsig_Hino_135_Zee_stack_2->SetBinContent(2,8.563517);
   hsig_Hino_135_Zee_stack_2->SetBinContent(3,11.756);
   hsig_Hino_135_Zee_stack_2->SetBinContent(4,17.93921);
   hsig_Hino_135_Zee_stack_2->SetBinContent(5,20.64044);
   hsig_Hino_135_Zee_stack_2->SetBinContent(6,21.55551);
   hsig_Hino_135_Zee_stack_2->SetBinContent(7,35.89585);
   hsig_Hino_135_Zee_stack_2->SetBinContent(8,55.80571);
   hsig_Hino_135_Zee_stack_2->SetBinContent(9,85.69232);
   hsig_Hino_135_Zee_stack_2->SetBinError(1,0.7565281);
   hsig_Hino_135_Zee_stack_2->SetBinError(2,0.6656887);
   hsig_Hino_135_Zee_stack_2->SetBinError(3,1.100444);
   hsig_Hino_135_Zee_stack_2->SetBinError(4,1.999108);
   hsig_Hino_135_Zee_stack_2->SetBinError(5,2.875944);
   hsig_Hino_135_Zee_stack_2->SetBinError(6,2.181885);
   hsig_Hino_135_Zee_stack_2->SetBinError(7,2.652559);
   hsig_Hino_135_Zee_stack_2->SetBinError(8,3.732598);
   hsig_Hino_135_Zee_stack_2->SetBinError(9,15.74429);
   hsig_Hino_135_Zee_stack_2->SetEntries(9);

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

   ci = TColor::GetColor("#0000ff");
   hsig_Hino_135_Zee_stack_2->SetMarkerColor(ci);
   hsig_Hino_135_Zee_stack_2->SetMarkerStyle(26);
   hsig_Hino_135_Zee_stack_2->SetMarkerSize(2);
   hsig_Hino_135_Zee_stack_2->GetXaxis()->SetTitle("|z_{origin}| [mm]");
   hsig_Hino_135_Zee_stack_2->GetXaxis()->SetLabelFont(42);
   hsig_Hino_135_Zee_stack_2->GetXaxis()->SetTitleOffset(1);
   hsig_Hino_135_Zee_stack_2->GetXaxis()->SetTitleFont(42);
   hsig_Hino_135_Zee_stack_2->GetYaxis()->SetTitle(" Pointing Resolution [mm]");
   hsig_Hino_135_Zee_stack_2->GetYaxis()->SetLabelFont(42);
   hsig_Hino_135_Zee_stack_2->GetYaxis()->SetTitleFont(42);
   hsig_Hino_135_Zee_stack_2->GetZaxis()->SetLabelFont(42);
   hsig_Hino_135_Zee_stack_2->GetZaxis()->SetTitleOffset(1);
   hsig_Hino_135_Zee_stack_2->GetZaxis()->SetTitleFont(42);
   hs->Add(hsig_Hino_135_Zee_stack_2,"");
   Double_t xAxis4[10] = {0, 25, 50, 75, 100, 125, 150, 250, 400, 700}; 
   
   TH1F *hsig_Zee_stack_3 = new TH1F("hsig_Zee_stack_3","",9, xAxis4);
   hsig_Zee_stack_3->SetBinContent(1,8.287868);
   hsig_Zee_stack_3->SetBinContent(2,8.399764);
   hsig_Zee_stack_3->SetBinContent(3,9.768364);
   hsig_Zee_stack_3->SetBinContent(4,13.20532);
   hsig_Zee_stack_3->SetBinContent(5,16.96514);
   hsig_Zee_stack_3->SetBinContent(6,19.79158);
   hsig_Zee_stack_3->SetBinError(1,0.160437);
   hsig_Zee_stack_3->SetBinError(2,0.1934505);
   hsig_Zee_stack_3->SetBinError(3,0.3436717);
   hsig_Zee_stack_3->SetBinError(4,0.8746367);
   hsig_Zee_stack_3->SetBinError(5,3.174802);
   hsig_Zee_stack_3->SetBinError(6,4.653077);
   hsig_Zee_stack_3->SetEntries(6);

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

   ci = TColor::GetColor("#ff0000");
   hsig_Zee_stack_3->SetMarkerColor(ci);
   hsig_Zee_stack_3->SetMarkerStyle(24);
   hsig_Zee_stack_3->SetMarkerSize(2);
   hsig_Zee_stack_3->GetXaxis()->SetTitle("|z_{origin}| [mm]");
   hsig_Zee_stack_3->GetXaxis()->SetLabelFont(42);
   hsig_Zee_stack_3->GetXaxis()->SetTitleOffset(1);
   hsig_Zee_stack_3->GetXaxis()->SetTitleFont(42);
   hsig_Zee_stack_3->GetYaxis()->SetTitle(" Pointing Resolution [mm]");
   hsig_Zee_stack_3->GetYaxis()->SetLabelFont(42);
   hsig_Zee_stack_3->GetYaxis()->SetTitleFont(42);
   hsig_Zee_stack_3->GetZaxis()->SetLabelFont(42);
   hsig_Zee_stack_3->GetZaxis()->SetTitleOffset(1);
   hsig_Zee_stack_3->GetZaxis()->SetTitleFont(42);
   hs->Add(hsig_Zee_stack_3,"");
   Double_t xAxis5[10] = {0, 25, 50, 75, 100, 125, 150, 250, 400, 700}; 
   
   TH1F *hsig_data_all_stack_4 = new TH1F("hsig_data_all_stack_4","",9, xAxis5);
   hsig_data_all_stack_4->SetBinContent(1,8.64034);
   hsig_data_all_stack_4->SetBinContent(2,8.733828);
   hsig_data_all_stack_4->SetBinContent(3,9.991446);
   hsig_data_all_stack_4->SetBinContent(4,13.36151);
   hsig_data_all_stack_4->SetBinContent(5,14.70785);
   hsig_data_all_stack_4->SetBinContent(6,18.96221);
   hsig_data_all_stack_4->SetBinError(1,0.05157228);
   hsig_data_all_stack_4->SetBinError(2,0.07092891);
   hsig_data_all_stack_4->SetBinError(3,0.1413403);
   hsig_data_all_stack_4->SetBinError(4,0.3769031);
   hsig_data_all_stack_4->SetBinError(5,1.172572);
   hsig_data_all_stack_4->SetBinError(6,3.461457);
   hsig_data_all_stack_4->SetEntries(6);
   hsig_data_all_stack_4->SetMarkerStyle(20);
   hsig_data_all_stack_4->SetMarkerSize(2);
   hsig_data_all_stack_4->GetXaxis()->SetTitle("|z_{origin}| [mm]");
   hsig_data_all_stack_4->GetXaxis()->SetLabelFont(42);
   hsig_data_all_stack_4->GetXaxis()->SetTitleOffset(1);
   hsig_data_all_stack_4->GetXaxis()->SetTitleFont(42);
   hsig_data_all_stack_4->GetYaxis()->SetTitle(" Pointing Resolution [mm]");
   hsig_data_all_stack_4->GetYaxis()->SetLabelFont(42);
   hsig_data_all_stack_4->GetYaxis()->SetTitleFont(42);
   hsig_data_all_stack_4->GetZaxis()->SetLabelFont(42);
   hsig_data_all_stack_4->GetZaxis()->SetTitleOffset(1);
   hsig_data_all_stack_4->GetZaxis()->SetTitleFont(42);
   hs->Add(hsig_data_all_stack_4,"");
   hs->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("hsig_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("hsig_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("hsig_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("hsig_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->Modified();
   c->cd();
   c->SetSelected(c);
}