W-boson production in TMD factorization

At hadron colliders, the differential cross section for W production can be factorized and it is sensitive transverse momentum dependent distributions (TMD) for low boson transverse momentum. While, often, the corresponding non-perturbative QCD contributions are extrapolated from Z boson production, here we use an existing extraction (based on the code Artemide) of TMD which includes data coming from Drell–Yan and semi-inclusive deep inelastic scattering, to provide checks and predictions for the W case. Including fiducial cuts with different configurations and kinematical power corrections, we consider transverse momentum dependent cross sections within several intervals of the vector boson transverse mass. We perform the same study for the \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$p_T^{W^-}/p_T^{W^+}$$\end{document}pTW-/pTW+ and \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$p_T^Z/p_T^W$$\end{document}pTZ/pTW distributions. We compare our predictions with recent extractions of these quantities at ATLAS and CMS and results from TeVatron. The results encourage a broader experimental and phenomenological work, and a deeper study of TMD for the W case.


Introduction
Vector boson production is very relevant at hadron colliders and its measurement has achieved an increasing precision in the latest years . The experimental results have allowed the extraction of important Standard Model quantities like collinear parton densities and the mass of the W -boson [26][27][28][29]. The precise measurements of W cross section and similar observables need a good determination of QCD nonperturbative inputs which can be partially extracted from neutral boson mediated processes. Because we are interested in transverse momentum dependent observables it is mandatory to consider the extraction of transverse momentum dependent a e-mail: dangut01@ucm.es b e-mail: sergiol95@univie.ac.at c e-mail: ignazios@ucm.es (corresponding author) parton distribution functions (TMD) that have used the data of the neutral bosons at low transverse momentum [30][31][32][33][34][35].
In [36] it is pointed out that the low-transverse momentum spectrum for Z -boson production at LHC cannot be fully understood without the introduction of TMD, and further study is considered as necessary (see also [37]). The Wproduction is a natural test for the TMD factorization, and an explicit evaluation of the cross section within this approach is yet missing in the literature to our knowledge. On the other side there exist codes and predictions  for the W spectrum that include non-perturbative TMD effects only in a limited way.
The W -spectrum is also interesting in order to establish some properties of the TMD like their flavor dependence. In fact, out of all the extractions cited above, only [35] can be sensitive to some non-trivial flavor structure of TMD because it includes also data from semi-inclusive deep inelastic scattering (SIDIS) from Hermes [63] and Compass experiments [64], despite the fact that it analyzes only processes mediated by neutral vector bosons. Some preliminary work in this sense has been performed in [65] and later in [66] where the authors conclude that LHC data on W -boson production must be sensitive to the non-trivial flavor structure of TMD. This statement can partially be tested in the present work. We also consider of major importance to well establish the flavor dependence of the TMD before the starting of the Electron Ion Collider (EIC).
The results that we present want to provide a set of predictions for W -production which include the latest information from TMD studies. We revise the kinematics of W -boson transverse momentum distribution in order to include kinematical power suppressed terms as in [35] and to implement the fiducial cuts typical of these processes. We pay particular attention to errors in our predictions as coming from scale variations, PDF, TMD parameterization whose present knowledge is here described in detail. We recall that the TMD factorization applies to Drell-Yan (DY) type processes for values of the transverse momentum of the vector boson (q T ) much lower then its virtual mass (Q). The factorization theorem concludes that the non-perturbative parts of the q T -differential cross-sections of boson production when q T Q are included in the transverse momentum dependent (TMD) distributions and, separately, their evolution kernel [67][68][69][70][71][72][73][74][75][76][77]. Phenomenologically it has been found that the TMD factorization at leading order works for q T /Q 0.2 [32,34,35] and we consider this range of validity also here. A theoretical estimate of this range has been provided in [78].
The theoretical perturbative calculations for TMD distributions at small-b performed in recent years is highly significative and recently it has reached the N 3 LO precision 1 [79][80][81][82][83][84][85][86][87]. Results of the same N 3 LO already exist for the universal QCD anomalous dimensions [88][89][90][91][92], so that finally one has an extremely accurate perturbative input. At present the non-perturbative parts of the TMD are extracted at NNLO (which means N 3 LL in the evolution kernel and NNLO in the matching of the TMD to PDF) so that for consistency we will use the perturbative results up to this order which are included in the code Artemide [93] and the non-perturbative parameters as extracted in [35]. The theoretical settings of this work are very similar to the ones of [35] and we make explicit use of the ζ -prescription [77,94]. The main difference with respect to the neutral boson case is represented by the particular kinematics of the W -production, which we study in detail in the next sections. Thanks to our explicit study of the leptonic tensor and implementation of fiducial cuts we can compare our predictions with existing results from LHC (ATLAS, CMS) and TeVatron (CDF, D / 0). We explore explicitly new regions of the measured phase space where TMD effects can be significative and experimentally testable, especially for di-lepton masses below the W -mass peak.
In order to establish some notation we start writing the W mediated reaction where h 1,2 are hadrons (typically protons and antiprotons) and l ± = e ± , μ ± and ν l (ν l ) are their corresponding (anti)neutrinos. We consider a cross section differential in the vector boson transverse mass (m T ) and transverse momentum (q T ). The typical TMD condition q T Q becomes in this case q T m T , while the relation among m T and lepton momenta is highly non-trivial as outlined in the next sections. For our predictions, we consider the case of LHC experiments at √ s = 13 TeV, assuming that lepton cuts and fiducial cross sections are similar to the case of √ s = 7 TeV [25]. In principle this discussion can be extended to other center of mass energies. We consider several intervals of m T and q T which are relevant both for the TMD flavor determination and the mass of the W -boson. In this sense we discover that several intervals of these variables can be interesting for QCD studies and we investigate them in detail. We consider also the p Z T / p W T and p W − T / p W + T observables, in different ranges of the transverse invariant mass of the W .
The paper is prepared as in the following. In Sect. 2 we establish the notation and write the cross section of the W distributions with TMD factorization and we explain how fiducial cuts are implemented. The kinematical relations are here described in detail (and we are not aware of a similar detailed description in the literature). In Sect. 3 we list the source of errors for the various observables examined in this paper. In Sect. 4 we provide our prediction for transverse momentum differential distributions. In Sect. 5 we compare the result of our code with existing theoretical and experimental results. We conclude in Sect. 6 and we provide some details of our calculations in the Appendix.

