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;).

Transverse momentum distribution for the second leading lepton for events in the SR2 signal region for DATA and SM predictions.

Transverse momentum distribution for the third leading lepton for events in the SR1 signal region for DATA and SM predictions.

Transverse momentum distribution for the third leading lepton for events in the SR2 signal region for DATA and SM predictions.

Missing transverse energy for events in the SR1 signal region for DATA and SM predictions.

Missing transverse energy for events in the SR2 signal region for DATA and SM predictions.

Invariant mass of the same-flavour-opposite-sign (SFOS) lepton pair for events in the SR1 signal region for DATA and SM predictions.

The Cross Section in signal region SR1 for the SUSY pMSSM model with M1=100 GeV grid.

The Cross Section in signal region SR1 for the SUSY simplified model grid.

The Number of generated Events in signal region SR1 for the SUSY pMSSM model with M1=100 GeV grid.

The Number of generated Events in signal region SR1 for the SUSY simplified model grid.

The Efficiency in signal region SR1 for the SUSY pMSSM model with M1=100 GeV grid.

The Efficiency in signal region SR1 for the SUSY simplified model grid.

The Acceptance in signal region SR1 for the SUSY pMSSM model with M1=100 GeV grid.

The Acceptance in signal region SR1 for the SUSY simplified model grid.

The Acceptance*Efficiency in signal region SR1 for the SUSY pMSSM model with M1=100 GeV grid.

The Acceptance*Efficiency in signal region SR1 for the SUSY simplified model grid.

The Systematic Uncertainty of the data (excluding the Monte Carlo) in signal region SR1 for the SUSY pMSSM model with M1=100 GeV grid.

The Systematic Uncertainty of the data (excluding the Monte Carlo) in signal region SR1 for the SUSY simplified model grid.

CL values for the pMSSM with M1=100 GeV model grid for the SR1 signal region.

CL values for the simplified model model grid for the SR1 signal region.

Post-fit SR yields of the high-mass channel. The upper panel shows the observed number of events, as well the post-fit background predictions in each region. The bottom panel shows the ratio of the observed data and the total background prediction. The shaded areas correspond to the total statistical and systematic uncertainties obtained after the fit and described in Section 6.

Post-fit SR yields of the high-mass channel. The upper panel shows the observed number of events, as well the post-fit background predictions in each region. The bottom panel shows the ratio of the observed data and the total background prediction. The shaded areas correspond to the total statistical and systematic uncertainties obtained after the fit and described in Section 6.

Post-fit SR yields of the high-mass channel. The upper panel shows the observed number of events, as well the post-fit background predictions in each region. The bottom panel shows the ratio of the observed data and the total background prediction. The shaded areas correspond to the total statistical and systematic uncertainties obtained after the fit and described in Section 6.

Post-fit SR yields of the high-mass channel. The upper panel shows the observed number of events, as well the post-fit background predictions in each region. The bottom panel shows the ratio of the observed data and the total background prediction. The shaded areas correspond to the total statistical and systematic uncertainties obtained after the fit and described in Section 6.

Post-fit SR yields of the high-mass channel. The upper panel shows the observed number of events, as well the post-fit background predictions in each region. The bottom panel shows the ratio of the observed data and the total background prediction. The shaded areas correspond to the total statistical and systematic uncertainties obtained after the fit and described in Section 6.

Post-fit SR yields of the high-mass channel. The upper panel shows the observed number of events, as well the post-fit background predictions in each region. The bottom panel shows the ratio of the observed data and the total background prediction. The shaded areas correspond to the total statistical and systematic uncertainties obtained after the fit and described in Section 6.

Post-fit SR yields of the high-mass channel. The upper panel shows the observed number of events, as well the post-fit background predictions in each region. The bottom panel shows the ratio of the observed data and the total background prediction. The shaded areas correspond to the total statistical and systematic uncertainties obtained after the fit and described in Section 6.

Pre-fit data and background (reweighted $2b$) predictions for each $4b$ SR $E_\text{T}^\text{miss}$ and $m_\text{eff}$ bin of the low-mass channel for the 2016 data-taking period. The bottom panel shows the significance of any differences between the observed $4b$ data and the background prediction. The $1\sigma$ and $2\sigma$ bands are shown in green and yellow, respectively. All systematics are included except the background normalization, which is 2.3%.

Pre-fit data and background (reweighted $2b$) predictions for each $4b$ SR $E_\text{T}^\text{miss}$ and $m_\text{eff}$ bin of the low-mass channel for the 2017 data-taking period. The bottom panel shows the significance of any differences between the observed $4b$ data and the background prediction. The $1\sigma$ and $2\sigma$ bands are shown in green and yellow, respectively. All systematics are included except the background normalization, which is 3.7%.

Pre-fit data and background (reweighted $2b$) predictions for each $4b$ SR $E_\text{T}^\text{miss}$ and $m_\text{eff}$ bin of the low-mass channel for the 2018 data-taking period. The bottom panel shows the significance of any differences between the observed $4b$ data and the background prediction. The $1\sigma$ and $2\sigma$ bands are shown in green and yellow, respectively. All systematics are included except the background normalization, which is 1.8%.

Exclusion limits of the low-mass and high-mass channels. The low-mass channel is used for $m_{\tilde{H}}<250$ GeV while the high-mass channel is used for $m_{\tilde{H}}\ge250$ GeV. The plot shows the observed (solid) and expected (dashed) 95% CL upper limits on the cross section of higgsino pair production, assuming a higgsino decay branching ratio of $\mathcal{B}(\tilde{H}\rightarrow h + \tilde{G})=100\%$. The theory cross section and its uncertainty are shown by the solid red line and red shading. Results from a previous ATLAS search using 24.3-36.1 fb$^{-1}$ [13] are shown by the solid (observed) and dashed (expected) blue lines. The bottom panel shows the ratio of the limits to the theory cross section. The phase space above the lines is excluded.

Exclusion limits of the low-mass and high-mass channels. The low-mass channel is used for $m_{\tilde{H}}<250$ GeV while the high-mass channel is used for $m_{\tilde{H}}\ge250$ GeV. The plot shows the observed (solid) and expected (dashed) 95% CL upper limits on the cross section of higgsino pair production, assuming a higgsino decay branching ratio of $\mathcal{B}(\tilde{H}\rightarrow h + \tilde{G})=100\%$. The theory cross section and its uncertainty are shown by the solid red line and red shading. Results from a previous ATLAS search using 24.3-36.1 fb$^{-1}$ [13] are shown by the solid (observed) and dashed (expected) blue lines. The bottom panel shows the ratio of the limits to the theory cross section. The phase space above the lines is excluded.

Exclusion limits of the low-mass and high-mass channels. The low-mass channel is used for $m_{\tilde{H}}<250$ GeV while the high-mass channel is used for $m_{\tilde{H}}\ge250$ GeV. The plot shows the observed (solid) and expected (dashed) 95% CL upper limits on the cross section of higgsino pair production, assuming a higgsino decay branching ratio of $\mathcal{B}(\tilde{H}\rightarrow h + \tilde{G})=100\%$. The theory cross section and its uncertainty are shown by the solid red line and red shading. Results from a previous ATLAS search using 24.3-36.1 fb$^{-1}$ [13] are shown by the solid (observed) and dashed (expected) blue lines. The bottom panel shows the ratio of the limits to the theory cross section. The phase space above the lines is excluded.

Exclusion limits of the low-mass and high-mass channels. The low-mass channel is used for $m_{\tilde{H}}<250$ GeV while the high-mass channel is used for $m_{\tilde{H}}\ge250$ GeV. The plot shows the observed (solid) and expected (dashed) 95% CL upper limits on the cross section of higgsino pair production, assuming a higgsino decay branching ratio of $\mathcal{B}(\tilde{H}\rightarrow h + \tilde{G})=100\%$. The theory cross section and its uncertainty are shown by the solid red line and red shading. Results from a previous ATLAS search using 24.3-36.1 fb$^{-1}$ [13] are shown by the solid (observed) and dashed (expected) blue lines. The bottom panel shows the ratio of the limits to the theory cross section. The phase space above the lines is excluded.

Exclusion limits of the low-mass and high-mass channels. The low-mass channel is used for $m_{\tilde{H}}<250$ GeV while the high-mass channel is used for $m_{\tilde{H}}\ge250$ GeV. The plot shows the observed (solid) and expected (dashed) 95% CL upper limits on the cross section of higgsino pair production, assuming a higgsino decay branching ratio of $\mathcal{B}(\tilde{H}\rightarrow h + \tilde{G})=100\%$. The theory cross section and its uncertainty are shown by the solid red line and red shading. Results from a previous ATLAS search using 24.3-36.1 fb$^{-1}$ [13] are shown by the solid (observed) and dashed (expected) blue lines. The bottom panel shows the ratio of the limits to the theory cross section. The phase space above the lines is excluded.

Exclusion limits of the low-mass and high-mass channels. The low-mass channel is used for $m_{\tilde{H}}<250$ GeV while the high-mass channel is used for $m_{\tilde{H}}\ge250$ GeV. The plot shows the observed (solid) and expected (dashed) 95% CL upper limits on the cross section of higgsino pair production, assuming a higgsino decay branching ratio of $\mathcal{B}(\tilde{H}\rightarrow h + \tilde{G})=100\%$. The theory cross section and its uncertainty are shown by the solid red line and red shading. Results from a previous ATLAS search using 24.3-36.1 fb$^{-1}$ [13] are shown by the solid (observed) and dashed (expected) blue lines. The bottom panel shows the ratio of the limits to the theory cross section. The phase space above the lines is excluded.

Exclusion limits of the low-mass and high-mass channels. The low-mass channel is used for $m_{\tilde{H}}<250$ GeV while the high-mass channel is used for $m_{\tilde{H}}\ge250$ GeV. The plot shows the observed (solid) and expected (dashed) 95% CL upper limits on the cross section of higgsino pair production, assuming a higgsino decay branching ratio of $\mathcal{B}(\tilde{H}\rightarrow h + \tilde{G})=100\%$. The theory cross section and its uncertainty are shown by the solid red line and red shading. Results from a previous ATLAS search using 24.3-36.1 fb$^{-1}$ [13] are shown by the solid (observed) and dashed (expected) blue lines. The bottom panel shows the ratio of the limits to the theory cross section. The phase space above the lines is excluded.

Exclusion limits of the low-mass and high-mass channels. The low-mass channel is used for $m_{\tilde{H}}<250$ GeV while the high-mass channel is used for $m_{\tilde{H}}\ge250$ GeV. The plot shows the observed (solid) and expected (dashed) 95% CL upper limits on the cross section of higgsino pair production, assuming a higgsino decay branching ratio of $\mathcal{B}(\tilde{H}\rightarrow h + \tilde{G})=100\%$. The theory cross section and its uncertainty are shown by the solid red line and red shading. Results from a previous ATLAS search using 24.3-36.1 fb$^{-1}$ [13] are shown by the solid (observed) and dashed (expected) blue lines. The bottom panel shows the ratio of the limits to the theory cross section. The phase space above the lines is excluded.

Exclusion limits of the low-mass and high-mass channels. The low-mass channel is used for $m_{\tilde{H}}<250$ GeV while the high-mass channel is used for $m_{\tilde{H}}\ge250$ GeV. The plot shows the observed (solid) and expected (dashed) 95% CL upper limits on the cross section of higgsino pair production, assuming a higgsino decay branching ratio of $\mathcal{B}(\tilde{H}\rightarrow h + \tilde{G})=100\%$. The theory cross section and its uncertainty are shown by the solid red line and red shading. Results from a previous ATLAS search using 24.3-36.1 fb$^{-1}$ [13] are shown by the solid (observed) and dashed (expected) blue lines. The bottom panel shows the ratio of the limits to the theory cross section. The phase space above the lines is excluded.

Exclusion limits of the low-mass and high-mass channels. The low-mass channel is used for $m_{\tilde{H}}<250$ GeV while the high-mass channel is used for $m_{\tilde{H}}\ge250$ GeV. The plot shows the observed (solid) and expected (dashed) 95% CL upper limits on the cross section of higgsino pair production, assuming a higgsino decay branching ratio of $\mathcal{B}(\tilde{H}\rightarrow h + \tilde{G})=100\%$. The theory cross section and its uncertainty are shown by the solid red line and red shading. Results from a previous ATLAS search using 24.3-36.1 fb$^{-1}$ [13] are shown by the solid (observed) and dashed (expected) blue lines. The bottom panel shows the ratio of the limits to the theory cross section. The phase space above the lines is excluded.

Exclusion limits of the low-mass and high-mass channels. The low-mass channel is used for $m_{\tilde{H}}<250$ GeV while the high-mass channel is used for $m_{\tilde{H}}\ge250$ GeV. The plot shows the 95% CL observed (solid) and expected (dashed) upper limits on $\mathcal{B}(\tilde{H}\rightarrow h + \tilde{G})$, assuming the theory cross section for higgsino pair production. The higgsinos are assumed to decay as $\tilde{H}\rightarrow h + \tilde{G}$ or $\tilde{H}\rightarrow Z + \tilde{G}$. The phase space above the lines is excluded.

