“…By default, we use the MSTW2008 set of PDFs [133], and we fix the renormalization and factorization scales entering the gluon-fusion cross section to µ R = µ F = m φ /2 [19,134], where φ = {h, H, A} denotes the considered Higgs boson. For bottom-quark annihilation, the central values of the scales are chosen as µ R = m φ and µ F = m φ /4 [72][73][74]135]. In the calculation of the gluon-fusion cross section we relate the bottom Yukawa coupling to the pole mass M b , computed at the three-loop level [136,137] from the input value for the running mass, m b (m b ).…”