From 01c77997f4212085a1cedc049e6c6bca98a5c1b6 Mon Sep 17 00:00:00 2001 From: shivesh Date: Tue, 10 Apr 2018 11:12:44 -0500 Subject: updates --- test/test_NSI.png | Bin 107718 -> 117422 bytes test/test_NSI.py | 40 +++++++++++++++++++++++++++------------- 2 files changed, 27 insertions(+), 13 deletions(-) (limited to 'test') diff --git a/test/test_NSI.png b/test/test_NSI.png index 2eb4345..000ea61 100644 Binary files a/test/test_NSI.png and b/test/test_NSI.png differ diff --git a/test/test_NSI.py b/test/test_NSI.py index 617c353..9b49c93 100644 --- a/test/test_NSI.py +++ b/test/test_NSI.py @@ -57,24 +57,38 @@ ax.step(binning, np.concatenate([[exp[0]], exp]), alpha=1, print '0.1 mutau min_llh', golem.MinLLH().likelihood print '0.1 mutau expectation', exp -# npp = gf.NewPhysicsParams() -# npp.epsilon_mutau = 0 -# # npp.epsilon_prime = 0 +np.epsilon_mutau = 0.2 -# golem.SetNewPhysicsParams(npp) +golem.SetNewPhysicsParams(npp) -# exp = np.sum(golem.GetExpectation(fit_params), axis=(0, 1, 2, 3)) -# ax.step(binning, np.concatenate([[exp[0]], exp]), alpha=1, -# drawstyle='steps-pre', label='1e-10 LV', linestyle='--') +exp = np.sum(golem.GetExpectation(fit_params), axis=(0, 1, 2, 3)) +ax.step(binning, np.concatenate([[exp[0]], exp]), alpha=1, + drawstyle='steps-pre', label='0.2 mutau', linestyle='--') -# print '1e10 LV min_llh', golem.MinLLH().likelihood -# print '1e10 LV expectation', exp +print '0.2 mutau min_llh', golem.MinLLH().likelihood +print '0.2 mutau expectation', exp -# npp = gf.NewPhysicsParams() -# npp.epsilon_mutau = 0 -# # npp.epsilon_prime = 0 +np.epsilon_mutau = 0.3 + +golem.SetNewPhysicsParams(npp) + +exp = np.sum(golem.GetExpectation(fit_params), axis=(0, 1, 2, 3)) +ax.step(binning, np.concatenate([[exp[0]], exp]), alpha=1, + drawstyle='steps-pre', label='0.3 mutau', linestyle='--') + +print '0.3 mutau min_llh', golem.MinLLH().likelihood +print '0.3 mutau expectation', exp + +np.epsilon_mutau = 0.4 + +golem.SetNewPhysicsParams(npp) + +exp = np.sum(golem.GetExpectation(fit_params), axis=(0, 1, 2, 3)) +ax.step(binning, np.concatenate([[exp[0]], exp]), alpha=1, + drawstyle='steps-pre', label='0.4 mutau', linestyle='--') -# golem.SetNewPhysicsParams(npp) +print '0.4 mutau min_llh', golem.MinLLH().likelihood +print '0.4 mutau expectation', exp ax.tick_params(axis='x', labelsize=12) ax.tick_params(axis='y', labelsize=12) -- cgit v1.2.3