Invariant mass distribution of final state particles in SR1 HF category ($\text{H}\to\text{J}/\psi\gamma$ signal)

Invariant mass distribution of final state particles in SR1 Z-LP category ($\text{Z}\to\text{J}/\psi\gamma$ signal)

Invariant mass distribution of final state particles in SR1 Z-HP category ($\text{Z}\to\text{J}/\psi\gamma$ signal)

Invariant mass distribution of final state particles in SR2 category ($\text{H}\to\psi(\text{2S})\gamma$ signal)

Invariant mass distribution of final state particles in SR2 category ($\text{H}\to\psi(\text{2S})\gamma$ signal)

Invariant mass distribution of final state particles in control region category ($\text{H}\to\mu\mu\gamma$ background)

Observed and expected (with $\pm 1$, $\pm 2$ standard deviation bands) exclusion limits at 95% C.L. on the branching fraction $\mathcal{B}$ of the $(\text{H,Z}) \to \psi(\text{nS})\gamma$ decays

Observed (expected) upper limits at 95% C.L. on the normalized values with respect to the SM expectation, denoted as the signal strength parameter $\mu$, of the product of the cross section $\sigma$ and the branching fraction $\mathcal{B}$ of the $(\text{H,Z}) \to \psi(\text{nS})\gamma$ decays, and the branching fraction, assuming a SM Z and H boson cross section.

Post-fit agreement between data and MC prediction for the $SR_{\mathrm{loose}}^{2\ell3j}$ signal region, which uses the invariant mass of the c-tagged jet and the geometrically closest b-tagged jet, $m_{\mathit{cb}}^{\mathrm{min}\Delta R}$, as an observable. The hatched uncertainty bands include all uncertainties and their correlations. The last bins contain overflow events. "Other Top" includes single-top-quark production and associated production of $t\bar{t}$ and single top quarks with bosons. "Non-Top" includes W+jets, Z+jets, and diboson processes.

Post-fit agreement between data and MC prediction for the $SR_{\mathrm{loose}}^{1\ell6ji}$ signal region, which uses the jet multiplicity, $N_{\mathrm{jets}}$, as an observable. The hatched uncertainty bands include all uncertainties and their correlations. The last bins contain overflow events. "Other Top" includes single-top-quark production and associated production of $t\bar{t}$ and single top quarks with bosons. "Non-Top" includes W+jets, Z+jets, and diboson processes.

Post-fit agreement between data and MC prediction for the $SR_{\mathrm{tight}}^{1\ell6ji}$ signal region, which uses the jet multiplicity, $N_{\mathrm{jets}}$, as an observable. The hatched uncertainty bands include all uncertainties and their correlations. The last bins contain overflow events. "Other Top" includes single-top-quark production and associated production of $t\bar{t}$ and single top quarks with bosons. "Non-Top" includes W+jets, Z+jets, and diboson processes.

Post-fit agreement between data and MC prediction for the $SR_{\mathrm{loose}}^{2\ell4ji}$ signal region, which uses the jet multiplicity, $N_{\mathrm{jets}}$, as an observable. The hatched uncertainty bands include all uncertainties and their correlations. The last bins contain overflow events. "Other Top" includes single-top-quark production and associated production of $t\bar{t}$ and single top quarks with bosons. "Non-Top" includes W+jets, Z+jets, and diboson processes.

Post-fit agreement between data and MC prediction for the $SR_{\mathrm{tight}}^{2\ell4ji}$ signal region, which uses the jet multiplicity, $N_{\mathrm{jets}}$, as an observable. The hatched uncertainty bands include all uncertainties and their correlations. The last bins contain overflow events. "Other Top" includes single-top-quark production and associated production of $t\bar{t}$ and single top quarks with bosons. "Non-Top" includes W+jets, Z+jets, and diboson processes.

Jet multiplicity and tagging criteria for the SRs and CRs in the single-lepton and dilepton channels. All single-lepton regions are defined in the 5-jet-exclusive and 6-jet-inclusive jet selections, all dilepton regions in the 3-jet-exclusive and 4-jet-inclusive jet selections. By construction, $\mathrm{SR}^{2\ell}_{\mathrm{tight}}$ only exists for the 4-jet-inclusive selection. Identical requirements on inclusive and exclusive working points, e.g., one c@22% and one c@11% in the $\mathrm{CR}_1^{1\ell}$, indicate a veto of additional tags at the inclusive working point.

Breakdown of the fiducial $t\bar{t}+{\geq}2c$ and $t\bar{t}+1c$ cross-section uncertainties. Systematic uncertainties are grouped into signal and background modeling, instrumental, and MC statistics categories. The table lists the fractional uncertainty in the measured cross-sections in percent and is estimated from the covariance matrix of the profile likelihood fit, following EPJC 84 (2024) 593. Individual groups of uncertainties can be added in quadrature to obtain the total uncertainty. JES and JER denote the jet energy scale and resolution uncertainties, respectively. The fractional uncertainties in the fitted $t\bar{t}+{\geq}1c$ and $t\bar{t}+\text{light}$ normalization factors are included in the data statistics category. For presentation purposes, all uncertainties are symmetrized.

Measured fiducial cross-section values in comparison with various NLO+PS predictions. The uncertainties in the predictions include independent and simultaneous variations of the $\mu_R$ and $\mu_F$ scales and uncertainties in the choice of the PDF set, but no uncertainties in the parton shower or hadronization. The uncertainties in the measured values include all statistical and systematic uncertainties. For presentation purposes, the quoted measurement and prediction uncertainties are symmetrized.

