Page 1 of 1

Bug in the B0 implementation in AddLoopFunctions.f90

Posted: 14. Jul 2017, 12:17
by PhilipDiessner
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

Re: Bug in the B0 implementation in AddLoopFunctions.f90

Posted: 15. Jul 2017, 13:24
by FStaub
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