CMS logoCMS event Hgg
Compact Muon Solenoid
LHC, CERN

CMS-PAS-B2G-24-001
Search for a new scalar resonance decaying to a Higgs boson and a new scalar with two bottom quarks and two photons in the final state in proton-proton collisions at $ \sqrt{s}= $ 13 TeV
Abstract: A search for a new scalar resonance, X, decaying to a standard model Higgs boson and a new scalar particle, Y, is presented. The Higgs boson further decays to a pair of b quarks, while the Y particle decays to a pair of photons. The search is performed in the mass range 240-1000 GeV for the resonance X, and in the mass range 70-800 GeV for the particle Y, using proton-proton collision data collected by the CMS experiment at $ \sqrt{s}= $ 13 TeV, corresponding to an integrated luminosity of 132 fb$ ^{-1} $. Observed (expected) upper limits on the product of the production cross section and the relevant branching fraction at the 95% confidence level are extracted for the $ \text{X}\to\text{Y}\text{H} $ process, and are found to be within the range of 0.05 $ -2.69\,(0.08-1.94) $ fb depending on $ m_\text{X} $. The most significant deviation from the background-only hypothesis is observed for X and Y masses of 300 and 77 GeV, respectively, with a local (global) significance of 3.33 $ \,(0.65) $ standard deviations.
Figures & Tables Summary References CMS Publications
Figures

png pdf
Figure 1:
Feynman diagram for the production of the BSM resonance $ \mathrm{X} $ and its subsequent decay to two scalars, one SM Higgs boson and one BSM scalar $ \text{Y} $, with $ \mathrm{H}\to\mathrm{b}\overline{\mathrm{b}} $ and $ \text{Y}\to\gamma\gamma $

png pdf
Figure 2:
Distributions of the transformed pNN score for the signal hypotheses of $ m_\text{X}= $ 280 GeV, $ m_\text{Y}= $ 125 GeV (left) and $ m_\text{X}= $ 600 GeV, $ m_\text{Y}= $ 70 GeV (right) in their corresponding SRs. The bin boundaries correspond to the SR boundaries of each mass point. The distributions are inclusive in the $ m_{\gamma\gamma} $ distribution. The grey bands in the lower panels show the statistical uncertainty on the background estimation.

png pdf
Figure 2-a:
Distributions of the transformed pNN score for the signal hypotheses of $ m_\text{X}= $ 280 GeV, $ m_\text{Y}= $ 125 GeV (left) and $ m_\text{X}= $ 600 GeV, $ m_\text{Y}= $ 70 GeV (right) in their corresponding SRs. The bin boundaries correspond to the SR boundaries of each mass point. The distributions are inclusive in the $ m_{\gamma\gamma} $ distribution. The grey bands in the lower panels show the statistical uncertainty on the background estimation.

png pdf
Figure 2-b:
Distributions of the transformed pNN score for the signal hypotheses of $ m_\text{X}= $ 280 GeV, $ m_\text{Y}= $ 125 GeV (left) and $ m_\text{X}= $ 600 GeV, $ m_\text{Y}= $ 70 GeV (right) in their corresponding SRs. The bin boundaries correspond to the SR boundaries of each mass point. The distributions are inclusive in the $ m_{\gamma\gamma} $ distribution. The grey bands in the lower panels show the statistical uncertainty on the background estimation.

png pdf
Figure 3:
Parametric models of the signal process for $ m_\text{X}= $ 600 GeV, $ m_\text{Y}= $ 70 GeV (left), and for $ m_\text{X}= $ 1000 GeV, $ m_\text{Y}= $ 800 GeV (right) in their most sensitive signal region. The histograms are normalized to unity. The acronym ``dof" stands for the numbers of degrees of freedom of the parametric model.

png
Figure 3-a:
Parametric models of the signal process for $ m_\text{X}= $ 600 GeV, $ m_\text{Y}= $ 70 GeV (left), and for $ m_\text{X}= $ 1000 GeV, $ m_\text{Y}= $ 800 GeV (right) in their most sensitive signal region. The histograms are normalized to unity. The acronym ``dof" stands for the numbers of degrees of freedom of the parametric model.

png
Figure 3-b:
Parametric models of the signal process for $ m_\text{X}= $ 600 GeV, $ m_\text{Y}= $ 70 GeV (left), and for $ m_\text{X}= $ 1000 GeV, $ m_\text{Y}= $ 800 GeV (right) in their most sensitive signal region. The histograms are normalized to unity. The acronym ``dof" stands for the numbers of degrees of freedom of the parametric model.

