CMS logoCMS event Hgg
Compact Muon Solenoid
LHC, CERN

CMS-PAS-SUS-20-004
Search for higgsinos in channels with two Higgs bosons and missing transverse momentum in proton-proton collisions at $\sqrt{s}=$ 13 TeV
Abstract: Results are presented from a search for physics beyond the standard model in proton-proton collisions at $\sqrt{s}=$ 13 TeV in final states containing two Higgs bosons, each decaying via the process $\text{H}\rightarrow\mathrm{b\overline{b}}$, and large missing transverse momentum. The search uses a data sample accumulated by the CMS experiment at the CERN LHC corresponding to an integrated luminosity of 137 fb$^{-1}$. The search is motivated by models of supersymmetry that predict the production of neutralinos, the neutral partners of the electroweak gauge and Higgs bosons. The observed event yields in the signal regions are found to be consistent with the standard model background expectations based on control regions in the data. The results are interpreted using simplified models of supersymmetry. For the electroweak production of nearly degenerate higgsinos, each of whose decay chains yields a neutralino state $\tilde{\chi}^0_1$ that in turn decays to a massless goldstino and a Higgs boson, $\tilde{\chi}^0_1$ masses in the range 175 to 1025 GeV are excluded at the 95% confidence level. For the strong production of gluino pairs decaying via a slightly lighter $\tilde{\chi}^0_2$ to H and a light $\tilde{\chi}^0_1$, gluino masses below 2330 GeV are excluded.
Figures & Tables Summary References CMS Publications
Figures

png pdf
Figure 1:
Diagrams for (left) the TChiHH-G signal model, $\tilde{\chi}^0_1 \tilde{\chi}^0_1 \to \mathrm{H} \mathrm{H} \tilde{\mathrm{G}} \tilde{\mathrm{G}} $, in which the $\tilde{\chi}^0_1$ NLSPs are produced indirectly through the cascade decays of several combinations of neutralinos and charginos, as described in the text; (center) TChiHH, in which the electroweak production of two neutralinos leads to two Higgs bosons and two neutralinos ($\tilde{\chi}^0_1 $); (right) T5HH, the strong production of a pair of gluinos, each of which decays via a three-body process to quarks and a neutralino, the neutralino subsequently decaying to a Higgs boson and a $\tilde{\chi}^0_1$ LSP. In each diagram, the filled circle represents the sum of processes that can lead to the SUSY particles shown.

png pdf
Figure 1-a:
Diagrams for (left) the TChiHH-G signal model, $\tilde{\chi}^0_1 \tilde{\chi}^0_1 \to \mathrm{H} \mathrm{H} \tilde{\mathrm{G}} \tilde{\mathrm{G}} $, in which the $\tilde{\chi}^0_1$ NLSPs are produced indirectly through the cascade decays of several combinations of neutralinos and charginos, as described in the text; (center) TChiHH, in which the electroweak production of two neutralinos leads to two Higgs bosons and two neutralinos ($\tilde{\chi}^0_1 $); (right) T5HH, the strong production of a pair of gluinos, each of which decays via a three-body process to quarks and a neutralino, the neutralino subsequently decaying to a Higgs boson and a $\tilde{\chi}^0_1$ LSP. In each diagram, the filled circle represents the sum of processes that can lead to the SUSY particles shown.

png pdf
Figure 1-b:
Diagrams for (left) the TChiHH-G signal model, $\tilde{\chi}^0_1 \tilde{\chi}^0_1 \to \mathrm{H} \mathrm{H} \tilde{\mathrm{G}} \tilde{\mathrm{G}} $, in which the $\tilde{\chi}^0_1$ NLSPs are produced indirectly through the cascade decays of several combinations of neutralinos and charginos, as described in the text; (center) TChiHH, in which the electroweak production of two neutralinos leads to two Higgs bosons and two neutralinos ($\tilde{\chi}^0_1 $); (right) T5HH, the strong production of a pair of gluinos, each of which decays via a three-body process to quarks and a neutralino, the neutralino subsequently decaying to a Higgs boson and a $\tilde{\chi}^0_1$ LSP. In each diagram, the filled circle represents the sum of processes that can lead to the SUSY particles shown.

png pdf
Figure 1-c:
Diagrams for (left) the TChiHH-G signal model, $\tilde{\chi}^0_1 \tilde{\chi}^0_1 \to \mathrm{H} \mathrm{H} \tilde{\mathrm{G}} \tilde{\mathrm{G}} $, in which the $\tilde{\chi}^0_1$ NLSPs are produced indirectly through the cascade decays of several combinations of neutralinos and charginos, as described in the text; (center) TChiHH, in which the electroweak production of two neutralinos leads to two Higgs bosons and two neutralinos ($\tilde{\chi}^0_1 $); (right) T5HH, the strong production of a pair of gluinos, each of which decays via a three-body process to quarks and a neutralino, the neutralino subsequently decaying to a Higgs boson and a $\tilde{\chi}^0_1$ LSP. In each diagram, the filled circle represents the sum of processes that can lead to the SUSY particles shown.

png pdf
Figure 2:
Distributions in the primary variables of the resolved signature for events satisfying the baseline requirements plus $ {N_{\text {b}}} =$ 3 or 4, except for those on the variable plotted: (upper left) ${\Delta m_{\mathrm{b} \mathrm{b}}}$, (upper right) ${\Delta R_{\text {max}}}$, (middle left) ${< m_{\mathrm{b} \mathrm{b}}>}$, (middle right) ${N_{\text {b}}}$, and (lower) ${{p_{\mathrm {T}}} ^\text {miss}}$. The rightmost bin includes overflow entries. The data are shown as black markers with error bars, simulated SM backgrounds by the stacked histograms (scaled by factors of 0.93-0.97 to match the total data yield in each plot), and simulations of $\text {TChiHH-G}(m(\tilde{\chi}^0_1),m(\tilde{\mathrm{G}}))$ signals by the lines. Gray shading represents the statistical uncertainties of the simulation. Vertical dashed lines indicate the boundaries for SR selection (dashed), or for the SR binning discussed in Section 7 (dotted). The lower panel shows the ratio of data to (scaled) simulation.

png pdf
Figure 2-a:
Distributions in the primary variables of the resolved signature for events satisfying the baseline requirements plus $ {N_{\text {b}}} =$ 3 or 4, except for those on the variable plotted: (upper left) ${\Delta m_{\mathrm{b} \mathrm{b}}}$, (upper right) ${\Delta R_{\text {max}}}$, (middle left) ${< m_{\mathrm{b} \mathrm{b}}>}$, (middle right) ${N_{\text {b}}}$, and (lower) ${{p_{\mathrm {T}}} ^\text {miss}}$. The rightmost bin includes overflow entries. The data are shown as black markers with error bars, simulated SM backgrounds by the stacked histograms (scaled by factors of 0.93-0.97 to match the total data yield in each plot), and simulations of $\text {TChiHH-G}(m(\tilde{\chi}^0_1),m(\tilde{\mathrm{G}}))$ signals by the lines. Gray shading represents the statistical uncertainties of the simulation. Vertical dashed lines indicate the boundaries for SR selection (dashed), or for the SR binning discussed in Section 7 (dotted). The lower panel shows the ratio of data to (scaled) simulation.

png pdf
Figure 2-b:
Distributions in the primary variables of the resolved signature for events satisfying the baseline requirements plus $ {N_{\text {b}}} =$ 3 or 4, except for those on the variable plotted: (upper left) ${\Delta m_{\mathrm{b} \mathrm{b}}}$, (upper right) ${\Delta R_{\text {max}}}$, (middle left) ${< m_{\mathrm{b} \mathrm{b}}>}$, (middle right) ${N_{\text {b}}}$, and (lower) ${{p_{\mathrm {T}}} ^\text {miss}}$. The rightmost bin includes overflow entries. The data are shown as black markers with error bars, simulated SM backgrounds by the stacked histograms (scaled by factors of 0.93-0.97 to match the total data yield in each plot), and simulations of $\text {TChiHH-G}(m(\tilde{\chi}^0_1),m(\tilde{\mathrm{G}}))$ signals by the lines. Gray shading represents the statistical uncertainties of the simulation. Vertical dashed lines indicate the boundaries for SR selection (dashed), or for the SR binning discussed in Section 7 (dotted). The lower panel shows the ratio of data to (scaled) simulation.

