SusHi: A program for the calculation of Higgs production in gluon fusion   and bottom-quark annihilation in the Standard Model and the MSSM
  This article describes the code SusHi (for "Supersymmetric Higgs") which calculates the cross sections $pp/p\bar{p}\rightarrow\phi+X$ in gluon fusion and bottom-quark annihilation in the SM and the MSSM, where $\phi$ is any of the neutral Higgs bosons within these models. Apart from inclusive cross sections up to NNLO QCD, differential cross sections with respect to the Higgs' transverse momentum $p_T$ and (pseudo-)rapidity $y(\eta)$ can be calculated through NLO QCD. In case of gluon fusion, SusHi contains NLO QCD contributions from the third family of quarks and squarks, NNLO corrections due to top-quarks, approximate NNLO corrections due to top-squarks, and electro-weak effects. It supports various renormalization schemes for the sbottom sector and the bottom Yukawa coupling, as well as resummation effects of higher order $\tan\beta$-enhanced sbottom contributions. SusHi provides a link to FeynHiggs for the calculation of the Higgs masses. 
Robert V. Harlander, Stefan Liebler, Hendrik Mantler
hep-ph