Comparison of Deep Inelastic Electron-Photon Scattering Data with the Herwig and Phojet Monte Carlo Models

Deep inelastic electron-photon scattering is studied in the Q**2 range from 1.2 to 30 GeV**2 using the LEP1 data taken with the ALEPH, L3 and OPAL detectors at centre-of-mass energies close to the mass of the Z boson. Distributions of the measured hadronic final state are corrected to the hadron level and compared to the predictions of the HERWIG and PHOJET Monte Carlo models. For large regions in most of the distributions studied the results of the different experiments agree with one another. However, significant differences are found between the data and the models. Therefore the combined LEP data serve as an important input to improve on the Monte Carlo models.


Introduction
The measurement of the hadronic structure function F γ 2 crucially depends on the accurate description of the hadronic final state by Monte Carlo models.The available models do not properly account for all features observed in the data, and therefore, at present, the accuracy of the measurement of F γ 2 is mainly limited by the imperfect description of the hadronic final state by the Monte Carlo models.In previous analyses of the individual LEP experiments [1][2][3][4], it had been shown that there are discrepancies in several distributions of the hadronic final state between the various QCD models and the data.It has also be seen that the data are precise enough to further constrain the models.The purpose of this paper is to combine the ALEPH, L3 and OPAL data to establish a consistent and significant measurement, which can be used to optimise the models.
In this paper the reaction e + e − → e + e − γ ⋆ γ → e + e − hadrons, proceeding via the exchange of two photons, is studied in the single tag configuration, where one scattered electron 3 is detected.The differential cross-section for the deep inelastic electron-photon scattering reaction, e(k) γ(p) → e(k ′ ) γ(p) γ * (q) → e X, where the terms in brackets denote the four-momentum vectors of the particles, is proportional to F γ 2 (x, Q 2 ) [5].Here Q 2 = −q 2 = −(k − k ′ ) 2 and x = Q 2 /2p • q.Experimentally, in the single tag configuration, the value of x is obtained using where W 2 is the hadronic invariant mass squared, and P 2 = −p 2 is neglected in calculating x.
The hadronic structure function F γ 2 receives contributions both from the point-like part and from the hadron-like part of the photon structure.The point-like part can be calculated in perturbative QCD.At low Q 2 the hadron-like part is usually modelled based on the Vector Meson Dominance model.The combined contributions are evolved using the DGLAP evolution equation.
Combining the results of three of the LEP experiments not only reduces the statistical errors compared to the individual results, but the difference between the results also gives a reliable estimate of the systematic, detector dependent, errors.The experimental data are fully corrected for trigger inefficiencies and background has been subtracted.The experimental distributions presented are also corrected for detector effects using different Monte Carlo models, and can directly be compared to the model predictions based on generated quantities only, i.e., without the simulation of the detector response, provided that a well defined set of hadron level cuts defined below is applied.This experimental information serves as a basis for improvements on the models.
A set of variables is chosen to compare the corrected data to the hadron level predictions of the Monte Carlo models.The variables used are: • the reconstructed invariant hadronic mass, W res , within a restricted range in polar angles with respect to the beam axis, • the transverse energy out of the plane defined by the beam direction and the direction 3 Electrons and positrons are referred to as electrons.
1 of the tagged electron, E t,out , • the number of charged particles, N trk , • the transverse momenta of charged particles with respect to the beam axis, p t,trk , • and the hadronic energy flow, 1/N • d E/d η, as a function of the pseudorapidity η = − ln(tan(θ/2)) with respect to the beam axis, where N denotes the number of events.The complete definition of how these variables are calculated is given in Section 3.
For W res , E t,out and N trk the result is a differential event cross-section, which means, the distributions have one entry per event, whereas for p t,trk the distribution is a oneparticle inclusive cross-section.The hadronic energy flow is shown as an average energy flow per event, E/N, where the sum runs over all objects and over all events in a given bin of pseudorapidity.
The analysis presented here is based on data of the individual experiments taken at centre-of-mass energies close to the mass of the Z boson.The approximate ranges in Q 2 and x of the different experiments used in this analysis, and the integrated luminosities are listed in Table 1.The Q 2 ranges are calculated from the requirements on the energy and angle of the deeply inelastically scattered electron, and the x ranges are derived from Eq. 1, using the range in Q 2 and the approximate reach in hadronic invariant mass of 3 < W < 35 GeV.
In the HERWIG model the hard interaction is simulated as e q → e q scattering, where the incoming quark is generated according to a set of parton distribution functions of the photon.The incoming quark is subject to an initial state parton shower which contains the γ → q q vertex, and the outgoing partons undergo final state parton showers as in the case of e + e − annihilations.The hadronisation is based on the cluster model.The initial state parton shower is designed in such a way that the hardest emission is matched to the sum of the matrix elements for the resolved processes, g → q q, q → q g and the point-like γ → q q process.The parton shower uses an angular evolution parameter, and so it obeys angular ordering.For point-like events the transverse momentum of the partons with respect to the direction of the incoming photon is given by perturbation theory.In contrast, for hadron-like events, the photon remnant gets a transverse momentum k t with respect to the direction of the incoming photon, where originally the transverse momentum was generated from a Gaussian distribution.
The HERWIG version used for event simulation is HERWIG5.9+kt , where the label k t denotes that the k t distribution has been altered from the program default.This change is made to improve the agreement for the high-Q 2 region between the ALEPH and OPAL data and the original HERWIG prediction, where the original program is denoted as HERWIG default.The default Gaussian behaviour is replaced by a power-law function of the form dk 2 t /(k 2 t + k 2 0 ) [9] with k 0 = 0.66 GeV.The change is motivated by the observation made in photoproduction studies at HERA [10], where the power-law function gave a better description of the data.This is a good example of how information from two different, but related reactions can be used to improve on a general purpose Monte Carlo program.It is interesting to note that the same value of k 0 is chosen to describe both, photoproduction and deep inelastic electron-photon scattering events.The upper limit of k 2 t in HERWIG+k t is fixed at k 2 t,max = 25 GeV 2 , which is almost the upper limit of Q 2 for the Q 2 region studied, as shown in Table 1.
Based on the distributions presented in Section 5 the fixed cut-off has been replaced by dynamically adjusting the upper limit on k t on an event by event basis [11].In this scheme the maximum transverse momentum k 2 t,max is set to the hardest virtuality scale in the event, which is of order Q 2 .The distributions produced with this procedure are denoted as HERWIG+k t (dyn).
The PHOJET Monte Carlo is based on the Dual Parton Model [12].It is designed for hadron-hadron, photon-hadron and photon-photon collisions, where originally only real or quasi-real photons were considered.It has recently been extended to simulate the deep inelastic electron-photon scattering case, where one of the photons is highly virtual.For the case of deep inelastic scattering the program is not based on the DIS formula using F γ 2 , but the γ ⋆ γ cross-section is calculated from the γγ cross-section by extrapolating in Q 2 on the basis of the Generalised Vector Dominance model.The events are generated from soft and hard partonic processes, where a cut-off of 2.5 GeV on the transverse momentum of the scattered partons in the photon-photon centreof-mass system is used to separate the two classes of events.The hard processes are sub-divided into direct processes, where the photon as a whole takes part in the hard interactions, and resolved processes.In resolved processes either one or both photons fluctuate into a hadronic state, and a quark or gluon of one hadronic state interacts either with the other photon, or with a quark or gluon of the second hadronic state.Also virtual photons can interact as resolved states, however, the parton distribution functions of the photons are suppressed as a function of the photon virtualities.The sum of the processes is matched to the deep inelastic scattering cross-section.Initial state parton showers are simulated with a backward evolution algorithm using the transverse momentum as evolution scale.Final state parton showers are generated with the Lund code JETSET [13].Both satisfy angular ordering implied by coherence effects.The hadronisation is based on the Lund string model as implemented in JETSET.

