Measurement of Particle Production and Inclusive Differential Cross Sections in Collisions at TeV
Abstract
We report a set of measurements of particle production in inelastic collisions collected with a minimumbias trigger at the Tevatron Collider with the CDF II experiment. The inclusive charged particle transverse momentum differential cross section is measured, with improved precision, over a range about ten times wider than in previous measurements. The former modeling of the spectrum appears to be incompatible with the high particle momenta observed. The dependence of the charged particle transverse momentum on the event particle multiplicity is analyzed to study the various components of hadron interactions. This is one of the observable variables most poorly reproduced by the available Monte Carlo generators. A first measurement of the event transverse energy sum differential cross section is also reported. A comparison with a pythia prediction at the hadron level is performed. The inclusive charged particle differential production cross section is fairly well reproduced only in the transverse momentum range available from previous measurements. At higher momentum the agreement is poor. The transverse energy sum is poorly reproduced over the whole spectrum. The dependence of the charged particle transverse momentum on the particle multiplicity needs the introduction of more sophisticated particle production mechanisms, such as multiple parton interactions, in order to be better explained.
pacs:
13.60.Hb, 13.85.HdCDF Collaboration
I Introduction
In hadron collisions, hard interactions are theoretically well defined and described as collisions of two incoming partons along with softer interactions from the remaining partons. The socalled “minimumbias” (MB) interactions, on the contrary, can only be defined through a description of the experimental apparatus that triggers the collection of the data. Such a trigger is set up so as to collect, with uniform acceptance, events from all possible inelastic interactions. At the energy of the Tevatron, MB data consist largely of the softer inelastic interactions. In this paper, only the inelastic particle production in the central part of the region orthogonal to the beam axis is exploited. The diffractive interactions are neglected. An exhaustive description of inelastic nondiffractive events can only be accomplished by a nonperturbative phenomenological model such as that made available by the pythia Monte Carlo generator.
The understanding of softer physics is interesting not only in its own right, but is also important for precision measurements of hard interactions in which soft effects need to be accounted for. For example, an interesting discussion on how nonperturbative color reconnection effects between the underlying event and the hard scattering partons may affect the top quark mass measurement can be found in (1). Also, effects due to multiple partonparton interactions must be accounted for in MB measurements. A detailed understanding of MB interactions is especially important in very high luminosity environments (such as at the Large Hadron Collider) (2) where a large number of such interactions is expected in the same bunch crossing. MB physics offers a unique ground for studying both the theoretically poorly understood softer phenomena and the interplay between the soft and the hard perturbative interactions.
The observables that are experimentally accessible in the MB final state, namely the particle inclusive distributions and correlations, represent a complicated mixture of different physics effects such that most models could readily be tuned to give an acceptable description of each single observable, but not to describe simultaneously the entire set. The pythia Tune A (3) event generator is, to our knowledge, the first model that comes close to describing a wide range of MB experimental distributions.
In this paper three observables of the final state of antiprotonproton interactions measured with the CDF detector at TeV are presented: 1) the inclusive charged particle transverse momentum () (4) differential cross section, 2) the transverse energy sum () differential cross section, and 3) the dependence of the charged particle average transverse momentum on the charged particle multiplicity, .
The first two measurements address two of the basic features of inelastic inclusive particle production. The measurement of the event transverse energy sum is new to the field, and represents a first attempt at describing the full final state including neutral particles. In this regard, it is complementary to the charged particle measurement in describing the global features of the inelastic cross section. In this article, previous CDF measurements (5) (6) are widely extended in range and precision. The single particle spectrum now extends to over 100 GeV/, and enables verification of the empirical modeling (7) of minimumbias production up to the high production region spanning more than twelve orders of magnitude in cross section. The is one of the variables most sensitive to the combination of the physical effects present in MB collisions, and is also the variable most poorly reproduced by the available Monte Carlo generators. Other soft production mechanisms (8), different from a phenomenological extrapolation of QCD to the nonperturbative region, might show up in the high multiplicity region of . Should this be the case, we might expect to observe finalstate particle correlations similar to those observed in ionion collisions (9).
A comparison with the pythia Monte Carlo generator model (10) is carried out for all the distributions and correlations studied.
The rest of this paper is organized as follows: Sec. II describes the detector components most relevant to this analysis. Section III describes the triggers and the datasets used, including a short description of the Monte Carlo generator tuning, the event selection and the backgrounds. In Sec. IV the methods used to correct the data for detector inefficiency and acceptance are discussed. Section V is devoted to the discussion of the systematic uncertainties. In Sec. VI the results are presented and compared to model predictions.
Ii The CDF detector
CDF II is a general purpose detector that combines precision charged particle tracking with projective geometry calorimeter towers. A detailed description of the detector can be found elsewhere (11). Here we briefly describe the detector components that are relevant to this analysis: the tracking system, the central calorimeters, and the forward luminosity counters.
The tracking system is situated immediately outside the beam pipe and is composed of an inner set of silicon microstrip detectors and an outer drift chamber (COT). The silicon detectors are located between radii of cm, and provide precision measurements of the track’s impact parameter with respect to the primary vertex. The innermost layer (L00) (12) is single sided, and is attached directly on the beam pipe. Five layers of doublesided silicon microstrips (SVXII) (13) cover the pseudorapidity region: in each layer one side is oriented at a stereo angle with respect to the beam axis to provide three dimensional measurements. The ISL (14) is located outside SVXII. It consists of one layer of silicon microstrips covering the region and of two layers at where the COT coverage is incomplete or missing. The COT (15) is a cylindrical opencell drift chamber with 96 sense wire layers grouped into eight alternating superlayers of stereo and axial wires. Its active volume covers cm and cm, thus providing fiducial coverage up to to tracks originating within cm. Outside the COT, a solenoid provides a 1.4 T magnetic field that allows the particle momenta to be computed from the trajectory curvature. The transverse momentum resolution is GeV for the integrated tracking system and GeV for the COT tracking alone.
Located outside the solenoid, two layers of segmented sampling calorimeters (electromagnetic (16) and hadronic (17)) are used to measure the energy of the particles. In the central region, , the calorimeter elements are arranged in a projective tower geometry of granularity . The electromagnetic components use leadscintillator sampling. A multiwire proportional chamber (CES) is embedded at approximately the depth of the shower maximum. The hadron calorimeter uses iron absorbers and scintillators. At normal incidence the total depth corresponds to about 18 radiation lengths in the electromagnetic calorimeter and 4.5 interaction lengths in the hadronic calorimeter.
The energy resolution of the electromagnetic calorimeter is (GeV)) 2% for electromagnetic particles. It is (GeV)) 3% for single pions when using both calorimeters.
Two systems of gas Cherenkov counters (CLC) (18), covering the forward regions , are used to measure the number of inelastic collisions per bunch crossing and to determine the luminosity. For triggering purposes only, this analysis exploits a TimeofFlight detector (TOF) (19) located between the COT and the solenoid at a mean radius of 140 cm. The TOF consists of 216 scintillator bars with photomultipliers at each end and covers roughly .
Iii Data Sample and Event Selection
This analysis is based on an integrated luminosity of 506 pb collected with the CDF II detector between October 2002 and August 2004. The data were collected with a minimumbias trigger that operates as follows. An antiprotonproton bunch crossing, signalled by the Tevatron radio frequency, is defined to contain at least one interaction if there is a coincidence in time of signals in both forward and backward CLC modules. This required coincidence is the start gate of the firstlevel CDF trigger (Level 1) and is the socalled minimumbias trigger. CDF uses a threelevel trigger system that selects events to be recorded to tape at Hz from the bunch crossing rate of approximately 2.5 MHz. The minimumbias trigger is rate limited at Level 1 in order to keep the Level 3 output at 1 Hz. A total of about bunch crossings was recorded.
Part of the analysis also uses data collected with a high multiplicity trigger that selects events that passed the minimumbias trigger precondition and in addition have a large number of primary charged particles. It functions at Level 1 by selecting events with at least 14 hit bars in the TOF system, a hit being defined as the coincidence of two signals from the photomultipliers at the two ends of each bar. At Level 3 this trigger requires at least 22 reconstructed tracks converging to the event vertex. The threshold of 14 TOF signals was selected as the highest compatible with a fully efficient trigger for events with offline charged particle multiplicity . The latter threshold was dictated by the statistics available in Run I and that expected for Run II. This data sample consists of about 64000 triggered events.
For transverse energy measurements, only part of the MB sample was used. Only runs with initial instantaneous luminosity below cms have been kept in order to reduce the effects of event pileup in the calorimeters. The total number of bunch crossings accepted in this subsample is about . The average instantaneous luminosities of the two MB samples are roughly cms for the energy subsample and cms for the full sample.
An offline event selection is applied to the recorded sample of minimumbias triggered events. Events that contain cosmicray candidates, identified by the combination of tracking and calorimeter timing, are rejected. Only those events collected when all the detector components were working correctly are included in the final reduced data sample.
iii.1 Event Selection
Primary vertices are identified by the convergence of reconstructed tracks along the axis. All tracks with hits in at least two COT layers are accepted. No efficiency correction is applied to the tracks used for this task. Vertices are classified in several quality classes: the higher the number of tracks and their reconstruction quality (Sec. IV.1), the higher the class quality assigned to the vertex. For vertices of lowest quality (mainly vertices with one to three tracks) a requirement that they be symmetric is added, i.e. there must be at least one track in both the positive and negative rapidity regions for the vertex to be accepted as primary. In other words, the quantity , where is the number of tracks in the positive or negative hemisphere, cannot equal one.
Events are accepted that contain one, and only one, primary vertex in the fiducial region cm centered around the nominal CDF position. This fiducial interval is further restricted to cm when measurements with the calorimeter are involved.
The event selection described contains an unavoidable contamination due to multiple vertices when the separation between vertices is less than the vertex resolution in the coordinate, which is about 3 cm. A correction for this effect is discussed in Sec. VI.
iii.2 Trigger and Vertex Acceptance
Due to small inefficiencies in the response of the CLC detector, the minimumbias trigger is not 100% efficient. The efficiency has been evaluated by monitoring the trigger with several central high transverse energy triggers, such as those containing a high track, a central high electron, or a central high jet. The results show that the trigger efficiency increases with the increase of some global event variables such as central multiplicity and central sum .
On the other hand, the total acceptance (including the efficiency) of the trigger has been measured by comparing it with a sample of zerobias events collected during the same period. The zerobias data set is collected without any trigger requirements, simply by starting the data acquisition at the Tevatron radiofrequency signal. The results are in agreement with previous studies (20) and indicate that the efficiency depends on a number of variables, most of which in some way are related to the number of tracks present in the detector: number of beam interactions, number of tracks, instantaneous luminosity and the CLC calibration. We parametrized the dependence on these variables so that a correction can be applied on an eventbyevent basis.
The total MB trigger acceptance increases linearly with the instantaneous luminosity. As a function of the number of tracks, the acceptance is well represented by a typical turnon curve starting at about 20% (two tracks) and reaching its plateau with a value between 97 and 99% for about 15 tracks.
As stated above, the present analysis includes data collected with the high multiplicity trigger previously described. The offline selection for these data is the same as that for the minimumbias. Events from the high multiplicity trigger are accepted if they have reconstructed charged track multiplicity at Level 3 greater than or equal to 22. This value is a compromise between the desire for larger statistics in the multiplicity region where the cross section drops and the available trigger bandwidth. The trigger efficiency for this multiplicity is higher than 97%.
The primary vertex recognition efficiency for the MB data sample is evaluated in two ways: by comparing the number of expected vertices on the basis of the instantaneous luminosity and by using a Monte Carlo simulation with multiple interactions. This efficiency was studied as a function of various event variables and found to be roughly flat for 40 cm, but strongly dependent on the number of interactions in the bunch crossing and on the number of tracks available for vertex clustering. Therefore the efficiency has been parametrized as a function of the number of tracks and of the instantaneous luminosity.
Because of their dependence on the number of tracks in the bunch crossing, a variable closely related to the event particle multiplicity, both the trigger and the vertex efficiencies affect not only the total cross section but also the shape of inclusive distributions. The efficiency values are computed on an eventbyevent basis, and are common to all the distributions analyzed.
iii.3 Backgrounds
Diffractive events, with finalstate particles mostly confined in the forward regions, may have some activity in the central region that enters as a background in our sample. By assuming the following indicative values 44.4/10.3/7.0 mb for the centralinelastic, single, and doublediffractive cross sections (21), respectively, and knowing the relative CLC acceptances, we estimate their contribution to the MB cross section to be approximately 6%. Roughly the same conclusion was drawn by analyzing a sample of diffractive events generated with the pythia simulation and passed through a MB trigger simulation. Considering that in about half of the diffractive events no primary vertex is reconstructed, we estimate that diffractive production forms up to of our MB sample and is concentrated in the region of low charged particle multiplicity and low .
For the energy measurements, the presence of calorimeter towers with significant energy deposits not due to particles originating from the interaction was checked. In a sample of zerobias events, after requiring no reconstructed tracks and no signal in the CES, about 0.002 towers per event were found above the pedestal threshold. This number increases with the instantaneous luminosity and is attributed to real particles crossing the calorimeter, probably scattered back from the forward calorimeters. The resulting average energy per event was subtracted from the measurement of each event .
iii.4 The Monte Carlo Sample
A sample of simulated Monte Carlo (MC) events about twice the size of the data was generated with pythia version 6.216 (22), with parameters optimized for the best reproduction of minimumbias interactions. pythia Tune A (3) describes the MB interactions starting from a leading order QCD matrix element augmented by initial and finalstate showers and multiple parton interactions (23), folded in with CTEQ5L parton distribution functions (24) and the Lund string fragmentation model (25). To model the mixture of hard and soft interactions, pythia introduces a cut off parameter (26) that regulates the divergence of the 2to2 partonparton perturbative cross section at low momenta. This parameter is used also to regulate the additional partonparton scatterings that may occur in the same collision. Thus, fixing the amount of multipleparton interactions (i.e., setting the cutoff) allows the hard 2to2 partonparton scattering to be extended all the way down to , without hitting a divergence. The amount of hard scattering in simulated MB events is, therefore, related to the activity of the socalled underlying event in the hard scattering processes. The final state, likewise, is subject to several effects such as the treatments of the beam remnants and color (re)connection effects. The pythia Tune A results presented here are the predictions, not fits.
The MC sample used for all the efficiency and acceptance corrections was generated with Tune A and GeV/. This tuning was found to give a similar output as the default () with only slightly better reproduction of the high particles and a somewhat larger particle multiplicity distribution.
The definition of primary particles was to consider all particles with mean lifetime s produced promptly in the interaction, and the decay products of those with shorter mean lifetimes. With this definition strange hadrons are included among the primary particles, and those that are not reconstructed are corrected for. On the other hand, their decay products (mainly from decays) are excluded, while those from heavier flavor hadrons are included.
A rundependent simulation with a realistic distribution of multiple interactions was employed. Events were fully simulated through the detector and successively reconstructed with the standard CDF reconstruction chain. The simulation includes the CLC detectors used to trigger the MB sample.
Iv Tracking and Energy Corrections
This section describes the procedures adopted to correct the data for detector inefficiencies and limited acceptance, and for reconstruction errors. First, charged particle tracks are selected in such a way as to remove the main sources of background such as secondary particles and misidentified tracks (Sec. IV.1). The tracking efficiency is then computed for the selected tracks, and an appropriate correction is applied to the data distributions (Sec. IV.2). The measurement of requires a careful evaluation of the calorimeter intrinsic response and acceptance, and of other distorting effects, especially in the lower range. A correction for each of these effects is described in Sec. IV.3 and is applied to the data.
iv.1 Track Selection and Acceptance
Reconstructed tracks are accepted if they comply with a minimal set of quality selections including a minimum number of hits, both in axial and stereo layers of the COT. These requirements are made more stringent if no hits in the silicon detectors are used.
All tracks are required to originate in a fiducial region in the plane , where is the nearest distance, projected in the transverse plane, between the track extrapolation and the beam axis; is the distance between the point of closest approach of the track to the axis and the coordinate of the event vertex. The actual region selected in the plane depends on the track itself. Tracks reconstructed including the information from silicon detectors are selected within cm; those reconstructed with no information from the silicon detectors have worse resolution in , and are accepted if cm. A similar selection criterion is used along the beam axis: cm for tracks with silicon information and cm for the remaining tracks. These track selection criteria are used to select primary tracks, and were determined from MC simulation as the ones that maximize the ratio of primary to secondary particles.
As a further requirement, primary charged particles must have a transverse momentum greater than 0.4 GeV/ and pseudorapidity in order to optimize the efficiency and acceptance conditions. The track sample used in this analysis is therefore very different from the one used to reconstruct the event vertex.
The number of primary charged particles in the event after the above selection is defined as the event multiplicity .
iv.2 Tracking Efficiency
The detector acceptance and the tracker efficiency have been analyzed with the aim of estimating a correction to each inclusive distribution presented in the paper. For each track, the multiplicative correction is computed using MC as
(1) 
where is the number of tracks reconstructed as primary and the number of generated primary charged particles. This correction factor includes the track detection and reconstruction efficiency, the correction for the contamination of secondary particles (particle interaction, pair creation), particle decays and misidentified tracks (in MC, reconstructed tracks that do not match to a generated charged particle).
The tracking efficiency is strongly dependent on the number of tracks with a trajectory passing close to the event vertex. To avoid biases due to an incorrect multiplicity distribution in the MC generator, the correction factor was evaluated, as a function of , in ten different ranges of track multiplicity.
The tracking efficiency is the largest contribution to . It is about 70% at GeV/ and increases to about 92% at 5 GeV/, where it reaches a plateau.
The fraction of secondary and misidentified tracks ranges between 1 and 3% over the whole spectrum. The final correction is roughly flat in and , and shows two broad peaks in that correspond to the edges of the silicon detector barrels.
The total correction, as defined in Eq. 1, includes also the smearing correction for very high tracks, where the small curvature may be a source of high dispersion in the reconstructed value, and introduces a significant deviation with respect to the generated .
The measured track distribution is corrected by weighting each track that enters the distribution by the correction (computed at the and values corresponding to that specific track) and by the eventrelated acceptances (trigger and vertex efficiency and diffractive event subtraction described in Sec. III.2 and III.3).
To illustrate the effect of the convolution of all the corrections on the final distribution, the ratio of the fully corrected to the raw distributions is shown in Fig. 1. The correction decreases from 1.6 at GeV/ to 1.05 above 100 GeV/.
The dependence (presented in Sec. VI.2) requires a specific two step correction procedure. First, for each data point at fixed , the correction to the is evaluated and is corrected accordingly. In a second step, a correction is applied for the smearing of the multiplicity of the events. Using MC, a matrix is generated that contains the probability that an event with reconstructed tracks was actually generated with particles:
(2) 
where and refer to the reconstructed and generated multiplicity bin, respectively. In doing this it is assumed that, for all multiplicities, the average of events with generated tracks is the same as that of the events with reconstructed tracks. This is indeed the case after the absolute correction on is applied.
iv.3 Calorimeter Response and Correction of the Distribution
The transverse energy is computed in the limited region as the scalar sum over the calorimeter towers of the transverse energies in the electromagnetic and hadronic compartments:
(3) 
where is the polar angle measured with respect to the direction of the proton beam from the actual primary vertex position. Towers with less than 100 MeV deposition are not included in the sum.
CDF calorimetry is optimized for the measurement of high energy depositions and the analysis of its energy response is not usually performed below a few GeV. In this paper the total distribution is pushed down below this limit and a specific study of the energy correction extension had to be done.
The calorimeter response to single charged particles was checked to be well represented by the simulation down to a track of about 400 MeV/. The simulation of the energy deposition of neutral particles is assumed to be correct. Since the fraction of charged and neutral energy produced in data and in our MC sample agree fairly well, we rely on MC simulation to measure down to GeV the integrated calorimeter response to the total energy deposited.
The list of corrections applied to the data distribution is the following. All corrections are made after the calibration of the calorimeters.

