Precision diboson measurements at hadron colliders

We discuss the measurements of the anomalous triple gauge couplings at Large Hadron Collider focusing on the contribution of the ${\cal O}_{3W}$ and ${\cal O}_{3\tilde W}$ operators. These deviations were known to be particularly hard to measure due to their suppressed interference with the SM amplitudes in the inclusive processes, leading to approximate flat directions in the space of these Wilson coefficients. We present the prospects for the measurements of these interactions at HL-LHC and HE-LHC using exclusive variables sensitive to the interference terms and taking carefully into account effects appearing due to NLO QCD corrections.


Introduction
The Standard Model (SM) is an exceptionally powerful theory in describing most of the observed phenomena in particle physics. The experimental measurements at the Large Hadron Collider (LHC) following the discovery of the Higgs boson [1,2] continue to confirm it to be the valid theory for the larger and larger scales of energy. At the same time we know that the SM has to be just an effective field theory, since it cannot provide explanations to various experimental observations such as the existence of neutrino masses and oscillations, the existence of dark matter and the large baryon asymmetry present in the Universe. Moreover, fine tuning arguments such as the naturalness of the electroweak (EW) scale, the large hierarchy among the Yukawa couplings of the SM fermions and the non observation of charge-parity (CP) violation effects in strong interactions also seem to suggest that an ultraviolet completion of this theory is needed. In particular, in the SM the Higgs mass parameter responsible for EW symmetry breaking is quadratically sensitive to any heavy new physics (NP) contribution. This hints to a relative low energy scale where new dynamic degrees of freedom should be present, unless one is willing to accept a large level of fine tuning in the EW sector. This paradigm has motivated in the last decades a huge experimental effort for the direct search for NP at and beyond the EW scale. The null results so far obtained have however almost ruled out the vanilla new physics models, except in some corners of their parameter space.
In view of this, precision studies of all the possibile deformations from the SM due to new states not directly accessible at current collider energies become a crucial task for the present and future experimental program. The language of the SM Effective Field Theory (SMEFT) provides a well defined organizing principle for characterizing the various deviations from the SM Lagrangian, given the (at least moderate) mass gap existing between the EW scale and the NP scale Λ. As well known in this language the new interactions are expressed as a series of higher dimensional operators so that the effective Lagrangian can be written as L SMEFT = L SM + L 6 + ...  .1) is valid under the assumption of lepton number conservation, so that the first corrections to the SM only appear at the order of dimension six operators, order for which a complete basis has been identified long ago in [3,4]. One of the main goals of the current and High-Luminosity program of the LHC (HL-LHC), as well as of future High-Energy options (HE-LHC), is the precise determination of the c i coefficients. The main objective of this paper is the study for the measurement of two dimension six operators affecting the triple gauge coupling among EW gauge bosons, namely where g and W µν are the SU (2) L gauge coupling constant and the field strength tensor andW µν its dual,W µν = 1 2 αβµν W αβ . It is well known that the measurement of the Wilson coefficient of the two operators of Eq. (1.2) is extremely challenging since the interference between the SM and NP contributions to diboson production in 2 → 2 scattering is suppressed in the high-energy regime as a consequence of certain helicity selection rules [5,6]. This makes it hard to precisely determine the magnitude of the c 3W andc 3W Wilson coefficients, as well as to measure their sign and to differentiate amongst their two different contributions to the scattering amplitudes. Based on the fact that the helicity selection rules of [6] are only valid for 2 → 2 scattering, various observables built out from the decay products of the diboson final states have been recently proposed [7,8]. These observables help to overcome the non interference problem ensuring a larger sensitivity to the Wilson coefficient of the operators of Eq. (1.2).
In this paper we update previous analyses by considering both the pp → W Z and pp → W γ diboson processes with the inclusion of the O 3W operator and carefully treating QCD nextto-leading-order (NLO) effects, which we find to be important since they partially restore the interference between the SM and BSM amplitudes (for the previous studies of the NLO effects in the presence of the higher dimensional operators see [9][10][11]). By closely following existing experimental analyses targeting diboson final states, we provide combined bounds in the c 3W −c 3W plane showing the potentiality of the HL-and HE-LHC options in testing these higher dimensional operators. Interestingly we find that some of the selection cuts which are necessary to suppress reducible QCD background processes automatically lead to a partial restoration of the interference also at LO, an effect extremely relevant for experimental analyses and which was overlooked in the previous literature.
We also compare our findings for the O 3W operator with the limits arising from the non observation of neutron and electron electric dipole moments (EDM). We find that the HL-LHC sensitivity on the CP odd operator becomes stronger than the bounds from neutron EDM but not than the ones from the electron EDM, which are one order of magnitude stronger.
The paper is organized as follows. In Sec. 2 we give a review of the methods recently proposed in the literature to restore the interference of the operators of Eq. (1.2). We then present our analyses for the W Z and W γ processes in Sec. 3 and Sec. 4. Limits from EDMs are discussed in Sec. 5 while prospect for the HE option of the LHC at 27 TeV are shown in Sec. 6. We conclude in Sec. 7.

