Photonic events with missing energy in e+e- collisions at sqrt(s)=189 GeV

Photonic events with large missing energy have been observed in e+e- collisions at a centre-of-mass energy of 189GeV using the OPAL detector at LEP. Results are presented for event topologies consistent with a single photon or with an acoplanar photon pair. Cross-section measurements are performed within the kinematic acceptance of each selection, and the number of light neutrino species is measured. Cross-section results are compared with the expectations from the Standard Model process e+e- to nu nubar + photon(s). No evidence is observed for new physics contributions to these final states. Upper limits are derived on sigma(e+e- to XY).BR(X to Y gamma) and sigma(e+e- to XX).BR**2(X to Y gamma) for the case of stable and invisible Y. These limits apply to single and pair production of excited neutrinos (X=nu*, Y = nu), to neutralino production (X=neutralino_2, Y=neutralino_1) and to supersymmetric models in which X = neutralino_1 and Y = light gravitino. The case of macroscopic decay lengths of particle X is considered for e+e- to XX, X to Y gamma, when M_Y is of order zero. The single-photon results are also used to place upper limits on superlight gravitino pair production as well as graviton-photon production in the context of theories with additional space dimensions.


Introduction
We describe measurements and searches using a data sample of photonic events with large missing energy collected in 1998 with the OPAL detector at LEP.The events result from e + e − collisions at a centre-of-mass energy of 188.6 GeV with an integrated luminosity of 177.3 pb −1 .The present paper builds on publications from earlier data samples at lower centre-of-mass energies [1,2].This data-set at 189 GeV gives discovery potential in a new kinematic regime with about a four-fold increase in integrated luminosity.Measurements of photonic event production have also been made by the other LEP collaborations at centre-of-mass energies above the W pair threshold [3], including new results from L3 and DELPHI at √ s = 189 GeV [4].
The single-photon and acoplanar-photons search topologies presented here are designed to select events with one or more photons and significant missing transverse energy, indicating the presence of at least one neutrino-like invisible particle which interacts only weakly with matter.The event selections for these search topologies are similar to those used in our recent publication [1].The single-photon search topology is sensitive to events in which there are one or two photons and missing energy which, within the Standard Model, are expected from the e + e − → ννγ(γ) process 1 .The acoplanar-photons search topology is designed to select events with two or more photons and significant missing transverse energy which, within the Standard Model, are expected from the e + e − → ννγγ(γ) process.
The single photon topology provides a direct measurement of the invisible width of the Z 0 and can probe charged and neutral triple gauge couplings.The acoplanar-photons topology can probe WWγγ quartic couplings in the e + e − → ν e ν e γγ process.The neutral and quartic gauge coupling measurements will be described in forthcoming papers based on the event selections described herein.
These photonic final-state topologies are sensitive to several different new physics scenarios.A generic classification is e + e − → XY or e + e − → XX where X is neutral and can decay radiatively (X → Yγ) and Y is stable and only weakly interacting.For the general case of massive X and Y this includes conventional supersymmetric processes (X = χ0 2 , Y = χ0 1 ).These topologies also have particularly good sensitivity for the special case of M Y ≈ 0. This applies both to the production of excited neutrinos (X = ν * , Y = ν) and to supersymmetric models in which the lightest supersymmetric particle (LSP) is a light gravitino and χ0 1 is the next-to-lightest supersymmetric particle (NLSP) which decays to a gravitino and a photon (X = χ0 1 , Y = G).The neutralino life-time in such models is a free parameter and so we also address the possibility of neutralino-pair production with macroscopic decay lengths.One type of new physics which could be seen in the single-photon topology is the production of an invisible particle in association with a photon.An example of this is graviton-photon production, e + e − → Gγ [5,6].This can occur within string theory models which allow gravitons to propagate in a higher-dimensional space but restrict Standard Model particles to the usual four space-time dimensions [7].Another type of new physics is the production of invisible particles tagged by initial-state radiation.One example2 is production of a pair of gravitinos, e + e − → G Gγ, as in the superlight gravitino model [8].The acoplanar-photons search topology also has sensitivity to the production of two particles, one invisible, or with an invisible decay mode, and the other decaying into two photons.This paper will first describe the OPAL detector and the Monte Carlo samples used.A brief summary of the event selections will then be given, followed by cross-section measurements for e + e − → ννγ(γ) and e + e − → ννγγ(γ) and comparisons with Standard Model expectations.The new physics search results will then be discussed.

Detector and Monte Carlo Samples
The OPAL detector, which is described in detail in [9], contains a silicon micro-vertex detector surrounded by a pressurized central tracking system operating inside a solenoid with a magnetic field of 0.435 T. The region outside the solenoid (barrel) and the pressure bell (endcap) is instrumented with scintillation counters, presamplers and the lead-glass electromagnetic calorimeter (ECAL).The magnet return yoke is instrumented for hadron calorimetry and is surrounded by external muon chambers.Electromagnetic calorimeters close to the beam axis measure luminosity and complete the acceptance.
The measurements presented here are mainly based on the observation of clusters of energy deposited in the lead-glass electromagnetic calorimeter.This consists of an array of 9,440 lead-glass blocks in the barrel (| cos θ| < 0.82) with a quasi-pointing geometry and two dome-shaped endcap arrays, each of 1,132 lead-glass blocks, covering the polar angle3 range (0.81 < | cos θ| < 0.984).Fully hermetic electromagnetic calorimeter coverage is achieved beyond the end of the ECAL down to small polar angles with the use of the the gamma-catcher calorimeter, the forward calorimeter (FD) and the silicon-tungsten calorimeter (SW).
Scintillators in the barrel and endcap regions are used to reject backgrounds from cosmic ray interactions by providing time measurements for the large fraction (≈ 80%) of photons which convert in the material in front of the ECAL.The barrel time-of-flight (TOF) scintillator bars are located outside the solenoid in front of the barrel ECAL and match its geometrical acceptance (| cos θ| < 0.82).Tile endcap (TE) scintillator arrays [10] are located at 0.81 < | cos θ | < 0.955 behind the pressure bell and in front of the endcap ECAL.
The integrated luminosities of the data samples are determined to better than 1% from smallangle Bhabha scattering events in the SW calorimeter.Triggers [11] based on electromagnetic energy deposits in either the barrel or endcap electromagnetic calorimeters lead to full trigger efficiency for photonic events passing the event selection criteria described in the following section.
To simulate possible new physics processes of the type e + e − → XY and e + e − → XX where X decays to Yγ and Y escapes detection, a modified version of the SUSYGEN [20] Monte Carlo generator was used to produce neutralino pair events of the type e + e − → χ0 2 χ0 1 and e + e − → χ0 2 χ0 2 , χ0 2 → χ0 1 γ, with isotropic angular distributions for the production and decay of χ0 2 and including initial-state radiation.Monte Carlo events were generated at 48 (for XY production) and 42 (for XX production) points in the kinematically accessible region of the (M X , M Y ) plane.A Monte Carlo generator was written to simulate the superlight gravitino signature G Gγ, discussed in Section 5.1.4.The KORALZ ννγ(γ) sample was used to determine the efficiency for graviton-photon production in the context of additional space dimensions, by means of an event reweighting (see section 5.1.5).The same procedure was also used to calculate the efficiency for G Gγ production, and compared with the prediction of the direct Monte Carlo simulation.All the Monte Carlo samples described above were processed through the OPAL detector simulation [21].
Simulation of e + e − → XX, X → Yγ signal events with M Y ≈ 0 where M X has a finite nonzero lifetime τ X was implemented in the full simulation of the OPAL detector.In particular, the massive quasi-stable neutral particle X was assigned properties similar to a heavy neutrino for the purpose of propagation, and then propagated within the GEANT [22] framework, according to its initial kinematics and τ X before forcing the X → Yγ decay.

