Signal acceptance times efficiency of HVT ZH(qqbb) events as a function of the resonance mass at different cut stages. Auxiliary table attached for 4 TeV mass point.

Dijet mass distributions in the WH 1-tag signal region. Uncertainties in the background and signal histograms are set to zero. The signal corresponds to that expected for HVT model B with resonance mass of 2 TeV.

Dijet mass distributions in the WH 2-tag signal region. Uncertainties in the background and signal histograms are set to zero. The signal corresponds to that expected for HVT model B with resonance mass of 2 TeV.

Dijet mass distributions in the ZH 1-tag signal region. Uncertainties in the background and signal histograms are set to zero. The signal corresponds to that expected for HVT model B with resonance mass of 2 TeV.

Dijet mass distributions in the ZH 2-tag signal region. Uncertainties in the background and signal histograms are set to zero. The signal corresponds to that expected for HVT model B with resonance mass of 2 TeV.

Distribution of $E^{miss}_T$ in data and for the expected SM background in the Two-Electron CR after performing the 'simplified shape fit'. Overflows are included in the fourth bin of each distribution. The $E^{miss}_T$ calculation in this CR does not include electron contribution. The error bars are statistical, and the dashed band includes statistical and systematic uncertainties determined by the fit. The lower panel shows the ratio of data to expected background event yields.

The observed exclusion contours on the $m_{\chi} - m_{med}$ plane for a simplified DM model involving an axial-vector mediator with couplings $g_{\chi}$ = 1, $g_q$ = 0.25 and $g_l$ = 0.

The expected exclusion contours on the $m_{\chi} - m_{med}$ plane for a simplified DM model involving an axial-vector mediator with couplings $g_{\chi}$ = 1, $g_q$ = 0.25 and $g_l$ = 0.

The expected + 1 sigma exclusion contours on the $m_{\chi} - m_{med}$ plane for a simplified DM model involving an axial-vector mediator with couplings $g_{\chi}$ = 1, $g_q$ = 0.25 and $g_l$ = 0.

The expected - 1 sigma exclusion contours on the $m_{\chi} - m_{med}$ plane for a simplified DM model involving an axial-vector mediator with couplings $g_{\chi}$ = 1, $g_q$ = 0.25 and $g_l$ = 0.

The observed exclusion contours on the $m_{\chi} - m_{med}$ plane for a simplified DM model involving an axial-vector mediator with couplings $g_{\chi}$ = 1, $g_q$ = 0.1 and $g_l$ = 0.1.

The expected exclusion contours on the $m_{\chi} - m_{med}$ plane for a simplified DM model involving an axial-vector mediator with couplings $g_{\chi}$ = 1, $g_q$ = 0.1 and $g_l$ = 0.1.

The expected + 1 sigma exclusion contours on the $m_{\chi} - m_{med}$ plane for a simplified DM model involving an axial-vector mediator with couplings $g_{\chi}$ = 1, $g_q$ = 0.1 and $g_l$ = 0.1.

The expected - 1 sigma exclusion contours on the $m_{\chi} - m_{med}$ plane for a simplified DM model involving an axial-vector mediator with couplings $g_{\chi}$ = 1, $g_q$ = 0.1 and $g_l$ = 0.1.

The observed exclusion contours on the $m_{\chi} - m_{med}$ plane for a vector mediator with couplings $g_{\chi}$ = 1, $g_q$ = 0.25 and $g_l$ = 0.

The expected exclusion contours on the $m_{\chi} - m_{med}$ plane for a vector mediator with couplings $g_{\chi}$ = 1, $g_q$ = 0.25 and $g_l$ = 0.

The expected + 1 sigma exclusion contours on the $m_{\chi} - m_{med}$ plane for a vector mediator with couplings $g_{\chi}$ = 1, $g_q$ = 0.25 and $g_l$ = 0.

The expected - 1 sigma exclusion contours on the $m_{\chi} - m_{med}$ plane for a vector mediator with couplings $g_{\chi}$ = 1, $g_q$ = 0.25 and $g_l$ = 0.

The observed exclusion contours on the $m_{\chi} - m_{med}$ plane for a vector mediator with couplings and $g_{\chi}$ = 1, $g_q$ = 0.1 and $g_l$ = 0.01

The expected exclusion contours on the $m_{\chi} - m_{med}$ plane for a vector mediator with couplings $g_{\chi}$ = 1, $g_q$ = 0.1 and $g_l$ = 0.01.

The expected + 1 sigma exclusion contours on the $m_{\chi} - m_{med}$ plane for a vector mediator with couplings $g_{\chi}$ = 1, $g_q$ = 0.1 and $g_l$ = 0.01.

The expected - 1 sigma exclusion contours on the $m_{\chi} - m_{med}$ plane for a vector mediator with couplings $g_{\chi}$ = 1, $g_q$ = 0.1 and $g_l$ = 0.01.

