NMSSMTools
TOOLS FOR THE CALCULATION OF THE
HIGGS AND SPARTICLE SPECTRUM IN THE NMSSM:
NMHDECAY, NMSPEC and NMGMSB
NMHDECAY
The Fortran code NMHDECAY computes the masses, couplings and decay
widths of all Higgs bosons of the NMSSM, and the masses of all
sparticles, in terms of its parameters at the electroweak (or susy
breaking) scale: the Yukawa couplings lambda and kappa, the soft
trilinear terms A_lambda and A_kappa, and tan(beta) and mu_eff =
lambda*S. (Instead of A_lambda, the MSSMlike parameter M_A can also
be used as input.) The computation of the Higgs spectrum includes
leading electroweak corrections, two loop terms and propagator
corrections. The computation of the decay widths is carried out as in
HDECAY, but
momentarily without three body decays. Each point in parameter space is
checked against negative Higgs boson searches at LEP, and negative
sparticle searches at LEP and the Tevatron, including unconventional
channels relevant for the NMSSM. B physics constraints from b > s
gamma, Delta M_q, B > mu+mu and B+ > tau+ nu_tau are included as in
ref. [4] below and updated occasionally. The dark matter relic density,
direct and indirect detection cross sections can be computed via a
link to a NMSSM version of the MicrOMEGAs code [3].
SLHA
conventions for input and output are used. NMHDECAY is part of the
NMSSMTools package that can be downloaded below.
* Laboratoire de Physique Theorique, Universite de Paris XI,
Batiment 210, F91405 Orsay, France.
** Davis Institute for High Energy Physics, University of California,
Davis, California 95616, USA.
*** Laboratoire Physique Theorique et Astroparticules,
Univ. Montpellier II, F34095 Montpellier, France
References:
[3] G. Belanger, F. Boudjema, C. Hugonie, A. Pukhov, A. Semenov,
JCAP 0509:001 (2005),
arXiv:hepph/0505142
NMSPEC
NMSPEC allows to compute the sparticle and Higgs
masses, Higgs decay widths and couplings in the NMSSM, with soft Susy
breaking terms specified at the GUT scale. Exceptions are the soft
singlet mass and kappa, that are both determined in terms of the other
parameters through the minimization equations of the Higgs potential.
The soft Higgs and gaugino masses at the GUT scale can be chosen as
nonuniversal, if desired. NMSPEC is included in the NMSSMTools package
below.
* Laboratoire de Physique Theorique, Universite de Paris XI,
Batiment 210, F91405 Orsay, France.
** Laboratoire Physique Theorique et Astroparticules,
Univ. Montpellier II, F34095 Montpellier, France
References:
U. Ellwanger and C. Hugonie, "NMSPEC: A Fortran code for the
sparticle and Higgs masses in the NMSSM with GUT scale boundary
conditions", Comput.Phys.Commun.177:399407,2007,
arXiv:hepph/0612134
NMGMSB
NMGMSB allows to compute the sparticle and Higgs masses, Higgs decay
widths and couplings in the NMSSM, where soft Susy breaking terms are
specified at by gauge mediated supersymmetry breaking, plus terms
originating from couplings of the Singlet to messengers.
The boundary conditions at the messenger scale are taken from the paper
below, and described in the file README. Either the soft Singlet mass or
a Susy tadpole term XIS are determined in terms of the other parameters
through the minimization equations of the Higgs potential. Sample input
and output files are included. It is part of the NMSSMTools package
described below.
* Laboratoire de Physique Theorique, Universite de Paris XI,
Batiment 210, F91405 Orsay, France
References:
U. Ellwanger, C.C. JeanLouis, A.M. Teixeira,
"Phenomenology of the General NMSSM with Gauge Mediated
Supersymmetry Breaking",
JHEP 0805:044 (2008),
arXiv:0803.2962 [hepph].
NMSDECAY
NMSDECAY allows to compute sparticle widths and branching ratios.
It is based on a generalization of SDECAY, including
the corresponding QCD corrections and 3body decay modes. Slepton 3body
decays, possibly relevant in case of a singlinolike LSP, have been
added.
* Laboratoire de Physique Theorique, Universite de Paris XI,
Batiment 210, F91405 Orsay, France
** Laboratoire de Physique Corpusculaire, CNRS/IN2P3  UMR
6533, Campus des Cezeaux, 24 Av. des Landais, F63171 Aubiere Cedex,
France
References:
Debottam Das, Ulrich Ellwanger, Ana M. Teixeira,
"NMSDECAY: A Fortran Code for Supersymmetric Particle Decays
in the NexttoMinimal Supersymmetric Standard Model"
arXiv:1106.5633 [hepph], Comput.Phys.Commun. 183 (2012) 774779
SHOULD BE CITED TOGETHER WITH
M. Muhlleitner (Karlsruhe, Inst. Technol.), A. Djouadi (Orsay, LPT
& CERN, Theory Division), Y. Mambrini (Orsay, LPT),
"SDECAY: A Fortran code for the decays of the supersymmetric particles
in the MSSM", Comput.Phys.Commun.168:4670 (2005),
arXiv:hepph/0311167.
NMSSMTools
NMSSMTools contains NMHDECAY, NMSPEC, NMGMSB and NMSDECAY.
* Laboratoire de Physique Theorique, Universite de Paris XI,
Batiment 210, F91405 Orsay, France.
** Laboratoire Physique Theorique et Astroparticules,
Univ. Montpellier II, F34095 Montpellier, France
Downloads:
HOW TO USE NMSSMTOOLS: README (plain text)
Directory NMSSMTools_4.9.3 (gzipped
and tared), contains NMHDECAY, NMSDECAY, NMSPEC, NMGMSB,
micromegas_3.3, EXPCON for experimental constraints, and sample input
and output files.
How to use NMSSMTools:
After the download of NMSSMTools_x.tgz (where x is the version
number), type "tar zxvf NMSSMTools_x.tgz". This will create a directory
named NMSSMTools_x.
Go into this directory.
There you find the directories
"sources", "main" and "micromegas", which contain source files;
"EXPCON", which contains data files corresponding to experimental
constraints;
"SAMPLES", which contains sample input and output files;
COMPILATION:
To compile, type first "make init", then "make". A first compilation
may take a while, since most subroutines of micromegas are compiled.
The following eight executable routines are created in the directory
"main": nmhdecay, nmhdecay_rand, nmhdecay_grid, nmspec, nmspec_rand,
nmspec_grid, nmgmsb and nmgmsb_rand.
If a subroutine in the directory "sources" was modified, one has to
type "make init" and "make" again. If a routine in the directory "main"
was modified, it suffices to type "make" again.
To delete all the already compiled codes type "make clean".
INPUT FILES:
From the version 2.0.0 onwards, the task to be performed by an input
file is independent from its name. Any name is allowed, provided it
contains the three letters "inp"; it can be of the general form
PREFIXinpSUFFIX where PREFIX and SUFFIX can contain dots etc..
The input file can be located in any directory specified by a PATH.
To run any input file PREFIXinpSUFFIX, type "run PATH/PREFIXinpSUFFIX".
(PATH is optional; if absent, the input file has to be located in the
same directory as the script file "run".)
The output files are located in the directory specified by PATH.
They have the following format:
If one single point in the parameter space is evaluated:
PREFIXspectrSUFFIX (includes the previous
PREFIXdecaySUFFIX),
PREFIXomegaSUFFIX (if the relic density is computed)
If scans are performed:
PREFIXerrSUFFIX as well as PREFIXoutSUFFIX
However, the task to be performed by an input file must be specified in
the BLOCK MODSEL at the beginning (see the SLHA2 conventions in
B. Allanach et al., SUSY Les Houches Accord 2, Comp.Phys.Commun.180:825,2009,
arXiv:0801.0045.
For more details see the description in README .
Latest update: 30.06.2016
History:
30.06.2016: version 4.9.3
 Minor update of bsg.f, with thanks to Florian Domingo.
17.05.2016: version 4.9.2
 Constraints from CMSPASHIG13010 on H_125>AA>4mu corrected (the bounds were too low
by a factor 1/100), with thanks to Robin Cameron Aggleton.
18.03.2016: version 4.9.1
 Gravitino relic density for nmspec_rand and nmspec_mcmc and the output
of XIFGUT, XISGUT in nmspec_grid.f, nmspec_rand.f and nmspec_mcmc.f corrected, with thanks to YongGon Yi.
04.03.2016: version 4.9.0
 New tests of Kphysics: K+> Pi+ nu nubar, KL> Pi0 nu nubar, DMK, epsK, thanks to Florian Domingo.
 Bphysics: Updates of previous and addition of new processes:
B>Xs mu+mu, B>Xs mu+mu, b > d gamma, Bd>mu+mu, B> Xs nu nubar, B+> K+ nu nubar,
B> Ks nu nubar, BR[B+ > D tau+ nu_tau]/BR[B+ > D l+ nu_l], BR[B+ > D* tau+ nu_tau]/BR[B+ > D* l+ nu_l],
thanks to Florian Domingo.
 Tests of SM Higgs couplings: Updates of parameters (based on Lilith1.1.3 with database 15.09, thanks to J. Bernon).
 Decays of light scalars: Muon/strange quark loops added to H/A > gammagamma (but BRs not reliable for M_H/A < 1 GeV).
 msferm.f: (minor) bugs corrected, thanks to Florian Domingo.
<
 The gravitino mass and, for flag I7=2, branching ratios into the gravitino are computed.
 Gauge mediation: for flag I3=/=0, the nonthermal relic density due to the gravitino is computed by MicrOmegas from the gravitino/NLSP mass ratio.
01.10.2015: version 4.8.2
 In bsg.f: Delta M_d,s also upgraded to NLO, numerically unstable functions corrected
(thanks to Daniele Barducci for spotting), with thanks to Florian Domingo,
see arXiv:1512.02091 [hepph].
23.09.2015: version 4.8.1
 Added missing files to the update of bsg.f in 4.8.0.
 Missing (rare) HiggstoHiggs decays are printed in the output, with thanks to Marcin Badziak.
09.09.2015: version 4.8.0
 Major update of bsg.f, with thanks to Florian Domingo,
see arXiv:1512.02091 [hepph]:
Most Bphysics observables are upgraded to full NLO;
recent SM results of 1503.01789 (B > Xs gamma), 1311.0903 (Bs > mu+mu) and 1503.04849 (B > Xsl+l) have been taken into account. Experimental and lattice input from current HFAG and FLAG data.
13.08.2015: version 4.7.1
 Bug in lhchig.f concerning chi2gam/chi2zz in case of two scalars near 125 GeV corrected, with thanks to
NilsErik Bomark
02.07.2015: version 4.7.0
 For Higgs mass precision >=1, running stop/sbottom masses and trilin. couplings at the scale QSTSB are
computed without threshold corrections such that they are defined in DR_bar (in order to comply with SLHA).
 The SUSY scale can now be defined in BLOCK EXTPAR, item 0 (in order to comply with SLHA).
 Yukawa corrections ~ht, hb to the stop/sbottom pole masses added
 Lmunu in runpar takes into account the Z_3 violating contribution MUP to the singlino mass
(with thanks to F. Domingo).
 Chi2 of the signal rates of the SMlike Higgs boson are computed for the CPeven Higgs with mass
closest to 125.1 GeV, indep. of its mass
 Gauge and SM Yukawa couplings at the SUSY scale are computed in an extra routine "getsusycoup" attached to
runpar.f
22.04.2015: version 4.6.0
 NMSDECAY can be called for GMSB models, decays into a gravitino are included
 More general GMSB models of the type discussed by A. Delgado, G. F. Giudice and P. Slavich in
PLB 653, 424 (2007) arXiv:0706.3873 can be defined.
 Correction of bugs concerning selectrons/smuons in NMSDECAY
(thanks to Daniele Barducci and Nishita Desai for spotting)
 Bugs in decays for very light (< 1 GeV) Higgs bosons corrected
 NMSPEC improved for very large (> 10 TeV) masses at the GUT scale
 Optional output without privately defined SLHABLOCK's (and some additional BLOCK's);
