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.

Correlation matrix for the individual sources of systematic uncertainty considered in the analysis, for the combined opposite sign (OS) and same sign (SS) binned-template profile likelihood fit to data. The OS or SS refers to the charge signs of the primary lepton and the soft muon. The gamma parameters are NPs used to describe the effect of the limited statistics of the samples.

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 observed exclusion contour at 95% CL as a function of the $m_{H^{\pm}}$ vs. tan$\beta$ and assuming $m_a$ = 150 $\mathrm{GeV}$, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.7$. Masses that are within the contours are excluded. Only signals simulating the tW+DM final states are considered in this contour.

The expected exclusion contour at 95% CL as a function of the $m_{H^{\pm}}$ vs. tan$\beta$ and assuming $m_a$ = 150 $\mathrm{GeV}$, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.7$. Masses that are within the contours are excluded. Only signals simulating the tW+DM final states are considered in this contour.

The observed exclusion contour at 95% CL as a function of the $m_{H^{\pm}}$ vs. tan$\beta$ and assuming $m_a$ = 250 $\mathrm{GeV}$, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.7$. Masses that are within the contours are excluded. Only signals simulating the tW+DM final states are considered in this contour.

The expected exclusion contour at 95% CL as a function of the $m_{H^{\pm}}$ vs. tan$\beta$ and assuming $m_a$ = 250 $\mathrm{GeV}$, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.7$. Masses that are within the contours are excluded. Only signals simulating the tW+DM final states are considered in this contour.

The observed exclusion contour at 95% CL as a function of the $m_a$ vs. $m_{H^{\pm}}$ and assuming tan$\beta$ = 1, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.7$. Masses that are within the contours are excluded. Signals simulating the tW+DM + tt+DM final states are considered in this contour.

The expected exclusion contour at 95% CL as a function of the $m_a$ vs. $m_{H^{\pm}}$ and assuming tan$\beta$ = 1, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.7$. Masses that are within the contours are excluded. Signals simulating the tW+DM + tt+DM final states are considered in this contour.

The observed exclusion contour at 95% CL as a function of the $m_{H^{\pm}}$ vs. tan$\beta$ and assuming $m_a$ = 150 $\mathrm{GeV}$, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.7$. Masses that are within the contours are excluded. Signals simulating the tW+DM + tt+DM final states are considered in this contour.

The expected exclusion contour at 95% CL as a function of the $m_{H^{\pm}}$ vs. tan$\beta$ and assuming $m_a$ = 150 $\mathrm{GeV}$, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.7$. Masses that are within the contours are excluded. Signals simulating the tW+DM + tt+DM final states are considered in this contour.

The observed exclusion contour at 95% CL as a function of the $m_{H^{\pm}}$ vs. tan$\beta$ and assuming $m_a$ = 250 $\mathrm{GeV}$, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.7$. Masses that are within the contours are excluded. Signals simulating the tW+DM + tt+DM final states are considered in this contour.

The expected exclusion contour at 95% CL as a function of the $m_{H^{\pm}}$ vs. tan$\beta$ and assuming $m_a$ = 250 $\mathrm{GeV}$, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.7$. Masses that are within the contours are excluded. Signals simulating the tW+DM + tt+DM final states are considered in this contour.

The observed exclusion contour at 95% CL as a function of the $m_a$ vs. $m_{H^{\pm}}$ and assuming tan$\beta$ = 1, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.35$. Masses that are within the contours are excluded. Only signals simulating the tW+DM final states are considered in this contour.

The expected exclusion contour at 95% CL as a function of the $m_a$ vs. $m_{H^{\pm}}$ and assuming tan$\beta$ = 1, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.35$. Masses that are within the contours are excluded. Only signals simulating the tW+DM final states are considered in this contour.

The observed exclusion contour at 95% CL as a function of the $m_{H^{\pm}}$ vs. tan$\beta$ and assuming $m_a$ = 150 $\mathrm{GeV}$, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.35$. Masses that are within the contours are excluded. Only signals simulating the tW+DM final states are considered in this contour.

The expected exclusion contour at 95% CL as a function of the $m_{H^{\pm}}$ vs. tan$\beta$ and assuming $m_a$ = 150 $\mathrm{GeV}$, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.35$. Masses that are within the contours are excluded. Only signals simulating the tW+DM final states are considered in this contour.

The observed exclusion contour at 95% CL as a function of the $m_{H^{\pm}}$ vs. tan$\beta$ and assuming $m_a$ = 250 $\mathrm{GeV}$, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.35$. Masses that are within the contours are excluded. Only signals simulating the tW+DM final states are considered in this contour.

The expected exclusion contour at 95% CL as a function of the $m_{H^{\pm}}$ vs. tan$\beta$ and assuming $m_a$ = 250 $\mathrm{GeV}$, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.35$. Masses that are within the contours are excluded. Only signals simulating the tW+DM final states are considered in this contour.

The observed exclusion contour at 95% CL as a function of the $m_a$ vs. $m_{H^{\pm}}$ and assuming tan$\beta$ = 1, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.35$. Masses that are within the contours are excluded. Signals simulating the tW+DM + tt+DM final states are considered in this contour.

The expected exclusion contour at 95% CL as a function of the $m_a$ vs. $m_{H^{\pm}}$ and assuming tan$\beta$ = 1, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.35$. Masses that are within the contours are excluded. Signals simulating the tW+DM + tt+DM final states are considered in this contour.

The observed exclusion contour at 95% CL as a function of the $m_{H^{\pm}}$ vs. tan$\beta$ and assuming $m_a$ = 150 $\mathrm{GeV}$, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.35$. Masses that are within the contours are excluded. Signals simulating the tW+DM + tt+DM final states are considered in this contour.

The expected exclusion contour at 95% CL as a function of the $m_{H^{\pm}}$ vs. tan$\beta$ and assuming $m_a$ = 150 $\mathrm{GeV}$, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.35$. Masses that are within the contours are excluded. Signals simulating the tW+DM + tt+DM final states are considered in this contour.

The observed exclusion contour at 95% CL as a function of the $m_{H^{\pm}}$ vs. tan$\beta$ and assuming $m_a$ = 250 $\mathrm{GeV}$, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.35$. Masses that are within the contours are excluded. Signals simulating the tW+DM + tt+DM final states are considered in this contour.

The expected exclusion contour at 95% CL as a function of the $m_{H^{\pm}}$ vs. tan$\beta$ and assuming $m_a$ = 250 $\mathrm{GeV}$, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.35$. Masses that are within the contours are excluded. Signals simulating the tW+DM + tt+DM final states are considered in this contour.

The observed exclusion contour at 95% CL as a function of the $m_a$ vs. $m_{H^{\pm}}$ and assuming tan$\beta$ = 1, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.7$. Masses that are within the contours are excluded. Only signals simulating the tW+DM final states are considered in this contour. These exclusion contours are derived using the 0L channel only.

The expected exclusion contour at 95% CL as a function of the $m_a$ vs. $m_{H^{\pm}}$ and assuming tan$\beta$ = 1, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.7$. Masses that are within the contours are excluded. Only signals simulating the tW+DM final states are considered in this contour. These exclusion contours are derived using the 0L channel only.

The observed exclusion contour at 95% CL as a function of the $m_{H^{\pm}}$ vs. tan$\beta$ and assuming $m_a$ = 150 $\mathrm{GeV}$, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.7$. Masses that are within the contours are excluded. Only signals simulating the tW+DM final states are considered in this contour. These exclusion contours are derived using the 0L channel only.

The expected exclusion contour at 95% CL as a function of the $m_{H^{\pm}}$ vs. tan$\beta$ and assuming $m_a$ = 150 $\mathrm{GeV}$, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.7$. Masses that are within the contours are excluded. Only signals simulating the tW+DM final states are considered in this contour. These exclusion contours are derived using the 0L channel only.

The observed exclusion contour at 95% CL as a function of the $m_{H^{\pm}}$ vs. tan$\beta$ and assuming $m_a$ = 250 $\mathrm{GeV}$, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.7$. Masses that are within the contours are excluded. Only signals simulating the tW+DM final states are considered in this contour. These exclusion contours are derived using the 0L channel only.

The expected exclusion contour at 95% CL as a function of the $m_{H^{\pm}}$ vs. tan$\beta$ and assuming $m_a$ = 250 $\mathrm{GeV}$, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.7$. Masses that are within the contours are excluded. Only signals simulating the tW+DM final states are considered in this contour. These exclusion contours are derived using the 0L channel only.

The observed exclusion contour at 95% CL as a function of the $m_a$ vs. $m_{H^{\pm}}$ and assuming tan$\beta$ = 1, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.7$. Masses that are within the contours are excluded. Signals simulating the tW+DM + tt+DM final states are considered in this contour. These exclusion contours are derived using the 0L channel only.

The expected exclusion contour at 95% CL as a function of the $m_a$ vs. $m_{H^{\pm}}$ and assuming tan$\beta$ = 1, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.7$. Masses that are within the contours are excluded. Signals simulating the tW+DM + tt+DM final states are considered in this contour. These exclusion contours are derived using the 0L channel only.

The observed exclusion contour at 95% CL as a function of the $m_{H^{\pm}}$ vs. tan$\beta$ and assuming $m_a$ = 150 $\mathrm{GeV}$, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.7$. Masses that are within the contours are excluded. Signals simulating the tW+DM + tt+DM final states are considered in this contour. These exclusion contours are derived using the 0L channel only.

The expected exclusion contour at 95% CL as a function of the $m_{H^{\pm}}$ vs. tan$\beta$ and assuming $m_a$ = 150 $\mathrm{GeV}$, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.7$. Masses that are within the contours are excluded. Signals simulating the tW+DM + tt+DM final states are considered in this contour. These exclusion contours are derived using the 0L channel only.

The observed exclusion contour at 95% CL as a function of the $m_{H^{\pm}}$ vs. tan$\beta$ and assuming $m_a$ = 250 $\mathrm{GeV}$, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.7$. Masses that are within the contours are excluded. Signals simulating the tW+DM + tt+DM final states are considered in this contour. These exclusion contours are derived using the 0L channel only.

The expected exclusion contour at 95% CL as a function of the $m_{H^{\pm}}$ vs. tan$\beta$ and assuming $m_a$ = 250 $\mathrm{GeV}$, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.7$. Masses that are within the contours are excluded. Signals simulating the tW+DM + tt+DM final states are considered in this contour. These exclusion contours are derived using the 0L channel only.

The observed exclusion contour at 95% CL as a function of the $m_a$ vs. $m_{H^{\pm}}$ and assuming tan$\beta$ = 1, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.35$. Masses that are within the contours are excluded. Only signals simulating the tW+DM final states are considered in this contour. These exclusion contours are derived using the 0L channel only.

The expected exclusion contour at 95% CL as a function of the $m_a$ vs. $m_{H^{\pm}}$ and assuming tan$\beta$ = 1, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.35$. Masses that are within the contours are excluded. Only signals simulating the tW+DM final states are considered in this contour. These exclusion contours are derived using the 0L channel only.

The observed exclusion contour at 95% CL as a function of the $m_{H^{\pm}}$ vs. tan$\beta$ and assuming $m_a$ = 150 $\mathrm{GeV}$, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.35$. Masses that are within the contours are excluded. Only signals simulating the tW+DM final states are considered in this contour. These exclusion contours are derived using the 0L channel only.

The expected exclusion contour at 95% CL as a function of the $m_{H^{\pm}}$ vs. tan$\beta$ and assuming $m_a$ = 150 $\mathrm{GeV}$, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.35$. Masses that are within the contours are excluded. Only signals simulating the tW+DM final states are considered in this contour. These exclusion contours are derived using the 0L channel only.

The observed exclusion contour at 95% CL as a function of the $m_{H^{\pm}}$ vs. tan$\beta$ and assuming $m_a$ = 250 $\mathrm{GeV}$, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.35$. Masses that are within the contours are excluded. Only signals simulating the tW+DM final states are considered in this contour. These exclusion contours are derived using the 0L channel only.

The expected exclusion contour at 95% CL as a function of the $m_{H^{\pm}}$ vs. tan$\beta$ and assuming $m_a$ = 250 $\mathrm{GeV}$, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.35$. Masses that are within the contours are excluded. Only signals simulating the tW+DM final states are considered in this contour. These exclusion contours are derived using the 0L channel only.

The observed exclusion contour at 95% CL as a function of the $m_a$ vs. $m_{H^{\pm}}$ and assuming tan$\beta$ = 1, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.35$. Masses that are within the contours are excluded. Signals simulating the tW+DM + tt+DM final states are considered in this contour. These exclusion contours are derived using the 0L channel only.

The expected exclusion contour at 95% CL as a function of the $m_a$ vs. $m_{H^{\pm}}$ and assuming tan$\beta$ = 1, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.35$. Masses that are within the contours are excluded. Signals simulating the tW+DM + tt+DM final states are considered in this contour. These exclusion contours are derived using the 0L channel only.

The observed exclusion contour at 95% CL as a function of the $m_{H^{\pm}}$ vs. tan$\beta$ and assuming $m_a$ = 150 $\mathrm{GeV}$, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.35$. Masses that are within the contours are excluded. Signals simulating the tW+DM + tt+DM final states are considered in this contour. These exclusion contours are derived using the 0L channel only.

The expected exclusion contour at 95% CL as a function of the $m_{H^{\pm}}$ vs. tan$\beta$ and assuming $m_a$ = 150 $\mathrm{GeV}$, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.35$. Masses that are within the contours are excluded. Signals simulating the tW+DM + tt+DM final states are considered in this contour. These exclusion contours are derived using the 0L channel only.

