Searches for neutral Higgs bosons in e + e − collisions at centre-of-mass energies from 192 to 202 GeV

Searches for neutral Higgs bosons are performed with the 237 pb − 1 of data collected in 1999 by the ALEPH detector at LEP, for centre-of-mass energies between 191.6 and 201 : 6 GeV. These searches apply to Higgs bosons within the context of the Standard Model and its minimal supersymmetric extension (MSSM) as well as to invisibly decaying Higgs bosons. No evidence of a signal is seen. A lower limit on the mass of the Standard Model Higgs boson of 107.7 GeV =c 2 at 95% con(cid:12)dence level is set. In the MSSM, lower limits of 91.2 and 91.6 GeV =c 2 are derived for the masses of the neutral Higgs bosons h and A, respectively. For a Higgs boson decaying invisibly and produced with the Standard Model cross section, masses below 106.4 GeV =c 2 are excluded.


Introduction
Searches for neutral Higgs bosons of the Standard Model, its minimal supersymmetric extension (MSSM), and extensions allowing Higgs boson decays into invisible final states were performed using data collected by the ALEPH detector at LEP during 1999.The data sample was taken at four centre-of-mass energies, 191.6, 195.5, 199.5, and 201.6 GeV at which 28.9, 79.9, 86.3, and 41.9 pb −1 of data were collected, respectively.The total data sample corresponding to an integrated luminosity of 237 pb −1 was analysed to search for topologies arising from the e + e − → hZ Higgsstrahlung process supplemented by W and Z gauge boson fusion, and from the e + e − → hA associated pair-production process of the MSSM.The production cross section of the Higgsstrahlung process in the MSSM is reduced by a factor sin 2 (β − α), where tan β is the ratio of the vacuum expectation values of the two Higgs doublets and α is the mixing angle in the CP-even Higgs sector.The hA production cross section is proportional to cos 2 (β − α).For an invisibly decaying Higgs boson, the observable rate of the Higgsstrahlung process can be expressed as ξ 2 σ SM (e + e − → hZ), where ξ 2 is the product of the branching ratio to invisible decays and a model-dependent factor which reduces the cross section with respect to that in the Standard Model.
Searches for neutral Higgs bosons with the ALEPH detector have already been carried out up to a centre-of-mass energy of 188.6 GeV [1,2]; no evidence of a signal was found.A lower limit at 95% confidence level (CL) was set at 92.9 GeV/c 2 on the Standard Model Higgs boson mass.In the MSSM for the benchmark parameter set with maximal stop mixing [3], lower limits of 82.5 GeV/c 2 and 82.6 GeV/c 2 were derived on the masses of the h and A Higgs bosons, respectively.For an invisibly decaying Higgs boson, a lower limit of 95.4 GeV/c 2 was set, for a production cross section equal to that in the Standard Model.Similar searches have been performed by the other LEP experiments [4].
The higher centre-of-mass energies and integrated luminosity in the 1999 data substantially increase the experimental sensitivity for the detection of Higgs bosons with respect to previous results.Nevertheless, the background processes are the same as those described in Refs.[1,2].The theoretical framework, the event selections, the study of systematic uncertainties, and the result extraction are therefore very similar to those previously described.The differences mainly consist of reoptimization of event selections, the introduction of new discriminating variables, and improvements related to the simulation of signal processes.

