On-shell condition for loop masses in the THDM-CPV model

[AndersKvellestad]

Hi Florian,

Unless I'm mistaken there's a bug in the SPheno code generated for THDM-CPV (and possibly other models). In the subroutine OneLoophh in LoopMassesTHDMCPV.f90 I think there's supposed to be an iterative calculation to satisfy the on-shell condition for eigenvalues of the loop-corrected mass matrix. However, in the code the self-energy is always calculated at p^2 = 0 for all four scalar states:

Code: Select all

i_count = 0 
Do  
i_count = i_count + 1 
test_m2 = mass2 
Do i1=1,4
PiSf(i1,:,:) = ZeroC 
p2 =  mass2(i1) 
If (i1.eq.1) p2 = 0._dp 
If (i1.eq.2) p2 = 0._dp 
If (i1.eq.3) p2 = 0._dp 
If (i1.eq.4) p2 = 0._dp 
Call Pi1Loophh(p2,MFd,MFd2,MFe,MFe2,MFu,MFu2,Mhh,Mhh2,MVZ,MVZ2,MHm,MHm2,              & 
& MVWm,MVWm2,cplcFdFdUhhL,cplcFdFdUhhR,cplcFeFeUhhL,cplcFeFeUhhR,cplcFuFuUhhL,           & 
& cplcFuFuUhhR,cplcgWmgWmUhh,cplcgWpCgWpCUhh,cplcgZgZUhh,cplUhhhhhh,cplUhhhhVZ,          & 
& cplUhhHmcHm,cplUhhHmcVWm,cplUhhcVWmVWm,cplUhhVZVZ,cplUhhUhhhhhh,cplUhhUhhHmcHm,        & 
& cplUhhUhhcVWmVWm,cplUhhUhhVZVZ,kont,PiSf(i1,:,:))
End Do
(...calculation of loop-corrected mass matrix...)
End Do

Consequently the loop ends after only two iterations, since the results of the first and second iterations are identical.

Kind regards,
Anders

[FStaub]

Hi,

the patch is attached.

Cheers,
Florian

[AndersKvellestad]

Thanks for the quick reply. This works like a charm.

Cheers,
Anders