The total coherent $\mathrm{J}/\psi$ photoproduction cross section as a function of photon-nuclear center-of-mass energy per nucleon $W_{\gamma \mathrm{N}}^{\mathrm{Pb}}$, measured in PbPb ultra-peripheral collisions at $\sqrt{s_{\mathrm{NN}}}$ = 5.02 TeV. The $W_{\gamma \mathrm{N}}^{\mathrm{Pb}}$ values used correspond to the center of each rapidity range. The theoretical uncertainties is due to the uncertainties in the photon flux.

The nuclear gluon suppression factor $R_{\mathrm{g}}^{\mathrm{Pb}}$ as a function of Bjorken $x$ extracted from the CMS measurement of the coherent $\mathrm{J}/\psi$ photoproduction in PbPb ultra-peripheral collisions at $\sqrt{s_{\mathrm{NN}}}$ = 5.02 TeV. The $x$ values are evaluated at the centers of their corresponding rapidity ranges. The theoretical uncertainties are due to the uncertainties in the photon flux and the impulse approximation model.

The nuclear gluon suppression factor $R_{\mathrm{g}}^{\mathrm{Pb}}$ as a function of Bjorken $x$ extracted from the CMS measurement of the coherent $\mathrm{J}/\psi$ photoproduction in PbPb ultra-peripheral collisions at $\sqrt{s_{\mathrm{NN}}}$ = 5.02 TeV. The $x$ values are evaluated at the centers of their corresponding rapidity ranges. The theoretical uncertainties are due to the uncertainties in the photon flux and the impulse approximation model.

The total covariance matrix of the total coherent photoproduction cross section as a function of photon-nuclear center-of-mass energy per nucleon $W_{\gamma \mathrm{N}}^{\mathrm{Pb}}$. The covariance matrix includes both the experimental and theoretical (photon flux) uncertainties. The bins are ordered as increasing in $W_{\gamma \mathrm{N}}^{\mathrm{Pb}}$.

The photon flux values from STARLight as a function of rapidity, in different neutron multiplicity classes: 0n0n, 0nXn, and XnXn.

The experimental covariance matrix of the total coherent photoproduction cross section as a function of photon-nuclear center-of-mass energy per nucleon $W_{\gamma \mathrm{N}}^{\mathrm{Pb}}$. The bins are ordered as increasing in $W_{\gamma \mathrm{N}}^{\mathrm{Pb}}$.

The covariance matrix for the flux in the 0n0n neutron multiplicity class.

The theoretical (photon flux) covariance matrix of the total coherent photoproduction cross section as a function of photon-nuclear center-of-mass energy per nucleon $W_{\gamma \mathrm{N}}^{\mathrm{Pb}}$. The bins are ordered as increasing in $W_{\gamma \mathrm{N}}^{\mathrm{Pb}}$.

The covariance matrix for the flux in the 0nXn neutron multiplicity class.

The covariance matrix for the flux in the XnXn neutron multiplicity class.

The total covariance matrix of the total coherent photoproduction cross section as a function of photon-nuclear center-of-mass energy per nucleon $W_{\gamma \mathrm{N}}^{\mathrm{Pb}}$. The covariance matrix includes both the experimental and theoretical (photon flux) uncertainties. The bins are ordered as increasing in $W_{\gamma \mathrm{N}}^{\mathrm{Pb}}$.

The experimental covariance matrix of the total coherent photoproduction cross section as a function of photon-nuclear center-of-mass energy per nucleon $W_{\gamma \mathrm{N}}^{\mathrm{Pb}}$. The bins are ordered as increasing in $W_{\gamma \mathrm{N}}^{\mathrm{Pb}}$.

The theoretical (photon flux) covariance matrix of the total coherent photoproduction cross section as a function of photon-nuclear center-of-mass energy per nucleon $W_{\gamma \mathrm{N}}^{\mathrm{Pb}}$. The bins are ordered as increasing in $W_{\gamma \mathrm{N}}^{\mathrm{Pb}}$.

