Fiducial differential cross section of the electroweak $W^\pm W^\pm jj$ production as a function of $N_{\mathrm{gap}\,\mathrm{jets}}$. The correlation of uncertainties of the measured cross section across bins is presented in Table 14.

Fiducial differential cross section of the electroweak $W^\pm W^\pm jj$ production as a function of $\xi_{\mathrm{j}3}$. The correlation of uncertainties of the measured cross section across bins is presented in Table 15.

Fiducial differential cross section of the inclusive $W^\pm W^\pm jj$ production as a function of $m_{\ell\ell}$. The correlation of uncertainties of the measured cross section across bins is presented in Table 16.

Fiducial differential cross section of the inclusive $W^\pm W^\pm jj$ production as a function of $m_{\mathrm{T}}$. The correlation of uncertainties of the measured cross section across bins is presented in Table 17.

Fiducial differential cross section of the inclusive $W^\pm W^\pm jj$ production as a function of $m_{\mathrm{jj}}$. The correlation of uncertainties of the measured cross section across bins is presented in Table 18.

Fiducial differential cross section of the inclusive $W^\pm W^\pm jj$ production as a function of $N_{\mathrm{gap}\,\mathrm{jets}}$. The correlation of uncertainties of the measured cross section across bins is presented in Table 19.

Fiducial differential cross section of the inclusive $W^\pm W^\pm jj$ production as a function of $\xi_{\mathrm{j}3}$. The correlation of uncertainties of the measured cross section across bins is presented in Table 20.

Observed correlations between the bins of the LH-unfolded cross section of the electroweak $W^\pm W^\pm jj$ production as a function of $m_{\ell\ell}$. See Table 1.

Observed correlations between the bins of the LH-unfolded cross section of the electroweak $W^\pm W^\pm jj$ production as a function of $m_{\mathrm{T}}$. See Table 2.

Observed correlations between the bins of the LH-unfolded cross section of the electroweak $W^\pm W^\pm jj$ production as a function of $m_{\mathrm{jj}}$. See Table 3.

Observed correlations between the bins of the LH-unfolded cross section of the electroweak $W^\pm W^\pm jj$ production as a function of $N_{\mathrm{gap}\,\mathrm{jets}}$. See Table 4.

Observed correlations between the bins of the LH-unfolded cross section of the electroweak $W^\pm W^\pm jj$ production as a function of $\xi_{\mathrm{j}3}$. See Table 5.

Observed correlations between the bins of the LH-unfolded cross section of the inclusive $W^\pm W^\pm jj$ production as a function of $m_{\ell\ell}$. See Table 6.

Observed correlations between the bins of the LH-unfolded cross section of the inclusive $W^\pm W^\pm jj$ production as a function of $m_{\mathrm{T}}$. See Table 7.

Observed correlations between the bins of the LH-unfolded cross section of the inclusive $W^\pm W^\pm jj$ production as a function of $m_{\mathrm{jj}}$. See Table 8.

Observed correlations between the bins of the LH-unfolded cross section of the inclusive $W^\pm W^\pm jj$ production as a function of $N_{\mathrm{gap}\,\mathrm{jets}}$. See Table 9.

Observed correlations between the bins of the LH-unfolded cross section of the inclusive $W^\pm W^\pm jj$ production as a function of $\xi_{\mathrm{j}3}$. See Table 10.

Evolution of the one-dimensional expected and observed limits at 95% CL on the parameters corresponding to the quartic operators with label M0 as a function of the cut-off scale. The unitarity bounds as a function of the cut-off scale are defined for one non-zero Wilson coefficient.

Evolution of the one-dimensional expected and observed limits at 95% CL on the parameters corresponding to the quartic operators with label M1 as a function of the cut-off scale. The unitarity bounds as a function of the cut-off scale are defined for one non-zero Wilson coefficient.

Evolution of the one-dimensional expected and observed limits at 95% CL on the parameters corresponding to the quartic operators with label M7 as a function of the cut-off scale. The unitarity bounds as a function of the cut-off scale are defined for one non-zero Wilson coefficient. The limits on M7 were obtained without taking into account the SM-EFT interference for the EW W$^\pm$Zjj final state.

Evolution of the one-dimensional expected and observed limits at 95% CL on the parameters corresponding to the quartic operators with label S02 as a function of the cut-off scale. The unitarity bounds as a function of the cut-off scale are defined for one non-zero Wilson coefficient.

Evolution of the one-dimensional expected and observed limits at 95% CL on the parameters corresponding to the quartic operators with label S1 as a function of the cut-off scale. The unitarity bounds as a function of the cut-off scale are defined for one non-zero Wilson coefficient.

Evolution of the one-dimensional expected and observed limits at 95% CL on the parameters corresponding to the quartic operators with label T0 as a function of the cut-off scale. The unitarity bounds as a function of the cut-off scale are defined for one non-zero Wilson coefficient.

Evolution of the one-dimensional expected and observed limits at 95% CL on the parameters corresponding to the quartic operators with label T1 as a function of the cut-off scale. The unitarity bounds as a function of the cut-off scale are defined for one non-zero Wilson coefficient.

Evolution of the one-dimensional expected and observed limits at 95% CL on the parameters corresponding to the quartic operators with label T2 as a function of the cut-off scale. The unitarity bounds as a function of the cut-off scale are defined for one non-zero Wilson coefficient.

Expected and observed exclusion limits at 95% CL for $\sigma_{\mathrm{VBF}}(\mathrm{H}^{\pm\pm}_5) \times \mathcal{B}(\mathrm{H}^{\pm\pm}_5 \rightarrow W^{\pm}W^{\pm})$ as a function of the doubly-charged Higgs mass.

Expected and observed exclusion limits at 95% CL for $\sin\theta_{\mathrm{H}}$ as a function of the doubly-charged Higgs mass in the Georgi-Machacek model.

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.

$cos\theta^{*}$ distribution of $\phi$ meson w.r.t. Production Plane in transverse momentum range 0.5-0.8 GeV/$c$ measured in Pb-Pb collisions at $\sqrt{s_{NN}}$ = 2.76 TeV.

$cos\theta^{*}$ distribution of $\rm{K}^{*0}$ (average of particle and anti-particle) meson w.r.t. Random Plane in transverse momentum range 0.8-1.2 GeV/$c$ measured in Pb-Pb collisions at $\sqrt{s_{NN}}$ = 2.76 TeV.

$cos\theta^{*}$ distribution of $\phi$ meson w.r.t. Random Plane in transverse momentum range 0.5-0.7 GeV/$c$ measured in Pb-Pb collisions at $\sqrt{s_{NN}}$ = 2.76 TeV.

$cos\theta^{*}$ distribution of K$^{0}_{S}$ meson w.r.t. Event Plane in transverse momentum range 0.6-0.8 GeV/$c$ measured in Pb-Pb collisions at $\sqrt{s_{NN}}$ = 2.76 TeV.

$cos\theta^{*}$ distribution of K$^{0}_{S}$ meson w.r.t. Production Plane in transverse momentum range 0.6-0.8 GeV/$c$ measured in Pb-Pb collisions at $\sqrt{s_{NN}}$ = 2.76 TeV.

$cos\theta^{*}$ distribution of $\rm{K}^{*0}$ (average of particle and anti-particle) meson w.r.t. Production Plane in transverse momentum range 0-0.6 GeV/$c$ measured in pp collisions at $\sqrt{s}$ = 13 TeV.

$cos\theta^{*}$ distribution of $\phi$ meson w.r.t. Production Plane in transverse momentum range 0.5-0.8 GeV/$c$ measured in pp collisions at $\sqrt{s}$ = 13 TeV.

$\rho_{00}$ as a function of $p_{\rm T}$ w.r.t. Event Plane for $\rm{K}^{*0}$ (average of particle and anti-particle) meson measured in Pb-Pb collisions at $\sqrt{s_{NN}}$ = 2.76 TeV.

$\rho_{00}$ as a function of $p_{\rm T}$ w.r.t. Event Plane for $\rm{K}^{0}_{S}$ meson measured in Pb-Pb collisions at $\sqrt{s_{NN}}$ = 2.76 TeV.

$\rho_{00}$ as a function of $p_{\rm T}$ w.r.t. Event Plane for $\phi$ meson measured in Pb-Pb collisions at $\sqrt{s_{NN}}$ = 2.76 TeV.

$\rho_{00}$ as a function of $p_{\rm T}$ w.r.t. Production Plane for $\rm{K}^{*0}$ (average of particle and anti-particle) meson measured in Pb-Pb collisions at $\sqrt{s_{NN}}$ = 2.76 TeV.

$\rho_{00}$ as a function of $p_{\rm T}$ w.r.t. Production Plane for $\rm{K}^{*0}$ (average of particle and anti-particle) meson measured in pp collisions at $\sqrt{s}$ = 13 TeV.

$\rho_{00}$ as a function of $p_{\rm T}$ w.r.t. Production Plane for $\rm{K}^{0}_{S}$ meson measured in Pb-Pb collisions at $\sqrt{s_{NN}}$ = 2.76 TeV.

$\rho_{00}$ as a function of $p_{\rm T}$ w.r.t. Production Plane for $\phi$ meson measured in Pb-Pb collisions at $\sqrt{s_{NN}}$ = 2.76 TeV.

$\rho_{00}$ as a function of $p_{\rm T}$ w.r.t. Production Plane for $\phi$ meson measured in pp collisions at $\sqrt{s}$ = 13 TeV.

$\rho_{00}$ as a function of $p_{\rm T}$ w.r.t. Random Plane for $\rm{K}^{*0}$ (average of particle and anti-particle) meson measured in Pb-Pb collisions at $\sqrt{s_{NN}}$ = 2.76 TeV.

$\rho_{00}$ as a function of $p_{\rm T}$ w.r.t. Random Plane for $\phi$ meson measured in Pb-Pb collisions at $\sqrt{s_{NN}}$ = 2.76 TeV.

$\rho_{00}$ as a function of $\langle N_{\rm{part}} \rangle$ w.r.t. Production Plane in transverse momentum range 0.4-1.2 GeV/$c$ for $\rm{K}^{*0}$ (average of particle and anti-particle) meson measured in Pb-Pb collisions at $\sqrt{s_{NN}}$ = 2.76 TeV.

$\rho_{00}$ as a function of $\langle N_{\rm{part}} \rangle$ w.r.t. Production Plane in transverse momentum range 3.0-5.0 GeV/$c$ for $\rm{K}^{*0}$ (average of particle and anti-particle) meson measured in Pb-Pb collisions at $\sqrt{s_{NN}}$ = 2.76 TeV.

$\rho_{00}$ as a function of $\langle N_{\rm{part}} \rangle$ w.r.t. Production Plane in transverse momentum range 0.5-0.8 GeV/$c$ for $\phi$ meson measured in Pb-Pb collisions at $\sqrt{s_{NN}}$ = 2.76 TeV.

$\rho_{00}$ as a function of $\langle N_{\rm{part}} \rangle$ w.r.t. Production Plane in transverse momentum range 3.0-5.0 GeV/$c$ for $\phi$ meson measured in Pb-Pb collisions at $\sqrt{s_{NN}}$ = 2.76 TeV.

$\rho_{00}$ as a function of $\langle N_{\rm{part}} \rangle$ w.r.t. Event Plane in transverse momentum range 0.8-1.2 GeV/$c$ for $\rm{K}^{*0}$ (average of particle and anti-particle) meson measured in Pb-Pb collisions at $\sqrt{s_{NN}}$ = 2.76 TeV.

$\rho_{00}$ as a function of $\langle N_{\rm{part}} \rangle$ w.r.t. Event Plane in transverse momentum range 3.0-5.0 GeV/$c$ for $\rm{K}^{*0}$ (average of particle and anti-particle) meson measured in Pb-Pb collisions at $\sqrt{s_{NN}}$ = 2.76 TeV.

$\rho_{00}$ as a function of $\langle N_{\rm{part}} \rangle$ w.r.t. Event Plane in transverse momentum range 0.5-0.7 GeV/$c$ for $\phi$ meson measured in Pb-Pb collisions at $\sqrt{s_{NN}}$ = 2.76 TeV.

$\rho_{00}$ as a function of $\langle N_{\rm{part}} \rangle$ w.r.t. Event Plane in transverse momentum range 3.0-5.0 GeV/$c$ for $\phi$ meson measured in Pb-Pb collisions at $\sqrt{s_{NN}}$ = 2.76 TeV.

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

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.

Fig. 1 Left, data for jet radius R=0.2. Unfolded pp full jet cross-section at $\sqrt{s}$ = 5.02 TeV for R = 0.1 − 0.6. No leading track requirement is imposed.

Fig. 1 Left, data for jet radius R=0.3. Unfolded pp full jet cross-section at $\sqrt{s}$ = 5.02 TeV for R = 0.1 − 0.6. No leading track requirement is imposed.

Fig. 1 Left, data for jet radius R=0.3. Unfolded pp full jet cross-section at $\sqrt{s}$ = 5.02 TeV for R = 0.1 − 0.6. No leading track requirement is imposed.

Fig. 1 Left, data for jet radius R=0.4. Unfolded pp full jet cross-section at $\sqrt{s}$ = 5.02 TeV for R = 0.1 − 0.6. No leading track requirement is imposed.

Fig. 1 Left, data for jet radius R=0.4. Unfolded pp full jet cross-section at $\sqrt{s}$ = 5.02 TeV for R = 0.1 − 0.6. No leading track requirement is imposed.

Fig. 1 Left, data for jet radius R=0.5. Unfolded pp full jet cross-section at $\sqrt{s}$ = 5.02 TeV for R = 0.1 − 0.6. No leading track requirement is imposed.

Fig. 1 Left, data for jet radius R=0.5. Unfolded pp full jet cross-section at $\sqrt{s}$ = 5.02 TeV for R = 0.1 − 0.6. No leading track requirement is imposed.

Fig. 1 Left, data for jet radius R=0.6. Unfolded pp full jet cross-section at $\sqrt{s}$ = 5.02 TeV for R = 0.1 − 0.6. No leading track requirement is imposed.

Fig. 1 Left, data for jet radius R=0.6. Unfolded pp full jet cross-section at $\sqrt{s}$ = 5.02 TeV for R = 0.1 − 0.6. No leading track requirement is imposed.

Fig. 2, values for jet radius R=0.1. Non-perturbative correction factor applied to parton-level NLO+NLL predictions, obtained from PYTHIA 8 tune A14 as the ratio of the inclusive jet spectrum at hadron-level with MPI compared to parton-level without MPI.