Experimental Method
Large data sets are generated with the PHOJET and HERWIG+k t programs respectively, for √ s = M Z .All unstable particles with lifetimes of less than 1 ns are allowed to decay in the event generation.In this way the particles of the final state correspond approximately to those actually seen in the detectors.The corresponding integrated luminosities for the PHOJET and HERWIG+k t samples are 831 pb −1 and 683 pb −1 respectively4 .The definitions of the phase space and observables include cuts at generator level both on the events and on the particles within the events.The cuts are chosen such that the individual detectors have good acceptance and therefore detector related uncertainties are expected to be small.To select the events at this stage, the following cuts are applied to the generated hadron level quantities: 1.The energy of the deeply inelastically scattered electron has to be larger than 35 GeV. 2. The polar angle θ tag of the deeply inelastically scattered electron with respect to either beam direction has to be in the ranges 27 − 55 mrad (low-Q 2 region) or 60 − 120 mrad (high-Q 2 region).The two ranges in scattering angles studied correspond to Q 2 ranges of about 1.2 < Q 2 < 6.3 GeV 2 and 5.7 < Q 2 < 30 GeV 2 .3. The events are required to contain no electron with energy of more than 35 GeV and polar angle above 25 mrad with respect to the beam direction in the hemisphere opposite to the one containing the deeply inelastically scattered electron.4. The number of charged particles N trk , calculated by summing over all charged particles which have a transverse momentum p t with respect to the beam axis of more than 200 MeV and polar angles, θ, with respect to the beam axis, in the range 20 < θ < 160 • , has to be greater than or equal to 3. The values chosen closely resemble the acceptance of the tracking detectors.5.The invariant mass W res , calculated by summing over all charged and neutral particles fulfilling p t > 200 MeV and 20 < θ < 160 • , corresponding to |η| < 1.735, has to be larger than 3 GeV.This set of particles and cuts defines the hadron level and the data are corrected to this level.The hadron level predictions of W res and E t,out are calculated from the charged and neutral particles defined above, and the distributions of N trk and p t,trk from the charged particles alone.The only exception is the hadronic energy flow.For the hadronic energy flow, the same event selection has been applied, but all particles are included in the distribution without applying a cut on transverse momentum.Figure 1 shows the differential cross-section dσ/dQ 2 within the cuts listed above and corrected for detector effects.The vertical line roughly separates the low-Q 2 and high-Q 2 regions.Since Q 2 depends on the energy and angle of the deeply inelastically scattered electron there is a slight overlap in Q 2 , but due to the cut on θ tag the two samples are statistically independent.In the kinematic region studied the cross-section prediction of HERWIG+k t is about 40% higher than the prediction based on PHOJET as shown in Figure 1.
To study the experimentally observable distributions at the detector level samples of 60k HERWIG+k t and 120k PHOJET events, which are statistically independent from the samples mentioned above, are passed through the detector simulation programs of the ALEPH, L3 and OPAL collaborations, ensuring that all experiments use identical events.The objects reconstructed after the detector simulation are energy clusters, measured in the electromagnetic and hadronic calorimeters, and tracks, measured in the tracking devices.Identical event selection cuts at the detector level are applied by all experiments, closely reflecting the cuts applied at the hadron level as described above: 1.A cluster of at least 35 GeV is required in one of the small angle electromagnetic luminosity monitors.2. The polar angle with respect to either beam direction of the cluster has to be in the range from 27 − 55 mrad or 60 − 120 mrad.3. The most energetic cluster in the hemisphere opposite to the tagged electron has to have energy less than 35 GeV. 4. At least 3 tracks, fulfilling a set of quality criteria, are required to be observed in the tracking devices with 20 With this strategy it is ensured that the hadron level distributions, which are obtained from the large size samples without detector simulation, and the detector level distributions, which are obtained from the samples of smaller size with a detailed detector simulation for each individual experiment, are statistically independent.Both samples are used in the correction procedure applied to the data described in the next section.

