From e79e1fcf612c386f5f808bfa9099af723fc28a76 Mon Sep 17 00:00:00 2001 From: tlatorre Date: Mon, 6 Jul 2020 12:27:27 -0500 Subject: update plot-muons to plot energy distributions for each particle ID for stopping muons --- utils/plot-muons | 12 ++++++++++++ 1 file changed, 12 insertions(+) (limited to 'utils/plot-muons') diff --git a/utils/plot-muons b/utils/plot-muons index ab2be0e..d069d01 100755 --- a/utils/plot-muons +++ b/utils/plot-muons @@ -175,6 +175,18 @@ if __name__ == '__main__': print(stopping_muons[['run','gtid','ke','T_dx','dT','gtid_michel','r_michel','ftp_r_michel','id','r']]) + fig = plt.figure() + plot_hist2_data_mc(stopping_muons,stopping_muons_mc) + despine(fig,trim=True) + if len(muons): + plt.tight_layout() + fig.legend(handles,labels,loc='upper right') + if args.save: + plt.savefig("stopping_muons.pdf") + plt.savefig("stopping_muons.eps") + else: + plt.suptitle("Stopping Muons") + fig = plt.figure() plt.hist((stopping_muons['ke']-stopping_muons['T_dx'])*100/stopping_muons['T_dx'], bins=np.linspace(-100,100,200), histtype='step', color='C0', label="Data") plt.hist((stopping_muons_mc['ke']-stopping_muons_mc['T_dx'])*100/stopping_muons_mc['T_dx'], bins=np.linspace(-100,100,200), histtype='step', color='C1', label="Monte Carlo") -- cgit