$W$ boson polarization in vector boson scattering at the LHC

Measuring the scattering of longitudinally-polarized vector bosons represents a fundamental test of ElectroWeak Symmetry Breaking. In addition to the challenges provided by low rates and large backgrounds, there are conceptual issues which need to be clarified for the definition of a suitable signal. Since vector bosons are unstable and can only be observed through their decay products, the polarization states interfere among themselves. Moreover, already at tree level, there are diagrams which cannot be interpreted as production times decay of EW bosons but are necessary for gauge invariance. We discuss a possible way to define a cross section for polarized $W$'s, dropping all non resonant diagrams, and projecting on shell the resonant ones, thus, preserving gauge invariance. In most cases, the sum of polarized distributions reproduces the full results. In the absence of cuts, the ratios of the polarized cross sections to the full one agree with the results of a standard projection on Legendre polynomials. While the latter cannot be employed in a realistic environment, a comparison of the data with the shapes of the angular distributions for polarized vector bosons allows the extraction of the polarization fractions in the presence of selection cuts on the charged leptons.

In particular, it would be useful to isolate the longitudinal component where deviations from the SM are larger and easier to reveal.
The basic tool has been known for a long time: each of the three on shell polarizations results in a specific decay distribution of the charged leptons. Two obstructions, however, hamper progress along this path: • Since the W 's are unstable particles, the decays of the indidual polarizations interfere among themselves, even in the limit of narrow width. These interference contributions cancel exactly only when an integration over the full azimuth of the lepton is performed. Acceptance cuts, however, inhibit collecting data over the full angular range. In addition, cuts affect differently the three polarizations and modify the angular distributions. Notice that the effects discussed here are present in every W production channel, not only in VBS.
• With few exceptions, ElectroWeak boson production processses are described by amplitudes including non resonant diagrams, which cannot be interpreted as production times decay of any vector boson, as shown in fig. 1(c) for the VBS case. These diagrams are essential for gauge invariance and cannot be ignored. For them, separating polarizations is simply unfeasible. Furthermore, as it will be shown later, their contribution is, in general, not negligible.
The polarization fractions in the SM for W + jets processes, without cuts on the charged leptons, have been discussed in ref. [25]. The effects of selection cuts have been studied in ref. [26] for W + jets and a number of other W production mechanisms. The analysis in ref. [26] showed how the decay distributions of the charged leptons get distorted by the presence of cuts and how the simple methods that allow to measure the polarization components when no cut on the charged leptons is imposed, fail when cuts are introduced. The interplay between interference among W polarizations and selection cuts has been also examined in ref. [27].
A number of measurements of the polarization fractions of the W have been performed at the LHC, both by CMS [28] and ATLAS [29], in the W + jets channel. Both collaborations have also studied the polarization of the W 's [30,31] in top-antitop events.
It should be mentioned that the charged lepton decay angles cannot be measured exactly because of the difficulties in reconstructing the center of mass frame of the W . Therefore, in practice, other, directly observable, quantities are studies as proxies. Examples are L P [28], cos θ 2D [29] and R p T [32], which is mostly useful for the W + W + channel.
In this paper, we discuss under which conditions it is possible to define VBS cross sections for polarized W 's, and study their basic properties at the LHC. We also show how to overcome the difficulties related to the presence of acceptance cuts for the charged leptons. We believe these preliminary steps to be essential for this kind of measurement.
We are far from being fully realistic. Nonetheless, our results pave the way to future phenomenological analyses. We limit ourselves to the simple case of W + W − production in VBS when both vector bosons decay leptonically. We also consider only the case in which the two leptons are an e − µ + pair, which avoids the small complication of the Z -2 -  contribution which appears when the two leptons belong to the same family. A cut on the mass of the two leptons would suffice to completely eliminate this channel.
While the leptonic decay of the W 's leads to a cleaner environment, the presence of two neutrinos makes the estimate of the invariant mass of the boson boson system more involved. Some of these difficulties might be alleviated studying the semileptonic channel instead of the fully leptonic one. However, in the semileptonic case the QCD background is larger and it is difficult to separate W W from W Z production. In addition, one of the two vector bosons needs to be identified from its hadronic decay.
The present study could be easily extended to include NLO QCD corrections [33,34] since they do not modify the W decay. EW corrections, which have been recently calculated [35] for same sign W 's, potentially mix the production and decay part of the amplitudes and will require an additional effort.
We are confident that if separation between the different polarizations can be achieved in Monte Carlo simulations it will also be within reach of the experiments, given sufficient luminosity.
The structure of the paper is the following: in the next section we recall how the polarizations of the W 's enter the amplitudes when the decay is taken into account exactly and how interferences between polarizations arise. Next, we present the approximations which we propose in order to separate the different polarizations and show how they reproduce the full result in the absence of cuts on the leptons. In particular we discuss how the complete differential distributions are reproduced, in most cases, by the sum of singly polarized distributions with the exception of those variables, like the charged lepton transverse momentum, which constrain the available angular range for the decay. In sect. 7 we introduce acceptance leptonic cuts and discuss how they spoil the cancellation of the interference contributions and modify the simple form of the decay angular distribution. Finally, in sect. 8 we show that the shapes of the angular decay distributions are sufficiently -3 -universal to allow an almost model independent measurement of the polarization fractions. We extract the polarization components in the Higgsless model and in one instance of a Singlet extension of the SM, fitting the full distribution with a sum of singly polarized SM shapes.
2 W boson polarization and angular distribution of its decay products As already mentioned, a vector boson production tree level reaction, in general, receives contribution from different classes of diagrams, both resonant and non resonant. If we concentrate on hadronic processes at O(α 2 EM α n S ), involving a single external W + which decays leptonically, each diagram includes one ElectroWeak propagator of timelike momentum. The amplitude can be written, in the Unitary Gauge, as The polarization tensor can be expressed in terms of four polarization vectors [36] − In a frame in which the off shell W boson propagates along the z-axis, with momentum κ, energy E and invariant mass Q 2 = √ E 2 − κ 2 , the polarizations read: The longitudinal and transverse polarizations obey the standard constraints ε i · k = 0, On shell, the auxiliary polarization is zero and the longitudinal polarization reduces to the usual expression: ε µ 0 = (κ, 0, 0, E)/M W . The most general case, in which the W propagates along a generic direction, can easily be obtained by a rotation.
The decay amplitudes of the W , depend on its polarization. In the rest frame of the ν pair, they are: where (θ, φ) are the charged lepton angles relative to the boson direction in the laboratory frame. The decay amplitude for the auxiliary polarization is zero, for massless leptons, because ε µ A is proportional to the four-momentum of the virtual boson. Hence, each physical polarization is uniquely associated with a specific angular distribution of the charged lepton, even when the W boson is off mass shell.
Defining a polarized production amplitude the full amplitude can be written as: where M F λ is the full amplitude with a single polarization for the intermediate W . Notice that in each M F λ all correlations between production and decay are exact. The squared amplitude becomes: The interference terms in eq.(2.9) are not, in general, zero. They cancel only when the squared amplitude is integrated over the full range of the angle φ, or, equivalently, when the charged lepton can be observed for any value of φ. This remains true in the Narrow Width Approximation in which 1/((k 2 − M 2 ) 2 + Γ 2 w M 2 ) is replaced by π δ(k 2 − M 2 )/(ΓM ). With this substitution, the integration over the invariant mass of the intermediate state becomes trivial, but the angular integration is unaffected.
If we denote by dσ(θ, φ, X)/d Lips the fully differential cross section, where θ, φ are the W decay variables in the boson rest frame and X stands for all additional phase space variables, by dσ(θ, X)/d cos θ its integral over φ, dσ(θ, X) d cos θ = dφ dσ(θ, φ, X) dLips (2.10) and by σ(X) the integral of dσ(θ, X)/d cos θ over cos θ σ(X) = d cos θ dσ(θ, X) d cos θ (2.11) one can write, using eq.(2.5) and eq.(2.6), where the upper sign is for W + and the lower sign for W − . In general the f i depend on the variables X which are not integrated over. The normalizations are chosen so that 13) and f L , f 0 and f R represent the left, longitudinal and right polarization fractions, respectively. The steps which lead to eq.(2.12) can be repeated even when a partial or complete integration over the X variables is performed. If eq.(2.12) holds, the polarized components can be extracted from the differential angular distribution by a projection on the first three Legendre polinomials: This procedure is completely equivalent to the method used in refs. [25,26] based on the calculation of the first few moments of the angular distribution.
The polarization fractions f 0 , f L , f R can be obtained as: where the upper/lower sign refers to the W + /W − . We note that the sum f 0 +f R +f L = 2α 0 is bound to be one. A word of caution is necessary when acceptance cuts are imposed on the charged leptons, as is unavoidable in practice. While the cancellation of the interference terms in eq.(2.9) is a necessary condition for the validity of eq.(2.12), this is by no means sufficient. A generic cut, think of the lepton transverse momentum, will depend on both the angular variables and the variables X. Integrating over X, in the presence of cuts, results in different theta dependence. As a consequence, the measured lepton decay distribution is not described any more by the simple formula eq.(2.12) and the polarization fractions cannot be computed as in eq.(2.15).
In order to separate the polarized components in the data, it is necessary to compute the individual amplitudes M F λ in eq.(2.8), which requires making the substitution in the W propagator. For the present analysis, this possibility has been introduced in PHANTOM [37].

