WPHACT 2.0: A fully massive Monte Carlo generator for four fermion physics at ee colliders

WPHACT 2.0 is the new fully massive version of a MC program and unweighted event generator which computes all Standard Model processes with four fermions in the final state at e⁺e⁻ colliders. The program can now generate unweighted events for any subset of all four fermion final states in a single run, by making use of dedicated pre-samples which can cover the entire phase space. Improvements with respect to WPHACT 1.0 include the Imaginary Fermion Loop gauge restoring scheme, new phase space mappings, a new input system, the possibility to compute subsets of Feynman diagrams and options for including ISR via QEDPS, running α_QED, CKM mixing, resonances in channels.     

SOFTSUSY: A program for calculating supersymmetric spectra

SOFTSUSY is a program which accurately calculates the spectrum of superparticles in the Minimal Supersymmetric Standard Model (MSSM). The program solves the renormalisation group equations with theoretical constraints on soft supersymmetry breaking terms provided by the user. Weak-scale gauge coupling and fermion mass data (including one-loop finite MSSM corrections) are used as a boundary condition, as well as successful radiative electroweak symmetry breaking. The program can also calculate a measure of fine-tuning. The program structure has been designed to easily generalize to extensions of the MSSM . This article serves as a self-contained guide to prospective users, and indicates the conventions and approximations used. Sample results are compared with similar calculations in the literature.     

SDECAY: a Fortran code for the decays of the supersymmetric particles in the MSSM

We present the Fortran code SDECAY, which calculates the decay widths and branching ratios of all the supersymmetric particles in the Minimal Supersymmetric Standard Model, including higher order effects. Besides the usual two-body decays of sfermions and gauginos and the three-body decays of charginos, neutralinos and gluinos, we have also implemented the three-body decays of stops and sbottoms, and even the four-body decays of the stop; the important loop-induced decay modes are also included. The QCD corrections to the two-body decays involving strongly interacting particles and the dominant components of the electroweak corrections to all decay modes are implemented.

Program summary

Title of program: SDECAY Version 1.1a (March 2005)Catalogue identifier: ADVJProgram summary URL:http://cpc.cs.qub.ac.uk/summaries/ADVJProgram obtainable: CPC Program Library, Queen's University of Belfast, N. IrelandLicensing provisions: noneComputer for which the program is designed: Any with a Fortran77 systemOperating systems under which the program has been tested: Linux, UnixTypical running time: A few seconds on modern personal computers and workstationsProgramming language used: Fortran77No. of lines in distributed program, including test data, etc.: 59 621No. of bytes in distributed program, including test data, etc.: 338 478Distribution format: tar.gzMemory required to execute (with test data): 7.3 MBDistribution format: ASCIINature of physical problem: Numerical calculation of the decay widths and branching ratios of supersymmetric particles in the Minimal Supersymmetric Standard Model (MSSM). The program calculates two-, three- and four-body decays and loop decays. It includes the SUSY-QCD corrections to two-body decays involving strongly interacting particles. The top-quark decays within the MSSM are evaluated as well.Method of solution: Two-dimensional numerical integration of the analytic formulae for the double differential decay widths of the three-body decays. The other decay widths are calculated analytically.Restrictions on the complexity of the problem: In the higher order decay modes the total decay widths of the virtually exchanged (s)particles are not included in their respective propagators. The higher order decays are calculated when the two-body decays are kinematically closed.     

eett6f v. 1.0: A program for top quark pair production and decay into 6 fermions at linear colliders

The first version of a computer program eett6f for calculating cross sections of e⁺e⁻→6 fermions processes relevant for a -pair production and decay at centre of mass energies typical for linear colliders is presented. eett6f v. 1.0 allows for calculating both the total and differential cross sections at tree level of the Standard Model (SM). The program can be used as the Monte Carlo generator of unweighted events as well.     

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.     

EPW: A program for calculating the electron–phonon coupling using maximally localized Wannier functions

EPW (Electron–Phonon coupling using Wannier functions) is a program written in Fortran90 for calculating the electron–phonon coupling in periodic systems using density-functional perturbation theory and maximally localized Wannier functions. EPW can calculate electron–phonon interaction self-energies, electron–phonon spectral functions, and total as well as mode-resolved electron–phonon coupling strengths. The calculation of the electron–phonon coupling requires a very accurate sampling of electron–phonon scattering processes throughout the Brillouin zone, hence reliable calculations can be prohibitively time-consuming. EPW combines the Kohn–Sham electronic eigenstates and the vibrational eigenmodes provided by the Quantum ESPRESSO package (see Giannozzi et al., 2009 [1]) with the maximally localized Wannier functions provided by the wannier90 package (see Mostofi et al., 2008 [2]) in order to generate electron–phonon matrix elements on arbitrarily dense Brillouin zone grids using a generalized Fourier interpolation. This feature of EPW leads to fast and accurate calculations of the electron–phonon coupling, and enables the study of the electron–phonon coupling in large and complex systems.

Program summary

