SPheno, a program for calculating supersymmetric spectra, SUSY particle decays and SUSY particle production at e+e- colliders. W. Porod.

PROGRAM SUMMARY
Title of program: SPheno
Catalogue identifier: ADRV
Ref. in CPC: 153(2003)275
Distribution format: tar gzip file
Operating system: Linux, Unix
Number of lines in distributed program, including test data, etc: 59999
Keywords: Supersymmetric standard model, Renormalization group equations, MSSM spectra, Runge-Kutta, Decays of super- symmetric particles, Production, Elementary particle physics, Phenomenological model.
Programming language used: Fortran
Computer: PC , HP machine , ALPHA machine .

Nature of physical problem:
The first issue is the determination of the masses and couplings of supersymmetric particles in the R-parity conserved MSSM. Low energy data on Standard Model fermion masses, gauge couplings and electroweak gauge boson masses serve as constraints. Radiative corrections from supersymmetric particles to these inputs must be calculated. Theoretical constraints on the soft SUSY breaking parameters from a high scale theory are imposed and the parameters at the electroweak scale are obtained from the high scale parameters by evaluating the corresponding renormalization group equations. These parameters must be consistent with the requirements of correct electroweak symmetry breaking. The second issue is to use the obtained masses and couplings for calculating decay widths and branching ratios of supersymmetric particles as well as the cross sections for these particles in electron positron annihilation. The third issue is to calculate the following low energy constraints: the branching ratio for b ->sgamma, the SUSY contributions to the rho parameter as well as the SUSY contributions to the anomalous magnetic moment of the muon.

Method of solution:
The renormalization connecting a high scale and the electroweak scale is calculated by the Runge-Kutta method. Iteration provides a solution consistent with the multi-boundary conditions. In case of three-body decays and for the calculation of initial state radiation Gaussian quadrature is used for the numerical solution of the integrals.

Restrictions:
MSSM spectra with real parameters only neglecting the effect of mixing between (s)fermion generations.

Typical running time:
0.3 s on a Pentium III.