The observed exclusion contour at 95% CL as a function of the $m_{H^{\pm}}$ vs. tan$\beta$ and assuming $m_a$ = 250 $\mathrm{GeV}$, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.35$. Masses that are within the contours are excluded. Signals simulating the tW+DM + tt+DM final states are considered in this contour. These exclusion contours are derived using the 0L channel only.

The expected exclusion contour at 95% CL as a function of the $m_{H^{\pm}}$ vs. tan$\beta$ and assuming $m_a$ = 250 $\mathrm{GeV}$, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.35$. Masses that are within the contours are excluded. Signals simulating the tW+DM + tt+DM final states are considered in this contour. These exclusion contours are derived using the 0L channel only.

The observed exclusion contour at 95% CL as a function of the $m_a$ vs. $m_{H^{\pm}}$ and assuming tan$\beta$ = 1, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.7$. Masses that are within the contours are excluded. Only signals simulating the tW+DM final states are considered in this contour. These exclusion contours are derived using the 1L channel only.

The expected exclusion contour at 95% CL as a function of the $m_a$ vs. $m_{H^{\pm}}$ and assuming tan$\beta$ = 1, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.7$. Masses that are within the contours are excluded. Only signals simulating the tW+DM final states are considered in this contour. These exclusion contours are derived using the 1L channel only.

The observed exclusion contour at 95% CL as a function of the $m_{H^{\pm}}$ vs. tan$\beta$ and assuming $m_a$ = 150 $\mathrm{GeV}$, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.7$. Masses that are within the contours are excluded. Only signals simulating the tW+DM final states are considered in this contour. These exclusion contours are derived using the 1L channel only.

The expected exclusion contour at 95% CL as a function of the $m_{H^{\pm}}$ vs. tan$\beta$ and assuming $m_a$ = 150 $\mathrm{GeV}$, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.7$. Masses that are within the contours are excluded. Only signals simulating the tW+DM final states are considered in this contour. These exclusion contours are derived using the 1L channel only.

The observed exclusion contour at 95% CL as a function of the $m_{H^{\pm}}$ vs. tan$\beta$ and assuming $m_a$ = 250 $\mathrm{GeV}$, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.7$. Masses that are within the contours are excluded. Only signals simulating the tW+DM final states are considered in this contour. These exclusion contours are derived using the 1L channel only.

The expected exclusion contour at 95% CL as a function of the $m_{H^{\pm}}$ vs. tan$\beta$ and assuming $m_a$ = 250 $\mathrm{GeV}$, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.7$. Masses that are within the contours are excluded. Only signals simulating the tW+DM final states are considered in this contour. These exclusion contours are derived using the 1L channel only.

The observed exclusion contour at 95% CL as a function of the $m_a$ vs. $m_{H^{\pm}}$ and assuming tan$\beta$ = 1, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.7$. Masses that are within the contours are excluded. Signals simulating the tW+DM + tt+DM final states are considered in this contour. These exclusion contours are derived using the 1L channel only.

The expected exclusion contour at 95% CL as a function of the $m_a$ vs. $m_{H^{\pm}}$ and assuming tan$\beta$ = 1, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.7$. Masses that are within the contours are excluded. Signals simulating the tW+DM + tt+DM final states are considered in this contour. These exclusion contours are derived using the 1L channel only.

The observed exclusion contour at 95% CL as a function of the $m_{H^{\pm}}$ vs. tan$\beta$ and assuming $m_a$ = 150 $\mathrm{GeV}$, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.7$. Masses that are within the contours are excluded. Signals simulating the tW+DM + tt+DM final states are considered in this contour. These exclusion contours are derived using the 1L channel only.

The expected exclusion contour at 95% CL as a function of the $m_{H^{\pm}}$ vs. tan$\beta$ and assuming $m_a$ = 150 $\mathrm{GeV}$, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.7$. Masses that are within the contours are excluded. Signals simulating the tW+DM + tt+DM final states are considered in this contour. These exclusion contours are derived using the 1L channel only.

The observed exclusion contour at 95% CL as a function of the $m_{H^{\pm}}$ vs. tan$\beta$ and assuming $m_a$ = 250 $\mathrm{GeV}$, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.7$. Masses that are within the contours are excluded. Signals simulating the tW+DM + tt+DM final states are considered in this contour. These exclusion contours are derived using the 1L channel only.

The expected exclusion contour at 95% CL as a function of the $m_{H^{\pm}}$ vs. tan$\beta$ and assuming $m_a$ = 250 $\mathrm{GeV}$, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.7$. Masses that are within the contours are excluded. Signals simulating the tW+DM + tt+DM final states are considered in this contour. These exclusion contours are derived using the 1L channel only.

The observed exclusion contour at 95% CL as a function of the $m_a$ vs. $m_{H^{\pm}}$ and assuming tan$\beta$ = 1, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.35$. Masses that are within the contours are excluded. Only signals simulating the tW+DM final states are considered in this contour. These exclusion contours are derived using the 1L channel only.

The expected exclusion contour at 95% CL as a function of the $m_a$ vs. $m_{H^{\pm}}$ and assuming tan$\beta$ = 1, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.35$. Masses that are within the contours are excluded. Only signals simulating the tW+DM final states are considered in this contour. These exclusion contours are derived using the 1L channel only.

The observed exclusion contour at 95% CL as a function of the $m_{H^{\pm}}$ vs. tan$\beta$ and assuming $m_a$ = 150 $\mathrm{GeV}$, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.35$. Masses that are within the contours are excluded. Only signals simulating the tW+DM final states are considered in this contour. These exclusion contours are derived using the 1L channel only.

The expected exclusion contour at 95% CL as a function of the $m_{H^{\pm}}$ vs. tan$\beta$ and assuming $m_a$ = 150 $\mathrm{GeV}$, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.35$. Masses that are within the contours are excluded. Only signals simulating the tW+DM final states are considered in this contour. These exclusion contours are derived using the 1L channel only.

The expected exclusion contour at 95% CL as a function of the $m_{H^{\pm}}$ vs. tan$\beta$ and assuming $m_a$ = 250 $\mathrm{GeV}$, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.35$. Masses that are within the contours are excluded. Only signals simulating the tW+DM final states are considered in this contour. These exclusion contours are derived using the 1L channel only.

The observed exclusion contour at 95% CL as a function of the $m_a$ vs. $m_{H^{\pm}}$ and assuming tan$\beta$ = 1, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.35$. Masses that are within the contours are excluded. Signals simulating the tW+DM + tt+DM final states are considered in this contour. These exclusion contours are derived using the 1L channel only.

The expected exclusion contour at 95% CL as a function of the $m_a$ vs. $m_{H^{\pm}}$ and assuming tan$\beta$ = 1, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.35$. Masses that are within the contours are excluded. Signals simulating the tW+DM + tt+DM final states are considered in this contour. These exclusion contours are derived using the 1L channel only.

The observed exclusion contour at 95% CL as a function of the $m_{H^{\pm}}$ vs. tan$\beta$ and assuming $m_a$ = 150 $\mathrm{GeV}$, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.35$. Masses that are within the contours are excluded. Signals simulating the tW+DM + tt+DM final states are considered in this contour. These exclusion contours are derived using the 1L channel only.

The expected exclusion contour at 95% CL as a function of the $m_{H^{\pm}}$ vs. tan$\beta$ and assuming $m_a$ = 150 $\mathrm{GeV}$, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.35$. Masses that are within the contours are excluded. Signals simulating the tW+DM + tt+DM final states are considered in this contour. These exclusion contours are derived using the 1L channel only.

The observed exclusion contour at 95% CL as a function of the $m_{H^{\pm}}$ vs. tan$\beta$ and assuming $m_a$ = 250 $\mathrm{GeV}$, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.35$. Masses that are within the contours are excluded. Signals simulating the tW+DM + tt+DM final states are considered in this contour. These exclusion contours are derived using the 1L channel only.

The expected exclusion contour at 95% CL as a function of the $m_{H^{\pm}}$ vs. tan$\beta$ and assuming $m_a$ = 250 $\mathrm{GeV}$, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.35$. Masses that are within the contours are excluded. Signals simulating the tW+DM + tt+DM final states are considered in this contour. These exclusion contours are derived using the 1L channel only.

The observed exclusion contour at 95% CL as a function of the $m_a$ vs. $m_{H^{\pm}}$ and assuming tan$\beta$ = 1, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.7$. Masses that are within the contours are excluded. Only signals simulating the tW+DM final states are considered in this contour. These exclusion contours are derived using the 2L channel only.

The expected exclusion contour at 95% CL as a function of the $m_a$ vs. $m_{H^{\pm}}$ and assuming tan$\beta$ = 1, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.7$. Masses that are within the contours are excluded. Only signals simulating the tW+DM final states are considered in this contour. These exclusion contours are derived using the 2L channel only.

The observed exclusion contour at 95% CL as a function of the $m_{H^{\pm}}$ vs. tan$\beta$ and assuming $m_a$ = 150 $\mathrm{GeV}$, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.7$. Masses that are within the contours are excluded. Only signals simulating the tW+DM final states are considered in this contour. These exclusion contours are derived using the 2L channel only.

The expected exclusion contour at 95% CL as a function of the $m_{H^{\pm}}$ vs. tan$\beta$ and assuming $m_a$ = 150 $\mathrm{GeV}$, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.7$. Masses that are within the contours are excluded. Only signals simulating the tW+DM final states are considered in this contour. These exclusion contours are derived using the 2L channel only.

The observed exclusion contour at 95% CL as a function of the $m_{H^{\pm}}$ vs. tan$\beta$ and assuming $m_a$ = 250 $\mathrm{GeV}$, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.7$. Masses that are within the contours are excluded. Only signals simulating the tW+DM final states are considered in this contour. These exclusion contours are derived using the 2L channel only.

The expected exclusion contour at 95% CL as a function of the $m_{H^{\pm}}$ vs. tan$\beta$ and assuming $m_a$ = 250 $\mathrm{GeV}$, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.7$. Masses that are within the contours are excluded. Only signals simulating the tW+DM final states are considered in this contour. These exclusion contours are derived using the 2L channel only.

The observed exclusion contour at 95% CL as a function of the $m_a$ vs. $m_{H^{\pm}}$ and assuming tan$\beta$ = 1, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.7$. Masses that are within the contours are excluded. Signals simulating the tW+DM + tt+DM final states are considered in this contour. These exclusion contours are derived using the 2L channel only.

The expected exclusion contour at 95% CL as a function of the $m_a$ vs. $m_{H^{\pm}}$ and assuming tan$\beta$ = 1, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.7$. Masses that are within the contours are excluded. Signals simulating the tW+DM + tt+DM final states are considered in this contour. These exclusion contours are derived using the 2L channel only.

The observed exclusion contour at 95% CL as a function of the $m_{H^{\pm}}$ vs. tan$\beta$ and assuming $m_a$ = 150 $\mathrm{GeV}$, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.7$. Masses that are within the contours are excluded. Signals simulating the tW+DM + tt+DM final states are considered in this contour. These exclusion contours are derived using the 2L channel only.

The expected exclusion contour at 95% CL as a function of the $m_{H^{\pm}}$ vs. tan$\beta$ and assuming $m_a$ = 150 $\mathrm{GeV}$, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.7$. Masses that are within the contours are excluded. Signals simulating the tW+DM + tt+DM final states are considered in this contour. These exclusion contours are derived using the 2L channel only.

The observed exclusion contour at 95% CL as a function of the $m_{H^{\pm}}$ vs. tan$\beta$ and assuming $m_a$ = 250 $\mathrm{GeV}$, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.7$. Masses that are within the contours are excluded. Signals simulating the tW+DM + tt+DM final states are considered in this contour. These exclusion contours are derived using the 2L channel only.

The expected exclusion contour at 95% CL as a function of the $m_{H^{\pm}}$ vs. tan$\beta$ and assuming $m_a$ = 250 $\mathrm{GeV}$, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.7$. Masses that are within the contours are excluded. Signals simulating the tW+DM + tt+DM final states are considered in this contour. These exclusion contours are derived using the 2L channel only.

The expected exclusion contour at 95% CL as a function of the $m_a$ vs. $m_{H^{\pm}}$ and assuming tan$\beta$ = 1, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.35$. Masses that are within the contours are excluded. Only signals simulating the tW+DM final states are considered in this contour. These exclusion contours are derived using the 2L channel only.

The expected exclusion contour at 95% CL as a function of the $m_{H^{\pm}}$ vs. tan$\beta$ and assuming $m_a$ = 150 $\mathrm{GeV}$, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.35$. Masses that are within the contours are excluded. Only signals simulating the tW+DM final states are considered in this contour. These exclusion contours are derived using the 2L channel only.

The expected exclusion contour at 95% CL as a function of the $m_{H^{\pm}}$ vs. tan$\beta$ and assuming $m_a$ = 250 $\mathrm{GeV}$, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.35$. Masses that are within the contours are excluded. Only signals simulating the tW+DM final states are considered in this contour. These exclusion contours are derived using the 2L channel only.

The expected exclusion contour at 95% CL as a function of the $m_a$ vs. $m_{H^{\pm}}$ and assuming tan$\beta$ = 1, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.35$. Masses that are within the contours are excluded. Signals simulating the tW+DM + tt+DM final states are considered in this contour. These exclusion contours are derived using the 2L channel only.