Measured and predicted values for the $t\bar{t}+{\geq}1b$, $t\bar{t}+{\geq}2c$, and $t\bar{t}+1c$ cross-sections and for the cross-section ratios of these processes to total $t\bar{t}+\text{jets}$ production. The quoted uncertainties in the measurements include statistical and systematic uncertainties. The predictions are taken from the nominal setup using inclusive $t\bar{t}$ Powheg+Pythia8 predictions for the $t\bar{t}+{\geq}2c$, $t\bar{t}+1c$, and $t\bar{t}+\text{light}$ processes, and $t\bar{t}+b\bar{b}$ Powheg+Pythia8 predictions for the $t\bar{t}+{\geq}1b$ process. The first use the 5FS, the latter the 4FS, where the $b$-quarks are treated as massive particles and are included in the matrix element calculation. The uncertainties in the predictions include independent and simultaneous variations of the $\mu_R$ and $\mu_F$ scales and uncertainties in the choice of the PDF set.

Nuisance parameters with a large impact on the fiducial measurement: ranking of the ten most impactful nuisance parameters for the $t\bar{t}+{\geq}2c$ cross-section measurement. The colored boxes indicate the impact of the nuisance parameters, $\Delta\mu$ (top axis), on the corresponding parameter of interest, $\mu$. They are calculated from the covariance matrix using the post-fit up (down) uncertainties of the parameter of interest and the nuisance parameter, $\sigma_{\mathrm{POI}}^{\mathrm{up}}$ and $\sigma_{\mathrm{NP}}^{\mathrm{up}}$ ($\sigma_{\mathrm{POI}}^{\mathrm{down}}$ and $\sigma_{\mathrm{NP}}^{\mathrm{down}}$), and their linear correlation coefficient, $\rho$, following EPJC 84 (2024) 593. The data points and error bars indicate the post-fit values and uncertainties of the nuisance parameters relative to their pre-fit values, $\theta_0$, in units of their pre-fit uncertainties, $\Delta \theta$ (bottom axis). A hollow circle indicates that the nuisance parameter has a Poisson constraint.

Nuisance parameters with a large impact on the fiducial measurement: ranking of the ten most impactful nuisance parameters for the $t\bar{t}+1c$ cross-section measurement. The colored boxes indicate the impact of the nuisance parameters, $\Delta\mu$ (top axis), on the corresponding parameter of interest, $\mu$. They are calculated from the covariance matrix using the post-fit up (down) uncertainties of the parameter of interest and the nuisance parameter, $\sigma_{\mathrm{POI}}^{\mathrm{up}}$ and $\sigma_{\mathrm{NP}}^{\mathrm{up}}$ ($\sigma_{\mathrm{POI}}^{\mathrm{down}}$ and $\sigma_{\mathrm{NP}}^{\mathrm{down}}$), and their linear correlation coefficient, $\rho$, following EPJC 84 (2024) 593. The data points and error bars indicate the post-fit values and uncertainties of the nuisance parameters relative to their pre-fit values, $\theta_0$, in units of their pre-fit uncertainties, $\Delta \theta$ (bottom axis). A hollow circle indicates that the nuisance parameter has a Poisson constraint.

Nuisance parameters with a large impact on the fiducial measurement: ranking of the ten most impactful nuisance parameters for the $t\bar{t}+{\geq}2c$ ratio to total $t\bar{t}+\text{jets}$ production. The colored boxes indicate the impact of the nuisance parameters, $\Delta\mu$ (top axis), on the corresponding parameter of interest, $\mu$. They are calculated from the covariance matrix using the post-fit up (down) uncertainties of the parameter of interest and the nuisance parameter, $\sigma_{\mathrm{POI}}^{\mathrm{up}}$ and $\sigma_{\mathrm{NP}}^{\mathrm{up}}$ ($\sigma_{\mathrm{POI}}^{\mathrm{down}}$ and $\sigma_{\mathrm{NP}}^{\mathrm{down}}$), and their linear correlation coefficient, $\rho$, following EPJC 84 (2024) 593. The data points and error bars indicate the post-fit values and uncertainties of the nuisance parameters relative to their pre-fit values, $\theta_0$, in units of their pre-fit uncertainties, $\Delta \theta$ (bottom axis). A hollow circle indicates that the nuisance parameter has a Poisson constraint.

Nuisance parameters with a large impact on the fiducial measurement: ranking of the ten most impactful nuisance parameters for the $t\bar{t}+1c$ ratio to total $t\bar{t}+\text{jets}$ production. The colored boxes indicate the impact of the nuisance parameters, $\Delta\mu$ (top axis), on the corresponding parameter of interest, $\mu$. They are calculated from the covariance matrix using the post-fit up (down) uncertainties of the parameter of interest and the nuisance parameter, $\sigma_{\mathrm{POI}}^{\mathrm{up}}$ and $\sigma_{\mathrm{NP}}^{\mathrm{up}}$ ($\sigma_{\mathrm{POI}}^{\mathrm{down}}$ and $\sigma_{\mathrm{NP}}^{\mathrm{down}}$), and their linear correlation coefficient, $\rho$, following EPJC 84 (2024) 593. The data points and error bars indicate the post-fit values and uncertainties of the nuisance parameters relative to their pre-fit values, $\theta_0$, in units of their pre-fit uncertainties, $\Delta \theta$ (bottom axis). A hollow circle indicates that the nuisance parameter has a Poisson constraint.

Efficiencies and rejection rates for the $b$- and $c$-jet tagging working points of the $b/c$-tagger, as well as the selection cuts on the discriminants, $\mathcal{D}_b'$ and $\mathcal{D}_c'$ that define the tagging bins. The primes indicate that a standard logistic function is applied to the discriminants. The b@70% and c@22% working points are inclusive of their respective tighter working points. The values were estimated from single-lepton and dilepton $t\bar{t}$ events simulated with Powheg + Pythia 8.

Efficiencies and rejection rates for the $b$- and $c$-jet tagging working points of the $b/c$-tagger, as well as the selection cuts on the discriminants, $\mathcal{D}_b'$ and $\mathcal{D}_c'$ that define the tagging bins. The primes indicate that a standard logistic function is applied to the discriminants. The b@70% and c@22% working points are inclusive of their respective tighter working points. The values were estimated from single-lepton and dilepton $t\bar{t}$ events simulated with Powheg + Pythia 8.

