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.

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).

Covariance for full matrix measurement all bins from $m(t\bar{t})$ fit

Results for full matrix measurement inclusive from $p_{T}(t)$

Covariance for full matrix measurement inclusive from $p_{T}(t)$

Results for full matrix measurement all bins from $p_{T}(t)$ fit

Covariance for full matrix measurement all bins from $p_{T}(t)$ fit

Results for full matrix measurement differential in $m(t\bar{t})$

Covariance for full matrix measurement differential in $m(t\bar{t})$

Results for full matrix measurement differential in $m(t\bar{t})$ for $|cos(\theta)| < 0.4$

Covariance for full matrix measurement differential in $m(t\bar{t})$ for $|cos(\theta)| < 0.4$

Results for full matrix measurement differential in $p_{T}(t)$

Covariance for full matrix measurement differential in $p_{T}(t)$

Results for $D$ measurement inclusive from $m(t\bar{t})$

Covariance for $D$ measurement inclusive from $m(t\bar{t})$

Results for $D$ measurement inclusive from $p_{T}(t)$

Covariance for $D$ measurement inclusive from $p_{T}(t)$

Results for $\tilde{D}$ measurement inclusive in $m(t\bar{t})$

Covariance for $\tilde{D}$ measurement inclusive in $m(t\bar{t})$

Results for $\tilde{D}$ measurement inclusive in $p_{T}(t)$

Covariance for $\tilde{D}$ measurement inclusive in $p_{T}(t)$

Results for $D$ measurement all bins from from $m(t\bar{t})$ fit

Covariance for $D$ measurement all bins from from $m(t\bar{t})$ fit

Results for $D$ measurement all bins from from $p_{T}(t)$ fit

Covariance for $D$ measurement all bins from from $p_{T}(t)$ fit

Results for $\tilde{D}$ measurement all bins from from $m(t\bar{t})$ fit

Covariance for $\tilde{D}$ measurement all bins from from $m(t\bar{t})$ fit

Results for $\tilde{D}$ measurement all bins from from $p_{T}(t)$ fit

Covariance for $\tilde{D}$ measurement all bins from from $p_{T}(t)$ fit

Results for $D$ measurement differential in $m(t\bar{t})$

Covariance for $D$ measurement differential in $m(t\bar{t})$

Results for $\tilde{D}$ measurement differential in $m(t\bar{t})$

Covariance for $\tilde{D}$ measurement differential in $m(t\bar{t})$

Results for $\tilde{D}$ measurement differential in $m(t\bar{t})$ for $|cos(\theta)| < 0.4$

Covariance for $\tilde{D}$ measurement differential in $m(t\bar{t})$ for $|cos(\theta)| < 0.4$

Results for $D$ measurement differential in $p_{T}(t)$

Covariance for $D$ measurement differential in $p_{T}(t)$

Results for $\tilde{D}$ measurement differential in $p_{T}(t)$

Covariance for $\tilde{D}$ measurement differential in $p_{T}(t)$

The 95% CL upper limits on the branching fraction $\mathcal{B}(\mathrm{H \to SS})$ with $\mathrm{S}\to\tau\tau$ decay, for different LLP masses $m_{\mathrm{S}}$ and proper decay lengths $c\tau_{0}$.

The 95% CL limits on the LLP mass $m_{\mathrm{S}}$ for different proper decay lengths $c\tau_{0}$ assuming a branching fraction of 1% for the $\mathrm{H \to SS}$ decay, and with subsequent $\mathrm{S \to b\overline{b}}$ decay.

The 95% CL limits on the LLP mass $m_{\mathrm{S}}$ for different proper decay lengths $c\tau_{0}$ assuming a branching fraction of 1% for the $\mathrm{H \to SS}$ decay, and with subsequent $\mathrm{S \to d \overline{d}}$ decay.

The 95% CL limits on the dark-sector top quark partner mass $m_{T}$ for different hidden glueball masses $m_{0}$ in the fraternal twin Higgs model.