The 90% CL exclusion limit on $\chi$-proton spin-dependent scattering cross section in an axial-vector model with couplings $g_{\chi}$ = 1, $g_q$ = 0.25 and $g_l$ = 0 as a function of the dark-matter mass $m_{\chi}$. The 90% CL exclusion limit is displayed only in the regime where the theory is perturbative.

The 90% CL exclusion limit on $\chi$-neutron spin-dependent scattering cross section in an axial-vector model with couplings $g_{\chi}$ = 1, $g_q$ = 0.25 and $g_l$ = 0 as a function of the dark-matter mass $m_{\chi}$. The 90% CL exclusion limit is displayed only in the regime where the theory is perturbative.

The 90% CL exclusion limit on $\chi$-nucleon spin-independent scattering cross section in a vector model with couplings $g_{\chi}$ = 1, $g_q$ = 0.25 and $g_l$ = 0 as a function of the dark-matter mass $m_{\chi}$.

Observed and expected exclusion at 95% CL on the coupling $c_W$ as a function of the effective scale $f_a$ for an ALP mass of 1 MeV.

The observed 95% CL upper limits on signal strength $\mu$, the ratio of the experimental limit to the predicted signal cross section, for a simplified model of dark matter production involving an axial-vector operator and couplings $g_{\chi}$ = 1, $g_q$ = 0.25 and $g_l$ = 0 as a function of the mass of the dark matter, $m_{\chi}$ and the mass of the axial-vector mediator, $m_{med}$.

The observed 95% CL upper limits on cross section for a simplified model of dark matter production involving an axial-vector operator and couplings $g_{\chi}$ = 1, $g_q$ = 0.25 and $g_l$ = 0 as a function of the mass of the dark matter, $m_{\chi}$ and the mass of the axial-vector mediator, $m_{med}$.

Acceptance in % for a simplified model of dark matter production involving an axial-vector operator with couplings $g_{\chi}$ = 1, $g_q$ = 0.25 and $g_l$ = 0 for different values of $m_{\chi}$ and $m_{med}$ in the four inclusive signal regions.

Acceptance in % for a simplified model of dark matter production involving an axial-vector operator with couplings $g_{\chi}$ = 1, $g_q$ = 0.25 and $g_l$ = 0 for different values of $m_{\chi}$ and $m_{med}$ in the three exclusive signal regions.

Efficiencies for a simplified model of dark matter production involving an axial-vector operator with couplings $g_{\chi}$ = 1, $g_q$ = 0.25 and $g_l$ = 0 for different values of $m_{\chi}$ and $m_{med}$ in the four inclusive signal regions.

Efficiencies for a simplified model of dark matter production involving an axial-vector operator with couplings $g_{\chi}$ = 1, $g_q$ = 0.25 and $g_l$ = 0 for different values of $m_{\chi}$ and $m_{med}$ in the three exclusive signal regions.

The measured combined normalized differential cross-sections, divided by the bin-width, in the fiducial volume at Born level as well as a factor $k_{dressed}$ to translate from the Born particle level to the dressed particle level.

The measured combined normalized differential cross-sections, divided by the bin-width, in the fiducial volume at Born level as well as a factor $k_{dressed}$ to translate from the Born particle level to the dressed particle level.

Results of the normalized differential cross-section $1/\sigma_\mathrm{fid}\times \mathrm{d}\sigma_\mathrm{fid}/\mathrm{d}p_{T}^{ll}$ measured on born level for the $Z\rightarrow\mu\mu$ decay channel. The following naming convention is used: Stat.(Data), Stat.(MC) an Eff.(Uncor.), denote the statistical uncertainties due limited data and MC as well as the uncorrelated lepton efficiency uncertainties; Scale and Res. denote the muon momentum scale and resolution uncertainties; Muon Sag. denotes the uncertainty due to the muon sagitta bias; Eff. (Cor.), Isolation, Trigger and TTVA denote the uncertainties of the data/MC scale-factors for the correlated muon reconstruction, isolation, trigger and track-to-vertex matching efficiencies; the uncertainties due to the primary vertex z-distribution and pile-up reweighting are denoted as Z-Pos and Pile-Up, while the model and background uncertainties are summarized under Model and Bkg.. The sign-information is kept to track bin-to-bin changes.

Results of the normalized differential cross-section $1/\sigma_\mathrm{fid}\times \mathrm{d}\sigma_\mathrm{fid}/\mathrm{d}\phi_{\eta}^{*}$ measured on born level for the $Z\rightarrow\mu\mu$ decay channel. The following naming convention is used: Stat.(Data), Stat.(MC) an Eff.(Uncor.), denote the statistical uncertainties due limited data and MC as well as the uncorrelated lepton efficiency uncertainties; Scale and Res. denote the muon momentum scale and resolution uncertainties; Muon Sag. denotes the uncertainty due to the muon sagitta bias; Eff. (Cor.), Isolation, Trigger and TTVA denote the uncertainties of the data/MC scale-factors for the correlated muon reconstruction, isolation, trigger and track-to-vertex matching efficiencies; the uncertainties due to the primary vertex z-distribution and pile-up reweighting are denoted as Z-Pos and Pile-Up, while the model and background uncertainties are summarized under Model and Bkg.. The sign-information is kept to track bin-to-bin changes.

