Search for doubly and singly charged Higgs bosons decaying into vector bosons in multi-lepton final states with the ATLAS detector using proton-proton collisions at s\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ \sqrt{\mathrm{s}} $$\end{document} = 13 TeV

A search for charged Higgs bosons decaying into W±W± or W±Z bosons is performed, involving experimental signatures with two leptons of the same charge, or three or four leptons with a variety of charge combinations, missing transverse momentum and jets. A data sample of proton-proton collisions at a centre-of-mass energy of 13 TeV recorded with the ATLAS detector at the Large Hadron Collider between 2015 and 2018 is used. The data correspond to a total integrated luminosity of 139 fb−1. The search is guided by a type-II seesaw model that extends the scalar sector of the Standard Model with a scalar triplet, leading to a phenomenology that includes doubly and singly charged Higgs bosons. Two scenarios are explored, corresponding to the pair production of doubly charged H±± bosons, or the associated production of a doubly charged H±± boson and a singly charged H± boson. No significant deviations from the Standard Model predictions are observed. H±± bosons are excluded at 95% confidence level up to 350 GeV and 230 GeV for the pair and associated production modes, respectively.


Introduction
Experimental signatures with two leptons with the same electric charge (same-charge) or multi-lepton final states are extensively exploited in searches for physics beyond the Standard Model (BSM physics) at hadron colliders. In many models, heavy BSM particles may be produced in proton-proton collisions and decay into multiple massive Standard Model (SM) electroweak gauge bosons or top quarks. Subsequent decays of these gauge bosons into final states with leptons can occur with considerable branching ratios. These final states are favourable to the search for new phenomena since the yields predicted within the SM are generally low and the experimental effects are well understood. At the Large Hadron Collider (LHC) [1], signatures with two leptons of the same charge, or three or four leptons with a variety of charge combinations have been used by the ATLAS and CMS experiments to explore the landscape of possible SM extensions and their phenomenology [2][3][4][5][6][7][8][9][10][11][12][13][14][15][16][17][18][19][20][21]. Among these proposed extensions, the addition to the SM sector of one weak gauge triplet of scalar fields with a non-zero vacuum expectation value of the neutral component is a compelling way to account for neutrino masses through the type-II seesaw mechanism [22]. is achieved by requiring the neutral components of the SM scalar doublet H and scalar triplet ∆ to acquire vacuum expectation values, v d and v t respectively. After the EWSB, the mixing between these fields results in seven scalar bosons: H ±± , H ± , A 0 (CP odd), H 0 (CP even) and h 0 (CP even). The h 0 boson is the SM Higgs boson, and a small mixing between the CP-even scalars is naturally obtained in the model. In addition, the triplet-neutrino Yukawa term provides non-zero neutrino masses proportional to the vacuum expectation value of the triplet v t . In this type-II seesaw model, constraints from electroweak precision measurements lead to an upper bound on v t of approximately 2.5 GeV, which is significantly lower than the electroweak scale and matches the need for small values suggested by the natural association of v t with the neutrino masses. Additional, theoretical constraints are the absence of tachyonic modes, requirements for the vacuum structure and stability of the potential, and unitarity requirements. A detailed discussion can be found in ref. [29]. The theoretical and experimental constraints help to choose allowed values for the parameters of the model from which the predominant production and decay modes follow. Some type-II seesaw models also extend the weak sector by introducing two scalar triplets and righthanded gauge bosons W R (left-right symmetric models). The model studied in this analysis corresponds to the left-handed triplet discussed in refs. [30][31][32][33][34].
Two scenarios with v t equal to 100 MeV are explored in this paper. In the first scenario, the mass difference between H ±± and H ± bosons is > 100 GeV, the H ± boson being the heavier one. Large mass splitting could still be allowed despite the |m H ± −m H ±± | < 40 GeV constraint deduced from the electroweak precision data [35] if radiative corrections reduce the tree-level contribution coming from the triplet vacuum expectation value. Only the pair production of H ±± bosons via the diagram shown in figure 1a (pp → γ * /Z * → H ±± H ∓∓ ) is considered. The triplet vacuum expectation value is taken to be v t = 100 MeV such that only the doubly charged H ±± boson decay into a pair of W bosons with the same charge, H ±± → W ± W ± , is possible [28]. The leptonic decays H ±± → ± ± are suppressed with increasing v t [36, 37] and not considered in this paper.
In the second scenario, the mass of the H ± boson is chosen to be at most 5 GeV different from the H ±± boson mass. The low mass difference ensures that the charged Higgs bosons do not decay into each other with a non-negligible branching fraction, and complies with -3 -

JHEP06(2021)146
the additional |m H ± − m H ±± | < 40 GeV constraint deduced from the electroweak precision data. In addition, the chosen mass difference between the doubly and singly charged Higgs bosons maximises the signal amplitudes. Only the associated production of H ±± and H ± bosons via the diagram shown in figure 1b (pp → W ± * → H ±± H ∓ ) is considered. The production cross sections for H ±± pair production with the same mass settings as for associated production of H ±± and H ± bosons can be large. However, this production mode is exactly the one described in the first scenario and it is not considered in the second scenario. This choice is motivated by the objective of the search, which is to study a characteristic signature with a benchmark parameter point and report its cross section. The H ±± boson decays into a pair of W bosons with the same charge, with a branching ratio of 100%. Only the bosonic decays of the singly charged bosons (H ± → W ± Z) are considered and, depending on the H ± boson mass, the branching ratio varies between 40% and 60% (see table 1). The branching fraction for the m H ±± = 300 GeVmass hypothesis is lower than for the neighbouring points because of the high dependency of this quantity on the mass of the scalar triplet ∆ [29]. Depending on the H ± boson mass, other decay modes are H ± → tb and H ± → W ± h 0 [28]. Studies at Monte Carlo generator level show that after the selection of at least two same-charge leptons or at least three leptons and no b-jet (where the b-jets are selected with 70% efficiency), the contribution from the other possible H ± decays is negligible. Similar conclusions are also reached after examining the results obtained for the various control regions used in the analysis. The effect on the mass below which a charged Higgs boson is excluded is negligible.
Pair production of H ± bosons is also possible, albeit with a much smaller production cross section than for H ±± pair production. Therefore, H ± pair production is not considered in this paper. The cross section for single H ±± production via vector-boson fusion (pp → W ± * W ± * → H ±± ) is proportional to v t , and hence negligible.
For the H ±± pair production mode, the mixing between the CP-even scalars is taken to be 10 −4 and the remaining five couplings in the potential are adjusted to obtain a given H ±± boson mass hypothesis while requiring h 0 to have a mass of 125 GeV. Similar settings are also used for the associated production mode. The branching fraction times crosssection calculation for the pair production of H ±± bosons and the associated production of H ±± and H ± bosons is performed for on-shell W and Z bosons, and therefore only the region m H ±± > 200 GeVis considered in the present analysis.
Extensive searches for leptonic decays have been performed at various experiments [38][39][40][41][42][43], excluding doubly charged H ±± bosons with masses up to about 870 GeV. The CMS collaboration published results for the H ±± → W ± W ± decay mode in the context of single H ±± production through vector-boson fusion with v t values of a few tens of GeV, for a model with two Higgs triplets [44][45][46]. In contrast, the H ±± → W ± W ± decay mode is not often investigated for v t values around 100 MeV. In this case, the difference is that single H ±± production is suppressed and only H ±± pair production is sizeable. A direct search for H ±± pair production with decays to W ± W ± pairs has been performed on a smaller data set by the ATLAS collaboration [23], validating the principle of such an approach.

