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.

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.

The contribution from different systematic uncertainties to the measured $t\bar{t}t\bar{t}$ production cross section, grouped in categories.

The results of the fitted signal strength $\mu$ in the 1L/2LOS and 2LSS/3L combined channel.

The results of fitted inclusive ${t\bar{t}t\bar{t}}$ cross-section in the 1L/2LOS and 2LSS/3L combined channel.

Comparison between data and prediction for the distribution of the sum of the pseudo-continuous b-tagging score over the six jets with the highest score in the 1L,$\geq$9j,$\geq$3b region before the fit.

Comparison between data and prediction for the distribution of the sum of the pseudo-continuous b-tagging score over the six jets with the highest score in the 1L,$\geq$9j,$\geq$3b region after the fit.

Comparison between data and prediction for the distribution of the sum of the pseudo-continuous b-tagging score over the six jets with the highest score in the 2LOS,$\geq$7j,$\geq$3b region before the fit.

Comparison between data and prediction for the distribution of the sum of the pseudo-continuous b-tagging score over the six jets with the highest score in the 2LOS,$\geq$7j,$\geq$3b region after the fit.

Comparison between data and prediction for the distribution of the scalar sum of all jet and lepton pT in the event in the 1L,$\geq$8j,$\geq$3b region before the fit.

Comparison between data and prediction for the distribution of the scalar sum of all jet and lepton pT in the event in the 1L,$\geq$8j,$\geq$3b region after the fit.

Comparison between data and prediction for the distribution of the scalar sum of all jet and lepton pT in the event in the 2LOS,$\geq$6j,$\geq$3b region before the fit.

Comparison between data and prediction for the distribution of the scalar sum of all jet and lepton pT in the event in the 2LOS,$\geq$6j,$\geq$3b region after the fit.

Comparison between data and prediction for the distribution of number of jets in the 1L,$\geq$8j,$\geq$3b region before the fit.

Comparison between data and prediction for the distribution of number of jets in the 1L,$\geq$8j,$\geq$3b region after the fit.

Comparison between data and prediction for the distribution of number of jets in the 2LOS,$\geq$6j,$\geq$3b region before the fit.

Comparison between data and prediction for the distribution of number of jets in the 2LOS,$\geq$6j,$\geq$3b region after the fit.

Comparison between data and prediction for the distribution of b-jets multiplicity in the 1L,$\geq$8j,$\geq$3b region before the fit.

Comparison between data and prediction for the distribution of b-jets multiplicity in the 1L,$\geq$8j,$\geq$3b region after the fit.

Comparison between data and prediction for the distribution of b-jets multiplicity in the 2LOS,$\geq$6j,$\geq$3b region before the fit.

Comparison between data and prediction for the distribution of b-jets multiplicity in the 2LOS,$\geq$6j,$\geq$3b region after the fit.

Comparison between data and prediction for the distribution of the BDT score in the 1L,9j,4b signal region before the fit.

Comparison between data and prediction for the distribution of the BDT score in the 1L,9j,4b signal region after the fit.

Comparison between data and prediction for the distribution of the BDT score in the 1L,9j,$\geq$5b signal region before the fit.

Comparison between data and prediction for the distribution of the BDT score in the 1L,9j,$\geq$5b signal region after the fit.

Comparison between data and prediction for the distribution of the BDT score in the 1L,$\geq$10j,3bL signal region before the fit.

Comparison between data and prediction for the distribution of the BDT score in the 1L,$\geq$10j,3bL signal region after the fit.

Comparison between data and prediction for the distribution of the BDT score in the 1L,$\geq$10j,3bH signal region before the fit.

Comparison between data and prediction for the distribution of the BDT score in the 1L,$\geq$10j,3bH signal region after the fit.

Comparison between data and prediction for the distribution of the BDT score in the 1L,$\geq$10j,4b signal region before the fit.

Comparison between data and prediction for the distribution of the BDT score in the 1L,$\geq$10j,4b signal region after the fit.

Comparison between data and prediction for the distribution of the BDT score in the 1L,$\geq$10j,$\geq$5b signal region before the fit.

Comparison between data and prediction for the distribution of the BDT score in the 1L,$\geq$10j,$\geq$5b signal region after the fit.

