aboutsummaryrefslogtreecommitdiffstats
path: root/utils/fr.py
diff options
context:
space:
mode:
authorShivesh Mandalia <shivesh.mandalia@outlook.com>2020-01-08 13:08:22 -0600
committerShivesh Mandalia <shivesh.mandalia@outlook.com>2020-01-08 13:08:22 -0600
commitf6542d12ef1b4eec3ec5e8209cbe58fd5695bf06 (patch)
tree795dc4ea2dca9a378507d924b97c0cac136d65a4 /utils/fr.py
parentd1ce931f5ffa17a916c49cc287a4a423675bebe8 (diff)
downloadGolemFlavor-f6542d12ef1b4eec3ec5e8209cbe58fd5695bf06.tar.gz
GolemFlavor-f6542d12ef1b4eec3ec5e8209cbe58fd5695bf06.zip
Revert "make O traceless"
This reverts commit d1ce931f5ffa17a916c49cc287a4a423675bebe8.
Diffstat (limited to 'utils/fr.py')
-rw-r--r--utils/fr.py13
1 files changed, 7 insertions, 6 deletions
diff --git a/utils/fr.py b/utils/fr.py
index c62ce55..bf0fb56 100644
--- a/utils/fr.py
+++ b/utils/fr.py
@@ -364,15 +364,16 @@ def params_to_BSMu(bsm_angles, dim, energy, mass_eigenvalues=MASS_EIGENVALUES,
z = 0.+1e-9
if texture is Texture.OEU:
- np_s12_2, np_c13_4, np_s23_2, np_dcp, sc = 0.5, 1.0, z, z, bsm_angles
+ np_s12_2, np_c13_4, np_s23_2, np_dcp, sc2 = 0.5, 1.0, z, z, bsm_angles
elif texture is Texture.OET:
- np_s12_2, np_c13_4, np_s23_2, np_dcp, sc = z, 0.25, z, z, bsm_angles
+ np_s12_2, np_c13_4, np_s23_2, np_dcp, sc2 = z, 0.25, z, z, bsm_angles
elif texture is Texture.OUT:
- np_s12_2, np_c13_4, np_s23_2, np_dcp, sc = z, 1.0, 0.5, z, bsm_angles
+ np_s12_2, np_c13_4, np_s23_2, np_dcp, sc2 = z, 1.0, 0.5, z, bsm_angles
else:
- np_s12_2, np_c13_4, np_s23_2, np_dcp, sc = bsm_angles
+ np_s12_2, np_c13_4, np_s23_2, np_dcp, sc2 = bsm_angles
- sc = np.power(10., sc)
+ sc2 = np.power(10., sc2)
+ sc1 = sc2 / 100.
mass_matrix = np.array(
[[0, 0, 0], [0, mass_eigenvalues[0], 0], [0, 0, mass_eigenvalues[1]]]
@@ -383,7 +384,7 @@ def params_to_BSMu(bsm_angles, dim, energy, mass_eigenvalues=MASS_EIGENVALUES,
else:
NP_U = angles_to_u((np_s12_2, np_c13_4, np_s23_2, np_dcp))
SC_U = np.array(
- [[-sc, 0, 0], [0, 0, 0], [0, 0, sc]]
+ [[0, 0, 0], [0, sc1, 0], [0, 0, sc2]]
)
bsm_term = (energy**(dim-3)) * np.dot(NP_U, np.dot(SC_U, NP_U.conj().T))
bsm_ham = sm_ham + bsm_term