# Forward production at high energy: centrality dependence and mean transverse momentum

###### Abstract

Forward rapidity meson production in proton-nucleus collisions can be an important constraint of descriptions of the small- nuclear wavefunction. In an earlier work we studied this process using a dipole cross section satisfying the Balitsky-Kovchegov equation, fit to HERA inclusive data and consistently extrapolated to the nuclear case using a standard Woods-Saxon distribution. In this paper we present further calculations of these cross sections, studying the mean transverse momentum of the meson and the dependence on collision centrality. We also extend the calculation to backward rapidities using nuclear parton distribution functions. We show that the parametrization is overall rather consistent with the available experimental data. However, there is a tendency towards a too strong centrality dependence. This can be traced back to the rather small transverse area occupied by small- gluons in the nucleon that is seen in the HERA data, compared to the total inelastic nucleon-nucleon cross section.

###### pacs:

13.85.Ni, 14.40.Pq, 24.85.+p, 25.75.Cj## I Introduction

The production of mesons at forward rapidity in high energy proton-proton and proton-nucleus collisions can provide valuable information on gluon saturation. Indeed, the production of particles at forward rapidity probes the target at very small , where saturation effects should be enhanced. In particular, the charm quark mass being of the same order of magnitude as the saturation scale, production should be sensitive to these dynamics. The charm quark mass is also large enough to provide a hard scale, making a perturbative study of this process possible. In addition, production has been the subject of many experimental studies at the LHC, both in proton-proton Khachatryan:2010yr (); Aaij:2011jh (); Aad:2011sp (); Chatrchyan:2011kc (); Abelev:2014qha (); Aaij:2015rla (); Adam:2015rta () and in proton-nucleus Abelev:2013yxa (); Aaij:2013zxa (); Adam:2015iga (); Aad:2015ddl (); Adam:2015jsa () collisions. This provides a lot of data to confront with nuclear effects predicted by various theoretical models, both in the Color Glass Condensate (CGC) framework Fujii:2013gxa (); Ducloue:2015gfa (); Ma:2015sia (); Watanabe:2015yca (); Fujii:2015lld () and as constraints for nuclear parton distribution functions Albacete:2013ei () and energy loss in cold nuclear matter Vogt:2010aa (); Arleo:2012rs (); Arleo:2014oha (); Vogt:2015uba ().

In a recent work Ducloue:2015gfa () we studied, in the CGC framework, the production of forward mesons in proton-proton and proton-nucleus collisions at the LHC. We showed that, when using the Glauber approach to generalize the dipole cross section to nuclei, the nuclear suppression for minimum bias events is smaller than in previous CGC calculations such as Fujii:2013gxa (), and much closer to experimental data. In this paper we will study, in the same framework, other observables of interest in this process, such as production at backward rapidity, the centrality dependence in the optical Glauber model and the mean transverse momentum of the produced ’s.

## Ii Formalism

Let us briefly recall the main steps of the calculation. For more details we refer the reader to Ref. Ducloue:2015gfa (). We use the color evaporation model (CEM) to relate the cross section for production to the pair production cross section. In this model a fixed fraction of the pairs produced below the -meson threshold are assumed to hadronize into mesons:

(1) |

where , and are the transverse momentum, the rapidity and the invariant mass of the pair respectively, GeV is the D meson mass and is the charm quark mass that we will vary between 1.2 and 1.5 GeV. Note that in this work we will focus on ratios where cancels so we do not need to fix it to any specific value here.

The study of gluon and quark pair production in the dilute-dense limit of the CGC formalism was started some time ago Blaizot:2004wu (); Blaizot:2004wv () (see also Kharzeev:2012py ()) and used in several calculations such as Fujii:2006ab (); Fujii:2005rm (); Fujii:2013gxa (); Fujii:2013yja (). The physical picture is the following: an incoming gluon from the projectile can split into a quark-antiquark pair either before or after the interaction with the target. The partons propagating trough the target are assumed to interact eikonally with it, picking up a Wilson line factor in either the adjoint (for gluons) or the fundamental (for quarks) representation. Since we study production at forward rapidity, where the projectile is probed at large , we will use the collinear approximation in which the incoming gluon is assumed to have zero transverse momentum. In this approach the cross section for production reads, in the large- limit Fujii:2013gxa (),