Comparison between data and prediction for the distribution of the BDT score in the 2LOS,7j,$\geq$4b signal region before the fit.

Comparison between data and prediction for the distribution of the BDT score in the 2LOS,7j,$\geq$4b signal region after the fit.

Comparison between data and prediction for the distribution of the BDT score in the 2LOS,$\geq$8j,3bL signal region before the fit.

Comparison between data and prediction for the distribution of the BDT score in the 2LOS,$\geq$8j,3bL signal region after the fit.

Comparison between data and prediction for the distribution of the BDT score in the 2LOS,$\geq$8j,3bH signal region before the fit.

Comparison between data and prediction for the distribution of the BDT score in the 2LOS,$\geq$8j,3bH signal region after the fit.

Comparison between data and prediction for the distribution of the BDT score in the 2LOS,$\geq$8j,$\geq$4b signal region before the fit.

Comparison between data and prediction for the distribution of the BDT score in the 2LOS,$\geq$8j,$\geq$4b signal region after the fit.

Comparison between data and prediction for the distribution of the scalar sum of all jet and lepton pT in the event in the 1L,8j,3bV validation region before the fit.

Comparison between data and prediction for the distribution of the scalar sum of all jet and lepton pT in the event in the 1L,8j,3bV validation region after the fit.

Comparison between data and prediction for the distribution of the BDT score in the 1L,9j,3bV validation region before the fit.

Comparison between data and prediction for the distribution of the BDT score in the 1L,9j,3bV validation region after the fit.

Comparison between data and prediction for the distribution of the BDT score in the 1L,$\geq$10j,3bV validation region before the fit.

Comparison between data and prediction for the distribution of the BDT score in the 1L,$\geq$10j,3bV validation region after the fit.

Comparison between data and prediction for the distribution of the scalar sum of all jet and lepton pT in the event in the 2LOS,6j,3bV validation region before the fit.

Comparison between data and prediction for the distribution of the scalar sum of all jet and lepton pT in the event in the 2LOS,6j,3bV validation region after the fit.

Comparison between data and prediction for the distribution of the BDT score in the 2LOS,7j,3bV validation region before the fit.

Comparison between data and prediction for the distribution of the BDT score in the 2LOS,7j,3bV validation region after the fit.

Comparison between data and prediction for the distribution of the BDT score in the 2LOS,$\geq$8j,3bV validation region before the fit.

Comparison between data and prediction for the distribution of the BDT score in the 2LOS,$\geq$8j,3bV validation region after the fit.

Comparison between data and prediction for the distribution of the scalar sum of all jet and lepton pT in the event in the 1L,8j,3bL control region before the fit.

Comparison between data and prediction for the distribution of the scalar sum of all jet and lepton pT in the event in the 1L,8j,3bL control region after the fit.

Comparison between data and prediction for the distribution of the scalar sum of all jet and lepton pT in the event in the 1L,8j,3bH control region before the fit.

Comparison between data and prediction for the distribution of the scalar sum of all jet and lepton pT in the event in the 1L,8j,3bH control region after the fit.

Comparison between data and prediction for the distribution of the scalar sum of all jet and lepton pT in the event in the 1L,9j,3bL control region before the fit.

Comparison between data and prediction for the distribution of the scalar sum of all jet and lepton pT in the event in the 1L,9j,3bL control region after the fit.

Comparison between data and prediction for the distribution of the scalar sum of all jet and lepton pT in the event in the 1L,9j,3bH control region before the fit.

Comparison between data and prediction for the distribution of the scalar sum of all jet and lepton pT in the event in the 1L,9j,3bH control region after the fit.

Comparison between data and prediction for the distribution of the scalar sum of all jet and lepton pT in the event in the 1L,8j,4b control region before the fit.

Comparison between data and prediction for the distribution of the scalar sum of all jet and lepton pT in the event in the 1L,8j,4b control region after the fit.

Comparison between data and prediction for the distribution of the scalar sum of all jet and lepton pT in the event in the 1L,8j,$\geq$5b control region before the fit.

Comparison between data and prediction for the distribution of the scalar sum of all jet and lepton pT in the event in the 1L,8j,$\geq$5b control region after the fit.

