Searches for Neutral Higgs Bosons in e+e- Collisions around sqrt(s) = 189 GeV

Searches for neutral Higgs bosons in the Standard Model and the MSSM have been performed using data collected by the DELPHI experiment at a centre-of-mass energy of 188.7 GeV, corresponding to an integrated luminosity of 158 pb^{-1}. These analyses are used, in combination with our results from lower energies, to set new 95% confidence level lower mass bounds on the Standard Model Higgs boson (94.6 GeV/c^2) and on the lightest neutral scalar (82.6 GeV/c^2) and neutral pseudoscalar (84.1 GeV/c^2) Higgs bosons in a representative scan of the MSSM parameters. The results are also interpreted in the framework of a general two-Higgs doublet model.


Introduction
In the framework of the Standard Model (SM) there is one physical Higgs boson, H, which is a neutral CP-even scalar. At LEP II the most likely production process is through the s-channel, e + e − → Z * →HZ. The W + W − and ZZ fusion t-channel production processes in some of the channels described below are not considered here, but their contribution is typically below 10% in the range of masses considered in this study.
The results of the search for the SM Higgs are also interpreted in terms of the lightest scalar Higgs boson, h, in the Minimal Super-symmetric Standard Model (MSSM). This model predicts also a CP-odd pseudo-scalar, A, produced mostly in the e + e − → hA process at LEP II. This associated production is also considered in this paper.
With the data taken previously at √ s= 183 GeV DELPHI excluded a SM Higgs boson with mass less than 85.7 GeV/c 2 [1], and set limits on h and A of the MSSM of 74.4 GeV/c 2 and 75.3 GeV/c 2 respectively. The present analyses therefore concentrate on masses between these and the kinematic limit. Note that the LEP Higgs working group [2] has found mass limits of 89.7 GeV/c 2 for H, 80.1 GeV/c 2 for h and 80.6 GeV/c 2 for A, under assumptions generally referred to as the benchmark scan, when combining the data of the four LEP experiments from data taken up to 183 GeV.
In the HZ channel, all known decays of the Z boson have been taken into account (hadrons, charged leptons and neutrinos) while the analyses have been optimised either for decays of the Higgs into bb, making use of the expected high branching fraction of this mode, or for Higgs boson decays into a pair of τ 's. A dedicated search for the invisible Higgs boson decay modes will be reported separately. The hA production has been searched for in the 4b and bb τ + τ − channels.
There are separate analyses for the different decay modes of the Higgs and Z bosons. Some common features are discussed in Sect. 3, the Hµ + µ − and He + e − channels in Sect. 4, channels involving jets and τ 's in Sect. 5 and Hνν in Sect. 6. Purely hadronic final states are discussed in Sect. 7. The results are presented in Sect. 8.

Data samples overview and the DELPHI detector
For most of the data collected in 1998, LEP was running at energies around 189 GeV. DELPHI recorded an integrated luminosity of (158±1) pb −1 at a mean energy of 188.7 GeV.
Large numbers of background and signal events have been produced by Monte Carlo simulation using the DELPHI detector simulation program [3]. The size of these samples is typically about 100 times the luminosity of the collected data. Background was generated with PYTHIA [4] and KORALZ [5] for (e + e − → ff γ), PYTHIA and EXCALIBUR [6] for the four-fermion background and TWOGAM [7] and BDK [8] for two-photon processes. BABAMC [9] was used to simulate Bhabha events in the main acceptance region.
Signal events were produced using the HZHA [10] generator. For the SM process the Higgs mass was varied in 5 GeV/c 2 steps from 70 GeV/c 2 to 100 GeV/c 2 , while for hA , the A mass was varied between 70 and 90 GeV/c 2 with tan β (the ratio of the vacuum expectation values of the two doublets) either 2 or 20. This fixes the h mass, almost equal to m A for tan β = 20 and significantly lower than m A if tan β = 2.
The HZ simulated samples were classified according to the Higgs and Z boson decay modes. For He + e − , Hµ + µ − and Hνν the natural SM mix of H decay modes was permitted. In the Hqq channel the τ τ decay mode was removed, and we generated separately the two channels involving τ leptons for which one of the bosons is forced to decay to τ 's and the other hadronically. Finally, for the hA simulations final states involving either four b quarks or two b quarks and two τ 's were simulated. Efficiencies were defined relative to these states. The size of these samples varied from 2000 to 20,000 events.
The detector was unchanged from the previous data taking period. Thus we refer to our previous publication [1] for a short description. More details can be found in references [3,11].
3 Common features for all channels 3

