Expected (dashed line) value of $-2\ln\Lambda$ as a function of $\kappa_{\lambda}$ for the combined fit, when all other coupling modifiers are fixed to their SM predictions.

Observed (solid line) value of $-2\ln\Lambda$ as a function of $\kappa_{2V}$ for the combined fit, when all other coupling modifiers are fixed to their SM predictions.

Expected (dashed line) value of $-2\ln\Lambda$ as a function of $\kappa_{2V}$ for the combined fit, when all other coupling modifiers are fixed to their SM predictions.

Likelihood observed contours at 68% CL (solid line) and 95% CL (dashed line) in the ($\kappa_{\lambda}$, $\kappa_{2V}$) parameter space, when all other coupling modifiers are fixed to their SM predictions. The best-fit value, denoted by a cross in the plot, corresponds to the very first entry in the table.

Likelihood expected contours at 68% CL (teal shaded region) and 95% CL (yellow shaded region) in the ($\kappa_{\lambda}$, $\kappa_{2V}$) parameter space, when all other coupling modifiers are fixed to their SM predictions. In the plot, the SM prediction is indicated by the star.

Likelihood contours at 68% (solid line) and 95% (dashed line) CL in the ($c_{hhh}$, $c_{gghh}$) parameter space, when all other Wilson coefficients are fixed to their SM values. The corresponding expected contours are shown by the inner and outer shaded regions. The SM prediction is indicated by the star, while the observed best-fit value is denoted by the black cross. The best-fit value corresponds to the very first entry in the table.

Likelihood contours at 68% (solid line) and 95% (dashed line) CL in the ($c_{hhh}$, $c_{gghh}$) parameter space, when all other Wilson coefficients are fixed to their SM values. The corresponding expected contours are shown by the inner and outer shaded regions. The SM prediction is indicated by the star, while the observed best-fit value is denoted by the black cross.

Likelihood contours at 68% (solid line) and 95% (dashed line) CL in the ($c_{hhh}$, $c_{tthh}$) parameter space, when all other Wilson coefficients are fixed to their SM values. The corresponding expected contours are shown by the inner and outer shaded regions. The SM prediction is indicated by the star, while the observed best-fit value is denoted by the black cross. The best-fit value corresponds to the very first entry in the table.

Likelihood contours at 68% (solid line) and 95% (dashed line) CL in the ($c_{hhh}$, $c_{tthh}$) parameter space, when all other Wilson coefficients are fixed to their SM values. The corresponding expected contours are shown by the inner and outer shaded regions. The SM prediction is indicated by the star, while the observed best-fit value is denoted by the black cross.

Likelihood contours at 68% (solid line) and 95% (dashed line) CL in the ($c_{H}$, $c_{H◻}$) parameter space, when all other Wilson coefficients are fixed to their SM values. The corresponding expected contours are shown by the inner and outer shaded regions. The SM prediction is indicated by the star, while the observed best-fit value is denoted by the black cross. The best-fit value corresponds to the very first entry in the table.

Likelihood contours at 68% (solid line) and 95% (dashed line) CL in the ($c_{H}$, $c_{H◻}$) parameter space, when all other Wilson coefficients are fixed to their SM values. The corresponding expected contours are shown by the inner and outer shaded regions. The SM prediction is indicated by the star, while the observed best-fit value is denoted by the black cross.

The acceptance times efficiency [in %] for the signal ggF $HH$ process as a function of the coupling modifier $\kappa_{\lambda}$ in each analysis category for the $\tau_{\mathrm{had}}\tau_{\mathrm{had}}$ channel. The other coupling modifiers affecting ggF $HH$ production are fixed to their SM predictions.

The acceptance times efficiency [in %] for the signal ggF $HH$ process as a function of the coupling modifier $\kappa_{\lambda}$ in each analysis category for the $\tau_{\mathrm{lep}}\tau_{\mathrm{had}}$ SLT channel. The other coupling modifiers affecting ggF $HH$ production are fixed to their SM predictions.

The acceptance times efficiency [in %] for the signal ggF $HH$ process as a function of the coupling modifier $\kappa_{\lambda}$ in each analysis category for the $\tau_{\mathrm{lep}}\tau_{\mathrm{had}}$ LTT channel. The other coupling modifiers affecting ggF $HH$ production are fixed to their SM predictions.

The acceptance times efficiency [in %] for the signal VBF $HH$ process as a function of the coupling modifier $\kappa_{\lambda}$ in each analysis category for the $\tau_{\mathrm{had}}\tau_{\mathrm{had}}$ channel. The other coupling modifiers affecting VBF $HH$ production are fixed to their SM predictions.

The acceptance times efficiency [in %] for the signal VBF $HH$ process as a function of the coupling modifier $\kappa_{\lambda}$ in each analysis category for the $\tau_{\mathrm{lep}}\tau_{\mathrm{had}}$ SLT channel. The other coupling modifiers affecting VBF $HH$ production are fixed to their SM predictions.

The acceptance times efficiency [in %] for the signal VBF $HH$ process as a function of the coupling modifier $\kappa_{\lambda}$ in each analysis category for the $\tau_{\mathrm{lep}}\tau_{\mathrm{had}}$ LTT channel. The other coupling modifiers affecting VBF $HH$ production are fixed to their SM predictions.

The acceptance times efficiency [in %] for the signal VBF $HH$ process as a function of the coupling modifier $\kappa_{2V}$ in each analysis category for the $\tau_{\mathrm{had}}\tau_{\mathrm{had}}$ channel. The other coupling modifiers affecting VBF $HH$ production are fixed to their SM predictions.

The acceptance times efficiency [in %] for the signal VBF $HH$ process as a function of the coupling modifier $\kappa_{2V}$ in each analysis category for the $\tau_{\mathrm{lep}}\tau_{\mathrm{had}}$ SLT channel. The other coupling modifiers affecting VBF $HH$ production are fixed to their SM predictions.

The acceptance times efficiency [in %] for the signal VBF $HH$ process as a function of the coupling modifier $\kappa_{2V}$ in each analysis category for the $\tau_{\mathrm{lep}}\tau_{\mathrm{had}}$ LTT channel. The other coupling modifiers affecting VBF $HH$ production are fixed to their SM predictions.

Expected (dashed line) value of $-2\ln\Lambda$ as a function of $\kappa_{\lambda}$ for the fit to the $\tau_{\mathrm{had}}\tau_{\mathrm{had}}$ channel, when all other coupling modifiers are fixed to their SM predictions.

Expected (dashed line) value of $-2\ln\Lambda$ as a function of $\kappa_{\lambda}$ for the fit to the $\tau_{\mathrm{lep}}\tau_{\mathrm{had}}$ channel, when all other coupling modifiers are fixed to their SM predictions.

Observed (solid line) value of $-2\ln\Lambda$ as a function of $\kappa_{\lambda}$ for the fit to the $\tau_{\mathrm{had}}\tau_{\mathrm{had}}$ channel, when all other coupling modifiers are fixed to their SM predictions.

Observed (solid line) value of $-2\ln\Lambda$ as a function of $\kappa_{\lambda}$ for the fit to the $\tau_{\mathrm{lep}}\tau_{\mathrm{had}}$ channel, when all other coupling modifiers are fixed to their SM predictions.

Expected (dashed line) value of $-2\ln\Lambda$ as a function of $\kappa_{2V}$ for the fit to the $\tau_{\mathrm{had}}\tau_{\mathrm{had}}$ channel, when all other coupling modifiers are fixed to their SM predictions.

Expected (dashed line) value of $-2\ln\Lambda$ as a function of $\kappa_{2V}$ for the fit to the $\tau_{\mathrm{lep}}\tau_{\mathrm{had}}$ channel, when all other coupling modifiers are fixed to their SM predictions.

Observed (solid line) value of $-2\ln\Lambda$ as a function of $\kappa_{2V}$ for the fit to the $\tau_{\mathrm{had}}\tau_{\mathrm{had}}$ channel, when all other coupling modifiers are fixed to their SM predictions.

Observed (solid line) value of $-2\ln\Lambda$ as a function of $\kappa_{2V}$ for the fit to the $\tau_{\mathrm{lep}}\tau_{\mathrm{had}}$ channel, when all other coupling modifiers are fixed to their SM predictions.

Summary of the expected (blue) and observed (orange) one-dimensional constraints at 68% (solid error bars) and 95% (dashed error bars) CL on the coefficients of selected HEFT (top) and SMEFT (bottom) operators describing Higgs boson interactions affecting $HH$ production through gluon-gluon fusion. The results are obtained from one dimensional scans of the profile log-likelihood assuming that all other Wilson coefficients are fixed to their SM values. The contribution from VBF $HH$ production is neglected.

The performance of jet energy resolution (JER) for jets with |y| < 2.1 evaluated as a function of pT_truth in different centrality bins. Simulated hard scatter events were overlaid onto events from a dedicated sample of minimum-bias Xe+Xe data. The fit parameters are listed in a sperate table (Extras 1)

The relative magnitude of systematic uncertainties for per-pair normalized xJ distributions in 0-10% Xe+Xe centrality

The relative magnitude of systematic uncertainties for absolutely normalized xJ distributions in 0-10% Xe+Xe centrality

The relative magnitude of systematic uncertainties for rho distributions for leading jets in 0-10% Xe+Xe centrality

Per-pair normalized xJ distribution evaluated in four centrality intervals and given pT1 interval.

Per-pair normalized xJ distribution evaluated in four centrality intervals and given pT1 interval.

Per-pair normalized xJ distribution evaluated in four centrality intervals and given pT1 interval.

Absolutely normalized xJ distribution evaluated in four centrality intervals and given pT1 interval.

Absolutely normalized xJ distribution evaluated in four centrality intervals and given pT1 interval.

Absolutely normalized xJ distribution evaluated in four centrality intervals and given pT1 interval.

Per-pair normalized xJ distribution in Xe+Xe collisions.

Per-pair normalized xJ distribution in Pb+Pb collisions.

Per-pair normalized xJ distribution in Xe+Xe collisions.

Per-pair normalized xJ distribution in Pb+Pb collisions.

Per-pair normalized xJ distribution in Xe+Xe collisions.

Per-pair normalized xJ distribution in Pb+Pb collisions.

Absolutely normalized xJ distribution in Xe+Xe collisions.

Absolutely normalized xJ distribution in Xe+Xe collisions.

Absolutely normalized xJ distribution in Xe+Xe collisions.

Absolutely normalized xJ distribution in Xe+Xe collisions.

Absolutely normalized xJ distribution in Xe+Xe collisions.

Absolutely normalized xJ distribution in Xe+Xe collisions.

The ratios of Xe+Xe and Pb+Pb pair nuclear-modification factors, rho, evaluated as a function of leading jet pT in the same centrality intervals.

The ratios of Xe+Xe and Pb+Pb pair nuclear-modification factors, rho, evaluated as a function of leading jet pT in the same centrality intervals.

The ratios of Xe+Xe and Pb+Pb pair nuclear-modification factors, rho, evaluated as a function of leading jet pT in the same centrality intervals.

The ratios of Xe+Xe and Pb+Pb pair nuclear-modification factors, rho, evaluated as a function of subleading jet pT in the same centrality intervals.

The ratios of Xe+Xe and Pb+Pb pair nuclear-modification factors, rho, evaluated as a function of subleading jet pT in the same centrality intervals.

The ratios of Xe+Xe and Pb+Pb pair nuclear-modification factors, rho, evaluated as a function of subleading jet pT in the same centrality intervals.

The ratios of Xe+Xe and Pb+Pb pair nuclear-modification factors, rho, evaluated as a function of leading jet pT in the same SUM ETFCal intervals (selecting equivalent event activity)

The ratios of Xe+Xe and Pb+Pb pair nuclear-modification factors, rho, evaluated as a function of leading jet pT in the same SUM ETFCal intervals (selecting equivalent event activity)

The ratios of Xe+Xe and Pb+Pb pair nuclear-modification factors, rho, evaluated as a function of leading jet pT in the same SUM ETFCal intervals (selecting equivalent event activity)

The ratios of Xe+Xe and Pb+Pb pair nuclear-modification factors, rho, evaluated as a function of subleading jet pT in the same SUM ETFCal intervals (selecting equivalent event activity)

The ratios of Xe+Xe and Pb+Pb pair nuclear-modification factors, rho, evaluated as a function of subleading jet pT in the same SUM ETFCal intervals (selecting equivalent event activity)

The ratios of Xe+Xe and Pb+Pb pair nuclear-modification factors, rho, evaluated as a function of subleading jet pT in the same SUM ETFCal intervals (selecting equivalent event activity)

Parameter a,b, and c from JER fits in Figure 1b.

Per-pair normalized xJ distribution in Xe+Xe collisions for selected Pb+Pb centrality and pT1 bin.

Per-pair normalized xJ distribution in Xe+Xe collisions for selected Pb+Pb centrality and pT1 bin.

Per-pair normalized xJ distribution in Xe+Xe collisions for selected Pb+Pb centrality and pT1 bin.

Per-pair normalized xJ distribution in Xe+Xe collisions for selected Pb+Pb centrality and pT1 bin.

Per-pair normalized xJ distribution in Xe+Xe collisions for selected Pb+Pb centrality and pT1 bin.

Per-pair normalized xJ distribution in Xe+Xe collisions for selected Pb+Pb centrality and pT1 bin.

Per-pair normalized xJ distribution in Xe+Xe collisions for selected Pb+Pb centrality and pT1 bin.

Per-pair normalized xJ distribution in Xe+Xe collisions for selected Pb+Pb centrality and pT1 bin.

