The post-fit distribution of the $M(\ell\ell)$ variable is shown for the high-MET bin for the $\text{t}\bar{\text{t}}$ CR. Uncertainties include both the statistical and systematic components.

The post-fit distribution of the $M(\ell\ell)$ variable is shown for the low-MET bin for the WZ-enriched region. Uncertainties include both the statistical and systematic components.

The post-fit distribution of the $M(\ell\ell)$ variable is shown for the high-MET bin for the WZ-enriched region. Uncertainties include both the statistical and systematic components.

The post-fit distribution of the $M(\ell\ell)$ variable is shown for the high-MET bin for the SS CR. Uncertainties include both the statistical and systematic components.

The 2$\ell$-Ewk SR: the post-fit distribution of the $M(\ell\ell)$ variable is shown for the low-MET bin. Uncertainties include both the statistical and systematic components. The signal distributions overlaid on the plot are from the TChiWZ and the simplified higgsino models in the scenario where the product of $\widetilde{m}_{\tilde{\chi}^0_2}\widetilde{m}_{\tilde{\chi}^0_1}$ eigenvalues is positive and negative, respectively. The numbers after the model name in the legend indicate the mass of the NLSP and the mass splitting between the NLSP and LSP, in GeV.

The 2$\ell$-Ewk SR: the post-fit distribution of the $M(\ell\ell)$ variable is shown for the med-MET bin. Uncertainties include both the statistical and systematic components. The signal distributions overlaid on the plot are from the TChiWZ and the simplified higgsino models in the scenario where the product of $\widetilde{m}_{\tilde{\chi}^0_2}\widetilde{m}_{\tilde{\chi}^0_1}$ eigenvalues is positive and negative, respectively. The numbers after the model name in the legend indicate the mass of the NLSP and the mass splitting between the NLSP and LSP, in GeV.

The 2$\ell$-Ewk SR: the post-fit distribution of the $M(\ell\ell)$ variable is shown for the high-MET bin. Uncertainties include both the statistical and systematic components. The signal distributions overlaid on the plot are from the TChiWZ and the simplified higgsino models in the scenario where the product of $\widetilde{m}_{\tilde{\chi}^0_2}\widetilde{m}_{\tilde{\chi}^0_1}$ eigenvalues is positive and negative, respectively. The numbers after the model name in the legend indicate the mass of the NLSP and the mass splitting between the NLSP and LSP, in GeV.

The 2$\ell$-Ewk SR: the post-fit distribution of the $M(\ell\ell)$ variable is shown for the ultra-MET bin. Uncertainties include both the statistical and systematic components. The signal distributions overlaid on the plot are from the TChiWZ and the simplified higgsino models in the scenario where the product of $\widetilde{m}_{\tilde{\chi}^0_2}\widetilde{m}_{\tilde{\chi}^0_1}$ eigenvalues is positive and negative, respectively. The numbers after the model name in the legend indicate the mass of the NLSP and the mass splitting between the NLSP and LSP, in GeV.

The 3$\ell$-Ewk search regions: the post-fit distribution of the $M^{\text{min}}_{\text{SFOS}}(\ell\ell)$ variable is shown for the low-MET bin. Uncertainties include both the statistical and systematic components. The signal distributions overlaid on the plot are from the TChiWZ and the simplified higgsino models in the scenario where the product of $\widetilde{m}_{\tilde{\chi}^0_2}\widetilde{m}_{\tilde{\chi}^0_1}$ eigenvalues is positive and negative, respectively. The numbers after the model name in the legend indicate the mass of the NLSP and the mass splitting between the NLSP and LSP, in GeV.

The 3$\ell$-Ewk search regions: the post-fit distribution of the $M^{\text{min}}_{\text{SFOS}}(\ell\ell)$ variable is shown for the high-MET bin. Uncertainties include both the statistical and systematic components. The signal distributions overlaid on the plot are from the TChiWZ and the simplified higgsino models in the scenario where the product of $\widetilde{m}_{\tilde{\chi}^0_2}\widetilde{m}_{\tilde{\chi}^0_1}$ eigenvalues is positive and negative, respectively. The numbers after the model name in the legend indicate the mass of the NLSP and the mass splitting between the NLSP and LSP, in GeV.

The 2$\ell$-Stop SR: the post-fit distribution of the leading lepton $p_{T}$ variable is shown for the low-MET bin. Uncertainties include both the statistical and systematic components. The signal distributions overlaid on the plot are from the T2bff$\tilde{\chi}^0_1$ and the T2bW models. The numbers after the model name in the legend indicate the mass of the top squark and the mass splitting between the top squark and LSP, in GeV.

The 2$\ell$-Stop SR: the post-fit distribution of the leading lepton $p_{T}$ variable is shown for the med-MET bin. Uncertainties include both the statistical and systematic components. The signal distributions overlaid on the plot are from the T2bff$\tilde{\chi}^0_1$ and the T2bW models. The numbers after the model name in the legend indicate the mass of the top squark and the mass splitting between the top squark and LSP, in GeV.

The 2$\ell$-Stop SR: the post-fit distribution of the leading lepton $p_{T}$ variable is shown for the high-MET bin. Uncertainties include both the statistical and systematic components. The signal distributions overlaid on the plot are from the T2bff$\tilde{\chi}^0_1$ and the T2bW models. The numbers after the model name in the legend indicate the mass of the top squark and the mass splitting between the top squark and LSP, in GeV.

