The geometric acceptance multiplied by the efficiency of the $\it{p}_{T}^\text{miss} >$ 200 GeV selection as a function of the proper decay length c$\tau$ for a scalar particle S with a mass of 40 GeV. The acceptance region for DT is defined by requiring the LLP to decay in the region with |z| < 661 cm and 380 cm < r < 736 cm. The acceptance region for CSC is defined by requiring the LLP decay in the region with $|\eta| < 2.4$, r < 695.5 cm, and 661 cm < |z| < 1100 cm or in the region with $|\eta| < 2.4$, r < 270 cm, and 500 cm < |z| < 661 cm. Single CSC cluster requires exactly one LLP to decay in CSC; Single DT cluster requires exactly one LLP to decay in DT; Double cluster requires both LLP to decay in CSC or DT. The denominator in this plot includes all generated events. The nominator includes events that pass the acceptance requirements above and $\it{p}_{T}^\text{miss} >$ 200 GeV.

Distributions of the cluster time for signal, where S decaying to $d\bar{d}$ for a proper decay length c$\tau$ of 1 m and mass of 40 GeV, and for a background-enriched sample in data selected by inverting the $N_\text{hits}$ requirement.

The distributions of $N_\text{hits}$ for single CSC clusters are shown for signal (assuming B(H $\rightarrow$ SS) = 1%, S $\rightarrow d\bar{d}$, and c$\tau$ = 1 m).

The distributions of $N_\text{hits}$ for single CSC clusters are shown for compared to the OOT background ($t_\text{clusters} < 12.5$ ns). The OOT background is representative of the overall background shape, because the background passing all the selections described above is dominated by pileup and underlying events.

The distributions of $\Delta\phi(\it{p}_{T}^\text{miss} \text{,cluster)}$ for single CSC clusters are shown for signal (assuming B(H $\rightarrow$ SS) = 1%, S $\rightarrow d\bar{d}$, and c$\tau$ = 1 m), compared to the shape of background in a selection in which the cluster is not matched to any RPC hit and passes all other selections. The background is dominated by clusters from noise and low-pT particles.

The distributions of $N_\text{hits}$ for single CSC clusters are shown for signal (assuming B(H $\rightarrow$ SS) = 1%, S $\rightarrow d\bar{d}$, and c$\tau$ = 1 m), compared to the shape of background in a selection in which the cluster is not matched to any RPC hit and passes all other selections. The background is dominated by clusters from noise and low-pT particles.

The distributions of $\Delta\phi(\it{p}_{T}^\text{miss} \text{,cluster)}$ for single CSC clusters are shown for signal (assuming B(H $\rightarrow$ SS) = 1%, S $\rightarrow d\bar{d}$, and c$\tau$ = 1 m), compared to the shape of background in a selection in which the cluster is not matched to any RPC hit and passes all other selections. The background is dominated by clusters from noise and low-pT particles.

The signal (assuming B(H $\rightarrow$ SS) = 1%, S $\rightarrow d\bar{d}$, and c$\tau$ = 1 m), background, and data distributions of $N_\text{clusters}$ passing the $N_\text{hits}$ selection in the search region for CSC-CSC category.

The signal (assuming B(H $\rightarrow$ SS) = 1%, S $\rightarrow d\bar{d}$, and c$\tau$ = 1 m), background, and data distributions of $N_\text{clusters}$ passing the $N_\text{hits}$ selection in the search region for DT-DT category.

The signal (assuming B(H $\rightarrow$ SS) = 1%, S $\rightarrow d\bar{d}$, and c$\tau$ = 1 m), background, and data distributions of $N_\text{clusters}$ passing the $N_\text{hits}$ selection in the search region for DT-CSC category.

The signal (assuming B(H $\rightarrow$ SS) = 1%, S $\rightarrow d\bar{d}$, and c$\tau$ = 1 m), background, and data distributions of $N_\text{hits}$ in the search region of the single-CSC cluster category are shown. The $N_\text{hits}$ distribution includes only events in bins A and D. The right-hand bin in the Nhits distribution includes overflow events.

The signal (assuming B(H $\rightarrow$ SS) = 1%, S $\rightarrow d\bar{d}$, and c$\tau$ = 1 m), background, and data distributions of $\Delta\phi(\it{p}_{T}^\text{miss} \text{,cluster)}$ in the search region of the single-CSC cluster category are shown. The $\Delta\phi(\it{p}_{T}^\text{miss} \text{,cluster)}$ distribution includes only events in bins A and B.

The signal (assuming B(H $\rightarrow$ SS) = 1%, S $\rightarrow d\bar{d}$, and c$\tau$ = 1 m), background, and data distributions of $N_\text{hits}$ in the search region of the single-DT cluster category are shown. The $N_\text{hits}$ distribution includes only events in bins A and D. The right-hand bin in the Nhits distribution includes overflow events.

The signal (assuming B(H $\rightarrow$ SS) = 1%, S $\rightarrow d\bar{d}$, and c$\tau$ = 1 m), background, and data distributions of $\Delta\phi(\it{p}_{T}^\text{miss} \text{,cluster)}$ in the search region of the single-DT cluster category are shown. The $\Delta\phi(\it{p}_{T}^\text{miss} \text{,cluster)}$ distribution includes only events in bins A and B.

The 95% CL observed and expected limits on the branching fraction B(H $\rightarrow$ SS) for 3 GeV mass and $ S \rightarrow d\bar{d}$ decay mode.

The 95% CL observed and expected limits on the branching fraction B(H $\rightarrow$ SS) for 7 GeV mass and $ S \rightarrow d\bar{d}$ decay mode.

The 95% CL observed and expected limits on the branching fraction B(H $\rightarrow$ SS) for 15 GeV mass and $ S \rightarrow d\bar{d}$ decay mode.

The 95% CL observed and expected limits on the branching fraction B(H $\rightarrow$ SS) for 40 GeV mass and $ S \rightarrow d\bar{d}$ decay mode.

The 95% CL observed and expected limits on the branching fraction B(H $\rightarrow$ SS) for 55 GeV mass and $ S \rightarrow d\bar{d}$ decay mode.

The 95% CL observed and expected limits on the branching fraction B(H $\rightarrow$ SS) for 0.4 GeV mass and $ S \rightarrow \pi^{0} \pi^{0}$ decay mode.

The 95% CL observed and expected limits on the branching fraction B(H $\rightarrow$ SS) for 1 GeV mass and $ S \rightarrow \pi^{0} \pi^{0}$ decay mode.

The 95% CL observed and expected limits on the branching fraction B(H $\rightarrow$ SS) for 7 GeV mass and $ S \rightarrow \tau^{+} \tau^{-}$ decay mode.

The 95% CL observed and expected limits on the branching fraction B(H $\rightarrow$ SS) for 15 GeV mass and $ S \rightarrow \tau^{+} \tau^{-}$ decay mode.

The 95% CL observed and expected limits on the branching fraction B(H $\rightarrow$ SS) for 40 GeV mass and $ S \rightarrow \tau^{+} \tau^{-}$ decay mode.

The 95% CL observed and expected limits on the branching fraction B(H $\rightarrow$ SS) for 55 GeV mass and $ S \rightarrow \tau^{+} \tau^{-}$ decay mode.

The 95% CL observed and expected limits on the branching fraction B(H $\rightarrow$ SS) for 15 GeV mass and $ S \rightarrow b\bar{b}$ decay mode.

The 95% CL observed and expected limits on the branching fraction B(H $\rightarrow$ SS) for 40 GeV mass and $ S \rightarrow b\bar{b}$ decay mode.

The 95% CL observed and expected limits on the branching fraction B(H $\rightarrow$ SS) for 55 GeV mass and $ S \rightarrow b\bar{b}$ decay mode.

The 95% CL observed and expected limits on the branching fraction B(H $\rightarrow$ SS) for 0.4 GeV mass and $ S \rightarrow \pi^{+} \pi^{-}$ decay mode.

The 95% CL observed and expected limits on the branching fraction B(H $\rightarrow$ SS) for 1 GeV mass and $ S \rightarrow \pi^{+} \pi^{-}$ decay mode.

The 95% CL observed and expected limits on the branching fraction B(H $\rightarrow$ SS) for 1.5 GeV mass and $ S \rightarrow K^{+}K^{-}$ decay mode.

The 95% CL observed and expected limits on the branching fraction B(H $\rightarrow$ SS) for 1.5 GeV mass and $ S \rightarrow K^{0}K^{0}$ decay mode.

The 95% CL observed and expected limits on the branching fraction B(H $\rightarrow$ SS) for 0.4 GeV mass and $ S \rightarrow \gamma\gamma$ decay mode.

The 95% CL observed and expected limits on the branching fraction B(H $\rightarrow$ SS) for 0.4 GeV mass and $ S \rightarrow \e^{+} \e^{-}$ decay mode.

The 95% CL observed and expected limits on the branching fraction B(H $\rightarrow$ SS) as functions of mass and c$\tau$ assuming S inherits all couplings from the Higgs boson evaluated at the LLP mass.

The 95% CL observed and expected limits on the branching fraction B(H $\rightarrow$ SS) as functions of mass and c$\tau$ assuming S inherits all couplings from the Higgs boson evaluated at the LLP mass.

The 95% CL observed and expected limits on the branching fraction B(H $\rightarrow$ SS) as functions of mass and c$\tau$ assuming S inherits all couplings from the Higgs boson evaluated at the LLP mass.

The 95% CL observed and expected limits on the branching fraction B(H $\rightarrow$ SS) as functions of mass and c$\tau$ assuming S inherits all couplings from the Higgs boson evaluated at the LLP mass.

The 95% CL observed and expected limits on the branching fraction B(H $\rightarrow$ SS) as functions of mass and c$\tau$ assuming S inherits all couplings from the Higgs boson evaluated at the LLP mass.

The 95% CL observed and expected limits on the branching fraction B(H $\rightarrow$ SS) as functions of mass and c$\tau$ assuming S inherits all couplings from the Higgs boson evaluated at the LLP mass.

The 95% CL observed and expected limits on the branching fraction B(H $\rightarrow$ SS) as functions of mass and c$\tau$ assuming S inherits all couplings from the Higgs boson evaluated at the LLP mass.

The 95% CL observed and expected limits on the branching fraction B(H $\rightarrow$ SS) as functions of mass and c$\tau$ assuming S inherits all couplings from the Higgs boson evaluated at the LLP mass.

The 95% CL observed and expected limits on the branching fraction B(H $\rightarrow$ SS) as functions of mass and c$\tau$ assuming S inherits all couplings from the Higgs boson evaluated at the LLP mass.