Charged DY cross section with transverse variables
Let us consider a charged DY process as in Eq. (1.1) in which two hadrons with momenta P 1 and P 2 lead to a lepton with momentum l and an (anti)neutrino with momentum l . We approximate the interaction among protons and the leptonneutrino pair through the production of a W ± boson with momentum q = l+l . In this paper we do not consider masses of the initial hadrons nor masses of final-state particles, so The combination of these momenta leads to the definition of the relevant kinematical variables of the DY process where s represents the square of the center of mass energy, Q 2 is the invariant mass of the lepton-neutrino pair and y is the rapidity of the produced W boson. As the neutrino is a nondetected particle the quantity Q 2 is not measurable anymore and one introduces the transverse mass [21,25,60,95-97] as a measurable quantity in processes with invisible particles in final state as where l T and l T are the transverse parts of the lepton and neutrino momenta respectively and φ lν is the relative angle between both particles. The cross section of a general DY process mediated by a gauge-boson G can be written as in [35] where α em = e 2 /4π , being e the electron charge. This cross section is written in terms of a product of a lepton tensor L μν GG and a hadron tensor W GG μν defined as where J G μ is the current for the production of a gauge boson G.
On the other hand, G is the (Feynman) propagator of the gauge boson, in this case G = W and where M W and W are the mass and width of the W boson given in [98], respectively. We do not include electroweak (EW) corrections, which are however calculated in [99].
The W -production is usually expressed in terms of the transverse mass. The relation between the invariant and transverse mass in a general case is non-trivial, so we describe it here in detail. We have and In order to include this function in the cross section we introduce the identity Note that the Dirac delta in Eq. (2.9) is relocated inside the definition of the lepton tensor which simplifies the practical evaluation of the integral.
The next two sections are devoted to the factorization of the hadronic tensor and the final expression for the lepton tensor affected by transverse mass and fiducial cuts. This information will lead us to the final form of the cross section of charged DY within TMD factorization.