The 2$\ell$-Stop SR: the post-fit distribution of the leading lepton $p_{T}$ variable is shown for the ultra-MET bin. Uncertainties include both the statistical and systematic components. The signal distributions overlaid on the plot are from the T2bff$\tilde{\chi}^0_1$ and the T2bW models. The numbers after the model name in the legend indicate the mass of the top squark and the mass splitting between the top squark and LSP, in GeV.

The observed 95% CL exclusion contours (black curves) assuming the NLO+NLL cross sections, with the variations (thin lines) corresponding to the uncertainty in the cross section for the TChiWZ model. The red curves present the 95% CL expected limits with the band (thin lines) covering 68% of the limits in the absence of signal. Results are reported for the $\widetilde{m}_{\tilde{\chi}^0_2} \widetilde{m}_{\tilde{\chi}^0_1} > 0$ $M(\ell\ell)$ spectrum reweighting scenario. The range of luminosities of the analysis regions included in the fit is indicated on the plot.

The observed 95% CL exclusion contours (black curves) assuming the NLO+NLL cross sections, with the variations (thin lines) corresponding to the uncertainty in the cross section for the TChiWZ model. The red curves present the 95% CL expected limits with the band (thin lines) covering 68% of the limits in the absence of signal. Results are reported for the $\widetilde{m}_{\tilde{\chi}^0_2} \widetilde{m}_{\tilde{\chi}^0_1} < 0$ $M(\ell\ell)$ spectrum reweighting scenario. The range of luminosities of the analysis regions included in the fit is indicated on the plot.

The observed 95% CL exclusion contours (black curves) assuming the NLO+NLL cross sections, with the variations (thin lines) corresponding to the uncertainty in the cross section for the simplified higgsino model. The simplified model includes both neutralino pair and neutralino-chargino production modes. The red curves present the 95% CL expected limits with the band (thin lines) covering 68% of the limits in the absence of signal. The results are reported for the $\widetilde{m}_{\tilde{\chi}^0_2}\widetilde{m}_{\tilde{\chi}^0_1}$ $M(\ell\ell)$ spectrum reweighting scenario. The range of luminosities of the analysis regions included in the fit is indicated on the plot.

The observed 95% CL exclusion contours (black curves) assuming the NLO+NLL cross sections, with the variations (thin lines) corresponding to the uncertainty in the cross section for the pMSSM higgsino model. The pMSSM one includes all possible production modes. The red curves present the 95% CL expected limits with the band (thin lines) covering 68% of the limits in the absence of signal. The results are reported for the $\widetilde{m}_{\tilde{\chi}^0_2}\widetilde{m}_{\tilde{\chi}^0_1}$ $M(\ell\ell)$ spectrum reweighting scenario. The range of luminosities of the analysis regions included in the fit is indicated on the plot.

The observed 95% CL exclusion contours (black curves) assuming the NLO+NLL cross sections, with the variations (thin lines) corresponding to the uncertainty in the cross section for the T2bff$\tilde{\chi}^0_1$ simplified model. The red curves present the 95% CL expected limits with the band (thin lines) covering 68% of the limits in the absence of signal. The range of luminosities of the analysis regions included in the fit is indicated on the plot.

The observed 95% CL exclusion contours (black curves) assuming the NLO+NLL cross sections, with the variations (thin lines) corresponding to the uncertainty in the cross section for the T2bW simplified model. The red curves present the 95% CL expected limits with the band (thin lines) covering 68% of the limits in the absence of signal. The range of luminosities of the analysis regions included in the fit is indicated on the plot.

The background estimate and the observed data in the number of tagged jets >= 2 bin, for each of the seven validation samples (VS$_{1}$ through VS$_{7}$), along with the signal sample (Sig S). Signal model distributions for scalar masses of 15 and 55 GeV with a proper mean decay length of 20 mm are also shown. The Higgs boson branching fraction to long-lived scalars (B(H$\rightarrow $ SS)) is set to 20%.

Exclusion limits at 95% CL on the Higgs boson branching fraction to long-lived scalars B(H$\rightarrow$ SS) as a function of the mean proper decay length of the scalar for scalar mass 15 GeV and the S$\rightarrow d\bar{d}$ decay mode.

Exclusion limits at 95% CL on the Higgs boson branching fraction to long-lived scalars B(H$\rightarrow$ SS) as a function of the mean proper decay length of the scalar for scalar mass 40 GeV and the S$\rightarrow d\bar{d}$ decay mode.

Exclusion limits at 95% CL on the Higgs boson branching fraction to long-lived scalars B(H$\rightarrow$ SS) as a function of the mean proper decay length of the scalar for scalar mass 55 GeV and the S$\rightarrow d\bar{d}$ decay mode.

Exclusion limits at 95% CL on the Higgs boson branching fraction to long-lived scalars B(H$\rightarrow$ SS) as a function of the mean proper decay length of the scalar for scalar mass 15 GeV and the S$\rightarrow b\bar{b}$ decay mode.

Exclusion limits at 95% CL on the Higgs boson branching fraction to long-lived scalars B(H$\rightarrow$ SS) as a function of the mean proper decay length of the scalar for scalar mass 40 GeV and the S$\rightarrow b\bar{b}$ decay mode.

Exclusion limits at 95% CL on the Higgs boson branching fraction to long-lived scalars B(H$\rightarrow$ SS) as a function of the mean proper decay length of the scalar for scalar mass 55 GeV and the S$\rightarrow b\bar{b}$ decay mode.

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

