EFT triangles in the same-sign WW scattering process at the HL-LHC and HE-LHC

We investigate the Beyond Standard Model discovery potential in the framework of the effective field theory (EFT) for the same-sign WW scattering process in purely leptonic W decay modes at the High-Luminosity and High-Energy phases of the Large Hadron Collider (LHC). The goal of this paper is to examine the applicability of the EFT approach, with one dimension-8 operator varied at a time, to describe a hypothetical new physics signal in the WWWW quartic coupling. In the considered process there is no experimental handle on the WW invariant mass, and it has previously been shown that the discovery potential at 14 TeV is rather slim. In this paper we report the results calculated for a 27 TeV machine and compare them with the discovery potential obtained at 14 TeV. We find that while the respective discovery regions shift to lower values of the Wilson coefficients, the overall discovery potential of this procedure does not get significantly larger with a higher beam energy.


Introduction
The Large Hadron Collider (LHC) has completed data taking for Run II. While a lot of collected data still awaits to be analyzed, no physics beyond the standard model (BSM) has been announced until now. The lack of direct indications for the presence of new physics (NP) makes indirect searches more interesting. The High Luminosity LHC (HL-LHC) upgrade will eventually collect an integrated luminosity of 3 ab −1 of data in pp collisions at a center-of-mass (c.o.m.) energy of 14 TeV, which should maximize the LHC potential to uncover new phenomena. It may however well a e-mail: pkozow@fuw.edu.pl (corresponding author) be that the NP degrees of freedom are at higher masses making it difficult at the LHC to identify experimentally new particles, or new paradigms. These considerations have been driving, in the last few years, intense activity worldwide to assess the future of collider experiments beyond the HL-LHC. Several proposals and studies have been performed. The prospects of pushing the LHC program further with the LHC tunnel and the whole CERN infrastructure, together with future magnet technology, is an exciting possibility that could push the energy up into an unexplored region with the 27 TeV High Energy LHC (HE-LHC), that could collect an integrated luminosity of 15 ab −1 .
Precision measurements provide an important tool to search for heavy BSM dynamics, associated with mass scales beyond the LHC direct energy reach, exploiting the fact that such dynamics can still have an impact on processes at smaller energy, via virtual effects. In this context the well-established framework of effective field theories (EFTs) allows to systematically parameterize BSM effects and elucidate how they modify SM processes. The BSM contributions are effectively parametrized in terms of higher dimension operators O where the superscript n = 6, 8 indicates the dimensionality of the corresponding operator. Following the usual notation, we introduce a set of Wilson coefficients f (n) i , defined as f (6) i = C (6) i which are free parameters since neither C (n) i nor of the full theory are known. Equation (1.1) represents in principle an infinite and model-independent parametrization of possible BSM effects, valid by construction up to the cutoff value in the energy scale of the studied process. The EFT approach is based on the assumption that a limited number of chosen operators contribute to the given process in the studied energy range. In particular, vector boson scattering (VBS) processes are widely recognized as the best laboratory to study the operators which modify only the V V V V quartic couplings. In the SM EFT framework these operators start at dimension-8. The usual ATLAS and CMS procedure in VBS data analyses to date involves testing one such dimension-8 operator at a time. The goal of this paper is to investigate the sensitivity of this approach by identifying for each operator the corresponding reach in the ( f i , ) plane at the HE-LHC.
In a recent paper [1] the physics potential of the single dimension-8 operator EFT approach has been tested on a hypothetical new physics signal observed in the same-sign W W scattering process at the HL-LHC. The analysis was focused on the "gold-plated" purely leptonic W decay modes: where l and l stand for any combination of electrons and muons. In this process the M W W invariant mass cannot be reconstructed experimentally on the event-by-event basis, leading to a restricted space in the ( f i , ) plane (the "EFT triangle") for which the single-operator EFT description of the data is viable. In this note we extend such investigations to the HE-LHC energy and expected luminosity domain. The paper is organized as follows. In Sect. 2 we remind the concept of "EFT triangles". In Sect. 3. we detail our procedure of event simulation and subsequent treatment of generated events. In Sect. 4 we summarize our findings and conclude. Three appendices contain supplementary material. In the first appendix the definitions of dimension-8 operators are recalled. In Appendix B we argue that qualitative features of the full W W scattering process, including off-shell effects, can be inferred from considering the on-shell W W scattering amplitudes and we discuss helicity amplitudes after adding the higher dimension operators. In the third appendix we address the question of what values of BSM couplings can be drawn from the discovery regions.

