Inclusive differential cross section for two electrons and two muons (Max = 1710~GeV).

Mass of the first dilepton pair in the $Z$ mass region (Max = 88.1~GeV).

Mass of the first dilepton pair in the $H$ mass region (Max = 107~GeV).

Mass of the first dilepton pair in the offshell region (Max = 146~GeV).

Mass of the first dilepton pair in the $ZZ$ mass region (Max = 232~GeV).

Mass of the second dilepton pair in the $Z$ mass region (Max = 44.2~GeV).

Mass of the second dilepton pair in the $H$ mass region (Max = 59.6~GeV).

Mass of the second dilepton pair in the offhell region (Max = 88.3~GeV).

Mass of the second dilepton pair in the $ZZ$ mass region (Max = 1000~GeV).

$p_T$ of the first dilepton pair in the $Z$ mass region (Max = 291~GeV).

$p_T$ of the first dilepton pair in the $H$ mass region (Max = 245~GeV).

$p_T$ of the first dilepton pair in the offshell region (Max = 216~GeV).

$p_T$ of the first dilepton pair in the $ZZ$ mass region (Max = 555~GeV).

$p_T$ of the second dilepton pair in the $Z$ mass region (Max = 117~GeV).

$p_T$ of the second dilepton pair in the $H$ mass region (Max = 124~GeV).

$p_T$ of the second dilepton pair in the offhell region (Max = 192~GeV).

$p_T$ of the second dilepton pair in the $ZZ$ mass region (Max = 555~GeV).

$\cos\theta^*$ between first dilepton pair in the $Z$ mass region.

$\cos\theta^*$ between first dilepton pair in the offshell region.

$\cos\theta^*$ between first dilepton pair in the $H$ mass region.

$\cos\theta^*$ between first dilepton pair in the $ZZ$ mass region.

$\cos\theta^*$ between second dilepton pair in the $Z$ mass region.

$\cos\theta^*$ between second dilepton pair in the $H$ mass region.

$\cos\theta^*$ between second dilepton pair in the offshell mass region.

$\cos\theta^*$ between second dilepton pair in the $ZZ$ mass region.

$\Delta y$ between the dilepton pairs in the $Z$ mass region.

$\Delta y$ between the dilepton pairs in the $H$ mass region.

$\Delta y$ between the dilepton pairs in the offshell region.

$\Delta y$ between the dilepton pairs in the $ZZ$ mass region.

$\Delta \phi$ between the dilepton pairs in the $Z$ mass region.

$\Delta \phi$ between the dilepton pairs in the $H$ mass region.

$\Delta \phi$ between the dilepton pairs in the offshell region.

$\Delta \phi$ between the dilepton pairs in the $ZZ$ mass region.

$\Delta \phi$ between the leading lepton pairs in the $Z$ mass region.

$\Delta \phi$ between the leading lepton pairs in the $H$ mass region.

$\Delta \phi$ between the leding lepton pairs in the offshell region.

$\Delta \phi$ between the leading lepton pairs in the $ZZ$ mass region.

Four lepton mass in first $p_T^{4l}$ slice (Max = 726~GeV)

Four lepton mass in second $p_T^{4l}$ slice (Max = 1079~GeV)

Four lepton mass in third $p_T^{4l}$ slice (Max = 1153~GeV)

Four lepton mass in fourth $p_T^{4l}$ slice (Max = 1320~GeV)

Four lepton mass in fifth $p_T^{4l}$ slice (Max = 1413~GeV)

Four lepton mass in first $y^{4l}$ slice (Max = 1710~GeV)

Four lepton mass in second $y^{4l}$ slice (Max = 1116~GeV)

Four lepton mass in third $y^{4l}$ slice (Max = 1413~GeV)

Four lepton mass in fourth $y^{4l}$ slice (Max = 1044~GeV)

Four lepton mass in fifth $y^{4l}$ slice (Max = 918~GeV)

Covariance matrix for M4l

Covariance matrix for M4lvChannel

Covariance matrix for m12_vs_M4l

Covariance matrix for m34_vs_M4l

Covariance matrix for pt12_vs_M4l

Covariance matrix for pt34_vs_M4l

Covariance matrix for CTS12_vs_M4l

Covariance matrix for CTS34_vs_M4l

Covariance matrix for dYpairs_vs_M4l

Covariance matrix for dPhiPairs_vs_M4l

Covariance matrix for dPhiLeadLep_vs_M4l

Covariance matrix for M4lvPt4lbin

