CMS logoCMS event Hgg
Compact Muon Solenoid
LHC, CERN

CMS-PAS-HIG-16-032
Search for $\mathrm{ H( {b\bar{b}} )H( {\gamma\gamma} )}$ decays at $ \sqrt{s} = $ 13 TeV
Abstract: The search for the production of a pair of Higgs bosons, in the final states containing two photons and two bottom quarks, is presented. Both resonant and nonresonant production mechanisms of the Higgs bosons are investigated. The analyzed data correspond to an integrated luminosity of 2.70 fb$^{-1}$ of pp collisions at $ \sqrt{s} = $ 13 TeV collected with the CMS detector in 2015. Upper limits on the production cross section of spin-0 and spin-2 particles are set at 95% confidence level. In addition, upper limits on nonresonant Higgs pair production are also set and compared to the SM Higgs pair production cross section via gluon fusion. The observed data is in agreement with standard model predictions.
Figures & Tables Summary References CMS Publications
Figures

png pdf
Figure 1:
Feynman diagrams that contribute to Higgs boson pair production by gluon fusion at leading order. On the left, via a heavy quark loop (with dominant top quark contribution). On the right, via the Higgs self coupling.

png pdf
Figure 2-a:
Signal $A\times \epsilon $ for (a) spin-0 radion and (b) spin-2 graviton signals, after each step of selection.

png pdf
Figure 2-b:
Signal $A\times \epsilon $ for (a) spin-0 radion and (b) spin-2 graviton signals, after each step of selection.

png pdf
Figure 3-a:
Distributions of (a) $ {M(\gamma \gamma )} $ and (b) $ {M(\mathrm{jj})} $ for 2.7 fb$^{-1}$ of data (points), expected SM single Higgs signal and background (shaded histograms), and resonant and nonresonant signals (lines). The 300 GeV Graviton, 600 GeV Radion and the SM nonresonant HH signal shapes are displayed, normalized assuming $\sigma ( \mathrm{ pp \to X \to HH \to \mathrm{ b \bar{b} }\gamma \gamma } ) = \sigma ( \mathrm{ pp \to HH \to \mathrm{ b \bar{b} }\gamma \gamma } ) = $ 100 fb. Events are required to pass the photon and jet selection described in Sec. 4 and to have at least one jet passing loose b-tagging requirements.

png pdf
Figure 3-b:
Distributions of (a) $ {M(\gamma \gamma )} $ and (b) $ {M(\mathrm{jj})} $ for 2.7 fb$^{-1}$ of data (points), expected SM single Higgs signal and background (shaded histograms), and resonant and nonresonant signals (lines). The 300 GeV Graviton, 600 GeV Radion and the SM nonresonant HH signal shapes are displayed, normalized assuming $\sigma ( \mathrm{ pp \to X \to HH \to \mathrm{ b \bar{b} }\gamma \gamma } ) = \sigma ( \mathrm{ pp \to HH \to \mathrm{ b \bar{b} }\gamma \gamma } ) = $ 100 fb. Events are required to pass the photon and jet selection described in Sec. 4 and to have at least one jet passing loose b-tagging requirements.

png pdf
Figure 4:
Pictorial description of the b-tag categorization scheme, in which the axes represent the b-tag score of the leading jet (y-axis) and of the subleading jet (x-axis). The low mass analysis categorization is shown on the left, while the nonresonant analysis categorization is shown on the right. HPC is the high-purity category, while MPC is the medium-purity category, and JCR is the jet control region.

png pdf
Figure 5:
Different 4-body mass reconstruction methods for different radion signal points: $ {\tilde{M}_{\mathrm {X}}} $ (dotted line) and $ {M(\mathrm{jj}\gamma \gamma )}$ (full line). Mass points are normalized independently, such that $ {M(\mathrm{jj}\gamma \gamma )}$ peak heights are set to unit. The different hypothesis mass samples are simulated assuming 1 GeV widths.

png pdf
Figure 6-a:
Distributions of (a) $ {M(\mathrm{jj}\gamma \gamma )}$ and (b) $ {\tilde{M}_{\mathrm {X}}} $ for 2.7 fb$^{-1}$ of data (points), expected SM single Higgs signal and background (shaded histograms), and resonant and nonresonant signals (lines). The 300 GeV Graviton, 600 GeV Radion and the SM nonresonant HH signal shapes are displayed, normalized assuming $\sigma (\mathrm{ pp\to X \to HH \to \mathrm{ b \bar{b} }\gamma \gamma } ) = \sigma ( \mathrm{ pp \to HH \to \mathrm{ b \bar{b} }\gamma \gamma } ) = $ 100 fb. Events are required to pass the photon and jet selection described in Sec. 4 and to have at least one jet passing loose b-tagging requirements.