Tower relative correction. The response to the energy entering each calorimeter tower was measured with MC as a function of the of the tower and of the coordinate of the primary vertex and then normalized to the value obtained for the tower with the best response. This correction is introduced to make the calorimeter response flat in and vertex .

Absolute correction for the calorimeter response to the total energy released in each event. This is calculated, using MC, as the ratio of the reconstructed in the calorimeter and corrected for the tower relative response in , to the sum of the transverse energies of the generated primary particles in whose trajectory extrapolates to the same region. The calorimeter response as a function of is shown in Fig. 2.

Correction for the different geometrical acceptance of the calorimeter to events in different positions along the axis as a function of the coordinate of the event vertex. This correction ranges from 1 at to about 0.9 at cm.

Correction for undetected charged particles that curl in the magnetic field and do not reach the calorimeter. The average energy due to low charged particles, estimated from MC, as a function of the event , is added to each event.

Correction for unresolved event pileup. Our rundependent MC sample represents well the average number of multiple interactions. This was checked by plotting the ratios of the distributions at high luminosity to the low luminosity ones. A correction was applied by weighting each event by the ratio of the distribution of the events with only one generated interaction to the distribution of events with only one reconstructed interaction. The correction is done for five different ranges of instantaneous luminosity. This weight ranges from about 0.9 to about 1.1.
In terms of the calorimeter response (Fig. 2), the region below about 5 GeV is the most critical. The reliability of MC in evaluating the calorimeter response was checked – for charged particles – against the single particle response measured from data. A more detailed discussion can be found in Sec. V.4 and leads to a systematic uncertainty as high as 15% on the measurement in this region.
Finally, an unfolding correction for the spread of the events with due to the finite energy resolution is applied. The unfolding is carried out in three steps. (a) An unfolding factor defined as
(4) 
where and indicate respectively the generated and the reconstructed values, is extracted from MC; (b) in order to avoid biases due to the fact that the MC does not perfectly reproduce the data, pythia Tune A is reweighted until it accurately follows the data distribution; (c) a new unfolding factor is computed from the reweighted MC sample and is applied to the corrected data distribution.
The unfolding factor as a function of the event is shown in Fig. 3. The final corrected distribution is therefore obtained as
(5) 
where and refer to the number of events in the corrected and raw distributions respectively. refers to the th correction in the numeration given above.
V Systematic Uncertainties
The selection criteria applied to the dataset, as well as the procedures and the MC generator used to correct for the distortions of the apparatus, efficiency, acceptance limitation, etc. are sources of systematic uncertainties. Each source may affect the final distributions in different ways. A description of potential sources of uncertainty, and the methods used to calculate their contributions to the systematic uncertainties on the final results is presented in the following. Table 1 shows a summary of the systematic uncertainties.
Source/Distribution  ()  event  () 

