CMS logoCMS event Hgg
Compact Muon Solenoid
LHC, CERN

CMS-PAS-HIN-18-005
Nuclear modification of $\Upsilon$ states in pPb collisions at ${\sqrt {\smash [b]{s_{_{\mathrm {NN}}}}}} = $ 5.02 TeV
Abstract: Production cross sections of $\Upsilon(\mathrm{1S})$, $\Upsilon(\mathrm{2S})$, and $\Upsilon(\mathrm{3S})$ states decaying into $\mu^{+}\mu^{-}$ in proton-lead (pPb) collisions are reported using data collected by the CMS experiment at ${\sqrt {\smash [b]{s_{_{\mathrm {NN}}}}}}= $ 5.02 TeV. Nuclear modification factors $R_\mathrm{pPb}$ for all three $\Upsilon$ states are obtained using measured proton-proton (pp) cross sections at the same collision energy. All $\Upsilon$ states are found to be suppressed in pPb collisions compared to pp collisions. The $\Upsilon$ $R_\mathrm{pPb}$ show a sequential ordering, with $\Upsilon(\mathrm{1S})$ least suppressed and $\Upsilon(\mathrm{3S})$ most suppressed, indicating presence of final-state modification of $\Upsilon$ states in pPb collisions. When presented as a function of transverse momentum $p_\mathrm{T}$ and center-of-mass rapidity $y_\mathrm{CM}$, the $R_\mathrm{pPb}$ of individual $\Upsilon$ states are found to be consistent with constant values; although there is slight indication of higher separation of the suppression level of excited states for low-$p_\mathrm{T}$ $\Upsilon$ in the lead-going direction, where more nuclear matter is present. The final-state comover interaction model, which predicts sequential suppression of bottomonia in pPb, is found to be in better agreement with $R_\mathrm{pPb}$ versus $y_\mathrm{CM}$ than initial-state models. The nuclear modification observed in pPb collisions is less pronounced than the strong modification observed in lead-lead collisions, suggesting presence of additional quark gluon plasma effects in the latter. Forward-backward production ratios $R_\mathrm{FB}$ of $\Upsilon$ states are reported as functions of event activity variables, obtained from midrapidity as well as forward and backward rapidity regions. The $R_\mathrm{FB}$ for all $\Upsilon$ states are found to be consistent with unity and constant with increasing event activity, irrespective of the rapidity region used to measure activity.
Figures Summary References CMS Publications
Figures

png pdf
Figure 1:
Measured dimuon invariant mass distributions (closed circles) for pp collisions (left) and pPb collisions (right). The total unbinned maximum-likelihood fits to the data are shown as solid blue lines, with the background component indicated by dashed blue lines. The individual $\Upsilon (\mathrm {1S})$, $\Upsilon (\mathrm {2S})$, and $\Upsilon (\mathrm {3S})$ signal shapes in pp are depicted as dashed green lines in the left panel. The dashed red line in the right panel is obtained by scaling the $\Upsilon (\mathrm {1S})$, $\Upsilon (\mathrm {2S})$, and $\Upsilon (\mathrm {3S})$ signal shapes in pPb (solid blue line) by the inverse of the measured ${R_{\mathrm {pPb}}}$ for each state.

png pdf
Figure 1-a:
Measured dimuon invariant mass distribution (closed circles) for pp collisions. The total unbinned maximum-likelihood fits to the data are shown as solid blue lines, with the background component indicated by dashed blue lines. The individual $\Upsilon (\mathrm {1S})$, $\Upsilon (\mathrm {2S})$, and $\Upsilon (\mathrm {3S})$ signal shapes in pp are depicted as dashed green lines in the left panel. The dashed red line in the right panel is obtained by scaling the $\Upsilon (\mathrm {1S})$, $\Upsilon (\mathrm {2S})$, and $\Upsilon (\mathrm {3S})$ signal shapes in pPb (solid blue line) by the inverse of the measured ${R_{\mathrm {pPb}}}$ for each state.

