Interpretation of results in the T1tttt model. The colored regions show the upper limits (95\% CL) on the production cross section for $pp\rightarrow \tilde{g}\tilde{g},\tilde{g}\rightarrow t\bar{t}\tilde{\chi}^0_1$ in the $m_{\tilde{g}}$-$m_{\tilde{\chi}^0_1}$ plane.

Interpretation of results in the T1tttt model. The colored regions show the upper limits (95\% CL) on the production cross section for $pp\rightarrow \tilde{g}\tilde{g},\tilde{g}\rightarrow t\bar{t}\tilde{\chi}^0_1$ in the $m_{\tilde{g}}$-$m_{\tilde{\chi}^0_1}$ plane. The curve shows the observed limit on the corresponding SUSY particle masses obtained by comparing the excluded cross section with theoretical cross sections.

Interpretation of results in the T1tttt model. The colored regions show the upper limits (95\% CL) on the production cross section for $pp\rightarrow \tilde{g}\tilde{g},\tilde{g}\rightarrow t\bar{t}\tilde{\chi}^0_1$ in the $m_{\tilde{g}}$-$m_{\tilde{\chi}^0_1}$ plane. The curve shows the observed limit on the corresponding SUSY particle masses obtained by comparing the excluded cross section with $+1\sigma$ theoretical cross sections.

Interpretation of results in the T1tttt model. The colored regions show the upper limits (95\% CL) on the production cross section for $pp\rightarrow \tilde{g}\tilde{g},\tilde{g}\rightarrow t\bar{t}\tilde{\chi}^0_1$ in the $m_{\tilde{g}}$-$m_{\tilde{\chi}^0_1}$ plane. The curve shows the observed limit on the corresponding SUSY particle masses obtained by comparing the excluded cross section with $-1\sigma$ theoretical cross sections.

Interpretation of results in the T1tttt model. The colored regions show the upper limits (95\% CL) on the production cross section for $pp\rightarrow \tilde{g}\tilde{g},\tilde{g}\rightarrow t\bar{t}\tilde{\chi}^0_1$ in the $m_{\tilde{g}}$-$m_{\tilde{\chi}^0_1}$ plane. The curve shows the expected limit on the corresponding SUSY particle masses obtained by comparing the excluded cross section with theoretical cross sections.

Interpretation of results in the T1tttt model. The colored regions show the upper limits (95\% CL) on the production cross section for $pp\rightarrow \tilde{g}\tilde{g},\tilde{g}\rightarrow t\bar{t}\tilde{\chi}^0_1$ in the $m_{\tilde{g}}$-$m_{\tilde{\chi}^0_1}$ plane. The curve shows the expected limit on the corresponding SUSY particle masses corresponding to a $+1\sigma$ variation of the experimental uncertainty.

Interpretation of results in the T1tttt model. The colored regions show the upper limits (95\% CL) on the production cross section for $pp\rightarrow \tilde{g}\tilde{g},\tilde{g}\rightarrow t\bar{t}\tilde{\chi}^0_1$ in the $m_{\tilde{g}}$-$m_{\tilde{\chi}^0_1}$ plane. The curve shows the expected limit on the corresponding SUSY particle masses corresponding to a $-1\sigma$ variation of the experimental uncertainty.

Observed excluded region (95\% CL), shown in blue, in the $m_{\tilde{g}}$-$m_{\tilde{\chi}^0_1}$ plane for a model combining T5tttt, gluino pair production, followed by gluino decay to an on-shell top squark, together with a model for direct top squark pair production. The top squarks decay via the two-body process $\tilde{t}\rightarrow t\tilde{\chi}^0_1$. The neutralino and top squark masses are related by the constraint $m_{\tilde{t}} = m_{\tilde{\chi}^0_1} + 175$ GeV.

Expected excluded region (95\% CL), shown in blue, in the $m_{\tilde{g}}$-$m_{\tilde{\chi}^0_1}$ plane for a model combining T5tttt, gluino pair production, followed by gluino decay to an on-shell top squark, together with a model for direct top squark pair production. The top squarks decay via the two-body process $\tilde{t}\rightarrow t\tilde{\chi}^0_1$. The neutralino and top squark masses are related by the constraint $m_{\tilde{t}} = m_{\tilde{\chi}^0_1} + 175$ GeV.

Figure 9a (upper left). The 95% CL upper limit on the production cross section of the T1tttt simplified models as a function of the gluino and LSP masses.