.1 Particle selection
In all analyses, charged particles are selected if their momentum is greater than 100 MeV/c and if they originate from the interaction region (within 10 cm along the beam direction and within 4 cm in the transverse plane). Neutral particles are defined either as energy clusters in the calorimeters not associated to charged particle tracks, or as reconstructed vertices of photon conversions, interactions of neutral hadrons or decays of neutral particles in the tracking volume. All neutral clusters of energy greater than 200 MeV (electromagnetic) or 500 MeV (hadronic) are used; clusters in the range 100-500 MeV are considered with specific quality criteria in some analyses. The π ± mass is used for all charged particles except identified leptons, while zero mass is used for electromagnetic clusters and the K 0 mass is assigned to neutral hadronic clusters.

b-quark identification
The method of separation of b quarks from other flavours is described in [12], where the various differences between B-hadrons and other particles are accumulated in a single variable, hereafter denoted x b for an event and x i b for jet i. A major input to the combined variable is the probability that all tracks in a group originate from the interaction point. x b combines this probability with information from secondary vertices (the mass computed from the particles assigned to the secondary vertex, the rapidity of those particles, and the fraction of the jet momentum carried by them) and also the transverse momentum (with respect to its jet axis) of the leptons, using the likelihood ratio technique. Increasing values of x b correspond to increasingly 'b-like' events (or jets).
The procedure is calibrated on events recorded in the same experimental conditions at the Z resonance. The performance of the combined b-tagging is described in Ref. [13], and the impact parameter tagging in Ref. [14]. The overall performance is illustrated in Figure 1.
A careful study of possible systematic effects, including data versus simulation agreement at the Z pole (checked inclusively, per flavour and for multi-jet events), leads to an overall relative b-tagging uncertainty below 5%, varying slightly with the exact tagging value used. At high energy, an inclusive comparison of data with simulation confirms this number. The gluon splitting rates into bb and cc have been rescaled in the simulation according to the DELPHI [15] measurement. In addition, a 50% uncertainty on these splitting rates is applied to the qq (γ) background estimate.

Constrained fits
In most channels a constrained fit [16] is performed to extract the Higgs mass, and often to reject background processes as well. If only total energy and momentum conservation are imposed then the fit is referred to as '4-C', while some fits require the Z mass as well, either as a fixed value, or taking into account the Breit-Wigner shape of the Z resonance. In both cases such fits are referred to as '5-C'. In order to allow the removal of most of the radiative return to the Zevents, an algorithm has been developed [17] in order to estimate the effective energy of the e + e − collision. This algorithm makes use of a '3-C' kinematic fit in order to test the presence of an initial state photon parallel to the beam direction and then lost in the beam pipe. This effective centre-of-mass energy is called √ s ′ throughout the paper.

Confidence level calculations
The procedure used to compute the confidence levels is the same as that used in our previous publications [1,18] but the discriminant information is now two-dimensional: the first variable used is the reconstructed mass, the second one is either the x b (for the electron and muon channels) or the likelihood (for all other channels). As far as the mass information is concerned, the reconstructed Higgs boson mass is used in the hZ channels and the sum of the reconstructed h and A masses in the hA channels (for the pairing with minimal mass difference in the four-jet channel). In order to make full use of the information contained in the second variable the selections are looser than in the past: the method used for deriving the confidence levels ensures that adding regions of lower signal and higher background can only enhance the performance expected from a tighter selection. Since the distributions, represented as two-dimensional histograms, are derived from simulation samples, the limited statistics in some bins are a potential problem: statistical fluctuations can artificially increase the expected sensitivity. This possible systematic shift of the confidence level has been estimated comparing the expected results using the full simulation sample with those derived from fractions of this sample. The bin sizes were carefully chosen to keep full sensitivity while avoiding any significant bias caused by this effect.
4 Higgs boson searches in events with jets and electrons or muons

