Ratio of the $\eta/\pi^{0}$ ratio to the $m_{\rm T}$-scaling prediction as a function of $p_{\rm T}$ for pp collisions at $\sqrt{s}$ = 13 TeV.

$x_{\rm T}$ scaled $\pi^{0}$ spectrum for pp collisions at $\sqrt{s}$ = 0.9 TeV, with $n = 5.01$.

$x_{\rm T}$ scaled $\pi^{0}$ spectrum for pp collisions at $\sqrt{s}$ = 2.76 TeV, with $n = 5.01$.

$x_{\rm T}$ scaled $\pi^{0}$ spectrum for pp collisions at $\sqrt{s}$ = 7 TeV, with $n = 5.01$.

$x_{\rm T}$ scaled $\pi^{0}$ spectrum for pp collisions at $\sqrt{s}$ = 8 TeV, with $n = 5.01$.

$x_{\rm T}$ scaled $\pi^{0}$ spectrum for pp collisions at $\sqrt{s}$ = 13 TeV, with $n = 5.01$.

$x_{\rm T}$ scaled $\eta$ spectrum for pp collisions at $\sqrt{s}$ = 2.76 TeV, with $n = 4.91$.

$x_{\rm T}$ scaled $\eta$ spectrum for pp collisions at $\sqrt{s}$ = 7 TeV, with $n = 4.91$.

$x_{\rm T}$ scaled $\eta$ spectrum for pp collisions at $\sqrt{s}$ = 8 TeV, with $n = 4.91$.

$x_{\rm T}$ scaled $\eta$ spectrum for pp collisions at $\sqrt{s}$ = 13 TeV, with $n = 4.91$.

Invariant differential yields of the $\pi^{0}$ meson versus transverse momentum for pp collisions at $\sqrt{s}$ = 13 TeV in the INEL>0 event class in the V0M multiplicity class 0-100%. (Scaled for visibility in the figure)

Invariant differential yields of the $\pi^{0}$ meson versus transverse momentum for pp collisions at $\sqrt{s}$ = 13 TeV in the INEL>0 event class in the V0M multiplicity class 70-100%. (Scaled for visibility in the figure)

Invariant differential yields of the $\pi^{0}$ meson versus transverse momentum for pp collisions at $\sqrt{s}$ = 13 TeV in the INEL>0 event class in the V0M multiplicity class 50-70%. (Scaled for visibility in the figure)

Invariant differential yields of the $\pi^{0}$ meson versus transverse momentum for pp collisions at $\sqrt{s}$ = 13 TeV in the INEL>0 event class in the V0M multiplicity class 30-50%. (Scaled for visibility in the figure)

Invariant differential yields of the $\pi^{0}$ meson versus transverse momentum for pp collisions at $\sqrt{s}$ = 13 TeV in the INEL>0 event class in the V0M multiplicity class 20-30%. (Scaled for visibility in the figure)

Invariant differential yields of the $\pi^{0}$ meson versus transverse momentum for pp collisions at $\sqrt{s}$ = 13 TeV in the INEL>0 event class in the V0M multiplicity class 10-20%. (Scaled for visibility in the figure)

Invariant differential yields of the $\pi^{0}$ meson versus transverse momentum for pp collisions at $\sqrt{s}$ = 13 TeV in the INEL>0 event class in the V0M multiplicity class 5-10%. (Scaled for visibility in the figure)

Invariant differential yields of the $\pi^{0}$ meson versus transverse momentum for pp collisions at $\sqrt{s}$ = 13 TeV in the INEL>0 event class in the V0M multiplicity class 1-5%. (Scaled for visibility in the figure)

Invariant differential yields of the $\pi^{0}$ meson versus transverse momentum for pp collisions at $\sqrt{s}$ = 13 TeV in the INEL>0 event class in the V0M multiplicity class 0-1%. (Scaled for visibility in the figure)

Invariant differential yields of the $\pi^{0}$ meson versus transverse momentum for pp collisions at $\sqrt{s}$ = 13 TeV in the INEL>0 event class in the V0M multiplicity class 0.05-0.1%. (Scaled for visibility in the figure)

Invariant differential yields of the $\pi^{0}$ meson versus transverse momentum for pp collisions at $\sqrt{s}$ = 13 TeV in the INEL>0 event class in the V0M multiplicity class 0-0.1%. (Scaled for visibility in the figure)

Invariant differential yields of the $\pi^{0}$ meson versus transverse momentum for pp collisions at $\sqrt{s}$ = 13 TeV in the INEL>0 event class in the V0M multiplicity class 0.01-0.05%. (Scaled for visibility in the figure)

Invariant differential yields of the $\pi^{0}$ meson versus transverse momentum for pp collisions at $\sqrt{s}$ = 13 TeV in the INEL>0 event class in the V0M multiplicity class 0-0.01%. (Scaled for visibility in the figure)

Invariant differential yields of the $\eta$ meson versus transverse momentum for pp collisions at $\sqrt{s}$ = 13 TeV in the INEL>0 event class in the V0M multiplicity class 0-100%. (Scaled for visibility in the figure)

Invariant differential yields of the $\eta$ meson versus transverse momentum for pp collisions at $\sqrt{s}$ = 13 TeV in the INEL>0 event class in the V0M multiplicity class 70-100%. (Scaled for visibility in the figure)

Invariant differential yields of the $\eta$ meson versus transverse momentum for pp collisions at $\sqrt{s}$ = 13 TeV in the INEL>0 event class in the V0M multiplicity class 50-70%. (Scaled for visibility in the figure)

Invariant differential yields of the $\eta$ meson versus transverse momentum for pp collisions at $\sqrt{s}$ = 13 TeV in the INEL>0 event class in the V0M multiplicity class 30-50%. (Scaled for visibility in the figure)

Invariant differential yields of the $\eta$ meson versus transverse momentum for pp collisions at $\sqrt{s}$ = 13 TeV in the INEL>0 event class in the V0M multiplicity class 20-30%. (Scaled for visibility in the figure)

Invariant differential yields of the $\eta$ meson versus transverse momentum for pp collisions at $\sqrt{s}$ = 13 TeV in the INEL>0 event class in the V0M multiplicity class 10-20%. (Scaled for visibility in the figure)

Invariant differential yields of the $\eta$ meson versus transverse momentum for pp collisions at $\sqrt{s}$ = 13 TeV in the INEL>0 event class in the V0M multiplicity class 5-10%. (Scaled for visibility in the figure)

Invariant differential yields of the $\eta$ meson versus transverse momentum for pp collisions at $\sqrt{s}$ = 13 TeV in the INEL>0 event class in the V0M multiplicity class 1-5%. (Scaled for visibility in the figure)

Invariant differential yields of the $\eta$ meson versus transverse momentum for pp collisions at $\sqrt{s}$ = 13 TeV in the INEL>0 event class in the V0M multiplicity class 0-1%. (Scaled for visibility in the figure)

Invariant differential yields of the $\eta$ meson versus transverse momentum for pp collisions at $\sqrt{s}$ = 13 TeV in the INEL>0 event class in the V0M multiplicity class 0.05-0.1%. (Scaled for visibility in the figure)

Invariant differential yields of the $\eta$ meson versus transverse momentum for pp collisions at $\sqrt{s}$ = 13 TeV in the INEL>0 event class in the V0M multiplicity class 0-0.1%. (Scaled for visibility in the figure)

Invariant differential yields of the $\eta$ meson versus transverse momentum for pp collisions at $\sqrt{s}$ = 13 TeV in the INEL>0 event class in the V0M multiplicity class 0.01-0.05%. (Scaled for visibility in the figure)

Invariant differential yields of the $\eta$ meson versus transverse momentum for pp collisions at $\sqrt{s}$ = 13 TeV in the INEL>0 event class in the V0M multiplicity class 0-0.01%. (Scaled for visibility in the figure)

Ratios of the invariant differential yields of $\pi^{0}$ meson in the V0M multiplicity class 70-100% to the inclusive spectrum in the INEL>0 event class.

Ratios of the invariant differential yields of $\pi^{0}$ meson in the V0M multiplicity class 50-70% to the inclusive spectrum in the INEL>0 event class.

Ratios of the invariant differential yields of $\pi^{0}$ meson in the V0M multiplicity class 30-50% to the inclusive spectrum in the INEL>0 event class.

Ratios of the invariant differential yields of $\pi^{0}$ meson in the V0M multiplicity class 20-30% to the inclusive spectrum in the INEL>0 event class (data points not shown in figure 11).