ATLAS detector
The ATLAS experiment [47] at the LHC is a multipurpose particle detector with a forwardbackward symmetric cylindrical geometry and nearly 4π coverage in solid angle. 1 It consists of an inner tracking detector (ID) surrounded by a thin superconducting solenoid providing a 2 T axial magnetic field, electromagnetic (EM) and hadron calorimeters, and a muon spectrometer (MS). The inner tracking detector covers the pseudorapidity range |η| < 2.5. It consists of silicon pixel, silicon microstrip and transition radiation tracking detectors; the innermost layer is 33 mm from the beamline [48][49][50]. Lead/liquid-argon (LAr) sampling calorimeters provide electromagnetic energy measurements with high granularity. A steel/scintillator-tile hadron calorimeter covers the central pseudorapidity range (|η| < 1.7). The endcap and forward regions are instrumented with LAr calorimeters for EM and hadronic energy measurements up to |η| = 4.9. The muon spectrometer surrounds the calorimeters and is based on three large air-core toroidal superconducting magnets with eight coils each. The field integral of the toroids ranges between 2.0 and 6.0 T·m across most of the detector. The muon spectrometer includes a system of precision tracking chambers and fast detectors for triggering. A two-level trigger system [51] is used to select events. The first-level trigger is implemented in hardware and uses a subset of the detector information to keep the accepted rate below 100 kHz. This is followed by a software-based trigger that reduces the accepted event rate to 1 kHz on average depending on the data-taking conditions.

Monte Carlo event simulation
Monte Carlo (MC) event generators were used to simulate the signal and background events produced in the proton-proton collisions. For the H ±± pair production and the H ±± and H ± associated production signal processes, the events at particle level were generated with MadGraph5_aMC@NLO [53] using leading-order (LO) matrix elements (ME) and NNPDF3.0lo parton distribution functions (PDF) [54]. The events were subsequently showered using Pythia 8 [55] with a set of tuned parameters called the A14 tune [56]. Signal processes were simulated for different mass hypotheses between 200 and 600 GeV, and presented in table 1. For the pair and associated production, the cross section decreases rapidly with the charged Higgs boson's mass as shown in table 1. The event samples were normalised using calculations at next-to-leading order (NLO) in perturbative QCD [52]. The NLO K-factor increases the expected event yields by a factor of 1.25, independently of the mass of the charged Higgs bosons. 1 ATLAS uses a right-handed coordinate system with its origin at the nominal interaction point (IP) in the centre of the detector and the z-axis along the beam pipe. The x-axis points from the IP to the centre of the LHC ring, and the y-axis points upwards. Cylindrical coordinates (r, φ) are used in the transverse plane, φ being the azimuthal angle around the z-axis. The pseudorapidity is defined in terms of the polar angle θ as η ≡ − ln tan(θ/2), and is equal to the rapidity y ≡ 0.5 ln ((E + pz)/(E − pz)) in the relativistic limit. Angular distance is measured in units of ∆R ≡ (∆y) 2 + (∆φ) 2 . The magnitude of the momentum in the plane transverse to the beam axis is denoted by p T .  associated production modes at next-to-leading order in QCD [52]. The branching ratios (B) of the charged Higgs bosons to W ± W ± or W ± Z are included in the quoted values.

Process
Generator ME accuracy PDF Parton shower and hadronisation Parameter set  Table 2. Summary of the event generators, parton shower models, and PDF sets used for the simulation of the background event samples. The notation V is used to refer to an electroweak gauge boson W or Z/γ * . In the final column, 'default' refers the to default parameter set provided with the event generator.
A summary of the MC event generators used to simulate the background is presented in table 2 and further details are given below. The notation V is used to refer to an electroweak gauge boson W or Z/γ * .  [57] generator depending on the process, including offshell effects and Higgs boson contributions, where appropriate. The V V and V γ processes were simulated using matrix elements at NLO accuracy in QCD for up to one additional parton and at LO accuracy for up to three additional parton emissions. The electroweak production of dibosons in association with two jets (V V -EW jj) was generated using LO accuracy matrix elements. The matrix element calculations were matched and merged -6 -