png pdf
Figure 2-c:
Distributions in the primary variables of the resolved signature for events satisfying the baseline requirements plus $ {N_{\text {b}}} =$ 3 or 4, except for those on the variable plotted: (upper left) ${\Delta m_{\mathrm{b} \mathrm{b}}}$, (upper right) ${\Delta R_{\text {max}}}$, (middle left) ${< m_{\mathrm{b} \mathrm{b}}>}$, (middle right) ${N_{\text {b}}}$, and (lower) ${{p_{\mathrm {T}}} ^\text {miss}}$. The rightmost bin includes overflow entries. The data are shown as black markers with error bars, simulated SM backgrounds by the stacked histograms (scaled by factors of 0.93-0.97 to match the total data yield in each plot), and simulations of $\text {TChiHH-G}(m(\tilde{\chi}^0_1),m(\tilde{\mathrm{G}}))$ signals by the lines. Gray shading represents the statistical uncertainties of the simulation. Vertical dashed lines indicate the boundaries for SR selection (dashed), or for the SR binning discussed in Section 7 (dotted). The lower panel shows the ratio of data to (scaled) simulation.

png pdf
Figure 2-d:
Distributions in the primary variables of the resolved signature for events satisfying the baseline requirements plus $ {N_{\text {b}}} =$ 3 or 4, except for those on the variable plotted: (upper left) ${\Delta m_{\mathrm{b} \mathrm{b}}}$, (upper right) ${\Delta R_{\text {max}}}$, (middle left) ${< m_{\mathrm{b} \mathrm{b}}>}$, (middle right) ${N_{\text {b}}}$, and (lower) ${{p_{\mathrm {T}}} ^\text {miss}}$. The rightmost bin includes overflow entries. The data are shown as black markers with error bars, simulated SM backgrounds by the stacked histograms (scaled by factors of 0.93-0.97 to match the total data yield in each plot), and simulations of $\text {TChiHH-G}(m(\tilde{\chi}^0_1),m(\tilde{\mathrm{G}}))$ signals by the lines. Gray shading represents the statistical uncertainties of the simulation. Vertical dashed lines indicate the boundaries for SR selection (dashed), or for the SR binning discussed in Section 7 (dotted). The lower panel shows the ratio of data to (scaled) simulation.

png pdf
Figure 2-e:
Distributions in the primary variables of the resolved signature for events satisfying the baseline requirements plus $ {N_{\text {b}}} =$ 3 or 4, except for those on the variable plotted: (upper left) ${\Delta m_{\mathrm{b} \mathrm{b}}}$, (upper right) ${\Delta R_{\text {max}}}$, (middle left) ${< m_{\mathrm{b} \mathrm{b}}>}$, (middle right) ${N_{\text {b}}}$, and (lower) ${{p_{\mathrm {T}}} ^\text {miss}}$. The rightmost bin includes overflow entries. The data are shown as black markers with error bars, simulated SM backgrounds by the stacked histograms (scaled by factors of 0.93-0.97 to match the total data yield in each plot), and simulations of $\text {TChiHH-G}(m(\tilde{\chi}^0_1),m(\tilde{\mathrm{G}}))$ signals by the lines. Gray shading represents the statistical uncertainties of the simulation. Vertical dashed lines indicate the boundaries for SR selection (dashed), or for the SR binning discussed in Section 7 (dotted). The lower panel shows the ratio of data to (scaled) simulation.

png pdf
Figure 3:
Distributions in the most pertinent variables after the baseline requirements for the boosted signature, except for those on the variable shown: (upper left) ${{p_{\mathrm {T}}} ^\text {miss}}$, (upper right) jet ${p_{\mathrm {T}}}$, (lower left) ${D_{ {\mathrm{b} \mathrm{b}}}}$, and (lower right) ${m_{\text {J}}}$. Except for the ${{p_{\mathrm {T}}} ^\text {miss}}$ distribution, each plot contains two entries per event, for each of the two ${p_{\mathrm {T}}}$-leading wide-cone jets. The data are shown by the black markers with error bars, and SM backgrounds from simulation (scaled by 86% to match the data integral) by the filled histograms. Line histograms show simulations of the signals $\text {TChiHH-G}(m(\tilde{\chi}^0_1),m(\tilde{\mathrm{G}}))$ or $\text {T5HH}(m({\mathrm{\tilde{g}}}),m(\tilde{\chi}^0_1))$. In the lower-left plot, the vertical dashed line denotes the ${D_{ {\mathrm{b} \mathrm{b}}}} $ threshold for classification of the jet as a ${\mathrm{H} \to b {}\mathrm{\bar{b}}}$ candidate; in the lower-right plot, the vertical dashed lines denote the boundaries of the H mass window. The lower panel in each plot shows the ratio of observed to (scaled) simulated yields.

png pdf
Figure 3-a:
Distributions in the most pertinent variables after the baseline requirements for the boosted signature, except for those on the variable shown: (upper left) ${{p_{\mathrm {T}}} ^\text {miss}}$, (upper right) jet ${p_{\mathrm {T}}}$, (lower left) ${D_{ {\mathrm{b} \mathrm{b}}}}$, and (lower right) ${m_{\text {J}}}$. Except for the ${{p_{\mathrm {T}}} ^\text {miss}}$ distribution, each plot contains two entries per event, for each of the two ${p_{\mathrm {T}}}$-leading wide-cone jets. The data are shown by the black markers with error bars, and SM backgrounds from simulation (scaled by 86% to match the data integral) by the filled histograms. Line histograms show simulations of the signals $\text {TChiHH-G}(m(\tilde{\chi}^0_1),m(\tilde{\mathrm{G}}))$ or $\text {T5HH}(m({\mathrm{\tilde{g}}}),m(\tilde{\chi}^0_1))$. In the lower-left plot, the vertical dashed line denotes the ${D_{ {\mathrm{b} \mathrm{b}}}} $ threshold for classification of the jet as a ${\mathrm{H} \to b {}\mathrm{\bar{b}}}$ candidate; in the lower-right plot, the vertical dashed lines denote the boundaries of the H mass window. The lower panel in each plot shows the ratio of observed to (scaled) simulated yields.

png pdf
Figure 3-b:
Distributions in the most pertinent variables after the baseline requirements for the boosted signature, except for those on the variable shown: (upper left) ${{p_{\mathrm {T}}} ^\text {miss}}$, (upper right) jet ${p_{\mathrm {T}}}$, (lower left) ${D_{ {\mathrm{b} \mathrm{b}}}}$, and (lower right) ${m_{\text {J}}}$. Except for the ${{p_{\mathrm {T}}} ^\text {miss}}$ distribution, each plot contains two entries per event, for each of the two ${p_{\mathrm {T}}}$-leading wide-cone jets. The data are shown by the black markers with error bars, and SM backgrounds from simulation (scaled by 86% to match the data integral) by the filled histograms. Line histograms show simulations of the signals $\text {TChiHH-G}(m(\tilde{\chi}^0_1),m(\tilde{\mathrm{G}}))$ or $\text {T5HH}(m({\mathrm{\tilde{g}}}),m(\tilde{\chi}^0_1))$. In the lower-left plot, the vertical dashed line denotes the ${D_{ {\mathrm{b} \mathrm{b}}}} $ threshold for classification of the jet as a ${\mathrm{H} \to b {}\mathrm{\bar{b}}}$ candidate; in the lower-right plot, the vertical dashed lines denote the boundaries of the H mass window. The lower panel in each plot shows the ratio of observed to (scaled) simulated yields.

png pdf
Figure 3-c:
Distributions in the most pertinent variables after the baseline requirements for the boosted signature, except for those on the variable shown: (upper left) ${{p_{\mathrm {T}}} ^\text {miss}}$, (upper right) jet ${p_{\mathrm {T}}}$, (lower left) ${D_{ {\mathrm{b} \mathrm{b}}}}$, and (lower right) ${m_{\text {J}}}$. Except for the ${{p_{\mathrm {T}}} ^\text {miss}}$ distribution, each plot contains two entries per event, for each of the two ${p_{\mathrm {T}}}$-leading wide-cone jets. The data are shown by the black markers with error bars, and SM backgrounds from simulation (scaled by 86% to match the data integral) by the filled histograms. Line histograms show simulations of the signals $\text {TChiHH-G}(m(\tilde{\chi}^0_1),m(\tilde{\mathrm{G}}))$ or $\text {T5HH}(m({\mathrm{\tilde{g}}}),m(\tilde{\chi}^0_1))$. In the lower-left plot, the vertical dashed line denotes the ${D_{ {\mathrm{b} \mathrm{b}}}} $ threshold for classification of the jet as a ${\mathrm{H} \to b {}\mathrm{\bar{b}}}$ candidate; in the lower-right plot, the vertical dashed lines denote the boundaries of the H mass window. The lower panel in each plot shows the ratio of observed to (scaled) simulated yields.

