Energy loss effect of incoming gluons from J/ψ production in p-A collisions

The energy loss effect of incoming gluons from J/ψ production in p-A (or d-A) collisions is investigated by means of the E866, RHIC and LHC experimental data. The gluon mean energy loss per unit path length dE/dL = 2.18 ± 0.14 GeV/fm is extracted by fitting the E866 experimental data for J/ψ production cross section ratios RW(Fe)/Be(xF). The obtained result indicates that the incoming gluons lose more energy than the incident quarks. By comparing the theoretical results with E866, RHIC, and LHC experimental data, it is found that the nuclear suppression due to the incident gluon (quark) energy loss reduces (increases) with the increase of the kinematic variable xF (or y). The energy loss effect of incoming gluons plays an important role in the suppression of J/ψ production in a wide energy range from to , and the influence of incident quark energy loss can be ignored for high energies (such as at RHIC and LHC energy).


Introduction
In order to quantify the properties of the quark-gluon plasma (QGP) created in heavy-ion collisions, a solid understanding of the nuclear modification of particle spectra in cold nuclear matter is fundamentally important. J/ψ production in proton-nucleus collisions provides an ideal tool to test the microscopic dynamics of mediuminduced parton energy loss.
Drastic nuclear suppression effects are observed over a wide collision energy range for minimum bias p-A and d-A collisions, such as at the NA3 [1], E772 [2], E866 [3,4], NA50 [5], HEAR-B [6], LHC [7,8] and RHIC [9]experiments. However, it is striking that there is no consensus on the origin of J/ψ suppression in some kinematical conditions [10]. Some approaches attribute J/ψ suppression to an effective absorption cross section σ abs of the cc pair [11][12]; other models attribute J/ψ suppression to the increase of the cc pair invariant mass by multiple soft rescatterings through the nucleus, leading to a reduction of the overlap with the J/ψ wave function [13].
In the nucleus rest frame, a high-energy J/ψ is formed long after the nucleus, thus what actually propagates through the nucleus is the parent cc pair. Our previous works [14][15] support that the nuclear modification of the parton distribution functions and the incident proton energy loss owing to multiple scattering on the surrounding nucleon and gluon radiation are the main initial state effects which induce the J/ψ suppression, and the energy loss of the color octet cc is the dominant final state effect when the cc pair remains colored on its entire path through the medium. In Ref. [15], by using the EPS09 nuclear parton distributions [16] together with the energy loss of the proton beam in the initial state (the center-of-mass system energy loss per collision ∆ √ s = 0.18 GeV is determined from the nuclear Drell-Yan experimental data in the Glauber model [17]) and the linear quark energy loss in the final state, we extracted the charm quark mean energy loss per unit path length (dE/dL = 1.49±0.37GeV/fm with χ 2 /ndf = 0.91) by fitting the E866 experimental data [4] in the region 0.2 < x F < 0.65.
To further investigate the microscopic dynamics of medium-induced parton energy loss, the color charge of parton energy loss has received significant interest. This issue is of fundamental importance for accurately understanding the dynamics of modifying a hard probe and the dense QCD properties of what is probed. Previ-  [18][19][20][21][22] predicted that gluons lose more energy than quarks because of the stronger coupling to the medium.
In the J/ψ production for p-A (or d-A) collisions, the observed suppression induced by the incident parton energy loss effect can give a better way to discriminatingly identify the energy loss of incoming gluons and quarks. Following our previous work, in the present study we investigate the incoming gluon energy loss effect by means of the E866 [4], RHIC [9] and LHC [7,8] experimental data, and hope that our research can provide a useful reference for deep understanding of the microscopic dynamics of medium-induced parton energy loss.
The remainder of this paper is organized as follows. In Section 2, the theoretical framework of our study is introduced. Section 3 is devoted to the results and discussion. Finally, a summary is presented.

