void bhw_spectrum(const Char_t *fileName) {
  TFile *file = new TFile(fileName);

  TH1D *strue = (TH1D*)file->Get("total");
  TH1D *reco = (TH1D*)file->Get("xreco");

  reco->SetMarkerStyle(24);
  strue->SetLineColor(3);

  strue->DrawCopy("hist");
  reco->DrawCopy("samee1");
  gPad->SetLogy();
}

TH1D* bhw_mcfactor(const Char_t *fileName) {
  TFile *file = new TFile(fileName);

  TH1D *strue = (TH1D*)file->Get("total");
  TH1D *reco = (TH1D*)file->Get("xreco");

  reco->Divide(strue);
  reco->Draw("e1");
  gPad->SetLogy(kFALSE);

  return reco;
}