micrOMEGAs: Version 1.3

We present the latest version of micrOMEGAs, a code that calculates the relic density of the lightest supersymmetric particle (LSP) in the minimal supersymmetric standard model (MSSM). All tree-level processes for the annihilation of the LSP are included as well as all possible coannihilation processes with neutralinos, charginos, sleptons, squarks and gluinos. The cross-sections extracted from CalcHEP are calculated exactly using loop-corrected masses and mixings as specified in the SUSY Les Houches Accord. Relativistic formulae for the thermal average are used and care is taken to handle poles and thresholds by adopting specific integration routines. The input parameters can be either the soft SUSY parameters in a general MSSM or the parameters of a SUGRA model specified at the GUT scale. In the latter case, a link with Suspect, SOFTSUSY, Spheno and Isajet allows one to calculate the supersymmetric spectrum, Higgs masses, as well as mixing matrices. Higher-order corrections to Higgs couplings to quark pairs including QCD as well as some SUSY corrections (Δmb) are implemented. Routines calculating μ(g−2), b→sγ and Bs→μ⁺μ⁻ are also included. In particular the b→sγ routine includes an improved NLO for the SM and the charged Higgs while the SUSY large tanβ effects beyond leading-order are included. This new version also provides cross-sections for any 2→2 process as well as partial decay widths for two-body final states in the MSSM allowing for easy simulation at colliders.

Program summary

Program title:micrOMEGAs1.3Catalogue identifier:ADQR_v1_3Program summary URL:http://cpc.cs.qub.ac.uk/summaries/ADQR_v1_3Program obtainable from: CPC Program Library, Queen's University of Belfast, N. IrelandLicensing provisions:noneComputer:PC, Alpha, Silicon graphics, SunProgramming language:C and FortranOperating system:UNIX (Linux, OSF1, IRIX64, SunOS)RAM:20 MB depending on the number of processes requiredNo of lines in distributed program, including test data, etc.:78 314No. of bytes in distributed program, including test data, etc.:703 112Distribution format:tar.gzNumber of processors used:1External routines/libraries:Library of Fortran functions, for example, -lg2c (platform dependent)Catalogue identifier of previous version:ADQRJournal reference of previous version:Comput. Phys. Comm. 149 (2002) 103Does the new version supersede the previous version?:yesNature of problem:Calculation of the relic density of the lightest supersymmetric particle in the MSSM.Solution method:In numerically solving the evolution equation for the density of dark matter, relativistic formulae for the thermal average are used. All tree-level processes for annihilation and coannihilation of SUSY particles are included. The cross-sections for all processes are calculated exactly with CalcHEP. Higher-order corrections to Higgs masses and Higgs couplings to quark pairs including QCD as well as some SUSY corrections are implemented. The input parameters can be either the soft SUSY parameters in a general MSSM or the parameters of a SUGRA model specified at the GUT scale. In the latter case, a link with Suspect, SOFTSUSY, Spheno and Isajet allows to calculate the supersymmetric spectrum, Higgs masses, as well as mixing matrices.Reasons for the new version:This new version contains a more accurate calculation of the relic density of dark matter as well as many new features both for interface with codes that calculate the supersymmetric spectrum as well as for computation of cross-sections and decays relevant for collider physics.Summary of revisions:

•

Interface with the main codes to calculate the supersymmetric spectrum: Suspect, Isajet, Spheno and SOFTSUSY in models defined at some high scale.

•

Includes loop corrected sparticle masses and mixing matrices.

•

Includes loop-corrected Higgs masses and widths. QCD corrections to the Higgs couplings to fermion pairs are included as well as, via an effective Lagrangian, the Δmb correction relevant at large tanβ.

•

Provides exact numerical solution of the Boltzmann equation by Runge-Kutta.

•

Outputs the relative contribution of each channel to 1/Ω.

•

Computes cross-sections for any 2→2 process at the parton level.

•

Calculates decay widths for all particles at tree-level including all 1→2 decay modes.

•

Calculates constraints on MSSM: Bs→μ⁺μ⁻ and NLO corrections to b→sγ.

Unusual features:Depending on the parameters of the model, the program generates additional new code for matrix elements, compiles it and loads it dynamically.Running time:0.2 sec.     

CPsuperH2.0: An improved computational tool for Higgs phenomenology in the MSSM with explicit CP violation