Per-pair normalized xJ distribution in Xe+Xe collisions for selected Pb+Pb centrality and pT1 bin.

The unfolded differential charge asymmetry as a function of the transverse momentum of the top pair system. The measured values are given with statistical and systematic uncertainties. The SM theory predictions calculated at NNLO in QCD and NLO in EW theory are listed. The scale uncertainty is obtained by varying renormalisation and factorisation scales independently by a factor of 2 or 0.5 around $\mu_0$ to calculate the maximum and minimum value of the asymmetry, respectively. The nominal value $\mu_0$ is chosen as $H_T/4$. The variations in which one scale is multiplied by 2 while the other scale is divided by 2 are excluded. Finally, the scale and MC integration uncertainties are added in quadrature.

The unfolded differential charge asymmetry as a function of the longitudinal boost of the top pair system. The measured values are given with statistical and systematic uncertainties. The SM theory predictions calculated at NNLO in QCD and NLO in EW theory are listed. The scale uncertainty is obtained by varying renormalisation and factorisation scales independently by a factor of 2 or 0.5 around $\mu_0$ to calculate the maximum and minimum value of the asymmetry, respectively. The nominal value $\mu_0$ is chosen as $H_T/4$. The variations in which one scale is multiplied by 2 while the other scale is divided by 2 are excluded. Finally, the scale and MC integration uncertainties are added in quadrature.

The unfolded inclusive leptonic asymmetry. The unfolded $A_C^{\ell\bar{\ell}}$ is obtained in the reduced phase-space defined by the requirement $|\Delta |\eta_{\ell\bar{\ell}}||<2.5$. The measured values are given with statistical and systematic uncertainties. The SM theory predictions calculated at NLO in QCD and NLO in EW theory are listed. The theory uncertainty is obtained by varying both scales by a factor of 0.5 or 2.0 to calculate the minimum and maximum value of the asymmetry, respectively.

The unfolded differential leptonic asymmetry as a function of the invariant mass of the di-lepton pair. The unfolded $A_C^{\ell\bar{\ell}}$ is obtained in the reduced phase-space defined by the requirement $|\Delta |\eta_{\ell\bar{\ell}}||<2.5$. The measured values are given with statistical and systematic uncertainties. The SM theory predictions calculated at NLO in QCD and NLO in EW theory are listed. The theory uncertainty is obtained by varying both scales by a factor of 0.5 or 2.0 to calculate the minimum and maximum value of the asymmetry, respectively.

The unfolded differential leptonic asymmetry as a function of the transverse momentum of the di-lepton pair. The unfolded $A_C^{\ell\bar{\ell}}$ is obtained in the reduced phase-space defined by the requirement $|\Delta |\eta_{\ell\bar{\ell}}||<2.5$. The measured values are given with statistical and systematic uncertainties. The SM theory predictions calculated at NLO in QCD and NLO in EW theory are listed. The theory uncertainty is obtained by varying both scales by a factor of 0.5 or 2.0 to calculate the minimum and maximum value of the asymmetry, respectively.

The unfolded differential leptonic asymmetry as a function of the longitudinal boost of the di-lepton pair. The unfolded $A_C^{\ell\bar{\ell}}$ is obtained in the reduced phase-space defined by the requirement $|\Delta |\eta_{\ell\bar{\ell}}||<2.5$. The measured values are given with statistical and systematic uncertainties. The SM theory predictions calculated at NLO in QCD and NLO in EW theory are listed. The theory uncertainty is obtained by varying both scales by a factor of 0.5 or 2.0 to calculate the minimum and maximum value of the asymmetry, respectively.

Individual 68% and 95% CL bounds on the relevant Wilson coefficients of the SM Effective Field Theory in units of $\text{TeV}^{-2}$. The bounds are derived from the $A_C^{t\bar{t}}$ inclusive measurement. The experimental uncertainties are accounted for, in the form of the complete covariance matrix that keeps track of correlations between bins for the differential measurement. The theory uncertainty from the NNLO QCD + NLO EW calculation is included by explicitly varying the renormalization and factorization scales, or the parton density functions, in the calculation and registering the variations in the intervals.

Individual 68% and 95% CL bounds on the relevant Wilson coefficients of the SM Effective Field Theory in units of $\text{TeV}^{-2}$. The bounds are derived from the $A_C^{t\bar{t}}$ vs $m_{t\bar{t}}$ measurement. The experimental uncertainties are accounted for, in the form of the complete covariance matrix that keeps track of correlations between bins for the differential measurement. The theory uncertainty from the NNLO QCD + NLO EW calculation is included by explicitly varying the renormalization and factorization scales, or the parton density functions, in the calculation and registering the variations in the intervals.

Ranking of the systematic uncertainties with marginalisation for the $A_C^{t\bar{t}}$ inclusive measurement. The effect on unfolded $A_C$ for down and up variation of the systematic uncertainty is shown, respectively. The pulls and constraints of the ranked NPs are obtained from data.

Ranking of the systematic uncertainties with marginalisation for the $A_C^{t\bar{t}}$ vs $\beta_{z,t\bar{t}}$ measurement for $\beta_{z,t\bar{t}}$ $\in$ [0,0.3]. The effect on unfolded $A_C$ for down and up variation of the systematic uncertainty is shown, respectively. The pulls and constraints of the ranked NPs are obtained from data.

Ranking of the systematic uncertainties with marginalisation for the $A_C^{t\bar{t}}$ vs $\beta_{z,t\bar{t}}$ measurement for $\beta_{z,t\bar{t}}$ $\in$ [0.3,0.6]. The effect on unfolded $A_C$ for down and up variation of the systematic uncertainty is shown, respectively. The pulls and constraints of the ranked NPs are obtained from data.

Ranking of the systematic uncertainties with marginalisation for the $A_C^{t\bar{t}}$ vs $\beta_{z,t\bar{t}}$ measurement for $\beta_{z,t\bar{t}}$ $\in$ [0.6,0.8]. The effect on unfolded $A_C$ for down and up variation of the systematic uncertainty is shown, respectively. The pulls and constraints of the ranked NPs are obtained from data.

Ranking of the systematic uncertainties with marginalisation for the $A_C^{t\bar{t}}$ vs $\beta_{z,t\bar{t}}$ measurement for $\beta_{z,t\bar{t}}$ $\in$ [0.8,1]. The effect on unfolded $A_C$ for down and up variation of the systematic uncertainty is shown, respectively. The pulls and constraints of the ranked NPs are obtained from data.

Ranking of the systematic uncertainties with marginalisation for the $A_C^{t\bar{t}}$ vs $m_{t\bar{t}}$ measurement for $m_{t\bar{t}}$ < 500 GeV. The effect on unfolded $A_C$ for down and up variation of the systematic uncertainty is shown, respectively. The pulls and constraints of the ranked NPs are obtained from data.

Ranking of the systematic uncertainties with marginalisation for the $A_C^{t\bar{t}}$ vs $m_{t\bar{t}}$ measurement for $m_{t\bar{t}}$ $\in$ [500,750] GeV. The effect on unfolded $A_C$ for down and up variation of the systematic uncertainty is shown, respectively. The pulls and constraints of the ranked NPs are obtained from data.

Ranking of the systematic uncertainties with marginalisation for the $A_C^{t\bar{t}}$ vs $m_{t\bar{t}}$ measurement for $m_{t\bar{t}}$ $\in$ [750,1000] GeV. The effect on unfolded $A_C$ for down and up variation of the systematic uncertainty is shown, respectively. The pulls and constraints of the ranked NPs are obtained from data.

Ranking of the systematic uncertainties with marginalisation for the $A_C^{t\bar{t}}$ vs $m_{t\bar{t}}$ measurement for $m_{t\bar{t}}$ $\in$ [1000,1500] GeV. The effect on unfolded $A_C$ for down and up variation of the systematic uncertainty is shown, respectively. The pulls and constraints of the ranked NPs are obtained from data.

Ranking of the systematic uncertainties with marginalisation for the $A_C^{t\bar{t}}$ vs $m_{t\bar{t}}$ measurement for $m_{t\bar{t}}$ > 1500 GeV. The effect on unfolded $A_C$ for down and up variation of the systematic uncertainty is shown, respectively. The pulls and constraints of the ranked NPs are obtained from data.

Ranking of the systematic uncertainties with marginalisation for the $A_C^{t\bar{t}}$ vs $p_{T,t\bar{t}}$ measurement for $p_{T,t\bar{t}}$ < 30 GeV. The effect on unfolded $A_C$ for down and up variation of the systematic uncertainty is shown, respectively. The pulls and constraints of the ranked NPs are obtained from data.

Ranking of the systematic uncertainties with marginalisation for the $A_C^{t\bar{t}}$ vs $p_{T,t\bar{t}}$ measurement for $p_{T,t\bar{t}}$ $\in$ [30,120] GeV. The effect on unfolded $A_C$ for down and up variation of the systematic uncertainty is shown, respectively. The pulls and constraints of the ranked NPs are obtained from data.

Ranking of the systematic uncertainties with marginalisation for the $A_C^{t\bar{t}}$ vs $p_{T,t\bar{t}}$ measurement for $p_{T,t\bar{t}}$ > 120 GeV. The effect on unfolded $A_C$ for down and up variation of the systematic uncertainty is shown, respectively. The pulls and constraints of the ranked NPs are obtained from data.

Ranking of the systematic uncertainties with marginalisation for the $A_C^{\ell\ell}$ inclusive measurement. The effect on unfolded $A_C$ for down and up variation of the systematic uncertainty is shown, respectively. The pulls and constraints of the ranked NPs are obtained from data.

Ranking of the systematic uncertainties with marginalisation for the $A_C^{\ell\ell}$ vs $\beta_{z,\ell\bar{\ell}}$ measurement for $\beta_{z,\ell\bar{\ell}}$ $\in$[0,0.3]. The effect on unfolded $A_C$ for down and up variation of the systematic uncertainty is shown, respectively. The pulls and constraints of the ranked NPs are obtained from data.

Ranking of the systematic uncertainties with marginalisation for the $A_C^{\ell\ell}$ vs $\beta_{z,\ell\bar{\ell}}$ measurement for $\beta_{z,\ell\bar{\ell}}$ $\in$[0.3,0.6]. The effect on unfolded $A_C$ for down and up variation of the systematic uncertainty is shown, respectively. The pulls and constraints of the ranked NPs are obtained from data.

Ranking of the systematic uncertainties with marginalisation for the $A_C^{\ell\ell}$ vs $\beta_{z,\ell\bar{\ell}}$ measurement for $\beta_{z,\ell\bar{\ell}}$ $\in$[0.6,0.8]. The effect on unfolded $A_C$ for down and up variation of the systematic uncertainty is shown, respectively. The pulls and constraints of the ranked NPs are obtained from data.

Ranking of the systematic uncertainties with marginalisation for the $A_C^{\ell\ell}$ vs $\beta_{z,\ell\bar{\ell}}$ measurement for $\beta_{z,\ell\bar{\ell}}$ $\in$[0.8,1]. The effect on unfolded $A_C$ for down and up variation of the systematic uncertainty is shown, respectively. The pulls and constraints of the ranked NPs are obtained from data.

Ranking of the systematic uncertainties with marginalisation for the $A_C^{\ell\ell}$ vs $m_{\ell\bar{\ell}}$ measurement for $m_{\ell\bar{\ell}}$ < 200 GeV. The effect on unfolded $A_C$ for down and up variation of the systematic uncertainty is shown, respectively. The pulls and constraints of the ranked NPs are obtained from data.

Ranking of the systematic uncertainties with marginalisation for the $A_C^{\ell\ell}$ vs $m_{\ell\bar{\ell}}$ measurement for $m_{\ell\bar{\ell}}$ $\in$ [200,300] GeV. The effect on unfolded $A_C$ for down and up variation of the systematic uncertainty is shown, respectively. The pulls and constraints of the ranked NPs are obtained from data.

Ranking of the systematic uncertainties with marginalisation for the $A_C^{\ell\ell}$ vs $m_{\ell\bar{\ell}}$ measurement for $m_{\ell\bar{\ell}}$ $\in$ [300,400] GeV. The effect on unfolded $A_C$ for down and up variation of the systematic uncertainty is shown, respectively. The pulls and constraints of the ranked NPs are obtained from data.

Ranking of the systematic uncertainties with marginalisation for the $A_C^{\ell\ell}$ vs $m_{\ell\bar{\ell}}$ measurement for $m_{\ell\bar{\ell}}$ > 400 GeV. The effect on unfolded $A_C$ for down and up variation of the systematic uncertainty is shown, respectively. The pulls and constraints of the ranked NPs are obtained from data.

Ranking of the systematic uncertainties with marginalisation for the $A_C^{\ell\ell}$ vs $p_{T,\ell\bar{\ell}}$ measurement for $p_{T,\ell\bar{\ell}}$ < 20 GeV. The effect on unfolded $A_C$ for down and up variation of the systematic uncertainty is shown, respectively. The pulls and constraints of the ranked NPs are obtained from data.

Ranking of the systematic uncertainties with marginalisation for the $A_C^{\ell\ell}$ vs $p_{T,\ell\bar{\ell}}$ measurement for $p_{T,\ell\bar{\ell}}$ $\in$ [20, 70] GeV. The effect on unfolded $A_C$ for down and up variation of the systematic uncertainty is shown, respectively. The pulls and constraints of the ranked NPs are obtained from data.

