Partonic transverse momenta in soft collisions

The partonic transverse momentum, kt\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$k_t$$\end{document}, distribution plays a crucial role in driving high-energy hadron interactions. If kt\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$k_t$$\end{document} is limited we have old fashioned Regge theory. If kt\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$k_t$$\end{document} increases with energy the interaction may be described by perturbative QCD. We use BFKL diffusion in lnkt\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\ln k_t$$\end{document}, supplemented by a stronger absorption of low kt\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$k_t$$\end{document} partons, to estimate the growth of the mean transverse momenta ⟨kt⟩\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\langle k_t\rangle $$\end{document} with energy. This growth reveals itself in the distribution of secondaries produced at the collider energies. We present a simple, BFKL-based, model to demonstrate the possible size of the effect. Moreover, we propose a way to evaluate experimentally the shape of the parton transverse momenta distribution by studying the spectra of the (D\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$D$$\end{document} or B\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$B$$\end{document}) mesons which contain one heavy quark.


Introduction
Contrary to old Regge theory, where it was assumed that the transverse momenta of all the particles are limited, QCD is a logarithmic theory where there is a possibility that the parton's (quark, gluon) transverse momentum may increase during the evolution. In particular, already in leading order (LO) BFKL evolution there is diffusion in ln k t space [1,2]. From the experimental point of view, it is relevant to note that the growth of the mean transverse momenta, p t , of secondary hadrons with collider energy was observed at the Tevatron and at the LHC (see e.g. [3,4]). In order to describe this growth in DGLAP-based Monte Carlo generators [5,6] an additional infrared cutoff, k min , was introduced. Of course, in any case, we need a cutoff to avoid the infrared divergence of the amplitude of the hard (parton-parton interaction) subprocess. However, at first sight, we would expect this cutoff to have its origin in confinement. It should be less than 1 GeV and should not depend on energy. On the contrary, it turns out that to reproduce the energy dependence of the data, the value a e-mail: a.d.martin@durham.ac.uk of k min should increase as k min ∝ s 0.12 [5]; such that at the Tevatron energy k min 2 GeV, while at the LHC k min 3 GeV.
In Sect. 2 we present a simple model which accounts for BFKL ln k t diffusion, together with the absorptive effects which additionally suppress the low k t partons, since the absorptive cross section behaves as σ abs ∝ 1/k 2 t . That is, we now have a dynamically induced infrared cutoff. 1 In Sect. 3 we use this model to obtain the expected energy and rapidity dependence of k t distributions. In Sect. 4, we discuss the possibility to directly study these effects experimentally by measuring the p t spectra of D (and/or B) mesons. Due to the strong leading particle effect (see e.g. [11][12][13]) the transverse momentum of mesons which contain a heavy quark is close to the transverse momentum of the heavy quark. Moreover, final-state interactions and confinement do not appreciably distort the original distribution of these heavy mesons.