JHEP06(2021)146
with the Sherpa parton shower based on Catani-Seymour dipole factorisation [58,59] using the MEPS@NLO prescription [60][61][62][63]. The virtual QCD corrections were provided by the OpenLoops library [64,65]. The production of triboson (V V V ) events was simulated with the Sherpa v2.2.2 generator using factorised gauge-boson decays. All diboson and triboson processes were generated using the NNPDF3.0nnlo PDF set [54], along with the dedicated set of tuned parton-shower parameters developed by the Sherpa authors. The cross sections from the event generator were used for the normalisation. The production of an electroweak gauge boson or virtual photon in association with jets (V +jets) was simulated with the Sherpa v2.2.1 generator using NLO matrix elements for up to two partons, and LO matrix elements for up to four partons calculated with the Comix [58] and OpenLoops libraries. The event samples were generated using the NNPDF3.0nnlo PDF set, along with the dedicated set of tuned parton-shower parameters developed by the Sherpa authors. The V +jets event samples were normalised to a NNLO prediction.
Higgs boson production in association with a vector boson was simulated at LO with Pythia 8.186 and EvtGen [66] using the A14 tune, along with the NNPDF2.3lo PDF set. The Monte Carlo prediction was normalised to cross sections calculated at NNLO in QCD with NLO electroweak corrections for qq/qg → V H, and at next-to-leading-logarithm accuracy in QCD with NLO EW corrections for gg → ZH [67][68][69][70][71][72][73]. The production of ttH events was modelled using the Powheg Box v2 [74][75][76][77][78] generator at NLO with the NNPDF3.0nlo PDF set [54]. The events were interfaced to Pythia 8.230 using the A14 tune and the NNPDF2.3lo PDF set [54]. The decays of bottom and charm hadrons were performed by EvtGen v1.6.0. The production of ttV , tW Z and tZ events was modelled using the MadGraph5_aMC@NLO v2.3.3 generator at NLO with the NNPDF3.0nlo PDF. The events were interfaced to Pythia 8.210 (Pythia 8.212 for tW Z) using the A14 tune and the NNPDF2.3lo PDF set. The decays of bottom and charm hadrons were simulated using the EvtGen v1.2.0 program. The simulated ttH, ttV , tW Z and tZ events were normalised to cross-section calculations at NLO accuracy in QCD. In addition, NLO EW corrections were included for ttH, ttV and tW Z. Backgrounds from top-quark pair production (tt) and tW production were estimated at NLO accuracy in QCD using the hvq program [74] in Powheg Box v2. The event samples were generated using the A14 tune and NNPDF3.0nnlo PDF set. The interference between tt and tW production is neglected as it has a negligible impact on the analysis. The tt event sample was normalised to the cross-section prediction at NNLO in QCD including the resummation of next-to-next-to-leading logarithmic (NNLL) soft-gluon terms calculated using Top++2.0 [79][80][81][82][83][84][85]. The inclusive tW cross section was corrected to the theory prediction calculated at NLO in QCD with NNLL soft gluon corrections [86,87].
The production of tttt, ttt, ttW W and ttW Z events was modelled using the Mad-Graph5_aMC@NLO v2.3.3 generator at NLO with the NNPDF3.1nlo PDF [54]. The events were interfaced with Pythia 8.230 using the A14 tune and the NNPDF2.3lo PDF set. The tttt, ttW W and ttW Z contributions were normalised to theoretical cross sections at NLO in QCD including electroweak corrections [88]. The ttt production was normalised to a cross section calculated to LO in QCD.

JHEP06(2021)146
Backgrounds such as pp → W + W + W − W − or pp → W ± W ± ZW ∓ with 2 sc , 3 or 4 final states, which have the same signature as the considered signal, have very small production cross sections times branching ratio and their contribution in the preselection, control and signal regions was found to be negligible.
The signal and background events were passed through the GEANT4 [89] simulation of the ATLAS detector [90] and reconstructed using the same algorithms as are used for the data. The effect of multiple proton-proton interactions in the same or nearby bunch crossings (pile-up) is accounted for using inelastic proton-proton interactions generated by Pythia 8 [91], with the A3 tune [92] and the NNPDF2.3lo PDF set [93]. These inelastic proton-proton interactions were added to the signal and background event samples and weighted such that the distribution of the average number of proton-proton interactions in simulation matches that observed in the data.

Event reconstruction
The analysis is performed in pp collision data recorded by the ATLAS detector between 2015 and 2018. In this period, the LHC delivered colliding beams with a peak instantaneous luminosity up to L = 2.1 × 10 34 cm −2 s −1 , achieved in 2018, and an average number of pp interactions per bunch crossing of 34. After requirements on the stability of the beams, the operational status of all ATLAS detector components, and the quality of the recorded data, the total integrated luminosity of the data set corresponds to 139 fb −1 [94].
Proton-proton interaction vertices are reconstructed from charged-particle tracks with p T > 500 MeV [95,96] in the ID. The presence of at least one such vertex with a minimum of two associated tracks is required, and the vertex with the largest sum of p 2 T of associated tracks is chosen as the primary vertex.
The anti-k t algorithm [97,98] with radius parameter R = 0.4 is used to reconstruct jets up to |η| = 4.9. It uses as inputs particle-flow objects, combining tracking and calorimetric information [99]. The jets are calibrated as described in ref. [100]. Only jets with p T > 20 GeV and |η| < 2.5 are considered further. Events are removed if they contain jets induced by calorimeter noise or non-collision background, according to criteria similar to those described in ref. [101]. Jets with p T < 60 GeV and |η| < 2.4 from pile-up interactions are suppressed using a jet-vertex tagging discriminant [102].
Jets containing b-flavored hadrons are identified in the region |η| < 2.5 with the DL1r b-tagging algorithm based on a recurrent neural network [103,104]. It makes use of the impact parameters of tracks associated with the jet, the position of reconstructed secondary vertices and their compatibility with the decay chains of such hadrons. The b-tagging average efficiency is 70%, as measured in tt events.
Electron candidates are reconstructed as tracks in the ID matched to energy clusters in the EM calorimeter, within |η| < 2.47 [105]. Only electrons with p T > 10 GeV and not in the transition region between the barrel and endcap calorimeters (1.37 < |η| < 1.52) are considered. The electron identification is based on a multivariate likelihood-based discriminant that uses the shower shapes in the EM calorimeter and the associated track properties measured in the ID. The electron candidates must satisfy the 'Loose' identification criteria described in ref. [105]. Signal electrons are required to satisfy the 'Tight' -8 -