Photonic Event Selection
This section summarizes the criteria for selecting single-photon and acoplanar-photons events.The kinematic acceptance of each selection is defined in terms of the photon energy, E γ , and the photon polar angle, θ.In addition, the scaled energy, x γ , is defined as E γ /E beam , and the scaled transverse energy, x T , as x γ sin θ.
• The transverse momentum p γγ T of the two-photon system consisting of the two highest energy photons must satisfy p γγ T /E beam > 0.05.
In each of the two cases we retain acceptance for events with additional photons in which the resulting photonic system is still consistent with the presence of significant missing energy.This reduces the sensitivity of each measurement to the modelling of higher-order contributions.Consequently, a large fraction of the kinematic acceptance of the acoplanar-photons selection is also contained in the kinematic acceptance of the single-photon selection.

Single-Photon Event Selection
The single-photon selection criteria are explained in detail in a previous publication [1].After defining the kinematic acceptance, additional cuts on cluster quality, forward energy, muon chamber and hadron calorimeter information, and multiphoton kinematics are used to remove cosmic ray backgrounds, beam related backgrounds and standard model physics backgrounds that have no missing energy.Then events are classed as either having or not having photon conversion candidates based on tracking chamber information.Each of these two classes has somewhat different additional selection criteria based on charged track activity, time-of-flight information, and other background suppression cuts.
Several improvements have been made to the single-photon selection used in this analysis compared with the previous analysis [1].They are described in detail in Appendix A.

Acoplanar-Photons Event Selection
The acoplanar-photons selection has classes of selection requirements similar to those of the singlephoton selection.The details of the acoplanar-photons selection are described in our previous analysis [1]; changes with respect to the previous analysis are minimal and are given in Appendix A.
The requirement of a second photon in an event reduces many of the backgrounds which are otherwise a problem for the single-photon analysis.For this reason, the acoplanar photons selection can tolerate a lower energy threshold for the most energetic photon, as well as looser but more inclusive acceptance requirements for photon conversions.On the other hand, the single-photon selection has more acceptance for events having no time-of-flight information for the photons.In order to obtain the best overall acceptance for acoplanar-photons, we have added to the acoplanar-photons selection that part of the single-photon selection which contains two photons within the kinematic acceptance of the acoplanar-photons selection.This addition results in a relative increase in efficiency of 9.6% for Standard Model e + e − → ννγγ(γ) events.

Selection Results
The results of the single-photon and acoplanar-photons selections and the corresponding cross-section measurements and other measured event quantities are given below in sections 4.1 and 4.2.Results from both selections are summarized in Table 1.

Single-Photon
After applying the single-photon selection criteria to the data sample, 643 events are selected.The expected contribution from cosmic ray and beam-related backgrounds is 4.6 ± 1.5 events.These backgrounds have been estimated from events having out-of-time TOF or TE information, but passing all other selection criteria, and from events selected with looser criteria that have been visually scanned.Of the expected physics backgrounds from plausible sources, e + e − → e + e − γ, e + e − → ℓ + ℓ − ν ν(γ), e + e − → e + e − ℓ + ℓ − , e + e − → ν νqq, e + e − → µ + µ − γ and e + e − → τ + τ − γ have non-negligible contributions 4 .The total number of expected physics background events is 4.6 ± 0.5 and the contributing sources are summarized in Table 2.The number of events expected from the Standard Model process e + e − → ννγ(γ) as predicted by KORALZ is 679 ± 5 (stat) ± 14 (sys) where the systematic error is from experimental sources as discussed later.The number of events observed is consistent with the number expected from e + e − → ννγ(γ) plus the background.The efficiency for selecting e + e − → ννγ(γ) events within the kinematic acceptance of the single-photon selection is (82.1 ± 1.7)%.For both the single-photon and acoplanar-photons selections, efficiency losses due to vetoes on random detector occupancy range from about (3-5)%.Quoted efficiencies include these losses.The single-photon cross-section for √ s = 189 GeV, accounting for detector and selection efficiencies and subtracting the estimated background, is 4.35 ± 0.17 (stat) ± 0.09 (sys) pb.
The systematic error on the cross-section measurement from knowledge of the efficiency and normalisation is estimated to be 2.1%.The contributing uncertainties are summarized in Table 3.The dominant systematics are the selection efficiency uncertainty (1.5%) and the uncertainty on the detector occupancy estimate (1%).The event selection efficiency is controlled using a data sample of around 1200 e + e − → γγ events.In particular, the efficiency and time response of the TB and TE scintillators are measured from these data and small correction factors applied to the efficiency.The selection efficiency systematic includes the uncertainties on such corrections and an estimate of the residual uncertainty on the efficiency.The occupancy, estimated from random beam-crossing triggers, is typically 4% and we assign an error of 1%.We have compared the estimated efficiency using the NUNUGPV98 and KORALZ event generators.We find that the relative efficiency difference is (0.5 ± 0.4)% and assign an error of 0.5% to take into account the sensitivity to the modelling of the photon energy, angle and multiplicity distributions.This source of error is considerably reduced from our previous publications, largely because both generators that are now used are designed to model accurately ν νγγ events.Additional systematic errors on the efficiency arise from uncertainties in modelling the material close to or inside the beam-pipe which accounts for photons which convert early (20% relative uncertainty) and from systematic effects in the track reconstruction and parameter estimation which affect the association of tracks to ECAL clusters.Uncertainties in the photon energy scale, resolution and angular measurement give small additional systematic contributions.
The cross-section as a function of centre-of-mass energy is plotted in Figure 1.Cross-section results from the current analysis and from our earlier publications [1,2] are plotted.The curve shows the predicted cross-section from the KORALZ event generator for the Standard Model process e + e − → ννγ(γ).The data are consistent with the prediction.
Figure 2a shows the recoil mass distribution, where the recoil mass M recoil is defined as the mass recoiling against the photon (or against the two-photon system).The peak in the distribution at M Z is due to a large contribution from the decay Z 0 → νν. Figure 2b shows the polar angle distribution along with the e + e − → ννγ(γ) Monte Carlo expectation.In both distributions, there is consistency between data and Monte Carlo, although we note a deficit in the radiative return peak and a slight excess in the low energy, high recoil mass region.
The single-photon selection is designed to allow for the presence of a second photon in order to accept events from the e + e − → ννγγ process.Thirty-six observed events are considered to be two-photon events (i.e. have a second photon with deposited energy exceeding 300 MeV and with 15 • < θ < 165 • ), consistent with the expectation of 33.6 ± 1.5 events from the KORALZ Monte Carlo.Of these, 20 fall within the kinematic acceptance of the acoplanar photon selection, as compared to the KORALZ expectation of 23.5 ± 1.0.