The 95% CL limits on the dark-sector top quark partner mass $m_{T}$ for different hidden glueball masses $m_{0}$ in the folded SUSY model.

The 95% CL upper limits on $\sigma(\mathrm{pp}\to\mathrm{Z}')\mathcal{B}(\mathrm{Z}'\to\mathrm{Q_{dark}\overline{Q}_{dark}})$ when the dark pion mass $m_{\pi_{\mathrm{dark}}}$ is $20~\mathrm{GeV}$. For the dark sector parton showering and hadronization $N_{f}=7$, $N_{c}=3$, $\Lambda_{\mathrm{dark}}/m_{\pi_{\mathrm{dark}}}=2$, and $m_{\rho_{\mathrm{dark}}}/m_{\pi_{\mathrm{dark}}}=4$ are assumed.

The 95% CL upper limits on $\sigma(\mathrm{pp}\to\mathrm{Z}')\mathcal{B}(\mathrm{Z}'\to\mathrm{Q_{dark}\overline{Q}_{dark}})$ when the dark pion mass $m_{\pi_{\mathrm{dark}}}$ is $15~\mathrm{GeV}$. For the dark sector parton showering and hadronization $N_{f}=7$, $N_{c}=3$, $\Lambda_{\mathrm{dark}}/m_{\pi_{\mathrm{dark}}}=2$, and $m_{\rho_{\mathrm{dark}}}/m_{\pi_{\mathrm{dark}}}=4$ are assumed.

The 95% CL upper limits on $\sigma(\mathrm{pp}\to\mathrm{Z}')\mathcal{B}(\mathrm{Z}'\to\mathrm{Q_{dark}\overline{Q}_{dark}})$ when the dark pion mass $m_{\pi_{\mathrm{dark}}}$ is $10~\mathrm{GeV}$. For the dark sector parton showering and hadronization $N_{f}=7$, $N_{c}=3$, $\Lambda_{\mathrm{dark}}/m_{\pi_{\mathrm{dark}}}=2$, and $m_{\rho_{\mathrm{dark}}}/m_{\pi_{\mathrm{dark}}}=4$ are assumed.

The 95% CL upper limits on the branching fraction for the $125~\mathrm{GeV}$ Higgs boson to decay to a pair of dark quarks $\mathcal{B}(\mathrm{H}\to \mathrm{Q_{dark}\overline{Q}_{dark}})$ when the dark sector $N_{f}=1$, $N_{c}=3$. The dark pion $\pi_{\mathrm{dark}}$ is assumed to decay to SM particles through an offshell Higgs boson. The upper limits are shown for different dark pion masses $m_{\pi_{\mathrm{dark}}}$ and proper decay lengths $c\tau_{\pi_{\mathrm{dark}}}$, with $\Lambda_{\mathrm{dark}}/m_{\pi_{\mathrm{dark}}}=2.5$, $m_{\rho_{\mathrm{dark}}}/m_{\pi_{\mathrm{dark}}}=2.5$.

The 95% CL upper limits on the branching fraction for the $125~\mathrm{GeV}$ Higgs boson to decay to a pair of dark quarks $\mathcal{B}(\mathrm{H}\to \mathrm{Q_{dark}\overline{Q}_{dark}})$ when the dark sector $N_{f}=1$, $N_{c}=3$. The dark pion $\pi_{\mathrm{dark}}$ is assumed to decay to SM particles through an offshell Higgs boson. The upper limits are shown for different dark pion masses $m_{\pi_{\mathrm{dark}}}$ and proper decay lengths $c\tau_{\pi_{\mathrm{dark}}}$, with $\Lambda_{\mathrm{dark}}/m_{\pi_{\mathrm{dark}}}=1$, $m_{\rho_{\mathrm{dark}}}/m_{\pi_{\mathrm{dark}}}=2.5$.

