Skip to content

Commit

Permalink
use matching scale instead of mb in damping factor
Browse files Browse the repository at this point in the history
  • Loading branch information
RoyStegeman committed Oct 27, 2023
1 parent 2d49d63 commit c5927eb
Showing 1 changed file with 6 additions and 8 deletions.
14 changes: 6 additions & 8 deletions src/pineko/fonll.py
Original file line number Diff line number Diff line change
Expand Up @@ -120,15 +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."""
mc = theorycard_constituent_fks["mc"]
mb = theorycard_constituent_fks["mb"]
kcthr = theorycard_constituent_fks["kcThr"]
kbthr = theorycard_constituent_fks["kbThr"]
cmatching = theorycard_constituent_fks["kcThr"] * theorycard_constituent_fks["mc"]
bmatching = theorycard_constituent_fks["kbThr"] * theorycard_constituent_fks["mb"]
q2grid = fonll_info.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]
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]
damping_factor_charm *= step_function_charm
damping_factor_bottom *= step_function_bottom
return {"mc": damping_factor_charm, "mb": damping_factor_bottom}
Expand Down

0 comments on commit c5927eb

Please sign in to comment.