Fig. 2, values for jet radius R=0.1. Non-perturbative correction factor applied to parton-level NLO+NLL predictions, obtained from PYTHIA 8 tune A14 as the ratio of the inclusive jet spectrum at hadron-level with MPI compared to parton-level without MPI.

Fig. 2, values for jet radius R=0.2. Non-perturbative correction factor applied to parton-level NLO+NLL predictions, obtained from PYTHIA 8 tune A14 as the ratio of the inclusive jet spectrum at hadron-level with MPI compared to parton-level without MPI.

Fig. 2, values for jet radius R=0.2. Non-perturbative correction factor applied to parton-level NLO+NLL predictions, obtained from PYTHIA 8 tune A14 as the ratio of the inclusive jet spectrum at hadron-level with MPI compared to parton-level without MPI.

Fig. 2, values for jet radius R=0.3. Non-perturbative correction factor applied to parton-level NLO+NLL predictions, obtained from PYTHIA 8 tune A14 as the ratio of the inclusive jet spectrum at hadron-level with MPI compared to parton-level without MPI.

Fig. 2, values for jet radius R=0.3. Non-perturbative correction factor applied to parton-level NLO+NLL predictions, obtained from PYTHIA 8 tune A14 as the ratio of the inclusive jet spectrum at hadron-level with MPI compared to parton-level without MPI.

Fig. 2, values for jet radius R=0.4. Non-perturbative correction factor applied to parton-level NLO+NLL predictions, obtained from PYTHIA 8 tune A14 as the ratio of the inclusive jet spectrum at hadron-level with MPI compared to parton-level without MPI.

Fig. 2, values for jet radius R=0.4. Non-perturbative correction factor applied to parton-level NLO+NLL predictions, obtained from PYTHIA 8 tune A14 as the ratio of the inclusive jet spectrum at hadron-level with MPI compared to parton-level without MPI.

Fig. 2, values for jet radius R=0.5. Non-perturbative correction factor applied to parton-level NLO+NLL predictions, obtained from PYTHIA 8 tune A14 as the ratio of the inclusive jet spectrum at hadron-level with MPI compared to parton-level without MPI.

Fig. 2, values for jet radius R=0.5. Non-perturbative correction factor applied to parton-level NLO+NLL predictions, obtained from PYTHIA 8 tune A14 as the ratio of the inclusive jet spectrum at hadron-level with MPI compared to parton-level without MPI.

Fig. 2, values for jet radius R=0.6. Non-perturbative correction factor applied to parton-level NLO+NLL predictions, obtained from PYTHIA 8 tune A14 as the ratio of the inclusive jet spectrum at hadron-level with MPI compared to parton-level without MPI.

Fig. 2, values for jet radius R=0.6. Non-perturbative correction factor applied to parton-level NLO+NLL predictions, obtained from PYTHIA 8 tune A14 as the ratio of the inclusive jet spectrum at hadron-level with MPI compared to parton-level without MPI.

Fig. 3 Bottom, data for jet radius ratio R=0.1/R=0.2. Unfolded pp jet cross-section ratios for various R. With R = 0.1 in the numerator.

Fig. 3 Bottom, data for jet radius ratio R=0.1/R=0.2. Unfolded pp jet cross-section ratios for various R. With R = 0.1 in the numerator.

Fig. 3 Bottom, data for jet radius ratio R=0.1/R=0.3. Unfolded pp jet cross-section ratios for various R. With R = 0.1 in the numerator.

Fig. 3 Bottom, data for jet radius ratio R=0.1/R=0.3. Unfolded pp jet cross-section ratios for various R. With R = 0.1 in the numerator.

Fig. 3 Bottom, data for jet radius ratio R=0.1/R=0.4. Unfolded pp jet cross-section ratios for various R. With R = 0.1 in the numerator.

Fig. 3 Bottom, data for jet radius ratio R=0.1/R=0.4. Unfolded pp jet cross-section ratios for various R. With R = 0.1 in the numerator.

Fig. 3 Bottom, data for jet radius ratio R=0.1/R=0.5. Unfolded pp jet cross-section ratios for various R. With R = 0.1 in the numerator.

Fig. 3 Bottom, data for jet radius ratio R=0.1/R=0.5. Unfolded pp jet cross-section ratios for various R. With R = 0.1 in the numerator.

Fig. 3 Bottom, data for jet radius ratio R=0.1/R=0.6. Unfolded pp jet cross-section ratios for various R. With R = 0.1 in the numerator.

Fig. 3 Bottom, data for jet radius ratio R=0.1/R=0.6. Unfolded pp jet cross-section ratios for various R. With R = 0.1 in the numerator.

Fig. 3 Top, data for jet radius ratio R=0.2/R=0.3. Unfolded pp jet cross-section ratios for various R. With R = 0.2 in the numerator.

Fig. 3 Top, data for jet radius ratio R=0.2/R=0.3. Unfolded pp jet cross-section ratios for various R. With R = 0.2 in the numerator.

Fig. 3 Top, data for jet radius ratio R=0.2/R=0.4. Unfolded pp jet cross-section ratios for various R. With R = 0.2 in the numerator.

Fig. 3 Top, data for jet radius ratio R=0.2/R=0.4. Unfolded pp jet cross-section ratios for various R. With R = 0.2 in the numerator.

Fig. 3 Top, data for jet radius ratio R=0.2/R=0.5. Unfolded pp jet cross-section ratios for various R. With R = 0.2 in the numerator.

Fig. 3 Top, data for jet radius ratio R=0.2/R=0.5. Unfolded pp jet cross-section ratios for various R. With R = 0.2 in the numerator.

Fig. 3 Top, data for jet radius ratio R=0.2/R=0.6. Unfolded pp jet cross-section ratios for various R. With R = 0.2 in the numerator.

Fig. 3 Top, data for jet radius ratio R=0.2/R=0.6. Unfolded pp jet cross-section ratios for various R. With R = 0.2 in the numerator.

Fig. 4 Left, data for pp. Unfolded pp and Pb-Pb full jet spectra at $\sqrt{s}=5.02$ TeV for $R=0.2$, with 5 GeV leading track requirement. The pp data points correspond to $\frac{\mathrm{d}^{2}\sigma}{\mathrm{d}p_{\mathrm{T,jet}} \mathrm{d}\eta_{\mathrm{jet}}}$.

Fig. 4 Left, data for pp. Unfolded pp and Pb-Pb full jet spectra at $\sqrt{s}=5.02$ TeV for $R=0.2$, with 5 GeV leading track requirement. The pp data points correspond to $\frac{\mathrm{d}^{2}\sigma}{\mathrm{d}p_{\mathrm{T,jet}} \mathrm{d}\eta_{\mathrm{jet}}}$.

Fig. 4 Left, data for PbPb. Unfolded pp and Pb-Pb full jet spectra at $\sqrt{s}=5.02$ TeV for $R=0.2$, with 5 GeV leading track requirement. The pp data points correspond to $\frac{\mathrm{d}^{2}\sigma}{\mathrm{d}p_{\mathrm{T,jet}} \mathrm{d}\eta_{\mathrm{jet}}}$.

Fig. 4 Left, data for PbPb. Unfolded pp and Pb-Pb full jet spectra at $\sqrt{s}=5.02$ TeV for $R=0.2$, with 5 GeV leading track requirement. The pp data points correspond to $\frac{\mathrm{d}^{2}\sigma}{\mathrm{d}p_{\mathrm{T,jet}} \mathrm{d}\eta_{\mathrm{jet}}}$.

Fig. 4 Right, data for pp. Unfolded pp and Pb-Pb full jet spectra at $\sqrt{s}=5.02$ TeV for $R=0.4$, with 7 GeV leading track requirement.

Fig. 4 Right, data for pp. Unfolded pp and Pb-Pb full jet spectra at $\sqrt{s}=5.02$ TeV for $R=0.4$, with 7 GeV leading track requirement.

Fig. 4 Right, data for PbPb. Unfolded pp and Pb-Pb full jet spectra at $\sqrt{s}=5.02$ TeV for $R=0.4$, with 7 GeV leading track requirement.

Fig. 4 Right, data for PbPb. Unfolded pp and Pb-Pb full jet spectra at $\sqrt{s}=5.02$ TeV for $R=0.4$, with 7 GeV leading track requirement.

Fig. 5 Left, data for pp jet cross section with leading track bias ratio 5/0. Ratio of the pp jet cross-section with various leading charged particle requirements

Fig. 5 Left, data for pp jet cross section with leading track bias ratio 5/0. Ratio of the pp jet cross-section with various leading charged particle requirements

Fig. 5 Left, data for pp jet cross section with leading track bias ratio 7/0. Ratio of the pp jet cross-section with various leading charged particle requirements

Fig. 5 Left, data for pp jet cross section with leading track bias ratio 7/0. Ratio of the pp jet cross-section with various leading charged particle requirements

Fig. 5 Left, data for pp jet cross section with leading track bias ratio 7/5. Ratio of the pp jet cross-section with various leading charged particle requirements

Fig. 5 Left, data for pp jet cross section with leading track bias ratio 7/5. Ratio of the pp jet cross-section with various leading charged particle requirements

Fig. 5 Right. Ratio of the R = 0.2 Pb–Pb jet cross-section with a 7 GeV/c leading charged particle requirement compared to a 5 GeV/c leading charged particle requirement.

Fig. 5 Right. Ratio of the R = 0.2 Pb–Pb jet cross-section with a 7 GeV/c leading charged particle requirement compared to a 5 GeV/c leading charged particle requirement.

Fig. 6 Top Left. Jet $R_{\mathrm{AA}}$ at $\sqrt{s}=5.02$ TeV for $R=0.2$. Same leading track requirement of 5 GeV is imposed on the numerator as well as denominator.

Fig. 6 Top Left. Jet $R_{\mathrm{AA}}$ at $\sqrt{s}=5.02$ TeV for $R=0.2$. Same leading track requirement of 5 GeV is imposed on the numerator as well as denominator.

Fig. 6 Top Right. Jet $R_{\mathrm{AA}}$ at $\sqrt{s}=5.02$ TeV for $R=0.4$. Same leading track requirement of 7 GeV is imposed on the numerator as well as denominator.

Fig. 6 Top Right. Jet $R_{\mathrm{AA}}$ at $\sqrt{s}=5.02$ TeV for $R=0.4$. Same leading track requirement of 7 GeV is imposed on the numerator as well as denominator.

Fig. 6 Bottom Left. Jet $R_{\mathrm{AA}}$ at $\sqrt{s}=5.02$ TeV for $R=0.2$. PbPb with leading track requirement of 5 GeV. pp reference without a leading track requirement.

Fig. 6 Bottom Left. Jet $R_{\mathrm{AA}}$ at $\sqrt{s}=5.02$ TeV for $R=0.2$. PbPb with leading track requirement of 5 GeV. pp reference without a leading track requirement.

Fig. 6 Bottom Right. Jet $R_{\mathrm{AA}}$ at $\sqrt{s}=5.02$ TeV for $R=0.4$. PbPb with leading track requirement of 7 GeV. pp reference without a leading track requirement.

Fig. 6 Bottom Right. Jet $R_{\mathrm{AA}}$ at $\sqrt{s}=5.02$ TeV for $R=0.4$. PbPb with leading track requirement of 7 GeV. pp reference without a leading track requirement.

Charged-hadron spectrum in the centrality interval 5-10% for p+Pb, divided by 〈TPPB〉. 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 〈TPPB〉. 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 〈TPPB〉. 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 〈TPPB〉. 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 〈TPPB〉. 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 〈TPPB〉. 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 〈TPPB〉. 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 〈TAA〉. 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 〈TAA〉. 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 〈TAA〉. 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 〈TAA〉. 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 〈TAA〉. 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 〈TAA〉. 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 〈TAA〉. 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 〈TAA〉. 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 〈TAA〉. 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 〈TAA〉. 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 〈TAA〉. 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 〈TAA〉. 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 〈TAA〉. 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 〈TAA〉. 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 〈TAA〉. 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 〈TAA〉. 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.

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.

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 λ 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 >100 MeV at √s=13 TeV in the multiplicity interval 71 ≤ n_ch < 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 >100 MeV at √s=13 TeV in the multiplicity interval 231 ≤ n_ch < 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, λ(m_ch), on rescaled multiplicity, m_ch, obtained from the exponential fit of the R₂(Q) correlation functions for tracks with p_T > 100  MeV and p_T > 500 MeV at √s = 13  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 λ(m_ch) for p_T >100 MeV and p_T >500 MeV, respectively.

The dependence of the correlation strength, λ(m_ch), on rescaled multiplicity, m_ch, obtained from the exponential fit of the R₂(Q) correlation functions for tracks with p_T > 100  MeV and p_T > 500 MeV at √s = 13  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 λ(m_ch) for p_T >100 MeV and p_T >500 MeV, respectively.

The dependence of the correlation strength, λ(m_ch), on rescaled multiplicity, m_ch, obtained from the exponential fit of the R₂(Q) correlation functions for tracks with p_T > 100  MeV and p_T > 500 MeV at √s = 13  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 λ(m_ch) for p_T >100 MeV and p_T >500 MeV, respectively.

The dependence of the correlation strength, λ(m_ch), on rescaled multiplicity, m_ch, obtained from the exponential fit of the R₂(Q) correlation functions for tracks with p_T > 100  MeV and p_T > 500 MeV at √s = 13  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 λ(m_ch) for p_T >100 MeV and p_T >500 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 ∛m_ch < 1.2 for p_T >100 MeV and p_T >500 MeV, respectively. The black and blue dotted curves are extensions of the black and blue solid curves beyond ∛m_ch > 1.2, respectively. The black and brown dashed curves represent the saturation value of R(m_ch) for ∛m_ch > 1.45 with p_T >100 MeV and for ∛m_ch > 1.6 with p_T >500 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 ∛m_ch < 1.2 for p_T >100 MeV and p_T >500 MeV, respectively. The black and blue dotted curves are extensions of the black and blue solid curves beyond ∛m_ch > 1.2, respectively. The black and brown dashed curves represent the saturation value of R(m_ch) for ∛m_ch > 1.45 with p_T >100 MeV and for ∛m_ch > 1.6 with p_T >500 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 ∛m_ch < 1.2 for p_T >100 MeV and p_T >500 MeV, respectively. The black and blue dotted curves are extensions of the black and blue solid curves beyond ∛m_ch > 1.2, respectively. The black and brown dashed curves represent the saturation value of R(m_ch) for ∛m_ch > 1.45 with p_T >100 MeV and for ∛m_ch > 1.6 with p_T >500 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 ∛m_ch < 1.2 for p_T >100 MeV and p_T >500 MeV, respectively. The black and blue dotted curves are extensions of the black and blue solid curves beyond ∛m_ch > 1.2, respectively. The black and brown dashed curves represent the saturation value of R(m_ch) for ∛m_ch > 1.45 with p_T >100 MeV and for ∛m_ch > 1.6 with p_T >500 MeV, respectively.