Comparison between data and prediction for the distribution of the scalar sum of all jet and lepton pT in the event in the 2LOS,6j,3bL control region before the fit.

Comparison between data and prediction for the distribution of the scalar sum of all jet and lepton pT in the event in the 2LOS,6j,3bL control region after the fit.

Comparison between data and prediction for the distribution of the scalar sum of all jet and lepton pT in the event in the 2LOS,6j,3bH control region before the fit.

Comparison between data and prediction for the distribution of the scalar sum of all jet and lepton pT in the event in the 2LOS,6j,3bH control region after the fit.

Comparison between data and prediction for the distribution of the scalar sum of all jet and lepton pT in the event in the 2LOS,7j,3bL control region before the fit.

Comparison between data and prediction for the distribution of the scalar sum of all jet and lepton pT in the event in the 2LOS,7j,3bL control region after the fit.

Comparison between data and prediction for the distribution of the scalar sum of all jet and lepton pT in the event in the 2LOS,7j,3bH control region before the fit.

Comparison between data and prediction for the distribution of the scalar sum of all jet and lepton pT in the event in the 2LOS,7j,3bH control region after the fit.

Comparison between data and prediction for the distribution of the scalar sum of all jet and lepton pT in the event in the 2LOS,6j,$\geq$4b control region before the fit.

Comparison between data and prediction for the distribution of the scalar sum of all jet and lepton pT in the event in the 2LOS,6j,$\geq$4b control region after the fit.

Measured and predicted fiducial cross-sections in both the lllljj and ll$\nu\nu$jj channels for the inclusive ZZjj processes. Uncertainties due to different sources are presented.

Measured and predicted fiducial cross-sections in both the lllljj and ll$\nu\nu$jj channels for the inclusive ZZjj processes. Uncertainties due to different sources are presented.

Observed and expected multivariate discriminant distribution in the $\ell\ell\ell\ell jj$ QCD CR.

Observed and expected multivariate discriminant distribution in the $\ell\ell\ell\ell jj$ QCD CR.

Observed and expected multivariate discriminant distribution in the $\ell\ell\ell\ell jj$ SR.

Observed and expected multivariate discriminant distribution in the $\ell\ell\ell\ell jj$ SR.

Observed and expected multivariate discriminant distribution in the $\ell\ell\nu\nu jj$ SR.

Observed and expected multivariate discriminant distribution in the $\ell\ell\nu\nu jj$ SR.

The total inelastic cross section.

The measured differential elastic cross section. In addition to the statistical and total systematic uncertainties, the following 22 systematic shifts are given, which are included in the profile fit with their signs: -- Constraints: Beam optics uncertainty obtained by varying the ALFA constraints in the optics fit -- QScan: Variation by +/- 0.1 % of the quadrupole strength -- Q2: Fit of the strength of Q2 using the best value for the strength of Q1 and Q3 -- Q5Q6: Variation of the strength of Q5 and Q6 by -0.2% as indicated by machine constraints -- MadX: Uncertainty related to the beam transport replacing matrix transport by MadX PTC tracking -- Qmisal: Uncertainty due to the mis-alignment of the quadrupoles in the beam line -- Q1Q3: Propagation of the optics fit uncertainty in the strenght of Q1 and Q3 on the differential elastic cross section -- Aopt: Alignment uncertainty from the optimization procedure -- Offv: Alignment uncertainty related to the vertical beam center offset -- Offh: Alignment uncertainty related to the horizontal beam center offset -- Ang: Alignment uncertainty related to the detector rotation in the x-y plane -- BGn: Uncertainty from the background normalization -- BGs: Uncertainty from the background shape -- MCres: Error from modelling of the detector response -- Slope: Residual dependence on the physics model estimated by varying the nuclear slope in the simulation by +/- 1 GeV^-2 -- Emit: Uncertainty from the emittance used to calculate beam divergence in the simulation -- Unf: Unfolding uncertainty from the data-driven closure test -- Trac: Uncertainty from the variation of the track reconstruction selection cuts -- Xing: Uncertainty from residual crossing angle in the horizontal plane -- Eff: Uncertainty from the reconstruction efficiency -- Lumi: Luminosity uncertainty (+/- 1.5%) -- Ebeam: Uncertainty from the nominal beam energy (+/- 0.65%) Small differences in the values given here compared to the published version are related to insignificant rounding issues.