The observed exclusion contour at 95% CL as a function of the $m_{H^{\pm}}$ vs. tan$\beta$ and assuming $m_a$ = 150 $\mathrm{GeV}$, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.35$. Masses that are within the contours are excluded. Signals simulating the tW+DM + tt+DM final states are considered in this contour. These exclusion contours are derived using the 2L channel only.

The expected exclusion contour at 95% CL as a function of the $m_{H^{\pm}}$ vs. tan$\beta$ and assuming $m_a$ = 150 $\mathrm{GeV}$, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.35$. Masses that are within the contours are excluded. Signals simulating the tW+DM + tt+DM final states are considered in this contour. These exclusion contours are derived using the 2L channel only.

The observed exclusion contour at 95% CL as a function of the $m_{H^{\pm}}$ vs. tan$\beta$ and assuming $m_a$ = 250 $\mathrm{GeV}$, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.35$. Masses that are within the contours are excluded. Signals simulating the tW+DM + tt+DM final states are considered in this contour. These exclusion contours are derived using the 2L channel only.

The expected exclusion contour at 95% CL as a function of the $m_{H^{\pm}}$ vs. tan$\beta$ and assuming $m_a$ = 250 $\mathrm{GeV}$, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.35$. Masses that are within the contours are excluded. Signals simulating the tW+DM + tt+DM final states are considered in this contour. These exclusion contours are derived using the 2L channel only.

Model dependent upper limit on the cross section for the $m_a$ vs. $ m_{H^{\pm}}$ signal grid assuming tan$\beta$ = 1, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.7$. Only signals simulating the tW+DM final states are considered. Upper limits with large $\mu_{\mathrm{sig}}$ for the observed limit are capped at 500.

Model dependent upper limit on the cross section for the $m_{H^{\pm}}$ vs. tan$\beta$ signal grid assuming $m_a$ = 150 GeV, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.7$. Only signals simulating the tW+DM final states are considered. Upper limits with large $\mu_{\mathrm{sig}}$ for the observed limit are capped at 500.

Model dependent upper limit on the cross section for the $m_{H^{\pm}}$ vs. tan$\beta$ signal grid assuming $m_a$ = 250 GeV, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.7$. Only signals simulating the tW+DM final states are considered. Upper limits with large $\mu_{\mathrm{sig}}$ for the observed limit are capped at 500.

Model dependent upper limit on the cross section for the $m_a$ vs. $ m_{H^{\pm}}$ signal grid assuming tan$\beta$ = 1, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.7$. Signals simulating the tW+DM + tt+DM final states are considered. Upper limits with large $\mu_{\mathrm{sig}}$ for the observed limit are capped at 500.

Model dependent upper limit on the cross section for the $m_{H^{\pm}}$ vs. tan$\beta$ signal grid assuming $m_a$ = 150 GeV, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.7$. Signals simulating the tW+DM + tt+DM final states are considered. Upper limits with large $\mu_{\mathrm{sig}}$ for the observed limit are capped at 500.

Model dependent upper limit on the cross section for the $m_{H^{\pm}}$ vs. tan$\beta$ signal grid assuming $m_a$ = 250 GeV, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.7$. Signals simulating the tW+DM + tt+DM final states are considered. Upper limits with large $\mu_{\mathrm{sig}}$ for the observed limit are capped at 500.

Model dependent upper limit on the cross section for the $m_a$ vs. $ m_{H^{\pm}}$ signal grid assuming tan$\beta$ = 1, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.35$. Only signals simulating the tW+DM final states are considered. Upper limits with large $\mu_{\mathrm{sig}}$ for the observed limit are capped at 500.

Model dependent upper limit on the cross section for the $m_{H^{\pm}}$ vs. tan$\beta$ signal grid assuming $m_a$ = 150 GeV, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.35$. Only signals simulating the tW+DM final states are considered. Upper limits with large $\mu_{\mathrm{sig}}$ for the observed limit are capped at 500.

Model dependent upper limit on the cross section for the $m_{H^{\pm}}$ vs. tan$\beta$ signal grid assuming $m_a$ = 250 GeV, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.35$. Only signals simulating the tW+DM final states are considered. Upper limits with large $\mu_{\mathrm{sig}}$ for the observed limit are capped at 500.

Model dependent upper limit on the cross section for the $m_a$ vs. $ m_{H^{\pm}}$ signal grid assuming tan$\beta$ = 1, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.35$. Signals simulating the tW+DM + tt+DM final states are considered. Upper limits with large $\mu_{\mathrm{sig}}$ for the observed limit are capped at 500.

Model dependent upper limit on the cross section for the $m_{H^{\pm}}$ vs. tan$\beta$ signal grid assuming $m_a$ = 150 GeV, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.35$. Signals simulating the tW+DM + tt+DM final states are considered. Upper limits with large $\mu_{\mathrm{sig}}$ for the observed limit are capped at 500.

Model dependent upper limit on the cross section for the $m_{H^{\pm}}$ vs. tan$\beta$ signal grid assuming $m_a$ = 250 GeV, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.35$. Signals simulating the tW+DM + tt+DM final states are considered. Upper limits with large $\mu_{\mathrm{sig}}$ for the observed limit are capped at 500.

Model dependent upper limit on the cross section for the $m_a$ vs. $ m_{H^{\pm}}$ signal grid assuming tan$\beta$ = 1, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.7$. Only signals simulating the tW+DM final states are considered. Upper limits with large $\mu_{\mathrm{sig}}$ for the observed limit are capped at 500.

Model dependent upper limit on the cross section for the $m_{H^{\pm}}$ vs. tan$\beta$ signal grid assuming $m_a$ = 150 GeV, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.7$. Only signals simulating the tW+DM final states are considered. Upper limits with large $\mu_{\mathrm{sig}}$ for the observed limit are capped at 500.

Model dependent upper limit on the cross section for the $m_{H^{\pm}}$ vs. tan$\beta$ signal grid assuming $m_a$ = 250 GeV, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.7$. Only signals simulating the tW+DM final states are considered. Upper limits with large $\mu_{\mathrm{sig}}$ for the observed limit are capped at 500.

Model dependent upper limit on the cross section for the $m_a$ vs. $ m_{H^{\pm}}$ signal grid assuming tan$\beta$ = 1, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.7$. Signals simulating the tW+DM + tt+DM final states are considered. Upper limits with large $\mu_{\mathrm{sig}}$ for the observed limit are capped at 500.

Model dependent upper limit on the cross section for the $m_{H^{\pm}}$ vs. tan$\beta$ signal grid assuming $m_a$ = 150 GeV, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.7$. Signals simulating the tW+DM + tt+DM final states are considered. Upper limits with large $\mu_{\mathrm{sig}}$ for the observed limit are capped at 500.

Model dependent upper limit on the cross section for the $m_{H^{\pm}}$ vs. tan$\beta$ signal grid assuming $m_a$ = 250 GeV, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.7$. Signals simulating the tW+DM + tt+DM final states are considered. Upper limits with large $\mu_{\mathrm{sig}}$ for the observed limit are capped at 500.

Model dependent upper limit on the cross section for the $m_a$ vs. $ m_{H^{\pm}}$ signal grid assuming tan$\beta$ = 1, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.35$. Only signals simulating the tW+DM final states are considered. Upper limits with large $\mu_{\mathrm{sig}}$ for the observed limit are capped at 500.

Model dependent upper limit on the cross section for the $m_{H^{\pm}}$ vs. tan$\beta$ signal grid assuming $m_a$ = 150 GeV, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.35$. Only signals simulating the tW+DM final states are considered. Upper limits with large $\mu_{\mathrm{sig}}$ for the observed limit are capped at 500.

Model dependent upper limit on the cross section for the $m_{H^{\pm}}$ vs. tan$\beta$ signal grid assuming $m_a$ = 250 GeV, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.35$. Only signals simulating the tW+DM final states are considered. Upper limits with large $\mu_{\mathrm{sig}}$ for the observed limit are capped at 500.

Model dependent upper limit on the cross section for the $m_a$ vs. $ m_{H^{\pm}}$ signal grid assuming tan$\beta$ = 1, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.35$. Signals simulating the tW+DM + tt+DM final states are considered. Upper limits with large $\mu_{\mathrm{sig}}$ for the observed limit are capped at 500.

Model dependent upper limit on the cross section for the $m_{H^{\pm}}$ vs. tan$\beta$ signal grid assuming $m_a$ = 150 GeV, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.35$. Signals simulating the tW+DM + tt+DM final states are considered. Upper limits with large $\mu_{\mathrm{sig}}$ for the observed limit are capped at 500.

Model dependent upper limit on the cross section for the $m_{H^{\pm}}$ vs. tan$\beta$ signal grid assuming $m_a$ = 250 GeV, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.35$. Signals simulating the tW+DM + tt+DM final states are considered. Upper limits with large $\mu_{\mathrm{sig}}$ for the observed limit are capped at 500.

Model dependent upper limit on the cross section for the $m_a$ vs. $ m_{H^{\pm}}$ signal grid assuming tan$\beta$ = 1, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.7$. Only signals simulating the tW+DM final states are considered. Upper limits with large $\mu_{\mathrm{sig}}$ for the observed limit are capped at 500.

Model dependent upper limit on the cross section for the $m_{H^{\pm}}$ vs. tan$\beta$ signal grid assuming $m_a$ = 150 GeV, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.7$. Only signals simulating the tW+DM final states are considered. Upper limits with large $\mu_{\mathrm{sig}}$ for the observed limit are capped at 500.

Model dependent upper limit on the cross section for the $m_{H^{\pm}}$ vs. tan$\beta$ signal grid assuming $m_a$ = 250 GeV, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.7$. Only signals simulating the tW+DM final states are considered. Upper limits with large $\mu_{\mathrm{sig}}$ for the observed limit are capped at 500.

Model dependent upper limit on the cross section for the $m_a$ vs. $ m_{H^{\pm}}$ signal grid assuming tan$\beta$ = 1, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.7$. Signals simulating the tW+DM + tt+DM final states are considered. Upper limits with large $\mu_{\mathrm{sig}}$ for the observed limit are capped at 500.

Model dependent upper limit on the cross section for the $m_{H^{\pm}}$ vs. tan$\beta$ signal grid assuming $m_a$ = 150 GeV, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.7$. Signals simulating the tW+DM + tt+DM final states are considered. Upper limits with large $\mu_{\mathrm{sig}}$ for the observed limit are capped at 500.

Model dependent upper limit on the cross section for the $m_{H^{\pm}}$ vs. tan$\beta$ signal grid assuming $m_a$ = 250 GeV, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.7$. Signals simulating the tW+DM + tt+DM final states are considered. Upper limits with large $\mu_{\mathrm{sig}}$ for the observed limit are capped at 500.

Model dependent upper limit on the cross section for the $m_a$ vs. $ m_{H^{\pm}}$ signal grid assuming tan$\beta$ = 1, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.35$. Only signals simulating the tW+DM final states are considered. Upper limits with large $\mu_{\mathrm{sig}}$ for the observed limit are capped at 500.

Model dependent upper limit on the cross section for the $m_{H^{\pm}}$ vs. tan$\beta$ signal grid assuming $m_a$ = 150 GeV, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.35$. Only signals simulating the tW+DM final states are considered. Upper limits with large $\mu_{\mathrm{sig}}$ for the observed limit are capped at 500.

Model dependent upper limit on the cross section for the $m_{H^{\pm}}$ vs. tan$\beta$ signal grid assuming $m_a$ = 250 GeV, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.35$. Only signals simulating the tW+DM final states are considered. Upper limits with large $\mu_{\mathrm{sig}}$ for the observed limit are capped at 500.

Model dependent upper limit on the cross section for the $m_a$ vs. $ m_{H^{\pm}}$ signal grid assuming tan$\beta$ = 1, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.35$. Signals simulating the tW+DM + tt+DM final states are considered. Upper limits with large $\mu_{\mathrm{sig}}$ for the observed limit are capped at 500.

Model dependent upper limit on the cross section for the $m_{H^{\pm}}$ vs. tan$\beta$ signal grid assuming $m_a$ = 150 GeV, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.35$. Signals simulating the tW+DM + tt+DM final states are considered. Upper limits with large $\mu_{\mathrm{sig}}$ for the observed limit are capped at 500.

Model dependent upper limit on the cross section for the $m_{H^{\pm}}$ vs. tan$\beta$ signal grid assuming $m_a$ = 250 GeV, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.35$. Signals simulating the tW+DM + tt+DM final states are considered. Upper limits with large $\mu_{\mathrm{sig}}$ for the observed limit are capped at 500.

The distributions of $m_{\mathrm{b1},\mathrm{W-tagged}}$ in the 0L inclusive signal region. For each bin yields for the data and total SM prediction are provided. The SM prediction is provided with the total uncertainty, including the MC statistical uncertainty, detector-related systematic uncertainties and theoretical uncertainties. The rightmost bin includes overflow events.

The distributions of $m_{\mathrm{T}}^{\mathrm{b,E_{\mathrm{T}^{\mathrm{miss}}}}}$ in the 0L inclusive signal region. For each bin yields for the data and total SM prediction are provided. The SM prediction is provided with the total uncertainty, including the MC statistical uncertainty, detector-related systematic uncertainties and theoretical uncertainties. The rightmost bin includes overflow events.

The distributions of $N_{\mathrm{W-tagged}}$ in the 0L inclusive signal region. For each bin yields for the data and total SM prediction are provided. The SM prediction is provided with the total uncertainty, including the MC statistical uncertainty, detector-related systematic uncertainties and theoretical uncertainties. The rightmost bin includes overflow events.