Ratios of the invariant differential yields of $\pi^{0}$ meson in the V0M multiplicity class 10-20% to the inclusive spectrum in the INEL>0 event class.

Ratios of the invariant differential yields of $\pi^{0}$ meson in the V0M multiplicity class 5-10% to the inclusive spectrum in the INEL>0 event class (data points not shown in figure 11).

Ratios of the invariant differential yields of $\pi^{0}$ meson in the V0M multiplicity class 1-5% to the inclusive spectrum in the INEL>0 event class.

Ratios of the invariant differential yields of $\pi^{0}$ meson in the V0M multiplicity class 0-1% to the inclusive spectrum in the INEL>0 event class (data points not shown in figure 11).

Ratios of the invariant differential yields of $\pi^{0}$ meson in the V0M multiplicity class 0.05-0.1% to the inclusive spectrum in the INEL>0 event class (data points not shown in figure 11).

Ratios of the invariant differential yields of $\pi^{0}$ meson in the V0M multiplicity class 0-0.1% to the inclusive spectrum in the INEL>0 event class (data points not shown in figure 11).

Ratios of the invariant differential yields of $\pi^{0}$ meson in the V0M multiplicity class 0.01-0.05% to the inclusive spectrum in the INEL>0 event class (data points not shown in figure 11).

Ratios of the invariant differential yields of $\pi^{0}$ meson in the V0M multiplicity class 0-0.01% to the inclusive spectrum in the INEL>0 event class.

Ratios of the invariant differential yields of $\eta$ meson in the V0M multiplicity class 70-100% to the inclusive spectrum in the INEL>0 event class.

Ratios of the invariant differential yields of $\eta$ meson in the V0M multiplicity class 50-70% to the inclusive spectrum in the INEL>0 event class.

Ratios of the invariant differential yields of $\eta$ meson in the V0M multiplicity class 30-50% to the inclusive spectrum in the INEL>0 event class.

Ratios of the invariant differential yields of $\eta$ meson in the V0M multiplicity class 20-30% to the inclusive spectrum in the INEL>0 event class (data points not shown in figure 11).

Ratios of the invariant differential yields of $\eta$ meson in the V0M multiplicity class 10-20% to the inclusive spectrum in the INEL>0 event class.

Ratios of the invariant differential yields of $\eta$ meson in the V0M multiplicity class 5-10% to the inclusive spectrum in the INEL>0 event class (data points not shown in figure 11).

Ratios of the invariant differential yields of $\eta$ meson in the V0M multiplicity class 1-5% to the inclusive spectrum in the INEL>0 event class.

Ratios of the invariant differential yields of $\eta$ meson in the V0M multiplicity class 0-1% to the inclusive spectrum in the INEL>0 event class (data points not shown in figure 11).

Ratios of the invariant differential yields of $\eta$ meson in the V0M multiplicity class 0.05-0.1% to the inclusive spectrum in the INEL>0 event class (data points not shown in figure 11).

Ratios of the invariant differential yields of $\eta$ meson in the V0M multiplicity class 0-0.1% to the inclusive spectrum in the INEL>0 event class (data points not shown in figure 11).

Ratios of the invariant differential yields of $\eta$ meson in the V0M multiplicity class 0.01-0.05% to the inclusive spectrum in the INEL>0 event class (data points not shown in figure 11).

Ratios of the invariant differential yields of $\eta$ meson in the V0M multiplicity class 0-0.01% to the inclusive spectrum in the INEL>0 event class.

Slope parameter (n) of a power-law fit ($f(p_{\rm T}) = \alpha p_{\rm T}^{-n}$) for $5 < p_{\rm T} < 15$ GeV/$c$ as a function the charged-particle multiplicity density in units of the average normalized multiplicity density for the INEL>0 event class for the neutral pion.

Slope parameter (k) of an exponential fit ($f(p_{\rm T}) = \alpha e^{-p_{\rm T}/k}$) for $0.4 < p_{\rm T} < 2$ GeV/$c$ as a function the charged-particle multiplicity density in units of the average normalized multiplicity density for the INEL>0 event class for the neutral pion.

Slope parameter (n) of a power-law fit ($f(p_{\rm T}) = \alpha p_{\rm T}^{-n}$) for $5 < p_{\rm T} < 15$ GeV/$c$ as a function the charged-particle multiplicity density in units of the average normalized multiplicity density for the INEL>0 event class for the eta meson.

Slope parameter (k) of an exponential fit ($f(p_{\rm T}) = \alpha e^{-p_{\rm T}/k}$) for $1.5 < p_{\rm T} < 4$ GeV/$c$ as a function the charged-particle multiplicity density in units of the average normalized multiplicity density for the INEL>0 event class for the eta meson.

$\eta/\pi^{0}$ ratio as a function of $p_{\rm T}$ for pp collisions at $\sqrt{s}$ = 13 TeV in INEL>0 collisions in the V0M multiplicity class 0-100%.

$\eta/\pi^{0}$ ratio as a function of $p_{\rm T}$ for pp collisions at $\sqrt{s}$ = 13 TeV in INEL>0 collisions in the V0M multiplicity class 70-100% (data points not shown in figure 13).

$\eta/\pi^{0}$ ratio as a function of $p_{\rm T}$ for pp collisions at $\sqrt{s}$ = 13 TeV in INEL>0 collisions in the V0M multiplicity class 50-70%.

$\eta/\pi^{0}$ ratio as a function of $p_{\rm T}$ for pp collisions at $\sqrt{s}$ = 13 TeV in INEL>0 collisions in the V0M multiplicity class 30-50% (data points not shown in figure 13).

$\eta/\pi^{0}$ ratio as a function of $p_{\rm T}$ for pp collisions at $\sqrt{s}$ = 13 TeV in INEL>0 collisions in the V0M multiplicity class 20-30% (data points not shown in figure 13).

$\eta/\pi^{0}$ ratio as a function of $p_{\rm T}$ for pp collisions at $\sqrt{s}$ = 13 TeV in INEL>0 collisions in the V0M multiplicity class 10-20% (data points not shown in figure 13).

$\eta/\pi^{0}$ ratio as a function of $p_{\rm T}$ for pp collisions at $\sqrt{s}$ = 13 TeV in INEL>0 collisions in the V0M multiplicity class 5-10% (data points not shown in figure 13).

$\eta/\pi^{0}$ ratio as a function of $p_{\rm T}$ for pp collisions at $\sqrt{s}$ = 13 TeV in INEL>0 collisions in the V0M multiplicity class 1-5% (data points not shown in figure 13).

$\eta/\pi^{0}$ ratio as a function of $p_{\rm T}$ for pp collisions at $\sqrt{s}$ = 13 TeV in INEL>0 collisions in the V0M multiplicity class 0-1% (data points not shown in figure 13).

$\eta/\pi^{0}$ ratio as a function of $p_{\rm T}$ for pp collisions at $\sqrt{s}$ = 13 TeV in INEL>0 collisions in the V0M multiplicity class 0.05-0.1% (data points not shown in figure 13).

$\eta/\pi^{0}$ ratio as a function of $p_{\rm T}$ for pp collisions at $\sqrt{s}$ = 13 TeV in INEL>0 collisions in the V0M multiplicity class 0-0.1%.

$\eta/\pi^{0}$ ratio as a function of $p_{\rm T}$ for pp collisions at $\sqrt{s}$ = 13 TeV in INEL>0 collisions in the V0M multiplicity class 0.01-0.05% (data points not shown in figure 13).

$\eta/\pi^{0}$ ratio as a function of $p_{\rm T}$ for pp collisions at $\sqrt{s}$ = 13 TeV in INEL>0 collisions in the V0M multiplicity class 0-0.01% (data points not shown in figure 13).

Ratio of $\eta/\pi^{0}$ ratio in INEL>0 collisions in the V0M multiplicity class 70-100% to the inclusive $\eta/\pi^{0}$ ratio as a function of $p_{\rm T}$ for pp collisions at $\sqrt{s}$ = 13 TeV (data points not shown in figure 13).

Ratio of $\eta/\pi^{0}$ ratio in INEL>0 collisions in the V0M multiplicity class 50-70% to the inclusive $\eta/\pi^{0}$ ratio as a function of $p_{\rm T}$ for pp collisions at $\sqrt{s}$ = 13 TeV.

