diff options
| author | tlatorre <tlatorre@uchicago.edu> | 2020-07-06 12:27:27 -0500 |
|---|---|---|
| committer | tlatorre <tlatorre@uchicago.edu> | 2020-07-06 12:27:27 -0500 |
| commit | e79e1fcf612c386f5f808bfa9099af723fc28a76 (patch) | |
| tree | 75c93ff9990521b0cf5ac15b7104e5246054abb9 | |
| parent | d48be389e496801d48af305fad67c0d79bb7585a (diff) | |
| download | sddm-e79e1fcf612c386f5f808bfa9099af723fc28a76.tar.gz sddm-e79e1fcf612c386f5f808bfa9099af723fc28a76.tar.bz2 sddm-e79e1fcf612c386f5f808bfa9099af723fc28a76.zip | |
update plot-muons to plot energy distributions for each particle ID for stopping 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() |