png pdf
Figure 4:
Background-only fit (red line) and signal+background fit (blue line) for the mass point hypothesis of $ m_\text{X}= $ 280 GeV, $ m_\text{Y}= $ 90 GeV. The red line in the lower panel shows the background-only fit, which is by definition zero, and it is added as a visual aid. From left to right, the first and second most sensitive signal regions are shown. The choice of the background functional form is determined by the maximum likelihood fit.

png
Figure 4-a:
Background-only fit (red line) and signal+background fit (blue line) for the mass point hypothesis of $ m_\text{X}= $ 280 GeV, $ m_\text{Y}= $ 90 GeV. The red line in the lower panel shows the background-only fit, which is by definition zero, and it is added as a visual aid. From left to right, the first and second most sensitive signal regions are shown. The choice of the background functional form is determined by the maximum likelihood fit.

png
Figure 4-b:
Background-only fit (red line) and signal+background fit (blue line) for the mass point hypothesis of $ m_\text{X}= $ 280 GeV, $ m_\text{Y}= $ 90 GeV. The red line in the lower panel shows the background-only fit, which is by definition zero, and it is added as a visual aid. From left to right, the first and second most sensitive signal regions are shown. The choice of the background functional form is determined by the maximum likelihood fit.

png pdf
Figure 5:
Observed\,(expected) upper limits on the $ \sigma \times BR $ for the $ \text{X}\to\text{Y}(\gamma\gamma)\mathrm{H}(\mathrm{b}\overline{\mathrm{b}}) $ signal with the different mass hypotheses are shown with solid\,(dashed) lines. The green and the yellow bands define the $ \pm 68% $ and $ \pm 95% $ uncertainty bands, respectively. For visualization purposes, the upper limit for mass points with different $ m_\text{X} $ has been multiplied with a constant factor quoted on the right of each band. Mass points with $ m_\text{X} $ ranging from 240 to 550 GeV are shown.

png pdf
Figure 6:
Observed\,(expected) upper limits on the $ \sigma \times BR $ for the $ \text{X}\to\text{Y}(\gamma\gamma)\mathrm{H}(\mathrm{b}\overline{\mathrm{b}}) $ signal with the different mass hypotheses are shown with solid\,(dashed) lines. The green and the yellow bands define the $ \pm 68% $ and $ \pm 95% $ uncertainty bands, respectively. For visualization purposes, the upper limit for mass points with different $ m_\text{X} $ has been multiplied with a constant factor quoted on the right of each band. Mass points with $ m_\text{X} $ ranging from 600 to 1000 GeV are shown.

png pdf
Figure 7:
Observed\,(expected) upper limits on the $ \sigma \times BR $ for the $ \text{X}\to\text{Y}(\gamma\gamma)\mathrm{H}(\mathrm{b}\overline{\mathrm{b}}) $ signal with the different $ m_\text{Y} $ hypotheses are shown with solid\,(dashed) lines. The inner (green) band and the outer (yellow) band indicate the regions containing 68% and 95%, respectively, of the distribution of limits expected under the background-only hypothesis. The left plot shows the set of mass points with the lowest $ m_\text{X} = $ 240 GeV, and the right plot shows the set of mass points with the highest $ m_\text{X} = $ 1000 GeV.

png
Figure 7-a:
Observed\,(expected) upper limits on the $ \sigma \times BR $ for the $ \text{X}\to\text{Y}(\gamma\gamma)\mathrm{H}(\mathrm{b}\overline{\mathrm{b}}) $ signal with the different $ m_\text{Y} $ hypotheses are shown with solid\,(dashed) lines. The inner (green) band and the outer (yellow) band indicate the regions containing 68% and 95%, respectively, of the distribution of limits expected under the background-only hypothesis. The left plot shows the set of mass points with the lowest $ m_\text{X} = $ 240 GeV, and the right plot shows the set of mass points with the highest $ m_\text{X} = $ 1000 GeV.

