LU TP 1533
March 2016
The DrellYan process as a testing ground for parton distributions up to LHC
Eduardo Basso Department of Astronomy and Theoretical Physics,
Lund University, SE 22362 Lund, Sweden
Claude Bourrely Aix Marseille Université, Université de Toulon, CNRS, CPT, UMR 7332
13288 Marseille, Cedex 09, France
Roman Pasechnik Department of Astronomy and Theoretical Physics,
Lund University, SE 22362 Lund, Sweden
Jacques Soffer Physics Department, Temple University,
1925 N, 12th Street, Philadelphia, PA 191221801, USA
Abstract
The DrellYan massive dilepton production in hadronhadron collisions provides a unique tool, complementary to Deep Inelastic Scattering, for improving our understanding of hadronic substructure and in particular for testing parton distributions. We will consider measurements of the differential and doubledifferential DrellYan cross sections from FNAL Tevatron up to CERN LHC energies and they will be compared to the predictions of perturbative QCD calculations using most recent sets (CT14 and MMHT14) of parton distribution functions, as well as those provided by the statistical approach.
Key words: DrellYan process, parton distribution functions
PACS numbers: 12.40.Ee, 13.38.Bx, 14.70.Hp
1 Introduction
The Deep Inelastic Scattering (DIS) of leptons and nucleons is indeed our main source of information to study the internal nucleon structure in terms of parton distribution functions (PDF). However at hadron colliders, one has access to the DrellYan (DY) process [1] and the measurement of massive dilepton production, via exchange in the channel , is also an excellent, unique and clean observable for testing QCD and, in particular, the nucleon structure encoded in PDFs [2]. Needless to mention, the DY process is one of the standard candles for New Physics searches at the energy and luminosity frontiers which thus requires a significant reduction in theoretical uncertainties. This is the reason why several DY experiments were already performed in almost all highenergy hadronic facilities and they are now being intensively studied in all ongoing major experiments at the Large Hadron Collider (LHC) at CERN. This will allow us to test various theoretical models in a vast kinematical domain and to analyse the corresponding PDF uncertainties with this simple process whose mechanism is dominated by partonantiparton annihilation in collisions.
The rapidity and the invariant mass of the dilepton system produced in protonproton collisions are related, at leading order (LO), to the momentum fraction carried by the parton in the forwardgoing (backwardgoing) proton, according to the formula . Therefore, the rapidity and mass distributions are sensitive to the PDFs of the interacting partons. Besides, the transverse momentum distributions of the dilepton provides an additional information about the dynamics of proton collisions at high energy particularly sensitive to the higherorder QCD corrections. In this work, several most recent PDFs parameterizations at the nexttoleading order (NLO) such as MMHT14 [3] and CT14 [4] models, as well as the NLO statistical PDFs previously developed in Ref. [5], will be employed for description of the existing DY data for the differential , and distributions, both at low and high energies. For this purpose, we have selected the data from a limited number of experiments at FNAL Tevatron (E866, D0) and CERN LHC (CMS, ATLAS) keeping only those with the highest integrated luminosity, corresponding to the highest precision.
The paper is organized as follows. In Section 2, we review the main features of three sets of PDFs we have used for our calculations. In Section 3, we consider the invariant mass distribution at LHC (7/8 TeV) energies vs available data over a very broad dilepton mass range up to 2 TeV, and also in a much smaller mass range from a fixedtarget FNAL experiment. In Section 4, the differential cross section as a function of the rapidity is analyzed. In Section 5, we study the transverse momentum spectra for two different centreofmass energies. We give our final remarks and conclusions in Section 6.
2 PDFs selection
We will now summarize the essential properties of three sets of PDFs which will be tested in our analysis of the DY process.
First, let us recall the main features of the statistical approach [5] for building up the PDFs as opposed to the standard polynomial type parameterizations based on Regge theory at low and on counting rules at large . The fermion distributions are given by the sum of two terms, a quasi FermiDirac function and a helicity independent diffractive contribution
(1) 
(2) 
defined at the input energy scale . We note that the diffractive term is absent in the quark helicity distribution and in the quark valence contribution .
In Eqs. (1) and (2) the multiplicative factors and in the numerators of the nondiffractive parts of the ’s and ’s distributions, imply a modification of the quantum statistical form, which was proposed in order to agree with experimental data ^{1}^{1}1These factors were fully justified in the extension of the PDFs to include their transverse momentum dependence (TMD) [6].. The parameter plays the role of a universal temperature and are the two thermodynamical potentials of the quark , with helicity . They represent the fundamental characteristics of the model. Notice the change of sign of the potentials and helicity for the antiquarks^{2}^{2}2At variance with statistical mechanics where the distributions are expressed in terms of the energy, here one uses which is clearly the natural variable entering in all the sum rules of the parton model..
Although the statistical approach to the starting PDFs allows the simultaneous description of unpolarized cross sections and helicity asymmetries, a unique situation in the literature, in this work we will restrain ourselves to spinindependent DY observables. For a given flavor , the corresponding quark and antiquark distributions involve eight free parameters: , , , , , and . It reduces to since one of them is fixed by the valence sum rule, , where , respectively. In the light quark sector , the total number of free parameters is reduced to by applying additional constraints as was done in Ref. [7] (for a more detailed review, see e.g. Ref. [8])
For the gluon PDF at the input energy scale, we consider the blackbody inspired expression
(3) 
a quasi BoseEinstein function, with being the only free parameter, since is determined by the momentum sum rule. To summarize, this determination of all PDF sets^{3}^{3}3In Ref. [5] we have considered the helicity gluon distribution which is irrelevant in the present work involves a total of seventeen free parameters. Namely, in addition to the temperature and the exponent of the gluon distribution, we have eight free parameters for the light quarks , seven free parameters for the strange quarks as was outlined above. These parameters were determined in Ref. [5], from a NLO QCD fit of a large set of accurate DIS data . The resulting PDFs, denoted from now on as BS15, are illustrated for the energy scale in Fig. 1.
Next we consider an improved version of the socalled MSTW framework which was proposed seven years ago in Ref. [9]. The main shortcoming of the MSTW PDFs was an incorrect description of the lepton charge asymmetry from decays as a function of the lepton rapidity. The new input distributions result from major changes in the theoretical procedure with respect to original Ref. [9] since they involve Chebyshev polynomials. In this new version proposed in Ref. [3] the majority of the starting PDFs have the following form
(4) 
where is the input energy scale, and are Chebyshev polynomials in , with , where one takes . For each PDF, namely, one has to determine the parameters . Here, and denote the lightquark valence distributions and is the lightquark sea distribution. For one sets .
We still have to specify the parameterisations of the gluon and of the differences and . For one sets and one uses the following expression
(5) 
For the poorly determined strange quark difference one takes
(6) 
Finally, for the gluon distribution, as proposed long time ago [10], one needs a second term for the small behavior as shown below
(7) 
and we notice that it involves seven free parameters since is constrained by the momentum sum rule. This is a major difference with respect to the statistical approach. In total, there are thirty seven free parameters, a large number, while one should remember that there are three constraints from the valence sum rules as already mentioned above. In addition, there is also the strong coupling defined at the mass scale of the boson, i.e. , which is allowed to be free when determining the best fit. The authors claim that the advantage of using a parameterisation based on Chebyshev polynomials is the stability and a good convergence of the values found for the coeffients . In what follows, we refer to these PDFs as to the MMHT14 NLO model.
The parameters were determined from a global analysis of a variety of new data sets, from the LHC including some DY results, updated Tevatron data and HERA combined H1 and ZEUS data on the total and charm structure functions. For DIS data on deuterium targets deuteron corrections were taken into account, as well as nuclear corrections for neutrino data. The NLO and a nexttonexttoleading order (NNLO) QCD fits of the data were performed and we display the resulting PDFs in Fig. 2 at the same energy scale as in Fig. 1. The fit quality and the determination of the PDF uncertainties were carefully studied in Ref. [3]. It is interesting to observe that in most cases the new MMHT PDFs are within one standard deviation from those of the MSTW framework [9], a remarkable stability of the procedure.
We now turn to the third PDF set we will use in our calculations of the DY cross sections. It results from different versions [11] up to NNLO, from the CTEQTEA global analysis of QCD, including the use of data from HERA, Tevatron and LHC. New theoretical developments are associated to the achievement of an increasing data precision in DIS, vector boson production and singlejet production. It is interesting to note that DY data from colliders were not considered so far. Standard parametrizations for each flavor are of the form
(8) 
where the first two factors are suggested by the Regge theory and by counting rules, and the remaining factor is assumed to be slowly varying. In the previous CTEQ analyses, for each flavor was chosen as an exponential of a polynomial in or , for instance,
(9) 
for valence quarks .
In their most recent work refered to as CT14 [4], for the valence quarks they reexpress the polynomial as a linear combination of Bernstein polynomials in
(10) 
where . Then seven parameters for each flavor are reduced to just four by setting , and using the valence sum rule. So the valence quark and are determined in terms of a total of eight free parameters.
The CT14 model uses a similar parameterisation for the gluon but with a polynomial of a lower order since the data provide fewer constraints on the gluon distribution
(11) 
where . However, instead of , it employs . The momentum sum rule reduces the total number of parameters of the gluon distributions down to five. The sea quark distributions and were parametrized using fourthorder polynomials in with the same variable that was used for the gluon. They assumed at , which implies .
All in all, the CT14 model has eight free parameters associated with the valence quarks, five parameters associated with the gluon, and thirteen parameters associated with sea quarks, which in total amounts to twenty six fitting parameters. The initial energy scale =1.295 GeV has been used to perform the NLO and NNLO QCD fits of a large set of data. We display the resulting PDFs in Fig. 3 at as a reference.
All the PDFs considered in this work were obtained from their correspondent LHAPDF6 [12] grids, which were used in our MonteCarlo simulations described below^{4}^{4}4The LHADPF grids for the BS15 NLO model shall appear in online database soon..
3 Dilepton invariant mass distribution
Consider first, as displayed in Fig. 4, the measurement of DY cross section at the LHC energy TeV allowing to get dilepton mass up to 2 TeV.
The NLO cross section at LHC energies was calculated using the subtraction method implemented in the DYNNLO code [14]. This semianalytical tool employs the dipole subtraction formalism of Catani and Seymor [15] realised by using the MCFM event generator [16] and stands as one of the uptodate calculations of the DrellYan process up to NNLO. In the case of the invariant mass distributions the full phase space is integrated within each bin in the dilepton mass . Note, in all the calculations here and below the factorisation and renormalisation scale are taken to be equal to the transverse mass of the dilepton pair, by convention. One notices that the theoretical results with distinct PDFs behave similarly for a broad range of the invariant mass with minor differences.
The DY low mass region was specifically measured by the ATLAS Collaboration at luminosities of 1.6 fb and 35 pb [17], for which a comparison with three distinct NLO PDFs based on the DYNNLO implementation is shown in Fig. 5. The data used the invariant mass region of 26 – 66 GeV for the higher luminosity whereas less precise data extend this region down to 12 GeV. Here, one starts to observe a minor deviation between predictions of BS15 and the other two PDFs. While the results of MHHT and CT models are closer to data, the BS15 predictions are still within the experimental error bars.
We also selected the data for the high mass DY electronpositron pair production measured by the ATLAS Collaboration [18] which correspond to 4.9 fb luminocity within the mass range of GeV and in the kinematics defined by dielectron pairs with GeV and . In Fig. 6 the high mass distribution of pairs is compared with the corresponding NLO pQCD predictions. The CT14 and MMHT14 results are pretty close to each other and show a very good agreement with data, while the BS15 result is somewhat below the data but is still within the error bars.
At much lower energies at FNAL = 38.8 GeV the E866/NuSea Collaboration provided the DY differential distributions over a broad kinematics range [20]. In this case, we have implemented the formalism by Sutton et al. [19] where both NLO corrections for Compton and annihilation subprocesses were included. This formalism results in a fair description of the E866/NuSea data on DY differential distributions in the protonproton collisions for three PDF models (BS15, CT14 and MMHT14) at various and dilepton masses shown in Fig. 7. Here, an overall uncertainty w.r.t. scale variations within the interval is shown as filled error bands. Note, the uncertainties grow with similarly in all three models. Noticeable deviations between the predictions from distinct PDF sets emerge mostly at small invariant masses GeV where some deficiency of the predictions can be observed. Specifically for this case, one observes that the CT14 PDF is rather unstable at and the CT14/BS15 PDFs do not work well at , while MMHT14 demonstrates a better description in both normalisation and shape. For intermediate values of the Feynman variable and large all the considered PDFs provide a fairly good description of the low energy data.
4 boson rapidity distribution
The events and analyses for this and further sections were generated with Pythia v8 [21] and Rivet [22]. For the boson rapidity distribution the standard LO Pythia (with parton showers and NLO PDFs) was employed. The parton shower within Pythia effectively resums the higherorder real corrections due to emissions off the colored initial states. The latter corrections significantly affect the kinematics of the final states and is need for a proper treatment of the differential distributions. Therefore, this approach differs from the full NLO calculation by the virtual corrections only. However, these corrections can only affect the overall normalisation of the cross sections and are typically accounted for by a universal factor which is practically independent on kinematics of the final states. Since the experimental data for the DrellYan process at Tevatron and LHC are typically provided in terms of ratios of the differentialtototal cross sections, one does not need to explicitly compute the virtual corrections since they cancel out in these ratios and do not affect our conclusions. Besides, other model uncertainties such as the scale uncertainty mostly cancel in such ratios as well. So it is not surprising that this method leads to a very good agreement of various PDF models with the Tevatron and LHC data on the ratios in the whole kinematic region (see below).
For validation of the considered PDF models and theoretical methods, we display in Fig. 8 the results on the rapidity distribution from the CMS Collaboration at = 7 TeV [23] corresponding to the integrated luminosity of 4.5(4.8) fb. The dimuon channel was used for rapidity reconstruction within certain invariant mass intervals, from which we have selected the one corresponding to a vicinity of the peak, GeV. Analogically, the rapidity distributions from the D0 Collaboration at a lower energy = 1.96 TeV [24] in the dielectron channel and the mass range 71 111 GeV corresponding to an integrated luminosity of 0.4fb were compared to the theoretical predictions with BS, MHHT and CT models. The corresponding results are shown in Fig. 9. All the NLO pQCD predictions exhibit a fair agreement with the existing data from both LHC and Tevatron measurements.
5 boson transverse momentum spectrum
A highprecision data on the boson transverse momentum distribution up to 800 GeV has been presented by ATLAS in Ref. [25] at = 7 TeV corresponding to the integrated luminosity of 4.7 fb. In Fig. 10 we consider the full detector acceptance whereas the boson transverse momentum for three different rapidity bins is shown in Fig. 11. A good description was obtained as well including rather high region where in principle the higher order corrections could start to play an important role, as well as the statistical fluctuations for the forward region () could become noticeable due to a smaller number of events generated in the forward/high kinematics. A disagreement of the model predictions in the highest bins, besides poor statistics, could be due to the fact we employ the simplistic analysis with LO matrix elements and NLO PDFs (plus parton shower) and not the fullprecision higherorder matrix elements. Remarkably enough, one can see that the PDFs tested here and our simplistic LO analysis give a fair description of data for almost entire range of boson ’s.
It is worth to present similar results for the boson transverse momentum distribution from the D0 Collaboration at Tevatron [26] in collisions at =1.96 TeV in the dielectron channel corresponding to the integrated luminosity of 0.98 fb. In Fig. 12 the boson transverse momentum distribution is plotted in two regions, the one that covers the full detector acceptance and the one which is focused on forward production, i.e., for and GeV. Similarly to the LHC case, the theoretical modelling with the chosen PDFs exhibit rather good description of the Tevatron data within the experimental uncertainties.
6 Concluding remarks
To summarize, the NLO pQCD predictions with the MHHT14, CT14 and BS15 PDFs studied in this work have resulted in a fair description of a broad range of DY data from FNALNUSEA energies up to LHC energies while a few minor deviations, mainly at the edges of the respective phase space, have been observed. Given the fact that three very different PDF sets lead to rather similar results we conclude that the NLO BS15 model having a much fewer free parameters which are fitted to the DIS data only is not worse than the other most recent parameterisations and should be considered on the same footing as the current global PDF fits when it comes to the DY process at high energies.
Extrapolating to the DrellYan process at even larger energies of LHC Run II ( TeV) we do not expect more significant differences between the BS15 predictions and those of other PDF models since typical values we probe there get even smaller at higher energies. In Fig. 13 we present (anti)quark PDF CT14toBS15 ratios (central values) as functions of and GeV. This figure demonstrates that the differences between the considered PDFs becomes smaller at larger and typical for DrellYan at large at the LHC; inclusion of error bars does not change this picture. Similar situation holds for comparison of BS15 with the MHHT14 PDFs. As an illustration, in Fig. 14 (left) we show the transverse momentum distributions of the boson at TeV for BS15, CT14 and MHHT14 PDFs. A noticeable deviation of the BS15 prediction, which lies below the other PDFs, is seen only at low GeV, otherwise differences between the predictions for this and other observables such as and distributions are small. We conclude that the BS15 NLO PDF is a good tool to investigate pQCD physics in the second run of LHC measurements.
As is seen in Fig. 7, at low energies the PDF models exhibit more substantial differences in shapes of the invariant mass and distributions than those at high energies. This fact indicates the critical importance of low energy DrellYan measurements such as the FNAL E906/SeaQuest experiment [27]. In particular, the E906 experiment will determine the ratio at large (up to ). The statistical BS15 model predicts at large as illustrated in Fig. 14 (right) together with CT14 and MMHT14 predictions where the filled bands correspond to scale variations within GeV interval for each PDF model. This characteristic effect is based on the Pauli exclusion principle and the fact that in the proton there are two quarks and one quark such that more compared to . A high precision measurement of such a ratio would therefore become a crucial test of physics behind the PDF models.
Acknowledgments
Useful discussions with Gunnar Ingelman, Jesper Roy Christiansen and Torbjorn Sjöstrand are gratefully acknowledged. E. B. is funded by CAPES and CNPq (Brazil), contract numbers 2362/139 and 150674/20155. R. P. is supported by the Swedish Research Council, contract number 6212013428.
References