png pdf
Figure 3-d:
Distributions in the most pertinent variables after the baseline requirements for the boosted signature, except for those on the variable shown: (upper left) ${{p_{\mathrm {T}}} ^\text {miss}}$, (upper right) jet ${p_{\mathrm {T}}}$, (lower left) ${D_{ {\mathrm{b} \mathrm{b}}}}$, and (lower right) ${m_{\text {J}}}$. Except for the ${{p_{\mathrm {T}}} ^\text {miss}}$ distribution, each plot contains two entries per event, for each of the two ${p_{\mathrm {T}}}$-leading wide-cone jets. The data are shown by the black markers with error bars, and SM backgrounds from simulation (scaled by 86% to match the data integral) by the filled histograms. Line histograms show simulations of the signals $\text {TChiHH-G}(m(\tilde{\chi}^0_1),m(\tilde{\mathrm{G}}))$ or $\text {T5HH}(m({\mathrm{\tilde{g}}}),m(\tilde{\chi}^0_1))$. In the lower-left plot, the vertical dashed line denotes the ${D_{ {\mathrm{b} \mathrm{b}}}} $ threshold for classification of the jet as a ${\mathrm{H} \to b {}\mathrm{\bar{b}}}$ candidate; in the lower-right plot, the vertical dashed lines denote the boundaries of the H mass window. The lower panel in each plot shows the ratio of observed to (scaled) simulated yields.

png pdf
Figure 4:
Configuration of the search and control regions for the (left) resolved and (right) boosted signatures. The patterns shown are repeated in each of several bins in kinematic or topological variables for improved sensitivity, as discussed in the text.

png pdf
Figure 5:
The double ratio $\kappa = ({N_{\text {SR}}} / {N_{\text {SB}}})/({N_{\text {CSR}}} / {N_{\text {CSB}}})$ from the SM simulation for the 3b and 4b search samples for each $({{p_{\mathrm {T}}} ^\text {miss}}, {\Delta R_{\text {max}}})$ bin of the resolved signature. The value $\Delta $ gives the deviation of ${\kappa}$ from unity, and $\sigma $ its relative statistical uncertainty.

png pdf
Figure 6:
Distributions in (left) the ratio from the SM simulation of the yield in the ${m_{\text {J}}}$ signal region to the yield in the sideband for each of the 0H, 1H, and 2H analysis regions of the boosted signature, and (right) $ {{p_{\mathrm {T}}} ^\text {miss}} $ for the 2H and 1H SRs and the 0H+b control region from SM simulation (MC, solid points), and for the control region in data (open black points). In the right-hand plot the upper panel shows the distributions normalized to unit area. The blue points, and in one bin the open point, are hidden under the solid black points. The lower panel shows the (unnormalized) ratios to the 0H+b yields of the SRs 1H (blue) and 2H (red) for the simulation. The error bars show the statistical uncertainties.

png pdf
Figure 6-a:
Distributions in (left) the ratio from the SM simulation of the yield in the ${m_{\text {J}}}$ signal region to the yield in the sideband for each of the 0H, 1H, and 2H analysis regions of the boosted signature, and (right) $ {{p_{\mathrm {T}}} ^\text {miss}} $ for the 2H and 1H SRs and the 0H+b control region from SM simulation (MC, solid points), and for the control region in data (open black points). In the right-hand plot the upper panel shows the distributions normalized to unit area. The blue points, and in one bin the open point, are hidden under the solid black points. The lower panel shows the (unnormalized) ratios to the 0H+b yields of the SRs 1H (blue) and 2H (red) for the simulation. The error bars show the statistical uncertainties.

png pdf
Figure 6-b:
Distributions in (left) the ratio from the SM simulation of the yield in the ${m_{\text {J}}}$ signal region to the yield in the sideband for each of the 0H, 1H, and 2H analysis regions of the boosted signature, and (right) $ {{p_{\mathrm {T}}} ^\text {miss}} $ for the 2H and 1H SRs and the 0H+b control region from SM simulation (MC, solid points), and for the control region in data (open black points). In the right-hand plot the upper panel shows the distributions normalized to unit area. The blue points, and in one bin the open point, are hidden under the solid black points. The lower panel shows the (unnormalized) ratios to the 0H+b yields of the SRs 1H (blue) and 2H (red) for the simulation. The error bars show the statistical uncertainties.

png pdf
Figure 7:
Distributions in ${< m_{\mathrm{b} \mathrm{b}}>}$ for (above) 3b and (below) 4b data, denoted by black markers, with error bars indicating the statistical uncertainties. The overlaid cyan histogram shows the corresponding distribution of the 2b data multiplied by ${\kappa}$ and scaled in area to the 3b or 4b distribution, with statistical uncertainties indicated by the cyan shading. The ratio of these distributions appears in the lower panel. The red, green, and black histograms show simulations of representative signals, denoted $\text {TChiHH-G}(m(\tilde{\chi}^0_1),m(\tilde{\mathrm{G}}))$ in the legends. The left and right plots show the low- and high- ${\Delta R_{\text {max}}}$ plane, respectively.

png pdf
Figure 7-a:
Distributions in ${< m_{\mathrm{b} \mathrm{b}}>}$ for (above) 3b and (below) 4b data, denoted by black markers, with error bars indicating the statistical uncertainties. The overlaid cyan histogram shows the corresponding distribution of the 2b data multiplied by ${\kappa}$ and scaled in area to the 3b or 4b distribution, with statistical uncertainties indicated by the cyan shading. The ratio of these distributions appears in the lower panel. The red, green, and black histograms show simulations of representative signals, denoted $\text {TChiHH-G}(m(\tilde{\chi}^0_1),m(\tilde{\mathrm{G}}))$ in the legends. The left and right plots show the low- and high- ${\Delta R_{\text {max}}}$ plane, respectively.

png pdf
Figure 7-b:
Distributions in ${< m_{\mathrm{b} \mathrm{b}}>}$ for (above) 3b and (below) 4b data, denoted by black markers, with error bars indicating the statistical uncertainties. The overlaid cyan histogram shows the corresponding distribution of the 2b data multiplied by ${\kappa}$ and scaled in area to the 3b or 4b distribution, with statistical uncertainties indicated by the cyan shading. The ratio of these distributions appears in the lower panel. The red, green, and black histograms show simulations of representative signals, denoted $\text {TChiHH-G}(m(\tilde{\chi}^0_1),m(\tilde{\mathrm{G}}))$ in the legends. The left and right plots show the low- and high- ${\Delta R_{\text {max}}}$ plane, respectively.

png pdf
Figure 7-c:
Distributions in ${< m_{\mathrm{b} \mathrm{b}}>}$ for (above) 3b and (below) 4b data, denoted by black markers, with error bars indicating the statistical uncertainties. The overlaid cyan histogram shows the corresponding distribution of the 2b data multiplied by ${\kappa}$ and scaled in area to the 3b or 4b distribution, with statistical uncertainties indicated by the cyan shading. The ratio of these distributions appears in the lower panel. The red, green, and black histograms show simulations of representative signals, denoted $\text {TChiHH-G}(m(\tilde{\chi}^0_1),m(\tilde{\mathrm{G}}))$ in the legends. The left and right plots show the low- and high- ${\Delta R_{\text {max}}}$ plane, respectively.

png pdf
Figure 7-d:
Distributions in ${< m_{\mathrm{b} \mathrm{b}}>}$ for (above) 3b and (below) 4b data, denoted by black markers, with error bars indicating the statistical uncertainties. The overlaid cyan histogram shows the corresponding distribution of the 2b data multiplied by ${\kappa}$ and scaled in area to the 3b or 4b distribution, with statistical uncertainties indicated by the cyan shading. The ratio of these distributions appears in the lower panel. The red, green, and black histograms show simulations of representative signals, denoted $\text {TChiHH-G}(m(\tilde{\chi}^0_1),m(\tilde{\mathrm{G}}))$ in the legends. The left and right plots show the low- and high- ${\Delta R_{\text {max}}}$ plane, respectively.