Ranking of the systematic uncertainties with marginalisation for the $A_C^{\ell\ell}$ vs $p_{T,\ell\bar{\ell}}$ measurement for $p_{T,\ell\bar{\ell}}$ > 70 GeV. The effect on unfolded $A_C$ for down and up variation of the systematic uncertainty is shown, respectively. The pulls and constraints of the ranked NPs are obtained from data.

Post-marginalisation correlation coefficients $\rho_{ij}$ of nuisance parameters for the $A_C^{t\bar{t}}$ inclusive measurement. Only $|\rho_{ij}| > 0.05$ values are included.

Post-marginalisation correlation coefficients $\rho_{ij}$ of nuisance parameters for the $A_C^{t\bar{t}}$ vs $m_{t\bar{t}}$ measurement. Only $|\rho_{ij}| > 0.05$ values are included.

Post-marginalisation correlation coefficients $\rho_{ij}$ of nuisance parameters for the $A_C^{t\bar{t}}$ vs $p_{T,t\bar{t}}$ measurement. Only $|\rho_{ij}| > 0.05$ values are included.

Post-marginalisation correlation coefficients $\rho_{ij}$ of nuisance parameters for the $A_C^{t\bar{t}}$ vs $\beta_{z,t\bar{t}}$ measurement. Only $|\rho_{ij}| > 0.05$ values are included.

Post-marginalisation correlation coefficients $\rho_{ij}$ of nuisance parameters for the $A_C^{\ell\ell}$ inclusive measurement. Only $|\rho_{ij}| > 0.05$ values are included.

Post-marginalisation correlation coefficients $\rho_{ij}$ of nuisance parameters for the $A_C^{\ell\ell}$ vs $m_{\ell\bar{\ell}}$ measurement. Only $|\rho_{ij}| > 0.05$ values are included.

Post-marginalisation correlation coefficients $\rho_{ij}$ of nuisance parameters for the $A_C^{\ell\ell}$ vs $p_{T,\ell\bar{\ell}}$ measurement. Only $|\rho_{ij}| > 0.05$ values are included.

Post-marginalisation correlation coefficients $\rho_{ij}$ of nuisance parameters for the $A_C^{\ell\ell}$ vs $\beta_{z,\ell\bar{\ell}}$ measurement. Only $|\rho_{ij}| > 0.05$ values are included.

Covariance matrix for the $A_C^{t\bar{t}}$ vs $m_{t\bar{t}}$ measurement. The total (stat. + syst.) uncertainties are considered.

Covariance matrix for the $A_C^{t\bar{t}}$ vs $p_{T,t\bar{t}}$ measurement. The total (stat. + syst.) uncertainties are considered.

Covariance matrix for the $A_C^{t\bar{t}}$ vs $\beta_{z,t\bar{t}}$ measurement. The total (stat. + syst.) uncertainties are considered.

Covariance matrix for the $A_C^{\ell\ell}$ vs $m_{\ell\bar{\ell}}$ measurement. The total (stat. + syst.) uncertainties are considered.

Covariance matrix for the $A_C^{\ell\ell}$ vs $p_{T,\ell\bar{\ell}}$ measurement. The total (stat. + syst.) uncertainties are considered.

Covariance matrix for the $A_C^{\ell\ell}$ vs $\beta_{z,\ell\bar{\ell}}$ measurement. The total (stat. + syst.) uncertainties are considered.

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

The best-fit predicted and observed distributions of the combined NN output in the high-$p_\text{T}$-SR for the $\mu\tau_e$ channel. The first and last bin include underflow and overflow events, respectively.

The best-fit predicted and observed distributions of the collinear mass $m_\text{coll}$ in the high-$p_\text{T}$-SR for the $e\tau_\mu$ channel. The first and last bin include underflow and overflow events, respectively.

The best-fit predicted and observed distributions of the collinear mass $m_\text{coll}$ in the high-$p_\text{T}$-SR for the $\mu\tau_e$ channel. The first and last bin include underflow and overflow events, respectively.

Observed and expected upper limits on $\mathcal{B}(Z\rightarrow\ell\tau)$ with leptonically decaying $\tau$ at 95% confidence level.

Observed and expected upper limits on $\mathcal{B}(Z\rightarrow\ell\tau)$ at 95% confidence level, combination of hadronically and leptonically decaying $\tau$ channels.

The best-fit predicted and observed distributions of the collinear mass $m_\text{coll}$ in the CRZ$\tau\tau$ for the $e\tau_\mu$ channel. The first and last bin include underflow and overflow events, respectively.

The best-fit predicted and observed distributions of the collinear mass $m_\text{coll}$ in the CRZ$\tau\tau$ for the $\mu\tau_e$ channel. The first and last bin include underflow and overflow events, respectively.

The best-fit predicted and observed distributions of the combined NN output in the CRZ$\tau\tau$ for the $e\tau_\mu$ channel. The first and last bin include underflow and overflow events, respectively.

The best-fit predicted and observed distributions of the combined NN output in the CRZ$\tau\tau$ for the $\mu\tau_e$ channel. The first and last bin include underflow and overflow events, respectively.

The best-fit predicted and observed distributions of the leading lepton transverse momentum $p_\text{T}(e)$ in the high-$p_\text{T}$-SR for the $e\tau_\mu$ channel. The first and last bin include underflow and overflow events, respectively.

The best-fit predicted and observed distributions of the leading lepton transverse momentum $p_\text{T}(\mu)$ in the high-$p_\text{T}$-SR for the $\mu\tau_e$ channel. The first and last bin include underflow and overflow events, respectively.

The best-fit predicted and observed distributions of the missing transverse momentum $E^\text{miss}_\text{T}$ in the low-$p_\text{T}$-SR for the $e\tau_\mu$ channel. The first and last bin include underflow and overflow events, respectively.

The best-fit predicted and observed distributions of the missing transverse momentum $E^\text{miss}_\text{T}$ in the low-$p_\text{T}$-SR for the $\mu\tau_e$ channel. The first and last bin include underflow and overflow events, respectively.

The +1 $\sigma$ band of the observed exclusion contor for the $Z^{\prime}_{B}$ model in the $m_{\chi}$-$m_{Z^{\prime}_{B}}$ plane.

The -1 $\sigma$ band of the observed exclusion contor for the $Z^{\prime}_{B}$ model in the $m_{\chi}$-$m_{Z^{\prime}_{B}}$ plane.

A comparison of the inferred limits to the constraints from direct detection experiments on the spin-independent DM--nucleon cross section in the context of the $Z'_B$ simplified model with vector couplings. Limits are shown at 90% CL.

The observed exclusion contor for the $Z^{\prime}$-2HDM model in the $m_{A}$-$m_{Z^{\prime}}$ plane.

The expected exclusion contor for the $Z^{\prime}$-2HDM model in the $m_{A}$-$m_{Z^{\prime}}$ plane.

The +1 $\sigma$ band of the observed exclusion contor for the $m_{A}$-$Z^{\prime}$-2HDM model in the $m_{Z^{\prime}}$ plane.

The -1 $\sigma$ band of the observed exclusion contor for the $m_{A}$-$Z^{\prime}$-2HDM model in the $m_{Z^{\prime}}$ plane.

The observed exclusion contor for the 2HDM-a model in the $m_{A}$-$m_{a}$ plane.

The expected exclusion contor for the 2HDM-a model in the $m_{A}$-$m_{a}$ plane.

The +1 $\sigma$ band of the observed exclusion contor for the 2HDM-a model in the $m_{A}$-$m_{a}$ plane.

The -1 $\sigma$ band of the observed exclusion contor for the 2HDM-a model in the $m_{A}$-$m_{a}$ plane.

The observed exclusion contor for the 2HDM-a model in the $tan\beta$-$m_{a}$ plane.

The expected exclusion contor for the 2HDM-a model in the $tan\beta$-$m_{a}$ plane.

The +1 $\sigma$ band of the observed exclusion contor for the 2HDM-a model in the $tan\beta$-$m_{a}$ plane.

The -1 $\sigma$ band of the observed exclusion contor for the 2HDM-a model in the $tan\beta$-$m_{a}$ plane.

The exclusion limits at 95% CL for the 2HDM+a model as a function of $\sin \theta$ for $m_{A,H^{\pm},H}$= 600GeV, $m_a$ = 200GeV, $\tan \beta$ = 1.0$.

The exclusion limits at 95% CL for the 2HDM+a model as a function of $\sin \theta$ for $m_{A,H^{\pm},H}$= 1000GeV, $m_a$ = 350GeV, $\tan \beta$ = 1.0$.

Breakdown of the dominant systematic uncertainties.

Event yields in the range of 120 $<m_{\gamma\gamma}<$ 130 GeV for data, signal models, the SM Higgs boson background and non-resonant background in each analysis category, for an integrated luminosity of $139$fb$^{-1}$.

Detailed background contributions from the SM Higgs boson and continuum background for each cut

Detailed contributions from the signals for each cut.

Acceptance times efficiency for several signals in each category.

The definitions of the four tetralepton signal regions. The regions are defined to target different $b$-jet multiplicities and flavour combinations of the non-Z leptons.

The absolute particle-level differential cross-section measured in the fiducial phase-space as a function of the $p_{\textrm{T}}$ of the $Z$ boson in the 3$\ell$+4$\ell$ channel. The uncertainty is decomposed into four components which are the signal modelling uncertainty, the background modelling uncertainty, the experimental uncertainty, and the data statistical uncertainty.

The absolute parton-level differential cross-section measured in the fiducial phase-space as a function of the $p_{\textrm{T}}$ of the $Z$ boson in the 3$\ell$+4$\ell$ channel. The uncertainty is decomposed into four components which are the signal modelling uncertainty, the background modelling uncertainty, the experimental uncertainty, and the data statistical uncertainty.

The normalised particle-level differential cross-section measured in the fiducial phase-space as a function of the $p_{\textrm{T}}$ of the $Z$ boson in the 3$\ell$+4$\ell$ channel. The uncertainty is decomposed into four components which are the signal modelling uncertainty, the background modelling uncertainty, the experimental uncertainty, and the data statistical uncertainty.

The normalised parton-level differential cross-section measured in the fiducial phase-space as a function of the $p_{\textrm{T}}$ of the $Z$ boson in the 3$\ell$+4$\ell$ channel. The uncertainty is decomposed into four components which are the signal modelling uncertainty, the background modelling uncertainty, the experimental uncertainty, and the data statistical uncertainty.

The absolute parton-level differential cross-section measured in the fiducial phase-space as a function of the absolute value of rapidity of the $Z$ boson in the 3$\ell$+4$\ell$ channel. The uncertainty is decomposed into four components which are the signal modelling uncertainty, the background modelling uncertainty, the experimental uncertainty, and the data statistical uncertainty.

The absolute particle-level differential cross-section measured in the fiducial phase-space as a function of the number of jets in the 3$\ell$ channel. The uncertainty is decomposed into four components which are the signal modelling uncertainty, the background modelling uncertainty, the experimental uncertainty, and the data statistical uncertainty.

The absolute particle-level differential cross-section measured in the fiducial phase-space as a function of the number of jets in the 4$\ell$ channel. The uncertainty is decomposed into four components which are the signal modelling uncertainty, the background modelling uncertainty, the experimental uncertainty, and the data statistical uncertainty.

The absolute parton-level differential cross-section measured in the fiducial phase-space as a function of the $p_{\textrm{T}}^{l \textrm{non-}Z}$ in the 3$\ell$ channel. The uncertainty is decomposed into four components which are the signal modelling uncertainty, the background modelling uncertainty, the experimental uncertainty, and the data statistical uncertainty.

The absolute parton-level differential cross-section measured in the fiducial phase-space as a function of the $|\Delta \phi (Z, t_{\textrm{lep}})|/\pi$ in the 3$\ell$ channel. The uncertainty is decomposed into four components which are the signal modelling uncertainty, the background modelling uncertainty, the experimental uncertainty, and the data statistical uncertainty.

The absolute parton-level differential cross-section measured in the fiducial phase-space as a function of the $|\Delta y (Z, t_{\textrm{lep}})|$ in the 3$\ell$ channel. The uncertainty is decomposed into four components which are the signal modelling uncertainty, the background modelling uncertainty, the experimental uncertainty, and the data statistical uncertainty.

The absolute parton-level differential cross-section measured in the fiducial phase-space as a function of the $|\Delta \phi (l_{t}^{+}, l_{\bar{t}}^{-})|/\pi$ in the 4$\ell$ channel. The uncertainty is decomposed into four components which are the signal modelling uncertainty, the background modelling uncertainty, the experimental uncertainty, and the data statistical uncertainty.

The absolute parton-level differential cross-section measured in the fiducial phase-space as a function of the $|\Delta \phi (t\bar{t}, Z)|/\pi$ in the 4$\ell$ channel. The uncertainty is decomposed into four components which are the signal modelling uncertainty, the background modelling uncertainty, the experimental uncertainty, and the data statistical uncertainty.

The absolute parton-level differential cross-section measured in the fiducial phase-space as a function of the $p_{\textrm{T}}^{t\bar{t}}$ in the 4$\ell$ channel. The uncertainty is decomposed into four components which are the signal modelling uncertainty, the background modelling uncertainty, the experimental uncertainty, and the data statistical uncertainty.

The normalised parton-level differential cross-section measured in the fiducial phase-space as a function of the $p_{\textrm{T}}^{l \textrm{non-}Z}$ in the 3$\ell$ channel. The uncertainty is decomposed into four components which are the signal modelling uncertainty, the background modelling uncertainty, the experimental uncertainty, and the data statistical uncertainty.