Electron channel
The analysis follows what has already been published [1], with the following improvements in the selection cuts. To reinforce the Bhabha veto, the preselection described in [1] has been complemented by a rejection of electron candidate pairs having acoplanarity (defined as the supplement of the angle between the transverse momenta of the two electrons) below 3 degrees and energies above 40 GeV. To allow for the tau decays of the Higgs boson while keeping a good purity, the requirement on the minimum event charged multiplicity has been raised to 8 except if the recoiling system from the electron candidate pair is made of two jets, each with a charged multiplicity lower than or equal to 3 and with a mass below 2 GeV/c 2 . This defines the preselection.
The energy of the slower/faster electron is required to be above 15/20 GeV. Electron isolation angles with respect to the closest jet are required to be more than 20 • for the most isolated electron and more than 8 • for the other. Global kinematic fits [16] are performed, imposing total energy and momentum conservation and constraining the invariant mass of the e + e − system to m Z (5-C fits). If the fit probability is below 10 −8 , the fit procedure is redone for fixed values of the e + e − mass between 36.5 and 105 GeV/c 2 , in order to allow for the tails of the Z mass distribution. For each mass a combined probability is determined as the product of the χ 2 probability times the Breit-Wigner probability for the Z mass, and the mass giving the maximum combined probability is retained. Events with a combined probability below 10 −8 are rejected. As the search is restricted to high mass Higgs bosons produced in association with a Z particle, the sum of the masses of the electron pair and of the recoiling system as given by the kinematic fit is required to be above 150 GeV/c 2 and the difference in the range from -100 GeV/c 2 to 50 GeV/c 2 . The fitted recoil mass and the global b-tagging variable x b are used for the two-dimensional calculation of the confidence levels. Table 1 shows the effect of the cuts on data, simulated background and signal events. The agreement between data and background simulation is illustrated at preselection level in Figure 2 which shows the distributions of the main analysis variables, namely, the slow electron energy, the fitted mass of the jet system, the minimum electron isolation angle and the global event b-tagging variable x b .
The final background amounts to 6.63 ± 0.26 (stat.) +0.59 −0.93 (syst.) events, and is mainly due to e + e − qq (ZZ) events. Illustrations of the two-dimensional distribution, used as input for the confidence level computations are shown in figure 3 for data, simulated background and signal events (for m H =95 GeV/c 2 ). Table 2 shows the selection efficiency for different Higgs boson masses. The systematic uncertainties have been evaluated as described in [1]. Among the events selected in data, one has a high x b value and thus is kept in the final mass-plot, for which the supplementary cut x b > −1.8 was applied.

Muon channel
The analysis is based upon the same discriminant variables as in [1], but the selection criteria have been re-optimised, as explained in [1], to afford the best sensitivity to the expected signal at 188.7 GeV. The preselection remains unchanged except that events must now have at least nine charged particle tracks. Two muons are required with opposite charges, momenta greater than 34 GeV/c and 19 GeV/c, and an opening angle larger than 10 • . The muon isolation angle with respect to the closer jet must be greater than 16 • for the most isolated muon and greater than 8 • for the other one. A 5-C kinematic fit is then performed to test the compatibility of the di-muon mass with the Z mass. Events are kept only if the fit converges. The fitted recoiling mass is chosen as the first discriminant variable for the two-dimensional calculation of the confidence levels. The second variable is the global b-quark variable x b . Table 1 details the effect of the selections on data and simulated samples of background and signal events. The agreement of simulation with data is good. This can also be seen in Figure 4, which shows the total energy of the charged particles, the momentum of the fast muon candidate and the content in b-quark of the event after the preselection. The isolation of the muons with respect to the closest jet is also given after the lepton pair selection. At the end of the analysis, 5 events are selected in the data in good agreement with the expected background of 5.09 ± 0.19(stat.) ± 0.21(syst.) events coming mainly from ZZ. Finally the signal efficiencies for different Higgs boson masses are given in Table 2. The systematic uncertainties on background and efficiencies have been derived as explained in [1]. The two events with the largest values of x b are kept for the final mass-plot.

