$\Delta\phi$ projection of balance function in low $p_{T}$ in 0-10% centrality

$\Delta\phi$ projection of balance function in low $p_{T}$ in 30-40% centrality

1D $\Delta\eta$ projection of balance function in low $p_{T}$ in 70-80% centrality

$\Delta\eta$ projection of balance function in low $p_{T}$ in 0-40 $N_{trk}$

$\Delta\eta$ projection of balance function in low $p_{T}$ in 120-150 $N_{trk}$

$\Delta\eta$ projection of balance function in low $p_{T}$ in 270-300 $N_{trk}$

$\Delta\phi$ projection of balance function in low $p_{T}$ in 0-40 $N_{trk}$

$\Delta\phi$ projection of balance function in low $p_{T}$ in 120-150 $N_{trk}$

$\Delta\phi$ projection of balance function in low $p_{T}$ in 270-300 $N_{trk}$

The width of balance function in $<|\Delta\eta|>$ in low $p_{T}$ for PbPb collisions

Balance function width in low $p_{T}$ PbPb relative $<|\Delta\eta|>/<|\Delta\eta|>_{N_{ch}<65}$

The width of balance function in $<|\Delta\eta|>$ in low $p_{T}$ for pPb collisions

Balance function width in low $p_{T}$ in pPb relative $<|\Delta\eta|>/<|\Delta\eta|>_{N_{ch}<24}$

The width of the balance function in $<|\Delta\phi|>$ in low $p_{T}$ for PbPb collisions

Balance function width in low $p_{T}$ in PbPb for relative $<|\Delta\phi|>/<|\Delta\phi|>_{N_{ch}<65}$

The width of the balance function in $<|\Delta\phi|>$ in low $p_{T}$ for pPb collisions

Balance function width in low $p_{T}$ in pPb for relative $|\Delta\phi|/<|\Delta\phi|>_{N_{ch}<24}$

Balance function projection as a function of $\Delta\eta$ in intermediate $p_{T}$ in PbPb for 0-10% centrality

Balance function projection as a function of $\Delta\eta$ in intermediate $p_{T}$ in PbPb for 30-40% centrality

Balance function projection as a function of $\Delta\eta$ in intermediate $p_{T}$ in PbPb for 70-80% centrality

Balance function projection as a function of $\Delta\phi$ in intermediate $p_{T}$ in PbPb for 0-10% centrality

Balance function projection as a function of $\Delta\phi$ in intermediate $p_{T}$ in PbPb for 30-40% centrality

Balance function projection as a function of $\Delta\phi$ in intermediate $p_{T}$ in PbPb for 70-80% centrality

Balance function projection as a function of $\Delta\eta$ in intermediate $p_{T}$ in pPb for 0-40 $N_{trk}$

Balance function projection as a function of $\Delta\eta$ in intermediate $p_{T}$ in pPb for 120-150 $N_{trk}$

Balance function projection as a function of $\Delta\eta$ in intermediate $p_{T}$ in pPb for 270-300 $N_{trk}$

Balance function projection as a function of $\Delta\phi$ in intermediate $p_{T}$ in pPb for $N_{trk}$ 0-40

Balance function projection as a function of $\Delta\phi$ in intermediate $p_{T}$ in pPb for $N_{trk}$ 120-150

Balance function projection as a function of $\Delta\phi$ in intermediate $p_{T}$ in pPb for $N_{trk}$ 270-300

Balance function projection as a function of $\Delta\eta$ in high $p_{T}$ in PbPb for 0-10% centrality

Balance function projection as a function of $\Delta\eta$ in high $p_{T}$ in PbPb for 30-40% centrality

Balance function projection as a function of $\Delta\eta$ in high $p_{T}$ in PbPb for 70-80% centrality

Balance function projection as a function of $\Delta\phi$ in high $p_{T}$ in PbPb for 0-10% centrality

Balance function projection as a function of $\Delta\phi$ in high $p_{T}$ in PbPb for 30-40% centrality

Balance function projection as a function of $\Delta\phi$ in high $p_{T}$ in PbPb for 70-80% centrality

Balance function projection as a function of $\Delta\eta$ in high $p_{T}$ in pPb for $N_{trk}$ 0-40

Balance function projection as a function of $\Delta\eta$ in high $p_{T}$ in pPb for $N_{trk}$ 120-150

Balance function projection as a function of $\Delta\eta$ in high $p_{T}$ in pPb for $N_{trk}$ 270-300

Balance function projection as a function of $\Delta\phi$ in high $p_{T}$ in PbPb for $N_{trk}$ 0-40

Balance function projection as a function of $\Delta\phi$ in high $p_{T}$ in PbPb for $N_{trk}$ 120-150

Balance function projection as a function of $\Delta\phi$ in high $p_{T}$ in PbPb for $N_{trk}$ 270-300

Average width of the balance function in $\Delta \eta$ with $N_{ch}$ for low $p_{T}$

Average width of the balance function in $\Delta \eta$ with $N_{ch}$ for intermediate $p_{T}$

Average width of the balance function in $\Delta \eta$ with $N_{ch}$ for high $p_{T}$

Average width of the balance function in $\Delta \phi$ with $p_{T}$ for low $p_{T}$

Average width of the balance function in $\Delta \phi$ with $N_{ch}$ for intermediate $p_{T}$

Average width of the balance function in $\Delta \phi$ with $N_{ch}$ for high $p_{T}$

Average width of the balance function in $\Delta \eta$ with $N_{ch}$ for low $p_{T}$

Average width of the balance function in $\Delta \eta$ with $N_{ch}$ for intermediate $p_{T}$

Average width of the balance function in $\Delta \phi$ with $N_{ch}$ for high $p_{T}$

Average width of the balance function in $\Delta \phi$ with $N_{ch}$ for low $p_{T}$

Average width of the balance function in $\Delta \phi$ with $N_{ch}$ for intermediate $p_{T}$

Average width of the balance function in $\Delta \phi$ with $N_{ch}$ for high $p_{T}$

Expected exclusion limits at 95% CL from Fig 7(a) for $\tilde{g}$ decays into SM gauge bosons and $\tilde{\chi}^{0}_{1}$

One $\sigma$ band of expected exclusion limits at 95% CL from Fig 7(a) for $\tilde{g}$ decays into SM gauge bosons and $\tilde{\chi}^{0}_{1}$

Observed exclusion limits at 95% CL from Fig 7(c) for $\tilde{g}$ decays into sleptons and subsequently to SM leptons and $\tilde{\chi}^{0}_{1}$

Positive one $\sigma$ observed exclusion limits at 95% CL from Fig 7(c) for $\tilde{g}$ decays into sleptons and subsequently to SM leptons and $\tilde{\chi}^{0}_{1}$

Negative one $\sigma$ observed exclusion limits at 95% CL from Fig 7(c) for $\tilde{g}$ decays into sleptons and subsequently to SM leptons and $\tilde{\chi}^{0}_{1}$

Expected exclusion limits at 95% CL from Fig 7(c) for $\tilde{g}$ decays into sleptons and subsequently to SM leptons and $\tilde{\chi}^{0}_{1}$

One $\sigma$ band of expected exclusion limits at 95% CL from Fig 7(c) for $\tilde{g}$ decays into sleptons and subsequently to SM leptons and $\tilde{\chi}^{0}_{1}$

Observed exclusion limits at 95% CL from Fig 7(f) for $\tilde{g}$ decays into anti-top and $\tilde{t}$ and $\tilde{t}$ decays via a non-zero RPV coupling $\lambda''$

Positive one $\sigma$ observed exclusion limits at 95% CL from Fig 7(f) for $\tilde{g}$ decays into anti-top and $\tilde{t}$ and $\tilde{t}$ decays via a non-zero RPV coupling $\lambda''$

Negative one $\sigma$ observed exclusion limits at 95% CL from Fig 7(f) for $\tilde{g}$ decays into anti-top and $\tilde{t}$ and $\tilde{t}$ decays via a non-zero RPV coupling $\lambda''$

Expected exclusion limits at 95% CL from Fig 7(f) for $\tilde{g}$ decays into anti-top and $\tilde{t}$ and $\tilde{t}$ decays via a non-zero RPV coupling $\lambda''$

One $\sigma$ band of expected exclusion limits at 95% CL from Fig 7(f) for $\tilde{g}$ decays into anti-top and $\tilde{t}$ and $\tilde{t}$ decays via a non-zero RPV coupling $\lambda''$

Observed exclusion limits at 95% CL from Fig 7(e) for direct $\tilde{\chi_{1}^{0}}$ decay into SM leptons and quarks via a non-zero RPV coupling $\lambda'$

Positive one $\sigma$ observed exclusion limits at 95% CL from Fig 7(e) for direct $\tilde{\chi_{1}^{0}}$ decay into SM leptons and quarks via a non-zero RPV coupling $\lambda'$

Negative one $\sigma$ observed exclusion limits at 95% CL from Fig 7(e) for direct $\tilde{\chi_{1}^{0}}$ decay into SM leptons and quarks via a non-zero RPV coupling $\lambda'$

Expected exclusion limits at 95% CL from Fig 7(e) for direct $\tilde{\chi_{1}^{0}}$ decay into SM leptons and quarks via a non-zero RPV coupling $\lambda'$

One $\sigma$ band of expected exclusion limits at 95% CL from Fig 7(e) for direct $\tilde{\chi_{1}^{0}}$ decay into SM leptons and quarks via a non-zero RPV coupling $\lambda'$

Observed exclusion limits at 95% CL from Fig 7(b) for $\tilde{q}$ decays into SM gauge bosons and $\tilde{\chi}^{0}_{1}$

Positive one $\sigma$ observed exclusion limits at 95% CL from Fig 7(b) for $\tilde{q}$ decays into SM gauge bosons and $\tilde{\chi}^{0}_{1}$

Negative one $\sigma$ observed exclusion limits at 95% CL from Fig 7(b) for $\tilde{q}$ decays into SM gauge bosons and $\tilde{\chi}^{0}_{1}$

Expected exclusion limits at 95% CL from Fig 7(b) for $\tilde{q}$ decays into SM gauge bosons and $\tilde{\chi}^{0}_{1}$

One $\sigma$ band of expected exclusion limits at 95% CL from Fig 7(b) for $\tilde{q}$ decays into SM gauge bosons and $\tilde{\chi}^{0}_{1}$

Observed exclusion limits at 95% CL from Fig 7(d) for $\tilde{q}$ decays into sleptons and subsequently to SM leptons and $\tilde{\chi}^{0}_{1}$

Positive one $\sigma$ observed exclusion limits at 95% CL from Fig 7(d) for $\tilde{q}$ decays into sleptons and subsequently to SM leptons and $\tilde{\chi}^{0}_{1}$

Negative one $\sigma$ observed exclusion limits at 95% CL from Fig 7(d) for $\tilde{q}$ decays into sleptons and subsequently to SM leptons and $\tilde{\chi}^{0}_{1}$

Expected exclusion limits at 95% CL from Fig 7(d) for $\tilde{q}$ decays into sleptons and subsequently to SM leptons and $\tilde{\chi}^{0}_{1}$

One $\sigma$ band of expected exclusion limits at 95% CL from Fig 7(d) for $\tilde{q}$ decays into sleptons and subsequently to SM leptons and $\tilde{\chi}^{0}_{1}$

N-1 distribution for $m_{\mathrm{eff}}$of observed data and expected background in SRGGWZ-H.

N-1 distribution for $E_{\mathrm{T}}^{\mathrm{miss}}$of observed data and expected background in SRGGSlep-M.

N-1 distribution for $\sum{p_{\mathrm{T}}^\mathrm{jet}}$of observed data and expected background in SRUDD-ge2b.

N-1 distribution for $m_{\mathrm{eff}}$of observed data and expected background in SRLQD.

N-1 distribution for $m_{\mathrm{eff}}$of observed data and expected background in SRSSWZ-H.

N-1 distribution for $m_{\mathrm{eff}}$of observed data and expected background in SRSSSlep-H(loose).

Signal acceptance for SRGGWZ-H signal region from Fig 10(c) in a SUSY scenario where $\tilde{g}$ decays into SM gauge bosons and $\tilde{\chi}^{0}_{1}$

Signal efficiency for SRGGWZ-H signal region from Fig 15(c) in a SUSY scenario where $\tilde{g}$ decays into SM gauge bosons and $\tilde{\chi}^{0}_{1}$

Signal acceptance for SRGGWZ-M signal region from Fig 10(b) in a SUSY scenario where $\tilde{g}$ decays into SM gauge bosons and $\tilde{\chi}^{0}_{1}$

Signal efficiency for SRGGWZ-M signal region from Fig 15(b) in a SUSY scenario where $\tilde{g}$ decays into SM gauge bosons and $\tilde{\chi}^{0}_{1}$

Signal acceptance for SRGGWZ-L signal region from Fig 10(a) in a SUSY scenario where $\tilde{g}$ decays into SM gauge bosons and $\tilde{\chi}^{0}_{1}$

Signal efficiency for SRGGWZ-L signal region from Fig 15(a) in a SUSY scenario where $\tilde{g}$ decays into SM gauge bosons and $\tilde{\chi}^{0}_{1}$

Signal acceptance for SRGGSlep-L signal region from Fig 12(a) in a SUSY scenario where $\tilde{g}$ decays into sleptons and subsequently to SM leptons and $\tilde{\chi}^{0}_{1}$

Signal efficiency for SRGGSlep-L signal region from Fig 17(a) in a SUSY scenario where $\tilde{g}$ decays into sleptons and subsequently to SM leptons and $\tilde{\chi}^{0}_{1}$

Signal acceptance for SRGGSlep-M signal region from Fig 12(b) in a SUSY scenario where $\tilde{g}$ decays into sleptons and subsequently to SM leptons and $\tilde{\chi}^{0}_{1}$

Signal efficiency for SRGGSlep-M signal region from Fig 17(b) in a SUSY scenario where $\tilde{g}$ decays into sleptons and subsequently to SM leptons and $\tilde{\chi}^{0}_{1}$

Signal acceptance for SRGGSlep-H signal region from Fig 12(c) in a SUSY scenario where $\tilde{g}$ decays into sleptons and subsequently to SM leptons and $\tilde{\chi}^{0}_{1}$

Signal efficiency for SRGGSlep-H signal region from Fig 17(c) in a SUSY scenario where $\tilde{g}$ decays into sleptons and subsequently to SM leptons and $\tilde{\chi}^{0}_{1}$

Signal acceptance for SRUDD-1b signal region from Fig 14(b) in a SUSY scenario where $\tilde{g}$ decays into anti-top and $\tilde{t}$ and $\tilde{t}$ decays via a non-zero RPV coupling $\lambda''$

Signal efficiency for SRUDD-1b signal region from Fig 19(b) in a SUSY scenario where $\tilde{g}$ decays into anti-top and $\tilde{t}$ and $\tilde{t}$ decays via a non-zero RPV coupling $\lambda''$

Signal acceptance for SRUDD-2b signal region from Fig 14(c) in a SUSY scenario where $\tilde{g}$ decays into anti-top and $\tilde{t}$ and $\tilde{t}$ decays via a non-zero RPV coupling $\lambda''$

Signal efficiency for SRUDD-2b signal region from Fig 19(c) in a SUSY scenario where $\tilde{g}$ decays into anti-top and $\tilde{t}$ and $\tilde{t}$ decays via a non-zero RPV coupling $\lambda''$

Signal acceptance for SRUDD-ge2b signal region from Fig 14(d) in a SUSY scenario where $\tilde{g}$ decays into anti-top and $\tilde{t}$ and $\tilde{t}$ decays via a non-zero RPV coupling $\lambda''$

Signal efficiency for SRUDD-ge2b signal region from Fig 19(d) in a SUSY scenario where $\tilde{g}$ decays into anti-top and $\tilde{t}$ and $\tilde{t}$ decays via a non-zero RPV coupling $\lambda''$

Signal acceptance for SRUDD-ge3b signal region from Fig 14(e) in a SUSY scenario where $\tilde{g}$ decays into anti-top and $\tilde{t}$ and $\tilde{t}$ decays via a non-zero RPV coupling $\lambda''$

Signal efficiency for SRUDD-ge3b signal region from Fig 19(e) in a SUSY scenario where $\tilde{g}$ decays into anti-top and $\tilde{t}$ and $\tilde{t}$ decays via a non-zero RPV coupling $\lambda''$

Signal acceptance for SRLQD signal region from Fig 14(a) in a SUSY scenario where direct $\tilde{\chi_{1}^{0}}$ decay into SM leptons and quarks via a non-zero RPV coupling $\lambda'$

Signal efficiency for SRLQD signal region from Fig 19(a) in a SUSY scenario where direct $\tilde{\chi_{1}^{0}}$ decay into SM leptons and quarks via a non-zero RPV coupling $\lambda'$

Signal acceptance for SRSSWZ-L signal region from Fig 11(a) in a SUSY scenario where $\tilde{q}$ decays into SM gauge bosons and $\tilde{\chi}^{0}_{1}$

Signal efficiency for SRSSWZ-L signal region from Fig 16(a) in a SUSY scenario where $\tilde{q}$ decays into SM gauge bosons and $\tilde{\chi}^{0}_{1}$

Signal acceptance for SRSSWZ-ML signal region from Fig 11(b) in a SUSY scenario where $\tilde{q}$ decays into SM gauge bosons and $\tilde{\chi}^{0}_{1}$

Signal efficiency for SRSSWZ-ML signal region from Fig 16(b) in a SUSY scenario where $\tilde{q}$ decays into SM gauge bosons and $\tilde{\chi}^{0}_{1}$

Signal acceptance for SRSSWZ-MH signal region from Fig 11(c) in a SUSY scenario where $\tilde{q}$ decays into SM gauge bosons and $\tilde{\chi}^{0}_{1}$

Signal efficiency for SRSSWZ-MH signal region from Fig 16(c) in a SUSY scenario where $\tilde{q}$ decays into SM gauge bosons and $\tilde{\chi}^{0}_{1}$

Signal acceptance for SRSSWZ-H signal region from Fig 11(d) in a SUSY scenario where $\tilde{q}$ decays into SM gauge bosons and $\tilde{\chi}^{0}_{1}$

Signal efficiency for SRSSWZ-H signal region from Fig 16(d) in a SUSY scenario where $\tilde{q}$ decays into SM gauge bosons and $\tilde{\chi}^{0}_{1}$

Signal acceptance for SRSSSlep-H signal region from Fig 13(d) in a SUSY scenario where $\tilde{q}$ decays into sleptons and subsequently to SM leptons and $\tilde{\chi}^{0}_{1}$

Signal efficiency for SRSSSlep-H signal region from Fig 18(d) in a SUSY scenario where $\tilde{q}$ decays into sleptons and subsequently to SM leptons and $\tilde{\chi}^{0}_{1}$

Signal acceptance for SRSSSlep-MH signal region from Fig 13(c) in a SUSY scenario where $\tilde{q}$ decays into sleptons and subsequently to SM leptons and $\tilde{\chi}^{0}_{1}$

Signal efficiency for SRSSSlep-MH signal region from Fig 18(c) in a SUSY scenario where $\tilde{q}$ decays into sleptons and subsequently to SM leptons and $\tilde{\chi}^{0}_{1}$

Signal acceptance for SRSSSlep-L signal region from Fig 13(a) in a SUSY scenario where $\tilde{q}$ decays into sleptons and subsequently to SM leptons and $\tilde{\chi}^{0}_{1}$

Signal efficiency for SRSSSlep-L signal region from Fig 18(a) in a SUSY scenario where $\tilde{q}$ decays into sleptons and subsequently to SM leptons and $\tilde{\chi}^{0}_{1}$

Signal acceptance for SRSSSlep-ML signal region from Fig 13(b) in a SUSY scenario where $\tilde{q}$ decays into sleptons and subsequently to SM leptons and $\tilde{\chi}^{0}_{1}$

Signal efficiency for SRSSSlep-ML signal region from Fig 18(b) in a SUSY scenario where $\tilde{q}$ decays into sleptons and subsequently to SM leptons and $\tilde{\chi}^{0}_{1}$

Signal acceptance for SRSSSlep-H(loose) signal region from Fig 13(e) in a SUSY scenario where $\tilde{q}$ decays into sleptons and subsequently to SM leptons and $\tilde{\chi}^{0}_{1}$

Signal efficiency for SRSSSlep-H(loose) signal region from Fig 18(e) in a SUSY scenario where $\tilde{q}$ decays into sleptons and subsequently to SM leptons and $\tilde{\chi}^{0}_{1}$

Number of signal events expected for 139 fb$^{-1}$ at different stages of the event selection for the signal region SRGGWZ-H in a susy scenario where $\tilde{g}$ decays into SM gauge bosons and $\tilde{\chi}^{0}_{1}$. The masses of the superpartners involved in the process are set to $m(\tilde{g})$ = 1400 GeV, $m(\tilde{\chi_{1}^{0}})$ = 1000 GeV. Only statistical uncertainties are shown.

Number of signal events expected for 139 fb$^{-1}$ at different stages of the event selection for the signal region SRGGWZ-M in a susy scenario where $\tilde{g}$ decays into SM gauge bosons and $\tilde{\chi}^{0}_{1}$. The masses of the superpartners involved in the process are set to $m(\tilde{g})$ = 1400 GeV, $m(\tilde{\chi_{1}^{0}})$ = 1000 GeV. Only statistical uncertainties are shown.

Number of signal events expected for 139 fb$^{-1}$ at different stages of the event selection for the signal region SRGGWZ-L in a susy scenario where $\tilde{g}$ decays into SM gauge bosons and $\tilde{\chi}^{0}_{1}$. The masses of the superpartners involved in the process are set to $m(\tilde{g})$ = 1400 GeV, $m(\tilde{\chi_{1}^{0}})$ = 1000 GeV. Only statistical uncertainties are shown.

Number of signal events expected for 139 fb$^{-1}$ at different stages of the event selection for the signal region SRGGSlep-L in a susy scenario where $\tilde{g}$ decays into sleptons and subsequently to SM leptons and $\tilde{\chi}^{0}_{1}$. The masses of the superpartners involved in the process are set to $m(\tilde{g})$ = 2000 GeV, $m(\tilde{\chi_{1}^{0}})$ = 500 GeV. Only statistical uncertainties are shown.

Number of signal events expected for 139 fb$^{-1}$ at different stages of the event selection for the signal region SRGGSlep-M in a susy scenario where $\tilde{g}$ decays into sleptons and subsequently to SM leptons and $\tilde{\chi}^{0}_{1}$. The masses of the superpartners involved in the process are set to $m(\tilde{g})$ = 2000 GeV, $m(\tilde{\chi_{1}^{0}})$ = 500 GeV. Only statistical uncertainties are shown.

Number of signal events expected for 139 fb$^{-1}$ at different stages of the event selection for the signal region SRGGSlep-H in a susy scenario where $\tilde{g}$ decays into sleptons and subsequently to SM leptons and $\tilde{\chi}^{0}_{1}$. The masses of the superpartners involved in the process are set to $m(\tilde{g})$ = 2000 GeV, $m(\tilde{\chi_{1}^{0}})$ = 500 GeV. Only statistical uncertainties are shown.

Number of signal events expected for 139 fb$^{-1}$ at different stages of the event selection for the signal region SRUDD-1b in a susy scenario where $\tilde{g}$ decays into anti-top and $\tilde{t}$ and $\tilde{t}$ decays via a non-zero RPV coupling $\lambda''$. The masses of the superpartners involved in the process are set to $m(\tilde{g})$ = 1600 GeV, $m(\tilde{t})$ = 600 GeV. Only statistical uncertainties are shown.

Number of signal events expected for 139 fb$^{-1}$ at different stages of the event selection for the signal region SRUDD-2b in a susy scenario where $\tilde{g}$ decays into anti-top and $\tilde{t}$ and $\tilde{t}$ decays via a non-zero RPV coupling $\lambda''$. The masses of the superpartners involved in the process are set to $m(\tilde{g})$ = 1600 GeV, $m(\tilde{t})$ = 600 GeV. Only statistical uncertainties are shown.

Number of signal events expected for 139 fb$^{-1}$ at different stages of the event selection for the signal region SRUDD-ge2b in a susy scenario where $\tilde{g}$ decays into anti-top and $\tilde{t}$ and $\tilde{t}$ decays via a non-zero RPV coupling $\lambda''$. The masses of the superpartners involved in the process are set to $m(\tilde{g})$ = 1600 GeV, $m(\tilde{t})$ = 600 GeV. Only statistical uncertainties are shown.

Number of signal events expected for 139 fb$^{-1}$ at different stages of the event selection for the signal region SRUDD-ge3b in a susy scenario where $\tilde{g}$ decays into anti-top and $\tilde{t}$ and $\tilde{t}$ decays via a non-zero RPV coupling $\lambda''$. The masses of the superpartners involved in the process are set to $m(\tilde{g})$ = 1600 GeV, $m(\tilde{t})$ = 600 GeV. Only statistical uncertainties are shown.

Number of signal events expected for 139 fb$^{-1}$ at different stages of the event selection for the signal region SRLQD in a susy scenario where direct $\tilde{\chi_{1}^{0}}$ decay into SM leptons and quarks via a non-zero RPV coupling $\lambda'$. The masses of the superpartners involved in the process are set to $m(\tilde{g})$ = 2200 GeV, $m(\tilde{\chi_{1}^{0}})$ = 1870 GeV. Only statistical uncertainties are shown.

Number of signal events expected for 139 fb$^{-1}$ at different stages of the event selection for the signal region SRSSWZ-L in a susy scenario where $\tilde{q}$ decays into SM gauge bosons and $\tilde{\chi}^{0}_{1}$. The masses of the superpartners involved in the process are set to $m(\tilde{q})$ = 800 GeV, $m(\tilde{\chi_{1}^{0}})$ = 600 GeV. Only statistical uncertainties are shown.

Number of signal events expected for 139 fb$^{-1}$ at different stages of the event selection for the signal region SRSSWZ-ML in a susy scenario where $\tilde{q}$ decays into SM gauge bosons and $\tilde{\chi}^{0}_{1}$. The masses of the superpartners involved in the process are set to $m(\tilde{q})$ = 800 GeV, $m(\tilde{\chi_{1}^{0}})$ = 600 GeV. Only statistical uncertainties are shown.

Number of signal events expected for 139 fb$^{-1}$ at different stages of the event selection for the signal region SRSSWZ-MH in a susy scenario where $\tilde{q}$ decays into SM gauge bosons and $\tilde{\chi}^{0}_{1}$. The masses of the superpartners involved in the process are set to $m(\tilde{q})$ = 800 GeV, $m(\tilde{\chi_{1}^{0}})$ = 600 GeV. Only statistical uncertainties are shown.

Number of signal events expected for 139 fb$^{-1}$ at different stages of the event selection for the signal region SRSSWZ-H in a susy scenario where $\tilde{q}$ decays into SM gauge bosons and $\tilde{\chi}^{0}_{1}$. The masses of the superpartners involved in the process are set to $m(\tilde{q})$ = 800 GeV, $m(\tilde{\chi_{1}^{0}})$ = 600 GeV. Only statistical uncertainties are shown.

Number of signal events expected for 139 fb$^{-1}$ at different stages of the event selection for the signal region SRSSSlep-H in a susy scenario where $\tilde{q}$ decays into sleptons and subsequently to SM leptons and $\tilde{\chi}^{0}_{1}$. The masses of the superpartners involved in the process are set to $m(\tilde{q})$ = 1000 GeV, $m(\tilde{\chi_{1}^{0}})$ = 800 GeV. Only statistical uncertainties are shown.

Number of signal events expected for 139 fb$^{-1}$ at different stages of the event selection for the signal region SRSSSlep-MH in a susy scenario where $\tilde{q}$ decays into sleptons and subsequently to SM leptons and $\tilde{\chi}^{0}_{1}$. The masses of the superpartners involved in the process are set to $m(\tilde{q})$ = 1000 GeV, $m(\tilde{\chi_{1}^{0}})$ = 800 GeV. Only statistical uncertainties are shown.

Number of signal events expected for 139 fb$^{-1}$ at different stages of the event selection for the signal region SRSSSlep-L in a susy scenario where $\tilde{q}$ decays into sleptons and subsequently to SM leptons and $\tilde{\chi}^{0}_{1}$. The masses of the superpartners involved in the process are set to $m(\tilde{q})$ = 1000 GeV, $m(\tilde{\chi_{1}^{0}})$ = 800 GeV. Only statistical uncertainties are shown.

Number of signal events expected for 139 fb$^{-1}$ at different stages of the event selection for the signal region SRSSSlep-ML in a susy scenario where $\tilde{q}$ decays into sleptons and subsequently to SM leptons and $\tilde{\chi}^{0}_{1}$. The masses of the superpartners involved in the process are set to $m(\tilde{q})$ = 1000 GeV, $m(\tilde{\chi_{1}^{0}})$ = 800 GeV. Only statistical uncertainties are shown.

Number of signal events expected for 139 fb$^{-1}$ at different stages of the event selection for the signal region SRSSSlep-H(loose) in a susy scenario where $\tilde{q}$ decays into sleptons and subsequently to SM leptons and $\tilde{\chi}^{0}_{1}$. The masses of the superpartners involved in the process are set to $m(\tilde{q})$ = 1000 GeV, $m(\tilde{\chi_{1}^{0}})$ = 800 GeV. Only statistical uncertainties are shown.

Cross-section upper limits at 95% CL from Fig1(a) for $\tilde{g}$ decays into SM gauge bosons and $\tilde{\chi}^{0}_{1}$

Cross-section upper limits at 95% CL from Fig1(c) for $\tilde{g}$ decays into sleptons and subsequently to SM leptons and $\tilde{\chi}^{0}_{1}$

Cross-section upper limits at 95% CL from Fig1(f) for $\tilde{g}$ decays into anti-top and $\tilde{t}$ and $\tilde{t}$ decays via a non-zero RPV coupling $\lambda''$

Cross-section upper limits at 95% CL from Fig1(e) for direct $\tilde{\chi_{1}^{0}}$ decay into SM leptons and quarks via a non-zero RPV coupling $\lambda'$

Cross-section upper limits at 95% CL from Fig1(b) for $\tilde{q}$ decays into SM gauge bosons and $\tilde{\chi}^{0}_{1}$

Cross-section upper limits at 95% CL from Fig1(d) for $\tilde{q}$ decays into sleptons and subsequently to SM leptons and $\tilde{\chi}^{0}_{1}$

Unfolded $M(l^{-}\gamma)$ distribution for $Z \to \mu\mu\gamma$ process with dressed leptons and bkg subtraction. $M_{ll}>20$ GeV. Nexp.un f. = 64822.9 $\pm$ 254.6, NPowHeg truth =634276.

Unfolded dR distribution for $Z \to ee\gamma$ process with dressed leptons and bkg subtraction. $M_{ll}>20$ GeV. Nexp.un f. = 64273.2 $\pm$ 253.5, NPowHeg truth =338671.

Unfolded dR distribution for $Z \to \mu\mu\gamma$ process with dressed leptons and bkg subtraction. $M_{ll}>20$ GeV. Nexp.un f. = 65362.4 $\pm$ 255.7 , NPowHeg truth =634214.

Unfolded $P_{T}^{\gamma}$ distribution for $Z \to ee\gamma$ process with dressed leptons and bkg subtraction. $M_{ll}>20$ GeV. Nexp.un f. = 63486.8 $\pm$ 252.0, NPowHeg truth =333348.

Unfolded $P_{T}^{\gamma}$ distribution for $Z \to \mu\mu\gamma$ process with dressed leptons and bkg subtraction. $M_{ll}>20$ GeV. Nexp.un f. = 64432.6 $\pm$ 253.8, NPowHeg truth =624059.

Unfolded $M(l^{+}\gamma)$ distribution for $Z \to ee\gamma$ process with dressed leptons and bkg subtraction. $M_{ll}>45$ GeV. Nexp.un f. = 57095.9 $\pm$ 254 , NPowHeg truth =301666, N Sherpa truth =76327.