useful as param_card.dat for MadGraph as indicated in README
25.02.2015: version 4.5.1
 Correction of a bug in the compilation of micrOMEGAs on Mac (with thanks to A. Pukhov)
 Improvement of the convergence algorithm in NMSPEC.
 Contributions of two distinct Higgs bosons to the gammagamma/ZZ signal rates are summed only if the mass difference is below 3 GeV.
 Constraint on the invisible Zwidth reduced to 0.5 MeV due to update of the SM contribution in A. Freitas, 1401.2447.
18.12.2014: version 4.5.0
 Calculation of relic density in MicrOmegas made possible for LSP masses below 1 GeV (with thanks to A. Pukhov)
 Scalar and pseudoscalar Higgs masses below 1 GeV made possible (branching ratios not quite reliable, however)
 The switch 1 in BLOCK MODSEL allows to insist on the Z_3invariant NMSSM for the values 3,4,5
21.11.2014: version 4.4.1
 Hadronic parameters for Bphysics are updated, with thanks to F. Domingo. For details see the file bsg.f.
 Calculation of the relic density in MicrOmegas corrected for nonZ_3invariant terms,
with thanks to J. Gainer et al. for raising the issue, and A. Semenov and A. Pukhov for the corrections.
 Calculation of the relic density in MicrOmegas improved in cases of schannel resonances, with thanks to A. Pukhov.
 Minor bug (sign) in mgluino.f corrected, impact on the gluino mass <~2 GeV.