Corrections for Detector Effects
Before a measured quantity can be compared to theoretical predictions or to the results of other measurements it must first be corrected for various detector related effects, such as geometrical acceptance, detector inefficiency and resolution, decays, particle interactions with the material of the detector and the effects of the event and track selections.Figure 2 shows the energy flow, 1/N •d E/d η, predicted by the HERWIG+k t model at the hadron level as well as at the detector level, for the three detectors.The events are entered in the figure such that the deeply inelastically scattered electron is always at negative rapidities, but not shown.Also shown in the figure is the coverage in η of the various sub-detectors used in this analysis.A detailed description of the ALEPH, L3 and OPAL detectors can be found in [16], [17], and [18] respectively.The central region of all the detectors with −1.735 < η < 1.735 is covered by tracking and calorimetry.The forward η > 1.735, and backward η < −1.735 regions are covered by electromagnetic calorimeters for luminosity measurements and typically extend out to η = ±4.3.The electromagnetic calorimeter of L3 between 1.735 < |η| < 3.411 is not used in the present analysis.
For the analysis presented here, the correction of the data to the hadron level is done with multiplicative factors, f , relating the measured value X meas data of a quantity X, such as a bin content, to the corrected value, X corr data using the relation: For distributions, the correction factors are computed bin by bin, e.g. for the energy flow, f is the ratio of the hadron level (lightly shaded) and the detector level (darkly shaded) distributions shown in Figure 2. In this way of correcting the data the assumption is made that within the restricted angular range there is little smearing of the variables between bins, hence a simple correction factor is justified, and therefore no attempt to use an unfolding procedure has been made.The Monte Carlo was used to verify the accuracy of this assumption.Application of this correction results in measurements corrected to a well-defined kinematical region and particle composition, as defined in Section 3.
The correction factors for the energy flow in the low-Q 2 region are shown in Figure 3.The correction factors are near one in the central region of pseudorapidity where identical cuts have been applied and they are similar for the three experiments.However, there is a much larger spread in the region of larger pseudorapidity η > 1.735, where the experiments have different sub-detectors, different angular coverage, and apply different cuts.For ALEPH and OPAL the clusters in the forward detectors are required to have an energy of at least 1 GeV, while for L3 this requirement is at least 4 GeV.This leads to larger correction factors for L3 in that region.In addition, for L3, the clusters measured in the forward detectors on the side of the tagged electron, i.e. at η < −2, are not considered in this study.Another difference in the treatment of clusters measured in the forward detectors is that OPAL uses a correction function, obtained from Monte Carlo studies of the detector response to hadrons, to correct for losses in the measurement of the hadronic energy in the forward detector, thereby also reducing the correction factors.It should be stressed that for most of the distributions used for the comparisons in the following section, only the central detector part for polar angles with respect to the beam axis in the range 20 < θ < 160 • , that is |η| < 1.735, is exploited.The only exception is the energy flow where no cuts on the angle are applied.
The correction factors in the low-Q 2 region for the other chosen variables and for the ALEPH, L3 and OPAL experiments using the HERWIG+k t and PHOJET models, are shown in Figures 4 and 5.The quoted errors of the correction factors are the combined statistical uncertainties of the generated hadron level and the simulated detector level quantities.While the differences between the correction factors obtained from HERWIG+k t and PHOJET are very small for the energy flow, they can vary significantly for other variables.For example, in the case of OPAL, for W res the HERWIG+k t correction factors are on average about 20% higher than the factors obtained with PHOJET.The correction factors for the low-Q 2 and high-Q 2 regions typically differ by less than 10%.

