Inclusive Lambda Production in Two-Photon Collisions at LEP

The reactions e^+e^- ->e^+e^- Lambda X and e^+e^- ->e^+e^- Lambda X are studied using data collected at LEP with the L3 detector at centre-of-mass energies between 189 and 209 GeV. Inclusive differential cross sections are measured as a function of the lambda transverse momentum, p_t, and pseudo-rapidity, eta, in the ranges 0.4 GeV<p_t<2.5 GeV and |\eta|<1.2. The data are compared to Monte Carlo predictions. The differential cross section as a function of p_t is well described by an exponential of the form A exp (- p_t /<p_t>)$.


Introduction
Two-photon collisions are the main source of hadron production in high-energy e + e − interactions at LEP via the process e + e − → e + e − γ * γ * → e + e − hadrons. The outgoing electron and positron carry almost the full beam energy and their transverse momenta are usually so small that they escape, undetected, along the beam pipe. In this process, the negative fourmomentum squared of the photons, Q 2 , has an average value Q 2 ≃ 0.2 GeV 2 and they can therefore be considered as quasi-real. In the Vector Dominance Model (VDM), each photon can transform into a vector meson with the same quantum numbers, thus initiating a strong interaction process with characteristics similar to hadron-hadron interactions. This process dominates in the "soft" interaction region, where hadrons are produced with a low transverse momentum, p t . Hadrons with high p t may also be produced by the QED process γγ → qq (direct process) or by QCD processes originating from the partonic content of the photon (resolved processes).
The L3 Collaboration has previously studied inclusive π 0 , K 0 S [1] and charged hadron [2] production in quasi-real two-photon collisions. In this Letter, the inclusive Λ and Λ production 1) from quasi-real photons is studied for a centre-of-mass energy of the two interacting photons, W γγ , greater than 5 GeV. The results are expressed in bins of transverse momentum, p t , and pseudo-rapidity, η, in the ranges 0.4 GeV < p t < 2.5 GeV and |η| < 1.2. The η range is further divided in two subsamples with 0.4 GeV < p t ≤ 1 GeV and 1 GeV < p t < 2.5 GeV. The data sample corresponds to a total integrated luminosity of 610 pb −1 collected with the L3 detector [3], at e + e − centre-of-mass energies √ s = 189 − 209 GeV, with a luminosity weighted average value of √ s = 198 GeV . Results on inclusive Λ production for a smaller data sample at lower √ s were previously reported by the TOPAZ Collaboration in the range 0.75 GeV < p t < 2.75 GeV [4]. The H1 Collaboration investigated the Λ photoproduction at HERA [5].

Monte Carlo simulation
The process e + e − → e + e − hadrons is modelled with the PYTHIA [6] and PHOJET [7] event generators with two times more luminosity than the data. Both generators include the VDM, direct and resolved processes. In PYTHIA, including these processes for each photon leads to six classes of interactions. A smooth transition between these classes is then obtained by introducing a transverse momentum scale. The two-photon luminosity function is based on the EPA approximation [8] with a cut at the mass of the rho meson. The SaS-1D parametrisation is used for the photon structure [9]. PHOJET relies on the Dual Parton Model [10], with a soft and a hard component. The two-photon luminosity functions is calculated in the formalism of Reference 8. The leading-order GRV parametrisation is used for the photon structure [11]. For both programs, matrix elements are calculated at the leading order and higher-order terms are simulated by parton shower in the leading-log approximation. The fragmentation is performed using the Lund string fragmentation scheme as implemented in JETSET [6], which is also used to simulate the hadronisation process. The strangeness suppression factor in JETSET is chosen as 0.3, while a value α S (m Z ) = 0.12 is used for the strong coupling constrant. From a study of Monte Carlo events, the hard component is found to be larger than the soft component. Their ratio goes from around two at low values of p t to around three at high values of p t .
The following Monte Carlo generators are used to simulate the background processes: KK2f [12] for the annihilation process e + e − → qq (γ); KORALZ [13] for e + e − → τ + τ − (γ); KORALW [14] for e + e − → W + W − and DIAG36 [15] for e + e − → e + e − τ + τ − . The response of the L3 detector is simulated using the GEANT [16] and GHEISHA [17] programs. Time dependent detector inefficiencies, as monitored during each data taking period, are included in the simulations. All simulated events are passed through the same reconstruction program as the data.

