Single vector leptoquark (LQv) k = 1 total selection efficiency, accounting for both the decay branching fraction and the event selection, for events that pass the signal region requirements and any of the top quark or b jet categories defined in the search.

Neutron multiplicity dependence of acoplanarity ($\alpha$) from process $\gamma\gamma$ to $\mu^+\mu^-$ in ultraperipheral PbPb at $\sqrt{s_{NN}}=5.02$ TeV.

Neutron multiplicity dependence of acoplanarity ($\alpha$) from process $\gamma\gamma$ to $\mu^+\mu^-$ in ultraperipheral PbPb at $\sqrt{s_{NN}}=5.02$ TeV.

Neutron multiplicity dependence of acoplanarity ($\alpha$) from process $\gamma\gamma$ to $\mu^+\mu^-$ in ultraperipheral PbPb at $\sqrt{s_{NN}}=5.02$ TeV.

Neutron multiplicity dependence of average core acoplanarity ($<\alpha^{core}>$) distributions from process $\gamma\gamma$ to $\mu^+\mu^-$ in ultraperipheral PbPb at $\sqrt{s_{NN}}=5.02$ TeV.

Neutron multiplicity dependence of average invariant mass ($<m_{\mu\mu}>$) from process $\gamma\gamma$ to $\mu^+\mu^-$ in ultraperipheral PbPb at $\sqrt{s_{NN}}=5.02$ TeV.

Acoplanarity ($\alpha$) distributions from process $\gamma\gamma$ to $\mu^+\mu^-$ for 0n1n class in ultraperipheral PbPb at $\sqrt{s_{NN}}=5.02$ TeV. Dimuon rapidity is in the hemisphere containing larger neutron multiplicity.

Acoplanarity ($\alpha$) distributions from process $\gamma\gamma$ to $\mu^+\mu^-$ for 0nXn class in ultraperipheral PbPb at $\sqrt{s_{NN}}=5.02$ TeV. Dimuon rapidity is in the hemisphere containing larger neutron multiplicity.

Acoplanarity ($\alpha$) distributions from process $\gamma\gamma$ to $\mu^+\mu^-$ for 1nXn class in ultraperipheral PbPb at $\sqrt{s_{NN}}=5.02$ TeV. Dimuon rapidity is in the hemisphere containing larger neutron multiplicity.

Acoplanarity ($\alpha$) distributions from process $\gamma\gamma$ to $\mu^+\mu^-$ for 0n1n class in ultraperipheral PbPb at $\sqrt{s_{NN}}=5.02$ TeV. Dimuon rapidity is in the hemisphere containing smaller neutron multiplicity.

Acoplanarity ($\alpha$) distributions from process $\gamma\gamma$ to $\mu^+\mu^-$ for 0nXn class in ultraperipheral PbPb at $\sqrt{s_{NN}}=5.02$ TeV. Dimuon rapidity is in the hemisphere containing smaller neutron multiplicity.

Acoplanarity ($\alpha$) distributions from process $\gamma\gamma$ to $\mu^+\mu^-$ for 1nXn class in ultraperipheral PbPb at $\sqrt{s_{NN}}=5.02$ TeV. Dimuon rapidity is in the hemisphere containing smaller neutron multiplicity.

Summary of the sources of systematic and statistical uncertainties and their impact on the measurement of the ttH and tH signal rates, and the measured value of the unconstrained nuisance parameters. The quantity $\Delta r_{x}/r_{x}$ corresponds to the change in uncertainty when fixing the nuisance parameters associated with that uncertainty in the fit. Under the label "MC and sideband statistical uncertainty" are the uncertainties associated with the limited number of simulated MC events and the amount of data events in the application region of the MP method.

Number of events selected in the 2$\ell$ss+0$\tau_{\mathrm{h}}$, 3$\ell$+0$\tau_{\mathrm{h}}$, and 2$\ell$ss+1$\tau_{\mathrm{h}}$ analysis channels compared to the event yields expected from the $ tH$ and $ H$ signals and from background processes. The expected event yields are computed for the values of nuisance parameters and of the POI obtained from the ML fit. The best fit values of the POI amount to $\hat{r_{ttH}} = 0.92$ and $\hat{r_{tH}} = 5.7$. Quoted uncertainties represent the sum of statistical and systematic components.

