{'fraction_non_alphanumeric': 0.032247133901846727, 'fraction_numerical': 0.018862596713164287, 'mean_word_length': 4.633052330205185, 'pattern_counts': {'":': 0, '<': 12, '<?xml version=': 0, '>': 2, 'https://': 0, 'lorem ipsum': 0, 'www.': 0, 'xml': 0}, 'pii_count': 0, 'substrings_counts': 0, 'word_list_counts': {'cursed_substrings.json': 0, 'profanity_word_list.json': 0, 'sexual_word_list.json': 0, 'zh_pornsignals.json': 0}}
Production of pizero and eta mesons at large transverse momenta in pp and pBe interactions at 530 and 800 GeV/c

We present results on the production of high transverse momentum pizero and eta mesons in pp and pBe interactions at 530 and 800 GeV/c. The data span the kinematic ranges: 1<p_T<10 GeV/c in transverse momentum and 1.5 units in rapidity. The inclusive pizero cross sections are compared with next-to-leading order QCD calculations and to expectations based on a phenomenological parton-k_T model.


I. INTRODUCTION
The study of inclusive single-hadron production at large transverse momentum (p T ) has been a useful probe in the development of perturbative quantum chromodynamics (PQCD) [1,2]. Early in the evolution of the parton model, a departure from an exponential dependence of particle production at lower p T was interpreted in terms of the onset of interactions between pointlike constituents (partons) contained in hadrons. Large p T is a regime where perturbative methods have been applied to QCD to provide quantitative comparisons with data. Such comparisons yield information on the validity of the PQCD description, and on parton distribution functions of hadrons and fragmentation functions of partons.
This paper reports high-precision measurements of the production of π 0 and η mesons with large p T . The π 0 production cross sections are compared with next-to-leading order (NLO) PQCD calculations [3]. As illustrated in a previous publication [4], our data (for both inclusive π 0 and direct-photon production) are not described satisfactorily by the available NLO PQCD calculations, using standard choices of parameters. Similar discrepancies have been observed [5] between conventional PQCD calculations and other measurements of highp T π 0 and direct-photon cross sections (see also [6][7][8]). The origin of these discrepancies may be attributed to the effects of initial-state soft-gluon radiation. Such radiation generates transverse components of initial-state parton momenta, referred to below as k T [9]. Evidence of significant k T in various processes, and a phenomenological model for incorporating its effect on calculated high-p T cross sections, have been extensively discussed in Refs. [4,5]; recent studies of the photoproduction of direct photons at HERA may provide additional insights [10][11][12][13][14]. The inadequacy of NLO PQCD in describing k T -sensitive distributions has been discussed in Ref. [15]. In this paper, we follow the phenomenological prescription of Ref. [5] when comparing calculations with our π 0 data. We also present cross sections for η meson production at large-p T . As might have been expected from earlier measurements (see, e.g., [16]), η production, relative to π 0 production, shows little dependence on p T or on center of mass rapidity (y cm ).

II. THE EXPERIMENTAL SETUP
The E706 experiment at Fermilab was designed to measure direct-photon production at high-p T , and to investigate the structure of events containing direct photons. The data collection phase of the experiment spanned three fixed-target running periods, and included a relatively low statistics commissioning run in 1987-88 [16][17][18][19], and primary data runs in