The distributions of $m_{\mathrm{b1},\mathrm{\cancel{b1}}}$ in the hadronic top inclusive signal region. For each bin yields for the data and total SM prediction are provided. The SM prediction is provided with the total uncertainty, including the MC statistical uncertainty, detector-related systematic uncertainties and theoretical uncertainties. The rightmost bin includes overflow events.

The distributions of $m_{\mathrm{b1},\mathrm{\cancel{b1}}}$ in the leptonic top inclusive signal region. For each bin yields for the data and total SM prediction are provided. The SM prediction is provided with the total uncertainty, including the MC statistical uncertainty, detector-related systematic uncertainties and theoretical uncertainties. The rightmost bin includes overflow events.

The distributions of $m_{\mathrm{b1},\mathrm{\cancel{b1}}}$ in the leptonic top inclusive signal region. For each bin yields for the data and total SM prediction are provided. The SM prediction is provided with the total uncertainty, including the MC statistical uncertainty, detector-related systematic uncertainties and theoretical uncertainties. The rightmost bin includes overflow events.

Cutflow for the reference point $(\it{m}_{\mathrm{H^{\pm}}}, \it{m}_{a}, tan\beta, sin\theta )=$ (500,100,1,0.7) , (800,150,20,0.7), (600,250,30,0.7), (1000,400,1,0.7) in 0L regions. Results are shown including all correction factors applied to simulation, and is normalised to 139 fb$^{-1}$.

Cutflow for the reference point $(\it{m}_{\mathrm{H^{\pm}}}, \it{m}_{a}, tan\beta, sin\theta )=$ (500,100,1,0.7) , (800,150,20,0.7), (600,250,30,0.7), (1000,400,1,0.7) in 1L leptonic top regions. Results are shown including all correction factors applied to simulation, and is normalised to 139 fb$^{-1}$.

Cutflow for the reference point $(\it{m}_{\mathrm{H^{\pm}}}, \it{m}_{a}, tan\beta, sin\theta )=$ (500,100,1,0.7) , (800,150,20,0.7), (600,250,30,0.7), (1000,400,1,0.7) in 1L hadronic top regions. Results are shown including all correction factors applied to simulation, and is normalised to 139 fb$^{-1}$.

Signal acceptance in the 0L region for 2HDM+a model DM signals on the plane defined by m$_a$--m$_{H^{\pm}}$ assuming tan$\beta$ = 1, m$_{\chi}$= 10 GeV and sin$\theta$ = 0.7. Please mind that the acceptance given in the table is multiplied by factor of $10^{3}$

Signal efficiency in the 0L region for 2HDM+a model DM signals on the plane defined by m$_a$--m$_{H^{\pm}}$ assuming tan$\beta$ = 1, m$_{\chi}$= 10 GeV and sin$\theta$ = 0.7. Please mind that the efficiency given in the table is multiplied by factor of $10^{2}$

Signal acceptance in the 0L region for 2HDM+a model DM signals on the plane defined by m$_{H^{\pm}}$--tan$\beta$ assuming m$_a$ = 150 GeV, m$_{\chi}$= 10 GeV and sin$\theta$ = 0.7. Please mind that the acceptance given in the table is multiplied by factor of $10^{3}$

Signal efficiency in the 0L region for 2HDM+a model DM signals on the plane defined by m$_{H^{\pm}}$--tan$\beta$ assuming m$_a$ = 150 GeV, m$_{\chi}$= 10 GeV and sin$\theta$ = 0.7. Please mind that the efficiency given in the table is multiplied by factor of $10^{2}$

Signal acceptance in the 0L region for 2HDM+a model DM signals on the plane defined by m$_{H^{\pm}}$--tan$\beta$ assuming m$_a$ = 250 GeV, m$_{\chi}$= 10 GeV and sin$\theta$ = 0.7. Please mind that the acceptance given in the table is multiplied by factor of $10^{3}$

Signal efficiency in the 0L region for 2HDM+a model DM signals on the plane defined by m$_{H^{\pm}}$--tan$\beta$ assuming m$_a$ = 250 GeV, m$_{\chi}$= 10 GeV and sin$\theta$ = 0.7. Please mind that the efficiency given in the table is multiplied by factor of $10^{2}$

Signal acceptance in the 1L region for 2HDM+a model DM signals on the plane defined by m$_a$--m$_{H^{\pm}}$ assuming tan$\beta$ = 1, m$_{\chi}$= 10 GeV and sin$\theta$ = 0.7. Please mind that the acceptance given in the table is multiplied by factor of $10^{3}$

Signal efficiency in the 1L region for 2HDM+a model DM signals on the plane defined by m$_a$--m$_{H^{\pm}}$ assuming tan$\beta$ = 1, m$_{\chi}$= 10 GeV and sin$\theta$ = 0.7. Please mind that the efficiency given in the table is multiplied by factor of $10^{2}$

Signal acceptance in the 1L region for 2HDM+a model DM signals on the plane defined by m$_{H^{\pm}}$--tan$\beta$ assuming m$_a$ = 150 GeV, m$_{\chi}$= 10 GeV and sin$\theta$ = 0.7. Please mind that the acceptance given in the table is multiplied by factor of $10^{3}$

Signal efficiency in the 1L region for 2HDM+a model DM signals on the plane defined by m$_{H^{\pm}}$--tan$\beta$ assuming m$_a$ = 150 GeV, m$_{\chi}$= 10 GeV and sin$\theta$ = 0.7. Please mind that the efficiency given in the table is multiplied by factor of $10^{2}$

Signal acceptance in the 1L region for 2HDM+a model DM signals on the plane defined by m$_{H^{\pm}}$--tan$\beta$ assuming m$_a$ = 250 GeV, m$_{\chi}$= 10 GeV and sin$\theta$ = 0.7. Please mind that the acceptance given in the table is multiplied by factor of $10^{3}$

Signal efficiency in the 1L region for 2HDM+a model DM signals on the plane defined by m$_{H^{\pm}}$--tan$\beta$ assuming m$_a$ = 250 GeV, m$_{\chi}$= 10 GeV and sin$\theta$ = 0.7. Please mind that the efficiency given in the table is multiplied by factor of $10^{2}$

Signal acceptance in the 0L region for 2HDM+a model DM signals on the plane defined by m$_a$--m$_{H^{\pm}}$ assuming tan$\beta$ = 1, m$_{\chi}$= 10 GeV and sin$\theta$ = 0.35. Please mind that the acceptance given in the table is multiplied by factor of $10^{3}$

Signal efficiency in the 0L region for 2HDM+a model DM signals on the plane defined by m$_a$--m$_{H^{\pm}}$ assuming tan$\beta$ = 1, m$_{\chi}$= 10 GeV and sin$\theta$ = 0.35. Please mind that the efficiency given in the table is multiplied by factor of $10^{2}$

Signal acceptance in the 0L region for 2HDM+a model DM signals on the plane defined by m$_{H^{\pm}}$--tan$\beta$ assuming m$_a$ = 150 GeV, m$_{\chi}$= 10 GeV and sin$\theta$ = 0.35. Please mind that the acceptance given in the table is multiplied by factor of $10^{3}$

Signal efficiency in the 0L region for 2HDM+a model DM signals on the plane defined by m$_{H^{\pm}}$--tan$\beta$ assuming m$_a$ = 150 GeV, m$_{\chi}$= 10 GeV and sin$\theta$ = 0.35. Please mind that the efficiency given in the table is multiplied by factor of $10^{2}$

Signal acceptance in the 0L region for 2HDM+a model DM signals on the plane defined by m$_{H^{\pm}}$--tan$\beta$ assuming m$_a$ = 250 GeV, m$_{\chi}$= 10 GeV and sin$\theta$ = 0.35. Please mind that the acceptance given in the table is multiplied by factor of $10^{3}$

Signal efficiency in the 0L region for 2HDM+a model DM signals on the plane defined by m$_{H^{\pm}}$--tan$\beta$ assuming m$_a$ = 250 GeV, m$_{\chi}$= 10 GeV and sin$\theta$ = 0.35. Please mind that the efficiency given in the table is multiplied by factor of $10^{2}$

Signal acceptance in the 1L region for 2HDM+a model DM signals on the plane defined by m$_a$--m$_{H^{\pm}}$ assuming tan$\beta$ = 1, m$_{\chi}$= 10 GeV and sin$\theta$ = 0.35. Please mind that the acceptance given in the table is multiplied by factor of $10^{3}$

Signal efficiency in the 1L region for 2HDM+a model DM signals on the plane defined by m$_a$--m$_{H^{\pm}}$ assuming tan$\beta$ = 1, m$_{\chi}$= 10 GeV and sin$\theta$ = 0.35. Please mind that the efficiency given in the table is multiplied by factor of $10^{2}$

Signal acceptance in the 1L region for 2HDM+a model DM signals on the plane defined by m$_{H^{\pm}}$--tan$\beta$ assuming m$_a$ = 150 GeV, m$_{\chi}$= 10 GeV and sin$\theta$ = 0.35. Please mind that the acceptance given in the table is multiplied by factor of $10^{3}$

Signal efficiency in the 1L region for 2HDM+a model DM signals on the plane defined by m$_{H^{\pm}}$--tan$\beta$ assuming m$_a$ = 150 GeV, m$_{\chi}$= 10 GeV and sin$\theta$ = 0.35. Please mind that the efficiency given in the table is multiplied by factor of $10^{2}$

Signal acceptance in the 1L region for 2HDM+a model DM signals on the plane defined by m$_{H^{\pm}}$--tan$\beta$ assuming m$_a$ = 250 GeV, m$_{\chi}$= 10 GeV and sin$\theta$ = 0.35. Please mind that the acceptance given in the table is multiplied by factor of $10^{3}$

Signal efficiency in the 1L region for 2HDM+a model DM signals on the plane defined by m$_{H^{\pm}}$--tan$\beta$ assuming m$_a$ = 250 GeV, m$_{\chi}$= 10 GeV and sin$\theta$ = 0.35. Please mind that the efficiency given in the table is multiplied by factor of $10^{2}$

The EW-Zy jj differential cross-section in the Signal Region as a function of the leading lepton pT. The lower panels show the ratios of the MC predictions to the data. The band around unfolded data represents the total uncertainty (including statistical uncertainty) and takes into account the correlations as obtained from the fit. The hatched area represents the uncertainty in the prediction. Events beyond the upper limit of the histogram are included in the last bin.

The EW-Zy jj differential cross-section in the Signal Region as a function of the leading photon pT. The lower panels show the ratios of the MC predictions to the data. The band around unfolded data represents the total uncertainty (including statistical uncertainty) and takes into account the correlations as obtained from the fit. The hatched area represents the uncertainty in the prediction. Events beyond the upper limit of the histogram are included in the last bin.

The EW-Zy jj differential cross-section in the Signal Region as a function of the leading jet pT. The lower panels show the ratios of the MC predictions to the data. The band around unfolded data represents the total uncertainty (including statistical uncertainty) and takes into account the correlations as obtained from the fit. The hatched area represents the uncertainty in the prediction. Events beyond the upper limit of the histogram are included in the last bin.

The EW-Zy jj differential cross-section in the Signal Region as a function of the Zy system pT. The lower panels show the ratios of the MC predictions to the data. The band around unfolded data represents the total uncertainty (including statistical uncertainty) and takes into account the correlations as obtained from the fit. The hatched area represents the uncertainty in the prediction. Events beyond the upper limit of the histogram are included in the last bin.

The EW-Zy jj differential cross-section in the Signal Region as a function of the dijet invariant mass. The lower panels show the ratios of the MC predictions to the data. The band around unfolded data represents the total uncertainty (including statistical uncertainty) and takes into account the correlations as obtained from the fit. The hatched area represents the uncertainty in the prediction. Events beyond the upper limit of the histogram are included in the last bin.

The EW-Zy jj differential cross-section in the Signal Region as a function of the dijet rapidity difference. The lower panels show the ratios of the MC predictions to the data. The band around unfolded data represents the total uncertainty (including statistical uncertainty) and takes into account the correlations as obtained from the fit. The hatched area represents the uncertainty in the prediction. Events beyond the upper limit of the histogram are included in the last bin.

The EW-Zy jj differential cross-section in the Signal Region as a function of the Zy and dijet azimuthal difference. The lower panels show the ratios of the MC predictions to the data. The band around unfolded data represents the total uncertainty (including statistical uncertainty) and takes into account the correlations as obtained from the fit. The hatched area represents the uncertainty in the prediction. Events beyond the upper limit of the histogram are included in the last bin.

The total Zy jj differential cross-section in the Extended SR as a function of the leading lepton pT. The lower panels show the ratios of the MC predictions to the data. The band around unfolded data represents the total uncertainty (including statistical uncertainty) and takes into account the correlations as obtained from the fit. The hatched area represents the uncertainty in the prediction. Events beyond the upper limit of the histogram are included in the last bin.

The total Zy jj differential cross-section in the Extended SR as a function of the leading photon pT. The lower panels show the ratios of the MC predictions to the data. The band around unfolded data represents the total uncertainty (including statistical uncertainty) and takes into account the correlations as obtained from the fit. The hatched area represents the uncertainty in the prediction. Events beyond the upper limit of the histogram are included in the last bin.

The total Zy jj differential cross-section in the Extended SR as a function of the leading jet pT. The lower panels show the ratios of the MC predictions to the data. The band around unfolded data represents the total uncertainty (including statistical uncertainty) and takes into account the correlations as obtained from the fit. The hatched area represents the uncertainty in the prediction. Events beyond the upper limit of the histogram are included in the last bin.