The 95% CL upper limits on the branching fraction for the $125~\mathrm{GeV}$ Higgs boson to decay to a pair of dark quarks $\mathcal{B}(\mathrm{H}\to \mathrm{Q_{dark}\overline{Q}_{dark}})$ when the dark sector $N_{f}=1$, $N_{c}=3$. The dark pion $\pi_{\mathrm{dark}}$ is assumed to decay to SM particles through an offshell Higgs boson. The upper limits are shown for different dark pion masses $m_{\pi_{\mathrm{dark}}}$ and proper decay lengths $c\tau_{\pi_{\mathrm{dark}}}$, with $\Lambda_{\mathrm{dark}}/m_{\pi_{\mathrm{dark}}}=1$, $m_{\rho_{\mathrm{dark}}}/m_{\pi_{\mathrm{dark}}}=1$.

The 95% CL upper limits on the branching fraction for the $125~\mathrm{GeV}$ Higgs boson to decay to a pair of dark quarks $\mathcal{B}(\mathrm{H}\to \mathrm{Q_{dark}\overline{Q}_{dark}})$ when the dark sector $N_{f}=1$, $N_{c}=3$. The dark pion $\pi_{\mathrm{dark}}$ is assumed to decay to a pair of gluons. The upper limits are shown for different dark pion masses $m_{\pi_{\mathrm{dark}}}$ and proper decay lengths $c\tau_{\pi_{\mathrm{dark}}}$, with $\Lambda_{\mathrm{dark}}/m_{\pi_{\mathrm{dark}}}=2.5$, $m_{\rho_{\mathrm{dark}}}/m_{\pi_{\mathrm{dark}}}=2.5$.

The 95% CL upper limits on the branching fraction for the $125~\mathrm{GeV}$ Higgs boson to decay to a pair of dark quarks $\mathcal{B}(\mathrm{H}\to \mathrm{Q_{dark}\overline{Q}_{dark}})$ when the dark sector $N_{f}=1$, $N_{c}=3$. The dark pion $\pi_{\mathrm{dark}}$ is assumed to decay to a pair of gluons. The upper limits are shown for different dark pion masses $m_{\pi_{\mathrm{dark}}}$ and proper decay lengths $c\tau_{\pi_{\mathrm{dark}}}$, with $\Lambda_{\mathrm{dark}}/m_{\pi_{\mathrm{dark}}}=1$, $m_{\rho_{\mathrm{dark}}}/m_{\pi_{\mathrm{dark}}}=2.5$.

The 95% CL upper limits on the branching fraction for the $125~\mathrm{GeV}$ Higgs boson to decay to a pair of dark quarks $\mathcal{B}(\mathrm{H}\to \mathrm{Q_{dark}\overline{Q}_{dark}})$ when the dark sector $N_{f}=1$, $N_{c}=3$. The dark pion $\pi_{\mathrm{dark}}$ is assumed to decay to a pair of gluons. The upper limits are shown for different dark pion masses $m_{\pi_{\mathrm{dark}}}$ and proper decay lengths $c\tau_{\pi_{\mathrm{dark}}}$, with $\Lambda_{\mathrm{dark}}/m_{\pi_{\mathrm{dark}}}=1$, $m_{\rho_{\mathrm{dark}}}/m_{\pi_{\mathrm{dark}}}=1$.

The 95% CL upper limits on the branching fraction for the $125~\mathrm{GeV}$ Higgs boson to decay to a pair of dark quarks $\mathcal{B}(\mathrm{H}\to \mathrm{Q_{dark}\overline{Q}_{dark}})$ when the dark sector $N_{f}=7$, $N_{c}=3$. The dark pion $\pi_{\mathrm{dark}}$ is assumed to decay to a pair of down quarks. The upper limits are shown for different dark pion masses $m_{\pi_{\mathrm{dark}}}$ and proper decay lengths $c\tau_{\pi_{\mathrm{dark}}}$, with $\Lambda_{\mathrm{dark}}/m_{\pi_{\mathrm{dark}}}=2$, $m_{\rho_{\mathrm{dark}}}/m_{\pi_{\mathrm{dark}}}=4$.

