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

The data collected at centre-of-mass energies of 188.6 GeV by ALEPH at LEP, corresponding to an integrated luminosity of 176.2 pb-1, are analysed in a search for pair-produced charged Higgs bosons H+/-. Three analyses are employed to select the taunutaunu, taunucs and cscs final states. No evidence for a signal is found. Upper limits are set on the production cross section as a function of the branching fraction BR(H+ to tau nu) and of the mass M(H+), assuming that the sum of the branching ratios is equal to one. In the framework of a two-Higgs-doublet model, charged Higgs bosons with masses below 65.4 GeV/c2 are excluded at 95% confidence level independently of the decay mode.


Introduction
In the Standard Model, particle masses are generated via the Higgs mechanism implemented using one doublet of complex scalar fields.In this process one physical state remains in the spectrum, known as the Standard Model Higgs boson.The most important phenomenological consequence of an extended Higgs sector is the appearance of additional spin-0 states, both neutral and charged.For example, with the addition of one more doublet, five physical states remain after the spontaneous breaking of the SU(2) L × U(1) Y symmetry: three neutral and a pair of charged Higgs bosons.Among the possible extensions of the Higgs sector, those obtained by adding more doublets are preferred because they naturally lead at tree level to M W ≃ M Z cos θ W , a relation very well verified by experiment.
The ALEPH data collected at centre-of-mass energies up to 184 GeV have been used in Refs.[1,2] to search for pair production of charged Higgs bosons predicted in models with two Higgs doublets.The negative result of the search was translated into a lower limit on the H ± mass m H ± of 59 GeV/c 2 at 95% confidence level.In this paper an update of the search based on the data collected at √ s = 188.6GeV (hereafter referred to as the 189 GeV data) is presented.The theoretical framework and underlying assumptions are the same as detailed in Refs.[1,2].The H + is assumed to decay predominantly into cs or τ + ν τ final states (and respective charge conjugates for the H − ).Other decay modes are not considered and B(H + → τ + ν τ ) + B(H + → cs) = 1 is assumed, but the analysis is equally sensitive to other hadronic decay modes.As a consequence, H + H − pair production leads to three final states (cssc, csτ − ντ /csτ + ν τ and τ + ν τ τ − ντ ) for which separate searches are performed.
This letter is organized as follows.After a brief description of the ALEPH detector in Section 2, the event selections are described in Section 3. The results and the conclusions are given in Sections 4 and 5.

The ALEPH detector
Only a brief description of the ALEPH subdetectors relevant for this analysis are given here.A more comprehensive description of the detector components is given in Ref. [3] and of the reconstruction algorithms in Ref. [4].
The trajectories of charged particles are measured with a silicon vertex detector, a cylindrical drift chamber, and a large time projection chamber (TPC).These are immersed in a 1.5 T axial field provided by a superconducting solenoidal coil.This system yields a resolution of δp T /p T = 6 ×10 −4 p T ⊕0.005 (p T in GeV/c).Hereafter, charged particle tracks reconstructed with at least four hits in the TPC, and originating from within a cylinder of length 20 cm and radius 2 cm coaxial with the beam and centred at the nominal collision point, are referred to as good tracks.
The electromagnetic calorimeter, placed between the tracking system and the coil, is a highly segmented sampling calorimeter which is used to identify electrons and photons 1 and to measure their energies.It has a total thickness of 22 radiation lengths at normal incidence and provides a relative energy resolution of 0.18/ √ E + 0.009 (E in GeV).The luminosity monitors extend the calorimetric coverage down to 34 mrad from the beam axis.
Muons are identified by their penetration in the hadron calorimeter, a 1.2 m thick iron yoke instrumented with 23 layers of streamer tubes, together with two surrounding layers of muon chambers.The hadron calorimeter also provides a measurement of the energy of charged and neutral hadrons with a relative resolution of 85%/ √ E (E in GeV).
The calorimetry and tracking information are combined in an energy flow algorithm, classifying a set of energy flow "particles" as photons, neutral hadrons and charged particles.From these objects, jets are reconstructed with an energy resolution of (0.60 √ E + 0.60) × (1 + cos 2 θ) where E in GeV and θ are the jet energy and polar angle, respectively.

Monte Carlo samples
Fully simulated Monte Carlo event samples reconstructed with the same program as the data have been used for background estimates, design of selections and cut optimization.Samples of all background sources corresponding to at least 20 times the collected luminosity were generated.The most important background sources are e + e − → τ + τ − , qq, fourfermion processes and two-photon collisions, simulated with the KORALZ [5], PYTHIA [6], KORALW [7] and PHOT02 [8] generators.
The signal Monte Carlo events were generated using the HZHA [9] generator.Samples of at least 2000 signal events were simulated for each of the various final states for charged Higgs boson masses between 50 and 75 GeV/c 2 .