C. Calorimetry
The central element of the E706 apparatus was the finely-segmented liquid argon electromagnetic calorimeter (EMLAC) used to detect and measure electromagnetic showers [27,28].
The EMLAC had a cylindrical geometry with an inner radius of 20 cm and an outer radius of 160 cm. It was divided into four mechanically independent quadrants, which were further subdivided electronically to create octants. The calorimeter had 33 longitudinal cells read out in two sections: an 11 cell front section (≈8.5 radiation lengths) and a 22 cell back section (≈18 radiation lengths). This front/back split was used for measuring the direction of incidence of showering particles, for discriminating between electromagnetic and hadronic showers and for resolving closely separated electromagnetic showers. The longitudinal cells consisted of 2 mm thick lead cathodes (the first cathode was constructed of aluminum), double-sided copper-clad G-10 radial (R) anode boards, followed by 2 mm thick lead cathodes and double-sided copper-clad G-10 azimuthal (Φ) anode boards. There were 2.5 mm argon gaps between each of these layers in a cell. The physical layout is illustrated in Fig. 3.
The copper-cladding on the anode boards was cut to form strips. Signals from corresponding strips from all R (or Φ) anode boards in the front (or back) section were jumpered together. The copper-cladding on the radial anode boards was cut into concentric strips centered on the nominal beam axis. The width of the strips on the first R board was 5.5 mm. The width of the R strips on the following R boards increased slightly so that the radial geometry was projective relative to the target, which was located 9 m upstream of the EMLAC. The azimuthal readout was subdivided at a radius of 40 cm into inner and outer segments, with each inner Φ strip subtending an azimuthal angle of π/192 radians, and outer Φ strips covering π/384 radians. Subdivision of the azimuthal strips in the outer portion of the detector improved both the position and energy resolution for showers reconstructed in this region. It also reduced R-Φ correlation ambiguities from multiple showers in the same octant of the calorimeter.
Data acquisition and trigger-signal processing for the EMLAC was based upon the FNAL RABBIT system [29]. To achieve tolerable deadtimes, the zero suppression features of this system were used extensively during the experiment's 1987-88 commissioning run. However, zero suppression limited our ability to understand the effects of out-of-time photoninduced showers and the tails of hadron-induced showers, thereby compromising efforts to understand the detailed response of the detector. Consequently, FASTBUS modules (the ICBM-Intelligent Control and Buffering Module-and the Wolf interface [30,31]) were developed by E706 to replace the original, CDF-designed, MX readout controllers. These FASTBUS modules enabled us to eliminate zero suppression during the experiment's two primary data runs.
The apparatus also included two other calorimeters: a hadronic calorimeter (HALAC) located downstream of the EMLAC in the same cryostat, and a steel and scintillator calorimeter (FCAL), positioned further downstream, to increase coverage in the very forward region.
The HALAC had 53 longitudinal cells read out in 2 sections: a 14 cell front section (≈2 interaction lengths), and a 39 cell back section (≈6 interaction lengths). Each cell consisted of read-out planes separated by 3 mm argon gaps and a 2.5 cm thick steel plate.
The FCAL acceptance covered the beam hole region of the EMLAC. It was split into three longitudinally similar sections. Each section was composed of alternating layers of 1.9 cm thick steel absorber plates and 4.8 mm thick acrylic scintillator sheets. The distance between steel plates was 6.9 mm. The downstream module contained 32 steel absorber plates and 33 scintillator sheets; the other two modules were comprised of 28 steel absorber plates and 29 scintillator sheets. Together, the three modules constituted ≈10.5 interaction lengths of material.

D. Muon identification
The E672 muon spectrometer, consisting of a toroidal magnet, shielding, scintillators, and proportional wire chambers, was deployed immediately downstream of the FCAL. The combined Meson West Spectrometer was triggered on high-mass muon pairs in order to investigate the hadroproduction of J/ψ, ψ(2s), χ c , and B mesons [32][33][34][35]. E706 and E672 collected data simultaneously and shared trigger logic, the data acquisition system, and event reconstruction programs. Data collected with the dimuon trigger were also used for several technical studies-for example, the J/ψ −→ µ + µ − signal was used to calibrate the momentum scale of the tracking system.

