Neutral Current Forward-Backward Asymmetry: from $\theta_W$ to PDF Determinations

Measurements of the forward-backward asymmetry in neutral-current Drell-Yan di-lepton production have primarily been used for determinations of the weak mixing angle $\theta_W$. We observe that, unlike the case of Run-I of the Large Hadron Collider (LHC Run-I), for the first time at the LHC Run-II the reconstructed forward-backward asymmetry has the capability of placing useful constraints on the determination of the parton distribution functions (PDFs). By examining the statistical and the PDF uncertainties on the reconstructed forward-backward asymmetry, we investigate its potential for disentangling the flavour content of quark and antiquark PDFs. Access to the valence/sea $u$-quark and to the sea up-type antiquark PDFs, in particular, may be gained by the appropriate use of selection cuts in the rapidity of the emerging lepton pair in regions both near the $Z$-boson peak and away from it, in a manner complementary, though more indirect, to the case of the charged-current asymmetry. We study the extension of these results for the planned high-luminosity (HL) LHC.


I. INTRODUCTION
The forward-backward asymmetry in neutral-current (NC) Drell-Yan (DY) lepton pair production has traditionally been used at hadron-hadron colliders in the context of precision electroweak measurements, in particular for determinations of the weak mixing angle θ W [1][2][3]. However, mostly due to the theoretical systematics coming from the parton distribution functions (PDFs), the accuracy reached in θ W determinations from measurements of the forwardbackward asymmetry at the Large Hadron Collider (LHC) is currently an order of magnitude lower than in the case of θ W determinations from lepton colliders [4].
In Ref. [5] we have proposed that measurements of the NC forward-backward asymmetry (henceforth denoted as A FB ) at the LHC Run-II (as well as in future high-luminosity runs at the HL-LHC) can usefully be employed in PDF determinations, in a fashion complementary to measurements of the charged-current (CC) asymmetry, which have long been used in PDF fits to access the d/u quark PDF ratio (see [6][7][8][9][10][11] and references therein). While CC and NC doubly differential cross sections in the DY di-lepton mass and rapidity have extensively been used for PDF determinations, Ref. [5] investigates the role of the angular information encoded in the A FB , which is related to the single-lepton pseudorapidity and, once combined with di-lepton mass and rapidity, would qualitatively correspond to triply differential cross sections. Precision measurements of such cross sections [12] have not yet been exploited systematically for PDF determinations. It is noted in [5] that, once mapped in the invariant mass of the di-lepton final state, both the NC differential cross section and the A FB display statistical errors which, with the luminosity reachable at Run-II, become competitive with those assigned to quark PDFs in a broad kinematic range around the Z-boson peak.
In this paper we follow the proposal [5] and explore the NC asymmetry focusing on its potential for disentangling the flavour content of PDFs. Specifically, we find that it is possible to probe the quark PDF flavour structure with the A FB through a judicious use of selection cuts in the rapidity of the emerging di-lepton pair, whose four-momentum is completely accessible in NC DY, unlike the case of the final state of CC DY which contains a neutrino. Although access to the d/u quark PDF ratio is more indirect with the NC A FB than with the CC asymmetry, the PDF kinematical region in longitudinal momentum fraction x and evolution scale Q 2 is larger in the NC case, as both on-peak and off-peak cross sections can be used, whilst in the CC case the escaping neutrino prevents the reconstruction of the W ± mass so that the W ± peak region necessarily dominates the cross section. We observe that selecting high dilepton rapidities enhances the contribution of higher-order radiative QCD effects [13,14] and subleading production channels [15,16] so that this region also becomes sensitive to the gluon PDF.
The paper is organised as follows. In Sec. II we briefly summarise the set-up of the calculation of Ref. [5] (which was applied to the CT14NNLO [6] and NNPDF3.1 [7] PDF sets) and we extend its results to the case of the ABMP16 [8], HERA2.0 [10] and MMHT2014 [9] PDF sets. In Sec. III we address the question of flavour content and investigate parton luminosities in the kinematic regions relevant for the A FB around the Z-boson peak. In Sec. IV we analyse the A FB in the high-rapidity region and we study the role of rapidity selection cuts to disentangle the PDF flavour structure. We give conclusions in Sec. V.

