# The production of primordial $J/\psi$ in p+p and relativistic heavy-ion collisions

## TL;DR

This paper models the production of primordial J/ψ mesons in proton-proton and heavy-ion collisions using PYTHIA and Wigner density formalism, successfully reproducing experimental data and exploring QGP effects.

## Contribution

It introduces a novel approach combining PYTHIA and Wigner density formalism to accurately predict J/ψ production and its distributions in high-energy collisions.

## Key findings

- Reproduces experimental J/ψ yield, rapidity, and pT distributions.
- Finds significant enhancement of primordial J/ψ in heavy-ion collisions without QGP effects.
- Validates the modeling approach against available data.

## Abstract

$J/\psi$ observables of p+p collisions at $\sqrt{s}=$ 200 GeV and 2.76 TeV are calculated using the PYTHIA event generator for the production of $c$ and $\bar c$ and the Wigner density formalism. The charm quark momentum from PYTHIA is tuned in such a way that the transverse momentum and rapidity distribution correspond to the Fixed-Order Next-to-Leading Logarithm (FONLL) calculations and hence to the experimental data for open charm mesons. Using the Wigner density of the charmonium wave function we calculate for each charm quark pair the probability to form a charmonium. As a result, the experimental data on the total $J/\psi$ yield, its rapidity and transverse momentum distribution as well as the fractions of feed-down from excited states of charmonium are reproduced. Applying the same approach to relativistic heavy-ion collisions, we find, if quark gluon plasma (QGP) effects are ignored, a considerable enhancement of the primordial $J/\psi$'s.