CMS logoCMS event Hgg
Compact Muon Solenoid
LHC, CERN

CMS-PAS-B2G-23-006
Search for heavy neutral Higgs bosons A and H in the $ \mathrm{t\bar{t}} $Z final state
Abstract: A direct search for heavy neutral Higgs bosons A and H in the $ \mathrm{t\bar{t}} $Z final state is presented, targeting the process $ \mathrm{A}\rightarrow\mathrm{Z}\mathrm{H} $ with $ \mathrm{H}\rightarrow\mathrm{t\bar{t}} $. For the first time, the final state with decays of the Z boson to electrons or muons together with all-jet decays of the $ \mathrm{t\bar{t}} $ system is considered. The analysis is conducted with proton-proton collision data collected at the CERN LHC with the CMS experiment at $ \sqrt{s}= $ 13 TeV, which correspond to an integrated luminosity of 138 fb$ ^{-1} $. No signal is observed. Model-independent upper limits on the cross section times branching ratio are derived for A and H boson masses up to 2100 GeV and 1900 GeV, respectively, assuming narrow width for both A and H. The results are also interpreted within two-Higgs-doublet models. Regions in parameter space with A boson masses between 550 and up to 1500 GeV, and H boson masses between 350 and up to 700 GeV, depending on the model parameters, are excluded at 95% confidence level, complementing and substantially extending the reach of previous searches.
Figures Summary References CMS Publications
Figures

png pdf
Figure 1:
Definition of the event categories per lepton-flavour channel ($ \mathrm{e}^+\mathrm{e}^- $ or $ \mu^{+}\mu^{-} $) and jet multiplicity $ n $ (5 or $ \geq $6).

png pdf
Figure 2:
Expected $ (\Delta m,p^{\mathrm{Z}}_{\text{T}}) $ distribution for signal events with $ m_{\mathrm{A}}= $ 1000 GeV and $ m_{\mathrm{H}}= $ 600 GeV (blue shading), and boundaries of the elliptical bins (solid lines) used to construct the final observable.

png pdf
Figure 3:
Distribution of events in the $ p^{\mathrm{Z}}_{\text{T}}\times\Delta m $ analysis bins in the CR after the fit to data described in Section 6 with a $ (m_{\mathrm{A}},m_{\mathrm{H}}) $ signal hypothesis of (1000, 350) GeV (left) and (1000,850) GeV (right). The six bins in each region correspond to the six quantiles of the $ p^{\mathrm{Z}}_{\text{T}}\times\Delta m $ distribution, as described in the text. The hatched area represents the total post-fit uncertainty. The signal (solid red line) and background (coloured histograms) distributions are shown with their best fit normalisations from the simultaneous fit to the data (``post-fit''), where the yields in the $ \mathrm{e}\mathrm{e} $ and $ \mu\mu $ channels are added for representation purposes. In the left plot the post-fit signal cross section is found to be 0.0 $ \pm $ 0.1 fb, as indicated in the legend, and therefore the corresponding histogram is not displayed. The signal is also shown for a normalisation to 25 fb (``pre-fit'') in the dashed red line.

png pdf
Figure 3-a:
Distribution of events in the $ p^{\mathrm{Z}}_{\text{T}}\times\Delta m $ analysis bins in the CR after the fit to data described in Section 6 with a $ (m_{\mathrm{A}},m_{\mathrm{H}}) $ signal hypothesis of (1000, 350) GeV (left) and (1000,850) GeV (right). The six bins in each region correspond to the six quantiles of the $ p^{\mathrm{Z}}_{\text{T}}\times\Delta m $ distribution, as described in the text. The hatched area represents the total post-fit uncertainty. The signal (solid red line) and background (coloured histograms) distributions are shown with their best fit normalisations from the simultaneous fit to the data (``post-fit''), where the yields in the $ \mathrm{e}\mathrm{e} $ and $ \mu\mu $ channels are added for representation purposes. In the left plot the post-fit signal cross section is found to be 0.0 $ \pm $ 0.1 fb, as indicated in the legend, and therefore the corresponding histogram is not displayed. The signal is also shown for a normalisation to 25 fb (``pre-fit'') in the dashed red line.