The dependence of the 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 ∛m_ch < 1.2 for p_T >100 MeV and p_T >500 MeV, respectively. The black and blue dotted curves are extensions of the black and blue solid curves beyond ∛m_ch > 1.2, respectively. The black and brown dashed curves represent the saturation value of R(m_ch) for ∛m_ch > 1.45 with p_T >100 MeV and for ∛m_ch > 1.6 with p_T >500 MeV, respectively

The dependence of the 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 ∛m_ch < 1.2 for p_T >100 MeV and p_T >500 MeV, respectively. The black and blue dotted curves are extensions of the black and blue solid curves beyond ∛m_ch > 1.2, respectively. The black and brown dashed curves represent the saturation value of R(m_ch) for ∛m_ch > 1.45 with p_T >100 MeV and for ∛m_ch > 1.6 with p_T >500 MeV, respectively

The dependence of the 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 ∛m_ch < 1.2 for p_T >100 MeV and p_T >500 MeV, respectively. The black and blue dotted curves are extensions of the black and blue solid curves beyond ∛m_ch > 1.2, respectively. The black and brown dashed curves represent the saturation value of R(m_ch) for ∛m_ch > 1.45 with p_T >100 MeV and for ∛m_ch > 1.6 with p_T >500 MeV, respectively

The dependence of the 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 ∛m_ch < 1.2 for p_T >100 MeV and p_T >500 MeV, respectively. The black and blue dotted curves are extensions of the black and blue solid curves beyond ∛m_ch > 1.2, respectively. The black and brown dashed curves represent the saturation value of R(m_ch) for ∛m_ch > 1.45 with p_T >100 MeV and for ∛m_ch > 1.6 with p_T >500 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 TeV (black circles) and 7 TeV (open blue circles), and the ratio of C₂^7 TeV (Q) to C₂^13 TeV (Q) (bottom panel). Comparison of C₂^data (Q) for representative multiplicity region 3.09 < m_ch ≤ 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  TeV (black circles) and 7  TeV (open blue circles), and the ratio of C₂^7 TeV (Q) to C₂^13 TeV (Q) (bottom panel). Comparison of C₂^data (Q) for representative k_T region 400 < k_T ≤500  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, λ(k_T), obtained from the exponential fit to the R₂(Q) correlation functions for events with multiplicity n_ch ≥ 2 and transfer momentum of tracks with p_T >100 MeV and p_T >500 MeV at √s=13 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 λ(k_T).

The k_T dependence of the correlation strength, λ(k_T), obtained from the exponential fit to the R₂(Q) correlation functions for events with multiplicity n_ch ≥ 2 and transfer momentum of tracks with p_T >100 MeV and p_T >500 MeV at √s=13 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 λ(k_T).

The k_T dependence of the correlation strength, λ(k_T), obtained from the exponential fit to the R₂(Q) correlation functions for events with multiplicity n_ch ≥ 2 and transfer momentum of tracks with p_T >100 MeV and p_T >500 MeV at √s=13 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 λ(k_T).