Factorization of the hadronic tensor
In this paper we use the unpolarized part of the hadron tensor from [35] neglecting the masses of both hadrons. Thus, for a generic gauge boson G we omit power suppressed highertwist TMD contributions and we write where C V is the matching coefficient for vector current to collinear/anti-collinear vector and z GG f f are the EW factors that will be defined later. On the other hand f 1 is the unpolarized TMDPDF defined as where p is the momentum of the hadron, W n (x) is a Wilson line rooted at x and pointing along vector n to infinity. The momentum fraction x is defined in DY kinematics as (2.13) In the small-b limit, the TMDPDF can be re-factorized in terms of matching coefficients (calculated up to NNLO in [82]) and integrated PDFs. In order to write a complete TMD, we consider also a function modeling non-perturbative effects that are not included into collinear PDFs (2.14) where the symbol ⊗ represents the Mellin convolution in the x variable. We use the ansatz for f N P suggested in [35] where the parameters λ 1 , . . . , λ 5 are extracted from a combined DY+SIDIS fit to data in [35]. In the TMD parameterization of Eq. (2.14) the flavor dependence is in principle contained in the parameters λ i and the PDF. The extraction of [35] has found that, given actual SIDIS data it is sufficient to count flavor dependence in PDF only. While we make predictions using the values of λ i as in [35], it is possible that the W -boson measurements that we propose will be able to achieve a better control of these parameters, including a flavor dependence. The angular part of the Fourier integral in Eq. (2.11) can be solved analytically and we can rewrite the hadron tensor where

Lepton tensor and fiducial cuts
The lepton tensor that enters in the cross section is integrated over the lepton and neutrino momenta, thus contracting it with the Lorentz structure that comes from the hadron tensor in Eq. (2.11) we obtain, where the extra label cuts represents the fiducial cuts introduced in different experiments. These cuts are implemented over the momenta of lepton and neutrino and the rapidity of the lepton On the other hand note that the difference of this lepton tensor with the one for neutral DY production [32,35] is the extra dependence on the transverse mass through the delta function introduced in Eq. (2.9). The integral in Eq. (2.18) cannot be solved analytically and only numerical results can be obtained. A detailed discussion about this integral can be found in Appendix A.

Final expression of the cross section
Once we have precise definitions for hadron and lepton tensor and taking into account that we get the desired differential cross section as where, in the case of W boson production, the product of EW factors for hadronic and leptonic part and propagators can be written as where V f f are the elements of CKM matrix that mixes flavors, e f are the quark charges in terms of the electron charge and (2.24) Thus, the final expression for the desired cross section is (2.25) The plots of next sections are done using the PDF set of NNPDF31_nnlo_as_0118 [28]. The numerical inputs that we have used are

Sources of errors
In the course of the paper we consider several sources of errors that come from the way we write and parameterize the cross section. We list them here as a reference for the following sections.

Scale variations
In the present work we use the ζprescription defined in [77] and we include the nonperturbative fixing of the ζ -scale as defined in [35,94]. As a result the scale variation is done changing the parameters c 2,4 ∈ [0.5, 2]. These parameters define the uncertainty in the hard matching scale and in the matching of the TMD on the collinear PDF respectively. Error from a reference PDF set using replicas The central value of our predictions is deduced using the PDF set NNPDF31_nnlo_as_0118 [28], which includes also LHC data. The error on each bin is evaluated taking the variance over 1000 replicas.
Predictions from different PDF sets In [35] the authors have performed fits of DY data at low and high energy with different sets of PDF and we report their results in Appendix B. We consider here the same sets of PDF and we show how the predictions change with these different sets.
Error from TMD parametrization In [35] the authors provide a set of replicas of their TMD parameterization and we provide the main estimation of this error using them. Alternatively, one could use a more standard error propagation coming from the estimate of the non-perturbative parameters. We have used this second method as a check, and we provide more details in Appendix C.

