Search for charged Higgs bosons in e+e- collisions at energies up to sqrt(s) = 209GeV

A search for charged Higgs bosons produced in pairs is performed with data collected at centre-of-mass energies ranging from 189 to 209 GeV by ALEPH at LEP, corresponding to a total luminosity of 629 invpb. The three final states taunutaunu, taunucs and cscs are considered. No evidence for a signal is found and lower limits are set on the mass M_H+ as a function of the branching fraction B(H to taunu). In the framework of a two-Higgs-doublet model, and assuming B(H+ to taunu + B(H+ to cs) = 1 charged Higgs bosons with masses below 79.3 Gev/c2 are excluded at 95% confidence level independently of the branching ratios.


Introduction
The Standard Model of electroweak interactions requires only one doublet of complex scalar fields, resulting in a single neutral Higgs particle.The simplest extensions of the Standard Model assume two complex scalar-field doublets, with a total of eight degrees of freedom.As in the Standard Model, three of the degrees of freedom are associated with the longitudinal components of the W ± and Z bosons.The remaining five degrees of freedom appear as five physical scalar Higgs states: three neutral Higgs bosons and the charged Higgs bosons H ± .
In the two-Higgs-doublet case, the charged Higgs boson couplings are completely specified in terms of the electric charge and the weak mixing angle θ W .The production cross-section thus depends only on the mass m H ± .For masses accessible at LEP 2 energies, the charged Higgs boson decays with negligible lifetime and width into either cs/c b or τ + ν τ final states.Because the analyses are not sensitive to the quark flavour, and because the cs decay mode dominates over c b, cs stands for either cs or c b in the following.Therefore, B(H + → τ + ν τ )+B(H + → cs)=1 is assumed and H + H − pair production leads to three final states (τ + ν τ τ − ντ , csτ − ντ /csτ + ν τ and cssc) for which separate searches are performed.
The ALEPH data collected at energies up to 189 GeV have already been analysed and the search results published in Refs.[1,2,3].The negative result of the search, under the hypotheses specified above, was translated into a lower limit on the H ± mass of 65.5 GeV/c 2 at 95% confidence level (C.L.).Results from other experiments are given in Ref. [4].The present letter describes the search for pair-produced charged Higgs bosons using the data collected up to the end of data taking.An improved analysis has been designed for the fully leptonic channel.In the semileptonic search, the rejection of the W + W − background has been refined with a method based on a combination of the charge-tagged boson production angle and a τ polarization estimator.For the four-jet event selection, the linear discriminant analysis (LDA) has been re-optimized to account for the additional integrated luminosity collected at increased centre-of-mass energies.

The ALEPH detector and event samples
A complete and detailed description of the ALEPH detector and its performance, as well as of the standard reconstruction and analysis algorithms can be found in Refs.[5,6].Only those items relevant for the final states under study in this letter are summarized below.
The trajectories of the charged particles (called charged tracks in the following) are measured with the central tracking system, formed by a silicon vertex detector, an inner drift chamber and a large time projection chamber, all immersed in the 1.5 T axial magnetic field from a superconducting solenoidal coil.Electrons and photons are identified in the electromagnetic calorimeter, a highly segmented sampling calorimeter placed between the tracking device and the coil.Muons are identified in the hadron calorimeter, a 1.2 m thick iron yoke instrumented with 23 layers of streamer tubes, surrounded with two double layers of muon chambers.Together with the luminometers, the hermetic calorimetric coverage extends down to 34 mrad of the beam axis.The missing energy and momentum from, e.g., tau charged Higgs boson decays, are determined with an energy-flow algorithm which combines particle identification, tracking and calorimetry information into a set of energyflow particles, used in the present analyses.
The data analysed in this letter were collected at LEP between 1998 and 2000 at e + e − centre-of-mass energies ranging from 189 to 209 GeV, corresponding to a total integrated luminosity of 629 pb −1 .The details for each sample are given in Table 1.Fully simulated samples of events reconstructed with the same programs as the data were used for the background estimates, the design of the selections and the optimization of the selection cuts.The most important background sources are (i) difermion events (e + e − → τ + τ − and qq) simulated with the KORALZ [7] generator; and (ii) e + e − → W + W − and other four-fermion processes simulated with the KORALW [8] and PYTHIA [9] generators.Event samples of these background processes, corresponding to at least 20 times the collected luminosity, were generated.The W + W − cross sections predicted by RACOONWW [10] and YFSWW [11] were used as discussed in Ref. [12].Finally, the two-photon interactions (γγ → leptons) were simulated with the PHOT02 [13] generator.Samples of these events with at least six times the collected luminosity were generated.
The signal events generated with the HZHA [14] program were simulated for each of the final states and centre-of-mass energies (Table 1), and for charged Higgs boson masses between 45 and 100 GeV/c 2 .

Analyses
An event selection has been defined for each of the τ + ν τ τ − ντ , csτ − ντ /csτ + ν τ (hereafter referred to as csτ − ντ ) and cssc channels, and was optimized for B(H + → τ + ν τ ) = 100%, 50% and 0%, respectively.The selection criteria were chosen to achieve the highest 95% confidence level expected limit on the charged Higgs boson mass in the absence of signal.

The τ + ν τ τ − ντ final state
Events with two to six charged tracks (at least one and at most four of each sign) are considered.Leptonic events W + W − → ℓνℓ ′ ν (ℓ, ℓ ′ = e or µ) are rejected by requiring that the momentum of any identified electron or muon be less than 0.1 √ s.The events are then forced to form two jets with the JADE algorithm [15].An event is selected if both jet polar angles θ 1,2 satisfy | cos θ 1,2 | < 0.96, if their reconstructed masses are less than 3 GeV/c 2 and if each jet contains at least one charged track.To suppress the high cross section γγ → f f processes, the total visible mass is required to be in excess of 0.075 √ s, the momentum transverse to the beam is required to be greater than 10 GeV/c, and there must be no energy deposited in a cone of 12 • around the beam axis.The signal selection efficiency of the latter cut is corrected for the effect of the beam-related background, not included in the simulation, and is estimated from events triggered at random beam crossings.The relative loss of signal efficiency is about 7%.
Nearly coplanar tau pairs from e + e − → τ + τ − (γ) are rejected by requiring that the angle α between the two tau jets be less than 170 • and the angle between the projections of their momenta onto the plane transverse to the beam axis be less than 165 • .The missing energy is required to be greater than 80 GeV and the missing mass greater than 70 GeV/c 2 .In order to improve the W + W − background rejection, an LDA has been used to construct a discriminant variable D 0 from a combination of the following four quantities: • a charge-tagged angular variable calculated from the polar angles of the τ jets and their charges as • the angle α between the two tau jets; • the missing transverse momentum of the event P miss T ; • the value y 23 of the jet-clustering resolution parameter for which the transition from two to three jets occurs.
The optimal discriminant variable was found to be where α is in radians and P miss T in GeV/c.The distribution of D 0 is displayed in Fig. 1.This quantity is used as a discriminant variable in the derivation of the mass limit.
The signal event selection efficiencies, parametrized as a function of m H ± , are given in Table 2 for √ s = 206 GeV.The selection efficiencies are almost independent of the centreof-mass energy and increase only slightly with m H ± .For a signal with m H ± =85 GeV/c 2 and B(H + → τ + ν τ )=1, a total of 16.5 events is expected in the data taken at centre-ofmass energies between 189 GeV and 209 GeV.The numbers of events selected are given in Table 3, compared to the expectations from the Standard Model backgrounds, dominated by W + W − production.The systematic uncertainty on the number of expected signal events is estimated to be 3.1%, dominated by the effect of limited Monte Carlo statistics (2.4%) and the uncertainty on the cross section for charged Higgs boson production (2%).The systematic error on the background level is estimated to be 1.5%, dominated by the effects of limited Monte Carlo statistics (1.3%), by the uncertainty on the cross section for the W + W − process (0.5%) and the uncertainty on the cross section for two-photon production (5%).

The csτ − ντ final state
The mixed final state csτ − ντ is characterized by two jets originating from the hadronic decay of one of the charged Higgs bosons and a τ jet with missing energy due to the prompt neutrino as well as to the neutrino(s) from the subsequent τ decay.
The preselection is the same as that described in Ref. [3].In order to identify the τ jet an algorithm based on "minijets" is used as described in Ref. [16].If a minijet satisfies the τ -jet selection criteria, the rest of the event is clustered into two jets using the Durham [17] clustering algorithm.A kinematic fit is performed with the constraints of energy and momentum conservation and equality of the cs and τ + ν τ masses.If there is more than one τ candidate the combination with the lowest χ 2 is taken.In order to reject background from W + W − → (e/µ)νqq ′ , the measured energy of the τ jet boosted into the Higgs rest frame is required to be less than 0.175 √ s.The boost is performed using the information from the hadronic side of the event.
After this procedure the following four variables are chosen to further suppress the background: • the total missing transverse momentum of the event, P miss T ; • the isolation angle θ iso of the τ , defined as the half-angle of the cone around the τ jet direction containing 5% of the total energy of the rest of the event; • the χ 2 from the kinematic fit; • the decay angle θ ch τ , defined as the angle between the τ momentum in the Higgs boson centre-of-mass frame and the Higgs boson flight direction, charge-tagged with the charge of the τ , to exploit the asymmetry in the W system, absent for scalars.
The four variables are linearly combined into one variable, D 1 , defined as where P miss T is in GeV/c, and θ iso and θ ch τ are in radians.Events are selected by requiring that D 1 > −0.1.The background consists primarily of W + W − → ℓνqq ′ events.
Due to the scalar nature of the H + , the τ + from its decay is produced in a lefthanded helicity state, in contrast to the τ + 's from W + decays.Variables designed for the measurement of the τ polarization at LEP 1 [18] have been used to form an event-by-event helicity estimator, E τ .This variable, together with the charge-tagged production angle θ ch prod [3], is used to discriminate further between W + W − → τ νqq ′ and H + H − → csτ − ντ events.The two variables are combined into another variable, D 2 , defined as where θ ch prod is expressed in radians.The distribution of D 2 is shown in Fig. 2a.The cut optimization yields D 2 > −0.3 for m H ± =75 GeV/c 2 .The selection efficiencies are given in Table 2 as a function of the Higgs boson mass for √ s = 206 GeV.They are only weakly dependent on √ s.In the data collected between √ s = 189 and 209 GeV, the numbers of selected events are compared with the background expectations in Table 3.The fitted-mass distribution of the Higgs boson candidates is shown in Fig. 2b.For m H ± =77 GeV/c 2 , close to the sensitivity of this search, and for B(H + → τ + ν τ )=0.5, a total of 21.2 signal events is expected.
The systematic uncertainty on the number of expected signal events is estimated to be 3.0%.The main contributions are the finite size of the simulated event samples (2.2%), calorimeter calibration uncertainties (0.5%) and the uncertainty on the cross section for charged Higgs boson production (2%).The systematic error on the background level was estimated to be 3.9%.The main contributions are from limited statistics of the simulated event samples (2.5%), uncertainty on the cross section for the W + W − process (0.5%) and calibration uncertainties (3%).

The cssc final state
The hadronic decays of pair-produced charged Higgs bosons lead to a four-jet final state with equal mass dijet systems.The preselection remains unchanged with respect to Ref. [3].
A five-constraint kinematic fit is performed with energy-momentum conservation and equal dijet-mass constraints.In this fitting procedure, the errors on the jet energies and angles are parametrized as for the W mass measurement in the four-jet channel [19].The pairing is chosen as the dijet combination giving the minimum χ 2 .
To evaluate the mass difference between the two dijet invariant masses, momentum and energy conservation is imposed to rescale the energies of the four jets, fixing the jet velocities at their measured values.The mass difference ∆m between the two rescaled dijets is required to be smaller than 30 GeV/c 2 .
To improve the background rejection a linear discriminant D 3 is constructed, combining the following five variables: • the production polar angle θ prod , i.e. the angle between the Higgs boson momentum direction and the beam axis; • the difference ∆m between the two rescaled dijet masses; • the χ 2 of the 5C kinematic fit; • the product of the minimum jet energy E min and the minimum jet-jet angle θ qq ; • the logarithm of the QCD four-jet matrix element squared M QCD [20].
The optimized LDA coefficients were determined at √ s = 206 GeV with a cocktail of five charged Higgs boson masses ranging between 80 and 88 GeV/c 2 , leading to: For m H ± =75 GeV/c 2 and B(H + → τ + ν τ )=0, a total of 101.9 events is expected in the data.The efficiency does not depend on √ s.
After the complete selection, the comparison between data and simulation is displayed in Fig. 3b for the dijet invariant mass.The numbers of events observed in the data are compared in Table 3 to the expected background from Standard Model processes, dominated by W + W − production.An overall 2.4 standard deviation deficit with respect to expectation is observed.It is correlated with the deficit observed in the measurement of the W + W − hadronic cross section [12], which was ascribed to a statistical fluctuation.
The systematic error on the number of expected signal events is estimated to be 2.5%.The main contributions are from limited sample statistics (1.3%), uncertainty on the cross section for charged Higgs production (2%) and accuracy of the simulation (0.5%).The systematic error on the expected background, dominated by W + W − and qq production, is estimated to be 2.0%.The main contributions are from the simulated sample statistics (0.4% for W + W − and 1.6% for qq), the uncertainty on the cross section (0.5% for W + W − and 5% for qq), and the adequacy of the simulation (1.4% for W + W − and 2.1% for qq).

Results
No evidence for a signal is observed in the data.The results of the three selections have been combined to set a 95% C.L. lower limit on the mass of charged Higgs bosons.
Full background subtraction has been performed in setting the limit with the likelihood ratio test statistic [21].Systematic uncertainties are taken into account according to Ref. [22].To improve the sensitivity of the analysis, the charged Higgs boson mass has been used as a discriminating variable for the cssc and csτ − ντ channels.In the previous publications [1,2,3], only event counting was used in the τ + ν τ τ − ντ channel.In this analysis, the discriminant variable D 0 has been introduced in the limit setting procedure.
The result of the combination of the three analyses is shown in Fig. 4. Charged Higgs bosons with mass lower than 79.3 GeV/c 2 are excluded at the 95% C.L. independently of B(H + → τ + ν τ ).The corresponding expected exclusion is 77.1 GeV/c 2 .For the values B(H + → τ + ν τ ) = 0 and 1, 95% C.L. lower limits on m H ± are set at 80.4 GeV/c 2 (with 78.2 GeV/c 2 expected) and 87.8 GeV/c 2 (with 89.2 GeV/c 2 expected) respectively.The expected (dash-dotted) and observed (solid) exclusion curves are shown for the combination of the three analyses, using the full 189-209 GeV data set.
Upper limits can also be derived on the H + H − cross section at √ s = 200 GeV, as a function of the Higgs boson mass, for B(H + → τ + ν τ )=0, 50 and 100%.To combine the data at different centre-of-mass energies, the limit on the cross section was extrapolated to 200 GeV with the expected √ s dependence for the production of a charged scalar particle pair.The result is shown in Fig. 5 as a function of m H ± .

Conclusions
Pair-produced charged Higgs bosons have been searched for in the three final states τ + ν τ τ − ντ , csτ − ντ and cssc, with 629 pb −1 of data collected at centre-of-mass energies from 189 to 209 GeV.No evidence for Higgs boson production was found and lower limits were set on m H ± as a function of B(H + → τ + ν τ ), within the framework of two-Higgs-doublet models.Assuming B(H + → τ + ν τ )+B(H + → cs)=1, charged Higgs bosons with mass below 79.3 GeV/c 2 are excluded at 95% C.L., independent of B(H + → τ + ν τ ).

Figure 1 :
Figure 1: The distribution of the discriminant variable D 0 described in the text for the fullyleptonic channel.The points are the data, the open histogram is the Standard Model background and the hatched histogram represents the Higgs signal expectation, absolutely normalized, with m H ± = 85 GeV/c 2 .

D 1 2 )Figure 2 :
Figure 2: (a) The distribution of the discriminant variable D 2 described in the text for the semileptonic channel.(b) The distribution of the fitted mass of the Higgs boson candidates after the cut on D 2 .The points are the data, the open histogram is the Standard Model background and the hatched histogram represents the Higgs boson signal expectation with m H ± = 75 GeV/c 2 .The signal is arbitrarily normalized.

Figure 3 :
Figure 3: (a) The distribution of the discriminating variable D 3 .(b) The distribution of the reconstructed mass of the Higgs boson candidates after the cut on the discriminating variable.The points are the data, the open histograms are the Standard Model backgrounds and the hatched histogram represents the Higgs signal expectation for m H ± = 75 GeV/c 2 .The signal is arbitrarily normalized.

Table 1 :
Integrated luminosities, centre-of-mass energy ranges and mean centre-of-mass energy values for the data collected with the ALEPH detector from 1998 to 2000.

Table 2 :
The signal event selection efficiencies ǫ (in %), parametrized as a function of the charged Higgs boson mass m H ± , at √ s= 206 GeV.

Table 3 :
Numbers of candidate events and background expected from Standard Model processes, for each of the three years of data taking.