Corrected Data Comparisons
The discussion of the comparison is subdivided into three parts.First the corrected data from the individual experiments are compared to each other and to the Monte Carlo models.In this comparison only statistical errors are used and no attempt has been made to obtain estimates of systematic errors for the individual experiments.Based on the above comparison a modified version of the HERWIG+k t model has been developed which is described next.Finally the data are combined and compared to the Monte Carlo models.In the combination of the data the spread of the experiments is used as an estimate of the systematic uncertainty of the measured distributions.The numerical results are listed in Tables 2-12 and can be obtained in electronic form [19].

Data from individual experiments
The  show the corrected differential cross-sections, calculated in the kinematical range described in Section 3, for the data compared with the HERWIG+k t and PHOJET predictions.Figures 6,8,10 and 12 show the data on a linear and logarithmic scale, corrected with the HERWIG+k t model, while Figures 7,9,11 and 13 show the same data corrected with PHOJET.For all distributions the errors shown are the quadratic sum of the statistical errors of the measured quantity and the statistical errors of the correction factors.As an example, Table 2 shows the results for the W res distribution for the three experiments, listing the statistical error on the data and the statistical error on the correction factors, f , separately.
The experimental results from ALEPH, L3 and OPAL agree with each other within errors for large regions in most of the variables studied.However, there are also regions where they significantly differ from each other, for example, in the region of W res < 10 GeV, E t,out < 5 GeV, for low charged multiplicities and low p t,trk in the low-Q 2 region.The level of agreement also depends on the Monte Carlo model used for correction.The agreement between the experiments is slightly better for the data corrected with PHOJET, than for the data corrected with HERWIG+k t .In the combination of the data, the differences between the measured distributions of the different experiments will be used as an estimate of the systematic error.
There are significant differences between the Monte Carlo distributions and the data, particularly in the low-Q 2 region.The PHOJET distributions lie consistently below HERWIG+k t , especially at the low end of the distributions, while the agreement of the tails in the high-Q 2 region is reasonable.In general the PHOJET predictions agree reasonably well with the data, corrected with PHOJET, for large values of the variables.However, there are large differences in the low-Q 2 region between the data, corrected with HERWIG+k t , and the HERWIG+k t predictions in all the distributions.For the W res , E t,out and p t,trk distributions (except in the region of low values, where the data are inconsistent) the differences between the experiments are much less than the HERWIG+k t −PHOJET differences.
Figures 14 and 15 show the corrected energy flow as a function of pseudorapidity in the low-Q 2 and high-Q 2 region for the individual experiments compared to the HERWIG+k t and PHOJET predictions.Following the energy flow analyses at HERA, the statistical errors for the energy flow are taken as √ E 2 /N.In Figure 14 the data are corrected with the HERWIG+k t model and in Figure 15 with the PHOJET model.In the central region of the detector, between −1.5 < η < 1.5, the data lie between the HERWIG+k t and the PHOJET predictions in the low-Q 2 region, whereas in the high-Q 2 region the two MC predictions lie closer to each other than the data of OPAL and ALEPH.These corrected energy flows in the central region are stable against changes in the requirements for the hadronic final state, e.g.changes to the quality cuts of the accepted tracks and clusters.
For η > 1.5 the data from the various experiments vary much more than the statistical errors.Furthermore, it has been found that the corrected energy flows are not very stable against variations of the event selection cuts such as the anti-tag condition or a maximum energy cut on the energy deposited in the forward detectors.The observed changes are much larger at the detector level than at the hadron level, leading to the conclusion that the modelling of the energy response to the hadronic energy in the forward region is poor.This can be understood since the sub-detectors covering this region have poor energy resolution for hadronic energies and no particle identification.Therefore the uncertainty of the hadronic energy flow in the forward region is larger than is indicated by the spread of the experiments and it is difficult to draw firm conclusions on the description by the Monte Carlo models.However, the data appear to lie consistently below the prediction of either generator.