ALEPH detector
The components of the ALEPH detector that are most relevant for the analyses presented here are summarized in this section.A more detailed description of the detector can be found in Ref. [5] and its performance in Ref. [6].
Three coaxial tracking devices are located inside a solenoidal superconducting coil which produces an axial magnetic field of 1.5 T. The vertex detector (VDET) [7] consists of two cylindrical layers of silicon wafers situated at average radii of 6.3 and 11.0 cm.Charged particles with a polar angle in the range |cos θ| < 0.88 (0.95) traverse at least two (one) VDET layers.The VDET is surrounded by an inner tracking wire chamber (ITC) and by a large time projection chamber (TPC), which measures up to 21 three-dimensional coordinates per 1 charged particle between radii of 31 and 180 cm.The tracking achieves a momentum resolution σ(p T )/p T of 6 × 10 −4 p T ⊕ 0.005, with p T in GeV/c.The resolution on the three-dimensional impact parameter of tracks can be parametrized as (34 + 70/p)(1 + 1.6 cos 4 θ) µm, with p in GeV/c.
The electromagnetic calorimeter (ECAL) is also situated inside the coil.It is segmented into projective towers of typically 0.9 • × 0.9 • , which allows electrons and photons to be identified within jets.Luminosity calorimeters of similar construction to ECAL are installed between the endcaps and the beam pipe and are treated as an extension of the ECAL.A silicon-tungsten sampling calorimeter extends the electromagnetic calorimeter coverage down to 34 mrad.
Outside the coil, a hadron calorimeter (HCAL) measures the hadronic energy, acts as a filter for the identification of muons, and serves as a return yoke for the magnetic field.The outermost detectors are two double layers of muon chambers.
The measurements of charged particle tracks and of energy deposits in the calorimeters, combined with the identification of photons, electrons, and muons, are used to produce a list of charged and neutral energy flow particles.Hadron jets, formed by clustering these particles, have an energy resolution of σ(E) = (0.60 , where E is the energy in GeV and θ the polar angle of the jet.The resolution on the jet angles is approximately 20 mrad in both θ and φ.
The tagging of b quark jets is accomplished by combining several discriminating variables together using neural networks.These variables include an impact-parameter-based probability [8], a displaced vertex χ 2 [9], and the transverse momentum of identified leptons with respect to the jet axis [10].Three jet-shape quantities supplement these variables.The resulting six variable neural network b-tag is described in Ref. [11]; a more specialized four variable b-tag used for the e + e − → hA→ b bb b final state is described in Ref. [1].

Signal and background simulation
For each of the four centre-of-mass energies, fully simulated samples of signal and background events were generated.The hzha program [12], used to generate the signal events, was also used to compute the signal cross sections, the Higgs boson decay branching fractions, and the radiative corrections to the Higgs boson masses in the MSSM.The most recent version contains, in particular, the latest refinements on these radiative corrections [3,13,14] and a more complete simulation of the he + e − and hν e νe final states which now includes the interference between the boson fusion and the Higgsstrahlung processes [15].
The e + e − → qq(γ) background was generated with koralz [16], instead of pythia [17], to benefit from a more accurate treatment of multiple initial state photon radiation.The e + e − → Weν events, generated with pythia, were compared to events generated with grace4f [18], and the differences were taken into account in the systematic uncertainties.The other Standard Model four-fermion processes were simulated in a manner similar to that described in Refs.[1,2].The sizes of the simulated signal and background samples correspond to at least 50 times the collected luminosity used in the analyses.

Event selections
Analyses had already been developed [1,2] for most final states relevant for the searches for neutral Higgs bosons.For the visible e + e − → hZ search, these final states are the leptonic final state (h + − where denotes an electron or a muon), the missing energy final state (hν ν), the four-jet final state (hqq), and the tau final states (hτ + τ − and h→ τ + τ − , Z→ qq).For the e + e − → hA search, they are the four-b final state (b bb b) and the tau final state (b bτ + τ − ).For the invisible e + e − → hZ search, they are the acoplanar lepton pair final state (h + − ) and the acoplanar jet pair final state (hqq).As in Ref. [1], the visible Higgs boson search was conducted in two alternative analysis streams, a neural-network-based stream (NN) and a cut-based stream (cut), although both streams share the same leptonic and tau-final-state selections.This section highlights changes with respect to Refs.[1,2].Most of the hZ and hA selections were optimized for Higgs boson masses of 107 and 90 GeV/c 2 , respectively, i.e., close to the search sensitivity in both cases.
A reoptimization of the selection criteria and retraining of the neural networks in all channels were performed to account for the higher centre-of-mass energies and integrated luminosity.In some cases, several new discriminating variables were adopted and refinements to analysis techniques implemented, as described below.
• In the missing energy channel, the reconstructed Higgs boson mass was removed from the variables of the single neural network selection (referred to as "A") in order to increase the selection efficiency for low Higgs boson masses.The structure of the three-neural-network selection (referred to as "B") remains unchanged.
• In the tau final state selection, the overlap with leptonic-final-state events (he + e − or hµ + µ − ) was reduced by requiring that the measured invariant mass of the two tau jets be smaller than 75 GeV/c 2 when a decay particle of at least one of the two taus is identified as an electron or a muon.
• In the four-jet channels, the following improvements were made.
i) A four-constraint fit, identical to that used in Ref. [19], was adopted in all selections, instead of a simple energy rescaling.It allows the jet angles to be fitted as well as the jet energies, thus improving the dijet mass resolution, particularly near the kinematic threshold.
ii) In the four-b channel, a new variable was introduced to characterize the unbalanced topology of three jets recoiling against one.For each of the possible three-jet combinations, the energy-weighted angular dispersion is computed as where the indices i, j, and k refer to jets, θ l, u is the angle between the momentum direction of the object l contained in jet i, j, k and the vectorial sum u of the three jet momenta, and E l is the energy of object l.The variable ∆θ 3 is defined as the minimum value of this dispersion for all four possible three-jet systems.iii) In the neural network selection, the former 17 variables were supplemented with the decay angles relative to the direction of flight of the Z and Higgs boson candidates, α 12 and α 34 .The output of the neural network for the four-jet final state is shown in Fig. 2.