Acoplanar-Photons
The acoplanar-photons selection applied to the data sample yields 24 events, in good agreement with the KORALZ prediction of 26.9 ± 1.2 events for the Standard Model e + e − → ννγγ(γ) contribution.The expected contribution from other Standard Model processes and from cosmic ray and beam related backgrounds is 0.11 ± 0.04 events.The selection efficiency for e + e − → ννγγ(γ) events within the kinematic acceptance of the selection is (66.4 ± 2.9)%.The corresponding cross-section is 0.204 ± 0.043 pb, compared to a KORALZ prediction of 0.228 ± 0.002 pb.The OPAL measurements of the cross-sections at √ s = 130, 136, 161, 172, 183 and 189 GeV are summarized in Table 4. Results for √ s < 189 GeV have been taken from our previous publications [1,2].
The dominant source of systematic uncertainties is modelling of the reconstruction efficiency, especially the simulation of the detector material and consequent photon conversion probabilities.Other sources arise from uncertainties on the electromagnetic calorimeter energy scale and resolution, on the integrated luminosity measurement, on detector occupancy estimates and from comparisons of different Monte Carlo event generators for the process e + e − → ννγγ(γ).The relative systematic error from all sources is 4.3%.
The kinematic properties of the selected events are displayed in Figure 3 where they are compared with the predicted distributions for e + e − → ννγγ(γ) obtained using the KORALZ generator normalized to the integrated luminosity of the data.Plot (a) shows the recoil mass distribution of the selected acoplanar-photon pairs (or of the two most-energetic photons in the case of events with three or more photons).The distribution is peaked near the mass of the Z 0 as is expected for contributions from e + e − → ννγγ(γ).The resolution of the recoil mass is typically 3-5 GeV for M recoil ≈ M Z .Plot (b) shows the distribution of the scaled energy of the second most energetic photon.Plot (c) shows the γγ invariant-mass distribution for which the mass resolution is typically 0.7-1.9GeV.Plot (d) shows the distribution in scaled transverse momentum of the selected two-photon system.There is 1 selected event having a third photon with deposited energy above 300 MeV and within the polar-angle acceptance of the selection.The corresponding expectation from KORALZ is 1.20 ± 0.08 events.

Data Interpretation
The results of the single-photon and acoplanar-photons selections are used to test the Standard Model and search for new physics contributions.
For the XY and XX searches, we test the range of the following product branching ratios which are consistent with the data: σ(e + e − → XY) • BR(X → Yγ) and σ(e + e − → XX) • BR 2 (X → Yγ).Given lack of evidence for signal, we then set 95% CL upper limits on these quantities.This is done both for the general case of massive X and Y, and also separately for the special case of M Y ≈ 0. All efficiencies are first evaluated under the assumption that the decay length of X is zero.For the XX search in the special case of M Y ≈ 0, we evaluate the efficiency as a function of the non-zero lifetime of the X particle, thus quantifying the sensitivity of the search for the general case of non-prompt decay.For both the XY and XX searches, Monte Carlo samples were generated for a variety of mass points in the kinematically accessible region of the (M X , M Y ) plane.To set limits for arbitrary M X and M Y , the efficiency over the entire (M X , M Y ) plane is parameterized using the efficiencies calculated at the generated mass points.As justified previously [1], we restrict the searches to M X + M Y > M Z for the single-photon topology and to M X values larger than about M Z /2 for the acoplanar-photons topology.

Single-Photon
With the single-photon topology we measure the number of light neutrino species and observe a rise in the cross-section at low photon energies which is consistent with the additional cross-section expected from charged current contributions to e + e − → ν e ν e γ(γ).We also give results of searches for XY production, superlight gravitino pair production, and graviton-photon production in the context of models with extra space dimensions.