Figure 9a (upper left). Observed exclusion region at 95% CL assuming 100% branching fraction.

Figure 9a (upper left). Expected exclusion region at 95% CL assuming 100% branching fraction.

Figure 9b (upper right). The 95% CL upper limit on the production cross section of the T1ttbb simplified models as a function of the gluino and LSP masses.

Figure 9b (upper right). Observed exclusion region at 95% CL.

Figure 9b (upper right). Expected exclusion region at 95%.

Figure 9c (bottom left). The 95% CL upper limit on the production cross section of the T5tttt simplified models as a function of the gluino and LSP masses. imits are not given for the T5tttt model for mχ~01< 50 GeV for the reason stated in the text.

Figure 9c (bottom left). Observed exclusion region at 95% CL assuming 100% branching fraction.

Figure 9c (bottom left). Expected exclusion region at 95% CL assuming 100% branching fraction.

Figure 9d (bottom right). The 95% CL upper limit on the production cross section of the T5ttcc simplified models as a function of the gluino and LSP masses.

Figure 9d (bottom right). Observed exclusion region at 95% CL.

Figure 9d (bottom right). Expected exclusion region at 95% CL.

The observed number of events and the total background prediction for search regions with Nt = 1 and Nb = 1.

The observed number of events and the total background prediction for search regions with Nt = 1 and Nb ≥ 2.

The observed number of events and the total background prediction for search regions with Nt ≥ 2.

Selected T2tt signal yields for the aggregate search bins.

Selected gluino mediated signal yields for the aggregate search bins.

Figure 2. -1 sigma excluded gluino and neutralino masses at 95% CL for the T1tttt.

Figure 2. Expected exclusion for the T1tttt model.

Figure 2. +1 sigma expected exclusion for the T1tttt model.

Figure 2. -1 sigma expected exclusion for the T1tttt model.

Figure 2. Excluded gluino and neutralino masses at 95% CL for the T5tttt.

Figure 2. Expected exclusion for the T5tttt model.

Table 1. Observed event yields and mean background yields from the predictive fit in the 18 bins of the signal region R4.

Table 2. Observed event yields and mean background yields from the predictive fit in four aggregate search bins.

Observed exclusion contour at 95% CL for the T6gg model.

Expected exclusion contour at 95% CL for the T6gg model.

Exclusion limits on the SUSY cross section at 95% CL for the T6Wg model.

Observed exclusion contour at 95% CL for the T6Wg model.

Expected exclusion contour at 95% CL for the T6Wg model.

Exclusion limits on the SUSY cross section at 95% CL for the T5gg model.

Observed exclusion contour at 95% CL for the T5gg model.

Expected exclusion contour at 95% CL for the T5gg model.

Exclusion limits on the SUSY cross section at 95% CL for the T5Wg model.

Observed exclusion contour at 95% CL for the T5Wg model.

Expected exclusion contour at 95% CL for the T5Wg model.

Distribution of $m_{jj}$ in data and simulation for events in the signal regions of the low-mass search. The points show the data while the filled histograms show the background contributions. The gray band shows the statistical and systematic uncertainties in the background, while the dashed vertical region ("Higgs") shows the expected SM Higgs boson mass range, which is excluded from this analysis. A 600 GeV bulk graviton signal prediction is represented by the black dashed histogram; for visibility, the signal cross-section is increased by a factor of 50. The background normalization is derived from the final fit to the $m_{VZ}$ observable in data.

The m$_{j}$ distributions of the events in data, compared to the expected background shape, for the high-mass analysis in the electron channel, for the high-purity category. The expected background shape is extracted from a fit to the data sidebands (Z+jets) or derived from simulation ("top quark" and "VV"). The dashed region around the background sum represents the uncertainty in the Z+jets distribution, while the dashed vertical region ("Higgs") shows the expected SM Higgs boson mass range, excluded from the analysis.

The m$_{j}$ distributions of the events in data, compared to the expected background shape, for the high-mass analysis in the electron channel, for the low-purity category. The expected background shape is extracted from a fit to the data sidebands (Z+jets) or derived from simulation ("top quark" and "VV"). The dashed region around the background sum represents the uncertainty in the Z+jets distribution, while the dashed vertical region ("Higgs") shows the expected SM Higgs boson mass range, excluded from the analysis.

The m$_{j}$ distributions of the events in data, compared to the expected background shape, for the high-mass analysis in the muon channel, for the high-purity category. The expected background shape is extracted from a fit to the data sidebands (Z+jets) or derived from simulation ("top quark" and "VV"). The dashed region around the background sum represents the uncertainty in the Z+jets distribution, while the dashed vertical region ("Higgs") shows the expected SM Higgs boson mass range, excluded from the analysis.