E. Triggering
The E706 trigger selected events yielding high transverse momentum showers in the EM-LAC. The selection process involved four stages: beam and interaction definitions, pretrigger requirements, and the final trigger requirements. Beam particles were detected using a hodoscope consisting of three planes (arranged in X, Y and U views) of scintillator, located ≈2 m upstream of the target region. Each plane contained 12 scintillator strips. The widths of the strips varied from 1 mm in the central region, to 5 mm along the edges. The edges of individual strips in each plane were overlapped to avoid gaps in acceptance. A beam signal was generated by a coincidence of signals from counters in at least two of the three hodoscope planes. A beam1 signal required that less than two hodoscope planes detected two or more isolated clusters of hits in coincidence with beam. This beam1 requirement rejected events with multiple beam particles incident upon the target. A plane of four scintillation counters (referred to as beam hole counters) arranged to produce a 0.95 cm diameter central hole was located downstream of the beam hodoscope and used to reject interactions initiated by particles in the beam halo.
Two pairs of scintillation counters were mounted on the dipole analysis magnet, one pair upstream and the other downstream of the magnet. Each pair had a central hole that allowed non-interacting beam particles to pass through undetected. An interaction was defined as a coincidence between signals from at least two of these four interaction counters.
To minimize potential confusion in the EMLAC due to out-of-time interactions, a filter was used to reject interactions that occurred within 60 ns of one other.
For those interactions that satisfied both the beam1 and interaction definitions, the p T deposited in various regions of the EMLAC was evaluated by weighting the energy signals from the fast outputs of the EMLAC R channel amplifiers by ≈sin θ i , where θ i is the polar angle that the i th strip subtends with respect to the nominal beam axis. The pretrigger hi requirement was satisfied when the p T detected either in the inner 128 R channels or the outer R channels of any octant was greater than the threshold value of ≈2 GeV/c. A pretrigger signal was issued only if the signals from a given octant satisfied the pretrigger requirement, there was no evidence in that octant of substantial noise or significant p T attributable to an earlier interaction, and there was no incident beam halo muon detected by the walls of scintillation counters surrounding the hadron shield upstream of the spectrometer. The pretrigger signal latched data from the various subsystems until a final trigger decision was made.
Localized trigger groups were formed for each octant by clustering the R channels into 32 groups of 8 channels. Each of the adjacent pairs of 8 channel groups (groups 1 and 2, 2 and 3, . . . , 31 and 32) formed a local group of 16 strips. If the p T detected in any of these groups of 16 was above a specified high (or low) threshold, then a local hi (or local lo) signal was generated for that octant. The single local hi (single local lo) trigger required a local hi (local lo) signal from an octant that also satisfied the pretrigger hi. The local hi (local lo) threshold was ≈3.5 GeV/c (≈2 GeV/c). The single local lo trigger was prescaled by a factor of ≈200.
In addition to these high-p T triggers, prescaled samples of low-bias beam, interaction, and pretrigger events were also recorded. The prescale factors for these triggers were typically set at 15 6 , 15 5 , and 15 3 , respectively. These low-bias triggers constituted ≈10% of the events recorded.
III. π 0 AND η ANALYSIS The data sample used in this analysis corresponds to an integrated luminosity of 6.8 (1.2) events/pb for 530 GeV/c pBe (pp) interactions, and 6.5 (1.2) events/pb at 800 GeV/c.
The following subsections describe the analysis procedures and methods used to correct the data for losses resulting from inefficiencies and selection biases.

A. Event reconstruction
The charged-track reconstruction algorithm produced track segments upstream of the dipole magnet using information from the SSDs, and downstream of the magnet using information from the PWCs and STDCs. These track segments were projected to the center of the magnet, and linked to form the final tracks, and the interaction vertex. The chargedtrack reconstruction and vertex-finding methodology are described in Ref. [20]. To reduce backgrounds from hadronic showers, only showers with at least 20% of their shower energy deposited in the front part of EMLAC (E FRONT /E TOTAL > 0.2) were accepted as photon candidates. To correct the cross sections for this requirement, as well as for other larger effects including resolution-smearing and reconstruction losses, a full simulation of the showers in the calorimeter (described below) was employed. In addition, γγ combinations were considered as π 0 or η candidates only when the two photons were detected in the same octant (to simplify subsequent analysis of the trigger response), and only those combinations with energy asymmetry [A γγ ≡ |E γ 1 − E γ 2 |/(E γ 1 + E γ 2 )] less than 0.75 were considered (to reduce uncertainties due to low energy photons).
The γγ invariant mass distribution in the π 0 and η regions for photon pairs that satisfied the above requirements is shown in Fig. 4 for representative low p T intervals. A π 0 candidate was defined as a combination of two photons with invariant mass, M γγ [38], in the range 100 MeV/c 2 < M γγ < 180 MeV/c 2 . An η candidate was defined as a two-photon combination in the range 450 MeV/c 2 < M γγ < 650 MeV/c 2 . To evaluate production cross sections, combinatorial background in the π 0 and η regions was evaluated as follows. Sideband regions were defined to cover an equivalent mass range of the π 0 and η peak regions (using the same acceptance criteria as for the peak regions). The p T and rapidity distributions from these side bands were then subtracted from the corresponding distributions within the π 0 and η mass ranges to obtain the respective signals. This technique is appropriate as long as the combinatorial background depends approximately linearly dependence on M γγ . At low p T (below ≈ 2 GeV/c), the shape of the combinatorial background in the signal regions is not linear, and a more sophisticated fitting procedure was used to evaluate the background.
The γγ mass distributions were fitted to Gaussians for signal plus second-order polynomials in M γγ to represent the background. The combinatorial background in the peak regions was then defined using the resultant parameters of the fit, and the signals defined as the differences between the totals and the fitted backgrounds.
For the cross section measurements, the signals have been corrected for losses due to the energy asymmetry cut and the branching fractions [39] for the γγ decay modes. The correction for losses due to the conversion of one or both of the photons into e + e − pairs was evaluated by projecting each reconstructed photon from the event vertex to the reconstructed position in the EMLAC. The number of radiation lengths of material traversed along the photon path was calculated on the basis of a detailed description of the detector. The photon conversion probability was then evaluated and used to account for losses from conversion.