EFT triangles
Since the truncation of the expansion in Eq. (1.1) introduces model dependence, in Ref. [1] the concept of "EFT models" has been introduced where they are defined by the choice of operators O i . The EFT description is valid up to a cutoff energy at which new states are expected to appear; the cutoff value is unknown a priori. However, in the presence of higher dimension operators the scattering amplitudes grow with energy and eventually break the perturbative unitarity limit M U . The condition < M U ( f i ) defines the upper bound on the range of possible values of as a function of f i .
A specific feature of the process in Eq. (1.3) is that the scale M ≡ M W W , i.e., the invariant mass of the scattered W W bosons, is not experimentally accessible, making it impossible to properly apply the cutoff on the data. Any BSM signal, S, is defined as the deviation from the SM prediction in the differential distributions dσ/dx i of some observable Any collected data sample will in general be a sum of the contributions from M < and M > (unless happens to be out of kinematic reach). To cope with events with M > , different solutions have been advocated in the literature, e.g.: from discarding these events at the level of simulation, to invoking unitarization procedures (usually assuming = M U ) [2,3], to ignoring the cutoff altogether. Any of the above prescriptions is related to additional arbitrariness of choices and therefore affects the physics interpretation of the results. Genuine data interpretation in the EFT language requires its successful description without any additional assumptions as to the nature of BSM physics at the scale above . This is only possible if the bulk of the total observed BSM signal originates indeed from the EFT-controlled range.
The EFT-controlled signal reads:

2)
Here M max is the kinematic limit of the W W invariant mass. Eq. (2.2) defines signal coming uniquely from the "EFT model" in its range of validity and assumes only the SM contribution in the region M > .
The additional contribution from the region above may enhance the signal, but it may also preclude proper description of the data within the EFT. The total BSM signal can be estimated without detailed knowledge of the UV completion from the expected asymptotic behavior for M → ∞, i.e., by assuming that all the helicity amplitudes above remain constant at their respective values they reach at (hence superscript A = const): For every value of , BSM observability imposes some minimum value of f i for which the total BSM signal defined by Eq. (2.3) has enough statistical significance, 5σ in our example. Successful description in the EFT framework imposes some maximum value of f i such that signal estimates computed from Eqs. (2.2) and (2.3) remain statistically consistent, e.g., within 2σ . The "EFT triangle" is the region in the ( f i , ) plane for which a statistically significant BSM signal can be successfully described with a chosen higher dimension operator O i . It is bounded from three sides: • from above by the unitarity limit M U ( f i ), In the HL-LHC case, for all the individual dimension-8 operators that affect the W W W W quartic coupling such triangles were found to be rather narrow or even entirely empty (for O S1 ) [1]. In this paper we extend this analysis to the HE-LHC case, in an attempt to verify if an increased beam energy and integrated luminosity will translate into larger EFT triangles.
Throughout this work we follow the MadGraph convention for the definition of dimension-8 operators (implemented therein via public UFO files [4]), in which the field strength Such conversion factors are equivalent to absorbing the electroweak coupling constants g, each one explicitly factored out for each occurrence of the field stress tensor, in the effective couplings C i . For the reader's convenience, Appendix A lists dimension-8 operators; more details can be found in Ref. [5].