The m$_{j}$ distributions of the events in data, compared to the expected background shape, for the high-mass analysis in the muon channel, for the low-purity category. The expected background shape is extracted from a fit to the data sidebands (Z+jets) or derived from simulation ("top quark" and "VV"). The dashed region around the background sum represents the uncertainty in the Z+jets distribution, while the dashed vertical region ("Higgs") shows the expected SM Higgs boson mass range, excluded from the analysis.

Expected and observed distributions of the resonance candidate mass $m_{\text{VZ}}$ in the high-mass analysis, in the electron channel, high-purity category. The shaded area represents the post-fit uncertainty in the background. The expected contribution from W' signal candidates with mass $m_{\text{X}} = 2000$ GeV, normalized to a cross section of 100 fb$^{-1}$, is also shown.

Expected and observed distributions of the resonance candidate mass $m_{\text{VZ}}$ in the high-mass analysis, in the electron channel, low-purity category. The shaded area represents the post-fit uncertainty in the background. The expected contribution from W' signal candidates with mass $m_{\text{X}} = 2000$ GeV, normalized to a cross section of 100 fb$^{-1}$, is also shown.

Expected and observed distributions of the resonance candidate mass $m_{\text{VZ}}$ in the high-mass analysis, in the muon channel, high-purity category. The shaded area represents the post-fit uncertainty in the background. The expected contribution from W' signal candidates with mass $m_{\text{X}} = 2000$ GeV, normalized to a cross section of 100 fb$^{-1}$, is also shown.

Expected and observed distributions of the resonance candidate mass $m_{\text{VZ}}$ in the high-mass analysis, in the muon channel, low-purity category. The shaded area represents the post-fit uncertainty in the background. The expected contribution from W' signal candidates with mass $m_{\text{X}} = 2000$ GeV, normalized to a cross section of 100 fb$^{-1}$, is also shown.

Sideband $m_{ZV}$ distribution for the low-mass search in the merged V, untagged category, after fitting the sideband data alone. The points show the data while the filled histograms show the background contributions. Electron and muon categories are combined. The gray band indicates the statistical and post-fit systematic uncertainties in the normalization and shape of the background. Larger bin widths are used at higher values of $m_{ZV}$; the bin widths are indicated by the horizontal error bars.

Sideband $m_{ZV}$ distribution for the low-mass search in the merged V, tagged category, after fitting the sideband data alone. The points show the data while the filled histograms show the background contributions. Electron and muon categories are combined. The gray band indicates the statistical and post-fit systematic uncertainties in the normalization and shape of the background. Larger bin widths are used at higher values of $m_{ZV}$; the bin widths are indicated by the horizontal error bars.

Sideband $m_{ZV}$ distribution for the low-mass search in the resolved V, untagged category, after fitting the sideband data alone. The points show the data while the filled histograms show the background contributions. Electron and muon categories are combined. The gray band indicates the statistical and post-fit systematic uncertainties in the normalization and shape of the background. Larger bin widths are used at higher values of $m_{ZV}$; the bin widths are indicated by the horizontal error bars.

Sideband $m_{ZV}$ distribution for the low-mass search in the resolved V, tagged category, after fitting the sideband data alone. The points show the data while the filled histograms show the background contributions. Electron and muon categories are combined. The gray band indicates the statistical and post-fit systematic uncertainties in the normalization and shape of the background. Larger bin widths are used at higher values of $m_{ZV}$; the bin widths are indicated by the horizontal error bars.

The signal region $m_{VZ}$ distribution for the low-mass search, in the the merged V, untagged category, after fitting the signal and sideband regions. Electron and muon categories are combined. A 600 GeV bulk graviton signal prediction is represented by the black dashed histogram. The gray band indicates the statistical and post-fit systematic uncertainties in the normalization and shape of the background. Larger bin widths are used at higher values of $m_{VZ}$; the bin widths are indicated by the horizontal error bars.

The signal region $m_{VZ}$ distribution for the low-mass search, in the the merged V, tagged category, after fitting the signal and sideband regions. Electron and muon categories are combined. A 600 GeV bulk graviton signal prediction is represented by the black dashed histogram. The gray band indicates the statistical and post-fit systematic uncertainties in the normalization and shape of the background. Larger bin widths are used at higher values of $m_{VZ}$; the bin widths are indicated by the horizontal error bars.

