NMSSMCALCEW

Authors: Thi Nhung Dao, Julien Baglio, Margarete Mühlleitner, Shruti Patel and Kodai Sakurai


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 and charged Higgs boson decays into fermion and massive gauge boson pairs, into gauge and Higgs boson final states, chargino and neutralino pairs, and into third generation sfermion pairs. The decays into quarks and squarks also include the SUSY-QCD corrections. The one-loop corrections to the electroweakino, stop, sbottom, stau and sneutrino masses include the one-loop corrections in the OS and the $\DRb$ scheme, where for the electroweakino sector two different OS renormalization schemes are provided. The code is based on the Fortran package NMSSMCALC version 3.0. The code has been tested on a Linux machine.


Description:

The Fortran code NMSSMCALCEW is derived from the code NMSSMCALC. The original code NMSSMCALC 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 scheme. NMSSMCALC calculates the Higgs decay widths and branching ratios including the state-of-the-art higher-order QCD corrections. In NMSSMCALC, in the neutral Higgs boson 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. Analogously, for the charged Higgs boson the higher-order SUSY corrections have been implemented through effective couplings for the decays into fermion pairs. 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 momentum dependence 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, sbottom, stau and sneutrino masses are computed for different renormalization schemes, OS and DR_bar. In the electroweakino sector two different OS schemes have been implemented in addition to the DR_bar scheme. 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. In the decays of the charged Higgs boson we have included the SUSY-EW and, where applicable, the SUSY-QCD corrections to the decays into fermion pairs, a charged W and neutral Higgs boson final state pair, electroweakino and third generation sfermion pairs. The renormalization schemes are the same as the ones applied in the neutral Higgs boson decays.

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.


Files:

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 the 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.
  • 'ChargedHiggs_Bosons.F' for the computation of the NLO corrections to the charged Higgs boson decays into W H_i,stau-sneutrino pairs and stop-sbottom pairs
  • 'ChargedHiggs_Fermions.F' for the computation of the NLO corrections to the charged Higgs boson decays into top bottom, tau neutrino and chargino-neutralino pairs
  • 'RenConsts.F' for the computation of the renormalization constants for the stop/sbottom and sneutrino/stau sectors

    • 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 HdecayNLO1.F, HdecayNLO2.F, HdecayNLO3.F, HdecayNLOSelf.F, RenConsts.F, ChargedHiggs_Bosons.F, ChargedHiggs_Fermions.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 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 see http://arxiv.org/abs/1903.11358.)

    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).
  • 'ischstau' to choose the renormalization scheme for the stau/sneutrino sector. The options are 1 (OS) and 3 (DRbar).
  • '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 corrections as 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 work in.
    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 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.
    Type 'make clean' to remove the executable as well as the object files generated in the 'obj' directory.