Within the standard approach of effective field theory of weak interactions for ΔB = 1 transitions, we look for possibly unexpected subtle New Physics effects, here dubbed "flavourful Easter eggs". We perform a Bayesian global fit using the publicly available HEPfit package, taking into account state-of-the-art experimental information concerning these processes, including the suggestive measurements from LHCb of R K and R K * , the latter available only very recently. We parametrise New Physics contributions to b → s transitions in terms of shifts of Wilson coefficients of the electromagnetic dipole and semileptonic operators, assuming CP-conserving effects, but allowing in general for violation of lepton flavour universality. We show how optimistic/conservative hadronic estimates can impact quantitatively the size of New Physics extracted from the fit. With a conservative approach to hadronic uncertainties we find nonzero New Physics contributions to Wilson coefficients at the level of ∼ 3σ , depending on the model chosen. Furthermore, given the interplay between hadronic contributions and New Physics effects in the leptonic vector current, a scenario with nonstandard leptonic axial currents is comparable to the more widely advocated one with New Physics in the leptonic vector current.

In light of the very recent updates on the R K and R K * measurements from the LHCb and Belle collaborations, we systematically explore here imprints of New Physics in b → s + − transitions using the language of effective field theories. We focus on effects that violate Lepton Flavour Universality both in the Weak Effective Theory and in the Standard Model Effective Field Theory. In the Weak Effective Theory we find a preference for scenarios with the simultaneous presence of two operators, a left-handed quark current with vector muon coupling and a right-handed quark current with axial muon coupling, irrespective of the treatment of hadronic uncertainties. In the Standard Model Effective Field Theory we select different scenarios according to the treatment of hadronic effects: while an aggressive estimate of hadronic uncertainties points to the simultaneous presence of two operators, one with left-handed quark and muon couplings and one with left-handed quark and right-handed muon couplings, a more conservative treatment of hadronic matrix elements leaves room for a broader set of scenarios, including the one involving only the purely left-handed operator with muon coupling.

HEPfit is a flexible open-source tool which, given the Standard Model or any of its extensions, allows to (i) fit the model parameters to a given set of experimental observables; (ii) obtain predictions for observables. HEPfit can be used either in Monte Carlo mode, to perform a Bayesian Markov Chain Monte Carlo analysis of a given model, or as a library, to obtain predictions of observables for a given point in the parameter space of the model, allowing HEPfit to be used in any statistical framework. In the present version, around a thousand observables have been implemented in the Standard Model and in several new physics scenarios. In this paper, we describe the general structure of the code as well as models and observables implemented in the current release.