png pdf
Figure 3-b:
Distribution of events in the $ p^{\mathrm{Z}}_{\text{T}}\times\Delta m $ analysis bins in the CR after the fit to data described in Section 6 with a $ (m_{\mathrm{A}},m_{\mathrm{H}}) $ signal hypothesis of (1000, 350) GeV (left) and (1000,850) GeV (right). The six bins in each region correspond to the six quantiles of the $ p^{\mathrm{Z}}_{\text{T}}\times\Delta m $ distribution, as described in the text. The hatched area represents the total post-fit uncertainty. The signal (solid red line) and background (coloured histograms) distributions are shown with their best fit normalisations from the simultaneous fit to the data (``post-fit''), where the yields in the $ \mathrm{e}\mathrm{e} $ and $ \mu\mu $ channels are added for representation purposes. In the left plot the post-fit signal cross section is found to be 0.0 $ \pm $ 0.1 fb, as indicated in the legend, and therefore the corresponding histogram is not displayed. The signal is also shown for a normalisation to 25 fb (``pre-fit'') in the dashed red line.

png pdf
Figure 4:
Distributions of $ p^{\mathrm{Z}}_{\text{T}}\times\Delta m $ in the SR after the fit to data with a $ (m_{\mathrm{A}},m_{\mathrm{H}}) $ signal hypothesis of (1000, 350) GeV (left) and (1000, 850) GeV (right). The hatched area represents the total uncertainty. The signal (solid red line) and background (coloured histograms) distributions are shown with their best fit normalisations from the simultaneous fit to the data (``post-fit''), where the yields in the $ \mathrm{e}\mathrm{e} $ and $ \mu\mu $ channels are added for representation purposes. In the left plot the post-fit signal cross section is found to be 0.0 $ \pm $ 0.1 fb, as indicated in the legend, and therefore the corresponding histogram is not displayed. The signal is also shown for a normalisation to 25 fb (``pre-fit'') in the dashed red line.

png pdf
Figure 4-a:
Distributions of $ p^{\mathrm{Z}}_{\text{T}}\times\Delta m $ in the SR after the fit to data with a $ (m_{\mathrm{A}},m_{\mathrm{H}}) $ signal hypothesis of (1000, 350) GeV (left) and (1000, 850) GeV (right). The hatched area represents the total uncertainty. The signal (solid red line) and background (coloured histograms) distributions are shown with their best fit normalisations from the simultaneous fit to the data (``post-fit''), where the yields in the $ \mathrm{e}\mathrm{e} $ and $ \mu\mu $ channels are added for representation purposes. In the left plot the post-fit signal cross section is found to be 0.0 $ \pm $ 0.1 fb, as indicated in the legend, and therefore the corresponding histogram is not displayed. The signal is also shown for a normalisation to 25 fb (``pre-fit'') in the dashed red line.

png pdf
Figure 4-b:
Distributions of $ p^{\mathrm{Z}}_{\text{T}}\times\Delta m $ in the SR after the fit to data with a $ (m_{\mathrm{A}},m_{\mathrm{H}}) $ signal hypothesis of (1000, 350) GeV (left) and (1000, 850) GeV (right). The hatched area represents the total uncertainty. The signal (solid red line) and background (coloured histograms) distributions are shown with their best fit normalisations from the simultaneous fit to the data (``post-fit''), where the yields in the $ \mathrm{e}\mathrm{e} $ and $ \mu\mu $ channels are added for representation purposes. In the left plot the post-fit signal cross section is found to be 0.0 $ \pm $ 0.1 fb, as indicated in the legend, and therefore the corresponding histogram is not displayed. The signal is also shown for a normalisation to 25 fb (``pre-fit'') in the dashed red line.