JHEP06(2021)146
identification [105] for better rejection of non-prompt electrons. The electron identification efficiency varies with increasing p T in Z → ee events, from 65% at p T = 10 GeV to 88% at p T = 100 GeV for the Tight operating point, and from 85% at p T = 20 GeV to 95% at p T = 100 GeV for the Loose operating point.
The longitudinal impact parameter of the electron track, z 0 , is required to satisfy |z 0 sin θ| < 0.5 mm, where θ is the polar angle of the track. The transverse impact parameter divided by its uncertainty, |d 0 |/σ(d 0 ), is required to be less than five. For all signal electrons there must be no association with a vertex from a reconstructed photon conversion [105] in the detector material. To further reduce the photon conversion background in the 2 sc channel, additional requirements are applied to the signal electrons [4]: i) the candidate must not have a reconstructed displaced vertex with radius r > 20 mm whose reconstruction uses the track associated with the electron, ii) the invariant mass of the system formed by the track associated with the electron and the closest track at the primary or a conversion vertex is required to be larger than 100 MeV; this selection is referred to as the photon conversion veto.
For the signal electrons, the identification criteria are complemented by an isolation requirement, which is based on the energy in a cone around the electron candidate calculated using either reconstructed tracks or energy clusters. Electrons with wrongly reconstructed charge (charge-flip) are suppressed using a boosted decision tree discriminant exploiting additional tracks in the vicinity of the electron and track-to-cluster matching variables [105].
Muon candidates are reconstructed in the region |η| < 2.5 from MS tracks matching ID tracks. The analysis only considers muons with p T > 10 GeV satisfying the 'Medium' quality requirements defined in ref.
[106]. The muon reconstruction efficiency is approximately 98% in Z → µµ events. The longitudinal impact parameter of the muon track must satisfy |z 0 sin θ| < 0.5 mm and the transverse impact parameter must satisfy |d 0 |/σ(d 0 ) < 3. For signal muons the candidate must satisfy calorimeter-and track-based isolation requirements.
Non-prompt electrons and muons from the decays of b-and c-flavored hadrons are further rejected using a boosted decision tree discriminant based on isolation and secondary vertex information, referred to as the non-prompt-lepton veto [107].
To avoid cases where the detector response to a single physical object is reconstructed as two different final-state objects, e.g. an electron reconstructed as both an electron and a jet, several steps are followed to remove such overlaps, as described in ref.
[13]. The overlap removal procedure is performed using candidate leptons.
The different lepton selections used in the analysis are summarised in table 3. Three types of signal lepton requirements are used for both the electrons and muons: 'tight', 'loose', and 'loose and minimally isolated'. The tight leptons, and the loose and minimally isolated leptons, are a subset of the loose leptons. The photon conversion veto is applied on top of the tight and loose electron selection requirements, only in the 2 sc channel. The Loose and FixedCutLoose isolation criteria applied for the loose and minimally isolated electrons and muons are described in refs.
The missing transverse momentum, with magnitude E miss T , is defined as the negative vector sum of the transverse momenta of all identified objects (muon candidates, electron -9 -

JHEP06(2021)146
Electrons Muons  Table 3. The requirements applied to define the categories of candidate leptons: loose (L), loose and minimally isolated (L * ) and tight (T) leptons. The overlap removal procedure is not applied for the candidate leptons. In the 2 sc channel, the photon conversion veto is required in addition to the loose and tight criteria.
candidates and jets) and an additional soft term [108, 109]. The soft term is constructed from all tracks that are matched to the primary vertex and are not associated with any other objects.

Event selection
Candidate events are selected for read-out using lepton triggers that require the one electron or muon to satisfy identification criteria similar to those used in the offline reconstruction, isolation criteria, and a transverse momentum requirement of p T > 26 GeV [110,111]. With increasing p T the requirements on identification and isolation become less stringent. All events must contain at least one offline tight lepton with p T > 30 GeVthat triggered the event. The event selection proceeds in two steps: the preselection and the signal regions (SRs) selection. The preselection is defined in table 4; the three channels (2 sc , 3 and 4 ) are defined to be mutually exclusive. Events are selected only if the absolute value of the sum of charges of the leptons is two, one, and two or zero for the 2 sc , 3 and 4 channels, respectively. In the 2 sc channel, the second-highest-p T lepton is required to have p T > 20 GeVand both leptons are required to be tight. Similarly in the 3 channel, each lepton in the pair of leptons of the same charge is required to have p T > 20 GeVand be tight.
Further preselection requirements are based on E miss T , the jet multiplicity N jets and the number of b-jets, N b-jet . The lowest E miss T value is 30 GeV in the 3 and 4 channels, and 70 GeV in the 2 sc channel; this selection helps to reduce the non-prompt lepton, electron charge-flip and V V backgrounds. In the 2 sc (3 ) channel only events with at least three (two) jets are considered. The background from top-quark production is highly reduced by requiring zero b-jets in the event. In order to reduce the background from Drell-Yan processes and neutral mesons, the invariant mass of same-flavour opposite-charge lepton pairs is required to be greater than 15 GeV for the 3 and 4 channels. The invariant mass of the same-flavour opposite-charge lepton pair must differ from the nominal Zboson mass by 10 GeV. For the 2 sc channel, the Z-boson invariant mass veto is also -10 -