15.09.2014: version 4.4.0
 Constraints from Bphysics are updated using values from the HFAG web page. For details see the file bsg.f.
 The default value of the SMlike Higgs mass is 125.1 GeV. Constraints on signal strengths are updated, with
thanks to J. Bernon, B. Dumont and S. Kraml.
 Constraints on H>tautau from ATLASCONF2014049 added to the ones from CMS.
 Constraints from CMSPASHIG13010 on H_125>AA>4mu added.
 Regularisation of IR logarithm in mhiggs.f for very a light singlino and Higgs precision 0, with thanks to M. CahillRowley, J. Gainer, J. Hewett and T. Rizzo.
 Regularisation of some IR logarithms for a very light pseudoscalar, and of some potential divergencies for very heavy squarks for Higgs precision 12.
 The use of Higgs precision 12 is blocked in the presence of Z_3 violating terms, since this part of the code
does not include Z_3 violating terms.
26.06.2014: version 4.3.0
 Subroutines for Markov Chain Monte Carlo scans in parameter space are
provided, separately for the general NMSSM (NMHDECAY_MCMC), the NMSSM
with parameters defined at the GUT scale (NMSPEC_MCMC) and the NMSSM
with gauge mediated supersymmetry breaking (NMGMSB_MCMC). For
applications see README and the files provided in the directory
SAMPLES.
 corrections in the Higgs mass calculation (Higgs mass precision 1 and 2):
 correction of the two loop term in the MSbarDRbar transition of