png pdf
Figure 5:
Expected (left) and observed (right) 95% CL upper limits on the production cross section times branching ratio of the $ \mathrm{A}\to\mathrm{Z}\mathrm{H}\to\mathrm{Z}{\mathrm{t}\overline{\mathrm{t}}} $ process in the $ (m_{\mathrm{A}}, m_{\mathrm{H}}) $ plane.

png pdf
Figure 5-a:
Expected (left) and observed (right) 95% CL upper limits on the production cross section times branching ratio of the $ \mathrm{A}\to\mathrm{Z}\mathrm{H}\to\mathrm{Z}{\mathrm{t}\overline{\mathrm{t}}} $ process in the $ (m_{\mathrm{A}}, m_{\mathrm{H}}) $ plane.

png pdf
Figure 5-b:
Expected (left) and observed (right) 95% CL upper limits on the production cross section times branching ratio of the $ \mathrm{A}\to\mathrm{Z}\mathrm{H}\to\mathrm{Z}{\mathrm{t}\overline{\mathrm{t}}} $ process in the $ (m_{\mathrm{A}}, m_{\mathrm{H}}) $ plane.

png pdf
Figure 6:
Median expected (dashed lines) and observed (filled contours) 95% CL exclusion regions in the ($ m_{\mathrm{H}} $,$ m_{\mathrm{A}} $) parameter space of the type-II 2HDM for $ \tan\beta = $ 0.5 (blue), 1 (orange), and 2 (green). The enclosed regions are excluded. The dash-dotted lines bound the regions where the width of the A boson is up to 25% relative to $ m_{\mathrm{A}} $.

png pdf
Figure 7:
Median expected (dashed lines) and observed (filled contours) 95% CL exclusion regions in the type-II 2HDM ($ \tan\beta $, $ m_{\mathrm{A}} $) parameter space at $ m_{\mathrm{H}}= $ 400 GeV (left) and in the ($ \tan\beta $, $ \cos(\beta-\alpha) $) parameter space at $ m_{\mathrm{A}}= $ 600 GeV and $ m_{\mathrm{H}}= $ 400 GeV (right). The enclosed regions are excluded. The dash-dotted lines bound the regions where the width of the A boson is up to 25% relative to $ m_{\mathrm{A}} $.

png pdf
Figure 7-a:
Median expected (dashed lines) and observed (filled contours) 95% CL exclusion regions in the type-II 2HDM ($ \tan\beta $, $ m_{\mathrm{A}} $) parameter space at $ m_{\mathrm{H}}= $ 400 GeV (left) and in the ($ \tan\beta $, $ \cos(\beta-\alpha) $) parameter space at $ m_{\mathrm{A}}= $ 600 GeV and $ m_{\mathrm{H}}= $ 400 GeV (right). The enclosed regions are excluded. The dash-dotted lines bound the regions where the width of the A boson is up to 25% relative to $ m_{\mathrm{A}} $.