Pre-fit event yields in the twelve single-lepton and dilepton control regions (CRs) of the analysis. The quoted uncertainties include both statistical and systematic components, but all normalization factors are fixed to their pre-fit values with no uncertainties. All uncertainties are assumed to be uncorrelated. "Other Top" includes single-top-quark production and associated production of $t\bar{t}$ and single top quarks with bosons. "Non-Top" includes W+jets, Z+jets, and diboson processes. The "Fakes" category includes events from misidentified leptons and non-prompt leptons, estimated through the data-driven matrix method in the single-lepton channel and through MC predictions in the dilepton channel.

Pre-fit event yields in the seven signal regions (SRs) of the analysis. The quoted uncertainties include both statistical and systematic components, but all normalization factors are fixed to their pre-fit values with no uncertainties. All uncertainties are assumed to be uncorrelated. "Other Top" includes single-top-quark production and associated production of $t\bar{t}$ and single top quarks with bosons. "Non-Top" includes W+jets, Z+jets, and diboson processes. The "Fakes" category includes events from misidentified leptons and non-prompt leptons, estimated through the data-driven matrix method in the single-lepton channel and through MC predictions in the dilepton channel.

Post-fit event yields in the twelve single-lepton and dilepton control regions (CRs) of the analysis. The quoted uncertainties include both statistical and systematic components, and post-fit correlations between uncertainties are taken into account. "Other Top" includes single-top-quark production and associated production of $t\bar{t}$ and single top quarks with bosons. "Non-Top" includes W+jets, Z+jets, and diboson processes. The "Fakes" category includes events from misidentified leptons and non-prompt leptons, estimated through the data-driven matrix method in the single-lepton channel and through MC predictions in the dilepton channel.

Post-fit event yields in the seven signal regions (SRs) of the analysis. The quoted uncertainties include both statistical and systematic components, and post-fit correlations between uncertainties are taken into account. "Other Top" includes single-top-quark production and associated production of $t\bar{t}$ and single top quarks with bosons. "Non-Top" includes W+jets, Z+jets, and diboson processes. The "Fakes" category includes events from misidentified leptons and non-prompt leptons, estimated through the data-driven matrix method in the single-lepton channel and through MC predictions in the dilepton channel.

Observed 95% CL upper limits on the production cross section times branching ratio of the A → ZH → Ztt process in the (mA , mH) plane.

Median expected (dashed lines) and observed (filled contours) 95% CL exclusion regions in the (mH, mA) parameter space of the type-II 2HDM for tan β = 0.5 (blue), 1 (orange), and 1.5 (red), derived assuming narrow resonances. The enclosed regions are excluded. The dotted lines indicate points of constant natural line width of the A boson relative to its mass in the 2HDM.

Median expected (dashed lines) and observed (filled contours) 95% CL exclusion regions in the (mH, mA) parameter space of the type-II 2HDM for tan β = 0.5 (blue), 1 (orange), and 1.5 (red), derived assuming narrow resonances. The enclosed regions are excluded. The dotted lines indicate points of constant natural line width of the A boson relative to its mass in the 2HDM.

Median expected (dashed lines) and observed (filled contours) 95% CL exclusion regions in the (mH, mA) parameter space of the type-II 2HDM for tan β = 0.5 (blue), 1 (orange), and 1.5 (red), derived assuming narrow resonances. The enclosed regions are excluded. The dotted lines indicate points of constant natural line width of the A boson relative to its mass in the 2HDM.

Median expected (dashed lines) and observed (filled contours) 95% CL exclusion regions in the (mH, mA) parameter space of the type-II 2HDM for tan β = 0.5 (blue), 1 (orange), and 1.5 (red), derived assuming narrow resonances. The enclosed regions are excluded. The dotted lines indicate points of constant natural line width of the A boson relative to its mass in the 2HDM.

Median expected (dashed lines) and observed (filled contours) 95% CL exclusion regions in the (mH, mA) parameter space of the type-II 2HDM for tan β = 0.5 (blue), 1 (orange), and 1.5 (red), derived assuming narrow resonances. The enclosed regions are excluded. The dotted lines indicate points of constant natural line width of the A boson relative to its mass in the 2HDM.

Median expected (dashed lines) and observed (filled contours) 95% CL exclusion regions in the (mH, mA) parameter space of the type-II 2HDM for tan β = 0.5 (blue), 1 (orange), and 1.5 (red), derived assuming narrow resonances. The enclosed regions are excluded. The dotted lines indicate points of constant natural line width of the A boson relative to its mass in the 2HDM.

Median expected (dashed lines) and observed (filled contours) 95% CL exclusion regions in the type-II 2HDM (tan β, mA) parameter space at mH = 400 GeV, derived assuming narrow resonances. The enclosed regions are excluded. The dotted lines indicate points of constant natural line width of the A boson relative to its mass in the 2HDM.

Median expected (dashed lines) and observed (filled contours) 95% CL exclusion regions in the type-II 2HDM (tan β, mA) parameter space at mH = 400 GeV, derived assuming narrow resonances. The enclosed regions are excluded. The dotted lines indicate points of constant natural line width of the A boson relative to its mass in the 2HDM.

Median expected (dashed lines) and observed (filled contours) 95% CL exclusion regions in the type-II 2HDM (tanβ, cos(β−α)) parameter space at mA = 600 GeV and mH = 400 GeV, derived assuming narrow resonances. The enclosed regions are excluded. The dotted lines indicate points of constant natural line width of the A boson relative to its mass in the 2HDM.

Median expected (dashed lines) and observed (filled contours) 95% CL exclusion regions in the type-II 2HDM (tanβ, cos(β−α)) parameter space at mA = 600 GeV and mH = 400 GeV, derived assuming narrow resonances. The enclosed regions are excluded. The dotted lines indicate points of constant natural line width of the A boson relative to its mass in the 2HDM.