We describe the Fortran code CPsuperH2.0, which contains several improvements and extensions of its predecessor CPsuperH. It implements improved calculations of the Higgs-boson pole masses, notably a full treatment of the 4×4 neutral Higgs propagator matrix including the Goldstone boson and a more complete treatment of threshold effects in self-energies and Yukawa couplings, improved treatments of two-body Higgs decays, some important three-body decays, and two-loop Higgs-mediated contributions to electric dipole moments. CPsuperH2.0 also implements an integrated treatment of several B-meson observables, including the branching ratios of Bs→μ⁺μ⁻, Bd→τ⁺τ⁻, Bu→τν, B→Xsγ and the latter's CP-violating asymmetry ACP, and the supersymmetric contributions to the mass differences. These additions make CPsuperH2.0 an attractive integrated tool for analyzing supersymmetric CP and flavour physics as well as searches for new physics at high-energy colliders such as the Tevatron, LHC and linear colliders.¹

Program summary

Program title: CPsuperH2.0Catalogue identifier: ADSR_v2_0Program summary URL:http://cpc.cs.qub.ac.uk/summaries/ADSR_v2_0.htmlProgram obtainable from: CPC Program Library, Queen's University, Belfast, N. IrelandLicensing provisions: Standard CPC licence, http://cpc.cs.qub.ac.uk/licence/licence.htmlNo. of lines in distributed program, including test data, etc.: 13 290No. of bytes in distributed program, including test data, etc.: 89 540Distribution format: tar.gzProgramming language: Fortran 77Computer: PC running under Linux and computers in Unix environmentOperating system: LinuxRAM: 32 MbytesClassification: 11.1Catalogue identifier of the previous version: ADSR_v1_0Journal reference of the previous version: CPC 156 (2004) 283Does the new version supersede the previous version?: YesNature of problem: The calculations of mass spectrum, decay widths and branching ratios of the neutral and charged Higgs bosons in the Minimal Supersymmetric Standard Model with explicit CP violation have been improved. The program is based on recent renormalization-group-improved diagrammatic calculations that include dominant higher-order logarithmic and threshold corrections, b-quark Yukawa-coupling resummation effects and improved treatment of Higgs-boson pole-mass shifts. The couplings of the Higgs bosons to the Standard Model gauge bosons and fermions, to their supersymmetric partners and all the trilinear and quartic Higgs-boson self-couplings are also calculated. The new implementations include a full treatment of the 4×4(2×2) neutral (charged) Higgs propagator matrix together with the center-of-mass dependent Higgs-boson couplings to gluons and photons, two-loop Higgs-mediated contributions to electric dipole moments, and an integrated treatment of several B-meson observables.Solution method: One-dimensional numerical integration for several Higgs-decay modes, iterative treatment of the threshold corrections and Higgs-boson pole masses, and the numerical diagonalization of the neutralino mass matrix.Reasons for new version: Mainly to provide a coherent numerical framework which calculates consistently observables for both low- and high-energy experiments.Summary of revisions: Improved treatment of Higgs-boson masses and propagators. Improved treatment of Higgs-boson couplings and decays. Higgs-mediated two-loop electric dipole moments. B-meson observables.Running time: Less than 0.1 seconds.     

SuSpect: A Fortran code for the Supersymmetric and Higgs particle spectrum in the MSSM

We present the Fortran code SuSpect version 2.3, which calculates the Supersymmetric and Higgs particle spectrum in the Minimal Supersymmetric Standard Model (MSSM). The calculation can be performed in constrained models with universal boundary conditions at high scales such as the gravity (mSUGRA), anomaly (AMSB) or gauge (GMSB) mediated supersymmetry breaking models, but also in the non-universal MSSM case with R-parity and CP conservation. Care has been taken to treat important features such as the renormalization group evolution of parameters between low and high energy scales, the consistent implementation of radiative electroweak symmetry breaking and the calculation of the physical masses of the Higgs bosons and supersymmetric particles taking into account the dominant radiative corrections. Some checks of important theoretical and experimental features, such as the absence of non-desired minima, large fine-tuning in the electroweak symmetry breaking condition, as well as agreement with precision measurements can be performed. The program is simple to use, self-contained and can easily be linked to other codes; it is rather fast and flexible, thus allowing scans of the parameter space with several possible options and choices for model assumptions and approximations.

Program summary