The 95% CL observed and expected limits on the branching fraction B(H $\rightarrow$ SS) as functions of mass and c$\tau$ assuming S inherits all couplings from the Higgs boson evaluated at the LLP mass.

The 95% CL observed and expected limits on the branching fraction B(H $\rightarrow \Psi\Psi$) for vector portal, assuming $(\xi_{\omega}$, $\xi_{\Lambda}) = (1, 1)$ and 2 GeV LLP mass.

The 95% CL observed and expected limits on the branching fraction B(H $\rightarrow \Psi\Psi$) for vector portal, assuming $(\xi_{\omega}$, $\xi_{\Lambda}) = (1, 1)$ and 5 GeV LLP mass.

The 95% CL observed and expected limits on the branching fraction B(H $\rightarrow \Psi\Psi$) for vector portal, assuming $(\xi_{\omega}$, $\xi_{\Lambda}) = (1, 1)$ and 10 GeV LLP mass.

The 95% CL observed and expected limits on the branching fraction B(H $\rightarrow \Psi\Psi$) for vector portal, assuming $(\xi_{\omega}$, $\xi_{\Lambda}) = (1, 1)$ and 15 GeV LLP mass.

The 95% CL observed and expected limits on the branching fraction B(H $\rightarrow \Psi\Psi$) for vector portal, assuming $(\xi_{\omega}$, $\xi_{\Lambda}) = (1, 1)$ and 20 GeV LLP mass.

The 95% CL observed and expected limits on the branching fraction B(H $\rightarrow \Psi\Psi$) for gluon portal, assuming $(\xi_{\omega}$, $\xi_{\Lambda}) = (2.5, 1)$ and 5 GeV LLP mass.

The 95% CL observed and expected limits on the branching fraction B(H $\rightarrow \Psi\Psi$) for gluon portal, assuming $(\xi_{\omega}$, $\xi_{\Lambda}) = (2.5, 1)$ and 10 GeV LLP mass.

The 95% CL observed and expected limits on the branching fraction B(H $\rightarrow \Psi\Psi$) for gluon portal, assuming $(\xi_{\omega}$, $\xi_{\Lambda}) = (2.5, 1)$ and 15 GeV LLP mass.

The 95% CL observed and expected limits on the branching fraction B(H $\rightarrow \Psi\Psi$) for gluon portal, assuming $(\xi_{\omega}$, $\xi_{\Lambda}) = (2.5, 1)$ and 20 GeV LLP mass.

The 95% CL observed and expected limits on the branching fraction B(H $\rightarrow \Psi\Psi$) for gluon portal, assuming $(\xi_{\omega}$, $\xi_{\Lambda}) = (2.5, 2.5)$ and 5 GeV LLP mass.

The 95% CL observed and expected limits on the branching fraction B(H $\rightarrow \Psi\Psi$) for gluon portal, assuming $(\xi_{\omega}$, $\xi_{\Lambda}) = (2.5, 2.5)$ and 10 GeV LLP mass.

The 95% CL observed and expected limits on the branching fraction B(H $\rightarrow \Psi\Psi$) for gluon portal, assuming $(\xi_{\omega}$, $\xi_{\Lambda}) = (2.5, 2.5)$ and 15 GeV LLP mass.

The 95% CL observed and expected limits on the branching fraction B(H $\rightarrow \Psi\Psi$) for gluon portal, assuming $(\xi_{\omega}$, $\xi_{\Lambda}) = (2.5, 2.5)$ and 20 GeV LLP mass.

The 95% CL observed and expected limits on the branching fraction B(H $\rightarrow \Psi\Psi$) for gluon portal, assuming $(\xi_{\omega}$, $\xi_{\Lambda}) = (1, 1)$ and 5 GeV LLP mass.

The 95% CL observed and expected limits on the branching fraction B(H $\rightarrow \Psi\Psi$) for gluon portal, assuming $(\xi_{\omega}$, $\xi_{\Lambda}) = (1, 1)$ and 10 GeV LLP mass.

The 95% CL observed and expected limits on the branching fraction B(H $\rightarrow \Psi\Psi$) for gluon portal, assuming $(\xi_{\omega}$, $\xi_{\Lambda}) = (1, 1)$ and 15 GeV LLP mass.

The 95% CL observed and expected limits on the branching fraction B(H $\rightarrow \Psi\Psi$) for gluon portal, assuming $(\xi_{\omega}$, $\xi_{\Lambda}) = (1, 1)$ and 20 GeV LLP mass.

The 95% CL observed and expected limits on the branching fraction B(H $\rightarrow \Psi\Psi$) for photon portal, assuming $(\xi_{\omega}$, $\xi_{\Lambda}) = (2.5, 1)$ and 2 GeV LLP mass.

The 95% CL observed and expected limits on the branching fraction B(H $\rightarrow \Psi\Psi$) for photon portal, assuming $(\xi_{\omega}$, $\xi_{\Lambda}) = (2.5, 1)$ and 5 GeV LLP mass.

The 95% CL observed and expected limits on the branching fraction B(H $\rightarrow \Psi\Psi$) for photon portal, assuming $(\xi_{\omega}$, $\xi_{\Lambda}) = (2.5, 1)$ and 10 GeV LLP mass.

The 95% CL observed and expected limits on the branching fraction B(H $\rightarrow \Psi\Psi$) for photon portal, assuming $(\xi_{\omega}$, $\xi_{\Lambda}) = (2.5, 1)$ and 15 GeV LLP mass.

The 95% CL observed and expected limits on the branching fraction B(H $\rightarrow \Psi\Psi$) for photon portal, assuming $(\xi_{\omega}$, $\xi_{\Lambda}) = (2.5, 1)$ and 20 GeV LLP mass.

The 95% CL observed and expected limits on the branching fraction B(H $\rightarrow \Psi\Psi$) for photon portal, assuming $(\xi_{\omega}$, $\xi_{\Lambda}) = (2.5, 2.5)$ and 2 GeV LLP mass.

The 95% CL observed and expected limits on the branching fraction B(H $\rightarrow \Psi\Psi$) for photon portal, assuming $(\xi_{\omega}$, $\xi_{\Lambda}) = (2.5, 2.5)$ and 5 GeV LLP mass.

The 95% CL observed and expected limits on the branching fraction B(H $\rightarrow \Psi\Psi$) for photon portal, assuming $(\xi_{\omega}$, $\xi_{\Lambda}) = (2.5, 2.5)$ and 10 GeV LLP mass.

The 95% CL observed and expected limits on the branching fraction B(H $\rightarrow \Psi\Psi$) for photon portal, assuming $(\xi_{\omega}$, $\xi_{\Lambda}) = (2.5, 2.5)$ and 15 GeV LLP mass.

The 95% CL observed and expected limits on the branching fraction B(H $\rightarrow \Psi\Psi$) for photon portal, assuming $(\xi_{\omega}$, $\xi_{\Lambda}) = (2.5, 2.5)$ and 20 GeV LLP mass.

The 95% CL observed and expected limits on the branching fraction B(H $\rightarrow \Psi\Psi$) for photon portal, assuming $(\xi_{\omega}$, $\xi_{\Lambda}) = (1, 1)$ and 2 GeV LLP mass.

The 95% CL observed and expected limits on the branching fraction B(H $\rightarrow \Psi\Psi$) for photon portal, assuming $(\xi_{\omega}$, $\xi_{\Lambda}) = (1, 1)$ and 5 GeV LLP mass.

The 95% CL observed and expected limits on the branching fraction B(H $\rightarrow \Psi\Psi$) for photon portal, assuming $(\xi_{\omega}$, $\xi_{\Lambda}) = (1, 1)$ and 10 GeV LLP mass.

The 95% CL observed and expected limits on the branching fraction B(H $\rightarrow \Psi\Psi$) for photon portal, assuming $(\xi_{\omega}$, $\xi_{\Lambda}) = (1, 1)$ and 15 GeV LLP mass.

The 95% CL observed and expected limits on the branching fraction B(H $\rightarrow \Psi\Psi$) for photon portal, assuming $(\xi_{\omega}$, $\xi_{\Lambda}) = (1, 1)$ and 20 GeV LLP mass.

The 95% CL observed and expected limits on the branching fraction B(H $\rightarrow \Psi\Psi$) for higgs portal, assuming $(\xi_{\omega}$, $\xi_{\Lambda}) = (2.5, 1)$ and 4 GeV LLP mass.

The 95% CL observed and expected limits on the branching fraction B(H $\rightarrow \Psi\Psi$) for higgs portal, assuming $(\xi_{\omega}$, $\xi_{\Lambda}) = (2.5, 1)$ and 5 GeV LLP mass.

The 95% CL observed and expected limits on the branching fraction B(H $\rightarrow \Psi\Psi$) for higgs portal, assuming $(\xi_{\omega}$, $\xi_{\Lambda}) = (2.5, 1)$ and 10 GeV LLP mass.

The 95% CL observed and expected limits on the branching fraction B(H $\rightarrow \Psi\Psi$) for higgs portal, assuming $(\xi_{\omega}$, $\xi_{\Lambda}) = (2.5, 1)$ and 15 GeV LLP mass.

The 95% CL observed and expected limits on the branching fraction B(H $\rightarrow \Psi\Psi$) for higgs portal, assuming $(\xi_{\omega}$, $\xi_{\Lambda}) = (2.5, 1)$ and 20 GeV LLP mass.

The 95% CL observed and expected limits on the branching fraction B(H $\rightarrow \Psi\Psi$) for higgs portal, assuming $(\xi_{\omega}$, $\xi_{\Lambda}) = (2.5, 2.5)$ and 4 GeV LLP mass.

The 95% CL observed and expected limits on the branching fraction B(H $\rightarrow \Psi\Psi$) for higgs portal, assuming $(\xi_{\omega}$, $\xi_{\Lambda}) = (2.5, 2.5)$ and 5 GeV LLP mass.

The 95% CL observed and expected limits on the branching fraction B(H $\rightarrow \Psi\Psi$) for higgs portal, assuming $(\xi_{\omega}$, $\xi_{\Lambda}) = (2.5, 2.5)$ and 10 GeV LLP mass.

The 95% CL observed and expected limits on the branching fraction B(H $\rightarrow \Psi\Psi$) for higgs portal, assuming $(\xi_{\omega}$, $\xi_{\Lambda}) = (2.5, 2.5)$ and 15 GeV LLP mass.

The 95% CL observed and expected limits on the branching fraction B(H $\rightarrow \Psi\Psi$) for higgs portal, assuming $(\xi_{\omega}$, $\xi_{\Lambda}) = (2.5, 2.5)$ and 20 GeV LLP mass.