png pdf
Figure 7-b:
Median expected (dashed lines) and observed (filled contours) 95% CL exclusion regions in the type-II 2HDM ($ \tan\beta $, $ m_{\mathrm{A}} $) parameter space at $ m_{\mathrm{H}}= $ 400 GeV (left) and in the ($ \tan\beta $, $ \cos(\beta-\alpha) $) parameter space at $ m_{\mathrm{A}}= $ 600 GeV and $ m_{\mathrm{H}}= $ 400 GeV (right). The enclosed regions are excluded. The dash-dotted lines bound the regions where the width of the A boson is up to 25% relative to $ m_{\mathrm{A}} $.
Summary
A direct search for heavy neutral Higgs bosons A and H in the $ {\mathrm{t}\overline{\mathrm{t}}} \mathrm{Z} $ final state has been conducted, utilising 138 fb$ ^{-1} $ of proton-proton collision data collected by the CMS experiment at a centre-of-mass energy of 13 TeV. The search targets the process $ \mathrm{A}\to\mathrm{Z}\mathrm{H} $ with $ \mathrm{H}\to{\mathrm{t}\overline{\mathrm{t}}} $. For the first time, the final state with decays of the Z boson to electrons or muons together with all-jet decays of the $ \mathrm{t} \overline{\mathrm{t}} $ system is considered. No signal is observed, and stringent upper limits are set on the cross section times branching ratio for A and H boson masses up to 2100 GeV and 1900 GeV, respectively, assuming narrow resonances. The results are further used to constrain the parameter space of two-Higgs-doublet models. Regions with A boson masses between 550 and up to 1500 GeV, and H boson masses between 350 and up to 700 GeV, are excluded at 95% confidence level, depending on the value of $ \tan\beta $. The results presented in this note do not confirm the excess reported by the ATLAS Collaboration in the region around $ (m_{\mathrm{A}},m_{\mathrm{H}})= $ (650, 450) GeV with a local significance of 2.85 SD [16]. The results complement and substantially extend the reach of previous searches, constraining parameter regions relevant in the context of baryogenesis.
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 Collaboration A detailed map of Higgs boson interactions by the ATLAS experiment ten years after the discovery Nature 607 (2022) 52 2207.00092
5 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
6 H. E. Haber and G. L. Kane The search for Supersymmetry: Probing physics beyond the standard model Phys. Rept. 117 (1985) 75
7 J. E. Kim Light pseudoscalars, particle physics and cosmology Phys. Rept. 150 (1987) 1
8 M. Bauer, U. Haisch, and F. Kahlhoefer Simplified dark matter models with two Higgs doublets: I. pseudoscalar mediators JHEP 05 (2017) 138 1701.07427
9 M. Trodden Electroweak baryogenesis: A brief review in 33rd Rencontres de Moriond: Electroweak Interactions and Unified Theories, 1998 hep-ph/9805252
10 T. Biekötter et al. The trap in the early universe: impact on the interplay between gravitational waves and LHC physics in the 2HDM JCAP 03 (2023) 031 2208.14466
11 T. D. Lee A theory of spontaneous T violation PRD 8 (1973) 1226
12 G. C. Branco et al. Theory and phenomenology of two-Higgs-doublet models Phys. Rept. 516 (2012) 1 1106.0034
13 U. Haisch and G. Polesello Searching for heavy Higgs bosons in the $ \mathrm{t}\overline{\mathrm{t}}\mathrm{Z} $ and $ \mathrm{t}\mathrm{b}\mathrm{W} $ final states JHEP 09 (2018) 151 1807.07734
14 G. C. Dorsch, S. J. Huber, and J. M. No A strong electroweak phase transition in the 2HDM after LHC8 JHEP 10 (2013) 029 1305.6610
15 P. Basler et al. Strong first order electroweak phase transition in the CP-conserving 2HDM revisited JHEP 02 (2017) 121 1612.04086
16 ATLAS Collaboration Search for a CP-odd Higgs boson decaying to a heavy CP-even Higgs boson and a Z boson in the $ \ell\ell{\mathrm{t}\overline{\mathrm{t}}} $ and $ \nu\overline{\nu}\mathrm{b}\overline{\mathrm{b}} $ final states using 140 fb$^{-1}$ of data collected with the ATLAS detector ATLAS conference note ATLAS-CONF-2023-034, 2023
link
17 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
18 CMS Collaboration CMS luminosity measurement for the 2017 data-taking period at $ \sqrt{s}= $ 13 TeV CMS Physics Analysis Summary
CMS-PAS-LUM-17-004
CMS-PAS-LUM-17-004
19 CMS Collaboration CMS luminosity measurement for the 2018 data-taking period at $ \sqrt{s}= $ 13 TeV CMS Physics Analysis Summary, 2018
CMS-PAS-LUM-18-002
CMS-PAS-LUM-18-002
20 CMS Collaboration The CMS experiment at the CERN LHC JINST 3 (2008) S08004
21 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
22 CMS Collaboration The CMS trigger system JINST 12 (2017) P01020 CMS-TRG-12-001
1609.02366
23 CMS Collaboration Particle-flow reconstruction and global event description with the CMS detector JINST 12 (2017) P10003 CMS-PRF-14-001
1706.04965
24 CMS Collaboration Technical proposal for the phase-2 upgrade of the Compact Muon Solenoid CMS Technical Proposal CERN-LHCC-2015-010, CMS-TDR-15-02, 2015
CDS
25 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
26 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
27 CMS Collaboration ECAL 2016 refined calibration and Run2 summary plots CMS Detector Performance Summary CMS-DP-2020-021, 2020
CDS
28 CMS Collaboration Pileup mitigation at CMS in 13 TeV data JINST 15 (2020) P09018 CMS-JME-18-001
2003.00503
29 M. Cacciari, G. P. Salam, and G. Soyez The anti-$ k_{\mathrm{T}} $ jet clustering algorithm JHEP 04 (2008) 063 0802.1189
30 M. Cacciari, G. P. Salam, and G. Soyez FastJet user manual EPJC 72 (2012) 1896 1111.6097
31 M. Cacciari, G. P. Salam, and G. Soyez The Catchment Area of Jets JHEP 04 (2008) 005 0802.1188
32 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
33 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
34 E. Bols et al. Jet flavour classification using DeepJet JINST 15 (2020) P12012 2008.10519
35 CMS Collaboration Performance summary of AK4 jet b tagging with data from proton-proton collisions at 13 TeV with the CMS detector CMS Detector Performance Summary CMS-DP-2023-005, 2023
CDS
36 GEANT4 Collaboration GEANT 4---a simulation toolkit NIM A 506 (2003) 250
37 NNPDF Collaboration Parton distributions from high-precision collider data EPJC 77 (2017) 663 1706.00428
38 T. Sjöstrand et al. An introduction to PYTHIA 8.2 Comput. Phys. Commun. 191 (2015) 159 1410.3012
39 CMS Collaboration Extraction and validation of a new set of CMS PYTHIA 8 tunes from underlying-event measurements EPJC 80 (2020) 4 CMS-GEN-17-001
1903.12179
40 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
41 C. Degrande Automatic evaluation of UV and R2 terms for beyond the standard model Lagrangians: a proof-of-principle Comput. Phys. Commun. 197 (2015) 239 1406.3030
42 G. C. Branco, L. Lavoura, and J. P. Silva CP Violation volume 103. Clarandon Press, Oxford, 1999
43 C. Degrande et al. UFO - the Universal FeynRules Output Comput. Phys. Commun. 183 (2012) 1201 1108.2040
44 P. Nason A new method for combining NLO QCD with shower Monte Carlo algorithms JHEP 11 (2004) 040 hep-ph/0409146
45 S. Frixione, P. Nason, and C. Oleari Matching NLO QCD computations with parton shower simulations: the POWHEG method JHEP 11 (2007) 070 0709.2092
46 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
47 T. Je \v z o and P. Nason On the treatment of resonances in next-to-leading order calculations matched to a parton shower JHEP 12 (2015) 065 1509.09071
48 M. Beneke, P. Falgari, S. Klein, and C. Schwinn Hadronic top-quark pair production with NNLL threshold resummation NPB 855 (2012) 695 1109.1536
49 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
50 P. B ä rnreuther, M. Czakon, and A. Mitov Percent-level-precision physics at the tevatron: next-to-next-to-leading order QCD corrections to $ \mathrm{q}\overline{\mathrm{q}}\to\mathrm{t}\overline{\mathrm{t}}\text{+X} $ PRL 109 (2012) 132001 1204.5201
51 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
52 M. Czakon and A. Mitov NNLO corrections to top pair production at hadron colliders: the quark-gluon reaction JHEP 01 (2013) 080 1210.6832
53 M. Czakon, P. Fiedler, and A. Mitov Total top-quark pair-production cross section at hadron colliders through $ o(\alpha\frac{4}{S}) $ PRL 110 (2013) 252004 1303.6254
54 M. Czakon and A. Mitov Top++: a program for the calculation of the top-pair cross-section at hadron colliders Comput. Phys. Commun. 185 (2014) 2930 1112.5675
55 CMS Collaboration Measurement of differential cross sections for the production of top quark pairs and of additional jets in lepton+jets events from pp collisions at $ \sqrt{s}= $ 13 TeV PRD 97 (2018) 112003 CMS-TOP-17-002
1803.08856
56 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
57 J. M. Lindert et al. Precise predictions for V+jets dark matter backgrounds EPJC 77 (2017) 829 1705.04664
58 R. Frederix and S. Frixione Merging meets matching in MC@NLO JHEP 12 (2012) 061 1209.6215
59 A. Kulesza et al. Associated production of a top quark pair with a heavy electroweak gauge boson at NLO+NNLL accuracy EPJC 79 (2019) 249 1812.08622
60 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
61 E. Re Single-top Wt-channel production matched with parton showers using the POWHEG method EPJC 71 (2011) 1547 1009.2450
62 M. Aliev et al. HATHOR: HAdronic Top and Heavy quarks crOss section calculatoR Comput. Phys. Commun. 182 (2011) 1034 1007.1327
63 P. Kant et al. HatHor for single top-quark production: Updated predictions and uncertainty estimates for single top-quark production in hadronic collisions Comput. Phys. Commun. 191 (2015) 74 1406.4403
64 N. Kidonakis Two-loop soft anomalous dimensions for single top quark associated production with $ \mathrm{W^-} $ or $ \mathrm{H}^{-} $ PRD 82 (2010) 054018 1005.4451
65 T. Gehrmann et al. $ \mathrm{W^+}\mathrm{W^-} $ production at hadron colliders in next to next to leading order QCD PRL 113 (2014) 212001 1408.5243
66 J. M. Campbell, R. K. Ellis, and C. Williams Vector boson pair production at the LHC JHEP 07 (2011) 018 1105.0020
67 R. J. Barlow and C. Beeston Fitting using finite Monte Carlo samples Comput. Phys. Commun. 77 (1993) 219
68 J. S. Conway Incorporating nuisance parameters in likelihoods for multisource spectra in Proceedings, workshop on statistical issues related to discovery claims in search experiments and unfolding, 2011
CERN-2011-006.115
1103.0354
69 T. Junk Confidence level computation for combining searches with small statistics NIM A 434 (1999) 435 hep-ex/9902006
70 A. L. Read Presentation of search results: The $ \text{CL}_\text{s} $ technique JPG 28 (2002) 2693
71 G. Cowan, K. Cranmer, E. Gross, and O. Vitells Asymptotic formulae for likelihood-based tests of new physics EPJC 71 (2011) 1554 1007.1727
72 ATLAS and CMS Collaborations, and the LHC Higgs Combination Group Procedure for the LHC Higgs boson search combination in summer 2011 ATL-PHYS-PUB-2011-011, CMS NOTE-2011/005, 2011
link
73 R. V. Harlander, S. Liebler, and H. Mantler SusHi Bento: Beyond NNLO and the heavy-top limit Comput. Phys. Commun. 212 (2017) 239 1605.03190
74 D. Eriksson, J. Rathsman, and O. Stal 2HDMC: Two-Higgs-Doublet Model Calculator physics and manual Comput. Phys. Commun. 181 (2010) 189 0902.0851
75 J. Bernon et al. Scrutinizing the alignment limit in two-Higgs-doublet models: m$ _h $=125\,GeV PRD 92 (2015) 075004 1507.00933
Compact Muon Solenoid
LHC, CERN