The statistical covariance matrix for the differential elastic cross section.

Cross-section upper limit in the $m(\tilde b_1)$-$m(\tilde\chi^0_1)$ plane for the best expected signal region.

Cross-section upper limit in the $m(\tilde b_1)$-$m(\tilde\chi^0_1)$ plane for the SRA250 signal region.

Cross-section upper limit in the $m(\tilde b_1)$-$m(\tilde\chi^0_1)$ plane for the SRA350 signal region.

Cross-section upper limit in the $m(\tilde b_1)$-$m(\tilde\chi^0_1)$ plane for the SRA450 signal region.

Cross-section upper limit in the $m(\tilde b_1)$-$m(\tilde\chi^0_1)$ plane for the SRB signal region.

Expected CLs values in the $m(\tilde b_1)$-$m(\tilde\chi^0_1)$ plane for the best expected signal region.

Expected CLs values in the $m(\tilde b_1)$-$m(\tilde\chi^0_1)$ plane for the SRA250 signal region.

Expected exclusion limit at 95% CL in the $m(\tilde b_1)$-$m(\tilde\chi^0_1)$ plane for the sbottom pair production scenario, for signal region SRA250.

Observed exclusion limit at 95% CL in the $m(\tilde b_1)$-$m(\tilde\chi^0_1)$ plane for the sbottom pair production scenario, for signal region SRA250.

Expected CLs values in the $m(\tilde b_1)$-$m(\tilde\chi^0_1)$ plane for the SRA350 signal region.

Expected exclusion limit at 95% CL in the $m(\tilde b_1)$-$m(\tilde\chi^0_1)$ plane for the sbottom pair production scenario, for signal region SRA350.

Observed exclusion limit at 95% CL in the $m(\tilde b_1)$-$m(\tilde\chi^0_1)$ plane for the sbottom pair production scenario, for signal region SRA350.

Expected CLs values in the $m(\tilde b_1)$-$m(\tilde\chi^0_1)$ plane for the SRA450 signal region.

Expected exclusion limit at 95% CL in the $m(\tilde b_1)$-$m(\tilde\chi^0_1)$ plane for the sbottom pair production scenario, for signal region SRA450.

Observed exclusion limit at 95% CL in the $m(\tilde b_1)$-$m(\tilde\chi^0_1)$ plane for the sbottom pair production scenario, for signal region SRA450.

Expected CLs values in the $m(\tilde b_1)$-$m(\tilde\chi^0_1)$ plane for the SRB signal region.

Expected exclusion limit at 95% CL in the $m(\tilde b_1)$-$m(\tilde\chi^0_1)$ plane for the sbottom pair production scenario, for signal region SRB.

Observed exclusion limit at 95% CL in the $m(\tilde b_1)$-$m(\tilde\chi^0_1)$ plane for the sbottom pair production scenario, for signal region SRB.

Observed CLs values in the $m(\tilde b_1)$-$m(\tilde\chi^0_1)$ plane for the best expected signal region.

Observed CLs values in the $m(\tilde b_1)$-$m(\tilde\chi^0_1)$ plane for the SRA250 signal region.

Observed CLs values in the $m(\tilde b_1)$-$m(\tilde\chi^0_1)$ plane for the SRA350 signal region.

Observed CLs values in the $m(\tilde b_1)$-$m(\tilde\chi^0_1)$ plane for the SRA450 signal region.

Observed CLs values in the $m(\tilde b_1)$-$m(\tilde\chi^0_1)$ plane for the SRB signal region.

Signal efficiency (in %) in the $m(\tilde b_1)$-$m(\tilde\chi^0_1)$ plane for the sbottom pair production scenario, for the best expected signal region.

Signal efficiency (in %) in the $m(\tilde b_1)$-$m(\tilde\chi^0_1)$ plane for the sbottom pair production scenario, for the SRA250 signal region.

Signal efficiency (in %) in the $m(\tilde b_1)$-$m(\tilde\chi^0_1)$ plane for the sbottom pair production scenario, for the SRA350 signal region.