Unfolded $M(l^{-}\gamma)$ distribution for $Z \to ee\gamma$ process with dressed leptons and bkg subtraction. $M_{ll}>45$ GeV. Nexp.un f. = 57130 $\pm$ 254 , NPowHeg truth =301666, N Sherpa truth =763273.

Unfolded $M(l^{+}\gamma)$ distribution for $Z \to \mu\mu\gamma$ process with dressed leptons and bkg subtraction. $M_{ll}>45$ GeV. Nexp.un f. = 58310.9 $\pm$ 245, NPowHeg truth =564672, N Sherpa truth =814016.

Unfolded $M(l^{-}\gamma)$ distribution for $Z \to \mu\mu\gamma$ process with dressed leptons and bkg subtraction. $M_{ll}>45$ GeV. Nexp.un f. = 58320.9 $\pm$ 245, NPowHeg truth =564672, N Sherpa truth =814016.

Unfolded dR distribution for $Z \to ee\gamma$ process with dressed leptons and bkg subtraction. $M_{ll}>45$ GeV. Nexp.un f. = 57705.4 $\pm$ 254 , NPowHeg truth =301655, N Sherpa truth =763259.

Unfolded dR distribution for $Z \to \mu\mu\gamma$ process with dressed leptons and bkg subtraction. $M_{ll}>45$ GeV. Nexp.un f. = 58869.2 $\pm$ 256, NPowHeg truth =564660, NS herpa truth =813995.

Unfolded $P_{T}^{\gamma}$ distribution for $Z \to ee\gamma$ process with dressed leptons and bkg subtraction. $M_{ll}>45$ GeV. Nexp.un f. = 57358.6 $\pm$ 254 , NPowHeg truth =298662, N Sherpa truth =756609.

Unfolded $P_{T}^{\gamma}$ distribution for $Z \to \mu\mu\gamma$ process with dressed leptons and bkg subtraction. $M_{ll}>45$ GeV. Nexp.un f. = 58226.4 $\pm$ 254, NPowHeg truth =558949, N Sherpa truth =806988.

Unfolded $M(l^{+}\gamma)$ distribution for $Z \to ee\gamma$ process with bare leptons and bkg subtraction. $M_{ll}>20$ GeV. Nexp.un f. = 63717.4 $\pm$ 252.4, NPowHeg truth =338714.

Unfolded $M(l^{-}\gamma)$ distribution for $Z \to ee\gamma$ process with bare leptons and bkg subtraction. $M_{ll}>20$ GeV. Nexp.un f. = 63855.8 $\pm$ 252.7 , NPowHeg truth =338708.

Unfolded $M(l^{+}\gamma)$ distribution for $Z \to \mu\mu\gamma$ process with bare leptons and bkg subtraction. $M_{ll}>20$ GeV. Nexp.un f. = 64809.8 $\pm$ 254.6, NPowHeg truth =634285.

Unfolded $M(l^{-}\gamma)$ distribution for $Z \to \mu\mu\gamma$ process with bare leptons and bkg subtraction. $M_{ll}>20$ GeV. Nexp.un f. = 64822.9 $\pm$ 254.6, NPowHeg truth =634276.

Unfolded dR distribution for $Z \to ee\gamma$ process with bare leptons and bkg subtraction. $M_{ll}>20$ GeV. Nexp.un f. = 64273.2 $\pm$ 253.5, NPowHeg truth =338671.

Unfolded dR distribution for $Z \to \mu\mu\gamma$ process with bare leptons and bkg subtraction. $M_{ll}>20$ GeV. Nexp.un f. = 65362.4 $\pm$ 255.7 , NPowHeg truth =634214.

Unfolded $P_{T}^{\gamma}$ distribution for $Z \to ee\gamma$ process with bare leptons and bkg subtraction. $M_{ll}>20$ GeV. Nexp.un f. = 63486.8 $\pm$ 252.0, NPowHeg truth =333348.

Unfolded $P_{T}^{\gamma}$ distribution for $Z \to \mu\mu\gamma$ process with bare leptons and bkg subtraction. $M_{ll}>20$ GeV. Nexp.un f. = 64432.6 $\pm$ 253.8, NPowHeg truth =624059.

Unfolded $M(l^{+}\gamma)$ distribution for $Z \to ee\gamma$ process with bare leptons and bkg subtraction. $M_{ll}>45$ GeV. Nexp.un f. = 57095.9 $\pm$ 254 , NPowHeg truth =301666, N Sherpa truth =76327.

Unfolded $M(l^{-}\gamma)$ distribution for $Z \to ee\gamma$ process with bare leptons and bkg subtraction. $M_{ll}>45$ GeV. Nexp.un f. = 57130 $\pm$ 254 , NPowHeg truth =301666, N Sherpa truth =763273.

Unfolded $M(l^{+}\gamma)$ distribution for $Z \to \mu\mu\gamma$ process with bare leptons and bkg subtraction. $M_{ll}>45$ GeV. Nexp.un f. = 58310.9 $\pm$ 245, NPowHeg truth =564672, N Sherpa truth =814016.

Unfolded $M(l^{-}\gamma)$ distribution for $Z \to \mu\mu\gamma$ process with bare leptons and bkg subtraction. $M_{ll}>45$ GeV. Nexp.un f. = 58320.9 $\pm$ 245, NPowHeg truth =564672, N Sherpa truth =814016.

Unfolded dR distribution for $Z \to ee\gamma$ process with bare leptons and bkg subtraction. $M_{ll}>45$ GeV. Nexp.un f. = 57705.4 $\pm$ 254 , NPowHeg truth =301655, N Sherpa truth =763259.

Unfolded dR distribution for $Z \to \mu\mu\gamma$ process with bare leptons and bkg subtraction. $M_{ll}>45$ GeV. Nexp.un f. = 58869.2 $\pm$ 256, NPowHeg truth =564660, NS herpa truth =813995.

Unfolded $P_{T}^{\gamma}$ distribution for $Z \to ee\gamma$ process with bare leptons and bkg subtraction. $M_{ll}>45$ GeV. Nexp.un f. = 57358.6 $\pm$ 254 , NPowHeg truth =298662, N Sherpa truth =756609.

Unfolded $P_{T}^{\gamma}$ distribution for $Z \to \mu\mu\gamma$ process with bare leptons and bkg subtraction. $M_{ll}>45$ GeV. Nexp.un f. = 58226.4 $\pm$ 254, NPowHeg truth =558949, N Sherpa truth =806988.

Truth $P_{T}^{Z}$ distribution after reweghting to data.

Combined Covariance Matrix for Unfolded $M(l^{+}\gamma)$ distribution for $Z \to ee\gamma$ process with dressed leptons and bkg subtraction. $M_{ll}>20$ GeV. Nexp.un f. = 63717.4 $\pm$ 252.4, NPowHeg truth =338714.

Combined Covariance Matrix for Unfolded $M(l^{-}\gamma)$ distribution for $Z \to ee\gamma$ process with dressed leptons and bkg subtraction. $M_{ll}>20$ GeV. Nexp.un f. = 63855.8 $\pm$ 252.7 , NPowHeg truth =338708.

Combined Covariance Matrix for Unfolded $M(l^{+}\gamma)$ distribution for $Z \to \mu\mu\gamma$ process with dressed leptons and bkg subtraction. $M_{ll}>20$ GeV. Nexp.un f. = 64809.8 $\pm$ 254.6, NPowHeg truth =634285.

Combined Covariance Matrix for Unfolded $M(l^{-}\gamma)$ distribution for $Z \to \mu\mu\gamma$ process with dressed leptons and bkg subtraction. $M_{ll}>20$ GeV. Nexp.un f. = 64822.9 $\pm$ 254.6, NPowHeg truth =634276.

Combined Covariance Matrix for Unfolded dR distribution for $Z \to ee\gamma$ process with dressed leptons and bkg subtraction. $M_{ll}>20$ GeV. Nexp.un f. = 64273.2 $\pm$ 253.5, NPowHeg truth =338671.

Combined Covariance Matrix for Unfolded dR distribution for $Z \to \mu\mu\gamma$ process with dressed leptons and bkg subtraction. $M_{ll}>20$ GeV. Nexp.un f. = 65362.4 $\pm$ 255.7 , NPowHeg truth =634214.

Combined Covariance Matrix for Unfolded $P_{T}^{\gamma}$ distribution for $Z \to ee\gamma$ process with dressed leptons and bkg subtraction. $M_{ll}>20$ GeV. Nexp.un f. = 63486.8 $\pm$ 252.0, NPowHeg truth =333348.

Combined Covariance Matrix for Unfolded $P_{T}^{\gamma}$ distribution for $Z \to \mu\mu\gamma$ process with dressed leptons and bkg subtraction. $M_{ll}>20$ GeV. Nexp.un f. = 64432.6 $\pm$ 253.8, NPowHeg truth =624059.

Combined Covariance Matrix for Unfolded $M(l^{+}\gamma)$ distribution for $Z \to ee\gamma$ process with dressed leptons and bkg subtraction. $M_{ll}>45$ GeV. Nexp.un f. = 57095.9 $\pm$ 254 , NPowHeg truth =301666, N Sherpa truth =76327.

Combined Covariance Matrix for Unfolded $M(l^{-}\gamma)$ distribution for $Z \to ee\gamma$ process with dressed leptons and bkg subtraction. $M_{ll}>45$ GeV. Nexp.un f. = 57130 $\pm$ 254 , NPowHeg truth =301666, N Sherpa truth =763273.

Combined Covariance Matrix for Unfolded $M(l^{+}\gamma)$ distribution for $Z \to \mu\mu\gamma$ process with dressed leptons and bkg subtraction. $M_{ll}>45$ GeV. Nexp.un f. = 58310.9 $\pm$ 245, NPowHeg truth =564672, N Sherpa truth =814016.

Combined Covariance Matrix for Unfolded $M(l^{-}\gamma)$ distribution for $Z \to \mu\mu\gamma$ process with dressed leptons and bkg subtraction. $M_{ll}>45$ GeV. Nexp.un f. = 58320.9 $\pm$ 245, NPowHeg truth =564672, N Sherpa truth =814016.

Combined Covariance Matrix for Unfolded dR distribution for $Z \to ee\gamma$ process with dressed leptons and bkg subtraction. $M_{ll}>45$ GeV. Nexp.un f. = 57705.4 $\pm$ 254 , NPowHeg truth =301655, N Sherpa truth =763259.

Combined Covariance Matrix for Unfolded dR distribution for $Z \to \mu\mu\gamma$ process with dressed leptons and bkg subtraction. $M_{ll}>45$ GeV. Nexp.un f. = 58869.2 $\pm$ 256, NPowHeg truth =564660, NS herpa truth =813995.

Combined Covariance Matrix for Unfolded $P_{T}^{\gamma}$ distribution for $Z \to ee\gamma$ process with dressed leptons and bkg subtraction. $M_{ll}>45$ GeV. Nexp.un f. = 57358.6 $\pm$ 254 , NPowHeg truth =298662, N Sherpa truth =756609.

Combined Covariance Matrix for Unfolded $P_{T}^{\gamma}$ distribution for $Z \to \mu\mu\gamma$ process with dressed leptons and bkg subtraction. $M_{ll}>45$ GeV. Nexp.un f. = 58226.4 $\pm$ 254, NPowHeg truth =558949, N Sherpa truth =806988.

Combined Covariance Matrix for Unfolded $M(l^{+}\gamma)$ distribution for $Z \to ee\gamma$ process with bare leptons and bkg subtraction. $M_{ll}>20$ GeV. Nexp.un f. = 63717.4 $\pm$ 252.4, NPowHeg truth =338714.

Combined Covariance Matrix for Unfolded $M(l^{-}\gamma)$ distribution for $Z \to ee\gamma$ process with bare leptons and bkg subtraction. $M_{ll}>20$ GeV. Nexp.un f. = 63855.8 $\pm$ 252.7 , NPowHeg truth =338708.

Combined Covariance Matrix for Unfolded $M(l^{+}\gamma)$ distribution for $Z \to \mu\mu\gamma$ process with bare leptons and bkg subtraction. $M_{ll}>20$ GeV. Nexp.un f. = 64809.8 $\pm$ 254.6, NPowHeg truth =634285.

Combined Covariance Matrix for Unfolded $M(l^{-}\gamma)$ distribution for $Z \to \mu\mu\gamma$ process with bare leptons and bkg subtraction. $M_{ll}>20$ GeV. Nexp.un f. = 64822.9 $\pm$ 254.6, NPowHeg truth =634276.

Combined Covariance Matrix for Unfolded dR distribution for $Z \to ee\gamma$ process with bare leptons and bkg subtraction. $M_{ll}>20$ GeV. Nexp.un f. = 64273.2 $\pm$ 253.5, NPowHeg truth =338671.

Combined Covariance Matrix for Unfolded dR distribution for $Z \to \mu\mu\gamma$ process with bare leptons and bkg subtraction. $M_{ll}>20$ GeV. Nexp.un f. = 65362.4 $\pm$ 255.7 , NPowHeg truth =634214.

Combined Covariance Matrix for Unfolded $P_{T}^{\gamma}$ distribution for $Z \to ee\gamma$ process with bare leptons and bkg subtraction. $M_{ll}>20$ GeV. Nexp.un f. = 63486.8 $\pm$ 252.0, NPowHeg truth =333348.

Combined Covariance Matrix for Unfolded $P_{T}^{\gamma}$ distribution for $Z \to \mu\mu\gamma$ process with bare leptons and bkg subtraction. $M_{ll}>20$ GeV. Nexp.un f. = 64432.6 $\pm$ 253.8, NPowHeg truth =624059.

Combined Covariance Matrix for Unfolded $M(l^{+}\gamma)$ distribution for $Z \to ee\gamma$ process with bare leptons and bkg subtraction. $M_{ll}>45$ GeV. Nexp.un f. = 57095.9 $\pm$ 254 , NPowHeg truth =301666, N Sherpa truth =76327.

Combined Covariance Matrix for Unfolded $M(l^{-}\gamma)$ distribution for $Z \to ee\gamma$ process with bare leptons and bkg subtraction. $M_{ll}>45$ GeV. Nexp.un f. = 57130 $\pm$ 254 , NPowHeg truth =301666, N Sherpa truth =763273.

Combined Covariance Matrix for Unfolded $M(l^{+}\gamma)$ distribution for $Z \to \mu\mu\gamma$ process with bare leptons and bkg subtraction. $M_{ll}>45$ GeV. Nexp.un f. = 58310.9 $\pm$ 245, NPowHeg truth =564672, N Sherpa truth =814016.

Combined Covariance Matrix for Unfolded $M(l^{-}\gamma)$ distribution for $Z \to \mu\mu\gamma$ process with bare leptons and bkg subtraction. $M_{ll}>45$ GeV. Nexp.un f. = 58320.9 $\pm$ 245, NPowHeg truth =564672, N Sherpa truth =814016.

Combined Covariance Matrix for Unfolded dR distribution for $Z \to ee\gamma$ process with bare leptons and bkg subtraction. $M_{ll}>45$ GeV. Nexp.un f. = 57705.4 $\pm$ 254 , NPowHeg truth =301655, N Sherpa truth =763259.

Combined Covariance Matrix for Unfolded dR distribution for $Z \to \mu\mu\gamma$ process with bare leptons and bkg subtraction. $M_{ll}>45$ GeV. Nexp.un f. = 58869.2 $\pm$ 256, NPowHeg truth =564660, NS herpa truth =813995.

Combined Covariance Matrix for Unfolded $P_{T}^{\gamma}$ distribution for $Z \to ee\gamma$ process with bare leptons and bkg subtraction. $M_{ll}>45$ GeV. Nexp.un f. = 57358.6 $\pm$ 254 , NPowHeg truth =298662, N Sherpa truth =756609.

Combined Covariance Matrix for Unfolded $P_{T}^{\gamma}$ distribution for $Z \to \mu\mu\gamma$ process with bare leptons and bkg subtraction. $M_{ll}>45$ GeV. Nexp.un f. = 58226.4 $\pm$ 254, NPowHeg truth =558949, N Sherpa truth =806988.

Unfolded $M(l^{+}\gamma\gamma)$ distribution for $Z \to ll\gamma\gamma$ process with dressed leptons.

Unfolded $M(l^{+}\gamma\gamma)$ distribution for $Z \to ll\gamma\gamma$ process with bare leptons.

Unfolded dR distribution for $Z \to ll\gamma\gamma$ process with dressed leptons for leading photon.

Unfolded dR distribution for $Z \to ll\gamma\gamma$ process with bare leptons for leading photon.

Unfolded dR distribution for $Z \to ll\gamma\gamma$ process with dressed leptons for second photon.

Unfolded dR distribution for $Z \to ll\gamma\gamma$ process with bare leptons for second photon.

Unfolded dRgg distribution for $Z \to ll\gamma\gamma$ process with dressed leptons.

Unfolded dRgg distribution for $Z \to ll\gamma\gamma$ process with bare leptons.

Unfolded $P_{T}^{\gamma 1}$ distribution for $Z \to ll\gamma\gamma$ process with dressed leptons for leading photon.

Unfolded $P_{T}^{\gamma 1}$ distribution for $Z \to ll\gamma\gamma$ process with bare leptons for leading photon.

Unfolded $P_{T}^{\gamma 2}$ distribution for $Z \to ll\gamma\gamma$ process with dressed leptons for second photon.

Unfolded $P_{T}^{\gamma 2}$ distribution for $Z \to ll\gamma\gamma$ process with bare leptons for second photon.

The JEC, JER, and total systematic uncertainties in unfolded cross sections as functions of transverse momentum, for 1.5<|y|<2. The total systematic uncertainty includes also the luminosity, jet identification and trigger efficiency uncertainties.

The unfolded measured particle-level inclusive jet cross section as functions of jet pT (markers), for |y|<0.5, compared to the NLO perturbative QCD prediction (histogram), using the CT14NLO PDF set, with muR = muF = HT, and corrected for the NP effects. The experimental and theoretical systematic uncertainties are shown.

The unfolded measured particle-level inclusive jet cross section as functions of jet pT (markers), for 0.5<|y|<1, compared to the NLO perturbative QCD prediction (histogram), using the CT14NLO PDF set, with muR = muF = HT, and corrected for the NP effects. The experimental and theoretical systematic uncertainties are shown.

The unfolded measured particle-level inclusive jet cross section as functions of jet pT (markers), for 1<|y|<1.5, compared to the NLO perturbative QCD prediction (histogram), using the CT14NLO PDF set, with muR = muF = HT, and corrected for the NP effects. The experimental and theoretical systematic uncertainties are shown.

The unfolded measured particle-level inclusive jet cross section as functions of jet pT (markers), for 1.5<|y|<2, compared to the NLO perturbative QCD prediction (histogram), using the CT14NLO PDF set, with muR = muF = HT, and corrected for the NP effects. The experimental and theoretical systematic uncertainties are shown.

Ratios (points) of the unfolded measured cross sections to the NLO theoretical predictions, using the CT14NLO PDF set, with mu = pT, for |y|<0.5. The vertical error bars show the statistical experimental uncertainty. The systematic experimental uncertainty, the total theoretical uncertainty, and the individual sources of theoretical uncertainty are shown.

Ratios (points) of the unfolded measured cross sections to the NLO theoretical predictions, using the CT14NLO PDF set, with mu = pT, for 0.5<|y|<1. The vertical error bars show the statistical experimental uncertainty. The systematic experimental uncertainty, the total theoretical uncertainty, and the individual sources of theoretical uncertainty are shown.

Ratios (points) of the unfolded measured cross sections to the NLO theoretical predictions, using the CT14NLO PDF set, with mu = pT, for 1<|y|<1.5. The vertical error bars show the statistical experimental uncertainty. The systematic experimental uncertainty, the total theoretical uncertainty, and the individual sources of theoretical uncertainty are shown.

Ratios (points) of the unfolded measured cross sections to the NLO theoretical predictions, using the CT14NLO PDF set, with mu = pT, for 1.5<|y|<2. The vertical error bars show the statistical experimental uncertainty. The systematic experimental uncertainty, the total theoretical uncertainty, and the individual sources of theoretical uncertainty are shown.

Ratios (points) of the unfolded measured cross sections to the NLO theoretical predictions, using the CT14NLO PDF set, with mu = HT, for |y|<0.5. The vertical error bars show the statistical experimental uncertainty. The systematic experimental uncertainty, the total theoretical uncertainty, and the individual sources of theoretical uncertainty are shown.

Ratios (points) of the unfolded measured cross sections to the NLO theoretical predictions, using the CT14NLO PDF set, with mu = HT, for 0.5<|y|<1. The vertical error bars show the statistical experimental uncertainty. The systematic experimental uncertainty, the total theoretical uncertainty, and the individual sources of theoretical uncertainty are shown.

Ratios (points) of the unfolded measured cross sections to the NLO theoretical predictions, using the CT14NLO PDF set, with mu = HT, for 1<|y|<1.5. The vertical error bars show the statistical experimental uncertainty. The systematic experimental uncertainty, the total theoretical uncertainty, and the individual sources of theoretical uncertainty are shown.

Ratios (points) of the unfolded measured cross sections to the NLO theoretical predictions, using the CT14NLO PDF set, with mu = HT, for 1.5<|y|<2. The vertical error bars show the statistical experimental uncertainty. The systematic experimental uncertainty, the total theoretical uncertainty, and the individual sources of theoretical uncertainty are shown.

Ratios (points) of the unfolded measured cross sections to the NNLO theoretical predictions, using the CT14NNLO PDF set, with mu = HT, for |y|<0.5. The vertical error bars show the statistical experimental uncertainty. The systematic experimental uncertainty, the total theoretical uncertainty, and the individual sources of theoretical uncertainty are shown.

Ratios (points) of the unfolded measured cross sections to the NNLO theoretical predictions, using the CT14NNLO PDF set, with mu = HT, for 0.5<|y|<1. The vertical error bars show the statistical experimental uncertainty. The systematic experimental uncertainty, the total theoretical uncertainty, and the individual sources of theoretical uncertainty are shown.

Ratios (points) of the unfolded measured cross sections to the NNLO theoretical predictions, using the CT14NNLO PDF set, with mu = HT, for 1<|y|<1.5. The vertical error bars show the statistical experimental uncertainty. The systematic experimental uncertainty, the total theoretical uncertainty, and the individual sources of theoretical uncertainty are shown.

Ratios (points) of the unfolded measured cross sections to the NNLO theoretical predictions, using the CT14NNLO PDF set, with mu = HT, for 1.5<|y|<2. The vertical error bars show the statistical experimental uncertainty. The systematic experimental uncertainty, the total theoretical uncertainty, and the individual sources of theoretical uncertainty are shown.

Ratios (points) of the unfolded measured cross sections to the NNLO theoretical predictions, using the NNPDF31NNLO PDF set, with mu = HT, for |y|<0.5. The vertical error bars show the statistical experimental uncertainty. The systematic experimental uncertainty, the total theoretical uncertainty, and the individual sources of theoretical uncertainty are shown.

Ratios (points) of the unfolded measured cross sections to the NNLO theoretical predictions, using the NNPDF31NNLO PDF set, with mu = HT, for 0.5<|y|<1. The vertical error bars show the statistical experimental uncertainty. The systematic experimental uncertainty, the total theoretical uncertainty, and the individual sources of theoretical uncertainty are shown.

Ratios (points) of the unfolded measured cross sections to the NNLO theoretical predictions, using the NNPDF31NNLO PDF set, with mu = HT, for 1<|y|<1.5. The vertical error bars show the statistical experimental uncertainty. The systematic experimental uncertainty, the total theoretical uncertainty, and the individual sources of theoretical uncertainty are shown.

Ratios (points) of the unfolded measured cross sections to the NNLO theoretical predictions, using the NNPDF31NNLO PDF set, with mu = HT, for 1.5<|y|<2. The vertical error bars show the statistical experimental uncertainty. The systematic experimental uncertainty, the total theoretical uncertainty, and the individual sources of theoretical uncertainty are shown.

The effect of aS(MZ) variation, for |y|<0.5. The NNLO theoretical cross section predictions using the NNPDF31NNLO PDF with m = HT, calculated for different choices of aS (0.108, 0.110, 0.112, 0.114, 0.116, 0.117, 0.118, 0.119, 0.120, 0.122, and 0.124), are divided by the benchmark NNLO prediction for aS = 0.118 and the same choice of PDF set, muR, and muF. Also shown is the experimental unfolded measurement divided by the same benchmark prediction. The width of the unity line corresponds to the statistical uncertainty from the MC integration for the determination of the NNLO prediction. The error bars on the unfolded data correspond to the total experimental statistical and systematic uncertainty added in quadrature.

The effect of aS(MZ) variation, for 0.5<|y|<1. The NNLO theoretical cross section predictions using the NNPDF31NNLO PDF with m = HT, calculated for different choices of aS (0.108, 0.110, 0.112, 0.114, 0.116, 0.117, 0.118, 0.119, 0.120, 0.122, and 0.124), are divided by the benchmark NNLO prediction for aS = 0.118 and the same choice of PDF set, muR, and muF. Also shown is the experimental unfolded measurement divided by the same benchmark prediction. The width of the unity line corresponds to the statistical uncertainty from the MC integration for the determination of the NNLO prediction. The error bars on the unfolded data correspond to the total experimental statistical and systematic uncertainty added in quadrature.

The effect of aS(MZ) variation, for 1<|y|<1.5. The NNLO theoretical cross section predictions using the NNPDF31NNLO PDF with m = HT, calculated for different choices of aS (0.108, 0.110, 0.112, 0.114, 0.116, 0.117, 0.118, 0.119, 0.120, 0.122, and 0.124), are divided by the benchmark NNLO prediction for aS = 0.118 and the same choice of PDF set, muR, and muF. Also shown is the experimental unfolded measurement divided by the same benchmark prediction. The width of the unity line corresponds to the statistical uncertainty from the MC integration for the determination of the NNLO prediction. The error bars on the unfolded data correspond to the total experimental statistical and systematic uncertainty added in quadrature.

The effect of aS(MZ) variation, for 1.5<|y|<2. The NNLO theoretical cross section predictions using the NNPDF31NNLO PDF with m = HT, calculated for different choices of aS (0.108, 0.110, 0.112, 0.114, 0.116, 0.117, 0.118, 0.119, 0.120, 0.122, and 0.124), are divided by the benchmark NNLO prediction for aS = 0.118 and the same choice of PDF set, muR, and muF. Also shown is the experimental unfolded measurement divided by the same benchmark prediction. The width of the unity line corresponds to the statistical uncertainty from the MC integration for the determination of the NNLO prediction. The error bars on the unfolded data correspond to the total experimental statistical and systematic uncertainty added in quadrature.

Unfolded correlation matrix for |y|<0.5.

Unfolded correlation matrix for 0.5<|y|<1.

Unfolded correlation matrix for 1<|y|<1.5.

Unfolded correlation matrix for 1.5<|y|<2.

Summary of results for cross-section of non-prompt $\psi(2S)$ decaying to a muon pair for 13 TeV data in fb/GeV. Uncertainties are statistical and systematic, respectively.

Summary of results for non-prompt fraction of $J/\psi$ decaying to a muon pair as a percentage for 13 TeV data. Uncertainties are statistical and systematic, respectively.

Summary of results for non-prompt fraction of $\psi(2S)$ decaying to a muon pair as a percentage for 13 TeV data. Uncertainties are statistical and systematic, respectively.

Summary of results for prompt $\psi(2S)$ over $J/\psi$ ratio as a percentage for 13 TeV data. Uncertainties are statistical and systematic, respectively.

Summary of results for non-prompt $\psi(2S)$ over $J/\psi$ ratio as a percentage for 13 TeV data. Uncertainties are statistical and systematic, respectively.

Correction factors for an assumed angular dependence $\propto 1+\lambda_{\theta} \cos^2\theta $ in the helicity frame, for several values of the parameter $\lambda_{\theta}$. The correction factors were found to be the same (within 1% - 2%) for $J\psi$ and $\psi(2S)$ mesons, for prompt and non-prompt production mechanisms, and for the three rapidity intervals considered in this paper.

Observed excluded cross-section at 95% CL for the WZ-mediated simplified model of wino $\tilde{\chi}^{\pm}_{1}/\tilde{\chi}^{0}_{2}$ production from Fig 8(aux).

Expected exclusion limits at 95% CL for the WZ-mediated simplified model of wino $\tilde{\chi}^{\pm}_{1}/\tilde{\chi}^{0}_{2}$ production.

Observed exclusion limits at 95% CL for the Wh-mediated simplified model of wino $\tilde{\chi}^{\pm}_{1}/\tilde{\chi}^{0}_{2}$ production from from Fig 13(a) and from Fig 7 and Fig 10(aux).

Observed excluded cross-section at 95% CL for the Wh-mediated simplified model of wino $\tilde{\chi}^{\pm}_{1}/\tilde{\chi}^{0}_{2}$ production from Fig 7(aux) and Fig 10(aux).

positive one $\sigma$ observed exclusion limits at 95% CL for the Wh-mediated simplified model of wino $\tilde{\chi}^{\pm}_{1}/\tilde{\chi}^{0}_{2}$ production from from Fig 13(a) and from Fig 7 and Fig 10(aux).

negative one $\sigma$ observed exclusion limits at 95% CL for the Wh-mediated simplified model of wino $\tilde{\chi}^{\pm}_{1}/\tilde{\chi}^{0}_{2}$ production from from Fig 13(a) and from Fig 7 and Fig 10(aux).

Expected exclusion limits at 95% CL for the Wh-mediated simplified model of wino $\tilde{\chi}^{\pm}_{1}/\tilde{\chi}^{0}_{2}$ production.

Expected exclusion limits at 95% CL for the Wh-mediated simplified model of wino $\tilde{\chi}^{\pm}_{1}/\tilde{\chi}^{0}_{2}$ production.

Expected exclusion limits at 95% CL for the Wh-mediated simplified model of wino $\tilde{\chi}^{\pm}_{1}/\tilde{\chi}^{0}_{2}$ production.

Number of signal events expected for 139 fb$^{-1}$ at different stages of the event selection for the signal region $SR^{bRPV}_{2l-SS}$. in a susy scenario where $\tilde{\chi}^{0} _{1} and \tilde{\chi}^{0} _{2}$ are produced in pairs and decay to all possible allowed bRPV decays. The masses of the superpartners involved in the process are set to $m(\tilde{\chi}^{0} _{1}/\tilde{\chi}^{0} _{2})$ = 200 GeV, tan$\beta$=5. Only statistical uncertainties are shown.

Number of signal events expected for 139 fb$^{-1}$ at different stages of the event selection for the signal region $SR^{bRPV}_{3l}$. in a susy scenario where $\tilde{\chi}^{0} _{1} and \tilde{\chi}^{0} _{2}$ are produced in pairs and decay to all possible allowed bRPV decays. The masses of the superpartners involved in the process are set to $m(\tilde{\chi}^{0} _{1}/\tilde{\chi}^{0} _{2})$ = 200 GeV, tan$\beta$=5. Only statistical uncertainties are shown.

Number of signal events expected for 139 fb$^{-1}$ at different stages of the event selection for the signal region $SR^{WZ}_{high-m_{T2}}$. The wino-like doublet pair ($\tilde{\chi}^{\pm} _{1} and \tilde{\chi}^{0} _{2}$) were produced and then decays into $bino-like \tilde{\chi}^{0} _{1}$ which is the lightest SUSY particle (LSP) accompanied by mass on-shell or mass off-shell W and Z bosons. The masses of the superpartners involved in the process are set to $m(\tilde{\chi}^{\pm} _{1}/\tilde{\chi}^{0} _{2})$ = 150 GeV, $m(\tilde{\chi}^{0} _{1})$ = 50 GeV. Only statistical uncertainties are shown.

