flav-io / flavio

A Python package for flavour physics phenomenology in the Standard model and beyond
http://flav-io.github.io/
MIT License
71 stars 62 forks source link

NP contributions to R(D(*)) from scalar / tensor operators #157

Closed whereforebound closed 3 years ago

whereforebound commented 3 years ago

This issue is related to the implementation of the lepton-flavour universality ratios R(D) and R(D) in the newest version of flavio. There are analytic expressions for these observables used in the literature, given in terms of the SM values for R(D()) and the new-physics contributions which comprise the relevant Wilson coefficients and numerical factors capturing the form factors and kinematics. While the contributions from the vector operator seem to agree very well between Flavio and the analytics, there appear to be sizeable discrepancies for the contributions from the scalar and tensor operator. I appended a pdf file for further illustration.

flaviovsanalytics.pdf

I am wondering if these discrepancies could arise from a problem with the implementation in Flavio, or if they simply reflect the uncertainties coming along with different form-factor parametrisations? With the potentially related issue #156 regarding the implementation of the SM uncertainties for R(D(*)) having been raised here recently, I assume it might be worth considering the first option.

Many thanks, Tobias

peterstangl commented 3 years ago

The issue with the scalar Wilson coefficients turned out to be a bug and was solved by PR #164.

The differences seen in case of non-zero tensor Wilson coefficients on the other hand seem to just reflects different form factor implementations. In fact, a comparison of flavio and eos shows good agreement: CT

@whereforebound do you agree and can this issue be closed?

whereforebound commented 3 years ago

@peterstangl yes, I agree that these discrepancies are most likely due to different form-factor input. I think the issue can be closed.