JHEP06(2021)146
Selection criteria 2 sc 3 4 At least one offline tight lepton with p T > 30 GeVthat triggered the event GeV   Table 4. The preselection criteria for the 2 sc , 3 and 4 analysis channels. The leptons are required to pass the loose (L), loose and minimally isolated (L * ) or tight (T) requirements. The leptons are ordered by decreasing p T ( 1 , 2 , . . .) in the 2 sc and 4 channels, while for the 3 channel 1 and 2 denote the two same-charge leptons and 0 denotes the lepton with a charge opposite to the total lepton charge. Q denotes the charge of each lepton. SFOC refers to same-flavour opposite-charge lepton pairs. The symbol "-" means no requirement is made. The equal sign (=) is used to emphasise that the selection criterion has to be exactly the given number.
applied to e ± e ± events, in order to reduce the contributions originating from electron charge misidentification.
In addition to E miss T the following variables are used to define SRs: • The invariant mass of all selected leptons in the event, m x , where x can be 2, 3 or 4 corresponding to the 2 sc , 3 or 4 channels.
• The invariant mass of all jets in the event, m jets . When there are more than four jets in the event, only the leading four jets are used. This variable is only used for the 2 sc channel.
• The distance in η-φ between two same-charge leptons, ∆R ± ± . It is used for the 2 sc and 3 channels. In the 4 channel, two such variables can be calculated per event, ∆R min ± ± and ∆R max ± ± , denoting the minimum and maximum values, respectively.
• The transverse momentum of the highest-p T jet, p leading jet T . This variable is used in the 3 channel.
• The transverse momentum of the highest-p T lepton, p 1 T . This variable is used in the 4 channel.
• The azimuthal distance between the dilepton system and E miss T , ∆φ ,E miss T . It is only used in the 2 sc channel.

JHEP06(2021)146
• The smallest distance in η-φ between any lepton and its closest jet, ∆R jet . This variable is used in the 3 channel.
• The variable S, used for the 2 sc channel to describe the event topology in the transverse plane, and defined using the spread of the φ angles of the leptons, E miss T , and jets as follows: where the R is the root mean square that quantifies the spread, These variables are found to discriminate between the signal and the background. They exploit both the boosted decay topology of the charged Higgs bosons and the high energy of the decay products. Signal regions are defined for each channel, as summarised in table 5. The selection criteria defining the signal regions result from a scan of the multidimensional parameter space of the discriminating variables mentioned above. The SRs were designed by optimising the sensitivity for the H ±± pair production mode, using the m H ±± = 200, 300, 400 and 500 GeV mass hypotheses. The same SRs are used to study the H ± associated production mode; even though this approach is not optimal, it is still preferred instead of significantly increasing the number of signal regions. Additional mass hypotheses are used in the signal region defined for the lower mass hypotheses, since the signal discrimination power does not vary significantly in this regime. In particular, for the H ±± pair production mode the m H ±± = 300 GeVsignal regions are also used for m H ±± = 350 GeV. For the H ±± and H ± associated production mode, the m H ±± = 200 GeV, 400 GeV and 500 GeV signal regions are also used for m H ±± = 220 GeV, 450 GeV and 550 GeV, respectively. The SRs defined for the 2 sc channel are further divided into ee, eµ and µµ final states. Events in the 3 SRs are separated into two categories according to whether or not a same-flavour opposite-charge lepton pair exists in the event. This procedure further improves the expected significance by exploring differences in background composition and lepton-flavour composition between signal and background. The number of events observed in data is shown together with the expected signal and estimated background yields in section 9.

Background estimation
The background sources can be divided into two main categories. One category is populated by the SM events which contain only reconstructed charged (prompt) leptons originating from leptonic decays of W and Z bosons. The second category is formed by non-prompt leptons and charge-flip electrons. The non-prompt-lepton category refers to the leptons that originate from decays of b-and c-hadrons, or single pions that mimic electron signatures. The electrons from hadron decays into photons which convert into pairs of electrons in the beam pipe or detector material also enter this category. Lepton candidates reconstructed from these different sources share the properties of being less isolated, having  larger impact parameters relative to the primary vertex and being less likely to satisfy the lepton identification criteria. The backgrounds from SM processes with prompt leptons are estimated with MC simulations, except for background from W Z production, for which the normalisation is corrected using data in a dedicated control region. Background from non-prompt leptons and electron charge-flip are estimated using data-based methods. Background from V γ production can contribute in the SRs if the photon converts to an electron-positron pair, and is estimated using MC simulations. Background from W W production is estimated from simulation if the two W bosons have the same electric charge and from data if the two W bosons have opposite electric charge, since the latter only contributes through electron charge-flip and non-prompt leptons.

Background from W Z production
The W Z process is a dominant source of background in the 2 sc and 3 SRs. To correct a mismodelling seen in the jet multiplicity distribution [112], a normalisation factor is computed and then applied to the W Z background events containing two or more jets. The normalisation factor is measured in a dedicated W Z control region. It is selected by requiring exactly three tight leptons with p T > 20 GeV, at least two jets and no b-jet in the event. Finally, there must be at least one pair of same-flavour opposite-charge leptons with a mass compatible with the Z-boson mass (|m oc − m Z | < 10 GeV). The latter criterion ensures that the W Z control region and the 3 SRs are orthogonal.
The normalisation factor is derived in a fit of a first-order polynomial, as a function of jet multiplicity, to ratios of the data event yield (subtracting all non-W Z contributions) to the W Z event yield predicted in MC simulation. The value of the polynomial at N jets = 0 is used to scale the predicted W Z yields and found is to be 0.83 ± 0.07, where the uncertainty includes the statistical and systematic components. The different sources of systematic uncertainty are discussed in section 8. The jet multiplicity distribution in the W Z control region is shown in figure 2a, and for illustration the distribution in the 3 preselection region is also shown in figure 2b. The normalisation factor is applied to the W Z MC contribution. The sum of estimated backgrounds agrees with data within the assigned systematic uncertainties.

Electron charge-flip background
In the 2 sc channel, a background contribution is expected from events with oppositecharge lepton pairs where the charge of one of the leptons is misidentified. The charge-flip background is only significant for electrons and is mainly due to interactions of the electron with the ID material.
The misidentification rate is measured using a large data sample of dielectron events originating mainly from Z → e + e − decays selected by requiring two tight electrons with an invariant mass between 80 and 100 GeV. The sample contains mostly opposite-charge dielectron pairs, with a small fraction of same-charge dielectron pairs. The fraction of samecharge dielectron events is used to extract the charge misidentification rate as a function of the electron p T and η, using the method described in ref.
[105]. This rate is found to range between 0.01% and 4%, where higher values are obtained at large rapidities due to the larger amount of material traversed by the electrons. The statistical uncertainty of this estimate varies between 2% and 26% and is taken as a systematic uncertainty in the charge misidentification rate. The background in both the opposite-charge and same-charge samples is estimated in a cubic polynomial function fit of the high (100-120 GeV) and low (60-80 GeV) m sidebands. The uncertainty in the background is estimated by varying the sidebands and propagates to a 3% uncertainty in the estimated charge misidentification rate. The final systematic uncertainty combines all the sources mentioned.
The charge-flip background in a given region is estimated by selecting data events with opposite-charge dilepton pairs, but otherwise identical selection, and weighting them by the probability that the charge of the electrons is misidentified. Another source of systematic uncertainty is estimated by comparing, in simulated V +jets, tt, and W W events, the number of same-charge events estimated from opposite-charge events with the prediction; it accounts for differences between the charge misidentification rates in different processes and is found to be approximately 10%. For this test, the misidentification rate was measured using Z → e + e − MC simulations, using the same method as for the measurement performed with data.
To estimate the electron charge-flip background in the regions defined with loose electrons the methodology described above is applied. For this estimation, dedicated electrons misidentification rates measured with a Z → e + e − sample selected by requiring one tight and one loose electrons are used.

