From 35e2ebfbf1313fa753d5b18787b2d551ec243a33 Mon Sep 17 00:00:00 2001 From: tlatorre Date: Mon, 6 Jul 2020 09:41:33 -0500 Subject: add first draft of plot-muons This commit adds a first draft of a script to plot the michel energy distribution and particle id histograms for data and Monte Carlo and to plot the energy bias and resolution for stopping muons. --- utils/plot-energy | 2 -- 1 file changed, 2 deletions(-) (limited to 'utils/plot-energy') diff --git a/utils/plot-energy b/utils/plot-energy index ef34f97..4568bcc 100755 --- a/utils/plot-energy +++ b/utils/plot-energy @@ -295,8 +295,6 @@ if __name__ == '__main__': stopping_muons['T_dx'] = dx_to_energy(stopping_muons.dx) stopping_muons['dT'] = stopping_muons['energy1'] - stopping_muons['T_dx'] - print(stopping_muons[['r','r_psup','dx','T_dx','energy1','r_michel','x','y','z','x_michel','y_michel','z_michel','ftp_r_michel']]) - fig = plt.figure() plt.hist((stopping_muons['energy1']-stopping_muons['T_dx'])*100/stopping_muons['T_dx'], bins=np.linspace(-100,100,200), histtype='step') despine(fig,trim=True) -- cgit