Exclusion limits of the low-mass and high-mass channels. The low-mass channel is used for $m_{\tilde{H}}<250$ GeV while the high-mass channel is used for $m_{\tilde{H}}\ge250$ GeV. The plot shows the 95% CL observed (solid) and expected (dashed) upper limits on $\mathcal{B}(\tilde{H}\rightarrow h + \tilde{G})$, assuming the theory cross section for higgsino pair production. The higgsinos are assumed to decay as $\tilde{H}\rightarrow h + \tilde{G}$ or $\tilde{H}\rightarrow Z + \tilde{G}$. The phase space above the lines is excluded.

Exclusion limits of the low-mass and high-mass channels. The low-mass channel is used for $m_{\tilde{H}}<250$ GeV while the high-mass channel is used for $m_{\tilde{H}}\ge250$ GeV. The plot shows the 95% CL observed (solid) and expected (dashed) upper limits on $\mathcal{B}(\tilde{H}\rightarrow h + \tilde{G})$, assuming the theory cross section for higgsino pair production. The higgsinos are assumed to decay as $\tilde{H}\rightarrow h + \tilde{G}$ or $\tilde{H}\rightarrow Z + \tilde{G}$. The phase space above the lines is excluded.

Exclusion limits of the low-mass and high-mass channels. The low-mass channel is used for $m_{\tilde{H}}<250$ GeV while the high-mass channel is used for $m_{\tilde{H}}\ge250$ GeV. The plot shows the 95% CL observed (solid) and expected (dashed) upper limits on $\mathcal{B}(\tilde{H}\rightarrow h + \tilde{G})$, assuming the theory cross section for higgsino pair production. The higgsinos are assumed to decay as $\tilde{H}\rightarrow h + \tilde{G}$ or $\tilde{H}\rightarrow Z + \tilde{G}$. The phase space above the lines is excluded.

Exclusion limits of the low-mass and high-mass channels. The low-mass channel is used for $m_{\tilde{H}}<250$ GeV while the high-mass channel is used for $m_{\tilde{H}}\ge250$ GeV. The plot shows the 95% CL observed (solid) and expected (dashed) upper limits on $\mathcal{B}(\tilde{H}\rightarrow h + \tilde{G})$, assuming the theory cross section for higgsino pair production. The higgsinos are assumed to decay as $\tilde{H}\rightarrow h + \tilde{G}$ or $\tilde{H}\rightarrow Z + \tilde{G}$. The phase space above the lines is excluded.

Exclusion limits of the low-mass channel. The plot shows the observed (solid) and expected (dashed) 95% CL upper limits on the cross section of higgsino pair production, assuming a higgsino decay branching ratio of $\mathcal{B}(\tilde{H}\rightarrow h + \tilde{G})=100\%$. The theory cross section and its uncertainty are shown by the solid red line and red shading. The bottom panel shows the ratio of the limits to the theory cross section. The phase space above the lines is excluded.

Exclusion limits of the low-mass channel. The plot shows the observed (solid) and expected (dashed) 95% CL upper limits on the cross section of higgsino pair production, assuming a higgsino decay branching ratio of $\mathcal{B}(\tilde{H}\rightarrow h + \tilde{G})=100\%$. The theory cross section and its uncertainty are shown by the solid red line and red shading. The bottom panel shows the ratio of the limits to the theory cross section. The phase space above the lines is excluded.

Exclusion limits of the high-mass channel. The plot shows the observed (solid) and expected (dashed) 95% CL upper limits on the cross section of higgsino pair production, assuming a higgsino decay branching ratio of $\mathcal{B}(\tilde{H}\rightarrow h + \tilde{G})=100\%$. The theory cross section and its uncertainty are shown by the solid red line and red shading. The bottom panel shows the ratio of the limits to the theory cross section. The phase space above the lines is excluded.

Exclusion limits of the high-mass channel. The plot shows the observed (solid) and expected (dashed) 95% CL upper limits on the cross section of higgsino pair production, assuming a higgsino decay branching ratio of $\mathcal{B}(\tilde{H}\rightarrow h + \tilde{G})=100\%$. The theory cross section and its uncertainty are shown by the solid red line and red shading. The bottom panel shows the ratio of the limits to the theory cross section. The phase space above the lines is excluded.

Results of the background-only fit in the low-mass channel discovery region SR_LM_150. Both pre-fit and post-fit values are shown.

Results of the background-only fit in the low-mass channel discovery region SR_LM_300. Both pre-fit and post-fit values are shown.

The experimental efficiency of the low-mass channel for the exclusion and discovery signal regions as a function of higgsino mass. The experimental efficiency is defined as the number of events passing the detector-level event selections divided by the number of events passing the event selections for a perfect detector. The denominator is obtained by implementing particle-level event selections that emulate the detector-level selections. This treats the lack of availability of $b$-jet triggers as an inefficiency.

The particle-level acceptance for the low-mass exclusion and discovery signal regions, shown as a function of higgsino mass. The acceptance is defined as the fraction of signal events passing the particle-level event selection that emulates the detector-level selection. The acceptance calculation considers only those signal events where both higgsinos decay to Higgs bosons.

The experimental efficiency of the high-mass channel discovery regions as a function of higgsino mass. For each higgsino mass, the efficiency is shown for the SR-1 region corresponding to the mass. For masses above 1100 GeV, SR-1-1100 is used. The experimental efficiency is defined as the number of events passing the detector-level event selections divided by the number of events passing the event selections for a perfect detector. The denominator is obtained by implementing particle-level event selections that emulate the detector-level selections. The efficiency calculation considers only those signal events where both higgsinos decay to Higgs bosons.

The particle-level acceptance for the high-mass signal regions, shown as a function of higgsino mass. For each higgsino mass, the acceptance is shown for the SR-1 region corresponding to the mass. For masses above 1100 GeV, SR-1-1100 is used. The acceptance is defined as the fraction of signal events passing the particle-level event selection that emulates the detector-level selection. The acceptance calculation considers only those signal events where both higgsinos decay to Higgs bosons.

Cutflow for the low-mass channel for a representative 130 GeV signal. The preselection requires 4 or more jets with $p_\text{T}>25$ GeV and 2 or more $b$-jets with $p_\text{T}>25$ GeV. The $b$-jet cut requires 4 or more $b$-jets with $p_\text{T}>40$ GeV. As the samples are generated with $\mathcal{B}(\tilde{H}\rightarrow h\tilde{G})$=50%, $\mathcal{B}(\tilde{H}\rightarrow Z\tilde{G})$=50% to allow for both decays to be studied, the $hh$ events selection is used to select the events where each of the higgsinos decays to a Higgs boson. Expected yields are normalized to a luminosity of 139 fb$^{-1}$, with the availability of $b$-jet triggers lowering the luminosity to 126 fb$^{-1}$. All selections are cumulative, with the exception of the SR cuts, which are each applied separately.

Cutflow for the low-mass channel for a representative 150 GeV signal. The preselection requires 4 or more jets with $p_\text{T}>25$ GeV and 2 or more $b$-jets with $p_\text{T}>25$ GeV. The $b$-jet cut requires 4 or more $b$-jets with $p_\text{T}>40$ GeV. As the samples are generated with $\mathcal{B}(\tilde{H}\rightarrow h\tilde{G})$=50%, $\mathcal{B}(\tilde{H}\rightarrow Z\tilde{G})$=50% to allow for both decays to be studied, the $hh$ events selection is used to select the events where each of the higgsinos decays to a Higgs boson. Expected yields are normalized to a luminosity of 139 fb$^{-1}$, with the availability of $b$-jet triggers lowering the luminosity to 126 fb$^{-1}$. All selections are cumulative, with the exception of the SR cuts, which are each applied separately.

Cutflow for the low-mass channel for a representative 200 GeV signal. The preselection requires 4 or more jets with $p_\text{T}>25$ GeV and 2 or more $b$-jets with $p_\text{T}>25$ GeV. The $b$-jet cut requires 4 or more $b$-jets with $p_\text{T}>40$ GeV. As the samples are generated with $\mathcal{B}(\tilde{H}\rightarrow h\tilde{G})$=50%, $\mathcal{B}(\tilde{H}\rightarrow Z\tilde{G})$=50% to allow for both decays to be studied, the $hh$ events selection is used to select the events where each of the higgsinos decays to a Higgs boson. Expected yields are normalized to a luminosity of 139 fb$^{-1}$, with the availability of $b$-jet triggers lowering the luminosity to 126 fb$^{-1}$. All selections are cumulative, with the exception of the SR cuts, which are each applied separately.

Cutflow for the low-mass channel for a representative 250 GeV signal. The preselection requires 4 or more jets with $p_\text{T}>25$ GeV and 2 or more $b$-jets with $p_\text{T}>25$ GeV. The $b$-jet cut requires 4 or more $b$-jets with $p_\text{T}>40$ GeV. As the samples are generated with $\mathcal{B}(\tilde{H}\rightarrow h\tilde{G})$=50%, $\mathcal{B}(\tilde{H}\rightarrow Z\tilde{G})$=50% to allow for both decays to be studied, the $hh$ events selection is used to select the events where each of the higgsinos decays to a Higgs boson. Expected yields are normalized to a luminosity of 139 fb$^{-1}$, with the availability of $b$-jet triggers lowering the luminosity to 126 fb$^{-1}$. All selections are cumulative, with the exception of the SR cuts, which are each applied separately.

Cutflow for the low-mass channel for a representative 300 GeV signal. The preselection requires 4 or more jets with $p_\text{T}>25$ GeV and 2 or more $b$-jets with $p_\text{T}>25$ GeV. The $b$-jet cut requires 4 or more $b$-jets with $p_\text{T}>40$ GeV. As the samples are generated with $\mathcal{B}(\tilde{H}\rightarrow h\tilde{G})$=50%, $\mathcal{B}(\tilde{H}\rightarrow Z\tilde{G})$=50% to allow for both decays to be studied, the $hh$ events selection is used to select the events where each of the higgsinos decays to a Higgs boson. Expected yields are normalized to a luminosity of 139 fb$^{-1}$, with the availability of $b$-jet triggers lowering the luminosity to 126 fb$^{-1}$. All selections are cumulative, with the exception of the SR cuts, which are each applied separately.

Cutflow for the low-mass channel for a representative 400 GeV signal. The preselection requires 4 or more jets with $p_\text{T}>25$ GeV and 2 or more $b$-jets with $p_\text{T}>25$ GeV. The $b$-jet cut requires 4 or more $b$-jets with $p_\text{T}>40$ GeV. As the samples are generated with $\mathcal{B}(\tilde{H}\rightarrow h\tilde{G})$=50%, $\mathcal{B}(\tilde{H}\rightarrow Z\tilde{G})$=50% to allow for both decays to be studied, the $hh$ events selection is used to select the events where each of the higgsinos decays to a Higgs boson. Expected yields are normalized to a luminosity of 139 fb$^{-1}$, with the availability of $b$-jet triggers lowering the luminosity to 126 fb$^{-1}$. All selections are cumulative, with the exception of the SR cuts, which are each applied separately.

Cutflow for the low-mass channel for a representative 500 GeV signal. The preselection requires 4 or more jets with $p_\text{T}>25$ GeV and 2 or more $b$-jets with $p_\text{T}>25$ GeV. The $b$-jet cut requires 4 or more $b$-jets with $p_\text{T}>40$ GeV. As the samples are generated with $\mathcal{B}(\tilde{H}\rightarrow h\tilde{G})$=50%, $\mathcal{B}(\tilde{H}\rightarrow Z\tilde{G})$=50% to allow for both decays to be studied, the $hh$ events selection is used to select the events where each of the higgsinos decays to a Higgs boson. Expected yields are normalized to a luminosity of 139 fb$^{-1}$, with the availability of $b$-jet triggers lowering the luminosity to 126 fb$^{-1}$. All selections are cumulative, with the exception of the SR cuts, which are each applied separately.

Cutflow for the low-mass channel for a representative 600 GeV signal. The preselection requires 4 or more jets with $p_\text{T}>25$ GeV and 2 or more $b$-jets with $p_\text{T}>25$ GeV. The $b$-jet cut requires 4 or more $b$-jets with $p_\text{T}>40$ GeV. As the samples are generated with $\mathcal{B}(\tilde{H}\rightarrow h\tilde{G})$=50%, $\mathcal{B}(\tilde{H}\rightarrow Z\tilde{G})$=50% to allow for both decays to be studied, the $hh$ events selection is used to select the events where each of the higgsinos decays to a Higgs boson. Expected yields are normalized to a luminosity of 139 fb$^{-1}$, with the availability of $b$-jet triggers lowering the luminosity to 126 fb$^{-1}$. All selections are cumulative, with the exception of the SR cuts, which are each applied separately.

Cutflow for the low-mass channel for a representative 700 GeV signal. The preselection requires 4 or more jets with $p_\text{T}>25$ GeV and 2 or more $b$-jets with $p_\text{T}>25$ GeV. The $b$-jet cut requires 4 or more $b$-jets with $p_\text{T}>40$ GeV. As the samples are generated with $\mathcal{B}(\tilde{H}\rightarrow h\tilde{G})$=50%, $\mathcal{B}(\tilde{H}\rightarrow Z\tilde{G})$=50% to allow for both decays to be studied, the $hh$ events selection is used to select the events where each of the higgsinos decays to a Higgs boson. Expected yields are normalized to a luminosity of 139 fb$^{-1}$, with the availability of $b$-jet triggers lowering the luminosity to 126 fb$^{-1}$. All selections are cumulative, with the exception of the SR cuts, which are each applied separately.

