void make_plots(TFile *event_file, TTree *T, TFile *out_file)
{
  TH1F *e_ke = new TH1F("e_ke", "e- (or e+) Kinetic Energy",1000,0,3.0);
  e_ke->SetYTitle("Events per bin");
  e_ke->SetXTitle("Energy (MeV)");

  TH1F *gamma_ke = new TH1F("gamma_ke", "Gamma Kinetic Energy",1000,0,3.0);
  gamma_ke->SetYTitle("Events per bin");
  gamma_ke->SetXTitle("Gamma Energy (MeV)");

  TH1F *alpha_ke = new TH1F("alpha_ke", "Alpha Kinetic Energy",1000,0,10.0);
  alpha_ke->SetYTitle("Events per bin");
  alpha_ke->SetXTitle("Alpha Energy (MeV)");

  TH1F *time = new TH1F("time", "Initial Time of Particles",1000,0,10000.0);
  time->SetYTitle("Events per bin");
  time->SetXTitle("time (ns)");

  T->Draw("ds.mc.particles.kineticEnergy>>e_ke", "ds.mc.particles.pdgCode==11||ds.mc.particles.pdgCode==-11","goff");
  T->Draw("ds.mc.particles.kineticEnergy>>gamma_ke", "ds.mc.particles.pdgCode==22","goff");
  T->Draw("ds.mc.particles.kineticEnergy>>alpha_ke", "ds.mc.particles.pdgCode==1000020040","goff");
  T->Draw("ds.mc.particles.mcTime>>time","","goff");

  e_ke->Write();
  gamma_ke->Write();
  alpha_ke->Write();
  time->Write();
}