SOFTSUSY: by Ben Allanach

This code provides a SUSY spectrum in the MSSM consistent with input low energy data, and a user supplied high energy constraint (eg minmal SUGRA). It is written in C++ with an emphasis on easy generalisability. Full three-family couplings and renormalisation group equations are employed, as well as one-loop finite corrections a la Bagger, Matchev, Pierce and Zhang. It can be used in conjunction with other programs for many different particle physics calculations: see a SUSY tools review.
SOFTSUSY has been used to generate points in the CMSSM that fit current indirect data: you can obtain such points with their properties at the KISMET page.
If you use SOFTSUSY to write a paper, please cite

B.C. Allanach, Comput. Phys. Commun. 143 (2002) 305-331, hep-ph/0104145,
which is the SOFTSUSY manual. The version on the electronic hep-ph/ archive will be updated with more recent versions (since 1.2).
You can run SOFTSUSY for one parameter point checking it against other codes from this web-page. There are detailed comparisons between SOFTSUSY and other publicly available codes in

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

Installation

To run SOFTSUSY, you should only need standard C++ libraries.

The following releases contain a test program (main.cpp) and the code for the SOFTSUSY library (libsoft.a, link with -lsoft). There is also a Makefile: in linux, just unpack the files with (eg for softsusy2.0)
> gunzip softsusy2.0.tar.gz
> tar -xvf softsusy2.0.tar
then type
> make
The C++ test program is then run by
> ./softpoint.x

> ./softsusy.x

Release 2.0.18 21/05/08

BUG-fixes: Gzipped tarred [source 2.0.18 | manual | documentation]: bug-fixed Z radiative correction to neutralino mass. Per mille differences resulting. Thanks to R Ruiz.
05/03/08: [source 2.0.17] More bugs in form of charged Higgs loops introduced by recent "bug fixes" - it should all be fixed in this version.
03/03/08: [source 2.0.16] More bugs in sign of charged Higgs loops in Higgs masses introduced in 2.0.15: thanks to Roberto Ruiz. 2.0.15 could have Higgs masses out by 3-4 GeV for heavy points and other sfermions out by a percent or so.
11/02/08: [source 2.0.15] Bug-fixed sign of charged Higgs corrections, introduced in v2.0.1. Could change slepton masses by several percent in some regions of parameter space. Thanks to J Kersten and W Porod.
09/08/07: [source 2.0.14] Bug-fixed tan beta and Higgs VEV RGEs: added Feynman gauge dependent pieces. Numerically, this makes a negligible difference (smaller than 10^-5 fractional error in all masses/mixing angles in MSUGRA). Thanks to A Shelpnyakov
09/07/07 [source 2.0.13] now works with negative m0 values in msugra - request from A Pukhov
28/06/07 [source 2.0.12] minor bug-fix in slha output file - spaces added so that SUSYHIT can read it. Thanks to S Sekmen
13/02/07 [source 2.0.11] small bug-fix in higgs mass squared parameters in the EXTPAR option input of SLHA files. Thanks to Steve Martin.
23/01/07 [source 2.0.10] Bug-fixed D0 Passarino-Veltman function - special cases now handled correctly (these are very unlikely to affect your results). Thanks to Shandou Tetsuo.
15/01/07 [source 2.0.9] bug-fixed loop corrections to neutralino mass. Negligible change to neutralino_1 mass, but heavier neutralino masses change by a percent or so at SPS1a. Thanks to Steve Kom.
07/09/06 [source 2.0.8] Bug-fixed 1-loop threshold correction to sneutrino mass. Resulting pole mass changes very slightly. Thanks to J Kersten.
22/06/06 [source 2.0.7] bug-fixed SLHA implentation of GMSB in softpoint.cpp. No change to rest of calculation. Thanks to Peter Skands.
25/04/06 [source 2.0.6] at high masses, there were sometimes convergence problems associated with numerical rounding errors in the p=0 limit of Passarino-Veltman functions. This has now been fixed by adding more tolerance before going to p=0.
20/02/06 [source 2.0.5 ] small error in two-loop term of right-handed charged slepton mass squared RGE fixed. Thanks to R Ruiz.
17/01/06 [ source 2.0.4 ]: mistake introduced in 2.0.3 in charged Higgs loop corrections to neutral higgs masses fixed. Thanks to R Ruiz for flagging the error.
14/11/05 [ source 2.0.3 ]: sbottom/stop mass ordering debugged:
program would sometimes crash due to inconsistent mass ordering. No changes to calculation, except of fine-tuning. ht fine-tuning now not included by default.
27/10/05 [ source 2.0.2 ]: grid of Q-values now conforms to SLHA definition. No changes to calculation.
04/10/05 [ source 2.0.1 ]: sfermions' threshold corrections debugged.
26/08/05 [ source 2.0 ]
Upgraded accuracy: full 2-loop MSSM RGEs default, full 1-loop MSSM threshold corrections.
Simpler user interface: massIn has been done away with. To change Standard Model inputs, use either a main program like main.cpp for a scan, or the SUSY les Houches Accord Input format. softpoint.x otherwise uses default Standard Model inputs. Deleted confusing general boundary conditions option in SOFTSUSY format - use the SLHA.
SUSY Les Houches Accord updated and corrected. Isajet7.64 interface debugged.
mgluino now calculated at Q=MSUSY
number of loops in scalar RGEs promoted to a changeable global variable