The normalised parton-level differential cross-section measured in the fiducial phase-space as a function of the $|\Delta \phi (Z, t_{\textrm{lep}})|/\pi$ in the 3$\ell$ channel. The uncertainty is decomposed into four components which are the signal modelling uncertainty, the background modelling uncertainty, the experimental uncertainty, and the data statistical uncertainty.

The normalised parton-level differential cross-section measured in the fiducial phase-space as a function of the $|\Delta y (Z, t_{\textrm{lep}})|$ in the 3$\ell$ channel. The uncertainty is decomposed into four components which are the signal modelling uncertainty, the background modelling uncertainty, the experimental uncertainty, and the data statistical uncertainty.

The normalised parton-level differential cross-section measured in the fiducial phase-space as a function of the $|\Delta \phi (l_{t}^{+}, l_{\bar{t}}^{-})|/\pi$ in the 4$\ell$ channel. The uncertainty is decomposed into four components which are the signal modelling uncertainty, the background modelling uncertainty, the experimental uncertainty, and the data statistical uncertainty.

The normalised parton-level differential cross-section measured in the fiducial phase-space as a function of the $|\Delta \phi (t\bar{t}, Z)|/\pi$ in the 4$\ell$ channel. The uncertainty is decomposed into four components which are the signal modelling uncertainty, the background modelling uncertainty, the experimental uncertainty, and the data statistical uncertainty.

The normalised parton-level differential cross-section measured in the fiducial phase-space as a function of the $p_{\textrm{T}}^{t\bar{t}}$ in the 4$\ell$ channel. The uncertainty is decomposed into four components which are the signal modelling uncertainty, the background modelling uncertainty, the experimental uncertainty, and the data statistical uncertainty.

The normalised parton-level differential cross-section measured in the fiducial phase-space as a function of the absolute value of rapidity of the $Z$ boson in the 3$\ell$+4$\ell$ channel. The uncertainty is decomposed into four components which are the signal modelling uncertainty, the background modelling uncertainty, the experimental uncertainty, and the data statistical uncertainty.

The absolute particle-level differential cross-section measured in the fiducial phase-space as a function of the $p_{\textrm{T}}^{l \textrm{non-}Z}$ in the 3$\ell$ channel. The uncertainty is decomposed into four components which are the signal modelling uncertainty, the background modelling uncertainty, the experimental uncertainty, and the data statistical uncertainty.

The absolute particle-level differential cross-section measured in the fiducial phase-space as a function of the $|\Delta \phi (Z, t_{\textrm{lep}})|/\pi$ in the 3$\ell$ channel. The uncertainty is decomposed into four components which are the signal modelling uncertainty, the background modelling uncertainty, the experimental uncertainty, and the data statistical uncertainty.

The absolute particle-level differential cross-section measured in the fiducial phase-space as a function of the $|\Delta y (Z, t_{\textrm{lep}})|$ in the 3$\ell$ channel. The uncertainty is decomposed into four components which are the signal modelling uncertainty, the background modelling uncertainty, the experimental uncertainty, and the data statistical uncertainty.

The absolute particle-level differential cross-section measured in the fiducial phase-space as a function of the $|\Delta \phi (l_{t}^{+}, l_{\bar{t}}^{-})|/\pi$ in the 4$\ell$ channel. The uncertainty is decomposed into four components which are the signal modelling uncertainty, the background modelling uncertainty, the experimental uncertainty, and the data statistical uncertainty.

The absolute particle-level differential cross-section measured in the fiducial phase-space as a function of the $|\Delta \phi (t\bar{t}, Z)|/\pi$ in the 4$\ell$ channel. The uncertainty is decomposed into four components which are the signal modelling uncertainty, the background modelling uncertainty, the experimental uncertainty, and the data statistical uncertainty.

The absolute particle-level differential cross-section measured in the fiducial phase-space as a function of the $p_{\textrm{T}}^{t\bar{t}}$ in the 4$\ell$ channel. The uncertainty is decomposed into four components which are the signal modelling uncertainty, the background modelling uncertainty, the experimental uncertainty, and the data statistical uncertainty.

The absolute particle-level differential cross-section measured in the fiducial phase-space as a function of the absolute value of rapidity of the $Z$ boson in the 3$\ell$+4$\ell$ channel. The uncertainty is decomposed into four components which are the signal modelling uncertainty, the background modelling uncertainty, the experimental uncertainty, and the data statistical uncertainty.

The normalised particle-level differential cross-section measured in the fiducial phase-space as a function of the $p_{\textrm{T}}^{l \textrm{non-}Z}$ in the 3$\ell$ channel. The uncertainty is decomposed into four components which are the signal modelling uncertainty, the background modelling uncertainty, the experimental uncertainty, and the data statistical uncertainty.

The normalised particle-level differential cross-section measured in the fiducial phase-space as a function of the $|\Delta \phi (Z, t_{\textrm{lep}})|/\pi$ in the 3$\ell$ channel. The uncertainty is decomposed into four components which are the signal modelling uncertainty, the background modelling uncertainty, the experimental uncertainty, and the data statistical uncertainty.

The normalised particle-level differential cross-section measured in the fiducial phase-space as a function of the $|\Delta y (Z, t_{\textrm{lep}})|$ in the 3$\ell$ channel. The uncertainty is decomposed into four components which are the signal modelling uncertainty, the background modelling uncertainty, the experimental uncertainty, and the data statistical uncertainty.

The normalised particle-level differential cross-section measured in the fiducial phase-space as a function of the $|\Delta \phi (l_{t}^{+}, l_{\bar{t}}^{-})|/\pi$ in the 4$\ell$ channel. The uncertainty is decomposed into four components which are the signal modelling uncertainty, the background modelling uncertainty, the experimental uncertainty, and the data statistical uncertainty.

The normalised particle-level differential cross-section measured in the fiducial phase-space as a function of the $|\Delta \phi (t\bar{t}, Z)|/\pi$ in the 4$\ell$ channel. The uncertainty is decomposed into four components which are the signal modelling uncertainty, the background modelling uncertainty, the experimental uncertainty, and the data statistical uncertainty.

The normalised particle-level differential cross-section measured in the fiducial phase-space as a function of the $p_{\textrm{T}}^{t\bar{t}}$ in the 4$\ell$ channel. The uncertainty is decomposed into four components which are the signal modelling uncertainty, the background modelling uncertainty, the experimental uncertainty, and the data statistical uncertainty.

The normalised particle-level differential cross-section measured in the fiducial phase-space as a function of the absolute value of rapidity of the $Z$ boson in the 3$\ell$+4$\ell$ channel. The uncertainty is decomposed into four components which are the signal modelling uncertainty, the background modelling uncertainty, the experimental uncertainty, and the data statistical uncertainty.

The normalised particle-level differential cross-section measured in the fiducial phase-space as a function of the number of jets in the 3$\ell$ channel. The uncertainty is decomposed into four components which are the signal modelling uncertainty, the background modelling uncertainty, the experimental uncertainty, and the data statistical uncertainty.

The normalised particle-level differential cross-section measured in the fiducial phase-space as a function of the number of jets in the 4$\ell$ channel. The uncertainty is decomposed into four components which are the signal modelling uncertainty, the background modelling uncertainty, the experimental uncertainty, and the data statistical uncertainty.

The total correlation matrix of the absolute parton-level differential cross-section measured in the fiducial phase-space as a function of the $p_{\textrm{T}}$ of the $Z$ boson in the 3$\ell$+4$\ell$ channel.

The total correlation matrix of the absolute parton-level differential cross-section measured in the fiducial phase-space as a function of the absolute value of rapidity of the $Z$ boson in the 3$\ell$+4$\ell$ channel.

The total correlation matrix of the absolute parton-level differential cross-section measured in the fiducial phase-space as a function of the $p_{\textrm{T}}^{l \textrm{non-}Z}$ in the 3$\ell$ channel.

The total correlation matrix of the absolute parton-level differential cross-section measured in the fiducial phase-space as a function of the $|\Delta \phi (Z, t_{\textrm{lep}})|/\pi$ in the 3$\ell$ channel.

The total correlation matrix of the absolute parton-level differential cross-section measured in the fiducial phase-space as a function of the $|\Delta y (Z, t_{\textrm{lep}})|/\pi$ in the 3$\ell$ channel.

The total correlation matrix of the absolute parton-level differential cross-section measured in the fiducial phase-space as a function of the $|\Delta \phi (l_{t}^{+}, l_{\bar{t}}^{-})|/\pi$ in the 4$\ell$ channel.

The total correlation matrix of the absolute parton-level differential cross-section measured in the fiducial phase-space as a function of the $|\Delta \phi (t\bar{t}, Z)|/\pi$ in the 4$\ell$ channel.

The total correlation matrix of the absolute parton-level differential cross-section measured in the fiducial phase-space as a function of the $p_{\textrm{T}}^{t\bar{t}}$ in the 4$\ell$ channel.

The total correlation matrix of the normalised parton-level differential cross-section measured in the fiducial phase-space as a function of the $p_{\textrm{T}}$ of the $Z$ boson in the 3$\ell$+4$\ell$ channel.

The total correlation matrix of the normalised parton-level differential cross-section measured in the fiducial phase-space as a function of the absolute value of rapidity of the $Z$ boson in the 3$\ell$+4$\ell$ channel.

The total correlation matrix of the normalised parton-level differential cross-section measured in the fiducial phase-space as a function of the $p_{\textrm{T}}^{l \textrm{non-}Z}$ in the 3$\ell$ channel.

The total correlation matrix of the normalised parton-level differential cross-section measured in the fiducial phase-space as a function of the $|\Delta \phi (Z, t_{\textrm{lep}})|/\pi$ in the 3$\ell$ channel.

The total correlation matrix of the normalised parton-level differential cross-section measured in the fiducial phase-space as a function of the $|\Delta y (Z, t_{\textrm{lep}})|/\pi$ in the 3$\ell$ channel.

The total correlation matrix of the normalised parton-level differential cross-section measured in the fiducial phase-space as a function of the $|\Delta \phi (l_{t}^{+}, l_{\bar{t}}^{-})|/\pi$ in the 4$\ell$ channel.

The total correlation matrix of the normalised parton-level differential cross-section measured in the fiducial phase-space as a function of the $|\Delta \phi (t\bar{t}, Z)|/\pi$ in the 4$\ell$ channel.

The total correlation matrix of the normalised parton-level differential cross-section measured in the fiducial phase-space as a function of the $p_{\textrm{T}}^{t\bar{t}}$ in the 4$\ell$ channel.

The total correlation matrix of the absolute particle-level differential cross-section measured in the fiducial phase-space as a function of the $p_{\textrm{T}}$ of the $Z$ boson in the 3$\ell$+4$\ell$ channel.

The total correlation matrix of the absolute particle-level differential cross-section measured in the fiducial phase-space as a function of the absolute value of rapidity of the $Z$ boson in the 3$\ell$+4$\ell$ channel.

The total correlation matrix of the absolute particle-level differential cross-section measured in the fiducial phase-space as a function of the $p_{\textrm{T}}^{l \textrm{non-}Z}$ in the 3$\ell$ channel.

The total correlation matrix of the absolute particle-level differential cross-section measured in the fiducial phase-space as a function of the $|\Delta \phi (Z, t_{\textrm{lep}})|/\pi$ in the 3$\ell$ channel.

The total correlation matrix of the absolute particle-level differential cross-section measured in the fiducial phase-space as a function of the $|\Delta y (Z, t_{\textrm{lep}})|/\pi$ in the 3$\ell$ channel.

The total correlation matrix of the absolute particle-level differential cross-section measured in the fiducial phase-space as a function of the number of jets in the 3$\ell$ channel.

The total correlation matrix of the absolute particle-level differential cross-section measured in the fiducial phase-space as a function of the $|\Delta \phi (l_{t}^{+}, l_{\bar{t}}^{-})|/\pi$ in the 4$\ell$ channel.

The total correlation matrix of the absolute particle-level differential cross-section measured in the fiducial phase-space as a function of the $|\Delta \phi (t\bar{t}, Z)|/\pi$ in the 4$\ell$ channel.

The total correlation matrix of the absolute particle-level differential cross-section measured in the fiducial phase-space as a function of the $p_{\textrm{T}}^{t\bar{t}}$ in the 4$\ell$ channel.

The total correlation matrix of the absolute particle-level differential cross-section measured in the fiducial phase-space as a function of the number of jets in the 4$\ell$ channel.

The total correlation matrix of the normalised particle-level differential cross-section measured in the fiducial phase-space as a function of the $p_{\textrm{T}}$ of the $Z$ boson in the 3$\ell$+4$\ell$ channel.

The total correlation matrix of the normalised particle-level differential cross-section measured in the fiducial phase-space as a function of the absolute value of rapidity of the $Z$ boson in the 3$\ell$+4$\ell$ channel.

The total correlation matrix of the normalised particle-level differential cross-section measured in the fiducial phase-space as a function of the $p_{\textrm{T}}^{l \textrm{non-}Z}$ in the 3$\ell$ channel.

The total correlation matrix of the normalised particle-level differential cross-section measured in the fiducial phase-space as a function of the $|\Delta \phi (Z, t_{\textrm{lep}})|/\pi$ in the 3$\ell$ channel.

The total correlation matrix of the normalised particle-level differential cross-section measured in the fiducial phase-space as a function of the $|\Delta y (Z, t_{\textrm{lep}})|/\pi$ in the 3$\ell$ channel.

The total correlation matrix of the normalised particle-level differential cross-section measured in the fiducial phase-space as a function of the number of jets in the 3$\ell$ channel.