Results of the normalized differential cross-section $1/\sigma_\mathrm{fid}\times \mathrm{d}\sigma_\mathrm{fid}/\mathrm{d}p_{T}^{ll}$ measured on born level for the $Z\rightarrow ee$ decay channel. The following naming convention is used: Stat.(Data), Stat.(MC) and Eff.(Uncor.), denote the statistical uncertainties due limited data and MC as well as the uncorrelated lepton efficiency uncertainties; Scale and Res. for the electron momentum scale and resolution uncertainties; Elec. (Reco), Elec. (ID), Isolation, Trigger and Charge-ID denote the correlated uncertainties of the data/MC scale-factors for the electron reconstruction, identification, isolation, trigger and charge-identification efficiencies; The uncertainties due to the primary vertex z-distribution and pile-up reweighting are denoted as Z-Pos and Pile-Up, while the model and background uncertainties are summarized under Model and Bkg. The sign-information is kept to track bin-to-bin changes.

Results of the normalized differential cross-section $1/\sigma_\mathrm{fid} \times \mathrm{d}\sigma_\mathrm{fid}/\mathrm{d}\phi_{\eta}^{*}$ measured on born level for the $Z\rightarrow ee$ decay channel. The following naming convention is used: Stat.(Data), Stat.(MC) and Eff.(Uncor.), denote the statistical uncertainties due limited data and MC as well as the uncorrelated lepton efficiency uncertainties; Scale and Res. for the electron momentum scale and resolution uncertainties; Elec. (Reco), Elec. (ID), Isolation, Trigger and Charge-ID denote the correlated uncertainties of the data/MC scale-factors for the electron reconstruction, identification, isolation, trigger and charge-identification efficiencies; The uncertainties due to the primary vertex z-distribution and pile-up reweighting are denoted as Z-Pos and Pile-Up, while the model and background uncertainties are summarized under Model and Bkg. The sign-information is kept to track bin-to-bin changes.

Results of the normalized differential cross-section $1/\sigma_\mathrm{fid} \times \mathrm{d}\sigma_\mathrm{fid}/\mathrm{d}p_{T}^{ll}$ measured on bare level for the $Z\rightarrow\mu\mu$ decay channel. The following naming convention is used: Stat.(Data), Stat.(MC) an Eff.(Uncor.), denote the statistical uncertainties due limited data and MC as well as the uncorrelated lepton efficiency uncertainties; Scale and Res. denote the muon momentum scale and resolution uncertainties; Muon Sag. denotes the uncertainty due to the muon sagitta bias; Eff. (Cor.), Isolation, Trigger and TTVA denote the uncertainties of the data/MC scale-factors for the correlated muon reconstruction, isolation, trigger and track-to-vertex matching efficiencies; the uncertainties due to the primary vertex z-distribution and pile-up reweighting are denoted as Z-Pos and Pile-Up, while the model and background uncertainties are summarized under Model and Bkg.. The sign-information is kept to track bin-to-bin changes.

Results of the normalized differential cross-section $1/\sigma_\mathrm{fid} \times \mathrm{d}\sigma_\mathrm{fid}/\mathrm{d}\phi_{\eta}^{*}$ measured on bare level for the $Z\rightarrow\mu\mu$ decay channel. The following naming convention is used: Stat.(Data), Stat.(MC) an Eff.(Uncor.), denote the statistical uncertainties due limited data and MC as well as the uncorrelated lepton efficiency uncertainties; Scale and Res. denote the muon momentum scale and resolution uncertainties; Muon Sag. denotes the uncertainty due to the muon sagitta bias; Eff. (Cor.), Isolation, Trigger and TTVA denote the uncertainties of the data/MC scale-factors for the correlated muon reconstruction, isolation, trigger and track-to-vertex matching efficiencies; the uncertainties due to the primary vertex z-distribution and pile-up reweighting are denoted as Z-Pos and Pile-Up, while the model and background uncertainties are summarized under Model and Bkg.. The sign-information is kept to track bin-to-bin changes.

Measured combined normalized differential cross-section in the fiducial volume at Born level as well as a factor $k_{dressed}$ to translate from the Born particle level to the dressed particle.

Measured combined normalized differential cross-section in the fiducial volume at Born level as well as a factor $k_{dressed}$ to translate from the Born particle level to the dressed particle.

The distribution of events passing the selection requirements in the electron channel as a function of dilepton invariant mass $m_{ll}$ , the latter with one entry for each lepton per event. The MC signal sample is simulated using Powheg+Pythia8. The predictions of the MC signal sample together with the MC background samples are normalized to the integral of the data and the total experimental uncertainty of the predicted values is shown as a grey band in the ratio of the prediction to data.

