2HDMC – two-Higgs-doublet model calculator

• Published in: Computer physics communications 2010, Vol.181 (1), p.189-205
• Main Authors: Eriksson, David, Rathsman, Johan, Stål, Oscar
• Format: Article
• Language: English
• Subjects: Fysik; Higgs physics; NATURAL SCIENCES; NATURVETENSKAP; Physics; Two-Higgs-doublet model
• ISSN: 0010-4655; 1879-2944; 1879-2944
• DOI: 10.1016/j.cpc.2009.09.011

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 2 -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 title: 2HDMC Catalogue identifier: AEFI_v1_0 Program summary URL: http://cpc.cs.qub.ac.uk/summaries/AEFI_v1_0.html Program obtainable from: CPC Program Library, Queen's University, Belfast, N. Ireland Licensing provisions: GNU GPL No. of lines in distributed program, including test data, etc.: 12 032 No. of bytes in distributed program, including test data, etc.: 90 699 Distribution format: tar.gz Programming language: C++ Computer: Any computer running Linux Operating system: Linux RAM: 5 Mb Classification: 11.1 External 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