The Bs meson $p_{\rm{T}}$-dependent RAA in PpPp. The measurment was carried out inside a fiducial region respecting ($p_{\rm{T}}$<10 & 1.5<|y|<2.4) and ($p_{\rm{T}}$>10 & |y|<2.4).

Distribution of the RF output for events in the 2j1b region. The number of observed events (points) and estimated signal and background events (filled histograms) before the maximum likelihood fit are shown. The vertical bars on the points represent the statistical uncertainty in the data, and the hatched band the total uncertainty in the estimated events before the fit. The lower panels display the ratio of the data to the sum of the estimated events (points) before the fit, with the bands giving the corresponding uncertainties.

The distribution of the Subleading jet $p_{T}$ for events in the 2j2b region. The number of observed events (points) and estimated signal and background events (filled histograms) from the maximum likelihood fit are shown. The vertical bars on the points represent the statistical uncertainty in the data, and the hatched band the total uncertainty in the estimated events after the fit. The lower panels display the ratio of the data to the sum of the estimated events (points) after the fit, with the bands giving the corresponding uncertainties.

Distribution of the Subleading jet $p_{T}$ for events in the 2j2b region. The number of observed events (points) and estimated signal and background events (filled histograms) before the maximum likelihood fit are shown. The vertical bars on the points represent the statistical uncertainty in the data, and the hatched band the total uncertainty in the estimated events before the fit. The lower panels display the ratio of the data to the sum of the estimated events (points) before the fit, with the bands giving the corresponding uncertainties.

Observed and theoretical cross sections. In the observed, the first uncertainty is statistical, the second is the systematic, and the third the luminosity. In the theoretical, the first uncertainty is due to scale variations, the second due to the choice of PDF. The theoretical cross section has been computed at approximate third-order in QCD, with the addition of third-order corrections of soft gluon emission terms, assuming a top quark mass of 172.5 GeV and using the PDF4LHC21 PDF set.

The distribution of the RF output for events in the 1j1b region. The number of observed events (points) and estimated signal and background events (filled histograms) from the maximum likelihood fit are shown. The vertical bars on the points represent the statistical uncertainty in the data, and the hatched band the total uncertainty in the estimated events after the fit. The lower panels display the ratio of the data to the sum of the estimated events (points) after the fit, with the bands giving the corresponding uncertainties.

The number of estimated signal and background events after the fit in the 1j1b, 2j1b, and 2j2b regions compared to the observed number of events. The total uncertainties in the estimated events after the fit are given.

The distribution of the RF output for events in the 2j1b region. The number of observed events (points) and estimated signal and background events (filled histograms) from the maximum likelihood fit are shown. The vertical bars on the points represent the statistical uncertainty in the data, and the hatched band the total uncertainty in the estimated events after the fit. The lower panels display the ratio of the data to the sum of the estimated events (points) after the fit, with the bands giving the corresponding uncertainties.

Normalised fiducial differential tW production cross section as a function of the $Leading$ $lepton$ $p_{T}$. The horizontal bars on the points show the bin width. Predictions from POWHEG (PH) DR and DS + PYTHIA 8 (P8), POWHEG DR + HERWIG 7 (H7), MADGRAPH5_aMC@NLO (aMC) DR, DR2, DS and DS with a dynamic factor + PYTHIA 8 are also shown. The grey band represents the statistical uncertainty and the yellow band the total uncertainty. In the lower panels, the ratio of the predictions to the data is shown.

The distribution of the Subleading jet $p_{T}$ for events in the 2j2b region. The number of observed events (points) and estimated signal and background events (filled histograms) from the maximum likelihood fit are shown. The vertical bars on the points represent the statistical uncertainty in the data, and the hatched band the total uncertainty in the estimated events after the fit. The lower panels display the ratio of the data to the sum of the estimated events (points) after the fit, with the bands giving the corresponding uncertainties.

