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.

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

The model-independent 95\% \CL expected and observed upper limits set on ${\sigma(\PP\to 2\Pa+\PX)\mathcal{B}^2(\Pa\to 2\PGm)\alphaGen}$ over the range $0.21 < \MPa < 9\GeV$ for the combined 2016, 2017, and 2018 analyses. Mass ranges that overlap with \JPsi and \PgU resonances are excluded from the search

The 95\% \CL observed upper limits on the effective coupling $\CahLambda$ of the ALP to the SM Higgs assuming the branching fraction of the ALP to muons is 1 (blue) and 0.1 (orange), for both the expected (dashed) and observed (solid) limits

The 95\% \CL observed upper limits on the effective coupling $\cll/\Lambda$ of the ALP to the SM leptons. The orange shaded region represents the parameter space excluded by this search under three choices of $\CahLambda$ $1\TeV^{-2}$ (solid), $0.1\TeV^{-2}$ (dashed), and $0.01\TeV^{-2}$ (dotted)

The 95\% \CL observed upper limits on $\varepsilon^{2} \mathcal{B}(\ZDark \to \sDark \overline{\sDark}) \mathcal{B}^{2}(\sDark \rightarrow 2\PGm)$ for the vector portal model as a function of the dark scalar mass $\MsDark$ and dark vector boson mass $\MZDark$. The yellow and green bands indicate one and two standard deviation values of the limit, respectively. Because the model-independent limits are calculated only up to a dimuon mass of 60\GeV, the \sDark mass considered for these limits is below 60\GeV

The 95\% \CL observed upper limits for $\sigma(\PP\to\PhOneTwo\to 2\PaOne) \mathcal{B}^2(\PaOne\to 2\PGm)$ for the NMSSM as a function of $\MPaOne$ for $\MPhOneTwo$ = 90~GEV. The data presented here reflect the results of the 2017 and 2018 datasets combined with the previously published results of the 2016 dataset in \cite{CMS:2018jid}

The 95\% \CL observed upper limits for $\sigma(\PP\to\PhOneTwo\to 2\PaOne) \mathcal{B}^2(\PaOne\to 2\PGm)$ for the NMSSM as a function of $\MPaOne$ for $\MPhOneTwo$ = 125~GEV. The data presented here reflect the results of the 2017 and 2018 datasets combined with the previously published results of the 2016 dataset in \cite{CMS:2018jid}

The 95\% \CL observed upper limits for $\sigma(\PP\to\PhOneTwo\to 2\PaOne) \mathcal{B}^2(\PaOne\to 2\PGm)$ for the NMSSM as a function of $\MPaOne$ for $\MPhOneTwo$ = 150~GEV. The data presented here reflect the results of the 2017 and 2018 datasets combined with the previously published results of the 2016 dataset in \cite{CMS:2018jid}

The 95\% \CL observed upper limits (black solid curves) from this search as interpreted in the dark SUSY scenario for the process $\PP \to \Ph \to 2\nOne \to 2\gammaDark + 2\nDark \to 4\PGm + \PX$, with \mbox{$\MnOne = 60\GeV$} and $\MnDark = 1\GeV$. The limits are presented in the plane of $\varepsilon$ and $\MgammaDark$. The color gradient represents different branching fraction assumptions for \mbox{$\mathcal{B}(\Ph \to 2\nOne \to 2\gammaDark + 2\nDark)$}, ranging from dark orange (0.05\%) to light orange (10\%). The degradation of the limit around 1\GeV is attributed to the drop in the dimuon branching fraction $\mathcal{B}(\gammaDark \to 2\mu)$ due to the dimuon resonance of the $\phi$ meson~\cite{Batell:2009yf}

Transverse momentum spectrum of charged hadrons measured at midrapidity ($|y|<0.5$) in INEL>0 pp collisions at $\sqrt{s}$ = 13 TeV for different flattenicity event classes selected with the V0M estimator at forward rapidity (top figure, upper panel)