Modification of the HERWIG model
As discussed above, the distribution of transverse momentum k t of the photon remnant with respect to the direction of the incoming photon has been altered, motivated by the findings in photoproduction at HERA.At LEP, the modification was initially studied as a possible improvement of the agreement between the HERWIG prediction and the high-Q 2 data of OPAL and ALEPH.While HERWIG+k t seems to reasonably describe the data in the high-Q 2 region for low and high values of E t,out , it overestimates the distribution for medium values of E t,out , as shown in Figure 8.Even though the description of the energy flow is improved with the HERWIG+k t generator, it fails to accurately describe the data over a wide range of the pseudorapidity η.
While the fixed limit of k 2 t,max = 25 GeV 2 is adequate for the high-Q 2 region, in the low-Q 2 region it introduces too much transverse momentum, which is most clearly seen in the transverse momentum distribution of the tracks in Figure 12.Therefore the dynamical adjustment, HERWIG+k t (dyn), as discussed in Section 2, has been introduced.As will be seen in the next section, when comparing with the combined LEP data, this change leads to an improved description of the data also for the low-Q 2 region.

Combined data
In this section the results from the individual experiments discussed in Section 5.1 are combined.In large ranges of the phase space the individual results agree within the statistical precision quoted, however there are also significant differences as discussed above.These differences are expected because the previous analyses of the individual LEP experiments [1][2][3][4] showed that the systematic errors, which are not included above, dominate.Because the Monte Carlo models do not sufficiently well resemble the data, evaluating the experimental systematic errors of the measurements based on these models will not be very reliable.On the other hand a combined result is desirable in order to serve as a guidance for the model builders to improve on their Monte Carlo programs.Therefore, in a first attempt to make a combined measurement, the experimental systematic error is taken from the spread of the measured results, wherever they are significantly different on the basis of the statistical error alone.In this case the purely statistical error of the combined result is inflated as discussed below.
The combined distributions from ALEPH, L3 and OPAL are shown in comparison to PHOJET and various predictions from HERWIG in .The measured values are listed in Tables 3-12.The combination of the data follows the procedure recommended by the Particle Data Group in section 4.2.2 of [20].Since this is a crucial point of the analysis and because specific choices have to be made in the combination, the procedure is briefly discussed below.
The combined value for a given bin is calculated as the weighted average of the measurements of the individual experiments using the statistical errors to calculate both the weights and the error of the combined value.To calculate the average bin content x and its error σ x from the individual contents x i and their statistical errors σ i the following procedure is applied.The average content is calculated from x = w i • x i / w i , using the weights w i = 1/σ 2 i .The sum runs over all experiments i = 1, . . ., N ex with non-zero entries in that bin, and the error σ x is taken to be 1/ √ w i .The χ 2 of the average is calculated from If in the tail of a distribution, for example, some experiments measure zero in a particular bin, then x and χ 2 are scaled by the ratio of the number of experiments with nonzero entries to the number of experiments being averaged for that bin.If χ 2 is larger than N ex − 1, the error is increased by a factor S = χ 2 /(N ex − 1), thereby taking the spread of the experiments as an estimate of the experimental systematic uncertainty.Finally, to obtain the errors quoted the uncertainties due to the correction factors are added in quadrature.Since the same data sets are used by each experiment to calculate the correction factors, the corresponding errors are strongly correlated between the experiments.To take this into account, the error on the correction factors is included by taking the smallest quoted error of the individual experiments as an estimate of this systematic error, which is assumed to be 100% correlated amongst the experiments.
Also in the combination of the results of the individual experiments the 1/N •d E/d η distribution is treated slightly differently from the other distributions.There is a large scatter in the measured values of the different experiments, especially in the forward region, as can be seen from Figures 14 and 15.As a consequence there is also a large scatter in the scale factors listed in Tables 11 and 12.To avoid the combination procedure manufacturing artificially small errors for bins where the measurements happen to coincide, the scale factor applied to obtain the combined measurement is taken as the average of the individual scale factors from three neighbouring bins centered around the bin under study.
It is apparent from Figures 16,18,20,22 and 24, that the new HERWIG+k t (dyn) with the dynamical cut-off lies much closer to the low-Q 2 averaged data than the version of HERWIG+k t using the fixed cut-off.However, the difference between HERWIG+k t (dyn) and HERWIG default is rather small.
In the case of the energy flow none of the various HERWIG models is able to accurately describe the data.This suggests that even though the new HERWIG+k t (dyn) better describes most of the data distributions, the energy flow is still not well understood.
The PHOJET model,Figures 17,19,21,23 and 25, describes the data reasonably well in both regions of Q 2 , but underestimates the cross-section near the lower limit of the distributions.This is understood as a consequence of the high transverse momentum cut-off of 2.5 GeV for the scattered partons in the hard scattering matrix element.Below this cut-off only soft events are generated.Due to the different Q 2 dependences of the soft and hard components in PHOJET, this leads to a suppression of the low-W events [21].
The energy flow corrected with HERWIG+k t and with PHOJET Figures 24 and 25, mostly agree with each other within errors, except in a few bins in the forward region.In this region the data corrected with PHOJET lie below the data corrected with HERWIG+k t in the region of the peak at η ≃ 2 and above the flow corrected with HERWIG+k t in the region of η > 2.5.