Title of program:SuSpectCatalogue identifier:ADYR_v1_0Program summary URL:http://cpc.cs.qub.ac.uk/summaries/ADYR_v1_0Program obtainable from: CPC Program Library, Queen's University of Belfast, N. IrelandLicensing provisions:noneProgramming language used:FORTRAN 77Computer:Unix machines, PCNo. of lines in distributed program, including test data, etc.:21 821No. of bytes in distributed program, including test data, etc.:249 657Distribution format:tar.gzOperating system:Unix (or Linux)RAM:approximately 2500 KbytesNumber of processors used:1 processorNature of problem:SuSpect calculates the supersymmetric and Higgs particle spectrum (masses and some other relevant parameters) in the unconstrained Minimal Supersymmetric Standard Model (MSSM), as well as in constrained models (cMSSMs) such as the minimal Supergravity (mSUGRA), the gauge mediated (GMSB) and anomaly mediated (AMSB) Supersymmetry breaking scenarii. The following features and ingredients are included: renormalization group evolution between low and high energy scales, consistent implementation of radiative electroweak symmetry breaking, calculation of the physical particle masses with radiative corrections at the one- and two-loop level.Solution method:The main methods used in the code are: (1) an (adaptative fourth-order) Runge-Kutta type algorithm (following a standard algorithm described in “Numerical Recipes”), used to solve numerically a set of coupled differential equations resulting from the renormalization group equations at the two-loop level of the perturbative expansions; (2) diagonalizations of mass matrices; (3) some mathematical (Spence, etc) functions resulting from the evaluation of one and two-loop integrals using the Feynman graphs techniques for radiative corrections to the particle masses; (4) finally, some fixed-point iterative algorithms to solve non-linear equations for some of the relevant output parameters.Restrictions:(1) The code is limited at the moment to real input parameters. (2) It also does not include flavor non-diagonal terms which are possible in the most general soft supersymmetry breaking Lagrangian. (3) There are some (mild) limitations on the possible range of values of input parameter, i.e. not any arbitrary values of some input parameters are allowed: these limitations are essentially based on physical rather than algorithmic issues, and warning flags and other protections are installed to avoid as much as possible execution failure if unappropriate input values are used.Running time:between 1 and 3 seconds depending on options, with a 1 GHz processor.     

HFOLD - A program package for calculating two-body MSSM Higgs decays at full one-loop level

HFOLD (Higgs Full One Loop Decays) is a Fortran program package for calculating all MSSM Higgs two-body decay widths and the corresponding branching ratios at full one-loop level. The package is done in the SUSY Parameter Analysis convention and supports the SUSY Les Houches Accord input and output format.

Program summary

Program title: HFOLDCatalogue identifier: AEJG_v1_0Program summary URL:http://cpc.cs.qub.ac.uk/summaries/AEJG_v1_0.htmlProgram obtainable from: CPC Program Library, Queen?s University, Belfast, N. IrelandLicensing provisions: Standard CPC licence, http://cpc.cs.qub.ac.uk/licence/licence.htmlNo. of lines in distributed program, including test data, etc.: 340 621No. of bytes in distributed program, including test data, etc.: 1 760 051Distribution format: tar.gzProgramming language: Fortran 77Computer: Workstation, PCOperating system: LinuxRAM: 524 288 000 BytesClassification: 11.1External routines: LoopTools 2.2 (http://www.feynarts.de/looptools/), SLHALib 2.2 (http://www.feynarts.de/slha/). The LoopTools code is included in the distribution package.Nature of problem: A future high-energy e⁺e⁻ linear collider will be the best environment for the precise measurements of masses, cross sections, branching ratios, etc. Experimental accuracies are expected at the per-cent down to the per-mile level. These must be matched from the theoretical side. Therefore higher order calculations are mandatory.Solution method: This program package calculates all MSSM Higgs two-body decay widths and the corresponding branching ratios at full one-loop level. The renormalization is done in the DR scheme following the SUSY Parameter Analysis convention. The program supports the SUSY Les Houches Accord input and output format.Running time: The example provided takes only a few seconds to run.     

2HDMC - two-Higgs-doublet model calculator

We describe the public C++ code 2HDMC which can be used to perform calculations in a general, CP-conserving, two-Higgs-doublet model (2HDM). The program features simple conversion between different parametrizations of the 2HDM potential, a flexible Yukawa sector specification with choices of different Z₂-symmetries or more general couplings, a decay library including all two-body - and some three-body - decay modes for the Higgs bosons, and the possibility to calculate observables of interest for constraining the 2HDM parameter space, as well as theoretical constraints from positivity and unitarity. The latest version of the 2HDMC code and full documentation is available from: http://www.isv.uu.se/thep/MC/2HDMC.

Program summary

Program title:2HDMCCatalogue identifier: AEFI_v1_0Program summary URL:http://cpc.cs.qub.ac.uk/summaries/AEFI_v1_0.htmlProgram obtainable from: CPC Program Library, Queen's University, Belfast, N. IrelandLicensing provisions: GNU GPLNo. of lines in distributed program, including test data, etc.: 12 032No. of bytes in distributed program, including test data, etc.: 90 699Distribution format: tar.gzProgramming language: C++Computer: Any computer running LinuxOperating system: LinuxRAM: 5 MbClassification: 11.1External routines: GNU Scientific Library (http://www.gnu.org/software/gsl/)Nature of problem: Determining properties of the potential, calculation of mass spectrum, couplings, decay widths, oblique parameters, muon g−2, and collider constraints in a general two-Higgs-doublet model.Solution method: From arbitrary potential and Yukawa sector, tree-level relations are used to determine Higgs masses and couplings. Decay widths are calculated at leading order, including FCNC decays when applicable. Decays to off-shell vector bosons are obtained by numerical integration. Observables are computed (analytically or numerically) as function of the input parameters.Restrictions: CP-violation is not treated.Running time: Less than 0.1 s on a standard PC     

