  A novel computation of the fully-differential cross section for neutral
Higgs-boson production through gluon fusion in the CP-conserving NMSSM is
presented. Based on the calculation of NLO corrections to the total cross
section [1], we implemented the NMSSM amplitudes in three codes, applying
different resummation techniques: analytic transverse-momentum resummation at
NLO+NLL, and two fully-differential NLO+PS Monte-Carlo approaches using the
MC@NLO and POWHEG matching procedures, respectively. We study phenomenological
predictions for distributions in the NMSSM with a special emphasis on the Higgs
transverse-momentum spectrum. Reasonable agreement among the various approaches
is found, once well-motivated choices for the unphysical matching scales and
the determination of the related uncertainties are made.
