Correct c and bmatching to be squared

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}