A Guidebook to Hunting Charged Higgs Bosons at the LHC

We perform a comprehensive global analysis in the Minimal Supersymmetric Standard Model (MSSM) as well as in the 2-Higgs Doublet Model (2HDM) of the production and decay mechanisms of charged Higgs bosons (H±) at the Large Hadron Collider (LHC). We start from accounting for the most recent experimental results (SM-like Higgs boson signal strengths and search limits for new Higgs boson states obtained at Run-1 and -2 of the LHC and previous colliders), from (both direct and indirect) searches for supersymmetric particles as well as from flavor observables (from both e+e− factories and hadron colliders). We then present precise predictions for H± cross sections and decay rates in different reference scenarios of the two aforementioned models in terms of the parameter space currently available, specifically, mapped over the customary (mA,H±,tanβ) planes. These include the mhmod+ and hMSSM configurations of the MSSM and the 2HDM Type-I, -II, -X, and -Y for which we also enforce theoretical constraints, such as vacuum stability, perturbativity, and unitarity. We also define specific Benchmark Points (BPs) which are always close to (or coinciding with) the best fits of the theoretical scenarios to experimental data. We finally briefly discuss the ensuing phenomenology for the purpose of aiding future searches for such charged Higgs boson states.


INTRODUCTION
The Higgs boson discovery of 2012 [1][2][3][4] at the CERN Large Hadron Collider (LHC) has led to the confirmation of the Standard Model (SM) as the proper theory of the Electro-Weak (EW) scale. However, there is much evidence that the SM is not appropriate at all scales, rather it should be viewed as an effective low-energy realization of a more complete and fundamental theory on setting beyond the EW regime. Among the many proposals for the latter, one can list theories with some new symmetry, e.g., Supersymmetry (SUSY), or an enlarged particle content (e.g., in the Higgs sector), or both. Following the aforementioned discovery, no new particle has however been seen at the LHC, implying that new physics at the EW scale should be weakly interacting or that strongly interacting particles, if present, should lead to signatures involving soft decay products or in channels with overwhelming (ir)reducible backgrounds. We shall adopt here the first assumption.
Many SM extensions possess in their spectra additional neutral and/or charged Higgs states. Amongst these, SUSY [5] is indeed considered the most appealing one as it addresses several shortcomings of the SM, including the problem of the large hierarchy between the EW and Planck scales as well as the dark matter puzzle. While the search for SUSY was unsuccessful during the first LHC run, the increase in the Center-of-Mass (CM) energy of the machine from 8 to 13 TeV plus the additional luminosity of the second run are improving greatly the sensitivity to the new superparticles which are predicted. While the jury is still out on these, we remind here the reader that SUSY also requires at least two Higgs doublets for a successful EW Symmetry Breaking (EWSB) pattern. For exactly two such fields, yielding the so-called Minimal Supersymmetric Standard Model (MSSM), also having the same gauge group structure of the SM, one obtains four physical Higgs particles, in addition to the discovered SM-like one with the observed mass of 125 GeV. In fact, the same Higgs mass spectrum also belongs to a generic 2-Higgs Doublet Model (2HDM), i.e., one not originating from SUSY. In neither case, though, there exists a precise prediction of the typical masses of the new Higgs states, though we know already that the MSSM allows for one to be lighter than the 125 GeV state in a very small region of parameter space [6], whereas the 2HDM generally does so over a significantly larger expanse of it [7,8]. Either way, the presence of extra physical Higgs boson states alongside the SMlike one is thus one of the characteristics of Beyond the SM (BSM) physics, whether within SUSY or otherwise. Hence, looking for these additional states in various production and decay channels over a wide range of kinematic regimes is an important part of the physics programme of the multi-purpose LHC experiments, ATLAS, and CMS. Specifically, the discovery of a (singly) charged Higgs boson would point to a likely additional Higgs doublet (or a Higgs field with higher representation, such as triplet). Hence, we concentrate on this Higgs state here.
The two Higgs doublet fields pertaining to the MSSM are required to break the EW symmetry and to generate the isospinup and -down type fermions as well as the W ± and Z boson masses [9][10][11]. The Higgs spectrum herein is given by the following states: two charged H ± 's, a CP-odd A and two CP-even Higgses h and H, with m h < m H (conventionally, wherein h is the SM-like Higgs state). The tree-level phenomenology of the Higgs sector of the MSSM is described entirely by two input parameters, one Higgs mass (that can be taken to be that of the CP-odd Higgs state, m A ) and the ratio tan β of the Vacuum Expectation Values (VEVs) of the two Higgs doublet fields. Note that one of the most powerful prediction of SUSY is the existence of a light Higgs boson that could be produced at colliders. In the MSSM, at the tree-level, the light CP-even h is predicted to be lighter than the Z boson. However, it is well-known that loop effects could raise the h mass upper bound to 135 GeV for a large soft breaking trilinear parameter, A t , and/or a heavy scalar top [12][13][14][15][16][17][18][19][20][21][22]. After the Higgs boson discovery at the LHC, MSSM benchmark scenarios have been refined to match the experimental data and to reveal characteristic features of certain regions of the parameter space [23][24][25]. Of the many MSSM frameworks presented in literature, we consider in this work the so-called m mod+ h [24] and hMSSM [26][27][28][29][30] ones, which will be described in the coming section. As for the 2HDM, one ought to specify the Yukawa sector, in order to proceed to study phenomenologically its manifestations. While SUSY enforces this in the form of a so-called Type-II, this is only one of four Ultra-Violet (UV) complete realizations of a generic 2HDM, the others been termed Type-I, -X, and -Y. The difference between these four scenarios is the way the fermionic masses are generated. We define as Type-I the model where only  one doublet couples to all fermions, Type-II is the scenario where  one doublet couples to up-type quarks and the other to downtype quarks and leptons, the Type-X is the model where one doublet couples to all quarks and the other to all leptons while a Type-Y is built such that one doublet couples to up-type quarks and to leptons and the other to down-type quarks. In all such cases, the number of free parameters at tree-level is seven to start with, hence it becomes more cumbersome than in SUSY to map experimental results onto theoretical constraints. Yet, in virtue of the fact that a 2HDM is the simplest realization of a BSM scenario based solely on doublet Higgs fields, its study is vigorously being pursued experimentally.
So far, the non-observation of any Higgs signal events in direct searches above and beyond those of the SM-like Higgs state constrains the parameter space of the underlying physics model. Specifically, in the case of the H ± boson, wherein the relevant phenomenological parameters are m H ± and tan β in whatever scenario, one can pursue the study of its production and decay modes in a model independent way, which results can a posteriori be translated to exclude the relevant parameter space in a given scenario (whether it be the MSSM, 2HDM, or something else). This recasting is conveniently done on the (m A , tan β) and (m H ± , tan β) planes for the MSSM and 2HDM, respectively, so that we will map our findings in the same way.
At hadron colliders, there exists many production modes for charged Higgs bosons which are rather similar in the MSSM and 2HDM. For a light charged Higgs, i.e., with mass m H ± +m b < m t , its production comes mainly from top decay. At the LHC, the production of top quark pairs proceeds via Quantum Chromo-Dynamics (QCD) interactions and, when kinematically allowed, one top could decay into a charged Higgs state and a bottom quark in a competition with the SM decay into a W ± boson and again a bottom quark. Therefore, the complete H ± production mechanism qq, gg → tt → tbH − provides the main source of light charged Higgs bosons at the LHC and offers a much more copious signature than any other form of direct production. After crossing the top-bottom threshold, i.e., when m H ± + m b > m t , a charged Higgs (pseudo)scalar can be produced through the process gb → tH − [31][32][33][34]. In fact, these two mechanisms can be simultaneously captured via the process gg → tbH − [35,36], which again makes it clear that one should expect large H ± cross sections induced by QCD interactions also in the heavy H ± mass range 1 .
In the MSSM, and also in a variety of 2HDM Types, light charged Higgs bosons would decay almost exclusively into a (hadronic or leptonic) τ lepton and its associated neutrino for tan β > ∼ 1. When the top-bottom channel is kinematically open, then H + → tb would compete with H ± → hW ± , HW ± , AW ± decays as well as various SUSY channels in the MSSM. In the latter, H + → tb → bbW + is the dominant channel and the bosonic decays H ± → hW ± , HW ± , AW ± (also yielding bbW + are explored in the mass range m H ± = 200-1,000 GeV, in the latter mode using multi-jet final states with one electron or muon from the top quark decay. No significant excess above the backgroundonly hypothesis has been observed and upper limits are set on the pp → tbH ± production cross section times BR(H ± → tb). Several interpretations of these limits have eventually been given in benchmark scenarios of the MSSM, including those mentioned above. Note that current ATLAS and CMS bounds are significantly weakened in the 2HDM once the exotic decay channels into a lighter neutral Higgs, e.g., H ± → hW ± or H ± → AW ± , are open. This scenario could also happen in the MSSM if one of the SUSY decay channels of charged Higgs bosons are open (such as into chargino-neutralino pairs). In the 2HDM, the possibility of producing a light charged Higgs boson from top decay with a subsequent step H ± → hW ± or H ± → AW ± was studied in Arhrib et al. [45] and it was shown that it can lead to sizable cross sections at low tan β. We stress here that there exist several recent analyses dedicated to 2HDM phenomenology [34,46,47] that we consulted. However, unlike Akeroyd et al. [34], Arbey et al. [46], and Bernon et al. [47] only concentrates on neutral Higgs phenomenology and discuss the charged Higgs contribution only to flavor physics observables without singling out the relevant charged Higgs production and decay channels at the LHC, which is indeed one of the aims of this analysis.
In this paper, we analyze the allowed σ (pp → tbH + + c.c.) × BR(H ± → anything) rates by taking into account both theoretical and experiments constraints on the underlying BSM model, the latter including the latest ATLAS and CMS results for SM-Higgs (h) and other Higgs (H, A, H + ) searches with the full set of 36.5 fb −1 data collected in the second LHC phase. We will then interpret these results under the proposed scenarios to quantify the magnitude of the available parameter space to be covered by future LHC analyses. In doing so, we will extract several Benchmarks Points (BPs) that could lead to detectable signals, all of which are consistent with the best fit regions in both the MSSM and 2HDM.
The paper is organized as follows. In the second section we review the MSSM and introduce the benchmark scenarios that we will discuss. The 2HDM, with its parameterizations and Yukawa textures, is described in the third section. The fourth section is devoted to a discussion of the theoretical and experimental constraints used in our study. Results and discussions for the MSSM and 2HDM are presented in the fifth section and we finish with our conclusions.

THE MSSM
In the MSSM, due to the holomorphy of the superpotential, one introduce two Higgs doublets 1,2 in order to give masses to up-type quarks as well as down-type quarks and leptons, respectively. Both Higgs fields acquire VEVs, denoted by v 1,2 . After EWSB takes place, the spectrum of the model contains the aforementioned Higgs states: h, H, A and H ± . The MSSM Higgs sector is parameterized at tree-level by tan β = v 2 /v 1 , and e.g., the CP-odd mass m A . One of the interesting features of the MSSM is the prediction, at the tree-level, of a light CP-even Higgs h with a mass m h ≤ m Z . However, such tree-level prediction is strongly modified by radiative corrections at one-and two-loop level [12][13][14][15][16][17][18][19][20][21]. It has been shown in Degrassi et al. [22], on the one hand, that the loop effects can make the light CP-even mass m h reach a value of 135 GeV and, on the other hand, that the theoretical uncertainties due to unknown high order effects can be of the order of 3 GeV. In fact, these large loop effects are welcome in order to shift the light CP-even Higgs mass to the measured experimental value m h ≈ 125 GeV. Note also that the loop effects will modify not only the tree-level Higgs mass relations but also the Higgs self-couplings and the Higgs coupling to SUSY particles. Therefore, beside the tree-level parameters tan β and m A , the top quark mass and the associated squark masses and their soft SUSY breaking parameters enter through radiative corrections [12-21, 48, 49]. In fact, when trying to push the light CP-even mass from m Z to 125 GeV through loop effects, one needs to introduce a large SUSY scale with large soft trilinear parameter A t . Such a large SUSY scale puts automatically the SUSY spectrum at the TeV scale, which is consistent with negative searches for SUSY particles at the LHC.
To compute the masses and couplings of Higgs bosons in a given point of the MSSM parameter space we use the code FeynHiggs [50,51] for the m mod+ h scenario and the program HDECAY for the hMSSM case [52,53]. Both codes include the full one-loop and a large part of the dominant two-loop corrections to the neutral Higgs masses. Since the theoretical uncertainty on the Higgs mass calculation in the FeynHiggs code has been estimated to be of the order of 3 GeV, we consider as phenomenologically acceptable the points in the MSSM parameter space where FeynHiggs predicts the existence of a scalar state with mass between 122.5 and 128.5 GeV and with approximately SM-like couplings to gauge bosons and fermions. In addition to the tree-level scalar potential parameters, tan β and m A , when taking into account high order corrections, the MSSM parameters most relevant to the prediction of the masses and production cross sections of the Higgs bosons are: the soft SUSY-breaking masses for the stop and sbottom squarks, which, for simplicity, we assume all equal to a common mass parameter M SUSY , the soft SUSY-breaking gluino mass mg, the soft SUSY-breaking Higgs-squark-squark couplings A t and A b , the superpotential Higgs(ino)-mass parameter µ and the leftright mixing terms in the stop and sbottom mass matrices (divided by m t and m b ) respectively. Since the (approximate) two-loop calculation of the Higgs masses implemented in FeynHiggs and the Next-to-Leading Order (NLO) calculation of QCD corrections to the production cross section implemented in SusHi [54,55] employ the same renormalization (on-shell) scheme, the input values of the soft SUSY-breaking parameters can be passed seamlessly from the Higgs mass to the cross section calculations.
In the light of the latest LHC data on the discovered Higgslike boson, and given the fact that the MSSM contains many independent parameters which makes it a fastidious task to perform a full scan, there have been many studies which lead to several benchmarks that could fit the observed Higgs boson as well as be tested at the future LHC with higher luminosity [24,[26][27][28][29]. As intimated already, in this study, we will concentrate on two of these benchmark scenarios: the m mod+ h and hMSSM ones, which we will describe hereafter.

Scenario
The m mod+ h scenario is a modification of the time-honored m max h scenario (also called maximal mixing scenario), which was originally defined to give conservative exclusion bounds on tan β in the context of Higgs boson searches at LEP [56]. The m max h scenario was introduced in order to maximize the value of m h by incorporating large radiative correction effects for a large m A >> m Z mass, fixed value of tan β > 8 and large SUSY scale of the order 1 TeV. However, this scenario predicts m h to be much higher than the observed Higgs boson mass, due to the large mixing in the scalar top sector.
Hence, the maximal mixing scenario has been modified, by reducing the amount of scalar top mixing, such that the mass of the lightest Higgs state, m h , is compatible with the mass of the observed Higgs boson within ±3 GeV in a large fraction of the considered parameter space. In fact, modifications of the m max h scenario can be done in two ways depending on the sign of (A t − µ cot β)/M SUSY , leading to an m mod− h and m mod+ h [24]. It has been demonstrated in Carena et al. [24] that when m mod+ h is confronted with LHC data, there is a substantial region in the plane (m A , tan β) with tan β > 7 for which the light CP-even Higgs mass is in a good agreement with the measured one at the LHC, hence our choice of this scenario.
The SUSY input parameters in this scenario are fixed as 2 where M SUSY is the aforementioned SUSY mass scale.

The hMSSM Scenario
In the previous scenario, one need to input tan β, m A and also the other SUSY parameters to get the Higgs and SUSY (mass and coupling) spectrum. Taking into account the theoretical uncertainty of the order 3 GeV on the Higgs mass, which could originate from unknown high order loop effects, a light CP-even Higgs boson with a mass in the range [122, 128] GeV would be an MSSM candidate for the observed Higgs-like particle. However, 2 Notice that this m mod+ h configuration is compliant with the theoretical and experimental constraints discussed below, including Dark Matter (DM) ones. So is the case for the hMSSM configuration below. plenty of points on the (m A , tan β) plane would correspond to one configuration of m h mass. To avoid this situation, the hMSSM benchmark was introduced [26][27][28][29]. In this scenario, the light CP-even Higgs state is enforced to be ∼125 GeV while setting the SUSY mass scale M SUSY to be rather high (i.e., > 1 TeV) in order to explain the non-observation of any SUSY particle at colliders. The hMSSM setup thus describes the MSSM Higgs sector in terms of just m A and tan β, exactly like for tree-level predictions, given the experimental knowledge of m Z and m h . In this scenario, therefore, the dominant radiative corrections would be fixed by the measured experimental value of m h which in turn fixes the SUSY scale [26][27][28][29]. It defines a largely model-independent scenario, because the predictions of the properties of the MSSM Higgs bosons do not depend on the details of the SUSY sector, somewhat unlike the previous case, wherein squark masses are fine-tuned to obtain m h ≈ 125 GeV.
The SUSY input parameters in this scenario are similar to the previous one, Equation (2), except that we take X t = 2M SUSY = 2 TeV.

Setup
Both scenarios introduced above are characterized by relatively large values of the ratio X t /M SUSY , ensuring that the MSSM mass of the SM-like Higgs state falls within the required range without the need for an extremely heavy stop. In addition, the masses of the gluino and first two generation squarks are set to 1.5 TeV, large enough to evade the current ATLAS and CMS limits stemming from SUSY searches [57][58][59][60][61]. We vary the parameters tan β and m A within the following ranges: The soft trilinear term A t is set to be equal to A b . Due to the smallness of the light quarks masses, the left-right mixing of the first two generation squarks is neglected. The gaugino mass parameters M 1 , M 2 and the soft SUSY-breaking gluino mass mg are all related through Renormalization Group Equation (RGE) running to some common high scale m 1/2 soft term which yields the relations mg ≈ 3.5M 2 and M 1 ≈ 0.5M 2 . In our analysis, we assume Grand Unified Theory (GUT) relations only between M 1 and M 2 while M 2 and mg are taken independent from each other. Finally, the soft SUSY-breaking parameters in the slepton sector have a very small impact on the predictions for the Higgs masses and production cross sections, therefore we do not report on them here.

THE 2HDM
In this section, we define the scalar potential and the Yukawa sector of the general 2HDM. The most general scalar potential which is SU(2) L ⊗ U(1) Y invariant is given by [62,63] V( 1 , 2 ) = m 2 Frontiers in Physics | www.frontiersin.org The complex (pseudo)scalar doublets i (i = 1, 2) can be parameterized as with v 1,2 ≥ 0 being the VEVs satisfying 22 GeV. Hermiticity of the potential forces λ 1,2,3,4 to be real while λ 5 and m 2 12 can be complex. In this work we choose to work in a CP-conserving potential where both VEVs are real and so are also λ 5 and m 2 12 . After EWSB, three of the eight degrees of freedom in the Higgs sector of the 2HDM are eaten by the Goldstone bosons (G ± and G) to give masses to the longitudinal gauge bosons (W ± and Z). The remaining five degrees of freedom become the aforementioned physical Higgs bosons. After using the minimization conditions for the potential together with the W ± boson mass requirement, we end up with seven independent parameters which will be taken as where, as usual, tan β ≡ v 2 /v 1 and β is also the angle that diagonalizes the mass matrices of both the CP-odd and charged Higgs sector while the angle α does so in the CP-even Higgs sector. The most commonly used versions of a CP-conserving 2HDM are the ones that satisfy a discrete Z 2 symmetry i → (−1) i+1 i (i = 1, 2), that, when extended to the Yukawa sector, guarantees the absence of Flavor Changing Neutral Currents (FCNCs). Such a symmetry would also require m 2 12 = 0, unless we tolerate a soft violation of this by the dimension two term m 2 12 (as we do here). The Yukawa Lagrangian can then be written as where Q T L = (u L , d L ) and L T L = (l L , l L ) are the left-handed quark doublet and lepton doublet, respectively, the Y f k 's (k = 1, 2 and f = u, d, l) denote the 3 × 3 Yukawa matrices and˜ k = iσ 2 * k (k = 1, 2). The mass matrices of the quarks and leptons are a linear combination of Y f 1 and Y f 2 , Y d,l 1,2 and Y u 1,2 . Since they cannot be diagonalized simultaneously in general, neutral Higgs Yukawa couplings with flavor violation appear at tree-level and contribute significantly to FCNC processes, such as M K,B,D as well as B d,s → µ + µ − mediated by neutral Higgs exchanges. To avoid having those large FCNC processes, one known solution is to extend the Z 2 symmetry to the Yukawa sector. When doing so, we ended up with the already discussed four possibilities regarding the Higgs bosons couplings to fermions [63].
Here, the shorthand notation cx ≡ cosx and sx ≡ sinx is used.
After EWSB, the Yukawa Lagrangian can be expressed in the mass eigenstate basis as follows [64,65]: We give in Table 1 the values of the Yukawa couplings κ φ f (φ = h, H, A), i.e., the Higgs boson interactions normalized to the SM vertices introduced in David et al. [66], in the four 2HDM Types. The couplings of h and H to gauge bosons W ± , Z are proportional to sin(β − α) and cos(β − α), respectively. Since these are gauge couplings, they are the same for all Yukawa types. As we are considering the scenario where the lightest neutral Higgs state is the 125 GeV scalar, the SM-like Higgs boson h is recovered when cos(β − α) ≈ 0. As one can see from Table 1, for all 2HDM Types, this is also the limit where the Yukawa couplings of the discovered Higgs boson become SM-like. The limit cos(β − α) ≈ 0 seems to be favored by LHC data, except for the possibility of a wrong sign limit [67,68], where the couplings to down-type quarks can have a relative sign to the gauge bosons ones, thus oppositely to those of the SM. Our benchmarks will focus on the SM-like limit where indeed cos(β − α) ≈ 0. We end this section by noticing that we have used the public program 2HDMC [69] to evaluate the 2HDM spectrum as well as the decay rates and BRs of all Higgs particles. We have used 2HDMC to also enforce the aforementioned theoretical constraints onto both BSM scenarios considered here.

THEORETICAL AND EXPERIMENTAL CONSTRAINTS
In order to perform a systematic scan over the parameter space of the two MSSM configurations and the four 2HDM Types, we take into account the following theoretical 3 and experimental constraints.

Theoretical Constraints
We list these here as itemized entries.
• Vacuum stability: To ensure that the scalar potential is bounded from below, it is enough to assume that the quartic couplings should satisfy the following relations [70]: We also impose that the potential has a minimum that is compatible with EWSB. If this minimum is CP-conserving, any other possible charged or CP-violating stationary points will be a saddle point above the minimum [71]. However, there is still the possibility of having two coexisting CP-conserving minima. In order to force the minimum compatible with EWSB, one need to impose the following simple condition [72]: Writing the minimum conditions as allows us to express m 2 11 and m 2 22 in terms of the soft Z 2 breaking term m 2 12 and the quartic couplings λ 1−5 . • Perturbative unitarity: Another important theoretical constraint on the (pseudo)scalar sector of the 2HDM is the perturbative unitarity requirement. We require that the S-wave component of the various (pseudo)scalar scattering amplitudes of Goldstone and Higgs states remains unitary. Such a condition implies a set of constraints that have to be fulfilled and are given by [73,74] |a ± |, |b ± |, |c ± |, |f ± |, |e 1,2 |, |f 1 |, |p 1 | < 8π, where • EW precision observables (EWPOs): The additional neutral and charged (pseudo)scalars, beyond the SM-like Higgs state, contribute to the gauge bosons vacuum polarization through their coupling to gauge bosons. In particular, the universal parameters S, T, and U provide constraints on the mass splitting between the heavy states m H , m H ± , and m A in the scenario in which h is identified with the SM-like Higgs state. The general expressions for the parameters S, T, and U in 2HDMs can be found in Kanemura et al. [75]. To derive constraints on the scalar spectrum we consider the following values for S, T, and U: while using the corresponding covariance matrix given in Baak et al. [76]. The χ 2 function is then expressed as with correlation factor +0.91. The aforementioned 2HDMC program allows us to check most of the above theoretical constraints, such as perturbative unitarity, boundedness from below of the scalar potential as well as EWPOs (S, T, and U), which are all turned on during the calculation, and can be adapted to the MSSM as well.

Experimental Constraints
The parameter space of our benchmark scenarios is already partially constrained by the limits obtained from various searches for additional Higgs bosons at the LHC and elsewhere as well as the requirement that one of the neutral scalar states should match the properties of the observed SM-like Higgs boson. We evaluate the former constraints with the code HiggsBounds [77][78][79][80] and the latter with the code HiggsSignals [81]. We stress, however, that our study of the existing constraints cannot truly replace a dedicated analysis of the proposed benchmark scenarios by ATLAS and CMS, which alone would be able to combine the results of different searches taking into account all correlations. In this section we briefly summarize the relevant features of HiggsBounds and HiggsSignals used in our study.

Collider Constraints
The code HiggsBounds tests each parameter point for 95% Confidence Level (CL) exclusion from Higgs searches at the LHC as well as LEP and Tevatron. First, the code determines the most sensitive experimental search available, as judged by the expected limit, for each additional Higgs boson in the model. Then, only the selected channels are applied to the model, i.e., the predicted signal rate for the most sensitive search of each additional Higgs boson is compared to the observed upper limit. In the case the prediction exceeds the limit, the parameter point is regarded as excluded. For more details on the procedure, the reader can see Bechtle et al. [80].
Among the searches that are relevant in constraining our scenarios for charged Higgs studies, the version we have used, 5.2.0beta, of HiggsBounds includes the following.     uncertainty estimates for the gluon-fusion and bottomquark annihilation cross sections. The most conservative (i.e., weakest) determination of the exclusion region is obtained by taking simultaneously the lowest values in the uncertainty range for both production processes of each of the heavier Higgs bosons, while the least conservative (i.e., strongest) determination is obtained by taking simultaneously the highest values in the uncertainty range.
With the use of the code HiggsSignals, we test the compatibility of our scenarios with the observed SM-like Higgs signals, by comparing the predictions of SusHi, FeynHiggs, and 2HDMC for the signal strengths of Higgs production and decay in a variety of channels against ATLAS and CMS measurements.

DM Constraints
These have naturally been enforced only in the MSSM case, by using the program micrOMEGAs version 5.0.9 [114] 4 . Such a code calculates the properties of DM in terms of its relic density as well as its direct and indirect detection rates. For the two MSSM scenarios considered here, the DM candidate, i.e., the Lightest Supersymmetric Particle (LSP), is the lightest neutralino. We require that the outcome of the calculation of the relic density should be in agreement with the latest Planck measurement [115].

NUMERICAL RESULTS
In this section, we present our findings for the MSSM and 2HDM in turn. 4 https://lapth.cnrs.fr/micromegas/

MSSM Results
In the hMSSM scenario, all superparticles are chosen to be rather heavy so that production and decays of the MSSM Higgs bosons are only mildly affected by their presence due to the decoupling properties of SUSY. In particular, the loopinduced SUSY contributions to the couplings of the light CPeven scalars are small and the heavy Higgs bosons with the masses even up to 2 TeV decay only to SM particles. Therefore, the phenomenology of this scenario at the LHC resembles that of a 2HDM Type-II with MSSM-inspired Higgs couplings and mass relations.
The masses of the third generation squarks and that of the gluino are safely above the current bounds from direct searches at the LHC, as intimated. Specifically, we refer to CMS Collaboration [116,117] [129] for the gluino. The value chosen for X t is close to the one for which the maximal value of m h is obtained. The m mod+ h scenario is very similar to the hMSSM one except the fact that we take X t = 2M SUSY = 1 TeV.
In Figure 1 the allowed regions on the (m A , tan β) plane are depicted for various χ 2 , wherein the left and right panel are, respectively, for the hMSSM and m mod+ h scenarios. For the hMSSM and χ 2 ≤ 12, one can see that m A should be heavier than about 400 GeV. In the case of m A ≈ 400 − 600 GeV, tan β should be in the range [1,9] while for m A around 1 TeV tan β is in the range [1,15]. The dashed (solid) line represents the 95% (68%) CL obtained by the HiggsSignals fit and the best fit point is located at m A ≈ 1 TeV and tan β ≈ 2. For the m mod+ h scenario, the situation is quite different. In order to accommodate m h ≈ 125 GeV, one needs tan β > 10. Similarly to the left panel, also in the right one the dashed(solid) line represents the 95% (68%) CL obtained by the HiggsSignals fit and the best fit point is located at m A ≈ 1 TeV and tan β ≈ 20. For this scenario and for χ 2 < 12, all tan β ≤ 10 are excluded. Note that after imposing DM constraints onto the best fit analysis, we observe that the best fit points for both the hMSSM and m mod+ h scenarios move to somewhat lighter values of the charged Higgs boson mass. The BPs given in Table 2 account for this effect.
In Figure 2 we present the total width of the charged Higgs boson, again, over the (m A , tan β) plane, for both hMSSM (left) and m mod+ h (right). As one can see from the left panel, the total width for the hMSSM case is largest for tan β ≤ 3, which is when Ŵ H ± ≈ 7 − 10 GeV, while for tan β ≥ 5 the width drops to 1-3 GeV. This effect can be attributed to the fact that the total width is fully dominated by H + → tb, whenever this channel is open, in which the top effect is more pronounced for low tan β. In this case, H + → τ ν is subleading and also the decay modes H + →χ + iχ 0 j are suppressed. In the case of m mod+ h , since small tan β is not allowed, the total charged Higgs boson width is generally smaller than in the hMSSM case, as a consequence of the fact that H + → tb is therefore smaller in this scenario. The maximal total width is here obtained for m A ≈ 1 TeV and a large tan β ≈ 20. In the m mod+ h scenario, the decay H + →χ + 2χ 0 2 could have a significant BR, reaching 30%. Hence, the H ± is always rather narrow, whichever its mass. In fact, owing to the degeneracy between m A and m H ± in the MSSM, as dictated by h data, a remarkable result is that in the minimal SUSY scenario a charged Higgs boson is essentially always heavier than the top quark.
In Figure 3 we show the production cross section for single charged Higgs boson production in association with a top quark (as appropriate for the m H ± > m t case) times the BR of H + into a specific final state for both the hMSSM and m mod+ h scenarios using Prospino [130][131][132]. In fact, as we have seen previously, the total width of the charged Higgs state is rather small in both cases, in relation to the mass, so that one can use the Narrow   As for BR(B 0 q → µ + µ − ), the combined results from the LHCb and CMS collaborations are shown as in Archilli [136]. Width Approximation (NWA) to estimate such a cross section (which we have done here). In the top-left(top-right) panel of Figure 3, we show the size of the cross section of σ (pp → tH − + c.c.) × BR(H ± → tb)(σ (pp → tH − + c.c.) × BR(H ± → τ ν)), given in pb.
For the hMSSM scenario, one can see that in the tb channel the largest cross section (more than 0.1 pb) is reached for small tan β < 3. There is also a wide region with m H ± ∈ [400, 600] GeV and tan β < 10 where the cross section is still rather important: between 10 −3 and 0.1 pb. As for the τ ν channel, the cross section is maximized when tan β is in the range [4,9] and the largest cross section is seen around 10 −3 pb. However, amongst the bosonic channels, H ± → W ± A is hopeless because BR(H ± → W ± A) is very suppressed while H ± → W ± h can have a rate that is close to 10 −2 pb for small tan β ≈ 1. Note that, for completeness, we have also drawn the exclusion region due to BR(B → X s γ ), even though we can always assume some kind of flavor violation that takes place in the MSSM and can bring the BR(B → X s γ ) to a correct value. In terms of σ (pp → tH − + c.c.) × BR(H ± → XY) for the m mod+ h scenario, the situation is worse. The best channels are H + → tb and H + →χ + 1χ 0 1 with the maximum cross section in the allowed region being between 10 −3 and 10 −2 pb for charged Higgs boson masses in the range 400-600 GeV, as can be seen from Figure 4.
We conclude this section by presenting in Table 2 two BPs, one each for the m mod+ h and hMSSM scenarios, to aid future analyses of Run-2 (and possibly Run-3) data from the LHC. Notice that these BPs do not correspond to the best fit points in these two MSSM configurations, as the latter would yield too small cross sections 5 , owing to the very large charged Higgs mass involved (of order 1 TeV). Yet, the BPs presented correspond to rather large values of m H ± , as dictated by the compatibility tests of the m mod+ h and hMSSM scenarios with current datasets, still giving production and decay rates (in one or more channels) potentially testable in the near future.

2HDM Results
We now move on to discuss the 2HDM. In this scenario, we consider h as being again the 125 GeV SM-like Higgs and vary the other six parameters as indicated in Table 3. When performing the scan over the 2HDM parameter space, other than taking into account the usual LHC, Tevatron and LEP bounds (as implemented in HiggsBounds and HiggsSignals) as well as the theoretical ones (as implemented in 2HDMC), we also have to consider flavor observables. In fact, unlike the MSSM, where potentially significant contributions to (especially) B-physics due to the additional Higgs states entering the 2HDM beyond the SMlike one can be canceled by the corresponding sparticle effects (and besides, are generally small because of the rather heavy H, A, and H ± masses), the 2HDM has to be tested against a variety of such data. The B-physics observables that we have considered to that effect are listed in Table 4. We have computed the 2HDM predictions for these in all 2HDM Types using our 5 Probably accessible only at the High-Luminosity LHC [133]. own implementation, which output in fact agrees with the one from SuperIso [137] (when run in 2HDM mode).
Based on such constrained scans, we first illustrate in Figure 5, on the (cos(β − α), tan β) plane, the best fit points for the four 2HDM Types. Herein, are also shown the compatibility regions with the observed Higgs signal at the 1σ (green) and 2σ (yellow) level. The details of the best fit points herein (red stars) are given in Table 5 together with the values of the following observables: the total charged Higgs width Ŵ H ± , σ (pp → tH − + c.c.), BR(H ± → τ ν), BR(H ± → tb), BR(H ± → AW ± ), and BR(H ± → hW ± ). Note that in the 2HDM Type-II and -Y, the best fit point is located at a charged Higgs mass around 600 GeV because of the B → X s γ constraints while in the 2HDM Type-I and-X one can fit data with a rather light charged Higgs state.
In Figures 6-9, we show (in gray) over the (m H ± , tan β) plane the 95% CL exclusion region from the non-observation of the additional Higgs states for 2HDM Type-I(-II)[-X]{-Y}. In all these plots, we also draw (as a solid yellow line) the 95% CL exclusion from BR(B → X s γ ) together with a solid green line representing the 1σ compatibility with the Higgs signals observed at the LHC. As a green star, we also give the best fit point to these data over the available parameter space for all Types (these are the same as the red stars in the previous figure). It is clear from these plots that, in the 2HDM-I and -X, one can still have relatively light charged Higgs states (of the order 100-200 GeV in mass) that are consistent with all aforementioned data, crucially including B-physics observables. In addition, such light charged Higgs does not affect too much the rate of h → γ γ which is strongly dominated by the W ± loops while the charged Higgs loops are subleading. In the case of the 2HDM Type-II and -Y, the BR(B → X s γ ) constraint pushes the charged Higgs boson mass to be higher than 580 GeV. (Note that, in the 2HDM Type-II, it is clear that, like for the MSSM case, large tan β is excluded mainly from H, A → τ + τ − as well as from H + → τ ν searches at the LHC). However, for 2HDM Type-X, one can see that light charged Higgs states, with m H ± ≤ 170 GeV, are excluded for all tan β's and this is due to charged Higgs searches failing to detect H ± → τ ν.
We now discuss the size of the charged Higgs production cross section times its BRs in decay channels, such as FIGURE 6 | The σ (pp → tH − + c.c.) × BR(H ± → XY ) rate (in pb) at √ s=14 TeV in the 2HDM Type-I, for XY ≡ tb (top left), XY ≡ τ ν (top right), XY ≡ AW ± (bottom left), and XY ≡ hW ± (bottom right). Exclusion bounds at 95% CL from the non-observation of the additional Higgs states are overlaid in gray. The green contour indicates compatibility with the observed Higgs signal at 68% CL and the best fit (benchmark) points are marked by green stars. The solid yellow line contours are the boundary of 95% CL exclusion from B → X s γ measurements. The maximum of the cross section σ (pp → tH − + c.c.) × BR(H ± → XY ) is 3.1 and 1.83 pb for XY ≡ AW ± , tb, respectively.
H + → τ ν, tb, AW ± and hW ± . In Figure 6 (top-left panel) we illustrate the values of σ (pp → tH − + c.c.) × BR(H ± → tb) (in pb) where we can see that it is possible tohave a production times decay rate in the range 0.01 to 0.2 pb for 1 ≤ tan β ≤ 6 and 180 GeV < m H ± < 300 GeV. This could lead to more than thousands raw ttb signal events for 100 fb −1 luminosity. In the case of H ± → τ ν and H ± → hW ± , which are suppressed, respectively, by 1/ tan β and cos(β − α) ≈ 0, the rate is much smaller than for the tb mode. In contrast, since the coupling H ± W ∓ A is a gauge coupling without any suppression factor, when H ± → AW ± is open, it may dominate over the H ± → tb channel. One can see from Figure 6 (bottom-left panel) that, for 100 GeV < m H ± < 220 GeV and for 1 ≤ tan β ≤ 14, the corresponding rate for σ (pp → tH − + c.c.) × BR(H ± → AW ± ) ≥ 0.01 pb. This could lead to an interesting final state bW + W − A where one W ± could decay leptonically, hence offering a clean trigger.
In the case of 2HDM Type-II and -Y, as one can see from Figures 7, 9, respectively, there is a wide region over the (m H ± , tan β) plane where the rate for σ (pp → tH − + c.c.) × BR(H ± → tb) is rather sizable for both moderate (m H ± ≤ 300 GeV) and heavy (otherwise) charged Higgs masses (top-left panel). However, if one takes into account the B → X s γ constraint, then m H ± is required to be much heavier than 580 GeV (as already discussed), which makes FIGURE 7 | The σ (pp → tH − + c.c.) × BR(H ± → XY ) rate (in pb) at √ s=14 TeV in the 2HDM Type-II, for XY ≡ tb (top left), XY ≡ τ ν (top right), XY ≡ AW ± (bottom left), and XY ≡ hW ± (bottom right). The maximum of the cross section σ (pp → tH − + c.c.) × BR(H ± → XY ) is 2.3 and 1.38 pb for XY ≡ AW ± , tb, respectively. The color coding is the same as in Figure 6.
the rate σ (pp → tH − + c.c.) × BR(H ± → tb) ≥ 0.1 pb only for tan β < 1.5. All the other channels (in the three remaining panels) have smaller production times decay rates.
The 2HDM Type-X is depicted in Figure 8, wherein the usual production times BR rates are shown. The top-right panel is again for the H + → tb channel, which exhibits a potentially interesting cross section (≥ 1 fb) in the H + → tb channel for both a light charged Higgs mass (around 200 GeV) and a heavy one (around 420 GeV). In the case of the τ ν channel (top-right panel), one can get sizable rates for σ (pp → tH − + c.c.) × BR(H ± → τ ν) for a charged Higgs mass around 200 GeV and tan β ≥ 2.
In all 2HDM Types, we elect the best fit points to also be the BPs amenable to experimental tests by ATLAS and CMS.

CONCLUSIONS
We have studied charged Higgs boson phenomenology in both the MSSM and 2HDM, the purpose being to define BPs amenable to phenomenological investigation already with the full Run-1 and -2 datasets and certainly accessible with the Run-3 one of the LHC. They have been singled out following the enforcement of the latest theoretical and experimental constraints, so as to be entirely up-to-date. Furthermore, they have been defined with the intent of increasing sensitivity of dedicated (modeldependent) H ± searches to some of the most probable parameter space configurations of either scenario. With this in mind, we have listed in two tables their input and output values, the former in terms of the fundamental parameters of the model concerned and the latter in terms of key observables (like, TeV in the 2HDM Type-X, for XY ≡ tb (top left), XY ≡ τ ν (top right), XY ≡ AW ± (bottom left), and XY ≡ hW ± (bottom right). The maximum of the cross section σ (pp → tH − + c.c.) × BR(H ± → XY ) is 2.3 and 1.23 pb for XY ≡ AW ± , tb, respectively. The color coding is the same as in Figure 6. e.g., physical masses and couplings, production cross sections and decay BRs). We have also specified which numerical tools we have used to produce all such an information, including their settings.
For the MSSM we have concentrated on two popular scenarios, i.e., the hMSSM and m mod+ h ones. It was found that the hMSSM case still possesses a rather large available parameter space, here mapped over the (m A , tan β) plane, while the m mod+ h one is instead much more constrained. In terms of the largest production and decay rates, in the hMSSM scenario one finds that the most copious channels, assuming pp → tH − + c.c. production, are via the decay H + → tb followed by H + → τ ν whereas for the m mod+ h scenario the decay modes H + → tb and H + →χ + 1χ 0 offer the largest rates. In both cases, only m H ± > m t values are truly admissible by current data.
Within the 2HDM, we have looked at the four standard Yukawa setups, known as Type-I, -II, -X, and -Y. Because of B → X s γ constraints, the profile of a charged Higgs in the 2HDM Type-II and -Y is a rather heavy one, with a mass required to be more than 580 GeV. While this puts an obvious limit to LHC sensitivity owing to a large phase space suppression in production, we have emphasized that H ± → bbW ± channels should be searched for, with intermediate contributions from the AW ± and tb modes (including their interference [138]), alongside H ± → τ ν. In the case of the 2HDM Type-I and -X, a much lighter charged Higgs state is still allowed by data, in fact, even with a mass below that of the top quark. While the configuration m H ± < m t is best probed by using tt production and decays into τ ν, the complementary mass region, i.e., m H ± > m t (wherein pp → tH − + c.c. is the production mode), may well FIGURE 9 | The σ (pp → tH − + c.c.) × BR(H ± → XY ) rate (in pb) at √ s=14 TeV in the 2HDM Type-Y, for XY ≡ tb (top left), XY ≡ τ ν (top right), XY ≡ AW ± (bottom left), and XY ≡ hW ± (bottom right). The maximum of the cross section σ (pp → tH − + c.c.) × BR(H ± → XY ) is 3.54 and 1.85 pb for XY ≡ AW ± , tb, respectively. The color coding is the same as in Figure 6.
be accessible via a combination of H + → tb and H ± → AW ± (in Type-I) plus H ± → τ ν (in Type-X).

DATA AVAILABILITY STATEMENT
The datasets analysed in this study can be found in the HiggsBounds repository: https://higgsbounds.hepforge.org.

AUTHOR'S NOTE
Since the original submission of this paper, several new experimental analyses have been carried out by ATLAS and CMS using the full Run-2 data sample of ≈ 139 fb −1 .
Some of these, covering both measurements of the SM-like Higgs Boson and the search for new (pseudo)scalar Higgs states, both charged and neutral, have been captured by the latest versions of HiggsBounds and HiggsSignals, HiggsBounds-5.3.2beta and HiggsSignals-2.2.3beta, respectively. Likewise, further analyses by LHCb of flavor observables have been carried out since and most of these have been captured by the latest version of SuperIso. Hence, we have repeated our scans using all such tools and found negligible differences between our original results and the new ones. Further, we have investigated which ones of the full Run-2 data set analyses were not incorporated in the above codes and found that their ad-hoc application to our analysis did not change our results either.