png pdf
Figure 1-b:
Measured dimuon invariant mass distribution (closed circles) for pPb collisions. The total unbinned maximum-likelihood fits to the data are shown as solid blue lines, with the background component indicated by dashed blue lines. The individual $\Upsilon (\mathrm {1S})$, $\Upsilon (\mathrm {2S})$, and $\Upsilon (\mathrm {3S})$ signal shapes in pp are depicted as dashed green lines in the left panel. The dashed red line in the right panel is obtained by scaling the $\Upsilon (\mathrm {1S})$, $\Upsilon (\mathrm {2S})$, and $\Upsilon (\mathrm {3S})$ signal shapes in pPb (solid blue line) by the inverse of the measured ${R_{\mathrm {pPb}}}$ for each state.

png pdf
Figure 2:
Cross section times dimuon branching fraction of $\Upsilon (\mathrm {1S})$ (red circles), $\Upsilon (\mathrm {2S})$ (blue squares), and $\Upsilon (\mathrm {3S})$ (green diamonds) as a function of $p_\mathrm {T}$ (left) and rapidity (right) in pPb collisions. Error bars on the points represent statistical and fit uncertainties and filled boxes represent systematic uncertainties.

png pdf
Figure 2-a:
Cross section times dimuon branching fraction of $\Upsilon (\mathrm {1S})$ (red circles), $\Upsilon (\mathrm {2S})$ (blue squares), and $\Upsilon (\mathrm {3S})$ (green diamonds) as a function of $p_\mathrm {T}$ in pPb collisions. Error bars on the points represent statistical and fit uncertainties and filled boxes represent systematic uncertainties.

png pdf
Figure 2-b:
Cross section times dimuon branching fraction of $\Upsilon (\mathrm {1S})$ (red circles), $\Upsilon (\mathrm {2S})$ (blue squares), and $\Upsilon (\mathrm {3S})$ (green diamonds) as a function of rapidity in pPb collisions. Error bars on the points represent statistical and fit uncertainties and filled boxes represent systematic uncertainties.

png pdf
Figure 3:
${R_{\mathrm {pPb}}}$ of $\Upsilon (\mathrm {1S})$ (red circles), $\Upsilon (\mathrm {2S})$ (blue squares), and $\Upsilon (\mathrm {3S})$ (green diamonds) as a function of $p_\mathrm {T}$ for $| {y_{\mathrm {CM}}}| < $ 1.93 (left) and versus ${y_{\mathrm {CM}}}$ for $p_\mathrm {T} < $ 30 GeV/$c$ (right). Error bars on the points represent statistical and fit uncertainties and filled boxes represent systematic uncertainties. The gray box around the line at unity represents the global uncertainty due to luminosity normalization. All three $\Upsilon $ states are suppressed in pPb collisions compared to pp collisions throughout the kinematic region explored. For each $\Upsilon $ state, the measured ${R_{\mathrm {pPb}}}$ is consistent with a constant value across the kinematic range. The $\Upsilon $ states show a sequential pattern of suppression, with $\Upsilon (\mathrm {1S})$ the least suppressed.

png pdf
Figure 3-a:
${R_{\mathrm {pPb}}}$ of $\Upsilon (\mathrm {1S})$ (red circles), $\Upsilon (\mathrm {2S})$ (blue squares), and $\Upsilon (\mathrm {3S})$ (green diamonds) as a function of $p_\mathrm {T}$ for $| {y_{\mathrm {CM}}}| < $ 1.93. Error bars on the points represent statistical and fit uncertainties and filled boxes represent systematic uncertainties. The gray box around the line at unity represents the global uncertainty due to luminosity normalization. All three $\Upsilon $ states are suppressed in pPb collisions compared to pp collisions throughout the kinematic region explored. For each $\Upsilon $ state, the measured ${R_{\mathrm {pPb}}}$ is consistent with a constant value across the kinematic range. The $\Upsilon $ states show a sequential pattern of suppression, with $\Upsilon (\mathrm {1S})$ the least suppressed.

