  The NLO cross section of the gluon fusion process is matched to parton
showers in the MC@NLO approach. We work in the framework of MadGraph5_aMC@NLO
and document the inclusion of the full quark-mass dependence in the SM as well
as the state-of-the-art squark and gluino effects within the MSSM embodied in
the program SusHi. The combination of the two programs is realized by a script
which is publicly available and whose usage is detailed. We discuss the input
cards and the relevant parameter switches. One of our focuses is on the shower
scale which is specifically important for gluon-induced Higgs production,
particularly in models with enhanced Higgs-bottom Yukawa coupling.