Separating the resonant contribution: On Shell Projection
After our discussion of W boson polarization in processes with a single W , we now turn to reactions which contain non resonant contributions, in particular to those in which two lepton pairs, e −ν e and µ + ν µ , are produced. For the set of diagrams in which each leptonic line is connected to a single intermediate W , as in figs. 1(a), 1(b), which we will call double resonant or just resonant for short, one can proceed as in the previous section. However, -6 -there are many diagrams, like the one shown in fig. 1(c), which cannot be expanded in a similar fashion. The only way to proceed, in order to define amplitudes with definite W polarization, is to devise an approximation to the full result that only involves double resonant diagrams.
There are, obviously, several conceivable approximations. The simplest one is simply to drop all non resonant diagrams, possibly restricting the mass of the lepton pair in the decay to lie close to the W mass. Since this procedure clearly violates gauge invariance, one expects it would produce distinctly wrong cross sections, at least in particular regions of phase space. For this reason we have chosen an On Shell Projection (OSP) method, similar to the procedure employed for the calculation of EW radiative corrections to W + W − production in refs. [38][39][40][41]. This can be realized as a completely gauge invariant approximation.
It consists in substituting the momentum of the ν pair in the calculation of M P λ and M D λ in eq.(2.8) with a momentum on the W mass shell, while the denominator in the W propagator is left untouched. However, this projection is not uniquely defined. In order to have an unambiguous prescription we have chosen to conserve: 1. the total four-momentum of the W W system (thus, also M W W is conserved); 2. the direction of the two W bosons in the W W center of mass frame; 3. the angles of each charged lepton, in the corresponding W center of mass frame, relative to the boson direction in the lab.
Since the original W momenta are typically only slightly off shell, the modification of the kinematics is expected, in most cases, to be minimal. The modified momenta affect only the calculation of the weight of the event. In the LHA event file [42] all particles are assigned the original, unprojected momenta. This procedure can only be applied for M 2 2ν > 2M W . In the following we will refer to the total invariant mass of the four leptons as M W W , for simplicity.
The OSP requires a further adjustment in the computation of the projected amplitudes. In the full calculation, the presence of unstable intermediate, timelike, W and Z bosons, forces the introduction of an imaginary part, tipically −iΓM in tree level processes, into their propagators. This corresponds to the partial resummation of a particular class of higher order contributions, effectively mixing different perturbative orders, and has the collateral effect of complicating the issue of gauge invariance [43][44][45].
When massive vector bosons appear only as virtual states, a simple and effective way of preserving gauge invariance is the Complex Mass Scheme (CMSc) [40,46]. In the CMSc, the squared masses in all massive vector bosons propagators become complex, M 2 → M 2 − iΓM , and the cosine of the Weinberg mixing angle is defined in terms of complex masses as cos 2 , becoming itself complex. The CMSc, however, is not gauge invariant for amplitudes with external Weak gauge bosons, which are implicitely regarded as stable. The simplest process which exemplifies these issues is e + ν e → W + γ. If the polarization vector of the photon is substituted by its momentum, the amplitude should become zero, because of the electromagnetic Ward identity. It is easily verified that this happens only if Γ W = 0. In the end, gauge invariance forces all widths of intermediate vector bosons in the OSP amplitudes to be set to zero. As a consequence all weak bosons whose momentum can get on their mass shell must be projected on shell simultaneously.
In Appendix A we discuss, in a more quantitative and detailed way, the pitfalls of the non fully gauge invariant approximations mentioned in this section. In the Appendix we focus on the high transverse momentum, high mass region where gauge violating effects are expected to be magnified, particularly when longitudinally polarized vector bosons are involved.
All results, in the following, which are labeled as full or complete have been obtained in the CMSc. All OSP results have been obtained projecting on mass shell all resonant vector bosons, setting to zero all widths in non resonant propagators and using real couplings.

