Open hoodyn opened 3 years ago
Summary of the issues being solved and proposed changes:
Issues:
1: Call for a prediction of 'BR(KS->emu,mue)'
was giving an error due to a bad implementation.
2: BR(KL->emu,mue)
was calculated incorrectly, since simplified formulas were used which hold only for LF conserving cases.
Suggested updates
K='KS', l1='e', l2='mu'
is no longer an invalid combination of options in the function amplitudes_eff
.
Amplitudes for LFV decays (i.e. l1!=l2
) were not calculated correctly. Fixing the physics required some nontrivial extension of the code:
the original function function amplitudes
has been renamed to amplitudes_weak_eigst
as it actually does not yield the relevant amplitudes for KL -> l1+l2- nor KS->l1+l2-, but for the decay of weak eigenstate K0bar -> l1+ l2-. Thus the usage:
S_K0bar, P_K0bar = amplitudes_weak_eigst(par, wc, l1, l2)
Also the amplitude for K0->l1+l2- can be obtained from this function as
S_aux, P_aux = amplitudes_weak_eigst(par, wc, l2, l1)
S_K0 = -S_aux.conjugate()
P_K0 = P_aux.conjugate()
I will attach some updated notes explaining these matters in detail.
An entirely new function amplitudes(par, wc, K, l1, l2)
has been implemented. It simply makes proper linear combinations of amplitudes for K0 and K0bar which correspond to the approximate mass eigenstates KL, KS (in fact, CP eigenstates are used instead). In the special case l1 == l2
, this can be simplified to the expressions with imaginary or real parts of S,P, which had been there originally, adopted from ArXiv:1711.11030 (I guess). I implemented this special case separately inside the amplitudes, for the reason of both effectivity and better possibility to check that the new formulae are consistent with the old ones. If the maintainers find this separate treatment redundant, it can be removed, of course. The new amplitudes function now takes K
as a mandatory argument (to be 'KL'
or 'KS'
).
There is no longer a need to take real or imaginary parts of the amplitudes within the amplitudes_eff
function, which now serves purely to sum the short-distance amplitudes
with amplitudes_LD
in the LF-conserving cases. It also seems more consistent now to me. Concerning LF-conserving decays: Due to different sign conventions, amplitudes_eff
subtracts the LD and SD contributions in this version, but in fact the interference works the same way as it did before. (For 'KL'
, an extra minus sign had been present in original implementation of P
in amplitudes_LD
. For 'KS'
, the sign is compensated with S=-S_K0bar.real
inside amplitudes
function yields .
Furthermore, I made an additional small commit:
amplitudes_LD(par, K, l)
was originally giving long-distance contributions for both 'KL'
and 'KS'
, although it takes K
as an argument. Eventually, everything had been calculated correctly but the implementation seemed little awkward to me. So I've changed it in order to make things clearer.Here's my notes about the leptonic kaon decays in EFT. (I haven't found any suitable text in the literature.) Might be helpful to understand the changes in the implementation. kaon_decays.pdf
@hoodyn thank you again for all the work you have put into this! Since your results differ from those in published papers, I think it would be good to do some cross checks before merging this PR. Unfortunately, I don't have enough time at the moment to do that. So I will keep this PR open until I find time for this or somebody else can do some cross checks. Did you actually contact the authors of the paper that seems to contain mistakes?
@peterstangl Hi, I haven't contacted the authors of 1602.00881. I think that the mistakes in that paper can be qualified just as "big typos". This paper does not give formulae for KL,KS -> emu anyway, so I don't think it is important for this PR.
Btw, I used my branch of flavio to successfully crosscheck predictions of Ref. 1801.02895. The same crosscheck fails for the current master branch. This, however, says nothing about the problems I see in Ref. 1602.00881.
@hoodyn thank you again for all the work you have put into this! Since your results differ from those in published papers, I think it would be good to do some cross checks before merging this PR. Unfortunately, I don't have enough time at the moment to do that. So I will keep this PR open until I find time for this or somebody else can do some cross checks. Did you actually contact the authors of the paper that seems to contain mistakes?
Here's my notes about the leptonic kaon decays in EFT. (I haven'tfound any suitable text in the literature.) Might be helpful to understand the changes in the implementation. kaon_decays.pdf
@hoodyn I cross-checked your notes and found overall agreement. I found a few disagreements, which are, however, probably just typos. Nevertheless, I list them in the following:
m_K0
missing.S_sdmue S_sdemu - P_sdmue P_sdemu
).Concerning the implementation, I will comment directly on the code.
This is the (second) pull request aiming to fix the problems in calculations od LFV kaon decays to leptons, based on discussion in #138. (I'm aborting my original PR #140 which I polluted by bugs and wrong argumentation.)