Luminosity and Trigger  6%  —  6% 
Vertex  0 – 0.6%  0 – 0.5%  0 – 2% 
Diffractive events  0 – 0.5%  0 – 1%  0 – 8% 
MC tuning  1 – 4 %  %  5 – 15% 
Method  1%  —  — 
Lost  —  —  1% 
PileUp  —  —  0 – 3% 
v.1 Integrated Luminosity, Trigger Efficiency
There is an overall global 6% systematic uncertainty on the effective timeintegrated luminosity measurement (27) that is to be added to all the cross section measurements.
Since the trigger uses the same subdetectors as the luminosity measurement, the uncertainty on the trigger efficiency is already included in the systematic uncertainty on the integrated luminosity measurement.
v.2 Vertex Selection and Efficiency
The final cross sections depend on the correction for vertex reconstruction inefficiency that was evaluated with MC. This correction, applied to the MC sample itself, returns a number of reconstructed vertices that differs by from the number of generated ones. The variation on the track distribution from this effect is minor: it has a maximum of at GeV/ and is negligible above 5 GeV/. On the event the variation is about in the multiplicity region between 1 and 5. On the distribution it is larger: from at GeV to a negligible value above 6 GeV.
v.3 Background of Diffractive Events
There are two possible uncertainties on the correction for the contamination of diffractive events: the value of the diffractive cross section with respect to the inelastic nondiffractive one, and the average number of diffractive particles in the COT region. We let the contribution of diffractive events in MB vary from 5 to 7% and the average multiplicity from 1.0 to 1.4 tracks per event. These values are estimates of the contribution of diffractive processes to the inelastic central production. We take as the uncertainty the maximum variation obtained, which is about 30% of the correction itself. The correction piles up in the low multiplicity region. This uncertainty affects the track cross section by at GeV/, the event by less than in the first two multiplicity bins, and the cross section by 8 to 1% in GeV.
v.4 Uncertainties Related to the MC Generator
The Monte Carlo modeling of any of the kinematic distributions of particles always introduces an uncertainty on the corrections when the data distributions are not well reproduced. To evaluate this uncertainty, a second sample of events was simulated with the same Monte Carlo generator but different tuning (tune DW (28)). This tuning, when employed for MB production, yields less energy per event than both data and Tune A.
The track reconstruction has a small, but non zero, inefficiency in any kinematic variable. The difference produced by different pythia configurations on the final corrected distributions is taken as a systematic uncertainty. We find that the corrected track distribution varies by to % and the dependence varies by less than 1%. To avoid biases due to an incorrect multiplicity distribution in the MC generator, the correction was evaluated in different multiplicity bins. We compare the distributions corrected inclusively (integrating over all particle multiplicities) and differentially with respect to the multiplicity, and we find a relative difference of about over the whole spectrum.
Another uncertainty is due to the contamination of secondary particles. To address this effect, our selection (track and ) is varied both in data and MC and the resulting average number of tracks is compared. No significant variations were observed, after correction, on the average multiplicity.
For the energy measurement, the largest uncertainty is due to the simulation of neutral particles, including the detector simulation and the particle generator. There is no way to disentangle these effects, but their combination may be reflected by a different fraction of neutral energy in MC and in data. This, in turn, may affect the global correction since the energy from neutral particles has a higher calorimeter response than the energy from charged particles. The observed difference in neutral fraction from 0.42 to 0.48 (average values) in data, with respect to MC, corresponds to a variation in the calorimeter response to by 2%.
We take the difference between the distributions corrected with different MC tunings as the uncertainty due to the generator. The uncertainty is about 15% at GeV, drops to about 5% at 10 GeV and then remains roughly constant. Note that, at least in part, this uncertainty includes the previous one concerning the simulation of neutral particles.
The uncertainty on the amount of energy per event due to low looping charged particles depends directly on the generator because the region of lower momenta is difficult to compare to data. The two pythia tunings that we employ give a difference of about 1% in over the whole spectrum, which corresponds to about the same uncertainty on the distribution shape.
v.5 Uncertainties Originating from Event PileUp
Finally, there is an uncertainty due to unresolved pileup of events within 3 cm to each other along the beam line. None of the algorithms that we tried was able to separate these overlaps efficiently.
The impact on was estimated by comparing the average multiplicity at different instantaneous luminosities and it was found to be tracks per event, this being the difference in multiplicity between lowest and highest luminosity regions (Fig. 4). For the uncertainty on the total number of particles in the whole MB sample, we take the difference in multiplicity between the lower and the average luminosity: about 0.04 tracks per event, corresponding to of the average raw multiplicity.
The contribution from such events has been taken into account when counting the number of events that enter the cross section calculation (Sec. VI.1), but an uncertainty on the correction remains. It amounts to 0.005 tracks per event, which corresponds to a variation of 0.1% of the total MB cross section.
The impact on the average track is negligible; the maximum variation observed when varying the luminosity is about 0.004 GeV/. The uncertainty on the shape of the distribution is therefore negligible.
The effect on is also negligible. This becomes clearer when considering that since the effect on the is almost zero, any variation could only be due to the reallocation of events along the multiplicity axis. The ratio of two plots from samples of high and low luminosities shows negligible variation.
In the case of the energy measurement, the effect of undetected pileup is much larger and was corrected for (Sec. VI.3), but a small uncertainty still remains on the correction itself due to the uncertainty on the calibration of the MC pileup process.
We may assume that there is no pileup below a given luminosity (e.g., cm) and use this low luminosity sample to compare to our distribution. The ratio of the two is compatible with unity. However, although the pileup probability in the low luminosity sample is small (), it is not negligible. We may then assume an uncertainty proportional to that of the MB inelastic nondiffractive cross section used by the MC generator. By assuming conservatively an uncertainty of the MB inelastic nondiffractive cross section used by the MC generator of 6 mb, we calculate that this is equivalent to a variation in the sample average luminosity of cm, which would be reflected as a of GeV. This, in turn, corresponds to an uncertainty on the distribution of % at GeV and negligible at GeV.
range (GeV/)  