Observables in W production with TMD
In this section we concentrate on some observables that can be, in principle, sensitive to TMD effects and/or are relevant to establish some important properties of the W as its mass. We consider the W -boson differential cross section, the ratios for p Z T / p W T , and p W − T / p W + T all of them as a function of the boson transverse momentum distribution and different intervals of the transverse mass. The fiducial cuts have been set as the ATLAS experiment [25] but for √ s = 13 TeV. We have identified three interesting intervals of m T , namely [50,66] GeV, [66,99] GeV, [99,120] GeV corresponding to m T values below, around, above the W -mass, with typical cuts on lepton momenta. The m T -interval which is mostly studied in the literature is the one with m T ∈ [66,99] GeV, however this is not the only interesting one from the perspective of studying the TMD properties. This is because non-perturbative TMD effect are expected to be more relevant for low values of m T , similarly to Drell-Yan processes.
The q T interval that we have considered, that is q T /m T 0.2 is consistent with the application for TMD factorization. In order to have a reference setup for PDF we have chosen NNPDF31_nnlo_as_0118 [28] as in [35] as it is one of the sets that include LHC data. Nevertheless we have also examined the results taking into account different sets. We have considered a q T binning of 1 GeV, which much probably exceeds the possibilities of current experiments. We motivate this choice because it has been used also by other authors, see f.i. [60], and also because it is allows to better evidence the TMD effects.

Spectra of the W ±
We have considered the spectra of the W ± and in Fig. 1 we show the case for the W − (similar plots are obtained for W + ). In the upper panels of Fig. 1 we have considered three possible cuts for lepton momenta as reported in the same figures. In the upper panels of Fig In the lower panels of Fig. 1 we normalize the cross section to its value integrated in q T in the interval shown in the figure. The shape of the curves now changes and in the left and central panels the more peaked cross sections is the one obtained for m T ∈ [50,66] GeV, which is partly due to the fact that we always select q T /m T ≤ 0.2.
In order to drive a conclusion from these plots we have to recall the previous experience of the fit of [34,35] and also the results of [36]. In these works it shown that in TMD analyses increasing the value of Q the non-perturbative QCD effects are washed out. As a a result having data below the W, Z boson peaks is extremely useful for this kind of research and it can provide valuable information. The plots shown in Fig. 1 actually show that within the current facility it is possible to achieve this goal. The bands shown in Fig. 1   (3) the uncertainty as coming from different sets of PDF; (4) the uncertainty due to nonperturbative parameters: in this case the central line is given by the the central replica of NNPDF31_nnlo_as_0118 [28].
There is a series of observations that one can make on these plots. As a first the smallest variations are provided by replicas and non-perturbative parameters and both of them are below 1%. The scale variation is certainly the biggest source of error: it is in range 2-9% for q T ≤ 4 GeV, and about 1-2% for q T ≥ 4 GeV. We recall that the present analysis uses N 3 LL for TMD evolution and NNLO for all the rest. In the introduction it was pointed out that recently higher order perturbative calculations have been performed. It will be interesting to observe how the scale variation changes when including this higher order term. On the other side before doing this one should have a TMD extraction of the same order.
In exploring different PDF extractions we have considered the ones for which a TMD extraction has been provided in [35] . We recall that the HERA PDF [100] actually provided the best fit in that case. Observing Figs. 2 and 3 one sees that actually this set behaves differently with respect to the others, however the difference with NNPDF31 is below 1-2 % an all over the range of qT that we have considered. For all the other sets, there is a major difference with NNPDF31 for q T ≤4 GeV while they agree within 1% for greater values of q T .
The uncertainties due to non-perturbative parameters result to be almost negligible. The non-perturbative param-eters have uncertainties coming from their extraction from data, that is from a χ 2 analysis of fits. Despite the fact that some of these uncertainties may be significative on each single parameter the impact of each single variation is just a fraction of the impact of the whole non-perturbative TMD structure because their value is significantly different from zero. It is also possible that the errors on each parameter is under-estimated as suggested in [35].