The total Zy jj differential cross-section in the Extended SR as a function of the Z boson pT. The lower panels show the ratios of the MC predictions to the data. The band around unfolded data represents the total uncertainty (including statistical uncertainty) and takes into account the correlations as obtained from the fit. The hatched area represents the uncertainty in the prediction. Events beyond the upper limit of the histogram are included in the last bin.

The total Zy jj differential cross-section in the Extended SR as a function of the Zy system pT. The lower panels show the ratios of the MC predictions to the data. The band around unfolded data represents the total uncertainty (including statistical uncertainty) and takes into account the correlations as obtained from the fit. The hatched area represents the uncertainty in the prediction. Events beyond the upper limit of the histogram are included in the last bin.

The total Zy jj differential cross-section in the Extended SR as a function of the dijet invariant mass. The lower panels show the ratios of the MC predictions to the data. The band around unfolded data represents the total uncertainty (including statistical uncertainty) and takes into account the correlations as obtained from the fit. The hatched area represents the uncertainty in the prediction. Events beyond the upper limit of the histogram are included in the last bin.

The total Zy jj differential cross-section in the Extended SR as a function of the dijet rapidity difference. The lower panels show the ratios of the MC predictions to the data. The band around unfolded data represents the total uncertainty (including statistical uncertainty) and takes into account the correlations as obtained from the fit. The hatched area represents the uncertainty in the prediction. Events beyond the upper limit of the histogram are included in the last bin.

The total Zy jj differential cross-section in the Extended SR as a function of the Zy and dijet azimuthal difference. The lower panels show the ratios of the MC predictions to the data. The band around unfolded data represents the total uncertainty (including statistical uncertainty) and takes into account the correlations as obtained from the fit. The hatched area represents the uncertainty in the prediction. Events beyond the upper limit of the histogram are included in the last bin.

The total Zy jj differential cross-section in the Extended SR as a function of the centrality of the system ζ(Zy). The lower panels show the ratios of the MC predictions to the data. The band around unfolded data represents the total uncertainty (including statistical uncertainty) and takes into account the correlations as obtained from the fit. The hatched area represents the uncertainty in the prediction. Events beyond the upper limit of the histogram are included in the last bin.

Expected and observed 95% CL for the 1SFH mu Majorana model.

Expected and observed 95% CL for the 2QDH NH Dirac model.

Expected and observed 95% CL for the 2QDH NH Majorana model.

Expected and observed 95% CL for the 2QDH IH Dirac model.

Expected and observed 95% CL for the 2QDH IH Majorana model.

Cutflow for six simulated signal channels showing the weighted number of expected events based on the single-flavour mixing model in the Majorana limit. Each column uses the generated signal sample with the mass hypothesis $m_N = 10$ GeV and proper decay length $c\tau_N = 10$ mm.

Cutflow for the six channels in data showing the number of events passing each successive signal selection for Majorana HNLs.

The event selection efficiency for each mass-lifetime point in all six studied channels. Shown is the fraction of the produced MC simulation events that pass all signal region selections. An entry of 0 indicates no events were selected.

The dominant signal uncertainty is due to differences in reconstruction of displaced vertices and tracks between data and MC. This is evaluated by comparing $K^{0}_{S} \rightarrow \pi^+\pi^-$ event yields in the VR and in MC produced with Pythia8.186 in bins of $p_\mathrm{T}$ and $r_\mathrm{DV}$. The data/MC ratio is normalized to the bin nearest the IP where the tracking and vertexing reconstruction algorithms are expected to be most robust. The symmetrized difference from 1.0 is applied to each signal vertex as a per-event systematic variation.

Expected and observed yields in the different analysis regions (prefit) for the 1SFH e Dirac model (10 GeV, 10mm).

Expected and observed yields in the different analysis regions (postfit) for the 1SFH e Dirac model (10 GeV, 10mm).

Expected and observed yields in the different analysis regions (prefit) for the 1SFH e Majorana model (10 GeV, 10mm).

Expected and observed yields in the different analysis regions (postfit) for the 1SFH e Majorana model (10 GeV, 10mm).

Expected and observed yields in the different analysis regions (prefit) for the 1SFH u Dirac model (10 GeV, 10mm).

Expected and observed yields in the different analysis regions (postfit) for the 1SFH u Dirac model (10 GeV, 10mm).

Expected and observed yields in the different analysis regions (prefit) for the 1SFH u Majorana model (10 GeV, 10mm).

Expected and observed yields in the different analysis regions (postfit) for the 1SFH u Majorana model (10 GeV, 10mm).

Expected and observed yields in the different analysis regions (prefit) for the 2QDH (NH) Dirac model (10 GeV, 10mm).

Expected and observed yields in the different analysis regions (postfit) for the 2QDH (NH) Dirac model (10 GeV, 10mm).

Expected and observed yields in the different analysis regions (prefit) for the 2QDH (NH) Majorana model (10 GeV, 10mm).

Expected and observed yields in the different analysis regions (postfit) for the 2QDH (NH) Majorana model (10 GeV, 10mm).

Expected and observed yields in the different analysis regions (prefit) for the 2QDH (IH) Dirac model (10 GeV, 10mm).

Expected and observed yields in the different analysis regions (postfit) for the 2QDH (IH) Dirac model (10 GeV, 10mm).

Expected and observed yields in the different analysis regions (prefit) for the 2QDH (IH) Majorana model (10 GeV, 10mm).

Expected and observed yields in the different analysis regions (postfit) for the 2QDH (IH) Majorana model (10 GeV, 10mm).

The total displaced vertexing efficiency as a function of $r_{DV}$ for the custom configuration used in this analysis. The definition of the secondary vertex efficiency can be found in defined in \cite{ATL-PHYS-PUB-2019-013}. The efficiency is shown for $\mu-\mu\mu$, $\mu-\mu e$ and $\mu-ee$ signals with $m_N=10$~GeV and $c\tau_N=10$~mm.

Measured fiducial differential cross sections for the jet multiplicity in the high-p$_{\mathrm{T}}$ region in Z($\to \ell^{+} \ell^{-}$) + high p$_{\mathrm{T}}$ jets events. The statistical, systematic, and luminosity uncertainties are given.

Measured fiducial differential cross sections for the $\Delta R_{Z,j}^{min}$ in the high-p$_{\mathrm{T}}$ region in Z($\to \ell^{+} \ell^{-}$) + high p$_{\mathrm{T}}$ jets events. The statistical, systematic, and luminosity uncertainties are given.

Measured fiducial differential cross sections for the $r_{Z,j}$ in the high-p$_{\mathrm{T}}$ region in Z($\to \ell^{+} \ell^{-}$) + high p$_{\mathrm{T}}$ jets events. The statistical, systematic, and luminosity uncertainties are given.

Measured fiducial differential cross sections for the $r_{Z,j}$ in the collinear region in Z($\to \ell^{+} \ell^{-}$) + high p$_{\mathrm{T}}$ jets events. The statistical, systematic, and luminosity uncertainties are given.

Measured fiducial differential cross sections for the $r_{Z,j}$ in the back-to-back region in Z($\to \ell^{+} \ell^{-}$) + high p$_{\mathrm{T}}$ jets events. The statistical, systematic, and luminosity uncertainties are given.

Measured fiducial differential cross sections for the jet multiplicity in the collinear region in Z($\to \ell^{+} \ell^{-}$) + high p$_{\mathrm{T}}$ jets events. The statistical, systematic, and luminosity uncertainties are given.

Measured fiducial differential cross sections for the jet multiplicity in the back-to-back region in Z($\to \ell^{+} \ell^{-}$) + high p$_{\mathrm{T}}$ jets events. The statistical, systematic, and luminosity uncertainties are given.

Measured fiducial differential cross sections for the H$_{\mathrm{T}}$ in Z($\to \ell^{+} \ell^{-}$) + high p$_{\mathrm{T}}$ jets events. The statistical, systematic, and luminosity uncertainties are given.

Measured fiducial differential cross sections for the $\Delta R_{Z,j}^{min}$ in the high-S$_{\mathrm{T}}$ region in Z($\to \ell^{+} \ell^{-}$) + high p$_{\mathrm{T}}$ jets events. The statistical, systematic, and luminosity uncertainties are given.

Measured fiducial differential cross sections for the jet multiplicity in the high-S$_{\mathrm{T}}$ region in Z($\to \ell^{+} \ell^{-}$) + high p$_{\mathrm{T}}$ jets events. The statistical, systematic, and luminosity uncertainties are given.

Systematic uncertainties for the Z boson p$_{\mathrm{T}}$ in Z($\to \ell^{+} \ell^{-}$) + high p$_{\mathrm{T}}$ jets events. The uncertainties are presented as a percentage of the measured cross-section for the upward variation of each source of uncertainty in each bin.

Systematic uncertainties for the leading jet p$_{\mathrm{T}}$ in Z($\to \ell^{+} \ell^{-}$) + high p$_{\mathrm{T}}$ jets events. The uncertainties are presented as a percentage of the measured cross-section for the upward variation of each source of uncertainty in each bin.

Systematic uncertainties for the jet multiplicity in Z($\to \ell^{+} \ell^{-}$) + high p$_{\mathrm{T}}$ jets events. The uncertainties are presented as a percentage of the measured cross-section for the upward variation of each source of uncertainty in each bin.

Systematic uncertainties for the jet multiplicity in the high-p$_{\mathrm{T}}$ region in Z($\to \ell^{+} \ell^{-}$) + high p$_{\mathrm{T}}$ jets events. The uncertainties are presented as a percentage of the measured cross-section for the upward variation of each source of uncertainty in each bin.

Systematic uncertainties for the $\Delta R_{Z,j}^{min}$ in the high-p$_{\mathrm{T}}$ region in Z($\to \ell^{+} \ell^{-}$) + high p$_{\mathrm{T}}$ jets events. The uncertainties are presented as a percentage of the measured cross-section for the upward variation of each source of uncertainty in each bin.

Systematic uncertainties for the $r_{Z,j}$ in the high-p$_{\mathrm{T}}$ region in Z($\to \ell^{+} \ell^{-}$) + high p$_{\mathrm{T}}$ jets events. The uncertainties are presented as a percentage of the measured cross-section for the upward variation of each source of uncertainty in each bin.

Systematic uncertainties for the $r_{Z,j}$ in the collinear region in Z($\to \ell^{+} \ell^{-}$) + high p$_{\mathrm{T}}$ jets events. The uncertainties are presented as a percentage of the measured cross-section for the upward variation of each source of uncertainty in each bin.

Systematic uncertainties for the $r_{Z,j}$ in the back-to-back region in Z($\to \ell^{+} \ell^{-}$) + high p$_{\mathrm{T}}$ jets events. The uncertainties are presented as a percentage of the measured cross-section for the upward variation of each source of uncertainty in each bin.

Systematic uncertainties for the jet multiplicity in the collinear region in Z($\to \ell^{+} \ell^{-}$) + high p$_{\mathrm{T}}$ jets events. The uncertainties are presented as a percentage of the measured cross-section for the upward variation of each source of uncertainty in each bin.

Systematic uncertainties for the jet multiplicity in the back-to-back region in Z($\to \ell^{+} \ell^{-}$) + high p$_{\mathrm{T}}$ jets events. The uncertainties are presented as a percentage of the measured cross-section for the upward variation of each source of uncertainty in each bin.

Systematic uncertainties for the H$_{\mathrm{T}}$ in Z($\to \ell^{+} \ell^{-}$) + high p$_{\mathrm{T}}$ jets events. The uncertainties are presented as a percentage of the measured cross-section for the upward variation of each source of uncertainty in each bin.

Systematic uncertainties for the $\Delta R_{Z,j}^{min}$ in the high-S$_{\mathrm{T}}$ region in Z($\to \ell^{+} \ell^{-}$) + high p$_{\mathrm{T}}$ jets events. The uncertainties are presented as a percentage of the measured cross-section for the upward variation of each source of uncertainty in each bin.

Systematic uncertainties for the jet multiplicity in the high-S$_{\mathrm{T}}$ region in Z($\to \ell^{+} \ell^{-}$) + high p$_{\mathrm{T}}$ jets events. The uncertainties are presented as a percentage of the measured cross-section for the upward variation of each source of uncertainty in each bin.

Correction scale factor from the cross section calculated with leptons at the Born-level to the cross section calculated with dressed leptons as a function of the Z boson p$_{\mathrm{T}}$ in Z($\to \ell^{+} \ell^{-}$) + high p$_{\mathrm{T}}$ jets events, averaging the electron and muon channels, derived with Sherpa2.2.11. The systematic uncertainty is obtained with an enveloppe around scale factors computed from Sherpa2.2.1 and MG5_aMC+Py8 CKKWL.

Correction scale factor from the cross section calculated with leptons at the Born-level to the cross section calculated with dressed leptons as a function of the leading jet p$_{\mathrm{T}}$ in Z($\to \ell^{+} \ell^{-}$) + high p$_{\mathrm{T}}$ jets events, averaging the electron and muon channels, derived with Sherpa2.2.11. The systematic uncertainty is obtained with an enveloppe around scale factors computed from Sherpa2.2.1 and MG5_aMC+Py8 CKKWL.

Correction scale factor from the cross section calculated with leptons at the Born-level to the cross section calculated with dressed leptons as a function of the jet multiplicity in Z($\to \ell^{+} \ell^{-}$) + high p$_{\mathrm{T}}$ jets events, averaging the electron and muon channels, derived with Sherpa2.2.11. The systematic uncertainty is obtained with an enveloppe around scale factors computed from Sherpa2.2.1 and MG5_aMC+Py8 CKKWL.

