Bug in the B0 implementation in AddLoopFunctions.f90

Report the bugs you found
Post Reply
PhilipDiessner
Posts: 4
Joined: 14. Jul 2017, 12:07

Bug in the B0 implementation in AddLoopFunctions.f90

Post by PhilipDiessner » 14. Jul 2017, 12:17

Hi Florian,
The triviality check for the B0 function misses the condition for zero momentum:

Code: Select all

Complex(dp) Function SA_B0(xp, xm1, xm2)
 Implicit None
  Real(dp), Intent(in) :: xp, xm1, xm2
  
  If ((xm1.eq.0._dp).and.(xm2.eq.0._dp)) Then 
    SA_B0=0._dp
    Return
  End if
  ...
This screws up the neutrino contribution to Z self energy and probably some other things.

Cheers,
Philip

FStaub
Site Admin
Posts: 822
Joined: 13. Apr 2016, 14:05

Re: Bug in the B0 implementation in AddLoopFunctions.f90

Post by FStaub » 15. Jul 2017, 13:24

Hi Philipp,

good question why these lines are there.. this check for zero masses (together with the one if p2 is zero or not), is done on the SPheno side. Thanks for pointing this out.

Cheers,
Florian

Post Reply