diff options
| author | tlatorre <tlatorre@uchicago.edu> | 2020-06-23 08:36:58 -0500 |
|---|---|---|
| committer | tlatorre <tlatorre@uchicago.edu> | 2020-06-23 08:36:58 -0500 |
| commit | bbc07b5fcd3022d171bef9ced727d4cf67cf35ae (patch) | |
| tree | 14b3fa2acf438d9a1dccf31efb35264fc5f2781e /utils/plot-energy | |
| parent | a3b4f83d38c6c8aa6eb4ba7c0944fc31cb66c7ab (diff) | |
| download | sddm-bbc07b5fcd3022d171bef9ced727d4cf67cf35ae.tar.gz sddm-bbc07b5fcd3022d171bef9ced727d4cf67cf35ae.tar.bz2 sddm-bbc07b5fcd3022d171bef9ced727d4cf67cf35ae.zip | |
first draft of a script to do the chi2 analysis
This commit adds a first draft of a script called chi2. This script calculates
a chi2 for the null hypothesis test if the events in the energy range 20 MeV -
10 GeV match what we expect from atmospheric neutrino events.
Diffstat (limited to 'utils/plot-energy')
| -rwxr-xr-x | utils/plot-energy | 2 |
1 files changed, 2 insertions, 0 deletions
diff --git a/utils/plot-energy b/utils/plot-energy index 4568bcc..ef34f97 100755 --- a/utils/plot-energy +++ b/utils/plot-energy @@ -295,6 +295,8 @@ 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) |