From 4f63d9b3bc6d15dd8962b5de81ccfc78a63f3b75 Mon Sep 17 00:00:00 2001
From: tlatorre <tlatorre@uchicago.edu>
Date: Mon, 6 Jul 2020 11:20:00 -0500
Subject: 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
---
 utils/chi2 | 3 ---
 1 file changed, 3 deletions(-)

(limited to 'utils/chi2')

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)
-- 
cgit