Signal efficiency (in %) in the $m(\tilde b_1)$-$m(\tilde\chi^0_1)$ plane for the sbottom pair production scenario, for the SRA450 signal region.

Signal efficiency (in %) in the $m(\tilde b_1)$-$m(\tilde\chi^0_1)$ plane for the sbottom pair production scenario, for the SRB signal region.

Signal acceptance (in %) in the $m(\tilde b_1)$-$m(\tilde\chi^0_1)$ plane for the sbottom pair production scenario, for the best expected signal region.

Signal acceptance (in %) in the $m(\tilde b_1)$-$m(\tilde\chi^0_1)$ plane for the sbottom pair production scenario, for the SRA250 signal region.

Signal acceptance (in %) in the $m(\tilde b_1)$-$m(\tilde\chi^0_1)$ plane for the sbottom pair production scenario, for the SRA350 signal region.

Signal acceptance (in %) in the $m(\tilde b_1)$-$m(\tilde\chi^0_1)$ plane for the sbottom pair production scenario, for the SRA450 signal region.

Signal acceptance (in %) in the $m(\tilde b_1)$-$m(\tilde\chi^0_1)$ plane for the sbottom pair production scenario, for the SRB signal region.

Total experimental systematic uncertainty in percent on the signal efficiency times acceptance in the $m(\tilde b_1)$-$m(\tilde\chi^0_1)$ plane. The best expected signal region selection is used per point.

Data (bold boxes) and background estimates (colour fill) for m_beta vs. m_betagamma for the gluino R-hadron search (1000 GeV). The blue thin-line boxes illustrate the expected signal (on top of background) for the given R-hadron mass hypothesis. The black dashed vertical/horizontal lines at 500 GeV show the mass selection (signal region in the top-right). Two events pass this selection.

Expected (dashed black line) and observed (solid red line) 95% CL upper limits on the cross section as a function of mass for the production of long-lived gluino R-hadrons. The theory prediction along with its +-1sigma uncertainty is show as a black line and a blue band, respectively. The observed 8 TeV Run-1 limit and theory prediction [arXiv:1411.6795] are shown in dash-dotted and dotted lines, respectively.

Expected (dashed black line) and observed (solid red line) 95% CL upper limits on the cross section as a function of mass for the production of bottom-squark R-hadrons. The theory prediction along with its +-1sigma uncertainty is show as a black line and a blue band, respectively. The observed 8 TeV Run-1 limit and theory prediction [arXiv:1411.6795] are shown in dash-dotted and dotted lines, respectively.

Expected (dashed black line) and observed (solid red line) 95% CL upper limits on the cross section as a function of mass for the production of top-squark R-hadrons. The theory prediction along with its +-1sigma uncertainty is show as a black line and a blue band, respectively. The observed 8 TeV Run-1 limit and theory prediction [arXiv:1411.6795] are shown in dash-dotted and dotted lines, respectively.

Final selection requirements as a function of the simulated R-hadron mass.

Summary of all studied systematic uncertainties. Ranges indicate a dependency on the R-hadron mass hypothesis (from low to high masses).

Expected signal yield (Nsig) and efficiency (eff.), estimated background (Nbkg) and observed number of events in data (Nobs) for the full mass range after the final selection using 3.2/fb of data. The stated uncertainties include both the statistical and systematic contribution.

Distribution of the truth-level beta for gluino R-hadrons in exemplary signal MC samples and muons in a Zmumu MC sample. All distributions have been normalised to one. The last bin contains the overflow of the histograms. The distributions illustrate the good discriminating power of the variables.

Distribution of the truth-level betagamma for gluino R-hadrons in exemplary signal MC samples and muons in a Zmumu MC sample. All distributions have been normalised to one. The last bin contains the overflow of the histograms. The distributions illustrate the good discriminating power of the variables.

Expected (dashed black line) and observed (solid red line) 95% confidence level upper limits on the cross section as a function of mass for the production of long-lived gluino R-hadrons. The theory prediction along with its +-1sigma uncertainty is show as a black line and a blue band, respectively. For meta-stable gluinos with a lifetime of 50 ns. (mass exclusion: about 1660 GeV expected, 1520 GeV observed).

