Measurement of the Cross-Section for the Process GammaGamma ->Ppbar at sqrt(s)_ee=183-189 GeV at LEP

The exclusive production of proton-antiproton pairs in the collisions of two quasi-real photons had been studied using data taken at sqrt(s)_ee=183 GeV and 189 GeV with the OPAL detector at LEP. Results are presented for Ppbar invariant masses, W, in the range 2.15 W<<3.95 GeV. The cross-section measurements are compared with previous data and with recent analytic calculations based on the quark-diquark model.


Introduction
The exclusive production of proton-antiproton (pp) pairs in the collision of two quasi-real photons can be used to test predictions of QCD.At LEP the photons are emitted by the beam electrons 1 and the pp pairs are produced in the process e + e − → e + e − γγ → e + e − pp.
The application of QCD to exclusive photon-photon reactions is based on the work of Brodsky and Lepage [1].According to their formalism the process is factorized into a nonperturbative part, which is the hadronic wave function of the final state, and a perturbative part.Calculations based on this ansatz [2,3] use a specific model of the proton's three-quark wave function by Chernyak and Zhitnitsky [4].This calculation yields cross-sections about one order of magnitude smaller than the existing experimental results [5][6][7][8][9][10], for pp centre-of-mass energies W greater than 2.5 GeV.
To model non-perturbative effects, the introduction of quark-diquark systems has been proposed [11].Within this model, baryons are viewed as a combination of a quark and a diquark rather than a three-quark system.The composite nature of the diquark is taken into account by form factors.
Recent studies [12] have extended the systematic investigation of hard exclusive reactions within the quark-diquark model to photon-photon processes [13][14][15][16].In these studies the crosssections have been calculated down to W values of 2.2 GeV below which the quark-diquark model is no longer expected to be valid.Most of the experimental data, however, have been taken at such low energies.
The calculations of the integrated cross-section for the process γγ → pp in the angular range | cos θ * | < 0.6 (where θ * is the angle between the proton's momentum and the electron beam direction in the pp centre-of-mass system) and for W > 2.5 GeV are in good agreement with experimental results [9,10], whereas the pure quark model predicts much smaller crosssections [2,3].
In this paper, we present a measurement of the cross-section for the exclusive process e + e − → e + e − pp in the range 2.15 < W < 3.95 GeV, using data taken with the OPAL detector at √ s ee = 183 GeV and 189 GeV at LEP.The integrated luminosities for the two energies are 62.8 pb −1 and 186.2 pb −1 .

The OPAL detector
The OPAL detector and trigger system are described in detail elsewhere [17].We briefly describe only those features particularly relevant to this analysis.The tracking system for charged particles is inside a solenoid that provides a uniform axial magnetic field of 0.435 T. The system consists of a silicon micro-vertex detector, a high-resolution vertex drift chamber, a large-volume jet chamber and surrounding z-chambers.The micro-vertex detector surrounds the beam pipe covering the angular range of | cos θ| < 0.9 and provides tracking information in the r-φ and z directions2 .The jet chamber records the momentum and energy loss of charged particles over 98% of the solid angle.In the range of | cos θ| < 0.73, up to 159 points are measured along each track.The energy loss, dE/dx, of a charged particle in the chamber gas is measured from the integrated charges of each hit at both ends of each signal wire with a resolution of about 3.5% for isolated tracks with the maximum of 159 points.The z-chambers are used to improve the track measurement in the z direction.
The barrel time-of-flight (TOF) scintillation counters are located immediately outside the solenoid at a mean radius of 2.36 m, covering the polar angle range | cos θ| < 0.82.The outer parts of the detector, in the barrel and endcaps, consist of lead-glass electromagnetic calorimetry (ECAL) followed by an instrumented iron yoke for hadron calorimetry and four layers of external muon chambers.Forward electromagnetic calorimeters complete the acceptance for electromagnetically interacting particles down to polar angles of about 24 mrad.
The trigger signatures required for this analysis are based on a combination of time-of-flight and track trigger information.

Kinematics
The production of proton-antiproton pairs in photon-photon interactions proceeds via the process where q i , p i denote the four-momenta (i = 1, 2).Each of the two incoming electrons emits a photon and the final state produced by the two colliding photons consists of one proton (p) and one antiproton (p).The four-momentum squared of the two photons is (i = 1, 2) Since the electrons are scattered at small angles and they remain undetected, the four-momenta squared of each of the two photons are small, i.e. the photons are quasi-real.In this case the transverse 3 component of the momentum sum of the proton and the antiproton in the laboratory system is expected to be small whereas the longitudinal component of the momentum sum can be large.The photon-photon centre-of-mass system is generally boosted along the beam axis.The larger the boost, the closer the produced (anti) protons are to the beam direction.This feature, combined with the typically low mass of the final state pp, and the low efficiency for tracking at small angles, leads to significant acceptance losses for these types of events.

Monte Carlo generators
The e + e − → e + e − pp events are simulated with the PC Monte Carlo generator which has been developed to study exclusive photon-photon processes [18].The PC Monte Carlo generator has been expanded for use in this analysis to simulate the kinematics of exclusive baryonantibaryon final states, e + e − → e + e − pp, e + e − → e + e − ΛΛ, and e + e − → e + e − Ξ − Ξ + .A total of 20 000 e + e − → e + e − pp Monte Carlo events have been generated as a control sample for the event selection described in the next section.Similarly, for the trigger and detection efficiency determination, a total of 360 000 e + e − → e + e − pp Monte Carlo events have been generated in the range of 2.15 < W < 3.95 GeV.The background coming from e + e − → e + e − π + π − is generated with the PC Monte Carlo program, as is the feed-down background from protonantiproton pairs coming from the reaction e + e − → e + e − ΛΛ → ppπ + π − where the pions are not detected (200 000 events).
The leptonic photon-photon background processes e + e − → e + e − µ + µ − , e + e − → e + e − e + e − and e + e − → e + e − τ + τ − are simulated with the Vermaseren generator [19].The KORALZ generator [20] is used to simulate the background processes e + e − → µ + µ − , and e + e − → τ + τ − .The e + e − → e + e − background process is simulated with the BHWIDE [21] generator.Table 1 lists all the generated background Monte Carlo samples.Monte Carlo events are generated at √ s ee = 189 GeV only, since the change in acceptance between √ s ee = 183 and 189 GeV is small compared to the statistical uncertainty of the measurement.They have been processed through a full simulation of the OPAL detector [22] and have been analysed using the same reconstruction algorithms that are used for the data.

Event selection
The e + e − → e + e − pp events are selected by the following set of cuts: 1.The sum of the energies measured in the barrel and endcap sections of the electromagnetic calorimeter must be less than half the beam energy.
2. Exactly two oppositely charged tracks are required with each track having at least 20 hits in the central jet chamber to ensure a reliable determination of the specific energy loss dE/dx.The point of closest approach to the interaction point must be less than 1 cm in the rφ plane and less than 50 cm in the z direction.
3. For each track the polar angle must be in the range | cos θ| < 0.75 and the transverse momentum p ⊥ must be larger than 0.4 GeV.These cuts ensure a high trigger efficiency and good particle identification.
4. The invariant mass W of the pp final state must be in the range 2.15 < W < 3.95 GeV.The invariant mass is determined from the measured momenta of the two tracks using the proton mass.
5. The events are boosted into the rest system of the measured pp final state.The scattering angle of the tracks in this system has to satisfy | cos θ * | < 0.6.
6.All events must fulfil the trigger conditions defined in Section 6.
7. The large background from other exclusive processes, mainly the production of e + e − , µ + µ − , and π + π − pairs, is reduced by particle identification using the specific energy loss dE/dx in the jet chamber and the energy in the electromagnetic calorimeter.The dE/dx probabilities of the tracks must be consistent with the p and p hypothesis.
-Events where the ratio E/p for each track lies in the range 0.4 < E/p < 1.84 are regarded as possible e + e − → e + e − e + e − candidates.These events are rejected if the dE/dx probabilities of the two tracks are consistent with the electron hypothesis.
-Events where the ratio E/p for each track is less than 0.8, as expected for a minimum ionizing particle, are regarded as possible background from e + e − → e + e − µ + µ − events.This background is reduced by rejecting events where the dE/dx probability for both tracks is consistent with the muon hypothesis.This cut is also effective in reducing the π + π − background.
-The dE/dx probability for the proton hypothesis has to be greater than 0.1% for each track and it has to be larger than the probabilities for the pion and kaon hypotheses.
The product of the dE/dx probabilities for both tracks to be (anti) protons has to be larger than the product of the dE/dx probabilities for both tracks to be electrons.
8. Cosmic ray background is eliminated by applying a muon veto [23].
9. Exclusive two-particle final states are selected by requiring the transverse component of the momentum sum squared of the two tracks, | p ⊥ | 2 , to be smaller than 0.04 GeV 2 .By restricting the maximum value of Q 2 i , this cut also ensures that the interacting photons are quasi-real.Therefore no further cut rejecting events with scattered electrons in the detector needs to be applied.Fig. 1 shows the | p ⊥ | 2 distribution obtained after applying all cuts except the cut on After all cuts 163 data events are selected, 35 events at √ s ee = 183 GeV and 128 events at √ s ee = 189 GeV.The distribution of measured dE/dx values versus the particle momentum for the selected data events is shown in Fig. 2a.Background from events containing particles other than (anti-)protons is negligible due to the good rejection power of the dE/dx cuts.Since no event remains after applying the event selection to the background Monte Carlo samples, the 95% confidence level upper limits for the number of background events expected to contribute to the selected data sample are given in Table 1.
Since the pp final state is fully reconstructed, the experimental resolution for W (determined with Monte Carlo simulation) is better than 1%.The experimental resolution for | cos θ * | is about 0.014.Fig. 2b shows the W distribution for data and Monte Carlo signal events after the final selection.The Monte Carlo distributions agree well with the data.

Trigger and detection efficiencies
The e + e − → e + e − pp events contain only two tracks with momenta in the range 0.4 GeV to 2 GeV.Special triggers are required to select such low multiplicity events with only low momentum particles and the efficiencies of these triggers must be well known.The e + e − → e + e − pp events are mainly triggered by a combination of triggers using hits in the time-of-flight counters and tracks.The track trigger takes data on the z coordinate of hits from the vertex drift chamber, and from three groups of 12 wires in the jet chamber.The selected events must satisfy at least one of the following trigger conditions: A Two tracks in the barrel region from the track-trigger (TT).This corresponds to an angular acceptance of approximately | cos θ| < 0.75.Condition A is highly efficient for events within its geometrical acceptance but it triggers on two tracks.Condition B is used to trigger on a single track and to measure the trigger efficiency in combination with condition A.
It was checked that the trigger efficiency does not depend on φ.Under this condition and assuming the efficiency of each trigger component for each track to be independent, the combined event trigger can be written as: Here ε TT is the efficiency for one track to be triggered by the track trigger and ε TT−TOF is the efficiency for one track to be triggered by the TT-TOF coincidence.The first term in (3) gives the efficiency for both tracks to be triggered by the track trigger (condition A) and the second term gives the efficiency for the events not triggered by condition A to be triggered by condition B for either of the two tracks.
Data events are used to calculate the trigger efficiency.The trigger efficiency is determined by considering events in which one track detected in one half of the rφ plane (e.g.0 < φ < 180 degrees) satisfies the required track trigger and time-of-flight θ − φ matching hit while the other track in the other half plane (e.g. 180 < φ < 360 degrees) is used to measure the efficiency of the track and time-of-flight triggers.
To study the efficiency with a large event sample, events from photon-photon processes with two tracks in the final state such as e + e − and µ + µ − are used together with the pp events to determine the track trigger efficiency as a function of p ⊥ .This efficiency exceeds 90% for p ⊥ > 0.4 GeV and the observed dependence of the track trigger efficiency on p ⊥ is found to be consistent between electrons and muons.
In a second step, e + e − → e + e − pp Monte Carlo events in the range 2.15 < W < 3.95 GeV are used to obtain the trigger efficiencies as a function of W , ε TRIG (W ).The Monte Carlo events are reweighted here according to the trigger efficiency which has been determined as a function of the transverse momentum p ⊥ using (3).These reweighted Monte Carlo events have been used only to determine the trigger efficiency as function of W . Finally the values obtained for ε TRIG (W ) are normalized to give the overall trigger efficiency ε TRIG = 94.7% above the threshold p ⊥ > 0.4 GeV as already determined from data only.The region W < 2.15 GeV is excluded from the analysis because the trigger efficiency drops rapidly below 70%.
The detection efficiency is determined by comparing the number of Monte Carlo events passing all cuts with the total number of Monte Carlo events generated within a polar angle | cos θ * GEN | < 0.6 in the pp centre-of-mass system: where | cos θ * | refers to the reconstructed polar angle in the pp centre-of-mass system.The detection efficiency is about 2% at high W and about 4% at low W .To be able to compare the measured cross-section with any given model in bins of | cos θ * | and W , the detection efficiency has been determined from the signal Monte Carlo in bins of W and | cos θ * |.

Cross-section measurements
The differential cross-section for the process e + e − → e + e − pp is given by where N ev is the number of events selected in each (W, | cos θ * |) bin, ε TRIG is the trigger efficiency, ε DET is the detection efficiency, L e + e − is the measured integrated luminosity, and ∆W and ∆| cos θ * | are the bin widths in W and in | cos θ * |.
The total cross-section σ(γγ → pp) for a given value of √ s ee is obtained from the differential cross-section dσ(e + e − → e + e − pp)/dW using the luminosity function dL γγ /dW [24]: The luminosity function dL γγ /dW is calculated by the Galuga program [25].The resulting differential cross-sections for the process γγ → pp in bins of W and | cos θ * | are then summed over | cos θ * | to obtain the total cross-section as a function of W for | cos θ * | < 0.6.

Systematic uncertainties
The following sources of systematic uncertainties have been taken into account (Table 5): Luminosity function: The accuracy of the photon-photon luminosity function used in (6) has been estimated by taking into account different models such as the ρ-pole model [26], the Equivalent Photon Approximation (EPA) [27], the Generalized Vector Dominance Model (GVMD) [25], the luminosity functions given in [28], and in [29,30].Each of the resulting numbers has then been compared with the VDM form-factor model [25] which is used for the final result.We take the largest deviation resulting from these comparisons as the uncertainty in the luminosity functions which is about ±5%.
Trigger efficiency: The trigger efficiency model of (3) has been tested by making a comparison between the measured relative frequencies of the four different sub-combinations, (N TT , N TT−TOF ) = (1,1), (2,0), (2,1), (2,2), which can trigger the events, and the predicted values, by using fitted efficiencies for TT and TT-TOF (Section 6).Here N TT denotes the number of tracks triggered by the TT and N TT−TOF the number of tracks triggered by the TT-TOF trigger.Table 2 gives the measured and predicted fractions for events with W > 2.15 GeV.The fit of the two efficiencies ε TT and ε TT−TOF yields a χ 2 of 9.6 over 2 degrees of freedom.Although the fit is poor the overall efficiency result is consistent with the other determination.An additional systematic uncertainty of 5.0% is therefore assigned to the trigger efficiency.This value is obtained by increasing the uncertainties to obtain a normalized χ 2 of one.

Monte Carlo statistics:
The statistical uncertainty on the detection efficiency due to the number of simulated Monte Carlo events varies from 4.5% at low W to 6% at high W .
dE/dx cuts: The systematic uncertainties due to the dE/dx cuts are determined by recalculating the Monte Carlo detection efficiency after varying the dE/dx values: -The measured values are shifted by ±1σ, where σ is the theoretical uncertainty of the measurement.-The measured values are smeared with a Gaussian distribution of width 0.1σ, which is the typical systematic uncertainty of each individual dE/dx measurement, as found for pions from K 0 decays [31].
The modified dE/dx values are transformed into weights and the new detection efficiency is calculated by applying the event selection on the modified Monte Carlo events.The systematic uncertainty assigned to each (W, | cos θ * |) bin is the quadratic sum of the full deviation of the detection efficiency with smearing and the average absolute deviation value of the ±1σ shifted values with respect to the original values.
The systematic uncertainties due to the variation of the dE/dx cuts are larger at high values of W .They vary between 0.1% for W < 2.6 GeV and up to 5% for W > 2.6 GeV.
Residual background: Residual background can come from non-exclusive production of pp pairs in processes like e + e − → e + e − ppππ if the pions are not detected.The | p ⊥ | 2 cut eliminates most of this background.The data events are almost coplanar, i.e. no data event has an acoplanarity of more than 0.262 rad.To estimate the contribution of this background, the shape of the acoplanarity distribution for pp pairs in Monte Carlo e + e − → e + e − ΛΛ and e + e − → e + e − pp events has been fitted to the data.This yields an upper limit of 10 events for this background contribution.An additional systematic uncertainty of 6% is therefore taken into account.
Additional uncertainties due to the measured integrated e + e − luminosity, the track reconstruction efficiency and the momentum resolution for the protons are negligible.The total systematic uncertainty is obtained by adding all systematic uncertainties in quadrature.

Results and Discussion
The measured cross-sections in bins of W are given in Table 3.The average W in each bin has been determined by applying the procedure described in [32].The measured cross-sections σ(γγ → pp) for 2.15 < W < 3.95 GeV and for | cos θ * | < 0.6 are compared with the results obtained by ARGUS [8], CLEO [9] and VENUS [10] in Fig. 3a and to the results obtained by TASSO [5], JADE [6] and TPC/2γ [7] in Fig. 3b.The quark-diquark model predictions [12] are also shown.Reasonable agreement is found between this measurement and the results obtained by other experiments for W > 2.3 GeV.At lower W our measurements agree with the measurements by JADE [6] and ARGUS [8], but lie below the results obtained by CLEO [9], and VENUS [10].The cross-section measurements reported here extend towards higher values of W than previous results.
Fig. 3a-b show the measured γγ → pp cross-section as a function of W together with some predictions based on the quark-diquark model [11,12,15].There is good agreement between our results and the older quark-diquark model predictions [11,15].The most recent calculations [12] lie above the data, but within the estimated theoretical uncertainties the predictions are in agreement with the measurement.
An important consequence of the pure quark hard scattering picture is the power law which follows from the dimensional counting rules [33,34].The dimensional counting rules state that an exclusive cross-section at fixed angle has an energy dependence connected with the number of hadronic constituents participating in the process under investigation.We expect that for asymptotically large W and fixed where n = 8 is the number of elementary fields and The introduction of diquarks modifies the power law by decreasing n to n = 6.This power law is compared to the data in Fig. 3c with σ(γγ → pp) ∼ W −2(n−3) using three values of the exponent n: fixed values n = 8, n = 6, and the fitted value n = 7.5 ± 0.8 obtained by taking into account statistical uncertainties only.More data covering a wider range of W would be required to determine the exponent n more precisely.
The measured differential cross-sections dσ(γγ → pp)/d| cos θ * | in different W ranges and for | cos θ * | < 0.6 are given in Table 4 and in Fig. 4. The differential cross-section in the range 2.15 < W < 2.55 GeV lies below the results reported by VENUS [10] and CLEO [9] (Fig. 4a).Since the CLEO measurements are given for the lower W range 2.0 < W < 2.5 GeV, we rescale their results by a factor 0.635 which is the ratio of the two CLEO total cross-section measurements integrated over the W ranges 2.0 < W < 2.5 GeV and 2.15 < W < 2.55 GeV.This leads to a better agreement between the two measurements but the OPAL results are still consistently lower.The shapes of the | cos θ * | dependence of all measurements are consistent apart from the highest | cos θ * | bin, where the OPAL measurement is significantly lower than the measurements of the other two experiments.
In Fig. 4b-c the differential cross-sections dσ(γγ → pp)/d| cos θ * | in the W ranges 2.35 < W < 2.85 GeV and 2.55 < W < 2.95 GeV are compared to the measurements by TASSO, VENUS and CLEO in similar W ranges.The measurements are consistent within the uncertainties.
The comparison of the differential cross-section as a function of | cos θ * | for 2.55 < W < 2.95 GeV with the calculation of [12] at W = 2.8 GeV for different distribution amplitudes (DA) is shown in Fig. 5a.The pure quark model [2,3] and the quark-diquark model predictions lie below the data, but the shapes of the curves are consistent with those of the data.
In Fig. 5b the differential cross-section dσ(γγ → pp)/d| cos θ * | is shown versus | cos θ * | for 2.15 < W < 2.55 GeV.The cross-section decreases at large | cos θ * |; the shape of the angular distribution is different from that at higher W values.This indicates that for low W the perturbative calculations of [2,3] are not valid.
Another important consequence of the hard scattering picture is the hadron helicity conservation rule.For each exclusive reaction like γγ → pp the sum of the two initial helicities equals the sum of the two final ones [35].According to the simplification used in [11], only scalar diquarks are considered, and the (anti) proton carries the helicity of the single (anti) quark.Neglecting quark masses, quark and antiquark and hence proton and antiproton have to be in opposite helicity states.If the (anti) proton is considered as a point-like particle, simple QED rules determine the angular dependence of the unpolarized γγ → pp differential cross-section [29]: This expression is compared to the data in two W ranges, 2.55 < W < 2.95 GeV (Fig. 5a) and 2.15 < W < 2.55 GeV (Fig. 5b).The normalisation in each case is determined by the best fit to the data.In the higher W range, the prediction ( 8) is in agreement with the data within the experimental uncertainties.In the lower W range this simple model does not describe the data.At low W soft processes such as meson exchange are expected to introduce other partial waves, so that the approximations leading to (8) become invalid [36].

B
A coincidence of at least one barrel track from the track trigger and a θ − φ coincidence of a track from the track trigger with hits from the time-of-flight detector (TOF).The barrel track and the track forming the θ − φ coincidence are not necessarily identical.The angular acceptance of the track trigger is approximately | cos θ| < 0.75, whereas the acceptance of the θ − φ coincidence is | cos θ| < 0.82.

Figure 1 : 6 WFigure 2 :
Figure 1: Distribution of the transverse momentum balance, | p ⊥ | 2 , for the e + e − → e + e − pp events in data (black points with error bars) compared to the Monte Carlo simulation (histogram).The distributions are obtained after applying all cuts except the | p ⊥ | 2 cut.The arrow indicates the cut value.The error bars are statistical only.The Monte Carlo distribution is normalized to the number of selected data events.

Figure 3 :
Figure 3: Cross-sections σ(γγ → pp) as a function of W .The data and the theoretical predictions cover a range of | cos θ * | < 0.6.The data points are plotted at the value of W . a,b)The data are compared to other experimental results[7][8][9][10] and to the quark-diquark model prediction[12].The error bars include statistical and systematic uncertainties, except for TASSO[5] where the uncertainties are statistical only.c) The data are compared to the quark-diquark model of[11] (dash-dotted line), and of[12] (solid line), using the standard distribution amplitude (DA) with and without neglecting the mass m p of the proton, and with the predictions of the power law with fixed and with fitted exponent n.The inner error bars are the statistical uncertainties and the outer error bars are the total uncertainties.