Formalism for J/ψ production differential cross sections
In the the color evaporation model (CEM) [23], for J/ψ production in p-A collisions, quarkonium production is treated identically to open heavy-quark production except that the invariant mass of the heavy quark pair is restricted to be less than twice the mass of the lightest meson that can be formed with one heavy constituent quark. For charmonium the upper limit on the cc pair mass is then 2m D . The hadroproduction of a heavy quark at leading order (LO) in perturbative QCD is the sum of contributions from qq annihilation and gg fusion. The charmonium production cross section dσ p−p /dx F is a convolution of the qq and gg partonic cross sections with the parton distribution functions f i in the incident proton and f i in the target proton, and is expressed as [24]: Here, in the rest frame of the target nucleus, x 1(2) is the projectile proton (target) parton momentum fractions, √ s is the center of mass energy of the hadronic collision, m 2 = x 1 x 2 s, m c = 1.2 GeV and m D = 1.87 GeV are respectively the charm quark and D meson mass, σ gg (σ qq ) is the LO cc partonic production cross section from the gluon fusion (quark-antiquark annihilation), and ρ J/ψ is the fraction of the cc pair which produces the J/ψ state. In J/ψ production from p-A (or d-A) collisions, owing to multiple scattering on the surrounding nucleon and gluon radiation while the incident parton propagates through the nucleus, the incoming gluon (quark) can lose its energy ∆E g (∆E q ). The energy loss of the incoming gluon (quark) results in an average change in its momentum fraction prior to the collision, According to the parametrization for parton energy loss [25][26], the mean energy loss of an incoming gluon (quark) can be expressed as: Here, , and α, β are the parameters that can be extracted from the experimental data by adopting the χ 2 analysis method. When the J/ψ hadronization occurs outside the nucleus, nuclear absorption should play little or no role and the energy loss of the color octet cc is the dominant final state effect. In view of the shift in x F due to the energy loss of color octet cc (∆E cc ), the momentum fraction of the incident gluon (quark) is actually: [15], τ = m 2 /s. The J/ψ differential production cross section in p-A collisions dσ p−A /dx F is written as: dm 2m Here, in consideration of the shift in x F due to the energy loss of color octet cc, the target parton momentum fraction is actually ]. Further, considering the energy loss of the incident gluon, incoming quark and the color octet cc, the leading order for J/ψ production cross section as a function of y should be written as: Here, with and