Systematic uncertainties
Systematic uncertainties caused by inaccuracies of the simulation were evaluated as in Refs.[1,2].Whenever possible, they were extracted from 3.5 pb −1 of data taken at the Z peak during the same year.A discrepancy between data and simulation, affecting the acceptance for signal and for most of the background processes, was observed in the impact-parameterbased b tagging quantities.To correct for this effect, a smearing of the track parameters was performed on the simulated events to bring them in better agreement with the data.Half of the correction was conservatively taken as a systematic uncertainty for all channels.The systematic uncertainties on signal efficiencies are typically 5%, and those on the background yields are summarized in Table 1.

Selection overlaps
Events selected by more than one analysis were treated with a well-defined prescription to remove overlaps.This procedure prevents any given event from affecting the confidence level calculations more than once.The potential overlaps of the searches for different topologies arising from the Higgsstrahlung process were suppressed by assigning a precedence order to the four selections, and by rejecting any candidate event already selected by a higher precedence search.The order of the precedence was chosen to be 1) the missing energy final state; 2) the tau final states; 3) the leptonic final state; and 4) the four-jet final state.In the data, no event was affected by this procedure.
The significant overlaps of the hZ and hA searches for tau final states as well as for four-jet final states in the cut stream were treated in a manner similar to that used in Ref. [1] for the missing energy channel.This method divides the analyses into three statistically independent branches: one branch containing the events selected only by the hZ search; a second branch containing the events selected only by the hA search; and a third branch containing the events selected simultaneously by both searches.The use of different reconstructed Higgs boson mass definitions in the hZ searches (with the Z mass constraint) and the hA searches (with an equal mass constraint) introduces an additional difficulty, solved in the following manner.The hZ (hA) mass definition is used for the hZ(hA)-only branch.For the third branch, the mass definition depends on the MSSM parameter set and is chosen such that it maximizes the combined expected 95% CL sensitivity.Typically, the hZ mass definition is chosen for small Table 1: The numbers of signal and background events expected, and the numbers of candidate events observed in the data for each channel.The signal expectation is computed for the hZ process with a Higgs boson mass of 107 GeV/c 2 , and that for the hA process with a common Higgs boson mass of 90 GeV/c 2 .The background yield is subdivided into ZZ (including Zee and Zν ν), WW (including Weν), and f f (including γγ → f f).The f f composition is entirely qq for the visible Higgs boson and hadronic invisible Higgs boson decay channels, and + − and γγ → + − for the leptonic invisible Higgs boson decay channel.Systematic uncertainties are given for each background.The three-branch subdivision in the NN stream four-jet final state was replaced by a simpler and similarly sensitive two-branch treatment: the first branch (4b) containing the events selected by the four-b final state search, and the second branch (2b) containing the events not selected by the 4b search.