Covariance matrix for M4lvRapiditybin

EFT results (linear)

EFT results (full)

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.

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 number of data and fitted simulated events in each bin of the $|d_{0}^{\mu}|$ distribution in the $10<p_{\textrm{T}}^{\mu}<20$ GeV selection in the $\mu-\mu$ channel.

The number of data and fitted simulated events in each bin of the $|d_{0}^{\mu}|$ distribution in the $20<p_{\textrm{T}}^{\mu}<250$ GeV selection in the $e-\mu$ channel.

The number of data and fitted simulated events in each bin of the $|d_{0}^{\mu}|$ distribution in the $20<p_{\textrm{T}}^{\mu}<250$ GeV selection in the $\mu-\mu$ channel.

The measurement of the ratio of the rate of decay of W bosons to τ-leptons and muons, $R(\tau/\mu)=B(W\rightarrow\tau\nu_\tau)/B(W\rightarrow \mu\nu_\mu)$.

Acceptances on TRUTH level of the DMt samples in the tW1L channel for all bins in the SR. The acceptance is defined as the number of weighted TRUTH events in the SR over the number of expected events without any selections. The maps include all samples in the $m_H - tan\beta$ plane with $m_a = 250~GeV$.

Efficiencies of the DMt samples in the tW2L SR. The efficiency is defined as the number of weighted reconstructed events over the number of weighted TRUTH events in the SR. The maps include all samples in the $m_a - m_H$ plane with $tan\beta = 1$.

Acceptances on TRUTH level of the DMt samples in the tW2L SR. The acceptance is defined as the number of weighted TRUTH events in the SR over the number of expected events without any selections. The maps include all samples in the $m_a - m_H$ plane with $tan\beta = 1$.

Efficiencies of the DMt samples in the tW2L SR. The efficiency is defined as the number of weighted reconstructed events over the number of weighted TRUTH events in the SR. The maps include all samples in the $m_H - tan\beta$ plane with $m_a = 250~GeV$.

Acceptances on TRUTH level of the DMt samples in the tW2L SR. The acceptance is defined as the number of weighted TRUTH events in the SR over the number of expected events without any selections. The maps include all samples in the $m_H - tan\beta$ plane with $m_a = 250~GeV$.

Efficiencies of the DMt samples in the tj1L channel for all bins in the SR. The efficiency is defined as the number of weighted reconstructed events over the number of weighted TRUTH events in the SR. The map includes all used samples in the $m_H - tan\beta$ plane with $m_a = 250~GeV$.

Acceptances on TRUTH level of the DMt samples in the tj1L channel for all bins in the SR. The acceptance is defined as the number of weighted TRUTH events in the SR over the number of expected events without any selections. The map includes all used samples in the $m_H - tan\beta$ plane with $m_a = 250~GeV$.

Upper limits on signal strength (excluded cross section over theoretical cross section) of the tW1L analysis considering only DMt signal.

Upper limits on excluded cross sections of the tW1L analysis considering only the DMt signal.

The expected exclusion contours as a function of $(m_a, m_{H^{\pm}})$, assuming only $tW$+DM contributions, for the tW1L analysis channel.

The observed exclusion contours as a function of $(m_a, m_{H^{\pm}})$, assuming only $tW$+DM contributions, for the tW1L analysis channel.

Upper limits on signal strength (excluded cross section over theoretical cross section) of the tW1L analysis considering only DMt signal.

Upper limits on excluded cross sections of the tW1L analysis considering only the DMt signal.

The expected exclusion contours as a function of $(m_{H^{\pm}}, \tan\beta)$, assuming only $tW$+DM contributions, for the tW1L analysis channel.

The observed exclusion contours as a function of $(m_{H^{\pm}}, \tan\beta)$, assuming only $tW$+DM contributions, for the tW1L analysis channel.

Upper limits on signal strength (excluded cross section over theoretical cross section) of the tW2L analysis considering only DMt signal.

Upper limits on excluded cross sections of the tW2L analysis considering only the DMt signal.

The expected exclusion contours as a function of $(m_a, m_{H^{\pm}})$, assuming only $tW$+DM contributions, for the tW2L analysis channel.

The observed exclusion contours as a function of $(m_a, m_{H^{\pm}})$, assuming only $tW$+DM contributions, for the tW2L analysis channel.

Upper limits on signal strength (excluded cross section over theoretical cross section) of the tW2L analysis considering only DMt signal.

Upper limits on excluded cross sections of the tW2L analysis considering only the DMt signal.

