  We include finite top- and bottom-mass effects in the next-to-next-to-leading
order parton shower (NNLOPS) event generator for inclusive Higgs boson
production in gluon fusion based upon the POWHEG+MiNLO approach. Since
fixed-order results for quark-mass effects only reach NLO accuracy, we add them
to the NNLOPS generator at that accuracy. We explore uncertainties related to
the unknown all-order logarithmic structure of bottom-mass effects by comparing
the assumption of full exponentiation to no exponentiation at all.
Phenomenological results showing the effects of finite quark-masses in the
NNLOPS simulation are presented. These suggest that the aforementioned
uncertainty is well contained within the envelope of plain renormalization and
factorization scale uncertainties.
