Probability of finding at least one TAR jet, where the p_T-leading TAR jet passes the m_Wcand and D₂^β=1 requirements, as a function of m_s. The probability is determined in a sample of signal events with m_Z'=2100 GeV, with the preselections applied.

Probability of finding a W_had candidate reconstructed as a pair of R=0.4 PFlow jets, as a function of m_s. The probability is determined in a sample of signal events with m_Z'=500 GeV, with the preselections applied that do not pass the requirements of the merged category.

Probability of finding a W_had candidate reconstructed as a pair of R=0.4 PFlow jets, as a function of m_s. The probability is determined in a sample of signal events with m_Z'=1000 GeV, with the preselections applied that do not pass the requirements of the merged category.

Probability of finding a W_had candidate reconstructed as a pair of R=0.4 PFlow jets, as a function of m_s. The probability is determined in a sample of signal events with m_Z'=1700 GeV, with the preselections applied that do not pass the requirements of the merged category.

Probability of finding a W_had candidate reconstructed as a pair of R=0.4 PFlow jets, as a function of m_s. The probability is determined in a sample of signal events with m_Z'=2100 GeV, with the preselections applied that do not pass the requirements of the merged category.

Observed exclusion contour at 95% C.L. for the dark Higgs model in the (m_Z', m_s) plane for g_q=0.25, g_χ=1, m_χ=200 GeV, and sinθ=0.01.

Expected exclusion contour at 95% C.L. for the dark Higgs model in the (m_Z', m_s) plane for g_q=0.25, g_χ=1, m_χ=200 GeV, and sinθ=0.01.

Expected+1σ exclusion contour at 95% C.L. for the dark Higgs model in the (m_Z', m_s) plane for g_q=0.25, g_χ=1, m_χ=200 GeV, and sinθ=0.01.

Expected-1σ exclusion contour at 95% C.L. for the dark Higgs model in the (m_Z', m_s) plane for g_q=0.25, g_χ=1, m_χ=200 GeV, and sinθ=0.01.

Expected+2σ exclusion contour at 95% C.L. for the dark Higgs model in the (m_Z', m_s) plane for g_q=0.25, g_χ=1, m_χ=200 GeV, and sinθ=0.01.

Expected-2σ exclusion contour at 95% C.L. for the dark Higgs model in the (m_Z', m_s) plane for g_q=0.25, g_χ=1, m_χ=200 GeV, and sinθ=0.01.

Observed upper limits at 95% C.L. on σ(pp → s χχ) × B(s → W^± W^∓) for m_Z'=0.5 TeV as a function of m_s. The expected limits, varied up and down by one and two standard deviations, are shown as green and yellow bands, respectively. The observed and expected limits are compared to the theoretical LO cross section for the σ(pp → s χχ) × B(s → W^± W^∓) process for m_Z'=0.5 TeV, shown in dashed blue.

Observed upper limits at 95% C.L. on σ(pp → s χχ) × B(s → W^± W^∓) for m_Z'=1 TeV as a function of m_s. The expected limits, varied up and down by one and two standard deviations, are shown as green and yellow bands, respectively. The observed and expected limits are compared to the theoretical LO cross section for the σ(pp → s χχ) × B(s → W^± W^∓) process for m_Z'=1 TeV, shown in dashed blue.

Observed upper limits at 95% C.L. on σ(pp → s χχ) × B(s → W^± W^∓) for m_Z'=1.7 TeV as a function of m_s. The expected limits, varied up and down by one and two standard deviations, are shown as green and yellow bands, respectively. The observed and expected limits are compared to the theoretical LO cross section for the σ(pp → s χχ) × B(s → W^± W^∓) process for m_Z'=1.7 TeV, shown in dashed blue.

Observed upper limits at 95% C.L. on σ(pp → s χχ) × B(s → W^± W^∓) for m_Z'=2.1 TeV as a function of m_s. The expected limits, varied up and down by one and two standard deviations, are shown as green and yellow bands, respectively. The observed and expected limits are compared to the theoretical LO cross section for the σ(pp → s χχ) × B(s → W^± W^∓) process for m_Z'=2.1 TeV, shown in dashed blue.

Data overlaid on SM background yields stacked in each SR and CR category after the fit to data ('Post-fit'). The yields in the SR are broken down into their contributions to the individual bins. The maximum-likelihood estimators are set to the conditional values of the CR-only fit, and propagated to SR and CRs.

Dominant sources of uncertainty for three dark Higgs scenarios after the fit to data. The uncertainties are quantified in terms of their contribution to the fitted signal uncertainty that is expressed relative to the theory prediction. Three representative dark Higgs signal scenarios with g_q=0.25, g_χ=1.0, sinθ=0.01 and m_χ=200 GeV are considered, which are indicated using the (m_Z', m_s) format in units of GeV in the table columns.

Cumulative efficiencies in the merged category for three representative dark Higgs signal scenarios with g_q=0.25, g_χ=1.0, sinθ=0.01, m_Z' = 1 TeV, and m_χ=200 GeV considering s→W(ℓν)W(qq) decays only.

Cumulative efficiencies in the resolved category for three representative dark Higgs signal scenarios with g_q=0.25, g_χ=1.0, sinθ=0.01, m_Z' = 1 TeV, and m_χ=200 GeV considering s→W(ℓν)W(qq) decays only.

Theoretical cross section for &sigma;(pp &rarr; s&chi;&chi;) &times; B(s &rarr; W^&plusmn;W^&#8723;) for each of the dark Higgs signal points at m_Z′ ={300, 350, 400, 500, 750} GeV, with g_q = 0.25, g<sub>&chi; = 1.0, sin&theta; = 0.01, m_Z′ = 1 TeV , and m_&chi; = 200 GeV. Also shown are experimentally excluded cross sections of &sigma;(pp &rarr; s&chi;&chi;) &times; B(s &rarr; W^&plusmn;W^&#8723;) (Obs.) together with the expectations (Exp.) varied up and down by one standard deviation (&plusmn;1&sigma;).

Theoretical cross section for &sigma;(pp &rarr; s&chi;&chi;) &times; B(s &rarr; W^&plusmn;W^&#8723;) for each of the dark Higgs signal points at m_Z′ ={1000, 1700} GeV, with g_q = 0.25, g<sub>&chi; = 1.0, sin&theta; = 0.01, m_Z′ = 1 TeV , and m_&chi; = 200 GeV. Also shown are experimentally excluded cross sections of &sigma;(pp &rarr; s&chi;&chi;) &times; B(s &rarr; W^&plusmn;W^&#8723;) (Obs.) together with the expectations (Exp.) varied up and down by one standard deviation (&plusmn;1&sigma;).

Theoretical cross section for &sigma;(pp &rarr; s&chi;&chi;) &times; B(s &rarr; W^&plusmn;W^&#8723;) for each of the dark Higgs signal points at m_Z′ ={2100, 2500, 2900, 3300} GeV, with g_q = 0.25, g<sub>&chi; = 1.0, sin&theta; = 0.01, m_Z′ = 1 TeV , and m_&chi; = 200 GeV. Also shown are experimentally excluded cross sections of &sigma;(pp &rarr; s&chi;&chi;) &times; B(s &rarr; W^&plusmn;W^&#8723;) (Obs.) together with the expectations (Exp.) varied up and down by one standard deviation (&plusmn;1&sigma;).

Expected 95% CL exclusion contours in the $(m(W_{R}), m(N_{R}))$ plane in the muon channel for boosted.

Observed 95% CL exclusion contours for the assumption of Dirac type in the $(m(W_{R}), m(N_{R}))$ plane in the electron channel for boosted.

Expected 95% CL exclusion contours for the assumption of Dirac type in the $(m(W_{R}), m(N_{R}))$ plane in the electron channel for boosted.

Observed 95% CL exclusion contours for the assumpition of Dirac type in the $(m(W_{R}), m(N_{R}))$ plane in the muon channel for boosted.

Expected 95% CL exclusion contours for the assumption of Dirac type in the $(m(W_{R}), m(N_{R}))$ plane in the muon channel for boosted.

Observed 95% CL exclusion contours for the assumption of Majorana type in the $(m(W_{R}), m(N_{R}))$ plane in the electron channel for resolved.

Expected 95% CL exclusion contours for the assumption of Majorana type in the $(m(W_{R}), m(N_{R}))$ plane in the electron channel for resolved.

Observed 95% CL exclusion contours for the assumption of Majorana type in the $(m(W_{R}), m(N_{R}))$ plane in the muon channel for resolved.

Expected 95% CL exclusion contours for the assumption of Majorana type in the $(m(W_{R}), m(N_{R}))$ plane in the muon channel for resolved.

Observed 95% CL exclusion contours for the assumption of Dirac type in the $(m(W_{R}), m(N_{R}))$ plane in the electron channel for resolved.

Expected 95% CL exclusion contours for the assumption of Dirac type in the $(m(W_{R}), m(N_{R}))$ plane in the electron channel for resolved.

Observed 95% CL exclusion contours for the assumption of Dirac type in the $(m(W_{R}), m(N_{R}))$ plane in the muon channel for resolved.

Expected 95% CL exclusion contours for the assumption of Dirac type in the $(m(W_{R}), m(N_{R}))$ plane in the muon channel for resolved.

The $m_{eJ}$ distribution in the boosted one electron channel.

The $m_{eeJ}$ distribution in the boosted two electron channel.

The $m_{\mu\mu J}$ distribution in the boosted two muon channel.

The $m_{jj}$ distribution in the resolved electron channel.

The $m_{jj}$ distribution in the resolved muon channel.

The $m_{eejj}$ distribution in the resolved electron channel.

The $m_{\mu\mu jj}$ distribution in the resolved muon channel.

The $h_{T}$ distribution in the resolved electron channel.

The $h_{T}$ distribution in the resolved muon channel.

Acceptance $\times$ efficiency for bSR1e.

Acceptance $\times$ efficiency for bSR2e.

Acceptance $\times$ efficiency for bSR2mu.

Acceptance $\times$ efficiency for rSROS2e.

Acceptance $\times$ efficiency for rSROS2mu.

Acceptance $\times$ efficiency for rSRSS2e.

Acceptance $\times$ efficiency for rSRSS2mu.

Expected and observed 95% CL upper limits for $\sigma(pp\to W_{R})\times B(W_{R}\to eeqq)$ as a function of $m(W_{R})$ for the special case of $m(W_{R}) = 2\times m(N_{R})$. The theoretical cross section corresponds to either the Dirac or Majorana neutrino interpretation by ignoring the charge of the leptons.

Expected and observed 95% CL upper limits for $\sigma(pp\to W_{R})\times B(W_{R} \to \mu\mu qq)$ as a function of $m(W_{R})$ for the special case of $m(W_{R}) = 2\times m(N_{R})$. The theoretical cross section corresponds to either the Dirac or Majorana neutrino interpretation by ignoring the charge of the leptons.

Expected and observed 95% CL upper limits for $\sigma(pp\to W_{R})\times B(W_{R} \to eeqq)$ as a function of $m(W_{R})$ for the special case of $m(N_{R})=200$ GeV. The theoretical cross section corresponds to either the Dirac or Majorana neutrino interpretation by ignoring the charge of the leptons.

Expected and observed 95% CL upper limits for $\sigma(pp\to W_{R})\times B(W_{R} \to \mu \mu qq)$ as a function of $m(W_{R})$ for the special case of $m(N_{R})=200$ GeV. The theoretical cross section corresponds to either the Dirac or Majorana neutrino interpretation by ignoring the charge of the leptons.

Expected and observed 95% CL upper limits for $\sigma(pp \to W_{R}) \times B(W_{R} \to eeqq)$ as a function of $m(W_{R})$ for the Majorana neutrino interpretations for the special case of $m(W_{R}) = 2\times m(N_{R})$.

Expected and observed 95% CL upper limits for $\sigma(pp \to W_{R}) \times B(W_{R} \to \mu\mu qq)$ as a function of $m(W_{R})$ for the Majorana neutrino interpretations for the special case of $m(W_{R}) = 2 \times m(N_{R})$.

Expected and observed 95% CL upper limits for $\sigma(pp \to W_{R}) \times B(W_{R} \to eeqq)$ as a function of $m(W_{R})$ for the Dirac neutrino interpretations for the special case of $m(W_{R}) = 2 \times m(N_{R})$.

Expected and observed 95% CL upper limits for $\sigma(pp \to W_{R}) \times B(W_{R} \to \mu \mu qq)$ as a function of $m(W_{R})$ for the Dirac neutrino interpretations for the special case of $m(W_{R}) = 2 \times m(N_{R})$.

Expected CLs values in the High mA SR, mA vs ma signal grid.

Observed CLs values in the High mA SR, mA vs ma signal grid.

Expected CLs values in the Low mA SR, mA vs ma signal grid.

Observed CLs values in the High mA SR, mA vs ma signal grid.

CLs+b values in the Low mA SR, mA vs tanB signal grid.

CLs+b values in the High mA SR, mA vs ma signal grid.

CLs+b values in the Low mA SR, mA vs ma signal grid.

Cut flow of the 2HDM+a signal points, gluon–gluon fusion production, Low mA SR. tanB = 1, $sin\theta$ = 0.35. The first two entries in the tables are number of raw MC events, third entry is theoretical prediction, and all other lines include the correct weights. Note that during the generation Higgs boson’s branching ratio to taus has been set to 1. An additional factor of 0.0627 is used to account for that, starting from the ‘Initial’ entry.

Cut flow of the 2HDM+a signal points, gluon–gluon fusion production, High mA SR. tanB = 1, $sin\theta$ = 0.35. The first two entries in the tables are number of raw MC events, third entry is theoretical prediction, and all other lines include the correct weights. Note that during the generation Higgs boson’s branching ratio to taus has been set to 1. An additional factor of 0.0627 is used to account for that, starting from the ‘Initial’ entry.

Cut flow of the 2HDM+a signal points, bb annihilation production, Low mA SR. tanB = 1, $sin\theta$ = 0.35. The first two entries in the tables are number of raw MC events, third entry is theoretical prediction, and all other lines include the correct weights. Note that during the generation Higgs boson’s branching ratio to taus has been set to 1. An additional factor of 0.0627 is used to account for that, starting from the ‘Initial’ entry.

Cut flow of the 2HDM+a signal points, bb annihilation production, High mA SR. tanB = 1, $sin\theta$ = 0.35. The first two entries in the tables are number of raw MC events, third entry is theoretical prediction, and all other lines include the correct weights. Note that during the generation Higgs boson’s branching ratio to taus has been set to 1. An additional factor of 0.0627 is used to account for that, starting from the ‘Initial’ entry.

A comparison of the observed and expected yields in the four bins of the Low mA SR.

A comparison of the observed and expected yields in the two bins of the High mA SR.

Expected exclusion contours at 95% CL as a function of mA and ma.

Observed exclusion contours at 95% CL as a function of mA and ma.

Expected +- 1sigma exclusion contours at 95% CL as a function of mA and ma.

Expected exclusion contours at 95% CL as a function of mA and ma.

Observed exclusion contours at 95% CL as a function of mA and ma.

Expected +- 1sigma exclusion contours at 95% CL as a function of mA and ma.

Acceptance, High mA SR, mA vs tanB grid, 400-750 GeV, bb prod

Acceptance, High mA SR, mA vs tanB grid, >750 GeV, bb prod

Acceptance, Low mA SR, mA vs tanB grid, 100-250 GeV, bb prod

Acceptance, Low mA SR, mA vs tanB grid, 250-400 GeV, bb prod

Acceptance, Low mA SR, mA vs tanB grid, 400-550 GeV, bb prod

Acceptance, Low mA SR, mA vs tanB grid, >550 GeV, bb prod

Acceptance, High mA SR, mA vs ma grid, 400-750 GeV, bb prod

Acceptance, High mA SR, mA vs ma grid, >750 GeV, bb prod

Acceptance, Low mA SR, mA vs ma grid, 100-250 GeV, bb prod

Acceptance, Low mA SR, mA vs ma grid, 250-400 GeV, bb prod

Acceptance, Low mA SR, mA vs ma grid, 400-550 GeV, bb prod

Acceptance, Low mA SR, mA vs ma grid, >550 GeV, bb prod

Acceptance, High mA SR, mA vs tanB grid, 400-750 GeV, ggF prod

Acceptance, High mA SR, mA vs tanB grid, >750 GeV, ggF prod

Acceptance, Low mA SR, mA vs tanB grid, 100-250 GeV, ggF prod

Acceptance, Low mA SR, mA vs tanB grid, 250-400 GeV, ggF prod

Acceptance, Low mA SR, mA vs tanB grid, 400-550 GeV, ggF prod

Acceptance, Low mA SR, mA vs tanB grid, >550 GeV, ggF prod

Acceptance, High mA SR, mA vs ma grid, 400-750 GeV, ggF prod

Acceptance, High mA SR, mA vs ma grid, >750 GeV, ggF prod

Acceptance, Low mA SR, mA vs ma grid, 100-250 GeV, ggF prod

Acceptance, Low mA SR, mA vs ma grid, 250-400 GeV, ggF prod

Acceptance, Low mA SR, mA vs ma grid, 400-550 GeV, ggF prod

Acceptance, Low mA SR, mA vs ma grid, >550 GeV, ggF prod

Efficiency, High mA SR, mA vs tanB grid, 400-750 GeV, bb prod

Efficiency, High mA SR, mA vs tanB grid, >750 GeV, bb prod

Efficiency, Low mA SR, mA vs tanB grid, 100-250 GeV, bb prod

Efficiency, Low mA SR, mA vs tanB grid, 250-400 GeV, bb prod

Efficiency, Low mA SR, mA vs tanB grid, 400-550 GeV, bb prod

Efficiency, Low mA SR, mA vs tanB grid, >550 GeV, bb prod

Efficiency, High mA SR, mA vs ma grid, 400-750 GeV, bb prod

Efficiency, High mA SR, mA vs ma grid, >750 GeV, bb prod

Efficiency, Low mA SR, mA vs ma grid, 100-250 GeV, bb prod

Efficiency, Low mA SR, mA vs ma grid, 250-400 GeV, bb prod

Efficiency, Low mA SR, mA vs ma grid, 400-550 GeV, bb prod

Efficiency, Low mA SR, mA vs ma grid, >550 GeV, bb prod

Efficiency, High mA SR, mA vs tanB grid, 400-750 GeV, ggF prod

Efficiency, High mA SR, mA vs tanB grid, >750 GeV, ggF prod

Efficiency, Low mA SR, mA vs tanB grid, 100-250 GeV, ggF prod

Efficiency, Low mA SR, mA vs tanB grid, 250-400 GeV, ggF prod

Efficiency, Low mA SR, mA vs tanB grid, 400-550 GeV, ggF prod

Efficiency, Low mA SR, mA vs tanB grid, >550 GeV, ggF prod

Efficiency, High mA SR, mA vs ma grid, 400-750 GeV, ggF prod

Efficiency, High mA SR, mA vs ma grid, >750 GeV, ggF prod

Efficiency, Low mA SR, mA vs ma grid, 100-250 GeV, ggF prod

Efficiency, Low mA SR, mA vs ma grid, 250-400 GeV, ggF prod

Efficiency, Low mA SR, mA vs ma grid, 400-550 GeV, ggF prod

Efficiency, Low mA SR, mA vs ma grid, >550 GeV, ggF prod

Figure 4a

Figure 4b

Figure 4c

Figure 5a

Figure 5b

Figure 5c

Figure 7

Figure 8a

Figure 8b

Figure 8c

Figure 9a

Figure 9b

Figure 9c

Figure 10a

Figure 10b

Figure 10c

$P_{\rm H}$ measurements are shown as a function of hyperon $p_{\rm T}$ at $\sqrt{s_{\rm NN}}$=19.6 and 27 GeV. Statistical uncertainties are represented as lines while systematic uncertainties are represented as boxes. There is no observed dependence of $P_{\rm H}$ on $p_{\rm T}$ at $\sqrt{s_{\rm NN}}$=19.6 or 27 GeV, consistent with previous observations.

$P_{\rm H}$ measurements are shown as a function of hyperon $y$ at $\sqrt{s_{\rm NN}}$=19.6 and 27 GeV. Statistical uncertainties are represented as lines while systematic uncertainties are represented as boxes. The data set at $\sqrt{s_{\rm NN}}$=19.6 GeV takes advantage of STAR upgrades to reach larger $|y|$.

Observed excluded cross-section at 95% CL for the WZ-mediated simplified model of wino $\tilde{\chi}^{\pm}_{1}/\tilde{\chi}^{0}_{2}$ production from Fig 8(aux).

Expected exclusion limits at 95% CL for the WZ-mediated simplified model of wino $\tilde{\chi}^{\pm}_{1}/\tilde{\chi}^{0}_{2}$ production.

Observed exclusion limits at 95% CL for the Wh-mediated simplified model of wino $\tilde{\chi}^{\pm}_{1}/\tilde{\chi}^{0}_{2}$ production from from Fig 13(a) and from Fig 7 and Fig 10(aux).

Observed excluded cross-section at 95% CL for the Wh-mediated simplified model of wino $\tilde{\chi}^{\pm}_{1}/\tilde{\chi}^{0}_{2}$ production from Fig 7(aux) and Fig 10(aux).

positive one $\sigma$ observed exclusion limits at 95% CL for the Wh-mediated simplified model of wino $\tilde{\chi}^{\pm}_{1}/\tilde{\chi}^{0}_{2}$ production from from Fig 13(a) and from Fig 7 and Fig 10(aux).

negative one $\sigma$ observed exclusion limits at 95% CL for the Wh-mediated simplified model of wino $\tilde{\chi}^{\pm}_{1}/\tilde{\chi}^{0}_{2}$ production from from Fig 13(a) and from Fig 7 and Fig 10(aux).

Expected exclusion limits at 95% CL for the Wh-mediated simplified model of wino $\tilde{\chi}^{\pm}_{1}/\tilde{\chi}^{0}_{2}$ production.

Expected exclusion limits at 95% CL for the Wh-mediated simplified model of wino $\tilde{\chi}^{\pm}_{1}/\tilde{\chi}^{0}_{2}$ production.

Expected exclusion limits at 95% CL for the Wh-mediated simplified model of wino $\tilde{\chi}^{\pm}_{1}/\tilde{\chi}^{0}_{2}$ production.

Number of signal events expected for 139 fb$^{-1}$ at different stages of the event selection for the signal region $SR^{bRPV}_{2l-SS}$. in a susy scenario where $\tilde{\chi}^{0} _{1} and \tilde{\chi}^{0} _{2}$ are produced in pairs and decay to all possible allowed bRPV decays. The masses of the superpartners involved in the process are set to $m(\tilde{\chi}^{0} _{1}/\tilde{\chi}^{0} _{2})$ = 200 GeV, tan$\beta$=5. Only statistical uncertainties are shown.

Number of signal events expected for 139 fb$^{-1}$ at different stages of the event selection for the signal region $SR^{bRPV}_{3l}$. in a susy scenario where $\tilde{\chi}^{0} _{1} and \tilde{\chi}^{0} _{2}$ are produced in pairs and decay to all possible allowed bRPV decays. The masses of the superpartners involved in the process are set to $m(\tilde{\chi}^{0} _{1}/\tilde{\chi}^{0} _{2})$ = 200 GeV, tan$\beta$=5. Only statistical uncertainties are shown.

Number of signal events expected for 139 fb$^{-1}$ at different stages of the event selection for the signal region $SR^{WZ}_{high-m_{T2}}$. The wino-like doublet pair ($\tilde{\chi}^{\pm} _{1} and \tilde{\chi}^{0} _{2}$) were produced and then decays into $bino-like \tilde{\chi}^{0} _{1}$ which is the lightest SUSY particle (LSP) accompanied by mass on-shell or mass off-shell W and Z bosons. The masses of the superpartners involved in the process are set to $m(\tilde{\chi}^{\pm} _{1}/\tilde{\chi}^{0} _{2})$ = 150 GeV, $m(\tilde{\chi}^{0} _{1})$ = 50 GeV. Only statistical uncertainties are shown.

Number of signal events expected for 139 fb$^{-1}$ at different stages of the event selection for the signal region $SR^{WZ}_{low-m_{T2}}$. The wino-like doublet pair ($\tilde{\chi}^{\pm} _{1} and \tilde{\chi}^{0} _{2}$) were produced and then decays into $bino-like \tilde{\chi}^{0} _{1}$ which is the lightest SUSY particle (LSP) accompanied by mass on-shell or mass off-shell W and Z bosons. The masses of the superpartners involved in the process are set to $m(\tilde{\chi}^{\pm} _{1}/\tilde{\chi}^{0} _{2})$ = 150 GeV, $m(\tilde{\chi}^{0} _{1})$ = 50 GeV. Only statistical uncertainties are shown.

Number of signal events expected for 139 fb$^{-1}$ at different stages of the event selection for the low mass $SR^{RPV}_{2l1b}$, where the $\tilde{\chi}^{0} _{1} and \tilde{\chi}^{0} _{2}$ are directly produced and undergoes prompt RPV decays. The masses of the superpartners involved in the process are set to $m(\tilde{\chi}^{0} _{1}/\tilde{\chi}^{0} _{2})$ = 200 GeV. Only statistical uncertainties are shown.

Number of signal events expected for 139 fb$^{-1}$ at different stages of the event selection for the medium mass $SR^{RPV}_{2l1b}$, where the $\tilde{\chi}^{0} _{1} and \tilde{\chi}^{0} _{2}$ are directly produced and undergoes prompt RPV decays. The masses of the superpartners involved in the process are set to $m(\tilde{\chi}^{0} _{1}/\tilde{\chi}^{0} _{2})$ = 200 GeV. Only statistical uncertainties are shown.

Number of signal events expected for 139 fb$^{-1}$ at different stages of the event selection for the low mass $SR^{RPV}_{2l2b}$, where the $\tilde{\chi}^{0} _{1} and \tilde{\chi}^{0} _{2}$ are directly produced and undergoes prompt RPV decays. The masses of the superpartners involved in the process are set to $m(\tilde{\chi}^{0} _{1}/\tilde{\chi}^{0} _{2})$ = 200 GeV. Only statistical uncertainties are shown.

Number of signal events expected for 139 fb$^{-1}$ at different stages of the event selection for the medium mass $SR^{RPV}_{2l2b}$, where the $\tilde{\chi}^{0} _{1} and \tilde{\chi}^{0} _{2}$ are directly produced and undergoes prompt RPV decays. The masses of the superpartners involved in the process are set to $m(\tilde{\chi}^{0} _{1}/\tilde{\chi}^{0} _{2})$ = 200 GeV. Only statistical uncertainties are shown.

Number of signal events expected for 139 fb$^{-1}$ at different stages of the event selection for the high mass $SR^{RPV}_{2l2b}$, where the $\tilde{\chi}^{0} _{1} and \tilde{\chi}^{0} _{2}$ are directly produced and undergoes prompt RPV decays. The masses of the superpartners involved in the process are set to $m(\tilde{\chi}^{0} _{1}/\tilde{\chi}^{0} _{2})$ = 200 GeV. Only statistical uncertainties are shown.

Number of signal events expected for 139 fb$^{-1}$ at different stages of the event selection for the low mass $SR^{RPV}_{2l3b}$, where the $\tilde{\chi}^{0} _{1} and \tilde{\chi}^{0} _{2}$ are directly produced and undergoes prompt RPV decays. The masses of the superpartners involved in the process are set to $m(\tilde{\chi}^{0} _{1}/\tilde{\chi}^{0} _{2})$ = 200 GeV. Only statistical uncertainties are shown.

Number of signal events expected for 139 fb$^{-1}$ at different stages of the event selection for the medium mass $SR^{RPV}_{2l3b}$, where the $\tilde{\chi}^{0} _{1} and \tilde{\chi}^{0} _{2}$ are directly produced and undergoes prompt RPV decays. The masses of the superpartners involved in the process are set to $m(\tilde{\chi}^{0} _{1}/\tilde{\chi}^{0} _{2})$ = 200 GeV. Only statistical uncertainties are shown.

Number of signal events expected for 139 fb$^{-1}$ at different stages of the event selection for the high mass $SR^{RPV}_{2l3b}$, where the $\tilde{\chi}^{0} _{1} and \tilde{\chi}^{0} _{2}$ are directly produced and undergoes prompt RPV decays. The masses of the superpartners involved in the process are set to $m(\tilde{\chi}^{0} _{1}/\tilde{\chi}^{0} _{2})$ = 200 GeV. Only statistical uncertainties are shown.

Number of signal events expected for 139 fb$^{-1}$ at different stages of the event selection for the $SR^{Wh}_{low-m_{T2} }$. The wino-like doublet pair ($\tilde{\chi}^{\pm} _{1} and \tilde{\chi}^{0} _{2}$) were produced and then decays into $bino-like \tilde{\chi}^{0} _{1}$ which is the lightest SUSY particle (LSP) accompanied by mass on-shell or mass off-shell W and Higgs bosons. The masses of the superpartners involved in the process are set to $m(\tilde{\chi}^{\pm} _{1}/\tilde{\chi}^{0} _{2})$ = 300 GeV, $m(\tilde{\chi}^{0} _{1})$ = 100 GeV. Only statistical uncertainties are shown.

Number of signal events expected for 139 fb$^{-1}$ at different stages of the event selection for the $SR^{Wh}_{high-m_{T2} }$. The wino-like doublet pair ($\tilde{\chi}^{\pm} _{1} and \tilde{\chi}^{0} _{2}$) were produced and then decays into $bino-like \tilde{\chi}^{0} _{1}$ which is the lightest SUSY particle (LSP) accompanied by mass on-shell or mass off-shell W and Higgs bosons. The masses of the superpartners involved in the process are set to $m(\tilde{\chi}^{\pm} _{1}/\tilde{\chi}^{0} _{2})$ = 300 GeV, $m(\tilde{\chi}^{0} _{1})$ = 100 GeV. Only statistical uncertainties are shown.

Signal Hepdataeptance for $SR^{bRPV}_{2l-SS}$ signal region from Fig 13(a)(aux) in a SUSY scenario where $\tilde{\chi}^{0} _{1} and \tilde{\chi}^{0} _{2}$ are produced in pairs and decay to all possible allowed bRPV decays.

Signal Hepdataeptance for $SR^{bRPV}_{3l}$ signal region from Fig 13(b)(aux) in a SUSY scenario where $\tilde{\chi}^{0} _{1} and \tilde{\chi}^{0} _{2}$ are produced in pairs and decay to all possible allowed bRPV decays.

Signal acceptance for $SR^{WZ}_{high-m_{T2}}$ in a SUSY scenario where the wino-like doublet pair ($\tilde{\chi}^{\pm} _{1} and \tilde{\chi}^{0} _{2}$) were produced and then decays into $bino-like \tilde{\chi}^{0} _{1}$ which is the lightest SUSY particle (LSP) accompanied by mass on-shell or mass off-shell W and Z bosons.

Signal acceptance for $SR^{WZ}_{low-m_{T2}}$ in a SUSY scenario where the wino-like doublet pair ($\tilde{\chi}^{\pm} _{1} and \tilde{\chi}^{0} _{2}$) were produced and then decays into $bino-like \tilde{\chi}^{0} _{1}$ which is the lightest SUSY particle (LSP) accompanied by mass on-shell or mass off-shell W and Z bosons.

Signal acceptance for $SR^{RPV}_{2l1b}-L$ signal region in a SUSY scenario where the $\tilde{\chi}^{0} _{1} and \tilde{\chi}^{0} _{2}$ are directly produced and undergoes prompt RPV decays.

Signal acceptance for $SR^{RPV}_{2l1b}-M$ signal region in a SUSY scenario where the $\tilde{\chi}^{0} _{1} and \tilde{\chi}^{0} _{2}$ are directly produced and undergoes prompt RPV decays.

Signal acceptance for $SR^{RPV}_{2l2b}-L$ signal region in a SUSY scenario where the $\tilde{\chi}^{0} _{1} and \tilde{\chi}^{0} _{2}$ are directly produced and undergoes prompt RPV decays.

Signal acceptance for $SR^{RPV}_{2l2b}-M$ signal region in a SUSY scenario where the $\tilde{\chi}^{0} _{1} and \tilde{\chi}^{0} _{2}$ are directly produced and undergoes prompt RPV decays.

Signal acceptance for $SR^{RPV}_{2l2b}-H$ signal region in a SUSY scenario where the $\tilde{\chi}^{0} _{1} and \tilde{\chi}^{0} _{2}$ are directly produced and undergoes prompt RPV decays.

Signal acceptance for $SR^{RPV}_{2l3b}-L$ signal region in a SUSY scenario where the $\tilde{\chi}^{0} _{1} and \tilde{\chi}^{0} _{2}$ are directly produced and undergoes prompt RPV decays.

Signal acceptance for $SR^{RPV}_{2l3b}-M$ signal region in a SUSY scenario where the $\tilde{\chi}^{0} _{1} and \tilde{\chi}^{0} _{2}$ are directly produced and undergoes prompt RPV decays.

Signal acceptance for $SR^{RPV}_{2l3b}-H$ signal region in a SUSY scenario where the $\tilde{\chi}^{0} _{1} and \tilde{\chi}^{0} _{2}$ are directly produced and undergoes prompt RPV decays.

Signal efficiency for $SR^{bRPV}_{2l-SS}$ signal region in a SUSY scenario where $\tilde{\chi}^{0} _{1} and \tilde{\chi}^{0} _{2}$ are produced in pairs and decay to all possible allowed bRPV decays.

Signal efficiency for $SR^{bRPV}_{3l}$ signal region in a SUSY scenario where $\tilde{\chi}^{0} _{1} and \tilde{\chi}^{0} _{2}$ are produced in pairs and decay to all possible allowed bRPV decays.

Signal efficiency for $SR^{WZ}_{high-m_{T2}}$ in a SUSY scenario where the wino-like doublet pair ($\tilde{\chi}^{\pm} _{1} and \tilde{\chi}^{0} _{2}$) were produced and then decays into $bino-like \tilde{\chi}^{0} _{1}$ which is the lightest SUSY particle (LSP) accompanied by mass on-shell or mass off-shell W and Z bosons.

Signal efficiency for $SR^{WZ}_{low-m_{T2}}$ in a SUSY scenario where the wino-like doublet pair ($\tilde{\chi}^{\pm} _{1} and \tilde{\chi}^{0} _{2}$) were produced and then decays into $bino-like \tilde{\chi}^{0} _{1}$ which is the lightest SUSY particle (LSP) accompanied by mass on-shell or mass off-shell W and Z bosons.

Signal efficiency for $SR^{RPV}_{2l1b}-L$ signal region in a SUSY scenario where the $\tilde{\chi}^{0} _{1} and \tilde{\chi}^{0} _{2}$ are directly produced and undergoes prompt RPV decays.

Signal efficiency for $SR^{RPV}_{2l1b}-M$ signal region in a SUSY scenario where the $\tilde{\chi}^{0} _{1} and \tilde{\chi}^{0} _{2}$ are directly produced and undergoes prompt RPV decays.

Signal efficiency for $SR^{RPV}_{2l2b}-L$ signal region in a SUSY scenario where the $\tilde{\chi}^{0} _{1} and \tilde{\chi}^{0} _{2}$ are directly produced and undergoes prompt RPV decays.

Signal efficiency for $SR^{RPV}_{2l2b}-M$ signal region in a SUSY scenario where the $\tilde{\chi}^{0} _{1} and \tilde{\chi}^{0} _{2}$ are directly produced and undergoes prompt RPV decays.

Signal efficiency for $SR^{RPV}_{2l2b}-H$ signal region in a SUSY scenario where the $\tilde{\chi}^{0} _{1} and \tilde{\chi}^{0} _{2}$ are directly produced and undergoes prompt RPV decays.

Signal efficiency for $SR^{RPV}_{2l3b}-L$ signal region in a SUSY scenario where the $\tilde{\chi}^{0} _{1} and \tilde{\chi}^{0} _{2}$ are directly produced and undergoes prompt RPV decays.

Signal efficiency for $SR^{RPV}_{2l3b}-M$ signal region in a SUSY scenario where the $\tilde{\chi}^{0} _{1} and \tilde{\chi}^{0} _{2}$ are directly produced and undergoes prompt RPV decays.

Signal efficiency for $SR^{RPV}_{2l3b}-H$ signal region in a SUSY scenario where the $\tilde{\chi}^{0} _{1} and \tilde{\chi}^{0} _{2}$ are directly produced and undergoes prompt RPV decays.

Signal acceptance for $SR^{Wh}_{high-m_{T2} }$ signal region from Fig 11(a)(aux) in a SUSY scenario where direct production of a lightest $\tilde{\chi}^{\pm} _{1} and \tilde{\chi}^{0} _{2}$ , decay with 100% branching ratio to a final state with a same sign light lepton (e or $\mu$) pair and two lightest neutralino1, via the on-shell emission of SM W and Higgs bosons,

Signal acceptance for $SR^{Wh}_{low-m_{T2} }$ signal region from Fig 11(b)(aux) in a SUSY scenario where direct production of a lightest $\tilde{\chi}^{\pm} _{1} and \tilde{\chi}^{0} _{2}$ , decay with 100% branching ratio to a final state with a same sign light lepton (e or $\mu$) pair and two lightest neutralino1, via the on-shell emission of SM W and Higgs bosons,

Signal efficiency for $SR^{Wh}_{high-m_{T2} }$ signal region from Fig 15(a)(aux) in a SUSY scenario where direct production of a lightest $\tilde{\chi}^{\pm} _{1} and \tilde{\chi}^{0} _{2}$ , decay with 100% branching ratio to a final state with a same sign light lepton (e or $\mu$) pair and two lightest neutralino1, via the on-shell emission of SM W and Higgs bosons,

Signal efficiency for $SR^{Wh}_{low-m_{T2} }$ signal region from Fig 15(b)(aux) in a SUSY scenario where direct production of a lightest $\tilde{\chi}^{\pm} _{1} and \tilde{\chi}^{0} _{2}$ , decay with 100% branching ratio to a final state with a same sign light lepton (e or $\mu$) pair and two lightest neutralino1, via the on-shell emission of SM W and Higgs bosons,

Observed 95% X-section upper limits as a function of higgsino $\tilde{\chi}^{\pm}_{1}/\tilde{\chi}^{0}_{1}/\tilde{\chi}^{0}_{2}$ mass in the bilinear RPV model from Fig 14.

Observed 95% X-section upper limits as a function of higgsino $\tilde{\chi}^{0}_{1}/\tilde{\chi}^{0}_{2}$ mass in the UDD RPV model from Fig 18.

Observed 95% X-section upper limits as a function of wino $\tilde{\chi}^{\pm}_{1}/\tilde{\chi}^{0}_{2}$ mass in the WZ-mediated simplified model of wino $\tilde{\chi}^{\pm}_{1}/\tilde{\chi}^{0}_{2}$ production from Fig 9(aux).

N-1 distributions for $m_{T2}$ of observed data and expected background towards $SR^{WZ}_{high-m_{T2}}$ from publication's Figure 11(a) . The last bin is inclusive.

N-1 distributions for $m_{T2}$ of observed data and expected background towards $SR^{WZ}_{low-m_{T2}}$ from publication's Figure 11(b) . The last bin is inclusive.

N-1 distributions for $m_{T2}$ of observed data and expected background towards $SR^{bRPV}_{2l-SS}$ from publication's Figure 11(c) . The last bin is inclusive.

N-1 distributions for $m_{T2}$ of observed data and expected background towards $SR^{bRPV}_{3l}$ from publication's Figure 11(d) . The last bin is inclusive.

N-1 distributions for $\sum p^{b-jet}_{T}/\sum p^{jet}_{T}$ of observed data and expected background towards $SR^{RPV}_{2l1b}-L$ from publication's Figure 16(a) . The last bin is inclusive.

N-1 distributions for $\sum p^{b-jet}_{T}/\sum p^{jet}_{T}$ of observed data and expected background towards $SR^{RPV}_{2l2b}-M$ from publication's Figure 16(b) . The last bin is inclusive.

N-1 distributions for $\sum p^{b-jet}_{T}/\sum p^{jet}_{T}$ of observed data and expected background towards $SR^{RPV}_{2l3b}-H$ from publication's Figure 16(c) . The last bin is inclusive.

N-1 distribution for $E_{T}^{miss}$ in $SR^{Wh}_{high-m_{T2} }$ in ee channel

N-1 distribution for $E_{T}^{miss}$ in $SR^{Wh}_{high-m_{T2} }$ in e$\mu$ channel

N-1 distribution for $E_{T}^{miss}$ in $SR^{Wh}_{high-m_{T2} }$ in $\mu\mu$ channel

N-1 distribution for $\mathcal{S}(E_{T}^{miss})$ in $SR^{Wh}_{low-m_{T2} }$ in ee channel

N-1 distribution for $\mathcal{S}(E_{T}^{miss})$ in $SR^{Wh}_{low-m_{T2} }$ in e$\mu$ channel

N-1 distribution for $\mathcal{S}(E_{T}^{miss})$ in $SR^{Wh}_{low-m_{T2} }$ in $\mu\mu$ channel

The values of $\langle \mathrm{d}N_\mathrm{ch} / \mathrm{d}\eta \rangle$ averaged over $|\eta|<0.5$ for the INEL, NSD, and INEL>0 event classes at $\sqrt{s} = 5.02$ TeV

Pseudorapidity density distributions of charged particles, $\mathrm{d}N_{\mathrm{ch}}/\mathrm{d}\eta$, in pp collisions at $\sqrt{s} = $ 5.02 TeV for the INEL>0 with momentum threshold $p_{T}>0.15 \mathrm{GeV}/c$ ($\mathrm{INEL>0}_{p_{\mathrm{T}}>0.15}^{|\eta|<0.8}$).

Pseudorapidity density distributions of charged particles, $\mathrm{d}N_{\mathrm{ch}}/\mathrm{d}\eta$, in pp collisions at $\sqrt{s} = $ 5.02 TeV for the INEL>0 with momentum threshold $p_{T}>0.5 \mathrm{GeV}/c$ ($\mathrm{INEL>0}_{p_{\mathrm{T}}>0.5}^{|\eta|<0.8}$).

Pseudorapidity density distributions of charged particles, $\mathrm{d}N_{\mathrm{ch}}/\mathrm{d}\eta$, in pp collisions at $\sqrt{s} = $ 5.02 TeV for the INEL>0 with momentum threshold $p_{T}>1 \mathrm{GeV}/c$ ($\mathrm{INEL>0}_{p_{\mathrm{T}}>1}^{|\eta|<0.8}$).

Pseudorapidity density distributions of charged particles, $\mathrm{d}N_{\mathrm{ch}}/\mathrm{d}\eta$, in pp collisions at $\sqrt{s} = $ 5.02 TeV for the INEL>0 with momentum threshold $p_{T}>2 \mathrm{GeV}/c$ ($\mathrm{INEL>0}_{p_{\mathrm{T}}>2}^{|\eta|<0.8}$).

Pseudorapidity density distributions of charged particles, $\mathrm{d}N_{\mathrm{ch}}/\mathrm{d}\eta$, in pp collisions at $\sqrt{s} = $ 13 TeV for the INEL>0 with momentum threshold $p_{T}>0.15 \mathrm{GeV}/c$ ($\mathrm{INEL>0}_{p_{\mathrm{T}}>0.15}^{|\eta|<0.8}$).

Pseudorapidity density distributions of charged particles, $\mathrm{d}N_{\mathrm{ch}}/\mathrm{d}\eta$, in pp collisions at $\sqrt{s} = $ 13 TeV for the INEL>0 with momentum threshold $p_{T}>0.5 \mathrm{GeV}/c$ ($\mathrm{INEL>0}_{p_{\mathrm{T}}>0.5}^{|\eta|<0.8}$).

Pseudorapidity density distributions of charged particles, $\mathrm{d}N_{\mathrm{ch}}/\mathrm{d}\eta$, in pp collisions at $\sqrt{s} = $ 13 TeV for the INEL>0 with momentum threshold $p_{T}>1 \mathrm{GeV}/c$ ($\mathrm{INEL>0}_{p_{\mathrm{T}}>1}^{|\eta|<0.8}$).

Pseudorapidity density distributions of charged particles, $\mathrm{d}N_{\mathrm{ch}}/\mathrm{d}\eta$, in pp collisions at $\sqrt{s} = $ 13 TeV for the INEL>0 with momentum threshold $p_{T}>2 \mathrm{GeV}/c$ ($\mathrm{INEL>0}_{p_{\mathrm{T}}>2}^{|\eta|<0.8}$).

HFE (electrons from semileptonic decays of heavy-flavor hadrons) invariant yields in different centrality intervals of Au+Au collisions at $\sqrt{s_{\rm NN}}=200$ GeV. The vertical bars and the boxes represent statistical and systematic uncertainties, respectively. The horizontal bars indicate the bin width.

HFE (electrons from semileptonic decays of heavy-flavor hadrons) $R_{\rm AA}$ (red circles) as a function of $p_{\rm T}$ in different centrality intervals of Au+Au collisions at $\sqrt{s_{\rm NN}}=200$ GeV, compared with STAR (yellow stars) and PHENIX (green squares) published results, and Duke ((modified Langevin transport model, blue line) and PHSD (parton-hadron-string dynamics model, orange line) model calculations. Vertical bars and boxes around data points represent combined statistical and systematic uncertainties from both Au+Au and $p$+$p$ measurements, respectively. Boxes at unity show the global uncertainties, which for this analysis include the 8% global uncertainty on $p$+$p$ reference and the $N_{\rm coll}$ uncertainties. The left box is for PHENIX and the right one for STAR.

HFE (electrons from semileptonic decays of heavy-flavor hadrons) $R_{\rm AA}$ (red circles) as a function of $N_{\rm part}$ in Au+Au collisions at $\sqrt{s_{\rm NN}}=200$ GeV, compared with PHENIX measurements (green squares), and Duke (blue line) and PHSD (orange line) model calculations. Vertical bars and boxes around data points represent statistical and systematic uncertainties from Au+Au measurements, respectively. The gray band represents the $N_{\rm coll}$ uncertainties. The boxes at unity show the global uncertainties including the total uncertainties of the $p$+$p$ reference.

Lévy scale parameter $R$ versus average $m_T$ of the pair. The graphical representation of statistical and systematic uncertainties is the same as in Fig. 4.

Lévy index parameter $\alpha$ versus average $m_T$ of the pair. Statistical and systematic uncertainties are indicated similarly to Fig. 4. The horizontal line, %\aplha$ = 1.207, represents the 0%-30% centrality average value of %\aplha$.

Inverse square of the Lévy scale parameter 1/$R^2$ versus average $m_T$ of the pair. Statistical and systematic uncertainties shown as bars and boxes, respectively.

New scale parameter $\hat{R}$ versus average $m_T$ of the pair, with a linear fit. Statistical and systematic uncertainties shown as bars and boxes, respectively.

Normalized correlation strength parameter $\lambda/\lambda_{max}$ versus average $m_T$ of the pair. The data are compared with parameter scans from Refs. [58, 59] using different values of in-medium $\eta^{\prime}$ mass $m^{*}_{\eta^{\prime}}$ and slope parameter $B^{-1}_{\eta^{\prime}}$ . A best fit with Eq. (63) and the resulting $H$ and $\sigma$ parameters are also shown.

Normalized correlation strength parameter $\lambda/\lambda_{max}$ versus average $m_T$ of the pair. The data are compared with parameter scans from Refs. [58, 59] using different values of in-medium $\eta^{\prime}$ mass $m^{*}_{\eta^{\prime}}$ and slope parameter $B^{-1}_{\eta^{\prime}}$ . A best fit with Eq. (63) and the resulting $H$ and $\sigma$ parameters are also shown.

Normalized correlation strength parameter $\lambda/\lambda_{max}$ versus average $m_T$ of the pair. The data are compared with parameter scans from Refs. [58, 59] using different values of in-medium $\eta^{\prime}$ mass $m^{*}_{\eta^{\prime}}$ and slope parameter $B^{-1}_{\eta^{\prime}}$ . A best fit with Eq. (63) and the resulting $H$ and $\sigma$ parameters are also shown.

Normalized correlation strength parameter $\lambda/\lambda_{max}$ versus average $m_T$ of the pair. The data are compared with parameter scans from Refs. [58, 59] using different values of in-medium $\eta^{\prime}$ mass $m^{*}_{\eta^{\prime}}$ and slope parameter $B^{-1}_{\eta^{\prime}}$ . A best fit with Eq. (63) and the resulting $H$ and $\sigma$ parameters are also shown.

Normalized correlation strength parameter $\lambda/\lambda_{max}$ versus average $m_T$ of the pair. The data are compared with parameter scans from Refs. [58, 59] using different values of in-medium $\eta^{\prime}$ mass $m^{*}_{\eta^{\prime}}$ and slope parameter $B^{-1}_{\eta^{\prime}}$ . A best fit with Eq. (63) and the resulting $H$ and $\sigma$ parameters are also shown.

Two dimensional observed likelihood profile scan as a function of mu_ggH and mu_qqH

Distributions of $\Delta\phi$ as obtained from simulation, requiring various $\textrm{t}$ tag multiplicities for the signal in two representative combinations of (gluino, neutralino) masses with large (2.2, 0.1)$\mathrm{TeV}$ and small (1.8, 1.3)$\mathrm{TeV}$ mass difference.

Results of fits to the $n_{\textrm{b}}$ multiplicity for control regions for the muon channel and with the requirements $3\leq n_{\textrm{jet}}\leq4$, $250<L_T<350~\mathrm{GeV}$, $500<H_T<750~\mathrm{GeV}$, $n_{\textrm{W}}\geq1$, $\Delta \phi<1$. The shaded area shows the fit uncertainty of the total background.

Results of fits to the $n_{\textrm{b}}$ multiplicity for control regions for the muon channel and with the requirements $3\leq n_{\textrm{jet}}\leq4$, $350<L_T<450~\mathrm{GeV}$, $H_T>1000~\mathrm{GeV}$, $n_{\textrm{W}}\geq0$, $\Delta \phi<1$. The shaded area shows the fit uncertainty of the total background.

Jet multiplicity distribution after the single-lepton baseline selection excluding the SRs for the multi-b analysis. The simulation is normalized to data with the SF mentioned in the plot.

Jet multiplicity distribution after the single-lepton baseline selection excluding the SRs for the the zero-b analysis (right). The simulation is normalized to data with the SF mentioned in the plot.

Jet multiplicity distribution in the dilepton CRs for the multi-b analysis. The simulation is normalized to data with the SF mentioned in the plot.

Jet multiplicity distribution in the dilepton CRs for the zero-b analysis. The simulation is normalized to data with the SF mentioned in the plot.

The double ratio of the single-lepton and dilepton ratio between data and simulation together with fit results and their uncertainties is shown for the multi-b analysis. The fits are performed for each data-taking year; 2018 is shown as an example.

The double ratio of the single-lepton and dilepton ratio between data and simulation together with fit results and their uncertainties is shown for the zero-b analysis. The fits are performed for each data-taking year; 2018 is shown as an example.

The prefit $L_{\mathrm{P}}$ distribution for selected electron candidates in the baseline QCD selection, with modified requirements of $n_{\text{jet}}\in[3,4]$ and $n_{\mathrm{b}}=0$.

The prefit $L_{\mathrm{P}}$ distribution for anti-selected electron candidates in the baseline QCD selection, with modified requirements of $n_{\text{jet}}\in[3,4]$ and $n_{\mathrm{b}}=0$.

Observed event yields in the MB SRs of the multi-b analysis compared to signal and background predictions. The relative fraction of the different SM EW background contributions determined in simulation is shown by the stacked, colored histograms, normalized so that their sum is equal to the background estimated using data control regions. The QCD background is predicted using the $L_p$ method. The signal is shown for two representative combinations of (gluino, neutralino) masses with large (2.2, 0.1) $\mathrm{TeV}$ and small (1.8, 1.3) $\mathrm{TeV}$ mass differences.

Observed event yields in the MB SRs of the zero-b analysis compared to signal and background predictions. The $\textrm{W}$+jets, $\textrm{t}\bar{\textrm{t}}$+jets, and QCD predictions are extracted from data control samples, while the other background contributions are estimated from simulation. The signal is shown for two representative combinations of (gluino, neutralino) masses with large (2.2, 0.1) $\mathrm{TeV}$ and small (1.8, 1.3) $\mathrm{TeV}$ mass differences.

Cross section limits at 95% CL for the T1tttt (left) model, as functions of the gluino and LSP masses, assuming a branching fraction of 100%. The mass of the intermediate chargino is taken to be halfway between the gluino and the neutralino masses. The solid black (dashed red) lines correspond to the observed (expected) mass limits, with the thicker lines representing the central values and the thinner lines representing the $\pm1\sigma$ uncertainty bands related to the theoretical (experimental) uncertainties.

Cross section limits at 95% CL for the T1tttt model, as functions of the gluino and LSP masses, assuming a branching fraction of 100%. The mass of the intermediate chargino is taken to be halfway between the gluino and the neutralino masses. The solid black (dashed red) lines correspond to the observed (expected) mass limits, with the thicker lines representing the central values and the thinner lines representing the $\pm1\sigma$ uncertainty bands related to the theoretical (experimental) uncertainties.

Cross section limits at 95% CL for the T1tttt model, as functions of the gluino and LSP masses, assuming a branching fraction of 100%. The mass of the intermediate chargino is taken to be halfway between the gluino and the neutralino masses. The solid black (dashed red) lines correspond to the observed (expected) mass limits, with the thicker lines representing the central values and the thinner lines representing the $\pm1\sigma$ uncertainty bands related to the theoretical (experimental) uncertainties.

Cross section limits at 95% CL for the T1tttt model, as functions of the gluino and LSP masses, assuming a branching fraction of 100%. The mass of the intermediate chargino is taken to be halfway between the gluino and the neutralino masses. The solid black (dashed red) lines correspond to the observed (expected) mass limits, with the thicker lines representing the central values and the thinner lines representing the $\pm1\sigma$ uncertainty bands related to the theoretical (experimental) uncertainties.

Cross section limits at 95% CL for the T1tttt model, as functions of the gluino and LSP masses, assuming a branching fraction of 100%. The mass of the intermediate chargino is taken to be halfway between the gluino and the neutralino masses. The solid black (dashed red) lines correspond to the experved (expected) mass limits, with the thicker lines representing the central values and the thinner lines representing the $\pm1\sigma$ uncertainty bands related to the theoretical (experimental) uncertainties.

Cross section limits at 95% CL for the T1tttt model, as functions of the gluino and LSP masses, assuming a branching fraction of 100%. The mass of the intermediate chargino is taken to be halfway between the gluino and the neutralino masses. The solid black (dashed red) lines correspond to the experved (expected) mass limits, with the thicker lines representing the central values and the thinner lines representing the $\pm1\sigma$ uncertainty bands related to the theoretical (experimental) uncertainties.

Cross section limits at 95% CL for the T1tttt model, as functions of the gluino and LSP masses, assuming a branching fraction of 100%. The mass of the intermediate chargino is taken to be halfway between the gluino and the neutralino masses. The solid black (dashed red) lines correspond to the experved (expected) mass limits, with the thicker lines representing the central values and the thinner lines representing the $\pm1\sigma$ uncertainty bands related to the theoretical (experimental) uncertainties.

Cross section limits at 95% CL for the T5qqqqWW model, as functions of the gluino and LSP masses, assuming a branching fraction of 100%. The mass of the intermediate chargino is taken to be halfway between the gluino and the neutralino masses. The solid black (dashed red) lines correspond to the observed (expected) mass limits, with the thicker lines representing the central values and the thinner lines representing the $\pm1\sigma$ uncertainty bands related to the theoretical (experimental) uncertainties.

Cross section limits at 95% CL for the T5qqqqWW model, as functions of the gluino and LSP masses, assuming a branching fraction of 100%. The mass of the intermediate chargino is taken to be halfway between the gluino and the neutralino masses. The solid black (dashed red) lines correspond to the observed (expected) mass limits, with the thicker lines representing the central values and the thinner lines representing the $\pm1\sigma$ uncertainty bands related to the theoretical (experimental) uncertainties.

Cross section limits at 95% CL for the T5qqqqWW model, as functions of the gluino and LSP masses, assuming a branching fraction of 100%. The mass of the intermediate chargino is taken to be halfway between the gluino and the neutralino masses. The solid black (dashed red) lines correspond to the observed (expected) mass limits, with the thicker lines representing the central values and the thinner lines representing the $\pm1\sigma$ uncertainty bands related to the theoretical (experimental) uncertainties.

Cross section limits at 95% CL for the T5qqqqWW model, as functions of the gluino and LSP masses, assuming a branching fraction of 100%. The mass of the intermediate chargino is taken to be halfway between the gluino and the neutralino masses. The solid black (dashed red) lines correspond to the observed (expected) mass limits, with the thicker lines representing the central values and the thinner lines representing the $\pm1\sigma$ uncertainty bands related to the theoretical (experimental) uncertainties.

Cross section limits at 95% CL for the T5qqqqWW model, as functions of the gluino and LSP masses, assuming a branching fraction of 100%. The mass of the intermediate chargino is taken to be halfway between the gluino and the neutralino masses. The solid black (dashed red) lines correspond to the experved (expected) mass limits, with the thicker lines representing the central values and the thinner lines representing the $\pm1\sigma$ uncertainty bands related to the theoretical (experimental) uncertainties.

Cross section limits at 95% CL for the T5qqqqWW model, as functions of the gluino and LSP masses, assuming a branching fraction of 100%. The mass of the intermediate chargino is taken to be halfway between the gluino and the neutralino masses. The solid black (dashed red) lines correspond to the experved (expected) mass limits, with the thicker lines representing the central values and the thinner lines representing the $\pm1\sigma$ uncertainty bands related to the theoretical (experimental) uncertainties.

Cross section limits at 95% CL for the T5qqqqWW model, as functions of the gluino and LSP masses, assuming a branching fraction of 100%. The mass of the intermediate chargino is taken to be halfway between the gluino and the neutralino masses. The solid black (dashed red) lines correspond to the experved (expected) mass limits, with the thicker lines representing the central values and the thinner lines representing the $\pm1\sigma$ uncertainty bands related to the theoretical (experimental) uncertainties.

Observed number of events in the MB SR bins of the multi-b analysis. All bins are defined with $\Delta\phi > 0.75$.

Observed number of events in the MB SR bins of the zero-b analysis.

The simplified likelihood can be constructed using this covariance matrix of nuisance parameters in the signal region (MB SR). The bin contents have been aggregated accross the three years to simplify the likelihood.

Selection efficiency for the multi-b signal for the SRs.

The simplified likelihood can be constructed using this covariance matrix of nuisance parameters in the signal region (MB SR). The bins were aggregated across HT to stabilize the simplified likelihood.

Selection efficiency for the zero-b signal for the SRs. The bins correspond to the aggregated bins that were used for the covariance matrix.

Dependence of $\overline{p}$-$\overline{d}$ correlation on pseudorapidity acceptance of $\overline{p}$ and $\overline{d}$ selection in Pb-Pb collisions at $\sqrt{s_\mathrm{NN}}$ = 5.02 TeV. Results are for 20.0--30.0$\%$ collision centrality.

Dependence of $\overline{p}$-$\overline{d}$ correlation on pseudorapidity acceptance of $\overline{p}$ and $\overline{d}$ selection in Pb-Pb collisions at $\sqrt{s_\mathrm{NN}}$ = 5.02 TeV. Results are for 50.0--60.0$\%$ collision centrality.

The obtained p-value from comparing data to background estimation across all bins, defined by windows in the X and Y particle masses, in the anomaly signal region.

The expected 95% CL limits on the cross-section $\sigma(pp \rightarrow Y \rightarrow XH \rightarrow q\bar{q}b\bar{b}$) in pb in the two-dimensional space of $m_Y$ versus $m_X$, obtained from a simultaneous fit of both merged and resolved two-prong signal regions with all statistical and systematic uncertainties.

The observed 95% CL limits on the cross-section $\sigma(pp \rightarrow Y \rightarrow XH \rightarrow q\bar{q}b\bar{b}$) in pb in the two-dimensional space of $m_Y$ versus $m_X$, obtained from a simultaneous fit of both merged and resolved two-prong signal regions with all statistical and systematic uncertainties.

Breakdown of systematic uncertainties in the measured fiducial cross-section. The impact of the top-quark mass on the cross-section is included in the table and not counted in the total uncertainty entry in the paper.

Data bootstrap post unfolding for the fiducial cross-section. The replicas are obtained by reweighting each observed data event by a random integer generated according to Poisson statistics, using the BootstrapGenerator software package (https://gitlab.cern.ch/atlas-physics/sm/StandardModelTools_BootstrapGenerator/BootstrapGenerator), which implements a technique described in ATL-PHYS-PUB-2021-011 (https://cds.cern.ch/record/2759945). The ATLAS event number and run number of each event are used as seed to uniquely but reproducibly initialise the random number generator for each event. All the provided numbers originate from pseudo-data, including the 0th entry, and are in units of [fb].

Data bootstrap post unfolding for the fiducial cross-section. The replicas are obtained by reweighting each observed data event by a random integer generated according to Poisson statistics, using the BootstrapGenerator software package (https://gitlab.cern.ch/atlas-physics/sm/StandardModelTools_BootstrapGenerator/BootstrapGenerator), which implements a technique described in ATL-PHYS-PUB-2021-011 (https://cds.cern.ch/record/2759945). The ATLAS event number and run number of each event are used as seed to uniquely but reproducibly initialise the random number generator for each event. All the provided numbers originate from pseudo-data, including the 0th entry, and are in units of [fb].

Breakdown of systematic uncertainties in the measured inclusive cross-section. The impact of the top-quark mass on the cross-section is included in the table and not counted in the total uncertainty entry in the paper.

Breakdown of systematic uncertainties in the measured inclusive cross-section. The impact of the top-quark mass on the cross-section is included in the table and not counted in the total uncertainty entry in the paper.

Data bootstrap post unfolding for the inclusive cross-section. The replicas are obtained by reweighting each observed data event by a random integer generated according to Poisson statistics, using the BootstrapGenerator software package (https://gitlab.cern.ch/atlas-physics/sm/StandardModelTools_BootstrapGenerator/BootstrapGenerator), which implements a technique described in ATL-PHYS-PUB-2021-011 (https://cds.cern.ch/record/2759945). The ATLAS event number and run number of each event are used as seed to uniquely but reproducibly initialise the random number generator for each event. All the provided numbers originate from pseudo-data, including the 0th entry, and are in units of [fb].

Data bootstrap post unfolding for the inclusive cross-section. The replicas are obtained by reweighting each observed data event by a random integer generated according to Poisson statistics, using the BootstrapGenerator software package (https://gitlab.cern.ch/atlas-physics/sm/StandardModelTools_BootstrapGenerator/BootstrapGenerator), which implements a technique described in ATL-PHYS-PUB-2021-011 (https://cds.cern.ch/record/2759945). The ATLAS event number and run number of each event are used as seed to uniquely but reproducibly initialise the random number generator for each event. All the provided numbers originate from pseudo-data, including the 0th entry, and are in units of [fb].

Differential absolute cross-section for $\textrm{p}_{\textrm{T}}^{l}$. The impact of the top-quark mass on the cross-section is included in the table and not counted in the total uncertainty entry in the paper. The last bin of the distribution contains the overflow.

Differential absolute cross-section for $\textrm{p}_{\textrm{T}}^{l}$. The impact of the top-quark mass on the cross-section is included in the table and not counted in the total uncertainty entry in the paper. The last bin of the distribution contains the overflow.

Covariance Matrix (statistical uncertainties only) for the differential absolute cross-section for $\textrm{p}_{\textrm{T}}^{l}$. The last bin of the distribution contains the overflow.

Covariance Matrix (statistical uncertainties only) for the differential absolute cross-section for $\textrm{p}_{\textrm{T}}^{l}$. The last bin of the distribution contains the overflow.

Data bootstrap post unfolding for the differential absolute cross-section for $\textrm{p}_{\textrm{T}}^{l}$. The replicas are obtained by reweighting each observed data event by a random integer generated according to Poisson statistics, using the BootstrapGenerator software package (https://gitlab.cern.ch/atlas-physics/sm/StandardModelTools_BootstrapGenerator/BootstrapGenerator), which implements a technique described in ATL-PHYS-PUB-2021-011 (https://cds.cern.ch/record/2759945). The ATLAS event number and run number of each event are used as seed to uniquely but reproducibly initialise the random number generator for each event. All the provided numbers originate from pseudo-data, including the 0th entry, and are in units of [fb/GeV]. The last bin of the distribution contains the overflow.

Data bootstrap post unfolding for the differential absolute cross-section for $\textrm{p}_{\textrm{T}}^{l}$. The replicas are obtained by reweighting each observed data event by a random integer generated according to Poisson statistics, using the BootstrapGenerator software package (https://gitlab.cern.ch/atlas-physics/sm/StandardModelTools_BootstrapGenerator/BootstrapGenerator), which implements a technique described in ATL-PHYS-PUB-2021-011 (https://cds.cern.ch/record/2759945). The ATLAS event number and run number of each event are used as seed to uniquely but reproducibly initialise the random number generator for each event. All the provided numbers originate from pseudo-data, including the 0th entry, and are in units of [fb/GeV]. The last bin of the distribution contains the overflow.

Differential absolute cross-section for $|\eta^{l}|$. The impact of the top-quark mass on the cross-section is included in the table and not counted in the total uncertainty entry in the paper.

Differential absolute cross-section for $|\eta^{l}|$. The impact of the top-quark mass on the cross-section is included in the table and not counted in the total uncertainty entry in the paper.

Covariance Matrix (statistical uncertainties only) for the differential absolute cross-section for $|\eta^{l}|$.

Covariance Matrix (statistical uncertainties only) for the differential absolute cross-section for $|\eta^{l}|$.

Data bootstrap post unfolding for the differential absolute cross-section for $|\eta^{l}|$. The replicas are obtained by reweighting each observed data event by a random integer generated according to Poisson statistics, using the BootstrapGenerator software package (https://gitlab.cern.ch/atlas-physics/sm/StandardModelTools_BootstrapGenerator/BootstrapGenerator), which implements a technique described in ATL-PHYS-PUB-2021-011 (https://cds.cern.ch/record/2759945). The ATLAS event number and run number of each event are used as seed to uniquely but reproducibly initialise the random number generator for each event. All the provided numbers originate from pseudo-data, including the 0th entry, and are in units of [fb].

Data bootstrap post unfolding for the differential absolute cross-section for $|\eta^{l}|$. The replicas are obtained by reweighting each observed data event by a random integer generated according to Poisson statistics, using the BootstrapGenerator software package (https://gitlab.cern.ch/atlas-physics/sm/StandardModelTools_BootstrapGenerator/BootstrapGenerator), which implements a technique described in ATL-PHYS-PUB-2021-011 (https://cds.cern.ch/record/2759945). The ATLAS event number and run number of each event are used as seed to uniquely but reproducibly initialise the random number generator for each event. All the provided numbers originate from pseudo-data, including the 0th entry, and are in units of [fb].

Differential absolute cross-section for $\textrm{E}^{e} + \textrm{E}^{\mu}$. The impact of the top-quark mass on the cross-section is included in the table and not counted in the total uncertainty entry in the paper. The last bin of the distribution contains the overflow.

Differential absolute cross-section for $\textrm{E}^{e} + \textrm{E}^{\mu}$. The impact of the top-quark mass on the cross-section is included in the table and not counted in the total uncertainty entry in the paper. The last bin of the distribution contains the overflow.

Covariance Matrix (statistical uncertainties only) for the differential absolute cross-section for $\textrm{E}^{e} + \textrm{E}^{\mu}$. The last bin of the distribution contains the overflow.

Covariance Matrix (statistical uncertainties only) for the differential absolute cross-section for $\textrm{E}^{e} + \textrm{E}^{\mu}$. The last bin of the distribution contains the overflow.

Data bootstrap post unfolding for the differential absolute cross-section for $\textrm{E}^{e} + \textrm{E}^{\mu}$. The replicas are obtained by reweighting each observed data event by a random integer generated according to Poisson statistics, using the BootstrapGenerator software package (https://gitlab.cern.ch/atlas-physics/sm/StandardModelTools_BootstrapGenerator/BootstrapGenerator), which implements a technique described in ATL-PHYS-PUB-2021-011 (https://cds.cern.ch/record/2759945). The ATLAS event number and run number of each event are used as seed to uniquely but reproducibly initialise the random number generator for each event. All the provided numbers originate from pseudo-data, including the 0th entry, and are in units of [fb/GeV]. The last bin of the distribution contains the overflow.

Data bootstrap post unfolding for the differential absolute cross-section for $\textrm{E}^{e} + \textrm{E}^{\mu}$. The replicas are obtained by reweighting each observed data event by a random integer generated according to Poisson statistics, using the BootstrapGenerator software package (https://gitlab.cern.ch/atlas-physics/sm/StandardModelTools_BootstrapGenerator/BootstrapGenerator), which implements a technique described in ATL-PHYS-PUB-2021-011 (https://cds.cern.ch/record/2759945). The ATLAS event number and run number of each event are used as seed to uniquely but reproducibly initialise the random number generator for each event. All the provided numbers originate from pseudo-data, including the 0th entry, and are in units of [fb/GeV]. The last bin of the distribution contains the overflow.

Differential absolute cross-section for $m^{e\mu}$. The impact of the top-quark mass on the cross-section is included in the table and not counted in the total uncertainty entry in the paper. The last bin of the distribution contains the overflow.

Differential absolute cross-section for $m^{e\mu}$. The impact of the top-quark mass on the cross-section is included in the table and not counted in the total uncertainty entry in the paper. The last bin of the distribution contains the overflow.

Covariance Matrix (statistical uncertainties only) for the differential absolute cross-section for $m^{e\mu}$. The last bin of the distribution contains the overflow.

Covariance Matrix (statistical uncertainties only) for the differential absolute cross-section for $m^{e\mu}$. The last bin of the distribution contains the overflow.

Data bootstrap post unfolding for the differential absolute cross-section for $m^{e\mu}$. The replicas are obtained by reweighting each observed data event by a random integer generated according to Poisson statistics, using the BootstrapGenerator software package (https://gitlab.cern.ch/atlas-physics/sm/StandardModelTools_BootstrapGenerator/BootstrapGenerator), which implements a technique described in ATL-PHYS-PUB-2021-011 (https://cds.cern.ch/record/2759945). The ATLAS event number and run number of each event are used as seed to uniquely but reproducibly initialise the random number generator for each event. All the provided numbers originate from pseudo-data, including the 0th entry, and are in units of [fb/GeV]. The last bin of the distribution contains the overflow.

Data bootstrap post unfolding for the differential absolute cross-section for $m^{e\mu}$. The replicas are obtained by reweighting each observed data event by a random integer generated according to Poisson statistics, using the BootstrapGenerator software package (https://gitlab.cern.ch/atlas-physics/sm/StandardModelTools_BootstrapGenerator/BootstrapGenerator), which implements a technique described in ATL-PHYS-PUB-2021-011 (https://cds.cern.ch/record/2759945). The ATLAS event number and run number of each event are used as seed to uniquely but reproducibly initialise the random number generator for each event. All the provided numbers originate from pseudo-data, including the 0th entry, and are in units of [fb/GeV]. The last bin of the distribution contains the overflow.

Differential absolute cross-section for $\textrm{p}_{\textrm{T}}^{e} + \textrm{p}_{\textrm{T}}^{\mu}$. The impact of the top-quark mass on the cross-section is included in the table and not counted in the total uncertainty entry in the paper. The last bin of the distribution contains the overflow.

Differential absolute cross-section for $\textrm{p}_{\textrm{T}}^{e} + \textrm{p}_{\textrm{T}}^{\mu}$. The impact of the top-quark mass on the cross-section is included in the table and not counted in the total uncertainty entry in the paper. The last bin of the distribution contains the overflow.

Covariance Matrix (statistical uncertainties only) for the differential absolute cross-section for $\textrm{p}_{\textrm{T}}^{e} + \textrm{p}_{\textrm{T}}^{\mu}$. The last bin of the distribution contains the overflow.

Covariance Matrix (statistical uncertainties only) for the differential absolute cross-section for $\textrm{p}_{\textrm{T}}^{e} + \textrm{p}_{\textrm{T}}^{\mu}$. The last bin of the distribution contains the overflow.

Data bootstrap post unfolding for the differential absolute cross-section for $\textrm{p}_{\textrm{T}}^{e} + \textrm{p}_{\textrm{T}}^{\mu}$. The replicas are obtained by reweighting each observed data event by a random integer generated according to Poisson statistics, using the BootstrapGenerator software package (https://gitlab.cern.ch/atlas-physics/sm/StandardModelTools_BootstrapGenerator/BootstrapGenerator), which implements a technique described in ATL-PHYS-PUB-2021-011 (https://cds.cern.ch/record/2759945). The ATLAS event number and run number of each event are used as seed to uniquely but reproducibly initialise the random number generator for each event. All the provided numbers originate from pseudo-data, including the 0th entry, and are in units of [fb/GeV]. The last bin of the distribution contains the overflow.

Data bootstrap post unfolding for the differential absolute cross-section for $\textrm{p}_{\textrm{T}}^{e} + \textrm{p}_{\textrm{T}}^{\mu}$. The replicas are obtained by reweighting each observed data event by a random integer generated according to Poisson statistics, using the BootstrapGenerator software package (https://gitlab.cern.ch/atlas-physics/sm/StandardModelTools_BootstrapGenerator/BootstrapGenerator), which implements a technique described in ATL-PHYS-PUB-2021-011 (https://cds.cern.ch/record/2759945). The ATLAS event number and run number of each event are used as seed to uniquely but reproducibly initialise the random number generator for each event. All the provided numbers originate from pseudo-data, including the 0th entry, and are in units of [fb/GeV]. The last bin of the distribution contains the overflow.

Differential absolute cross-section for $\textrm{p}_{\textrm{T}}^{e\mu}$. The impact of the top-quark mass on the cross-section is included in the table and not counted in the total uncertainty entry in the paper. The last bin of the distribution contains the overflow.

Differential absolute cross-section for $\textrm{p}_{\textrm{T}}^{e\mu}$. The impact of the top-quark mass on the cross-section is included in the table and not counted in the total uncertainty entry in the paper. The last bin of the distribution contains the overflow.

Covariance Matrix (statistical uncertainties only) for the differential absolute cross-section for $\textrm{p}_{\textrm{T}}^{e\mu}$. The last bin of the distribution contains the overflow.

Covariance Matrix (statistical uncertainties only) for the differential absolute cross-section for $\textrm{p}_{\textrm{T}}^{e\mu}$. The last bin of the distribution contains the overflow.

Data bootstrap post unfolding for the differential absolute cross-section for $\textrm{p}_{\textrm{T}}^{e\mu}$. The replicas are obtained by reweighting each observed data event by a random integer generated according to Poisson statistics, using the BootstrapGenerator software package (https://gitlab.cern.ch/atlas-physics/sm/StandardModelTools_BootstrapGenerator/BootstrapGenerator), which implements a technique described in ATL-PHYS-PUB-2021-011 (https://cds.cern.ch/record/2759945). The ATLAS event number and run number of each event are used as seed to uniquely but reproducibly initialise the random number generator for each event. All the provided numbers originate from pseudo-data, including the 0th entry, and are in units of [fb/GeV]. The last bin of the distribution contains the overflow.

Data bootstrap post unfolding for the differential absolute cross-section for $\textrm{p}_{\textrm{T}}^{e\mu}$. The replicas are obtained by reweighting each observed data event by a random integer generated according to Poisson statistics, using the BootstrapGenerator software package (https://gitlab.cern.ch/atlas-physics/sm/StandardModelTools_BootstrapGenerator/BootstrapGenerator), which implements a technique described in ATL-PHYS-PUB-2021-011 (https://cds.cern.ch/record/2759945). The ATLAS event number and run number of each event are used as seed to uniquely but reproducibly initialise the random number generator for each event. All the provided numbers originate from pseudo-data, including the 0th entry, and are in units of [fb/GeV]. The last bin of the distribution contains the overflow.

Differential absolute cross-section for $|\Delta\phi^{e\mu}|$. The impact of the top-quark mass on the cross-section is included in the table and not counted in the total uncertainty entry in the paper.

Differential absolute cross-section for $|\Delta\phi^{e\mu}|$. The impact of the top-quark mass on the cross-section is included in the table and not counted in the total uncertainty entry in the paper.

Covariance Matrix (statistical uncertainties only) for the differential absolute cross-section for $|\Delta\phi^{e\mu}|$.

Covariance Matrix (statistical uncertainties only) for the differential absolute cross-section for $|\Delta\phi^{e\mu}|$.

Data bootstrap post unfolding for the differential absolute cross-section for $|\Delta\phi^{e\mu}|$. The replicas are obtained by reweighting each observed data event by a random integer generated according to Poisson statistics, using the BootstrapGenerator software package (https://gitlab.cern.ch/atlas-physics/sm/StandardModelTools_BootstrapGenerator/BootstrapGenerator), which implements a technique described in ATL-PHYS-PUB-2021-011 (https://cds.cern.ch/record/2759945). The ATLAS event number and run number of each event are used as seed to uniquely but reproducibly initialise the random number generator for each event. All the provided numbers originate from pseudo-data, including the 0th entry, and are in units of [fb/rad].

Data bootstrap post unfolding for the differential absolute cross-section for $|\Delta\phi^{e\mu}|$. The replicas are obtained by reweighting each observed data event by a random integer generated according to Poisson statistics, using the BootstrapGenerator software package (https://gitlab.cern.ch/atlas-physics/sm/StandardModelTools_BootstrapGenerator/BootstrapGenerator), which implements a technique described in ATL-PHYS-PUB-2021-011 (https://cds.cern.ch/record/2759945). The ATLAS event number and run number of each event are used as seed to uniquely but reproducibly initialise the random number generator for each event. All the provided numbers originate from pseudo-data, including the 0th entry, and are in units of [fb/rad].

Differential absolute cross-section for $|y^{e\mu}|$. The impact of the top-quark mass on the cross-section is included in the table and not counted in the total uncertainty entry in the paper.

Differential absolute cross-section for $|y^{e\mu}|$. The impact of the top-quark mass on the cross-section is included in the table and not counted in the total uncertainty entry in the paper.

Covariance Matrix (statistical uncertainties only) for the differential absolute cross-section for $|y^{e\mu}|$.

Covariance Matrix (statistical uncertainties only) for the differential absolute cross-section for $|y^{e\mu}|$.

Data bootstrap post unfolding for the differential absolute cross-section for $|y^{e\mu}|$. The replicas are obtained by reweighting each observed data event by a random integer generated according to Poisson statistics, using the BootstrapGenerator software package (https://gitlab.cern.ch/atlas-physics/sm/StandardModelTools_BootstrapGenerator/BootstrapGenerator), which implements a technique described in ATL-PHYS-PUB-2021-011 (https://cds.cern.ch/record/2759945). The ATLAS event number and run number of each event are used as seed to uniquely but reproducibly initialise the random number generator for each event. All the provided numbers originate from pseudo-data, including the 0th entry, and are in units of [fb].

Data bootstrap post unfolding for the differential absolute cross-section for $|y^{e\mu}|$. The replicas are obtained by reweighting each observed data event by a random integer generated according to Poisson statistics, using the BootstrapGenerator software package (https://gitlab.cern.ch/atlas-physics/sm/StandardModelTools_BootstrapGenerator/BootstrapGenerator), which implements a technique described in ATL-PHYS-PUB-2021-011 (https://cds.cern.ch/record/2759945). The ATLAS event number and run number of each event are used as seed to uniquely but reproducibly initialise the random number generator for each event. All the provided numbers originate from pseudo-data, including the 0th entry, and are in units of [fb].

Double differential absolute cross-section for $|y^{e\mu}|$ x $m^{e\mu}$. The impact of the top-quark mass on the cross-section is included in the table and not counted in the total uncertainty entry in the paper. The highest invariant mass bin of the two distributions contains the overflows.

Double differential absolute cross-section for $|y^{e\mu}|$ x $m^{e\mu}$. The impact of the top-quark mass on the cross-section is included in the table and not counted in the total uncertainty entry in the paper. The highest invariant mass bin of the two distributions contains the overflows.

Covariance Matrix (statistical uncertainties only) for the double differential absolute cross-section for $|y^{e\mu}|$ x $m^{e\mu}$. The highest invariant mass bin of the two distributions contains the overflows.

Covariance Matrix (statistical uncertainties only) for the double differential absolute cross-section for $|y^{e\mu}|$ x $m^{e\mu}$. The highest invariant mass bin of the two distributions contains the overflows.

Data bootstrap post unfolding for the double differential absolute cross-section for $|y^{e\mu}|$ x $m^{e\mu}$. The replicas are obtained by reweighting each observed data event by a random integer generated according to Poisson statistics, using the BootstrapGenerator software package (https://gitlab.cern.ch/atlas-physics/sm/StandardModelTools_BootstrapGenerator/BootstrapGenerator), which implements a technique described in ATL-PHYS-PUB-2021-011 (https://cds.cern.ch/record/2759945). The ATLAS event number and run number of each event are used as seed to uniquely but reproducibly initialise the random number generator for each event. All the provided numbers originate from pseudo-data, including the 0th entry, and are in units of [fb/GeV]. The highest invariant mass bin of the two distributions contains the overflows.

Data bootstrap post unfolding for the double differential absolute cross-section for $|y^{e\mu}|$ x $m^{e\mu}$. The replicas are obtained by reweighting each observed data event by a random integer generated according to Poisson statistics, using the BootstrapGenerator software package (https://gitlab.cern.ch/atlas-physics/sm/StandardModelTools_BootstrapGenerator/BootstrapGenerator), which implements a technique described in ATL-PHYS-PUB-2021-011 (https://cds.cern.ch/record/2759945). The ATLAS event number and run number of each event are used as seed to uniquely but reproducibly initialise the random number generator for each event. All the provided numbers originate from pseudo-data, including the 0th entry, and are in units of [fb/GeV]. The highest invariant mass bin of the two distributions contains the overflows.

Double differential absolute cross-section for $|\Delta\phi^{e\mu}|$ x $m^{e\mu}$. The impact of the top-quark mass on the cross-section is included in the table and not counted in the total uncertainty entry in the paper. The highest invariant mass bin of the two distributions contains the overflows.

Double differential absolute cross-section for $|\Delta\phi^{e\mu}|$ x $m^{e\mu}$. The impact of the top-quark mass on the cross-section is included in the table and not counted in the total uncertainty entry in the paper. The highest invariant mass bin of the two distributions contains the overflows.

Covariance Matrix (statistical uncertainties only) for the double differential absolute cross-section for $|\Delta\phi^{e\mu}|$ x $m^{e\mu}$. The highest invariant mass bin of the two distributions contains the overflows.

Covariance Matrix (statistical uncertainties only) for the double differential absolute cross-section for $|\Delta\phi^{e\mu}|$ x $m^{e\mu}$. The highest invariant mass bin of the two distributions contains the overflows.

Data bootstrap post unfolding for the double differential absolute cross-section for $|\Delta\phi^{e\mu}|$ x $m^{e\mu}$. The replicas are obtained by reweighting each observed data event by a random integer generated according to Poisson statistics, using the BootstrapGenerator software package (https://gitlab.cern.ch/atlas-physics/sm/StandardModelTools_BootstrapGenerator/BootstrapGenerator), which implements a technique described in ATL-PHYS-PUB-2021-011 (https://cds.cern.ch/record/2759945). The ATLAS event number and run number of each event are used as seed to uniquely but reproducibly initialise the random number generator for each event. All the provided numbers originate from pseudo-data, including the 0th entry, and are in units of [fb/GeV/rad]. The highest invariant mass bin of the two distributions contains the overflows.

Data bootstrap post unfolding for the double differential absolute cross-section for $|\Delta\phi^{e\mu}|$ x $m^{e\mu}$. The replicas are obtained by reweighting each observed data event by a random integer generated according to Poisson statistics, using the BootstrapGenerator software package (https://gitlab.cern.ch/atlas-physics/sm/StandardModelTools_BootstrapGenerator/BootstrapGenerator), which implements a technique described in ATL-PHYS-PUB-2021-011 (https://cds.cern.ch/record/2759945). The ATLAS event number and run number of each event are used as seed to uniquely but reproducibly initialise the random number generator for each event. All the provided numbers originate from pseudo-data, including the 0th entry, and are in units of [fb/GeV/rad]. The highest invariant mass bin of the two distributions contains the overflows.

Double differential absolute cross-section for $|\Delta\phi^{e\mu}|$ x $\textrm{p}_{\textrm{T}}^{e\mu}$. The impact of the top-quark mass on the cross-section is included in the table and not counted in the total uncertainty entry in the paper. The highest $\textrm{p}_{\textrm{T}}^{e\mu}$ bin contains the overflows.

Double differential absolute cross-section for $|\Delta\phi^{e\mu}|$ x $\textrm{p}_{\textrm{T}}^{e\mu}$. The impact of the top-quark mass on the cross-section is included in the table and not counted in the total uncertainty entry in the paper. The highest $\textrm{p}_{\textrm{T}}^{e\mu}$ bin contains the overflows.

Covariance Matrix (statistical uncertainties only) for the double differential absolute cross-section for $|\Delta\phi^{e\mu}|$ x $\textrm{p}_{\textrm{T}}^{e\mu}$. The highest $\textrm{p}_{\textrm{T}}^{e\mu}$ bin contains the overflows.

Covariance Matrix (statistical uncertainties only) for the double differential absolute cross-section for $|\Delta\phi^{e\mu}|$ x $\textrm{p}_{\textrm{T}}^{e\mu}$. The highest $\textrm{p}_{\textrm{T}}^{e\mu}$ bin contains the overflows.

Data bootstrap post unfolding for the double differential absolute cross-section for $|\Delta\phi^{e\mu}|$ x $\textrm{p}_{\textrm{T}}^{e\mu}$. The replicas are obtained by reweighting each observed data event by a random integer generated according to Poisson statistics, using the BootstrapGenerator software package (https://gitlab.cern.ch/atlas-physics/sm/StandardModelTools_BootstrapGenerator/BootstrapGenerator), which implements a technique described in ATL-PHYS-PUB-2021-011 (https://cds.cern.ch/record/2759945). The ATLAS event number and run number of each event are used as seed to uniquely but reproducibly initialise the random number generator for each event. All the provided numbers originate from pseudo-data, including the 0th entry, and are in units of [fb/GeV/rad]. The highest $\textrm{p}_{\textrm{T}}^{e\mu}$ bin contains the overflows.

Data bootstrap post unfolding for the double differential absolute cross-section for $|\Delta\phi^{e\mu}|$ x $\textrm{p}_{\textrm{T}}^{e\mu}$. The replicas are obtained by reweighting each observed data event by a random integer generated according to Poisson statistics, using the BootstrapGenerator software package (https://gitlab.cern.ch/atlas-physics/sm/StandardModelTools_BootstrapGenerator/BootstrapGenerator), which implements a technique described in ATL-PHYS-PUB-2021-011 (https://cds.cern.ch/record/2759945). The ATLAS event number and run number of each event are used as seed to uniquely but reproducibly initialise the random number generator for each event. All the provided numbers originate from pseudo-data, including the 0th entry, and are in units of [fb/GeV/rad]. The highest $\textrm{p}_{\textrm{T}}^{e\mu}$ bin contains the overflows.

Double differential absolute cross-section for $|\Delta\phi^{e\mu}|$ x $\textrm{E}^{e} + \textrm{E}^{\mu}$. The impact of the top-quark mass on the cross-section is included in the table and not counted in the total uncertainty entry in the paper. The highest $\textrm{E}^{e} + \textrm{E}^{\mu}$ bin contains the overflows.

Double differential absolute cross-section for $|\Delta\phi^{e\mu}|$ x $\textrm{E}^{e} + \textrm{E}^{\mu}$. The impact of the top-quark mass on the cross-section is included in the table and not counted in the total uncertainty entry in the paper. The highest $\textrm{E}^{e} + \textrm{E}^{\mu}$ bin contains the overflows.

Covariance Matrix (statistical uncertainties only) for the double differential absolute cross-section for $|\Delta\phi^{e\mu}|$ x $\textrm{E}^{e} + \textrm{E}^{\mu}$. The highest $\textrm{E}^{e} + \textrm{E}^{\mu}$ bin contains the overflows.

Covariance Matrix (statistical uncertainties only) for the double differential absolute cross-section for $|\Delta\phi^{e\mu}|$ x $\textrm{E}^{e} + \textrm{E}^{\mu}$. The highest $\textrm{E}^{e} + \textrm{E}^{\mu}$ bin contains the overflows.

Data bootstrap post unfolding for the double differential absolute cross-section for $|\Delta\phi^{e\mu}|$ x $\textrm{E}^{e} + \textrm{E}^{\mu}$. The replicas are obtained by reweighting each observed data event by a random integer generated according to Poisson statistics, using the BootstrapGenerator software package (https://gitlab.cern.ch/atlas-physics/sm/StandardModelTools_BootstrapGenerator/BootstrapGenerator), which implements a technique described in ATL-PHYS-PUB-2021-011 (https://cds.cern.ch/record/2759945). The ATLAS event number and run number of each event are used as seed to uniquely but reproducibly initialise the random number generator for each event. All the provided numbers originate from pseudo-data, including the 0th entry, and are in units of [fb/GeV/rad]. The highest $\textrm{E}^{e} + \textrm{E}^{\mu}$ bin contains the overflows.

Data bootstrap post unfolding for the double differential absolute cross-section for $|\Delta\phi^{e\mu}|$ x $\textrm{E}^{e} + \textrm{E}^{\mu}$. The replicas are obtained by reweighting each observed data event by a random integer generated according to Poisson statistics, using the BootstrapGenerator software package (https://gitlab.cern.ch/atlas-physics/sm/StandardModelTools_BootstrapGenerator/BootstrapGenerator), which implements a technique described in ATL-PHYS-PUB-2021-011 (https://cds.cern.ch/record/2759945). The ATLAS event number and run number of each event are used as seed to uniquely but reproducibly initialise the random number generator for each event. All the provided numbers originate from pseudo-data, including the 0th entry, and are in units of [fb/GeV/rad]. The highest $\textrm{E}^{e} + \textrm{E}^{\mu}$ bin contains the overflows.

Differential normalised cross-section for $\textrm{p}_{\textrm{T}}^{l}$. The impact of the top-quark mass on the cross-section is included in the table and not counted in the total uncertainty entry in the paper. The last bin of the distribution contains the overflow.

Differential normalised cross-section for $\textrm{p}_{\textrm{T}}^{l}$. The impact of the top-quark mass on the cross-section is included in the table and not counted in the total uncertainty entry in the paper. The last bin of the distribution contains the overflow.

Covariance Matrix (statistical uncertainties only) for the differential normalised cross-section for $\textrm{p}_{\textrm{T}}^{l}$. The last bin of the distribution contains the overflow.

Covariance Matrix (statistical uncertainties only) for the differential normalised cross-section for $\textrm{p}_{\textrm{T}}^{l}$. The last bin of the distribution contains the overflow.

Data bootstrap post unfolding for the differential normalised cross-section for $\textrm{p}_{\textrm{T}}^{l}$. The replicas are obtained by reweighting each observed data event by a random integer generated according to Poisson statistics, using the BootstrapGenerator software package (https://gitlab.cern.ch/atlas-physics/sm/StandardModelTools_BootstrapGenerator/BootstrapGenerator), which implements a technique described in ATL-PHYS-PUB-2021-011 (https://cds.cern.ch/record/2759945). The ATLAS event number and run number of each event are used as seed to uniquely but reproducibly initialise the random number generator for each event. All the provided numbers originate from pseudo-data, including the 0th entry, and are in units of [1/GeV]. The last bin of the distribution contains the overflow.

Data bootstrap post unfolding for the differential normalised cross-section for $\textrm{p}_{\textrm{T}}^{l}$. The replicas are obtained by reweighting each observed data event by a random integer generated according to Poisson statistics, using the BootstrapGenerator software package (https://gitlab.cern.ch/atlas-physics/sm/StandardModelTools_BootstrapGenerator/BootstrapGenerator), which implements a technique described in ATL-PHYS-PUB-2021-011 (https://cds.cern.ch/record/2759945). The ATLAS event number and run number of each event are used as seed to uniquely but reproducibly initialise the random number generator for each event. All the provided numbers originate from pseudo-data, including the 0th entry, and are in units of [1/GeV]. The last bin of the distribution contains the overflow.

Differential normalised cross-section for $|\eta^{l}|$. The impact of the top-quark mass on the cross-section is included in the table and not counted in the total uncertainty entry in the paper.

Differential normalised cross-section for $|\eta^{l}|$. The impact of the top-quark mass on the cross-section is included in the table and not counted in the total uncertainty entry in the paper.

Covariance Matrix (statistical uncertainties only) for the differential normalised cross-section for $|\eta^{l}|$.

Covariance Matrix (statistical uncertainties only) for the differential normalised cross-section for $|\eta^{l}|$.

Data bootstrap post unfolding for the differential normalised cross-section for $|\eta^{l}|$. The replicas are obtained by reweighting each observed data event by a random integer generated according to Poisson statistics, using the BootstrapGenerator software package (https://gitlab.cern.ch/atlas-physics/sm/StandardModelTools_BootstrapGenerator/BootstrapGenerator), which implements a technique described in ATL-PHYS-PUB-2021-011 (https://cds.cern.ch/record/2759945). The ATLAS event number and run number of each event are used as seed to uniquely but reproducibly initialise the random number generator for each event. All the provided numbers originate from pseudo-data, including the 0th entry.

Data bootstrap post unfolding for the differential normalised cross-section for $|\eta^{l}|$. The replicas are obtained by reweighting each observed data event by a random integer generated according to Poisson statistics, using the BootstrapGenerator software package (https://gitlab.cern.ch/atlas-physics/sm/StandardModelTools_BootstrapGenerator/BootstrapGenerator), which implements a technique described in ATL-PHYS-PUB-2021-011 (https://cds.cern.ch/record/2759945). The ATLAS event number and run number of each event are used as seed to uniquely but reproducibly initialise the random number generator for each event. All the provided numbers originate from pseudo-data, including the 0th entry.

Differential normalised cross-section for $\textrm{E}^{e} + \textrm{E}^{\mu}$. The impact of the top-quark mass on the cross-section is included in the table and not counted in the total uncertainty entry in the paper. The last bin of the distribution contains the overflow.

Differential normalised cross-section for $\textrm{E}^{e} + \textrm{E}^{\mu}$. The impact of the top-quark mass on the cross-section is included in the table and not counted in the total uncertainty entry in the paper. The last bin of the distribution contains the overflow.

Covariance Matrix (statistical uncertainties only) for the differential normalised cross-section for $\textrm{E}^{e} + \textrm{E}^{\mu}$. The last bin of the distribution contains the overflow.

Covariance Matrix (statistical uncertainties only) for the differential normalised cross-section for $\textrm{E}^{e} + \textrm{E}^{\mu}$. The last bin of the distribution contains the overflow.

Data bootstrap post unfolding for the differential normalised cross-section for $\textrm{E}^{e} + \textrm{E}^{\mu}$. The replicas are obtained by reweighting each observed data event by a random integer generated according to Poisson statistics, using the BootstrapGenerator software package (https://gitlab.cern.ch/atlas-physics/sm/StandardModelTools_BootstrapGenerator/BootstrapGenerator), which implements a technique described in ATL-PHYS-PUB-2021-011 (https://cds.cern.ch/record/2759945). The ATLAS event number and run number of each event are used as seed to uniquely but reproducibly initialise the random number generator for each event. All the provided numbers originate from pseudo-data, including the 0th entry. The last bin of the distribution contains the overflow.

Data bootstrap post unfolding for the differential normalised cross-section for $\textrm{E}^{e} + \textrm{E}^{\mu}$. The replicas are obtained by reweighting each observed data event by a random integer generated according to Poisson statistics, using the BootstrapGenerator software package (https://gitlab.cern.ch/atlas-physics/sm/StandardModelTools_BootstrapGenerator/BootstrapGenerator), which implements a technique described in ATL-PHYS-PUB-2021-011 (https://cds.cern.ch/record/2759945). The ATLAS event number and run number of each event are used as seed to uniquely but reproducibly initialise the random number generator for each event. All the provided numbers originate from pseudo-data, including the 0th entry. The last bin of the distribution contains the overflow.

Differential normalised cross-section for $m^{e\mu}$. The impact of the top-quark mass on the cross-section is included in the table and not counted in the total uncertainty entry in the paper. The last bin of the distribution contains the overflow.

Differential normalised cross-section for $m^{e\mu}$. The impact of the top-quark mass on the cross-section is included in the table and not counted in the total uncertainty entry in the paper. The last bin of the distribution contains the overflow.

Covariance Matrix (statistical uncertainties only) for the differential normalised cross-section for $m^{e\mu}$. The last bin of the distribution contains the overflow.

Covariance Matrix (statistical uncertainties only) for the differential normalised cross-section for $m^{e\mu}$. The last bin of the distribution contains the overflow.

Data bootstrap post unfolding for the differential normalised cross-section for $m^{e\mu}$. The replicas are obtained by reweighting each observed data event by a random integer generated according to Poisson statistics, using the BootstrapGenerator software package (https://gitlab.cern.ch/atlas-physics/sm/StandardModelTools_BootstrapGenerator/BootstrapGenerator), which implements a technique described in ATL-PHYS-PUB-2021-011 (https://cds.cern.ch/record/2759945). The ATLAS event number and run number of each event are used as seed to uniquely but reproducibly initialise the random number generator for each event. All the provided numbers originate from pseudo-data, including the 0th entry. The last bin of the distribution contains the overflow.

Data bootstrap post unfolding for the differential normalised cross-section for $m^{e\mu}$. The replicas are obtained by reweighting each observed data event by a random integer generated according to Poisson statistics, using the BootstrapGenerator software package (https://gitlab.cern.ch/atlas-physics/sm/StandardModelTools_BootstrapGenerator/BootstrapGenerator), which implements a technique described in ATL-PHYS-PUB-2021-011 (https://cds.cern.ch/record/2759945). The ATLAS event number and run number of each event are used as seed to uniquely but reproducibly initialise the random number generator for each event. All the provided numbers originate from pseudo-data, including the 0th entry. The last bin of the distribution contains the overflow.

Differential normalised cross-section for $\textrm{p}_{\textrm{T}}^{e} + \textrm{p}_{\textrm{T}}^{\mu}$. The impact of the top-quark mass on the cross-section is included in the table and not counted in the total uncertainty entry in the paper. The last bin of the distribution contains the overflow.

Differential normalised cross-section for $\textrm{p}_{\textrm{T}}^{e} + \textrm{p}_{\textrm{T}}^{\mu}$. The impact of the top-quark mass on the cross-section is included in the table and not counted in the total uncertainty entry in the paper. The last bin of the distribution contains the overflow.

Covariance Matrix (statistical uncertainties only) for the differential normalised cross-section for $\textrm{p}_{\textrm{T}}^{e} + \textrm{p}_{\textrm{T}}^{\mu}$. The last bin of the distribution contains the overflow.

Covariance Matrix (statistical uncertainties only) for the differential normalised cross-section for $\textrm{p}_{\textrm{T}}^{e} + \textrm{p}_{\textrm{T}}^{\mu}$. The last bin of the distribution contains the overflow.

Data bootstrap post unfolding for the differential normalised cross-section for $\textrm{p}_{\textrm{T}}^{e} + \textrm{p}_{\textrm{T}}^{\mu}$. The replicas are obtained by reweighting each observed data event by a random integer generated according to Poisson statistics, using the BootstrapGenerator software package (https://gitlab.cern.ch/atlas-physics/sm/StandardModelTools_BootstrapGenerator/BootstrapGenerator), which implements a technique described in ATL-PHYS-PUB-2021-011 (https://cds.cern.ch/record/2759945). The ATLAS event number and run number of each event are used as seed to uniquely but reproducibly initialise the random number generator for each event. All the provided numbers originate from pseudo-data, including the 0th entry, and are in units of [1/GeV]. The last bin of the distribution contains the overflow.

Data bootstrap post unfolding for the differential normalised cross-section for $\textrm{p}_{\textrm{T}}^{e} + \textrm{p}_{\textrm{T}}^{\mu}$. The replicas are obtained by reweighting each observed data event by a random integer generated according to Poisson statistics, using the BootstrapGenerator software package (https://gitlab.cern.ch/atlas-physics/sm/StandardModelTools_BootstrapGenerator/BootstrapGenerator), which implements a technique described in ATL-PHYS-PUB-2021-011 (https://cds.cern.ch/record/2759945). The ATLAS event number and run number of each event are used as seed to uniquely but reproducibly initialise the random number generator for each event. All the provided numbers originate from pseudo-data, including the 0th entry, and are in units of [1/GeV]. The last bin of the distribution contains the overflow.

Differential normalised cross-section for $\textrm{p}_{\textrm{T}}^{e\mu}$. The impact of the top-quark mass on the cross-section is included in the table and not counted in the total uncertainty entry in the paper. The last bin of the distribution contains the overflow.

Differential normalised cross-section for $\textrm{p}_{\textrm{T}}^{e\mu}$. The impact of the top-quark mass on the cross-section is included in the table and not counted in the total uncertainty entry in the paper. The last bin of the distribution contains the overflow.

Covariance Matrix (statistical uncertainties only) for the differential normalised cross-section for $\textrm{p}_{\textrm{T}}^{e\mu}$. The last bin of the distribution contains the overflow.

Covariance Matrix (statistical uncertainties only) for the differential normalised cross-section for $\textrm{p}_{\textrm{T}}^{e\mu}$. The last bin of the distribution contains the overflow.

Data bootstrap post unfolding for the differential normalised cross-section for $\textrm{p}_{\textrm{T}}^{e\mu}$. The replicas are obtained by reweighting each observed data event by a random integer generated according to Poisson statistics, using the BootstrapGenerator software package (https://gitlab.cern.ch/atlas-physics/sm/StandardModelTools_BootstrapGenerator/BootstrapGenerator), which implements a technique described in ATL-PHYS-PUB-2021-011 (https://cds.cern.ch/record/2759945). The ATLAS event number and run number of each event are used as seed to uniquely but reproducibly initialise the random number generator for each event. All the provided numbers originate from pseudo-data, including the 0th entry, and are in units of [1/GeV]. The last bin of the distribution contains the overflow.

Data bootstrap post unfolding for the differential normalised cross-section for $\textrm{p}_{\textrm{T}}^{e\mu}$. The replicas are obtained by reweighting each observed data event by a random integer generated according to Poisson statistics, using the BootstrapGenerator software package (https://gitlab.cern.ch/atlas-physics/sm/StandardModelTools_BootstrapGenerator/BootstrapGenerator), which implements a technique described in ATL-PHYS-PUB-2021-011 (https://cds.cern.ch/record/2759945). The ATLAS event number and run number of each event are used as seed to uniquely but reproducibly initialise the random number generator for each event. All the provided numbers originate from pseudo-data, including the 0th entry, and are in units of [1/GeV]. The last bin of the distribution contains the overflow.

Differential normalised cross-section for $|\Delta\phi^{e\mu}|$. The impact of the top-quark mass on the cross-section is included in the table and not counted in the total uncertainty entry in the paper.

Differential normalised cross-section for $|\Delta\phi^{e\mu}|$. The impact of the top-quark mass on the cross-section is included in the table and not counted in the total uncertainty entry in the paper.

Covariance Matrix (statistical uncertainties only) for the differential normalised cross-section for $|\Delta\phi^{e\mu}|$.

Covariance Matrix (statistical uncertainties only) for the differential normalised cross-section for $|\Delta\phi^{e\mu}|$.

Data bootstrap post unfolding for the differential normalised cross-section for $|\Delta\phi^{e\mu}|$. The replicas are obtained by reweighting each observed data event by a random integer generated according to Poisson statistics, using the BootstrapGenerator software package (https://gitlab.cern.ch/atlas-physics/sm/StandardModelTools_BootstrapGenerator/BootstrapGenerator), which implements a technique described in ATL-PHYS-PUB-2021-011 (https://cds.cern.ch/record/2759945). The ATLAS event number and run number of each event are used as seed to uniquely but reproducibly initialise the random number generator for each event. All the provided numbers originate from pseudo-data, including the 0th entry, and are in units of [1/rad].

Data bootstrap post unfolding for the differential normalised cross-section for $|\Delta\phi^{e\mu}|$. The replicas are obtained by reweighting each observed data event by a random integer generated according to Poisson statistics, using the BootstrapGenerator software package (https://gitlab.cern.ch/atlas-physics/sm/StandardModelTools_BootstrapGenerator/BootstrapGenerator), which implements a technique described in ATL-PHYS-PUB-2021-011 (https://cds.cern.ch/record/2759945). The ATLAS event number and run number of each event are used as seed to uniquely but reproducibly initialise the random number generator for each event. All the provided numbers originate from pseudo-data, including the 0th entry, and are in units of [1/rad].

Differential normalised cross-section for $|y^{e\mu}|$. The impact of the top-quark mass on the cross-section is included in the table and not counted in the total uncertainty entry in the paper.

Differential normalised cross-section for $|y^{e\mu}|$. The impact of the top-quark mass on the cross-section is included in the table and not counted in the total uncertainty entry in the paper.

Covariance Matrix (statistical uncertainties only) for the differential normalised cross-section for $|y^{e\mu}|$.

Covariance Matrix (statistical uncertainties only) for the differential normalised cross-section for $|y^{e\mu}|$.

Data bootstrap post unfolding for the differential normalised cross-section for $|y^{e\mu}|$. The replicas are obtained by reweighting each observed data event by a random integer generated according to Poisson statistics, using the BootstrapGenerator software package (https://gitlab.cern.ch/atlas-physics/sm/StandardModelTools_BootstrapGenerator/BootstrapGenerator), which implements a technique described in ATL-PHYS-PUB-2021-011 (https://cds.cern.ch/record/2759945). The ATLAS event number and run number of each event are used as seed to uniquely but reproducibly initialise the random number generator for each event. All the provided numbers originate from pseudo-data, including the 0th entry.

Data bootstrap post unfolding for the differential normalised cross-section for $|y^{e\mu}|$. The replicas are obtained by reweighting each observed data event by a random integer generated according to Poisson statistics, using the BootstrapGenerator software package (https://gitlab.cern.ch/atlas-physics/sm/StandardModelTools_BootstrapGenerator/BootstrapGenerator), which implements a technique described in ATL-PHYS-PUB-2021-011 (https://cds.cern.ch/record/2759945). The ATLAS event number and run number of each event are used as seed to uniquely but reproducibly initialise the random number generator for each event. All the provided numbers originate from pseudo-data, including the 0th entry.

Double differential normalised cross-section for $|y^{e\mu}|$ x $m^{e\mu}$. The impact of the top-quark mass on the cross-section is included in the table and not counted in the total uncertainty entry in the paper. The highest invariant mass bin of the two distributions contains the overflows.

Double differential normalised cross-section for $|y^{e\mu}|$ x $m^{e\mu}$. The impact of the top-quark mass on the cross-section is included in the table and not counted in the total uncertainty entry in the paper. The highest invariant mass bin of the two distributions contains the overflows.

Covariance Matrix (statistical uncertainties only) for the double differential normalised cross-section for $|y^{e\mu}|$ x $m^{e\mu}$. The highest invariant mass bin of the two distributions contains the overflows.

Covariance Matrix (statistical uncertainties only) for the double differential normalised cross-section for $|y^{e\mu}|$ x $m^{e\mu}$. The highest invariant mass bin of the two distributions contains the overflows.

Data bootstrap post unfolding for the double differential normalised cross-section for $|y^{e\mu}|$ x $m^{e\mu}$. The replicas are obtained by reweighting each observed data event by a random integer generated according to Poisson statistics, using the BootstrapGenerator software package (https://gitlab.cern.ch/atlas-physics/sm/StandardModelTools_BootstrapGenerator/BootstrapGenerator), which implements a technique described in ATL-PHYS-PUB-2021-011 (https://cds.cern.ch/record/2759945). The ATLAS event number and run number of each event are used as seed to uniquely but reproducibly initialise the random number generator for each event. All the provided numbers originate from pseudo-data, including the 0th entry, and are in units of [1/GeV]. The highest invariant mass bin of the two distributions contains the overflows.

Data bootstrap post unfolding for the double differential normalised cross-section for $|y^{e\mu}|$ x $m^{e\mu}$. The replicas are obtained by reweighting each observed data event by a random integer generated according to Poisson statistics, using the BootstrapGenerator software package (https://gitlab.cern.ch/atlas-physics/sm/StandardModelTools_BootstrapGenerator/BootstrapGenerator), which implements a technique described in ATL-PHYS-PUB-2021-011 (https://cds.cern.ch/record/2759945). The ATLAS event number and run number of each event are used as seed to uniquely but reproducibly initialise the random number generator for each event. All the provided numbers originate from pseudo-data, including the 0th entry, and are in units of [1/GeV]. The highest invariant mass bin of the two distributions contains the overflows.

Double differential normalised cross-section for $|\Delta\phi^{e\mu}|$ x $m^{e\mu}$. The impact of the top-quark mass on the cross-section is included in the table and not counted in the total uncertainty entry in the paper. The highest invariant mass bin of the two distributions contains the overflows.

Double differential normalised cross-section for $|\Delta\phi^{e\mu}|$ x $m^{e\mu}$. The impact of the top-quark mass on the cross-section is included in the table and not counted in the total uncertainty entry in the paper. The highest invariant mass bin of the two distributions contains the overflows.

Covariance Matrix (statistical uncertainties only) for the double differential normalised cross-section for $|\Delta\phi^{e\mu}|$ x $m^{e\mu}$. The highest invariant mass bin of the two distributions contains the overflows.

Covariance Matrix (statistical uncertainties only) for the double differential normalised cross-section for $|\Delta\phi^{e\mu}|$ x $m^{e\mu}$. The highest invariant mass bin of the two distributions contains the overflows.

Data bootstrap post unfolding for the double differential normalised cross-section for $|\Delta\phi^{e\mu}|$ x $m^{e\mu}$. The replicas are obtained by reweighting each observed data event by a random integer generated according to Poisson statistics, using the BootstrapGenerator software package (https://gitlab.cern.ch/atlas-physics/sm/StandardModelTools_BootstrapGenerator/BootstrapGenerator), which implements a technique described in ATL-PHYS-PUB-2021-011 (https://cds.cern.ch/record/2759945). The ATLAS event number and run number of each event are used as seed to uniquely but reproducibly initialise the random number generator for each event. All the provided numbers originate from pseudo-data, including the 0th entry, and are in units of [1/GeV/rad]. The highest invariant mass bin of the two distributions contains the overflows.

Data bootstrap post unfolding for the double differential normalised cross-section for $|\Delta\phi^{e\mu}|$ x $m^{e\mu}$. The replicas are obtained by reweighting each observed data event by a random integer generated according to Poisson statistics, using the BootstrapGenerator software package (https://gitlab.cern.ch/atlas-physics/sm/StandardModelTools_BootstrapGenerator/BootstrapGenerator), which implements a technique described in ATL-PHYS-PUB-2021-011 (https://cds.cern.ch/record/2759945). The ATLAS event number and run number of each event are used as seed to uniquely but reproducibly initialise the random number generator for each event. All the provided numbers originate from pseudo-data, including the 0th entry, and are in units of [1/GeV/rad]. The highest invariant mass bin of the two distributions contains the overflows.

Double differential normalised cross-section for $|\Delta\phi^{e\mu}|$ x $\textrm{p}_{\textrm{T}}^{e\mu}$. The impact of the top-quark mass on the cross-section is included in the table and not counted in the total uncertainty entry in the paper. The highest $\textrm{p}_{\textrm{T}}^{e\mu}$ bin contains the overflows.

Double differential normalised cross-section for $|\Delta\phi^{e\mu}|$ x $\textrm{p}_{\textrm{T}}^{e\mu}$. The impact of the top-quark mass on the cross-section is included in the table and not counted in the total uncertainty entry in the paper. The highest $\textrm{p}_{\textrm{T}}^{e\mu}$ bin contains the overflows.

Covariance Matrix (statistical uncertainties only) for the double differential normalised cross-section for $|\Delta\phi^{e\mu}|$ x $\textrm{p}_{\textrm{T}}^{e\mu}$. The highest $\textrm{p}_{\textrm{T}}^{e\mu}$ bin contains the overflows.

Covariance Matrix (statistical uncertainties only) for the double differential normalised cross-section for $|\Delta\phi^{e\mu}|$ x $\textrm{p}_{\textrm{T}}^{e\mu}$. The highest $\textrm{p}_{\textrm{T}}^{e\mu}$ bin contains the overflows.

Data bootstrap post unfolding for the double differential normalised cross-section for $|\Delta\phi^{e\mu}|$ x $\textrm{p}_{\textrm{T}}^{e\mu}$. The replicas are obtained by reweighting each observed data event by a random integer generated according to Poisson statistics, using the BootstrapGenerator software package (https://gitlab.cern.ch/atlas-physics/sm/StandardModelTools_BootstrapGenerator/BootstrapGenerator), which implements a technique described in ATL-PHYS-PUB-2021-011 (https://cds.cern.ch/record/2759945). The ATLAS event number and run number of each event are used as seed to uniquely but reproducibly initialise the random number generator for each event. All the provided numbers originate from pseudo-data, including the 0th entry, and are in units of [1/GeV/rad]. The highest $\textrm{p}_{\textrm{T}}^{e\mu}$ bin contains the overflows.

Data bootstrap post unfolding for the double differential normalised cross-section for $|\Delta\phi^{e\mu}|$ x $\textrm{p}_{\textrm{T}}^{e\mu}$. The replicas are obtained by reweighting each observed data event by a random integer generated according to Poisson statistics, using the BootstrapGenerator software package (https://gitlab.cern.ch/atlas-physics/sm/StandardModelTools_BootstrapGenerator/BootstrapGenerator), which implements a technique described in ATL-PHYS-PUB-2021-011 (https://cds.cern.ch/record/2759945). The ATLAS event number and run number of each event are used as seed to uniquely but reproducibly initialise the random number generator for each event. All the provided numbers originate from pseudo-data, including the 0th entry, and are in units of [1/GeV/rad]. The highest $\textrm{p}_{\textrm{T}}^{e\mu}$ bin contains the overflows.

Double differential normalised cross-section for $|\Delta\phi^{e\mu}|$ x $\textrm{E}^{e} + \textrm{E}^{\mu}$. The impact of the top-quark mass on the cross-section is included in the table and not counted in the total uncertainty entry in the paper. The highest $\textrm{E}^{e} + \textrm{E}^{\mu}$ bin contains the overflows.

Double differential normalised cross-section for $|\Delta\phi^{e\mu}|$ x $\textrm{E}^{e} + \textrm{E}^{\mu}$. The impact of the top-quark mass on the cross-section is included in the table and not counted in the total uncertainty entry in the paper. The highest $\textrm{E}^{e} + \textrm{E}^{\mu}$ bin contains the overflows.

Covariance Matrix (statistical uncertainties only) for the double differential normalised cross-section for $|\Delta\phi^{e\mu}|$ x $\textrm{E}^{e} + \textrm{E}^{\mu}$. The highest $\textrm{E}^{e} + \textrm{E}^{\mu}$ bin contains the overflows.

Covariance Matrix (statistical uncertainties only) for the double differential normalised cross-section for $|\Delta\phi^{e\mu}|$ x $\textrm{E}^{e} + \textrm{E}^{\mu}$. The highest $\textrm{E}^{e} + \textrm{E}^{\mu}$ bin contains the overflows.

Data bootstrap post unfolding for the double differential normalised cross-section for $|\Delta\phi^{e\mu}|$ x $\textrm{E}^{e} + \textrm{E}^{\mu}$. The replicas are obtained by reweighting each observed data event by a random integer generated according to Poisson statistics, using the BootstrapGenerator software package (https://gitlab.cern.ch/atlas-physics/sm/StandardModelTools_BootstrapGenerator/BootstrapGenerator), which implements a technique described in ATL-PHYS-PUB-2021-011 (https://cds.cern.ch/record/2759945). The ATLAS event number and run number of each event are used as seed to uniquely but reproducibly initialise the random number generator for each event. All the provided numbers originate from pseudo-data, including the 0th entry, and are in units of [1/GeV/rad]. The highest $\textrm{E}^{e} + \textrm{E}^{\mu}$ bin contains the overflows.

Data bootstrap post unfolding for the double differential normalised cross-section for $|\Delta\phi^{e\mu}|$ x $\textrm{E}^{e} + \textrm{E}^{\mu}$. The replicas are obtained by reweighting each observed data event by a random integer generated according to Poisson statistics, using the BootstrapGenerator software package (https://gitlab.cern.ch/atlas-physics/sm/StandardModelTools_BootstrapGenerator/BootstrapGenerator), which implements a technique described in ATL-PHYS-PUB-2021-011 (https://cds.cern.ch/record/2759945). The ATLAS event number and run number of each event are used as seed to uniquely but reproducibly initialise the random number generator for each event. All the provided numbers originate from pseudo-data, including the 0th entry, and are in units of [1/GeV/rad]. The highest $\textrm{E}^{e} + \textrm{E}^{\mu}$ bin contains the overflows.

Inclusive cross section predictions up to QCD NNLO accuracy using the NNPDF3.1 PDF set

Measured production cross sections ratio $\sigma(W^+ + \overline{c})$ / $\sigma(W^- + c)$

Predictions for the cross sections ratio $\sigma(W^+ + \overline{c})$ / $\sigma(W^- + c)$ at QCD NLO accuracy from MCFM using different PDF sets

Predictions for the cross sections ratio $\sigma(W^+ + \overline{c})$ / $\sigma(W^- + c)$ up to QCD NNLO accuracy using the NNPDF3.1 PDF set

Measured diferential cross sections $\sigma(W^- + c) + \sigma(W^+ + \overline{c})$ as a function of the absolute value of the pseudorapidity of the lepton from the W decay, unfolded to the particle level.

Measured diferential cross sections as a function of the transverse momentum of the lepton from the W decay, unfolded to the particle level.

Measured diferential cross sections $\sigma(W^- + c) + \sigma(W^+ + \overline{c})$ as a function of the absolute value of the pseudorapidity of the lepton from the W decay, unfolded to the parton level.

Measured diferential cross sections as a function of the transverse momentum of the lepton from the W decay, unfolded to the parton level.

Measured diferential cross sections ratio $R=\sigma(W^+ + \overline{c}) / \sigma(W^- + c)$ as a function of the absolute value of the pseudorapidity of the lepton from the W decay.

Measured diferential cross sections ratio $R=\sigma(W^+ + \overline{c}) / \sigma(W^- + c)$ as a function of the transverse momentum of the lepton from the W decay.

The rejection factor (inverse of the efficiency) for jets that have the other origins (background jets) are shown. The background jet rejection factor is calculated using pre-selected large-$R$ jets in the sample of simulated $Z\rightarrow\nu\nu$ + jets events, dominated by initial state radiation (ISR) jets. The efficiency correction factors are applied on the background rejection for the $W_{qq}$/$Z_{qq}$-tagging. The systematic uncertainty is represented by the hashed bands.

The boson-tagging efficiency for jets arising from $Z/h$ bosons decaying into $b\bar{b}$ (signal jets). The signal jet efficiency of $Z_{bb}$/$h_{bb}$-tagging is evaluated using a sample of pre-selected large-$R$ jets ($p_{\textrm{T}}>200~\textrm{GeV}, |\eta|<2.0, m_{J} > 40~\textrm{GeV}$) in the simulated $(\tilde{W},\tilde{B})$ simplified model signal events with $\Delta m(\tilde{\chi}_{\textrm{heavy}},~\tilde{\chi}_{\textrm{light}}) \ge 400~\textrm{GeV}$. The jets are matched with generator-level $Z/h$-bosons by $\Delta R<1.0$ which decay into $b\bar{b}$. The systematic uncertainty is represented by the hashed bands.

The rejection factor (inverse of the efficiency) for jets that have the other origins (background jets) are shown. The background jet rejection factor is calculated using pre-selected large-$R$ jets in the sample of simulated $Z\rightarrow\nu\nu$ + jets events, dominated by initial state radiation (ISR) jets. As for the $Z_{bb}$/$h_{bb}$-tagging, the rejection is shown as the function of number of $b$- or $c$-quarks contained in the large-$R$ jet within $\Delta R<1.0$. The systematic uncertainty is represented by the hashed bands.

The signal jet efficiency and background jet rejection for the $W_{qq}$-/$Z_{qq}$-tagging working point used in the analysis in comparison with those for the nominal working points (''50% efficiency W/Z tagger'' in ATL-PHYS-PUB-2020-017). The definition of the efficiency and rejection follow the caption of Figure 4. Note that the efficiency for the nominal working point is not 50% since they are tuned using jet samples with different pre-selection.

The signal jet efficiency and background jet rejection for the $W_{qq}$-/$Z_{qq}$-tagging working point used in the analysis in comparison with those for the nominal working points (''50% efficiency W/Z tagger'' in ATL-PHYS-PUB-2020-017). The definition of the efficiency and rejection follow the caption of Figure 4. Note that the efficiency for the nominal working point is not 50% since they are tuned using jet samples with different pre-selection.

The signal jet efficiency and background jet rejection for the $W_{qq}$-/$Z_{qq}$-tagging working point used in the analysis in comparison with those for the nominal working points (''50% efficiency W/Z tagger'' in ATL-PHYS-PUB-2020-017). The definition of the efficiency and rejection follow the caption of Figure 4. Note that the efficiency for the nominal working point is not 50% since they are tuned using jet samples with different pre-selection.

The signal jet efficiency and background jet rejection for the $W_{qq}$-/$Z_{qq}$-tagging working point used in the analysis in comparison with those for the nominal working points (''50% efficiency W/Z tagger'' in ATL-PHYS-PUB-2020-017). The definition of the efficiency and rejection follow the caption of Figure 4. Note that the efficiency for the nominal working point is not 50% since they are tuned using jet samples with different pre-selection.

The total post-fit uncertainty in each of the SRs and VRs.

Summary of the observed data and predicted SM background in all SRs. The background prediction in SR-4Q (SR-2B2Q) is obtained by a background-only fit to CR0L-4Q (CR0L-2B2Q). The total systematic uncertainty on the background prediction is shown by the hatched area. Distributions of a few representative signals are overlaid. A few representative signals are overlaid. For the $(\tilde{W},\tilde{B})$ simplified model models, the label $(x,y)$ indicates $(m(\tilde{\chi}_{1}^{\pm}), m(\tilde{\chi}_{1}^{0}))=(x,y)~\textrm{GeV}$. The bottom panel shows the statistical significance of the discrepancy between the observed number of events and the SM expectation.

Number of observed data events and the SM backgrounds in the SRs and the CR0L bins. The SM backgrounds are predicted by the background-only fits. Note that the relative uncertainty on the expected yield in the CRs between the reducible backgrounds are identical, given that a common normalization factor is assigned for all of them in the fit. Yields for negligible backgrounds are indicated by "0.0000001±0.0000001" in the table.

Number of observed data events and the post-fit SM background prediction in the VR1L(1Y) bins and the corresponding CR1L(1Y) bins. "-" indicates negligibly small contribution. Note that the relative uncertainty on the expected yield in the CRs between the reducible backgrounds are identical, given that a common normalization factor is assigned for all of them in the fit. Yields for negligible backgrounds are indicated by "0.0000001±0.0000001" in the table.

$m_{\textrm{eff}}$ distribution in SR-4Q-VV. The post-fit SM background expectation using the background-only fit is shown in a histogram stack. Distributions of a few representative signals are overlaid. The bottom panels show the ratio of the observed data to the background prediction. The selection criteria on the variable shown by each plot is removed, while the arrow indicates the cut value to define the region. A few representative signals are overlaid. For the $(\tilde{W},\tilde{B})$ simplified model models, the label $(x,y)$ indicates $(m(\tilde{\chi}_{1}^{\pm}), m(\tilde{\chi}_{1}^{0}))=(x,y)~\text{GeV}$.

Distribution of $m(J_{1})$ in SR-4Q-VV. The post-fit SM background expectation using the background-only fit is shown in a histogram stack. Distributions of relevant signals are overlaid. The bottom panels show the ratio of the observed data to the background prediction. The selection criteria on the variable shown by each plot is removed, while the arrow indicates the cut value to define the region. A few representative signals are overlaid. For the $(\tilde{W},\tilde{B})$ simplified model models, the label $(x,y)$ indicates $(m(\tilde{\chi}_{1}^{\pm}), m(\tilde{\chi}_{1}^{0}))=(x,y)~\text{GeV}$.

Distribution of $D_2(J_{1})$ in SR-4Q-VV. For $D_2$, the cut value applied for $V_{qq}$-tagging ($D_{2,\text{cut}}$) is subtracted as the off-set so that $D_2-D_{2,\text{cut}}(J)<0$ represents $J$ passing the $D_2$ selection. The post-fit SM background expectation using the background-only fit is shown in a histogram stack. Distributions of relevant signals are overlaid. The bottom panels show the ratio of the observed data to the background prediction. The selection criteria on the variable shown by each plot is removed, while the arrow indicates the cut value to define the region. A few representative signals are overlaid. For the $(\tilde{W},\tilde{B})$ simplified model models, the label $(x,y)$ indicates $(m(\tilde{\chi}_{1}^{\pm}), m(\tilde{\chi}_{1}^{0}))=(x,y)~\text{GeV}$.

Distribution of $m(J_{2})$ in SR-4Q-VV. The post-fit SM background expectation using the background-only fit is shown in a histogram stack. Distributions of relevant signals are overlaid. The bottom panels show the ratio of the observed data to the background prediction. The selection criteria on the variable shown by each plot is removed, while the arrow indicates the cut value to define the region. A few representative signals are overlaid. For the $(\tilde{W},\tilde{B})$ simplified model models, the label $(x,y)$ indicates $(m(\tilde{\chi}_{1}^{\pm}), m(\tilde{\chi}_{1}^{0}))=(x,y)~\text{GeV}$.

Distribution of $D_2(J_{2})$ in SR-4Q-VV. For $D_2$, the cut value applied for $V_{qq}$-tagging ($D_{2,\text{cut}}$) is subtracted as the off-set so that $D_2-D_{2,\text{cut}}(J)<0$ represents $J$ passing the $D_2$ selection. The post-fit SM background expectation using the background-only fit is shown in a histogram stack. Distributions of relevant signals are overlaid. The bottom panels show the ratio of the observed data to the background prediction. The selection criteria on the variable shown by each plot is removed, while the arrow indicates the cut value to define the region. A few representative signals are overlaid. For the $(\tilde{W},\tilde{B})$ simplified model models, the label $(x,y)$ indicates $(m(\tilde{\chi}_{1}^{\pm}), m(\tilde{\chi}_{1}^{0}))=(x,y)~\text{GeV}$.

$m_{\textrm{T2}}$ distributions in SR-2B2Q-VZ. The post-fit SM background expectation using the background-only fit is shown in a histogram stack. Distributions of a few representative signals are overlaid. The bottom panels show the ratio of the observed data to the background prediction. The selection criteria on the variable shown by each plot is removed, while the arrow indicates the cut value to define the region. A few representative signals are overlaid. For the $(\tilde{W},\tilde{B})$ simplified model models, the label $(x,y)$ indicates $(m(\tilde{\chi}_{1}^{\pm}), m(\tilde{\chi}_{1}^{0}))=(x,y)~\text{GeV}$.

Distribution of $m(J_{bb})$ in SR-2B2Q-VZ. The post-fit SM background expectation using the background-only fit is shown in a histogram stack. Distributions of relevant signals are overlaid. The bottom panels show the ratio of the observed data to the background prediction. The selection criteria on the variable shown by each plot is removed, while the arrow indicates the cut value to define the region. A few representative signals are overlaid. For the $(\tilde{W},\tilde{B})$ simplified model models, the label $(x,y)$ indicates $(m(\tilde{\chi}_{1}^{\pm}), m(\tilde{\chi}_{1}^{0}))=(x,y)~\text{GeV}$.

Distribution of $m_{\textrm{eff}} $ in SR-2B2Q-VZ. The post-fit SM background expectation using the background-only fit is shown in a histogram stack. Distributions of relevant signals are overlaid. The bottom panels show the ratio of the observed data to the background prediction. The selection criteria on the variable shown by each plot is removed, while the arrow indicates the cut value to define the region. A few representative signals are overlaid. For the $(\tilde{W},\tilde{B})$ simplified model models, the label $(x,y)$ indicates $(m(\tilde{\chi}_{1}^{\pm}), m(\tilde{\chi}_{1}^{0}))=(x,y)~\text{GeV}$.

$m_{\textrm{T2}}$ distributions in SR-2B2Q-Vh. The post-fit SM background expectation using the background-only fit is shown in a histogram stack. Distributions of a few representative signals are overlaid. The bottom panels show the ratio of the observed data to the background prediction. The selection criteria on the variable shown by each plot is removed, while the arrow indicates the cut value to define the region. A few representative signals are overlaid. For the $(\tilde{W},\tilde{B})$ simplified model models, the label $(x,y)$ indicates $(m(\tilde{\chi}_{1}^{\pm}), m(\tilde{\chi}_{1}^{0}))=(x,y)~\text{GeV}$.

Distribution of $m(J_{bb})$ in SR-2B2Q-Vh. The post-fit SM background expectation using the background-only fit is shown in a histogram stack. Distributions of relevant signals are overlaid. The bottom panels show the ratio of the observed data to the background prediction. The selection criteria on the variable shown by each plot is removed, while the arrow indicates the cut value to define the region. A few representative signals are overlaid. For the $(\tilde{W},\tilde{B})$ simplified model models, the label $(x,y)$ indicates $(m(\tilde{\chi}_{1}^{\pm}), m(\tilde{\chi}_{1}^{0}))=(x,y)~\text{GeV}$.

Distribution of $m_{\textrm{eff}} $ in SR-2B2Q-Vh. The post-fit SM background expectation using the background-only fit is shown in a histogram stack. Distributions of relevant signals are overlaid. The bottom panels show the ratio of the observed data to the background prediction. The selection criteria on the variable shown by each plot is removed, while the arrow indicates the cut value to define the region. A few representative signals are overlaid. For the $(\tilde{W},\tilde{B})$ simplified model models, the label $(x,y)$ indicates $(m(\tilde{\chi}_{1}^{\pm}), m(\tilde{\chi}_{1}^{0}))=(x,y)~\text{GeV}$.

Exclusion limits for $(\tilde{W},\tilde{B})$ simplifiec model as a function of the produced wino mass ($m(\tilde{\chi}_{1}^{\pm}/\tilde{\chi}_{2}^{0})$) and the bino LSP mass ($m(\tilde{\chi}_{1}^{0})$). Expected (dashed) and observed (solid) 95% CL exclusion limits on $(\tilde{W},\tilde{B})$ simplified models; C1C1-WW. The expected limit for the 1$\sigma$ down variation is not shown as no mass points could be excluded.

Exclusion limits for $(\tilde{W},\tilde{B})$ simplifiec model as a function of the produced wino mass ($m(\tilde{\chi}_{1}^{\pm}/\tilde{\chi}_{2}^{0})$) and the bino LSP mass ($m(\tilde{\chi}_{1}^{0})$). Expected (dashed) and observed (solid) 95% CL exclusion limits on $(\tilde{W},\tilde{B})$ simplified models; C1C1-WW. The expected limit for the 1$\sigma$ down variation is not shown as no mass points could be excluded.

Exclusion limits for $(\tilde{W},\tilde{B})$ simplifiec model as a function of the produced wino mass ($m(\tilde{\chi}_{1}^{\pm}/\tilde{\chi}_{2}^{0})$) and the bino LSP mass ($m(\tilde{\chi}_{1}^{0})$). Expected (dashed) and observed (solid) 95% CL exclusion limits on $(\tilde{W},\tilde{B})$ simplified models; C1C1-WW. The expected limit for the 1$\sigma$ down variation is not shown as no mass points could be excluded.

Exclusion limits for $(\tilde{W},\tilde{B})$ simplifiec model as a function of the produced wino mass ($m(\tilde{\chi}_{1}^{\pm}/\tilde{\chi}_{2}^{0})$) and the bino LSP mass ($m(\tilde{\chi}_{1}^{0})$). Expected (dashed) and observed (solid) 95% CL exclusion limits on $(\tilde{W},\tilde{B})$ simplified models; C1C1-WW. The expected limit for the 1$\sigma$ down variation is not shown as no mass points could be excluded.

Exclusion limits for $(\tilde{W},\tilde{B})$ simplifiec model as a function of the produced wino mass ($m(\tilde{\chi}_{1}^{\pm}/\tilde{\chi}_{2}^{0})$) and the bino LSP mass ($m(\tilde{\chi}_{1}^{0})$). Expected (dashed) and observed (solid) 95% CL exclusion limits on $(\tilde{W},\tilde{B})$ simplified models; C1C1-WW. The expected limit for the 1$\sigma$ down variation is not shown as no mass points could be excluded.

Exclusion limits for $(\tilde{W},\tilde{B})$ simplifiec model as a function of the produced wino mass ($m(\tilde{\chi}_{1}^{\pm}/\tilde{\chi}_{2}^{0})$) and the bino LSP mass ($m(\tilde{\chi}_{1}^{0})$). Expected (dashed) and observed (solid) 95% CL exclusion limits on $(\tilde{W},\tilde{B})$ simplified models; C1N2-WZ.

Exclusion limits for $(\tilde{W},\tilde{B})$ simplifiec model as a function of the produced wino mass ($m(\tilde{\chi}_{1}^{\pm}/\tilde{\chi}_{2}^{0})$) and the bino LSP mass ($m(\tilde{\chi}_{1}^{0})$). Expected (dashed) and observed (solid) 95% CL exclusion limits on $(\tilde{W},\tilde{B})$ simplified models; C1N2-WZ.

Exclusion limits for $(\tilde{W},\tilde{B})$ simplifiec model as a function of the produced wino mass ($m(\tilde{\chi}_{1}^{\pm}/\tilde{\chi}_{2}^{0})$) and the bino LSP mass ($m(\tilde{\chi}_{1}^{0})$). Expected (dashed) and observed (solid) 95% CL exclusion limits on $(\tilde{W},\tilde{B})$ simplified models; C1N2-WZ.

Exclusion limits for $(\tilde{W},\tilde{B})$ simplifiec model as a function of the produced wino mass ($m(\tilde{\chi}_{1}^{\pm}/\tilde{\chi}_{2}^{0})$) and the bino LSP mass ($m(\tilde{\chi}_{1}^{0})$). Expected (dashed) and observed (solid) 95% CL exclusion limits on $(\tilde{W},\tilde{B})$ simplified models; C1N2-WZ.

Exclusion limits for $(\tilde{W},\tilde{B})$ simplifiec model as a function of the produced wino mass ($m(\tilde{\chi}_{1}^{\pm}/\tilde{\chi}_{2}^{0})$) and the bino LSP mass ($m(\tilde{\chi}_{1}^{0})$). Expected (dashed) and observed (solid) 95% CL exclusion limits on $(\tilde{W},\tilde{B})$ simplified models; C1N2-WZ.

Exclusion limits for $(\tilde{W},\tilde{B})$ simplifiec model as a function of the produced wino mass ($m(\tilde{\chi}_{1}^{\pm}/\tilde{\chi}_{2}^{0})$) and the bino LSP mass ($m(\tilde{\chi}_{1}^{0})$). Expected (dashed) and observed (solid) 95% CL exclusion limits on $(\tilde{W},\tilde{B})$ simplified models; C1N2-WZ.

Exclusion limits for $(\tilde{W},\tilde{B})$ simplifiec model as a function of the produced wino mass ($m(\tilde{\chi}_{1}^{\pm}/\tilde{\chi}_{2}^{0})$) and the bino LSP mass ($m(\tilde{\chi}_{1}^{0})$). Expected (dashed) and observed (solid) 95% CL exclusion limits on $(\tilde{W},\tilde{B})$ simplified models; C1N2-Wh.

Exclusion limits for $(\tilde{W},\tilde{B})$ simplifiec model as a function of the produced wino mass ($m(\tilde{\chi}_{1}^{\pm}/\tilde{\chi}_{2}^{0})$) and the bino LSP mass ($m(\tilde{\chi}_{1}^{0})$). Expected (dashed) and observed (solid) 95% CL exclusion limits on $(\tilde{W},\tilde{B})$ simplified models; C1N2-Wh.