pith. machine review for the scientific record. sign in

arxiv: 2604.07968 · v1 · submitted 2026-04-09 · ✦ hep-ex · nucl-ex

Recognition: 2 theorem links

· Lean Theorem

Multiplicity dependence of prompt and non-prompt J/psi production at midrapidity in pp collisions at sqrt{s} = 13 TeV

Authors on Pith no claims yet

Pith reviewed 2026-05-10 18:03 UTC · model grok-4.3

classification ✦ hep-ex nucl-ex
keywords J/ψ productionmultiplicity dependenceprompt non-prompt separationpp collisionsdielectron decayazimuthal regionsself-normalized yieldsheavy quarkonium
0
0 comments X

The pith

Prompt and non-prompt J/ψ production in proton-proton collisions rises stronger than linearly with charged-particle multiplicity.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The paper measures the yields of J/ψ mesons produced directly in the collision and those originating from beauty-hadron decays, both normalized to their average values in inelastic events. These self-normalized yields are examined as a function of the number of charged particles produced at midrapidity in 13 TeV pp collisions, using the dielectron decay channel. Both prompt and non-prompt components exhibit a stronger-than-linear increase with multiplicity, and the trends remain similar between the two. When the multiplicity is divided into azimuthal sectors relative to the J/ψ direction, the strong rise persists in the toward region but weakens in the transverse and away regions. The ratio of prompt J/ψ to D0 yields is also reported as varying with multiplicity.

Core claim

Self-normalized yields of both prompt and non-prompt J/ψ increase more steeply than linearly with rising midrapidity charged-particle multiplicity in pp collisions at √s = 13 TeV. The increase is comparable for the two components and is strongest when multiplicity is measured in the azimuthal toward region aligned with the J/ψ, while it is weaker in the transverse and away regions. The multiplicity dependence of the prompt J/ψ to D0 yield ratio is additionally quantified.

What carries the argument

Self-normalized yields of J/ψ mesons separated into prompt and non-prompt components via the dielectron decay channel, evaluated in bins of midrapidity charged-particle multiplicity and further subdivided into three azimuthal regions (toward, transverse, away) relative to the J/ψ momentum direction.

If this is right

  • Both direct and beauty-decay J/ψ production mechanisms are enhanced in high-multiplicity pp events beyond simple linear scaling with particle density.
  • The similar multiplicity dependence for prompt and non-prompt J/ψ points to a common underlying driver affecting heavy-quark production.
  • The azimuthal variation, with stronger rise toward the J/ψ direction, indicates that the enhancement correlates with the direction of the leading particle.
  • The reported multiplicity dependence of the prompt J/ψ to D0 ratio constrains models of charm and beauty production in small systems.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

  • The observed trends in pp collisions may share a common origin with collective-like behavior reported in high-multiplicity small systems.
  • If the enhancement is universal across quarkonia, it could affect the baseline used to interpret suppression signals in heavy-ion collisions.
  • Measurements with other heavy-flavor species or at different collision energies would test whether the stronger-than-linear scaling is a general feature of high-multiplicity pp events.

Load-bearing premise

Prompt and non-prompt J/ψ components can be accurately separated in the dielectron channel without multiplicity-dependent biases from efficiency corrections or event selection.

What would settle it

A new measurement that finds only linear or weaker growth in self-normalized J/ψ yields with multiplicity, or that obtains inconsistent prompt versus non-prompt trends when using an independent separation technique such as decay-vertex displacement.

Figures

Figures reproduced from arXiv: 2604.07968 by ALICE Collaboration.

