void PlotPointingRes()
{
//=========Macro generated from canvas: c/c
//=========  (Mon May 16 10:17:38 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,119.8484);
   hsig_Hino_135_Hyy_stack_1->SetBinContent(2,79.53299);
   hsig_Hino_135_Hyy_stack_1->SetBinContent(3,61.96319);
   hsig_Hino_135_Hyy_stack_1->SetBinContent(4,54.35619);
   hsig_Hino_135_Hyy_stack_1->SetBinContent(5,37.97338);
   hsig_Hino_135_Hyy_stack_1->SetBinContent(6,31.50521);
   hsig_Hino_135_Hyy_stack_1->SetBinContent(7,23.70501);
   hsig_Hino_135_Hyy_stack_1->SetBinContent(8,20.25397);
   hsig_Hino_135_Hyy_stack_1->SetBinContent(9,13.85109);
   hsig_Hino_135_Hyy_stack_1->SetBinContent(10,13.33894);
   hsig_Hino_135_Hyy_stack_1->SetBinContent(11,10.3361);
   hsig_Hino_135_Hyy_stack_1->SetBinContent(12,11.80116);
   hsig_Hino_135_Hyy_stack_1->SetBinContent(13,13.36342);
   hsig_Hino_135_Hyy_stack_1->SetBinContent(14,12.80615);
   hsig_Hino_135_Hyy_stack_1->SetBinContent(15,19.45242);
   hsig_Hino_135_Hyy_stack_1->SetBinContent(16,22.37934);
   hsig_Hino_135_Hyy_stack_1->SetBinContent(17,29.01382);
   hsig_Hino_135_Hyy_stack_1->SetBinContent(18,38.52811);
   hsig_Hino_135_Hyy_stack_1->SetBinContent(19,53.58908);
   hsig_Hino_135_Hyy_stack_1->SetBinContent(20,64.15059);
   hsig_Hino_135_Hyy_stack_1->SetBinContent(21,81.9718);
   hsig_Hino_135_Hyy_stack_1->SetBinContent(22,119.4573);
   hsig_Hino_135_Hyy_stack_1->SetBinError(1,4.257597);
   hsig_Hino_135_Hyy_stack_1->SetBinError(2,2.641256);
   hsig_Hino_135_Hyy_stack_1->SetBinError(3,3.009805);
   hsig_Hino_135_Hyy_stack_1->SetBinError(4,2.324481);
   hsig_Hino_135_Hyy_stack_1->SetBinError(5,1.631565);
   hsig_Hino_135_Hyy_stack_1->SetBinError(6,1.169908);
   hsig_Hino_135_Hyy_stack_1->SetBinError(7,1.127506);
   hsig_Hino_135_Hyy_stack_1->SetBinError(8,0.9025109);
   hsig_Hino_135_Hyy_stack_1->SetBinError(9,0.7176992);
   hsig_Hino_135_Hyy_stack_1->SetBinError(10,0.528372);
   hsig_Hino_135_Hyy_stack_1->SetBinError(11,0.4469939);
   hsig_Hino_135_Hyy_stack_1->SetBinError(12,0.4718177);
   hsig_Hino_135_Hyy_stack_1->SetBinError(13,0.5247985);
   hsig_Hino_135_Hyy_stack_1->SetBinError(14,0.5646459);
   hsig_Hino_135_Hyy_stack_1->SetBinError(15,0.983002);
   hsig_Hino_135_Hyy_stack_1->SetBinError(16,1.109804);
   hsig_Hino_135_Hyy_stack_1->SetBinError(17,1.589436);
   hsig_Hino_135_Hyy_stack_1->SetBinError(18,1.330116);
   hsig_Hino_135_Hyy_stack_1->SetBinError(19,2.211967);
   hsig_Hino_135_Hyy_stack_1->SetBinError(20,2.613463);
   hsig_Hino_135_Hyy_stack_1->SetBinError(21,3.062751);
   hsig_Hino_135_Hyy_stack_1->SetBinError(22,3.919309);
   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,105.7434);
   hsig_Hino_135_Zee_stack_2->SetBinContent(2,68.82827);
   hsig_Hino_135_Zee_stack_2->SetBinContent(3,61.19965);
   hsig_Hino_135_Zee_stack_2->SetBinContent(4,48.86045);
   hsig_Hino_135_Zee_stack_2->SetBinContent(5,40.29806);
   hsig_Hino_135_Zee_stack_2->SetBinContent(6,29.92201);
   hsig_Hino_135_Zee_stack_2->SetBinContent(7,27.04185);
   hsig_Hino_135_Zee_stack_2->SetBinContent(8,20.34592);
   hsig_Hino_135_Zee_stack_2->SetBinContent(9,14.54985);
   hsig_Hino_135_Zee_stack_2->SetBinContent(10,12.96849);
   hsig_Hino_135_Zee_stack_2->SetBinContent(11,10.46032);
   hsig_Hino_135_Zee_stack_2->SetBinContent(12,11.80038);
   hsig_Hino_135_Zee_stack_2->SetBinContent(13,11.78617);
   hsig_Hino_135_Zee_stack_2->SetBinContent(14,15.71257);
   hsig_Hino_135_Zee_stack_2->SetBinContent(15,18.58191);
   hsig_Hino_135_Zee_stack_2->SetBinContent(16,21.41976);
   hsig_Hino_135_Zee_stack_2->SetBinContent(17,29.96891);
   hsig_Hino_135_Zee_stack_2->SetBinContent(18,39.54059);
   hsig_Hino_135_Zee_stack_2->SetBinContent(19,53.11871);
   hsig_Hino_135_Zee_stack_2->SetBinContent(20,59.55907);
   hsig_Hino_135_Zee_stack_2->SetBinContent(21,71.80996);
   hsig_Hino_135_Zee_stack_2->SetBinContent(22,112.0623);
   hsig_Hino_135_Zee_stack_2->SetBinError(1,3.374766);
   hsig_Hino_135_Zee_stack_2->SetBinError(2,2.235111);
   hsig_Hino_135_Zee_stack_2->SetBinError(3,2.788368);
   hsig_Hino_135_Zee_stack_2->SetBinError(4,1.686733);
   hsig_Hino_135_Zee_stack_2->SetBinError(5,1.867264);
   hsig_Hino_135_Zee_stack_2->SetBinError(6,1.737683);
   hsig_Hino_135_Zee_stack_2->SetBinError(7,1.253575);
   hsig_Hino_135_Zee_stack_2->SetBinError(8,1.010188);
   hsig_Hino_135_Zee_stack_2->SetBinError(9,0.7205655);
   hsig_Hino_135_Zee_stack_2->SetBinError(10,0.6534386);
   hsig_Hino_135_Zee_stack_2->SetBinError(11,0.4770041);
   hsig_Hino_135_Zee_stack_2->SetBinError(12,0.7059391);
   hsig_Hino_135_Zee_stack_2->SetBinError(13,0.5734357);
   hsig_Hino_135_Zee_stack_2->SetBinError(14,0.7656096);
   hsig_Hino_135_Zee_stack_2->SetBinError(15,0.8835917);
   hsig_Hino_135_Zee_stack_2->SetBinError(16,1.032765);
   hsig_Hino_135_Zee_stack_2->SetBinError(17,1.824513);
   hsig_Hino_135_Zee_stack_2->SetBinError(18,1.47599);
   hsig_Hino_135_Zee_stack_2->SetBinError(19,2.270567);
   hsig_Hino_135_Zee_stack_2->SetBinError(20,2.615326);
   hsig_Hino_135_Zee_stack_2->SetBinError(21,2.233279);
   hsig_Hino_135_Zee_stack_2->SetBinError(22,3.564514);
   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,28.29477);
   hsig_Zee_stack_3->SetBinContent(7,22.45319);
   hsig_Zee_stack_3->SetBinContent(8,17.75502);
   hsig_Zee_stack_3->SetBinContent(9,15.24209);
   hsig_Zee_stack_3->SetBinContent(10,13.4921);
   hsig_Zee_stack_3->SetBinContent(11,12.76801);
   hsig_Zee_stack_3->SetBinContent(12,12.4687);
   hsig_Zee_stack_3->SetBinContent(13,13.62934);
   hsig_Zee_stack_3->SetBinContent(14,14.73515);
   hsig_Zee_stack_3->SetBinContent(15,17.69694);
   hsig_Zee_stack_3->SetBinContent(16,20.18929);
   hsig_Zee_stack_3->SetBinContent(17,24.68773);
   hsig_Zee_stack_3->SetBinError(6,3.58892);
   hsig_Zee_stack_3->SetBinError(7,0.7957209);
   hsig_Zee_stack_3->SetBinError(8,0.3503433);
   hsig_Zee_stack_3->SetBinError(9,0.1916545);
   hsig_Zee_stack_3->SetBinError(10,0.1111746);
   hsig_Zee_stack_3->SetBinError(11,0.09343322);
   hsig_Zee_stack_3->SetBinError(12,0.08791612);
   hsig_Zee_stack_3->SetBinError(13,0.1150452);
   hsig_Zee_stack_3->SetBinError(14,0.1751972);
   hsig_Zee_stack_3->SetBinError(15,0.3482794);
   hsig_Zee_stack_3->SetBinError(16,0.8206779);
   hsig_Zee_stack_3->SetBinError(17,2.430733);
   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,"");
   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);
   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);
}