void makePlot(int iShow){
  gROOT->Reset();
  gROOT->SetStyle("Plain");
  gStyle->SetOptStat(0);

  TFile *hFile1 = new TFile("mcOut.root");

  TH1D *hDen; hFile1->GetObject("hEthrown",hDen);
  hDen->SetName("hDen");

  TH1D *hNum; hFile1->GetObject("hEout40",hNum);
  hNum->SetName("hNum");

  hNum->Sumw2();

  TH1D *hDiv = (TH1D*) hDen->Clone("hDiv");
  
  hDiv->Divide(hNum,hDen);
  //Fancy plot stuff
  hDiv->GetXaxis()->SetTitle("Energy (MeV)");
  //hDiv->GetYaxis()->SetTitle("Fraction of #gamma Me_{Lower}^{Upper}");
  hDiv->GetYaxis()->SetTitle("Fraction of #gamma");

  //See https://root.cern.ch/doc/master/classTAttMarker.html
  hDiv->SetMarkerStyle(20);
  hDiv->SetMarkerColor(2);
  hDiv->SetLineColor(2);
  hDiv->GetXaxis()->SetRangeUser(0.0,500.0);
  hDiv->SetMinimum(0.6);
  hDiv->SetMaximum(1.1);
  //hDiv->Fit("pol0");

  hDiv->Draw("e1");

}