Abstract We analyze the relative yields of different bottomonia and charmonia states produced in Pb-Pb, p−Pb, and high multiplicity p−p collisions at the CERN Large Hadron Collider, within a semiclassical grand canonical ensemble approach. The underlying assumption is the early thermalization and subsequent freeze-out of these heavy hadrons resulting in their chemical freeze-out at a temperature of approximately 230MeV, significantly higher than that of light and strange hadrons.