png pdf
Figure 3-b:
${R_{\mathrm {pPb}}}$ of $\Upsilon (\mathrm {1S})$ (red circles), $\Upsilon (\mathrm {2S})$ (blue squares), and $\Upsilon (\mathrm {3S})$ (green diamonds) versus ${y_{\mathrm {CM}}}$ for $p_\mathrm {T} < $ 30 GeV/$c$. Error bars on the points represent statistical and fit uncertainties and filled boxes represent systematic uncertainties. The gray box around the line at unity represents the global uncertainty due to luminosity normalization. All three $\Upsilon $ states are suppressed in pPb collisions compared to pp collisions throughout the kinematic region explored. For each $\Upsilon $ state, the measured ${R_{\mathrm {pPb}}}$ is consistent with a constant value across the kinematic range. The $\Upsilon $ states show a sequential pattern of suppression, with $\Upsilon (\mathrm {1S})$ the least suppressed.

png pdf
Figure 4:
${R_{\mathrm {pPb}}}$ of $\Upsilon (\mathrm {1S})$ versus ${y_{\mathrm {CM}}}$ with shadowing predictions from R. Vogt [23] (left) and energy loss with and without shadowing corrections from F. Arleo and S. Peigne [24] (right). These initial- and intermediate-state effects modify all $\Upsilon $ states similarly. Error bars on the points represent statistical and fit uncertainties and filled boxes represent systematic uncertainties. The gray box around the line at unity represents the global uncertainty due to luminosity normalization.

png pdf
Figure 4-a:
${R_{\mathrm {pPb}}}$ of $\Upsilon (\mathrm {1S})$ versus ${y_{\mathrm {CM}}}$ with shadowing predictions from R. Vogt [23]. Error bars on the points represent statistical and fit uncertainties and filled boxes represent systematic uncertainties. The gray box around the line at unity represents the global uncertainty due to luminosity normalization.

png pdf
Figure 4-b:
${R_{\mathrm {pPb}}}$ of $\Upsilon (\mathrm {1S})$ versus ${y_{\mathrm {CM}}}$ with energy loss with and without shadowing corrections from F. Arleo and S. Peigne [24]. Error bars on the points represent statistical and fit uncertainties and filled boxes represent systematic uncertainties. The gray box around the line at unity represents the global uncertainty due to luminosity normalization.

png pdf
Figure 5:
${R_{\mathrm {pPb}}}$ versus ${y_{\mathrm {CM}}}$ with comover effect predictions from E. Ferreiro and J. Lansberg [26] with shadowing corrections using nCTEQ15 and EPS09 for $\Upsilon (\mathrm {1S})$ (top left), $\Upsilon (\mathrm {2S})$ (top right) and $\Upsilon (\mathrm {3S})$ (bottom). The final-state comover effect is seen to modify the $\Upsilon $ states sequentially. Error bars on the points represent statistical and fit uncertainties and filled boxes represent systematic uncertainties. The gray box around the line at unity represents the global uncertainty due to luminosity normalization.

png pdf
Figure 5-a:
${R_{\mathrm {pPb}}}$ versus ${y_{\mathrm {CM}}}$ with comover effect predictions from E. Ferreiro and J. Lansberg [26] with shadowing corrections using nCTEQ15 and EPS09 for $\Upsilon (\mathrm {1S})$. Error bars on the points represent statistical and fit uncertainties and filled boxes represent systematic uncertainties. The gray box around the line at unity represents the global uncertainty due to luminosity normalization.