The 95% CL observed and expected limits on the branching fraction B(H $\rightarrow \Psi\Psi$) for higgs portal, assuming $(\xi_{\omega}$, $\xi_{\Lambda}) = (1, 1)$ and 4 GeV LLP mass.

The 95% CL observed and expected limits on the branching fraction B(H $\rightarrow \Psi\Psi$) for higgs portal, assuming $(\xi_{\omega}$, $\xi_{\Lambda}) = (1, 1)$ and 5 GeV LLP mass.

The 95% CL observed and expected limits on the branching fraction B(H $\rightarrow \Psi\Psi$) for higgs portal, assuming $(\xi_{\omega}$, $\xi_{\Lambda}) = (1, 1)$ and 10 GeV LLP mass.

The 95% CL observed and expected limits on the branching fraction B(H $\rightarrow \Psi\Psi$) for higgs portal, assuming $(\xi_{\omega}$, $\xi_{\Lambda}) = (1, 1)$ and 15 GeV LLP mass.

The 95% CL observed and expected limits on the branching fraction B(H $\rightarrow \Psi\Psi$) for higgs portal, assuming $(\xi_{\omega}$, $\xi_{\Lambda}) = (1, 1)$ and 20 GeV LLP mass.

The 95% CL observed and expected limits on the branching fraction B(H $\rightarrow \Psi\Psi$) for darkphoton portal, assuming $(\xi_{\omega}$, $\xi_{\Lambda}) = (2.5, 1)$ and 5 GeV LLP mass.

The 95% CL observed and expected limits on the branching fraction B(H $\rightarrow \Psi\Psi$) for darkphoton portal, assuming $(\xi_{\omega}$, $\xi_{\Lambda}) = (2.5, 1)$ and 10 GeV LLP mass.

The 95% CL observed and expected limits on the branching fraction B(H $\rightarrow \Psi\Psi$) for darkphoton portal, assuming $(\xi_{\omega}$, $\xi_{\Lambda}) = (2.5, 1)$ and 15 GeV LLP mass.

The 95% CL observed and expected limits on the branching fraction B(H $\rightarrow \Psi\Psi$) for darkphoton portal, assuming $(\xi_{\omega}$, $\xi_{\Lambda}) = (2.5, 1)$ and 20 GeV LLP mass.

The 95% CL observed and expected limits on the branching fraction B(H $\rightarrow \Psi\Psi$) for darkphoton portal, assuming $(\xi_{\omega}$, $\xi_{\Lambda}) = (2.5, 2.5)$ and 5 GeV LLP mass.

The 95% CL observed and expected limits on the branching fraction B(H $\rightarrow \Psi\Psi$) for darkphoton portal, assuming $(\xi_{\omega}$, $\xi_{\Lambda}) = (2.5, 2.5)$ and 10 GeV LLP mass.

The 95% CL observed and expected limits on the branching fraction B(H $\rightarrow \Psi\Psi$) for darkphoton portal, assuming $(\xi_{\omega}$, $\xi_{\Lambda}) = (2.5, 2.5)$ and 15 GeV LLP mass.

The 95% CL observed and expected limits on the branching fraction B(H $\rightarrow \Psi\Psi$) for darkphoton portal, assuming $(\xi_{\omega}$, $\xi_{\Lambda}) = (2.5, 2.5)$ and 20 GeV LLP mass.

The 95% CL observed and expected limits on the branching fraction B(H $\rightarrow \Psi\Psi$) for darkphoton portal, assuming $(\xi_{\omega}$, $\xi_{\Lambda}) = (1, 1)$ and 5 GeV LLP mass.

The 95% CL observed and expected limits on the branching fraction B(H $\rightarrow \Psi\Psi$) for darkphoton portal, assuming $(\xi_{\omega}$, $\xi_{\Lambda}) = (1, 1)$ and 10 GeV LLP mass.

The 95% CL observed and expected limits on the branching fraction B(H $\rightarrow \Psi\Psi$) for darkphoton portal, assuming $(\xi_{\omega}$, $\xi_{\Lambda}) = (1, 1)$ and 15 GeV LLP mass.

The 95% CL observed and expected limits on the branching fraction B(H $\rightarrow \Psi\Psi$) for darkphoton portal, assuming $(\xi_{\omega}$, $\xi_{\Lambda}) = (1, 1)$ and 20 GeV LLP mass.

Expected (dashed black line) and observed (solid red line) 95% CL exclusion limits on the higgsino simplified model being considered. These are shown with $\pm 1\sigma_{\mathrm{exp}}$ (yellow band) from experimental systematic and statistical uncertainties, and with $\pm 1\sigma_{\mathrm{theory}}^{\mathrm{SUSY}}$ (red dotted lines) from signal cross-section uncertainties, respectively. The limits set by the latest ATLAS searches using the soft lepton and disappearing track signatures are illustrated by the blue and green regions, respectively, while the limit imposed by the LEP experiments is shown in gray. The dot-dashed gray line indicates the predicted mass-splitting for the pure higgsino scenario.

Expected (dashed black line) and observed (solid red line) 95% CL exclusion limits on the higgsino simplified model being considered. These are shown with $\pm 1\sigma_{\mathrm{exp}}$ (yellow band) from experimental systematic and statistical uncertainties, and with $\pm 1\sigma_{\mathrm{theory}}^{\mathrm{SUSY}}$ (red dotted lines) from signal cross-section uncertainties, respectively. The limits set by the latest ATLAS searches using the soft lepton and disappearing track signatures are illustrated by the blue and green regions, respectively, while the limit imposed by the LEP experiments is shown in gray. The dot-dashed gray line indicates the predicted mass-splitting for the pure higgsino scenario.

Expected (dashed black line) and observed (solid red line) 95% CL exclusion limits on the higgsino simplified model being considered. These are shown with $\pm 1\sigma_{\mathrm{exp}}$ (yellow band) from experimental systematic and statistical uncertainties, and with $\pm 1\sigma_{\mathrm{theory}}^{\mathrm{SUSY}}$ (red dotted lines) from signal cross-section uncertainties, respectively. The limits set by the latest ATLAS searches using the soft lepton and disappearing track signatures are illustrated by the blue and green regions, respectively, while the limit imposed by the LEP experiments is shown in gray. The dot-dashed gray line indicates the predicted mass-splitting for the pure higgsino scenario.

Expected (dashed black line) and observed (solid red line) 95% CL exclusion limits on the higgsino simplified model being considered. These are shown with $\pm 1\sigma_{\mathrm{exp}}$ (yellow band) from experimental systematic and statistical uncertainties, and with $\pm 1\sigma_{\mathrm{theory}}^{\mathrm{SUSY}}$ (red dotted lines) from signal cross-section uncertainties, respectively. The limits set by the latest ATLAS searches using the soft lepton and disappearing track signatures are illustrated by the blue and green regions, respectively, while the limit imposed by the LEP experiments is shown in gray. The dot-dashed gray line indicates the predicted mass-splitting for the pure higgsino scenario.

Expected and observed CLs values per signal point represented by the grey numbers. The expected (dashed) and observed (solid) 95% CL exclusion limits are overlaid along with $\pm 1\sigma_{\mathrm{exp}}$ (yellow band) from experimental systematic and statistical uncertainties, and with $\pm 1\sigma_{\mathrm{theory}}^{\mathrm{SUSY}}$ (red dotted lines) from signal cross-section uncertainties, respectively.

Expected and observed CLs values per signal point represented by the grey numbers. The expected (dashed) and observed (solid) 95% CL exclusion limits are overlaid along with $\pm 1\sigma_{\mathrm{exp}}$ (yellow band) from experimental systematic and statistical uncertainties, and with $\pm 1\sigma_{\mathrm{theory}}^{\mathrm{SUSY}}$ (red dotted lines) from signal cross-section uncertainties, respectively.

Expected and observed cross-section upper-limit per signal point represented by the grey numbers. The expected (dashed) and observed (solid) 95% CL exclusion limits are overlaid along with $\pm 1\sigma_{\mathrm{exp}}$ (yellow band) from experimental systematic and statistical uncertainties, and with $\pm 1\sigma_{\mathrm{theory}}^{\mathrm{SUSY}}$ (red dotted lines) from signal cross-section uncertainties, respectively.

Expected and observed cross-section upper-limit per signal point represented by the grey numbers. The expected (dashed) and observed (solid) 95% CL exclusion limits are overlaid along with $\pm 1\sigma_{\mathrm{exp}}$ (yellow band) from experimental systematic and statistical uncertainties, and with $\pm 1\sigma_{\mathrm{theory}}^{\mathrm{SUSY}}$ (red dotted lines) from signal cross-section uncertainties, respectively.

Truth-level signal acceptances for each production process ($\tilde{\chi}_1^\pm \tilde{\chi}_1^0$, $\tilde{\chi}_1^\pm \tilde{\chi}_2^0$, $\tilde{\chi}_1^+ \tilde{\chi}_1^-$, and $\tilde{\chi}_2^0 \tilde{\chi}_1^0$) in a SR with the $S(d_0)$ requirement removed. The acceptance is defined as the fraction of accepted events divided by the total number of events in the generator-level signal Monte Carlo simulation, where the signal candidate track is identified as the charged particle with the largest distance between the interaction vertex and the secondary vertex of the higgsino decays.

Truth-level signal acceptances for each production process ($\tilde{\chi}_1^\pm \tilde{\chi}_1^0$, $\tilde{\chi}_1^\pm \tilde{\chi}_2^0$, $\tilde{\chi}_1^+ \tilde{\chi}_1^-$, and $\tilde{\chi}_2^0 \tilde{\chi}_1^0$) in a SR with the $S(d_0)$ requirement removed. The acceptance is defined as the fraction of accepted events divided by the total number of events in the generator-level signal Monte Carlo simulation, where the signal candidate track is identified as the charged particle with the largest distance between the interaction vertex and the secondary vertex of the higgsino decays.

Truth-level signal acceptances for each production process ($\tilde{\chi}_1^\pm \tilde{\chi}_1^0$, $\tilde{\chi}_1^\pm \tilde{\chi}_2^0$, $\tilde{\chi}_1^+ \tilde{\chi}_1^-$, and $\tilde{\chi}_2^0 \tilde{\chi}_1^0$) in a SR with the $S(d_0)$ requirement removed. The acceptance is defined as the fraction of accepted events divided by the total number of events in the generator-level signal Monte Carlo simulation, where the signal candidate track is identified as the charged particle with the largest distance between the interaction vertex and the secondary vertex of the higgsino decays.