Ratio of $\eta/\pi^{0}$ ratio in INEL>0 collisions in the V0M multiplicity class 30-50% to the inclusive $\eta/\pi^{0}$ ratio as a function of $p_{\rm T}$ for pp collisions at $\sqrt{s}$ = 13 TeV (data points not shown in figure 13).

Ratio of $\eta/\pi^{0}$ ratio in INEL>0 collisions in the V0M multiplicity class 20-30% to the inclusive $\eta/\pi^{0}$ ratio as a function of $p_{\rm T}$ for pp collisions at $\sqrt{s}$ = 13 TeV (data points not shown in figure 13).

Ratio of $\eta/\pi^{0}$ ratio in INEL>0 collisions in the V0M multiplicity class 10-20% to the inclusive $\eta/\pi^{0}$ ratio as a function of $p_{\rm T}$ for pp collisions at $\sqrt{s}$ = 13 TeV (data points not shown in figure 13).

Ratio of $\eta/\pi^{0}$ ratio in INEL>0 collisions in the V0M multiplicity class 5-10% to the inclusive $\eta/\pi^{0}$ ratio as a function of $p_{\rm T}$ for pp collisions at $\sqrt{s}$ = 13 TeV (data points not shown in figure 13).

Ratio of $\eta/\pi^{0}$ ratio in INEL>0 collisions in the V0M multiplicity class 1-5% to the inclusive $\eta/\pi^{0}$ ratio as a function of $p_{\rm T}$ for pp collisions at $\sqrt{s}$ = 13 TeV (data points not shown in figure 13).

Ratio of $\eta/\pi^{0}$ ratio in INEL>0 collisions in the V0M multiplicity class 0-1% to the inclusive $\eta/\pi^{0}$ ratio as a function of $p_{\rm T}$ for pp collisions at $\sqrt{s}$ = 13 TeV (data points not shown in figure 13).

Ratio of $\eta/\pi^{0}$ ratio in INEL>0 collisions in the V0M multiplicity class 0.05-0.1% to the inclusive $\eta/\pi^{0}$ ratio as a function of $p_{\rm T}$ for pp collisions at $\sqrt{s}$ = 13 TeV (data points not shown in figure 13).

Ratio of $\eta/\pi^{0}$ ratio in INEL>0 collisions in the V0M multiplicity class 0-0.1% to the inclusive $\eta/\pi^{0}$ ratio as a function of $p_{\rm T}$ for pp collisions at $\sqrt{s}$ = 13 TeV.

Ratio of $\eta/\pi^{0}$ ratio in INEL>0 collisions in the V0M multiplicity class 0.01-0.05% to the inclusive $\eta/\pi^{0}$ ratio as a function of $p_{\rm T}$ for pp collisions at $\sqrt{s}$ = 13 TeV (data points not shown in figure 13).

Ratio of $\eta/\pi^{0}$ ratio in INEL>0 collisions in the V0M multiplicity class 0-0.01% to the inclusive $\eta/\pi^{0}$ ratio as a function of $p_{\rm T}$ for pp collisions at $\sqrt{s}$ = 13 TeV (data points not shown in figure 13).

Mean values of the ratio of the $\eta/\pi^{0}$ ratio in a certain multiplicity interval to the $\eta/\pi^{0}$ ratio in the INEL>0 event class for $1 < p_{\rm{ T}} < 4 \rm{~GeV}/c$ as a function of the normalized mean charged-particle multiplicity.

Mean values of the ratio of the $\eta/\pi^{0}$ ratio in a certain multiplicity interval to the $\eta/\pi^{0}$ ratio in the INEL>0 event class for $5 < p_{\rm{ T}} < 15 \rm{~GeV}/c$ as a function of the normalized mean charged-particle multiplicity.

Correlation functions $P_2^{\rm LS}$ of charged hadrons measured in minimum bias pp collisions at $\sqrt{s}=13\;\text{TeV}$. Charged hadrons are selected in the range $0.2 < p_{\rm T} < 2.0$ GeV/$c$ and with pseudorapidity $|\eta| < 0.8$.

Correlation functions $R_2^{\rm CI}$ of charged hadrons measured in minimum bias pp collisions at $\sqrt{s}=13\;\text{TeV}$. Charged hadrons are selected in the range $0.2 < p_{\rm T} < 2.0$ GeV/$c$ and with pseudorapidity $|\eta| < 0.8$.

Correlation functions $P_2^{\rm CI}$ of charged hadrons measured in minimum bias pp collisions at $\sqrt{s}=13\;\text{TeV}$. Charged hadrons are selected in the range $0.2 < p_{\rm T} < 2.0$ GeV/$c$ and with pseudorapidity $|\eta| < 0.8$.

Projections onto $\Delta\eta$ of $R_2^{\rm CI}$ correlation functions of charged hadrons calculated in the selected $p_{\rm T}$ range in pp collisions at $\sqrt{s}=13\;\text{TeV}$. The $\Delta\eta$ projection is calculated as averages of the two-dimensional correlations in the range $|\Delta\varphi| \leq \pi$ and $|\Delta\eta| \leq 1.6$, respectively. Simulations using PYTHIA8 with the Monash 2013 tune and color reconnection (CR) enabled are added.

Projections onto $\Delta\varphi$ of $R_2^{\rm CI}$ correlation functions of charged hadrons calculated in the selected $p_{\rm T}$ range in pp collisions at $\sqrt{s}=13\;\text{TeV}$. The $\Delta\varphi$ projection is calculated as averages of the two-dimensional correlations in the range $|\Delta\varphi| \leq \pi$ and $|\Delta\eta| \leq 1.6$, respectively. Simulations using PYTHIA8 with the Monash 2013 tune and color reconnection (CR) enabled are added.

Projections onto $\Delta\eta$ of $P_2^{\rm CI}$ correlation functions of charged hadrons calculated in the selected $p_{\rm T}$ range in pp collisions at $\sqrt{s}=13\;\text{TeV}$. The $\Delta\eta$ projection is calculated as averages of the two-dimensional correlations in the range $|\Delta\varphi| \leq \pi$ and $|\Delta\eta| \leq 1.6$, respectively. Simulations using PYTHIA8 with the Monash 2013 tune and color reconnection (CR) enabled are added.

Projections onto $\Delta\varphi$ of $P_2^{\rm CI}$ correlation functions of charged hadrons calculated in the selected $p_{\rm T}$ range in pp collisions at $\sqrt{s}=13\;\text{TeV}$. The $\Delta\varphi$ projection is calculated as averages of the two-dimensional correlations in the range $|\Delta\varphi| \leq \pi$ and $|\Delta\eta| \leq 1.6$, respectively. Simulations using PYTHIA8 with the Monash 2013 tune and color reconnection (CR) enabled are added.

Correlation functions $R_2^{\rm CD}$ of charged hadrons measured in minimum bias pp collisions at $\sqrt{s}=13\;\text{TeV}$. Charged hadrons are selected in the range $0.2 < p_{\rm T} < 2.0$ GeV/$c$ and with pseudorapidity $|\eta| < 0.8$.

Correlation functions $P_2^{\rm CD}$ of charged hadrons measured in minimum bias pp collisions at $\sqrt{s}=13\;\text{TeV}$. Charged hadrons are selected in the range $0.2 < p_{\rm T} < 2.0$ GeV/$c$ and with pseudorapidity $|\eta| < 0.8$.

Projections onto $\Delta\eta$ of $R_2^{\rm CD}$ correlation functions of charged hadrons calculated in the selected $p_{\rm T}$ range in pp collisions at $\sqrt{s}=13\;\text{TeV}$. The $\Delta\eta$ projection is calculated as averages of the two-dimensional correlations in the range $|\Delta\varphi| \leq \pi$ and $|\Delta\eta| \leq 1.6$, respectively. Simulations using PYTHIA8 with the Monash 2013 tune and color reconnection (CR) enabled are added.

Projections onto $\Delta\varphi$ of $R_2^{\rm CD}$ correlation functions of charged hadrons calculated in the selected $p_{\rm T}$ range in pp collisions at $\sqrt{s}=13\;\text{TeV}$. The $\Delta\varphi$ projection is calculated as averages of the two-dimensional correlations in the range $|\Delta\varphi| \leq \pi$ and $|\Delta\eta| \leq 1.6$, respectively. Simulations using PYTHIA8 with the Monash 2013 tune and color reconnection (CR) enabled are added.

