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}$.

Event yields after a background-only fit and data in the three signal regions. All the backgrounds except the fake $\tau_{\mathrm{had}}$ are estimated using the real $\tau_{\mathrm{had}}$ contributions from MC. The fake $\tau_{\mathrm{had}}$ background contribution is data-driven. The top-quark, $V$+jets, and small backgrounds are included in "Others". The total uncertainties include both systematic and statistical. In the $WW1\ell2\tau_{\mathrm{had}}$ channel, the smaller uncertainty in the total background than the fake $\tau_{\mathrm{had}}$ is due to the anticorrelations between the nuisance parameters of the fake $\tau_{\mathrm{had}}$ and diboson backgrounds obtained in the fit.

Uncertainties on $\sigma(pp\to X\to SH)$ obtained from the combined fit to the data for the least and most sensitive mass points. The uncertainties are symmetrised and grouped into the categories described in the systematic section.

Summary of observed and expected 95\% CL upper limits on the production cross-section of $\sigma(pp\to X\to SH)$ in fb in the combination of $WW$+$ZZ$, and in the $WW$ and $ZZ$ individual channels. The best expected limit is set to 85 fb for $m_{X}$ = 1250 GeV and $_{S}$ = 300 GeV. The combination limits gain about 26--53% over the best individual limit from the $WW1\ell2\tau_{\mathrm{had}}$ channel.

Values of $\Delta Z$ for the discovery sensitivity, as defined in the text, as a function of the invariant mass $\textit{m}$. The j+b-jet invariant mass distribution is calculated in the 10 pb AR. Positive percentages indicate improvements in sensitivity. Horizontal dashed lines are drawn at 100% and 200% to guide the eye. The five benchmark BSM models are (1) charged Higgs boson production in association with a top quark, $tbH^{+}$ with $H^{+} \rightarrow t\bar{b}$; (2) a Kaluza-Klein gauge boson, $W_{KK}$, with the SM $W$ boson and a radion $\phi$; (3) a $Z'$ boson decaying to a composite lepton $E$ and $\ell$, with $E \rightarrow Z\ell$; (4) the sequential standard model $W' \rightarrow W Z' \rightarrow \ell\nu q\bar{q}$; (5) a simplified dark-matter model with an axial-vector mediator $Z' \rightarrow q\bar{q}$, where one of the quarks radiates a $W$ boson decaying to $\ell\nu$. The multiple markers shown for the composite-lepton model at the same invariant mass values correspond to different composite lepton ($E$) masses between 0.25 and 3.5 TeV. The center positions of the markers are set to the masses of the corresponding heavy particles.

Values of $\Delta Z$ for the discovery sensitivity, as defined in the text, as a function of the invariant mass $\textit{m}$. The 2b-jet invariant mass distribution is calculated in the 10 pb AR. Positive percentages indicate improvements in sensitivity. Horizontal dashed lines are drawn at 100% and 200% to guide the eye. The five benchmark BSM models are (1) charged Higgs boson production in association with a top quark, $tbH^{+}$ with $H^{+} \rightarrow t\bar{b}$; (2) a Kaluza-Klein gauge boson, $W_{KK}$, with the SM $W$ boson and a radion $\phi$; (3) a $Z'$ boson decaying to a composite lepton $E$ and $\ell$, with $E \rightarrow Z\ell$; (4) the sequential standard model $W' \rightarrow W Z' \rightarrow \ell\nu q\bar{q}$; (5) a simplified dark-matter model with an axial-vector mediator $Z' \rightarrow q\bar{q}$, where one of the quarks radiates a $W$ boson decaying to $\ell\nu$. The multiple markers shown for the composite-lepton model at the same invariant mass values correspond to different composite lepton ($E$) masses between 0.25 and 3.5 TeV. The center positions of the markers are set to the masses of the corresponding heavy particles.

Values of $\Delta Z$ for the discovery sensitivity, as defined in the text, as a function of the invariant mass $\textit{m}$. The j+e invariant mass distribution is calculated in the 10 pb AR. Positive percentages indicate improvements in sensitivity. Horizontal dashed lines are drawn at 100% and 200% to guide the eye. The five benchmark BSM models are (1) charged Higgs boson production in association with a top quark, $tbH^{+}$ with $H^{+} \rightarrow t\bar{b}$; (2) a Kaluza-Klein gauge boson, $W_{KK}$, with the SM $W$ boson and a radion $\phi$; (3) a $Z'$ boson decaying to a composite lepton $E$ and $\ell$, with $E \rightarrow Z\ell$; (4) the sequential standard model $W' \rightarrow W Z' \rightarrow \ell\nu q\bar{q}$; (5) a simplified dark-matter model with an axial-vector mediator $Z' \rightarrow q\bar{q}$, where one of the quarks radiates a $W$ boson decaying to $\ell\nu$. The multiple markers shown for the composite-lepton model at the same invariant mass values correspond to different composite lepton ($E$) masses between 0.25 and 3.5 TeV. The center positions of the markers are set to the masses of the corresponding heavy particles.

Values of $\Delta Z$ for the discovery sensitivity, as defined in the text, as a function of the invariant mass $\textit{m}$. The b-jet+e invariant mass distribution is calculated in the 10 pb AR. Positive percentages indicate improvements in sensitivity. Horizontal dashed lines are drawn at 100% and 200% to guide the eye. The five benchmark BSM models are (1) charged Higgs boson production in association with a top quark, $tbH^{+}$ with $H^{+} \rightarrow t\bar{b}$; (2) a Kaluza-Klein gauge boson, $W_{KK}$, with the SM $W$ boson and a radion $\phi$; (3) a $Z'$ boson decaying to a composite lepton $E$ and $\ell$, with $E \rightarrow Z\ell$; (4) the sequential standard model $W' \rightarrow W Z' \rightarrow \ell\nu q\bar{q}$; (5) a simplified dark-matter model with an axial-vector mediator $Z' \rightarrow q\bar{q}$, where one of the quarks radiates a $W$ boson decaying to $\ell\nu$. The multiple markers shown for the composite-lepton model at the same invariant mass values correspond to different composite lepton ($E$) masses between 0.25 and 3.5 TeV. The center positions of the markers are set to the masses of the corresponding heavy particles.

Values of $\Delta Z$ for the discovery sensitivity, as defined in the text, as a function of the invariant mass $\textit{m}$. The j+$\gamma$ invariant mass distribution is calculated in the 10 pb AR. Positive percentages indicate improvements in sensitivity. Horizontal dashed lines are drawn at 100% and 200% to guide the eye. The five benchmark BSM models are (1) charged Higgs boson production in association with a top quark, $tbH^{+}$ with $H^{+} \rightarrow t\bar{b}$; (2) a Kaluza-Klein gauge boson, $W_{KK}$, with the SM $W$ boson and a radion $\phi$; (3) a $Z'$ boson decaying to a composite lepton $E$ and $\ell$, with $E \rightarrow Z\ell$; (4) the sequential standard model $W' \rightarrow W Z' \rightarrow \ell\nu q\bar{q}$; (5) a simplified dark-matter model with an axial-vector mediator $Z' \rightarrow q\bar{q}$, where one of the quarks radiates a $W$ boson decaying to $\ell\nu$. The multiple markers shown for the composite-lepton model at the same invariant mass values correspond to different composite lepton ($E$) masses between 0.25 and 3.5 TeV. The center positions of the markers are set to the masses of the corresponding heavy particles.

Values of $\Delta Z$ for the discovery sensitivity, as defined in the text, as a function of the invariant mass $\textit{m}$. The j+$\mu$ invariant mass distribution is calculated in the 10 pb AR. Positive percentages indicate improvements in sensitivity. Horizontal dashed lines are drawn at 100% and 200% to guide the eye. The five benchmark BSM models are (1) charged Higgs boson production in association with a top quark, $tbH^{+}$ with $H^{+} \rightarrow t\bar{b}$; (2) a Kaluza-Klein gauge boson, $W_{KK}$, with the SM $W$ boson and a radion $\phi$; (3) a $Z'$ boson decaying to a composite lepton $E$ and $\ell$, with $E \rightarrow Z\ell$; (4) the sequential standard model $W' \rightarrow W Z' \rightarrow \ell\nu q\bar{q}$; (5) a simplified dark-matter model with an axial-vector mediator $Z' \rightarrow q\bar{q}$, where one of the quarks radiates a $W$ boson decaying to $\ell\nu$. The multiple markers shown for the composite-lepton model at the same invariant mass values correspond to different composite lepton ($E$) masses between 0.25 and 3.5 TeV. The center positions of the markers are set to the masses of the corresponding heavy particles.

Values of $\Delta Z$ for the discovery sensitivity, as defined in the text, as a function of the invariant mass $\textit{m}$. The b-jet+$\mu$ invariant mass distribution is calculated in the 10 pb AR. Positive percentages indicate improvements in sensitivity. Horizontal dashed lines are drawn at 100% and 200% to guide the eye. The five benchmark BSM models are (1) charged Higgs boson production in association with a top quark, $tbH^{+}$ with $H^{+} \rightarrow t\bar{b}$; (2) a Kaluza-Klein gauge boson, $W_{KK}$, with the SM $W$ boson and a radion $\phi$; (3) a $Z'$ boson decaying to a composite lepton $E$ and $\ell$, with $E \rightarrow Z\ell$; (4) the sequential standard model $W' \rightarrow W Z' \rightarrow \ell\nu q\bar{q}$; (5) a simplified dark-matter model with an axial-vector mediator $Z' \rightarrow q\bar{q}$, where one of the quarks radiates a $W$ boson decaying to $\ell\nu$. The multiple markers shown for the composite-lepton model at the same invariant mass values correspond to different composite lepton ($E$) masses between 0.25 and 3.5 TeV. The center positions of the markers are set to the masses of the corresponding heavy particles.

Values of $\Delta Z$ for the discovery sensitivity, as defined in the text, as a function of the invariant mass $\textit{m}$. The b-jet+$\gamma$ invariant mass distribution is calculated in the 10 pb AR. Positive percentages indicate improvements in sensitivity. Horizontal dashed lines are drawn at 100% and 200% to guide the eye. The five benchmark BSM models are (1) charged Higgs boson production in association with a top quark, $tbH^{+}$ with $H^{+} \rightarrow t\bar{b}$; (2) a Kaluza-Klein gauge boson, $W_{KK}$, with the SM $W$ boson and a radion $\phi$; (3) a $Z'$ boson decaying to a composite lepton $E$ and $\ell$, with $E \rightarrow Z\ell$; (4) the sequential standard model $W' \rightarrow W Z' \rightarrow \ell\nu q\bar{q}$; (5) a simplified dark-matter model with an axial-vector mediator $Z' \rightarrow q\bar{q}$, where one of the quarks radiates a $W$ boson decaying to $\ell\nu$. The multiple markers shown for the composite-lepton model at the same invariant mass values correspond to different composite lepton ($E$) masses between 0.25 and 3.5 TeV. The center positions of the markers are set to the masses of the corresponding heavy particles.

The 95% CL upper limits on the cross section times acceptance ($A$), efficiency ($\epsilon$), and branching ratio ($B$) for Gaussian-shaped signals with different signal widths. The limits are calculated for events with at least one lepton with pT > 60 GeV in the 10 pb AR. Two width hypotheses are shown; $\sigma X / mX = 0$ and $\sigma X / mX = 0.15$. In both cases, the detector resolution for jets is included in the simulation of signal samples. The $\pm1 \sigma$ and $\pm2 \sigma$ bands around the expected limit are shown for $\sigma X / mX = 0$ signals. Mass points are spaced 5% apart, relative to the preceding point, starting at 0.3 TeV.

Distributions of the anomaly score from the AE for data and five benchmark BSM models. Their legends, from top to bottom, are; (1) charged Higgs boson production in association with a top quark, $tbH^{+}$ with $H^{+} \rightarrow t\bar{b}$, with the $H^{+}$ mass of 2 TeV; (2) a 6 TeV Kaluza-Klein gauge boson, $W_{KK}$, with the SM $W$ boson and a radion $\phi$; (3) a 4 TeV $Z'$ boson decaying to a composite lepton $E$ with a mass of 0.5 TeV; (4) the sequential standard model $W' \rightarrow WZ' \rightarrow \ell\nu q\bar{q}$, with the $W'$ mass of 6.2 TeV and the $Z'$ mass of 6 TeV; (5) a simplified dark-matter model with an axial-vector mediator $Z' \rightarrow q\bar{q}$, where one of the quarks radiates a $W$ boson decaying to $\ell\nu$ with a mass of 6 TeV. The distributions for the BSM models are smoothed to remove fluctuations due to low MC event counts. The vertical lines indicate the start of the three anomaly regions (ARs). The labels of the three ARs indicate the visible cross section for hypothetical processes yielding the same number of events as observed in the 140 $fb^{-1}$ dataset. The AE is applied to preselected events without any requirements on invariant mass distributions.

