search

wilson-eft / wilson

A Python package for the running and matching of Wilson coefficients above and below the electroweak scale
https://wilson-eft.github.io
MIT License
26 stars 18 forks source link

running wrong #55

Open jasonaebischerGIT opened 3 years ago

jasonaebischerGIT commented 3 years ago

when running lequ1_2233 = 1/160**2 in the warsaw basis from from 160 GeV to 60 GeV the output of lequ1_2233 is wrong by several orders of magnitude

jackypheno commented 3 years ago

@jasonaebischerGIT

Looks like its an issue related to the size of WCs, for lequ1_2233 I get sensible output for

'wc_in = Wilson({'lequ1_2233':1e-6}, 160, 'SMEFT', 'Warsaw') wc = wc_in.match_run(60, 'SMEFT', 'Warsaw') print(wc.dict['lequ1_2233'])'

(1.0333888075747854e-06+2.115327309804727e-18j)

But the moment we change the initial condition to 1e-5 it gives

(-1.498326492711292e-11-1.39516319299177e-22j)

The problem persists even when using the 'leading_log' approximation.

DavidMStraub commented 3 years ago

I would check what happens before flavour rotations. Perhaps it's not a bug but just an effect of flavour rotations (i.e. that the output basis is not the same as the input basis) like in #63.

jackypheno commented 3 years ago

@DavidMStraub Good idea! Let me check that.

peterstangl commented 2 years ago

@jackypheno did you actually find out if the results was just due to a change of the flavour basis or if there is actually an issue?

peterstangl commented 2 years ago

I changed the label since it is not clear that this is actually a bug.