png pdf
Figure 6-b:
Distributions of (a) $ {M(\mathrm{jj}\gamma \gamma )}$ and (b) $ {\tilde{M}_{\mathrm {X}}} $ for 2.7 fb$^{-1}$ of data (points), expected SM single Higgs signal and background (shaded histograms), and resonant and nonresonant signals (lines). The 300 GeV Graviton, 600 GeV Radion and the SM nonresonant HH signal shapes are displayed, normalized assuming $\sigma (\mathrm{ pp\to X \to HH \to \mathrm{ b \bar{b} }\gamma \gamma } ) = \sigma ( \mathrm{ pp \to HH \to \mathrm{ b \bar{b} }\gamma \gamma } ) = $ 100 fb. Events are required to pass the photon and jet selection described in Sec. 4 and to have at least one jet passing loose b-tagging requirements.

png pdf
Figure 7:
$ {M(\gamma \gamma )} $:$ {\tilde{M}_{\mathrm {X}}} $ distribution for the dominant background, modeled with diphoton+jets MC, after the full selection. The red and gray lines represent the mass window regions for resonances with masses of 250 GeV and 300 GeV, respectively.

png pdf
Figure 8-a:
Signal fits for the resonant 320 GeV Radion (spin-0) sample after selection for the high (a) and medium (b) purity b-tag categories. Upper plots show fits to $ {M(\gamma \gamma )} $, while lower plots show fits to $ {M(\mathrm{jj})} $.

png pdf
Figure 8-b:
Signal fits for the resonant 320 GeV Radion (spin-0) sample after selection for the high (a) and medium (b) purity b-tag categories. Upper plots show fits to $ {M(\gamma \gamma )} $, while lower plots show fits to $ {M(\mathrm{jj})} $.

png pdf
Figure 8-c:
Signal fits for the resonant 320 GeV Radion (spin-0) sample after selection for the high (a) and medium (b) purity b-tag categories. Upper plots show fits to $ {M(\gamma \gamma )} $, while lower plots show fits to $ {M(\mathrm{jj})} $.

png pdf
Figure 8-d:
Signal fits for the resonant 320 GeV Radion (spin-0) sample after selection for the high (a) and medium (b) purity b-tag categories. Upper plots show fits to $ {M(\gamma \gamma )} $, while lower plots show fits to $ {M(\mathrm{jj})} $.

png pdf
Figure 9-a:
Signal fits for the SM HH nonresonant sample after selection for the high (a) and medium (b) purity b-tag categories. Upper plots show fits to $ {M(\gamma \gamma )} $, while lower plots show fits to $ {M(\mathrm{jj})} $.

png pdf
Figure 9-b:
Signal fits for the SM HH nonresonant sample after selection for the high (a) and medium (b) purity b-tag categories. Upper plots show fits to $ {M(\gamma \gamma )} $, while lower plots show fits to $ {M(\mathrm{jj})} $.

png pdf
Figure 9-c:
Signal fits for the SM HH nonresonant sample after selection for the high (a) and medium (b) purity b-tag categories. Upper plots show fits to $ {M(\gamma \gamma )} $, while lower plots show fits to $ {M(\mathrm{jj})} $.

png pdf
Figure 9-d:
Signal fits for the SM HH nonresonant sample after selection for the high (a) and medium (b) purity b-tag categories. Upper plots show fits to $ {M(\gamma \gamma )} $, while lower plots show fits to $ {M(\mathrm{jj})} $.

png pdf
Figure 10-a:
Background-only fits for the resonant search with the 320 GeV mass window selection, for (a) $ {M(\gamma \gamma )} $ and (b) $ {M(\mathrm{jj})} $ , in the high purity b-tag category.

png pdf
Figure 10-b:
Background-only fits for the resonant search with the 320 GeV mass window selection, for (a) $ {M(\gamma \gamma )} $ and (b) $ {M(\mathrm{jj})} $ , in the high purity b-tag category.