Program title: EPWCatalogue identifier: AEHA_v1_0Program summary URL:http://cpc.cs.qub.ac.uk/summaries/AEHA_v1_0.htmlProgram obtainable from: CPC Program Library, Queen's University, Belfast, N. IrelandLicensing provisions: GNU Public LicenseNo. of lines in distributed program, including test data, etc.: 304 443No. of bytes in distributed program, including test data, etc.: 1 487 466Distribution format: tar.gzProgramming language: Fortran 90Computer: Any architecture with a Fortran 90 compilerOperating system: Any environment with a Fortran 90 compilerHas the code been vectorized or parallelized?: Yes, optimized for 1 to 64 processorsRAM: Heavily system dependent, as small as a few MBSupplementary material: A copy of the “EPW/examples” directory containing the phonon binary files can be downloadedClassification: 7External routines: MPI, Quantum-ESPRESSO package [1], BLAS, LAPACK, FFTW. (The necessary Blas, Lapack and FFTW routines are included in the Quantum-ESPRESSO package [1].)Nature of problem: The calculation of the electron–phonon coupling from first-principles requires a very accurate sampling of electron–phonon scattering processes throughout the Brillouin zone; hence reliable calculations can be prohibitively timeconsuming.Solution method: EPW makes use of a real-space formulation and combines the Kohn–Sham electronic eigenstates and the vibrational eigenmodes provided by the Quantum-ESPRESSO package with the maximally localized Wannier functions provided by the wannier90 package in order to generate electron–phonon matrix elements on arbitrarily dense Brillouin zone grids using a generalized Fourier interpolation.Running time: Single processor examples typically take 5–10 minutes.References:

• [1] 

  P. Giannozzi, et al., J. Phys. Condens. Matter 21 (2009), 395502, http://www.quantum-espresso.org/.

     

TSIL: a program for the calculation of two-loop self-energy integrals

TSIL is a library of utilities for the numerical calculation of dimensionally regularized two-loop self-energy integrals. A convenient basis for these functions is given by the integrals obtained at the end of O.V. Tarasov's recurrence relation algorithm. The program computes the values of all of these basis functions, for arbitrary input masses and external momentum. When analytical expressions in terms of polylogarithms are available, they are used. Otherwise, the evaluation proceeds by a Runge-Kutta integration of the coupled first-order differential equations for the basis integrals, using the external momentum invariant as the independent variable. The starting point of the integration is provided by known analytic expressions at (or near) zero external momentum. The code is written in C, and may be linked from C/C++ or Fortran. A Fortran interface is provided. We describe the structure and usage of the program, and provide a simple example application. We also compute two new cases analytically, and compare all of our notations and conventions for the two-loop self-energy integrals to those used by several other groups.

Program summary

Title of program:TSILVersion number: 1.0Catalogue identifier: ADWSProgram summary URL:http://cpc.cs.qub.ac.uk/summaries/ADWSProgram obtainable from: CPC Program Library, Queen's University of Belfast, N. IrelandProgramming language: CPlatform: Any platform supporting the GNU Compiler Collection (gcc), the Intel C compiler (icc), or a similar C compiler with support for complex mathematicsNo. of lines in distributed program, including test data, etc.: 42 730No. of bytes in distributed program, including test data, etc.: 297 101Distribution format: tar.gzNature of physical problem: Numerical evaluation of dimensionally regulated Feynman integrals needed in two-loop self-energy calculations in relativistic quantum field theory in four dimensions.Method of solution: Analytical evaluation in terms of polylogarithms when possible, otherwise through Runge-Kutta solution of differential equations.Limitations: Loss of accuracy in some unnatural threshold cases that do not have vanishing masses.Typical running time: Less than a second.     

RacoonWW1.3: A Monte Carlo program for four-fermion production at ee colliders

We present the Monte Carlo generator RacoonWW that computes cross sections to all processes e⁺e⁻→4f and e⁺e⁻→4fγ and calculates the complete electroweak radiative corrections to e⁺e⁻→WW→4f in the electroweak Standard Model in double-pole approximation. The calculation of the tree-level processes e⁺e⁻→4f and e⁺e⁻→4fγ is based on the full matrix elements for massless (polarized) fermions. When calculating radiative corrections to e⁺e⁻→WW→4f, the complete virtual doubly-resonant electroweak corrections are included, i.e. the factorizable and non-factorizable virtual corrections in double-pole approximation, and the real corrections are based on the full matrix elements for e⁺e⁻→4fγ. The matching of soft and collinear singularities between virtual and real corrections is done alternatively in two different ways, namely by using a subtraction method or by applying phase-space slicing. Higher-order initial-state photon radiation and naive QCD corrections are taken into account. RacoonWW also provides anomalous triple gauge-boson couplings for all processes e⁺e⁻→4f and anomalous quartic gauge-boson couplings for all processes e⁺e⁻→4fγ.     

OptaDOS: A tool for obtaining density of states,core-level and optical spectra from electronic structure codes