Cutflow for the low-mass channel for a representative 800 GeV signal. The preselection requires 4 or more jets with $p_\text{T}>25$ GeV and 2 or more $b$-jets with $p_\text{T}>25$ GeV. The $b$-jet cut requires 4 or more $b$-jets with $p_\text{T}>40$ GeV. As the samples are generated with $\mathcal{B}(\tilde{H}\rightarrow h\tilde{G})$=50%, $\mathcal{B}(\tilde{H}\rightarrow Z\tilde{G})$=50% to allow for both decays to be studied, the $hh$ events selection is used to select the events where each of the higgsinos decays to a Higgs boson. Expected yields are normalized to a luminosity of 139 fb$^{-1}$, with the availability of $b$-jet triggers lowering the luminosity to 126 fb$^{-1}$. All selections are cumulative, with the exception of the SR cuts, which are each applied separately.

Cutflow for the low-mass channel for a representative 900 GeV signal. The preselection requires 4 or more jets with $p_\text{T}>25$ GeV and 2 or more $b$-jets with $p_\text{T}>25$ GeV. The $b$-jet cut requires 4 or more $b$-jets with $p_\text{T}>40$ GeV. As the samples are generated with $\mathcal{B}(\tilde{H}\rightarrow h\tilde{G})$=50%, $\mathcal{B}(\tilde{H}\rightarrow Z\tilde{G})$=50% to allow for both decays to be studied, the $hh$ events selection is used to select the events where each of the higgsinos decays to a Higgs boson. Expected yields are normalized to a luminosity of 139 fb$^{-1}$, with the availability of $b$-jet triggers lowering the luminosity to 126 fb$^{-1}$. All selections are cumulative, with the exception of the SR cuts, which are each applied separately.

Cutflow for the low-mass channel for a representative 1000 GeV signal. The preselection requires 4 or more jets with $p_\text{T}>25$ GeV and 2 or more $b$-jets with $p_\text{T}>25$ GeV. The $b$-jet cut requires 4 or more $b$-jets with $p_\text{T}>40$ GeV. As the samples are generated with $\mathcal{B}(\tilde{H}\rightarrow h\tilde{G})$=50%, $\mathcal{B}(\tilde{H}\rightarrow Z\tilde{G})$=50% to allow for both decays to be studied, the $hh$ events selection is used to select the events where each of the higgsinos decays to a Higgs boson. Expected yields are normalized to a luminosity of 139 fb$^{-1}$, with the availability of $b$-jet triggers lowering the luminosity to 126 fb$^{-1}$. All selections are cumulative, with the exception of the SR cuts, which are each applied separately.

Cutflow for the low-mass channel for a representative 1100 GeV signal. The preselection requires 4 or more jets with $p_\text{T}>25$ GeV and 2 or more $b$-jets with $p_\text{T}>25$ GeV. The $b$-jet cut requires 4 or more $b$-jets with $p_\text{T}>40$ GeV. As the samples are generated with $\mathcal{B}(\tilde{H}\rightarrow h\tilde{G})$=50%, $\mathcal{B}(\tilde{H}\rightarrow Z\tilde{G})$=50% to allow for both decays to be studied, the $hh$ events selection is used to select the events where each of the higgsinos decays to a Higgs boson. Expected yields are normalized to a luminosity of 139 fb$^{-1}$, with the availability of $b$-jet triggers lowering the luminosity to 126 fb$^{-1}$. All selections are cumulative, with the exception of the SR cuts, which are each applied separately.

Cutflow for the high-mass channel for a representative 200 GeV signal. The preselection requires 4 or more jets with $p_\text{T}>25$ GeV and 2 or more $b$-jets with $p_\text{T}>25$ GeV. As the samples are generated with $\mathcal{B}(\tilde{H}\rightarrow h\tilde{G})$=50%, $\mathcal{B}(\tilde{H}\rightarrow Z\tilde{G})$=50% to allow for both decays to be studied, the $hh$ events selection is used to select the events where each of the higgsinos decays to a Higgs boson. Expected yields are normalized to a luminosity of 139 fb$^{-1}$. All selections are cumulative, with the exception of the SR cuts, which are each applied separately.

Cutflow for the high-mass channel for a representative 250 GeV signal. The preselection requires 4 or more jets with $p_\text{T}>25$ GeV and 2 or more $b$-jets with $p_\text{T}>25$ GeV. As the samples are generated with $\mathcal{B}(\tilde{H}\rightarrow h\tilde{G})$=50%, $\mathcal{B}(\tilde{H}\rightarrow Z\tilde{G})$=50% to allow for both decays to be studied, the $hh$ events selection is used to select the events where each of the higgsinos decays to a Higgs boson. Expected yields are normalized to a luminosity of 139 fb$^{-1}$. All selections are cumulative, with the exception of the SR cuts, which are each applied separately.

Cutflow for the high-mass channel for a representative 300 GeV signal. The preselection requires 4 or more jets with $p_\text{T}>25$ GeV and 2 or more $b$-jets with $p_\text{T}>25$ GeV. As the samples are generated with $\mathcal{B}(\tilde{H}\rightarrow h\tilde{G})$=50%, $\mathcal{B}(\tilde{H}\rightarrow Z\tilde{G})$=50% to allow for both decays to be studied, the $hh$ events selection is used to select the events where each of the higgsinos decays to a Higgs boson. Expected yields are normalized to a luminosity of 139 fb$^{-1}$. All selections are cumulative, with the exception of the SR cuts, which are each applied separately.

Cutflow for the high-mass channel for a representative 400 GeV signal. The preselection requires 4 or more jets with $p_\text{T}>25$ GeV and 2 or more $b$-jets with $p_\text{T}>25$ GeV. As the samples are generated with $\mathcal{B}(\tilde{H}\rightarrow h\tilde{G})$=50%, $\mathcal{B}(\tilde{H}\rightarrow Z\tilde{G})$=50% to allow for both decays to be studied, the $hh$ events selection is used to select the events where each of the higgsinos decays to a Higgs boson. Expected yields are normalized to a luminosity of 139 fb$^{-1}$. All selections are cumulative, with the exception of the SR cuts, which are each applied separately.

Cutflow for the high-mass channel for a representative 500 GeV signal. The preselection requires 4 or more jets with $p_\text{T}>25$ GeV and 2 or more $b$-jets with $p_\text{T}>25$ GeV. As the samples are generated with $\mathcal{B}(\tilde{H}\rightarrow h\tilde{G})$=50%, $\mathcal{B}(\tilde{H}\rightarrow Z\tilde{G})$=50% to allow for both decays to be studied, the $hh$ events selection is used to select the events where each of the higgsinos decays to a Higgs boson. Expected yields are normalized to a luminosity of 139 fb$^{-1}$. All selections are cumulative, with the exception of the SR cuts, which are each applied separately.

Cutflow for the high-mass channel for a representative 600 GeV signal. The preselection requires 4 or more jets with $p_\text{T}>25$ GeV and 2 or more $b$-jets with $p_\text{T}>25$ GeV. As the samples are generated with $\mathcal{B}(\tilde{H}\rightarrow h\tilde{G})$=50%, $\mathcal{B}(\tilde{H}\rightarrow Z\tilde{G})$=50% to allow for both decays to be studied, the $hh$ events selection is used to select the events where each of the higgsinos decays to a Higgs boson. Expected yields are normalized to a luminosity of 139 fb$^{-1}$. All selections are cumulative, with the exception of the SR cuts, which are each applied separately.

Cutflow for the high-mass channel for a representative 700 GeV signal. The preselection requires 4 or more jets with $p_\text{T}>25$ GeV and 2 or more $b$-jets with $p_\text{T}>25$ GeV. As the samples are generated with $\mathcal{B}(\tilde{H}\rightarrow h\tilde{G})$=50%, $\mathcal{B}(\tilde{H}\rightarrow Z\tilde{G})$=50% to allow for both decays to be studied, the $hh$ events selection is used to select the events where each of the higgsinos decays to a Higgs boson. Expected yields are normalized to a luminosity of 139 fb$^{-1}$. All selections are cumulative, with the exception of the SR cuts, which are each applied separately.

Cutflow for the high-mass channel for a representative 800 GeV signal. The preselection requires 4 or more jets with $p_\text{T}>25$ GeV and 2 or more $b$-jets with $p_\text{T}>25$ GeV. As the samples are generated with $\mathcal{B}(\tilde{H}\rightarrow h\tilde{G})$=50%, $\mathcal{B}(\tilde{H}\rightarrow Z\tilde{G})$=50% to allow for both decays to be studied, the $hh$ events selection is used to select the events where each of the higgsinos decays to a Higgs boson. Expected yields are normalized to a luminosity of 139 fb$^{-1}$. All selections are cumulative, with the exception of the SR cuts, which are each applied separately.

Cutflow for the high-mass channel for a representative 900 GeV signal. The preselection requires 4 or more jets with $p_\text{T}>25$ GeV and 2 or more $b$-jets with $p_\text{T}>25$ GeV. As the samples are generated with $\mathcal{B}(\tilde{H}\rightarrow h\tilde{G})$=50%, $\mathcal{B}(\tilde{H}\rightarrow Z\tilde{G})$=50% to allow for both decays to be studied, the $hh$ events selection is used to select the events where each of the higgsinos decays to a Higgs boson. Expected yields are normalized to a luminosity of 139 fb$^{-1}$. All selections are cumulative, with the exception of the SR cuts, which are each applied separately.

Cutflow for the high-mass channel for a representative 1000 GeV signal. The preselection requires 4 or more jets with $p_\text{T}>25$ GeV and 2 or more $b$-jets with $p_\text{T}>25$ GeV. As the samples are generated with $\mathcal{B}(\tilde{H}\rightarrow h\tilde{G})$=50%, $\mathcal{B}(\tilde{H}\rightarrow Z\tilde{G})$=50% to allow for both decays to be studied, the $hh$ events selection is used to select the events where each of the higgsinos decays to a Higgs boson. Expected yields are normalized to a luminosity of 139 fb$^{-1}$. All selections are cumulative, with the exception of the SR cuts, which are each applied separately.

Cutflow for the high-mass channel for a representative 1100 GeV signal. The preselection requires 4 or more jets with $p_\text{T}>25$ GeV and 2 or more $b$-jets with $p_\text{T}>25$ GeV. As the samples are generated with $\mathcal{B}(\tilde{H}\rightarrow h\tilde{G})$=50%, $\mathcal{B}(\tilde{H}\rightarrow Z\tilde{G})$=50% to allow for both decays to be studied, the $hh$ events selection is used to select the events where each of the higgsinos decays to a Higgs boson. Expected yields are normalized to a luminosity of 139 fb$^{-1}$. All selections are cumulative, with the exception of the SR cuts, which are each applied separately.

Cutflow for the high-mass channel for a representative 1200 GeV signal. The preselection requires 4 or more jets with $p_\text{T}>25$ GeV and 2 or more $b$-jets with $p_\text{T}>25$ GeV. As the samples are generated with $\mathcal{B}(\tilde{H}\rightarrow h\tilde{G})$=50%, $\mathcal{B}(\tilde{H}\rightarrow Z\tilde{G})$=50% to allow for both decays to be studied, the $hh$ events selection is used to select the events where each of the higgsinos decays to a Higgs boson. Expected yields are normalized to a luminosity of 139 fb$^{-1}$. All selections are cumulative, with the exception of the SR cuts, which are each applied separately.

Cutflow for the high-mass channel for a representative 1300 GeV signal. The preselection requires 4 or more jets with $p_\text{T}>25$ GeV and 2 or more $b$-jets with $p_\text{T}>25$ GeV. As the samples are generated with $\mathcal{B}(\tilde{H}\rightarrow h\tilde{G})$=50%, $\mathcal{B}(\tilde{H}\rightarrow Z\tilde{G})$=50% to allow for both decays to be studied, the $hh$ events selection is used to select the events where each of the higgsinos decays to a Higgs boson. Expected yields are normalized to a luminosity of 139 fb$^{-1}$. All selections are cumulative, with the exception of the SR cuts, which are each applied separately.

Cutflow for the high-mass channel for a representative 1400 GeV signal. The preselection requires 4 or more jets with $p_\text{T}>25$ GeV and 2 or more $b$-jets with $p_\text{T}>25$ GeV. As the samples are generated with $\mathcal{B}(\tilde{H}\rightarrow h\tilde{G})$=50%, $\mathcal{B}(\tilde{H}\rightarrow Z\tilde{G})$=50% to allow for both decays to be studied, the $hh$ events selection is used to select the events where each of the higgsinos decays to a Higgs boson. Expected yields are normalized to a luminosity of 139 fb$^{-1}$. All selections are cumulative, with the exception of the SR cuts, which are each applied separately.

Cutflow for the high-mass channel for a representative 1500 GeV signal. The preselection requires 4 or more jets with $p_\text{T}>25$ GeV and 2 or more $b$-jets with $p_\text{T}>25$ GeV. As the samples are generated with $\mathcal{B}(\tilde{H}\rightarrow h\tilde{G})$=50%, $\mathcal{B}(\tilde{H}\rightarrow Z\tilde{G})$=50% to allow for both decays to be studied, the $hh$ events selection is used to select the events where each of the higgsinos decays to a Higgs boson. Expected yields are normalized to a luminosity of 139 fb$^{-1}$. All selections are cumulative, with the exception of the SR cuts, which are each applied separately.

