Merge pull request #134 from NNPDF/fix_thresholds_nfonll

account for k{fl}Thr in damping step function fonll

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."""
-    mc = theorycard_constituent_fks["mc"]
-    mb = theorycard_constituent_fks["mb"]
+    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 = mc**2 < q2grid
-    step_function_bottom = 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}