95% CL expected and observed frequentist upper limit for the LbL cross section within the $m_{\gamma\gamma} > 350$ GeV, $0.070 < \xi^+ < 0.111$, and $0.070 < \xi^- < 0.138$ search region.

95% expected and observed one-dimensional limits on $\zeta_1$ and $\zeta_2$ anomalous LbyL production parameters, when the other parameter is set to zero. This corresponds to a search region of $m_{\gamma\gamma} > 350$ GeV, $0.070 < \xi^+ < 0.111$, and $0.070 < \xi^- < 0.138$.

95% expected and observed limits on $\zeta_1$ and $\zeta_2$ anomalous LbyL production parameters. The parametric elliptic form is assumed: $a_0\zeta_1^2+a_1\zeta_1\zeta_2+a_2\zeta_2^2=1$. This corresponds to a search region of $m_{\gamma\gamma} > 350$ GeV, $0.070 < \xi^+ < 0.111$, and $0.070 < \xi^- < 0.138$.

The distributions of the transfer factor ($\alpha$) versus $m_T$ in the HP VBF category is shown. The last bin corresponds to the value obtained by integrating events above the penultimate bin.

The distributions of the transfer factor ($\alpha$) versus $m_T$ in the HP ggF/DY category is shown. The last bin corresponds to the value obtained by integrating events above the penultimate bin.

The distributions of the transfer factor ($\alpha$) versus $m_T$ in the LP VBF category is shown. The last bin corresponds to the value obtained by integrating events above the penultimate bin.

The distributions of the transfer factor ($\alpha$) versus $m_T$ in the LP ggF/DY category is shown. The last bin corresponds to the value obtained by integrating events above the penultimate bin.

Distributions of mT for high-purity ggF/DY CR events after performing background-only fits. The last bin corresponds to the yields integrated above the penultimate bin. Both the ggF and VBF-produced 1 TeV graviton signals (normalized to 10 fb) are shown in the plot. Nonresonant background corresponds to the Z+jets and W+jets events. Resonant background corresponds to the ttbar, single top, and di/tri-boson events.

Distributions of mT for high-purity VBF CR events after performing background-only fits. The last bin corresponds to the yields integrated above the penultimate bin. Both the ggF and VBF-produced 1 TeV graviton signals (normalized to 10 fb) are shown in the plot. Nonresonant background corresponds to the Z+jets and W+jets events. Resonant background corresponds to the ttbar, single top, and di/tri-boson events.

Distributions of mT for low-purity ggF/DY CR events after performing background-only fits. The last bin corresponds to the yields integrated above the penultimate bin. Both the ggF and VBF-produced 1 TeV graviton signals (normalized to 10 fb) are shown in the plot. Nonresonant background corresponds to the Z+jets and W+jets events. Resonant background corresponds to the ttbar, single top, and di/tri-boson events.

Distributions of mT for low-purity VBF CR events after performing background-only fits. The last bin corresponds to the yields integrated above the penultimate bin. Both the ggF and VBF-produced 1 TeV graviton signals (normalized to 10 fb) are shown in the plot. Nonresonant background corresponds to the Z+jets and W+jets events. Resonant background corresponds to the ttbar, single top, and di/tri-boson events.

Distributions of mT for high-purity ggF/DY SR events after performing background-only fits. The last bin corresponds to the yields integrated above the penultimate bin. Both the ggF and VBF-produced 1 TeV graviton signals (normalized to 10 fb) are shown in the plot. Nonresonant background corresponds to the Z+jets and W+jets events. Resonant background corresponds to the ttbar, single top, and di/tri-boson events.

Distributions of mT for high-purity VBF SR events after performing background-only fits. The last bin corresponds to the yields integrated above the penultimate bin. Both the ggF and VBF-produced 1 TeV graviton signals (normalized to 10 fb) are shown in the plot. Nonresonant background corresponds to the Z+jets and W+jets events. Resonant background corresponds to the ttbar, single top, and di/tri-boson events.

Distributions of mT for low-purity ggF/DY SR events after performing background-only fits. The last bin corresponds to the yields integrated above the penultimate bin. Both the ggF and VBF-produced 1 TeV graviton signals (normalized to 10 fb) are shown in the plot. Nonresonant background corresponds to the Z+jets and W+jets events. Resonant background corresponds to the ttbar, single top, and di/tri-boson events.

Distributions of mT for low-purity VBF SR events after performing background-only fits. The last bin corresponds to the yields integrated above the penultimate bin. Both the ggF and VBF-produced 1 TeV graviton signals (normalized to 10 fb) are shown in the plot. Nonresonant background corresponds to the Z+jets and W+jets events. Resonant background corresponds to the ttbar, single top, and di/tri-boson events.

Expected and observed 95% CL upper limits on the product of the radion (R) production (gluon-gluon fusion) cross section and the R -> ZZ branching fraction versus the radion mass.

Expected and observed 95% CL upper limits on the product of the radion (R) production (vector boson fusion) cross section and the R -> ZZ branching fraction versus the radion mass.

Expected and observed 95% CL upper limits on the product of the W' (W') production (Drell-Yan) cross section and the W' -> WZ branching fraction versus the W' mass.

Expected and observed 95% CL upper limits on the product of the W' (W') production (vector boson fusion) cross section and the W' -> WZ branching fraction versus the W' mass.