II. THE DIFFERENTIAL FORWARD-BACKWARD ASYMMETRY
The forward-backward asymmetry A FB in the γ * , Z mediated DY production process pp → γ * , Z → + − (with = e, µ) is given by where M ( + − ) is the DY lepton pair invariant mass, η( + ) is charged lepton pseudorapidity while the identification of the Forward (F) and Backward (B) hemispheres via the restrictions η( + ) > 0 and η( + ) < 0 is obtained by exploiting the di-lepton boost variable as discussed in [17][18][19][20][21]. Hence, this is effectively a 'di-lepton FB asymmetry'. This asymmetry is to be contrasted with the 'lepton charge asymmetry' in the CC channel, defined as which has traditionally been used to constrain PDFs, particularly the ratio of d-quark to u-quark distributions. From the charged current data, different observables can be actually reconstructed. The two cross sections σ(ud → + ν) and σ(dū → −n u) give rise to different distributions owing to the different valence u and d-quarks that dominate their production. The distribution in the rapidity of the mediator, W ± , is governed by a sea quark and a sea anti-quark at low values. Only a tiny difference is thus expected between the two above-mentioned cross sections. Instead, at large values of y(W ± ), the process is mainly initiated by a valence quark and a sea anti-quark. A much bigger discrepancy is therefore expected between the two charged current cross sections. This feature should enable one to distinguish between valence/sea u/d quark PDFs. This highlights the potential role of a binning (even a cut) on the rapidity of the di-lepton system for charged current processes and in general. Starting from the well known relations TeV and x i = | p i |/E beam , the W ± rapidity distributions can be directly related to the quark and anti-quark fractional momenta thus constraining the corresponding PDFs. For |y(W ± )| 2.5, the measurable region would correspond to x 1 0.1 and x 2 3 · 10 −4 (low-x anti-quark and large-x quark as the charged lepton is likely to be emitted in the same direction as the valence quark). This range could be used to constrain the PDFs and determine the parton luminosities.
For the CC leptonic channel, the W ± -rapidity can be derived using the reconstruction of the neutrino on the Wboson peak, √ŝ M W . This reconstruction algorithm however suffers from uncertainties. A better observable is the charged lepton pseudo-rapidity that is directly related to the scattering angle, η(l ± ) = − ln(tan θ/2). The difference between the two processes, qq → W ± → l ± ν l , is then enhanced at high values of |η(l ± )| of course compatible with the detector acceptance.
In order to reduce the systematics, including the theoretical uncertainty due to the higher order QCD and EW corrections, it is a standard practise to construct cross section ratios. That is the reason to define the lepton charge asymmetry in Eq. II.2. This strategy is also highly advisable to reduce the theoretical error due to the PDFs determination. The PDFs uncertainty indeed cancels out, at a large extent, when considering cross section ratios. This observation is the kernel of our analysis, which addresses the neutral current DY processes qq → Z, γ * → l + l − . This process is used to analyse the Q 2 dependence of the PDFs. While the CC process is sensitive to the sea anti-quark and the valence quark, the NC process encodes informations about the sum of the sead,ū anti-quarks and the valence u, d quarks. Oppositely to the CC differential cross sections that appear in ratios, the NC ones are not normalised usually. Recently, in the literature, the LHC NC data have been released in the form of triple differential cross sections in the di-lepton rapidity, di-lepton invariant mass and lepton scattering angle. This powerful informations however suffer from experimental and theoretical systematical errors. In particular, the theoretical PDF's uncertainty can be so large to bury the potential differences coming from the comparison of the NC cross sections obtained using different PDF's sets. This would of course weaken enormously the sensitivity of the LHC measurements to the q andq PDFs. At present, the discrepancy between the predictions from any two different PDF sets is compared to the statistical precision. At the W ± or Z-boson peaks, the cross sections are very high. The signal is clear and the background is rather low. That is why a sensitivity to small differences in the PDF's parametrisation is claimed to be at hand.
However, in our view, the PDF's uncertainty could still represent a real obstacle. This is the main reason for proposing the continuation of the use of the Forward-Backward asymmetry, which is a ratio of (differential) cross sections in the NC channel. It was observed in [5], using the method of [19,20] to define the reconstructed FB asymmetry, A * FB , and to evaluate its PDF uncertainty, that whilst no gain in quark PDF determinations should be expected from the reconstruction of this observable with the data collected at √ s = 8 TeV, the scenario changes at the LHC Run-II and at the planned HL-LHC. For instance, for the projected luminosity L = 3 ab −1 of the HL LHC [22], and for the CT14NNLO [6] and NNPDF3.1 [7] PDF sets, the statistical error is found [5] to be lower than the PDF error in the low to intermediate mass region M ∼ < 400 GeV.
Furthermore, as previously mentioned, the reconstructed A * FB experiences a partial cancellation of systematical uncertainties. As the uncorrelated and correlated components of the uncertainty are determined for A * FB and added up in quadrature, the result is that the global error is significantly reduced compared to the case of the pure differential cross sections. From this point of view, the measured A * FB could be highly competitive with the single-differential cross section and even with the recently measured triple-differential cross section [12] from which it can be obtained.
We here consider the differential distribution of the A * FB observable in the di-electron final state, focusing on the invariant mass region around the Z-boson pole, where the high statistics due to its resonant production allows very precise experimental measurements. The results are presented taking into account the resolution and the experimental acceptance and efficiency factors of this channel [23]. When computing the statistical error on the observable, we include the QCD next-to-next-to-leading order (NNLO) corrections to the cross section [24,25] to estimate the expected number of events.
Following up on the results we already presented in [5], here we further explore the effects of imposing a rapidity cut on the di-lepton system on the observable A * FB . As it is well known from the literature, such a kinematical cut improves the efficiency of the reconstruction so that a closer shape to the original A FB is recovered. Conversely, as we mentioned before for the CC channel, the discrepancy in the A * FB predictions obtained with any two different PDF sets are enhanced by the effect of the rapidity cut. We extend here the results of [5] by including additional PDF sets.
In Fig. 1, the A * FB distributions in the di-lepton invariant mass have been obtained with three different PDF sets: ABMP16 [8], HERA2.0 [10] and MMHT2014 [9]. We have imposed a lower rapidity cut of |y| > 0 (first row), |y| > 0.8 (second row) and |y| > 1.5 (third row) on the di-electron final state. The three distributions are accompanied by the PDF (plots on the right) and the statistical (plots on the left) error bands. The latter have been computed assuming an integrated luminosity of 300 fb −1 , that is, the value which is supposed to be achieved by the end of the LHC Run-II programme.
The increasing discrepancies between the predictions as the rapidity cut gets stronger are a consequence of the different parametrisations in the partonic content of the proton, characteristic of each PDF's set. Quarks and antiquarks PDFs suffer from larger uncertainties in the regions of high-x and low-x, respectively: the rapidity cut indeed selects these kinematical regions. Future experimental measurements of A * FB in these invariant mass and high rapidity regions have the potential of reaching sufficient precision to discriminate between different PDF parametrisations in the corresponding x range. If using A * FB , which is a cross section ratio, this discrimination can be remarkably outside the PDF's uncertainty bands, as shown by the right-hand side plots. This is true in the invariant mass interval just before the Z-boson peak, 60 GeV≤ M (ll) ≤ 80 GeV, and soon after it, 100 GeV≤ M (ll) ≤ 180 GeV, for a high enough y(ll) cut. This condition is not at all guaranteed by the triple differential cross section that contains the same informations but is not normalised. This powerful feature of the A * FB observable is not taken into account, generally, and it should.
From Fig. 1, one can clearly distinguish the HERA2.0 predictions from the ABMP16 and MMHT2014 ones already for y(ll) ≥ 0.8. Analogous results can be found in the comparison between the NNPDF3.1 and the CT14NNLO predictions as displayed in Ref. [5].