Cutflow for two representative signal samples used in this analysis. The excited tau mass $m_{\tau^\ast} = 2.75$ TeV and the contact interaction scale $\Lambda=10$ TeV. The LQ mass $m_\mathrm{LQ} = 1.3$ TeV. The event yields include all correction factors applied to simulation, and is normalised to 139 fb$^{-1}$.

Acceptance x efficiency of the $\tau^\ast$ signal SR selection

Acceptance x efficiency of the LQ signal SR selection

Background-only fit to the observed numbers of events (black points) as a function of the diphoton invariant mass $m_{\gamma\gamma}$.

Expected and observed 95% CL upper limits on the signal fiducial cross section $\sigma_{\textrm{fid}}$, assuming 100% branching ratio for ALP decay into two photons, as a function of the hypothetical ALP mass $m_{\textrm{X}}$. The 1$\sigma$ and 2$\sigma$ confidence intervals are shown by the coloured bands.

Expected and observed 95% CL upper limits on the the ALP coupling constant, assuming 100% branching ratio for ALP decay into two photons, as a function of the hypothetical ALP mass $m_{\textrm{X}}$. The 1$\sigma$ and 2$\sigma$ confidence intervals are shown by the coloured bands. Contours of the ALP natural width $\Gamma$ are illustrated by the smooth blue solid lines.

Fiducial region acceptance of signal events for coupling constant f$^{-1}=0.05$ TeV$^{-1}$. EL, SD, and DD stand for the exclusive, single-dissociative, and double-dissociative processes. The acceptances are parameterised as $p_0 + p^{i}_{1} e^{p_{2}^{i} m_{\textrm{X}}} + p_{3}^{i} e^{p_{4}^{i} m_{\textrm{X}}} $ , where $i$ specifies the process type, exclusive, single-dissociative, or double-dissociative.

Correction factors $A_{\textrm{X}}/A_{0}$ and $C_{\textrm{X}}$ for the exclusive signal events with ALP coupling constant f$^{-1}$=0.05 TeV$^{-1}$, where acceptance $A_{\textrm{X}}$ and efficiency $C_{\textrm{X}}$ are defined in arXiv:2102.13405. $A_{\textrm{X}}$/$A_{0}$ and $C_{\textrm{X}}$ are parameterised as $p_0 + p^{i}_{1} e^{p_{2}^{i} m_{\textrm{X}}} + p_{3}^{i} e^{p_{4}^{i} m_{\textrm{X}}}$, where $i$ specifies the process type.

Simulated ALP signal cutflow for $m_{\rm{X}}$=300 GeV and $m_{\rm{X}}$=1200 GeV. The yields are normalized to a luminosity of 14.6 fb$^{-1}$. Selections are applied for each side after the selection on acoplanarity. The union of the sets of selected events is taken finally.

Expected exclusion limits at 95% CL from Fig 7(a) for $\tilde{g}$ decays into SM gauge bosons and $\tilde{\chi}^{0}_{1}$

One $\sigma$ band of expected exclusion limits at 95% CL from Fig 7(a) for $\tilde{g}$ decays into SM gauge bosons and $\tilde{\chi}^{0}_{1}$

Observed exclusion limits at 95% CL from Fig 7(c) for $\tilde{g}$ decays into sleptons and subsequently to SM leptons and $\tilde{\chi}^{0}_{1}$

Positive one $\sigma$ observed exclusion limits at 95% CL from Fig 7(c) for $\tilde{g}$ decays into sleptons and subsequently to SM leptons and $\tilde{\chi}^{0}_{1}$

Negative one $\sigma$ observed exclusion limits at 95% CL from Fig 7(c) for $\tilde{g}$ decays into sleptons and subsequently to SM leptons and $\tilde{\chi}^{0}_{1}$

Expected exclusion limits at 95% CL from Fig 7(c) for $\tilde{g}$ decays into sleptons and subsequently to SM leptons and $\tilde{\chi}^{0}_{1}$

One $\sigma$ band of expected exclusion limits at 95% CL from Fig 7(c) for $\tilde{g}$ decays into sleptons and subsequently to SM leptons and $\tilde{\chi}^{0}_{1}$

Observed exclusion limits at 95% CL from Fig 7(f) for $\tilde{g}$ decays into anti-top and $\tilde{t}$ and $\tilde{t}$ decays via a non-zero RPV coupling $\lambda''$

Positive one $\sigma$ observed exclusion limits at 95% CL from Fig 7(f) for $\tilde{g}$ decays into anti-top and $\tilde{t}$ and $\tilde{t}$ decays via a non-zero RPV coupling $\lambda''$

Negative one $\sigma$ observed exclusion limits at 95% CL from Fig 7(f) for $\tilde{g}$ decays into anti-top and $\tilde{t}$ and $\tilde{t}$ decays via a non-zero RPV coupling $\lambda''$

Expected exclusion limits at 95% CL from Fig 7(f) for $\tilde{g}$ decays into anti-top and $\tilde{t}$ and $\tilde{t}$ decays via a non-zero RPV coupling $\lambda''$

One $\sigma$ band of expected exclusion limits at 95% CL from Fig 7(f) for $\tilde{g}$ decays into anti-top and $\tilde{t}$ and $\tilde{t}$ decays via a non-zero RPV coupling $\lambda''$

Observed exclusion limits at 95% CL from Fig 7(e) for direct $\tilde{\chi_{1}^{0}}$ decay into SM leptons and quarks via a non-zero RPV coupling $\lambda'$

Positive one $\sigma$ observed exclusion limits at 95% CL from Fig 7(e) for direct $\tilde{\chi_{1}^{0}}$ decay into SM leptons and quarks via a non-zero RPV coupling $\lambda'$

Negative one $\sigma$ observed exclusion limits at 95% CL from Fig 7(e) for direct $\tilde{\chi_{1}^{0}}$ decay into SM leptons and quarks via a non-zero RPV coupling $\lambda'$

Expected exclusion limits at 95% CL from Fig 7(e) for direct $\tilde{\chi_{1}^{0}}$ decay into SM leptons and quarks via a non-zero RPV coupling $\lambda'$

One $\sigma$ band of expected exclusion limits at 95% CL from Fig 7(e) for direct $\tilde{\chi_{1}^{0}}$ decay into SM leptons and quarks via a non-zero RPV coupling $\lambda'$

Observed exclusion limits at 95% CL from Fig 7(b) for $\tilde{q}$ decays into SM gauge bosons and $\tilde{\chi}^{0}_{1}$

Positive one $\sigma$ observed exclusion limits at 95% CL from Fig 7(b) for $\tilde{q}$ decays into SM gauge bosons and $\tilde{\chi}^{0}_{1}$

Negative one $\sigma$ observed exclusion limits at 95% CL from Fig 7(b) for $\tilde{q}$ decays into SM gauge bosons and $\tilde{\chi}^{0}_{1}$

Expected exclusion limits at 95% CL from Fig 7(b) for $\tilde{q}$ decays into SM gauge bosons and $\tilde{\chi}^{0}_{1}$

One $\sigma$ band of expected exclusion limits at 95% CL from Fig 7(b) for $\tilde{q}$ decays into SM gauge bosons and $\tilde{\chi}^{0}_{1}$

Observed exclusion limits at 95% CL from Fig 7(d) for $\tilde{q}$ decays into sleptons and subsequently to SM leptons and $\tilde{\chi}^{0}_{1}$

Positive one $\sigma$ observed exclusion limits at 95% CL from Fig 7(d) for $\tilde{q}$ decays into sleptons and subsequently to SM leptons and $\tilde{\chi}^{0}_{1}$

Negative one $\sigma$ observed exclusion limits at 95% CL from Fig 7(d) for $\tilde{q}$ decays into sleptons and subsequently to SM leptons and $\tilde{\chi}^{0}_{1}$

Expected exclusion limits at 95% CL from Fig 7(d) for $\tilde{q}$ decays into sleptons and subsequently to SM leptons and $\tilde{\chi}^{0}_{1}$

One $\sigma$ band of expected exclusion limits at 95% CL from Fig 7(d) for $\tilde{q}$ decays into sleptons and subsequently to SM leptons and $\tilde{\chi}^{0}_{1}$

N-1 distribution for $m_{\mathrm{eff}}$of observed data and expected background in SRGGWZ-H.

N-1 distribution for $E_{\mathrm{T}}^{\mathrm{miss}}$of observed data and expected background in SRGGSlep-M.

N-1 distribution for $\sum{p_{\mathrm{T}}^\mathrm{jet}}$of observed data and expected background in SRUDD-ge2b.

N-1 distribution for $m_{\mathrm{eff}}$of observed data and expected background in SRLQD.

N-1 distribution for $m_{\mathrm{eff}}$of observed data and expected background in SRSSWZ-H.

N-1 distribution for $m_{\mathrm{eff}}$of observed data and expected background in SRSSSlep-H(loose).

Signal acceptance for SRGGWZ-H signal region from Fig 10(c) in a SUSY scenario where $\tilde{g}$ decays into SM gauge bosons and $\tilde{\chi}^{0}_{1}$

Signal efficiency for SRGGWZ-H signal region from Fig 15(c) in a SUSY scenario where $\tilde{g}$ decays into SM gauge bosons and $\tilde{\chi}^{0}_{1}$

Signal acceptance for SRGGWZ-M signal region from Fig 10(b) in a SUSY scenario where $\tilde{g}$ decays into SM gauge bosons and $\tilde{\chi}^{0}_{1}$

Signal efficiency for SRGGWZ-M signal region from Fig 15(b) in a SUSY scenario where $\tilde{g}$ decays into SM gauge bosons and $\tilde{\chi}^{0}_{1}$

Signal acceptance for SRGGWZ-L signal region from Fig 10(a) in a SUSY scenario where $\tilde{g}$ decays into SM gauge bosons and $\tilde{\chi}^{0}_{1}$

Signal efficiency for SRGGWZ-L signal region from Fig 15(a) in a SUSY scenario where $\tilde{g}$ decays into SM gauge bosons and $\tilde{\chi}^{0}_{1}$

Signal acceptance for SRGGSlep-L signal region from Fig 12(a) in a SUSY scenario where $\tilde{g}$ decays into sleptons and subsequently to SM leptons and $\tilde{\chi}^{0}_{1}$

Signal efficiency for SRGGSlep-L signal region from Fig 17(a) in a SUSY scenario where $\tilde{g}$ decays into sleptons and subsequently to SM leptons and $\tilde{\chi}^{0}_{1}$

Signal acceptance for SRGGSlep-M signal region from Fig 12(b) in a SUSY scenario where $\tilde{g}$ decays into sleptons and subsequently to SM leptons and $\tilde{\chi}^{0}_{1}$

Signal efficiency for SRGGSlep-M signal region from Fig 17(b) in a SUSY scenario where $\tilde{g}$ decays into sleptons and subsequently to SM leptons and $\tilde{\chi}^{0}_{1}$

Signal acceptance for SRGGSlep-H signal region from Fig 12(c) in a SUSY scenario where $\tilde{g}$ decays into sleptons and subsequently to SM leptons and $\tilde{\chi}^{0}_{1}$

Signal efficiency for SRGGSlep-H signal region from Fig 17(c) in a SUSY scenario where $\tilde{g}$ decays into sleptons and subsequently to SM leptons and $\tilde{\chi}^{0}_{1}$

Signal acceptance for SRUDD-1b signal region from Fig 14(b) in a SUSY scenario where $\tilde{g}$ decays into anti-top and $\tilde{t}$ and $\tilde{t}$ decays via a non-zero RPV coupling $\lambda''$

Signal efficiency for SRUDD-1b signal region from Fig 19(b) in a SUSY scenario where $\tilde{g}$ decays into anti-top and $\tilde{t}$ and $\tilde{t}$ decays via a non-zero RPV coupling $\lambda''$

Signal acceptance for SRUDD-2b signal region from Fig 14(c) in a SUSY scenario where $\tilde{g}$ decays into anti-top and $\tilde{t}$ and $\tilde{t}$ decays via a non-zero RPV coupling $\lambda''$

Signal efficiency for SRUDD-2b signal region from Fig 19(c) in a SUSY scenario where $\tilde{g}$ decays into anti-top and $\tilde{t}$ and $\tilde{t}$ decays via a non-zero RPV coupling $\lambda''$

Signal acceptance for SRUDD-ge2b signal region from Fig 14(d) in a SUSY scenario where $\tilde{g}$ decays into anti-top and $\tilde{t}$ and $\tilde{t}$ decays via a non-zero RPV coupling $\lambda''$

Signal efficiency for SRUDD-ge2b signal region from Fig 19(d) in a SUSY scenario where $\tilde{g}$ decays into anti-top and $\tilde{t}$ and $\tilde{t}$ decays via a non-zero RPV coupling $\lambda''$

Signal acceptance for SRUDD-ge3b signal region from Fig 14(e) in a SUSY scenario where $\tilde{g}$ decays into anti-top and $\tilde{t}$ and $\tilde{t}$ decays via a non-zero RPV coupling $\lambda''$