Distributions of $m_T$ in the $W^{+}$ signal selection for mu final states for the pp collisions at $\sqrt{s}=$ 13TeV after the maximum likelihood fit. The EW backgrounds include the contributions from DY, $W\to\tau\nu$, and diboson processes.

Distributions of $m_T$ in the $W^{-}$ signal selection for e final states for the pp collisions at $\sqrt{s}=$ 5TeV after the maximum likelihood fit. The EW backgrounds include the contributions from DY, $W\to\tau\nu$, and diboson processes.

Distributions of $m_T$ in the $W^{-}$ signal selection for mu final states for the pp collisions at $\sqrt{s}=$ 5TeV after the maximum likelihood fit. The EW backgrounds include the contributions from DY, $W\to\tau\nu$, and diboson processes.

Distributions of $m_T$ in the $W^{-}$ signal selection for e final states for the pp collisions at $\sqrt{s}=$ 13TeV after the maximum likelihood fit. The EW backgrounds include the contributions from DY, $W\to\tau\nu$, and diboson processes.

Distributions of $m_T$ in the $W^{-}$ signal selection for mu final states for the pp collisions at $\sqrt{s}=$ 13TeV after the maximum likelihood fit. The EW backgrounds include the contributions from DY, $W\to\tau\nu$, and diboson processes.

Distributions of $m_{ll}$ in the Z signal selection for ee final states for the pp collisions at $\sqrt{s}=$ 5TeV after the maximum likelihood fit. The EW backgrounds include the contributions from DY, $W\to l\nu$, and diboson processes.

Distributions of $m_{ll}$ in the Z signal selection for mumu final states for the pp collisions at $\sqrt{s}=$ 5TeV after the maximum likelihood fit. The EW backgrounds include the contributions from DY, $W\to l\nu$, and diboson processes.

Distributions of $m_{ll}$ in the Z signal selection for ee final states for the pp collisions at $\sqrt{s}=$ 13TeV after the maximum likelihood fit. The EW backgrounds include the contributions from DY, $W\to l\nu$, and diboson processes.

Distributions of $m_{ll}$ in the Z signal selection for mumu final states for the pp collisions at $\sqrt{s}=$ 13TeV after the maximum likelihood fit. The EW backgrounds include the contributions from DY, $W\to l\nu$, and diboson processes.

Figure 4 Fiducial cross sections and ratios for $W\to\ell\nu$ and $Z\to\ell\ell$ processes at 13TeV

Figure 5 Fiducial cross sections and ratios for $W\to\ell\nu$ and $Z\to\ell\ell$ processes at 5.02TeV

Figure 6 Fiducial cross sections and ratios for $W\to\ell\nu$ and $Z\to\ell\ell$ processes between 13TeV and 5.02TeV

Figure 7 Total cross sections and ratios for $W\to\ell\nu$ and $Z\to\ell\ell$ processes at 13TeV

Figure 8 Total cross sections and ratios for $W\to\ell\nu$ and $Z\to\ell\ell$ processes at 5.02TeV

Figure 9 Total cross sections and ratios for $W\to\ell\nu$ and $Z\to\ell\ell$ processes between 13TeV and 5.02TeV

Summary of theoretical predictions of the total inclusive cross sections for W and Z boson production times leptonic branching fractions at ppbar collisions.

Summary of experimental measurements of the total inclusive cross sections for W and Z boson production times leptonic branching fractions at UA1 and UA2.

Summary of experimental measurements of the total inclusive cross sections for W and Z boson production times leptonic branching fractions at D0.

Summary of experimental measurements of the total inclusive cross sections for W and Z boson production times leptonic branching fractions at CDF.

Summary of theoretical predictions of the total inclusive cross sections for W boson production times leptonic branching fractions at pp collisions.

Summary of experimental measurements of the total inclusive cross sections for Z boson production times leptonic branching fractions at pp collisions.

Summary of experimental measurements of the total inclusive cross sections for W and Z boson production times leptonic branching fractions at CMS.

Observed two-dimensional likelihood scans for $c_{Hq}^{(1)}$ vs. $c_{Hq}^{(3)}$ while other coefficients are fixed at their SM values.

Observed two-dimensional likelihood scans for $c_{Hq}^{(1)}$ vs. $c_{Hu}$ while allowing the other coefficients to float freely at each point of the sca.

Observed two-dimensional likelihood scans for $c_{Hq}^{(1)}$ vs. $c_{Hu}$ while other coefficients are fixed at their SM values.

Observed two-dimensional likelihood scans for $c_{Hq}^{(1)}$ vs. $c_{Hd}$ while allowing the other coefficients to float freely at each point of the sca.

Observed two-dimensional likelihood scans for $c_{Hq}^{(1)}$ vs. $c_{Hd}$ while other coefficients are fixed at their SM values.

Observed two-dimensional likelihood scans for $c_{Hq}^{(3)}$ vs. $c_{Hu}$ while allowing the other coefficients to float freely at each point of the sca.

Observed two-dimensional likelihood scans for $c_{Hq}^{(3)}$ vs. $c_{Hu}$ while other coefficients are fixed at their SM values.

Observed two-dimensional likelihood scans for $c_{Hq}^{(3)}$ vs. $c_{Hd}$ while allowing the other coefficients to float freely at each point of the sca.

Observed two-dimensional likelihood scans for $c_{Hq}^{(3)}$ vs. $c_{Hd}$ while other coefficients are fixed at their SM values.

Observed two-dimensional likelihood scans for $c_{Hu}$ vs. $c_{Hd}$ while allowing the other coefficients to float freely at each point of the sca.

Observed two-dimensional likelihood scans for $c_{Hu}$ vs. $c_{Hd}$ while other coefficients are fixed at their SM values.