The expected exclusion contours as a function of $(m_{H^{\pm}}, \tan\beta)$, assuming only $tW$+DM contributions, for the tW2L analysis channel.

The observed exclusion contours as a function of $(m_{H^{\pm}}, \tan\beta)$, assuming only $tW$+DM contributions, for the tW2L analysis channel.

Upper limits on signal strength (excluded cross section over theoretical cross section) of the combined tW1L and tW2L analyses considering only the DMt signal.

Upper limits on excluded cross sections of the combined tW1L and tW2L analyses considering only the DMt signal.

The expected exclusion contours as a function of $(m_a, m_{H^{\pm}})$, assuming only $tW$+DM contributions, for the statistical combination of the tW1L and tW2L analysis channel.

The observed exclusion contours as a function of $(m_a, m_{H^{\pm}})$, assuming only $tW$+DM contributions, for the statistical combination of the tW1L and tW2L analysis channel.

Upper limits on signal strength (excluded cross section over theoretical cross section) of the combined tW1L and tW2L analyses considering only the DMt signal.

Upper limits on excluded cross sections of the combined tW1L and tW2L analyses considering only the DMt signal.

The expected exclusion contours as a function of $(m_{H^{\pm}}, \tan\beta)$, assuming only $tW$+DM contributions, for the statistical combination of the tW1L and tW2L analysis channel.

The observed exclusion contours as a function of $(m_{H^{\pm}}, \tan\beta)$, assuming only $tW$+DM contributions, for the statistical combination of the tW1L and tW2L analysis channel.

Upper limits on signal strength (excluded cross section over theoretical cross section) of the tW1L analysis considering the DMt$\bar{t}$+DMt signal.

The expected exclusion contours as a function of $(m_a, m_{H^{\pm}})$, assuming DM$t\bar{t}$ and DM$t$ contributions, for the tW1L analysis channel.

The observed exclusion contours as a function of $(m_a, m_{H^{\pm}})$, assuming DM$t\bar{t}$ and DM$t$ contributions, for the tW1L analysis channel.

Upper limits on signal strength (excluded cross section over theoretical cross section) of the tW1L analysis considering the DMt$\bar{t}$+DMt signal.

The expected exclusion contours as a function of $(m_{H^{\pm}}, \tan\beta)$, assuming DM$t\bar{t}$ and DM$t$ contributions, for the tW1L analysis channel.

The observed exclusion contours as a function of $(m_{H^{\pm}}, \tan\beta)$, assuming DM$t\bar{t}$ and DM$t$ contributions, for the tW1L analysis channel.

Upper limits on signal strength (excluded cross section over theoretical cross section) of the tW2L analysis considering the DMt$\bar{t}$+DMt signal.

The expected exclusion contours as a function of $(m_a, m_{H^{\pm}})$, assuming DM$t\bar{t}$ and DM$t$ contributions, for the tW2L analysis channel.

The observed exclusion contours as a function of $(m_a, m_{H^{\pm}})$, assuming DM$t\bar{t}$ and DM$t$ contributions, for the tW2L analysis channel.

Upper limits on signal strength (excluded cross section over theoretical cross section) of the tW2L analysis considering the DMt$\bar{t}$+DMt signal.

The expected exclusion contours as a function of $(m_{H^{\pm}}, \tan\beta)$, assuming DM$t\bar{t}$ and DM$t$ contributions, for the tW2L analysis channel.

The observed exclusion contours as a function of $(m_{H^{\pm}}, \tan\beta)$, assuming DM$t\bar{t}$ and DM$t$ contributions, for the tW2L analysis channel.

Upper limits on signal strength (excluded cross section over theoretical cross section) of the combined tW1L and tW2L analyses considering the DMt$\bar{t}$+DMt signal.

The expected exclusion contours as a function of $(m_a, m_{H^{\pm}})$, assuming DM$t\bar{t}$ and DM$t$ contributions, for the statistical combination of the tW1L and tW2L analysis channel.

The observed exclusion contours as a function of $(m_a, m_{H^{\pm}})$, assuming DM$t\bar{t}$ and DM$t$ contributions, for the statistical combination of the tW1L and tW2L analysis channel.

Upper limits on signal strength (excluded cross section over theoretical cross section) of the combined tW1L and tW2L analyses considering the DMt$\bar{t}$+DMt signal.