Expected and observed 95% CL upper limits on the product of the graviton (G) production (gluon-gluon fusion) cross section and the G -> ZZ branching fraction versus the graviton mass.

Expected and observed 95% CL upper limits on the product of the graviton (G) production (vector boson fusion) cross section and the G -> ZZ branching fraction versus the graviton mass.

Post-fit $m_\text{T}(\gamma, E_\text{T}^\text{miss})$ distribution in the inclusive signal region for the dark-photon search with the 125 GeV mass $\mathcal{B}(H \to \gamma \gamma_\text{d})$ signal normalization set to zero. A $H \to \gamma \gamma_\text{d}$ decay signal is shown for two different mass hypotheses, 125 GeV and 500 GeV, and scaled to a $\mathcal{B}(H \to \gamma \gamma_\text{d})$ of 2% and 1%, respectively. Events with $m_\text{T}(\gamma, E_\text{T}^\text{miss})$ larger than the rightmost bin boundary are added to that bin.

The 95% CL upper limit on the Higgs boson production cross-section times branching ratio to $\gamma \gamma_\text{d}$ is shown for different VBF-produced scalar-mediator-mass hypotheses in the NWA. The theoretically predicted cross-section of a Higgs boson produced via VBF and with the $\mathcal{B}(H \to \gamma \gamma_\text{d}) =$ 5% is superimposed on the $\pm 1\sigma$ and $\pm 2\sigma$ NNLO QCD + NLO EW uncertainty band of the expected production cross-section limit.

Post-fit $m_\text{jj}$ distribution in the inclusive signal region. The Higgs boson invisible decay signal is scaled to a $\mathcal{B}_\text{inv}$ of 37%. Events with $m_\text{jj}$ larger than the rightmost bin boundary are added to that bin.

Post-fit $m_\text{jj}$ distribution in the one-lepton control region $W_{\ell \nu}^\gamma$ CR. Events with $m_\text{jj}$ larger than the rightmost bin boundary are added to that bin.

Post-fit $m_\text{T}$ distribution in the one lepton control region. Events with $m_\text{T}$ larger than the rightmost bin boundary are added to that bin.

Post-fit photon centrality distribution in the zero lepton signal plus control region with the $\mathcal{B}_\text{inv}$ signal normalization set to zero in the fit.

Post-fit photon $E_\text{T}$ distribution in the zero lepton signal region with the $\mathcal{B}_\text{inv}$ signal normalization set to zero in the fit.

Post-fit photon centrality distribution in the zero lepton signal plus control region resulting from the fit to the $m_\text{jj}$ distribution for EW $Z \gamma + \text{jets}$. The post-fit uncertainties include statistical, experimental, and theory contributions.

Post-fit photon $E_\text{T}$ distribution in the zero lepton signal region resulting from the fit to the $m_\text{jj}$ distribution for EW $Z \gamma + \text{jets}$. The post-fit uncertainties include statistical, experimental, and theory contributions.

Post-fit DNN output score distribution in the one lepton control region.

Yields for the EW $Z \gamma + \text{jets}$ process are shown after each selection along with relative and absolute signal acceptance efficiencies.

Yields for the 125 GeV Higgs boson with $\mathcal{B}_\text{inv.} =$ 1 signal produced by the vector boson fusion process in association with a final state photon are shown after each selection along with relative and absolute signal acceptance efficiencies.

Yields for the 125 GeV Higgs boson with $\mathcal{B}(H \to \gamma \gamma_\text{d}) =$ 1 signal produced by the vector boson fusion process are shown after each selection along with relative and absolute signal acceptance efficiencies.

Expected +- 1 sigma 95% CL exclusion limit for the Zprime-2HDM model.

Expected +- 2 sigma 95% CL exclusion limit for the Zprime-2HDM model.

Observed 95% CL exclusion limit for the 2HDM+a model ggF production.

Expected 95% CL exclusion limit for the 2HDM+a model ggF production.

Expected +- 1 sigma 95% CL exclusion limit for the 2HDM+a model ggF production.

Expected +- 2 sigma 95% CL exclusion limit for the 2HDM+a model ggF production.

Observed 95% CL exclusion limit for the 2HDM+a model bbA production.

Expected 95% CL exclusion limit for the 2HDM+a model bbA production.

Expected +- 1 sigma 95% CL exclusion limit for the 2HDM+a model bbA production.

Expected +- 2 sigma 95% CL exclusion limit for the 2HDM+a model bbA production.

Observed 95% CL exclusion limit for the Zprime-2HDM model with the benchmark used in arXiv:1707.01302.

Expected 95% CL exclusion limit for the Zprime-2HDM model with the benchmark used in arXiv:1707.01302.

Expected +- 1 sigma 95% CL exclusion limit for the Zprime-2HDM model with the benchmark used in arXiv:1707.01302.

Expected +- 2 sigma 95% CL exclusion limit for the Zprime-2HDM model with the benchmark used in arXiv:1707.01302.

Expected and observed upper limits at 95% CL on cross-section for Zprime-2HDM model.

Expected and observed upper limits at 95% CL on cross-section for ggF producton in the 2HDM+a model.

Expected and observed upper limits at 95% CL on cross-section for bbA producton in the 2HDM+a model.

Model-independent upper limits on the visible cross-section $σ_{vis, $h(\bar{b})+DM} ≡ σ_{h+DM} \times B(h \to b\bar{b}) \times \mathcal{A} \times \epsilon$ in the different signal regions.