Results and discussion
In order to determine the value of the incoming gluon energy loss parameter α, we give a phenomenological analysis at leading order for the J/ψ production cross section ratios R W(Fe)/Be (x F ): for the E866 experimental data (49 points) by using the EPS09 nuclear parton distributions [16] together with the energy loss parameter of the incident quark (β = 1.21 ± 0.09 GeV/fm) determined from the nuclear Drell-Yan experimental data [26] and the color octet cc energy loss (α = 2.97 GeV/fm) determined in our previous work [15]. By minimizing χ 2 with the CERN subroutine MINUIT [28] the value of parameter α is extracted: α = 2.18 ± 0.14 GeV/fm. One standard deviation of the optimum parameter corresponds to an increase of χ 2 by 1 unit from its minimum χ 2 min . The result indicates that the incoming gluons lose more energy than the incident quarks in J/ψ production from p-A collisions, which is in accord with the prediction that gluons lose more energy than quarks because of the stronger coupling to the medium [18][19][20][21][22]. In addition, due to the effects of the modification of the gluon parton distribution functions on the nucleus leading to an additional J/ψ suppression in p-A collisions, the EPS09 uncertainties can be the main source of uncertainty associated with our results.
To identify the energy loss effect of the incoming gluon and incident quark on the J/ψ suppression, the theoretical results are compared with E866 experimental data [4] at √ s = 38.7 GeV in Figs. 1 and 2, RHIC experimental data [9] at √ s = 200 GeV in Fig. 3, and LHC experimental data [7,8] at √ s = 5.0 TeV in Fig.   4, respectively. The dotted, dashed and solid lines correspond to the results given without initial state energy loss, by considering the incident quark energy loss effect, and the energy loss of the incident quark together with incoming gluon energy loss.  As can be seen in Figs. 1 and 2, the nuclear suppression due to the incident quark energy loss is negligible in the region x F < 0.3, increases gradually for x F < 0.8, and becomes steeper for x F > 0.8. The suppression from the energy loss effect of the incoming gluon, on the other hand, is much steeper in the region x F < 0.3, reduces gradually for x F < 0.8, and becomes negligible for x F > 0.8. It is clear that the incident gluon energy loss plays an important role in the suppression of the J/ψ production cross section ratios R W(Fe)/Be (x F ) in the small x F region (especially for x F < 0.3), and the energy loss effect of incoming quarks is obvious in the large x F region (especially for x F > 0.8). We can see that the experimental data on J/ψ production at E866 energy ( √ s = 38.7 GeV) can give a good test for the identity of the incident parton which loses its energy in the nuclear medium.
In Figs. 3 and 4, the theoretical results for J/ψ production cross section ratios R Au(Pb)/p as a function of y are compared with RHIC [9] and LHC [7,8] experimental data, respectively. From Fig. 3 we can see that the dotted line and the dashed line appear to overlap, which indicates that the energy loss effect due to the incoming quark plays no role in the J/ψ production at RHIC energy. In contrast, the nuclear suppression due to the incident gluon energy loss is obvious, especially in the range y < −1.5, reduces gradually with the increase of y, and becomes negligible for y > 2.0. From Fig. 4, the incident quark energy loss effect has little impact on the J/ψ production cross section ratio R Pb/p (y) at LHC energy, and the energy loss due to the incoming gluon plays an important role in the nuclear suppression, especially in the range y < −3.5, reduces gradually with the increase −2.0 The solid triangles are the RHIC experimental data [9]. The solid triangles and filled circles correspond to the experimental data from the ALICE Collaboration [7] and LHCb Collaboration [8] at the LHC, respectively.
of y, and becomes negligible for y > −1.5. In the present work, it is found that the energy loss of the incoming gluon plays an important role in the suppression of J/ψ production in a wide energy range from √ s = 38.7 GeV to √ s = 5.0 TeV, and the influence of incident quark energy loss can be ignored for high energies, such as RHIC energy and LHC energy.

Summary
Following our previous work [14,15], we have studied the energy loss effect of incoming gluons from J/ψ production in p-A (or d-A) collisions. By means of the EPS09 nuclear parton distributions [16] together with the energy loss of the incident quark (dE/dL = 1.21±0.09 GeV/fm determined in our work [26]) and color octet cc (dE/dL = 2.97 ± 0.74 GeV/fm determined in our study [15]), we have given a phenomenological analysis at leading order for the J/ψ production cross section ratios R W(Fe)/Be (x F ) for the E866 experimental data (49 points) and extracted the gluon mean energy loss per unit path length dE/dL = 2.18 ± 0.14 GeV/fm by minimizing χ 2 with the CERN subroutine MINUIT [28]. This result indicates that the incoming gluons lose more energy than the incident quarks, which supports the prediction that gluons lose more energy than quarks because of the stronger coupling to the medium [18][19][20][21][22]. In addition, the EPS09 uncertainties can be the main source of uncertainty associated with our results, owing to the effects of the modification of the gluon parton distribution functions on the nucleus leading to an additional J/ψ suppression in p-A collisions. To identify the energy loss effect of the incoming gluon and quark on the J/ψ suppression, the theoretical results were compared with E866 [4], RHIC [9], and LHC [7,8] experimental data. We find that the nuclear suppression due to the incident gluon (quark) energy loss reduces (increases) with the increase of the kinematic variable x F (or y). The energy loss of the incoming gluon plays an important role in the suppression of J/ψ production in a wide energy range from √ s = 38.7 GeV to √ s = 5.0 TeV, and the influence of incident quark energy loss can be ignored for high energy, such as at RHIC energy and LHC energy.