CMS logoCMS event Hgg
Compact Muon Solenoid
LHC, CERN

CMS-PAS-HIN-24-001
Observation of the multiplicity dependence of $ \sigma_{\psi(2S)}/\sigma_{J/\psi} $ in pPb collisions at 8.16 TeV
Abstract: The ratio of $ \psi(2S) $ and J/$ \psi $ meson production cross sections is measured as a function of the charged-particle multiplicity in pPb collisions at center-of-mass energy per nucleon pair of $ \sqrt {\smash [b]{s_{_{\mathrm {NN}}}}} = $ 8.16 TeV using the CMS detector. The data corresponds to an integrated luminosity of 175 fb$^{-1}$. Both charmonium states are measured using the dimuon decay channel. The ratio is measured for both prompt and nonprompt $ \psi(2S) $ and J/$ \psi $ mesons separately. A multiplicity-dependent modification of the ratio is observed for prompt mesons. The ratio of the nonprompt $ \psi(2S) $ to J/$ \psi $ cross sections show no dependency on multiplicity. No significant rapidity-dependence is observed within the charmonium rapidity range $ -$2.865 $ < y_{CM} < $ 1.935 and transverse momentum range probed, 3 $ < p_{T} < $ 30 GeV. These measurements constrain hadronization models of heavy quarks in nuclear collisions and support a picture in which co-moving particles may preferentially dissociate weakly-bound charmonia excited states as compared to the ground state during the evolution of the system.
Figures Summary References CMS Publications
Figures

png pdf
Figure 1:
Prompt and nonprompt normalized $ \sigma_{$ \psi (\text{2S}) $} / \sigma_{\mathrm{J}/\psi} $ for $ 1 < y_{CM} < 1.935 $ (upper left), $ -1 < y_{CM} < $ 1 (upper right), $ -2 < y_{CM} < - $1 (lower left), and $ -2.865 < y_{CM} < - 2$ (lower right) for 6.5 $ < p_{\mathrm{T}} < $ 30 GeV.

png pdf
Figure 1-a:
Prompt and nonprompt normalized $ \sigma_{$ \psi (\text{2S}) $} / \sigma_{\mathrm{J}/\psi} $ for $ 1 < y_{CM} < 1.935 $ (upper left), $ -1 < y_{CM} < $ 1 (upper right), $ -2 < y_{CM} < - $1 (lower left), and $ -2.865 < y_{CM} < - 2$ (lower right) for 6.5 $ < p_{\mathrm{T}} < $ 30 GeV.

png pdf
Figure 1-b:
Prompt and nonprompt normalized $ \sigma_{$ \psi (\text{2S}) $} / \sigma_{\mathrm{J}/\psi} $ for $ 1 < y_{CM} < 1.935 $ (upper left), $ -1 < y_{CM} < $ 1 (upper right), $ -2 < y_{CM} < - $1 (lower left), and $ -2.865 < y_{CM} < - 2$ (lower right) for 6.5 $ < p_{\mathrm{T}} < $ 30 GeV.

png pdf
Figure 1-c:
Prompt and nonprompt normalized $ \sigma_{$ \psi (\text{2S}) $} / \sigma_{\mathrm{J}/\psi} $ for $ 1 < y_{CM} < 1.935 $ (upper left), $ -1 < y_{CM} < $ 1 (upper right), $ -2 < y_{CM} < - $1 (lower left), and $ -2.865 < y_{CM} < - 2$ (lower right) for 6.5 $ < p_{\mathrm{T}} < $ 30 GeV.

png pdf
Figure 1-d:
Prompt and nonprompt normalized $ \sigma_{$ \psi (\text{2S}) $} / \sigma_{\mathrm{J}/\psi} $ for $ 1 < y_{CM} < 1.935 $ (upper left), $ -1 < y_{CM} < $ 1 (upper right), $ -2 < y_{CM} < - $1 (lower left), and $ -2.865 < y_{CM} < - 2$ (lower right) for 6.5 $ < p_{\mathrm{T}} < $ 30 GeV.