(2) |

where and denote the transverse momenta of the quarks, and their rapidities, , and is the dimension of the adjoint representation of SU(). The longitudinal momentum fractions and probed in the projectile and the target respectively are given by

(3) |

The explicit expression for the “hard matrix element” is given in Ref. Fujii:2013gxa (). In Eq. (2), is the gluon density in the probe and is described, in the collinear approximation that we use here, in terms of usual parton distribution functions (PDFs). In the following we use, unless otherwise stated, the MSTW 2008 Martin:2009iq () parametrization for . For consistency we use the leading other (LO) parton distribution since also the rest of the calculation is made at this order. To estimate the uncertainty associated with the choice of the factorization scale , we will vary it between and .

The function describes the propagation of a pair in the color field of the target and reads

(4) |

where is the impact parameter. In this expression, the function contains all the information about the target. It is the fundamental representation dipole correlator:

(5) |

with

(6) |

where is a fundamental representation Wilson line in the color field of the target. The dipole correlator is obtained by solving numerically the running coupling Balitsky-Kovchegov (rcBK) equation Balitsky:1995ub (); Kovchegov:1999ua (); Balitsky:2006wa (). For the initial condition we use the ’MV’ parametrization introduced in Ref. Lappi:2013zma (), which reads, in the case of a proton target,

(7) |

where . The running coupling is taken as:

(8) |

where parametrizes the uncertainty related to the scale of the strong coupling in the transverse coordinate space.

The free parameters in these expressions are obtained by fitting the combined inclusive HERA DIS cross section data Aaron:2009aa () for GeV and . Their best fit values (with ) are GeV, , and mb. In the case of a proton target, the dipole amplitude does not have an explicit impact parameter dependence and we thus make the replacement

(9) |

in Eq. (4), corresponding to the effective proton transverse area measured in DIS experiments. The fit only includes light quarks. In particular this leaves the charm quark mass that would be consistent with the DIS data in this model still uncertain, which is why we vary it in a rather large range for the uncertainty estimate in this work.

To generalize this proton dipole correlator to the case of a nuclear target we use, as in Lappi:2013zma (), the optical Glauber model. In this model the gluons at the initial rapidity are localized in the individual nucleons of the nucleus. The nucleons are then taken to be distributed randomly and independently in the transverse plane according to the standard Woods-Saxon nuclear density profile. An analytical average over the positions of the nucleons leads to the following initial condition for the rcBK evolution of a nuclear target:

(10) |

Here is the standard Woods-Saxon transverse thickness function of the nucleus:

(11) |

with and , and is defined such that is normalized to unity. All the other parameters in the initial condition (10), which is evolved using the rcBK equation for each , are the same as in the proton case.

In this model the dipole amplitude depends on the impact parameter and we need to integrate explicitly over it. The impact parameter dependence, which carries over to the centrality dependence, thus appears naturally in this model. Our practical procedure for carrying out this comparison will be discussed in more detail in Sec. IV.

## Iii Backward rapidity

In our previous work Ducloue:2015gfa () we only considered production in proton-proton and proton-nucleus collisions at forward rapidity. In these kinematics the process can be seen as the collision of a dilute proton probed at large , which can be described using well known parton distribution functions (PDFs), and a dense target described in terms of classical color fields. The nuclear modification of production was also measured at backward rapidity by ALICE Abelev:2013yxa (); Adam:2015iga () and LHCb Aaij:2013zxa (). In this case the produced is moving in the direction of the incoming nucleus and the physical picture is the same as at forward rapidity, but with the roles of the projectile and the target interchanged, i.e. a dilute proton or nucleus interacts with a dense proton target. The latter is described in the same way as in proton-proton collisions at forward rapidity, while the projectile is described either by a PDF in the case of a proton or by a nuclear PDF (nPDF) in the case of a nucleus. Therefore the calculation is very similar to the case of proton-proton collisions in our previous study. Note, however, that in general nuclear PDFs are less tightly constrained by experimental data than usual proton PDFs. In the following we will use the leading order EPS09 nPDF parametrization Eskola:2009uj () which provides additional error sets to estimate this uncertainty. The LO EPS09 analysis uses the CTEQ6L1 Pumplin:2002vw () proton PDFs as a reference, therefore for consistency we use same parametrization of proton PDFs when computing the nuclear modification factor at backward rapidity. Nevertheless, while it could be sizeable for the cross section, the difference compared to the MSTW 2008 parametrization is very small for the nuclear modification factor.