Number of signal events expected for 139 fb$^{-1}$ at different stages of the event selection for the signal region $SR^{WZ}_{low-m_{T2}}$. The wino-like doublet pair ($\tilde{\chi}^{\pm} _{1} and \tilde{\chi}^{0} _{2}$) were produced and then decays into $bino-like \tilde{\chi}^{0} _{1}$ which is the lightest SUSY particle (LSP) accompanied by mass on-shell or mass off-shell W and Z bosons. The masses of the superpartners involved in the process are set to $m(\tilde{\chi}^{\pm} _{1}/\tilde{\chi}^{0} _{2})$ = 150 GeV, $m(\tilde{\chi}^{0} _{1})$ = 50 GeV. Only statistical uncertainties are shown.

Number of signal events expected for 139 fb$^{-1}$ at different stages of the event selection for the low mass $SR^{RPV}_{2l1b}$, where the $\tilde{\chi}^{0} _{1} and \tilde{\chi}^{0} _{2}$ are directly produced and undergoes prompt RPV decays. The masses of the superpartners involved in the process are set to $m(\tilde{\chi}^{0} _{1}/\tilde{\chi}^{0} _{2})$ = 200 GeV. Only statistical uncertainties are shown.

Number of signal events expected for 139 fb$^{-1}$ at different stages of the event selection for the medium mass $SR^{RPV}_{2l1b}$, where the $\tilde{\chi}^{0} _{1} and \tilde{\chi}^{0} _{2}$ are directly produced and undergoes prompt RPV decays. The masses of the superpartners involved in the process are set to $m(\tilde{\chi}^{0} _{1}/\tilde{\chi}^{0} _{2})$ = 200 GeV. Only statistical uncertainties are shown.

Number of signal events expected for 139 fb$^{-1}$ at different stages of the event selection for the low mass $SR^{RPV}_{2l2b}$, where the $\tilde{\chi}^{0} _{1} and \tilde{\chi}^{0} _{2}$ are directly produced and undergoes prompt RPV decays. The masses of the superpartners involved in the process are set to $m(\tilde{\chi}^{0} _{1}/\tilde{\chi}^{0} _{2})$ = 200 GeV. Only statistical uncertainties are shown.

Number of signal events expected for 139 fb$^{-1}$ at different stages of the event selection for the medium mass $SR^{RPV}_{2l2b}$, where the $\tilde{\chi}^{0} _{1} and \tilde{\chi}^{0} _{2}$ are directly produced and undergoes prompt RPV decays. The masses of the superpartners involved in the process are set to $m(\tilde{\chi}^{0} _{1}/\tilde{\chi}^{0} _{2})$ = 200 GeV. Only statistical uncertainties are shown.

Number of signal events expected for 139 fb$^{-1}$ at different stages of the event selection for the high mass $SR^{RPV}_{2l2b}$, where the $\tilde{\chi}^{0} _{1} and \tilde{\chi}^{0} _{2}$ are directly produced and undergoes prompt RPV decays. The masses of the superpartners involved in the process are set to $m(\tilde{\chi}^{0} _{1}/\tilde{\chi}^{0} _{2})$ = 200 GeV. Only statistical uncertainties are shown.

Number of signal events expected for 139 fb$^{-1}$ at different stages of the event selection for the low mass $SR^{RPV}_{2l3b}$, where the $\tilde{\chi}^{0} _{1} and \tilde{\chi}^{0} _{2}$ are directly produced and undergoes prompt RPV decays. The masses of the superpartners involved in the process are set to $m(\tilde{\chi}^{0} _{1}/\tilde{\chi}^{0} _{2})$ = 200 GeV. Only statistical uncertainties are shown.

Number of signal events expected for 139 fb$^{-1}$ at different stages of the event selection for the medium mass $SR^{RPV}_{2l3b}$, where the $\tilde{\chi}^{0} _{1} and \tilde{\chi}^{0} _{2}$ are directly produced and undergoes prompt RPV decays. The masses of the superpartners involved in the process are set to $m(\tilde{\chi}^{0} _{1}/\tilde{\chi}^{0} _{2})$ = 200 GeV. Only statistical uncertainties are shown.

Number of signal events expected for 139 fb$^{-1}$ at different stages of the event selection for the high mass $SR^{RPV}_{2l3b}$, where the $\tilde{\chi}^{0} _{1} and \tilde{\chi}^{0} _{2}$ are directly produced and undergoes prompt RPV decays. The masses of the superpartners involved in the process are set to $m(\tilde{\chi}^{0} _{1}/\tilde{\chi}^{0} _{2})$ = 200 GeV. Only statistical uncertainties are shown.

Number of signal events expected for 139 fb$^{-1}$ at different stages of the event selection for the $SR^{Wh}_{low-m_{T2} }$. The wino-like doublet pair ($\tilde{\chi}^{\pm} _{1} and \tilde{\chi}^{0} _{2}$) were produced and then decays into $bino-like \tilde{\chi}^{0} _{1}$ which is the lightest SUSY particle (LSP) accompanied by mass on-shell or mass off-shell W and Higgs bosons. The masses of the superpartners involved in the process are set to $m(\tilde{\chi}^{\pm} _{1}/\tilde{\chi}^{0} _{2})$ = 300 GeV, $m(\tilde{\chi}^{0} _{1})$ = 100 GeV. Only statistical uncertainties are shown.

Number of signal events expected for 139 fb$^{-1}$ at different stages of the event selection for the $SR^{Wh}_{high-m_{T2} }$. The wino-like doublet pair ($\tilde{\chi}^{\pm} _{1} and \tilde{\chi}^{0} _{2}$) were produced and then decays into $bino-like \tilde{\chi}^{0} _{1}$ which is the lightest SUSY particle (LSP) accompanied by mass on-shell or mass off-shell W and Higgs bosons. The masses of the superpartners involved in the process are set to $m(\tilde{\chi}^{\pm} _{1}/\tilde{\chi}^{0} _{2})$ = 300 GeV, $m(\tilde{\chi}^{0} _{1})$ = 100 GeV. Only statistical uncertainties are shown.

Signal Hepdataeptance for $SR^{bRPV}_{2l-SS}$ signal region from Fig 13(a)(aux) in a SUSY scenario where $\tilde{\chi}^{0} _{1} and \tilde{\chi}^{0} _{2}$ are produced in pairs and decay to all possible allowed bRPV decays.

Signal Hepdataeptance for $SR^{bRPV}_{3l}$ signal region from Fig 13(b)(aux) in a SUSY scenario where $\tilde{\chi}^{0} _{1} and \tilde{\chi}^{0} _{2}$ are produced in pairs and decay to all possible allowed bRPV decays.

Signal acceptance for $SR^{WZ}_{high-m_{T2}}$ in a SUSY scenario where the wino-like doublet pair ($\tilde{\chi}^{\pm} _{1} and \tilde{\chi}^{0} _{2}$) were produced and then decays into $bino-like \tilde{\chi}^{0} _{1}$ which is the lightest SUSY particle (LSP) accompanied by mass on-shell or mass off-shell W and Z bosons.

Signal acceptance for $SR^{WZ}_{low-m_{T2}}$ in a SUSY scenario where the wino-like doublet pair ($\tilde{\chi}^{\pm} _{1} and \tilde{\chi}^{0} _{2}$) were produced and then decays into $bino-like \tilde{\chi}^{0} _{1}$ which is the lightest SUSY particle (LSP) accompanied by mass on-shell or mass off-shell W and Z bosons.

Signal acceptance for $SR^{RPV}_{2l1b}-L$ signal region in a SUSY scenario where the $\tilde{\chi}^{0} _{1} and \tilde{\chi}^{0} _{2}$ are directly produced and undergoes prompt RPV decays.

Signal acceptance for $SR^{RPV}_{2l1b}-M$ signal region in a SUSY scenario where the $\tilde{\chi}^{0} _{1} and \tilde{\chi}^{0} _{2}$ are directly produced and undergoes prompt RPV decays.

Signal acceptance for $SR^{RPV}_{2l2b}-L$ signal region in a SUSY scenario where the $\tilde{\chi}^{0} _{1} and \tilde{\chi}^{0} _{2}$ are directly produced and undergoes prompt RPV decays.

Signal acceptance for $SR^{RPV}_{2l2b}-M$ signal region in a SUSY scenario where the $\tilde{\chi}^{0} _{1} and \tilde{\chi}^{0} _{2}$ are directly produced and undergoes prompt RPV decays.

Signal acceptance for $SR^{RPV}_{2l2b}-H$ signal region in a SUSY scenario where the $\tilde{\chi}^{0} _{1} and \tilde{\chi}^{0} _{2}$ are directly produced and undergoes prompt RPV decays.

Signal acceptance for $SR^{RPV}_{2l3b}-L$ signal region in a SUSY scenario where the $\tilde{\chi}^{0} _{1} and \tilde{\chi}^{0} _{2}$ are directly produced and undergoes prompt RPV decays.

Signal acceptance for $SR^{RPV}_{2l3b}-M$ signal region in a SUSY scenario where the $\tilde{\chi}^{0} _{1} and \tilde{\chi}^{0} _{2}$ are directly produced and undergoes prompt RPV decays.

Signal acceptance for $SR^{RPV}_{2l3b}-H$ signal region in a SUSY scenario where the $\tilde{\chi}^{0} _{1} and \tilde{\chi}^{0} _{2}$ are directly produced and undergoes prompt RPV decays.

Signal efficiency for $SR^{bRPV}_{2l-SS}$ signal region in a SUSY scenario where $\tilde{\chi}^{0} _{1} and \tilde{\chi}^{0} _{2}$ are produced in pairs and decay to all possible allowed bRPV decays.

Signal efficiency for $SR^{bRPV}_{3l}$ signal region in a SUSY scenario where $\tilde{\chi}^{0} _{1} and \tilde{\chi}^{0} _{2}$ are produced in pairs and decay to all possible allowed bRPV decays.

Signal efficiency for $SR^{WZ}_{high-m_{T2}}$ in a SUSY scenario where the wino-like doublet pair ($\tilde{\chi}^{\pm} _{1} and \tilde{\chi}^{0} _{2}$) were produced and then decays into $bino-like \tilde{\chi}^{0} _{1}$ which is the lightest SUSY particle (LSP) accompanied by mass on-shell or mass off-shell W and Z bosons.

Signal efficiency for $SR^{WZ}_{low-m_{T2}}$ in a SUSY scenario where the wino-like doublet pair ($\tilde{\chi}^{\pm} _{1} and \tilde{\chi}^{0} _{2}$) were produced and then decays into $bino-like \tilde{\chi}^{0} _{1}$ which is the lightest SUSY particle (LSP) accompanied by mass on-shell or mass off-shell W and Z bosons.

Signal efficiency for $SR^{RPV}_{2l1b}-L$ signal region in a SUSY scenario where the $\tilde{\chi}^{0} _{1} and \tilde{\chi}^{0} _{2}$ are directly produced and undergoes prompt RPV decays.

Signal efficiency for $SR^{RPV}_{2l1b}-M$ signal region in a SUSY scenario where the $\tilde{\chi}^{0} _{1} and \tilde{\chi}^{0} _{2}$ are directly produced and undergoes prompt RPV decays.

Signal efficiency for $SR^{RPV}_{2l2b}-L$ signal region in a SUSY scenario where the $\tilde{\chi}^{0} _{1} and \tilde{\chi}^{0} _{2}$ are directly produced and undergoes prompt RPV decays.

Signal efficiency for $SR^{RPV}_{2l2b}-M$ signal region in a SUSY scenario where the $\tilde{\chi}^{0} _{1} and \tilde{\chi}^{0} _{2}$ are directly produced and undergoes prompt RPV decays.

Signal efficiency for $SR^{RPV}_{2l2b}-H$ signal region in a SUSY scenario where the $\tilde{\chi}^{0} _{1} and \tilde{\chi}^{0} _{2}$ are directly produced and undergoes prompt RPV decays.

Signal efficiency for $SR^{RPV}_{2l3b}-L$ signal region in a SUSY scenario where the $\tilde{\chi}^{0} _{1} and \tilde{\chi}^{0} _{2}$ are directly produced and undergoes prompt RPV decays.

Signal efficiency for $SR^{RPV}_{2l3b}-M$ signal region in a SUSY scenario where the $\tilde{\chi}^{0} _{1} and \tilde{\chi}^{0} _{2}$ are directly produced and undergoes prompt RPV decays.

Signal efficiency for $SR^{RPV}_{2l3b}-H$ signal region in a SUSY scenario where the $\tilde{\chi}^{0} _{1} and \tilde{\chi}^{0} _{2}$ are directly produced and undergoes prompt RPV decays.

Signal acceptance for $SR^{Wh}_{high-m_{T2} }$ signal region from Fig 11(a)(aux) in a SUSY scenario where direct production of a lightest $\tilde{\chi}^{\pm} _{1} and \tilde{\chi}^{0} _{2}$ , decay with 100% branching ratio to a final state with a same sign light lepton (e or $\mu$) pair and two lightest neutralino1, via the on-shell emission of SM W and Higgs bosons,

Signal acceptance for $SR^{Wh}_{low-m_{T2} }$ signal region from Fig 11(b)(aux) in a SUSY scenario where direct production of a lightest $\tilde{\chi}^{\pm} _{1} and \tilde{\chi}^{0} _{2}$ , decay with 100% branching ratio to a final state with a same sign light lepton (e or $\mu$) pair and two lightest neutralino1, via the on-shell emission of SM W and Higgs bosons,

Signal efficiency for $SR^{Wh}_{high-m_{T2} }$ signal region from Fig 15(a)(aux) in a SUSY scenario where direct production of a lightest $\tilde{\chi}^{\pm} _{1} and \tilde{\chi}^{0} _{2}$ , decay with 100% branching ratio to a final state with a same sign light lepton (e or $\mu$) pair and two lightest neutralino1, via the on-shell emission of SM W and Higgs bosons,

Signal efficiency for $SR^{Wh}_{low-m_{T2} }$ signal region from Fig 15(b)(aux) in a SUSY scenario where direct production of a lightest $\tilde{\chi}^{\pm} _{1} and \tilde{\chi}^{0} _{2}$ , decay with 100% branching ratio to a final state with a same sign light lepton (e or $\mu$) pair and two lightest neutralino1, via the on-shell emission of SM W and Higgs bosons,

Observed 95% X-section upper limits as a function of higgsino $\tilde{\chi}^{\pm}_{1}/\tilde{\chi}^{0}_{1}/\tilde{\chi}^{0}_{2}$ mass in the bilinear RPV model from Fig 14.

Observed 95% X-section upper limits as a function of higgsino $\tilde{\chi}^{0}_{1}/\tilde{\chi}^{0}_{2}$ mass in the UDD RPV model from Fig 18.

Observed 95% X-section upper limits as a function of wino $\tilde{\chi}^{\pm}_{1}/\tilde{\chi}^{0}_{2}$ mass in the WZ-mediated simplified model of wino $\tilde{\chi}^{\pm}_{1}/\tilde{\chi}^{0}_{2}$ production from Fig 9(aux).

N-1 distributions for $m_{T2}$ of observed data and expected background towards $SR^{WZ}_{high-m_{T2}}$ from publication's Figure 11(a) . The last bin is inclusive.

N-1 distributions for $m_{T2}$ of observed data and expected background towards $SR^{WZ}_{low-m_{T2}}$ from publication's Figure 11(b) . The last bin is inclusive.

N-1 distributions for $m_{T2}$ of observed data and expected background towards $SR^{bRPV}_{2l-SS}$ from publication's Figure 11(c) . The last bin is inclusive.

N-1 distributions for $m_{T2}$ of observed data and expected background towards $SR^{bRPV}_{3l}$ from publication's Figure 11(d) . The last bin is inclusive.

N-1 distributions for $\sum p^{b-jet}_{T}/\sum p^{jet}_{T}$ of observed data and expected background towards $SR^{RPV}_{2l1b}-L$ from publication's Figure 16(a) . The last bin is inclusive.

N-1 distributions for $\sum p^{b-jet}_{T}/\sum p^{jet}_{T}$ of observed data and expected background towards $SR^{RPV}_{2l2b}-M$ from publication's Figure 16(b) . The last bin is inclusive.

N-1 distributions for $\sum p^{b-jet}_{T}/\sum p^{jet}_{T}$ of observed data and expected background towards $SR^{RPV}_{2l3b}-H$ from publication's Figure 16(c) . The last bin is inclusive.

N-1 distribution for $E_{T}^{miss}$ in $SR^{Wh}_{high-m_{T2} }$ in ee channel

N-1 distribution for $E_{T}^{miss}$ in $SR^{Wh}_{high-m_{T2} }$ in e$\mu$ channel

N-1 distribution for $E_{T}^{miss}$ in $SR^{Wh}_{high-m_{T2} }$ in $\mu\mu$ channel

N-1 distribution for $\mathcal{S}(E_{T}^{miss})$ in $SR^{Wh}_{low-m_{T2} }$ in ee channel

N-1 distribution for $\mathcal{S}(E_{T}^{miss})$ in $SR^{Wh}_{low-m_{T2} }$ in e$\mu$ channel

N-1 distribution for $\mathcal{S}(E_{T}^{miss})$ in $SR^{Wh}_{low-m_{T2} }$ in $\mu\mu$ channel

Exclusion limits at 95% CL for direct gluino pair production, where the gluinos decay to light-flavor quarks and a $\tilde{\chi}_1^\pm$ that decays to $W^\pm\tilde{\chi}_1^0$. Signal cross sections are calculated at approximately NNLO+NNLL order in $\alpha_S$, assuming unity branching fraction to $q_i\bar{q}_j W\pm\tilde{\chi}_1^0$.

Exclusion limits at 95% CL for direct gluino pair production, where the gluinos decay to bottom quarks ($\tilde{g}\to b\bar{b}\tilde{\chi}_1^0$). Signal cross sections are calculated at approximately NNLO+NNLL order in $\alpha_S$, assuming unity branching fraction for the given decay.

Exclusion limits at 95% CL for direct gluino pair production, where the gluinos decay to top quarks ($\tilde{g}\to t\bar{t}\tilde{\chi}_1^0$). Signal cross sections are calculated at approximately NNLO+NNLL order in $\alpha_S$, assuming unity branching fraction for the given decay.

Exclusion limits at 95% CL for light-flavor squark pair production, where the squarks decay to $q\tilde{\chi}_1^0$. Signal cross sections are calculated at approximately NNLO+NNLL order in $\alpha_S$, assuming unity branching fraction for the given decay and 1-fold degeneracy in the light-flavor squarks (corresponding to the inner set of curves in the limit plot). To get the theory cross section for other N-fold degeneracy assumptions (e.g. 8-fold for the outer curves in the limit plot), just multiply by N.

Exclusion limits at 95% CL for bottom squark pair production, where the squarks decay to $b\tilde{\chi}_1^0$. Signal cross sections are calculated at approximately NNLO+NNLL order in $\alpha_S$, assuming unity branching fraction for the given decay.

Exclusion limits at 95% CL for top squark pair production, where the squarks decay to $t\tilde{\chi}_1^0$. Signal cross sections are calculated at approximately NNLO+NNLL order in $\alpha_S$, assuming unity branching fraction for the given decay.

Exclusion limits at 95% CL for top squark pair production, where the squarks decay to $b\tilde{\chi}_1^\pm$ and the $\tilde{\chi}_1^0$ decay to $W^\pm\tilde{\chi}_1^0$. Signal cross sections are calculated at approximately NNLO+NNLL order in $\alpha_S$, assuming unity branching fraction for the given decay.

Exclusion limits at 95% CL for top squark pair production, where the squarks decay either to $b\tilde{\chi}_1^\pm\to bW^\pm\tilde{\chi}_1^0$ or to $t\tilde{\chi}_1^0$. Signal cross sections are calculated at approximately NNLO+NNLL order in $\alpha_S$, assuming unity branching fraction for the given decay.

Exclusion limits at 95% CL for top squark pair production, where the squarks decay to $c\tilde{\chi}_1^0$. Signal cross sections are calculated at approximately NNLO+NNLL order in $\alpha_S$, assuming unity branching fraction for the given decay.

Exclusion limits at 95% CL for the mono-$\phi$ model, in which a resonantly-produced colored scalar decays to a massive Dirac fermion and a quark. Signal cross sections are calculated at leading order in $\alpha_S$, assuming unity branching fraction for the given decay.

Cross section limits for $\mathrm{LQ}\to\mathrm{q}\nu$, where $q=u,\,d,\,s,\,\mathrm{or}\,c$. Limits are at the 95% confidence level. Theory cross sections are LO for vector LQ, and NLO for scalar LQ. Branching ratio is assumed to be 100% to $\mathrm{q}\nu$.

Cross section limits for $\mathrm{LQ}\to\mathrm{b}\nu$. Limits are at the 95% confidence level. Theory cross sections are LO for vector LQ, and NLO for scalar LQ. Branching ratio is assumed to be 100% to $\mathrm{b}\nu$.

Cross section limits for $\mathrm{LQ}\to\mathrm{t}\nu$. Limits are at the 95% confidence level. Theory cross sections are LO for vector LQ, and NLO for scalar LQ. Branching ratios are assumed to be $\mathcal{B}(\mathrm{LQ}\to\mathrm{t}\nu)=1-\beta$, and $\mathcal{B}(\mathrm{LQ}\to\mathrm{b}\tau)=\beta$.

Predictions and observations for monojet signal regions

Predictions and observations for signal regions with $250 \leq H_\mathrm{T} < 450$ GeV

Predictions and observations for signal regions with $450 \leq H_\mathrm{T} < 575$ GeV and $N_\mathrm{j}<7$

Predictions and observations for signal regions with $450 \leq H_\mathrm{T} < 575$ GeV and $N_\mathrm{j}\geq7$

Predictions and observations for signal regions with $575 \leq H_\mathrm{T} < 1200$ GeV and $N_\mathrm{j}^\mathrm{hi}<4$

Predictions and observations for signal regions with $575 \leq H_\mathrm{T} < 1200$ GeV and $4\leq N_\mathrm{j}^\mathrm{hi}<7$

Predictions and observations for signal regions with $575 \leq H_\mathrm{T} < 1200$ GeV and $N_\mathrm{j}\geq7$

Predictions and observations for signal regions with $1200 \leq H_\mathrm{T} < 1500$ GeV and $N_\mathrm{j}^\mathrm{hi}<4$

Predictions and observations for signal regions with $1200 \leq H_\mathrm{T} < 1500$ GeV and $4\leq N_\mathrm{j}^\mathrm{hi}<7$

Predictions and observations for signal regions with $1200 \leq H_\mathrm{T} < 1500$ GeV and $N_\mathrm{j}\geq7$

Predictions and observations for signal regions with $H_\mathrm{T} \geq 1500$ GeV and $N_\mathrm{j}<7$

Predictions and observations for signal regions with $H_\mathrm{T} \geq 1500$ GeV and $N_\mathrm{j}\geq7$

Covariance matrix for the 282 signal regions of the inclusive $M_\mathrm{T2}$ search

Correlation matrix for the 282 signal regions of the inclusive $M_\mathrm{T2}$ search

Bin number definitions for the $M_\mathrm{T2}$ covariance and correlation matrices

Exclusion limits at 95% CL for direct gluino pair production, where the gluinos decay to light-flavor quarks and either the lightest neutralino, or the lightest chargino, and the chargino is long-lived with $c\tau_0 = 10$ cm and mass O(100) MeV greater than the neutralino's mass. Signal cross sections are calculated at approximately NNLO+NNLL order in $\alpha_S$.

Exclusion limits at 95% CL for direct gluino pair production, where the gluinos decay to light-flavor quarks and either the lightest neutralino, or the lightest chargino, and the chargino is long-lived with $c\tau_0 = 50$ cm and mass O(100) MeV greater than the neutralino's mass. Signal cross sections are calculated at approximately NNLO+NNLL order in $\alpha_S$.

Exclusion limits at 95% CL for direct gluino pair production, where the gluinos decay to light-flavor quarks and either the lightest neutralino, or the lightest chargino, and the chargino is long-lived with $c\tau_0 = 200$ cm and mass O(100) MeV greater than the neutralino's mass. Signal cross sections are calculated at approximately NNLO+NNLL order in $\alpha_S$.

Exclusion limits at 95% CL for direct light squark pair production, where the squarks decay to light-flavor quarks and either the lightest neutralino, or the lightest chargino, and the chargino is long-lived with $c\tau_0 = 10$ cm and mass O(100) MeV greater than the neutralino's mass. Signal cross sections are calculated at approximately NNLO+NNLL order in $\alpha_S$, for a single light squark.

Exclusion limits at 95% CL for direct light squark pair production, where the squarks decay to light-flavor quarks and either the lightest neutralino, or the lightest chargino, and the chargino is long-lived with $c\tau_0 = 50$ cm and mass O(100) MeV greater than the neutralino's mass. Signal cross sections are calculated at approximately NNLO+NNLL order in $\alpha_S$, for a single light squark.

Exclusion limits at 95% CL for direct light squark pair production, where the squarks decay to light-flavor quarks and either the lightest neutralino, or the lightest chargino, and the chargino is long-lived with $c\tau_0 = 200$ cm and mass O(100) MeV greater than the neutralino's mass. Signal cross sections are calculated at approximately NNLO+NNLL order in $\alpha_S$, for a single light squark.

Exclusion limits at 95% CL for direct stop pair production, where the stops decay to either a top and the lightest neutralino, or a bottom and the lightest chargino, and the chargino is long-lived with $c\tau_0 = 10$ cm and mass O(100) MeV greater than the neutralino's mass. Signal cross sections are calculated at approximately NNLO+NNLL order in $\alpha_S$.

Exclusion limits at 95% CL for direct stop pair production, where the stops decay to either a top and the lightest neutralino, or a bottom and the lightest chargino, and the chargino is long-lived with $c\tau_0 = 50$ cm and mass O(100) MeV greater than the neutralino's mass. Signal cross sections are calculated at approximately NNLO+NNLL order in $\alpha_S$.

Exclusion limits at 95% CL for direct stop pair production, where the stops decay to either a top and the lightest neutralino, or a bottom and the lightest chargino, and the chargino is long-lived with $c\tau_0 = 200$ cm and mass O(100) MeV greater than the neutralino's mass. Signal cross sections are calculated at approximately NNLO+NNLL order in $\alpha_S$.

The maximum chargino mass excluded at 95% CL for direct gluino pair production, where the gluinos decay to light-flavor quarks and either the lightest neutralino, or the lightest chargino, and the chargino mass is O(100) MeV greater than the neutralino's mass. The chargino's lifetime is varied from $c\tau_{0} = 1$ to 2000 cm while the gluino mass is fixed to 1900 GeV. Signal cross sections are calculated at approximately NNLO+NNLL order in $\alpha_S$. If all kinematically allowed chargino masses are excluded, the curves, including 68 and 95% expected, tend to overlap. At short decay lengths, horizontal exclusion lines are obtained from the inclusive analysis, as this is not affected by track reconstruction inefficiencies, which may arise when the chargino decays before the CMS tracker, and therefore shows better sensitivity to scenarios with very small lifetime compared to the disappearing track search, based on median expected limits.

The maximum chargino mass excluded at 95% CL for direct squark pair production, where the squarks decay to light-flavor quarks and either the lightest neutralino, or the lightest chargino, and the chargino mass is O(100) MeV greater than the neutralino's mass. The chargino's lifetime is varied from $c\tau_{0} = 1$ to 2000 cm while the squark mass is fixed to 900 GeV. Signal cross sections are calculated at approximately NNLO+NNLL order in $\alpha_S$, for a single light squark. If all kinematically allowed chargino masses are excluded, the curves, including 68 and 95% expected, tend to overlap. At short decay lengths, horizontal exclusion lines are obtained from the inclusive analysis, as this is not affected by track reconstruction inefficiencies, which may arise when the chargino decays before the CMS tracker, and therefore shows better sensitivity to scenarios with very small lifetime compared to the disappearing track search, based on median expected limits.

The maximum chargino mass excluded at 95% CL for direct squark pair production, where the squarks decay to light-flavor quarks and either the lightest neutralino, or the lightest chargino, and the chargino mass is O(100) MeV greater than the neutralino's mass. The chargino's lifetime is varied from $c\tau_{0} = 1$ to 2000 cm while the squark mass is fixed to 1500 GeV. Signal cross sections are calculated at approximately NNLO+NNLL order in $\alpha_S$, and the eight light squarks' masses are assumed to be degenerate. If all kinematically allowed chargino masses are excluded, the curves, including 68 and 95% expected, tend to overlap. At short decay lengths, horizontal exclusion lines are obtained from the inclusive analysis, as this is not affected by track reconstruction inefficiencies, which may arise when the chargino decays before the CMS tracker, and therefore shows better sensitivity to scenarios with very small lifetime compared to the disappearing track search, based on median expected limits.

Exclusion limits at 95% CL for direct stop pair production, where the stops decay to either a top and the lightest neutralino, or a bottom and the lightest chargino, and the chargino mass is O(100) MeV greater than the neutralino's mass. The chargino's lifetime is varied from $c\tau_{0} = 1$ to 2000 cm while the stop mass is fixed to 1000 GeV. Signal cross sections are calculated at approximately NNLO+NNLL order in $\alpha_S$. If all kinematically allowed chargino masses are excluded, the curves, including 68 and 95% expected, tend to overlap. At short decay lengths, horizontal exclusion lines are obtained from the inclusive analysis, as this is not affected by track reconstruction inefficiencies, which may arise when the chargino decays before the CMS tracker, and therefore shows better sensitivity to scenarios with very small lifetime compared to the disappearing track search, based on median expected limits.

Exclusion limits at 95% CL for direct gluino pair production, where the gluinos decay to light-flavor quarks and either the lightest neutralino, or the lightest chargino, and the chargino mass is O(100) MeV greater than the neutralino's mass. The chargino's lifetime is varied from $c\tau_{0} = 5$ to 1000 cm while the gluino mass is fixed to 1600 GeV and the neutralino's mass is fixed to 1575 GeV. Signal cross sections are calculated at approximately NNLO+NNLL order in $\alpha_S$.

Exclusion limits at 95% CL for direct squark pair production, where the squarks decay to light-flavor quarks and either the lightest neutralino, or the lightest chargino, and the chargino mass is O(100) MeV greater than the neutralino's mass. The chargino's lifetime is varied from $c\tau_{0} = 5$ to 1000 cm while the squark mass is fixed to 2000 GeV and the neutralino's mass is fixed to 1000 GeV. Signal cross sections are calculated at approximately NNLO+NNLL order in $\alpha_S$, and the eight light squarks' masses are assumed to be degenerate.

Exclusion limits at 95% CL for direct stop pair production, where the stops decay to either a top and the lightest neutralino, or a bottom and the lightest chargino, and the chargino mass is O(100) MeV greater than the neutralino's mass. The chargino's lifetime is varied from $c\tau_{0} = 5$ to 1000 cm while the stop mass is fixed to 1100 GeV and the neutralino's mass is fixed to 1000 GeV. Signal cross sections are calculated at approximately NNLO+NNLL order in $\alpha_S$.

Predictions and observations for 2016 disappearing track signal regions

Predictions and observations for 2017-2018 pixel track signal regions

Predictions and observations for 2017-2018 medium (M) and long (L) length track signal regions

Covariance matrix for the 68 signal regions of the disappearing tracks $M_\mathrm{T2}$ search

Correlation matrix for the 68 signal regions of the disappearing tracks $M_\mathrm{T2}$ search

Self-normalised yield of electrons from heavy-flavour hadron decays as a function of normalised charged-particle pseudorapidity density at midrapidity computed in pp collisions at $\sqrt{s}$ = 13 TeV in different $p_{\rm T}$ intervals

Self-normalised yield of electrons from heavy-flavour hadron decays as a function of normalised charged-particle pseudorapidity density at midrapidity computed in p--pPb collisions at $\sqrt{s_{\rm NN}}$ = 8.16 TeV in different $p_{\rm T}$ intervals