Number of events selected in the 2$\ell$ss+0$\tau_{\mathrm{h}}$, 3$\ell$+0$\tau_{\mathrm{h}}$, and 2$\ell$ss+1$\tau_{\mathrm{h}}$ analysis channels compared to the event yields expected from the $ tH$ and $ H$ signals and from background processes. The expected event yields are computed for the values of nuisance parameters and of the POI obtained from the ML fit. The best fit values of the POI amount to $\hat{r_{ttH}} = 0.92$ and $\hat{r_{tH}} = 5.7$. Quoted uncertainties represent the sum of statistical and systematic components.

Number of events selected in the 1$\ell$+1$\tau_{\mathrm{h}}$, 0$\ell$+2$\tau_{\mathrm{h}}$, 2$\ell$os+1$\tau_{\mathrm{h}}$, and 1$\ell$+2$\tau_{\mathrm{h}}$ analysis channels compared to the event yields expected from the $ tH$ and $ H$ signals and from background processes. The expected event yields are computed for the values of nuisance parameters and of the POI obtained from the ML fit. The best fit values of the POI amount to $\hat{r_{ttH}} = 0.92$ and $\hat{r_{tH}} = 5.7$. Quoted uncertainties represent the sum of statistical and systematic components.

Number of events selected in the 1$\ell$+1$\tau_{\mathrm{h}}$, 0$\ell$+2$\tau_{\mathrm{h}}$, 2$\ell$os+1$\tau_{\mathrm{h}}$, and 1$\ell$+2$\tau_{\mathrm{h}}$ analysis channels compared to the event yields expected from the $ tH$ and $ H$ signals and from background processes. The expected event yields are computed for the values of nuisance parameters and of the POI obtained from the ML fit. The best fit values of the POI amount to $\hat{r_{ttH}} = 0.92$ and $\hat{r_{tH}} = 5.7$. Quoted uncertainties represent the sum of statistical and systematic components.

Number of events selected in the 4$\ell$+0$\tau_{\mathrm{h}}$, 3$\ell$+1$\tau_{\mathrm{h}}$, and 2$\ell$+2$\tau_{\mathrm{h}}$ analysis channels compared to the event yields expected from the $ tH$ and $ H$ signals and from background processes. The expected event yields are computed for the values of nuisance parameters and of the POI obtained from the ML fit. The best fit values of the POI amount to $\hat{r_{ttH}} = 0.92$ and $\hat{r_{tH}} = 5.7$. Quoted uncertainties represent the sum of statistical and systematic components.

Number of events selected in the 4$\ell$+0$\tau_{\mathrm{h}}$, 3$\ell$+1$\tau_{\mathrm{h}}$, and 2$\ell$+2$\tau_{\mathrm{h}}$ analysis channels compared to the event yields expected from the $ tH$ and $ H$ signals and from background processes. The expected event yields are computed for the values of nuisance parameters and of the POI obtained from the ML fit. The best fit values of the POI amount to $\hat{r_{ttH}} = 0.92$ and $\hat{r_{tH}} = 5.7$. Quoted uncertainties represent the sum of statistical and systematic components.

Dependence of the likelihood function L as a function of $\kappa_{\mathrm{t}}$ and $\kappa_{\mathrm{V}}$.

Dependence of the likelihood function L as a function of $\kappa_{\mathrm{t}}$ and $\kappa_{\mathrm{V}}$.

$-2\Delta log( L)$ as a function of the production rate of ttH and tH.

$-2\Delta log( L)$ as a function of the production rate of ttH and tH.

$-2\Delta log( L)$ as a function of the production rate of ttH and ttW.

$-2\Delta log( L)$ as a function of the production rate of ttH and ttW.

$-2\Delta log( L)$ as a function of the production rate of ttH and ttZ.

$-2\Delta log( L)$ as a function of the production rate of ttH and ttZ.

