Probing the $H^\pm W^\mp Z$ interaction at the high energy upgrade of the LHC

An $H^\pm W^\mp Z$ interaction at the tree level is common feature of new physics models that feature scalar triplets. In this study, we aim to probe the strength of the aforementioned interaction in a model-agnostic fashion at the futuristic 27 TeV proton-proton collider. We assume that the $H^\pm$ couples dominantly to ($W^\pm,Z$) and ($t,b$). We specifically study the processes that involve the $H^\pm W^\mp Z$ vertex at the production level, that is, $p p \to H^\pm j j$ and $p p \to Z H^\pm$. Moreover, we look into both $H^\pm \to W^\pm Z,~t b$ decays for either production process. Our investigations reveal that the $H^\pm j j$ production process has a greater reach compared to $Z H^\pm$. Moreover, the discovery potential of a charged Higgs improves markedly with respect to the earlier studies corresponding to lower centre-of-mass energies. Finally, we recast our results in the context of the popular Georgi-Machacek model.


I. INTRODUCTION
With the discovery of a Higgs boson of mass 125 GeV [1,2], the particle spectrum of the Standard Model (SM) is complete. Moreover, the properties of the discovered boson are found to be increasingly consistent with that of the SM Higgs. Despite this success, the SM remains far from being the complete framework. Certain issues on both theoretical and experimental fronts advocate for additional dynamics beyond the SM. Interestingly, extending just the scalar sector of the SM can suffice to address such issues. This puts forth extended Higgs sectors as prototypes of new physics (NP).
A subset of extended Higgs sectors features additional SU (2) L scalar multiplets and all of them predict at least one singly charged scalar H + . While the tree level fermionic couplings of H + are generally proportional to the fermion mass for all SU (2) L representations, it is the coupling to the gauge bosons that can differ. For instance, the H + W − Z interaction vanishes at the tree level for a two-Higgs doublet model (2HDM) [3][4][5][6]. This coupling is generated at one-loop [7][8][9][10] and hence is generally small in magnitude. And this is true regardless the number of such scalar doublets. On the contrary, the aforementioned coupling is non-zero at the tree level itself in case of SU (2) L scalar triplets [11][12][13][14] and is proportional to the vacuum expectation value (VEV) acquired by the neutral component of the triplet. The triplet VEV is however tightly constrained by measurements of the ρ-parameter for the simplistic Higgs triplet model (HTM) (that employs a complex scalar triplet over and above the doublet) thereby obscuring the observability of the H + W − Z coupling. That said, there exist nonminimal extensions of the HTM where this problem can be circumvented. These involve combining a real scalar triplets with a complex ones in a custodially symmetric manner such that the ρ-parameter is unity for all values of the triplet VEV. The most popular of such extensions arguably is the Georgi-Machacek (GM) model (see [15][16][17][18][19][20][21][22][23][24][25][26][27][28][29][30][31][32][33][34] for a partial list of references) that features an SU (2) L × SU (2) R global symmetry in the scalar potential. This framework predicts two charged Higgses both of which couple to W + , Z. In fact, they also get to interact with the SM fermions by virtue of mixings induced by the scalar potential. Given that the strengths of such Yukawa interactions are proportional to the corresponding fermion masses, the one with the (t, b) pair is of foremost importance here. The charged Higgses emerging from the more non-minimal extensions also naturally feature sizeable couplings to (W, Z) and (t, b). Therefore, owing to the fact that the H + W − Z interaction features in a plethora of new physics (NP) scenarios, there is the enticing possibility of probing the same in a model-independent fashion.
A singly charged Higgs lighter than the t-quark has been looked for at the LHC via the t → H + b decay. On the other hand, a heavier one is searched in processes where it is produced in association with a t and subsequently decays to a pair of fermions. Such channels include H + → tb [35][36][37], H + → τ + ν τ [38,39], H + → cb [40] and H + → cs [41]. These searches mainly rely on the strength of the Yukawa interactions of H + . In addition, bosonic decays of the form H + → W + X, where X denotes a scalar have also been of substantial interest [42][43][44][45][45][46][47]. On the other hand, the H + W − Z interaction has been probed in pp → H ± jj, H ± → W ± Z [48,49] in context of the GM model and scalar triplet extensions of the minimal supersymmetric standard model (MSSM) [50][51][52][53]. A recent review on charged Higgs phenomenology in 2HDM is [54].
In this work, we set out to study the H + W − Z interaction at the proposed 27 TeV energy upgrade [55] of the Large Hadron Collider (LHC), the so-called high-energy LHC (HE-LHC). We adopt a generic but simplistic framework comprising H + W − Z and H + tb interactions only. We look for the H + in processes that involve the H + W − Z coupling at the production vertex. Two such important processes are W Z fusion and the W -mediated pp → ZH ± . Following production, the H ± → W ± Z, tb decay modes are looked at and bounds on the corresponding branching fractions are derived corresponding to 15 ab −1 integrated luminosity. The generality of the adopted framework makes the obtained results applicable to a wide variety of models that consist of at least one charged Higgs.
The study is organised as follows. We propose the theoretical framework in section II.
The existing LHC constraints and our analysis strategy are discussed in section III. Section IV contains the analyses for the proposed channels. The results so obtained are recasted as favoured regions in the model parameter space in section V. Finally we summarise and conclude in section VI.

