update plot-muons

- add get_multinomial_prob() function to stats.py
- add plot_hist2_data_mc() function to do the normal particle id plot
  but also print p values
- other small bug fixes

1 files changed, 0 insertions, 3 deletions

diff --git a/utils/chi2 b/utils/chi2
index 72ad89b..f1ae974 100755
--- a/utils/chi2
+++ b/utils/chi2
@@ -121,9 +121,6 @@ def make_nll(data, mc_hists):
         return nll - norm.logpdf(x[1],ENERGY_SCALE_MEAN,ENERGY_SCALE_UNCERTAINTY)
     return nll
 
-def chi2(samples,expected):
-    return np.sum((samples-expected)**2/expected,axis=-1)
-
 def get_mc_hist(data,x,bins):
     hist = np.histogram(data,bins=bins)[0]
     return sample_mc_hist(hist,norm=x[0]/100.0)