CMS logoCMS event Hgg
Compact Muon Solenoid
LHC, CERN

CMS-B2G-23-007 ; CERN-EP-2026-001
Search for heavy scalar resonances decaying to Lorentz-boosted Higgs and Higgs-like bosons in the $ \mathrm{b}\overline{\mathrm{b}}4\mathrm{q} $ final state at $ \sqrt{s} = $ 13 TeV
Submitted to the Journal of High Energy Physics
Abstract: A search is performed for a heavy scalar resonance ($\mathrm{X}$) decaying to a Higgs boson ($\mathrm{H}$) and a Higgs-like scalar boson ($\mathrm{Y}$) in the two bottom quark ($ \mathrm{H}\to\mathrm{b}\overline{\mathrm{b}} $) and four quark ($ \mathrm{Y} \to\mathrm{V}\mathrm{V}\to4\mathrm{q} $) final state, where $\mathrm{V}$ denotes a $\mathrm{W}$ or $\mathrm{Z}$ boson. Masses of the $\mathrm{X}$ between 900 and 4000 GeV and the $\mathrm{Y}$ between 60 and 2800 GeV are considered. The search is performed in data collected by the CMS experiment at the CERN LHC from proton-proton collisions at 13 TeV center-of-mass energy, with a data set corresponding to a total integrated luminosity of 138 fb$ ^{-1} $. It targets the Lorentz-boosted regime, in which the products of the $ \mathrm{H}\to\mathrm{b}\overline{\mathrm{b}} $ decay can be reconstructed as a single large-area jet, and those from the $ \mathrm{Y} \to\mathrm{V}\mathrm{V}\to4\mathrm{q} $ decay as either one $ \mathrm{Y} \to4\mathrm{q} $ or two $ \mathrm{V}\to\mathrm{q}\overline{\mathrm{q}} $ jets. Jet identification and mass reconstruction exploit machine-learning tools, including a novel attention-based ``particle transformer'' for $ \mathrm{Y} \to4\mathrm{q} $ identification. No significant excess is observed in the data above the standard model background expectation. Upper limits on the product of production cross section and branching fraction as low as 0.2fb$^{-1}$ are derived at 95% confidence level for various mass hypotheses. This is the first search at the LHC for scalar resonances in the all-hadronic $ \mathrm{b}\overline{\mathrm{b}}\mathrm{V}\mathrm{V} $ decay channel.
Figures & Tables Summary References CMS Publications
Figures

png pdf
Figure 1:
Schematic diagrams of the $ \mathrm{X} \to \mathrm{HY} $ decay in the $ \mathrm{H}\to\mathrm{b}\overline{\mathrm{b}} $ and $ {{Y}} \to\mathrm{V}\mathrm{V}\to4\mathrm{q} $ all-hadronic final state, with $ \mathrm{V}=\mathrm{W} $ or $\mathrm{Z} $. Left: The fully merged regime ($ {m_{{Y}} } \lesssim 0.1 {m_{{X}} } $) in which {Y} is highly boosted, with its decay products merged and reconstructed as a single large-area jet. Right: The semimerged regime (0.1 $ {m_{{X}} } \lesssim {m_{{Y}} } < {m_{{X}} } - {m_\mathrm{H}} $) in which {Y} is not as boosted, and each daughter V boson is instead reconstructed individually as a large-area jet.

png pdf
Figure 2:
Schematic of the event categorization into the signal pass (SP), signal fail (SF), validation pass (VP), and validation fail (VF) regions for the fully merged (left) and semimerged (right) categories. The pass and fail tagging regions of the fully merged and semimerged categories are determined by application of the high purity (HP) working point (WP) of the $ T_\text{YWW} $ and $ T_\text{Xbb} $ taggers, respectively.

png pdf
Figure 2-a:
Schematic of the event categorization into the signal pass (SP), signal fail (SF), validation pass (VP), and validation fail (VF) regions for the fully merged (left) and semimerged (right) categories. The pass and fail tagging regions of the fully merged and semimerged categories are determined by application of the high purity (HP) working point (WP) of the $ T_\text{YWW} $ and $ T_\text{Xbb} $ taggers, respectively.

png pdf
Figure 2-b:
Schematic of the event categorization into the signal pass (SP), signal fail (SF), validation pass (VP), and validation fail (VF) regions for the fully merged (left) and semimerged (right) categories. The pass and fail tagging regions of the fully merged and semimerged categories are determined by application of the high purity (HP) working point (WP) of the $ T_\text{YWW} $ and $ T_\text{Xbb} $ taggers, respectively.