Theory cross-section for Zprime-2HDM model.

Theory cross-section for bbA production in the 2HDM+a model.

Theory cross-section for ggF production in the 2HDM+a model.

Distribution of Higgs boson candidate mass in 2b region with MET=150-200 GeV.

Distribution of Higgs boson candidate mass in 2b region with MET=200-350 GeV.

Distribution of Higgs boson candidate mass in 2b region with MET=350-500 GeV.

Distribution of Higgs boson candidate mass in 2b region with MET=500-750 GeV.

Distribution of Higgs boson candidate mass in 2b region with MET > 750 GeV.

Distribution of Higgs boson candidate mass in 3b region with MET=150-200 GeV.

Distribution of Higgs boson candidate mass in 3b region with MET=200-350 GeV.

Distribution of Higgs boson candidate mass in 3b region with MET=350-500 GeV.

Distribution of Higgs boson candidate mass in 3b region with MET > 500 GeV.

Yields in 1-lepton control region.

Yields in 2-lepton control region.

MET distribution in 2b region of the 0-lepton channel.

MET distribution in 3b region of the 0-lepton channel.

Expected signal yields after certain selection cuts in 2b region with MET=150-200 GeV.

Expected signal yields after certain selection cuts in 2b region with MET=200-350 GeV.

Expected signal yields after certain selection cuts in 2b region with MET=350-500 GeV.

Expected signal yields after certain selection cuts in 2b region with MET=500-750 GeV.

Expected signal yields after certain selection cuts in 2b region with MET > 750 GeV.

Expected signal yields after certain selection cuts in 3b region with MET=150-200 GeV.

Expected signal yields after certain selection cuts in 3b region with MET=200-350 GeV.

Expected signal yields after certain selection cuts in 3b region with MET=350-500 GeV.

Expected signal yields after certain selection cuts in 3b region with MET > 500 GeV.

Acceptance times efficiency for bbA production in the 2HDM+a model - 2b region with MET=150-200 GeV.

Acceptance times efficiency for bbA production in the 2HDM+a model - 2b region with MET=200-350 GeV.

Acceptance times efficiency for bbA production in the 2HDM+a model - 2b region with MET=350-500 GeV.

Acceptance times efficiency for bbA production in the 2HDM+a model - 2b region with MET=500-750 GeV.

Acceptance times efficiency for bbA production in the 2HDM+a model - 2b region with MET > 750 GeV.

Acceptance times efficiency for bbA production in the 2HDM+a model - 3b region with MET=150-200 GeV.

Acceptance times efficiency for bbA production in the 2HDM+a model - 3b region with MET=200-350 GeV.

Acceptance times efficiency for bbA production in the 2HDM+a model - 3b region with MET=350-500 GeV.

Acceptance times efficiency for bbA production in the 2HDM+a model - 3b region with MET>500 GeV.

Acceptance times efficiency for ggF production in the 2HDM+a model - 2b region with MET=150-200 GeV.

Acceptance times efficiency for ggF production in the 2HDM+a model - 2b region with MET=200-350 GeV.

Acceptance times efficiency for ggF production in the 2HDM+a model - 2b region with MET=350-500 GeV.

Acceptance times efficiency for ggF production in the 2HDM+a model - 2b region with MET=500-750 GeV.

Acceptance times efficiency for ggF production in the 2HDM+a model - 2b region with MET > 750 GeV.

Acceptance times efficiency for ggF production in the 2HDM+a model - 3b region with MET=150-200 GeV.

Acceptance times efficiency for ggF production in the 2HDM+a model - 3b region with MET=200-350 GeV.

Acceptance times efficiency for ggF production in the 2HDM+a model - 3b region with MET=350-500 GeV.

Acceptance times efficiency for ggF production in the 2HDM+a model - 3b region with MET > 500 GeV.

Acceptance times efficiency for ggF production in the Zprime-2HDM model - 2b region with MET=150-200 GeV.

Acceptance times efficiency for ggF production in the Zprime-2HDM model - 2b region with MET=200-350 GeV.

Acceptance times efficiency for ggF production in the Zprime-2HDM model - 2b region with MET=350-500 GeV.

Acceptance times efficiency for ggF production in the Zprime-2HDM model - 2b region with MET=500-750 GeV.

Acceptance times efficiency for ggF production in the Zprime-2HDM model - 2b region with MET > 750 GeV.

Acceptance times efficiency for ggF production in the Zprime-2HDM model - 3b region with MET=150-200 GeV.

Acceptance times efficiency for ggF production in the Zprime-2HDM model - 3b region with MET=200-350 GeV.

Acceptance times efficiency for ggF production in the Zprime-2HDM model - 3b region with MET=350-500 GeV.

Acceptance times efficiency for ggF production in the Zprime-2HDM model - 3b region with MET > 500 GeV.

The measured differential branching fractions as a function of the invariant hadronic mass squared of the $X_u$ system ($M_X^2$).

The measured differential branching fractions as a function of the light-cone momenta of the hadronic $X_u$ system [$P_+ = (E_X^B - |old{p}_X^B|)$].

The measured differential branching fractions as a function of the light-cone momenta of the hadronic $X_u$ system [$P_- = (E_X^B + |old{p}_X^B|)$].

The background-subtracted yields as a function of $E_\ell^B$.

The background-subtracted yields as a function of $q^2$.

The background-subtracted yields as a function of $M_X$.

The background-subtracted yields as a function of $M_X^2$.