The expected exclusion contours as a function of $(m_{H^{\pm}}, \tan\beta)$, assuming DM$t\bar{t}$ and DM$t$ contributions, for the statistical combination of the tW1L and tW2L analysis channel.

The observed exclusion contours as a function of $(m_{H^{\pm}}, \tan\beta)$, assuming DM$t\bar{t}$ and DM$t$ contributions, for the statistical combination of the tW1L and tW2L analysis channel.

Upper limits on signal strength (excluded cross section over theoretical cross section) of the tj1L analysis considering only the DMt signal.

Upper limits on upper limits on excluded cross sections of the tj1L analysis considering only the DMt signal.

The expected and observed cross section exclusion limits as a function of $m_{H^{\pm}}$ in the tj1L analysis channel for signal models with $m_a = 250~GeV$, and $\tan\beta=0.3$. The $\sigma^{}_\mathrm{BSM}$ is the cross section of the $t$-channel DM production process.

The expected and observed cross section exclusion limits as a function of $m_{H^{\pm}}$ in the tj1L analysis channel for signal models with $m_a = 250~GeV$, and $\tan\beta=0.5$. The $\sigma^{}_\mathrm{BSM}$ is the cross section of the $t$-channel DM production process.

Cross sections of the DMt samples in the tW1L channel. The maps include all samples in the $m_a - m_H$ plane with $tan\beta = 1$.

Cross sections of the DMt samples in the tW1L channel. The maps include all samples in the $m_H - tan\beta$ plane with $m_a = 250~GeV$.

Cross sections times branching ratio of the DMt samples in the tW2L channel. The maps include all samples in the $m_a - m_H$ plane with $tan\beta = 1$.

Cross sections times branching ratio of the DMt samples in the tW2L channel. The maps include all samples in the $m_H - tan\beta$ plane with $m_a = 250~GeV$.

Cross sections of the DMt samples in the tj1L channel. The map includes all samples in the $m_H - tan\beta$ plane with $m_a = 250~GeV$.

MC generator filter efficiencies of the DMt samples in the tW1L channel. The maps include all samples in the $m_a - m_H$ plane with $tan\beta = 1$.

MC generator filter efficiencies of the DMt samples in the tW1L channel. The maps include all samples in the $m_H - tan\beta$ plane with $m_a = 250~GeV$.

MC generator filter efficiencies of the DMt samples in the tW2L channel. The maps include all samples in the $m_a - m_H$ plane with $tan\beta = 1$.

MC generator filter efficiencies of the DMt samples in the tW2L channel. The maps include all samples in the $m_H - tan\beta$ plane with $m_a = 250~GeV$.

MC generator filter efficiencies of the DMt samples in the tj1L channel. The map includes all samples in the $m_H - tan\beta$ plane with $m_a = 250~GeV$.

Background-only fit results for the tW1L and tW2L signal regions. The backgrounds which contribute only a small amount (rare processes such as triboson, Higgs boson production processes, $t\bar{t}t\bar{t}$, $t\bar{t}WW$ and non-prompt or misidentified leptons background) are grouped and labelled as ``Others´´. The quoted uncertainties on the fitted SM background include both the statistical and systematic uncertainties.

Background-only fit results for the tj1L signal regions. The backgrounds which contribute only a small amount ($Z$+jets, rare processes such as $tWZ$, triboson, Higgs boson production processes, ,$t\bar{t}t\bar{t}$, $t\bar{t}WW$) are grouped and labelled as ``Others´´. The quoted uncertainties on the fitted SM background include both the statistical and systematic uncertainties.

Cutflow of the weighted events with statistical uncertainties for two DMt samples in all bins of the tW1L channel. The PreSelection includes at least 1 lepton in the event, at least 1 $b$-jet with $p_{\mathrm{T}} > 50~GeV$, $m\mathrm{_{T}^{lep}} > 30~GeV$, $\Delta\phi\mathrm{_{4jets, MET}^{min}} > 0.5$ and $E\mathrm{_{T}^{miss}} > 200~GeV$.

Cutflow of the weighted events with statistical uncertainties for two DMt samples in the tW2L channel. The PreSelection includes at least 2 leptons in the event, at least 1 $b$-jet with $p_{\mathrm{T}} > 40~GeV$, $m_{ll} > 40~GeV$, $m\mathrm{_{T2}} > 40~GeV$, $\Delta\phi\mathrm{_{4jets, MET}^{min}} > 0.5$ and $E\mathrm{_{T}^{miss}} > 200~GeV$.