png pdf
Figure 3:
Measured transfer functions $ R^\text{sim}({{m_{{X}} },{m_{{Y}} }}) $ (left) and $ R({{m_{{X}} },{m_{{Y}} }}) $ (right) between the signal pass (SP) and signal fail (SF) regions of the fully merged category used to estimate the QCD multijet background contribution.

png pdf
Figure 3-a:
Measured transfer functions $ R^\text{sim}({{m_{{X}} },{m_{{Y}} }}) $ (left) and $ R({{m_{{X}} },{m_{{Y}} }}) $ (right) between the signal pass (SP) and signal fail (SF) regions of the fully merged category used to estimate the QCD multijet background contribution.

png pdf
Figure 3-b:
Measured transfer functions $ R^\text{sim}({{m_{{X}} },{m_{{Y}} }}) $ (left) and $ R({{m_{{X}} },{m_{{Y}} }}) $ (right) between the signal pass (SP) and signal fail (SF) regions of the fully merged category used to estimate the QCD multijet background contribution.

png pdf
Figure 4:
Projected distributions of the $ {m_{{X}} ^\text{rec}} $ (left) and $ {m_{{Y}} ^\text{rec}} $ (right) observables in the validation pass (VP) region of the fully merged (FM) category, after a background-only maximum likelihood fit to the data. The upper panels show the data and fitted background estimates, as well as the $ \pm $1 standard deviation uncertainty in the total background estimate ($ {\sigma_\text{syst}} $) per bin, divided by the respective bin widths, in gray. The vertical error bars for the data represent the statistical uncertainty, while the horizontal bars represent the bin widths. The lower panels show the pull per bin, defined as $\frac{ (\mathrm{data} - \mathrm{Bkg.})}{{\sigma_\text{stat}}} $, as well as a gray band indicating the systematic uncertainty in the pulls, estimated by dividing the corresponding systematic uncertainty in the main panel by $ {\sigma_\text{stat}} $.

png pdf
Figure 4-a:
Projected distributions of the $ {m_{{X}} ^\text{rec}} $ (left) and $ {m_{{Y}} ^\text{rec}} $ (right) observables in the validation pass (VP) region of the fully merged (FM) category, after a background-only maximum likelihood fit to the data. The upper panels show the data and fitted background estimates, as well as the $ \pm $1 standard deviation uncertainty in the total background estimate ($ {\sigma_\text{syst}} $) per bin, divided by the respective bin widths, in gray. The vertical error bars for the data represent the statistical uncertainty, while the horizontal bars represent the bin widths. The lower panels show the pull per bin, defined as $\frac{ (\mathrm{data} - \mathrm{Bkg.})}{{\sigma_\text{stat}}} $, as well as a gray band indicating the systematic uncertainty in the pulls, estimated by dividing the corresponding systematic uncertainty in the main panel by $ {\sigma_\text{stat}} $.

png pdf
Figure 4-b:
Projected distributions of the $ {m_{{X}} ^\text{rec}} $ (left) and $ {m_{{Y}} ^\text{rec}} $ (right) observables in the validation pass (VP) region of the fully merged (FM) category, after a background-only maximum likelihood fit to the data. The upper panels show the data and fitted background estimates, as well as the $ \pm $1 standard deviation uncertainty in the total background estimate ($ {\sigma_\text{syst}} $) per bin, divided by the respective bin widths, in gray. The vertical error bars for the data represent the statistical uncertainty, while the horizontal bars represent the bin widths. The lower panels show the pull per bin, defined as $\frac{ (\mathrm{data} - \mathrm{Bkg.})}{{\sigma_\text{stat}}} $, as well as a gray band indicating the systematic uncertainty in the pulls, estimated by dividing the corresponding systematic uncertainty in the main panel by $ {\sigma_\text{stat}} $.

png pdf
Figure 5:
Projected distributions of the $ {m_{{X}} ^\text{rec}} $ (left) and $ {m_{{Y}} ^\text{rec}} $ (right) observables in the validation pass (VP) region of the semimerged (SM) category, after a background-only maximum likelihood fit to the data. The upper panels show the data and fitted background estimates, as well as the $ \pm $1 standard deviation uncertainty in the total background estimate ($ {\sigma_\text{syst}} $) per bin, divided by the respective bin widths, in gray. The vertical error bars for the data represent the statistical uncertainty, while the horizontal bars represent the bin widths. The lower panels show the pull per bin, defined as $\frac{ (\mathrm{data} - \mathrm{Bkg.})}{{\sigma_\text{stat}}} $ as well as a gray band indicating the systematic uncertainty in the pulls, estimated by dividing the corresponding systematic uncertainty in the main panel by $ {\sigma_\text{stat}} $.

