Authors: Thi Nhung Dao, Julien Baglio and Margarete
A Fortran package for the calculation of the loop corrected NMSSM
Higgs boson masses and their decay widths and branching
ratios in the CP-conserving and in the CP-violating NMSSM. The already available state-of-the-art higher-order QCD corrections are supplemented by the complete next-to-leading order (SUSY-)electroweak
corrections to the neutral NMSSM Higgs boson decays into fermion and
massive gauge boson pairs, into gauge and Higgs boson final states,
chargino and neutralino pairs, and into sbottom and stop pairs. The latter also include the SUSY-QCD corrections. The one-loop corrections to the
electroweakino, stop and sbottom masses include the one-loop corrections
in two different OS and the $\DRb$ scheme. The code is based on the
Fortran package NMSSMCALC version 3.0. The code has been tested on a Linux
The Fortran code NMSSMCALCEW is derived from the code NMSSMCALC. The original
calculates the one-loop corrected and the two-loop O(alpha_t alpha_s+alpha_t^2) corrected NMSSM Higgs
boson masses in a mixed renormalisation scheme based on on-shell (OS) and DR_bar conditions of the Higgs
sector together with two possible renormalization schemes of the top/stop sector, the OS and DR_bar
NMSSMCALC calculates the Higgs decay widths and branching ratios including the state-of-the-art
higher-order QCD corrections. In NMSSMCALC, in the decays into a pair of SM down type quarks (ss_bar and
bb_bar), the dominant QCD corrections together with the large-tan(beta) corrections are resummed into
effective couplings and the running quark masses. Also the dominant SUSY-EW corrections are taken into
account through effective couplings in these decays and the decays into tau leptons. Off-shell decays into
two massive gauge bosons, gauge and Higgs boson and Higgs pair final states as well as into heavy quark
pairs are included as well.
The new code NMSSMCALCEW provides the presently
highest level of precision in the calculation of the NMSSM Higgs boson decays.
In contrast to NMSSMCALC, we have included the Z^H factor with full
in order to render the loop-corrected Higgs boson masses on-shell.
In the new code NMSSMCALCEW, the one-loop corrections to the chargino, neutralino, stop and sbottom
masses are computed for different renormalization schemes, OS and DR_bar. The NLO EW corrections
to the decays of the neutral NMSSM Higgs bosons into a pair of OS massive gauge bosons (WW,
ZZ), a Z and Higgs boson final state pair, into Higgs boson pairs and
into chargino and neutralino pairs are computed. For the Higgs boson
pair final states we include the complete one-loop corrections
(see, D.T. Nhung et al., JHEP 1311 (2013) 181 and
M. Mühlleitner et al., JHEP 1512 (2015) 034). In the decays into bb_bar and ss_bar and into up-type quarks, the remaining NLO SUSY-EW and SUSY-QCD corrections are included. For the tau taubar final states the remaining NLO SUSY-EW corrections are computed. In the decays into a stop or sbottom pair both the NLO-EW and NLO-QCD corrections are included.
The loop-induced decays into gamma gamma, Z gamma and gluon gluon are not changed with respect to NMSSMCALC 3.00. And the here described higher-order corrections in NMSSMCALCEW are only included for on-shell final states so that the off-shell decays into two massive gauge bosons, gauge and Higgs bosons and Higgs pair final states are kept the same as in NMSSMCALC 3.00. Decays of the charged Higgs boson are also the same as in NMSSMCALC 3.00.
Approximations: Both for the decays into ZZ and into WW
we include the respective one-loop squared terms in case the one-loop
contribution is larger than 80% of the tree-level decay width. For the
decays into neutralinos, in case the NLO width is negative, the improved
tree-level width, with the Z^H factor, is computed. In both cases a
warning is printed on the screen.
The program package NMSSMCALCEW is distributed in two folders. The folder
'src' contains the source files and the folder 'inc' include files. Furthermore, the file makefile is provided for the creation of the excecutable.
The folder 'src' contains the following files:
a wrap file 'nmssmcalc.f'.
'CalcMasses.F' is the main file for the computation of the loop
corrected NMSSM Higgs boson masses up to two-loop order
O(alphat * alphas + alphat^2).
'CalcMasses2loop.F' , containing the subroutines and functions necessary for the computation of the higher-order corrected Higgs boson masses. It returns the renormalised self-energy at O(alpha_t alpha_s) in the gaugeless limit at vanishing external momentum.
'CalcMasses_atat.F' for the computation of the O(alpha_t^2) corrections to the Higgs boson masses in the gaugeless limit at vanishing external momentum.
'conversion_polesmass_DRbarmass.f' for the conversion of the top-quark pole to the running DRbar mass.
'Zfactors_New.F' for the computation of the Z-factor and one-loop Higgs self-energies with complex momentum.
'bhdecay.f' ('bhdecay_c.f') for the calculation of the NMSSM Higgs boson decay widths and branching ratios in the real (complex) NMSSM in the original NMSSMCALC 3.00 implementation.
'dmb.f', 'hgaga.f', 'hsqsq_nmssm.f' (the latter only in the real NMSSM) contain additionnal decay routines.
'Xvegas.f' and 'susylha_nmssm.f' for the numerical integration and for the SUSY Les Houches Accord routines.
'EDMs.F' and 'EDMloopfunc.F' for the computation of the EDMs (in the complex NMSSM only).
It also contains the folder 'NLOEW' containing all the routines for the calculation of the NLO SUSY-EW corrections to the on-shell two-body decays of the neutral Higgs bosons. In this folder, we have th files
'bhdecay_rEWN.F' (bhdecay_cEWN.F) for the calculation of the NMSSM
Higgs boson decay widths and branching ratios in the real (complex)
NMSSM including the NLO SUSY-EW and the dominant
higher-order QCD corrections as already
implemented in NMSSMCALC for the
on-shell two-body decay widths of the
neutral Higgs bosons.
'HdecayNLO_hhh.F' for the full one-loop corrections to the Higgs-to-Higgs decays widths.
'HdecayNLO1.F' for the one-loop corrections to the masses into
electroweakinos and stops and
sbottoms (SUSY-EW and also SUSY-QCD for
the stops and sbottoms).
'HdecayNLO2.F' for the NLO-EW (NLO-QCD) corrections to decays into a pair of fermions, a Z boson and Higgs boson pair, a chargino pair, a stop and sbottom pair.
'HdecayNLO3.F' for the NLO-EW correction to the decays into a pair of neutralinos and a pair of massive gauge bosons (VV=WW, ZZ).
'HdecayNLOLibs.F' contains the bremstrahlung functions and the numerical calculation of the loop-corrected electroweakino masses.
'HdecayNLOself.F' for the computation of the self-energies for top/stop and bottom/sbottom.
The folder 'inc' contains all files for the parameter declarations which
are used in many subroutines:
'bh_EWNMSSM.in' is used in bhdecay_cEWN.F, HdecayNLO1.F,HdecayNLO2.F.
'diag-f.h' is used in CalcMasses.F.
'TLMT4.h' is used in CalcMasses_atat.F.
The folder 'inc' also contains the folder 'NLOEW' with the files:
'HdecayCTs.inc' and 'HdecayTCLs.in' are used in
The program package needs two input files named 'inp.dat' and
'bhdecay.in'. The file 'inp.dat' is read in by nmssmcalc.f. It must be
provided in the SLHA format and has to contain the blocks MODSEL, SMINPUTS
(extended by an input for the W boson pole mass, to be given as the 9th
entry), MINPAR and EXTPAR. In the complex case also IMEXPTAR and CMPLX
have to be given. The latter block is not foreseen in the original SLHA2
format. We have added it to provide the complex phase of the vacuum
expectation value for H_u. For the computation of the loop-corrected Higgs
spectrum, there are three loop levels that can be chosen by the user in
the file 'inp.dat': '1' for one-loop level, '2' for two-loop level at
O(alpha_t alpha_s), '3' for two-loop level at O(alpha_t alpha_s +
alpha_t^2). There are furthermore three options for the renormalization of
the top/stop sector: '1' and '2' for DRbar renormalization without and
with inclusion of the g1 and g2 gauge coupling contributions in the
conversion and the running. '3' for on-shell renormalization. (For details
The file 'bhdecay.in' is read in by bhdecay(_c).f and contains the CKM mixing matrix elements and several flags for the decay calculation. They are explained at the beginning of bhdecay(_c).f, in the file itself, and in the manual of NMSSMCALC. After the computation of the loop-corrected NMSSM Higgs boson masses CalcMasses.F writes out all necessary parameters, masses and mixing angles in an SLHA format file slha.in, which is read in by bhdecay(_c).f. This Fortran routine then computes the decay widths and branching ratios and writes these out in an SLHA format file slha_decay.out. Furthermore it writes out all parameters and mixing angles.
In 'bhdecay.in' that is used here by NMSSMCALCEW, new options have been included. They are
'ischhXX' to choose the renormalization scheme for the loop corrected electroweakino masses. The options are 1
(OS1), 2 (OS2) and 3 (DRbar). The two OS schemes are specified in arXiv:1907.12060.
'ischhst' to choose the renormalization scheme for the stop sector. The options are 1 (OS NLO-EW, OS NLO-QCD), 2 (OS NLO-EW, DRbar NLO-QCD), 3 (DRbar NLO-EW, DRbar NLO-QCD).
'ischhsb' to choose the renormalization scheme for the sbottom sector. The options are 1 (OS NLO-EW, OS NLO-QCD), 2 (OS NLO-EW, DRbar NLO-QCD), 3 (DRbar NLO-EW, DRbar NLO-QCD).
'iewh' to choose the levels of NLO SUSY-EW (SUSY-QCD) corrections that are
included. The options are 0 (as in NMSSMCALC 3.00), 1 (decays as
implemented in NMSSMCALC, but with the Z^H factor), 2 (full NLO
described in arXiv:1907.12060. Both for
option 1 and 2, the loop-corrected electroweakino and stop/sbottom masses
are used in the phase space factor.
How to compile and run the program package:
NMSSMCALCEW requires LoopTools 2.14 (or higher) to work with the EW
corrections in the decay part. If LoopTools is not installed on your
computer yet, go to http://www.feynarts.de/looptools/ and download the package LoopTools-2.n.tar.gz there (n>=14). In order to install LoopTools, proceed as follows:
1/ unpack the LoopTools-2.n.tar.gz archive,
2/ change into the LoopTools-2.n directory,
3/ run "./configure",
4/ run "make",
5/ run "make install".
Now download and unpack the tar file 'nmssmcalcew.tar.gz'. It
contains two subdirectories called 'nmssmcalc_rew_alphat2-master' and
'nmssmcalc_cew_alphat2-master' for the real and complex NMSSM,
respectively. Go to the subdirectory of the version in which you want to
Open in a text editor the file 'makefile' and in line 31 provide the
absolute path to your own LoopTools binary directory, located in the
main directory of LoopTools, for example
LT = /myhome/LoopTools-2.15/x86_64-Linux
Modify also the line 66 to make sure it refers to the correct 'lib'
sub-directory within the LoopTools binary directory, for example 'lib64',
STDLIBS = -L$(LT)/lib64 -looptools
In case you want to compile the package without the EW corrections in the
decay widths, you can switch the flag yesEW to FALSE on line 19 and LoopTools is not needed any more.
Now compile all files by typing 'make'. An executable 'run' is
Type './run' to run the executable. The user has to provide the input files for CalcMasses.F (default name 'inp.dat') and for bhdecay(_c).f (to be named 'bhdecay.in'). The user also has the choice to provide in the command line the names of the input and output files for CalcMasses.F (first and second argument) and the name of the output file provided by the decay routine (third argument). Hence the command will be 'run name_file1 name_file2 name_file3' in this case. Sample input files 'inp.dat' and 'bhdecay.in' are included in the .tar files.
Type 'make clean' to remove the executable as well as the object files generated in the 'obj' directory.