Towards the ultimate differential SMEFT analysis

We obtain SMEFT bounds using an approach that utilises the complete multi-dimensional differential information of a process. This approach is based on the fact that at a given EFT order, the full angular distribution in the most important electroweak processes can be expressed as a sum of a fixed number of basis functions. The coefficients of these basis functions — the so-called angular moments — and their energy dependance, thus form an ideal set of experimental observables that encapsulates the complete multi-dimensional differential information of the process. This approach is generic and the observables constructed allow to avoid blind directions in the SMEFT parameter space. While this method is applicable to many of the important electroweak processes, as a first example we study the pp → V(ℓℓ)h(bb) process (V ≡ Z/W±, ℓℓ ≡ ℓ+ℓ−/ℓ±ν), including QCD NLO effects, differentially. We show that using the full differential data in this way plays a crucial role in simultaneously and maximally constraining the different vertex structures of the Higgs coupling to gauge bosons. In particular, our method yields bounds on the $$ {hV}_{\mu \nu}{V}^{\mu \nu},{hV}_{\mu \nu}{\tilde{V}}^{\mu \nu} $$ hV μν V μν , hV μν V ˜ μν and $$ hVff\left( ff\equiv f\overline{f}/f\overline{f}^{\prime}\right) $$ hVff ff ≡ f f ¯ / f f ¯ ′ couplings, stronger than projected bounds reported in any other process. This matrix-element-based method can provide a transparent alternative to complement machine learning techniques that also aim to disentangle correlations in the SMEFT parameter space.