Cutflow of the weighted events with the statistical uncertainties (except for the first cuts) for two DMt samples in all bins off the tj1L channel. The PreSelection includes at least 1 lepton in the event, at least 1 $b$-jet with $p_{\mathrm{T}} > 50~GeV$, $m\mathrm{_{T}^{lep}} > 30~GeV$, $\Delta\phi\mathrm{_{4jets, MET}^{min}} > 0.5$ and $E\mathrm{_{T}^{miss}} > 200~GeV$.

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.

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

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

Here THETA is the angle (in c.m. system) between the plane containing partons 1 and 3 and the plane containing partons 4 and 5: cos(theta) = [p1*p3]*[p4*p ]/(|[p1*p3]|*|[p4*p5]|). The estimated systematic uncertainty is 6 PCT.

Here M is scaled invariant mass of jet pairs. The estimated systematic uncertainty is 6 PCT.

Here M is scaled invariant mass of jet pairs. The estimated systematic uncertainty is 6 PCT.

Here M is scaled invariant mass of jet pairs. The estimated systematic uncertainty is 6 PCT.

The estimated systematic uncertainty is 6 PCT. 4-JET events in their center-of-mass system.

The estimated systematic uncertainty is 6 PCT. 4-JET events in their center-of-mass system.

The estimated systematic uncertainty is 6 PCT. 4-JET events in their center-of-mass system.

The estimated systematic uncertainty is 6 PCT. 4-JET events in their center-of-mass system.

Errors are statistical only. The estimated systematic uncertainty is 6 PCT. 4-JET events in their center-of-mass system.

Errors are statistical only. The estimated systematic uncertainty is 6 PCT. 4-JET events in their center-of-mass system.

Errors are statistical only. The estimated systematic uncertainty is 6 PCT. 4-JET events in their center-of-mass system.

Errors are statistical only. The estimated systematic uncertainty is 6 PCT. 4-JET events in their center-of-mass system.

Errors are statistical only. The estimated systematic uncertainty is 6 PCT. The measured distribution of cosine of space angles between pairs of jets for the 4-JET events in their center-of-mass system.

Errors are statistical only. The estimated systematic uncertainty is 6 PCT. The measured distribution of cosine of space angles between pairs of jets for the 4-JET events in their center-of-mass system.

Errors are statistical only. The estimated systematic uncertainty is 6 PCT. The measured distribution of cosine of space angles between pairs of jets for the 4-JET events in their center-of-mass system.

Errors are statistical only. The estimated systematic uncertainty is 6 PCT. The measured distribution of cosine of space angles between pairs of jets for the 4-JET events in their center-of-mass system.

Errors are statistical only. The estimated systematic uncertainty is 6 PCT. The measured distribution of cosine of space angles between pairs of jets for the 4-JET events in their center-of-mass system.

Errors are statistical only. The estimated systematic uncertainty is 6 PCT. The measured distribution of cosine of space angles between pairs of jets for the 4-JET events in their center-of-mass system.

Errors are statistical only. The estimated systematic uncertainty is 6 PCT. The measured distribution of scaled jets pair masses for the 4-JET events in their center-of-mass system.

Errors are statistical only. The estimated systematic uncertainty is 6 PCT. The measured distribution of scaled jets pair masses for the 4-JET events in their center-of-mass system.

Errors are statistical only. The estimated systematic uncertainty is 6 PCT. The measured distribution of scaled jets pair masses for the 4-JET events in their center-of-mass system.

Errors are statistical only. The estimated systematic uncertainty is 6 PCT. The measured distribution of scaled jets pair masses for the 4-JET events in their center-of-mass system.

Errors are statistical only. The estimated systematic uncertainty is 6 PCT. The measured distribution of scaled jets pair masses for the 4-JET events in their center-of-mass system.

Errors are statistical only. The estimated systematic uncertainty is 6 PCT. The measured distribution of scaled jets pair masses for the 4-JET events in their center-of-mass system.

Errors are statistical only. The estimated systematic uncertainty is 6 PCT. Here THETA is 'Bengtsson-Zerwas' (BZ) angle defined as COS(THETA(N=BZ)) = [p3*p4]*[p5*p6]/(|p3*p4|*|p5*p6|), where P(i) are jet momentum in the center-o f-mass frame.

Errors are statistical only. The estimated systematic uncertainty is 6 PCT. Here THETA is 'Nachtmann-Reiter' (NR) angle defined as COS(THETA(N=NR))= (p3-p4)*(p5-p6)/(|p3-p4|*|p5-p6|).