Decay of super-heavy particles: user guide of the SHdecay program

I give here a detailed user guide for the C++ program SHdecay, which has been developed for computing the final spectra of stable particles (protons, photons, LSPs, electrons, neutrinos of the three species and their antiparticles) arising from the decay of a super-heavy X particle. It allows to compute in great detail the complete decay cascade for any given decay mode into particles of the Minimal Supersymmetric Standard Model (MSSM). In particular, it takes into account all interactions of the MSSM during the perturbative cascade (including not only QCD, but also the electroweak and 3rd generation Yukawa interactions), and includes a detailed treatment of the SUSY decay cascade (for a given set of parameters) and of the non-perturbative hadronization process. All these features allow us to ensure energy conservation over the whole cascade up to a numerical accuracy of a few per mille. Yet, this program also allows to restrict the computation to QCD or SUSY-QCD frameworks. I detail the input and output files, describe the role of each part of the program, and include some advice for using it best.

Program summary

Title of program: SHdecayCatalogue identifier:ADSLProgram summary URL:http://cpc.cs.qub.ac.uk/summaries/ADSLProgram obtainable from: CPC Program Library, Queen's University of Belfast, N. IrelandComputer and operating system: Program tested on PC running Linux KDE and Suse 8.1Programming language used: C with STL C++ library and using the standard gnu g++ compilerNo. lines in distributed program: 14 955No. of bytes in distributed program, including test data, etc.: 624 487Distribution format: tar gzip fileKeywords: Super-heavy particles, fragmentation functions, DGLAP equations, supersymmetry, MSSM, UHECRNature of physical problem: Obtaining the energy spectra of the final stable decay products (protons, photons, electrons, the three species of neutrinos and the LSPs) of a decaying super-heavy X particle, within the framework of the Minimal Supersymmetric Standard Model (MSSM). It can be done numerically by solving the full set of DGLAP equations in the MSSM for the perturbative evolution of the fragmentation functions Dp2p1(x,Q) of any particle p1 into any other p2 (x is the energy fraction carried by the particle p2 and Q its virtuality), and by treating properly the different decay cascades of all unstable particles and the final hadronization of quarks and gluons. In order to obtain proper results at very low values of x (up to x∼10−13), NLO color coherence effects have been included by using the Modified Leading Log Approximation (MLLA).Method of solution: the DGLAP equations are solved by a four order Runge-Kutta method with a fixed step.Typical running time: Around 35 hours for the first run, but the most time consuming sub-programs can be run only once for most applications.     

SuperIso Relic: A program for calculating relic density and flavor physics observables in Supersymmetry

We describe SuperIso Relic, a public program for evaluation of relic density and flavor physics observables in the minimal supersymmetric extension of the Standard Model (MSSM). SuperIso Relic is an extension of the SuperIso program which adds to the flavor observables of SuperIso the computation of all possible annihilation and coannihilation processes of the LSP which are required for the relic density calculation. All amplitudes have been generated at the tree level with FeynArts/FormCalc, and widths of the Higgs bosons are computed with FeynHiggs at the two-loop level. SuperIso Relic also provides the possibility to modify the assumptions of the cosmological model, and to study their consequences on the relic density.

Program summary

Program title: SuperIso RelicCatalogue identifier: AEGD_v1_0Program summary URL:http://cpc.cs.qub.ac.uk/summaries/AEGD_v1_0.htmlProgram obtainable from: CPC Program Library, Queen's University, Belfast, N. IrelandLicensing provisions: yesNo. of lines in distributed program, including test data, etc.: 2 274 720No. of bytes in distributed program, including test data, etc.: 6 735 649Distribution format: tar.gzProgramming language: C (C99 Standard compliant) and FortranComputer: 32- or 64-bit PC, MacOperating system: Linux, MacOSRAM: 100 MbClassification: 1.9, 11.6External routines: ISASUGRA/ISAJET and/or SOFTSUSY, FeynHiggsDoes the new version supersede the previous version?: No (AEAN_v2_0)Nature of problem: Calculation of the lightest supersymmetric particle relic density, as well as flavor physics observables, in order to derive constraints on the supersymmetric parameter space.Solution method: SuperIso Relic uses a SUSY Les Houches Accord file, which can be either generated automatically via a call to SOFTSUSY or ISAJET, or provided by the user. This file contains the masses and couplings of the supersymmetric particles. SuperIso Relic then computes the lightest supersymmetric particle relic density as well as the most constraining flavor physics observables. To do so, it calculates first the widths of the Higgs bosons with FeynHiggs, and then it evaluates the squared amplitudes of the diagrams needed for the relic density calculation. These thousands of diagrams have been previously generated with the FeynArts/FormCalc package. SuperIso Relic is able to perform the calculations in different supersymmetry breaking scenarios, such as mSUGRA, NUHM, AMSB and GMSB.Reasons for new version: This version incorporates the calculation of the relic density, which is often used to constrain Supersymmetry.Summary of revisions:

•

Addition of the relic density calculation

•

Replacement of "float" type by "double".

Unusual features: SuperIso Relic includes the possibility of altering the underlying cosmological model and testing the influence of the cosmological assumptions.Additional comments: This program is closely associated with the "SuperIso" program - CPC Program Library, Catalogue Id. AEAN.Running time:Compilation time: a couple of hours for the statically linked version, a few minutes for the dynamically linked version. Running time: about 1 second, or a few seconds if libraries need to be compiled on the fly.     

MC-TESTER: a universal tool for comparisons of Monte Carlo predictions for particle decays in high energy physics

Theoretical predictions in high energy physics are routinely provided in the form of Monte Carlo generators. Comparisons of predictions from different programs and/or different initialization set-ups are often necessary. MC-TESTER can be used for such tests of decays of intermediate states (particles or resonances) in a semi-automated way. Our test consists of two steps. Different Monte Carlo programs are run; events with decays of a chosen particle are searched, decay trees are analyzed and appropriate information is stored. Then, at the analysis step, a list of all found decay modes is defined and branching ratios are calculated for both runs. Histograms of all scalar Lorentz-invariant masses constructed from the decay products are plotted and compared for each decay mode found in both runs. For each plot a measure of the difference of the distributions is calculated and its maximal value over all histograms for each decay channel is printed in a summary table. As an example of MC-TESTER application, we include a test with the τ lepton decay Monte Carlo generators, TAUOLA and PYTHIA. The HEPEVT (or LUJETS) common block is used as exclusive source of information on the generated events.

Program summary

Title of the program:MC-TESTER, version 1.1Catalogue identifier: ADSMProgram summary URL:http://cpc.cs.qub.ac.uk/summaries/ADSMProgram obtainable from: CPC Program Library, Queen's University of Belfast, N. IrelandComputer: PC, two Intel Xeon 2.0 GHz processors, 512MB RAMOperating system: Linux Red Hat 6.1, 7.2, and also 8.0Programming language used:C++, FORTRAN77: gcc 2.96 or 2.95.2 (also 3.2) compiler suite with g++ and g77Size of the package: 7.3 MB directory including example programs (2 MB compressed distribution archive), without ROOT libraries (additional 43 MB).No. of bytes in distributed program, including test data, etc.: 2 024 425Distribution format: tar gzip fileAdditional disk space required: Depends on the analyzed particle: 40 MB in the case of τ lepton decays (30 decay channels, 594 histograms, 82-pages booklet).Keywords: particle physics, decay simulation, Monte Carlo methods, invariant mass distributions, programs comparisonNature of the physical problem: The decays of individual particles are well defined modules of a typical Monte Carlo program chain in high energy physics. A fast, semi-automatic way of comparing results from different programs is often desirable, for the development of new programs, to check correctness of the installations or for discussion of uncertainties.Method of solution: A typical HEP Monte Carlo program stores the generated events in the event records such as HEPEVT or PYJETS. MC-TESTER scans, event by event, the contents of the record and searches for the decays of the particle under study. The list of the found decay modes is successively incremented and histograms of all invariant masses which can be calculated from the momenta of the particle decay products are defined and filled. The outputs from the two runs of distinct programs can be later compared. A booklet of comparisons is created: for every decay channel, all histograms present in the two outputs are plotted and parameter quantifying shape difference is calculated. Its maximum over every decay channel is printed in the summary table.Restrictions on the complexity of the problem: For a list of limitations see Section 6.Typical running time: Varies substantially with the analyzed decay particle. On a PC/Linux with 2.0 GHz processors MC-TESTER increases the run time of the τ-lepton Monte Carlo program TAUOLA by 4.0 seconds for every 100000 analyzed events (generation itself takes 26 seconds). The analysis step takes 13 seconds; processing takes additionally 10 seconds. Generation step runs may be executed simultaneously on multi-processor machines.Accessibility: web page: http://cern.ch/Piotr.Golonka/MC/MC-TESTER e-mails: Piotr.Golonka@CERN.CH, T.Pierzchala@friend.phys.us.edu.pl, Zbigniew.Was@CERN.CH.     