In Fig. 1 we show the nuclear modification factor , defined as

(12) |

as a function of and obtained in this way at negative rapidity compared with data from ALICE Abelev:2013yxa (); Adam:2015iga () and LHCb Aaij:2013zxa () experiments. The uncertainty in our calculation is significantly larger than at forward rapidity Ducloue:2015gfa (). This is due to the fact that we include in our uncertainty band, in addition to the variation of the charm quark mass and of the factorization scale, the nuclear PDF uncertainty obtained following the procedure described in Ref. Eskola:2009uj (). In particular, our lower bound for the factorization scale is which can reach values of less than 2 GeV at small transverse momentum. The nuclear PDFs are not well constrained at such small scales at the moment. Nevertheless the general agreement with data is quite good taking into account the rather large theoretical and experimental uncertainties. In this case deviations of from unity are entirely due to the nuclear PDFs.

Now that we have computed the nuclear modification factor both at forward Ducloue:2015gfa () and backward rapidities, we have access to the forward to backward ratio , defined as

(13) |

This ratio can be interesting to study because there may be an additional cancellation of some uncertainties common to the numerator and the denominator. In particular, when determining the nuclear modification factor, experimental studies such as Abelev:2013yxa (); Aaij:2013zxa () have to use an interpolation for the reference proton-proton cross section since there is no data at TeV. This interpolation is not needed to study , but the final statistical uncertainty may be larger if the coverage in rapidity by the detector is not symmetric with respect to 0. Concerning our calculation, we have seen that at negative rapidity the uncertainty on due to nuclear PDFs is rather large. This error will remain in since the computation of at positive rapidity does not involve nPDFs. Indeed, we see from Fig. 2, where we show the forward to backward ratio as a function of and , that the uncertainty on this quantity is still quite large. Nevertheless, within this error band the agreement with data is reasonable, although the variation at low seems to be steeper than in the data. In Fig. 2 (R) we only show our results for as a function of integrated over the same range as ALICE data Abelev:2013yxa (), which is slightly smaller than for LHCb data Aaij:2013zxa (), but our results for would be very similar and ALICE and LHCb data are compatible with each other.

Centrality class | [fm] | ||
---|---|---|---|

2–10% | 14.7 | 4.14 | |

10–20% | 13.6 | Â | 4.44 |

20–40% | 11.4 | 4.94 | |

40–60% | 7.7 | 5.64 | |

60–80% | 3.7 | 6.29 | |

80–100% | 1.5 | 6.91 |

## Iv Centrality dependence

We have seen that the optical Glauber model contains an explicit impact parameter dependence which can be related to centrality determinations at experiments. In this section we will compare this centrality dependence with its recent measurement at the LHC by the ALICE collaboration in the range Adam:2015jsa ().

### iv.1 Centrality in the optical Glauber model

We still need to relate the explicit impact parameter dependence in our model to the definition of centrality used by experiments. The experimental data are usually presented in terms of centrality classes. In the optical Glauber model these classes would be defined by calculating the impact parameter range corresponding to the centrality class using the relation

(14) |

Here is the total inelastic proton-nucleus cross section, given by

(15) |

and the scattering probability at impact parameter is

(16) |

where is the total inelastic nucleon-nucleon cross section. The particle yield in each centrality class is then given by

(17) |

where the values of and are obtained from Eq. (14) and corresponds to the expression of the cross section before integrating over .

In practice, however, this straightforward procedure cannot be used for comparing our calculation with the experimental centrality classes. This can be seen most clearly by calculating the average number of binary nucleon-nucleon collisions for each centrality class. In the optical Glauber model this is computed using the relation

(18) |

with

(19) |

Table 1 shows for the optical Glauber centrality classes compared to the values given by ALICE Adam:2015jsa () for the experimental ones. For central collisions the average number of binary collisions estimated by ALICE is smaller than in the optical Glauber model, while the opposite is true for peripheral collisions.