Conclusion
For the first time the results of deep inelastic electron-photon scattering from three of the LEP experiments have been combined and compared to predictions from the PHOJET and the HERWIG models.It is found that the data from the ALEPH, L3 and OPAL experiments agree within statistical errors except near the edges of the distributions.Where the spread is larger than expected from the statistical errors, as, for example, for low charged multiplicities, this difference is taken as an estimate of the detector dependent systematic uncertainty of the measurement.
In the comparison of the data with the HERWIG+k t model the most striking discrepancy is seen in the distributions of the low-Q 2 region, where the HERWIG+k t model systematically overestimates the data.This discrepancy is found to be mainly due to the fixed cut-off for the intrinsic transverse momentum of the quarks in the photon in the HERWIG+k t model.By dynamically adjusting the cut-off according to the kinematics of the individual event in the HERWIG+k t (dyn) model the description of the data is significantly improved, particularly in the low-Q 2 region.
The PHOJET model describes the data reasonably well in both regions of Q 2 , but underestimate the cross-section near the lower limit of the distributions, due to the high transverse momentum cut-off for the scattered partons in the hard scattering matrix element.
The energy flow of the data lies between the predictions of the HERWIG and PHOJET Monte Carlo models in the central regions of the detectors.In the forward region the Monte Carlo predictions lie systematically above the data.It should be noted, however, that it is difficult to assess the systematic errors in this region because of the poor resolution of the hadronic energy measured in the electromagnetic luminosity monitors.
The method of combining the data of several of the LEP experiments has proven useful to detect shortcomings of Monte Carlo models in the description of these data.For the HERWIG Monte Carlo this investigation already demonstrated that the changed distribution of transverse momentum k t of the photon remnant with respect to the direction of the incoming photon, the HERWIG+k t (dyn) model, gives a better description of the LEP data.As the data distributions are corrected to the hadron level, they can be directly compared to the predictions of the Monte Carlo models, without the need of detector simulation, and thus can be used more easily by Monte Carlo model builders.0.181 0.147 ± 0.045 1.854 0.151 0.138 ± 0.031 0.710 9.0 − 10.0 0.114 0.083 ± 0.015 0.596 0.088 0.134 ± 0.028 0.000 10.0 − 11.0 0.084 0.072 ± 0.022 0.280 0.073 0.018 ± 0.012 0.000 11.0 − 12.0 0.051 0.018 ± 0.009 0.515 0.057 0.055 ± 0.016 0.000 12.0 − 13.0 0.048 0.013 ± 0.013 1.410 0.034 0.000 ± 0.000 0.000 13.0 − 14.0 0.024 0.036 ± 0.018 0.509 0.019 0.000 ± 0.000 0.000 14.0 − 15.0 0.030 0.012 ± 0.017 0.174 0.019 0.000 ± 0.000 0.000  10: The combined differential cross-section dσ/dp t,trk calculated in the kinematical range defined in the text for the low-Q 2 and high-Q 2 regions.The data are corrected with the PHOJET model.The PHOJET prediction (MC) and the scale factor S are shown in addition.HERWIG corrected low Q       HERWIG corrected low Q