Truth-level signal acceptances for each production process ($\tilde{\chi}_1^\pm \tilde{\chi}_1^0$, $\tilde{\chi}_1^\pm \tilde{\chi}_2^0$, $\tilde{\chi}_1^+ \tilde{\chi}_1^-$, and $\tilde{\chi}_2^0 \tilde{\chi}_1^0$) in a SR with the $S(d_0)$ requirement removed. The acceptance is defined as the fraction of accepted events divided by the total number of events in the generator-level signal Monte Carlo simulation, where the signal candidate track is identified as the charged particle with the largest distance between the interaction vertex and the secondary vertex of the higgsino decays.

Truth-level signal acceptances for each production process ($\tilde{\chi}_1^\pm \tilde{\chi}_1^0$, $\tilde{\chi}_1^\pm \tilde{\chi}_2^0$, $\tilde{\chi}_1^+ \tilde{\chi}_1^-$, and $\tilde{\chi}_2^0 \tilde{\chi}_1^0$) in a SR with the $S(d_0)$ requirement removed. The acceptance is defined as the fraction of accepted events divided by the total number of events in the generator-level signal Monte Carlo simulation, where the signal candidate track is identified as the charged particle with the largest distance between the interaction vertex and the secondary vertex of the higgsino decays.

Truth-level signal acceptances for each production process ($\tilde{\chi}_1^\pm \tilde{\chi}_1^0$, $\tilde{\chi}_1^\pm \tilde{\chi}_2^0$, $\tilde{\chi}_1^+ \tilde{\chi}_1^-$, and $\tilde{\chi}_2^0 \tilde{\chi}_1^0$) in a SR with the $S(d_0)$ requirement removed. The acceptance is defined as the fraction of accepted events divided by the total number of events in the generator-level signal Monte Carlo simulation, where the signal candidate track is identified as the charged particle with the largest distance between the interaction vertex and the secondary vertex of the higgsino decays.

Signal efficiencies in SR-Low for each production process ($\tilde{\chi}_1^\pm \tilde{\chi}_1^0$, $\tilde{\chi}_1^\pm \tilde{\chi}_2^0$, $\tilde{\chi}_1^+ \tilde{\chi}_1^-$, and $\tilde{\chi}_2^0 \tilde{\chi}_1^0$), defined by the number of events of reconstructed-level signal simulation divided by the number of events obtained at generator level, where the $S(d_0)$ selecton efficiency has the largest impact. The higgsino decay products from $\Delta \mathrm{m}(\tilde{\chi}_1^\pm,\tilde{\chi}_1^0) < 0.4$ GeV signal have $p_{\mathrm{T}}$ too low to be reconstructed as the signal candidate tracks, and therefore the identified signal candidate tracks are typically from pile-up collisions or underlying events similar to the QCD track background, causing a low $S(d_0)$ selection efficiency in these plots.

Signal efficiencies in SR-Low for each production process ($\tilde{\chi}_1^\pm \tilde{\chi}_1^0$, $\tilde{\chi}_1^\pm \tilde{\chi}_2^0$, $\tilde{\chi}_1^+ \tilde{\chi}_1^-$, and $\tilde{\chi}_2^0 \tilde{\chi}_1^0$), defined by the number of events of reconstructed-level signal simulation divided by the number of events obtained at generator level, where the $S(d_0)$ selecton efficiency has the largest impact. The higgsino decay products from $\Delta \mathrm{m}(\tilde{\chi}_1^\pm,\tilde{\chi}_1^0) < 0.4$ GeV signal have $p_{\mathrm{T}}$ too low to be reconstructed as the signal candidate tracks, and therefore the identified signal candidate tracks are typically from pile-up collisions or underlying events similar to the QCD track background, causing a low $S(d_0)$ selection efficiency in these plots.

Signal efficiencies in SR-Low for each production process ($\tilde{\chi}_1^\pm \tilde{\chi}_1^0$, $\tilde{\chi}_1^\pm \tilde{\chi}_2^0$, $\tilde{\chi}_1^+ \tilde{\chi}_1^-$, and $\tilde{\chi}_2^0 \tilde{\chi}_1^0$), defined by the number of events of reconstructed-level signal simulation divided by the number of events obtained at generator level, where the $S(d_0)$ selecton efficiency has the largest impact. The higgsino decay products from $\Delta \mathrm{m}(\tilde{\chi}_1^\pm,\tilde{\chi}_1^0) < 0.4$ GeV signal have $p_{\mathrm{T}}$ too low to be reconstructed as the signal candidate tracks, and therefore the identified signal candidate tracks are typically from pile-up collisions or underlying events similar to the QCD track background, causing a low $S(d_0)$ selection efficiency in these plots.

Signal efficiencies in SR-Low for each production process ($\tilde{\chi}_1^\pm \tilde{\chi}_1^0$, $\tilde{\chi}_1^\pm \tilde{\chi}_2^0$, $\tilde{\chi}_1^+ \tilde{\chi}_1^-$, and $\tilde{\chi}_2^0 \tilde{\chi}_1^0$), defined by the number of events of reconstructed-level signal simulation divided by the number of events obtained at generator level, where the $S(d_0)$ selecton efficiency has the largest impact. The higgsino decay products from $\Delta \mathrm{m}(\tilde{\chi}_1^\pm,\tilde{\chi}_1^0) < 0.4$ GeV signal have $p_{\mathrm{T}}$ too low to be reconstructed as the signal candidate tracks, and therefore the identified signal candidate tracks are typically from pile-up collisions or underlying events similar to the QCD track background, causing a low $S(d_0)$ selection efficiency in these plots.

Signal efficiencies in SR-Low for each production process ($\tilde{\chi}_1^\pm \tilde{\chi}_1^0$, $\tilde{\chi}_1^\pm \tilde{\chi}_2^0$, $\tilde{\chi}_1^+ \tilde{\chi}_1^-$, and $\tilde{\chi}_2^0 \tilde{\chi}_1^0$), defined by the number of events of reconstructed-level signal simulation divided by the number of events obtained at generator level, where the $S(d_0)$ selecton efficiency has the largest impact. The higgsino decay products from $\Delta \mathrm{m}(\tilde{\chi}_1^\pm,\tilde{\chi}_1^0) < 0.4$ GeV signal have $p_{\mathrm{T}}$ too low to be reconstructed as the signal candidate tracks, and therefore the identified signal candidate tracks are typically from pile-up collisions or underlying events similar to the QCD track background, causing a low $S(d_0)$ selection efficiency in these plots.

Signal efficiencies in SR-Low for each production process ($\tilde{\chi}_1^\pm \tilde{\chi}_1^0$, $\tilde{\chi}_1^\pm \tilde{\chi}_2^0$, $\tilde{\chi}_1^+ \tilde{\chi}_1^-$, and $\tilde{\chi}_2^0 \tilde{\chi}_1^0$), defined by the number of events of reconstructed-level signal simulation divided by the number of events obtained at generator level, where the $S(d_0)$ selecton efficiency has the largest impact. The higgsino decay products from $\Delta \mathrm{m}(\tilde{\chi}_1^\pm,\tilde{\chi}_1^0) < 0.4$ GeV signal have $p_{\mathrm{T}}$ too low to be reconstructed as the signal candidate tracks, and therefore the identified signal candidate tracks are typically from pile-up collisions or underlying events similar to the QCD track background, causing a low $S(d_0)$ selection efficiency in these plots.

Signal efficiencies in SR-High for each production process ($\tilde{\chi}_1^\pm \tilde{\chi}_1^0$, $\tilde{\chi}_1^\pm \tilde{\chi}_2^0$, $\tilde{\chi}_1^+ \tilde{\chi}_1^-$, and $\tilde{\chi}_2^0 \tilde{\chi}_1^0$), defined by the number of events of reconstructed-level signal simulation divided by the number of events obtained at generator level, where the $S(d_0)$ selecton efficiency has the largest impact. The higgsino decay products from $\Delta \mathrm{m}(\tilde{\chi}_1^\pm,\tilde{\chi}_1^0) < 0.4$ GeV signal have $p_{\mathrm{T}}$ too low to be reconstructed as the signal candidate tracks, and therefore the identified signal candidate tracks are typically from pile-up collisions or underlying events similar to the QCD track background, causing a low $S(d_0)$ selection efficiency in these plots.

Signal efficiencies in SR-High for each production process ($\tilde{\chi}_1^\pm \tilde{\chi}_1^0$, $\tilde{\chi}_1^\pm \tilde{\chi}_2^0$, $\tilde{\chi}_1^+ \tilde{\chi}_1^-$, and $\tilde{\chi}_2^0 \tilde{\chi}_1^0$), defined by the number of events of reconstructed-level signal simulation divided by the number of events obtained at generator level, where the $S(d_0)$ selecton efficiency has the largest impact. The higgsino decay products from $\Delta \mathrm{m}(\tilde{\chi}_1^\pm,\tilde{\chi}_1^0) < 0.4$ GeV signal have $p_{\mathrm{T}}$ too low to be reconstructed as the signal candidate tracks, and therefore the identified signal candidate tracks are typically from pile-up collisions or underlying events similar to the QCD track background, causing a low $S(d_0)$ selection efficiency in these plots.

Signal efficiencies in SR-High for each production process ($\tilde{\chi}_1^\pm \tilde{\chi}_1^0$, $\tilde{\chi}_1^\pm \tilde{\chi}_2^0$, $\tilde{\chi}_1^+ \tilde{\chi}_1^-$, and $\tilde{\chi}_2^0 \tilde{\chi}_1^0$), defined by the number of events of reconstructed-level signal simulation divided by the number of events obtained at generator level, where the $S(d_0)$ selecton efficiency has the largest impact. The higgsino decay products from $\Delta \mathrm{m}(\tilde{\chi}_1^\pm,\tilde{\chi}_1^0) < 0.4$ GeV signal have $p_{\mathrm{T}}$ too low to be reconstructed as the signal candidate tracks, and therefore the identified signal candidate tracks are typically from pile-up collisions or underlying events similar to the QCD track background, causing a low $S(d_0)$ selection efficiency in these plots.

Signal efficiencies in SR-High for each production process ($\tilde{\chi}_1^\pm \tilde{\chi}_1^0$, $\tilde{\chi}_1^\pm \tilde{\chi}_2^0$, $\tilde{\chi}_1^+ \tilde{\chi}_1^-$, and $\tilde{\chi}_2^0 \tilde{\chi}_1^0$), defined by the number of events of reconstructed-level signal simulation divided by the number of events obtained at generator level, where the $S(d_0)$ selecton efficiency has the largest impact. The higgsino decay products from $\Delta \mathrm{m}(\tilde{\chi}_1^\pm,\tilde{\chi}_1^0) < 0.4$ GeV signal have $p_{\mathrm{T}}$ too low to be reconstructed as the signal candidate tracks, and therefore the identified signal candidate tracks are typically from pile-up collisions or underlying events similar to the QCD track background, causing a low $S(d_0)$ selection efficiency in these plots.