Projections onto $\Delta\eta$ of $P_2^{\rm CD}$ correlation functions of charged hadrons calculated in the selected $p_{\rm T}$ range in pp collisions at $\sqrt{s}=13\;\text{TeV}$. The $\Delta\eta$ projection is calculated as averages of the two-dimensional correlations in the range $|\Delta\varphi| \leq \pi$ and $|\Delta\eta| \leq 1.6$, respectively. Simulations using PYTHIA8 with the Monash 2013 tune and color reconnection (CR) enabled are added.

Projections onto $\Delta\varphi$ of $P_2^{\rm CD}$ correlation functions of charged hadrons calculated in the selected $p_{\rm T}$ range in pp collisions at $\sqrt{s}=13\;\text{TeV}$. The $\Delta\varphi$ projection is calculated as averages of the two-dimensional correlations in the range $|\Delta\varphi| \leq \pi$ and $|\Delta\eta| \leq 1.6$, respectively. Simulations using PYTHIA8 with the Monash 2013 tune and color reconnection (CR) enabled are added.

Width of $R_2^{\rm CI}(\Delta\eta)$ and $P_2^{\rm CI}(\Delta\eta)$ correlation functions along $\Delta\eta$ measured within $|\Delta\varphi| < \pi$ in pp collisions and compared with results from $p-Pb$ and $Pb-Pb$ collisions as a function of $\langle \mathrm{d}N_\mathrm{ch}/\mathrm{d}\eta \rangle_{|\eta|<0.5}$.

Width of $R_2^{\rm CD}(\Delta\eta)$ and $P_2^{\rm CD}(\Delta\eta)$ correlation functions along $\Delta\eta$ measured within $|\Delta\varphi| < \pi$ in pp collisions and compared with results from p$-$Pb and Pb$-$Pb collisions as a function of $\langle \mathrm{d}N_\mathrm{ch}/\mathrm{d}\eta \rangle_{|\eta|<0.5}$.

Width of $R_2^{\rm CD}(\Delta\varphi)$ and $P_2^{\rm CD}(\Delta\varphi)$ correlation functions along $\Delta\varphi$ measured within $|\Delta\eta| < 1.6$ in pp collisions and compared with results from p$-$Pb and Pb$-$Pb collisions as a function of $\langle \mathrm{d}N_\mathrm{ch}/\mathrm{d}\eta \rangle_{|\eta|<0.5}$.

Balance function of charged hadrons, in the selected $p_{\rm T}$ range, obtained using ALICE data in pp collisions at $\sqrt{s}=13\;\text{TeV}$.

Integral of $B$ results for charged hadrons in pp collisions at $\sqrt{s}=13\;\text{TeV}$ compared with the integral of $B$ for identified hadron ($\pi \pi$, $\pi$K, and $\pi$p) pairs in Pb$-$Pb collisions at $\sqrt{s_{\rm NN}}=2.76\;\text{TeV}$.

HF Fraction at $1.2 < |\eta| < 2.0$ Top Muon Rich Data

HF Fraction at $1.2 < |\eta| < 2.0$ Bottom Ratio

Charged Hadron $v_2$; $1.2 < |\eta| < 2.0$

Charged Hadron; $1.2 < |\eta| < 2.0$

HF -> $\mu$; $1.2 < |\eta| < 2.0$

HF -> $e$; $|\eta| < 0.35$

JPSI v2 in Au+Au collisions as a function of pT (GeV/c) for 10%--40% centrality with pT binned by [0, 2], [2, 5], and [5, 10] GeV/c.

Isolated-photon cross section ratio of pp collisions at $\sqrt{s}=13~\mathrm{TeV}$ over $\sqrt{s}=5.02~\mathrm{TeV}$ in data (left column) and NLO pQCD calculation from JETPHOX (right column) for $R=0.4$. Data statistical and systematic uncertainties are provided. The theory scale and PDF uncertainties are provided. The data normalisation uncertainty is provided in the paper.

Isolated-photon cross section ratio of pp collisions at $\sqrt{s}=7~\mathrm{TeV}$ over $\sqrt{s}=5.02~\mathrm{TeV}$ in data (left column) and NLO pQCD calculation from JETPHOX (right column) for $R=0.4$. Data statistical and systematic uncertainties are provided. The theory scale and PDF uncertainties are provided. The data normalisation uncertainty is provided in the paper.

Nuclear modification factor ($R_{\rm AA}$) of isolated photons for Pb$-$Pb and pp collisions at $\sqrt{s_{\mathrm{NN}}}=5.02~\mathrm{TeV}$. $R_{\rm AA}$ obtained with a cone size value of $R=0.2$. Each column for each Pb$-$Pb collisions centrality class. The last column for the NLO pQCD JETPHOX calculations, Pb$-$Pb (nPDF) over pp (PDF). Data statistical and systematic uncertainties are provided. The theory scale and PDF uncertainties are provided.

Nuclear modification factor ($R_{\rm AA}$) of isolated photons for Pb$-$Pb and pp collisions at $\sqrt{s_{\mathrm{NN}}}=5.02~\mathrm{TeV}$. $R_{\rm AA}$ obtained with a cone size value of $R=0.4$. Each column for each Pb$-$Pb collisions centrality class. The last column for the NLO pQCD JETPHOX calculations, Pb$-$Pb (nPDF) over pp (PDF). Data statistical and systematic uncertainties are provided. The theory scale and PDF uncertainties are provided.

Ratios $R_{\rm CSP}$ and $R_{\rm CSC}$ of isolated photons cross sections in data for Pb$-$Pb collisions at $\sqrt{s_{\mathrm{NN}}}=5.02~\mathrm{TeV}$. $R_{\rm CSP}$, ratio central over semi-peripheral collisions 50$-$70%. $R_{\rm CSC}$, ratio central over semi-central collisions 30$-$50%. Ratios obtained with a cone size value of $R=0.2$. Data statistical and systematic uncertainties are provided.

Ratios $R_{\rm CSP}$ and $R_{\rm CSC}$ of isolated photons cross sections in data for Pb$-$Pb collisions at $\sqrt{s_{\mathrm{NN}}}=5.02~\mathrm{TeV}$. $R_{\rm CSP}$, ratio central over semi-peripheral collisions 50$-$70%. $R_{\rm CSC}$, ratio central over semi-central collisions 30$-$50%. Ratios obtained with a cone size value of $R=0.4$. Data statistical and systematic uncertainties are provided.

Self-normalised ZN energy as a function of the self-normalised charged-particle-multiplicity in pp collisions at $\sqrt{s}$ = 13 TeV in the high-ZN-energy classification (V0M+SPDClusters event classes).

Self-normalised ZN energy as a function of the self-normalised charged-particle-multiplicity in pp collisions at $\sqrt{s}$ = 13 TeV in the low-ZN-energy classification (V0M+SPDClusters event classes).

$K^{0}_{S}$ transverse momentum spectrum in pp collisions at $\sqrt{s}$ = 13 TeV in the high-ZN-energy classification (V0M+SPDClusters event classes).

$\Lambda+\overline{\Lambda}$ transverse momentum spectrum in pp collisions at $\sqrt{s}$ = 13 TeV in the high-ZN-energy classification (V0M+SPDClusters event classes).

$\Xi^{-}+\overline{\Xi}^{+}$ transverse momentum spectrum in pp collisions at $\sqrt{s}$ = 13 TeV in the high-ZN-energy classification (V0M+SPDClusters event classes).

$K^{0}_{S}$ transverse momentum yield ratios to class III in the high-ZN-energy classification (V0M+SPDClusters event classes) in pp collisions at $\sqrt{s}$ = 13 TeV.

$\Lambda+\overline{\Lambda}$ transverse momentum yield ratios to class III in the high-ZN-energy classification (V0M+SPDClusters event classes) in pp collisions at $\sqrt{s}$ = 13 TeV.

$\Xi^{-}+\overline{\Xi}^{+}$ transverse momentum yield ratios to class III in the high-ZN-energy classification (V0M+SPDClusters event classes) in pp collisions at $\sqrt{s}$ = 13 TeV.

$K^{0}_{S}$ transverse momentum spectrum in pp collisions at $\sqrt{s}$ = 13 TeV in the high-local-multiplicity classification (V0M+SPDClusters event classes).

$\Lambda+\overline{\Lambda}$ transverse momentum spectrum in pp collisions at $\sqrt{s}$ = 13 TeV in the high-local-multiplicity classification (V0M+SPDClusters event classes).

$\Xi^{-}+\overline{\Xi}^{+}$ transverse momentum spectrum in pp collisions at $\sqrt{s}$ = 13 TeV in the high-local-multiplicity classification (V0M+SPDClusters event classes).