C. Trigger response
As mentioned previously, the single local hi and single local lo trigger decisions were based upon depositions of p T in the EMLAC within groups of 16 contiguous radial strips. For each π 0 or η candidate, a probability to satisfy the trigger was defined, based upon energy deposition in the entire octant, as P = 1 − (1 − p i ), where p i is the efficiency of the i th trigger group in the octant containing the candidate [40]. The inverse of this probability was applied as a trigger weight to each meson candidate. To avoid excessively large trigger weights, meson candidates with trigger probabilities of < 0.1 were excluded from further consideration. The correction for losses from this requirement was determined from Monte Carlo, and absorbed into the reconstruction efficiency.
The efficiency of the pretrigger hi was determined in a manner similar to that used for the local triggers [40]. Additional details on the trigger can be found in Refs. [20,41,42].
The cross sections presented in this paper reflect composite measurements, utilizing a combination of results from the interaction, pretrigger hi, single local lo, and single local hi triggers. The p T spectra (corrected for only prescale factors) from a representative sample of these triggers are shown in Fig. 5. The transition points chosen between the high and low threshold triggers were determined by comparing the fully corrected results from each trigger, and were different for π 0 and η mesons and also depended on rapidity.

D. Rejection of beam halo muons
Spurious triggers were produced by muons in the beam halo that radiated energy in the electromagnetic calorimeter in random coincidence with an interaction in the target. Particularly in the outer regions of the EMLAC, this energy was recorded as a high-p T deposition that satisfied the local trigger requirements. This occurred much more frequently in data from the 530 GeV/c secondary beam than for the primary 800 GeV/c beam, due to the absence of an upstream interaction target in the latter case. To reduce this background, the pretrigger logic relied on signals from the veto walls of scintillator counters to reject events associated with such muons in the beam halo. In the off-line analysis, we employed expanded requirements on the latched veto-wall signals, the direction of reconstructed showers [28], the shower shape (halo muon-induced showers have a different shape than photon or electron-induced showers), and the total p T (balance) in the event. For the latter, we calculated the net p T of the photons and charged particles which, based upon their initial trajectories, would intercept the EMLAC in the transverse plane within the 120 • sector opposite the meson candidate. In interactions which generate a high-p T meson, the ratio of this "away-side" p T (P away T ) to that of the meson p T should be near unity. However, for events triggered by showers from beam halo muons, P away T /p T should be near zero, since the interaction in random coincidence with the beam halo muon is typically a soft (low p T ) interaction. Candidates with P away T /p T < 0.3 were considered likely to be due to beam halo muons and were rejected.
To illustrate the effect of the above off-line requirements, the γγ invariant mass distribution, both before and after application of the rejection criteria, is shown in Fig. 6 for the 530 and 800 GeV/c data, for p T > 7.0 GeV/c. The large muon-induced background at low γγ mass values in the 530 GeV/c data is due to the occasional transverse splitting of the muon-induced showers into two closely-separated photon candidates. This happens because the reconstruction software assumes that the showers originate from the target-region rather than from the beam halo. The 800 GeV/c data have very few muon-induced triggers, and is consequently not affected very strongly by these rejection criteria.
The impact of these rejection criteria on the physics signal was checked using more restrictive selection criteria to define a pure sample of γγ pairs. The fraction of signal lost by the application of each of the muon-rejection requirements determined a correction to the cross section. The product of the correction factors for muon rejection corresponds to an increase of ≈8% in the cross secton at p T = 4 GeV/c, and ≈10% at p T = 7 GeV/c, for the 530 GeV/c beam data. (The corrections were smaller for the 800 GeV/c data because of cleaner beam conditions.)