SPICE: Simulation Package for Including Flavor in Collider Events

We describe SPICE: Simulation Package for Including Flavor in Collider Events. SPICE takes as input two ingredients: a standard flavor-conserving supersymmetric spectrum and a set of flavor-violating slepton mass parameters, both of which are specified at some high “mediation” scale. SPICE then combines these two ingredients to form a flavor-violating model, determines the resulting low-energy spectrum and branching ratios, and outputs HERWIG and SUSY Les Houches files, which may be used to generate collider events. The flavor-conserving model may be any of the standard supersymmetric models, including minimal supergravity, minimal gauge-mediated supersymmetry breaking, and anomaly-mediated supersymmetry breaking supplemented by a universal scalar mass. The flavor-violating contributions may be specified in a number of ways, from specifying charges of fields under horizontal symmetries to completely specifying all flavor-violating parameters. SPICE is fully documented and publicly available, and is intended to be a user-friendly aid in the study of flavor at the Large Hadron Collider and other future colliders.

Program summary

Program title: SPICECatalogue identifier: AEFL_v1_0Program summary URL:http://cpc.cs.qub.ac.uk/summaries/AEFL_v1_0.htmlProgram obtainable from: CPC Program Library, Queen's University, Belfast, N. IrelandLicensing provisions: Standard CPC licence, http://cpc.cs.qub.ac.uk/licence/licence.htmlNo. of lines in distributed program, including test data, etc.: 8153No. of bytes in distributed program, including test data, etc.: 67 291Distribution format: tar.gzProgramming language: C++Computer: Personal computerOperating system: Tested on Scientific Linux 4.xClassification: 11.1External routines: SOFTSUSY [1,2] and SUSYHIT [3]Nature of problem: Simulation programs are required to compare theoretical models in particle physics with present and future data at particle colliders. SPICE determines the masses and decay branching ratios of supersymmetric particles in theories with lepton flavor violation. The inputs are the parameters of any of several standard flavor-conserving supersymmetric models, supplemented by flavor-violating parameters determined, for example, by horizontal flavor symmetries. The output are files that may be used for detailed simulation of supersymmetric events at particle colliders.Solution method: Simpson's rule integrator, basic algebraic computation.Additional comments: SPICE interfaces with SOFTSUSY and SUSYHIT to produce the low energy sparticle spectrum. Flavor mixing for sleptons and sneutrinos is fully implemented; flavor mixing for squarks is not included.Running time: <1 minute. Running time is dominated by calculating the possible and relevant three-body flavor-violating decays of sleptons, which is usually 10-15 seconds per slepton.References:

[1]

B.C. Allanach, Comput. Phys. Commun. 143 (2002) 305, arXiv:hep-ph/0104145.

[2]

B.C. Allanach, M.A. Bernhardt, arXiv:0903.1805 [hep-ph].

[3]

A. Djouadi, M.M. Muhlleitner, M. Spira, Acta Phys. Pol. B 38 (2007) 635, arXiv:hep-ph/0609292.

     

2HDMC — two-Higgs-doublet model calculator

We describe version 1.0.6 of the public C++ code 2HDMC, which can be used to perform calculations in a general, CP-conserving, two-Higgs-doublet model (2HDM). The program features simple conversion between different parametrizations of the 2HDM potential, a flexible Yukawa sector specification with choices of different Z₂-symmetries or more general couplings, a decay library including all two-body — and some three-body — decay modes for the Higgs bosons, and the possibility to calculate observables of interest for constraining the 2HDM parameter space, as well as theoretical constraints from positivity and unitarity. The latest version of the 2HDMC code and full documentation is available from: http://www.isv.uu.se/thep/MC/2HDMC.

New version program summary

Program title: 2HDMCCatalogue identifier: AEFI_v1_1Program summary URL:http://cpc.cs.qub.ac.uk/summaries/AEFI_v1_1.htmlProgram obtainable from: CPC Program Library, Queen's University, Belfast, N. IrelandLicensing provisions: GNU GPLNo. of lines in distributed program, including test data, etc.: 12 110No. of bytes in distributed program, including test data, etc.: 92 731Distribution format: tar.gzProgramming language: C++Computer: Any computer running LinuxOperating system: LinuxRAM: 5 MbCatalogue identifier of previous version: AEFI_v1_0Journal reference of previous version: Comput. Phys. Comm. 180 (2010) 189Classification: 11.1External routines: GNU Scientific Library (http://www.gnu.org/software/gsl/)Does the new version supersede the previous version?: YesNature of problem: Determining properties of the potential, calculation of mass spectrum, couplings, decay widths, oblique parameters, muon g−2, and collider constraints in a general two-Higgs-doublet model.Solution method: From arbitrary potential and Yukawa sector, tree-level relations are used to determine Higgs masses and couplings. Decay widths are calculated at leading order, including FCNC decays when applicable. Decays to off-shell vector bosons are obtained by numerical integration. Observables are computed (analytically or numerically) as function of the input parameters.Reasons for new version: Improved calculation of the oblique parameters.Summary of revisions: The computation of the oblique parameters has been improved to give reliable results in the case of degenerate masses for the Higgs bosons. Another issue in the oblique parameter calculation, affecting the numerical values of S, U, V, and X (independently of the Higgs boson masses), has been corrected.Restrictions: CP-violation is not treated.Running time: Less than 0.1 s on a standard PC.     

