From bc28b9e2a31666839e837e6f0e49161713527085 Mon Sep 17 00:00:00 2001
From: shivesh <s.p.mandalia@qmul.ac.uk>
Date: Wed, 11 Apr 2018 13:56:39 -0500
Subject: GOLEMFIT takes in Haar measure params, fix. Add Bayes Factor
 calculation

---
 utils/gf.py | 11 +++++++++--
 1 file changed, 9 insertions(+), 2 deletions(-)

(limited to 'utils/gf.py')

diff --git a/utils/gf.py b/utils/gf.py
index ebc8538..e575094 100644
--- a/utils/gf.py
+++ b/utils/gf.py
@@ -44,14 +44,21 @@ def set_up_as(fitter, params):
 
 def get_llh(fitter, params):
     fitparams = gf.FitParameters(gf.sampleTag.HESE)
-    # print params
     for parm in params:
         fitparams.__setattr__(parm.name, parm.value)
     llh = -fitter.EvalLLH(fitparams)
-    # print '=== llh = {0}'.format(llh)
     return llh
 
 
+def get_llh_freq(fitter, params):
+    print 'setting to {0}'.format(params)
+    fitparams = gf.FitParameters(gf.sampleTag.HESE)
+    for parm in params:
+        fitparams.__setattr__(parm.name, parm.value)
+    fitter.SetFitParametersSeed(fitparams)
+    return -fitter.MinLLH().likelihood
+
+
 def data_distributions(fitter):
     hdat = fitter.GetDataDistribution()
     binedges = np.asarray([fitter.GetZenithBinsData(), fitter.GetEnergyBinsData()])
-- 
cgit v1.2.3