On-shell condition for loop masses in the THDM-CPV model
Posted: 3. May 2016, 15:15
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:
Consequently the loop ends after only two iterations, since the results of the first and second iterations are identical.
Kind regards,
Anders
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
Kind regards,
Anders