[1]
S. D. Drell and T.M. Yan, Phys. Rev. Lett. 25, 316 (1970);
Ann. Physics 66, 578 (1971).  [2] J. C. Peng and J. W. Qiu, Prog. Part. Nucl. Phys. 76, 43 (2014).
 [3] L. A. HarlandLang, A. D. Martin, P. Motylinsky and R. S. Thorne, Eur. Phys. J. 75 5, 204 (2015).
 [4] S. Dulat et al., arXiv:1506.07443v2 [hepph].
 [5] C. Bourrely and J. Soffer, Nucl. Phys. A 941, 307 (2015).
 [6] C. Bourrely, F. Buccella and J. Soffer, Int. J. Mod. Phys. A 28, 1350026 (2013).
 [7] C. Bourrely, F. Buccella and J. Soffer, Eur. Phys. J. C 23, 487 (2002).
 [8] C. Bourrely, J. Soffer and F. Buccella, Eur. Phys. J. C 41, 327 (2005).
 [9] A. D. Martin, W. J. Stirling, R. S. Thorne and G. Watt, Eur. Phys. J. 63, 189 (2009).
 [10] A. D. Martin, R. Roberts, W. J. Stirling and R. S. Thorne, Eur. Phys. J. 23, 73 (2002).

[11]
P. M. Nadolski et al., Phys. Rev. D 78, 013004 (2008);
J. Pumplin et al., Phys. Rev. D 80, 014019 (2009);
H. L. Lai et al., Phys. Rev. D 82, 074024 (2010).  [12] A. Buckley et al., Eur. Phys. J. C 75 132 (2015).
 [13] V. Khachatryan et al. (CMS Collaboration), Eur. Phys. J. 75, 147 (2015).

