Search for charged Higgs bosons produced in vector boson fusion processes and decaying into vector boson pairs in proton–proton collisions at \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\sqrt{s} = 13\,{\text {TeV}} $$\end{document}s=13TeV

A search for charged Higgs bosons produced in vector boson fusion processes and decaying into vector bosons, using proton–proton collisions at \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\sqrt{s}=13\,{\text {TeV}} $$\end{document}s=13TeV at the LHC, is reported. The data sample corresponds to an integrated luminosity of 137\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\,{\text {fb}}^{-1}$$\end{document}fb-1 collected with the CMS detector. Events are selected by requiring two or three electrons or muons, moderate missing transverse momentum, and two jets with a large rapidity separation and a large dijet mass. No excess of events with respect to the standard model background predictions is observed. Model independent upper limits at 95% confidence level are reported on the product of the cross section and branching fraction for vector boson fusion production of charged Higgs bosons as a function of mass, from 200 to 3000\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\,{\text {GeV}}$$\end{document}GeV. The results are interpreted in the context of the Georgi–Machacek model.

At the LHC, interactions from VBS are characterized by the presence of two gauge bosons in association with two forward jets with a large pseudorapidity separation (| η jj |) and a large dijet invariant mass (m jj ). An excess of events with respect to the SM predictions could indicate the presence of new resonances, such as singly or doubly charged e-mail: cms-publication-committee-chair@cern.ch Higgs bosons. Extended Higgs sectors with additional SU (2) isotriplet scalars give rise to charged Higgs bosons with couplings to W and Z bosons at the tree-level [19]. Specifically, the Georgi-Machacek (GM) model [18,20], with both real and complex triplets, preserves a global symmetry SU(2) L ×SU(2) R , which is broken by the Higgs vacuum expectation value to the diagonal subgroup SU(2) L+R . Thus, the tree-level ratio of the W and Z boson masses is protected against large radiative corrections. In this model, singly (doubly) charged Higgs bosons that decay to W and Z bosons (same-sign W boson pairs) are produced via vector boson fusion (VBF).
The charged Higgs bosons H ± and H ±± in the GM model are degenerate in mass (denoted as m H 5 ) at tree level and transform as a quintuplet under the SU(2) L+R symmetry. The H ± and H ±± bosons are also collectively referred to as H 5 in the context of the GM model. Production and decays of the H 5 states depend on the two parameters m H 5 and s H , where s 2 H characterizes the fraction of the W boson mass squared generated by the vacuum expectation value of the triplet fields. The H 5 states are fermiophobic and are assumed to decay to vector boson pairs with branching fraction of 100% [21]. Figure 1 shows representative Feynman diagrams for the production and decay of the charged Higgs bosons. There are additional charged Higgs bosons H ± predicted in the GM model that transform as a triplet under the SU(2) L+R symmetry. These H ± bosons have only fermionic couplings and are not considered here. This paper presents a search for H ± and H ±± that are produced via VBF and decay to WZ and W ± W ± boson pairs, respectively, using proton-proton (pp) collisions at √ s = 13 TeV. The data sample corresponds to an integrated luminosity of 137 ± 2 fb −1 [22][23][24], collected with the CMS detector [25] in three separate LHC operating periods during 2016, 2017, and 2018. The three data sets are analyzed independently, with appropriate calibrations and corrections, to account for the various LHC running conditions and the performance of the CMS detector. The W ± W ± and WZ channels are simultaneously studied by performing a binned maximum-likelihood fit of distributions sensitive to these processes, following the methods described in Ref. [26]. The searches for H ± and H ±± are performed in the leptonic decay modes W ± Z → ± ν ± ∓ and W ± W ± → ± ν ± ν, where , = e, μ. Candidate events contain either two identified leptons of the same charge or three identified charged leptons with the total charge of ±1, moderate missing transverse momentum ( p miss T ), and two jets with large values of | η jj | and m jj .
Model independent upper limits at 95% confidence level (CL) are reported on the product of the cross section and branching fraction for vector boson fusion production of the H ± and H ±± bosons individually. The results are also interpreted in the context of the GM model including the simultaneous contributions of the H ± and H ±± bosons. Searches for charged Higgs bosons in these topologies have been performed by the CMS Collaboration at 13 TeV using the data sample collected during 2016 [27][28][29]. The ATLAS and CMS Collaborations have also set constraints on the GM model by performing searches for charged Higgs bosons in semileptonic final states at 8 TeV [30] and 13 TeV [31], respectively.