Integrated $D^0$ $R_{\rm AA}$ as a function of $N_{\mathrm{part}}$ in different $p_{\rm T}$ regions, 0$-$8 ${\mathrm{GeV/}}c$ (squares), 0.7$-$2.2 ${\mathrm{GeV/}}c$ (diamonds) and 3$-$8 ${\mathrm{GeV/}}c$ (circles). Open symbols are for the 0$-$80% minimum bias events. The vertical bar around unity denotes the overall normalization uncertainty from p+p reference.

Near-side ratios of V^{0}-triggered yields to charged-hadron triggered (h$-$h) ones, ${{Y}^{K^{0}_{S}-\mathrm{h}}_{\Delta\varphi}}/{{Y}^{\mathrm{h-h}}_{\Delta\varphi}}$ and ${{Y}^{(\Lambda+\overline{\Lambda})-\mathrm{h}}_{\Delta\varphi}}/{{Y}^{\mathrm{h-h}}_{\Delta\varphi}}$ in Pb-Pb collisions.

Away-side ratios of V^{0}-triggered yields to charged-hadron triggered (h$-$h) ones, ${{Y}^{K^{0}_{S}-\mathrm{h}}_{\Delta\varphi}}/{{Y}^{\mathrm{h-h}}_{\Delta\varphi}}$ and ${{Y}^{(\Lambda+\overline{\Lambda})-\mathrm{h}}_{\Delta\varphi}}/{{Y}^{\mathrm{h-h}}_{\Delta\varphi}}$ in Pb-Pb collisions.

Near-side ratios of V^{0}-triggered yields to charged-hadron triggered (h--h) ones, ${{Y}^{K^{0}_{S}-\mathrm{h}}_{\Delta\varphi}}/{{Y}^{\mathrm{h-h}}_{\Delta\varphi}}$ and ${{Y}^{(\Lambda+\overline{\Lambda})-\mathrm{h}}_{\Delta\varphi}}/{{Y}^{\mathrm{h-h}}_{\Delta\varphi}}$ in pp collisions.

Away-side ratios of V^{0}-triggered yields to charged-hadron triggered (h--h) ones, ${{Y}^{K^{0}_{S}-\mathrm{h}}_{\Delta\varphi}}/{{Y}^{\mathrm{h-h}}_{\Delta\varphi}}$ and ${{Y}^{(\Lambda+\overline{\Lambda})-\mathrm{h}}_{\Delta\varphi}}/{{Y}^{\mathrm{h-h}}_{\Delta\varphi}}$ in Pb-Pb collisions.

$p_{T}$-dependent nuclear modification factor of $\rm{K}^{*0}$ (average of particle and anti-particle) meson measured in 30-40% centrality class for Pb-Pb collisions at $\sqrt{s_{NN}}$ = 5.02 TeV.

$p_{T}$-dependent nuclear modification factor of $\rm{K}^{*0}$ (average of particle and anti-particle) meson measured in 40-50% centrality class for Pb-Pb collisions at $\sqrt{s_{NN}}$ = 5.02 TeV.

$p_{T}$-dependent nuclear modification factor of $\rm{K}^{*0}$ (average of particle and anti-particle) meson measured in 50-60% centrality class for Pb-Pb collisions at $\sqrt{s_{NN}}$ = 5.02 TeV.

$p_{T}$-dependent nuclear modification factor of $\rm{K}^{*0}$ (average of particle and anti-particle) meson measured in 60-70% centrality class for Pb-Pb collisions at $\sqrt{s_{NN}}$ = 5.02 TeV.

$p_{T}$-dependent nuclear modification factor of $\rm{K}^{*0}$ (average of particle and anti-particle) meson measured in 70-80% centrality class for Pb-Pb collisions at $\sqrt{s_{NN}}$ = 5.02 TeV.

$p_{T}$-dependent nuclear modification factor of $\rm{K}^{*0}$ (average of particle and anti-particle) meson measured in 20-40% centrality class for Pb-Pb collisions at $\sqrt{s_{NN}}$ = 5.02 TeV.

$p_{T}$-dependent nuclear modification factor of $\rm{K}^{*0}$ (average of particle and anti-particle) meson measured in 40-60% centrality class for Pb-Pb collisions at $\sqrt{s_{NN}}$ = 5.02 TeV.

$p_{T}$-dependent nuclear modification factor of $\rm{K}^{*0}$ (average of particle and anti-particle) meson measured in 60-80% centrality class for Pb-Pb collisions at $\sqrt{s_{NN}}$ = 5.02 TeV.

$p_{T}$-dependent nuclear modification factor of $\phi$ meson measured in 0-10% centrality class for Pb-Pb collisions at $\sqrt{s_{NN}}$ = 5.02 TeV.

$p_{T}$-dependent nuclear modification factor of $\phi$ meson measured in 10-20% centrality class for Pb-Pb collisions at $\sqrt{s_{NN}}$ = 5.02 TeV.

$p_{T}$-dependent nuclear modification factor of $\phi$ meson measured in 20-30% centrality class for Pb-Pb collisions at $\sqrt{s_{NN}}$ = 5.02 TeV.

$p_{T}$-dependent nuclear modification factor of $\phi$ meson measured in 30-40% centrality class for Pb-Pb collisions at $\sqrt{s_{NN}}$ = 5.02 TeV.

$p_{T}$-dependent nuclear modification factor of $\phi$ meson measured in 40-50% centrality class for Pb-Pb collisions at $\sqrt{s_{NN}}$ = 5.02 TeV.

$p_{T}$-dependent nuclear modification factor of $\phi$ meson measured in 50-60% centrality class for Pb-Pb collisions at $\sqrt{s_{NN}}$ = 5.02 TeV.

$p_{T}$-dependent nuclear modification factor of $\phi$ meson measured in 60-70% centrality class for Pb-Pb collisions at $\sqrt{s_{NN}}$ = 5.02 TeV.

$p_{T}$-dependent nuclear modification factor of $\phi$ meson measured in 70-80% centrality class for Pb-Pb collisions at $\sqrt{s_{NN}}$ = 5.02 TeV.

$p_{T}$-dependent nuclear modification factor of $\phi$ meson measured in 20-40% centrality class for Pb-Pb collisions at $\sqrt{s_{NN}}$ = 5.02 TeV.

$p_{T}$-dependent nuclear modification factor of $\phi$ meson measured in 40-60% centrality class for Pb-Pb collisions at $\sqrt{s_{NN}}$ = 5.02 TeV.

$p_{T}$-dependent nuclear modification factor of $\phi$ meson measured in 60-80% centrality class for Pb-Pb collisions at $\sqrt{s_{NN}}$ = 5.02 TeV.

Differential polarized inclusive jet cross sections as a function of jet transverse energy.

Differential polarized inclusive jet cross sections as a function of photon virtuality, Q**2.

Differential polarized inclusive jet cross sections as a function of photon virtuality, Q**2.

Differential polarized inclusive jet cross sections as a function of X. Update: contrary to the publication, this is D(SIG)/DLOG10(X) not D(SIG)/DX.

Differential polarized inclusive jet cross sections as a function of X. Update: contrary to the publication, this is D(SIG)/DLOG10(X) not D(SIG)/DX.

Integrated polarized inclusive-jet cross section in the given kinematical range for the E- beam.

Integrated polarized inclusive-jet cross section in the given kinematical range for the E+ beam.

Differential unpolarized cross section for single jet production as a function of the jet pseudorapidity.

Differential unpolarized cross section for two jet production as a function of the mean jet pseudorapidity.

Differential unpolarized cross section for three jet production as a function of the mean jet pseudorapidity.

Differential unpolarized cross section for single jet production as a function of the jet transverse energy.

Differential unpolarized cross section for two jet production as a function of the mean jet transverse energy.

Differential unpolarized cross section for three jet production as a function of the mean jet transverse energy.

Differential unpolarized cross section for single jet production as a function of photon virtuality, Q**2.

Differential unpolarized cross section for two jet production as a function of photon virtuality, Q**2.

Differential unpolarized cross section for three jet production as a function of photon virtuality, Q**2.

Differential unpolarized inclusive-jet cross section as a function of X. Update: contrary to the publication, this is D(SIG)/DLOG10(X) not D(SIG)/DX.

Integrated unpolarized jet cross section in the given kinemetical region.

Differential unpolarized dijet cross section as a function of the dijet mass.

Differential unpolarized three-jet cross section as a function of the three-jet mass.

Charged-hadron spectrum in the centrality interval 5-10% for p+Pb, divided by &#9001;TPPB&#9002;. The systematic uncertainties are described in the section 7 of the paper. The total systematic uncertainties are determined by adding the contributions from all relevant sources in quadrature.

Charged-hadron spectrum in the centrality interval 10-20% for p+Pb, divided by &#9001;TPPB&#9002;. The systematic uncertainties are described in the section 7 of the paper. The total systematic uncertainties are determined by adding the contributions from all relevant sources in quadrature.

Charged-hadron spectrum in the centrality interval 20-30% for p+Pb, divided by &#9001;TPPB&#9002;. The systematic uncertainties are described in the section 7 of the paper. The total systematic uncertainties are determined by adding the contributions from all relevant sources in quadrature.

Charged-hadron spectrum in the centrality interval 30-40% for p+Pb, divided by &#9001;TPPB&#9002;. The systematic uncertainties are described in the section 7 of the paper. The total systematic uncertainties are determined by adding the contributions from all relevant sources in quadrature.

Charged-hadron spectrum in the centrality interval 40-60% for p+Pb, divided by &#9001;TPPB&#9002;. The systematic uncertainties are described in the section 7 of the paper. The total systematic uncertainties are determined by adding the contributions from all relevant sources in quadrature.

Charged-hadron spectrum in the centrality interval 60-90% for p+Pb, divided by &#9001;TPPB&#9002;. The systematic uncertainties are described in the section 7 of the paper. The total systematic uncertainties are determined by adding the contributions from all relevant sources in quadrature.

Charged-hadron spectrum in the centrality interval 0-90% for p+Pb, divided by &#9001;TPPB&#9002;. The systematic uncertainties are described in the section 7 of the paper. The total systematic uncertainties are determined by adding the contributions from all relevant sources in quadrature.

Charged-hadron cross-section in pp collisions. The systematic uncertainties are described in the section 7 of the paper. The total systematic uncertainties are determined by adding the contributions from all relevant sources in quadrature.

Charged-hadron spectrum in the centrality interval 0-5% for Pb+Pb, divided by &#9001;TAA&#9002;. The systematic uncertainties are described in the section 7 of the paper. The total systematic uncertainties are determined by adding the contributions from all relevant sources in quadrature. The systematic uncertainty on momentum bias is negligible at low pT; in such cases, it is omitted in the table below.

Charged-hadron spectrum in the centrality interval 5-10% for Pb+Pb, divided by &#9001;TAA&#9002;. The systematic uncertainties are described in the section 7 of the paper. The total systematic uncertainties are determined by adding the contributions from all relevant sources in quadrature. The systematic uncertainty on momentum bias is negligible at low pT; in such cases, it is omitted in the table below.

Charged-hadron spectrum in the centrality interval 10-20% for Pb+Pb, divided by &#9001;TAA&#9002;. The systematic uncertainties are described in the section 7 of the paper. The total systematic uncertainties are determined by adding the contributions from all relevant sources in quadrature. The systematic uncertainty on momentum bias is negligible at low pT; in such cases, it is omitted in the table below.

Charged-hadron spectrum in the centrality interval 20-30% for Pb+Pb, divided by &#9001;TAA&#9002;. The systematic uncertainties are described in the section 7 of the paper. The total systematic uncertainties are determined by adding the contributions from all relevant sources in quadrature. The systematic uncertainty on momentum bias is negligible at low pT; in such cases, it is omitted in the table below.

Charged-hadron spectrum in the centrality interval 30-40% for Pb+Pb, divided by &#9001;TAA&#9002;. The systematic uncertainties are described in the section 7 of the paper. The total systematic uncertainties are determined by adding the contributions from all relevant sources in quadrature. The systematic uncertainty on momentum bias is negligible at low pT; in such cases, it is omitted in the table below.

Charged-hadron spectrum in the centrality interval 40-50% for Pb+Pb, divided by &#9001;TAA&#9002;. The systematic uncertainties are described in the section 7 of the paper. The total systematic uncertainties are determined by adding the contributions from all relevant sources in quadrature. The systematic uncertainty on momentum bias is negligible at low pT; in such cases, it is omitted in the table below.

Charged-hadron spectrum in the centrality interval 50-60% for Pb+Pb, divided by &#9001;TAA&#9002;. The systematic uncertainties are described in the section 7 of the paper. The total systematic uncertainties are determined by adding the contributions from all relevant sources in quadrature. The systematic uncertainty on momentum bias is negligible at low pT; in such cases, it is omitted in the table below.

Charged-hadron spectrum in the centrality interval 60-80% for Pb+Pb, divided by &#9001;TAA&#9002;. The systematic uncertainties are described in the section 7 of the paper. The total systematic uncertainties are determined by adding the contributions from all relevant sources in quadrature. The systematic uncertainty on momentum bias is negligible at low pT; in such cases, it is omitted in the table below.

Charged-hadron cross-section in pp collisions. The systematic uncertainties are described in the section 7 of the paper. The total systematic uncertainties are determined by adding the contributions from all relevant sources in quadrature.

Charged-hadron spectrum in the centrality interval 0-5% for Xe+Xe, divided by &#9001;TAA&#9002;. The systematic uncertainties are described in the section 7 of the paper. The total systematic uncertainties are determined by adding the contributions from all relevant sources in quadrature.

Charged-hadron spectrum in the centrality interval 5-10% for Xe+Xe, divided by &#9001;TAA&#9002;. The systematic uncertainties are described in the section 7 of the paper. The total systematic uncertainties are determined by adding the contributions from all relevant sources in quadrature.

Charged-hadron spectrum in the centrality interval 10-20% for Xe+Xe, divided by &#9001;TAA&#9002;. The systematic uncertainties are described in the section 7 of the paper. The total systematic uncertainties are determined by adding the contributions from all relevant sources in quadrature.

Charged-hadron spectrum in the centrality interval 20-30% for Xe+Xe, divided by &#9001;TAA&#9002;. The systematic uncertainties are described in the section 7 of the paper. The total systematic uncertainties are determined by adding the contributions from all relevant sources in quadrature.

Charged-hadron spectrum in the centrality interval 30-40% for Xe+Xe, divided by &#9001;TAA&#9002;. The systematic uncertainties are described in the section 7 of the paper. The total systematic uncertainties are determined by adding the contributions from all relevant sources in quadrature.

Charged-hadron spectrum in the centrality interval 40-50% for Xe+Xe, divided by &#9001;TAA&#9002;. The systematic uncertainties are described in the section 7 of the paper. The total systematic uncertainties are determined by adding the contributions from all relevant sources in quadrature.

Charged-hadron spectrum in the centrality interval 50-60% for Xe+Xe, divided by &#9001;TAA&#9002;. The systematic uncertainties are described in the section 7 of the paper. The total systematic uncertainties are determined by adding the contributions from all relevant sources in quadrature.

Charged-hadron spectrum in the centrality interval 60-80% for Xe+Xe, divided by &#9001;TAA&#9002;. The systematic uncertainties are described in the section 7 of the paper. The total systematic uncertainties are determined by adding the contributions from all relevant sources in quadrature.

Nuclear modification factor in centrality interval 0-5% for p+Pb. The systematic uncertainties are described in the section 7 of the paper. The total systematic uncertainties are determined by adding the contributions from all relevant sources in quadrature.

Nuclear modification factor in centrality interval 5-10% for p+Pb. The systematic uncertainties are described in the section 7 of the paper. The total systematic uncertainties are determined by adding the contributions from all relevant sources in quadrature.

Nuclear modification factor in centrality interval 10-20% for p+Pb. The systematic uncertainties are described in the section 7 of the paper. The total systematic uncertainties are determined by adding the contributions from all relevant sources in quadrature.

Nuclear modification factor in centrality interval 20-30% for p+Pb. The systematic uncertainties are described in the section 7 of the paper. The total systematic uncertainties are determined by adding the contributions from all relevant sources in quadrature.

Nuclear modification factor in centrality interval 30-40% for p+Pb. The systematic uncertainties are described in the section 7 of the paper. The total systematic uncertainties are determined by adding the contributions from all relevant sources in quadrature.

Nuclear modification factor in centrality interval 40-60% for p+Pb. The systematic uncertainties are described in the section 7 of the paper. The total systematic uncertainties are determined by adding the contributions from all relevant sources in quadrature.

Nuclear modification factor in centrality interval 60-90% for p+Pb. The systematic uncertainties are described in the section 7 of the paper. The total systematic uncertainties are determined by adding the contributions from all relevant sources in quadrature.

Nuclear modification factor in centrality interval 0-90% for p+Pb. The systematic uncertainties are described in the section 7 of the paper. The total systematic uncertainties are determined by adding the contributions from all relevant sources in quadrature.

Nuclear modification factor in centrality interval 0-5% for Pb+Pb. The systematic uncertainties are described in the section 7 of the paper. The total systematic uncertainties are determined by adding the contributions from all relevant sources in quadrature. The systematic uncertainty on momentum bias is negligible at low pT; in such cases, it is omitted in the table below.

Nuclear modification factor in centrality interval 5-10% for Pb+Pb. The systematic uncertainties are described in the section 7 of the paper. The total systematic uncertainties are determined by adding the contributions from all relevant sources in quadrature. The systematic uncertainty on momentum bias is negligible at low pT; in such cases, it is omitted in the table below.

Nuclear modification factor in centrality interval 10-20% for Pb+Pb. The systematic uncertainties are described in the section 7 of the paper. The total systematic uncertainties are determined by adding the contributions from all relevant sources in quadrature. The systematic uncertainty on momentum bias is negligible at low pT; in such cases, it is omitted in the table below.

Nuclear modification factor in centrality interval 20-30% for Pb+Pb. The systematic uncertainties are described in the section 7 of the paper. The total systematic uncertainties are determined by adding the contributions from all relevant sources in quadrature. The systematic uncertainty on momentum bias is negligible at low pT; in such cases, it is omitted in the table below.

Nuclear modification factor in centrality interval 30-40% for Pb+Pb. The systematic uncertainties are described in the section 7 of the paper. The total systematic uncertainties are determined by adding the contributions from all relevant sources in quadrature. The systematic uncertainty on momentum bias is negligible at low pT; in such cases, it is omitted in the table below.

Nuclear modification factor in centrality interval 40-50% for Pb+Pb. The systematic uncertainties are described in the section 7 of the paper. The total systematic uncertainties are determined by adding the contributions from all relevant sources in quadrature. The systematic uncertainty on momentum bias is negligible at low pT; in such cases, it is omitted in the table below.

Nuclear modification factor in centrality interval 50-60% for Pb+Pb. The systematic uncertainties are described in the section 7 of the paper. The total systematic uncertainties are determined by adding the contributions from all relevant sources in quadrature. The systematic uncertainty on momentum bias is negligible at low pT; in such cases, it is omitted in the table below.

Nuclear modification factor in centrality interval 60-80% for Pb+Pb. The systematic uncertainties are described in the section 7 of the paper. The total systematic uncertainties are determined by adding the contributions from all relevant sources in quadrature. The systematic uncertainty on momentum bias is negligible at low pT; in such cases, it is omitted in the table below.

Nuclear modification factor in centrality interval 0-5% for Xe+Xe. The systematic uncertainties are described in the section 7 of the paper. The total systematic uncertainties are determined by adding the contributions from all relevant sources in quadrature.

Nuclear modification factor in centrality interval 5-10% for Xe+Xe. The systematic uncertainties are described in the section 7 of the paper. The total systematic uncertainties are determined by adding the contributions from all relevant sources in quadrature.

Nuclear modification factor in centrality interval 10-20% for Xe+Xe. The systematic uncertainties are described in the section 7 of the paper. The total systematic uncertainties are determined by adding the contributions from all relevant sources in quadrature.

Nuclear modification factor in centrality interval 20-30% for Xe+Xe. The systematic uncertainties are described in the section 7 of the paper. The total systematic uncertainties are determined by adding the contributions from all relevant sources in quadrature.

Nuclear modification factor in centrality interval 30-40% for Xe+Xe. The systematic uncertainties are described in the section 7 of the paper. The total systematic uncertainties are determined by adding the contributions from all relevant sources in quadrature.

Nuclear modification factor in centrality interval 40-50% for Xe+Xe. The systematic uncertainties are described in the section 7 of the paper. The total systematic uncertainties are determined by adding the contributions from all relevant sources in quadrature.

Nuclear modification factor in centrality interval 50-60% for Xe+Xe. The systematic uncertainties are described in the section 7 of the paper. The total systematic uncertainties are determined by adding the contributions from all relevant sources in quadrature.

Nuclear modification factor in centrality interval 60-80% for Xe+Xe. The systematic uncertainties are described in the section 7 of the paper. The total systematic uncertainties are determined by adding the contributions from all relevant sources in quadrature.

Nuclear modification factor in centrality interval 0-5% for p+Pb. The systematic uncertainties are described in the section 7 of the paper. The total systematic uncertainties are determined by adding the contributions from all relevant sources in quadrature.

Nuclear modification factor in centrality interval 0-5% for p+Pb. The systematic uncertainties are described in the section 7 of the paper. The total systematic uncertainties are determined by adding the contributions from all relevant sources in quadrature.

Nuclear modification factor in centrality interval 0-5% for p+Pb. The systematic uncertainties are described in the section 7 of the paper. The total systematic uncertainties are determined by adding the contributions from all relevant sources in quadrature.

Nuclear modification factor in centrality interval 0-5% for p+Pb. The systematic uncertainties are described in the section 7 of the paper. The total systematic uncertainties are determined by adding the contributions from all relevant sources in quadrature.

Nuclear modification factor in centrality interval 5-10% for p+Pb. The systematic uncertainties are described in the section 7 of the paper. The total systematic uncertainties are determined by adding the contributions from all relevant sources in quadrature.

Nuclear modification factor in centrality interval 5-10% for p+Pb. The systematic uncertainties are described in the section 7 of the paper. The total systematic uncertainties are determined by adding the contributions from all relevant sources in quadrature.

Nuclear modification factor in centrality interval 5-10% for p+Pb. The systematic uncertainties are described in the section 7 of the paper. The total systematic uncertainties are determined by adding the contributions from all relevant sources in quadrature.

Nuclear modification factor in centrality interval 5-10% for p+Pb. The systematic uncertainties are described in the section 7 of the paper. The total systematic uncertainties are determined by adding the contributions from all relevant sources in quadrature.

Nuclear modification factor in centrality interval 10-20% for p+Pb. The systematic uncertainties are described in the section 7 of the paper. The total systematic uncertainties are determined by adding the contributions from all relevant sources in quadrature.

Nuclear modification factor in centrality interval 10-20% for p+Pb. The systematic uncertainties are described in the section 7 of the paper. The total systematic uncertainties are determined by adding the contributions from all relevant sources in quadrature.

Nuclear modification factor in centrality interval 10-20% for p+Pb. The systematic uncertainties are described in the section 7 of the paper. The total systematic uncertainties are determined by adding the contributions from all relevant sources in quadrature.

Nuclear modification factor in centrality interval 10-20% for p+Pb. The systematic uncertainties are described in the section 7 of the paper. The total systematic uncertainties are determined by adding the contributions from all relevant sources in quadrature.

Nuclear modification factor in centrality interval 20-30% for p+Pb. The systematic uncertainties are described in the section 7 of the paper. The total systematic uncertainties are determined by adding the contributions from all relevant sources in quadrature.

Nuclear modification factor in centrality interval 20-30% for p+Pb. The systematic uncertainties are described in the section 7 of the paper. The total systematic uncertainties are determined by adding the contributions from all relevant sources in quadrature.

Nuclear modification factor in centrality interval 20-30% for p+Pb. The systematic uncertainties are described in the section 7 of the paper. The total systematic uncertainties are determined by adding the contributions from all relevant sources in quadrature.

Nuclear modification factor in centrality interval 20-30% for p+Pb. The systematic uncertainties are described in the section 7 of the paper. The total systematic uncertainties are determined by adding the contributions from all relevant sources in quadrature.

Nuclear modification factor in centrality interval 30-40% for p+Pb. The systematic uncertainties are described in the section 7 of the paper. The total systematic uncertainties are determined by adding the contributions from all relevant sources in quadrature.

Nuclear modification factor in centrality interval 30-40% for p+Pb. The systematic uncertainties are described in the section 7 of the paper. The total systematic uncertainties are determined by adding the contributions from all relevant sources in quadrature.

Nuclear modification factor in centrality interval 30-40% for p+Pb. The systematic uncertainties are described in the section 7 of the paper. The total systematic uncertainties are determined by adding the contributions from all relevant sources in quadrature.

Nuclear modification factor in centrality interval 30-40% for p+Pb. The systematic uncertainties are described in the section 7 of the paper. The total systematic uncertainties are determined by adding the contributions from all relevant sources in quadrature.

Nuclear modification factor in centrality interval 40-60% for p+Pb. The systematic uncertainties are described in the section 7 of the paper. The total systematic uncertainties are determined by adding the contributions from all relevant sources in quadrature.

Nuclear modification factor in centrality interval 40-60% for p+Pb. The systematic uncertainties are described in the section 7 of the paper. The total systematic uncertainties are determined by adding the contributions from all relevant sources in quadrature.

Nuclear modification factor in centrality interval 40-60% for p+Pb. The systematic uncertainties are described in the section 7 of the paper. The total systematic uncertainties are determined by adding the contributions from all relevant sources in quadrature.

Nuclear modification factor in centrality interval 40-60% for p+Pb. The systematic uncertainties are described in the section 7 of the paper. The total systematic uncertainties are determined by adding the contributions from all relevant sources in quadrature.

Nuclear modification factor in centrality interval 60-90% for p+Pb. The systematic uncertainties are described in the section 7 of the paper. The total systematic uncertainties are determined by adding the contributions from all relevant sources in quadrature.

Nuclear modification factor in centrality interval 60-90% for p+Pb. The systematic uncertainties are described in the section 7 of the paper. The total systematic uncertainties are determined by adding the contributions from all relevant sources in quadrature.

Nuclear modification factor in centrality interval 60-90% for p+Pb. The systematic uncertainties are described in the section 7 of the paper. The total systematic uncertainties are determined by adding the contributions from all relevant sources in quadrature.

Nuclear modification factor in centrality interval 60-90% for p+Pb. The systematic uncertainties are described in the section 7 of the paper. The total systematic uncertainties are determined by adding the contributions from all relevant sources in quadrature.

Nuclear modification factor in centrality interval 0-90% for p+Pb. The systematic uncertainties are described in the section 7 of the paper. The total systematic uncertainties are determined by adding the contributions from all relevant sources in quadrature.

Nuclear modification factor in centrality interval 0-90% for p+Pb. The systematic uncertainties are described in the section 7 of the paper. The total systematic uncertainties are determined by adding the contributions from all relevant sources in quadrature.

Nuclear modification factor in centrality interval 0-90% for p+Pb. The systematic uncertainties are described in the section 7 of the paper. The total systematic uncertainties are determined by adding the contributions from all relevant sources in quadrature.

Nuclear modification factor in centrality interval 0-90% for p+Pb. The systematic uncertainties are described in the section 7 of the paper. The total systematic uncertainties are determined by adding the contributions from all relevant sources in quadrature.

Nuclear modification factor in centrality interval 0-5% for Pb+Pb. The systematic uncertainties are described in the section 7 of the paper. The total systematic uncertainties are determined by adding the contributions from all relevant sources in quadrature.

Nuclear modification factor in centrality interval 0-5% for Pb+Pb. The systematic uncertainties are described in the section 7 of the paper. The total systematic uncertainties are determined by adding the contributions from all relevant sources in quadrature.

Nuclear modification factor in centrality interval 0-5% for Pb+Pb. The systematic uncertainties are described in the section 7 of the paper. The total systematic uncertainties are determined by adding the contributions from all relevant sources in quadrature.

Nuclear modification factor in centrality interval 0-5% for Pb+Pb. The systematic uncertainties are described in the section 7 of the paper. The total systematic uncertainties are determined by adding the contributions from all relevant sources in quadrature.

Nuclear modification factor in centrality interval 5-10% for Pb+Pb. The systematic uncertainties are described in the section 7 of the paper. The total systematic uncertainties are determined by adding the contributions from all relevant sources in quadrature.

Nuclear modification factor in centrality interval 5-10% for Pb+Pb. The systematic uncertainties are described in the section 7 of the paper. The total systematic uncertainties are determined by adding the contributions from all relevant sources in quadrature.

Nuclear modification factor in centrality interval 5-10% for Pb+Pb. The systematic uncertainties are described in the section 7 of the paper. The total systematic uncertainties are determined by adding the contributions from all relevant sources in quadrature.

Nuclear modification factor in centrality interval 5-10% for Pb+Pb. The systematic uncertainties are described in the section 7 of the paper. The total systematic uncertainties are determined by adding the contributions from all relevant sources in quadrature.

Nuclear modification factor in centrality interval 10-20% for Pb+Pb. The systematic uncertainties are described in the section 7 of the paper. The total systematic uncertainties are determined by adding the contributions from all relevant sources in quadrature.

Nuclear modification factor in centrality interval 10-20% for Pb+Pb. The systematic uncertainties are described in the section 7 of the paper. The total systematic uncertainties are determined by adding the contributions from all relevant sources in quadrature.

Nuclear modification factor in centrality interval 10-20% for Pb+Pb. The systematic uncertainties are described in the section 7 of the paper. The total systematic uncertainties are determined by adding the contributions from all relevant sources in quadrature.

Nuclear modification factor in centrality interval 10-20% for Pb+Pb. The systematic uncertainties are described in the section 7 of the paper. The total systematic uncertainties are determined by adding the contributions from all relevant sources in quadrature.

Nuclear modification factor in centrality interval 20-30% for Pb+Pb. The systematic uncertainties are described in the section 7 of the paper. The total systematic uncertainties are determined by adding the contributions from all relevant sources in quadrature.

Nuclear modification factor in centrality interval 20-30% for Pb+Pb. The systematic uncertainties are described in the section 7 of the paper. The total systematic uncertainties are determined by adding the contributions from all relevant sources in quadrature.

Nuclear modification factor in centrality interval 20-30% for Pb+Pb. The systematic uncertainties are described in the section 7 of the paper. The total systematic uncertainties are determined by adding the contributions from all relevant sources in quadrature.

Nuclear modification factor in centrality interval 20-30% for Pb+Pb. The systematic uncertainties are described in the section 7 of the paper. The total systematic uncertainties are determined by adding the contributions from all relevant sources in quadrature.

Nuclear modification factor in centrality interval 30-40% for Pb+Pb. The systematic uncertainties are described in the section 7 of the paper. The total systematic uncertainties are determined by adding the contributions from all relevant sources in quadrature.

Nuclear modification factor in centrality interval 30-40% for Pb+Pb. The systematic uncertainties are described in the section 7 of the paper. The total systematic uncertainties are determined by adding the contributions from all relevant sources in quadrature.

Nuclear modification factor in centrality interval 30-40% for Pb+Pb. The systematic uncertainties are described in the section 7 of the paper. The total systematic uncertainties are determined by adding the contributions from all relevant sources in quadrature.

Nuclear modification factor in centrality interval 30-40% for Pb+Pb. The systematic uncertainties are described in the section 7 of the paper. The total systematic uncertainties are determined by adding the contributions from all relevant sources in quadrature.

Nuclear modification factor in centrality interval 40-50% for Pb+Pb. The systematic uncertainties are described in the section 7 of the paper. The total systematic uncertainties are determined by adding the contributions from all relevant sources in quadrature.

Nuclear modification factor in centrality interval 40-50% for Pb+Pb. The systematic uncertainties are described in the section 7 of the paper. The total systematic uncertainties are determined by adding the contributions from all relevant sources in quadrature.

