CPsuperH: a computational tool for Higgs phenomenology in the minimal supersymmetric standard model with explicit CP violation

Abstract

We provide a detailed description of the Fortran code CPsuperH, a newly-developed computational package that calculates the mass spectrum and decay widths of the neutral and charged Higgs bosons in the Minimal Supersymmetric Standard Model with explicit CP violation. 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 Higgs-boson pole-mass shifts. The code CPsuperH is self-contained (with all subroutines included), is easy and fast to run, and is organized to allow further theoretical developments to be easily implemented. 1 1 The program may be obtained from http://theory.ph.man.ac.uk/~jslee/CPsuperH.html. The fact that the masses and couplings of the charged and neutral Higgs bosons are computed at a similar high-precision level makes it an attractive tool for Tevatron, LHC and LC studies, also in the CP-conserving case. Program summary Title of program: CPsuperH Catalogue number: ADSR Program summary URL: http://cpc.cs.qub.ac.uk/summaries/ADSR Program obtainable from: CPC Program Library, Queen's University of Belfast, N. Ireland Computer for which the program is designed and others on which it is operable: PC running under Linux and computers in Unix environment Programming Language used: Fortran77 High-speed storage required: No No. of bytes in distributed program, including test data, etc.: 49 729 Distribution format: tar gzip file No. of lines in combined program and test case: 7417 Separate documentation available: The program and the paper may be obtainable from http://theory.ph.man.ac.uk/~jslee/CPsuperH.html Keywords: Higgs bosons, supersymmetry, CP Nature of physical problem: The mass spectrum, decay widths and branching ratios of the neutral and charged Higgs bosons in the Minimal Supersymmetric Standard Model with explicit CP violation are calculated. 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 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. Method of solution: One-dimensional numerical integration for the Higgs-decay modes involving more than one massive gauge boson, iterative treatment of the threshold corrections, and the numerical diagonalization of the neutralino mass matrix. Restrictions on the complexity of the problem: Only the dominant leading-order QCD corrections which remain unaffected by CP violation have been included in this version. Typical running time: Less than 0.1 s Unusual features of the program: No