The impact parameter is not directly observable, so the experimental centrality selection has to use some other observable as a proxy for it. In the case of the ALICE analysis Adam:2015jsa () the observable used is the energy of the lead-going side Zero-Degree Calorimeter (ZDC). Due to the large fluctuations in the signal for a fixed impact parameter, the values of vary less strongly with centrality in the experimental classes than in the optical Glauber ones. Thus, for example, a relatively peripheral smaller- event can end up in the most central class in the case of a fluctuation in the ZDC signal.

### iv.2 Fixed impact parameter approximation

Developing a detailed Monte Carlo Glauber model that would enable us to exactly match the experimental centrality classes would be beyond the scope of this work. Therefore we start from the assumption that the ALICE Glauber model is correct and produces a reliable estimate for the in each experimental class. We then assume that the mapping between this and the impact parameter is accurately enough described by our optical Glauber model, and use the relation to determine a mean impact parameter corresponding to the experimental centrality class. The values of resulting from this procedure are shown in Table 1. We then calculate the nuclear modification factor for the centrality class using this fixed value of . We note that following this procedure for the 80-100% centrality class would lead to a value of for which our calculation would not be applicable because the saturation scale of the nucleus falls below the one of the proton. Therefore in the following we will only show results for the five most central classes considered by ALICE.

In Fig. 3 we show a comparison of our calculation with ALICE data for the nuclear modification factor in different centrality classes, , defined as

(20) |

as a function of . We observe that the description of experimental data is generally satisfactory in the first four bins. For the fifth bin the value of obtained in the optical Glauber model is almost constant and very close to one, while the data still shows a significant variation with .

The discrepancy with experimental data for peripheral collisions comes from the fact that in our model the saturation scale of the lead nucleus falls below the one of the proton for a value of of the order of 6.3 fm, as shown in Fig. 4 (L). We see from Fig. 4 (R), where we show the number of binary collisions as a function of , that this corresponds to . This is the point where, by definition, reaches 1 in our calculation and beyond which the validity of the framework we have used is questionable. On the other hand, ALICE data shows that is still significantly smaller than 1 down to (see Fig. 6).

The strong centrality dependence is caused by the value of extracted from DIS fits being much smaller than the total inelastic nucleon-nucleon cross section. The HERA data, both the inclusive cross section fitted in Ref. Lappi:2013zma () and data on exclusive vector meson production (see e.g. Chekanov:2004mw ()), lead to a picture where the small- gluons that can participate in a hard process in a proton are concentrated in a rather small area in the transverse plane (see also the discussion in Ref. Frankfurt:2010ea ()). This small- gluon “hot spot” is then surrounded by a larger “cloud” that only participates in soft interactions, contibuting to the total nucleon-nucleon inelastic cross section. Our model takes this picture to the extreme, by assuming that at the initial rapidity the gluons contributing to production are concetrated in the area inside the target nucleons. Thus, for peripheral collisions, the probe proton can overlap with the soft cloud of target nucleons while still seeing on average only one small- gluon hot spot in the target, thus leaving a hard process like production approximately unmodified.

This centrality dependence could probably be mildened by using a larger value for of the order of (as is effectively done in Fujii:2013gxa ()). This would, however, lose the consistency of our description of the nucleon from HERA to the LHC. Also, as discussed in Refs. Lappi:2013zma (); Ducloue:2015gfa (), varying these parameters in an uncontrolled way could very easily, depending on how exactly it is done, lead to an excessive suppression for minimum bias collisions, or to an that is very far from unity even at high transverse momentum.

Note that the saturation scale (or gluon density) at the edge of the nucleus being smaller than that of the proton is an artefact of averaging over the transverse locations of the dense but small gluon hot spots of the nucleons in the target nucleus. Therefore we do not use the optical Glauber parametrization in this region, but explicitly set to unity. Nevertheless, even an explicit Monte Carlo Glauber procedure with the same parameters would not change the ordering that leads to the absence of nuclear effects for peripheral collisions with .

### iv.3 Explicit integration over the impact parameter