The τ + ν τ τ − ντ final state
The final state produced by leptonic decays of both charged Higgs bosons consists of two acoplanar τ 's and missing energy carried away by the neutrinos.As this topology is identical to that expected from stau pair production with massless neutralinos, the "Large ∆M" selection described in Ref. [10] is used here to search for charged Higgs bosons in this channel.Efficiencies to select events from H + H − → τ + ν τ τ − ντ are of the The total expected background amounts to 15.5 events, mainly consisting of irreducible background from W + W − → τ + ν τ τ − ντ .In the data, 20 events are selected, in agreement with the expectation.The systematic uncertainty on the number of expected signal events is estimated to be 3.0%, dominated by the effect of limited Monte Carlo statistics (2.7%) and uncertainties on the cross-section for charged Higgs production (1.0%).The systematic error on the background is estimated to be 8%.This is dominated by the effect of limited Monte Carlo statistics (4%), uncertainties on the cross section for W pair production (2%), and uncertainities on the cross section for two-photon processes (7%).The systematic error on the luminosity is estimated to be 0.5%.

The csτ − ντ final state
The mixed final state csτ − ντ is characterised by two jets originating from the hadronic decay of one of the charged Higgs bosons and a thin jet plus missing energy due to the neutrinos from the subsequent decay of the charged Higgs and of the τ .
As a first step in the analysis the thrust of the event is required to be less than 0.96 and the total number of good charged tracks greater than 7. To reduce background from twophoton processes and the contribution of beam related background which is not simulated, the energy deposited in a 12 • cone around the beam axis is required to be less than 2.5% of the centre-of-mass energy.Background from e + e − → qq(γ) events is reduced by demanding that the polar angle θ miss of the missing momentum vector point away from the beam axis such that |cosθ miss | < 0.9.To reduce the background from W + W − → qq ′ ℓν ℓ where ℓ corresponds to an electron or muon, it is required that the events contain no identified lepton with a momentum greater than 10% of the centre-of-mass energy.
At this point the events are clustered into three jets using the JADE algorithm [11].The y cut value where the transition from two to three jets occurs is required to be greater than 0.001.The jet with the lowest charged track multiplicity is taken as the τ jet candidate.If at least two jets have the same multiplicity the τ jet candidate is taken to be the lowest momentum jet.The following quality cuts are then applied on the τ jet.It is required that the charged multiplicity of the τ jet be between one and three.The angle between the τ jet and the closest quark jet candidate is required to be between 30 • and 125 • .The energy of the τ jet boosted into the charged Higgs boson rest frame, defined as the frame recoiling against the hadronic system, is required to be less than 40 GeV.
To further suppress backgrounds the following four variables are used: The angle θ qq ′ between the two hadronic jets.
• The total transverse momentum p tot t divided by the total visible energy E vis .
• The production angle of the events, reconstructed from the sum of the momenta of the two quark jets.This production angle is then charge tagged with the tracks from the τ jet to form the variable θ ch prod .In the case of two tracks in the τ jet the charge of the highest momentum track is used.
• The χ 2 per number of degrees of freedom from a kinematic fit to the events using the constraints of energy and momentum conservation and the equality of the two masses produced in each side of the event.
The four variables are linearly combined in a discriminating variable D, displayed in Fig. 1a.The cut optimisation leads to D > 0.26.In the data collected at √ s = 189 GeV, 20 events are selected, in agreement with the background expectation of 22.6.The efficiencies of selection for a range of masses are given in Table 1.The mass of the hadronic system, and hence the mass of the charged Higgs candidates, is rescaled such that the energy of the two jets is equal to the beam energy in order to improve the resolution.The reconstructed masses of the candidate events are displayed in Fig. 1b.The cutoff near 80 GeV/c 2 is due to the influence of the θ qq ′ variable.
The systematic uncertainty on the number of expected signal events is estimated to be 4.1%.The main contributions are from Monte Carlo statistics (3.5%), calibration uncertainties (1.5%) and uncertainty on the cross section for charged Higgs production (1.0%).The systematic error on the background level is estimated to be 12%.The main contributions are from Monte Carlo statistics (3%), uncertainty on the cross section for the W + W − process (2%) and from the statistics of data/Monte Carlo comparisons (10%).The systematic error on the luminosity is estimated to be 0.5%.

