void PlotPointingBiasVrho()
{
//=========Macro generated from canvas: c_m/c
//=========  (Thu Jun 16 09:39:03 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(-250,-250,2250,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 Resolution [mm];");
   hs_m->SetMinimum(-200);
   hs_m->SetMaximum(200);
   Double_t xAxis1[10] = {0, 20, 40, 60, 80, 100, 120, 300, 500, 2000}; 
   
   TH1F *hs_m_stack_1 = new TH1F("hs_m_stack_1","",9, xAxis1);
   hs_m_stack_1->SetMinimum(-200);
   hs_m_stack_1->SetMaximum(200);
   hs_m_stack_1->SetDirectory(0);
   hs_m_stack_1->SetStats(0);

   Int_t ci;      // for color index setting
   TColor *color; // for color definition with alpha
   ci = TColor::GetColor("#000099");
   hs_m_stack_1->SetLineColor(ci);
   hs_m_stack_1->GetXaxis()->SetTitle("|z_{origin}| [mm]");
   hs_m_stack_1->GetXaxis()->SetLabelFont(42);
   hs_m_stack_1->GetXaxis()->SetTitleOffset(1);
   hs_m_stack_1->GetXaxis()->SetTitleFont(42);
   hs_m_stack_1->GetYaxis()->SetTitle(" Pointing Resolution [mm]");
   hs_m_stack_1->GetYaxis()->SetLabelFont(42);
   hs_m_stack_1->GetYaxis()->SetTitleFont(42);
   hs_m_stack_1->GetZaxis()->SetLabelFont(42);
   hs_m_stack_1->GetZaxis()->SetTitleOffset(1);
   hs_m_stack_1->GetZaxis()->SetTitleFont(42);
   hs_m->SetHistogram(hs_m_stack_1);
   
   Double_t xAxis2[10] = {0, 20, 40, 60, 80, 100, 120, 300, 500, 2000}; 
   
   TH1F *hmean_Hino_135_2ns_HyyHSM_stack_1 = new TH1F("hmean_Hino_135_2ns_HyyHSM_stack_1","",9, xAxis2);
   hmean_Hino_135_2ns_HyyHSM_stack_1->SetBinContent(1,8.317298);
   hmean_Hino_135_2ns_HyyHSM_stack_1->SetBinContent(2,-37.77903);
   hmean_Hino_135_2ns_HyyHSM_stack_1->SetBinContent(3,-10.09165);
   hmean_Hino_135_2ns_HyyHSM_stack_1->SetBinContent(4,20.36626);
   hmean_Hino_135_2ns_HyyHSM_stack_1->SetBinContent(5,-47.60487);
   hmean_Hino_135_2ns_HyyHSM_stack_1->SetBinContent(6,9.88545);
   hmean_Hino_135_2ns_HyyHSM_stack_1->SetBinContent(7,23.28899);
   hmean_Hino_135_2ns_HyyHSM_stack_1->SetBinContent(8,59.35661);
   hmean_Hino_135_2ns_HyyHSM_stack_1->SetBinContent(9,170.4682);
   hmean_Hino_135_2ns_HyyHSM_stack_1->SetBinError(1,12.33144);
   hmean_Hino_135_2ns_HyyHSM_stack_1->SetBinError(2,133.4347);
   hmean_Hino_135_2ns_HyyHSM_stack_1->SetBinError(3,198.4193);
   hmean_Hino_135_2ns_HyyHSM_stack_1->SetBinError(4,98.21629);
   hmean_Hino_135_2ns_HyyHSM_stack_1->SetBinError(5,225.8613);
   hmean_Hino_135_2ns_HyyHSM_stack_1->SetBinError(6,35.10566);
   hmean_Hino_135_2ns_HyyHSM_stack_1->SetBinError(7,5.400688);
   hmean_Hino_135_2ns_HyyHSM_stack_1->SetBinError(8,13.63337);
   hmean_Hino_135_2ns_HyyHSM_stack_1->SetBinError(9,213.4588);
   hmean_Hino_135_2ns_HyyHSM_stack_1->SetEntries(9);

   ci = TColor::GetColor("#ff66ff");
   hmean_Hino_135_2ns_HyyHSM_stack_1->SetLineColor(ci);

   ci = TColor::GetColor("#ff66ff");
   hmean_Hino_135_2ns_HyyHSM_stack_1->SetMarkerColor(ci);
   hmean_Hino_135_2ns_HyyHSM_stack_1->SetMarkerStyle(26);
   hmean_Hino_135_2ns_HyyHSM_stack_1->SetMarkerSize(2);
   hmean_Hino_135_2ns_HyyHSM_stack_1->GetXaxis()->SetLabelFont(42);
   hmean_Hino_135_2ns_HyyHSM_stack_1->GetXaxis()->SetTitleOffset(1);
   hmean_Hino_135_2ns_HyyHSM_stack_1->GetXaxis()->SetTitleFont(42);
   hmean_Hino_135_2ns_HyyHSM_stack_1->GetYaxis()->SetLabelFont(42);
   hmean_Hino_135_2ns_HyyHSM_stack_1->GetYaxis()->SetTitleFont(42);
   hmean_Hino_135_2ns_HyyHSM_stack_1->GetZaxis()->SetLabelFont(42);
   hmean_Hino_135_2ns_HyyHSM_stack_1->GetZaxis()->SetTitleOffset(1);
   hmean_Hino_135_2ns_HyyHSM_stack_1->GetZaxis()->SetTitleFont(42);
   hs_m->Add(hmean_Hino_135_2ns_HyyHSM_stack_1,"");
   Double_t xAxis3[10] = {0, 20, 40, 60, 80, 100, 120, 300, 500, 2000}; 
   
   TH1F *hmean_Hino_135_2ns_ZeeZSM_stack_2 = new TH1F("hmean_Hino_135_2ns_ZeeZSM_stack_2","",9, xAxis3);
   hmean_Hino_135_2ns_ZeeZSM_stack_2->SetBinContent(1,-1.811006);
   hmean_Hino_135_2ns_ZeeZSM_stack_2->SetBinContent(2,5.869547);
   hmean_Hino_135_2ns_ZeeZSM_stack_2->SetBinContent(3,-59.9592);
   hmean_Hino_135_2ns_ZeeZSM_stack_2->SetBinContent(4,-41.93364);
   hmean_Hino_135_2ns_ZeeZSM_stack_2->SetBinContent(5,-5.187752);
   hmean_Hino_135_2ns_ZeeZSM_stack_2->SetBinContent(6,15.54259);
   hmean_Hino_135_2ns_ZeeZSM_stack_2->SetBinContent(7,-4.068924);
   hmean_Hino_135_2ns_ZeeZSM_stack_2->SetBinContent(8,28.65927);
   hmean_Hino_135_2ns_ZeeZSM_stack_2->SetBinContent(9,6.22108);
   hmean_Hino_135_2ns_ZeeZSM_stack_2->SetBinError(1,115.707);
   hmean_Hino_135_2ns_ZeeZSM_stack_2->SetBinError(2,228.004);
   hmean_Hino_135_2ns_ZeeZSM_stack_2->SetBinError(3,227.8956);
   hmean_Hino_135_2ns_ZeeZSM_stack_2->SetBinError(4,222.2318);
   hmean_Hino_135_2ns_ZeeZSM_stack_2->SetBinError(5,172.0877);
   hmean_Hino_135_2ns_ZeeZSM_stack_2->SetBinError(6,43.34871);
   hmean_Hino_135_2ns_ZeeZSM_stack_2->SetBinError(7,6.300182);
   hmean_Hino_135_2ns_ZeeZSM_stack_2->SetBinError(8,21.64894);
   hmean_Hino_135_2ns_ZeeZSM_stack_2->SetBinError(9,88.07653);
   hmean_Hino_135_2ns_ZeeZSM_stack_2->SetEntries(9);

   ci = TColor::GetColor("#009900");
   hmean_Hino_135_2ns_ZeeZSM_stack_2->SetLineColor(ci);

   ci = TColor::GetColor("#009900");
   hmean_Hino_135_2ns_ZeeZSM_stack_2->SetMarkerColor(ci);
   hmean_Hino_135_2ns_ZeeZSM_stack_2->SetMarkerStyle(26);
   hmean_Hino_135_2ns_ZeeZSM_stack_2->SetMarkerSize(2);
   hmean_Hino_135_2ns_ZeeZSM_stack_2->GetXaxis()->SetLabelFont(42);
   hmean_Hino_135_2ns_ZeeZSM_stack_2->GetXaxis()->SetTitleOffset(1);
   hmean_Hino_135_2ns_ZeeZSM_stack_2->GetXaxis()->SetTitleFont(42);
   hmean_Hino_135_2ns_ZeeZSM_stack_2->GetYaxis()->SetLabelFont(42);
   hmean_Hino_135_2ns_ZeeZSM_stack_2->GetYaxis()->SetTitleFont(42);
   hmean_Hino_135_2ns_ZeeZSM_stack_2->GetZaxis()->SetLabelFont(42);
   hmean_Hino_135_2ns_ZeeZSM_stack_2->GetZaxis()->SetTitleOffset(1);
   hmean_Hino_135_2ns_ZeeZSM_stack_2->GetZaxis()->SetTitleFont(42);
   hs_m->Add(hmean_Hino_135_2ns_ZeeZSM_stack_2,"");
   Double_t xAxis4[10] = {0, 20, 40, 60, 80, 100, 120, 300, 500, 2000}; 
   
   TH1F *hmean_Hino_135_10ns_HyyHSM_stack_3 = new TH1F("hmean_Hino_135_10ns_HyyHSM_stack_3","",9, xAxis4);
   hmean_Hino_135_10ns_HyyHSM_stack_3->SetBinContent(1,130.1398);
   hmean_Hino_135_10ns_HyyHSM_stack_3->SetBinContent(2,7.669408);
   hmean_Hino_135_10ns_HyyHSM_stack_3->SetBinContent(3,31.41259);
   hmean_Hino_135_10ns_HyyHSM_stack_3->SetBinContent(4,27.97955);
   hmean_Hino_135_10ns_HyyHSM_stack_3->SetBinContent(5,-4.663884);
   hmean_Hino_135_10ns_HyyHSM_stack_3->SetBinContent(6,95.10497);
   hmean_Hino_135_10ns_HyyHSM_stack_3->SetBinContent(7,49.27845);
   hmean_Hino_135_10ns_HyyHSM_stack_3->SetBinContent(8,99.51669);
   hmean_Hino_135_10ns_HyyHSM_stack_3->SetBinContent(9,-40.66927);
   hmean_Hino_135_10ns_HyyHSM_stack_3->SetBinError(1,1845.438);
   hmean_Hino_135_10ns_HyyHSM_stack_3->SetBinError(2,994.0897);
   hmean_Hino_135_10ns_HyyHSM_stack_3->SetBinError(3,247.3444);
   hmean_Hino_135_10ns_HyyHSM_stack_3->SetBinError(4,965.5057);
   hmean_Hino_135_10ns_HyyHSM_stack_3->SetBinError(5,2248.891);
   hmean_Hino_135_10ns_HyyHSM_stack_3->SetBinError(6,313.3614);
   hmean_Hino_135_10ns_HyyHSM_stack_3->SetBinError(7,33.13779);
   hmean_Hino_135_10ns_HyyHSM_stack_3->SetBinError(8,242.8156);
   hmean_Hino_135_10ns_HyyHSM_stack_3->SetBinError(9,128.7137);
   hmean_Hino_135_10ns_HyyHSM_stack_3->SetEntries(9);

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

   ci = TColor::GetColor("#ffcc00");
   hmean_Hino_135_10ns_HyyHSM_stack_3->SetMarkerColor(ci);
   hmean_Hino_135_10ns_HyyHSM_stack_3->SetMarkerStyle(26);
   hmean_Hino_135_10ns_HyyHSM_stack_3->SetMarkerSize(2);
   hmean_Hino_135_10ns_HyyHSM_stack_3->GetXaxis()->SetLabelFont(42);
   hmean_Hino_135_10ns_HyyHSM_stack_3->GetXaxis()->SetTitleOffset(1);
   hmean_Hino_135_10ns_HyyHSM_stack_3->GetXaxis()->SetTitleFont(42);
   hmean_Hino_135_10ns_HyyHSM_stack_3->GetYaxis()->SetLabelFont(42);
   hmean_Hino_135_10ns_HyyHSM_stack_3->GetYaxis()->SetTitleFont(42);
   hmean_Hino_135_10ns_HyyHSM_stack_3->GetZaxis()->SetLabelFont(42);
   hmean_Hino_135_10ns_HyyHSM_stack_3->GetZaxis()->SetTitleOffset(1);
   hmean_Hino_135_10ns_HyyHSM_stack_3->GetZaxis()->SetTitleFont(42);
   hs_m->Add(hmean_Hino_135_10ns_HyyHSM_stack_3,"");
   Double_t xAxis5[10] = {0, 20, 40, 60, 80, 100, 120, 300, 500, 2000}; 
   
   TH1F *hmean_Hino_135_10ns_ZeeZSM_stack_4 = new TH1F("hmean_Hino_135_10ns_ZeeZSM_stack_4","",9, xAxis5);
   hmean_Hino_135_10ns_ZeeZSM_stack_4->SetBinContent(1,26.94438);
   hmean_Hino_135_10ns_ZeeZSM_stack_4->SetBinContent(2,39.49816);
   hmean_Hino_135_10ns_ZeeZSM_stack_4->SetBinContent(3,66.64555);
   hmean_Hino_135_10ns_ZeeZSM_stack_4->SetBinContent(4,-135.6647);
   hmean_Hino_135_10ns_ZeeZSM_stack_4->SetBinContent(5,-73.14718);
   hmean_Hino_135_10ns_ZeeZSM_stack_4->SetBinContent(6,18.01078);
   hmean_Hino_135_10ns_ZeeZSM_stack_4->SetBinContent(7,4.185076);
   hmean_Hino_135_10ns_ZeeZSM_stack_4->SetBinContent(8,65.1384);
   hmean_Hino_135_10ns_ZeeZSM_stack_4->SetBinContent(9,8.507764);
   hmean_Hino_135_10ns_ZeeZSM_stack_4->SetBinError(1,2785.248);
   hmean_Hino_135_10ns_ZeeZSM_stack_4->SetBinError(2,2823.797);
   hmean_Hino_135_10ns_ZeeZSM_stack_4->SetBinError(3,1724.284);
   hmean_Hino_135_10ns_ZeeZSM_stack_4->SetBinError(4,3978.891);
   hmean_Hino_135_10ns_ZeeZSM_stack_4->SetBinError(5,3376.098);
   hmean_Hino_135_10ns_ZeeZSM_stack_4->SetBinError(6,1912.823);
   hmean_Hino_135_10ns_ZeeZSM_stack_4->SetBinError(7,20.18882);
   hmean_Hino_135_10ns_ZeeZSM_stack_4->SetBinError(8,262.2686);
   hmean_Hino_135_10ns_ZeeZSM_stack_4->SetBinError(9,2490.366);
   hmean_Hino_135_10ns_ZeeZSM_stack_4->SetEntries(9);

   ci = TColor::GetColor("#0099ff");
   hmean_Hino_135_10ns_ZeeZSM_stack_4->SetLineColor(ci);

   ci = TColor::GetColor("#0099ff");
   hmean_Hino_135_10ns_ZeeZSM_stack_4->SetMarkerColor(ci);
   hmean_Hino_135_10ns_ZeeZSM_stack_4->SetMarkerStyle(26);
   hmean_Hino_135_10ns_ZeeZSM_stack_4->SetMarkerSize(2);
   hmean_Hino_135_10ns_ZeeZSM_stack_4->GetXaxis()->SetLabelFont(42);
   hmean_Hino_135_10ns_ZeeZSM_stack_4->GetXaxis()->SetTitleOffset(1);
   hmean_Hino_135_10ns_ZeeZSM_stack_4->GetXaxis()->SetTitleFont(42);
   hmean_Hino_135_10ns_ZeeZSM_stack_4->GetYaxis()->SetLabelFont(42);
   hmean_Hino_135_10ns_ZeeZSM_stack_4->GetYaxis()->SetTitleFont(42);
   hmean_Hino_135_10ns_ZeeZSM_stack_4->GetZaxis()->SetLabelFont(42);
   hmean_Hino_135_10ns_ZeeZSM_stack_4->GetZaxis()->SetTitleOffset(1);
   hmean_Hino_135_10ns_ZeeZSM_stack_4->GetZaxis()->SetTitleFont(42);
   hs_m->Add(hmean_Hino_135_10ns_ZeeZSM_stack_4,"");
   Double_t xAxis6[10] = {0, 20, 40, 60, 80, 100, 120, 300, 500, 2000}; 
   
   TH1F *hmean_Hino_325_2ns_HyyHSM_stack_5 = new TH1F("hmean_Hino_325_2ns_HyyHSM_stack_5","",9, xAxis6);
   hmean_Hino_325_2ns_HyyHSM_stack_5->SetBinContent(1,-36.68216);
   hmean_Hino_325_2ns_HyyHSM_stack_5->SetBinContent(2,1.540303);
   hmean_Hino_325_2ns_HyyHSM_stack_5->SetBinContent(3,12.86671);
   hmean_Hino_325_2ns_HyyHSM_stack_5->SetBinContent(4,-88.32486);
   hmean_Hino_325_2ns_HyyHSM_stack_5->SetBinContent(5,30.02132);
   hmean_Hino_325_2ns_HyyHSM_stack_5->SetBinContent(6,117.9762);
   hmean_Hino_325_2ns_HyyHSM_stack_5->SetBinContent(7,1.133477);
   hmean_Hino_325_2ns_HyyHSM_stack_5->SetBinContent(8,11.76073);
   hmean_Hino_325_2ns_HyyHSM_stack_5->SetBinContent(9,141.9276);
   hmean_Hino_325_2ns_HyyHSM_stack_5->SetBinError(1,106.595);
   hmean_Hino_325_2ns_HyyHSM_stack_5->SetBinError(2,48.02097);
   hmean_Hino_325_2ns_HyyHSM_stack_5->SetBinError(3,159.1241);
   hmean_Hino_325_2ns_HyyHSM_stack_5->SetBinError(4,177.1929);
   hmean_Hino_325_2ns_HyyHSM_stack_5->SetBinError(5,78.68913);
   hmean_Hino_325_2ns_HyyHSM_stack_5->SetBinError(6,299.8126);
   hmean_Hino_325_2ns_HyyHSM_stack_5->SetBinError(7,8.013911);
   hmean_Hino_325_2ns_HyyHSM_stack_5->SetBinError(8,19.49306);
   hmean_Hino_325_2ns_HyyHSM_stack_5->SetBinError(9,307.6325);
   hmean_Hino_325_2ns_HyyHSM_stack_5->SetEntries(9);

   ci = TColor::GetColor("#ff9900");
   hmean_Hino_325_2ns_HyyHSM_stack_5->SetLineColor(ci);

   ci = TColor::GetColor("#ff9900");
   hmean_Hino_325_2ns_HyyHSM_stack_5->SetMarkerColor(ci);
   hmean_Hino_325_2ns_HyyHSM_stack_5->SetMarkerStyle(26);
   hmean_Hino_325_2ns_HyyHSM_stack_5->SetMarkerSize(2);
   hmean_Hino_325_2ns_HyyHSM_stack_5->GetXaxis()->SetLabelFont(42);
   hmean_Hino_325_2ns_HyyHSM_stack_5->GetXaxis()->SetTitleOffset(1);
   hmean_Hino_325_2ns_HyyHSM_stack_5->GetXaxis()->SetTitleFont(42);
   hmean_Hino_325_2ns_HyyHSM_stack_5->GetYaxis()->SetLabelFont(42);
   hmean_Hino_325_2ns_HyyHSM_stack_5->GetYaxis()->SetTitleFont(42);
   hmean_Hino_325_2ns_HyyHSM_stack_5->GetZaxis()->SetLabelFont(42);
   hmean_Hino_325_2ns_HyyHSM_stack_5->GetZaxis()->SetTitleOffset(1);
   hmean_Hino_325_2ns_HyyHSM_stack_5->GetZaxis()->SetTitleFont(42);
   hs_m->Add(hmean_Hino_325_2ns_HyyHSM_stack_5,"");
   Double_t xAxis7[10] = {0, 20, 40, 60, 80, 100, 120, 300, 500, 2000}; 
   
   TH1F *hmean_Hino_625_2ns_ZeeZSM_stack_6 = new TH1F("hmean_Hino_625_2ns_ZeeZSM_stack_6","",9, xAxis7);
   hmean_Hino_625_2ns_ZeeZSM_stack_6->SetBinContent(1,-432.1705);
   hmean_Hino_625_2ns_ZeeZSM_stack_6->SetBinContent(2,-26.47761);
   hmean_Hino_625_2ns_ZeeZSM_stack_6->SetBinContent(3,-90.72552);
   hmean_Hino_625_2ns_ZeeZSM_stack_6->SetBinContent(4,-80.87433);
   hmean_Hino_625_2ns_ZeeZSM_stack_6->SetBinContent(5,-57.13115);
   hmean_Hino_625_2ns_ZeeZSM_stack_6->SetBinContent(6,20.36743);
   hmean_Hino_625_2ns_ZeeZSM_stack_6->SetBinContent(7,2.683105);
   hmean_Hino_625_2ns_ZeeZSM_stack_6->SetBinContent(8,-36.97381);
   hmean_Hino_625_2ns_ZeeZSM_stack_6->SetBinContent(9,-97.2748);
   hmean_Hino_625_2ns_ZeeZSM_stack_6->SetBinError(1,5971.367);
   hmean_Hino_625_2ns_ZeeZSM_stack_6->SetBinError(2,204.0347);
   hmean_Hino_625_2ns_ZeeZSM_stack_6->SetBinError(3,219.366);
   hmean_Hino_625_2ns_ZeeZSM_stack_6->SetBinError(4,486.7974);
   hmean_Hino_625_2ns_ZeeZSM_stack_6->SetBinError(5,381.9382);
   hmean_Hino_625_2ns_ZeeZSM_stack_6->SetBinError(6,376.394);
   hmean_Hino_625_2ns_ZeeZSM_stack_6->SetBinError(7,11.42846);
   hmean_Hino_625_2ns_ZeeZSM_stack_6->SetBinError(8,51.57036);
   hmean_Hino_625_2ns_ZeeZSM_stack_6->SetBinError(9,141.0084);
   hmean_Hino_625_2ns_ZeeZSM_stack_6->SetEntries(9);

   ci = TColor::GetColor("#9966ff");
   hmean_Hino_625_2ns_ZeeZSM_stack_6->SetLineColor(ci);

   ci = TColor::GetColor("#9966ff");
   hmean_Hino_625_2ns_ZeeZSM_stack_6->SetMarkerColor(ci);
   hmean_Hino_625_2ns_ZeeZSM_stack_6->SetMarkerStyle(26);
   hmean_Hino_625_2ns_ZeeZSM_stack_6->SetMarkerSize(2);
   hmean_Hino_625_2ns_ZeeZSM_stack_6->GetXaxis()->SetLabelFont(42);
   hmean_Hino_625_2ns_ZeeZSM_stack_6->GetXaxis()->SetTitleOffset(1);
   hmean_Hino_625_2ns_ZeeZSM_stack_6->GetXaxis()->SetTitleFont(42);
   hmean_Hino_625_2ns_ZeeZSM_stack_6->GetYaxis()->SetLabelFont(42);
   hmean_Hino_625_2ns_ZeeZSM_stack_6->GetYaxis()->SetTitleFont(42);
   hmean_Hino_625_2ns_ZeeZSM_stack_6->GetZaxis()->SetLabelFont(42);
   hmean_Hino_625_2ns_ZeeZSM_stack_6->GetZaxis()->SetTitleOffset(1);
   hmean_Hino_625_2ns_ZeeZSM_stack_6->GetZaxis()->SetTitleFont(42);
   hs_m->Add(hmean_Hino_625_2ns_ZeeZSM_stack_6,"");
   Double_t xAxis8[10] = {0, 20, 40, 60, 80, 100, 120, 300, 500, 2000}; 
   
   TH1F *hmean_data_skimmed_stack_7 = new TH1F("hmean_data_skimmed_stack_7","",9, xAxis8);
   hmean_data_skimmed_stack_7->SetBinContent(1,-2.239078);
   hmean_data_skimmed_stack_7->SetBinContent(2,-8.123496);
   hmean_data_skimmed_stack_7->SetBinContent(3,-14.69265);
   hmean_data_skimmed_stack_7->SetBinContent(4,-18.39138);
   hmean_data_skimmed_stack_7->SetBinContent(5,-22.64304);
   hmean_data_skimmed_stack_7->SetBinContent(6,-31.59721);
   hmean_data_skimmed_stack_7->SetBinContent(7,-34.37787);
   hmean_data_skimmed_stack_7->SetBinError(1,0.2873326);
   hmean_data_skimmed_stack_7->SetBinError(2,0.3403423);
   hmean_data_skimmed_stack_7->SetBinError(3,0.4703144);
   hmean_data_skimmed_stack_7->SetBinError(4,0.760187);
   hmean_data_skimmed_stack_7->SetBinError(5,1.48342);
   hmean_data_skimmed_stack_7->SetBinError(6,2.991224);
   hmean_data_skimmed_stack_7->SetBinError(7,5.540205);
   hmean_data_skimmed_stack_7->SetEntries(9);

   ci = TColor::GetColor("#ff6666");
   hmean_data_skimmed_stack_7->SetLineColor(ci);

   ci = TColor::GetColor("#ff6666");
   hmean_data_skimmed_stack_7->SetMarkerColor(ci);
   hmean_data_skimmed_stack_7->SetMarkerStyle(20);
   hmean_data_skimmed_stack_7->SetMarkerSize(2);
   hmean_data_skimmed_stack_7->GetXaxis()->SetLabelFont(42);
   hmean_data_skimmed_stack_7->GetXaxis()->SetTitleOffset(1);
   hmean_data_skimmed_stack_7->GetXaxis()->SetTitleFont(42);
   hmean_data_skimmed_stack_7->GetYaxis()->SetLabelFont(42);
   hmean_data_skimmed_stack_7->GetYaxis()->SetTitleFont(42);
   hmean_data_skimmed_stack_7->GetZaxis()->SetLabelFont(42);
   hmean_data_skimmed_stack_7->GetZaxis()->SetTitleOffset(1);
   hmean_data_skimmed_stack_7->GetZaxis()->SetTitleFont(42);
   hs_m->Add(hmean_data_skimmed_stack_7,"");
   Double_t xAxis9[10] = {0, 20, 40, 60, 80, 100, 120, 300, 500, 2000}; 
   
   TH1F *hmean_Zee_truZ_stack_8 = new TH1F("hmean_Zee_truZ_stack_8","",9, xAxis9);
   hmean_Zee_truZ_stack_8->SetBinContent(1,-1.88591);
   hmean_Zee_truZ_stack_8->SetBinContent(2,-8.363108);
   hmean_Zee_truZ_stack_8->SetBinContent(3,-14.67068);
   hmean_Zee_truZ_stack_8->SetBinContent(4,-19.33361);
   hmean_Zee_truZ_stack_8->SetBinContent(5,-20.65642);
   hmean_Zee_truZ_stack_8->SetBinContent(6,-30.82584);
   hmean_Zee_truZ_stack_8->SetBinContent(7,-19.39207);
   hmean_Zee_truZ_stack_8->SetBinError(1,0.5725419);
   hmean_Zee_truZ_stack_8->SetBinError(2,0.6567366);
   hmean_Zee_truZ_stack_8->SetBinError(3,0.8137838);
   hmean_Zee_truZ_stack_8->SetBinError(4,1.2035);
   hmean_Zee_truZ_stack_8->SetBinError(5,1.923899);
   hmean_Zee_truZ_stack_8->SetBinError(6,3.948596);
   hmean_Zee_truZ_stack_8->SetBinError(7,9.798736);
   hmean_Zee_truZ_stack_8->SetEntries(9);
   hmean_Zee_truZ_stack_8->SetMarkerStyle(20);
   hmean_Zee_truZ_stack_8->SetMarkerSize(2);
   hmean_Zee_truZ_stack_8->GetXaxis()->SetLabelFont(42);
   hmean_Zee_truZ_stack_8->GetXaxis()->SetTitleOffset(1);
   hmean_Zee_truZ_stack_8->GetXaxis()->SetTitleFont(42);
   hmean_Zee_truZ_stack_8->GetYaxis()->SetLabelFont(42);
   hmean_Zee_truZ_stack_8->GetYaxis()->SetTitleFont(42);
   hmean_Zee_truZ_stack_8->GetZaxis()->SetLabelFont(42);
   hmean_Zee_truZ_stack_8->GetZaxis()->SetTitleOffset(1);
   hmean_Zee_truZ_stack_8->GetZaxis()->SetTitleFont(42);
   hs_m->Add(hmean_Zee_truZ_stack_8,"");
   hs_m->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("hmean_Hino_135_2ns_HyyHSM_stack_1","Hino_135_2ns_HyyHSM","lpf");
   entry->SetFillStyle(1001);

   ci = TColor::GetColor("#ff66ff");
   entry->SetLineColor(ci);
   entry->SetLineStyle(1);
   entry->SetLineWidth(1);

   ci = TColor::GetColor("#ff66ff");
   entry->SetMarkerColor(ci);
   entry->SetMarkerStyle(26);
   entry->SetMarkerSize(2);
   entry->SetTextFont(42);
   entry=leg->AddEntry("hmean_Hino_135_2ns_ZeeZSM_stack_2","Hino_135_2ns_ZeeZSM","lpf");
   entry->SetFillStyle(1001);

   ci = TColor::GetColor("#009900");
   entry->SetLineColor(ci);
   entry->SetLineStyle(1);
   entry->SetLineWidth(1);

   ci = TColor::GetColor("#009900");
   entry->SetMarkerColor(ci);
   entry->SetMarkerStyle(26);
   entry->SetMarkerSize(2);
   entry->SetTextFont(42);
   entry=leg->AddEntry("hmean_Hino_135_10ns_HyyHSM_stack_3","Hino_135_10ns_HyyHSM","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("hmean_Hino_135_10ns_ZeeZSM_stack_4","Hino_135_10ns_ZeeZSM","lpf");
   entry->SetFillStyle(1001);

   ci = TColor::GetColor("#0099ff");
   entry->SetLineColor(ci);
   entry->SetLineStyle(1);
   entry->SetLineWidth(1);

   ci = TColor::GetColor("#0099ff");
   entry->SetMarkerColor(ci);
   entry->SetMarkerStyle(26);
   entry->SetMarkerSize(2);
   entry->SetTextFont(42);
   entry=leg->AddEntry("hmean_Hino_325_2ns_HyyHSM_stack_5","Hino_325_2ns_HyyHSM","lpf");
   entry->SetFillStyle(1001);

   ci = TColor::GetColor("#ff9900");
   entry->SetLineColor(ci);
   entry->SetLineStyle(1);
   entry->SetLineWidth(1);

   ci = TColor::GetColor("#ff9900");
   entry->SetMarkerColor(ci);
   entry->SetMarkerStyle(26);
   entry->SetMarkerSize(2);
   entry->SetTextFont(42);
   entry=leg->AddEntry("hmean_Hino_625_2ns_ZeeZSM_stack_6","Hino_625_2ns_ZeeZSM","lpf");
   entry->SetFillStyle(1001);

   ci = TColor::GetColor("#9966ff");
   entry->SetLineColor(ci);
   entry->SetLineStyle(1);
   entry->SetLineWidth(1);

   ci = TColor::GetColor("#9966ff");
   entry->SetMarkerColor(ci);
   entry->SetMarkerStyle(26);
   entry->SetMarkerSize(2);
   entry->SetTextFont(42);
   entry=leg->AddEntry("hmean_data_skimmed_stack_7","data_skimmed","lpf");
   entry->SetFillStyle(1001);

   ci = TColor::GetColor("#ff6666");
   entry->SetLineColor(ci);
   entry->SetLineStyle(1);
   entry->SetLineWidth(1);

   ci = TColor::GetColor("#ff6666");
   entry->SetMarkerColor(ci);
   entry->SetMarkerStyle(20);
   entry->SetMarkerSize(2);
   entry->SetTextFont(42);
   entry=leg->AddEntry("hmean_Zee_truZ_stack_8","Zee_truZ","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);
}