Setup of the simulations
We have studied pp → jje −ν e µ + ν µ at parton level. All events have been generated with PHANTOM [37], using the NNPDF30 lo as 0130 PDF set [47] with scale Q = M W W . We consider only ElectroWeak processes at O(α 6 EM ). We neglect O(α 4 EM α 2 S ) processes and exclude the top-antitop background assuming perfect b quark veto.
The OSP requires the mass of the four lepton system to be larger than twice the mass of the W , thus, only events with M W W > 300 GeV have been retained. Since we are mainly interested in VBS at large invariant masses, this is not a limitation.

Validating the approximation in the absence of cuts on the charged leptons
In this section we compare the differential distributions of a number of kinematic variables obtained from the exact matrix element with the incoherent sum of three OSP distributions in which the negatively charged intermediate W boson is polarized while the positively charged one remains unpolarized. Only doubly resonant diagrams are included for polarized processes. No lepton cut is applied and, as a consequence, the cross section from the e-θ cos  incoherent sum of polarized cross sections coincides with the cross section obtained from the coherent sum. Our aim is: • to demonstrate that the W polarization fractions obtained with the OSP are in agreement with those obtained from a standard expansion in Legendre polynomials.
• to show that the OSP reproduces well the rigorous cross section and distributions.
• to explore the differences between the individual polarized distributions as a tool to separate them in the data.
The full cross section is 1.748 fb while the incoherent sum of singly polarized cross sections is 1.731 fb, which differs from the complete result by about 1%.
In figs. 2, 3 the black curves refers to the rigorous result and the violet ones to the sum of polarized distributions. The individual distributions are shown in red (longitudinal), blue (left) and green (right). We have also computed the coherent sum of the three polarized contributions, which is, in all cases, almost indistinguishable from the full result and, therefore, is not shown.
In fig. 2, on the left, we show the distribution of the decay angle of the electron in the reference frame of the e − ν e pair integrated over the full range M W W > 300 GeV. On the right we present the polarization fractions as a function of the invariant mass of the four leptons, which provides additional information since the distributions in fig. 2(a) are dominated by events with relatively small M W W .
The polarization components obtained by expanding the full angular distribution on Legendre polynomials, as discussed in sect. 2, are shown as lighter shade smooth lines in fig. 2(a) and show that the distributions from the polarized generations, in the absence of cuts on the charged leptons, have the expected functional form, and that the polarization -9 -    The incoherent sum of three OSP distributions agrees very well with the full result for the W W invariant mass, fig. 3(a), the mass of the two tag jets, fig. 3(b), the transverse momentum of the e − ν e pair, fig. 3(c), the rapidity of the e − ν e pair, fig. 3(d) and the mass of the eµ pair, fig. 3(e) .
The agreement is less satisfactory for the transverse momentum of the electron, fig. 3(f). Since, as already mentioned, the coherent sum of the three contributions agrees with the complete result, the discrepancy between the black and the violet lines in fig. 3(f) is due to the interference among the different polarizations. The interference is non zero because forcing the lepton p t to a single bin restricts the angular range of the leptons, spoiling the complete cancellation of the interferences among polarizations.
The distribution, shown in fig. 3(g), of the azimuthal angle of the electron, φ e , which is defined following ref. [25], stands apart. Each of the three singly polarized contributions is isotropic, since the azimuthal angle enters the decay amplitudes only as a phase. The coherent combination of the three amplitudes produces a non trivial modulation of the differential cross section. The incoherent sum of the polarized results reproduces well only the average value of the full distribution.
The singly polarized distributions of the transverse momentum of the W , of its pseudorapidity, of the transverse momentum of the negatively charged lepton and of the invariant mass of the two charged leptons depend significantly on the polarization of the W . Fig. 3(c) shows that longitudinally polarized W 's have a markedly softer p t spectrum in comparison with transversely polarized ones. In fact they dominate for p W t < 50 GeV. Fig. 3(d) shows that, while transversely polarized W 's are predominantly produced at small rapidities, the distribution for longitudinally polarized one presents a dip at zero and peaks at η W − = ± 2.
The p e − t and M ll distributions are harder for left polarized W − 's than for longitudinal polarized ones. The right handed W − 's have the softest spectrum. This behaviour is -11 - 6 Joint polarization fractions for the two W 's The projection described in eqs.(2.14-2.15) can be readily generalized to a simultaneous expansion in products of Legendre polynomials of the two variables cos θ e and cos θ µ . Similarly, the substitution in eq.(2.16) can be performed for each of the two final state W 's. The outcome is shown in fig. 4. For ease of presentation, the right and left handed contributions are summed together in the transverse component W T = W R + W L . The polarization components obtained by expanding, in each bin, the full angular distribution on Legendre polynomials are shown in lighter colors. The darker histograms are obtained integrating the amplitudes squared with definite polarization for each W . The two independent determinations of the joint polarization fractions agree extremely well over the full range in M W W . This implies that the method we propose can be relied on for analyzing double polarized cross sections. Fig. 4 shows that the W + T W − T fraction is always the largest one and dominates at large invariant masses, comprising about 70% of the total cross section. The W + T W − 0 and W + 0 W − T components are essentially equal and almost constant at about 18%. The longitudinal-longitudinal fraction is the smallest one, of the order of a few percent. This implies that measuring the scattering with two longitudinally polarized W 's in the final state, will require determining the polarization of both vector bosons, since a longitudinal W is expected in most cases to be produced in association with a transversely polarized companion.