png pdf
Figure 5-b:
${R_{\mathrm {pPb}}}$ versus ${y_{\mathrm {CM}}}$ with comover effect predictions from E. Ferreiro and J. Lansberg [26] with shadowing corrections using nCTEQ15 and EPS09 for $\Upsilon (\mathrm {2S})$. Error bars on the points represent statistical and fit uncertainties and filled boxes represent systematic uncertainties. The gray box around the line at unity represents the global uncertainty due to luminosity normalization.

png pdf
Figure 5-c:
${R_{\mathrm {pPb}}}$ versus ${y_{\mathrm {CM}}}$ with comover effect predictions from E. Ferreiro and J. Lansberg [26] with shadowing corrections using nCTEQ15 and EPS09 for $\Upsilon (\mathrm {3S})$. Error bars on the points represent statistical and fit uncertainties and filled boxes represent systematic uncertainties. The gray box around the line at unity represents the global uncertainty due to luminosity normalization.

png pdf
Figure 6:
${R_{\mathrm {pPb}}}$ of $\Upsilon (\mathrm {1S})$ (red circles), $\Upsilon (\mathrm {2S})$ (blue squares), and $\Upsilon (\mathrm {3S})$ (green diamonds) at forward and backward rapidity for 0 $ < {p_{\mathrm {T}}} < $ 6 GeV/$c$ (left) and 6 $ < {p_{\mathrm {T}}} < $ 30 GeV/$c$ (right). The points are shifted horizontally for better visibility. Error bars on the points represent statistical and fit uncertainties and filled boxes represent systematic uncertainties. The gray box around the line at unity represents the global uncertainty due to luminosity normalization.

png pdf
Figure 6-a:
${R_{\mathrm {pPb}}}$ of $\Upsilon (\mathrm {1S})$ (red circles), $\Upsilon (\mathrm {2S})$ (blue squares), and $\Upsilon (\mathrm {3S})$ (green diamonds) at forward and backward rapidity for 0 $ < {p_{\mathrm {T}}} < $ 6 GeV/$c$. The points are shifted horizontally for better visibility. Error bars on the points represent statistical and fit uncertainties and filled boxes represent systematic uncertainties. The gray box around the line at unity represents the global uncertainty due to luminosity normalization.

png pdf
Figure 6-b:
${R_{\mathrm {pPb}}}$ of $\Upsilon (\mathrm {1S})$ (red circles), $\Upsilon (\mathrm {2S})$ (blue squares), and $\Upsilon (\mathrm {3S})$ (green diamonds) at forward and backward rapidity for 6 $ < {p_{\mathrm {T}}} < $ 30 GeV/$c$. The points are shifted horizontally for better visibility. Error bars on the points represent statistical and fit uncertainties and filled boxes represent systematic uncertainties. The gray box around the line at unity represents the global uncertainty due to luminosity normalization.

png pdf
Figure 7:
${R_{\mathrm {FB}}}$ vs. midrapidity $N_\mathrm {tracks}$ (left) and forward/backward rapidity $E_{T}^{HF}$ (right) of $\Upsilon (\mathrm {1S})$ (red circles), $\Upsilon (\mathrm {2S})$ (blue squares), and $\Upsilon (\mathrm {3S})$ (green diamonds) for $p_\mathrm {T} < $ 30 GeV/$c$ and $| {y_{\mathrm {CM}}}| < $ 1.93. Error bars on the points represent statistical and fit uncertainties and filled boxes represent systematic uncertainties.

png pdf
Figure 7-a:
${R_{\mathrm {FB}}}$ vs. midrapidity $N_\mathrm {tracks}$ of $\Upsilon (\mathrm {1S})$ (red circles), $\Upsilon (\mathrm {2S})$ (blue squares), and $\Upsilon (\mathrm {3S})$ (green diamonds) for $p_\mathrm {T} < $ 30 GeV/$c$ and $| {y_{\mathrm {CM}}}| < $ 1.93. Error bars on the points represent statistical and fit uncertainties and filled boxes represent systematic uncertainties.

