SusHi Bento: Beyond NNLO and the heavy-top limit

Abstract Version 1.6.0 of the code SusHi is presented. Concerning inclusive CP-even Higgs production in gluonfusion, 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 ˆ → 1 ; N 3 LO QCD corrections in terms of the soft expansion; top-quarkmass 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 b b → 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.doi.org/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@NLO can be established. Summary of revisions: Inclusion of 2HDM, NMSSM; Improvements in the prediction of the gluon-fusion cross section: Top-quarkmass 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-quarkeffective 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 gluonfusion and heavy-quark annihilation in the SM and extended Higgs sectorsthrough 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