A novel experimental search channel for very light higgs bosons in the 2HDM type I

We present a reinterpretation study of existing results from the CMS Collaboration, speciﬁcally, searches for light Beyond the Standard Model (BSM) Higgs pairs produced in the chain decay pp → H SM → hh ( aa ) into a variety of ﬁnal states, in the context of the CP-conserving 2-Higgs Doublet Model (2HDM) Type-I. Through this, we test the Large Hadron Collider (LHC) sensitivity to a possible new signature, pp → H SM → Z A → Z Zh , with Z Z → j j μ + μ − and h → b ¯ b . We perform a systematic scan over the 2HDM Type-I parameter space, by taking into account all available theoretical and experimental constraints, in order to ﬁnd a region with a potentially visible signal. We investigate the signiﬁcance of it through a full Monte Carlo simulation down to the parametrised detector level. We show that such a signal is an alternative promising channel to standard four-body searches for light BSM Higgs bosons at the LHC already with an integrated luminosity of L = 300 fb − 1 . For a tenfold increase of the latter, discovery should be possible over most of the allowed parameter space.


Introduction
One of the main goals of the LHC machine is to investigate the individual properties (mass, width, spin, CP quantum numbers) of the Higgs boson as well as its interactions (with both matter and forces) and to look into evidence for new physics.These features have been probed by the ATLAS and CMS collaborations in different production and decay channels, using proton-proton ( pp) collision data collected at centreof-mass energies of 7, 8 [1] and 13 TeV [2,3].More recently, there have been observations involving the rare decay of the Higgs boson into a Z -boson and a photon (H → Z γ ) a e-mail: souad.semlali@soton.ac.uk (corresponding author) with a significance of 3.4σ [4], which could provide insights into the Higgs boson's coupling to both the electroweak gauge bosons and photons, contributing to our comprehension of the dynamics and interplays within the Higgs sector.Although the current measurements of the Higgs mass, spin, width and couplings to SM fermions and vector bosons [2,3] are all indeed in good agreement with the SM theoretical predictions, the small deviations in the SM-like Higgs couplings probed in various production modes for the five key decay channels H → γ γ, Z Z * , W W * , ττ and bb [2,3] still provide signs of possible potential BSM contributions to the total Higgs width and hints of new physics through the invisible and/or undetected decays.The indirect constraints from the current fit of couplings measurements and direct searches for H → inv (i.e., to'invisible' final states) performed by ATLAS and CMS collaborations have placed upper limits on the Branching Ratio (BR) of Higgs boson to invisible particles and undetected BSM particles at 95% C.L. (Confidence Level) [5][6][7][8].
The quest to uncover the trilinear couplings also holds a prominent position on the LHC's research agenda.Such interactions can be probed with sufficient luminosity, although, at present (i.e., at the end of Run 2), they are not determined yet.A measurement of this interaction is one of the highest priority goals during, possibly, Run 3 and, certainly, at the High Luminosity LHC (HL-LHC), both of which would, therefore, start shedding light on the nature of the Higgs boson and the shape of the Higgs potential, which in turn has implications for the vacuum metastability, the hierarchy problem as well as the strength of the Electro-Weak (EW) phase transition.However, probing Higgs's self-interactions, both trilinear and quartic couplings, in multi-Higgs production is experimentally very challenging.From the theory side on the other hand, many models with an extended scalar sector, like the 2-Higgs Doublet model (2HDM) [9][10][11][12][13], the Next-to-2HDM (N2HDM) [14][15][16][17] and a variety of both minimal and nonminimal Supersymmetric (SUSY) models [18][19][20][21][22][23][24] motivate additional features of new di-Higgs final states, as they all present with additional CP-even and/or -odd Higgs states, which can be accessible by the LHC experiments in a variety of signatures.Among the prominent processes are searches for exotic Higgs decays to a pair of light scalars or pseudoscalars, e.g.H → aa, hh, which then decay to SM particles.
This paper focuses on the popular 2HDM.After EW Symmetry Breaking (EWSB), the scalar sector of the 2HDM predicts five physical Higgs states, two CP-even Higgs bosons (h, H , with m h < m H ), one CP-odd one (a) and a pair of charged ones (H ± ).The rich (pseudo)scalar sector of the 2HDM and the different sets of Yukawa couplings that can be realised then offer a very interesting production and decay phenomenology of neutral and charged Higgs states at the LHC, even after scrutinising the 2HDM parameter space by considering different theoretical (vacuum stability, perturbativity, unitarity, etc.) and experimental (from SM-like Higgs data and null searches for companion states, flavour physics and low energy observables, etc.) constraints.Furthermore, the 2HDM is also attractive because one can impose a simple Z 2 discrete symmetry to the Yukawa sector in order to suppress Flavour Changing Neutral Currents (FCNCs) at tree level [25,26], which then forces one doublet to couple to a given type of fermions and leading as a result to four Yukawa interactions (termed, Type-I, Type-II, Type-X and Type-Y).In fact, in order to realise EWSB in such a way that the 2HDM is compliant with all experimental data, it is finally customary to allow for a soft breaking of this Z 2 symmetry.Herein, we will use the latter setup with a Type-I Yukawa structure.This particular scenario can accommodate both light neutral (10-100 GeV) and charged scalars (100-200 GeV), while one can not obtain such light Higgs states within the 2HDM Type-II or Type-Y configuration due ot the significant constraints arising from flavour physics, specifically from B → X s γ which enforces a lower bound of 800 GeV on the charged Higgs mass [27].
In the present study, we plan to take advantage of the direct access to some trilinear Higgs couplings that the LHC can access, entering multi-Higgs processes such as H → hh(aa) and H ± → W ± a, to test the sensitivity of the parameter space combination of 2HDM Type-I to light Higgs searches in cascade (or chain) decays.As the analysis progress, we aim to explore the scope of a new search for light Higgs bosons at the LHC Run 3 (with an integrated luminosity of 300 fb −1 ) as well as the HL-LHC (with an integrated luminosity of 3000 fb −1 ), on the basis of the knowledge acquired from the study of the aforementioned signatures.We focus mainly on the configuration of Type-I with inverted scenario, which in turn offers an alternative and new promising signal, in the form of the following cascade decays H → Z * a → Z * Z * h → bbμ − μ + j j.The main Higgs production process is via gluon fusion gg → H .