png
Figure 7-b:
Observed\,(expected) upper limits on the $ \sigma \times BR $ for the $ \text{X}\to\text{Y}(\gamma\gamma)\mathrm{H}(\mathrm{b}\overline{\mathrm{b}}) $ signal with the different $ m_\text{Y} $ hypotheses are shown with solid\,(dashed) lines. The inner (green) band and the outer (yellow) band indicate the regions containing 68% and 95%, respectively, of the distribution of limits expected under the background-only hypothesis. The left plot shows the set of mass points with the lowest $ m_\text{X} = $ 240 GeV, and the right plot shows the set of mass points with the highest $ m_\text{X} = $ 1000 GeV.
Tables

png pdf
Table 1:
Additional photon requirements, as a function of $ |\eta| $ and $ R_{9} $. The variable $ \rho $ is the median of the transverse energy density per unit area in the event.

png pdf
Table 2:
The training variables included as input to the pNN used for the final selection of this search.
Summary
A search for the production of a new scalar resonance, $ \mathrm{X} $, which decays to a standard model Higgs boson and a new scalar resonance, $ \text{Y} $, has been presented. The final state involves a pair of b quarks from the Higgs boson decay, and a pair of photons, from the $ \text{Y} $ particle decay. This search is the first one targeting this final state combination. The analysis probes the mass range of 240--1000 GeV for the resonance X and of 70--800 GeV for the particle $ \text{Y} $, and uses proton-proton collision data collected by the CMS experiment at $ \sqrt{s} = $ 13 TeV, corresponding to an integrated luminosity of 132 fb$^{-1}$. Upper limits on the $ \mathrm{X}\to\text{Y}\mathrm{H} $ cross section have been derived at the 95% confidence level as functions of the masses of the $ \mathrm{X} $ and the $ \text{Y} $ particles. The observed\,(expected) upper limits on the product of the production cross section of $ \mathrm{X} $ and the branching fraction to the $ \mathrm{b}\overline{\mathrm{b}}\gamma\gamma $ final state are found to be between 0.05--2.69\,(0.08--1.94)\,fb depending on the specific signal masses hypothesis, and are compatible with the background-only expectation. A local\,(global) significance of 3.33\,(0.65) standard deviations has been observed for the most significant deviation from the background-only hypothesis, corresponding to $ m_\text{X}= $ 300 GeV and $ m_\text{Y}= $ 77 GeV.
References
1 ATLAS Collaboration Observation of a new particle in the search for the standard model Higgs boson with the detector at the LHC PLB 716 (2012) 1 1207.7214
2 CMS Collaboration Observation of a new boson at a mass of 125 GeV with the CMS experiment at the LHC PLB 716 (2012) 30 CMS-HIG-12-028
1207.7235
3 ATLAS Collaboration A detailed map of higgs boson interactions by the ATLAS experiment ten years after the discovery Nature 607 (2022) 52 2207.00092
4 CMS Collaboration A portrait of the Higgs boson by the cms experiment ten years after the discovery Nature 607 (2022) 60 CMS-HIG-22-001
2207.00043
5 U. Ellwanger, C. Hugonie, and A. M. Teixeira The Next-to-Minimal Supersymmetric Standard Model Physics Reports 496 (2010) 1 0910.1785
6 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
7 H. Abouabid et al. Benchmarking di-Higgs production in various extended Higgs sector models JHEP 2022 (2022) 011 2112.12515
8 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 Physics Letters B 842 (2023) 137392 2204.12413
9 CMS Collaboration Search for a heavy Higgs boson decaying into two lighter Higgs bosons in the $ \tau\tau $bb final state at 13 TeV JHEP 2021 (2021) 057 CMS-HIG-20-014
2106.10361
10 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 2024 (2024) 316 CMS-HIG-21-011
2310.01643
11 CMS Collaboration Search for a standard model-like Higgs boson in the mass range between 70 and 110 GeV in the diphoton final state in proton-proton collisions at $ \sqrt{s} $ = 13 TeV \hrefhttp://www.arXiv.org/abs/arXiv:2405.18149\textttarXiv:arXiv:2405.18149, 2025
Physics Letters B 860 (2025) 139067
12 CMS Collaboration The CMS experiment at the CERN LHC JINST 3 (2008) S08004
13 CMS Collaboration Development of the CMS detector for the CERN LHC Run 3 JINST 19 (2024) P05064 CMS-PRF-21-001
2309.05466
14 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
15 CMS Collaboration The CMS trigger system JINST 12 (2017) P01020 CMS-TRG-12-001
1609.02366
16 CMS Collaboration Performance of the CMS high-level trigger during LHC Run 2 JINST 19 (2024) P11021 CMS-TRG-19-001
2410.17038
17 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
18 CMS Collaboration Particle-flow reconstruction and global event description with the CMS detector JINST 12 (2017) P10003 CMS-PRF-14-001
1706.04965
19 CMS Collaboration Performance of photon reconstruction and identification with the CMS detector in proton-proton collisions at sqrt(s) = 8 TeV JINST 10 (2015) P08010 CMS-EGM-14-001
1502.02702
20 CMS Collaboration A measurement of the Higgs boson mass in the diphoton decay channel PLB 805 (2020) 135425 CMS-HIG-19-004
2002.06398
21 CMS Collaboration Electron and photon reconstruction and identification with the CMS experiment at the CERN LHC JINST 16 (2021) P05014 CMS-EGM-17-001
2012.06888
22 CMS Collaboration ECAL 2016 refined calibration and Run2 summary plots CMS Detector Performance Summary CMS-DP-2020-021, 2020
CDS
23 CMS Collaboration Performance of the CMS muon detector and muon reconstruction with proton-proton collisions at $ \sqrt{s}= $ 13 TeV JINST 13 (2018) P06015 CMS-MUO-16-001
1804.04528
24 M. Cacciari, G. P. Salam, and G. Soyez The anti-$ k_{\mathrm{T}} $ jet clustering algorithm JHEP 04 (2008) 063 0802.1189
25 M. Cacciari, G. P. Salam, and G. Soyez Fastjet user manual EPJC 72 (2012) 1896 1111.6097
26 CMS Collaboration Pileup mitigation at CMS in 13 TeV data JINST 15 (2020) P09018 CMS-JME-18-001
2003.00503
27 CMS Collaboration Jet energy scale and resolution in the CMS experiment in pp collisions at 8 TeV JINST 12 (2017) P02014 CMS-JME-13-004
1607.03663
28 E. Bols et al. Jet flavour classification using DeepJet JINST 15 (2020) P12012 2008.10519
29 CMS Collaboration Performance of the DeepJet b tagging algorithm using 41.9 fb of data from proton-proton collisions at 13 TeV with Phase 1 CMS detector CMS Detector Performance Note CMS-DP-2018-058, 2018
CDS
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 2014 (2014) 079 1405.0301
31 R. Frederix et al. The automation of next-to-leading order electroweak calculations JHEP 2018 (2018) 185 1804.10017
32 D. de Florian et al. Handbook of LHC Higgs cross sections: 4. Deciphering the nature of the Higgs sector CERN Report CERN-2017-002-M, 2016
link
1610.07922
33 S. Frixione, P. Nason, and C. Oleari Matching NLO QCD computations with parton shower simulations: the POWHEG method JHEP 11 (2007) 070 0709.2092
34 S. Alioli, P. Nason, C. Oleari, and E. Re NLO single-top production matched with shower in POWHEG: s- and t-channel contributions [Erratum: JHEP 02, 011 ()], 2009
JHEP 09 (2009) 111
0907.4076
35 T. Melia, P. Nason, R. Rontsch, and G. Zanderighi W+W-, WZ and ZZ production in the POWHEG BOX JHEP 11 (2011) 078 1107.5051
36 P. Nason and G. Zanderighi $ W^+ W^- $, $ W Z $ and $ Z Z $ production in the POWHEG-BOX-V2 EPJC 74 (2014) 2702 1311.1365
37 E. Bagnaschi, G. Degrassi, P. Slavich, and A. Vicini Higgs production via gluon fusion in the POWHEG approach in the SM and in the MSSM JHEP 02 (2012) 088 1111.2854
38 G. Heinrich, S. P. Jones, M. Kerner, and L. Scyboz A non-linear EFT description of $ gg\to hh $ at NLO interfaced to POWHEG JHEP 10 (2020) 021 2006.16877
39 G. Heinrich et al. Probing the trilinear Higgs boson coupling in di-Higgs production at NLO QCD including parton shower effects JHEP 06 (2019) 066 1903.08137
40 S. Jones and S. Kuttimalai Parton shower and NLO-matching uncertainties in Higgs boson pair production JHEP 02 (2018) 176 1711.03319
41 G. Heinrich et al. NLO predictions for Higgs boson pair production with full top quark mass dependence matched to parton showers JHEP 08 (2017) 088 1703.09252
42 G. Buchalla et al. Higgs boson pair production in non-linear effective field theory with full $ m_t $-dependence at NLO QCD JHEP 09 (2018) 057 1806.05162
43 T. Gleisberg et al. Event generation with SHERPA 1.1 JHEP 02 (2009) 007 0811.4622
44 T. Sjöstrand et al. An introduction to PYTHIA 8.2 Comput. Phys. Commun. 191 (2015) 159 1410.3012
45 CMS Collaboration Event generator tunes obtained from underlying event and multiparton scattering measurements EPJC 76 (2016) 155 CMS-GEN-14-001
1512.00815
46 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
47 NNPDF Collaboration Parton distributions from high-precision collider data EPJC 77 (2017) 663 1706.00428
48 \GEANTfour Collaboration GEANT 4 --- a simulation toolkit NIM A 506 (2003) 250
49 CMS Collaboration Performance of missing transverse momentum reconstruction in proton-proton collisions at $ \sqrt{s} = $ 13\,TeV using the CMS detector JINST 14 (2019) P07004 CMS-JME-17-001
1903.06078
50 CMS Collaboration Measurements of $ t\overline{t}H $ production and the $ CP $ structure of the Yukawa interaction between the Higgs boson and top quark in the diphoton decay channel PRL 125 (2020) 061801 CMS-HIG-19-013
2003.10866
51 P. Baldi et al. Parameterized neural networks for high-energy physics EPJC 76 (2016) 1601.07913
52 M. J. Oreglia A study of the reactions $ \psi^\prime \to \gamma \gamma \psi $ PhD thesis, Stanford University, . SLAC Report SLAC-R-236, see Appendix D, 1980
link
53 J. E. Gaiser Charmonium spectroscopy from radiative decays of the $ J/\psi $ and $ \psi^\prime $ PhD thesis, Stanford University, . SLAC Report SLAC-R-255, 1982
link
54 S. Choi and H. Oh Improved extrapolation methods of data-driven background estimations in high energy physics EPJC 81 (2021)
55 P. D. Dauncey, M. Kenzie, N. Wardle, and G. J. Davies Handling uncertainties in background shapes JINST 10 (2015) P04015 1408.6865
56 J. Baglio et al. $ gg \rightarrow HH $: Combined uncertainties PRD 103 (2021)
57 J. Butterworth et al. PDF4LHC recommendations for LHC Run II JPG 43 (2016) 023001 1510.03865
58 S. Heinemeyer et al. Handbook of LHC Higgs cross sections: 3. Higgs properties \hrefhttp://www.arXiv.org/abs/arXiv:1307.1347\textttarXiv:arXiv:1307.1347
link
59 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
60 CMS Collaboration CMS luminosity measurement for the 2017 data-taking period at $ \sqrt{s} $ = 13 TeV CMS Physics Analysis Summary, 2018
link
CMS-PAS-LUM-17-004
61 CMS Collaboration CMS luminosity measurement for the 2018 data-taking period at $ \sqrt{s} $ = 13 TeV CMS Physics Analysis Summary, 2019
link
CMS-PAS-LUM-18-002
62 CMS Collaboration Measurement of the inclusive $ W $ and $ Z $ production cross sections in $ pp $ collisions at $ \sqrt{s}= $ 7 TeV JHEP 10 (2011) 132 CMS-EWK-10-005
1107.4789
63 CMS Collaboration The CMS statistical analysis and combination tool: Combine Comput. Softw. Big Sci. 8 (2024) 19 CMS-CAT-23-001
2404.06614
64 W. Verkerke and D. Kirkby The \textscRooFit toolkit for data modeling in th International Conference on Computing in High Energy and Nuclear Physics (CHEP ): La Jolla CA, United States, March 24--28, . . . [eConf C0303241 MOLT007], 2003
Proc. 1 (2003) 3
physics/0306116
65 L. Moneta et al. The \textscRooStats project in th International Workshop on Advanced Computing and Analysis Techniques in Physics Research (ACAT ): Jaipur, India, February 22--27, . . . [PoS (ACAT) 057], 2010
Proc. 1 (2010) 3
1009.1003
66 E. Gross and O. Vitells Trial factors for the look elsewhere effect in high energy physics EPJC 70 (2010) 1005.1891
Compact Muon Solenoid
LHC, CERN