SuperIso v3.0, flavor physics observables calculations: extension to NMSSM

SuperIso v3.0 is a public program for evaluation of flavor physics observables in the minimal supersymmetric extension of the Standard Model (MSSM) and the next to minimal supersymmetric extension of the Standard Model (NMSSM). SuperIso v3.0 incorporates many flavor observables such as the inclusive branching ratio of B→Xsγ, the isospin asymmetry of B→K^∗γ, the branching ratio of Bs→μ⁺μ⁻, the branching ratio of B→τντ, the branching ratio of B→Dτντ, the branching ratio of K→μνμ and the branching ratios of Ds→τντ and Ds→μνμ. The calculation of the branching ratio of B→Xsγ includes NNLO Standard Model contributions. The program also computes the muon anomalous magnetic moment (g−2). Seven sample models are included in the package, namely mSUGRA, NUHM, AMSB and GMSB for the MSSM, and CNMSSM, NGMSB and NNUHM for the NMSSM. SuperIso uses a SUSY Les Houches Accord file (SLHA1 or SLHA2) as input, which can be either generated automatically by the program via a call to external spectrum calculators (SOFTSUSY, ISAJET or NMSSMTools), or provided by the user.

New version program summary

Program title:SuperIso v3.0Catalogue identifier: AEAN_v3_0Program summary URL:http://cpc.cs.qub.ac.uk/summaries/AEAN_v3_0.htmlProgram obtainable from: CPC Program Library, Queen's University, Belfast, N. IrelandLicensing provisions: GNU General Public LicenceNo. of lines in distributed program, including test data, etc.: 6869No. of bytes in distributed program, including test data, etc.: 42 627Distribution format: tar.gzProgramming language: C (C99 Standard compliant)Computer: 32- or 64-bit PC, MacOperating system: Linux, MacOSRAM: less than 1 MBClassification: 11.6External routines: ISASUGRA/ISAJET, SOFTSUSY and/or NMSSMToolsDoes the new version supersede the previous version?: YesNature of problem: Calculation of flavor physics observables as well as the muon anomalous magnetic moment in the Minimal Supersymmetric Standard Model with minimal flavor violation and in the Next to Minimal Supersymmetric Standard Model, in order to derive constraints on the supersymmetric parameter spaces.Solution method:SuperIso uses a SUSY Les Houches Accord (SLHA1 or SLHA2) file, which can be either generated automatically via a call to SOFTSUSY, ISAJET or NMSSMTools, or provided by the user. This file contains the masses, mixings and couplings of the supersymmetric particles. SuperIso then computes the most constraining flavor physics observables and the muon (g−2). SuperIso is able to perform the calculations in different supersymmetry breaking scenarios, such as mSUGRA, NUHM, AMSB and GMSB, as well as constrained NMSSM scenarios such as CNMSSM, NNUHM and NGMSB.Reasons for new version:SuperIso has been extended to the next to minimal supersymmetric extension of the Standard Model (NMSSM). The implemented routines are therefore extensively modified.Summary of revisions:

•

Improvement of the SLHA2 reader.

•

Replacement of “float” variables by “double”.

•

Implementation of an interface with NMSSMTools.

•

Extension of the calculation of flavor observables as well as the muon anomalous magnetic moment to NMSSM.

•

Addition of three different NMSSM scenarios: CNMSSM, NGMSB and NNUHM.

•

Three sample main programs have been added: cnmssm.c, ngmsb.c and nnuhm.c. Additional instructions to use them are given when running them without arguments.

Unusual features: The code is very flexible, and new observables can be added easily.Running time: Less than 1 sec     

FeynHiggs: A program for the calculation of MSSM Higgs-boson observables - Version 2.6.5