Figure 1 :
Figure1: The differential cross-section, dσ/dQ 2 , observed in the two ranges in scattering angles θ tag studied, compared to the predictions of the HERWIG+k t and PHOJET models.The cross-sections are given in the kinematical ranges described in the text.The errors are statistical only.

Figure 2 :Figure 3 :
Figure 2: The HERWIG+k t energy flow, 1/N • d E/d η, at the hadron level (lightly shaded) as well as on detector level (darkly shaded), as measured by the three detectors.The coverage of the different detector components are indicated by dot-dashed lines.In the case of L3 the ALR is not used in this analysis.

Figure 4 :Figure 5 :
Figure 4: The HERWIG+k t correction factors, f , for the ALEPH, L3 and OPAL W res , E t,out , N trk and p t,trk distributions for the low-Q 2 region.The symbols are slightly displaced for better visibility.

Figure 6 :
Figure 6: The W res distributions from ALEPH, L3 and OPAL for the low-Q 2 (left) and high-Q 2 region (right), corrected with the HERWIG+k t model on a linear scale (top) and on a log scale (bottom).

Figure 7 :
Figure 7: The W res distributions from ALEPH, L3 and OPAL for the low-Q 2 (left) and high-Q 2 region (right), corrected with the PHOJET model on a linear scale (top) and on a log scale (bottom).

Figure 10 :
Figure 10: The N trk distributions from ALEPH, L3 and OPAL for the low-Q 2 (left) and high-Q 2 region (right), corrected with the HERWIG+k t model on a linear scale (top) and on a log scale (bottom).

Figure 11 :
Figure 11: The N trk distributions from ALEPH, L3 and OPAL for the low-Q 2 (left) and high-Q 2 region (right), corrected with the PHOJET model on a linear scale (top) and on a log scale (bottom).

Figure 12 :Figure 13 :
Figure 12: The p t,trk distributions from ALEPH, L3 and OPAL for the low-Q 2 (left) and high-Q 2 region (right), corrected with the HERWIG+k t model on a linear scale (top) and on a log scale (bottom).

Figure 14 :Figure 15 :
Figure 14: The hadronic energy flow from ALEPH, L3 and OPAL for the low-Q 2 (left) and high-Q 2 region (right), corrected with the HERWIG+k t model.The shaded band indicates the forward region of the experiments.

Figure 16 :Figure 17 :
Figure 16: The combined W res distribution from ALEPH, L3 and OPAL for the low-Q 2 (left) and high-Q 2 region (right), corrected with the HERWIG+k t model on a linear scale (top) and on a log scale (bottom).The data are compared to three different model assumptions of the HERWIG+k t model.

