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 | ||
CMS Collaboration | ||
6 April 2025 | ||
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. | ||
Links: CDS record (PDF) ; CADI line (restricted) ; |
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 |
![]() |
![]() |
![]() |
![]() |
![]() |
![]() |