Selection results
For each analysis, the numbers of signal and background events expected, and the observed number of events are summarized in Table 1.For the Standard Model Higgs boson searches, a total of 76 events are selected in the data with the NN stream, while 97.8 events are expected from Standard Model processes.Figure 3 shows the reconstructed Higgs boson mass distributions for the data and the background expectations for the hZ selections.The hA selections yield 5 events in the data with the NN stream, while 9.4 events are expected from Standard Model processes.In the search for invisibly decaying Higgs bosons, a total of 58 candidate events are selected with 60.8 events expected from Standard Model processes; the distributions of their reconstructed Higgs boson mass are shown in Fig. 4. Fewer events were selected in the data than were expected.This deficit is mostly due to the four-jet analysis of the NN stream where 46.4 events were expected, but only 30 events were selected in the data.In the four-jet cut analysis 27 events are selected in data with 36.3 expected.As can be seen in Fig. 3, the deficit of events is mostly near the mass region dominated by e + e − → ZZ events.
At the preselection level, the background expectation is dominated by four-jet e + e − → W + W − events which are kinematically similar to signal events, and no significant discrepancies were found between data and Monte Carlo simulation for all relevant kinematic variables.Figure 2 shows the distribution of the neural network output distribution at the preselection level.Overall, the agreement between the simulated expectation and the data is very good with 2572 events expected and 2523 seen in the data.The deficit is, however, apparent in Fig. 2 at high values of the neural network output, which is dominated by b tagged e + e − → ZZ events.
As mentioned in Section 3.3, a slight disagreement was found and corrected for in the b tagging distributions.The corrections were determined from Z peak data taken during the same year.These were verified with e + e − → qqγ radiative return events and semileptonic e + e − → W + W − → qq ν events collected at high centre-of-mass energies.The qqγ events provide a known b b content while the semileptonic W + W − events provide a control sample with negligible b-jet content.The impact of these b-tagging corrections is propagated to the simulation, but the effect is too small to account for the observed deficit of events.The conclusion of these studies is that the cause of the deficit is most likely due to a statistical

Exclusion limits
The 176.2 pb −1 of data collected in 1998 with a centre-of-mass energy of 188.6 GeV were combined with the 1999 data to determine the final results.As was done in Refs.[1,2], confidence levels on the signal hypothesis were drawn using the Signal Estimator method [20] rather than with the method used by the other LEP experiments [21], thus increasing the expected mass limit by ∼ 0.5 GeV/c 2 .Expected confidence levels were calculated as the median confidence level in the absence of signal.The tau final state analyses and the cut analyses for the four-jet, four-b, and missing energy topologies use only the reconstructed Higgs boson mass as a discriminant in the confidence level calculations.The other analyses use an additional discriminant variable: the neural network output in the case of the 2b branch of the NN stream four-jet analysis and a b-tag sensitive variable in the 4b branch as well as in the h + − final state analysis.
The expected and observed confidence levels are shown in Fig. 5 for the Standard Model Higgs boson searches.The low confidence levels for the background hypothesis at small Higgs boson masses are a reflection of the observed deficit of events discussed in Section 4.2.The lower limits at 95% confidence level on the mass of the Standard Model Higgs boson with systematic uncertainties taken into account [22] are summarized in Table 2. Since the NN stream gives an expected Higgs boson mass lower limit higher by about 0.5 GeV/c 2 , it is used to set the final results.
Figure 6 shows the hZ, the hA, and their combined Higgs boson mass exclusions as a function of sin 2 (β − α).These combined exclusions are interpreted within two benchmark parameter sets, one which maximizes the radiative corrections to the lighter CP-even Higgs boson mass as a function of tan β (referred to as the m max h scenario) [14] and one in which no mixing in the scalar top sector is assumed.The excluded domains in the [m h ,tan β] plane are shown in Fig. 7 as obtained using the NN stream.The results in the m max h scenario are summarized in Table 2.The tan β exclusion range is sensitive to the top quark mass assumption which was taken to be 175 GeV/c 2 .For a top quark mass of 180 GeV/c 2 , the tan β exclusion range reduces to [1.0,1.5].
For the invisibly decaying Higgs boson search, no discriminant variable is used in the confidence level calculations, as was done in the previous search at 188.6 GeV [2], because the sliding analysis technique effectively acts as if the reconstructed mass were used as a discriminant variable.The results are interpreted as an exclusion domain in the (m h ,ξ 2 ) plane, presented in Fig. 8.For ξ 2 = 1, invisibly decaying Higgs bosons with masses below 106.4 GeV/c 2 are excluded at 95% CL with an expected limit of 105.6 GeV/c 2 .