Neutrino Counting
Single-photon events are expected within the Standard Model from the process e + e − → ννγ(γ).At the tree level, the cross-section for muon-type and tau-type neutrinos is attributable solely to s-channel Z production with initial-state radiation.For electron-type neutrinos, the single photon cross-section arises from Feynman diagrams corresponding to s-channel Z production with initial-state radiation and t-channel W exchange with radiation from the initial-state or the exchanged W. Higher order electroweak processes such as WW boxes are expected to give a negligible contribution to the crosssection.
The results of the single-photon selection are used to measure the size of the s-channel Z production contributions and the W-related contributions (W amplitude squared plus the W-Z interference), which we parameterize in terms of N ν and f W .The measurement of the pure s-channel Z production contributions is a direct measurement of the Z invisible width, which is related to the effective number of light neutrino generations, N ν , defined as the ratio of the Z invisible width to the expected width in the Standard Model for one neutrino generation.The W-related contributions are parameterized by a multiplicative scale factor, f W , defined to be 1 for the Standard Model expectation.
In order to measure N ν and/or f W , we perform a binned fit to both the overall event rate and to the shape of the photon energy distribution by minimising the negative log-likelihood: where n exp is given by Here n obs is the observed number of events, n exp is the expected number of events as a function of N ν and f W derived by reweighting fully simulated ννγ(γ) events, and P represents the Poisson probability for observing n obs events given an expectation of n exp .n W exp is the number of events for Standard Model W contributions 5 and n Z exp is the number of s-channel events expected per neutrino generation; both n W exp and n Z exp were evaluated using the NUNUGPV98 generator6 .
The dominant systematic errors arise from the uncertainty on the selection efficiency and from the theoretical uncertainty on the expected number of events, with a minor contribution from the limited Monte Carlo statistics.The relative uncertainty on the selection efficiency is 2.1% and was discussed in Section 4. We assign a value of 2% to the theoretical uncertainty based on comparisons and estimated precisions of the NUNUGPV98, KORALZ and grcννγ [23] event generators.In addition, an uncertainty related to the modelling of the photon energy spectrum is assigned, again estimated by comparing the NUNUGPV98 and KORALZ event generators.Other sources of systematic error, such as the uncertainty on the centre-of-mass energy and the expected background, result in negligible contributions.
In a first step, we check the consistency with the Standard Model predictions by fitting for both N ν and f W .The results are f W = 1.12 ± 0.13 (stat) ± 0.12 (sys) and N ν = 2.63 ± 0.15 (stat) ± 0.11 (sys) with a correlation coefficient of −41%. Figure 4 shows the 70%, 95% and 99% probability contours in the space of the two parameters, while Figure 5 shows the photon energy distribution for the events selected at 189 GeV, compared to the expectation with f W = 1.12,N ν = 2.63.The data are seen to be in fair agreement with the Standard Model prediction for N ν = 3 and f W = 1.In particular, the W contributions are observed with a high degree of significance, even when allowing N ν to be unconstrained.
In a second step, we assume the W contributions to be as predicted by the Standard Model and fit for N ν .The result is: N ν = 2.69 ± 0.13 (stat) ± 0.11 (sys).
Alternatively, one can assume N ν = 3 as indicated by the precise, but less direct, measurements of the Z lineshape [24], and fit for the relative size of the W-contributions.The result is f W = 0.99 ± 0.11 (stat) ± 0.12 (sys), establishing that the W-contributions are observed and are consistent with the expectations from the Standard Model.We search for evidence of new physics processes of type e + e − → XY, X → Yγ.To select candidate new physics events, kinematic-consistency and degraded-resolution cuts are applied to events in the singlephoton event sample.The kinematic-consistency cuts require that the energy of the most energetic photon be within the range kinematically consistent with mass values (M X , M Y ) after accounting for energy resolution effects.The degraded-resolution cuts reject events in which the most energetic photon is in the angular regions 0.78 < | cos θ| < 0.82 or | cos θ| > 0.95; energy resolution is significantly degraded in these regions and the Monte Carlo simulation is less reliable.These cuts are approximately 94% efficient in selecting signal events within the kinematic acceptance of the single-photon selection, assuming uniform energy and cos θ signal distributions.
Two methods are used to search for, and place upper limits on, contributions from e + e − → XY, X → Yγ.The first is an event-counting method, based simply on the total number of events selected as candidates.The kinematic cuts, described above, are sufficiently loose that the results are relatively insensitive to the shapes of signal distributions in photon energy or cos θ that result from the specifics of a particular model.Thus this method gives the more generally applicable crosssection limits for XY production.In many models, however, the distributions for the production and decay angles in the process e + e − → XY, X → Yγ are approximately isotropic.We therefore also perform a likelihood-based analysis to search for new physics contributions under the assumption of isotropically distributed production and decay angles.This results in a substantial increase in experimental sensitivity, since energy and angular distributions can now be used to differentiate signal from background.
In both methods, Standard Model background is assumed to be from e + e − → ννγ(γ) only.Background from other sources, including the estimated cosmic ray and beam-related background, is small and is neglected in the limit calculations.Since uncertainties due to energy scale and resolution effects in the low M X − M Y region lead to large relative uncertainties in the estimated efficiencies, we restrict our search to the region in which M X − M Y > 5 GeV.

Event-Counting Method:
Upper limits are calculated based on the total numbers of observed and expected background events using the method described in reference [25].This procedure is similar to the method used in our previous publication [1], with the exception that only the kinematic-consistency and degraded-resolution cuts are applied in order to improve the generality of the results.
The total signal efficiency within the single-photon kinematic acceptance is approximately 85%, varying by less than 1% over most of the (M X , M Y ) plane.However, the efficiency decreases significantly when the mass difference M X − M Y becomes small, due to low photon energies resulting from the X → Yγ decay.Contributions to the systematic error on the efficiency for selecting events from potential new physics sources are similar in nature and size to those discussed in our previous publication.The total relative systematic error varies from 2.5 to 5.5%, depending on M X and M Y .These systematic errors have been treated according to the method in reference [26]; the effect on the upper limits is small.Uncertainties in the ννγ(γ) background estimate have a more significant effect on the upper limits, and are treated by means of a convolution within the limit calculations.A relative 3% uncertainty has been assigned to the expected background contribution, based on factors considered in the cross-section measurement and assigning a 2% theoretical uncertainty.
There are 552 events selected after the degraded-resolution cuts are applied; the KORALZ expectation is 601 ± 14.The photon energy distribution for these events is shown in Figure 6.The number of selected events after further application of the kinematic-consistency cuts depends on the values for M X and M Y ; ranges for observed and expected events are displayed as contour plots in Figure 7.The consistency between the number of observed events (n obs ) and the number of expected background events (n b ) is given in Figure 8. Plotted is P fluct , the probability for observing at least n obs events given a background expectation of n b ± σ b events, defined by: This figure is intended to highlight regions in the (M X , M Y ) plane with an excess of events that could potentially be indicative of new physics.Such an excess would appear as a very low probability (P fluct ≪ 1%).However, a deficit of events is actually observed for most of the (M X , M Y ) plane.This leads to high values for P fluct ; for example, for n obs ≫ 1, a value P fluct = 0.97 indicates a deficit consistent at approximately the 3% level.The calculated values of P fluct range from 0.71% to 99.94%.
Values with rather low and high probability occur partly because we consider several thousand points in the (M X , M Y ) plane; each of these has a different set of kinematic-consistency cuts which select different parts of the photon energy spectrum.The net effect is that almost every possible energy range is selected for some (M X , M Y ) value.For example, as seen in Figure 6, there is a deficit in the energy range from approximately 26 to 38 GeV; this range corresponds roughly to M X = 113 GeV and M Y = 74 GeV.
The general structure of Figure 8 can also be understood from the photon energy spectrum.There is an overall deficit of approximately 49 events (∼ 1.7σ).As seen in Figure 6, most of the deficit lies in the region of the radiative return peak, while there is an excess in the low-energy region.For regions in the (M X , M Y ) plane for which low-energy photons are not kinematically allowed, events populating the low-energy region of the energy spectrum are rejected and the significance of the remaining deficit is enhanced.In particular, this occurs when M X − M Y is relatively large, and accounts for the large P fluct values in the low-M Y region of Figure 8. Conversely, regions with a small mass difference M X − M Y correspond to low photon energies, which in turn correspond to the region of the photon energy spectrum in which there is an excess.This accounts for the small P fluct values in the low M X − M Y region of Figure 8.
Cross-section times branching-ratio upper limits calculated using the event-counting method are shown in Figure 9 and range from 43 to 409 fb.Results for one M X , M Y pair (M X = 125, M Y = 63 GeV) are shown as an addition to the expected ννγ(γ) background in Figure 6.The added signal contribution is that which would be expected from a cross-section equal to the 95% CL upper limit.