The distribution of events passing the selection requirements in the muon channel as a function of dilepton invariant mass $m_{ll}$, the latter with one entry for each lepton per event. The MC signal sample is simulated using Powheg+Pythia8. The predictions of the MC signal sample together with the MC background samples are normalized to the integral of the data and the total experimental uncertainty of the predicted values is shown as a grey band in the ratio of the prediction to data.

The distribution of events passing the selection requirements in the electron channel as a function of lepton pseudorapidity $\eta$, the latter with one entry for each lepton per event. The MC signal sample is simulated using Powheg+Pythia8. The predictions of the MC signal sample together with the MC background samples are normalized to the integral of the data and the total experimental uncertainty of the predicted values is shown as a grey band in the ratio of the prediction to data.

The distribution of events passing the selection requirements in the muon channel as a function of lepton pseudorapidity $\eta$, the latter with one entry for each lepton per event. The MC signal sample is simulated using Powheg+Pythia8. The predictions of the MC signal sample together with the MC background samples are normalized to the integral of the data and the total experimental uncertainty of the predicted values is shown as a grey band in the ratio of the prediction to data.

The distribution of events passing the selection requirements in the electron channel as a function of dilepton transverse momentum. The MC signal sample is simulated using Powheg+Pythia8. The predictions are normalized to the integral of the data and the total experimental uncertainty of the predicted values is shown as a grey band in the ratio of the prediction to data.

The distribution of events passing the selection requirements in the muon channel as a function of dilepton transverse momentum. The MC signal sample is simulated using Powheg+Pythia8. The predictions are normalized to the integral of the data and the total experimental uncertainty of the predicted values is shown as a grey band in the ratio of the prediction to data.

The distribution of events passing the selection requirements in the electron channel as a function of $\phi_{\eta}^{*}$. The MC signal sample is simulated using Powheg+Pythia8. The predictions are normalized to the integral of the data and the total experimental uncertainty of the predicted values is shown as a grey band in the ratio of the prediction to data.

The distribution of events passing the selection requirements in the muon channel as a function of $\phi_{\eta}^{*}$. The MC signal sample is simulated using Powheg+Pythia8. The predictions are normalized to the integral of the data and the total experimental uncertainty of the predicted values is shown as a grey band in the ratio of the prediction to data.

The measured normalized cross section as a function of $p_{ll}$ for the electron and muon channels and the combined result as well as their ratio together with the total uncertainties, shown as a blue band. The pull distribution between the electron and muon channels, defined as the difference between the two channels divided by the combined uncorrelated uncertainty, is also shown. The $p_{ll}$ distribution is split into linear and logarithmic scales at 30 GeV.

The measured normalized cross section as a function of $\phi_{\eta}^{*}$ for the electron and muon channels and the combined result as well as their ratio together with the total uncertainties, shown as a blue band. The pull distribution between the electron and muon channels, defined as the difference between the two channels divided by the combined uncorrelated uncertainty, is also shown.

Comparison of the normalized $p_{ll}$ distributions predicted by different computations: Pythia8 with the AZ tune, Powheg+Pythia8 with the AZNLO tune, Sherpa v2.2.1 and RadISH with the Born level combined measurement. The uncertainties of the measurement are shown as vertical bars and uncertainties of the Sherpa and RadISH predictions are indicated by the coloured bands.

Comparison of the normalized $\phi_{\eta}^{*}$ distributions predicted by different computations: Pythia8 with the AZ tune, Powheg+Pythia8 with the AZNLO tune, Sherpa v2.2.1 and RadISH with the Born level combined measurement. The uncertainties of the measurement are shown as vertical bars and uncertainties of the Sherpa and RadISH predictions are indicated by the coloured bands.

Comparison of the normalized $p_{ll}$ distribution in the range $p_{ll}$ > 10 GeV. The Born level combined measurement is compared with predictions by Sherpa v2.2.1, fixed-order NNLOjet and NNLOjet supplied with NLO electroweak corrections. The uncertainties in the measurement are shown as vertical bars and the uncertainties in the predictions are indicated by the coloured bands.

The measured combined normalized differential cross-sections, divided by the bin-width, in the fiducial volume at dressed level.

The measured combined normalized differential cross-sections, divided by the bin-width, in the fiducial volume at dressed level.

Distribution of $S$, which is one of the discriminating variables used to define the $2\ell^{sc}$ SRs. The events are selected with the preselection requirements listed in Table 4 in the paper. The data (dots) are compared with the expected contributions from the relevant background sources (histograms). The expected signal distributions for $m_{H^{\pm\pm}} = 300~GeV$ are also shown, scaled to the observed number of events. The last bin includes overflows.

Distribution of $E_{T}^{miss}$, which is one of the discriminating variables used to define the $3\ell$ SRs. The events are selected with the preselection requirements listed in Table 4 in the paper. The data (dots) are compared with the expected contributions from the relevant background sources (histograms). The expected signal distributions for $m_{H^{\pm\pm}} = 300~GeV$ are also shown, scaled to the observed number of events. The last bin includes overflows.

