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.

Distribution of the three leading leptons flavour in the SR-WZ with uncertainties evaluated after the inclusive cross section fit

Efficiency, acceptance, and proportion of events with leptonic tau decays in WZ production

WZ fiducial cross section in the four flavour exclusive and the flavour inclusive channels

WZ total cross section extrapolated from the four flavour exclusive and the flavour inclusive channels

Distribution of the total lepton charge in the SR-WZ with uncertainties evaluated after the inclusive cross section fit

W$^{+}$Z fiducial cross section in the four flavour exclusive and the flavour inclusive channels

W$^{-}$Z fiducial cross section in the four flavour exclusive and the flavour inclusive channels

WZ charge asymmetry ratio measured on each of the four flavour exclusive and the flavour inclusive channels

Distribution of the cosine of the W polarization angle times total lepton charge in the SR-WZ with uncertainties evaluated after the W polarization fit

Distribution of the cosine of the Z polarization angle in the SR-WZ with uncertainties evaluated after the Z polarization fit

Best fits to the W and Z polarization fractions

2D confidence regions at the 68, 95, and 99% CL in the $f_O^W$-$f_{L}^W-f_R^W$ plane

2D confidence regions at the 68, 95, and 99% CL in the $f_O^Z$-$f_{L}^Z-f_R^Z$ plane

Distribution of the invariant mass of the WZ system in the SR-WZ with uncertainties evaluated after the inclusive cross section fit

Best fit values and one dimensional confidence regions in several EFT coefficients obtained from the EFT fit considering both the SM interferences and purely BSM (order $\Lambda^{-2}$ and $\Lambda^{-4}$) terms

2D confidence regions at the 68, 95, and 99% CL in the $c_{www}$-$c_{w}$ plane

2D confidence regions at the 68, 95, and 99% CL in the $c_{w}$-$c_{b}$ plane

2D confidence regions at the 68, 95, and 99% CL in the $c_{www}$-$c_{w}$ plane

Best fit values and one dimensional confidence regions in several EFT coefficients obtained from the EFT fit considering only the SM-EFT interference (order $\Lambda^{-2}$) terms

Evolution of the best fit and expected and observed 95% CI for the $c_{w}$ parameter as a function of the cutoff scale

Evolution of the best fit and expected and observed 95% CI for the $c_{b}$ parameter as a function of the cutoff scale

Evolution of the best fit and expected and observed 95% CI for the $c_{www}$ parameter as a function of the cutoff scale

Evolution of the best fit and expected and observed 95% CI for the $\tilde{c}_{www}$ parameter as a function of the cutoff scale

Evolution of the best fit and expected and observed 95% CI for the $\tilde{c}_{w}$ parameter as a function of the cutoff scale

Differential cross section with respect to the transverse momentum of the Z boson

Correlation matrix for the unfolded results obtained using NNLO bias, area con-straint, and no additional regularization for the $p_{T}$ of the Z boson

Response matrix for the $p_{T}$ of the Z boson obtained with POWHEG

Differential cross section with respect to the transverse momentum of the leading jet

Correlation matrix for the unfolded results obtained using NNLO bias, area con-straint, and no additional regularization for the $p_{T}$ of the leading jet

Response matrix for the $p_{T}$ of the leading jet obtained with POWHEG

Differential cross section with respect to the jet multiplicity

Correlation matrix for the unfolded results obtained using NNLO bias, area con-straint, and no additional regularization for the jet multiplicity

Response matrix for the jet multiplicity obtained with POWHEG

Differential cross section with respect to the invariant mass of the WZ system

Correlation matrix for the unfolded results obtained using NNLO bias, area con-straint, and no additional regularization for the invariant mass of the WZ system

Response matrix for the invariant mass of the WZ system obtained with POWHEG

Differential cross section with respect to the transverse momentum of the lepton associated to the W boson

Correlation matrix for the unfolded results obtained using NNLO bias, area con-straint, and no additional regularization for the $p_{T}$ of the lepton associated to the W boson

Response matrix for the $p_{T}$ of the lepton associated to the W boson obtained with POWHEG

Differential cross section with respect to the transverse momentum of the lepton associated to the W boson, W$^{+}$Z only

Correlation matrix for the unfolded results obtained using NNLO bias, area con-straint, and no additional regularization for the $p_{T}$ of the lepton associated to the W boson, W$^{+}$Z only