Normalised fiducial differential tW production cross section as a function of the $p_{Z}(e^{\pm}, \mu^{\pm}, j)$. The horizontal bars on the points show the bin width. Predictions from POWHEG (PH) DR and DS + PYTHIA 8 (P8), POWHEG DR + HERWIG 7 (H7), MADGRAPH5_aMC@NLO (aMC) DR, DR2, DS and DS with a dynamic factor + PYTHIA 8 are also shown. The grey band represents the statistical uncertainty and the yellow band the total uncertainty. In the lower panels, the ratio of the predictions to the data is shown.

Observed and theoretical cross sections. In the observed, the first uncertainty is statistical, the second is the systematic, and the third the luminosity. In the theoretical, the first uncertainty is due to scale variations, the second due to the choice of PDF. The theoretical cross section has been computed at approximate third-order in QCD, with the addition of third-order corrections of soft gluon emission terms, assuming a top quark mass of 172.5 GeV and using the PDF4LHC21 PDF set.

Normalised fiducial differential tW production cross section as a function of the $Jet$ $p_{T}$. The horizontal bars on the points show the bin width. Predictions from POWHEG (PH) DR and DS + PYTHIA 8 (P8), POWHEG DR + HERWIG 7 (H7), MADGRAPH5_aMC@NLO (aMC) DR, DR2, DS and DS with a dynamic factor + PYTHIA 8 are also shown. The grey band represents the statistical uncertainty and the yellow band the total uncertainty. In the lower panels, the ratio of the predictions to the data is shown.

The number of estimated signal and background events after the fit in the 1j1b, 2j1b, and 2j2b regions compared to the observed number of events. The total uncertainties in the estimated events after the fit are given.

Normalised fiducial differential tW production cross section as a function of the $m(e^{\pm}, \mu^{\pm}, j)$. The horizontal bars on the points show the bin width. Predictions from POWHEG (PH) DR and DS + PYTHIA 8 (P8), POWHEG DR + HERWIG 7 (H7), MADGRAPH5_aMC@NLO (aMC) DR, DR2, DS and DS with a dynamic factor + PYTHIA 8 are also shown. The grey band represents the statistical uncertainty and the yellow band the total uncertainty. In the lower panels, the ratio of the predictions to the data is shown.

Normalised fiducial differential tW production cross section as a function of the $Leading$ $lepton$ $p_{T}$. The horizontal bars on the points show the bin width. Predictions from POWHEG (PH) DR and DS + PYTHIA 8 (P8), POWHEG DR + HERWIG 7 (H7), MADGRAPH5_aMC@NLO (aMC) DR, DR2, DS and DS with a dynamic factor + PYTHIA 8 are also shown. The grey band represents the statistical uncertainty and the yellow band the total uncertainty. In the lower panels, the ratio of the predictions to the data is shown.

Normalised fiducial differential tW production cross section as a function of the $\Delta\phi(e^{\pm}, \mu^{\pm})/\pi$. The horizontal bars on the points show the bin width. Predictions from POWHEG (PH) DR and DS + PYTHIA 8 (P8), POWHEG DR + HERWIG 7 (H7), MADGRAPH5_aMC@NLO (aMC) DR, DR2, DS and DS with a dynamic factor + PYTHIA 8 are also shown. The grey band represents the statistical uncertainty and the yellow band the total uncertainty. In the lower panels, the ratio of the predictions to the data is shown.

Normalised fiducial differential tW production cross section as a function of the $p_{Z}(e^{\pm}, \mu^{\pm}, j)$. The horizontal bars on the points show the bin width. Predictions from POWHEG (PH) DR and DS + PYTHIA 8 (P8), POWHEG DR + HERWIG 7 (H7), MADGRAPH5_aMC@NLO (aMC) DR, DR2, DS and DS with a dynamic factor + PYTHIA 8 are also shown. The grey band represents the statistical uncertainty and the yellow band the total uncertainty. In the lower panels, the ratio of the predictions to the data is shown.