Correction scale factor from the cross section calculated with leptons at the Born-level to the cross section calculated with dressed leptons as a function of the jet multiplicity in the high-p$_{\mathrm{T}}$ region in Z($\to \ell^{+} \ell^{-}$) + high p$_{\mathrm{T}}$ jets events, averaging the electron and muon channels, derived with Sherpa2.2.11. The systematic uncertainty is obtained with an enveloppe around scale factors computed from Sherpa2.2.1 and MG5_aMC+Py8 CKKWL.

Correction scale factor from the cross section calculated with leptons at the Born-level to the cross section calculated with dressed leptons as a function of the $\Delta R_{Z,j}^{min}$ in the high-p$_{\mathrm{T}}$ region in Z($\to \ell^{+} \ell^{-}$) + high p$_{\mathrm{T}}$ jets events, averaging the electron and muon channels, derived with Sherpa2.2.11. The systematic uncertainty is obtained with an enveloppe around scale factors computed from Sherpa2.2.1 and MG5_aMC+Py8 CKKWL.

Correction scale factor from the cross section calculated with leptons at the Born-level to the cross section calculated with dressed leptons as a function of the $r_{Z,j}$ in the high-p$_{\mathrm{T}}$ region in Z($\to \ell^{+} \ell^{-}$) + high p$_{\mathrm{T}}$ jets events, averaging the electron and muon channels, derived with Sherpa2.2.11. The systematic uncertainty is obtained with an enveloppe around scale factors computed from Sherpa2.2.1 and MG5_aMC+Py8 CKKWL.

Correction scale factor from the cross section calculated with leptons at the Born-level to the cross section calculated with dressed leptons as a function of the $r_{Z,j}$ in the collinear region in Z($\to \ell^{+} \ell^{-}$) + high p$_{\mathrm{T}}$ jets events, averaging the electron and muon channels, derived with Sherpa2.2.11. The systematic uncertainty is obtained with an enveloppe around scale factors computed from Sherpa2.2.1 and MG5_aMC+Py8 CKKWL.

Correction scale factor from the cross section calculated with leptons at the Born-level to the cross section calculated with dressed leptons as a function of the $r_{Z,j}$ in the back-to-back region in Z($\to \ell^{+} \ell^{-}$) + high p$_{\mathrm{T}}$ jets events, averaging the electron and muon channels, derived with Sherpa2.2.11. The systematic uncertainty is obtained with an enveloppe around scale factors computed from Sherpa2.2.1 and MG5_aMC+Py8 CKKWL.

Correction scale factor from the cross section calculated with leptons at the Born-level to the cross section calculated with dressed leptons as a function of the jet multiplicity in the collinear region in Z($\to \ell^{+} \ell^{-}$) + high p$_{\mathrm{T}}$ jets events, averaging the electron and muon channels, derived with Sherpa2.2.11. The systematic uncertainty is obtained with an enveloppe around scale factors computed from Sherpa2.2.1 and MG5_aMC+Py8 CKKWL.

Correction scale factor from the cross section calculated with leptons at the Born-level to the cross section calculated with dressed leptons as a function of the jet multiplicity in the back-to-back region in Z($\to \ell^{+} \ell^{-}$) + high p$_{\mathrm{T}}$ jets events, averaging the electron and muon channels, derived with Sherpa2.2.11. The systematic uncertainty is obtained with an enveloppe around scale factors computed from Sherpa2.2.1 and MG5_aMC+Py8 CKKWL.

Correction scale factor from the cross section calculated with leptons at the Born-level to the cross section calculated with dressed leptons as a function of the H$_{\mathrm{T}}$ in Z($\to \ell^{+} \ell^{-}$) + high p$_{\mathrm{T}}$ jets events, averaging the electron and muon channels, derived with Sherpa2.2.11. The systematic uncertainty is obtained with an enveloppe around scale factors computed from Sherpa2.2.1 and MG5_aMC+Py8 CKKWL.

Correction scale factor from the cross section calculated with leptons at the Born-level to the cross section calculated with dressed leptons as a function of the $\Delta R_{Z,j}^{min}$ in the high-S$_{\mathrm{T}}$ region in Z($\to \ell^{+} \ell^{-}$) + high p$_{\mathrm{T}}$ jets events, averaging the electron and muon channels, derived with Sherpa2.2.11. The systematic uncertainty is obtained with an enveloppe around scale factors computed from Sherpa2.2.1 and MG5_aMC+Py8 CKKWL.

Correction scale factor from the cross section calculated with leptons at the Born-level to the cross section calculated with dressed leptons as a function of the jet multiplicity in the high-S$_{\mathrm{T}}$ region in Z($\to \ell^{+} \ell^{-}$) + high p$_{\mathrm{T}}$ jets events, averaging the electron and muon channels, derived with Sherpa2.2.11. The systematic uncertainty is obtained with an enveloppe around scale factors computed from Sherpa2.2.1 and MG5_aMC+Py8 CKKWL.

Correction scale factor from the cross section calculated with an overlap removal with jets of pT greater than 100 GeV to the cross section calculated with an overlap removal with jets of pT greater than 30 GeV as a function of the Z boson p$_{\mathrm{T}}$ in Z($\to \ell^{+} \ell^{-}$) + high p$_{\mathrm{T}}$ jets events, averaging the electron and muon channels, derived with Sherpa2.2.11. The systematic uncertainty is obtained with an enveloppe around scale factors computed from Sherpa2.2.1 and MG5_aMC+Py8 CKKWL.

Correction scale factor from the cross section calculated with an overlap removal with jets of pT greater than 100 GeV to the cross section calculated with an overlap removal with jets of pT greater than 30 GeV as a function of the leading jet p$_{\mathrm{T}}$ in Z($\to \ell^{+} \ell^{-}$) + high p$_{\mathrm{T}}$ jets events, averaging the electron and muon channels, derived with Sherpa2.2.11. The systematic uncertainty is obtained with an enveloppe around scale factors computed from Sherpa2.2.1 and MG5_aMC+Py8 CKKWL.

Correction scale factor from the cross section calculated with an overlap removal with jets of pT greater than 100 GeV to the cross section calculated with an overlap removal with jets of pT greater than 30 GeV as a function of the jet multiplicity in Z($\to \ell^{+} \ell^{-}$) + high p$_{\mathrm{T}}$ jets events, averaging the electron and muon channels, derived with Sherpa2.2.11. The systematic uncertainty is obtained with an enveloppe around scale factors computed from Sherpa2.2.1 and MG5_aMC+Py8 CKKWL.

Correction scale factor from the cross section calculated with an overlap removal with jets of pT greater than 100 GeV to the cross section calculated with an overlap removal with jets of pT greater than 30 GeV as a function of the jet multiplicity in the high-p$_{\mathrm{T}}$ region in Z($\to \ell^{+} \ell^{-}$) + high p$_{\mathrm{T}}$ jets events, averaging the electron and muon channels, derived with Sherpa2.2.11. The systematic uncertainty is obtained with an enveloppe around scale factors computed from Sherpa2.2.1 and MG5_aMC+Py8 CKKWL.

Correction scale factor from the cross section calculated with an overlap removal with jets of pT greater than 100 GeV to the cross section calculated with an overlap removal with jets of pT greater than 30 GeV as a function of the $\Delta R_{Z,j}^{min}$ in the high-p$_{\mathrm{T}}$ region in Z($\to \ell^{+} \ell^{-}$) + high p$_{\mathrm{T}}$ jets events, averaging the electron and muon channels, derived with Sherpa2.2.11. The systematic uncertainty is obtained with an enveloppe around scale factors computed from Sherpa2.2.1 and MG5_aMC+Py8 CKKWL.

Correction scale factor from the cross section calculated with an overlap removal with jets of pT greater than 100 GeV to the cross section calculated with an overlap removal with jets of pT greater than 30 GeV as a function of the $r_{Z,j}$ in the high-p$_{\mathrm{T}}$ region in Z($\to \ell^{+} \ell^{-}$) + high p$_{\mathrm{T}}$ jets events, averaging the electron and muon channels, derived with Sherpa2.2.11. The systematic uncertainty is obtained with an enveloppe around scale factors computed from Sherpa2.2.1 and MG5_aMC+Py8 CKKWL.

Correction scale factor from the cross section calculated with an overlap removal with jets of pT greater than 100 GeV to the cross section calculated with an overlap removal with jets of pT greater than 30 GeV as a function of the $r_{Z,j}$ in the collinear region in Z($\to \ell^{+} \ell^{-}$) + high p$_{\mathrm{T}}$ jets events, averaging the electron and muon channels, derived with Sherpa2.2.11. The systematic uncertainty is obtained with an enveloppe around scale factors computed from Sherpa2.2.1 and MG5_aMC+Py8 CKKWL.

Correction scale factor from the cross section calculated with an overlap removal with jets of pT greater than 100 GeV to the cross section calculated with an overlap removal with jets of pT greater than 30 GeV as a function of the $r_{Z,j}$ in the back-to-back region in Z($\to \ell^{+} \ell^{-}$) + high p$_{\mathrm{T}}$ jets events, averaging the electron and muon channels, derived with Sherpa2.2.11. The systematic uncertainty is obtained with an enveloppe around scale factors computed from Sherpa2.2.1 and MG5_aMC+Py8 CKKWL.

Correction scale factor from the cross section calculated with an overlap removal with jets of pT greater than 100 GeV to the cross section calculated with an overlap removal with jets of pT greater than 30 GeV as a function of the jet multiplicity in the collinear region in Z($\to \ell^{+} \ell^{-}$) + high p$_{\mathrm{T}}$ jets events, averaging the electron and muon channels, derived with Sherpa2.2.11. The systematic uncertainty is obtained with an enveloppe around scale factors computed from Sherpa2.2.1 and MG5_aMC+Py8 CKKWL.

Correction scale factor from the cross section calculated with an overlap removal with jets of pT greater than 100 GeV to the cross section calculated with an overlap removal with jets of pT greater than 30 GeV as a function of the jet multiplicity in the back-to-back region in Z($\to \ell^{+} \ell^{-}$) + high p$_{\mathrm{T}}$ jets events, averaging the electron and muon channels, derived with Sherpa2.2.11. The systematic uncertainty is obtained with an enveloppe around scale factors computed from Sherpa2.2.1 and MG5_aMC+Py8 CKKWL.

Correction scale factor from the cross section calculated with an overlap removal with jets of pT greater than 100 GeV to the cross section calculated with an overlap removal with jets of pT greater than 30 GeV as a function of the H$_{\mathrm{T}}$ in Z($\to \ell^{+} \ell^{-}$) + high p$_{\mathrm{T}}$ jets events, averaging the electron and muon channels, derived with Sherpa2.2.11. The systematic uncertainty is obtained with an enveloppe around scale factors computed from Sherpa2.2.1 and MG5_aMC+Py8 CKKWL.

Correction scale factor from the cross section calculated with an overlap removal with jets of pT greater than 100 GeV to the cross section calculated with an overlap removal with jets of pT greater than 30 GeV as a function of the $\Delta R_{Z,j}^{min}$ in the high-S$_{\mathrm{T}}$ region in Z($\to \ell^{+} \ell^{-}$) + high p$_{\mathrm{T}}$ jets events, averaging the electron and muon channels, derived with Sherpa2.2.11. The systematic uncertainty is obtained with an enveloppe around scale factors computed from Sherpa2.2.1 and MG5_aMC+Py8 CKKWL.

Correction scale factor from the cross section calculated with an overlap removal with jets of pT greater than 100 GeV to the cross section calculated with an overlap removal with jets of pT greater than 30 GeV as a function of the jet multiplicity in the high-S$_{\mathrm{T}}$ region in Z($\to \ell^{+} \ell^{-}$) + high p$_{\mathrm{T}}$ jets events, averaging the electron and muon channels, derived with Sherpa2.2.11. The systematic uncertainty is obtained with an enveloppe around scale factors computed from Sherpa2.2.1 and MG5_aMC+Py8 CKKWL.

Measured fiducial differential cross sections for the Z boson p$_{\mathrm{T}}$ in Z($\to \ell^{+} \ell^{-}$) + high p$_{\mathrm{T}}$ jets events, where the EW Zjj contribution is treated as signal and not subtracted as background. The statistical, systematic, and luminosity uncertainties are given.

Measured fiducial differential cross sections for the leading jet p$_{\mathrm{T}}$ in Z($\to \ell^{+} \ell^{-}$) + high p$_{\mathrm{T}}$ jets events, where the EW Zjj contribution is treated as signal and not subtracted as background. The statistical, systematic, and luminosity uncertainties are given.

Measured fiducial differential cross sections for the jet multiplicity in Z($\to \ell^{+} \ell^{-}$) + high p$_{\mathrm{T}}$ jets events, where the EW Zjj contribution is treated as signal and not subtracted as background. The statistical, systematic, and luminosity uncertainties are given.

Measured fiducial differential cross sections for the jet multiplicity in the high-p$_{\mathrm{T}}$ region in Z($\to \ell^{+} \ell^{-}$) + high p$_{\mathrm{T}}$ jets events, where the EW Zjj contribution is treated as signal and not subtracted as background. The statistical, systematic, and luminosity uncertainties are given.

Measured fiducial differential cross sections for the $\Delta R_{Z,j}^{min}$ in the high-p$_{\mathrm{T}}$ region in Z($\to \ell^{+} \ell^{-}$) + high p$_{\mathrm{T}}$ jets events, where the EW Zjj contribution is treated as signal and not subtracted as background. The statistical, systematic, and luminosity uncertainties are given.