The signal region $m_{VZ}$ distribution for the low-mass search, in the the resolved V, untagged category, after fitting the signal and sideband regions. Electron and muon categories are combined. A 600 GeV bulk graviton signal prediction is represented by the black dashed histogram. The gray band indicates the statistical and post-fit systematic uncertainties in the normalization and shape of the background. Larger bin widths are used at higher values of $m_{VZ}$; the bin widths are indicated by the horizontal error bars.

The signal region $m_{VZ}$ distribution for the low-mass search, in the the resolved V, tagged category, after fitting the signal and sideband regions. Electron and muon categories are combined. A 600 GeV bulk graviton signal prediction is represented by the black dashed histogram. The gray band indicates the statistical and post-fit systematic uncertainties in the normalization and shape of the background. Larger bin widths are used at higher values of $m_{VZ}$; the bin widths are indicated by the horizontal error bars.

Observed and expected 95% CL upper limit on $\sigma_{W'} Br(W' \to ZW)$ as a function of the resonance mass, taking into account all statistical and systematic uncertainties. The electron and muon channels and the various categories used in the analysis are combined together. The green (inner) and yellow (outer) bands represent the 68% and 95% coverage of the expected limit in the background-only hypothesis.

Observed and expected 95% CL upper limit on $\sigma_{G} Br(G \to ZZ)$ as a function of the resonance mass, taking into account all statistical and systematic uncertainties. The electron and muon channels and the various categories used in the analysis are combined together. The green (inner) and yellow (outer) bands represent the 68% and 95% coverage of the expected limit in the background-only hypothesis.

Observed local p-values for W' narrow resonances as a function of the resonance mass.

Observed local p-values for G narrow resonances as a function of the resonance mass.

Figure 8. Gluino pair production cross section.

Table 1. Post-fit yields for the background-only fit, observed data, and expected yields for mgluino = 1600 GeV in each search bin.

Number of predicted signal events in the $\mu\mu X$ final states, for several benchmark signal hypotheses for Majorana HNL. The quoted uncertainties include statistical and systematic uncertainties.

The 95% CL limits on $|V_{Ne}|^2$ as a function of the HNL mass for a Majorana HNL. Values of $-1$ indicate that no limit is available for the mass point.

The 95% CL limits on $|V_{N\mu}|^2$ as a function of the HNL mass for a Majorana HNL. Values of $-1$ indicate that no limit is available for the mass point.

The 95% CL limits on $|V_{Ne}|^2$ as a function of the HNL mass for a Dirac HNL. Values of $-1$ indicate that no limit is available for the mass point.

The 95% CL limits on $|V_{N\mu}|^2$ as a function of the HNL mass for a Dirac HNL. Values of $-1$ indicate that no limit is available for the mass point.

Covariance matrix for the background prediction in the signal region bins for the $eeX$ final states. The bin labels specify the channel, the $m(\ell_2,\ell_3)$ bin, and the $\Delta_{\mathrm{2D}}$ bin.

Covariance matrix for the background prediction in the signal region bins for the $\mu\mu X$ final states. The bin labels specify the channel, the $m(\ell_2,\ell_3)$ bin, and the $\Delta_{\mathrm{2D}}$ bin.

Correlation matrix for the background prediction the signal region bins for the $eeX$ final states. The bin labels specify the channel, the $m(\ell_2,\ell_3)$ bin, and the $\Delta_{\mathrm{2D}}$ bin.

Correlation matrix for the background prediction the signal region bins for the $\mu\mu X$ final states. The bin labels specify the channel, the $m(\ell_2,\ell_3)$ bin, and the $\Delta_{\mathrm{2D}}$ bin.

Mean lifetime and production cross section for HNLs of Majorana nature, assuming the decay to three charged leptons and coupling to only electrons or muons, as evaluated from simulated event samples.

Distribution of $m_{\mathrm{jjj}}$ for preselected events with $\mathrm{N}_{j}$ = 3

Distribution of $m_{\mathrm{j}}$ for preselected events with $\mathrm{N}_{j}$ = 3

Distribution of the deep-WH value of the highest-mass jet with 60 < $m_{\mathrm{j}}$ < 100 GeV for preselected events with $\mathrm{N}_{j}$ = 3

scale factors (SFs) for W, $t^{2}$, and q/g matched jets in the low-$m_{\mathrm{j}}$ and low-$p_{\mathrm{T}}$ (LL) bin, as functions of the deep-W discriminant value.