Nuclear modification factor in centrality interval 40-50% for Pb+Pb. The systematic uncertainties are described in the section 7 of the paper. The total systematic uncertainties are determined by adding the contributions from all relevant sources in quadrature.

Nuclear modification factor in centrality interval 40-50% for Pb+Pb. The systematic uncertainties are described in the section 7 of the paper. The total systematic uncertainties are determined by adding the contributions from all relevant sources in quadrature.

Nuclear modification factor in centrality interval 50-60% for Pb+Pb. The systematic uncertainties are described in the section 7 of the paper. The total systematic uncertainties are determined by adding the contributions from all relevant sources in quadrature.

Nuclear modification factor in centrality interval 50-60% for Pb+Pb. The systematic uncertainties are described in the section 7 of the paper. The total systematic uncertainties are determined by adding the contributions from all relevant sources in quadrature.

Nuclear modification factor in centrality interval 50-60% for Pb+Pb. The systematic uncertainties are described in the section 7 of the paper. The total systematic uncertainties are determined by adding the contributions from all relevant sources in quadrature.

Nuclear modification factor in centrality interval 50-60% for Pb+Pb. The systematic uncertainties are described in the section 7 of the paper. The total systematic uncertainties are determined by adding the contributions from all relevant sources in quadrature.

Nuclear modification factor in centrality interval 60-80% for Pb+Pb. The systematic uncertainties are described in the section 7 of the paper. The total systematic uncertainties are determined by adding the contributions from all relevant sources in quadrature.

Nuclear modification factor in centrality interval 60-80% for Pb+Pb. The systematic uncertainties are described in the section 7 of the paper. The total systematic uncertainties are determined by adding the contributions from all relevant sources in quadrature.

Nuclear modification factor in centrality interval 60-80% for Pb+Pb. The systematic uncertainties are described in the section 7 of the paper. The total systematic uncertainties are determined by adding the contributions from all relevant sources in quadrature.

Nuclear modification factor in centrality interval 60-80% for Pb+Pb. The systematic uncertainties are described in the section 7 of the paper. The total systematic uncertainties are determined by adding the contributions from all relevant sources in quadrature.

Nuclear modification factor in centrality interval 0-5% for Xe+Xe. The systematic uncertainties are described in the section 7 of the paper. The total systematic uncertainties are determined by adding the contributions from all relevant sources in quadrature.

Nuclear modification factor in centrality interval 0-5% for Xe+Xe. The systematic uncertainties are described in the section 7 of the paper. The total systematic uncertainties are determined by adding the contributions from all relevant sources in quadrature.

Nuclear modification factor in centrality interval 0-5% for Xe+Xe. The systematic uncertainties are described in the section 7 of the paper. The total systematic uncertainties are determined by adding the contributions from all relevant sources in quadrature.

Nuclear modification factor in centrality interval 5-10% for Xe+Xe. The systematic uncertainties are described in the section 7 of the paper. The total systematic uncertainties are determined by adding the contributions from all relevant sources in quadrature.

Nuclear modification factor in centrality interval 5-10% for Xe+Xe. The systematic uncertainties are described in the section 7 of the paper. The total systematic uncertainties are determined by adding the contributions from all relevant sources in quadrature.

Nuclear modification factor in centrality interval 5-10% for Xe+Xe. The systematic uncertainties are described in the section 7 of the paper. The total systematic uncertainties are determined by adding the contributions from all relevant sources in quadrature.

Nuclear modification factor in centrality interval 10-20% for Xe+Xe. The systematic uncertainties are described in the section 7 of the paper. The total systematic uncertainties are determined by adding the contributions from all relevant sources in quadrature.

Nuclear modification factor in centrality interval 10-20% for Xe+Xe. The systematic uncertainties are described in the section 7 of the paper. The total systematic uncertainties are determined by adding the contributions from all relevant sources in quadrature.

Nuclear modification factor in centrality interval 10-20% for Xe+Xe. The systematic uncertainties are described in the section 7 of the paper. The total systematic uncertainties are determined by adding the contributions from all relevant sources in quadrature.

Nuclear modification factor in centrality interval 20-30% for Xe+Xe. The systematic uncertainties are described in the section 7 of the paper. The total systematic uncertainties are determined by adding the contributions from all relevant sources in quadrature.

Nuclear modification factor in centrality interval 20-30% for Xe+Xe. The systematic uncertainties are described in the section 7 of the paper. The total systematic uncertainties are determined by adding the contributions from all relevant sources in quadrature.

Nuclear modification factor in centrality interval 20-30% for Xe+Xe. The systematic uncertainties are described in the section 7 of the paper. The total systematic uncertainties are determined by adding the contributions from all relevant sources in quadrature.

Nuclear modification factor in centrality interval 30-40% for Xe+Xe. The systematic uncertainties are described in the section 7 of the paper. The total systematic uncertainties are determined by adding the contributions from all relevant sources in quadrature.

Nuclear modification factor in centrality interval 30-40% for Xe+Xe. The systematic uncertainties are described in the section 7 of the paper. The total systematic uncertainties are determined by adding the contributions from all relevant sources in quadrature.

Nuclear modification factor in centrality interval 30-40% for Xe+Xe. The systematic uncertainties are described in the section 7 of the paper. The total systematic uncertainties are determined by adding the contributions from all relevant sources in quadrature.

Nuclear modification factor in centrality interval 40-50% for Xe+Xe. The systematic uncertainties are described in the section 7 of the paper. The total systematic uncertainties are determined by adding the contributions from all relevant sources in quadrature.

Nuclear modification factor in centrality interval 40-50% for Xe+Xe. The systematic uncertainties are described in the section 7 of the paper. The total systematic uncertainties are determined by adding the contributions from all relevant sources in quadrature.

Nuclear modification factor in centrality interval 40-50% for Xe+Xe. The systematic uncertainties are described in the section 7 of the paper. The total systematic uncertainties are determined by adding the contributions from all relevant sources in quadrature.

Nuclear modification factor in centrality interval 50-60% for Xe+Xe. The systematic uncertainties are described in the section 7 of the paper. The total systematic uncertainties are determined by adding the contributions from all relevant sources in quadrature.

Nuclear modification factor in centrality interval 50-60% for Xe+Xe. The systematic uncertainties are described in the section 7 of the paper. The total systematic uncertainties are determined by adding the contributions from all relevant sources in quadrature.

Nuclear modification factor in centrality interval 50-60% for Xe+Xe. The systematic uncertainties are described in the section 7 of the paper. The total systematic uncertainties are determined by adding the contributions from all relevant sources in quadrature.

Nuclear modification factor in centrality interval 60-80% for Xe+Xe. The systematic uncertainties are described in the section 7 of the paper. The total systematic uncertainties are determined by adding the contributions from all relevant sources in quadrature.

Nuclear modification factor in centrality interval 60-80% for Xe+Xe. The systematic uncertainties are described in the section 7 of the paper. The total systematic uncertainties are determined by adding the contributions from all relevant sources in quadrature.

Nuclear modification factor in centrality interval 60-80% for Xe+Xe. The systematic uncertainties are described in the section 7 of the paper. The total systematic uncertainties are determined by adding the contributions from all relevant sources in quadrature.

Measured unfolded differential cross-section as a function of the dilepton transverse momentum $p^{ll}_{\mathrm{T}}$. NLO predictions from Sherpa 2.2.10 and MadGraph5_aMC@NLO 2.7.3 are also shown. The uncertainty in the predictions is divided into statistical and theoretical uncertainties (scale and PDF+$\alpha_{s}$).

Measured unfolded differential cross-section as a function of the the four-body transverse momentum $p^{ll\gamma\gamma}_{\mathrm{T}}$. NLO predictions from Sherpa 2.2.10 and MadGraph5_aMC@NLO 2.7.3 are also shown. The uncertainty in the predictions is divided into statistical and theoretical uncertainties (scale and PDF+$\alpha_{s}$).

Measured unfolded differential cross-section as a function of the diphoton invariant mass $m_{\gamma\gamma}$. NLO predictions from Sherpa 2.2.10 and MadGraph5_aMC@NLO 2.7.3 are also shown. The uncertainty in the predictions is divided into statistical and theoretical uncertainties (scale and PDF+$\alpha_{s}$).

Measured unfolded differential cross-section as a function of the four-body invariant mass $m_{ll\gamma\gamma}$. NLO predictions from Sherpa 2.2.10 and MadGraph5_aMC@NLO 2.7.3 are also shown. The uncertainty in the predictions is divided into statistical and theoretical uncertainties (scale and PDF+$\alpha_{s}$).

Expected and observed $95\%$ confidence intervals for the coupling parameters $f_{T,j}/\Lambda^{4}$ of transverse dimension-8 operators. All parameter values outside of the stated range are excluded at the chosen confidence level. No unitarity constraints are applied.

Expected and observed unitarised $95\%$ confidence intervals for the coupling parameter $f_{T,8}/\Lambda^{4}$ in the clipping energy range between 1.1 and 5 TeV. The non-unitarised limits ($E_c = \infty$) are also shown. All parameter values outside of the stated range are excluded at the chosen confidence level.

Expected and observed unitarised $95\%$ confidence intervals for the coupling parameter $f_{T,0}/\Lambda^{4}$ in the clipping energy range between 1.1 and 5 TeV. The non-unitarised limits ($E_c = \infty$) are also shown. All parameter values outside of the stated range are excluded at the chosen confidence level.

Expected and observed unitarised $95\%$ confidence intervals for the coupling parameter $f_{T,1}/\Lambda^{4}$ in the clipping energy range between 1.1 and 5 TeV. The non-unitarised limits ($E_c = \infty$) are also shown. All parameter values outside of the stated range are excluded at the chosen confidence level.

Expected and observed unitarised $95\%$ confidence intervals for the coupling parameter $f_{T,2}/\Lambda^{4}$ in the clipping energy range between 1.1 and 5 TeV. The non-unitarised limits ($E_c = \infty$) are also shown. All parameter values outside of the stated range are excluded at the chosen confidence level.

Expected and observed unitarised $95\%$ confidence intervals for the coupling parameter $f_{T,5}/\Lambda^{4}$ in the clipping energy range between 1.1 and 5 TeV. The non-unitarised limits ($E_c = \infty$) are also shown. All parameter values outside of the stated range are excluded at the chosen confidence level.

Expected and observed unitarised $95\%$ confidence intervals for the coupling parameter $f_{T,6}/\Lambda^{4}$ in the clipping energy range between 1.1 and 5 TeV. The non-unitarised limits ($E_c = \infty$) are also shown. All parameter values outside of the stated range are excluded at the chosen confidence level.

Expected and observed unitarised $95\%$ confidence intervals for the coupling parameter $f_{T,7}/\Lambda^{4}$ in the clipping energy range between 1.1 and 5 TeV. The non-unitarised limits ($E_c = \infty$) are also shown. All parameter values outside of the stated range are excluded at the chosen confidence level.

Expected and observed unitarised $95\%$ confidence intervals for the coupling parameter $f_{T,9}/\Lambda^{4}$ in the clipping energy range between 1.1 and 5 TeV. The non-unitarised limits ($E_c = \infty$) are also shown. All parameter values outside of the stated range are excluded at the chosen confidence level.

Ratio of measured PSI(2S) over J/PSI(1S) cross section in charged-particle multiplicity intervals and integrated in multiplicity.

Ratio of measured PSI(2S) over J/PSI(1S) cross section in charged-particle multiplicity intervals and integrated in multiplicity.

Ratio of measured PSI(2S) over J/PSI(1S) cross section in charged-particle multiplicity intervals and integrated in multiplicity.

Ranking of the systematic uncertainties on the measured cross-section, normalised to the predicted value, in the inclusive fit to data. The impact of each nuisance parameter, $\Delta \sigma_{\text{inc}}/\sigma^{\text{pred.}}_{\text{inc}}$, is computed by comparing the nominal best-fit value of $\sigma_{\text{inc}}/\sigma^{\text{pred}}_{\text{inc}}$ with the result of the fit when fixing the considered nuisance parameter to its best-fit value, $\theta$, shifted by its pre-fit (post-fit) uncertainties $\pm \Delta \theta$ ($\pm \Delta \hat{\theta}$). The figure shows the effect of the ten most significant uncertainties.

Ranking of the systematic uncertainties on the measured cross-section, normalised to the predicted value, in the fiducial fit to data. The impact of each nuisance parameter, $\Delta \sigma_{\text{fid}}/\sigma^{\text{pred.}}_{\text{fid}}$, is computed by comparing the nominal best-fit value of $\sigma_{\text{fid}}/\sigma^{\text{pred}}_{\text{fid}}$ with the result of the fit when fixing the considered nuisance parameter to its best-fit value, $\theta$, shifted by its pre-fit (post-fit) uncertainties $\pm \Delta \theta$ ($\pm \Delta \hat{\theta}$). The figure shows the effect of the ten most significant uncertainties.

Ranking of the systematic uncertainties on the measured cross-section, normalised to the predicted value, in the fiducial fit to data. The impact of each nuisance parameter, $\Delta \sigma_{\text{fid}}/\sigma^{\text{pred.}}_{\text{fid}}$, is computed by comparing the nominal best-fit value of $\sigma_{\text{fid}}/\sigma^{\text{pred}}_{\text{fid}}$ with the result of the fit when fixing the considered nuisance parameter to its best-fit value, $\theta$, shifted by its pre-fit (post-fit) uncertainties $\pm \Delta \theta$ ($\pm \Delta \hat{\theta}$). The figure shows the effect of the ten most significant uncertainties.

Impact of different categories of systematic uncertainties on the fiducial and inclusive measurements. The quoted values are obtained by repeating the fit, fixing a set of nuisance parameters of the sources corresponding to the considered category, and subtracting in quadrature the resulting uncertainty from the total uncertainty of the nominal fit. The total uncertainty is different from the sum in quadrature of the different components due to correlations between nuisance parameters built by the fit.

Impact of different categories of systematic uncertainties on the fiducial and inclusive measurements. The quoted values are obtained by repeating the fit, fixing a set of nuisance parameters of the sources corresponding to the considered category, and subtracting in quadrature the resulting uncertainty from the total uncertainty of the nominal fit. The total uncertainty is different from the sum in quadrature of the different components due to correlations between nuisance parameters built by the fit.

Fiducial region definition

Fiducial region definition

The rho-parameter, i.e. the ratio of the real to imaginary part of the elastic scattering amplitude extrapolated to t=0. The systematic uncertainty includes experimental and theoretical uncerainties.

The nuclear slope parameter B from a fit of the form exp(-Bt-Ct^2-Dt^3). The systematic uncertainty includes experimental and theoretical uncerainties.

The nuclear slope parameter B from a fit of the form exp(-Bt-Ct^2-Dt^3). The systematic uncertainty includes experimental and theoretical uncerainties.

The nuclear slope parameter C from a fit of the form exp(-Bt-Ct^2-Dt^3). The systematic uncertainty includes experimental and theoretical uncerainties.

The nuclear slope parameter C from a fit of the form exp(-Bt-Ct^2-Dt^3). The systematic uncertainty includes experimental and theoretical uncerainties.

The nuclear slope parameter D from a fit of the form exp(-Bt-Ct^2-Dt^3). The systematic uncertainty includes experimental and theoretical uncerainties.

The nuclear slope parameter D from a fit of the form exp(-Bt-Ct^2-Dt^3). The systematic uncertainty includes experimental and theoretical uncerainties.

The total elastic cross section measured inside the fiducial volume. The systematic uncertainty includes experimental uncertainties.

The total elastic cross section measured inside the fiducial volume. The systematic uncertainty includes experimental uncertainties.

The total elastic cross section obtained from the fitted parameters, extrapolated to full phase space using only the nuclear amplitude.

The total elastic cross section obtained from the fitted parameters, extrapolated to full phase space using only the nuclear amplitude.

The total inelastic cross section.

The total inelastic cross section.

The ratio of elastic to total cross section.

The ratio of elastic to total cross section.

The differential elastic cross section as function of t with statistical and systematic uncertainties. The systematic uncertainties are given as signed relative change for 20 sources of experimental uncertainty associated to nuisance parameters used in the fit for the extraction of physics parameters.

The differential elastic cross section as function of t with statistical and systematic uncertainties. The systematic uncertainties are given as signed relative change for 20 sources of experimental uncertainty associated to nuisance parameters used in the fit for the extraction of physics parameters.

Statistical covariance matrix for the measurement of the differential elastic cross section as function of t.

Statistical covariance matrix for the measurement of the differential elastic cross section as function of t.

Jet multiplicity measured for a leading-pT jet ($p_{T1}$) with 400 < $p_{T1}$ < 800 GeV and for an azimuthal separation between the two leading jets of $0 < \Delta\Phi_{1,2} < 150^{\circ}$. The full breakdown of the uncertainties is displayed, with PU corresponding to Pileup, PREF to Trigger Prefering, PTHAT to the hard-scale (renormalization and factorization scales), MISS and FAKE to the inefficienties and background, LUMI to integrated luminosity. With JES, JER and stat. unc. following the notation in the paper.

Jet multiplicity measured for a leading-pT jet ($p_{T1}$) with 400 < $p_{T1}$ < 800 GeV and for an azimuthal separation between the two leading jets of $150 < \Delta\Phi_{1,2} < 170^{\circ}$. The full breakdown of the uncertainties is displayed, with PU corresponding to Pileup, PREF to Trigger Prefering, PTHAT to the hard-scale (renormalization and factorization scales), MISS and FAKE to the inefficienties and background, LUMI to integrated luminosity. With JES, JER and stat. unc. following the notation in the paper.

Jet multiplicity measured for a leading-pT jet ($p_{T1}$) with 400 < $p_{T1}$ < 800 GeV and for an azimuthal separation between the two leading jets of $170 < \Delta\Phi_{1,2} < 180^{\circ}$. The full breakdown of the uncertainties is displayed, with PU corresponding to Pileup, PREF to Trigger Prefering, PTHAT to the hard-scale (renormalization and factorization scales), MISS and FAKE to the inefficienties and background, LUMI to integrated luminosity. With JES, JER and stat. unc. following the notation in the paper.

Jet multiplicity measured for a leading-pT jet ($p_{T1}$) with $p_{T1}$ > 800 and for an azimuthal separation between the two leading jets of $0 < \Delta\Phi_{1,2} < 150^{\circ}$. The full breakdown of the uncertainties is displayed, with PU corresponding to Pileup, PREF to Trigger Prefering, PTHAT to the hard-scale (renormalization and factorization scales), MISS and FAKE to the inefficienties and background, LUMI to integrated luminosity. With JES, JER and stat. unc. following the notation in the paper.

Jet multiplicity measured for a leading-pT jet ($p_{T1}$) with $p_{T1}$ > 800 and for an azimuthal separation between the two leading jets of $150 < \Delta\Phi_{1,2} < 170^{\circ}$. The full breakdown of the uncertainties is displayed, with PU corresponding to Pileup, PREF to Trigger Prefering, PTHAT to the hard-scale (renormalization and factorization scales), MISS and FAKE to the inefficienties and background, LUMI to integrated luminosity. With JES, JER and stat. unc. following the notation in the paper.

Jet multiplicity measured for a leading-pT jet ($p_{T1}$) with $p_{T1}$ > 800 and for an azimuthal separation between the two leading jets of $170 < \Delta\Phi_{1,2} < 180^{\circ}$. The full breakdown of the uncertainties is displayed, with PU corresponding to Pileup, PREF to Trigger Prefering, PTHAT to the hard-scale (renormalization and factorization scales), MISS and FAKE to the inefficienties and background, LUMI to integrated luminosity. With JES, JER and stat. unc. following the notation in the paper.

Measured transverse momentum of the leading $p_{T}$ jet ($p_{T1}$). The full breakdown of the uncertainties is displayed, with PU corresponding to Pileup, PREF to Trigger Prefering, PTHAT to the hard-scale (renormalization and factorization scales), MISS and FAKE to the inefficienties and background, LUMI to integrated luminosity. With JES, JER and stat. unc. following the notation in the paper.

Measured transverse momentum of the subleading $p_{T}$ jet ($p_{T2}$). The full breakdown of the uncertainties is displayed, with PU corresponding to Pileup, PREF to Trigger Prefering, PTHAT to the hard-scale (renormalization and factorization scales), MISS and FAKE to the inefficienties and background, LUMI to integrated luminosity. With JES, JER and stat. unc. following the notation in the paper.

Measured transverse momentum of the third leading $p_{T}$ jet ($p_{T3}$). The full breakdown of the uncertainties is displayed, with PU corresponding to Pileup, PREF to Trigger Prefering, PTHAT to the hard-scale (renormalization and factorization scales), MISS and FAKE to the inefficienties and background, LUMI to integrated luminosity. With JES, JER and stat. unc. following the notation in the paper.

Measured transverse momentum of the fourth leading $p_{T}$ jet ($p_{T4}$). The full breakdown of the uncertainties is displayed, with PU corresponding to Pileup, PREF to Trigger Prefering, PTHAT to the hard-scale (renormalization and factorization scales), MISS and FAKE to the inefficienties and background, LUMI to integrated luminosity. With JES, JER and stat. unc. following the notation in the paper.

Jet multiplicity distribution in TUnfold binning ($N_{jets},\Delta\phi_{1,2},p_{T1}$) as indicated in the XML file provided as additional resource. The uncertainties follow the notation of Table 1.

Correlation matrix at particle level for the measured jet multiplicity in TUnfold binning ($N_{jets},\Delta\phi_{1,2},p_{T1}$) as indicated in the XML file provided as additional resource.

Jet $p_{T}$ distributions in TUnfold binning ($p_{T1},p_{T2},p_{T3},p_{T4}$) as indicated in the XML file provided as additional resource. The uncertainties follow the notation of Table 1.

Correlation matrix at particle level for the measured jet $p_{T}$ distributions in TUnfold binning ($p_{T1},p_{T2},p_{T3},p_{T4}$) as indicated in the XML file provided as additional resource.

The extrapolated ($\tau > 30$ ps) average charged-particle muliplicity per unit of rapidity for ETARAP=0 as a function of the centre-of-mass energy.

Charged-particle multiplicities in proton-proton collisions at a centre-of mass energy of 13000 GeV as a function of pseudorapidity for events with the number of charged particles >=2 having transverse momentum >100 MeV and absolute(pseudorapidity) <2.5.

Charged-particle multiplicities in proton-proton collisions at a centre-of mass energy of 13000 GeV as a function of pseudorapidity for events with the number of charged particles >=2 having transverse momentum >100 MeV and absolute(pseudorapidity) <2.5.

Charged-particle multiplicities in proton-proton collisions at a centre-of mass energy of 13000 GeV as a function of transverse momentum for events with the number of charged particles >=2 having transverse momentum >100 MeV and absolute(pseudorapidity) <2.5.

Charged-particle multiplicities in proton-proton collisions at a centre-of mass energy of 13000 GeV as a function of transverse momentum for events with the number of charged particles >=2 having transverse momentum >100 MeV and absolute(pseudorapidity) <2.5.

Charged-particle multiplicity distribution in proton-proton collisions at a centre-of mass energy of 13000 GeV for events with the number of charged particles >=2 having transverse momentum >100 MeV and absolute(pseudorapidity) <2.5.

Charged-particle multiplicity distribution in proton-proton collisions at a centre-of mass energy of 13000 GeV for events with the number of charged particles >=2 having transverse momentum >100 MeV and absolute(pseudorapidity) <2.5.

Average transverse momentum in proton-proton collisions at a centre-of mass energy of 13000 GeV as a function of the number of charged particles in the event for events with the number of charged particles >=2 having transverse momentum >100 MeV and absolute(pseudorapidity) <2.5.

Average transverse momentum in proton-proton collisions at a centre-of mass energy of 13000 GeV as a function of the number of charged particles in the event for events with the number of charged particles >=2 having transverse momentum >100 MeV and absolute(pseudorapidity) <2.5.

Extrapolated ($\tau > 30$ ps) charged-particle multiplicities in proton-proton collisions at a centre-of mass energy of 13000 GeV as a function of pseudorapidity for events with the number of charged particles >=2 having transverse momentum >100 MeV and absolute(pseudorapidity) <2.5.

Extrapolated ($\tau > 30$ ps) charged-particle multiplicities in proton-proton collisions at a centre-of mass energy of 13000 GeV as a function of pseudorapidity for events with the number of charged particles >=2 having transverse momentum >100 MeV and absolute(pseudorapidity) <2.5.

Extrapolated ($\tau > 30$ ps) charged-particle multiplicities in proton-proton collisions at a centre-of mass energy of 13000 GeV as a function of transverse momentum for events with the number of charged particles >=2 having transverse momentum >100 MeV and absolute(pseudorapidity) <2.5.

Extrapolated ($\tau > 30$ ps) charged-particle multiplicities in proton-proton collisions at a centre-of mass energy of 13000 GeV as a function of transverse momentum for events with the number of charged particles >=2 having transverse momentum >100 MeV and absolute(pseudorapidity) <2.5.

Extrapolated ($\tau > 30$ ps) charged-particle multiplicity distributions in proton-proton collisions at a centre-of mass energy of 13000 GeV for events with the number of charged particles >=2 having transverse momentum >100 MeV and absolute(pseudorapidity) <2.5.

Extrapolated ($\tau > 30$ ps) charged-particle multiplicity distributions in proton-proton collisions at a centre-of mass energy of 13000 GeV for events with the number of charged particles >=2 having transverse momentum >100 MeV and absolute(pseudorapidity) <2.5.

Extrapolated ($\tau > 30$ ps) average transverse momentum in proton-proton collisions at a centre-of mass energy of 13000 GeV as a function of the number of charged particles in the event for events with the number of charged particles >=2 having transverse momentum >100 MeV and absolute(pseudorapidity) <2.5.

Extrapolated ($\tau > 30$ ps) average transverse momentum in proton-proton collisions at a centre-of mass energy of 13000 GeV as a function of the number of charged particles in the event for events with the number of charged particles >=2 having transverse momentum >100 MeV and absolute(pseudorapidity) <2.5.

The Bose-Einstein correlation (BEC) parameter R as a function of n_ch for HMT events using different MC generators in the calculation of R₂(Q). The uncertainties shown are statistical. The lower panel of each plot shows the ratio of the BEC parameters obtained using EPOS LHC (red circles), Pythia 8 Monash (blue squares) and Herwig++ UE-EE-5 (green triangles) compared with the parameters obtained using Pythia 8 A2. The gray band in the lower panels is the MC systematic uncertainty, obtained as explained in the text.

The Bose-Einstein correlation (BEC) parameter R as a function of k_T for MB events using different MC generators in the calculation of R₂(Q). The uncertainties shown are statistical. The lower panel of each plot shows the ratio of the BEC parameters obtained using EPOS LHC (red circles), Pythia 8 Monash (blue squares) and Herwig++ UE-EE-5 (green triangles) compared with the parameters obtained using Pythia 8 A2. The gray band in the lower panels is the MC systematic uncertainty, obtained as explained in the text.

The Bose-Einstein correlation (BEC) parameter &lambda; as a function of k_T for MB events using different MC generators in the calculation of R₂(Q). The uncertainties shown are statistical. The lower panel of each plot shows the ratio of the BEC parameters obtained using EPOS LHC (red circles), Pythia 8 Monash (blue squares) and Herwig++ UE-EE-5 (green triangles) compared with the parameters obtained using Pythia 8 A2. The gray band in the lower panels is the MC systematic uncertainty, obtained as explained in the text.

The two-particle double-ratio correlation function, R₂(Q), for pp collisions for track p_T &gt;100&nbsp;MeV at &radic;s=13&nbsp;TeV in the multiplicity interval 71 &le; n_ch &lt; 80 for the minimum-bias (MB) events. The blue dashed and red solid lines show the results of the exponential and Gaussian fits, respectively. The region excluded from the fits is shown. The statistical uncertainty and the systematic uncertainty for imperfections in the data reconstruction procedure are added in quadrature.

The two-particle double-ratio correlation function, R₂(Q), for pp collisions for track p_T &gt;100&nbsp;MeV at &radic;s=13&nbsp;TeV in the multiplicity interval 231 &le; n_ch &lt; 300 for the high-multiplicity track (HMT) events. The blue dashed and red solid lines show the results of the exponential and Gaussian fits, respectively. The region excluded from the fits is shown. The statistical uncertainty and the systematic uncertainty for imperfections in the data reconstruction procedure are added in quadrature.

The dependence of the correlation strength, &lambda;(m_ch), on rescaled multiplicity, m_ch, obtained from the exponential fit of the R₂(Q) correlation functions for tracks with p_T &gt; 100 &nbsp;MeV and p_T &gt; 500&nbsp;MeV at &radic;s = 13 &nbsp;TeV for the minimum-bias (MB) and high multiplicity track (HMT) data. The uncertainties represent the sum in quadrature of the statistical and asymmetric systematic contributions. The black and blue solid curves represent the exponential fit of &lambda;(m_ch) for p_T &gt;100&nbsp;MeV and p_T &gt;500&nbsp;MeV, respectively.

The dependence of the correlation strength, &lambda;(m_ch), on rescaled multiplicity, m_ch, obtained from the exponential fit of the R₂(Q) correlation functions for tracks with p_T &gt; 100 &nbsp;MeV and p_T &gt; 500&nbsp;MeV at &radic;s = 13 &nbsp;TeV for the minimum-bias (MB) and high multiplicity track (HMT) data. The uncertainties represent the sum in quadrature of the statistical and asymmetric systematic contributions. The black and blue solid curves represent the exponential fit of &lambda;(m_ch) for p_T &gt;100&nbsp;MeV and p_T &gt;500&nbsp;MeV, respectively.

The dependence of the correlation strength, &lambda;(m_ch), on rescaled multiplicity, m_ch, obtained from the exponential fit of the R₂(Q) correlation functions for tracks with p_T &gt; 100 &nbsp;MeV and p_T &gt; 500&nbsp;MeV at &radic;s = 13 &nbsp;TeV for the minimum-bias (MB) and high multiplicity track (HMT) data. The uncertainties represent the sum in quadrature of the statistical and asymmetric systematic contributions. The black and blue solid curves represent the exponential fit of &lambda;(m_ch) for p_T &gt;100&nbsp;MeV and p_T &gt;500&nbsp;MeV, respectively.

The dependence of the correlation strength, &lambda;(m_ch), on rescaled multiplicity, m_ch, obtained from the exponential fit of the R₂(Q) correlation functions for tracks with p_T &gt; 100 &nbsp;MeV and p_T &gt; 500&nbsp;MeV at &radic;s = 13 &nbsp;TeV for the minimum-bias (MB) and high multiplicity track (HMT) data. The uncertainties represent the sum in quadrature of the statistical and asymmetric systematic contributions. The black and blue solid curves represent the exponential fit of &lambda;(m_ch) for p_T &gt;100&nbsp;MeV and p_T &gt;500&nbsp;MeV, respectively.

The dependence of the source radius, R(m_ch), on m_ch. The uncertainties represent the sum in quadrature of the statistical and asymmetric systematic contributions. The black and blue solid curves represent the fit of R(m_ch) for &#8731;m_ch &lt; 1.2 for p_T &gt;100&nbsp;MeV and p_T &gt;500&nbsp;MeV, respectively. The black and blue dotted curves are extensions of the black and blue solid curves beyond &#8731;m_ch &gt; 1.2, respectively. The black and brown dashed curves represent the saturation value of R(m_ch) for &#8731;m_ch &gt; 1.45 with p_T &gt;100&nbsp;MeV and for &#8731;m_ch &gt; 1.6 with p_T &gt;500&nbsp;MeV, respectively.