Ratio p Z T / p W T
It has been pointed out since a long time that the observable can show a very small uncertainty because the PDF contributions tend to cancel and it can be useful for the measurements of the W mass (some recent phenomenological work can be found in [101,102]). The region of q T where this observable is actually measured corresponds to the region of validity of the TMD factorization theorem, so it is interesting to observe what is the impact of TMD. In formula 4.1 the numerator and the denominator are weighted by σ W,Z that is the cross section integrated in the q T interval under study. This is preferable to normalizing to the total cross section σ W,Z which should be extracted elsewhere. The result of our prediction is shown in Figs. 4 and 5 for W − and W + respectively. We have considered two intervals of m T , being the lower one useful to control better TMD effects and several sources of error as in the previous section. The scale uncertainty now is  [66,99] GeV (middle column), [99,120] GeV (right column). In the first line we report the theoretical error from scale uncertainties as explained in the text. In the second line we have the error calculated as a variance in each bin of 100 replicas of the set NNPDF31_nnlo_as_0118 [28]. The uncertainty is referred to the average value of each bin (red line). The value of the cross section given by the central replica is represented by the green line. On the third line we represent the value of each been with different sets of PDF. On the fourth row we have the uncertainty due to nonperturbative parameters. The central value is given by the the central replica of NNPDF31_nnlo_as_0118 [28] treated as in [60], considering separately the correlated and uncorrelated cases.
The error from scale variations can be considered as correlated or uncorrelated as in [60], see second and third row in Fig. 4. The uncorrelated uncertainty is well below 1% for the central interval and most of the low interval of m T . On the other side the correlated uncertainty is below 2% only for q T ≥ 4 GeV. A similar trend is shown when different sets of PDF are used. In this case the spread of the results can be at most 2% for q T ≤ 4 GeV, and at most 1% for q T ≥ 4 GeV. We find remarkable that the PDF error results slightly inferior for the low m T interval, which confirms that this is an interesting case to study. Another interesting observation is that the difference between the HERA and NNPDF31 sets is always well below 1%. The other sources of error that we have considered, give uncertainties less the 1% on all over the intervals that we have studied.

Ratio p W
The ratio of W ± cross sections shown in Fig. 6. We have considered again two intervals of m T and it is nice to observe the similarity of these curves in the two intervals. The theoretical errors have been estimated and scale variations represent In the first line we report the theoretical error from scale uncertainties as explained in the text. In the second line we have the error calculated as a variance in each bin of 100 replicas of the set NNPDF31_nnlo_as_0118 [28]. The uncertainty is referred to the average value of each bin (red line). The value of the cross section given by the central replica is represented by the green line. On the third line we represent the value of each been with different sets of PDF. On the fourth row we have the uncertainty due to non-perturbative parameters. The central value is given by the the central replica of NNPDF31_nnlo_as_0118 [28]. The cross section for the central replica is shown in Fig. 7 the biggest error. We have considered both correlated and uncorrelated scale variations and their value is very similar in this case. The scale variations keeps below 2% for q T > 4 GeV and grows below this value. All other uncertainties keep below the 1%.