Normalised fiducial differential tW production cross section as a function of the $m_{T}(e^{\pm}, \mu^{\pm}, j, p_{T}^{miss})$. The horizontal bars on the points show the bin width. Predictions from POWHEG (PH) DR and DS + PYTHIA 8 (P8), POWHEG DR + HERWIG 7 (H7), MADGRAPH5_aMC@NLO (aMC) DR, DR2, DS and DS with a dynamic factor + PYTHIA 8 are also shown. The grey band represents the statistical uncertainty and the yellow band the total uncertainty. In the lower panels, the ratio of the predictions to the data is shown.

Normalised fiducial differential tW production cross section as a function of the $Jet$ $p_{T}$. The horizontal bars on the points show the bin width. Predictions from POWHEG (PH) DR and DS + PYTHIA 8 (P8), POWHEG DR + HERWIG 7 (H7), MADGRAPH5_aMC@NLO (aMC) DR, DR2, DS and DS with a dynamic factor + PYTHIA 8 are also shown. The grey band represents the statistical uncertainty and the yellow band the total uncertainty. In the lower panels, the ratio of the predictions to the data is shown.

The p-values from the goodness-of-fit tests comparing the six differential cross section measurements with the predictions from POWHEG (PH) DR and DS + PYTHIA 8 (P8), POWHEG DR + HERWIG 7 (H7), MADGRAPH5 aMC@NLO (aMC) DR, DR2, DS, and DS with a dynamic factor + PYTHIA 8. The complete covariance matrix from the results and the statistical uncertainties in the predictions are taken into account.

Normalised fiducial differential tW production cross section as a function of the $m(e^{\pm}, \mu^{\pm}, j)$. The horizontal bars on the points show the bin width. Predictions from POWHEG (PH) DR and DS + PYTHIA 8 (P8), POWHEG DR + HERWIG 7 (H7), MADGRAPH5_aMC@NLO (aMC) DR, DR2, DS and DS with a dynamic factor + PYTHIA 8 are also shown. The grey band represents the statistical uncertainty and the yellow band the total uncertainty. In the lower panels, the ratio of the predictions to the data is shown.

Distribution of the RF output for events in the 1j1b region. The number of observed events (points) and estimated signal and background events (filled histograms) before the maximum likelihood fit are shown. The vertical bars on the points represent the statistical uncertainty in the data, and the hatched band the total uncertainty in the estimated events before the fit. The lower panels display the ratio of the data to the sum of the estimated events (points) before the fit, with the bands giving the corresponding uncertainties.

Normalised fiducial differential tW production cross section as a function of the $\Delta\phi(e^{\pm}, \mu^{\pm})/\pi$. The horizontal bars on the points show the bin width. Predictions from POWHEG (PH) DR and DS + PYTHIA 8 (P8), POWHEG DR + HERWIG 7 (H7), MADGRAPH5_aMC@NLO (aMC) DR, DR2, DS and DS with a dynamic factor + PYTHIA 8 are also shown. The grey band represents the statistical uncertainty and the yellow band the total uncertainty. In the lower panels, the ratio of the predictions to the data is shown.

Distribution of the RF output for events in the 2j1b region. The number of observed events (points) and estimated signal and background events (filled histograms) before the maximum likelihood fit are shown. The vertical bars on the points represent the statistical uncertainty in the data, and the hatched band the total uncertainty in the estimated events before the fit. The lower panels display the ratio of the data to the sum of the estimated events (points) before the fit, with the bands giving the corresponding uncertainties.

Normalised fiducial differential tW production cross section as a function of the $m_{T}(e^{\pm}, \mu^{\pm}, j, p_{T}^{miss})$. The horizontal bars on the points show the bin width. Predictions from POWHEG (PH) DR and DS + PYTHIA 8 (P8), POWHEG DR + HERWIG 7 (H7), MADGRAPH5_aMC@NLO (aMC) DR, DR2, DS and DS with a dynamic factor + PYTHIA 8 are also shown. The grey band represents the statistical uncertainty and the yellow band the total uncertainty. In the lower panels, the ratio of the predictions to the data is shown.

