From dd852d8112fb27a037b8f27a08c2e8c61ac4749b Mon Sep 17 00:00:00 2001 From: Andrea Barontini Date: Fri, 1 Dec 2023 12:32:37 +0100 Subject: [PATCH] Correct c and bmatching to be squared --- src/pineko/fonll.py | 12 ++++++------ 1 file changed, 6 insertions(+), 6 deletions(-) diff --git a/src/pineko/fonll.py b/src/pineko/fonll.py index 213de023..aa511c29 100644 --- a/src/pineko/fonll.py +++ b/src/pineko/fonll.py @@ -120,13 +120,13 @@ def update_fk_theorycard(combined_fk, input_theorycard_path): def produce_dampings(theorycard_constituent_fks, fonll_info, damp): """Return the damping factors for each of the relevant masses.""" - cmatching = theorycard_constituent_fks["kcThr"] * theorycard_constituent_fks["mc"] - bmatching = theorycard_constituent_fks["kbThr"] * theorycard_constituent_fks["mb"] + cmatching2 = (theorycard_constituent_fks["kcThr"] * theorycard_constituent_fks["mc"]) **2 + bmatching2 = (theorycard_constituent_fks["kbThr"] * theorycard_constituent_fks["mb"]) **2 q2grid = fonll_info.Q2grid - step_function_charm = (cmatching) ** 2 < q2grid - step_function_bottom = (bmatching) ** 2 < q2grid - damping_factor_charm = (1 - cmatching / q2grid) ** damp[1] - damping_factor_bottom = (1 - bmatching / q2grid) ** damp[1] + step_function_charm = cmatching2 < q2grid + step_function_bottom = bmatching2 < q2grid + damping_factor_charm = (1 - cmatching2 / q2grid) ** damp[1] + damping_factor_bottom = (1 - bmatching2 / q2grid) ** damp[1] damping_factor_charm *= step_function_charm damping_factor_bottom *= step_function_bottom return {"mc": damping_factor_charm, "mb": damping_factor_bottom}