Proton-Antiproton Pair Production in Two-Photon Collisions at LEP

The reaction e^+e^- ->e^+e^- proton antiproton is studied with the L3 detector at LEP. The analysis is based on data collected at e^+e^- center-of-mass energies from 183 GeV to 209 GeV, corresponding to an integrated luminosity of 667 pb^-1. The gamma gamma ->proton antiproton differential cross section is measured in the range of the two-photon center-of-mass energy from 2.1 GeV to 4.5 GeV. The results are compared to the predictions of the three-quark and quark-diquark models.


Introduction
Electron-positron colliders are a suitable place for the study of two-photon interactions via the process e + e − → e + e − γ * γ * → e + e − X, where γ * denotes a virtual photon. The outgoing electron and positron carry almost the full beam energy and are usually undetected, due to their small transverse momenta. The final state X has, therefore, a low mass as compared to the e + e − center-of-mass energy, √ s, and its transverse momentum is almost zero. The small photon virtuality allows the extraction of the cross section σ(γγ → X) for real photon collisions, once the photon flux is calculated by QED [1]. Calculations of the cross section γγ → baryon antibaryon were performed using the hard scattering approach of Brodsky and Lepage [2]. In this formalism, the process is factorized into a perturbative γγ → qq amplitude and a non-perturbative part described by the quark distribution functions of the baryon. Such calculations with three-quark distribution functions [3,4] yielded results about one order of magnitude below early γγ → pp measurements [5] for γγ center-of-mass energies, W γγ , less than 3 GeV. The quark-diquark model [6] was proposed as a possible way to model non-perturbative effects. Within this model, the partonic structure of the baryon is described by a quark-diquark system rather than three quarks. The composite nature of diquarks is taken into account by form factors ensuring that the three-quark model is recovered at an asymptotically large momentum transfer.
This Letter presents a study of the differential cross section of the reaction γγ → pp in the W γγ range from 2.1 GeV to 4.5 GeV. The data sample corresponds to a total integrated luminosity of 667 pb −1 collected with the L3 detector [7] at √ s = 183 − 209 GeV. The analysis is based on the central tracking system and the high resolution BGO electromagnetic calorimeter. The events are selected by the track triggers [8].
The pp final state in two-photon collisions was previously studied at lower √ s [5,9] and recently at LEP [10] with lower statistics. The present study extends the W γγ range. Our results are compared to these experiments and to recent theoretical predictions of the quark-diquark model [11].

Monte Carlo simulation
Monte Carlo events for the γγ → pp reaction, as well as for the background processes γγ → K + K − and γγ → π + π − , are generated with the PC Monte Carlo program [12], for each beam energy, within the formalism of Reference 1. The e + e − → e + e − τ + τ − , e + e − → e + e − µ + µ − and e + e − → e + e − e + e − background processes are simulated with the DIAG36 [13] generator which includes the full set of O(α 4 ) QED diagrams. The generated events are passed through the full L3 detector simulation using the GEANT [14] and GEISHA [15] programs and are reconstructed with the same programs as the data. Time dependent detector inefficiencies, as monitored during the data taking period, are taken into account.

Proton-antiproton event selection
In order to match the experimental acceptance with the range of the theoretical predictions, the data are analysed only in the region | cos θ * | < 0.6, where θ * is the center-of-mass production angle of the proton. Events are first selected by requiring two well-reconstructed tracks of opposite charge. The track selection criteria are: • a distance of closest approach to the interaction point less than 3 mm in the plane transverse to the beam direction, • at least 30 hits, out of a maximum of 62, in the tracking chamber, • a matched energy cluster in the electromagnetic calorimeter, • a transverse momentum, p t , greater than 400 MeV, to ensure a high trigger efficiency and electron rejection.
The identification of γγ → pp events is mainly based on three artificial neural networks, used to separate antiprotons from e − , µ − and h − , where h − represents either a π − or K − . Each neural network consists of five input nodes, a single layer of five hidden neurons and one output neuron. The following measured quantities are associated with the five input nodes: • the momentum of the antiproton, • the probability for the proton mass hypothesis based on the mean energy loss dE/dx measured in the tracker, • the ratio E t /p t , where E t is the transverse energy measured in the electromagnetic calorimeter, • the number of BGO crystals in the calorimetric cluster associated with the antiproton, • the ratio between the energy deposited in the central crystal of the cluster and the sum of the energies deposited in the three-by-three matrix of crystals around it.
The two last variables exploit the typical signature of an antiproton annihilating in the BGO electromagnetic calorimeter: a broad shower spanning several crystals, as opposed to narrow showers for electrons, minimum ionisation deposition for muons and compact showers for low energy hadrons. Each neural network is trained with the corresponding sample from Monte Carlo simulations so that its output value is close to one for antiprotons.
A particle is identified as an antiproton after a cut on the output of the three neural networks. The efficiency for antiproton detection is found to be 74%. Even though the e − neural network rejects more than 98% of the electrons, the e + e − → e + e − e + e − cross section is more than two orders of magnitude greater than the e + e − → e + e − pp cross section and the electron contamination still remains important. To reduce this background, the ratio E t /p t of the proton candidate is required to be less than 0.6, as shown in Figure 1. In addition, the confidence level that the dE/dx measurement is consistent with a proton must be more than 5%. These cuts eliminate more than 95% of the remaining e + e − → e + e − e + e − background.

