From 2d49d6360bdba38b72e2fa8d529638e29a73a70a Mon Sep 17 00:00:00 2001 From: RoyStegeman Date: Fri, 27 Oct 2023 09:36:59 +0100 Subject: [PATCH] account for k{fl}Thr in damping step function fonll --- src/pineko/fonll.py | 6 ++++-- 1 file changed, 4 insertions(+), 2 deletions(-) diff --git a/src/pineko/fonll.py b/src/pineko/fonll.py index 3075e227..ba087d08 100644 --- a/src/pineko/fonll.py +++ b/src/pineko/fonll.py @@ -122,9 +122,11 @@ def produce_dampings(theorycard_constituent_fks, fonll_info, damp): """Return the damping factors for each of the relevant masses.""" mc = theorycard_constituent_fks["mc"] mb = theorycard_constituent_fks["mb"] + kcthr = theorycard_constituent_fks["kcThr"] + kbthr = theorycard_constituent_fks["kbThr"] q2grid = fonll_info.Q2grid - step_function_charm = mc**2 < q2grid - step_function_bottom = mb**2 < q2grid + step_function_charm = (kcthr * mc) ** 2 < q2grid + step_function_bottom = (kbthr * mb) ** 2 < q2grid damping_factor_charm = (1 - mc / q2grid) ** damp[1] damping_factor_bottom = (1 - mb / q2grid) ** damp[1] damping_factor_charm *= step_function_charm