Run 0, 1800 GeV (Eq.8)  1.290.02  8.260.08  –  0.4  10.  102/64 
Run 0, 1800 GeV (Eq.8)  1.290.02  8.260.07  –  0.5  10.  90/62 
Run 0, 1800 GeV (Eq.8)  1.3 fixed  8.280.02  –  0.4  10.  103/65 
Run II, 1960 GeV (Eq.8)  1.2300.004  8.130.01  –  0.4  10.  352/192 
Run II, 1960 GeV (Eq.8)  1.2230.005  8.110.01  –  0.5  10.  258/182 
Run II, 1960 GeV (Eq.9)  1.290.02  8.300.07  4.30.1  0.4  150.  94/233 
Run II, 1960 GeV (Eq.9)  1.360.04  8.470.09  4.640.07  0.5  150.  80/223 
v.6 Total Systematic Uncertainties
All the sources of uncertainty mentioned in Sec.V add up to the total systematic uncertainty that we attribute to each distribution as shown in the relative plots. Those originating from MC are added linearly, and their sum is added in quadrature with the others. Uncertainties arising due to the finite MC statistics used to calculate the corrections are represented in the error bars on the data points; their contribution is about 50%. For the track distribution, the summed systematic uncertainties range between 3 and 6%, for the correlation from negligible values up to 1.5%, and for the distribution from 5% to 25%. These numbers do not include the 6% uncertainty on the integrated luminosity.
It is worth noting that in this paper the measurements of and spectra are pushed down to very low particle energies. CDF II has limited sensitivity in these regions, so that the correction must necessarily rely heavily on simulation.
Vi Results
vi.1 Track Cross Section
The single particle invariant cross section per unit phasespace element is defined as
(6) 
where , , and are the particle energy, momentum, and rapidity, respectively. The charged particle distributions in bins of and have the same shape and mean values. Therefore the cross section factorizes in and and we may write the invariant differential form as
(7) 
where is the raw number of charged particles that is to be corrected for all efficiencies, , and acceptance . is the effective timeintegrated luminosity of the sample.
The accepted region in is calculated from the for each charged track, always assuming the charged pion mass. To obtain a number of tracks per unit rapidity interval, each track is weighted by evaluated at . This procedure introduces a bias that could be avoided only by assigning the correct particle mass to all the reconstructed tracks, which is not possible experimentally. Using MC, it was estimated that this bias is at most 5% at GeV/, and becomes negligible above 5 GeV/. This estimate has in turn an uncertainty that is difficult to estimate due to the lack of measurements of the relative abundance of particles in MB data.
The acceptance takes into account the limited region and the rejection of crossings with event pileup. In the latter case the number of undetected events was estimated indirectly by plotting the average as a function of the instantaneous luminosity (Fig. 4). In this plot, the increase in is due to the increase in number of pileup events. We assume that virtually no pileup is present at a luminosity of cms. The difference with respect to the at the average luminosity of the sample yields the estimated number of events that went unobserved. The final acceptance within of our event selections for this event sample is .
The differential cross section is shown in Fig. 5. The same measurement was discussed in (29) and last published by the CDF collaboration in 1988 (5). For historical reasons, the data published in 1988 were based on the average of positive plus negative tracks, i.e. only half of the total tracks were included, which explains most of the scale factor of about 2 between the two measurements. Besides this, the new measurement shows a cross section about 4% higher than the previous one. At least part of this difference may be explained by the increased centerofmass energy of the collisions from 1800 to 1960 GeV. It should be noted, however, that in 1988 the integrated luminosity was determined indirectly from the UA4 cross section (30) and from the number of events selected. In the region where the 1800 GeV data are available, the distributions have the same shape.
We observe that modeling the particle spectrum with the powerlaw form used in 1988 to fit the distribution (Eq. 8), does not account for the high tail observed in this measurement (Fig. 5). The form in Eq. 8 is merely empirical, and the s of the 1988 data fits were already quite poor. Nevertheless, in the limited region up to GeV/, we obtain, for the present data, a set of fit parameters compatible with those published in 1988 (Table 2).
(8) 
In our measurement, the tail of the distribution is at least three orders of magnitude higher than what could be expected by simply extrapolating to high the function that fits the low region. In order to fit the whole spectrum, we introduced a more sophisticated parametrization (Eq.9):
(9) 
With this new function, we obtain a good (see table 2) but the data are still not well reproduced above about 100 GeV/.
Figure 6 shows the ratio of data over pythia at hadron level. Also in this case, the data show a larger cross section at high starting from about 20 GeV/. The MC generator does not produce any particles at all beyond 50 GeV/.
vi.2 Mean vs Event Multiplicity
The dependence of on multiplicity is computed as the average of all charged particles in events with the same charged multiplicity , as a function of :
(10) 
The rate of change of versus is a measure of the amount of hard versus soft processes contributing to minimumbias collisions; in simulation the rate is sensitive to the modeling of the multipleparton interactions (MPI) (1). The model that currently best reproduces the correlation, pythia Tune A, was tuned to fit the activity in the socalled underlying event in high transverse momentum jet production (31). However, it uses the same cutoff parameter to regulate the divergence of the primary 2to2 partonparton scattering and the number of additional partonparton interactions in the same collision. In addition, in pythia the final state is subject to color (re)connection effects between different parton interactions of the same collision.
The naive expectation from an uncorrelated system of strings decaying to hadrons would be that the should be independent of . However, already at the ISR and at the (32), and more recently at RHIC and at the Tevatron (29) (33), such flat behavior was convincingly ruled out. A study of the dependence of the mean transverse momentum on the charged multiplicity was already performed by CDF in Run I and published in (6). In the analysis presented here an extension to higher multiplicities, well over 40 particles in the central rapidity region, is presented. The precision greatly benefits from the larger statistics obtained with a dedicated trigger (Sec. III). Data from the high multiplicity trigger are included by merging them into the MB sample. Comparison with Run I data (Fig. 8) suggests that there is no faster rise of at the higher multiplicities. Such a rise could have been considered as an indication of a thermodynamic behavior of an expanding initial state of hadronic matter (34).
If only two processes contribute to the MB final state, one soft, and one hard (the hard 2to2 partonparton scattering), then demanding large would preferentially select the hard process and lead to a high . However, we see from Fig. 8 (Tune A, no MPI) that with these two processes alone, the average increases much too rapidly. MPI provide another mechanism for producing large multiplicities that are harder than the beambeam remnants, but not as hard as the primary 2to2 hard scattering. By introducing this mechanism, pythia in the Tune A configuration gives a fairly good description of and, although the data are quantitatively not exactly reproduced, there is great progress over fits to Run I data (6). Note that the systematic uncertainty is always within 2%, a value significantly smaller than the discrepancy with data. pythia Tune A does a better job at describing the data than the ATLAS tune as described in (35). Both include MPI, but with different choices for the color connections (1). In Fig. 8, the ATLAS, no MPI and Tune A distributions do not reach multiplicities greater than about 35 solely due to the limited statistics of the generated samples.
vi.3 Cross Section
The differential cross section is given by
(11) 
where is the timeintegrated luminosity for this subsample of events and is the corresponding corrected number of events. The efficiency includes all trigger and vertex efficiencies and the acceptance takes into account the limited region ( cm for this analysis) and the rejection of crossings with event pileup.
The differential cross section in for is shown in Fig. 9. The raw and corrected event average transverse energies are and GeV, respectively. This measurement, which represents the total inelastic nondiffractive cross section for events of given , is not comparable with previous results since it is the first of its kind at the Tevatron energies.
Figure 10 shows a comparison with the pythia Tune A simulation at hadron level. The simulation does not closely reproduce the data over the whole spectrum. In particular, we observe that the peak of the MC distribution is slightly shifted to higher energies with respect to the data.
Vii Conclusions
Minimumbias collisions are a mixture of hard processes (perturbative QCD) and soft processes (nonperturbative QCD) and, therefore, are very difficult to simulate. They contain soft beambeam remnants, hard QCD 2to2 partonparton scattering, and multiple parton interactions (soft and hard). To simulate such collisions correctly, the appropriate combination of all the processes involved must be known.
This paper provides a set of high precision measurements of the final state in minimumbias interactions and compares them to the best available MC model. The following observations may be made:

The former powerlaw modeling of the particle spectrum is not compatible with the high momentum tail ( GeV/) observed in data. The change of slope confirms that the MB spectrum is modeled by the mixing of soft and hard interactions. This distribution may be seen as an indirect measurement of such compositeness. The continuity of the spectrum and of the dependence, and the absence of threshold effects on such a large scale, indicate that there is no clear separation of hard and soft processes other than an arbitrary experimental choice. The more recent tunings of the pythia MC generator (Tune A) reproduce the inclusive charged particle distribution in data within 10% up to GeV/ but the prediction lies below the data at high . This may mean that the tune does not have exactly the right fraction of hard 2to2 partonparton scattering and, also, that there is more energy from soft processes in the data than predicted.

The cross section represents the first attempt to measure the neutral particle activity in MB at CDF. The MC generator tuned to reproduce charged particle production does not closely reproduce the shape of the distribution. This might be related to the observation that there is an excess of energy in the underlying event in high transverse momentum jet production over the prediction of pythia Tune A.

Among the observables in MB collisions, the dependence of the chargedparticle momentum on the event multiplicity seems to be one of the most sensitive variables to the relative contributions by several components of MB interactions. This correlation is reproduced fairly well only with pythia Tune A: the mechanism of multiple parton interactions (with strong finalstate correlations among them) has been shown to be very useful in order to reproduce high multiplicity final states with the correct particle transverse momenta. In fact, the data very much disfavor models without MPI, and put strong constraints on multipleparton interaction models.
The results presented here can be used to improve QCD Monte Carlo models for minimumbias collisions and further our understanding of multiple parton interactions.
Viii Acknowledgments
We thank the Fermilab staff and the technical staffs of the participating institutions for their vital contributions. This work was supported by the U.S. Department of Energy and National Science Foundation; the Italian Istituto Nazionale di Fisica Nucleare; the Ministry of Education, Culture, Sports, Science and Technology of Japan; the Natural Sciences and Engineering Research Council of Canada; the National Science Council of the Republic of China; the Swiss National Science Foundation; the A.P. Sloan Foundation; the Bundesministerium für Bildung und Forschung, Germany; the Korean Science and Engineering Foundation and the Korean Research Foundation; the Science and Technology Facilities Council and the Royal Society, UK; the Institut National de Physique Nucleaire et Physique des Particules/CNRS; the Russian Foundation for Basic Research; the Ministerio de Ciencia e Innovación, and Programa ConsoliderIngenio 2010, Spain; the Slovak R&D Agency; and the Academy of Finland.
Ix Appendix: Data Tables
(GeV/)  (mb/(GeV/))  stat. err.  (GeV/)  (mb/(GeV/))  stat. err. 