$K^{0}_{S}$ transverse momentum yield ratios to class IV in the high-local-multiplicity classification (V0M+SPDClusters event classes) in pp collisions at $\sqrt{s}$ = 13 TeV.

$\Lambda+\overline{\Lambda}$ transverse momentum yield ratios to class IV in the high-local-multiplicity classification (V0M+SPDClusters event classes) in pp collisions at $\sqrt{s}$ = 13 TeV.

$\Xi^{-}+\overline{\Xi}^{+}$ transverse momentum yield ratios to class IV in the high-local-multiplicity classification (V0M+SPDClusters event classes) in pp collisions at $\sqrt{s}$ = 13 TeV.

Self-normalised yield ratios of $K^{0}_{S}$, $\Lambda+\overline{\Lambda}$, and $\Xi^{-}+\overline{\Xi}^{+}$ as a function of the self-normalised charged-particle multiplicity in pp collisions at $\sqrt{s}$ = 13 TeV in the standalone classification (V0M event classes).

Self-normalised yield ratios of $K^{0}_{S}$, $\Lambda+\overline{\Lambda}$, and $\Xi^{-}+\overline{\Xi}^{+}$ as a function of the self-normalised $ZN$ signal in pp collisions at $\sqrt{s}$ = 13 TeV in the standalone classification (V0M event classes).

Self-normalised yield ratios of $K^{0}_{S}$, $\Lambda+\overline{\Lambda}$, and $\Xi^{-}+\overline{\Xi}^{+}$ as a function of the self-normalised charged-particle multiplicity in pp collisions at $\sqrt{s}$ = 13 TeV in the high-local-multiplicity classification (V0M+SPDClusters event classes).

Self-normalised yield ratios of $K^{0}_{S}$, $\Lambda+\overline{\Lambda}$, and $\Xi^{-}+\overline{\Xi}^{+}$ as a function of the self-normalised charged-particle multiplicity in pp collisions at $\sqrt{s}$ = 13 TeV in the low-local-multiplicity classification (V0M+SPDClusters event classes).

Self-normalised yield ratios of $K^{0}_{S}$, $\Lambda+\overline{\Lambda}$, and $\Xi^{-}+\overline{\Xi}^{+}$ as a function of the self-normalised charged-particle multiplicity in pp collisions at $\sqrt{s}$ = 13 TeV in the high-ZN-energy classification (V0M+SPDClusters event classes).

Self-normalised yield ratios of $K^{0}_{S}$, $\Lambda+\overline{\Lambda}$, and $\Xi^{-}+\overline{\Xi}^{+}$ as a function of the self-normalised charged-particle multiplicity in pp collisions at $\sqrt{s}$ = 13 TeV in the low-ZN-energy classification (V0M+SPDClusters event classes).

Self-normalised yield ratios of $K^{0}_{S}$, $\Lambda+\overline{\Lambda}$, and $\Xi^{-}+\overline{\Xi}^{+}$ as a function of the self-normalised $ZN$ signal in pp collisions at $\sqrt{s}$ = 13 TeV in the high-local-multiplicity classification (V0M+SPDClusters event classes).

Self-normalised yield ratios of $K^{0}_{S}$, $\Lambda+\overline{\Lambda}$, and $\Xi^{-}+\overline{\Xi}^{+}$ as a function of the self-normalised $ZN$ signal in pp collisions at $\sqrt{s}$ = 13 TeV in the low-local-multiplicity classification (V0M+SPDClusters event classes).

Self-normalised yield ratios of $K^{0}_{S}$, $\Lambda+\overline{\Lambda}$, and $\Xi^{-}+\overline{\Xi}^{+}$ as a function of the self-normalised $ZN$ signal in pp collisions at $\sqrt{s}$ = 13 TeV in the high-ZN-energy classification (V0M+SPDClusters event classes).

Self-normalised yield ratios of $K^{0}_{S}$, $\Lambda+\overline{\Lambda}$, and $\Xi^{-}+\overline{\Xi}^{+}$ as a function of the self-normalised $ZN$ signal in pp collisions at $\sqrt{s}$ = 13 TeV in the low-ZN-energy classification (V0M+SPDClusters event classes).

Near side asociated charged particle yield per trigger in the range $0.15<p_T<1$ GeV/$c$ for prompt J/$\psi$, as a function of $p_T$, using the HM event samples.

Near side asociated charged particle yield per trigger in the range $0.15<p_T<1$ GeV/$c$ for non-prompt J/$\psi$, as a function of $p_T$, using the MB and EMCal event samples.

Near side asociated charged particle yield per trigger in the range $0.15<p_T<1$ GeV/$c$ for non-prompt J/$\psi$, as a function of $p_T$, using the HM event samples.

Near side asociated charged particle yield per trigger in the range $1<p_T<3$ GeV/$c$ for inclusive J/$\psi$, as a function of $p_T$, using the MB and EMCal event samples.

Near side asociated charged particle yield per trigger in the range $1<p_T<3$ GeV/$c$ for inclusive J/$\psi$, as a function of $p_T$, using the HM event samples.

Near side asociated charged particle yield per trigger in the range $1<p_T<3$ GeV/$c$ for prompt J/$\psi$, as a function of $p_T$, using the MB and EMCal event samples.

Near side asociated charged particle yield per trigger in the range $1<p_T<3$ GeV/$c$ for prompt J/$\psi$, as a function of $p_T$, using the HM event samples.

Near side asociated charged particle yield per trigger in the range $1<p_T<3$ GeV/$c$ for non prompt J/$\psi$, as a function of $p_T$, using the MB and EMCal event samples.

Near side asociated charged particle yield per trigger in the range $1<p_T<3$ GeV/$c$ for non prompt J/$\psi$, as a function of $p_T$, using the HM event samples.

Near side asociated charged particle yield per trigger in the range $3<p_T<10$ GeV/$c$ for inclusive J/$\psi$, as a function of $p_T$, using the MB and EMCal event samples.

Near side asociated charged particle yield per trigger in the range $3<p_T<10$ GeV/$c$ for inclusive J/$\psi$, as a function of $p_T$, using the HM event samples.

Near side asociated charged particle yield per trigger in the range $3<p_T<10$ GeV/$c$ for prompt J/$\psi$, as a function of $p_T$, using the MB and EMCal event samples.

Near side asociated charged particle yield per trigger in the range $3<p_T<10$ GeV/$c$ for prompt J/$\psi$, as a function of $p_T$, using the HM event samples.

Near side asociated charged particle yield per trigger in the range $3<p_T<10$ GeV/$c$ for non prompt J/$\psi$, as a function of $p_T$, using the MB and EMCal event samples.

Near side asociated charged particle yield per trigger in the range $3<p_T<10$ GeV/$c$ for non prompt J/$\psi$, as a function of $p_T$, using the HM event samples.

Pythia calculations for near side asociated charged particle yield per trigger in the range $0.15<p_T<1$ GeV/$c$ for inclusive J/$\psi$, as a function of $p_T$, using the MB and EMCal event samples.

Pythia calculations for near side asociated charged particle yield per trigger in the range $0.15<p_T<1$ GeV/$c$ for inclusive J/$\psi$, as a function of $p_T$, using the HM event samples.

Pythia calculations for near side asociated charged particle yield per trigger in the range $0.15<p_T<1$ GeV/$c$ for prompt J/$\psi$, as a function of $p_T$, using the MB and EMCal event samples.

Pythia calculations for near side asociated charged particle yield per trigger in the range $0.15<p_T<1$ GeV/$c$ for prompt J/$\psi$, as a function of $p_T$, using the HM event samples.

Pythia calculations for near side asociated charged particle yield per trigger in the range $0.15<p_T<1$ GeV/$c$ for non prompt J/$\psi$, as a function of $p_T$, using the MB and EMCal event samples.

Pythia calculations for near side asociated charged particle yield per trigger in the range $0.15<p_T<1$ GeV/$c$ for non prompt J/$\psi$, as a function of $p_T$, using the HM event samples.

Pythia calculations for near side asociated charged particle yield per trigger in the range $1<p_T<3$ GeV/$c$ for inclusive J/$\psi$, as a function of $p_T$, using the MB and EMCal event samples.

Pythia calculations for near side asociated charged particle yield per trigger in the range $1<p_T<3$ GeV/$c$ for inclusive J/$\psi$, as a function of $p_T$, using the HM event samples.