Differential cross section with respect to the transverse momentum of the lepton associated to the W boson, W$^{-}$Z only

Correlation matrix for the unfolded results obtained using NNLO bias, area con-straint, and no additional regularization for the $p_{T}$ of the lepton associated to the W boson, W$^{-}$Z only

Differential cross section with respect to the cosine of the W polarization angle times total lepton charge

Correlation matrix for the unfolded results obtained using NNLO bias, area con-straint, and no additional regularization for the cosine of the W polarization angle times total lepton charge

Response matrix for the cosine of the W polarization angle times total lepton charge obtained with POWHEG

Differential cross section with respect to the cosine of the W polarization angle times total lepton charge, W$^{+}$Z only

Correlation matrix for the unfolded results obtained using NNLO bias, area con-straint, and no additional regularization for the cosine of the W polarization angle times total lepton charge, W$^{+}$Z only

Differential cross section with respect to the cosine of the W polarization angle times total lepton charge, W$^{-}$Z only

Correlation matrix for the unfolded results obtained using NNLO bias, area con-straint, and no additional regularization for the cosine of the W polarization angle times total lepton charge, W$^{-}$Z only

Differential cross section with respect to the cosine of the Z polarization angle

Correlation matrix for the unfolded results obtained using NNLO bias, area con-straint, and no additional regularization for the cosine of the Z polarization angle

Response matrix for the cosine of the Z polarization angle obtained with POWHEG

Differential cross section with respect to the cosine of the Z polarization angle, W$^{+}$Z only

Correlation matrix for the unfolded results obtained using NNLO bias, area con-straint, and no additional regularization for the cosine of the Z polarization angle, W$^{+}$Z only

Differential cross section with respect to the cosine of the Z polarization angle, W$^{-}$Z only

Correlation matrix for the unfolded results obtained using NNLO bias, area con-straint, and no additional regularization for the cosine of the Z polarization angle, W$^{-}$Z only

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

Theoretical cross-sections for heavy composite majorana neutrino signals for $\Lambda$=13 TeV.

Exclusion limits on the product of cross section and branching fraction for the electron channel, as a function of the resonance mass hypothesis.

Exclusion limits on the product of cross section and branching fraction for the muon channel, as a function of the resonance mass hypothesis.

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.

Observed (black solid) and expected (black dotted) limits at 95% confidence level on the product of the the cross section for bbA prodcution and the branching fraction $\mathrm{A} \rightarrow \tau\tau$, obtained from the combination of the $\mathrm{e}\tau_\mathrm{h}$ and $\mu\tau_\mathrm{h}$ channels in the 1 b tag category, as a function of the A boson mass. The green and yellow bands represent the one and two standard deviation uncertainties in the expected limits.

Observed (black solid) and expected (black dotted) limits at 95% confidence level on the product of the the cross section for bbA prodcution and the branching fraction $\mathrm{A} \rightarrow \tau\tau$, obtained in the $\mathrm{e}\tau_\mathrm{h}$ channel of the 1 b tag category, as a function of the A boson mass. The green and yellow bands represent the one and two standard deviation uncertainties in the expected limits.

Observed (black solid) and expected (black dotted) limits at 95% confidence level on the product of the the cross section for bbA prodcution and the branching fraction $\mathrm{A} \rightarrow \tau\tau$, obtained in the $\mu\tau_\mathrm{h}$ channel of the 1 b tag category, as a function of the A boson mass. The green and yellow bands represent the one and two standard deviation uncertainties in the expected limits.

The product of acceptance, efficiency, and branching fraction of the qbX signal with $\mathrm{X} \rightarrow \tau\tau$ and $m_{\mathrm{B}}=170~\mathrm{GeV}$ in $\mathrm{e}\tau_\mathrm{h}$ and $\mu\tau_\mathrm{h}$ channels of the 1b1f and 1b1c event categories, as a function of the X boson mass. The selections are as described in the paper. The uncertainty refers to the statistical uncertainty only.

Observed $m_{\tau\tau}$ distribution in the $\mathrm{e}\tau_\mathrm{h}$ channel of the 1b1f category, compared to the expected SM background contributions. The distributions of the qbX signal with $m_{\mathrm{B}}=170~\mathrm{GeV}$, and X boson mass 40 and 60 GeV are overlaid to illustrate the sensitivity. They are normalized to the cross section times branching fraction of 20 pb. The uncertainty band represents the sum in quadrature of statistical and systematic uncertainties obtained from the fit. The lower panel shows the ratio between the observed and expected events in each bin.