The total correlation matrix of the normalised particle-level differential cross-section measured in the fiducial phase-space as a function of the $|\Delta \phi (l_{t}^{+}, l_{\bar{t}}^{-})|/\pi$ in the 4$\ell$ channel.

The total correlation matrix of the normalised particle-level differential cross-section measured in the fiducial phase-space as a function of the $|\Delta \phi (t\bar{t}, Z)|/\pi$ in the 4$\ell$ channel.

The total correlation matrix of the normalised particle-level differential cross-section measured in the fiducial phase-space as a function of the $p_{\textrm{T}}^{t\bar{t}}$ in the 4$\ell$ channel.

The total correlation matrix of the normalised particle-level differential cross-section measured in the fiducial phase-space as a function of the number of jets in the 4$\ell$ channel.

The $m_{\mathrm{eff}}$ distribution in SR0$_{\mathrm{bveto}}^{\mathrm{loose}}$ and SR0$_{\mathrm{bveto}}^{\mathrm{tight}}$ for events passing the signal region requirements except the $m_{\mathrm{eff}}$ requirement. Distributions for data, the estimated SM backgrounds after the background-only fit, and an example SUSY scenario are shown. "Other" is the sum of the $tWZ$, $t\bar{t}WW$, $t\bar{t} ZZ$, $t\bar{t} WH$, $t\bar{t} HH$, $t\bar{t} tW$, and $t\bar{t}t\bar{t}$ backgrounds. The last bin captures the overflow events. The lower panel shows the ratio of the observed data to the expected SM background yield in each bin. Both the statistical and systematic uncertainties in the SM background are included in the shaded band. The red arrows indicate the $m_{\mathrm{eff}}$ selections in the signal regions.

The $m_{\mathrm{eff}}$ distribution in SR1$_{\mathrm{bveto}}^{\mathrm{loose}}$ and SR1$_{\mathrm{bveto}}^{\mathrm{tight}}$ for events passing the signal region requirements except the $m_{\mathrm{eff}}$ requirement. Distributions for data, the estimated SM backgrounds after the background-only fit, and an example SUSY scenario are shown. "Other" is the sum of the $tWZ$, $t\bar{t}WW$, $t\bar{t} ZZ$, $t\bar{t} WH$, $t\bar{t} HH$, $t\bar{t} tW$, and $t\bar{t}t\bar{t}$ backgrounds. The last bin captures the overflow events. The lower panel shows the ratio of the observed data to the expected SM background yield in each bin. Both the statistical and systematic uncertainties in the SM background are included in the shaded band. The red arrows indicate the $m_{\mathrm{eff}}$ selections in the signal regions.

The $m_{\mathrm{eff}}$ distribution in SR2$_{\mathrm{bveto}}^{\mathrm{loose}}$ and SR2$_{\mathrm{bveto}}^{\mathrm{tight}}$ for events passing the signal region requirements except the $m_{\mathrm{eff}}$ requirement. Distributions for data, the estimated SM backgrounds after the background-only fit, and an example SUSY scenario are shown. "Other" is the sum of the $tWZ$, $t\bar{t}WW$, $t\bar{t} ZZ$, $t\bar{t} WH$, $t\bar{t} HH$, $t\bar{t} tW$, and $t\bar{t}t\bar{t}$ backgrounds. The last bin captures the overflow events. The lower panel shows the ratio of the observed data to the expected SM background yield in each bin. Both the statistical and systematic uncertainties in the SM background are included in the shaded band. The red arrows indicate the $m_{\mathrm{eff}}$ selections in the signal regions.

The $m_{\mathrm{eff}}$ distribution in SR0$_{\mathrm{breq}}$ for events passing the signal region requirements except the $m_{\mathrm{eff}}$ requirement. Distributions for data, the estimated SM backgrounds after the background-only fit, and an example SUSY scenario are shown. "Other" is the sum of the $tWZ$, $t\bar{t}WW$, $t\bar{t} ZZ$, $t\bar{t} WH$, $t\bar{t} HH$, $t\bar{t} tW$, and $t\bar{t}t\bar{t}$ backgrounds. The last bin captures the overflow events. The lower panel shows the ratio of the observed data to the expected SM background yield in each bin. Both the statistical and systematic uncertainties in the SM background are included in the shaded band. The red arrows indicate the $m_{\mathrm{eff}}$ selections in the signal regions.

The $m_{\mathrm{eff}}$ distribution in SR1$_{\mathrm{breq}}$ for events passing the signal region requirements except the $m_{\mathrm{eff}}$ requirement. Distributions for data, the estimated SM backgrounds after the background-only fit, and an example SUSY scenario are shown. "Other" is the sum of the $tWZ$, $t\bar{t}WW$, $t\bar{t} ZZ$, $t\bar{t} WH$, $t\bar{t} HH$, $t\bar{t} tW$, and $t\bar{t}t\bar{t}$ backgrounds. The last bin captures the overflow events. The lower panel shows the ratio of the observed data to the expected SM background yield in each bin. Both the statistical and systematic uncertainties in the SM background are included in the shaded band. The red arrows indicate the $m_{\mathrm{eff}}$ selections in the signal regions.

The $m_{\mathrm{eff}}$ distribution in SR2$_{\mathrm{breq}}$ for events passing the signal region requirements except the $m_{\mathrm{eff}}$ requirement. Distributions for data, the estimated SM backgrounds after the background-only fit, and an example SUSY scenario are shown. "Other" is the sum of the $tWZ$, $t\bar{t}WW$, $t\bar{t} ZZ$, $t\bar{t} WH$, $t\bar{t} HH$, $t\bar{t} tW$, and $t\bar{t}t\bar{t}$ backgrounds. The last bin captures the overflow events. The lower panel shows the ratio of the observed data to the expected SM background yield in each bin. Both the statistical and systematic uncertainties in the SM background are included in the shaded band. The red arrows indicate the $m_{\mathrm{eff}}$ selections in the signal regions.

Expected 95% CL exclusion limits on the higgsino GGM models. The limits are set using the statistical combination of disjoint signal regions. Where two (or more) signal regions overlap, the signal region contributing its observed $\mathrm{CL}_{\mathrm{s}}$ value to the combination is the one with the better (best) expected $\mathrm{CL}_{\mathrm{s}}$ value.

$+1\sigma$ expected 95% CL exclusion limits on the higgsino GGM models. The limits are set using the statistical combination of disjoint signal regions. Where two (or more) signal regions overlap, the signal region contributing its observed $\mathrm{CL}_{\mathrm{s}}$ value to the combination is the one with the better (best) expected $\mathrm{CL}_{\mathrm{s}}$ value.

$-1\sigma$ expected 95% CL exclusion limits on the higgsino GGM models. The limits are set using the statistical combination of disjoint signal regions. Where two (or more) signal regions overlap, the signal region contributing its observed $\mathrm{CL}_{\mathrm{s}}$ value to the combination is the one with the better (best) expected $\mathrm{CL}_{\mathrm{s}}$ value.

Observed 95% CL exclusion limits on the higgsino GGM models. The limits are set using the statistical combination of disjoint signal regions. Where two (or more) signal regions overlap, the signal region contributing its observed $\mathrm{CL}_{\mathrm{s}}$ value to the combination is the one with the better (best) expected $\mathrm{CL}_{\mathrm{s}}$ value.

$+1\sigma$ observed 95% CL exclusion limits on the higgsino GGM models. The limits are set using the statistical combination of disjoint signal regions. Where two (or more) signal regions overlap, the signal region contributing its observed $\mathrm{CL}_{\mathrm{s}}$ value to the combination is the one with the better (best) expected $\mathrm{CL}_{\mathrm{s}}$ value.

$-1\sigma$ observed 95% CL exclusion limits on the higgsino GGM models. The limits are set using the statistical combination of disjoint signal regions. Where two (or more) signal regions overlap, the signal region contributing its observed $\mathrm{CL}_{\mathrm{s}}$ value to the combination is the one with the better (best) expected $\mathrm{CL}_{\mathrm{s}}$ value.

Expected 95% CL exclusion limits on wino NLSP pair production with RPV LSP decays via $\lambda_{12k}$, where $k \in{1,2}$. The limits are set using the statistical combination of disjoint signal regions. Where two (or more) signal regions overlap, the signal region contributing its observed $\mathrm{CL}_{\mathrm{s}}$ value to the combination is the one with the better (best) expected $\mathrm{CL}_{\mathrm{s}}$ value.

$+1\sigma$ expected 95% CL exclusion limits on wino NLSP pair production with RPV LSP decays via $\lambda_{12k}$, where $k \in{1,2}$. The limits are set using the statistical combination of disjoint signal regions. Where two (or more) signal regions overlap, the signal region contributing its observed $\mathrm{CL}_{\mathrm{s}}$ value to the combination is the one with the better (best) expected $\mathrm{CL}_{\mathrm{s}}$ value.

$-1\sigma$ expected 95% CL exclusion limits on wino NLSP pair production with RPV LSP decays via $\lambda_{12k}$, where $k \in{1,2}$. The limits are set using the statistical combination of disjoint signal regions. Where two (or more) signal regions overlap, the signal region contributing its observed $\mathrm{CL}_{\mathrm{s}}$ value to the combination is the one with the better (best) expected $\mathrm{CL}_{\mathrm{s}}$ value.

Observed 95% CL exclusion limits on wino NLSP pair production with RPV LSP decays via $\lambda_{12k}$, where $k \in{1,2}$. The limits are set using the statistical combination of disjoint signal regions. Where two (or more) signal regions overlap, the signal region contributing its observed $\mathrm{CL}_{\mathrm{s}}$ value to the combination is the one with the better (best) expected $\mathrm{CL}_{\mathrm{s}}$ value.

$+1\sigma$ bserved 95% CL exclusion limits on wino NLSP pair production with RPV LSP decays via $\lambda_{12k}$, where $k \in{1,2}$. The limits are set using the statistical combination of disjoint signal regions. Where two (or more) signal regions overlap, the signal region contributing its observed $\mathrm{CL}_{\mathrm{s}}$ value to the combination is the one with the better (best) expected $\mathrm{CL}_{\mathrm{s}}$ value.

$-1\sigma$ observed 95% CL exclusion limits on wino NLSP pair production with RPV LSP decays via $\lambda_{12k}$, where $k \in{1,2}$. The limits are set using the statistical combination of disjoint signal regions. Where two (or more) signal regions overlap, the signal region contributing its observed $\mathrm{CL}_{\mathrm{s}}$ value to the combination is the one with the better (best) expected $\mathrm{CL}_{\mathrm{s}}$ value.

Expected 95% CL exclusion limits on wino NLSP pair production with RPV LSP decays via $\lambda_{i33}$, where $i \in{1,2}$. The limits are set using the statistical combination of disjoint signal regions. Where two (or more) signal regions overlap, the signal region contributing its observed $\mathrm{CL}_{\mathrm{s}}$ value to the combination is the one with the better (best) expected $\mathrm{CL}_{\mathrm{s}}$ value.

$+1\sigma$ expected 95% CL exclusion limits on wino NLSP pair production with RPV LSP decays via $\lambda_{i33}$, where $i \in{1,2}$. The limits are set using the statistical combination of disjoint signal regions. Where two (or more) signal regions overlap, the signal region contributing its observed $\mathrm{CL}_{\mathrm{s}}$ value to the combination is the one with the better (best) expected $\mathrm{CL}_{\mathrm{s}}$ value.

$-1\sigma$ expected 95% CL exclusion limits on wino NLSP pair production with RPV LSP decays via $\lambda_{i33}$, where $i \in{1,2}$. The limits are set using the statistical combination of disjoint signal regions. Where two (or more) signal regions overlap, the signal region contributing its observed $\mathrm{CL}_{\mathrm{s}}$ value to the combination is the one with the better (best) expected $\mathrm{CL}_{\mathrm{s}}$ value.

Observed 95% CL exclusion limits on wino NLSP pair production with RPV LSP decays via $\lambda_{i33}$, where $i \in{1,2}$. The limits are set using the statistical combination of disjoint signal regions. Where two (or more) signal regions overlap, the signal region contributing its observed $\mathrm{CL}_{\mathrm{s}}$ value to the combination is the one with the better (best) expected $\mathrm{CL}_{\mathrm{s}}$ value.

$+1\sigma$ observed 95% CL exclusion limits on wino NLSP pair production with RPV LSP decays via $\lambda_{i33}$, where $i \in{1,2}$. The limits are set using the statistical combination of disjoint signal regions. Where two (or more) signal regions overlap, the signal region contributing its observed $\mathrm{CL}_{\mathrm{s}}$ value to the combination is the one with the better (best) expected $\mathrm{CL}_{\mathrm{s}}$ value.

$-1\sigma$ observed 95% CL exclusion limits on wino NLSP pair production with RPV LSP decays via $\lambda_{i33}$, where $i \in{1,2}$. The limits are set using the statistical combination of disjoint signal regions. Where two (or more) signal regions overlap, the signal region contributing its observed $\mathrm{CL}_{\mathrm{s}}$ value to the combination is the one with the better (best) expected $\mathrm{CL}_{\mathrm{s}}$ value.

Expected 95% CL exclusion limits on slepton/sneutrino NLSP pair production with RPV LSP decays via $\lambda_{12k}$, where $k \in{1,2}$. The limits are set using the statistical combination of disjoint signal regions. Where two (or more) signal regions overlap, the signal region contributing its observed $\mathrm{CL}_{\mathrm{s}}$ value to the combination is the one with the better (best) expected $\mathrm{CL}_{\mathrm{s}}$ value.