Pythia calculations for near side asociated charged particle yield per trigger in the range $1<p_T<3$ GeV/$c$ for prompt J/$\psi$, as a function of $p_T$, using the MB and EMCal event samples.

Pythia calculations for near side asociated charged particle yield per trigger in the range $1<p_T<3$ GeV/$c$ for prompt J/$\psi$, as a function of $p_T$, using the HM event samples.

Pythia calculations for near side asociated charged particle yield per trigger in the range $1<p_T<3$ GeV/$c$ for non prompt J/$\psi$, as a function of $p_T$, using the MB and EMCal event samples.

Pythia calculations for near side asociated charged particle yield per trigger in the range $1<p_T<3$ GeV/$c$ for non prompt J/$\psi$, as a function of $p_T$, using the HM event samples.

Pythia calculations for near side asociated charged particle yield per trigger in the range $3<p_T<10$ GeV/$c$ for inclusive J/$\psi$, as a function of $p_T$, using the MB and EMCal event samples.

Pythia calculations for near side asociated charged particle yield per trigger in the range $3<p_T<10$ GeV/$c$ for inclusive J/$\psi$, as a function of $p_T$, using the HM event samples.

Pythia calculations for near side asociated charged particle yield per trigger in the range $3<p_T<10$ GeV/$c$ for prompt J/$\psi$, as a function of $p_T$, using the MB and EMCal event samples.

Pythia calculations for near side asociated charged particle yield per trigger in the range $3<p_T<10$ GeV/$c$ for prompt J/$\psi$, as a function of $p_T$, using the HM event samples.

Pythia calculations for near side asociated charged particle yield per trigger in the range $3<p_T<10$ GeV/$c$ for non prompt J/$\psi$, as a function of $p_T$, using the MB and EMCal event samples.

Pythia calculations for near side asociated charged particle yield per trigger in the range $3<p_T<10$ GeV/$c$ for non prompt J/$\psi$, as a function of $p_T$, using the HM event samples.

Away side asociated charged particle yield per trigger in the range $0.15<p_T<1$ GeV/$c$ for inclusive J/$\psi$, as a function of $p_T$, using the MB and EMCal event samples.

Away side asociated charged particle yield per trigger in the range $0.15<p_T<1$ GeV/$c$ for inclusive J/$\psi$, as a function of $p_T$, using the HM event samples.

Away side asociated charged particle yield per trigger in the range $0.15<p_T<1$ GeV/$c$ for prompt J/$\psi$, as a function of $p_T$, using the MB and EMCal event samples.

Away side asociated charged particle yield per trigger in the range $0.15<p_T<1$ GeV/$c$ for prompt J/$\psi$, as a function of $p_T$, using the HM event samples.

Away side asociated charged particle yield per trigger in the range $0.15<p_T<1$ GeV/$c$ for non prompt J/$\psi$, as a function of $p_T$, using the MB and EMCal event samples.

Away side asociated charged particle yield per trigger in the range $0.15<p_T<1$ GeV/$c$ for non prompt J/$\psi$, as a function of $p_T$, using the HM event samples.

Away side asociated charged particle yield per trigger in the range $1<p_T<3$ GeV/$c$ for inclusive J/$\psi$, as a function of $p_T$, using the MB and EMCal event samples.

Away side asociated charged particle yield per trigger in the range $1<p_T<3$ GeV/$c$ for inclusive J/$\psi$, as a function of $p_T$, using the HM event samples.

Away side asociated charged particle yield per trigger in the range $1<p_T<3$ GeV/$c$ for prompt J/$\psi$, as a function of $p_T$, using the MB and EMCal event samples.

Away side asociated charged particle yield per trigger in the range $1<p_T<3$ GeV/$c$ for prompt J/$\psi$, as a function of $p_T$, using the HM event samples.

Away side asociated charged particle yield per trigger in the range $1<p_T<3$ GeV/$c$ for non prompt J/$\psi$, as a function of $p_T$, using the MB and EMCal event samples.

Away side asociated charged particle yield per trigger in the range $1<p_T<3$ GeV/$c$ for non prompt J/$\psi$, as a function of $p_T$, using the HM event samples.

Away side asociated charged particle yield per trigger in the range $3<p_T<10$ GeV/$c$ for inclusive J/$\psi$, as a function of $p_T$, using the MB and EMCal event samples.

Away side asociated charged particle yield per trigger in the range $3<p_T<10$ GeV/$c$ for inclusive J/$\psi$, as a function of $p_T$, using the HM event samples.

Away side asociated charged particle yield per trigger in the range $3<p_T<10$ GeV/$c$ for prompt J/$\psi$, as a function of $p_T$, using the MB and EMCal event samples.

Away side asociated charged particle yield per trigger in the range $3<p_T<10$ GeV/$c$ for prompt J/$\psi$, as a function of $p_T$, using the HM event samples.

Away side asociated charged particle yield per trigger in the range $3<p_T<10$ GeV/$c$ for non prompt J/$\psi$, as a function of $p_T$, using the MB and EMCal event samples.

Away side asociated charged particle yield per trigger in the range $3<p_T<10$ GeV/$c$ for non prompt J/$\psi$, as a function of $p_T$, using the HM event samples.

Pythia calculations for away side asociated charged particle yield per trigger in the range $0.15<p_T<1$ GeV/$c$ for inclusive J/$\psi$, as a function of $p_T$, using the MB and EMCal event samples.

Pythia calculations for away side asociated charged particle yield per trigger in the range $0.15<p_T<1$ GeV/$c$ for inclusive J/$\psi$, as a function of $p_T$, using the HM event samples.

Pythia calculations for away side asociated charged particle yield per trigger in the range $0.15<p_T<1$ GeV/$c$ for prompt J/$\psi$, as a function of $p_T$, using the MB and EMCal event samples.

Pythia calculations for away side asociated charged particle yield per trigger in the range $0.15<p_T<1$ GeV/$c$ for prompt J/$\psi$, as a function of $p_T$, using the HM event samples.

Pythia calculations for away side asociated charged particle yield per trigger in the range $0.15<p_T<1$ GeV/$c$ for non prompt J/$\psi$, as a function of $p_T$, using the MB and EMCal event samples.

Pythia calculations for away side asociated charged particle yield per trigger in the range $0.15<p_T<1$ GeV/$c$ for non prompt J/$\psi$, as a function of $p_T$, using the HM event samples.

Pythia calculations for away side asociated charged particle yield per trigger in the range $1<p_T<3$ GeV/$c$ for inclusive J/$\psi$, as a function of $p_T$, using the MB and EMCal event samples.

Pythia calculations for away side asociated charged particle yield per trigger in the range $1<p_T<3$ GeV/$c$ for inclusive J/$\psi$, as a function of $p_T$, using the HM event samples.

Pythia calculations for away side asociated charged particle yield per trigger in the range $1<p_T<3$ GeV/$c$ for prompt J/$\psi$, as a function of $p_T$, using the MB and EMCal event samples.

Pythia calculations for away side asociated charged particle yield per trigger in the range $1<p_T<3$ GeV/$c$ for prompt J/$\psi$, as a function of $p_T$, using the HM event samples.

Pythia calculations for away side asociated charged particle yield per trigger in the range $1<p_T<3$ GeV/$c$ for non prompt J/$\psi$, as a function of $p_T$, using the MB and EMCal event samples.

Pythia calculations for away side asociated charged particle yield per trigger in the range $1<p_T<3$ GeV/$c$ for non prompt J/$\psi$, as a function of $p_T$, using the HM event samples.

Pythia calculations for away side asociated charged particle yield per trigger in the range $3<p_T<10$ GeV/$c$ for inclusive J/$\psi$, as a function of $p_T$, using the MB and EMCal event samples.

Pythia calculations for away side asociated charged particle yield per trigger in the range $3<p_T<10$ GeV/$c$ for inclusive J/$\psi$, as a function of $p_T$, using the HM event samples.

Pythia calculations for away side asociated charged particle yield per trigger in the range $3<p_T<10$ GeV/$c$ for prompt J/$\psi$, as a function of $p_T$, using the MB and EMCal event samples.

Pythia calculations for away side asociated charged particle yield per trigger in the range $3<p_T<10$ GeV/$c$ for prompt J/$\psi$, as a function of $p_T$, using the HM event samples.

Pythia calculations for away side asociated charged particle yield per trigger in the range $3<p_T<10$ GeV/$c$ for non prompt J/$\psi$, as a function of $p_T$, using the MB and EMCal event samples.

