void PlotPointingResVrho()
{
//=========Macro generated from canvas: c/c
//=========  (Tue May 31 10:43:13 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[21] = {-750, -400, -300, -140, -120, -100, -80, -60, -40, -20, 0, 20, 40, 60, 80, 100, 120, 140, 300, 400, 750}; 
   
   TH1F *hs_stack_1 = new TH1F("hs_stack_1","",20, 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[21] = {-750, -400, -300, -140, -120, -100, -80, -60, -40, -20, 0, 20, 40, 60, 80, 100, 120, 140, 300, 400, 750}; 
   
   TH1F *hsig_Hino_2ns_stack_1 = new TH1F("hsig_Hino_2ns_stack_1","",20, xAxis2);
   hsig_Hino_2ns_stack_1->SetBinContent(1,151.9577);
   hsig_Hino_2ns_stack_1->SetBinContent(2,86.16159);
   hsig_Hino_2ns_stack_1->SetBinContent(3,46.25172);
   hsig_Hino_2ns_stack_1->SetBinContent(4,28.30871);
   hsig_Hino_2ns_stack_1->SetBinContent(5,23.79232);
   hsig_Hino_2ns_stack_1->SetBinContent(6,17.77771);
   hsig_Hino_2ns_stack_1->SetBinContent(7,14.99352);
   hsig_Hino_2ns_stack_1->SetBinContent(8,11.90906);
   hsig_Hino_2ns_stack_1->SetBinContent(9,9.516913);
   hsig_Hino_2ns_stack_1->SetBinContent(10,9.88975);
   hsig_Hino_2ns_stack_1->SetBinContent(11,8.837668);
   hsig_Hino_2ns_stack_1->SetBinContent(12,9.48929);
   hsig_Hino_2ns_stack_1->SetBinContent(13,11.34502);
   hsig_Hino_2ns_stack_1->SetBinContent(14,13.49197);
   hsig_Hino_2ns_stack_1->SetBinContent(15,18.34609);
   hsig_Hino_2ns_stack_1->SetBinContent(16,21.73999);
   hsig_Hino_2ns_stack_1->SetBinContent(17,25.16727);
   hsig_Hino_2ns_stack_1->SetBinContent(18,45.04312);
   hsig_Hino_2ns_stack_1->SetBinContent(19,87.60326);
   hsig_Hino_2ns_stack_1->SetBinContent(20,147.8604);
   hsig_Hino_2ns_stack_1->SetBinError(1,1.878389);
   hsig_Hino_2ns_stack_1->SetBinError(2,1.109038);
   hsig_Hino_2ns_stack_1->SetBinError(3,0.3565159);
   hsig_Hino_2ns_stack_1->SetBinError(4,0.46465);
   hsig_Hino_2ns_stack_1->SetBinError(5,0.4679214);
   hsig_Hino_2ns_stack_1->SetBinError(6,0.3102414);
   hsig_Hino_2ns_stack_1->SetBinError(7,0.2420095);
   hsig_Hino_2ns_stack_1->SetBinError(8,0.1666201);
   hsig_Hino_2ns_stack_1->SetBinError(9,0.1310653);
   hsig_Hino_2ns_stack_1->SetBinError(10,0.1466994);
   hsig_Hino_2ns_stack_1->SetBinError(11,0.1026816);
   hsig_Hino_2ns_stack_1->SetBinError(12,0.1324708);
   hsig_Hino_2ns_stack_1->SetBinError(13,0.2376289);
   hsig_Hino_2ns_stack_1->SetBinError(14,0.2578978);
   hsig_Hino_2ns_stack_1->SetBinError(15,0.2617832);
   hsig_Hino_2ns_stack_1->SetBinError(16,0.3594349);
   hsig_Hino_2ns_stack_1->SetBinError(17,0.4771407);
   hsig_Hino_2ns_stack_1->SetBinError(18,0.3653861);
   hsig_Hino_2ns_stack_1->SetBinError(19,1.152753);
   hsig_Hino_2ns_stack_1->SetBinError(20,1.790982);
   hsig_Hino_2ns_stack_1->SetEntries(20);

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

   ci = TColor::GetColor("#00ff00");
   hsig_Hino_2ns_stack_1->SetMarkerColor(ci);
   hsig_Hino_2ns_stack_1->SetMarkerStyle(26);
   hsig_Hino_2ns_stack_1->SetMarkerSize(2);
   hsig_Hino_2ns_stack_1->GetXaxis()->SetTitle("Pointing Resolution [mm]");
   hsig_Hino_2ns_stack_1->GetXaxis()->SetLabelFont(42);
   hsig_Hino_2ns_stack_1->GetXaxis()->SetTitleOffset(1);
   hsig_Hino_2ns_stack_1->GetXaxis()->SetTitleFont(42);
   hsig_Hino_2ns_stack_1->GetYaxis()->SetTitle(" Pointing Resolution [mm]");
   hsig_Hino_2ns_stack_1->GetYaxis()->SetLabelFont(42);
   hsig_Hino_2ns_stack_1->GetYaxis()->SetTitleFont(42);
   hsig_Hino_2ns_stack_1->GetZaxis()->SetLabelFont(42);
   hsig_Hino_2ns_stack_1->GetZaxis()->SetTitleOffset(1);
   hsig_Hino_2ns_stack_1->GetZaxis()->SetTitleFont(42);
   hs->Add(hsig_Hino_2ns_stack_1,"");
   Double_t xAxis3[21] = {-750, -400, -300, -140, -120, -100, -80, -60, -40, -20, 0, 20, 40, 60, 80, 100, 120, 140, 300, 400, 750}; 
   
   TH1F *hsig_Hino_2ns_stack_2 = new TH1F("hsig_Hino_2ns_stack_2","",20, xAxis3);
   hsig_Hino_2ns_stack_2->SetBinContent(1,335.3139);
   hsig_Hino_2ns_stack_2->SetBinContent(2,263.2471);
   hsig_Hino_2ns_stack_2->SetBinContent(3,45.00698);
   hsig_Hino_2ns_stack_2->SetBinContent(4,32.54061);
   hsig_Hino_2ns_stack_2->SetBinContent(5,28.51133);
   hsig_Hino_2ns_stack_2->SetBinContent(6,20.30662);
   hsig_Hino_2ns_stack_2->SetBinContent(7,14.82661);
   hsig_Hino_2ns_stack_2->SetBinContent(8,11.95373);
   hsig_Hino_2ns_stack_2->SetBinContent(9,9.279762);
   hsig_Hino_2ns_stack_2->SetBinContent(10,9.143309);
   hsig_Hino_2ns_stack_2->SetBinContent(11,8.884949);
   hsig_Hino_2ns_stack_2->SetBinContent(12,10.52467);
   hsig_Hino_2ns_stack_2->SetBinContent(13,11.2731);
   hsig_Hino_2ns_stack_2->SetBinContent(14,13.49437);
   hsig_Hino_2ns_stack_2->SetBinContent(15,18.3974);
   hsig_Hino_2ns_stack_2->SetBinContent(16,23.41066);
   hsig_Hino_2ns_stack_2->SetBinContent(17,30.65145);
   hsig_Hino_2ns_stack_2->SetBinContent(18,45.33886);
   hsig_Hino_2ns_stack_2->SetBinContent(19,229.7699);
   hsig_Hino_2ns_stack_2->SetBinContent(20,323.0261);
   hsig_Hino_2ns_stack_2->SetBinError(1,256.1049);
   hsig_Hino_2ns_stack_2->SetBinError(2,222.9619);
   hsig_Hino_2ns_stack_2->SetBinError(3,4.533006);
   hsig_Hino_2ns_stack_2->SetBinError(4,3.623187);
   hsig_Hino_2ns_stack_2->SetBinError(5,2.613176);
   hsig_Hino_2ns_stack_2->SetBinError(6,1.290545);
   hsig_Hino_2ns_stack_2->SetBinError(7,0.6956964);
   hsig_Hino_2ns_stack_2->SetBinError(8,0.4241089);
   hsig_Hino_2ns_stack_2->SetBinError(9,0.2687505);
   hsig_Hino_2ns_stack_2->SetBinError(10,0.4767509);
   hsig_Hino_2ns_stack_2->SetBinError(11,0.2163525);
   hsig_Hino_2ns_stack_2->SetBinError(12,0.2489547);
   hsig_Hino_2ns_stack_2->SetBinError(13,0.576822);
   hsig_Hino_2ns_stack_2->SetBinError(14,0.7699426);
   hsig_Hino_2ns_stack_2->SetBinError(15,1.126713);
   hsig_Hino_2ns_stack_2->SetBinError(16,2.245219);
   hsig_Hino_2ns_stack_2->SetBinError(17,2.659292);
   hsig_Hino_2ns_stack_2->SetBinError(18,3.85844);
   hsig_Hino_2ns_stack_2->SetBinError(19,122.1499);
   hsig_Hino_2ns_stack_2->SetBinError(20,192.2908);
   hsig_Hino_2ns_stack_2->SetEntries(20);

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

   ci = TColor::GetColor("#ffcc00");
   hsig_Hino_2ns_stack_2->SetMarkerColor(ci);
   hsig_Hino_2ns_stack_2->SetMarkerStyle(26);
   hsig_Hino_2ns_stack_2->SetMarkerSize(2);
   hsig_Hino_2ns_stack_2->GetXaxis()->SetTitle("Pointing Resolution [mm]");
   hsig_Hino_2ns_stack_2->GetXaxis()->SetLabelFont(42);
   hsig_Hino_2ns_stack_2->GetXaxis()->SetTitleOffset(1);
   hsig_Hino_2ns_stack_2->GetXaxis()->SetTitleFont(42);
   hsig_Hino_2ns_stack_2->GetYaxis()->SetTitle(" Pointing Resolution [mm]");
   hsig_Hino_2ns_stack_2->GetYaxis()->SetLabelFont(42);
   hsig_Hino_2ns_stack_2->GetYaxis()->SetTitleFont(42);
   hsig_Hino_2ns_stack_2->GetZaxis()->SetLabelFont(42);
   hsig_Hino_2ns_stack_2->GetZaxis()->SetTitleOffset(1);
   hsig_Hino_2ns_stack_2->GetZaxis()->SetTitleFont(42);
   hs->Add(hsig_Hino_2ns_stack_2,"");
   Double_t xAxis4[21] = {-750, -400, -300, -140, -120, -100, -80, -60, -40, -20, 0, 20, 40, 60, 80, 100, 120, 140, 300, 400, 750}; 
   
   TH1F *hsig_Hino_2ns_stack_3 = new TH1F("hsig_Hino_2ns_stack_3","",20, xAxis4);
   hsig_Hino_2ns_stack_3->SetBinContent(1,383.7872);
   hsig_Hino_2ns_stack_3->SetBinContent(2,329.4429);
   hsig_Hino_2ns_stack_3->SetBinContent(3,97.19864);
   hsig_Hino_2ns_stack_3->SetBinContent(4,190.2463);
   hsig_Hino_2ns_stack_3->SetBinContent(5,30.80784);
   hsig_Hino_2ns_stack_3->SetBinContent(6,18.95185);
   hsig_Hino_2ns_stack_3->SetBinContent(7,13.13874);
   hsig_Hino_2ns_stack_3->SetBinContent(8,12.20965);
   hsig_Hino_2ns_stack_3->SetBinContent(9,11.9781);
   hsig_Hino_2ns_stack_3->SetBinContent(10,10.10868);
   hsig_Hino_2ns_stack_3->SetBinContent(11,8.941812);
   hsig_Hino_2ns_stack_3->SetBinContent(12,9.009976);
   hsig_Hino_2ns_stack_3->SetBinContent(13,11.22417);
   hsig_Hino_2ns_stack_3->SetBinContent(14,14.354);
   hsig_Hino_2ns_stack_3->SetBinContent(15,15.26798);
   hsig_Hino_2ns_stack_3->SetBinContent(16,60.47858);
   hsig_Hino_2ns_stack_3->SetBinContent(17,46.60015);
   hsig_Hino_2ns_stack_3->SetBinContent(18,147.5839);
   hsig_Hino_2ns_stack_3->SetBinContent(19,230.2374);
   hsig_Hino_2ns_stack_3->SetBinContent(20,1686.242);
   hsig_Hino_2ns_stack_3->SetBinError(1,697.4059);
   hsig_Hino_2ns_stack_3->SetBinError(2,885.5746);
   hsig_Hino_2ns_stack_3->SetBinError(3,64.91308);
   hsig_Hino_2ns_stack_3->SetBinError(4,158.5467);
   hsig_Hino_2ns_stack_3->SetBinError(5,11.05758);
   hsig_Hino_2ns_stack_3->SetBinError(6,3.337204);
   hsig_Hino_2ns_stack_3->SetBinError(7,1.336888);
   hsig_Hino_2ns_stack_3->SetBinError(8,1.151457);
   hsig_Hino_2ns_stack_3->SetBinError(9,0.7825452);
   hsig_Hino_2ns_stack_3->SetBinError(10,0.4323735);
   hsig_Hino_2ns_stack_3->SetBinError(11,0.4902937);
   hsig_Hino_2ns_stack_3->SetBinError(12,0.5537802);
   hsig_Hino_2ns_stack_3->SetBinError(13,0.8500611);
   hsig_Hino_2ns_stack_3->SetBinError(14,1.808121);
   hsig_Hino_2ns_stack_3->SetBinError(15,3.229491);
   hsig_Hino_2ns_stack_3->SetBinError(16,43.6526);
   hsig_Hino_2ns_stack_3->SetBinError(17,52.03905);
   hsig_Hino_2ns_stack_3->SetBinError(18,114.5854);
   hsig_Hino_2ns_stack_3->SetBinError(19,514.9638);
   hsig_Hino_2ns_stack_3->SetBinError(20,843.1483);
   hsig_Hino_2ns_stack_3->SetEntries(20);

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

   ci = TColor::GetColor("#ff00ff");
   hsig_Hino_2ns_stack_3->SetMarkerColor(ci);
   hsig_Hino_2ns_stack_3->SetMarkerStyle(26);
   hsig_Hino_2ns_stack_3->SetMarkerSize(2);
   hsig_Hino_2ns_stack_3->GetXaxis()->SetTitle("Pointing Resolution [mm]");
   hsig_Hino_2ns_stack_3->GetXaxis()->SetLabelFont(42);
   hsig_Hino_2ns_stack_3->GetXaxis()->SetTitleOffset(1);
   hsig_Hino_2ns_stack_3->GetXaxis()->SetTitleFont(42);
   hsig_Hino_2ns_stack_3->GetYaxis()->SetTitle(" Pointing Resolution [mm]");
   hsig_Hino_2ns_stack_3->GetYaxis()->SetLabelFont(42);
   hsig_Hino_2ns_stack_3->GetYaxis()->SetTitleFont(42);
   hsig_Hino_2ns_stack_3->GetZaxis()->SetLabelFont(42);
   hsig_Hino_2ns_stack_3->GetZaxis()->SetTitleOffset(1);
   hsig_Hino_2ns_stack_3->GetZaxis()->SetTitleFont(42);
   hs->Add(hsig_Hino_2ns_stack_3,"");
   Double_t xAxis5[21] = {-750, -400, -300, -140, -120, -100, -80, -60, -40, -20, 0, 20, 40, 60, 80, 100, 120, 140, 300, 400, 750}; 
   
   TH1F *hsig_Zee_truZ_stack_4 = new TH1F("hsig_Zee_truZ_stack_4","",20, xAxis5);
   hsig_Zee_truZ_stack_4->SetBinContent(4,20.37015);
   hsig_Zee_truZ_stack_4->SetBinContent(5,22.11153);
   hsig_Zee_truZ_stack_4->SetBinContent(6,16.67682);
   hsig_Zee_truZ_stack_4->SetBinContent(7,15.04218);
   hsig_Zee_truZ_stack_4->SetBinContent(8,12.26145);
   hsig_Zee_truZ_stack_4->SetBinContent(9,10.66822);
   hsig_Zee_truZ_stack_4->SetBinContent(10,11.24329);
   hsig_Zee_truZ_stack_4->SetBinContent(11,10.52154);
   hsig_Zee_truZ_stack_4->SetBinContent(12,11.53468);
   hsig_Zee_truZ_stack_4->SetBinContent(13,12.11427);
   hsig_Zee_truZ_stack_4->SetBinContent(14,14.08047);
   hsig_Zee_truZ_stack_4->SetBinContent(15,16.6785);
   hsig_Zee_truZ_stack_4->SetBinContent(16,17.57574);
   hsig_Zee_truZ_stack_4->SetBinContent(17,21.90544);
   hsig_Zee_truZ_stack_4->SetBinContent(18,18.66935);
   hsig_Zee_truZ_stack_4->SetBinError(4,1.424691);
   hsig_Zee_truZ_stack_4->SetBinError(5,1.402546);
   hsig_Zee_truZ_stack_4->SetBinError(6,0.4830379);
   hsig_Zee_truZ_stack_4->SetBinError(7,0.3219565);
   hsig_Zee_truZ_stack_4->SetBinError(8,0.1891684);
   hsig_Zee_truZ_stack_4->SetBinError(9,0.1621675);
   hsig_Zee_truZ_stack_4->SetBinError(10,0.1043363);
   hsig_Zee_truZ_stack_4->SetBinError(11,0.1390061);
   hsig_Zee_truZ_stack_4->SetBinError(12,0.125853);
   hsig_Zee_truZ_stack_4->SetBinError(13,0.182362);
   hsig_Zee_truZ_stack_4->SetBinError(14,0.2710565);
   hsig_Zee_truZ_stack_4->SetBinError(15,0.4973943);
   hsig_Zee_truZ_stack_4->SetBinError(16,1.005759);
   hsig_Zee_truZ_stack_4->SetBinError(17,2.273173);
   hsig_Zee_truZ_stack_4->SetBinError(18,4.333411);
   hsig_Zee_truZ_stack_4->SetEntries(15);

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

   ci = TColor::GetColor("#0000ff");
   hsig_Zee_truZ_stack_4->SetMarkerColor(ci);
   hsig_Zee_truZ_stack_4->SetMarkerStyle(20);
   hsig_Zee_truZ_stack_4->SetMarkerSize(2);
   hsig_Zee_truZ_stack_4->GetXaxis()->SetTitle("Pointing Resolution [mm]");
   hsig_Zee_truZ_stack_4->GetXaxis()->SetLabelFont(42);
   hsig_Zee_truZ_stack_4->GetXaxis()->SetTitleOffset(1);
   hsig_Zee_truZ_stack_4->GetXaxis()->SetTitleFont(42);
   hsig_Zee_truZ_stack_4->GetYaxis()->SetTitle(" Pointing Resolution [mm]");
   hsig_Zee_truZ_stack_4->GetYaxis()->SetLabelFont(42);
   hsig_Zee_truZ_stack_4->GetYaxis()->SetTitleFont(42);
   hsig_Zee_truZ_stack_4->GetZaxis()->SetLabelFont(42);
   hsig_Zee_truZ_stack_4->GetZaxis()->SetTitleOffset(1);
   hsig_Zee_truZ_stack_4->GetZaxis()->SetTitleFont(42);
   hs->Add(hsig_Zee_truZ_stack_4,"");
   hs->Draw("nostack p");
   
   TLegend *leg = new TLegend(0.45,0.7,0.62,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_2ns_stack_1","Hino_2ns:#tau_{Hino} < 100.00 ns","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("hsig_Hino_2ns_stack_2","Hino_2ns:#tau_{Hino} < 0.20 ns","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("hsig_Hino_2ns_stack_3","Hino_2ns:#tau_{Hino} < 0.05 ns","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_Zee_truZ_stack_4","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->Modified();
   c->cd();
   c->SetSelected(c);
}