Studies of Dimension-Six EFT effects in Vector Boson Scattering

We discuss the implications of dimension-six operators of the Effective Field Theory (EFT) framework in the study of Vector Boson Scattering (VBS) in the $pp \to Z Z j j $ channel. We show that operators of dimension-six should not be neglected in favour of those of dimension-eight. We observe that this process is very sensitive to some of the operators commonly fit using LEP and Higgs data, and that it can be used to improve the bounds on the former. Further we show that other operators than the ones generating anomalous triple and quartic gauge couplings (aTGCs/aQCGs) can have a non-negligible impact on the total and differential rates and their shapes. For this reason, a correct interpretation of the experimental results can only be achieved by including all the relevant, bosonic and fermionic operators; we finally discuss how such an interpretation of experimental measurements can be done.


Introduction
Effective Field Theories (EFT) have become extremely popular in the last few years, , proving to be a robust tool for New Physics (NP) searches and BSM physics studies. In general, one can use the EFT to find the low-energy behaviour of a given UV theory (see for example Refs. [36][37][38]). Alternatively, one can use EFT in an almost model-independent way: as a generalised SM extension that can be used to parametrise small deviations observed on experimental measurements of SM observables. The latter is known as the "bottom-up" approach and will be the one used in this work. The main underlying idea to the bottom-up approach is to add higher dimensional operators to the (dimension-four) Standard Model Lagrangian, in a way that is consistent with the known symmetries: SU (2) × SU (3) × U (1). Some relevant early works in this direction can be found in Refs. [39][40][41][42][43][44] and some interesting reviews on the topic are Refs. [45][46][47]. By adding new higher-dimensional terms to the SM Lagrangian, it is possible to parametrise small deviations from the original SM prediction. If such small deviations are found in the experimental data, it is possible to start mapping the space of "New Physics directions", ruling out some and focusing on others. While few questions remain unanswered in the current picture of fundamental interactions, some seem to be quite far from an answer: How does gravity relate with the other interactions? what are dark matter and dark energy? Other important questions, however, can be tackled at LHC and from an SM perspective. The most important ones concerning the origin of electroweak symmetry breaking (EWSB). The Higgs mechanism [48][49][50][51][52] is a very good description of the EWSB, but some details of the latter are still unknown to us: for example, the fact that the spontaneous symmetry breaking can be realised in a linear as well as a non-linear representation. The answer to this enigma may lay in the gauge couplings, which have only been partially studied at LEP: only some triple gauge couplings have been observed, in a very concrete energy regime and under a set of assumptions regarding the final-state radiation whereas the interactions between four gauge bosons will only be observed at LHC. To address the last question, the detailed study of the Vector Boson Scattering (VBS) process is fundamental. In this process is characterised at tree-level by the exchange of weak gauge bosons between two quarks or a quark and an anti-quark. This means, that it gives us direct access to both triple and quartic gauge couplings. Since this is a purely electroweak process, it has a relatively small cross-section at LHC where the large QCD backgrounds dominate everywhere. However, the family of VBS processes has very particular experimental signatures: two very energetic forward jets with a big rapidity gap between them. In this work we combine these two interesting fields, by applying EFT techniques to the study of VBS at LHC. In particular, we study the effects of different dimension-six (dim = 6 ) operators in the total cross section and differential distributions of the purely electroweak contribution to the process pp → ZZjj, more commonly known as VBS(ZZ). We perform here a numerical study of the aforementioned quantities, a full analytical description will follow in an upcoming publication. The outline of the paper is as follows: In section 2 we introduce our notation and conventions as well and define of the family of VBS processes. In section 3 we summarise the state of the art, both on theoretical and experimental aspects, and we discuss the issue of anomalous couplings. In section 4.1 we compare some of the published results for dimension-six operator fits with our predictions for the cross-section of this process. In sections 4 and 5 we study the impact of different operators of the Warsaw basis on the differential distributions. At first, we focus only on TGC/QGC operators. The full analysis is shown in section 5.1, where we take into account all the Warsaw basis operators, bosonic and fermionic. In section 6 we focus on the main signatures for VBS: dijet observables. Showing that the effects of certain operators (in particular the four-quark ones) can be enhanced on such observables. As anticipated, the family of VBS processes have a relatively small cross section at LHC. The situation for the ZZ final state is particularly dramatic, where the QCD background is very large, with a signal to background ratio of up to 1/20 in some phase space regions. For this reason, a rigorous treatment of the process demands also the study of the EFT effects on the corresponding background, which we perform in section 7. To finalise, in section 8, we discuss a possible strategy for a global analysis including all the dimension-six operators relevant to this process.