Non-prompt-lepton background
The composition of the non-prompt-lepton background in the SRs varies considerably among the analysis channels. Therefore, the methods to estimate these contributions are different for the 2 sc , 3 and 4 channels. In the 2 sc and 3 channels the non-prompt-lepton background is estimated using a fake-factor method [23], while the simulation is corrected with scale factors measured in data for the 4 channel.
Non-prompt-lepton background estimate for the 2 sc channel. The estimate of the background from non-prompt leptons assumes that these contributions can be extrapolated from a control region enriched in non-prompt leptons with a so-called fake factor.

JHEP06(2021)146
The control region is selected using the kinematic requirements of the preselection or the signal regions but alternative lepton selection criteria. The latter means that at least one of the selected leptons is required to satisfy the loose but not the tight lepton requirements.
The fake factors are calculated separately for electrons and muons in control regions with kinematic selections designed to enhance their content in non-prompt leptons. This is achieved by applying the preselection requirements of the 2 sc channel, except that E miss T must be lower than 70 GeV. Only events with electron and muon same-charge pairs are then used. The fake factor is defined as ratio of the number of events in the control region with all the selected leptons required to pass the tight signal requirements, to the number of events in the same region but with one of the selected leptons satisfying alternative lepton requirements. The measurement is performed as a function of the lepton p T . The fake factor dependency on the lepton η was also checked and found to be negligible. The SM processes with prompt leptons and the charge-flip contributions are subtracted in the control region. For the electron and muon fake-factor measurements, the lepton with the second highest p T is assumed to be the non-prompt one.
The measured fake factors are 0.03 ± 0.01 for electrons and muons up to p T = 40 GeV, and increase to 0.16 ± 0.05 and 0.09 ± 0.02, respectively, for electrons and muons with p T > 60 GeV. The uncertainties are statistical only. A systematic uncertainty of 20% (10%) in the electron (muon) fake factor is estimated by studying the variation of the fake factor with E miss T . For this measurement two E miss T bins are considered, < 70 GeV and > 70 GeV. This uncertainty accounts for the different compositions of the non-prompt leptons in the control region and the SRs. Another source of systematic uncertainty accounts for how often the non-prompt lepton is actually the one with the highest lepton p T and not the one with the second highest p T , as assumed. It is estimated from generatorlevel studies performed with MC simulations. This source is dominant in the region with the lepton p T > 60 GeV, where the uncertainty reaches 45% (80%) for electrons (muons). Uncertainties in the SM processes with prompt leptons (approximately 20%) and electron charge-flip (15%) background subtraction are also included. The overall uncertainty amounts to approximately 30% (20%) and 55% (85%) for the fake factors for electrons (muons) with 20 < p T < 60 GeV and p T > 60 GeV, respectively.
Non-prompt-lepton background estimate for the 3 channel. The same method as the one employed for the 2 sc channel is used in the 3 channel. Here the opposite-charge lepton 0 passes the loose selection, and it is assumed to be prompt, an assumption that is found to be valid in MC simulation. The control region used to calculate the fake factors uses the 3 preselection requirements, except that exactly one jet is required. Only events with electron and muon same-charge pairs are then used. The muon fake factor is found to be 0.03±0.01 and the electron fake factor is found to be 0.02±0.01, where the uncertainties are statistical only. A systematic uncertainty of 15% is estimated by measuring the lepton fake factor in a control region enriched in events from tt production. This uncertainty accounts for the different compositions of the non-prompt leptons in the control region and the SRs. Uncertainties in the subtraction of the SM processes with prompt leptons are found to be approximately 55% (45%) for the electron (muon) fake factor. Another 20%

JHEP06(2021)146
Sample Z+jets-enriched region tt-enriched region GeV - Table 6. The selection criteria that define the control regions enriched in non-prompt leptons used to determine the MC scale factors for the 4 channel. The leptons are required to pass the loose and minimally isolated (L * ) requirements. The transverse mass, m T , is calculated as the invariant mass of the vector sum of the transverse momentum of the non-prompt-lepton candidate and the missing transverse momentum. The symbol "-" means no requirement is made. The equal sign (=) is used to emphasise that the selection criterion has to be exactly the given number. uncertainty comes from generator-level studies performed with MC simulations to test the assumption that 0 is a prompt lepton. The fake factors' dependency on the lepton p T is also studied, and the deviations from the nominal fake factors are within the statistical uncertainty. When all sources of systematic uncertainty are combined, the total systematic uncertainty of the electron (muon) fake factor is found to be 60% (50%). -17 -

JHEP06(2021)146
The scale factors are measured to be λ e HF = 0.98 ± 0.18, λ e LF = 1.34 ± 0.17 and λ µ HF = 0.94 ± 0.04, where the uncertainties are statistical only. Systematic uncertainties of 5%, 4% and 2% for the λ e HF , λ e LF and λ µ HF scale factors, respectively, are estimated by varying the jet multiplicity and the lepton p T threshold in the nominal Z+jets-and ttenriched control regions. These uncertainties account for the different compositions of the non-prompt leptons in the control regions and the SRs. Uncertainties in the prompt-lepton subtraction are dominated by theoretical uncertainties. The total systematic uncertainties combine all the sources mentioned and are approximately 6%, 15% and 4% for the λ e HF , λ e LF and λ µ HF scale factors, respectively.