$\phi$ meson production cross section $\mathrm{d}^2\sigma/(\mathrm{d}y\mathrm{d}p_\mathrm{T})$ as a function of $p_\mathrm{T}$ at $\sqrt{s}=5.02$ TeV in pp collisions, in the $2.5<y<3$ rapidity interval.

$\phi$ meson production cross section $\mathrm{d}^2\sigma/(\mathrm{d}y\mathrm{d}p_\mathrm{T})$ as a function of $p_\mathrm{T}$ at $\sqrt{s}=5.02$ TeV in pp collisions, in the $3<y<3.25$ rapidity interval.

$\phi$ meson production cross section $\mathrm{d}^2\sigma/(\mathrm{d}y\mathrm{d}p_\mathrm{T})$ as a function of $p_\mathrm{T}$ at $\sqrt{s}=5.02$ TeV in pp collisions, in the $3.25<y<3.5$ rapidity interval.

$\phi$ meson production cross section $\mathrm{d}^2\sigma/(\mathrm{d}y\mathrm{d}p_\mathrm{T})$ as a function of $p_\mathrm{T}$ at $\sqrt{s}=5.02$ TeV in pp collisions, in the $3.5<y<4$ rapidity interval.

$\phi$ meson production cross section $\mathrm{d}^2\sigma/(\mathrm{d}y\mathrm{d}p_\mathrm{T})$ as a function of $p_\mathrm{T}$ at $\sqrt{s}=13$ TeV in pp collisions, in the $2.5<y<2.75$ rapidity interval.

$\phi$ meson production cross section $\mathrm{d}^2\sigma/(\mathrm{d}y\mathrm{d}p_\mathrm{T})$ as a function of $p_\mathrm{T}$ at $\sqrt{s}=13$ TeV in pp collisions, in the $2.75<y<3$ rapidity interval.

$\phi$ meson production cross section $\mathrm{d}^2\sigma/(\mathrm{d}y\mathrm{d}p_\mathrm{T})$ as a function of $p_\mathrm{T}$ at $\sqrt{s}=13$ TeV in pp collisions, in the $3<y<3.25$ rapidity interval.

$\phi$ meson production cross section $\mathrm{d}^2\sigma/(\mathrm{d}y\mathrm{d}p_\mathrm{T})$ as a function of $p_\mathrm{T}$ at $\sqrt{s}=13$ TeV in pp collisions, in the $3.25<y<3.5$ rapidity interval.

$\phi$ meson production cross section $\mathrm{d}^2\sigma/(\mathrm{d}y\mathrm{d}p_\mathrm{T})$ as a function of $p_\mathrm{T}$ at $\sqrt{s}=13$ TeV in pp collisions, in the $3.5<y<4$ rapidity interval.

$\phi$ meson production cross section $\mathrm{d}^2\sigma/(\mathrm{d}y\mathrm{d}p_\mathrm{T})$ as a function of rapidity at $\sqrt{s}=5.02$ TeV in pp collisions, in the $1<p_T<2$ GeV/$c$ interval.

$\phi$ meson production cross section $\mathrm{d}^2\sigma/(\mathrm{d}y\mathrm{d}p_\mathrm{T})$ as a function of rapidity at $\sqrt{s}=5.02$ TeV in pp collisions, in the $2<p_T<5$ GeV/$c$ interval.

$\phi$ meson production cross section $\mathrm{d}^2\sigma/(\mathrm{d}y\mathrm{d}p_\mathrm{T})$ as a function of rapidity at $\sqrt{s}=5.02$ TeV in pp collisions, in the $5<p_T<8$ GeV/$c$ interval.

$\phi$ meson production cross section $\mathrm{d}^2\sigma/(\mathrm{d}y\mathrm{d}p_\mathrm{T})$ as a function of rapidity at $\sqrt{s}=13$ TeV in pp collisions, in the $1<p_T<2$ GeV/$c$ interval.

$\phi$ meson production cross section $\mathrm{d}^2\sigma/(\mathrm{d}y\mathrm{d}p_\mathrm{T})$ as a function of rapidity at $\sqrt{s}=13$ TeV in pp collisions, in the $2<p_T<5$ GeV/$c$ interval.