Observed two-dimensional likelihood scans for $c_{Hq}^{(1)}$ vs. $g^{(2)}_{ZZ}$ while allowing the other coefficients to float freely at each point of the sca.

Observed two-dimensional likelihood scans for $c_{Hq}^{(1)}$ vs. $g^{(2)}_{ZZ}$ while other coefficients are fixed at their SM values.

Observed two-dimensional likelihood scans for $c_{Hq}^{(3)}$ vs. $g^{(2)}_{ZZ}$ while allowing the other coefficients to float freely at each point of the sca.

Observed two-dimensional likelihood scans for $c_{Hq}^{(3)}$ vs. $g^{(2)}_{ZZ}$ while other coefficients are fixed at their SM values.

Observed two-dimensional likelihood scans for $c_{Hu}$ vs. $g^{(2)}_{ZZ}$ while allowing the other coefficients to float freely at each point of the sca.

Observed two-dimensional likelihood scans for $c_{Hu}$ vs. $g^{(2)}_{ZZ}$ while other coefficients are fixed at their SM values.

Observed two-dimensional likelihood scans for $c_{Hq}^{(1)}$ vs. $g^{(4)}_{ZZ}$ while allowing the other coefficients to float freely at each point of the sca.

Observed two-dimensional likelihood scans for $c_{Hq}^{(1)}$ vs. $g^{(4)}_{ZZ}$ while other coefficients are fixed at their SM values.

Observed two-dimensional likelihood scans for $c_{Hq}^{(3)}$ vs. $g^{(4)}_{ZZ}$ while allowing the other coefficients to float freely at each point of the sca.

Observed two-dimensional likelihood scans for $c_{Hq}^{(3)}$ vs. $g^{(4)}_{ZZ}$ while other coefficients are fixed at their SM values.

Observed two-dimensional likelihood scans for $c_{Hu}$ vs. $g^{(4)}_{ZZ}$ while allowing the other coefficients to float freely at each point of the sca.

Observed two-dimensional likelihood scans for $c_{Hu}$ vs. $g^{(4)}_{ZZ}$ while other coefficients are fixed at their SM values.

Observed two-dimensional likelihood scans for $c_{Hd}$ vs. $g^{(2)}_{ZZ}$ while allowing the other coefficients to float freely at each point of the sca.

Observed two-dimensional likelihood scans for $c_{Hd}$ vs. $g^{(2)}_{ZZ}$ while other coefficients are fixed at their SM values.

Observed two-dimensional likelihood scans for $c_{Hd}$ vs. $g^{(4)}_{ZZ}$ while allowing the other coefficients to float freely at each point of the sca.

Observed two-dimensional likelihood scans for $c_{Hd}$ vs. $g^{(4)}_{ZZ}$ while other coefficients are fixed at their SM values.

Observed two-dimensional likelihood scans for $g^{(2)}_{ZZ}$ vs. $g^{(4)}_{ZZ}$ while allowing the other coefficients to float freely at each point of the sca.

Observed two-dimensional likelihood scans for $g^{(2)}_{ZZ}$ vs. $g^{(4)}_{ZZ}$ while other coefficients are fixed at their SM values.

Results for the $P_3$ angular observable. The first uncertainties are statistical and the second systematic.

Results for the $P_4^\prime$ angular observable. The first uncertainties are statistical and the second systematic.

Results for the $P_5^\prime$ angular observable. The first uncertainties are statistical and the second systematic.

Results for the $P_6^\prime$ angular observable. The first uncertainties are statistical and the second systematic.

Results for the $P_8^\prime$ angular observable. The first uncertainties are statistical and the second systematic.

Results for the CP averaged observables $F_L$ and $P_1$–$P_8^\prime$. The first uncertainties are statistical and the second systematic.

Correlation matrix of angular observables Fl and P1–P8', considering only the statistical uncertainties from the maximum-likelihood fit in the interval 1.1 < $q^2$ < 2 GeV$^2$

Correlation matrix of angular observables Fl and P1–P8', considering only the statistical uncertainties from the maximum-likelihood fit in the interval 2 < $q^2$ < 4.3 GeV$^2$

Correlation matrix of angular observables Fl and P1–P8', considering only the statistical uncertainties from the maximum-likelihood fit in the interval 4.3 < $q^2$ < 6 GeV$^2$

Correlation matrix of angular observables Fl and P1–P8', considering only the statistical uncertainties from the maximum-likelihood fit in the interval 6 < $q^2$ < 8.68 GeV$^2$

Correlation matrix of angular observables Fl and P1–P8', considering only the statistical uncertainties from the maximum-likelihood fit in the interval 10.09 < $q^2$ < 12.86 GeV$^2$

Correlation matrix of angular observables Fl and P1–P8', considering only the statistical uncertainties from the maximum-likelihood fit in the interval 14.18 < $q^2$ < 16 GeV$^2$

Distributions of $S_{\mathrm{ML}}$ for data, simulated background and signal events with $n_{\mathrm{track}}$ of 4. The distributions are shown for split-SUSY signals with a gluino mass of 2000 GeV and neutralino mass of 1800 GeV. Different gluino proper decay lengths are shown as $c\tau$ in the legend. All distributions are normalized to unity.

Distributions of $S_{\mathrm{ML}}$ for data, simulated background and signal events with $n_{\mathrm{track}}$ of $\geq$ 5. The distributions are shown for split-SUSY signals with a gluino mass of 2000 GeV and neutralino mass of 1900 GeV. Different gluino proper decay lengths are shown as $c\tau$ in the legend. All distributions are normalized to unity.

Distributions of $S_{\mathrm{ML}}$ for data, simulated background and signal events with $n_{\mathrm{track}}$ of $\geq$ 5. The distributions are shown for split-SUSY signals with a gluino mass of 2000 GeV and neutralino mass of 1800 GeV. Different gluino proper decay lengths are shown as $c\tau$ in the legend. All distributions are normalized to unity.