Analysis
In this section we present a generator-level study aimed at finding the EFT triangles for the individual n = 8 operators at the HE-LHC. Event samples of the process pp → j jμ + μ + νν at 27 TeV were generated for each n = 8 operator O i that modifies the W W W W quartic coupling, i = S0, S1 (so called scalar operators), T 0, T 1, T 2 (transverse), and M0, M1, M6, M7 (mixed 1 ones). Generation has been done at LO using MadGraph5_aMC@NLO v5.2.6.2 generator [6], with the appropriate UFO files containing additional vertices involving the desired n = 8 operators. A scan of f i values for each operator was made using the Mad-Graph reweight command, including f i = 0 to represent the SM case. The Pythia package v6.4.1.9 [7] was used for hadronization as well as initial and final state radiation processes. Unitarity limits were determined using the VBFNLO [8,9] calculator v1.4.0, after applying appropriate Wilson coefficients conversion factors. Cross sections at the output of MadGraph were multiplied by a factor 4 to account for all the lepton (electron and/or muon) combinations in the final state (in this work only positively charged leptons are taken into account, although the same analysis can be done for the negative charges). Only signal samples were generated and the SM case was treated as irreducible background in the study of possible BSM effects.
The applied analysis chain to the events generated for HE-LHC is a carbon copy of the one described in detail in Ref. [1], here we only briefly outline the main points, with special emphasis on the differences. Standard VBS selection criteria were applied, namely we require at least two reconstructed jets and exactly two leptons (muons or electrons) satisfying the following conditions: M j j > 500 GeV, GeV and |η l | < 2.5. Although we have not optimized our selection criteria for 27 TeV, one can safely suppose they will not be much different. A detailed optimization will be done with full knowledge of the future detector geometry. Like before, the total BSM signal is estimated according to Eq. (2.3) by suppressing the high-mass tail above the calculated value of . This is achieved by applying an additional weight of the form ( /M W W ) 4 to each generated event in this region. The EFT-controlled signal is calculated according to Eq. (2.2) by replacing the generated high-mass tail with the one expected in the SM (known as "clipping method"). Signal significances are computed as the square root of a χ 2 resulting from a binby-bin comparison of the event yields, with statistical errors such as expected from the data at 3 ab −1 , in the distribution of the most sensitive kinematic variables. Compared to the 14 TeV analysis, the binning of histograms was changed so that in the highest bin the SM prediction normalized to 3/ab is still between 2 and 3 events. In Fig. 1 as an example shown are distributions of four chosen variables: • invariant mass of two leptons M ll , • ratio of transverse momenta of leptons and jets R p T ≡ in the SM and in the case f M1 = 0.2 TeV −4 and = 4.9 TeV. As it was for 14 TeV, we found R p T to be the most sensitive kinematic variable for O S0 and O S1 , and M o1 for the remaining operators.
As the unitarity limit we take always the lower of the two values between on-shell W + W + and on-shell W + W − scattering, calculated from T-matrix diagonalization in the helicity space. Indeed, both processes probe the same quartic coupling and are governed by the same Wilson coefficients, as further explained in Appendix B.3. Figure 2 shows the results for the individual operators S0, S1, T 0, T 1, T 2, M0, M1 and M7, in comparison with results at 14 TeV (for positive f values). Not unexpectedly, all the triangles are shifted to lower f values compared Assumed is √ s = 27 TeV to 14 TeV, the shift being as large as almost an order of magnitude. However, the total area of the triangles does not get significantly larger as we increase the energy. This is because the EFT consistency criterion pushes the effective upper limits on f in a similar manner as does the BSM observability criterion for the lower limits. Overall, the shapes and sizes of all the EFT triangles are remarkably similar for 27 TeV as for 14 TeV, only their respective positions differ. Figure 3 shows the respective results for negative f values of S0, S1, T 0, T 1, T 2, M1 and M7. Here exactly the same observations can be made again. The negative f values of M0 look virtually identical to their positive counterparts, since for these operators the SM-BSM interference term in the total amplitude calculation is practically negligible (see the Appendix for details), and so we do not show them here. There is no triangle at all for S1, for which the overall lower limit for BSM observability is about 1.2 TeV −4 and the upper limit for EFT consistency is 1.4 TeV −4 . Here as well we observe a similar behavior as for 14 TeV.
Both analyses, at 14 and 27 TeV, were done at the generator level. Reducible backgrounds were not simulated, since they are known to be strongly detector dependent (for the different compositions of reducible background as measured at 13 TeV by CMS and ATLAS, see Figs. 2-a and 2-b in Ref. [10] and Fig. 3 in Ref. [11], respectively, where it has been found that 'non-prompt' leptons (from hadron decay and fakes) is the largest reducible background and WZ is small in CMS, while in ATLAS WZ is the largest background and 'nonprompt' is small). Such simulations have to be carried by each experiment for its specific detector performances. For a realistic estimate of the sensitivity limits to new physics effects the reader is referred to literature, e.g., Ref. [12]. Full detector simulation will move the EFT triangles to larger Wilson coefficients, but it cannot make them larger, hence all our conclusions still hold.
For the sake of a convenient comparison between the respective results at two different pp beam energies, in the bulk of this study we have always assumed the same integrated luminosity of 3 ab −1 for both cases. This number is appropriate for the HL-LHC stage, but underestimates the expected statistical power of the HE-LHC. However, it is trivial to recalculate all the results to 15/ab in order to get the true expected discovery reach of the HE-LHC, taking into account its actual expected luminosity. An increase of statistics by a factor 5 will lead to a further shift of all the EFT triangles by a factor close to √ 5, both in the 5σ discovery and the 2σ consistency curves (in fact, somewhat less than that because of non-linear dependence of the BSM signal on the value of the individual Wilson coefficients). It will not significantly change either the shape nor the size of triangles. A comparison of results calculated for the same pp beam energy of 27 TeV and two different integrated luminosities is exemplified for the M1 operator in Fig. 4.
Our simple procedure to suppress the high-mass tails by applying a ( /M W W ) 4 weight to events generated above the scale of works reasonably well in the vicinity of the unitarity limit. In this region it produces a tail falling approximately like 1/M 2 W W , which is the expected asymptotic (i.e., for M W W >> ) behavior of the total cross section after regularization. It nonetheless becomes too strong as we go to << M U , where the total cross section is still dominated by the SM contribution which does not require any further suppression. Moreover, for low values of the tail itself becomes large, leading to large uncertainties due to the details of its modeling. We have discontinued the curves in Figs. 2 and 3 below the values at which we find the method lead to the unphysical result of signal being suppressed below the SM level itself. For this reason the EFT triangles for T 0, T 1 and T 2 do not close. For the remaining cases, however, they are completely contained in the region where our simple method is still viable.
According to Ref. [13], full NLO (EW + QCD) corrections to the SM same-sign WW scattering process at 13 TeV lower the total cross section within a fiducial volume defined by standard VBS cuts by 17%; the relative effect increases in the high invariant mass region. Corresponding corrections at 27 TeV were reported in [14]. They are relatively a few percent larger than in the former case. Therefore, accounting for higher-order corrections is well motivated in further more dedicated studies.
Positivity constraints applied to VBS, derived e.g. in Ref. [15], require that certain linear combinations of Wilson coefficients of dimension-8 operators to be positive definite. Since we choose one operator at a time, the bounds would determine what signs of the Wilson coefficients should be considered for each operator. Therefore accounting for these bounds using our results is straightforward.

