CMS logoCMS event Hgg
Compact Muon Solenoid
LHC, CERN

CMS-SMP-25-005 ; CERN-EP-2026-183
Measurement of the fragmentation properties of jets containing $ \Upsilon{\textrm{(nS)}} $ mesons in proton-proton collisions at $ \sqrt{s}= $ 13 TeV
Submitted to the Journal of High Energy Physics
Abstract: A measurement of the fragmentation properties of jets containing $ \Upsilon{\textrm{(nS)}} $ mesons using proton-proton collision data at $ \sqrt{s}= $ 13 TeV, corresponding to an integrated luminosity of 138 fb$ ^{-1} $, is presented. The $ \Upsilon{\textrm{(nS)}} $ mesons associated with jets are reconstructed via their decays producing pairs of oppositely charged muons. The longitudinal and transverse projections of the momenta of the $ \Upsilon{\textrm{(nS)}} $ mesons along the momenta of the corresponding jets are studied. The results are compared with Monte Carlo predictions, including recent developments in the modelling of quarkonia production in parton showers. The description of the data by these Monte Carlo predictions is found to be unsatisfactory, leaving room for improvements in the modelling of such processes.
Figures Summary References CMS Publications
Figures

png pdf
Figure 1:
Example diagrams for the hard production of $ \Upsilon{\textrm{(nS)}} $ mesons in the colour singlet (left) and colour octet mechanisms (right).

png pdf
Figure 1-a:
Example diagrams for the hard production of $ \Upsilon{\textrm{(nS)}} $ mesons in the colour singlet (left) and colour octet mechanisms (right).

png pdf
Figure 1-b:
Example diagrams for the hard production of $ \Upsilon{\textrm{(nS)}} $ mesons in the colour singlet (left) and colour octet mechanisms (right).

png pdf
Figure 2:
Fits to the invariant mass distributions of $ \Upsilon{\textrm{(nS)}} $ candidates in two representative ranges of the jet $ p_{\mathrm{T}} $ and the fragmentation variables $ z $ (left) and $ p_{\mathrm{T}}^{\text{rel}} $ (right). The fit components, including signal and background functions, are also shown on the upper panels. The vertical error bars represent the statistical uncertainties on the data points. The lower panel shows the difference between the data and the fitted function, divided by the statistical uncertainty of the data. The model shows, in general, a good agreement with the data.

png pdf
Figure 2-a:
Fits to the invariant mass distributions of $ \Upsilon{\textrm{(nS)}} $ candidates in two representative ranges of the jet $ p_{\mathrm{T}} $ and the fragmentation variables $ z $ (left) and $ p_{\mathrm{T}}^{\text{rel}} $ (right). The fit components, including signal and background functions, are also shown on the upper panels. The vertical error bars represent the statistical uncertainties on the data points. The lower panel shows the difference between the data and the fitted function, divided by the statistical uncertainty of the data. The model shows, in general, a good agreement with the data.

png pdf
Figure 2-b:
Fits to the invariant mass distributions of $ \Upsilon{\textrm{(nS)}} $ candidates in two representative ranges of the jet $ p_{\mathrm{T}} $ and the fragmentation variables $ z $ (left) and $ p_{\mathrm{T}}^{\text{rel}} $ (right). The fit components, including signal and background functions, are also shown on the upper panels. The vertical error bars represent the statistical uncertainties on the data points. The lower panel shows the difference between the data and the fitted function, divided by the statistical uncertainty of the data. The model shows, in general, a good agreement with the data.

png pdf
Figure 3:
Detector-level longitudinal profile for $ \Upsilon{\textrm{(1S)}} $ mesons (left) and transverse profile for $ \Upsilon{\textrm{(3S)}} $ mesons (right) for two different ranges of the jet transverse momentum. The error bars on the data distributions represent the statistical uncertainties.

png pdf
Figure 3-a:
Detector-level longitudinal profile for $ \Upsilon{\textrm{(1S)}} $ mesons (left) and transverse profile for $ \Upsilon{\textrm{(3S)}} $ mesons (right) for two different ranges of the jet transverse momentum. The error bars on the data distributions represent the statistical uncertainties.

png pdf
Figure 3-b:
Detector-level longitudinal profile for $ \Upsilon{\textrm{(1S)}} $ mesons (left) and transverse profile for $ \Upsilon{\textrm{(3S)}} $ mesons (right) for two different ranges of the jet transverse momentum. The error bars on the data distributions represent the statistical uncertainties.

png pdf
Figure 4:
Transfer matrices for the longitudinal profile of jets containing $ \Upsilon{\textrm{(1S)}} $ mesons (left) and for the transverse profile of jets containing $ \Upsilon{\textrm{(3S)}} $ (right). Each transfer matrix parameterizes the migrations among the three bins on the jet $ p_{\mathrm{T}} $ used in the measurement, shown as boxes delimited by dashed lines, as well as among bins on $ z $ or $ p_{\mathrm{T}}^{\text{rel}} $.