$-2\Delta log( L)$ as a function of the production rate of ttW and ttZ.

$-2\Delta log( L)$ as a function of the production rate of ttW and ttZ.

Covariance matrix of the fitted parameters of interest.

Covariance matrix of the fitted parameters of interest.

Production rate of the ttH signal measured in each of the ten channels individually and for the combination of all channels

Production rate of the ttH signal measured in each of the ten channels individually and for the combination of all channels

Production rate of the tH signal measured in each of the ten channels individually and for the combination of all channels

Production rate of the tH signal measured in each of the ten channels individually and for the combination of all channels

Probability for ttH signal events to pass the event selection in each of the analysis categories employed.

Prompt $D^0$ meson $v_3$ as a function of $p_T$ in PbPb collisions at $\sqrt{s_{NN}}=5.02~TeV$.

Prompt $D^0$ meson $v_3$ as a function of $p_T$ in PbPb collisions at $\sqrt{s_{NN}}=5.02~TeV$.

Prompt $D^0$ meson $v_3$ as a function of $p_T$ in PbPb collisions at $\sqrt{s_{NN}}=5.02~TeV$.

Prompt $D^0$ meson $v_2$ as a function of $p_T$ in PbPb collisions at $\sqrt{s_{NN}}=5.02~TeV$.

Prompt $D^0$ meson $v_2$ as a function of $p_T$ in PbPb collisions at $\sqrt{s_{NN}}=5.02~TeV$.

Prompt $D^0$ meson $v_2$ as a function of $p_T$ in PbPb collisions at $\sqrt{s_{NN}}=5.02~TeV$.

Prompt $D^0$ meson $v_3$ as a function of $p_T$ in PbPb collisions at $\sqrt{s_{NN}}=5.02~TeV$.

Prompt $D^0$ meson $v_3$ as a function of $p_T$ in PbPb collisions at $\sqrt{s_{NN}}=5.02~TeV$.

Prompt $D^0$ meson $v_3$ as a function of $p_T$ in PbPb collisions at $\sqrt{s_{NN}}=5.02~TeV$.

Prompt $D^0$ meson $v_2$ as a function of centrality in PbPb collisions at $\sqrt{s_{NN}}=5.02~TeV$.

Prompt $D^0$ meson $v_2$ as a function of centrality in PbPb collisions at $\sqrt{s_{NN}}=5.02~TeV$.

Prompt $D^0$ meson $v_3$ as a function of centrality in PbPb collisions at $\sqrt{s_{NN}}=5.02~TeV$.

Prompt $D^0$ meson $v_3$ as a function of centrality in PbPb collisions at $\sqrt{s_{NN}}=5.02~TeV$.

Prompt $D^0$ meson $v_2$ as a function of rapidity in PbPb collisions at $\sqrt{s_{NN}}=5.02~TeV$.

Prompt $D^0$ meson $v_3$ as a function of rapidity in PbPb collisions at $\sqrt{s_{NN}}=5.02~TeV$.

Prompt $D^0$ meson $\Delta v_2$ as a function of rapidity in PbPb collisions at $\sqrt{s_{NN}}=5.02~TeV$.

Measured fiducial cross sections for the $\mathrm{W}^\pm_{\mathrm{L}}\mathrm{W}^\pm_{\mathrm{L}}$ and $\mathrm{W}^\pm_{\mathrm{X}}\mathrm{W}^\pm_{\mathrm{T}}$ processes, and for the $\mathrm{W}^\pm_{\mathrm{L}}\mathrm{W}^\pm_{\mathrm{X}}$ and $\mathrm{W}^\pm_{\mathrm{T}}\mathrm{W}^\pm_{\mathrm{T}}$ processes for the helicity eigenstates defined in the parton-parton center-of-mass frame. The combination of the statistical and systematic uncertainties is shown. $\mathcal{B}$ is the branching fraction for $\mathrm{W}\mathrm{W} \rightarrow \ell \nu \ell' \nu$. The fiducial region is defined by requiring two same-sign leptons with $p_{T}>20$, $|\eta|<2.5$, and $m_{ll}>20$, and two jets with $m_{jj}>500$ and $|\Delta \eta_{jj}|>2.5$. The jets at generator level are clustered from stable particles, excluding neutrinos, using the anti-kt clustering algorithm with R = 0.4, and are required to have $p_{T}>50$ and $|\eta|<4.7$. The jets within $\Delta R<0.4$ of the selected charged leptons are not included. The theoretical predictions including the $\mathcal{O}(\alpha_{s}\alpha^6)$ and $\mathcal{O}(\alpha^7)$ corrections to the \MGvATNLO LO cross sections, as described in arXiv:2009.09429, are also shown. The theoretical uncertainties include statistical, PDF, and LO scale uncertainties.

