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(); }