png pdf
Figure 2:
Prompt and nonprompt normalized $ \sigma_{$ \psi (\text{2S}) $} / \sigma_{\mathrm{J}/\psi} $ for $ -2.865 < y_{CM} < - 2$ (left) for 6.5 $ < p_{\mathrm{T}} < $ 30 GeV. The trend-line slope for both prompt and nonprompt double ratio is plotted (right) with respect to the $ y_{CM} $.

png pdf
Figure 2-a:
Prompt and nonprompt normalized $ \sigma_{$ \psi (\text{2S}) $} / \sigma_{\mathrm{J}/\psi} $ for $ -2.865 < y_{CM} < - 2$ (left) for 6.5 $ < p_{\mathrm{T}} < $ 30 GeV. The trend-line slope for both prompt and nonprompt double ratio is plotted (right) with respect to the $ y_{CM} $.

png pdf
Figure 2-b:
Prompt and nonprompt normalized $ \sigma_{$ \psi (\text{2S}) $} / \sigma_{\mathrm{J}/\psi} $ for $ -2.865 < y_{CM} < - 2$ (left) for 6.5 $ < p_{\mathrm{T}} < $ 30 GeV. The trend-line slope for both prompt and nonprompt double ratio is plotted (right) with respect to the $ y_{CM} $.

png pdf
Figure 3:
Dimuon invariant mass distribution showing the J/$ \psi $ and $ \psi $(2S) peaks in pPb data, after applying the $ l^{3D} \mathrm{J}/\psi $ selection, for $ -2.865 < y_{CM} < - 2$ (left) and $ 1 < y_{CM} < 1.935 $ (right) for 3 $ < p_{\mathrm{T}} < $ 30 GeV. The invariant mass fits are performed using the unbinned maximum extended likelihood method in the region of 2.5 $ < m_{\mu^{+} \mu^{-}} < $ 4.3 GeV. Both the J/$\psi$ and $\psi$(2S) nominal signal shapes are defined by the sum of two Crystal Ball (CB) functions. The background is described by a third-order Chebyshev polynomial. The total fit is shown with a solid blue line, and the dotted blue line represents the background.

png pdf
Figure 3-a:
Dimuon invariant mass distribution showing the J/$ \psi $ and $ \psi $(2S) peaks in pPb data, after applying the $ l^{3D} \mathrm{J}/\psi $ selection, for $ -2.865 < y_{CM} < - 2$ (left) and $ 1 < y_{CM} < 1.935 $ (right) for 3 $ < p_{\mathrm{T}} < $ 30 GeV. The invariant mass fits are performed using the unbinned maximum extended likelihood method in the region of 2.5 $ < m_{\mu^{+} \mu^{-}} < $ 4.3 GeV. Both the J/$\psi$ and $\psi$(2S) nominal signal shapes are defined by the sum of two Crystal Ball (CB) functions. The background is described by a third-order Chebyshev polynomial. The total fit is shown with a solid blue line, and the dotted blue line represents the background.

png pdf
Figure 3-b:
Dimuon invariant mass distribution showing the J/$ \psi $ and $ \psi $(2S) peaks in pPb data, after applying the $ l^{3D} \mathrm{J}/\psi $ selection, for $ -2.865 < y_{CM} < - 2$ (left) and $ 1 < y_{CM} < 1.935 $ (right) for 3 $ < p_{\mathrm{T}} < $ 30 GeV. The invariant mass fits are performed using the unbinned maximum extended likelihood method in the region of 2.5 $ < m_{\mu^{+} \mu^{-}} < $ 4.3 GeV. Both the J/$\psi$ and $\psi$(2S) nominal signal shapes are defined by the sum of two Crystal Ball (CB) functions. The background is described by a third-order Chebyshev polynomial. The total fit is shown with a solid blue line, and the dotted blue line represents the background.

png pdf
Figure 4:
Prompt and nonprompt normalized $ \sigma_{$ \psi (\text{2S}) $} / \sigma_{\mathrm{J}/\psi} $ for $ -2.865 < y_{CM} < - 2$ (left) and $ 1 < y_{CM} < 1.935 $ (right) for 3 $ < p_{\mathrm{T}} < $ 6.5 GeV are shown. Both prompt and nonprompt cross-section ratios are fitted with straight lines (shown in solid red for prompt and solid blue for nonprompt), considering statistical and systematic uncertainties added in quadrature. The color bands represent the 1$ \sigma $ confidence intervals. In both regions, the prompt ratio shows a dependence on multiplicity, whereas the nonprompt ratio remains independent across multiplicity. A stronger relative suppression in the Pb-going ($ -2.865 < y_{CM} < - 2$) direction compared to the p-going ($ 1 < y_{CM} < 1.935 $) direction is found for the prompt ratio, but any firm conclusion is limited by the larger uncertainty.

