Spheno Higgs Bounds output is missing a pseudo scalar higgs

[jackaraz]

Hi Florian

I'm using SARAh 4.9.3 (tried also 12.2 same issue). I realised that SPhenos higgsbounds output doesnot contain the correct number of pseudoscalar higgs. here is an example with secluded mssm where higgs sector from slha file is

Code: Select all

        25     0.00000000E+00  # hh_1
        35     0.00000000E+00  # hh_2
        45     1.23866041E+02  # hh_3
        55     1.79366743E+02  # hh_4
        65     2.37668646E+03  # hh_5
        75     2.42229613E+03  # hh_6
        36     2.44099648E+03  # Ah_2
        46     5.33754868E+00  # Ah_3
        56     7.74173876E+00  # Ah_4
        66     4.25065369E+01  # Ah_5
        76     2.37682768E+03  # Ah_6
        37     2.37889584E+03  # Hpm_2

However the higgsbounds subrutine is missing the Ah_2

Code: Select all

 Subroutine WriteHiggsBounds
Implicit None
Open(87,file="[i]MH_GammaTot.dat[/i]",status="unknown")
Open(88,file="MHplus_GammaTot.dat",status="unknown")
Open(89,file="effC.dat",status="unknown")
Open(90,file="BR_H_NP.dat",status="unknown")
Open(91,file="BR_Hplus.dat",status="unknown")
Open(92,file="BR_t.dat",status="unknown")
Open(93,file="LEP_HpHm_CS_ratios.dat",status="unknown")
Write(87,"(I1)",advance="No") 1
Write(88,"(I1)",advance="No") 1
Write(87,"(2e16.8)",advance="No") Mhh(1)
Write(87,"(2e16.8)",advance="No") Mhh(2)
Write(87,"(2e16.8)",advance="No") Mhh(3)
Write(87,"(2e16.8)",advance="No") Mhh(4)
Write(87,"(2e16.8)",advance="No") Mhh(5)
Write(87,"(2e16.8)",advance="No") Mhh(6)
Write(87,"(2e16.8)",advance="No") MAh(3)
Write(87,"(2e16.8)",advance="No") MAh(4)
Write(87,"(2e16.8)",advance="No") MAh(5)

in case of MH_GammaTot.dat it is easy to fix but other files,such as effC.dat, has the same issue as well. by any chance do you have a simple patch for it? (I added new pid's in higgsbounds PDGnumbering.f90)

Thank you very much

best regards

[FStaub]

Hi,

I actually think that there is a mistake in the model files: this model has two neutral Goldstones because two neutral gauge bosons are there, i.e. it should read in particles.m

Code: Select all

     {Ah   ,  {    Description -> "Pseudo-Scalar Higgs",
                 PDG -> {0, 0, 36, 46,56,66},
                 PDG.IX ->{0, 0, 102000001,102000002,102000003,102000004  } }},  

So, the appearance of Ah_2 in the mass block is wrong.

Cheers,
Florian