Likelihood Method:
The likelihood-based analysis is a straightforward extension of an extended maximum likelihood fit [27].Upper limits are calculated using information from the photon energy and angular distributions as well as the total number of observed events.The number of observed events, number of expected background events, and signal efficiencies are the same as in the event-counting method.
The likelihood function is given by where the product is over all selected events in the data, and • n b = expected number of background events, • n s = ǫ s σ s L = expected number of signal events, with ǫ s the efficiency for observing signal events, σ s the signal cross-section times branching ratio, and L the integrated luminosity, • n obs = number of candidate events observed in the data, • P((n b + n s ) → n obs ) = e −(n b +ns) (n b + n s ) n obs /n obs != the Poisson probability to observe n obs events given an expectation of (n b + n s ), • P b,s (E i , θ i ) = probability density (normalized to one) for a photon i resulting from a background (b) or signal (s) process to have an energy E i and polar angle θ i , and • f b,s = n b,s /(n b + n s ) = relative fraction of background or signal events.
The expected number of background events, n b , is determined from a Monte Carlo sample of events generated by KORALZ.The signal efficiency, ǫ s , is calculated by integrating an energy and angular distribution function over the region of kinematic acceptance (including kinematic-consistency and degraded-resolution cuts), and scaling the result to account for additional efficiency losses due to other cuts in the single-photon selection.This energy and angular distribution function serves as the the signal probability density function P s (E, θ).The background probability density function P b (E, θ) is obtained from a parameterization of the distributions formed by simulated ννγ(γ) events generated by KORALZ and NUNUGPV98.In the region | cos θ| < 0.72, the photon energy and angular distributions are independent, with the angular distribution given by 1/sin 2 θ and the energy distribution determined with a parameterization.In the region | cos θ| > 0.72, energy resolution is dependent on cos θ.In this region, energy distributions are parameterized separately for slices of width 0.1 in | cos θ|; when properly normalized, these parameterizations together form the 2-dimensional probability distribution.
The value for σ s is not fixed; it is instead treated as a free parameter in the likelihood function.When properly normalized, the likelihood function can be thought of as a probability density function for a hypothesized σ s consistent with the observed data.Restricting σ s to be non-negative, the 95% CL upper limit σ 95 is therefore determined from the following equation: The validity of the likelihood-based method was tested with Monte Carlo simulations.Two types of test were performed.The first test relies on the following definition of a 95% CL limit: if the true signal cross-section happened to be equal to the 95% CL limit, then at least 95% of a large number of identical experiments would result in data that are more signal-like (as defined by a likelihood comparison) than the experiment from which the 95% CL likelihood was derived.This type of test was performed for various combinations of large and small numbers of expected signal and background events, using events generated randomly according to distributions derived from the expected XY signal and ννγ(γ) background.The second type of test involved the calculation of a large number of 95% CL upper limits from samples generated with an identical signal cross-section.In this case, at least 95% of the calculated upper limits should be larger than the actual signal cross-section.Both types of test yielded agreement with expectations to within statistical errors (typically less than 0.5%).In addition, tests with fully simulated Monte Carlo events were used to test for any biases due to the signal and background parameterizations.No significant biases were found.
Sources of systematic errors are the same as for the event-counting method, and are treated in the same manner.However, when calculating limits using the likelihood method, an additional 4% relative uncertainty in experimental sensitivity has been added to account for effects due to uncertainties in the background and signal parameterizations.The 4% estimate is the result of tests using different parameterizations; the uncertainty is again treated according to the method given in reference [26].
The 95% CL upper limits resulting from the likelihood method are shown as a function of M X and M Y in Figure 10.The values range from 23 fb to 371 fb.

5.1.3
Search for e + e − → XY, X → Yγ ; Special case: M Y ≈ 0 The case M Y ≈ 0 is applicable to excited neutrino models and to some supersymmetric models mentioned earlier.The results presented above include this case and no separate analysis is performed, although the results are highlighted here.The upper limits on σ(e + e − → XY) • BR(X → Yγ) for M Y ≈ 0 as a function of M X range from 23 to 107 fb using the likelihood method, and from 43 to 170 fb using the event-counting method.The limits are shown in Figure 11 for both methods.Also shown is the average limit expected in the absence of signal for the event-counting method.The actual limit is significantly lower than the expected limit (see discussion in Section 5.1.2).Table 5 gives more details of the search results.

Gravitino Pair Production
A supersymmetric model has been proposed in which the gravitino is very light [8].This model predicts a new source of single-photon events from the process e + e − → G Gγ.We use the singlephoton topology to place constraints on the cross-section and therefore the gravitino mass in this model.The differential cross-section is given by where α is the fine structure constant, G N is the gravitational constant, m G is the gravitino mass, x γ is the photon scaled energy (E γ /E beam ) and θ is the polar angle, with leading to a soft photon energy spectrum.
Two methods were used to determine the efficiency for observing events from e + e − → G Gγ.The first is by reweighting simulated ννγ(γ) events generated by KORALZ.Each simulated event is weighted by f G Gγ (x, cos θ)/f KZ (x, cos θ), where f KZ is the 2-dimensional energy and angular distribution for the most energetic photon in events generated by KORALZ.We use a parameterization f KZ which is valid up to photon energies of about 60 GeV, and thus we establish a kinematic acceptance region for G Gγ production: x T > 0.05, 15 • < θ < 165 • , and E γ < 60 GeV.We note that the cross-section for G Gγ production with photon energies above 60 GeV is negligible.
In addition, a Monte Carlo generator was written to generate events according to the distribution given in Equation 8. Initial state radiation was treated in a manner identical to that used by the EXOTIC [28] generator.Efficiencies were calculated and found to be identical to those using the reweighting method to within statistical errors.
To achieve maximum sensitivity 7 , we place a maximum energy requirement on the observed photon energy: E γ < 30 GeV.With this requirement, we observe 195 candidates with an expected ν νγ background of 179.6 ± 5.4.The efficiency within the above kinematic acceptance region is estimated by the Monte Carlo to be 82.8%, and was assigned a relative 5% uncertainty.Using the event-counting method described earlier, we place a 95% CL cross-section upper limit of 293 fb, giving a lower limit8 on the gravitino mass of 8.7 µeV.