Signal efficiencies in SR-High for each production process ($\tilde{\chi}_1^\pm \tilde{\chi}_1^0$, $\tilde{\chi}_1^\pm \tilde{\chi}_2^0$, $\tilde{\chi}_1^+ \tilde{\chi}_1^-$, and $\tilde{\chi}_2^0 \tilde{\chi}_1^0$), defined by the number of events of reconstructed-level signal simulation divided by the number of events obtained at generator level, where the $S(d_0)$ selecton efficiency has the largest impact. The higgsino decay products from $\Delta \mathrm{m}(\tilde{\chi}_1^\pm,\tilde{\chi}_1^0) < 0.4$ GeV signal have $p_{\mathrm{T}}$ too low to be reconstructed as the signal candidate tracks, and therefore the identified signal candidate tracks are typically from pile-up collisions or underlying events similar to the QCD track background, causing a low $S(d_0)$ selection efficiency in these plots.

Signal efficiencies in SR-High for each production process ($\tilde{\chi}_1^\pm \tilde{\chi}_1^0$, $\tilde{\chi}_1^\pm \tilde{\chi}_2^0$, $\tilde{\chi}_1^+ \tilde{\chi}_1^-$, and $\tilde{\chi}_2^0 \tilde{\chi}_1^0$), defined by the number of events of reconstructed-level signal simulation divided by the number of events obtained at generator level, where the $S(d_0)$ selecton efficiency has the largest impact. The higgsino decay products from $\Delta \mathrm{m}(\tilde{\chi}_1^\pm,\tilde{\chi}_1^0) < 0.4$ GeV signal have $p_{\mathrm{T}}$ too low to be reconstructed as the signal candidate tracks, and therefore the identified signal candidate tracks are typically from pile-up collisions or underlying events similar to the QCD track background, causing a low $S(d_0)$ selection efficiency in these plots.

Event selection cutflows for signal samples with $m(\tilde{\chi}_{1}^0)$ = 150 GeV and $\Delta m(\tilde{\chi}_{1}^\pm, \tilde{\chi}_{1}^0)$ = 1.5, 1.0, and 0.75 GeV, including all six production processes ($\tilde{\chi}_1^\pm \tilde{\chi}_1^0$, $\tilde{\chi}_1^\pm \tilde{\chi}_2^0$, $\tilde{\chi}_1^+ \tilde{\chi}_1^-$, and $\tilde{\chi}_2^0 \tilde{\chi}_1^0$). The cross-section used to obtain the initial number of events ($\sigma(\mathrm{n}_{\mathrm{jets}}) \geq 1$) refers to an emission of at least one gluon or quark with $p_{\mathrm{T}} > 50$ GeV at the parton level.

Event selection cutflows for signal samples with $m(\tilde{\chi}_{1}^0)$ = 150 GeV and $\Delta m(\tilde{\chi}_{1}^\pm, \tilde{\chi}_{1}^0)$ = 0.5, 0.35, and 0.25 GeV, including all six production processes ($\tilde{\chi}_1^\pm \tilde{\chi}_1^0$, $\tilde{\chi}_1^\pm \tilde{\chi}_2^0$, $\tilde{\chi}_1^+ \tilde{\chi}_1^-$, and $\tilde{\chi}_2^0 \tilde{\chi}_1^0$). The cross-section used to obtain the initial number of events ($\sigma(\mathrm{n}_{\mathrm{jets}}) \geq 1$) refers to an emission of at least one gluon or quark with $p_{\mathrm{T}} > 50$ GeV at the parton level.

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

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

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

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

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

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

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

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

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

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

Measured min-$d_{xy}$ distribution in the 2-match category, after requiring the $d_{xy}$ muon to pass the isolation requirement $I_{\mathrm{rel}}^{\mathrm{PF}} <0.25$ (i. e., the B and D bins of the ABCD plane). Overlaid with a red histogram is the background predicted from the region of the ABCD plane failing the same requirement (the A and C bins), as well as three signal benchmark hypotheses (as defined in the legends), assuming $\alpha_D = \alpha_{\mathrm{EM}}$ (the fine-structure constant). The red hatched bands correspond to the background prediction uncertainty. The last bin includes the overflow.

Two-dimensional exclusion surfaces for $\Delta = 0.1 \, m_1$ as a function of the DM mass $m_1$ and the signal strength $y$, with $m_{A'} = 3 \, m_1$. Filled histograms denote observed limits on $\sigma(\mathrm{pp} \to A' \to \chi_2 \chi_1) \, \mathcal{B}(\chi_2 \to \chi_1 \mu^+ \mu^-)$. Solid (dashed) curves denote the observed (expected) exclusion limits at 95% CL, with 68% CL uncertainty bands around the expectation. Regions above the curves are excluded, depending on the $\alpha_D$ hypothesis: $\alpha_{\mathrm{D}} = \alpha_{\mathrm{EM}}$ (dark blue) or 0.1 (light magenta).

Two-dimensional exclusion surfaces for $\Delta = 0.4 \, m_1$ as a function of the DM mass $m_1$ and the signal strength $y$, with $m_{A'} = 3 \, m_1$. Filled histograms denote observed limits on $\sigma(\mathrm{pp} \to A' \to \chi_2 \chi_1) \, \mathcal{B}(\chi_2 \to \chi_1 \mu^+ \mu^-)$. Solid (dashed) curves denote the observed (expected) exclusion limits at 95% CL, with 68% CL uncertainty bands around the expectation. Regions above the curves are excluded, depending on the $\alpha_D$ hypothesis: $\alpha_{\mathrm{D}} = \alpha_{\mathrm{EM}}$ (dark blue) or 0.1 (light magenta).

Measured ratio of $\mathcal{B}_{4\mu}/\mathcal{B}_{2\mu}$

Measured branching fraction $\mathcal{B}_{4\mu}$

Measured branching fraction $\mathcal{B}_{4\mu}$

Comparison of the $m_{miss}$ shapes for the simulated signal events within the fiducial region and those outside it, after including the effect of PU protons as describe in the text, for a generated $m_{X}$ mass of 1000 GeV. The distributions are shown for single(+z)-single(−z) proton reconstruction categories.

Product of the acceptance times the combined reconstruction and identification efficiency, as a function of $m_{X}$ , for events generated inside the fiducial volume defined in Table 1. The curves shown panel display the different final states. The definition of the fiducial region and signal model used to estimate the acceptance is provided in the text.

Product of the acceptance times the combined reconstruction and identification efficiency, as a function of $m_{X}$ , for different proton reconstruction categories in the $Z \rightarrow \mu\mu $ analysis, and for events generated inside the fiducial volume defined in Table 1. The curves shown panel display the different final states. The definition of the fiducial region and signal model used to estimate the acceptance is provided in the text.

Distributions of the reconstructed proton $\xi$ in the negative arm from the proton mixing procedure with simulated MC events, are compared to data. The ee final states are shown.The ee events are displayed without the Z boson $p_T$ requirement.For illustration, the simulated signal distributions are superimposed for various choices of $m_{X}$ , normalised to a generated fiducial cross section of 100 pb.

Distributions of the reconstructed proton $\xi$ in the positive arm from the proton mixing procedure with simulated MC events, are compared to data. The ee final states are shown.The ee events are displayed without the Z boson $p_T$ requirement.For illustration, the simulated signal distributions are superimposed for various choices of $m_{X}$ , normalised to a generated fiducial cross section of 100 pb.

Distributions of the reconstructed di-proton rapidity from the proton mixing procedure with simulated MC events, are compared to data. The ee final states are shown.The ee events are displayed without the Z boson $p_T$ requirement.For illustration, the simulated signal distributions are superimposed for various choices of $m_{X}$ , normalised to a generated fiducial cross section of 100 pb.

Distributions of the reconstructed proton $\xi$ in the negative arm from the proton mixing procedure with simulated MC events, are compared to data. The $\mu\mu$ final states are shown.The $\mu\mu$ events are displayed without the Z boson $p_T$ requirement.For illustration, the simulated signal distributions are superimposed for various choices of $m_{X}$ , normalised to a generated fiducial cross section of 100 pb.

Distributions of the reconstructed proton $\xi$ in the positive arm from the proton mixing procedure with simulated MC events, are compared to data. The $\mu\mu$ final states are shown.The $\mu\mu$ events are displayed without the Z boson $p_T$ requirement.For illustration, the simulated signal distributions are superimposed for various choices of $m_{X}$ , normalised to a generated fiducial cross section of 100 pb.

Distributions of the reconstructed di-proton rapidity from the proton mixing procedure with simulated MC events, are compared to data. The $\mu\mu$ final states are shown.The $\mu\mu$ events are displayed without the Z boson $p_T$ requirement.For illustration, the simulated signal distributions are superimposed for various choices of $m_{X}$ , normalised to a generated fiducial cross section of 100 pb.

Distributions of the reconstructed proton $\xi$ in the negative arm from the proton mixing procedure with simulated MC events, are compared to data. The $\gamma$ final states are shown.For illustration, the simulated signal distributions are superimposed for various choices of $m_{X}$ , normalised to a generated fiducial cross section of 100 pb.

Distributions of the reconstructed proton $\xi$ in the positive arm from the proton mixing procedure with simulated MC events, are compared to data. The $\gamma$ final states are shown.For illustration, the simulated signal distributions are superimposed for various choices of $m_{X}$ , normalised to a generated fiducial cross section of 100 pb.

Distributions of the reconstructed di-proton rapidity from the proton mixing procedure with simulated MC events, are compared to data. The $\gamma$ final states are shown.For illustration, the simulated signal distributions are superimposed for various choices of $m_{X}$ , normalised to a generated fiducial cross section of 100 pb.

Validation of the background modelling method, using the $e\mu$ control sample. Selected $e\mu$ events are mixed with protons from $Z\rightarrow \mu \mu$ events with $p_{T}(Z)<10$ GeV to simulate the combinatorial background shape, while the data points are unaltered $e\mu$ events. The proton $\xi$ distribution for the positive CT-PPS arm is shown.

Validation of the background modelling method, using the $e\mu$ control sample. Selected $e\mu$ events are mixed with protons from $Z\rightarrow \mu \mu$ events with $p_{T}(Z)<10$ GeV to simulate the combinatorial background shape, while the data points are unaltered $e\mu$ events. The proton $\xi$ distribution for the negative CT-PPS arm is shown.