Measured fiducial differential cross sections for the $r_{Z,j}$ in the high-p$_{\mathrm{T}}$ region in Z($\to \ell^{+} \ell^{-}$) + high p$_{\mathrm{T}}$ jets events, where the EW Zjj contribution is treated as signal and not subtracted as background. The statistical, systematic, and luminosity uncertainties are given.

Measured fiducial differential cross sections for the $r_{Z,j}$ in the collinear region in Z($\to \ell^{+} \ell^{-}$) + high p$_{\mathrm{T}}$ jets events, where the EW Zjj contribution is treated as signal and not subtracted as background. The statistical, systematic, and luminosity uncertainties are given.

Measured fiducial differential cross sections for the $r_{Z,j}$ in the back-to-back region in Z($\to \ell^{+} \ell^{-}$) + high p$_{\mathrm{T}}$ jets events, where the EW Zjj contribution is treated as signal and not subtracted as background. The statistical, systematic, and luminosity uncertainties are given.

Measured fiducial differential cross sections for the jet multiplicity in the collinear region in Z($\to \ell^{+} \ell^{-}$) + high p$_{\mathrm{T}}$ jets events, where the EW Zjj contribution is treated as signal and not subtracted as background. The statistical, systematic, and luminosity uncertainties are given.

Measured fiducial differential cross sections for the jet multiplicity in the back-to-back region in Z($\to \ell^{+} \ell^{-}$) + high p$_{\mathrm{T}}$ jets events, where the EW Zjj contribution is treated as signal and not subtracted as background. The statistical, systematic, and luminosity uncertainties are given.

Measured fiducial differential cross sections for the H$_{\mathrm{T}}$ in Z($\to \ell^{+} \ell^{-}$) + high p$_{\mathrm{T}}$ jets events, where the EW Zjj contribution is treated as signal and not subtracted as background. The statistical, systematic, and luminosity uncertainties are given.

Measured fiducial differential cross sections for the $\Delta R_{Z,j}^{min}$ in the high-S$_{\mathrm{T}}$ region in Z($\to \ell^{+} \ell^{-}$) + high p$_{\mathrm{T}}$ jets events, where the EW Zjj contribution is treated as signal and not subtracted as background. The statistical, systematic, and luminosity uncertainties are given.

Measured fiducial differential cross sections for the jet multiplicity in the high-S$_{\mathrm{T}}$ region in Z($\to \ell^{+} \ell^{-}$) + high p$_{\mathrm{T}}$ jets events, where the EW Zjj contribution is treated as signal and not subtracted as background. The statistical, systematic, and luminosity uncertainties are given.

Systematic uncertainties for the Z boson p$_{\mathrm{T}}$ in Z($\to \ell^{+} \ell^{-}$) + high p$_{\mathrm{T}}$ jets events, where the EW Zjj contribution is treated as signal and not subtracted as background. The uncertainties are presented as a percentage of the measured cross-section for the upward variation of each source of uncertainty in each bin.

Systematic uncertainties for the leading jet p$_{\mathrm{T}}$ in Z($\to \ell^{+} \ell^{-}$) + high p$_{\mathrm{T}}$ jets events, where the EW Zjj contribution is treated as signal and not subtracted as background. The uncertainties are presented as a percentage of the measured cross-section for the upward variation of each source of uncertainty in each bin.

Systematic uncertainties for the jet multiplicity in Z($\to \ell^{+} \ell^{-}$) + high p$_{\mathrm{T}}$ jets events, where the EW Zjj contribution is treated as signal and not subtracted as background. The uncertainties are presented as a percentage of the measured cross-section for the upward variation of each source of uncertainty in each bin.

Systematic uncertainties for the jet multiplicity in the high-p$_{\mathrm{T}}$ region in Z($\to \ell^{+} \ell^{-}$) + high p$_{\mathrm{T}}$ jets events, where the EW Zjj contribution is treated as signal and not subtracted as background. The uncertainties are presented as a percentage of the measured cross-section for the upward variation of each source of uncertainty in each bin.

Systematic uncertainties for the $\Delta R_{Z,j}^{min}$ in the high-p$_{\mathrm{T}}$ region in Z($\to \ell^{+} \ell^{-}$) + high p$_{\mathrm{T}}$ jets events, where the EW Zjj contribution is treated as signal and not subtracted as background. The uncertainties are presented as a percentage of the measured cross-section for the upward variation of each source of uncertainty in each bin.

Systematic uncertainties for the $r_{Z,j}$ in the high-p$_{\mathrm{T}}$ region in Z($\to \ell^{+} \ell^{-}$) + high p$_{\mathrm{T}}$ jets events, where the EW Zjj contribution is treated as signal and not subtracted as background. The uncertainties are presented as a percentage of the measured cross-section for the upward variation of each source of uncertainty in each bin.

Systematic uncertainties for the $r_{Z,j}$ in the collinear region in Z($\to \ell^{+} \ell^{-}$) + high p$_{\mathrm{T}}$ jets events, where the EW Zjj contribution is treated as signal and not subtracted as background. The uncertainties are presented as a percentage of the measured cross-section for the upward variation of each source of uncertainty in each bin.

Systematic uncertainties for the $r_{Z,j}$ in the back-to-back region in Z($\to \ell^{+} \ell^{-}$) + high p$_{\mathrm{T}}$ jets events, where the EW Zjj contribution is treated as signal and not subtracted as background. The uncertainties are presented as a percentage of the measured cross-section for the upward variation of each source of uncertainty in each bin.

Systematic uncertainties for the jet multiplicity in the collinear region in Z($\to \ell^{+} \ell^{-}$) + high p$_{\mathrm{T}}$ jets events, where the EW Zjj contribution is treated as signal and not subtracted as background. The uncertainties are presented as a percentage of the measured cross-section for the upward variation of each source of uncertainty in each bin.

Systematic uncertainties for the jet multiplicity in the back-to-back region in Z($\to \ell^{+} \ell^{-}$) + high p$_{\mathrm{T}}$ jets events, where the EW Zjj contribution is treated as signal and not subtracted as background. The uncertainties are presented as a percentage of the measured cross-section for the upward variation of each source of uncertainty in each bin.

Systematic uncertainties for the H$_{\mathrm{T}}$ in Z($\to \ell^{+} \ell^{-}$) + high p$_{\mathrm{T}}$ jets events, where the EW Zjj contribution is treated as signal and not subtracted as background. The uncertainties are presented as a percentage of the measured cross-section for the upward variation of each source of uncertainty in each bin.

Systematic uncertainties for the $\Delta R_{Z,j}^{min}$ in the high-S$_{\mathrm{T}}$ region in Z($\to \ell^{+} \ell^{-}$) + high p$_{\mathrm{T}}$ jets events, where the EW Zjj contribution is treated as signal and not subtracted as background. The uncertainties are presented as a percentage of the measured cross-section for the upward variation of each source of uncertainty in each bin.

Systematic uncertainties for the jet multiplicity in the high-S$_{\mathrm{T}}$ region in Z($\to \ell^{+} \ell^{-}$) + high p$_{\mathrm{T}}$ jets events, where the EW Zjj contribution is treated as signal and not subtracted as background. The uncertainties are presented as a percentage of the measured cross-section for the upward variation of each source of uncertainty in each bin.

The observed exclusion contours as a function of (m(med), m($\chi$)), with Vector mediator)

The expected exclusion contours as a function of (m(a), tan($\beta$)), with sin($\theta$) = 0.35)

The observed exclusion contours as a function of (m(a), tan($\beta$)), with sin($\theta$) = 0.35)

The expected exclusion contours as a function of (m(a), tan($\beta$)), with sin($\theta$) = 0.7)

The observed exclusion contours as a function of (m(a), tan($\beta$)), with sin($\theta$) = 0.7)

The expected exclusion contours as a function of (m(H), tan($\beta$)), with sin($\theta$) = 0.35)

The observed exclusion contours as a function of (m(H), tan($\beta$)), with sin($\theta$) = 0.35)

The expected exclusion contours as a function of (m(H), tan($\beta$)), with sin($\theta$) = 0.7)

The observed exclusion contours as a function of (m(H), tan($\beta$)), with sin($\theta$) = 0.7)

The expected exclusion contours as a function of (m(a), m(H)), with sin($\theta$) = 0.35)

The observed exclusion contours as a function of (m(a), m(H)), with sin($\theta$) = 0.35)

The expected exclusion contours as a function of (m(a), m(H)), with sin($\theta$) = 0.7)

The observed exclusion contours as a function of (m(a), m(H)), with sin($\theta$) = 0.7)

Expected lower limit on signal strength at 95% CL as a function of sin($\theta$), with m(a) = 200 GeV, m(H) = 600 GeV.

Observed lower limit on signal strength at 95% CL as a function of sin($\theta$), with m(a) = 200 GeV, m(H) = 600 GeV.

Expected lower limit on signal strength at 95% CL as a function of sin($\theta$), with m(a) = 350 GeV, m(H) = 1000 GeV.

Observed lower limit on signal strength at 95% CL as a function of sin($\theta$), with m(a) = 350 GeV, m(H) = 1000 GeV.

Observed lower limit on WIMP-nucleon cross section at 90% CL as a function of m(WIMP), assuming Higgs-portal scenario with Scalar WIMP.

Observed lower limit on WIMP-nucleon cross section at 90% CL as a function of m(WIMP), assuming Higgs-portal scenario with Majorana WIMP.

Observed lower limit on the spin-dependent WIMP–proton scattering cross-section.

Observed lower limit on the spin-independent WIMP–nucleon scattering cross-section.

Cutflow of unweighted and weighted events of $ZH$ signals in the electron channel.

Cutflow of unweighted and weighted events of $ZH$ signals in the muon channel.

Cutflow of unweighted and weighted events of DM signals (simplified DM axial-vector with m(med) = 900 GeV and m($\chi$) = 40 GeV, 2HDM+$a$ signal with m(A) = 600 GeV, m(a) = 400 GeV, tan($\beta$) = 1.0 and m($\chi$) = 10 GeV) in the electron channel.

Cutflow of unweighted and weighted events of DM signals (simplified DM axial-vector with m(med) = 900 GeV and m($\chi$) = 40 GeV, 2HDM+$a$ signal with m(A) = 600 GeV, m(a) = 400 GeV, tan($\beta$) = 1.0 and m($\chi$) = 10 GeV) in the muon channel.

Observed (points with error bars) and expected background (solid histograms) distributions for $E_{T}^{miss}$ in the signal region (a) SRL, (b) SRM and (c) SRH after the background-only fit applied to the CRs. The predicted signal distributions for the two models with a gluino mass of 2000 GeV and neutralino mass of 250 GeV (SRL), 1050 GeV (SRM) or 1950 GeV (SRH) are also shown for comparison. The uncertainties in the SM background are only statistical.

Observed and expected exclusion limit in the gluino-neutralino mass plane at 95% CL combined using the signal region with the best expected sensitivity at each point, for the full Run-2 dataset corresponding to an integrated luminosity of $139~\mathrm{fb}^{-1}$, for $\gamma/Z$ (a) and $\gamma/h$ (b) signal models. The black solid line corresponds to the expected limits at 95% CL, with the light (yellow) bands indicating the 1$\sigma$ exclusions due to experimental and background-theory uncertainties. The observed limits are indicated by medium (red) curves, the solid contour represents the nominal limit, and the dotted lines are obtained by varying the signal cross section by the theoretical scale and PDF uncertainties. For each point in the higgsino-bino parameter space, the labels indicate the best-expected signal region, where L, M and H mean SRL, SRM and SRH, respectively.

Observed and expected exclusion limit in the gluino-neutralino mass plane at 95% CL combined using the signal region with the best expected sensitivity at each point, for the full Run-2 dataset corresponding to an integrated luminosity of $139~\mathrm{fb}^{-1}$, for $\gamma/Z$ (a) and $\gamma/h$ (b) signal models. The black solid line corresponds to the expected limits at 95% CL, with the light (yellow) bands indicating the 1$\sigma$ exclusions due to experimental and background-theory uncertainties. The observed limits are indicated by medium (red) curves, the solid contour represents the nominal limit, and the dotted lines are obtained by varying the signal cross section by the theoretical scale and PDF uncertainties. For each point in the higgsino-bino parameter space, the labels indicate the best-expected signal region, where L, M and H mean SRL, SRM and SRH, respectively.

Acceptance (left) and efficiency (right) for the $\gamma/Z$ model signal grid for SRL (top), SRM (middle) and SRH (bottom).

Acceptance (left) and efficiency (right) for the $\gamma/Z$ model signal grid for SRL (top), SRM (middle) and SRH (bottom).

Acceptance (left) and efficiency (right) for the $\gamma/Z$ model signal grid for SRL (top), SRM (middle) and SRH (bottom).

Acceptance (left) and efficiency (right) for the $\gamma/Z$ model signal grid for SRL (top), SRM (middle) and SRH (bottom).

Acceptance (left) and efficiency (right) for the $\gamma/Z$ model signal grid for SRL (top), SRM (middle) and SRH (bottom).

Acceptance (left) and efficiency (right) for the $\gamma/Z$ model signal grid for SRL (top), SRM (middle) and SRH (bottom).

Acceptance (left) and efficiency (right) for the $\gamma/h$ model signal grid for SRL (top), SRM (middle) and SRH (bottom).

Acceptance (left) and efficiency (right) for the $\gamma/h$ model signal grid for SRL (top), SRM (middle) and SRH (bottom).