Interference suppression and its restoration
In this section we will review the main results regarding the interference suppression from the helicity selection rules [6] and the possible strategies to overcome it recently proposed in [7,8]. The reader familiar with the topic can directly skip to the next Section.
Generically, the scattering cross section for any 2 → 2 process in the presence of higher dimensional beyond the SM (BSM) operators can be written as where E is the typical energy of the scattering process, M is the mass of the SM particles and ellipses stand for the smaller terms in the M 2 E 2 expansion The interference terms between the SM and BSM as well as a pure BSM terms are indicated explicitly. In the high energy limit E M the leading contribution comes from the a SM,int,BSM 0 terms in the brackets corresponding to the zero mass limit of the SM particles. In [6] it was shown that a int 0 (the leading contribution to the interference term) is equal to zero for all of the processes containing transversely polarized vector bosons. This effect comes from the fact that the SM and NP amplitudes contain transverse vector bosons in the different helicity eigenstates, for which the interference vanishes. Dramatically, this interference suppressions implies that the high energy measurements of the Wilson coefficients will not benefit from the usual growth of the amplitudes with the energy expected from dimension six operators. This negatively affects the possibilities of high-energy hadron colliders, where the strongest bounds can usually be obtained by exploiting the relative enhancement of the NP contribution compared to the SM one in the high energy distribution tails [12][13][14][15][16][17][18][19][20][21].

Modulation from azimuthal angles: ideal case
For concreteness let us consider the process qq → V T V T , where V = W ± , Z, γ and we will always work in the high energy limit, E m V . In the SM then the only amplitudes that will be generated at leading order in energy are A SM (qq → V T,± V T,∓ ), where the helicities of the final state vector bosons are explicitly indicated. At the same time the dimension six operators in Eq. (1.2) generate only the amplitudes A BSM (qq → V T,± V T,± ). Clearly, there is no interference between the BSM and SM contributions. This is the core of the above mentioned helicity selection rules. However note that at least one of the vector bosons in the final sate is not stable. Hence the physical process is not a 2 → 2 but instead a 2 → 3 or 2 → 4 scattering. For simplicity let us consider the case of qq → W T γ with a leptonically decaying W . The differential cross section can then be schematically written as where the sum runs over the intermediate W polarizations, with the W assumed to be on-shell, and 3 is the Lorentz Invariant Phase Space. In the narrow width approximation the leading contribution to the interference, i.e. the cross term SM × BSM in Eq. (2.2), is given by: where we have ignored the contributions to the longitudinal polarizations in the SM. A simple calculation shows that where φ is the angle spanned by the plane of the W decay products and the W γ scattering plane. As explicitly shown in [8], the phase of the expression M SM can be identified using the optical theorem and its properties under CP transformations. Let's consider an arbitrary amplitude A(a → b). Then the optical theorem (if there are no strong phases, i.e. contributions of nearly on-shell particles) fixes (2.5) At the same time the transformation under CP implies where η CP = 1(−1) for interactions respecting (violating) CP symmetry. By combining Eq. (2.5) and Eq. (2.6) we can infer Applying this result to the qq → W γ process we obtain that Similar arguments can be applied to the case of W Z production. There, since only one pair of the intermediate vector bosons have opposite helicities, the modulation factorizes into a sum of two independent terms and reads The take home message is that by exploiting the modulations of Eq. (2.9) and Eq. (2.10) is possible to increase the precision on the determination of the Wilson coefficients associated with the O 3W and O 3W operators by overcoming the suppression of the interference terms of the cross section, suppression that is recovered with no ambiguity by performing a complete integration over the φ i angles.