Leptonic cuts and their effects
In this section we document how the distributions presented in sect. 5 are modified by the introduction of realistic acceptance cuts on the electron which is the decay product of the W − , whose polarization we wish to determine. Our results confirm that the polarization fractions of the W cannot be determined anymore by a projection on the first three Legendre polynomials. We show that, in the presence of standard leptonic cuts, the interference among the polarized amplitudes is small. Therefore, the incoherent sum of the three OSP results approximates fairly well, in most cases, the full distribution. We require: p e t > 20 GeV, |η e | < 2.5. (7.1) The full cross section is 1.411(1) fb, the coherent sum of OSP polarized amplitudes gives 1.401(1) fb, while the incoherent sum of singly polarized cross sections is 1.382(1) fb, which differs from the full result by about 2%.
In fig. 5 and fig. 6 we show the same set of results presented sect. 5. A number of features are worth noticing. The cross section from the coherent sum of polarized amplitudes is very close to the incoherent sum of cross sections, which indicates that the interference among polarizations is generally small. The rigorous cross section is well reproduced by both approximations.
In fig. 5, on the left, we show the distribution of the decay angle of the electron in the reference frame of the e − ν e pair. On the right we present the polarization fractions as a function of the invariant mass of the four leptons.
The full angular distribution is approximated within a few percent, over the full range, by the sum of the unpolarized results. The complete result, shown by the black histogram, however, is not of the form of eq.(2.12) and cannot be described in terms of the three first Legendre polynomials. This becomes clear expanding the full result as in eqs.(2.14-2.15), which yields the blue, green and orange smooth curves in fig. 5. Their sum is the smooth gray curve which fails to describe the correct distribution. As already noticed in ref. [26], There are, however, subtler effects into play, as can be seen from the polarization fractions as a function of M W W in fig. 5(b). One notices that the longitudinal component decreases faster than the right handed one for increasing M W W . The two curves, which are very close at M W W < 400 GeV, separate at large invariant masses, in contrast with the trend they display in fig. 2(b). This is due to the different distribution in W rapidity of the two polarized cross sections. With increasing M W W the average absolute value of the W rapidity with longitudinal polarization increases and a larger number of the corresponding leptons fail the rapidity acceptance cut.
Even in the presence of leptonic cuts, the fraction of longitudinally polarized W − 's is well above 10%.
The incoherent sum of three OSP distributions agrees well with the full result for the W W invariant mass, fig. 6(a), the mass of the two tag jets, fig. 6(b), the transverse momentum of the e − ν e pair, fig. 6(c), the rapidity of the e − ν e pair, fig. 6(d) and the mass of the eµ pair, fig. 6(e) .
As before, the transverse momentum of the electron, fig. 6(f), is affected by interferences, but the effect is not particularly enhanced by the presence of cuts.
The distribution of the azimuthal angle of the electron in fig. 6(g) displays again peculiar features. Each of the three singly polarized contributions develops a dip at φ e = 0 and φ e = ±π. The incoherent combination of the three amplitudes follows more closely the full result than when no cuts are applied but does not reproduce the correct curve.
The individual polarizations are not affected equally by the cuts. Typically, the cross section for right handed W 's is reduced the most, followed by the cross section for longitudinally polarized W 's. Left handed W bosons seem to be the least sensitive to acceptance cuts.