png pdf
Figure 7-b:
${R_{\mathrm {FB}}}$ vs. forward/backward rapidity $E_{T}^{HF}$ of $\Upsilon (\mathrm {1S})$ (red circles), $\Upsilon (\mathrm {2S})$ (blue squares), and $\Upsilon (\mathrm {3S})$ (green diamonds) for $p_\mathrm {T} < $ 30 GeV/$c$ and $| {y_{\mathrm {CM}}}| < $ 1.93. Error bars on the points represent statistical and fit uncertainties and filled boxes represent systematic uncertainties.

png pdf
Figure 8:
${R_{\mathrm {pPb}}}$ of $\Upsilon (\mathrm {1S})$, $\Upsilon (\mathrm {2S})$ and $\Upsilon (\mathrm {3S})$ (red circles) for the integrated kinematic range 0 $ < {p_{\mathrm {T}}} < $ 30 GeV/$c$ and $| {y_{\mathrm {CM}}}| < $ 1.93. The ${R_{\mathrm {pPb}}}$ results are compared to the CMS results on $\Upsilon $ ${R_{\mathrm {AA}}}$ (blue squares for $\Upsilon (\mathrm {1S})$ and $\Upsilon (\mathrm {2S})$ and blue arrow for $\Upsilon (\mathrm {3S})$ at 95% confidence level) for 0 $ < {p_{\mathrm {T}}} < $ 30 GeV/$c$ and $| {y_{\mathrm {CM}}}| < $ 2.4 at the same energy [17]. Error bars represent statistical and fit uncertainties and filled boxes around points represent systematic uncertainties. The gray and red boxes around the line at unity depict the uncertainty in the pp and pPb luminosity normalizations, respectively. The blue box around unity depicts the global uncertainty pertaining to PbPb data.

png pdf
Figure 9:
Cross section times dimuon branching fraction of $\Upsilon (\mathrm {1S})$ (red circles), $\Upsilon (\mathrm {2S})$ (blue squares), and $\Upsilon (\mathrm {3S})$ (green diamonds) as a function of $p_\mathrm {T}$ (left) and rapidity (right) in pp collisions. Error bars on the points represent statistical and fit uncertainties and filled boxes represent systematic uncertainties.

png pdf
Figure 9-a:
Cross section times dimuon branching fraction of $\Upsilon (\mathrm {1S})$ (red circles), $\Upsilon (\mathrm {2S})$ (blue squares), and $\Upsilon (\mathrm {3S})$ (green diamonds) as a function of $p_\mathrm {T}$ in pp collisions. Error bars on the points represent statistical and fit uncertainties and filled boxes represent systematic uncertainties.