$+1\sigma$ expected 95% CL exclusion limits on slepton/sneutrino NLSP pair production with RPV LSP decays via $\lambda_{12k}$, where $k \in{1,2}$. The limits are set using the statistical combination of disjoint signal regions. Where two (or more) signal regions overlap, the signal region contributing its observed $\mathrm{CL}_{\mathrm{s}}$ value to the combination is the one with the better (best) expected $\mathrm{CL}_{\mathrm{s}}$ value.

$-1\sigma$ expected 95% CL exclusion limits on slepton/sneutrino NLSP pair production with RPV LSP decays via $\lambda_{12k}$, where $k \in{1,2}$. The limits are set using the statistical combination of disjoint signal regions. Where two (or more) signal regions overlap, the signal region contributing its observed $\mathrm{CL}_{\mathrm{s}}$ value to the combination is the one with the better (best) expected $\mathrm{CL}_{\mathrm{s}}$ value.

Observed 95% CL exclusion limits on slepton/sneutrino NLSP pair production with RPV LSP decays via $\lambda_{12k}$, where $k \in{1,2}$. The limits are set using the statistical combination of disjoint signal regions. Where two (or more) signal regions overlap, the signal region contributing its observed $\mathrm{CL}_{\mathrm{s}}$ value to the combination is the one with the better (best) expected $\mathrm{CL}_{\mathrm{s}}$ value.

$+1\sigma$ observed 95% CL exclusion limits on slepton/sneutrino NLSP pair production with RPV LSP decays via $\lambda_{12k}$, where $k \in{1,2}$. The limits are set using the statistical combination of disjoint signal regions. Where two (or more) signal regions overlap, the signal region contributing its observed $\mathrm{CL}_{\mathrm{s}}$ value to the combination is the one with the better (best) expected $\mathrm{CL}_{\mathrm{s}}$ value.

$-1\sigma$ observed 95% CL exclusion limits on slepton/sneutrino NLSP pair production with RPV LSP decays via $\lambda_{12k}$, where $k \in{1,2}$. The limits are set using the statistical combination of disjoint signal regions. Where two (or more) signal regions overlap, the signal region contributing its observed $\mathrm{CL}_{\mathrm{s}}$ value to the combination is the one with the better (best) expected $\mathrm{CL}_{\mathrm{s}}$ value.

Expected 95% CL exclusion limits on slepton/sneutrino NLSP pair production with RPV LSP decays via $\lambda_{i33}$, where $i \in{1,2}$. The limits are set using the statistical combination of disjoint signal regions. Where two (or more) signal regions overlap, the signal region contributing its observed $\mathrm{CL}_{\mathrm{s}}$ value to the combination is the one with the better (best) expected $\mathrm{CL}_{\mathrm{s}}$ value.

$+1\sigma$ expected 95% CL exclusion limits on slepton/sneutrino NLSP pair production with RPV LSP decays via $\lambda_{i33}$, where $i \in{1,2}$. The limits are set using the statistical combination of disjoint signal regions. Where two (or more) signal regions overlap, the signal region contributing its observed $\mathrm{CL}_{\mathrm{s}}$ value to the combination is the one with the better (best) expected $\mathrm{CL}_{\mathrm{s}}$ value.

$-1\sigma$ expected 95% CL exclusion limits on slepton/sneutrino NLSP pair production with RPV LSP decays via $\lambda_{i33}$, where $i \in{1,2}$. The limits are set using the statistical combination of disjoint signal regions. Where two (or more) signal regions overlap, the signal region contributing its observed $\mathrm{CL}_{\mathrm{s}}$ value to the combination is the one with the better (best) expected $\mathrm{CL}_{\mathrm{s}}$ value.

Observed 95% CL exclusion limits on slepton/sneutrino NLSP pair production with RPV LSP decays via $\lambda_{i33}$, where $i \in{1,2}$. The limits are set using the statistical combination of disjoint signal regions. Where two (or more) signal regions overlap, the signal region contributing its observed $\mathrm{CL}_{\mathrm{s}}$ value to the combination is the one with the better (best) expected $\mathrm{CL}_{\mathrm{s}}$ value.

$+1\sigma$ observed 95% CL exclusion limits on slepton/sneutrino NLSP pair production with RPV LSP decays via $\lambda_{i33}$, where $i \in{1,2}$. The limits are set using the statistical combination of disjoint signal regions. Where two (or more) signal regions overlap, the signal region contributing its observed $\mathrm{CL}_{\mathrm{s}}$ value to the combination is the one with the better (best) expected $\mathrm{CL}_{\mathrm{s}}$ value.

$-1\sigma$ observed 95% CL exclusion limits on slepton/sneutrino NLSP pair production with RPV LSP decays via $\lambda_{i33}$, where $i \in{1,2}$. The limits are set using the statistical combination of disjoint signal regions. Where two (or more) signal regions overlap, the signal region contributing its observed $\mathrm{CL}_{\mathrm{s}}$ value to the combination is the one with the better (best) expected $\mathrm{CL}_{\mathrm{s}}$ value.

Expected 95% CL exclusion limits on gluino NLSP pair production with RPV LSP decays via $\lambda_{12k}$, where $k \in{1,2}$. The limits are set using the statistical combination of disjoint signal regions. Where two (or more) signal regions overlap, the signal region contributing its observed $\mathrm{CL}_{\mathrm{s}}$ value to the combination is the one with the better (best) expected $\mathrm{CL}_{\mathrm{s}}$ value.

$+1\sigma$ expected 95% CL exclusion limits on gluino NLSP pair production with RPV LSP decays via $\lambda_{12k}$, where $k \in{1,2}$. The limits are set using the statistical combination of disjoint signal regions. Where two (or more) signal regions overlap, the signal region contributing its observed $\mathrm{CL}_{\mathrm{s}}$ value to the combination is the one with the better (best) expected $\mathrm{CL}_{\mathrm{s}}$ value.

$-1\sigma$ expected 95% CL exclusion limits on gluino NLSP pair production with RPV LSP decays via $\lambda_{12k}$, where $k \in{1,2}$. The limits are set using the statistical combination of disjoint signal regions. Where two (or more) signal regions overlap, the signal region contributing its observed $\mathrm{CL}_{\mathrm{s}}$ value to the combination is the one with the better (best) expected $\mathrm{CL}_{\mathrm{s}}$ value.

Observed 95% CL exclusion limits on gluino NLSP pair production with RPV LSP decays via $\lambda_{12k}$, where $k \in{1,2}$. The limits are set using the statistical combination of disjoint signal regions. Where two (or more) signal regions overlap, the signal region contributing its observed $\mathrm{CL}_{\mathrm{s}}$ value to the combination is the one with the better (best) expected $\mathrm{CL}_{\mathrm{s}}$ value.

$+1\sigma$ observed 95% CL exclusion limits on gluino NLSP pair production with RPV LSP decays via $\lambda_{12k}$, where $k \in{1,2}$. The limits are set using the statistical combination of disjoint signal regions. Where two (or more) signal regions overlap, the signal region contributing its observed $\mathrm{CL}_{\mathrm{s}}$ value to the combination is the one with the better (best) expected $\mathrm{CL}_{\mathrm{s}}$ value.

$-1\sigma$ observed 95% CL exclusion limits on gluino NLSP pair production with RPV LSP decays via $\lambda_{12k}$, where $k \in{1,2}$. The limits are set using the statistical combination of disjoint signal regions. Where two (or more) signal regions overlap, the signal region contributing its observed $\mathrm{CL}_{\mathrm{s}}$ value to the combination is the one with the better (best) expected $\mathrm{CL}_{\mathrm{s}}$ value.

Expected 95% CL exclusion limits on gluino NLSP pair production with RPV LSP decays via $\lambda_{i33}$, where $i \in{1,2}$. The limits are set using the statistical combination of disjoint signal regions. Where two (or more) signal regions overlap, the signal region contributing its observed $\mathrm{CL}_{\mathrm{s}}$ value to the combination is the one with the better (best) expected $\mathrm{CL}_{\mathrm{s}}$ value.

$+1\sigma$ expected 95% CL exclusion limits on gluino NLSP pair production with RPV LSP decays via $\lambda_{i33}$, where $i \in{1,2}$. The limits are set using the statistical combination of disjoint signal regions. Where two (or more) signal regions overlap, the signal region contributing its observed $\mathrm{CL}_{\mathrm{s}}$ value to the combination is the one with the better (best) expected $\mathrm{CL}_{\mathrm{s}}$ value.

$-1\sigma$ expected 95% CL exclusion limits on gluino NLSP pair production with RPV LSP decays via $\lambda_{i33}$, where $i \in{1,2}$. The limits are set using the statistical combination of disjoint signal regions. Where two (or more) signal regions overlap, the signal region contributing its observed $\mathrm{CL}_{\mathrm{s}}$ value to the combination is the one with the better (best) expected $\mathrm{CL}_{\mathrm{s}}$ value.

Observed 95% CL exclusion limits on gluino NLSP pair production with RPV LSP decays via $\lambda_{i33}$, where $i \in{1,2}$. The limits are set using the statistical combination of disjoint signal regions. Where two (or more) signal regions overlap, the signal region contributing its observed $\mathrm{CL}_{\mathrm{s}}$ value to the combination is the one with the better (best) expected $\mathrm{CL}_{\mathrm{s}}$ value.

$+1\sigma$ observed 95% CL exclusion limits on gluino NLSP pair production with RPV LSP decays via $\lambda_{i33}$, where $i \in{1,2}$. The limits are set using the statistical combination of disjoint signal regions. Where two (or more) signal regions overlap, the signal region contributing its observed $\mathrm{CL}_{\mathrm{s}}$ value to the combination is the one with the better (best) expected $\mathrm{CL}_{\mathrm{s}}$ value.

$-1\sigma$ observed 95% CL exclusion limits on gluino NLSP pair production with RPV LSP decays via $\lambda_{i33}$, where $i \in{1,2}$. The limits are set using the statistical combination of disjoint signal regions. Where two (or more) signal regions overlap, the signal region contributing its observed $\mathrm{CL}_{\mathrm{s}}$ value to the combination is the one with the better (best) expected $\mathrm{CL}_{\mathrm{s}}$ value.

Observed upper limit on the signal cross section in fb for the wino NLSP models with RPV LSP decays via $\lambda_{12k}$ where $k \in{1,2}$. The interpolation between signal scenarios studied is included for illustration purposes only and may be subject to interpolation effects in sparsely populated areas.

Observed upper limit on the signal cross section in fb for the wino NLSP models with RPV LSP decays via $\lambda_{i33}$ where $i \in{1,2}$. The interpolation between signal scenarios studied is included for illustration purposes only and may be subject to interpolation effects in sparsely populated areas.

Observed upper limit on the signal cross section in fb for the slepton/sneutrino NLSP models with RPV LSP decays via $\lambda_{12k}$ where $k \in{1,2}$. The interpolation between signal scenarios studied is included for illustration purposes only and may be subject to interpolation effects in sparsely populated areas.

Observed upper limit on the signal cross section in fb for the slepton/sneutrino NLSP models with RPV LSP decays via $\lambda_{i33}$ where $i \in{1,2}$. The interpolation between signal scenarios studied is included for illustration purposes only and may be subject to interpolation effects in sparsely populated areas.

Observed upper limit on the signal cross section in fb for the gluino NLSP models with RPV LSP decays via $\lambda_{12k}$ where $k \in{1,2}$. The interpolation between signal scenarios studied is included for illustration purposes only and may be subject to interpolation effects in sparsely populated areas.

Observed upper limit on the signal cross section in fb for the gluino NLSP models with RPV LSP decays via $\lambda_{i33}$ where $i \in{1,2}$. The interpolation between signal scenarios studied is included for illustration purposes only and may be subject to interpolation effects in sparsely populated areas.

Observed upper limit on the signal cross section in fb for the higgsino GGM models. The interpolation between signal scenarios studied is included for illustration purposes only and may be subject to interpolation effects in sparsely populated areas.

Best expected SR for the wino NLSP models with RPV LSP decays via $\lambda_{12k}$ where $k \in{1,2}$. A value of 1 corresponds to SR0$_{\mathrm{bveto}}^{\mathrm{loose}}$SR1$_{\mathrm{bveto}}^{\mathrm{loose}}$SR2$_{\mathrm{bveto}}^{\mathrm{loose}}$, 2 corresponds to SR0$_{\mathrm{bveto}}^{\mathrm{loose}}$SR1$_{\mathrm{bveto}}^{\mathrm{loose}}$SR2$_{\mathrm{bveto}}^{\mathrm{tight}}$, 3 corresponds to SR0$_{\mathrm{bveto}}^{\mathrm{loose}}$SR1$_{\mathrm{bveto}}^{\mathrm{tight}}$SR2$_{\mathrm{bveto}}^{\mathrm{loose}}$, 4 corresponds to SR0$_{\mathrm{bveto}}^{\mathrm{loose}}$SR1$_{\mathrm{bveto}}^{\mathrm{tight}}$SR2$_{\mathrm{bveto}}^{\mathrm{tight}}$, and 5 corresponds to SR0$_{\mathrm{bveto}}^{\mathrm{tight}}$SR1$_{\mathrm{bveto}}^{\mathrm{tight}}$SR2$_{\mathrm{bveto}}^{\mathrm{tight}}$.