png pdf
Figure 8:
Distributions in ${< m_{\mathrm{b} \mathrm{b}}>}$ for (above) 3b and (below) 4b data, denoted by black markers, with error bars indicating the statistical uncertainties. The overlaid cyan histogram shows the corresponding distribution of the 2b data multiplied by ${\kappa}$ and scaled in area to the 3b or 4b distribution, with statistical uncertainties indicated by the cyan shading. The ratio of these distributions appears in the lower panel. The red, green, and black histograms show simulations of representative signals, denoted $\text {TChiHH-G}(m(\tilde{\chi}^0_1),m(\tilde{\mathrm{G}}))$ in the legends. The left and right plots represent the $ {{p_{\mathrm {T}}} ^\text {miss}} < $ 200 GeV and $ {{p_{\mathrm {T}}} ^\text {miss}} \geq $ 300 GeV plane, respectively.

png pdf
Figure 8-a:
Distributions in ${< m_{\mathrm{b} \mathrm{b}}>}$ for (above) 3b and (below) 4b data, denoted by black markers, with error bars indicating the statistical uncertainties. The overlaid cyan histogram shows the corresponding distribution of the 2b data multiplied by ${\kappa}$ and scaled in area to the 3b or 4b distribution, with statistical uncertainties indicated by the cyan shading. The ratio of these distributions appears in the lower panel. The red, green, and black histograms show simulations of representative signals, denoted $\text {TChiHH-G}(m(\tilde{\chi}^0_1),m(\tilde{\mathrm{G}}))$ in the legends. The left and right plots represent the $ {{p_{\mathrm {T}}} ^\text {miss}} < $ 200 GeV and $ {{p_{\mathrm {T}}} ^\text {miss}} \geq $ 300 GeV plane, respectively.

png pdf
Figure 8-b:
Distributions in ${< m_{\mathrm{b} \mathrm{b}}>}$ for (above) 3b and (below) 4b data, denoted by black markers, with error bars indicating the statistical uncertainties. The overlaid cyan histogram shows the corresponding distribution of the 2b data multiplied by ${\kappa}$ and scaled in area to the 3b or 4b distribution, with statistical uncertainties indicated by the cyan shading. The ratio of these distributions appears in the lower panel. The red, green, and black histograms show simulations of representative signals, denoted $\text {TChiHH-G}(m(\tilde{\chi}^0_1),m(\tilde{\mathrm{G}}))$ in the legends. The left and right plots represent the $ {{p_{\mathrm {T}}} ^\text {miss}} < $ 200 GeV and $ {{p_{\mathrm {T}}} ^\text {miss}} \geq $ 300 GeV plane, respectively.

png pdf
Figure 8-c:
Distributions in ${< m_{\mathrm{b} \mathrm{b}}>}$ for (above) 3b and (below) 4b data, denoted by black markers, with error bars indicating the statistical uncertainties. The overlaid cyan histogram shows the corresponding distribution of the 2b data multiplied by ${\kappa}$ and scaled in area to the 3b or 4b distribution, with statistical uncertainties indicated by the cyan shading. The ratio of these distributions appears in the lower panel. The red, green, and black histograms show simulations of representative signals, denoted $\text {TChiHH-G}(m(\tilde{\chi}^0_1),m(\tilde{\mathrm{G}}))$ in the legends. The left and right plots represent the $ {{p_{\mathrm {T}}} ^\text {miss}} < $ 200 GeV and $ {{p_{\mathrm {T}}} ^\text {miss}} \geq $ 300 GeV plane, respectively.

png pdf
Figure 8-d:
Distributions in ${< m_{\mathrm{b} \mathrm{b}}>}$ for (above) 3b and (below) 4b data, denoted by black markers, with error bars indicating the statistical uncertainties. The overlaid cyan histogram shows the corresponding distribution of the 2b data multiplied by ${\kappa}$ and scaled in area to the 3b or 4b distribution, with statistical uncertainties indicated by the cyan shading. The ratio of these distributions appears in the lower panel. The red, green, and black histograms show simulations of representative signals, denoted $\text {TChiHH-G}(m(\tilde{\chi}^0_1),m(\tilde{\mathrm{G}}))$ in the legends. The left and right plots represent the $ {{p_{\mathrm {T}}} ^\text {miss}} < $ 200 GeV and $ {{p_{\mathrm {T}}} ^\text {miss}} \geq $ 300 GeV plane, respectively.

png pdf
Figure 9:
Distributions in ${m_{\text {J}}}$ for the boosted signature, integrated in ${{p_{\mathrm {T}}} ^\text {miss}}$. The projections (left) contain two entries per event, with statistical uncertainties in the data and simulation given by the vertical bars and gray shading, respectively. The SM components are scaled to the data integral, by factors of 0.84, 0.92, and 1.13 in the 0H (upper), 1H (middle), and 2H (lower) plot, respectively. The data to simulation ratio appears in each lower panel. Simulated signals $\text {TChiHH-G}(m(\tilde{\chi}^0_1),m(\tilde{\mathrm{G}}))$ and $\text {T5HH}(m({\mathrm{\tilde{g}}}),m(\tilde{\chi}^0_1))$ are also shown. In the correlation plots (right) the continuous color scale represents the SM background, black dots the data, and red dots the expected signal. The dashed lines and box denote the boundaries of the SR.

png pdf
Figure 9-a:
Distributions in ${m_{\text {J}}}$ for the boosted signature, integrated in ${{p_{\mathrm {T}}} ^\text {miss}}$. The projections (left) contain two entries per event, with statistical uncertainties in the data and simulation given by the vertical bars and gray shading, respectively. The SM components are scaled to the data integral, by factors of 0.84, 0.92, and 1.13 in the 0H (upper), 1H (middle), and 2H (lower) plot, respectively. The data to simulation ratio appears in each lower panel. Simulated signals $\text {TChiHH-G}(m(\tilde{\chi}^0_1),m(\tilde{\mathrm{G}}))$ and $\text {T5HH}(m({\mathrm{\tilde{g}}}),m(\tilde{\chi}^0_1))$ are also shown. In the correlation plots (right) the continuous color scale represents the SM background, black dots the data, and red dots the expected signal. The dashed lines and box denote the boundaries of the SR.

png
Figure 9-b:
Distributions in ${m_{\text {J}}}$ for the boosted signature, integrated in ${{p_{\mathrm {T}}} ^\text {miss}}$. The projections (left) contain two entries per event, with statistical uncertainties in the data and simulation given by the vertical bars and gray shading, respectively. The SM components are scaled to the data integral, by factors of 0.84, 0.92, and 1.13 in the 0H (upper), 1H (middle), and 2H (lower) plot, respectively. The data to simulation ratio appears in each lower panel. Simulated signals $\text {TChiHH-G}(m(\tilde{\chi}^0_1),m(\tilde{\mathrm{G}}))$ and $\text {T5HH}(m({\mathrm{\tilde{g}}}),m(\tilde{\chi}^0_1))$ are also shown. In the correlation plots (right) the continuous color scale represents the SM background, black dots the data, and red dots the expected signal. The dashed lines and box denote the boundaries of the SR.

png pdf
Figure 9-c:
Distributions in ${m_{\text {J}}}$ for the boosted signature, integrated in ${{p_{\mathrm {T}}} ^\text {miss}}$. The projections (left) contain two entries per event, with statistical uncertainties in the data and simulation given by the vertical bars and gray shading, respectively. The SM components are scaled to the data integral, by factors of 0.84, 0.92, and 1.13 in the 0H (upper), 1H (middle), and 2H (lower) plot, respectively. The data to simulation ratio appears in each lower panel. Simulated signals $\text {TChiHH-G}(m(\tilde{\chi}^0_1),m(\tilde{\mathrm{G}}))$ and $\text {T5HH}(m({\mathrm{\tilde{g}}}),m(\tilde{\chi}^0_1))$ are also shown. In the correlation plots (right) the continuous color scale represents the SM background, black dots the data, and red dots the expected signal. The dashed lines and box denote the boundaries of the SR.

png
Figure 9-d:
Distributions in ${m_{\text {J}}}$ for the boosted signature, integrated in ${{p_{\mathrm {T}}} ^\text {miss}}$. The projections (left) contain two entries per event, with statistical uncertainties in the data and simulation given by the vertical bars and gray shading, respectively. The SM components are scaled to the data integral, by factors of 0.84, 0.92, and 1.13 in the 0H (upper), 1H (middle), and 2H (lower) plot, respectively. The data to simulation ratio appears in each lower panel. Simulated signals $\text {TChiHH-G}(m(\tilde{\chi}^0_1),m(\tilde{\mathrm{G}}))$ and $\text {T5HH}(m({\mathrm{\tilde{g}}}),m(\tilde{\chi}^0_1))$ are also shown. In the correlation plots (right) the continuous color scale represents the SM background, black dots the data, and red dots the expected signal. The dashed lines and box denote the boundaries of the SR.