Executable files: after installation

softpoint.x: command-line interface. GMSB, AMSB, mSUGRA and general boundary conditions possible.
softsusy.x: example C++ test program - loops over tan beta, printing out Higgs masses around SPS1a.

Source Files

Contained in the gzipped tarred archive above, but can be obtained seperately here:

Input and info files

[README] contains installation instructions and a list of included files
[outputTest] is the output from the test program in main.cpp
[lesHouchesInput] is an alternative input file in the SUSY Les Houches Accord format
[lesHouchesOutput] is the result of running with the above input file, ie
softpoint.x leshouches < lesHouchesInput > lesHouchesOutput

Code

[def.h] contains switches (options) and parameters such as default fermion masses, required accuracy etc
[main.cpp | main.h] contains a main C++ program to calculate Higgs masses as a function of tan beta
[softpoint.cpp | softpoint.h] is the main program for the command line interface
[mycomplex.cpp | mycomplex.h] contains Complex class - standard complex algebra
[linalg.cpp | linalg.h] contains Matrix and Vector classes - standard linear algebra and diagonalisation routines
[lowe.cpp | lowe.h] contains QedQcd class - contains fermion mass data, as well as QED x QCD gauge couplings. Also runs them to 1 loop x 3 loops.
[numerics.cpp | numerics.h] contains numerical routines - differential equation solver, differentiator and function minimiser
[rge.cpp | rge.h] a general RGE class - describes how to run an RGE object, and how many parameters the class has (and its scale) and how to access the beta functions etc
[softsusy.cpp | softsusy.h] the main part of the code. A MssmSoftsusy object contains the low energy data. Also contained are the iterating routines to do the main calculation and radiative corrections etc.
[softpars.cpp | softpars.h]. A SoftParsMssm objecy contains all soft breaking and SUSY (and physical) parameters,
[physpars.cpp | physpars.h]. Contains flags for potential problems in sProblem structure, and structure for containing physical MSSM parameters such as pole masses in sPhysical.
[susy.cpp | susy.h] MssmSusy contains all SUSY couplings and tan beta, as well as their beta functions
[utils.cpp | utils.h] contains a few handy little mathematical functions
[twoloophiggs.f | twoloophiggs.h] contains the Slavich-et-al two loop Higgs mass/REWSB corrections
[Makefile]

Extensions

Non-standard supersymmetry breaking

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

Ben Allanach's homepage:

Ben Allanach B2.04, DAMTP, CMS, University of Cambridge Wilberforce Road, Cambridge CB3 0WA, UK. Tel: +44 (0)1223 766889