Distribution of $\Delta R_{\ell^{\pm}\ell^{\pm}}$, which is one of the discriminating variables used to define the $3\ell$ SRs. The events are selected with the preselection requirements listed in Table 4 in the paper. The data (dots) are compared with the expected contributions from the relevant background sources (histograms). The expected signal distributions for $m_{H^{\pm\pm}} = 300~GeV$ are also shown, scaled to the observed number of events. The last bin includes overflows.

Distribution of $m_{x\ell}$ ($x$=3), which is one of the discriminating variables used to define the $3\ell$ SRs. The events are selected with the preselection requirements listed in Table 4 in the paper. The data (dots) are compared with the expected contributions from the relevant background sources (histograms). The expected signal distributions for $m_{H^{\pm\pm}} = 300~GeV$ are also shown, scaled to the observed number of events. The last bin includes overflows.

Distribution of $p_{T}^{leading jet}$, which is one of the discriminating variables used to define the $3\ell$ SRs. The events are selected with the preselection requirements listed in Table 4 in the paper. The data (dots) are compared with the expected contributions from the relevant background sources (histograms). The expected signal distributions for $m_{H^{\pm\pm}} = 300~GeV$ are also shown, scaled to the observed number of events. The last bin includes overflows.

Distribution of $E_{T}^{miss}$, which is one of the discriminating variables used to define the $4\ell$ SRs. The events are selected with the preselection requirements listed in Table 4 in the paper. The data (dots) are compared with the expected contributions from the relevant background sources (histograms). The expected signal distributions for $m_{H^{\pm\pm}} = 300~GeV$ are also shown, scaled to the observed number of events. The last bin includes overflows.

Distribution of $\Delta R_{\ell^{\pm}\ell^{\pm}}^{min}$, which is one of the discriminating variables used to define the $4\ell$ SRs. The events are selected with the preselection requirements listed in Table 4 in the paper. The data (dots) are compared with the expected contributions from the relevant background sources (histograms). The expected signal distributions for $m_{H^{\pm\pm}} = 300~GeV$ are also shown, scaled to the observed number of events. The last bin includes overflows.

Distribution of $m_{x\ell}$ ($x$=4), which is one of the discriminating variables used to define the $4\ell$ SRs. The events are selected with the preselection requirements listed in Table 4 in the paper. The data (dots) are compared with the expected contributions from the relevant background sources (histograms). The expected signal distributions for $m_{H^{\pm\pm}} = 300~GeV$ are also shown, scaled to the observed number of events. The last bin includes overflows.

Distribution of $p_{T}^{\ell_{1}}$, which is one of the discriminating variables used to define the $4\ell$ SRs. The events are selected with the preselection requirements listed in Table 4 in the paper. The data (dots) are compared with the expected contributions from the relevant background sources (histograms). The expected signal distributions for $m_{H^{\pm\pm}} = 300~GeV$ are also shown, scaled to the observed number of events. The last bin includes overflows.

The $E_{T}^{miss}$ distribution for the SRs of the $m_{H^{\pm\pm}} = 300~GeV$ signal mass hypothesis, where the selection requirement on $E_{T}^{miss}$ has been removed. In the attached plot, the signals are stacked on top of the backgrounds while individuals contributions of the $2\ell^{sc}$ channel are shown in HEPData. The last bin, isolated by a vertical red dashed line, is inclusive and corresponds to the SR.

The $E_{T}^{miss}$ distribution for the SRs of the $m_{H^{\pm\pm}} = 300~GeV$ signal mass hypothesis, where the selection requirement on $E_{T}^{miss}$ has been removed. In the attached plot, the signals are stacked on top of the backgrounds while individuals contributions of the $3\ell$ channel are shown in HEPData. The last bin, isolated by a vertical red dashed line, is inclusive and corresponds to the SR.

The $E_{T}^{miss}$ distribution for the SRs of the $m_{H^{\pm\pm}} = 300~GeV$ signal mass hypothesis, where the selection requirement on $E_{T}^{miss}$ has been removed. In the attached plot, the signals are stacked on top of the backgrounds while individuals contributions of the $4\ell$ channel are shown in HEPData. The last bin, isolated by a vertical red dashed line, is inclusive and corresponds to the SR.

Distribution of $v_{3}$ from MBT events plotted as a function of the A-particle $p_\mathrm{T}$ for 0-5% centrality.

Distribution of $v_{3}$ from $p_{T}^{jet}>75$ GeV events plotted as a function of the A-particle $p_\mathrm{T}$ for 0-5% centrality.

Distribution of $v_{3}$ from $p_{T}^{jet}>100$ GeV events plotted as a function of the A-particle $p_\mathrm{T}$ for 0-5% centrality.