Acceptance (left) and efficiency (right) for the $\gamma/h$ model signal grid for SRL (top), SRM (middle) and SRH (bottom).

Acceptance (left) and efficiency (right) for the $\gamma/h$ model signal grid for SRL (top), SRM (middle) and SRH (bottom).

Acceptance (left) and efficiency (right) for the $\gamma/h$ model signal grid for SRL (top), SRM (middle) and SRH (bottom).

Acceptance (left) and efficiency (right) for the $\gamma/h$ model signal grid for SRL (top), SRM (middle) and SRH (bottom).

Cutflow for the SRL selection, for two relevant signal points for both $\gamma/Z$ and $\gamma/h$ models, where the gluinos have mass of 2000 GeV and the neutralinos have a mass of 250 GeV (10000 generated events). The numbers are normalized to a luminosity of 139 $fb^{-1}$.

Cutflow for the SRM selection, for two relevant signal points for both $\gamma/Z$ and $\gamma/h$ models, where the gluinos have mass of 2000 GeV and the neutralinos have a mass of 1050 GeV (10000 generated events). The numbers are normalized to a luminosity of 139 $fb^{-1}$.

Cutflow for the SRH selection, for two relevant signal points for both $\gamma/Z$ and $\gamma/h$ models, where the gluinos have mass of 2000 GeV and the neutralinos have a mass of 1950 GeV (10000 generated events). The numbers are normalized to a luminosity of 139 $fb^{-1}$.

Observed and expected exclusion limits in the gluino–neutralino mass plane at 95% CL for the full Run-2 dataset corresponding to an integrated luminosity of 139 fb−1 , for the (a) $\gamma/Z$ and (b) $\gamma/h$ signal models. They are obtained by combining limits from the signal region with the best expected sensitivity at each point. The dashed (black) line corresponds to the expected limits at 95% CL, with the light (yellow) band indicating the $\pm 1\sigma$ excursions due to experimental and background-theory uncertainties. The observed limits are indicated by medium (red) curves: the solid contour represents the nominal limit, and the dotted lines are obtained by varying the signal cross section by the theoretical scale and PDF uncertainties.

Observed and expected exclusion limits in the gluino–neutralino mass plane at 95% CL for the full Run-2 dataset corresponding to an integrated luminosity of 139 fb−1 , for the (a) $\gamma/Z$ and (b) $\gamma/h$ signal models. They are obtained by combining limits from the signal region with the best expected sensitivity at each point. The dashed (black) line corresponds to the expected limits at 95% CL, with the light (yellow) band indicating the $\pm 1\sigma$ excursions due to experimental and background-theory uncertainties. The observed limits are indicated by medium (red) curves: the solid contour represents the nominal limit, and the dotted lines are obtained by varying the signal cross section by the theoretical scale and PDF uncertainties.

Observed and expected exclusion limits in the gluino–neutralino mass plane at 95% CL for the full Run-2 dataset corresponding to an integrated luminosity of 139 fb−1 , for the (a) $\gamma/Z$ and (b) $\gamma/h$ signal models. They are obtained by combining limits from the signal region with the best expected sensitivity at each point. The dashed (black) line corresponds to the expected limits at 95% CL, with the light (yellow) band indicating the $\pm 1\sigma$ excursions due to experimental and background-theory uncertainties. The observed limits are indicated by medium (red) curves: the solid contour represents the nominal limit, and the dotted lines are obtained by varying the signal cross section by the theoretical scale and PDF uncertainties.

Observed and expected exclusion limits in the gluino–neutralino mass plane at 95% CL for the full Run-2 dataset corresponding to an integrated luminosity of 139 fb−1 , for the (a) $\gamma/Z$ and (b) $\gamma/h$ signal models. They are obtained by combining limits from the signal region with the best expected sensitivity at each point. The dashed (black) line corresponds to the expected limits at 95% CL, with the light (yellow) band indicating the $\pm 1\sigma$ excursions due to experimental and background-theory uncertainties. The observed limits are indicated by medium (red) curves: the solid contour represents the nominal limit, and the dotted lines are obtained by varying the signal cross section by the theoretical scale and PDF uncertainties.

Observed and expected exclusion limits in the gluino–neutralino mass plane at 95% CL for the full Run-2 dataset corresponding to an integrated luminosity of 139 fb−1 , for the (a) $\gamma/Z$ and (b) $\gamma/h$ signal models. They are obtained by combining limits from the signal region with the best expected sensitivity at each point. The dashed (black) line corresponds to the expected limits at 95% CL, with the light (yellow) band indicating the $\pm 1\sigma$ excursions due to experimental and background-theory uncertainties. The observed limits are indicated by medium (red) curves: the solid contour represents the nominal limit, and the dotted lines are obtained by varying the signal cross section by the theoretical scale and PDF uncertainties.

Observed and expected exclusion limits in the gluino–neutralino mass plane at 95% CL for the full Run-2 dataset corresponding to an integrated luminosity of 139 fb−1 , for the (a) $\gamma/Z$ and (b) $\gamma/h$ signal models. They are obtained by combining limits from the signal region with the best expected sensitivity at each point. The dashed (black) line corresponds to the expected limits at 95% CL, with the light (yellow) band indicating the $\pm 1\sigma$ excursions due to experimental and background-theory uncertainties. The observed limits are indicated by medium (red) curves: the solid contour represents the nominal limit, and the dotted lines are obtained by varying the signal cross section by the theoretical scale and PDF uncertainties.

Observed and expected exclusion limits in the gluino–neutralino mass plane at 95% CL for the full Run-2 dataset corresponding to an integrated luminosity of 139 fb−1 , for the (a) $\gamma/Z$ and (b) $\gamma/h$ signal models. They are obtained by combining limits from the signal region with the best expected sensitivity at each point. The dashed (black) line corresponds to the expected limits at 95% CL, with the light (yellow) band indicating the $\pm 1\sigma$ excursions due to experimental and background-theory uncertainties. The observed limits are indicated by medium (red) curves: the solid contour represents the nominal limit, and the dotted lines are obtained by varying the signal cross section by the theoretical scale and PDF uncertainties.

Observed and expected exclusion limits in the gluino–neutralino mass plane at 95% CL for the full Run-2 dataset corresponding to an integrated luminosity of 139 fb−1 , for the (a) $\gamma/Z$ and (b) $\gamma/h$ signal models. They are obtained by combining limits from the signal region with the best expected sensitivity at each point. The dashed (black) line corresponds to the expected limits at 95% CL, with the light (yellow) band indicating the $\pm 1\sigma$ excursions due to experimental and background-theory uncertainties. The observed limits are indicated by medium (red) curves: the solid contour represents the nominal limit, and the dotted lines are obtained by varying the signal cross section by the theoretical scale and PDF uncertainties.

Observed and expected exclusion limits in the gluino–neutralino mass plane at 95% CL for the full Run-2 dataset corresponding to an integrated luminosity of 139 fb−1 , for the (a) $\gamma/Z$ and (b) $\gamma/h$ signal models. They are obtained by combining limits from the signal region with the best expected sensitivity at each point. The dashed (black) line corresponds to the expected limits at 95% CL, with the light (yellow) band indicating the $\pm 1\sigma$ excursions due to experimental and background-theory uncertainties. The observed limits are indicated by medium (red) curves: the solid contour represents the nominal limit, and the dotted lines are obtained by varying the signal cross section by the theoretical scale and PDF uncertainties.

Observed and expected exclusion limits in the gluino–neutralino mass plane at 95% CL for the full Run-2 dataset corresponding to an integrated luminosity of 139 fb−1 , for the (a) $\gamma/Z$ and (b) $\gamma/h$ signal models. They are obtained by combining limits from the signal region with the best expected sensitivity at each point. The dashed (black) line corresponds to the expected limits at 95% CL, with the light (yellow) band indicating the $\pm 1\sigma$ excursions due to experimental and background-theory uncertainties. The observed limits are indicated by medium (red) curves: the solid contour represents the nominal limit, and the dotted lines are obtained by varying the signal cross section by the theoretical scale and PDF uncertainties.

Observed and expected exclusion limits in the gluino–neutralino mass plane at 95% CL for the full Run-2 dataset corresponding to an integrated luminosity of 139 fb−1 , for the (a) $\gamma/Z$ and (b) $\gamma/h$ signal models. They are obtained by combining limits from the signal region with the best expected sensitivity at each point. The dashed (black) line corresponds to the expected limits at 95% CL, with the light (yellow) band indicating the $\pm 1\sigma$ excursions due to experimental and background-theory uncertainties. The observed limits are indicated by medium (red) curves: the solid contour represents the nominal limit, and the dotted lines are obtained by varying the signal cross section by the theoretical scale and PDF uncertainties.

Observed and expected exclusion limits in the gluino–neutralino mass plane at 95% CL for the full Run-2 dataset corresponding to an integrated luminosity of 139 fb−1 , for the (a) $\gamma/Z$ and (b) $\gamma/h$ signal models. They are obtained by combining limits from the signal region with the best expected sensitivity at each point. The dashed (black) line corresponds to the expected limits at 95% CL, with the light (yellow) band indicating the $\pm 1\sigma$ excursions due to experimental and background-theory uncertainties. The observed limits are indicated by medium (red) curves: the solid contour represents the nominal limit, and the dotted lines are obtained by varying the signal cross section by the theoretical scale and PDF uncertainties.

The $1\sigma_{\text{exp}}$ variation of the expected exclusion contour at 95% CL as a function of the smuon and neutralino masses, for $\lambda_{231}^{'}=1.0$.

The observed exclusion contour at 95% CL as a function of the smuon and neutralino masses, for $\lambda_{231}^{'}=0.1$.

The expected exclusion contour at 95% CL as a function of the smuon and neutralino masses, for $\lambda_{231}^{'}=0.1$.

The observed exclusion contour at 95% CL as a function of the smuon and neutralino masses, for $\lambda_{231}^{'}=0.15$.

The expected exclusion contour at 95% CL as a function of the smuon and neutralino masses, for $\lambda_{231}^{'}=0.15$.

The observed exclusion contour at 95% CL as a function of the smuon and neutralino masses, for $\lambda_{231}^{'}=0.2$.

The expected exclusion contour at 95% CL as a function of the smuon and neutralino masses, for $\lambda_{231}^{'}=0.2$.

The observed exclusion contour at 95% CL as a function of the smuon and neutralino masses, for $\lambda_{231}^{'}=0.4$.

The expected exclusion contour at 95% CL as a function of the smuon and neutralino masses, for $\lambda_{231}^{'}=0.4$.

The observed exclusion contour at 95% CL as a function of the smuon and neutralino masses, for $\lambda_{231}^{'}=0.6$.

The expected exclusion contour at 95% CL as a function of the smuon and neutralino masses, for $\lambda_{231}^{'}=0p6$.

The observed exclusion contour at 95% CL as a function of the smuon and neutralino masses, for $\lambda_{231}^{'}=1.5$.

The expected exclusion contour at 95% CL as a function of the smuon and neutralino masses, for $\lambda_{231}^{'}=1.5$.

The observed exclusion contour at 95% CL as a function of the leptoquark mass and coupling strength.

The expected exclusion contour at 95% CL as a function of the leptoquark mass and coupling strength.

The minus $1\sigma_{\text{theory}}$ variation of the observed exclusion contour at 95% CL as a function of the leptoquark mass and coupling strength.

The plus $1\sigma_{\text{theory}}$ variation of the observed exclusion contour at 95% CL as a function of the leptoquark mass and coupling strength.

The $1\sigma_{\text{exp}}$ variation of the expected exclusion contour at 95% CL as a function of the leptoquark mass and coupling strength.

Observed yields, and fake lepton background yields in the $e^{+}\mu^{-}$ and $e^{-}\mu^{+}$ channels of SR-MET, along with the results of the $e^{+}\mu^{-}/e^{-}\mu^{+}$ ratio measurement and 1-sided p-value in SR-MET, binned in $M_{T2}$.

Observed yields, and fake lepton background yields in the $e^{+}\mu^{-}$ and $e^{-}\mu^{+}$ channels of SR-JET, along with the results of the $e^{+}\mu^{-}/e^{-}\mu^{+}$ ratio measurement and 1-sided p-value in SR-JET, binned in $H_{\text{P}}$.

Observed and expected 95% CL upper limits on the total number of signal events entering the $e^{+}\mu^{-}$ and $e^{-}\mu^{+}$ channels of each bin of SR-MET. The regions are binned in the same way as the ratio $\rho$ measurement. The limits are shown for a selection of 'z' values, where 'z' is the fraction of the total signal events entering the $e^{+}\mu^{-}$ channel.

Observed and expected 95% CL upper limits on the total number of signal events entering the $e^{+}\mu^{-}$ and $e^{-}\mu^{+}$ channels of each bin of SR-JET. The regions are binned in the same way as the ratio $\rho$ measurement. The limits are shown for a selection of 'z' values, where 'z' is the fraction of the total signal events entering the $e^{+}\mu^{-}$ channel.

Signal yields following each cut in the analysis, for representative $R$-parity-violating supersymmetry and leptoquark signals. All yields are MC generator-weighted and normalised to $139~\text{fb}^{-1}$. The cut labelled `Preselection' includes trigger requirements, and requires exactly one Baseline electron and one Baseline muon. At this point, the muon charge-bias correction weights are also applied. The $R$-parity-violating supersymmetry models were generated by specifying a top-quark in the final state and applying a two-lepton filter, hence the first row also includes events where the top quark decays to an electron.