Distribution of an input to the BDT classifier in the $3\ell$ category: The linear discrimminant $p_{T}^{miss,LD}$. The discriminant is defined by the relation $p_{T}^{miss,LD} = 0.6 p_{T}^{miss} + 0.4 H_{T}^{miss}$, where $H_{T}^{miss}$ corresponds to the magnitude of the vector $p_{T}$ sum of all electrons, muons, taus and small radius jets passing the selection critertia, and $p_{T}^{miss}$ represents the missing transverse momentum vector computed as the negative vector $p_{T}$ sum of all the particles reconstructed by the PF algorithm in the event.

Distributions of the transverse mass $m_{T}$ of the lepton that is not originating from the Z boson decay and the missing transverse momentum in the $3\ell$/WZ control region. The distributions expected for the WZ and ZZ as well as for other background processes are shown for the values of nuisance parameters obtained from the ML fit used in the signal extraction.

Distributions of the four lepton invariant mass $m_{4\ell}$ in the $4\ell$/ZZ control region. The distributions expected for the ZZ as well as for other background processes are shown for the values of nuisance parameters obtained from the ML fit used in the signal extraction.

Distributions of the transverse mass $m_{T}$ of the leading lepton and the missing transverse momentum in the $2\ell$(ss) control region. The distributions expected for the misidentified $\ell$ background as well as for other background processes is shown for the values of nuisance parameters obtained from the background-only ML fit, in which the HH signal is constrained to be zero.

Distributions of the reconstructed HH mass $m_{HH}$ in the $2\ell+2 au_{h}$ control region. The distributions expected for the misidentified $\ell/ au_{h}$ background as well as for other background processes is shown for the values of nuisance parameters obtained from the background-only ML fit, in which the HH signal is constrained to be zero.

Distribution in the output of the BDT trained for nonresonant HH production and evaluated for the benchmark scenario JHEP04 BM7 for the $2\ell$(ss) category. The SM HH signal is shown for a cross section amounting to $30$ times the value predicted in the SM. The distributions expected for the background processes are shown for the values of nuisance parameters obtained from the ML fit of the signal+background hypothesis to the data. The binning is chosen to approximate signal quantiles.

Distribution in the output of the BDT trained for nonresonant HH production and evaluated for the benchmark scenario JHEP04 BM7 for the $3\ell$ category. The SM HH signal is shown for a cross section amounting to $30$ times the value predicted in the SM. The distributions expected for the background processes are shown for the values of nuisance parameters obtained from the ML fit of the signal+background hypothesis to the data. The binning is chosen to approximate signal quantiles.

Distribution in the output of the BDT trained for nonresonant HH production and evaluated for the benchmark scenario JHEP04 BM7 for the $4\ell$ category. The SM HH signal is shown for a cross section amounting to $30$ times the value predicted in the SM. The distributions expected for the background processes are shown for the values of nuisance parameters obtained from the ML fit of the signal+background hypothesis to the data. The binning is chosen to approximate signal quantiles.

Distribution in the output of the BDT trained for nonresonant HH production and evaluated for the benchmark scenario JHEP04 BM7 for the $3\ell+1 au_{h}$ category. The SM HH signal is shown for a cross section amounting to $30$ times the value predicted in the SM. The distributions expected for the background processes are shown for the values of nuisance parameters obtained from the ML fit of the signal+background hypothesis to the data. The binning is chosen to approximate background quantiles.

Distribution in the output of the BDT trained for nonresonant HH production and evaluated for the benchmark scenario JHEP04 BM7 for the $2\ell+2 au_{h}$ category. The SM HH signal is shown for a cross section amounting to $30$ times the value predicted in the SM. The distributions expected for the background processes are shown for the values of nuisance parameters obtained from the ML fit of the signal+background hypothesis to the data. The binning is chosen to approximate background quantiles.