The dependence of the source radius, R(m_ch), on m_ch. The uncertainties represent the sum in quadrature of the statistical and asymmetric systematic contributions. The black and blue solid curves represent the fit of R(m_ch) for &#8731;m_ch &lt; 1.2 for p_T &gt;100&nbsp;MeV and p_T &gt;500&nbsp;MeV, respectively. The black and blue dotted curves are extensions of the black and blue solid curves beyond &#8731;m_ch &gt; 1.2, respectively. The black and brown dashed curves represent the saturation value of R(m_ch) for &#8731;m_ch &gt; 1.45 with p_T &gt;100&nbsp;MeV and for &#8731;m_ch &gt; 1.6 with p_T &gt;500&nbsp;MeV, respectively.

The dependence of the source radius, R(m_ch), on m_ch. The uncertainties represent the sum in quadrature of the statistical and asymmetric systematic contributions. The black and blue solid curves represent the fit of R(m_ch) for &#8731;m_ch &lt; 1.2 for p_T &gt;100&nbsp;MeV and p_T &gt;500&nbsp;MeV, respectively. The black and blue dotted curves are extensions of the black and blue solid curves beyond &#8731;m_ch &gt; 1.2, respectively. The black and brown dashed curves represent the saturation value of R(m_ch) for &#8731;m_ch &gt; 1.45 with p_T &gt;100&nbsp;MeV and for &#8731;m_ch &gt; 1.6 with p_T &gt;500&nbsp;MeV, respectively.

The dependence of the source radius, R(m_ch), on m_ch. The uncertainties represent the sum in quadrature of the statistical and asymmetric systematic contributions. The black and blue solid curves represent the fit of R(m_ch) for &#8731;m_ch &lt; 1.2 for p_T &gt;100&nbsp;MeV and p_T &gt;500&nbsp;MeV, respectively. The black and blue dotted curves are extensions of the black and blue solid curves beyond &#8731;m_ch &gt; 1.2, respectively. The black and brown dashed curves represent the saturation value of R(m_ch) for &#8731;m_ch &gt; 1.45 with p_T &gt;100&nbsp;MeV and for &#8731;m_ch &gt; 1.6 with p_T &gt;500&nbsp;MeV, respectively.

The dependence of the R(m_ch) on &#8731;m_ch. The uncertainties represent the sum in quadrature of the statistical and asymmetric systematic contributions. The black and blue solid curves represent the fit of R(m_ch) for &#8731;m_ch &lt; 1.2 for p_T &gt;100&nbsp;MeV and p_T &gt;500&nbsp;MeV, respectively. The black and blue dotted curves are extensions of the black and blue solid curves beyond &#8731;m_ch &gt; 1.2, respectively. The black and brown dashed curves represent the saturation value of R(m_ch) for &#8731;m_ch &gt; 1.45 with p_T &gt;100&nbsp;MeV and for &#8731;m_ch &gt; 1.6 with p_T &gt;500&nbsp;MeV, respectively

The dependence of the R(m_ch) on &#8731;m_ch. The uncertainties represent the sum in quadrature of the statistical and asymmetric systematic contributions. The black and blue solid curves represent the fit of R(m_ch) for &#8731;m_ch &lt; 1.2 for p_T &gt;100&nbsp;MeV and p_T &gt;500&nbsp;MeV, respectively. The black and blue dotted curves are extensions of the black and blue solid curves beyond &#8731;m_ch &gt; 1.2, respectively. The black and brown dashed curves represent the saturation value of R(m_ch) for &#8731;m_ch &gt; 1.45 with p_T &gt;100&nbsp;MeV and for &#8731;m_ch &gt; 1.6 with p_T &gt;500&nbsp;MeV, respectively

The dependence of the R(m_ch) on &#8731;m_ch. The uncertainties represent the sum in quadrature of the statistical and asymmetric systematic contributions. The black and blue solid curves represent the fit of R(m_ch) for &#8731;m_ch &lt; 1.2 for p_T &gt;100&nbsp;MeV and p_T &gt;500&nbsp;MeV, respectively. The black and blue dotted curves are extensions of the black and blue solid curves beyond &#8731;m_ch &gt; 1.2, respectively. The black and brown dashed curves represent the saturation value of R(m_ch) for &#8731;m_ch &gt; 1.45 with p_T &gt;100&nbsp;MeV and for &#8731;m_ch &gt; 1.6 with p_T &gt;500&nbsp;MeV, respectively

The dependence of the R(m_ch) on &#8731;m_ch. The uncertainties represent the sum in quadrature of the statistical and asymmetric systematic contributions. The black and blue solid curves represent the fit of R(m_ch) for &#8731;m_ch &lt; 1.2 for p_T &gt;100&nbsp;MeV and p_T &gt;500&nbsp;MeV, respectively. The black and blue dotted curves are extensions of the black and blue solid curves beyond &#8731;m_ch &gt; 1.2, respectively. The black and brown dashed curves represent the saturation value of R(m_ch) for &#8731;m_ch &gt; 1.45 with p_T &gt;100&nbsp;MeV and for &#8731;m_ch &gt; 1.6 with p_T &gt;500&nbsp;MeV, respectively

Comparison of single-ratio two-particle correlation functions, using the unlike-charge particle (UCP) pair reference sample, for minimum-bias (MB) events, showing C₂^data(Q) (top panel) at 13&nbsp;TeV (black circles) and 7&nbsp;TeV (open blue circles), and the ratio of C₂^7&nbsp;TeV (Q) to C₂^13&nbsp;TeV (Q) (bottom panel). Comparison of C₂^data (Q) for representative multiplicity region 3.09 &lt; m_ch &le; 3.86. The statistical and systematic uncertainties, combined in quadrature, are presented. The systematic uncertainties include track efficiency, Coulomb correction, non-closure and multiplicity-unfolding uncertainties.

Comparison of single-ratio two-particle correlation functions, using the unlike-charge particle (UCP) pair reference sample, for minimum-bias (MB) events, showing C₂^data(Q) (top panel) at 13 &nbsp;TeV (black circles) and 7 &nbsp;TeV (open blue circles), and the ratio of C₂^7&nbsp;TeV (Q) to C₂^13&nbsp;TeV (Q) (bottom panel). Comparison of C₂^data (Q) for representative k_T region 400 &lt; k_T &le;500 &nbsp;MeV. The statistical and systematic uncertainties, combined in quadrature, are presented. The systematic uncertainties include track efficiency, Coulomb correction, non-closure and multiplicity-unfolding uncertainties.

The k_T dependence of the correlation strength, &lambda;(k_T), obtained from the exponential fit to the R₂(Q) correlation functions for events with multiplicity n_ch &ge; 2 and transfer momentum of tracks with p_T &gt;100&nbsp;MeV and p_T &gt;500&nbsp;MeV at &radic;s=13&nbsp;TeV for the minimum-bias (MB) and high-multiplicity track (HMT) events. The uncertainties represent the sum in quadrature of the statistical and systematic contributions. The curves represent the exponential fits to &lambda;(k_T).

The k_T dependence of the correlation strength, &lambda;(k_T), obtained from the exponential fit to the R₂(Q) correlation functions for events with multiplicity n_ch &ge; 2 and transfer momentum of tracks with p_T &gt;100&nbsp;MeV and p_T &gt;500&nbsp;MeV at &radic;s=13&nbsp;TeV for the minimum-bias (MB) and high-multiplicity track (HMT) events. The uncertainties represent the sum in quadrature of the statistical and systematic contributions. The curves represent the exponential fits to &lambda;(k_T).

The k_T dependence of the correlation strength, &lambda;(k_T), obtained from the exponential fit to the R₂(Q) correlation functions for events with multiplicity n_ch &ge; 2 and transfer momentum of tracks with p_T &gt;100&nbsp;MeV and p_T &gt;500&nbsp;MeV at &radic;s=13&nbsp;TeV for the minimum-bias (MB) and high-multiplicity track (HMT) events. The uncertainties represent the sum in quadrature of the statistical and systematic contributions. The curves represent the exponential fits to &lambda;(k_T).

The k_T dependence of the correlation strength, &lambda;(k_T), obtained from the exponential fit to the R₂(Q) correlation functions for events with multiplicity n_ch &ge; 2 and transfer momentum of tracks with p_T &gt;100&nbsp;MeV and p_T &gt;500&nbsp;MeV at &radic;s=13&nbsp;TeV for the minimum-bias (MB) and high-multiplicity track (HMT) events. The uncertainties represent the sum in quadrature of the statistical and systematic contributions. The curves represent the exponential fits to &lambda;(k_T).

The k_T dependence of the source radius, R(k_T), obtained from the exponential fit to the R₂(Q) correlation functions for events with multiplicity n_ch &ge; 2 and transfer momentum of tracks with p_T &gt;100&nbsp;MeV and p_T &gt;500&nbsp;MeV at &radic;s=13&nbsp;TeV for the minimum-bias (MB) and high-multiplicity track (HMT) events. The uncertainties represent the sum in quadrature of the statistical and systematic contributions. The curves represent the exponential fits to R(k_T).

The k_T dependence of the source radius, R(k_T), obtained from the exponential fit to the R₂(Q) correlation functions for events with multiplicity n_ch &ge; 2 and transfer momentum of tracks with p_T &gt;100&nbsp;MeV and p_T &gt;500&nbsp;MeV at &radic;s=13&nbsp;TeV for the minimum-bias (MB) and high-multiplicity track (HMT) events. The uncertainties represent the sum in quadrature of the statistical and systematic contributions. The curves represent the exponential fits to R(k_T).

The k_T dependence of the source radius, R(k_T), obtained from the exponential fit to the R₂(Q) correlation functions for events with multiplicity n_ch &ge; 2 and transfer momentum of tracks with p_T &gt;100&nbsp;MeV and p_T &gt;500&nbsp;MeV at &radic;s=13&nbsp;TeV for the minimum-bias (MB) and high-multiplicity track (HMT) events. The uncertainties represent the sum in quadrature of the statistical and systematic contributions. The curves represent the exponential fits to R(k_T).

The k_T dependence of the source radius, R(k_T), obtained from the exponential fit to the R₂(Q) correlation functions for events with multiplicity n_ch &ge; 2 and transfer momentum of tracks with p_T &gt;100&nbsp;MeV and p_T &gt;500&nbsp;MeV at &radic;s=13&nbsp;TeV for the minimum-bias (MB) and high-multiplicity track (HMT) events. The uncertainties represent the sum in quadrature of the statistical and systematic contributions. The curves represent the exponential fits to R(k_T).

The two-dimensional dependence on m_ch and k_T for p_T &gt; 100&nbsp;MeV for the correlation strength, &lambda;, obtained from the exponential fit to the R₂(Q) correlation functions using the MB sample for m_ch &le; 3.08 and the HMT sample for m_ch &gt; 3.08.

The two-dimensional dependence on m_ch and k_T for p_T &gt; 100&nbsp;MeV for the source radius, R, obtained from the exponential fit to the R₂(Q) correlation functions using the MB sample for m_ch &le; 3.08 and the HMT sample for m_ch &gt; 3.08.

The parameter &lambda; for p_T &gt; 100&nbsp;MeV as a function of k_T in selected low m_ch intervals. The error bars and boxes represent the statistical and systematic contributions, respectively.

The parameter &lambda; for p_T &gt; 100&nbsp;MeV as a function of k_T in selected low m_ch intervals. The error bars and boxes represent the statistical and systematic contributions, respectively.

The parameter &lambda; for p_T &gt; 100&nbsp;MeV as a function of k_T in selected high m_ch intervals. The error bars and boxes represent the statistical and systematic contributions, respectively.

The parameter &lambda; for p_T &gt; 100&nbsp;MeV as a function of k_T in selected high m_ch intervals. The error bars and boxes represent the statistical and systematic contributions, respectively.

The parameter &lambda; for p_T &gt; 100&nbsp;MeV as a function of m_ch in k_T intervals between 0.1 and 0.5&nbsp;GeV. The error bars and boxes represent the statistical and systematic contributions, respectively.

The parameter &lambda; for p_T &gt; 100&nbsp;MeV as a function of m_ch in k_T intervals between 0.1 and 0.5&nbsp;GeV. The error bars and boxes represent the statistical and systematic contributions, respectively.

The parameter &lambda; for p_T &gt; 100&nbsp;MeV as a function of m_ch in k_T intervals between 0.5 and 1.5&nbsp;GeV. The error bars and boxes represent the statistical and systematic contributions, respectively.

The parameter &lambda; for p_T &gt; 100&nbsp;MeV as a function of m_ch in k_T intervals between 0.5 and 1.5&nbsp;GeV. The error bars and boxes represent the statistical and systematic contributions, respectively.

The parameter R for p_T &gt; 100&nbsp;MeV as a function of k_T in selected low m_ch intervals. The error bars and boxes represent the statistical and systematic contributions, respectively.

The parameter R for p_T &gt; 100&nbsp;MeV as a function of k_T in selected low m_ch intervals. The error bars and boxes represent the statistical and systematic contributions, respectively.

The parameter R for p_T &gt; 100&nbsp;MeV as a function of k_T in selected high m_ch intervals. The error bars and boxes represent the statistical and systematic contributions, respectively.

The parameter R for p_T &gt; 100&nbsp;MeV as a function of k_T in selected high m_ch intervals. The error bars and boxes represent the statistical and systematic contributions, respectively.

The parameter R for p_T &gt; 100&nbsp;MeV as a function of m_ch in k_T intervals between 0.1 and 0.5&nbsp;GeV. The error bars and boxes represent the statistical and systematic contributions, respectively.

The parameter R for p_T &gt; 100&nbsp;MeV as a function of m_ch in k_T intervals between 0.1 and 0.5&nbsp;GeV. The error bars and boxes represent the statistical and systematic contributions, respectively.

The parameter R for p_T &gt; 100&nbsp;MeV as a function of m_ch in k_T intervals between 0.5 and 1.5&nbsp;GeV. The error bars and boxes represent the statistical and systematic contributions, respectively.

The parameter R for p_T &gt; 100&nbsp;MeV as a function of m_ch in k_T intervals between 0.5 and 1.5&nbsp;GeV. The error bars and boxes represent the statistical and systematic contributions, respectively.

The fit parameter &mu; describing the dependence of the correlation strength, &lambda;, on charged-particle scaled multiplicity, for track p_T&gt;100&nbsp;MeV and track p_T&gt;500&nbsp;MeV in the minimum-bias (MB) and high-multiplicity track (HMT) samples at &radic;s = 13&nbsp;TeV. The error bars and boxes represent the statistical and systematic contributions, respectively. The black solid (blue dashed) curve represents the exponential fit of the dependence of parameter &mu; on m_ch for tracks with p_T &gt;100&nbsp;MeV (p_T &gt;500&nbsp;MeV).

The fit parameter &mu; describing the dependence of the correlation strength, &lambda;, on charged-particle scaled multiplicity, for track p_T&gt;100&nbsp;MeV and track p_T&gt;500&nbsp;MeV in the minimum-bias (MB) and high-multiplicity track (HMT) samples at &radic;s = 13&nbsp;TeV. The error bars and boxes represent the statistical and systematic contributions, respectively. The black solid (blue dashed) curve represents the exponential fit of the dependence of parameter &mu; on m_ch for tracks with p_T &gt;100&nbsp;MeV (p_T &gt;500&nbsp;MeV).

The fit parameter &mu; describing the dependence of the correlation strength, &lambda;, on charged-particle scaled multiplicity, for track p_T&gt;100&nbsp;MeV and track p_T&gt;500&nbsp;MeV in the minimum-bias (MB) and high-multiplicity track (HMT) samples at &radic;s = 13&nbsp;TeV. The error bars and boxes represent the statistical and systematic contributions, respectively. The black solid (blue dashed) curve represents the exponential fit of the dependence of parameter &mu; on m_ch for tracks with p_T &gt;100&nbsp;MeV (p_T &gt;500&nbsp;MeV).

The fit parameter &mu; describing the dependence of the correlation strength, &lambda;, on charged-particle scaled multiplicity, for track p_T&gt;100&nbsp;MeV and track p_T&gt;500&nbsp;MeV in the minimum-bias (MB) and high-multiplicity track (HMT) samples at &radic;s = 13&nbsp;TeV. The error bars and boxes represent the statistical and systematic contributions, respectively. The black solid (blue dashed) curve represents the exponential fit of the dependence of parameter &mu; on m_ch for tracks with p_T &gt;100&nbsp;MeV (p_T &gt;500&nbsp;MeV).

The fit parameter &nu; describing the dependence of the correlation strength, &lambda;, on charged-particle scaled multiplicity, for track p_T&gt;100&nbsp;MeV and track p_T&gt;500&nbsp;MeV in the minimum-bias (MB) and high-multiplicity track (HMT) samples at &radic;s = 13&nbsp;TeV. The error bars and boxes represent the statistical and systematic contributions, respectively. The black solid (blue dashed) curve represents the exponential fit of the dependence of parameter &nu; on m_ch for tracks with p_T &gt;100&nbsp;MeV (p_T &gt;500&nbsp;MeV).

The fit parameter &nu; describing the dependence of the correlation strength, &lambda;, on charged-particle scaled multiplicity, for track p_T&gt;100&nbsp;MeV and track p_T&gt;500&nbsp;MeV in the minimum-bias (MB) and high-multiplicity track (HMT) samples at &radic;s = 13&nbsp;TeV. The error bars and boxes represent the statistical and systematic contributions, respectively. The black solid (blue dashed) curve represents the exponential fit of the dependence of parameter &nu; on m_ch for tracks with p_T &gt;100&nbsp;MeV (p_T &gt;500&nbsp;MeV).

The fit parameter &nu; describing the dependence of the correlation strength, &lambda;, on charged-particle scaled multiplicity, for track p_T&gt;100&nbsp;MeV and track p_T&gt;500&nbsp;MeV in the minimum-bias (MB) and high-multiplicity track (HMT) samples at &radic;s = 13&nbsp;TeV. The error bars and boxes represent the statistical and systematic contributions, respectively. The black solid (blue dashed) curve represents the exponential fit of the dependence of parameter &nu; on m_ch for tracks with p_T &gt;100&nbsp;MeV (p_T &gt;500&nbsp;MeV).

The fit parameter &nu; describing the dependence of the correlation strength, &lambda;, on charged-particle scaled multiplicity, for track p_T&gt;100&nbsp;MeV and track p_T&gt;500&nbsp;MeV in the minimum-bias (MB) and high-multiplicity track (HMT) samples at &radic;s = 13&nbsp;TeV. The error bars and boxes represent the statistical and systematic contributions, respectively. The black solid (blue dashed) curve represents the exponential fit of the dependence of parameter &nu; on m_ch for tracks with p_T &gt;100&nbsp;MeV (p_T &gt;500&nbsp;MeV).

The parameter &xi; describing the dependence of the source radius, R, on charged-particle scaled multiplicity, m_ch, for track p_T&gt;100&nbsp;MeV and track p_T&gt;500&nbsp;MeV in the minimum-bias (MB) and high-multiplicity track (HMT) samples at &radic;s = 13&nbsp;TeV. The error bars and boxes represent the statistical and systematic contributions, respectively. The black solid and blue dashed curves represent the saturated value of the parameter &xi; for m_ch &gt; 3.0 for tracks with p_T &gt;100&nbsp;MeV and for m_ch &gt; 2.8 for tracks with p_T &gt;500&nbsp;MeV, respectively.

The parameter &xi; describing the dependence of the source radius, R, on charged-particle scaled multiplicity, m_ch, for track p_T&gt;100&nbsp;MeV and track p_T&gt;500&nbsp;MeV in the minimum-bias (MB) and high-multiplicity track (HMT) samples at &radic;s = 13&nbsp;TeV. The error bars and boxes represent the statistical and systematic contributions, respectively. The black solid and blue dashed curves represent the saturated value of the parameter &xi; for m_ch &gt; 3.0 for tracks with p_T &gt;100&nbsp;MeV and for m_ch &gt; 2.8 for tracks with p_T &gt;500&nbsp;MeV, respectively.

The parameter &xi; describing the dependence of the source radius, R, on charged-particle scaled multiplicity, m_ch, for track p_T&gt;100&nbsp;MeV and track p_T&gt;500&nbsp;MeV in the minimum-bias (MB) and high-multiplicity track (HMT) samples at &radic;s = 13&nbsp;TeV. The error bars and boxes represent the statistical and systematic contributions, respectively. The black solid and blue dashed curves represent the saturated value of the parameter &xi; for m_ch &gt; 3.0 for tracks with p_T &gt;100&nbsp;MeV and for m_ch &gt; 2.8 for tracks with p_T &gt;500&nbsp;MeV, respectively.

The parameter &xi; describing the dependence of the source radius, R, on charged-particle scaled multiplicity, m_ch, for track p_T&gt;100&nbsp;MeV and track p_T&gt;500&nbsp;MeV in the minimum-bias (MB) and high-multiplicity track (HMT) samples at &radic;s = 13&nbsp;TeV. The error bars and boxes represent the statistical and systematic contributions, respectively. The black solid and blue dashed curves represent the saturated value of the parameter &xi; for m_ch &gt; 3.0 for tracks with p_T &gt;100&nbsp;MeV and for m_ch &gt; 2.8 for tracks with p_T &gt;500&nbsp;MeV, respectively.

The parameter &kappa; describing the dependence of the source radius, R, on charged-particle scaled multiplicity, m_ch, for track p_T&gt;100&nbsp;MeV and track p_T&gt;500&nbsp;MeV in the minimum-bias (MB) and high-multiplicity track (HMT) samples at &radic;s = 13&nbsp;TeV. The error bars and boxes represent the statistical and systematic contributions, respectively. The black solid and blue dashed curves represent the exponential fit to the parameter &kappa; for tracks with p_T &gt;100&nbsp;MeV and for tracks with p_T &gt;500&nbsp;MeV, respectively.

The parameter &kappa; describing the dependence of the source radius, R, on charged-particle scaled multiplicity, m_ch, for track p_T&gt;100&nbsp;MeV and track p_T&gt;500&nbsp;MeV in the minimum-bias (MB) and high-multiplicity track (HMT) samples at &radic;s = 13&nbsp;TeV. The error bars and boxes represent the statistical and systematic contributions, respectively. The black solid and blue dashed curves represent the exponential fit to the parameter &kappa; for tracks with p_T &gt;100&nbsp;MeV and for tracks with p_T &gt;500&nbsp;MeV, respectively.

The parameter &kappa; describing the dependence of the source radius, R, on charged-particle scaled multiplicity, m_ch, for track p_T&gt;100&nbsp;MeV and track p_T&gt;500&nbsp;MeV in the minimum-bias (MB) and high-multiplicity track (HMT) samples at &radic;s = 13&nbsp;TeV. The error bars and boxes represent the statistical and systematic contributions, respectively. The black solid and blue dashed curves represent the exponential fit to the parameter &kappa; for tracks with p_T &gt;100&nbsp;MeV and for tracks with p_T &gt;500&nbsp;MeV, respectively.

The parameter &kappa; describing the dependence of the source radius, R, on charged-particle scaled multiplicity, m_ch, for track p_T&gt;100&nbsp;MeV and track p_T&gt;500&nbsp;MeV in the minimum-bias (MB) and high-multiplicity track (HMT) samples at &radic;s = 13&nbsp;TeV. The error bars and boxes represent the statistical and systematic contributions, respectively. The black solid and blue dashed curves represent the exponential fit to the parameter &kappa; for tracks with p_T &gt;100&nbsp;MeV and for tracks with p_T &gt;500&nbsp;MeV, respectively.

The two-dimensional dependence on m_ch and k_T for p_T &gt; 500&nbsp;MeV for the correlation strength, &lambda;, obtained from the exponential fit to the R₂(Q) correlation functions using the MB sample for m_ch &le; 3.08 and the HMT sample for m_ch &gt; 3.08.

The two-dimensional dependence on m_ch and k_T for p_T &gt; 500&nbsp;MeV for the source radius, R, obtained from the exponential fit to the R₂(Q) correlation functions using the MB sample for m_ch &le; 3.08 and the HMT sample for m_ch &gt; 3.08.

The parameter &lambda; for p_T &gt; 500&nbsp;MeV as a function of k_T in selected low m_ch intervals. The error bars and boxes represent the statistical and systematic contributions, respectively.

The parameter &lambda; for p_T &gt; 500&nbsp;MeV as a function of k_T in selected low m_ch intervals. The error bars and boxes represent the statistical and systematic contributions, respectively.

The parameter &lambda; for p_T &gt; 500&nbsp;MeV as a function of k_T in selected high m_ch intervals. The error bars and boxes represent the statistical and systematic contributions, respectively.

The parameter &lambda; for p_T &gt; 500&nbsp;MeV as a function of k_T in selected high m_ch intervals. The error bars and boxes represent the statistical and systematic contributions, respectively.

The parameter &lambda; for p_T &gt; 500&nbsp;MeV as a function of m_ch in k_T intervals between 0.5 and 1.5&nbsp;GeV. The error bars and boxes represent the statistical and systematic contributions, respectively.

The parameter &lambda; for p_T &gt; 500&nbsp;MeV as a function of m_ch in k_T intervals between 0.5 and 1.5&nbsp;GeV. The error bars and boxes represent the statistical and systematic contributions, respectively.

The parameter R for p_T &gt; 500&nbsp;MeV as a function of k_T in selected low m_ch intervals. The error bars and boxes represent the statistical and systematic contributions, respectively.

The parameter R for p_T &gt; 500&nbsp;MeV as a function of k_T in selected low m_ch intervals. The error bars and boxes represent the statistical and systematic contributions, respectively.

The parameter R for p_T &gt; 500&nbsp;MeV as a function of k_T in selected high m_ch intervals. The error bars and boxes represent the statistical and systematic contributions, respectively.

The parameter R for p_T &gt; 500&nbsp;MeV as a function of k_T in selected high m_ch intervals. The error bars and boxes represent the statistical and systematic contributions, respectively.

The parameter R for p_T &gt; 500&nbsp;MeV as a function of m_ch in k_T intervals between 0.5 and 1.5&nbsp;GeV. The error bars and boxes represent the statistical and systematic contributions, respectively.

The parameter R for p_T &gt; 500&nbsp;MeV as a function of m_ch in k_T intervals between 0.5 and 1.5&nbsp;GeV. The error bars and boxes represent the statistical and systematic contributions, respectively.

ATLAS and CMS results for the source radius R as a function of n_ch in pp interactions at 13&nbsp;TeV. The CMS results (open circles) have been adjusted (by the CMS collaboration) to the ATLAS kinematic region&#8758; p_T &gt; 100&nbsp;MeV and |&eta;|&lt;2.5. The ATLAS uncertainties are the sum in quadrature of the statistical and asymmetric systematic uncertainties. For CMS, only the systematic uncertainties are shown since the statistical uncertainties are smaller than the marker size. The dashed blue (ATLAS) and black (CMS) lines represent the fit to &#8731;n_ch at low multiplicity, continued as dotted lines beyond the fit range. The solid green (ATLAS) and broken black (CMS) lines indicate the plateau level at high multiplicity.

ATLAS and CMS results for the source radius R as a function of n_ch in pp interactions at 13&nbsp;TeV. The CMS results (open circles) have been adjusted (by the CMS collaboration) to the ATLAS kinematic region&#8758; p_T &gt; 100&nbsp;MeV and |&eta;|&lt;2.5. The ATLAS uncertainties are the sum in quadrature of the statistical and asymmetric systematic uncertainties. For CMS, only the systematic uncertainties are shown since the statistical uncertainties are smaller than the marker size. The dashed blue (ATLAS) and black (CMS) lines represent the fit to &#8731;n_ch at low multiplicity, continued as dotted lines beyond the fit range. The solid green (ATLAS) and broken black (CMS) lines indicate the plateau level at high multiplicity.

ATLAS and CMS results for the source radius R as a function of n_ch in pp interactions at 13&nbsp;TeV. The CMS results (open circles) have been adjusted (by the CMS collaboration) to the ATLAS kinematic region&#8758; p_T &gt; 100&nbsp;MeV and |&eta;|&lt;2.5. The ATLAS uncertainties are the sum in quadrature of the statistical and asymmetric systematic uncertainties. For CMS, only the systematic uncertainties are shown since the statistical uncertainties are smaller than the marker size. The dashed blue (ATLAS) and black (CMS) lines represent the fit to &#8731;n_ch at low multiplicity, continued as dotted lines beyond the fit range. The solid green (ATLAS) and broken black (CMS) lines indicate the plateau level at high multiplicity.

ATLAS and CMS results for the source radius R as a function of &#8731;n_ch in pp interactions at 13&nbsp;TeV. The CMS results (open circles) have been adjusted (by the CMS collaboration) to the ATLAS kinematic region&#8758; p_T &gt; 100&nbsp;MeV and |&eta;|&lt;2.5. The ATLAS uncertainties are the sum in quadrature of the statistical and asymmetric systematic uncertainties. For CMS, only the systematic uncertainties are shown since the statistical uncertainties are smaller than the marker size. The dashed blue (ATLAS) and black (CMS) lines represent the fit to &#8731;n_ch at low multiplicity, continued as dotted lines beyond the fit range. The solid green (ATLAS) and broken black (CMS) lines indicate the plateau level at high multiplicity.

ATLAS and CMS results for the source radius R as a function of &#8731;n_ch in pp interactions at 13&nbsp;TeV. The CMS results (open circles) have been adjusted (by the CMS collaboration) to the ATLAS kinematic region&#8758; p_T &gt; 100&nbsp;MeV and |&eta;|&lt;2.5. The ATLAS uncertainties are the sum in quadrature of the statistical and asymmetric systematic uncertainties. For CMS, only the systematic uncertainties are shown since the statistical uncertainties are smaller than the marker size. The dashed blue (ATLAS) and black (CMS) lines represent the fit to &#8731;n_ch at low multiplicity, continued as dotted lines beyond the fit range. The solid green (ATLAS) and broken black (CMS) lines indicate the plateau level at high multiplicity.

ATLAS and CMS results for the source radius R as a function of &#8731;n_ch in pp interactions at 13&nbsp;TeV. The CMS results (open circles) have been adjusted (by the CMS collaboration) to the ATLAS kinematic region&#8758; p_T &gt; 100&nbsp;MeV and |&eta;|&lt;2.5. The ATLAS uncertainties are the sum in quadrature of the statistical and asymmetric systematic uncertainties. For CMS, only the systematic uncertainties are shown since the statistical uncertainties are smaller than the marker size. The dashed blue (ATLAS) and black (CMS) lines represent the fit to &#8731;n_ch at low multiplicity, continued as dotted lines beyond the fit range. The solid green (ATLAS) and broken black (CMS) lines indicate the plateau level at high multiplicity.

Systematic uncertainties (in percent) in the correlation strength, &lambda;, and source radius, R, for the exponential fit of the two-particle double-ratio correlation functions, R₂(Q), for p_T &gt; 100&nbsp;MeV at &radic;s= 13&nbsp;TeV for the MB and HMT events. The choice of MC generator gives rise to asymmetric uncertainties, denoted by uparrow and downarrow. This asymmetry propagates through to the cumulative uncertainty. The columns under ‘Uncertainty range’ show the range of systematic uncertainty from the fits in the various n_ch intervals.

The results of the fits to the dependencies of the correlation strength, &lambda;, and source radius, R, on the average rescaled charged-particle multiplicity, m_ch, for |&eta;| &lt; 2.5 and both p_T &gt; 100&nbsp;MeV and p_T &gt; 500&nbsp;MeV at &radic;s = 13&nbsp;TeV for the minimum-bias (MB) and the high-multiplicity track (HMT) events. The parameters &gamma; and &delta; resulting from a joint fit to the MB and HMT data are presented. The total uncertainties are shown.

The results of the fits to the dependencies of the correlation strength, &lambda;, and source radius, R, on the pair average transverse momentum, k_T, for various functional forms and for minimum-bias (MB) and high-multiplicity track (HMT) events for p_T &gt; 100&nbsp;MeV and p_T &gt; 500&nbsp;MeV at &radic;s = 13&nbsp;TeV. The total uncertainties are shown.