Exclusive γγ → pp event selection
To select exclusive γγ → pp events the following further cuts are applied: • Events with a photon candidate outside a cone of 36 degrees half-opening angle with respect to the antiproton direction are rejected. A photon candidate is defined as a shower in the electromagnetic calorimeter with at least two adjacent crystals, an energy greater than 100 MeV and no charged tracks within a cone of 200 mrad half-opening angle.
• The square of the total transverse momentum of the proton-antiproton pair ( p t ) 2 must be less than 0.1 GeV 2 . This cut reduces inclusive background events of the type e + e − → e + e − ppX, where X represents one or more unobserved particles.
These cuts yield a total number of 989 selected events. Figure 2 shows the distribution of the effective mass of the pp pair, identified with W γγ . It is obtained by assigning the proton mass to the two tracks.
The exclusive background is estimated by processing the corresponding Monte Carlo events with the same analysis cuts. The γγ → π + π − and γγ → K + K − cross sections given in Reference 16 are used. This background is found to be negligible in the region W γγ < 2.5 GeV and increases up to 25% above 2.5 GeV. The e + e − → e + e − e + e − background varies from 5% in the low mass region to 30% in the high mass region. The e + e − → e + e − µ + µ − contamination is less than 0.2% and neglected. The contamination of inclusive channels and of e + e − → e + e − τ + τ − is estimated by fitting the tail of the ( p t ) 2 distribution. The extrapolation to the region ( p t ) 2 < 0.1 GeV 2 gives an average background level of (2 ± 1)%. The sum of all backgrounds is subtracted bin by bin for the cross section determination. The background composition of each bin is detailed in Table 1.

Cross section measurements
The detection efficiency is determined by Monte Carlo simulation in bins of W γγ and | cos θ * |. It takes into account the track acceptance and selection criteria, the exclusive proton-antiproton identification criteria and the track trigger efficiency. This trigger efficiency is 91% at W γγ = 2.1 GeV and rises to 98% for W γγ > 3 GeV. The higher level trigger efficiencies are estimated from the data themselves, using prescaled events, and range from 88% at W γγ = 2.1 GeV up to 99% for W γγ > 3 GeV. The total selection efficiency is found to be maximum, about 6%, at W γγ ≃ 2.5 GeV and | cos θ * | < 0.1. It decreases both at larger | cos θ * | and W γγ , due respectively to limited angular acceptance and antiproton identification efficiency. The average efficiencies as a function of W γγ for | cos θ * | < 0.6 and as a function of | cos θ * | for W γγ < 4.5 GeV are presented in Figure 3.
The cross sections are evaluated in bins of W γγ and | cos θ * | for different values of √ s. Due to the limited selection efficiency near threshold (Figure 3), the measurement is restricted to the range W γγ > 2.1 GeV. A total of 938 events are selected after this cut. The differential cross sections are integrated to obtain the production cross sections e + e − → e + e − pp for | cos θ * | < 0.6 and 2.1 GeV < W γγ < 4.5 GeV as a function of √ s. Since the results show no significant √ s dependence they are combined into a single measurement at √ s = 197 GeV: σ(e + e − → e + e − pp) = 26.7 ± 0.9 ± 2.7 pb where the first uncertainty is statistical, the second systematic.
The following systematic uncertainties are considered. The uncertainty due to the selection procedure is evaluated by varying the selection cuts. The uncertainty on the neural network outputs is estimated by varying the input variables according to their resolution. An uncertainty of 50% on the γγ → π + π − and γγ → K + K − cross sections is propagated in the background subtraction. The uncertainty on the e + e − → e + e − e + e − background is estimated to be 30%. The uncertainty of 50% on the level of inclusive background is also included. For low values of W γγ the dominant systematic uncertainty is about 7% and is due to selection cut variation. It remains the most important systematic uncertainty for W γγ < 3.2 GeV, as large as 30%. For W γγ > 3.2 GeV, the background subtraction uncertainty becomes the most important, rising to 55% in the last W γγ bin. The systematic uncertainties shown in Table 2 are the quadratic sum of the different sources.
The differential cross section dσ(γγ → pp)/d| cos θ * | for real photon-photon collisions is extracted as a function of W γγ by dividing out the two-photon luminosity function and extrapolating to Q 2 = 0 with a GVDM form factor [17]. The luminosity functions are evaluated for each √ s and W γγ interval. The measured differential cross sections are integrated to obtain the cross sections σ(γγ → pp). An additional uncertainty of 5%, due to the choice of the photon form factor, is included in the systematics. The results are reported in Table 2. The present measurement is of higher statistical precision and extends towards higher values of W γγ than the previous results [5,9,10]. Agreement is observed within the quoted uncertainties, except for the measurements of Reference 10 which lie below our data in the low mass region.