We present OptaDOS, a program for calculating core-electron and low-loss electron energy loss spectra (EELS) and optical spectra along with total-, projected- and joint-density of electronic states (DOS) from single-particle eigenenergies and dipole transition coefficients. Energy-loss spectroscopy is an important tool for probing bonding within a material. Interpreting these spectra can be aided by first principles calculations. The spectra are generated from the eigenenergies through integration over the Brillouin zone. An important feature of this code is that this integration is performed using a choice of adaptive or linear extrapolation broadening methods which we show produces higher accuracy spectra than standard fixed-width Gaussian broadening. OptaDOS  may be straightforwardly interfaced to any electronic structure code. OptaDOS  is freely available under the GNU General Public licence from http://www.optados.org.     

carlomat, version 2 of the program for automatic computation of lowest order cross sections

Version 2 of carlomat, a program for automatic computation of the lowest order cross sections of multiparticle reactions, is described. The substantial modifications with respect to version 1 of the program include: generation of a single phase space parameterization for the Feynman diagrams of the same topology, an interface to parton density functions, improvement of the color matrix computation, the Cabibbo–Kobayashi–Maskawa mixing in the quark sector, the effective models including scalar electrodynamics, the Wtb$Wtb$ interaction with operators of dimension up to 5 and a general top–Higgs coupling. Moreover, some minor modifications have been made and several bugs in the program have been corrected.     

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.     

HypExp, a Mathematica package for expanding hypergeometric functions around integer-valued parameters

We present the Mathematica package HypExp which allows to expand hypergeometric functions around integer parameters to arbitrary order. At this, we apply two methods, the first one being based on an integral representation, the second one on the nested sums approach. The expansion works for both symbolic argument z and unit argument. We also implemented new classes of integrals that appear in the first method and that are, in part, yet unknown to Mathematica.

Program summary

Title of program:HypExpCatalogue identifier:ADXF_v1_0Program summary URL:http://cpc.cs.qub.ac.uk/summaries/ADXF_v1_0Program obtainable from: CPC Program Library, Queen's University of Belfast, N. IrelandLicence:noneComputers:Computers running Mathematica under Linux or WindowsOperating system:Linux, WindowsProgram language:MathematicaNo. of bytes in distributed program, including test data, etc.:739 410No. of lines in distributed program, including test data, etc.:89 747Distribution format:tar.gzOther package needed:the package HPL, included in the distributionExternal file required:noneNature of the physical problem:Expansion of hypergeometric functions around integer-valued parameters. These are needed in the context of dimensional regularization for loop and phase space integrals.Method of solution:Algebraic manipulation of nested sums and integral representation.Restrictions on complexity of the problem:Limited by the memory availableTypical running time:Strongly depending on the problem and the availability of libraries.     

CROSSV,a simple FORTRAN 77 program for calculating 2-dimensional experimental cross-variograms

A simple portable FORTRAN 77 program is proposed for computing experimental cross-variograms. The program is suitable especially for use on personal computers and can be translated easily into other programming languages.     

Analytical ray tracing system: Introducing art, a C-library designed for seismic applications

Ray tracing technique is an important tool not only to forward but also for inverse problems in Geophysics, which most of the seismic processing steps depend on. However, implementing ray tracing codes can be very time consuming. This article presents a computer library to trace rays in 2.5D media composed by a stack of layers. The velocity profile inside each layer is such that the eikonal equation can be analytically solved. Therefore, the ray tracing within such profile is made fast and accurate. The great advantage of an analytical ray tracing library is the numerical precision of the quantities computed and the fast execution of the implemented codes. Even though ray tracing programs exist for a long time, for example the seis88 package by ?ervený, most of those programs use a numerical approach to compute the ray. Regardless of the fact that numerical methods can solve more general problems, the analytical ones could be part of a more sophisticated simulation process, where the ray tracing time is completely relevant. We demonstrate the feasibility of our codes using several examples (Miqueles et al., 2013)  [1]. The library can also be used for other applications besides seismic, e.g., optics and tomography.     

CONFREG: a BASIC program for calculating and plotting confidence regions based on correlational analyses

Many observations encountered in biological and medical research are randomly distributed in bivariate scales, and thus not susceptible to simple regression analyses. Since such data are depicted by ellipses in scatter diagrams, a computer program to calculate the confidence regions for the means or the total data of bivariate samples was written in BASIC for correlational analyses. The program, based on the principal axes algorithm, plots the calculated confidence regions as an elliptic area, using the fitted equations for its major and minor axes. The program displays the sample parameters required to perform comparisons between different groups of experimental conditions.     

ISICS: A program for calculating K-, L- and M-shell cross sections from ECPSSR theory using a personal computer

A versatile and fast C-language program has been developed to exactly calculate ionization and x-ray production cross sections using the ECPSSR theory. All elements and any projectile can be selected for any range of energies and energy increments. Coulomb excitation form factors are numerically calculated in the program. Comparisons are presented between program calculations and those obtained by standard procedures from other sources. The program is available in compiled form and in source code for modification by the user. User options include selecting a relativistic calculation of the projectile ion velocity. Electron capture by the projectile is not incorporated into the program at this time.