diff options
| author | shivesh <s.p.mandalia@qmul.ac.uk> | 2018-04-11 13:56:39 -0500 |
|---|---|---|
| committer | shivesh <s.p.mandalia@qmul.ac.uk> | 2018-04-11 13:56:39 -0500 |
| commit | bc28b9e2a31666839e837e6f0e49161713527085 (patch) | |
| tree | eef2f71d6fcc4b4bd60b71b744b33c2d94622293 /utils/gf.py | |
| parent | 326ff3bacfe0c2925afde031aa6287ebe0af0b33 (diff) | |
| download | GolemFlavor-bc28b9e2a31666839e837e6f0e49161713527085.tar.gz GolemFlavor-bc28b9e2a31666839e837e6f0e49161713527085.zip | |
GOLEMFIT takes in Haar measure params, fix. Add Bayes Factor calculation
Diffstat (limited to 'utils/gf.py')
| -rw-r--r-- | utils/gf.py | 11 |
1 files changed, 9 insertions, 2 deletions
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()]) |