SMEFT: notations and conventions
In this work, the bottom-up approach to EFT is used. The SM Lagrangian is extended with higher dimensional operators, consistent with the known SM symmetries. Further, we assume a linear representation for the physical Higgs field, in the form an SU (2) doublet. Such a theory is commonly known as SMEFT: (2.1) At dimension-five, dim = 5, there is only one possible operator, from Ref. [53], which doesn't enter the process studied in this work. At dim = 6 , the complete basis has 59 operators in the flavour universal case and 2499 in the most general one. In this work we will use a parametrisation of the former commonly known as the Warsaw basis, from Ref. [54]. A general method to construct higher dimensional bases using Hilbert series was proposed in Ref. [55]. In the context of VBS, some subsets of dim = 8 operators affecting quartic gauge couplings have been proposed in Refs. [56,57] Other EFT bases There are additional dimension-six bases, other than the Warsaw basis. It is quite common to use the SILH basis, from Ref. [1] in Higgs phenomenology, however it is not optimised for multiboson processes. Instead, there is a VBS-dedicated basis, typically known as the HISZ basis, from Ref. [58].
Parameter Shifts Adding higher dimensional terms to the SM Lagrangian has three consequences: firstly, new vertices appear. For example those with four-fermions. Secondly, the SM vertices get modified with an additional EFT contribution of the form: V SM = a · g + b · g · c i /Λ 2 , where g is the SM coupling and c i is the Wilson coefficient associated with the i th dim = 6 operator. Thirdly, there are shifts on the other SM parameters. Namely the masses, vev, weak mixing angle and gauge fixing parameters. For a detailed discussion on the parameter shifts and gauge fixing in SMEFT see Refs. [6,[59][60][61].
The easiest example to understand parameter-shifts is that of the Higgs field: if we add the Warsaw basis operators to the SM Lagrangian, the Higgs part of the Lagrangian gets modified as: where: and Φ is the Higgs doublet, in the Feynman and unitary gauges respectively. Expanding Equation (2.2), the kinetic term gets modified, as well as the potential (i.e. its minimum gets shifted). To restore the correct vacuum expectation value and canonical normalization of the kinetic term, the corresponding SM parameter and field have to be redefined as: After some simple algebra one finds: As is customary in the literature, we re-define the Wilson coefficients as:c i = v 2 Λ 2 c i . The consequence of the field and parameter redefinitions in Equation(2.5) is that each and every vertex containing the Higgs field will be dependent on the Wilson coefficients {c H , c HD }. This effect is nothing else than a wave-function renormalization in the more classical sense.
The same phenomenon occurs for all the other fields and parameters. The case of the EW sector has to be handled with great care, specially when working on the Feynman gauge, since new linear transitions between gauge and Goldstone bosons appear. The gauge piece of the Lagrangian gets modified as: As a consequence, the gauge fields get shifted in the same fashion as shown for the Higgs field. The neutral gauge boson for example 1 : where θ w is the weak mixing angle. For this reason, when one studies a concrete process, it is important to also take into account these shifts and not only the EFT effects on single vertices. For example, the operator O HB does not directly modify triple or quartic gauge vertices, but it enters the Z field normalization and hence, any vertex containing the former. Additionally to the field shifts, also get shifted the parameters: {M W , M z , θ w , g, G F }. It is well known that these parameters are not all independent. The relation M z = M W cos(θw) , is valid at leading-order in the SM, but has to be corrected at NLO. The same holds for the SMEFT, where this and similar relations need to be re-examined (e.g θ w (g, g ), M H (G F , v)). For this reason, the input parameter set chosen for a calculation will involve different EFT parameters depending on the choice.
Importance of the Input Parameter Set (IPS) It is important to recall that different IPS lead to different predictions, already at tree level. While the "α-scheme": {M z , α, G F } might be convenient at lower energies, or when using experimental measurements of EWPD, for the energy scale of our process (E ≈ 2M z ) the "M W -scheme": {M z , M W , G F } is a better choice. The modifications to the vev are shown in Equation(2.5), and analogously, G F gets modified by the {c , c Hl } operators. The one-loop renormalization of G F in SMEFT has very recently been calculated in Ref. [62].
The Wilson coefficients of the effective theory are not observable quantities, and the M S renormalization scheme is adopted. The SM masses and the electroweak coupling, on the other hand, can be related to experimental quantities and the on-shell renormalization scheme can be adopted.

SMEFT amplitudes and cross-sections
In order to derive physical quantities, it is convenient to look further than the EFT Lagrangian, to S−Matrix elements, since those are the gauge invariant objects that can be projected into experimental quantities. The most general EFT amplitude can be expressed as: The first EFT term corresponds one dim = 6 insertion in the original tree-level diagram, the second term represents two dim = 6 insertions, then one insertion and one loop, two insertions and one loop, then one dim = 8 insertion at tree level. And so on, as illustrated in figure 1. The leading order EFT corresponds to the first term, A (1,1) 6 . The term following that, A (1,2) 6 , with two insertions at tree level, is tricky, since it is of order (1/Λ 2 ) 2 and the EFT expansion is only renormalizable order by order in the expansion parameter 1/Λ. For this reason it is normally not included in the EFT nor NLO-EFT calculations, unless the corresponding counterterms are available. The term A The perturbative EFT expansion, then, grows in different directions: higher dimensional operators or higher loops. If the tree level dim = 8 correction is larger than the NLO-dim = 6 term will depend on the energy of the next new physics scale Λ which, up to now, is unknown. In this work, the leading order for an EFT amplitude is defined as: . Still, some ambiguities appear when squaring the amplitude in (2.8), The term ∼ḡ 2 |A (1,1) 6 | 2 , is commonly called "quadratic EFT" in the literature. It can be used as an estimate for the theoretical uncertainty; it is interesting to study it in depth too, since it is positive definite and hence it can have relevant effects on the differential distributions and on the unitarity bounds of the EFT expansion. This term will be studied in the context of VBS in section 4.2.

Off-shell effects
Higher dimensional operators are always suppressed by a power of the new physics cut-off, Λ. This means that the expansion parameter in the perturbative EFT expansion is E/Λ, where E is the energy scale of the process under study 2 . I.e. near the Z-pole we can think of the EFT expansion in terms of M z 2 /Λ 2 whereas away from the peak the EFT effects are more accurately parametrised by p 2 T /Λ 2 . For this reason, the high-energy regions (tails of the p T distributions) are the ones where the EFT effects are expected to be largest (i.e. E 1 /Λ >> E 2 /Λ, for E 2 on the pole and E 1 on the tail).
The VBS process is defined by two very energetic jets. This means that the p T (j) and m jj distributions are privileged kinematic variables where to expect EFT effects.

VBS: Definition of the process
The family of vector boson scattering processes is very interesting, since it lays at the heart of electroweak symmetry breaking. Some works describing the details of these processes are Refs. [63][64][65][66][67]. Unitarity and gauge-invariance are conserved in this process thanks to a series of cancellations between Feynman diagrams, and fundamentally thanks to the introduction of the Higgs boson, see Ref. [68,69], for these reasons, the VBS represent a set of privileged channels for New Physics studies. For some applications of NP searches to VBS see for example Refs. [70][71][72][73][74][75][76][77][78].
The possible definitions of the VBS process are multiple. Typically, there are substantial differences between the theoretical definitions, in terms of initial and final states, and the experimental definitions, that constrain the phase space of the final states as well. In particular, it is common that the experimental analyses impose certain cuts to try and decouple the vector-boson fusion (VBF) process, where a Higgs boson in the s-channel is produced from the exchange of weak bosons between a quark and an antiquark. This way, the VBF channel is studied in dedicated Higgs analyses, whereas the VBS channels belong to the multiboson analyses.
Regarding such cuts, there are different VBS-regions that are widely accepted, but they obviously lead to different results. In the most recent LHC results, in Ref. [79], the VBS(ZZ) is defined as the purely electroweak component of pp → ZZjj → ¯ ¯ jj, measured in the region defined by the following cuts: One can define a VBS-enriched region, with the additional cuts: which was also used in some parts of that analysis. In this work we chose a compromise between both regions, and apply the cuts: As anticipated, further cuts impose the two Z bosons to be on-shell, to remove the VBF contamination. Nevertheless it is important to keep in mind that the experimental cut defining "on-shell" Z's (M ¯ ∈ [60, 120] GeV ) is not strictly the same thing as the theory definition for the Z pole.
For a matter of clarity, in this work we focus only the process pp → ZZjj before the on-shell decay. The difference between the process pp → ZZjj followed by the on-shell decay: Z → ¯ and the process pp → ¯ ¯ jj in the aforementioned fiducial region, is negligible. In general when applying multivariate-analysis techniques the first definition is preferred, since it populates the phase space in a more effective way.
For a rigorous EFT treatment, the same study should be performed, to make sure that the difference between both options is also negligible in SMEFT. It is clear that new operators will come into play, mainly the ones connecting quarks and leptons in the final state, however intuition and experience tell us that the VBS-cuts will most likely remove the bulk of that contribution.

Anomalous couplings
Anomalous gauge couplings were introduced in Ref. [80], at a time when the EWSB mechanism had not been thoroughly tested yet, and before the Higgs boson was discovered. Such couplings, defined in terms of ad-hoc variations on the Lagrangian parameters, might be good in a first approximation, but present serious theoretical inconsistencies. The main problem being that they violate gauge invariance and unitarity beyond the leading order.
The EFT approach aims to parametrise small deviations from the SM predictions, which are currently being tested with unprecedented precision at LHC. In that regard, a more consistent approach than the anomalous couplings is needed. In particular, one that is consistent at next-to-leading order. The SMEFT approach considered in this work, in terms of dimension-six operators, represents an optimal solution to the anomalous-coupling problem: it can be understood as a SM Lagrangian where all the parameters are anomalous, but in a way consistent with the QFT rules.
Numerous works regarding anomalous gauge couplings in SMEFT can be found in the literature. Some or the earliest studies are Refs. [58,[81][82][83][84], and more current ones can be found for example in Refs. [85][86][87][88][89][90]. In the upcoming sections we will investigate and discuss the differences between allowing the EFT operators only on the weak-couplings (in the spirit of the anomalous-coupling approach) or allowing them to occur anywhere.

SMEFT for triple and quartic gauge couplings
As part of the legacy of the LEP experiment, it is common in the experimental analysis to study triple gauge couplings (TGCs) in multiboson production channels. For some LEP/LEP2 results, see Refs. [91][92][93][94][95] and for some LHC analysis see Refs. [96][97][98]. Quartic gauge couplings (QCGs), on the contrary, were more difficultly accessed by LEP/LEP2 (for example in e + e − → γγνν and e + e − → γγqq were studied in Refs. [99,100]) and not all QGCs where accessible.
This fact makes the QGCs an interesting goal for the LHC experiments, where QGCs are studied in the VBS channels, for example in the analysis [79]. Still, it is important to emphasise that this approach is extremely misleading, since it implies identifying a collection of thousands of Feynman diagrams with a single (off-shell!) vertex, and more importantly, it implies the assumption that triple and quartic gauge couplings are not originated simultaneously through the EWSB. For this reason an analysis in terms of crosssections, differential distributions or (pseudo-)observables should always be preferred.
At the LEP experiment, the s-channel production of Z/W -bosons was very well under control, as well as their decays, since the electroweak radiation could be deconvoluted pretty accurately from the process. This made it posible to treat triple-gauge couplings as pseudoobservables, which could be measured by the experiments. This is not the case any more at the LHC, and hence, one should not aim at "measuring" triple or quartic electroweak couplings.
Pseudo-observables are well-defined theoretical quantities that can be measured in the experiment, a classic example is that of the set of EWPD. The most promising alternative for LHC physics relies on the study of the residues of S-Matrix poles, which are by default gauge-invariant quantities. For some work in this direction, see Refs. [101][102][103][104][105], and the reviews [46,106].
The impact of dim = 6 operators on triple gauge couplings has been studied in Refs. [87,88]. And the set of dim = 8 operators affecting quartic gauge couplings has been studied in depth in Ref. [56,57,107], and other dim = 8 subset, relevant for diboson studies has been shown in Ref. [108]. Similar studies to the one presented here, tackling vector boson scattering in the dim = 8 basis have been presented in Refs. [109,110], as well as works on VBS in the context of the electroweak chiral Lagrangian, in Refs. [111,112]. However there is no study, to our knowledge, addressing dim = 6 effects in the context of VBS.

Subsets of operators and gauge invariance
Each of the operators in the Warsaw basis is independently gauge-invariant and in principle, it is possible in a tree-level study to select a subset of operators without breaking this gauge invariance. This situation however, will not hold beyond tree-level, where different operators enter through the dim = 6 counterterms, and the full basis is needed for UVrenormalization. For a further discussion on the SMEFT renormalization see Refs. [6,7,9,59].
Gauge invariance is also broken if the effects of certain dim = 6 or dim = 8 operators are only included on a certain vertex and not in other vertices or wave function normalizations. Hl , that enter as corrections to G F .

EFT for the Gauge Couplings
To allow a straightforward comparison with the existing literature, in this section we study the impact of a handful of EFT operators. In particular, the operators that directly affect triple and quartic gauge vertices, and are CP-even. 3 For this preliminary study we generated 3 · 10 5 events for the process defined in section 2.3 using a modified version of the SMEFTsim package [104], interfaced with Madgraph5_aMC@NLO [113] via FeynRules [114] and MadAnalysis5 [115]. We studied the impact of each of the three TGC/QGC operators separately, as well as the sum of them, following the definition of leading-order EFT given in section 2.1.
In Figure 2 we see the impact of each of the three operators individually and the sum of them, for four different observables: the invariant mass of the two final Z and that of the two final jets, and the transverse momentum of the leading Z and leading jet. For the numerical values of the coefficients. In this section, we choose the democratic values c W =c HW =c HW B = 0.06, which correspond to c W = c HW = c HW B = 1 with Λ = 1TeV. In section 4.1 we will discuss the case of the available best-fit values.
We observe that the EFT effects on the invariant mass distributions are relatively homogeneous, in particular in the two-jet case. The VBS signature is characterised by two very energetic jets and the high energy phase space is quite well populated. On the contrary, for the ZZ invariant mass we find that at very high energies we reach the limit where the SM production tends to zero, and the EFT effects become sizeable. This effect is of course dominated by the Monte Carlo uncertainty, but it points us to a region that should be studied in detail. Such regions where the SM production becomes negligible are also those where the quadratic-EFT of equation (2.9) will have a dominant role. This will be discussed in section 4.2.
The case transverse momentum distributions are very interesting themselves. We find out that the EFT effects get enhanced on the tails of such distributions, which is something that we would have expected a priori for all four observables, but is not so pronounced as expected for the invariant masses. For both cases: p T (Z 1 ) and p T (j 1 ), we also reach the regime where the SM production is negligible but the EFT effects remain.

Figure 2:
Here we show the impact of the 3 Warsaw basis operators that affect the triple and quartic gauge couplings. We set the valuesc = c v 2 Λ 2 = 0.06 which correspond to c = 1 for Λ = 1TeV. However it is important to recall that one of the main assumptions of the EFT is that there are no new light resonances, and hence in a histogram like this one we are implicitly assuming Λ > 3TeV, while keepingc = c v 2 Λ 2 = 0.06.

Comparison with LEP and Higgs bounds
Several works have appeared in the last years, where SMEFT predictions are compared with LEP data (in Refs. [4,116,117]) and LHC data (the SILH basis in Refs. [118,119] and more recently the Warsaw Basis in Refs. [120][121][122]).
In order to be consistent, a global fit of the full Warsaw basis would be desirable. For this reason it is necessary to include as many measurements as possible, including fiducial cross-sections as well as differential distributions.
In this section we study how do the published best-fit values enter the VBS(ZZ) fiducial cross-section. In Figure 3 we show the signal strength µ = σ EF T σ SM for the central values given by the profile fit of the operators (black, points) and their 95% confidence level bounds (red, error bars) 4 .
For the LEP case, we find that the values for some operators are off by a large amount (≥ 200%) which means that this channel could be an interesting one to constrain such a fit better. This is not surprising since, as precise as it was, LEP is still a "low-energy" experiment 5 , compared with the energies considered here and throughout LHC Run-2. Moreover, the treatment of the radiation in LEP raises some doubts on the applicability of such measurements for EFT searches, as discussed in section I.5 of Ref. [46].
In the LHC case, we find that the fit is more precise and its predictions are compatible with the ones here. The total contribution including its error bar, however, still departs from the SM expected value. Further, it is important to recall that only a subset of EFT operators is included in the former fits, and hence some directions in the phase space have not been tested. The results shown on Figure 3 further show that only one of the fourfermion operators (O , the one entering G F ) is well constrained, whereas other 11 of these operators enter the VBS cross-section. Contrary to the common belief that four-fermion operators are maximally constrained from LEP measurements, we see that, while this might be the case for most of the leptonic operators, it is not the case for four-quark ones.  [117] and on the right from the LHC fit Ref. [121] mentioned in the text. 4 The data are taken from the aforementioned papers. For the LEP fit we took the values corresponding to the MW IPS and to a 1% theoretical uncertainty. For the Higgs fit we took the values in Table 4 of Ref. [121], but it was not possible to clarify some of the details with the authors 5 An important step of any EFT calculation requires the matching of the EFT calculation, at a low energy scale E0 with the "new physics" scale Λ, to account for the running of the EFT coupling.

Linear and Quadratic contributions
Another issue that should be handled with care is the treatment of quadratic dim = 6 terms. This effect was discussed in section 2.1, and it was thoroughly studied in section I.4.7 of Ref. [46] and in Refs. [123,124]. The quadratic contribution to the cross-section is not included in the SMEFT predictions, as a matter of consistency: its perturbative order is higher than the linear term, equivalent to the dim = 8 operators that are not included either. However, this quadratic term is, by definition, always positive, and hence can have a large impact on the behaviour of the final distributions. There are some situations where the quadratic contribution should be carefully studied: • In the regions when the SM prediction becomes very small. In that case the interference term SM × EF T 6 is dominated by the EF T contribution, which is expected to be large. If this is the case, the quadratic term must also be quite large and should be calculated as part of the theoretical uncertainty.
• In the cases where the linear dim = 8 interference term is included (SM × EF T 8 ), since they are both of the same perturbative order, O(Λ − 4).
The latter situation, where dim = 8 operators where compared with dim = 6 quadratic terms, was recently studied in Ref. [125]. In figure 4 we show an example of the impact of including the quadratic contribution in some key distributions, for the canonical value Λ = 1 TeV. Some good news is that, for higher values of the New-Physics scale, Λ, the ratio between quadratic and linear contributions is expected to smoothen. I.e. r = (1/Λ 2 ) (1/Λ) = 1 Λ , hence r 2 << r 1 for Λ 2 >> Λ 1 . Figure 4: Quadratic effects on some kinematic distributions. The red lines represent the linear contribution to the cross section (dim = 6 interference with the SM) and the blue lines represent the previous contribution plus the purely dim = 6 term. We take the canonical valuesc W = 0.06 and Λ = 1 TeV. The quadratic effects are relevant in general on the tails of the distributions, and in particular on the bins where the SM production vanishes. They also play a very important role in the bins where the interference with the SM is negative, since they may restore unitarity. Λ = 1 TeV is the worst-case-scenario, for higher values of Λ the difference between the linear and quadratic terms decreases.

EFT for the full process
In this section, we investigate the effects of including all the Warsaw-basis operators on the computation of the VBS(ZZ) cross-section. Some examples of the Feynman diagrams that contribute to the process are shown in figure 5. In particular we found, numerically, the following expressions for the bosonic contribution to the total cross-section: And for the fermionic contribution: Both of this expressions have been extracted from simple numerical analysis of relatively small Monte Carlo samples and hence are dominated by the MC uncertainty. Which we estimate of the order of 10% for each of the interference terms. The only purpose of displaying them here is to give an impression of the relative sensitivity of this process to the different EFT operators. It is interesting to observe that the bosonic interference is generally positive, while the fermionic one is generally negative. This means that in the case that all the Wilson coefficients would have the same sign, both interferences could extensively cancel, giving rise to a very small SM deviation in the total cross section. For this reason it is fundamental to define observables and regions where the EFT effects are maximised.
To understand the impact of the full dim = 6 basis, we defined different benchmark scenarios where we study the differential distributions for the most interesting VBS observables.

Benchmarks 1 and 2
We consider all the bosonic, CP-even operators in the Warsaw basis, and given the relative contribution of each of them to the total (linear) cross section, given by equation (5.1), we find two possible solutions to the equation: one that gives an O(10%) enhancement and one that gives a O(10%) decrease (negative interference) to the SM total cross section. This is approximately the sensitivity we have to the process in LHC Run-2, and the order of the EW corrections to VBS at very high energies, see Refs. [65,[126][127][128]. In Figure 6 we see the effects of these two benchmarks for four different VBS observables, and in table 1 we give the values used for the numerical simulation. For this part of the study we generated 9 · 10 5 events for each benchmark scenario, using the tools described in section 4. Figure 6: Bosonic benchmarks B1 and B2. The "local" effects on individual observables and bins are very different from the global enhancement or decrease of the cross-section. The observables related with transverse momentum seem to be more discriminating than the ones related to invariant masses. As anticipated, the EFT effects are larger on the tails of the distributions.

Four-fermion operators
In this section, we repeat the same procedure for purely fermionic operators. An interesting feature of four-fermion operators is that they are always generated at tree-level in the UVcompletion, whereas the rest of the dimension-six operators can be generated either at tree or loop level in the UV theory. This means that the effects of four-fermion operators will be often enhanced by a factor of 16π 2 with respect to the rest of the basis. This is particularly relevant for the case of study here, since purely gauge operators, those built of three field strength tensors in the dim = 6 case of four field strength tensors in the dim = 8 basis, are always generated from loops in the UV-completion, and hence suppressed by 16π 2 = 157.914. For the original works on the PTG/LG classification see Refs. [44]. For more recent discussions see Refs. [10,129] In this study, we consider all the fermionic operators. There are 12 operators that dominate the EFT contribution to this process, plus 3 additional ones which are very colour suppressed. Out of these, only O can be constrained from the available fits, since it enters G F in the way described in section 2. The rest of the four-fermion operators remain unconstrained, to our knowledge.

Benchmarks 3, 4 and 5
Equivalently to the previous section, we find the solutions that enhance (B4) or diminish (B3) the SM cross-section by approximately 10%. In particular we chose the values given in table 2. It is important to understand that these becnhmarks are solutions to the equation (5.2), that conspire to enhance or decrease the SM cross-section, but the single entries of Table 2 have no physical meaning by themselves. Additionally we show the benchmark B5 which gives rise to the same total cross section as B3, but has very different kinematics. This study can be seen in figures 7 and 8.

Di-jet Observables
Di-jet observables (m jj , ∆η jj ) are characteristic quantities to any VBS study, since the main signature of such processes is given by their jets. Di-jet data in general represent a great opportunity to constrain four-quark operators at LHC, some attempts in this direction have been presented in Refs. [130,131]. In this section we show the rapidity distributions corresponding to the different set-ups that were studied. This variable seems to be a very good pointer for NP effects in VBS. Figure 9: Top: Bososnic benchmarks (left) and fermionic benchmarks (right). The ∆η variable seems to be a very good flag for new physics effects in VBS, where two rapidity peaks are always observed. Bottom: Effect of the B1 scenario in the background process, discussed in section 7. The effects in the background are more subtle. Still it is important to notice that the number of background events is about one order of magnitude larger per bin, even in the VBS enriched region studied in this work.

EFT in the Background
The main background at LHC for the previously studied process is the QCD induced pp → zzjj. Which has a very large cross-section compared with the VBS one and it is mainly discriminated from the signal thanks to the jet signature of the latter. This background was analysed in LHC in Refs. [96,132].
Moreover, in the available VBS(ZZ) analysis, the S/B discrimination was only achieved by means of a boosted decision tree (BDT) and a Matrix Element (ME) discriminator, described in Refs. [79,133].
Given this situation, it wouldn't be sensible to add the EFT effects in the signal and not in the much larger background. In the following we show a preliminary study of the EFT effects in the background and we discuss which regions and observables are best for the observation of EFT effects.

Characterization of the background
For this study, we generated a sample of events of the QCD component, i.e. O(g 2 s ), to the pp → ZZjj process. The cross-section for this component is much larger than the purely EW one, although this discrepancy can be minimised in the VBS fiducial region, yielding: The sensitivity to the different dim = 6 operators can be extracted numerically in the same way as in equations (5.1)-(5.2), Hl + 0.018c (1) This background is sensitive to less operators than the signal, and it is sensitive to other operators {O HG , O HG } to which the signal was blind. Some examples of Feynman diagrams for the dim = 6 background can be seen in Fig.10. In order to see the effects of the dim = 6 in the background, we generated a sample of background events, including the EFT operators corresponding to the benchmark scenario B1 in table 1, and leaving the new operators equal to zero 6 . The first interesting observation is that this benchmark, that produced a negative interference of 20% in the signal cross-section, gives rise to an enhancement in the background cross section, of approximately ≈ 1%. Given the fact that the background cross-section is 4 times larger than the signal one, eq. (7.1), a few per-cent number of events in the background can be equivalent to an O(20%) number of signal events. This means that the EFT effects in the background should always be well studied and understood. Else, an effect like this: background enhancement and signal reduction, could lead us to miss some interesting EFT effects if we only look at signals and total rates.
In figure 11 we show the corresponding plots for the background study.  : Effect of the B1 scenario in the background process. The effect is in principle very small, however it is important to remark that the number of background events is about one order of magnitude larger per bin, than that for the signal.

General Strategy for EFT studies at LHC
In this work we have shown that at least 11 bosonic and 12 fermionic dimension-six operators enter the VBS(ZZ) signal cross-section. If we try to isolate the effects of each of these operators using a single cross-section measurement, the solutions to this system define a 22 dimensional manifold. In order to further constrain the space of solutions it is necessary to add further experimental measurements.
In this work we have proposed different VBS observables 7 : invariant mass of the two gauge bosons, invariant mass and rapidity of the di-jet system, transverse momentum of the leading gauge boson and jet and the rapidity separation between jets. As we have seen throughout section 5, the EFT scales differently in different regions of phase space (with the high energy "tails" being privileged). This means that a differential measurement with n-bins could be used to add n-equations to the aforementioned 23-variable system, reducing the number of directions in the solution-hyperplane.
Moreover, if we assume factorisation of the EFT effects, by defining the amplitude as in section 2.1, with one insertion per diagram and at tree-level, the fiducial cross-section can be written as, σ EF T = σ SM + A · (F(c 1 ,c 2 , . . . ,c 23 )), (8.1) where F is a linear function of the Wilson coefficientsc i , and A is a global factor modulating the EFT contribution. In that case, the most important task is to find the relationship between the coefficients, F, since the global factor A can be constrained with the measurement of the fiducial cross-section. The best way to do that includes using multivariate analysis (MVA) techniques, like the ones currently used in the experimental VBS searches. For example in the analysis discussed here, or in the analysis of the VBS(WZ) in ATLAS, in Ref. [134], where a 15-variable BDT was used. It is always possible to add different measurements (from LEP, Higgs physics, etc.), in order to reduce the number of unknowns in the system of equations. But this should be dealt-with with great care: a bound set for a certain Wilson coefficient is only valid in the energy regime where the calculation is performed. In order to know the value of this coefficient at a different energy scale, the renormalization group evolution of the coefficient has to be calculated. In that regard, it is safer to include more observables and more bins for signal and background measurements than to mix VBS-signal, Higgs-signal and LEP-signal measurements.
Ideally, gauge boson polarization measurements, not studied in detail here, could also be useful. Such polarization studies happen to be very interesting already at the SM level,as pointed out in recent works like Ref. [77,135]. Further, the study of the parton shower effects on the EFT distributions, might also shed new light on the problem. There are very few references in this direction so far, for example Ref. [122]. Last but not least, the study of the backgrounds, as shown in section 7, is also necessary for a correct determination of the EFT effects.
At the current level of precision of LHC measurements, and the high level of agreement of the former with the SM predictions, it is advisable to perform any EFT analysis with a great deal of accuracy. For this purpose EFT operators cannot be studied on a case-by-case basis, and a global study of the set of dimension-six operators is necessary in the first place. In a second stage, dimension-eight as well as NLO and quadratic dimension-six effects need to be studied in order to improve the EFT theoretical uncertainties.