SUSY-HIT Version 1.5 - Changes with respect to Version 1.3

Changes in SDECAY (now version 1.5a)

  • Changes by Margarete Mühlleitner: The 3-body decays of the neutralinos and the charginos now include the mass effects of the muon, the charm quark and the strange quark.

    • Changes in SDECAY (now version 1.5)

  • Changes by Ramona Gröber: implementation of the stop 3-body decay into W b neutralino1 with FV couplings; changes to treat the transition region between the 4-body and this 3-body decay correctly. Changes have been made in sdecay.f and lightst4bod.f. New are the files lighst3bod.f und Xvegas.f. (22 Feb 2015)

    • Changes in SDECAY (now version 1.4)

  • Change in definition of scaleofewsb (now not read in as extval(0) any more, as this leads to conflicts). Thanks to R. Gröber. (20 August 2014)
  • Changes by Ramona Gröber in order to implement the FCNC stop decays. (20 Aug 2014)

    • New Fortran file lightst4bod.f

  • A new Fortran file lightst4bod.f has been added to the package. It calculates the FCNC stop 4-body decays. The makefile has been changed accordingly.

    • Changes in HDECAY (version 3.4)

  • In order to support the gfortran compiler change of 'COMPLEX FUNCTION func*16' to 'COMPLEX*16 FUNCTION func'. (20 Aug 2014)

    • Changes in hgaga.f

  • In order to support the gfortran compiler change of 'COMPLEX FUNCTION func*16' to 'COMPLEX*16 FUNCTION func'. (20 Aug 2014)

    • Changes in suspect2.f (version 2.41)

  • Several changes in order to support the gfortran compiler. (20 Aug 2014)

    • SUSY-HIT Version 1.3 - Changes with respect to Version 1.2

    The newest versions of HDECAY and SuSpect have been linked in the SUSY-HIT Version 1.3. SuSpect by default now reads in the parameters from an input file in the SLHA format.

    Changes in SDECAY (now version 1.3b)

  • Typo fix in the function SD_ftotqcd, which plays a role for the QCD corrections to the light squark decays. The fixed typo in particular changes the QCD corrections for compressed spectra. Thanks to W. Hollik, J. Lindert, D. Pagani. (14 July 2012)
  • Some small changes such that in case SuSpect is taken as spectrum calculator, SuSpect now reads in the input file in a SLHA format. Therefore suspect2.in has been taken off and replaced by the input file suspect_slha.in which is in the SLHA format.

    • Changes in HDECAY (now version 3.4)

  • The old HDECAY 3.303 version has been replaced with the new version 3.4.

    • SUSY-HIT Version 1.2 - Changes with respect to Version 1.1

    Changes in SDECAY (now version 1.3a)

  • The new version can now also treat negative chargino masses. Thanks to Valery Zukhov for pointing out the problem.
  • Bugs fixed in the QCD corrections. This concerns the subroutines SD_ftfunctions and SD_fbfunctions.

    • Changes in HDECAY (now version 3.303)

  • The old HDECAY version has been replaced with the new version 3.303. Apart from the replacement of hdecay.f this implicated also the replacement of hgaga.f
  • Bug fixed in gaugino subroutine. In case the input is provided by an input file in the SLHA format, the Higgs-neutralino-neutralino couplings were not defined correctly.
  • This HDECAY version can now also treat negative chargino masses.

    • Changes in SuSpect (now version 2.41)

  • The old Suspect2 version has been replaced with the new version 2.41. Apart from the replacement of suspect2.f this implicated also the replacement of the files twoloophiggs.f, bsg.f and suspect2.in

    • SUSY-HIT Version 1.1 - Changes with respect to Version 1.0

    Changes in SDECAY (now version 1.3)

  • Some changes so that the decays into an LSP gravitino within GMSB models are calculated, also in case SuSpect is not linked. Some further decays within GMSB models have been added. Existing ones have been corrected. Thanks to CMS SUSYBSM.
  • For initialization, all partial and total widths are set to zero. Thanks to CMS SUSYBSM.
  • Correction of some PDG codes of leptons and sleptons in the SLHA output of the MSSM Higgs particle decays. Thanks to CMS SUSYBSM.
  • Correction of the PDG codes for stau1+, stau2+ in the Snu_tauL decay block of the SLHA output. Thanks to CMS SUSYBSM.
  • Some corrections in the output of SUSY-HIT, so that now also parameter values, which are 0 (and are allowed to be 0), are given out. Thanks to CMS SUSYBSM.
  • The value for extval(24) is now correctly read in. It is defined as M_A^2 in the SLHA. Thanks to Johan Rathsman.
  • Changes in the intitialization of the variables hmixval, minval, extval, m_softval of the corresponding common blocks defined in the SLHA input in order to jump over them in case they are not defined in the SLHA input file.
  • Added a comma at the end of the first line of the common block SD_sup2body. Added a comma at the end of the first line of the common block SD_sel2body. Put a comma in the common block SD_snel2body between brsnellneut and brsnellchar. Thanks to Rémi Lafaye.
  • Bug corrected: The PDG code for "h" is set to 25. (Was 26 before, which is wrong.) Thanks to Rémi Lafaye.
  • Correction in the SLHA input/output concerning soft SUSY breaking parameters. Thanks to Phillip Grajek.

    Changes in HDECAY

  • HDECAY version 3.200 has been replaced by version 3.302. This version now also needs the file hgaga.f which has been added to the SUSYHIT files on this homepage.
    The new HDECAY version includes the fully massive (SUSY-)QCD corrections to Higgs -> gamma gamma and the delta_mb effects in charged Higgs decays. Further changes are:
  • Correction of some PDG codes of leptons and sleptons in the SLHA output of the SM/MSSM Higgs particle decays. Thanks to CMS SUSYBSM.
  • The value for extval(24) is now correctly read in. It is defined as M_A^2 in the SLHA. Thanks to Johan Rathsman.
  • Changes in the initialization of the variables hmixval, minval, extval, m_softval of the corresponding common blocks defined in the SLHA input in order to jump over them in case they are not defined in the SLHA input file.
  • Correction in the setting of the stop, sbottom, stau mixing angles.
  • Further changes/corrections in the SLHA input and output.

    • Changes in susylha.f, the file to read in an SLHA input in HDECAY

  • The values minval, extval, hmixval, msoftval of the corresponding SLHA format blocks have been initialized differently than in the previous version in order to jump over undefined variables.

    • Changes in SuSpect

  • The value for extval(24) is now correctly read in. It is defined as M_A^2 in the SLHA. Thanks to Johan Rathsman.
  • The variable sgnmu0 is now correctly related to extval(23).
  • The variable extval(26) is now correctly related to M_A.
  • The value imod(2) is now correctly related to the various minimal models according to the SLHA format.

    • HDECAY within SUSY-HIT: Changes with respect to the original version HDECAY 3.200

  • HDECAY has become a subroutine.
  • The name of the SLHA input file to be read in is now "slhaspectrum.in". In the original version it is "slha.in".
  • There is no input file hdecay.in any more. HDECAY gets the necessary parameters from the SLHA input file. The parameter values not given in the SLHA input file are read in by SDECAY from the susyhit.in input file and passed on to HDECAY. The HDECAY flags are hard-coded: The decays of all MSSM Higgs bosons are calculated by default at higher order, as well as off-shell decays are included and the decays into SUSY particles.
  • For the sake of passing on parameters read in by SDECAY from susyhit.in the following common block has been created: COMMON/SUSYHITIN/flagshsin,amsin,amcin,ammuonin,alphin,gamwin,gamzin,vusin,vcbin,rvubin
  • HDECAY does not provide an output file within the SUSY-HIT package.
  • In the subroutine susycp_hdec the call for the subroutines haber and feynhiggs has been taken off. They are not needed within the program package SUSY-HIT.
  • Some commenting lines added.
  • The lines cw2calc = amw**2/amz**2
    sw2calc = 1-cw2calc
    cwcalc = dsqrt(cw2calc)
    swcalc = dsqrt(sw2calc)
    vewsb = 1.D0/dsqrt(dsqrt(2.D0)*gf)
    g2ew = 2*amw/vewsb
    g1ew = g2ew*swcalc/cwcalc
    and the command STOP have been taken out.
  • The lines
    write(nout,52) 1,amudrbar,'# mu(Q)'
    write(nout,52) 2,tgbet,'# tanbeta(Q)'
    have been changed to
    write(nout,52) 1,amudrbar,'mu(Q)'
    write(nout,52) 2,tgbet,'tanbeta(Q)'
  • The lines
    write(nout,53) 1,1,au,'A_u(Q) DRbar'
    write(nout,53) 2,2,au,'A_c(Q) DRbar'
    write(nout,53) 3,3,au,'A_t(Q) DRbar'
    have been changed to
    write(nout,53) 1,1,auval(1,1),'A_u(Q) DRbar'
    write(nout,53) 2,2,auval(2,2),'A_c(Q) DRbar'
    write(nout,53) 3,3,auval(3,3),'A_t(Q) DRbar'
    and this
    write(nout,53) 1,1,ad,'A_d(Q) DRbar'
    write(nout,53) 2,2,ad, 'A_s(Q) DRbar'
    write(nout,53) 3,3,ad,'A_b(Q) DRbar'
    to
    write(nout,53) 1,1,adval(1,1),'A_d(Q) DRbar'
    write(nout,53) 2,2,adval(2,2),'A_s(Q) DRbar'
    write(nout,53) 3,3,adval(3,3),'A_b(Q) DRbar'
    and this
    write(nout,53) 1,1,al,'A_e(Q) DRbar'
    write(nout,53) 2,2,al, 'A_mu(Q) DRbar'
    write(nout,53) 3,3,al,'A_tau(Q) DRbar'
    to
    write(nout,53) 1,1,aeval(1,1),'A_e(Q) DRbar'
    write(nout,53) 2,2,aeval(2,2),'A_mu(Q) DRbar'
    write(nout,53) 3,3,aeval(3,3),'A_tau(Q) DRbar'

    SDECAY within SUSY-HIT: Changes with respect to the original version SDECAY 1.2

  • As for the spectrum, SDECAY does not read in any more the input file "sdecay.in", but the newly created input file "susyhit.in".
  • Since there is no input file "sdecay.in" any more, the former flags for the original program SDECAY have been hard-coded. In detail, they have been set to
    flagoutput = 1.D0, flagqcd = 1.D0, flagmulti = 1.D0, flagloop = 1.D0, flagsusytop = 1.D0, flagnlspgmsb = 1.D0, idrbar = 1, scalechoice = 1.D0, scaleofcoupling = 100.D0
    This means that the QCD corrections to 2-body decays involving coloured particles will be calculated by default as well as multi-body and loop-induced decays, the decays of the NLSP in GMSB models and the top decays. The running strong coupling constant and quark masses will be evaluated in the DRbar scheme at the EWSB scale.
  • As for the mass spectrum and soft SUSY breaking parameters, SDECAY reads in by default an SLHA input file, which has to be called "slhaspectrum.in". It writes out by default an SLHA output file which is called "susyhit_slha.out".
  • To get the total widths and branching ratios of the MSSM Higgs bosons, SDECAY calls now HDECAY. To pass the needed results from HDECAY to SDECAY, some common blocks of HDECAY have been implemented in SDECAY: COMMON/WIDTHA_HDEC, /WIDTHHL_HDEC, /WIDTHHH_HDEC, /WIDTHHC_HDEC, /WISUSY_HDEC, /WISFER_HDEC, /HD_golddec. The MSSM Higgs boson total widths and branching ratios are given out in "susyhit_slha.out".
  • A new common block
    COMMON/SUSYHITIN/flagshsin,amsin,amcin,ammmuonin,alphin,gamwin,gamzin,vusin,vcb
    has been created. It passes on the flag "flagshsin" which decides whether SuSpect is linked or not and which is read in from "susyhit.in". And it passes on the values for the strange and charm quark masses, for the muon mass, for alpha, the total W and Z boson widths and some CKM matrix elements to HDECAY which are needed there. All theses parameters are read in from "susyhit.in".
  • Change of the common block COMMON/SD_INPUT/flaginput,ninlha to COMMON/SD_INPUT/ninlha. Since the flag "flaginput" is not needed any more.
  • Several comment lines have been changed according to the needs and also the comments in the output file.

    • SuSpect2 within SUSY-HIT: Changes with respect to the original version SuSpect2 2.34

  • The SuSpect output file in the SLHA format is now called "slhaspectrum.in". It will be read in by HDECAY and SDECAY.
  • The functions FO_HDEC and SUBH_HDEC have been renamed F0_TLH and SUBH_TLH in order to avoid clashes with the same names appearing in hdecay.f.
  • In the subroutine twoloophiggs.f the suffix _HDEC of the following functions has been changed to _TLH in order to avoid clashes with the same functions also appearing in HDECAY: f0_HDEC, t_HDEC, gfun_HDEC, delmb_HDEC, SUBH_HDEC.
  • The charm quark mass value mc=1.42d0 has been changed to mc=1.40d0 whereever it appears in order to be in accordance with the usually set value in the input file.
  • The lines
    write(nout,61) 1,'SuSpect # RGE +Spectrum calculator'
    write(nout,61) 2,'2.33 # version number'
    have been changed to
    write(nout,61) 1,'SuSpect # RGE +Spectrum calculator'
    write(nout,61) 2,'2.34 # version number'
    Note that the version number has been adapted!
    Margarete Mühlleitner
    Last modified: Wed Sept 28 10:34:12 MET 2016