Pythia calculations for away side asociated charged particle yield per trigger in the range $3<p_T<10$ GeV/$c$ for non prompt J/$\psi$, as a function of $p_T$, using the HM event samples.

Multiplicity dependence of the near-side width $\sigma$ in pp collisions at $\sqrt{s_{\rm NN}} = 13$ TeV. Obtained in transverse momentum intervals $3.0 < p_\mathrm{T, trig} < 4.0$ GeV/$c$ and $1.0 < p_\mathrm{T, assoc} < 2.0$ GeV/$c$. The mulitplicity is estimated with midrapidity multiplicity estimator ($|\eta|<1.0,\,p_\mathrm{T}>0.2$ GeV/$c$).

Multiplicity dependence of the near-side width $\sigma$ in pp collisions at $\sqrt{s_{\rm NN}} = 13$ TeV. Obtained in transverse momentum intervals $3.0 < p_\mathrm{T, trig} < 4.0$ GeV/$c$ and $2.0 < p_\mathrm{T, assoc} < 3.0$ GeV/$c$. The mulitplicity is estimated with midrapidity multiplicity estimator ($|\eta|<1.0,\,p_\mathrm{T}>0.2$ GeV/$c$).

Multiplicity dependence of the near-side width $\sigma$ in pp collisions at $\sqrt{s_{\rm NN}} = 13$ TeV. Obtained in transverse momentum intervals $3.0 < p_\mathrm{T, assoc} < p_\mathrm{T, trig} < 4.0$ GeV/$c$. The multiplicity is estimated with midrapidity multiplicity estimator ($|\eta|<1.0,\,p_\mathrm{T}>0.2$ GeV/$c$).

Multiplicity dependence of the near-side width $\sigma$ in pp collisions at $\sqrt{s_{\rm NN}} = 13$ TeV. Obtained in transverse momentum intervals $4.0 < p_\mathrm{T, trig} < 8.0$ GeV/$c$ and $1.0 < p_\mathrm{T, assoc} < 2.0$ GeV/$c$. The mulitplicity is estimated with midrapidity multiplicity estimator ($|\eta|<1.0,\,p_\mathrm{T}>0.2$ GeV/$c$).

Multiplicity dependence of the near-side width $\sigma$ in pp collisions at $\sqrt{s_{\rm NN}} = 13$ TeV. Obtained in transverse momentum intervals $4.0 < p_\mathrm{T, trig} < 8.0$ GeV/$c$ and $2.0 < p_\mathrm{T, assoc} < 3.0$ GeV/$c$. The mulitplicity is estimated with midrapidity multiplicity estimator ($|\eta|<1.0,\,p_\mathrm{T}>0.2$ GeV/$c$).

Multiplicity dependence of the near-side width $\sigma$ in pp collisions at $\sqrt{s_{\rm NN}} = 13$ TeV. Obtained in transverse momentum intervals $4.0 < p_\mathrm{T, trig} < 8.0$ GeV/$c$ and $3.0 < p_\mathrm{T, assoc} < 4.0$ GeV/$c$. The mulitplicity is estimated with midrapidity multiplicity estimator ($|\eta|<1.0,\,p_\mathrm{T}>0.2$ GeV/$c$).

Multiplicity dependence of the near-side width $\sigma$ in pp collisions at $\sqrt{s_{\rm NN}} = 13$ TeV. Obtained in transverse momentum intervals $4.0 < p_\mathrm{T, assoc} < p_\mathrm{T, trig} < 8.0$ GeV/$c$. The multiplicity is estimated with midrapidity multiplicity estimator ($|\eta|<1.0,\,p_\mathrm{T}>0.2$ GeV/$c$).

Multiplicity dependence of the near-side width $\sigma$ in pp collisions at $\sqrt{s_{\rm NN}} = 13$ TeV. Obtained in transverse momentum intervals $1.0 < p_\mathrm{T, assoc} < p_\mathrm{T, trig} < 2.0$ GeV/$c$. The multiplicity is estimated with forward multiplicity estimator.

Multiplicity dependence of the near-side width $\sigma$ in pp collisions at $\sqrt{s_{\rm NN}} = 13$ TeV. Obtained in transverse momentum intervals $2.0 < p_\mathrm{T, trig} < 3.0$ GeV/$c$ and $1.0 < p_\mathrm{T, assoc} < 2.0$ GeV/$c$. The mulitplicity is estimated with forward multiplicity estimator.

Multiplicity dependence of the near-side width $\sigma$ in pp collisions at $\sqrt{s_{\rm NN}} = 13$ TeV. Obtained in transverse momentum intervals $2.0 < p_\mathrm{T, assoc} < p_\mathrm{T, trig} < 3.0$ GeV/$c$. The multiplicity is estimated with forward multiplicity estimator.

Multiplicity dependence of the near-side width $\sigma$ in pp collisions at $\sqrt{s_{\rm NN}} = 13$ TeV. Obtained in transverse momentum intervals $3.0 < p_\mathrm{T, trig} < 4.0$ GeV/$c$ and $1.0 < p_\mathrm{T, assoc} < 2.0$ GeV/$c$. The mulitplicity is estimated with forward multiplicity estimator.

Multiplicity dependence of the near-side width $\sigma$ in pp collisions at $\sqrt{s_{\rm NN}} = 13$ TeV. Obtained in transverse momentum intervals $3.0 < p_\mathrm{T, trig} < 4.0$ GeV/$c$ and $2.0 < p_\mathrm{T, assoc} < 3.0$ GeV/$c$. The mulitplicity is estimated with forward multiplicity estimator.

Multiplicity dependence of the near-side width $\sigma$ in pp collisions at $\sqrt{s_{\rm NN}} = 13$ TeV. Obtained in transverse momentum intervals $3.0 < p_\mathrm{T, assoc} < p_\mathrm{T, trig} < 4.0$ GeV/$c$. The multiplicity is estimated with forward multiplicity estimator.

Multiplicity dependence of the near-side width $\sigma$ in pp collisions at $\sqrt{s_{\rm NN}} = 13$ TeV. Obtained in transverse momentum intervals $4.0 < p_\mathrm{T, trig} < 8.0$ GeV/$c$ and $1.0 < p_\mathrm{T, assoc} < 2.0$ GeV/$c$. The mulitplicity is estimated with forward multiplicity estimator.

Multiplicity dependence of the near-side width $\sigma$ in pp collisions at $\sqrt{s_{\rm NN}} = 13$ TeV. Obtained in transverse momentum intervals $4.0 < p_\mathrm{T, trig} < 8.0$ GeV/$c$ and $2.0 < p_\mathrm{T, assoc} < 3.0$ GeV/$c$. The mulitplicity is estimated with forward multiplicity estimator.

Multiplicity dependence of the near-side width $\sigma$ in pp collisions at $\sqrt{s_{\rm NN}} = 13$ TeV. Obtained in transverse momentum intervals $4.0 < p_\mathrm{T, trig} < 8.0$ GeV/$c$ and $3.0 < p_\mathrm{T, assoc} < 4.0$ GeV/$c$. The mulitplicity is estimated with forward multiplicity estimator.

Multiplicity dependence of the near-side width $\sigma$ in pp collisions at $\sqrt{s_{\rm NN}} = 13$ TeV. Obtained in transverse momentum intervals $4.0 < p_\mathrm{T, assoc} < p_\mathrm{T, trig} < 8.0$ GeV/$c$. The multiplicity is estimated with forward multiplicity estimator.

Multiplicity dependence of the fragmentation yield $Y^{frag}$ in pp collisions at $\sqrt{s_{\rm NN}} = 13$ TeV. Obtained in transverse momentum intervals $1.0 < p_\mathrm{T, assoc} < p_\mathrm{T, trig} < 2.0$ GeV/$c$. The mulitplicity is estimated with midrapidity multiplicity estimator ($|\eta|<1.0,\,p_\mathrm{T}>0.2$ GeV/$c$).

Multiplicity dependence of the fragmentation yield $Y^{frag}$ in pp collisions at $\sqrt{s_{\rm NN}} = 13$ TeV. Obtained in transverse momentum intervals $2.0 < p_\mathrm{T, trig} < 3.0$ GeV/$c$ and $1.0 < p_\mathrm{T, assoc} < 2.0$ GeV/$c$. The mulitplicity is estimated with midrapidity multiplicity estimator ($|\eta|<1.0,\,p_\mathrm{T}>0.2$ GeV/$c$).

