Distributions of the number of hits in the cluster (Nhits) for the CSC category in the out-of-time (OOT) region. The last histogram bin contains all overflow events.

The 95% CL observed and expected upper limits on the VLL production cross section as a function of the VLL mass for the pseudoscalar mean proper decay length c$\tau$ = 0.025 m. The pseudoscalar mass is 2 GeV.

The 95% CL observed and expected upper limits on the VLL production cross section as a function of the pseudoscalar mean proper decay length c$\tau$ for VLL Mass of 700 GeV. The pseudoscalar mass is 2 GeV.

The 95% CL observed upper limits on the VLL production cross section as a function of the VLL mass and the pseudoscalar mean proper decay length c$\tau$. The pseudoscalar mass is 2 GeV.

The 95% CL observed exclusion contour on the VLL production cross section as a function of the VLL mass and the pseudoscalar mean proper decay length c$\tau$. The pseudoscalar mass is 2 GeV.

Selection efficiencies in the CSC category signal region (SR) for different signal hypothesis. The pseudoscalar mass is 2 GeV.

Selection efficiencies in the DT category signal region (SR) for different signal hypothesis. The pseudoscalar mass is 2 GeV.

The product of detector acceptance and analysis selection efficiency in the muon channel as functions of the particle X mass. Three analysis requirements are applied consecutively: event reconstruction, HLT, and final signal selection. The product of detector acceptance and analysis selection efficiencies are shown at each stage in red, blue, and orange, respectively.

Background-only fit to data. The lower panel contains the pull distribution, defined as the difference between the data yield and the background prediction divided by their combined uncertainty (electron channel).

Background-only fit to data. The lower panel contains the pull distribution, defined as the difference between the data yield and the background prediction divided by their combined uncertainty (muon channel).

The systematic uncertainties affecting signal normalization (electron channel).

The systematic uncertainties affecting signal normalization (muon channel).

The 95% CL expected and observed limits on P P --> X --> W gamma as a function of the resonant mass with leptonic decays of W bosons (narrow width).

The 95% CL expected and observed limits on P P --> X --> W gamma as a function of the resonant mass with leptonic decays of W bosons (broad width).

The 95% CL expected and observed limits on P P --> X --> W gamma as a function of the resonant mass with both leptonic and hadronic decays of W bosons (narrow width).

The 95% CL expected and observed limits on P P --> X --> W gamma as a function of the resonant mass with both leptonic and hadronic decays of W bosons (broad width).

The observed local p-values for narrow resonance width hypotheses.

The observed local p-values for broad resonance width hypotheses.

Fit results for the region 2 of pair produced resolved three jet mass ($m_{jjj}$) distribution, after $H_{\rm T}>600$ GeV, $|\eta|<2.4$, sixth jet $p_{\rm T}>50$ GeV, $D^2_{[(6,3)+(3,2)]}<1.0$, $A_m <0.175$, $\Delta>180$ GeV, $D^2_{[3,2]}<0.175$

Fit results for the region 3 of pair produced resolved three jet mass ($m_{jjj}$) distribution, after $H_{\rm T}>600$ GeV, $|\eta|<2.4$, sixth jet $p_{\rm T}>100$ GeV, $D^2_{[(6,3)+(3,2)]}<0.9$, $A_m <0.15$, $\Delta>150$ GeV, $D^2_{[3,2]}<0.2$

Expected and observed upper limits on the signal strength for the pair produced merged trijets. Acceptance after the selection $p_{\rm T}>300$ GeV, $|\eta|<2.4$, and $\tau_{32,\mathrm{DDT}}<0$ on both leading and subleading jet and $A_m<0.15$ between the two leading jets.

Expected and observed upper limits on the signal strength for the pair produced merged dijets. Acceptance after the selection $p_{\rm T}>300$ GeV, $|\eta|<2.4$, and $N^1_{2,\mathrm{DDT}}<0$ on both leading and subleading jet and $A_m<0.15$ between the two leading jets.

Expected and observed upper limits on the signal strength for the pair produced resolved trijets (including overlaps). Acceptance for three regions search requions described in table 1 of the paper.

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.

