Distributions for neutral Higgs production in the NMSSM

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.