Multiplicity dependence of the fragmentation yield $Y^{frag}$ in pp collisions at $\sqrt{s_{\rm NN}} = 13$ TeV. Obtained in transverse momentum intervals $2.0 < p_\mathrm{T, assoc} < p_\mathrm{T, trig} < 3.0$ GeV/$c$. The mulitplicity is estimated with midrapidity multiplicity estimator ($|\eta|<1.0,\,p_\mathrm{T}>0.2$ GeV/$c$).

Multiplicity dependence of the fragmentation yield $Y^{frag}$ in pp collisions at $\sqrt{s_{\rm NN}} = 13$ TeV. Obtained in transverse momentum intervals $3.0 < p_\mathrm{T, trig} < 4.0$ GeV/$c$ and $1.0 < p_\mathrm{T, assoc} < 2.0$ GeV/$c$. The mulitplicity is estimated with midrapidity multiplicity estimator ($|\eta|<1.0,\,p_\mathrm{T}>0.2$ GeV/$c$).

Multiplicity dependence of the fragmentation yield $Y^{frag}$ in pp collisions at $\sqrt{s_{\rm NN}} = 13$ TeV. Obtained in transverse momentum intervals $3.0 < p_\mathrm{T, trig} < 4.0$ GeV/$c$ and $2.0 < p_\mathrm{T, assoc} < 3.0$ GeV/$c$. The mulitplicity is estimated with midrapidity multiplicity estimator ($|\eta|<1.0,\,p_\mathrm{T}>0.2$ GeV/$c$).

Multiplicity dependence of the fragmentation yield $Y^{frag}$ in pp collisions at $\sqrt{s_{\rm NN}} = 13$ TeV. Obtained in transverse momentum intervals $3.0 < p_\mathrm{T, assoc} < p_\mathrm{T, trig} < 4.0$ GeV/$c$. The mulitplicity is estimated with midrapidity multiplicity estimator ($|\eta|<1.0,\,p_\mathrm{T}>0.2$ GeV/$c$).

Multiplicity dependence of the fragmentation yield $Y^{frag}$ in pp collisions at $\sqrt{s_{\rm NN}} = 13$ TeV. Obtained in transverse momentum intervals $4.0 < p_\mathrm{T, trig} < 8.0$ GeV/$c$ and $1.0 < p_\mathrm{T, assoc} < 2.0$ GeV/$c$. The mulitplicity is estimated with midrapidity multiplicity estimator ($|\eta|<1.0,\,p_\mathrm{T}>0.2$ GeV/$c$).

Multiplicity dependence of the fragmentation yield $Y^{frag}$ in pp collisions at $\sqrt{s_{\rm NN}} = 13$ TeV. Obtained in transverse momentum intervals $4.0 < p_\mathrm{T, trig} < 8.0$ GeV/$c$ and $2.0 < p_\mathrm{T, assoc} < 3.0$ GeV/$c$. The mulitplicity is estimated with midrapidity multiplicity estimator ($|\eta|<1.0,\,p_\mathrm{T}>0.2$ GeV/$c$).

Multiplicity dependence of the fragmentation yield $Y^{frag}$ in pp collisions at $\sqrt{s_{\rm NN}} = 13$ TeV. Obtained in transverse momentum intervals $4.0 < p_\mathrm{T, trig} < 8.0$ GeV/$c$ and $3.0 < p_\mathrm{T, assoc} < 4.0$ GeV/$c$. The mulitplicity is estimated with midrapidity multiplicity estimator ($|\eta|<1.0,\,p_\mathrm{T}>0.2$ GeV/$c$).

Multiplicity dependence of the fragmentation yield $Y^{frag}$ in pp collisions at $\sqrt{s_{\rm NN}} = 13$ TeV. Obtained in transverse momentum intervals $4.0 < p_\mathrm{T, assoc} < p_\mathrm{T, trig} < 8.0$ GeV/$c$. The mulitplicity is estimated with midrapidity multiplicity estimator ($|\eta|<1.0,\,p_\mathrm{T}>0.2$ GeV/$c$).

Multiplicity dependence of the fragmentation yield $Y^{frag}$ in pp collisions at $\sqrt{s_{\rm NN}} = 13$ TeV. Obtained in transverse momentum intervals $1.0 < p_\mathrm{T, assoc} < p_\mathrm{T, trig} < 2.0$ GeV/$c$. The mulitplicity is estimated with forward multiplicity estimator.

Multiplicity dependence of the fragmentation yield $Y^{frag}$ in pp collisions at $\sqrt{s_{\rm NN}} = 13$ TeV. Obtained in transverse momentum intervals $2.0 < p_\mathrm{T, trig} < 3.0$ GeV/$c$ and $1.0 < p_\mathrm{T, assoc} < 2.0$ GeV/$c$. The mulitplicity is estimated with forward multiplicity estimator.

Multiplicity dependence of the fragmentation yield $Y^{frag}$ in pp collisions at $\sqrt{s_{\rm NN}} = 13$ TeV. Obtained in transverse momentum intervals $2.0 < p_\mathrm{T, assoc} < p_\mathrm{T, trig} < 3.0$ GeV/$c$. The mulitplicity is estimated with forward multiplicity estimator.

Multiplicity dependence of the fragmentation yield $Y^{frag}$ in pp collisions at $\sqrt{s_{\rm NN}} = 13$ TeV. Obtained in transverse momentum intervals $3.0 < p_\mathrm{T, trig} < 4.0$ GeV/$c$ and $1.0 < p_\mathrm{T, assoc} < 2.0$ GeV/$c$. The mulitplicity is estimated with forward multiplicity estimator.

Multiplicity dependence of the fragmentation yield $Y^{frag}$ in pp collisions at $\sqrt{s_{\rm NN}} = 13$ TeV. Obtained in transverse momentum intervals $3.0 < p_\mathrm{T, trig} < 4.0$ GeV/$c$ and $2.0 < p_\mathrm{T, assoc} < 3.0$ GeV/$c$. The mulitplicity is estimated with forward multiplicity estimator.

Multiplicity dependence of the fragmentation yield $Y^{frag}$ in pp collisions at $\sqrt{s_{\rm NN}} = 13$ TeV. Obtained in transverse momentum intervals $3.0 < p_\mathrm{T, assoc} < p_\mathrm{T, trig} < 4.0$ GeV/$c$. The mulitplicity is estimated with forward multiplicity estimator.

Multiplicity dependence of the fragmentation yield $Y^{frag}$ in pp collisions at $\sqrt{s_{\rm NN}} = 13$ TeV. Obtained in transverse momentum intervals $4.0 < p_\mathrm{T, trig} < 8.0$ GeV/$c$ and $1.0 < p_\mathrm{T, assoc} < 2.0$ GeV/$c$. The mulitplicity is estimated with forward multiplicity estimator.

Multiplicity dependence of the fragmentation yield $Y^{frag}$ in pp collisions at $\sqrt{s_{\rm NN}} = 13$ TeV. Obtained in transverse momentum intervals $4.0 < p_\mathrm{T, trig} < 8.0$ GeV/$c$ and $2.0 < p_\mathrm{T, assoc} < 3.0$ GeV/$c$. The mulitplicity is estimated with forward multiplicity estimator.

Multiplicity dependence of the fragmentation yield $Y^{frag}$ in pp collisions at $\sqrt{s_{\rm NN}} = 13$ TeV. Obtained in transverse momentum intervals $4.0 < p_\mathrm{T, trig} < 8.0$ GeV/$c$ and $3.0 < p_\mathrm{T, assoc} < 4.0$ GeV/$c$. The mulitplicity is estimated with forward multiplicity estimator.

Multiplicity dependence of the fragmentation yield $Y^{frag}$ in pp collisions at $\sqrt{s_{\rm NN}} = 13$ TeV. Obtained in transverse momentum intervals $4.0 < p_\mathrm{T, assoc} < p_\mathrm{T, trig} < 8.0$ GeV/$c$. The mulitplicity is estimated with forward multiplicity estimator.

Self-Normalized $J/\psi$ yields in the same arms before and after dimuon subtraction and the comparison with PYTHIA 8 simulations

Self-normalized $\psi(2S)$ to $J/\psi$ yields ratio as a function of Self-normalized $N_{ch}$ in the same arms

Self-normalized $\psi(2S)$ to $J/\psi$ yields ratio as a function of Self-normalized $N_{ch}$ in opposite arms