png pdf
Figure 4-a:
Transfer matrices for the longitudinal profile of jets containing $ \Upsilon{\textrm{(1S)}} $ mesons (left) and for the transverse profile of jets containing $ \Upsilon{\textrm{(3S)}} $ (right). Each transfer matrix parameterizes the migrations among the three bins on the jet $ p_{\mathrm{T}} $ used in the measurement, shown as boxes delimited by dashed lines, as well as among bins on $ z $ or $ p_{\mathrm{T}}^{\text{rel}} $.

png pdf
Figure 4-b:
Transfer matrices for the longitudinal profile of jets containing $ \Upsilon{\textrm{(1S)}} $ mesons (left) and for the transverse profile of jets containing $ \Upsilon{\textrm{(3S)}} $ (right). Each transfer matrix parameterizes the migrations among the three bins on the jet $ p_{\mathrm{T}} $ used in the measurement, shown as boxes delimited by dashed lines, as well as among bins on $ z $ or $ p_{\mathrm{T}}^{\text{rel}} $.

png pdf
Figure 5:
Relative values of the systematic uncertainties discussed in the text, together with the total uncertainty band, for the longitudinal and transverse profiles of $ \Upsilon{\textrm{(1S)}} $ (top), $ \Upsilon{\textrm{(2S)}} $ (middle) and $ \Upsilon{\textrm{(3S)}} $ mesons (bottom) in jets with 60 $ < p_{\mathrm{T}} < $ 75 GeV.

png pdf
Figure 5-a:
Relative values of the systematic uncertainties discussed in the text, together with the total uncertainty band, for the longitudinal and transverse profiles of $ \Upsilon{\textrm{(1S)}} $ (top), $ \Upsilon{\textrm{(2S)}} $ (middle) and $ \Upsilon{\textrm{(3S)}} $ mesons (bottom) in jets with 60 $ < p_{\mathrm{T}} < $ 75 GeV.

png pdf
Figure 5-b:
Relative values of the systematic uncertainties discussed in the text, together with the total uncertainty band, for the longitudinal and transverse profiles of $ \Upsilon{\textrm{(1S)}} $ (top), $ \Upsilon{\textrm{(2S)}} $ (middle) and $ \Upsilon{\textrm{(3S)}} $ mesons (bottom) in jets with 60 $ < p_{\mathrm{T}} < $ 75 GeV.

png pdf
Figure 5-c:
Relative values of the systematic uncertainties discussed in the text, together with the total uncertainty band, for the longitudinal and transverse profiles of $ \Upsilon{\textrm{(1S)}} $ (top), $ \Upsilon{\textrm{(2S)}} $ (middle) and $ \Upsilon{\textrm{(3S)}} $ mesons (bottom) in jets with 60 $ < p_{\mathrm{T}} < $ 75 GeV.

png pdf
Figure 5-d:
Relative values of the systematic uncertainties discussed in the text, together with the total uncertainty band, for the longitudinal and transverse profiles of $ \Upsilon{\textrm{(1S)}} $ (top), $ \Upsilon{\textrm{(2S)}} $ (middle) and $ \Upsilon{\textrm{(3S)}} $ mesons (bottom) in jets with 60 $ < p_{\mathrm{T}} < $ 75 GeV.

png pdf
Figure 5-e:
Relative values of the systematic uncertainties discussed in the text, together with the total uncertainty band, for the longitudinal and transverse profiles of $ \Upsilon{\textrm{(1S)}} $ (top), $ \Upsilon{\textrm{(2S)}} $ (middle) and $ \Upsilon{\textrm{(3S)}} $ mesons (bottom) in jets with 60 $ < p_{\mathrm{T}} < $ 75 GeV.

png pdf
Figure 5-f:
Relative values of the systematic uncertainties discussed in the text, together with the total uncertainty band, for the longitudinal and transverse profiles of $ \Upsilon{\textrm{(1S)}} $ (top), $ \Upsilon{\textrm{(2S)}} $ (middle) and $ \Upsilon{\textrm{(3S)}} $ mesons (bottom) in jets with 60 $ < p_{\mathrm{T}} < $ 75 GeV.

png pdf
Figure 6:
Particle-level results for the longitudinal (left) and transverse profiles (right) of jets containing $ \Upsilon{\textrm{(1S)}} $ mesons in three different $ p_{\mathrm{T}} $ regions. The data (black points) are compared to the expectations of PYTHIA8.240 with the CP1 (blue) and CP5 tunes (red), as well as the PYTHIA8.310 predictions with the CP1 (orange) and CP5 tunes (violet). The vertical error bars represent the statistical uncertainties on the data points, while the grey bands represent the total experimental uncertainties.

png pdf
Figure 6-a:
Particle-level results for the longitudinal (left) and transverse profiles (right) of jets containing $ \Upsilon{\textrm{(1S)}} $ mesons in three different $ p_{\mathrm{T}} $ regions. The data (black points) are compared to the expectations of PYTHIA8.240 with the CP1 (blue) and CP5 tunes (red), as well as the PYTHIA8.310 predictions with the CP1 (orange) and CP5 tunes (violet). The vertical error bars represent the statistical uncertainties on the data points, while the grey bands represent the total experimental uncertainties.