Modulation from azimuthal angles: real case
In the previous Section we have discussed an ideal situation, assuming that the azimuthal angles between the plane spanned by the vector bosons decaying products and the scattering plane can be exactly determined. However the azimuthal angle determination suffers from a twofold degeneracy as pointed out in [8]. Let us recall the definitions of the φ angles which can be used experimentally. First we define two normalsn i decay where the index i refers to the first or the second vector boson, l i are the leptons from its decay and ± indicate the lepton helicities. The azimuthal angle φ between the two planes orthogonal to the normals is thus defined as Note that in the case of the Z boson, since its coupling to left-and right-handed charged leptons are approximately equal, we cannot unanbiguosly identify the helicities of the final state leptons. As a consequence the normal vectorn Z decay is defined only up to an overall sign. By using the definition of Eq. (2.12), this translates into an ambiguity (2.13) None of the modulations of the Eq. (2.10) are however affected by this ambiguity, since they are functions of 2φ Z . Now let us look at the azimuthal angle of the leptons from the W boson decay. Differently than for the Z boson, in this case the helicities of the final state leptons are fixed by the pure left-handed nature of the EW interactions. However in this case the azimuthal angle determination suffers from a twofold ambiguity on the determination of the longitudinal momentum of the invisible neutrino, arising from the quadratic equation determining the on-shellness of the W boson. All together for boosted W bosons this leads to the approximated ambiguity This is illustrated in Fig. 1, where we plot the φ W angle constructed assuming a randomly chosen p ν z solution against the same angle where the real, but experimentally unaccessible, value of p ν z has been used. The ambiguity of Eq. (2.14) clearly washes away the sin 2φ W modulations of Eq. (2.9) and Eq. (2.10).
Before concluding this Section a comment is in order. Our definition of the diboson scattering plane of Eq. (2.11) strictly assumes a 2 → 2 scattering process, where the two vector bosons are produced back to back in the center of mass frame of the initial partons. In the case of real radiation emission, as in the case of the presence of initial state radiation jets, the diboson scattering plane has to be defined directly through the momenta of the two vector boson. However in the case of the W Z and W γ processes the determination of this plane will be again affected by the neutrino reconstruction ambiguities. We then decide to use the definition of Eq. (2.11) when building the azimuthal angles of Eq. (2.12) throughout our analysis. This is a good approximation, since only processes with a hard jet emissions, which are kinematically suppressed, can lead to a significant differences between the planes orientations.

pp → W Z process
We begin by studying the fully leptonic pp → W Z process at the LHC. Before doing so we wish to describe to simulation environment which will also be used for the analysis of the fully leptonic pp → W γ process discussed in Sec. 4.