Graviton-Photon: Search for Extra Space Dimensions
There has been recent interest in string theory models which postulate the existence of additional compactified space dimensions; these models allow gravitons to propagate freely in the higher-dimensional space while restricting Standard Model particles to a 3+1 dimensional hypersurface [7].The fundamental mass scale in this class of theories, M D , governs the rate of graviton production; it is possible that direct graviton-photon production could occur at significant rates at LEP2 energies with an experimental signature of a single photon with missing energy [5,6].The single-photon search topology can therefore be used to place constraints on the fundamental mass scale M D (or, equivalently, on the size of the extra dimensions) by placing limits on the graviton-photon cross-section.
The differential cross-section is given by where δ is the number of extra dimensions and S δ−1 is the surface area of a δ-dimensional sphere of unit radius, with We use the event-counting method to place limits on graviton-photon production in the cases 2 ≤ δ ≤ 7. We use the same kinematic acceptance region as described in Section 5.1.4:x T > 0.05, E γ < 60 GeV, and 15 • < θ < 165 • .The expected photon energy spectrum is soft, so in order to improve sensitivity, we require observed photon energies to be less than 34 GeV for all values of δ.
The efficiency for observing events from e + e − → Gγ is determined by reweighting simulated ννγ(γ) events generated by KORALZ in a manner equivalent to that described in Section 5.1.4:simulated events are weighted by f Gγ (x, cos θ)/f KZ (x, cos θ), where f KZ is the same KORALZ parameterization used earlier.The distributions f KZ and f Gγ have roughly similar shapes, and the calculated efficiency is insensitive to small differences in the KORALZ parameterization.
We observe 208 candidates with an expected ν νγ contribution of 196.0 ± 5.9 events.Assuming a 5% relative uncertainty in the signal efficiency, upper limits on the cross-section and corresponding lower limits on M D for 2 ≤ δ ≤ 7 have been calculated 9 .The results are given in Table 6 based on the convention 10 of Equation 2 in reference [5].The searches for e + e − → XX, X → Yγ, both for the general case discussed here and the special case of M Y ≈ 0 discussed in 5.2.2, use the methods described in our previous publication [1].Selected events are classified as consistent with a given value of M X and M Y if the energy of each of the photons falls within the region kinematically accessible to photons from the process e + e − → XX, X → Yγ, including resolution effects.The selection efficiencies at each generated grid point for the e + e − → XX, X → Yγ Monte Carlo events are shown in Table 7.These values include the efficiency of the kinematic consistency requirement which is higher than 95% at each generated point in the (M X , M Y ) plane for which M X − M Y > 5 GeV.Events from e + e − → ννγγ(γ) are typically characterized by a high-energy photon from the radiative return to the Z 0 and a second lower energy photon.The kinematic consistency requirement is such that the two photons must have energies within the same (kinematically accessible) region.As M X and M Y increase, the allowed range of energy for the photons narrows, and fewer ννγγ(γ) events will be accepted.For the 24 selected events, the distribution of the number of events consistent with a given mass point (M X ,M Y ) is consistent with the expectation from e + e − → ννγγ(γ) Monte Carlo, over the full (M X ,M Y ) plane.

Acoplanar-Photons
In our previous publication, because of concerns about the modelling of the Standard Model process e + e − → ννγγ(γ), all limits derived from the acoplanar photons analysis were obtained without accounting for the expected background.The theoretical situation is now greatly improved, with two event generators [12,13] agreeing to better than 1% for the total cross-section for this process within the kinematic acceptance of this analysis.For that reason, in this paper, all limits derived from this selection have been calculated taking the (KORALZ) background estimate into account.Figure 12 shows the 95% CL exclusion regions for σ(e + e − → XX) • BR 2 (X → Yγ).The limits vary from 33 fb to 103 fb for M X > 45 GeV and M X − M Y > 5 GeV.In the region 2.5 GeV ≤ M X − M Y < 5.0 GeV, the efficiency falls off rapidly (see Table 7).As this rapid fall increases the associated uncertainty in the efficiency, no limits have been set in this region.
Systematic errors are due primarily to limited Monte Carlo statistics at the generated (M X , M Y ) points and the uncertainty on the efficiency parameterization across the (M X , M Y ) plane.The combined relative uncertainty on the efficiency varies from about (3-6)% across the plane (for M X −M Y > 5 GeV).All systematic uncertainties are accounted for in the manner advocated in reference [26].This also applies to the limits for the M Y ≈ 0 case, presented in the next section. 9Evaluated also with α = 1 128 10 The convention of reference [6] differs.

5.2.2
Search for e + e − → XX, X → Yγ ; Special case: M Y ≈ 0 For the special case of M Y ≈ 0 the kinematic consistency requirements differ from those used for the general case.One can calculate [29] the maximum mass, M max X , which is consistent with the measured three-momenta of the two photons, assuming a massless Y.A cut on M max X provides further suppression of the ννγγ(γ) background while retaining high efficiency for the signal hypothesis.This is discussed in more detail in reference [2].We require that the maximum kinematically allowed mass be greater than M X − 5 GeV, which retains (95.5 +2.0 −1.0 )% relative efficiency for signal at all values of M X while suppressing much of the remaining ννγγ(γ) background.In our previous publication we also applied a recoil-mass cut at 80 GeV.In the case where background is not accounted for in the limit calculation, such a cut improves the expected sensitivity of the analysis.With background subtraction this is no longer the case.Therefore that cut has been removed.
Figure 13 shows the expected M max X distribution for signal Monte Carlo events with M X = 90 GeV and for e + e − → ννγγ(γ) Monte Carlo events.Also shown is the distribution of the selected data events.For the M Y ≈ 0 case, the efficiencies calculated from Monte Carlo events are shown in Table 8 after application of the event selection criteria and then after the cut on M max X .Also shown in Table 8 are the number of selected events consistent with each value of M X as well as the expected number of events from e + e − → ννγγ(γ).The number of selected events consistent with a given value of M X varies from 14, for M X ≥ 45 GeV, to 3 events at the kinematic limit.The expected number of events decreases from 15.8 ± 0.7 at M X ≥ 45 GeV to 1.34 ± 0.07 consistent with M X ≥ 94 GeV.
Based on the efficiencies and the number of selected events, we calculate a 95% CL upper limit on σ(e + e − → XX) • BR 2 (X → Yγ) for M Y ≈ 0 as a function of M X .This is shown as the solid line in Figure 14.The limit is between 50 and 80 fb for M X values from 45 GeV up to the kinematic limit.Also shown as a dashed line is the expected limit, defined as the average limit one would expect to set in the absence of signal.The limits can be used to set model-dependent limits on the mass of the lightest neutralino in supersymmetric models in which the NLSP is the lightest neutralino and the LSP is a light gravitino (X = χ0 1 , Y = G).Shown in Figure 14 as a dotted line is the (Bornlevel) cross-section prediction from a specific light gravitino LSP model [30] in which the neutralino composition is purely bino, with m ẽR = 1.35m χ0 1 and m ẽL = 2.7m χ0 1 .Within the framework of this model, χ0 1 masses between 45 and 88.3 GeV are excluded at 95% CL.
As described in section 2, the efficiencies over the full angular range have been calculated using isotropic angular distributions for production and decay of X.The validity of this model has been examined based on the angular distributions calculated for photino pair production in reference [31].For models proposed in reference [32], the production angular distributions are more central and so this procedure is conservative.For a 1 + cos 2 θ production angular distribution expected for t-channel exchange of a very heavy particle according to reference [31], the relative efficiency reduction would be less than 2% at all points in the (M X , M Y ) plane.