Validation of the background modelling method, using the $e\mu$ control sample. Selected $e\mu$ events are mixed with protons from $Z\rightarrow \mu \mu$ events with $p_{T}(Z)<10$ GeV to simulate the combinatorial background shape, while the data points are unaltered $e\mu$ events. The di-proton invariant mass is shown.

Validation of the background modelling method, using the $e\mu$ control sample. Selected $e\mu$ events are mixed with protons from $Z\rightarrow \mu \mu$ events with $p_{T}(Z)<10$ GeV to simulate the combinatorial background shape, while the data points are unaltered $e\mu$ events. The $m_{miss}$ is shown.

$m_{miss}$ distributions in the $pp\rightarrow ppZeeX$, with the protons reconstructed with the multi(+z)-multi(-z)method. The background distributions are shown after the fit. The expectations for a signal with $m_{X} = 1000$ GeV are superimposed, where the fiducial production cross section is normalised to 1 pb.

$m_{miss}$ distributions in the $pp\rightarrow ppZeeX$, with the protons reconstructed with the multi(+z)-single(-z)method. The background distributions are shown after the fit. The expectations for a signal with $m_{X} = 1000$ GeV are superimposed, where the fiducial production cross section is normalised to 1 pb.

$m_{miss}$ distributions in the $pp\rightarrow ppZeeX$, with the protons reconstructed with the single(+z)-multi(-z)method. The background distributions are shown after the fit. The expectations for a signal with $m_{X} = 1000$ GeV are superimposed, where the fiducial production cross section is normalised to 1 pb.

$m_{miss}$ distributions in the $pp\rightarrow ppZeeX$, with the protons reconstructed with the single(+z)-single(-z)method. The background distributions are shown after the fit. The expectations for a signal with $m_{X} = 1000$ GeV are superimposed, where the fiducial production cross section is normalised to 1 pb.

$m_{miss}$ distributions in the $pp\rightarrow ppZ\mu\mu X$, with the protons reconstructed with the multi(+z)-multi(-z)method. The background distributions are shown after the fit. The expectations for a signal with $m_{X} = 1000$ GeV are superimposed, where the fiducial production cross section is normalised to 1 pb.

$m_{miss}$ distributions in the $pp\rightarrow ppZ\mu\mu X$, with the protons reconstructed with the multi(+z)-single(-z)method. The background distributions are shown after the fit. The expectations for a signal with $m_{X} = 1000$ GeV are superimposed, where the fiducial production cross section is normalised to 1 pb.

$m_{miss}$ distributions in the $pp\rightarrow ppZ\mu\mu X$, with the protons reconstructed with the single(+z)-multi(-z)method. The background distributions are shown after the fit. The expectations for a signal with $m_{X} = 1000$ GeV are superimposed, where the fiducial production cross section is normalised to 1 pb.

$m_{miss}$ distributions in the $pp\rightarrow ppZ\mu\mu X$, with the protons reconstructed with the single(+z)-single(-z)method. The background distributions are shown after the fit. The expectations for a signal with $m_{X} = 1000$ GeV are superimposed, where the fiducial production cross section is normalised to 1 pb.

$m_{miss}$ distributions in the $pp\rightarrow pp\gamma X$, with the protons reconstructed with the multi(+z)-multi(-z)method. The background distributions are shown after the fit. The expectations for a signal with $m_{X} = 1000$ GeV are superimposed, where the fiducial production cross section is normalised to 10 pb.

$m_{miss}$ distributions in the $pp\rightarrow pp\gamma X$, with the protons reconstructed with the multi(+z)-single(-z)method. The background distributions are shown after the fit. The expectations for a signal with $m_{X} = 1000$ GeV are superimposed, where the fiducial production cross section is normalised to 10 pb.

$m_{miss}$ distributions in the $pp\rightarrow pp\gamma X$, with the protons reconstructed with the single(+z)-multi(-z)method. The background distributions are shown after the fit. The expectations for a signal with $m_{X} = 1000$ GeV are superimposed, where the fiducial production cross section is normalised to 10 pb.

$m_{miss}$ distributions in the $pp\rightarrow pp\gamma X$, with the protons reconstructed with the single(+z)-single(-z)method. The background distributions are shown after the fit. The expectations for a signal with $m_{X} = 1000$ GeV are superimposed, where the fiducial production cross section is normalised to 10 pb.

Upper limits on the $pp\rightarrow ppZ/\gamma X$ cross section at 95% CL, as a function of $m_X$. The 68 and 95% central intervals of the expected limits are represented by the green and yellow bands, respectively, while the observed limit is superimposed as a curve. The plot corresponds to the $Z\rightarrow ee$ analysis.

Upper limits on the $pp\rightarrow ppZ/\gamma X$ cross section at 95% CL, as a function of $m_X$. The 68 and 95% central intervals of the expected limits are represented by the green and yellow bands, respectively, while the observed limit is superimposed as a curve. The plot corresponds to the $Z\rightarrow \mu\mu$ analysis.

Upper limits on the $pp\rightarrow ppZ/\gamma X$ cross section at 95% CL, as a function of $m_X$. The 68 and 95% central intervals of the expected limits are represented by the green and yellow bands, respectively, while the observed limit is superimposed as a curve. The plot corresponds to the $Z$ analysis.

Upper limits on the $pp\rightarrow ppZ/\gamma X$ cross section at 95% CL, as a function of $m_X$. The 68 and 95% central intervals of the expected limits are represented by the green and yellow bands, respectively, while the observed limit is superimposed as a curve. The plot corresponds to the $\gamma$ analysis.

Two-dimensional distribution of the invariant mass $m_{DV}$ and the track multiplicity in the High-pT jet SR for observed data events

Two-dimensional distribution of the invariant mass $m_{DV}$ and the track multiplicity in the High-pT jet SR for expected signal events in the strong gluino pair pair production model with m(gluino)=1.8 TeV, m(chi0)=0.2 TeV, tau(chi0)=0.1 ns

Two-dimensional distribution of the invariant mass $m_{DV}$ and the track multiplicity in the Trackless jet SR for observed data events

Two-dimensional distribution of the invariant mass $m_{DV}$ and the track multiplicity in the Trackless jet SR for expected signal events in the electroweak pair production model

Expected exclusion limits at 95% CL on the lifetime and mass of the neutralino in electroweakino pair production models

Expected (+1 sigma) exclusion limits at 95% CL on the lifetime and mass of the neutralino in electroweakino pair production models

Expected (-1 sigma) exclusion limits at 95% CL on the lifetime and mass of the neutralino in electroweakino pair production models

Observed exclusion limits at 95% CL on the lifetime and mass of the neutralino in electroweakino pair production models

Observed (+1 sigma) exclusion limits at 95% CL on the lifetime and mass of the neutralino in electroweakino pair production models

Observed (-1 sigma) exclusion limits at 95% CL on the lifetime and mass of the neutralino in electroweakino pair production models

Expected exclusion limits at 95% CL on the lifetime and mass of the neutralino in strong gluino pair production models and m(gluino)=2.4 TeV

Expected (+1 sigma) exclusion limits at 95% CL on the lifetime and mass of the neutralino in strong gluino pair production models and m(gluino)=2.4 TeV

Expected (-1 sigma) exclusion limits at 95% CL on the lifetime and mass of the neutralino in strong gluino pair production models and m(gluino)=2.4 TeV

Observed exclusion limits at 95% CL on the lifetime and mass of the neutralino in strong gluino pair production models and m(gluino)=2.4 TeV

Observed (+1 sigma) exclusion limits at 95% CL on the lifetime and mass of the neutralino in strong gluino pair production models and m(gluino)=2.4 TeV

Observed (-1 sigma) exclusion limits at 95% CL on the lifetime and mass of the neutralino in strong gluino pair production models and m(gluino)=2.4 TeV

Exclusion limits at 95% CL on the production cross section in the electroweak pair production model.

Exclusion limits at 95% CL on the production cross section in the strong gluino pair production models and m(gluino)=2.4 TeV

Expected exclusion limits at 95% CL on the lifetime and mass of the neutralino in strong gluino pair production models and m(gluino)=2.0 TeV

Expected (+1 sigma) exclusion limits at 95% CL on the lifetime and mass of the neutralino in strong gluino pair production models and m(gluino)=2.0 TeV

Expected (-1 sigma) exclusion limits at 95% CL on the lifetime and mass of the neutralino in strong gluino pair production models and m(gluino)=2.0 TeV

Observed exclusion limits at 95% CL on the lifetime and mass of the neutralino in strong gluino pair production models and m(gluino)=2.0 TeV

Observed (+1 sigma) exclusion limits at 95% CL on the lifetime and mass of the neutralino in strong gluino pair production models and m(gluino)=2.0 TeV

Observed (-1 sigma) exclusion limits at 95% CL on the lifetime and mass of the neutralino in strong gluino pair production models and m(gluino)=2.0 TeV

Expected exclusion limits at 95% CL on the lifetime and mass of the neutralino in strong gluino pair production models and m(gluino)=2.2 TeV

Expected (+1 sigma) exclusion limits at 95% CL on the lifetime and mass of the neutralino in strong gluino pair production models and m(gluino)=2.2 TeV

Expected (-1 sigma) exclusion limits at 95% CL on the lifetime and mass of the neutralino in strong gluino pair production models and m(gluino)=2.2 TeV

Observed exclusion limits at 95% CL on the lifetime and mass of the neutralino in strong gluino pair production models and m(gluino)=2.2 TeV

Observed (+1 sigma) exclusion limits at 95% CL on the lifetime and mass of the neutralino in strong gluino pair production models and m(gluino)=2.2 TeV

Observed (-1 sigma) exclusion limits at 95% CL on the lifetime and mass of the neutralino in strong gluino pair production models and m(gluino)=2.2 TeV

Expected exclusion limits at 95% CL on the lifetime and mass of the gluino in strong gluino pair production models and m(chi0)=50 GeV

Expected (+1 sigma) exclusion limits at 95% CL on the lifetime and mass of the gluino in strong gluino pair production models and m(chi0)=50 GeV

Expected (-1 sigma) exclusion limits at 95% CL on the lifetime and mass of the gluino in strong gluino pair production models and m(chi0)=50 GeV

Observed exclusion limits at 95% CL on the lifetime and mass of the gluino in strong gluino pair production models and m(chi0)=50 GeV

Observed (+1 sigma) exclusion limits at 95% CL on the lifetime and mass of the gluino in strong gluino pair production models and m(chi0)=50 GeV