Event selection
Two-photon interaction events are mainly collected by the track triggers [18], with a low transverse momentum threshold of about 150 MeV, and the calorimetric energy trigger [19]. The selection of e + e − → e + e − hadrons events [20] is based on information from the central tracking detectors and the electromagnetic and hadronic calorimeters. It consists of: • A multiplicity cut. To select hadronic final states, at least six objects must be detected, where an object can be a track satisfying minimal quality requirements or an isolated calorimetric cluster of energy greater than 100 MeV.
• Energy cuts. The total energy deposited in the calorimeters must be less than 40% of √ s to suppress events from the e + e − → qq(γ) and e + e − → τ + τ − (γ) processes. In addition, the total energy in the electromagnetic calorimeter is required to be greater than 500 MeV in order to suppress beam-gas and beam-wall interactions and less than 50 GeV to remove events from the annihilation process e + e − → qq(γ).
• An anti-tag condition. Events with a cluster in the luminosity monitor, which covers the angular region 31 mrad < θ < 62 mrad, with an electromagnetic shower shape and energy greater than 30 GeV are excluded from the analysis. In addition, events with an electron scattered above 62 mrad are rejected by the energy cuts.
• A mass cut. The invariant mass of all visible particles, W vis , must be greater than 5 GeV. In this calculation, the pion mass is attributed to tracks while isolated electromagnetic clusters are treated as massless. The distribution of W vis for data and Monte Carlo after all other cuts are applied is shown in Figure 1. Values of W vis up to 200 GeV are accessible.
About 3 million hadronic events are selected by these criteria with an overall efficiency of 45%. The background level of this sample is less than 1% and is mainly due to the e + e − → qq(γ) and e + e − → e + e − τ + τ − processes. The background from beam-gas and beam-wall interactions is found to be negligible.
The Λ baryons are reconstructed using the decay Λ → pπ. Secondary decay vertices are selected which are formed by two oppositely charged tracks. The secondary vertices must satisfy the following criteria: • The distance, in the plane transverse to the beam direction, between the secondary vertex and the primary e + e − interaction point must be greater than 3 mm.
• The angle between the total transverse momentum vector of the two outgoing tracks and the direction in the transverse plane between the primary interaction point and the secondary vertex must be less than 100 mrad.
The distributions of these variables are presented in Figure 2. The proton is identified as the track with the largest momentum. Monte Carlo studies show that this association is correct for more than 99% of the vertices. In addition, the dE/dx measurement of both proton and pion candidates must be consistent with this hypothesis with a confidence level greater than 0.01.
After these cuts, about 70000 events are selected. The distribution of the invariant mass of the pπ system, m(pπ), shows a clear Λ peak, as shown in Figures 3 and 4 for the different p t bins listed in Table 1. The resolution of m(pπ) is found to be around 3 MeV and is well reproduced by Monte Carlo simulations.