The results we have shown in the previous section were obtained with a fixed impact parameter chosen so that the number of binary collisions in the optical Glauber model is equal to the average number of binary collisions estimated by ALICE in each centrality class. However, the nuclear modification ratio in a given centrality bin receives contributions from a distribution of different impact parameters. One could therefore argue that having a profile in the impact parameter space and integrating over could lead to different results. To quantify this effect we will here use two different kinds of distributions to obtain distributions in the impact parameter space. The first one is provided by ALICE alicecent () and is obtained from the Slow Nucleon Model (SNM) Adam:2014qja (). Since in this model the average number of binary collisions is not the same as the one obtained in the hybrid method used in Ref. Adam:2015jsa (), we shift the distributions so that matches the one in the third column of Table 1. It should be noted that, contrary to the hybrid method, this method is biased Adam:2014qja (). In addition, this is only one possible way of extracting distributions at experiments. Other methods could yield significantly different distributions. To try to quantify the dependence of our results on the particular shape of the distributions, we will also use, for the 60-80% centrality bin which is the most sensitive to fluctuations, a simple linearly decreasing distribution. The two parameters of this distribution, its height at the origin () and the value at which it vanishes (), are determined by imposing that it is normalized to unity and that : , .

In Fig. 5 we show the values obtained for the nuclear modification factor as a function of in the 60-80% centrality bin, both when using a fixed impact parameter and when integrating explicitly over using the two distributions described previously (SNM and linear). The explicit integration over leads to a smaller at small transverse momentum. The effect is more pronounced with the linear distribution. Similar results are obtained when looking at integrated over as a function of , as show in Fig. 6. In particular, we see that the value of in the 60-80% centrality bin obtained with the linear distribution is significantly closer to the ALICE data point. On the other hand, the value obtained with the SNM distribution is very close to the fixed impact parameter result.

In conclusion, it is not possible for now to directly compare our impact parameter dependent results with the centrality dependent measurement performed by ALICE. Indeed, for this one would need to have access to an unbiased determination of the distributions in each centrality bin, which does not exist at the moment. Here we have tried to estimate the importance of this effect by using distributions obtained in two models. The fact that these two models lead to significantly different results for peripheral collisions while the central bins are much less sensitive to fluctuations means that the variation of as a function of centrality is too model dependent to have a reliable comparison with experimental data.

## V Mean transverse momentum

In Ref. Ducloue:2015gfa () we found that the uncertainty on cross sections both in proton-proton and proton-nucleus collisions was rather large. This uncertainty mostly affects the normalization and therefore quantities such as the nuclear modification factor show a smaller uncertainty. Besides the nuclear modification factor, another observable which is not sensitive to the absolute normalization of the cross section is the mean transverse momentum of the produced meson. In Fig. 7 (L) we show this quantity as a function of the rapidity in proton-proton collisions at a center of mass energy of 7 TeV and compare with LHCb data Aaij:2011jh (). We observe that our calculation is compatible with the data but it is still affected by a relatively large uncertainty. In the collinear approximation on the proton side that we are using here, the mean transverse momentum increases slightly with rapidity, a trend not seen in the data. One must, however, keep in mind that towards central rapidities the intrinsic transverse momentum from also the proton should increase, leading to the opposite behavior. A matching between the collinear and -factorized approximations required to fully quantify this effect is beyond the scope of this paper. In Fig. 7 (R) we show the same quantity in proton-lead collisions at a center of mass energy TeV.

The ALICE collaboration has also presented results for as a function of . On Fig. 8 (L) we see that our calculation agrees with this measurement within the rather large uncertainty band (except for most peripheral collisions, where our calculation is not applicable as explained previously). When one considers the difference between in proton-lead and in proton-proton collisions, as shown on Fig. 8 (R), the uncertainty on our calculation shrinks and shows a too strong variation as a function of , both when using a fixed impact parameter and when integrating explicitly over using the distributions obtained in the Slow Nucleon Model. However, as in section IV.3, the results obtained when integrating over depend strongly on the exact shape of the distributions used. In particular, one can see that using a linear distribution leads to a better agreement with data for peripheral collisions.

## Vi Dependence on the center of mass energy

### vi.1 Proton-proton collisions

