void PlotPointingRes()
{
//=========Macro generated from canvas: c/c
//=========  (Mon May 23 10:16:29 2022) by ROOT version 6.20/06
   TCanvas *c = new TCanvas("c", "c",0,0,1600,1200);
   gStyle->SetOptStat(0);
   c->Range(-937.5,-18.75,937.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[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_stack_1 = new TH1F("hs_stack_1","",22, 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[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 *hsig_Hino_135_Hyy_stack_1 = new TH1F("hsig_Hino_135_Hyy_stack_1","",22, xAxis2);
   hsig_Hino_135_Hyy_stack_1->SetBinContent(1,98.64541);
   hsig_Hino_135_Hyy_stack_1->SetBinContent(2,69.07349);
   hsig_Hino_135_Hyy_stack_1->SetBinContent(3,49.48498);
   hsig_Hino_135_Hyy_stack_1->SetBinContent(4,48.10503);
   hsig_Hino_135_Hyy_stack_1->SetBinContent(5,33.04204);
   hsig_Hino_135_Hyy_stack_1->SetBinContent(6,27.2696);
   hsig_Hino_135_Hyy_stack_1->SetBinContent(7,20.99718);
   hsig_Hino_135_Hyy_stack_1->SetBinContent(8,17.476);
   hsig_Hino_135_Hyy_stack_1->SetBinContent(9,11.54739);
   hsig_Hino_135_Hyy_stack_1->SetBinContent(10,10.4617);
   hsig_Hino_135_Hyy_stack_1->SetBinContent(11,8.814338);
   hsig_Hino_135_Hyy_stack_1->SetBinContent(12,8.895816);
   hsig_Hino_135_Hyy_stack_1->SetBinContent(13,11.26875);
   hsig_Hino_135_Hyy_stack_1->SetBinContent(14,11.57821);
   hsig_Hino_135_Hyy_stack_1->SetBinContent(15,16.12786);
   hsig_Hino_135_Hyy_stack_1->SetBinContent(16,20.91427);
   hsig_Hino_135_Hyy_stack_1->SetBinContent(17,25.021);
   hsig_Hino_135_Hyy_stack_1->SetBinContent(18,35.46091);
   hsig_Hino_135_Hyy_stack_1->SetBinContent(19,46.98877);
   hsig_Hino_135_Hyy_stack_1->SetBinContent(20,51.57434);
   hsig_Hino_135_Hyy_stack_1->SetBinContent(21,68.66442);
   hsig_Hino_135_Hyy_stack_1->SetBinContent(22,92.37844);
   hsig_Hino_135_Hyy_stack_1->SetBinError(1,4.484064);
   hsig_Hino_135_Hyy_stack_1->SetBinError(2,2.468139);
   hsig_Hino_135_Hyy_stack_1->SetBinError(3,2.899005);
   hsig_Hino_135_Hyy_stack_1->SetBinError(4,2.886529);
   hsig_Hino_135_Hyy_stack_1->SetBinError(5,1.704821);
   hsig_Hino_135_Hyy_stack_1->SetBinError(6,1.524902);
   hsig_Hino_135_Hyy_stack_1->SetBinError(7,1.232663);
   hsig_Hino_135_Hyy_stack_1->SetBinError(8,0.9397301);
   hsig_Hino_135_Hyy_stack_1->SetBinError(9,0.6321346);
   hsig_Hino_135_Hyy_stack_1->SetBinError(10,0.5512);
   hsig_Hino_135_Hyy_stack_1->SetBinError(11,0.4024501);
   hsig_Hino_135_Hyy_stack_1->SetBinError(12,0.4317606);
   hsig_Hino_135_Hyy_stack_1->SetBinError(13,0.4962507);
   hsig_Hino_135_Hyy_stack_1->SetBinError(14,0.5834724);
   hsig_Hino_135_Hyy_stack_1->SetBinError(15,0.8835087);
   hsig_Hino_135_Hyy_stack_1->SetBinError(16,1.029426);
   hsig_Hino_135_Hyy_stack_1->SetBinError(17,1.203293);
   hsig_Hino_135_Hyy_stack_1->SetBinError(18,2.009265);
   hsig_Hino_135_Hyy_stack_1->SetBinError(19,2.596437);
   hsig_Hino_135_Hyy_stack_1->SetBinError(20,3.035047);
   hsig_Hino_135_Hyy_stack_1->SetBinError(21,3.433821);
   hsig_Hino_135_Hyy_stack_1->SetBinError(22,3.829138);
   hsig_Hino_135_Hyy_stack_1->SetEntries(22);

   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[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 *hsig_Hino_135_Zee_stack_2 = new TH1F("hsig_Hino_135_Zee_stack_2","",22, xAxis3);
   hsig_Hino_135_Zee_stack_2->SetBinContent(1,84.02205);
   hsig_Hino_135_Zee_stack_2->SetBinContent(2,62.31391);
   hsig_Hino_135_Zee_stack_2->SetBinContent(3,53.96534);
   hsig_Hino_135_Zee_stack_2->SetBinContent(4,43.27344);
   hsig_Hino_135_Zee_stack_2->SetBinContent(5,35.04352);
   hsig_Hino_135_Zee_stack_2->SetBinContent(6,25.59093);
   hsig_Hino_135_Zee_stack_2->SetBinContent(7,22.83068);
   hsig_Hino_135_Zee_stack_2->SetBinContent(8,16.33397);
   hsig_Hino_135_Zee_stack_2->SetBinContent(9,13.57303);
   hsig_Hino_135_Zee_stack_2->SetBinContent(10,10.4972);
   hsig_Hino_135_Zee_stack_2->SetBinContent(11,9.782319);
   hsig_Hino_135_Zee_stack_2->SetBinContent(12,10.7781);
   hsig_Hino_135_Zee_stack_2->SetBinContent(13,10.9332);
   hsig_Hino_135_Zee_stack_2->SetBinContent(14,13.23772);
   hsig_Hino_135_Zee_stack_2->SetBinContent(15,18.27172);
   hsig_Hino_135_Zee_stack_2->SetBinContent(16,20.31099);
   hsig_Hino_135_Zee_stack_2->SetBinContent(17,26.93551);
   hsig_Hino_135_Zee_stack_2->SetBinContent(18,31.25374);
   hsig_Hino_135_Zee_stack_2->SetBinContent(19,45.26451);
   hsig_Hino_135_Zee_stack_2->SetBinContent(20,54.52573);
   hsig_Hino_135_Zee_stack_2->SetBinContent(21,66.01048);
   hsig_Hino_135_Zee_stack_2->SetBinContent(22,90.20917);
   hsig_Hino_135_Zee_stack_2->SetBinError(1,3.142331);
   hsig_Hino_135_Zee_stack_2->SetBinError(2,2.672663);
   hsig_Hino_135_Zee_stack_2->SetBinError(3,2.704861);
   hsig_Hino_135_Zee_stack_2->SetBinError(4,2.038169);
   hsig_Hino_135_Zee_stack_2->SetBinError(5,1.989608);
   hsig_Hino_135_Zee_stack_2->SetBinError(6,1.2349);
   hsig_Hino_135_Zee_stack_2->SetBinError(7,1.338115);
   hsig_Hino_135_Zee_stack_2->SetBinError(8,0.9013076);
   hsig_Hino_135_Zee_stack_2->SetBinError(9,0.7286341);
   hsig_Hino_135_Zee_stack_2->SetBinError(10,0.5239369);
   hsig_Hino_135_Zee_stack_2->SetBinError(11,0.4865178);
   hsig_Hino_135_Zee_stack_2->SetBinError(12,0.6030272);
   hsig_Hino_135_Zee_stack_2->SetBinError(13,0.5519032);
   hsig_Hino_135_Zee_stack_2->SetBinError(14,0.70834);
   hsig_Hino_135_Zee_stack_2->SetBinError(15,1.173099);
   hsig_Hino_135_Zee_stack_2->SetBinError(16,1.109912);
   hsig_Hino_135_Zee_stack_2->SetBinError(17,1.586313);
   hsig_Hino_135_Zee_stack_2->SetBinError(18,1.402527);
   hsig_Hino_135_Zee_stack_2->SetBinError(19,2.422382);
   hsig_Hino_135_Zee_stack_2->SetBinError(20,3.55739);
   hsig_Hino_135_Zee_stack_2->SetBinError(21,2.933567);
   hsig_Hino_135_Zee_stack_2->SetBinError(22,4.020144);
   hsig_Hino_135_Zee_stack_2->SetEntries(22);

   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[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 *hsig_Zee_stack_3 = new TH1F("hsig_Zee_stack_3","",22, xAxis4);
   hsig_Zee_stack_3->SetBinContent(6,20.06786);
   hsig_Zee_stack_3->SetBinContent(7,18.50749);
   hsig_Zee_stack_3->SetBinContent(8,14.5961);
   hsig_Zee_stack_3->SetBinContent(9,11.55928);
   hsig_Zee_stack_3->SetBinContent(10,10.02466);
   hsig_Zee_stack_3->SetBinContent(11,9.342149);
   hsig_Zee_stack_3->SetBinContent(12,9.116945);
   hsig_Zee_stack_3->SetBinContent(13,10.07881);
   hsig_Zee_stack_3->SetBinContent(14,11.3309);
   hsig_Zee_stack_3->SetBinContent(15,13.60473);
   hsig_Zee_stack_3->SetBinContent(16,15.97229);
   hsig_Zee_stack_3->SetBinContent(17,20.85042);
   hsig_Zee_stack_3->SetBinError(6,1.287907);
   hsig_Zee_stack_3->SetBinError(7,1.077778);
   hsig_Zee_stack_3->SetBinError(8,0.5190371);
   hsig_Zee_stack_3->SetBinError(9,0.2127911);
   hsig_Zee_stack_3->SetBinError(10,0.1482353);
   hsig_Zee_stack_3->SetBinError(11,0.1256035);
   hsig_Zee_stack_3->SetBinError(12,0.1180328);
   hsig_Zee_stack_3->SetBinError(13,0.1504221);
   hsig_Zee_stack_3->SetBinError(14,0.2420328);
   hsig_Zee_stack_3->SetBinError(15,0.4479278);
   hsig_Zee_stack_3->SetBinError(16,1.14396);
   hsig_Zee_stack_3->SetBinError(17,3.61993);
   hsig_Zee_stack_3->SetEntries(12);

   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[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 *hsig_data_all_stack_4 = new TH1F("hsig_data_all_stack_4","",22, xAxis5);
   hsig_data_all_stack_4->SetBinContent(6,20.57784);
   hsig_data_all_stack_4->SetBinContent(7,16.94414);
   hsig_data_all_stack_4->SetBinContent(8,14.25091);
   hsig_data_all_stack_4->SetBinContent(9,11.68749);
   hsig_data_all_stack_4->SetBinContent(10,10.29713);
   hsig_data_all_stack_4->SetBinContent(11,9.848927);
   hsig_data_all_stack_4->SetBinContent(12,9.852605);
   hsig_data_all_stack_4->SetBinContent(13,10.18299);
   hsig_data_all_stack_4->SetBinContent(14,11.40652);
   hsig_data_all_stack_4->SetBinContent(15,14.40001);
   hsig_data_all_stack_4->SetBinContent(16,16.58342);
   hsig_data_all_stack_4->SetBinContent(17,21.5459);
   hsig_data_all_stack_4->SetBinError(6,2.052405);
   hsig_data_all_stack_4->SetBinError(7,0.6088256);
   hsig_data_all_stack_4->SetBinError(8,0.1764385);
   hsig_data_all_stack_4->SetBinError(9,0.07711715);
   hsig_data_all_stack_4->SetBinError(10,0.04824938);
   hsig_data_all_stack_4->SetBinError(11,0.03410348);
   hsig_data_all_stack_4->SetBinError(12,0.03504824);
   hsig_data_all_stack_4->SetBinError(13,0.05080615);
   hsig_data_all_stack_4->SetBinError(14,0.08190694);
   hsig_data_all_stack_4->SetBinError(15,0.2303701);
   hsig_data_all_stack_4->SetBinError(16,0.6134965);
   hsig_data_all_stack_4->SetBinError(17,2.902713);
   hsig_data_all_stack_4->SetEntries(12);
   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.55,0.7,0.72,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);
}