Details on the event simulation
We simulate the hard scattering fully leptonic pp → W Z process via the MadGraph5 aMCNLO platform [22] using the HELatNLO UFO model that have been implemented in the FeynRules package [23] and exported under the UFO format [24] by the authors of [25] 1 . We perform our study at NLO in QCD commenting, when relevant, the differences with respect to the results obtained at leading-order as well as NLO with an extra jet radiation in the matrix-element (hereafter NLO+j), which partially mimics the next-to-next-to-leading-order (NNLO) accuracy. Parton showering and hadronization of partonic events has been performed with PYTHIA8 [26]. Matching and merging between hard-scattering and parton shower have been performed through PYTHIA8 for the NLO case and PYTHIA8 + FxFx algorithm as described in [27] for the NLO+j case. We report in a compact way in Tab. 1 the summary of the tools used for each level of the perturbative expansions. When analyzing the events, jets have been reconstructed via the anti-κ T algorithm [28] with ∆R = 0.4 and a p T threshold of 20 GeV through the MadAnalysis5 package [29] as implemented in MadGraph5 aMCNLO.
While our event generation has been performed at the partonic level, we wish to mimic (at least partially) detector smearing effects when building the angular variables used for our analysis, without performing a dedicated detector simulation for all our event samples. We do so as follows. We choose one event sample and compare, on an event by event basis, the values of the φ Z and φ W variables before and after having applied detector effects, which we have evaluated through the Delphes 3 package [30]. We build the distributions of the ∆φ smear Z,W = φ parton Z,W − φ detector Z,W difference and construct the corresponding probability distribution function. We approximate the latter with a three rectangles shape and dress the parton level values of the azimuthal angles with a ∆φ smear Z,W evaluated with the computed probability. For concreteness we use the following functions 2 [0, π/4] with probability 0.66 [π/4, π/2] with probability 0.34.
(3.1) 1 We thank the authors of [25] for sharing their NLO model files with the addition of the CP-even and CP-odd operators of Eq. (1.2) previous to the publication of their paper. 2 We define the change of the angle due to the smearing to be in the interval [− π 2 , π 2 ] due to the 2φ periodicity of the modulation terms of Eq. (2.10) in the cross section.

Order
Hard-scattering Parton Shower Jet Merging LO MadGraph5 aMCNLO PYTHIA8  In the lower panels we show the NLO/LO and NLO+j/NLO differential cross section ratios.