Observed $m_{\tau\tau}$ distribution in the $\mathrm{e}\tau_\mathrm{h}$ channel of the 1b1c category, compared to the expected SM background contributions. The distributions of the qbX signal with $m_{\mathrm{B}}=170~\mathrm{GeV}$, and X boson mass 40 and 60 GeV are overlaid to illustrate the sensitivity. They are normalized to the cross section times branching fraction of 20 pb. The uncertainty band represents the sum in quadrature of statistical and systematic uncertainties obtained from the fit. The lower panel shows the ratio between the observed and expected events in each bin.

Observed $m_{\tau\tau}$ distribution in the $\mu\tau_\mathrm{h}$ channel of the 1b1f category, compared to the expected SM background contributions. The distributions of the qbX signal with $m_{\mathrm{B}}=170~\mathrm{GeV}$, and X boson mass 40 and 60 GeV are overlaid to illustrate the sensitivity. They are normalized to the cross section times branching fraction of 20 pb. The uncertainty band represents the sum in quadrature of statistical and systematic uncertainties obtained from the fit. The lower panel shows the ratio between the observed and expected events in each bin.

Observed $m_{\tau\tau}$ distribution in the $\mu\tau_\mathrm{h}$ channel of the 1b1c category, compared to the expected SM background contributions. The distributions of the qbX signal with $m_{\mathrm{B}}=170~\mathrm{GeV}$, and X boson mass 40 and 60 GeV are overlaid to illustrate the sensitivity. They are normalized to the cross section times branching fraction of 20 pb. The uncertainty band represents the sum in quadrature of statistical and systematic uncertainties obtained from the fit. The lower panel shows the ratio between the observed and expected events in each bin.

Observed (black solid) and expected (black dotted) limits at 95% confidence level on the product of the the cross section for the prodcution of the qbX signal with mB = 170 GeV and the branching fraction $\mathrm{X} \rightarrow \tau\tau$, obtained from the combination of the $\mathrm{e}\tau_\mathrm{h}$ and $\mu\tau_\mathrm{h}$ channels, and the 1b1f and 1b1c categories, as a function of the X boson mass. The green and yellow bands represent the one and two standard deviation uncertainties in the expected limits.

Observed (black solid) and expected (black dotted) limits at 95% confidence level on the product of the the cross section for the prodcution of the qbX signal with mB = 300 GeV and the branching fraction $\mathrm{X} \rightarrow \tau\tau$, obtained from the combination of the $\mathrm{e}\tau_\mathrm{h}$ and $\mu\tau_\mathrm{h}$ channels, and the 1b1f and 1b1c categories, as a function of the X boson mass. The green and yellow bands represent the one and two standard deviation uncertainties in the expected limits.

Observed (black solid) and expected (black dotted) limits at 95% confidence level on the product of the the cross section for the prodcution of the qbX signal with mB = 450 GeV and the branching fraction $\mathrm{X} \rightarrow \tau\tau$, obtained from the combination of the $\mathrm{e}\tau_\mathrm{h}$ and $\mu\tau_\mathrm{h}$ channels, and the 1b1f and 1b1c categories, as a function of the X boson mass. The green and yellow bands represent the one and two standard deviation uncertainties in the expected limits.

Observed (black solid) and expected (black dotted) limits at 95% confidence level on the product of the the cross section for the prodcution of the qbX signal with mB = 170 GeV and the branching fraction $\mathrm{X} \rightarrow \tau\tau$, obtained in the $\mathrm{e}\tau_\mathrm{h}$ channel of the 1b1f category, as a function of the X boson mass. The green and yellow bands represent the one and two standard deviation uncertainties in the expected limits.

Observed (black solid) and expected (black dotted) limits at 95% confidence level on the product of the the cross section for the prodcution of the qbX signal with mB = 170 GeV and the branching fraction $\mathrm{X} \rightarrow \tau\tau$, obtained in the $\mathrm{e}\tau_\mathrm{h}$ channel of the 1b1c category, as a function of the X boson mass. The green and yellow bands represent the one and two standard deviation uncertainties in the expected limits.

