CMS-PAS-HIG-17-017

CMS-PAS-HIG-17-017
Search for non-resonant Higgs pair-production in the $\mathrm{b}\overline{\mathrm{b}}\mathrm{b}\overline{\mathrm{b}}$ final state with the CMS detector
CMS Collaboration
July 2018

Abstract: Results of a search for non-resonant production of Higgs boson pairs decaying to two $\rm{b \bar{b}}$ pairs are presented. This search uses data from proton-proton collisions at a centre-of-mass energy of 13 TeV, collected by the CMS experiment at the LHC in 2016 and corresponding to an integrated luminosity of 35.9 fb$^{-1}$. No signal is observed, and an upper limit of 847fb is set on the standard model cross section for non-resonant Higgs boson pair-production, where the Higgs bosons decay to $\rm{b \bar{b}}$ pairs. The same signature is studied in the context of theories of physics beyond the standard model that predict modified couplings of the Higgs boson, and used to extract limits on Higgs boson anomalous couplings.
Links: CDS record (PDF) ; inSPIRE record ; CADI line (restricted) ; These preliminary results are superseded in this paper, JHEP 04 (2019) 112. The superseded preliminary plots can be found here.

Figures & Tables	Summary	References	CMS Publications

Figures
png pdf	Figure 1: Feynman diagrams that contribute to HH production by gluon-gluon fusion at leading order. Diagrams (a) and (b) correspond to SM-like processes, while diagrams (c), (d), and (e) correspond to pure BSM effects: (c) and (d) describe contact interactions between the H and gluons, and (e) exploits the contact interaction of two Higgs bosons with top quarks.
png pdf	Figure 2: Left: comparison of the distribution of BDT output for real data (histogram) selected in a region of the leading versus trailing Higgs mass plane that excludes a 60-GeV-wide box around the most probable values of the di-jet masses of signal events, with the corresponding output on an artificial sample obtained from the same data set by hemisphere mixing. Right: bin-by-bin differences between data and model, in standard deviation units before (top) and after (bottom) bias correction; pull distribution for the differences, fit to a Gaussian distribution. Bias correction uncertainty is increased to take the sigma of the residuals to one.
png pdf	Figure 3: Results of the fit to the BDT distribution for the SM HH production signal. In the bottom panel a comparison is shown between the best-fit signal and best-fit background subtracted from measured data. The band, centred at zero, shows the total uncertainty.
png pdf	Figure 4: Cross section limits on the 13 BSM models investigated. See Table xxxxx for their respective parameter values.
png pdf	Figure 5: Cross section limits on $\rm {H({{\mathrm {b}} {\overline {\mathrm {b}}}})H({{\mathrm {b}} {\overline {\mathrm {b}}}})}$ for values of $k_\lambda $ in the [-20;20] range, assuming $k_t = $ 1; the theoretical prediction with $k_t = $ 1 is also shown.

Tables
png pdf	Table 1: The values of the anomalous coupling parameter for the the 13 benchmark models studied\nobreakspace {} [19]. For reference the values of the parameters in the SM are also included.
png pdf	Table 2: Cut-flow efficiency for the SM like signal $\rm {pp} \rightarrow \rm {HH} \rightarrow \rm {b \bar{b} }$; the efficiency and the relative reduction of each successive selection is shown. The number of expected events for an integrated luminosity of 1fb$^{-1}$ is also reported.
png pdf	Table 3: List of BDT input variables.
png pdf	Table 4: Systematic uncertainties considered in the analysis and relative impact on the expected limit for SM HH production. The relative impact is obtained by fixing the nuisance parameters corresponding to each source and recalculating the expected limit.
png pdf	Table 5: The expected upper limit on $\sigma (\rm {pp \rightarrow H({{\mathrm {b}} {\overline {\mathrm {b}}}})H({{\mathrm {b}} {\overline {\mathrm {b}}}})})$ in the SM at 95% confidence level using 35.9 fb$^{-1}$ of data considering the whole set of nuisance parameters.
png pdf	Table 6: Expected and observed limits on the HH production cross section for the 13 BSM benchmark models.

Summary