Distribution of $v_{2}$ from MBT events plotted as a function of the event centrality and A-particle $p_\mathrm{T}$ in (0.5-2 GeV).

Distribution of $v_{2}$ from $p_{T}^{jet}>75$ GeV events plotted as a function of the event centrality and A-particle $p_\mathrm{T}$ in (0.5-2 GeV).

Distribution of $v_{2}$ from $p_{T}^{jet}>100$ GeV events plotted as a function of the event centrality and A-particle $p_\mathrm{T}$ in (0.5-2 GeV).

Distribution of $v_{2}$ from MBT events plotted as a function of the event centrality and A-particle $p_\mathrm{T}$ in (2-9 GeV).

Distribution of $v_{2}$ from $p_{T}^{jet}>75$ GeV events plotted as a function of the event centrality and A-particle $p_\mathrm{T}$ in (2-9 GeV).

Distribution of $v_{2}$ from $p_{T}^{jet}>100$ GeV events plotted as a function of the event centrality and A-particle $p_\mathrm{T}$ in (2-9 GeV).

Distribution of $v_{2}$ from MBT events plotted as a function of the event centrality and A-particle $p_\mathrm{T}$ in (9-100 GeV).

Distribution of $v_{2}$ from $p_{T}^{jet}>75$ GeV events plotted as a function of the event centrality and A-particle $p_\mathrm{T}$ in (9-100 GeV).

Distribution of $v_{2}$ from $p_{T}^{jet}>100$ GeV events plotted as a function of the event centrality and A-particle $p_\mathrm{T}$ in (9-100 GeV).

Distribution of relative UE-UE pair fractions from MBT events plotted as a function of A-particle $p_\mathrm{T}$ for 0-5% centrality.

Distribution of relative HS-HS pair fractions from MBT events plotted as a function of A-particle $p_\mathrm{T}$ for 0-5% centrality.

Distribution of relative HS-UE pair fractions from MBT events plotted as a function of A-particle $p_\mathrm{T}$ for 0-5% centrality.

Distribution of relative UE-HS pair fractions from MBT events plotted as a function of A-particle $p_\mathrm{T}$ for 0-5% centrality.

Distribution of relative UE-UE pair fractions from $p_{T}^{jet}>75$ GeV events plotted as a function of A-particle $p_\mathrm{T}$ for 0-5% centrality.

Distribution of relative HS-HS pair fractions from $p_{T}^{jet}>75$ GeV events plotted as a function of A-particle $p_\mathrm{T}$ for 0-5% centrality.

Distribution of relative HS-UE pair fractions from $p_{T}^{jet}>75$ GeV events plotted as a function of A-particle $p_\mathrm{T}$ for 0-5% centrality.

Distribution of relative UE-HS pair fractions from $p_{T}^{jet}>75$ GeV events plotted as a function of A-particle $p_\mathrm{T}$ for 0-5% centrality.

Distribution of relative UE-UE pair fractions from $p_{T}^{jet}>100$ GeV events plotted as a function of A-particle $p_\mathrm{T}$ for 0-5% centrality.

Distribution of relative HS-HS pair fractions from $p_{T}^{jet}>100$ GeV events plotted as a function of A-particle $p_\mathrm{T}$ for 0-5% centrality.

Distribution of relative HS-UE pair fractions from $p_{T}^{jet}>100$ GeV events plotted as a function of A-particle $p_\mathrm{T}$ for 0-5% centrality.

Distribution of relative UE-HS pair fractions from $p_{T}^{jet}>100$ GeV events plotted as a function of A-particle $p_\mathrm{T}$ for 0-5% centrality.

Distribution of relative UE-UE pair fractions from MBT events plotted as a function of centrality A-particle $p_\mathrm{T}$ in (0.5-2 GeV).

Distribution of relative HS-UE pair fractions from MBT events plotted as a function of centrality A-particle $p_\mathrm{T}$ in (0.5-2 GeV).

Distribution of relative UE-UE pair fractions from $p_{T}^{jet}>75$ GeV events plotted as a function of centrality A-particle $p_\mathrm{T}$ in (0.5-2 GeV).

Distribution of relative HS-UE pair fractions from $p_{T}^{jet}>75$ GeV events plotted as a function of centrality A-particle $p_\mathrm{T}$ in (0.5-2 GeV).

Distribution of relative UE-UE pair fractions from $p_{T}^{jet}>100$ GeV events plotted as a function of centrality A-particle $p_\mathrm{T}$ in (0.5-2 GeV).

Distribution of relative HS-UE pair fractions from $p_{T}^{jet}>100$ GeV events plotted as a function of centrality A-particle $p_\mathrm{T}$ in (0.5-2 GeV).

Distribution of relative UE-UE pair fractions from MBT events plotted as a function of centrality A-particle $p_\mathrm{T}$ in (2-9 GeV).

Distribution of relative HS-UE pair fractions from MBT events plotted as a function of centrality A-particle $p_\mathrm{T}$ in (2-9 GeV).

