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 DR_bar conditions. In the two-loop O(alpha_t alpha_s) corrections the user can choose between two 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 the decays into stop and sbottom pairs, respectively, contain the SUSY-QCD corrections. 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 are: 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 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.

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 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.

Last modified: Thur Aug 13 22:45:31 CET 2015