aboutsummaryrefslogtreecommitdiff
path: root/utils/plot-energy
diff options
context:
space:
mode:
authortlatorre <tlatorre@uchicago.edu>2020-07-06 09:41:33 -0500
committertlatorre <tlatorre@uchicago.edu>2020-07-06 09:41:33 -0500
commit35e2ebfbf1313fa753d5b18787b2d551ec243a33 (patch)
tree08dad1cf8c3d9033d77790fbfd5a7cedeafb1177 /utils/plot-energy
parentb7eda016951273419aed1f0c68dd0150a6026a3e (diff)
downloadsddm-35e2ebfbf1313fa753d5b18787b2d551ec243a33.tar.gz
sddm-35e2ebfbf1313fa753d5b18787b2d551ec243a33.tar.bz2
sddm-35e2ebfbf1313fa753d5b18787b2d551ec243a33.zip
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.
Diffstat (limited to 'utils/plot-energy')
-rwxr-xr-xutils/plot-energy2
1 files changed, 0 insertions, 2 deletions
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)