Validation
The data-based methods employed to estimate the backgrounds are validated by comparing the event yields in data with the combined predictions for these backgrounds, added to MC predictions for SM processes with prompt signal leptons. Distributions of selected variables in the 2 sc , 3 and 4 channels are shown in figures 3, 4 and 5, respectively, after the preselection requirements from table 4 are applied. Good agreement is observed in both normalisation and shape, demonstrating that the background contributions are well estimated. The expected distributions of both signal models are shown for m H ±± = 300 GeV to illustrate the discrimination power of the selected variables.
The signal contamination at the preselection level was studied for m H ±± = 300 GeV and m H ±± = 200 GeV mass hypotheses, corresponding to the H ±± pair production mode and the H ±± and H ± associated production mode, respectively. These mass hypotheses were selected because they are close to the sensitive mass range of the current analysis. A maximum value of 8.5% (2%) H ±± H ∓∓ (H ±± H ∓ ) signal contamination was found in all the individual 2 sc , 3 and 4 channels at preselection level.

Systematic uncertainties
Uncertainties in the signal and background yields arise from experimental uncertainties and from the theoretical accuracy of the prediction of the SM background yields. The experimental uncertainties arise from the luminosity determination, modelling of pile-up interactions, the reconstruction and identification of electrons, muons and jets and from the uncertainties associated with the data-based methods that are used to estimate the non-prompt lepton and charge-flip electron backgrounds.
The uncertainty in the combined 2015-2018 integrated luminosity is 1.7% [113], obtained using the LUCID-2 detector [114] for the primary luminosity measurements. Uncertainties in the modelling of pile-up interactions are estimated by reweighting the distribution of the average number of pp interactions per bunch crossing in the simulation, such that the average number of interactions changes by ±9%. The impact of these uncertainties on the background event yields estimated from MC simulation is found to be lower than 2%.
The uncertainties related to event reconstruction include the lepton [106, 115] and jet [100]   The overall impact on the background from SM processes with prompt leptons and the signal yields in the SRs from these systematic uncertainties is found to be lower than 10%. The dominant contribution comes from the jet energy scale component.
The uncertainties in the efficiencies of the electron [116] and muon [106] reconstruction, identification and trigger are also included. Their impact on the estimated yields in the SRs for the signal and the background from SM processes with prompt leptons is lower than 5%. The uncertainties in the b-jet identification are found to be negligible.     rections are evaluated by varying the renormalisation and factorisation scales independently by factors of two and one-half, and removing combinations where the variations differ by a factor of four. The uncertainties due to the PDF and the α s value used in the PDF determination are evaluated using the PDF4LHC prescription [117]. In the SRs defined for m H ±± = 300 GeV, the theory uncertainty in the background yields from ZZ and W W processes varies between 15% and 40%. For the background from V V V it is approximately 10%, for the background from V γ it is approximately 35% and for the tZ and ttX processes it is approximately 14%. The total uncertainty in the estimated W Z yields is 9% and includes the statistical uncertainty as well as two sources of theoretical uncertainty. Systematic uncertainties due to higher-order QCD corrections are evaluated using the same prescription as for the other diboson processes, and are found to be 3%. The second source is evaluated by comparing the results obtained with the linear fit model (section 7.1) in the W Z control region, and in a region defined with the W Z control region requirements except that three jets must be in the event. This uncertainty is found to be 8.4%. To validate the assigned uncertainty, several checks are performed. The linear fit function is changed to a quadratic one to check the quality of the fit model. The fit parameterisation is studied in the W Z-enriched region and in a region defined with the W Z control region requirements except that three jets must be in the event. The differences between the obtained results are found to be covered by the total uncertainty. The choice of parton shower model is studied in the W Z control region by comparing event samples simulated with Sherpa v2.2.1 and MadGraph5_aMC@NLO, and the differences are found to be covered by the total uncertainty. For the W Z background in regions with lower jet multiplicities, the uncertainty is estimated with the same methodology as for the other diboson processes described above.
An uncertainty of 50% is assigned to the other backgrounds (ttt, tttt, ttW W and V H), This large value is assigned to cover uncertainties from missing higher-order corrections and the PDF sets. Since these processes produce a larger number of jets at the first order of the perturbative expansion, they are less sensitive to parton shower modelling uncertainties.
The relative uncertainty in the background yields obtained in a fit of background to the observed data is shown in figure 6 for all SRs. The statistical uncertainties originate from the limited number of preselected and opposite-charge data events used in the fake-

JHEP06(2021)146
factor method and the charge-flip electron background estimate, respectively, as well as the effect of the limited number of simulated events for SM processes with prompt leptons. The total uncertainty is computed using all sources of uncertainty, and they are treated as uncorrelated. The uncertainties range from 10% to 30% and are dominated by the statistical uncertainties in the non-prompt-lepton estimate and the theory uncertainties. An exception is the 2 sc SR defined for m H ±± = 300 GeV, where the uncertainties from most sources are of similar size. The uncertainties associated with the charge-flip background are small in all 2 sc SRs. In the 4 channel, the statistical uncertainty of the non-prompt leptons comes from the limited number of events in the Z+jets, tt and tW MC simulation samples.
The theoretical uncertainties in the predicted signal yields arise from the parton shower model, missing higher-order corrections, and parton distribution functions. The systematic uncertainty due to the parton shower model is evaluated by comparing event samples simulated by Pythia 8 using the A14 tune with samples from Herwig using the H7-uemmht underlying-event tune [118], and is found to be less than 5%. The uncertainties due to the PDFs are found to be less than 5%. The uncertainty due to missing higher-order corrections is less than 10% [52]. When those uncertainties are combined in quadrature, an overall uncertainty in the signal yields of approximately 10% is obtained for the signal normalisation.