$\phi$ meson production cross section $\mathrm{d}^2\sigma/(\mathrm{d}y\mathrm{d}p_\mathrm{T})$ as a function of rapidity at $\sqrt{s}=13$ TeV in pp collisions, in the $5<p_T<10$ GeV/$c$ interval.

$\phi$ meson production integrated cross section as a function of $\sqrt{s}$ for $1.5 < p_\mathrm{T} <5$ GeV/$c$ at forward rapidity in pp collisions.

The results of fits of the two-particle double-ratio correlation function $R_{2} (Q)$ for $\sqrt{s} = 0.9\ TeV$ events with the unlike-charge reference sample for various multiplicity intervals $n_{ch}$ for the exponential parametrization $\Omega^{(E)}$. The uncertainties for parameters $\lambda$ and $R$ are square root from the quadratic sum of statistical and systematic errors. Where only one error is shown for parameters $C_{0}$ and $\epsilon$, this represents the statistical uncertainty only.

The results of fits of the two-particle double-ratio correlation function $R_{2} (Q)$ for $\sqrt{s} = 7\ TeV$ events with the unlike-charge reference sample for various multiplicity intervals $n_{ch}$ for the exponential parametrization $\Omega^{(E)}$. The uncertainties for parameters $\lambda$ and $R$ are square root from the quadratic sum of statistical and systematic errors. The statistical uncertainties at $7\ TeV$ were corrected on the $\sqrt{\chi^{2}/ndf}$ and they more smaller than systematic uncertainties. Where only one error is shown for parameters $C_{0}$ and $\epsilon$, this represents the statistical uncertainty only.

The results of fits of the two-particle double-ratio correlation function $R_{2} (Q)$ for $\sqrt{s} = 7\ TeV\ HM$ events with the unlike-charge reference sample for various multiplicity intervals $n_{ch}$ for the exponential parametrization $\Omega^{(E)}$. The uncertainties for parameters $\lambda$ and $R$ are square root from the quadratic sum of statistical and systematic errors. Where only one error is shown for parameters $C_{0}$ and $\epsilon$, this represents the statistical uncertainty only.

The results of fits of the two-particle double-ratio correlation function $R_{2} (Q)$ for $\sqrt{s} = 0.9\ TeV$ ($n_{ch} \ge 2$) events with the unlike-charge reference sample for various $k_{T}$ intervals for the exponential parametrization $\Omega^{(E)}$. The uncertainties for parameters $\lambda$ and $R$ are square root from the quadratic sum of statistical and systematic errors. The statistical uncertainties were corrected on the $\sqrt{\chi^{2}/ndf}$ when $\chi^{2}/ndf > 1$. Where only one error is shown for parameters $C_{0}$ and $\epsilon$, this represents the statistical uncertainty only.

The results of fits of the two-particle double-ratio correlation function $R_{2} (Q)$ for $\sqrt{s} = 7\ TeV$ ($n_{ch} \ge 2$) events with the unlike-charge reference sample for various $k_{T}$ intervals for the exponential parametrization $\Omega^{(E)}$. The uncertainties for parameters $\lambda$ and $R$ are square root from the quadratic sum of statistical and systematic errors. Where only one error is shown for parameters $C_{0}$ and $\epsilon$ , this represents the statistical uncertainty only. The statistical uncertainties were corrected on the $\sqrt{\chi^{2}/ndf}$ when $\chi^{2}/ndf > 1$. The statistical uncertainties at $7\ TeV$ more smaller than systematic uncertainties.

The results of fits of the two-particle double-ratio correlation function $R_{2} (Q)$ for $\sqrt{s} = 7\ TeV\ HM$ ($n_{ch} \ge 150$) events with the unlike-charge reference sample for various $k_{T}$ intervals for the exponential parametrization $\Omega^{(E)}$. The uncertainties for parameters $\lambda$ and $R$ are square root from the quadratic sum of statistical and systematic errors. The statistical uncertainties were corrected on the $\sqrt{\chi^{2}/ndf}$ when $\chi^{2}/ndf > 1$.Where only one error is shown for parameters $C_{0}$ and $\epsilon$, this represents the statistical uncertainty only.