Results of elliptic flow, corrected for short range correlations, for prompt neutral D mesons, as a function of multiplicity for $|y_{lab}|< 1$, with 6$ < p_{T} < $8 GeV in pp collisions at 13 TeV.

Results of elliptic flow, corrected for short range correlations, for prompt neutral D mesons, as a function of multiplicity for $|y_{lab}|< 1$, with 2$ < p_{T} < $4 GeV in pPb collisions at 8.16 TeV.

Results of elliptic flow, corrected for short range correlations, for prompt neutral D mesons, as a function of multiplicity for $|y_{lab}|< 1$, with 4$ < p_{T} < $6 GeV in pPb collisions at 8.16 TeV.

Results of elliptic flow, corrected for short range correlations, for prompt neutral D mesons, as a function of multiplicity for $|y_{lab}|< 1$, with 6$ < p_{T} < $8 GeV in pPb collisions at 8.16 TeV.

Results of elliptic flow, corrected for short range correlations, for nonprompt neutral D mesons, as a function of transverse momenta for $|y_{lab}|< 1$, with $185 \leq N^{offline}_{trk} < 250$ in pPb collisions at 8.16 TeV.

The observed DNN output distribution for data collected in 2016 in the VBF-SB region compared to the post-fit background estimate from SM processes. The predicted backgrounds are obtained from a S+B fit performed across analysis regions and years. The total post-fit and pre-fit uncertainties on the background prediction are also reported.

The observed DNN output distribution for data collected in 2017 in the VBF-SB region compared to the post-fit background estimate from SM processes. The predicted backgrounds are obtained from a S+B fit performed across analysis regions and years. The total post-fit and pre-fit uncertainties on the background prediction are also reported.

The observed DNN output distribution for data collected in 2018 in the VBF-SB region compared to the post-fit background estimate from SM processes. The predicted backgrounds are obtained from a S+B fit performed across analysis regions and years. The total post-fit and pre-fit uncertainties on the background prediction are also reported.

The observed DNN output distribution in the VBF-SB regions for the combination of 2016, 2017, and 2018 data, compared to the post-fit prediction from SM processes. The total post-fit and pre-fit uncertainties on the background prediction are also reported.

The observed DNN output distribution in the VBF-SR region for the combination of 2016, 2017, and 2018 data, compared to the post-fit prediction from SM processes. The post-fit distributions for the Higgs boson signal produced via ggH and VBF modes with mass of 125.38 GeV, along with the total post-fit and pre-fit uncertainties on the background prediction, are also reported.

The observed BDT output distribution compared to the prediction from the simulation of various SM background processes. Dimuon events passing the event selection requirements of the ggH category, with dimuon mass between 110–150 GeV, are considered. The expected distributions for ggH, VBF, and other signal processes are also reported. The background uncertainty is obtained by summing in quadrature the uncertainty due to the limited size of the simulated samples with those arising from experimental and theoretical sources of systematic uncertainty.

Comparison between the data and the total background extracted from a S+B fit performed in the ggH-cat1 subcategory. The one and two standard deviation bands include the uncertainties in the background component of the fit. The best fit signal contribution with mass of 125.38 GeV is also reported.

Comparison between the data and the total background extracted from a S+B fit performed in the ggH-cat2 subcategory. The one and two standard deviation bands include the uncertainties in the background component of the fit. The best fit signal contribution with mass of 125.38 GeV is also reported.