Measuring the polarization fractions
The fact that the full distribution is well described by the incoherent sum of the polarized differential distributions allows for a measurement of the polarization fractions within a single model, even in the presence of cuts on the charged leptons. The measured angular distribution can be fitted to a linear combination of the normalized shapes obtained from -15 -Montecarlo simulations of VBS events with final state W 's of definite polarization. We have verified that this procedure works well in the SM for one polarized vector boson. The results shown in sect. 6 suggest that it should be straightforward to extend the method to double polarized events.
However, it would be inconvenient to have to generate templates for all extensions of the SM, being unknown which specific model is realized in nature, should the SM need extending. A natural question is whether this method can provide a model independent, practical way of extracting the polarization fractions from the data, comparable with the Legendre expansion, which is applicable when the W decay is unrestricted by cuts. The essential condition for this approach, is that the shapes of the polarized distributions are sufficiently universal. On one hand, once a polarized W transverse momentum and rapidity are fixed, its decay is completely specified. Therefore, in any sufficiently small bin in phase space, the angular distribution of the charged lepton will be the same, irrespective of the underlying dynamics. On the other hand, the influence of cuts depends on the transverse momentum and rapidity of the W . Different models yield different distributions of these quantities. Hence, for finite size bins, the details of the averaging over phase space will vary, and some degree of model dependence is to be expected.
In this section we pursue this possibility comparing the normalized distributions produced in the SM with those obtained in different models.
To this purpose, we have studied the Higgsless model obtained from the SM by sending the Higgs mass to infinity. Such model, while not experimentally viable any more, can be viewed as the most extreme modification of the SM in the VBS domain, in which the unitarity violating growth of the scattering of longitudinally polarized vector bosons is completely unmitigated by Higgs exchange.
In fig. 7 we present some results for the Higgsless model. On the left, we show the distribution of the decay angle of the electron in the reference frame of the e − ν e pair. On the right we present the polarization fractions as a function of the invariant mass of the four leptons.  Comparing fig. 7(a) with fig. 5(a) it is immediately evident that the two angular distributions of the charged lepton, in the center of mass of the W , for M W W > 300 GeV, are markedly different. The region around cos θ = 0 is clearly more populated in the Higgsless model as a consequence of a larger fraction of longitudinally polarized bosons. This is confirmed by the comparison betwen polarization fractions on the right hand side of the two figures. In fig. 7(b) the growth of the longitudinal component is quite prominent. Notice, however, that, even in this extreme model, the left handed fraction is the largest one for all W W invariant masses up to about 1800 GeV.
We have also examined a Z 2 -symmetric Singlet extension of the SM with an additional heavy scalar [48][49][50][51][52][53][54][55][56][57][58][59][60][61][62]. In this model the couplings of the two Higgses to SM particles are proportional to the SM couplings of the Higgs, multiplied by universal factors: where xx represents a pair of SM fermions or vectors, and α is the mixing angle. As a consequence, since only small values of sin α are allowed [63][64][65], the width of the heavy Higgs is much smalller than the width of a SM Higgs of the same mass. We have taken sin α = 0.2, m H = 600 GeV and tan β = 0.3, where tan β is the ratio of the vacuum expectation values of the two neutral scalar fields, which yields Γ H = 6.45 GeV.
With the exception of the mass window in the vicinity of the heavy Higgs resonance all results for the Singlet model follow closely those of the SM.
We have checked that, for each W polarization, the shape of charged lepton angular distributions, 1/σ ·dσ/d cos θ, in both the singlet and the noHiggs models, are in reasonable agreement with the corresponding SM shape in all 100 GeV intervals in the W W invariant mass, even though the normalizations can be quite different. We have also verified the agreement between the sum of singly polarized distributions and the complete result in both cases.
The results are shown in fig. 8 for the Higgsless model, and in fig. 9 for the Singlet one. In both figures the black histogram shows the result of the exact calculation. The lighter red, green and blue curves are the singly polarized distributions in the Higgsless and Singlet models. The darker red, green and blue curves show the results for the longitudinal, right and left handed polarizations, respectively, of the fit of the full results using the SM templates. The longitudinal and right handed components have been multiplied by three and two, respectively, to improve the overall readability of the plots. Fig. 8(a) and fig. 9(a) show that, when a large range of invariant masses is taken into account, the shapes of the polarized components are indeed very similar in the three models. This in not surprising since the distributions are dominated by low invariant mass events, for which the differences between the various models are expected to be small.   fig. 8. p e t > 20 GeV, |η e | < 2.5. The longitudinal and right handed components have been multiplied by three and two, respectively. Fig. 8(b), fig. 8(c) and fig. 9(b)  and Singlet model appears slightly shifted to large values of cos θ in comparison with the SM distribution. The effect seems to increase for increasing M W W , in fact, it is more pronounced in fig. 8(b), which contains events of higher average mass, than in fig. 8(c). Despite these differences, the fit of both models using SM templates is quite good. We have fitted the full result for each model, by means of a χ 2 minimization, with a linear combination of four normalized functions extracted from the SM distributions: the shapes for the longitudinal, right and left handed components together with the shape of the SM interference. The latter is defined as the normalized difference between the unpolarized, complete, distribution and the sum of the three singly polarized ones.
In tab. 1, we show a selection of results. The rows labeled SM, no Higgs and Singlet show the percent ratio of the singly polarized cross sections to the full, unpolarized cross section for each model. The rows labeled Fit give the coefficients, in percent, of the four SM shapes obtained from the fit of the angular distribution of the electron in the corresponding model.
When the fit is extended to the full range of invariant masses, M W W > 300 GeV, the polarization fractions returned by the fit reproduce the results obtained from the ratio of the OSP polarized cross sections to the complete one.
For the Higgsless model, we have investigated the agreement between these two determinations of the polarization fractions in each 100 GeV interval in the W W invariant mass. The accord is quite good in each bin, even though the normalizations become progressively -19 -different as the invariant mass grows.
The M W W > 1000 GeV range, where the longitudinal polarization fractions in the Higgsless model is more than twice the SM one, gives the poorest agreement. However, the discrepancy between the polarization fractions from the fit and those from the ratio of cross sections is of the order of 10% at most.
The fitted polarization fractions for the Singlet model in the heavy Higgs peak region, 580 GeV < M W W < 620 GeV, are in very good agreement with those estimated from the ratio of the singly polarized cross sections to the full result in that bin.
The results presented in this section suggest that it is possible to extract the polarization fractions of the W in an almost model-independent way, exploiting the similarity of the shapes of polarized cos θ e distributions for different underlying theories, in most of kinematic regions.
The angular distributions we have discussed, are difficult to measure, particularly when both W 's decay leptonically. However, the same approach can be applied to other kinematic variables whose distributions discriminate among the different polarizations, as for instance the invariant mass of the charged lepton pair and the transverse momentum of the charged lepton, as shown in fig. 6. In order to enhance the accuracy of the obtained results, more refined multivariate fit methods could be investigated.

