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 2 H /ŝ → 1; N 3 LO QCD corrections in terms of the soft expansion; top-quark mass suppressed terms through NNLO; matching to the cross section at x → 0 through N 3 LO. 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 bb → 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. Compared to version 1.0.0 the newest SusHi version 1.6.0 now supports the 2-Higgs-DoubletModel (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@NLO can be established.
Summary of revisions:Inclusion of 2HDM, NMSSM; Improvements in the prediction of the gluon-fusion cross section: Topquark 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 heavyquark 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