Expected (dashed black line) and observed (solid red line) 95% confidence level upper limits on the cross section as a function of mass for the production of long-lived gluino R-hadrons. The theory prediction along with its +-1sigma uncertainty is show as a black line and a blue band, respectively. For meta-stable gluinos with a lifetime of 30 ns. (mass exclusion: about 1660 GeV expected, 1520 GeV observed).

Expected (dashed black line) and observed (solid red line) 95% confidence level upper limits on the cross section as a function of mass for the production of long-lived gluino R-hadrons. The theory prediction along with its +-1sigma uncertainty is show as a black line and a blue band, respectively. For meta-stable gluinos with a lifetime of 10 ns. (mass exclusion: about 1660 GeV expected, 1520 GeV observed).

Object-quality selection cut-flow with observed data and exemplary expected events (scaled to 3.2/fb for MC) in the gluino R-hadron search.

Object-quality selection cut-flow with observed data and exemplary expected events (scaled to 3.2/fb for MC) in the squark R-hadron search.

Expected signal yield (Nsig) and efficiency (eff.), estimated background (Nbkg) and observed number of events in data (Nobs) for the full mass range in the meta-stable gluino R-hadron search using 3.2/fb of data. The stated uncertainties include both the statistical and systematic contribution.

Comparison of data with estimated backgrounds in the $\tilde{E}_\text{T}^\text{miss}$ distribution with the TZCR1 event selection except for the requirement on $\tilde{E}_\text{T}^\text{miss}$. The variables $\tilde{E}_\text{T}^\text{miss}$ and $\tilde{m}_\text{T}$ are constructed in the same way as $E_\text{T}^\text{miss}$ and $m_\text{T}$ but treating the leading photon transverse momentum as invisible. The predicted backgrounds are scaled with normalization factors. The uncertainty band includes statistical and all experimental systematic uncertainties. The last bin includes overflow.

Comparison of data with estimated backgrounds in the $\tilde{m}_\text{T}$ distribution with the TZCR1 event selection except for the requirement on $\tilde{m}_\text{T}$. The variables $\tilde{E}_\text{T}^\text{miss}$ and $\tilde{m}_\text{T}$ are constructed in the same way as $E_\text{T}^\text{miss}$ and $m_\text{T}$ but treating the leading photon transverse momentum as invisible. The predicted backgrounds are scaled with normalization factors. The uncertainty band includes statistical and all experimental systematic uncertainties. The last bin includes overflow.

Comparison of the observed data ($n_\text{obs}$) with the predicted background ($n_\text{exp}$) in the validation and signal regions. The background predictions are obtained using the background-only fit configuration. The bottom panel shows the significance of the difference between data and predicted background, where the significance is based on the total uncertainty ($\sigma_\text{tot}$).

Jet multiplicity distributions for events where exactly two signal leptons are selected. No correction factors are included in the background normalizations. The uncertainty band includes statistical and all experimental systematic uncertainties. The last bin includes overflow.

Jet multiplicity distributions for events where exactly one lepton plus one $\tau$ candidate are selected. No correction factors are included in the background normalizations. The uncertainty band includes statistical and all experimental systematic uncertainties. The last bin includes overflow.

The $E_\text{T}^\text{miss}$ distribution in SR1. In the plot, the full event selection in the corresponding signal region is applied, except for the requirement on $E_\text{T}^\text{miss}$. The predicted backgrounds are scaled with normalization factors. The uncertainty band includes statistical and all experimental systematic uncertainties. The last bin contains the overflow. Benchmark signal models are overlaid for comparison. The benchmark models are specified by the gluino and stop masses, given in TeV in the table.

The $m_\text{T}$ distribution in SR1. In the plot, the full event selection in the corresponding signal region is applied, except for the requirement on $m_\text{T}$. The predicted backgrounds are scaled with normalization factors. The uncertainty band includes statistical and all experimental systematic uncertainties. The last bin contains the overflow. Benchmark signal models are overlaid for comparison. The benchmark models are specified by the gluino and stop masses, given in TeV in the table.

Expected (black dashed) 95% excluded regions in the plane of $m_{\tilde{g}}$ versus $m_{\tilde{t}_1}$ for gluino-mediated stop production.