Higgs boson searches in events with jets and taus
Three channels are covered by this analysis, two for the SM, depending of which boson decays into τ + τ − , and one for the MSSM. Hadronic events are selected by requiring at least ten charged particles, a total reconstructed energy greater than 0.4 √ s , a reconstructed charged energy above 0.2 √ s and an effective centre-of-mass energy, √ s ′ , greater than 120 GeV. This defines the preselection.
A search for τ lepton candidates is then performed using a likelihood technique. Clusters of one or three charged particles are first preselected if they are isolated from all other particles by more than 10 • , if the cluster momentum is above 2 GeV/c and if all particles in a 10 • cone around the cluster direction make an invariant mass below 2 GeV/c 2 . The likelihood variable is calculated for the preselected clusters using distributions of the cluster momentum, of its isolation angle and of the probability that the tracks forming the cluster come from the primary vertex. Pairs of τ candidates with opposite charges and an opening angle of at least 90 • are selected using a cut on the product of their likelihoods, considering both the 1-1 (where the two τ leptons decay to one prong) and 1-3 (with at least one τ decaying to three prongs) topologies. As an example, Figure 5a shows the τ selection likelihood distribution for the selected events in the 1-1 topology.
Two slim jets are then reconstructed with all particles (charged and neutral) inside a 10 • cone around the cluster directions. The rest of the event is forced into two jets using the DURHAM algorithm. The slim jets are constrained to be in the 20 • ≤ θ τ ≤ 160 • polar angle region to reduce the Ze + e − background, while the hadronic dijet invariant mass is required to be between 20 and 110 GeV/c 2 in order to reduce the qq(γ) and Zγ * backgrounds. The jet energies and masses are then rescaled, imposing energy and momentum conservation, in order to improve the estimation of the masses of the dijets (τ + τ − and qq ). Both dijets are required to have a rescaled mass above 20 GeV/c 2 and below √ s , and each hadronic jet must have a rescaling factor in the range 0.4 to 1.5.
The remaining background comes from genuine ℓ + ℓ − qq events. In order to reject the qq e + e − and qq µ + µ − backgrounds the measured mass of the leptonic dijet is required to be between 10 and 80 GeV/c 2 and its electromagnetic energy to be below 60 GeV (see Figure 5c). The effect of the selections on data, simulated background and signal events is given in Table 1, while the selection efficiencies are summarised in Tables 2 (SM) and 3 (hA). Systematic uncertainties have been estimated by moving each selection cut according to the resolution of the corresponding variable. The main contributions are due to the τ + τ − invariant mass and electromagnetic energy. The total systematic uncertainties amount to ±6% on signal efficiencies and ±11% on the background.
At the end of the above selections, the reconstructed Higgs boson mass is estimated from the sum of the rescaled dijet masses in the hA channel and by subtracting the nominal Z mass in the hZ channels ( Figure 5e). Besides this reconstructed mass, the twodimensional calculation of the confidence levels makes use of a discriminating variable, again using a likelihood technique. This variable is built from the distributions of the rescaling factors of the τ jets, the τ momenta and the global b-tagging x b variable (see Figure 5g). Since the three possible τ + τ − qq signals are analysed in the same way, the confidence level computation uses only one global τ + τ − qq channel. At each test point, the signal expectation and distribution in this channel are obtained by summing the contributions from the three signals weighted by their expected rates.