png pdf
Figure 5-a:
Projected distributions of the $ {m_{{X}} ^\text{rec}} $ (left) and $ {m_{{Y}} ^\text{rec}} $ (right) observables in the validation pass (VP) region of the semimerged (SM) category, after a background-only maximum likelihood fit to the data. The upper panels show the data and fitted background estimates, as well as the $ \pm $1 standard deviation uncertainty in the total background estimate ($ {\sigma_\text{syst}} $) per bin, divided by the respective bin widths, in gray. The vertical error bars for the data represent the statistical uncertainty, while the horizontal bars represent the bin widths. The lower panels show the pull per bin, defined as $\frac{ (\mathrm{data} - \mathrm{Bkg.})}{{\sigma_\text{stat}}} $ as well as a gray band indicating the systematic uncertainty in the pulls, estimated by dividing the corresponding systematic uncertainty in the main panel by $ {\sigma_\text{stat}} $.

png pdf
Figure 5-b:
Projected distributions of the $ {m_{{X}} ^\text{rec}} $ (left) and $ {m_{{Y}} ^\text{rec}} $ (right) observables in the validation pass (VP) region of the semimerged (SM) category, after a background-only maximum likelihood fit to the data. The upper panels show the data and fitted background estimates, as well as the $ \pm $1 standard deviation uncertainty in the total background estimate ($ {\sigma_\text{syst}} $) per bin, divided by the respective bin widths, in gray. The vertical error bars for the data represent the statistical uncertainty, while the horizontal bars represent the bin widths. The lower panels show the pull per bin, defined as $\frac{ (\mathrm{data} - \mathrm{Bkg.})}{{\sigma_\text{stat}}} $ as well as a gray band indicating the systematic uncertainty in the pulls, estimated by dividing the corresponding systematic uncertainty in the main panel by $ {\sigma_\text{stat}} $.

png pdf
Figure 6:
Projected distributions of the $ {m_{{X}} ^\text{rec}} $ (left) and $ {m_{{Y}} ^\text{rec}} $ (right) observables in the signal pass (SP) region of the fully merged (FM) category, after a maximum likelihood fit to the data with an $ {{{X}} \to\mathrm{H}{{Y}} \to\mathrm{b}\overline{\mathrm{b}}\mathrm{W}\mathrm{W}} $ signal with $ {m_{{X}} } = $ 900 GeV and $ {m_{{Y}} } = $ 80 GeV and product of best-fit production cross section and branching fraction of 15fb$^{-1}$. The signal distribution has been scaled by a factor of two in the upper panels to help distinguish it from the background distributions. The lower panels show the pull per bin, defined as $\frac{ (\mathrm{data} - \mathrm{Bkg.})}{{\sigma_\text{stat}}} $, as well as a gray band indicating the systematic uncertainty in the pulls, estimated by dividing the corresponding systematic uncertainty in the main panel by $ {\sigma_\text{stat}} $.

png pdf
Figure 6-a:
Projected distributions of the $ {m_{{X}} ^\text{rec}} $ (left) and $ {m_{{Y}} ^\text{rec}} $ (right) observables in the signal pass (SP) region of the fully merged (FM) category, after a maximum likelihood fit to the data with an $ {{{X}} \to\mathrm{H}{{Y}} \to\mathrm{b}\overline{\mathrm{b}}\mathrm{W}\mathrm{W}} $ signal with $ {m_{{X}} } = $ 900 GeV and $ {m_{{Y}} } = $ 80 GeV and product of best-fit production cross section and branching fraction of 15fb$^{-1}$. The signal distribution has been scaled by a factor of two in the upper panels to help distinguish it from the background distributions. The lower panels show the pull per bin, defined as $\frac{ (\mathrm{data} - \mathrm{Bkg.})}{{\sigma_\text{stat}}} $, as well as a gray band indicating the systematic uncertainty in the pulls, estimated by dividing the corresponding systematic uncertainty in the main panel by $ {\sigma_\text{stat}} $.