Comparisons with other groups and experiments
The study of TMD in W -production has not been explored in its full potential in the literature, because the formulation of the TMD factorization theorem with full details is very recent. The main difficulty in establishing this kind of studies is that the TMD factorization only holds for q T Q, and in this W -case q T m T . In this way we are interested in describing only a (relevant) part of the spectrum, the one which contains the peak of the distribution.
In this section we would like to compare the outcome of the results as coming from TMD factorization and recent TMD fits with theoretical and experimental results.
A theoretical prediction of the LHC case can be found in [60] and the comparison with us is done in Fig. 7. In Fig. 4 Ratio of Z /W − spectrum for m T ∈ [50,66] GeV (left column) , and m T ∈ [66,99] GeV (right column). On the first row we have the Z /W − spectrum using the central replica of NNPDF31_nnlo_as_0118 [28]. Uncorrelated and correlated theoretical uncertainties are given in second and third row respectively. In the third line we have the error calculated as a variance in each bin of 100 replicas of the set NNPDF31_nnlo_as_0118 [28]. The uncertainty is referred to the average value of each bin (red line). The value of the observable given by the central replica is represented by the green line. On the fourth line we represent the value of each been with different sets of PDF. On the fifth row we have the uncertainty due to non-perturbative parameters . On the first row we have the Z /W − spectrum using the central replica of NNPDF31_nnlo_as_0118 [28]. Uncorrelated and correlated theoretical uncertainties are given in second and third row respectively. In the third line we have the error calculated as a variance in each bin of 100 replicas of the set NNPDF31_nnlo_as_0118 [28]. The uncertainty is referred to the average value of each bin (red line). The value of the observable given by the central replica is represented by the green line. On the fourth line we represent the value of each been with different sets of PDF. On the fifth row we have the uncertainty due to non-perturbative parameters Fig. 6 Ratio of W − /W + spectrum for m T ∈ [50,66] GeV (left column) , and m T ∈ [66,99] GeV (right column). On the first row we have the W − /W + spectrum using the central replica of NNPDF31_nnlo_as_0118 [28]. Uncorrelated and correlated theoretical uncertainties are given in second and third row respectively. In the third line we have the error calculated as a variance in each bin of 100 replicas of the set NNPDF31_nnlo_as_0118 [28]. The uncertainty is referred to the average value of each bin (red line). The value of the observable given by the central replica is undistinguishable from the red line. On the fourth line we represent the value of each been with different sets of PDF. On the fifth row we have the uncertainty due to non-perturbative parameters Fig. 7 Comparison of our prediction with [60] including errors coming from scales variation. On the right the two cross sections are normalized to the central value of the prediction from Artemide shown. The left, central, right panels correspond respectively to the cases of 50 GeV < m T < 66; GeV, 66 GeV < m T < 99; GeV, 99 GeV < m T < 120; GeV this plot we show the cross section for the two groups in the interval 0 GeV < q T < 20 GeV. The cross section is normalized to the total cross section as provided by [60], in order to have a reference value. The scale uncertainty in the two curves is similar, although slightly reduced in our case. There can be multiple reasons for this, like a different choice of scales or the usage of Monte Carlo calculation in [60]. One observes also a remarkable difference in the central values of the cross section for very low values of q T ≤ 5 GeV. The origin of the difference may have multiple motives starting from the parametrization of the non-perturbative effects (which is based on a broader data analysis in our case) to the set of prescription and scale fixing used by the two groups. In general this difference is expected where non-perturbative effects are significative as it is shown by the figure, while it is much less significative for q T ≥ 5 GeV. In order to explore this we have considered a modified version of the TMD where the effect of the constants λ i in Eq. (2.15) is nullified. This is achieved considering f N P (x, b) = exp(−λb 2 ) with λ 10 −3 . The results are shown in Fig. 8 and they show to be consistent with Fig. 7 where the TMD non-perturbative part was not considered, despite the fact that the kinematical cuts considered in the two figure are slightly different. We postpone anyhow a deeper study to a different work.
We have considered also a comparison with Pythia 8.3 [103], with AZ tune as adapted by ATLAS experiment. We have found two version of this tune one as used in [60] 2 In the comparison of the cross sections Fig. 9 we observe a general shift of the Pythia result with respect to ours, whose sign depends on the value of q T , although the two estimates are mostly compatible within the errors (more for 66 GeV < m T < 99 GeV than in the less energetic interval). In the ratio of W ± cross section (Fig. 10) instead the agreement is complete, both considering correlate and uncorrelated errors, providing Artemide a smaller error band. Similar conclusion come from the p Z T / p W T in Fig. 11, although with more difference for transverse momentum less than 5 GeV.
Current experiments at LHC at the moment have a limited number of data points. For ATLAS we have found that for the low q T region, one has only two points in [95], shown in Fig. 12. The cross section is normalized to its integrated value in the plotted q T interval. The theoretical and experimental values are compatible within the errors. However it is evident that such a large binning does not allow a precise statement. 2 We thank Pier Monni for communicating us that the AZ tune as provided by the Pythia collaboration does not coincide exactly with the one used in ATLAS experiment and used actually in [60] and one as provided by ATLAS experiment. We have included both in our plots which are shown in Figs. 9, 10, 11. The pictures show the differential cross section normalized to their integration over the shown interval. The errors from Pythia come from statistical uncertainty and we have checked that they are similar to the one obtained by the variation of the parameters of the tune. Fig. 9 Comparison of Artemide cross section with Pythia 8.3 AZ tune as in [60] (blue band) and as in the original ATLAS release (green band) for W + (top panels) and W − (low panels). The Artemide error comes from scale variations, the Pythia errors are commented in the text  Fig. 13 we consider the cases of the W , Z /W and W − /W + transverse momentum spectrum comparing our results with [21]. The large bin for q T > 18 GeV is the most sensitive to power corrections which however do not seem to create particular problems. The first thing that one can observe is that there is a good agreement between data and our prediction when theoretical errors are included, which suggests the fact that most of the QCD corrections in the CMS experiment are due to the TMD region.
Going to TeVatron experiments we have found data for D / 0 and CDF and we show a comparison with our prediction in Fig. 14. In the central plot of this figure we do not show the experimental error, because the uncertainty of the bins are correlated through non-trivial correlation matrices. The agreement with these experiments is in general greater than We have considered the cross sections χ 2 for the different experiments with different sets of PDF and the results are shown in Table 1 (we also report the number of relevant points for each experiment). The CDF result are the ones with larger errors and TMD predictions agrees with them while for the D / 0 at √ s = 1.8 GeV, the agreement is worse. We have not considered the case of D / 0 at √ s = 1.96 GeV, because the computation of the error in this case involves the knowledge of the W-spectrum up to q T = 600 GeV, which we do not have. In the case of LHC, in general we have a very limited number of points for each experiment. ATLAS has a remarkable agreement with the HERA20 PDF sets, definitely better than with other sets. The CMS case is more elaborate. For the case of electronic decay of W at CMS, the high value of χ 2 is basically driven by just one point out of four, as can be seen Fig. 13, left panel (it is the point in the bin 7.5 GeV< q T < 12.5 GeV). Removing this point, the χ 2 is very similar to the ATLAS case. For the muon channel instead we cannot find a particular point which is responsible for the high χ 2 . As a final remark we recall that the CMS observable is totally inclusive on m T , which does not allow a perfect control of factorization hypothesis at the theoretical level. We have checked the impact on the χ 2 of the theoretical errors (scale uncertainty and PDF replicas) for the case of NNPDF and HERA20 and we have reported it in Table 2. In both cases the errors on each bin are considered uncorrelated to the rest. The scale uncertainty error is dominated by the variation of c 4 , the parameter which associated to the scale at which TMD are matched onto PDF. We observe that including this error we have a big reduction of the χ 2 on all experiments. In order to understand better this issue it would help to have experimental results with a definite interval of m T and also a TMD extraction that includes the W processes, that we postpone to a future work. Concerning this last point we have compared the χ 2 as coming The different extractions are obtained with the NNPDF31, PDF set, and using different data sets as shown in Table 3. We find a substantial agreement among all extractions, which suggests a mild flavor dependence of the TMD. In Table 2, the muon channel in CMS is still not fully agreeing with the theoretical prediction. This fact, which is put in evidence here for the first time, needs further study beyond the present work.

