{
//=========Macro generated from canvas: Systematics/Systematics
//=========  (Thu Feb  4 15:51:00 2016) by ROOT version5.27/06b
   TCanvas *Systematics = new TCanvas("Systematics", "Systematics",0,0,900,900);
   gStyle->SetOptFit(1);
   gStyle->SetOptStat(0);
   gStyle->SetOptTitle(0);
   Systematics->SetHighLightColor(2);
   Systematics->Range(0,0,1,1);
   Systematics->SetFillColor(0);
   Systematics->SetBorderMode(0);
   Systematics->SetBorderSize(2);
   Systematics->SetTickx(1);
   Systematics->SetTicky(1);
   Systematics->SetLeftMargin(0.25);
   Systematics->SetRightMargin(0.05);
   Systematics->SetTopMargin(0.07);
   Systematics->SetBottomMargin(0.14);
   Systematics->SetFrameFillStyle(0);
   Systematics->SetFrameBorderMode(0);
  
// ------------>Primitives in pad: Absolute_Values
   TPad *Absolute_Values = new TPad("Absolute_Values", "Absolute_Values",0,0.4,1,1);
   Absolute_Values->Draw();
   Absolute_Values->cd();
   Absolute_Values->Range(0,0,1,1);
   Absolute_Values->SetFillColor(0);
   Absolute_Values->SetBorderMode(0);
   Absolute_Values->SetBorderSize(2);
   Absolute_Values->SetLogy();
   Absolute_Values->SetTickx(1);
   Absolute_Values->SetTicky(1);
   Absolute_Values->SetLeftMargin(0.2);
   Absolute_Values->SetRightMargin(0.05);
   Absolute_Values->SetTopMargin(0.05);
   Absolute_Values->SetBottomMargin(0);
   Absolute_Values->SetFrameFillStyle(0);
   Absolute_Values->SetFrameBorderMode(0);
   
   TGraphAsymmErrors *grae = new TGraphAsymmErrors(24);
   grae->SetName("Graph");
   grae->SetTitle("Graph");

   Int_t ci;   // for color index setting
   ci = TColor::GetColor("#c0c0c0");
   grae->SetFillColor(ci);
   grae->SetLineWidth(2);
   grae->SetMarkerStyle(20);
   grae->SetPoint(0,150,0.486747);
   grae->SetPointError(0,18.75,18.75,0.04176004,0.4699325);
   grae->SetPoint(1,187.5,0.486747);
   grae->SetPointError(1,18.75,18.75,0.04176004,0.4699325);
   grae->SetPoint(2,187.5,0.3365233);
   grae->SetPointError(2,23.4375,23.4375,0.06815357,0.09651989);
   grae->SetPoint(3,234.375,0.3365233);
   grae->SetPointError(3,23.4375,23.4375,0.06815357,0.09651989);
   grae->SetPoint(4,234.375,0.2056066);
   grae->SetPointError(4,29.29688,29.29688,0.06580265,0.2591217);
   grae->SetPoint(5,292.9688,0.2056066);
   grae->SetPointError(5,29.29688,29.29688,0.06580265,0.2591217);
   grae->SetPoint(6,292.9688,0.1122084);
   grae->SetPointError(6,36.62109,36.62109,0.04776169,0.05706626);
   grae->SetPoint(7,366.2109,0.1122084);
   grae->SetPointError(7,36.62109,36.62109,0.04776169,0.05706626);
   grae->SetPoint(8,366.2109,0.05253398);
   grae->SetPointError(8,45.77637,45.77637,0.02635883,0.02338781);
   grae->SetPoint(9,457.7637,0.05253398);
   grae->SetPointError(9,45.77637,45.77637,0.02635883,0.02338781);
   grae->SetPoint(10,457.7637,0.02196398);
   grae->SetPointError(10,57.22046,57.22046,0.01247788,0.01263618);
   grae->SetPoint(11,572.2046,0.02196398);
   grae->SetPointError(11,57.22046,57.22046,0.01247788,0.01263618);
   grae->SetPoint(12,572.2046,0.007710425);
   grae->SetPointError(12,71.52557,71.52557,0.004726953,0.007340367);
   grae->SetPoint(13,715.2557,0.007710425);
   grae->SetPointError(13,71.52557,71.52557,0.004726953,0.007340367);
   grae->SetPoint(14,715.2557,0.00235306);
   grae->SetPointError(14,89.40697,89.40697,0.001625346,0.003116022);
   grae->SetPoint(15,894.0697,0.00235306);
   grae->SetPointError(15,89.40697,89.40697,0.001625346,0.003116022);
   grae->SetPoint(16,894.0697,0.000779402);
   grae->SetPointError(16,111.7587,111.7587,0.0005828685,0.0009863261);
   grae->SetPoint(17,1117.587,0.000779402);
   grae->SetPointError(17,111.7587,111.7587,0.0005828685,0.0009863261);
   grae->SetPoint(18,1117.587,0.0002567834);
   grae->SetPointError(18,139.6984,139.6984,0.0001994997,0.0003785163);
   grae->SetPoint(19,1396.984,0.0002567834);
   grae->SetPointError(19,139.6984,139.6984,0.0001994997,0.0003785163);
   grae->SetPoint(20,1396.984,6.096639e-05);
   grae->SetPointError(20,174.623,174.623,4.874606e-05,9.503296e-05);
   grae->SetPoint(21,1746.23,6.096639e-05);
   grae->SetPointError(21,174.623,174.623,4.874606e-05,9.503296e-05);
   grae->SetPoint(22,1746.23,9.085633e-06);
   grae->SetPointError(22,218.2787,218.2787,8.870719e-06,3.48777e-05);
   grae->SetPoint(23,2182.787,9.085633e-06);
   grae->SetPointError(23,218.2787,218.2787,8.870719e-06,3.48777e-05);
   
   TH1F *Graph1 = new TH1F("Graph1","Graph",100,0,2628.048);
   Graph1->SetMinimum(7.5e-06);
   Graph1->SetMaximum(0.5354217);
   Graph1->SetDirectory(0);
   Graph1->SetStats(0);
   Graph1->SetLineStyle(0);
   Graph1->SetMarkerStyle(20);
   Graph1->GetXaxis()->SetTitle("E (GeV)");
   Graph1->GetXaxis()->SetRange(13,67);
   Graph1->GetXaxis()->SetLabelFont(42);
   Graph1->GetXaxis()->SetLabelOffset(0.007);
   Graph1->GetXaxis()->SetLabelSize(0.05);
   Graph1->GetXaxis()->SetTitleSize(0.06);
   Graph1->GetXaxis()->SetTitleOffset(1.1);
   Graph1->GetXaxis()->SetTitleFont(42);
   Graph1->GetYaxis()->SetTitle("#frac{d#sigma}{dE} [mb/GeV]");
   Graph1->GetYaxis()->SetLabelFont(42);
   Graph1->GetYaxis()->SetLabelOffset(0.007);
   Graph1->GetYaxis()->SetLabelSize(0.05);
   Graph1->GetYaxis()->SetTitleSize(0.06);
   Graph1->GetYaxis()->SetTitleOffset(1.25);
   Graph1->GetYaxis()->SetTitleFont(42);
   Graph1->GetZaxis()->SetLabelFont(42);
   Graph1->GetZaxis()->SetLabelOffset(0.007);
   Graph1->GetZaxis()->SetLabelSize(0.05);
   Graph1->GetZaxis()->SetTitleSize(0.06);
   Graph1->GetZaxis()->SetTitleFont(42);
   grae->SetHistogram(Graph1);
   
   grae->Draw("ae3");
   
   grae = new TGraphAsymmErrors(24);
   grae->SetName("Graph");
   grae->SetTitle("Graph");
   grae->SetFillColor(13);
   grae->SetLineWidth(2);
   grae->SetMarkerStyle(20);
   grae->SetPoint(0,150,0.486747);
   grae->SetPointError(0,18.75,18.75,0.01084962,0.009572623);
   grae->SetPoint(1,187.5,0.486747);
   grae->SetPointError(1,18.75,18.75,0.01084962,0.009572623);
   grae->SetPoint(2,187.5,0.3365233);
   grae->SetPointError(2,23.4375,23.4375,0.04069026,0.03197586);
   grae->SetPoint(3,234.375,0.3365233);
   grae->SetPointError(3,23.4375,23.4375,0.04069026,0.03197586);
   grae->SetPoint(4,234.375,0.2056066);
   grae->SetPointError(4,29.29688,29.29688,0.04309834,0.04511864);
   grae->SetPoint(5,292.9688,0.2056066);
   grae->SetPointError(5,29.29688,29.29688,0.04309834,0.04511864);
   grae->SetPoint(6,292.9688,0.1122084);
   grae->SetPointError(6,36.62109,36.62109,0.0305342,0.03590714);
   grae->SetPoint(7,366.2109,0.1122084);
   grae->SetPointError(7,36.62109,36.62109,0.0305342,0.03590714);
   grae->SetPoint(8,366.2109,0.05253398);
   grae->SetPointError(8,45.77637,45.77637,0.01540801,0.02129915);
   grae->SetPoint(9,457.7637,0.05253398);
   grae->SetPointError(9,45.77637,45.77637,0.01540801,0.02129915);
   grae->SetPoint(10,457.7637,0.02196398);
   grae->SetPointError(10,57.22046,57.22046,0.006671538,0.009961976);
   grae->SetPoint(11,572.2046,0.02196398);
   grae->SetPointError(11,57.22046,57.22046,0.006671538,0.009961976);
   grae->SetPoint(12,572.2046,0.007710425);
   grae->SetPointError(12,71.52557,71.52557,0.002208281,0.00375652);
   grae->SetPoint(13,715.2557,0.007710425);
   grae->SetPointError(13,71.52557,71.52557,0.002208281,0.00375652);
   grae->SetPoint(14,715.2557,0.00235306);
   grae->SetPointError(14,89.40697,89.40697,0.0008437036,0.0012594);
   grae->SetPoint(15,894.0697,0.00235306);
   grae->SetPointError(15,89.40697,89.40697,0.0008437036,0.0012594);
   grae->SetPoint(16,894.0697,0.000779402);
   grae->SetPointError(16,111.7587,111.7587,0.0002888264,0.0003669074);
   grae->SetPoint(17,1117.587,0.000779402);
   grae->SetPointError(17,111.7587,111.7587,0.0002888264,0.0003669074);
   grae->SetPoint(18,1117.587,0.0002567834);
   grae->SetPointError(18,139.6984,139.6984,9.231423e-05,8.613581e-05);
   grae->SetPoint(19,1396.984,0.0002567834);
   grae->SetPointError(19,139.6984,139.6984,9.231423e-05,8.613581e-05);
   grae->SetPoint(20,1396.984,6.096639e-05);
   grae->SetPointError(20,174.623,174.623,2.446374e-05,4.417039e-05);
   grae->SetPoint(21,1746.23,6.096639e-05);
   grae->SetPointError(21,174.623,174.623,2.446374e-05,4.417039e-05);
   grae->SetPoint(22,1746.23,9.085633e-06);
   grae->SetPointError(22,218.2787,218.2787,6.591629e-06,1.376454e-05);
   grae->SetPoint(23,2182.787,9.085633e-06);
   grae->SetPointError(23,218.2787,218.2787,6.591629e-06,1.376454e-05);
   
   TH1F *Graph2 = new TH1F("Graph2","Graph",100,0,2628.048);
   Graph2->SetMinimum(7.5e-06);
   Graph2->SetMaximum(0.5354217);
   Graph2->SetDirectory(0);
   Graph2->SetStats(0);
   Graph2->SetLineStyle(0);
   Graph2->SetMarkerStyle(20);
   Graph2->GetXaxis()->SetTitle("E (GeV)");
   Graph2->GetXaxis()->SetRange(12,67);
   Graph2->GetXaxis()->SetLabelFont(42);
   Graph2->GetXaxis()->SetLabelOffset(0.007);
   Graph2->GetXaxis()->SetLabelSize(0.05);
   Graph2->GetXaxis()->SetTitleSize(0.06);
   Graph2->GetXaxis()->SetTitleOffset(1.1);
   Graph2->GetXaxis()->SetTitleFont(42);
   Graph2->GetYaxis()->SetTitle("#frac{dN}{dE}");
   Graph2->GetYaxis()->SetLabelFont(42);
   Graph2->GetYaxis()->SetLabelOffset(0.007);
   Graph2->GetYaxis()->SetLabelSize(0.05);
   Graph2->GetYaxis()->SetTitleSize(0.06);
   Graph2->GetYaxis()->SetTitleOffset(1.25);
   Graph2->GetYaxis()->SetTitleFont(42);
   Graph2->GetZaxis()->SetLabelFont(42);
   Graph2->GetZaxis()->SetLabelOffset(0.007);
   Graph2->GetZaxis()->SetLabelSize(0.05);
   Graph2->GetZaxis()->SetTitleSize(0.06);
   Graph2->GetZaxis()->SetTitleFont(42);
   grae->SetHistogram(Graph2);
   
   grae->Draw("pe3");
   
   grae = new TGraphAsymmErrors(12);
   grae->SetName("Graph");
   grae->SetTitle("Graph");
   grae->SetFillColor(1);
   grae->SetLineWidth(2);
   grae->SetMarkerStyle(20);
   grae->SetMarkerSize(1.25);
   grae->SetPoint(0,168.3176,0.486747);
   grae->SetPointError(0,18.31763,19.18237,0.0007333285,0.0007333285);
   grae->SetPoint(1,210.2622,0.3365233);
   grae->SetPointError(1,22.76216,24.11284,0.0006157064,0.0006157064);
   grae->SetPoint(2,262.6172,0.2056066);
   grae->SetPointError(2,28.24225,30.3515,0.0003555667,0.0003555667);
   grae->SetPoint(3,327.9434,0.1122084);
   grae->SetPointError(3,34.97468,38.26751,0.0002332533,0.0002332533);
   grae->SetPoint(4,409.4183,0.05253398);
   grae->SetPointError(4,43.20734,48.3454,0.0001374996,0.0001374996);
   grae->SetPoint(5,510.9785,0.02196398);
   grae->SetPointError(5,53.21482,61.2261,7.665124e-05,7.665124e-05);
   grae->SetPoint(6,637.4918,0.007710425);
   grae->SetPointError(6,65.28724,77.7639,4.191671e-05,4.191671e-05);
   grae->SetPoint(7,794.9645,0.00235306);
   grae->SetPointError(7,79.70878,99.10515,2.068978e-05,2.068978e-05);
   grae->SetPoint(8,990.7922,0.000779402);
   grae->SetPointError(8,96.72253,126.7949,1.202693e-05,1.202693e-05);
   grae->SetPoint(9,1242.253,0.0002567834);
   grae->SetPointError(9,124.666,154.7308,5.374283e-06,5.374283e-06);
   grae->SetPoint(10,1548.298,6.096639e-05);
   grae->SetPointError(10,151.3137,197.9323,3.221391e-06,3.221391e-06);
   grae->SetPoint(11,1928.509,9.085633e-06);
   grae->SetPointError(11,182.279,254.2784,4.729672e-07,4.729672e-07);
   grae->Draw("ep");
   Double_t xAxis9[13] = {150, 187.5, 234.375, 292.9688, 366.2109, 457.7637, 572.2046, 715.2557, 894.0697, 1117.587, 1396.984, 1746.23, 2182.787}; 
   
   TH1D *hGenJet_energy_noCuts = new TH1D("hGenJet_energy_noCuts","GenJet energy distribution",12, xAxis9);
   hGenJet_energy_noCuts->SetBinContent(0,0.4588835);
   hGenJet_energy_noCuts->SetBinContent(1,0.4588835);
   hGenJet_energy_noCuts->SetBinContent(2,0.2834228);
   hGenJet_energy_noCuts->SetBinContent(3,0.1599712);
   hGenJet_energy_noCuts->SetBinContent(4,0.08396457);
   hGenJet_energy_noCuts->SetBinContent(5,0.04049441);
   hGenJet_energy_noCuts->SetBinContent(6,0.01890777);
   hGenJet_energy_noCuts->SetBinContent(7,0.008182113);
   hGenJet_energy_noCuts->SetBinContent(8,0.003238356);
   hGenJet_energy_noCuts->SetBinContent(9,0.00123091);
   hGenJet_energy_noCuts->SetBinContent(10,0.000404966);
   hGenJet_energy_noCuts->SetBinContent(11,0.0001026833);
   hGenJet_energy_noCuts->SetBinContent(12,1.46982e-05);
   hGenJet_energy_noCuts->SetBinContent(13,0.001212026);
   hGenJet_energy_noCuts->SetBinError(0,16.20901);
   hGenJet_energy_noCuts->SetBinError(1,0.1106205);
   hGenJet_energy_noCuts->SetBinError(2,0.0777583);
   hGenJet_energy_noCuts->SetBinError(3,0.05225107);
   hGenJet_energy_noCuts->SetBinError(4,0.03385847);
   hGenJet_energy_noCuts->SetBinError(5,0.02103109);
   hGenJet_energy_noCuts->SetBinError(6,0.01285374);
   hGenJet_energy_noCuts->SetBinError(7,0.007562877);
   hGenJet_energy_noCuts->SetBinError(8,0.004255608);
   hGenJet_energy_noCuts->SetBinError(9,0.0023467);
   hGenJet_energy_noCuts->SetBinError(10,0.001203923);
   hGenJet_energy_noCuts->SetBinError(11,0.0005422307);
   hGenJet_energy_noCuts->SetBinError(12,0.0001834896);
   hGenJet_energy_noCuts->SetBinError(13,0.03481416);
   hGenJet_energy_noCuts->SetEntries(4443496);
   hGenJet_energy_noCuts->SetStats(0);
   hGenJet_energy_noCuts->SetLineColor(2);
   hGenJet_energy_noCuts->SetLineWidth(3);
   hGenJet_energy_noCuts->SetMarkerStyle(20);
   hGenJet_energy_noCuts->GetXaxis()->SetTitle("E_{det} (GeV)");
   hGenJet_energy_noCuts->GetXaxis()->SetLabelFont(42);
   hGenJet_energy_noCuts->GetXaxis()->SetLabelOffset(0.007);
   hGenJet_energy_noCuts->GetXaxis()->SetLabelSize(0.05);
   hGenJet_energy_noCuts->GetXaxis()->SetTitleSize(0.06);
   hGenJet_energy_noCuts->GetXaxis()->SetTitleOffset(0.9);
   hGenJet_energy_noCuts->GetXaxis()->SetTitleFont(42);
   hGenJet_energy_noCuts->GetYaxis()->SetTitle("#frac{dN}{dE}");
   hGenJet_energy_noCuts->GetYaxis()->SetLabelFont(42);
   hGenJet_energy_noCuts->GetYaxis()->SetLabelOffset(0.007);
   hGenJet_energy_noCuts->GetYaxis()->SetLabelSize(0.05);
   hGenJet_energy_noCuts->GetYaxis()->SetTitleSize(0.06);
   hGenJet_energy_noCuts->GetYaxis()->SetTitleOffset(1.05);
   hGenJet_energy_noCuts->GetYaxis()->SetTitleFont(42);
   hGenJet_energy_noCuts->GetZaxis()->SetLabelFont(42);
   hGenJet_energy_noCuts->GetZaxis()->SetLabelOffset(0.007);
   hGenJet_energy_noCuts->GetZaxis()->SetLabelSize(0.05);
   hGenJet_energy_noCuts->GetZaxis()->SetTitleSize(0.06);
   hGenJet_energy_noCuts->GetZaxis()->SetTitleFont(42);
   hGenJet_energy_noCuts->Draw("histsame");
   Double_t xAxis10[13] = {150, 187.5, 234.375, 292.9688, 366.2109, 457.7637, 572.2046, 715.2557, 894.0697, 1117.587, 1396.984, 1746.23, 2182.787}; 
   
   TH1D *hGenJet_energy_noCuts = new TH1D("hGenJet_energy_noCuts","GenJet energy distribution",12, xAxis10);
   hGenJet_energy_noCuts->SetBinContent(0,0.5103658);
   hGenJet_energy_noCuts->SetBinContent(1,0.5103658);
   hGenJet_energy_noCuts->SetBinContent(2,0.3089233);
   hGenJet_energy_noCuts->SetBinContent(3,0.1722346);
   hGenJet_energy_noCuts->SetBinContent(4,0.08951848);
   hGenJet_energy_noCuts->SetBinContent(5,0.04339886);
   hGenJet_energy_noCuts->SetBinContent(6,0.02009845);
   hGenJet_energy_noCuts->SetBinContent(7,0.008757462);
   hGenJet_energy_noCuts->SetBinContent(8,0.003638224);
   hGenJet_energy_noCuts->SetBinContent(9,0.001334577);
   hGenJet_energy_noCuts->SetBinContent(10,0.0004235205);
   hGenJet_energy_noCuts->SetBinContent(11,0.0001171166);
   hGenJet_energy_noCuts->SetBinContent(12,2.136865e-05);
   hGenJet_energy_noCuts->SetBinContent(13,0.001291658);
   hGenJet_energy_noCuts->SetBinError(0,17.10806);
   hGenJet_energy_noCuts->SetBinError(1,0.1166608);
   hGenJet_energy_noCuts->SetBinError(2,0.08118106);
   hGenJet_energy_noCuts->SetBinError(3,0.05421688);
   hGenJet_energy_noCuts->SetBinError(4,0.03496034);
   hGenJet_energy_noCuts->SetBinError(5,0.02177226);
   hGenJet_energy_noCuts->SetBinError(6,0.01325228);
   hGenJet_energy_noCuts->SetBinError(7,0.007824263);
   hGenJet_energy_noCuts->SetBinError(8,0.004510701);
   hGenJet_energy_noCuts->SetBinError(9,0.002443522);
   hGenJet_energy_noCuts->SetBinError(10,0.001231194);
   hGenJet_energy_noCuts->SetBinError(11,0.0005790867);
   hGenJet_energy_noCuts->SetBinError(12,0.0002212421);
   hGenJet_energy_noCuts->SetBinError(13,0.03593965);
   hGenJet_energy_noCuts->SetEntries(4899740);
   hGenJet_energy_noCuts->SetStats(0);

   ci = TColor::GetColor("#ff9933");
   hGenJet_energy_noCuts->SetLineColor(ci);
   hGenJet_energy_noCuts->SetLineStyle(2);
   hGenJet_energy_noCuts->SetLineWidth(3);
   hGenJet_energy_noCuts->SetMarkerStyle(20);
   hGenJet_energy_noCuts->GetXaxis()->SetTitle("E_{det} (GeV)");
   hGenJet_energy_noCuts->GetXaxis()->SetLabelFont(42);
   hGenJet_energy_noCuts->GetXaxis()->SetLabelOffset(0.007);
   hGenJet_energy_noCuts->GetXaxis()->SetLabelSize(0.05);
   hGenJet_energy_noCuts->GetXaxis()->SetTitleSize(0.06);
   hGenJet_energy_noCuts->GetXaxis()->SetTitleOffset(0.9);
   hGenJet_energy_noCuts->GetXaxis()->SetTitleFont(42);
   hGenJet_energy_noCuts->GetYaxis()->SetTitle("#frac{dN}{dE}");
   hGenJet_energy_noCuts->GetYaxis()->SetLabelFont(42);
   hGenJet_energy_noCuts->GetYaxis()->SetLabelOffset(0.007);
   hGenJet_energy_noCuts->GetYaxis()->SetLabelSize(0.05);
   hGenJet_energy_noCuts->GetYaxis()->SetTitleSize(0.06);
   hGenJet_energy_noCuts->GetYaxis()->SetTitleOffset(1.05);
   hGenJet_energy_noCuts->GetYaxis()->SetTitleFont(42);
   hGenJet_energy_noCuts->GetZaxis()->SetLabelFont(42);
   hGenJet_energy_noCuts->GetZaxis()->SetLabelOffset(0.007);
   hGenJet_energy_noCuts->GetZaxis()->SetLabelSize(0.05);
   hGenJet_energy_noCuts->GetZaxis()->SetTitleSize(0.06);
   hGenJet_energy_noCuts->GetZaxis()->SetTitleFont(42);
   hGenJet_energy_noCuts->Draw("histsame");
   Double_t xAxis11[13] = {150, 187.5, 234.375, 292.9688, 366.2109, 457.7637, 572.2046, 715.2557, 894.0697, 1117.587, 1396.984, 1746.23, 2182.787}; 
   
   TH1D *hGenJet_energy_noCuts = new TH1D("hGenJet_energy_noCuts","GenJet energy distribution",12, xAxis11);
   hGenJet_energy_noCuts->SetBinContent(0,0.6260885);
   hGenJet_energy_noCuts->SetBinContent(1,0.6260885);
   hGenJet_energy_noCuts->SetBinContent(2,0.4108947);
   hGenJet_energy_noCuts->SetBinContent(3,0.2493341);
   hGenJet_energy_noCuts->SetBinContent(4,0.1396299);
   hGenJet_energy_noCuts->SetBinContent(5,0.07161827);
   hGenJet_energy_noCuts->SetBinContent(6,0.03283741);
   hGenJet_energy_noCuts->SetBinContent(7,0.01350987);
   hGenJet_energy_noCuts->SetBinContent(8,0.004765999);
   hGenJet_energy_noCuts->SetBinContent(9,0.001455141);
   hGenJet_energy_noCuts->SetBinContent(10,0.0003507694);
   hGenJet_energy_noCuts->SetBinContent(11,6.034743e-05);
   hGenJet_energy_noCuts->SetBinContent(12,6.207164e-06);
   hGenJet_energy_noCuts->SetBinContent(13,0.0001505435);
   hGenJet_energy_noCuts->SetBinError(0,17.46973);
   hGenJet_energy_noCuts->SetBinError(1,0.1292118);
   hGenJet_energy_noCuts->SetBinError(2,0.09362561);
   hGenJet_energy_noCuts->SetBinError(3,0.06523268);
   hGenJet_energy_noCuts->SetBinError(4,0.04366249);
   hGenJet_energy_noCuts->SetBinError(5,0.02796896);
   hGenJet_energy_noCuts->SetBinError(6,0.01693923);
   hGenJet_energy_noCuts->SetBinError(7,0.009718068);
   hGenJet_energy_noCuts->SetBinError(8,0.005162693);
   hGenJet_energy_noCuts->SetBinError(9,0.002551507);
   hGenJet_energy_noCuts->SetBinError(10,0.00112047);
   hGenJet_energy_noCuts->SetBinError(11,0.0004156844);
   hGenJet_energy_noCuts->SetBinError(12,0.0001192411);
   hGenJet_energy_noCuts->SetBinError(13,0.01226962);
   hGenJet_energy_noCuts->SetEntries(5132262);
   hGenJet_energy_noCuts->SetStats(0);

   ci = TColor::GetColor("#0000ff");
   hGenJet_energy_noCuts->SetLineColor(ci);
   hGenJet_energy_noCuts->SetLineStyle(3);
   hGenJet_energy_noCuts->SetLineWidth(3);
   hGenJet_energy_noCuts->SetMarkerStyle(20);
   hGenJet_energy_noCuts->GetXaxis()->SetTitle("E_{det} (GeV)");
   hGenJet_energy_noCuts->GetXaxis()->SetLabelFont(42);
   hGenJet_energy_noCuts->GetXaxis()->SetLabelOffset(0.007);
   hGenJet_energy_noCuts->GetXaxis()->SetLabelSize(0.05);
   hGenJet_energy_noCuts->GetXaxis()->SetTitleSize(0.06);
   hGenJet_energy_noCuts->GetXaxis()->SetTitleOffset(0.9);
   hGenJet_energy_noCuts->GetXaxis()->SetTitleFont(42);
   hGenJet_energy_noCuts->GetYaxis()->SetTitle("#frac{dN}{dE}");
   hGenJet_energy_noCuts->GetYaxis()->SetLabelFont(42);
   hGenJet_energy_noCuts->GetYaxis()->SetLabelOffset(0.007);
   hGenJet_energy_noCuts->GetYaxis()->SetLabelSize(0.05);
   hGenJet_energy_noCuts->GetYaxis()->SetTitleSize(0.06);
   hGenJet_energy_noCuts->GetYaxis()->SetTitleOffset(1.05);
   hGenJet_energy_noCuts->GetYaxis()->SetTitleFont(42);
   hGenJet_energy_noCuts->GetZaxis()->SetLabelFont(42);
   hGenJet_energy_noCuts->GetZaxis()->SetLabelOffset(0.007);
   hGenJet_energy_noCuts->GetZaxis()->SetLabelSize(0.05);
   hGenJet_energy_noCuts->GetZaxis()->SetTitleSize(0.06);
   hGenJet_energy_noCuts->GetZaxis()->SetTitleFont(42);
   hGenJet_energy_noCuts->Draw("histsame");
   Double_t xAxis12[13] = {150, 187.5, 234.375, 292.9688, 366.2109, 457.7637, 572.2046, 715.2557, 894.0697, 1117.587, 1396.984, 1746.23, 2182.787}; 
   
   TH1D *hGenJet_Pythia84C_MPI_off = new TH1D("hGenJet_Pythia84C_MPI_off","Genjet energy Pythia84C, MPI off",12, xAxis12);
   hGenJet_Pythia84C_MPI_off->SetBinContent(0,0.3586622);
   hGenJet_Pythia84C_MPI_off->SetBinContent(1,0.3586622);
   hGenJet_Pythia84C_MPI_off->SetBinContent(2,0.1998119);
   hGenJet_Pythia84C_MPI_off->SetBinContent(3,0.09939794);
   hGenJet_Pythia84C_MPI_off->SetBinContent(4,0.04530314);
   hGenJet_Pythia84C_MPI_off->SetBinContent(5,0.0189721);
   hGenJet_Pythia84C_MPI_off->SetBinContent(6,0.007799408);
   hGenJet_Pythia84C_MPI_off->SetBinContent(7,0.003291711);
   hGenJet_Pythia84C_MPI_off->SetBinContent(8,0.001366796);
   hGenJet_Pythia84C_MPI_off->SetBinContent(9,0.0005820164);
   hGenJet_Pythia84C_MPI_off->SetBinContent(10,0.0002087848);
   hGenJet_Pythia84C_MPI_off->SetBinContent(11,6.480762e-05);
   hGenJet_Pythia84C_MPI_off->SetBinContent(12,1.242998e-05);
   hGenJet_Pythia84C_MPI_off->SetBinContent(13,0.0009282);
   hGenJet_Pythia84C_MPI_off->SetBinError(0,0.134095);
   hGenJet_Pythia84C_MPI_off->SetBinError(1,0.0008263733);
   hGenJet_Pythia84C_MPI_off->SetBinError(2,0.0005516823);
   hGenJet_Pythia84C_MPI_off->SetBinError(3,0.0003480264);
   hGenJet_Pythia84C_MPI_off->SetBinError(4,0.0002101515);
   hGenJet_Pythia84C_MPI_off->SetBinError(5,0.0001216385);
   hGenJet_Pythia84C_MPI_off->SetBinError(6,6.975723e-05);
   hGenJet_Pythia84C_MPI_off->SetBinError(7,4.053351e-05);
   hGenJet_Pythia84C_MPI_off->SetBinError(8,2.336148e-05);
   hGenJet_Pythia84C_MPI_off->SetBinError(9,1.363518e-05);
   hGenJet_Pythia84C_MPI_off->SetBinError(10,7.304456e-06);
   hGenJet_Pythia84C_MPI_off->SetBinError(11,3.639959e-06);
   hGenJet_Pythia84C_MPI_off->SetBinError(12,1.425816e-06);
   hGenJet_Pythia84C_MPI_off->SetBinError(13,0.0002574364);
   hGenJet_Pythia84C_MPI_off->SetEntries(4024689);

   ci = TColor::GetColor("#cc00cc");
   hGenJet_Pythia84C_MPI_off->SetLineColor(ci);
   hGenJet_Pythia84C_MPI_off->SetLineStyle(4);
   hGenJet_Pythia84C_MPI_off->SetLineWidth(3);
   hGenJet_Pythia84C_MPI_off->SetMarkerStyle(25);
   hGenJet_Pythia84C_MPI_off->GetXaxis()->SetTitle("E_{det} (GeV)");
   hGenJet_Pythia84C_MPI_off->GetXaxis()->SetLabelFont(42);
   hGenJet_Pythia84C_MPI_off->GetXaxis()->SetLabelSize(0.035);
   hGenJet_Pythia84C_MPI_off->GetXaxis()->SetTitleSize(0.035);
   hGenJet_Pythia84C_MPI_off->GetXaxis()->SetTitleFont(42);
   hGenJet_Pythia84C_MPI_off->GetYaxis()->SetTitle("#frac{dN}{dE}");
   hGenJet_Pythia84C_MPI_off->GetYaxis()->SetLabelFont(42);
   hGenJet_Pythia84C_MPI_off->GetYaxis()->SetLabelSize(0.035);
   hGenJet_Pythia84C_MPI_off->GetYaxis()->SetTitleSize(0.035);
   hGenJet_Pythia84C_MPI_off->GetYaxis()->SetTitleFont(42);
   hGenJet_Pythia84C_MPI_off->GetZaxis()->SetLabelFont(42);
   hGenJet_Pythia84C_MPI_off->GetZaxis()->SetLabelSize(0.035);
   hGenJet_Pythia84C_MPI_off->GetZaxis()->SetTitleSize(0.035);
   hGenJet_Pythia84C_MPI_off->GetZaxis()->SetTitleFont(42);
   hGenJet_Pythia84C_MPI_off->Draw("plsame");
   Double_t xAxis13[13] = {150, 187.5, 234.375, 292.9688, 366.2109, 457.7637, 572.2046, 715.2557, 894.0697, 1117.587, 1396.984, 1746.23, 2182.787}; 
   
   TH1D *hGenJet_Pythia84C_MPI_on = new TH1D("hGenJet_Pythia84C_MPI_on","Genjet energy Pythia84C, MPI on",12, xAxis13);
   hGenJet_Pythia84C_MPI_on->SetBinContent(0,0.5051282);
   hGenJet_Pythia84C_MPI_on->SetBinContent(1,0.5051282);
   hGenJet_Pythia84C_MPI_on->SetBinContent(2,0.30601);
   hGenJet_Pythia84C_MPI_on->SetBinContent(3,0.1700922);
   hGenJet_Pythia84C_MPI_on->SetBinContent(4,0.08738557);
   hGenJet_Pythia84C_MPI_on->SetBinContent(5,0.04234163);
   hGenJet_Pythia84C_MPI_on->SetBinContent(6,0.01961488);
   hGenJet_Pythia84C_MPI_on->SetBinContent(7,0.008597379);
   hGenJet_Pythia84C_MPI_on->SetBinContent(8,0.003495452);
   hGenJet_Pythia84C_MPI_on->SetBinContent(9,0.001343429);
   hGenJet_Pythia84C_MPI_on->SetBinContent(10,0.0004475201);
   hGenJet_Pythia84C_MPI_on->SetBinContent(11,0.0001128511);
   hGenJet_Pythia84C_MPI_on->SetBinContent(12,1.877581e-05);
   hGenJet_Pythia84C_MPI_on->SetBinContent(13,0.00089964);
   hGenJet_Pythia84C_MPI_on->SetBinError(0,0.06452997);
   hGenJet_Pythia84C_MPI_on->SetBinError(1,0.0004385804);
   hGenJet_Pythia84C_MPI_on->SetBinError(2,0.0003053242);
   hGenJet_Pythia84C_MPI_on->SetBinError(3,0.0002036014);
   hGenJet_Pythia84C_MPI_on->SetBinError(4,0.0001305279);
   hGenJet_Pythia84C_MPI_on->SetBinError(5,8.126662e-05);
   hGenJet_Pythia84C_MPI_on->SetBinError(6,4.947277e-05);
   hGenJet_Pythia84C_MPI_on->SetBinError(7,2.929554e-05);
   hGenJet_Pythia84C_MPI_on->SetBinError(8,1.670764e-05);
   hGenJet_Pythia84C_MPI_on->SetBinError(9,9.264368e-06);
   hGenJet_Pythia84C_MPI_on->SetBinError(10,4.782552e-06);
   hGenJet_Pythia84C_MPI_on->SetBinError(11,2.148084e-06);
   hGenJet_Pythia84C_MPI_on->SetBinError(12,7.836872e-07);
   hGenJet_Pythia84C_MPI_on->SetBinError(13,0.000113344);
   hGenJet_Pythia84C_MPI_on->SetEntries(2.448937e+07);

   ci = TColor::GetColor("#ff9999");
   hGenJet_Pythia84C_MPI_on->SetLineColor(ci);
   hGenJet_Pythia84C_MPI_on->SetLineStyle(5);
   hGenJet_Pythia84C_MPI_on->SetLineWidth(3);
   hGenJet_Pythia84C_MPI_on->SetMarkerStyle(26);
   hGenJet_Pythia84C_MPI_on->GetXaxis()->SetTitle("E_{det} (GeV)");
   hGenJet_Pythia84C_MPI_on->GetXaxis()->SetLabelFont(42);
   hGenJet_Pythia84C_MPI_on->GetXaxis()->SetLabelSize(0.035);
   hGenJet_Pythia84C_MPI_on->GetXaxis()->SetTitleSize(0.035);
   hGenJet_Pythia84C_MPI_on->GetXaxis()->SetTitleFont(42);
   hGenJet_Pythia84C_MPI_on->GetYaxis()->SetTitle("#frac{dN}{dE}");
   hGenJet_Pythia84C_MPI_on->GetYaxis()->SetLabelFont(42);
   hGenJet_Pythia84C_MPI_on->GetYaxis()->SetLabelSize(0.035);
   hGenJet_Pythia84C_MPI_on->GetYaxis()->SetTitleSize(0.035);
   hGenJet_Pythia84C_MPI_on->GetYaxis()->SetTitleFont(42);
   hGenJet_Pythia84C_MPI_on->GetZaxis()->SetLabelFont(42);
   hGenJet_Pythia84C_MPI_on->GetZaxis()->SetLabelSize(0.035);
   hGenJet_Pythia84C_MPI_on->GetZaxis()->SetTitleSize(0.035);
   hGenJet_Pythia84C_MPI_on->GetZaxis()->SetTitleFont(42);
   hGenJet_Pythia84C_MPI_on->Draw("plsame");
   Double_t xAxis14[13] = {150, 187.5, 234.375, 292.9688, 366.2109, 457.7637, 572.2046, 715.2557, 894.0697, 1117.587, 1396.984, 1746.23, 2182.787}; 
   
   TH1D *hGenJet_energy = new TH1D("hGenJet_energy","hGenJet_energy",12, xAxis14);
   hGenJet_energy->SetBinContent(0,0.6685529);
   hGenJet_energy->SetBinContent(1,0.6685529);
   hGenJet_energy->SetBinContent(2,0.4307103);
   hGenJet_energy->SetBinContent(3,0.2496474);
   hGenJet_energy->SetBinContent(4,0.1322579);
   hGenJet_energy->SetBinContent(5,0.06357385);
   hGenJet_energy->SetBinContent(6,0.02774823);
   hGenJet_energy->SetBinContent(7,0.01089358);
   hGenJet_energy->SetBinContent(8,0.003973309);
   hGenJet_energy->SetBinContent(9,0.001309663);
   hGenJet_energy->SetBinContent(10,0.0003556107);
   hGenJet_energy->SetBinContent(11,7.619856e-05);
   hGenJet_energy->SetBinContent(12,1.021563e-05);
   hGenJet_energy->SetBinContent(13,0.00065799);
   hGenJet_energy->SetBinError(0,0.1479347);
   hGenJet_energy->SetBinError(1,0.00114167);
   hGenJet_energy->SetBinError(2,0.0008196159);
   hGenJet_energy->SetBinError(3,0.0005581183);
   hGenJet_energy->SetBinError(4,0.0003633444);
   hGenJet_energy->SetBinError(5,0.000225316);
   hGenJet_energy->SetBinError(6,0.0001331421);
   hGenJet_energy->SetBinError(7,7.461534e-05);
   hGenJet_energy->SetBinError(8,4.030546e-05);
   hGenJet_energy->SetBinError(9,2.069724e-05);
   hGenJet_energy->SetBinError(10,9.646397e-06);
   hGenJet_energy->SetBinError(11,3.99389e-06);
   hGenJet_energy->SetBinError(12,1.307978e-06);
   hGenJet_energy->SetBinError(13,0.00021933);
   hGenJet_energy->SetEntries(5205900);

   ci = TColor::GetColor("#cccc00");
   hGenJet_energy->SetLineColor(ci);
   hGenJet_energy->SetLineStyle(6);
   hGenJet_energy->SetLineWidth(3);
   hGenJet_energy->SetMarkerStyle(27);
   hGenJet_energy->GetXaxis()->SetTitle("E_{det} (GeV)");
   hGenJet_energy->GetXaxis()->SetLabelFont(42);
   hGenJet_energy->GetXaxis()->SetLabelSize(0.035);
   hGenJet_energy->GetXaxis()->SetTitleSize(0.035);
   hGenJet_energy->GetXaxis()->SetTitleFont(42);
   hGenJet_energy->GetYaxis()->SetTitle("#frac{dN}{dE}");
   hGenJet_energy->GetYaxis()->SetLabelFont(42);
   hGenJet_energy->GetYaxis()->SetLabelSize(0.035);
   hGenJet_energy->GetYaxis()->SetTitleSize(0.035);
   hGenJet_energy->GetYaxis()->SetTitleFont(42);
   hGenJet_energy->GetZaxis()->SetLabelFont(42);
   hGenJet_energy->GetZaxis()->SetLabelSize(0.035);
   hGenJet_energy->GetZaxis()->SetTitleSize(0.035);
   hGenJet_energy->GetZaxis()->SetTitleFont(42);
   hGenJet_energy->Draw("plsame");
   Double_t xAxis15[13] = {150, 187.5, 234.375, 292.9688, 366.2109, 457.7637, 572.2046, 715.2557, 894.0697, 1117.587, 1396.984, 1746.23, 2182.787}; 
   
   TH1D *hGenJet_energy = new TH1D("hGenJet_energy","hGenJet_energy",12, xAxis15);
   hGenJet_energy->SetBinContent(0,0.6628924);
   hGenJet_energy->SetBinContent(1,0.6628924);
   hGenJet_energy->SetBinContent(2,0.4416972);
   hGenJet_energy->SetBinContent(3,0.2723311);
   hGenJet_energy->SetBinContent(4,0.1503449);
   hGenJet_energy->SetBinContent(5,0.07268219);
   hGenJet_energy->SetBinContent(6,0.03097814);
   hGenJet_energy->SetBinContent(7,0.01107);
   hGenJet_energy->SetBinContent(8,0.00334412);
   hGenJet_energy->SetBinContent(9,0.0008430193);
   hGenJet_energy->SetBinContent(10,0.0002174646);
   hGenJet_energy->SetBinContent(11,4.184129e-05);
   hGenJet_energy->SetBinContent(12,3.523477e-06);
   hGenJet_energy->SetBinContent(13,0.00038455);
   hGenJet_energy->SetBinError(0,0.3378028);
   hGenJet_energy->SetBinError(1,0.002607248);
   hGenJet_energy->SetBinError(2,0.001903567);
   hGenJet_energy->SetBinError(3,0.001336901);
   hGenJet_energy->SetBinError(4,0.0008884645);
   hGenJet_energy->SetBinError(5,0.0005525286);
   hGenJet_energy->SetBinError(6,0.0003226364);
   hGenJet_energy->SetBinError(7,0.0001725062);
   hGenJet_energy->SetBinError(8,8.480406e-05);
   hGenJet_energy->SetBinError(9,3.808373e-05);
   hGenJet_energy->SetBinError(10,1.730055e-05);
   hGenJet_energy->SetBinError(11,6.787553e-06);
   hGenJet_energy->SetBinError(12,1.761738e-06);
   hGenJet_energy->SetBinError(13,0.00038455);
   hGenJet_energy->SetEntries(993165);

   ci = TColor::GetColor("#33ff00");
   hGenJet_energy->SetLineColor(ci);
   hGenJet_energy->SetLineStyle(7);
   hGenJet_energy->SetLineWidth(3);
   hGenJet_energy->SetMarkerStyle(28);
   hGenJet_energy->GetXaxis()->SetTitle("E_{det} (GeV)");
   hGenJet_energy->GetXaxis()->SetLabelFont(42);
   hGenJet_energy->GetXaxis()->SetLabelSize(0.035);
   hGenJet_energy->GetXaxis()->SetTitleSize(0.035);
   hGenJet_energy->GetXaxis()->SetTitleFont(42);
   hGenJet_energy->GetYaxis()->SetTitle("#frac{dN}{dE}");
   hGenJet_energy->GetYaxis()->SetLabelFont(42);
   hGenJet_energy->GetYaxis()->SetLabelSize(0.035);
   hGenJet_energy->GetYaxis()->SetTitleSize(0.035);
   hGenJet_energy->GetYaxis()->SetTitleFont(42);
   hGenJet_energy->GetZaxis()->SetLabelFont(42);
   hGenJet_energy->GetZaxis()->SetLabelSize(0.035);
   hGenJet_energy->GetZaxis()->SetTitleSize(0.035);
   hGenJet_energy->GetZaxis()->SetTitleFont(42);
   hGenJet_energy->Draw("plsame");
   
   TLegend *leg = new TLegend(1210,3.562213e-321,1.130634e-312,2.697616e-316,NULL,"brNDC");
   leg->SetTextFont(62);
   leg->SetLineColor(1);
   leg->SetLineStyle(1);
   leg->SetLineWidth(2);
   leg->SetFillColor(0);
   leg->SetFillStyle(1001);
   TLegendEntry *entry=leg->AddEntry("response","Unfolded data, 30 it.","lp");
   entry->SetLineColor(1);
   entry->SetLineStyle(1);
   entry->SetLineWidth(1);
   entry->SetMarkerColor(1);
   entry->SetMarkerStyle(21);
   entry->SetMarkerSize(1);
   entry=leg->AddEntry("Graph","Syst. errors","f");
   entry->SetLineColor(1);
   entry->SetLineStyle(1);
   entry->SetLineWidth(1);
   entry->SetMarkerColor(1);
   entry->SetMarkerStyle(21);
   entry->SetMarkerSize(1);
   entry=leg->AddEntry("Graph","Syst. errors (no JES)","f");
   entry->SetLineColor(1);
   entry->SetLineStyle(1);
   entry->SetLineWidth(1);
   entry->SetMarkerColor(1);
   entry->SetMarkerStyle(21);
   entry->SetMarkerSize(1);
   entry=leg->AddEntry("hGenJet_energy_noCuts","Pythia6 (Z2*)","l");
   entry->SetLineColor(1);
   entry->SetLineStyle(1);
   entry->SetLineWidth(1);
   entry->SetMarkerColor(1);
   entry->SetMarkerStyle(21);
   entry->SetMarkerSize(1);
   entry=leg->AddEntry("hGenJet_energy_noCuts","Pythia8 (4C)","l");
   entry->SetLineColor(1);
   entry->SetLineStyle(1);
   entry->SetLineWidth(1);
   entry->SetMarkerColor(1);
   entry->SetMarkerStyle(21);
   entry->SetMarkerSize(1);
   entry=leg->AddEntry("hGenJet_energy_noCuts","EPOS","l");
   entry->SetLineColor(1);
   entry->SetLineStyle(1);
   entry->SetLineWidth(1);
   entry->SetMarkerColor(1);
   entry->SetMarkerStyle(21);
   entry->SetMarkerSize(1);
   entry=leg->AddEntry("hGenJet_Pythia84C_MPI_off","Pythia8.212 4C (no MPI)","lp");
   entry->SetLineColor(1);
   entry->SetLineStyle(1);
   entry->SetLineWidth(1);
   entry->SetMarkerColor(1);
   entry->SetMarkerStyle(21);
   entry->SetMarkerSize(1);
   entry=leg->AddEntry("hGenJet_Pythia84C_MPI_on","Pythia8.212 4C (MPI)","lp");
   entry->SetLineColor(1);
   entry->SetLineStyle(1);
   entry->SetLineWidth(1);
   entry->SetMarkerColor(1);
   entry->SetMarkerStyle(21);
   entry->SetMarkerSize(1);
   entry=leg->AddEntry("hGenJet_energy","QGSJet","lp");
   entry->SetLineColor(1);
   entry->SetLineStyle(1);
   entry->SetLineWidth(1);
   entry->SetMarkerColor(1);
   entry->SetMarkerStyle(21);
   entry->SetMarkerSize(1);
   entry=leg->AddEntry("hGenJet_energy","Sibyll","lp");
   entry->SetLineColor(1);
   entry->SetLineStyle(1);
   entry->SetLineWidth(1);
   entry->SetMarkerColor(1);
   entry->SetMarkerStyle(21);
   entry->SetMarkerSize(1);
   leg->Draw();
   Absolute_Values->Modified();
   Systematics->cd();
  
// ------------>Primitives in pad: Ratios
   Ratios = new TPad("Ratios", "Ratios",0,0,1,0.4);
   Ratios->Draw();
   Ratios->cd();
   Ratios->Range(0,0,1,1);
   Ratios->SetFillColor(0);
   Ratios->SetBorderMode(0);
   Ratios->SetBorderSize(2);
   Ratios->SetTickx(1);
   Ratios->SetTicky(1);
   Ratios->SetLeftMargin(0.2);
   Ratios->SetRightMargin(0.05);
   Ratios->SetTopMargin(0);
   Ratios->SetBottomMargin(0.15);
   Ratios->SetFrameFillStyle(0);
   Ratios->SetFrameBorderMode(0);
   
   grae = new TGraphAsymmErrors(24);
   grae->SetName("Graph");
   grae->SetTitle("TGraphAsymmErrors All");

   ci = TColor::GetColor("#c0c0c0");
   grae->SetFillColor(ci);
   grae->SetLineWidth(2);

   ci = TColor::GetColor("#c0c0c0");
   grae->SetMarkerColor(ci);
   grae->SetMarkerStyle(21);
   grae->SetPoint(0,150,1);
   grae->SetPointError(0,18.75,18.75,0.08579413,0.9654553);
   grae->SetPoint(1,187.5,1);
   grae->SetPointError(1,18.75,18.75,0.08579413,0.9654553);
   grae->SetPoint(2,187.5,1);
   grae->SetPointError(2,23.4375,23.4375,0.2025226,0.2868149);
   grae->SetPoint(3,234.375,1);
   grae->SetPointError(3,23.4375,23.4375,0.2025226,0.2868149);
   grae->SetPoint(4,234.375,1);
   grae->SetPointError(4,29.29688,29.29688,0.3200415,1.260279);
   grae->SetPoint(5,292.9688,1);
   grae->SetPointError(5,29.29688,29.29688,0.3200415,1.260279);
   grae->SetPoint(6,292.9688,1);
   grae->SetPointError(6,36.62109,36.62109,0.4256518,0.5085741);
   grae->SetPoint(7,366.2109,1);
   grae->SetPointError(7,36.62109,36.62109,0.4256518,0.5085741);
   grae->SetPoint(8,366.2109,1);
   grae->SetPointError(8,45.77637,45.77637,0.5017481,0.445194);
   grae->SetPoint(9,457.7637,1);
   grae->SetPointError(9,45.77637,45.77637,0.5017481,0.445194);
   grae->SetPoint(10,457.7637,1);
   grae->SetPointError(10,57.22046,57.22046,0.5681069,0.5753137);
   grae->SetPoint(11,572.2046,1);
   grae->SetPointError(11,57.22046,57.22046,0.5681069,0.5753137);
   grae->SetPoint(12,572.2046,1);
   grae->SetPointError(12,71.52557,71.52557,0.61306,0.9520055);
   grae->SetPoint(13,715.2557,1);
   grae->SetPointError(13,71.52557,71.52557,0.61306,0.9520055);
   grae->SetPoint(14,715.2557,1);
   grae->SetPointError(14,89.40697,89.40697,0.6907372,1.324243);
   grae->SetPoint(15,894.0697,1);
   grae->SetPointError(15,89.40697,89.40697,0.6907372,1.324243);
   grae->SetPoint(16,894.0697,1);
   grae->SetPointError(16,111.7587,111.7587,0.7478407,1.265491);
   grae->SetPoint(17,1117.587,1);
   grae->SetPointError(17,111.7587,111.7587,0.7478407,1.265491);
   grae->SetPoint(18,1117.587,1);
   grae->SetPointError(18,139.6984,139.6984,0.7769182,1.474068);
   grae->SetPoint(19,1396.984,1);
   grae->SetPointError(19,139.6984,139.6984,0.7769182,1.474068);
   grae->SetPoint(20,1396.984,1);
   grae->SetPointError(20,174.623,174.623,0.7995562,1.558776);
   grae->SetPoint(21,1746.23,1);
   grae->SetPointError(21,174.623,174.623,0.7995562,1.558776);
   grae->SetPoint(22,1746.23,1);
   grae->SetPointError(22,218.2787,218.2787,0.9763457,3.838775);
   grae->SetPoint(23,2182.787,1);
   grae->SetPointError(23,218.2787,218.2787,0.9763457,3.838775);
   
   TH1F *Graph3 = new TH1F("Graph3","TGraphAsymmErrors All",100,0,2628.048);
   Graph3->SetMinimum(0);
   Graph3->SetMaximum(3);
   Graph3->SetDirectory(0);
   Graph3->SetStats(0);
   Graph3->SetLineStyle(0);
   Graph3->SetMarkerStyle(20);
   Graph3->GetXaxis()->SetTitle("E (GeV)");
   Graph3->GetXaxis()->SetRange(13,67);
   Graph3->GetXaxis()->SetLabelFont(42);
   Graph3->GetXaxis()->SetLabelOffset(0.007);
   Graph3->GetXaxis()->SetLabelSize(0.05);
   Graph3->GetXaxis()->SetTitleSize(0.06);
   Graph3->GetXaxis()->SetTitleOffset(1.1);
   Graph3->GetXaxis()->SetTitleFont(42);
   Graph3->GetYaxis()->SetTitle("Ratio");
   Graph3->GetYaxis()->SetLabelFont(42);
   Graph3->GetYaxis()->SetLabelOffset(0.007);
   Graph3->GetYaxis()->SetLabelSize(0.05);
   Graph3->GetYaxis()->SetTitleSize(0.06);
   Graph3->GetYaxis()->SetTitleOffset(1.25);
   Graph3->GetYaxis()->SetTitleFont(42);
   Graph3->GetZaxis()->SetLabelFont(42);
   Graph3->GetZaxis()->SetLabelOffset(0.007);
   Graph3->GetZaxis()->SetLabelSize(0.05);
   Graph3->GetZaxis()->SetTitleSize(0.06);
   Graph3->GetZaxis()->SetTitleFont(42);
   grae->SetHistogram(Graph3);
   
   grae->Draw("ae3");
   
   grae = new TGraphAsymmErrors(24);
   grae->SetName("Graph");
   grae->SetTitle("TGraphAsymmErrors Example");
   grae->SetFillColor(13);
   grae->SetLineWidth(2);
   grae->SetMarkerColor(13);
   grae->SetMarkerStyle(21);
   grae->SetPoint(0,150,1);
   grae->SetPointError(0,18.75,18.75,0.02229006,0.01966653);
   grae->SetPoint(1,187.5,1);
   grae->SetPointError(1,18.75,18.75,0.02229006,0.01966653);
   grae->SetPoint(2,187.5,1);
   grae->SetPointError(2,23.4375,23.4375,0.1209137,0.09501827);
   grae->SetPoint(3,234.375,1);
   grae->SetPointError(3,23.4375,23.4375,0.1209137,0.09501827);
   grae->SetPoint(4,234.375,1);
   grae->SetPointError(4,29.29688,29.29688,0.2096155,0.2194416);
   grae->SetPoint(5,292.9688,1);
   grae->SetPointError(5,29.29688,29.29688,0.2096155,0.2194416);
   grae->SetPoint(6,292.9688,1);
   grae->SetPointError(6,36.62109,36.62109,0.2721206,0.3200042);
   grae->SetPoint(7,366.2109,1);
   grae->SetPointError(7,36.62109,36.62109,0.2721206,0.3200042);
   grae->SetPoint(8,366.2109,1);
   grae->SetPointError(8,45.77637,45.77637,0.2932961,0.4054358);
   grae->SetPoint(9,457.7637,1);
   grae->SetPointError(9,45.77637,45.77637,0.2932961,0.4054358);
   grae->SetPoint(10,457.7637,1);
   grae->SetPointError(10,57.22046,57.22046,0.3037491,0.4535598);
   grae->SetPoint(11,572.2046,1);
   grae->SetPointError(11,57.22046,57.22046,0.3037491,0.4535598);
   grae->SetPoint(12,572.2046,1);
   grae->SetPointError(12,71.52557,71.52557,0.2864019,0.4872001);
   grae->SetPoint(13,715.2557,1);
   grae->SetPointError(13,71.52557,71.52557,0.2864019,0.4872001);
   grae->SetPoint(14,715.2557,1);
   grae->SetPointError(14,89.40697,89.40697,0.3585559,0.5352182);
   grae->SetPoint(15,894.0697,1);
   grae->SetPointError(15,89.40697,89.40697,0.3585559,0.5352182);
   grae->SetPoint(16,894.0697,1);
   grae->SetPointError(16,111.7587,111.7587,0.3705744,0.470755);
   grae->SetPoint(17,1117.587,1);
   grae->SetPointError(17,111.7587,111.7587,0.3705744,0.470755);
   grae->SetPoint(18,1117.587,1);
   grae->SetPointError(18,139.6984,139.6984,0.3595023,0.3354415);
   grae->SetPoint(19,1396.984,1);
   grae->SetPointError(19,139.6984,139.6984,0.3595023,0.3354415);
   grae->SetPoint(20,1396.984,1);
   grae->SetPointError(20,174.623,174.623,0.4012661,0.7245039);
   grae->SetPoint(21,1746.23,1);
   grae->SetPointError(21,174.623,174.623,0.4012661,0.7245039);
   grae->SetPoint(22,1746.23,1);
   grae->SetPointError(22,218.2787,218.2787,0.7255003,1.514979);
   grae->SetPoint(23,2182.787,1);
   grae->SetPointError(23,218.2787,218.2787,0.7255003,1.514979);
   
   TH1F *Graph4 = new TH1F("Graph4","TGraphAsymmErrors Example",100,0,2628.048);
   Graph4->SetMinimum(0);
   Graph4->SetMaximum(2);
   Graph4->SetDirectory(0);
   Graph4->SetStats(0);
   Graph4->SetLineStyle(0);
   Graph4->SetMarkerStyle(20);
   Graph4->GetXaxis()->SetTitle("E (GeV)");
   Graph4->GetXaxis()->SetRange(12,67);
   Graph4->GetXaxis()->SetLabelFont(42);
   Graph4->GetXaxis()->SetLabelOffset(0.007);
   Graph4->GetXaxis()->SetLabelSize(0.05);
   Graph4->GetXaxis()->SetTitleSize(0.06);
   Graph4->GetXaxis()->SetTitleOffset(1.1);
   Graph4->GetXaxis()->SetTitleFont(42);
   Graph4->GetYaxis()->SetLabelFont(42);
   Graph4->GetYaxis()->SetLabelOffset(0.007);
   Graph4->GetYaxis()->SetLabelSize(0.05);
   Graph4->GetYaxis()->SetTitleSize(0.06);
   Graph4->GetYaxis()->SetTitleOffset(1.25);
   Graph4->GetYaxis()->SetTitleFont(42);
   Graph4->GetZaxis()->SetLabelFont(42);
   Graph4->GetZaxis()->SetLabelOffset(0.007);
   Graph4->GetZaxis()->SetLabelSize(0.05);
   Graph4->GetZaxis()->SetTitleSize(0.06);
   Graph4->GetZaxis()->SetTitleFont(42);
   grae->SetHistogram(Graph4);
   
   grae->Draw("pe3");
   
   grae = new TGraphAsymmErrors(12);
   grae->SetName("Graph");
   grae->SetTitle("Graph");
   grae->SetFillColor(1);
   grae->SetLineWidth(2);
   grae->SetMarkerStyle(20);
   grae->SetMarkerSize(1.25);
   grae->SetPoint(0,168.3176,1);
   grae->SetPointError(0,18.31763,19.18237,0.001506591,0.001506591);
   grae->SetPoint(1,210.2622,1);
   grae->SetPointError(1,22.76216,24.11284,0.00182961,0.00182961);
   grae->SetPoint(2,262.6172,1);
   grae->SetPointError(2,28.24225,30.3515,0.001729354,0.001729354);
   grae->SetPoint(3,327.9434,1);
   grae->SetPointError(3,34.97468,38.26751,0.002078752,0.002078752);
   grae->SetPoint(4,409.4183,1);
   grae->SetPointError(4,43.20734,48.3454,0.002617346,0.002617346);
   grae->SetPoint(5,510.9785,1);
   grae->SetPointError(5,53.21482,61.2261,0.003489862,0.003489862);
   grae->SetPoint(6,637.4918,1);
   grae->SetPointError(6,65.28724,77.7639,0.005436368,0.005436368);
   grae->SetPoint(7,794.9645,1);
   grae->SetPointError(7,79.70878,99.10515,0.008792714,0.008792714);
   grae->SetPoint(8,990.7922,1);
   grae->SetPointError(8,96.72253,126.7949,0.01543097,0.01543097);
   grae->SetPoint(9,1242.253,1);
   grae->SetPointError(9,124.666,154.7308,0.02092925,0.02092925);
   grae->SetPoint(10,1548.298,1);
   grae->SetPointError(10,151.3137,197.9323,0.0528388,0.0528388);
   grae->SetPoint(11,1928.509,1);
   grae->SetPointError(11,182.279,254.2784,0.05205661,0.05205661);
   grae->Draw("ep");
   Double_t xAxis16[13] = {150, 187.5, 234.375, 292.9688, 366.2109, 457.7637, 572.2046, 715.2557, 894.0697, 1117.587, 1396.984, 1746.23, 2182.787}; 
   
   TH1D *hGen_ratio = new TH1D("hGen_ratio","GenJet energy distribution",12, xAxis16);
   hGen_ratio->SetBinContent(1,0.9427556);
   hGen_ratio->SetBinContent(2,0.8422087);
   hGen_ratio->SetBinContent(3,0.7780449);
   hGen_ratio->SetBinContent(4,0.7482916);
   hGen_ratio->SetBinContent(5,0.7708231);
   hGen_ratio->SetBinContent(6,0.8608539);
   hGen_ratio->SetBinContent(7,1.061175);
   hGen_ratio->SetBinContent(8,1.376232);
   hGen_ratio->SetBinContent(9,1.579301);
   hGen_ratio->SetBinContent(10,1.577072);
   hGen_ratio->SetBinContent(11,1.68426);
   hGen_ratio->SetBinContent(12,1.617741);
   hGen_ratio->SetBinError(1,0.2272693);
   hGen_ratio->SetBinError(2,0.2310688);
   hGen_ratio->SetBinError(3,0.2541349);
   hGen_ratio->SetBinError(4,0.3017505);
   hGen_ratio->SetBinError(5,0.4003382);
   hGen_ratio->SetBinError(6,0.5852269);
   hGen_ratio->SetBinError(7,0.9808809);
   hGen_ratio->SetBinError(8,1.808583);
   hGen_ratio->SetBinError(9,3.010996);
   hGen_ratio->SetBinError(10,4.688591);
   hGen_ratio->SetBinError(11,8.894373);
   hGen_ratio->SetBinError(12,20.19575);
   hGen_ratio->SetEntries(0.3661618);
   hGen_ratio->SetStats(0);
   hGen_ratio->SetLineColor(2);
   hGen_ratio->SetLineWidth(3);
   hGen_ratio->SetMarkerStyle(20);
   hGen_ratio->GetXaxis()->SetTitle("E_{det} (GeV)");
   hGen_ratio->GetXaxis()->SetLabelFont(42);
   hGen_ratio->GetXaxis()->SetLabelOffset(0.007);
   hGen_ratio->GetXaxis()->SetLabelSize(0.05);
   hGen_ratio->GetXaxis()->SetTitleSize(0.06);
   hGen_ratio->GetXaxis()->SetTitleOffset(0.9);
   hGen_ratio->GetXaxis()->SetTitleFont(42);
   hGen_ratio->GetYaxis()->SetTitle("#frac{dN}{dE}");
   hGen_ratio->GetYaxis()->SetLabelFont(42);
   hGen_ratio->GetYaxis()->SetLabelOffset(0.007);
   hGen_ratio->GetYaxis()->SetLabelSize(0.05);
   hGen_ratio->GetYaxis()->SetTitleSize(0.06);
   hGen_ratio->GetYaxis()->SetTitleOffset(1.05);
   hGen_ratio->GetYaxis()->SetTitleFont(42);
   hGen_ratio->GetZaxis()->SetLabelFont(42);
   hGen_ratio->GetZaxis()->SetLabelOffset(0.007);
   hGen_ratio->GetZaxis()->SetLabelSize(0.05);
   hGen_ratio->GetZaxis()->SetTitleSize(0.06);
   hGen_ratio->GetZaxis()->SetTitleFont(42);
   hGen_ratio->Draw("histsame");
   Double_t xAxis17[13] = {150, 187.5, 234.375, 292.9688, 366.2109, 457.7637, 572.2046, 715.2557, 894.0697, 1117.587, 1396.984, 1746.23, 2182.787}; 
   
   TH1D *hGen_ratio = new TH1D("hGen_ratio","GenJet energy distribution",12, xAxis17);
   hGen_ratio->SetBinContent(1,1.048524);
   hGen_ratio->SetBinContent(2,0.917985);
   hGen_ratio->SetBinContent(3,0.83769);
   hGen_ratio->SetBinContent(4,0.797788);
   hGen_ratio->SetBinContent(5,0.8261102);
   hGen_ratio->SetBinContent(6,0.9150641);
   hGen_ratio->SetBinContent(7,1.135795);
   hGen_ratio->SetBinContent(8,1.546167);
   hGen_ratio->SetBinContent(9,1.71231);
   hGen_ratio->SetBinContent(10,1.64933);
   hGen_ratio->SetBinContent(11,1.921003);
   hGen_ratio->SetBinContent(12,2.351916);
   hGen_ratio->SetBinError(1,0.2396797);
   hGen_ratio->SetBinError(2,0.2412404);
   hGen_ratio->SetBinError(3,0.2636963);
   hGen_ratio->SetBinError(4,0.3115707);
   hGen_ratio->SetBinError(5,0.4144471);
   hGen_ratio->SetBinError(6,0.6033728);
   hGen_ratio->SetBinError(7,1.014783);
   hGen_ratio->SetBinError(8,1.917);
   hGen_ratio->SetBinError(9,3.135235);
   hGen_ratio->SetBinError(10,4.794803);
   hGen_ratio->SetBinError(11,9.499);
   hGen_ratio->SetBinError(12,24.35107);
   hGen_ratio->SetEntries(0.3398603);
   hGen_ratio->SetStats(0);

   ci = TColor::GetColor("#ff9933");
   hGen_ratio->SetLineColor(ci);
   hGen_ratio->SetLineStyle(2);
   hGen_ratio->SetLineWidth(3);
   hGen_ratio->SetMarkerStyle(20);
   hGen_ratio->GetXaxis()->SetTitle("E_{det} (GeV)");
   hGen_ratio->GetXaxis()->SetLabelFont(42);
   hGen_ratio->GetXaxis()->SetLabelOffset(0.007);
   hGen_ratio->GetXaxis()->SetLabelSize(0.05);
   hGen_ratio->GetXaxis()->SetTitleSize(0.06);
   hGen_ratio->GetXaxis()->SetTitleOffset(0.9);
   hGen_ratio->GetXaxis()->SetTitleFont(42);
   hGen_ratio->GetYaxis()->SetTitle("#frac{dN}{dE}");
   hGen_ratio->GetYaxis()->SetLabelFont(42);
   hGen_ratio->GetYaxis()->SetLabelOffset(0.007);
   hGen_ratio->GetYaxis()->SetLabelSize(0.05);
   hGen_ratio->GetYaxis()->SetTitleSize(0.06);
   hGen_ratio->GetYaxis()->SetTitleOffset(1.05);
   hGen_ratio->GetYaxis()->SetTitleFont(42);
   hGen_ratio->GetZaxis()->SetLabelFont(42);
   hGen_ratio->GetZaxis()->SetLabelOffset(0.007);
   hGen_ratio->GetZaxis()->SetLabelSize(0.05);
   hGen_ratio->GetZaxis()->SetTitleSize(0.06);
   hGen_ratio->GetZaxis()->SetTitleFont(42);
   hGen_ratio->Draw("histsame");
   Double_t xAxis18[13] = {150, 187.5, 234.375, 292.9688, 366.2109, 457.7637, 572.2046, 715.2557, 894.0697, 1117.587, 1396.984, 1746.23, 2182.787}; 
   
   TH1D *hGen_ratio = new TH1D("hGen_ratio","GenJet energy distribution",12, xAxis18);
   hGen_ratio->SetBinContent(1,1.286271);
   hGen_ratio->SetBinContent(2,1.220999);
   hGen_ratio->SetBinContent(3,1.212676);
   hGen_ratio->SetBinContent(4,1.244381);
   hGen_ratio->SetBinContent(5,1.363275);
   hGen_ratio->SetBinContent(6,1.495058);
   hGen_ratio->SetBinContent(7,1.752157);
   hGen_ratio->SetBinContent(8,2.025447);
   hGen_ratio->SetBinContent(9,1.866997);
   hGen_ratio->SetBinContent(10,1.366013);
   hGen_ratio->SetBinContent(11,0.9898474);
   hGen_ratio->SetBinContent(12,0.6831845);
   hGen_ratio->SetBinError(1,0.265467);
   hGen_ratio->SetBinError(2,0.2782233);
   hGen_ratio->SetBinError(3,0.3172763);
   hGen_ratio->SetBinError(4,0.3891284);
   hGen_ratio->SetBinError(5,0.5324094);
   hGen_ratio->SetBinError(6,0.7712457);
   hGen_ratio->SetBinError(7,1.260416);
   hGen_ratio->SetBinError(8,2.194106);
   hGen_ratio->SetBinError(9,3.2738);
   hGen_ratio->SetBinError(10,4.363576);
   hGen_ratio->SetBinError(11,6.818455);
   hGen_ratio->SetBinError(12,13.12418);
   hGen_ratio->SetEntries(1.063562);
   hGen_ratio->SetStats(0);

   ci = TColor::GetColor("#0000ff");
   hGen_ratio->SetLineColor(ci);
   hGen_ratio->SetLineStyle(3);
   hGen_ratio->SetLineWidth(3);
   hGen_ratio->SetMarkerStyle(20);
   hGen_ratio->GetXaxis()->SetTitle("E_{det} (GeV)");
   hGen_ratio->GetXaxis()->SetLabelFont(42);
   hGen_ratio->GetXaxis()->SetLabelOffset(0.007);
   hGen_ratio->GetXaxis()->SetLabelSize(0.05);
   hGen_ratio->GetXaxis()->SetTitleSize(0.06);
   hGen_ratio->GetXaxis()->SetTitleOffset(0.9);
   hGen_ratio->GetXaxis()->SetTitleFont(42);
   hGen_ratio->GetYaxis()->SetTitle("#frac{dN}{dE}");
   hGen_ratio->GetYaxis()->SetLabelFont(42);
   hGen_ratio->GetYaxis()->SetLabelOffset(0.007);
   hGen_ratio->GetYaxis()->SetLabelSize(0.05);
   hGen_ratio->GetYaxis()->SetTitleSize(0.06);
   hGen_ratio->GetYaxis()->SetTitleOffset(1.05);
   hGen_ratio->GetYaxis()->SetTitleFont(42);
   hGen_ratio->GetZaxis()->SetLabelFont(42);
   hGen_ratio->GetZaxis()->SetLabelOffset(0.007);
   hGen_ratio->GetZaxis()->SetLabelSize(0.05);
   hGen_ratio->GetZaxis()->SetTitleSize(0.06);
   hGen_ratio->GetZaxis()->SetTitleFont(42);
   hGen_ratio->Draw("histsame");
   Double_t xAxis19[13] = {150, 187.5, 234.375, 292.9688, 366.2109, 457.7637, 572.2046, 715.2557, 894.0697, 1117.587, 1396.984, 1746.23, 2182.787}; 
   
   TH1D *hGen_ratio = new TH1D("hGen_ratio","Genjet energy Pythia84C, MPI off",12, xAxis19);
   hGen_ratio->SetBinContent(1,0.7368554);
   hGen_ratio->SetBinContent(2,0.5937534);
   hGen_ratio->SetBinContent(3,0.4834375);
   hGen_ratio->SetBinContent(4,0.4037412);
   hGen_ratio->SetBinContent(5,0.3611396);
   hGen_ratio->SetBinContent(6,0.3551);
   hGen_ratio->SetBinContent(7,0.4269169);
   hGen_ratio->SetBinContent(8,0.5808591);
   hGen_ratio->SetBinContent(9,0.7467474);
   hGen_ratio->SetBinContent(10,0.8130774);
   hGen_ratio->SetBinContent(11,1.063006);
   hGen_ratio->SetBinContent(12,1.368092);
   hGen_ratio->SetBinError(1,0.002028486);
   hGen_ratio->SetBinError(2,0.001966628);
   hGen_ratio->SetBinError(3,0.001887889);
   hGen_ratio->SetBinError(4,0.002052322);
   hGen_ratio->SetBinError(5,0.00250093);
   hGen_ratio->SetBinError(6,0.003409196);
   hGen_ratio->SetBinError(7,0.0057465);
   hGen_ratio->SetBinError(8,0.01116479);
   hGen_ratio->SetBinError(9,0.02094839);
   hGen_ratio->SetBinError(10,0.03314748);
   hGen_ratio->SetBinError(11,0.08197223);
   hGen_ratio->SetBinError(12,0.172335);
   hGen_ratio->SetEntries(1649.592);

   ci = TColor::GetColor("#cc00cc");
   hGen_ratio->SetLineColor(ci);
   hGen_ratio->SetLineStyle(4);
   hGen_ratio->SetLineWidth(3);
   hGen_ratio->SetMarkerStyle(25);
   hGen_ratio->GetXaxis()->SetTitle("E_{det} (GeV)");
   hGen_ratio->GetXaxis()->SetLabelFont(42);
   hGen_ratio->GetXaxis()->SetLabelSize(0.035);
   hGen_ratio->GetXaxis()->SetTitleSize(0.035);
   hGen_ratio->GetXaxis()->SetTitleFont(42);
   hGen_ratio->GetYaxis()->SetTitle("#frac{dN}{dE}");
   hGen_ratio->GetYaxis()->SetLabelFont(42);
   hGen_ratio->GetYaxis()->SetLabelSize(0.035);
   hGen_ratio->GetYaxis()->SetTitleSize(0.035);
   hGen_ratio->GetYaxis()->SetTitleFont(42);
   hGen_ratio->GetZaxis()->SetLabelFont(42);
   hGen_ratio->GetZaxis()->SetLabelSize(0.035);
   hGen_ratio->GetZaxis()->SetTitleSize(0.035);
   hGen_ratio->GetZaxis()->SetTitleFont(42);
   hGen_ratio->Draw("plsame");
   Double_t xAxis20[13] = {150, 187.5, 234.375, 292.9688, 366.2109, 457.7637, 572.2046, 715.2557, 894.0697, 1117.587, 1396.984, 1746.23, 2182.787}; 
   
   TH1D *hGen_ratio = new TH1D("hGen_ratio","Genjet energy Pythia84C, MPI on",12, xAxis20);
   hGen_ratio->SetBinContent(1,1.037763);
   hGen_ratio->SetBinContent(2,0.9093278);
   hGen_ratio->SetBinContent(3,0.8272702);
   hGen_ratio->SetBinContent(4,0.7787795);
   hGen_ratio->SetBinContent(5,0.8059855);
   hGen_ratio->SetBinContent(6,0.8930477);
   hGen_ratio->SetBinContent(7,1.115033);
   hGen_ratio->SetBinContent(8,1.485492);
   hGen_ratio->SetBinContent(9,1.723667);
   hGen_ratio->SetBinContent(10,1.742792);
   hGen_ratio->SetBinContent(11,1.851038);
   hGen_ratio->SetBinContent(12,2.066539);
   hGen_ratio->SetBinError(1,0.001804539);
   hGen_ratio->SetBinError(2,0.001895026);
   hGen_ratio->SetBinError(3,0.001739922);
   hGen_ratio->SetBinError(4,0.001993485);
   hGen_ratio->SetBinError(5,0.002615946);
   hGen_ratio->SetBinError(6,0.003845362);
   hGen_ratio->SetBinError(7,0.007154059);
   hGen_ratio->SetBinError(8,0.01486669);
   hGen_ratio->SetBinError(9,0.02913305);
   hGen_ratio->SetBinError(10,0.04095527);
   hGen_ratio->SetBinError(11,0.1039594);
   hGen_ratio->SetBinError(12,0.1378871);
   hGen_ratio->SetEntries(7109.603);

   ci = TColor::GetColor("#ff9999");
   hGen_ratio->SetLineColor(ci);
   hGen_ratio->SetLineStyle(5);
   hGen_ratio->SetLineWidth(3);
   hGen_ratio->SetMarkerStyle(26);
   hGen_ratio->GetXaxis()->SetTitle("E_{det} (GeV)");
   hGen_ratio->GetXaxis()->SetLabelFont(42);
   hGen_ratio->GetXaxis()->SetLabelSize(0.035);
   hGen_ratio->GetXaxis()->SetTitleSize(0.035);
   hGen_ratio->GetXaxis()->SetTitleFont(42);
   hGen_ratio->GetYaxis()->SetTitle("#frac{dN}{dE}");
   hGen_ratio->GetYaxis()->SetLabelFont(42);
   hGen_ratio->GetYaxis()->SetLabelSize(0.035);
   hGen_ratio->GetYaxis()->SetTitleSize(0.035);
   hGen_ratio->GetYaxis()->SetTitleFont(42);
   hGen_ratio->GetZaxis()->SetLabelFont(42);
   hGen_ratio->GetZaxis()->SetLabelSize(0.035);
   hGen_ratio->GetZaxis()->SetTitleSize(0.035);
   hGen_ratio->GetZaxis()->SetTitleFont(42);
   hGen_ratio->Draw("plsame");
   Double_t xAxis21[13] = {150, 187.5, 234.375, 292.9688, 366.2109, 457.7637, 572.2046, 715.2557, 894.0697, 1117.587, 1396.984, 1746.23, 2182.787}; 
   
   TH1D *hGen_ratio = new TH1D("hGen_ratio","hGenJet_energy",12, xAxis21);
   hGen_ratio->SetBinContent(1,1.373512);
   hGen_ratio->SetBinContent(2,1.279883);
   hGen_ratio->SetBinContent(3,1.214199);
   hGen_ratio->SetBinContent(4,1.178681);
   hGen_ratio->SetBinContent(5,1.210147);
   hGen_ratio->SetBinContent(6,1.263352);
   hGen_ratio->SetBinContent(7,1.412838);
   hGen_ratio->SetBinContent(8,1.688571);
   hGen_ratio->SetBinContent(9,1.680343);
   hGen_ratio->SetBinContent(10,1.384866);
   hGen_ratio->SetBinContent(11,1.249845);
   hGen_ratio->SetBinContent(12,1.124372);
   hGen_ratio->SetBinError(1,0.003127859);
   hGen_ratio->SetBinError(2,0.003378661);
   hGen_ratio->SetBinError(3,0.003431846);
   hGen_ratio->SetBinError(4,0.004060646);
   hGen_ratio->SetBinError(5,0.005331736);
   hGen_ratio->SetBinError(6,0.007495632);
   hGen_ratio->SetBinError(7,0.01235482);
   hGen_ratio->SetBinError(8,0.022668);
   hGen_ratio->SetBinError(9,0.03711486);
   hGen_ratio->SetBinError(10,0.04744797);
   hGen_ratio->SetBinError(11,0.09302067);
   hGen_ratio->SetBinError(12,0.1554049);
   hGen_ratio->SetEntries(6927.831);

   ci = TColor::GetColor("#cccc00");
   hGen_ratio->SetLineColor(ci);
   hGen_ratio->SetLineStyle(6);
   hGen_ratio->SetLineWidth(3);
   hGen_ratio->SetMarkerStyle(27);
   hGen_ratio->GetXaxis()->SetTitle("E_{det} (GeV)");
   hGen_ratio->GetXaxis()->SetLabelFont(42);
   hGen_ratio->GetXaxis()->SetLabelSize(0.035);
   hGen_ratio->GetXaxis()->SetTitleSize(0.035);
   hGen_ratio->GetXaxis()->SetTitleFont(42);
   hGen_ratio->GetYaxis()->SetTitle("#frac{dN}{dE}");
   hGen_ratio->GetYaxis()->SetLabelFont(42);
   hGen_ratio->GetYaxis()->SetLabelSize(0.035);
   hGen_ratio->GetYaxis()->SetTitleSize(0.035);
   hGen_ratio->GetYaxis()->SetTitleFont(42);
   hGen_ratio->GetZaxis()->SetLabelFont(42);
   hGen_ratio->GetZaxis()->SetLabelSize(0.035);
   hGen_ratio->GetZaxis()->SetTitleSize(0.035);
   hGen_ratio->GetZaxis()->SetTitleFont(42);
   hGen_ratio->Draw("plsame");
   Double_t xAxis22[13] = {150, 187.5, 234.375, 292.9688, 366.2109, 457.7637, 572.2046, 715.2557, 894.0697, 1117.587, 1396.984, 1746.23, 2182.787}; 
   
   TH1D *hGen_ratio = new TH1D("hGen_ratio","hGenJet_energy",12, xAxis22);
   hGen_ratio->SetBinContent(1,1.361883);
   hGen_ratio->SetBinContent(2,1.312531);
   hGen_ratio->SetBinContent(3,1.324525);
   hGen_ratio->SetBinContent(4,1.339873);
   hGen_ratio->SetBinContent(5,1.383527);
   hGen_ratio->SetBinContent(6,1.410407);
   hGen_ratio->SetBinContent(7,1.435719);
   hGen_ratio->SetBinContent(8,1.421179);
   hGen_ratio->SetBinContent(9,1.081623);
   hGen_ratio->SetBinContent(10,0.8468794);
   hGen_ratio->SetBinContent(11,0.6863008);
   hGen_ratio->SetBinContent(12,0.3878075);
   hGen_ratio->SetBinError(1,0.005736);
   hGen_ratio->SetBinError(2,0.006145207);
   hGen_ratio->SetBinError(3,0.006893891);
   hGen_ratio->SetBinError(4,0.008393583);
   hGen_ratio->SetBinError(5,0.01112347);
   hGen_ratio->SetBinError(6,0.01549207);
   hGen_ratio->SetBinError(7,0.02369548);
   hGen_ratio->SetBinError(8,0.0381448);
   hGen_ratio->SetBinError(9,0.0516347);
   hGen_ratio->SetBinError(10,0.06966656);
   hGen_ratio->SetBinError(11,0.1170897);
   hGen_ratio->SetBinError(12,0.1949518);
   hGen_ratio->SetEntries(3167.774);

   ci = TColor::GetColor("#33ff00");
   hGen_ratio->SetLineColor(ci);
   hGen_ratio->SetLineStyle(7);
   hGen_ratio->SetLineWidth(3);
   hGen_ratio->SetMarkerStyle(28);
   hGen_ratio->GetXaxis()->SetTitle("E_{det} (GeV)");
   hGen_ratio->GetXaxis()->SetLabelFont(42);
   hGen_ratio->GetXaxis()->SetLabelSize(0.035);
   hGen_ratio->GetXaxis()->SetTitleSize(0.035);
   hGen_ratio->GetXaxis()->SetTitleFont(42);
   hGen_ratio->GetYaxis()->SetTitle("#frac{dN}{dE}");
   hGen_ratio->GetYaxis()->SetLabelFont(42);
   hGen_ratio->GetYaxis()->SetLabelSize(0.035);
   hGen_ratio->GetYaxis()->SetTitleSize(0.035);
   hGen_ratio->GetYaxis()->SetTitleFont(42);
   hGen_ratio->GetZaxis()->SetLabelFont(42);
   hGen_ratio->GetZaxis()->SetLabelSize(0.035);
   hGen_ratio->GetZaxis()->SetTitleSize(0.035);
   hGen_ratio->GetZaxis()->SetTitleFont(42);
   hGen_ratio->Draw("plsame");
   Ratios->Modified();
   Systematics->cd();
   Systematics->Modified();
   Systematics->cd();
   Systematics->SetSelected(Systematics);
}