eHDECAY: An implementation of the Higgs effective Lagrangian into HDECAY

We present eHDECAY, a modified version of the program HDECAY which includes the full list of leading bosonic operators of the Higgs effective Lagrangian with a linear or non-linear realization of the electroweak symmetry and implements two benchmark composite Higgs models.Program summaryProgram title: eHDECAYCatalogue identifier: AEUJ_v1_0Program summary URL:http://cpc.cs.qub.ac.uk/summaries/AEUJ_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.: 71512No. of bytes in distributed program, including test data, etc.: 666827Distribution format: tar.gzProgramming language: Fortran77.Computer: Any with a Fortran77 system.Operating system: Linux, Unix.RAM: 0.5 MBClassification: 11.1.Nature of problem: Numerical calculation of the decay widths and branching ratios of a Higgs-like boson within four different parametrizations: the non-linear Lagrangian, the Strongly-Interacting Light Higgs (SILH) Lagrangian and the MCHM4 and MCHM5 Lagrangians. The Fortran program eHDECAY includes the most important higher-order QCD effects and in case of the SILH and composite Higgs parametrization also the electroweak (EW) higher order corrections. The user furthermore has the possibility to turn off these EW corrections.Solution method: The necessary input values, the choice of the parametrization and the values of the various couplings are set in the input file ehdecay.in. These are read in by the main routine ehdecay.f. The main routine calculates the decay widths and branching ratios through analytical formulae, by using several help routines (dmb.f, elw.f, feynhiggs.f, haber.f, hgaga.f, hgg.f, hsqsq.f, susylha.f). The calculated branching ratios and total width are given out in the files br.eff1 and br.eff2.Restrictions: The EW corrections are included only in the SILH and the composite Higgs parametrization in an approximate way. They are consistently not included in the non-linear case. This would require the explicit calculation of the EW higher order corrections in this framework. The program does not provide any distributions.Running time: Less than one second per point