Transverse momentum spectrum of $\pi^{+} + \pi^{-}$ measured at midrapidity ($|y|<0.5$) in INEL>0 pp collisions at $\sqrt{s}$ = 13 TeV for different flattenicity event classes and HM events (0--1% V0M) in the same flattenicity event classes (bottom figure, upper panel)

Transverse momentum spectrum of $K^{+} + K^{-}$ measured at midrapidity ($|y|<0.5$) in INEL>0 pp collisions at $\sqrt{s}$ = 13 TeV for different flattenicity event classes and HM events (0--1% V0M) in the same flattenicity event classes (bottom figure, upper panel)

Transverse momentum spectrum of $p + \overline{p}$ measured at midrapidity ($|y|<0.5$) in INEL>0 pp collisions at $\sqrt{s}$ = 13 TeV for different flattenicity event classes and HM events (0--1% V0M) in the same flattenicity event classes (bottom figure, upper panel)

Transverse momentum spectrum of charged hadrons measured at midrapidity ($|y|<0.5$) in INEL>0 pp collisions at $\sqrt{s}$ = 13 TeV for different flattenicity event classes and HM events (0--1% V0M) in the same flattenicity event classes (bottom figure, upper panel)

The $p_{\rm T}$-differential $(K^{+} + K^{-})/(\pi^{+}+\pi^{-})$ particle ratio for two extremes of flattenicity event classes measured in multiplicity-integrated events in pp collisions at $\sqrt{s}$ = 13 TeV

The $p_{\rm T}$-differential $(p + \overline{p})$/($\pi^{+}+\pi^{-}$) particle ratio for two extremes of flattenicity event classes measured in multiplicity-integrated events in pp collisions at $\sqrt{s}$ = 13 TeV

The $p_{\rm T}$-differential $(K^{+} + K^{-})/(\pi^{+}+\pi^{-})$ particle ratio for two extremes of flattenicity event classes measured in the 0-1% V0M event class in pp collisions at $\sqrt{s}$ = 13 TeV

The $p_{\rm T}$-differential $(p + \overline{p})$/($\pi^{+}+\pi^{-}$) particle ratio for two extremes of flattenicity event classes measured in the 0-1% V0M event class in pp collisions at $\sqrt{s}$ = 13 TeV

Transverse momentum-integrated $(K^{+} + K^{-})/(\pi^{+}+\pi^{-})$ particle ratio as a function of the average charged-particle density calculated for different flattenicity event classes in pp collisions at $\sqrt{s}$ = 13 TeV

Transverse momentum-integrated $(p + \overline{p})$/($\pi^{+}+\pi^{-}$) particle ratio as a function of the average charged-particle density calculated for different flattenicity event classes in pp collisions at $\sqrt{s}$ = 13 TeV

Average transverse momentum of $\pi^{+}+\pi^{-}$ as a function of the average charged-particle density calculated for different flattenicity event classes in pp collisions at $\sqrt{s}$ = 13 TeV

Average transverse momentum of $K^{+} + K^{-}$ as a function of the average charged-particle density calculated for different flattenicity event classes in pp collisions at $\sqrt{s}$ = 13 TeV

Average transverse momentum of $p + \overline{p}$ as a function of the average charged-particle density calculated for different flattenicity event classes in pp collisions at $\sqrt{s}$ = 13 TeV

$\chi_{c0}(4500) \to J/\psi(\to \mu^+ \mu^-)\phi(\to K^+ K^-)$ diffractive production cross-section, multiplied by $\mathcal{B}(J/\psi \to \mu^+ \mu^-)$ and $\mathcal{B}(\phi \to K^+ K^-)$ branching ratios.

$\chi_{c1}(4685) \ \mathrm{or} \ \chi_{c0}(4700) \to J/\psi(\to \mu^+ \mu^-)\phi(\to K^+ K^-)$ diffractive production cross-section, multiplied by $\mathcal{B}(J/\psi \to \mu^+ \mu^-)$ and $\mathcal{B}(\phi \to K^+ K^-)$ branching ratios.