III. PARTON LUMINOSITIES AT HIGH RAPIDITY
In this section we want to further investigate the origin of the behaviour of A * FB that we have observed. We'll show that such a behaviour is strictly related to the flavour content of the proton and can give access to the up-type PDFs in a unique way. To this aim, we compare the parton luminosities from the u −ū and the d −d interactions at some fixed value of the centre-of-mass energy around the Z-boson peak focusing on their dependence on the partonic rapidity, y qq . We follow the formulae of Ref. [26] and in Fig. 2 we show their predictions for the CT14NNLO set.
As visible, the d −d initiated process gets more suppressed at high rapidity in comparison to the u −ū one. This effect gets enhanced with increasing the centre-of-mass energy of the hard quark-antiquark scattering. It's indeed particularly evident on the righthand-side of the Z-boson peak. It is then interesting to visualise the ratio of the u −ū and d −d partonic luminosities as a function of the c.o.m. rapidity for the various PDF sets that we are considering. This is shown in Figs. 3(a)-(e) for the four values of the centre-of-mass energy of the hard quark-antiquark scattering that, on purpose, have been selected before and after the Z-boson peak (see legend). While at low rapidity, y qq , the two luminosities are comparable, with increasing y qq the d −d component gets depleted by roughly an order of magnitude. The level of suppression of the d −d contribution is different within each PDF set. In order to estimate such differences in the predictions, one needs to include the PDF errors from each set.
In Fig. 3(f), the ratio of the u −ū and d −d partonic luminosities as a function of the c.o.m. rapidity is displayed with the corresponding PDF error band as calculated within the various PDF sets at fixed √ŝ = M Z . In order to detail the displayed differences in the ratio predicted by the various collaborations, we take as reference point the rapidity value y = 4.5. Here, the relative contribution from the down-quark initiated process is 12% with a relative error of 82% with the NNPDF3.1 set, 21% with a relative error of 36% with the CT14NNLO set, 12% with a relative error of 17% with the ABMP16, 18% with a relative error of 23% with the HERA2.0 set and 20% with a relative error of 20% with the MMHT2014 set.
This has two consequences. First of all, there are differences expected in the Z-boson production cross sections obtained with the various PDF sets that are well outside the PDF's error band at large rapidity. A good sensitivity could be reached to small differences in the u, d quark and anti-quark parametrisation. Second, by imposing a sufficiently high rapidity cut on the data sample, all the considered PDF sets agree within errors that the dd component of the luminosity gets reduced to at most 20%. One could therefore have a direct handle on the processes that are initiated by the u −ū interactions. This implies that one could have direct access to the valence u-quark and the sea u quark and anti-quark PDFs. Up to some degree of purity, the up-quark flavour could be extracted from the neutral Drell-Yan data and the corresponding PDFs constrained (or tuned) by a direct comparison. Using the valence u-quark PDF already well measured, one could adjust the sea u andū distributions to fit the data. The kinematical region of large rapidity corresponds to select a low xū and a high x u in the uū scattering process. These are the ranges least known for the quarks PDFs. Therefore, at large rapidity the neutral Drell-Yan data could be used to constrain the u andū PDFs and determine the corresponding uū parton luminosity in these two extreme regions of x.