E. Detector simulation
The Meson West spectrometer was modeled using a detailed geant [43] simulation.
Because the full simulation of electromagnetic showers requires extensive computing time, we developed a hybrid approach using geant-tracking through the magnetic spectrometer and in the initial stages of the shower development in the calorimeter. We used the standard geant algorithms for tracking particles with energies above 10 MeV, below which we relied on an empirical parameterization for the deposition of energy in the EMLAC [44]. This cutoff was selected to be at the point at which bremsstrahlung still dominates the energy loss in lead, and led to significant improvement in processing speed. In doing this, we took advantage of the steady advances in computational power of the FNAL UNIX farms [45] to reach the desired level of statistical accuracy.
As inputs to the geant simulation, we employed single particle distributions, reconstructed data events, and two physics event generators: herwig [46] and pythia [47]. For the analyses described herein, we chose herwig as the principal Monte Carlo event generator based on a better match of particle multiplicities between data and Monte Carlo using the default parameters. Over 5.5 million herwig events were passed through the geant simulation. We weighted the herwig π 0 and η spectra (in p T and rapidity) to match our measured results, so that the corrections obtained from the Monte Carlo were based on the data distributions rather than on the behavior of the physics generator.
The calibration of the energy response of the EMLAC was based on the reconstructed masses of π 0 mesons in the γγ decay mode [28]. The steeply falling p T spectrum for π 0 production, combined with the calorimeter's resolution, produced a small offset (≈1%) in the mean reconstructed photon energies. Using the same calibration procedure in the simulated EMLAC as in the detector, we corrected this offset and minimized any potential biases in the calibration. We also employed the geant Monte Carlo simulation to evaluate the mean correction (as a function of photon energy) for energy deposited in the material upstream of the EMLAC. The impact of detector resolution on the energy scale and on the π 0 and η production spectra was incorporated in the overall reconstruction efficiency corrections.
To ensure that the Monte Carlo simulation reproduced the data, a special preprocessor was used to convert geant information into signals and strip energies as measured in the various detectors, and to simulate hardware effects, such as channel noise and gain variations.
The generated Monte Carlo events were then processed through the same reconstruction software used for the analysis of data, and thereby provided measures of the inefficiencies and biases of the reconstruction algorithms.  Figure 11 shows a comparison between the Monte Carlo and the data for the sideband-subtracted π 0 energy asymmetry distribution, and the agreement indicates that the Monte Carlo simulation describes accurately the losses of low-energy photons.
Reconstruction inefficiencies for π 0 and η mesons (which satisfied the A γγ , energy asymmetry, requirement) were relatively small over most of the kinematic range. Figure 12 shows the probability for a π 0 to pass the selection requirements imposed on the Monte Carlo events at 530 GeV/c, as a function of p T , for different rapidity intervals. This probability includes losses due to the reconstruction algorithm, the E FRONT /E TOTAL requirement, and the 10% minimum trigger probability requirement. The inefficiency at forward rapidities and high-p T is attributable to the increased difficulty in separating two photons from π 0 decays (coalescence) in this kinematic region.

