Probing the gluon density of the proton in the exclusive photoproduction of vector mesons at the LHC: A phenomenological analysis

The current uncertainty on the gluon density extracted from the global parton analysis is large in the kinematical range of small values of the Bjorken - $x$ variable and low values of the hard scale $Q^2$. An alternative to reduces this uncertainty is the analysis of the exclusive vector meson photoproduction in photon - hadron and hadron - hadron collisions. This process offers a unique opportunity to constrain the gluon density of the proton, since its cross section is proportional to the gluon density squared. In this paper we consider current parametrizations for the gluon distribution and estimate the exclusive vector meson photoproduction cross section at HERA and LHC using the leading logarithmic formalism. We perform a fit of the normalization of the $\gamma h$ cross section and the value of the hard scale for the process and demonstrate that the current LHCb experimental data are better described by models that assume a slow increasing of the gluon distribution at small - $x$ and low $Q^2$.

in the global analysis of the experimental data. We present the results for two different values of the hard scale: Q = M J/ψ /2 (left panel) and MΥ/2 (right panel). M V is the mass of the vector meson and √ s is the center -of -mass energy of the hadron -hadron collision. This property was the main motivation for the proposition presented in Ref. [28], which was improved by several authors in the last fourteen years [29][30][31][32][33][34]. Our goal in this paper is to complement these previous studies by performing a phenomenological analysis of the exclusive vector meson photoproduction which can illuminate several aspects of the formalism and about the current parametrizations for the gluon distribution obtained in the global parton analysis. Basically, we consider different models for xg and estimate the cross section for the photoproduction of vector mesons in γp interactions using the leading logarithmic formalism [14,15]. We demonstrate that assuming the usual value of the factorization scaleQ = M V /2, the HERA and LHCb data for the J/Ψ photoproduction are not described by these models. Taking into account that this formalism have been derived at leading order, which implies a uncertainty in the choice ofQ, we determine its value by fitting the γh data. Using these best fit parameters we estimate the rapidity distributions for the J/Ψ, Ψ(2S) and Υ production in pp collisions at the LHC and compare with the recent LHCb data [22][23][24]. Our results indicate that the current experimental data for the exclusive vector meson production in pp collisions can only be described by models which predict a slow increasing of the gluon distribution at small -x.
Lets initially present a brief review of the formalism for the calculation of the vector meson production in pp collisions. The process is represented in Fig. 2, with its rapidity distribution being given by where we have taken into account that the two incident protons can be the source of the photon and ⊗ characterizes the presence of a rapidity gap in the final state. The photon -hadron center of mass energies squared W 2 ± and the photon energies ω ± are given by where √ s h1h2 is the hadron -hadron center -of -mass energy, M V is the mass of the vector mass and y its rapidity. Moreover, the factors S 2 (W ± ) characterizes absorptive corrections which can destroy the rapidity gaps generated in exclusive processes [33,34]. The elastic photon flux for the proton can be expressed by [54] with the notation Ω = 1 + [ (0.71 GeV 2 )/Q 2 min ] and Q 2 The main input for the calculation of the rapidity distribution [ Eq. (1)] is the photon -hadron cross section for the vector meson production σ γh→V h . In the leading logarithmic approximation, the cross section for the vector meson production off any hadronic target, including a nucleus, at small-x and for a sufficiently hard scale, is proportional to the square of the gluon parton density of the target. To lowest order, the γh → V h (h = p, A) amplitude can be factorized into the product of the γ → qq transition (q = c, b), the scattering of the qq system on the target via (colorless) two-gluon exchange, and finally the formation of the quarkonium from the outgoing qq pair. The heavy meson mass M V ensures that perturbative QCD can be applied to photoproduction. The calculation was performed some years ago to leading logarithmic (log(Q 2 )) approximation, assuming the produced vector meson quarkonium system to be nonrelativistic [14,15] and improved in distinct aspects [55,56]. Assuming a non-relativistic wave function for the vector meson one have that the t = 0 differential cross section of photoproduction of heavy vector mesons in leading order approximation is given by [14,15] dσ γh→V h dt where xg h is the target gluon distribution and x = 4Q 2 /W 2 , with W the γh center of mass energy, Q 2 the characteristic hard scale of the processes and N = 1 at leading order (see below). Moreover, Γ ee is the leptonic decay width of the vector meson. In Refs. [55][56][57][58][59]  . It is important to emphasize that magnitude of these corrections is still a subject of discussions [60][61][62]. In particular, the value of the hard scaleQ is not fixed reliably at leading order. Such corrections have direct impact in the normalization of the cross section and in its energy dependence. In what follows we take into account of the contributions associated to real part of the scattering amplitude and to the skewness effect, which is equivalent to multiply the Eq. (4) by a factor (1 + β)R 2 g , where In order to estimate the total cross section we assume an exponential parametrization for the small |t| behaviour of the amplitude, which implies with b V being given by the following parametrization [33] b V (W ) = 4.9 + 0.24 ln(W/90 GeV), which is compatible with the HERA data. Moreover, we consider the values of mass and electronic decay widths of the vector mesons as given in Ref. [63]. Finally, in our calculations of the rapidity distribution we will assume that the absorptive corrections S 2 (W ± ) are given by the gap survival probability computed in [33,34] using the model proposed in Ref. [64].
In Fig. 1 we present the predictions for the gluon distribution obtained at leading order by the different groups [5][6][7][8][9] that perform the global parton analysis of the experimental data to extract the parton distributions. We show the results for two different values for the hard scale: Q = M J/Ψ /2 (left panel) and Q = M Υ /2 (right panel). These values are usually assumed as beingQ in the calculations of the exclusive J/Ψ and Υ photoproduction cross sections, respectively. We have that the distinct predictions differ significantly, mainly at small -x and low Q, which demonstrate that the global analysis do not reliably determine the gluon in this kinematical range. Basically, these distinct behaviours of the gluon distribution for low values of the hard scale are directly connected with the different assumptions for the initial conditions of the DGLAP equations considered by the distinct groups. With the increasing of the hard scale, the predictions becomes less dependent of these assumptions, and the distinct predictions becomes similar at very large Q. These results demonstrate the importance of probes of the gluon distribution at small -x and low Q. In Fig. 3 we present the resulting predictions for the energy dependence of the exclusive J/Ψ photoproduction cross section obtained from Eq. (4) assumingQ = M J/Ψ /2. The experimental data from Refs. [21][22][23][65][66][67][68] are presented for comparison. We obtain that the different models for the gluon distribution are not able to describe the normalization and/or the energy dependence of the data. The differences present in Fig. 1 are amplified in the exclusive J/Ψ photoproduction cross section due to its quadratic dependence on xg. In particular, the distinct xbehaviour predicted by the NNPDF and CT10 parametrizations for Q = M J/Ψ /2, with xg decreasing in the range 10 −4 ≤ x ≤ 10 −2 before to increase for x < 10 −4 , implies the anomalous behaviour in σ γp→J/Ψp observed in Fig. 3.
A possible interpretation of the results presented in Fig. 3 is that the different models fail to describe the data since our predictions for σ γp→J/Ψp have been obtained at leading order, which can be strongly affected by theoretical uncertainties associated to higher order corrections [60,61]. It is generally believed that these higher order corrections should be important, but a full calculation still remains a challenge. In order to analyse the possible impacts of the higher order corrections in our phenomenological analysis, in what follows we will assume that the general behaviour of the cross section, Eq. (4), will remain unchanged after the inclusion of the higher order corrections and that they  can be effectively incorporated by taking the values of the normalization and the hard scale from a fit to the γh data. In other words, we will assume that the main effect of the higher order corrections will be the modification of the normalization of the cross section and that they can be taken into account by an appropriate choice of the hard scale (For a similar approach see [32]). Basically we will assumeQ = ξM V /2 and will take ξ and N in Eq.
(4) as free parameters to be determined by fitting the experimental data for the σ γp→V p cross section. We consider the experimental data from Fermilab, HERA and LHC for the J/Ψ, Ψ(2S) and Υ production and performed the minimization of χ 2 /d.o.f. in order to determine ξ and N . These new parameters are then used to calculate the rapidity distribution for the vector meson photoproduction in pp collisions.
In Table I we present the fitted parameters for the different processes and a comparison between the predictions and the HERA data are shown in the left panel of the Figs. 4, 5 and 6. In general, the results gives reasonable values for the χ 2 /d.o.f., inside the 98% confidence level. We obtain smaller values of χ 2 /d.o.f. for heavier states, mainly due the small number of experimental data. In contrast, for the J/ψ data, a small χ 2 /d.o.f. is only obtained using CT10 and NNPDF parametrizations, which describe the data in the full energy range, with the other parametrizations being able to describe the data in a restrict range. Moreover, as can be verified from the analysis of the Figs. 4, 5 and 6, the predictions of the distinct parametrizations differ significantly in the kinematical range beyond the HERA data. It is important to emphasize that the results for ξ and N obtained using the CT10 and NNPDF parametrizations indicate that in order to describe the data for J/Ψ and Ψ(2S) production a large amount of higher order corrections are necessary. In contrast, these corrections are small for the Υ case.
In the right panel of the Figs. 4, 5 and 6 we present the corresponding predictions for the rapidity distributions for the vector meson photoproduction in pp collisions at √ s = 7 TeV. We compare our predictions with the recent data from the LHCb Collaboration [22][23][24]. As the rapidity distributions probe a large range of γp center of mass energies, the differences present in the predictions for σ γp→V p are amplified in dσ/dy, particularly for large rapidities in the case of the lighter mesons. For J/Ψ and Ψ(2S) production the models characterized by a strong increasing of the gluon distribution at small -x and small hard scales (Alekhin, MMHT and GJR) predict a double peak structure in the distributions. In contrast, the CT10 and NNPDF parametrizations predict a flat y distribution in a large rapidity range. These parametrizations describe quite well the LHCb data for the J/Ψ production and overestimate the Ψ(2S) data. It is important to emphasize that we have obtained these predictions using the model for S 2 (W ± ) proposed in [33,34]. Another aspect is the description of the Ψ(2S) wave function used in our calculations, which is currently a subject of intense debate. In the case of the Υ production, we obtain that the different models are not able to describe the current LHCb data, with the predictions at central rapidities being largely distinct. These results can indicate e.g. smaller values for S 2 (W ± ) than those proposed in [33]. Certainly, the future CMS data for y ≈ 0, which are currently under analysis, will be important to get more definitive conclusions.
A comment is in order. In our study we have considered leading order gluon distributions in order to be theoretically consistent with the fact that our expression for the total cross section have been derived at the leading logarithmic approximation. However, we have verified that our main conclusions are not modified if next -to -leading -order (NLO) gluon distributions are used as input in the calculations. Basically, the NLO corrections in the global analysis implies that the associated gluon distributions are not so steep at small -x and low values of the hard scale. As a consequence, a better description of the HERA data for the vector meson production is feasible. For example, in the case of the GJR parametrization, the value of xg for x = 10 −4 and µ = M J/Ψ /2 decreases by a factor 2. The resulting values for the χ 2 /d.o.f. for the description of the HERA data becomes 1.5 / 2.2 / 1.9 for the J/Ψ/Ψ(2S)/Υ production, respectively. However, these values still are larger than those obtained using the LO and NLO CT10 parametrizations. Moreover, as the difference between the GJR and CT10 gluon distributions increases at smaller values of x, probed at larger values of |y|, the resulting GJR (NLO) predictions for the rapidity distributions are not able to describe the LHC data. Similar results are obtained using the NLO gluon distributions derived in Refs. [5,8].
In the case of the CT10 and NNPDF parametrizations, we have that the impact of the NLO corrections is small, with  the resulting gluon distributions being similar to those obtained at leading order. As a consequence, the associated predictions for the rapidity distributions are similar to those presented in the Figs. 4, 5 and 6. These results indicate that the HERA and LHC data are better described when the gluon distribution present a slow increasing at smallx and low Q, in agreement with conclusion obtained at leading order. Finally, lets summarize our main conclusions. Exclusive vector meson photoproduction offers a unique opportunity to constrain the gluon density of the proton in the kinematical range of small -x and low Q 2 , which is the range in which the global analysis do not reliable determine the gluon distribution. In this paper we have performed a  phenomenological study of this process using the leading logarithmic formalism and different models for xg predicted by distinct groups that perform the global analysis of the experimental data in order to obtain the PDFs. We demonstrated that the combination of anyone of these models and the leading logarithmic formalism is not able to describe the experimental data for the J/Ψ photoproduction. Assuming that the main modifications in the formalism due to higher order corrections are the change in the normalization of the cross section and in the value of the hard scale probed in the process, we have fitted the current experimental data for the J/Ψ, Ψ(2S) and Υ photoproduction and used this new set parameters as input to calculate the rapidity distributions for the vector meson photoproduction in pp collisions. We have demonstrated that a better description of the γh data is obtained when the gluon distribution presents a slow increasing at small -x and low Q. Moreover, our results indicated that the LHCb data also are better described by these models. However, results for Ψ(2S) and Υ shown that the normalizations are not perfectly described, which can indicated the a more detailed analysis about the magnitude of the factor S 2 (W ± ) for the different final states is important in future. Future experimental data, particularly at central rapidities, will be important to constrain the different models. A final comment is in order. In this paper, in order to estimate the exclusive photoproduction of vector mesons, we have assumed the leading logarithmic formalism and different solutions of the linear DGLAP equation. Our goal was to try to describe the current experimental data for the exclusive vector photoproduction by improving this formalism and assuming that non linear effects in the QCD dynamics can be disregarded in this process. Our results indicate that it is not an easy task. Certainly, more definitive conclusions could be obtained using the full NLO expression for the cross section, which still is in progress. However, it is important to emphasize that the vector meson photoproduction at HERA and LHC is quite well described using the color dipole formalism, taking into account the non linear effects in the QCD dynamics and assuming S 2 (W ± ) = 1 (See e.g. Refs [41,42]). Such aspects demonstrated that future experimental results for the exclusive vector meson photoproduction in hadronic colliders are fundamental to understanding the QCD dynamics at high energies.