Testing the parton evolution with the use of two-body final states

We consider the production of bb¯\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$b{\bar{b}}$$\end{document} quarks and Drell–Yan lepton pairs under LHC conditions focusing attention on the total transverse momentum of the produced pair and on the azimuthal angle between the momenta of the outgoing particles. Plotting the corresponding distributions in bins of the final-state invariant mass, one can reconstruct the full map of the transverse momentum dependent parton densities in a proton. We give examples of how these distributions can look like at the LHC energies.

Experiments of new generation running at the LHC yield plenty of high precision data. In order to properly interpret these data we need the parton distribution functions to be known with adequately good accuracy. This, in turn, raises the question of a detailed measurement of the parton distributions. In this note we focus attention on two important kinematic observables which enable us to reconstruct the full map of the transverse momentum dependent (TMD), or unintegrated, parton densities. We address the LHC conditions ( pp collisions at √ s = 7 TeV), for which we give a number of illustrations.
The evolution of TMD gluon densities can be explored with the production of bb pairs. At the LHC energies, this process is dominated by the direct leading-order (LO) offshell gluon-gluon fusion subprocess while the contribution from the quark-antiquark annihilation is of almost no importance because of the comparatively low quark densities. The four-momenta of corresponding particles are given in the parentheses. The present calculation of the process (1) is fully identical to that performed previously [1]. The evolution of TMD quark densities can be explored with the production of Drell-Yan lepton pairs. This a e-mail: lipatov@theory.sinp.msu.ru process is dominated by the off-shell quark-antiquark annihilation subprocess where q is for the valence and sea quarks andq stands for the sea anti-quarks. The present calculation of the process (2) is fully identical to that from [2]. We do not consider here higher-order corrections q +q → l + + l − + g since they are already taken into account in the k T -factorization approach [3][4][5][6][7] as a part of the evolution of TMD quark densities. The final states of the processes (1) and (2) are represented by two-body systems with fully reconstructible kinematics where the transverse momentum p T of the bb or lepton pair measures the net transverse momentum of the initial gluons or quarks, the invariant mass of the pair measures the product of longitudinal momentum fractions, M 2 = x 1 x 2 s, and the rapidity of the pairs measures the ratio of the momentum fractions, y = (1/2) ln(x 1 /x 2 ). A useful complementary observable is the difference between the azimuthal angles of produced particles φ. In the LO of collinear QCD factorization, the p T and φ distributions degenerate into delta functions at p T = 0 and φ = 0, and the continuous spectra can only be obtained by including higher-order corrections. In the k T -factorization approach, these radiative corrections are automatically taken into account in the form of TMD parton densities. Comparing the p T and φ spectra at varying gluon momentum fraction x we watch the evolution of parton distributions.
To simulate the bb pair production we used the latest JH'2013 parametrization [8] for the TMD gluon densities in a proton. The input parameters of this gluon distribution were fitted to describe the proton structure function F 2 . To simulate the production of Drell-Yan lepton pairs we applied complementary TMD valence quark distributions from the same set [8]. The necessary TMD sea quark densities are calculated from the gluon ones in the approximation where the sea quarks occur in the last gluon-to-quark splitting [9]. The results of our calculations are displayed in Figs. 1, 2, 3, 4, 5, 6, and 7. Shown in Figs. 1 and 2 are the spectra of bb pair and dilepton transverse momentum p T and the azimuthal angle φ plotted for several different intervals of their invariant mass M. Here, to make the changes in shape more easily recognizable, we show the normalized differential cross sections. We see that with increasing M the maximum in the p T spectrum shifts gradually to higher values, and the whole distribution becomes more flat. The φ distribution moves toward φ π , which is due to the inequality M p T .
The latter becomes even stronger at high M (see Fig. 3). As one can see from Besides the restrictions on the invariant mass, the special kinematical cuts on the final state give us further possibilities to find the region of x and/or partonic transverse momenta we desire. It is illustrated in Figs. 6 and 7, where we plot the normalized differential cross sections of the considered subprocesses calculated as functions of x and k 2 T (the longitudinal momentum fraction and transverse momentum of one of the colliding partons) with the additional cuts applied to the rapidity y of the final-state quark or lepton pair. As an example, we used y < 1 and 3 < y < 4. We show that under these cuts one can probe different x and/or k 2 T regions and extract information on the TMD parton distributions at the scale given by M. Note that the different k 2 T regions can be obtained under additional restrictions on the quark or lepton pair transverse momentum p T and/or azimuthal angle φ.
The TMD parton distributions in a proton can be calculated using the different approaches (see, for example, [10] and the references therein). Of course, their different behavior as a function of x and/or k 2 T is reflected in the description of the LHC data. For example, we compared the rapidity distribution of Drell-Yan lepton pair production calculated regarding the kinematical conditions imposed by the LHCb Collaboration [11] using the JH'2013 parton densities (as above) and the ones obtained from the Kimber-Martin-Ryskin (KMR) prescription [12,13] (see Fig. 8, left panel). The latter is a formalism to construct the TMD parton densities from the known conventional parton distributions. The key assumption is that the k T dependence enters at the last evolution step, so that usual DGLAP evolution [14][15][16][17] can be used up to this step. The rapidity distribution of Drell-Yan pair production is sensitive to the x-behavior of TMD partons. One can clearly see that the difference between the TMD parton densities applied (Fig. 8, right panel) leads to a different description of recent LHCb data [11]. Therefore, the LHC experimental data for the processes considered can be used to constrain the TMD parton distributions in a proton. The experimental data are from LHCb [11]. Right panel the TMD up quark distributions calculated as a function of x at k 2 T = 10 GeV and μ 2 = 100 GeV 2 . The solid and dashed curves on both panels correspond to the JH'2013 and KMR sets, respectively Additionally, we investigate the dependence of estimates presented above on the parton shower effects using the Monte Carlo event generator cascade [18]. As expected, we observe only a very small contribution of the initial-state parton shower, since in the k T -factorization approach it does not influence the transverse momentum of the gluons (because it is determined from the TMD gluon density).
Thus, we conclude that one can map the evolution of parton distributions at the scale M from high values of proton longitudinal momentum fraction x to low ones by applying different cuts on the final states. This is important for further precise determination of the TMD quark and gluon densities in a proton from the LHC data.