scale factors (SFs) for W, $t^{2}$, and q/g matched jets in the low-$m_{\mathrm{j}}$ and high-$p_{\mathrm{T}}$ (LH) bin, as functions of the deep-W discriminant value.

scale factors (SFs) for $t^{2}$, $t^{3,4}$, and q/g matched jets in the high-$m_{\mathrm{j}}$ and low-$p_{\mathrm{T}}$ (HL) bin, as functions of the deep-WH discriminant value.

scale factors (SFs) for $t^{2}$, $t^{3,4}$, and q/g matched jets in the high-$m_{\mathrm{j}}$ and high-$p_{\mathrm{T}}$ (HH) bin, as functions of the deep-WH discriminant value.

The deep-W discriminant of the jet with highest mass in the single-lepton sideband for LL samples.

The deep-W discriminant of the jet with highest mass in the single-lepton sideband for LH samples.

The deep-WH discriminant of the jet with highest mass in the single-lepton sideband for HL samples.

The deep-WH discriminant of the jet with highest mass in the single-lepton sideband for HH samples.

Comparison of the distribution of data and simulated backgrounds, as a function of the deep-W discriminant value for the highest-mass jet in CR1, after the scale factors have been applied.

Comparison of the distribution of data and simulated backgrounds, as a function of the deep-WH discriminant value for the highest-mass jet in CR2, after the scale factors have been applied.

Comparison of the distribution of data and simulated backgrounds, as a function of the deep-WH discriminant value for the highest-mass jet in CR3, after the scale factors have been applied.

Comparison of the distribution of data and simulated backgrounds, as a function of the deep-W discriminant value for the highest-mass jet in CR45, after the scale factors have been applied.

Comparison of the distribution of data and simulated backgrounds, as a function of the deep-W discriminant value for the highest-mass jet in CR6, after the scale factors have been applied.

The $m_{\mathrm{j}^{\mathrm{max}}}$ distributions for different jet types for SR1–3 events of the signal with $m_{\mathrm{W}_{\mathrm{KK}}}$ = 2.5 TeV, $m_{\mathrm{R}}$ = 0.2 TeV$ without deep-W (WH) constraints.

The deep-W distribution normalized to unity for the shown components of of the signal with $m_{\mathrm{W}_{\mathrm{KK}}}$ = 2.5 TeV, $m_{\mathrm{R}}$ = 0.2 TeV. The $t^{3,4}$ jets from the preselected sample, normalized to unity, are superimposed to compare shapes with the $R^{3q}$ and $R^{4q}$ distributions.

The deep-WH distribution normalized to unity for the shown components of of the signal with $m_{\mathrm{W}_{\mathrm{KK}}}$ = 2.5 TeV, $m_{\mathrm{R}}$ = 0.2 TeV. The $t^{3,4}$ jets from the preselected sample, normalized to unity, are superimposed to compare shapes with the $R^{3q}$ and $R^{4q}$ distributions.

The $m_{\mathrm{jj}}$ distribution for CR1 for data and simulation.

The $m_{\mathrm{jj}}$ distribution for CR2 for data and simulation.

The $m_{\mathrm{jj}}$ distribution for CR3 for data and simulation.

The $m_{\mathrm{jjj}}$ distribution for CR45 for data and simulation.

The $m_{\mathrm{jjj}}$ distribution for CR6 for data and simulation.

Post-fit distributions of the reconstructed triboson system ($m_{\mathrm{jj}}$) in data and simulation for SR1.

Post-fit distributions of the reconstructed triboson system ($m_{\mathrm{jj}}$) in data and simulation for SR2.

Post-fit distributions of the reconstructed triboson system ($m_{\mathrm{jj}}$) in data and simulation for SR3.

Post-fit distributions of the reconstructed triboson system ($m_{\mathrm{jjj}}$) in data and simulation for SR4.

Post-fit distributions of the reconstructed triboson system ($m_{\mathrm{jjj}}$) in data and simulation for SR5.

Post-fit distributions of the reconstructed triboson system ($m_{\mathrm{jjj}}$) in data and simulation for SR6.

Observed upper limits at $95\%$ CL on the signal cross section $\times$ branching fraction as functions of the $m_{\mathrm{W}_{\mathrm{KK}}}$ and $m_{\mathrm{R}}$ resonance masses.

Expected median lower limit contour on the $m_{\mathrm{W}_{\mathrm{KK}}}$ and $m_{\mathrm{R}}$ plane.