Results
The observed data event yields and the corresponding estimates for the backgrounds in the SRs are shown in figure 7. More details of the event yields in the signal region defined for m H ±± = 300 GeV are given in table 7. No significant excess over the expected yields is observed in any of the SRs. Table 7 includes the acceptance for the signal from pair production of H ±± bosons, A PP , and from associated production of H ±± and H ± bosons, A AP . It is defined as the number of selected reconstructed events divided by the number of events at the event generation stage and represents the signal reduction due to phase-space acceptance, branching ratio and detector efficiency. The results in all the SRs are shown in the appendix A.
The E miss T distribution is shown in figure 8 for the SRs of the m H ±± = 300 GeV signal mass hypothesis, where the selection requirement on E miss T has been removed. The other criteria summarised in table 5 used to define the SR are applied. Good agreement between data and the expected event yields is observed for low values of E miss T , demonstrating that the background contributions are also well estimated close to the SRs, where the requirements on the kinematic variables are tighter than at preselection.
The statistical interpretation is based on a likelihood ratio test [119] using the CL s method [120]. The signal strength, a free parameter in the fit, modifies the cross section of the signal hypothesis under investigation. A separate likelihood function is constructed for every signal hypothesis as the product of the Poisson probability distributions of the six channels of the corresponding signal region. Gaussian distributions are used to constrain the nuisance parameters associated with the systematic uncertainties. The widths of the    Table 7. The expected background and the observed data event yields in the signal region defined for the m H ±± = 300 GeV mass hypothesis. The signal yield is for the corresponding mass point and is normalised to the luminosity of 139 fb −1 . The displayed numbers include all sources of statistical and systematic uncertainties. The overall signal acceptances A PP and A AP and the observed upper limit on extra contributions to each signal region at 95% confidence level, n 95 , are also presented. Selections with µµ, three or four leptons are not affected by the electron charge-flip background, so these contributions are denoted by "-". Gaussian distributions correspond to the magnitudes of these uncertainties. Statistical uncertainties in the backgrounds are estimated using Poisson distributions.
The expected and observed upper limits on the H ±± pair production and the H ±± and H ± associated production cross sections times branching fraction at 95% confidence level (CL) for the five mass hypotheses are shown in figures 9a and 9b, respectively. They are obtained from the combination of 2 sc , 3 and 4 SRs. Assuming a linear dependence of the cross-section limit between neighbouring mass points, the observed 95% CL lower limit on the mass of the H ±± boson is 350 GeV for the pair production mode and 230 GeV for the associated production mode. To confirm the validity of the linear extrapolations close to the exclusion threshold, upper limits on the H ±± pair production cross section times branching fraction are also computed for the complementary mass hypothesis of 350 GeV, as shown in figure 9a, and are found to match the extrapolated values well. A similar test was also done for the H ±± and H ± associated production mode, and the upper limits on the cross section times branching fraction are computed for the 220 GeV, 450 GeV and 550 GeV complementary mass hypotheses. The results are shown in figure 9b, and again good matching is found.
A tighter limit on the H ±± boson mass is obtained for the H ±± pair production mode mainly due to the cross sections times branching fraction for H ±± pair production being higher than those for associated production of H ±± and H ± . Other important reasons are the different branching ratios of H ±± H ∓ (≈ 16%) and H ±± H ∓∓ (≈ 26%) decays into 2 sc , 3 and 4 leptons, and the signal acceptance, which is higher for the pair production mode than for the associated production mode (see tables 7 and 8). The definition of common signal regions also plays a role. As discussed in section 6, the SRs are not optimal for associated production of H ±± and H ± . For similar reasons, the upper limits on the and H ± associated production cross section times branching fraction at 95% CL obtained from the combination of 2 sc , 3 and 4 channels. The region above the observed limit is excluded by the measurement. The bands represent the expected exclusion curves within one and two standard deviations. The theoretical prediction [29] including the NLO QCD corrections [52] is also shown.
charged Higgs boson production cross sections times branching fraction are weaker for associated production of H ±± and H ± than for H ±± pair production.
For associated production of H ±± and H ± , the limit on the H ±± boson mass implies a constraint on the H ∓ boson mass which is at most 5 GeV different from the H ±± mass. Upper limits on the H ±± pair (H ±± and H ± associated) production cross section times branching fraction range from 15 fb to 5 fb (40 fb to 10 fb) when the assumed H ±± mass varies from 200 GeV to 600 GeV. Observed upper limits at 95% CL on the number of BSM events for each signal region are derived using the CL s prescription and the results are shown in table 7 for m H ±± = 300 GeV and for all SRs in the appendix A.

Conclusion
A search for pair production of H ±± bosons, and for associated production of H ±± and H ± bosons, in the context of a type-II seesaw model, is presented. For the H ±± pair production mode the H ±± boson mass is at least 100 GeV lower than the H ± boson mass, while for the associated production of H ±± and H ± bosons the mass difference between the doubly and singly charged Higgs bosons is at most 5 GeV. Only the bosonic decays H ±± → W ± W ± and H ± → W ± Z are studied. Dedicated signal regions in the 2 sc , 3 and 4 channels are defined as a function of different mass hypotheses within the model to look for evidence of H ±± and H ± bosons in the √ s = 13 TeV proton-proton collision data sample collected between 2015 and 2018 with the ATLAS detector at the LHC. The data sample corresponds to an integrated luminosity of 139 fb −1 . The increase in the total integrated luminosity allowed significant improvements in the lepton fake factors measurement and the estimation of the associated uncertainties. The data are found to be in good agreement with the estimated background for all channels investigated. Combining those channels, the model considered is excluded at 95% confidence level for H ±± boson masses below 350 GeV for the pair production of H ±± bosons, and below 230 GeV for the associated -26 -

JHEP06(2021)146
production of H ±± and H ± bosons. Upper limits on the pair (associated) production cross section times branching fraction range from 15 fb to 5 fb (40 fb to 10 fb) when the assumed H ±± mass varies from 200 GeV to 600 GeV. The results obtained for the H ±± pair production scenario raise the exclusion limits beyond those from a previous, similar search by ATLAS using a smaller 13 TeV data set by approximately 130 GeV. In addition, the associated production of H ±± and H ± bosons is also explored, extending the previous search.

A Supplementary results
The