Comparison between the data and the total background extracted from a S+B fit performed in the ggH-cat3 subcategory. The one and two standard deviation bands include the uncertainties in the background component of the fit. The best fit signal contribution with mass of 125.38 GeV is also reported.

Comparison between the data and the total background extracted from a S+B fit performed in the ggH-cat4 subcategory. The one and two standard deviation bands include the uncertainties in the background component of the fit. The best fit signal contribution with mass of 125.38 GeV is also reported.

Comparison between the data and the total background extracted from a S+B fit performed in the ggH-cat5 subcategory. The one and two standard deviation bands include the uncertainties in the background component of the fit. The best fit signal contribution with mass of 125.38 GeV is also reported.

The observed BDT output distribution compared to the prediction from the simulation of various SM background processes. Dimuon events passing the event selection requirements of the ttH hadronic category, with dimuon mass between 110–150 GeV, are considered. The expected distributions from different production modes of the Higgs boson are also reported. The background uncertainty is obtained by summing in quadrature the uncertainty due to the limited size of the simulated samples with those arising from experimental and theoretical sources of systematic uncertainty.

The observed BDT output distribution compared to the prediction from the simulation of various SM background processes. Dimuon events passing the event selection requirements of the ttH leptonic category, with dimuon mass between 110–150 GeV, are considered. The expected distributions from different production modes of the Higgs boson are also reported. The background uncertainty is obtained by summing in quadrature the uncertainty due to the limited size of the simulated samples with those arising from experimental and theoretical sources of systematic uncertainty.

Comparison between the data and the total background extracted from a S+B fit performed in the ttHhad-cat1 subcategory. The one and two standard deviation bands include the uncertainties in the background component of the fit. The best fit signal contribution with mass of 125.38 GeV is also reported.

Comparison between the data and the total background extracted from a S+B fit performed in the ttHhad-cat2 subcategory. The one and two standard deviation bands include the uncertainties in the background component of the fit. The best fit signal contribution with mass of 125.38 GeV is also reported.

Comparison between the data and the total background extracted from a S+B fit performed in the ttHhad-cat3 subcategory. The one and two standard deviation bands include the uncertainties in the background component of the fit. The best fit signal contribution with mass of 125.38 GeV is also reported.

Comparison between the data and the total background extracted from a S+B fit performed in the ttHlep-cat1 subcategory. The one and two standard deviation bands include the uncertainties in the background component of the fit. The best fit signal contribution with mass of 125.38 GeV is also reported.

Comparison between the data and the total background extracted from a S+B fit performed in the ttHlep-cat2 subcategory. The one and two standard deviation bands include the uncertainties in the background component of the fit. The best fit signal contribution with mass of 125.38 GeV is also reported.

The observed BDT output distribution compared to the prediction from the simulation of various SM background processes. Dimuon events passing the event selection requirements of the WH leptonic category, with dimuon mass between 110–150 GeV, are considered. The expected distributions from different production modes of the Higgs boson are also reported. The background uncertainty is obtained by summing in quadrature the uncertainty due to the limited size of the simulated samples with those arising from experimental and theoretical sources of systematic uncertainty.

The observed BDT output distribution compared to the prediction from the simulation of various SM background processes. Dimuon events passing the event selection requirements of the ZH leptonic category, with dimuon mass between 110–150 GeV, are considered. The expected distributions from different production modes of the Higgs boson are also reported. The background uncertainty is obtained by summing in quadrature the uncertainty due to the limited size of the simulated samples with those arising from experimental and theoretical sources of systematic uncertainty.

Comparison between the data and the total background extracted from a S+B fit performed in the WH-cat1 subcategory. The one and two standard deviation bands include the uncertainties in the background component of the fit. The best fit signal contribution with mass of 125.38 GeV is also reported.

Comparison between the data and the total background extracted from a S+B fit performed in the WH-cat2 subcategory. The one and two standard deviation bands include the uncertainties in the background component of the fit. The best fit signal contribution with mass of 125.38 GeV is also reported.