png pdf
Figure 6-b:
Particle-level results for the longitudinal (left) and transverse profiles (right) of jets containing $ \Upsilon{\textrm{(1S)}} $ mesons in three different $ p_{\mathrm{T}} $ regions. The data (black points) are compared to the expectations of PYTHIA8.240 with the CP1 (blue) and CP5 tunes (red), as well as the PYTHIA8.310 predictions with the CP1 (orange) and CP5 tunes (violet). The vertical error bars represent the statistical uncertainties on the data points, while the grey bands represent the total experimental uncertainties.

png pdf
Figure 6-c:
Particle-level results for the longitudinal (left) and transverse profiles (right) of jets containing $ \Upsilon{\textrm{(1S)}} $ mesons in three different $ p_{\mathrm{T}} $ regions. The data (black points) are compared to the expectations of PYTHIA8.240 with the CP1 (blue) and CP5 tunes (red), as well as the PYTHIA8.310 predictions with the CP1 (orange) and CP5 tunes (violet). The vertical error bars represent the statistical uncertainties on the data points, while the grey bands represent the total experimental uncertainties.

png pdf
Figure 6-d:
Particle-level results for the longitudinal (left) and transverse profiles (right) of jets containing $ \Upsilon{\textrm{(1S)}} $ mesons in three different $ p_{\mathrm{T}} $ regions. The data (black points) are compared to the expectations of PYTHIA8.240 with the CP1 (blue) and CP5 tunes (red), as well as the PYTHIA8.310 predictions with the CP1 (orange) and CP5 tunes (violet). The vertical error bars represent the statistical uncertainties on the data points, while the grey bands represent the total experimental uncertainties.

png pdf
Figure 6-e:
Particle-level results for the longitudinal (left) and transverse profiles (right) of jets containing $ \Upsilon{\textrm{(1S)}} $ mesons in three different $ p_{\mathrm{T}} $ regions. The data (black points) are compared to the expectations of PYTHIA8.240 with the CP1 (blue) and CP5 tunes (red), as well as the PYTHIA8.310 predictions with the CP1 (orange) and CP5 tunes (violet). The vertical error bars represent the statistical uncertainties on the data points, while the grey bands represent the total experimental uncertainties.

png pdf
Figure 6-f:
Particle-level results for the longitudinal (left) and transverse profiles (right) of jets containing $ \Upsilon{\textrm{(1S)}} $ mesons in three different $ p_{\mathrm{T}} $ regions. The data (black points) are compared to the expectations of PYTHIA8.240 with the CP1 (blue) and CP5 tunes (red), as well as the PYTHIA8.310 predictions with the CP1 (orange) and CP5 tunes (violet). The vertical error bars represent the statistical uncertainties on the data points, while the grey bands represent the total experimental uncertainties.

png pdf
Figure 7:
Particle-level results for the longitudinal (left) and transverse profiles (right) of jets containing $ \Upsilon{\textrm{(2S)}} $ mesons in three different $ p_{\mathrm{T}} $ regions. The data (black points) are compared to the expectations of PYTHIA8.240 with the CP1 (blue) and CP5 tunes (red), as well as the PYTHIA8.310 predictions with the CP1 (orange) and CP5 tunes (violet). The vertical error bars represent the statistical uncertainties on the data points, while the grey bands represent the total experimental uncertainties.

png pdf
Figure 7-a:
Particle-level results for the longitudinal (left) and transverse profiles (right) of jets containing $ \Upsilon{\textrm{(2S)}} $ mesons in three different $ p_{\mathrm{T}} $ regions. The data (black points) are compared to the expectations of PYTHIA8.240 with the CP1 (blue) and CP5 tunes (red), as well as the PYTHIA8.310 predictions with the CP1 (orange) and CP5 tunes (violet). The vertical error bars represent the statistical uncertainties on the data points, while the grey bands represent the total experimental uncertainties.

png pdf
Figure 7-b:
Particle-level results for the longitudinal (left) and transverse profiles (right) of jets containing $ \Upsilon{\textrm{(2S)}} $ mesons in three different $ p_{\mathrm{T}} $ regions. The data (black points) are compared to the expectations of PYTHIA8.240 with the CP1 (blue) and CP5 tunes (red), as well as the PYTHIA8.310 predictions with the CP1 (orange) and CP5 tunes (violet). The vertical error bars represent the statistical uncertainties on the data points, while the grey bands represent the total experimental uncertainties.

png pdf
Figure 7-c:
Particle-level results for the longitudinal (left) and transverse profiles (right) of jets containing $ \Upsilon{\textrm{(2S)}} $ mesons in three different $ p_{\mathrm{T}} $ regions. The data (black points) are compared to the expectations of PYTHIA8.240 with the CP1 (blue) and CP5 tunes (red), as well as the PYTHIA8.310 predictions with the CP1 (orange) and CP5 tunes (violet). The vertical error bars represent the statistical uncertainties on the data points, while the grey bands represent the total experimental uncertainties.

