Released by Julien Baglio, Ramona Gröber, Margarete
Mühlleitner, Dao Thi Nhung, Heidi Rzehak, Michael Spira,
Juraj Streicher and Kathrin Walz
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 Fortran package NMSSMCALC computes the one-loop corrected and the
two-loop O(alpha_t alpha_s) corrected NMSSM
Higgs boson masses in a mixed renormalisation scheme based on on-shell and
In the two-loop O(alpha_t alpha_s) corrections the user can choose
renormalisation schemes of the top/stop-sector: a DR_bar scheme and an OS
scheme. NMSSMCALC calculates the Higgs decay widths and branching
ratios including the dominant higher order QCD corrections. The decays
into fermion pairs with down-type final state leptons include the higher
order SUSY-electroweak corrections, the ones into down-type final state
quarks the SUSY-QCD corrections. In case of b-quark(s) in the final state
both SUSY-QCD and SUSY-electroweak corrections have been implemented.
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. In the real NMSSM
decays into stop and sbottom pairs, respectively, contain the
The input and output files feature the SUSY Les Houches Accord (SLHA).
For the complex NMSSM, since the version 2.00 (13 August 2015) also the
EDMs are calculated. The user can choose to calculate the EDMs via the
switch 10 in the BLOCK MODSEL. The meanings of the second switch there
0 - no EDMs are calculated, this is the default setting; 1 - the EDMs are
calculated; 2 - the individual contributions of the EDMs are given out.
NMSSMCALC consists of a wrap file
nmssmcalc.f and three main files: CalcMasses.F for the computation of the loop
corrected NMSSM Higgs boson masses up to two-loop O(alpha_t alpha_s) in
the real and in the complex NMSSM; CalcMasses2loop.F which contains necessary
subroutines and functions and returns the renormalised self-energy at
O(alpha_t alpha_s); EDMs.F and EDMloopfunc.F for the computation of the EDMs in
the complex NMSSM;
bhdecay.f for the calculation of the NMSSM
Higgs boson decay widths and branching ratios in the real NMSSM; bhdecay_c.f for the calculation of the NMSSM
Higgs boson decay widths and branching ratios in the complex NMSSM.
Further linked files are
dmb.f, hgaga.f, hsqsq.f (only in the real NMSSM),
Xvegas.f and susylha_nmssm.f.
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 added it to provide the
complex phase of the vacuum expectation value for H_u. 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 and in the manual. After the
computation of the loop-corrected NMSSM Higgs boson masses CalcMasses.F
out all necessary parameters, masses and mixing angles in an SLHA format
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.
How to compile and run the program package:
Download and unpack the tar file nmssmcalc.tar.gz for the real NMSSM or
nmssmcalc_c.tar.gz for the complex NMSSM. Compile all files by typing
'make'. An exectuable 'run' is created. Type './run' to run the
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.
Last modified: Thur Aug 13 22:45:31 CET 2015