[14]
S. Catani, L. Cieri, G. Ferrera, D. de Florian, M. Grazzini, Phys. Rev. Lett. 103 082001 (2009);
S. Catani, M. Grazzini, Phys. Rev. Lett. 98 222002 (2007). 
[15]
S. Catani and M. H. Seymour, Phys. Lett. B 378 (1996) 287;
Nucl. Phys. B 485 (1997) 291 [Erratumibid. B 510 (1998) 503].  [16] J. Campbell, R.K. Ellis, MCFM  Monte Carlo for FeMtobarn processes, http://mcfm.fnal.gov.
 [17] G. Aad et al. (ATLAS Collaboration), JHEP 1406, 112 (2014).
 [18] G. Aad et al. (ATLAS Collaboration), Phys. Lett. B 725 223 (2013).
 [19] P. J. Sutton, A. D. Martin, R. G. Roberts and W. J. Stirling, Phys. Rev. D 45, 2349 (1992).
 [20] J. C. Webb et al. (FNAL E866/NuSea Collaboration), arXiv:hepex/0302019; arXiv:hepex/0301031.
 [21] T. Sjöstrand et al., Comput. Phys. Commun. 191, 159 (2015).
 [22] A. Buckley et al., Comput. Phys. Commun. 184, 2803 (2013).
 [23] S. Chatrchyan et al. (CMS Collaboration), JHEP 1312, 030 (2013).
 [24] V. Abazov et al. (D0 Collaboration), Phys. Rev. D 76, 012003 (2007).
 [25] G. Aad et al. (ATLAS Collaboration), JHEP 1409, 145 (2014).
 [26] V. Abazov et al. (D0 Collaboration), Phys. Rev. Lett. 100, 102002 (2008).
 [27] P. E. Reimer [Fermilab SeaQuest Collaboration], J. Phys. Conf. Ser. 295, 012011 (2011).