Higgs boson searches in events with jets and missing energy
The analysis starts with a preselection which is done in two steps. The first step aims at reducing the γγ contamination and requires a total charged multiplicity greater than 10 (with at least one charged particle with a transverse momentum above 2 GeV/c), a total charged energy greater than 30 GeV and the sum of the transverse energies of the charged particles with respect to the beam axis greater than 28 GeV. The total transverse momentum has to be greater than 2 GeV/c. Furthermore, events where both the total transverse momentum and the largest single charged particle transverse momentum are less than 5 GeV/c have also been rejected. After these cuts, the γγ contamination is reduced to 1.5% of the total background, which is now dominated by qq(γ) events.
Then, jets are reconstructed from the event particles using the LUCLUS [19] algorithm with the DURHAM distance (y cut = 0.005). The results will be hereafter referred to as "free-jet clustering". Events are also forced in a two-jet topology using the same algorithm (with a result referred to as "two-jet clustering") in order to tag specifically qq(γ) events with the photon emitted along the beam axis. The rest of the preselection is designed to remove a large fraction of the remaining background without affecting the signal efficiency using the good discrimination between background and signal in the distributions of many analysis variables. The selection requires the most energetic electromagnetic cluster associated to a charged particle to be lower than 25 GeV (or 10 GeV if the charged particle associated to the cluster failed the charged particle selection criteria), the effective centreof-mass energy, √ s ′ , to be greater than 100 GeV, the absolute value of the cosine of the polar angle of the missing momentum to be lower than 0.98, the most forward jet in the free-jet clustering to be more than 16 • from the beam axis, the fraction of electromagnetic energy per jet in the free-jet clustering to be lower than 0.8, the energy deposited in the forward region within 30 • around the beam axis to be smaller than 20 GeV and the energy of the more (less) energetic jet in the two-jet clustering to be between 30 and 90 GeV (10 and 60 GeV). This defines the preselection.
The final discrimination between signal and background is achieved through a multidimensional variable built using the likelihood ratio method. The input variables are the effective centre-of-mass energy √ s ′ , the global b-tagging variable x b , the missing momentum P mis , and the cosine of its polar angle, the charged multiplicities of the jets in the free-jet clustering, the energy of the most energetic jet in the two-jet clustering, E jet1 , the acoplanarity, defined as the supplement of the angle between the transverse momenta of the two jets in the two-jet clustering, the maximal (over all particles in an event) transverse momentum with respect to the axis of the closest jet in the two-jet clustering, P tmax , and the output of a veto algorithm based on the response of the lead-scintillator counters installed at polar angles around 40 • to detect photons crossing this insensitive region of the electromagnetic calorimeters.
The distributions at preselection level of some of the input variables are shown in Figure 6 while that of the discriminant variable is given in Figure 7. Figure 8 shows the expected background rate as a function of the efficiency for a Higgs mass m h =95 GeV/c 2 when varying the cut on the discriminant variable. The final selection yields a total background of 27.8 ± 1.0(stat.) at 54.3% efficiency for the signal. The number of observed events is 27. Table 1 details the effects of the selections on data and simulated samples of background and signal events while the selection efficiencies as a function of the Higgs boson mass are summarised in Table 2.
The systematic uncertainty is dominated by the imperfect modelling of the energy flow. The corresponding error has been estimated by comparing data and simulation in test samples of Zγ events at high energy and data taken at the peak. It amounts to 11.0% relative to the background. Other sources include the imperfect modelling of b-tagging and jet angular resolutions, the dependence on the jet algorithm and uncertainties in cross-sections. This yields a total systematic error of ±12.4%.

Higgs boson searches in pure hadronic events
The aim of the four-jet preselection is to eliminate radiative and γγ events and to reduce the QCD and Zγ * background. This preselection, common to HZ and hA analyses, has not changed with respect to last year's analysis [1], except that the number of charged particles for the di-jet recognised as the Z is required to be greater than or equal to five, in order to remove events in which the Z decays into charged leptons.

The HZ four-jet channel
The present analysis is an update of the method used by DELPHI at 183 GeV [1]. Events are selected using a discriminant variable which is defined as the ratio of likelihood products for signal and background hypotheses for a set of quantities having a different behaviour in the two cases. These variables can be divided into two categories related respectively to the shape and to the b-content of the events. The six shape variables are those defined previously [1] and a new quantity: the fitted mass of the di-jet assigned to the Z. The agreement between data and background simulation is illustrated at preselection level in Figure 9 which shows the distributions of four analysis variables.
In the previous version of this analysis the Z boson mass was fixed at its nominal value and the "best" pairing of jets to select the Higgs and Z candidates was found by maximising an expression in which the b-content of the different jets and the χ 2 probability of the five-constraint fit contribute. When the production of the Higgs boson is close to the kinematic limit, the Z mass distribution no longer has a Breit-Wigner shape centred on m Z . The previous procedure has been generalised using the Z mass distribution given by the simulation for a fixed Higgs mass equal to 95 GeV/c 2 . For values of the Higgs mass which differ from 95 GeV/c 2 , the same Z mass distribution has been used, even if not optimal, in order to be independent of the assumed value for the Higgs mass. The resulting mass distribution, for preselected events in data and simulation, is shown in Figure 10.
Events originating from the signal and from the background are separated using the value of a discriminant variable. This variable combines the information from shape variables to reduce the QCD background and from the b-tagging variable to reduce the contribution of W pairs. The most effective variable against the W + W − background was found to be x i b , the combined b-tagging variable [12] measured for each jet. The likelihoods that each event is of HZ, W + W − or QCD origin are evaluated. The final discriminant variable is obtained as the ratio between signal and background likelihoods. Figure 11 shows the number of expected Monte-Carlo and observed data events, as a function of the efficiency on the signal when cutting on this discriminant variable.
The two-dimensional distribution obtained by combining the Higgs mass estimate and the likelihood ratio is used for the final limit calculation. The most background-like events are suppressed by demanding that the (log)likelihood ratio be greater than -1.0.
The systematic uncertainties have been evaluated by considering a b-tagging, a QCD related and a 4-fermion part, independently, together with an uncertainty on the crosssections of all processes, resulting in an overall relative uncertainty of 7.5% at the final selection level.