Observed (black solid) and expected (black dotted) limits at 95% confidence level on the product of the the cross section for the prodcution of the qbX signal with mB = 170 GeV and the branching fraction $\mathrm{X} \rightarrow \tau\tau$, obtained in the $\mu\tau_\mathrm{h}$ channel of the 1b1f category, as a function of the X boson mass. The green and yellow bands represent the one and two standard deviation uncertainties in the expected limits.

Observed (black solid) and expected (black dotted) limits at 95% confidence level on the product of the the cross section for the prodcution of the qbX signal with mB = 170 GeV and the branching fraction $\mathrm{X} \rightarrow \tau\tau$, obtained in the $\mu\tau_\mathrm{h}$ channel of the 1b1c category, as a function of the X boson mass. The green and yellow bands represent the one and two standard deviation uncertainties in the expected limits.

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

Upper limits at 95% CL on $\sigma\times\mathcal{B}^2(\mathrm{H}_1\rightarrow b\bar{b})$ as a function of $m_{\mathrm{H_1}}$ for $m_{\mathrm{SUSY}}=1600\;\mathrm{GeV}$.

Upper limits at 95% CL on $\sigma\times\mathcal{B}^2(\mathrm{H}_1\rightarrow b\bar{b})$ as a function of $m_{\mathrm{H_1}}$ for $m_{\mathrm{SUSY}}=2000\;\mathrm{GeV}$.

Upper limits at 95% CL on $\sigma\times\mathcal{B}^2(\mathrm{H}_1\rightarrow b\bar{b})$ as a function of $m_{\mathrm{H_1}}$ for $m_{\mathrm{SUSY}}=2200\;\mathrm{GeV}$.

Upper limits at 95% CL on $\sigma\times\mathcal{B}^2(\mathrm{H}_1\rightarrow b\bar{b})$ as a function of $m_{\mathrm{H_1}}$ for $m_{\mathrm{SUSY}}=2400\;\mathrm{GeV}$.

Upper limits at 95% CL on $\sigma\times\mathcal{B}^2(\mathrm{H}_1\rightarrow b\bar{b})$ as a function of $m_{\mathrm{H_1}}$ for $m_{\mathrm{SUSY}}=2600\;\mathrm{GeV}$.

Upper limits at 95% CL on $\sigma\times\mathcal{B}^2(\mathrm{H}_1\rightarrow b\bar{b})$ as a function of $m_{\mathrm{H_1}}$ for $m_{\mathrm{SUSY}}=2800\;\mathrm{GeV}$.

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 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 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 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 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 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 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 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 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 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 55 GeV mass and $ S \rightarrow b\bar{b}$ decay mode.

The cluster efficiency in bins of hadronic and EM energy in region A. Region A is defined as 391 cm $< r <$ 695.5 cm and 400 cm $< |z| <$ 671 cm. The cluster efficiency is estimated with LLPs decaying to $\tau^{+} \tau^{-}$. The sample contains equal fractions of events with LLP mass of 7, 15, 40, and 55 GeV and LLP lifetime of 0.1, 1, 10, and 100m. The first hadronic energy bins correspond to LLPs that decayed leptonically with 0 hadronic energy. The cluster efficiency includes all cluster-level selections described in the paper, except for the jet veto, time cut, and $\Delta\phi$ cut. The full simulation signal yield prediction for samples with various LLP mass between 7 - 55 GeV, lifetime between 0.1 - 100 m, and decay mode to $d\bar{d}$ and $\tau^{+} \tau^{-}$ can be reproduced using this parameterization to within 35% and 20% for region A and B, respectively.

The cluster efficiency in bins of hadronic and EM energy in region A. Region A is defined as 391 cm $< r <$ 695.5 cm and 400 cm $< |z| <$ 671 cm. The cluster efficiency is estimated with LLPs decaying to $\tau^{+} \tau^{-}$. The sample contains equal fractions of events with LLP mass of 7, 15, 40, and 55 GeV and LLP lifetime of 0.1, 1, 10, and 100m. The first hadronic energy bins correspond to LLPs that decayed leptonically with 0 hadronic energy. The cluster efficiency includes all cluster-level selections described in the paper, except for the jet veto, time cut, and $\Delta\phi$ cut. The full simulation signal yield prediction for samples with various LLP mass between 7 - 55 GeV, lifetime between 0.1 - 100 m, and decay mode to $d\bar{d}$ and $\tau^{+} \tau^{-}$ can be reproduced using this parameterization to within 35% and 20% for region A and B, respectively.