Distribution of the Subleading jet $p_{T}$ for events in the 2j2b region. The number of observed events (points) and estimated signal and background events (filled histograms) before the maximum likelihood fit are shown. The vertical bars on the points represent the statistical uncertainty in the data, and the hatched band the total uncertainty in the estimated events before the fit. The lower panels display the ratio of the data to the sum of the estimated events (points) before the fit, with the bands giving the corresponding uncertainties.

The p-values from the goodness-of-fit tests comparing the six differential cross section measurements with the predictions from POWHEG (PH) DR and DS + PYTHIA 8 (P8), POWHEG DR + HERWIG 7 (H7), MADGRAPH5 aMC@NLO (aMC) DR, DR2, DS, and DS with a dynamic factor + PYTHIA 8. The complete covariance matrix from the results and the statistical uncertainties in the predictions are taken into account.

Response matrix between detector and particle level for the $Leading$ $lepton$ $p_{T}$.

Response matrix between detector and particle level for the $Leading$ $lepton$ $p_{T}$.

Response matrix between detector and particle level for the $p_{Z}(e^{\pm}, \mu^{\pm}, j)$.

Response matrix between detector and particle level for the $p_{Z}(e^{\pm}, \mu^{\pm}, j)$.

Response matrix between detector and particle level for the $Jet$ $p_{T}$.

Response matrix between detector and particle level for the $Jet$ $p_{T}$.

Response matrix between detector and particle level for the $m(e^{\pm}, \mu^{\pm}, j)$.

Response matrix between detector and particle level for the $m(e^{\pm}, \mu^{\pm}, j)$.

Response matrix between detector and particle level for the $\Delta\phi(e^{\pm}, \mu^{\pm})/\pi$.

Response matrix between detector and particle level for the $\Delta\phi(e^{\pm}, \mu^{\pm})/\pi$.

Response matrix between detector and particle level for the $m_{T}(e^{\pm}, \mu^{\pm}, j, p_{T}^{miss})$.

Response matrix between detector and particle level for the $m_{T}(e^{\pm}, \mu^{\pm}, j, p_{T}^{miss})$.

Covariance matrix including all uncertainties of the normalised differential cross section in bins of $Leading$ $lepton$ $p_{T}$ in units of $1/GeV^{2}$.

Covariance matrix including all uncertainties of the normalised differential cross section in bins of $Leading$ $lepton$ $p_{T}$ in units of $1/GeV^{2}$.

Covariance matrix including all uncertainties of the normalised differential cross section in bins of $p_{Z}(e^{\pm}, \mu^{\pm}, j)$ in units of $1/GeV^{2}$.

Covariance matrix including all uncertainties of the normalised differential cross section in bins of $p_{Z}(e^{\pm}, \mu^{\pm}, j)$ in units of $1/GeV^{2}$.

Covariance matrix including all uncertainties of the normalised differential cross section in bins of $Jet$ $p_{T}$ in units of $1/GeV^{2}$.

Covariance matrix including all uncertainties of the normalised differential cross section in bins of $Jet$ $p_{T}$ in units of $1/GeV^{2}$.

Covariance matrix including all uncertainties of the normalised differential cross section in bins of $m(e^{\pm}, \mu^{\pm}, j)$ in units of $1/GeV^{2}$.

Covariance matrix including all uncertainties of the normalised differential cross section in bins of $m(e^{\pm}, \mu^{\pm}, j)$ in units of $1/GeV^{2}$.

Covariance matrix including all uncertainties of the normalised differential cross section in bins of $\Delta\phi(e^{\pm}, \mu^{\pm})/\pi$ in units of 1.

Covariance matrix including all uncertainties of the normalised differential cross section in bins of $\Delta\phi(e^{\pm}, \mu^{\pm})/\pi$ in units of 1.

Covariance matrix including all uncertainties of the normalised differential cross section in bins of $m_{T}(e^{\pm}, \mu^{\pm}, j, p_{T}^{miss})$ in units of $1/GeV^{2}$.

Covariance matrix including all uncertainties of the normalised differential cross section in bins of $m_{T}(e^{\pm}, \mu^{\pm}, j, p_{T}^{miss})$ in units of $1/GeV^{2}$.

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