png pdf
Figure 9-e:
Distributions in ${m_{\text {J}}}$ for the boosted signature, integrated in ${{p_{\mathrm {T}}} ^\text {miss}}$. The projections (left) contain two entries per event, with statistical uncertainties in the data and simulation given by the vertical bars and gray shading, respectively. The SM components are scaled to the data integral, by factors of 0.84, 0.92, and 1.13 in the 0H (upper), 1H (middle), and 2H (lower) plot, respectively. The data to simulation ratio appears in each lower panel. Simulated signals $\text {TChiHH-G}(m(\tilde{\chi}^0_1),m(\tilde{\mathrm{G}}))$ and $\text {T5HH}(m({\mathrm{\tilde{g}}}),m(\tilde{\chi}^0_1))$ are also shown. In the correlation plots (right) the continuous color scale represents the SM background, black dots the data, and red dots the expected signal. The dashed lines and box denote the boundaries of the SR.

png
Figure 9-f:
Distributions in ${m_{\text {J}}}$ for the boosted signature, integrated in ${{p_{\mathrm {T}}} ^\text {miss}}$. The projections (left) contain two entries per event, with statistical uncertainties in the data and simulation given by the vertical bars and gray shading, respectively. The SM components are scaled to the data integral, by factors of 0.84, 0.92, and 1.13 in the 0H (upper), 1H (middle), and 2H (lower) plot, respectively. The data to simulation ratio appears in each lower panel. Simulated signals $\text {TChiHH-G}(m(\tilde{\chi}^0_1),m(\tilde{\mathrm{G}}))$ and $\text {T5HH}(m({\mathrm{\tilde{g}}}),m(\tilde{\chi}^0_1))$ are also shown. In the correlation plots (right) the continuous color scale represents the SM background, black dots the data, and red dots the expected signal. The dashed lines and box denote the boundaries of the SR.

png pdf
Figure 10:
Observed and predicted yields in the search regions identified by the legend text. The points with error bars represent the observed yields, the magenta outline bands the background yields derived from the data sidebands with their total uncertainty, and blue shaded bands the values determined by the background-only fit.

png pdf
Figure 11:
Observed and expected upper limits at 95% CL on the cross section for the gauge-mediated symmetry breaking signal model $\tilde{\chi}^0_1 \tilde{\chi}^0_1 \to \mathrm{H} \mathrm{H} \tilde{\mathrm{G}} \tilde{\mathrm{G}} $. The thin dashed black line with green and yellow bands shows the expected limit with its 1- and 2-standard deviation uncertainties, while the solid black line shows the observed limit. The theoretical cross section is indicated by the dashed red line under the assumption that the decay chains leading to the $\tilde{\chi}^0_1 \tilde{\chi}^0_1 $ intermediate state include a degenerate set of all charginos and neutralinos, and by the dotted magenta lines under the assumption that only the combination $\tilde{\chi}^0_1$ $\tilde{\chi}^{0}_2$ is accessible.

png pdf
Figure 12:
Limits at 95% CL on the cross section for the TChiHH signal model in which production of the intermediate state $\tilde{\chi}^{0}_2$ $\tilde{\chi}^{0}_2$ (assumed mass degenerate) is followed by the decay of each to $\tilde{\chi}^0_1$ $\mathrm{H}$. The color scale gives the cross section limit as a function of $(m(\tilde{\chi}^{0}_2),m(\tilde{\chi}^0_1))$. The red solid and dashed contours show the expected limit with its 1-standard deviation uncertainty. At its central value the observed cross section surface does not exclude any part of the plane for this model; only the $+1$-standard deviation (dashed black) surface intersects the theoretical cross section, near $m(\tilde{\chi}^{0}_2) = $ 300 GeV for the smallest values of $m(\tilde{\chi}^0_1)$.

png pdf
Figure 13:
Observed and expected upper limits at 95% CL on the cross section for the simplified model T5HH, the strong production of a pair of gluinos each of which decays via a three-body process to quarks and a neutralino, the neutralino subsequently decaying to a Higgs boson and a $\tilde{\chi}^0_1$ LSP. The thin dashed black line with green and yellow bands shows the expected limit with its 1- and 2-standard deviation uncertainties; the solid black line shows the observed limit, and the dashed red line the theoretical cross section.
Tables

png pdf
Table 1:
Summary of background uncertainties for the resolved signature. The sources are statistical uncertainties in the determination of ${\kappa}$ and the uncertainties $\Delta (\text {data, MC})$ derived from the comparison of simulation with data in the single-lepton, dilepton, and low-$ {\Delta \phi}$ control samples. The unit of ${{p_{\mathrm {T}}} ^\text {miss}}$ is GeV, and the last column gives the bin-to-bin correlation $\rho $.

png pdf
Table 2:
Summary of systematic uncertainties in the background estimation for the boosted signature alongside the data sideband statistical uncertainties. The values from the ABCD measurement with the ${p_{\mathrm {T}}}$-integrated sample appear in the rows labeled $ {{p_{\mathrm {T}}} ^\text {miss}} > $ 300 GeV. The row labeled $\Delta (\text {data, MC})$ gives the contribution derived from the comparison of simulation with data in the one-lepton control sample, and those labeled "MC corrections'' refer to the corrections to simulation discussed in Section 5.

png pdf
Table 3:
Sources of systematic uncertainties and their typical impact on the signal yields obtained from simulation. The range is reported as the median 68% among all signal regions for every signal mass point considered. Entries reported as 0 correspond to values smaller than 0.5%.

png pdf
Table 4:
For each SR of the resolved signature, the MC closure factor ${\kappa}$, predicted background yield $N_{\text {bkg}}$ (pre-fit), background yield from the background-only fit $N_{\text {bkg}}$ (post-fit), and observed yield $N_{\text {obs}}$. The first and second uncertainties in the ${\kappa}$ factors are statistical and systematic, respectively. The uncertainties in the predictions, extracted from the maximum-likelihood fit described in the text, include both statistical and systematic contributions. The interpretation of results in the bin corresponding to $ {\Delta R_{\text {max}}} < $ 1.1, $ {N_{\text {b}}} =$ 3, 300 $ < {{p_{\mathrm {T}}} ^\text {miss}} < $ 400 GeV is discussed in the text.

png pdf
Table 5:
For each ${N_{{\mathrm{H}}}}$ SR of the boosted signature, the total predicted background yield ${{N_{\text {SR}}}}$, and for each ${{p_{\mathrm {T}}} ^\text {miss}}$ bin the fraction ${f_{\text {bkg}}}$, both with their statistical uncertainties, the predicted background yield $N_{\text {bkg}}$-pre-fit, the yield from the background-only fit $N_{\text {bkg}}$-post-fit, and the observed yield $N_{\text {obs}}$. Both ${{N_{\text {SR}}}}$ and ${f_{\text {bkg}}}$ are determined from data control samples. The $N_{\text {bkg}}$ values are extracted from the maximum-likelihood fit described in the text, with uncertainties that include both statistical and systematic contributions.
Summary
A search is presented for physics beyond the standard model in final states with pairs of Higgs bosons and an imbalance of transverse momentum, produced in proton-proton collisions at $\sqrt{s} = $ 13 TeV. The data, collected by the CMS experiment at the LHC, correspond to an integrated luminosity of 137 fb$^{-1}$. The Higgs bosons are reconstructed via their decay to a pair of b quarks, observed either as distinct b quark jets or as wide-cone jets each containing the pair of b quarks. No significant excess of events beyond those predicted by standard model processes is observed. This finding is used to set limits on the cross sections for the production of SUSY particles, considering both the direct production of neutralinos and their production through intermediate states with gluinos.