Feynman diagrams of $\mathrm{Z'}\to\mu^{-}\mu^{+}$ with a $\mathrm{Z'}$ boson produced via $\mathrm{s}\overline{\mathrm{b}}\to\mathrm{Z'}$, with one $\mathrm{b}$ quark and one $\overline{\mathrm{s}}$ quark in the final state.

Distribution of $\min(m_{\mu\mathrm{b}})$ as obtained from simulation in events with $N_{\mathrm{b}} \geq 1$ passing all the other selection requirements. In this search, we require $\min(m_{\mu\mathrm{b}})>175~\mathrm{GeV}$. The stacked histogram displays the expected distribution from the simulation of the SM backgrounds, while the overlaid open histograms illustrate the size and shape of the $\mathrm{Z'}$ contribution from the LFU model described in Eq. (1), for several $\mathrm{Z'}$ mass hypotheses. For illustrative purposes, we choose couplings $|g_{\ell}| = |g_{\nu}| = |g_{\mathrm{b}}| = 0.03$ and $\delta_{\mathrm{bs}}=0$. The contribution of background processes other than DY and $\mathrm{t}\overline{\mathrm{t}}$ is so small that it is only barely visible at the bottom of the stacked histogram. The hatched region indicates the statistical uncertainty arising from the limited size of the SM simulated samples. Histograms are normalized to unit area.

Distributions of $m_{\mu\mu}$ in the $N_{\mathrm{b}}=1$ event category. The stacked histogram displays the expected distribution from the SM background simulation. The overlaid open distributions illustrate the $\mathrm{Z'}$ contribution from the LFU model at $|g_{\ell}| = |g_{\nu}| = |g_{\mathrm{b}}| = 0.03$ and $\delta_{\mathrm{bs}}=0$ for a variety of $\mathrm{Z'}$ mass hypotheses. The observed data are shown as black points with statistical error bars. The hatched region indicates the statistical uncertainty arising from the limited size of the SM simulated samples. The size of the bins increases as a function of $m_{\mu\mu}$. In extracting the results of the search, the background is estimated directly from data, so the SM background simulation is only illustrative in these distributions.

Distributions of $m_{\mu\mu}$ in the $N_{\mathrm{b}}\geq 2$ event category. The stacked histogram displays the expected distribution from the SM background simulation. The overlaid open distributions illustrate the $\mathrm{Z'}$ contribution from the LFU model at $|g_{\ell}| = |g_{\nu}| = |g_{\mathrm{b}}| = 0.03$ and $\delta_{\mathrm{bs}}=0$ for a variety of $\mathrm{Z'}$ mass hypotheses. The observed data are shown as black points with statistical error bars. The hatched region indicates the statistical uncertainty arising from the limited size of the SM simulated samples. The size of the bins increases as a function of $m_{\mu\mu}$. In extracting the results of the search, the background is estimated directly from data, so the SM background simulation is only illustrative in these distributions.

Invariant mass $m_{\mu\mu}$ distributions in the $N_{\mathrm{b}} = 1$ category, shown together with the corresponding selected background functional forms used as input to the $discrete~profiling$ method when probing the $m_{\mathrm{Z'}}=500~\mathrm{GeV}$ hypothesis. The expected signal distribution for the LFU model described in Eq. (1), with couplings $|g_{\ell}| = |g_{\nu}| = |g_{\mathrm{b}}| = 0.03$ and $\delta_{\mathrm{bs}}=0$, is overlaid. The displayed mass range corresponds to the fit window used for this $m_{\mathrm{Z'}}$ hypothesis, which is $\pm 10\, \sigma_{\mathrm{mass}}$ around the probed $m_{\mathrm{Z'}}$ value. While the likelihood fits are performed on unbinned data, here we present the data in binned histograms with binning chosen to reflect the size of $\sigma_{\mathrm{mass}}$.

Invariant mass $m_{\mu\mu}$ distributions in the $N_{\mathrm{b}} \geq 2$ category, shown together with the corresponding selected background functional forms used as input to the $discrete~profiling$ method when probing the $m_{\mathrm{Z'}}=500~\mathrm{GeV}$ hypothesis. The expected signal distribution for the LFU model described in Eq. (1), with couplings $|g_{\ell}| = |g_{\nu}| = |g_{\mathrm{b}}| = 0.03$ and $\delta_{\mathrm{bs}}=0$, is overlaid. The displayed mass range corresponds to the fit window used for this $m_{\mathrm{Z'}}$ hypothesis, which is $\pm 10\, \sigma_{\mathrm{mass}}$ around the probed $m_{\mathrm{Z'}}$ value. While the likelihood fits are performed on unbinned data, here we present the data in binned histograms with binning chosen to reflect the size of $\sigma_{\mathrm{mass}}$.

Exclusion limits at $95\%$ CL on the number of selected BSM events with $N_{\mathrm{b}} \geq 1$ as functions of $m_{\mathrm{Z'}}$ for $f_{2\mathrm{b}}=0$. The quantity $f_{2\mathrm{b}}$ is the fraction of BSM events passing the analysis selection that have at least two $\mathrm{b}$ quark jets. The solid black (dashed red) curve represents the observed (median expected) exclusion. The inner green (outer yellow) band indicates the region containing $68~(95)\%$ of the distribution of limits expected under the background-only hypothesis.

Exclusion limits at $95\%$ CL on the number of selected BSM events with $N_{\mathrm{b}} \geq 1$ as functions of $m_{\mathrm{Z'}}$ for $f_{2\mathrm{b}}=0.25$. The quantity $f_{2\mathrm{b}}$ is the fraction of BSM events passing the analysis selection that have at least two $\mathrm{b}$ quark jets. The solid black (dashed red) curve represents the observed (median expected) exclusion. The inner green (outer yellow) band indicates the region containing $68~(95)\%$ of the distribution of limits expected under the background-only hypothesis.

Exclusion limits at $95\%$ CL on the number of selected BSM events with $N_{\mathrm{b}} \geq 1$ as functions of $m_{\mathrm{Z'}}$ for $f_{2\mathrm{b}}=0.75$. The quantity $f_{2\mathrm{b}}$ is the fraction of BSM events passing the analysis selection that have at least two $\mathrm{b}$ quark jets. The solid black (dashed red) curve represents the observed (median expected) exclusion. The inner green (outer yellow) band indicates the region containing $68~(95)\%$ of the distribution of limits expected under the background-only hypothesis.

Exclusion limits at $95\%$ CL on the number of selected BSM events with $N_{\mathrm{b}} \geq 1$ as functions of $m_{\mathrm{Z'}}$ for $f_{2\mathrm{b}}=1$. The quantity $f_{2\mathrm{b}}$ is the fraction of BSM events passing the analysis selection that have at least two $\mathrm{b}$ quark jets. The solid black (dashed red) curve represents the observed (median expected) exclusion. The inner green (outer yellow) band indicates the region containing $68~(95)\%$ of the distribution of limits expected under the background-only hypothesis.

Observed (solid) and median expected (dashed) exclusion limits at $95\%$ CL in the $|g_{\mathrm{b}}|$--$|g_{\ell}|$ plane for the LFU model, with $\delta_{\mathrm{bs}}=0$ and $|g_{\nu}|=|g_{\ell}|$. The exclusion limits are given up to coupling values at which the $\mathrm{Z'}$ width is equal to half of the $\mu\mu$ invariant mass resolution, marked by the dotted curves. Beyond these coupling values, the narrow width approximation intrinsic to the search strategy is not considered valid. The enclosed regions are excluded.

Observed exclusion limits at $95\%$ CL in the $|g_{\mathrm{b}}|$--$|g_{\ell}|$ plane for the LFU model, with $\delta_{\mathrm{bs}}=0$ and $|g_{\nu}|=|g_{\ell}|$, for $m_{\mathrm{Z'}}=350~\mathrm{GeV}$. The exclusion limits are given up to coupling values at which the $\mathrm{Z'}$ width is equal to half of the $\mu\mu$ invariant mass resolution. Beyond these coupling values, the narrow width approximation intrinsic to the search strategy is not considered valid. The enclosed regions are excluded.

Expected exclusion limits at $95\%$ CL in the $|g_{\mathrm{b}}|$--$|g_{\ell}|$ plane for the LFU model, with $\delta_{\mathrm{bs}}=0$ and $|g_{\nu}|=|g_{\ell}|$, for $m_{\mathrm{Z'}}=350~\mathrm{GeV}$. The exclusion limits are given up to coupling values at which the $\mathrm{Z'}$ width is equal to half of the $\mu\mu$ invariant mass resolution. Beyond these coupling values, the narrow width approximation intrinsic to the search strategy is not considered valid. The enclosed regions are excluded.

Exclusion limits at $95\%$ CL in the $|g_{\mathrm{b}}|$--$|g_{\ell}|$ plane for the LFU model, with $\delta_{\mathrm{bs}}=0$ and $|g_{\nu}|=|g_{\ell}|$, for $m_{\mathrm{Z'}}=350~\mathrm{GeV}$. The coupling values are shown at which the $\mathrm{Z'}$ width is equal to half of the $\mu\mu$ invariant mass resolution. Beyond these coupling values, the narrow width approximation intrinsic to the search strategy is not considered valid.

Observed exclusion limits at $95\%$ CL in the $|g_{\mathrm{b}}|$--$|g_{\ell}|$ plane for the LFU model, with $\delta_{\mathrm{bs}}=0$ and $|g_{\nu}|=|g_{\ell}|$, for $m_{\mathrm{Z'}}=500~\mathrm{GeV}$. The exclusion limits are given up to coupling values at which the $\mathrm{Z'}$ width is equal to half of the $\mu\mu$ invariant mass resolution. Beyond these coupling values, the narrow width approximation intrinsic to the search strategy is not considered valid. The enclosed regions are excluded.

Expected exclusion limits at $95\%$ CL in the $|g_{\mathrm{b}}|$--$|g_{\ell}|$ plane for the LFU model, with $\delta_{\mathrm{bs}}=0$ and $|g_{\nu}|=|g_{\ell}|$, for $m_{\mathrm{Z'}}=500~\mathrm{GeV}$. The exclusion limits are given up to coupling values at which the $\mathrm{Z'}$ width is equal to half of the $\mu\mu$ invariant mass resolution. Beyond these coupling values, the narrow width approximation intrinsic to the search strategy is not considered valid. The enclosed regions are excluded.

Exclusion limits at $95\%$ CL in the $|g_{\mathrm{b}}|$--$|g_{\ell}|$ plane for the LFU model, with $\delta_{\mathrm{bs}}=0$ and $|g_{\nu}|=|g_{\ell}|$, for $m_{\mathrm{Z'}}=500~\mathrm{GeV}$. The coupling values are shown at which the $\mathrm{Z'}$ width is equal to half of the $\mu\mu$ invariant mass resolution. Beyond these coupling values, the narrow width approximation intrinsic to the search strategy is not considered valid.

Observed exclusion limits at $95\%$ CL in the $|g_{\mathrm{b}}|$--$|g_{\ell}|$ plane for the LFU model, with $\delta_{\mathrm{bs}}=0$ and $|g_{\nu}|=|g_{\ell}|$, for $m_{\mathrm{Z'}}=700~\mathrm{GeV}$. The exclusion limits are given up to coupling values at which the $\mathrm{Z'}$ width is equal to half of the $\mu\mu$ invariant mass resolution. Beyond these coupling values, the narrow width approximation intrinsic to the search strategy is not considered valid. The enclosed regions are excluded.

Expected exclusion limits at $95\%$ CL in the $|g_{\mathrm{b}}|$--$|g_{\ell}|$ plane for the LFU model, with $\delta_{\mathrm{bs}}=0$ and $|g_{\nu}|=|g_{\ell}|$, for $m_{\mathrm{Z'}}=700~\mathrm{GeV}$. The exclusion limits are given up to coupling values at which the $\mathrm{Z'}$ width is equal to half of the $\mu\mu$ invariant mass resolution. Beyond these coupling values, the narrow width approximation intrinsic to the search strategy is not considered valid. The enclosed regions are excluded.

Exclusion limits at $95\%$ CL in the $|g_{\mathrm{b}}|$--$|g_{\ell}|$ plane for the LFU model, with $\delta_{\mathrm{bs}}=0$ and $|g_{\nu}|=|g_{\ell}|$, for $m_{\mathrm{Z'}}=700~\mathrm{GeV}$. The coupling values are shown at which the $\mathrm{Z'}$ width is equal to half of the $\mu\mu$ invariant mass resolution. Beyond these coupling values, the narrow width approximation intrinsic to the search strategy is not considered valid.

Observed exclusion limits at $95\%$ CL in the $|g_{\mathrm{b}}|$--$|g_{\ell}|$ plane for the LFU model, with $\delta_{\mathrm{bs}}=0$ and $|g_{\nu}|=|g_{\ell}|$, for $m_{\mathrm{Z'}}=1~\mathrm{TeV}$. The exclusion limits are given up to coupling values at which the $\mathrm{Z'}$ width is equal to half of the $\mu\mu$ invariant mass resolution. Beyond these coupling values, the narrow width approximation intrinsic to the search strategy is not considered valid. The enclosed regions are excluded.

Expected exclusion limits at $95\%$ CL in the $|g_{\mathrm{b}}|$--$|g_{\ell}|$ plane for the LFU model, with $\delta_{\mathrm{bs}}=0$ and $|g_{\nu}|=|g_{\ell}|$, for $m_{\mathrm{Z'}}=1~\mathrm{TeV}$. The exclusion limits are given up to coupling values at which the $\mathrm{Z'}$ width is equal to half of the $\mu\mu$ invariant mass resolution. Beyond these coupling values, the narrow width approximation intrinsic to the search strategy is not considered valid. The enclosed regions are excluded.

Exclusion limits at $95\%$ CL in the $|g_{\mathrm{b}}|$--$|g_{\ell}|$ plane for the LFU model, with $\delta_{\mathrm{bs}}=0$ and $|g_{\nu}|=|g_{\ell}|$, for $m_{\mathrm{Z'}}=1~\mathrm{TeV}$. The coupling values are shown at which the $\mathrm{Z'}$ width is equal to half of the $\mu\mu$ invariant mass resolution. Beyond these coupling values, the narrow width approximation intrinsic to the search strategy is not considered valid.

Observed exclusion limits at $95\%$ CL in the $|g_{\mathrm{b}}|$--$|g_{\ell}|$ plane for the LFU model, with $\delta_{\mathrm{bs}}=0$ and $|g_{\nu}|=|g_{\ell}|$, for $m_{\mathrm{Z'}}=1.5~\mathrm{TeV}$. The exclusion limits are given up to coupling values at which the $\mathrm{Z'}$ width is equal to half of the $\mu\mu$ invariant mass resolution. Beyond these coupling values, the narrow width approximation intrinsic to the search strategy is not considered valid. The enclosed regions are excluded.

Expected exclusion limits at $95\%$ CL in the $|g_{\mathrm{b}}|$--$|g_{\ell}|$ plane for the LFU model, with $\delta_{\mathrm{bs}}=0$ and $|g_{\nu}|=|g_{\ell}|$, for $m_{\mathrm{Z'}}=1.5~\mathrm{TeV}$. The exclusion limits are given up to coupling values at which the $\mathrm{Z'}$ width is equal to half of the $\mu\mu$ invariant mass resolution. Beyond these coupling values, the narrow width approximation intrinsic to the search strategy is not considered valid. The enclosed regions are excluded.

Exclusion limits at $95\%$ CL in the $|g_{\mathrm{b}}|$--$|g_{\ell}|$ plane for the LFU model, with $\delta_{\mathrm{bs}}=0$ and $|g_{\nu}|=|g_{\ell}|$, for $m_{\mathrm{Z'}}=1.5~\mathrm{TeV}$. The coupling values are shown at which the $\mathrm{Z'}$ width is equal to half of the $\mu\mu$ invariant mass resolution. Beyond these coupling values, the narrow width approximation intrinsic to the search strategy is not considered valid.

Observed exclusion limits at $95\%$ CL in the $|g_{\mathrm{b}}|$--$|g_{\ell}|$ plane for the LFU model, with $\delta_{\mathrm{bs}}=0$ and $|g_{\nu}|=|g_{\ell}|$, for $m_{\mathrm{Z'}}=2~\mathrm{TeV}$. The exclusion limits are given up to coupling values at which the $\mathrm{Z'}$ width is equal to half of the $\mu\mu$ invariant mass resolution. Beyond these coupling values, the narrow width approximation intrinsic to the search strategy is not considered valid. The enclosed regions are excluded.

Expected exclusion limits at $95\%$ CL in the $|g_{\mathrm{b}}|$--$|g_{\ell}|$ plane for the LFU model, with $\delta_{\mathrm{bs}}=0$ and $|g_{\nu}|=|g_{\ell}|$, for $m_{\mathrm{Z'}}=2~\mathrm{TeV}$. The exclusion limits are given up to coupling values at which the $\mathrm{Z'}$ width is equal to half of the $\mu\mu$ invariant mass resolution. Beyond these coupling values, the narrow width approximation intrinsic to the search strategy is not considered valid. The enclosed regions are excluded.

Exclusion limits at $95\%$ CL in the $|g_{\mathrm{b}}|$--$|g_{\ell}|$ plane for the LFU model, with $\delta_{\mathrm{bs}}=0$ and $|g_{\nu}|=|g_{\ell}|$, for $m_{\mathrm{Z'}}=2~\mathrm{TeV}$. The coupling values are shown at which the $\mathrm{Z'}$ width is equal to half of the $\mu\mu$ invariant mass resolution. Beyond these coupling values, the narrow width approximation intrinsic to the search strategy is not considered valid.

Observed (solid) and median expected (dashed) exclusion limits at $95\%$ CL in the $|g_{\mathrm{b}}|$--$|g_{\ell}|$ plane for the LFU model, with $|\delta_{\mathrm{bs}}|=0.25$ and $|g_{\nu}|=|g_{\ell}|$. The exclusion limits are given up to coupling values at which the $\mathrm{Z'}$ width is equal to half of the $\mu\mu$ invariant mass resolution, marked by the dotted curves. Beyond these coupling values, the narrow width approximation intrinsic to the search strategy is not considered valid. The enclosed regions are excluded.

Observed exclusion limits at $95\%$ CL in the $|g_{\mathrm{b}}|$--$|g_{\ell}|$ plane for the LFU model, with $|\delta_{\mathrm{bs}}|=0.25$ and $|g_{\nu}|=|g_{\ell}|$, for $m_{\mathrm{Z'}}=350~\mathrm{GeV}$. The exclusion limits are given up to coupling values at which the $\mathrm{Z'}$ width is equal to half of the $\mu\mu$ invariant mass resolution. Beyond these coupling values, the narrow width approximation intrinsic to the search strategy is not considered valid. The enclosed regions are excluded.

Expected exclusion limits at $95\%$ CL in the $|g_{\mathrm{b}}|$--$|g_{\ell}|$ plane for the LFU model, with $|\delta_{\mathrm{bs}}|=0.25$ and $|g_{\nu}|=|g_{\ell}|$, for $m_{\mathrm{Z'}}=350~\mathrm{GeV}$. The exclusion limits are given up to coupling values at which the $\mathrm{Z'}$ width is equal to half of the $\mu\mu$ invariant mass resolution. Beyond these coupling values, the narrow width approximation intrinsic to the search strategy is not considered valid. The enclosed regions are excluded.

Exclusion limits at $95\%$ CL in the $|g_{\mathrm{b}}|$--$|g_{\ell}|$ plane for the LFU model, with $|\delta_{\mathrm{bs}}|=0.25$ and $|g_{\nu}|=|g_{\ell}|$, for $m_{\mathrm{Z'}}=350~\mathrm{GeV}$. The coupling values are shown at which the $\mathrm{Z'}$ width is equal to half of the $\mu\mu$ invariant mass resolution. Beyond these coupling values, the narrow width approximation intrinsic to the search strategy is not considered valid.

Observed exclusion limits at $95\%$ CL in the $|g_{\mathrm{b}}|$--$|g_{\ell}|$ plane for the LFU model, with $|\delta_{\mathrm{bs}}|=0.25$ and $|g_{\nu}|=|g_{\ell}|$, for $m_{\mathrm{Z'}}=500~\mathrm{GeV}$. The exclusion limits are given up to coupling values at which the $\mathrm{Z'}$ width is equal to half of the $\mu\mu$ invariant mass resolution. Beyond these coupling values, the narrow width approximation intrinsic to the search strategy is not considered valid. The enclosed regions are excluded.

Expected exclusion limits at $95\%$ CL in the $|g_{\mathrm{b}}|$--$|g_{\ell}|$ plane for the LFU model, with $|\delta_{\mathrm{bs}}|=0.25$ and $|g_{\nu}|=|g_{\ell}|$, for $m_{\mathrm{Z'}}=500~\mathrm{GeV}$. The exclusion limits are given up to coupling values at which the $\mathrm{Z'}$ width is equal to half of the $\mu\mu$ invariant mass resolution. Beyond these coupling values, the narrow width approximation intrinsic to the search strategy is not considered valid. The enclosed regions are excluded.

Exclusion limits at $95\%$ CL in the $|g_{\mathrm{b}}|$--$|g_{\ell}|$ plane for the LFU model, with $|\delta_{\mathrm{bs}}|=0.25$ and $|g_{\nu}|=|g_{\ell}|$, for $m_{\mathrm{Z'}}=500~\mathrm{GeV}$. The coupling values are shown at which the $\mathrm{Z'}$ width is equal to half of the $\mu\mu$ invariant mass resolution. Beyond these coupling values, the narrow width approximation intrinsic to the search strategy is not considered valid.

Observed exclusion limits at $95\%$ CL in the $|g_{\mathrm{b}}|$--$|g_{\ell}|$ plane for the LFU model, with $|\delta_{\mathrm{bs}}|=0.25$ and $|g_{\nu}|=|g_{\ell}|$, for $m_{\mathrm{Z'}}=700~\mathrm{GeV}$. The exclusion limits are given up to coupling values at which the $\mathrm{Z'}$ width is equal to half of the $\mu\mu$ invariant mass resolution. Beyond these coupling values, the narrow width approximation intrinsic to the search strategy is not considered valid. The enclosed regions are excluded.

Expected exclusion limits at $95\%$ CL in the $|g_{\mathrm{b}}|$--$|g_{\ell}|$ plane for the LFU model, with $|\delta_{\mathrm{bs}}|=0.25$ and $|g_{\nu}|=|g_{\ell}|$, for $m_{\mathrm{Z'}}=700~\mathrm{GeV}$. The exclusion limits are given up to coupling values at which the $\mathrm{Z'}$ width is equal to half of the $\mu\mu$ invariant mass resolution. Beyond these coupling values, the narrow width approximation intrinsic to the search strategy is not considered valid. The enclosed regions are excluded.

Exclusion limits at $95\%$ CL in the $|g_{\mathrm{b}}|$--$|g_{\ell}|$ plane for the LFU model, with $|\delta_{\mathrm{bs}}|=0.25$ and $|g_{\nu}|=|g_{\ell}|$, for $m_{\mathrm{Z'}}=700~\mathrm{GeV}$. The coupling values are shown at which the $\mathrm{Z'}$ width is equal to half of the $\mu\mu$ invariant mass resolution. Beyond these coupling values, the narrow width approximation intrinsic to the search strategy is not considered valid.

Observed exclusion limits at $95\%$ CL in the $|g_{\mathrm{b}}|$--$|g_{\ell}|$ plane for the LFU model, with $|\delta_{\mathrm{bs}}|=0.25$ and $|g_{\nu}|=|g_{\ell}|$, for $m_{\mathrm{Z'}}=1~\mathrm{TeV}$. The exclusion limits are given up to coupling values at which the $\mathrm{Z'}$ width is equal to half of the $\mu\mu$ invariant mass resolution. Beyond these coupling values, the narrow width approximation intrinsic to the search strategy is not considered valid. The enclosed regions are excluded.

Expected exclusion limits at $95\%$ CL in the $|g_{\mathrm{b}}|$--$|g_{\ell}|$ plane for the LFU model, with $|\delta_{\mathrm{bs}}|=0.25$ and $|g_{\nu}|=|g_{\ell}|$, for $m_{\mathrm{Z'}}=1~\mathrm{TeV}$. The exclusion limits are given up to coupling values at which the $\mathrm{Z'}$ width is equal to half of the $\mu\mu$ invariant mass resolution. Beyond these coupling values, the narrow width approximation intrinsic to the search strategy is not considered valid. The enclosed regions are excluded.

Exclusion limits at $95\%$ CL in the $|g_{\mathrm{b}}|$--$|g_{\ell}|$ plane for the LFU model, with $|\delta_{\mathrm{bs}}|=0.25$ and $|g_{\nu}|=|g_{\ell}|$, for $m_{\mathrm{Z'}}=1~\mathrm{TeV}$. The coupling values are shown at which the $\mathrm{Z'}$ width is equal to half of the $\mu\mu$ invariant mass resolution. Beyond these coupling values, the narrow width approximation intrinsic to the search strategy is not considered valid.

Observed exclusion limits at $95\%$ CL in the $|g_{\mathrm{b}}|$--$|g_{\ell}|$ plane for the LFU model, with $|\delta_{\mathrm{bs}}|=0.25$ and $|g_{\nu}|=|g_{\ell}|$, for $m_{\mathrm{Z'}}=1.5~\mathrm{TeV}$. The exclusion limits are given up to coupling values at which the $\mathrm{Z'}$ width is equal to half of the $\mu\mu$ invariant mass resolution. Beyond these coupling values, the narrow width approximation intrinsic to the search strategy is not considered valid. The enclosed regions are excluded.

Expected exclusion limits at $95\%$ CL in the $|g_{\mathrm{b}}|$--$|g_{\ell}|$ plane for the LFU model, with $|\delta_{\mathrm{bs}}|=0.25$ and $|g_{\nu}|=|g_{\ell}|$, for $m_{\mathrm{Z'}}=1.5~\mathrm{TeV}$. The exclusion limits are given up to coupling values at which the $\mathrm{Z'}$ width is equal to half of the $\mu\mu$ invariant mass resolution. Beyond these coupling values, the narrow width approximation intrinsic to the search strategy is not considered valid. The enclosed regions are excluded.

Exclusion limits at $95\%$ CL in the $|g_{\mathrm{b}}|$--$|g_{\ell}|$ plane for the LFU model, with $|\delta_{\mathrm{bs}}|=0.25$ and $|g_{\nu}|=|g_{\ell}|$, for $m_{\mathrm{Z'}}=1.5~\mathrm{TeV}$. The coupling values are shown at which the $\mathrm{Z'}$ width is equal to half of the $\mu\mu$ invariant mass resolution. Beyond these coupling values, the narrow width approximation intrinsic to the search strategy is not considered valid.

Observed exclusion limits at $95\%$ CL in the $|g_{\mathrm{b}}|$--$|g_{\ell}|$ plane for the LFU model, with $|\delta_{\mathrm{bs}}|=0.25$ and $|g_{\nu}|=|g_{\ell}|$, for $m_{\mathrm{Z'}}=2~\mathrm{TeV}$. The exclusion limits are given up to coupling values at which the $\mathrm{Z'}$ width is equal to half of the $\mu\mu$ invariant mass resolution. Beyond these coupling values, the narrow width approximation intrinsic to the search strategy is not considered valid. The enclosed regions are excluded.

Exclusion limits at $95\%$ CL in the $|g_{\mathrm{b}}|$--$|g_{\ell}|$ plane for the LFU model, with $|\delta_{\mathrm{bs}}|=0.25$ and $|g_{\nu}|=|g_{\ell}|$, for $m_{\mathrm{Z'}}=2~\mathrm{TeV}$. The coupling values are shown at which the $\mathrm{Z'}$ width is equal to half of the $\mu\mu$ invariant mass resolution. Beyond these coupling values, the narrow width approximation intrinsic to the search strategy is not considered valid.

Exclusion limits at $95\%$ CL in the $|\theta_{23}|$--$|g_{\mathrm{Z'}}|$ plane for the $B_{3}\!-\!L_{2}$ model from JHEP 04 (2023) 033, for $m_{\mathrm{Z'}}=500~\mathrm{GeV}$. The solid black (dashed red) curves represent the observed (median expected) exclusions. The dotted curves denote the coupling values at which the $\mathrm{Z'}$ width equals one half of the $\mu\mu$ invariant mass resolution. For larger values of the couplings, the narrow width approximation intrinsic to the search strategy is not considered valid. The region enclosed between the solid black (dashed red) and dotted curves is (expected to be) excluded. The dotted curve for $m_{\mathrm{Z'}}=500~\mathrm{GeV}$ lies beyond the displayed $|g_{\mathrm{Z'}}|$ range and is, therefore, not shown. The shaded blue area represents the region preferred from the global fit in JHEP 04 (2023) 033 at $95\%$ CL. The region above the green dash-dotted curve is incompatible at $95\%$ CL with the measurement of the mass difference between the mass eigenstates of the neutral $\mathrm{B}_\mathrm{s}$ mesons.

Observed exclusion limits at $95\%$ CL in the $|\theta_{23}|$--$|g_{\mathrm{Z'}}|$ plane for the $B_{3}\!-\!L_{2}$ model from JHEP 04 (2023) 033, for $m_{\mathrm{Z'}}=500~\mathrm{GeV}$.

Expected exclusion limits at $95\%$ CL in the $|\theta_{23}|$--$|g_{\mathrm{Z'}}|$ plane for the $B_{3}\!-\!L_{2}$ model from JHEP 04 (2023) 033, for $m_{\mathrm{Z'}}=500~\mathrm{GeV}$.

Exclusion limits at $95\%$ CL in the $|\theta_{23}|$--$|g_{\mathrm{Z'}}|$ plane for the $B_{3}\!-\!L_{2}$ model from JHEP 04 (2023) 033, for $m_{\mathrm{Z'}}=1~\mathrm{TeV}$. The solid black (dashed red) curves represent the observed (median expected) exclusions. The dotted curves denote the coupling values at which the $\mathrm{Z'}$ width equals one half of the $\mu\mu$ invariant mass resolution. For larger values of the couplings, the narrow width approximation intrinsic to the search strategy is not considered valid. The region enclosed between the solid black (dashed red) and dotted curves is (expected to be) excluded. The shaded blue area represents the region preferred from the global fit in JHEP 04 (2023) 033 at $95\%$ CL. The region above the green dash-dotted curve is incompatible at $95\%$ CL with the measurement of the mass difference between the mass eigenstates of the neutral $\mathrm{B}_\mathrm{s}$ mesons.

Observed exclusion limits at $95\%$ CL in the $|\theta_{23}|$--$|g_{\mathrm{Z'}}|$ plane for the $B_{3}\!-\!L_{2}$ model from JHEP 04 (2023) 033, for $m_{\mathrm{Z'}}=1~\mathrm{TeV}$.

Expected exclusion limits at $95\%$ CL in the $|\theta_{23}|$--$|g_{\mathrm{Z'}}|$ plane for the $B_{3}\!-\!L_{2}$ model from JHEP 04 (2023) 033, for $m_{\mathrm{Z'}}=1~\mathrm{TeV}$.

Exclusion limits at $95\%$ CL in the $|\theta_{23}|$--$|g_{\mathrm{Z'}}|$ plane for the $B_{3}\!-\!L_{2}$ model from JHEP 04 (2023) 033, for $m_{\mathrm{Z'}}=1~\mathrm{TeV}$. The coupling values are shown at which the $\mathrm{Z'}$ width equals one half of the $\mu\mu$ invariant mass resolution. For larger values of the couplings, the narrow width approximation intrinsic to the search strategy is not considered valid.

Exclusion limits at $95\%$ CL in the $|\theta_{23}|$--$|g_{\mathrm{Z'}}|$ plane for the $B_{3}\!-\!L_{2}$ model from JHEP 04 (2023) 033, for $m_{\mathrm{Z'}}=1.5~\mathrm{TeV}$. The solid black (dashed red) curves represent the observed (median expected) exclusions. The dotted curves denote the coupling values at which the $\mathrm{Z'}$ width equals one half of the $\mu\mu$ invariant mass resolution. For larger values of the couplings, the narrow width approximation intrinsic to the search strategy is not considered valid. The region enclosed between the solid black (dashed red) and dotted curves is (expected to be) excluded. The shaded blue area represents the region preferred from the global fit in JHEP 04 (2023) 033 at $95\%$ CL. The region above the green dash-dotted curve is incompatible at $95\%$ CL with the measurement of the mass difference between the mass eigenstates of the neutral $\mathrm{B}_\mathrm{s}$ mesons.

Observed exclusion limits at $95\%$ CL in the $|\theta_{23}|$--$|g_{\mathrm{Z'}}|$ plane for the $B_{3}\!-\!L_{2}$ model from JHEP 04 (2023) 033, for $m_{\mathrm{Z'}}=1.5~\mathrm{TeV}$.

Expected exclusion limits at $95\%$ CL in the $|\theta_{23}|$--$|g_{\mathrm{Z'}}|$ plane for the $B_{3}\!-\!L_{2}$ model from JHEP 04 (2023) 033, for $m_{\mathrm{Z'}}=1.5~\mathrm{TeV}$.

Exclusion limits at $95\%$ CL in the $|\theta_{23}|$--$|g_{\mathrm{Z'}}|$ plane for the $B_{3}\!-\!L_{2}$ model from JHEP 04 (2023) 033, for $m_{\mathrm{Z'}}=1.5~\mathrm{TeV}$. The coupling values are shown at which the $\mathrm{Z'}$ width equals one half of the $\mu\mu$ invariant mass resolution. For larger values of the couplings, the narrow width approximation intrinsic to the search strategy is not considered valid.

Exclusion limits at $95\%$ CL in the $|\theta_{23}|$--$|g_{\mathrm{Z'}}|$ plane for the $B_{3}\!-\!L_{2}$ model from JHEP 04 (2023) 033, for $m_{\mathrm{Z'}}=2~\mathrm{TeV}$. The solid black (dashed red) curves represent the observed (median expected) exclusions. The dotted curves denote the coupling values at which the $\mathrm{Z'}$ width equals one half of the $\mu\mu$ invariant mass resolution. For larger values of the couplings, the narrow width approximation intrinsic to the search strategy is not considered valid. The region enclosed between the solid black (dashed red) and dotted curves is (expected to be) excluded. The shaded blue area represents the region preferred from the global fit in JHEP 04 (2023) 033 at $95\%$ CL. The region above the green dash-dotted curve is incompatible at $95\%$ CL with the measurement of the mass difference between the mass eigenstates of the neutral $\mathrm{B}_\mathrm{s}$ mesons.

Observed exclusion limits at $95\%$ CL in the $|\theta_{23}|$--$|g_{\mathrm{Z'}}|$ plane for the $B_{3}\!-\!L_{2}$ model from JHEP 04 (2023) 033, for $m_{\mathrm{Z'}}=2~\mathrm{TeV}$.

Expected exclusion limits at $95\%$ CL in the $|\theta_{23}|$--$|g_{\mathrm{Z'}}|$ plane for the $B_{3}\!-\!L_{2}$ model from JHEP 04 (2023) 033, for $m_{\mathrm{Z'}}=2~\mathrm{TeV}$.

Exclusion limits at $95\%$ CL in the $|\theta_{23}|$--$|g_{\mathrm{Z'}}|$ plane for the $B_{3}\!-\!L_{2}$ model from JHEP 04 (2023) 033, for $m_{\mathrm{Z'}}=2~\mathrm{TeV}$. The coupling values are shown at which the $\mathrm{Z'}$ width equals one half of the $\mu\mu$ invariant mass resolution. For larger values of the couplings, the narrow width approximation intrinsic to the search strategy is not considered valid.

Exclusion limits at $95\%$ CL in the $|g_{\mathrm{Z'}}|$--$m_{\mathrm{Z'}}$ plane for the $B_{3}\!-\!L_{2}$ model from JHEP 04 (2023) 033, for a fixed value of $\theta_{23}=0$. The solid black (dashed red) curve represents the observed (median expected) exclusion. The dotted curve denotes the coupling values at which the $\mathrm{Z'}$ width equals one half of the $\mu\mu$ invariant mass resolution. For larger values of the couplings, the narrow width approximation intrinsic to the search strategy is not considered valid. The region enclosed between the solid black (dashed red) and the dotted curves is (expected to be) excluded. The shaded blue area represents the $|g_{\mathrm{Z'}}|$ range preferred from a global fit in JHEP 04 (2023) 033 at $95\%$ CL.

Observed exclusion limits at $95\%$ CL in the $|g_{\mathrm{Z'}}|$--$m_{\mathrm{Z'}}$ plane for the $B_{3}\!-\!L_{2}$ model from JHEP 04 (2023) 033, for a fixed value of $\theta_{23}=0$.

Expected exclusion limits at $95\%$ CL in the $|g_{\mathrm{Z'}}|$--$m_{\mathrm{Z'}}$ plane for the $B_{3}\!-\!L_{2}$ model from JHEP 04 (2023) 033, for a fixed value of $\theta_{23}=0$.

Exclusion limits at $95\%$ CL in the $|g_{\mathrm{Z'}}|$--$m_{\mathrm{Z'}}$ plane for the $B_{3}\!-\!L_{2}$ model from JHEP 04 (2023) 033, for a fixed value of $\theta_{23}=0$. The coupling values are shown at which the $\mathrm{Z'}$ width equals one half of the $\mu\mu$ invariant mass resolution. For larger values of the couplings, the narrow width approximation intrinsic to the search strategy is not considered valid.

The product of acceptance and selection efficiency within the fiducial region for the VBF $\mathrm{H}^{\pm\pm}\rightarrow\mathrm{W}^{\pm}\mathrm{W}^{\pm}\rightarrow 2\ell 2v$ and $\mathrm{H}^{\pm}\rightarrow\mathrm{W}^{\pm}\mathrm{Z}\rightarrow 3\ell v$ processes, as a function of $m_{\mathrm{H}_{5}}$. The combination of the statistical and systematic uncertainties is shown. The theoretical uncertainties in the acceptance are also included.

Expected and observed exclusion limits at 95\% confidence level for $\sigma_\mathrm{VBF}(\mathrm{H}^{\pm\pm}) \mathcal{B}(\mathrm{H}^{\pm\pm} \rightarrow \mathrm{W}^{\pm}\mathrm{W}^{\pm})$ as functions of $m_{\mathrm{H}^{\pm\pm}}$. The contribution of the $\mathrm{H}^{\pm}$ boson signal is set to zero for the derivation of the exclusion limits on the $\sigma_\mathrm{VBF}(\mathrm{H}^{\pm\pm}) \mathcal{B}(\mathrm{H}^{\pm\pm} \rightarrow \mathrm{W}^{\pm}\mathrm{W}^{\pm})$. Values above the curves are excluded.

Expected and observed exclusion limits at 95\% confidence level for $\sigma_\mathrm{VBF}(\mathrm{H}^{\pm}) \mathcal{B}(\mathrm{H}^{\pm} \rightarrow \mathrm{W}^{\pm}\mathrm{Z})$ as functions of $m_{\mathrm{H}^{\pm}}$. The contribution of the $\mathrm{H}^{\pm\pm}$ boson signal is set to zero for the derivation of the exclusion limits on the $\sigma_\mathrm{VBF}(\mathrm{H}^{\pm}) \mathcal{B}(\mathrm{H}^{\pm} \rightarrow \mathrm{W}^{\pm}\mathrm{Z})$. Values above the curves are excluded.

Expected and observed exclusion limits at 95\% confidence level for $s_{\mathrm{H}}$ as functions of $m_{\mathrm{H}_{5}}$ in the Georgi-Machacek model. Values above the curves are excluded. The exclusion limits for $s_{\mathrm{H}}$ are shown up to $m_{\mathrm{H}_{5}}=2000$ GeV, given the low sensitivity in the Georgi-Machacek model for values above that mass.

Covariance matrix for all the bins used in Figure 4. The covariance matrix from the background-only fit is shown.

Measured fiducial cross sections for the $\mathrm{W}^\pm_{\mathrm{L}}\mathrm{W}^\pm_{\mathrm{L}}$ and $\mathrm{W}^\pm_{\mathrm{X}}\mathrm{W}^\pm_{\mathrm{T}}$ processes, and for the $\mathrm{W}^\pm_{\mathrm{L}}\mathrm{W}^\pm_{\mathrm{X}}$ and $\mathrm{W}^\pm_{\mathrm{T}}\mathrm{W}^\pm_{\mathrm{T}}$ processes for the helicity eigenstates defined in the parton-parton center-of-mass frame. The combination of the statistical and systematic uncertainties is shown. $\mathcal{B}$ is the branching fraction for $\mathrm{W}\mathrm{W} \rightarrow \ell \nu \ell' \nu$. The fiducial region is defined by requiring two same-sign leptons with $p_{T}>20$, $|\eta|<2.5$, and $m_{ll}>20$, and two jets with $m_{jj}>500$ and $|\Delta \eta_{jj}|>2.5$. The jets at generator level are clustered from stable particles, excluding neutrinos, using the anti-kt clustering algorithm with R = 0.4, and are required to have $p_{T}>50$ and $|\eta|<4.7$. The jets within $\Delta R<0.4$ of the selected charged leptons are not included. The theoretical predictions including the $\mathcal{O}(\alpha_{s}\alpha^6)$ and $\mathcal{O}(\alpha^7)$ corrections to the \MGvATNLO LO cross sections, as described in arXiv:2009.09429, are also shown. The theoretical uncertainties include statistical, PDF, and LO scale uncertainties.

The measured inclusive fiducial cross section measurements. The WW fiducial region is defined by requiring two same-sign leptons with $p_{T}>20$, $|\eta|<2.5$, and $m_{ll}>20$, and two jets with $m_{jj}>500$ and $|\Delta \eta_{jj}|>2.5$. The jets at generator level are clustered from stable particles, excluding neutrinos, using the anti-kt clustering algorithm with R = 0.4, and are required to have $p_{T}>50$ and $|\eta|<4.7$. The jets within $\Delta R<0.4$ of the selected charged leptons are not included. The WZ fiducial region is defined by requiring three leptons with $p_{T}>20$, $|\eta|<2.5$, a pair of opposite charge same-flavor lepton pair with $|m_{ll}-m_{Z}|<15$, and two jets with $m_{jj}>500$ and $|\Delta \eta_{jj}|>2.5$.

Distributions of $m_{jj}$ in the WW signal region.

Distributions of $m_{ll}$ in the WW signal region.

Distributions of $m_{jj}$ in the WZ signal region.

Distributions of BDT score in the WZ signal region.

Absolute WW cross section in $m_{jj}$ bins.

Normalized WW cross section in $m_{jj}$ bins.

Absolute WW cross section in $m_{ll}$ bins.

Normalized WW cross section in $m_{ll}$ bins.

Absolute WW cross section in $p_{T}^{l max}$ bins.

Normalized WW cross section in $p_{T}^{l max}$ bins.

Absolute WZ cross section in $m_{jj}$ bins.

Normalized WZ cross section in $m_{jj}$ bins.

Distributions of $m_{T}^{WW}$ in the WW signal region.

Distributions of $m_{T}^{WZ}$ in the WZ signal region.

Observed and expected lower and upper 95\% confidence level limits in TeV$^{-4}$ on the parameters of the quartic, obtained without using any unitarization procedure.

Observed and expected lower and upper 95\% confidence level limits in TeV$^{-4}$ on the parameters of the quartic by cutting the EFT expansion at the unitarity limit.

The relative statistical and systematic uncertainties in % from various sources for the absolute cross section measurements in bins of Z pt on dimuons.

The relative statistical and systematic uncertainties in % from various sources for the absolute cross section measurements in bins of Z pt on dielectrons.

The relative statistical and systematic uncertainties in % from various sources for the absolute cross section measurements in bins of |y(Z)| on dimuons.

The relative statistical and systematic uncertainties in % from various sources for the absolute cross section measurements in bins of |y(Z)| on dielectrons.

The relative statistical and systematic uncertainties in % from various sources for the absolute cross section measurements in bins of $\phi^{\scriptscriptstyle *}_\eta$ on dimuons.

The relative statistical and systematic uncertainties in % from various sources for the absolute cross section measurements in bins of $\phi^{\scriptscriptstyle *}_\eta$ on dielectrons.

The relative statistical and systematic uncertainties in % from various sources for the normalized cross section measurements in bins of Z pt on dimuons.

The relative statistical and systematic uncertainties in % from various sources for the normalized cross section measurements in bins of Z pt on dielectrons.

The relative statistical and systematic uncertainties in % from various sources for the normalized cross section measurements in bins of |y(Z)| on dimuons.

The relative statistical and systematic uncertainties in % from various sources for the normalized cross section measurements in bins of |y(Z)| on dielectrons.

The relative statistical and systematic uncertainties in % from various sources for the normalized cross section measurements in bins of $\phi^{\scriptscriptstyle *}_\eta$ on dimuons.

The relative statistical and systematic uncertainties in % from various sources for the normalized cross section measurements in bins of $\phi^{\scriptscriptstyle *}_\eta$ on dielectrons.

The relative statistical and systematic uncertainties in % from various sources for the absolute double-differential cross section measurements in bins of Z pt for the 0 < |y(Z)| < 0.4 bin in the dimuon final state.

The relative statistical and systematic uncertainties in % from various sources for the absolute double-differential cross section measurements in bins of Z pt for the 0.4 < |y(Z)| < 0.8 bin in the dimuon final state.

The relative statistical and systematic uncertainties in % from various sources for the absolute double-differential cross section measurements in bins of Z pt for the 0.8 < |y(Z)| < 1.2 bin in the dimuon final state.

The relative statistical and systematic uncertainties in % from various sources for the absolute double-differential cross section measurements in bins of Z pt for the 1.2 < |y(Z)| < 1.6 bin in the dimuon final state.

The relative statistical and systematic uncertainties in % from various sources for the absolute double-differential cross section measurements in bins of Z pt for the 1.6 < |y(Z)| < 2.4 bin in the dimuon final state.

The relative statistical and systematic uncertainties in % from various sources for the absolute double-differential cross section measurements in bins of Z pt for the 0 < |y(Z)| < 0.4 bin in the dielectron final state.

The relative statistical and systematic uncertainties in % from various sources for the absolute double-differential cross section measurements in bins of Z pt for the 0.4 < |y(Z)| < 0.8 bin in the dielectron final state.

The relative statistical and systematic uncertainties in % from various sources for the absolute double-differential cross section measurements in bins of Z pt for the 0.8 < |y(Z)| < 1.2 bin in the dielectron final state.

The relative statistical and systematic uncertainties in % from various sources for the absolute double-differential cross section measurements in bins of Z pt for the 1.2 < |y(Z)| < 1.6 bin in the dielectron final state.

The relative statistical and systematic uncertainties in % from various sources for the absolute double-differential cross section measurements in bins of Z pt for the 1.6 < |y(Z)| < 2.4 bin in the dielectron final state.

The measured absolute cross sections in bins of |y(Z)|, using dressed level leptons. The cross sections are normalized by the bin width.

The measured absolute cross sections in bins of Z pt, using dressed level leptons. The cross sections are normalized by the bin width.

The measured absolute cross sections in bins of $\phi^{\scriptscriptstyle *}_\eta$, using dressed level leptons. The cross sections are normalized by the bin width.

The measured absolute cross sections in bins of Z pt different |y(Z)| bins, using dressed level leptons. The cross sections are normalized by the bin width.

The measured normalized cross sections in bins of Z pt, using dressed level leptons. The cross sections are divided by the bin width and then the distributions are normalized to 1.

The measured normalized cross sections in bins of $\phi^{\scriptscriptstyle *}_\eta$, using dressed level leptons. The cross sections are divided by the bin width and then the distributions are normalized to 1.

The measured normalized cross sections in bins of |y(Z)|, using dressed level leptons. The cross sections are divided by the bin width and then the distributions are normalized to 1.

The measured normalized cross sections in bins of Z pt in |y(Z)| bins, using dressed level leptons. The cross sections are divided by the bin width and then the distributions are normalized to 1.

Covariance matrix using dressed level leptons for all bins used in bins of Z pt in the dimuon final state.

Covariance matrix using dressed level leptons for all bins used in bins of |y(Z)| in the dimuon final state.

Covariance matrix using dressed level leptons for all bins used in bins of $\phi^{\scriptscriptstyle *}_\eta$ in the dimuon final state.

Covariance matrix using dressed level leptons for all bins used in bins of Z pt for the 0 < |y(Z)| < 0.4 bin in the dimuon final state.

Covariance matrix using dressed level leptons for all bins used in bins of Z pt for the 0.4 < |y(Z)| < 0.8 bin in the dimuon final state.

Covariance matrix using dressed level leptons for all bins used in bins of Z pt for the 0.8 < |y(Z)| < 1.2 bin in the dimuon final state.

Covariance matrix using dressed level leptons for all bins used in bins of Z pt for the 1.2 < |y(Z)| < 1.6 bin in the dimuon final state.

Covariance matrix using dressed level leptons for all bins used in bins of Z pt for the 1.6 < |y(Z)| < 2.4 bin in the dimuon final state.

Covariance matrix using dressed level leptons for all bins used in bins of Z pt in the dielectron final state.

Covariance matrix using dressed level leptons for all bins used in bins of |y(Z)| in the dielectron final state.

Covariance matrix using dressed level leptons for all bins used in bins of $\phi^{\scriptscriptstyle *}_\eta$ in the dielectron final state.

Covariance matrix using dressed level leptons for all bins used in bins of Z pt for the 0 < |y(Z)| < 0.4 bin in the dielectron final state.

Covariance matrix using dressed level leptons for all bins used in bins of Z pt for the 0.4 < |y(Z)| < 0.8 bin in the dielectron final state.

Covariance matrix using dressed level leptons for all bins used in bins of Z pt for the 0.8 < |y(Z)| < 1.2 bin in the dielectron final state.

Covariance matrix using dressed level leptons for all bins used in bins of Z pt for the 1.2 < |y(Z)| < 1.6 bin in the dielectron final state.

Covariance matrix using dressed level leptons for all bins used in bins of Z pt for the 1.6 < |y(Z)| < 2.4 bin in the dielectron final state.

Covariance matrix using born level leptons for all bins used in bins of Z pt in the dimuon final state.

Covariance matrix using born level leptons for all bins used in bins of |y(Z)| in the dimuon final state.

Covariance matrix using born level leptons for all bins used in bins of $\phi^{\scriptscriptstyle *}_\eta$ in the dimuon final state.

Covariance matrix using born level leptons for all bins used in bins of Z pt for the 0 < |y(Z)| < 0.4 bin in the dimuon final state.

Covariance matrix using born level leptons for all bins used in bins of Z pt for the 0.4 < |y(Z)| < 0.8 bin in the dimuon final state.

Covariance matrix using born level leptons for all bins used in bins of Z pt for the 0.8 < |y(Z)| < 1.2 bin in the dimuon final state.

Covariance matrix using born level leptons for all bins used in bins of Z pt for the 1.2 < |y(Z)| < 1.6 bin in the dimuon final state.

Covariance matrix using born level leptons for all bins used in bins of Z pt for the 1.6 < |y(Z)| < 2.4 bin in the dimuon final state.

Covariance matrix using born level leptons for all bins used in bins of Z pt in the dielectron final state.

Covariance matrix using born level leptons for all bins used in bins of |y(Z)| in the dielectron final state.

Covariance matrix using born level leptons for all bins used in bins of $\phi^{\scriptscriptstyle *}_\eta$ in the dielectron final state.

Covariance matrix using born level leptons for all bins used in bins of Z pt for the 0 < |y(Z)| < 0.4 bin in the dielectron final state.

Covariance matrix using born level leptons for all bins used in bins of Z pt for the 0.4 < |y(Z)| < 0.8 bin in the dielectron final state.

Covariance matrix using born level leptons for all bins used in bins of Z pt for the 0.8 < |y(Z)| < 1.2 bin in the dielectron final state.

Covariance matrix using born level leptons for all bins used in bins of Z pt for the 1.2 < |y(Z)| < 1.6 bin in the dielectron final state.

Covariance matrix using born level leptons for all bins used in bins of Z pt for the 1.6 < |y(Z)| < 2.4 bin in the dielectron final state.

Covariance matrix for normalized cross sections using dressed level leptons for all bins used in bins of Z pt in the dimuon final state.

Covariance matrix for normalized cross sections using dressed level leptons for all bins used in bins of |y(Z)| in the dimuon final state.

Covariance matrix for normalized cross sections using dressed level leptons for all bins used in bins of $\phi^{\scriptscriptstyle *}_\eta$ in the dimuon final state.

Covariance matrix for normalized cross sections using dressed level leptons for all bins used in bins of Z pt for the 0 < |y(Z)| < 0.4 bin in the dimuon final state.

Covariance matrix for normalized cross sections using dressed level leptons for all bins used in bins of Z pt for the 0.4 < |y(Z)| < 0.8 bin in the dimuon final state.

Covariance matrix for normalized cross sections using dressed level leptons for all bins used in bins of Z pt for the 0.8 < |y(Z)| < 1.2 bin in the dimuon final state.

Covariance matrix for normalized cross sections using dressed level leptons for all bins used in bins of Z pt for the 1.2 < |y(Z)| < 1.6 bin in the dimuon final state.

Covariance matrix for normalized cross sections using dressed level leptons for all bins used in bins of Z pt for the 1.6 < |y(Z)| < 2.4 bin in the dimuon final state.

Covariance matrix for normalized cross sections using dressed level leptons for all bins used in bins of Z pt in the dielectron final state.

Covariance matrix for normalized cross sections using dressed level leptons for all bins used in bins of |y(Z)| in the dielectron final state.

Covariance matrix for normalized cross sections using dressed level leptons for all bins used in bins of $\phi^{\scriptscriptstyle *}_\eta$ in the dielectron final state.

Covariance matrix for normalized cross sections using dressed level leptons for all bins used in bins of Z pt for the 0 < |y(Z)| < 0.4 bin in the dielectron final state.

Covariance matrix for normalized cross sections using dressed level leptons for all bins used in bins of Z pt for the 0.4 < |y(Z)| < 0.8 bin in the dielectron final state.

Covariance matrix for normalized cross sections using dressed level leptons for all bins used in bins of Z pt for the 0.8 < |y(Z)| < 1.2 bin in the dielectron final state.

Covariance matrix for normalized cross sections using dressed level leptons for all bins used in bins of Z pt for the 1.2 < |y(Z)| < 1.6 bin in the dielectron final state.

Covariance matrix for normalized cross sections using dressed level leptons for all bins used in bins of Z pt for the 1.6 < |y(Z)| < 2.4 bin in the dielectron final state.

Covariance matrix for normalized cross sections using born level leptons for all bins used in bins of Z pt in the dimuon final state.

Covariance matrix for normalized cross sections using born level leptons for all bins used in bins of |y(Z)| in the dimuon final state.

Covariance matrix for normalized cross sections using born level leptons for all bins used in bins of $\phi^{\scriptscriptstyle *}_\eta$ in the dimuon final state.

Covariance matrix for normalized cross sections using born level leptons for all bins used in bins of Z pt for the 0 < |y(Z)| < 0.4 bin in the dimuon final state.

Covariance matrix for normalized cross sections using born level leptons for all bins used in bins of Z pt for the 0.4 < |y(Z)| < 0.8 bin in the dimuon final state.

Covariance matrix for normalized cross sections using born level leptons for all bins used in bins of Z pt for the 0.8 < |y(Z)| < 1.2 bin in the dimuon final state.

Covariance matrix for normalized cross sections using born level leptons for all bins used in bins of Z pt for the 1.2 < |y(Z)| < 1.6 bin in the dimuon final state.

Covariance matrix for normalized cross sections using born level leptons for all bins used in bins of Z pt for the 1.6 < |y(Z)| < 2.4 bin in the dimuon final state.

Covariance matrix for normalized cross sections using born level leptons for all bins used in bins of Z pt in the dielectron final state.

Covariance matrix for normalized cross sections using born level leptons for all bins used in bins of |y(Z)| in the dielectron final state.

Covariance matrix for normalized cross sections using born level leptons for all bins used in bins of $\phi^{\scriptscriptstyle *}_\eta$ in the dielectron final state.

Covariance matrix for normalized cross sections using born level leptons for all bins used in bins of Z pt for the 0 < |y(Z)| < 0.4 bin in the dielectron final state.

Covariance matrix for normalized cross sections using born level leptons for all bins used in bins of Z pt for the 0.4 < |y(Z)| < 0.8 bin in the dielectron final state.

Covariance matrix for normalized cross sections using born level leptons for all bins used in bins of Z pt for the 0.8 < |y(Z)| < 1.2 bin in the dielectron final state.

Covariance matrix for normalized cross sections using born level leptons for all bins used in bins of Z pt for the 1.2 < |y(Z)| < 1.6 bin in the dielectron final state.

Covariance matrix for normalized cross sections using born level leptons for all bins used in bins of Z pt for the 1.6 < |y(Z)| < 2.4 bin in the dielectron final state.

The measured absolute cross sections in bins of |y(Z)|, using born level leptons. The cross sections are normalized by the bin width.

The measured absolute cross sections in bins of Z pt, using born level leptons. The cross sections are normalized by the bin width. The first bin (0-1 GeV) is not shown as large differences were observed in aMC@NLO and POWHEG predictions at Born level.

The measured absolute cross sections in bins of $\phi^{\scriptscriptstyle *}_\eta$, using born level leptons. The cross sections are normalized by the bin width.

The measured absolute cross sections in bins of Z pt different |y(Z)| bins, using born level leptons. The cross sections are normalized by the bin width. The first bin (0-1 GeV) is not shown as large differences were observed in aMC@NLO and POWHEG predictions at Born level.

The measured normalized cross sections in bins of Z pt, using born level leptons. The first bin (0-1 GeV) is not shown as large differences were observed in aMC@NLO and POWHEG predictions at Born level. The cross sections are divided by the bin width and then the distributions are normalized to 1.

The measured normalized cross sections in bins of $\phi^{\scriptscriptstyle *}_\eta$, using born level leptons. The cross sections are divided by the bin width and then the distributions are normalized to 1.

The measured normalized cross sections in bins of |y(Z)|, using born level leptons.

The measured normalized cross sections (left) in bins of Z pt in |y(Z)| bins, using born level leptons. The first bin (0-1 GeV) is not shown as large differences were observed in aMC@NLO and POWHEG predictions at Born level. The cross sections are divided by the bin width and then the distributions are normalized to 1.

Covariance matrix using dressed level leptons for all bins used in bins of Z pt in the dilepton final state.

Covariance matrix using dressed level leptons for all bins used in bins of |y(Z)| in the dilepton final state.

Covariance matrix using dressed level leptons for all bins used in bins of $\phi^{\scriptscriptstyle *}_\eta$ in the dilepton final state.

Covariance matrix using dressed level leptons for all bins used in bins of Z pt for the 0 < |y(Z)| < 0.4 bin in the dilepton final state.

Covariance matrix using dressed level leptons for all bins used in bins of Z pt for the 0.4 < |y(Z)| < 0.8 bin in the dilepton final state.

Covariance matrix using dressed level leptons for all bins used in bins of Z pt for the 0.8 < |y(Z)| < 1.2 bin in the dilepton final state.

Covariance matrix using dressed level leptons for all bins used in bins of Z pt for the 1.2 < |y(Z)| < 1.6 bin in the dilepton final state.

Covariance matrix using dressed level leptons for all bins used in bins of Z pt for the 1.6 < |y(Z)| < 2.4 bin in the dilepton final state.

Covariance matrix using born level leptons for all bins used in bins of Z pt in the dilepton final state.

Covariance matrix using born level leptons for all bins used in bins of |y(Z)| in the dilepton final state.

Covariance matrix using born level leptons for all bins used in bins of $\phi^{\scriptscriptstyle *}_\eta$ in the dilepton final state.

Covariance matrix using born level leptons for all bins used in bins of Z pt for the 0 < |y(Z)| < 0.4 bin in the dilepton final state.

Covariance matrix using born level leptons for all bins used in bins of Z pt for the 0.4 < |y(Z)| < 0.8 bin in the dilepton final state.

Covariance matrix using born level leptons for all bins used in bins of Z pt for the 0.8 < |y(Z)| < 1.2 bin in the dilepton final state.

Covariance matrix using born level leptons for all bins used in bins of Z pt for the 1.2 < |y(Z)| < 1.6 bin in the dilepton final state.

Covariance matrix using born level leptons for all bins used in bins of Z pt for the 1.6 < |y(Z)| < 2.4 bin in the dilepton final state.

Covariance matrix for normalized cross sections using dressed level leptons for all bins used in bins of Z pt in the dilepton final state.

Covariance matrix for normalized cross sections using dressed level leptons for all bins used in bins of |y(Z)| in the dilepton final state.

Covariance matrix for normalized cross sections using dressed level leptons for all bins used in bins of $\phi^{\scriptscriptstyle *}_\eta$ in the dilepton final state.

Covariance matrix for normalized cross sections using dressed level leptons for all bins used in bins of Z pt for the 0 < |y(Z)| < 0.4 bin in the dilepton final state.

Covariance matrix for normalized cross sections using dressed level leptons for all bins used in bins of Z pt for the 0.4 < |y(Z)| < 0.8 bin in the dilepton final state.

Covariance matrix for normalized cross sections using dressed level leptons for all bins used in bins of Z pt for the 0.8 < |y(Z)| < 1.2 bin in the dilepton final state.

Covariance matrix for normalized cross sections using dressed level leptons for all bins used in bins of Z pt for the 1.2 < |y(Z)| < 1.6 bin in the dilepton final state.

Covariance matrix for normalized cross sections using dressed level leptons for all bins used in bins of Z pt for the 1.6 < |y(Z)| < 2.4 bin in the dilepton final state.

Covariance matrix for normalized cross sections using born level leptons for all bins used in bins of Z pt in the dilepton final state.

Covariance matrix for normalized cross sections using born level leptons for all bins used in bins of |y(Z)| in the dilepton final state.

Covariance matrix for normalized cross sections using born level leptons for all bins used in bins of $\phi^{\scriptscriptstyle *}_\eta$ in the dilepton final state.

Covariance matrix for normalized cross sections using born level leptons for all bins used in bins of Z pt for the 0 < |y(Z)| < 0.4 bin in the dilepton final state.

Covariance matrix for normalized cross sections using born level leptons for all bins used in bins of Z pt for the 0.4 < |y(Z)| < 0.8 bin in the dilepton final state.

Covariance matrix for normalized cross sections using born level leptons for all bins used in bins of Z pt for the 0.8 < |y(Z)| < 1.2 bin in the dilepton final state.

Covariance matrix for normalized cross sections using born level leptons for all bins used in bins of Z pt for the 1.2 < |y(Z)| < 1.6 bin in the dilepton final state.

Covariance matrix for normalized cross sections using born level leptons for all bins used in bins of Z pt for the 1.6 < |y(Z)| < 2.4 bin in the dilepton final state.

$\mathrm{M_{\mathrm{T}}}$ distribution for the $\mathrm{e}\mu$ control region.

$\mathrm{M_{\mathrm{T}}}$ distribution for the WZ control region.

$\mathrm{M_{\mathrm{T}}}$ distribution for the ZZ control region.

$\mathrm{M_{\mathrm{T}}}$ distribution for the signal region, events with $|\eta_{\gamma}|<1$.

$\mathrm{M_{\mathrm{T}}}$ distribution for the signal region, events with $|\eta_{\gamma}|>1$.

Muon reconstruction, identification, and isolation efficiency as a function of $|\eta|$ and $\mathrm{p_\mathrm{T}}$.

Electron reconstruction, identification, and isolation efficiency as a function of $|\eta|$ and $\mathrm{p_\mathrm{T}}$.

Photon efficiency as a function of $|\eta|$ and $\mathrm{p_\mathrm{T}}$.

Covariance matrix for all bins used in the analysis. There are 45 bins in total, 15 for every data-taking year. For every year, the first bin corresponds to events in the $\mathrm{e}\mu$ control region, the following five bins correspond to events with $|\eta^\gamma|< 1$ in the signal region, the next five bins correspond to events with $|\eta^\gamma|> 1$ in the signal region, the next two bins correspond to events in the WZ control region, and finally the last two bins correspond to events in the ZZ control region.