Nonresonant $J/\psi(\to \mu^+ \mu^-)\phi(\to K^+ K^-)$ diffractive production cross-section, multiplied by $\mathcal{B}(J/\psi \to \mu^+ \mu^-)$ and $\mathcal{B}(\phi \to K^+ K^-)$ branching ratios.

Slices of 2D distributions of observed events and the post-fit templates in the TT pass region, projected onto the plane of leading jet mass mJ1, including expected radion signal at 1.5 TeV.

Slices of 2D distributions of observed events and the post-fit templates in the TT pass region, projected onto the plane of leading jet mass mJ1, including expected radion signal at 1.5 TeV.

Slices of 2D distributions of observed events and the post-fit templates in the TT pass region, projected onto the plane of leading jet mass mJ1, including expected radion signal at 1.5 TeV.

Slices of 2D distributions of observed events and the post-fit templates in the LL pass region, projected onto the plane of corrected HH mass (mHH) mHH, including expected radion signal at 1.5 TeV.

Slices of 2D distributions of observed events and the post-fit templates in the LL pass region, projected onto the plane of corrected HH mass (mHH) mHH, including expected radion signal at 1.5 TeV.

Slices of 2D distributions of observed events and the post-fit templates in the LL pass region, projected onto the plane of corrected HH mass (mHH) mHH, including expected radion signal at 1.5 TeV.

Slices of 2D distributions of observed events and the post-fit templates in the TT pass region, projected onto the plane of corrected HH mass (mHH) mHH, including expected radion signal at 1.5 TeV.

Slices of 2D distributions of observed events and the post-fit templates in the TT pass region, projected onto the plane of corrected HH mass (mHH) mHH, including expected radion signal at 1.5 TeV.

Slices of 2D distributions of observed events and the post-fit templates in the TT pass region, projected onto the plane of corrected HH mass (mHH) mHH, including expected radion signal at 1.5 TeV.

Slices of 2D distributions of observed events and the post-fit templates in the semi-resolved pass region, projected onto the plane of corrected HH mass (mHH) mHH, including expected radion signal at 1.5 TeV.

Slices of 2D distributions of observed events and the post-fit templates in the semi-resolved pass region, projected onto the plane of corrected HH mass (mHH) mHH, including expected radion signal at 1.5 TeV.

Slices of 2D distributions of observed events and the post-fit templates in the semi-resolved pass region, projected onto the plane of corrected HH mass (mHH) mHH, including expected radion signal at 1.5 TeV.

Slices of 2D distributions of observed events and the post-fit templates in the semi-resolved pass region, projected onto the plane of leading jet mass mJ1, including expected radion signal at 1.5 TeV.

Slices of 2D distributions of observed events and the post-fit templates in the semi-resolved pass region, projected onto the plane of leading jet mass mJ1, including expected radion signal at 1.5 TeV.

Slices of 2D distributions of observed events and the post-fit templates in the semi-resolved pass region, projected onto the plane of leading jet mass mJ1, including expected radion signal at 1.5 TeV.

The observed (solid black line) and expected (dashed black line) upper limits at 95% CL on $\sigma$(pp->X)B(X->HH->4b) for the narrow spin-0 radion model. The green (yellow) bands represent one (two) standard deviations from the expected limit. The predicted theoretical cross sections are also shown.

The observed (solid black line) and expected (dashed black line) upper limits at 95% CL on $\sigma$(pp->X)B(X->HH->4b) for the narrow width spin-2 bulk graviton model. The green (yellow) bands represent one (two) standard deviations from the expected limit. The predicted theoretical cross sections are also shown.

The observed (solid black line) and expected (dashed black line) upper limits at 95% CL on $\sigma$(pp->X)B(X->HH->4b) for the 10%-width spin-0 radion model. The green (yellow) bands represent one (two) standard deviations from the expected limit. The predicted theoretical cross sections are also shown.

The observed (solid black line) and expected (dashed black line) upper limits at 95% CL on $\sigma$(pp->X)B(X->HH->4b) for the 10%-width spin-2 bulk graviton model. The green (yellow) bands represent one (two) standard deviations from the expected limit. The predicted theoretical cross sections are also shown.