II. THEORETICAL FRAMEWORK
The following Lagrangian describes the dimension ≤ 4 interactions of an H ± in a generic fashion [56][57][58] 1 : In the above, F, A tb and B tb are dimensionless parameters quantifying the strengths of the H ± W ∓ Z and H ± tb interactions respectively. To generate some perspective, we comment here on the typical values for F , A tb and B tb for two specific classes of models. In scalar sectors containing SU (2) L doublets alone, F is generated radiatively at the one-loop level.
The largest |F | for a 2HDM is reported to be 0.01 [59] with further enhancement expected 1 Momentum dependent interactions borne out of higher dimensional operators [58] have not been considered in this study.
with increasing the number of scalar doublets. For example, addition of an additional inert (color-octet) doublet elevates the maximum value of |F | to 0.03 (0.025) [59,60]. While this trend is encouraging, one must note that arbitrarily increasing the number of scalar doublets not only is limited by the experimental constraints such as that of the electroweak T -parameter and the diphoton signal strength, the model becomes aesthetically unattractive and loses predictive power. Therefore, |F | is in the O(10 −2 ) ballpark for realistic multidoublet extensions of the SM. The other parameter of interest, A tb equals cotβ for a 2HDM and its inert doublet and color-octet doublet extensions. Here, tanβ = v 2 v 1 with v 1 and v 2 denoting the VEVs of the two active doublets. Since tanβ < 1 is ruled out by flavour constraints, one has A tb < 1 for a 2HDM. In fact, tanβ ∈ [1,10] are typical values allowed for the pure 2HDM (see [61] for a global analysis of the 2HDM parameter space) leading to On the other hand, B tb equals −cotβ (tanβ) for a Type-I (Type-II) 2HDM.
The other class of models are characterised by SU (2) L scalar triplets where F is non-zero at the tree level itself. One derives F gv ∆ c W M W for a single complex triplet. Here, v ∆ , M W , c W and g denote the triplet VEV, W -mass, cosine of the Weinberg mixing angle and the SU (2) L gauge coupling respectively. In this case however, the ρ-parameter constraint dictates v ∆ 5 GeV [62] leading to F 0.017. The singly charged scalar H + couples to quarks by virtue of doublet-triplet mixing due to which A tb = B tb 2v ∆ v φ where v φ denotes the VEV of the SM-like doublet. And the maximum value permitted by the ρ-parameter constraint is A tb 0.030. The small v ∆ for this model therefore renders both |F |, |A tb | O(0.01). The same is not true in case of the GM model where v ∆ is allowed to be larger by virtue of the global SU (2) L × SU (2) R symmetry of the scalar potential. A global analysis of the minimal GM model in [63] gives the bound v ∆ 35 GeV thereby indicating F 0.32 and A tb 0.56.
Other possible bosonic and fermionic interactions of H ± are not included in the present study. Motivated by the some NP scenarios such as the Type-I 2HDM and the GM model, we further take M b B tb << M t A tb in this work for the subsequent analysis. Therefore, with the framework being describable by the two interactions, the total decay width of H + reads The branching ratios read Finally, we give the expressions for the corresponding partial widths below.
Here, x P =