Signal efficiency for SRUDD-ge3b signal region from Fig 19(e) in a SUSY scenario where $\tilde{g}$ decays into anti-top and $\tilde{t}$ and $\tilde{t}$ decays via a non-zero RPV coupling $\lambda''$

Signal acceptance for SRLQD signal region from Fig 14(a) in a SUSY scenario where direct $\tilde{\chi_{1}^{0}}$ decay into SM leptons and quarks via a non-zero RPV coupling $\lambda'$

Signal efficiency for SRLQD signal region from Fig 19(a) in a SUSY scenario where direct $\tilde{\chi_{1}^{0}}$ decay into SM leptons and quarks via a non-zero RPV coupling $\lambda'$

Signal acceptance for SRSSWZ-L signal region from Fig 11(a) in a SUSY scenario where $\tilde{q}$ decays into SM gauge bosons and $\tilde{\chi}^{0}_{1}$

Signal efficiency for SRSSWZ-L signal region from Fig 16(a) in a SUSY scenario where $\tilde{q}$ decays into SM gauge bosons and $\tilde{\chi}^{0}_{1}$

Signal acceptance for SRSSWZ-ML signal region from Fig 11(b) in a SUSY scenario where $\tilde{q}$ decays into SM gauge bosons and $\tilde{\chi}^{0}_{1}$

Signal efficiency for SRSSWZ-ML signal region from Fig 16(b) in a SUSY scenario where $\tilde{q}$ decays into SM gauge bosons and $\tilde{\chi}^{0}_{1}$

Signal acceptance for SRSSWZ-MH signal region from Fig 11(c) in a SUSY scenario where $\tilde{q}$ decays into SM gauge bosons and $\tilde{\chi}^{0}_{1}$

Signal efficiency for SRSSWZ-MH signal region from Fig 16(c) in a SUSY scenario where $\tilde{q}$ decays into SM gauge bosons and $\tilde{\chi}^{0}_{1}$

Signal acceptance for SRSSWZ-H signal region from Fig 11(d) in a SUSY scenario where $\tilde{q}$ decays into SM gauge bosons and $\tilde{\chi}^{0}_{1}$

Signal efficiency for SRSSWZ-H signal region from Fig 16(d) in a SUSY scenario where $\tilde{q}$ decays into SM gauge bosons and $\tilde{\chi}^{0}_{1}$

Signal acceptance for SRSSSlep-H signal region from Fig 13(d) in a SUSY scenario where $\tilde{q}$ decays into sleptons and subsequently to SM leptons and $\tilde{\chi}^{0}_{1}$

Signal efficiency for SRSSSlep-H signal region from Fig 18(d) in a SUSY scenario where $\tilde{q}$ decays into sleptons and subsequently to SM leptons and $\tilde{\chi}^{0}_{1}$

Signal acceptance for SRSSSlep-MH signal region from Fig 13(c) in a SUSY scenario where $\tilde{q}$ decays into sleptons and subsequently to SM leptons and $\tilde{\chi}^{0}_{1}$

Signal efficiency for SRSSSlep-MH signal region from Fig 18(c) in a SUSY scenario where $\tilde{q}$ decays into sleptons and subsequently to SM leptons and $\tilde{\chi}^{0}_{1}$

Signal acceptance for SRSSSlep-L signal region from Fig 13(a) in a SUSY scenario where $\tilde{q}$ decays into sleptons and subsequently to SM leptons and $\tilde{\chi}^{0}_{1}$

Signal efficiency for SRSSSlep-L signal region from Fig 18(a) in a SUSY scenario where $\tilde{q}$ decays into sleptons and subsequently to SM leptons and $\tilde{\chi}^{0}_{1}$

Signal acceptance for SRSSSlep-ML signal region from Fig 13(b) in a SUSY scenario where $\tilde{q}$ decays into sleptons and subsequently to SM leptons and $\tilde{\chi}^{0}_{1}$

Signal efficiency for SRSSSlep-ML signal region from Fig 18(b) in a SUSY scenario where $\tilde{q}$ decays into sleptons and subsequently to SM leptons and $\tilde{\chi}^{0}_{1}$

Signal acceptance for SRSSSlep-H(loose) signal region from Fig 13(e) in a SUSY scenario where $\tilde{q}$ decays into sleptons and subsequently to SM leptons and $\tilde{\chi}^{0}_{1}$

Signal efficiency for SRSSSlep-H(loose) signal region from Fig 18(e) in a SUSY scenario where $\tilde{q}$ decays into sleptons and subsequently to SM leptons and $\tilde{\chi}^{0}_{1}$

Number of signal events expected for 139 fb$^{-1}$ at different stages of the event selection for the signal region SRGGWZ-H in a susy scenario where $\tilde{g}$ decays into SM gauge bosons and $\tilde{\chi}^{0}_{1}$. The masses of the superpartners involved in the process are set to $m(\tilde{g})$ = 1400 GeV, $m(\tilde{\chi_{1}^{0}})$ = 1000 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 SRGGWZ-M in a susy scenario where $\tilde{g}$ decays into SM gauge bosons and $\tilde{\chi}^{0}_{1}$. The masses of the superpartners involved in the process are set to $m(\tilde{g})$ = 1400 GeV, $m(\tilde{\chi_{1}^{0}})$ = 1000 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 SRGGWZ-L in a susy scenario where $\tilde{g}$ decays into SM gauge bosons and $\tilde{\chi}^{0}_{1}$. The masses of the superpartners involved in the process are set to $m(\tilde{g})$ = 1400 GeV, $m(\tilde{\chi_{1}^{0}})$ = 1000 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 SRGGSlep-L in a susy scenario where $\tilde{g}$ decays into sleptons and subsequently to SM leptons and $\tilde{\chi}^{0}_{1}$. The masses of the superpartners involved in the process are set to $m(\tilde{g})$ = 2000 GeV, $m(\tilde{\chi_{1}^{0}})$ = 500 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 SRGGSlep-M in a susy scenario where $\tilde{g}$ decays into sleptons and subsequently to SM leptons and $\tilde{\chi}^{0}_{1}$. The masses of the superpartners involved in the process are set to $m(\tilde{g})$ = 2000 GeV, $m(\tilde{\chi_{1}^{0}})$ = 500 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 SRGGSlep-H in a susy scenario where $\tilde{g}$ decays into sleptons and subsequently to SM leptons and $\tilde{\chi}^{0}_{1}$. The masses of the superpartners involved in the process are set to $m(\tilde{g})$ = 2000 GeV, $m(\tilde{\chi_{1}^{0}})$ = 500 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 SRUDD-1b in a susy scenario where $\tilde{g}$ decays into anti-top and $\tilde{t}$ and $\tilde{t}$ decays via a non-zero RPV coupling $\lambda''$. The masses of the superpartners involved in the process are set to $m(\tilde{g})$ = 1600 GeV, $m(\tilde{t})$ = 600 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 SRUDD-2b in a susy scenario where $\tilde{g}$ decays into anti-top and $\tilde{t}$ and $\tilde{t}$ decays via a non-zero RPV coupling $\lambda''$. The masses of the superpartners involved in the process are set to $m(\tilde{g})$ = 1600 GeV, $m(\tilde{t})$ = 600 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 SRUDD-ge2b in a susy scenario where $\tilde{g}$ decays into anti-top and $\tilde{t}$ and $\tilde{t}$ decays via a non-zero RPV coupling $\lambda''$. The masses of the superpartners involved in the process are set to $m(\tilde{g})$ = 1600 GeV, $m(\tilde{t})$ = 600 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 SRUDD-ge3b in a susy scenario where $\tilde{g}$ decays into anti-top and $\tilde{t}$ and $\tilde{t}$ decays via a non-zero RPV coupling $\lambda''$. The masses of the superpartners involved in the process are set to $m(\tilde{g})$ = 1600 GeV, $m(\tilde{t})$ = 600 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 SRLQD in a susy scenario where direct $\tilde{\chi_{1}^{0}}$ decay into SM leptons and quarks via a non-zero RPV coupling $\lambda'$. The masses of the superpartners involved in the process are set to $m(\tilde{g})$ = 2200 GeV, $m(\tilde{\chi_{1}^{0}})$ = 1870 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 SRSSWZ-L in a susy scenario where $\tilde{q}$ decays into SM gauge bosons and $\tilde{\chi}^{0}_{1}$. The masses of the superpartners involved in the process are set to $m(\tilde{q})$ = 800 GeV, $m(\tilde{\chi_{1}^{0}})$ = 600 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 SRSSWZ-ML in a susy scenario where $\tilde{q}$ decays into SM gauge bosons and $\tilde{\chi}^{0}_{1}$. The masses of the superpartners involved in the process are set to $m(\tilde{q})$ = 800 GeV, $m(\tilde{\chi_{1}^{0}})$ = 600 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 SRSSWZ-MH in a susy scenario where $\tilde{q}$ decays into SM gauge bosons and $\tilde{\chi}^{0}_{1}$. The masses of the superpartners involved in the process are set to $m(\tilde{q})$ = 800 GeV, $m(\tilde{\chi_{1}^{0}})$ = 600 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 SRSSWZ-H in a susy scenario where $\tilde{q}$ decays into SM gauge bosons and $\tilde{\chi}^{0}_{1}$. The masses of the superpartners involved in the process are set to $m(\tilde{q})$ = 800 GeV, $m(\tilde{\chi_{1}^{0}})$ = 600 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 SRSSSlep-H in a susy scenario where $\tilde{q}$ decays into sleptons and subsequently to SM leptons and $\tilde{\chi}^{0}_{1}$. The masses of the superpartners involved in the process are set to $m(\tilde{q})$ = 1000 GeV, $m(\tilde{\chi_{1}^{0}})$ = 800 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 SRSSSlep-MH in a susy scenario where $\tilde{q}$ decays into sleptons and subsequently to SM leptons and $\tilde{\chi}^{0}_{1}$. The masses of the superpartners involved in the process are set to $m(\tilde{q})$ = 1000 GeV, $m(\tilde{\chi_{1}^{0}})$ = 800 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 SRSSSlep-L in a susy scenario where $\tilde{q}$ decays into sleptons and subsequently to SM leptons and $\tilde{\chi}^{0}_{1}$. The masses of the superpartners involved in the process are set to $m(\tilde{q})$ = 1000 GeV, $m(\tilde{\chi_{1}^{0}})$ = 800 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 SRSSSlep-ML in a susy scenario where $\tilde{q}$ decays into sleptons and subsequently to SM leptons and $\tilde{\chi}^{0}_{1}$. The masses of the superpartners involved in the process are set to $m(\tilde{q})$ = 1000 GeV, $m(\tilde{\chi_{1}^{0}})$ = 800 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 SRSSSlep-H(loose) in a susy scenario where $\tilde{q}$ decays into sleptons and subsequently to SM leptons and $\tilde{\chi}^{0}_{1}$. The masses of the superpartners involved in the process are set to $m(\tilde{q})$ = 1000 GeV, $m(\tilde{\chi_{1}^{0}})$ = 800 GeV. Only statistical uncertainties are shown.

Cross-section upper limits at 95% CL from Fig1(a) for $\tilde{g}$ decays into SM gauge bosons and $\tilde{\chi}^{0}_{1}$

Cross-section upper limits at 95% CL from Fig1(c) for $\tilde{g}$ decays into sleptons and subsequently to SM leptons and $\tilde{\chi}^{0}_{1}$

Cross-section upper limits at 95% CL from Fig1(f) for $\tilde{g}$ decays into anti-top and $\tilde{t}$ and $\tilde{t}$ decays via a non-zero RPV coupling $\lambda''$

Cross-section upper limits at 95% CL from Fig1(e) for direct $\tilde{\chi_{1}^{0}}$ decay into SM leptons and quarks via a non-zero RPV coupling $\lambda'$

Cross-section upper limits at 95% CL from Fig1(b) for $\tilde{q}$ decays into SM gauge bosons and $\tilde{\chi}^{0}_{1}$

Cross-section upper limits at 95% CL from Fig1(d) for $\tilde{q}$ decays into sleptons and subsequently to SM leptons and $\tilde{\chi}^{0}_{1}$

Observed 95% CL upper limits on the cross section for all masses and charges of photon-fusion pair-produced spin-1/2 monopoles as a function of mass for magnetic charges $|g|=1g_D$ and $|g|=2g_D$.

Observed 95% CL upper limits on the cross section for all masses and charges of Drell-Yan spin-0 HECOs production as a function of mass for various values of electric charge in the range $20\le|z|\le100$.

Observed 95% CL upper limits on the cross section for all masses and charges of Drell-Yan spin-1/2 HECOs production as a function of mass for various values of electric charge in the range $20\le|z|\le100$.

Observed 95% CL upper limits on the cross section for all masses and charges of photon-fusion pair-produced spin-0 HECOs as a function of mass for various values of electric charge in the range $20\le|z|\le100$.

Observed 95% CL upper limits on the cross section for all masses and charges of photon-fusion pair-produced spin-1/2 HECOs as a function of mass for various values of electric charge in the range $20\le|z|\le100$.

Selection efficiency as a function of transverse kinetic energy $E^\text{kin}_\text{T}=E_\text{kin}\sin\theta$ and pseudorapidity $|\eta|$ for $g=1g_\textrm{D}$ monopoles of mass 200 GeV.