png pdf
Figure 6-b:
Projected distributions of the $ {m_{{X}} ^\text{rec}} $ (left) and $ {m_{{Y}} ^\text{rec}} $ (right) observables in the signal pass (SP) region of the fully merged (FM) category, after a maximum likelihood fit to the data with an $ {{{X}} \to\mathrm{H}{{Y}} \to\mathrm{b}\overline{\mathrm{b}}\mathrm{W}\mathrm{W}} $ signal with $ {m_{{X}} } = $ 900 GeV and $ {m_{{Y}} } = $ 80 GeV and product of best-fit production cross section and branching fraction of 15fb$^{-1}$. The signal distribution has been scaled by a factor of two in the upper panels to help distinguish it from the background distributions. The lower panels show the pull per bin, defined as $\frac{ (\mathrm{data} - \mathrm{Bkg.})}{{\sigma_\text{stat}}} $, as well as a gray band indicating the systematic uncertainty in the pulls, estimated by dividing the corresponding systematic uncertainty in the main panel by $ {\sigma_\text{stat}} $.

png pdf
Figure 7:
Projected distributions of the $ {m_{{X}} ^\text{rec}} $ (left) and $ {m_{{Y}} ^\text{rec}} $ (right) observables in the signal pass (SP) region of the semimerged (SM) category, after a maximum likelihood fit to the data with an $ {{{X}} \to\mathrm{H}{{Y}} \to\mathrm{b}\overline{\mathrm{b}}\mathrm{W}\mathrm{W}} $ signal with $ {m_{{X}} } = $ 1200 GeV and $ {m_{{Y}} } = $ 900 GeV and product of best-fit production cross section and branching fraction of 250fb$^{-1}$. The lower panels show the pull per bin, defined as $\frac{ (\mathrm{data} - \mathrm{Bkg.})}{{\sigma_\text{stat}}} $, as well as a gray band indicating the systematic uncertainty in the pulls, estimated by dividing the corresponding systematic uncertainty in the main panel by $ {\sigma_\text{stat}} $.

png pdf
Figure 7-a:
Projected distributions of the $ {m_{{X}} ^\text{rec}} $ (left) and $ {m_{{Y}} ^\text{rec}} $ (right) observables in the signal pass (SP) region of the semimerged (SM) category, after a maximum likelihood fit to the data with an $ {{{X}} \to\mathrm{H}{{Y}} \to\mathrm{b}\overline{\mathrm{b}}\mathrm{W}\mathrm{W}} $ signal with $ {m_{{X}} } = $ 1200 GeV and $ {m_{{Y}} } = $ 900 GeV and product of best-fit production cross section and branching fraction of 250fb$^{-1}$. The lower panels show the pull per bin, defined as $\frac{ (\mathrm{data} - \mathrm{Bkg.})}{{\sigma_\text{stat}}} $, as well as a gray band indicating the systematic uncertainty in the pulls, estimated by dividing the corresponding systematic uncertainty in the main panel by $ {\sigma_\text{stat}} $.

png pdf
Figure 7-b:
Projected distributions of the $ {m_{{X}} ^\text{rec}} $ (left) and $ {m_{{Y}} ^\text{rec}} $ (right) observables in the signal pass (SP) region of the semimerged (SM) category, after a maximum likelihood fit to the data with an $ {{{X}} \to\mathrm{H}{{Y}} \to\mathrm{b}\overline{\mathrm{b}}\mathrm{W}\mathrm{W}} $ signal with $ {m_{{X}} } = $ 1200 GeV and $ {m_{{Y}} } = $ 900 GeV and product of best-fit production cross section and branching fraction of 250fb$^{-1}$. The lower panels show the pull per bin, defined as $\frac{ (\mathrm{data} - \mathrm{Bkg.})}{{\sigma_\text{stat}}} $, as well as a gray band indicating the systematic uncertainty in the pulls, estimated by dividing the corresponding systematic uncertainty in the main panel by $ {\sigma_\text{stat}} $.