III. ANALYSIS STRATEGY AND EXISTING LHC LIMITS
Two kinematically distinct topologies for H ± production in pp-collisions that involve the H ± W ∓ Z interaction are pp → H ± jj and pp → ZH ± . In the former, an H ± is produced by the fusion of W ± and Z and two forward jets are emitted. This is essentially vector boson-fusion (VBF) that has close semblance to similar processes for h-production, where, h denotes the SM-like Higgs of mass 125 GeV. On the other hand, the ZH ± process is a W -mediated s-channel topology. In this work, we aim to probe the H + → W + Z, tb decays for both the aforementioned production processes. Leptonic decays of t, W ± and Z gives rise to the following four cascades: pp → H ± Z → tbZ → 2b + 3l + / E T , where l = e, µ. The corresponding Feynman diagrams are shown in Fig.1.
We reiterate for clarity that both processes (a) and (b) describe H ± -production via VBF with subsequent decays to W ± Z and tb modes respectively. Similarly, processes (c) and (d) correspond to the W ± Z and tb decays of H ± produced via pp → ZH ± .
As mentioned before, the framework can be described by M H + , and the couplings F and A tb . We aim to explore the discovery potential of a 27 TeV pp-machine with integrated luminosity L = 15 ab −1 through a detailed signal and background analysis of the aforementioned signal channels. Further, we supplement the conventional cut-based analyses (CBA) with the more sophisticated multivariate analyses (MVA) using the Boosted Decision Tree Decorrelated (BDTD) algorithm [64]. An overview of BDTD analysis will be given in subsection The performances of the proposed channels can be compared by drawing the exclusion and discovery contours corresponding to each in the F − A tb plane for a given M H + . It is imperative to discuss possible exclusions on {M H + , F, A tb } from previous collider searches. The proposed process (a) itself has been searched at the LHC by the ATLAS and pp → H ± Z → tbZ → 3l + 2b + / E T . In the diagram, j represents the a light quark jet, while l ± denote e, µ.
CMS collaborations [48,49] can be obtained as, Here σ |F |=1 | M H + denotes the pp → H ± jj cross section in our framework for |F | = 1 for a given M H + in absence of kinematical cuts.
Another pertinent search by ATLAS at 13 TeV and 36.2 fb −1 of integrated luminosity is pp → tbH ± , H ± → tb [35]. A more stringent bound comes from the recent search at 139 fb −1 data [65]. Adopting the more recent bound, a limit on σ tbH ± × BR(H + → tb) is derived as : where, σ |A tb |=1 | M H + denotes the pp → tbH ± cross section for |A tb | = 1 for a given M H + in absence of kinematical cuts.  [68]. The multivariate analysis is done using the TMVA package [69]. The signal significance S is derived using [70], with N S (N B ) denoting the number of signal (background) events surviving the kinematical cuts 2 .
In order to estimate the effects on the final signal significance by assuming a systematic uncertainty of σ sys_un , the signal significance formula changes as: where σ B = σ sys_un × B.
We set out to analyse the various channels following the aforementioned strategy. For processes (a),(c) and (b),(d) in Fig.1, (F ,BR(H + → W + Z)) and (F ,BR(H + → tb)) are respectively treated as the free parameters. That said, for a given signal, it is important to 2 The number of signal (N S ) and the background (N B ) events can be calculated as: where σ S(B) , L, S ( B ) denote the signal (background) cross section, integrated luminosity and signal (background) cut-efficiency respectively.
compare the results from CBA to those from MVA for specific reference values of F and the branching ratios. Therefore, we take (F ,BR(H + → W + Z)) = (0.4,0.4) and (F ,BR(H + → tb)) = (0.4,0.4) as the reference for processes (a),(c) and (b),(d) respectively for the ensuing analysis.
A. The 3l + 2j + / E T channel This subsection discusses the VBF production of H ± followed by the H ± → W ± Z decay.
The VBF topology always leads to a couple of light jets in the forward and backward directions with negligible hadronic activity in the intervening rapidity gap. Since the two jets (j 1 , j 2 ) reside in different η-hemisphere for VBF production, the forward (backward) jet is identified by η > 0 (η < 0). We therefore demand η j 1 η j 2 < 0 to tag the two leading light jets as forward and backward ones. The H ± → W ± Z decay ultimately leads to three leptons, two of which are of same flavor and opposite sign (SFOS) that come from Z and the third one comes from W ± . In addition, a neutrino also originates from W ± . We choose to treat the H ± → W ± Z branching ratio as a free parameter at this level. The signal cross sections for different M H + are given in Table I.
The most dominant SM backgrounds corresponding to signal come from W ± Zjj and ZZjj production followed by leptonic cascades. These two background processes involve jets at the production level as does the signal. Therefore, we demanded p j T > 20 GeV, M j 1 j 2 > 500 GeV, |η j | < 5 while generating these in order to pick out the phase space kinematically relevant to the VBF topology. The resulting LO cross sections of these two background processes along with their leptonic decay chains are given in Table I. Sub-leading contributions to the backgrounds come from W + W − Z, W ± ZZ, ZZZ, ttW ± , ttZ. We do not impose cuts at the generation level for the sub-leading ones since they do not involve jets at the production level, and, the corresponding cross sections are again listed in Table I. A couple of observations then emerge. First, W ± ZZ → 4l + 2j, ZZZ have negligible cross sections compared to the leading backgrounds. Secondly, though it might naively appear that the other sub-leading channels offer rates comparable to, or, even higher than ZZjj, one must remember that their cross sections are computed in the absence of kinematical cuts. Therefore, stringent VBF selection cuts almost nearly obliterate ttW ± , ttZ, W ± ZZ → 3l + 2j + / E T . In all, it is only the W ± Zjj and ZZjj channels that are kinematically relevant to the analysis.
Events are selected by demanding at least two light jets (N j ≥ 2) and three leptons after vetoing b-jets and τ -jets in the final state. Of the two hardest light jets, one is demanded to be forward and the other as the backward one as mentioned earlier. We first initiate the cut-based approach. The following trigger-level cuts are applied: Signal / Backgrounds Process Cross section σ (LO) (fb) To start with, we shall compute the signal significance using the cut-based method. In addition to the trigger cuts, certain specific kinematic variables are identified to extract the signal with higher efficiency. They are We denote the corresponding cuts to be A 1 , A 2 , A 3 , A 4 respectively and they are optimised by looping over a few configurations and selecting the one that yields the maximum statistical significance.
The cut-flow for the signal and background processes corresponding to the chosen BPs are displayed in Table III-VI. The following discussion motivates these variables and compares the signal and backgrounds in their terms.
• A 1 and A 2 : A VBF topology demands a high invariant mass of the two forward jet system (M j 1 j 2 ) and a large separation between their pseudo rapidities (|∆η j 1 j 2 |). Thus high invariant mass cuts (given in Table II) are used for all the benchmarks. Strong cuts on |∆η j 1 j 2 | for enhancing the significance for all benchmark points can be found in Table II. Corresponding normalised distributions for both signal and backgrounds are drawn in Fig.2(d) and Fig.2(e). It can be observed that the peak of M j 1 j 2 distribution shifts towards higher value with increase in M H + . From the distributions, it is evident that these cuts reduce the backgrounds to a large extent. The same is also reflected in the first two columns of Table III, Table IV and Table V. • A 3 : Next comes a cut on the invariant mass of the 3l + / Comparing these two variables with M H + , the closest one is chosen to be M W Z inv . For clarity we have plotted distributions of (M W Z inv ) 1 and (M W Z inv ) 2 for both signal and dominant backgrounds before applying any cut. Whereas, for backgrounds since the source of 3l + / E T is not a single particle, the corresponding Fig.2 • A 4 : Another potentially important variable for this signal is the transverse mass of W Z-system. It can be constructed as: where p T (W ) is computed by adding p miss T and p T of the lepton (coming from W ± ) vectorially and E T (W ) is obtained from the scalar sum of p miss T and the lepton transverse energy. p T (Z) is constructed by adding the p T of the two leptons of SFOS that peak near the Z-mass and E T (Z) is calculated using the following relation.
As the W ± , Z coming from heavy H ± are highly boosted, peak of the normalised distribution of M W Z T will be shifted to the higher end starting from BP1 to BP4 in  Table II.
Overall, the CBA shows that M H + = 500 GeV has the highest observability. Next we perform an MVA to improve over the results obtained from cut-based analysis. Before going . / E T,x and / E T,y are the x and y component of the missing transverse momentum. 4 While drawing the distributions, we only have shown the distribution for two dominant backgrounds, i.e. W ± Zjj and ZZjj.
Signal 862 654 491 262   First the signal and the background-like events are distinguished using decision trees as classifier. To segregate the signal and background-like events depending on its purity for chosen benchmark points for 3l + 2j + / E T channel. 5 , each node of the decision tree is associated with one discerning variable on which an optimised cut value is imposed. By modifying the BDTD variable NCuts, one can perform the aforementioned task within TMVA. Starting from the zeroth node or root node, the training of the decision trees is carried on unless a particular depth (MaxDepth) specified by the user is achieved. Finally from the final leaf nodes, an event can be tagged as signal or background according to the purity.
The decision trees being prone to statistical fluctuations of the training sample, are said to be weak classifiers. This problem can be alleviated by combining a set of weak classifiers into a stronger one through the modification of the weight of the events. One can thus create new decision trees using this particular procedure termed as Boosting. Several kinematic variables might be highly correlated. We select a reasonable number of top ranking variables among them for providing best possible separation between the signal and backgrounds. Thus for this channel following kinematic variables with maximum discerning abilities are proposed for MVA: Here, ∆R ij denotes the distance in the η − φ plane between the ith and jth particle. While is the scalar (vector) sum of the transverse momentum of the system of ith, jth and kth particle, M j 1 j 2 l 1 l 2 l 3 represents invariant mass of "j 1 , j 2 , l 1 , l 2 , l 3 "-system. The rest of the variables have been defined earlier in the text. According to the BDT ranking, among these variables, the five most important variables to differentiate the signal from backgrounds are :  Table VII for convenience.
To maximise the significance, one can adjust BDT cut value or BDT score. The BDT cut values adjusted for four benchmarks can be found in Table VII. In Fig.3(a), we have presented the variation of the signal significances with BDT score. Here it can be easily inferred that the significance attains a maximum value at a particular BDT score 6 . Next we draw Receiver's Operative Characteristic (ROC) curves 7 for all benchmarks in Fig.3 to estimate the degree of rejecting the backgrounds with respect to the signal. From Fig.3, it can be inferred that the degree of rejecting backgrounds reduces with decrease in M H + .
But owing to large signal cross section of BP1 compared to the rest of the three BPs, BP1 fares the best in probing this channel.
The yields for signal (for M H + = 200 GeV, 300 GeV, 500 GeV, 1 TeV) and corresponding backgrounds after optimisation through BDTD-analysis are given in Total 3632 ZZjj 33 Total 1152  In the previous subsection we explored VBF production of the charged Higgs H ± where one forward and one backward light jet is produced along with H ± and H ± further decays into W ± and Z. In this subsection, the H ± → tb decay is probed with the t decaying leptonically leading to a final state containing two b-tagged jets, two light jets, one isolated The cross sections for the chosen signal benchmarks and background processes are given in Table IX. We impose the following generation-level cuts while generating these backgrounds to better handle the statistics for large backgrounds :  In addition to the generation-level cuts we add stronger trigger-level cuts on |η j,b,l | as : • B 4 : We construct an observable called M ef f which is defined as the scalar sum of transverse momenta of all visible particles and / E T . This variable helps us separate signal from background especially in the high M H + region, because the decay products of a heavy charged Higgs will naturally be boosted rendering higher M ef f for the whole system.
• B 5 : We also utilise the contransverse mass variable, M CT (a, b) [71] which is an invariant quantity for two objects a and b having contra-linear (opposite in direction) and equal magnitude boost. This variable is defined as, Here, M (a) and p T (a) corresponds to the invariant mass and transverse momentum vector respectively for the object a, whereas the transverse energy, E T is defined as such as M 2blν , we have utilised this variable as an input for our BDTD analysis which will be discussed shortly.
The optimised cuts on the aforementioned variables for each of the benchmark points are given in Table. X. These cuts are applied to select the events over and above the basic generation-level and trigger-level cuts discussed above.
Tables XII-XV show the cut-flow for the signal and the background processes, yielding a fair indication of the efficiency of each cut. We have also calculated the projected significance for each benchmark point at 27 TeV LHC with 15 ab −1 luminosity.
We can see in Tables XII-XV that BP2 fares best in terms of signal significance, followed by BP1, BP3 and BP4. Therefore, BP1 which has the lowest mass of charged Higgs, has lowest discriminating capacity between the signal and background. However, its large crosssection enables us to achieve substantial signal significance. On the other hand, BP4, which 200 GeV |∆η j 1 j 2 | > 3 M j 1 j 2 > 500 GeV 100 GeV < M H + < 300 GeV -300 GeV |∆η j 1 j 2 | > 3 M j 1 j 2 > 500 GeV 100 GeV < M H + < 300 GeV -   corresponds to the heaviest charged Higgs has the best discriminating power among all the BPs. But its meagre production cross-section makes it difficult to obtain a good significance.
However, as we can rely on various good separation variables in case of BP4, we expect significant improvement over the cut-based analysis when we perform multivariate analysis.
Evidently, for BP2 and BP3 there are good optimisations between production cross section and signal-background separation and consequently they yield considerable significance even Signal yield = 24870 Total background yield = 1.8×10 7 Signal significance (0% σ sys_un ) = 5.9 Signal significance (2% (5%) σ sys_un ) = 0.07 (0.03)     Fig.19 in Appendix A. From this figure it can be inferred that the best possible signal background separation occurs for BP4.
The variation of significance with BDT cut value or BDT score, depicted in Fig.6(a) shows that the significance becomes highest for a particular value of BDT score. While computing the final signal significance, we use the BDT cut values for which the significances are maximised for four benchmarks. Fig.6(b) represents the ROC curves for four benchmarks.
it can be clearly seen that the background rejection efficiency for BP4 is best compare to the rest. Inspite of having the best background rejection efficiency, BP4 has the lowest signal significance owing to small signal cross section (Table IX). According to the ability to differentiate between the signal and background, let us choose following thirteen kinematic variables: where η j 1 j 2 and |∆η j 1 j 2 | are the pseudorapidity of the system of two leading light jets and ∆η j 1 j 2 between the two leading jets respectively. ∆φ bb,ljj is the azimuthal angle separation between the two b-jet system and the system of lepton with two leading light jets, and the rest of the variables have their usual meaning. The most important variables to distinguish between the signal and backgrounds are ∆R bb , M j 1 j 2 , η j 1 j 2 , M ef f , |∆η j 1 j 2 |, M bbj 1 , M bbj 2 , M 2blν .
While doing the analysis, we have set F = 0.4 and we fix the branching ratio of the decay H + → tb to be 40%. The normalised distributions of the best kinematic variables for different M H + are shown in Fig. 5. After doing the BDTD analysis, we summarise the results in Table XVI.
Comparing the results obtained from BDTD analysis (Table. XVI) and those from cutbased analysis (Table XII-XV) we can see that BP1 is hardly improved by the BDTD analysis. Since the kinematic variables do not offer a good separation in this case, the BDTD will perform marginally better than an optimized cut-based analysis. The results improve slightly in case of BP2. The improvement becomes striking in case of BP3 and BP4, because these benchmarks correspond to observables with significant discriminating power.
Considering 2% and 5% systematic uncertainties, the significances has also been calculated and tabulated in Table XVI.
With the results of the VBF analyses in hand, we turn our attention to the pp → ZH ± → 5l + / E T channel. To elaborate a bit, the H ± here is produced in association with a Z-boson via an s-channel exchange of the W -boson and it subsequently leads to the hadronically quiet H ± → W ± Z → 3l + / E T decay cascade. The other Z-boson also decays leptonically. The Feynman diagram corresponding to this signal can be found in Fig.1(c). Throughout this analysis, we remain agnostic to the fact that the prospects of observing an H + in this channel could be low and this would be attributed to a small signal cross section(s).
Background processes leading to an exactly 5l + / E T final state are pp → W ± ZZ and pp → hW ± . Subleading backgrounds come from pp → hZ → 6l and pp → tth → 6l+2b+ / E T . Mis-identification of a lepton as missing energy can lead to a 5l + / E T state originating from the former. As for the latter, obtaining the same final state as the signal would entail misidentification of one lepton and 2 b-jets. We show the cross sections for the signal and the backgrounds in Table XVII. BP4 is excluded from the subsequent analysis owing to the small signal cross section. We start by vetoing the b-jets and τ -jets in the final state, so that the final state (comprising of five leptons) is completely leptonic. Here by lepton we mean e, µ and have taken all possible combinations of those five leptons leading to two same flavour opposite sign (SFOS) pairs and one isolated lepton.
Since the MVA was found to outperform the CBA for the two VBF signals, we choose to omit the CBA for this process. Different BDT parameters used to carry out the MVA, are given in Table XVIII. Fig.20 in Appendix A depicts the distribution for signal and backgrounds along with the KS-scores for this channel. In Fig.9(a) and Fig.9(b) we show the variation of significance with BDT score and ROC curves for three benchmarks respectively.
It can be seen that the background rejection efficiency increases with increasing M H + . The Here we reconstruct the two Z-boson invariant masses (M ll ) 1 and (M ll ) 2 by combinatorially identifying the corresponding two SFOS lepton pairs. Any cut on these is expected to suppress the hW ± and tth backgrounds appropriately that do not involve a pair of on-shell Z-bosons. Next, the invariant mass of the two pairs of SFOS lepton system (M llll ), originating from the decay of two Z-bosons, is constructed. Since for the hW ± background, h decays into ZZ * , the distribution of M llll peaks around 125 GeV in Fig.8. Thus a veto cut on M llll around 125 GeV for all values of M H + would help suppressing the aforementioned background to a large extent. Similarly to what was done for the H ± jj → W ± Zjj channel, we construct the invariant mass of the lllν system (M W Z inv ) that comes from H ± . For W ± ZZ background, W ± Z are not originated from the decay of a single mother particle as in the signal. Therefore one can rely on this variable which has substantial discriminatory power between the signal and the backgrounds. Here ∆φ l i / E T is the azimuthal angle between l i and the missing transverse energy vector. Rest of the aforementioned kinematic variables have been defined earlier.
Among these twenty variables, (M ll ) 1 , (M ll ) 2 , M llll , M W Z inv , ∆R l 1 l 2 , ∆R l 1 l 3 , / E T turn out to be the best performing having high discerning ability to separate signal and backgrounds. Corresponding yields for signal and backgrounds at L = 15 ab −1 for M H + = 200 GeV, 300 GeV and 500 GeV along with the signal significances are given in Table XIX. The maximum significance of 4.24 is predicted for M H + = 200 GeV following the BDTD analysis. The significance declines for the higher charged Higgs masses. Significances with 2% and 5% systematic uncertainties are presented in Table XIX.