In this work we use the simple color evaporation model to describe the hadronization of pairs into mesons. The normalization of cross sections then depends on a non perturbative constant , see (1). The uncertainty associated with this parameter can be eliminated by studying the ratio of cross sections at different center of mass energies. In addition, from the experimental point of view, systematic uncertainties can cancel to some extent in this ratio. Such a measurement has been made possible at the LHC for proton-proton collisions thanks to the recent increase of from 8 to 13 TeV. In particular the ratio was studied as a function of and by the LHCb collaboration Aaij:2015rla (). In Fig. 9 we compare these data with the results that we obtain for this ratio in our model. The resulting uncertainty is rather small and the agreement with data is quite good, in particular at large rapidity and relatively low transverse momentum which is the kinematical domain where our calculation is expected to be the most reliable.

### vi.2 Proton-nucleus collisions

Thanks to its recent upgrade, the LHC may also perform proton-lead collisions at a higher center of mass energy in the future. Here we study how our results would be affected by a change of from 5 to 8 TeV. In Fig. 10 we show the nuclear modification factor at forward rapidity as a function of and at these two energies, as well as existing LHC data at TeV. The values at TeV shown here differ slightly from the ones shown in Figs. 8 and 10 of Ref. Ducloue:2015gfa () because we corrected a numerical problem which was causing the region of large impact parameters (where the saturation scale of the lead nucleus falls below the one of the proton, see Fig. 4 (L)) to be neglected. Here we impose in this region, as in Ref. Lappi:2013zma (). As one could expect, the higher center of mass energy leads to a stronger suppression due to the higher densities reached in the target. However, the effect is quite small, in particular compared to the size of the uncertainties. For this reason we only show, for TeV, our results for the “central” values of the parameters ( GeV and ). For TeV we show both the central value and the uncertainty band corresponding to the variation of and . Therefore, while measuring forward production in 8 TeV proton-lead collisions could help reduce experimental uncertainties by getting rid of the interpolation needed for the proton-proton reference, we do not expect significantly stronger nuclear effects at this energy. This is not surprising since we use the same dipole cross sections as in Ref. Lappi:2013zma (), where a weak energy dependence of the nuclear modification factor was found in single inclusive particle production.

## Vii Conclusions

In this paper we have extended our study of forward production in proton-nucleus collisions in the Color Glass Condensate framework to new kinematics and observables. In particular we have studied the nuclear suppression at negative rapidities by describing the nucleus probed at large in terms of nuclear parton distribution functions. We achieved a quite good description of experimental measurements of this quantity, even if the uncertainty is larger at backward than at forward rapidities because the nuclear PDFs are not yet very strongly constrained by data. This allowed us to compute the forward to backward ratio, again with a good agreement with data within the rather large uncertainties. We have also studied the centrality dependence of our calculation. While using the optical Glauber model to extend the description from a proton target to a nucleus leads to a better agreement with experimental data for minimum bias observables than previous calculations in the same framework, it is difficult to compare directly the resulting centrality dependence to experimental data. Indeed, using a fixed impact parameter obtained in the optical Glauber model from the average number of binary collisions estimated by ALICE leads to a too strong centrality dependence. On the other hand, to integrate explicitly over the impact parameter one has to use distributions based on various assumptions and our calculation is very sensitive to the exact shape of these distributions. Finally we have studied how the nuclear modification factor would be affected by an increase of the center of mass energy achievable at the LHC. As expected nuclear effects are stronger but the change is too small to be significant given the size of theoretical uncertainties.

## Acknowledgments

T. L. and B. D. are supported by the Academy of Finland, projects 267321 and 273464. H. M. is supported under DOE Contract No. DE-SC0012704. This research used computing resources of CSC – IT Center for Science in Espoo, Finland. We would like to thank C. Hadjidakis and I. Lakomov for discussions on the ALICE data.

## References