The Bose-Einstein correlation (BEC) parameters &lambda; and R as a function of n_ch and k_T using different MC generators in the calculation of R₂(Q). (a) &lambda; versus n_ch for MB events, (b) &lambda; versus n_ch for HMT events, (c) &lambda; versus k_T and (d) R versus k_T for MB events. The uncertainties shown are statistical. The lower panel of each plot shows the ratio of the BEC parameters obtained using EPOS LHC (red circles), Pythia 8 Monash (blue squares) and Herwig++ UE-EE-5 (green triangles) compared with the parameters obtained using Pythia 8 A2. The gray band in the lower panels is the MC systematic uncertainty, obtained as explained in the text.

The Bose-Einstein correlation (BEC) parameters &lambda; and R as a function of n_ch and k_T using different MC generators in the calculation of R₂(Q). (a) &lambda; versus n_ch for MB events, (b) &lambda; versus n_ch for HMT events, (c) &lambda; versus k_T and (d) R versus k_T for MB events. The uncertainties shown are statistical. The lower panel of each plot shows the ratio of the BEC parameters obtained using EPOS LHC (red circles), Pythia 8 Monash (blue squares) and Herwig++ UE-EE-5 (green triangles) compared with the parameters obtained using Pythia 8 A2. The gray band in the lower panels is the MC systematic uncertainty, obtained as explained in the text.

The Bose-Einstein correlation (BEC) parameters &lambda; and R as a function of n_ch and k_T using different MC generators in the calculation of R₂(Q). (a) &lambda; versus n_ch for MB events, (b) &lambda; versus n_ch for HMT events, (c) &lambda; versus k_T and (d) R versus k_T for MB events. The uncertainties shown are statistical. The lower panel of each plot shows the ratio of the BEC parameters obtained using EPOS LHC (red circles), Pythia 8 Monash (blue squares) and Herwig++ UE-EE-5 (green triangles) compared with the parameters obtained using Pythia 8 A2. The gray band in the lower panels is the MC systematic uncertainty, obtained as explained in the text.

The Bose-Einstein correlation (BEC) parameters &lambda; and R as a function of n_ch and k_T using different MC generators in the calculation of R₂(Q). (a) &lambda; versus n_ch for MB events, (b) &lambda; versus n_ch for HMT events, (c) &lambda; versus k_T and (d) R versus k_T for MB events. The uncertainties shown are statistical. The lower panel of each plot shows the ratio of the BEC parameters obtained using EPOS LHC (red circles), Pythia 8 Monash (blue squares) and Herwig++ UE-EE-5 (green triangles) compared with the parameters obtained using Pythia 8 A2. The gray band in the lower panels is the MC systematic uncertainty, obtained as explained in the text.

The single-ratio two-particle correlation functions, C₂^data(Q), for the minimum-bias (MB) events using the unlike-charge particle (UCP) pairs reference sample for n_ch - intervals&#8758; (a) 2 &lt; n_ch &le; 10, (b) 11 &lt; n_ch &le; 20, (c) 21 &lt; n_ch &le; 30, (d) 31 &lt; n_ch &le; 40, (e) 41 &lt; n_ch &le; 50, (f) 51 &lt; n_ch &le; 60, (g) 61 &lt; n_ch &le; 70, (h) 71 &lt; n_ch &le; 80 and (i) 81 &lt; n_ch &le; 90. The error bars represent the statistical uncertainties. The boxes represent the systematic uncertainties, which are the sum in quadrature of a variation of the Coulomb correction, the track reconstruction efficiency and the unfolding matrix.

The single-ratio two-particle correlation functions, C₂^data(Q), for the minimum-bias (MB) events using the unlike-charge particle (UCP) pairs reference sample for n_ch - intervals&#8758; (a) 2 &lt; n_ch &le; 10, (b) 11 &lt; n_ch &le; 20, (c) 21 &lt; n_ch &le; 30, (d) 31 &lt; n_ch &le; 40, (e) 41 &lt; n_ch &le; 50, (f) 51 &lt; n_ch &le; 60, (g) 61 &lt; n_ch &le; 70, (h) 71 &lt; n_ch &le; 80 and (i) 81 &lt; n_ch &le; 90. The error bars represent the statistical uncertainties. The boxes represent the systematic uncertainties, which are the sum in quadrature of a variation of the Coulomb correction, the track reconstruction efficiency and the unfolding matrix.

The single-ratio two-particle correlation functions, C₂^data(Q), for the minimum-bias (MB) events using the unlike-charge particle (UCP) pairs reference sample for n_ch - intervals&#8758; (a) 2 &lt; n_ch &le; 10, (b) 11 &lt; n_ch &le; 20, (c) 21 &lt; n_ch &le; 30, (d) 31 &lt; n_ch &le; 40, (e) 41 &lt; n_ch &le; 50, (f) 51 &lt; n_ch &le; 60, (g) 61 &lt; n_ch &le; 70, (h) 71 &lt; n_ch &le; 80 and (i) 81 &lt; n_ch &le; 90. The error bars represent the statistical uncertainties. The boxes represent the systematic uncertainties, which are the sum in quadrature of a variation of the Coulomb correction, the track reconstruction efficiency and the unfolding matrix.

The single-ratio two-particle correlation functions, C₂^data(Q), for the minimum-bias (MB) events using the unlike-charge particle (UCP) pairs reference sample for n_ch - intervals&#8758; (a) 2 &lt; n_ch &le; 10, (b) 11 &lt; n_ch &le; 20, (c) 21 &lt; n_ch &le; 30, (d) 31 &lt; n_ch &le; 40, (e) 41 &lt; n_ch &le; 50, (f) 51 &lt; n_ch &le; 60, (g) 61 &lt; n_ch &le; 70, (h) 71 &lt; n_ch &le; 80 and (i) 81 &lt; n_ch &le; 90. The error bars represent the statistical uncertainties. The boxes represent the systematic uncertainties, which are the sum in quadrature of a variation of the Coulomb correction, the track reconstruction efficiency and the unfolding matrix.

The single-ratio two-particle correlation functions, C₂^data(Q), for the minimum-bias (MB) events using the unlike-charge particle (UCP) pairs reference sample for n_ch - intervals&#8758; (a) 2 &lt; n_ch &le; 10, (b) 11 &lt; n_ch &le; 20, (c) 21 &lt; n_ch &le; 30, (d) 31 &lt; n_ch &le; 40, (e) 41 &lt; n_ch &le; 50, (f) 51 &lt; n_ch &le; 60, (g) 61 &lt; n_ch &le; 70, (h) 71 &lt; n_ch &le; 80 and (i) 81 &lt; n_ch &le; 90. The error bars represent the statistical uncertainties. The boxes represent the systematic uncertainties, which are the sum in quadrature of a variation of the Coulomb correction, the track reconstruction efficiency and the unfolding matrix.

The single-ratio two-particle correlation functions, C₂^data(Q), for the minimum-bias (MB) events using the unlike-charge particle (UCP) pairs reference sample for n_ch - intervals&#8758; (a) 2 &lt; n_ch &le; 10, (b) 11 &lt; n_ch &le; 20, (c) 21 &lt; n_ch &le; 30, (d) 31 &lt; n_ch &le; 40, (e) 41 &lt; n_ch &le; 50, (f) 51 &lt; n_ch &le; 60, (g) 61 &lt; n_ch &le; 70, (h) 71 &lt; n_ch &le; 80 and (i) 81 &lt; n_ch &le; 90. The error bars represent the statistical uncertainties. The boxes represent the systematic uncertainties, which are the sum in quadrature of a variation of the Coulomb correction, the track reconstruction efficiency and the unfolding matrix.

The single-ratio two-particle correlation functions, C₂^data(Q), for the minimum-bias (MB) events using the unlike-charge particle (UCP) pairs reference sample for n_ch - intervals&#8758; (a) 2 &lt; n_ch &le; 10, (b) 11 &lt; n_ch &le; 20, (c) 21 &lt; n_ch &le; 30, (d) 31 &lt; n_ch &le; 40, (e) 41 &lt; n_ch &le; 50, (f) 51 &lt; n_ch &le; 60, (g) 61 &lt; n_ch &le; 70, (h) 71 &lt; n_ch &le; 80 and (i) 81 &lt; n_ch &le; 90. The error bars represent the statistical uncertainties. The boxes represent the systematic uncertainties, which are the sum in quadrature of a variation of the Coulomb correction, the track reconstruction efficiency and the unfolding matrix.

The single-ratio two-particle correlation functions, C₂^data(Q), for the minimum-bias (MB) events using the unlike-charge particle (UCP) pairs reference sample for n_ch - intervals&#8758; (a) 2 &lt; n_ch &le; 10, (b) 11 &lt; n_ch &le; 20, (c) 21 &lt; n_ch &le; 30, (d) 31 &lt; n_ch &le; 40, (e) 41 &lt; n_ch &le; 50, (f) 51 &lt; n_ch &le; 60, (g) 61 &lt; n_ch &le; 70, (h) 71 &lt; n_ch &le; 80 and (i) 81 &lt; n_ch &le; 90. The error bars represent the statistical uncertainties. The boxes represent the systematic uncertainties, which are the sum in quadrature of a variation of the Coulomb correction, the track reconstruction efficiency and the unfolding matrix.

The single-ratio two-particle correlation functions, C₂^data(Q), for the minimum-bias (MB) events using the unlike-charge particle (UCP) pairs reference sample for n_ch - intervals&#8758; (a) 2 &lt; n_ch &le; 10, (b) 11 &lt; n_ch &le; 20, (c) 21 &lt; n_ch &le; 30, (d) 31 &lt; n_ch &le; 40, (e) 41 &lt; n_ch &le; 50, (f) 51 &lt; n_ch &le; 60, (g) 61 &lt; n_ch &le; 70, (h) 71 &lt; n_ch &le; 80 and (i) 81 &lt; n_ch &le; 90. The error bars represent the statistical uncertainties. The boxes represent the systematic uncertainties, which are the sum in quadrature of a variation of the Coulomb correction, the track reconstruction efficiency and the unfolding matrix.

The single-ratio two-particle correlation functions, C₂^data(Q), for the minimum-bias (MB) events using the unlike-charge particle (UCP) pairs reference sample for n_ch - intervals&#8758; (a) 91 &lt; n_ch &le; 100, (b) 101 &lt; n_ch &le; 125, (c) 126 &lt; n_ch &le; 150, (d) 151 &lt; n_ch &le; 200, (e) 201 &lt; n_ch &le; 250. The error bars represent the statistical uncertainties. The boxes represent the systematic uncertainties, which are the sum in quadrature of a variation of the Coulomb correction, the track reconstruction efficiency and the unfolding matrix.

The single-ratio two-particle correlation functions, C₂^data(Q), for the minimum-bias (MB) events using the unlike-charge particle (UCP) pairs reference sample for n_ch - intervals&#8758; (a) 91 &lt; n_ch &le; 100, (b) 101 &lt; n_ch &le; 125, (c) 126 &lt; n_ch &le; 150, (d) 151 &lt; n_ch &le; 200, (e) 201 &lt; n_ch &le; 250. The error bars represent the statistical uncertainties. The boxes represent the systematic uncertainties, which are the sum in quadrature of a variation of the Coulomb correction, the track reconstruction efficiency and the unfolding matrix.

The single-ratio two-particle correlation functions, C₂^data(Q), for the minimum-bias (MB) events using the unlike-charge particle (UCP) pairs reference sample for n_ch - intervals&#8758; (a) 91 &lt; n_ch &le; 100, (b) 101 &lt; n_ch &le; 125, (c) 126 &lt; n_ch &le; 150, (d) 151 &lt; n_ch &le; 200, (e) 201 &lt; n_ch &le; 250. The error bars represent the statistical uncertainties. The boxes represent the systematic uncertainties, which are the sum in quadrature of a variation of the Coulomb correction, the track reconstruction efficiency and the unfolding matrix.

The single-ratio two-particle correlation functions, C₂^data(Q), for the minimum-bias (MB) events using the unlike-charge particle (UCP) pairs reference sample for n_ch - intervals&#8758; (a) 91 &lt; n_ch &le; 100, (b) 101 &lt; n_ch &le; 125, (c) 126 &lt; n_ch &le; 150, (d) 151 &lt; n_ch &le; 200, (e) 201 &lt; n_ch &le; 250. The error bars represent the statistical uncertainties. The boxes represent the systematic uncertainties, which are the sum in quadrature of a variation of the Coulomb correction, the track reconstruction efficiency and the unfolding matrix.

The single-ratio two-particle correlation functions, C₂^data(Q), for the minimum-bias (MB) events using the unlike-charge particle (UCP) pairs reference sample for n_ch - intervals&#8758; (a) 91 &lt; n_ch &le; 100, (b) 101 &lt; n_ch &le; 125, (c) 126 &lt; n_ch &le; 150, (d) 151 &lt; n_ch &le; 200, (e) 201 &lt; n_ch &le; 250. The error bars represent the statistical uncertainties. The boxes represent the systematic uncertainties, which are the sum in quadrature of a variation of the Coulomb correction, the track reconstruction efficiency and the unfolding matrix.

The single-ratio two-particle correlation functions, C₂^data(Q), for the high-multiplicity track (HMT) events using the unlike-charge particle (UCP) pairs reference sample for n_ch - intervals&#8758; (a) 101 &lt; n_ch &le; 110, (b) 111 &lt; n_ch &le; 120, (c) 121 &lt; n_ch &le; 130, (d) 131 &lt; n_ch &le; 140, (e) 141 &lt; n_ch &le; 155, (f) 156 &lt; n_ch &le; 175, (g) 176 &lt; n_ch &le; 200, (h) 201 &lt; n_ch &le; 230 and (i) 231 &lt; n_ch &le; 300. The error bars represent the statistical uncertainties. The boxes represent the systematic uncertainties, which are the sum in quadrature of a variation of the Coulomb correction, the track reconstruction efficiency and the unfolding matrix.

The single-ratio two-particle correlation functions, C₂^data(Q), for the high-multiplicity track (HMT) events using the unlike-charge particle (UCP) pairs reference sample for n_ch - intervals&#8758; (a) 101 &lt; n_ch &le; 110, (b) 111 &lt; n_ch &le; 120, (c) 121 &lt; n_ch &le; 130, (d) 131 &lt; n_ch &le; 140, (e) 141 &lt; n_ch &le; 155, (f) 156 &lt; n_ch &le; 175, (g) 176 &lt; n_ch &le; 200, (h) 201 &lt; n_ch &le; 230 and (i) 231 &lt; n_ch &le; 300. The error bars represent the statistical uncertainties. The boxes represent the systematic uncertainties, which are the sum in quadrature of a variation of the Coulomb correction, the track reconstruction efficiency and the unfolding matrix.

The single-ratio two-particle correlation functions, C₂^data(Q), for the high-multiplicity track (HMT) events using the unlike-charge particle (UCP) pairs reference sample for n_ch - intervals&#8758; (a) 101 &lt; n_ch &le; 110, (b) 111 &lt; n_ch &le; 120, (c) 121 &lt; n_ch &le; 130, (d) 131 &lt; n_ch &le; 140, (e) 141 &lt; n_ch &le; 155, (f) 156 &lt; n_ch &le; 175, (g) 176 &lt; n_ch &le; 200, (h) 201 &lt; n_ch &le; 230 and (i) 231 &lt; n_ch &le; 300. The error bars represent the statistical uncertainties. The boxes represent the systematic uncertainties, which are the sum in quadrature of a variation of the Coulomb correction, the track reconstruction efficiency and the unfolding matrix.

The single-ratio two-particle correlation functions, C₂^data(Q), for the high-multiplicity track (HMT) events using the unlike-charge particle (UCP) pairs reference sample for n_ch - intervals&#8758; (a) 101 &lt; n_ch &le; 110, (b) 111 &lt; n_ch &le; 120, (c) 121 &lt; n_ch &le; 130, (d) 131 &lt; n_ch &le; 140, (e) 141 &lt; n_ch &le; 155, (f) 156 &lt; n_ch &le; 175, (g) 176 &lt; n_ch &le; 200, (h) 201 &lt; n_ch &le; 230 and (i) 231 &lt; n_ch &le; 300. The error bars represent the statistical uncertainties. The boxes represent the systematic uncertainties, which are the sum in quadrature of a variation of the Coulomb correction, the track reconstruction efficiency and the unfolding matrix.

The single-ratio two-particle correlation functions, C₂^data(Q), for the high-multiplicity track (HMT) events using the unlike-charge particle (UCP) pairs reference sample for n_ch - intervals&#8758; (a) 101 &lt; n_ch &le; 110, (b) 111 &lt; n_ch &le; 120, (c) 121 &lt; n_ch &le; 130, (d) 131 &lt; n_ch &le; 140, (e) 141 &lt; n_ch &le; 155, (f) 156 &lt; n_ch &le; 175, (g) 176 &lt; n_ch &le; 200, (h) 201 &lt; n_ch &le; 230 and (i) 231 &lt; n_ch &le; 300. The error bars represent the statistical uncertainties. The boxes represent the systematic uncertainties, which are the sum in quadrature of a variation of the Coulomb correction, the track reconstruction efficiency and the unfolding matrix.

The single-ratio two-particle correlation functions, C₂^data(Q), for the high-multiplicity track (HMT) events using the unlike-charge particle (UCP) pairs reference sample for n_ch - intervals&#8758; (a) 101 &lt; n_ch &le; 110, (b) 111 &lt; n_ch &le; 120, (c) 121 &lt; n_ch &le; 130, (d) 131 &lt; n_ch &le; 140, (e) 141 &lt; n_ch &le; 155, (f) 156 &lt; n_ch &le; 175, (g) 176 &lt; n_ch &le; 200, (h) 201 &lt; n_ch &le; 230 and (i) 231 &lt; n_ch &le; 300. The error bars represent the statistical uncertainties. The boxes represent the systematic uncertainties, which are the sum in quadrature of a variation of the Coulomb correction, the track reconstruction efficiency and the unfolding matrix.

The single-ratio two-particle correlation functions, C₂^data(Q), for the high-multiplicity track (HMT) events using the unlike-charge particle (UCP) pairs reference sample for n_ch - intervals&#8758; (a) 101 &lt; n_ch &le; 110, (b) 111 &lt; n_ch &le; 120, (c) 121 &lt; n_ch &le; 130, (d) 131 &lt; n_ch &le; 140, (e) 141 &lt; n_ch &le; 155, (f) 156 &lt; n_ch &le; 175, (g) 176 &lt; n_ch &le; 200, (h) 201 &lt; n_ch &le; 230 and (i) 231 &lt; n_ch &le; 300. The error bars represent the statistical uncertainties. The boxes represent the systematic uncertainties, which are the sum in quadrature of a variation of the Coulomb correction, the track reconstruction efficiency and the unfolding matrix.

The single-ratio two-particle correlation functions, C₂^data(Q), for the high-multiplicity track (HMT) events using the unlike-charge particle (UCP) pairs reference sample for n_ch - intervals&#8758; (a) 101 &lt; n_ch &le; 110, (b) 111 &lt; n_ch &le; 120, (c) 121 &lt; n_ch &le; 130, (d) 131 &lt; n_ch &le; 140, (e) 141 &lt; n_ch &le; 155, (f) 156 &lt; n_ch &le; 175, (g) 176 &lt; n_ch &le; 200, (h) 201 &lt; n_ch &le; 230 and (i) 231 &lt; n_ch &le; 300. The error bars represent the statistical uncertainties. The boxes represent the systematic uncertainties, which are the sum in quadrature of a variation of the Coulomb correction, the track reconstruction efficiency and the unfolding matrix.

The single-ratio two-particle correlation functions, C₂^data(Q), for the high-multiplicity track (HMT) events using the unlike-charge particle (UCP) pairs reference sample for n_ch - intervals&#8758; (a) 101 &lt; n_ch &le; 110, (b) 111 &lt; n_ch &le; 120, (c) 121 &lt; n_ch &le; 130, (d) 131 &lt; n_ch &le; 140, (e) 141 &lt; n_ch &le; 155, (f) 156 &lt; n_ch &le; 175, (g) 176 &lt; n_ch &le; 200, (h) 201 &lt; n_ch &le; 230 and (i) 231 &lt; n_ch &le; 300. The error bars represent the statistical uncertainties. The boxes represent the systematic uncertainties, which are the sum in quadrature of a variation of the Coulomb correction, the track reconstruction efficiency and the unfolding matrix.

The single-ratio two-particle correlation functions, C₂^data(Q), for the minimum-bias (MB) events using the unlike-charge particle (UCP) pairs reference sample k_T - intervals&#8758; (a) 100 &lt; k_T &le; 200&nbsp;MeV, (b) 200 &lt; k_T &le; 300&nbsp;MeV, (c) 300 &lt; k_T &le; 400&nbsp;MeV, (d) 400 &lt; k_T &le; 500&nbsp;MeV, (e) 500 &lt; k_T &le; 600&nbsp;MeV, (f) 600 &lt; k_T &le; 700&nbsp;MeV, (g) 700 &lt; k_T &le; 1000&nbsp;MeV, and (h) 1000 &lt; k_T &le; 1500&nbsp;MeV. The error bars represent the statistical uncertainties. The boxes represent the systematic uncertainties, which are the sum in quadrature of a variation of the Coulomb correction, the track reconstruction efficiency and the unfolding matrix.

The single-ratio two-particle correlation functions, C₂^data(Q), for the minimum-bias (MB) events using the unlike-charge particle (UCP) pairs reference sample k_T - intervals&#8758; (a) 100 &lt; k_T &le; 200&nbsp;MeV, (b) 200 &lt; k_T &le; 300&nbsp;MeV, (c) 300 &lt; k_T &le; 400&nbsp;MeV, (d) 400 &lt; k_T &le; 500&nbsp;MeV, (e) 500 &lt; k_T &le; 600&nbsp;MeV, (f) 600 &lt; k_T &le; 700&nbsp;MeV, (g) 700 &lt; k_T &le; 1000&nbsp;MeV, and (h) 1000 &lt; k_T &le; 1500&nbsp;MeV. The error bars represent the statistical uncertainties. The boxes represent the systematic uncertainties, which are the sum in quadrature of a variation of the Coulomb correction, the track reconstruction efficiency and the unfolding matrix.

The single-ratio two-particle correlation functions, C₂^data(Q), for the minimum-bias (MB) events using the unlike-charge particle (UCP) pairs reference sample k_T - intervals&#8758; (a) 100 &lt; k_T &le; 200&nbsp;MeV, (b) 200 &lt; k_T &le; 300&nbsp;MeV, (c) 300 &lt; k_T &le; 400&nbsp;MeV, (d) 400 &lt; k_T &le; 500&nbsp;MeV, (e) 500 &lt; k_T &le; 600&nbsp;MeV, (f) 600 &lt; k_T &le; 700&nbsp;MeV, (g) 700 &lt; k_T &le; 1000&nbsp;MeV, and (h) 1000 &lt; k_T &le; 1500&nbsp;MeV. The error bars represent the statistical uncertainties. The boxes represent the systematic uncertainties, which are the sum in quadrature of a variation of the Coulomb correction, the track reconstruction efficiency and the unfolding matrix.

The single-ratio two-particle correlation functions, C₂^data(Q), for the minimum-bias (MB) events using the unlike-charge particle (UCP) pairs reference sample k_T - intervals&#8758; (a) 100 &lt; k_T &le; 200&nbsp;MeV, (b) 200 &lt; k_T &le; 300&nbsp;MeV, (c) 300 &lt; k_T &le; 400&nbsp;MeV, (d) 400 &lt; k_T &le; 500&nbsp;MeV, (e) 500 &lt; k_T &le; 600&nbsp;MeV, (f) 600 &lt; k_T &le; 700&nbsp;MeV, (g) 700 &lt; k_T &le; 1000&nbsp;MeV, and (h) 1000 &lt; k_T &le; 1500&nbsp;MeV. The error bars represent the statistical uncertainties. The boxes represent the systematic uncertainties, which are the sum in quadrature of a variation of the Coulomb correction, the track reconstruction efficiency and the unfolding matrix.

The single-ratio two-particle correlation functions, C₂^data(Q), for the minimum-bias (MB) events using the unlike-charge particle (UCP) pairs reference sample k_T - intervals&#8758; (a) 100 &lt; k_T &le; 200&nbsp;MeV, (b) 200 &lt; k_T &le; 300&nbsp;MeV, (c) 300 &lt; k_T &le; 400&nbsp;MeV, (d) 400 &lt; k_T &le; 500&nbsp;MeV, (e) 500 &lt; k_T &le; 600&nbsp;MeV, (f) 600 &lt; k_T &le; 700&nbsp;MeV, (g) 700 &lt; k_T &le; 1000&nbsp;MeV, and (h) 1000 &lt; k_T &le; 1500&nbsp;MeV. The error bars represent the statistical uncertainties. The boxes represent the systematic uncertainties, which are the sum in quadrature of a variation of the Coulomb correction, the track reconstruction efficiency and the unfolding matrix.

The single-ratio two-particle correlation functions, C₂^data(Q), for the minimum-bias (MB) events using the unlike-charge particle (UCP) pairs reference sample k_T - intervals&#8758; (a) 100 &lt; k_T &le; 200&nbsp;MeV, (b) 200 &lt; k_T &le; 300&nbsp;MeV, (c) 300 &lt; k_T &le; 400&nbsp;MeV, (d) 400 &lt; k_T &le; 500&nbsp;MeV, (e) 500 &lt; k_T &le; 600&nbsp;MeV, (f) 600 &lt; k_T &le; 700&nbsp;MeV, (g) 700 &lt; k_T &le; 1000&nbsp;MeV, and (h) 1000 &lt; k_T &le; 1500&nbsp;MeV. The error bars represent the statistical uncertainties. The boxes represent the systematic uncertainties, which are the sum in quadrature of a variation of the Coulomb correction, the track reconstruction efficiency and the unfolding matrix.

The single-ratio two-particle correlation functions, C₂^data(Q), for the minimum-bias (MB) events using the unlike-charge particle (UCP) pairs reference sample k_T - intervals&#8758; (a) 100 &lt; k_T &le; 200&nbsp;MeV, (b) 200 &lt; k_T &le; 300&nbsp;MeV, (c) 300 &lt; k_T &le; 400&nbsp;MeV, (d) 400 &lt; k_T &le; 500&nbsp;MeV, (e) 500 &lt; k_T &le; 600&nbsp;MeV, (f) 600 &lt; k_T &le; 700&nbsp;MeV, (g) 700 &lt; k_T &le; 1000&nbsp;MeV, and (h) 1000 &lt; k_T &le; 1500&nbsp;MeV. The error bars represent the statistical uncertainties. The boxes represent the systematic uncertainties, which are the sum in quadrature of a variation of the Coulomb correction, the track reconstruction efficiency and the unfolding matrix.

The single-ratio two-particle correlation functions, C₂^data(Q), for the minimum-bias (MB) events using the unlike-charge particle (UCP) pairs reference sample k_T - intervals&#8758; (a) 100 &lt; k_T &le; 200&nbsp;MeV, (b) 200 &lt; k_T &le; 300&nbsp;MeV, (c) 300 &lt; k_T &le; 400&nbsp;MeV, (d) 400 &lt; k_T &le; 500&nbsp;MeV, (e) 500 &lt; k_T &le; 600&nbsp;MeV, (f) 600 &lt; k_T &le; 700&nbsp;MeV, (g) 700 &lt; k_T &le; 1000&nbsp;MeV, and (h) 1000 &lt; k_T &le; 1500&nbsp;MeV. The error bars represent the statistical uncertainties. The boxes represent the systematic uncertainties, which are the sum in quadrature of a variation of the Coulomb correction, the track reconstruction efficiency and the unfolding matrix.

The single-ratio two-particle correlation functions, C₂^data(Q), for the high-multiplicity track (HMT) events using the unlike-charge particle (UCP) pairs reference sample for k_T - intervals&#8758; (a) 100 &lt; k_T &le; 200&nbsp;MeV, (b) 200 &lt; k_T &le; 300&nbsp;MeV, (c) 300 &lt; k_T &le; 400&nbsp;MeV, (d) 400 &lt; k_T &le; 500&nbsp;MeV, (e) 500 &lt; k_T &le; 600&nbsp;MeV, (f) 600 &lt; k_T &le; 700&nbsp;MeV, (g) 700 &lt; k_T &le; 1000&nbsp;MeV, and (h) 1000 &lt; k_T &le; 1500&nbsp;MeV. The error bars represent the statistical uncertainties. The boxes represent the systematic uncertainties, which are the sum in quadrature of a variation of the Coulomb correction, the track reconstruction efficiency and the unfolding matrix.

The single-ratio two-particle correlation functions, C₂^data(Q), for the high-multiplicity track (HMT) events using the unlike-charge particle (UCP) pairs reference sample for k_T - intervals&#8758; (a) 100 &lt; k_T &le; 200&nbsp;MeV, (b) 200 &lt; k_T &le; 300&nbsp;MeV, (c) 300 &lt; k_T &le; 400&nbsp;MeV, (d) 400 &lt; k_T &le; 500&nbsp;MeV, (e) 500 &lt; k_T &le; 600&nbsp;MeV, (f) 600 &lt; k_T &le; 700&nbsp;MeV, (g) 700 &lt; k_T &le; 1000&nbsp;MeV, and (h) 1000 &lt; k_T &le; 1500&nbsp;MeV. The error bars represent the statistical uncertainties. The boxes represent the systematic uncertainties, which are the sum in quadrature of a variation of the Coulomb correction, the track reconstruction efficiency and the unfolding matrix.

The single-ratio two-particle correlation functions, C₂^data(Q), for the high-multiplicity track (HMT) events using the unlike-charge particle (UCP) pairs reference sample for k_T - intervals&#8758; (a) 100 &lt; k_T &le; 200&nbsp;MeV, (b) 200 &lt; k_T &le; 300&nbsp;MeV, (c) 300 &lt; k_T &le; 400&nbsp;MeV, (d) 400 &lt; k_T &le; 500&nbsp;MeV, (e) 500 &lt; k_T &le; 600&nbsp;MeV, (f) 600 &lt; k_T &le; 700&nbsp;MeV, (g) 700 &lt; k_T &le; 1000&nbsp;MeV, and (h) 1000 &lt; k_T &le; 1500&nbsp;MeV. The error bars represent the statistical uncertainties. The boxes represent the systematic uncertainties, which are the sum in quadrature of a variation of the Coulomb correction, the track reconstruction efficiency and the unfolding matrix.

The single-ratio two-particle correlation functions, C₂^data(Q), for the high-multiplicity track (HMT) events using the unlike-charge particle (UCP) pairs reference sample for k_T - intervals&#8758; (a) 100 &lt; k_T &le; 200&nbsp;MeV, (b) 200 &lt; k_T &le; 300&nbsp;MeV, (c) 300 &lt; k_T &le; 400&nbsp;MeV, (d) 400 &lt; k_T &le; 500&nbsp;MeV, (e) 500 &lt; k_T &le; 600&nbsp;MeV, (f) 600 &lt; k_T &le; 700&nbsp;MeV, (g) 700 &lt; k_T &le; 1000&nbsp;MeV, and (h) 1000 &lt; k_T &le; 1500&nbsp;MeV. The error bars represent the statistical uncertainties. The boxes represent the systematic uncertainties, which are the sum in quadrature of a variation of the Coulomb correction, the track reconstruction efficiency and the unfolding matrix.

The single-ratio two-particle correlation functions, C₂^data(Q), for the high-multiplicity track (HMT) events using the unlike-charge particle (UCP) pairs reference sample for k_T - intervals&#8758; (a) 100 &lt; k_T &le; 200&nbsp;MeV, (b) 200 &lt; k_T &le; 300&nbsp;MeV, (c) 300 &lt; k_T &le; 400&nbsp;MeV, (d) 400 &lt; k_T &le; 500&nbsp;MeV, (e) 500 &lt; k_T &le; 600&nbsp;MeV, (f) 600 &lt; k_T &le; 700&nbsp;MeV, (g) 700 &lt; k_T &le; 1000&nbsp;MeV, and (h) 1000 &lt; k_T &le; 1500&nbsp;MeV. The error bars represent the statistical uncertainties. The boxes represent the systematic uncertainties, which are the sum in quadrature of a variation of the Coulomb correction, the track reconstruction efficiency and the unfolding matrix.