png pdf
Figure 4-a:
Prompt and nonprompt normalized $ \sigma_{$ \psi (\text{2S}) $} / \sigma_{\mathrm{J}/\psi} $ for $ -2.865 < y_{CM} < - 2$ (left) and $ 1 < y_{CM} < 1.935 $ (right) for 3 $ < p_{\mathrm{T}} < $ 6.5 GeV are shown. Both prompt and nonprompt cross-section ratios are fitted with straight lines (shown in solid red for prompt and solid blue for nonprompt), considering statistical and systematic uncertainties added in quadrature. The color bands represent the 1$ \sigma $ confidence intervals. In both regions, the prompt ratio shows a dependence on multiplicity, whereas the nonprompt ratio remains independent across multiplicity. A stronger relative suppression in the Pb-going ($ -2.865 < y_{CM} < - 2$) direction compared to the p-going ($ 1 < y_{CM} < 1.935 $) direction is found for the prompt ratio, but any firm conclusion is limited by the larger uncertainty.

png pdf
Figure 4-b:
Prompt and nonprompt normalized $ \sigma_{$ \psi (\text{2S}) $} / \sigma_{\mathrm{J}/\psi} $ for $ -2.865 < y_{CM} < - 2$ (left) and $ 1 < y_{CM} < 1.935 $ (right) for 3 $ < p_{\mathrm{T}} < $ 6.5 GeV are shown. Both prompt and nonprompt cross-section ratios are fitted with straight lines (shown in solid red for prompt and solid blue for nonprompt), considering statistical and systematic uncertainties added in quadrature. The color bands represent the 1$ \sigma $ confidence intervals. In both regions, the prompt ratio shows a dependence on multiplicity, whereas the nonprompt ratio remains independent across multiplicity. A stronger relative suppression in the Pb-going ($ -2.865 < y_{CM} < - 2$) direction compared to the p-going ($ 1 < y_{CM} < 1.935 $) direction is found for the prompt ratio, but any firm conclusion is limited by the larger uncertainty.
Summary
In summary, this note presents a measurement of the normalized production cross-section ratio of $ \psi $(2S) over J/$ \psi $ mesons as a function of charged-particle multiplicity in pPb collisions at center-of-mass energy per nucleon pair of $ \sqrt{\smash[b]{s_{_{\mathrm{NN}}}}} = $ 8.16 TeV with a data sample collected by the CMS detector, corresponding to an integrated luminosity of 175 nb$^{-1}$. The ratio is measured for both prompt and nonprompt $ \psi $(2S) and J/$ \psi $ mesons separately. A multiplicity-dependent modification of the ratio is observed for prompt mesons, indicating a decreasing trend with increasing multiplicity. The ratio of the nonprompt $ \psi $(2S) to J/$ \psi $ cross sections show no dependency on multiplicity. No significant rapidity-dependence is observed within the charmonium rapidity range $ -2.865 < y_{CM} < 1.935 $ and transverse momentum range, 3 $ < p_{\mathrm{T}} < $ 30 GeV probed. The observations in this study suggest that co-moving particles may dissociate weakly-bound excited state of charmonia more than the ground state. These findings provide valuable new constraints for understanding the mechanisms of hadronization in nuclear collisions involving heavy quarks.
References
1 G. T. Bodwin, E. Braaten, and G. P. Lepage Rigorous QCD analysis of inclusive annihilation and production of heavy quarkonium PRD 51 (1995) 1125 hep-ph/9407339
2 W. Busza, K. Rajagopal, and W. van der Schee Heavy ion collisions: The big picture, and the big questions Ann. Rev. Nucl. Part. Sci. 68 (2018) 339 1802.04801
3 H. Khanpour and S. Atashbar Tehrani Global analysis of nuclear parton distribution functions and their uncertainties at next-to-next-to-leading order PRD 93 (2016) 014026 1601.00939
4 F. Arleo and S. Peign é Quarkonium suppression in heavy-ion collisions from coherent energy loss in cold nuclear matter JHEP 10 (2014) 073 1407.5054
5 T. Matsui and H. Satz $ J/\psi $ suppression by quark-gluon plasma formation PLB 178 (1986) 416
6 ALICE Collaboration Centrality and transverse momentum dependence of inclusive J/\ensuremath\psi production at midrapidity in Pb-Pb collisions at $ \sqrt{\smash[b]{s_{_{\mathrm{NN}}}}} = $ 5.02 TeV PLB 805 (2020) 135434 1910.14404
7 CMS Collaboration Event Activity Dependence of Y(nS) Production in $ \sqrt {\smash [b]{s_{_{\mathrm {NN}}}}} $=5.02 TeV pPb and $ \sqrt{s} $=2.76 TeV pp Collisions JHEP 04 (2014) 103 CMS-HIN-13-003
1312.6300
8 NA50 Collaboration $ \psi' $ production in Pb-Pb collisions at 158-GeV/nucleon EPJC 49 (2007) 559 nucl-ex/0612013
9 PHENIX Collaboration Measurement of $ \psi(2S) $ nuclear modification at backward and forward rapidity in $ \mathrm{p}+\mathrm{p} $, $ \mathrm{p}+\mathrm{Al} $, and $ \mathrm{p}+\mathrm{Au} $ collisions at $ \sqrt{\smash[b]{s_{_{\mathrm{NN}}}}}= $ 200 GeV Phys. Rev. C 105 (2022) 064912 2202.03863
10 LHCb Collaboration Study of $ \Upsilon $ production in $ \mathrm{p} $Pb collisions at $ \sqrt{\smash[b]{s_{_{\mathrm{NN}}}}}= $ 8.16 TeV [Erratum: JHEP 02, 093 ()], 2018
JHEP 11 (2018) 194
1810.07655
11 ALICE Collaboration Measurement of nuclear effects on $ \psi\rm{(2S)} $ production in p-Pb collisions at $ \sqrt{\smash[b]{s_{_{\mathrm{NN}}}}} = $ 8.16 TeV JHEP 07 (2020) 237 2003.06053
12 LHCb Collaboration Study of $ \psi(2S) $ production and cold nuclear matter effects in pPb collisions at $ \sqrt{\smash[b]{s_{_{\mathrm{NN}}}}} = $ 5 TeV JHEP 03 (2016) 133 1601.07878
13 CMS Collaboration Nuclear modification of $ \Upsilon $ states in pPb collisions at $ \sqrt{\smash[b]{s_{_{\mathrm{NN}}}}} = $ 5.02 TeV PLB 835 (2022) 137397 CMS-HIN-18-005
2202.11807
14 ATLAS Collaboration Measurement of quarkonium production in proton-lead and proton-proton collisions at 5.02 $ \mathrm {TeV} $ with the ATLAS detector EPJC 78 (2018) 171 1709.03089
15 E. G. Ferreiro Charmonium dissociation and recombination at LHC: Revisiting comovers PLB 731 (2014) 57 1210.3209
16 E. G. Ferreiro Excited charmonium suppression in proton--nucleus collisions as a consequence of comovers PLB 749 (2015) 98 1411.0549
17 S. Gavin and R. Vogt Charmonium suppression by comover scattering in Pb + Pb collisions PRL 78 (1997) 1006 hep-ph/9606460
18 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? [Erratum: JHEP 03, 063 ()], 2018
JHEP 10 (2018) 094
1804.04474
19 LHCb Collaboration Multiplicity dependence of $ \sigma_{\psi(2S)}/\sigma_{J/\psi} $ in pp collisions at $ \sqrt{s}= $ 13 TeV JHEP 05 (2024) 243 2312.15201
20 ALICE Collaboration Measurement of $ \psi $(2S) production as a function of charged-particle pseudorapidity density in pp collisions at $ \sqrt{s} $ = 13 TeV and p-Pb collisions at $ \sqrt{\smash[b]{s_{_{\mathrm{NN}}}}} = $ 8.16 TeV with ALICE at the LHC JHEP 06 (2023) 147 2204.10253
21 LHCb Collaboration Prompt and nonprompt \ensuremath\psi(2S) production in pPb collisions at $ \sqrt{s_{\textrm{NN}}} $ = 8.16 TeV JHEP 04 (2024) 111 2401.11342
22 CMS Collaboration The CMS experiment at the CERN LHC JINST 3 (2008) S08004
23 CMS Collaboration Development of the CMS detector for the CERN LHC Run 3 JINST 19 (2024) P05064 CMS-PRF-21-001
2309.05466
24 CMS Collaboration Performance of the CMS muon detector and muon reconstruction with proton-proton collisions at $ \sqrt{s}= $ 13 TeV JINST 13 (2018) P06015 CMS-MUO-16-001
1804.04528
25 CMS Collaboration Performance of the CMS Level-1 trigger in proton-proton collisions at $ \sqrt{s} = $ 13 TeV JINST 15 (2020) P10017 CMS-TRG-17-001
2006.10165
26 CMS Collaboration The CMS trigger system JINST 12 (2017) P01020 CMS-TRG-12-001
1609.02366
27 CMS Collaboration Observation of correlated azimuthal anisotropy fourier harmonics in pp and $ \mathrm{p}+\mathrm{Pb} $ collisions at the LHC PRL 120 (2018) 092301 CMS-HIN-16-022
1709.09189
28 CMS Collaboration Constraints on the chiral magnetic effect using charge-dependent azimuthal correlations in $ \mathrm{p}\mathrm{Pb} $ and PbPb collisions at the CERN Large Hadron Collider Phys. Rev. C 97 (2018) 044912 CMS-HIN-17-001
1708.01602
29 CMS Collaboration Description and performance of track and primary-vertex reconstruction with the CMS tracker JINST 9 (2014) P10009 CMS-TRK-11-001
1405.6569
30 C. Bierlich et al. A comprehensive guide to the physics and usage of PYTHIA 8.3 SciPost Phys. Codeb. 2022 (2022) 8 2203.11601
31 CMS Collaboration Event generator tunes obtained from underlying event and multiparton scattering measurements EPJC 76 (2016) 155 CMS-GEN-14-001
1512.00815
32 D. J. Lange The EvtGen particle decay simulation package NIM A 462 (2001) 152
33 N. Davidson, T. Przedzinski, and Z. Was PHOTOS interface in C++: Technical and physics documentation Comput. Phys. Commun. 199 (2016) 86 1011.0937
34 T. Pierog et al. EPOS LHC: Test of collective hadronization with data measured at the CERN Large Hadron Collider Phys. Rev. C 92 (2015) 034906 1306.0121
35 GEANT4 Collaboration GEANT 4---a simulation toolkit NIM A 506 (2003) 250
36 CMS Collaboration Measurement of the inclusive W and Z production cross sections in pp collisions at $ \sqrt{s}= $ 7 TeV JHEP 10 (2011) 132 CMS-EWK-10-005
1107.4789
37 CMS Collaboration Performance of CMS muon reconstruction in heavy ion collisions CMS Physics Analysis Summary, 2023
link
CMS-PAS-MUO-21-001
38 Particle Data Group Collaboration Review of particle physics PTEP 2022 (2022) 083C01
39 CMS Collaboration Measurement of prompt $ \psi(2S) \to J/\psi $ yield ratios in Pb-Pb and $ \mathrm{p}-\mathrm{p} $ collisions at $ \sqrt{\smash[b]{s_{_{\mathrm{NN}}}}} = $ 2.76 TeV PRL 113 (2014) 262301 CMS-HIN-12-007
1410.1804
40 CMS Collaboration Relative modification of prompt \ensuremath\psi(2S) and J/\ensuremath\psi yields from pp to PbPb collisions at $ \sqrt{\smash[b]{s_{_{\mathrm{NN}}}}}= $ 5.02 TeV PRL 118 (2017) 162301 CMS-HIN-16-004
1611.01438
41 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 PLB 790 (2019) 509 CMS-HIN-16-015
1805.02248
42 M. Oreglia A study of the reactions $ \psi' \to \gamma\gamma\psi $ Phd thesis, Stanford University, 1980
Compact Muon Solenoid
LHC, CERN