Conclusion and outlook
Although an increase of the LHC energy vastly improves the sensitivity to new physics effects in VBS processes, the question of EFT applicability is a different one and cannot be solved by changing the energy. The same-sign W W process with its purely leptonic W decays is often considered "goldplated" due to its relatively good signal to background ratio, but the lack of experimental access to the W W invariant mass poses a severe problem in describing the data in terms of the EFT. Despite reasonable sensitivity to BSM effects, such effects if observed will most likely not be possible to interpret using the data from this process alone and applying the usual framework of testing one dimension-8 operator at a time. This conclusion holds regardless of the actual proton beam energy.
The present results reinforce the former conclusion that future VBS data analysis, both at the LHC experiments as well as future proton-proton colliders, should evolve in the direction of multidimensional fits with many higher dimension operators varied at a time. This in turn may require global simultaneous fits to many processes (including W Z, Z Z and semi-leptonic W V , if not other processes) to help disentangle the correlations between signals originating from different operators. Independently of the goal and results of this paper, one always has to remember that the approach of varying one operator at-a-time is not the most appropriate in any case -such a procedure breaks the model independence of the EFT description and therefore EFT requires ultimately a global approach.
Helpful in disentangling the effects of different operators may be also the polarizations of the outgoing W bosons, as different operator subsets, S, T and M, affect different polarizations. There are new theoretical ideas how to project the total VBS cross sections onto individual polarizations without invoking the rather crude W on-shell approximation [16]. W W polarizations can be extracted from the data by fitting simulated templates of the corresponding polarized distributions. Unfortunately, purely leptonic W W decays do not offer the possibility to reconstruct the W decay angle, which is the only strictly model-independent signature of W polarization. While many other distributions exhibit qualitative differences between the different polarizations, they are usually also strongly model-dependent. Consequently, SM templates cannot be used in the BSM case without the risk of losing sensitivity to the BSM signal. If, however, a set of observables is identified for which sufficiently modelindependent templates for W L W L , W T W T and W T W L can be constructed, it could vastly improve the perspectives of future VBS data analysis in the framework of the EFT. is supported by the Spanish MINECO Project FPA2016-78220-C3-1-P (Fondos FEDER) and by National Science Centre, Poland, the PRE-LUDIUM project under Contract 2018/29/N/ST2/01153. JK was supported in part by the National Science Centre, Poland, the HARMONIA project under contract UMO-2015/18/M/ST2/00518 (2016-2020). JK and MS were partly supported by the COST Action CA15108 and are grateful to all the members of the Action for inspiring discussions. ST is supported by Fermi Research Alliance, LLC under Contract No. De-AC02-07CH11359 with the United States Department of Energy.

