Top: Self-normalized yield of $\Upsilon$ as a function of self-normalized charged-particle multiplicity, compared to inclusive J/$\psi$ measured in the forward rapidity region at 5.02~TeV~, 7~TeV~, and 13~TeV~, and to inclusive J/$\psi$ measured in the central rapidity region at 13~TeV~. The error bars represent the statistical uncertainty on the quarkonium yields, while the quadratic sum of the point-by-point systematic uncertainties on the quarkonium yields as well as on $\textrm{d}N_{\mathrm{ch}} / \textrm{d}\eta$ / $\langle \textrm{d}N_{\mathrm{ch}}/\textrm{d}\eta \rangle$ is depicted as boxes. Bottom: Self-normalized yield ratio of $\Upsilon$(1S) over J/$\psi$ as a function of self-normalized charged-particle multiplicity, compared to model computations from 3-pomeron CGC approach~, PYTHIA 8.2~, comovers~, and CPP~..