Distribution of relative UE-UE pair fractions from $p_{T}^{jet}>75$ GeV events plotted as a function of centrality A-particle $p_\mathrm{T}$ in (2-9 GeV).

Distribution of relative HS-UE pair fractions from $p_{T}^{jet}>75$ GeV events plotted as a function of centrality A-particle $p_\mathrm{T}$ in (2-9 GeV).

Distribution of relative UE-UE pair fractions from $p_{T}^{jet}>100$ GeV events plotted as a function of centrality A-particle $p_\mathrm{T}$ in (2-9 GeV).

Distribution of relative HS-UE pair fractions from $p_{T}^{jet}>100$ GeV events plotted as a function of centrality A-particle $p_\mathrm{T}$ in (2-9 GeV).

Distribution of relative UE-UE pair fractions from MBT events plotted as a function of centrality A-particle $p_\mathrm{T}$ in (9-100 GeV).

Distribution of relative HS-UE pair fractions from MBT events plotted as a function of centrality A-particle $p_\mathrm{T}$ in (9-100 GeV).

Distribution of relative UE-UE pair fractions from $p_{T}^{jet}>75$ GeV events plotted as a function of centrality A-particle $p_\mathrm{T}$ in (9-100 GeV).

Distribution of relative HS-UE pair fractions from $p_{T}^{jet}>75$ GeV events plotted as a function of centrality A-particle $p_\mathrm{T}$ in (9-100 GeV).

Distribution of relative UE-UE pair fractions from $p_{T}^{jet}>100$ GeV events plotted as a function of centrality A-particle $p_\mathrm{T}$ in (9-100 GeV).

Distribution of relative HS-UE pair fractions from $p_{T}^{jet}>100$ GeV events plotted as a function of centrality A-particle $p_\mathrm{T}$ in (9-100 GeV).

Post-fit background and signal yields in the full signal region. The total systematic uncertainty differs from the sum in quadrature of the different uncertainties due to correlations. Q mis-id refers to the charge misassignment background. Mat. Conv. and Low $m_{\gamma^*}$ refer respectively to events with one non-prompt electron originating from photon conversion in the detector material and to events with a virtual photon leading to an $e^+e^-$ pair. HF $e$ (HF $\mu$) refers to events with one non-prompt electron (muon) from heavy-flavour hadron decay, LF refers to events with a lepton originating from light-meson decay, and 'Other $t\bar{t}X$' includes events coming from $t\bar{t}WZ$, $t\bar{t}Z$+jets, $t\bar{t}WH$, $t\bar{t}HH$.

Ranking of the nuisance parameters included in the fit according to their impact on the signal strength $\mu$. Only the 20 most highly ranked nuisance parameters are shown in the thumbnail figure. The empty blue rectangles correspond to the pre-fit impact on $\mu$ and the filled blue ones to the post-fit impact on $\mu$, both referring to the upper $x$-axis scale. The impact of each nuisance parameter, $\Delta\mu$, is computed by comparing the nominal best-fit value of $\mu$ with the result of the fit when fixing the nuisance parameter to its best-fit value, $\hat{\theta}$, shifted by its pre-fit (post-fit) uncertainties $\pm \Delta\theta$ ($\pm \Delta\hat{\theta}$). The black points show the pulls of the parameter relative to their nominal values, $\theta_0$. The nominal value for all parameters is $\theta_0 = 0$. The pulls of the most impactful parameters and their relative post-fit errors, $\Delta\hat{\theta}/\Delta\theta$, are referenced to the lower $x$-axis scale. In the table, POI means Parameter of Interest.

Comparison between data and prediction before the fit ('Pre-Fit') for the 'HF e' control region. Mat. Conv. and Low $m_{\gamma^*}$ refer respectively to events with one non-prompt electron originating from photon conversion in the detector material and to events with a virtual photon leading to an $e^+e^-$ pair.

Comparison between data and prediction before the fit ('Pre-Fit') for the 'HF mu' control region. Mat. Conv. and Low $m_{\gamma^*}$ refer respectively to events with one non-prompt electron originating from photon conversion in the detector material and to events with a virtual photon leading to an $e^+e^-$ pair.

Comparison between data and prediction before the fit ('Pre-Fit') for the $m_{ee}^{PV}$ variable in the 'Conv.' control region. The first and last bins contain underflow and overflow events, respectively. Q mis-id refers to the charge misassignment background. Mat. Conv. and Low $m_{\gamma^*}$ refer respectively to events with one non-prompt electron originating from photon conversion in the detector material and to events with a virtual photon leading to an $e^+e^-$ pair.

Comparison between data and prediction before the fit ('Pre-Fit') for the $\Sigma p_{T}^{l}$ variable in the 'ttW' control region. The first and last bins contain underflow and overflow events, respectively. Q mis-id refers to the charge misassignment background. Mat. Conv. and Low $m_{\gamma^*}$ refer respectively to events with one non-prompt electron originating from photon conversion in the detector material and to events with a virtual photon leading to an $e^+e^-$ pair.