5.2.3
Search for e + e − → XX, X → Yγ ; Special case: M Y ≈ 0 and macroscopic decay length As an extension to the special case of M Y ≈ 0, we consider the sensitivity of the acoplanar-photons search to e + e − → XX, X → Yγ when X has a macroscopic decay length.This extension evaluates the selection efficiency using signal Monte Carlo samples with various lifetime values τ X .The X → Yγ decay is treated by a modified OPAL detector simulation package designed to handle delayed decays (as described in section 2).The range of lifetimes considered extends from the near-zero lifetime value of τ X = 10 −15 s to τ X = 10 −7 s at which point cτ is 30 m.Such long lifetimes lead to decays outside the sensitive volume of the OPAL detector (determined by the outer radius of the ECAL), which implies a natural cutoff in sensitivity.
As with the zero lifetime case, the (M max X > M X − 5 GeV) requirement is applied.The relative efficiency of this cut decreases with increasing lifetime, but exceeds 90% for all (M X , τ X ) combinations with high (greater than 50%) selection efficiencies.For (M X , τ X ) combinations with lower selection efficiencies, the relative efficiency falls more rapidly with increasing τ X , but does not drop below 60% for the (M X , τ X ) combinations considered.The loss in relative efficiency of the M max X cut is expected due to the definition of M max X , but the cut is maintained to allow comparison between the prompt decay and macroscopic decay length selection efficiencies.The resulting selection efficiencies for macroscopic decay lengths are given in Table 9.Also listed in Table 9 are the number of events observed in the data and the expected e + e − → ννγγ(γ) background.The efficiencies indicate that sensitivity to macroscopic decay lengths is maintained up to lifetimes of 10 −9 s but falls rapidly for longer lifetimes.This is as expected, since for τ X = 10 −9 s the decay length ranges from 30 to 3 cm for the values of M X considered, so that the X → Yγ decays typically occur well within the confines of the detector.For lifetimes greater than 10 −9 s, the decrease in efficiency is two-fold.The primary loss is due to the increased decay length, implying fewer X → Yγ decays within the sensitive detector volume.Selection efficiency is also lost with increasing τ X since larger lifetimes result in a delayed arrival time of the photon at the ECAL, which, if sufficiently late, causes the event to be vetoed due to the timing cuts imposed by the analysis.The effect of the event kinematics is also seen in the increase in selection efficiency for a given lifetime as M X approaches its threshold value, corresponding to a drop in β from 0.84 at M X =50 GeV to 0.10 at M X =94 GeV.
Contributions to the systematic error in the selection efficiency are (in order of significance) the limited Monte Carlo statistics, the efficiency parameterisation in the (M X , τ X ) plane and the modelling of the timing cuts.The total contributions are given in Table 9.
Given the efficiencies and the number of selected events, the 95% CL upper limit on σ(e + e − → XX)•BR 2 (X → Yγ) for M Y ≈ 0 as a function of M X and τ X is calculated, and is shown in figure 15.The background-subtracted limits from the prompt decay analysis are used to define the profile as a function of M X for τ X = 10 −15 s.The interpolation of the limit to larger lifetimes is done by factoring the selection efficiencies of Table 9 into the 95% confidence level limits on σ(e + e − → XX) • BR 2 (X → Yγ) for M Y ≈ 0. As in Section 5.2.2 the Born-level cross-section of the neutralino NLSP Gravitino LSP model of [30] is used to give an exclusion region in the (neutralino mass, lifetime) plane.This exclusion region is superimposed on Figure 15 and gives the excluded domain in the (M X ,τ X ) plane within the same specific light gravitino LSP model [30] discussed above.

Conclusions
We have searched for photonic events with large missing energy in two topologies in data taken with the OPAL detector at LEP, at a centre-of-mass energy of 189 GeV.
In the single-photon selection, which requires at least one photon with x T > 0.05 in the region 15 • < θ < 165 • (| cos θ | < 0.966), 643 events are observed in the data.The background-subtracted crosssection measurement of 4.35 ± 0.17 (stat) ± 0.09 (sys) pb is consistent with the KORALZ prediction of 4.66 pb from the Standard Model ννγ(γ) process.Interpreting the results as a measurement of the effective number of light neutrino species, we measure N ν = 2.69 ± 0.13(stat) ± 0.11(sys).We also observe significant W contributions to the cross-section with a rate consistent with the Standard Model expectation.
We calculate upper limits on the cross-section times branching ratio for the process e + e − → XY, X → Yγ using two methods: an event-counting method, which is insensitive to energy or angular distribution shapes and is therefore relatively model-independent, and a likelihood method, which assumes isotropic production and decay angular distributions and has greater sensitivity.In the region of interest in the (M X , M Y ) plane, the limits vary from 43 to 409 fb using the event-counting method and from 23 to 371 fb using the likelihood method.These limits include the special case of M Y ≈ 0, where the limit varies between 43 and 170 fb using the event-counting method and from 23 to 107 fb using the likelihood method.We note that some of these limits are much more stringent than would be expected on average in the absence of signal contributions.
We set a 95% CL cross-section upper limit on e + e − → G Gγ production of 293 fb implying a lower limit on the gravitino mass of 8.7 µeV in the superlight gravitino model of reference [8].In the context of string theory models with extra space dimensions, upper limits on the cross-section for graviton-photon production of between 309 and 271 fb at 95% CL are set, giving lower limits on the fundamental mass scale varying between 1086 and 470 GeV for between 2 and 7 additional dimensions respectively.
The acoplanar-photons selection requires at least two photons with scaled energy x γ > 0.05 within the polar angle region 15 • < θ < 165 • or at least two photons with energy E γ > 1.75 GeV with one satisfying | cos θ | < 0.8 and the other satisfying 15 • < θ < 165 • .In each case, the requirement p γγ T /E beam > 0.05 is also applied.There are 24 events selected.The KORALZ prediction for the contribution from e + e − → ννγγ(γ) is 26.9 ± 1.2 events; contribution from other sources is 0.11 events.The number of events observed in the data and their kinematic distributions are consistent with Standard Model expectations.We derive 95% CL upper limits on σ(e + e − → XX) • BR 2 (X → Yγ) ranging from 33 to 103 fb for the general case of massive X and Y.For the special case of M Y ≈ 0, the 95% CL upper limits on σ(e + e − → XX) • BR 2 (X → Yγ) range from 50 to 80 fb.For the case of macroscopic decay lengths, these values range from 50 fb at short lifetimes to 8.8 pb for lifetimes as long as 10 −7 s.
The results of the acoplanar-photons search are used to place model-dependent lower limits on the χ0 1 mass in a specific light gravitino LSP model [30].Masses between 45 and 88.3 GeV are excluded at 95% CL for promptly decaying neutralinos, while for a χ0 1 lifetime of 10 −8 s (cτ = 3m), masses between 45 and 81 GeV are excluded.
• Improved redundancy in the rejection of beam-related backgrounds for conversion candidates with x T < 0.1 and | cos θ| > 0.90 is introduced, based on the measured z-difference of the track's point of closest approach to the interaction point.
• An additional beam-halo veto was implemented; events are rejected if ECAL clusters exceeding 100 MeV in deposited energy are found within 20 cm in the r − φ plane but on opposite ends of the ECAL endcap.This veto provides better rejection of beam-halo background and has negligible effect on signal efficiencies.
• A minor addition was made to the special background vetoes used for non-conversion candidates.
Events are now explicitly rejected if there is a large amount of hadronic energy in an HCAL cluster associated with the ECAL cluster of the primary photon candidate.This provides better redundancy with existing cuts in terms of background rejection, with negligible inefficiency even for beam energy photons.