png pdf
Figure 11-a:
Background-only fits for the resonant search with the 320 GeV mass window selection, for (a) $ {M(\gamma \gamma )} $ and (b) $ {M(\mathrm{jj})} $ , in the medium purity b-tag category.

png pdf
Figure 11-b:
Background-only fits for the resonant search with the 320 GeV mass window selection, for (a) $ {M(\gamma \gamma )} $ and (b) $ {M(\mathrm{jj})} $ , in the medium purity b-tag category.

png pdf
Figure 12-a:
Background fits for the nonresonant analysis, after selection, for (a) $ {M(\gamma \gamma )} $ and (b) $ {M(\mathrm{jj})} $ , high purity b-tag category.

png pdf
Figure 12-b:
Background fits for the nonresonant analysis, after selection, for (a) $ {M(\gamma \gamma )} $ and (b) $ {M(\mathrm{jj})} $ , high purity b-tag category.

png pdf
Figure 13-a:
Background fits for the nonresonant analysis, after selection, for (a) $ {M(\gamma \gamma )} $ and (b) $ {M(\mathrm{jj})} $ , medium purity b-tag category.

png pdf
Figure 13-b:
Background fits for the nonresonant analysis, after selection, for (a) $ {M(\gamma \gamma )} $ and (b) $ {M(\mathrm{jj})} $ , medium purity b-tag category.

png pdf
Figure 14-a:
Expected and observed upper limits on the production of $\mathrm{ pp \to X\to { {\mathrm {H}} {\mathrm {H}} } \to {\mathrm{ b \bar{b} }\gamma \gamma } } $ for the resonant analysis, assuming (a) spin-0 and (b) spin-2 for the resonance. Cross sections predicted by the RS WED model are shown (dashed red lines). The SM branching for $ { {\mathrm {H}} {\mathrm {H}} } \to {\mathrm{ b \bar{b} }\gamma \gamma } $ is assumed.

png pdf
Figure 14-b:
Expected and observed upper limits on the production of $\mathrm{ pp \to X\to { {\mathrm {H}} {\mathrm {H}} } \to {\mathrm{ b \bar{b} }\gamma \gamma } } $ for the resonant analysis, assuming (a) spin-0 and (b) spin-2 for the resonance. Cross sections predicted by the RS WED model are shown (dashed red lines). The SM branching for $ { {\mathrm {H}} {\mathrm {H}} } \to {\mathrm{ b \bar{b} }\gamma \gamma } $ is assumed.

png pdf
Figure 15:
Expected and observed upper limits on the nonresonant production of the SM process $\mathrm{ pp } \to { {\mathrm {H}} {\mathrm {H}} } \to {\mathrm{ b \bar{b} }\gamma \gamma } $, separated in analysis categories: high purity, medium purity and the final combined result.
Tables

png pdf
Table 1:
Summary of the search analysis methods used in this search.