2HDM type-I
The 2HDM is one of the simplest well-motivated extensions of the SM.In this section, we briefly review the theoretical structure of this model.The scalar sector of the 2HDM consists of two complex SU (2) L doublets, Φ 1 and Φ 2 , with hypercharge Y = +1.The most general SU (2) L × U (1) Y invariant scalar potential can be written as follows [11][12][13]: Assuming CP-conservation in the 2HDM and following the hermiticity of the scalar potential, m 2 11 , m 2 22 , m 2 12 , λ 1,2,3,4,5,6 are real parameters.Invoking the described Z 2 symmetry, to avoid tree-level Higgs-mediated FCNCs at tree level, implies that λ 6 = λ 7 = 0. Also notice that the bilinear term proportional to m 2 12 breaks the Z 2 symmetry softly.Using the two minimisation conditions of the scalar potential and the combination , one can then trade the Lagrangian parameters of the 2HDM for a more convenient set of variables, , where α is the CP-even mixing angle, v 1 and v 2 are the Vaccum Expectations Values (VEVs) of the two Higgs doublets Φ 1 and Φ 2 , respectively.

Yukawa couplings
The general structure of the Yukawa Lagrangian when both Higgs fields couple to all fermions is given by [13]: where Q L and L L are the weak isospin quark and lepton doublets, u R and d R denote the right-handed quark singlets while Y u 1,2 , Y d 1,2 and Y l 1,2 are couplings matrices in flavour space.This form of Yukawa interaction gives rise to large FCNCs at tree level, which is strongly constrained by B-physics observables.Implementing Z 2 symmetry [25,26] allows only one doublet to couple to a given right-handed fermion field.
Depending on the Z 2 assignment, one can have the four types of models previously refered to as Type-I, Type-II, Type-X and Type-Y.In the mass-eigenstate basis, they can be unified and expressed as follows: where P L ,R = (1 ± γ 5 )/2 are the projection operators and V denotes the Cabibbo-Kobayashi-Maskawa (CKM) matrix.
Here, we focus only on Type-I, where only one doublet Φ 2 couples to all fermions, and thus the Higgs-fermion couplings are flavour diagonal in the fermion mass basis and depend only on the mixing angles, α and β, as follows: where we have used the short-hand notation c and s for cos and sin, respectively.