png pdf
Figure 7-d:
Particle-level results for the longitudinal (left) and transverse profiles (right) of jets containing $ \Upsilon{\textrm{(2S)}} $ mesons in three different $ p_{\mathrm{T}} $ regions. The data (black points) are compared to the expectations of PYTHIA8.240 with the CP1 (blue) and CP5 tunes (red), as well as the PYTHIA8.310 predictions with the CP1 (orange) and CP5 tunes (violet). The vertical error bars represent the statistical uncertainties on the data points, while the grey bands represent the total experimental uncertainties.

png pdf
Figure 7-e:
Particle-level results for the longitudinal (left) and transverse profiles (right) of jets containing $ \Upsilon{\textrm{(2S)}} $ mesons in three different $ p_{\mathrm{T}} $ regions. The data (black points) are compared to the expectations of PYTHIA8.240 with the CP1 (blue) and CP5 tunes (red), as well as the PYTHIA8.310 predictions with the CP1 (orange) and CP5 tunes (violet). The vertical error bars represent the statistical uncertainties on the data points, while the grey bands represent the total experimental uncertainties.

png pdf
Figure 7-f:
Particle-level results for the longitudinal (left) and transverse profiles (right) of jets containing $ \Upsilon{\textrm{(2S)}} $ mesons in three different $ p_{\mathrm{T}} $ regions. The data (black points) are compared to the expectations of PYTHIA8.240 with the CP1 (blue) and CP5 tunes (red), as well as the PYTHIA8.310 predictions with the CP1 (orange) and CP5 tunes (violet). The vertical error bars represent the statistical uncertainties on the data points, while the grey bands represent the total experimental uncertainties.

png pdf
Figure 8:
Particle-level results for the longitudinal (left) and transverse profiles (right) of jets containing $ \Upsilon{\textrm{(3S)}} $ mesons in three different $ p_{\mathrm{T}} $ regions. The data (black points) are compared to the expectations of PYTHIA8.240 with the CP1 (blue) and CP5 tunes (red), as well as the PYTHIA8.310 predictions with the CP1 (orange) and CP5 tunes (violet). The vertical error bars represent the statistical uncertainties on the data points, while the grey bands represent the total experimental uncertainties.

png pdf
Figure 8-a:
Particle-level results for the longitudinal (left) and transverse profiles (right) of jets containing $ \Upsilon{\textrm{(3S)}} $ mesons in three different $ p_{\mathrm{T}} $ regions. The data (black points) are compared to the expectations of PYTHIA8.240 with the CP1 (blue) and CP5 tunes (red), as well as the PYTHIA8.310 predictions with the CP1 (orange) and CP5 tunes (violet). The vertical error bars represent the statistical uncertainties on the data points, while the grey bands represent the total experimental uncertainties.

png pdf
Figure 8-b:
Particle-level results for the longitudinal (left) and transverse profiles (right) of jets containing $ \Upsilon{\textrm{(3S)}} $ mesons in three different $ p_{\mathrm{T}} $ regions. The data (black points) are compared to the expectations of PYTHIA8.240 with the CP1 (blue) and CP5 tunes (red), as well as the PYTHIA8.310 predictions with the CP1 (orange) and CP5 tunes (violet). The vertical error bars represent the statistical uncertainties on the data points, while the grey bands represent the total experimental uncertainties.

png pdf
Figure 8-c:
Particle-level results for the longitudinal (left) and transverse profiles (right) of jets containing $ \Upsilon{\textrm{(3S)}} $ mesons in three different $ p_{\mathrm{T}} $ regions. The data (black points) are compared to the expectations of PYTHIA8.240 with the CP1 (blue) and CP5 tunes (red), as well as the PYTHIA8.310 predictions with the CP1 (orange) and CP5 tunes (violet). The vertical error bars represent the statistical uncertainties on the data points, while the grey bands represent the total experimental uncertainties.

png pdf
Figure 8-d:
Particle-level results for the longitudinal (left) and transverse profiles (right) of jets containing $ \Upsilon{\textrm{(3S)}} $ mesons in three different $ p_{\mathrm{T}} $ regions. The data (black points) are compared to the expectations of PYTHIA8.240 with the CP1 (blue) and CP5 tunes (red), as well as the PYTHIA8.310 predictions with the CP1 (orange) and CP5 tunes (violet). The vertical error bars represent the statistical uncertainties on the data points, while the grey bands represent the total experimental uncertainties.

png pdf
Figure 8-e:
Particle-level results for the longitudinal (left) and transverse profiles (right) of jets containing $ \Upsilon{\textrm{(3S)}} $ mesons in three different $ p_{\mathrm{T}} $ regions. The data (black points) are compared to the expectations of PYTHIA8.240 with the CP1 (blue) and CP5 tunes (red), as well as the PYTHIA8.310 predictions with the CP1 (orange) and CP5 tunes (violet). The vertical error bars represent the statistical uncertainties on the data points, while the grey bands represent the total experimental uncertainties.