Comparison of perturbative expansions for the SM
Before proceeding to the study of the sensitivity on the O 3W and O 3W operators we wish to validate our simulation framework against the existing literature for the case of the SM. We consider the pp → W Z processes separately for the two W boson charge signs at LO, NLO and NLO+j order for the LHC with a center of mass energy of 14 TeV. We force the Z boson to decay into a muon pair and the W boson into an electron and the associated neutrino. By applying only a 20 GeV cut on the transverse momenta of all visible leptons we obtain a cross section value at NLO and LO of 37.8 fb and 18.0 fb for the W + case and of 26.7 fb and 11.2 fb for the W − case. The addition of an extra jet in the matrix element increases these value of an extra ∼ 10%. These findings nicely agree with the latest results of [31], computed for √ s = 13 TeV. For the same processes we then compare the differential cross sections in function of the transverse momentum of the charged lepton from the W decay, also reported in [31] for √ s = 8 TeV. By taking into account the parton luminosity rescaling factor between our and their center of mass energy (which is ∼ 2 for the qq scattering of proton's valence quarks for √ŝ = 300 GeV) we find an overall good agreement in the distributions shapes between our LO and NLO results and the ones of [31], thus further validating our simulation framework. Again, we observe that there is a small difference between the NLO and NLO+j calculations. Given the larger computation time needed for the latter simulation, we will present our results only at NLO accuracy in QCD commenting however, where relevant, what could be the effect of the extra real radiation on the processes under consideration.

Sensitivity to the BSM operators
We now turn on the BSM operator O 3W and O 3W defined in Eq. (1.2) and simulate LO and NLO events with the same strategy as for the SM case described in Sec. 3.2, and applying a final combinatorial factor to take into account all possibile final state flavor configurations involving the first two generation lepton flavor. We generate events with only the CP-even or CP-odd operator different from zero, as well as events where both the operators are present, so as to determine the contribution to the cross section due to the interference of the two deformations. We closely follow the ATLAS experimental analysis of [32] and we define our signal region imposing the following sets of cuts: p e T > 20 GeV, p µ T > 15 GeV, |η µ,e | < 2.5, ∆R( , ) > 0.2, ∆R( , j) > 0.4, where = e, µ and the p T threshold for jets is 20 GeV. We further require that the same flavor opposite charge lepton pair reconstruct the Z boson asking |m µ + µ − − m Z | < 20 GeV and we impose a cut of 30 GeV on the W boson transverse mass 3 . We then bin our events with respect to the W Z system transverse mass, which we define as 2) where / p i is the i-th component of the missing transverse momentum of the event. We finally build the φ Z and φ W azimuthal angles as defined in Eq. (2.12) and categorize the events with respect to φ Z and φ W , both defined in the range 0 to π. Now we can proceed to the analysis of the various BSM contributions. Generically the production cross section in the presence of the operators of Eq. (1.2) is given by We firstly compare in Fig. 3 the LO, NLO and NLO+j interference, σ int , (left) and quadratic, σ BSM 1 , (right) terms of the cross section in presence of the CP-even operator O 3W in the angular region φ Z ∈ [ π 4 , 3π 4 ] in function of the m T W Z . We observe that for the pure BSM term the κ-factor between NLO and LO is ∼ 1.3, only mildly growing with the partonic energy of the process, and that the addition of an extra jet in the matrix element only provide a small increase, around 5%,   with respect to the NLO process, similarly to what has been found for the inclusive process in the SM case. On the other side for the interference case, the κ-factor shows a slightly decreasing pattern with the energy of the system, reaching a value of ∼ 2 for m T W Z ∼ 1 TeV. Furthermore helicity selection rules are not applicable at NLO level leading to a mild restoration of the interference effects between the SM and BSM contributions [7,11].
We now proceed in setting the bounds on the c 3W andc 3W Wilson coefficients as follows. We categorize our events with respect to four angular φ Z and two φ W bins, equally spaced in the range 0 to π, and with respect to the W Z system transverse mass, with m T W Z bins between [0,1000] GeV in steps of 100 GeV, [1000,1200] GeV and [1200,1500] GeV. We consider only the SM irreducible W Z background and we impose a global efficiency of 0.6 for reconstructing the final state for all lepton flavor combinations. Then, by assuming a Poissonian distributed statistics, we perform a Bayesian statistical analysis estimating the systematical error through one nuisance parameter (see [7] for more details). We find that the binning in φ W has a marginal impact on the limits determination, which is due to the large smearing on the φ W variable with respect to Z decay products azimuthal angles. Binning our events with respect to φ Z , φ W and m T W Z , we obtain the 95% posterior probability limits 4 on c 3W andc 3W shown in Fig. 4. The limits are shown in function of the maximum m T W Z bin value used for the computation of the bounds and for an integrated luminosity of 3000 fb −1 , i.e. at the end of the high luminosity phase of the LHC, assuming a systematic error of 5%.
We then fix a maximum value of 1500 GeV for the m T W Z bin considered and we show in Fig. 5 the 68% and 95% limits in the c 3W −c 3W plane assuming the SM (left panel) or a signal injection with c 3W =c 3W = 0.4 TeV −2 , again with a systematic uncertainty of 5% and an integrated luminosity of 3000 fb −1 . There the black and red curves correspond to the probability contours with and without the binning in the φ Z and φ W angles and the shaded areas in the left panel correspond to the bounds derived from the non observation of a neutron (dark blue) and electron (light blu) EDM, discussed in Sec. 5.
We observe that the use of the azimuthal variables marginally improves on the limits when   the SM is assumed. This comes out from the combination of three different effects. Firstly, we are considering both the linear and the quadratic term in the EFT expansion, where the latter is not affected from the helicity selection rules cancellation. Secondly the helicity selection rules are violated by QCD NLO effects. Lastly, the imposition of kinematic cuts to select the analysis signal region have also the effect of restoring the interference between the SM and the BSM amplitude. Indeed, we have checked that some of the cuts lead to a partial selection of the azimuthal angles. We postpone the discussion of this effect to Sec. 4.1 when we discuss the W γ process, since the effect is much stronger and the smaller number of final state particles makes it easier to understand the kinematic origin of this behavior. We notice however that the use of the azimuthal angles is crucial in the case of a signal discovery at the LHC. As illustrated in the right panel of Fig. 5  4 pp → W ± γ process We next turn to another process which can be used to test the CP-Even and CP-Odd operators of Eq. (1.2): pp → W ± γ. As for the W Z case, also here we consider a fully leptonic final state which, despite having a smaller branching ratio and the presence of an invisible neutrino, is generally a cleaner channel with respect to the hadronic counterpart. Having validated our simulation framework for the W Z case, we do not perform a comparison of the LO and higher orders samples for the W γ process, and we consider from the beginning of our discussion the event samples generated at NLO accuracy.