png pdf
Table 2:
Summary of systematic uncertainties. The uncertainty in the b-tagging efficiency is anticorrelated between the b tag categories.
Summary
A search is performed by the CMS collaboration for resonant and nonresonant production of two Higgs bosons in the decay channel ${\mathrm{ H }\mathrm{ H }} \to \mathrm{ b \bar{b} } \gamma\gamma $, based on an integrated luminosity of 2.70 fb$^{-1}$ of $pp$ collisions collected at $ \sqrt{s} = $ 13 TeV in 2015. Resonances are sought in the mass range between 250 and 900 GeV. Expected and observed upper limits at a 95% CL are measured on the cross sections for the production of new particles decaying to ${\mathrm{ H }\mathrm{ H }} \to \mathrm{ b \bar{b} } \gamma \gamma $. The limits are compared to BSM predictions, based on the assumption of the existence of a warped extra dimension. No statistically significant deviations from the null hypothesis are found. The observed limits exclude radion (spin-0) mass points below 750 GeV, except in the vicinity of $M_{X} = $ 650 GeV, assuming $\Lambda_{R} = $ 1 TeV. For nonresonant production with SM-like kinematics, a 95% CL upper limit is set on $\sigma( \mathrm{ pp \to HH \to {b\bar{b} } \gamma \gamma } )$ at 7.90 fb.
References
1 S. Dawson, A. Ismail, and I. Low What is in the loop? The anatomy of double Higgs production PRD91 (2015), no. 11, 115008 1504.05596
2 Z. Heng, L. Shang, Y. Zhang, and J. Zhu Pair production of 125 GeV Higgs boson in the SM extension with color-octet scalars at the LHC JHEP 02 (2014) 083 1312.4260
3 R. Grober and M. Muhlleitner Composite Higgs Boson Pair Production at the LHC JHEP 06 (2011), no. KA-TP-37-2010, 020 1012.1562
4 M. Moretti et al. Higgs boson self-couplings at the LHC as a probe of extended Higgs sectors JHEP 02 (2005), no. FNT-T-2004-03, SHEP-03-12, CAFPE-26-03, UG-FT-156-03, LAPTH-1028-04, BA-TH-495-04, 024 hep-ph/0410334
5 W. D. Goldberger and M. B. Wise Modulus stabilization with bulk fields Phys.Rev.Lett. 83 (1999), no. CALT-68-2232, 4922--4925 hep-ph/9907447
6 O. DeWolfe, D. Freedman, S. Gubser, and A. Karch Modeling the fifth-dimension with scalars and gravity Phys.Rev. D62 (2000), no. HUTP-99-A048, MIT-CTP-2903, 046008 hep-th/9909134
7 C. Csaki, M. Graesser, L. Randall, and J. Terning Cosmology of brane models with radion stabilization Phys.Rev. D62 (2000), no. SCIPP-99-49, HUTP-A061, NSF-ITP-99-130, 045015 hep-ph/9911406
8 H. Davoudiasl, J. L. Hewett, and T. G. Rizzo Phenomenology of the Randall-Sundrum Gauge Hierarchy Model PRL 84 (2000), no. SLAC-PUB-8241 hep-ph/9909255
9 C. Csaki, M. L. Graesser, and G. D. Kribs Radion dynamics and electroweak physics Phys.Rev. D63 (2001), no. SCIPP-00-27, 065002 hep-th/0008151
10 K. Agashe, H. Davoudiasl, G. Perez, and A. Soni Warped Gravitons at the LHC and Beyond Phys.Rev. D76 (2007), no. SU-4252-843, BNL-HET-07-3, YITP-SB-07-02, 036006 hep-ph/0701186
11 A. Djouadi The Anatomy of electro-weak symmetry breaking. II. The Higgs bosons in the minimal supersymmetric model Phys.Rept. 459 (2008), no. LPT-ORSAY-05-18, 1--241 hep-ph/0503173
12 R. Barbieri et al. One or more Higgs bosons? 1307.4937
13 G. Branco et al. Theory and phenomenology of two-Higgs-doublet models Phys.Rept. 516 (2012) 1--102 1106.0034
14 M. Gouzevitch et al. Scale-invariant resonance tagging in multijet events and new physics in Higgs pair production JHEP 07 (2013) 148 1303.6636
15 B. Mellado Garcia, P. Musella, M. Grazzini, and R. Harlander CERN Report 4: Part I Standard Model Predictions Technical Report LHCHXSWG-DRAFT-INT-2016-008, May
16 D. de Florian and J. Mazzitelli Higgs pair production at next-to-next-to-leading logarithmic accuracy at the LHC 1505.07122
17 D. de Florian and J. Mazzitelli Higgs Boson Pair Production at Next-to-Next-to-Leading Order in QCD PRL 111 (2013) 201801 1309.6594
18 S. Borowka et al. Higgs boson pair production in gluon fusion at NLO with full top-quark mass dependence 1604.06447
19 S. Dawson, S. Dittmaier, and M. Spira Neutral Higgs boson pair production at hadron colliders: QCD corrections PRD58 (1998), no. DESY-98-028, BNL-HET-98-15, CERN-TH-98-135, 115012 hep-ph/9805244
20 J. Grigo, K. Melnikov, and M. Steinhauser Virtual corrections to Higgs boson pair production in the large top quark mass limit Nucl. Phys. B888 (2014), no. TTP14-024, LPN14-101, 17--29 1408.2422
21 CMS Collaboration Search for two Higgs bosons in final states containing two photons and two bottom quarks Submitted to: PRD (2016) CMS-HIG-13-032
1603.06896
22 ATLAS Collaboration Searches for Higgs boson pair production in the $ hh\to bb\tau\tau, \gamma\gamma WW^*, \gamma\gamma bb, bbbb $ channels with the ATLAS detector PRD92 (2015) 092004 1509.04670
23 ATLAS Collaboration Search for Higgs boson pair production in the $ b\bar{b}\gamma\gamma $ final state using pp collision data at $ \sqrt{s}=13 $ TeV with the ATLAS detector ATLAS-CONF-2016-004
24 CMS Collaboration The CMS experiment at the CERN LHC JINST 3 (2008) S08004 CMS-00-001
25 J. Alwall et al. MadGraph 5 : Going Beyond JHEP 06 (2011), no. FERMILAB-PUB-11-448-T, 128 1106.0522
26 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 1407 (2014), no. CERN-PH-TH-2014-064, CP3-14-18, LPN14-066, MCNET-14-09, ZU-TH-14-14, 079 1405.0301
27 A. Buckley et al. LHAPDF6: parton density access in the LHC precision era EPJC75 (2015), no. GLAS-PPE-2014-05, MCNET-14-29, IPPP-14-111, DCPT-14-222, 132 1412.7420
28 S. Carrazza, J. I. Latorre, J. Rojo, and G. Watt A compression algorithm for the combination of PDF sets EPJC 75 (2015), no. TIF-UNIMI-2015-2, OUTP-15-01P, IPPP-15-22, DCPT-15-44, 474 1504.06469
29 J. Butterworth et al. PDF4LHC recommendations for LHC Run II 1510.03865
30 S. Dulat et al. The CT14 Global Analysis of Quantum Chromodynamics 1506.07443
31 L. A. Harland-Lang, A. D. Martin, P. Motylinski, and R. S. Thorne Parton distributions in the LHC era: MMHT 2014 PDFs EPJC 75 (2015), no. LCTS-2014-47, IPPP-14-97, DCPT-14-194, 204 1412.3989
32 NNPDF Collaboration Parton distributions for the LHC Run II JHEP 04 (2015), no. EDINBURGH-2014-15, IFUM-1034-FT, CERN-PH-TH-2013-253, OUTP-14-11P, CAVENDISH-HEP-14-11, 040 1410.8849
33 T. Sjostrand, S. Mrenna, and P. Z. Skands A Brief Introduction to PYTHIA 8.1 CPC 178 (2008), no. CERN-LCGAPP-2007-04, LU-TP-07-28, FERMILAB-PUB-07-512-CD-T, 852--867 0710.3820
34 J. Allison et al. Geant4 developments and applications IEEE Trans. Nucl. Sci. 53 (2006) 270
35 A. Oliveira and R. Rosenfeld Hidden sector effects on double higgs production near threshold at the LHC Phys.Lett. B702 (2011) 201--204 1009.4497
36 P. de Aquino, K. Hagiwara, Q. Li, and F. Maltoni Simulating graviton production at hadron colliders JHEP 06 (2011) 132 1101.5499
37 Alexandra Oliveira Models for collider study di-Higgs channel twiki
38 Claude Duhr The Higgs effective theory web page
39 A. Alloul et al. FeynRules 2.0 - A complete toolbox for tree-level phenomenology CPC 185 (2014), no. CERN-PH-TH-2013-239, MCNET-13-14, IPPP-13-71, DCPT-13-142, --PITT-PACC-1308, 2250--2300 1310.1921
40 J. Baglio et al. The measurement of the Higgs self-coupling at the LHC: theoretical status JHEP 1304 (2013), no. KA-TP-44-2012, SFB-CPP-12-102, LPT-ORSAY-12-124, PSI-PR-12-10, 151 1212.5581
41 T. Gleisberg et al. Event generation with SHERPA 1.1 JHEP 02 (2009), no. FERMILAB-PUB-08-477-T, SLAC-PUB-13420, ZU-TH-17-08, DCPT-08-138, IPPP-08-69, EDINBURGH-2008-30, MCNET-08-14, 007 0811.4622
42 CMS Collaboration First results on Higgs to gammagamma at 13 TeV Technical Report CMS-PAS-HIG-15-005
43 CMS Collaboration Electron and Photon performance using data collected by CMS at $ \sqrt{s} $ = 13 TeV and 25ns CDS
44 CMS Collaboration Search for the standard model Higgs boson produced through vector boson fusion and decaying to bb with proton-proton collisions at sqrt(s) = 13 TeV CMS-PAS-HIG-16-003 CMS-PAS-HIG-16-003
45 CMS Collaboration Identification of b quark jets at the CMS Experiment in the LHC Run 2 Technical Report CMS-PAS-BTV-15-001
Compact Muon Solenoid
LHC, CERN