Comparison between data and prediction after the fit ('Post-Fit') for the 'HF e' control region. Mat. Conv. and Low $m_{\gamma^*}$ refer respectively to events with one non-prompt electron originating from photon conversion in the detector material and to events with a virtual photon leading to an $e^+e^-$ pair.

Comparison between data and prediction after the fit ('Post-Fit') for the 'HF mu' control region. Mat. Conv. and Low $m_{\gamma^*}$ refer respectively to events with one non-prompt electron originating from photon conversion in the detector material and to events with a virtual photon leading to an $e^+e^-$ pair.

Comparison between data and prediction after the fit ('Post-Fit') for the $m_{ee}^{PV}$ variable in the 'Conv.' control region. The first and last bins contain underflow and overflow events, respectively. Q mis-id refers to the charge misassignment background. Mat. Conv. and Low $m_{\gamma^*}$ refer respectively to events with one non-prompt electron originating from photon conversion in the detector material and to events with a virtual photon leading to an $e^+e^-$ pair.

Comparison between data and prediction after the fit ('Post-Fit') for the $\Sigma p_{T}^{l}$ variable in the in 'ttW' control region. The first and last bins contain underflow and overflow events, respectively. Q mis-id refers to the charge misassignment background. Mat. Conv. and Low $m_{\gamma^*}$ refer respectively to events with one non-prompt electron originating from photon conversion in the detector material and to events with a virtual photon leading to an $e^+e^-$ pair.

Comparison between data and prediction before the fit ('Pre-Fit') for the BDT score in the signal region. The first and last bins contain underflow and overflow events, respectively. Q mis-id refers to the charge misassignment background. Mat. Conv. and Low $m_{\gamma^*}$ refer respectively to events with one non-prompt electron originating from photon conversion in the detector material and to events with a virtual photon leading to an $e^+e^-$ pair.

Comparison between data and prediction after the fit ('Post-Fit') for the BDT score in the signal region. The first and last bins contain underflow and overflow events, respectively. Q mis-id refers to the charge misassignment background. Mat. Conv. and Low $m_{\gamma^*}$ refer respectively to events with one non-prompt electron originating from photon conversion in the detector material and to events with a virtual photon leading to an $e^+e^-$ pair.

Expected and observed 95% CL upper limits on the production of a heavy, narrow-width, scalar resonance decaying to a pair of Higgs bosons ($X\rightarrow HH$). The final state used in the search consists of a boosted $b\bar{b}$ pair and a boosted hadronically decaying $\tau^{+}\tau^{-}$ pair, and the SM braching ratio of the Higgs boson are assumed. The $\pm 1\sigma$ and $\pm 2\sigma$ variations about the expected limit are indicated by the error bands. Two different requirements are applied on the visible mass of the two boosted Higgs boson candidates for the resonance mass hypotheses of 1.6 TeV and 2.5 TeV, leading to discontinuities in the limits (at 1.6 TeV, the difference between imposing no requirement and $m^{vis}_{HH}>900$ GeV is less than 1% though).

Event yields for a simulated $ZH(\rightarrow aa\rightarrow (b\bar{b})(b\bar{b}))$ sample with $m_a = 17.5\,\text{GeV}$. The signal sample is produced with cross section equals to the standard model $pp\to ZH$, i.e. $0.88\,\text{pb}$. Cut 0 corresponds to the initial number of events. Cut 1 requires the single lepton trigger. Cut 2 requires 2 identified leptons. Cut 3 requires the Z-boson mass window. Cut 4 requires 2 reconstructed $a\to b\bar{b}$ candidates. Cut 5a requires 2 identified $a\to b\bar{b}$ candidates in the 1HPC1LPC region. Cut 6a requires the 2 $a\to b\bar{b}$ candidates in the 1HPC1LPC region to be inside the Higgs mass window. Cut 5b requires 2 identified $a\to b\bar{b}$ candidates in the 2HPC region. Cut 6b requires the 2 $a\to b\bar{b}$ candidates in the 2HPC region to be inside the Higgs mass window.

Background yield table for Z+jets, $t\bar{t}$, and rare sources. Observed data yield. Signal $ZH(\rightarrow aa\rightarrow (b\bar{b})(b\bar{b}))$ yield with $m_a = 20\,\text{GeV}$. The signal sample is produced with cross section equals to the standard model $pp\to ZH$, i.e. $0.88\,\text{pb}$. The table includes the yields in two signal regions with leptons consistent with an on-shell Z-boson decay, one with 2 $a\to b\bar{b}$ candidates in the 2HPC region and one with 2 $a\to b\bar{b}$ candidates in the 1HPC1LPC region. The table also includes the yields in four control regions, one with leptons consistent with an on-shell Z-boson decay and 2 $a\to b\bar{b}$ candidates in the Low Purity Category (LPC), and three others where the leptons are not consistent an on-shell Z-boson decay.