The hA four-b channel
A likelihood method has been applied to search for hA production in the four-jet channel. After the common four-jet preselection, tighter cuts were applied to the remaining events, namely, a cut in the parameter of the DURHAM algorithm (y cut ≥ 0.003) is imposed as well as the requirement of at least two charged particles per jet. Finally, an event is rejected if its maximum inter-jet energy difference is greater than 70 GeV. The resulting number of expected events and the signal efficiencies after this preselection are given in Table 1.
The following eight variables are combined in the likelihood: the event thrust, the second and fourth Fox-Wolfram moments, H 2 and H 4 , the minimal (among the three possible pairings of jets) di-jet-masses difference, the production angle of the candidate bosons, the sum (over the four jets) of the b-tag jet variable, the minimum di-jet btag and the number of secondary vertices. For each event, the measured value of each of these discriminant variables is compared with probability density functions obtained from simulated events. Figure 12 shows the resulting efficiency versus the total background varying the likelihood cut (assuming m A = 80 GeV/c 2 and tan β = 20). The final cut value is chosen depending on the efficiency-background point desired; for the derivation of the limits, a cut on the likelihood output at 2.0 has been chosen. A total of 13 events is observed while 11.41 are expected. A more stringent cut on the likelihood output corresponding to a requirement at 3.05 is used for the mass plot (see Fig 13). This yields a total background of 4.72 ± 0.22 coming from qq(γ)( 2.11 ± 0.14) and 4-fermion processes (2.61 ± 0.17) and 3 events selected in data with a sum of their di-jet masses of 147, 180 and 176 GeV/c 2 . Efficiencies obtained for different masses and tan β are summarised in Table 3.
To check systematic uncertainties on the total background due to the modeling of the shape of the probability density functions, the training and validation sample were exchanged and the analysis repeated. The uncertainty on the total background due to this effect has been estimated at the level of 5% and has been added (quadratically) to the other sources of errors, in particular the one coming from the b-tagging estimation, resulting in a final relative systematic uncertainty of 8%.

Summary and results
The results of the searches presented in the previous sections can be translated into exclusion limits on the masses of the neutral Higgs bosons in the SM and MSSM.

Summary
For each analysis of the HZ and hA channels at 188.7 GeV, the integrated luminosity, the expected backgrounds and their errors, and the number of observed events at various levels of the analyses are summarised in Table 1