The cluster efficiency in bins of hadronic and EM energy in region B. Region B is defined as 671 cm $< |z| <$ 1100 cm, $r <$ 695.5 cm, and $|\eta| <$ 2. The cluster efficiency is estimated with LLPs decaying to $\tau^{+} \tau^{-}$. The sample contains equal fractions of events with LLP mass of 7, 15, 40, and 55 GeV and LLP lifetime of 0.1, 1, 10, and 100m. The first hadronic energy bins correspond to LLPs that decayed leptonically with 0 hadronic energy. The cluster efficiency includes all cluster-level selections described in the paper, except for the jet veto, time cut, and $\Delta\phi$ cut. The full simulation signal yield prediction for samples with various LLP mass between 7 - 55 GeV, lifetime between 0.1 - 100 m, and decay mode to $d\bar{d}$ and $\tau^{+} \tau^{-}$ can be reproduced using this parameterization to within 35% and 20% for region A and B, respectively.

The cluster efficiency in bins of hadronic and EM energy in region B. Region B is defined as 671 cm $< |z| <$ 1100 cm, $r <$ 695.5 cm, and $|\eta| <$ 2. The cluster efficiency is estimated with LLPs decaying to $\tau^{+} \tau^{-}$. The sample contains equal fractions of events with LLP mass of 7, 15, 40, and 55 GeV and LLP lifetime of 0.1, 1, 10, and 100m. The first hadronic energy bins correspond to LLPs that decayed leptonically with 0 hadronic energy. The cluster efficiency includes all cluster-level selections described in the paper, except for the jet veto, time cut, and $\Delta\phi$ cut. The full simulation signal yield prediction for samples with various LLP mass between 7 - 55 GeV, lifetime between 0.1 - 100 m, and decay mode to $d\bar{d}$ and $\tau^{+} \tau^{-}$ can be reproduced using this parameterization to within 35% and 20% for region A and B, respectively.

The efficiency of $N_{station} > 1$ requirement in bins of hadronic energy in region B. Region B is defined as 671 cm $< |z| <$ 1100 cm, $r <$ 695.5 cm, and $|\eta| <$ 2. The cluster efficiency is estimated with LLPs decaying to $\tau^{+} \tau^{-}$. The sample contains equal fractions of events with LLP mass of 7, 15, 40, and 55 GeV and LLP lifetime of 0.1, 1, 10, and 100m. The first hadronic energy bin corresponds to LLPs that decayed leptonically with 0 hadronic energy. The efficiency is calculated with respect to clusters that pass all cluster-level cuts described in the paper, except for the jet veto, time cut, and $\Delta\phi$ cut. The full simulation signal yield prediction for samples with various LLP mass between 7 - 55 GeV, lifetime between 0.1 - 100 m, and decay mode to $d\bar{d}$ and $\tau^{+} \tau^{-}$ can be reproduced using this parameterization to within 10%.

The efficiency of $N_{station} > 1$ requirement in bins of hadronic energy in region B. Region B is defined as 671 cm $< |z| <$ 1100 cm, $r <$ 695.5 cm, and $|\eta| <$ 2. The cluster efficiency is estimated with LLPs decaying to $\tau^{+} \tau^{-}$. The sample contains equal fractions of events with LLP mass of 7, 15, 40, and 55 GeV and LLP lifetime of 0.1, 1, 10, and 100m. The first hadronic energy bin corresponds to LLPs that decayed leptonically with 0 hadronic energy. The efficiency is calculated with respect to clusters that pass all cluster-level cuts described in the paper, except for the jet veto, time cut, and $\Delta\phi$ cut. The full simulation signal yield prediction for samples with various LLP mass between 7 - 55 GeV, lifetime between 0.1 - 100 m, and decay mode to $d\bar{d}$ and $\tau^{+} \tau^{-}$ can be reproduced using this parameterization to within 10%.

The geometric signal acceptance as a function of $c\tau$. The acceptance is shown for four different LLP mass hypotheses: 7, 15, 40, and 55 GeV. The acceptance is defined by requiring at least 1 LLP to decay in the region defined as 400 cm $< |z| <$ 1100 cm, $r < $ 695.5 cm, and $|\eta| <$ 2.4.