0.40  0.41  1.0145 E+02  6.3 E01  2.45  2.50  1.395 E01  1.5 E03 
0.41  0.42  1.0215 E+02  6.4 E01  2.50  2.55  1.243 E01  1.4 E03 
0.42  0.43  9.685 E+01  6.2 E01  2.55  2.60  1.111 E01  1.3 E03 
0.43  0.44  9.245 E+01  6.0 E01  2.60  2.65  1.004 E01  1.1 E03 
0.44  0.45  8.811 E+01  5.8 E01  2.65  2.70  8.97 E02  1.0 E03 
0.45  0.46  8.403 E+01  5.6 E01  2.70  2.75  8.232 E02  9.8 E04 
0.46  0.47  8.007 E+01  5.4 E01  2.75  2.80  7.325 E02  8.8 E04 
0.47  0.48  7.688 E+01  5.2 E01  2.80  2.85  6.656 E02  8.0 E04 
0.48  0.49  7.360 E+01  5.0 E01  2.85  2.90  5.952 E02  7.4 E04 
0.49  0.50  7.021 E+01  4.8 E01  2.90  2.95  5.390 E02  6.8 E04 
0.50  0.51  6.701 E+01  4.6 E01  2.95  3.00  4.949 E02  6.2 E04 
0.51  0.52  6.404 E+01  4.4 E01  3.00  3.05  4.475 E02  5.7 E04 
0.52  0.53  6.126 E+01  4.3 E01  3.05  3.10  4.070 E02  5.3 E04 
0.53  0.54  5.846 E+01  4.1 E01  3.10  3.15  3.698 E02  5.2 E04 
0.54  0.55  5.563 E+01  3.9 E01  3.15  3.20  3.345 E02  4.6 E04 
0.55  0.56  5.318 E+01  3.8 E01  3.20  3.25  2.994 E02  4.2 E04 
0.56  0.57  5.077 E+01  3.6 E01  3.25  3.30  2.824 E02  4.1 E04 
0.57  0.58  4.851 E+01  3.5 E01  3.30  3.35  2.549 E02  3.7 E04 
0.58  0.59  4.634 E+01  3.3 E01  3.35  3.40  2.349 E02  3.5 E04 
0.59  0.60  4.412 E+01  3.2 E01  3.40  3.45  2.123 E02  3.2 E04 
0.60  0.61  4.233 E+01  3.1 E01  3.45  3.50  1.932 E02  3.0 E04 
0.61  0.62  4.029 E+01  3.0 E01  3.50  3.55  1.808 E02  2.9 E04 
0.62  0.63  3.858 E+01  2.8 E01  3.55  3.60  1.634 E02  2.6 E04 
0.63  0.64  3.681 E+01  2.7 E01  3.60  3.65  1.532 E02  2.5 E04 
0.64  0.65  3.528 E+01  2.6 E01  3.65  3.70  1.402 E02  2.4 E04 
0.65  0.66  3.375 E+01  2.5 E01  3.70  3.75  1.282 E02  2.1 E04 
0.66  0.67  3.228 E+01  2.4 E01  3.75  3.80  1.193 E02  2.1 E04 
0.67  0.68  3.091 E+01  2.3 E01  3.80  3.85  1.092 E02  1.9 E04 
0.68  0.69  2.967 E+01  2.2 E01  3.85  3.90  1.009 E02  1.8 E04 
0.69  0.70  2.829 E+01  2.1 E01  3.90  3.95  9.30 E03  1.7 E04 
0.70  0.71  2.715 E+01  2.0 E01  3.95  4.00  8.53 E03  1.6 E04 
0.71  0.72  2.601 E+01  2.0 E01  4.00  4.05  8.07 E03  1.5 E04 
0.72  0.73  2.499 E+01  1.9 E01  4.05  4.10  7.46 E03  1.5 E04 
0.73  0.74  2.392 E+01  1.8 E01  4.10  4.15  6.72 E03  1.4 E04 
0.74  0.75  2.293 E+01  1.8 E01  4.15  4.20  6.41 E03  1.3 E04 
0.75  0.76  2.204 E+01  1.7 E01  4.20  4.25  5.93 E03  1.2 E04 
0.76  0.77  2.115 E+01  1.6 E01  4.25  4.30  5.39 E03  1.1 E04 
0.77  0.78  2.027 E+01  1.6 E01  4.30  4.35  5.04 E03  1.1 E04 
0.78  0.79  1.943 E+01  1.5 E01  4.35  4.40  4.61 E03  1.0 E04 
0.79  0.80  1.871 E+01  1.5 E01  4.40  4.45  4.353 E03  9.8 E05 
0.80  0.81  1.803 E+01  1.4 E01  4.45  4.50  4.067 E03  9.6 E05 
0.81  0.82  1.727 E+01  1.3 E01  4.50  4.55  3.693 E03  9.2 E05 
0.82  0.83  1.655 E+01  1.3 E01  4.55  4.60  3.522 E03  8.4 E05 
0.83  0.84  1.594 E+01  1.3 E01  4.60  4.65  3.165 E03  8.1 E05 
0.84  0.85  1.533 E+01  1.2 E01  4.65  4.70  3.119 E03  7.8 E05 
0.85  0.86  1.469 E+01  1.2 E01  4.70  4.75  2.919 E03  7.4 E05 
0.86  0.87  1.415 E+01  1.1 E01  4.75  4.80  2.705 E03  7.1 E05 
0.87  0.88  1.361 E+01  1.1 E01  4.80  4.85  2.404 E03  6.5 E05 
0.88  0.89  1.313 E+01  1.0 E01  4.85  4.90  2.314 E03  6.3 E05 
0.89  0.90  1.258 E+01  1.0 E01  4.90  4.95  2.155 E03  6.0 E05 
0.90  0.91  1.212 E+01  1.0 E02  4.95  5.00  2.038 E03  5.8 E05 
0.91  0.92  1.1678 E+01  9.6 E02  5.00  5.20  1.784 E03  3.4 E05 
0.92  0.93  1.1216 E+01  9.4 E02  5.20  5.40  1.339 E03  2.8 E05 
0.93  0.94  1.0829 E+01  9.8 E02  5.40  5.60  1.105 E03  2.3 E05 
0.94  0.95  1.0396 E+01  9.3 E02  5.60  5.80  8.392 E04  1.9 E05 
0.95  0.96  1.0021 E+01  9.1 E02  5.80  6.00  6.59 E04  1.7 E05 
0.96  0.97  9.713 E+00  7.9 E02  6.00  6.20  5.54 E04  1.5 E05 
0.97  0.98  9.325 E+00  7.6 E02  6.20  6.40  4.32 E04  1.2 E05 
0.98  0.99  9.024 E+00  7.4 E02  6.40  6.60  3.58 E04  1.1 E05 
0.99  1.00  8.664 E+00  7.0 E02  6.60  6.80  2.979 E04  9.4 E06 
1.00  1.02  8.227 E+00  6.4 E02  6.80  7.00  2.361 E04  8.4 E06 
1.02  1.04  7.662 E+00  6.0 E02  7.00  7.20  1.999 E04  7.2 E06 
1.04  1.06  7.129 E+00  5.6 E02  7.20  7.40  1.655 E04  6.4 E06 
1.06  1.08  6.635 E+00  5.3 E02  7.40  7.60  1.422 E04  5.7 E06 
1.08  1.10  6.188 E+00  4.9 E02  7.60  7.80  1.276 E04  5.3 E06 
1.10  1.12  5.777 E+00  4.6 E02  7.80  8.00  9.60 E05  4.8 E06 
1.12  1.14  5.404 E+00  4.4 E02  8.00  8.20  9.44 E05  4.4 E06 
1.14  1.16  5.057 E+00  4.0 E02  8.20  8.40  7.05 E05  3.7 E06 
1.16  1.18  4.707 E+00  3.8 E02  8.40  8.60  5.97 E05  3.3 E06 
1.18  1.20  4.412 E+00  3.6 E02  8.60  8.80  5.02 E05  3.0 E06 
1.20  1.22  4.127 E+00  3.4 E02  8.80  9.00  4.69 E05  2.8 E06 
1.22  1.24  3.858 E+00  3.2 E02  9.00  9.20  3.98 E05  2.6 E06 
1.24  1.26  3.614 E+00  3.0 E02  9.20  9.40  3.47 E05  2.4 E06 
1.26  1.28  3.409 E+00  2.8 E02  9.40  9.60  3.23 E05  2.2 E06 
1.28  1.30  3.188 E+00  2.7 E02  9.60  9.80  2.18 E05  1.8 E06 
1.30  1.32  2.985 E+00  2.5 E02  9.80  10.00  2.25 E05  1.8 E06 
1.32  1.34  2.809 E+00  2.4 E02  10.00  10.50  1.89 E05  1.1 E06 
1.34  1.36  2.6298 E+00  2.2 E02  10.50  11.00  1.307 E05  8.6 E07 
1.36  1.38  2.476 E+00  2.1 E02  11.00  11.50  1.085 E05  7.5 E07 
1.38  1.40  2.325 E+00  2.0 E02  11.50  12.00  7.29 E06  5.9 E07 
1.40  1.42  2.192 E+00  1.9 E02  12.00  12.50  6.85 E06  5.7 E07 
1.42  1.44  2.053 E+00  1.8 E02  12.50  13.00  4.56 E06  4.4 E07 
1.44  1.46  1.939 E+00  1.7 E02  13.00  13.50  2.99 E06  3.5 E07 
1.46  1.48  1.822 E+00  1.6 E02  13.50  14.00  2.77 E06  3.3 E07 
1.48  1.50  1.725 E+00  1.5 E02  14.00  14.50  2.39 E06  3.0 E07 
1.50  1.52  1.624 E+00  1.5 E02  14.50  15.00  1.73 E06  2.5 E07 
1.52  1.54  1.536 E+00  1.4 E02  15.00  15.50  1.34 E06  2.1 E07 
1.54  1.56  1.441 E+00  1.3 E02  15.50  16.00  1.20 E06  2.0 E07 
1.56  1.58  1.358 E+00  1.2 E02  16.00  16.50  7.1 E07  1.5 E07 
1.58  1.60  1.287 E+00  1.2 E02  16.50  17.00  1.11 E06  1.8 E07 
1.60  1.62  1.212 E+00  1.1 E02  17.00  17.50  5.9 E07  1.3 E07 
1.62  1.64  1.153 E+00  1.1 E02  17.50  18.00  4.2 E07  1.1 E07 
1.64  1.66  1.084 E+00  1.0 E02  18.00  18.50  4.6 E07  1.1 E07 
1.66  1.68  1.0273 E+00  9.7 E03  18.50  19.00  5.5 E07  1.2 E07 
1.68  1.70  9.741 E01  9.3 E03  19.00  19.50  4.2 E07  1.0 E07 
1.70  1.72  9.176 E01  8.8 E03  19.50  20.00  3.84 E07  9.8 E08 
1.72  1.74  8.649 E01  8.3 E03  20.00  21.00  2.61 E07  5.8 E08 
1.74  1.76  8.238 E01  8.0 E03  21.00  22.00  1.45 E07  4.1 E08 
1.76  1.78  7.822 E01  7.6 E03  22.00  23.00  2.27 E07  5.1 E08 
1.78  1.80  7.389 E01  7.2 E03  23.00  24.00  1.45 E07  3.9 E08 
1.80  1.82  6.992 E01  6.9 E03  24.00  25.00  1.16 E07  3.5 E08 
1.82  1.84  6.612 E01  6.5 E03  25.00  26.00  1.00 E07  3.1 E08 
1.84  1.86  6.290 E01  6.3 E03  26.00  27.00  1.48 E07  3.8 E08 
1.86  1.88  5.963 E01  6.1 E03  27.00  28.00  6.20 E08  2.4 E08 
1.88  1.90  5.642 E01  5.7 E03  28.00  29.00  1.08 E07  3.1 E08 
1.90  1.92  5.382 E01  5.5 E03  29.00  30.00  9.2 E09  8.9 E09 
1.92  1.94  5.081 E01  5.3 E03  30.00  32.00  2.28 E08  9.6 E09 
1.94  1.96  4.864 E01  5.1 E03  32.00  34.00  1.77 E08  8.1 E09 
1.96  1.98  4.631 E01  4.8 E03  34.00  36.00  3.07 E08  1.0 E08 
1.98  2.00  4.358 E01  4.7 E03  36.00  38.00  2.69 E08  9.5 E09 
2.00  2.05  4.021 E01  3.8 E03  38.00  40.00  6.8 E09  4.8 E09 
2.05  2.10  3.533 E01  3.3 E03  40.00  42.00  1.69 E08  7.2 E09 
2.10  2.15  3.125 E01  3.0 E03  42.00  44.00  1.77 E08  6.9 E09 
2.15  2.20 