From bbc07b5fcd3022d171bef9ced727d4cf67cf35ae Mon Sep 17 00:00:00 2001 From: tlatorre Date: Tue, 23 Jun 2020 08:36:58 -0500 Subject: 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. --- utils/plot-energy | 2 ++ 1 file changed, 2 insertions(+) (limited to 'utils/plot-energy') 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) -- cgit