The k_T dependence of the correlation strength, λ(k_T), obtained from the exponential fit to the R₂(Q) correlation functions for events with multiplicity n_ch ≥ 2 and transfer momentum of tracks with p_T >100 MeV and p_T >500 MeV at √s=13 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 λ(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 ≥ 2 and transfer momentum of tracks with p_T >100 MeV and p_T >500 MeV at √s=13 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 ≥ 2 and transfer momentum of tracks with p_T >100 MeV and p_T >500 MeV at √s=13 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 ≥ 2 and transfer momentum of tracks with p_T >100 MeV and p_T >500 MeV at √s=13 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 ≥ 2 and transfer momentum of tracks with p_T >100 MeV and p_T >500 MeV at √s=13 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 > 100 MeV for the correlation strength, λ, obtained from the exponential fit to the R₂(Q) correlation functions using the MB sample for m_ch ≤ 3.08 and the HMT sample for m_ch > 3.08.

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

The parameter λ for p_T > 100 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 λ for p_T > 100 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 λ for p_T > 100 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 λ for p_T > 100 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 λ for p_T > 100 MeV as a function of m_ch in k_T intervals between 0.1 and 0.5 GeV. The error bars and boxes represent the statistical and systematic contributions, respectively.

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

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

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

The parameter R for p_T > 100 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 > 100 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 > 100 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 > 100 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 > 100 MeV as a function of m_ch in k_T intervals between 0.1 and 0.5 GeV. The error bars and boxes represent the statistical and systematic contributions, respectively.

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

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

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

The fit parameter μ describing the dependence of the correlation strength, λ, on charged-particle scaled multiplicity, for track p_T>100 MeV and track p_T>500 MeV in the minimum-bias (MB) and high-multiplicity track (HMT) samples at √s = 13 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 μ on m_ch for tracks with p_T >100 MeV (p_T >500 MeV).

The fit parameter μ describing the dependence of the correlation strength, λ, on charged-particle scaled multiplicity, for track p_T>100 MeV and track p_T>500 MeV in the minimum-bias (MB) and high-multiplicity track (HMT) samples at √s = 13 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 μ on m_ch for tracks with p_T >100 MeV (p_T >500 MeV).

The fit parameter μ describing the dependence of the correlation strength, λ, on charged-particle scaled multiplicity, for track p_T>100 MeV and track p_T>500 MeV in the minimum-bias (MB) and high-multiplicity track (HMT) samples at √s = 13 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 μ on m_ch for tracks with p_T >100 MeV (p_T >500 MeV).

The fit parameter μ describing the dependence of the correlation strength, λ, on charged-particle scaled multiplicity, for track p_T>100 MeV and track p_T>500 MeV in the minimum-bias (MB) and high-multiplicity track (HMT) samples at √s = 13 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 μ on m_ch for tracks with p_T >100 MeV (p_T >500 MeV).

The fit parameter ν describing the dependence of the correlation strength, λ, on charged-particle scaled multiplicity, for track p_T>100 MeV and track p_T>500 MeV in the minimum-bias (MB) and high-multiplicity track (HMT) samples at √s = 13 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 ν on m_ch for tracks with p_T >100 MeV (p_T >500 MeV).

The fit parameter ν describing the dependence of the correlation strength, λ, on charged-particle scaled multiplicity, for track p_T>100 MeV and track p_T>500 MeV in the minimum-bias (MB) and high-multiplicity track (HMT) samples at √s = 13 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 ν on m_ch for tracks with p_T >100 MeV (p_T >500 MeV).

The fit parameter ν describing the dependence of the correlation strength, λ, on charged-particle scaled multiplicity, for track p_T>100 MeV and track p_T>500 MeV in the minimum-bias (MB) and high-multiplicity track (HMT) samples at √s = 13 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 ν on m_ch for tracks with p_T >100 MeV (p_T >500 MeV).

The fit parameter ν describing the dependence of the correlation strength, λ, on charged-particle scaled multiplicity, for track p_T>100 MeV and track p_T>500 MeV in the minimum-bias (MB) and high-multiplicity track (HMT) samples at √s = 13 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 ν on m_ch for tracks with p_T >100 MeV (p_T >500 MeV).

The parameter ξ describing the dependence of the source radius, R, on charged-particle scaled multiplicity, m_ch, for track p_T>100 MeV and track p_T>500 MeV in the minimum-bias (MB) and high-multiplicity track (HMT) samples at √s = 13 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 ξ for m_ch > 3.0 for tracks with p_T >100 MeV and for m_ch > 2.8 for tracks with p_T >500 MeV, respectively.

The parameter ξ describing the dependence of the source radius, R, on charged-particle scaled multiplicity, m_ch, for track p_T>100 MeV and track p_T>500 MeV in the minimum-bias (MB) and high-multiplicity track (HMT) samples at √s = 13 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 ξ for m_ch > 3.0 for tracks with p_T >100 MeV and for m_ch > 2.8 for tracks with p_T >500 MeV, respectively.

The parameter ξ describing the dependence of the source radius, R, on charged-particle scaled multiplicity, m_ch, for track p_T>100 MeV and track p_T>500 MeV in the minimum-bias (MB) and high-multiplicity track (HMT) samples at √s = 13 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 ξ for m_ch > 3.0 for tracks with p_T >100 MeV and for m_ch > 2.8 for tracks with p_T >500 MeV, respectively.

The parameter ξ describing the dependence of the source radius, R, on charged-particle scaled multiplicity, m_ch, for track p_T>100 MeV and track p_T>500 MeV in the minimum-bias (MB) and high-multiplicity track (HMT) samples at √s = 13 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 ξ for m_ch > 3.0 for tracks with p_T >100 MeV and for m_ch > 2.8 for tracks with p_T >500 MeV, respectively.

The parameter κ describing the dependence of the source radius, R, on charged-particle scaled multiplicity, m_ch, for track p_T>100 MeV and track p_T>500 MeV in the minimum-bias (MB) and high-multiplicity track (HMT) samples at √s = 13 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 κ for tracks with p_T >100 MeV and for tracks with p_T >500 MeV, respectively.

The parameter κ describing the dependence of the source radius, R, on charged-particle scaled multiplicity, m_ch, for track p_T>100 MeV and track p_T>500 MeV in the minimum-bias (MB) and high-multiplicity track (HMT) samples at √s = 13 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 κ for tracks with p_T >100 MeV and for tracks with p_T >500 MeV, respectively.

The parameter κ describing the dependence of the source radius, R, on charged-particle scaled multiplicity, m_ch, for track p_T>100 MeV and track p_T>500 MeV in the minimum-bias (MB) and high-multiplicity track (HMT) samples at √s = 13 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 κ for tracks with p_T >100 MeV and for tracks with p_T >500 MeV, respectively.

The parameter κ describing the dependence of the source radius, R, on charged-particle scaled multiplicity, m_ch, for track p_T>100 MeV and track p_T>500 MeV in the minimum-bias (MB) and high-multiplicity track (HMT) samples at √s = 13 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 κ for tracks with p_T >100 MeV and for tracks with p_T >500 MeV, respectively.

The two-dimensional dependence on m_ch and k_T for p_T > 500 MeV for the correlation strength, λ, obtained from the exponential fit to the R₂(Q) correlation functions using the MB sample for m_ch ≤ 3.08 and the HMT sample for m_ch > 3.08.

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

The parameter λ for p_T > 500 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 λ for p_T > 500 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 λ for p_T > 500 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 λ for p_T > 500 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 λ for p_T > 500 MeV as a function of m_ch in k_T intervals between 0.5 and 1.5 GeV. The error bars and boxes represent the statistical and systematic contributions, respectively.

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

The parameter R for p_T > 500 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 > 500 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 > 500 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 > 500 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 > 500 MeV as a function of m_ch in k_T intervals between 0.5 and 1.5 GeV. The error bars and boxes represent the statistical and systematic contributions, respectively.

The parameter R for p_T > 500 MeV as a function of m_ch in k_T intervals between 0.5 and 1.5 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 TeV. The CMS results (open circles) have been adjusted (by the CMS collaboration) to the ATLAS kinematic region∶ p_T > 100 MeV and |η|<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 ∛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 TeV. The CMS results (open circles) have been adjusted (by the CMS collaboration) to the ATLAS kinematic region∶ p_T > 100 MeV and |η|<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 ∛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 TeV. The CMS results (open circles) have been adjusted (by the CMS collaboration) to the ATLAS kinematic region∶ p_T > 100 MeV and |η|<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 ∛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 TeV. The CMS results (open circles) have been adjusted (by the CMS collaboration) to the ATLAS kinematic region∶ p_T > 100 MeV and |η|<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 ∛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 TeV. The CMS results (open circles) have been adjusted (by the CMS collaboration) to the ATLAS kinematic region∶ p_T > 100 MeV and |η|<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 ∛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 TeV. The CMS results (open circles) have been adjusted (by the CMS collaboration) to the ATLAS kinematic region∶ p_T > 100 MeV and |η|<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 ∛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, λ, and source radius, R, for the exponential fit of the two-particle double-ratio correlation functions, R₂(Q), for p_T > 100 MeV at √s= 13 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, λ, and source radius, R, on the average rescaled charged-particle multiplicity, m_ch, for |η| < 2.5 and both p_T > 100 MeV and p_T > 500 MeV at √s = 13 TeV for the minimum-bias (MB) and the high-multiplicity track (HMT) events. The parameters γ and δ 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, λ, 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 > 100 MeV and p_T > 500 MeV at √s = 13 TeV. The total uncertainties are shown.

The Bose-Einstein correlation (BEC) parameters λ and R as a function of n_ch and k_T using different MC generators in the calculation of R₂(Q). (a) λ versus n_ch for MB events, (b) λ versus n_ch for HMT events, (c) λ 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 λ and R as a function of n_ch and k_T using different MC generators in the calculation of R₂(Q). (a) λ versus n_ch for MB events, (b) λ versus n_ch for HMT events, (c) λ 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 λ and R as a function of n_ch and k_T using different MC generators in the calculation of R₂(Q). (a) λ versus n_ch for MB events, (b) λ versus n_ch for HMT events, (c) λ 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 λ and R as a function of n_ch and k_T using different MC generators in the calculation of R₂(Q). (a) λ versus n_ch for MB events, (b) λ versus n_ch for HMT events, (c) λ 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∶ (a) 2 < n_ch ≤ 10, (b) 11 < n_ch ≤ 20, (c) 21 < n_ch ≤ 30, (d) 31 < n_ch ≤ 40, (e) 41 < n_ch ≤ 50, (f) 51 < n_ch ≤ 60, (g) 61 < n_ch ≤ 70, (h) 71 < n_ch ≤ 80 and (i) 81 < n_ch ≤ 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∶ (a) 2 < n_ch ≤ 10, (b) 11 < n_ch ≤ 20, (c) 21 < n_ch ≤ 30, (d) 31 < n_ch ≤ 40, (e) 41 < n_ch ≤ 50, (f) 51 < n_ch ≤ 60, (g) 61 < n_ch ≤ 70, (h) 71 < n_ch ≤ 80 and (i) 81 < n_ch ≤ 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∶ (a) 2 < n_ch ≤ 10, (b) 11 < n_ch ≤ 20, (c) 21 < n_ch ≤ 30, (d) 31 < n_ch ≤ 40, (e) 41 < n_ch ≤ 50, (f) 51 < n_ch ≤ 60, (g) 61 < n_ch ≤ 70, (h) 71 < n_ch ≤ 80 and (i) 81 < n_ch ≤ 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∶ (a) 2 < n_ch ≤ 10, (b) 11 < n_ch ≤ 20, (c) 21 < n_ch ≤ 30, (d) 31 < n_ch ≤ 40, (e) 41 < n_ch ≤ 50, (f) 51 < n_ch ≤ 60, (g) 61 < n_ch ≤ 70, (h) 71 < n_ch ≤ 80 and (i) 81 < n_ch ≤ 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∶ (a) 2 < n_ch ≤ 10, (b) 11 < n_ch ≤ 20, (c) 21 < n_ch ≤ 30, (d) 31 < n_ch ≤ 40, (e) 41 < n_ch ≤ 50, (f) 51 < n_ch ≤ 60, (g) 61 < n_ch ≤ 70, (h) 71 < n_ch ≤ 80 and (i) 81 < n_ch ≤ 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∶ (a) 2 < n_ch ≤ 10, (b) 11 < n_ch ≤ 20, (c) 21 < n_ch ≤ 30, (d) 31 < n_ch ≤ 40, (e) 41 < n_ch ≤ 50, (f) 51 < n_ch ≤ 60, (g) 61 < n_ch ≤ 70, (h) 71 < n_ch ≤ 80 and (i) 81 < n_ch ≤ 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∶ (a) 2 < n_ch ≤ 10, (b) 11 < n_ch ≤ 20, (c) 21 < n_ch ≤ 30, (d) 31 < n_ch ≤ 40, (e) 41 < n_ch ≤ 50, (f) 51 < n_ch ≤ 60, (g) 61 < n_ch ≤ 70, (h) 71 < n_ch ≤ 80 and (i) 81 < n_ch ≤ 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∶ (a) 2 < n_ch ≤ 10, (b) 11 < n_ch ≤ 20, (c) 21 < n_ch ≤ 30, (d) 31 < n_ch ≤ 40, (e) 41 < n_ch ≤ 50, (f) 51 < n_ch ≤ 60, (g) 61 < n_ch ≤ 70, (h) 71 < n_ch ≤ 80 and (i) 81 < n_ch ≤ 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∶ (a) 2 < n_ch ≤ 10, (b) 11 < n_ch ≤ 20, (c) 21 < n_ch ≤ 30, (d) 31 < n_ch ≤ 40, (e) 41 < n_ch ≤ 50, (f) 51 < n_ch ≤ 60, (g) 61 < n_ch ≤ 70, (h) 71 < n_ch ≤ 80 and (i) 81 < n_ch ≤ 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∶ (a) 91 < n_ch ≤ 100, (b) 101 < n_ch ≤ 125, (c) 126 < n_ch ≤ 150, (d) 151 < n_ch ≤ 200, (e) 201 < n_ch ≤ 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∶ (a) 91 < n_ch ≤ 100, (b) 101 < n_ch ≤ 125, (c) 126 < n_ch ≤ 150, (d) 151 < n_ch ≤ 200, (e) 201 < n_ch ≤ 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∶ (a) 91 < n_ch ≤ 100, (b) 101 < n_ch ≤ 125, (c) 126 < n_ch ≤ 150, (d) 151 < n_ch ≤ 200, (e) 201 < n_ch ≤ 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∶ (a) 91 < n_ch ≤ 100, (b) 101 < n_ch ≤ 125, (c) 126 < n_ch ≤ 150, (d) 151 < n_ch ≤ 200, (e) 201 < n_ch ≤ 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∶ (a) 91 < n_ch ≤ 100, (b) 101 < n_ch ≤ 125, (c) 126 < n_ch ≤ 150, (d) 151 < n_ch ≤ 200, (e) 201 < n_ch ≤ 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∶ (a) 101 < n_ch ≤ 110, (b) 111 < n_ch ≤ 120, (c) 121 < n_ch ≤ 130, (d) 131 < n_ch ≤ 140, (e) 141 < n_ch ≤ 155, (f) 156 < n_ch ≤ 175, (g) 176 < n_ch ≤ 200, (h) 201 < n_ch ≤ 230 and (i) 231 < n_ch ≤ 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∶ (a) 101 < n_ch ≤ 110, (b) 111 < n_ch ≤ 120, (c) 121 < n_ch ≤ 130, (d) 131 < n_ch ≤ 140, (e) 141 < n_ch ≤ 155, (f) 156 < n_ch ≤ 175, (g) 176 < n_ch ≤ 200, (h) 201 < n_ch ≤ 230 and (i) 231 < n_ch ≤ 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∶ (a) 101 < n_ch ≤ 110, (b) 111 < n_ch ≤ 120, (c) 121 < n_ch ≤ 130, (d) 131 < n_ch ≤ 140, (e) 141 < n_ch ≤ 155, (f) 156 < n_ch ≤ 175, (g) 176 < n_ch ≤ 200, (h) 201 < n_ch ≤ 230 and (i) 231 < n_ch ≤ 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∶ (a) 101 < n_ch ≤ 110, (b) 111 < n_ch ≤ 120, (c) 121 < n_ch ≤ 130, (d) 131 < n_ch ≤ 140, (e) 141 < n_ch ≤ 155, (f) 156 < n_ch ≤ 175, (g) 176 < n_ch ≤ 200, (h) 201 < n_ch ≤ 230 and (i) 231 < n_ch ≤ 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∶ (a) 101 < n_ch ≤ 110, (b) 111 < n_ch ≤ 120, (c) 121 < n_ch ≤ 130, (d) 131 < n_ch ≤ 140, (e) 141 < n_ch ≤ 155, (f) 156 < n_ch ≤ 175, (g) 176 < n_ch ≤ 200, (h) 201 < n_ch ≤ 230 and (i) 231 < n_ch ≤ 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∶ (a) 101 < n_ch ≤ 110, (b) 111 < n_ch ≤ 120, (c) 121 < n_ch ≤ 130, (d) 131 < n_ch ≤ 140, (e) 141 < n_ch ≤ 155, (f) 156 < n_ch ≤ 175, (g) 176 < n_ch ≤ 200, (h) 201 < n_ch ≤ 230 and (i) 231 < n_ch ≤ 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∶ (a) 101 < n_ch ≤ 110, (b) 111 < n_ch ≤ 120, (c) 121 < n_ch ≤ 130, (d) 131 < n_ch ≤ 140, (e) 141 < n_ch ≤ 155, (f) 156 < n_ch ≤ 175, (g) 176 < n_ch ≤ 200, (h) 201 < n_ch ≤ 230 and (i) 231 < n_ch ≤ 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∶ (a) 101 < n_ch ≤ 110, (b) 111 < n_ch ≤ 120, (c) 121 < n_ch ≤ 130, (d) 131 < n_ch ≤ 140, (e) 141 < n_ch ≤ 155, (f) 156 < n_ch ≤ 175, (g) 176 < n_ch ≤ 200, (h) 201 < n_ch ≤ 230 and (i) 231 < n_ch ≤ 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∶ (a) 101 < n_ch ≤ 110, (b) 111 < n_ch ≤ 120, (c) 121 < n_ch ≤ 130, (d) 131 < n_ch ≤ 140, (e) 141 < n_ch ≤ 155, (f) 156 < n_ch ≤ 175, (g) 176 < n_ch ≤ 200, (h) 201 < n_ch ≤ 230 and (i) 231 < n_ch ≤ 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∶ (a) 100 < k_T ≤ 200 MeV, (b) 200 < k_T ≤ 300 MeV, (c) 300 < k_T ≤ 400 MeV, (d) 400 < k_T ≤ 500 MeV, (e) 500 < k_T ≤ 600 MeV, (f) 600 < k_T ≤ 700 MeV, (g) 700 < k_T ≤ 1000 MeV, and (h) 1000 < k_T ≤ 1500 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∶ (a) 100 < k_T ≤ 200 MeV, (b) 200 < k_T ≤ 300 MeV, (c) 300 < k_T ≤ 400 MeV, (d) 400 < k_T ≤ 500 MeV, (e) 500 < k_T ≤ 600 MeV, (f) 600 < k_T ≤ 700 MeV, (g) 700 < k_T ≤ 1000 MeV, and (h) 1000 < k_T ≤ 1500 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∶ (a) 100 < k_T ≤ 200 MeV, (b) 200 < k_T ≤ 300 MeV, (c) 300 < k_T ≤ 400 MeV, (d) 400 < k_T ≤ 500 MeV, (e) 500 < k_T ≤ 600 MeV, (f) 600 < k_T ≤ 700 MeV, (g) 700 < k_T ≤ 1000 MeV, and (h) 1000 < k_T ≤ 1500 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∶ (a) 100 < k_T ≤ 200 MeV, (b) 200 < k_T ≤ 300 MeV, (c) 300 < k_T ≤ 400 MeV, (d) 400 < k_T ≤ 500 MeV, (e) 500 < k_T ≤ 600 MeV, (f) 600 < k_T ≤ 700 MeV, (g) 700 < k_T ≤ 1000 MeV, and (h) 1000 < k_T ≤ 1500 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∶ (a) 100 < k_T ≤ 200 MeV, (b) 200 < k_T ≤ 300 MeV, (c) 300 < k_T ≤ 400 MeV, (d) 400 < k_T ≤ 500 MeV, (e) 500 < k_T ≤ 600 MeV, (f) 600 < k_T ≤ 700 MeV, (g) 700 < k_T ≤ 1000 MeV, and (h) 1000 < k_T ≤ 1500 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∶ (a) 100 < k_T ≤ 200 MeV, (b) 200 < k_T ≤ 300 MeV, (c) 300 < k_T ≤ 400 MeV, (d) 400 < k_T ≤ 500 MeV, (e) 500 < k_T ≤ 600 MeV, (f) 600 < k_T ≤ 700 MeV, (g) 700 < k_T ≤ 1000 MeV, and (h) 1000 < k_T ≤ 1500 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∶ (a) 100 < k_T ≤ 200 MeV, (b) 200 < k_T ≤ 300 MeV, (c) 300 < k_T ≤ 400 MeV, (d) 400 < k_T ≤ 500 MeV, (e) 500 < k_T ≤ 600 MeV, (f) 600 < k_T ≤ 700 MeV, (g) 700 < k_T ≤ 1000 MeV, and (h) 1000 < k_T ≤ 1500 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∶ (a) 100 < k_T ≤ 200 MeV, (b) 200 < k_T ≤ 300 MeV, (c) 300 < k_T ≤ 400 MeV, (d) 400 < k_T ≤ 500 MeV, (e) 500 < k_T ≤ 600 MeV, (f) 600 < k_T ≤ 700 MeV, (g) 700 < k_T ≤ 1000 MeV, and (h) 1000 < k_T ≤ 1500 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∶ (a) 100 < k_T ≤ 200 MeV, (b) 200 < k_T ≤ 300 MeV, (c) 300 < k_T ≤ 400 MeV, (d) 400 < k_T ≤ 500 MeV, (e) 500 < k_T ≤ 600 MeV, (f) 600 < k_T ≤ 700 MeV, (g) 700 < k_T ≤ 1000 MeV, and (h) 1000 < k_T ≤ 1500 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∶ (a) 100 < k_T ≤ 200 MeV, (b) 200 < k_T ≤ 300 MeV, (c) 300 < k_T ≤ 400 MeV, (d) 400 < k_T ≤ 500 MeV, (e) 500 < k_T ≤ 600 MeV, (f) 600 < k_T ≤ 700 MeV, (g) 700 < k_T ≤ 1000 MeV, and (h) 1000 < k_T ≤ 1500 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∶ (a) 100 < k_T ≤ 200 MeV, (b) 200 < k_T ≤ 300 MeV, (c) 300 < k_T ≤ 400 MeV, (d) 400 < k_T ≤ 500 MeV, (e) 500 < k_T ≤ 600 MeV, (f) 600 < k_T ≤ 700 MeV, (g) 700 < k_T ≤ 1000 MeV, and (h) 1000 < k_T ≤ 1500 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∶ (a) 100 < k_T ≤ 200 MeV, (b) 200 < k_T ≤ 300 MeV, (c) 300 < k_T ≤ 400 MeV, (d) 400 < k_T ≤ 500 MeV, (e) 500 < k_T ≤ 600 MeV, (f) 600 < k_T ≤ 700 MeV, (g) 700 < k_T ≤ 1000 MeV, and (h) 1000 < k_T ≤ 1500 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∶ (a) 100 < k_T ≤ 200 MeV, (b) 200 < k_T ≤ 300 MeV, (c) 300 < k_T ≤ 400 MeV, (d) 400 < k_T ≤ 500 MeV, (e) 500 < k_T ≤ 600 MeV, (f) 600 < k_T ≤ 700 MeV, (g) 700 < k_T ≤ 1000 MeV, and (h) 1000 < k_T ≤ 1500 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∶ (a) 100 < k_T ≤ 200 MeV, (b) 200 < k_T ≤ 300 MeV, (c) 300 < k_T ≤ 400 MeV, (d) 400 < k_T ≤ 500 MeV, (e) 500 < k_T ≤ 600 MeV, (f) 600 < k_T ≤ 700 MeV, (g) 700 < k_T ≤ 1000 MeV, and (h) 1000 < k_T ≤ 1500 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∶ (a) 100 < k_T ≤ 200 MeV, (b) 200 < k_T ≤ 300 MeV, (c) 300 < k_T ≤ 400 MeV, (d) 400 < k_T ≤ 500 MeV, (e) 500 < k_T ≤ 600 MeV, (f) 600 < k_T ≤ 700 MeV, (g) 700 < k_T ≤ 1000 MeV, and (h) 1000 < k_T ≤ 1500 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∶ (a) 100 < k_T ≤ 200 MeV, (b) 200 < k_T ≤ 300 MeV, (c) 300 < k_T ≤ 400 MeV, (d) 400 < k_T ≤ 500 MeV, (e) 500 < k_T ≤ 600 MeV, (f) 600 < k_T ≤ 700 MeV, (g) 700 < k_T ≤ 1000 MeV, and (h) 1000 < k_T ≤ 1500 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∶ (a) 2 < n_ch ≤ 9, (b) 10 < n_ch ≤ 18, (c) 19 < n_ch ≤ 27, (d) 28 < n_ch ≤ 36, (e) 37 < n_ch ≤ 45, (f) 46 < n_ch ≤ 54, (g) 55 < n_ch ≤ 63, (h) 64 < n_ch ≤ 72, (i) 73 < n_ch ≤ 81, (j) 82 < n_ch ≤ 90, (k) 91 < n_ch ≤ 113, and (l) 114 < n_ch ≤ 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∶ (a) 2 < n_ch ≤ 9, (b) 10 < n_ch ≤ 18, (c) 19 < n_ch ≤ 27, (d) 28 < n_ch ≤ 36, (e) 37 < n_ch ≤ 45, (f) 46 < n_ch ≤ 54, (g) 55 < n_ch ≤ 63, (h) 64 < n_ch ≤ 72, (i) 73 < n_ch ≤ 81, (j) 82 < n_ch ≤ 90, (k) 91 < n_ch ≤ 113, and (l) 114 < n_ch ≤ 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∶ (a) 2 < n_ch ≤ 9, (b) 10 < n_ch ≤ 18, (c) 19 < n_ch ≤ 27, (d) 28 < n_ch ≤ 36, (e) 37 < n_ch ≤ 45, (f) 46 < n_ch ≤ 54, (g) 55 < n_ch ≤ 63, (h) 64 < n_ch ≤ 72, (i) 73 < n_ch ≤ 81, (j) 82 < n_ch ≤ 90, (k) 91 < n_ch ≤ 113, and (l) 114 < n_ch ≤ 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∶ (a) 2 < n_ch ≤ 9, (b) 10 < n_ch ≤ 18, (c) 19 < n_ch ≤ 27, (d) 28 < n_ch ≤ 36, (e) 37 < n_ch ≤ 45, (f) 46 < n_ch ≤ 54, (g) 55 < n_ch ≤ 63, (h) 64 < n_ch ≤ 72, (i) 73 < n_ch ≤ 81, (j) 82 < n_ch ≤ 90, (k) 91 < n_ch ≤ 113, and (l) 114 < n_ch ≤ 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∶ (a) 2 < n_ch ≤ 9, (b) 10 < n_ch ≤ 18, (c) 19 < n_ch ≤ 27, (d) 28 < n_ch ≤ 36, (e) 37 < n_ch ≤ 45, (f) 46 < n_ch ≤ 54, (g) 55 < n_ch ≤ 63, (h) 64 < n_ch ≤ 72, (i) 73 < n_ch ≤ 81, (j) 82 < n_ch ≤ 90, (k) 91 < n_ch ≤ 113, and (l) 114 < n_ch ≤ 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∶ (a) 2 < n_ch ≤ 9, (b) 10 < n_ch ≤ 18, (c) 19 < n_ch ≤ 27, (d) 28 < n_ch ≤ 36, (e) 37 < n_ch ≤ 45, (f) 46 < n_ch ≤ 54, (g) 55 < n_ch ≤ 63, (h) 64 < n_ch ≤ 72, (i) 73 < n_ch ≤ 81, (j) 82 < n_ch ≤ 90, (k) 91 < n_ch ≤ 113, and (l) 114 < n_ch ≤ 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∶ (a) 2 < n_ch ≤ 9, (b) 10 < n_ch ≤ 18, (c) 19 < n_ch ≤ 27, (d) 28 < n_ch ≤ 36, (e) 37 < n_ch ≤ 45, (f) 46 < n_ch ≤ 54, (g) 55 < n_ch ≤ 63, (h) 64 < n_ch ≤ 72, (i) 73 < n_ch ≤ 81, (j) 82 < n_ch ≤ 90, (k) 91 < n_ch ≤ 113, and (l) 114 < n_ch ≤ 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∶ (a) 2 < n_ch ≤ 9, (b) 10 < n_ch ≤ 18, (c) 19 < n_ch ≤ 27, (d) 28 < n_ch ≤ 36, (e) 37 < n_ch ≤ 45, (f) 46 < n_ch ≤ 54, (g) 55 < n_ch ≤ 63, (h) 64 < n_ch ≤ 72, (i) 73 < n_ch ≤ 81, (j) 82 < n_ch ≤ 90, (k) 91 < n_ch ≤ 113, and (l) 114 < n_ch ≤ 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∶ (a) 2 < n_ch ≤ 9, (b) 10 < n_ch ≤ 18, (c) 19 < n_ch ≤ 27, (d) 28 < n_ch ≤ 36, (e) 37 < n_ch ≤ 45, (f) 46 < n_ch ≤ 54, (g) 55 < n_ch ≤ 63, (h) 64 < n_ch ≤ 72, (i) 73 < n_ch ≤ 81, (j) 82 < n_ch ≤ 90, (k) 91 < n_ch ≤ 113, and (l) 114 < n_ch ≤ 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∶ (a) 2 < n_ch ≤ 9, (b) 10 < n_ch ≤ 18, (c) 19 < n_ch ≤ 27, (d) 28 < n_ch ≤ 36, (e) 37 < n_ch ≤ 45, (f) 46 < n_ch ≤ 54, (g) 55 < n_ch ≤ 63, (h) 64 < n_ch ≤ 72, (i) 73 < n_ch ≤ 81, (j) 82 < n_ch ≤ 90, (k) 91 < n_ch ≤ 113, and (l) 114 < n_ch ≤ 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∶ (a) 2 < n_ch ≤ 9, (b) 10 < n_ch ≤ 18, (c) 19 < n_ch ≤ 27, (d) 28 < n_ch ≤ 36, (e) 37 < n_ch ≤ 45, (f) 46 < n_ch ≤ 54, (g) 55 < n_ch ≤ 63, (h) 64 < n_ch ≤ 72, (i) 73 < n_ch ≤ 81, (j) 82 < n_ch ≤ 90, (k) 91 < n_ch ≤ 113, and (l) 114 < n_ch ≤ 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∶ (a) 2 < n_ch ≤ 9, (b) 10 < n_ch ≤ 18, (c) 19 < n_ch ≤ 27, (d) 28 < n_ch ≤ 36, (e) 37 < n_ch ≤ 45, (f) 46 < n_ch ≤ 54, (g) 55 < n_ch ≤ 63, (h) 64 < n_ch ≤ 72, (i) 73 < n_ch ≤ 81, (j) 82 < n_ch ≤ 90, (k) 91 < n_ch ≤ 113, and (l) 114 < n_ch ≤ 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∶ (a) 100 < k_T ≤ 200 MeV, (b) 200 < k_T ≤ 300 MeV, (c) 300 < k_T ≤ 400 MeV, (d) 400 < k_T ≤ 500 MeV, (e) 500 < k_T ≤ 600 MeV, (f) 600 < k_T ≤ 700 MeV, (g) 700 < k_T ≤ 1000 MeV, and (h) 1000 < k_T ≤ 1500 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∶ (a) 100 < k_T ≤ 200 MeV, (b) 200 < k_T ≤ 300 MeV, (c) 300 < k_T ≤ 400 MeV, (d) 400 < k_T ≤ 500 MeV, (e) 500 < k_T ≤ 600 MeV, (f) 600 < k_T ≤ 700 MeV, (g) 700 < k_T ≤ 1000 MeV, and (h) 1000 < k_T ≤ 1500 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∶ (a) 100 < k_T ≤ 200 MeV, (b) 200 < k_T ≤ 300 MeV, (c) 300 < k_T ≤ 400 MeV, (d) 400 < k_T ≤ 500 MeV, (e) 500 < k_T ≤ 600 MeV, (f) 600 < k_T ≤ 700 MeV, (g) 700 < k_T ≤ 1000 MeV, and (h) 1000 < k_T ≤ 1500 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∶ (a) 100 < k_T ≤ 200 MeV, (b) 200 < k_T ≤ 300 MeV, (c) 300 < k_T ≤ 400 MeV, (d) 400 < k_T ≤ 500 MeV, (e) 500 < k_T ≤ 600 MeV, (f) 600 < k_T ≤ 700 MeV, (g) 700 < k_T ≤ 1000 MeV, and (h) 1000 < k_T ≤ 1500 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∶ (a) 100 < k_T ≤ 200 MeV, (b) 200 < k_T ≤ 300 MeV, (c) 300 < k_T ≤ 400 MeV, (d) 400 < k_T ≤ 500 MeV, (e) 500 < k_T ≤ 600 MeV, (f) 600 < k_T ≤ 700 MeV, (g) 700 < k_T ≤ 1000 MeV, and (h) 1000 < k_T ≤ 1500 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∶ (a) 100 < k_T ≤ 200 MeV, (b) 200 < k_T ≤ 300 MeV, (c) 300 < k_T ≤ 400 MeV, (d) 400 < k_T ≤ 500 MeV, (e) 500 < k_T ≤ 600 MeV, (f) 600 < k_T ≤ 700 MeV, (g) 700 < k_T ≤ 1000 MeV, and (h) 1000 < k_T ≤ 1500 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∶ (a) 100 < k_T ≤ 200 MeV, (b) 200 < k_T ≤ 300 MeV, (c) 300 < k_T ≤ 400 MeV, (d) 400 < k_T ≤ 500 MeV, (e) 500 < k_T ≤ 600 MeV, (f) 600 < k_T ≤ 700 MeV, (g) 700 < k_T ≤ 1000 MeV, and (h) 1000 < k_T ≤ 1500 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∶ (a) 100 < k_T ≤ 200 MeV, (b) 200 < k_T ≤ 300 MeV, (c) 300 < k_T ≤ 400 MeV, (d) 400 < k_T ≤ 500 MeV, (e) 500 < k_T ≤ 600 MeV, (f) 600 < k_T ≤ 700 MeV, (g) 700 < k_T ≤ 1000 MeV, and (h) 1000 < k_T ≤ 1500 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, λ, 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 > 100 MeV at √s = 13 TeV. Statistical uncertainties for √χ²/ndf>1 are corrected by the √χ²/ndf. The total uncertainties are shown.

The correlation strength, λ, 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 > 500 MeV at √s = 13 TeV. Statistical uncertainties for √χ²/ndf>1 are corrected by the √χ²/ndf. The total uncertainties are shown.

The correlation strength, λ, 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 > 100 MeV at √s = 13 TeV. Statistical uncertainties for √χ²/ndf>1 are corrected by the √χ²/ndf. The total uncertainties are shown.

The correlation strength, λ, 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 > 500 MeV at √s = 13 TeV. Statistical uncertainties for √χ²/ndf>1 are corrected by the √χ²/ndf. The total uncertainties are shown.

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.

Pre-fit distribution of the reconstructed Higgs boson candidate $p_T^H$ for the dilepton $SR^{\geq 4j}_{\geq 4b}$ signal region. The dashed line shows the ${t\bar {t}H}$ signal distribution normalised to the total background prediction. The uncertainty band includes all uncertainties and their correlations, except for the uncertainty in the $k({t\bar {t}+{\geq }1b})$ normalisation factor which is not defined pre-fit. The last bin includes the overflow.

Pre-fit distribution of the reconstructed Higgs boson candidate $p_T^H$ for the single-lepton resolved $SR^{\geq 6j}_{\geq 4b}$ signal region. The dashed line shows the ${t\bar {t}H}$ signal distribution normalised to the total background prediction. The uncertainty band includes all uncertainties and their correlations, except for the uncertainty in the $k({t\bar {t}+{\geq }1b})$ normalisation factor which is not defined pre-fit. The last bin includes the overflow.

Pre-fit distribution of the reconstructed Higgs boson candidate $p_T^H$ for the single-lepton boosted ${{SR}_{{boosted}}}$ signal region. The dashed line shows the ${t\bar {t}H}$ signal distribution normalised to the total background prediction. The uncertainty band includes all uncertainties and their correlations, except for the uncertainty in the $k({t\bar {t}+{\geq }1b})$ normalisation factor which is not defined pre-fit. The last bin includes the overflow.

Comparison of predicted and observed event yields in each of the control and signal regions in the dilepton channel after the fit to the data. The uncertainty band includes all uncertainties and their correlations.

Comparison of predicted and observed event yields in each of the control and signal regions in the single-lepton channels after the fit to the data. The uncertainty band includes all uncertainties and their correlations.

Comparison between data and prediction for the BDT discriminant in the dilepton SRs after the inclusive fit to the data for $0\le p_T^H<120$ GeV. The ${t\bar {t}H}$ signal yield (solid red) is normalised to the fitted $\mu $ value from the inclusive fit. The dashed line shows the ${t\bar {t}H}$ signal distribution normalised to the total background prediction. The uncertainty band includes all uncertainties and their correlations.

Comparison between data and prediction for the BDT discriminant in the dilepton SRs after the inclusive fit to the data for $120\le p_T^H<200$ GeV. The ${t\bar {t}H}$ signal yield (solid red) is normalised to the fitted $\mu $ value from the inclusive fit. The dashed line shows the ${t\bar {t}H}$ signal distribution normalised to the total background prediction. The uncertainty band includes all uncertainties and their correlations.

Comparison between data and prediction for the BDT discriminant in the dilepton SRs after the inclusive fit to the data for $200\le p_T^H<300$ GeV. The ${t\bar {t}H}$ signal yield (solid red) is normalised to the fitted $\mu $ value from the inclusive fit. The dashed line shows the ${t\bar {t}H}$ signal distribution normalised to the total background prediction. The uncertainty band includes all uncertainties and their correlations.

Comparison between data and prediction for the BDT discriminant in the dilepton SRs after the inclusive fit to the data for $p_{{T}}^{H}\ge 300$ GeV. The ${t\bar {t}H}$ signal yield (solid red) is normalised to the fitted $\mu $ value from the inclusive fit. The dashed line shows the ${t\bar {t}H}$ signal distribution normalised to the total background prediction. The uncertainty band includes all uncertainties and their correlations.

Comparison between data and prediction for the BDT discriminant in the single-lepton resolved SRs after the inclusive fit to the data for $0\le p_T^H<120$ GeV. The ${t\bar {t}H}$ signal yield (solid red) is normalised to the fitted $\mu $ value from the inclusive fit. The dashed line shows the ${t\bar {t}H}$ signal distribution normalised to the total background prediction. The uncertainty band includes all uncertainties and their correlations.

Comparison between data and prediction for the BDT discriminant in the single-lepton resolved SRs after the inclusive fit to the data for $120\le p_T^H<200$ GeV. The ${t\bar {t}H}$ signal yield (solid red) is normalised to the fitted $\mu $ value from the inclusive fit. The dashed line shows the ${t\bar {t}H}$ signal distribution normalised to the total background prediction. The uncertainty band includes all uncertainties and their correlations.

Comparison between data and prediction for the BDT discriminant in the single-lepton resolved SRs after the inclusive fit to the data for $200\le p_T^H<300$ GeV. The ${t\bar {t}H}$ signal yield (solid red) is normalised to the fitted $\mu $ value from the inclusive fit. The dashed line shows the ${t\bar {t}H}$ signal distribution normalised to the total background prediction. The uncertainty band includes all uncertainties and their correlations.

Comparison between data and prediction for the BDT discriminant in the single-lepton resolved SRs after the inclusive fit to the data for $300\le p_T^H<450$ GeV. The ${t\bar {t}H}$ signal yield (solid red) is normalised to the fitted $\mu $ value from the inclusive fit. The dashed line shows the ${t\bar {t}H}$ signal distribution normalised to the total background prediction. The uncertainty band includes all uncertainties and their correlations.

Comparison between data and prediction for the BDT discriminant in the single-lepton resolved SRs after the inclusive fit to the data for $p_{{T}}^{H}\ge 450$ GeV (yield only). The ${t\bar {t}H}$ signal yield (solid red) is normalised to the fitted $\mu $ value from the inclusive fit. The dashed line shows the ${t\bar {t}H}$ signal distribution normalised to the total background prediction. The uncertainty band includes all uncertainties and their correlations.

Comparison between data and prediction for the BDT discriminant in the single-lepton boosted SRs after the inclusive fit to the data for $300\le p_T^H<450$ GeV. The ${t\bar {t}H}$ signal yield (solid red) is normalised to the fitted $\mu $ value from the inclusive fit. The dashed line shows the ${t\bar {t}H}$ signal distribution normalised to the total background prediction. The uncertainty band includes all uncertainties and their correlations.

Comparison between data and prediction for the BDT discriminant in the single-lepton boosted SRs after the inclusive fit to the data for $p_{{T}}^{H}\ge 450$ GeV. The ${t\bar {t}H}$ signal yield (solid red) is normalised to the fitted $\mu $ value from the inclusive fit. The dashed line shows the ${t\bar {t}H}$ signal distribution normalised to the total background prediction. The uncertainty band includes all uncertainties and their correlations.

Comparison between data and prediction for ${\Delta R^{{avg}}_{bb}}$ after the inclusive fit to the data in the single-lepton $CR^{5j}_{{\geq}4b\ lo}$ control region. The ${t\bar {t}H}$ signal yield (solid red) is normalised to the fitted $\mu $ value from the inclusive fit. The dashed line shows the ${t\bar {t}H}$ signal distribution normalised to the total background prediction. The uncertainty band includes all uncertainties and their correlations. The first (last) bin includes the underflow (overflow).

Comparison between data and prediction for ${\Delta R^{{avg}}_{bb}}$ after the inclusive fit to the data in the single-lepton $CR^{5j}_{{\geq}4b\ hi}$ control region. The ${t\bar {t}H}$ signal yield (solid red) is normalised to the fitted $\mu $ value from the inclusive fit. The dashed line shows the ${t\bar {t}H}$ signal distribution normalised to the total background prediction. The uncertainty band includes all uncertainties and their correlations. The first (last) bin includes the underflow (overflow).

Post-fit yields of signal ($S$) and total background ($B$) as a function of $\log (S/B)$, compared with data. Final-discriminant bins in all dilepton and single-lepton analysis regions are combined into bins of $\log (S/B)$, with the signal normalised to the SM prediction used for the computation of $\log (S/B)$. The signal is then shown normalised to the best-fit value and the SM prediction. The lower frame reports the ratio of data to background, and this is compared with the expected ${t\bar {t}H}$-signal-plus-background yield divided by the background-only yield for the best-fit signal strength (solid red line) and the SM prediction (dashed orange line).

Comparison between data and prediction for the reconstruction BDT score for the Higgs boson candidate identified using Higgs boson information, after the inclusive fit to the data in the dilepton resolved channel for $0\le p_T^H<120$ GeV. The ${t\bar {t}H}$ signal yield (solid red) is normalised to the fitted $\mu $ value from the inclusive fit. The dashed line shows the ${t\bar {t}H}$ signal distribution normalised to the total background prediction. The uncertainty band includes all uncertainties and their correlations.

Comparison between data and prediction for the average $\Delta \eta $ between $b$-tagged jets, after the inclusive fit to the data in the dilepton resolved channel for $0\le p_T^H<120$ GeV. The ${t\bar {t}H}$ signal yield (solid red) is normalised to the fitted $\mu $ value from the inclusive fit. The dashed line shows the ${t\bar {t}H}$ signal distribution normalised to the total background prediction. The uncertainty band includes all uncertainties and their correlations.

Comparison between data and prediction for the likelihood discriminant, after the inclusive fit to the data in the single-lepton resolved channel for $0\le p_T^H<120$ GeV. The ${t\bar {t}H}$ signal yield (solid red) is normalised to the fitted $\mu $ value from the inclusive fit. The dashed line shows the ${t\bar {t}H}$ signal distribution normalised to the total background prediction. The uncertainty band includes all uncertainties and their correlations.

Comparison between data and prediction for the average $\Delta R$ for all possible combinations of $b$-tagged jet pairs, after the inclusive fit to the data in the single-lepton resolved channel for $0\le p_T^H<120$ GeV. The ${t\bar {t}H}$ signal yield (solid red) is normalised to the fitted $\mu $ value from the inclusive fit. The dashed line shows the ${t\bar {t}H}$ signal distribution normalised to the total background prediction. The uncertainty band includes all uncertainties and their correlations.

Comparison between data and prediction for the DNN $P(H)$ output for the Higgs boson candidate after the inclusive fit to the data in the single-lepton boosted channel for $300\le p_T^H<450$ GeV. The ${t\bar {t}H}$ signal yield (solid red) is normalised to the fitted $\mu $ value from the inclusive fit. The dashed line shows the ${t\bar {t}H}$ signal distribution normalised to the total background prediction. The uncertainty band includes all uncertainties and their correlations.

Comparison between data and prediction for the DNN $P(H)$ output for the Higgs boson candidate after the inclusive fit to the data in the single-lepton boosted channel for $p_{{T}}^{H}\ge 450$ GeV. The ${t\bar {t}H}$ signal yield (solid red) is normalised to the fitted $\mu $ value from the inclusive fit. The dashed line shows the ${t\bar {t}H}$ signal distribution normalised to the total background prediction. The uncertainty band includes all uncertainties and their correlations.

Post-fit distribution of the reconstructed Higgs boson candidate mass for the dilepton $SR^{\geq 4j}_{\geq 4b}$ signal region. The ${t\bar {t}H}$ signal yield (solid red) is normalised to the fitted $\mu $ value from the inclusive fit. The dashed line shows the ${t\bar {t}H}$ signal distribution normalised to the total background prediction. The uncertainty band includes all uncertainties and their correlations. The first (last) bin includes the underflow (overflow).

Post-fit distribution of the reconstructed Higgs boson candidate mass for the single-lepton resolved $SR^{\geq 6j}_{\geq 4b}$ signal region. The ${t\bar {t}H}$ signal yield (solid red) is normalised to the fitted $\mu $ value from the inclusive fit. The dashed line shows the ${t\bar {t}H}$ signal distribution normalised to the total background prediction. The uncertainty band includes all uncertainties and their correlations. The first (last) bin includes the underflow (overflow).

Post-fit distribution of the reconstructed Higgs boson candidate mass for the single-lepton boosted ${{SR}_{{boosted}}}$ signal region. The ${t\bar {t}H}$ signal yield (solid red) is normalised to the fitted $\mu $ value from the inclusive fit. The dashed line shows the ${t\bar {t}H}$ signal distribution normalised to the total background prediction. The uncertainty band includes all uncertainties and their correlations. The first (last) bin includes the underflow (overflow).

Fitted values of the ${t\bar {t}H}$ signal strength parameter in the individual channels and in the inclusive signal-strength measurement.

Ranking of the 20 nuisance parameters with the largest post-fit impact on $\mu $ in the fit. Nuisance parameters corresponding to statistical uncertainties in the simulated event samples are not included. The empty blue rectangles correspond to the pre-fit impact on $\mu $ and the filled blue ones to the post-fit impact on $\mu $, both referring to the upper scale. The impact of each nuisance parameter, $\Delta \mu $, is computed by comparing the nominal best-fit value of $\mu $ with the result of the fit when fixing the considered nuisance parameter to its best-fit value, $\hat{\theta }$, shifted by its pre-fit (post-fit) uncertainties $\pm \Delta \theta $ ($\pm \Delta \hat{\theta }$). The black points show the pulls of the nuisance parameters relative to their nominal values, $\theta _0$. These pulls and their relative post-fit errors, $\Delta \hat{\theta }/\Delta \theta $, refer to the lower scale. The `ljets' (`dilep') label refers to the single-lepton (dilepton) channel.

Pre-fit distribution of the number of jets in the dilepton $SR^{\geq 4j}_{\geq 4b}$ signal region. The ${t\bar {t}H}$ signal yield (solid red) is normalised to the Standard Model expectation. The dashed line shows the ${t\bar {t}H}$ signal distribution normalised to the total background prediction. The uncertainty band includes all uncertainties and their correlations, except the uncertainty in the $k({t\bar {t}+{\geq }1b})$ normalisation factor that is not defined pre-fit.

Pre-fit distribution of the number of jets in the single-lepton resolved $SR^{\geq 6j}_{\geq 4b}$ signal region. The ${t\bar {t}H}$ signal yield (solid red) is normalised to the Standard Model expectation. The dashed line shows the ${t\bar {t}H}$ signal distribution normalised to the total background prediction. The uncertainty band includes all uncertainties and their correlations, except the uncertainty in the $k({t\bar {t}+{\geq }1b})$ normalisation factor that is not defined pre-fit.

Pre-fit distribution of the number of jets in the single-lepton boosted ${{SR}_{{boosted}}}$ signal region. The ${t\bar {t}H}$ signal yield (solid red) is normalised to the Standard Model expectation. The dashed line shows the ${t\bar {t}H}$ signal distribution normalised to the total background prediction. The uncertainty band includes all uncertainties and their correlations, except the uncertainty in the $k({t\bar {t}+{\geq }1b})$ normalisation factor that is not defined pre-fit.

Post-fit distribution of the number of jets in the dilepton $SR^{\geq 4j}_{\geq 4b}$ signal region. The ${t\bar {t}H}$ signal yield (solid red) is normalised to the fitted $\mu $ value from the inclusive fit. The dashed line shows the ${t\bar {t}H}$ signal distribution normalised to the total background prediction. The uncertainty band includes all uncertainties and their correlations.

Post-fit distribution of the number of jets in the single-lepton resolved $SR^{\geq 6j}_{\geq 4b}$ signal region. The ${t\bar {t}H}$ signal yield (solid red) is normalised to the fitted $\mu $ value from the inclusive fit. The dashed line shows the ${t\bar {t}H}$ signal distribution normalised to the total background prediction. The uncertainty band includes all uncertainties and their correlations.

Post-fit distribution of the number of jets in the single-lepton boosted ${{SR}_{{boosted}}}$ signal region. The ${t\bar {t}H}$ signal yield (solid red) is normalised to the fitted $\mu $ value from the inclusive fit. The dashed line shows the ${t\bar {t}H}$ signal distribution normalised to the total background prediction. The uncertainty band includes all uncertainties and their correlations.

Post-fit distribution of the reconstructed Higgs boson candidate $p_T^H$ for the dilepton $SR^{\geq 4j}_{\geq 4b}$ signal region. The ${t\bar {t}H}$ signal yield (solid red) is normalised to the fitted $\mu $ value from the inclusive fit. The dashed line shows the ${t\bar {t}H}$ signal distribution normalised to the total background prediction. The uncertainty band includes all uncertainties and their correlations. The last bin includes the overflow.

Post-fit distribution of the reconstructed Higgs boson candidate $p_T^H$ for the single-lepton resolved $SR^{\geq 6j}_{\geq 4b}$ signal region. The ${t\bar {t}H}$ signal yield (solid red) is normalised to the fitted $\mu $ value from the inclusive fit. The dashed line shows the ${t\bar {t}H}$ signal distribution normalised to the total background prediction. The uncertainty band includes all uncertainties and their correlations. The last bin includes the overflow.

Post-fit distribution of the reconstructed Higgs boson candidate $p_T^H$ for the single-lepton boosted ${{SR}_{{boosted}}}$ signal region. The ${t\bar {t}H}$ signal yield (solid red) is normalised to the fitted $\mu $ value from the inclusive fit. The dashed line shows the ${t\bar {t}H}$ signal distribution normalised to the total background prediction. The uncertainty band includes all uncertainties and their correlations. The last bin includes the overflow.

Signal-strength measurements in the individual STXS ${\hat{p}_{{T}}^{H}}$ bins, as well as the inclusive signal strength.

95% CL simplified template cross-section upper limits in the individual STXS ${\hat{p}_{{T}}^{H}}$ bins, as well as the inclusive limit. The observed limits are shown (solid black lines), together with the expected limits both in the background-only hypothesis (dotted black lines) and in the SM hypothesis (dotted red lines). In the case of the expected limits in the background-only hypothesis, one- and two-standard-deviation uncertainty bands are also shown. The hatched uncertainty bands correspond to the theory uncertainty in the fiducial cross-section prediction in each bin.

The ratios $S/B$ (black solid line, referring to the vertical axis on the left) and $S/\sqrt{B}$ (red dashed line, referring to the vertical axis on the right) for each category in the inclusive analysis in the dilepton channel (left) and in the single-lepton channels (right), where $S$ ($B$) is the number of selected signal (background) events predicted by the simulation and normalised to a luminosity of 139 fb$^{-1}$ .

Comparison between data and prediction for the $\Delta R$ between the Higgs candidate and the ${t\bar {t}}$ candidate system, after the inclusive fit to the data in the dilepton resolved channel for $0\le p_T^H<120$ GeV. The ${t\bar {t}H}$ signal yield (solid red) is normalised to the fitted $\mu $ value from the inclusive fit. The dashed line shows the ${t\bar {t}H}$ signal distribution normalised to the total background prediction. The uncertainty band includes all uncertainties and their correlations.

Comparison between data and prediction for the number of $b$-tagged jet pairs with an invariant mass within 30 GeV of 125 GeV, after the inclusive fit to the data in the dilepton resolved channel for $0\le p_T^H<120$ GeV. The ${t\bar {t}H}$ signal yield (solid red) is normalised to the fitted $\mu $ value from the inclusive fit. The dashed line shows the ${t\bar {t}H}$ signal distribution normalised to the total background prediction. The uncertainty band includes all uncertainties and their correlations.

Comparison between data and prediction for the reconstruction BDT score for the Higgs boson candidate identified using Higgs boson information, after the inclusive fit to the data in the single-lepton resolved channel for $0\le p_T^H<120$ GeV. The ${t\bar {t}H}$ signal yield (solid red) is normalised to the fitted $\mu $ value from the inclusive fit. The dashed line shows the ${t\bar {t}H}$ signal distribution normalised to the total background prediction. The uncertainty band includes all uncertainties and their correlations.

Comparison between data and prediction for the $\Delta R$ between the two highest ${p_{{T}}}$ $b$-tagged jets, after the inclusive fit to the data in the single-lepton resolved channel for $0\le p_T^H<120$ GeV. The ${t\bar {t}H}$ signal yield (solid red) is normalised to the fitted $\mu $ value from the inclusive fit. The dashed line shows the ${t\bar {t}H}$ signal distribution normalised to the total background prediction. The uncertainty band includes all uncertainties and their correlations.

Comparison between data and prediction for the sum of $b$-tagging discriminants of jets from Higgs, hadronic top and leptonic top candidates, after the inclusive fit to the data in the single-lepton boosted channel for $300\le p_T^H<450$ GeV. The ${t\bar {t}H}$ signal yield (solid red) is normalised to the fitted $\mu $ value from the inclusive fit. The dashed line shows the ${t\bar {t}H}$ signal distribution normalised to the total background prediction. The uncertainty band includes all uncertainties and their correlations. The first (last) bin includes the underflow (overflow).

Comparison between data and prediction for the sum of $b$-tagging discriminants of jets from Higgs, hadronic top and leptonic top candidates, after the inclusive fit to the data in the single-lepton boosted channel for $p_{{T}}^{H}\ge 450$ GeV. The ${t\bar {t}H}$ signal yield (solid red) is normalised to the fitted $\mu $ value from the inclusive fit. The dashed line shows the ${t\bar {t}H}$ signal distribution normalised to the total background prediction. The uncertainty band includes all uncertainties and their correlations. The first (last) bin includes the underflow (overflow).

Comparison between data and prediction for the hadronic top candidate invariant mass, after the inclusive fit to the data in the single-lepton boosted channel for $300\le p_T^H<450$ GeV. The ${t\bar {t}H}$ signal yield (solid red) is normalised to the fitted $\mu $ value from the inclusive fit. The dashed line shows the ${t\bar {t}H}$ signal distribution normalised to the total background prediction. The uncertainty band includes all uncertainties and their correlations. The first (last) bin includes the underflow (overflow).

Comparison between data and prediction for the hadronic top candidate invariant mass, after the inclusive fit to the data in the single-lepton boosted channel for $p_{{T}}^{H}\ge 450$ GeV. The ${t\bar {t}H}$ signal yield (solid red) is normalised to the fitted $\mu $ value from the inclusive fit. The dashed line shows the ${t\bar {t}H}$ signal distribution normalised to the total background prediction. The uncertainty band includes all uncertainties and their correlations. The first (last) bin includes the underflow (overflow).

Comparison between data and prediction for the fraction of the sum of $b$-tagging discriminants of all jets not associated to the Higgs or hadronic top candidates, after the inclusive fit to the data in the single-lepton boosted channel for $300\le p_T^H<450$ GeV. The ${t\bar {t}H}$ signal yield (solid red) is normalised to the fitted $\mu $ value from the inclusive fit. The dashed line shows the ${t\bar {t}H}$ signal distribution normalised to the total background prediction. The uncertainty band includes all uncertainties and their correlations. The first (last) bin includes the underflow (overflow).

Comparison between data and prediction for the fraction of the sum of $b$-tagging discriminants of all jets not associated to the Higgs or hadronic top candidates, after the inclusive fit to the data in the single-lepton boosted channel for $p_{{T}}^{H}\ge 450$ GeV. The ${t\bar {t}H}$ signal yield (solid red) is normalised to the fitted $\mu $ value from the inclusive fit. The dashed line shows the ${t\bar {t}H}$ signal distribution normalised to the total background prediction. The uncertainty band includes all uncertainties and their correlations. The first (last) bin includes the underflow (overflow).

Ranking of the 20 nuisance parameters with the largest post-fit impact on $\mu $ in the STXS fit for $0\le {\hat{p}_{{T}}^{H}}<120$ GeV. Nuisance parameters corresponding to statistical uncertainties in the simulated event samples are not included. The empty blue rectangles correspond to the pre-fit impact on $\mu $ and the filled blue ones to the post-fit impact on $\mu $, both referring to the upper scale. The impact of each nuisance parameter, $\Delta \mu $, is computed by comparing the nominal best-fit value of $\mu $ with the result of the fit when fixing the considered nuisance parameter to its best-fit value, $\hat{\theta }$, shifted by its pre-fit (post-fit) uncertainties $\pm \Delta \theta $ ($\pm \Delta \hat{\theta }$). The black points show the pulls of the nuisance parameters relative to their nominal values, $\theta _0$. These pulls and their relative post-fit errors, $\Delta \hat{\theta }/\Delta \theta $, refer to the lower scale. For experimental uncertainties that are decomposed into several independent sources, NP X corresponds to the X$^{th}$ nuisance parameter, ordered by their impact on $\mu $. The `ljets' (`dilep') label refers to the single-lepton (dilepton) channel.

Ranking of the 20 nuisance parameters with the largest post-fit impact on $\mu $ in the STXS fit for $120\le {\hat{p}_{{T}}^{H}}<200$ GeV. Nuisance parameters corresponding to statistical uncertainties in the simulated event samples are not included. The empty blue rectangles correspond to the pre-fit impact on $\mu $ and the filled blue ones to the post-fit impact on $\mu $, both referring to the upper scale. The impact of each nuisance parameter, $\Delta \mu $, is computed by comparing the nominal best-fit value of $\mu $ with the result of the fit when fixing the considered nuisance parameter to its best-fit value, $\hat{\theta }$, shifted by its pre-fit (post-fit) uncertainties $\pm \Delta \theta $ ($\pm \Delta \hat{\theta }$). The black points show the pulls of the nuisance parameters relative to their nominal values, $\theta _0$. These pulls and their relative post-fit errors, $\Delta \hat{\theta }/\Delta \theta $, refer to the lower scale. For experimental uncertainties that are decomposed into several independent sources, NP X corresponds to the X$^{th}$ nuisance parameter, ordered by their impact on $\mu $. The `ljets' (`dilep') label refers to the single-lepton (dilepton) channel.

Ranking of the 20 nuisance parameters with the largest post-fit impact on $\mu $ in the STXS fit for $200\le {\hat{p}_{{T}}^{H}}<300$ GeV. Nuisance parameters corresponding to statistical uncertainties in the simulated event samples are not included. The empty blue rectangles correspond to the pre-fit impact on $\mu $ and the filled blue ones to the post-fit impact on $\mu $, both referring to the upper scale. The impact of each nuisance parameter, $\Delta \mu $, is computed by comparing the nominal best-fit value of $\mu $ with the result of the fit when fixing the considered nuisance parameter to its best-fit value, $\hat{\theta }$, shifted by its pre-fit (post-fit) uncertainties $\pm \Delta \theta $ ($\pm \Delta \hat{\theta }$). The black points show the pulls of the nuisance parameters relative to their nominal values, $\theta _0$. These pulls and their relative post-fit errors, $\Delta \hat{\theta }/\Delta \theta $, refer to the lower scale. For experimental uncertainties that are decomposed into several independent sources, NP X corresponds to the X$^{th}$ nuisance parameter, ordered by their impact on $\mu $. The `ljets' (`dilep') label refers to the single-lepton (dilepton) channel.

Ranking of the 20 nuisance parameters with the largest post-fit impact on $\mu $ in the STXS fit for $300\le {\hat{p}_{{T}}^{H}}<450$ GeV. Nuisance parameters corresponding to statistical uncertainties in the simulated event samples are not included. The empty blue rectangles correspond to the pre-fit impact on $\mu $ and the filled blue ones to the post-fit impact on $\mu $, both referring to the upper scale. The impact of each nuisance parameter, $\Delta \mu $, is computed by comparing the nominal best-fit value of $\mu $ with the result of the fit when fixing the considered nuisance parameter to its best-fit value, $\hat{\theta }$, shifted by its pre-fit (post-fit) uncertainties $\pm \Delta \theta $ ($\pm \Delta \hat{\theta }$). The black points show the pulls of the nuisance parameters relative to their nominal values, $\theta _0$. These pulls and their relative post-fit errors, $\Delta \hat{\theta }/\Delta \theta $, refer to the lower scale. For experimental uncertainties that are decomposed into several independent sources, NP X corresponds to the X$^{th}$ nuisance parameter, ordered by their impact on $\mu $. The `ljets' (`dilep') label refers to the single-lepton (dilepton) channel.

Ranking of the 20 nuisance parameters with the largest post-fit impact on $\mu $ in the STXS fit for ${\hat{p}_{{T}}^{H}}\ge 450$ GeV. Nuisance parameters corresponding to statistical uncertainties in the simulated event samples are not included. The empty blue rectangles correspond to the pre-fit impact on $\mu $ and the filled blue ones to the post-fit impact on $\mu $, both referring to the upper scale. The impact of each nuisance parameter, $\Delta \mu $, is computed by comparing the nominal best-fit value of $\mu $ with the result of the fit when fixing the considered nuisance parameter to its best-fit value, $\hat{\theta }$, shifted by its pre-fit (post-fit) uncertainties $\pm \Delta \theta $ ($\pm \Delta \hat{\theta }$). The black points show the pulls of the nuisance parameters relative to their nominal values, $\theta _0$. These pulls and their relative post-fit errors, $\Delta \hat{\theta }/\Delta \theta $, refer to the lower scale. For experimental uncertainties that are decomposed into several independent sources, NP X corresponds to the X$^{th}$ nuisance parameter, ordered by their impact on $\mu $. The `ljets' (`dilep') label refers to the single-lepton (dilepton) channel.

95% confidence level upper limits on signal-strength measurements in the individual STXS ${\hat{p}_{{T}}^{H}}$ bins, as well as the inclusive signal-strength limit, after the fit used to extract multiple signal-strength parameters. The observed limits are shown (solid black lines), together with the expected limits both in the background-only hypothesis (dotted black lines) and in the SM hypothesis (dotted red lines). In the case of the expected limits in the background-only hypothesis, one- and two-standard-deviation uncertainty bands are also shown.

Post-fit correlation matrix (in percentages) between the $\mu $ values obtained in the STXS bins.

Performance of the Higgs boson reconstruction algorithms. For each row of `truth' ${\hat{p}_{{T}}^{H}}$, the matrix shows (in percentages) the fraction of Higgs boson candidates which are truth-matched to ${b\bar {b}}$ decays, with reconstructed $p_T^H$ in the various bins of the dilepton (left), single lepton resolved (middle) and boosted (right) channels.

Pre-fit event yields in the dilepton signal regions and control regions. All uncertainties are included except the $k({t\bar {t}+{\geq }1b})$ uncertainty that is not defined pre-fit. For the ${t\bar {t}H}$ signal, the pre-fit yield values correspond to the theoretical prediction and corresponding uncertainties. `Other sources' refers to s-channel, t-channel, $tW$, $tWZ$, $tZq$, $Z+$ jets and diboson events.

Post-fit event yields in the dilepton signal regions and control regions, after the inclusive fit in all channels. All uncertainties are included, taking into account correlations. For the ${t\bar {t}H}$ signal, the post-fit yield and uncertainties correspond to those in the inclusive signal-strength measurement. `Other sources' refers to s-channel, t-channel, $tW$, $tWZ$, $tZq$, $Z+$ jets and diboson events.

Pre-fit event yields in the single-lepton resolved and boosted signal regions and control regions. All uncertainties are included except the $k({t\bar {t}+{\geq }1b})$ uncertainty that is not defined pre-fit. For the ${t\bar {t}H}$ signal, the pre-fit yield values correspond to the theoretical prediction and corresponding uncertainties. `Other top sources' refers to s-channel, t-channel, $tWZ$ and $tZq$ events.

Post-fit event yields in the single-lepton resolved and boosted signal regions and control regions, after the inclusive fit in all channels. All uncertainties are included, taking into account correlations. For the ${t\bar {t}H}$ signal, the post-fit yield and uncertainties correspond to those in the inclusive signal-strength measurement. `Other top sources' refers to s-channel, t-channel, $tWZ$ and $tZq$ events.

Breakdown of the contributions to the uncertainties in $\mu$. The contributions from the different sources of uncertainty are evaluated after the fit. The $\Delta \mu $ values are obtained by repeating the fit after having fixed a certain set of nuisance parameters corresponding to a group of systematic uncertainties, and then evaluating $(\Delta \mu)^2$ by subtracting the resulting squared uncertainty of $\mu $ from its squared uncertainty found in the full fit. The same procedure is followed when quoting the effect of the ${t\bar {t}+{\geq }1b}$ normalisation. The total uncertainty is different from the sum in quadrature of the different components due to correlations between nuisance parameters existing in the fit.

Fraction (in percentages) of signal events, after SR and CR selections, originating from $b\bar {b}$, $WW$ and other remaining Higgs boson decay modes in the dilepton channel.

Fraction (in percentages) of signal events, after SR and CR selections, originating from $b\bar {b}$, $WW$ and other remaining Higgs boson decay modes in the single-lepton channels.

Predicted SM ${t\bar {t}H}$ cross-section in each of the five STXS ${\hat{p}_{{T}}^{H}}$ bins and signal acceptance times efficiency (including all event selection criteria) in each STXS bin as well as for the inclusive ${\hat{p}_{{T}}^{H}}$ range.

Number of expected signal events before the fit, after each selection requirement applied to enter the dilepton channel $SR^{\geq 4j}_{\geq 4b}$ region. All ${t\bar {t}H}$ signal events are included, regardless of the $H$ or ${t\bar {t}H}$ decay mode. All object corrections are applied, except for the initial number of events which is calculated using the NLO QCD+EW theoretical prediction. All quoted numbers are rounded to unity. More details on the selection criteria can be found in the text.

Number of expected signal events before the fit, after each selection requirement applied to enter the single-lepton channel resolved $SR^{\geq 6j}_{\geq 4b}$ region. All ${t\bar {t}H}$ signal events are included, regardless of the $H$ or ${t\bar {t}H}$ decay mode. All object corrections are applied, except for the initial number of events which is calculated using the NLO QCD+EW theoretical prediction. All quoted numbers are rounded to unity. More details on the selection criteria can be found in the text.

Number of expected signal events before the fit, after each selection requirement applied to enter the single-lepton channel boosted $SR_{boosted}$ region. All ${t\bar {t}H}$ signal events are included, regardless of the $H$ or ${t\bar {t}H}$ decay mode. All object corrections are applied, except for the initial number of events which is calculated using the NLO QCD+EW theoretical prediction. All quoted numbers are rounded to unity. More details on the selection criteria can be found in the text.

$K^-$ (a), invariant yields as a function of $m_T-m_0$ for various rapidity regions in 10--40\% central Au+Au collisions at ${\sqrt{s_{\mathrm{NN}}} = \mathrm{3\,GeV}}$. Statistics and systematic uncertainties are added quadratic here for plotting. Solid and dashed black lines depict $m_T$ exponential function fits to the measured data points with arbitrate scaling factors in each rapidity windows.

$\phi$ meson (b) invariant yields as a function of $m_T-m_0$ for various rapidity regions in 10--40\% central Au+Au collisions at ${\sqrt{s_{\mathrm{NN}}} = \mathrm{3\,GeV}}$. Statistics and systematic uncertainties are added quadratic here for plotting. Solid and dashed black lines depict $m_T$ exponential function fits to the measured data points with arbitrate scaling factors in each rapidity windows.

$\Xi^-$ (c) invariant yields as a function of $m_T-m_0$ for various rapidity regions in 10--40\% central Au+Au collisions at ${\sqrt{s_{\mathrm{NN}}} = \mathrm{3\,GeV}}$. Statistics and systematic uncertainties are added quadratic here for plotting. Solid and dashed black lines depict $m_T$ exponential function fits to the measured data points with arbitrate scaling factors in each rapidity windows.

$K^-$ (a), invariant yields as a function of $m_T-m_0$ for various rapidity regions in 40--60\% central Au+Au collisions at ${\sqrt{s_{\mathrm{NN}}} = \mathrm{3\,GeV}}$. Statistics and systematic uncertainties are added quadratic here for plotting. Solid and dashed black lines depict $m_T$ exponential function fits to the measured data points with arbitrate scaling factors in each rapidity windows.

$\phi$ meson (b) invariant yields as a function of $m_T-m_0$ for various rapidity regions in 40--60\% central Au+Au collisions at ${\sqrt{s_{\mathrm{NN}}} = \mathrm{3\,GeV}}$. Statistics and systematic uncertainties are added quadratic here for plotting. Solid and dashed black lines depict $m_T$ exponential function fits to the measured data points with arbitrate scaling factors in each rapidity windows.

Rapidity density distributions of $K^-$ (squares), $p_T$-integrated yields $dN/dy$ in 0--10\% (a), 10--40\% (b) and 40--60\% central (c) Au+Au collisions at ${\sqrt{s_{\mathrm{NN}}} = \mathrm{3\,GeV}}$. %The full symbols show the measured data, while the open ones are data reflected with respect to $y=0$ in the center-of-mass frame. Solid lines depict Gaussian function fits to the data points.

Rapidity density distributions of $\phi$ meson (circles) $p_T$-integrated yields $dN/dy$ in 0--10\% (a), 10--40\% (b) and 40--60\% central (c) Au+Au collisions at ${\sqrt{s_{\mathrm{NN}}} = \mathrm{3\,GeV}}$. %The full symbols show the measured data, while the open ones are data reflected with respect to $y=0$ in the center-of-mass frame. Solid lines depict Gaussian function fits to the data points.

Rapidity density distributions of $\Xi^-$ (diamonds) $p_T$-integrated yields $dN/dy$ in 0--10\% (a), 10--40\% (b) and 40--60\% central (c) Au+Au collisions at ${\sqrt{s_{\mathrm{NN}}} = \mathrm{3\,GeV}}$. %The full symbols show the measured data, while the open ones are data reflected with respect to $y=0$ in the center-of-mass frame. Solid lines depict Gaussian function fits to the data points.

$\phi/K^-$ (a) and $\phi/\Xi^-$ (b) ratio as a function of collision energy, $\sqrt{s_{\mathrm{NN}}}$. The solid black circles show the measurements presented here in 0-10\% centrality bin, while empty markers in black or grey are used for data from various other energies and/or collision systems. The grey solid line represents a THERMUS calculation based on the Grand Canonical Ensemble (GCE) while the dotted lines depict calculations based on the Canonical Ensemble (CE) with different parameters of strangeness correlation radius ($r_c$). The green dashed line, green shaded band and the solid red line show transport model calculations from the public versions $\mathrm{UrQMD}^{1}$, modified $\mathrm{UrQMD}^{2}$ and SMASH, respectively.

J/PSI(1S)-->MU+MU- nuclear modification in p+Al collisions as a function of rapidity. The statistical and systematic uncertainties vary point-to-point and are listed for each measured value. An additional global systematic uncertainty is provided in each column heading, which applies to all data points per column.

PSI(2S)-->MU+MU- nuclear modification in p+Au collisions as a function of rapidity. The statistical and systematic uncertainties vary point-to-point and are listed for each measured value. An additional global systematic uncertainty is provided in each column heading, which applies to all data points per column.

J/PSI(1S)-->MU+MU- nuclear modification in p+Au collisions as a function of rapidity. The statistical and systematic uncertainties vary point-to-point and are listed for each measured value. An additional global systematic uncertainty is provided in each column heading, which applies to all data points per column.

PSI(2S)-->MU+MU- nuclear modification in p+Au collisions as a function of centrality at forward rapidity. The statistical and systematic uncertainties vary point-to-point and are listed for each measured value. An additional global systematic uncertainty is provided in each column heading, which applies to all data points per column.

PSI(2S)-->MU+MU- nuclear modification in p+Au collisions as a function of centrality at backward rapidity. The statistical and systematic uncertainties vary point-to-point and are listed for each measured value. An additional global systematic uncertainty is provided in each column heading, which applies to all data points per column.

J/PSI(1S)-->MU+MU- nuclear modification in p+Au collisions as a function of centrality. The statistical and systematic uncertainties vary point-to-point and are listed for each measured value. An additional global systematic uncertainty is provided in each column heading, which applies to all data points per column.

Invariant yield ratio of PSI(2S)-->MU+MU- with respect to J/PSI(1S)-->MU+MU- in p+p collisions at rapidity 1.2 < |y| < 2.2. The statistical and systematic uncertainties vary point-to-point and are listed for each measured value. The global systematic uncertainty cancels when taking the ratio.

Invariant yield ratio of PSI(2S)-->MU+MU- with respect to J/PSI(1S)-->MU+MU- in p+Au collisions as a function of centrality at forward rapidity. The statistical and systematic uncertainties vary point-to-point and are listed for each measured value. The global systematic uncertainty cancels when taking the ratio.

Invariant yield ratio of PSI(2S)-->MU+MU- with respect to J/PSI(1S)-->MU+MU- in p+Au collisions as a function of centrality at backward rapidity. The statistical and systematic uncertainties vary point-to-point and are listed for each measured value. The global systematic uncertainty cancels when taking the ratio.

$p_{\rm T}$-differential yields of $\Lambda$(1520) (sum of particle and anti-particle states) in p--Pb collisions at $\sqrt{s_{\mathrm{NN}}}$ $\mathrm{=}$ 5.02 TeV in multiplicity interval 20--40\%. The uncertainty 'sys,$p_{\rm T}$-correlated' indicates the systematic uncertainty after removing the contributions of $p_{\rm T}$-uncorrelated uncertainty.

$p_{\rm T}$-differential yields of $\Lambda$(1520) (sum of particle and anti-particle states) in p--Pb collisions at $\sqrt{s_{\mathrm{NN}}}$ $\mathrm{=}$ 5.02 TeV in multiplicity interval 40--60\%. The uncertainty 'sys,$p_{\rm T}$-correlated' indicates the systematic uncertainty after removing the contributions of $p_{\rm T}$-uncorrelated uncertainty.

$p_{\rm T}$-differential yields of $\Lambda$(1520) (sum of particle and anti-particle states) in p--Pb collisions at $\sqrt{s_{\mathrm{NN}}}$ $\mathrm{=}$ 5.02 TeV in multiplicity interval 60--100\%. The uncertainty 'sys,$p_{\rm T}$-correlated' indicates the systematic uncertainty after removing the contributions of $p_{\rm T}$-uncorrelated uncertainty.

Ratio of $p_{\rm T}$-integrated yields of $\Lambda$(1520) (sum of particle and anti-particle states) to charged $\pi$ ($\pi^{+} + \pi^{-}$) at midrapidity as a function of $\langle {\rm d}N_{\rm ch}/{\rm d} \eta_{\rm lab} \rangle$ in inelastic pp collisions at $\sqrt{s}$ $\mathrm{=}$ 7 TeV.

Ratio of $p_{\rm T}$-integrated yields of $\Lambda$(1520) (sum of particle and anti-particle states) to charged $\pi$ ($\pi^{+} + \pi^{-}$) as a function of $\langle {\rm d}N_{\rm ch}/{\rm d} \eta_{\rm lab} \rangle$ in p--Pb collisions at $\sqrt{s}$ $\mathrm{=}$ 5.02 TeV. The uncertainty 'sys,uncorrelated' indicates the multiplicity uncorrelated systematic uncertainty.

Ratio of $p_{\rm T}$-integrated yields of $\Lambda$(1520) (sum of particle and anti-particle states) to K$^{-}$ at midrapidity as a function of $\langle {\rm d}N_{\rm ch}/{\rm d} \eta_{\rm lab} \rangle$ in inelastic pp collisions at $\sqrt{s}$ $\mathrm{=}$ 7 TeV.

Ratio of $p_{\rm T}$-integrated yields of $\Lambda$(1520) (sum of particle and anti-particle states) to K$^{-}$ as a function of $\langle {\rm d}N_{\rm ch}/{\rm d} \eta_{\rm lab} \rangle$ in p--Pb collisions at $\sqrt{s}$ $\mathrm{=}$ 5.02 TeV. The uncertainty 'sys,uncorrelated' indicates the multiplicity uncorrelated systematic uncertainty.

Ratio of $p_{\rm T}$-integrated yields of $\Lambda$(1520) (sum of particle and anti-particle states) to its ground state particle $\Lambda$ ($\Lambda + \bar{\Lambda}$) at midrapidity as a function of $\langle {\rm d}N_{\rm ch}/{\rm d} \eta_{\rm lab} \rangle$ in inelastic pp collisions at $\sqrt{s}$ $\mathrm{=}$ 7 TeV.

Ratio of $p_{\rm T}$-integrated yields of $\Lambda$(1520) (sum of particle and anti-particle states) to its ground state particle $\Lambda$ ($\Lambda + \bar{\Lambda}$) as a function of $\langle {\rm d}N_{\rm ch}/{\rm d} \eta_{\rm lab} \rangle$ in p--Pb collisions at $\sqrt{s}$ $\mathrm{=}$ 5.02 TeV. The uncertainty 'sys,uncorrelated' indicates the multiplicity uncorrelated systematic uncertainty.

Azimuthal correlation distribution between heavy-flavour decay electrons and charged hadrons, scaled by the number of electrons in EMCal triggered events in the electron transverse momentum range 4.5-6 GeV/c.

Relative beauty contribution to the heavy-flavour electron yield obtained with the method based on the track impact parameter.

Relative beauty contribution to the heavy-flavour electron yield obtained with the method based on the track impact parameter.

Relative beauty contribution to the heavy-flavour electron yield obtained from the method based on the azimuthal correlation of heavy-flavour decay electrons and charged tracks.

Relative beauty contribution to the heavy-flavour electron yield obtained from the method based on the azimuthal correlation of heavy-flavour decay electrons and charged tracks.

The $p_{\rm T}$-differential inclusive cross section of electrons (multiplied by $10^3$ for readability) from beauty hadron decays using the electron-hadron azimuthal correlation method.

The pT-differential inclusive cross section of electrons from beauty hadron decays using the electron-hadron azimuthal correlation method.

The $p_{\rm T}$-differential inclusive cross section of electrons (multiplied by $10^3$ for readability) from beauty hadron decays. A 1.9% normalization uncertainty not shown.

The pT-differential inclusive cross section of electrons from beauty hadron decays. A 1.9% normalization uncertainty not shown.

Inclusive bbbar production cross section per rapidity unit measured in mid-rapidity.

Inclusive bbbar production cross section per rapidity unit measured.

Total bbbar production cross section.

Total bbbar production cross section.

The difference in $v_{2}$ between protons and antiprotons as a function of $\sqrt{s_{NN}}$ for 0%-10%, 10%-40% and 40%-80% central Au + Au collisions. The systematic errors are shown by the hooked error bars. The dashed lines in the plot are fits with a power-law function.

No description provided.

The relative difference. The systematic errors are shown by the hooked error bars. The dashed lines in the plot are fits with a power-law function.

No description provided.

The $v_{2}$ difference between protons and antiprotons (and between $\pi^{+}$ and $pi^{-}$) for 10%-40% centrality Au+Au collisions at 7.7, 11.5, 14.5, and 19.6 GeV. The $v_{2}{BBC} results were slightly shifted horizontally.

No description provided.

No description provided.

No description provided.

No description provided.

No description provided.

No description provided.

No description provided.

No description provided.

No description provided.

No description provided.

No description provided.

No description provided.

No description provided.

No description provided.

No description provided.

No description provided.

No description provided.

No description provided.

No description provided.

No description provided.

No description provided.

No description provided.

No description provided.

No description provided.

No description provided.

No description provided.

No description provided.

No description provided.

No description provided.

No description provided.

No description provided.

No description provided.

No description provided.

No description provided.

No description provided.

No description provided.

No description provided.

No description provided.

No description provided.

No description provided.

No description provided.

No description provided.

No description provided.

No description provided.

No description provided.

No description provided.

No description provided.

No description provided.

No description provided.

No description provided.

No description provided.

No description provided.

No description provided.

No description provided.

No description provided.

No description provided.

No description provided.

No description provided.

No description provided.

No description provided.

No description provided.

No description provided.

No description provided.

No description provided.

No description provided.

No description provided.

No description provided.

No description provided.

No description provided.

No description provided.

No description provided.

No description provided.

No description provided.

No description provided.

No description provided.

No description provided.

No description provided.

No description provided.

No description provided.

No description provided.

No description provided.

No description provided.

No description provided.

No description provided.

No description provided.

No description provided.

No description provided.

No description provided.

No description provided.

No description provided.

No description provided.

No description provided.

No description provided.

No description provided.

No description provided.

No description provided.

No description provided.

No description provided.

No description provided.

No description provided.

No description provided.

No description provided.

No description provided.

No description provided.

No description provided.

No description provided.

No description provided.

No description provided.

No description provided.

No description provided.

No description provided.

No description provided.

No description provided.

No description provided.

No description provided.

No description provided.

No description provided.

No description provided.