4.24
Signal significance (2% (5%) σ sys_un ) 4.2 (4.0) Significance (0% σ sys_un ) 1.5 Signal significance (2% (5%) σ sys_un ) 1.5 (1.48) Significance (0% σ sys_un ) 0.7 Signal significance (2% (5%) σ sys_un ) 0.7 (0.7)    Fig. 4, one concludes that for an H + decaying to the massive gauge bosons, ZH + production has a much lower reach compared to VBF production. To quote a few numbers, an H + of mass 500 GeV that corresponds to |F | = 0.4 can be discovered at 5σ for BR(H + → W + Z) 3% via the 3l + 2j + / E T final state. In contrast, the 5l + / E T final state demands F 0.8 and BR(H + → W + Z) 100% to attain the same discovery potential. This subsection contains the analysis for the signal where following production in association with a Z, the charged boson decays to tb. In addition to Z → ll decay, we have the t → blν leptonic decay in this case. This gives rise to three lepton final state with two btagged jets and / E T . The dominant contribution to the background in this final state comes from the ttX type production where X = h, Z, W ± . The ZZZ and W ± ZZ production, Zh, W ± h and bbh are the sub-dominant backgrounds in this channel owing to their smaller cross section compared to the previous ones. We present the cross sections of chosen signal benchmarks and background processes in Table XX.
We generate all of these backgrounds with the application of following generation level  cuts : Further, exactly two b-tagged jets and three isolated leptons with the aforementioned transverse momenta satisfying |η b,l | < 2.5 are demanded. Out of these three leptons we first construct the invariant mass of the SFOS lepton pair, which must peak at M Z . Then the remaining lepton along with the / E T and b-jets is used to reconstruct the charged Higgs mass. We apply the same technique to reconstruct M H + as described in section IV B. Next, we perform the MVA by adjusting the BDT parameters as mentioned in Table XXI. The signal and background distributions and the KS-scores are depicted in Fig.21 in Appendix A. The significance vs. BDT score plot has been presented in Fig.12(a). Degree of rejecting the backgrounds can be identified from the ROC curves for three benchmarks in Fig.12(b).  Table. XXII. In addition, the signal significances considering 2% and 5% systematic uncertainty have been relegated in Table   XXII. Also, We show the normalised distributions of the kinematic variables having best discriminatory power between the signal and the backgrounds in Fig. 11.   Signal 249 Signal significance (0% σ sys_un ) 4.5 Signal significance (2% (5%) σ sys_un ) 3.0 (1.5) GeV and 1 TeV along with signal significances for the 3l +2b+ / E T channel after the BDTD analysis.
40%, we present the contours of constant signal significance in the |F |-BR(H + → tb) plane in Fig. 13. A straightforward comparison with Fig. 7 establishes that for H ± → tb too, the prospects of the ZH ± production process is much weaker compared to the VBF process in searching for an H + . We quote a few numbers to put this into perspective. For F = 0.4, the VBF process offers 5σ observability for BR(H + → tb) 15%. This correspondingly weakens to BR(H + → tb) 70% in case of ZH ± .