Summary and outlook
In this work we have set the status of the current knowledge of the W -boson spectrum within the TMD factorization formulation of its cross section and using the latest extraction of these distributions in [35]. We have considered the W -production differential cross section, the p Z T / p W T and the p W + T / p W − T distributions as functions of the boson transverse momentum q T in different intervals of m T . The kinematic of the processes has been deeply studied, providing details that were not treated in the literature (to our knowledge),   and showing explicitly how fiducial cuts are implemented together with kinematical power corrections. The kinematical description so achieved has put us in the condition to use the TMD extraction made in [35] with the same level of precision. The perturbative inputs are N 3 LL for the evolution and NNLO for all the rest. Because of our experience with Z -boson production we have given particular attention to values of m T below ( [50,66] GeV) and around the Wmass ( [66,99] GeV). In fact one expects to better control the non-perturbative QCD effects once the details of the observables below the W -mass are also known. It is remarkable that observables at low values for m T have been poorly studied despite the fact that there is no particular suppression of the cross section and all errors look very similar. An important chapter is the error analysis. We have considered as sources of error, scale variations, replicas uncertainty within the NNPDF31 PDF set, uncertainties due to different sets of PDF and the ones due to TMD non-perturbative parameterization. This analysis is interesting because it gives a first step in our understanding the errors. We consider it an important starting point that can certainly be improved with future work within TMD studies. For instance scale uncertainties can be checked including one more perturbative order (including very recent results at N 3 LO, [83][84][85][86][87]), which how-ever require a new TMD extraction. Then, the sets of PDF that we have checked are the ones analyzed also in [35], which is still limited. More PDF sets should be used and included in the TMD extraction. The actual impact of each source of error depends on the observable, and we have described it in Sect. 4.
At present the available data provide already some information on the observables that we have studied. We have considered the cases of TeVatron and LHC. The agreement in the cross section is reasonable for CDF, D / 0 at √ s = 1.8 TeV, ATLAS, while it is not completely satisfactory CMS. For all these cross sections we have done a χ 2 analysis. The message is that inferring and estimating QCD non-perturbative contributions to W processes from neutral boson mediated ones has some effectiveness and that a more global treatment of data is also worth studying. This raises also the question of whether a more sophisticated TMD flavor dependence can improve the agreement with data. Other effects like QED contributions [104,105] will also be addressed in future studies.
In order to perform this study we have proposed to consider several intervals of m T , for which we do not observe a particular suppression of the cross section. Comparing experiments, we find a striking difference between LHC and TeVatron data in the q T binning of the final result. Definitely the 2 GeV binning of TeVatron allows a much better understanding of the QCD effects, despite the LHC precision. A similar binning would be also desirable for LHC. This problem has already been discussed also in the case of Drell-Yan in [36]. valuable  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 .
For simplicity we can rewrite the term coming from the Jacobian in Eq. (A.5) as Finally the term that only depends on momenta in Eq. (A.2) is rewritten as Thus, the final result is written in term of the auxiliary functions J, A, M as a one-dimensional integral T , q T , cos φ)θ(cuts), (A.13) that should be done numerically.

