SusHi Bento: Beyond NNLO and the heavy-top limit

Version 1.6.0 of the code SusHi is presented. Concerning inclusive CP-even Higgs production in gluon fusion, the following new features with respect to previous versions have been implemented: expansion of the partonic cross section in the soft limit, i.e. around x = M-H(2)(s) over cap -> 1; (NLO)-L-3 QCD corrections in terms of the soft expansion; top-quark mass suppressed terms through NNLO; matching to the cross section at x -> 0 through (NLO)-L-3. For CP-even and-odd scalars, an efficient evaluation of the renormalization-scale dependence is included, and effects of dimension-5 operators can be studied, which we demonstrate for the SM Higgs boson and for a CP-even scalar with a mass of 750 GeV. In addition, as a generalization of the previously available b (b) over bar -> H cross section, SusHi_1. 6. 0 provides the cross section for charged and neutral Higgs production in the annihilation of arbitrary heavy quarks. At fixed order in perturbation theory, SusHi thus allows to obtain Higgs cross-section predictions in different models to the highest precision known today. For the SM Higgs boson of M-H = 125 GeV, SusHi yields 48.28 pb for the gluon fusion cross section at the LHC at 13 TeV. Simultaneously, SusHi provides the renormalization-scale uncertainty of +/- 1.97 pb. New version program summary Program title: SusHi. Program Files doi: http://dx.doLorg/10.17632/3bb3z4k9hp.1 Licensing provisions: GNU General Public License 3 (GPL). Programming language: Fortran 77. Journal Reference of previous version: Comp. Phys. Comm. 184, Issue 6 (2013) 1605-1617 Does the new version supersede the previous version?: Yes. The new version also includes all features of previous versions. Reasons for the new version: Compared to version 1.0.0 the newest SusHi version 1.6.0 now supports the 2-Higgs-Doublet-Model (2HDM) and the next-to-minimal supersymmetric standard model (NMSSM). The effects of dimension-5 operators in the calculation of the gluon-fusion cross section can be studied. It allows to calculate the Higgs production cross section from the annihilation of heavy quarks and includes various new features which improve the gluon-fusion cross-section prediction and the associated uncertainty estimate. Links to external codes 2HDMC, MoRe SusHi and MadGraph5_aMC@NL0 can be established. Summary of revisions: Inclusion of 2HDM, NMSSM; Improvements in the prediction of the gluon-fusion cross section: Top quark mass terms up to next-to-next-to leading order, soft expansion and next-to-next-to-next-to leading order corrections in the heavy top-quark effective theory, top squark corrections up to next-to-next-to leading order; dimension-5 operators; analytic determination of the renormalization scale dependence. Inclusion of heavy-quark annihilation cross sections. Link to MoRe SusHi for the calculation of resummed transverse-momentum distributions. Nature of problem: Calculation of inclusive and exclusive Higgs production cross sections in gluon fusion and heavy-quark annihilation in the SM and extended Higgs sectors through next-to-leading order QCD, including (next to -)next-to-next-to-leading order top-(s)quark contributions and electroweak effects. Solution method: Numerical Monte Carlo integration. References: http://sushi.hepforge.org. (C) 2016 Elsevier B.V. All rights reserved.