The distribution of $n_{\mathrm{track}}$ in different $S_{\mathrm{ML}}$ regions for simulated background events. Events with 0 $ < S_{\mathrm{ML}} < $ 0.2 (blue), 0.2 $ < S_{\mathrm{ML}} < $ 0.6 (red), and 0.6 $ < S_{\mathrm{ML}} < $ 1.0 (green) are compared. All distributions are normalized to unity. The similar $n_{\mathrm{track}}$ distributions demonstrate that $n_{\mathrm{track}}$ and $S_{\mathrm{ML}}$ are decorrelated.

The distribution of $d_{\mathrm{BV}}$ in $K_{\mathrm{S}}^{\mathrm{0}}$ vertices in data (black) and simulation (purple). The lower panel shows the ratio between data and simulation.

The vertex reconstruction efficiency for artificially displaced vertices in data (black) and simulation (red). In this example, the artificially displaced vertices are corrected to mimic split-SUSY signal events with gluino mass of 2000 GeV and neutralino mass of 1800 GeV.

The ML tagging efficiency for artificially displaced vertices in data (black) and simulation (red). In this example, the artificially displaced vertices are corrected to mimic split-SUSY signal events with gluino mass of 2000 GeV and neutralino mass of 1800 GeV.

Number of predicted and observed events in the control, validation, and search regions. Predictions are calculated using Eqs. (2) and (3) and fitting the data under the background-only hypothesis. Regions are organized by $S_{\mathrm{ML}}$ and $n_{\mathrm{track}}$ values, and region names corresponding with Fig. 7 are given in parentheses. The predicted number of events that pass the $S_{\mathrm{ML}}$ selection and the observed number of events that pass or fail the $S_{\mathrm{ML}}$ selection are shown in seperate rows.

The 95% CL upper limit on the product of the cross section and branching fraction squared for the split-SUSY signal model with a mass splitting of 100 GeV, shown as a function of gluino mass and $c\tau$. The observed (solid black) and expected (dashed red) exclusion curves are overlaid on the limit plot.

The 95% CL upper limit on the product of the cross section and branching fraction squared for the split-SUSY signal model with a mass splitting of 100 GeV, shown as a function of gluino mass and $c\tau$. The observed (solid black) and expected (dashed red) exclusion curves are overlaid on the limit plot.

The 95% CL upper limit on the product of the cross section and branching fraction squared for the split-SUSY model with a $c\tau$ of 10 mm, shown as a function of gluino mass and mass splitting. The observed (solid black) and expected (dashed red) exclusion curves are overlaid on the limit plot.

The 95% CL upper limit on the product of the cross section and branching fraction squared for the split-SUSY model with a $c\tau$ of 10 mm, shown as a function of gluino mass and mass splitting. The observed (solid black) and expected (dashed red) exclusion curves are overlaid on the limit plot.

The 95% CL upper limit on the product of the cross section and branching fraction squared for the GMSB SUSY signal model, shown as a function of gluino mass and $c\tau$. The observed (solid black) and expected (dashed red) exclusion curves are overlaid on the limit plot.

The 95% CL upper limit on the product of the cross section and branching fraction squared for the GMSB SUSY signal model, shown as a function of gluino mass and $c\tau$. The observed (solid black) and expected (dashed red) exclusion curves are overlaid on the limit plot.

Efficiency of signal region A for all signal samples interperted in this search. The cross sections of all signal models are taken to be 1 fb in the table.

The cutflow table for $N1$ and $N2$ signal samples for different mass points. The samples consist of both $ee$ and $\mu\mu$ final states in equal ratio. Each sample is normalized to a total cross-section of 0.1 pb. The abbreviations used in the table are, SC = Same Charge, SF = Same Flavor.

The cutflow table for $N3$ signal samples for different mass points. The samples consist of both $ee$ and $\mu\mu$ final states in equal ratio. Each sample is normalized to a total cross-section of 10 pb. The abbreviations used in the table are, SS = Same Charge, SF = Same Flavor.

Expected and observed upper limits on the strength of HNL mixing with electron neutrinos

Expected and observed upper limits on the strength of HNL mixing with muon neutrinos

Expected and observed upper limits on the strength of HNL mixing with tau neutrinos

Distributions of the $\mathrm{NN^{SvsB}}$ output for data and the expected background after the likelihood fit in the $SR_{cQu --}$ signal region. The post-fit background expectations are shown as filled histograms, the combined pre-fit background expectations are shown as dashed lines. The signal distribution using the Wilson coefficient values $c_{tu}^{(1)}=0.04$, $c_{Qu}^{(1)}=0.1$, $c_{Qu}^{(8)}=0.1$ is shown with a dotted line, normalized to the same number of events as the background.

Summary of the observed and predicted number of events in the five 2$\ell$ control regions. The background prediction is shown after the combined likelihood fit to data under the signal-plus-background hypothesis across all control regions and signal regions. The uncertainties in the total yields are smaller than the sum in quadrature of the uncertainties in the individual contributions due to anti-correlations resulting from the likelihood fit. Int Conv refers to events with photon conversions in the inner detector, while Mat Conv refers to events with photon conversions in the detector material. QMisID refers to events with charge-misidentified lepton.

Summary of the observed and predicted number of events in the four 3$\ell$ control regions. The background prediction is shown after the combined likelihood fit to data under the signal-plus-background hypothesis across all control regions and signal regions. The uncertainties in the total yields are smaller than the sum in quadrature of the uncertainties in the individual contributions due to anti-correlations resulting from the likelihood fit. Int Conv refers to events with photon conversions in the inner detector, while Mat Conv refers to events with photon conversions in the detector material. QMisID refers to events with charge-misidentified lepton.