IV. SENSITIVITY TO THE u ANDū PDFS
We now want to exploit this idea as far as the LHC sensitivity would allow. For a meaningful analysis we cannot push the rapidity cut too high due to the reduction of the data sample though. In Fig. 4, we show the expected number of events as a function of the di-lepton invariant mass obtained for different PDF choices when imposing a rapidity cut |y| > 4.5. For this value, the measurable region would correspond to x 1 0.6 and x 2 7 · 10 −4 (low-x anti-quark and high-x quark). This range could be used to constrain the valence/sea up-quarks and sea up-type antiquarks PDFs, where they are least known, and determine the corresponding parton luminosities. In the plot on the left we can observe the statistical error band evaluated for an integrated luminosity of 300 fb −1 . As visible, in the region around the Z-boson peak the statistics is very high despite the stiff kinematical cut, however, the distribution rapidly falls moving away from the resonant production. Up to M ll ≤ 120 GeV, the various PDF sets could be distinguished if relying only on the statistical error. Conversely, the plot on the right-hand side shows that with the current PDF uncertainties the predictions obtained with the various PDF sets cannot be disentangled. They are indeed all within the errors. The spectrum does not seem to be able to constrain the PDF of the valence/sea up-quark and the sea up-type antiquark. In Fig. 5 we show the equivalent picture for the case of the A * FB observable. In the left-hand side plot, we observe that the precision allowed by the available statistics is sufficient for a precise experimental measurement of the observable, reflecting what shown for the predicted di-lepton spectrum. In the plot on the right-hand side, we can appreciate the good features of the A * FB in terms of reduction of the PDF uncertainties, such that for this observable it will be possible to distinguish the predictions obtained with some of the PDF sets in specific invariant mass regions. In particular, the NNPDF prediction cannot be disentangled from the others just by observing the differential cross section in the di-lepton invariant mass, because of the large PDF error which the spectrum is affected by. In the A * FB distribution instead the NNPDF central value and error band come out to be incompatible with those of the other PDF sets. A new sensitivity to the PDFs is then gained by associating the measurement of A * FB to the spectrum.

V. CONCLUSIONS
In summary, building upon the recent observation made in [5] that the forward-backward asymmetry of di-lepton pairs produced via NC DY processes at the LHC can be used to constrain the PDFs inside the protons, we emphasise here that its exploitation in conjunction with a cut on the rapidity of the final state gives access also to the flavour content of the proton. In the presence of such a cut, one could indeed separate the valence/sea up-quark and the sea up-type antiquark, as the y-cut enforces a suppression of the dd-quark luminosity. Depending on the size of the data sample, the statistical and PDF errors for the different PDF sets could be not overlapping, thus offering the possibility to constrain or even eventually measure the up-type quark/antiquark PDFs.
This should be intended as a preliminary analysis, quite promising though. Higher order EW and QCD corrections can become quite important in the kinematical region of high rapidity. They should be properly included in order to make solid statements. Precise and reliable Monte Carlo Event Generators are available to this purpose and are regularly used by the LHC experimental collaborations for data interpretation. Having these tools at hand, it would be highly advisable to perform a new PDFs fitting procedure aimed at establishing, in a firm way and by direct comparison, whether the use of the reconstructed forward-backward asymmetry can give better results in terms of PDFs constraints with respect to the triple differential cross sections now in use.