For the electroweak production of nearly-degenerate higgsinos each of whose decay cascades yield a $\tilde{\chi}^0_1$ that in turn decays through a Higgs boson to the lightest SUSY particle (LSP), a massless goldstino, $\tilde{\chi}^0_1$ masses in the range 175 to 1025 GeV are excluded at the 95% confidence level. For a model with a mass splitting between the directly produced higgsinos and a bino LSP, the cross section upper limit lies slightly above the model cross section over the entire mass-parameter region where the experiment has expected sensitivity. For the strong production of gluino pairs decaying via a slightly lighter $\tilde{\chi}^0_2$ to a Higgs boson and a light $\tilde{\chi}^0_1$ LSP, gluino masses below 2330 GeV are excluded.
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 ATLAS and CMS Collaborations Combined measurement of the Higgs boson mass in pp collisions at $ \sqrt{s}= $ 7 and 8 TeV with the ATLAS and CMS experiments PRL 114 (2015) 191803 1503.07589
5 P. Ramond Dual theory for free fermions PRD 3 (1971) 2415
6 Y. A. Golfand and E. P. Likhtman Extension of the algebra of Poincar$ \'e $ group generators and violation of P invariance JEPTL 13 (1971)323
7 A. Neveu and J. H. Schwarz Factorizable dual model of pions NPB 31 (1971) 86
8 D. V. Volkov and V. P. Akulov Possible universal neutrino interaction JEPTL 16 (1972)438
9 J. Wess and B. Zumino A Lagrangian model invariant under supergauge transformations PLB 49 (1974) 52
10 J. Wess and B. Zumino Supergauge transformations in four dimensions NPB 70 (1974) 39
11 P. Fayet Supergauge invariant extension of the Higgs mechanism and a model for the electron and its neutrino NPB 90 (1975) 104
12 P. Fayet and S. Ferrara Supersymmetry PR 32 (1977) 249
13 H. P. Nilles Supersymmetry, supergravity and particle physics Phys. Rep. 110 (1984) 1
14 S. P. Martin A supersymmetry primer Adv. Ser. Direct. High Energy Phys. 21 (2010) 1 hep-ph/9709356
15 G. 't Hooft Naturalness, chiral symmetry, and spontaneous chiral symmetry breaking NATO Sci. Ser. B 59 (1980) 135
16 E. Witten Dynamical breaking of supersymmetry NPB 188 (1981) 513
17 M. Dine, W. Fischler, and M. Srednicki Supersymmetric technicolor NPB 189 (1981) 575
18 S. Dimopoulos and S. Raby Supercolor NPB 192 (1981) 353
19 S. Dimopoulos and H. Georgi Softly broken supersymmetry and SU(5) NPB 193 (1981) 150
20 R. K. Kaul and P. Majumdar Cancellation of quadratically divergent mass corrections in globally supersymmetric spontaneously broken gauge theories NPB 199 (1982) 36
21 S. Dimopoulos and G. F. Giudice Naturalness constraints in supersymmetric theories with nonuniversal soft terms PLB 357 (1995) 573 hep-ph/9507282
22 R. Barbieri and D. Pappadopulo S-particles at their naturalness limits JHEP 10 (2009) 061 0906.4546
23 M. Papucci, J. T. Ruderman, and A. Weiler Natural SUSY endures JHEP 09 (2012) 035 1110.6926
24 G. R. Farrar and P. Fayet Phenomenology of the production, decay, and detection of new hadronic states associated with supersymmetry PLB 76 (1978) 575
25 Particle Data Group, P. A. Zyla et al. Review of particle physics Prog. Theor. Exp. Phys. 2020 (2020) 083C01
26 S. Dimopoulos, M. Dine, S. Raby, and S. D. Thomas Experimental signatures of low-energy gauge mediated supersymmetry breaking PRL 76 (1996) 3494 hep-ph/9601367
27 K. T. Matchev and S. D. Thomas Higgs and $ Z $ boson signatures of supersymmetry PRD 62 (2000) 077702 hep-ph/9908482
28 N. Arkani-Hamed et al. MARMOSET: The path from LHC data to the new standard model via on-shell effective theories hep-ph/0703088
29 J. Alwall, M.-P. Le, M. Lisanti, and J. G. Wacker Model-independent jets plus missing energy searches PRD 79 (2009) 015005 0809.3264
30 J. Alwall, P. Schuster, and N. Toro Simplified models for a first characterization of new physics at the LHC PRD 79 (2009) 075020 0810.3921
31 D. Alves et al. Simplified models for LHC new physics searches JPG 39 (2012) 105005 1105.2838
32 J. T. Ruderman and D. Shih General Neutralino NLSPs at the Early LHC JHEP 08 (2012) 159 1103.6083
33 B. Fuks, M. Klasen, D. R. Lamprea, and M. Rothering Gaugino production in proton-proton collisions at a center-of-mass energy of 8 TeV JHEP 10 (2012) 081 1207.2159
34 B. Fuks, M. Klasen, D. R. Lamprea, and M. Rothering Precision predictions for electroweak superpartner production at hadron colliders with $ \sc $ Resummino EPJC 73 (2013) 2480 1304.0790
35 H. Baer et al. Natural SUSY with a bino- or wino-like LSP PRD 91 (2015) 075005 1501.06357
36 Z. Kang, P. Ko, and J. Li New Physics Opportunities in the Boosted Di-Higgs-Boson Plus Missing Transverse Energy Signature PRL 116 (2016) 131801 1504.04128
37 G. D. Kribs, A. Martin, T. S. Roy, and M. Spannowsky Discovering Higgs Bosons of the MSSM using Jet Substructure PRD 82 (2010) 095012 1006.1656
38 S. Gori, P. Schwaller, and C. E. M. Wagner Search for Higgs Bosons in SUSY Cascade Decays and Neutralino Dark Matter PRD 83 (2011) 115022 1103.4138
39 LHC Higgs Cross Section Working Group Handbook of LHC Higgs cross sections: 4. deciphering the nature of the Higgs sector CERN (2016) 1610.07922
40 ATLAS Collaboration Search for supersymmetry in events with photons, bottom quarks, and missing transverse momentum in proton-proton collisions at a centre-of-mass energy of 7 TeV with the ATLAS detector PLB 719 (2013) 261 1211.1167
41 ATLAS Collaboration Search for direct pair production of a chargino and a neutralino decaying to the 125 GeV Higgs boson in $ \sqrt{s} = 8 TeV {pp} $ collisions with the ATLAS detector EPJC 75 (2015) 208 1501.07110
42 ATLAS Collaboration Search for supersymmetry in final states with two same-sign or three leptons and jets using 36 fb$ ^{-1} $ of $ \sqrt{s} = $ 13 TeV pp collision data with the ATLAS detector JHEP 09 (2017) 084 1706.03731
43 ATLAS Collaboration Search for new phenomena with large jet multiplicities and missing transverse momentum using large-radius jets and flavour-tagging at ATLAS in 13 TeV pp collisions JHEP 12 (2017) 034 1708.02794
44 ATLAS Collaboration Search for squarks and gluinos in events with an isolated lepton, jets, and missing transverse momentum at $ \sqrt{s}= $ 13 TeV with the ATLAS detector PRD 96 (2017) 112010 1708.08232
45 ATLAS Collaboration Search for squarks and gluinos in final states with jets and missing transverse momentum using 36 fb$ ^{-1} $ of $ \sqrt{s}= $ 13 TeV pp collision data with the ATLAS detector PRD 97 (2018) 112001 1712.02332
46 ATLAS Collaboration Search for electroweak production of supersymmetric states in scenarios with compressed mass spectra at $ \sqrt{s}= $ 13 TeV with the ATLAS detector PRD 97 (2018) 052010 1712.08119
47 ATLAS Collaboration Search for photonic signatures of gauge-mediated supersymmetry in 13 TeV pp collisions with the ATLAS detector PRD 97 (2018) 092006 1802.03158
48 ATLAS Collaboration Search for electroweak production of supersymmetric particles in final states with two or three leptons at $ \sqrt{s}= $ 13 TeV with the ATLAS detector EPJC 78 (2018) 995 1803.02762
49 ATLAS Collaboration Search for chargino-neutralino production using recursive jigsaw reconstruction in final states with two or three charged leptons in proton-proton collisions at $ \sqrt{s}= $ 13 TeV with the ATLAS detector PRD 98 (2018) 092012 1806.02293
50 ATLAS Collaboration Search for pair production of higgsinos in final states with at least three $ b $-tagged jets in $ \sqrt{s} = $ 13 TeV pp collisions using the ATLAS detector PRD 98 (2018) 092002 1806.04030
51 ATLAS Collaboration Search for chargino and neutralino production in final states with a Higgs boson and missing transverse momentum at $ \sqrt{s} = $ 13 TeV with the ATLAS detector PRD 100 (2019) 012006 1812.09432
52 ATLAS Collaboration Search for electroweak production of charginos and sleptons decaying into final states with two leptons and missing transverse momentum in $ \sqrt{s} = $ 13 TeV pp collisions using the ATLAS detector EPJC 80 (2020) 123 1908.08215
53 ATLAS Collaboration Search for direct production of electroweakinos in final states with one lepton, missing transverse momentum and a Higgs boson decaying into two $ b $-jets in $ pp $ collisions at $ \sqrt{s}= $ 13 TeV with the ATLAS detector EPJC 80 (2020) 691 1909.09226
54 ATLAS Collaboration Searches for electroweak production of supersymmetric particles with compressed mass spectra in $ \sqrt{s} = $ 13 TeV pp collisions with the ATLAS detector PRD 101 (2020) 052005 1911.12606
55 ATLAS Collaboration Search for chargino-neutralino production with mass splittings near the electroweak scale in three-lepton final states in $ \sqrt{s} = $ 13 TeV pp collisions with the ATLAS detector PRD 101 (2020) 072001 1912.08479
56 ATLAS Collaboration Search for direct production of electroweakinos in final states with missing transverse momentum and a Higgs boson decaying into photons in pp collisions at $ \sqrt{s} = $ 13 TeV with the ATLAS detector JHEP 10 (2020) 005 2004.10894
57 ATLAS Collaboration Search for supersymmetry in events with four or more charged leptons in 139 fb$^{-1} $ of $ \sqrt{s} = $ 13 TeV pp collisions with the ATLAS detector 2103.11684
58 CMS Collaboration Searches for electroweak production of charginos, neutralinos, and sleptons decaying to leptons and W, Z, and Higgs bosons in pp collisions at 8 TeV EPJC 74 (2014) 3036 CMS-SUS-13-006
1405.7570
59 CMS Collaboration Search for top squark and higgsino production using diphoton Higgs boson decays PRL 112 (2014) 161802 CMS-SUS-13-014
1312.3310
60 CMS Collaboration Searches for electroweak neutralino and chargino production in channels with Higgs, Z, and W bosons in pp collisions at 8 TeV PRD 90 (2014) 092007 CMS-SUS-14-002
1409.3168
61 CMS Collaboration Search for physics beyond the standard model in events with two leptons of same sign, missing transverse momentum, and jets in proton-proton collisions at $ \sqrt{s} = $ 13 TeV EPJC 77 (2017) 578 CMS-SUS-16-035
1704.07323
62 CMS Collaboration Search for electroweak production of charginos and neutralinos in WH events in proton-proton collisions at $ \sqrt{s}= $ 13 TeV JHEP 11 (2017) 029 CMS-SUS-16-043
1706.09933
63 CMS Collaboration Search for supersymmetry in events with at least one photon, missing transverse momentum, and large transverse event activity in proton-proton collisions at $ \sqrt{s}= $ 13 TeV JHEP 12 (2017) 142 CMS-SUS-16-047
1707.06193
64 CMS Collaboration Search for supersymmetry with Higgs boson to diphoton decays using the razor variables at $ \sqrt{s} = $ 13 TeV PLB 779 (2018) 166 CMS-SUS-16-045
1709.00384
65 CMS Collaboration Search for Higgsino pair production in $ pp $ collisions at $ \sqrt{s} = $ 13 TeV in final states with large missing transverse momentum and two Higgs bosons decaying via $ H \to b \bar b $ PRD 97 (2018) 032007 CMS-SUS-16-044
1709.04896
66 CMS Collaboration Combined search for electroweak production of charginos and neutralinos in proton-proton collisions at $ \sqrt{s} = $ 13 TeV JHEP 03 (2018) 160 CMS-SUS-17-004
1801.03957
67 CMS Collaboration Search for new phenomena in final states with two opposite-charge, same-flavor leptons, jets, and missing transverse momentum in pp collisions at $ \sqrt{s}= $ 13 TeV JHEP 03 (2018) 076 CMS-SUS-16-034
1709.08908
68 CMS Collaboration Search for supersymmetry in events with one lepton and multiple jets exploiting the angular correlation between the lepton and the missing transverse momentum in proton-proton collisions at $ \sqrt{s} = $ 13 TeV PLB 780 (2018) 384 CMS-SUS-16-042
1709.09814
69 CMS Collaboration Search for supersymmetry in events with at least three electrons or muons, jets, and missing transverse momentum in proton-proton collisions at $ \sqrt{s}= $ 13 TeV JHEP 02 (2018) 067 CMS-SUS-16-041
1710.09154
70 CMS Collaboration Search for gauge-mediated supersymmetry in events with at least one photon and missing transverse momentum in pp collisions at $ \sqrt{s} = $ 13 TeV PLB 780 (2018) 118 CMS-SUS-16-046
1711.08008
71 CMS Collaboration Search for Physics Beyond the Standard Model in Events with High-Momentum Higgs Bosons and Missing Transverse Momentum in Proton-Proton Collisions at 13 TeV PRL 120 (2018) 241801 CMS-SUS-17-006
1712.08501
72 CMS Collaboration Search for new physics in events with two soft oppositely charged leptons and missing transverse momentum in proton-proton collisions at $ \sqrt{s}= $ 13 TeV PLB 782 (2018) 440 CMS-SUS-16-048
1801.01846
73 CMS Collaboration Searches for pair production of charginos and top squarks in final states with two oppositely charged leptons in proton-proton collisions at $ \sqrt{s}= $ 13 TeV JHEP 11 (2018) 079 CMS-SUS-17-010
1807.07799
74 CMS Collaboration Search for supersymmetry in events with a photon, a lepton, and missing transverse momentum in proton-proton collisions at $ \sqrt{s} = $ 13 TeV JHEP 01 (2019) 154 CMS-SUS-17-012
1812.04066
75 CMS Collaboration Search for supersymmetry in events with a photon, jets, $ \mathrm {b} $ -jets, and missing transverse momentum in proton-proton collisions at 13 TeV EPJC 79 (2019) 444 CMS-SUS-18-002
1901.06726
76 CMS Collaboration Search for supersymmetry in final states with photons and missing transverse momentum in proton-proton collisions at 13 TeV JHEP 06 (2019) 143 CMS-SUS-17-011
1903.07070
77 CMS Collaboration Combined search for supersymmetry with photons in proton-proton collisions at $ \sqrt{s}= $ 13 TeV PLB 801 (2020) 135183 CMS-SUS-18-005
1907.00857
78 CMS Collaboration Search for supersymmetry using Higgs boson to diphoton decays at $ \sqrt{s} = $ 13 TeV JHEP 11 (2019) 109 CMS-SUS-18-007
1908.08500
79 CMS Collaboration Searches for physics beyond the standard model with the $ M_\mathrm{T2} $ variable in hadronic final states with and without disappearing tracks in proton-proton collisions at $ \sqrt{s}= $ 13 TeV EPJC 80 (2020) 3 CMS-SUS-19-005
1909.03460
80 CMS Collaboration Search for physics beyond the standard model in events with jets and two same-sign or at least three charged leptons in proton-proton collisions at $ \sqrt{s}= $ 13 TeV EPJC 80 (2020) 752 CMS-SUS-19-008
2001.10086
81 CMS Collaboration Search for supersymmetry in proton-proton collisions at $ \sqrt{s} = $ 13 TeV in events with high-momentum Z bosons and missing transverse momentum JHEP 09 (2020) 149 CMS-SUS-19-013
2008.04422
82 CMS Collaboration Search for supersymmetry in final states with two oppositely charged same-flavor leptons and missing transverse momentum in proton-proton collisions at $ \sqrt{s} = $ 13 TeV JHEP 04 (2021) 123 CMS-SUS-20-001
2012.08600
83 CMS Collaboration Search for electroweak production of charginos and neutralinos in proton-proton collisions at $ \sqrt{s} = $ 13 TeV CMS-SUS-19-012
2106.14246
84 CMS Collaboration The CMS experiment at the CERN LHC JINST 3 (2008) S08004 CMS-00-001
85 CMS Collaboration Particle-flow reconstruction and global event description with the CMS detector JINST 12 (2017) P10003 CMS-PRF-14-001
1706.04965
86 M. Cacciari, G. P. Salam, and G. Soyez The anti-$ k_{\rm t} $ jet clustering algorithm JHEP 04 (2008) 063 0802.1189
87 M. Cacciari, G. P. Salam, and G. Soyez FastJet User Manual EPJC 72 (2012) 1896 1111.6097
88 CMS Collaboration Jet performance in pp collisions at $ \sqrt{s}= $ 7 TeV CMS-PAS-JME-10-003
89 CMS Collaboration Jet algorithms performance in 13 TeV data CMS-PAS-JME-16-003 CMS-PAS-JME-16-003
90 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
91 M. Cacciari and G. P. Salam Pileup subtraction using jet areas PLB 659 (2008) 119 0707.1378
92 D. Bertolini, P. Harris, M. Low, and N. Tran Pileup Per Particle Identification JHEP 10 (2014) 059 1407.6013
93 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
94 CMS Collaboration Performance of deep tagging algorithms for boosted double quark jet topology inproton-proton collisions at 13 tev with the phase-0 cms detector CDS
95 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
96 CMS Collaboration Electron and photon reconstruction and identification with the CMS experiment at the CERN LHC CMS-EGM-17-001
2012.06888
97 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
98 K. Rehermann and B. Tweedie Efficient Identification of Boosted Semileptonic Top Quarks at the LHC JHEP 03 (2011) 059 1007.2221
99 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
100 J. Alwall et al. Comparative study of various algorithms for the merging of parton showers and matrix elements in hadronic collisions EPJC 53 (2008) 473 0706.2569
101 R. Frederix and S. Frixione Merging meets matching in MC@NLO JHEP 12 (2012) 061 1209.6215
102 P. Nason A new method for combining NLO QCD with shower Monte Carlo algorithms JHEP 11 (2004) 040 hep-ph/0409146
103 S. Frixione, P. Nason, and C. Oleari Matching NLO QCD computations with parton shower simulations: the POWHEG method JHEP 11 (2007) 070 0709.2092
104 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
105 S. Alioli, P. Nason, C. Oleari, and E. Re NLO single-top production matched with shower in POWHEG: $ s $- and $ t $-channel contributions JHEP 09 (2009) 111 0907.4076
106 E. Re Single-top Wt-channel production matched with parton showers using the POWHEG method EPJC 71 (2011) 1547 1009.2450
107 GEANT4 Collaboration GEANT4--a simulation toolkit NIMA 506 (2003) 250
108 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
109 M. Beneke, P. Falgari, S. Klein, and C. Schwinn Hadronic top-quark pair production with NNLL threshold resummation NPB 855 (2012) 695 1109.1536
110 M. Cacciari et al. Top-pair production at hadron colliders with next-to-next-to-leading logarithmic soft-gluon resummation PLB 710 (2012) 612 1111.5869
111 P. Barnreuther, M. Czakon, and A. Mitov Percent level precision physics at the Tevatron: First genuine NNLO QCD corrections to $ \mathrm{q\bar{q}}\to\mathrm{t\bar{t}} +X $ PRL 109 (2012) 132001 1204.5201
112 M. Czakon and A. Mitov NNLO corrections to top-pair production at hadron colliders: the all-fermionic scattering channels JHEP 12 (2012) 054 1207.0236
113 M. Czakon and A. Mitov NNLO corrections to top pair production at hadron colliders: the quark-gluon reaction JHEP 01 (2013) 080 1210.6832
114 M. Czakon, P. Fiedler, and A. Mitov Total top-quark pair-production cross section at hadron colliders through $ O({\alpha_S}^4) $ PRL 110 (2013) 252004 1303.6254
115 R. Gavin, Y. Li, F. Petriello, and S. Quackenbush W physics at the LHC with FEWZ 2.1 CPC 184 (2013) 208 1201.5896
116 R. Gavin, Y. Li, F. Petriello, and S. Quackenbush FEWZ 2.0: A code for hadronic Z production at next-to-next-to-leading order CPC 182 (2011) 2388 1011.3540
117 W. Beenakker, R. Hopker, M. Spira, and P. M. Zerwas Squark and gluino production at hadron colliders NPB 492 (1997) 51 hep-ph/9610490
118 A. Kulesza and L. Motyka Threshold resummation for squark-antisquark and gluino-pair production at the LHC PRL 102 (2009) 111802 0807.2405
119 A. Kulesza and L. Motyka Soft gluon resummation for the production of gluino-gluino and squark-antisquark pairs at the LHC PRD 80 (2009) 095004 0905.4749
120 W. Beenakker et al. Soft-gluon resummation for squark and gluino hadroproduction JHEP 12 (2009) 041 0909.4418
121 W. Beenakker et al. Squark and gluino hadroproduction Int. J. Mod. Phys. A 26 (2011) 2637 1105.1110
122 W. Beenakker et al. NNLL-fast: predictions for coloured supersymmetric particle production at the LHC with threshold and Coulomb resummation JHEP 12 (2016) 133 1607.07741
123 W. Beenakker et al. NNLL resummation for squark-antisquark pair production at the LHC JHEP 01 (2012) 076 1110.2446
124 W. Beenakker et al. Towards NNLL resummation: hard matching coefficients for squark and gluino hadroproduction JHEP 10 (2013) 120 1304.6354
125 W. Beenakker et al. NNLL resummation for squark and gluino production at the LHC JHEP 12 (2014) 023 1404.3134
126 W. Beenakker et al. Stop production at hadron colliders NPB 515 (1998) 3 hep-ph/9710451
127 W. Beenakker et al. Supersymmetric top and bottom squark production at hadron colliders JHEP 08 (2010) 098 1006.4771
128 W. Beenakker et al. NNLL resummation for stop pair-production at the LHC JHEP 05 (2016) 153 1601.02954
129 S. Abdullin et al. The fast simulation of the CMS detector at LHC J. Phys. Conf. Ser. 331 (2011) 032049
130 A. Giammanco The fast simulation of the CMS experiment J. Phys. Conf. Ser. 513 (2014) 022012
131 T. Sjostrand et al. An Introduction to PYTHIA 8.2 CPC 191 (2015) 159 1410.3012
132 CMS Collaboration Event generator tunes obtained from underlying event and multiparton scattering measurements EPJC 76 (2016) 155 CMS-GEN-14-001
1512.00815
133 CMS Collaboration Extraction and validation of a new set of CMS PYTHIA8 tunes from underlying-event measurements Submitted to EPJC CMS-GEN-17-001
1903.12179
134 NNPDF Collaboration Parton distributions with QED corrections NPB 877 (2013) 290 1308.0598
135 NNPDF Collaboration Parton distributions from high-precision collider data EPJC 77 (2017) 663 1706.00428
136 CMS Collaboration Search for Top-Squark Pair Production in the Single-Lepton Final State in pp Collisions at $ \sqrt{s} = $ 8 TeV EPJC 73 (2013) 2677 CMS-SUS-13-011
1308.1586
137 CMS Collaboration The CMS trigger system JINST 12 (2017) P01020 CMS-TRG-12-001
1609.02366
138 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
139 A. J. Larkoski, S. Marzani, G. Soyez, and J. Thaler Soft Drop JHEP 05 (2014) 146 1402.2657
140 M. Dasgupta, A. Fregoso, S. Marzani, and G. P. Salam Towards an understanding of jet substructure JHEP 09 (2013) 029 1307.0007
141 A. Kalogeropoulos and J. Alwall The SysCalc code: A tool to derive theoretical systematic uncertainties 1801.08401
142 S. Catani, D. de Florian, M. Grazzini, and P. Nason Soft gluon resummation for Higgs boson production at hadron colliders JHEP 07 (2003) 028 hep-ph/0306211
143 M. Cacciari et al. The $ \mathrm{t\bar{t}} $ cross-section at 1.8 TeV and 1.96 TeV: a study of the systematics due to parton densities and scale dependence JHEP 04 (2004) 068 hep-ph/0303085
144 CMS Collaboration Measurement of the inelastic proton-proton cross section at $ \sqrt{s}= $ 13 TeV JHEP 07 (2018) 161 CMS-FSQ-15-005
1802.02613
145 CMS Collaboration Precision luminosity measurement in proton-proton collisions at $ \sqrt{s} = $ 13 TeV in 2015 and 2016 at CMS Submitted to EPJC CMS-LUM-17-003
2104.01927
146 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
147 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
148 G. Cowan, K. Cranmer, E. Gross, and O. Vitells Asymptotic formulae for likelihood-based tests of new physics EPJC 71 (2011) 1554 1007.1727
149 T. Junk Confidence level computation for combining searches with small statistics NIMA 434 (1999) 435 hep-ex/9902006
150 A. L. Read Presentation of search results: the CLs technique JPG 28 (2002) 2693
Compact Muon Solenoid
LHC, CERN