Comparison of the observed limits on the three WCs to the limit obtained by the ATLAS Run~1 analysis of events with $b$-tagged jets and a same-sign lepton pair [JHEP10(2015)150]. The cross-section limits from that analysis are converted to the limits on the WCs by using the fitted EFT parametrization. The bounds on the WCs are shown at the 68$\%$ CL (solid) and/or 95$\%$ CL (dashed) levels. For the ATLAS Run~1 analysis only the bounds at 95$\%$ CL are available. The vertical bar represents the SM prediction. The new-physics scale is set to $\Lambda = 1 \, \text{TeV}$.

Observed lower limits at 95$\%$ confidence level on the scale of new physics $\Lambda$ for Wilson coefficient values of 0.01, 1 and $4\pi^2$. The limits are compared with the limits obtained by the ATLAS Run~1 analysis of events with $b$-tagged jets and a same-sign lepton pair.

Observed and expected limits at 95$\%$ confidence level on the cross-section of $\sigma(pp \rightarrow tt)$.

2D likelihood scans of the Wilson coefficients $c_{tu}^{(1)}$ versus $c_{Qu}^{(1)}$. The new-physics scale is set to $\Lambda = 1 \, \text{TeV}$. In the table the $2\Delta$NLL values are shown which correspond to the minimized negative log likelihood for a given Wilson coefficient pair.

2D likelihood scans of the Wilson coefficients $c_{tu}^{(1)}$ versus $c_{Qu}^{(8)}$. The new-physics scale is set to $\Lambda = 1 \, \text{TeV}$. In the table the $2\Delta$NLL values are shown which correspond to the minimized negative log likelihood for a given Wilson coefficient pair.

2D likelihood scans of the Wilson coefficients $c_{Qu}^{(1)}$ versus $c_{Qu}^{(8)}$. The new-physics scale is set to $\Lambda = 1 \, \text{TeV}$. In the table the $2\Delta$NLL values are shown which correspond to the minimized negative log likelihood for a given Wilson coefficient pair.

1D-likelihood scan for $c_{tu}^{(1)}$. The observed and expected limits at 68\% and 95\% confidence level (CL) can be read off by finding the intersection of the likelihood scan curve with the 1$\sigma$ and 2$\sigma$ lines. The new-physics scale is set to $\Lambda = 1 \, \text{TeV}$. In the table the $2\Delta$NLL values are shown which correspond to the minimized negative log likelihood for a given Wilson coefficient value.

1D-likelihood scan for $c_{Qu}^{(1)}$. The observed and expected limits at 68\% and 95\% confidence level (CL) can be read off by finding the intersection of the likelihood scan curve with the 1$\sigma$ and 2$\sigma$ lines. The new-physics scale is set to $\Lambda = 1 \, \text{TeV}$. In the table the $2\Delta$NLL values are shown which correspond to the minimized negative log likelihood for a given Wilson coefficient value.

1D-likelihood scan for $c_{Qu}^{(8)}$. The observed and expected limits at 68\% and 95\% confidence level (CL) can be read off by finding the intersection of the likelihood scan curve with the 1$\sigma$ and 2$\sigma$ lines. The new-physics scale is set to $\Lambda = 1 \, \text{TeV}$. In the table the $2\Delta$NLL values are shown which correspond to the minimized negative log likelihood for a given Wilson coefficient value.

Expected (dashed line) value of $-2\ln\Lambda$ as a function of $\kappa_{\lambda}$ for the combined fit, when all other coupling modifiers are fixed to their SM predictions.

Observed (solid line) value of $-2\ln\Lambda$ as a function of $\kappa_{2V}$ for the combined fit, when all other coupling modifiers are fixed to their SM predictions.

Expected (dashed line) value of $-2\ln\Lambda$ as a function of $\kappa_{2V}$ for the combined fit, when all other coupling modifiers are fixed to their SM predictions.

Likelihood observed contours at 68% CL (solid line) and 95% CL (dashed line) in the ($\kappa_{\lambda}$, $\kappa_{2V}$) parameter space, when all other coupling modifiers are fixed to their SM predictions. The best-fit value, denoted by a cross in the plot, corresponds to the very first entry in the table.

Likelihood expected contours at 68% CL (teal shaded region) and 95% CL (yellow shaded region) in the ($\kappa_{\lambda}$, $\kappa_{2V}$) parameter space, when all other coupling modifiers are fixed to their SM predictions. In the plot, the SM prediction is indicated by the star.

Likelihood contours at 68% (solid line) and 95% (dashed line) CL in the ($c_{hhh}$, $c_{gghh}$) parameter space, when all other Wilson coefficients are fixed to their SM values. The corresponding expected contours are shown by the inner and outer shaded regions. The SM prediction is indicated by the star, while the observed best-fit value is denoted by the black cross. The best-fit value corresponds to the very first entry in the table.

Likelihood contours at 68% (solid line) and 95% (dashed line) CL in the ($c_{hhh}$, $c_{gghh}$) parameter space, when all other Wilson coefficients are fixed to their SM values. The corresponding expected contours are shown by the inner and outer shaded regions. The SM prediction is indicated by the star, while the observed best-fit value is denoted by the black cross.

Likelihood contours at 68% (solid line) and 95% (dashed line) CL in the ($c_{hhh}$, $c_{tthh}$) parameter space, when all other Wilson coefficients are fixed to their SM values. The corresponding expected contours are shown by the inner and outer shaded regions. The SM prediction is indicated by the star, while the observed best-fit value is denoted by the black cross. The best-fit value corresponds to the very first entry in the table.

Likelihood contours at 68% (solid line) and 95% (dashed line) CL in the ($c_{hhh}$, $c_{tthh}$) parameter space, when all other Wilson coefficients are fixed to their SM values. The corresponding expected contours are shown by the inner and outer shaded regions. The SM prediction is indicated by the star, while the observed best-fit value is denoted by the black cross.

Likelihood contours at 68% (solid line) and 95% (dashed line) CL in the ($c_{H}$, $c_{H◻}$) parameter space, when all other Wilson coefficients are fixed to their SM values. The corresponding expected contours are shown by the inner and outer shaded regions. The SM prediction is indicated by the star, while the observed best-fit value is denoted by the black cross. The best-fit value corresponds to the very first entry in the table.