The single-ratio two-particle correlation functions, C₂^data(Q), for the high-multiplicity track (HMT) events using the unlike-charge particle (UCP) pairs reference sample for k_T - intervals&#8758; (a) 100 &lt; k_T &le; 200&nbsp;MeV, (b) 200 &lt; k_T &le; 300&nbsp;MeV, (c) 300 &lt; k_T &le; 400&nbsp;MeV, (d) 400 &lt; k_T &le; 500&nbsp;MeV, (e) 500 &lt; k_T &le; 600&nbsp;MeV, (f) 600 &lt; k_T &le; 700&nbsp;MeV, (g) 700 &lt; k_T &le; 1000&nbsp;MeV, and (h) 1000 &lt; k_T &le; 1500&nbsp;MeV. The error bars represent the statistical uncertainties. The boxes represent the systematic uncertainties, which are the sum in quadrature of a variation of the Coulomb correction, the track reconstruction efficiency and the unfolding matrix.

The single-ratio two-particle correlation functions, C₂^data(Q), for the high-multiplicity track (HMT) events using the unlike-charge particle (UCP) pairs reference sample for k_T - intervals&#8758; (a) 100 &lt; k_T &le; 200&nbsp;MeV, (b) 200 &lt; k_T &le; 300&nbsp;MeV, (c) 300 &lt; k_T &le; 400&nbsp;MeV, (d) 400 &lt; k_T &le; 500&nbsp;MeV, (e) 500 &lt; k_T &le; 600&nbsp;MeV, (f) 600 &lt; k_T &le; 700&nbsp;MeV, (g) 700 &lt; k_T &le; 1000&nbsp;MeV, and (h) 1000 &lt; k_T &le; 1500&nbsp;MeV. The error bars represent the statistical uncertainties. The boxes represent the systematic uncertainties, which are the sum in quadrature of a variation of the Coulomb correction, the track reconstruction efficiency and the unfolding matrix.

The single-ratio two-particle correlation functions, C₂^data(Q), for the high-multiplicity track (HMT) events using the unlike-charge particle (UCP) pairs reference sample for k_T - intervals&#8758; (a) 100 &lt; k_T &le; 200&nbsp;MeV, (b) 200 &lt; k_T &le; 300&nbsp;MeV, (c) 300 &lt; k_T &le; 400&nbsp;MeV, (d) 400 &lt; k_T &le; 500&nbsp;MeV, (e) 500 &lt; k_T &le; 600&nbsp;MeV, (f) 600 &lt; k_T &le; 700&nbsp;MeV, (g) 700 &lt; k_T &le; 1000&nbsp;MeV, and (h) 1000 &lt; k_T &le; 1500&nbsp;MeV. The error bars represent the statistical uncertainties. The boxes represent the systematic uncertainties, which are the sum in quadrature of a variation of the Coulomb correction, the track reconstruction efficiency and the unfolding matrix.

The single-ratio two-particle correlation functions, C₂^data(Q), at 7 TeV for the minimum-bias (MB) events using the unlike-charge particle (UCP) pairs reference sample for n_ch - intervals&#8758; (a) 2 &lt; n_ch &le; 9, (b) 10 &lt; n_ch &le; 18, (c) 19 &lt; n_ch &le; 27, (d) 28 &lt; n_ch &le; 36, (e) 37 &lt; n_ch &le; 45, (f) 46 &lt; n_ch &le; 54, (g) 55 &lt; n_ch &le; 63, (h) 64 &lt; n_ch &le; 72, (i) 73 &lt; n_ch &le; 81, (j) 82 &lt; n_ch &le; 90, (k) 91 &lt; n_ch &le; 113, and (l) 114 &lt; n_ch &le; 136. The error bars represent the statistical uncertainties. The boxes represent the systematic uncertainties, which are the sum in quadrature of a variation of the Coulomb correction, the track reconstruction efficiency and the unfolding matrix.

The single-ratio two-particle correlation functions, C₂^data(Q), at 7 TeV for the minimum-bias (MB) events using the unlike-charge particle (UCP) pairs reference sample for n_ch - intervals&#8758; (a) 2 &lt; n_ch &le; 9, (b) 10 &lt; n_ch &le; 18, (c) 19 &lt; n_ch &le; 27, (d) 28 &lt; n_ch &le; 36, (e) 37 &lt; n_ch &le; 45, (f) 46 &lt; n_ch &le; 54, (g) 55 &lt; n_ch &le; 63, (h) 64 &lt; n_ch &le; 72, (i) 73 &lt; n_ch &le; 81, (j) 82 &lt; n_ch &le; 90, (k) 91 &lt; n_ch &le; 113, and (l) 114 &lt; n_ch &le; 136. The error bars represent the statistical uncertainties. The boxes represent the systematic uncertainties, which are the sum in quadrature of a variation of the Coulomb correction, the track reconstruction efficiency and the unfolding matrix.

The single-ratio two-particle correlation functions, C₂^data(Q), at 7 TeV for the minimum-bias (MB) events using the unlike-charge particle (UCP) pairs reference sample for n_ch - intervals&#8758; (a) 2 &lt; n_ch &le; 9, (b) 10 &lt; n_ch &le; 18, (c) 19 &lt; n_ch &le; 27, (d) 28 &lt; n_ch &le; 36, (e) 37 &lt; n_ch &le; 45, (f) 46 &lt; n_ch &le; 54, (g) 55 &lt; n_ch &le; 63, (h) 64 &lt; n_ch &le; 72, (i) 73 &lt; n_ch &le; 81, (j) 82 &lt; n_ch &le; 90, (k) 91 &lt; n_ch &le; 113, and (l) 114 &lt; n_ch &le; 136. The error bars represent the statistical uncertainties. The boxes represent the systematic uncertainties, which are the sum in quadrature of a variation of the Coulomb correction, the track reconstruction efficiency and the unfolding matrix.

The single-ratio two-particle correlation functions, C₂^data(Q), at 7 TeV for the minimum-bias (MB) events using the unlike-charge particle (UCP) pairs reference sample for n_ch - intervals&#8758; (a) 2 &lt; n_ch &le; 9, (b) 10 &lt; n_ch &le; 18, (c) 19 &lt; n_ch &le; 27, (d) 28 &lt; n_ch &le; 36, (e) 37 &lt; n_ch &le; 45, (f) 46 &lt; n_ch &le; 54, (g) 55 &lt; n_ch &le; 63, (h) 64 &lt; n_ch &le; 72, (i) 73 &lt; n_ch &le; 81, (j) 82 &lt; n_ch &le; 90, (k) 91 &lt; n_ch &le; 113, and (l) 114 &lt; n_ch &le; 136. The error bars represent the statistical uncertainties. The boxes represent the systematic uncertainties, which are the sum in quadrature of a variation of the Coulomb correction, the track reconstruction efficiency and the unfolding matrix.

The single-ratio two-particle correlation functions, C₂^data(Q), at 7 TeV for the minimum-bias (MB) events using the unlike-charge particle (UCP) pairs reference sample for n_ch - intervals&#8758; (a) 2 &lt; n_ch &le; 9, (b) 10 &lt; n_ch &le; 18, (c) 19 &lt; n_ch &le; 27, (d) 28 &lt; n_ch &le; 36, (e) 37 &lt; n_ch &le; 45, (f) 46 &lt; n_ch &le; 54, (g) 55 &lt; n_ch &le; 63, (h) 64 &lt; n_ch &le; 72, (i) 73 &lt; n_ch &le; 81, (j) 82 &lt; n_ch &le; 90, (k) 91 &lt; n_ch &le; 113, and (l) 114 &lt; n_ch &le; 136. The error bars represent the statistical uncertainties. The boxes represent the systematic uncertainties, which are the sum in quadrature of a variation of the Coulomb correction, the track reconstruction efficiency and the unfolding matrix.

The single-ratio two-particle correlation functions, C₂^data(Q), at 7 TeV for the minimum-bias (MB) events using the unlike-charge particle (UCP) pairs reference sample for n_ch - intervals&#8758; (a) 2 &lt; n_ch &le; 9, (b) 10 &lt; n_ch &le; 18, (c) 19 &lt; n_ch &le; 27, (d) 28 &lt; n_ch &le; 36, (e) 37 &lt; n_ch &le; 45, (f) 46 &lt; n_ch &le; 54, (g) 55 &lt; n_ch &le; 63, (h) 64 &lt; n_ch &le; 72, (i) 73 &lt; n_ch &le; 81, (j) 82 &lt; n_ch &le; 90, (k) 91 &lt; n_ch &le; 113, and (l) 114 &lt; n_ch &le; 136. The error bars represent the statistical uncertainties. The boxes represent the systematic uncertainties, which are the sum in quadrature of a variation of the Coulomb correction, the track reconstruction efficiency and the unfolding matrix.

The single-ratio two-particle correlation functions, C₂^data(Q), at 7 TeV for the minimum-bias (MB) events using the unlike-charge particle (UCP) pairs reference sample for n_ch - intervals&#8758; (a) 2 &lt; n_ch &le; 9, (b) 10 &lt; n_ch &le; 18, (c) 19 &lt; n_ch &le; 27, (d) 28 &lt; n_ch &le; 36, (e) 37 &lt; n_ch &le; 45, (f) 46 &lt; n_ch &le; 54, (g) 55 &lt; n_ch &le; 63, (h) 64 &lt; n_ch &le; 72, (i) 73 &lt; n_ch &le; 81, (j) 82 &lt; n_ch &le; 90, (k) 91 &lt; n_ch &le; 113, and (l) 114 &lt; n_ch &le; 136. The error bars represent the statistical uncertainties. The boxes represent the systematic uncertainties, which are the sum in quadrature of a variation of the Coulomb correction, the track reconstruction efficiency and the unfolding matrix.

The single-ratio two-particle correlation functions, C₂^data(Q), at 7 TeV for the minimum-bias (MB) events using the unlike-charge particle (UCP) pairs reference sample for n_ch - intervals&#8758; (a) 2 &lt; n_ch &le; 9, (b) 10 &lt; n_ch &le; 18, (c) 19 &lt; n_ch &le; 27, (d) 28 &lt; n_ch &le; 36, (e) 37 &lt; n_ch &le; 45, (f) 46 &lt; n_ch &le; 54, (g) 55 &lt; n_ch &le; 63, (h) 64 &lt; n_ch &le; 72, (i) 73 &lt; n_ch &le; 81, (j) 82 &lt; n_ch &le; 90, (k) 91 &lt; n_ch &le; 113, and (l) 114 &lt; n_ch &le; 136. The error bars represent the statistical uncertainties. The boxes represent the systematic uncertainties, which are the sum in quadrature of a variation of the Coulomb correction, the track reconstruction efficiency and the unfolding matrix.

The single-ratio two-particle correlation functions, C₂^data(Q), at 7 TeV for the minimum-bias (MB) events using the unlike-charge particle (UCP) pairs reference sample for n_ch - intervals&#8758; (a) 2 &lt; n_ch &le; 9, (b) 10 &lt; n_ch &le; 18, (c) 19 &lt; n_ch &le; 27, (d) 28 &lt; n_ch &le; 36, (e) 37 &lt; n_ch &le; 45, (f) 46 &lt; n_ch &le; 54, (g) 55 &lt; n_ch &le; 63, (h) 64 &lt; n_ch &le; 72, (i) 73 &lt; n_ch &le; 81, (j) 82 &lt; n_ch &le; 90, (k) 91 &lt; n_ch &le; 113, and (l) 114 &lt; n_ch &le; 136. The error bars represent the statistical uncertainties. The boxes represent the systematic uncertainties, which are the sum in quadrature of a variation of the Coulomb correction, the track reconstruction efficiency and the unfolding matrix.

The single-ratio two-particle correlation functions, C₂^data(Q), at 7 TeV for the minimum-bias (MB) events using the unlike-charge particle (UCP) pairs reference sample for n_ch - intervals&#8758; (a) 2 &lt; n_ch &le; 9, (b) 10 &lt; n_ch &le; 18, (c) 19 &lt; n_ch &le; 27, (d) 28 &lt; n_ch &le; 36, (e) 37 &lt; n_ch &le; 45, (f) 46 &lt; n_ch &le; 54, (g) 55 &lt; n_ch &le; 63, (h) 64 &lt; n_ch &le; 72, (i) 73 &lt; n_ch &le; 81, (j) 82 &lt; n_ch &le; 90, (k) 91 &lt; n_ch &le; 113, and (l) 114 &lt; n_ch &le; 136. The error bars represent the statistical uncertainties. The boxes represent the systematic uncertainties, which are the sum in quadrature of a variation of the Coulomb correction, the track reconstruction efficiency and the unfolding matrix.

The single-ratio two-particle correlation functions, C₂^data(Q), at 7 TeV for the minimum-bias (MB) events using the unlike-charge particle (UCP) pairs reference sample for n_ch - intervals&#8758; (a) 2 &lt; n_ch &le; 9, (b) 10 &lt; n_ch &le; 18, (c) 19 &lt; n_ch &le; 27, (d) 28 &lt; n_ch &le; 36, (e) 37 &lt; n_ch &le; 45, (f) 46 &lt; n_ch &le; 54, (g) 55 &lt; n_ch &le; 63, (h) 64 &lt; n_ch &le; 72, (i) 73 &lt; n_ch &le; 81, (j) 82 &lt; n_ch &le; 90, (k) 91 &lt; n_ch &le; 113, and (l) 114 &lt; n_ch &le; 136. The error bars represent the statistical uncertainties. The boxes represent the systematic uncertainties, which are the sum in quadrature of a variation of the Coulomb correction, the track reconstruction efficiency and the unfolding matrix.

The single-ratio two-particle correlation functions, C₂^data(Q), at 7 TeV for the minimum-bias (MB) events using the unlike-charge particle (UCP) pairs reference sample for n_ch - intervals&#8758; (a) 2 &lt; n_ch &le; 9, (b) 10 &lt; n_ch &le; 18, (c) 19 &lt; n_ch &le; 27, (d) 28 &lt; n_ch &le; 36, (e) 37 &lt; n_ch &le; 45, (f) 46 &lt; n_ch &le; 54, (g) 55 &lt; n_ch &le; 63, (h) 64 &lt; n_ch &le; 72, (i) 73 &lt; n_ch &le; 81, (j) 82 &lt; n_ch &le; 90, (k) 91 &lt; n_ch &le; 113, and (l) 114 &lt; n_ch &le; 136. The error bars represent the statistical uncertainties. The boxes represent the systematic uncertainties, which are the sum in quadrature of a variation of the Coulomb correction, the track reconstruction efficiency and the unfolding matrix.

The single-ratio two-particle correlation functions, C₂^data(Q), at 7 TeV for the minimum-bias (MB) events using the unlike-charge particle (UCP) pairs reference sample k_T - intervals&#8758; (a) 100 &lt; k_T &le; 200&nbsp;MeV, (b) 200 &lt; k_T &le; 300&nbsp;MeV, (c) 300 &lt; k_T &le; 400&nbsp;MeV, (d) 400 &lt; k_T &le; 500&nbsp;MeV, (e) 500 &lt; k_T &le; 600&nbsp;MeV, (f) 600 &lt; k_T &le; 700&nbsp;MeV, (g) 700 &lt; k_T &le; 1000&nbsp;MeV, and (h) 1000 &lt; k_T &le; 1500&nbsp;MeV. The error bars represent the statistical uncertainties. The boxes represent the systematic uncertainties, which are the sum in quadrature of a variation of the Coulomb correction, the track reconstruction efficiency and the unfolding matrix.

The single-ratio two-particle correlation functions, C₂^data(Q), at 7 TeV for the minimum-bias (MB) events using the unlike-charge particle (UCP) pairs reference sample k_T - intervals&#8758; (a) 100 &lt; k_T &le; 200&nbsp;MeV, (b) 200 &lt; k_T &le; 300&nbsp;MeV, (c) 300 &lt; k_T &le; 400&nbsp;MeV, (d) 400 &lt; k_T &le; 500&nbsp;MeV, (e) 500 &lt; k_T &le; 600&nbsp;MeV, (f) 600 &lt; k_T &le; 700&nbsp;MeV, (g) 700 &lt; k_T &le; 1000&nbsp;MeV, and (h) 1000 &lt; k_T &le; 1500&nbsp;MeV. The error bars represent the statistical uncertainties. The boxes represent the systematic uncertainties, which are the sum in quadrature of a variation of the Coulomb correction, the track reconstruction efficiency and the unfolding matrix.

The single-ratio two-particle correlation functions, C₂^data(Q), at 7 TeV for the minimum-bias (MB) events using the unlike-charge particle (UCP) pairs reference sample k_T - intervals&#8758; (a) 100 &lt; k_T &le; 200&nbsp;MeV, (b) 200 &lt; k_T &le; 300&nbsp;MeV, (c) 300 &lt; k_T &le; 400&nbsp;MeV, (d) 400 &lt; k_T &le; 500&nbsp;MeV, (e) 500 &lt; k_T &le; 600&nbsp;MeV, (f) 600 &lt; k_T &le; 700&nbsp;MeV, (g) 700 &lt; k_T &le; 1000&nbsp;MeV, and (h) 1000 &lt; k_T &le; 1500&nbsp;MeV. The error bars represent the statistical uncertainties. The boxes represent the systematic uncertainties, which are the sum in quadrature of a variation of the Coulomb correction, the track reconstruction efficiency and the unfolding matrix.

The single-ratio two-particle correlation functions, C₂^data(Q), at 7 TeV for the minimum-bias (MB) events using the unlike-charge particle (UCP) pairs reference sample k_T - intervals&#8758; (a) 100 &lt; k_T &le; 200&nbsp;MeV, (b) 200 &lt; k_T &le; 300&nbsp;MeV, (c) 300 &lt; k_T &le; 400&nbsp;MeV, (d) 400 &lt; k_T &le; 500&nbsp;MeV, (e) 500 &lt; k_T &le; 600&nbsp;MeV, (f) 600 &lt; k_T &le; 700&nbsp;MeV, (g) 700 &lt; k_T &le; 1000&nbsp;MeV, and (h) 1000 &lt; k_T &le; 1500&nbsp;MeV. The error bars represent the statistical uncertainties. The boxes represent the systematic uncertainties, which are the sum in quadrature of a variation of the Coulomb correction, the track reconstruction efficiency and the unfolding matrix.

The single-ratio two-particle correlation functions, C₂^data(Q), at 7 TeV for the minimum-bias (MB) events using the unlike-charge particle (UCP) pairs reference sample k_T - intervals&#8758; (a) 100 &lt; k_T &le; 200&nbsp;MeV, (b) 200 &lt; k_T &le; 300&nbsp;MeV, (c) 300 &lt; k_T &le; 400&nbsp;MeV, (d) 400 &lt; k_T &le; 500&nbsp;MeV, (e) 500 &lt; k_T &le; 600&nbsp;MeV, (f) 600 &lt; k_T &le; 700&nbsp;MeV, (g) 700 &lt; k_T &le; 1000&nbsp;MeV, and (h) 1000 &lt; k_T &le; 1500&nbsp;MeV. The error bars represent the statistical uncertainties. The boxes represent the systematic uncertainties, which are the sum in quadrature of a variation of the Coulomb correction, the track reconstruction efficiency and the unfolding matrix.

The single-ratio two-particle correlation functions, C₂^data(Q), at 7 TeV for the minimum-bias (MB) events using the unlike-charge particle (UCP) pairs reference sample k_T - intervals&#8758; (a) 100 &lt; k_T &le; 200&nbsp;MeV, (b) 200 &lt; k_T &le; 300&nbsp;MeV, (c) 300 &lt; k_T &le; 400&nbsp;MeV, (d) 400 &lt; k_T &le; 500&nbsp;MeV, (e) 500 &lt; k_T &le; 600&nbsp;MeV, (f) 600 &lt; k_T &le; 700&nbsp;MeV, (g) 700 &lt; k_T &le; 1000&nbsp;MeV, and (h) 1000 &lt; k_T &le; 1500&nbsp;MeV. The error bars represent the statistical uncertainties. The boxes represent the systematic uncertainties, which are the sum in quadrature of a variation of the Coulomb correction, the track reconstruction efficiency and the unfolding matrix.

The single-ratio two-particle correlation functions, C₂^data(Q), at 7 TeV for the minimum-bias (MB) events using the unlike-charge particle (UCP) pairs reference sample k_T - intervals&#8758; (a) 100 &lt; k_T &le; 200&nbsp;MeV, (b) 200 &lt; k_T &le; 300&nbsp;MeV, (c) 300 &lt; k_T &le; 400&nbsp;MeV, (d) 400 &lt; k_T &le; 500&nbsp;MeV, (e) 500 &lt; k_T &le; 600&nbsp;MeV, (f) 600 &lt; k_T &le; 700&nbsp;MeV, (g) 700 &lt; k_T &le; 1000&nbsp;MeV, and (h) 1000 &lt; k_T &le; 1500&nbsp;MeV. The error bars represent the statistical uncertainties. The boxes represent the systematic uncertainties, which are the sum in quadrature of a variation of the Coulomb correction, the track reconstruction efficiency and the unfolding matrix.

The single-ratio two-particle correlation functions, C₂^data(Q), at 7 TeV for the minimum-bias (MB) events using the unlike-charge particle (UCP) pairs reference sample k_T - intervals&#8758; (a) 100 &lt; k_T &le; 200&nbsp;MeV, (b) 200 &lt; k_T &le; 300&nbsp;MeV, (c) 300 &lt; k_T &le; 400&nbsp;MeV, (d) 400 &lt; k_T &le; 500&nbsp;MeV, (e) 500 &lt; k_T &le; 600&nbsp;MeV, (f) 600 &lt; k_T &le; 700&nbsp;MeV, (g) 700 &lt; k_T &le; 1000&nbsp;MeV, and (h) 1000 &lt; k_T &le; 1500&nbsp;MeV. The error bars represent the statistical uncertainties. The boxes represent the systematic uncertainties, which are the sum in quadrature of a variation of the Coulomb correction, the track reconstruction efficiency and the unfolding matrix.

The correlation strength, &lambda;, and source radius, R, of the exponential fits to the two-particle double-ratio correlation functions, R₂(Q), in dependence on the multiplicity, m_ch, intervals for the minimum-bias (MB) and the high-multiplicity track (HMT) events for p_T &gt; 100&nbsp;MeV at &radic;s = 13&nbsp;TeV. Statistical uncertainties for &radic;&chi;²/ndf&gt;1 are corrected by the &radic;&chi;²/ndf. The total uncertainties are shown.

The correlation strength, &lambda;, and source radius, R, of the exponential fits to the two-particle double-ratio correlation functions, R₂(Q), in dependence on the multiplicity, m_ch, intervals for the minimum-bias (MB) and the high-multiplicity track (HMT) events for p_T &gt; 500&nbsp;MeV at &radic;s = 13&nbsp;TeV. Statistical uncertainties for &radic;&chi;²/ndf&gt;1 are corrected by the &radic;&chi;²/ndf. The total uncertainties are shown.

The correlation strength, &lambda;, and source radius, R, of the exponential fits to the two-particle double-ratio correlation functions, R₂(Q), in dependence on the pair transverse momentum, k_T, intervals for the minimum-bias (MB) and the high-multiplicity track (HMT) events for p_T &gt; 100&nbsp;MeV at &radic;s = 13&nbsp;TeV. Statistical uncertainties for &radic;&chi;²/ndf&gt;1 are corrected by the &radic;&chi;²/ndf. The total uncertainties are shown.

The correlation strength, &lambda;, and source radius, R, of the exponential fits to the two-particle double-ratio correlation functions, R₂(Q), in dependence on the pair transverse momentum, k_T, intervals for the minimum-bias (MB) and the high-multiplicity track (HMT) events for p_T &gt; 500&nbsp;MeV at &radic;s = 13&nbsp;TeV. Statistical uncertainties for &radic;&chi;²/ndf&gt;1 are corrected by the &radic;&chi;²/ndf. The total uncertainties are shown.

$p_{\mathrm T}$-differential yield of $\frac{\mathrm{K^{*0}} + \overline{\mathrm{K^{*0}}}}{2}$ in p-Pb collisions at $\sqrt{s_{\mathrm{NN}}}~=~$5.02 TeV ($-0.9 < y < -0.6$).

$p_{\mathrm T}$-differential yield of $\frac{\mathrm{K^{*0}} + \overline{\mathrm{K^{*0}}}}{2}$ in p-Pb collisions at $\sqrt{s_{\mathrm{NN}}}~=~$5.02 TeV ($-1.2 < y < -0.9$).

$p_{\mathrm T}$-differential yield of $\phi$ in p-Pb collisions at $\sqrt{s_{\mathrm{NN}}}~=~$5.02 TeV ($0.0 < y < 0.3$).

$p_{\mathrm T}$-differential yield of $\phi$ in p-Pb collisions at $\sqrt{s_{\mathrm{NN}}}~=~$5.02 TeV ($-0.3 < y < 0.0$).

$p_{\mathrm T}$-differential yield of $\phi$ in p-Pb collisions at $\sqrt{s_{\mathrm{NN}}}~=~$5.02 TeV ($-0.6 < y < -0.3$).

$p_{\mathrm T}$-differential yield of $\phi$ in p-Pb collisions at $\sqrt{s_{\mathrm{NN}}}~=~$5.02 TeV ($-0.9 < y < -0.6$).

$p_{\mathrm T}$-differential yield of $\phi$ in p-Pb collisions at $\sqrt{s_{\mathrm{NN}}}~=~$5.02 TeV ($-1.2 < y < -0.9$).

d$N$/d$y$ of $\frac{\mathrm{K^{*0}} + \overline{\mathrm{K^{*0}}}}{2}$ in p-Pb collisions at $\sqrt{s_{\mathrm{NN}}}~=~$5.02 TeV (V0A multiplicity classes).

d$N$/d$y$ of $\phi$ in p-Pb collisions at $\sqrt{s_{\mathrm{NN}}}~=~$5.02 TeV (V0A multiplicity classes).

$\langle p_\mathrm{T}\rangle$ of $\frac{\mathrm{K^{*0}} + \overline{\mathrm{K^{*0}}}}{2}$ in p-Pb collisions at $\sqrt{s_{\mathrm{NN}}}~=~$5.02 TeV (V0A multiplicity classes).

$\langle p_\mathrm{T}\rangle$ of $\phi$ in p-Pb collisions at $\sqrt{s_{\mathrm{NN}}}~=~$5.02 TeV (V0A multiplicity classes).

(d$N$/$\mathrm{d}y$)/(d$N$/$\mathrm{d}y)_{y=0}$ of $\frac{\mathrm{K^{*0}} + \overline{\mathrm{K^{*0}}}}{2}$ in p-Pb collisions at $\sqrt{s_{\mathrm{NN}}}~=~$5.02 TeV (0-100 %).

(d$N$/$\mathrm{d}y$)/(d$N$/$\mathrm{d}y)_{y=0}$ of $\phi$ in p-Pb collisions at $\sqrt{s_{\mathrm{NN}}}~=~$5.02 TeV (0-100 %).

$\langle p_\mathrm{T}\rangle_\mathrm{y}$/$\langle p_\mathrm{T}\rangle_\mathrm{y=0}$ of $\frac{\mathrm{K^{*0}} + \overline{\mathrm{K^{*0}}}}{2}$ in p-Pb collisions at $\sqrt{s_{\mathrm{NN}}}~=~$5.02 TeV (0-100 %).

$\langle p_\mathrm{T}\rangle_\mathrm{y}$/$\langle p_\mathrm{T}\rangle_\mathrm{y=0}$ of $\phi$ in p-Pb collisions at $\sqrt{s_{\mathrm{NN}}}~=~$5.02 TeV (0-100%).

$y_{asym}$ of $\frac{\mathrm{K^{*0}} + \overline{\mathrm{K^{*0}}}}{2}$ in p-Pb collisions at $\sqrt{s_{\mathrm{NN}}}~=~$5.02 TeV (V0A multiplicity classes).

$y_{asym}$ of $\phi$ in p-Pb collisions at $\sqrt{s_{\mathrm{NN}}}~=~$5.02 TeV (V0A multiplicity classes).

$Q_{CP}$ of $\frac{\mathrm{K^{*0}} + \overline{\mathrm{K^{*0}}}}{2}$ in p-Pb collisions at $\sqrt{s_{\mathrm{NN}}}~=~$5.02 TeV (0-10 %/40-100 %).

$Q_{CP}$ of $\phi$ in p-Pb collisions at $\sqrt{s_{\mathrm{NN}}}~=~$5.02 TeV (0-10 %/40-100%).

$Q_{CP}$ of $\frac{\mathrm{K^{*0}} + \overline{\mathrm{K^{*0}}}}{2}$ in p-Pb collisions at $\sqrt{s_{\mathrm{NN}}}~=~$5.02 TeV (10-40 %/40-100 %).

$Q_{CP}$ of $\phi$ in p-Pb collisions at $\sqrt{s_{\mathrm{NN}}}~=~$5.02 TeV (10-40 %/40-100%).

Fig. 6b: Number density $N_{\rm ch}$ (left) and $\Sigma p_{\rm T}$ (right) distributions as a function of $p_{\rm T}^{\rm trig}$ in Away and Toward regions after the subtraction of Number density $N_{\rm ch}$ and $\Sigma p_{\rm T}$ distributions in the transverse region for p-Pb collisions for $p_{\rm T} >$ 0.5 GeV/$c$. The shaded areas and the error bars around the data points represent the systematic and statistical uncertainties, respectively.

Fig. 7a: $<p_{\rm T}>$ as a function of $p_{\rm T}^{\rm trig}$ in Toward (left) and Away (right) regions after the subtraction of transverse region for pp collisions for $p_{\rm T} >$ 0.5 GeV/$c$. The shaded areas and the error bars around the data points represent the systematic and statistical uncertainties, respectively.

Fig. 7b: $<p_{\rm T}>$ as a function of $p_{\rm T}^{\rm trig}$ in Toward (left) and Away (right) regions after the subtraction of transverse region for p-Pb collisions for $p_{\rm T} >$ 0.5 GeV/$c$. The shaded areas and the error bars around the data points represent the systematic and statistical uncertainties, respectively.

Fig. A1: Number density $N_{\rm ch}$ (left) and $\Sigma p_{\rm T}$ (right) distributions as a function of $p_{\rm T}^{\rm trig}$ in Transverse, Away, and Toward regions for $p_{\rm T} >$ 0.15 GeV/$c$. The shaded areas and the error bars around the data points represent the systematic and statistical uncertainties, respectively.

Fig. A2: Number density $N_{\rm ch}$ (left) and $\Sigma p_{\rm T}$ (right) distributions as a function of $p_{\rm T}^{\rm trig}$ in Transverse, Away, and Toward regions for $p_{\rm T} >$ 1.0 GeV/$c$. The shaded areas and the error bars around the data points represent the systematic and statistical uncertainties, respectively.

Fig. A3: Number density $N_{\rm ch}$ (left) and $\Sigma p_{\rm T}$ (right) distributions as a function of $p_{\rm T}^{\rm trig}$ in Transverse, Away, and Toward regions for $p_{\rm T} >$ 0.15 GeV/$c$. The shaded areas and the error bars around the data points represent the systematic and statistical uncertainties, respectively.

Fig. A4: Number density $N_{\rm ch}$ (left) and $\Sigma p_{\rm T}$ (right) distributions as a function of $p_{\rm T}^{\rm trig}$ in Transverse, Away, and Toward regions for $p_{\rm T} >$ 1.0 GeV/$c$. The shaded areas and the error bars around the data points represent the systematic and statistical uncertainties, respectively.