Wilson coefficients of hermitian operators that should be real still become complex numbers after running. As far as I see, this is only a numerical issue and the imaginary parts are extremely small and basically zero.
However, the fact that the type of these coefficients becomes complex has the effect that some observables in flavio become complex numbers, too. In particular, the observables defined in https://github.com/flav-io/flavio/blob/master/flavio/physics/scattering/ee_ww.py in terms of linear combinations of Wilson coefficients become complex numbers.
Since the ee->WW scattering observables are used in smelli, finally also the global likelihood in smelli becomes a complex number.
Possible solutions:
The first solution is in some sense a bit cleaner, the second one is probably much easier to implement (and maybe faster since there is no need to check which Wilson coefficients should be real).
@DavidMStraub what do you think?
Wilson coefficients of hermitian operators that should be real still become complex numbers after running. As far as I see, this is only a numerical issue and the imaginary parts are extremely small and basically zero.
However, the fact that the type of these coefficients becomes complex has the effect that some observables in flavio become complex numbers, too. In particular, the observables defined in https://github.com/flav-io/flavio/blob/master/flavio/physics/scattering/ee_ww.py in terms of linear combinations of Wilson coefficients become complex numbers.
Since the ee->WW scattering observables are used in smelli, finally also the global likelihood in smelli becomes a complex number.
Possible solutions:
The first solution is in some sense a bit cleaner, the second one is probably much easier to implement (and maybe faster since there is no need to check which Wilson coefficients should be real).
@DavidMStraub what do you think?