Figure 1
Figure 1. Figure 1: Regions in azimuthal angle defined with respect to a J/ψ candidate. In events containing a J/ψ candidate, the multiplicity can be separated into several azimuthal regions relative to the J/ψ emission direction. These azimuthal regions are as labelled in [PITH_FULL_IMAGE:figures/full_fig_p005_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: Raw inclusive J/ψ yield as a function of the selection on the BDT non-prompt output probability for 39 ≤ Ntrks < 46, and pT > 1 GeV/c. The error bars indicate statistical uncertainties. The lines correspond to the prompt (dashed), non-prompt (dot-dashed) and total (continuous) J/ψ scaled efficiencies. See text for more details. estimated from MC, either in the specific Ntrks interval (for A i and ε i ) or … view at source ↗
Figure 3
Figure 3. Figure 3: Fraction of non-prompt J/ψ mesons, fB, for pT > 1 GeV/c as a function of the self-normalized charged￾particle multiplicity measured within |η| < 0.9 for the INEL>0 (full markers) and V0M 0–0.1% (open markers) event classes. The vertical bars and boxes indicate the statistical and systematic uncertainties, respectively. The results are compared to PYTHIA 8 simulations using different tunes, as indicated by … view at source ↗
Figure 4
Figure 4. Figure 4: Self-normalized prompt (left) and non-prompt (right) J/ψ yields with pT > 1 GeV/c as a function of self-normalized charged-particle multiplicity at midrapidity. The data are shown separately for the INEL>0 and V0M 0–0.1% event classes. The yields are compared to PYTHIA 8 using different settings [69, 88], as well as EPOS4 [48–51] and the 3-Pomeron CGC model [40, 41]. Lower panels show the ratios to the exp… view at source ↗
Figure 5
Figure 5. Figure 5: Prompt J/ψ yield for 1 < pT < 4 GeV/c (left), 4 < pT < 8 GeV/c (middle) and pT > 8 GeV/c (right) as a function of self-normalized charged-particle multiplicity within |η| < 0.9 for the INEL>0 and V0M 0–0.1% event classes. The data are compared to PYTHIA 8 using different settings [69, 88], as well as EPOS4HQ [48–51], the 3-Pomeron CGC model [40, 41] and the CGC+ICEM model [38]. Lower panels show the ratios… view at source ↗
Figure 6
Figure 6. Figure 6: Non-prompt J/ψ yield for 1 < pT < 4 GeV/c (left), 4 < pT < 8 GeV/c (middle) and pT > 8 GeV/c (right) as a function of self-normalized charged-particle multiplicity within |η| < 0.9 for the INEL>0 and V0M 0–0.1% event classes. The data are compared to PYTHIA 8 using different settings [69, 88], as well as EPOS4HQ with and without hydrodynamics activated [48, 92]. Lower panels show the ratios to the expectat… view at source ↗
Figure 7
Figure 7. Figure 7: The fractions of J/ψ mesons from beauty feed-down with pT > 1 GeV/c as a function of the self￾normalized charged-particle multiplicity in |η| < 0.9 in the toward (top-left), transverse (top-right), and away (bottom left) azimuthal regions with respect to the J/ψ momentum direction. The data are shown separately for the INEL>0 and V0M 0–0.1% event classes and are compared to results from PYTHIA 8 simulation… view at source ↗
Figure 8
Figure 8. Figure 8: Self-normalized prompt J/ψ yields with pT > 1 GeV/c as a function of the self-normalized charged￾particle multiplicity in |η| < 0.9 measured in the toward (left), transverse (middle), and away (right) azimuthal regions with respect to the J/ψ momentum direction. The data are shown separately for the INEL>0 and V0M 0– 0.1% event classes and are compared to PYTHIA 8 simulations using different settings [69, … view at source ↗
Figure 9
Figure 9. Figure 9: Self-normalized non-prompt J/ψ yields with pT > 1 GeV/c as a function of the self-normalized charged-particle multiplicity in |η| < 0.9 measured in the toward (left), transverse (middle), and away (right) azimuthal regions with respect to the J/ψ momentum direction. The data are shown separately for the INEL>0 and V0M 0–0.1% event classes and are compared to PYTHIA 8 simulations using different settings [6… view at source ↗
Figure 10
Figure 10. Figure 10: Self-normalized prompt J/ψ yields in the pT intervals 1 < pT < 4 GeV/c (left column), 4 < pT < 8 GeV/c (middle column), and pT > 8 GeV/c (right column) as a function of the self-normalized charged-particle multiplicity in |η| < 0.9 measured in toward (top row), transverse (middle row), and away (bottom row) az￾imuthal regions with respect to the J/ψ momentum direction. The data are shown separately for th… view at source ↗
Figure 11
Figure 11. Figure 11: Self-normalized non-prompt J/ψ yields in the pT intervals 1 < pT < 4 GeV/c (left column), 4 < pT < 8 GeV/c (middle column), and pT > 8 GeV/c (right column) as a function of the self-normalized charged￾particle multiplicity in |η| < 0.9 measured in toward (top row), transverse (middle row), and away (bottom row) azimuthal regions with respect to the J/ψ momentum direction. The data are shown separately for… view at source ↗
Figure 12
Figure 12. Figure 12: Ratio between the pT-integrated prompt J/ψ and prompt D0 yields as a function of charged-particle multiplicity within |η| < 0.5 in several collision systems [52, 96, 97]. This is compared to the prediction from the Statistical Hadronization Model (SHMc) [98]. The vertical bars, empty and shaded boxes indicate the statistical, systematic and extrapolation uncertainties, respectively. 4.3 J/ψ-to-D0 ratio Th… view at source ↗
read the original abstract

The yields of prompt and non-prompt J/$\psi$ and the fraction of non-prompt J/$\psi$ are measured at midrapidity ($|y| < 0.9$) via the dielectron decay channel as a function of the midrapidity charged-particle multiplicity ($|\eta| < 0.9$) in pp collisions at $\sqrt{s} = 13$ TeV. The J/$\psi$ yields and the multiplicity are normalized by their average value in inelastic collisions. The multiplicity-dependent yield ratio between prompt J/$\psi$ and D$^0$ is reported. The multiplicity is further divided into three azimuthal regions with respect to the J/$\psi$ momentum: toward the J/$\psi$ emission direction, transverse, or opposite to it. A stronger-than-linear increase of the self-normalized yields is observed for both prompt and non-prompt J/$\psi$ production, with similar trends. This behaviour is also observed in the toward region, while a weaker increase is observed in the transverse and away regions.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

2 major / 2 minor

Summary. The manuscript reports measurements of prompt and non-prompt J/ψ yields at midrapidity (|y| < 0.9) in pp collisions at √s = 13 TeV via the dielectron channel, as a function of midrapidity charged-particle multiplicity (|η| < 0.9). Both yields and multiplicity are self-normalized to their inelastic averages. Prompt/non-prompt separation uses template fits to the pseudo-proper decay length. The multiplicity is further subdivided into toward, transverse, and away azimuthal regions relative to the J/ψ direction. The central result is a stronger-than-linear rise in self-normalized yields for both prompt and non-prompt J/ψ (with similar trends), more pronounced in the toward region than in the transverse and away regions; the prompt J/ψ to D⁰ yield ratio versus multiplicity is also presented.

Significance. If the trends are robust, the work supplies new experimental input on charmonium production mechanisms in high-multiplicity pp events. The simultaneous treatment of prompt and non-prompt components together with azimuthal-region dependence helps discriminate between multiple-parton-interaction effects and possible collective phenomena in small systems. Self-normalization reduces certain systematic uncertainties and facilitates direct comparison with models and other experiments.

major comments (2)
  1. [Analysis procedure for prompt/non-prompt separation] The prompt/non-prompt separation via pseudo-proper decay length template fits is load-bearing for the reported trends. The manuscript should explicitly demonstrate that the fit quality, template shapes, and any multiplicity-dependent biases remain under control in the highest multiplicity bins (analysis section describing the separation procedure).
  2. [Yield extraction and efficiency corrections] Self-normalized yields rely on efficiency corrections and event selection being multiplicity-independent after normalization. The paper must quantify any residual multiplicity-dependent efficiency bias (e.g., via data-driven closure tests) because this directly affects the claimed stronger-than-linear behavior (results section on yield extraction).
minor comments (2)
  1. Figure captions and legends should explicitly label which curves correspond to prompt versus non-prompt J/ψ and to the three azimuthal regions for immediate readability.
  2. The abstract states that the prompt J/ψ–D⁰ ratio is reported; the main text should clarify the D⁰ selection criteria and whether the same multiplicity estimator is used for both species.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for the positive assessment of our manuscript and the recommendation for minor revision. We address the major comments point by point below and will revise the manuscript to incorporate the requested demonstrations and quantifications.

read point-by-point responses
  1. Referee: [Analysis procedure for prompt/non-prompt separation] The prompt/non-prompt separation via pseudo-proper decay length template fits is load-bearing for the reported trends. The manuscript should explicitly demonstrate that the fit quality, template shapes, and any multiplicity-dependent biases remain under control in the highest multiplicity bins (analysis section describing the separation procedure).

    Authors: We agree that explicit validation of the separation procedure is essential. In the revised manuscript we will expand the relevant analysis section with additional figures showing fit quality (including residuals and goodness-of-fit metrics) specifically for the highest multiplicity bins. We will also present comparisons of template shapes across multiplicity classes to demonstrate their stability and summarize Monte Carlo studies that bound any multiplicity-dependent biases to well below the reported uncertainties. revision: yes

  2. Referee: [Yield extraction and efficiency corrections] Self-normalized yields rely on efficiency corrections and event selection being multiplicity-independent after normalization. The paper must quantify any residual multiplicity-dependent efficiency bias (e.g., via data-driven closure tests) because this directly affects the claimed stronger-than-linear behavior (results section on yield extraction).

    Authors: We acknowledge the importance of quantifying residual biases. The revised manuscript will include, in the yield-extraction section, data-driven closure tests (based on embedded simulations) that explicitly evaluate multiplicity-dependent efficiency variations after self-normalization. These tests will provide quantitative upper limits on any residual bias, confirming that it does not affect the observed stronger-than-linear rise. revision: yes

Circularity Check

0 steps flagged

No significant circularity: direct experimental measurement

full rationale

The paper reports empirical measurements of self-normalized prompt and non-prompt J/ψ yields versus charged-particle multiplicity in pp collisions at 13 TeV, obtained via dielectron reconstruction, pseudo-proper decay length template fits for component separation, and normalization to inelastic averages. No theoretical derivation chain exists; the reported stronger-than-linear trends and azimuthal-region dependence are direct data outputs. Efficiency corrections, event selection, and multiplicity binning follow standard ALICE procedures without any fitted parameter being relabeled as a prediction or any result reducing to its own inputs by construction. Self-citations are limited to prior experimental methods and do not bear the central claim.

Axiom & Free-Parameter Ledger

0 free parameters · 0 axioms · 0 invented entities

The paper reports experimental measurements without introducing new theoretical parameters, axioms, or entities; it relies on standard particle physics analysis techniques and detector-specific corrections.

pith-pipeline@v0.9.0 · 5495 in / 1287 out tokens · 55488 ms · 2026-05-10T18:03:28.842413+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Lean theorems connected to this paper

Citations machine-checked in the Pith Canon. Every link opens the source theorem in the public Lean library.

What do these tags mean?
matches
The paper's claim is directly supported by a theorem in the formal canon.
supports
The theorem supports part of the paper's argument, but the paper may add assumptions or extra steps.
extends
The paper goes beyond the formal theorem; the theorem is a base layer rather than the whole result.
uses
The paper appears to rely on the theorem as machinery.
contradicts
The paper's claim conflicts with a theorem or certificate in the canon.
unclear
Pith found a possible connection, but the passage is too broad, indirect, or ambiguous to say the theorem truly supports the claim.

Reference graph

Works this paper leans on

44 extracted references · 42 canonical work pages · 1 internal anchor

  1. [1]

    Quarkonium Production in an Improved Color Evaporation Model

    Y .-Q. Ma and R. V ogt, “Quarkonium Production in an Improved Color Evaporation Model”, Phys. Rev. D94(2016) 114029,arXiv:1609.06042 [hep-ph]

  2. [2]

    Rigorous QCD Analysis of Inclusive Annihilation and Production of Heavy Quarkonium

    G. T. Bodwin, E. Braaten, and G. Lepage, “Rigorous QCD analysis of inclusive annihilation and production of heavy quarkonium”,Phys. Rev. D51(1995) 1125–1171, arXiv:hep-ph/9407339. [Erratum:Phys. Rev. D55(1997) 5853 ]. 23 Multiplicity-dependent prompt and non-prompt J/ψproduction at √s=13 TeV ALICE Collaboration

  3. [3]

    Comprehensive description of J/ψproduction in proton-proton collisions at collider energies

    Y .-Q. Ma and R. Venugopalan, “Comprehensive description of J/ψproduction in proton-proton collisions at collider energies”,Phys. Rev. Lett.113(2014) 192301,arXiv:1408.4075 [hep-ph]

  4. [4]

    ReconcilingJ/ψproduction at HERA, RHIC, Tevatron, and LHC with NRQCD factorization at next-to-leading order

    M. Butenschoen and B. A. Kniehl, “ReconcilingJ/ψproduction at HERA, RHIC, Tevatron, and LHC with NRQCD factorization at next-to-leading order”,Phys. Rev. Lett.106(2011) 022003, arXiv:1009.5662 [hep-ph]

  5. [5]

    Production and polarization of promptJ/ψin the improved color evaporation model using thek T -factorization approach

    V . Cheung and R. V ogt, “Production and polarization of promptJ/ψin the improved color evaporation model using thek T -factorization approach”,Phys. Rev. D98(2018) 114029, arXiv:1808.02909 [hep-ph]

  6. [6]

    Are there any challenges in the charmonia production and polarization at the LHC?

    S. P. Baranov and A. V . Lipatov, “Are there any challenges in the charmonia production and polarization at the LHC?”,Phys. Rev. D100(2019) 114021,arXiv:1906.07182 [hep-ph]

  7. [7]

    On relevance of triple gluon fusion inJ/ψhadroproduction

    L. Motyka and M. Sadzikowski, “On relevance of triple gluon fusion inJ/ψhadroproduction”, Eur . Phys. J. C75(2015) 213,arXiv:1501.04915 [hep-ph]

  8. [8]

    J/ψproduction in hadron scattering: three-pomeron contribution

    E. Levin and M. Siddikov, “J/ψproduction in hadron scattering: three-pomeron contribution”, Eur . Phys. J. C79(2019) 376,arXiv:1812.06783 [hep-ph]. [9]ALICECollaboration, S. Acharyaet al., “Charm production and fragmentation fractions at midrapidity in pp collisions at √s= 13 TeV”,JHEP12(2023) 086,arXiv:2308.04877 [hep-ex]. [10]ALICECollaboration, S. Acha...

  9. [9]

    A Decade of Collectivity in Small Systems,

    J. F. Grosse-Oetringhaus and U. A. Wiedemann, “A Decade of Collectivity in Small Systems”, arXiv:2407.07484 [hep-ex]. [15]ALICECollaboration, J. Adamet al., “Enhanced production of multi-strange hadrons in high-multiplicity proton-proton collisions”,Nature Phys.13(2017) 535–539, arXiv:1606.07424 [nucl-ex]. [16]ALICECollaboration, S. Acharyaet al., “The AL...

  10. [10]

    Event engineering studies for heavy flavor production and hadronization in high multiplicity hadron-hadron and hadron-nucleus collisions

    Y .-Q. Ma, P. Tribedy, R. Venugopalan, and K. Watanabe, “Event engineering studies for heavy flavor production and hadronization in high multiplicity hadron-hadron and hadron-nucleus collisions”,Phys. Rev. D98(2018) 074025,arXiv:1803.11093 [hep-ph]

  11. [11]

    Accessing subnuclear fluctuations and saturation with multiplicity dependent J/ψproduction in p+p and p+Pb collisions

    F. Salazar, B. Schenke, and A. Soto-Ontoso, “Accessing subnuclear fluctuations and saturation with multiplicity dependent J/ψproduction in p+p and p+Pb collisions”,Phys. Lett. B827(2022) 136952,arXiv:2112.04611 [hep-ph]

  12. [12]

    Multiplicity dependence of quarkonia production in the CGC approach

    E. Levin, I. Schmidt, and M. Siddikov, “Multiplicity dependence of quarkonia production in the CGC approach”,Eur . Phys. J. C80(2020) 560,arXiv:1910.13579 [hep-ph]

  13. [13]

    High energy QCD: multiplicity dependence of quarkonia production

    E. Gotsman and E. Levin, “High energy QCD: multiplicity dependence of quarkonia production”,Eur . Phys. J. C81(2021) 99,arXiv:2008.10911 [hep-ph]

  14. [14]

    Charmonium production in high multiplicity pp collisions and the structure of the proton

    R. Terra and F. S. Navarra, “Charmonium production in high multiplicity pp collisions and the structure of the proton”,Phys. Rev. D108(2023) 054002,arXiv:2306.14298 [hep-ph]

  15. [15]

    J/ψin high-multiplicity pp collisions: Lessons from pA collisions

    B. Z. Kopeliovich, H. J. Pirner, I. K. Potashnikova, K. Reygers, and I. Schmidt, “J/ψin high-multiplicity pp collisions: Lessons from pA collisions”,Phys. Rev. D88(2013) 116002, arXiv:1308.3638 [hep-ph]

  16. [16]

    Heavy quarkonium in the saturated environment of high-multiplicity pp collisions

    B. Z. Kopeliovich, H. J. Pirner, I. K. Potashnikova, K. Reygers, and I. Schmidt, “Heavy quarkonium in the saturated environment of high-multiplicity pp collisions”,Phys. Rev. D101 (2020) 054023,arXiv:1910.09682 [hep-ph]

  17. [17]

    High multiplicityppevents andJ/ψproduction at LHC

    E. G. Ferreiro and C. Pajares, “High multiplicityppevents andJ/ψproduction at LHC”,Phys. Rev. C86(2012) 034903,arXiv:1203.5936 [hep-ph]

  18. [18]

    A comprehensive guide to the physics and usage of PYTHIA 8.3

    C. Bierlichet al., “A comprehensive guide to the physics and usage of PYTHIA 8.3”,SciPost Phys. Codeb.(2022) 8,arXiv:2203.11601 [hep-ph]

  19. [19]

    String Formation Beyond Leading Colour

    J. R. Christiansen and P. Z. Skands, “String Formation Beyond Leading Colour”,JHEP08 (2015) 003,arXiv:1505.01681 [hep-ph]. 26 Multiplicity-dependent prompt and non-prompt J/ψproduction at √s=13 TeV ALICE Collaboration

  20. [20]

    Werner, Phys

    K. Werner, “Revealing a deep connection between factorization and saturation: New insight into modeling high-energy proton-proton and nucleus-nucleus scattering in the EPOS4 framework”, Phys. Rev. C108(2023) 064903,arXiv:2301.12517 [hep-ph]

  21. [21]

    Perturbative QCD concerning light and heavy flavor in the EPOS4 framework

    K. Werner and B. Guiot, “Perturbative QCD concerning light and heavy flavor in the EPOS4 framework”,Phys. Rev. C108(2023) 034904,arXiv:2306.02396 [hep-ph]

  22. [22]

    Core-corona procedure and microcanonical hadronization to understand strangeness enhancement in proton-proton and heavy ion collisions in the EPOS4 framework

    K. Werner, “Core-corona procedure and microcanonical hadronization to understand strangeness enhancement in proton-proton and heavy ion collisions in the EPOS4 framework”,Phys. Rev. C 109(2024) 014910,arXiv:2306.10277 [hep-ph]

  23. [23]

    Measurements of inclusive J/ψproduction at midrapidity and forward rapidity in Pb–Pb collisions at √sNN = 5.02 TeV

    K. Werner, “Parallel scattering, saturation, and generalized Abramovskii-Gribov-Kancheli (AGK) theorem in the EPOS4 framework, with applications for heavy-ion collisions at √sNN of 5.02 TeV and 200 GeV”,Phys. Rev. C109(2024) 034918,arXiv:2310.09380 [hep-ph]. [52]ALICECollaboration, S. Acharyaet al., “Measurements of inclusive J/ψproduction at midrapidity ...

  24. [24]

    Probing Collective Effects in Hadronisation with the Extremes of the Underlying Event

    T. Martin, P. Skands, and S. Farrington, “Probing Collective Effects in Hadronisation with the Extremes of the Underlying Event”,Eur . Phys. J. C76(2016) 299,arXiv:1603.05298 [hep-ph]. [55]ALICECollaboration, S. Acharyaet al., “Charged-particle production as a function of the relative transverse activity classifier in pp, p–Pb, and Pb–Pb collisions at the...

  25. [25]

    Elucidating the multiplicity dependence of J/ψproduction in proton-proton collisions with PYTHIA8

    S. G. Weber, A. Dubla, A. Andronic, and A. Morsch, “Elucidating the multiplicity dependence of J/ψproduction in proton-proton collisions with PYTHIA8”,Eur . Phys. J. C79(2019) 36, arXiv:1811.07744 [nucl-th]. [58]ALICECollaboration, S. Acharyaet al., “J/ψ-hadron correlations at midrapidity in pp collisions at √s= 13 TeV”,JHEP07(2025) 023,arXiv:2409.04364 [...

  26. [26]

    NRQCD Confronts LHCb Data on Quarkonium Production within Jets

    R. Bain, L. Dai, A. Leibovich, Y . Makris, and T. Mehen, “NRQCD Confronts LHCb Data on Quarkonium Production within Jets”,Phys. Rev. Lett.119(2017) 032002,arXiv:1702.05525 [hep-ph]. 27 Multiplicity-dependent prompt and non-prompt J/ψproduction at √s=13 TeV ALICE Collaboration [63]ALICECollaboration, B. Abelevet al., “The ALICE experiment at the CERN LHC”,...

  27. [27]

    The ALICE TPC, a large 3-dimensional tracking device with fast readout for ultra-high multiplicity events

    J. Almeet al., “The ALICE TPC, a large 3-dimensional tracking device with fast readout for ultra-high multiplicity events”,Nucl. Instrum. Meth. A622(2010) 316–367,arXiv:1001.1950 [physics.ins-det]. [67]ALICECollaboration, E. Abbaset al., “Performance of the ALICE VZERO system”,JINST8 (2013) P10016,arXiv:1306.3130 [nucl-ex]. [68]ALICECollaboration, S. Acha...

  28. [28]

    Tuning PYTHIA 8.1: the Monash 2013 Tune

    P. Skands, S. Carrazza, and J. Rojo, “Tuning PYTHIA 8.1: the Monash 2013 Tune”,Eur . Phys. J. C74(2014) 3024,arXiv:1404.5630 [hep-ph]

  29. [29]

    Geant detector description and simulation tool

    R. Brunet al., “Geant detector description and simulation tool”, (1994) CERN-W-5013. [71]ALICECollaboration, S. Acharyaet al., “Multiplicity dependence of charged-particle production in pp, p–Pb, Xe–Xe and Pb–Pb collisions at the LHC”,Phys. Lett. B845(2023) 138110,arXiv:2211.15326 [nucl-ex]

  30. [30]

    PYTHIA 6.4 Physics and Manual

    T. Sjostrand, S. Mrenna, and P. Z. Skands, “PYTHIA 6.4 Physics and Manual”,JHEP05(2006) 026,arXiv:hep-ph/0603175

  31. [31]

    The EvtGen particle decay simulation package

    D. Lange, “The EvtGen particle decay simulation package”,Nucl. Instrum. Meth. A462(2001) 152–155

  32. [32]

    PHOTOS: A Universal Monte Carlo for QED radiative corrections. Version 2.0

    E. Barberio and Z. Was, “PHOTOS: A Universal Monte Carlo for QED radiative corrections. Version 2.0”,Comput. Phys. Commun.79(1994) 291–308. [75]ALICECollaboration, “The ALICE definition of primary particles”, ALICE-PUBLIC-2017-005. [76]ALICECollaboration, S. Acharyaet al., “Transverse momentum spectra and nuclear modification factors of charged particles ...

  33. [33]

    A multidimensional unfolding method based on bayes’ theorem

    G. D’Agostini, “A multidimensional unfolding method based on bayes’ theorem”,Nucl. Instrum. Meth. A362(1995) 487–498. [78]ALICECollaboration, S. Acharyaet al., “Inclusive J/ψproduction at midrapidity in pp collisions at √s=13 TeV”,Eur . Phys. J. C81(2021) 1121,arXiv:2108.01906 [nucl-ex]. [79]TMV ACollaboration, A. Hockeret al., “TMV A - Toolkit for Multiv...

  34. [34]

    Malaescu, (2009), arXiv:0907.3791 [physics.data-an]

    B. Malaescu, “An Iterative, dynamically stabilized method of data unfolding”, arXiv:0907.3791 [physics.data-an]

  35. [35]

    Cacciari, M

    M. Cacciari, M. Greco, and P. Nason, “Thep T spectrum in heavy-flavour hadroproduction.”, JHEP05(1998) 007,arXiv:hep-ph/9803400

  36. [36]

    Cacciari, S

    M. Cacciari, S. Frixione, N. Houdeau, M. L. Mangano, P. Nason, and G. Ridolfi, “Theoretical predictions for charm and bottom production at the LHC”,JHEP10(2012) 137, arXiv:1205.6344 [hep-ph]

  37. [37]

    Non-relativistic quantum chromodynamics in parton showers

    N. Cooke, P. Ilten, L. Lönnblad, and S. Mrenna, “Non-relativistic quantum chromodynamics in parton showers”,Eur . Phys. J. C84(2024) 432,arXiv:2312.05203 [hep-ph]

  38. [38]

    Hadronic Collisions: A Quarkonium Factory

    R. Baier and R. Ruckl, “Hadronic Collisions: A Quarkonium Factory”,Z. Phys. C19(1983) 251

  39. [39]

    Color octet quarkonia production. 2

    P. L. Cho and A. K. Leibovich, “Color octet quarkonia production. 2.”,Phys. Rev. D53(1996) 6203–6217,arXiv:hep-ph/9511315

  40. [40]

    Nasonet al., Bottom production, In Workshop on Standard Model Physics (and more) at the LHC (First Plenary Meeting), pp

    P. Nasonet al., “Bottom production”, inWorkshop on Standard Model Physics (and more) at the LHC (First Plenary Meeting), pp. 231–304. 5, 1999.arXiv:hep-ph/0003142

  41. [41]

    Heavy flavor as a probe of hot QCD matter produced in proton-proton collisions

    J. Zhao, J. Aichelin, P. B. Gossiaux, and K. Werner, “Heavy flavor as a probe of hot QCD matter produced in proton-proton collisions”,Phys. Rev. D109(2024) 054011,arXiv:2310.08684 [hep-ph]

  42. [42]

    Heavy-flavor hadron production in relativistic heavy ion collisions at energies available at BNL RHIC and at the CERN LHC in the EPOS4HQ framework

    J. Zhao, J. Aichelin, P. B. Gossiaux, V . Ozvenchuk, and K. Werner, “Heavy-flavor hadron production in relativistic heavy ion collisions at energies available at BNL RHIC and at the CERN LHC in the EPOS4HQ framework”,Phys. Rev. C110(2024) 024909, arXiv:2401.17096 [hep-ph]

  43. [43]

    Correlations between heavy mesons and the creation of the charmonia, bottomonia, andB c mesons in high energyppcollisions

    J. Zhao, J. Aichelin, P. B. Gossiaux, and K. Werner, “Correlations between heavy mesons and the creation of the charmonia, bottomonia, andB c mesons in high energyppcollisions”, arXiv:2511.08796 [hep-ph]. [95]ALICECollaboration, S. Acharyaet al., “Search for jet quenching effects in high-multiplicity pp collisions at √s= 13 TeV via di-jet acoplanarity”,JH...

  44. [44]

    The multiple-charm hierarchy in the statistical hadronization model,

    A. Andronic, P. Braun-Munzinger, M. K. Köhler, A. Mazeliauskas, K. Redlich, J. Stachel, and V . Vislavicius, “The multiple-charm hierarchy in the statistical hadronization model”,JHEP07 (2021) 035,arXiv:2104.12754 [hep-ph]. [99]ALICECollaboration, S. Acharyaet al., “Pseudorapidity distributions of charged particles as a function of mid- and forward rapidi...