A.2 Acoplanar-Photons Selection
The acoplanar-photons selection is almost identical to that described in our previous publication.
The only difference involves the timing requirements for candidates with associated TOF hits.These were relaxed for certain periods of the data-taking in the same manner as described above for the single-photon selection.1: Summary of results for the single and acoplanar photons selections.Shown are the number of events observed in the data; the number of expected e + e − → ννγ(γ) (single-photon) or e + e − → ννγγ(γ) (acoplanar-photons) events based on the KORALZ event generator; the number of events expected from backgrounds; the selection efficiency within the kinematic acceptance of the selections (defined in section 3); the background-subtracted, measured cross-section within the kinematic acceptance; and the expected cross-section based on the KORALZ generator.Errors when not marked are statistical and systematic added in quadrature.When two errors are shown, the first is statistical and the second is systematic.6: Results for graviton-photon production in the context of extra dimensions.Shown are the number of extra dimensions δ, the efficiency ǫ, the average upper limit on the signal cross-section σ 95 exp expected in the absence of signal, the measured 95% CL upper limit on signal cross-section σ 95 , and the corresponding lower limit on the mass scale M D .Efficiencies and cross-section limits are evaluated within the restricted kinematic acceptance described in Section 5.1.5.      Figure 8: Consistency between the number of observed events and the number of events expected from e + e − → ννγ(γ) as estimated by KORALZ.Plotted is the probability for the number of expected events to fluctuate to the number of observed events or more, as described in the text.Values greater than 50% indicate a downward fluctuation where the number of observed events is fewer than expected.The boundaries and delineated regions are as defined for Figure 7. of M X , assuming M Y ≈ 0, calculated using the likelihood-based method and the event-counting method as described in the text.Also shown is the average limit expected in the absence of signal, calculated using the event-counting method.

Figure 1 :Figure 2 :
Figure 1: The measured value of σ(e + e − → γ(γ) + invisible particle(s)), within the kinematic acceptance of the single-photon selection, as a function of √ s.The data points with error bars are OPAL measurements at √ s = 130, 136, 161, 172, 183 and 189 GeV.The curve is the prediction for the Standard Model process e + e − → ννγ(γ) from the KORALZ generator.

Figure 3 :Figure 4 :Figure 5 :
Figure 3: Plots of kinematic quantities for the selected acoplanar-photons events for √ s = 189 GeV.a) Recoil-mass distribution.b) Distribution of the scaled energy of the second photon (x 2 ).c) Distribution of the invariant mass of the γγ system.d) Scaled transverse momentum distribution for the γγ system.The data points with error bars represent the selected OPAL data events.In each case the histogram shows the expected contribution from e + e − → ννγγ(γ) events, from KORALZ, normalized to the integrated luminosity of the data.

Figure 6 :
Figure 6: The observed single-photon energy distribution after application of the degraded-resolution cuts.The data are shown as the points with error bars and are compared with the expectation from the Standard Model process e + e − → ννγ(γ), evaluated with KORALZ.The extra hatched contribution around 45GeV illustrates the expected additional contribution from XY production which is excluded at 95% CL using the event-counting method for particular X and Y masses.The hatched contribution corresponds to M X =125 GeV and M Y =63 GeV with a cross-section times branching ratio of 142 fb.

Figure 7 :
Figure 7: Contour plots for the √ s = 189 GeV data sample search for e + e − → XY, X → Yγ, after application of both degraded-resolution and kinematic-consistency cuts.For each set of mass values (M X , M Y ), a) shows the range of numbers of observed single-photon candidate events between each contour and b) shows the range of numbers of expected events from e + e − → ννγ(γ) between each contour.The ννγ(γ) expectation was derived using KORALZ.Lines are drawn around the boundaries defined by M X + M Y = 189 GeV, M X = M Y , and M X + M Y = M Z .The region in which M X − M Y < 5 GeV is not considered.

Table 2 :
Number of eventse + e − → e + e − γ 0.3 ± 0.2 e + e − → ℓ + ℓ − ν ν(γ)1.80 ± 0.14 e + e − → e + e − ℓ + ℓ − Numbers of events expected from various background processes contributing to the singlephoton event sample.The Standard Model background contributions are given by process.Also shown are the expected other backgrounds from cosmic rays and beam related sources.The errors shown are statistical.

Table 3 :
Summary of the experimental systematic errors related to the efficiency and normalisation of the single photon cross-section measurement.

Table 4 :
[1,2]easured cross-section for the process e + e − → ννγγ(γ), within the kinematic acceptance defined in section 3, for different centre-of-mass energies.For √ s = 130 and 136 GeV the measurements are the weighted average of the results obtained from the 1997 data and the results obtained from the 1995 data.Results for √ s ≤ 183 GeV are taken from our previous publications[1,2].The final column shows the cross-section predictions from KORALZ.The quoted errors are statistical.M X E min E max N obs N expected

Table 5 :
Results for e + e − → XY, X → Yγ, with M Y ≈ 0. Shown are values for the mass M X , the minimum and maximum energies (E min , E max ) allowed by the kinematic consistency cuts, the number of observed events (N obs ), the number of events expected from e + e − → ννγ(γ) (N expected ), the value of P fluct as defined in Equation4, the calculated 95% CL upper limit on the number of new-physics signal events (N 95 ), and the average value of N 95 expected in the absence of signal ( N 95 ).Upper limits are calculated using the event-counting method.Masses and energies are in GeV.

Table 7 :
Acoplanar-photons selection efficiencies (%) for the process e + e − → XX, X → Yγ at √ s = 189 GeV for various M X and M Y (in GeV), after application of kinematic-consistency cuts.These efficiency values are used to perform the efficiency parameterization across the (M X , M Y ) plane.The errors shown are due to Monte Carlo statistics only.Efficiencies for the generated points at M Y = 20 and M Y = M X − 15 are not shown, but are similar to those for M Y = 0 and M Y = M X /2.

Table 8 :
X .The last column shows the expected number of events from the process e + e − → ννγγ(γ) (KORALZ).The errors shown due to from Monte Carlo statistics only.
Acoplanar-photons event selection efficiencies (%), as a function of mass, for the process e + e − → XX, X → Yγ, for M Y ≈ 0 at √ s = 189 GeV.The first column shows the efficiency of the selection described in section 3.2.The second column shows the efficiency (%) after the additional cut on M max X .The third column shows the number of selected events consistent with the mass value M