png pdf
Figure 8-f:
Particle-level results for the longitudinal (left) and transverse profiles (right) of jets containing $ \Upsilon{\textrm{(3S)}} $ mesons in three different $ p_{\mathrm{T}} $ regions. The data (black points) are compared to the expectations of PYTHIA8.240 with the CP1 (blue) and CP5 tunes (red), as well as the PYTHIA8.310 predictions with the CP1 (orange) and CP5 tunes (violet). The vertical error bars represent the statistical uncertainties on the data points, while the grey bands represent the total experimental uncertainties.
Summary
A measurement of the fragmentation properties of jets containing $ \Upsilon{\textrm{(nS)}} $ quarkonia has been presented for the first time. The measurement is performed using the dimuon decays of the $ \Upsilon $ mesons and the jets containing the dimuon system. The measurements include the longitudinal profile $ z $, the distribution of the fraction of the jet momentum carried by the $ \Upsilon $ meson in the direction of the jet axis, and the transverse profile $ p_{\mathrm{T}}^{\text{rel}} $, the distribution of the momentum of the $ \Upsilon $ meson in the direction perpendicular to the jet axis. The results are compared with the predictions from PYTHIA8.240, in which $ \Upsilon $ mesons are generated directly from the matrix element, as well as with the PYTHIA8.310 predictions, using an alternative approach for the generation of quarkonia from the parton shower. Both PYTHIA versions make use of the CP1 and CP5 tunes. These measurements provide new and stringent tests of non-relativistic QCD (NRQCD). Compared with charmonium, bottomonium production is expected to be more reliably described within the NRQCD framework, because of the larger mass of the b quark, which improves the separation between perturbative and non-perturbative scales. Moreover, the inclusion of the transverse momentum profile of the $ \Upsilon $ meson relative to the jet axis further extends the sensitivity of these measurements to transverse-momentum-dependent aspects of quarkonium production, providing additional constraints on the dynamics of quarkonium formation beyond collinear fragmentation approaches. The results show important discrepancies between the data and the predictions, pointing to $ \Upsilon $ mesons carrying a smaller momentum fraction of the jet than what is predicted in the PYTHIA implementations of NRQCD. These results are consistent in magnitude with previous measurements of jets containing $ \mathrm{J}/\psi $ mesons by the LHCb and CMS Collaborations. While the PYTHIA8.310 improvements tend to go in the correct direction, showing less isolated quarkonia, these results suggest that further theoretical progress is needed to achieve a satisfactory description. Overall, the present measurement is expected to provide additional input to the so-called quarkonium polarization puzzle by constraining the long-distance matrix elements relevant to $ \Upsilon{\textrm{(nS)}} $ production in NRQCD. In particular, the simultaneous sensitivity to the longitudinal and transverse structure of $ \Upsilon{\textrm{(nS)}} $ production in jets may help to disentangle the interplay between colour-singlet and colour-octet production mechanisms in this context, as well as the contributions from higher angular momentum Fock states.
References
1 W. E. Caswell and G. P. Lepage Effective Lagrangians for bound state problems in QED, QCD, and other field theories PLB 167 (1986) 437
2 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
3 P. Cho and A. K. Leibovich Color-octet quarkonia production PRD 53 (1996) 150 hep-ph/9505329
4 P. Cho and A. K. Leibovich Color-octet quarkonia production. II PRD 53 (1996) 6203 hep-ph/9511315
5 CDF Collaboration $ \mathrm{J}/\psi $ and \PGyP2S production in $ {\mathrm{p}\overline{\mathrm{p}}} $ collisions at $ \sqrt{s}= $ 1.8 TeV PRL 79 (1997) 572
6 CDF Collaboration $ \Upsilon $ production in $ {\mathrm{p}\overline{\mathrm{p}}} $ collisions at $ \sqrt{s}= $ 1.8 TeV PRL 75 (1995) 4358
7 \DZERO Collaboration $ \mathrm{J}/\psi $ production in $ {\mathrm{p}\overline{\mathrm{p}}} $ collisions at $ \sqrt{s}= $ 1.8 TeV PLB 370 (1996) 239
8 \DZERO Collaboration Measurement of inclusive differential cross sections for \PGUP1S production in $ {\mathrm{p}\overline{\mathrm{p}}} $ collisions at $ \sqrt{s}= $ 1.96 TeV PRL 94 (2005) 232001 hep-ex/0502030
9 CMS Collaboration Measurement of quarkonium production cross sections in $ {\mathrm{p}\mathrm{p}} $ collisions at $ \sqrt{s}= $ 13 TeV PLB 780 (2018) 251 CMS-BPH-15-005
1710.11002
10 ATLAS Collaboration Measurement of the production cross-section of $ \mathrm{J}/\psi $ and \PGyP2S mesons in $ {\mathrm{p}\mathrm{p}} $ collisions at $ \sqrt{s}= $ 13 TeV with the ATLAS detector EPJC 84 (2024) 169 2309.17177
11 LHCb Collaboration Measurement of forward $ \mathrm{J}/\psi $ production cross-sections in $ {\mathrm{p}\mathrm{p}} $ collisions at $ \sqrt{s}= $ 13 TeV JHEP 10 (2015) 172 1509.00771
12 ATLAS Collaboration Measurement of upsilon production in 7 TeV $ {\mathrm{p}\mathrm{p}} $ collisions at ATLAS PRD 87 (2013) 052004 1211.7255
13 CMS Collaboration Measurement of the \PGUP1S, \PGUP2S, and \PGUP3S differential cross sections in $ {\mathrm{p}\mathrm{p}} $ collisions at $ \sqrt{s}= $ 13.6 TeV Submitted to JHEP, 2026 CMS-BPH-24-004
2601.20023
14 LHCb Collaboration Measurement of $ \Upsilon $ production in $ {\mathrm{p}\mathrm{p}} $ collisions at $ \sqrt{s}= $ 13 TeV JHEP 07 (2018) 134 1804.09214
15 E. Braaten, B. A. Kniehl, and J. Lee Polarization of prompt $ \mathrm{J}/\psi $ at the Fermilab Tevatron PRD 62 (2000) 094005 hep-ph/9911436
16 M. Beneke and M. Kr ä mer Direct $ \mathrm{J}/\psi $ and $ {\psi^\prime} $ polarization and cross sections at the Fermilab Tevatron PRD 55 (1997) 5269 hep-ph/9611218
17 CDF Collaboration Polarizations of $ \mathrm{J}/\psi $ and \PGyP2S mesons produced in $ {\mathrm{p}\overline{\mathrm{p}}} $ collisions at $ \sqrt{s}= $ 1.96 TeV PRL 99 (2007) 132001 0704.0638
18 \DZERO Collaboration Measurement of the polarization of the \PGUP1S and \PGUP2S states in $ {\mathrm{p}\overline{\mathrm{p}}} $ collisions at $ \sqrt{s}= $ 1.96 TeV PRL 101 (2008) 182004 0804.2799
19 CDF Collaboration Measurements of the angular distributions of muons from $ \Upsilon $ decays in $ {\mathrm{p}\overline{\mathrm{p}}} $ collisions at $ \sqrt{s}= $ 1.96 TeV PRL 108 (2012) 151802 1112.1591
20 CMS Collaboration Measurement of the prompt $ \mathrm{J}/\psi $ and \PGyP2S polarizations in $ {\mathrm{p}\mathrm{p}} $ collisions at $ \sqrt{s}= $ 7 TeV PLB 727 (2013) 381 CMS-BPH-13-003
1307.6070
21 CMS Collaboration Measurement of the polarizations of prompt and non-prompt $ \mathrm{J}/\psi $ and \PGyP2s mesons produced in $ {\mathrm{p}\mathrm{p}} $ collisions at $ \sqrt{s}= $ 13 TeV PLB 858 (2024) 139044 CMS-BPH-22-009
2406.14409
22 CMS Collaboration Measurement of the \PGUP1S, \PGUP2S and \PGUP3S polarizations in $ {\mathrm{p}\mathrm{p}} $ collisions at $ \sqrt{s}= $ 7 TeV PRL 110 (2013) 081802 CMS-BPH-11-023
1209.2922
23 LHCb Collaboration Measurement of $ \mathrm{J}/\psi $ polarization in $ {\mathrm{p}\mathrm{p}} $ collisions at $ \sqrt{s}= $ 7 TeV EPJC 73 (2013) 2631 1307.6379
24 LHCb Collaboration Measurement of the $ \Upsilon{\textrm{(nS)}} $ polarizations in $ {\mathrm{p}\mathrm{p}} $ collisions at $ \sqrt{s}= $ 7 and 8 TeV JHEP 12 (2017) 110 1709.01301
25 B. Mele and P. Nason The fragmentation function for heavy quarks in QCD [Corrigendum: \DOI10./j.nuclphysb..05.005], 1991
NPB 361 (1991) 626
26 M. Cacciari and E. Gardi Heavy-quark fragmentation NPB 664 (2003) 299 hep-ph/0301047
27 ATLAS Collaboration Measurement of $ {\mathrm{D}^{\ast}(2010)^{\pm}} $ meson production in jets from $ {\mathrm{p}\mathrm{p}} $ collisions at $ \sqrt{s}= $ 7 TeV with the ATLAS detector PRD 85 (2012) 052005 1112.4432
28 ATLAS Collaboration Measurement of b-quark fragmentation properties in jets using the decay $ {{\mathrm{B}^{\pm}}\to\mathrm{J}/\psi \mathrm{K^{\pm}}} $ in $ {\mathrm{p}\mathrm{p}} $ collisions at $ \sqrt{s}= $ 13 TeV with the ATLAS detector JHEP 12 (2021) 131 2108.11650
29 ALICE Collaboration Measurement of the production of charm jets tagged with $ \mathrm{D^0} $ mesons in $ {\mathrm{p}\mathrm{p}} $ collisions at $ \sqrt{s}= $ 5.02 and 13 TeV JHEP 06 (2023) 133 2204.10167
30 LHCb Collaboration $ {\mathrm{B}} $-jet fragmentation with $ {{\mathrm{B}^{\pm}}\to\mathrm{J}/\psi \mathrm{K^{\pm}}} $ decays in $ \sqrt{s}= $ 13 TeV $ {\mathrm{p}\mathrm{p}} $ collisions at LHCb PRD 113 (2026) 112006 2603.23663
31 M. Baumgart, A. K. Leibovich, T. Mehen, and I. Z. Rothstein Probing quarkonium production mechanisms with jet substructure JHEP 11 (2014) 003 1406.2295
32 Z.-B. Kang et al. $ \mathrm{J}/\psi $ production and polarization within a jet PRL 119 (2017) 032001 1702.03287
33 LHCb Collaboration Study of $ \mathrm{J}/\psi $ production in jets PRL 118 (2017) 192001 1701.05116
34 LHCb Collaboration Measurements of \PGyP2S and \PGccDoP3872 production within fully reconstructed jets EPJC 85 (2025) 562 2410.18018
35 CMS Collaboration Study of $ \mathrm{J}/\psi $ meson production inside jets in $ {\mathrm{p}\mathrm{p}} $ collisions at $ \sqrt{s}= $ 8 TeV PLB 804 (2020) 135409 CMS-BPH-15-003
1910.01686
36 CMS Collaboration Fragmentation of jets containing a prompt $ \mathrm{J}/\psi $ meson in PbPb and $ {\mathrm{p}\mathrm{p}} $ collisions at $ {\sqrt{\smash[b]{s_{_{\mathrm{NN}}}}}= $ 5.02 TeV PLB 825 (2021) 136842 CMS-HIN-19-007
2106.13235
37 T. Sjöstrand et al. An introduction to PYTHIA8.2 Comput. Phys. Commun. 191 (2015) 159 1410.3012
38 C. Bierlich et al. A comprehensive guide to the physics and usage of PYTHIA8.3 SciPost Phys. Codeb. 8 (2022) 2203.11601
39 R. Bain et al. NRQCD confronts LHCb data on quarkonium production within jets PRL 119 (2017) 032002 1702.05525
40 Y. Wang, D. Kang, and H. S. Chung NRQCD reconfronts LHCb data on quarkonium production within jets PRD 112 (2025) L071504 2507.19022
41 N. Cooke, P. Ilten, L. Lönnblad, and S. Mrenna Non-relativistic quantum chromodynamics in parton showers EPJC 84 (2024) 432 2312.05203
42 F. G. Celiberto On the quarkonium-in-jet collinear fragmentation at moderate-to-large transverse momentum in Proc. Diffraction and Low-$ x $: Palermo, Italy, September 08--14,.. [-A22], 2025
Acta Phys. Polon. Supp. 18 (2025) 1
2412.05661
43 F. G. Celiberto and F. Lonigro Pseudoscalar heavy-quarkonium hadroproduction from nonrelativistic fragmentation at NLL/NLO\textsuperscript+ PRD 112 (2025) 114040 2510.10593
44 M. G. Echevarria, Y. Makris, and I. Scimemi Quarkonium TMD fragmentation functions in NRQCD JHEP 10 (2020) 164 2007.05547
45 M. G. Echevarria, S. F. Romera, and I. Scimemi Gluon TMD fragmentation function into quarkonium JHEP 12 (2023) 181 2308.12356
46 M. Copeland et al. Polarized TMD fragmentation functions for $ \mathrm{J}/\psi $ production PRD 109 (2024) 054017 2308.08605
47 E. Chapon et al. Prospects for quarkonium studies at the high-luminosity LHC Prog. Part. Nucl. Phys. 122 (2022) 103906 2012.14161
48 CMS Collaboration HEPData record for this analysis link
49 CMS Collaboration The CMS experiment at the CERN LHC JINST 3 (2008) S08004
50 CMS Collaboration Development of the CMS detector for the CERN LHC \mboxRun 3 JINST 19 (2024) P05064 CMS-PRF-21-001
2309.05466
51 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
52 CMS Collaboration The CMS trigger system JINST 12 (2017) P01020 CMS-TRG-12-001
1609.02366
53 CMS Collaboration Performance of the CMS high-level trigger during LHC \mboxRun 2 JINST 19 (2024) P11021 CMS-TRG-19-001
2410.17038
54 CMS Collaboration Electron and photon reconstruction and identification with the CMS experiment at the CERN LHC JINST 16 (2021) P05014 CMS-EGM-17-001
2012.06888
55 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
56 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
57 CMS Collaboration Particle-flow reconstruction and global event description with the CMS detector JINST 12 (2017) P10003 CMS-PRF-14-001
1706.04965
58 M. Cacciari, G. P. Salam, and G. Soyez The anti-$ k_{\mathrm{T}} $ jet clustering algorithm JHEP 04 (2008) 063 0802.1189
59 M. Cacciari, G. P. Salam, and G. Soyez FASTJET user manual EPJC 72 (2012) 1896 1111.6097
60 CMS Collaboration Pileup mitigation at CMS in 13 TeV data JINST 15 (2020) P09018 CMS-JME-18-001
2003.00503
61 CMS Collaboration Jet energy scale and resolution in the CMS experiment in $ {\mathrm{p}\mathrm{p}} $ collisions at 8 TeV JINST 12 (2017) P02014 CMS-JME-13-004
1607.03663
62 CMS Collaboration Extraction and validation of a new set of CMS PYTHIA8 tunes from underlying-event measurements EPJC 80 (2020) 4 CMS-GEN-17-001
1903.12179
63 NNPDF Collaboration Parton distributions from high-precision collider data EPJC 77 (2017) 663 1706.00428
64 R. Baier and R. R \"u ckl Hadronic collisions: A quarkonium factory Z. Phys. C 19 (1983) 251
65 R. Gastmans, W. Troost, and T. T. Wu Production of heavy quarkonia from gluons NPB 291 (1987) 731
66 D. J. Lange The EVTGEN particle decay simulation package NIM A 462 (2001) 152
67 GEANT4 Collaboration GEANT 4---a simulation toolkit NIM A 506 (2003) 250
68 A. Bodek et al. Extracting muon momentum scale corrections for hadron collider experiments EPJC 72 (2012) 2194 1208.3710
69 CMS Collaboration Measurement of the inelastic proton-proton cross section at $ \sqrt{s}= $ 13 TeV JHEP 07 (2018) 161 CMS-FSQ-15-005
1802.02613
70 CMS Collaboration Performance of the CMS muon trigger system in proton-proton collisions at $ \sqrt{s}= $ 13 TeV JINST 16 (2021) P07001 CMS-MUO-19-001
2102.04790
71 E. Braaten and T. C. Yuan Gluon fragmentation into heavy quarkonium PRL 71 (1993) 1673 hep-ph/9303205
72 E. Braaten, K. Cheung, and T. C. Yuan $ \mathrm{Z^0} $ decay into charmonium via charm quark fragmentation PRD 48 (1993) 4230 hep-ph/9302307
73 E. Braaten and T. C. Yuan Gluon fragmentation into $ {P} $-wave heavy quarkonium PRD 50 (1994) 3176 hep-ph/9403401
74 E. Braaten and T. C. Yuan Gluon fragmentation into spin-triplet $ {S} $-wave quarkonium PRD 52 (1995) 6627 hep-ph/9507398
75 T. C. Yuan Perturbative QCD fragmentation functions for production of $ {P} $-wave charm and bottom mesons PRD 50 (1994) 5664 hep-ph/9405348
76 CMS Collaboration Precision luminosity measurement in proton-proton collisions at $ \sqrt{s}= $ 13 TeV in 2015 and 2016 at CMS EPJC 81 (2021) 800 CMS-LUM-17-003
2104.01927
77 CMS Collaboration CMS luminosity measurement for the 2017 data-taking period at $ \sqrt{s}= $ 13 TeV CMS Physics Analysis Summary, 2018
CMS-PAS-LUM-17-004
CMS-PAS-LUM-17-004
78 CMS Collaboration CMS luminosity measurement for the 2018 data-taking period at $ \sqrt{s}= $ 13 TeV CMS Physics Analysis Summary, 2019
CMS-PAS-LUM-18-002
CMS-PAS-LUM-18-002
79 M. J. Oreglia A study of the reactions $ {\psi^\prime\to\gamma\gamma\psi} $ PhD thesis, Stanford University, SLAC-R-236, 1980
link
80 J. E. Gaiser Charmonium spectroscopy from radiative decays of the $ \mathrm{J}/\psi $ and $ {\psi^\prime} $ PhD thesis, Stanford University, SLAC-0255, 1982
link
81 Particle Data Group , S. Navas et al. Review of particle physics PRD 110 (2024) 030001
82 G. D'Agostini A multidimensional unfolding method based on Bayes' theorem NIM A 362 (1995) 487
83 T. Adye Unfolding algorithms and tests using RooUnfold in orkshop on Statistical Issues Related to Discovery Claims in Search Experiments and Unfolding (PHYSTAT ): Geneva, Switzerland, January 17--20,, 2011
Proc. 2011 (2011) W
1105.1160
84 CMS Collaboration Measurement of the inclusive W and Z production cross sections in $ {\mathrm{p}\mathrm{p}} $ collisions at $ \sqrt{s}= $ 7 TeV with the CMS experiment JHEP 10 (2011) 132 CMS-EWK-10-005
1107.4789
Compact Muon Solenoid
LHC, CERN