The CMS detector
The central feature of the CMS apparatus is a superconducting solenoid of 6 m internal diameter, providing a magnetic field of 3.8 T. Within the solenoid volume are a silicon pixel and strip tracker, a lead-tungstate crystal electromagnetic calorimeter (ECAL), and a brass and scintillator hadron calorimeter, each composed of a barrel and two endcap sections. Forward calorimeters extend the η coverage provided by the barrel and endcap detectors up to |η| < 5. Muons are detected in gas-ionization chambers embedded in the steel magnetic flux-return yoke outside the solenoid. A more detailed description of the CMS detector, together with a definition of the coordinate system and the relevant kinematic variables, is reported in Ref. [25]. Events of interest are selected using a two-tiered trigger system [32]. The first level, composed of custom hardware processors, uses information from the calorimeters and muon detectors to select events at a rate of around 100 kHz within a fixed latency of 4 μs. The second level, known as the high-level trigger, consists of a farm of processors running a version of the full event reconstruction software optimized for fast processing, and reduces the event rate to around 1 kHz before data storage.

Signal and background simulation
Processes characterized by the presence of two gauge bosons in association with two forward jets are an important background contribution. The processes contributing to diboson plus two jets production that proceeds via the EW interaction are referred to as EW-induced diboson production, leading to tree-level contributions at O(α 4 ), where α is the EW coupling. Figure 2 shows representative Feynman diagrams of EW-induced diboson production involving quartic vertices. An additional contribution to the diboson plus two jets production arises via quantum chromodynamics (QCD) radiation, leading to tree-level contributions at O(α 2 α 2 S ), where α S is the strong coupling. This class of processes is referred to as QCD-induced diboson production. Figure 3 shows representative Feynman diagrams of the QCD-induced production. The associated production of a Z boson and a single top quark, referred to as tZq production, is also an important background contribution. Additional background contributions arise from the tt, tW, ttW, ttZ, ttγ , triple vector boson (VVV, V = W, Z), and double parton scattering processes.
Multiple Monte Carlo (MC) event generators are used to simulate the signal and background contributions. The signal and background processes are produced with on-shell particles. Three sets of simulated events for each process are needed to match the data taking conditions in the three years. The charged Higgs boson signal samples are simulated using MadGraph5_amc@nlo 2.4.2 [33,34] at leading order (LO) accuracy. The predicted signal cross sections are taken at next-to-next-to-LO (NNLO) accuracy from the GM model [21].
The SM EW W ± W ± and WZ processes, where both bosons decay leptonically, are simulated using MadGraph5_amc@nlo at LO accuracy with six EW (O(α 6 )) and zero QCD vertices. The same generator is also Fig. 2 Representative Feynman diagrams of a VBS process contributing to the EW-induced production of events containing W ± W ± (left) and WZ (right) boson pairs decaying to leptons, and two forward jets Fig. 3 Representative Feynman diagrams of the QCD-induced production of W ± W ± (left) and WZ (right) boson pairs decaying to leptons, and two jets used to simulate the QCD-induced W ± W ± process with four EW and two QCD vertices. Contributions with an initialstate b quark are excluded from the EW WZ simulation because they are considered part of the tZq background process. Triboson processes, where the WZ boson pair is accompanied by a third vector boson that decays into jets, are included in the EW WZ simulation. The QCD-induced WZ process is simulated at LO with up to three additional partons in the matrix element calculations using the Mad-Graph5_amc@nlo generator with at least one QCD vertex at tree level. The different jet multiplicities are merged using the MLM scheme [35] to match matrix element and parton shower jets, and the inclusive contribution is normalized to NNLO predictions [36]. The interference between the EW and QCD diagrams is also accounted for with Mad-Graph5_amc@nlo.
A complete set of NLO QCD and EW corrections for the leptonic W ± W ± scattering process has been computed [37,38] and they reduce the LO cross section of the EW W ± W ± process by 10-15%, with the correction increasing in magnitude with increasing dilepton and dijet invariant masses. Similarly, the NLO QCD and EW corrections for the leptonic WZ scattering process have been computed at the orders of O(α S α 6 ) and O(α 7 ) [39], reducing the cross sections for the EW WZ process by 10%. The SM EW W ± W ± and WZ processes are normalized by applying these O(α S α 6 ) and O(α 7 ) corrections to MadGraph5_amc@nlo LO cross sections.
The powheg v2 [40-44] generator is used to simulate the tt, tW, ZZ, and W ± W ∓ processes at NLO accuracy in QCD. Production of ttW, ttZ, ttγ , and VVV events is simulated at NLO accuracy in QCD using the MadGraph5_amc@nlo 2.2.2 (2.4.2) generator for the 2016 (2017 and 2018) samples. The tZq process is simulated in the four-flavor scheme using MadGraph5_amc@nlo 2.3.3 at next-to-LO (NLO). Events in which two hard partonparton interactions occur within a single pp collision, referred to as double parton scattering W ± W ± production, are generated at LO using pythia 8.226 (8. All MC generated events are processed through a simulation of the CMS detector based on Geant4 [52] and are reconstructed with the same algorithms used for data. The simulated samples include additional interactions in the same and neighboring bunch crossings, referred to as pileup. The additional inelastic events are generated using pythia with the same underlying event tune as the main interaction and superimposed on the hard-scattering events. The distribution of the number of pileup interactions in the simulation is adjusted to match the one observed in the data. The average number of interactions per bunch crossing was 23 (32) in 2016 (2017 and 2018) corresponding to an inelastic pp cross-section of 69.2 mb.

Event reconstruction
The primary vertex (PV) is defined as the vertex with the largest value of summed physics-object p 2 T . The physics objects are the jets, clustered using the jet finding algorithm [53,54] with the tracks assigned to candidate vertices as inputs, and the associated missing transverse momentum, taken as the negative vector sum of the p T of those jets.
The CMS particle-flow (PF) algorithm [55] is used to combine the information from the tracker, calorimeters, and muon systems to reconstruct and identify charged and neutral hadrons, photons, muons, and electrons (PF candidates). The missing transverse momentum vector p miss T is defined as the projection onto the plane perpendicular to the beam axis of the negative vector momentum sum of all reconstructed PF candidates in an event. Its magnitude is referred to as p miss T . Jets are reconstructed by clustering PF candidates using the anti-k T algorithm [53] with a distance parameter of 0.4. Additional proton-proton interactions within the same or nearby bunch crossings can contribute additional tracks and calorimetric energy depositions, increasing the apparent jet momentum. To mitigate this effect, tracks identified to be originating from pileup vertices are discarded and an offset correction is applied to correct for remaining contributions [56]. Jet energy corrections are derived from simulation studies so that the average measured energy of jets becomes identical to that of particle level jets. In situ measurements of the momentum balance in dijet, photon+jet, Z+jet, and multijet events are used to determine any residual differences between the jet energy scale in data and in simulation, and appropriate corrections are made [57]. Corrections to jet energies to account for the detector response are propagated to p miss T [58]. Jets with transverse momentum p T > 30 GeV and |η| < 4.7 are included in the analysis.
Events with at least one jet with p T > 20 GeV and |η| < 2.4 that is consistent with the fragmentation of a bottom quark are rejected to reduce the number of top quark background events. The DeepCSV btagging algorithm [59] is used for this selection. For the chosen working point, the efficiency of the algorithm to select bquark jets is about 72% and the rate for incorrectly tagging jets originating from the hadronization of gluons or u, d, s quarks is about 1%. The rate for incorrectly tagging jets originating from the hadronization of c quarks is about 10%.
Events with at least one reconstructed hadronic decay of a τ lepton, denoted as τ h , with p T > 18 GeV and |η| < 2.3, are rejected to reduce the contribution of diboson processes with τ h decays. The τ h decays are reconstructed using the hadrons-plus-strips algorithm [60].
Electrons and muons are reconstructed by associating a track reconstructed in the tracking detectors with either a cluster of energy deposits in the ECAL [61,62] or a track in the muon system [63]. Electrons (muons) must pass loose identification criteria with p T > 10 GeV and |η| < 2.5 (2.

Event selection
Collision events are collected using single-electron and single-muon triggers that require the presence of an isolated lepton with p T > 27 and 24 GeV, respectively [67]. In addition, a set of dilepton triggers with lower p T thresholds is used, ensuring a trigger efficiency above 99% for events that satisfy the subsequent offline selection [67].
Several selection requirements are used to isolate the W ± W ± and WZ topologies defining the signal regions (SRs), while reducing the contributions from background processes [26]. Candidate events must contain exactly two isolated same-sign charged leptons or exactly three isolated charged leptons with p T > 10 GeV, and at least two jets with |η| < 4.7 and the leading jet p j T > 50 GeV. To exclude the selected electrons and muons from the jet sample, the jets are required to be separated from the identified leptons by where φ is the azimuthal angle in radians.
For the WZ candidate events, one of the oppositely charged same-flavor leptons from the Z boson candidate is required to have p T > 25 GeV and the other p T > 10 GeV with the invariant mass of the dilepton pair m satisfying |m − m Z | < 15 GeV. For candidate events with three same-flavor leptons, the oppositely charged lepton pair with the invariant mass closest to the world-average Z boson mass m Z [68] is selected as the Z boson candidate. The third lepton associated with the W boson is required to have p T > 20 GeV. In addition, the trilepton invariant mass m is required to exceed 100 GeV to exclude a region where production of Z bosons with final-state photon radiation is expected to contribute.
One of the leptons in the same-sign W ± W ± candidate events is required to have p T > 25 GeV and the other p T > 20 GeV. The value of m must be greater than 20 GeV. Candidate events in the dielectron final state with |m − m Z | < 15 GeV are rejected to reduce the number of Z boson background events where the sign of one of the electron candidates is misidentified.
The VBF topology is targeted by requiring the two highest p T jets to have a mass m jj > 500 GeV and a pseudorapidity separation | η jj | > 2.5. The W and Z bosons in the VBF topologies are mostly produced in the central rapidity region with respect to the two selected jets. The candidate W ± W ± (WZ) events are required to satisfy max(z for one of the selected leptons. Here η is the pseudorapidity of the lepton, and η j 1 and η j 2 are the pseudorapidities of the two candidates VBF jets.
The p miss T is required to exceed 30 GeV for both SRs. The selection requirements used to define the same-sign W ± W ± and WZ SRs are summarized in Table 1.

Background estimation
A combination of methods based on simulation and on control regions (CRs) in data is used to estimate background contributions. By inverting some of the requirements in Table 1 we select background-enriched CRs. Uncertainties related to the theoretical and experimental predictions are estimated as described in Sect. 8.
The nonprompt lepton backgrounds originating from leptonic decays of heavy quarks, hadrons misidentified as leptons, and electrons from photon conversions are suppressed by the identification and isolation requirements imposed on leptons. The remaining contribution from the nonprompt lepton background is dominant in the W ± W ± SR and is estimated directly from data following the technique described in Ref.
[70], using events selected by the final selection criteria, except for one of the leptons, which is requested to pass a looser criterion having failed the nominal selection.
The yield in this sample is extrapolated to the signal region using the efficiencies for such loosely identified leptons to pass the standard lepton selection criteria. This efficiency is calculated in a sample of events dominated by dijet production. An uncertainty of 20% is assigned for the nonprompt lepton background normalization to include possible differences in the composition of jets between the data sample used to derive these efficiencies and the data samples in the W ± W ± and WZ SRs [64]. The background contribution from the electron sign mismeasurement is estimated from the simulation by applying a data-to-simulation efficiency correction due to electrons with sign mismeasurement. These corrections are determined using Z → ee events in the Z boson peak region that were recorded with independent triggers. These corrections amount to 40% for data collected in 2017 and 2018, while they are negligible for 2016 data. The electron sign mismeasurement rate is about 0.01 (0.3)% in the barrel (endcap) region [61,62].
Three CRs are used to select nonprompt lepton, tZq, and ZZ background-enriched events to further estimate the normalization of these background processes from data. The nonprompt lepton CR is defined by requiring the same selection as for the W ± W ± SR, but with the bjet veto requirement inverted. The selected events are enriched in the nonprompt lepton background coming mostly from semileptonic tt events. Similarly, the tZq CR is defined by requiring the same selection as for the WZ SR, but with the bjet veto requirement inverted. The selected events are dominated by the tZq background process. Finally, the ZZ CR selects events with two opposite-sign same-flavor lepton pairs with the same VBS-like requirements. The three CRs are used The m jj distributions after requiring the same selection as for the WW (upper) and WZ (lower) SRs, but with a requirement of 200 < m jj < 500 GeV. The predicted yields are shown with their best fit normalizations from the simultaneous fit (described in Sect. 7) for the background-only hypothesis i.e., assuming no contributions from the H ± and H ±± processes. Vertical bars on data points represent the statistical uncertainty in the data. The histograms for tVx backgrounds include the contributions from ttV and tZq processes. The histograms for other backgrounds include the contributions from double parton scattering, VVV, and from oppositely charged dilepton final states from tt, tW, W + W − , and Drell-Yan processes. The overflow is included in the last bin. The lower panels show the ratio of the number of events observed in data to that of the total SM prediction. The hatched gray bands represent the uncertainties in the predicted yields. The solid lines show the signal predictions for values of s H = 1.0 and m H 5 = 500 GeV in the GM model together with the SRs to constrain the normalization of the nonprompt lepton, tZq, and ZZ background processes from data. All other background processes are estimated from simulation after applying corrections to account for the small differences between data and simulation. The shapes of the tZq and ZZ background processes are estimated from simulation as well.
The prediction for the QCD WZ background process is validated in a CR defined by requiring the same selection as for the WZ SR, but with a requirement of 200 < m jj < 500 GeV. The predicted yields are shown with their best fit normalizations from the simultaneous fit (described in Sect. 7) for the background-only hypothesis i.e., assuming no contributions from the H ± and H ±± processes. Good agreement between the data and post-fit predicted yields is observed in this CR as can be seen in Fig. 4.

Signal extraction
A binned maximum-likelihood fit is performed using the W ± W ± and WZ SRs, and the nonprompt lepton, tZq, and ZZ CRs to discriminate between the signal and the remaining backgrounds. Signal contributions with electrons and muons produced in the decay of a τ lepton are included. The normalization factors for the tZq and ZZ background processes, affecting both the SRs and CRs, are included as free param-eters in the maximum-likelihood fit together with the signal strength. The SM W ± W ± (WZ) contribution is obtained from the sum of the EW W ± W ± (WZ), QCD W ± W ± (WZ), and the interference contributions according to the SM predictions [26] and allowed to vary within the uncertainties.
The diboson transverse mass (m VV T ) is constructed from the four-momentum of the selected charged leptons and the p miss T . The four-momentum of the neutrino system is defined using the p miss T , assuming that the values of the longitudinal component of the momentum and the mass are zero. The value of m VV T , defined as where E i and p z,i are the energies and longitudinal components of the momenta of the leptons and neutrino system from the decay of the gauge bosons in the event, is effective in discriminating between the resonant signal and nonresonant background processes. The value of m jj is effective in discriminating between all non-VBS processes and the signal (plus EW VV) processes because VBF and VBS topologies typically exhibit large values for the dijet mass. A profile likelihood technique is used where systematic uncertainties are represented by nuisance parameters [71]. For each individual bin, a Poisson likelihood term describes the fluctuation of the data around the expected central value, which is given by the sum of the contributions from signal and background processes. The systematic uncertainties are treated as nuisance parameters and are profiled with the shape and normalization of each distribution varying within the respective uncertainties in the fit. The normalization uncertainties are treated as log-normal nuisance parameters. Correlation across bins is taken into account. The uncertainties affecting the shapes of the distributions are modeled in the fit as nuisance parameters with external Gaussian constraints. The dominant nuisance parameters are not significantly constrained by the data, i.e., the normalized nuisance parameter uncertainties are close to unity.

Systematic uncertainties
Several sources of systematic uncertainty are taken into account in the signal extraction procedure. For each source of uncertainty, the effects on the signal and background distributions are considered to be correlated.
The total Run 2 (2016-2018) integrated luminosity has an uncertainty of 1.8%, the improvement in precision relative to Refs. [22-24] reflecting the (uncorrelated) time evolution of some systematic effects.
The simulation of pileup events assumes an inelastic pp cross section of 69.2 mb, with an associated uncertainty of 5% [72], which has an impact on the expected signal and background yields of about 1%.
Discrepancies in the lepton reconstruction and identification efficiencies between data and simulation are corrected by applying scale factors to all simulation samples. These scale factors, which depend on the p T and η for both electrons and muons, are determined using Z → events in the Z boson peak region that were recorded with independent triggers [61,63,73]. The uncertainty in the determination of the trigger efficiency leads to an uncertainty smaller than 1% in the expected signal yield. The trigger efficiency in the simulation is corrected to account for the effect of a gradual time shifts in the forward region in the ECAL endcaps for the 2016 and 2017 data [74]. The uncertainty in this correction is included in the trigger efficiency uncertainty. The lepton momentum scale uncertainty is computed by varying the lepton momenta in simulation with their uncertainties, and repeating the analysis selection. The resulting uncertain-ties in the yields are ≈1% for both electrons and muons. These uncertainties are assumed to be correlated across the three data sets.
The uncertainty in the calibration of the jet energy scale (JES) directly affects the acceptance of the jet multiplicity requirement and the p miss T measurement. These effects are estimated by shifting the JES in the simulated samples up and down by one standard deviation. The uncertainty in the jet energy resolution (JER) smearing applied to simulated samples to match the p T resolution measured in data causes both a change in the normalization and in the shape of the distributions. The overall uncertainty in the JES and JER is 2-5%, depending on p T and η [57,75], and its impact on the expected signal and background yields is about 3%.
The b tagging efficiency in the simulation is corrected using scale factors determined from data [59]. These values are estimated separately for correctly and incorrectly tagged jets. Each set of values results in uncertainties in the b tagging efficiency of about 1-4% depending on p T and η, and the impact on the expected signal and background yields is about 1%. The uncertainties in the JER, JES and b tagging are treated as uncorrelated across the three data taking years, since the detector conditions have changed among the three years.
The theoretical uncertainties associated with the choice of the renormalization and factorization scales are estimated     78]. The uncertainty in the yields due to missing higher-order EW corrections in the GM model is estimated to be 7% [21]. These theoretical uncertainties may affect both the estimated signal and background rates. The statistical uncertainties that are associated with the limited number of simulated events and data events used to estimate the nonprompt lepton background are also considered as systematic uncertainties.
A summary of the impact of the systematic uncertainties on the signal strength, μ, defined as the ratio of the observed charged Higgs signal yield to the expected yield,  Table 2 for the case of a background-only simulated data set, i.e., assuming no contributions from the H ± and H ±± processes. Table 2 also shows systematic uncertainties including a charged Higgs boson signal for values of s H = 1.0 and m H 5 = 500 GeV in the GM model. The impacts shown in Table 2 result from a fit to two simulated samples: background-only (first column, expected μ = 0) and signalplus-background (second column, expected μ = 1). They differ from the impacts in percent on the expected signal and background yields given above, which are estimated before the fit. The total systematic uncertainty is smaller for the background-only simulated data set because the uncertainties partially cancel out between the SRs and the CRs for the background processes.

Results
The distributions of m jj and m VV T in the WW and WZ SRs are shown in Fig. 5. The m jj distributions in the WW and WZ SRs are shown with finer binning compared to the binning used in the two-dimensional distribution in the fit. Distributions for signal, backgrounds, and data for the bins used in the simultaneous fit are shown in Fig. 6. The data yields, together with the background expectations with the best fit normalizations for the background-only hypothesis, i.e., assuming no contributions from the H ± and H ±± processes, are shown in Table 3. The product of kinematic acceptance and selection efficiency within the fiducial region for the H ±± → W ± W ± → 2 2ν and H ± → WZ → 3 ν processes, as a function of m H 5 , is shown in Fig. 7. The drop of selection efficiency for the H ± → WZ → 3 ν process for masses above 1000 GeV is coming from the lepton isolation requirement as the leptons from high-momentum Z boson decay are produced with a small angular separation.  . 7 The product of acceptance and selection efficiency within the fiducial region for the VBF H ±± → W ± W ± → 2 2ν and H ± → WZ → 3 ν processes, as a function of m H 5 . The combination of the statistical and systematic uncertainties is shown. The theoretical uncertainties in the acceptance are also included No significant excess of events above the expectation from the SM background predictions is found. The 95% CL upper limits on the charged Higgs production cross sections are calculated using the modified frequentist approach with the CL s criterion [79,80] and asymptotic method for the test statistic [71,81].
Constraints on resonant charged Higgs boson production are derived. The exclusion limits on the product of the doubly charged Higgs boson cross section and branching fraction σ VBF (H ±± ) B(H ±± → W ± W ± ) at 95% CLas a function of m H ±± are shown in Fig. 8 (upper left). The exclusion limits on the product of the charged Higgs boson cross section and branching fraction σ VBF (H ± ) B(H ± → WZ) at 95% CLas a function of m H ± are shown in Fig. 8 Fig. 8 (lower). The blue shaded region shows the parameter space for which the H 5 total width exceeds 10% of m(H 5 ), where the model is not applicable because of perturbativity and vacuum stability requirements [21]. For the probed parameter space and m VV T distribution used for signal extraction, the varying width as a function of s H is assumed to have negligible effect on the result. The observed limit excludes s H values greater than 0.20-0.35 for the m H 5 range from 200 to 1500 GeV. The limit improves the sensitivity of the previous CMS results at 13

Summary
A search for charged Higgs bosons produced in vector boson fusion processes and decaying into vector bosons, using proton-proton collisions at √ s = 13 TeV at the LHC, is reported. The data sample corresponds to an integrated luminosity of 137 fb −1 , collected with the CMS detector between 2016 and 2018. The search is performed in the leptonic decay modes W ± W ± → ± ν ± ν and W ± Z → ± ν ± ∓ , where , = e, μ. The W ± W ± and WZ channels are simultaneously studied by performing a binned maximum-likelihood fit using the transverse mass m T and dijet invariant mass m jj distributions. No excess of events with respect to the standard model background predictions is observed. Model independent upper limits at 95% confidence level are reported on the product of the cross section and branching fraction for vector boson fusion production of charged Higgs bosons decaying into vector bosons as a function of mass from 200 to 3000 GeV. The results are interpreted in the Georgi-Machacek (GM) model for which the most stringent limits to date are derived. The observed 95% confidence level limits exclude GM s H parameter values greater than 0.20-0.35 for the mass range from 200 to 1500 GeV. accelerator departments for the excellent performance of the LHC and thank the technical and administrative staffs at CERN and at other CMS institutes for their contributions to the success of the CMS effort. In addition, we gratefully acknowledge the computing centers and personnel of the Worldwide LHC Computing Grid and other centres for delivering so effectively the computing infrastructure essential to our analyses. Finally, we acknowledge the enduring support for the construction and operation of the LHC, the CMS detector, and the supporting computing infrastructure provided by the following funding agencies:

Data Availability Statement
This manuscript has no associated data or the data will not be deposited. [Authors' comment: Release and preservation of data used by the CMS Collaboration as the basis for publications is guided by the CMS policy as written in its document "CMS data preservation, re-use and open access policy" (https://cms-docdb.cern. ch/cgi-bin/PublicDocDB/RetrieveFile?docid=6032&filename= CMSDataPolicyV1.2.pdf&version=2).] Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permit-ted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecomm ons.org/licenses/by/4.0/. Funded by SCOAP 3 .