Compact Muon Solenoid
LHC, CERN
Visit us: CMS Public Website, CMS Physics ;
Contact us: CMS Publications Committee
| CMS-PAS-HIG-22-007 | ||
| Search for exotic Higgs boson decays to a pair of pseudoscalars in the $ \mu\mu\mathrm{bb} $ and $ \tau\tau\mathrm{bb} $ final states in proton-proton collisions with the CMS experiment | ||
| CMS Collaboration | ||
| 23 March 2023 | ||
| Abstract: A search for the exotic decays of a standard model like Higgs boson (H) with a 125 GeV mass to a pair of light pseudoscalars $ \mathrm{a}_1 $ is performed in final states where one pseudoscalar decays to two b quarks and the other to two $ \tau $ leptons or muons. A data sample of proton-proton collisions at $ \sqrt{s}= $ 13 TeV corresponding to an integrated luminosity of 138 fb$ ^{-1} $ recorded with the CMS detector is exploited. No statistically significant excess is observed over the standard model backgrounds. Upper limits are set, at 95% confidence level (CL), on the Higgs boson branching fraction to $ \ell\ell\mathrm{b}\mathrm{b} $ via a pair of $ \mathrm{a}_1 $, where $ \ell $ stands for muons and tau leptons. The limits depend on the pseudoscalar mass $ m_{\mathrm{a}_1} $. The observed limits are in the range (0.17-3.3) $\times$ 10$^{-4} $ and (1.7-7.6) $\times$ 10$^{-2} $ in the $ \mu\mu\mathrm{b}\mathrm{b} $ and $ \tau\tau\mathrm{b}\mathrm{b} $ final states, respectively. The two final states are combined to obtain exclusion limits of the branching fraction $ \mathrm{B}(\mathrm{H}\to\mathrm{a}_1\mathrm{a}_1\to\ell\ell\mathrm{b}\mathrm{b}) $ at 95% CL for a broad class of models of a two Higgs doublet extended with a scalar singlet (2HDM+S). Upper bounds on the Higgs boson branching fraction $ \mathrm{B}(\mathrm{H}\to\mathrm{a}_1\mathrm{a}_1) $ are also extracted from the combination. $ \mathrm{B}(\mathrm{H}\to\mathrm{a}_1\mathrm{a}_1) $ values above 0.23 are excluded at 95% CL for most Type-II 2HDM+S models for $ m_{\mathrm{a}_1} $ values between 15 and 60 GeV. | ||
|
Links:
CDS record (PDF) ;
CADI line (restricted) ;
These preliminary results are superseded in this paper, Submitted to EPJC. The superseded preliminary plots can be found here. |
||
| Figures | |
png pdf |
Figure 1:
The distribution of (top) leading and subleading muon $ p_{\mathrm{T}} $ and (bottom) leading and subleading jet $ p_{\mathrm{T}} $ in selected events. The uncertainty band in the lower panel represents the limited size of simulated samples together with a 30% uncertainty on the low mass Drell-Yan cross section. Simulated samples are normalized to 138 fb$^{-1}$ with the corresponding theoretical cross sections. The benchmark criteria explained in the text are used for signal normalization. |
png pdf |
Figure 1-a:
The distribution of (top) leading and subleading muon $ p_{\mathrm{T}} $ and (bottom) leading and subleading jet $ p_{\mathrm{T}} $ in selected events. The uncertainty band in the lower panel represents the limited size of simulated samples together with a 30% uncertainty on the low mass Drell-Yan cross section. Simulated samples are normalized to 138 fb$^{-1}$ with the corresponding theoretical cross sections. The benchmark criteria explained in the text are used for signal normalization. |
png pdf |
Figure 1-b:
The distribution of (top) leading and subleading muon $ p_{\mathrm{T}} $ and (bottom) leading and subleading jet $ p_{\mathrm{T}} $ in selected events. The uncertainty band in the lower panel represents the limited size of simulated samples together with a 30% uncertainty on the low mass Drell-Yan cross section. Simulated samples are normalized to 138 fb$^{-1}$ with the corresponding theoretical cross sections. The benchmark criteria explained in the text are used for signal normalization. |
png pdf |
Figure 1-c:
The distribution of (top) leading and subleading muon $ p_{\mathrm{T}} $ and (bottom) leading and subleading jet $ p_{\mathrm{T}} $ in selected events. The uncertainty band in the lower panel represents the limited size of simulated samples together with a 30% uncertainty on the low mass Drell-Yan cross section. Simulated samples are normalized to 138 fb$^{-1}$ with the corresponding theoretical cross sections. The benchmark criteria explained in the text are used for signal normalization. |
png pdf |
Figure 1-d:
The distribution of (top) leading and subleading muon $ p_{\mathrm{T}} $ and (bottom) leading and subleading jet $ p_{\mathrm{T}} $ in selected events. The uncertainty band in the lower panel represents the limited size of simulated samples together with a 30% uncertainty on the low mass Drell-Yan cross section. Simulated samples are normalized to 138 fb$^{-1}$ with the corresponding theoretical cross sections. The benchmark criteria explained in the text are used for signal normalization. |
png pdf |
Figure 2:
The $ p_{\mathrm{T}} $ distributions of the (left) dimuon systems, (right) di-b-jets system. The uncertainty band in the lower panel represents the limited size of simulated samples together with a 30% uncertainty on the low mass Drell-Yan cross section. Simulated samples are normalized to 138 fb$^{-1}$ with the corresponding theoretical cross sections. The benchmark criteria explained in the text are used for signal normalization. |
png pdf |
Figure 2-a:
The $ p_{\mathrm{T}} $ distributions of the (left) dimuon systems, (right) di-b-jets system. The uncertainty band in the lower panel represents the limited size of simulated samples together with a 30% uncertainty on the low mass Drell-Yan cross section. Simulated samples are normalized to 138 fb$^{-1}$ with the corresponding theoretical cross sections. The benchmark criteria explained in the text are used for signal normalization. |
png pdf |
Figure 2-b:
The $ p_{\mathrm{T}} $ distributions of the (left) dimuon systems, (right) di-b-jets system. The uncertainty band in the lower panel represents the limited size of simulated samples together with a 30% uncertainty on the low mass Drell-Yan cross section. Simulated samples are normalized to 138 fb$^{-1}$ with the corresponding theoretical cross sections. The benchmark criteria explained in the text are used for signal normalization. |
png pdf |
Figure 3:
The distribution of $ \chi_{\mathrm{b}\mathrm{b}} $ versus $ \chi_{\mathrm{H}} $ as defined in Eq. (1) for (left) simulated background processes, and (right) the signal process with $ m_{\rm a_1} = $ 40 GeV. The contours encircle the area with $ \chi_{\rm tot} $ below an arbitrary value. The grey scale represents the expected yields at 138 fb$^{-1}$. |
png pdf |
Figure 3-a:
The distribution of $ \chi_{\mathrm{b}\mathrm{b}} $ versus $ \chi_{\mathrm{H}} $ as defined in Eq. (1) for (left) simulated background processes, and (right) the signal process with $ m_{\rm a_1} = $ 40 GeV. The contours encircle the area with $ \chi_{\rm tot} $ below an arbitrary value. The grey scale represents the expected yields at 138 fb$^{-1}$. |
png pdf |
Figure 3-b:
The distribution of $ \chi_{\mathrm{b}\mathrm{b}} $ versus $ \chi_{\mathrm{H}} $ as defined in Eq. (1) for (left) simulated background processes, and (right) the signal process with $ m_{\rm a_1} = $ 40 GeV. The contours encircle the area with $ \chi_{\rm tot} $ below an arbitrary value. The grey scale represents the expected yields at 138 fb$^{-1}$. |
png pdf |
Figure 4:
Signal ($ m_{\rm a_1}= $ 40 GeV) versus background efficiency for different thresholds on $ \chi_{\rm tot}^2 $ (gray) and $ \chi_{\rm d}^2 $ (red) variables. The black star indicates signal efficiency versus that of background for the optimized $ \chi_{\rm d}^2 $ requirement. |
png pdf |
Figure 5:
Pre-fit distributions of the transformed DNN score for the $ \mu\,\tau_\mathrm{h} $ \ channel divided into events with one (left) or at least two (right) b jets. The shape of the $ \mathrm{H}\to{\rm a_1}{\rm a_1} $ signal, where $ m_{\rm a_1} = $ 35 GeV, is indicated assuming $ \mathrm{B}(\mathrm{H}\to{\rm a_1}{\rm a_1}\to\tau\tau{\mathrm{b}}{\mathrm{b}}) $ to be 100%. The lower panel shows the ratio of the observed data to the expected yields. The grey band represents the constrained statistical and systematic uncertainties. |
png pdf |
Figure 5-a:
Pre-fit distributions of the transformed DNN score for the $ \mu\,\tau_\mathrm{h} $ \ channel divided into events with one (left) or at least two (right) b jets. The shape of the $ \mathrm{H}\to{\rm a_1}{\rm a_1} $ signal, where $ m_{\rm a_1} = $ 35 GeV, is indicated assuming $ \mathrm{B}(\mathrm{H}\to{\rm a_1}{\rm a_1}\to\tau\tau{\mathrm{b}}{\mathrm{b}}) $ to be 100%. The lower panel shows the ratio of the observed data to the expected yields. The grey band represents the constrained statistical and systematic uncertainties. |
png pdf |
Figure 5-b:
Pre-fit distributions of the transformed DNN score for the $ \mu\,\tau_\mathrm{h} $ \ channel divided into events with one (left) or at least two (right) b jets. The shape of the $ \mathrm{H}\to{\rm a_1}{\rm a_1} $ signal, where $ m_{\rm a_1} = $ 35 GeV, is indicated assuming $ \mathrm{B}(\mathrm{H}\to{\rm a_1}{\rm a_1}\to\tau\tau{\mathrm{b}}{\mathrm{b}}) $ to be 100%. The lower panel shows the ratio of the observed data to the expected yields. The grey band represents the constrained statistical and systematic uncertainties. |
png pdf |
Figure 6:
The best-fit background models together with 68% CL uncertainty band from the fit to the data under the background-only hypothesis for the (top left) TT category, (top right) TM, (middle left) TL category, (middle right) Low $ p_{\mathrm{T}} $ category, and (bottom) VBF category. |
png pdf |
Figure 6-a:
The best-fit background models together with 68% CL uncertainty band from the fit to the data under the background-only hypothesis for the (top left) TT category, (top right) TM, (middle left) TL category, (middle right) Low $ p_{\mathrm{T}} $ category, and (bottom) VBF category. |
png pdf |
Figure 6-b:
The best-fit background models together with 68% CL uncertainty band from the fit to the data under the background-only hypothesis for the (top left) TT category, (top right) TM, (middle left) TL category, (middle right) Low $ p_{\mathrm{T}} $ category, and (bottom) VBF category. |
png pdf |
Figure 6-c:
The best-fit background models together with 68% CL uncertainty band from the fit to the data under the background-only hypothesis for the (top left) TT category, (top right) TM, (middle left) TL category, (middle right) Low $ p_{\mathrm{T}} $ category, and (bottom) VBF category. |
png pdf |
Figure 6-d:
The best-fit background models together with 68% CL uncertainty band from the fit to the data under the background-only hypothesis for the (top left) TT category, (top right) TM, (middle left) TL category, (middle right) Low $ p_{\mathrm{T}} $ category, and (bottom) VBF category. |
png pdf |
Figure 6-e:
The best-fit background models together with 68% CL uncertainty band from the fit to the data under the background-only hypothesis for the (top left) TT category, (top right) TM, (middle left) TL category, (middle right) Low $ p_{\mathrm{T}} $ category, and (bottom) VBF category. |
png pdf |
Figure 7:
Post-fit distributions of the $ m_{\tau\tau} $ for the $ \mu\,\tau_\mathrm{h} $ \ channel signal regions in events with exactly one b-tagged jet: SR1 (top left), SR2 (top right), and SR3 (bottom). The shape of the $ \mathrm{H}\to{\rm a_1}{\rm a_1} $ signal, where $ m_{\rm a_1} = $ 35 GeV, is indicated assuming $ \mathrm{B}(\mathrm{H}\to{\rm a_1}{\rm a_1}\to\tau\tau{\mathrm{b}}{\mathrm{b}}) $ to be 10%. |
png pdf |
Figure 7-a:
Post-fit distributions of the $ m_{\tau\tau} $ for the $ \mu\,\tau_\mathrm{h} $ \ channel signal regions in events with exactly one b-tagged jet: SR1 (top left), SR2 (top right), and SR3 (bottom). The shape of the $ \mathrm{H}\to{\rm a_1}{\rm a_1} $ signal, where $ m_{\rm a_1} = $ 35 GeV, is indicated assuming $ \mathrm{B}(\mathrm{H}\to{\rm a_1}{\rm a_1}\to\tau\tau{\mathrm{b}}{\mathrm{b}}) $ to be 10%. |
png pdf |
Figure 7-b:
Post-fit distributions of the $ m_{\tau\tau} $ for the $ \mu\,\tau_\mathrm{h} $ \ channel signal regions in events with exactly one b-tagged jet: SR1 (top left), SR2 (top right), and SR3 (bottom). The shape of the $ \mathrm{H}\to{\rm a_1}{\rm a_1} $ signal, where $ m_{\rm a_1} = $ 35 GeV, is indicated assuming $ \mathrm{B}(\mathrm{H}\to{\rm a_1}{\rm a_1}\to\tau\tau{\mathrm{b}}{\mathrm{b}}) $ to be 10%. |
png pdf |
Figure 7-c:
Post-fit distributions of the $ m_{\tau\tau} $ for the $ \mu\,\tau_\mathrm{h} $ \ channel signal regions in events with exactly one b-tagged jet: SR1 (top left), SR2 (top right), and SR3 (bottom). The shape of the $ \mathrm{H}\to{\rm a_1}{\rm a_1} $ signal, where $ m_{\rm a_1} = $ 35 GeV, is indicated assuming $ \mathrm{B}(\mathrm{H}\to{\rm a_1}{\rm a_1}\to\tau\tau{\mathrm{b}}{\mathrm{b}}) $ to be 10%. |
png pdf |
Figure 8:
Post-fit distributions of the $ m_{\tau\tau} $ for the $ \mu\,\tau_\mathrm{h} $ \ channel signal regions in events with at least two b-tagged jets: SR1 (left) and SR2 (right). The shape of the $ \mathrm{H}\to{\rm a_1}{\rm a_1} $ signal, where $ m_{\rm a_1} = $ 35 GeV, is indicated assuming $ \mathrm{B}(\mathrm{H}\to{\rm a_1}{\rm a_1}\to\tau\tau{\mathrm{b}}{\mathrm{b}}) $ to be 10%. |
png pdf |
Figure 8-a:
Post-fit distributions of the $ m_{\tau\tau} $ for the $ \mu\,\tau_\mathrm{h} $ \ channel signal regions in events with at least two b-tagged jets: SR1 (left) and SR2 (right). The shape of the $ \mathrm{H}\to{\rm a_1}{\rm a_1} $ signal, where $ m_{\rm a_1} = $ 35 GeV, is indicated assuming $ \mathrm{B}(\mathrm{H}\to{\rm a_1}{\rm a_1}\to\tau\tau{\mathrm{b}}{\mathrm{b}}) $ to be 10%. |
png pdf |
Figure 8-b:
Post-fit distributions of the $ m_{\tau\tau} $ for the $ \mu\,\tau_\mathrm{h} $ \ channel signal regions in events with at least two b-tagged jets: SR1 (left) and SR2 (right). The shape of the $ \mathrm{H}\to{\rm a_1}{\rm a_1} $ signal, where $ m_{\rm a_1} = $ 35 GeV, is indicated assuming $ \mathrm{B}(\mathrm{H}\to{\rm a_1}{\rm a_1}\to\tau\tau{\mathrm{b}}{\mathrm{b}}) $ to be 10%. |
png pdf |
Figure 9:
Post-fit distributions of the $ m_{\tau\tau} $ for the $ \mu\,\tau_\mathrm{h} $ \ channel control regions in events with exactly one b-tagged jet (left) and at least two b-tagged jets (right). The shape of the $ \mathrm{H}\to{\rm a_1}{\rm a_1} $ signal, where $ m_{\rm a_1} = $ 35 GeV, is indicated assuming $ \mathrm{B}(\mathrm{H}\to{\rm a_1}{\rm a_1}\to\tau\tau{\mathrm{b}}{\mathrm{b}}) $ to be 10%. |
png pdf |
Figure 9-a:
Post-fit distributions of the $ m_{\tau\tau} $ for the $ \mu\,\tau_\mathrm{h} $ \ channel control regions in events with exactly one b-tagged jet (left) and at least two b-tagged jets (right). The shape of the $ \mathrm{H}\to{\rm a_1}{\rm a_1} $ signal, where $ m_{\rm a_1} = $ 35 GeV, is indicated assuming $ \mathrm{B}(\mathrm{H}\to{\rm a_1}{\rm a_1}\to\tau\tau{\mathrm{b}}{\mathrm{b}}) $ to be 10%. |
png pdf |
Figure 9-b:
Post-fit distributions of the $ m_{\tau\tau} $ for the $ \mu\,\tau_\mathrm{h} $ \ channel control regions in events with exactly one b-tagged jet (left) and at least two b-tagged jets (right). The shape of the $ \mathrm{H}\to{\rm a_1}{\rm a_1} $ signal, where $ m_{\rm a_1} = $ 35 GeV, is indicated assuming $ \mathrm{B}(\mathrm{H}\to{\rm a_1}{\rm a_1}\to\tau\tau{\mathrm{b}}{\mathrm{b}}) $ to be 10%. |
png pdf |
Figure 10:
Observed and expected upper limits at 95% CL on $ \mathrm{B}(\mathrm{H}\to{\rm a_1}{\rm a_1}\to \mu\mu{\mathrm{b}}{\mathrm{b}}) $ as a function of $ m_{\rm a_1} $. The inner and outer bands indicate the regions containing the distribution of limits located within 68 and 95% confidence intervals, respectively, of the expectation under the background-only hypothesis. |
png pdf |
Figure 11:
Observed and expected 95% CL exclusion limits on $ \mathrm{B}(\mathrm{H}\to{\rm a_1}{\rm a_1}\to\tau\tau{\mathrm{b}}{\mathrm{b}}) $ in %, for the combination of all years with an integrated luminosity of 138 fb$^{-1}$ per channel and the combination. Top left: $ \mu\,\tau_\mathrm{h} $, top right: $ \mathrm{e}\,\tau_\mathrm{h} $, bottom left: $ \mathrm{e}\,\mu $, and bottom right: combination of all channels is shown. |
png pdf |
Figure 11-a:
Observed and expected 95% CL exclusion limits on $ \mathrm{B}(\mathrm{H}\to{\rm a_1}{\rm a_1}\to\tau\tau{\mathrm{b}}{\mathrm{b}}) $ in %, for the combination of all years with an integrated luminosity of 138 fb$^{-1}$ per channel and the combination. Top left: $ \mu\,\tau_\mathrm{h} $, top right: $ \mathrm{e}\,\tau_\mathrm{h} $, bottom left: $ \mathrm{e}\,\mu $, and bottom right: combination of all channels is shown. |
png pdf |
Figure 11-b:
Observed and expected 95% CL exclusion limits on $ \mathrm{B}(\mathrm{H}\to{\rm a_1}{\rm a_1}\to\tau\tau{\mathrm{b}}{\mathrm{b}}) $ in %, for the combination of all years with an integrated luminosity of 138 fb$^{-1}$ per channel and the combination. Top left: $ \mu\,\tau_\mathrm{h} $, top right: $ \mathrm{e}\,\tau_\mathrm{h} $, bottom left: $ \mathrm{e}\,\mu $, and bottom right: combination of all channels is shown. |
png pdf |
Figure 11-c:
Observed and expected 95% CL exclusion limits on $ \mathrm{B}(\mathrm{H}\to{\rm a_1}{\rm a_1}\to\tau\tau{\mathrm{b}}{\mathrm{b}}) $ in %, for the combination of all years with an integrated luminosity of 138 fb$^{-1}$ per channel and the combination. Top left: $ \mu\,\tau_\mathrm{h} $, top right: $ \mathrm{e}\,\tau_\mathrm{h} $, bottom left: $ \mathrm{e}\,\mu $, and bottom right: combination of all channels is shown. |
png pdf |
Figure 11-d:
Observed and expected 95% CL exclusion limits on $ \mathrm{B}(\mathrm{H}\to{\rm a_1}{\rm a_1}\to\tau\tau{\mathrm{b}}{\mathrm{b}}) $ in %, for the combination of all years with an integrated luminosity of 138 fb$^{-1}$ per channel and the combination. Top left: $ \mu\,\tau_\mathrm{h} $, top right: $ \mathrm{e}\,\tau_\mathrm{h} $, bottom left: $ \mathrm{e}\,\mu $, and bottom right: combination of all channels is shown. |
png pdf |
Figure 12:
Observed and expected 95% CL upper limits on $ \mathrm{B}(\mathrm{H}\to{\rm a_1}{\rm a_1}\to\ell\ell\mathrm{b}\mathrm{b}) $ in %, where $ \ell $ stands for muons or tau leptons, obtained from the combination of $ \mu\mu{\mathrm{b}}{\mathrm{b}} $ and $ \tau\tau{\mathrm{b}}{\mathrm{b}} $ channels using the full Run 2 integrated luminosity of 138 fb$^{-1}$. The results are obtained as functions of $ m_{\rm a_1} $ for 2HDM+S models, independent of the type and $ \tan\beta $ parameter. |
png pdf |
Figure 13:
Observed 95% CL upper limits on $ \mathrm{B}(\mathrm{H}\to{\rm a_1}{\rm a_1}) $ in %, for the combination of $ \mu\mu{\mathrm{b}}{\mathrm{b}} $ and $ \tau\tau{\mathrm{b}}{\mathrm{b}} $ channels using the full Run 2 integrated luminosity of 138 fb$^{-1}$ for Type III (left) and Type IV (right) 2HDM+S in the $ \tan\beta $ vs. $ m_{\rm a_1} $ phase space. The contours corresponding to branching fractions of 100% and 16% are drawn using dashed lines where 16% refers to the combined upper limit on Higgs to BSM particle decays from previous Run 2 results [10]. Linear extrapolation has been used between different points on the figures. |
png pdf |
Figure 13-a:
Observed 95% CL upper limits on $ \mathrm{B}(\mathrm{H}\to{\rm a_1}{\rm a_1}) $ in %, for the combination of $ \mu\mu{\mathrm{b}}{\mathrm{b}} $ and $ \tau\tau{\mathrm{b}}{\mathrm{b}} $ channels using the full Run 2 integrated luminosity of 138 fb$^{-1}$ for Type III (left) and Type IV (right) 2HDM+S in the $ \tan\beta $ vs. $ m_{\rm a_1} $ phase space. The contours corresponding to branching fractions of 100% and 16% are drawn using dashed lines where 16% refers to the combined upper limit on Higgs to BSM particle decays from previous Run 2 results [10]. Linear extrapolation has been used between different points on the figures. |
png pdf |
Figure 13-b:
Observed 95% CL upper limits on $ \mathrm{B}(\mathrm{H}\to{\rm a_1}{\rm a_1}) $ in %, for the combination of $ \mu\mu{\mathrm{b}}{\mathrm{b}} $ and $ \tau\tau{\mathrm{b}}{\mathrm{b}} $ channels using the full Run 2 integrated luminosity of 138 fb$^{-1}$ for Type III (left) and Type IV (right) 2HDM+S in the $ \tan\beta $ vs. $ m_{\rm a_1} $ phase space. The contours corresponding to branching fractions of 100% and 16% are drawn using dashed lines where 16% refers to the combined upper limit on Higgs to BSM particle decays from previous Run 2 results [10]. Linear extrapolation has been used between different points on the figures. |
| Tables | |
png pdf |
Table 1:
The electron, muon, and $ \tau_\mathrm{h} p_{\mathrm{T}} $ thresholds in GeV at trigger level for the $ \tau\tau{\mathrm{b}}{\mathrm{b}} $ and $ \mu\mu{\mathrm{b}}{\mathrm{b}} $ final states. |
png pdf |
Table 2:
Event yields for simulated processes and the number of observed events in data after applying $ \chi_{\rm d}^2 < $ 1.5. The expected number of simulated events is normalized to the integrated luminosity of 138 fb$^{-1}$. The type-III parametrization of 2HDM+S with $ \tan\beta= $ 2 is used to evaluate Br($ {\rm a_1}\to ff $). |
png pdf |
Table 3:
Summary of categorization requirements. Events in these categories contain two muons and two b-jets. As stated in the text, L, M, and T respectively stand for the loose, medium, and tight b-tag criteria. |
png pdf |
Table 4:
The expected yields for backgrounds and different signal hypotheses in each category. The entries are rounded to first decimal place. |
png pdf |
Table 5:
Event categories and subregions for the $ \tau\tau{\mathrm{b}}{\mathrm{b}} $ channel. The values correspond to the transformed DNN score used to define the signal (SRn) and control (CR) regions. |
| Summary |
| A search for an exotic decay of the 125 GeV Higgs boson to a pair of light pseudoscalar bosons in the final state with two b quarks and two $ \tau $ leptons or muons has been presented. The results are based on a data sample of proton-proton collisions corresponding to an integrated luminosity of 138 fb$^{-1}$, accumulated by the CMS experiment during LHC Run 2 at a center-of-mass energy of 13 TeV. Final states with at least one leptonic $ \tau $ decay are studied in $ \tau\tau{\mathrm{b}}{\mathrm{b}} $, excluding those with two muons or two electrons. The results show significant improvement with respect to the earlier CMS analyses at 13\, TeV, beyond what is merely expected from the increase in the size of the data sample. The new $ \tau\tau{\mathrm{b}}{\mathrm{b}} $ analysis gains from the deep neural network based signal categorization, while a more thorough analysis of the signal properties using a single discrimination variable improves $ \mu\mu{\mathrm{b}}{\mathrm{b}} $. In the absence of any significant excess in the data over the standard model backgrounds, upper limits are set, at 95% confidence level on $ \mathrm{B}(\mathrm{H}\to{\rm a_1}{\rm a_1}\to\mu\mu{\mathrm{b}}{\mathrm{b}}) $ and $ \mathrm{B}(\mathrm{H}\to{\rm a_1}{\rm a_1}\to\tau\tau{\mathrm{b}}{\mathrm{b}}) $, in the $ \mu\mu{\mathrm{b}}{\mathrm{b}} $ and $ \tau\tau{\mathrm{b}}{\mathrm{b}} $ analyses respectively. Both analyses provide the most stringent expected limits to date. In $ \mu\mu{\mathrm{b}}{\mathrm{b}} $, the observed limits are in the range (0.35-2.6) $\times$ 10$^{-4} $ for a pseudoscalar mass, $ m_{\rm a_1} $, between 15 and 62.5 GeV. Combining all final states in $ \tau\tau{\mathrm{b}}{\mathrm{b}} $, limits are observed to be in range (1.8--7.7)% for $ m_{\rm a_1} $ between 12 and 60 GeV. In the context of 2HDM+S models, the allowed values of the branching fraction $ \mathrm{B}(\mathrm{H}\to{\rm a_1}{\rm a_1}) $ are obtained by combining the $ \mu\mu{\mathrm{b}}{\mathrm{b}} $ and $ \tau\tau{\mathrm{b}}{\mathrm{b}} $ channels. For $ m_{\rm a_1} $ values between 15 and 60 GeV, $ \mathrm{B}(\mathrm{H}\to{\rm a_1}{\rm a_1}) $ above 23% are excluded, at 95% confidence level, in most of the Type-II models. In Type III and IV, upper limits as low as about 2% and 3% are obtained, respectively, for $ \tan\beta= $ 2.0 and 0.5. |
| 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 Symmetry and the Mass of Gauge Vector Mesons | PRL 13 (1964) 321 | |
| 5 | G. C. Branco et al. | Theory and phenomenology of two-Higgs-doublet models | Phys. Rept. 516 (2012) 1 | 1106.0034 |
| 6 | A. Djouadi | The Anatomy of electro-weak symmetry breaking. II. The Higgs bosons in the minimal supersymmetric model | Phys. Rept. 459 (2008) 1 | hep-ph/0503173 |
| 7 | T. Robens and T. Stefaniak | Status of the Higgs Singlet Extension of the Standard Model after LHC Run 1 | EPJC 75 (2015) 104 | 1501.02234 |
| 8 | T. Robens, T. Stefaniak, and J. Wittbrodt | Two-real-scalar-singlet extension of the SM: LHC phenomenology and benchmark scenarios | EPJC 80 (2020) 151 | 1908.08554 |
| 9 | ATLAS Collaboration | A detailed map of Higgs boson interactions by the ATLAS experiment ten years after the discovery | Nature 607 (2022) 52 | 2207.00092 |
| 10 | CMS Collaboration | A portrait of the Higgs boson by the CMS experiment ten years after the discovery | Nature 607 (2022) 60 | CMS-HIG-22-001 2207.00043 |
| 11 | B. Grzadkowski and P. Osland | Tempered Two-Higgs-Doublet Model | PRD 82 (2010) 125026 | 0910.4068 |
| 12 | A. Drozd, B. Grzadkowski, J. F. Gunion, and Y. Jiang | Extending two-Higgs-doublet models by a singlet scalar field - the Case for Dark Matter | JHEP 11 (2014) 105 | 1408.2106 |
| 13 | S. Ramos-Sanchez | The $ \mu $-problem, the NMSSM and string theory | Fortschritte der Physik 58 (2010) 748 | 1003.1307 |
| 14 | D. Curtin et al. | Exotic decays of the 125 GeV Higgs boson | PRD 90 (2014) 075004 | 1312.4992 |
| 15 | LHC Higgs Cross Section Working Group Collaboration | Handbook of LHC Higgs Cross Sections: 4. Deciphering the Nature of the Higgs Sector | technical report, 2016 link |
1610.07922 |
| 16 | ATLAS Collaboration | Search for Higgs boson decays into a pair of pseudoscalar particles in the $ bb\mu\mu $ final state with the ATLAS detector in pp collisions at $ \sqrt{s}= $ 13 TeV | PRD 105 (2022) | 2110.00313 |
| 17 | ATLAS Collaboration | Search for Higgs boson decays into a pair of light bosons in the $ \mathrm{b}\mathrm{b}\mu\mu $ final state in pp collision at $ \sqrt{s} = $ 13 TeV with the ATLAS detector | PLB 790 (2019) 1 | 1807.00539 |
| 18 | CMS Collaboration | Search for an exotic decay of the Higgs boson to a pair of light pseudoscalars in the final state with two muons and two b quarks in pp collisions at 13 TeV | PLB 795 (2019) | CMS-HIG-18-011 1812.06359 |
| 19 | CMS Collaboration | Search for light bosons in decays of the 125 GeV Higgs boson in proton-proton collisions at $ \sqrt{s}= $ 8 TeV | JHEP 10 (2017) 076 | CMS-HIG-16-015 1701.02032 |
| 20 | CMS Collaboration | Search for an exotic decay of the Higgs boson to a pair of light pseudoscalars in the final state with two b quarks and two $ \tau $ leptons in proton-proton collisions at $ \sqrt{s}= $ 13 TeV | PLB 785 (2018) 462 | CMS-HIG-17-024 1805.10191 |
| 21 | CMS Collaboration | The CMS trigger system | JINST 12 (2017) P01020 | CMS-TRG-12-001 1609.02366 |
| 22 | CMS Collaboration | The CMS Experiment at the CERN LHC | JINST 3 (2008) S08004 | |
| 23 | T. Sjöstrand et al. | An Introduction to PYTHIA 8.2 | Comput. Phys. Commun. 191 (2015) 159 | 1410.3012 |
| 24 | CMS Collaboration | Measurements of differential production cross sections for a Z boson in association with jets in pp collisions at $ \sqrt{s} = $ 8 TeV | JHEP 04 (2017) 022 | CMS-SMP-14-013 1611.03844 |
| 25 | R. Frederix and S. Frixione | Merging meets matching in MC@NLO | JHEP 12 (2012) 061 | 1209.6215 |
| 26 | J. Alwall et al. | The automated computation of tree-level and next-to-leading order differential cross sections, and their matching to parton shower simulations | JHEP 07 (2014) 079 | 1405.0301 |
| 27 | P. Skands, S. Carrazza, and J. Rojo | Tuning PYTHIA 8.1: the Monash 2013 Tune | EPJC 74 (2014) 3024 | 1404.5630 |
| 28 | 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 |
| 29 | J. Allison et al. | GEANT4 developments and applications | IEEE Trans. Nucl. Sci. 53 (2006) 270 | |
| 30 | 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 |
| 31 | P. Nason | A New method for combining NLO QCD with shower Monte Carlo algorithms | JHEP 11 (2004) 040 | hep-ph/0409146 |
| 32 | S. Frixione, P. Nason, and C. Oleari | Matching NLO QCD computations with Parton Shower simulations: the POWHEG method | JHEP 11 (2007) 070 | 0709.2092 |
| 33 | 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 |
| 34 | S. Alioli et al. | Jet pair production in POWHEG | JHEP 04 (2011) 081 | 1012.3380 |
| 35 | M. Czakon et al. | Top-pair production at the LHC through NNLO QCD and NLO EW | JHEP 10 (2017) 186 | 1705.04105 |
| 36 | S. Alioli, P. Nason, C. Oleari, and E. Re | NLO Higgs boson production via gluon fusion matched with shower in POWHEG | JHEP 04 (2009) 002 | 0812.0578 |
| 37 | E. Bagnaschi, G. Degrassi, P. Slavich, and A. Vicini | Higgs production via gluon fusion in the POWHEG approach in the SM and in the MSSM | JHEP 02 (2012) 088 | 1111.2854 |
| 38 | P. Nason and C. Oleari | NLO Higgs boson production via vector-boson fusion matched with shower in POWHEG | JHEP 02 (2010) 037 | 0911.5299 |
| 39 | G. Luisoni, P. Nason, C. Oleari, and F. Tramontano | $ \mathrm{H}\mathrm{W}^{\pm} $/$HZ$ + 0 and 1 jet at NLO with the POWHEG BOX interfaced to GoSam and their merging within MiNLO | JHEP 10 (2013) 083 | 1306.2542 |
| 40 | H. B. Hartanto, B. Jager, L. Reina, and D. Wackeroth | Higgs boson production in association with top quarks in the POWHEG BOX | PRD 91 (2015) 094003 | 1501.04498 |
| 41 | CMS Collaboration | Particle-flow reconstruction and global event description with the CMS detector | JINST 12 (2017) P10003 | CMS-PRF-14-001 1706.04965 |
| 42 | CMS Collaboration | Technical proposal for the Phase-II upgrade of the Compact Muon Solenoid | CMS Technical Proposal CERN-LHCC-2015-010, CMS-TDR-15-02, 2015 CDS |
|
| 43 | 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 |
| 44 | CMS Collaboration | Electron and photon reconstruction and identification with the CMS experiment at the CERN LHC | JINST 16 (2021) P05014 | CMS-EGM-17-001 2012.06888 |
| 45 | M. Cacciari, G. P. Salam, and G. Soyez | The anti-$ k_{\mathrm{T}} $ jet clustering algorithm | JHEP 04 (2008) 063 | 0802.1189 |
| 46 | M. Cacciari, G. P. Salam, and G. Soyez | FastJet User Manual | EPJC 72 (2012) 1896 | 1111.6097 |
| 47 | CMS Collaboration | Jet algorithms performance in 13 TeV data | CMS Physics Analysis Summary, 2017 CMS-PAS-JME-16-003 |
CMS-PAS-JME-16-003 |
| 48 | CMS Collaboration | Determination of Jet Energy Calibration and Transverse Momentum Resolution in CMS | JINST 6 (2011) P11002 | CMS-JME-10-011 1107.4277 |
| 49 | E. Bols et al. | Jet Flavour Classification Using DeepJet | JINST 15 (2020) P12012 | 2008.10519 |
| 50 | 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 |
| 51 | CMS Collaboration | Performance of the DeepJet b tagging algorithm using 41.9/fb of data from proton-proton collisions at 13 TeV with Phase 1 CMS detector | CMS Detector Performance Note CMS-DP-2018-058, 2018 CDS |
|
| 52 | CMS Collaboration | Performance of reconstruction and identification of $ \tau $ leptons decaying to hadrons and $ \nu_\tau $ in pp collisions at $ \sqrt{s}= $ 13 TeV | JINST 13 (2018) P10005 | CMS-TAU-16-003 1809.02816 |
| 53 | CMS Collaboration | Identification of hadronic tau lepton decays using a deep neural network | JINST 17 (2022) P07023 | CMS-TAU-20-001 2201.08458 |
| 54 | 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 |
| 55 | CDF Collaboration | Search for neutral Higgs bosons of the minimal supersymmetric standard model decaying to $ \tau $ pairs in $ \mathrm{p}\bar{\mathrm{p}} $ collisions at $ \sqrt{s} = $ 1.96 TeV | PRL 96 (2006) 011802 | hep-ex/0508051 |
| 56 | L. Bianchini et al. | Reconstruction of the Higgs mass in events with Higgs bosons decaying into a pair of tau leptons using matrix element technique | NIM A 862 (2017) 54 | 1603.05910 |
| 57 | CMS Collaboration | An embedding technique to determine $ \tau\tau $ backgrounds in proton-proton collision data | JINST 14 (2019) P06032 | CMS-TAU-18-001 1903.01216 |
| 58 | 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 |
| 59 | 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 |
| 60 | ATLAS, CMS Collaboration | 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 |
| 61 | ATLAS and CMS Collaborations, and LHC Higgs Combination Group | Procedure for the LHC Higgs boson search combination in Summer 2011 | Technical Report CMS-NOTE-2011-005, ATL-PHYS-PUB-2011-11 | |
| 62 | T. Junk | Confidence level computation for combining searches with small statistics | NIM A 434 (1999) 435 | hep-ex/9902006 |
| 63 | A. L. Read | Presentation of search results: The CL$ _{\text{s}} $ technique | JPG 28 (2002) 2693 | |
| 64 | G. Cowan, K. Cranmer, E. Gross, and O. Vitells | Asymptotic formulae for likelihood-based tests of new physics | EPJC 71 (2011) 1554 | 1007.1727 |
| 65 | U. Haisch, J. F. Kamenik, A. Malinauskas, and M. Spira | Collider constraints on light pseudoscalars | JHEP 03 (2018) 178 | 1802.02156 |
| 66 | CMS Collaboration | Precision luminosity measurement in proton-proton collisions at $ \sqrt{s} = $ 13 TeV in 2015 and 2016 at CMS | EPJC 81 (2021) 800 | CMS-LUM-17-003 2104.01927 |
| 67 | CMS Collaboration | CMS luminosity measurement for the 2017 data-taking period at $ \sqrt{s}= $ 13 TeV | CMS Physics Analysis Summary, 2018 CMS-PAS-LUM-17-004 |
CMS-PAS-LUM-17-004 |
| 68 | CMS Collaboration | CMS luminosity measurement for the 2018 data-taking period at $ \sqrt{s}= $ 13 TeV | CMS Physics Analysis Summary, 2019 CMS-PAS-LUM-18-002 |
CMS-PAS-LUM-18-002 |
| 69 | 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 |
| 70 | 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 |
| 71 | CMS Collaboration | Measurement of the inclusive and differential Higgs boson production cross sections in the decay mode to a pair of $ \tau $ leptons in pp collisions at $ \sqrt{s} = $ 13 TeV | PRL 128 (2022) 081805 | CMS-HIG-20-015 2107.11486 |
| 72 | CMS Collaboration | Performance of electron reconstruction and selection with the CMS detector in proton-proton collisions at $ \sqrt{s} = $ 8 TeV | JINST 10 (2015) P06005 | CMS-EGM-13-001 1502.02701 |
| 73 | R. J. Barlow and C. Beeston | Fitting using finite Monte Carlo samples | Comput. Phys. Commun. 77 (1993) 219 | |
| 74 | J. Butterworth et al. | PDF4LHC recommendations for LHC Run II | JPG 43 (2016) 023001 | 1510.03865 |
|
|
Compact Muon Solenoid LHC, CERN |
|
|
|
|
|
|