- (1) CMS collaboration, V. Khachatryan et. al., Prompt and non-prompt production in collisions at TeV, Eur. Phys. J. C71 (2011) 1575 [arXiv:1011.4193 [hep-ex]].
- (2) LHCb collaboration, R. Aaij et. al., Measurement of production in collisions at , Eur. Phys. J. C71 (2011) 1645 [arXiv:1103.0423 [hep-ex]].
- (3) ATLAS collaboration, G. Aad et. al., Measurement of the differential cross-sections of inclusive, prompt and non-prompt production in proton-proton collisions at TeV, Nucl. Phys. B850 (2011) 387 [arXiv:1104.3038 [hep-ex]].
- (4) CMS collaboration, S. Chatrchyan et. al., and (2S) production in collisions at TeV, JHEP 02 (2012) 011 [arXiv:1111.1557 [hep-ex]].
- (5) ALICE collaboration, B. Abelev et. al., Measurement of quarkonium production at forward rapidity in collisions at TeV, Eur. Phys. J. C74 (2014) 2974 [arXiv:1403.3648 [nucl-ex]].
- (6) LHCb collaboration, R. Aaij et. al., Measurement of forward production cross-sections in collisions at TeV, JHEP 10 (2015) 172 [arXiv:1509.00771 [hep-ex]].
- (7) ALICE collaboration, J. Adam et. al., Inclusive quarkonium production at forward rapidity in pp collisions at TeV, Eur. Phys. J. C76 (2016) 184 [arXiv:1509.08258 [hep-ex]].
- (8) ALICE collaboration, B. Abelev et. al., production and nuclear effects in p-Pb collisions at = 5.02 TeV, JHEP 02 (2014) 073 [arXiv:1308.6726 [nucl-ex]].
- (9) LHCb collaboration, R. Aaij et. al., Study of production and cold nuclear matter effects in collisions at TeV, JHEP 02 (2014) 072 [arXiv:1308.6729 [nucl-ex]].
- (10) ALICE collaboration, J. Adam et. al., Rapidity and transverse-momentum dependence of the inclusive J/ nuclear modification factor in p-Pb collisions at 5.02 TeV, JHEP 06 (2015) 055 [arXiv:1503.07179 [nucl-ex]].
- (11) ATLAS collaboration, G. Aad et. al., Measurement of differential production cross sections and forward-backward ratios in pÂ +Â Pb collisions with the ATLAS detector, Phys. Rev. C92 (2015) 034904 [arXiv:1505.08141 [hep-ex]].
- (12) ALICE collaboration, J. Adam et. al., Centrality dependence of inclusive production in p-Pb collisions at TeV, JHEP 11 (2015) 127 [arXiv:1506.08808 [nucl-ex]].
- (13) H. Fujii and K. Watanabe, Heavy quark pair production in high energy pA collisions: Quarkonium, Nucl. Phys. A915 (2013) 1 [arXiv:1304.2221 [hep-ph]].
- (14) B. Ducloué, T. Lappi and H. Mäntysaari, Forward production in proton-nucleus collisions at high energy, Phys. Rev. D91 (2015) 114005 [arXiv:1503.02789 [hep-ph]].
- (15) Y.-Q. Ma, R. Venugopalan and H.-F. Zhang, production and suppression in high energy proton-nucleus collisions, Phys. Rev. D92 (2015) 071901 [arXiv:1503.07772 [hep-ph]].
- (16) K. Watanabe and B.-W. Xiao, Forward Heavy Quarkonium Productions at the LHC, Phys. Rev. D92 (2015) 111502 [arXiv:1507.06564 [hep-ph]].
- (17) H. Fujii and K. Watanabe, Leptons from heavy-quark semileptonic decay in pA collisions within the CGC framework, Nucl. Phys. A951 (2016) 45 [arXiv:1511.07698 [hep-ph]].
- (18) J. L. Albacete et. al., Predictions for Pb Collisions at TeV, Int. J. Mod. Phys. E22 (2013) 1330007 [arXiv:1301.3395 [hep-ph]].
- (19) R. Vogt, Cold Nuclear Matter Effects on and Production at the LHC, Phys. Rev. C81 (2010) 044903 [arXiv:1003.3497 [hep-ph]].
- (20) F. Arleo and S. Peigné, Heavy-quarkonium suppression in p-A collisions from parton energy loss in cold QCD matter, JHEP 03 (2013) 122 [arXiv:1212.0434 [hep-ph]].
- (21) F. Arleo and S. Peigné, Quarkonium suppression in heavy-ion collisions from coherent energy loss in cold nuclear matter, JHEP 10 (2014) 73 [arXiv:1407.5054 [hep-ph]].
- (22) R. Vogt, Shadowing effects on and production at energies available at the CERN Large Hadron Collider, Phys. Rev. C92 (2015) 034909 [arXiv:1507.04418 [hep-ph]].
- (23) J. P. Blaizot, F. Gelis and R. Venugopalan, High-energy pA collisions in the color glass condensate approach. 1. Gluon production and the Cronin effect, Nucl. Phys. A743 (2004) 13 [arXiv:hep-ph/0402256 [hep-ph]].
- (24) J. P. Blaizot, F. Gelis and R. Venugopalan, High-energy pA collisions in the color glass condensate approach. 2. Quark production, Nucl. Phys. A743 (2004) 57 [arXiv:hep-ph/0402257 [hep-ph]].
- (25) D. E. Kharzeev, E. M. Levin and K. Tuchin, Nuclear modification of the J/ transverse momentum distributions in high energy pA and AA collisions, Nucl. Phys. A924 (2014) 47 [arXiv:1205.1554 [hep-ph]].
- (26) H. Fujii, F. Gelis and R. Venugopalan, Quark pair production in high energy pA collisions: General features, Nucl. Phys. A780 (2006) 146 [arXiv:hep-ph/0603099 [hep-ph]].
- (27) H. Fujii, F. Gelis and R. Venugopalan, Quark production in high energy proton-nucleus collisions, Eur. Phys. J. C43 (2005) 139 [arXiv:hep-ph/0502204 [hep-ph]].
- (28) H. Fujii and K. Watanabe, Heavy quark pair production in high energy pA collisions: Open heavy flavors, Nucl. Phys. A920 (2013) 78 [arXiv:1308.1258 [hep-ph]].
- (29) A. D. Martin, W. J. Stirling, R. S. Thorne and G. Watt, Parton distributions for the LHC, Eur. Phys. J. C63 (2009) 189 [arXiv:0901.0002 [hep-ph]].
- (30) I. Balitsky, Operator expansion for high-energy scattering, Nucl. Phys. B463 (1996) 99 [arXiv:hep-ph/9509348 [hep-ph]].
- (31) Y. V. Kovchegov, Unitarization of the BFKL pomeron on a nucleus, Phys. Rev. D61 (2000) 074018 [arXiv:hep-ph/9905214 [hep-ph]].
- (32) I. Balitsky, Quark contribution to the small-x evolution of color dipole, Phys. Rev. D75 (2007) 014001 [arXiv:hep-ph/0609105 [hep-ph]].
- (33) T. Lappi and H. Mäntysaari, Single inclusive particle production at high energy from HERA data to proton-nucleus collisions, Phys. Rev. D88 (2013) 114020 [arXiv:1309.6963 [hep-ph]].
- (34) H1 and ZEUS collaborations, F. D. Aaron et. al., Combined Measurement and QCD Analysis of the Inclusive Scattering Cross Sections at HERA, JHEP 01 (2010) 109 [arXiv:0911.0884 [hep-ex]].
- (35) K. J. Eskola, H. Paukkunen and C. A. Salgado, EPS09 — A New Generation of NLO and LO Nuclear Parton Distribution Functions, JHEP 04 (2009) 065 [arXiv:0902.4154 [hep-ph]].
- (36) J. Pumplin, D. R. Stump, J. Huston, H. L. Lai, P. M. Nadolsky and W. K. Tung, New generation of parton distributions with uncertainties from global QCD analysis, JHEP 07 (2002) 012 [arXiv:hep-ph/0201195 [hep-ph]].
- (37) ZEUS collaboration, S. Chekanov et. al., Exclusive electroproduction of mesons at HERA, Nucl. Phys. B695 (2004) 3 [arXiv:hep-ex/0404008 [hep-ex]].
- (38) L. Frankfurt, M. Strikman and C. Weiss, Transverse nucleon structure and diagnostics of hard parton-parton processes at LHC, Phys. Rev. D83 (2011) 054012 [arXiv:1009.2559 [hep-ph]].
- (39) ALICE collaboration. Public note, to appear.
- (40) ALICE collaboration, J. Adam et. al., Centrality dependence of particle production in p-Pb collisions at = 5.02 TeV, Phys. Rev. C91 (2015) 064905 [arXiv:1412.6828 [nucl-ex]].