F. Normalization
Electronic scalers that counted signals from the beam hodoscope, interaction counters, and beam hole counters were used to determine the number of beam particles incident on the target. Other scalers logged the state of the trigger and of components of the data acquisition system. Information from these scalers was used to determine the number of beam particles that traversed the spectrometer when it was ready to record data. This number was corrected for multiple occupancy in the beam hodoscope (beyond that excluded via the beam1 requirement, and for absorption of beam in the target material.
The normalization of the low p T π 0 cross section was independently verified using events from the prescaled beam and interaction trigger samples. In these samples, the absolute normalization can be obtained just by counting events. For these low p T events, the normalization as determined via the scalers and via event counting techniques agree to better than 3%.
Based upon the good agreement between results from these independent normalization methods, combined with the stability of the cross section results from various sections of the run, an evaluation of the internal consistency of the scalers, and a detailed analysis of the design, implementation and performance of the trigger, the net systematic uncertainty in the overall normalization is ≈8%.

G. Secondary Beam Contamination
The 530 GeV/c secondary beam cross section measurements were corrected for the small admixture of pions and kaons present in the beam, estimated as 2.75% π + and 0.5% K + , respectively [21]. Although the percentage contamination is small, it's effect at high p T is enhanced by the presence of two, rather than three, valence quarks in incident mesons.
The effect of π + contamination was estimated using our high statistics study of π 0 and η meson production by 515 GeV/c π − beam [4,48]. This is justified because: neutral meson production by π + and π − beams is expected to be similar, based on isospin arguments and earlier measurements [49]; the π + component of the nominally 530 GeV/c positive secondary beam had a mean momentum of 515 GeV/c; and the ratio of the measured high p T π 0 cross section from theČerenkov-tagged π + component of the positive secondary beam to the corresponding cross section from the 515 GeV/c π − beam was consistent with unity [21].
The effect of K + contamination was assumed to be half that of the same amount of π + contamination, consistent with earlier, lower energy measurements [50] and with our own, more statistically limited, data [21]. After correcting for beam contamination, the 530 GeV/c cross sections were reduced by ≈ 2% at low p T and by ≈ 10% at high p T .

H. Summary of systematic uncertainties
The principal contributions to the systematic uncertainty arose from the following sources: calibration of photon energy response (5-9%), π 0 and η reconstruction efficiency and detector-resolution unsmearing (5%), the overall normalization (8%), and, for the 530 GeV/c secondary beam, the beam contamination (0-7%). A more complete list of systematic uncertainties is presented in Table I. Some of these uncertainties (e.g. normalization) are strongly correlated between bins. The systematic uncertainties, combined in quadrature, are quoted with cross sections in the appropriate tables. Combining all systematic uncertainties yields a net uncertainty that varies from 11% (13%) at low p T for π 0 (η) mesons to 15% at high-p T .
The secondary proton beam was determined to have a mean momentum of 530 ± 2 GeV/c with an estimated halfwidth of ≈ 30 GeV/c. This momentum spread introduces a small uncertainty (≈ 5%) in comparisons of theory with data. (For the 800 GeV primary beam, the momentum bite and the corresponding uncertainty are negligible.)

IV. INCLUSIVE CROSS SECTIONS
A. π 0 production The inclusive π 0 cross sections per nucleon versus p T are shown in Figs. 13 and 14 for 530 and 800 GeV/c proton beams, respectively, incident upon beryllium and liquid hydrogen targets. Because of the steeply falling spectra, the data are plotted at abscissa values that correspond to the average values of the cross section in each p T bin (assuming local exponential p T dependence) [51]. These cross sections are also tabulated in Tables II   and III These discrepancies have been interpreted [4][5][6] as arising from additional soft-gluon emission in the initial state that is not included in the NLO calculation, and which results in sizeable parton k T before the hard collision (for a different perspective, see the discus-sion in Ref. [8]). Soft-gluon (or k T ) effects are expected in all hard-scattering processes, such as the inclusive production of jets, high-p T mesons, and direct photons [55][56][57][58]. The Collins-Soper-Sterman resummation formalism [59] provides a rigorous basis for understanding these radiation effects, and there have been several recent efforts to derive resummation descriptions for the inclusive direct-photon [60][61][62][63] and dijet cross sections [64][65][66]. The calculation of Ref. [60] for inclusive direct-photon production, which includes the effects of soft-gluon resummation near the kinematic threshold limit (x T = 2p T / √ s −→ 1), has a far smaller sensitivity to scale, compared to NLO calculations, and provides cross sections close to those of NLO calculations with a scale of µ = p T /2. Also, for our energies, the calculation of Ref. [67] using k T -resummation, and of Ref. [63], which treats simultaneously threshold and recoil effects in direct-photon production, yield a substantially larger cross section than the NLO result. However, no such calculations are available for inclusive meson production. In their absence, we use a PQCD-based model that incorporates transverse kinematics of initial-state partons to study the principal consequences of additional k T for high-p T production processes.
Because the unmodified PQCD cross sections fall rapidly with increasing p T , the net effect of the "k T smearing" is to increase the expected yield at higher p T 's. An exact treatment of the modified parton kinematics has been implemented in a Monte Carlo calculation of the leading-order (LO) cross sections for high-p T particle production [68], with the k T distribution for each of the incoming partons represented by a Gaussian with one adjustable parameter ( k T ). Unfortunately, no such program is available for NLO calculations, and so we approximate the effect of k T smearing by multiplying the NLO cross sections by the corresponding LO k T -enhancement factors. Admittedly, this procedure involves a risk of double-counting since some of the k T -enhancement may already be contained in the NLO calculation. However, we expect such double-counting effects to be small. Such comparisons are shown in Fig. 19 for pBe interactions at 530 and 800 GeV/c. The k T -enhanced calculations using the KKP fragmentation functions are seen to reproduce the shape of the π 0 cross section better than calculations using the earlier BKK fragmentation functions, but small discrepancies are still present. Similar conclusions can be drawn from comparisons of NLO QCD with our data on pp collisions (Fig. 20), although in this case the fractional differences using the KKP fragmentation functions are systematically greater than zero when we employ the same k T values used in Fig. 19. It is, however, worth noting that the hydrogen nucleus is not typical, having an up to down quark ratio of two rather than a value of less than one as in the case of standard nuclei, such as Be (Cu) for which The results discussed in this section are not very sensitive to the reasonably well-known parton distribution functions [4] (quark distributions being of primary importance here).
Methods similar to the ones described in this paper have been applied to analyze high-p T hadron spectra in pp, pA, and AA collisions [71,72], and k T effects have been found important for describing data on inclusive production of charged mesons.

B. η production
Cross sections for inclusive η production are tabulated in Tables VIII through XIII. Theoretical descriptions of η-meson production differ from the π 0 case primarily because of differences in the fragmentation of partons into the particles of interest. To investigate this aspect, we present η/π 0 relative production rates as functions of p T and y cm (for two p T ranges) in Figs. 24 and 25 -the data average to a value of 0.45 ± 0.01 at 530 GeV/c, and 0.42 ± 0.01 at 800 GeV/c (for 3 < p T < 8 GeV/c).

V. SUMMARY
The invariant cross sections for π 0 and η production have been measured for pp and pBe collisions at 530 and 800 GeV/c as functions of p T and y cm , over the kinematic range 1 < p T < 10 GeV/c, and 1.5 units in rapidity. Results from k T -enhanced NLO QCD calculations are in reasonable agreement with our measured π 0 cross sections. Employing the recent KKP fragmentation functions in the NLO QCD calculations was found to improve the description of the detailed shape of our π 0 cross sections (as a function of p T ) relative to theoretical results obtained using the earlier BKK fragmentation functions. The measured η/π 0 production ratios, which provide information about the relative fragmentation of partons into these mesons, are 0.45 ± 0.01 at 530 GeV/c and 0.42 ± 0.01 at 800 GeV/c, averaged over p T and rapidity. [26] These hybrid planes were installed after the commissioning run to provide improved vertex resolution.
[ [36] The fiducial volume of the hydrogen target was defined to exclude interactions with the mylar walls of the target flask.
[37] Binned in p T , rapidity, and vertex position, simulated π 0 and η mesons were allowed to decay isotropically to two photons. Photon pairs that satisfied the energy asymmetry requirement (as defined in the text) were projected to the EMLAC, and selected on the basis of the standard fiducial acceptance criteria.
[38] The calculations of invariant mass and p T assumed that the photons originated from the reconstructed interaction vertex. [40] The response of each trigger group was evaluated as a function of p T reconstructed within the group, using data samples that were not biased in that region. The sample used for obtaining the efficiency of the single local hi trigger was based on information from octants that fired the single local lo trigger. The sample used to evaluate the response of the single local lo trigger used information from the octants opposite those that satisfied the single local hi trigger. In the case of the pretrigger hi, octants located opposite those that satisfied the single local hi trigger were used to measure the trigger efficiency for the trigger groups comprised of the inner and outer R channels in each octant.
[44] This is in contrast to the analysis of data from the commissioning run where we used just empirical parameterizations of the electromagnetic showers in simulating the EMLAC response. GeV/c data. The clearly-resolved individual target elements are, from the left: two copper foils, the liquid hydrogen target (a mylar flask enclosed between two thin Be windows), followed by the main beryllium target. The two silicon strip detectors (SSDs) immediately upstream and four downstream of the target can also be easily distinguished.               II. Invariant differential cross sections Ed 3 σ/dp 3 per nucleon for the inclusive reaction pBe → π 0 X at 530 and 800 GeV/c, averaged over the rapidity intervals −0.75 ≤ y cm ≤ 0.75 and −1.0 ≤ y cm ≤ 0.5, respectively. 335.0 ± 6.3 ± 37 811 ± 18 ± 91 5. 10 -5.20 254.0 ± 5.2 ± 28 658 ± 17 ± 74 5. 20 -5.30 197.3 ± 4.8 ± 22 486 ± 16 ± 55 5.30 -5.40