Distribution in the output of the BDT trained for nonresonant HH production and evaluated for the benchmark scenario JHEP04 BM7 for the $1\ell+3 au_{h}$ category. The SM HH signal is shown for a cross section amounting to $30$ times the value predicted in the SM. The distributions expected for the background processes are shown for the values of nuisance parameters obtained from the ML fit of the signal+background hypothesis to the data. The binning is chosen to approximate background quantiles.

Distribution in the output of the BDT trained for nonresonant HH production and evaluated for the benchmark scenario JHEP04 BM7 for the $4 au_{h}$ category. The SM HH signal is shown for a cross section amounting to $30$ times the value predicted in the SM. The distributions expected for the background processes are shown for the values of nuisance parameters obtained from the ML fit of the signal+background hypothesis to the data. The binning is chosen to approximate background quantiles.

Observed and expected $95\%$ CL upper limits on the SM HH production cross section, obtained for both individual search categories and from a simultaneous fit of all seven categories combined.

Observed and expected $95\%$ CL upper limits on the HH production cross section as a function of the Higgs boson self-coupling strength modifier $\kappa_\lambda$ for the combination of all seven categories. All Higgs boson couplings other than $\lambda$ are assumed to have the values predicted in the SM.

Observed and expected $95\%$ CL upper limits on the HH production cross section as a function of the Higgs boson self-coupling strength modifier $\kappa_\lambda$ for the seven different categories as well as their combination. All Higgs boson couplings other than $\lambda$ are assumed to have the values predicted in the SM.

Observed and expected $95\%$ CL upper limits on the HH production cross section for the twelve benchmark scenarios from doi:10.1007/JHEP04(2016)126, the additional benchmark scenario 8a from doi:10.1007/JHEP09(2018)057, the seven benchmark scenarios from doi:10.1007/JHEP03(2020)091, and for the SM.

Observed and expected $95\%$ CL upper limits on the HH production cross section for the twelve benchmark scenarios from doi:10.1007/JHEP04(2016)126, the additional benchmark scenario 8a from doi:10.1007/JHEP09(2018)057, the seven benchmark scenarios from doi:10.1007/JHEP03(2020)091, and for the SM. The limits are shown for all seven search categories as well as for their combination.

Observed and expected $95\%$ CL upper limits on the HH production cross section as a function of the effective coupling $c_{2}$ for the combination of all seven categories. All Higgs boson couplings other than $c_{2}$ are assumed to have the values predicted in the SM.

Observed and expected $95\%$ CL upper limits on the HH production cross section as a function of the effective coupling $c_{2}$ and the top Yukawa coupling modifier $\kappa_{t}$ for the combination of all seven categories. All Higgs boson couplings other than $c_{2}$ and $\kappa_{t}$ are assumed to have the values predicted in the SM. The position of the SM in the ($c_{2}-\kappa_{t}$) plane, as well as the best fit value of $(c_{2},\kappa_{t})=(1.05, 1.74)$ together with contours for the theory cross section are shown as well.

Observed and expected $95\%$ CL upper limits on the HH production cross section as a function of the Higgs boson self-coupling modifier $\kappa_{\lambda}$ and the top Yukawa coupling modifier $\kappa_{t}$ for the combination of all seven categories. All Higgs boson couplings other than $\kappa_{\lambda}$ and $\kappa_{t}$ are assumed to have the values predicted in the SM. The position of the SM in the ($\kappa_{\lambda}-\kappa_{t}$) plane, as well as the best fit value of $(\kappa_{\lambda},\kappa_{t})=(-3.54, -1.70)$ together with contours for the theory cross section are shown as well.

Observed and expected $95\%$ CL upper limits on the HH production cross section as a function of the effective coupling $c_{2}$ and the Higgs boson self-coupling modifier $\kappa_{\lambda}$ for the combination of all seven categories. All Higgs boson couplings other than $c_{2}$ and $\kappa_{\lambda}$ are assumed to have the values predicted in the SM. The position of the SM in the ($c_{2}-\kappa_{\lambda}$) plane, as well as the best fit value of $(c_{2},\kappa_{\lambda})=(-0.66, 5.33)$ together with contours for the theory cross section are shown as well.