V. COMPARING THE VBF SENSITIVITIES
In this section, we compare the sensitivities of the various channels in the |F |-|A tb | plane.
The magenta (sky blue) area in Fig.14 is excluded by the H + jj (tbH ± ) search. Since the pp → jjH ± process statistically outperforms the pp → ZH ± process, we choose it to compare the sensitivities of the H ± → tb, W ± Z decay channels corresponding to the former.
The dashed (solid) red, green and blue lines in Fig.14 respectively correspond to the 2σexclusion, 3σ-discovery and 5σ-discovery contours for the 2b + 1l + 2j + / E T (3l + 2j + / E T ) channel. These contours are drawn by appropriately scaling the results of the multivariate analysis of either channel. An inspection of Fig. 14  plane. This is concurred by an inspection of Fig. 14. An H + of mass 500 GeV is seen to offer a 5σ observability when |A tb | 0.3 and |F | 0.2.
On the other hand, the 3l + 2j + / E T final state comes handy in probing simultaneously low values of |F | and |A tb |. The reason to this can be traced back to the fact that |F | 2 +|A tb | 2 and therefore increasing (decreasing) |F | while keeping the significance fixed requires increasing (decreasing) |A tb | too. This is how the lower left corner of the |F | − |A tb | plane gets accessible, a region the 2b + 1l + 2j + / E T signal is blind to. For example, |F | |A tb | 0.1 for M H + = 200 GeV prospects a 5σ-discovery potential. A complementarity between the two search channels is thus seen.
In order to see how the present analyses for a 27 TeV pp collider fares against the 14 TeV HL-LHC, we extrapolate the 13 TeV, 3 ab −1 analyses of the two VBF channels as reported in [57], to 14 TeV. The extrapolation procedure makes use of the fact that kinematical distributions, and hence, the cut-efficiencies do not change appreciably upon changing from √ s = 13 TeV to 14 TeV. Keeping the integrated luminosity fixed, one can write region is ruled out at 95% confidence level by the pp → H ± jj, H ± → W ± Z (pp → tbH ± , H ± → tb) search. and 27 TeV in Fig.15. The 27 TeV collider at 15 ab −1 is seen to improve the sensitivity of the channel by a great margin. For example, the lowest value of |A tb | the HL-LHC can probe at 5σ is 0.67 and one requires a large |F | = 1 for it. The senstivity for higher M H + is also seen to be much better for 27 TeV.  The sensitivity we obtain in terms of the generic parameters |F | and |A tb | can be straightforwardly translated to a realistic model containing singly charged boson(s) coupling dominantly to (t, b) and (W ± , Z). As discussed before, a prime example is the Georgi-Machacek (GM) model [15] whose scalar sector comprises a real triplet ξ and a complex triplet ∆ over and above the regular doublet. The triplets can be expressed as, Here v ∆ (v ξ ) refers to the vacuum expectation value (VEV) of the triplet ∆(ξ). With v d denoting the VEV of the scalar doublet, one has v 2 d + 2v 2 ∆ + 4v 2 ξ = 246 GeV. A custodial- symmetric scalar potential such as in the GM model demands v ξ = v ∆ / √ 2 in which case there is no constraint on v ∆ from the ρ-parameter.
Counting the singly charged boson coming from the scalar doublet, the total number of singly charged scalars in the GM model becomes three. The corresponding 3×3 mass matrix is diagonalised by the action of a mixing angle α giving rise to a charged goldstone G + , and, the physical scalars H + 3 and H + 5 . Most importantly, both couple to W ± , Z and the third generation quarks for an arbitrary α. More precisely, c W being the cosine of Weinberg angle.
That the H + W − Z interaction is proportional to the triplet VEV is a generic feature of models featuring scalar triplets. The sensitivity contours in the |F |-|A tb | plane can therefore be mapped to the v ∆ -|sinα| plane for the GM model for a given charged boson.   9 . It is worthwhile here to comment on how the choice of α impacts the observability. The cross section of the pp → H ± 5 jj, H ± 5 → tb process is proportional to α . Therefore, the cross section becomes maximum around maximal mixing, i.e, α = π 4 thereby also maximising the observability. We find that a 27 TeV hadronic collider with an integrated luminosity of 15 ab −1 excludes an H + 5 of mass 200 GeV at 95% confidence level for v ∆ 32 GeV in the vicinity of maximal mixing, a significant improvement over an earlier 13 TeV study [63] that reports v ∆ 70 GeV to be the corresponding exclusion limit. The discovery potential also markedly improves upon switching from 13 TeV to 27 TeV. We read from Fig.17  TeV and 36 fb −1 stands at v ∆ 50 GeV, we deem our analysis a considerable improvement over the former and demonstrative of the potential of a 27 TeV pp-collider. A 5σ-discovery is predicted for v ∆ 9 GeV. The sensitivity of this channel for the higher masses remains approximately the same in this corner of the v ∆ -|sinα| plane as can be read from Fig. 17.
Prior to closing this section, we add that the results we obtained in this study can be interpreted in context of any NP scenario that predicts singly charged Higgs bosons dominantly coupling to the third generation quarks and W ± , Z. The non-minimal variants of the GM model can straightforwardly similarly be examined for observability of the singly charged scalars. However, we would like to comment on the observability in case of multi-doublet models where F is radiatively suppressed. As stated in the introduction, the Type-I 2HDM 8 There exist other GM model parameters [25,74] such as scalar couplings, masses and mixing angles that directly do not enter a parameter region shown in the v ∆ − |sinα| plane. It can therefore be inferred that the entire v ∆ − |sinα| plane gets allowed by electroweak precision measurements by a judicious choice of the other model parameters. Hence, the impact of electroweak precision constraints is not separately investigated for this mostly model-independent study. 9 We assume that H + 3 is heavier than H + 5 so that the H + 5 → H + 3 Z decay is disallowed kinematically. and its inert or color octet extensions predict F O(10 −2 ) which is below the sensitivity threshold of the present analysis. That said, there is still the option of arbitarily increasing the number of scalar doublets to enhance F and bring it within the sensitivity reach. However, such models are theoretically contorted and tightly constrained by experimental data.
Overall, it can therefore be concluded that despite such a huge enhancement in sensitivity compared to the HL-LHC, the 27 TeV pp machine with 15 ab −1 integrated luminosity cannot probe the radiative H + W − Z vertex, at least in the minimal multi-doublet models such as 2HDMs and 3HDMs.