Conclusions
In this paper we have investigated the possibility of defining cross sections for processes with polarized W bosons, including off shell and non resonant effects.
We have proposed a method which is based on the observation that the set of doubly resonant diagrams, suitably projected on shell to preserve gauge invariance, approximates well the full result and allows an expansion in terms of amplitudes in which each final state W has a definite polarization. This procedure agrees with the standard approach based on Legendre polynomials in the absence of cuts on the decay leptons. When acceptance cuts are imposed on the leptons, and the Legendre polynomials procedure fails, it is possible to extract the polarization fractions using singly polarized SM Monte Carlo templates. These same templates can be used to measure the polarization of the W with reasonable accuracy even if new physics is present. If, however, the comparison between the full calculation and the naive approximation is pushed to the high invariant mass, high p t , region, the latter fails to reproduce the full result. This is clearly demonstrated in fig. 11 where the transverse momentum of longitudinally polarized W − is shown for events with M W W > 1000 GeV. While the OSP prediction agrees with the full result over the whole range, the curve obtained simply requiring |M ν − M W | < 30 GeV overshoots the true result by a large factor for transverse momenta above 300 GeV.
We now move to the second example, namely computing the OSP amplitudes in the CMSc scheme, instead of setting the weak bosons widths to zero and using real couplings. Typically, for an amplitude without internal vector bosons which can go on mass shell, the difference between the CMSc and setting all widths to zero is small, of order O(Γ/M ) = O(α). However, when the energy of the vector bosons grows large, these differences can become significant, because any deviation from a fully gauge invariant treatment of the boson widths spoils the cancellation of terms which are proportional to the energy.  In order to examine quantitatively the relevance of these gauge violating effects in VBS, we study the single process uu → uu e −ν e µ + ν µ , for large invariant masses of the four lepton system, M W W > 1000 GeV, large transverse momentum of the e −ν e pair, p W − t > 800 GeV. The mass of the pair is required to be close to the W mass, |M v − M W | < 0.2 GeV. No cuts on the charged leptons are applied.
A complete calculation in the Complex Mass Scheme gives a cross section of 0.641(3) · 10 −8 pb, which we regard as the true result. Restricting the calculation to the doubly resonant diagrams and projecting them on mass shell, as discussed in sect. 3, using complex masses and couplings, yields 0.754(2) · 10 −8 pb, which is not compatible with the correct value. However, when all masses and couplings are taken to be real, the on shell projected result is 0.650(2) · 10 −8 pb, in excellent agreement with the cross section obtained from the full calculation. This proves, in our opinion, that failing to maintain gauge invariance could indeed lead to significantly wrong results. Obviously, it does not prove that any gauge invariant calculation would reproduce the correct answer. In particular, the correctness and usefulness of the projection procedure needs to be verified in each case.
Gauge effects are particularly strong for longitudinally polarized vector bosons. Therefore, we have also computed the on shell projected cross sections for a longitudinally polarized W − : using complex masses and couplings we obtain 1.282(5) · 10 −9 pb, while, taking them to be real gives 0.202(1) · 10 −9 pb, which differ by a factor of about six. Notice that the difference between these two cross sections is equal, within errors, to the difference between the corresponding unpolarized results, suggesting that the discrepancy between the results for longitudinally polarized W − accounts for the totality of the difference in the unpolarized case.
Since no cuts on the charged lepton is imposed, the cross section for a longitudinally polarized W boson can also be extracted from the full distribution, projecting on Legendre polynomials. This yields 0.211(5) · 10 −9 pb, confirming the OSP result with all internal widths set to zero.