The cssc final state
The hadronic decays of two charged Higgs bosons lead to a final state of four well separated jets, which can be combined into two dijets with equal masses.With respect to Ref. [2] the preselection and jet pairing method remain unchanged.
The following five variables are used: • The χ 2 per degree of freedom from a five-constraint kinematic fit.The constraints in the fit are from conservation of energy and momentum and the equality of the two charged Higgs masses in the event.
• The production polar angle θ prod , i.e., the angle between the charged Higgs boson momentum direction and the beam axis.
• The difference between the largest and the smallest jet energies, E max − E min .• The product of the minimum angle between any two jets, and the smallest jet energy, E min × θ qq ′ .
The variables are linearly combined into one discriminating variable D, shown in Fig. 2a.Events are accepted if D > 4.4.Including in D the charm tag used in Ref. [2] does not increase the discriminating power.
The analysis selects 263 events from the data for masses between 50 and 80 GeV/c 2 , corresponding to a background of 294.2 events expected from Standard Model processes.The fitted-mass distribution of the selected candidates can be seen in Fig. 2b.Efficiencies are of the order of 50% as shown in Table 1.
The systematic error on the number of events expected is estimated to be 3.0%.The main contributions are from Monte Carlo statistics (2.4%), statistics of data/Monte Carlo comparison (1.2%), and knowledge of signal cross sections (1.0%).The systematic error on the background level is estimated to be 3%.The main contribution is from knowledge of the W + W − cross section (2%).The contribution of the luminosity is estimated to be 0.5%.

Results
The numbers of candidates observed in the data collected at a centre-of-mass energy of 189 GeV are consistent with those expected from Standard Model processes for each of the three channels.Since, in addition, the mass distributions in the cssc and csτ − ντ channels do not show any significant accumulation outside the W + W − region (Figs. 1 and 2), the results of the three selections described in this letter are combined with those obtained using √ s = 172-183 GeV data to set a 95% confidence level upper limit on the cross section for pair production of charged Higgs bosons.
In setting the limits several new features with respect to Refs.[1,2] are to be noted.Full background subtraction is performed according to Ref. [13].The likelihood ratio test statistic is used.The confidence levels are calculated using the semi-analytical approach described in Ref. [14].Systematic errors are conservatively taken into account by reducing the efficiencies and subtracted backgrounds by one standard deviation.The reconstructed mass of the charged Higgs boson is used as discriminating variable for the cssc and csτ − ντ channels.
The upper limit on the H + H − production cross section at 188.6 GeV as a function of m H ± is shown in Fig. 3 for three values of B(H + → τ + ν τ ).The results from lower centreof-mass energies have been scaled to √ s = 188.6GeV according to the dependence of the cross section on the centre-of-mass energy.
In two-Higgs-doublet models the production cross section for H + H − depends, at lowest order, only on m H ± .The expected cross section at 188.6 GeV, corrected for initial state radiation, is shown in Fig. 3. Upper limits on the production cross section therefore translate into excluded domains for m H ± .
The result of the combination of the three analyses is displayed in Fig. 4. Charged Higgs bosons with masses less than 65.4 GeV/c 2 are excluded at 95% confidence level independently of B(H + → τ + ν τ ).The corresponding expected exclusion is 69.1 GeV/c 2 .For the values B(H + → τ + ν τ ) = 0, 0.5 and 1, 95% C.L. lower limits on m H ± are set at 69.9, 65.4 and 76.3 GeV/c 2 , respectively.

Conclusions
The search for pair-produced charged Higgs bosons in the final states τ + ν τ τ − ντ , csτ − ντ and cssc has been updated using 176.2 pb −1 of data collected at √ s = 188.6GeV.
No evidence of charged Higgs boson production has been found and upper limits have been set on the production cross section as a function of B(H + → τ + ν τ ) and of m H ± .Within the framework of two-Higgs-doublet models these results exclude at 95% confidence level charged Higgs bosons with masses below 65.4 GeV/c 2 independently of B(H + → τ + ν τ ) and assuming B(H + → τ + ν τ ) + B(H + → cs) = 1.Similar results have been reported by L3 [15].

Figure 1 :
Figure 1: (a) The distribution of the discriminating variable D described in the text for the semi leptonic channel at the level of preselection.(b) The reconstructed masses of the charged Higgs candidates after the cut on the discriminating variable.In both plots the points are the data, the open histograms are the Standard Model background and the hatched areas represent the charged Higgs signal expectation with m H ± = 70 GeV/c 2 .The signal is normalised arbitrarily in both plots.

Figure 2 :
Figure 2: (a) The distribution of the discriminating variable D described in the text for the cssc channel at the level of preselection.(b) The distribution of the Higgs candidate masses after the cut on the discriminating variable.In both plots the points are the data, the open histograms are the Standard Model background and the hatched areas represent the Higgs signal expectation with m H ± = 70 GeV/c 2 .The signal is normalised arbitrarily in both plots.

Figure 3 :
Figure 3: Upper limits at 95% confidence level on the H + H − production cross section at √ s = 188.6GeV for three values of B(H + → τ + ν τ ).The charged Higgs boson production cross section is shown as a solid curve.

Table 1 :
Selection efficiencies ǫ (in as a function of the charged Higgs boson mass m H ± ., as shown in Table 1 for a representative set of charged Higgs boson masses.