Observed (-1 sigma) exclusion limits at 95% CL on the lifetime and mass of the gluino in strong gluino pair production models and m(chi0)=50 GeV

Expected exclusion limits at 95% CL on the lifetime and mass of the gluino in strong gluino pair production models and m(chi0)=450 GeV

Expected (+1 sigma) exclusion limits at 95% CL on the lifetime and mass of the gluino in strong gluino pair production models and m(chi0)=450 GeV

Expected (-1 sigma) exclusion limits at 95% CL on the lifetime and mass of the gluino in strong gluino pair production models and m(chi0)=450 GeV

Observed exclusion limits at 95% CL on the lifetime and mass of the gluino in strong gluino pair production models and m(chi0)=450 GeV

Observed (+1 sigma) exclusion limits at 95% CL on the lifetime and mass of the gluino in strong gluino pair production models and m(chi0)=450 GeV

Observed (-1 sigma) exclusion limits at 95% CL on the lifetime and mass of the gluino in strong gluino pair production models and m(chi0)=450 GeV

Expected exclusion limits at 95% CL on the mass of the gluino and neutralino in strong gluino pair production models and tau(chi0)=0.01 ns

Expected (+1 sigma) exclusion limits at 95% CL on the mass of the gluino and neutralino in strong gluino pair production models and tau(chi0)=0.01 ns

Expected (-1 sigma) exclusion limits at 95% CL on the mass of the gluino and neutralino in strong gluino pair production models and tau(chi0)=0.01 ns

Observed exclusion limits at 95% CL on the mass of the gluino and neutralino in strong gluino pair production models and tau(chi0)=0.01 ns

Observed (+1 sigma) exclusion limits at 95% CL on the mass of the gluino and neutralino in strong gluino pair production models and tau(chi0)=0.01 ns

Observed (-1 sigma) exclusion limits at 95% CL on the mass of the gluino and neutralino in strong gluino pair production models and tau(chi0)=0.01 ns

Expected exclusion limits at 95% CL on the mass of the gluino and neutralino in strong gluino pair production models and tau(chi0)=0.1 ns

Expected (+1 sigma) exclusion limits at 95% CL on the mass of the gluino and neutralino in strong gluino pair production models and tau(chi0)=0.1 ns

Expected (-1 sigma) exclusion limits at 95% CL on the mass of the gluino and neutralino in strong gluino pair production models and tau(chi0)=0.1 ns

Observed exclusion limits at 95% CL on the mass of the gluino and neutralino in strong gluino pair production models and tau(chi0)=0.1 ns

Observed (+1 sigma) exclusion limits at 95% CL on the mass of the gluino and neutralino in strong gluino pair production models and tau(chi0)=0.1 ns

Observed (-1 sigma) exclusion limits at 95% CL on the mass of the gluino and neutralino in strong gluino pair production models and tau(chi0)=0.1 ns

Expected exclusion limits at 95% CL on the mass of the gluino and neutralino in strong gluino pair production models and tau(chi0)=1 ns

Expected (+1 sigma) exclusion limits at 95% CL on the mass of the gluino and neutralino in strong gluino pair production models and tau(chi0)=1 ns

Expected (-1 sigma) exclusion limits at 95% CL on the mass of the gluino and neutralino in strong gluino pair production models and tau(chi0)=1 ns

Observed exclusion limits at 95% CL on the mass of the gluino and neutralino in strong gluino pair production models and tau(chi0)=1 ns

Observed (+1 sigma) exclusion limits at 95% CL on the mass of the gluino and neutralino in strong gluino pair production models and tau(chi0)=1 ns

Observed (-1 sigma) exclusion limits at 95% CL on the mass of the gluino and neutralino in strong gluino pair production models and tau(chi0)=1 ns

Expected exclusion limits at 95% CL on the mass of the gluino and neutralino in strong gluino pair production models and tau(chi0)=10 ns

Expected (+1 sigma) exclusion limits at 95% CL on the mass of the gluino and neutralino in strong gluino pair production models and tau(chi0)=10 ns

Expected (-1 sigma) exclusion limits at 95% CL on the mass of the gluino and neutralino in strong gluino pair production models and tau(chi0)=10 ns

Observed exclusion limits at 95% CL on the mass of the gluino and neutralino in strong gluino pair production models and tau(chi0)=10 ns

Observed (+1 sigma) exclusion limits at 95% CL on the mass of the gluino and neutralino in strong gluino pair production models and tau(chi0)=10 ns

Observed (-1 sigma) exclusion limits at 95% CL on the mass of the gluino and neutralino in strong gluino pair production models and tau(chi0)=10 ns

Exclusion limits at 95% CL on the production cross section in the strong gluino pair production models and m($ ilde{\chi}^0_1$)=1.25 TeV

Acceptance cutflow for the High-pT SR for representative points in the strong gluino pair production model. See additional resources for more information.

Acceptance cutflow for the Trackless SR for representative points in the electroweak pair production model. See additional resources for more information.

Acceptance cutflow for the Trackless SR for representative points in the electroweak pair production model with heavy-flavor quarks final state. See additional resources for more information.

Acceptance cutflow for the High-pT SR for representative points in the electroweak pair production model with heavy-flavor quarks final state. See additional resources for more information.

Reinterpretation Material: Event-level Efficiency for HighPt SR selections, R < 1150 mm

Reinterpretation Material: Event-level Efficiency for HighPt SR selections, R [1150, 3870] mm

Reinterpretation Material: Event-level Efficiency for HighPt SR selections, R > 3870 mm

Reinterpretation Material: Event-level Efficiency for Trackless SR selections, R < 1150 mm

Reinterpretation Material: Event-level Efficiency for Trackless SR selections, R [1150, 3870] mm

Reinterpretation Material: Event-level Efficiency for Trackless SR selections, R > 3870 mm

Reinterpretation Material: Vertex-level Efficiency for R < 22 mm

Reinterpretation Material: Vertex-level Efficiency for R [22, 25] mm

Reinterpretation Material: Vertex-level Efficiency for R [25, 29] mm

Reinterpretation Material: Vertex-level Efficiency for R [29, 38] mm

Reinterpretation Material: Vertex-level Efficiency for R [38, 46] mm

Reinterpretation Material: Vertex-level Efficiency for R [46, 73] mm

Reinterpretation Material: Vertex-level Efficiency for R [73, 84] mm

Reinterpretation Material: Vertex-level Efficiency for R [84, 111] mm

Reinterpretation Material: Vertex-level Efficiency for R [111, 120] mm

Reinterpretation Material: Vertex-level Efficiency for R [120, 145] mm

Reinterpretation Material: Vertex-level Efficiency for R [145, 180] mm

Reinterpretation Material: Vertex-level Efficiency for R [180, 300] mm

Cutflow (acceptance x efficiency) for the High-pT SR for representative points in the strong gluino pair production model. See additional resources for more information.

Cutflow (acceptance x efficiency) for the Trackless SR for representative points in the electroweak pair production model. See additional resources for more information.

Cutflow (acceptance x efficiency) for the Trackless SR for representative points in the electroweak pair production model with heavy-flavor quarks. See additional resources for more information.

Cutflow (acceptance x efficiency) for the High-pT SR for representative points in the electroweak pair production model with heavy-flavor quarks. See additional resources for more information.

Electron $p_{T}$ distribution using 2016-2018 data sets. The complete set of selection criteria were applied. Two examples of SSM WPRIME boson signal samples with masses of 3.8 and 5.6 TeV are shown. The lower panel shows the ratio of data to SM prediction, and the shaded band represents the systematic uncertainties.

$p_{T}^{miss}$ distribution in the electron channel using 2016-2018 data sets. The complete set of selection criteria were applied. Two examples of SSM WPRIME boson signal samples with masses of 3.8 and 5.6 TeV are shown. The lower panel shows the ratio of data to SM prediction, and the shaded band represents the systematic uncertainties.

Muon $p_T$ distribution using 2016-2018 data sets. The complete set of selection criteria were applied. Two examples of SSM WPRIME boson signal samples with masses of 3.8 and 5.6 TeV are shown. The lower panel shows the ratio of data to SM prediction, and the shaded band represents the systematic uncertainties.

$p_T^{miss}$ distribution in the muon channel using 2016-2018 data sets. The complete set of selection criteria were applied. Two examples of SSM WPRIME boson signal samples with masses of 3.8 and 5.6 TeV are shown. The lower panel shows the ratio of data to SM prediction, and the shaded band represents the systematic uncertainties.

$M_T$ distribution for the E+INVISIBLE system using 2016-2018 data sets. The complete set of selection criteria were applied. Two examples of SSM WPRIME boson signal samples with masses of 3.8 and 5.6 TeV are shown. The lower panel shows the ratio of data to SM prediction, and the shaded band represents the systematic uncertainties.

$M_T$ distribution for the MU+INVISIBLE system using 2016-2018 data sets. The complete set of selection criteria were applied. Two examples of SSM WPRIME boson signal samples with masses of 3.8 and 5.6 TeV are shown. The lower panel shows the ratio of data to SM prediction, and the shaded band represents the systematic uncertainties.

Theoretical prediction for the production and decay cross-section of SSM WPRIME bosons (W --> LEPTON + NU) as a function of the WPRIME mass, for the range 200-6400 GeV and inclusive in phase-space (no cuts). The uncertainty shown covers PDF and ALPHAS uncertainties. These theoretical predictions are obtained with PYTHIA (LO) and corrected to NNLO in QCD using FEWZ program.

Observed (solid line) and expected (dashed line) upper limits at 95% CL on the product of WPRIME cross-section production and BF(WPRIME --> LEPTON + NU) for an SSM WPRIME boson, as a function of the boson mass, for the electron decay channel. The shaded bands represent the one (green) and two (yellow) standard deviation uncertainty bands for the expected limits. The theoretical predictions for the SSM at NNLO precision in QCD are shown with narrow gray bands corresponding to the PDF and ALPHAS uncertainties.

Observed (solid line) and expected (dashed line) upper limits at 95% CL on the product of WPRIME cross-section production and BF(WPRIME --> LEPTON + NU) for an SSM WPRIME boson, as a function of the boson mass, for the muon decay channel. The shaded bands represent the one (green) and two (yellow) standard deviation uncertainty bands for the expected limits. The theoretical predictions for the SSM at NNLO precision in QCD are shown with narrow gray bands corresponding to the PDF and ALPHAS uncertainties.

Observed (solid line) and expected (dashed line) upper limits at 95% CL on the product of WPRIME cross-section production and BF(WPRIME --> LEPTON + NU) for an SSM WPRIME boson, as a function of the boson mass, for the combined electron and muon decay channels. The shaded bands represent the one (green) and two (yellow) standard deviation uncertainty bands for the expected limits. The theoretical predictions for the SSM at NNLO precision in QCD are shown with narrow gray bands corresponding to the PDF and ALPHAS uncertainties.