Invariant mass distributions of $j+Y$ for $M_{jY}$ > 0.3 TeV after preselection along with the fit from Eq.(1). The fit is represented by red lines, and the associated uncertainties are indicated by yellow bands. The lower panels show the bin-by-bin significances of deviations from the fit, calculated as $(d_{\textit{i}} - f_{\textit{i}}) / \delta_{\textit{i}}$, where $d_{\textit{i}}$ is the data yield, $f_{\textit{i}}$ is the fit value, and $\delta_{i}$ is the uncertainty of data in the $\textit{i}$-th bin

Distributions of the anomaly score for data and several anomaly scenarios. The example BSM model (shown with the dashed blue lines) is the sequential standard model $W' \rightarrow WZ' \rightarrow \ell\nu q\bar{q}$. The mass of $W'$ is set to 2.2 TeV and the mass of the $Z'$ is set to 2 TeV. This model leads to the final state of one lepton, two jets, and small missing transverse energy that is similar to the SM backgrounds. The other histograms represent events from the same model with artificial modifications to represent anomalous events; (1) "Anomaly 1" is the case where all jets beyond the second jet ($N_{j}$ > 2) are replaced with photons; (2) "Anomaly 2" is the case where all jets beyond the second jet are replaced with $b$-jets; (3) "Anomaly 3" is the case of low-multiplicity events where all jets beyond the second jet are removed. The histograms are normalized to unit. The left peak near log(Loss) = -9 visible in Anomaly 1 and 2 is from the events without a third jet.

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 and expected upper limits at 95% CL on the production cross-section times branching fraction as a function of $m_{A'}$ at dark Higgs boson mass of 50 GeV

Observed and expected upper limits at 95% CL on the production cross-section times branching fraction as a function of $m_{A'}$ at dark Higgs boson mass of 60 GeV

Observed and expected upper limits at 95% CL on the production cross-section times branching fraction as a function of $m_{A'}$ at dark Higgs boson mass of 70 GeV

Observed 90% CL upper limits on $\alpha_{D}\epsilon^{2}$ as a function of $m_{A'}$ at dark Higgs boson mass of 20 GeV

Observed 90% CL upper limits on $\alpha_{D}\epsilon^{2}$ as a function of $m_{A'}$ at dark Higgs boson mass of 30 GeV

Observed 90% CL upper limits on $\alpha_{D}\epsilon^{2}$ as a function of $m_{A'}$ at dark Higgs boson mass of 40 GeV

Observed 90% CL upper limits on $\alpha_{D}\epsilon^{2}$ as a function of $m_{A'}$ at dark Higgs boson mass of 50 GeV

Observed 90% CL upper limits on $\alpha_{D}\epsilon^{2}$ as a function of $m_{A'}$ at dark Higgs boson mass of 60 GeV

Observed 90% CL upper limits on $\alpha_{D}\epsilon^{2}$ as a function of $m_{A'}$ at dark Higgs boson mass of 70 GeV

Observed 90% CL upper limits on $\epsilon^{2}$, assuming of $\alpha_{D}=0.1$, as a function of $m_{A'}$ at dark Higgs boson mass of 20 GeV

Observed 90% CL upper limits on $\epsilon^{2}$, assuming of $\alpha_{D}=0.1$, as a function of $m_{A'}$ at dark Higgs boson mass of 30 GeV

Observed 90% CL upper limits on $\epsilon^{2}$, assuming of $\alpha_{D}=0.1$, as a function of $m_{A'}$ at dark Higgs boson mass of 40 GeV

Observed 90% CL upper limits on $\epsilon^{2}$, assuming of $\alpha_{D}=0.1$, as a function of $m_{A'}$ at dark Higgs boson mass of 50 GeV

Observed 90% CL upper limits on $\epsilon^{2}$, assuming of $\alpha_{D}=0.1$, as a function of $m_{A'}$ at dark Higgs boson mass of 60 GeV

Observed 90% CL upper limits on $\epsilon^{2}$, assuming of $\alpha_{D}=0.1$, as a function of $m_{A'}$ at dark Higgs boson mass of 70 GeV

Observed and expected upper limits at 95% CL on the branching fraction as a function of $m_{A'}$ at dark Higgs boson mass of 20 GeV

Observed and expected upper limits at 95% CL on the branching fraction as a function of $m_{A'}$ at dark Higgs boson mass of 30 GeV

Observed and expected upper limits at 95% CL on the branching fraction as a function of $m_{A'}$ at dark Higgs boson mass of 40 GeV

Observed and expected upper limits at 95% CL on the branching fraction as a function of $m_{A'}$ at dark Higgs boson mass of 50 GeV

Observed and expected upper limits at 95% CL on the branching fraction as a function of $m_{A'}$ at dark Higgs boson mass of 60 GeV

Observed and expected upper limits at 95% CL on the branching fraction as a function of $m_{A'}$ at dark Higgs boson mass of 70 GeV

Observed and expected upper limits at 95% CL on $\alpha_{D}\epsilon^{2}$ as a function of $m_{A'}$ at dark Higgs boson mass of 20 GeV

Observed and expected upper limits at 95% CL on $\alpha_{D}\epsilon^{2}$ as a function of $m_{A'}$ at dark Higgs boson mass of 30 GeV

Observed and expected upper limits at 95% CL on $\alpha_{D}\epsilon^{2}$ as a function of $m_{A'}$ at dark Higgs boson mass of 40 GeV

Observed and expected upper limits at 95% CL on $\alpha_{D}\epsilon^{2}$ as a function of $m_{A'}$ at dark Higgs boson mass of 50 GeV

Observed and expected upper limits at 95% CL on $\alpha_{D}\epsilon^{2}$ as a function of $m_{A'}$ at dark Higgs boson mass of 60 GeV

Observed and expected upper limits at 95% CL on $\alpha_{D}\epsilon^{2}$ as a function of $m_{A'}$ at dark Higgs boson mass of 70 GeV

Acceptance times efficiency for signal grid in merged two-prong signal region.

Acceptance times efficiency for signal grid in resolved two-prong signal region.

Acceptance times efficiency for signal grid in resolved two-prong signal region.

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

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

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

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

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

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

2 sigma band of expected combination limits at 95% CL in the $(m_{a},m_{A})$ plane under the assumption of $sin\theta$ = 0.35.

Observed combination limits at 95% CL in the $(m_{a},m_{A})$ plane under the assumption of $sin\theta$ = 0.7.

Expected combination limits at 95% CL in the $(m_{a},m_{A})$ plane under the assumption of $sin\theta$ = 0.7.

1 sigma band of expected combination limits at 95% CL in the $(m_{a},m_{A})$ plane under the assumption of $sin\theta$ = 0.7.

2 sigma band of expected combination limits at 95% CL in the $(m_{a},m_{A})$ plane under the assumption of $sin\theta$ = 0.7.

Observed limits of H(inv) at 95% CL in the $(m_{a},m_{\chi})$ plane under the assumption of $m_{A}$ = 1.2 TeV, $tan\beta$ = 1.0 and $sin\theta$ = 0.35.

Expected limits of H(inv) at 95% CL in the $(m_{a},m_{\chi})$ plane under the assumption of $m_{A}$ = 1.2 TeV, $tan\beta$ = 1.0 and $sin\theta$ = 0.35.

Observed limits of monoh(bb) at 95% CL in the $(m_{a},m_{\chi})$ plane under the assumption of $m_{A}$ = 1.2 TeV, $tan\beta$ = 1.0 and $sin\theta$ = 0.35.

Expected limits of monoh(bb) at 95% CL in the $(m_{a},m_{\chi})$ plane under the assumption of $m_{A}$ = 1.2 TeV, $tan\beta$ = 1.0 and $sin\theta$ = 0.35.

Observed limits of $h\to aa \to \mu\mu\tau\tau$ at 95% CL in the $(m_{a},m_{\chi})$ plane under the assumption of $m_{A}$ = 1.2 TeV, $tan\beta$ = 1.0 and $sin\theta$ = 0.35. Data points included are corresponding to the lower bound of the exclusion contour. The region above the bound was excluded.

Observed limits of $h\to aa \to bbbb$ at 95% CL in the $(m_{a},m_{\chi})$ plane under the assumption of $m_{A}$ = 1.2 TeV, $tan\beta$ = 1.0 and $sin\theta$ = 0.35. Data points included are corresponding to the lower bound of the exclusion contour. The region above the bound was excluded.

Observed limits of $h\to aa \to bb\mu\mu$ at 95% CL in the $(m_{a},m_{\chi})$ plane under the assumption of $m_{A}$ = 1.2 TeV, $tan\beta$ = 1.0 and $sin\theta$ = 0.35. Data points included are corresponding to the lower bound of the exclusion contour. The region above the bound was excluded.

Observed limits of $h\to aa \to 4 \mu$ at 95% CL in the $(m_{a},m_{\chi})$ plane under the assumption of $m_{A}$ = 1.2 TeV, $tan\beta$ = 1.0 and $sin\theta$ = 0.35. Data points included are corresponding to the lower bound of the exclusion contour. The region above the bound was excluded.

Observed limits of $h\to aa \to 4 \mu$ at 95% CL in the $(m_{a},m_{\chi})$ plane under the assumption of $m_{A}$ = 1.2 TeV, $tan\beta$ = 1.0 and $sin\theta$ = 0.35. Data points included are corresponding to the lower bound of the exclusion contour. The region above the bound was excluded.

Observed limits of $h\to aa \to 4 \mu$ at 95% CL in the $(m_{a},m_{\chi})$ plane under the assumption of $m_{A}$ = 1.2 TeV, $tan\beta$ = 1.0 and $sin\theta$ = 0.35. Data points included are corresponding to the lower bound of the exclusion contour. The region above the bound was excluded.

Observed combination limits at 95% CL in the $(m_{A},tan\beta)$ plane under the assumption of $sin\theta$ = 0.35.

Expected combination limits at 95% CL in the $(m_{A},tan\beta)$ plane under the assumption of $sin\theta$ = 0.35.

1 sigma band of expected combination limits at 95% CL in the $(m_{A},tan\beta)$ plane under the assumption of $sin\theta$ = 0.35.

2 sigma band of expected combination limits at 95% CL in the $(m_{A},tan\beta)$ plane under the assumption of $sin\theta$ = 0.35.

Observed combination limits at 95% CL in the $(m_{A},tan\beta)$ plane under the assumption of $sin\theta$ = 0.7.

Expected combination limits at 95% CL in the $(m_{A},tan\beta)$ plane under the assumption of $sin\theta$ = 0.7.

1 sigma band of expected combination limits at 95% CL in the $(m_{A},tan\beta)$ plane under the assumption of $sin\theta$ = 0.7.

2 sigma band of expected combination limits at 95% CL in the $(m_{A},tan\beta)$ plane under the assumption of $sin\theta$ = 0.7.

Observed combination limits at 95% CL in the $(m_{a},tan\beta)$ plane under the assumption of $sin\theta$ = 0.35.

Expected combination limits at 95% CL in the $(m_{a},tan\beta)$ plane under the assumption of $sin\theta$ = 0.35.

1 sigma band of expected combination limits at 95% CL in the $(m_{a},tan\beta)$ plane under the assumption of $sin\theta$ = 0.35.

2 sigma band of expected combination limits at 95% CL in the $(m_{a},tan\beta)$ plane under the assumption of $sin\theta$ = 0.35.

Observed combination limits at 95% CL in the $(m_{a},tan\beta)$ plane under the assumption of $sin\theta$ = 0.7.

Expected combination limits at 95% CL in the $(m_{a},tan\beta)$ plane under the assumption of $sin\theta$ = 0.7.

1 sigma band of expected combination limits at 95% CL in the $(m_{a},tan\beta)$ plane under the assumption of $sin\theta$ = 0.7.

2 sigma band of expected combination limits at 95% CL in the $(m_{a},tan\beta)$ plane under the assumption of $sin\theta$ = 0.7.

Observed and expected combination limits at 95% CL as a function of $sin\theta$ for the low-mass hypothesis $m_{A}$ = 0.6 TeV, $m_{a}$ = 200 GeV.

Observed and expected combination limits at 95% CL as a function of $sin\theta$ for the high-mass hypothesis $m_{A}$ = 1 TeV, $m_{a}$ = 350 GeV.

Observed and expected combination limits at 95% CL as a function of $m_{\chi}$ following the parameter choices of $m_{A}$ = 1 TeV, $m_{a}$ = 400 GeV.

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

The post-fit distributions of $p_\textrm{T}^{\textrm{bal}}$ for the SR. Data are compared with background predictions, and six signal predictions covering a representative mediator mass and invisible fraction range are overlaid. The uncertainties include all systematic and statistical components. The last bin contains the overflow.

The post-fit distributions of |$\phi_{\textrm{max}} - \phi_{\textrm{min}}$| for the SR. Data are compared with background predictions, and six signal predictions covering a representative mediator mass and invisible fraction range are overlaid. The uncertainties include all systematic and statistical components. The last bin contains the overflow.

Observed 95% CL limits on the t-channel production cross-section for semi-visible jets for R$_{inv}$ = 0.1.

Observed 95% CL limits on the t-channel production cross-section for semi-visible jets for R$_{inv}$ = 0.2.

Observed 95% CL limits on the t-channel production cross-section for semi-visible jets for R$_{inv}$ = 0.4.

Observed 95% CL limits on the t-channel production cross-section for semi-visible jets for R$_{inv}$ = 0.6.

Observed 95% CL limits on the t-channel production cross-section for semi-visible jets for R$_{inv}$ = 0.8.

Observed 95% CL limits on the t-channel production cross-section for semi-visible jets for R$_{inv}$ = 0.9.

The LO theory prediction (for λ=1).

Cutflow table for 6 signals of different mediator masses and R$_{inv}$ fractions.

IRing2 for HT2>=500 GeV, NJets>=5

IRing2 for HT2>=1000 GeV, NJets>=2

IRing2 for HT2>=1000 GeV, NJets>=3

IRing2 for HT2>=1000 GeV, NJets>=4

IRing2 for HT2>=1000 GeV, NJets>=5

IRing2 for HT2>=1500 GeV, NJets>=2

IRing2 for HT2>=1500 GeV, NJets>=3

IRing2 for HT2>=1500 GeV, NJets>=4

IRing2 for HT2>=1500 GeV, NJets>=5

IRing128 for HT2>=500 GeV, NJets>=2

IRing128 for HT2>=500 GeV, NJets>=3

IRing128 for HT2>=500 GeV, NJets>=4

IRing128 for HT2>=500 GeV, NJets>=5

IRing128 for HT2>=1000 GeV, NJets>=2

IRing128 for HT2>=1000 GeV, NJets>=3

IRing128 for HT2>=1000 GeV, NJets>=4

IRing128 for HT2>=1000 GeV, NJets>=5

IRing128 for HT2>=1500 GeV, NJets>=2

IRing128 for HT2>=1500 GeV, NJets>=3

IRing128 for HT2>=1500 GeV, NJets>=4

IRing128 for HT2>=1500 GeV, NJets>=5

ICyl16 for HT2>=500 GeV, NJets>=2

ICyl16 for HT2>=500 GeV, NJets>=3

ICyl16 for HT2>=500 GeV, NJets>=4

ICyl16 for HT2>=500 GeV, NJets>=5

ICyl16 for HT2>=1000 GeV, NJets>=2

ICyl16 for HT2>=1000 GeV, NJets>=3

ICyl16 for HT2>=1000 GeV, NJets>=4

ICyl16 for HT2>=1000 GeV, NJets>=5

ICyl16 for HT2>=1500 GeV, NJets>=2

ICyl16 for HT2>=1500 GeV, NJets>=3

ICyl16 for HT2>=1500 GeV, NJets>=4

ICyl16 for HT2>=1500 GeV, NJets>=5

IRing2 covariance for HT2>=500 GeV, NJets>=2 (Table 1)

IRing2 covariance for HT2>=500 GeV, NJets>=3 (Table 2)

IRing2 covariance for HT2>=500 GeV, NJets>=4 (Table 3)

IRing2 covariance for HT2>=500 GeV, NJets>=5 (Table 4)

IRing2 covariance for HT2>=1000 GeV, NJets>=2 (Table 5)

IRing2 covariance for HT2>=1000 GeV, NJets>=3 (Table 6)

IRing2 covariance for HT2>=1000 GeV, NJets>=4 (Table 7)

IRing2 covariance for HT2>=1000 GeV, NJets>=5 (Table 8)

IRing2 covariance for HT2>=1500 GeV, NJets>=2 (Table 9)

IRing2 covariance for HT2>=1500 GeV, NJets>=3 (Table 10)

IRing2 covariance for HT2>=1500 GeV, NJets>=4 (Table 11)

IRing2 covariance for HT2>=1500 GeV, NJets>=5 (Table 12)

IRing128 covariance for HT2>=500 GeV, NJets>=2 (Table 13)

IRing128 covariance for HT2>=500 GeV, NJets>=3 (Table 14)

IRing128 covariance for HT2>=500 GeV, NJets>=4 (Table 15)

IRing128 covariance for HT2>=500 GeV, NJets>=5 (Table 16)

IRing128 covariance for HT2>=1000 GeV, NJets>=2 (Table 17)

IRing128 covariance for HT2>=1000 GeV, NJets>=3 (Table 18)

IRing128 covariance for HT2>=1000 GeV, NJets>=4 (Table 19)

IRing128 covariance for HT2>=1000 GeV, NJets>=5 (Table 20)

IRing128 covariance for HT2>=1500 GeV, NJets>=2 (Table 21)

IRing128 covariance for HT2>=1500 GeV, NJets>=3 (Table 22)

IRing128 covariance for HT2>=1500 GeV, NJets>=4 (Table 23)

IRing128 covariance for HT2>=1500 GeV, NJets>=5 (Table 24)

ICyl16 covariance for HT2>=500 GeV, NJets>=2 (Table 25)

ICyl16 covariance for HT2>=500 GeV, NJets>=3 (Table 26)

ICyl16 covariance for HT2>=500 GeV, NJets>=4 (Table 27)

ICyl16 covariance for HT2>=500 GeV, NJets>=5 (Table 28)

ICyl16 covariance for HT2>=1000 GeV, NJets>=2 (Table 29)

ICyl16 covariance for HT2>=1000 GeV, NJets>=3 (Table 30)

ICyl16 covariance for HT2>=1000 GeV, NJets>=4 (Table 31)

ICyl16 covariance for HT2>=1000 GeV, NJets>=5 (Table 32)

ICyl16 covariance for HT2>=1500 GeV, NJets>=2 (Table 33)

ICyl16 covariance for HT2>=1500 GeV, NJets>=3 (Table 34)

ICyl16 covariance for HT2>=1500 GeV, NJets>=4 (Table 35)

ICyl16 covariance for HT2>=1500 GeV, NJets>=5 (Table 36)

IRing2 covariance, complete

1-IRing128 covariance, complete

1-ICyl16 covariance, complete

Events in bins of the $E_{miss}^T$ in the SR.

Observed (solid line) and expected (dashed line) 95% CL upper limits on the cross-section of single plus non-resonant plus pair vector LQ production from the combination of the high b-jet $p_{T}$ signal region for the $\tau_\text{lep}\tau_\text{had}$ and $\tau_\text{had}\tau_\text{had}$ channels. [$U_1^{MIN}$ model ($\kappa$ = 1) with $\lambda$ = 1.0]

Observed (solid line) and expected (dashed line) 95% CL upper limits on the cross-section of single plus non-resonant plus pair vector LQ production from the combination of the high b-jet $p_{T}$ signal region for the $\tau_\text{lep}\tau_\text{had}$ and $\tau_\text{had}\tau_\text{had}$ channels. [$U_1^{MIN}$ model ($\kappa$ = 1) with $\lambda$ = 1.7]

Observed (solid line) and expected (dashed line) 95% CL upper limits on the cross-section of single plus non-resonant plus pair vector LQ production from the combination of the high b-jet $p_{T}$ signal region for the $\tau_\text{lep}\tau_\text{had}$ and $\tau_\text{had}\tau_\text{had}$ channels. [$U_1^{MIN}$ model ($\kappa$ = 1) with $\lambda$ = 2.5]

The two-dimensional 95% CL exclusion limits in the $\lambda$-$m_{LQ}$ plane for singly plus non-resonant produced vector LQ (green lines) and for the sum, referred as Total, of single plus non-resonant plus pair vector LQ production (blue lines). Regions to the left of the lines are excluded. The dotted area shows the preferred region where the chosen LQ model can explain observed B anomalies. [$\kappa$ = 0]

The two-dimensional 95% CL exclusion limits in the $\lambda$-$m_{LQ}$ plane for singly plus non-resonant produced vector LQ (green lines) and for the sum, referred as Total, of single plus non-resonant plus pair vector LQ production (blue lines). Regions to the left of the lines are excluded. The dotted area shows the preferred region where the chosen LQ model can explain observed B anomalies. [$\kappa$ = 1]

Observed (solid line) and expected (dashed line) 95% CL upper limits for $\lambda$ = 1.0 on the cross-section of single plus non-resonant plus pair $\widetilde{S_{1}}$ production hypotheses from the combination of the high b-jet $p_{T}$ category for the $\tau_\text{lep}\tau_\text{had}$ and $\tau_\text{had}\tau_\text{had}$ channels.

Observed (solid line) and expected (dashed line) 95% CL upper limits for $\lambda$ = 1.7 on the cross-section of single plus non-resonant plus pair $\widetilde{S_{1}}$ production hypotheses from the combination of the high b-jet $p_{T}$ category for the $\tau_\text{lep}\tau_\text{had}$ and $\tau_\text{had}\tau_\text{had}$ channels.

Observed (solid line) and expected (dashed line) 95% CL upper limits for $\lambda$ = 2.5 on the cross-section of single plus non-resonant plus pair $\widetilde{S_{1}}$ production hypotheses from the combination of the high b-jet $p_{T}$ category for the $\tau_\text{lep}\tau_\text{had}$ and $\tau_\text{had}\tau_\text{had}$ channels.

The two-dimensional 95% CL exclusion limits in the $\lambda$-$m_{LQ}$ plane for singly plus non-resonant produced $\widetilde{S_{1}}$ (green lines) and for the sum, referred as Total, of single plus non-resonant plus pair vector LQ production (blue lines). Regions to the left of the lines are excluded.

Observed (solid line) and expected (dashed line) 95% upper limits on the visible cross-section, $\sigma_{\mathrm{vis}}$, obtained from model-independent search by a signal-plus-background fit in the high and low b-jet $p_{T}$ category for the $\tau_\text{lep}\tau_\text{had}$ and $\tau_\text{had}\tau_\text{had}$ channels.

Observed (solid line) and expected (dashed line) 95% CL upper limits on the cross-section of singly plus non-resonant produced vector LQ signals from the combination of the high b-jet $p_{T}$ category for the $\tau_\text{lep}\tau_\text{had}$ and $\tau_\text{had}\tau_\text{had}$ channels. [$U_1^{YM}$ model ($\kappa$ = 0) with $\lambda$ = 1.0]

Observed (solid line) and expected (dashed line) 95% CL upper limits on the cross-section of singly plus non-resonant produced vector LQ signals from the combination of the high b-jet $p_{T}$ category for the $\tau_\text{lep}\tau_\text{had}$ and $\tau_\text{had}\tau_\text{had}$ channels. [$U_1^{YM}$ model ($\kappa$ = 0) with $\lambda$ = 1.7]

Observed (solid line) and expected (dashed line) 95% CL upper limits on the cross-section of singly plus non-resonant produced vector LQ signals from the combination of the high b-jet $p_{T}$ category for the $\tau_\text{lep}\tau_\text{had}$ and $\tau_\text{had}\tau_\text{had}$ channels. [$U_1^{YM}$ model ($\kappa$ = 0) with $\lambda$ = 2.5]

Observed (solid line) and expected (dashed line) 95% CL upper limits on the cross-section of singly plus non-resonant produced vector LQ signals from the combination of the high b-jet $p_{T}$ category for the $\tau_\text{lep}\tau_\text{had}$ and $\tau_\text{had}\tau_\text{had}$ channels. [$U_1^{MIN}$ model ($\kappa$ = 1) with $\lambda$ = 1.0]

Observed (solid line) and expected (dashed line) 95% CL upper limits on the cross-section of singly plus non-resonant produced vector LQ signals from the combination of the high b-jet $p_{T}$ category for the $\tau_\text{lep}\tau_\text{had}$ and $\tau_\text{had}\tau_\text{had}$ channels. [$U_1^{MIN}$ model ($\kappa$ = 1) with $\lambda$ = 1.7]

Observed (solid line) and expected (dashed line) 95% CL upper limits on the cross-section of singly plus non-resonant produced vector LQ signals from the combination of the high b-jet $p_{T}$ category for the $\tau_\text{lep}\tau_\text{had}$ and $\tau_\text{had}\tau_\text{had}$ channels. [$U_1^{MIN}$ model ($\kappa$ = 1) with $\lambda$ = 2.5]

Observed (solid line) and expected (dashed line) 95% CL upper limits for $\lambda$ = 1.0 on the cross-section of singly produced $\widetilde{S_{1}}$ signal hypotheses from the combination of the high b-jet $p_{T}$ category for the $\tau_\text{lep}\tau_\text{had}$ and $\tau_\text{had}\tau_\text{had}$ channels.

Observed (solid line) and expected (dashed line) 95% CL upper limits for $\lambda$ = 1.7 on the cross-section of singly produced $\widetilde{S_{1}}$ signal hypotheses from the combination of the high b-jet $p_{T}$ category for the $\tau_\text{lep}\tau_\text{had}$ and $\tau_\text{had}\tau_\text{had}$ channels.

Observed (solid line) and expected (dashed line) 95% CL upper limits for $\lambda$ = 2.5 on the cross-section of singly produced $\widetilde{S_{1}}$ signal hypotheses from the combination of the high b-jet $p_{T}$ category for the $\tau_\text{lep}\tau_\text{had}$ and $\tau_\text{had}\tau_\text{had}$ channels.

The signal acceptance times efficiency for single plus non-resonant produced LQ in the $\widetilde{S_{1}}$ model from the combination of the high and low b-jet $p_{T}$ signal regions for the $\tau_\text{lep}\tau_\text{had}$ and $\tau_\text{had}\tau_\text{had}$ channels. [$\lambda$ = 1.0]

The signal acceptance times efficiency for single plus non-resonant produced LQ in the $\widetilde{S_{1}}$ model from the combination of the high and low b-jet $p_{T}$ signal regions for the $\tau_\text{lep}\tau_\text{had}$ and $\tau_\text{had}\tau_\text{had}$ channels. [$\lambda$ = 1.7]

The signal acceptance times efficiency for single plus non-resonant produced LQ in the $\widetilde{S_{1}}$ model from the combination of the high and low b-jet $p_{T}$ signal regions for the $\tau_\text{lep}\tau_\text{had}$ and $\tau_\text{had}\tau_\text{had}$ channels. [$\lambda$ = 2.5]

The signal acceptance times efficiency for single plus non-resonant produced LQ in the $U_1^{MIN}$ model from the combination of the high and low b-jet $p_{T}$ signal regions for the $\tau_\text{lep}\tau_\text{had}$ and $\tau_\text{had}\tau_\text{had}$ channels. [$\lambda$ = 1.0]

The signal acceptance times efficiency for single plus non-resonant produced LQ in the $U_1^{MIN}$ model from the combination of the high and low b-jet $p_{T}$ signal regions for the $\tau_\text{lep}\tau_\text{had}$ and $\tau_\text{had}\tau_\text{had}$ channels. [$\lambda$ = 1.7]

The signal acceptance times efficiency for single plus non-resonant produced LQ in the $U_1^{MIN}$ model from the combination of the high and low b-jet $p_{T}$ signal regions for the $\tau_\text{lep}\tau_\text{had}$ and $\tau_\text{had}\tau_\text{had}$ channels. [$\lambda$ = 2.5]

The signal acceptance times efficiency for single plus non-resonant produced LQ in the $U_1^{YM}$ model from the combination of the high and low b-jet $p_{T}$ signal regions for the $\tau_\text{lep}\tau_\text{had}$ and $\tau_\text{had}\tau_\text{had}$ channels. [$\lambda$ = 1.0]

The signal acceptance times efficiency for single plus non-resonant produced LQ in the $U_1^{YM}$ model from the combination of the high and low b-jet $p_{T}$ signal regions for the $\tau_\text{lep}\tau_\text{had}$ and $\tau_\text{had}\tau_\text{had}$ channels. [$\lambda$ = 1.7]

The signal acceptance times efficiency for single plus non-resonant produced LQ in the $U_1^{YM}$ model from the combination of the high and low b-jet $p_{T}$ signal regions for the $\tau_\text{lep}\tau_\text{had}$ and $\tau_\text{had}\tau_\text{had}$ channels. [$\lambda$ = 2.5]

Observed (solid line) and expected (dashed line) 95% CL upper limits on the cross-section of singly produced vector LQ signals from the combination of the high and low b-jet $p_{T}$ category for the $\tau_\text{lep}\tau_\text{had}$ and $\tau_\text{had}\tau_\text{had}$ channels, neglecting the interference between the LQ non-resonant production processes and the SM processes. However, the interference between the non-resonant LQ production and SM diagrams, such as those from Z+jets, could be non-negligible for events in the low b-jet $p_{T}$ signal region. [$U_1^{YM}$ model ($\kappa$ = 0) with $\lambda$ = 1.0]

Observed (solid line) and expected (dashed line) 95% CL upper limits on the cross-section of singly produced vector LQ signals from the combination of the high and low b-jet $p_{T}$ category for the $\tau_\text{lep}\tau_\text{had}$ and $\tau_\text{had}\tau_\text{had}$ channels, neglecting the interference between the LQ non-resonant production processes and the SM processes. However, the interference between the non-resonant LQ production and SM diagrams, such as those from Z+jets, could be non-negligible for events in the low b-jet $p_{T}$ signal region. [$U_1^{YM}$ model ($\kappa$ = 0) with $\lambda$ = 1.7]

Observed (solid line) and expected (dashed line) 95% CL upper limits on the cross-section of singly produced vector LQ signals from the combination of the high and low b-jet $p_{T}$ category for the $\tau_\text{lep}\tau_\text{had}$ and $\tau_\text{had}\tau_\text{had}$ channels, neglecting the interference between the LQ non-resonant production processes and the SM processes. However, the interference between the non-resonant LQ production and SM diagrams, such as those from Z+jets, could be non-negligible for events in the low b-jet $p_{T}$ signal region. [$U_1^{YM}$ model ($\kappa$ = 0) with $\lambda$ = 2.5]

Observed (solid line) and expected (dashed line) 95% CL upper limits on the cross-section of singly produced vector LQ signals from the combination of the high and low b-jet $p_{T}$ category for the $\tau_\text{lep}\tau_\text{had}$ and $\tau_\text{had}\tau_\text{had}$ channels, neglecting the interference between the LQ non-resonant production processes and the SM processes. However, the interference between the non-resonant LQ production and SM diagrams, such as those from Z+jets, could be non-negligible for events in the low b-jet $p_{T}$ signal region. [$U_1^{MIN}$ model ($\kappa$ = 1) with $\lambda$ = 1.0]

Observed (solid line) and expected (dashed line) 95% CL upper limits on the cross-section of singly produced vector LQ signals from the combination of the high and low b-jet $p_{T}$ category for the $\tau_\text{lep}\tau_\text{had}$ and $\tau_\text{had}\tau_\text{had}$ channels, neglecting the interference between the LQ non-resonant production processes and the SM processes. However, the interference between the non-resonant LQ production and SM diagrams, such as those from Z+jets, could be non-negligible for events in the low b-jet $p_{T}$ signal region. [$U_1^{MIN}$ model ($\kappa$ = 1) with $\lambda$ = 1.7]

Observed (solid line) and expected (dashed line) 95% CL upper limits on the cross-section of singly produced vector LQ signals from the combination of the high and low b-jet $p_{T}$ category for the $\tau_\text{lep}\tau_\text{had}$ and $\tau_\text{had}\tau_\text{had}$ channels, neglecting the interference between the LQ non-resonant production processes and the SM processes. However, the interference between the non-resonant LQ production and SM diagrams, such as those from Z+jets, could be non-negligible for events in the low b-jet $p_{T}$ signal region. [$U_1^{MIN}$ model ($\kappa$ = 1) with $\lambda$ = 2.5]

Observed (solid line) and expected (dashed line) 95% CL upper limits on the cross-section of the single plus non-resonant plus pair production of vector LQ from the combination of the high and low b-jet $p_{T}$ category for the $\tau_\text{lep}\tau_\text{had}$ and $\tau_\text{had}\tau_\text{had}$ channels, neglecting the interference between the LQ non-resonant production processes and the SM processes. However, the interference between the non-resonant LQ production and SM diagrams, such as those from Z+jets, could be non-negligible for events in the low b-jet $p_{T}$ signal region. [$U_1^{YM}$ model ($\kappa$ = 0) with $\lambda$ = 1.0]

Observed (solid line) and expected (dashed line) 95% CL upper limits on the cross-section of the single plus non-resonant plus pair production of vector LQ from the combination of the high and low b-jet $p_{T}$ category for the $\tau_\text{lep}\tau_\text{had}$ and $\tau_\text{had}\tau_\text{had}$ channels, neglecting the interference between the LQ non-resonant production processes and the SM processes. However, the interference between the non-resonant LQ production and SM diagrams, such as those from Z+jets, could be non-negligible for events in the low b-jet $p_{T}$ signal region. [$U_1^{YM}$ model ($\kappa$ = 0) with $\lambda$ = 1.7]

Observed (solid line) and expected (dashed line) 95% CL upper limits on the cross-section of the single plus non-resonant plus pair production of vector LQ from the combination of the high and low b-jet $p_{T}$ category for the $\tau_\text{lep}\tau_\text{had}$ and $\tau_\text{had}\tau_\text{had}$ channels, neglecting the interference between the LQ non-resonant production processes and the SM processes. However, the interference between the non-resonant LQ production and SM diagrams, such as those from Z+jets, could be non-negligible for events in the low b-jet $p_{T}$ signal region. [$U_1^{YM}$ model ($\kappa$ = 0) with $\lambda$ = 2.5]

Observed (solid line) and expected (dashed line) 95% CL upper limits on the cross-section of the single plus non-resonant plus pair production of vector LQ from the combination of the high and low b-jet $p_{T}$ category for the $\tau_\text{lep}\tau_\text{had}$ and $\tau_\text{had}\tau_\text{had}$ channels, neglecting the interference between the LQ non-resonant production processes and the SM processes. However, the interference between the non-resonant LQ production and SM diagrams, such as those from Z+jets, could be non-negligible for events in the low b-jet $p_{T}$ signal region. [$U_1^{MIN}$ model ($\kappa$ = 1) with $\lambda$ = 1.0]

Observed (solid line) and expected (dashed line) 95% CL upper limits on the cross-section of the single plus non-resonant plus pair production of vector LQ from the combination of the high and low b-jet $p_{T}$ category for the $\tau_\text{lep}\tau_\text{had}$ and $\tau_\text{had}\tau_\text{had}$ channels, neglecting the interference between the LQ non-resonant production processes and the SM processes. However, the interference between the non-resonant LQ production and SM diagrams, such as those from Z+jets, could be non-negligible for events in the low b-jet $p_{T}$ signal region. [$U_1^{MIN}$ model ($\kappa$ = 1) with $\lambda$ = 1.7]

Observed (solid line) and expected (dashed line) 95% CL upper limits on the cross-section of the single plus non-resonant plus pair production of vector LQ from the combination of the high and low b-jet $p_{T}$ category for the $\tau_\text{lep}\tau_\text{had}$ and $\tau_\text{had}\tau_\text{had}$ channels, neglecting the interference between the LQ non-resonant production processes and the SM processes. However, the interference between the non-resonant LQ production and SM diagrams, such as those from Z+jets, could be non-negligible for events in the low b-jet $p_{T}$ signal region. [$U_1^{MIN}$ model ($\kappa$ = 1) with $\lambda$ = 2.5]

The two-dimensional 95% CL exclusion limits in the $\lambda$-$m_{LQ}$ plane for singly plus non-resonant produced vector LQ (green lines) and for the sum, referred as Total, of single plus non-resonant plus pair vector LQ production (blue lines).Regions to the left of the lines are excluded. The dotted area shows the preferred region where the chosen LQ model can explain observed B anomalies. Results are extracted from the combination of the high and low b-jet $p_{T}$ signal region for the $\tau_\text{lep}\tau_\text{had}$ and $\tau_\text{had}\tau_\text{had}$ channels, neglecting the interference between the LQ non-resonant production processes and the SM processes. However, the interference between the non-resonant LQ production and SM diagrams, such as those from Z+jets, could be non-negligible for events in the low b-jet $p_{T}$ signal region. [$\kappa$ = 0]

The two-dimensional 95% CL exclusion limits in the $\lambda$-$m_{LQ}$ plane for singly plus non-resonant produced vector LQ (green lines) and for the sum, referred as Total, of single plus non-resonant plus pair vector LQ production (blue lines).Regions to the left of the lines are excluded. The dotted area shows the preferred region where the chosen LQ model can explain observed B anomalies. Results are extracted from the combination of the high and low b-jet $p_{T}$ signal region for the $\tau_\text{lep}\tau_\text{had}$ and $\tau_\text{had}\tau_\text{had}$ channels, neglecting the interference between the LQ non-resonant production processes and the SM processes. However, the interference between the non-resonant LQ production and SM diagrams, such as those from Z+jets, could be non-negligible for events in the low b-jet $p_{T}$ signal region. [$\kappa$ = 1]

Observed (solid line) and expected (dashed line) 95% CL upper limits for $\lambda$ = 1.0 on the cross-section of singly-produced $\widetilde{S_{1}}$ signals from the combination of the high and low b-jet $p_{T}$ categories for the $\tau_\text{lep}\tau_\text{had}$ and $\tau_\text{had}\tau_\text{had}$ channels, neglecting the interference between the LQ non-resonant production processes and the SM processes. However, the interference between the non-resonant LQ production and SM diagrams, such as those from Z+jets, could be non-negligible for events in the low b-jet $p_{T}$ signal region.

Observed (solid line) and expected (dashed line) 95% CL upper limits for $\lambda$ = 1.7 on the cross-section of singly-produced $\widetilde{S_{1}}$ signals from the combination of the high and low b-jet $p_{T}$ categories for the $\tau_\text{lep}\tau_\text{had}$ and $\tau_\text{had}\tau_\text{had}$ channels, neglecting the interference between the LQ non-resonant production processes and the SM processes. However, the interference between the non-resonant LQ production and SM diagrams, such as those from Z+jets, could be non-negligible for events in the low b-jet $p_{T}$ signal region.

Observed (solid line) and expected (dashed line) 95% CL upper limits for $\lambda$ = 2.5 on the cross-section of singly-produced $\widetilde{S_{1}}$ signals from the combination of the high and low b-jet $p_{T}$ categories for the $\tau_\text{lep}\tau_\text{had}$ and $\tau_\text{had}\tau_\text{had}$ channels, neglecting the interference between the LQ non-resonant production processes and the SM processes. However, the interference between the non-resonant LQ production and SM diagrams, such as those from Z+jets, could be non-negligible for events in the low b-jet $p_{T}$ signal region.

Observed (solid line) and expected (dashed line) 95% CL upper limits for $\lambda$ = 1.0 on the cross-section of single plus non-resonant plus pair $\widetilde{S_{1}}$ production signal from the combination of the high and low b-jet $p_{T}$ categories for the $\tau_\text{lep}\tau_\text{had}$ and $\tau_\text{had}\tau_\text{had}$ channels, neglecting the interference between the LQ non-resonant production processes and the SM processes. However, the interference between the non-resonant LQ production and SM diagrams, such as those from Z+jets, could be non-negligible for events in the low b-jet $p_{T}$ signal region.

Observed (solid line) and expected (dashed line) 95% CL upper limits for $\lambda$ = 1.7 on the cross-section of single plus non-resonant plus pair $\widetilde{S_{1}}$ production signal from the combination of the high and low b-jet $p_{T}$ categories for the $\tau_\text{lep}\tau_\text{had}$ and $\tau_\text{had}\tau_\text{had}$ channels, neglecting the interference between the LQ non-resonant production processes and the SM processes. However, the interference between the non-resonant LQ production and SM diagrams, such as those from Z+jets, could be non-negligible for events in the low b-jet $p_{T}$ signal region.

Observed (solid line) and expected (dashed line) 95% CL upper limits for $\lambda$ = 2.5 on the cross-section of single plus non-resonant plus pair $\widetilde{S_{1}}$ production signal from the combination of the high and low b-jet $p_{T}$ categories for the $\tau_\text{lep}\tau_\text{had}$ and $\tau_\text{had}\tau_\text{had}$ channels, neglecting the interference between the LQ non-resonant production processes and the SM processes. However, the interference between the non-resonant LQ production and SM diagrams, such as those from Z+jets, could be non-negligible for events in the low b-jet $p_{T}$ signal region.

The two-dimensional 95% CL exclusion limits in the $\lambda$-$m_{LQ}$ plane for singly plus non-resonant produced $\widetilde{S_{1}}$ (green lines) and for the sum, referred as Total, of single plus non-resonant plus pair $\widetilde{S_{1}}$ production (blue lines). Regions to the left of the lines are excluded. Results are extracted from the combination of the high and low b-jet $p_{T}$ signal region for the $\tau_\text{lep}\tau_\text{had}$ and $\tau_\text{had}\tau_\text{had}$ channels, neglecting the interference between the LQ non-resonant production processes and the SM processes. However, the interference between the non-resonant LQ production and SM diagrams, such as those from Z+jets, could be non-negligible for events in the low b-jet $p_{T}$ signal region.

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

The expected minus one standard deviation exclusion limit at 95% CL for the clockwork gravity model projected in the $k–M_{5}$ parameter space for the $ee$ channel for the case with mass thresholds.

The expected plus two standard deviation exclusion limit at 95% CL for the clockwork gravity model projected in the $k–M_{5}$ parameter space for the $ee$ channel for the case with mass thresholds.

The expected minus two standard deviation exclusion limit at 95% CL for the clockwork gravity model projected in the $k–M_{5}$ parameter space for the $ee$ channel for the case with mass thresholds.

The observed exclusion limit at 95% CL for the clockwork gravity model projected in the $k–M_{5}$ parameter space for the $\gamma\gamma$ channel for the case with mass thresholds.

The median expected exclusion limit at 95% CL for the clockwork gravity model projected in the $k–M_{5}$ parameter space for the $\gamma\gamma$ channel for the case with mass thresholds.

The expected plus one standard deviation exclusion limit at 95% CL for the clockwork gravity model projected in the $k–M_{5}$ parameter space for the $\gamma\gamma$ channel for the case with mass thresholds.

The expected minus one standard deviation exclusion limit at 95% CL for the clockwork gravity model projected in the $k–M_{5}$ parameter space for the $\gamma\gamma$ channel for the case with mass thresholds.

The expected plus two standard deviation exclusion limit at 95% CL for the clockwork gravity model projected in the $k–M_{5}$ parameter space for the $\gamma\gamma$ channel for the case with mass thresholds.

The expected minus two standard deviation exclusion limit at 95% CL for the clockwork gravity model projected in the $k–M_{5}$ parameter space for the $\gamma\gamma$ channel for the case with mass thresholds.

The observed exclusion limit at 95% CL for the clockwork gravity model projected in the $k–M_{5}$ parameter space for the $ee$ channel for the case without mass thresholds.

The median expected exclusion limit at 95% CL for the clockwork gravity model projected in the $k–M_{5}$ parameter space for the $ee$ channel for the case without mass thresholds.

The expected plus one standard deviation exclusion limit at 95% CL for the clockwork gravity model projected in the $k–M_{5}$ parameter space for the $ee$ channel for the case without mass thresholds.

The expected minus one standard deviation exclusion limit at 95% CL for the clockwork gravity model projected in the $k–M_{5}$ parameter space for the $ee$ channel for the case without mass thresholds.

The expected plus two standard deviation exclusion limit at 95% CL for the clockwork gravity model projected in the $k–M_{5}$ parameter space for the $ee$ channel for the case without mass thresholds.

The expected minus two standard deviation exclusion limit at 95% CL for the clockwork gravity model projected in the $k–M_{5}$ parameter space for the $ee$ channel for the case without mass thresholds.

The observed exclusion limit at 95% CL for the clockwork gravity model projected in the $k–M_{5}$ parameter space for the $\gamma\gamma$ channel for the case without mass thresholds.

The median expected exclusion limit at 95% CL for the clockwork gravity model projected in the $k–M_{5}$ parameter space for the $\gamma\gamma$ channel for the case without mass thresholds.

The expected plus one standard deviation exclusion limit at 95% CL for the clockwork gravity model projected in the $k–M_{5}$ parameter space for the $\gamma\gamma$ channel for the case without mass thresholds.

The expected minus one standard deviation exclusion limit at 95% CL for the clockwork gravity model projected in the $k–M_{5}$ parameter space for the $\gamma\gamma$ channel for the case without mass thresholds.

The expected plus two standard deviation exclusion limit at 95% CL for the clockwork gravity model projected in the $k–M_{5}$ parameter space for the $\gamma\gamma$ channel for the case without mass thresholds.

The expected minus two standard deviation exclusion limit at 95% CL for the clockwork gravity model projected in the $k–M_{5}$ parameter space for the $\gamma\gamma$ channel for the case without mass thresholds.

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

Di-muon mass distribution for the $\mu\mu\mu$ signal region for data and expected background. The expected signal distribution for $m_a = 35$ GeV is shown assuming $\sigma(t\bar{t}a)\times \text{Br}(a\rightarrow\mu\mu) = 4$ fb. Rare background processes include $ZZ+$jets, $WWZ$, $WZZ$, $ZZZ$, and $t\bar{t}t\bar{t}$.

Fit of the Double Crystal Ball width as a function of the $a$-boson mass in the $e\mu\mu$ channel.

Expected and observed 95% CL upper limit on signal cross section as a function of $m_a$ for $t\bar{t}a$.

Expected and observed 95% CL upper limit on signal cross section as a function of $m_a$ for $H^{+}\rightarrow aW$, considering a top pair cross section of $\sigma(pp\rightarrow t\bar{t})= 833$ pb for the charged Higgs boson mass, $m_{H^{+}} = 120$ GeV.

Expected and observed 95% CL upper limit on signal cross section as a function of $m_a$ for $H^{+}\rightarrow aW$, considering a top pair cross section of $\sigma(pp\rightarrow t\bar{t})= 833$ pb for the charged Higgs boson mass, $m_{H^{+}} = 140$ GeV.

Expected and observed 95% CL upper limit on signal cross section as a function of $m_a$ for $H^{+}\rightarrow aW$, considering a top pair cross section of $\sigma(pp\rightarrow t\bar{t})= 833$ pb for the charged Higgs boson mass, $m_{H^{+}} = 160$ GeV.

Comparison of the expected background with data in a region dominated by non-prompt background. The region is defined exactly like the $e\mu\mu$ signal region but requiring the two selected muons to have the same electrical charge.

Scan of the observed p-value as a function of $m_a$ for the background-only hypothesis.

Scan of the observed p-value as a function of $m_a$ for the background-only hypothesis for $e\mu\mu$ channel.

Scan of the observed p-value as a function of $m_a$ for the background-only hypothesis for $\mu\mu\mu$ channel.

Acceptance times efficiency for the $t\bar{t}a$, $a\rightarrow \mu\mu$ signal as a function of the $m_a$ hypothesis.

Acceptance times efficiency for the $t \rightarrow bH^+$, $H^+ \rightarrow W^+ a$, $a\rightarrow \mu\mu$ signal as a function of the $m_a$ and $m_{H^+}$ hypotheses, for the $e\mu\mu$ signal region. The upper left region of the plot corresponds to signals for which on-shell decays of $H^+ \rightarrow Wa$ are not possible. The CP conjugate process is also included.

Acceptance times efficiency for the $t \rightarrow bH^+ $, $H^+ \rightarrow W^+ a$, $a\rightarrow \mu\mu$ signal as a function of the $m_a$ and $m_{H^+}$ hypotheses, for the $\mu\mu\mu$ signal region. The upper left region of the plot corresponds to signals for which on-shell decays of $H^+ \rightarrow Wa$ are not possible. The CP conjugate process is also included.

Expected lower limits on $\tan\beta$ for the $t\rightarrow H^{+}b, H^{+}\rightarrow aW, a\rightarrow\mu\mu$ signal in the context of the 2HDM as a function of $m_{H^{+}}$ and $m_a$.

Observed lower limits on $\tan\beta$ for the $t\rightarrow H^{+}b, H^{+}\rightarrow aW, a\rightarrow\mu\mu$ signal in the context of the 2HDM as a function of $m_{H^{+}}$ and $m_a$.

Cutflow for two signal mass points for $t\bar{t}a$, $m_{a} = 20$ GeV and $m_{a} = 60$ GeV, and one signal mass point for $H^+ \rightarrow Wa$ for the electron channel $e\mu\mu$. For the signal yields, a cross section times branching ratio of 1 fb is assumed. The yields are presented before the profile likelihood fit.

Cutflow for the dominant backgrounds estimated from simulation ($t\bar{t}Z$, $WZ$, $t\bar{t}H$) and data, for the electron channel $e\mu\mu$. The yields are presented before the profile likelihood fit.

Cutflow for the signal mass point $t\bar{t}a$, $m_{a} = 20$ GeV, as well as the dominant backgrounds estimated from simulation ($t\bar{t}Z$, $WZ$, $t\bar{t}H$) and data, for the corresponding signal mass hypothesis, for the muon channel $\mu\mu\mu$. For the signal yields, a cross section times branching ratio of 1 fb is assumed. The yields are presented before the profile likelihood fit.

Cutflow for the signal mass point $t\bar{t}a$, $m_{a} = 60$ GeV, as well as the dominant backgrounds estimated from simulation ($t\bar{t}Z$, $WZ$, $t\bar{t}H$) and data, for the corresponding signal mass hypothesis, for the muon channel $\mu\mu\mu$. For the signal yields, a cross section times branching ratio of 1 fb is assumed. The yields are presented before the profile likelihood fit.

Cutflow for $H^+ \rightarrow W a$ with $m_{H^{+}} = 120$ GeV, $m_{a} = 30$ GeV for the muon channel $\mu\mu\mu$. A cross section times branching ratio of 1 fb is assumed. The yields are presented before the profile likelihood fit.

Total yields in the $e\mu\mu$ and $\mu\mu\mu$ signal regions for the expected background and data after the profile likelihood fit for a $m_a = 30$ GeV signal mass hypothesis.

Average timing distributions for SR data and the estimated background as determined by the background-only fit, in each of the five exclusive $\rho$ categories. For comparison, the expected timing shapes for a few different signal models are superimposed, with each model labeled by the values of the $\tilde\chi^0_1$ mass and lifetime, as well as decay mode. To provide some indication of the variations in signal yield and shape, three signal models are shown for each of the $\tilde\chi^0_1$ decay modes, namely $\tilde\chi^0_1$ $\rightarrow$ $H \tilde G$ and $\tilde\chi^0_1$ $\rightarrow$ $Z \tilde G$. The models shown include a rather low $\tilde\chi^0_1$ mass value of 135 GeV for lifetimes of either 2 ns or 10 ns, and a higher $\tilde\chi^0_1$ mass value which is near the 95% CL exclusion limit for each decay mode for a lifetime of 2 ns. Each signal model is shown with the signal normalization corresponding to a BR value of unity for the decay mode in question.

Average timing distributions for SR data and the estimated background as determined by the background-only fit, in each of the five exclusive $\rho$ categories. For comparison, the expected timing shapes for a few different signal models are superimposed, with each model labeled by the values of the $\tilde\chi^0_1$ mass and lifetime, as well as decay mode. To provide some indication of the variations in signal yield and shape, three signal models are shown for each of the $\tilde\chi^0_1$ decay modes, namely $\tilde\chi^0_1$ $\rightarrow$ $H \tilde G$ and $\tilde\chi^0_1$ $\rightarrow$ $Z \tilde G$. The models shown include a rather low $\tilde\chi^0_1$ mass value of 135 GeV for lifetimes of either 2 ns or 10 ns, and a higher $\tilde\chi^0_1$ mass value which is near the 95% CL exclusion limit for each decay mode for a lifetime of 2 ns. Each signal model is shown with the signal normalization corresponding to a BR value of unity for the decay mode in question.

The 95% CL limits on $\sigma(pp \rightarrow \tilde\chi^0_1 \tilde\chi^0_1$) in fb as a function of $\tilde\chi^0_1$ mass (left) and $\tilde\chi^0_1$ lifetime (right), for the different decay modes of $\mathcal{B}$($\tilde\chi^0_1$ $\rightarrow H + \tilde{G}$) = 1 (top) and $\mathcal{B}$($\tilde\chi^0_1$ $\rightarrow Z +\tilde{G}$) = 1 (bottom). For the limits as a function of mass (lifetime), several signal models with varying lifetime (mass) are overlaid for comparison. Included are the theoretical expectations from higgsino production for each mass hypothesis, calculated from a GMSB SUSY model that assumes nearly degenerate $\tilde\chi^0_1$, $\tilde\chi^\pm_1$, and $\tilde\chi^0_2$.

The 95% CL limits on $\sigma(pp \rightarrow \tilde\chi^0_1 \tilde\chi^0_1$) in fb as a function of $\tilde\chi^0_1$ mass (left) and $\tilde\chi^0_1$ lifetime (right), for the different decay modes of $\mathcal{B}$($\tilde\chi^0_1$ $\rightarrow H + \tilde{G}$) = 1 (top) and $\mathcal{B}$($\tilde\chi^0_1$ $\rightarrow Z +\tilde{G}$) = 1 (bottom). For the limits as a function of mass (lifetime), several signal models with varying lifetime (mass) are overlaid for comparison. Included are the theoretical expectations from higgsino production for each mass hypothesis, calculated from a GMSB SUSY model that assumes nearly degenerate $\tilde\chi^0_1$, $\tilde\chi^\pm_1$, and $\tilde\chi^0_2$.

The 95% CL limits on $\sigma(pp \rightarrow \tilde\chi^0_1 \tilde\chi^0_1$) in fb as a function of $\tilde\chi^0_1$ mass (left) and $\tilde\chi^0_1$ lifetime (right), for the different decay modes of $\mathcal{B}$($\tilde\chi^0_1$ $\rightarrow H + \tilde{G}$) = 1 (top) and $\mathcal{B}$($\tilde\chi^0_1$ $\rightarrow Z +\tilde{G}$) = 1 (bottom). For the limits as a function of mass (lifetime), several signal models with varying lifetime (mass) are overlaid for comparison. Included are the theoretical expectations from higgsino production for each mass hypothesis, calculated from a GMSB SUSY model that assumes nearly degenerate $\tilde\chi^0_1$, $\tilde\chi^\pm_1$, and $\tilde\chi^0_2$.

The 95% CL limits on $\sigma(pp \rightarrow \tilde\chi^0_1 \tilde\chi^0_1$) in fb as a function of $\tilde\chi^0_1$ mass (left) and $\tilde\chi^0_1$ lifetime (right), for the different decay modes of $\mathcal{B}$($\tilde\chi^0_1$ $\rightarrow H + \tilde{G}$) = 1 (top) and $\mathcal{B}$($\tilde\chi^0_1$ $\rightarrow Z +\tilde{G}$) = 1 (bottom). For the limits as a function of mass (lifetime), several signal models with varying lifetime (mass) are overlaid for comparison. Included are the theoretical expectations from higgsino production for each mass hypothesis, calculated from a GMSB SUSY model that assumes nearly degenerate $\tilde\chi^0_1$, $\tilde\chi^\pm_1$, and $\tilde\chi^0_2$.

The 95% CL limits on $\sigma(pp \rightarrow \tilde\chi^0_1 \tilde\chi^0_1$) in fb as a function of $\tilde\chi^0_1$ mass (left) and $\tilde\chi^0_1$ lifetime (right), for the different decay modes of $\mathcal{B}$($\tilde\chi^0_1$ $\rightarrow H + \tilde{G}$) = 1 (top) and $\mathcal{B}$($\tilde\chi^0_1$ $\rightarrow Z +\tilde{G}$) = 1 (bottom). For the limits as a function of mass (lifetime), several signal models with varying lifetime (mass) are overlaid for comparison. Included are the theoretical expectations from higgsino production for each mass hypothesis, calculated from a GMSB SUSY model that assumes nearly degenerate $\tilde\chi^0_1$, $\tilde\chi^\pm_1$, and $\tilde\chi^0_2$.

The 95% CL limits on $\sigma(pp \rightarrow \tilde\chi^0_1 \tilde\chi^0_1$) in fb as a function of $\tilde\chi^0_1$ mass (left) and $\tilde\chi^0_1$ lifetime (right), for the different decay modes of $\mathcal{B}$($\tilde\chi^0_1$ $\rightarrow H + \tilde{G}$) = 1 (top) and $\mathcal{B}$($\tilde\chi^0_1$ $\rightarrow Z +\tilde{G}$) = 1 (bottom). For the limits as a function of mass (lifetime), several signal models with varying lifetime (mass) are overlaid for comparison. Included are the theoretical expectations from higgsino production for each mass hypothesis, calculated from a GMSB SUSY model that assumes nearly degenerate $\tilde\chi^0_1$, $\tilde\chi^\pm_1$, and $\tilde\chi^0_2$.

The 95% CL limits on $\sigma(pp \rightarrow \tilde\chi^0_1 \tilde\chi^0_1$) in fb as a function of $\tilde\chi^0_1$ mass (left) and $\tilde\chi^0_1$ lifetime (right), for the different decay modes of $\mathcal{B}$($\tilde\chi^0_1$ $\rightarrow H + \tilde{G}$) = 1 (top) and $\mathcal{B}$($\tilde\chi^0_1$ $\rightarrow Z +\tilde{G}$) = 1 (bottom). For the limits as a function of mass (lifetime), several signal models with varying lifetime (mass) are overlaid for comparison. Included are the theoretical expectations from higgsino production for each mass hypothesis, calculated from a GMSB SUSY model that assumes nearly degenerate $\tilde\chi^0_1$, $\tilde\chi^\pm_1$, and $\tilde\chi^0_2$.

The 95% CL limits on $\sigma(pp \rightarrow \tilde\chi^0_1 \tilde\chi^0_1$) in fb as a function of $\tilde\chi^0_1$ mass (left) and $\tilde\chi^0_1$ lifetime (right), for the different decay modes of $\mathcal{B}$($\tilde\chi^0_1$ $\rightarrow H + \tilde{G}$) = 1 (top) and $\mathcal{B}$($\tilde\chi^0_1$ $\rightarrow Z +\tilde{G}$) = 1 (bottom). For the limits as a function of mass (lifetime), several signal models with varying lifetime (mass) are overlaid for comparison. Included are the theoretical expectations from higgsino production for each mass hypothesis, calculated from a GMSB SUSY model that assumes nearly degenerate $\tilde\chi^0_1$, $\tilde\chi^\pm_1$, and $\tilde\chi^0_2$.

The 95% CL limits on $\sigma(pp \rightarrow \tilde\chi^0_1 \tilde\chi^0_1$) in fb as a function of $\tilde\chi^0_1$ mass (left) and $\tilde\chi^0_1$ lifetime (right), for the different decay modes of $\mathcal{B}$($\tilde\chi^0_1$ $\rightarrow H + \tilde{G}$) = 1 (top) and $\mathcal{B}$($\tilde\chi^0_1$ $\rightarrow Z +\tilde{G}$) = 1 (bottom). For the limits as a function of mass (lifetime), several signal models with varying lifetime (mass) are overlaid for comparison. Included are the theoretical expectations from higgsino production for each mass hypothesis, calculated from a GMSB SUSY model that assumes nearly degenerate $\tilde\chi^0_1$, $\tilde\chi^\pm_1$, and $\tilde\chi^0_2$.

The 95% CL limits on $\sigma(pp \rightarrow \tilde\chi^0_1 \tilde\chi^0_1$) in fb as a function of $\tilde\chi^0_1$ mass (left) and $\tilde\chi^0_1$ lifetime (right), for the different decay modes of $\mathcal{B}$($\tilde\chi^0_1$ $\rightarrow H + \tilde{G}$) = 1 (top) and $\mathcal{B}$($\tilde\chi^0_1$ $\rightarrow Z +\tilde{G}$) = 1 (bottom). For the limits as a function of mass (lifetime), several signal models with varying lifetime (mass) are overlaid for comparison. Included are the theoretical expectations from higgsino production for each mass hypothesis, calculated from a GMSB SUSY model that assumes nearly degenerate $\tilde\chi^0_1$, $\tilde\chi^\pm_1$, and $\tilde\chi^0_2$.

The 95% CL limits on $\sigma(pp \rightarrow \tilde\chi^0_1 \tilde\chi^0_1$) in fb as a function of $\tilde\chi^0_1$ mass (left) and $\tilde\chi^0_1$ lifetime (right), for the different decay modes of $\mathcal{B}$($\tilde\chi^0_1$ $\rightarrow H + \tilde{G}$) = 1 (top) and $\mathcal{B}$($\tilde\chi^0_1$ $\rightarrow Z +\tilde{G}$) = 1 (bottom). For the limits as a function of mass (lifetime), several signal models with varying lifetime (mass) are overlaid for comparison. Included are the theoretical expectations from higgsino production for each mass hypothesis, calculated from a GMSB SUSY model that assumes nearly degenerate $\tilde\chi^0_1$, $\tilde\chi^\pm_1$, and $\tilde\chi^0_2$.

The 95% CL limits on $\sigma(pp \rightarrow \tilde\chi^0_1 \tilde\chi^0_1$) in fb as a function of $\tilde\chi^0_1$ mass (left) and $\tilde\chi^0_1$ lifetime (right), for the different decay modes of $\mathcal{B}$($\tilde\chi^0_1$ $\rightarrow H + \tilde{G}$) = 1 (top) and $\mathcal{B}$($\tilde\chi^0_1$ $\rightarrow Z +\tilde{G}$) = 1 (bottom). For the limits as a function of mass (lifetime), several signal models with varying lifetime (mass) are overlaid for comparison. Included are the theoretical expectations from higgsino production for each mass hypothesis, calculated from a GMSB SUSY model that assumes nearly degenerate $\tilde\chi^0_1$, $\tilde\chi^\pm_1$, and $\tilde\chi^0_2$.

The 95% CL limits on $\sigma(pp \rightarrow \tilde\chi^0_1 \tilde\chi^0_1$) in fb as a function of $\tilde\chi^0_1$ mass (left) and $\tilde\chi^0_1$ lifetime (right), for the different decay modes of $\mathcal{B}$($\tilde\chi^0_1$ $\rightarrow H + \tilde{G}$) = 1 (top) and $\mathcal{B}$($\tilde\chi^0_1$ $\rightarrow Z +\tilde{G}$) = 1 (bottom). For the limits as a function of mass (lifetime), several signal models with varying lifetime (mass) are overlaid for comparison. Included are the theoretical expectations from higgsino production for each mass hypothesis, calculated from a GMSB SUSY model that assumes nearly degenerate $\tilde\chi^0_1$, $\tilde\chi^\pm_1$, and $\tilde\chi^0_2$.

The 95% CL limits on $\sigma(pp \rightarrow \tilde\chi^0_1 \tilde\chi^0_1$) in fb as a function of $\tilde\chi^0_1$ mass (left) and $\tilde\chi^0_1$ lifetime (right), for the different decay modes of $\mathcal{B}$($\tilde\chi^0_1$ $\rightarrow H + \tilde{G}$) = 1 (top) and $\mathcal{B}$($\tilde\chi^0_1$ $\rightarrow Z +\tilde{G}$) = 1 (bottom). For the limits as a function of mass (lifetime), several signal models with varying lifetime (mass) are overlaid for comparison. Included are the theoretical expectations from higgsino production for each mass hypothesis, calculated from a GMSB SUSY model that assumes nearly degenerate $\tilde\chi^0_1$, $\tilde\chi^\pm_1$, and $\tilde\chi^0_2$.

The 95% CL limits on $\sigma(pp \rightarrow \tilde\chi^0_1 \tilde\chi^0_1$) in fb as a function of $\tilde\chi^0_1$ mass (left) and $\tilde\chi^0_1$ lifetime (right), for the different decay modes of $\mathcal{B}$($\tilde\chi^0_1$ $\rightarrow H + \tilde{G}$) = 1 (top) and $\mathcal{B}$($\tilde\chi^0_1$ $\rightarrow Z +\tilde{G}$) = 1 (bottom). For the limits as a function of mass (lifetime), several signal models with varying lifetime (mass) are overlaid for comparison. Included are the theoretical expectations from higgsino production for each mass hypothesis, calculated from a GMSB SUSY model that assumes nearly degenerate $\tilde\chi^0_1$, $\tilde\chi^\pm_1$, and $\tilde\chi^0_2$.

The 95% CL limits on $\sigma(pp \rightarrow \tilde\chi^0_1 \tilde\chi^0_1$) in fb as a function of $\tilde\chi^0_1$ mass (left) and $\tilde\chi^0_1$ lifetime (right), for the different decay modes of $\mathcal{B}$($\tilde\chi^0_1$ $\rightarrow H + \tilde{G}$) = 1 (top) and $\mathcal{B}$($\tilde\chi^0_1$ $\rightarrow Z +\tilde{G}$) = 1 (bottom). For the limits as a function of mass (lifetime), several signal models with varying lifetime (mass) are overlaid for comparison. Included are the theoretical expectations from higgsino production for each mass hypothesis, calculated from a GMSB SUSY model that assumes nearly degenerate $\tilde\chi^0_1$, $\tilde\chi^\pm_1$, and $\tilde\chi^0_2$.

The 95% CL limits on $\sigma(pp \rightarrow \tilde\chi^0_1 \tilde\chi^0_1$) in fb as a function of $\tilde\chi^0_1$ branching ratio to the SM Higgs boson, where the assumed cross-section is for higgsino production, and $\mathcal{B}$($\tilde\chi^0_1$ $\rightarrow Z +\tilde{G}$) = 1 - $\mathcal{B}$($\tilde\chi^0_1$ $\rightarrow H + \tilde{G}$). Several signal hypotheses are overlaid that are labelled by the $\tilde\chi^0_1$ mass, all with a fixed $\tilde\chi^0_1$ lifetime of 2 ns.

The 95% CL limits on $\sigma(pp \rightarrow \tilde\chi^0_1 \tilde\chi^0_1$) in fb as a function of $\tilde\chi^0_1$ branching ratio to the SM Higgs boson, where the assumed cross-section is for higgsino production, and $\mathcal{B}$($\tilde\chi^0_1$ $\rightarrow Z +\tilde{G}$) = 1 - $\mathcal{B}$($\tilde\chi^0_1$ $\rightarrow H + \tilde{G}$). Several signal hypotheses are overlaid that are labelled by the $\tilde\chi^0_1$ mass, all with a fixed $\tilde\chi^0_1$ lifetime of 2 ns.

The 95% CL limits on $\sigma(pp \rightarrow \tilde\chi^0_1 \tilde\chi^0_1$) in fb as a function of $\tilde\chi^0_1$ branching ratio to the SM Higgs boson, where the assumed cross-section is for higgsino production, and $\mathcal{B}$($\tilde\chi^0_1$ $\rightarrow Z +\tilde{G}$) = 1 - $\mathcal{B}$($\tilde\chi^0_1$ $\rightarrow H + \tilde{G}$). Several signal hypotheses are overlaid that are labelled by the $\tilde\chi^0_1$ mass, all with a fixed $\tilde\chi^0_1$ lifetime of 2 ns.

The 95% CL limits on $\sigma(pp \rightarrow \tilde\chi^0_1 \tilde\chi^0_1$) in fb as a function of $\tilde\chi^0_1$ branching ratio to the SM Higgs boson, where the assumed cross-section is for higgsino production, and $\mathcal{B}$($\tilde\chi^0_1$ $\rightarrow Z +\tilde{G}$) = 1 - $\mathcal{B}$($\tilde\chi^0_1$ $\rightarrow H + \tilde{G}$). Several signal hypotheses are overlaid that are labelled by the $\tilde\chi^0_1$ mass, all with a fixed $\tilde\chi^0_1$ lifetime of 2 ns.

The 95% CL exclusion limits on the target signal hypothesis, for $\tilde\chi^0_1$ lifetime in ns as a function of $\tilde\chi^0_1$ mass in GeV. The overlaid curves correspond to different decay hypotheses, where the assumed cross-section is for higgsino production, and the $\tilde\chi^0_1$ decays to $H + \tilde{G}$ or $Z + \tilde{G}$ such that $\mathcal{B}(H + \tilde{G}) + \mathcal{B}(Z + \tilde{G})$ = 100%. The curve shown in red represents the decay hypothesis where the $\tilde\chi^0_1$ decays to $Z + \tilde{G}$ with 100% branching ratio. The curve shown in blue represents the decay hypothesis where the $\tilde\chi^0_1$ decays to $H + \tilde{G}$ with 100% branching ratio.

The 95% CL exclusion limits on the target signal hypothesis, for $\tilde\chi^0_1$ lifetime in ns as a function of $\tilde\chi^0_1$ mass in GeV. The overlaid curves correspond to different decay hypotheses, where the assumed cross-section is for higgsino production, and the $\tilde\chi^0_1$ decays to $H + \tilde{G}$ or $Z + \tilde{G}$ such that $\mathcal{B}(H + \tilde{G}) + \mathcal{B}(Z + \tilde{G})$ = 100%. The curve shown in red represents the decay hypothesis where the $\tilde\chi^0_1$ decays to $Z + \tilde{G}$ with 100% branching ratio. The curve shown in blue represents the decay hypothesis where the $\tilde\chi^0_1$ decays to $H + \tilde{G}$ with 100% branching ratio.

The 95% CL exclusion limits on the target signal hypothesis, for $\tilde\chi^0_1$ lifetime in ns as a function of $\tilde\chi^0_1$ mass in GeV. The overlaid curves correspond to different decay hypotheses, where the assumed cross-section is for higgsino production, and the $\tilde\chi^0_1$ decays to $H + \tilde{G}$ or $Z + \tilde{G}$ such that $\mathcal{B}(H + \tilde{G}) + \mathcal{B}(Z + \tilde{G})$ = 100%. The curve shown in red represents the decay hypothesis where the $\tilde\chi^0_1$ decays to $Z + \tilde{G}$ with 100% branching ratio. The curve shown in blue represents the decay hypothesis where the $\tilde\chi^0_1$ decays to $H + \tilde{G}$ with 100% branching ratio.

The 95% CL exclusion limits on the target signal hypothesis, for $\tilde\chi^0_1$ lifetime in ns as a function of $\tilde\chi^0_1$ mass in GeV. The overlaid curves correspond to different decay hypotheses, where the assumed cross-section is for higgsino production, and the $\tilde\chi^0_1$ decays to $H + \tilde{G}$ or $Z + \tilde{G}$ such that $\mathcal{B}(H + \tilde{G}) + \mathcal{B}(Z + \tilde{G})$ = 100%. The curve shown in red represents the decay hypothesis where the $\tilde\chi^0_1$ decays to $Z + \tilde{G}$ with 100% branching ratio. The curve shown in blue represents the decay hypothesis where the $\tilde\chi^0_1$ decays to $H + \tilde{G}$ with 100% branching ratio.

The 95% CL exclusion limits on the target signal hypothesis, for $\tilde\chi^0_1$ lifetime in ns as a function of $\tilde\chi^0_1$ mass in GeV. The overlaid curves correspond to different decay hypotheses, where the assumed cross-section is for higgsino production, and the $\tilde\chi^0_1$ decays to $H + \tilde{G}$ or $Z + \tilde{G}$ such that $\mathcal{B}(H + \tilde{G}) + \mathcal{B}(Z + \tilde{G})$ = 100%. The curve shown in red represents the decay hypothesis where the $\tilde\chi^0_1$ decays to $Z + \tilde{G}$ with 100% branching ratio. The curve shown in blue represents the decay hypothesis where the $\tilde\chi^0_1$ decays to $H + \tilde{G}$ with 100% branching ratio.

The 95% CL exclusion limits on the target signal hypothesis, for $\tilde\chi^0_1$ lifetime in ns as a function of $\tilde\chi^0_1$ mass in GeV. The overlaid curves correspond to different decay hypotheses, where the assumed cross-section is for higgsino production, and the $\tilde\chi^0_1$ decays to $H + \tilde{G}$ or $Z + \tilde{G}$ such that $\mathcal{B}(H + \tilde{G}) + \mathcal{B}(Z + \tilde{G})$ = 100%. The curve shown in red represents the decay hypothesis where the $\tilde\chi^0_1$ decays to $Z + \tilde{G}$ with 100% branching ratio. The curve shown in blue represents the decay hypothesis where the $\tilde\chi^0_1$ decays to $H + \tilde{G}$ with 100% branching ratio.

The 95% CL exclusion limits on the target signal hypothesis, for $\tilde\chi^0_1$ lifetime in ns as a function of $\tilde\chi^0_1$ mass in GeV. The overlaid curves correspond to different decay hypotheses, where the assumed cross-section is for higgsino production, and the $\tilde\chi^0_1$ decays to $H + \tilde{G}$ or $Z + \tilde{G}$ such that $\mathcal{B}(H + \tilde{G}) + \mathcal{B}(Z + \tilde{G})$ = 100%. The curve shown in red represents the decay hypothesis where the $\tilde\chi^0_1$ decays to $Z + \tilde{G}$ with 100% branching ratio. The curve shown in blue represents the decay hypothesis where the $\tilde\chi^0_1$ decays to $H + \tilde{G}$ with 100% branching ratio.

The 95% CL exclusion limits on the target signal hypothesis, for $\tilde\chi^0_1$ lifetime in ns as a function of $\tilde\chi^0_1$ mass in GeV. The overlaid curves correspond to different decay hypotheses, where the assumed cross-section is for higgsino production, and the $\tilde\chi^0_1$ decays to $H + \tilde{G}$ or $Z + \tilde{G}$ such that $\mathcal{B}(H + \tilde{G}) + \mathcal{B}(Z + \tilde{G})$ = 100%. The curve shown in red represents the decay hypothesis where the $\tilde\chi^0_1$ decays to $Z + \tilde{G}$ with 100% branching ratio. The curve shown in blue represents the decay hypothesis where the $\tilde\chi^0_1$ decays to $H + \tilde{G}$ with 100% branching ratio.

The 95% CL exclusion limits on the target signal hypothesis, for $\tilde\chi^0_1$ lifetime in ns as a function of $\tilde\chi^0_1$ mass in GeV. The overlaid curves correspond to different decay hypotheses, where the assumed cross-section is for higgsino production, and the $\tilde\chi^0_1$ decays to $H + \tilde{G}$ or $Z + \tilde{G}$ such that $\mathcal{B}(H + \tilde{G}) + \mathcal{B}(Z + \tilde{G})$ = 100%. The curve shown in red represents the decay hypothesis where the $\tilde\chi^0_1$ decays to $Z + \tilde{G}$ with 100% branching ratio. The curve shown in blue represents the decay hypothesis where the $\tilde\chi^0_1$ decays to $H + \tilde{G}$ with 100% branching ratio.

The 95% CL exclusion limits on the target signal hypothesis, for $\tilde\chi^0_1$ lifetime in ns as a function of $\tilde\chi^0_1$ mass in GeV. The overlaid curves correspond to different decay hypotheses, where the assumed cross-section is for higgsino production, and the $\tilde\chi^0_1$ decays to $H + \tilde{G}$ or $Z + \tilde{G}$ such that $\mathcal{B}(H + \tilde{G}) + \mathcal{B}(Z + \tilde{G})$ = 100%. The curve shown in red represents the decay hypothesis where the $\tilde\chi^0_1$ decays to $Z + \tilde{G}$ with 100% branching ratio. The curve shown in blue represents the decay hypothesis where the $\tilde\chi^0_1$ decays to $H + \tilde{G}$ with 100% branching ratio.

The 95% CL exclusion limits on the target signal hypothesis, for $\tilde\chi^0_1$ lifetime in ns as a function of $\tilde\chi^0_1$ mass in GeV. The overlaid curves correspond to different decay hypotheses, where the assumed cross-section is for higgsino production, and the $\tilde\chi^0_1$ decays to $H + \tilde{G}$ or $Z + \tilde{G}$ such that $\mathcal{B}(H + \tilde{G}) + \mathcal{B}(Z + \tilde{G})$ = 100%. The curve shown in red represents the decay hypothesis where the $\tilde\chi^0_1$ decays to $Z + \tilde{G}$ with 100% branching ratio. The curve shown in blue represents the decay hypothesis where the $\tilde\chi^0_1$ decays to $H + \tilde{G}$ with 100% branching ratio.

The 95% CL exclusion limits on the target signal hypothesis, for $\tilde\chi^0_1$ lifetime in ns as a function of $\tilde\chi^0_1$ mass in GeV. The overlaid curves correspond to different decay hypotheses, where the assumed cross-section is for higgsino production, and the $\tilde\chi^0_1$ decays to $H + \tilde{G}$ or $Z + \tilde{G}$ such that $\mathcal{B}(H + \tilde{G}) + \mathcal{B}(Z + \tilde{G})$ = 100%. The curve shown in red represents the decay hypothesis where the $\tilde\chi^0_1$ decays to $Z + \tilde{G}$ with 100% branching ratio. The curve shown in blue represents the decay hypothesis where the $\tilde\chi^0_1$ decays to $H + \tilde{G}$ with 100% branching ratio.

Cutflow for an example higgsino signal with mass 225 GeV and lifetime 10 ns, in the H decay mode. Acceptance is defined at truth level, and efficiency compares the events passing at reconstruction level with respect to truth.

Cutflow for an example higgsino signal with mass 225 GeV and lifetime 10 ns, in the Z decay mode. Acceptance is defined at truth level, and efficiency compares the events passing at reconstruction level with respect to truth.

Cutflow for an example higgsino signal with mass 225 GeV and lifetime 2 ns, in the H decay mode. Acceptance is defined at truth level, and efficiency compares the events passing at reconstruction level with respect to truth.

Cutflow for an example higgsino signal with mass 225 GeV and lifetime 2 ns, in the Z decay mode. Acceptance is defined at truth level, and efficiency compares the events passing at reconstruction level with respect to truth.

Acceptance across the H decay mode signal grid, calculated using truth information. The selection applied corresponds to the model-independent signal region (i.e. the standard SR with $t_{\text{avg}$ > 0.9 ns).

Acceptance across the Z decay mode signal grid, calculated using truth information. The selection applied corresponds to the model-independent signal region (i.e. the standard SR with $t_{\text{avg}$ > 0.9 ns).

Efficiency across the H decay mode signal grid, calculated using reco information. The selection applied corresponds to the model-independent signal region (i.e. the standard SR with $t_{\text{avg}$ > 0.9 ns). Here, the numerator is the signal yield passing the reco selection and the denominator is the signal yield passing the truth selection.

Efficiency across the Z decay mode signal grid, calculated using reco information. The selection applied corresponds to the model-independent signal region (i.e. the standard SR with $t_{\text{avg}$ > 0.9 ns). Here, the numerator is the signal yield passing the reco selection and the denominator is the signal yield passing the truth selection.

Total efficiencies for the selection of signal events after categorization and application of the tighter photon $E_{\mathrm{T}}^{\gamma}$ selection used to optimize the signal significance spin-1 $q\bar{q'}\to X^{\pm} \to W^{\pm}\gamma$. In addition to the total efficiency, contributions to the signal selection from each of the separate event categories are shown. The efficiencies calculated from MC samples with $W/Z$ hadronic decays are shown as the points on each curve. The line presents interpolated results.

Total efficiencies for the selection of signal events after categorization and application of the tighter photon $E_{\mathrm{T}}^{\gamma}$ selection used to optimize the signal significance spin-2 $q\bar{q}\to X^0 \to Z\gamma$. In addition to the total efficiency, contributions to the signal selection from each of the separate event categories are shown. The efficiencies calculated from MC samples with $W/Z$ hadronic decays are shown as the points on each curve. The line presents interpolated results.

Total efficiencies for the selection of signal events after categorization and application of the tighter photon $E_{\mathrm{T}}^{\gamma}$ selection used to optimize the signal significance spin-2 $gg\to X^0 \to Z\gamma$. In addition to the total efficiency, contributions to the signal selection from each of the separate event categories are shown. The efficiencies calculated from MC samples with $W/Z$ hadronic decays are shown as the points on each curve. The line presents interpolated results.

The $m_{J\gamma}$ distributions of data events selected for the spin-1 $q\bar{q'}\to X^{\pm} \to W^{\pm}\gamma$ search in the D2 category. The background-only fit function shape is shown as the solid curve overlaid with a shaded band corresponding to statistical uncertainties in background parameters. Various signal shapes with cross-sections corresponding to expected limits obtained in this analysis are shown as dashed lines.

The $m_{J\gamma}$ distributions of data events selected for the spin-1 $q\bar{q'}\to X^{\pm} \to W^{\pm}\gamma$ search in the WMASS category. The background-only fit function shape is shown as the solid curve overlaid with a shaded band corresponding to statistical uncertainties in background parameters. Various signal shapes with cross-sections corresponding to expected limits obtained in this analysis are shown as dashed lines.

The $m_{J\gamma}$ distributions of data events selected for the spin-0 $gg \to X^{0} \to Z\gamma$ search in the BTAG category. The background-only fit function shape is shown as the solid curve overlaid with a shaded band corresponding to statistical uncertainties in background parameters. Various signal shapes with cross-sections corresponding to expected limits obtained in this analysis are shown as dashed lines.

The $m_{J\gamma}$ distributions of data events selected for the spin-0 $gg \to X^{0} \to Z\gamma$ search in the D2 category. The background-only fit function shape is shown as the solid curve overlaid with a shaded band corresponding to statistical uncertainties in background parameters. Various signal shapes with cross-sections corresponding to expected limits obtained in this analysis are shown as dashed lines.

The $m_{J\gamma}$ distributions of data events selected for the spin-0 $gg \to X^{0} \to Z\gamma$ search in the ZMASS category. The background-only fit function shape is shown as the solid curve overlaid with a shaded band corresponding to statistical uncertainties in background parameters. Various signal shapes with cross-sections corresponding to expected limits obtained in this analysis are shown as dashed lines.

The $m_{J\gamma}$ distributions of data events selected for the spin-2 $gg \to X^{0} \to Z\gamma$ search in the BTAG category. The background-only fit function shape is shown as the solid curve overlaid with a shaded band corresponding to statistical uncertainties in background parameters. Various signal shapes with cross-sections corresponding to expected limits obtained in this analysis are shown as dashed lines.

The $m_{J\gamma}$ distributions of data events selected for the spin-2 $gg \to X^{0} \to Z\gamma$ search in the D2 category. The background-only fit function shape is shown as the solid curve overlaid with a shaded band corresponding to statistical uncertainties in background parameters. Various signal shapes with cross-sections corresponding to expected limits obtained in this analysis are shown as dashed lines.

The $m_{J\gamma}$ distributions of data events selected for the spin-2 $gg \to X^{0} \to Z\gamma$ search in the ZMASS category. The background-only fit function shape is shown as the solid curve overlaid with a shaded band corresponding to statistical uncertainties in background parameters. Various signal shapes with cross-sections corresponding to expected limits obtained in this analysis are shown as dashed lines.

The $m_{J\gamma}$ distributions of data events selected for the spin-2 $q\bar{q} \to X^{0} \to Z\gamma$ search in the BTAG category. The background-only fit function shape is shown as the solid curve overlaid with a shaded band corresponding to statistical uncertainties in background parameters. Various signal shapes with cross-sections corresponding to expected limits obtained in this analysis are shown as dashed lines.

The $m_{J\gamma}$ distributions of data events selected for the spin-2 $q\bar{q} \to X^{0} \to Z\gamma$ search in the D2 category. The background-only fit function shape is shown as the solid curve overlaid with a shaded band corresponding to statistical uncertainties in background parameters. Various signal shapes with cross-sections corresponding to expected limits obtained in this analysis are shown as dashed lines.

The $m_{J\gamma}$ distributions of data events selected for the spin-2 $q\bar{q} \to X^{0} \to Z\gamma$ search in the ZMASS category. The background-only fit function shape is shown as the solid curve overlaid with a shaded band corresponding to statistical uncertainties in background parameters. Various signal shapes with cross-sections corresponding to expected limits obtained in this analysis are shown as dashed lines.

The 95% CL upper limits on $\sigma(pp\to X)\times B(X\to W/Z\gamma)$ as a function of $m_{X}$ for spin-0 $gg \to X^{0} \to Z\gamma$. The observed limits are shown as a solid black line and the expected ones are shown as a dashed line with the 1$\sigma$ (2$\sigma$) uncertainty band presented as the green (yellow) band.

The 95% CL upper limits on $\sigma(pp\to X)\times B(X\to W/Z\gamma)$ as a function of $m_{X}$ for spin-2 $gg \to X^{0} \to Z\gamma$. The observed limits are shown as a solid black line and the expected ones are shown as a dashed line with the 1$\sigma$ (2$\sigma$) uncertainty band presented as the green (yellow) band.

The 95% CL upper limits on $\sigma(pp\to X)\times B(X\to W/Z\gamma)$ as a function of $m_{X}$ for spin-2 $q\bar{q} \to X^{0} \to Z\gamma$. The observed limits are shown as a solid black line and the expected ones are shown as a dashed line with the 1$\sigma$ (2$\sigma$) uncertainty band presented as the green (yellow) band.

The 95% CL upper limits on $\sigma(pp\to X)\times B(X\to W/Z\gamma)$ as a function of $m_{X}$ for spin-1 $q\bar{q'}\to X^{\pm} \to W^{\pm}\gamma$. The observed limits are shown as a solid black line and the expected ones are shown as a dashed line with the 1$\sigma$ (2$\sigma$) uncertainty band presented as the green (yellow) band.

The jet mass distribution of large-$R$ jets originating from the hadronic decay of $W$ and $Z$ bosons produced from the decay of BSM bosons with mass $m_X = 2000$ GeV. The decays simulated are for the production models $q\bar{q'}\to X^{\pm} \to W^{\pm}\gamma$ with a spin-1 resonance $X^{\pm}$ and $gg\to X^0 \to Z\gamma$ with a spin-0 resonance $X^{0}$.

The photon $E_{T}$ selection criteria as a function of resonance mass for D2 categories. Selections for four different signal channels are shown together.

The photon $E_{T}$ selection criteria as a function of resonance mass for ZMASS or WMASS categories. Selections for four different signal channels are shown together.

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 95% CL exclusion contours in the $(m(W_{R}), m(N_{R}))$ plane in the muon channel for boosted.

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

Acceptance $\times$ efficiency for bSR1e.

Acceptance $\times$ efficiency for bSR2e.

Acceptance $\times$ efficiency for bSR2mu.

Acceptance $\times$ efficiency for rSROS2e.

Acceptance $\times$ efficiency for rSROS2mu.

Acceptance $\times$ efficiency for rSRSS2e.

Acceptance $\times$ efficiency for rSRSS2mu.

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

The $\Delta p_{\mathrm{T}}$ of photon-tagged jets as a function of jet $p_{\mathrm{T}}$.

The double $R_{\mathrm{AA}}$ ratio of $R_{\mathrm{AA}}^{\gamma\text{-jet}}/R_{\mathrm{AA}}^{\text{inclusive jet}}$ as a function of jet $p_{\mathrm{T}}$.

The differential cross-section of photon-tagged jets as a function of jet $p_{\mathrm{T}}$ in pp collisions.

The yields of photon-tagged jets as a function of jet $p_{\mathrm{T}}$ in Pb+Pb collisions for different centrality intervals.

The nuclear modification factor of photon-tagged jets as a function of jet $p_{\mathrm{T}}$ for different centrality intervals.

The $\Delta p_{\mathrm{T}}$ of photon-tagged jets as a function of jet $p_{\mathrm{T}}$.

The $R_{\mathrm{CP}}$ of photon-tagged jets as a function of jet $p_{\mathrm{T}}$.

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