Figure 4 :
Figure 4: Differential cross-sections for γγ → pp as a function of | cos θ * | in different ranges of W ; a, c) compared with CLEO [9] and VENUS [10] data with statistical (inner error bars) and systematic errors (outer bars), and b) compared with TASSO [5].The TASSO error bars are statistical only.The data points are slightly displaced for clarity.

Figure 5 :
Figure 5: Measured differential cross-section, dσ(γγ → pp)/d| cos θ * |, with statistical (inner bars) and total uncertainties (outer bars) for a) 2.55 < W < 2.95 GeV and b) 2.15 < W < 2.55 GeV.The data are compared with the point-like approximation for the proton (8) scaled to fit the data.The other curves show the pure quark model[2], the diquark model of[11] with the Dziembowski distribution amplitudes (DZ-DA), and the diquark model of[12] using standard and asymptotic distribution amplitudes.

Table 1 :
Number of generated events and corresponding integrated luminosities for the different background processes.Since no events pass the selection cuts, 95% confidence level (CL) upper limits are given for the number of background events expected to contribute to the selected data sample.

Table 3 :
Number of events and cross-sections for | cos θ * | < 0.6.Statistical and systematic uncertainties are also given.A Poisson asymmetric statistical uncertainty, determined within a 68% confidence interval, has been calculated for the last two bins of W with only one selected event.

Table 4 :
Number of data events selected in each bin of | cos θ * | and the measured differential cross-sections dσ(γγ → pp)/d| cos θ * | for three W ranges.The statistical and systematic uncertainties are also given.

Table 5 :
Estimates of relative systematic uncertainties.The total systematic uncertainty is obtained by adding them in quadrature.