SOFTSUSY
by: B C Allanach, P Athron, A Bednyakov, M Bernhardt, T Cridge, D Grellscheid, M Hanussek, C H Kom, S Martin, D Robertson, R Ruiz de Austri, P Slavich, L Tunstall, A Voigt and A G Williams
manuals [main | decays | NMSSM | RPV | RPV nu | 2 loop gauge/Yukawa | 2 loop squark/gluino]
Summary
This program provides a SUSY spectrum in the NMSSM, or the MSSM including flavour violation and with or without R-parity consistent with input Standard Model fermion mass/mixings and electroweak/strong coupling data. The R-parity violating mode can calculate neutrino masses and mixings to 1 loop. SOFTSUSY can be used in conjunction with other programs for many different particle physics calculations: see a SUSY tools review. See graphs of [ gauge unification | (whole range) | MSSM RGE evolution ] (feel free to use them). SOFTSUSY now has a mode with 3 loop RGEs and some 2-loop threshold correction and 2-loop SUSY QCD corrections to gluino and squark pole masses. SOFTSUSY now computes decay branching ratios for the MSSM and NMSSM.
Latest Release
If you want to receive an email every time there's a new version, send an email to Ben Allanach, asking to be added to the mailing list.- source 4.0.3 (27/07/17)
- Removed GiNaC and CLN dependencies when two-loop corrections to gauge and Yukawa coupling threshold corrections are switched on (see High Accuracy Mode below): output same to 1.0e-5, but much faster running: BCA, A Bednyakov and A Voigt
- source 4.0.2 (12/07/17)
- higgs->Z gamma corrected in NMSSM case - it was missing a contribution
- higgs->mu mu decays implemented with pole muon mass when the Yukawa coupling is zero
- H->WW* H->ZZ* decays have prefactors altered to hopefully include more higher order effects
- Fixed problem that wasn't flagging the non-convergence of the fixed point iteration algorithm: thanks to A Ibrahim for finding it
- Made 0 mixing default (ie no quark mixing, but non-zero masses for second/third family)
- Fixed erroneous muon Yukawa in delta_vB calculation: less than per mille effects. Thanks to M Back and A Voigt.
- Fixed 2-loop mt threshold corrections: boils down to using running mt in definition of stop mixing in one-loop corrections. Only relevant when --enable-two-loop-gauge-yukawa-compilation is switched on.
- Major upgrade: added MSSM/NMSSM decays
References
These reference are collected here in a bibTeX file for handy inclusion into latex documents.
If you use SOFTSUSY to write a paper, please cite (see MCnet guidelines)
which is the SOFTSUSY manual for the R-parity conserving MSSM.
If you use the decay calculations, please cite [Allanach:2001kg] and
If you calculate in the NMSSM, please cite [Allanach:2001kg] and
If you use the R-parity violating aspects, please cite [Allanach:2001kg] and
If you use it to calculate neutrino masses and mixings, please cite [Allanach:2001kg], [Allanach:2009bv] and
If you use the three-loop RGEs or two-loop threshold corrections to gauge/Yukawa couplings, please cite [Allanach:2001kg] and
If you use the two-loop SUSY QCD corrections to squark and gluino pole masses, please cite [Allanach:2001kg] and
Installation
For ultra-basic instructions, see the introduction video. Otherwise, see a quick tutorial given at BUSSTEPP 2012.To run SOFTSUSY, you should need only standard C++ and fortran libraries.
The following releases contain a test program (main.cpp) and the
SOFTSUSY library (libsoft.a, link with -L.libs -lsoft). In linux, just unpack the files with (eg for softsusy-4.0)
> gunzip softsusy-4.0.tar.gz > tar -xvf softsusy-4.0.tar > cd softsusy-4.0
> ./configure > make
> ./softsusy.x > inOutFiles/outputTest > ./rpvsoftsusy.x > inOutFiles/rpvOutputTest > ./rpvneut.x > inOutFiles/neutOutputTest > ./softsusy-nmssm.x > inOutFiles/outputTest-nmssm > ./higher.x > ./softpoint.x
> ./softpoint.x leshouches < inOutFiles/lesHouchesInput > inOutFiles/lesHouchesOutput > ./softpoint.x leshouches < inOutFiles/nmssmSLHAnoZ3Input > inOutFiles/nmssmSLHAnoZ3Output > ./softpoint.x leshouches < inOutFiles/nmssmSLHAZ3Input > inOutFiles/nmssmSLHAZ3Output > ./softpoint.x leshouches < inOutFiles/rpvHouchesInput > inOutFiles/rpvHouchesOutput > ./softpoint.x leshouches < inOutFiles/slha2Input > inOutFiles/slha2Output
SOFTSUSY executables use no input or output files except from standard input or standard output.
Note that the executables are actually wrapper scripts, the "true" executables lie in the directory .libs/ (from the directory where you installed softsusy).
If you have trouble with compilation, try this page.
Particle Decays
An example point including the calculation of sparticle decays, neglecting modes with a branching ratio of less than 1.0e-5, and outputting the partial widths in the comments:./softpoint.x gmsb --n5=2 --mMess=1.0e6 --LAMBDA=5.0e5 --tanBeta=10 --sgnMu=1 --decays --minBR=1.0e-5 --outputPartialWidths
High Orders Mode
If the 2-loop SUSY QCD corrections to squark and gluino masses are required, do> ./configure --enable-two-loop-sparticle-mass-compilation > make
> ./softpoint.x sugra --tol=1.0e-4 --m0=1000 --m12=1000 --a0=0 --tanBeta=10 --sgnMu=1 --two-loop-sparticle-masses --two-loop-sparticle-mass-method=1
High Accuracy Mode
If the high accuracy mode with 3-loop RGEs and some 2-loop threshold corrections is required, use the configure options:> ./configure --enable-two-loop-gauge-yukawa-compilation --enable-three-loop-rge-compilation > make
> ./softpoint.x sugra --tol=1.0e-5 --m0=7240 --m12=800 --a0=-6000 --tanBeta=50 --sgnMu=1 --mt=173.2 --alpha_s=0.1187 --mbmb=4.18 --two-loop-gauge-yukawa --three-loop-rges
Executable files: after installation
- softpoint.x: command-line interface. SLHA2, NMSSM, RPV, GMSB, AMSB, mSUGRA and general boundary conditions possible. Main program: src/softpoint.cpp
- softsusy.x: example C++ test program - loops over tan beta, printing out Higgs masses around CMSSM10.1. Main program: src/main.cpp
- softsusy-nmssm.x: example NMSSM test program - loops over tan beta. Main program: src/main-nmssm.cpp
- rpvsoftsusy.x: example R-parity violating test program - loops over l''_{323} for CMSSM10.1, printing out the RH stop mass. Main program: src/rpvmain.cpp
- rpvneut.x: example neutrino mass calculating R-parity violating test program. Main program: src/rpvneut.cpp
- higher.x: example program to illustrate 2-loop QCD corrections to squark and gluino pole masses (note - you must install NLLFAST) to compile this and use ./configure --enable-two-loop-sparticle-mass-compilation before the make command.
SOFTSUSY-specific input for SUSY Les Houches Accord input files
Block SOFTSUSY # SOFTSUSY specific inputs
0 <calcDecays> # Set to 1 to calculate decays/0 (default) for no decays
1 <TOLERANCE> # desired fractional accuracy in output
2 <MIXING> # quark mixing option
3 <PRINTOUT> # gives additional verbose output during calculation
4 <QEWSB> # change electroweak symmetry breaking scale
5 <INCLUDE_2_LOOP_SCALAR_CORRECTIONS> # Full 2-loop running in RGEs
6 <PRECISION> # number of significant figures in SLHA output
7 <numHiggsLoops> # number of loops in REWSB/mh calculation
8 <susyRpvBCatMSUSY> # Switch MSUSY-scale RPV boundary conditions ON
9 <invertedOutput> # RPV neutrino output uses normal hierarchy (=0.0) or inverted (=1.0)
10 <forceSlha1> # if =1, tries to force output into SLHA *1* format
11 <m32> # sets gravitino mass to m32
12 <printSpectra> # Prints spectrum even when point is theoretically excluded if=1
13 <mAFlag> # If=0 (default), sets tachyonic mA=0, otherwise mA=sqrt(|mA|^2)
15 <NMSSMTools> # If=1, enables NMSSMTools compatible SLHA2 output
16 <MICROMEGAS> # Micromegas options for NMSSMTools use: 1=RD, 2=DD, 3=ID, 4=both
17 <NMSDECAY> # If=1, flags for sparticle decays to be calculated via NMSDECAY
18 <SoftHiggsOut> # If=1, then the EWSB conditions output soft Higgs masses in NMSSM
# 19 Only works with ./configure --enable-three-loop-rge-compilation
19 <threeLoopRGEs> # If=1, then 3-loop MSSM RGEs included (default of 0 to disable)
# 20 Only works with ./configure --enable-two-loop-gauge-yukawa-compilation
20 <gyThresholds> # If>0, switch on gauge/Yukawa two-loop thresholds (see manual [5] for details).
If=31, they all are switched on (default 0 to disable).
# 22,23 only work with --enable-two-loop-sparticle-mass-compilation
22 <2-loop squark/gluino> # Include 2-loop terms in gluino/squark masses (1 for on, 0 for off)
23 <expandAroundGluinoPole> # sets expandAroundGluinoPole parameter (default 3)
24 <minBR> # If decay BR is below this number, don't output
25 <threeBodyDecays> # If set to 0, don't calculate 3-body decays (1=default)
26 <outputPartialWidths> # If set to 1, output partial widths (0=default)
Source Files
Contained in the gzipped tarred archive above, but the details of some of them are listed here:Input, output and information files
These input and output files are contained in the inOutFiles/ subdirectory:- [slha2Input] is an alternative input file in the SUSY Les Houches Accord 2 format for CMSSM10.1
- [slha2Output] is the result of running with the above input file and includes flavour violation, for inclusion into codes like SusyBsg which include flavour corrections
- [outputTest] is the output from the test program in main.cpp
- [rpvOutputTest] is the output from the test program in rpvmain.cpp
- [lesHouchesInput] is an input file in the SUSY Les Houches Accord format for CMSSM10.1
- [lesHouchesOutput] is the result of
running with the above input file, ie
- [rpvHouchesInput] is an alternative input file in the SUSY Les Houches Accord 2 format
- [rpvHouchesOutput] is the result of
running with the above input file, ie
softpoint.x leshouches < lesHouchesInput > lesHouchesOutput
Code
The source code is contained in the src/ subdirectory.- [def.h] contains switches (options) and parameters such as default gauge boson masses, required accuracy etc
- List and description of files
- List and description of classes
Comparisons With Other SUSY Spectrum Generators
There are detailed comparisons between SOFTSUSY and other R-parity conserving MSSM publicly available programs in- B.C. Allanach and T. Cridge, arxiv:1703.09717
- Precise Determination of the Neutral Higgs Boson Masses in the MSSM, B.C. Allanach, A. Djouadi, J.L. Kneur, W. Porod, P. Slavich, JHEP 0409 (2004) 044, hep-ph/0406166,
- Theoretical uncertainties in sparticle mass predictions from computational tools, B.C. Allanach, S. Kraml, W. Porod, JHEP 03 (2003) 045, hep-ph/0302102
- B.C. Allanach and T. Cridge, arxiv:1703.09717
- Higgs mass predictions of public NMSSM spectrum generators, Staub et al, Comp. Phys. Comm. 202 (2016) 113, arXiv:1507.05093
- Improved predictions for intermediate and heavy Supersymmetry in the MSSM and beyond Staub and Porod, arXiv:1703.03267
License
SOFTSUSY Copyright (C) 2007 B.C. Allanach
This program is free software: you can redistribute it and/or modify
it under the terms of the GNU General Public License as published by
the Free Software Foundation, either version 3 of the License, or
(at your option) any later version.
This program is distributed in the hope that it will be useful,
but WITHOUT ANY WARRANTY; without even the implied warranty of
MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
GNU General Public License for more details.
See http://www.gnu.org/licenses/.