The results of fits of the two-particle double-ratio correlation function $R_{2} (Q)$ for $\sqrt{s} = 7\ TeV$ events with the unlike-charge reference sample for various $k_{T}$ intervals of different multiplicity regions $n_{ch}$ for the exponential parametrization $\Omega^{(E)}$. The statistical uncertainties were corrected on the $\sqrt{\chi^{2}/ndf}$ when $\chi^{2}/ndf > 1$. The statistical uncertainties smaller than systematic uncertainties. The uncertainties for parameters $\lambda$ and $R$ are square root from the quadratic sum of statistical and systematic errors. Where only one error is shown for parameters $C_{0}$ and $\epsilon$, this represents the statistical uncertainty only.

The $Q$ distribution measured at $0.9\ TeV$ for unlike-sign pairs for $p_{T} > 100\ MeV$ and $|\eta|<2.5$, renormalized to the total number of charged-particle pairs. The error bars represents the statistical uncertainties.

The $Q$ distribution measured at $0.9\ TeV$ for like-sign pairs for $p_{T} > 100\ MeV$ and $|\eta|<2.5$, renormalized to the total number of charged-particle pairs. The error bars represents the statistical uncertainties.

The $Q$ distribution measured at $7\ TeV$ for unlike-sign pairs for $p_{T} > 100\ MeV$ and $|\eta|<2.5$, renormalized to the total number of charged-particle pairs. The error bars represents the statistical uncertainties.

The $Q$ distribution measured at $7\ TeV$ for like-sign pairs for $p_{T} > 100\ MeV$ and $|\eta|<2.5$, renormalized to the total number of charged-particle pairs. The error bars represents the statistical uncertainties.

The $R_{2}(Q)$ correlation function measured at $0.9\ TeV$ using unlike-charge particle reference sample for different $n_{ch}$ intervals $p_{T} > 100\ MeV$ and $|\eta|<2.5$. The error bars represents only the statistical uncertainties.

The $R_{2}(Q)$ correlation function measured at $0.9\ TeV$ using unlike-charge particle reference sample for different $k_{T}$ intervals with multiplicity $n_{ch} \ge 2$. The error bars represents only the statistical uncertainties.

The $R_{2}(Q)$ correlation function measured at $7\ TeV$ using unlike-charge particle reference sample for different multiplicity intervals $n_{ch}$ with $p_{T} > 100\ MeV$ and $|\eta| < 2.5$. The error bars represents only the statistical uncertainties.

The $R_{2}(Q)$ correlation function measured at $7\ TeV$ using unlike-charge particle reference sample for different $k_{T}$ intervals with multiplicity $n_{ch} \ge 2$. The error bars represents only the statistical uncertainties.

The $R_{2}(Q)$ correlation function measured at $7\ TeV\ HM$ using unlike-charge particle reference sample for different $n_{ch}$ intervals $p_{T} > 100\ MeV$ and $|\eta|<2.5$. The error bars represents only the statistical uncertainties.

The $R_{2}(Q)$ correlation function measured with $7\ TeV\ HM$ data using unlike-charge particle reference sample for different $k_{T}$ intervals with multiplicity $n_{ch} \ge 150$. The error bars represents only the statistical uncertainties.

The $R_{2}(Q)$ correlation function measured at $7\ TeV$ using unlike-charge particle reference sample for different $k_{T}$ intervals within multiplicity interval $n_{ch} = 2-9$. The error bars represents only the statistical uncertainties.

The $R_{2}(Q)$ correlation function measured at $7\ TeV$ using unlike-charge particle reference sample for different $k_{T}$ intervals within multiplicity interval $n_{ch} = 10-24$. The error bars represents only the statistical uncertainties.

The $R_{2}(Q)$ correlation function measured at $7\ TeV$ using unlike-charge particle reference sample for different $k_{T}$ intervals within multiplicity interval $n_{ch} = 25-80$. The error bars represents only the statistical uncertainties.

The $R_{2}(Q)$ correlation function measured at $7\ TeV$ using unlike-charge particle reference sample for different $k_{T}$ intervals within multiplicity interval $n_{ch} = 81-125$. The error bars represents only the statistical uncertainties.