diff options
| author | shivesh <s.p.mandalia@qmul.ac.uk> | 2018-04-10 11:12:44 -0500 |
|---|---|---|
| committer | shivesh <s.p.mandalia@qmul.ac.uk> | 2018-04-10 11:12:44 -0500 |
| commit | 01c77997f4212085a1cedc049e6c6bca98a5c1b6 (patch) | |
| tree | 670789202e71876abf46677ed91bba28dd642b17 /test/test_NSI.py | |
| parent | e1fa7270eeb9219865446cb8ceb4a5762f6aab9b (diff) | |
| download | GolemFlavor-01c77997f4212085a1cedc049e6c6bca98a5c1b6.tar.gz GolemFlavor-01c77997f4212085a1cedc049e6c6bca98a5c1b6.zip | |
updates
Diffstat (limited to 'test/test_NSI.py')
| -rw-r--r-- | test/test_NSI.py | 40 |
1 files changed, 27 insertions, 13 deletions
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) |