FeynHiggs is a Fortran code for the calculation of physical observables in the field of high-energy physics. FeynHiggs calculates various observables in the Higgs sector of the Minimal Supersymmetric Standard Model (MSSM) for real or complex parameters. These observables comprise Higgs-boson masses, mixing angles, couplings, Tevatron/LHC production cross-sections, branching ratios, as well as some additional observables such as Δρ, MW, the effective leptonic weak mixing angle, μ(g−2), BR(b→sγ), electric dipole moments.

New version program summary

Program title: FeynHiggs v2.6.5Catalogue identifier: ADKT_v2_0Program summary URL:http://cpc.cs.qub.ac.uk/summaries/ADKT_v2_0.htmlProgram obtainable from: CPC Program Library, Queen's University, Belfast, N. IrelandLicensing provisions: GPLNo. of lines in distributed program, including test data, etc.: 156 455No. of bytes in distributed program, including test data, etc.: 1 058 758Distribution format: tar.gzProgramming language: Fortran 77, C, MathematicaComputer: Intel/AMD, PowerPC, AlphaOperating system: Linux, Windows (Cygwin), Mac OS, Tru64 UnixRAM: insignificant (typically a few MB)Classification: 11.5Catalogue identifier of previous version: ADKT_v1_0Journal reference of previous version: Comput. Phys. Comm. 124 (2000) 76Does the new version supersede the previous version?: YesNature of problem: The experimental searches for Higgs bosons have to be compared with theory predictions at a high level of accuracy. Radiative corrections are especially important in the Minimal Supersymmetric Standard Model (MSSM).Solution method: High-precision calculations (mostly based on the Feynman-diagrammatic approach) for various Higgs-boson observables in the MSSM are implemented. The program takes the basic model parameters as input and evaluates many MSSM Higgs-boson observables relevant for experimental Higgs-boson physics.Reasons for new version: Existing calculations have been refined (e.g., by the inclusion of complex parameters); many new observables have been implemented.Summary of revisions: Input parameters are now allowed to be complex, and non-minimal flavor violation is permitted. Besides the Higgs-boson masses also branching ratios and hadron-collider production cross-sections are calculated.Restrictions: Several observables are still missing and have to be implemented.Running time: The command-line frontend, which invokes essentially all FeynHiggs functions, takes 200 ms a single point in parameter space On a 3 GHz AMD CPU.     

HiggsBounds: Confronting arbitrary Higgs sectors with exclusion bounds from LEP and the Tevatron

HiggsBounds is a computer code that tests theoretical predictions of models with arbitrary Higgs sectors against the exclusion bounds obtained from the Higgs searches at LEP and the Tevatron. The included experimental information comprises exclusion bounds at 95% C.L. on topological cross sections. In order to determine which search topology has the highest exclusion power, the program also includes, for each topology, information from the experiments on the expected exclusion bound, which would have been observed in case of a pure background distribution. Using the predictions of the desired model provided by the user as input, HiggsBounds determines the most sensitive channel and tests whether the considered parameter point is excluded at the 95% C.L. HiggsBounds is available as a Fortran 77 and Fortran 90 code. The code can be invoked as a command line version, a subroutine version and an online version. Examples of exclusion bounds obtained with HiggsBounds are discussed for the Standard Model, for a model with a fourth generation of quarks and leptons and for the Minimal Supersymmetric Standard Model with and without CP-violation. The experimental information on the exclusion bounds currently implemented in HiggsBounds will be updated as new results from the Higgs searches become available.

Program summary

Program title: HiggsBoundsCatalogue identifier: AEFF_v1_0Program summary URL:http://cpc.cs.qub.ac.uk/summaries/AEFF_v1_0.htmlProgram obtainable from: CPC Program Library, Queen's University, Belfast, N. IrelandLicensing provisions: Standard CPC licence, http://cpc.cs.qub.ac.uk/licence/licence.htmlNo. of lines in distributed program, including test data, etc.: 55 733No. of bytes in distributed program, including test data, etc.: 1 986 213Distribution format: tar.gzProgramming language: Fortran 77, Fortran 90 (two code versions are offered).Computer: HiggsBounds can be built with any compatible Fortran 77 or Fortran 90 compiler. The program has been tested on x86 CPUs running under Linux (Ubuntu 8.04) and with the following compilers: The Portland Group Inc. Fortran compilers (pgf77, pgf90), the GNU project Fortran compilers (g77, gfortran).Operating system: LinuxRAM: minimum of about 6000 kbytes (dependent on the code version)Classification: 11.1External routines: HiggsBounds requires no external routines/libraries. Some sample programs in the distribution require the programs FeynHiggs 2.6.x or CPsuperH2 to be installed (see “Subprograms used”).Subprograms used:

Cat Id	Title	Reference
ADKT_v2_0	FeynHiggsv2.6.5	CPC 180(2009)1426
ADSR_v2_0	CPsuperH2.0	CPC 180(2009)312

Full-size table