For calculating the cross-sections for these processes, we use MadGraph5_aMC@NLO 2.7.0 [1]} with Pythia 8.244 [2]} as parton shower (PS) employing the A14 set of tuned parameters [3]}. As model file, we use the “Higgs charaterization model” [4]}, [5]}, [6]}. We derive the cross-sections at LO in the five-flavor scheme. For the parton-distribution functions (PDF) we use the MSTW2008LO [7]} fit evaluated through the LHAPDF interface [8]}. The derived results are then rescaled to the state-of-the-art SM predictions reported in Ref. [9]}.In the case of Higgs production via gluon fusion, this procedure could be improved by taking into account the full NLO predictions available for scalar and pseudoscalar production. However, we cross-checked our predictions with the ones of the code SusHi [10]}, [11]}, finding differences which are negligible within the theoretical and experimental uncertainties.