Selection efficiency as a function of transverse kinetic energy $E^\text{kin}_\text{T}=E_\text{kin}\sin\theta$ and pseudorapidity $|\eta|$ for $g=1g_\textrm{D}$ monopoles of mass 500 GeV.

Selection efficiency as a function of transverse kinetic energy $E^\text{kin}_\text{T}=E_\text{kin}\sin\theta$ and pseudorapidity $|\eta|$ for $g=1g_\textrm{D}$ monopoles of mass 1000 GeV.

Selection efficiency as a function of transverse kinetic energy $E^\text{kin}_\text{T}=E_\text{kin}\sin\theta$ and pseudorapidity $|\eta|$ for $g=1g_\textrm{D}$ monopoles of mass 1500 GeV.

Selection efficiency as a function of transverse kinetic energy $E^\text{kin}_\text{T}=E_\text{kin}\sin\theta$ and pseudorapidity $|\eta|$ for $g=1g_\textrm{D}$ monopoles of mass 2000 GeV.

Selection efficiency as a function of transverse kinetic energy $E^\text{kin}_\text{T}=E_\text{kin}\sin\theta$ and pseudorapidity $|\eta|$ for $g=1g_\textrm{D}$ monopoles of mass 2500 GeV.

Selection efficiency as a function of transverse kinetic energy $E^\text{kin}_\text{T}=E_\text{kin}\sin\theta$ and pseudorapidity $|\eta|$ for $g=1g_\textrm{D}$ monopoles of mass 3000 GeV.

Selection efficiency as a function of transverse kinetic energy $E^\text{kin}_\text{T}=E_\text{kin}\sin\theta$ and pseudorapidity $|\eta|$ for $g=1g_\textrm{D}$ monopoles of mass 4000 GeV.

Selection efficiency as a function of transverse kinetic energy $E^\text{kin}_\text{T}=E_\text{kin}\sin\theta$ and pseudorapidity $|\eta|$ for $g=2g_\textrm{D}$ monopoles of mass 200 GeV.

Selection efficiency as a function of transverse kinetic energy $E^\text{kin}_\text{T}=E_\text{kin}\sin\theta$ and pseudorapidity $|\eta|$ for $g=2g_\textrm{D}$ monopoles of mass 500 GeV.

Selection efficiency as a function of transverse kinetic energy $E^\text{kin}_\text{T}=E_\text{kin}\sin\theta$ and pseudorapidity $|\eta|$ for $g=2g_\textrm{D}$ monopoles of mass 1000 GeV.

Selection efficiency as a function of transverse kinetic energy $E^\text{kin}_\text{T}=E_\text{kin}\sin\theta$ and pseudorapidity $|\eta|$ for $g=2g_\textrm{D}$ monopoles of mass 1500 GeV.

Selection efficiency as a function of transverse kinetic energy $E^\text{kin}_\text{T}=E_\text{kin}\sin\theta$ and pseudorapidity $|\eta|$ for $g=2g_\textrm{D}$ monopoles of mass 2000 GeV.

Selection efficiency as a function of transverse kinetic energy $E^\text{kin}_\text{T}=E_\text{kin}\sin\theta$ and pseudorapidity $|\eta|$ for $g=2g_\textrm{D}$ monopoles of mass 2500 GeV.

Selection efficiency as a function of transverse kinetic energy $E^\text{kin}_\text{T}=E_\text{kin}\sin\theta$ and pseudorapidity $|\eta|$ for $g=2g_\textrm{D}$ monopoles of mass 3000 GeV.

Selection efficiency as a function of transverse kinetic energy $E^\text{kin}_\text{T}=E_\text{kin}\sin\theta$ and pseudorapidity $|\eta|$ for $g=2g_\textrm{D}$ monopoles of mass 4000 GeV.

Selection efficiency as a function of transverse kinetic energy $E^\text{kin}_\text{T}=E_\text{kin}\sin\theta$ and pseudorapidity $|\eta|$ for HECOs of charge $|z|=20$ of mass 200 GeV.

Selection efficiency as a function of transverse kinetic energy $E^\text{kin}_\text{T}=E_\text{kin}\sin\theta$ and pseudorapidity $|\eta|$ for HECOs of charge $|z|=20$ of mass 500 GeV.

Selection efficiency as a function of transverse kinetic energy $E^\text{kin}_\text{T}=E_\text{kin}\sin\theta$ and pseudorapidity $|\eta|$ for HECOs of charge $|z|=20$ of mass 1000 GeV.

Selection efficiency as a function of transverse kinetic energy $E^\text{kin}_\text{T}=E_\text{kin}\sin\theta$ and pseudorapidity $|\eta|$ for HECOs of charge $|z|=20$ of mass 1500 GeV.

Selection efficiency as a function of transverse kinetic energy $E^\text{kin}_\text{T}=E_\text{kin}\sin\theta$ and pseudorapidity $|\eta|$ for HECOs of charge $|z|=20$ of mass 2000 GeV.

Selection efficiency as a function of transverse kinetic energy $E^\text{kin}_\text{T}=E_\text{kin}\sin\theta$ and pseudorapidity $|\eta|$ for HECOs of charge $|z|=20$ of mass 2500 GeV.

Selection efficiency as a function of transverse kinetic energy $E^\text{kin}_\text{T}=E_\text{kin}\sin\theta$ and pseudorapidity $|\eta|$ for HECOs of charge $|z|=20$ of mass 3000 GeV.

Selection efficiency as a function of transverse kinetic energy $E^\text{kin}_\text{T}=E_\text{kin}\sin\theta$ and pseudorapidity $|\eta|$ for HECOs of charge $|z|=20$ of mass 4000 GeV.

Selection efficiency as a function of transverse kinetic energy $E^\text{kin}_\text{T}=E_\text{kin}\sin\theta$ and pseudorapidity $|\eta|$ for HECOs of charge $|z|=40$ of mass 200 GeV.

Selection efficiency as a function of transverse kinetic energy $E^\text{kin}_\text{T}=E_\text{kin}\sin\theta$ and pseudorapidity $|\eta|$ for HECOs of charge $|z|=40$ of mass 500 GeV.

Selection efficiency as a function of transverse kinetic energy $E^\text{kin}_\text{T}=E_\text{kin}\sin\theta$ and pseudorapidity $|\eta|$ for HECOs of charge $|z|=40$ of mass 1000 GeV.

Selection efficiency as a function of transverse kinetic energy $E^\text{kin}_\text{T}=E_\text{kin}\sin\theta$ and pseudorapidity $|\eta|$ for HECOs of charge $|z|=40$ of mass 1500 GeV.

Selection efficiency as a function of transverse kinetic energy $E^\text{kin}_\text{T}=E_\text{kin}\sin\theta$ and pseudorapidity $|\eta|$ for HECOs of charge $|z|=40$ of mass 2000 GeV.

Selection efficiency as a function of transverse kinetic energy $E^\text{kin}_\text{T}=E_\text{kin}\sin\theta$ and pseudorapidity $|\eta|$ for HECOs of charge $|z|=40$ of mass 2500 GeV.

Selection efficiency as a function of transverse kinetic energy $E^\text{kin}_\text{T}=E_\text{kin}\sin\theta$ and pseudorapidity $|\eta|$ for HECOs of charge $|z|=40$ of mass 3000 GeV.

Selection efficiency as a function of transverse kinetic energy $E^\text{kin}_\text{T}=E_\text{kin}\sin\theta$ and pseudorapidity $|\eta|$ for HECOs of charge $|z|=40$ of mass 4000 GeV.

Selection efficiency as a function of transverse kinetic energy $E^\text{kin}_\text{T}=E_\text{kin}\sin\theta$ and pseudorapidity $|\eta|$ for HECOs of charge $|z|=60$ of mass 200 GeV.

Selection efficiency as a function of transverse kinetic energy $E^\text{kin}_\text{T}=E_\text{kin}\sin\theta$ and pseudorapidity $|\eta|$ for HECOs of charge $|z|=60$ of mass 500 GeV.

Selection efficiency as a function of transverse kinetic energy $E^\text{kin}_\text{T}=E_\text{kin}\sin\theta$ and pseudorapidity $|\eta|$ for HECOs of charge $|z|=60$ of mass 1000 GeV.

Selection efficiency as a function of transverse kinetic energy $E^\text{kin}_\text{T}=E_\text{kin}\sin\theta$ and pseudorapidity $|\eta|$ for HECOs of charge $|z|=60$ of mass 1500 GeV.

Selection efficiency as a function of transverse kinetic energy $E^\text{kin}_\text{T}=E_\text{kin}\sin\theta$ and pseudorapidity $|\eta|$ for HECOs of charge $|z|=60$ of mass 2000 GeV.

Selection efficiency as a function of transverse kinetic energy $E^\text{kin}_\text{T}=E_\text{kin}\sin\theta$ and pseudorapidity $|\eta|$ for HECOs of charge $|z|=60$ of mass 2500 GeV.

Selection efficiency as a function of transverse kinetic energy $E^\text{kin}_\text{T}=E_\text{kin}\sin\theta$ and pseudorapidity $|\eta|$ for HECOs of charge $|z|=60$ of mass 3000 GeV.

Selection efficiency as a function of transverse kinetic energy $E^\text{kin}_\text{T}=E_\text{kin}\sin\theta$ and pseudorapidity $|\eta|$ for HECOs of charge $|z|=60$ of mass 4000 GeV.

Selection efficiency as a function of transverse kinetic energy $E^\text{kin}_\text{T}=E_\text{kin}\sin\theta$ and pseudorapidity $|\eta|$ for HECOs of charge $|z|=80$ of mass 200 GeV.

Selection efficiency as a function of transverse kinetic energy $E^\text{kin}_\text{T}=E_\text{kin}\sin\theta$ and pseudorapidity $|\eta|$ for HECOs of charge $|z|=80$ of mass 500 GeV.

Selection efficiency as a function of transverse kinetic energy $E^\text{kin}_\text{T}=E_\text{kin}\sin\theta$ and pseudorapidity $|\eta|$ for HECOs of charge $|z|=80$ of mass 1000 GeV.

Selection efficiency as a function of transverse kinetic energy $E^\text{kin}_\text{T}=E_\text{kin}\sin\theta$ and pseudorapidity $|\eta|$ for HECOs of charge $|z|=80$ of mass 1500 GeV.

Selection efficiency as a function of transverse kinetic energy $E^\text{kin}_\text{T}=E_\text{kin}\sin\theta$ and pseudorapidity $|\eta|$ for HECOs of charge $|z|=80$ of mass 2000 GeV.

Selection efficiency as a function of transverse kinetic energy $E^\text{kin}_\text{T}=E_\text{kin}\sin\theta$ and pseudorapidity $|\eta|$ for HECOs of charge $|z|=80$ of mass 2500 GeV.

Selection efficiency as a function of transverse kinetic energy $E^\text{kin}_\text{T}=E_\text{kin}\sin\theta$ and pseudorapidity $|\eta|$ for HECOs of charge $|z|=80$ of mass 3000 GeV.

Selection efficiency as a function of transverse kinetic energy $E^\text{kin}_\text{T}=E_\text{kin}\sin\theta$ and pseudorapidity $|\eta|$ for HECOs of charge $|z|=80$ of mass 4000 GeV.

Selection efficiency as a function of transverse kinetic energy $E^\text{kin}_\text{T}=E_\text{kin}\sin\theta$ and pseudorapidity $|\eta|$ for HECOs of charge $|z|=100$ of mass 200 GeV.

Selection efficiency as a function of transverse kinetic energy $E^\text{kin}_\text{T}=E_\text{kin}\sin\theta$ and pseudorapidity $|\eta|$ for HECOs of charge $|z|=100$ of mass 500 GeV.

Selection efficiency as a function of transverse kinetic energy $E^\text{kin}_\text{T}=E_\text{kin}\sin\theta$ and pseudorapidity $|\eta|$ for HECOs of charge $|z|=100$ of mass 1000 GeV.

Selection efficiency as a function of transverse kinetic energy $E^\text{kin}_\text{T}=E_\text{kin}\sin\theta$ and pseudorapidity $|\eta|$ for HECOs of charge $|z|=100$ of mass 1500 GeV.

Selection efficiency as a function of transverse kinetic energy $E^\text{kin}_\text{T}=E_\text{kin}\sin\theta$ and pseudorapidity $|\eta|$ for HECOs of charge $|z|=100$ of mass 2000 GeV.

Selection efficiency as a function of transverse kinetic energy $E^\text{kin}_\text{T}=E_\text{kin}\sin\theta$ and pseudorapidity $|\eta|$ for HECOs of charge $|z|=100$ of mass 2500 GeV.

Selection efficiency as a function of transverse kinetic energy $E^\text{kin}_\text{T}=E_\text{kin}\sin\theta$ and pseudorapidity $|\eta|$ for HECOs of charge $|z|=100$ of mass 3000 GeV.

Selection efficiency as a function of transverse kinetic energy $E^\text{kin}_\text{T}=E_\text{kin}\sin\theta$ and pseudorapidity $|\eta|$ for HECOs of charge $|z|=100$ of mass 4000 GeV.

The average tetrajet invariant mass distributions in data, along with the fitted background estimates for 0.16 < $\alpha$ < 0.18.

The average tetrajet invariant mass distributions in data, along with the fitted background estimates for 0.18 < $\alpha$ < 0.20.