B Inputs from fits
In this section we report some results from previous fit [35] that we have used in this paper. The different PDF sets that have been reported in that fit and that we have used are listed in Table 4. The results for the TMD constants are in Table 5.

C Errors propagation from TMD parameters
An alternative way to present the errors coming from nonperturbative parametrization of the TMD is the following. We take from [35] the different values for the non-perturbative parameters coming from fits using several PDF sets. We have a data array of seven non-pertubative parameters and five values for each one with its respective uncertainty as it is shown in Table 5. We calculate the mean value of each parameter and its uncertainty associated doing the square sum of the uncertainties and taking the square-root of the sum. We show the values in Table 6. In order to quantify the uncertainty in our cross section due to the uncertainties in our non-pertubative parameters we use Eq. (C.1) where θ i = {B NP , c 0 , λ 1 , λ 2 , λ 3 , λ 4 , λ 5 } are the nonperturbative parameters andθ i their mean values collected in Table 5; V i j and C i j are the covariance and correlation matrices respectively of our non-perturbative parameters and δ dσ dq T i is the difference of the differential cross section evaluated in the extremes of the interval of the non-perturbative parameter i and keeping the rest of nonperturbative parameters fixed at their mean value. In order to go from the first line to the second line, we have approximated the derivative of the cross section with respect to θ i by where δθ i is the uncertainty associated to the mean value of the non-perturbative parameter θ i . Using the relation between the correlation and covariance matrices C i j = V i j δθ i δθ j we arrive to Eq. (C.1).   The correlation matrix at NNLO for the λ i is provided in [35] and, as a reference, we have considered the central PDF of NNPDF31_nnlo_as_0118 [28].
The result of this error estimate is very similar to the one described in the main part of the paper, which is a confirmation of the Gaussian distribution of TMD parameter errors for this case.