The background-subtracted yields as a function of $P_+$. [$P_+ = (E_X^B - |P_X^B|)$]

The background-subtracted yields as a function of $P_-$. [$P_- = (E_X^B + |P_X^B|)$].

The full experimental correlation matrix for the background-subtracted spectrum of the lepton energy in the $B$ rest frame ($E_\ell^B$).

The full experimental correlation matrix for the background-subtracted spectrum of $q^2$).

The full experimental correlation matrix for the background-subtracted spectrum of $M_X$.

The full experimental correlation matrix for the background-subtracted spectrum of $M_X^2$.

The full experimental correlation matrix for the background-subtracted spectrum of $P_+$.

The full experimental correlation matrix for the background-subtracted spectrum of $P_-$.

Migration matrix of $E_\ell^B$ [in %].

Migration matrix of $q^2$ [in %].

Migration matrix of $M_X$ [in %].

Migration matrix of $M_X^2$ [in %].

Migration matrix of $P_+$ [in %].

Migration matrix of $P_-$ [in %].

The total correction factors and phase space acceptance as a function of the invariant hadronic mass of the $X_u$ system ($M_X$).

The total correction factors and phase space acceptance as a function of the invariant hadronic mass squared of the $X_u$ system ($M_X^2$).

The total correction factors and phase space acceptance as a function of the lepton energy in the $B$ rest frame ($E_\ell^B$).

The total correction factors and phase space acceptance as a function of the four-momentum-transfer squared of the $B$ to the $X_u$ system ($q^2$).

The total correction factors and phase space acceptance as a function of $P_+$.

The total correction factors and phase space acceptance as a function of $P_-$.

The statistical correlations of the differential branching fractions are shown. The matrix elements are ordered for $M_{X}$, $M_{X}^{2}$, $q^{2}$, $E_{\ell}^{B}$, $P_{+}$ and $P_{-}$.

The full experimental (statistical and systematical) correlations of the differential branching fractions are shown. The matrix elements are ordered for $M_{X}$, $M_{X}^{2}$, $q^{2}$, $E_{\ell}^{B}$, $P_{+}$ and $P_{-}$.

The total partial branching fraction with $E_{\ell}^{B} > 1$ GeV. The order of entries is: the 2D fit result of Phys. Rev. D 104, 012008 (2021), the results calcuated with the differential branching fractions of $M_{X}$, $M_{X}^{2}$, $q^{2}$, $E_{\ell}^{B}$, $P_{+}$ and $P_{-}$

The full experimental correlations between the total partial branching fractions from summing the individual bins are shown. The matrix elements are saved in the order of $M_{X}$, $M_{X}^{2}$, $q^{2}$, $E_{\ell}^{B}$, $P_{+}$ and $P_{-}$.

The first moment of the measured differential branching fraction of the lepton energy in the $B$ rest frame ($E_\ell^B$).

The second moment of the measured differential branching fraction of the lepton energy in the $B$ rest frame ($E_\ell^B$).

The third moment of the measured differential branching fraction of the lepton energy in the $B$ rest frame ($E_\ell^B$).

The first moment of the measured differential branching fraction of $q^2$.

The second moment of the measured differential branching fraction of $q^2$.

The third moment of the measured differential branching fraction of $q^2$.

The first moment of the measured differential branching fraction of $M_{X}$.

The second moment of the measured differential branching fraction of $M_{X}$.

The third moment of the measured differential branching fraction of $M_{X}$.

The first moment of the measured differential branching fraction of $P_{+}$.

The second moment of the measured differential branching fraction of $P_{+}$.

The third moment of the measured differential branching fraction of $P_{+}$.

The first moment of the measured differential branching fraction of $P_{-}$.

The second moment of the measured differential branching fraction of $P_{-}$.

The third moment of the measured differential branching fraction of $P_{-}$.

The matrix elements are saved in the order of $M_{X}$, $q^{2}$, $E_{\ell}^{B}$, $P_{+}$ and $P_{-}$. For each variable, the entries is ordered for the first, second and third moments.

Differential signal yields for various signal hypotheses.

Differential signal yields for various signal hypotheses.

Differential signal yields for various signal hypotheses.

Differential signal yields for various signal hypotheses.

Differential signal yields for various signal hypotheses.

Background and data yields in the control and signal region bins. The prediction before ("prefit") and after the background only fit ("b-only") are given separately.

Background and data yields in the control and signal region bins. The prediction before ("prefit") and after the background only fit ("b-only") are given separately.

Background and data yields in the control and signal region bins. The prediction before ("prefit") and after the background only fit ("b-only") are given separately.

Matrix of covariance coefficients between signal region bins. The coefficients are obtained from the background-only fit to the control regions, and serve as input to the simplified likelihood reinterpretation scheme.

Background prediction and observed data yields in the signal region bins. The background yields are obtained from the background-only fit to the corresponding control regions, and serve as input to the simplified likelihood reinterpretation scheme.

Matrix of covariance coefficients between signal region bins. The coefficients are obtained from the background-only fit to the control regions, and serve as input to the simplified likelihood reinterpretation scheme.

Background prediction and observed data yields in the signal region bins. The background yields are obtained from the background-only fit to the corresponding control regions, and serve as input to the simplified likelihood reinterpretation scheme.

Matrix of covariance coefficients between signal region bins. The coefficients are obtained from the background-only fit to the control regions, and serve as input to the simplified likelihood reinterpretation scheme.