The average tetrajet invariant mass distributions in data, along with the fitted background estimates for 0.20 < $\alpha$ < 0.22.

The average tetrajet invariant mass distributions in data, along with the fitted background estimates for 0.22 < $\alpha$ < 0.24.

The average tetrajet invariant mass distributions in data, along with the fitted background estimates for 0.24 < $\alpha$ < 0.26.

The average tetrajet invariant mass distributions in data, along with the fitted background estimates for 0.26 < $\alpha$ < 0.28.

The average tetrajet invariant mass distributions in data, along with the fitted background estimates for 0.28 < $\alpha$ < 0.30.

The average tetrajet invariant mass distributions in data, along with the fitted background estimates for 0.30 < $\alpha$ < 0.32.

The average tetrajet invariant mass distributions in data, along with the fitted background estimates for 0.32 < $\alpha$ < 0.34.

The average dijet invariant mass distributions in data, along with the fitted background estimates for 0.10 < $\alpha$ < 0.12.

The average dijet invariant mass distributions in data, along with the fitted background estimates for 0.12 < $\alpha$ < 0.14.

The average dijet invariant mass distributions in data, along with the fitted background estimates for 0.14 < $\alpha$ < 0.16.

The average dijet invariant mass distributions in data, along with the fitted background estimates for 0.16 < $\alpha$ < 0.18.

The average dijet invariant mass distributions in data, along with the fitted background estimates for 0.18 < $\alpha$ < 0.20.

The average dijet invariant mass distributions in data, along with the fitted background estimates for 0.20 < $\alpha$ < 0.22.

The average dijet invariant mass distributions in data, along with the fitted background estimates for 0.22 < $\alpha$ < 0.24.

The average dijet invariant mass distributions in data, along with the fitted background estimates for 0.24 < $\alpha$ < 0.26.

The average dijet invariant mass distributions in data, along with the fitted background estimates for 0.26 < $\alpha$ < 0.28.

The average dijet invariant mass distributions in data, along with the fitted background estimates for 0.28 < $\alpha$ < 0.30.

The average dijet invariant mass distributions in data, along with the fitted background estimates for 0.30 < $\alpha$ < 0.32.

The average dijet invariant mass distributions in data, along with the fitted background estimates for 0.32 < $\alpha$ < 0.34.

The expected and observed limits on $m_{4j}$ for signal templates using a 4-parameter fit function for 0.10 < $\alpha$ < 0.12.

The expected and observed limits on $m_{4j}$ for signal templates using a 4-parameter fit function for 0.12 < $\alpha$ < 0.14.

The expected and observed limits on $m_{4j}$ for signal templates using a 4-parameter fit function for 0.14 < $\alpha$ < 0.16.

The expected and observed limits on $m_{4j}$ for signal templates using a 4-parameter fit function for 0.16 < $\alpha$ < 0.18.

The expected and observed limits on $m_{4j}$ for signal templates using a 4-parameter fit function for 0.18 < $\alpha$ < 0.20.

The expected and observed limits on $m_{4j}$ for signal templates using a 4-parameter fit function for 0.20 < $\alpha$ < 0.22.

The expected and observed limits on $m_{4j}$ for signal templates using a 4-parameter fit function for 0.22 < $\alpha$ < 0.24.

The expected and observed limits on $m_{4j}$ for signal templates using a 4-parameter fit function for 0.24 < $\alpha$ < 0.26.

The expected and observed limits on $m_{4j}$ for signal templates using a 4-parameter fit function for 0.26 < $\alpha$ < 0.28.

The expected and observed limits on $m_{4j}$ for signal templates using a 4-parameter fit function for 0.28 < $\alpha$ < 0.30.

The expected and observed limits on $m_{4j}$ for signal templates using a 4-parameter fit function for 0.30 < $\alpha$ < 0.32.

The expected and observed limits on $m_{4j}$ for signal templates using a 4-parameter fit function for 0.32 < $\alpha$ < 0.34.

The expected and observed limits on $m_{\langle 2j \rangle}$ for signal templates using a 4-parameter fit function for 0.10 < $\alpha$ < 0.12.

The expected and observed limits on $m_{\langle 2j \rangle}$ for signal templates using a 4-parameter fit function for 0.12 < $\alpha$ < 0.14.

The expected and observed limits on $m_{\langle 2j \rangle}$ for signal templates using a 4-parameter fit function for 0.14 < $\alpha$ < 0.16.

The expected and observed limits on $m_{\langle 2j \rangle}$ for signal templates using a 4-parameter fit function for 0.16 < $\alpha$ < 0.18.

The expected and observed limits on $m_{\langle 2j \rangle}$ for signal templates using a 4-parameter fit function for 0.18 < $\alpha$ < 0.20.

The expected and observed limits on $m_{\langle 2j \rangle}$ for signal templates using a 4-parameter fit function for 0.20 < $\alpha$ < 0.22.

The expected and observed limits on $m_{\langle 2j \rangle}$ for signal templates using a 4-parameter fit function for 0.22 < $\alpha$ < 0.24.

The expected and observed limits on $m_{\langle 2j \rangle}$ for signal templates using a 4-parameter fit function for 0.24 < $\alpha$ < 0.26.

The expected and observed limits on $m_{\langle 2j \rangle}$ for signal templates using a 4-parameter fit function for 0.26 < $\alpha$ < 0.28.

The expected and observed limits on $m_{\langle 2j \rangle}$ for signal templates using a 4-parameter fit function for 0.28 < $\alpha$ < 0.30.

The expected and observed limits on $m_{\langle 2j \rangle}$ for signal templates using a 4-parameter fit function for 0.30 < $\alpha$ < 0.32.

The expected and observed limits on $m_{\langle 2j \rangle}$ for signal templates using a 4-parameter fit function for 0.32 < $\alpha$ < 0.34.

The expected and observed limits on $m_{4j}$ for Gaussian signal templates using a 4-parameter fit function for 0.10 < $\alpha$ < 0.12.

The expected and observed limits on $m_{4j}$ for Gaussian signal templates using a 4-parameter fit function for 0.12 < $\alpha$ < 0.14.

The expected and observed limits on $m_{4j}$ for Gaussian signal templates using a 4-parameter fit function for 0.14 < $\alpha$ < 0.16.

The expected and observed limits on $m_{4j}$ for Gaussian signal templates using a 4-parameter fit function for 0.16 < $\alpha$ < 0.18.

The expected and observed limits on $m_{4j}$ for Gaussian signal templates using a 4-parameter fit function for 0.18 < $\alpha$ < 0.20.

The expected and observed limits on $m_{4j}$ for Gaussian signal templates using a 4-parameter fit function for 0.20 < $\alpha$ < 0.22.

The expected and observed limits on $m_{4j}$ for Gaussian signal templates using a 4-parameter fit function for 0.22 < $\alpha$ < 0.24.

The expected and observed limits on $m_{4j}$ for Gaussian signal templates using a 4-parameter fit function for 0.24 < $\alpha$ < 0.26.

The expected and observed limits on $m_{4j}$ for Gaussian signal templates using a 4-parameter fit function for 0.26 < $\alpha$ < 0.28.

The expected and observed limits on $m_{4j}$ for Gaussian signal templates using a 4-parameter fit function for 0.28 < $\alpha$ < 0.30.

The expected and observed limits on $m_{4j}$ for Gaussian signal templates using a 4-parameter fit function for 0.30 < $\alpha$ < 0.32.

The expected and observed limits on $m_{4j}$ for Gaussian signal templates using a 4-parameter fit function for 0.32 < $\alpha$ < 0.34.

The expected and observed limits on $m_{\langle 2j \rangle}$ for Gaussian signal templates using a 4-parameter fit function for 0.10 < $\alpha$ < 0.12.

The expected and observed limits on $m_{\langle 2j \rangle}$ for Gaussian signal templates using a 4-parameter fit function for 0.12 < $\alpha$ < 0.14.

The expected and observed limits on $m_{\langle 2j \rangle}$ for Gaussian signal templates using a 4-parameter fit function for 0.14 < $\alpha$ < 0.16.

The expected and observed limits on $m_{\langle 2j \rangle}$ for Gaussian signal templates using a 4-parameter fit function for 0.16 < $\alpha$ < 0.18.

The expected and observed limits on $m_{\langle 2j \rangle}$ for Gaussian signal templates using a 4-parameter fit function for 0.18 < $\alpha$ < 0.20.

The expected and observed limits on $m_{\langle 2j \rangle}$ for Gaussian signal templates using a 4-parameter fit function for 0.20 < $\alpha$ < 0.22.

The expected and observed limits on $m_{\langle 2j \rangle}$ for Gaussian signal templates using a 4-parameter fit function for 0.22 < $\alpha$ < 0.24.

The expected and observed limits on $m_{\langle 2j \rangle}$ for Gaussian signal templates using a 4-parameter fit function for 0.24 < $\alpha$ < 0.26.

The expected and observed limits on $m_{\langle 2j \rangle}$ for Gaussian signal templates using a 4-parameter fit function for 0.26 < $\alpha$ < 0.28.

The expected and observed limits on $m_{\langle 2j \rangle}$ for Gaussian signal templates using a 4-parameter fit function for 0.28 < $\alpha$ < 0.30.

The expected and observed limits on $m_{\langle 2j \rangle}$ for Gaussian signal templates using a 4-parameter fit function for 0.30 < $\alpha$ < 0.32.

The expected and observed limits on $m_{\langle 2j \rangle}$ for Gaussian signal templates using a 4-parameter fit function for 0.32 < $\alpha$ < 0.34.

The product of the analysis acceptance and selection efficiency for all analysis selection criteria is shown as a function of $m_Y$ and $m_{X}/m_{Y}$.

The cross-section of the signal samples is shown as a function of $m_Y$ and $m_{X}/m_{Y}$.

Distribution of $m_{\mathrm{JJ}}$ in the SR (black points) compared to the background-only expectation, where the shape is derived from the CR and the normalisation is fitted in the SR.

Observed and expected limits at 95% CL on the cross-section times branching ratio for the production of dark quarks through a $Z^\prime$ as a function of the $Z^\prime$ mass for model A (see the text for details). The darker and lighter shaded bands around the expected limits represent the $\pm 1 \sigma$ and $\pm 2 \sigma$ uncertainty range, respectively.

The predicted theory cross-section times branching ratio for the production of dark quarks through a $Z^\prime$ as a function of the $Z^\prime$ mass for model A (see the text for details). In the theory model, a universal coupling of the $Z^\prime$ to each SM quark $g_q$ is assumed and the dark quarks are all represented by one species with an effective coupling $g_{q_d}$ to the $Z^\prime$. For model A, the predicted cross-section is shown for $g_q=0.05$ and $g_{q_d} = 0.2$.

Observed and expected limits at 95% CL on the cross-section times branching ratio for the production of dark quarks through a $Z^\prime$ as a function of the $Z^\prime$ mass for model B (see the text for details). The darker and lighter shaded bands around the expected limits represent the $\pm 1 \sigma$ and $\pm 2 \sigma$ uncertainty range, respectively.

The predicted theory cross-section times branching ratio for the production of dark quarks through a $Z^\prime$ as a function of the $Z^\prime$ mass for model B (see the text for details). In the theory model, a universal coupling of the $Z^\prime$ to each SM quark $g_q$ is assumed and the dark quarks are all represented by one species with an effective coupling $g_{q_d}$ to the $Z^\prime$. For model B, the predicted cross-section is shown for $g_q=0.15$ and $g_{q_d} = 0.5$.

Observed and expected limits at 95% CL on the cross-section times branching ratio for the production of dark quarks through a $Z^\prime$ as a function of the $Z^\prime$ mass for model C (see the text for details). The darker and lighter shaded bands around the expected limits represent the $\pm 1 \sigma$ and $\pm 2 \sigma$ uncertainty range, respectively.

The predicted theory cross-section times branching ratio for the production of dark quarks through a $Z^\prime$ as a function of the $Z^\prime$ mass for model C (see the text for details). In the theory model, a universal coupling of the $Z^\prime$ to each SM quark $g_q$ is assumed and the dark quarks are all represented by one species with an effective coupling $g_{q_d}$ to the $Z^\prime$. For model C, the predicted cross-section is shown for $g_q=0.15$ and $g_{q_d}= 0.5$.

Observed and expected limits at 95% CL on the cross-section times branching ratio for the production of dark quarks through a $Z^\prime$ as a function of the $Z^\prime$ mass for model D (see the text for details). The darker and lighter shaded bands around the expected limits represent the $\pm 1 \sigma$ and $\pm 2 \sigma$ uncertainty range, respectively.

The predicted theory cross-section times branching ratio for the production of dark quarks through a $Z^\prime$ as a function of the $Z^\prime$ mass for model D (see the text for details). In the theory model, a universal coupling of the $Z^\prime$ to each SM quark $g_q$ is assumed and the dark quarks are all represented by one species with an effective coupling $g_{q_d}$ to the $Z^\prime$. For model D, the predicted cross-section is shown for $g_q=0.05$ and $g_{q_d} = 0.2$.

Relative efficiency (in %) of each stage of the analysis selection with respect to the previous one for models A through D and $m_{Z^\prime}=2.5$ TeV. Approximately 90000 events were generated per signal mass hypothesis.