Data Availability Statement
This manuscript has no associated data or the data will not be deposited. [Authors' comment: The datasets generated during and/or analysed during the current study are available from the corresponding author on reasonable request.] Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecomm ons.org/licenses/by/4.0/. Funded by SCOAP 3 .

A Dimension 8 operators
The following dimension eight operators contribute to the W W W W vertex, without affecting tri-linear couplings: where is the Higgs doublet field, the covariant derivative D μ ≡ ∂ μ + i g 2 B μ + igW i μ τ i 2 and the field strength tensor

B W W scattering: off-shell versus on-shell
In this Appendix we investigate what can be said about the VBS subprocess in the full pp → j jll ν l ν l reaction from the analysis of the on-shell W W scattering process. We start with the discussion of the W W scattering in full pp process, then identify the helicity amplitudes that dominate the highenergy behavior in the presence of dimension-8 operators and discuss the question of determining the unitarity limits.

B.1 W W scattering in the full pp reaction
In the physical process pp → j jll νν l the W bosons are offshell. Nevertheless, in this subsection we would like to show that qualitative conclusions on the influence of dimension-8 operators on the full process can be drawn from the analysis of on-shell W W scattering. To this end, let us employ the identity [16]: to express the numerator of the off-shell vector boson as a sum over polarization vectors μ λ (k). In the frame in which the spatial component of k μ is in the z direction, k μ = (E, 0, 0, k), the explicit form of each polarization vector reads: where k 2 ≡ k μ k μ . In the on-shell limit k 2 → M 2 W the auxiliary polarization vanishes and 0 approaches the exact on-shell form of longitudinal polarization. With the help of Eq. (B.1) one can then rewrite each of the 4 W propagators in each of the diagram that has VBS topology, as Then the parton-level amplitude qq → qqll v l v l with VBS topology can be decomposed as follows The Therefore in the following subsections we will discuss in detail the high-energy behavior of on-shell W W scattering in the presence of contributions from dimension-8 operators and the unitarity bound.

B.2 The on-shell W W scattering and the helicity amplitudes
Let us consider the elastic on-shell W + W + → W + W + in the presence of BSM part represented by a single dimension-8 operator, as in an "EFT model". The scattering amplitude i M can be written as: where A SM denotes the SM part and A BSM represents the BSM part that depends on the Wilson coefficients f i . For the on-shell W bosons we choose to work in the helicity basis in which the polarizations are μ i with i = +, −, 0. There are in total 3 4 = 81 helicity amplitudes i M(i j → kl) corresponding to helicity configurations (i jkl) in the W W → W W scattering process. The total unpolarized onshell W W cross section can schematically be written as: Since there are orders of magnitude differences concerning contributions of different helicity amplitudes to the total cross section it is convenient, using discrete symmetries P and T and Bose statistics, to divide 81 polarization amplitudes into classes. Amplitudes from the same class yield the same contribution to (polarized) cross sections. Hence, in practice one can consider a reduced number of 13 independent polarization classes, taking into account their multiplicities when computing the cross section. It turns out that only a few helicity configurations contribute non-negligibly at high For the case of the SM the contribution from the saturating helicities to the total unpolarized cross-section is shown in Fig. 5. The four saturating helicity configurations are the only ones whose scattering amplitude is asymptotically constant in energy. The remaining helicity configurations behave asymptotically at most as 1/s, hence their contribution is strongly suppressed at large s = M 2 W W . All on-shell W W cross sections are computed with a 10 • cut in the forward and backward scattering regions (which is explicitly written in the plots).
In the presence of n = 8 operators some of the saturating helicities grow with s, maximally as s 2 . The corresponding case for each "EFT model" studied is shown in Figs. 6 and 7 for f i > 0 and f i < 0, respectively. In fact, for each "EFT model" there is at least one polarization configuration with the asymptotic s 2 energy dependence providing dominant contribution to the unpolarized cross section at M W W = M U . In particular, in the case of "EFT models" with scalar operators (S) only the amplitude with all W bosons polarized longitudinally grows as s 2 . In the case of transverse operators (T ) some amplitudes with all W polarized transversally grow as s 2 , while for the case of mixed operators (M) it happens for amplitudes with two longitudinal and two transverse polarizations. It follows from D μ and W μν building blocks of BSM operators which project mostly on the longitudinal and transverse modes, respectively. It is interesting to notice, however, that for different S, T and M distinct polarization configurations of the outgoing W 's dominate the total cross section at large M W W . Measurement of final state W polarizations would give an insight to the dynamics of their interactions.
Since helicity is an observable for the on-shell W W scattering reaction, different helicity configurations do not interfere among themselves. The total unpolarized elastic onshell W W cross sections as a function of the center of mass W W energy and its dependence on the f i sign is shown in Fig. 8. The sign dependence of the total unpolarized cross section, most visible for T 0, T 2 and also present for T 1, M1, M7, S0, S1, is due to the interference terms in Eq. (B.6). More precisely, the dependence on the sign of the f i is determined by the magnitude of SM-BSM terms relative to the BSM 2 ones in the region E ≤ M U . While there are always BSM 2 terms that asymptotically behave as s 4 / 8 , the earlier growth as ∼ s 2 / 4 of the interference terms is not necessarily visible in each of the "EFT models". If the helicity configurations for which the amplitude depends on energy as s 2 are not among the saturating helicities of the SM, extra suppression factor(s) of v/ << 1 with respect to the opposite case, will be present in the SM-BSM terms. The latter means suppressed sign dependence of the unpolarized cross section, i.e. suppressed interference. It can be inferred from the polarization decomposition plots in Figs. 6, 7 that it is the case for the M0 operator, and indeed in Fig. 8 it is seen that the interference effect is practically invisible for this operator.
Although for the off-shell bosons the helicities are not observable and amplitudes with different helicity configurations may interfere, their interference will be dumped by different structures of fermionic currents to which they are coupled.

B.3 The unitarity bound
The dominating polarization configurations in the total unpolarized cross section can be read from Figs. 6, 7. However, the helicity combination that determines the M U , i.e. that yields strongest unitarity bound, is not necessarily among them. The reason is as follows. The partial wave expansion of helicity amplitude starts with J min = max{|λ 1 − λ 2 |, |λ 3 − λ 4 |}, where λ 1,2 and λ 3,4 correspond initial and final W polarizations, and it is the J = J min partial wave that yields the strongest unitarity limit. It has been checked that for the same-sign W W helicity amplitude that depends on energy as s 2 for the case of M operators J min = 1, while for the S and T operators J min = 0. It would imply then that the unitarity limit for the M operators would be weaker than for S and T , especially if only the J = 0 partial waves were considered. However, the same operators affect both the samesign and opposite-sign W W scattering processes, and both processes should be considered for the determination of the unitarity bounds. In the case of the latter reaction the number of independent helicity configurations is 17 resulting from the fact that one can use all three C, P and T symmetries; however one cannot use the Bose symmetry. In this case for each "EFT model", including the M ones, there exists a helicity configuration that depends on energy as s 2 and has J min = 0. As a result, for the M-type "EFT model" the unitarity limit is considerably stronger as compared to the limit color distinction between the signs: except for M1 and M7, upper cross section curves correspond to f < 0; in S0, S1 (T 0, T 2) stronger unitarity limits correspond to f < 0 ( f > 0) Table 1 The helicity combinations yielding the strongest helicity partial waves unitarity limits for each operator in case of each sign of f i . It is always a J = 0 partial wave that yields the strongest unitarity limits derived from same-sign W W partial wave expansion. This should be kept in mind in particular when using a VBFNLO calculator to determine the unitarity bounds that both sameand opposite-sign W W scattering processes are looked at. The helicity combination yielding the strongest helicity partial wave unitarity limits for each operator are summed up in Table 1.

C Extracting the BSM coupling from the discovery regions found
Assuming that the departure form the SM predictions is indeed observed at the HL/HE-LHC we turn to the question what can be said about the couplings of higher dimension operators that defined the "EFT model". While the probed scale can be read off directly from Figs. 2, 3, the matching between fundamental parameters C i of a deeper BSM physics and the Wilson coefficients f i of the low energy approximation is needed to extract the information about couplings.
Let us start with operators that contain the stress tensor W μν . The W bosons, being fundamental SU(2) gauge bosons, would couple to the to-be-integrated-out BSM states via gauge coupling g. Therefore from the corresponding f i one can factor out g 2 for each W μν . The Naive Dimensional Analysis (NDA) [17] suggests then the following matching The factor 2 follows from the relation TrŴ αβŴμν = 1 2 W i αβ W i μν since in NDA the stress tensor W i μν is used for counting purposes rather than the matrix formŴ μν . Extra factor 2 in case of M6, M7 operators is due to differences in the SU(2) structure, as can be seen from the relation O M7 = 1 2 O M1 + . . . , while 1 16π 2 in front of T operators is a single loop suppression factor suggested by the 4-th power of the electroweak coupling factored out.
The remaining dimensionless factors c i could be combinations of where g * , y * and λ * are some gauge, Yukawa and scalar couplings of the BSM sector of full theory, respectively. In the weakly-coupled theory one would expect that c i are naturally of order 1. Using (C.1) it is straightforward to find the range of c i corresponding to the discovery regions shown in Figs. 2, 3. Table 2 For each "EFT model" characterized by a choice of a single n = 8 operator shown are the ranges of the overall coefficients c i in Eq. C.1 that correspond to the discovery regions found in the 27 TeV study c min -c max 130-770 120-1300 670-2200 23-32 45-133 33-140 c min -c max 110-1500 140-2600 410-4500 48-130 54-270 The numerical values, presented in Tables 2 and 3 for 27  and 14 TeV case, respectively, are found to be roughly consistent. However, instead of being of order 1 they are much larger. It suggests that in case of linearly realized spontaneous breaking of SU (2) × U (1) symmetry our method of probing BSM physics is sensitive only to strong dynamics. Alternatively one could consider relaxed assumptions concerning the symmetry breaking mechanism. For the analysis on "EFT triangles" in the so-called Higgs Effective Field Theory [18], see [19].
We turn now to the discussion of the S0 and S1 operators. We assume these are generated at loop level in the BSM. Otherwise they would come associated with n = 6 operators, which are neglected in our analysis. Then, the NDA suggests Again, c i are some combinations of BSM couplings and naturally expected to be of order 1. If we set c i = 1 in Eq. (C. 3), then for f S0 > 0 we find that that g * ∈ (8.5; 10) and g * ∈ (8.4; 10) in the 14 and 27 TeV case respectively. For f S0 < 0 we find g * ∈ (6.2; 8.9) and g * ∈ (7.3; 8.8) for the 14 and 27 TeV case, respectively. The coupling is large, but interestingly it satisfies g * < 4π . For f S1 the discovery regions are empty for both 14 and 27 TeV cases for both signs.