Modulation from cuts
Before proceeding with the analysis, we comment here on the partial restoration of the interference between the SM and the BSM amplitudes arising from the imposition of certain kinematic cuts, which we anticipated in Sec. 3.3. Let's consider for example the cut on the W boson transverse mass which is imposed in the experimental analysis [33] and which is defined as where / p T ≈ p ν T . By looking at the dependence of the azimuthal angle φ W with respect to the transverse mass M T W illustrated in the two panels of Fig. 6 we observe that there is a strong correlation between the two variables. In the left panel all events within the detector kinematic acceptance are shown, while in the right panel we additionally impose p γ T > 100 GeV. In both plots we see that a small M T W is in correspondence with a value of 0 or π for φ W . On the other side, for large p γ T , a cut on the W boson transverse mass automatically selects events in the azimuthal bin [π/4, 3π/4]. These two behaviors can easily be understood analytically.
Let's consider first the M T W ∼ 0 case. In this limit the transverse momenta of the decay products of the W boson are parallel: whereâ is a unit vector in the transverse plane. The momenta of the W boson and the charged lepton can be decomposed in a transverse and longitudinal part as whereẑ is a unit vector parallel to the beam line and α W,e and β W,e are two real coefficients. Then Eq. parallel n decay ∝ p ν × p e ∝ p W × p e â ×ẑ (4.4) n scat. ∝ p W ×ẑ â × z so that the azimuthal angle can only take the values of 0 or π. In the high energy regime we can also understand the correlation shown in the right panel of Fig. 6 in the M W T ∼ M W limit. Indeed let us assume that the W boson is strictly on shell. Then the condition M W T = M W leads to (4.5) Let us consider the limit p W T p W z , which is equivalent to requiring p γ T p γ z . This limit in combination with the condition in Eq. (4.5) forces p e,ν T p e,ν z . Hence in this case the normal to the decay plane will be always along theẑ direction, so that the azimuthal angle will take a value equal to π/2. All together we see that a high M W T cut, together with the requirement of a large photon transverse momentum, lead to the automatic selection of a preferred azimuthal angle bin. In the analysis that we describe in the next Section we will bin the events in function of the transverse mass of the W γ system, for analogy with what has been done for the W Z case, where we have used the m T W Z variable of Eq. (3.2). However for a 2 → 2 scattering there is a one to one correlation between the W boson and the photon transverse momenta. Hence, by selecting bins with high M T W γ we automatically select events with high p γ T which, as shown above, lead to the selection of events where φ W ∼ π/2.
It is important to stress that a cut on the W boson transverse mass that we have discussed is imposed in the experimental analysis that we consider [33]. This kinematic selection is used to suppress backgrounds arising from processes without genuine missing transverse momentum, such as the overwhelming QCD γj background where a jet is misidentified as a lepton. Hence this modulation from cuts behavior is always present when performing a real experimental analysis. This is an important effect which has been overlooked in similar studies in the previous literature and that leads to an enhanced sensitivity with respect to what is naively expected. A similar effect also occurs in W Z channel process discussed in Sec. 3.3 and the plots on Fig. 5 reflect this property. However quantitatively we find it to be less important than in the W γ case.