Figure 20 :
Figure 20: The combined N trk distribution from ALEPH, L3 and OPAL for the low-Q 2 (left) and high-Q 2 region (right), corrected with the HERWIG+k t model on a linear scale (top) and on a log scale (bottom).The data are compared to three different model assumptions of the HERWIG+k t model.

Figure 21 :
Figure 21: The combined N trk distribution from ALEPH, L3 and OPAL for the low-Q 2 (left) and high-Q 2 region (right), corrected with and compared to the PHOJET model.

Figure 22 :Figure 23 :
Figure 22: The combined p t,trk distribution from ALEPH, L3 and OPAL for the low-Q 2 (left) and high-Q 2 region (right), corrected with the HERWIG+k t model on a linear scale (top) and on a log scale (bottom).The data are compared to three different model assumptions of the HERWIG+k t model.

Table 1 :
The approximate ranges in Q 2 and x of the different experiments and the integrated luminosities L int used in this analysis.
• < θ < 160 • , and with p t,trk of at least 200 MeV. 5.The invariant mass, W res , calculated from all tracks and clusters with p t > 200 MeV and 20 < θ < 160 • is required to be greater than 3 GeV.These objects define the detector level.They are used for the W res , E t,out , N trk and p t,trk distributions.The only exception is the energy flow 1/N • d E/d η, where again no cut on the transverse momentum has been applied.

Table 2 :
L3OPAL W res[GeV]dσ/dW res [pb/GeV] dσ/dW res [pb/GeV] dσ/dW res [pb/GeV] The individual differential cross-section dσ/dW res in the low-Q 2 region for the ALEPH, L3 and OPAL data at √ s = 91 GeV, calculated in the kinematical range defined in the text.The data have been corrected with HERWIG+k t .The first error listed is the statistical error on the data only, the second one is the statistical error arising from the correction factors, f .

Table 3 :
The combined differential cross-section dσ/dW res calculated in the kinematical range defined in the text for the low-Q 2 and high-Q 2 regions.The data are corrected with the HERWIG+k t model.The HERWIG+k t prediction (MC) and the scale factor S are shown in addition.

Table 4 :
The combined differential cross-section dσ/dW res calculated in the kinematical range defined in the text for the low-Q 2 and high-Q 2 regions.The data are corrected with the PHOJET model.The PHOJET prediction (MC) and the scale factor S are shown in addition.

Table 5 :
The combined differential cross-section dσ/dE t,out calculated in the kinematical range defined in the text for the low-Q 2 and high-Q 2 regions.The data are corrected with the HERWIG+k t model.The HERWIG+k t prediction (MC) and the scale factor S are shown in addition.

Table 6 :
The combined differential cross-section dσ/dE t,out calculated in the kinematical range defined in the text for the low-Q 2 and high-Q 2 regions.The data are corrected with the PHOJET model.The PHOJET prediction (MC) and the scale factor S are shown in addition.

Table 7 :
The combined differential cross-section dσ/dN trk calculated in the kinematical range defined in the text for the low-Q 2 and high-Q 2 regions.The data are corrected with the HERWIG+k t model.The HERWIG+k t prediction (MC) and the scale factor S are shown in addition.

Table 8 :
The combined differential cross-section dσ/dN trk calculated in the kinematical range defined in the text for the low-Q 2 and high-Q 2 regions.The data are corrected with the PHOJET model.The PHOJET prediction (MC) and the scale factor S are shown in addition.
Figure 8: The E t,out distributions from ALEPH, L3 and OPAL for the low-Q 2 (left) and high-Q 2 region (right), corrected with the HERWIG+k t model on a linear scale (top) and on a log scale (bottom).
Figure 9: The E t,out distributions from ALEPH, L3 and OPAL for the low-Q 2 (left) and high-Q 2 region (right), corrected with the PHOJET model on a linear scale (top) and on a log scale (bottom).
Figure 18: The combined E t,out distribution from ALEPH, L3 and OPAL for the low-Q 2 (left) and high-Q 2 region(right), corrected with the HERWIG+k t model on a linear scale (top) and on a log scale (bottom).The data are compared to three different model assumptions of the HERWIG+k t model.Figure19:The combined E t,out distribution from ALEPH, L3 and OPAL for the low-Q 2 (left) and high-Q 2 region(right), corrected with and compared to the PHOJET model.