CMS-PAS-HIG-19-018 | ||
Search for nonresonant Higgs boson pair production in final states with two bottom quarks and two photons in pp collisions at $\sqrt{s}=$ 13 TeV | ||
CMS Collaboration | ||
October 2020 | ||
Abstract: A search for nonresonant production of Higgs boson pairs via gluon-gluon fusion and vector boson fusion in final states with two bottom quarks and two photons is presented. This search uses data from proton-proton collisions at a center-of-mass energy of $\sqrt{s}=$ 13 TeV recorded by the CMS detector at the LHC from 2016 to 2018, corresponding to an integrated luminosity of 137 fb$^{-1}$. No signal is observed, and a 95% confidence level upper limit is set on the product of the inclusive Higgs boson pair production cross section and branching fraction into $\gamma\gamma\mathrm{b\bar{b}}$. The observed (expected) upper limit is determined to be 0.67 (0.45) fb, which corresponds to 7.7 (5.2) times the standard model prediction. Assuming all other Higgs boson couplings are equal to their values in the standard model, the coupling modifiers of the trilinear self-coupling $\kappa_{\lambda}$ and the coupling between a pair of Higgs bosons and a pair of vector bosons $c_{2V}$ are constrained within the ranges $-3.3 < \kappa_{\lambda} < 8.5$ and $-1.3 < c_{2V} < 3.5$ at 95% confidence level. Constraints on $\kappa_{\lambda}$ are also set by combining this analysis with a search for single Higgs bosons produced in association with top quark-antiquark pairs, and by performing a simultaneous fit of $\kappa_{\lambda}$ and the top Yukawa coupling modifier $\kappa_{\mathrm{t}}$. | ||
Links:
CDS record (PDF) ;
CADI line (restricted) ;
These preliminary results are superseded in this paper, JHEP 03 (2021) 257. The superseded preliminary plots can be found here. |
Figures | |
png pdf |
Figure 1:
Feynman diagrams of the processes contributing to the production of Higgs boson pairs via ggF at LO. The top diagrams correspond to SM processes, involving the top Yukawa coupling ${y_{\mathrm{t}}}$ and the trilinear Higgs coupling ${\lambda _{\mathrm{H} \mathrm{H} \mathrm{H}}}$, respectively. The bottom diagrams correspond to BSM processes: the diagram on the left involves the contact interaction of two Higgs bosons with two top quarks (${c_2}$), the middle diagram shows the quartic coupling between the Higgs bosons and two gluons (${c_{2\mathrm{g}}}$), and the diagram on the right describes the contact interactions between the Higgs boson and gluons (${c_{\mathrm{g}}}$). |
png pdf |
Figure 1-a:
Feynman diagram of a SM process contributing to the production of Higgs boson pairs via ggF at LO, involving the top Yukawa coupling ${y_{\mathrm{t}}}$ and the trilinear Higgs coupling ${\lambda _{\mathrm{H} \mathrm{H} \mathrm{H}}}$. |
png pdf |
Figure 1-b:
Feynman diagram of a SM process contributing to the production of Higgs boson pairs via ggF at LO, involving the top Yukawa coupling ${y_{\mathrm{t}}}$ and the trilinear Higgs coupling ${\lambda _{\mathrm{H} \mathrm{H} \mathrm{H}}}$. |
png pdf |
Figure 1-c:
Feynman diagram of a BSM process contributing to the production of Higgs boson pairs via ggF at LO. This process involves the contact interaction of two Higgs bosons with two top quarks (${c_2}$). |
png pdf |
Figure 1-d:
Feynman diagram of a BSM process contributing to the production of Higgs boson pairs via ggF at LO. This process involves the quartic coupling between the Higgs bosons and two gluons (${c_{2\mathrm{g}}}$). |
png pdf |
Figure 1-e:
Feynman diagram of a BSM process contributing to the production of Higgs boson pairs via ggF at LO. This process involves the contact interactions between the Higgs boson and gluons (${c_{\mathrm{g}}}$) and the trilinear Higgs coupling ${\lambda _{\mathrm{H} \mathrm{H} \mathrm{H}}}$. |
png pdf |
Figure 2:
The Feynman diagrams that contribute to the production of Higgs boson pairs via VBF at LO. On the left the diagram involving the HHH vertex (${\lambda _{\mathrm{H} \mathrm{H} \mathrm{H}}}$), in the middle the diagram with two HVV vertices ($ {c_{V}} $), and on the right the diagram with the HHVV vertex ($ {c_{2V}} $). |
png pdf |
Figure 2-a:
A Feynman diagram that contributes to the production of Higgs boson pairs via VBF at LO. The diagram involves the HHH vertex (${\lambda _{\mathrm{H} \mathrm{H} \mathrm{H}}}$). |
png pdf |
Figure 2-b:
A Feynman diagram that contributes to the production of Higgs boson pairs via VBF at LO. The diagram involves two HVV vertices ($ {c_{V}} $). |
png pdf |
Figure 2-c:
A Feynman diagram that contributes to the production of Higgs boson pairs via VBF at LO. The diagram involves the HHVV vertex ($ {c_{2V}} $). |
png pdf |
Figure 3:
The distribution of $ {m_{\gamma \gamma}} $ (left) and $ {m_\text {jj}}$ (right) in data and simulated events. Data, dominated by the $ {\gamma \gamma}$+jets and $ {\gamma}$+jets backgrounds, are compared to the SM ggF HH signal samples and single H samples (ttH, ggH, VBF-H, VH) after requiring the selection criteria described in Section 5. The error bars on the data points indicate statistical uncertainties. The HH signal has been scaled by a factor of 10$^{3}$ for display purposes. |
png pdf |
Figure 3-a:
The distribution of $ {m_{\gamma \gamma}} $ in data and simulated events. Data, dominated by the $ {\gamma \gamma}$+jets and $ {\gamma}$+jets backgrounds, are compared to the SM ggF HH signal samples and single H samples (ttH, ggH, VBF-H, VH) after requiring the selection criteria described in Section 5. The error bars on the data points indicate statistical uncertainties. The HH signal has been scaled by a factor of 10$^{3}$ for display purposes. |
png pdf |
Figure 3-b:
The distribution of $ {m_\text {jj}}$ in data and simulated events. Data, dominated by the $ {\gamma \gamma}$+jets and $ {\gamma}$+jets backgrounds, are compared to the SM ggF HH signal samples and single H samples (ttH, ggH, VBF-H, VH) after requiring the selection criteria described in Section 5. The error bars on the data points indicate statistical uncertainties. The HH signal has been scaled by a factor of 10$^{3}$ for display purposes. |
png pdf |
Figure 4:
Distributions of $ {\tilde{M}_{\mathrm {X}}} $. The SM ggF HH signal is compared with several BSM hypotheses listed in Table 1 (left), and the SM VBF HH signal is compared with different anomalous values of $ {c_{2V}} $ (right). All distributions are normalized to unity. |
png pdf |
Figure 4-a:
Distribution of $ {\tilde{M}_{\mathrm {X}}} $. The SM ggF HH signal is compared with several BSM hypotheses listed in Table 1. All distributions are normalized to unity. |
png pdf |
Figure 4-b:
Distribution of $ {\tilde{M}_{\mathrm {X}}} $. The SM VBF HH signal is compared with different anomalous values of $ {c_{2V}} $. All distributions are normalized to unity. |
png pdf |
Figure 5:
The distribution of the ttHScore (left) and MVA output (right) in data and simulated events. Data, dominated by $ {\gamma \gamma}$+jets and $ {\gamma}$+jets background, are compared to the SM ggF HH signal samples and single H samples (ttH, ggH, VBF-H, VH) after requiring the selection criteria described in Section 5. The error bars on the data points indicate statistical uncertainties. The HH signal has been scaled by a factor of 10$^{3}$ for display purposes. |
png pdf |
Figure 5-a:
The distribution of the ttHScore in data and simulated events. Data, dominated by $ {\gamma \gamma}$+jets and $ {\gamma}$+jets background, are compared to the SM ggF HH signal samples and single H samples (ttH, ggH, VBF-H, VH) after requiring the selection criteria described in Section 5. The error bars on the data points indicate statistical uncertainties. The HH signal has been scaled by a factor of 10$^{3}$ for display purposes. |
png pdf |
Figure 5-b:
The distribution of the MVA output in data and simulated events. Data, dominated by $ {\gamma \gamma}$+jets and $ {\gamma}$+jets background, are compared to the SM ggF HH signal samples and single H samples (ttH, ggH, VBF-H, VH) after requiring the selection criteria described in Section 5. The error bars on the data points indicate statistical uncertainties. The HH signal has been scaled by a factor of 10$^{3}$ for display purposes. |
png pdf |
Figure 6:
The distribution of the two MVA outputs is shown in data and simulated events in the two VBF $ {\tilde{M}_{\mathrm {X}}} $ regions: $ {\tilde{M}_{\mathrm {X}}} > $ 500 GeV (left) and $ {\tilde{M}_{\mathrm {X}}} < $ 500 GeV (right) regions. Data, dominated by the $ {\gamma \gamma}$+jets and $ {\gamma}$+jets backgrounds, are compared to the VBF HH signal samples with SM couplings and $ {c_{2V}} =$ 0, and single H samples (ttH, ggH, VBF-H, VH) after requiring the VBF selection criteria described in Section 5. The error bars on the data points indicate statistical uncertainties. The HH signal has been scaled by a factor of 10$^{3}$ for display purposes. |
png pdf |
Figure 6-a:
The distribution of the two MVA outputs is shown in data and simulated events in the VBF $ {\tilde{M}_{\mathrm {X}}} > $ 500 GeV region. Data, dominated by the $ {\gamma \gamma}$+jets and $ {\gamma}$+jets backgrounds, are compared to the VBF HH signal samples with SM couplings and $ {c_{2V}} =$ 0, and single H samples (ttH, ggH, VBF-H, VH) after requiring the VBF selection criteria described in Section 5. The error bars on the data points indicate statistical uncertainties. The HH signal has been scaled by a factor of 10$^{3}$ for display purposes. |
png pdf |
Figure 6-b:
The distribution of the two MVA outputs is shown in data and simulated events in the VBF $ {\tilde{M}_{\mathrm {X}}} < $ 500 GeV region. Data, dominated by the $ {\gamma \gamma}$+jets and $ {\gamma}$+jets backgrounds, are compared to the VBF HH signal samples with SM couplings and $ {c_{2V}} =$ 0, and single H samples (ttH, ggH, VBF-H, VH) after requiring the VBF selection criteria described in Section 5. The error bars on the data points indicate statistical uncertainties. The HH signal has been scaled by a factor of 10$^{3}$ for display purposes. |
png pdf |
Figure 7:
Parametrized signal shape for $ {m_{\gamma \gamma}} $ (left) and ${m_\text {jj}}$ (right) in the best resolution ggF (top) and VBF (bottom) categories. The open squares represent weighted simulated events and the blue lines are the corresponding models. Also shown are the $\sigma _{\text {eff}}$ value (half the width of the narrowest interval containing 68.3% of the invariant mass distribution) and the corresponding interval as a gray band, and the full width at half the maximum (FWHM) and the corresponding interval as a double arrow. |
png pdf |
Figure 7-a:
Parametrized signal shape for $ {m_{\gamma \gamma}} $ in the best resolution ggF category. The open squares represent weighted simulated events and the blue lines are the corresponding models. Also shown are the $\sigma _{\text {eff}}$ value (half the width of the narrowest interval containing 68.3% of the invariant mass distribution) and the corresponding interval as a gray band, and the full width at half the maximum (FWHM) and the corresponding interval as a double arrow. |
png pdf |
Figure 7-b:
Parametrized signal shape for ${m_\text {jj}}$ in the best resolution ggF category. The open squares represent weighted simulated events and the blue lines are the corresponding models. Also shown are the $\sigma _{\text {eff}}$ value (half the width of the narrowest interval containing 68.3% of the invariant mass distribution) and the corresponding interval as a gray band, and the full width at half the maximum (FWHM) and the corresponding interval as a double arrow. |
png pdf |
Figure 7-c:
Parametrized signal shape for $ {m_{\gamma \gamma}} $ in the best resolution VBF category. The open squares represent weighted simulated events and the blue lines are the corresponding models. Also shown are the $\sigma _{\text {eff}}$ value (half the width of the narrowest interval containing 68.3% of the invariant mass distribution) and the corresponding interval as a gray band, and the full width at half the maximum (FWHM) and the corresponding interval as a double arrow. |
png pdf |
Figure 7-d:
Parametrized signal shape for ${m_\text {jj}}$ in the best resolution VBF category. The open squares represent weighted simulated events and the blue lines are the corresponding models. Also shown are the $\sigma _{\text {eff}}$ value (half the width of the narrowest interval containing 68.3% of the invariant mass distribution) and the corresponding interval as a gray band, and the full width at half the maximum (FWHM) and the corresponding interval as a double arrow. |
png pdf |
Figure 8:
Invariant mass distributions $ {m_{\gamma \gamma}} $ (top row) and $ {m_\text {jj}}$ (bottom row) for the selected events in data (black points) in the best resolution ggF (CAT0) and VBF (CAT0) categories are shown. The solid red line shows the sum of the fitted signal and background, the solid blue line shows the background component from the single Higgs boson and the nonresonant processes, and the dashed black line shows the nonresonant background component. The one (green) and two (yellow) standard deviation bands include the uncertainties in the background component of the fit. The lower panel in each plot shows the residuals after the background subtraction. |
png pdf |
Figure 8-a:
The invariant mass distribution $ {m_{\gamma \gamma}} $ for the selected events in data (black points) in the best resolution ggF (CAT0) category is shown. The solid red line shows the sum of the fitted signal and background, the solid blue line shows the background component from the single Higgs boson and the nonresonant processes, and the dashed black line shows the nonresonant background component. The one (green) and two (yellow) standard deviation bands include the uncertainties in the background component of the fit. The lower panel in each plot shows the residuals after the background subtraction. |
png pdf |
Figure 8-b:
The invariant mass distribution $ {m_{\gamma \gamma}} $ for the selected events in data (black points) in the best resolution VBF (CAT0) category is shown. The solid red line shows the sum of the fitted signal and background, the solid blue line shows the background component from the single Higgs boson and the nonresonant processes, and the dashed black line shows the nonresonant background component. The one (green) and two (yellow) standard deviation bands include the uncertainties in the background component of the fit. The lower panel in each plot shows the residuals after the background subtraction. |
png pdf |
Figure 8-c:
The invariant mass distribution $ {m_\text {jj}}$ for the selected events in data (black points) in the best resolution ggF (CAT0) category is shown. The solid red line shows the sum of the fitted signal and background, the solid blue line shows the background component from the single Higgs boson and the nonresonant processes, and the dashed black line shows the nonresonant background component. The one (green) and two (yellow) standard deviation bands include the uncertainties in the background component of the fit. The lower panel in each plot shows the residuals after the background subtraction. |
png pdf |
Figure 8-d:
The invariant mass distribution $ {m_\text {jj}}$ for the selected events in data (black points) in the best resolution VBF (CAT0) category is shown. The solid red line shows the sum of the fitted signal and background, the solid blue line shows the background component from the single Higgs boson and the nonresonant processes, and the dashed black line shows the nonresonant background component. The one (green) and two (yellow) standard deviation bands include the uncertainties in the background component of the fit. The lower panel in each plot shows the residuals after the background subtraction. |
png pdf |
Figure 9:
Invariant mass distribution $ {m_{\gamma \gamma}} $ (left) and $ {m_\text {jj}}$ (right) for the selected events in data (black points) weighted by S/(S + B), where S (B) is the number of expected signal (background) events in a $\pm $1$\sigma _{\mathrm {eff}}$ mass window centered on ${m_{\mathrm{H}}}$. The variable $\sigma _{\mathrm {eff}}$ is defined as the smallest interval containing 68.3% of the distribution. The solid red line shows the sum of the fitted signal and background, the solid blue line shows the background component from the single Higgs boson and the nonresonant processes, and the dashed black line shows the nonresonant background component. The one (green) and two (yellow) standard deviation bands include the uncertainties in the background component of the fit. The lower panel shows the residuals after the background subtraction. |
png pdf |
Figure 9-a:
Invariant mass distribution $ {m_{\gamma \gamma}} $ for the selected events in data (black points) weighted by S/(S + B), where S (B) is the number of expected signal (background) events in a $\pm $1$\sigma _{\mathrm {eff}}$ mass window centered on ${m_{\mathrm{H}}}$. The variable $\sigma _{\mathrm {eff}}$ is defined as the smallest interval containing 68.3% of the distribution. The solid red line shows the sum of the fitted signal and background, the solid blue line shows the background component from the single Higgs boson and the nonresonant processes, and the dashed black line shows the nonresonant background component. The one (green) and two (yellow) standard deviation bands include the uncertainties in the background component of the fit. The lower panel shows the residuals after the background subtraction. |
png pdf |
Figure 9-b:
Invariant mass distribution $ {m_\text {jj}}$ for the selected events in data (black points) weighted by S/(S + B), where S (B) is the number of expected signal (background) events in a $\pm $1$\sigma _{\mathrm {eff}}$ mass window centered on ${m_{\mathrm{H}}}$. The variable $\sigma _{\mathrm {eff}}$ is defined as the smallest interval containing 68.3% of the distribution. The solid red line shows the sum of the fitted signal and background, the solid blue line shows the background component from the single Higgs boson and the nonresonant processes, and the dashed black line shows the nonresonant background component. The one (green) and two (yellow) standard deviation bands include the uncertainties in the background component of the fit. The lower panel shows the residuals after the background subtraction. |
png pdf |
Figure 10:
Expected and observed 95% CL upper limits on the product of the HH production cross section and $\mathcal {B}({{\mathrm{H} \mathrm{H}}} \to {\gamma \gamma \mathrm{b} {}\mathrm{\bar{b}}})$ obtained for different values of $ {\kappa _{\lambda}} $ assuming $ {\kappa _{\mathrm{t}}} = $1. The green and yellow bands represent, respectively, the one and two standard deviation extensions beyond the expected limit. The red line shows the theoretical prediction. |
png pdf |
Figure 11:
Negative log-likelihood as a function of $ {\kappa _{\lambda}} $ evaluated with an Asimov data set assuming the SM hypothesis (left) and the observed data (right) are shown. The 68 and 98% CL intervals are shown with the dashed gray lines. The two curves are shown for the HH (blue) and HH+ttH (orange) analysis categories. All other couplings are set to their SM values. |
png pdf |
Figure 11-a:
Negative log-likelihood as a function of $ {\kappa _{\lambda}} $ evaluated with an Asimov data set assuming the SM hypothesis is shown. The 68 and 98% CL intervals are shown with the dashed gray lines. The two curves are shown for the HH (blue) and HH+ttH (orange) analysis categories. All other couplings are set to their SM values. |
png pdf |
Figure 11-b:
Negative log-likelihood as a function of $ {\kappa _{\lambda}} $ evaluated with an Asimov data set assuming the observed data is shown. The 68 and 98% CL intervals are shown with the dashed gray lines. The two curves are shown for the HH (blue) and HH+ttH (orange) analysis categories. All other couplings are set to their SM values. |
png pdf |
Figure 12:
Negative log-likelihood contours at 68% and 95% CL in the ($ {\kappa _{\lambda}} $, $ {\kappa _{\mathrm{t}}} $) plane evaluated with an Asimov data set assuming the SM hypothesis (left) and the observed data (right). The contours obtained using the HH analysis categories only are shown in blue, and in orange when combined with the ttH categories. The best fit value for the HH categories only ($ {\kappa _{\lambda}} = $ 0.6, $ {\kappa _{\mathrm{t}}} = $ 1.2) is indicated by a blue circle, for the HH+ttH categories ($ {\kappa _{\lambda}} = $ 1.4, $ {\kappa _{\mathrm{t}}} = $ 1.3) by a orange diamond, and the SM prediction ($ {\kappa _{\lambda}} = $ 1.0, $ {\kappa _{\mathrm{t}}} = $ 1.0) by a black star. The regions of the 2D scan where the ${\kappa _{\mathrm{t}}}$ parametrization for anomalous values of ${\kappa _{\lambda}}$ at LO is not reliable are shown with a gray band. |
png pdf |
Figure 12-a:
Negative log-likelihood contours at 68% and 95% CL in the ($ {\kappa _{\lambda}} $, $ {\kappa _{\mathrm{t}}} $) plane evaluated with an Asimov data set assuming the SM hypothesis. The contours obtained using the HH analysis categories only are shown in blue, and in orange when combined with the ttH categories. The best fit value for the HH categories only ($ {\kappa _{\lambda}} = $ 0.6, $ {\kappa _{\mathrm{t}}} = $ 1.2) is indicated by a blue circle, for the HH+ttH categories ($ {\kappa _{\lambda}} = $ 1.4, $ {\kappa _{\mathrm{t}}} = $ 1.3) by a orange diamond, and the SM prediction ($ {\kappa _{\lambda}} = $ 1.0, $ {\kappa _{\mathrm{t}}} = $ 1.0) by a black star. The regions of the 2D scan where the ${\kappa _{\mathrm{t}}}$ parametrization for anomalous values of ${\kappa _{\lambda}}$ at LO is not reliable are shown with a gray band. |
png pdf |
Figure 12-b:
Negative log-likelihood contours at 68% and 95% CL in the ($ {\kappa _{\lambda}} $, $ {\kappa _{\mathrm{t}}} $) plane evaluated with an Asimov data set assuming the observed data. The contours obtained using the HH analysis categories only are shown in blue, and in orange when combined with the ttH categories. The best fit value for the HH categories only ($ {\kappa _{\lambda}} = $ 0.6, $ {\kappa _{\mathrm{t}}} = $ 1.2) is indicated by a blue circle, for the HH+ttH categories ($ {\kappa _{\lambda}} = $ 1.4, $ {\kappa _{\mathrm{t}}} = $ 1.3) by a orange diamond, and the SM prediction ($ {\kappa _{\lambda}} = $ 1.0, $ {\kappa _{\mathrm{t}}} = $ 1.0) by a black star. The regions of the 2D scan where the ${\kappa _{\mathrm{t}}}$ parametrization for anomalous values of ${\kappa _{\lambda}}$ at LO is not reliable are shown with a gray band. |
png pdf |
Figure 13:
Negative log-likelihood scan as a function of $ {\kappa _{\mathrm{t}}} $ evaluated with an Asimov data set assuming the SM hypothesis (left) and the observed data (right) are shown. The 68 and 98% CL intervals are shown with the dashed gray lines. The two curves are shown for the HH (blue) and the HH+ttH (orange) analysis categories. All other couplings are fixed to their SM values. |
png pdf |
Figure 13-a:
Negative log-likelihood scan as a function of $ {\kappa _{\mathrm{t}}} $ evaluated with an Asimov data set assuming the SM hypothesis is shown. The 68 and 98% CL intervals are shown with the dashed gray lines. The two curves are shown for the HH (blue) and the HH+ttH (orange) analysis categories. All other couplings are fixed to their SM values. |
png pdf |
Figure 13-b:
Negative log-likelihood scan as a function of $ {\kappa _{\mathrm{t}}} $ evaluated with an Asimov data set assuming the observed data is shown. The 68 and 98% CL intervals are shown with the dashed gray lines. The two curves are shown for the HH (blue) and the HH+ttH (orange) analysis categories. All other couplings are fixed to their SM values. |
png pdf |
Figure 14:
Expected and observed 95% CL upper limits on the product of the VBF HH production cross section and $\mathcal {B}({{\mathrm{H} \mathrm{H}}} \to {\gamma \gamma \mathrm{b} {}\mathrm{\bar{b}}})$ obtained for different values of $ {c_{2V}} $. The green and yellow bands represent, respectively, the one and two standard deviation extensions beyond the expected limit. The red line shows the theoretical prediction. |
png pdf |
Figure 15:
Negative log-likelihood contours at 68% and 95% CL in the ($ {\kappa _{\lambda}} $, $ {c_{2V}} $) plane evaluated with an Asimov data set assuming the SM hypothesis (left) and with the observed data (right). The contours are obtained using the HH analysis categories only. The best fit value ($ {\kappa _{\lambda}} = $ 0.0, $ {c_{2V}} = $ 0.3) is indicated by a blue circle, and the SM prediction ($ {\kappa _{\lambda}} = $ 1.0, $ {c_{2V}} = $ 1.0) by a black star. |
png pdf |
Figure 15-a:
Negative log-likelihood contours at 68% and 95% CL in the ($ {\kappa _{\lambda}} $, $ {c_{2V}} $) plane evaluated with an Asimov data set assuming the SM hypothesis. The contours are obtained using the HH analysis categories only. The best fit value ($ {\kappa _{\lambda}} = $ 0.0, $ {c_{2V}} = $ 0.3) is indicated by a blue circle, and the SM prediction ($ {\kappa _{\lambda}} = $ 1.0, $ {c_{2V}} = $ 1.0) by a black star. |
png pdf |
Figure 15-b:
Negative log-likelihood contours at 68% and 95% CL in the ($ {\kappa _{\lambda}} $, $ {c_{2V}} $) plane evaluated with an Asimov data set assuming the observed data. The contours are obtained using the HH analysis categories only. The best fit value ($ {\kappa _{\lambda}} = $ 0.0, $ {c_{2V}} = $ 0.3) is indicated by a blue circle, and the SM prediction ($ {\kappa _{\lambda}} = $ 1.0, $ {c_{2V}} = $ 1.0) by a black star. |
png pdf |
Figure 16:
Expected and observed 95% CL upper limits on the product of the ggF HH production cross section and $\mathcal {B}({{\mathrm{H} \mathrm{H}}} \to {\gamma \gamma \mathrm{b} {}\mathrm{\bar{b}}})$ obtained for different nonresonant benchmark models (defined in Table 1). The green and yellow bands represent, respectively, the one and two standard deviation extensions beyond the expected limit. |
png pdf |
Figure 17:
Expected and observed 95% CL upper limits on the product of the ggF HH production cross section and $\mathcal {B}({{\mathrm{H} \mathrm{H}}} \to {\gamma \gamma \mathrm{b} {}\mathrm{\bar{b}}})$ obtained for different values of the BSM coupling $ {c_2} $. The green and yellow bands represent, respectively, the one and two standard deviation extensions beyond the expected limit. The red line shows the theoretical prediction. |
Tables | |
png pdf |
Table 1:
Coupling parameter values in the SM and in twelve BSM benchmark hypotheses identified using the method described in Ref. [27]. |
png pdf |
Table 2:
Summary of the analysis categories. Two VBF and 12 ggF-enriched categories are defined based on the output of the MVA classifier and the mass of the Higgs boson pair system $ {\tilde{M}_{\mathrm {X}}} $. The VBF and ggF categories are mutually exclusive. |
Summary |
A search for nonresonant Higgs boson pair production (HH) has been presented, where one of the Higgs bosons decays to a pair of bottom quarks and the other to a pair of photons. This search uses proton-proton collision data collected at $\sqrt{s} = $ 13 TeV by the CMS experiment at the LHC, corresponding to a total integrated luminosity of 137 fb$^{-1}$. No signal has been observed. Upper limits at 95% confidence level (CL) on the product of the HH production cross section and the branching fraction into ${\mathrm{b\bar{b}}} gg$ are extracted for production in the standard model (SM) and in several scenarios beyond the standard model (BSM). The expected upper limit on ${\sigma_{\mathrm{H}\mathrm{H}}} {\mathcal{B}({{\mathrm{H}\mathrm{H}}} \to\mathrm{b\bar{b}}\gamma\gamma)}$ is 0.45 fb, corresponding to about 5.2 times the SM prediction, while the observed upper limit is 0.67 fb, corresponding to 7.7 times the expected value for the SM process. The presented result has the highest sensitivity to the SM HH production to date. Upper limits at 95% CL on the SM HH production cross section are also derived as a function of the Higgs boson self-coupling modifier ${\kappa_{\lambda}} \equiv {\lambda_{\mathrm{H}\mathrm{H}\mathrm{H}}} /{\lambda^\mathrm{SM}_{\mathrm{H}\mathrm{H}\mathrm{H}}} $ assuming that the top quark Yukawa coupling is SM-like. The coupling modifier ${\kappa_{\lambda}} $ is constrained within a range $-3.3 < \kappa_{\lambda} < 8.5$, while the expected constraint is within a range $-2.5 < \kappa_{\lambda} < 8.2$ at 95% CL. This search is combined with an analysis that targets top quark-antiquark associated production of a single Higgs boson decaying to a diphoton pair. In the scenario in which the HH signal has the properties predicted by the SM, the coupling modifier ${\kappa_{\lambda}} $ has been constrained. In addition, a simultaneous measurement of ${\kappa_{\lambda}} $ and the modifier of the coupling between the Higgs boson and the top quark ${\kappa_{\mathrm{t}}}$ is presented when both the HH and single Higgs boson processes are considered as signals. Limits are also set on the cross section of nonresonant HH production via vector boson fusion. The most stringent limit to date is set on the product of the vector boson fusion HH production cross section and the branching fraction into ${\mathrm{b\bar{b}}} \gamma\gamma$. The observed (expected) upper limit at 95% CL amounts to 1.02 (0.94) fb, corresponding to 225 (208) times the SM prediction. Limits are also set as a function of the modifier of the coupling between two vector bosons and two Higgs bosons, ${c_2} v$. The observed excluded region corresponds to ${c_2} v < -1.3$ and ${c_2} v > 3.5$, while the expected exclusion is ${c_2} v < -0.9$ and ${c_2} v > 3.1$. Numerous BSM hypotheses and coupling modifiers have been explored, both in the context of inclusive Higgs boson pair production and for HH production via gluon-gluon fusion and vector boson fusion. The production of Higgs boson pairs was also combined with the top quark-antiquark pair associated production of a single Higgs boson. Overall, all of the results are consistent with the SM predictions. |
References | ||||
1 | ATLAS Collaboration | Observation of a new particle in the search for the standard model Higgs boson with the ATLAS 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 | CMS Collaboration | Observation of a new boson with mass near 125 GeV in pp collisions at $ \sqrt{s} = $ 7 and 8 TeV | JHEP 06 (2013) 081 | CMS-HIG-12-036 1303.4571 |
4 | F. Englert and R. Brout | Broken symmetries and the masses of gauge bosons | PRL 13 (1964) 321 | |
5 | P. W. Higgs | Broken symmetries and the masses of gauge bosons | PRL 13 (1964) 508 | |
6 | M. Grazzini et al. | Higgs boson pair production at NNLO with top quark mass effects | JHEP 05 (2018) 059 | 1803.02463 |
7 | F. A. Dreyer and A. Karlberg | Vector-Boson Fusion Higgs Pair Production at N$ ^3 $LO | PRD 98 (2018) 114016 | 1811.07906 |
8 | LHC Higgs Cross Section Working Group Collaboration | Handbook of LHC Higgs Cross Sections: 4. Deciphering the Nature of the Higgs Sector | FERMILAB-FN-1025-T, CERN-2017-002-M, CERN | 1610.07922 |
9 | ATLAS Collaboration | Search For Higgs Boson Pair Production in the $ \gamma\gamma b\bar{b} $ Final State using $ pp $ Collision Data at $ \sqrt{s}= $ 8 TeV from the ATLAS Detector | PRL 114 (2015) 081802 | 1406.5053 |
10 | ATLAS Collaboration | 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 | EPJC 75 (2015) 412 | 1506.00285 |
11 | 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 | PRD 92 (2015) 092004 | 1509.04670 |
12 | ATLAS Collaboration | Search for pair production of Higgs bosons in the $ \mathrm{b\bar{b}}\mathrm{b\bar{b}} $ final state using proton--proton collisions at $ \sqrt{s} = $ 13 TeV with the ATLAS detector | PRD 94 (2016) 052002 | 1606.04782 |
13 | ATLAS Collaboration | Combination of searches for Higgs boson pairs in $ pp $ collisions at $ \sqrt{s} = $ 13 TeV with the ATLAS detector | PLB 800 (2020) 135103 | 1906.02025 |
14 | CMS Collaboration | Search for two Higgs bosons in final states containing two photons and two bottom quarks in proton-proton collisions at 8 TeV | PRD 94 (2016) 052012 | CMS-HIG-13-032 1603.06896 |
15 | ATLAS Collaboration | Search for Higgs boson pair production in the $ \gamma\gamma b\bar{b} $ final state with 13 $ TeV pp $ collision data collected by the ATLAS experiment | JHEP 11 (2018) 040 | 1807.04873 |
16 | CMS Collaboration | Search for Higgs boson pair production in the $ \mathrm{b\bar{b}}\tau\tau $ final state in proton-proton collisions at $ \sqrt{s}= $ 8 TeV | PRD 96 (2017) 072004 | CMS-HIG-15-013 1707.00350 |
17 | 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 = $ 13TeV | PLB 778 (2018) 101--127 | CMS-HIG-17-002 1707.02909 |
18 | CMS Collaboration | Search for resonant and nonresonant Higgs boson pair production in the $ \mathrm{b\bar{b}}\ell\nu\ell\nu $ final state in proton-proton collisions at $ \sqrt{s}= $ 13 TeV | JHEP 01 (2018) 054 | CMS-HIG-17-006 1708.04188 |
19 | CMS Collaboration | Search for Higgs boson pair production in the $ \gamma\gamma\mathrm{b\overline{b}} $ final state in pp collisions at $ \sqrt{s}= $ 13 TeV | PLB 788 (2019) 7 | CMS-HIG-17-008 1806.00408 |
20 | CMS Collaboration | Combination of searches for Higgs boson pair production in proton-proton collisions at $ \sqrt{s} = $ 13 TeV | PRL 122 (2019) 121803 | CMS-HIG-17-030 1811.09689 |
21 | ATLAS Collaboration | Search for the $ HH \rightarrow b \bar{b} b \bar{b} $ process via vector-boson fusion production using proton-proton collisions at $ \sqrt{s} = $ 13 TeV with the ATLAS detector | JHEP 07 (2020) 108 | 2001.05178 |
22 | CMS Collaboration | Measurements of $ \mathrm{t\bar{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 |
23 | F. Maltoni, D. Pagani, A. Shivaji, and X. Zhao | Trilinear Higgs coupling determination via single-Higgs differential measurements at the LHC | EPJC 77 (2017) 887 | 1709.08649 |
24 | CMS Collaboration | The CMS trigger system | JINST 12 (2017) P01020 | CMS-TRG-12-001 1609.02366 |
25 | CMS Collaboration | The CMS experiment at the CERN LHC | JINST 3 (2008) S08004 | CMS-00-001 |
26 | G. Giudice, C. Grojean, A. Pomarol, and R. Rattazzi | The Strongly-Interacting Light Higgs | JHEP 06 (2007) 045 | hep-ph/0703164 |
27 | A. Carvalho et al. | Higgs Pair Production: Choosing Benchmarks With Cluster Analysis | JHEP 04 (2016) 126 | 1507.02245 |
28 | CMS Collaboration | Measurements of Higgs boson properties in the diphoton decay channel in proton-proton collisions at $ \sqrt{s} = $ 13 TeV | JHEP 11 (2018) 185 | CMS-HIG-16-040 1804.02716 |
29 | 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 |
30 | 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 |
31 | 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 |
32 | 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 |
33 | B. Hespel, D. Lopez-Val, and E. Vryonidou | Higgs pair production via gluon fusion in the two-Higgs-doublet model | JHEP 09 (2014) 124 | 1407.0281 |
34 | R. Frederix et al. | Higgs pair production at the LHC with NLO and parton-shower effects | PLB 732 (2014) 142 | 1401.7340 |
35 | A. Buckley et al. | LHAPDF6: parton density access in the LHC precision era | EPJC 75 (2015) 132 | 1412.7420 |
36 | S. Carrazza, J. I. Latorre, J. Rojo, and G. Watt | A compression algorithm for the combination of PDF sets | EPJC 75 (2015) 474 | 1504.06469 |
37 | J. Butterworth et al. | PDF4LHC recommendations for LHC Run II | JPG 43 (2016) 023001 | 1510.03865 |
38 | S. Dulat et al. | New parton distribution functions from a global analysis of quantum chromodynamics | PRD 93 (2016) 033006 | 1506.07443 |
39 | 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) 204 | 1412.3989 |
40 | NNPDF Collaboration | Parton distributions for the LHC Run II | JHEP 04 (2015) 040 | 1410.8849 |
41 | T. Gleisberg, S. Hoeche, F. Krauss, M. Schonherr, S. Schumann, F. Siegert, and J. Winter | Event generation with SHERPA 1.1 | JHEP 02 (2009) 007 | 0811.4622 |
42 | T. Sjostrand et al. | An introduction to PYTHIA 8.2 | CPC 191 (2015) 159 | 1410.3012 |
43 | P. Nason | A New method for combining NLO QCD with shower Monte Carlo algorithms | JHEP 11 (2004) 040 | hep-ph/0409146 |
44 | S. Frixione, P. Nason, and C. Oleari | Matching NLO QCD computations with Parton Shower simulations: the POWHEG method | JHEP 11 (2007) 070 | 0709.2092 |
45 | 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 |
46 | CMS Collaboration | Event generator tunes obtained from underlying event and multiparton scattering measurements | EPJC 76 (2016) 155 | CMS-GEN-14-001 1512.00815 |
47 | 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 |
48 | NNPDF Collaboration | Parton distributions from high-precision collider data | EPJC 77 (2017) 663 | 1706.00428 |
49 | GEANT4 Collaboration | GEANT4---a simulation toolkit | NIMA 506 (2003) 250 | |
50 | B. Cabouat and T. Sjostrand | Some Dipole Shower Studies | EPJC 78 (2018) 226 | 1710.00391 |
51 | B. Jager et al. | Parton-shower effects in Higgs production via Vector-Boson Fusion | EPJC 80 (2020) 756 | 2003.12435 |
52 | CMS Collaboration | Observation of the diphoton decay of the Higgs boson and measurement of its properties | EPJC 74 (2014) 3076 | CMS-HIG-13-001 1407.0558 |
53 | CMS Collaboration | Particle-flow reconstruction and global event description with the CMS detector | JINST 12 (2017) P10003 | CMS-PRF-14-001 1706.04965 |
54 | M. Cacciari, G. P. Salam, and G. Soyez | The anti-$ {k_{\mathrm{T}}} $ jet clustering algorithm | JHEP 04 (2008) 063 | 0802.1189 |
55 | M. Cacciari, G. P. Salam, and G. Soyez | FastJet user manual | EPJC 72 (2012) 1896 | 1111.6097 |
56 | CMS Collaboration, H. Weber | The phase-1 upgrade of the CMS pixel detector | in Proceedings, 2016 IEEE Nuclear Science Symposium and Medical Imaging Conference: NSS/MIC 2016: Strasbourg, France, p. 8069719 2016 | |
57 | CMS Collaboration | Identification of heavy-flavour jets with the CMS detector in pp collisions at 13 TeV | JINST 13 (2018) P05011 | CMS-BTV-16-002 1712.07158 |
58 | CMS Collaboration | Determination of jet energy calibration and transverse momentum resolution in CMS | JINST 6 (2011) P11002 | CMS-JME-10-011 1107.4277 |
59 | CMS Collaboration | A deep neural network for simultaneous estimation of b jet energy and resolution | Submitted to Comput. Softw. Big Sci | CMS-HIG-18-027 1912.06046 |
60 | 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 |
61 | CMS Collaboration | Pileup mitigation at CMS in 13 TeV data | JINST 15 (2020) P09018 | CMS-JME-18-001 2003.00503 |
62 | N. Kumar and S. P. Martin | LHC search for di-Higgs decays of stoponium and other scalars in events with two photons and two bottom jets | PRD 90 (2014) 055007 | 1404.0996 |
63 | T. Q. Nguyen et al. | Topology classification with deep learning to improve real-time event selection at the LHC | Comput. Softw. Big Sci. 3 (2019) 12 | 1807.00083 |
64 | F. Chollet et al. | |||
65 | M. Abadi et al. | TensorFlow: Large-scale machine learning on heterogeneous systems | 2015 Software available from tensorflow.org. \url http://tensorflow.org/ | |
66 | S. Bolognesi et al. | On the spin and parity of a single-produced resonance at the LHC | PRD 86 (2012) 095031 | 1208.4018 |
67 | 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, p. 785 ACM, New York, NY, USA | |
68 | T. Hastie, R. Tibshirani, and J. Friedman | Springer-Verlag New York, 2nd edition | ||
69 | CMS Collaboration | Performance of quark/gluon discrimination in 8 $ TeV {\mathrm{p}}{\mathrm{p}} $ data | CMS-PAS-JME-13-002 | CMS-PAS-JME-13-002 |
70 | CMS Collaboration | Jet algorithms performance in 13 TeV data | CMS-PAS-JME-16-003 | CMS-PAS-JME-16-003 |
71 | M. J. Oreglia | A study of the reactions $\psi' \to \gamma\gamma \psi$ | PhD thesis, Stanford University, 1980 SLAC Report SLAC-R-236, see Appendix D | |
72 | P. D. Dauncey, M. Kenzie, N. Wardle, and G. J. Davies | Handling uncertainties in background shapes: the discrete profiling method | JINST 10 (2015) P04015 | 1408.6865 |
73 | 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 |
74 | CMS Collaboration | CMS Luminosity Measurements for the 2016 Data Taking Period | CMS-PAS-LUM-17-001 | CMS-PAS-LUM-17-001 |
75 | CMS Collaboration | CMS luminosity measurement for the 2017 data-taking period at $ \sqrt{s} = $ 13 TeV | CMS-PAS-LUM-17-004 | CMS-PAS-LUM-17-004 |
76 | CMS Collaboration | CMS luminosity measurement for the 2018 data-taking period at $ \sqrt{s} = $ 13 TeV | CMS-PAS-LUM-18-002 | CMS-PAS-LUM-18-002 |
77 | A. L. Read | Presentation of search results: the CL$ _s $ technique | JPG 28 (2002) 2693 | |
78 | T. Junk | Confidence level computation for combining searches with small statistics | NIMA 434 (1999) 435 | hep-ex/9902006 |
79 | G. Cowan, K. Cranmer, E. Gross, and O. Vitells | Asymptotic formulae for likelihood-based tests of new physics | EPJC 71 (2011) 1554 | 1007.1727 |
80 | ATLAS and CMS Collaborations | Procedure for the LHC Higgs boson search combination in summer 2011 | CMS-NOTE-2011-005 | |
81 | CMS Collaboration | A measurement of the Higgs boson mass in the diphoton decay channel | PLB 805 (2020) 135425 | CMS-HIG-19-004 2002.06398 |
Compact Muon Solenoid LHC, CERN |