Examination of thermalization of quarkonia at energies available at the CERN Large Hadron Collider

We analyze the relative yields of different bottomonia and charmonia states produced in Pb-Pb, p-Pb and high multiplicity p-p collisions at LHC, within a semi-classical grand canonical ensemble approach. The underlying assumption is the early thermalization and subsequent freezeout of these heavy hadrons resulting in their chemical freezeout at a temperature of approximately $230$~MeV, significantly higher than that of light and strange hadrons. The systematic dependence of the freezeout temperature on the collision centrality is also investigated in details.