Observed (red solid) 95% excluded regions in the plane of $m_{\tilde{g}}$ versus $m_{\tilde{t}_1}$ for gluino-mediated stop production.

Expected (black dashed) 95% excluded regions in the plane of $m_{\tilde{t}_1}$ versus $m_{\tilde{\chi}_1^0}$ for direct stop production.

Observed (red solid) 95% excluded regions in the plane of $m_{\tilde{t}_1}$ versus $m_{\tilde{\chi}_1^0}$ for direct stop production.

The expected upper limits on $T$ quark pair production times the squared branching ratio for $T \rightarrow tZ$ as a function of the $T$ quark mass.

The observed upper limits on $T$ quark pair production times the squared branching ratio for $T \rightarrow tZ$ as a function of the $T$ quark mass.

The expected limits on $T$ quarks as a function of the branching ratios $B\left(T \rightarrow bW\right)$ and $B\left(T \rightarrow tH\right)$ for a $T$ quark with a mass of 800 GeV. The $T$ is assumed to decay in three possible ways: $T \to tZ$, $T \to tH$, and $T \to bW$.

The observed limits on $T$ quarks as a function of the branching ratios $B\left(T \rightarrow bW\right)$ and $B\left(T \rightarrow tH\right)$ for a $T$ quark with a mass of 800 GeV. The $T$ is assumed to decay in three possible ways: $T \to tZ$, $T \to tH$, and $T \to bW$.

The $m_\text{T}$ distribution in the WVR2-tail validation region which has the same preselection and jet $p_\text{T}$ requirements as SR2.

The $am_\text{T2}$ distribution in the WVR2-tail validation region which has the same preselection and jet $p_\text{T}$ requirements as SR2.

Large-radius jet mass ($R=1.2$), decomposed into the number of small-radius jet constituents. The lower panel shows the ratio of the total data to the total prediction (summed over all jet multiplicities). Events are required to have one lepton, four jets with $p_\text{T}>80,50,40,40$ GeV, at least one $b$-tagged jet, $E_\text{T}^\text{miss}>200$ GeV, and $m_\text{T}>30$ GeV.

Distribution of $m_\text{T2}^\tau$ in data for a selection enriched in $t\bar{t}$ events with one hadronically decaying $\tau$. Events that have no hadronic $\tau$ candidate (that passes the Loose identification criteria, as well as other requirements) are not shown in the plot.

Upper limits on the model cross-section in units of pb for the gluino-mediated stop models.

Upper limits on the model cross-section in units of pb for the models with direct stop pair production.

Illustration of the best expected signal region per signal grid point for the gluino-mediated stop models. This mapping is used for the final combined exclusion limits.

Illustration of the best expected signal region per signal grid point for models with direct stop pair production. This mapping is used for the final combined exclusion limits.

Expected $CL_s$ values for the gluino-mediated stop models.

Observed $CL_s$ values for the gluino-mediated stop models.

Expected $CL_s$ values for the direct stop pair production models.

Observed $CL_s$ values for the direct stop pair production models.

Expected limit using SR1 for models with direct stop pair production and an unpolarized stop (and bino LSP).

Expected limit using SR1 for models with direct stop pair production with $\tilde{t}_1=\tilde{t}_L$ (and bino LSP).

Expected limit using SR1 for models with direct stop pair production with $\tilde{t}_1\sim\tilde{t}_R$ (and bino LSP).

Observed limit using SR1 for models with direct stop pair production and an unpolarized stop (and bino LSP).

Observed limit using SR1 for models with direct stop pair production with $\tilde{t}_1=\tilde{t}_L$ (and bino LSP).

Observed limit using SR1 for models with direct stop pair production with $\tilde{t}_1\sim\tilde{t}_R$ (and bino LSP).

Expected limit using SR2 for models with direct stop pair production and an unpolarized stop (and bino LSP).

Expected limit using SR2 for models with direct stop pair production with $\tilde{t}_1=\tilde{t}_L$ (and bino LSP).

Expected limit using SR2 for models with direct stop pair production with $\tilde{t}_1\sim\tilde{t}_R$ (and bino LSP).

Observed limit using SR2 for models with direct stop pair production and an unpolarized stop (and bino LSP).