Sensitivity to the BSM operators
We now proceed to the analysis of the W γ final state closely following the 7 TeV CMS results reported in [33], where a measurement of the W γ inclusive cross section has been performed. As a first step we generate fully leptonic W γ events for a center of mass energy of 7 TeV and we apply the same cuts enforced in the considered CMS search. In particular CMS required the presence of a lepton with p T > 35 GeV and |η| < 2.5 and of a photon with p T > 15 GeV and |η| < 2.1 and asked for a separation ∆R( , γ) > 0.7. A cut on M W T > 70 GeV is also applied that, as mentioned, strongly suppresses the backgrounds from processes without genuine missing transverse energy. Then by comparing our NLO predictions with the results of [33] we extract the efficiencies for reconstructing the γ final state, which we quantify to be 0.45 for the electron and 0.7 for the muon. We then use the same efficiency values for the case of the 14 TeV LHC 5 . In order to estimate the detector effects on the determination of the azimuthal angle we follow exactly the same procedure as for the W Z process (see Eq. 3.1) and we find the the following smearing function For our analysis we choose to multiply the rate of the SM W γ process by a factor two, see again [33], to take into account the presence of other background contributions, which we have not simulated. We then bin our events with respect to two angular φ W bins, defined as φ ∈ [π/4, 3π/4] and φ ∈ [0, π/4] ∪ [3π/4, π], and with respect to the W γ system transverse mass defined as with m T W γ bins between [0,1000] GeV in steps of 100 GeV, [1000,1200] GeV and [1200,1500] GeV. We have chosen this variable for the binning in order to make the comparison with the W Z analysis as clear as possible. By adopting this procedure we obtain the results illustrated in Fig. 7 and Fig. 8. In Fig. 7 the bounds are shown in function of the maximum m T W γ bin value used for the computation and for an integrated luminosity of 3000 fb −1 assuming a systematic error of 5%. We can see that the dependence on the maximum m T W γ is different for the CP-even and CP-odd operators. This is due to the fact that we can only restore the interference for the CP-even operator, due to the ambiguity in the W boson decay azimuthal angle, see Eq. (2.14). We have also checked that for the obtained bounds with m T W γ 1 TeV the yields for the CP-even operator are dominated by the interference terms. On the other side at higher energies the quadratic terms start to dominate and the constraints on both the CP-even and CP-odd operators become similar. Then in Fig. 8 we have fixed a maximum value of 1500 GeV for m T W γ and we show the 68% and 95% confidence level limits. There the black and red curve are computed by binning in the φ W angle or inclusive in it respectively and where the left and right hand plot correspond assuming the SM or a BSM signal with c 3W = −c 3W = 0.3 and we again show the bounds from the neutron and electron EDM non observation. As for the W Z case, we can see that for the SM like signal the binning in   the φ W angle practically does not change the results. This is a consequence of modulation from cuts effect described in the previous Section, since the hard cut on the M W T in combination with a high p T of the photon automatically select the value of the W decay azimuthal angle to be close to π/2. Moreover we can see on the right panel of Fig. 8 that even in the case of assuming an injected signal, the results remain the same with and without the azimuthal angle binning unlike in the W Z case. As expected from Eq. (2.9) and Eq. (2.14) the analysis can differentiate the sign of the CP-even interaction c 3W but it is insensitive to the sign of the CP-oddc 3W coupling.

Bounds from EDMs
The CP-odd operatorÕ 3W of Eq. (1.2) gives also a one-loop contribution to the neutron and electron EDMs. Since there are strong constraints from the non observation of EDMs of elementary particles, these null measurements can potentially lead to tight bounds onc 3W . In particular the effective operator generates the EDM operator for a fermion ψ [34,35]. For the case of the neutron we use the form factors of [36] and we obtain By using the latest result reported in the particle data group [37], namely |d n | < 0.3 × 10 −25 e cm at 90% CL, we obtain a limit |λ γ | 0.0023 (5.5) which translates in c 3W TeV 2 0.36 TeV 2 (5.6) which is of the same order with the bounds attainable at the end of the HL-LHC phase from the precision measurements of the W γ and W Z processes. On the other side the experimental limit on the electron EDM is much stronger than the one of the neutron, |d e | < 0.87 × 10 −28 e cm at 90% CL [37]. This leads to a much stronger constraint on the Wilson coefficient of the CP violating triple gauge coupling operator. Namely we obtain We stress however that these bounds can potentially be relaxed in presence of additional new physics contribution affecting the O EDM operator of Eq. (5.2) and cancelling against the one-loop contribution arising fromÕ 3W . We don't discuss this possibility any further, stressing again that the limits arising from the non observation of an electron EDM are potentially more constraining that the ones arising from direct LHC measurements.

