diff options
Diffstat (limited to 'utils/plot-muons')
| -rwxr-xr-x | utils/plot-muons | 12 |
1 files changed, 12 insertions, 0 deletions
diff --git a/utils/plot-muons b/utils/plot-muons index ab2be0e..d069d01 100755 --- a/utils/plot-muons +++ b/utils/plot-muons @@ -176,6 +176,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") plt.legend() |