We present the results of our recent study of the theoretical uncertainties that affect the predictions for the Higgs-boson transverse-momentum in gluon fusion when fixed- and all-order
results are matched. Our investigation consists of a twofold analysis: first we present a detailed comparison of two recently introduced prescriptions for the determination of the matching scale, then we apply the results of these methods to three widely used matching frameworks, namely the aMC@NLO and POWHEG Monte Carlo approaches and
analytic resummation. The results of our study are applied to the production of the SM Higgs boson and of the neutral Higgs bosons
of the Two-Higgs-Doublet Model in a variety of scenarios.