Observed and expected 95% CL upper limits on the product of the production cross section ($\sigma$) and the branching fraction, obtained after combining all categories with 138 $\mathrm{fb}^{−1}$ of data at $\sqrt{s}$ = 13 TeV for $\mathrm{V'}$ to $VV$ + $VH$ signal in HVT model C.

Total signal efficiency including all event categories as a function of the resonance mass $M_X$ for all DY/ggF signal models under consideration. The estimates from MC simulation (markers) are interpolated with a functional form (lines).

Total signal efficiency including all event categories as a function of the resonance mass $M_X$ for all VBF signal models under consideration. The estimates from MC simulation (markers) are interpolated with a functional form (lines).

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

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

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.

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.

For the tu$\gamma$ signal channel, the total number of observed selected events in the data (N$_{obs}$), the SM expectations (N$_{SM}$), the efficiencies ($\epsilon$), and the upper limits on the cross sections $\sigma_{fid}$ at the 95% CL in the fiducial region, without and with a requirement on the presence of a single accompanying b jet.

For the tc$\gamma$ signal channel, the total number of observed selected events in the data (N$_{obs}$), the SM expectations (N$_{SM}$), the efficiencies ($\epsilon$), and the upper limits on the cross sections $\sigma_{fid}$ at the 95% CL in the fiducial region, without and with a requirement on the presence of a single accompanying b jet.

Signal selection efficiency as a function of mA.

Expected +2 sigma uncertainty band exclusion limit for Type-I models, as a function of tan beta and cos(beta-alpha).

Expected +1 sigma uncertainty band exclusion limit for Type-I models, as a function of tan beta and cos(beta-alpha).

Expected exclusion limit for Type-I models, as a function of tan beta and cos(beta-alpha).

Expected -1 sigma uncertainty band exclusion limit for Type-I models, as a function of tan beta and cos(beta-alpha).

Expected -2 sigma uncertainty band exclusion limit for Type-I models, as a function of tan beta and cos(beta-alpha).

Observed exclusion limit for Type-I models, as a function of tan beta and cos(beta-alpha).

Expected +2 sigma uncertainty band exclusion limit for Type-II models, as a function of tan beta and cos(beta-alpha).

Expected +1 sigma uncertainty band exclusion limit for Type-II models, as a function of tan beta and cos(beta-alpha).

Expected exclusion limit for Type-II models, as a function of tan beta and cos(beta-alpha).

Expected -1 sigma uncertainty band exclusion limit for Type-II models, as a function of tan beta and cos(beta-alpha).

Expected -2 sigma uncertainty band exclusion limit for Type-II models, as a function of tan beta and cos(beta-alpha).

Observed exclusion limit for Type-II models, as a function of tan beta and cos(beta-alpha).

Observed event yields and estimated backgrounds in the low- and high-mass control regions. The uncertainty in the background yield is the sum in quadrature of the statistical and systematic uncertainties.

Summary of the relative systematic uncertainties in heavy Majorana neutrino signal yields and the background from prompt same-sign leptons, both estimated from simulation. The relative systematic uncertainties assigned to the data-driven backgrounds for the misidentified lepton background and mismeasured charge background are also shown. The uncertainties are given for the low-mass selections.

Summary of the relative systematic uncertainties in heavy Majorana neutrino signal yields and the background from prompt same-sign leptons, both estimated from simulation. The relative systematic uncertainties assigned to the data-driven backgrounds for the misidentified lepton background and mismeasured charge background are also shown. The uncertainties are given for the high-mass selections.

Summary of contributions to the systematic uncertainty related to the prompt same-sign leptons background, misidentified lepton background, and mismeasured charge background on the total background uncertainty for the case of $m_{\rm N}$ = 100 and 500 GeV.

Observed event yields and estimated backgrounds after the application of all selection, except for the final optimization. The background predictions from prompt same-sign leptons (Prompt bkgd.), misidentified leptons (Misid. bkgd.), mismeasured charge (Charge mismeas. bkgd.) and the total background (Total bkgd.) are shown together withthe number of events observed in data. The uncertainties shown are the statistical and systematic components, respectively.

Dielectron channel results after final optimization. The background predictions from prompt same-sign leptons, misidentified leptons, and mismeasured charge are shown along with the total background estimate and the number of events observed in data. The uncertainties shown are the statistical and systematic components, respectively.

Electron-muon channel results after final optimization. The background predictions from prompt same-sign leptons and misidentified leptons are shown along with the total background estimate and the number of events observed in data. The uncertainties shown are the statistical and systematic components, respectively.

The product of signal acceptance and efficiency in the 2l categories for the signal produced via vector boson fusion.

$m_{X}^{T}$ distribution in data in the 0l 2b non-VBF category. The data is shown as Events/10 GeV. The distribution is shown up 4000 GeV, which corresponds to the event with the highest $m_{X}^{T}$ observed in the SR.

$m_{X}^{T}$ distribution in data in the 0l $\leq$1b non-VBF category. The data is shown as Events/10 GeV. The distribution is shown up 4000 GeV, which corresponds to the event with the highest $m_{X}^{T}$ observed in the SR.

$m_{X}$ distribution in data in the 2e 2b non-VBF category. The data is shown as Events/10 GeV. The distribution is shown up 4000 GeV, which corresponds to the event with the highest $m_{X}$ observed in the SR.

$m_{X}$ distribution in data in the 2e $\leq$1b non-VBF category. The data is shown as Events/10 GeV. The distribution is shown up 4000 GeV, which corresponds to the event with the highest $m_{X}$ observed in the SR.

$m_{X}$ distribution in data in the 2$\mu$ 2b non-VBF category. The data is shown as Events/10 GeV. The distribution is shown up 4000 GeV, which corresponds to the event with the highest $m_{X}$ observed in the SR.

$m_{X}$ distribution in data in the 2$\mu$ $\leq$1b non-VBF category. The data is shown as Events/10 GeV. The distribution is shown up 4000 GeV, which corresponds to the event with the highest $m_{X}$ observed in the SR.

$m_{X}^{T}$ distribution in data in the 0l 2b VBF category. The data is shown as Events/10 GeV. The distribution is shown up 4000 GeV, which corresponds to the event with the highest $m_{X}^{T}$ observed in the SR.

$m_{X}^{T}$ distribution in data in the 0l $\leq$1b VBF category. The data is shown as Events/10 GeV. The distribution is shown up 4000 GeV, which corresponds to the event with the highest $m_{X}^{T}$ observed in the SR.

$m_{X}$ distribution in data in the 2e 2b VBF category. The data is shown as Events/10 GeV. The distribution is shown up 4000 GeV, which corresponds to the event with the highest $m_{X}$ observed in the SR.

$m_{X}$ distribution in data in the 2e $\leq$1b VBF category. The data is shown as Events/10 GeV. The distribution is shown up 4000 GeV, which corresponds to the event with the highest $m_{X}$ observed in the SR.

$m_{X}$ distribution in data in the 2$\mu$ 2b VBF category. The data is shown as Events/10 GeV. The distribution is shown up 4000 GeV, which corresponds to the event with the highest $m_{X}$ observed in the SR.

$m_{X}$ distribution in data in the 2$\mu$ $\leq$1b VBF category. The data is shown as Events/10 GeV. The distribution is shown up 4000 GeV, which corresponds to the event with the highest $m_{X}$ observed in the SR.

Observed and expected 95% CL upper limit on $\sigma \mathcal{B}$(Z'-> ZH) with all categories combined for the non-VBF signal, including all statistical and systematic uncertainties. The inner green band and the outer yellow band indicate the regions containing 68 and 95%, respectively, of the distribution of expected limits under the background-only hypothesis. The CMS search for a heavy resonance using 2016 data and the same final state [JHEP 11 (2018) 172] is shown as a comparison.

Observed and expected 95% CL upper limit on $\sigma \mathcal{B}$(Z'-> ZH) with all categories combined for the VBF signal, including all statistical and systematic uncertainties. The inner green band and the outer yellow band indicate the regions containing 68 and 95%, respectively, of the distribution of expected limits under the background-only hypothesis.

Observed exclusion limit in the space of the HVT model parameters [$g_{V}c_{H}$, $g^{2}c_{F}/g_{V}$] for mass hypotheses of 2 TeV for the non-VBF signal.

Observed exclusion limit in the space of the HVT model parameters [$g_{V}c_{H}$, $g^{2}c_{F}/g_{V}$] for mass hypotheses of 3 TeV for the non-VBF signal.

Observed exclusion limit in the space of the HVT model parameters [$g_{V}c_{H}$, $g^{2}c_{F}/g_{V}$] for mass hypotheses of 4 TeV for the non-VBF signal.