Observed and expected $95\%$ CL upper limits on the production of new particles X of spin $0$ and mass $m_{X}$ in the range $250$-$1000$ GeV, which decay to Higgs boson pairs.

Observed and expected $95\%$ CL upper limits on the production of new particles X of spin $0$ and mass $m_{X}$ in the range $250$-$1000$ GeV, which decay to Higgs boson pairs. The Limit is shown for the seven different search categories as well as their combination.

Observed and expected $95\%$ CL upper limits on the production of new particles X of spin $2$ and mass $m_{X}$ in the range $250$-$1000$ GeV, which decay to Higgs boson pairs.

Observed and expected $95\%$ CL upper limits on the production of new particles X of spin $2$ and mass $m_{X}$ in the range $250$-$1000$ GeV, which decay to Higgs boson pairs. The Limit is shown for the seven different search categories as well as their combination.

Distribution in the output of the BDT trained for resonances of spin 2 and mass $750\,$GeV production for the $2\ell$(ss) category. The resonant HH signal is shown for a cross section amounting to $1\,$pb. The distributions expected for the background processes are shown for the values of nuisance parameters obtained from the ML fit of the signal+background hypothesis to the data. The binning is chosen to approximate signal quantiles.

Distribution in the output of the BDT trained for resonances of spin 2 and mass $750\,$GeV production for the $3\ell$ category. The resonant HH signal is shown for a cross section amounting to $1\,$pb. The distributions expected for the background processes are shown for the values of nuisance parameters obtained from the ML fit of the signal+background hypothesis to the data. The binning is chosen to approximate signal quantiles.

No description provided.

Y defined as (qp)/(kp), where k and p are four-momenta of the incident electron and proton, respectively, q is four-momentum transfer.

Y defined as (qp)/(kp), where k and p are four-momenta of the incident electron and proton, respectively, q is four-momentum transfer.

No description provided.

No description provided.

Jet rates for E+ P data.

Jet rates for E- P data.

Cross section for large rapidity gap events.

Measured angular distribution of the Lambda-Lambdabar pairs in all hadron events. Statistical errors only are given but the systematic errors are negigible. Data are read from the plot.

Measured rapidity difference os the Lambda and Lambdabars in the Lambda-Lambdabar pairs in all hadron events. Statistical errors only are given but the systematic errors are negigible. Data are read from the plot.

Upper limits on the production cross section times branching fraction of the b* RH hypothesis at a 95% CL. Dashed colored lines show the expected limits from the l+jets and all-hadronic channels, where the latter start at resonance masses of 1.4 TeV. The observed and expected limits from the combination are shown as solid and dashed black lines, respectively. The green and yellow bands show the 68 and 95% confidence intervals on the combined expected limits.

Upper limits on the production cross section times branching fraction of the b* VL hypothesis at a 95% CL. Dashed colored lines show the expected limits from the l+jets and all-hadronic channels, where the latter start at resonance masses of 1.4 TeV. The observed and expected limits from the combination are shown as solid and dashed black lines, respectively. The green and yellow bands show the 68 and 95% confidence intervals on the combined expected limits.

Upper limits on the production cross section times branching fraction of the B+b hypothesis at a 95% CL. Dashed colored lines show the expected limits from the l+jets and all-hadronic channels, where the latter start at resonance masses of 1.4 TeV. The observed and expected limits from the combination are shown as solid and dashed black lines, respectively. The green and yellow bands show the 68 and 95% confidence intervals on the combined expected limits.

Upper limits on the production cross section times branching fraction of the B+t hypothesis at a 95% CL. Dashed colored lines show the expected limits from the l+jets and all-hadronic channels, where the latter start at resonance masses of 1.4 TeV. The observed and expected limits from the combination are shown as solid and dashed black lines, respectively. The green and yellow bands show the 68 and 95% confidence intervals on the combined expected limits.