The 95% CL observed (solid line) and expected (dashed line) model independent cross section upper limits as function of the $MT_{min}$ threshold, for the electron decay channel. The shaded bands represent the one (green) and two (yellow) standard deviation uncertainty bands for the expected limits.

The 95% CL observed (solid line) and expected (dashed line) model independent cross section upper limits as function of the $MT_{min}$ threshold, for the muon decay channel. The shaded bands represent the one (green) and two (yellow) standard deviation uncertainty bands for the expected limits.

The 95% CL observed (solid line) and expected (dashed line) model independent cross section upper limits as function of the $MT_{min}$ threshold, for the combined electron and muon decay channels. The shaded bands represent the one (green) and two (yellow) standard deviation uncertainty bands for the expected limits.

Observed (solid line) and expected (dashed line) upper limits at 95% CL on the coupling strength ratio, gWprime/gW as a function of the mass of the WPRIME boson, for the electron channel. The one (green) and two (yellow) standard deviation uncertainty bands for the expected limits are shown. The area above the limit curve is excluded. The dotted line represents the case of SSM coupling, gWprime, being equal to gW, the SM coupling.

Observed (solid line) and expected (dashed line) upper limits at 95% CL on the coupling strength ratio, gWprime/gW, as a function of the mass of the WPRIME boson, for the muon channel. The one (green) and two (yellow) standard deviation uncertainty bands for the expected limits are shown. The area above the limit curve is excluded. The dotted line represents the case of SSM coupling, gWprime , being equal to gW, the SM coupling.

Exclusion limits in the 2D plane (1/R, $\mu$) for the split-UED interpretation for the n = 2 case, for the electron channel for the 2016-2018 data sets. R is the radius of the additional spatial dimension and $\mu$ the bulk mass for the fermion field in five dimensions. The expected limit is depicted as a black dashed line. The one (green) and two (yellow) standard deviation uncertainty bands for the expected limits are shown. The experimentally excluded region is the entire area to the left of the solid black line.

Exclusion limits in the 2D plane (1/R, $\mu$) for the split-UED interpretation for the n = 2 case, for the muon channel for the 2016-2018 data sets. R is the radius of the additional spatial dimension and $\mu$ the bulk mass for the fermion field in five dimensions. The expected limit is depicted as a black dashed line. The one (green) and two (yellow) standard deviation uncertainty bands for the expected limits are shown. The experimentally excluded region is the entire area to the left of the solid black line.

Exclusion limits in the 2D plane (1/R, $\mu$) for the split-UED interpretation for the n = 2 case, for the combined electron and muon channels for the 2016-2018 data sets. R is the radius of the additional spatial dimension and $\mu$ the bulk mass for the fermion field in five dimensions. The expected limit is depicted as a black dashed line. The one (green) and two (yellow) standard deviation uncertainty bands for the expected limits are shown. The experimentally excluded region is the entire area to the left of the solid black line.

Observed (solid line) and expected (dashed line) upper limits at 95% CL on the ${\lambda}_{231}$ coupling in the RPV SUSY model with a STAU mediator, as a function of the mass of STAU, for the electron channel and for the production coupling, ${\lambda}^{'}_{3ij}$=0.5 . The couplings ${\lambda}^{'}_{3ij}$, ${\lambda}_{231}$, and ${\lambda}_{132}$ are defined in Fig. 1. The one (green) and two (yellow) standard deviation uncertainty bands for the expected limits are shown. The area above the limit curve is excluded. Limits for other values of ${\lambda}^{'}_{3ij}$ are obtained multiplying the ones shown by the ratio ${\lambda}^{'}_{3ij}$/0.5

Observed (solid line) and expected (dashed line) upper limits at 95% CL on the ${\lambda}_{132}$ coupling in the RPV SUSY model with a stau mediator, as a function of the mass of the stau, for the muon channel and for the production coupling, ${\lambda}^{'}_{3ij}$ = 0.5 . The couplings ${\lambda}^{'}_{3ij}$, ${\lambda}_{231}$, and ${\lambda}_{132}$ are defined in Fig. 1. The one (green) and two (yellow) standard deviation uncertainty bands for the expected limits are shown. The area above the limit curve is excluded. Limits for other values of ${\lambda}^{'}_{3ij}$ are obtained multiplying these ones by the ratio ${\lambda}^{'}_{3ij}$/0.5

Region in the oblique (W,Y) parameter phase space allowed by the current analysis at 95% CL. The combination of the electron and muon channel distributions from the 2017 and 2018 data sets at sqrt(s) = 13 TeV is used, having 101 fb-1 of luminosity.

Negative Log-Likelihood for the determination of the oblique W parameter, for the combination of the electron and muon channel and for 2017 and 2018 data sets (L=101 fb-1) at sqrt(s) = 13 TeV.

Total background estimate and observed data in all bins used in the limit calculation for the SSM WPRIME limits. The total event sample is split into the three years of data taking and the two lepton flavors.The distributions being shown are $M_T$ ones. The associated integrated luminosities for each year are of 35.9, 41.5, and 59.4 fb$̣^{-1}. The mass bins are numbered in ascending order.

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

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

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

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

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

Observed exclusion contour at 95&#37; C.L. for the dark Higgs model in the (m_Z', m_s) plane for g_q=0.25, g_&chi;=1, m_&chi;=200 GeV, and sin&theta;=0.01.

Expected exclusion contour at 95&#37; C.L. for the dark Higgs model in the (m_Z', m_s) plane for g_q=0.25, g_&chi;=1, m_&chi;=200 GeV, and sin&theta;=0.01.

Expected+1&sigma; exclusion contour at 95&#37; C.L. for the dark Higgs model in the (m_Z', m_s) plane for g_q=0.25, g_&chi;=1, m_&chi;=200 GeV, and sin&theta;=0.01.

Expected-1&sigma; exclusion contour at 95&#37; C.L. for the dark Higgs model in the (m_Z', m_s) plane for g_q=0.25, g_&chi;=1, m_&chi;=200 GeV, and sin&theta;=0.01.

Expected+2&sigma; exclusion contour at 95&#37; C.L. for the dark Higgs model in the (m_Z', m_s) plane for g_q=0.25, g_&chi;=1, m_&chi;=200 GeV, and sin&theta;=0.01.

Expected-2&sigma; exclusion contour at 95&#37; C.L. for the dark Higgs model in the (m_Z', m_s) plane for g_q=0.25, g_&chi;=1, m_&chi;=200 GeV, and sin&theta;=0.01.

Observed upper limits at 95&#37; C.L. on &sigma;(pp &rarr; s &chi;&chi;) &times; B(s &rarr; W^&pm; W^&mp;) for m_Z'=0.5 TeV as a function of m_s. The expected limits, varied up and down by one and two standard deviations, are shown as green and yellow bands, respectively. The observed and expected limits are compared to the theoretical LO cross section for the &sigma;(pp &rarr; s &chi;&chi;) &times; B(s &rarr; W^&pm; W^&mp;) process for m_Z'=0.5 TeV, shown in dashed blue.

Observed upper limits at 95&#37; C.L. on &sigma;(pp &rarr; s &chi;&chi;) &times; B(s &rarr; W^&pm; W^&mp;) for m_Z'=1 TeV as a function of m_s. The expected limits, varied up and down by one and two standard deviations, are shown as green and yellow bands, respectively. The observed and expected limits are compared to the theoretical LO cross section for the &sigma;(pp &rarr; s &chi;&chi;) &times; B(s &rarr; W^&pm; W^&mp;) process for m_Z'=1 TeV, shown in dashed blue.

Observed upper limits at 95&#37; C.L. on &sigma;(pp &rarr; s &chi;&chi;) &times; B(s &rarr; W^&pm; W^&mp;) for m_Z'=1.7 TeV as a function of m_s. The expected limits, varied up and down by one and two standard deviations, are shown as green and yellow bands, respectively. The observed and expected limits are compared to the theoretical LO cross section for the &sigma;(pp &rarr; s &chi;&chi;) &times; B(s &rarr; W^&pm; W^&mp;) process for m_Z'=1.7 TeV, shown in dashed blue.

Observed upper limits at 95&#37; C.L. on &sigma;(pp &rarr; s &chi;&chi;) &times; B(s &rarr; W^&pm; W^&mp;) for m_Z'=2.1 TeV as a function of m_s. The expected limits, varied up and down by one and two standard deviations, are shown as green and yellow bands, respectively. The observed and expected limits are compared to the theoretical LO cross section for the &sigma;(pp &rarr; s &chi;&chi;) &times; B(s &rarr; W^&pm; W^&mp;) process for m_Z'=2.1 TeV, shown in dashed blue.

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

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

Cumulative efficiencies in the merged category for three representative dark Higgs signal scenarios with g_q=0.25, g_&chi;=1.0, sin&theta;=0.01, m_Z' = 1 TeV, and m_&chi;=200 GeV considering s&rarr;W(&ell;&nu;)W(qq) decays only.

Cumulative efficiencies in the resolved category for three representative dark Higgs signal scenarios with g_q=0.25, g_&chi;=1.0, sin&theta;=0.01, m_Z' = 1 TeV, and m_&chi;=200 GeV considering s&rarr;W(&ell;&nu;)W(qq) decays only.

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

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

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

Parametrization of the $A_{X}$ factor, defined as the fraction of diphoton resonances satisfying the fiducial acceptance at the particle level, as function of $m_{X}$ obtained from the narrow width simulated signal samples produced in gluon fusion.

The correction factor, $C_{X}$, defined as the ratio between the number of reconstructed signal events passing the analysis cuts and the number of signal events at the particle level generated within the fiducial volume, and acceptance correction factor, $A_{X}$, defined as the fraction of diphoton resonances satisfying the fiducial acceptance at the particle level. Both are computed for NWA spin-0 models as a function of $m_{X}$.

Effect of event selections on a scalar MC signal sample generated for $m_{X}$ = 15 GeV and on the data. For the MC sample, the efficiencies are shown after applying event weights and a truth level filter that requires two photons with $p^{\gamma\gamma}_{T}>40$ GeV; for the data, the absolute yields are shown. The initial yields for data include a trigger preselection that is the OR of a list of single photon and diphoton triggers. The "2 $loose$ photons" step includes the kinematic acceptance cuts.

Parameterization of the Double Sided Crystal Ball function parameters describing the scalar mass resolution model as a function of $m_{X}$ [GeV].