Observed limit using SR2 for models with direct stop pair production with $\tilde{t}_1=\tilde{t}_L$ (and bino LSP).

Observed limit using SR2 for models with direct stop pair production with $\tilde{t}_1\sim\tilde{t}_R$ (and bino LSP).

Expected limit using SR1+SR2 (best expected) for models with direct stop pair production and an unpolarized stop (and bino LSP).

Expected limit using SR1+SR2 (best expected) for models with direct stop pair production with $\tilde{t}_1=\tilde{t}_L$ (and bino LSP).

Expected limit using SR1+SR2 (best expected) for models with direct stop pair production with $\tilde{t}_1\sim\tilde{t}_R$ (and bino LSP).

Observed limit using SR1+SR2 (best expected) for models with direct stop pair production and an unpolarized stop (and bino LSP).

Observed limit using SR1+SR2 (best expected) for models with direct stop pair production with $\tilde{t}_1=\tilde{t}_L$ (and bino LSP).

Observed limit using SR1+SR2 (best expected) for models with direct stop pair production with $\tilde{t}_1\sim\tilde{t}_R$ (and bino LSP).

Acceptance for SR1 in the gluino-mediated stop models. The acceptance is defined as the fraction of signal events that pass the analysis selection performed on generator-level objects, therefore emulating an ideal detector with perfect particle identification and no measurement resolution effects.

Acceptance for SR1 in the direct stop pair production. The acceptance is defined as the fraction of signal events that pass the analysis selection performed on generator-level objects, therefore emulating an ideal detector with perfect particle identification and no measurement resolution effects.

Acceptance for SR2 in the gluino-mediated stop models. The acceptance is defined as the fraction of signal events that pass the analysis selection performed on generator-level objects, therefore emulating an ideal detector with perfect particle identification and no measurement resolution effects.

Acceptance for SR2 in the direct stop pair production. The acceptance is defined as the fraction of signal events that pass the analysis selection performed on generator-level objects, therefore emulating an ideal detector with perfect particle identification and no measurement resolution effects.

Acceptance for SR3 in the gluino-mediated stop models. The acceptance is defined as the fraction of signal events that pass the analysis selection performed on generator-level objects, therefore emulating an ideal detector with perfect particle identification and no measurement resolution effects.

Acceptance for SR3 in the direct stop pair production. The acceptance is defined as the fraction of signal events that pass the analysis selection performed on generator-level objects, therefore emulating an ideal detector with perfect particle identification and no measurement resolution effects.

Efficiency for SR1 in the gluino-mediated stop models. The efficiency is the ratio between the expected signal rate calculated with simulated data passing all the reconstruction level cuts applied to reconstructed objects, and the signal rate for an ideal detector (with perfect particle identification and no measurement resolution effects).

Efficiency for SR1 in the direct stop pair production. The efficiency is the ratio between the expected signal rate calculated with simulated data passing all the reconstruction level cuts applied to reconstructed objects, and the signal rate for an ideal detector (with perfect particle identification and no measurement resolution effects).

Efficiency for SR2 in the gluino-mediated stop models. The efficiency is the ratio between the expected signal rate calculated with simulated data passing all the reconstruction level cuts applied to reconstructed objects, and the signal rate for an ideal detector (with perfect particle identification and no measurement resolution effects).

Efficiency for SR2 in the direct stop pair production. The efficiency is the ratio between the expected signal rate calculated with simulated data passing all the reconstruction level cuts applied to reconstructed objects, and the signal rate for an ideal detector (with perfect particle identification and no measurement resolution effects).

Efficiency for SR3 in the gluino-mediated stop models. The efficiency is the ratio between the expected signal rate calculated with simulated data passing all the reconstruction level cuts applied to reconstructed objects, and the signal rate for an ideal detector (with perfect particle identification and no measurement resolution effects).

Efficiency for SR3 in the direct stop pair production. The efficiency is the ratio between the expected signal rate calculated with simulated data passing all the reconstruction level cuts applied to reconstructed objects, and the signal rate for an ideal detector (with perfect particle identification and no measurement resolution effects).

Median expected and observed 95% CL upper limits on the cross-section times branching ratio (sigma*B) for W'_SSM production for the exclusive muon and electron channels, and for both channels combined.