h_top (with thanks to M. Muehlleitner and P. Slavich)
 the light Higgs contribution to the transition from the top pole
mass to the running top mass is expressed in terms of h_top at m_top and not M_SUSY;
increases the SMlike Higgs mass by 0.20.3 GeV (with thanks to P. Slavich)
08.01.2014: Commented checks for PROBs in nmhdecay.rand uncommented (no change of version number)
10.12.2013: version 4.2.1
 BLOCKs in the files *spectr* for input for HiggsBounds corrected
 (minor) bug in ftpar.f corrected, with thanks to S. Trojanowski
6.12.2013: version 4.2.0
 LUX constraints on the direct detection cross section implemented
 Constraints from CMSPASHIG12050 on H>tautau implemented (quite
conservative for MA<~250 GeV, but realistic for MA>~300 GeV)
 Higgs decays into two offshell Z/Wbosons can be included via a flag
(from Hdecay, http://people.web.psi.ch/spira/hdecay/)
 Bug in the reduced couplings of CPodd Higgses to gluons/photons
corrected (the branching fractions were not affected)
 The reduced couplings of Higgses to bquarks include now the DELMB corrections
(also used in loop induced couplings)
 The default value of the bquark pole mass in init.f is set to 4.54 GeV
 Higgs signal rates into invisible final states (pairs of neutralinos)
indicated in the output
 In the output *spectr*: BLOCKs with parameters at the GUT scale renamed
 The content of the previous output files *decay* is appended to the output
files *spectr*
 The output files *spectr* can be used as input for
HiggsBounds from http://higgsbounds.hepforge.org/
(BLOCKs added)
 For Mac users: new compiler (clang) used for compatibility with OSX 10.9.
18.10.2013: version 4.1.2
 bug in sources/fullhiggs.f corrected. Leads to an increase of the SMlike
Higgs mass for PFLAG=1 or 2 and MA >> mtop.
 minor corrections in MicrOmegas and the files sources/full* in order to avoid
NaNs for unphysical parameter sets.
27.08.2013: version 4.1.1
 bug in the total width of stau2 in NS_sleptons.f corrected, with thanks to L. Duflot.
 bug in ftpar.f ~MUP corrected, with thanks to Doyoun Kim.
21.07.2013: version 4.1.0
 constraints on a SMlike Higgs mass (125.7 \pm 3 GeV) and its signal rates
(from arXiv 1306.2941 ) added;
modifiable in the subroutine Higgs_CHI2 at the end of sources/lhchig.f.
 guess for the initial values of the parameters in nmspec_rand.f corrected; with
thanks to Jiang Yun.
23.05.2013: version 4.0.0
 The general NMSSM with nonZ_3 invariant parameters either at the
(low) SUSY scale or at the (high) GUT scale is made possible, as well as
the nMSSM with kappa=0.
(See README for the allowed options; now most
parameters not given in input files are assumed to be = 0).
 NMSSM specific radiative corrections to
CPodd and charged Higgs masses from chargino/neutralino loops have been
added for "Higgs precision = 0".
 Chargino/neutralino 1loop corrections to
the CPeven Higgs masses use parameters at the scale QSTSB (resumming
some large logarithms), and some NMSSM specific large logarithms have
been added for "Higgs precision = 0".
 The relic density and the direct detection cross sections are
computed with the new version
micromegas 3.0
(as of May 2013),
which includes virtual contributions to neutralino annihilation
cross sections (e.g. neutralinoneutralino > WW)
and loop contributions to neutralino+neutralino > gamma+gamma.
(Hence a first compilation of the code and a calculation of the relic
density takes longer than before. In relden.F, virtual decays can be
switched off setting setVVdecay(0) instead of the default setVVdecay(1).)
 The direct detection cross section is computed in terms of sigma_PiN
and sigma_S instead of sigma_PiN and sigma_0; now sigma_PiN and sigma_S
can be specified in init.f.
 The exp. constraint on dark matter relic density is now 0.1187 +/
0.0017 (arXiv:1303.5076),
but one has to allow for theoretical errors which are difficult to
estimate. The default window in init.f is 0.1187 +/ 0.0085, but can be
modified by the user.
 New interface between NMSSMTools and\A0
micromegas 3.0
using the full
Higgs effective lagrangian with parameters computed at the same
precision level as for Higgs masses and couplings. These
parameters can be found in the\A0BLOCK EFFECTIVE_COUPLINGS at the end of the output file
#spectr#.
 With "Higgs precision = 1 or 2", the relic density and the direct
detection cross sections are computed by micromegas with the correct
Higgs masses (previously micromegas used always the Higgs masses as
computed with "Higgs precision = 0").
 The diagonal chargino loop contributions to H/A>Z+gamma are included
(with thanks to Peiwen Wu for discussions).
30.01.2013: version 3.2.4
 Factor 2 in the charged Higgs loop contribution to H>gammagamma and
H>Zgamma corrected (with thanks to Peiwen Wu, L. Sluter and P.
Richardson; impact on the partial widths <2% for a charged Higgs mass >~
300 GeV)
11.01.2013: version 3.2.3
 XENON100 bounds corrected, with thanks to pwwu@itp.ac.cn
 lhchig.f updated: more reduced significances for Higgs bosons at the
LHC added, but LHC bounds omitted (left to the user); notation of the
reduced significances in the ouput files corrected (with thanks to C.
Beskidt)
21.12.2012: version 3.2.2
 Constraints from Bphysics updated
 Constraints from XENON100 updated
 Bugs in tdecay concerning t > H+ + b
(decays near the threshold, and exp. constraints) corrected
(with thanks to Neil Christensen and Zhen Liu)
 Convergence of the charged Higgs mass calculation in fullhig insured
(with thanks to Jinmian Li)
12.7.2012: version 3.2.1
 Provisional suppression of LHC constraints on Higgs signal cross
sections (to be updated)
23.3.2012: version 3.2.0
 In NMSPEC: improved convergence of RGEs with boundary conditions at
the Susy and GUT scales for large Yukawa couplings
 In NMSPEC: mu_eff and kappa at the Susy scale can be given as input,
then MHD and MHU at the GUT scale become outputs
 bounds from XENON100 on the direct detection cross section can be
applied
 LHCb constraints on BR(Bs>mu+mu) applied
 LHC constraints on Higgs signal cross sections updated
 code of M3H^2 in the output corrected (with thanks to B. Fuks)
16.1.2012: version 3.1.0
 Constraints from single (not combined) Higgs search channels
H>tautau, H>bb, H>ZZ, H>WW, top>bH+ >tau+nu (CMS) and
H>2gamma (ATLAS) added
 Compilation improved to be compatible with Fedora
 Improvement of the numerical integration of the Boltzmann equations
inside MicrOmegas
 CPodd Higgs mixing corrected (for PFLAG > 0 only)
 The finetuning with respect to the input parameters is computed
following arXiv:1107.2472 [hepph]
(The command "CALL FTPAR(PAR,...)" in MAIN is commented for scans)
7.10.2011: version 3.0.2
 NMSDECAY: a warning appears on the screen, of QCD corrections to colored
sparticle decay widths are negative and larger than the tree level
widths. Then the relative corrections "corr", leading to a negative
width when added according to (1+corr), are replaced by 1/(1corr).
 the measured (not computed) W and Zwidths are added to the output
of NMSDECAY
 the BLOCKs for parameters at the GUTscale in the SLHAoutput of
NMHDECAY and NMSPEC are renamed
 Higgssignificances at the LHC are computed only if a corresponding
flag "7" in the input file is switched on
12.9.2011: version 3.0.1
 NMSDECAY corrected, with special thanks to Oscar Stal:
negative chargino masses now considered correctly;
output corrected
 RGE for the tadpole term xi_S corrected, with special thanks to Renato
Fonseca
 constraints from B_s > mu mu updated (comb. LHCb + CMS)
 constraints from H/A > tau tau updated according to CMSPASHIG11020
29.6.2011: Release of version 3.0.0

NMSDECAY is added, which allows to compute sparticle widths and
branching ratios (see above). It is based on a generalization of
SDECAY, including
the corresponding QCD corrections and 3body decay modes.
With special thanks to M. Muhlleitner, A. Djouadi and Y. Mambrini.
 gfortran is set as default compiler for LINUX. f77 is no longer
recommended due to potential problems with complex functions in NMSDECAY.
17.5.2011: version 2.3.7
 charged Higgs mass with complete 1loop and Yukawa induced 2loop
corrections (if improved corrections to Higgs masses are switched on),
with special thanks to K.H. Phan and P. Slavich
11.5.2011: bug in H>bb corrected
22.4.2011: version 2.3.6
 LHC constraints on H/A > 2 tau included
 constraints from ALEPH on H>AA>4tau removed for MA > 9.4 GeV (due to Aeta_b mixing, see arXiv:1105.1722 [hepph])
 decays H/A > bb allowed only for MH/MA > 2*MB ~ 10.6 GeV
4.2.2011: version 2.3.5
 top quark decays into charged Higgs and top squarks + chi_0 included
 Tevatron constraints on t > b + H+, H+ > c+s or tau+nu_tau or A1+A1 included
4.1.2011: version 2.3.4
 updated constraints on B_s > mu+ mu and B+ > tau+ nu_tau
13.12.2010: version 2.3.3
 corrected constraints from ALEPH (JHEP 1005:049,2010, arXiv:1003.0705)
14.10.2010: version 2.3.2
 the improved corrections to the Higgs masses as in G. Degrassi, P. Slavich,
Nucl.Phys.B825:119150,2010,
arXiv 0907.4682 ,
can be switched on (without or with
pole mass corrections), with special thanks to P. Slavich
 improved calculation of h_top(Q_Susy)/Zfactors for rad. corrections
 bug in the 2loop RGEs for A_lambda corrected (with thanks to F. Staub)
 bug in the BRs(A > sfermions) corrected (with thanks to L. Suter)
 normalisation of the partial width H > W+H corrected (with thanks to L. Suter)
 slight correction of the function AQCD in decay.f
 constraints from ALEPH (JHEP 1005:049,2010, arXiv:1003.0705) on H>AA>4tau included
 corrections in the output SLHA: Higgs vev v(Q) properly normalised, smuons indexed by L,R
18.11.2009: version 2.3.1
 benchmark points corrected
 small correction in the radiative contributions in decay.f
(for M_Higgs > M_stop)
 output format of the GUT/messenger scale corrected
10.11.2009: Release of version 2.3.0
 constraints from B_s > X_s mu+ mu added (with thanks to F. Domingo)
 algorithm for the iterative determination of mu_eff in nmspec improved
 conventions for the Susy and soft singlet masses mu' and m_S^2' in
nmgmsb modified according to the new version of SLHA2
(additional factors 1/2)
 small corrections in the radiative contributions in decay.f, lowmuk.f,
checkmin.f and the two loop RGE for h_t (> slightly modified benchmark points)
31.3.2009: Release of version 2.2.0
 Two loop RGEs for all couplings and soft terms including the
NMSSMspecific contributions completed, pole mass corrections for the
CPodd and charged Higgs scalars improved.
 Corresponding small changes
of the numerical values of the benchmark points.
6.2.2009: version 2.1.1
 Version "grid" added to NMGMSB
 Benchmark point P1 in BMPGMSB corrected
10.12.2008: Release of version 2.1.0
 Significances for the detection of all neutral Higgs bosons at the
LHC (at low luminosity, 30 fb1, and at high luminosity, 300 fb1) can
be computed, following U. Ellwanger, J.F. Gunion and C. Hugonie, JHEP
{0507} (2005) 041,
arXiv:hepph/0503203, and G. Belanger, C. Hugonie and A. Pukhov, arXiv:0811.3224 [hepph].
 New version of micrOMEGAS (2.2) included, which allows to compute the
LSP annihilation cross section, the resulting photon spectrum for
indirect detection and spin dependent/spin independent LSPnucleon cross
sections for direct detection, following G. Belanger, F. Boudjema, A.
Pukhov and A. Semenov,
arXiv:0803.2360 [hepph], and G. Belanger, C. Hugonie and A. Pukhov,
arXiv:0811.3224 [hepph].
 Constraints from radiative Upsilon decays and the eta_b mass
measurement are implemented as in F. Domingo, U. Ellwanger, E. Fullana,
C. Hugonie and M.A. SanchisLozano, arXiv:0810.4736 [hepph].
 Initial guess for the parameters at the scale Q2 in the routines
NMGMSB corrected (with thanks to Hisano, Junji)
30.5.2008: version 2.0.1
 Bug in the rad. corr. to M2 in mneu.f corrected (with thanks to Junjie
Cao)
 Output of the BLOCK HMIX corrected in order to coincide with SLHA2
conventions (with thanks to A. Pukhov)
 Inappropriate floats in bsg.f corrected (with thanks to A. Pukhov)
 Warnings in case of inappropriate input files modified
13.05.2008: Release of version 2.0.0
 The architecture of the code has been modified: From the version
2.0.0 onwards, the task to be performed by an input file is independent
from its name. However, the task to be performed by an input file must
be specified in the BLOCK MODSEL at the beginning, see the new
description README.
 The routines NMGMSB have been added, which allow for
boundary conditions specified by gauge mediated supersymmetry breaking
(see above).
 SM error bars for the muon anomalous magnetic moment have been
updated. Still: the SM hadronic contribution to (g2)_muon is estimated
from e+e data only, data from tau decays is ignored; consequently the
exp. BNL result differs from the SM by ~3 sigma. If the line "CALL
MAGNMU(PAR,PROB)" is present in the main routines, it is checked
whether a point in parameter space differs by at most by 2 sigma from
the exp. result. Hence, points in parameter space for which the Susy
contribuion is too small (or vanishing or of the wrong sign) generate
warnings or are considered as phenomenologically excluded. If this
constraint is considered as too strong by the user it is enough to
comment the line "CALL MAGNMU(PAR,PROB)" in the main program.
 Exp. constraints from the Tevatron on sbottoms updated.
17.03.2008: Release of version 1.3.0
 The muon anomalous magnetic moment is computed (to two loops) according
to F. Domingo and U. Ellwanger, JHEP 0712:090 (2007),
arXiv:0710.3714 [hepph]; a
warning is issued if the susy contribution incl. theoret. errors is
less than ~10^{9} or larger than ~5*10^{9}
 The input of the soft parameters of the first two generations is
simplified for the versions nmhdecay
 The versions nmspec_grid and _rand allow to scan over nonuniversal
soft parameters at the GUT scale
 A script "run" is added which allows to run userdefined input
files
25.02.2008: version 1.2.5, 2 loop RGE for A_muon corrected
20.02.2008: version 1.2.4
 Bug in electroweak corrs. in bsg.f corrected (effect essentially
on B>mumu) > the
Benchmark point 3 in BMP.tgz above has been updated
 The Smuon masses (and the trilinear coupling A_muon) are treated
separately from the Selectron masses
 Small bug in the error management (if PROB(I) < 0) corrected
08.02.2008: version 1.2.3
 Bug in msferm corrected (slepton masses of the two first generations)
 Small bug in mhiggs corrected (bosonic loop corrs.; effect <~1 MeV on Mh_1)
 Floats rewritten as "D0"; Fortran compiler g77 used if available
05.02.2008: version 1.2.2
 Treatment of floats in bsg.f improved, consistency of
the output format with SLHA2 improved
 Benchmark points from the paper "Benchmark scenarios for the NMSSM",
A. Djouadi et al., arXiv:0801.4321 [hepph],
made available (see Downloads above)
03.12.2007: version 1.2.1, small bug in the RGE for lambda
corrected
30.10.2007: Release of version 1.2.0:
 inclusion of constraints from B physics according to F. Domingo and U.
Ellwanger, JHEP 0712:090 (2007), arXiv:0710.3714 [hepph]
 nonuniversal boundary conditions at the GUT scale in NMSPEC are made
possible
 improved convergence of NMSPEC for large lambda, improved
numerical agreement between NMHDECAY and NMSPEC
 the scale of mu_eff in the subroutine msferm.f corrected (with thanks
to Junjie Cao)
12.07.2007: version 1.1.1:
 improved convergence of nmspec near focus point regions
 bug in MAdependent threshold corrections for g1 and g2 in rges.f corrected
(affects only the calculation of MGUT)
11.06.2007: Bug in the large tan(beta) corrections to the
branching ratios (in decay.f) corrected (with thanks to I. Rottlaender)
30.03.2007: Bug in the error management in nmhdecay (in case of
negative Higgs masses squared) corrected
27.03.2007: Release of version 1.1.0:
 the MSSMlike parameter MA at M_SUSY can be used as input instead of
A_lambda (within nmhdecay, nmhdecay_scan or nmhdecay_grid)
 improved diagonalisation of Higgs mass matrices (avoids problems in
cases with huge mass differences)
 bug in the calculation of the Higgs wave function renormalisation
constants Z in the case MA >> M_susy corrected
 if soft Higgs masses  as computed internally in nmhdecay  are much
larger than the SUSY scale, a warning is indicated
 the micrOMEGAs version included in the NMSSMTools package has been
updated:
 bugs in the calculation of the Higgs widths in micrOMEGAs
corrected
 if the Higgs widths computed by micrOMEGAs are still unreliable
(which could affect the relic density if the LSP annihilates through a
Higgs resonance) a warning is indicated in the output file
31.01.2007: version 1.0.2: bug in the Higgs loop corrections to
M_higgs corrected,
08.01.2007: version 1.0.1: bug in the Susy threshold corrections in the running
of h_top and h_bottom, and the electroweak two loop corrections
corrected (with thanks to P. Slavich)
12.12.2006: Release of NMSSMTools, version 1.0.0