Bayesian upper exclusion limits on the ratio $g_{W'}/g_{W}$ for an SSM-like W$\prime$ boson are shown. The unity coupling ratio (blue dotted curve) corresponds to the SSM common benchmark. The 68 and 95% quantiles of the limits are represented by the green and yellow bands, respectively.

Bayesian lower exclusion limits on the NUGIM G(221) mixing angle $\cot(\theta_{E})$ are shown as a function of the W$\prime$ boson mass. The theoretically excluded region is shaded in grey. The 68 and 95% quantiles of the limits are represented by the green and yellow bands, respectively.

Bayesian upper exclusion limits at 95% CL on the product of the production cross section and branching fraction of a QBH in an associated $ au$ lepton and neutrino final state. Minimum threshold masses $m_{th}$ of up to 6.6 TeV are excluded at 95% CL. The observed limit (solid line) is obtained from the intersection with the LO QBH cross section (blue dotted curve). The 68 and 95% quantiles of the limits are represented by the green and yellow bands, respectively.

Bayesian upper limits at 95% CL on the cross section of the process $pp\rightarrow\tau\nu$ mediated via LQ exchange in the t-channel. The inner (green) band and the outer (yellow) band indicate the regions containing 68 and 95%, respectively, of the distribution of limits expected under the background-only hypothesis. The predicted LQ cross section at LO in the three coupling benchmark scenarios is depicted in different colors for $g_{U}=1$. The uncertainty bandy correspond to the sum in quadrature of PDF and scale variations. The first benchmark scenario considers only couplings to left-handed SM fermions (i.e. $\beta_{\text{R}}^{ij}=0$) and is referred to as "best fit LH". The second benchmark, referred to as "best fit LH+RH", considers $|\beta_{\text{R}}^{\text{b}\tau}|=1$ and all other $\beta_{\text{R}}^{ij}=0$. In the third "democratic" benchmark, equal couplings only to LH fermions are assumed, i.e. $\beta_{\text{L}}^{ij}=1$ and $\beta_{\text{R}}^{ij}=0$ for all $i$ and $j$.

Expected and observed lower limits of the LQ mass as a function of the coupling $g_{U}$ in the LH scenario. The blue band shows the 68% and 95% regions of $g_{U}$ preferred by the fit to the b anomalies data.

Expected and observed lower limits of the LQ mass as a function of the coupling $g_{U}$ in the LH+RH scenario. The blue band shows the 68% and 95% regions of $g_{U}$ preferred by the fit to the b anomalies data.

Expected and observed lower limits of the LQ mass as a function of the coupling $g_{U}$ in the democratic scenario. The blue band shows the 68% and 95% regions of $g_{U}$ preferred by the fit to the b anomalies data.

Bayesian upper exclusion limits at 95% CL on each of the Wilson coefficients described by the EFT model based on 2016-2018 data. The three different coupling types represent a left-handed vector coupling ($\epsilon^{cb}_{L}$), tensor-like coupling ($\epsilon^{cb}_{T}$), and scalar-tensor-like coupling ($\epsilon^{cb}_{S_{L}}$). The 68 and 95% quantiles of the limits are represented by the green and yellow bands, respectively.

Summary of exclusion limits (expected and observed) calculated at 95% CL for full Run-2 CMS data.

Background prediction and observed data yields in the signal region bins. The background yields are obtained from the background-only fit and serve as input to the simplified likelihood reinterpretation scheme. The naming of the bins is "year_binnumber", following the binning from Figure 4.

Matrix of covariance coefficients between signal region bins. The coefficients are obtained from the background-only fit and serve as input to the simplified likelihood reinterpretation scheme. The naming of the bins is "year_binnumber", following the binning used in Figure 4.

Predicted signal yields for the 2017 data-taking period, corresponding to 41.3 fb$^{-1}$, after the application of each search requirement (cumulative) for various signal hypothesis. The requirements listed are presented as total efficiencies w.r.t. the previous selection step.

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.

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