pp→WZ LHC 27 TeV
Solid: L=3 ab -1 Dashed: L=15 ab -1 Figure 9: 95% bound on the c 3W andc 3W Wilson coefficients computed with four and two equally spaced angular bins for φ Z and φ W respectively, in function of the largest W Z system transverse mass bin used for the 27 TeV LHC with 3 ab −1 (solid) and 15 ab −1 (dashed) of integrated luminosity. A systematic error of 5% has been assumed.

High Energy LHC
By the end of 2035 the LHC experiments ATLAS and CMS will have collected ∼ 3 ab −1 of integrated luminosity each, ending the HL phase of the CERN machine. Various collider prototypes have been proposed in the recent years for the post LHC era. These include leptonic machines such as ILC and CLIC ideal for performing precision measurements of the Higgs couplings, and hadronic machines, as the FCC-hh, a 100 TeV proton-proton collider, with huge potentiality for the discovery of resonant new physics above the TeV scale, that however requires enormous efforts, among which a new ∼ 100 Km tunnel. Hence in the last years a lot of attention has been given to the possibility of building a higher energy proton collider within the LHC tunnel. Thanks to new techniques with which it would be possible to build 16 T magnets, a flavor of mass energy of 27 TeV can be envisaged. This doubling of energy with respect to the LHC can offer great physics opportunities [38] both for on-shell particle productions, but also for indirect measurements as the ones discussed in this paper.
We then show in this Section the prospects for measuring the c 3W andc 3W Wilson coefficients by applying analyses similar to the ones discussed in Sec. 3 and Sec. 4. We focus on two benchmark of integrated luminosites: 3 ab −1 and 15 ab −1 . The results are shown in Fig. 9-Fig. 12, in complete analogies with the figures of the previous Sections.
For the W Z analysis we can see that the relative improvement from the binning in φ Z and φ W angles increases compared to the 14 TeV analysis, since we are getting closer to the values of the Wilson coefficients when the interference term dominates the cross section. Similar effects hold for the W γ process. The effect of the modulation from cuts becomes less important since for the same values of the m T W γ variable, larger values of the longitudinal momentum are expected at higher collision energies, so that the φ W ∼ π/2 bin selection becomes less strong, see discussion in the Section 4.    Figure 11: 95% bound on the c 3W andc 3W Wilson coefficients computed with two equally spaced angular φ W bins in function of the largest W γ system transverse mass bin used for the 27 TeV LHC with 3 ab −1 (solid) and 15 ab −1 (dashed) of integrated luminosity. A systematic error of 5% has been assumed.

Summary
We have analyzed the diboson production, pp → W Z and pp → W γ at NLO QCD order in the presence of the dimension six operators of Eq. (1.2), paying a particular attention to the effects related to the interference between the SM and BSM contributions. We have found that NLO QCD effects mildly affects the results of the analogous LO analysis [7], since the helicity selections rules do not apply at NLO. For both the pp → W Z and pp → W Z, W γ processes the observables related to the azimuthal angles lead to an enhancement of the interference providing a better sensitivity to the new physics interactions. In order to estimate the LHC possibilities on measuring these interaction we have closely followed available experimental studies of diboson  production [33,32]. Interestingly we have found that some of the kinematic selection cuts needed to suppress the reducible backgrounds in realistic analyses are partially performing an azimuthal angular bin selection. This effect turns out to be particularly important for the pp → W γ processes where the strong cut on the M T W forces the azimuthal angle to be close to π/2, making a further binning in the azimuthal angle φ W less important with respect to what is naively expected. The prospects of the bounds at the HL and HE phases of the LHC are presented. This leads to a sensitivity ∼ 10 −3 on the triple gauge couplings λ Z andλ Z at HL-LHC 6 . The HE phase of the CERN machine can further improve the bounds by factor of ∼ 2−5. These results are summarized in Tab Table 2: Summary of the results for the various channels in terms of the CP-even and CP-odd anomalous triple gauge couplings. Only events with m T W Z,W γ < 1.5 TeV are used. 6 The coefficient λ Z is normalized as µ3 . An analogous definition holds forλ Z .