png pdf
Figure 8:
The 2D pulls in the $ {m_{{X}} ^\text{rec}} $ and $ {m_{{Y}} ^\text{rec}} $ observables in the signal pass (SP) regions of the fully merged (FM, left) and semimerged (SM, right) categories, after a maximum likelihood fit to the data with an $ {{{X}} \to\mathrm{H}{{Y}} \to\mathrm{b}\overline{\mathrm{b}}\mathrm{W}\mathrm{W}} $ signal with $ {m_{{X}} } = $ 900 GeV and $ {m_{{Y}} } = $ 80 GeV, and $ {m_{{X}} } = $ 1200 GeV and $ {m_{{Y}} } = $ 900 GeV, respectively. Pulls are defined as $\frac{(\mathrm{data} - \mathrm{Bkg.})}{\sigma_\text{stat}} $ where $\sigma$ is the sum in quadrature of the background and data uncertainties.

png pdf
Figure 8-a:
The 2D pulls in the $ {m_{{X}} ^\text{rec}} $ and $ {m_{{Y}} ^\text{rec}} $ observables in the signal pass (SP) regions of the fully merged (FM, left) and semimerged (SM, right) categories, after a maximum likelihood fit to the data with an $ {{{X}} \to\mathrm{H}{{Y}} \to\mathrm{b}\overline{\mathrm{b}}\mathrm{W}\mathrm{W}} $ signal with $ {m_{{X}} } = $ 900 GeV and $ {m_{{Y}} } = $ 80 GeV, and $ {m_{{X}} } = $ 1200 GeV and $ {m_{{Y}} } = $ 900 GeV, respectively. Pulls are defined as $\frac{(\mathrm{data} - \mathrm{Bkg.})}{\sigma_\text{stat}} $ where $\sigma$ is the sum in quadrature of the background and data uncertainties.

png pdf
Figure 8-b:
The 2D pulls in the $ {m_{{X}} ^\text{rec}} $ and $ {m_{{Y}} ^\text{rec}} $ observables in the signal pass (SP) regions of the fully merged (FM, left) and semimerged (SM, right) categories, after a maximum likelihood fit to the data with an $ {{{X}} \to\mathrm{H}{{Y}} \to\mathrm{b}\overline{\mathrm{b}}\mathrm{W}\mathrm{W}} $ signal with $ {m_{{X}} } = $ 900 GeV and $ {m_{{Y}} } = $ 80 GeV, and $ {m_{{X}} } = $ 1200 GeV and $ {m_{{Y}} } = $ 900 GeV, respectively. Pulls are defined as $\frac{(\mathrm{data} - \mathrm{Bkg.})}{\sigma_\text{stat}} $ where $\sigma$ is the sum in quadrature of the background and data uncertainties.

png pdf
Figure 9:
Median expected (left) and observed (right) upper limits at 95% CL on the product of production cross section and branching fraction for $ {{{X}} \to\mathrm{H}{{Y}} \to\mathrm{b}\overline{\mathrm{b}}\mathrm{W}\mathrm{W}} $ on the 2D $ {m_{{X}} } $ and $ {m_{{Y}} } $ plane. The dashed line separates the mass points considered in the fully merged (FM) and semimerged (SM) categories.

png pdf
Figure 9-a:
Median expected (left) and observed (right) upper limits at 95% CL on the product of production cross section and branching fraction for $ {{{X}} \to\mathrm{H}{{Y}} \to\mathrm{b}\overline{\mathrm{b}}\mathrm{W}\mathrm{W}} $ on the 2D $ {m_{{X}} } $ and $ {m_{{Y}} } $ plane. The dashed line separates the mass points considered in the fully merged (FM) and semimerged (SM) categories.

png pdf
Figure 9-b:
Median expected (left) and observed (right) upper limits at 95% CL on the product of production cross section and branching fraction for $ {{{X}} \to\mathrm{H}{{Y}} \to\mathrm{b}\overline{\mathrm{b}}\mathrm{W}\mathrm{W}} $ on the 2D $ {m_{{X}} } $ and $ {m_{{Y}} } $ plane. The dashed line separates the mass points considered in the fully merged (FM) and semimerged (SM) categories.

png pdf
Figure 10:
Expected and observed upper limits at 95% CL on the product of production cross section and branching fraction for $ {{{X}} \to\mathrm{H}{{Y}} \to\mathrm{b}\overline{\mathrm{b}}\mathrm{W}\mathrm{W}} $ signals. The dashed and solid black lines represent median expected and observed limits, respectively. The inner (green) and outer (yellow) bands represent the $ \pm $ 1 and $ \pm $ 2 standard deviations for the expected limits, respectively. The left and right set of points in each panel correspond to those considered in the fully merged (FM) and semimerged (SM) categories, respectively.
Tables