Likelihood contours at 68% (solid line) and 95% (dashed line) CL in the ($c_{H}$, $c_{H◻}$) parameter space, when all other Wilson coefficients are fixed to their SM values. The corresponding expected contours are shown by the inner and outer shaded regions. The SM prediction is indicated by the star, while the observed best-fit value is denoted by the black cross.

The acceptance times efficiency [in %] for the signal ggF $HH$ process as a function of the coupling modifier $\kappa_{\lambda}$ in each analysis category for the $\tau_{\mathrm{had}}\tau_{\mathrm{had}}$ channel. The other coupling modifiers affecting ggF $HH$ production are fixed to their SM predictions.

The acceptance times efficiency [in %] for the signal ggF $HH$ process as a function of the coupling modifier $\kappa_{\lambda}$ in each analysis category for the $\tau_{\mathrm{lep}}\tau_{\mathrm{had}}$ SLT channel. The other coupling modifiers affecting ggF $HH$ production are fixed to their SM predictions.

The acceptance times efficiency [in %] for the signal ggF $HH$ process as a function of the coupling modifier $\kappa_{\lambda}$ in each analysis category for the $\tau_{\mathrm{lep}}\tau_{\mathrm{had}}$ LTT channel. The other coupling modifiers affecting ggF $HH$ production are fixed to their SM predictions.

The acceptance times efficiency [in %] for the signal VBF $HH$ process as a function of the coupling modifier $\kappa_{\lambda}$ in each analysis category for the $\tau_{\mathrm{had}}\tau_{\mathrm{had}}$ channel. The other coupling modifiers affecting VBF $HH$ production are fixed to their SM predictions.

The acceptance times efficiency [in %] for the signal VBF $HH$ process as a function of the coupling modifier $\kappa_{\lambda}$ in each analysis category for the $\tau_{\mathrm{lep}}\tau_{\mathrm{had}}$ SLT channel. The other coupling modifiers affecting VBF $HH$ production are fixed to their SM predictions.

The acceptance times efficiency [in %] for the signal VBF $HH$ process as a function of the coupling modifier $\kappa_{\lambda}$ in each analysis category for the $\tau_{\mathrm{lep}}\tau_{\mathrm{had}}$ LTT channel. The other coupling modifiers affecting VBF $HH$ production are fixed to their SM predictions.

The acceptance times efficiency [in %] for the signal VBF $HH$ process as a function of the coupling modifier $\kappa_{2V}$ in each analysis category for the $\tau_{\mathrm{had}}\tau_{\mathrm{had}}$ channel. The other coupling modifiers affecting VBF $HH$ production are fixed to their SM predictions.

The acceptance times efficiency [in %] for the signal VBF $HH$ process as a function of the coupling modifier $\kappa_{2V}$ in each analysis category for the $\tau_{\mathrm{lep}}\tau_{\mathrm{had}}$ SLT channel. The other coupling modifiers affecting VBF $HH$ production are fixed to their SM predictions.

The acceptance times efficiency [in %] for the signal VBF $HH$ process as a function of the coupling modifier $\kappa_{2V}$ in each analysis category for the $\tau_{\mathrm{lep}}\tau_{\mathrm{had}}$ LTT channel. The other coupling modifiers affecting VBF $HH$ production are fixed to their SM predictions.

Expected (dashed line) value of $-2\ln\Lambda$ as a function of $\kappa_{\lambda}$ for the fit to the $\tau_{\mathrm{had}}\tau_{\mathrm{had}}$ channel, when all other coupling modifiers are fixed to their SM predictions.

Expected (dashed line) value of $-2\ln\Lambda$ as a function of $\kappa_{\lambda}$ for the fit to the $\tau_{\mathrm{lep}}\tau_{\mathrm{had}}$ channel, when all other coupling modifiers are fixed to their SM predictions.

Observed (solid line) value of $-2\ln\Lambda$ as a function of $\kappa_{\lambda}$ for the fit to the $\tau_{\mathrm{had}}\tau_{\mathrm{had}}$ channel, when all other coupling modifiers are fixed to their SM predictions.

Observed (solid line) value of $-2\ln\Lambda$ as a function of $\kappa_{\lambda}$ for the fit to the $\tau_{\mathrm{lep}}\tau_{\mathrm{had}}$ channel, when all other coupling modifiers are fixed to their SM predictions.

Expected (dashed line) value of $-2\ln\Lambda$ as a function of $\kappa_{2V}$ for the fit to the $\tau_{\mathrm{had}}\tau_{\mathrm{had}}$ channel, when all other coupling modifiers are fixed to their SM predictions.

Expected (dashed line) value of $-2\ln\Lambda$ as a function of $\kappa_{2V}$ for the fit to the $\tau_{\mathrm{lep}}\tau_{\mathrm{had}}$ channel, when all other coupling modifiers are fixed to their SM predictions.

Observed (solid line) value of $-2\ln\Lambda$ as a function of $\kappa_{2V}$ for the fit to the $\tau_{\mathrm{had}}\tau_{\mathrm{had}}$ channel, when all other coupling modifiers are fixed to their SM predictions.

Observed (solid line) value of $-2\ln\Lambda$ as a function of $\kappa_{2V}$ for the fit to the $\tau_{\mathrm{lep}}\tau_{\mathrm{had}}$ channel, when all other coupling modifiers are fixed to their SM predictions.

Summary of the expected (blue) and observed (orange) one-dimensional constraints at 68% (solid error bars) and 95% (dashed error bars) CL on the coefficients of selected HEFT (top) and SMEFT (bottom) operators describing Higgs boson interactions affecting $HH$ production through gluon-gluon fusion. The results are obtained from one dimensional scans of the profile log-likelihood assuming that all other Wilson coefficients are fixed to their SM values. The contribution from VBF $HH$ production is neglected.