Comparison between the data and the total background extracted from a S+B fit performed in the WH-cat3 subcategory. The one and two standard deviation bands include the uncertainties in the background component of the fit. The best fit signal contribution with mass of 125.38 GeV is also reported.

Comparison between the data and the total background extracted from a S+B fit performed in the ZH-cat1 subcategory. The one and two standard deviation bands include the uncertainties in the background component of the fit. The best fit signal contribution with mass of 125.38 GeV is also reported.

Comparison between the data and the total background extracted from a S+B fit performed in the ZH-cat2 subcategory. The one and two standard deviation bands include the uncertainties in the background component of the fit. The best fit signal contribution with mass of 125.38 GeV is also reported.

Observed local p-values as a function of the Higgs boson mass hypothesis extracted from the combined fit across event categories as well as from each individual production category.

Expected local p-values as a function of the Higgs boson mass hypothesis extracted from the combined fit across event categories as well as from each individual production category. They are calculated using the background-only expectation obtained from the S+B fit and injecting a signal with mass of 125.38 GeV and strenght of one.

Signal strength modifiers measured for a Higgs boson with mass of 125.38 GeV in each production category as well as from the combined fit.

Scan of the profiled likelihood ratio as a function of $\mu_{ggH,t\overline{t}H}$ and $\mu_{VBF,VH}$.

The 68% CL contour for the profiled likelihood scan as a function of $\mu_{ggH,t\overline{t}H}$ and $\mu_{VBF,VH}$.

The 95% CL contour for the profiled likelihood scan as a function of $\mu_{ggH,t\overline{t}H}$ and $\mu_{VBF,VH}$.

The mass distribution for the weighted combination of VBF-SB and VBF-SR events. Each event is weighted proportionally to the S/(S+B) ratio, calculated as a function of the mass-decorrelated DNN output. The total post-fit and pre-fit uncertainties on the background prediction, along with the best fit signal contribution, are also reported.

Comparison between the data and the total background extracted for the S/(S+B) weighted combination of all event categories. The one and two standard deviation bands include the uncertainties in the background component of the fit. The best fit signal contribution with mass of 125.38 GeV is also reported.

Observed and expected local p-values as a function of the Higgs boson mass hypothesis extracted from the combined fit performed on data recorded at centre-of-mass energies of 7, 8, and 13 TeV. The expected p-values are calculated using the background-only expectation obtained from the S+B fit and injecting a signal with mass of 125.38 GeV and strenght of one.

The observed profile likelihood ratio as a function of $\kappa_{\mu}$ for a Higgs boson mass of 125.38 GeV, obtained from a combined fit with Ref.[10] in the resolved $\kappa$-framework. The best fit value for $\kappa_{\mu}$ is 1.07 and the corresponding observed 68% CL interval is $0.85 < \kappa_{\mu} < 1.29$.

The best fit estimates for the reduced coupling modifiers extracted for fermions and weak bosons from the resolved $\kappa$-framework compared to their corresponding prediction from the SM. The error bars represent 68% CL intervals for the measured parameters.

Differential cross sections normalized to the fiducial cross section for the combined 4e, 2e2µ, and 4µ decay channels as a function of the invariant mass of the ZZ system

Differential cross sections normalized to the fiducial cross section for the combined 4e, 2e2µ, and 4µ decay channels as a function of azimuthal separation of the two Z bosons

Differential cross sections normalized to the fiducial cross section for the combined 4e, 2e2µ, and 4µ decay channels as a function of ∆R separation of the two Z bosons

Measured fiducial cross section for each data sample and combined. The first uncertainty is statistical, the second is experimental systematic, and the third is associated with the integrated luminosity.

Measured total σ(pp → ZZ) cross section for each data sample and combined. The first uncertainty is statistical, the second is experimental systematic, the third is theoretical systematic. The fourth uncertainty is associated with the integrated luminosity.

Expected and observed one-dimensional 95% CL limits on aTGC parameters. The corresponding constrains on EFT parameters are estimated using the transformation from Ref. [65].

Expected and observed one-dimensional 95% CL limits on aTGC parameters. The corresponding constrains on EFT parameters are estimated using the transformation from Ref. [65].