Theoretical and experimental constraints
We now describe briefly a set of, in turn, theoretical and experimental constraints that must be satisfied by the parameter space of the 2HDM.
The above constraints have been implemented in 2HDMC−1.8.0 [32].This public code is then used to explore the 2HDM parameter space and to compute the different Higgs BRs in each point of it.(2HDMC also provides an interface to HiggsBounds and HiggsSignals, see below.) Experimental observations impose the following constraints: • EW precision observables, i.e., the oblique parameters, S, T and U [33,34] are required to be within 95% C.L. of their experimental measurements, the current fit values (with the correlation parameters) are given by [35]: S = −0.02± 0.10, T = 0.03 ± 0.12, U = 0.01 ± 0.11, ρ ST = 0.92, ρ SU = −0.80,ρ T U = −0.93,χ ST, SU, T U < 5.99.
• Consistency with the Z width measurement Γ Z = 2.4952 ± 0.0023 GeV from LEP [36,37] is required.• Constraints from LHC, Tevatron and LEP searches which failed to find companion Higgs states are taken into account via HiggsBounds−5.10.0 [38], which allows to test the exclusion limits at 95% C.L. • The code HiggsSignals−2.6.2 [39] is used to check the signal strength measurements of the SM-like Higgs boson discovered at the LHC in 2012.• Constraints from B meson decays are enforced by Superiso-v1.4[40], using the current measurements: - • Constraints from recent searches for light pseudoscalar states in the mass range [15, 62.5] GeV, in protonproton collision at However, lately, additional constraints from such Higgs cascade decays have emerged, not included in the numerical tool, so we had to deal with these separately.For example, the CMS group has reported a search for H → aa → 4γ [49], using the data collected at √ s = 13 TeV, with an integrated luminosity of 132 fb −1 .Upper limits can then be set on BR(H → aa → 4γ ) at 95% C.L, since no significant deviation is observed. 1The ATLAS group [53] has also recently searched for the exotic decay of the Higgs boson into two light pseudoscalars in μ + μ − bb final state at √ s = 13 TeV with an integrated luminosity of 137 fb −1 , in the range of masses varying from 15 to 60 GeV.An upper limit is placed on BR(H → aa → μ + μ − bb) 2 at 95%C.L. Table 1 summarises several searches for exotic decays of the Higgs bosons in various final states, performed by the two collaborations ATLAS and CMS at Run 2, targeting a different ranges of masses.