This document presents a search for non-resonant HH production with both Higgs bosons decaying into $\rm{b \bar{b}}$ pairs. The SM production has been studied along with 13 BSM benchmark models, using a data set of $\sqrt{s}= $ 13 TeV proton-proton collision events corresponding to an integrated luminosity of 35.9 fb$^{-1}$ collected by the CMS detector during the 2016 LHC run. The analysis of events acquired by a hadronic multijet trigger includes the selection of events with 4 b-tagged jets and a classification using boosted decision trees targeting each of the signals considered. Limits on the HH production cross section times the square of the branching fraction for the H decay to b-quark pairs are extracted for the SM and each BSM model considered using binned likelihood fits of the shape of the BDT classifier output. The background model is derived from a novel data-driven technique that provides a multi-dimensional representation of the QCD multijet background, which makes up the major part of the data, and also models well the overall background distribution. The expected upper limit on $\sigma(\rm{pp \rightarrow H(\mathrm{b\bar{b}})H(\mathrm{b\bar{b}})})_{SM}$ is 419 fb. The observed upper limit is 847 fb, which corresponds to approximately 74 times the expected value for the SM process. Anomalous couplings of the H are also investigated. The upper limits extracted for the HH production cross section in the 13 BSM benchmark models range from 508 to 3513 fb.