Distribution of the diphoton invariant mass in validation region VR2. The solid histograms are stacked to show the SM expectations after the 2&times;2D background estimation technique is applied. Background and signal predictions are normalised to the luminosity. The background category "h (other)" includes events originating from VBF, Vh, ggF, thq, thW and bb&#772;h, all subdominant in this signature. Statistical and systematic uncertainties are indicated by the shaded area. The lower panel of each plot shows the ratio of the data to the SM prediction for the respective bin. The first and last bins include the underflows and overflows respectively.

Distribution of the missing transverse momentum in validation region VR2. The solid histograms are stacked to show the SM expectations after the 2&times;2D background estimation technique is applied. Background and signal predictions are normalised to the luminosity. The background category "h (other)" includes events originating from VBF, Vh, ggF, thq, thW and bb&#772;h, all subdominant in this signature. Statistical and systematic uncertainties are indicated by the shaded area. The lower panel of each plot shows the ratio of the data to the SM prediction for the respective bin. The first and last bins include the underflows and overflows respectively.

Distribution of the diphoton invariant mass with all selections of the signal regions applied, except on m_{&gamma;&gamma;} itself, for signal region SR1h. The background estimation techniques described in the text are applied. The different backgrounds are stacked to add up to the total SM prediction in each bin. The predicted yields for signal benchmark models of varying &chi;&#771;⁰₁ mass are also overlaid (not stacked) assuming B(&chi;&#771;⁰₁ &rarr; hG&#771; ) to equal 100%. Background and signal predictions are normalised to the luminosity. The background category "h (other)" includes events originating from VBF, Vh, ggF, thq, thW and bb&#772;h, all subdominant in this signature. The sizes of the statistical and systematic uncertainties are indicated by the shaded areas. The lower panels show the ratio of the data to the SM prediction. Arrows indicate the borders of the signal region (|m_{&gamma;&gamma;}-125 GeV|&lt;5 GeV). The first and last bins include the underflows and overflows respectively.

Distribution of the diphoton invariant mass with all selections of the signal regions applied, except on m_{&gamma;&gamma;} itself, for signal region SR1Z. The background estimation techniques described in the text are applied. The different backgrounds are stacked to add up to the total SM prediction in each bin. The predicted yields for signal benchmark models of varying &chi;&#771;⁰₁ mass are also overlaid (not stacked) assuming B(&chi;&#771;⁰₁ &rarr; hG&#771; ) to equal 50%. Background and signal predictions are normalised to the luminosity. The background category "h (other)" includes events originating from VBF, Vh, ggF, thq, thW and bb&#772;h, all subdominant in this signature. The sizes of the statistical and systematic uncertainties are indicated by the shaded areas. The lower panels show the ratio of the data to the SM prediction. Arrows indicate the borders of the signal region (|m_{&gamma;&gamma;}-125 GeV|&lt;5 GeV). The first and last bins include the underflows and overflows respectively.

Distribution of the diphoton invariant mass with all selections of the signal regions applied, except on m_{&gamma;&gamma;} itself, for signal region SR2. The background estimation techniques described in the text are applied. The different backgrounds are stacked to add up to the total SM prediction in each bin. The predicted yields for signal benchmark models of varying &chi;&#771;⁰₁ mass are also overlaid (not stacked) assuming B(&chi;&#771;⁰₁ &rarr; hG&#771; ) to equal 100%. Background and signal predictions are normalised to the luminosity. The background category "h (other)" includes events originating from VBF, Vh, ggF, thq, thW and bb&#772;h, all subdominant in this signature. The sizes of the statistical and systematic uncertainties are indicated by the shaded areas. The lower panels show the ratio of the data to the SM prediction. Arrows indicate the borders of the signal region (|m_{&gamma;&gamma;}-125 GeV|&lt;5 GeV). The first and last bins include the underflows and overflows respectively.

Observed and expected limits on the pure higgsino cross-section at 95% CL assuming B(&chi;&#771;⁰₁ &rarr; hG&#771; )=100% for different &chi;&#771;⁰₁ masses, obtained by a statistical combination of the three signal regions SR1h, SR1Z and SR2. The inner and outer bands indicate the 1&sigma; and 2&sigma; variation on the expected limit respectively.

Observed and expected 95% CL limits on the pure-higgsino branching fraction to B(&chi;&#771;⁰₁ &rarr; hG&#771; ) as a function of the higgsino mass m(&chi;&#771;⁰₁) assuming it decays via either &chi;&#771;⁰₁&rarr; hG&#771; or &chi;&#771;⁰₁&rarr; ZG&#771;. Limits are obtained by performing a statistical combination of the three signal regions SR1h, SR1Z and SR2. The &plusmn; 1&sigma; variation on the expected limit is shown. The dotted lines indicate the observed limit obtained by a variation of theoretical prediction for the neutralino production cross-section by &plusmn;1 &sigma;. Values of B(&chi;&#771;⁰₁ &rarr; hG&#771; ) larger than the observed 95% CL limit are excluded, as indicated by the hatched area.

Observed and expected 95% CL limits on the pure-higgsino branching fraction to B(&chi;&#771;⁰₁ &rarr; hG&#771; ) as a function of the higgsino mass m(&chi;&#771;⁰₁) assuming it decays via either &chi;&#771;⁰₁&rarr; hG&#771; or &chi;&#771;⁰₁&rarr; ZG&#771;. Limits are obtained by performing a statistical combination of the three signal regions SR1h, SR1Z and SR2. The &plusmn; 1&sigma; variation on the expected limit is shown. The dotted lines indicate the observed limit obtained by a variation of theoretical prediction for the neutralino production cross-section by &plusmn;1 &sigma;. Values of B(&chi;&#771;⁰₁ &rarr; hG&#771; ) larger than the observed 95% CL limit are excluded, as indicated by the hatched area.

Observed and expected 95% CL limits on the pure-higgsino branching fraction to B(&chi;&#771;⁰₁ &rarr; hG&#771; ) as a function of the higgsino mass m(&chi;&#771;⁰₁) assuming it decays via either &chi;&#771;⁰₁&rarr; hG&#771; or &chi;&#771;⁰₁&rarr; ZG&#771;. Limits are obtained by performing a statistical combination of the three signal regions SR1h, SR1Z and SR2. The &plusmn; 1&sigma; variation on the expected limit is shown. The dotted lines indicate the observed limit obtained by a variation of theoretical prediction for the neutralino production cross-section by &plusmn;1 &sigma;. Values of B(&chi;&#771;⁰₁ &rarr; hG&#771; ) larger than the observed 95% CL limit are excluded, as indicated by the hatched area.

Observed and expected 95% CL limits on the pure-higgsino branching fraction to B(&chi;&#771;⁰₁ &rarr; hG&#771; ) as a function of the higgsino mass m(&chi;&#771;⁰₁) assuming it decays via either &chi;&#771;⁰₁&rarr; hG&#771; or &chi;&#771;⁰₁&rarr; ZG&#771;. Limits are obtained by performing a statistical combination of the three signal regions SR1h, SR1Z and SR2. The &plusmn; 1&sigma; variation on the expected limit is shown. The dotted lines indicate the observed limit obtained by a variation of theoretical prediction for the neutralino production cross-section by &plusmn;1 &sigma;. Values of B(&chi;&#771;⁰₁ &rarr; hG&#771; ) larger than the observed 95% CL limit are excluded, as indicated by the hatched area.

Observed and expected 95% CL limits on the pure-higgsino branching fraction to B(&chi;&#771;⁰₁ &rarr; hG&#771; ) as a function of the higgsino mass m(&chi;&#771;⁰₁) assuming it decays via either &chi;&#771;⁰₁&rarr; hG&#771; or &chi;&#771;⁰₁&rarr; ZG&#771;. Limits are obtained by performing a statistical combination of the three signal regions SR1h, SR1Z and SR2. The &plusmn; 1&sigma; variation on the expected limit is shown. The dotted lines indicate the observed limit obtained by a variation of theoretical prediction for the neutralino production cross-section by &plusmn;1 &sigma;. Values of B(&chi;&#771;⁰₁ &rarr; hG&#771; ) larger than the observed 95% CL limit are excluded, as indicated by the hatched area.

Observed and expected 95% CL limits on the pure-higgsino branching fraction to B(&chi;&#771;⁰₁ &rarr; hG&#771; ) as a function of the higgsino mass m(&chi;&#771;⁰₁) assuming it decays via either &chi;&#771;⁰₁&rarr; hG&#771; or &chi;&#771;⁰₁&rarr; ZG&#771;. Limits are obtained by performing a statistical combination of the three signal regions SR1h, SR1Z and SR2. The &plusmn; 1&sigma; variation on the expected limit is shown. The dotted lines indicate the observed limit obtained by a variation of theoretical prediction for the neutralino production cross-section by &plusmn;1 &sigma;. Values of B(&chi;&#771;⁰₁ &rarr; hG&#771; ) larger than the observed 95% CL limit are excluded, as indicated by the hatched area.

Numbers of signal and background events in the signal regions. The respective background estimation techniques are applied. The background category "h (other)" includes events originating from VBF, Vh, ggF, thq, thW and bb&#772;h, all subdominant in this signature. The different backgrounds are stacked to add up to the total Standard Model prediction in each bin. The predicted yields for signal benchmark models of varying &chi;&#771;⁰₁ mass are also plotted (not stacked), assuming B(&chi;&#771;⁰₁ &rarr; hG&#771; )=100% and a &chi;&#771;⁰₁ mass of 130 or 200 GeV. The statistical and systematic uncertainties are indicated by the shaded areas in the top plot. The bottom panel shows the statistical significance&nbsp;<a href="https://atlas.web.cern.ch/Atlas/GROUPS/PHYSICS/PUBNOTES/ATL-PHYS-PUB-2020-025/">[Ref]</a> of the difference between the SM prediction and the observed data in each region.

Observed 95% CL limits in pb on the pure higgsino cross-section, shown in the m(&chi;&#771;⁰₁)-B(&chi;&#771;⁰₁ &rarr; hG&#771; ) plane. Limits are obtained by a statistical combination of the three signal regions SR1h, SR1Z and SR2, assuming the neutralino to decay via either &chi;&#771;⁰₁&rarr; hG&#771; or &chi;&#771;⁰₁&rarr; ZG&#771;.

Acceptances for all signal model points considered in the analysis, shown in the m(&chi;&#771;⁰₁)-B(&chi;&#771;⁰₁ &rarr; hG&#771; ) plane. Acceptances are provided separately for signal region SR1h, assuming the neutralino to decay via either &chi;&#771;⁰₁&rarr; hG&#771; or &chi;&#771;⁰₁&rarr; ZG&#771;.

Acceptances for all signal model points considered in the analysis, shown in the m(&chi;&#771;⁰₁)-B(&chi;&#771;⁰₁ &rarr; hG&#771; ) plane. Acceptances are provided separately for signal region SR1Z, assuming the neutralino to decay via either &chi;&#771;⁰₁&rarr; hG&#771; or &chi;&#771;⁰₁&rarr; ZG&#771;.

Acceptances for all signal model points considered in the analysis, shown in the m(&chi;&#771;⁰₁)-B(&chi;&#771;⁰₁ &rarr; hG&#771; ) plane. Acceptances are provided separately for signal region SR2, assuming the neutralino to decay via either &chi;&#771;⁰₁&rarr; hG&#771; or &chi;&#771;⁰₁&rarr; ZG&#771;.

Efficiencies for all signal model points considered in the analysis, shown in the m(&chi;&#771;⁰₁)-B(&chi;&#771;⁰₁ &rarr; hG&#771; ) plane. Efficiencies are provided separately for signal region SR1h, assuming the neutralino to decay via either &chi;&#771;⁰₁&rarr; hG&#771; or &chi;&#771;⁰₁&rarr; ZG&#771;.

Efficiencies for all signal model points considered in the analysis, shown in the m(&chi;&#771;⁰₁)-B(&chi;&#771;⁰₁ &rarr; hG&#771; ) plane. Efficiencies are provided separately for signal region SR1Z, assuming the neutralino to decay via either &chi;&#771;⁰₁&rarr; hG&#771; or &chi;&#771;⁰₁&rarr; ZG&#771;.

Efficiencies for all signal model points considered in the analysis, shown in the m(&chi;&#771;⁰₁)-B(&chi;&#771;⁰₁ &rarr; hG&#771; ) plane. Efficiencies are provided separately for signal region SR2, assuming the neutralino to decay via either &chi;&#771;⁰₁&rarr; hG&#771; or &chi;&#771;⁰₁&rarr; ZG&#771;.

Observed and expected numbers of events in the three signal regions. The background category "h (other)" includes events originating from VBF, Vh, ggF, thq, thW and bb&#772;h, all subdominant in this signature. The table also includes model-independent 95% CL upper limits on the visible number of BSM events (S⁹⁵_obs), the number of BSM events given the expected number of background events (S⁹⁵_exp) and the visible BSM cross-section (&lang;&epsilon; &sigma;&rang;_obs⁹⁵), all calculated from pseudo-experiments. The discovery p-value (p₀) is also shown and its value is capped at 0.5 if the observed number of events is below the expected number of events.

Cutflows of two benchmark signal points assuming B(&chi;&#771;⁰₁ &rarr; hG&#771; )=100% for all three discovery signal regions. The initial selection includes the leptons veto. Only statistical uncertainties are included. Expected yields are normalised to a luminosity of 139 fb^-1.