Conclusions
Searches for neutral Higgs bosons in e + e − collisions at centre-of-mass energies from 191.6 to 201.6 GeV have been carried out with the ALEPH detector using an integrated luminosity of 237 pb −1 .The major event topologies expected from the hZ process for visibly as well as invisibly decaying Higgs bosons, the weak boson fusion process, and the hA process have been analysed.The data were combined with the 176 pb −1 collected in 1998 at a centre-of-mass energy of 188.6 GeV.
A 95% CL lower limit on the Standard Model Higgs boson mass of 107.7 GeV/c 2 is obtained with an expected sensitivity of 107.8 GeV/c 2 .Lower limits of 91.2 and 91.6 GeV/c 2 are obtained at 95% CL for the masses of the h and A neutral Higgs bosons of the MSSM.These limits, determined in the m max h scenario, are valid for any tan β > 0.5.For a top quark mass of 175 GeV/c 2 , the tan β range between 0.8 and 1.9 is excluded with 95% confidence.
An invisibly decaying Higgs boson with a production cross section equal to that in the Standard Model is excluded at 95% CL for masses below 106.4 GeV/c 2 .

Figure 1 :
Figure 1: Distribution, at preselection level, of the minimum value of the energy-weighted angular dispersion ∆θ 3 in the four-b channel for data (dots with error bars), simulated background (solid histogram), and the signal for m h = m A = 90 GeV/c 2 (dashed histogram).The signal histogram has an arbitrary normalization.

Figure 2 :
Figure 2: Distributions of the neural network output used to select four-jet hZ candidates for data (dots with error bars), simulated background (solid histograms), and simulated Higgs signal with an arbitrary normalization for m h = 107 GeV/c 2 (dashed histogram), after preselection criteria have been applied.

Figure 3 :
Figure 3: Distributions of the reconstructed Higgs boson mass M h for the data (dots with error bars) selected in the hZ searches by (a) the NN stream, and (b) the cut stream.The histograms show the Standard Model background expectation.

Figure 4 :
Figure 4: Distribution of the reconstructed Higgs boson mass M h for the data (dots with error bars) and the expected background (histogram) in the invisibly decaying Higgs boson searches for the acoplanar lepton and the acoplanar jet final states.

Figure 5 :Figure 6 :
Figure 5: Observed (solid) and expected (dashed) CL curves for (a) the background hypothesis and (b) the signal hypothesis for the NN stream, as functions of the hypothesized Standard Model Higgs boson mass m h .For the signal hypothesis (b), the intersections of the horizontal line at 5% with the curves define the observed and expected 95% CL lower limits on the Standard Model Higgs boson mass.

Figure 7 :Figure 8 :
Figure 7: The expected (dashed) and observed (solid) excluded regions at 95% CL for the NN stream in the [m h , tan β] plane of the MSSM for (a) the m max h parameter set and (b) the no mixing scenario.The densely hatched regions are not allowed theoretically.

Table 2 :
Summary of the 95% CL lower limits on Higgs boson masses in the Standard Model and in the MSSM, and their expected values in absence of signal.The tan β excluded ranges in the m max h scenario are also given for a top quark mass of 175 GeV/c 2 .