Differential cross sections
The differential cross sections for Λ production as a function of p t and |η| are measured for W γγ > 5 GeV, with a mean value of 30 GeV, and a photon virtuality Q 2 < 8 GeV 2 with Q 2 ≃ 0.2 GeV 2 . This phase space is defined by cuts at the generator level of the Monte Carlo.
The number of Λ baryons in each p t and |η| bin is evaluated by means of a fit to the m(pπ) spectrum in the interval 1.085 GeV < m(pπ) < 1.145 GeV. The signal is modelled with a Gaussian and the background by a fourth-degree Chebyshev polynomial. All parameters, including the mass and width of the peak, are left free. The results are given in Tables 1, 2  and 3.
The overall efficiencies, also listed in Tables 1, 2 and 3, include reconstruction and trigger efficiencies and take into account the 64% branching fraction of the decay Λ → pπ. The reconstruction efficiency, which includes effects of the acceptance and the selection cuts, is calculated with the PHOJET and PYTHIA Monte Carlo generators. As both generators reproduce well the shapes of the experimental distributions of hadronic two-photon production [20], their average is used to calculate the selection efficiency. The efficiency does not depend on the Q 2 cut-off. The track trigger efficiency is calculated for each data taking period by comparing the number of events accepted by the track and the calorimetric energy triggers. The efficiency of the higher level triggers is measured using prescaled events. The total trigger efficiency increases from 82% for p t < 0.4 GeV to 85% in the high p t region.
The differential cross section as a function of p t is first measured for the three different data samples collected in 1998, 1999 and 2000 at different values of √ s and for different trigger and machine background conditions. Good agreement is obtained between the three measurements. The different data samples are therefore combined into a single data sample.
Sources of systematic uncertainties on the cross section measurements are: background subtraction, scale and resolution uncertainties, Monte Carlo modelling and the accuracy of the trigger efficiency measurement. Their contributions are listed in Table 4. The dominant source of systematic uncertainty is due to background subtraction. It is estimated using different background parameterizations and fit intervals in the fitting procedure. The scale and resolution uncertainties are assessed by varying the selection cuts. The main contributions arises from the secondary vertex selection (3.2%) and the proton and pion identification criteria (2.5%). The uncertainty due to the selection of e + e − → e + e − hadrons events is 1%. The Monte Carlo modelling uncertainty, taken as half the relative difference between PHOJET and PYTHIA, increases with p t from 1% to 5%. A systematic uncertainty of 2% is assigned to the determination of the trigger efficiency, which takes into account the determination procedure and time stability.
The sum of the differential cross sections for the e + e − → e + e − ΛX and e + e − → e + e − ΛX processes as a function of p t for |η| < 1.2 is presented in Table 1 and Figure 5. Mass effects explain the lower value obtained in the first bin. The behaviour of the cross section for 0.75 GeV < p t < 2.5 GeV is well described by an exponential of the form A exp(−p t / p t ), as shown in Figure 5a, with a mean value p t = 368 ± 17 MeV. For comparison, the values p t ≃ 230 MeV and p t ≃ 290 MeV are obtained for inclusive π 0 and K 0 S production, respectively [1].
The data are compared to the PHOJET and PYTHIA Monte Carlo predictions in Figure  5b. The region p t < 0.6 GeV is well described by PYTHIA, whereas PHOJET overestimates the cross section. On the contrary, the region 0.6 GeV ≤ p t ≤ 1 GeV is better reproduced by PHOJET while PYTHIA lies below the data. For p t > 1 GeV both PYTHIA and PHOJET underestimate the data. This level of agreement is not unusual in two-photon physics.
The differential cross sections as a function of |η| for 0.4 GeV < p t ≤ 1 GeV and 1 GeV < p t < 2.5 GeV are given in Tables 2 and 3 and shown in Figure 6. Both Monte Carlo programs describe well the almost uniform η shape, while the size of the discrepancy on the absolute normalization depends on the p t range.  Table 1: The number of Λ baryons estimated by the fit, together with the Λ mass, the overall efficiency and the corresponding differential cross section as a function of p t for |η| < 1.2. The first uncertainty on the cross section is statistical and the second systematic. |η| 13.4 ± 0.2 56 ± 1 ± 3 0.6−0.9 2904 ± 70 13.0 ± 0.2 61 ± 1 ± 4 0.9−1.2 2774 ± 89 11.1 ± 0.1 68 ± 2 ± 8 Table 2: The number of Λ baryons estimated by the fit, together with the overall efficiency and the corresponding differential cross section as a function of pseudorapidity for 0.4 GeV < p t ≤ 1 GeV. The first uncertainty on the cross section is statistical and the second systematic.
|η| Table 3: The number of Λ baryons estimated by the fit, together with the overall efficiency and the corresponding differential cross section as a function of pseudorapidity for 1 GeV < p t < 2.5 GeV. The first uncertainty on the cross section is statistical and the second systematic.  Table 4: Systematic uncertainty on the cross section of the e + e − → e + e − ΛX and e + e − → e + e − ΛX processes due to background subtraction, scales and resolution uncertainties, Monte Carlo modelling and trigger efficiency. The total systematic uncertainty is the quadratic sum of the different contributions.   Figure 4: The invariant mass of the pπ system for a) 1.3 GeV < p t ≤ 1.6 GeV, b) 1.6 GeV < p t ≤ 2 GeV and c) 2 GeV < p t < 2.5 GeV. The signal is modelled with a Gaussian and the background by a fourth order Chebyshev polynomial.