Best expected SR for the wino NLSP models with RPV LSP decays via $\lambda_{i33}$ where $i \in{1,2}$. A value of 1 corresponds to SR0$_{\mathrm{bveto}}^{\mathrm{loose}}$SR1$_{\mathrm{bveto}}^{\mathrm{loose}}$SR2$_{\mathrm{bveto}}^{\mathrm{loose}}$, 2 corresponds to SR0$_{\mathrm{bveto}}^{\mathrm{loose}}$SR1$_{\mathrm{bveto}}^{\mathrm{loose}}$SR2$_{\mathrm{bveto}}^{\mathrm{tight}}$, 3 corresponds to SR0$_{\mathrm{bveto}}^{\mathrm{loose}}$SR1$_{\mathrm{bveto}}^{\mathrm{tight}}$SR2$_{\mathrm{bveto}}^{\mathrm{loose}}$, 4 corresponds to SR0$_{\mathrm{bveto}}^{\mathrm{loose}}$SR1$_{\mathrm{bveto}}^{\mathrm{tight}}$SR2$_{\mathrm{bveto}}^{\mathrm{tight}}$, and 5 corresponds to SR0$_{\mathrm{bveto}}^{\mathrm{tight}}$SR1$_{\mathrm{bveto}}^{\mathrm{tight}}$SR2$_{\mathrm{bveto}}^{\mathrm{tight}}$.

Best expected SR for the slepton/sneutrino NLSP models with RPV LSP decays via $\lambda_{12k}$ where $k \in{1,2}$. A value of 1 corresponds to SR0$_{\mathrm{bveto}}^{\mathrm{loose}}$SR1$_{\mathrm{bveto}}^{\mathrm{loose}}$SR2$_{\mathrm{bveto}}^{\mathrm{loose}}$, 2 corresponds to SR0$_{\mathrm{bveto}}^{\mathrm{loose}}$SR1$_{\mathrm{bveto}}^{\mathrm{loose}}$SR2$_{\mathrm{bveto}}^{\mathrm{tight}}$, 3 corresponds to SR0$_{\mathrm{bveto}}^{\mathrm{loose}}$SR1$_{\mathrm{bveto}}^{\mathrm{tight}}$SR2$_{\mathrm{bveto}}^{\mathrm{loose}}$, 4 corresponds to SR0$_{\mathrm{bveto}}^{\mathrm{loose}}$SR1$_{\mathrm{bveto}}^{\mathrm{tight}}$SR2$_{\mathrm{bveto}}^{\mathrm{tight}}$, and 5 corresponds to SR0$_{\mathrm{bveto}}^{\mathrm{tight}}$SR1$_{\mathrm{bveto}}^{\mathrm{tight}}$SR2$_{\mathrm{bveto}}^{\mathrm{tight}}$.

Best expected SR for the slepton/sneutrino NLSP models with RPV LSP decays via $\lambda_{i33}$ where $i \in{1,2}$. A value of 1 corresponds to SR0$_{\mathrm{bveto}}^{\mathrm{loose}}$SR1$_{\mathrm{bveto}}^{\mathrm{loose}}$SR2$_{\mathrm{bveto}}^{\mathrm{loose}}$, 2 corresponds to SR0$_{\mathrm{bveto}}^{\mathrm{loose}}$SR1$_{\mathrm{bveto}}^{\mathrm{loose}}$SR2$_{\mathrm{bveto}}^{\mathrm{tight}}$, 3 corresponds to SR0$_{\mathrm{bveto}}^{\mathrm{loose}}$SR1$_{\mathrm{bveto}}^{\mathrm{tight}}$SR2$_{\mathrm{bveto}}^{\mathrm{loose}}$, 4 corresponds to SR0$_{\mathrm{bveto}}^{\mathrm{loose}}$SR1$_{\mathrm{bveto}}^{\mathrm{tight}}$SR2$_{\mathrm{bveto}}^{\mathrm{tight}}$, and 5 corresponds to SR0$_{\mathrm{bveto}}^{\mathrm{tight}}$SR1$_{\mathrm{bveto}}^{\mathrm{tight}}$SR2$_{\mathrm{bveto}}^{\mathrm{tight}}$.

Best expected SR for the gluino NLSP models with RPV LSP decays via $\lambda_{12k}$ where $k \in{1,2}$. A value of 1 corresponds to SR0$_{\mathrm{bveto}}^{\mathrm{loose}}$SR1$_{\mathrm{bveto}}^{\mathrm{loose}}$SR2$_{\mathrm{bveto}}^{\mathrm{loose}}$, 2 corresponds to SR0$_{\mathrm{bveto}}^{\mathrm{loose}}$SR1$_{\mathrm{bveto}}^{\mathrm{loose}}$SR2$_{\mathrm{bveto}}^{\mathrm{tight}}$, 3 corresponds to SR0$_{\mathrm{bveto}}^{\mathrm{loose}}$SR1$_{\mathrm{bveto}}^{\mathrm{tight}}$SR2$_{\mathrm{bveto}}^{\mathrm{loose}}$, 4 corresponds to SR0$_{\mathrm{bveto}}^{\mathrm{loose}}$SR1$_{\mathrm{bveto}}^{\mathrm{tight}}$SR2$_{\mathrm{bveto}}^{\mathrm{tight}}$, and 5 corresponds to SR0$_{\mathrm{bveto}}^{\mathrm{tight}}$SR1$_{\mathrm{bveto}}^{\mathrm{tight}}$SR2$_{\mathrm{bveto}}^{\mathrm{tight}}$.

Best expected SR for the gluino NLSP models with RPV LSP decays via $\lambda_{i33}$ where $i \in{1,2}$. A value of 1 corresponds to SR0$_{\mathrm{bveto}}^{\mathrm{loose}}$SR1$_{\mathrm{bveto}}^{\mathrm{loose}}$SR2$_{\mathrm{bveto}}^{\mathrm{loose}}$, 2 corresponds to SR0$_{\mathrm{bveto}}^{\mathrm{loose}}$SR1$_{\mathrm{bveto}}^{\mathrm{loose}}$SR2$_{\mathrm{bveto}}^{\mathrm{tight}}$, 3 corresponds to SR0$_{\mathrm{bveto}}^{\mathrm{loose}}$SR1$_{\mathrm{bveto}}^{\mathrm{tight}}$SR2$_{\mathrm{bveto}}^{\mathrm{loose}}$, 4 corresponds to SR0$_{\mathrm{bveto}}^{\mathrm{loose}}$SR1$_{\mathrm{bveto}}^{\mathrm{tight}}$SR2$_{\mathrm{bveto}}^{\mathrm{tight}}$, and 5 corresponds to SR0$_{\mathrm{bveto}}^{\mathrm{tight}}$SR1$_{\mathrm{bveto}}^{\mathrm{tight}}$SR2$_{\mathrm{bveto}}^{\mathrm{tight}}$.

Best expected SR for the higgsino GGM models. A value of 6 corresponds to SR0-ZZ$^{\mathrm{loose}}$, 7 corresponds to SR0-ZZ$^{\mathrm{tight}}$, 8 corresponds to SR0-ZZ$^{\mathrm{loose}}_{\mathrm{bveto}}$, and 9 corresponds to SR0-ZZ$^{\mathrm{tight}}_{\mathrm{bveto}}$.

Acceptance across the wino NLSP $\lambda_{12k}\neq 0$ models for SR0$_{\mathrm{bveto}}^{\mathrm{loose}}$. The interpolation between signal scenarios studied is included for illustration purposes only and may be subject to interpolation effects in sparsely populated areas.

Efficiency across the wino NLSP $\lambda_{12k}\neq 0$ models for SR0$_{\mathrm{bveto}}^{\mathrm{loose}}$. The interpolation between signal scenarios studied is included for illustration purposes only and may be subject to interpolation effects in sparsely populated areas.

Acceptance across the wino NLSP $\lambda_{12k}\neq 0$ models for SR0$_{\mathrm{bveto}}^{\mathrm{tight}}$. The interpolation between signal scenarios studied is included for illustration purposes only and may be subject to interpolation effects in sparsely populated areas.

Efficiency across the wino NLSP $\lambda_{12k}\neq 0$ models for SR0$_{\mathrm{bveto}}^{\mathrm{tight}}$. The interpolation between signal scenarios studied is included for illustration purposes only and may be subject to interpolation effects in sparsely populated areas.

Acceptance across the wino NLSP $\lambda_{12k}\neq 0$ models for SR0$_{\mathrm{breq}}$. The interpolation between signal scenarios studied is included for illustration purposes only and may be subject to interpolation effects in sparsely populated areas.

Efficiency across the wino NLSP $\lambda_{12k}\neq 0$ models for SR0$_{\mathrm{breq}}$. The interpolation between signal scenarios studied is included for illustration purposes only and may be subject to interpolation effects in sparsely populated areas.

Acceptance across the wino NLSP $\lambda_{i33}\neq 0$ models for SR1$_{\mathrm{bveto}}^{\mathrm{loose}}$. The interpolation between signal scenarios studied is included for illustration purposes only and may be subject to interpolation effects in sparsely populated areas.

Efficiency across the wino NLSP $\lambda_{i33}\neq 0$ models for SR1$_{\mathrm{bveto}}^{\mathrm{loose}}$. The interpolation between signal scenarios studied is included for illustration purposes only and may be subject to interpolation effects in sparsely populated areas.

Acceptance across the wino NLSP $\lambda_{i33}\neq 0$ models for SR1$_{\mathrm{bveto}}^{\mathrm{tight}}$. The interpolation between signal scenarios studied is included for illustration purposes only and may be subject to interpolation effects in sparsely populated areas.

Efficiency across the wino NLSP $\lambda_{i33}\neq 0$ models for SR1$_{\mathrm{bveto}}^{\mathrm{tight}}$. The interpolation between signal scenarios studied is included for illustration purposes only and may be subject to interpolation effects in sparsely populated areas.

Acceptance across the wino NLSP $\lambda_{i33}\neq 0$ models for SR1$_{\mathrm{breq}}$. The interpolation between signal scenarios studied is included for illustration purposes only and may be subject to interpolation effects in sparsely populated areas.

Efficiency across the wino NLSP $\lambda_{i33}\neq 0$ models for SR1$_{\mathrm{breq}}$. The interpolation between signal scenarios studied is included for illustration purposes only and may be subject to interpolation effects in sparsely populated areas.

Acceptance across the wino NLSP $\lambda_{i33}\neq 0$ models for SR2$_{\mathrm{bveto}}^{\mathrm{loose}}$. The interpolation between signal scenarios studied is included for illustration purposes only and may be subject to interpolation effects in sparsely populated areas.

Efficiency across the wino NLSP $\lambda_{i33}\neq 0$ models for SR2$_{\mathrm{bveto}}^{\mathrm{loose}}$. The interpolation between signal scenarios studied is included for illustration purposes only and may be subject to interpolation effects in sparsely populated areas.

Acceptance across the wino NLSP $\lambda_{i33}\neq 0$ models for SR2$_{\mathrm{bveto}}^{\mathrm{tight}}$. The interpolation between signal scenarios studied is included for illustration purposes only and may be subject to interpolation effects in sparsely populated areas.

Efficiency across the wino NLSP $\lambda_{i33}\neq 0$ models for SR2$_{\mathrm{bveto}}^{\mathrm{tight}}$. The interpolation between signal scenarios studied is included for illustration purposes only and may be subject to interpolation effects in sparsely populated areas.

Acceptance across the wino NLSP $\lambda_{i33}\neq 0$ models for SR2$_{\mathrm{breq}}$. The interpolation between signal scenarios studied is included for illustration purposes only and may be subject to interpolation effects in sparsely populated areas.

Efficiency across the wino NLSP $\lambda_{i33}\neq 0$ models for SR2$_{\mathrm{breq}}$. The interpolation between signal scenarios studied is included for illustration purposes only and may be subject to interpolation effects in sparsely populated areas.

Acceptance across the GGM Higgsino grid for SR0-ZZ$^{\mathrm{loose}}$. The interpolation between signal scenarios studied is included for illustration purposes only and may be subject to interpolation effects in sparsely populated areas.

Efficiency across the GGM Higgsino grid for SR0-ZZ$^{\mathrm{loose}}$. The interpolation between signal scenarios studied is included for illustration purposes only and may be subject to interpolation effects in sparsely populated areas.

Acceptance across the GGM Higgsino grid for SR0-ZZ$^{\mathrm{tight}}$. The interpolation between signal scenarios studied is included for illustration purposes only and may be subject to interpolation effects in sparsely populated areas.

Efficiency across the GGM Higgsino grid for SR0-ZZ$^{\mathrm{tight}}$. The interpolation between signal scenarios studied is included for illustration purposes only and may be subject to interpolation effects in sparsely populated areas.

Acceptance across the GGM Higgsino grid for SR0-ZZ$_{\mathrm{bveto}}^{\mathrm{loose}}$. The interpolation between signal scenarios studied is included for illustration purposes only and may be subject to interpolation effects in sparsely populated areas.

Efficiency across the GGM Higgsino grid for SR0-ZZ$_{\mathrm{bveto}}^{\mathrm{loose}}$. The interpolation between signal scenarios studied is included for illustration purposes only and may be subject to interpolation effects in sparsely populated areas.

Acceptance across the GGM Higgsino grid for SR0-ZZ$_{\mathrm{bveto}}^{\mathrm{tight}}$. The interpolation between signal scenarios studied is included for illustration purposes only and may be subject to interpolation effects in sparsely populated areas.

Efficiency across the GGM Higgsino grid for SR0-ZZ$_{\mathrm{bveto}}^{\mathrm{tight}}$. The interpolation between signal scenarios studied is included for illustration purposes only and may be subject to interpolation effects in sparsely populated areas.