png pdf
Table 1:
Sources of systematic uncertainties and their effects on different processes in the fully merged (FM) and semimerged (SM) categories.
Summary
A search has been performed for a heavy scalar resonance $\mathrm{X}$ decaying to a Higgs boson (H) and Higgs-like scalar boson ($\mathrm{Y}$) in the two bottom quark ($ \mathrm{H}\to\mathrm{b}\overline{\mathrm{b}} $) and four quark ($ {{Y}} \to\mathrm{V}\mathrm{V}\to4\mathrm{q} $) final state, where V denotes a W or Z boson. The search targets Lorentz-boosted $ {{{X}} \to\mathrm{H}{{Y}} } $ production where the daughter quarks from the H and {Y} bosons are merged within large-area jets. A mass-decorrelated graph neural network tagger, known as PARTICLENET, was used to select $ \mathrm{H}\to\mathrm{b}\overline{\mathrm{b}} $ and $ \mathrm{V}\to\mathrm{q}\overline{\mathrm{q}} $ jets while the new ``particle transformer'' is used for fully merged $ \mathrm{Y} \to\mathrm{V}\mathrm{V} $ jets. The signal is extracted from a two-dimensional maximum likelihood fit to data binned in the reconstructed {Y} and {X} masses. The dominant quantum chromodynamics multijet background is estimated using control regions in data with the tagger score selections inverted. Other minor backgrounds including top quark and V plus jets are estimated using Monte Carlo simulations. No significant excess is observed in the data above the standard model background expectation. Upper limits at 95% confidence level as low as 0.2fb$^{-1}$ are set for the product of resonant $ \mathrm{X} \to\mathrm{H}{{Y}} \to\mathrm{b}\overline{\mathrm{b}}\mathrm{V}\mathrm{V}} $ production cross section and branching fraction for various {X} and {Y} masses. This is the first search for scalar resonances in the all-hadronic $ \mathrm{b}\overline{\mathrm{b}}\mathrm{V}\mathrm{V} $ final state at the LHC.
References
1 G. C. Branco et al. Theory and phenomenology of two-Higgs-doublet models Phys. Rept. 516 (2012) 1 1106.0034
2 N. Craig, J. Galloway, and S. Thomas Searching for signs of the second Higgs doublet 1305.2424
3 F. Domingo and S. Pa\ss ehr About the bosonic decays of heavy Higgs states in the (N)MSSM EPJC 82 (2022) 962 2207.05776
4 T. Robens, T. Stefaniak, and J. Wittbrodt Two-real-scalar-singlet extension of the SM: LHC phenomenology and benchmark scenarios EPJC 80 (2020) 151 1908.08554
5 CMS Collaboration Search for a massive scalar resonance decaying to a light scalar and a Higgs boson in the four b quarks final state with boosted topology PLB 842 (2023) 137392 2204.12413
6 CMS Collaboration Search for Higgs boson pair production in the $ \textrm{b}\overline{\textrm{b}}{\textrm{W}}^{+}{\textrm{W}}^{-} $ decay mode in proton-proton collisions at $ \sqrt{s} = $ 13 TeV JHEP 07 (2024) 293 CMS-HIG-21-005
2403.09430
7 CMS Collaboration Search for a new resonance decaying into two spin-0 bosons in a final state with two photons and two bottom quarks in proton-proton collisions at $ \sqrt{s} = $ 13 TeV JHEP 05 (2024) 316 CMS-HIG-21-011
2310.01643
8 CMS Collaboration Search for a heavy Higgs boson decaying into two lighter Higgs bosons in the $ \tau\tau\mathrm{b}\mathrm{b} $ final state at 13 TeV JHEP 11 (2021) 057 CMS-HIG-20-014
2106.10361
9 ATLAS Collaboration Search for a resonance decaying into a scalar particle and a Higgs boson in the final state with two bottom quarks and two photons with 199 fb$ ^{-1} $ of data collected at $ \sqrt{s}= $ 13 TeV and $ \sqrt{s}= $ 13.6 TeV with the ATLAS detector Submitted to \emphPhys. Lett. B, 2025 2510.02857
10 ATLAS Collaboration Search for a resonance decaying into a scalar particle and a Higgs boson in final states with leptons and two photons in proton-proton collisions at $ \sqrt{s} = $ 13 TeV with the ATLAS detector JHEP 10 (2024) 104 2405.20926
11 ATLAS Collaboration Search for a new heavy scalar particle decaying into a Higgs boson and a new scalar singlet in final states with one or two light leptons and a pair of $ \tau $-leptons with the ATLAS detector JHEP 10 (2023) 009 2307.11120
12 CMS Collaboration Precision luminosity measurement in proton-proton collisions at $ \sqrt{s} = $ 13 TeV in 2015 and 2016 at CMS EPJC 81 (2021) 800 CMS-LUM-17-003
2104.01927
13 CMS Collaboration Precision luminosity measurement in proton-proton collisions at $ \sqrt{s} = $ 13 TeV with the CMS detector CMS Physics Analysis Summary, 2025
CMS-PAS-LUM-20-001
CMS-PAS-LUM-20-001
14 CMS Collaboration HEPData record for this analysis link
15 CMS Collaboration The CMS experiment at the CERN LHC JINST 3 (2008) S08004
16 CMS Collaboration Development of the CMS detector for the CERN LHC Run 3 JINST 19 (2024) P05064 CMS-PRF-21-001
2309.05466
17 CMS Collaboration Performance of the CMS Level-1 trigger in proton-proton collisions at $ \sqrt{s} = $ 13 TeV JINST 15 (2020) P10017 CMS-TRG-17-001
2006.10165
18 CMS Collaboration The CMS trigger system JINST 12 (2017) P01020 CMS-TRG-12-001
1609.02366
19 CMS Collaboration Performance of the CMS high-level trigger during LHC Run 2 JINST 19 (2024) P11021 CMS-TRG-19-001
2410.17038
20 CMS Collaboration Particle-flow reconstruction and global event description with the CMS detector JINST 12 (2017) P10003 CMS-PRF-14-001
1706.04965
21 CMS Collaboration Technical proposal for the Phase-II upgrade of the Compact Muon Solenoid CMS Technical Proposal CERN-LHCC-2015-010, CMS-TDR-15-02, 2015
CDS
22 M. Cacciari, G. P. Salam, and G. Soyez The anti-$ k_{\mathrm{T}} $ jet clustering algorithm JHEP 04 (2008) 063 0802.1189
23 M. Cacciari, G. P. Salam, and G. Soyez FastJet user manual EPJC 72 (2012) 1896 1111.6097
24 D. Bertolini, P. Harris, M. Low, and N. Tran Pileup per particle identification JHEP 10 (2014) 059 1407.6013
25 CMS Collaboration Pileup mitigation at CMS in 13 TeV data JINST 15 (2020) P09018 CMS-JME-18-001
2003.00503
26 D. Krohn, J. Thaler, and L.-T. Wang Jet trimming JHEP 02 (2010) 084 0912.1342
27 H. Qu and L. Gouskos Jet tagging via particle clouds PRD 101 (2020) 056019 1902.08570
28 CMS Collaboration Mass regression of highly-boosted jets using graph neural networks CMS Detector Performance Summary CMS-DP-2021-017, 2021
CDS
29 CMS Collaboration Particle transformers for identifying Lorentz-boosted Higgs bosons decaying to a pair of W bosons CMS Physics Analysis Summary, 2025
CMS-PAS-JME-25-001
CMS-PAS-JME-25-001
30 J. Alwall et al. The automated computation of tree-level and next-to-leading order differential cross sections, and their matching to parton shower simulations JHEP 07 (2014) 079 1405.0301
31 T. Sjöstrand et al. An introduction to PYTHIA8.2 Comput. Phys. Commun. 191 (2015) 159 1410.3012
32 CMS Collaboration Extraction and validation of a new set of CMS PYTHIA8 tunes from underlying-event measurements EPJC 80 (2020) 4 CMS-GEN-17-001
1903.12179
33 P. Nason A new method for combining NLO QCD with shower Monte Carlo algorithms JHEP 11 (2004) 040 hep-ph/0409146
34 S. Frixione, P. Nason, and G. Ridolfi A positive-weight next-to-leading-order Monte Carlo for heavy flavour hadroproduction JHEP 09 (2007) 126 0707.3088
35 S. Frixione, P. Nason, and C. Oleari Matching NLO QCD computations with parton shower simulations: the POWHEG method JHEP 11 (2007) 070 0709.2092
36 S. Alioli, P. Nason, C. Oleari, and E. Re A general framework for implementing NLO calculations in shower Monte Carlo programs: the POWHEG box JHEP 06 (2010) 043 1002.2581
37 T. Je \v z o et al. An NLO+PS generator for $ {\mathrm{t}\overline{\mathrm{t}}} $ and $ \mathrm{W}\mathrm{t} $ production and decay including non-resonant and interference effects EPJC 76 (2016) 691 1607.04538
38 NNPDF Collaboration Parton distributions for the LHC Run II JHEP 04 (2015) 040 1410.8849
39 NNPDF Collaboration Parton distributions from high-precision collider data EPJC 77 (2017) 663 1706.00428
40 GEANT4 Collaboration GEANT 4 -- a simulation toolkit NIM A 506 (2003) 250
41 CMS Collaboration Identification of highly Lorentz-boosted heavy particles using graph neural networks and new mass decorrelation techniques CMS Detector Performance Summary CMS-DP-2020-002, 2020
CDS
42 CMS Collaboration Performance of heavy-flavour jet identification in Lorentz-boosted topologies in proton-proton collisions at $ \sqrt{s} = $ 13 TeV JINST 20 (2025) P11006 CMS-BTV-22-001
2510.10228
43 H. Qu, C. Li, and S. Qian Particle transformer for jet tagging in th Int. Conf. on Machine Learning, K. Chaudhuri et al., eds., volume 162, 2022
Proc. 3 (2022) 18281
2202.03772
44 CMS Collaboration Search for a massive resonance decaying to a pair of Higgs bosons in the four b quark final state in proton-proton collisions at $ \sqrt{s}= $ 13 TeV PLB 781 (2018) 244 1710.04960
45 CMS Collaboration Search for production of Higgs boson pairs in the four b quark final state using large-area jets in proton-proton collisions at $ \sqrt{s}= $ 13 TeV JHEP 01 (2019) 040 1808.01473
46 CMS Collaboration Inclusive search for highly boosted Higgs bosons decaying to bottom quark-antiquark pairs in proton-proton collisions at $ \sqrt{s} = $ 13 TeV JHEP 12 (2020) 85 CMS-HIG-19-003
2006.13251
47 S. N. Bernstein Démonstration du théor\`eme de Weierstrass fondée sur le calcul des probabilités Comm. Soc. Math. de Kharkov 13 (1912) 1
48 R. A. Fisher On the interpretation of $ \chi^2 $ from contingency tables, and the calculation of p J. Royal Stat. Soc. 85 (1922) 87
49 S. Baker and R. D. Cousins Clarification of the Use of Chi Square and Likelihood Functions in Fits to Histograms NIM A 221 (1984) 437
50 G. Cowan, K. Cranmer, E. Gross, and O. Vitells Asymptotic formulae for likelihood-based tests of new physics EPJC 71 (2011) 1554 1007.1727
51 F. A. Dreyer, G. P. Salam, and G. Soyez The Lund jet plane JHEP 12 (2018) 064 1807.04758
52 CMS Collaboration A method for correcting the substructure of multiprong jets using the Lund jet plane Submitted to JHEP, 2025 CMS-JME-23-001
2507.07775
53 CMS Collaboration Jet energy scale and resolution measurement with Run 2 legacy data collected by CMS at 13 TeV CMS Detector Performance Summary CMS-DP-2021-033, 2021
CDS
54 CMS Collaboration Search for nonresonant pair production of highly energetic Higgs bosons decaying to bottom quarks PRL 131 (2023) 041803 2205.06667
55 CMS Collaboration The CMS statistical analysis and combination tool: Combine Comput. Softw. Big Sci. 8 (2024) 19 CMS-CAT-23-001
2404.06614
56 W. Verkerke and D. P. Kirkby The RooFit toolkit for data modeling in th International Conference for Computing in High-Energy and Nuclear Physics (CHEP03)... [eConf C0303241, MOLT007], 2003
Proceedings of the 1 (2003) 3
physics/0306116
57 L. Moneta et al. The RooStats Project PoS ACAT 057, 2010
link
1009.1003
58 ATLAS and CMS Collaborations, and LHC Higgs Combination Group Procedure for the LHC Higgs boson search combination in Summer 2011 Technical Report CMS-NOTE-2011-005, ATL-PHYS-PUB-2011-11, 2011
59 T. Junk Confidence level computation for combining searches with small statistics NIM A 434 (1999) 435 hep-ex/9902006
60 A. L. Read Presentation of search results: The $ \text{CL}_\text{s} $ technique JPG 28 (2002) 2693
Compact Muon Solenoid
LHC, CERN