References
1	CMS Collaboration	Observation of a new boson at a mass of 125 GeV with the CMS experiment at the LHC	Physics Letters B 716 (2012), no. 1
2	ATLAS Collaboration	Observation of a new particle in the search for the standard model higgs boson with the ATLAS detector at the LHC	Physics Letters B 716 (2012), no. 1
3	ATLAS and C. Collaborations	Measurements of the Higgs boson production and decay rates and constraints on its couplings from a combined ATLAS and CMS analysis of the lhc pp collision data at $ \sqrt{s} = $ 7 and 8 TeV	JHEP 2016 (Aug, 2016) 45
4	LHC Higgs Cross Section Working Group Collaboration	Handbook of LHC Higgs Cross Sections: 4. Deciphering the Nature of the Higgs Sector		1610.07922
5	B. Mellado Garcia, P. Musella, M. Grazzini, and R. Harlander	CERN Report 4: Part I Standard Model Predictions	(May, 2016)
6	D. de Florian and J. Mazzitelli	Higgs Boson Pair Production at Next-to-Next-to-Leading Order in QCD	PRL 111 (2013)	1309.6594
7	S. Dawson, S. Dittmaier, and M. Spira	Neutral Higgs boson pair production at hadron colliders: QCD corrections	PRD58 (1998)	hep-ph/9805244
8	S. Borowka et al.	Higgs boson pair production in gluon fusion at next-to-leading order with full top-quark mass dependence	PRL 117 (Jun, 2016) 012001
9	D. de Florian and J. Mazzitelli	Higgs pair production at next-to-next-to-leading logarithmic accuracy at the LHC	JHEP 09 (2015)	1505.07122
10	CMS Collaboration	Search for Higgs boson pair production in events with two bottom quarks and two tau leptons in proton-proton collisions at $ \sqrt{s} = $ 13 TeV	Physics Letters B 778 (2018) 101–127
11	ATLAS Collaboration	Search for pair production of Higgs bosons in the $ b\bar{b}b\bar{b} $ final state using proton-proton collisions at $ \sqrt{s} = $ 13 TeV with the ATLAS detector	Submitted to JHEP	1804.06174
12	A. Carvalho et al.	Analytical parametrization and shape classification of anomalous HH production in the EFT approach		1608.06578
13	G. Aad et al.	Search for higgs boson pair production in the $ b\bar{b}b\bar{b} $ final state from pp collisions at $ \sqrt{s}= $ 8 TeV with the ATLAS detector	The European Physical Journal C 75 (Sep, 2015) 412
14	CMS Collaboration	Search for higgs boson pair production in the $ bb {\tau} {\tau} $ final state in proton-proton collisions at $ \sqrt{(}s)=8\text{}\text{}\mathrm{TeV} $	PRD 96 (Oct, 2017) 072004
15	A. M. Sirunyan et al.	Search for resonant and nonresonant Higgs boson pair production in the $ b\bar{b}\ell\nu\ell\nu $ final state in proton-proton collisions at $ \sqrt{s} = $ 13 TeV	JHEP 2018 (Jan, 2018) 54
16	CMS Collaboration	Search for higgs boson pair production in the $ \gamma\gamma\mathrm{b\bar{b}} $ final state in pp collisions at $ \sqrt{s} = $ 13 TeV	CDS
17	A. Falkowski	Effective field theory approach to LHC Higgs data	Pramana 87 (Aug, 2016) 39
18	T. Corbett et al.	The Higgs legacy of the LHC run I	JHEP 08 (2015)	1505.05516
19	A. Carvalho et al.	Higgs pair production: choosing benchmarks with cluster analysis	JHEP 04 (2016) 126
20	CMS Collaboration	The CMS experiment at the CERN LHC	JINST 3 (2008) S08004	CMS-00-001
21	CMS Collaboration	Identification of b quark jets at the CMS Experiment in the LHC Run 2	CMS-PAS-BTV-15-001	CMS-PAS-BTV-15-001
22	B. Hespel, D. Lopez-Val, and E. Vryonidou	Higgs pair production via gluon fusion in the Two-Higgs-Doublet Model	JHEP 1409 (2014) 124	1407.0281
23	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) 079	1405.0301
24	J. Butterworth et al.	PDF4LHC recommendations for LHC Run II	JPG43 (2016) 023001	1510.03865
25	A. Carvalho et al.	On the reinterpretation of non-resonant searches for Higgs boson pairs		1710.08261
26	T. Sjostrand, S. Mrenna, and P. Skands	PYTHIA 6.4 physics and manual	JHEP 05 (2006) 026	hep-ph/0603175
27	T. Sjostrand, S. Mrenna, and P. Skands	A brief introduction to PYTHIA 8.1	Computer Physics Communications 2008 (2008), no. 178, 852
28	M. L. Mangano, M. Moretti, F. Piccinini, and M. Treccani	Matching matrix elements and shower evolution for top-quark production in hadronic collisions	JHEP 2007 (2007), no. 0701, 013
29	S. Alioli, P. Nason, C. Oleari, and E. Re	Erratum: NLO single-top production matched with shower in POWHEG: s- and t-channel contributions	JHEP 2010 (2010), no. 2
30	J. M. Campbell, R. K. Ellis, P. Nason, and E. Re	Top-pair production and decay at NLO matched with parton showers	JHEP 04 (2015)	1412.1828
31	S. Agostinelli et al.	Geant4 - a simulation toolkit	Nucl. Instr. Meth. A506 (2003), no. 3, 250 -- 303
32	CMS Collaboration	Particle-flow reconstruction and global event description with the cms detector	JINST 12 (2017) P10003	CMS-PRF-14-001 1706.04965
33	M. Cacciari, G. P. Salam, and G. Soyez	The anti-$ k_t $ jet clustering algorithm	JHEP 04 (2008) 063	0802.1189
34	M. Cacciari, G. P. Salam, and G. Soyez	Fastjet user manual	The European Physical Journal C 72 (2012), no. 3
35	CMS Collaboration	Determination of jet energy calibration and transverse momentum resolution in CMS	Journal of Instrumentation 6 (2011), no. 11
36	CMS Collaboration	Pileup Jet Identification	CMS-PAS-JME-13-005	CMS-PAS-JME-13-005
37	CMS Collaboration	Performance of b-Tagging Algorithms in Proton Collisions at 13 TeV using the 2016 Data	CDS
38	T. Chen and C. Guestrin	Xgboost: A scalable tree boosting system	in Proceedings of the 22Nd ACM SIGKDD International Conference on Knowledge Discovery and Data Mining, KDD '16 ACM, New York, NY, USA
39	M. Dall'Osso et al.	Hemisphere mixing: A fully data-driven model of the QCD multijet background for the non-resonant higgs pair production search	in The European Physical Society Conference on High Energy Physics (EPS-HEP2017), volume 314 SISSA
40	CMS Collaboration	CMS Luminosity Measurements or the 2016 Data Taking Period	CMS-PAS-LUM-17-001	CMS-PAS-LUM-17-001
41	A. L. Read	Presentation of search results: The CL(s) technique	JPG28 (2002)
42	ATLAS and CMS Collaborations	Procedure for the LHC Higgs boson search combination in summer 2011	CMS-NOTE-2011-005