The $p_{\mathrm{T}}$ of the light leptons in distribution in SR0$_{\mathrm{bveto}}^{\mathrm{loose}}$. Distributions for data, the estimated SM backgrounds after the background-only fit, and an example SUSY scenario are shown. "Other" is the sum of the $tWZ$, $t\bar{t}WW$, $t\bar{t} ZZ$, $t\bar{t} WH$, $t\bar{t} HH$, $t\bar{t} tW$, and $t\bar{t}t\bar{t}$ backgrounds. The last bin captures the overflow events. The lower panel shows the ratio of the observed data to the expected SM background yield in each bin. Both the statistical and systematic uncertainties in the SM background are included in the shaded band.

The $p_{\mathrm{T}}$ of the light leptons in distribution in SR0-ZZ$^{\mathrm{loose}}$. Distributions for data, the estimated SM backgrounds after the background-only fit, and an example SUSY scenario are shown. "Other" is the sum of the $tWZ$, $t\bar{t}WW$, $t\bar{t} ZZ$, $t\bar{t} WH$, $t\bar{t} HH$, $t\bar{t} tW$, and $t\bar{t}t\bar{t}$ backgrounds. The last bin captures the overflow events. The lower panel shows the ratio of the observed data to the expected SM background yield in each bin. Both the statistical and systematic uncertainties in the SM background are included in the shaded band.

The $p_{\mathrm{T}}$ of the light leptons in distribution in SR0-ZZ$^{\mathrm{tight}}$. Distributions for data, the estimated SM backgrounds after the background-only fit, and an example SUSY scenario are shown. "Other" is the sum of the $tWZ$, $t\bar{t}WW$, $t\bar{t} ZZ$, $t\bar{t} WH$, $t\bar{t} HH$, $t\bar{t} tW$, and $t\bar{t}t\bar{t}$ backgrounds. The last bin captures the overflow events. The lower panel shows the ratio of the observed data to the expected SM background yield in each bin. Both the statistical and systematic uncertainties in the SM background are included in the shaded band.

The $p_{\mathrm{T}}$ of the light leptons in distribution in SR0-ZZ$_{\mathrm{bveto}}^{\mathrm{loose}}$. Distributions for data, the estimated SM backgrounds after the background-only fit, and an example SUSY scenario are shown. "Other" is the sum of the $tWZ$, $t\bar{t}WW$, $t\bar{t} ZZ$, $t\bar{t} WH$, $t\bar{t} HH$, $t\bar{t} tW$, and $t\bar{t}t\bar{t}$ backgrounds. The last bin captures the overflow events. The lower panel shows the ratio of the observed data to the expected SM background yield in each bin. Both the statistical and systematic uncertainties in the SM background are included in the shaded band.

The $p_{\mathrm{T}}$ of the light leptons in distribution in SR5L. Distributions for data, the estimated SM backgrounds after the background-only fit, and an example SUSY scenario are shown. "Other" is the sum of the $tWZ$, $t\bar{t}WW$, $t\bar{t} ZZ$, $t\bar{t} WH$, $t\bar{t} HH$, $t\bar{t} tW$, and $t\bar{t}t\bar{t}$ backgrounds. The last bin captures the overflow events. The lower panel shows the ratio of the observed data to the expected SM background yield in each bin. Both the statistical and systematic uncertainties in the SM background are included in the shaded band.

The $p_{\mathrm{T}}$ of the light leptons in distribution in SR1$_{\mathrm{bveto}}^{\mathrm{loose}}$. Distributions for data, the estimated SM backgrounds after the background-only fit, and an example SUSY scenario are shown. "Other" is the sum of the $tWZ$, $t\bar{t}WW$, $t\bar{t} ZZ$, $t\bar{t} WH$, $t\bar{t} HH$, $t\bar{t} tW$, and $t\bar{t}t\bar{t}$ backgrounds. The last bin captures the overflow events. The lower panel shows the ratio of the observed data to the expected SM background yield in each bin. Both the statistical and systematic uncertainties in the SM background are included in the shaded band.

The $p_{\mathrm{T}}$ of the taus leptons in distribution in SR1$_{\mathrm{bveto}}^{\mathrm{loose}}$. Distributions for data, the estimated SM backgrounds after the background-only fit, and an example SUSY scenario are shown. "Other" is the sum of the $tWZ$, $t\bar{t}WW$, $t\bar{t} ZZ$, $t\bar{t} WH$, $t\bar{t} HH$, $t\bar{t} tW$, and $t\bar{t}t\bar{t}$ backgrounds. The last bin captures the overflow events. The lower panel shows the ratio of the observed data to the expected SM background yield in each bin. Both the statistical and systematic uncertainties in the SM background are included in the shaded band.

The $p_{\mathrm{T}}$ of the light taus in distribution in SR2$_{\mathrm{bveto}}^{\mathrm{loose}}$. Distributions for data, the estimated SM backgrounds after the background-only fit, and an example SUSY scenario are shown. "Other" is the sum of the $tWZ$, $t\bar{t}WW$, $t\bar{t} ZZ$, $t\bar{t} WH$, $t\bar{t} HH$, $t\bar{t} tW$, and $t\bar{t}t\bar{t}$ backgrounds. The last bin captures the overflow events. The lower panel shows the ratio of the observed data to the expected SM background yield in each bin. Both the statistical and systematic uncertainties in the SM background are included in the shaded band.

The lepton flavour and multiplicities in events with four light leptons and a Z veto. Distributions for data, the estimated SM backgrounds after the background-only fit, and an example SUSY scenario are shown. "Other" is the sum of the $tWZ$, $t\bar{t}WW$, $t\bar{t} ZZ$, $t\bar{t} WH$, $t\bar{t} HH$, $t\bar{t} tW$, and $t\bar{t}t\bar{t}$ backgrounds. The last bin captures the overflow events. The lower panel shows the ratio of the observed data to the expected SM background yield in each bin. Both the statistical and systematic uncertainties in the SM background are included in the shaded band.

The lepton flavour and multiplicities in events with four light leptons and one Z candidate. Distributions for data, the estimated SM backgrounds after the background-only fit, and an example SUSY scenario are shown. "Other" is the sum of the $tWZ$, $t\bar{t}WW$, $t\bar{t} ZZ$, $t\bar{t} WH$, $t\bar{t} HH$, $t\bar{t} tW$, and $t\bar{t}t\bar{t}$ backgrounds. The last bin captures the overflow events. The lower panel shows the ratio of the observed data to the expected SM background yield in each bin. Both the statistical and systematic uncertainties in the SM background are included in the shaded band.

The lepton flavour and multiplicities in events with four light leptons and two Z candidates. Distributions for data, the estimated SM backgrounds after the background-only fit, and an example SUSY scenario are shown. "Other" is the sum of the $tWZ$, $t\bar{t}WW$, $t\bar{t} ZZ$, $t\bar{t} WH$, $t\bar{t} HH$, $t\bar{t} tW$, and $t\bar{t}t\bar{t}$ backgrounds. The last bin captures the overflow events. The lower panel shows the ratio of the observed data to the expected SM background yield in each bin. Both the statistical and systematic uncertainties in the SM background are included in the shaded band.

The lepton flavour and multiplicities in events with exactly five light leptons. Distributions for data, the estimated SM backgrounds after the background-only fit, and an example SUSY scenario are shown. "Other" is the sum of the $tWZ$, $t\bar{t}WW$, $t\bar{t} ZZ$, $t\bar{t} WH$, $t\bar{t} HH$, $t\bar{t} tW$, and $t\bar{t}t\bar{t}$ backgrounds. The last bin captures the overflow events. The lower panel shows the ratio of the observed data to the expected SM background yield in each bin. Both the statistical and systematic uncertainties in the SM background are included in the shaded band.

The lepton flavour and multiplicities in events with three light leptons and one tau and a Z veto. Distributions for data, the estimated SM backgrounds after the background-only fit, and an example SUSY scenario are shown. "Other" is the sum of the $tWZ$, $t\bar{t}WW$, $t\bar{t} ZZ$, $t\bar{t} WH$, $t\bar{t} HH$, $t\bar{t} tW$, and $t\bar{t}t\bar{t}$ backgrounds. The last bin captures the overflow events. The lower panel shows the ratio of the observed data to the expected SM background yield in each bin. Both the statistical and systematic uncertainties in the SM background are included in the shaded band.

The lepton flavour and multiplicities in events with three light leptons and one tau and one Z candidate. Distributions for data, the estimated SM backgrounds after the background-only fit, and an example SUSY scenario are shown. "Other" is the sum of the $tWZ$, $t\bar{t}WW$, $t\bar{t} ZZ$, $t\bar{t} WH$, $t\bar{t} HH$, $t\bar{t} tW$, and $t\bar{t}t\bar{t}$ backgrounds. The last bin captures the overflow events. The lower panel shows the ratio of the observed data to the expected SM background yield in each bin. Both the statistical and systematic uncertainties in the SM background are included in the shaded band.

The lepton flavour and multiplicities in events with two light leptons and two taus and a Z veto. Distributions for data, the estimated SM backgrounds after the background-only fit, and an example SUSY scenario are shown. "Other" is the sum of the $tWZ$, $t\bar{t}WW$, $t\bar{t} ZZ$, $t\bar{t} WH$, $t\bar{t} HH$, $t\bar{t} tW$, and $t\bar{t}t\bar{t}$ backgrounds. The last bin captures the overflow events. The lower panel shows the ratio of the observed data to the expected SM background yield in each bin. Both the statistical and systematic uncertainties in the SM background are included in the shaded band.

The lepton flavour and multiplicities in events with two light leptons and two taus and one Z candidate. Distributions for data, the estimated SM backgrounds after the background-only fit, and an example SUSY scenario are shown. "Other" is the sum of the $tWZ$, $t\bar{t}WW$, $t\bar{t} ZZ$, $t\bar{t} WH$, $t\bar{t} HH$, $t\bar{t} tW$, and $t\bar{t}t\bar{t}$ backgrounds. The last bin captures the overflow events. The lower panel shows the ratio of the observed data to the expected SM background yield in each bin. Both the statistical and systematic uncertainties in the SM background are included in the shaded band.

Cutflow event yields in regions SR0$_{\mathrm{bveto}}^{\mathrm{loose}}$, SR0$_{\mathrm{bveto}}^{\mathrm{tight}}$, SR0$_{\mathrm{breq}}$, and SR5L for RPV models with the $\lambda_{12k}\neq 0$ coupling. All yields correspond to weighted events, so that effects from lepton reconstruction efficiencies, trigger corrections, pileup reweighting, etc., are included. They are normalized to the integrated luminosity of the data sample, $\int L dt = 139\,\mbox{fb\(^{-1}\)}$. The preliminary event reduction is a centralized stage where at least two electrons/muons with uncalibrated $p_{\mathrm{T}} >$ 9 GeV are required.

Cutflow event yields in regions SR1$_{\mathrm{bveto}}^{\mathrm{loose}}$, SR1$_{\mathrm{bveto}}^{\mathrm{tight}}$, and SR1$_{\mathrm{breq}}$ for RPV models with the $\lambda_{i33}\neq 0$ coupling. All yields correspond to weighted events, so that effects from lepton reconstruction efficiencies, trigger corrections, pileup reweighting, etc., are included. They are normalized to the integrated luminosity of the data sample, $\int L dt = 139\,\mbox{fb\(^{-1}\)}$. The preliminary event reduction is a centralized stage where at least two electrons/muons with uncalibrated $p_{\mathrm{T}} >$ 9 GeV are required.

Cutflow event yields in regions SR2$_{\mathrm{bveto}}^{\mathrm{loose}}$, SR2$_{\mathrm{bveto}}^{\mathrm{tight}}$, and SR2$_{\mathrm{breq}}$ for RPV models with the $\lambda_{i33}\neq 0$ coupling. All yields correspond to weighted events, so that effects from lepton reconstruction efficiencies, trigger corrections, pileup reweighting, etc., are included. They are normalized to the integrated luminosity of the data sample, $\int L dt = 139\,\mbox{fb\(^{-1}\)}$. The preliminary event reduction is a centralized stage where at least two electrons/muons with uncalibrated $p_{\mathrm{T}} >$ 9 GeV are required.

Cutflow event yields in regions SR0-ZZ$^{\mathrm{loose}}$, SR0-ZZ$^{\mathrm{tight}}$, SR0-ZZ$_{\mathrm{bveto}}^{\mathrm{loose}}$, SR0-ZZ$_{\mathrm{bveto}}^{\mathrm{tight}}$, and SR5L the higgsino GGM RPC model with BR($\tilde{\chi}^{0}_1 \rightarrow Z \tilde{G}$) = 50% and higgsino masses of 200 GeV, or BR($\tilde{\chi}^{0}_1 \rightarrow Z \tilde{G}$) = 100% and higgsino masses of 300 GeV. All yields correspond to weighted events, so that effects from lepton reconstruction efficiencies, trigger corrections, pileup reweighting, etc., are included. They are normalized to the integrated luminosity of the data sample, $\int L dt = 139\,\mbox{fb\(^{-1}\)}$. The generator filter is a selection of $\geq$4e/$\mu$/$\tau_{\mathrm{had-vis}}$ leptons with $p_{\mathrm{T}}(e,\mu)>4$GeV, $p_{\mathrm{T}}(\tau_{\mathrm{had-vis}})>15$GeV and $|\eta|<2.8$ and is applied during the MC generation of the simulated events. The preliminary event reduction is a centralized stage where at least two electrons/muons with uncalibrated $p_{\mathrm{T}} > 9$ GeV are required.