Numerical analyses
The (pseudo)scalar sector of the 2HDM involves two CPeven Higgs bosons, h and H .One of these scalars can be identified as the 125 GeV state observed at the LHC.As mentioned, in this analysis, we will assume that the heaviest Higgs state H is the SM-like one with a mass of 125 GeV and that h and a are lighter than H .We the perform a scan over the following ranges, with Assuming m H = 125 GeV and m h,a < 90 GeV, the decay channels H → hh, aa, a Z * could be open, leading to invisible or undetected SMlike Higgs decays that are restricted by the current precision measurements of Higgs couplings.CMS performed a combination of searches, using data collected at √ s = 7, 8, 13 TeV [7], for Higgs bosons decaying into invisible parti-Footnote 1 continued procedure which was validated against the case of BR(H → hh(aa) 2 Corresponding search data and exclusion limits are available at the HEPData database.cles, which targets the following production channels: Vector Boson Fusion (VBF), Higgs-Strahlung (HS) and gluongluon Fusion (ggF) (allowing for initial state radiation).The combination of all the searches, assuming these SM-like production modes, yields an observed (expected) upper limit on BR(H → inv) of 0.19 (0.15) at 95% C.L. The ATLAS group reported a direct search for Higgs bosons produced via VBF with subsequent invisible decays, for 139 fb −1 of pp collision data at √ s = 13 TeV [8].An observed (expected) upper limit of 0.145 (0.103) is placed on BR(H → inv) at 95% C.L., as a function of the assumed production cross sections.As for now, both ATLAS and CMS have placed, respectively, an upper bound of 0.10 [6] and 0.15 [5] on BR(H → invisible) at 95% CL.In our analysis, we will assume that BR(H → inv) designates the sum of the following decay rates, BR(H → hh), BR(H → aa) and BR(H → a Z * ).
After performing a random scan over 2HDM Type-I parameters, we show in Fig. 1 the allowed regions by the- Each (coloured) point in the (m h , m a ) plane implies that there is a combination of the scanned model parameters, which obeys the aforementioned theoretical constraints and evades the current experimental limits in all searched modes, whereas the white space corresponds to the case where any possible mass combination is forbidden by an observed signature(s) in one or more existing experimental searches.
Within this region, the most sensitive channels for the model parameter points, as determined by HiggsBounds, are shown by coloured dots.Note that each coloured point in the parameter space of the 2HDM is allowed by current experimental searches.Obviously, there are two distinct regions in the figure.The one in the top left corner corresponds to low masses of h (m h < m H /2), and high masses of a (m a > m H /2), while the second one corresponds to the m a, h > m H /2 scenario.It is interesting to note that there are no acceptable points when 40 GeV < m h < m H /2 and m a > m H /2. This is due to the fact that this parameter combination is excluded by LEP searches for e + e − → ah → bbbb [54] and an ATLAS search for events with at least 3γ in Finally, it is noteworthy that the most sensitives searches for the region with low m h and high m a are the LEP searches for processes such as e + e − → ah → bbbb and e + e − → (h)Z → (bb)Z [54].Therefore, an update from the LHC during Run 3 is unlikely to rule out this mass combination over the plane (m h , m a ) of the 2HDM Type-I.We will be focusing on this region in the second part of our study.
We now turn to the reinterpretation of exotic Higgs decay searches, i.e., H → aa in τ + τ − bb, μ + μ − bb and μ + μ − τ + τ − final states in the framework of the 2HDM Type-I, while taking advantage of the parameter space discussed above.The recasting of τ + τ − bb, μ + μ − bb and μ + μ − τ + τ − searches for H → hh is also possible since these processes share similar kinematics (in the same spirit as in Ref. [56]).It is relevant to note that the constraints from the search for light pseudoscalars in the τ + τ − bb final state are much stronger than the ones from μ + μ − bb and μ + μ − τ + τ − searches.CMS has set an upper limit, between 1.7% and 7.6%, on BR(H → aa → τ + τ − bb) at 95% C.L. [51], assuming the SM production of primary Higgs boson.We show in Fig. 2 the outcome from reinterpreting the H → aa(hh) → τ + τ − bb search [51] in the 2HDM Type-I.The yellow and green bands represent the uncertainties at ±1σ and ±2σ associated with the expected exclusion limits.Grey points satisfy theoretical constraints described in Sect.2.2, whereas red points are excluded by null searches (i.e., by HiggsBounds [38]).The blue points satisfy both theoretical and experimental constraints.The area of sensitivity to H → aa(hh) → τ + τ − bb is already excluded by existing experimental searches (red points).In this connection, the BR of Higgs SM-like Higgs state decaying into hh and/or aa is very restricted and cannot exceed 9% at 95% C.L., again, in the 2HDM Type-I.
One can draw a similar conclusion form reinterpreting H → hh(aa) → μ + μ − bb [51] in our reference framework.Figure 3 shows that the parameter space with sensitivity to this search is excluded.One should keep in mind that the μ + μ − bb final state is well-balanced between large BR(h/a → bb) and a clean di-muon resonance that is easy to trigger on.This exercise emphasises that the 2HDM Type-I may not be a good framework for reinterpreting searches for exotic Higgs decays into light pseudoscalar in "traditional" final states such as μ + μ − bb, τ + τ − bb and μ + μ − τ + τ − .
We also address here light charged Higgs decay in the mass ranges where m H ± < m t − m b and m h,a < 90 GeV.In this configuration, the charged Higgs state can be be produced from top quark decays, i.e., t → bH + , followed by its bosonic decays to H ± → W ± h(a), instead of the standard fermionic decay modes like τ ν and cs.Many studies motivated these channels as alternative modes to search for light charged Higgs bosons that could dominate over the conventional fermionic channels, because of large BRs when they are kinematically allowed, in models such as our 2HDM Type-I [57][58][59][60].ATLAS [61] and CMS [62] have considered the ranges m a ∈ [15,75] GeV and m H ± < m t −m b to search for light charged Higgs bosons in pp → tt → bbH + W − with H + → W + a and a → μ + μ − at √ s = 13 TeV, since the μ + μ − finale state provides the aforementioned experimental advantages, which offset the suppressed rate of BR(a → μ + μ − ).Previously, both CDF and the LEP collaborations have searched for H ± → W ± a with a → bb [63], a → τ + τ − [64] and a → bb [65].In addition, LEP experiments [66] have set a lower bound on the charged Higgs boson mass of m H ± > 72.5 GeV in the 2HDM Type-I for m a > 12 GeV at 95% C.L.
Figure 4 shows the CMS observed and expected exclusion limits on the product of the BRs of t → bH ± , H ± → W ± a and a → μ + μ − [62] as a function of m a predicted by the 2HDM Type-I, with respect to several theoretical and experimental constraints.We adopt here m H ± = m a + 85 GeV [62], which enables us to consider H ± → W ±( * ) a, with W ±( * ) being on/off shell, by randomly sampling values of the charged Higgs mass between 100 GeV and 160 GeV (see Eq. ( 4)).A noteworthy observation is that the 2HDM Type-I offers sufficient sensitivity, when the prediction of We move now to discuss a new analysis, where we deploy the parameter space of the 2HDM Type-I following the outcomes of reinterpreting previous searches for light Higgs bosons, pp → H SM → hh(aa), in different final states, in order to search for a new signature.
Figure 5 shows the result of performing a scan over the parameter space of 2HDM Type I, wherein (recall) the heav-  Instead, the predominant channel in this region is a → Z * h with Z being off-shell, boasting a branching ratio that can reach 90%.Furthermore, the decay of H SM → a Z * is suppressed in the region with large m h and m a .In fact, H SM tends to decay to a Z * in the specific part of the region where the observation of a → Z * h is not possible.We show in Fig. 6 the gg → H → a Z * → h Z * Z * cross section, where h → bb.The process could yield a cross section of 0.006 pb.In the right panel of Fig. 6 we show the BR of h → bb in this region of the 2HDM Type-I parameter space.Obviously, the decay width of h is dominated by the decay mode h → bb.Thus, in what follows, we focus on the case where h decays to bb and Z ( * ) Z ( * ) → μ + μ − j j.Such a scenario could be an alternative channel to search for light Higgs bosons at Run 3 and the HL-LHC.

Signal vs. background analysis
We describe here the toolbox used to generate and analyse MC events.MadGraph-v.9.2.5 [70] is used to generate parton level configurations of both signal and background processes. 4The events are passed then to PYTHIA8 [72] to simulate parton showering, hadronisation and decays.Finally, we use Delphes−3.5.0 [73] with the standard CMS card 5 to perform detector simulation.We resort to MadAnalysis [74] to apply cuts and to conduct the analysis.
The major background processes are top pair production in association with 2 Initial State Radiation (ISR) jets, 6Z Z production with additional bb quarks, Z W + bb → μ + μ − j jbb and Drell-Yan plus jets (DY+jets).We show in Table 3 the corresponding cross sections at √ s = 13 TeV for the LHC energy.We have generated MC samples of (10 6 ) events.Unsurprisingly, the irreducible backgrounds pp → Z ( * ) Z ( * ) bb → bb j jμ + μ − (from both QCD and EW interactions) and pp → Z W bb → bb j jμ + μ − are negligible whereas pp → ggtt → ggμ + μ − j jbbν μ ν μ and DY+jets are large.
We considered a few BPs for the signal given by gg → H → a Z * → h Z * Z * → μ + μ − j j bb to perform the MC simulation.The input parameters of each BP are given in 4 Background and signal events are generated at LO. Higher order corrections are quantified through K -factors.The NLO QCD correction to top pair production in association with 2 jets computed at the LHC is about −27% [71], which we adopt here.The NLO corrections to gg → H are very large , about a facor of 2, due to the contributions from gg pairs to QCD radiation, whereas K NNLO/NLO ≈ is much smaller than K NLO/LO , signifying a convergence of the QCD expansion, so we renormalise the signal to the NNLO rates through the K -factor K NNLO/LO = σ NNLO /σ LO ∼ 2.6 − 2.7. 5 It adopts the anti-k T algorithm to cluster final particles into jets, with jet parameter ΔR = 0.5 and p min T, j = 20 GeV (for both light and b jets) Table 4.Note that the light Higgs width, Γ (h), is not small enough to lead to a large lifetime and hence, long-lived particles producing displaced vertices inside the detector.The proper decay length cτ 0 is in fact only a tiny fraction of micrometers. 7The different kinematic distributions at parton level in Fig. 7 show that the requirement of central pseudorapidity of the muons is generally satisfied however the p T of these can be rather small.To address this, one can invoke the di-muon scouting trigger, 8 which involves lowering the transverse momentum ( p μ T ) of muons to 4.5 GeV [77,78].Note that the p T threshold of such trigger is reduced to 3 GeV [79] in Run 3. Figure 8 shows the invariant mass distributions of the two b-jets, m bb , and that of the full final state, m j jμ + μ − bb , for the signal and the irreducible background processes at parton level, noting that m bb is close to light Higgs mass m h and m j jμ + μ − bb is close to SM-like Higgs mass m H (for the signal, unlike the irreducible backgrounds).We will clearly leverage these underlying partonic shapes in our detector level 9 analysis, to which we proceed next, in the presence of the following sequence of acceptance cuts: Additionally, we only consider events with two oppositely charged muons, two b-jets and two light jets in the final state.Events are further preselected after requiring the invariant mass of the two b-jets (m b b), the two light jets (m j j ) and the di-muon system 10 (m μμ ) to lie within the following mass ranges, as shown in Fig. 9, 5 GeV < m b b < 40 GeV, 10 GeV < m μμ < 50 GeV, 10 GeV < m j j < 50 GeV.
Figure 10 displays the distributions of the missing transverse energy (E / T ) and the highest p T 's of b-jets, light-jets and muons for signal and background processes at detector level.(As mentioned previously, the irreducible backgrounds stemming from Z Zbb and Z W b b processes are negligible, so we have not emulated these at detector level.) 7The proper decay length cτ 0 falls within the range from 0.06μm to 0.19μm, where τ 0 is the light Higgs lifetime at rest [75]. 8The CMS standard dimuon triggers require the highest p T muon to possess a transverse momentum within the range of 12-17 GeV, while the second-highest p T muon reconstructed should exhibit a p T between 5-8 GeV [76].This consequently leads to a significant loss of signal acceptance in searches for dimuon masses below 40 GeV [77,78]. 9The minimum p T requirements for muons has been lowered to 5 GeV to cover the entire signal region effectively.This modification is implemented in the Delphes card to simulate the CMS detector. 10Both μμ and j j originate from the decays of off-shell Z bosons (Z * Z * → μ + μ − j j).To enhance the signals and suppress the backgrounds arising from ggtt and DY+jets, we have adopted several kinematic cuts, which choice is based on comparing different distributions of the signal and background processes at the detector level.Specifically, this has been done through 2D distributions correlating the missing transverse momentum to a series of kinematic variables pertaining to some of the visible objects in the final state as illustrated in Figs. 12 and 13 for signal and background processes mentioned above.One can read that the signal and backgrounds distributions are anticorrelated.In fact, forcing the missing transverse energy to be below 30 GeV will strongly favour the signal over the backgrounds, specifically, ggt t. 11 Additionally, selecting events with p j T < 75 GeV and p μ T < 40 GeV would enhance the 11 Only 10% of t tgg events have survived after applying the missing transverse energy cut, as indicated in Table 5.
signal significance and suppress both background processes.Through similar reasoning, we require the invariant mass of the system to satisfy m H < 180 GeV.We show in Table 5 the event rates of the backgrounds after applying the cutflow discussed above.Notably, DY+jets events are significantly reduced after considering a final state with exactly 2 b-jets, 2 light jets and 2 muons, while adhering to the invariant masses requirements. 12Both backgrounds processes are completely removed after applying the kinematic cuts and restricting the invariant mass of the system to be below 180 GeV to further differentiate the signal from the remaining background processes.
We have then computed the significance (for √ s = 13 TeV and L = 300 fb −1 ), defined as 13where S(B) is the signal(background) yield after the discussed cutflow, for not only our three initial BPs (whose Σ rates are 3.05, 3.18 and 3.42 for BP1, BP2 and BP3, respectively), but also those appearing in Table 6.We have done so in order to be able to map the 2HDM Type-I parameter space in detail, so as to acquire a sense of the true portion of it that can be tested by forthcoming experiments.Note that   we have kept the same cutflow already illustrated for all such new BPs too.Also, it is at this stage that we take into account the aforementioned QCD K -factor for the signal.Many of the latter can have a significance larger than 3 and up to nearly 4, for Run 3 energy and luminosity.To observe their distribution over the (m h , m a ) plane, we have finally produced Fig. 14, indeed, assuming √ s = 13 TeV and L = 300 fb −1 , where both significance (Σ) and efficiency ( ) are mapped.Hence, at Run 3, we can conclude that a substantial portion of the 2HDM Type-I parameter space can offer some evidence of the signal we have pursued.Furthermore, we notice that a larger efficiency can be obtained for small m a : this is because the loss of efficiency with b-tagging is over-compensated by a simultaneous higher efficiency for both j-and μ-tagging.Needless to say, at the HL-LHC, where L = 3000 fb −1 , most of the sampled parameter space of the 2HDM Type-I would be discoverable.

Conclusions
In this paper, we have shown the outcome of performing some recasting over the parameter space of the 2HDM Type I, wherein the heaviest CP-even Higgs state H is identified with the discovered SM-like one, H SM , while h and a are lighter.After considering the available experimental data from searches for exotic Higgs decay into two light (pseudo)scalars, we have found that the corresponding parameter space for which there is sensitivity via H SM → hh(aa) → τ + τ − bb at Run 2 is already excluded by existing constraints from BSM Higgs searches.Furthermore, we have shown that there are regions of the 2HDM Type-I parameter space compliant with theoretical and experimental constraints yielding substantial BR(H ± → W ± a) and BR(H → Z * Z * h).The large size of the former has been exploited in other literature.Here, concerning the lat-  ter, we have made the case for looking at the process pp → H SM → Z * A → Z * Z * h via Z Z → μ + μ − j j and h → bb decays, specifically, in the region with large m a and small m h .After performing a full MC analysis down to detector level, we have proven that the overwhelming backgrounds arising from both top-quark pair production in association with 2 ISR jets and DY+jets can be suppressed after applying efficient kinematics cuts, leading to a large significance of this hitherto unexplored light Higgs signature already at Run 3 of the LHC, where evidence of it can be seen, further affording one with clear discovery potential at the HL-LHC.

Fig. 1
Fig. 1 Allowed parameter space in the 2HDM Type-I at 95% C.L. The solid red and green regions are excluded by the Z width constraint and the LEP search for h → aa [54], respectively.The most sensitive searches for the relevant (m h , m a ) regions are shown by coloured dots.Each coloured point in the (m h , m a ) plane satisfies all theoretical requirements and up-to-date experimental constraints.In the white space, any combination of masses is ruled out by current experimental searches

Fig. 2
Fig. 2 Observed and expected upper limits on B(H → aa(hh) → τ + τ − bb) [51] at 95% C.L. in the 2HDM Type-I.Grey points are allowed by theoretical constraints.As stated above, red colour indicates the com-

Fig. 12
Fig. 12 Correlation between p j T and E / T for signal (BP3) (red), ggt t (blue) and DY+jets (green) at detector level

Fig. 13
Fig. 13 Correlation between p μ T and E / T for signal (BP3) (red colour), ggt t (blue colour) and DY+jets (green colour) at detector level

Fig. 14
Fig. 14 Significance (left) and efficiency (right) of each BP produced in our analysis over the (m h , m a ) projection of the 2HDM Type-I parameter space, after the full cutflow described in the text

Table 4
Selected BPs with parton level cross section and other observables at LO. (All masses and widths are in GeV, with m H = 125 GeV.)

Table 6
Extended list of BPs used in the MC simulation for the 2HDM Type-I parameter scan, highlighting the h and a masses as well as the signal LO cross section and event rate after the full cutflow, together with its significance Σ and efficiency .Recall that NNLO QCD K -factor has been used for Higgs production.Here, √ s = 13 TeV and L = 300 fb −1