Hi Florian, Javier invited me to try to figure this out (we were in contact as the author of 1601.03737). I am not sure I understand this line: If (gt1.ne.gt2) gam = 2._dp*gam When calculating the width for a scalar decaying to two fermions, one simply writes down the feynman diagram for $\phi > f1 ...