BFKL-based model
The original BFKL equation [14][15][16][17] may be written as an integral equation for the gluon distribution unintegrated over k t , in the form where the kernel is evaluated as The first term in the kernel 2 can be understood as the effect of the emission of a daughter gluon with momentum (x, k t ) from a parent gluon with momentum (x = x/z, k t ). This generates the ladder structure of the pomeron sketched in Fig. 1a. The remaining two terms in the kernel (depending on f (x/z, k t )) account for the loop corrections which occur in the trajectory of t-channel reggeised gluons.
It is important to evolve in k t (as well as x) to be able to understand the origin and the behaviour of the dynamical infrared cutoff-that is, to see how the k t (s, y) distribution is generated within perturbative QCD. Here y = ln(1/x) is the rapidity of the parton. This dynamically generated cutoff affects (i) the p T distribution of secondary hadrons, (ii) the slope, α P , of the (QCD) pomeron trajectory and (iii) the values of the triple-and multi-pomeron couplings which control the predictions of the cross sections for diffractive dissociation.
Equation (2) can be solved numerically, starting, for example, from an input gluon with where we take x 0 = 0.2. Since the probability to have a large k t gluon should be suppressed by the small QCD coupling, we have included in (4) the factor α s (k t ). We use the oneloop running coupling α s (k 2 t ) with QCD = 0.15 GeV, and the number of light quarks to be n f = 4. Besides this, we account for the simple kinematical constraint-when the parton carries fraction x of the initial proton momenta, the transverse momentum k t cannot exceed the value k t,max = √ xs, where √ s is the initial energy. It is natural to approximate the input by taking x 0 = 0.2. The reasons are as follows. For BFKL evolution we have to consider small x, but we would like to cover the largest possible rapidity interval. Therefore we start with x 0 = 0.2, reserving a larger x interval for the valence quarks, and possible Good-Walker diffractive eigenstates [18], which describe low-mass diffractive dissociation. Moreover the typical DGLAP input gluon has a (1 − x) 5 type distribution corresponding to a mean of x of about 0.2.
To include the effects of absorption (that is, the rescattering of intermediate partons along the ladder) we follow [19] and multiply the BFKL kernel K of Eq. (3) by a canonical absorptive factor of the form exp(−λ (y, k t )/2) which depends on the rapidity, y = ln(1/x) and the k t of the current parton. Here is the optical density of the target gluon, while the factor λ accounts for the value of the triple-pomeron vertex, such that λ is the opacity of an incoming proton-'current' parton interaction. 3 However, we must account for the absorption by both the incoming beam (a) and the target (b) protons interacting with intermediate partons. That is, actually the absorptive factor reads where y (y ) is the rapidity difference between the beam (target) proton and the current, intermediate gluon in the BFKL evolution. Denoting the rapidity separation between the beam and the target protons by Y , we have y = Y − y.
The simplest absorptive effect comes from the triplepomeron diagram shown in Fig. 1b. As in [19], we use the Leading Log expression for the BFKL triple-pomeron vertex, that is [8,20,21], The -function reflects the fact that (after averaging over the azimuthal angle) the large-size pomeron (i.e. the ladder with small k t ) does not 'see' the small-size colourless object described by the BFKL pomeron component with k t > k t . Note that the suppression factor, written in the form (5), includes not just the triple-pomeron diagram, but also a series of the multi-pomeron contributions generated by the vertices, g n m , coupling n to m pomerons. Here we prefer to take the simple eikonal-like expression which satisfies the AGK cutting rules [22,23]. However, this means that we have to replace the exponents exp(−λ /2) in (5) by the factor So now the absorptive factor (5) becomes with the λ(k t , k t ) given by (6).
In terms of gluon density f (x, k t ), the 'differential' opac- where B g /2 is the t-slope of the initial "constituent gluon" form factor; we take 5 B g = 1 GeV −2 . To obtain the full opacity we take the integral where the lower limit reflects the (k t − k t ) function in (6).
Since the opacity a (y, k t ) is proportional to f a (y, k t ) we may write the evolution equation in rapidity y, just in terms of opacities. Thus, finally, we obtain a system of two evolution equations. One equation evolving for b up from the target (b) at y = 0, and one for a evolving down from the beam (a) at y = Y k = ln(s/k 2 t ). That is, 4 This equation follows after integrating Eq. (17) of [19] over the impact parameter, b, or from [20]. 5 There are several arguments in favour of the effective slope B g being of the order of 1 GeV −2 ; that is, in favour of the small-size 'hot-spot' transverse area occupied by our gluon amplitude. The first reason, is the small radius of the gluonic form factor of the proton calculated using QCD sum rules [24]. The next argument is the small value of the effective slope of the pomeron trajectory observed experimentally. Further evidence is the success of the additive quark model, σ (π p)/σ ( pp) 2/3. Finally, in the explicit calculation of our amplitude, following [25], we indeed found an almost constant slope B g 0.9 GeV −2 for the present collider energy interval.
for the evolution of gluon distributions from both the target and the beam initial hadrons (protons) in the absorptive (background) field of both hadrons. This system can be solved by iteration. In fact, it converges after just a few iterations.

The parton k t distribution
The transverse momentum distribution at rapidity y has the form dσ The system of Eq. (12) was solved numerically by iteration, introducing an infrared cutoff k 0 = 0.5 GeV; that is, assuming f (y, k t < k 0 ) = 0. The resulting transverse momentum distributions are presented in Fig. 2. The solid lines are the predictions for the gluon distribution in the central plateau region (with rapidity y = Y/2), while the dashed lines correspond to distributions shifted to the fragmentation region of the incoming proton (i.e. initial gluon) with y = Y/6. The y = Y/6 curves are steeper and the corresponding mean transverse momentum is smaller than that in the centre (y = Y/2). As expected the distributions become flatter when the energy increases. However, at the Tevatron (thick black curves) and even at the 8 TeV LHC (thin black curve) we are still far from true saturation. Only at √ s = 100 TeV do we predict an horizontal interval for k t < 2 GeV. For  TeV if one neglects the absorptive effects, that is, we have the case when the survival factor S ≡ 1 in (12). The distribution then decreases approximately linearly with increasing k t .
Next in Fig. 3 the rapidity dependence of mean transverse momenta, k t is shown for the Tevatron ( √ s = 2 TeV) and the LHC ( √ s = 8 TeV), and for √ s = 100 TeV. The value of mean k t increases with energy, and it decreases as the rapidity approaches the position of the initial hadron.