The SM Higgs boson
As an illustration, Figure 14 shows the distribution of the reconstructed Higgs boson mass found in the HZ channel after the tight selection (last lines in Table 1) for data, simulated background and signal events. The last cut in each channel has been chosen such that the signal-over-background ratio (for the reference mass) be almost equal for all channels (between .2 and .35). Within the 188.7 GeV data, the total number of events observed in all channels is 31, which is consistent with the 35.9 expected background. Since the two hypotheses (background only and background plus signal at 95 GeV/c 2 ) are almost indistinguishable, a possible signal at 90 GeV/c 2 has been superimposed in order to visualise our resolution in the mass of such a signal.
We proceed to set a limit on the SM Higgs boson mass, combining these data with those taken at lower energies, namely 161,172 GeV [14] and 183 GeV [1]. The expected cross-sections and branching ratios are taken from [20,21], with the top mass set to 175 GeV/c 2 .
The confidence levels CL b , CL sb and CL s are computed as described in [1]. CL b and CL sb are the confidence levels in these hypotheses (background only and signal + background), while CL s is conservatively taken as their ratio (CL sb /CL b ).
In the presence of a sizeable Higgs signal the value of the observed CL b (top of Figure 15) would approach one, because it measures the fraction of background-only experiments which are less signal-like than the observation. On the contrary here, the observation agrees well with the expectation (background only). Furthermore the curve for the signal hypothesis shows that the expected 5σ discovery limit (horizontal line at 1 − CL b = 5.7 × 10 −7 ) is at 88.6 GeV/c 2 . The confidence level in the signal is shown in Fig. 15 (bottom). The observed 95% CL lower limit on the mass is m H > 94.6 GeV/c 2 , while the expected mean is 94.4 GeV/c 2 and the expected median (50% exclusion potential) is 95.3 GeV/c 2 . If errors had not been allowed for, the observed (expected) limit would have been increased by 0.2 GeV/c 2 (0.4 GeV/c 2 ).
The effective ∆χ 2 (−2∆ ln L) with which the SM Higgs is excluded is shown in Fig. 16. In the event of a discovery the ∆χ 2 would be negative, and could be used to extract the mass and its error, as can be seen on the bottom plot of Fig. 16.
Finally the data can be used to set 95% CL upper bounds on the HZZ coupling in non-standard models which assume that the Higgs boson decay properties are identical to those in the SM but the cross-section may be different. Figure 17 shows the excluded cross-section as a function of the test mass.

Neutral Higgs bosons in the MSSM
The results in the hZ and hA processes are combined with the same statistical method as for the SM, using also earlier results at LEP2 energies [1,14,22].
In the MSSM, at tree level, the production cross-sections and the Higgs branching fractions depend on two free parameters, tan β and one Higgs boson mass, or, alternatively, two Higgs boson masses, eg m A and m h . The properties of the MSSM Higgs bosons are modified by radiative corrections which introduce additional parameters: the mass of the top quark, the Higgs mixing parameter, µ, the common sfermion mass term at the EW scale, M S , the common SU(2) gaugino mass term 1 at the EW scale, M 2 , and the common squark trilinear coupling at the EW scale, A. The interpretation of the experimental results depends on the values assumed for these parameters as well as on the order of the calculated radiative corrections.
The results described hereafter rely on leading-order two-loop calculations of the radiative corrections in the renormalization group approach [23], with recent modifications (about top threshold and gluino two-loop corrections) that make the computations agree with fully diagrammatic two-loop calculations [24]. After these improvements, the benchmark prescriptions for the parameters beyond tree-level have also been refined, leading to two extreme scenarii for the theoretical upper bound on m h as a function of tan β [25] which differ only by the value of X t = A − µ cot β, the parameter which controls the mixing in the stop sector. In the following, we adopt these new prescriptions which correspond to: 175 GeV/c 2 for the top mass, 1 TeV/c 2 for M S , 200 GeV/c 2 for M 2 and -200 GeV/c 2 for µ. Two values have been considered for the mixing in the stop sector: The signal expectations in each channel are derived from the cross-sections, the experimental luminosity and the efficiencies. A correction is applied to account for differing branching fractions of the Higgs bosons into bb and τ + τ − between the input point and the simulation (e.g. for the hZ process, the simulation is done in the SM framework). For the hA channels, as there can be a difference between the masses of the h and A bosons at low tan β, the set of hA efficiencies obtained from the simulation at tan β = 20 is applied at all points with tan β above 2.5, while the set of efficiencies derived from the tan β = 2 simulation is applied below. The same holds for the discriminant information. The signal expectations, expected backgrounds and numbers of candidates enter in the computation of the observed confidence level in the signal hypothesis at the input point, CL s . The expected confidence level in the signal hypothesis is also derived at each point. As there is a large overlap in the background selected by the two four-jet channels, only one channel is selected at each input point, on the basis of the best signal over background ratio. This ensures that the channels which are combined in the confidence level computations are independent.
The results translate into regions of the MSSM parameter space excluded at 95% CL or more. The excluded regions are presented in the (m h , tan β) plane in Fig. 18, in the (m A , tan β) plane in Fig. 19 and in the (m A , m h ) plane in Fig. 20. As illustrated in the latter, there is a small region of the parameter space where the decay h→AA opens, in which case it supplants the h→bb decay. But, due to the high luminosity collected at 188.7 GeV, the results in the h→bb channel alone cover most of the area which remained unexcluded at 183 GeV in this region [1].
Finally, the results shown in Figs The expected limits are 81.3 GeV/c 2 in m h and 82.3 GeV/c 2 in m A . In the low tan β region, in the no mixing case, all values of m A up to 1 TeV/c 2 are excluded, providing an excluded range in tan β between 0.6 and 2.2, in agreement with the expected excluded range. On the other hand, no limit can be set on tan β in the m max h scenario (see Fig. 19).