png pdf
Figure 9-b:
Cross section times dimuon branching fraction of $\Upsilon (\mathrm {1S})$ (red circles), $\Upsilon (\mathrm {2S})$ (blue squares), and $\Upsilon (\mathrm {3S})$ (green diamonds) as a function of rapidity in pp collisions. Error bars on the points represent statistical and fit uncertainties and filled boxes represent systematic uncertainties.
Summary
In summary, the $\Upsilon$ family has been studied in pPb collisions at 5.02 TeV and the production cross sections presented. Together with pp data obtained at the same center-of-mass energy, we measured the nuclear modification factors for the three $\Upsilon$ states decaying in the dimuon channel. We observe a suppression of the $\Upsilon$ yields relative to the hypothesis of $A$-scaling for all three states, in the full kinematic range studied. No significant trend is seen for the suppression as functions of $p_\mathrm{T}$ or $y_\mathrm{CM}$, although there is slight indication of higher separation of the suppression level of excited states in the lead-going direction for low-$p_\mathrm{T}$ $\Upsilon$. The forward-backward production ratios $R_\mathrm{FB}$ of $\Upsilon$ states were studied separately as functions of event activity recorded near to and far away from the rapidity region where $\Upsilon$ mesons were measured. The $R_\mathrm{FB}$ values are consistent with unity for all states, independent of the rapidity region used to measure activity. The integrated nuclear modification factors for $\Upsilon$ were compared to those obtained in PbPb collisions. The $R_\mathrm{AA}$ values are much smaller than the corresponding $R_\mathrm{pPb}$ value for each state, as expected in the presence of deconfinement effects in PbPb. A modest sequential suppression, consistent with predictions from hadronic comover effects, is observed in pPb, indicating the presence of final-state effects in pPb collisions. These results will help to elucidate the contributions of cold and hot nuclear matter effects to the modification of bottomonia in heavy ion collisions.
References
1 T. Matsui and H. Satz $ J/\psi $ Suppression by Quark-Gluon Plasma Formation PLB178 (1986) 416
2 J. W. Harris and B. Muller The search for the quark gluon plasma Ann.Rev.Nucl.Part.Sci 46 (1996) 71
3 J. F. Gunion and R. Vogt Determining the existence and nature of the quark-gluon plasma by Upsilon suppression at the LHC Nucl.Phys. B 492 (1997) 301
4 S. Kim, P. Petreczky, and A. Rothkopf In-medium quarkonium properties from a lattice QCD based effective field theory NP A956 (2016) 713 1512.05289
5 A. Emerick, X. Zhao, and R. Rapp Bottomonia in the Quark-Gluon Plasma and their Produciton at RHIC and LHC EPJA 48 (2012)
6 A. Rothkopf What lattice QCD spectral functions can tell us about heavy quarkonium in the QGP PoS ICHEP2016 (2016) 362 1611.06517
7 M. Laine, O. Philipsen, P. Romatschke, and M. Tassler Real-time static potential in hot QCD JHEP 03 (2007) 054 hep-ph/0611300
8 J.-P. Blaizot, D. De Boni, P. Faccioli, and G. Garberoglio Heavy quark bound states in a quark-gluon plasma: Dissociation and recombination NP A946 (2016) 49 1503.03857
9 S. Chen and M. He Gluo-dissociation of heavy quarkonium in the quark-gluon plasma reexamined PRC96 (2017) 034901 1705.10110
10 A. Andronic et al. Heavy-flavour and quarkonium production in the LHC era: from proton-proton to heavy-ion collisions Phys. J. C 76 (2016)
11 S. Digal, P. Petreczky, and H. Satz Quarkonium Feed-Down and Sequential Suppression Phys.Rev. D 64 (2001)
12 R. Vogt Cold nuclear matter effects on J/$ \psi $ and $ \Upsilon $ production at the LHC Phys.Rev.C 81 (2010)
13 M. Strickland Thermal Bottomonium Suppression AIP Conf. Proc. (2013) 179
14 F. Arleo and S. Peigne Heavy-quarkonium suppression in p-A collisions from parton energy loss in cold QCD matter JHEP 03 (2013) 122 1212.0434
15 CMS Collaboration Suppression of $ \Upsilon(1S), \Upsilon(2S) $ and $ \Upsilon(3S) $ production in PbPb collisions at $ {\sqrt {\smash [b]{s_{_{\mathrm {NN}}}}}} = $ 2.76 TeV PLB770 (2017) 357 CMS-HIN-15-001
1611.01510
16 CMS Collaboration Suppression of Excited $ \Upsilon $ States Relative to the Ground State in Pb-Pb Collisions at $ {\sqrt {\smash [b]{s_{_{\mathrm {NN}}}}}} = $ 5.02 TeV PRL 120 (2018) 142301 CMS-HIN-16-008
1706.05984
17 CMS Collaboration Measurement of nuclear modification factors of $ \Upsilon $(1S), $ \Upsilon $(2S), and $ \Upsilon $(3S) mesons in PbPb collisions at $ {\sqrt {\smash [b]{s_{_{\mathrm {NN}}}}}} = $ 5.02 TeV CMS-HIN-16-023
1805.09215
18 N. Brambilla et al. Heavy quarkonium: progress, puzzles and opportunities Eur.Phys.J. C 71 (2011)
19 R. Arnaldi Experimental overview on quarkonium production Nuclear Physics A 956 (2016) 128
20 A. Accardi et al. Parton Propagation and Fragmentation in QCD Matter Riv. Nuovo Cim. 32 (2010) 439 0907.3534
21 M. C. Chu and T. Matsui Dynamic Debye screening for a heavy-quark-anti-quark pair traversing a quark - gluon plasma PRD39 (1989) 1892
22 H. Fujii, F. Gelis, and R. Venugopalan Quark pair production in high energy pA collisions: General features NP A780 (2006) 146 hep-ph/0603099
23 R. Vogt Shadowing effects on $ J/\psi $ and $ \Upsilon $ production at energies available at the CERN Large Hadron Collider PRC92 (2015) 034909 1507.04418
24 F. Arleo and S. Peigne Quarkonium suppression in heavy-ion collisions from coherent energy loss in cold nuclear matter JHEP 10 (2014) 073 1407.5054
25 E. G. Ferreiro Excited charmonium suppression in proton-nucleus collisions as a consequence of comovers PLB749 (2015) 98 1411.0549
26 E. G. Ferreiro and J.-P. Lansberg Is bottomonium suppression in proton-nucleus and nucleus-nucleus collisions at LHC energies due to the same effects? 1804.04474
27 CMS Collaboration Performance of CMS muon reconstruction in $ pp $ collision events at $ \sqrt{s}= $ 7 TeV JINST 7 (2012) P10002 CMS-MUO-10-004
1206.4071
28 CMS Collaboration The CMS experiment at the CERN LHC JINST 3 (2008) S08004 CMS-00-001
29 R. Kalman A New Approach To Linear Filtering and Prediction Problems J. Basic Eng 82D (1960) 35
30 M. Tanabashi et al. Particle Data Group PRD 98 (2018) 030001
31 T. Sjostrand et al. An Introduction to PYTHIA 8.2 CPC 191 (2015) 159 1410.3012
32 GEANT4 Collaboration GEANT4: A Simulation toolkit NIMA506 (2003) 250
33 ATLAS Collaboration Measurement of quarkonium production in proton-lead and proton-proton collisions at 5.02 TeV with the ATLAS detector EPJC78 (2018) 171 1709.03089
34 CMS Collaboration Measurement of prompt $ \psi $(2S) production cross sections in proton-lead and proton-proton collisions at $ {\sqrt {\smash [b]{s_{_{\mathrm {NN}}}}}}= $ 5.02 TeV CMS-HIN-16-015
1805.02248
35 CMS Collaboration Measurement of prompt and nonprompt $ \mathrm{J}/{\psi} $ production in $ \mathrm {p}\mathrm {p} $ and $ \mathrm {p}\mathrm {Pb} $ collisions at $ {\sqrt {\smash [b]{s_{_{\mathrm {NN}}}}}} = $ 5.02 TeV EPJC77 (2017) 269 CMS-HIN-14-009
1702.01462
36 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 0902.4154
37 K. Kovarik et al. nCTEQ15 - Global analysis of nuclear parton distributions with uncertainties in the CTEQ framework PRD93 (2016) 085037 1509.00792
38 ALICE Collaboration Production of inclusive $ \Upsilon $(1S) and $ \Upsilon $(2S) in p-Pb collisions at $ {\sqrt {\smash [b]{s_{_{\mathrm {NN}}}}}} = $ 5.02 TeV PLB740 (2015) 105 1410.2234
39 LHCb Collaboration Study of $ \Upsilon $ production and cold nuclear matter effects in pPb collisions at $ {\sqrt {\smash [b]{s_{_{\mathrm {NN}}}}}} = $ 5 TeV JHEP 07 (2014) 094 1405.5152
Compact Muon Solenoid
LHC, CERN