VI. SUMMARY AND CONCLUSIONS
In this work, we have investigated the prospects of observing a gauge-phillic singly charged Higgs boson at the 27 TeV upgrade of the LHC attaining 15 ab −1 integrated luminosity.
Motivated by certain realistic scenarios, it is assumed that the charged scalar dominantly interacts with (W ± , Z) and (t, b) via the generic couplings F and A tb respectively. Two kinematically distinct topologies for H ± production, i.e., the VBF process pp → H ± jj and pp → ZH ± are proposed for detailed study. The H ± so produced decays to the tb and W ± Z pairs and complete leptonic cascades of the same are chosen for detailed analyses. That is, the following signals were analysed. (a) pp → H ± jj → W ± Zjj → 3l + 2j + / E T , (b) pp → H ± jj → tbjj → 2b + 1l + 2j + / E T , (c) pp → H ± Z → W ± ZZ → 5l + / E T and (d) pp → H ± Z → tbZ → 3l + 2b + / E T , where l = e, µ. We took M H + = 200 GeV, 300 GeV, 500 GeV and 1 TeV as benchmark mass points for the ensuing analysis. In addition to using the conventional cut-based method, multivariate techniques stemming from the boosted-decision-tree algorithm were also adopted. We summarise below our key results.
• The statistical yield of the VBF production process supersedes that of the pp → ZH ± process for both H ± → W ± Z and H ± → tb.
• The two VBF channels are seen to probe somewhat complementary regions of the |F |-|A tb | plane, which clearly is an upshot of the present analysis. The pp → H ± jj, H ± → W ± Z signal is found to be able to probe the lower left corner of the |F |-A tb plane.
• The generic nature of the study has enabled to interpret the obtained results in context of the Georgi-Machacek model. The relevant parameters for this case are the VEV v ∆ and the mixing angle α. The pp → H ± jj, H ± → tb signal predicts 5σ discovery of a 500 GeV charged scalar for a v ∆ 38 GeV. With the pp → H ± jj, H ± → W ± Z, the 5σ threshold lowers to v ∆ 8 GeV.
As a passing remark, our study highlights that the H + W − Z vertex is an important interaction to search for a singly charged Higgs boson. It also establishes that a √ s = 27 TeV pp collider with an integrated luminosity of 15 ab −1 can be a powerful tool to improve the discovery potential of a charged Higgs via such an interaction over the earlier analyses carried out at lower centre-of-mass energies. The present analysis therefore serves as an important case study on the efficacy of a 27 TeV hadronic collider.