Interpretation in a general Two Higgs Doublet Model
These results can also be translated to the framework of a general Two Higgs Doublet Model (2HDM) with one assumption (the decay of both h and A is dominated by bb and/or cc final states) and two options: CP conserving or CP violating.
In the CP-conserving two Higgs doublet model, the h and H bosons are mixtures of the real parts of the neutral Higgs fields, while A derives from the imaginary components not absorbed by the Z. The coupling strengths are: C hZ = sin (α−β) and C hA = cos (α−β). These couplings clearly indicate the complementarity of the two processes. Besides, if one of them is experimentally out of reach for a given set of masses, no exclusion is possible since mixing angles could always be such that the other process is suppressed below detectability.
The exclusion plot (left of Figure 21) is obtained in the following way: for each pair of m h and m A values the number of expected events for the channels hA and hZ is calculated using the cross-sections, integrated luminosities, branching ratios and efficiencies quoted in this paper. This number depends obviously on the factors C hZ and C hA . Thus a minimisation of the confidence level with respect to these factors has been performed, taking into account the sum rule C 2 hZ + C 2 hA = 1. It should be noted that in 2HDM the branching ratios of A and h into bb are proportional to tan β and sin α/ cos β respectively. Imposing the condition | sin α| > cos β, which is barely restrictive for medium or large values of tan β, leads to a dominant coupling to b quarks for both h and A bosons, (Zone I in Figure 21).
The case of non-b decays has also been studied, using the selections of this paper except those referring to b-tagging. This takes care of a scenario with tan β< 1 which would allow a dominant decay of the Higgs boson into cc [26]. This region (named Zone II in Figure 21) occurs for values of the α and β angles such that α ∼ β ∼ 0, resulting in a very restrictive and particular parameter set of the 2HDM which represents the most pessimistic scenario in this kind of search. Any other situation will lead to an intermediate excluded region, as for example when sin α = 0 but tan β > 1 which implies qq bb decays of the hA signal.
CP violation in the SUSY sector is an open possibility and may even be necessary in the electroweak baryogenesis scenario (see for instance [27] and references therein). The violation leads to three neutral Higgs bosons (noted h 1 , h 2 and h 3 , sorted by mass) with undefined CP properties. It is shown in [28] that the previous sum rule, valid in the CP-conserving case, can be extended to the CP-violating model, giving the relation: C 2 h 1 Z +C 2 h 2 Z +C 2 h 1 h 2 = 1, which together with the results described in this paper allows the minimum number of events for each point (m h 1 , m h 2 ) to be calculated, with an extended procedure similar to the one used previously, leading to the excluded region at 95% CL shown in the lower plot of Figure 21 (only the conditions for Zone I were used in this case).
This analysis shows that it is possible to exclude a large region of Higgs boson masses even when relaxing the standard assumptions made in the MSSM scenario.

Conclusions
The 158 pb −1 of data taken by DELPHI at 188.7 GeV, combined with our lower energy data, sets the lower limit at 95% CL on the mass of the Standard Model Higgs boson at: m H > 94. 6 GeV/c 2 .
The MSSM studies described above give for all values of tan β above 0.6, and assuming m A > 20 GeV/c 2 : Beyond the results described above, DELPHI performed a more complete scan of the parameters of the MSSM: this analysis, described in the addendum, shows clearly the robustness of the limits obtained in the benchmark scenarios.
Other LEP experiments, using their data sets collected concurrently with the ones used in this work, have reported similar results [29,30].