Expected $+ 1$ s.d. lower limit contour on the $m_{\mathrm{W}_{\mathrm{KK}}}$ and $m_{\mathrm{R}}$ plane.

Expected - 1 s.d. lower limit contour on the $m_{\mathrm{W}_{\mathrm{KK}}}$ and $m_{\mathrm{R}}$ plane.

Observed lower limit contour on the $m_{\mathrm{W}_{\mathrm{KK}}}$ and $m_{\mathrm{R}}$ plane.

Distribution of $m_{T}(W)$ in the $N_{jet}\geq 3$ signal region.

Distribution of $M_{3}$ in the $N_{jet}\geq 3$ signal region.

Distribution of $M_{3}$ in the $N_{jet}\geq 3$ signal region.

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

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(l,\gamma)$ in the $N_{jet}\geq 3$ signal region.

Distribution of $\Delta R(j,\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 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.

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

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.

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

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.

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.

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 $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 $|\eta(\gamma)|$ and the comparison to simulations.

The unfolded differential cross sections for $\Delta R(l,\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 $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 $|\eta(\gamma)|$.

The covariance matrix of systematic uncertainties for the unfolded differential measurement for $\Delta R(l,\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 $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 $|\eta(\gamma)|$.

The covariance matrix of statistic uncertainties for the unfolded differential measurement for $\Delta R(l,\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 $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 $|\eta(\gamma)|$.

The correlation matrix of statistical uncertainties for the unfolded differential measurement for $\Delta R(l,\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 $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 $|\eta(\gamma)|$.

The correlation matrix of systematic uncertainties for the unfolded differential measurement for $\Delta R(l,\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.

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

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_{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 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_{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 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}$.

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

SR1 diboson invariant mass distribution in the resolved tagged category after fitting the signal and sideband regions using a signal (ALP linear ZZ) plus background model. The last bin includes events with diboson invariant masses up to 3000 GeV. The signal is normalized to the 95% CL cross section limit at $f_a$ = 3 TeV (the scale factor used in the original figure for better visibility is not applied here).

SR2 diboson invariant mass distribution in the boosted untagged category after fitting the signal and sideband regions using a signal (ALP chiral ZH) plus background model. The last bin includes events with diboson invariant masses up to 3000 GeV. The signal is normalized to the 95% CL cross section limit at $f_a$ = 3 TeV (the scale factor used in the original figure for better visibility is not applied here).

SR2 diboson invariant mass distribution in the boosted tagged category after fitting the signal and sideband regions using a signal (ALP chiral ZH) plus background model. The last bin includes events with diboson invariant masses up to 3000 GeV. The signal is normalized to the 95% CL cross section limit at $f_a$ = 3 TeV.

SR2 diboson invariant mass distribution in the resolved untagged category after fitting the signal and sideband regions using a signal (ALP chiral ZH) plus background model. The last bin includes events with diboson invariant masses up to 3000 GeV. The signal is normalized to the 95% CL cross section limit at $f_a$ = 3 TeV (the scale factor used in the original figure for better visibility is not applied here).

SR2 diboson invariant mass distribution in the resolved tagged category after fitting the signal and sideband regions using a signal (ALP chiral ZH) plus background model. The last bin includes events with diboson invariant masses up to 3000 GeV. The signal is normalized to the 95% CL cross section limit at $f_a$ = 3 TeV (the scale factor used in the original figure for better visibility is not applied here).

Observed and expected 95% CL upper limit on $\sigma \mathcal{B}$(G -> ZZ) as a function of the resonance mass, taking into account all statistical and systematic uncertainties. The electron and muon channels and the various categories used in the analysis are combined. The inner green band and the outer yellow band represent the 68 and 95% coverage of the expected limits in the background-only hypothesis.

Observed and expected 95% CL upper limit on $\sigma \mathcal{B}$(W' -> ZW) as a function of the resonance mass, taking into account all statistical and systematic uncertainties. The electron and muon channels and the various categories used in the analysis are combined. The inner green band and the outer yellow band represent the 68 and 95% coverage of the expected limits in the background-only hypothesis.

Observed and expected 95% CL upper limit on the product of the ALP coupling to gluons and the coupling to Z bosons as a function of the mass scale $f_a$ for ALP masses smaller than 100 GeV.

Observed and expected 95% CL upper limit on the product of the ALP coupling to gluons and the chiral $a_\widetilde{2D}$ coupling as a function of the mass scale $f_a$ for ALP masses smaller than 100 GeV.