Background prediction and observed data yields in the signal region bins. The background yields are obtained from the background-only fit to the corresponding control regions, and serve as input to the simplified likelihood reinterpretation scheme.

Unweighted signal acceptance times efficiency at every cut stage. The requirements called "HCAL mitigation" refer to the requirements imposed in the 2018 data set in order to mitigate the localized failure of the HCAL detector.

Unweighted signal acceptance times efficiency at every cut stage. The requirements called "HCAL mitigation" refer to the requirements imposed in the 2018 data set in order to mitigate the localized failure of the HCAL detector.

Unweighted signal acceptance times efficiency at every cut stage. The requirements called "HCAL mitigation" refer to the requirements imposed in the 2018 data set in order to mitigate the localized failure of the HCAL detector.

Unweighted signal acceptance times efficiency at every cut stage. The requirements called "HCAL mitigation" refer to the requirements imposed in the 2018 data set in order to mitigate the localized failure of the HCAL detector.

Unweighted signal acceptance times efficiency at every cut stage. The requirements called "HCAL mitigation" refer to the requirements imposed in the 2018 data set in order to mitigate the localized failure of the HCAL detector.

Unweighted signal acceptance times efficiency at every cut stage. The requirements called "HCAL mitigation" refer to the requirements imposed in the 2018 data set in order to mitigate the localized failure of the HCAL detector.

Unweighted signal acceptance times efficiency at every cut stage. The requirements called "HCAL mitigation" refer to the requirements imposed in the 2018 data set in order to mitigate the localized failure of the HCAL detector.

Unweighted signal acceptance times efficiency at every cut stage. The requirements called "HCAL mitigation" refer to the requirements imposed in the 2018 data set in order to mitigate the localized failure of the HCAL detector.

Unweighted signal acceptance times efficiency at every cut stage. The requirements called "HCAL mitigation" refer to the requirements imposed in the 2018 data set in order to mitigate the localized failure of the HCAL detector.

Unweighted signal acceptance times efficiency at every cut stage. The requirements called "HCAL mitigation" refer to the requirements imposed in the 2018 data set in order to mitigate the localized failure of the HCAL detector.

Unweighted signal acceptance times efficiency at every cut stage. The requirements called "HCAL mitigation" refer to the requirements imposed in the 2018 data set in order to mitigate the localized failure of the HCAL detector.

Unweighted signal acceptance times efficiency at every cut stage. The requirements called "HCAL mitigation" refer to the requirements imposed in the 2018 data set in order to mitigate the localized failure of the HCAL detector.

Unweighted signal acceptance times efficiency at every cut stage. The requirements called "HCAL mitigation" refer to the requirements imposed in the 2018 data set in order to mitigate the localized failure of the HCAL detector.

Unweighted signal acceptance times efficiency at every cut stage. The requirements called "HCAL mitigation" refer to the requirements imposed in the 2018 data set in order to mitigate the localized failure of the HCAL detector.

Unweighted signal acceptance times efficiency at every cut stage. The requirements called "HCAL mitigation" refer to the requirements imposed in the 2018 data set in order to mitigate the localized failure of the HCAL detector.

Unweighted signal acceptance times efficiency at every cut stage. The requirements called "HCAL mitigation" refer to the requirements imposed in the 2018 data set in order to mitigate the localized failure of the HCAL detector.

Unweighted signal acceptance times efficiency at every cut stage. The requirements called "HCAL mitigation" refer to the requirements imposed in the 2018 data set in order to mitigate the localized failure of the HCAL detector.

Unweighted signal acceptance times efficiency at every cut stage. The requirements called "HCAL mitigation" refer to the requirements imposed in the 2018 data set in order to mitigate the localized failure of the HCAL detector.

Median Expected exclusion contour in the $m_{med}$-$m_{\chi}$ plane in the simplified model with axial couplings.

Observed exclusion contour in the $m_{med}$-$m_{\chi}$ plane in the simplified model with axial couplings.

Expected plus 1 s.d. exclusion contour in the $m_{med}$-$m_{\chi}$ plane in the simplified model with axial couplings.

Expected minus 1 s.d. exclusion contour in the $m_{med}$-$m_{\chi}$ plane in the simplified model with axial couplings.

Median Expected exclusion contour in the $m_{med}$-$m_{\chi}$ plane in the simplified model with vector couplings.

Observed exclusion contour in the $m_{med}$-$m_{\chi}$ plane in the simplified model with vector couplings.

Expected plus 1 s.d. exclusion contour in the $m_{med}$-$m_{\chi}$ plane in the simplified model with vector couplings.

Expected minus 1 s.d. exclusion contour in the $m_{med}$-$m_{\chi}$ plane in the simplified model with vector couplings.

Upper limits on the coupling $g_{\chi}$ in the simplified model with a axial mediator.

Upper limits on the coupling $g_{q}$ in the simplified model with a axial mediator.

Upper limits on the coupling $g_{\chi}$ in the simplified model with a vector mediator.

Upper limits on the coupling $g_{q}$ in the simplified model with a vector mediator.

Exclusion limits on the signal strength in the simplified model with scalar couplings.

Exclusion limits on the signal strength in the simplified model with pseudoscalar couplings.

Exclusion limits on the fundamental Planck scale $M_{D}$ as a function of the number of extra dimensions $d$.

Median Expected exclusion contour in the $m_{med}$-$m_{\chi}$ plane in the fermion portal model.