How to measure the k t distribution
Note that the predicted values of k t are of the same 'order of magnitude', but smaller than, the value of the k min cutoff used in the PYTHIA Monte Carlo, which is based on DGLAP evolution. However, we have to recall that (a) these are not exactly the same quantities and (b) here we have used a simplified model based on the LO BFKL kernel. 6 The advantage of this model is that it is sufficiently transparent and practically has no free parameters. The only exceptions are the starting values of x 0 = 0.2 and k 0 = 0.5 GeV and the slope B g = 1 GeV −2 of the initial 'constituent' gluon. The parameters are not chosen to describe the data, but simply taken to have physically reasonable values. Besides this, there may be some 'intrinsic' transverse momentum of the initial gluon which will enlarge the final value of k t .
We should emphasise that the partonic k t distribution, although not directly observable, drives all soft high-energy interactions. Clearly it would be interesting to measure the gluon's k g,t experimentally. Can this be done? The problem is that actually we never observe partons, but only the final secondary hadrons, which are mainly pions. Unfortunately the distributions of light hadrons (such as pions, kaons) are strongly affected by final-state interactions: that is, by hadronisation, confinement and the decay of resonances. In particular, the p t distribution of secondary pions strongly depends on the possible colour re-connection. Therefore it appears better to study the distributions of mesons which contain one heavy quark. Due to the strong leading particle effect [11][12][13], the p t distribution of these mesons is close to that of the heavy quark. Since heavy quarks are mainly produced by the gg → QQ subprocess, we may expect that (modulo some smearing due to hadronisation when the heavy quark picks up a light antiquark) the mean momentum of such a meson should carry the momentum of the parent gluon. Final-state interactions and resonance decays do not appreciably distort the p t distributions of these heavy mesons.
On one hand, it might be the best to measure the p t distributions of heavy B-mesons, where the leading particle effect is more pronounced. On the other hand, the b-quark already receives a rather large from the hard gg → bb subprocess and it may be hard to observe the variation of the incoming gluon k g,t on the top of this large 'background', k background . Therefore, it seems better to detect D-mesons where the value of k background ∼ m c is comparable with the expected gluon's k g,t . We would hope to observe the growth of p D,t with energy at fixed rapidity, and a decrease of p D,t with pseudorapidity 7 at a fixed energy. The last effect can be observed by a comparison of the CMS/ATLAS data at η = 0 with the LHCb data at η = 3-4. Moreover, note that it possible to do better. We could suppress the k background contribution generated into the 'hard' gg → QQ subprocess if the transverse momenta of both heavy mesons (D andD or B andB) are measured. In such a case the transverse momentum of the QQ pair is simply equal to the momentum of the parent gluon pair. Of course, we cannot avoid the smearing due to hadronisation, but it is not so large since it is controlled by the confinement scale and not by the heavy quark mass. So it would be good to measure the vector sum of the momenta of the two heavy mesons, or just the coplanarity between the two heavy mesons.
Non-complanarity should increase with energy, but decrease with η.
Another attractive measurement is to compare the p t of the secondaries produced in the diffractive dissociation with those from non-diffractive inelastic events. It is usually expected (see, for example, [26,27]), that the spectra of particles produced in proton diffractive dissociation into a highmass (M X ) state, are similar to that in normal inelastic events taken at an energy √ s = M X . That is in the situation when the energies of the final states are the same. On the contrary, in the picture described above, even in the case of dissociation, the p t distribution of secondaries should be driven by the parton's k t formed by the whole initial energy √ s M X . That is, it does not matter whether the events have a large rapidity gap (LRG) or not. One consequence (see, also, [28]) is that in proton diffractive dissociation to a large M X system (but still M X √ s) the dissociation events, especially near the edge of the LRG, are expected to have a larger p t than those in a normal inelastic pp-collision at √ s = M X ; modulo to possible hadronisation effects. Moreover, the rapidity dependence of the p t spectra in LRG events are also similar to that in the inelastic interaction at full proton-proton energy √ s, and not to the inelastic events with the proton-proton energy equal to M X . Again, to reduce the effects of hadronisation, it would be better to make the comparison by measuring the distributions of D-mesons both in inelastic and high-mass dissociation events.

Conclusions
The transverse momentum distribution of partons plays a pivotal underlying role both in the spectral shape of secondaries and in the asymptotic behaviour of high-energy proton-proton collisions. At first sight, just from dimensional arguments, we expect dσ/dk 2 t ∝ 1/k 4 t . That is, the major contribution should come from low k t , close to the cutoff ( 0.3 GeV) provided by confinement. On the contrary, to describe the data, it was necessary to introduce a much higher cutoff, k min , in the hard matrix element of the order of a few GeV, with a value that increases with collider energy, like s 0.12 . Actually such a k min was obtained by tuning the Monte Carlo generators [5,6], but clearly it should be of theoretical origin. Moreover, k min of the order of a few GeV should be explained in terms of perturbative QCD.
Here, we use a model based on the LO BFKL equation, supplemented by absorptive multi-pomeron corrections. The original BFKL equation includes diffusion in logk t , with, at each step of the evolution, the possibility that k t may increase or decrease with equal probabilities. However, strong absorption of low k t partons leads to a growth of k t with collider energy. We demonstrate that this effect naturally explains the observed energy behaviour of the effective cutoff, k min .
We did not perform a fit to the data, but show, at a qualitative level, that a simplified model based on leading order perturbative QCD with a few physically motivated parameters, produces a reasonable k t distribution of the partons. We present the expected k t distributions at different collider energies and the dependence of k t on the energy and rapidity of the parton.
Although the k t of the parton is not directly observable, we discuss the possibility to experimentally verify these predictions. One way, is to measure the p t distributions of mesons containing a heavy c or b quark, or better to measure DD or BB meson pairs. Another possibility is to compare the p t spectra of diffractive dissociation events with those of nondiffractive inelastic scattering.