Discussion of the results
The results are compared to the predictions of the three-quark [3] and the recent quark-diquark model [11] in Figure 4. The three-quark prediction is based on the leading order QCD calculations, using the distribution function of Chernyak and Zhitnisky [18]. The quark-diquark calculation is performed with the standard proton distribution amplitudes [11] and includes first order corrections due to the non-vanishing proton mass. While the shapes of the theoretical curves are quite similar, the normalisations are significantly different. The predictions of the three-quark model are about an order of magnitude below the measurement, whereas the quark-diquark predictions describe the data much better. The apparent change in the logarithmic slope of the cross section observed in data at W γγ around 3 GeV is, however, not reproduced by this model. In order to investigate this further, the differential cross sections are summed separately in three mass intervals: 2.1 GeV < W γγ < 2.5 GeV, 2.5 GeV < W γγ < 3.0 GeV and 3.0 GeV < W γγ < 4.5 GeV. The results are reported in Table 3 and plotted in Figure 5. A distinctive difference between the three distributions can be observed. No prediction is available for the diquark model for W γγ < 2.5 GeV, but it can be seen that the data has a qualitatively different behaviour to the diquark predictions, as it is strongly peaked at large angles. Assuming the presence of a single angular momentum state in the s-channel of γγ → pp production, a fit to a single spherical harmonic: Y 0 0 , Y 0 2 , Y 1 2 or Y 2 2 is performed. This does not give a satisfactory result. An acceptable fit is obtained using all of them, but only the Y 0 2 contribution is significantly different from zero. A satisfactory fit is also obtained using only the Y 0 0 and Y 0 2 harmonics with the fractions 8% and 92% respectively, as shown in Figure 5a. The intermediate region exhibits a rather flat dependence, which partially agrees with the model predictions.
The forward peaking behaviour of the differential cross section in the high mass interval is well reproduced by the quark-diquark model. Only in this region, then, the data can be described by the Brodsky-Lepage hard scattering approach. The presence of two distinct production mechanisms can also be seen in Figure 6a and Table 4 by considering separately the cross sections as a function of W γγ in a large angle region, | cos θ * | < 0.3, and in a region 0.3 < | cos θ * | < 0.6. The shape of the efficiency as a function of W γγ is similar for the two regions. The large angle cross section dominates in the low mass region and shows a steeper fall-off with W γγ than the small angle cross section. A fit of the form σ(W ) ∝ W −n to the small angle region shown in Figure 6b gives a reasonably good description with n = 9.8 ± 0.3. The data are also in agreement with the quark-diquark model prediction in this region. On the other hand, Figure 6c, the large angle cross section cannot be described by a simple power law behaviour for W γγ less than 3 GeV and does not follow the quark-diquark prediction. In this case, the change of shape at 3 GeV is similar to, but more pronounced than that in Figure 4.
In conclusion, an accurate study of the γγ → pp process is performed. The present results are of higher statistical precision and extend towards higher values of W γγ than those of previous experiments. Current models are only moderately successful in describing the observed data. These data can provide useful inputs for predictions of other γγ → baryon antibaryon channels [19].  distribution. The first uncertainty is statistical, the second systematic.
dσ(γγ → pp)/d| cos θ * | (nb) | cos θ * | 2.1 GeV < W γγ < 2.5 GeV 2.5 GeV < W γγ < 3.0 GeV 3.0 GeV < W γγ < 4.5 GeV 0 − 0.       Figure 4: The γγ → pp cross sections as a function of W γγ for | cos θ * | < 0.6 compared to the three-quark model calculation [3] and to the recent quark-diquark model prediction [11], available for W γγ > 2.5 GeV. Statistical and systematic uncertainties are added in quadrature. The differential cross section as a function of | cos θ * | for a) 2.1 < W γγ < 2.5 GeV, b) 2.5 < W γγ < 3.0 GeV and c) 3.0 < W γγ < 4.5 GeV. The data are compared to the predictions of the quark-diquark model. The fit in a) is described in the text. a) Figure 6: a) The γγ → pp cross section as a function of W γγ for the large angle region, | cos θ * | < 0.3 (full circle), and the small angle region, 0.3 < | cos θ * | < 0.6 (open circle). b) The small angle and c) the large angle cross section with the quark-diquark model predictions and the fit described in text.