Observed exclusion contour in the $m_{med}$-$m_{\chi}$ plane in the fermion portal model.

Expected plus 1 s.d. exclusion contour in the $m_{med}$-$m_{\chi}$ plane in the fermion portal model.

Expected minus 1 s.d. exclusion contour in the $m_{med}$-$m_{\chi}$ plane in the fermion portal model.

Tagging efficiency for AK8 jets. The efficiency includes the effect of the machine-learning based DeepAK8 tagger, as well as the application of the mass window requirement on the jet. The efficiency is split depending on the matching generator-level object. For reinterpretation purposes, this efficiency can directly be applied to any AK8 jet that is matched to a given type of generator-level object, with no prior selection on the jet mass. Other acceptance requirements on jet $\p_{T}$ and $\eta$ should still be applied.

Distribution of $m(l,\gamma)$ in the $N_{jet}\geq 3$ signal region.

Distribution of $\Delta R(l,\gamma)$ in the $N_{jet}\geq 3$ signal region.

Distribution of $\Delta R(j,\gamma)$ in the $N_{jet}\geq 3$ signal region.

Fit result of the multijet template obtained with loosely isolated leptons and the electroweak background to the measured $m_{T}(W)$ distribution with isolated leptons in the $N_{jet}=2$, $N_{b jet}=0$ selection for electrons.

Fit result of the multijet template obtained with loosely isolated leptons and the electroweak background to the measured $m_{T}(W)$ distribution with isolated leptons in the $N_{jet}=2$, $N_{b jet}=0$ selection for muons.

Distribution of the invariant mass of the lepton and the photon ($m(l,\gamma)$) in the $N_{jet}\geq 3$, $N_{b jet}=0$ selection for the e channel.

Distribution of the invariant mass of the lepton and the photon ($m(l,\gamma)$) in the $N_{jet}\geq 3$, $N_{b jet}=0$ selection for the $\mu$ channel.

Extracted scale factors for the contribution from misidentified electrons for the three data-taking periods, and the Z$\gamma$, W$\gamma$ simulations.

Predicted and observed yields in the control regions in the $N_{jet}= 3$ and $\geq 4$ seletions using the post-fit values of the nuisance parameters.

Predicted and observed yields in the signal regions in the $N_{jet}= 3$ and $\geq 4$ seletions using the post-fit values of the nuisance parameters.

The measured inclusive ttgamma cross section in the fiducial phase space compared to the prediction from simulation using Madgraph_aMC@NLO at a center-of-mass energy of 13 TeV.

Summary of the measured cross section ratios with respect to the NLO cross section prediction for signal regions binned in the electron channel, muon channel and the combined single lepton measurement.

The unfolded differential cross sections for $p_{T}(\gamma)$ and the comparison to simulations.

The unfolded differential cross sections for $|\eta(\gamma)|$ and the comparison to simulations.

The unfolded differential cross sections for $\Delta R(l,\gamma)$ and the comparison to simulations.

The covariance matrix of systematic uncertainties for the unfolded differential measurement for $p_{T}(\gamma)$.

The covariance matrix of systematic uncertainties for the unfolded differential measurement for $|\eta(\gamma)|$.

The covariance matrix of systematic uncertainties for the unfolded differential measurement for $\Delta R(l,\gamma)$.

The covariance matrix of statistic uncertainties for the unfolded differential measurement for $p_{T}(\gamma)$.

The covariance matrix of statistic uncertainties for the unfolded differential measurement for $|\eta(\gamma)|$.

The covariance matrix of statistic uncertainties for the unfolded differential measurement for $\Delta R(l,\gamma)$.

The correlation matrix of statistical uncertainties for the unfolded differential measurement for $p_{T}(\gamma)$.

The correlation matrix of statistical uncertainties for the unfolded differential measurement for $|\eta(\gamma)|$.

The correlation matrix of statistical uncertainties for the unfolded differential measurement for $\Delta R(l,\gamma)$.

The correlation matrix of systematic uncertainties for the unfolded differential measurement for $p_{T}(\gamma)$.

The correlation matrix of systematic uncertainties for the unfolded differential measurement for $|\eta(\gamma)|$.

The correlation matrix of systematic uncertainties for the unfolded differential measurement for $\Delta R(l,\gamma)$.

Summary of the one-dimensional intervals at 68 and 95% CL.

The observed and predicted post-fit yields for the combined Run 2 data set in the SR3 signal region for the electron channel.

The observed and predicted post-fit yields for the combined Run 2 data set in the SR3 signal region for the muon channel.

The observed and predicted post-fit yields for the combined Run 2 data set in the SR4p signal region for the electron channel.

The observed and predicted post-fit yields for the combined Run 2 data set in the SR4p signal region for the muon channel.

Negative log-likelihood ratio values with respect to the best fit value of the one-dimensional profiled scan for the Wilson coefficient $c_{tZ}$.

Negative log-likelihood ratio values with respect to the best fit value of the one-dimensional profiled scan for the Wilson coefficient $c^{I}_{tZ}$.

Negative log-likelihood ratio values with respect to the best fit value of the one-dimensional scan for the Wilson coefficient $c_{tZ}$.

Negative log-likelihood ratio values with respect to the best fit value of the one-dimensional scan for the Wilson coefficient $c^{I}_{tZ}$.

Negative log-likelihood ratio values with respect to the best fit value of the two-dimensional scan for the Wilson coefficients $c_{tZ}$ and $c^{I}_{tZ}$.