Measurement and QCD Analysis of Neutral and Charged Current Cross Sections at HERA

The inclusive e(+)p single and double differential cross sections for neutral and charged current processes axe measured with the H1 detector at HERA. The data were taken in 1999 and 2000 at a centre-of-mass energy of roots = 319 GeV and correspond to an integrated luminosity of 65.2 pb(-1). The cross sections are measured in the range of four-momentum transfer squared Q(2) between 100 and 30 000 GeV2 and Bjorken x between 0.0013 and 0.65. The neutral current analysis for the new e(+)p data and the earlier e-p data taken in 1998 and 1999 is extended to small energies of the scattered electron and therefore to higher values of inelasticity y, allowing a determination of the longitudinal structure function F-L at high Q(2) (110 - 700 GeV2). A new measurement of the structure function x (F) over tilde (3) is obtained using the new e(+)p and previously published e p neutral current cross section data at high Q(2). These data together with H1 low Q(2) precision data are further used to perform new next-to-leading order QCD analyses in the framework of the Standard Model to extract flavour separated parton distributions in the proton.


Introduction
In 1992 the HERA accelerator began operation, colliding lepton and proton beams within the H1 and ZEUS experiments. The phase space covered by HERA inclusive deep inelastic scattering (DIS) cross section measurements ranges from small Bjorken x at low Q 2 , the four-momentum transfer squared, to large x at Q 2 values larger than the squared masses of the W and Z gauge bosons. These measurements provide an insight into the partonic structure of matter and the dynamics of strong interactions and test quantum chromodynamics (QCD) over a huge kinematic range.
Both neutral current (NC) interactions, ep → eX via γ or Z 0 exchange, and charged current (CC) interactions, ep → νX via W exchange, can be observed at HERA, yielding complementary information on the QCD and electroweak (EW) parts of the Standard Model. The cross sections are defined in terms of three kinematic variables Q 2 , x and y, where y quantifies the inelasticity of the interaction. The kinematic variables are related via Q 2 = sxy, where s is the ep centre-of-mass energy squared.
Measurements of the NC and CC cross sections in e + p scattering have been made by H1 and ZEUS based on ≃ 40 pb −1 data sets taken between 1994 and 1997 [1,2] with protons of energy 820 GeV and positrons of energy 27.6 GeV, leading to a centre-of-mass energy √ s = 301 GeV.
Here, new e + p NC and CC cross section measurements, based on data taken at √ s = 319 GeV in 1999 and 2000, are presented with improved precision using a luminosity of 65.2 pb −1 .
The increased centre-of-mass energy stems from the change in the proton beam energy from 820 GeV to 920 GeV since 1998. These data sets together provide the most accurate neutral and charged current cross sections measured by H1 at high Q 2 (≥ 100 GeV 2 ) in the first phase of HERA operation (HERA-I).
The NC analysis is extended to higher y up to 0.9 for 100 GeV 2 ≤ Q 2 ≤ 800 GeV 2 . This extension of the kinematic range allows a determination of the longitudinal structure function, F L (x, Q 2 ), to be made at high Q 2 for the first time. This analysis is performed on both the 99 − 00 e + p data and the e − p data, taken in 1998 and 1999 with a luminosity of 16.4 pb −1 at √ s = 319 GeV. The extended high-y e − p analysis and F L extraction complement the inclusive cross section measurements published in [3]. The difference in NC cross sections between e + p and e − p scattering at high Q 2 is employed to update the measurement of xF 3 , superseding the earlier measurement [3].
The accuracy and kinematic coverage of the H1 neutral and charged current cross section data enable dedicated QCD analyses, which test the prediction of logarithmic scaling violations over four orders of magnitude in Q 2 and allow parton distribution functions (PDFs) of the proton to be deduced. This in turn allows predictions to be made for future facilities such as the LHC, deviations from which may be due to exotic phenomena beyond the Standard Model.
A next-to-leading order (NLO) QCD analysis of the H1 data alone is performed, using a novel decomposition of the quark species into the up-and down-type quark distributions to which the NC and CC cross section data are sensitive. The fit parameter space is narrowed using theoretical constraints adapted to the new ansatz and the experimental and phenomenological uncertainties are systematically approached. This leads to a description of the complete set of NC and CC data as well as to new determinations of the PDFs and their uncertainties. For comparison, the QCD analysis is further extended to include the accurate proton and deuteron data from the BCDMS muon scattering experiment [4].
This paper is organised as follows. In section 2 the definitions of the inclusive NC and CC cross sections are given. In section 3 the detector, simulation and measurement procedures are described. The QCD analysis method is explained in section 4, followed by the measurements and the QCD analysis results in section 5. The paper is summarised in section 6.

Neutral and Charged Current Cross Sections
After correction for QED radiative effects, the measured NC cross section for the process e ± p → e ± X with unpolarised beams is given by with where α ≡ α(Q 2 = 0) is the fine structure constant. The ∆ ±,weak N C corrections are defined in [5], with α and the Z and W boson masses (taken here as in [3] to be M Z = 91.187 GeV and M W = 80.41 GeV) as the main electroweak inputs. The weak corrections are typically less than 1% and never more than 3%. The NC structure function term φ ± N C was introduced in [1] and is expressed in terms of the generalised structure functionsF 2 , xF 3 andF L . The helicity dependences of the electroweak interaction are contained in Y ± ≡ 1 ± (1 − y) 2 . The generalised structure functionsF 2 and xF 3 can be further decomposed as [6] In the quark parton model (QPM) the structure functions F 2 , F γZ 2 and F Z 2 are related to the sum of the quark and anti-quark momentum distributions, xq(x, Q 2 ) and xq(x, Q 2 ), [e 2 q , 2e q v q , v 2 q + a 2 q ]{q + q} (5) and the structure functions xF γZ 3 and xF Z 3 to their difference, which determines the valence quark distributions, xq v (x, Q 2 ), In eqs. 5 and 6, e q is the electric charge of quark q and v q and a q are respectively the vector and axial-vector weak coupling constants of the quarks to the Z 0 . In the QPM the longitudinal structure functionF L ≡ 0.
For CC interactions the measured unpolarised ep scattering cross section corrected for QED radiative effects may be expressed as with where ∆ ±,weak CC represents the CC weak radiative corrections. In this analysis G F is defined [8] using the weak boson masses and is in very good agreement with G F determined from the measurement of the muon lifetime [7]. The CC structure function term φ ± CC [1] is expressed in terms of the CC structure functions W ± L , W ± 2 and xW ± 3 , defined in a similar manner to the NC structure functions [8]. In the QPM (where W ± L ≡ 0), they may be interpreted as lepton beamcharge dependent sums and differences of quark and anti-quark distributions and are given for an unpolarised lepton beam by Below the b quark mass threshold, xU, xD, xU and xD are defined respectively as the sum of up-type, of down-type and of their anti-quark-type distributions For the presentation of the subsequent measurements it is convenient to define the NC and CC "reduced cross sections" as 6 3 Experimental Technique

H1 Apparatus and Trigger
The H1 co-ordinate system is defined such that the positive z axis is in the direction of the outgoing proton beam (forward direction). The polar angle θ is then defined with respect to this axis. A full description of the H1 detector can be found in [9][10][11]. The detector components most relevant to this analysis are the Liquid Argon (LAr) calorimeter, which measures the positions and energies of particles over the range 4 • < θ < 154 • , a lead-fibre calorimeter (SPACAL) covering the range 153 • < θ < 177 • , the Plug calorimeter covering the range 0.7 • < θ < 3.3 • and the inner tracking detectors, which measure the angles and momenta of charged particles over the range 7 • < θ < 165 • . In the central region, 25 • < ∼ θ < ∼ 165 • , the central jet chamber (CJC) measures charged track trajectories in the (r, φ) plane and is supplemented by two z drift chambers to improve the θ measurement of reconstructed tracks. The forward tracking detector, θ < ∼ 30 • , is used to determine the vertex position of events when no reconstructed CJC track is found.
The ep luminosity is determined by measuring the QED bremsstrahlung (ep → epγ) event rate by tagging the photon in a photon detector located at z = −103 m. An electron tagger is placed at z = −33 m adjacent to the beam-pipe. It is used to check the luminosity measurement and to provide information on ep → eX events at very low Q 2 (photoproduction) where the electron scatters through a small angle (π − θ < 5 mrad).
NC events are triggered mainly using information from the LAr calorimeter. The calorimeter has a finely segmented geometry allowing the trigger to select localised energy deposits in the electromagnetic section of the calorimeter. For electrons with energy above 11 GeV this is 100% efficient as determined using an independently triggered sample of events. At lower energies the triggers based on LAr information are supplemented by using additional information from the tracking detectors. In 1998 the LAr calorimeter electronics were upgraded in order to trigger scattered electrons with energies as low as 6 GeV, the minimum value considered in this analysis. This gives access to the high y kinematic region. For electron energies of 6 GeV, the overall trigger efficiency is 96% for the e + p data set and 90% for the earlier e − p data set.
The characteristic feature of CC events is a large missing transverse momentum P miss T , which is identified at the trigger level using the LAr calorimeter vector sum of "trigger towers", i.e. groups of trigger regions with a projective geometry pointing to the nominal interaction vertex. At low P miss T , the efficiency is enhanced by making use of an additional trigger requiring calorimeter energy in association with track information from the inner tracking chambers. For the minimum P miss T of 12 GeV considered in the analysis the efficiency is 60%, rising to 90% for P miss

Simulation Programs
Simulated DIS events are used in order to determine acceptance corrections. DIS processes are generated using the DJANGO [12] Monte Carlo (MC) simulation program, which is based on LEPTO [13] for the hard interaction and HERACLES [14] for single photon emission off the lepton line and virtual EW corrections. LEPTO combines O(α s ) matrix elements with higher order QCD effects using the colour dipole model as implemented in ARIADNE [15]. The JETSET program is used to simulate the hadronisation process [16]. In the event generation the DIS cross section is calculated with the PDFs of [17]. The simulated cross section is reweighted using a NLO QCD fit (H1 97 PDF fit) to previous data [1].
The detector response to events produced by the generation programs is simulated in detail using a program based on GEANT [18]. These simulated events are then subjected to the same reconstruction and analysis chain as the data.
The dominant photoproduction background processes are simulated using the PYTHIA [19] generator with leading order PDFs for the proton and photon taken from [20]. Further background from QED-Compton scattering, lepton pair production via two-photon interactions, prompt photon production and heavy gauge boson (W ± , Z 0 ) production are included in the background simulation. Further details are given in [1].

Polar Angle Measurement and Energy Calibration
In the neutral current analysis the polar angle of the scattered electron (θ e ) is determined using the position of its energy deposit (cluster) in the LAr calorimeter, together with the interaction vertex reconstructed with tracks from charged particles in the event. The relative alignment of the calorimeter and tracking chambers is determined using a sample of events with a well measured electron track, using information from both the CJC and the z drift chambers. Minimisation of the spatial discrepancy between the electron track and cluster allows the LAr calorimeter and the inner tracking chambers to be aligned. The residual discrepancy then determines the systematic uncertainty on the measurement of θ e , which varies from 1 mrad for θ e > 135 • to 3 mrad for θ e < 120 • .
The calibration of the electromagnetic part of the LAr calorimeter is performed using the method described in [1]. Briefly, the redundancy of the detector information allows a prediction of the scattered electron energy (E ′ e ) to be made based on the electron beam energy (E e ), the polar angle measurement of the scattered electron and the inclusive polar angle (γ h ) [1] of the hadronic final state. This prediction of the double angle (DA) kinematic reconstruction method [21] is then compared with the measured electromagnetic energy, allowing local calibration factors to be determined in a finely segmented grid in z and φ. The calibration procedure is also performed on the simulated data. The final calibration is obtained by application of a further small correction determined from simulation, which accounts for small biases in the reconstruction of γ h . The calibration is cross checked using independent data samples from QED-Compton scattering and two-photon e + e − pair production processes. The total systematic uncertainty on the absolute electromagnetic energy scale varies from 0.7% in the backward part of the calorimeter to 3% in the forward region, where statistics are limited. The hadronic final state is measured using energy deposits in the LAr and SPACAL calorimeters supplemented by low momentum tracks. Isolated low energy calorimetric deposits are classified as noise and excluded from the analysis. The response of the detector to hadrons is calibrated by requiring transverse momentum conservation between the precisely calibrated scattered electron and the hadronic final state in NC events as described in [1]. The electron transverse momentum is defined as P T,e = p 2 x,e + p 2 y,e . The hadronic transverse momentum is determined from P T,h = ( i p x,i ) 2 + ( i p y,i ) 2 , where the summation is performed over all hadronic final state particles i, assuming particles of zero rest mass.
Detailed studies and cross checks of the hadronic response of the calorimeter using the enlarged data sample have led to an improved understanding of the hadronic energy measurement. The calibration procedure is cross checked by requiring energy-momentum conservation, E − P z ≡ (E ′ e − p z,e ) + (E h − P z,h ) = 2E e , with E h − P z,h being the contribution of all hadronic final state particles i (E i − p z,i ). In addition, the reference scale may be taken from the double angle method prediction rather than from the scattered electron. These studies have allowed the systematic uncertainty on the hadronic scale to be reduced with respect to previous measurements [1,3] in the region 12 GeV < P T,h < 50 GeV and γ h > 15 • . The uncorrelated part (see section 3.6) of the hadronic scale uncertainty is reduced to 1% from 1.7% previously. Fig. 1 demonstrates the quality of the hadronic calibration in the stated region of P T,h and γ h , showing the level of agreement between data and simulation after the calibration procedure. In fig. 1(a) the mean value of the ratio P T,h /P T,e is shown. In fig. 1(b) the inelasticity y h , defined from the hadron reconstruction method [22] as y h = (E h − P z,h )/2E e , is compared with the DA variable y DA . In this analysis, it is the relative difference between data and simulation that is relevant and good agreement is found to within 1%. In addition a 1% correlated uncertainty is considered, accounting for possible remaining biases in the reference scale used for the calibrations.

Neutral Current Measurement Procedure
Events from inelastic ep interactions are required to have a well defined interaction vertex to suppress beam-induced background. High Q 2 neutral current events are selected by requiring a compact and isolated energy deposit in the electromagnetic part of the LAr calorimeter 2 . The scattered electron is identified as the cluster of highest transverse momentum. In the central detector region, θ > 35 • , the cluster has to be associated with a track measured in the inner tracking chambers.
As mentioned earlier, energy-momentum conservation requires E − P z = 2E e . Restricting the measured E − P z to be greater than 35 GeV thus considerably reduces the radiative corrections due to initial state bremsstrahlung, where photons escape undetected in the backward direction. It also suppresses photoproduction background in which the scattered electron is lost in the backward beam-pipe and a hadron fakes the electron signal in the LAr calorimeter. Since the photoproduction background contribution increases with y, the analysis is separated into two distinct regions where different techniques are employed to suppress this background. The nominal analysis is restricted to y e < 0.63 for 90 GeV 2 < Q 2 e < 890 GeV 2 and y e < 0.90 for Q 2 e > 890 GeV 2 . This limits the minimum E ′ e to 11 GeV. The high-y analysis is performed for E ′ e > 6 GeV, 0.63 < y e < 0.90 and 90 GeV 2 < Q 2 e < 890 GeV 2 . Here y e and Q 2 e are reconstructed using the scattered electron energy and angle, the so-called electron method.
The NC kinematics in the nominal analysis are reconstructed using the eΣ method [23], which uses E ′ e , θ e and E h − P z,h and has good resolution and small sensitivity to QED radiative corrections over the accessible phase space. In the high-y analysis the electron method gives the best resolution and is used to define the event kinematics.
The nominal data sample consists of about 185 000 events. The comparison of the data and the simulation is shown in fig. 2 for the scattered electron energy and polar angle spectra and the distribution of E − P z , which are used in the reconstruction of x and Q 2 . All distributions are well described by the simulation, which is normalised to the luminosity of the data. In the nominal analysis the small photoproduction contribution is statistically subtracted using the background simulation.
In the high-y analysis, the photoproduction background plays an increasingly important role, as low energies of the scattered electron are accessed. For this analysis, the calorimeter cluster of the scattered electron is linked to a well measured track having the same charge as the electron beam. This requirement removes a sizeable part of the background, where π 0 → γγ decays give rise to fake scattered electron candidates. The remaining background from photoproduction was estimated from the number of data events in which the detected lepton candidate has opposite charge to the beam lepton. This background is statistically subtracted assuming charge symmetry. The charge symmetry is determined to be 0.99 ± 0.07 by measuring the ratio of wrongly charged fake scattered lepton candidates in e + p and e − p scattering, taking into account the difference in luminosity. The charge symmetry is cross checked using a sample of data events in which the scattered electron is detected in the electron tagger and a systematic uncertainty of 10% on the charge symmetry is assigned. Further details are given in [24,25]. In total about 24 000 e + p events and 5 000 e − p events are selected in the high-y analysis. Figs. 3(a)-(c) show the scattered lepton energy spectrum, the polar angle distribution and the E − P z spectrum for both the e + p and e − p data sets after background subtraction. The simulation, normalised to the luminosity of the data, provides a good description of these distributions. In fig. 3(d) the energy spectra of wrong charge lepton candidates in the data sets are shown. Good agreement is observed when the e + p data are normalised to the luminosity of the e − p data set.

Charged Current Measurement Procedure
The selection of charged current events requires a large missing transverse momentum, P miss T ≡ P T,h ≥ 12 GeV, assumed to be carried by an undetected neutrino. In addition the event must have a well reconstructed vertex as for the NC selection. The kinematic variables y h and Q 2 h are H1 Neutral Current e + p and ep High-y Data . In (c) the E − P z cut is not applied, but is indicated by the dashed vertical line. Shown in (d) is the energy distribution of wrongly charged lepton candidates in background events. In (d) the e + p data have been normalised to the luminosity of the e − p data set. determined using the hadron kinematic reconstruction method [22]. In order to restrict the measurement to a region with good kinematic resolution the events are required to have y h < 0.85. In addition the measurement is confined to the region where the trigger efficiency is 50% by demanding y h > 0.03. The ep background is dominantly due to photoproduction events in which the electron escapes undetected in the backward direction and missing transverse momentum is reconstructed due to fluctuations in the detector response or undetected particles. This background is suppressed using the ratio V ap /V p and the difference in azimuth between P T,h as measured in the main detector and the Plug calorimeter, ∆φ h,Plug [3]. The quantities V p and V ap are respectively the transverse energy flow parallel and anti-parallel to P T,h , the transverse momentum vector of the hadronic final state. The residual ep background is negligible for most of the measured kinematic domain, though it reaches 15% at the lowest Q 2 and the highest y. The simulation is used to estimate this contribution, which is subtracted statistically from the CC data sample with a systematic uncertainty of 30% on the number of subtracted events. The non-ep background is rejected as described in [1] by searching for event topologies typical of cosmic ray and beam-induced background. For further details see [26,27].
After all selection criteria are applied, the final CC data sample contains about 1 500 events. The data and simulation are compared in fig. 4 for the P T,h and E h − P z,h spectra, which are directly used in the reconstruction of the kinematic variables y and Q 2 . In both cases the simulation gives a reasonable description of the data.

Cross Section Measurement and Systematic Uncertainties
For both the NC and CC analyses the selected event samples are corrected for detector acceptance and migrations using the simulation and are converted to bin-centred cross sections. The bins in the (x, Q 2 ) plane are defined as in refs. [1,3], based on the consideration of the detector resolution and event statistics. The bins used in the measurement are required to have values of stability and purity 3 larger than 30%. This restricts the range of the NC measurements to y 0.005. The QED radiative corrections (∆ QED N C(CC) ) are defined in [1] and were calculated using the program HERACLES [14] as implemented in DJANGO [12].
The systematic uncertainties on the cross section measurements are presented in tables 6-13. They are split into bin-to-bin correlated and uncorrelated parts. All the correlated systematic errors are found to be symmetric to a good approximation and are assumed so in the following. The total systematic error is formed by adding the individual errors in quadrature.
The correlated and uncorrelated systematic errors are discussed briefly below (see refs. [24][25][26][27][28] for more details). In addition, there is a global uncertainty of 1.5% and 1.8% on the luminosity measurement for the e + p and e − p data respectively, of which 0.5% is common to both (see section 4.1).
• The total uncertainty on the electron energy scale is 1% if the z impact position of the electron at the calorimeter surface (z imp ) is in the backward part of the detector (z imp < −150 cm), 0.7% in the region −150 cm < z imp < 20 cm, 1.5% for 20 < z imp < 100 cm and 3% in the forward part (z imp > 100 cm). The correlated part of the total uncertainty comes mainly from the possible bias of the calibration method and is estimated to be 0.5% throughout the LAr calorimeter. It results in a correlated systematic error on the NC cross section which is typically below 1%, increasing at low y to ∼ 3% for Q 2 < ∼ 1 000 GeV 2 and ∼ 8% for larger Q 2 .
• The correlated uncertainty on the electron polar angle is 1 mrad, 2 mrad and 3 mrad for θ e > 135 • , 120 • < θ e < 135 • and θ e < 120 • , respectively. This leads to a typical uncertainty on the NC reduced cross section of less than 1%, increasing up to ∼ 5% at high x.
• A 0.5% (1%) uncorrelated error originates from the electron identification efficiency in the NC nominal (high-y) analysis for z imp < −5 cm. For z imp > −5 cm the uncertainty is increased to 2%. The precision of this efficiency is estimated using an independent track based electron identification algorithm, limited for z imp > −5 cm by statistics.
• A 0.5% uncorrelated error is assigned on the efficiency of the scattered electron trackcluster link requirement in the NC nominal analysis. In the high-y analysis this is increased to 1%.
• An uncorrelated 1% uncertainty on the hadronic energy measured in the LAr calorimeter is assigned for the region 12 GeV < P T,h < 50 GeV and γ h > 15 • . Outside this region the uncertainty is increased to 1.7%. In addition, a 1% correlated component to the uncertainty is added in quadrature, originating from the calibration method and from the uncertainty on the reference scale. This yields a total uncertainty of 1.4% and 2% for the two regions respectively.
• The uncertainty on the hadronic energy scale of the SPACAL calorimeter is 5%. The uncertainty on the hadronic energy measurement due to the inclusion of low momentum tracks is obtained by shifting their contribution by 3%. The influence on the cross section measurements from these sources is small compared with that from the correlated uncertainty from the LAr calorimeter energy scale. The three contributions (LAr, SPACAL, tracks) are thus combined, resulting in a single correlated hadronic error from the hadronic energy measurement, which is given in the tables. The corresponding error on the NC and CC cross sections is typically < ∼ 1%, but increases at low y to ∼ 5%.
• A 25% uncertainty is assigned on the amount of energy in the LAr calorimeter attributed to noise, which gives rise to a sizeable correlated systematic error at low y, reaching ≃ 10% at x = 0.65 and Q 2 < ∼ 2 000 GeV 2 in the NC measurements.
• In the CC analysis the correlated uncertainties due to the cuts against photoproduction on V ap /V p and ∆φ h,Plug are only significant at high y, low Q 2 and low P T,h , reaching a maximum of ∼ 7%.
• In the CC and the NC nominal analyses the photoproduction background is estimated from simulation. A 30% correlated uncertainty on the subtracted photoproduction background is determined from a comparison of data and simulation for a phase space region dominated by photoproduction background. This results in a systematic error of typically < ∼ 1%.
• In the NC high-y analysis the photoproduction background is estimated directly from the data by using wrongly charged (fake) scattered lepton candidates, which leads to a 10% correlated uncertainty on the subtracted photoproduction background. The resulting uncertainty on the measured cross sections is 1% or less.
• A 0.3% uncorrelated error is considered on the trigger efficiency in the NC nominal analysis and 2 − 6% in the CC analysis. For the NC high-y analysis the uncertainty on the cross section is ∼ 2% at low Q 2 , decreasing to 0.6% at the highest Q 2 in the analysis.
• An uncorrelated error of 1% (NC) and 3% (CC) is estimated on the QED radiative corrections by comparing the radiative corrections predicted by the Monte Carlo program (DJANGO) with those calculated from HECTOR and EPRC [29]. The error also accounts for a small missing correction in DJANGO due to the exchange of two or more photons between the electron and the quark lines.
• A 3% uncorrelated error is assigned on the event losses due to the non-ep background finders in the CC analysis, estimated from pseudo-CC data (see section 3.1).
• A 2% uncorrelated error (5% for y < 0.1) on the vertex finding efficiency for CC events is estimated using pseudo-CC data.
Overall a typical total systematic error of about 3% (6%) is reached for the NC (CC) double differential cross section. This precision has been achieved through detector improvements for triggering and a better understanding of the hadronic response of the detector, the electron identification and its angular measurement.

QCD Analysis
The cross section data presented here, together with the low Q 2 precision data [30] and high Q 2 e ± p data [1,3] previously published by the H1 Collaboration, cover a huge range in Q 2 and x. The improved accuracy now available allows the predictions of perturbative QCD to be tested over four orders of magnitude in Q 2 from about 1 GeV 2 to above 10 4 GeV 2 , and x from below 10 −4 to 0.65. The measurements of NC and CC e ± p scattering cross sections provide complementary sensitivity to different quark distributions and the gluon distribution, xg(x, Q 2 ). This is used to determine the sum of up-type xU, of down-type xD and of their anti-quark-type xU and xD distributions, employing only H1 inclusive cross section data.
With the current beam energies, the HERA collider data do not give access to the large x region of deep inelastic scattering at medium Q 2 (∼ 100 GeV 2 ). Complementary information on quark distributions in this region is provided by fixed-target lepton-proton data. Lepton-deuteron scattering data, which provide further constraints on the PDFs, are not yet available from HERA. Therefore, in the subsequent analysis, the H1 data are also combined with the precise BCDMS muon-proton and muon-deuteron scattering data and the results are compared for cross checks of the PDFs obtained from the analysis of the H1 data alone.

Ansatz
Traditionally, QCD analyses of inclusive deep inelastic scattering cross section data have used parameterisations of the valence quark distributions and of a sea quark distribution, imposing additional assumptions on the flavour decomposition of the sea [1,[31][32][33]. The neutral and charged current cross section data presented here, however, are sensitive to four combinations of up-and down-type (anti-)quark distributions which, for Q 2 less than the bottom quark production threshold, are given above in eq. 10. Working in terms of these combinations weakens the influence of necessary assumptions on the flavour decomposition of the sea in the fit. The valence quark distributions are obtained from and are not fitted directly.
In the QPM, the charged current structure function terms φ ± CC are superpositions of the distributions given in eq. 10 according to The neutral current structure function terms φ ± N C are dominated by the electromagnetic structure function F 2 , which can be written as In the high Q 2 neutral current data, complementary sensitivity is obtained from the interference structure function xF γZ ]/3, but still higher luminosity is required to exploit this for a dedicated determination of the valence quarks.
In the fit to the H1 and BCDMS data, the isoscalar nucleon structure function F N 2 is determined by the singlet combination of parton distributions and a small contribution from the difference of strange and charm quark distributions, The nucleon data obtained from the BCDMS muon-deuteron cross sections are measured for x ≥ 0.07. For these data nuclear corrections are applied following [34]. In eq. 15 the charm and strange quark distributions occur explicitly and may be constrained using experimental data as provided by H1 and ZEUS on the charm contribution to F 2 [35,36] and from NuTeV on the strangeness content of the nucleon [37]. The analysis of the H1 data, however, is rather insensitive to these quark distributions. They are assumed to be fixed fractions of the up-and down-type quark distributions respectively at the initial scale of the QCD evolution (see section 4.2).
The analysis is performed in the MS renormalisation scheme using the DGLAP evolution equations [38] at NLO [39]. The structure function formulae given here are thus replaced by integral convolutions of coefficient functions and PDFs. An approach is used whereby all quarks are taken to be massless, including the charm and bottom quarks, which provides an adequate description of the parton distributions in the high Q 2 kinematic range of the new data presented here. The bottom quark distribution, xb, is assumed to be zero for Fits are performed to the measured cross sections calculating the longitudinal structure functions to order α 2 s and assuming the strong coupling constant to be equal to α s (M 2 Z ) = 0.1185 [7]. All terms in eqs. 1 and 7 are calculated, including the weak corrections, ∆ ±,weak N C,CC . The analysis uses an x space program developed inside the H1 collaboration [40], with cross checks performed using an independent program [41]. In the fit procedure, a χ 2 function is minimised which is defined in [30]. The minimisation takes into account correlations of data points caused by systematic uncertainties allowing the error parameters (see table 2), including the relative normalisation of the various data sets, to be determined by the fit. The fit to only H1 inclusive cross section data, termed H1 PDF 2000, uses the data sets as specified in table 1. The table additionally lists the BCDMS data used in a further fit for comparison with the H1 PDF 2000 fit.  Table 1: Data sets from H1 used in the H1 PDF 2000 fit and from BCDMS µ-proton and µdeuteron scattering used in the H1+BCDMS fit. As for the previous H1 QCD analysis [30], the original BCDMS data are used at four different beam energies imposing the constraint y µ > 0.3. The inelasticity y µ was defined using BCDMS beam energies. The normalisation uncertainties of each data set (δ L ) are given as well as the kinematic ranges in x and Q 2 . The uncertainty δ L includes a common error of 0.5% for the H1 data sets (see text). The nominal analysis and high-y analysis do not overlap in kinematic coverage (see section 3.4).
The correlated systematic uncertainties for the H1 cross section measurements may be correlated across data sets as well as between data points, since they may arise from the same source. They are thus not treated independently in the QCD analysis presented here. The relationship between the error sources as used in the fitting procedure is summarised in table 2 for each of the eight correlated systematics considered. This leads to 18 independent error parameters. In addition, all H1 quoted luminosity uncertainties have a common contribution of 0.5% arising from the theoretical uncertainty on the Bethe-Heitler cross section. This common contribution has been taken into account in the QCD analysis. Table 2: Treatment of the correlated systematic error sources for the H1 data sets used in the fits. For each of the eight correlated systematic error sources, one or more parameters are included in the QCD fit procedure. The sources considered are due to the luminosity uncertainty (δ L ), the electron energy uncertainty (δ E ), the electron polar angle measurement uncertainty (δ θ ), the hadronic energy uncertainty (δ h ), the uncertainty due to noise subtraction (δ N ), the photoproduction dominated background simulation error (δ B ), the uncertainty due to the cuts against photoproduction in the CC analysis (δ V ) and the error on the charge symmetry in the high-y analysis (δ S ). The table entries indicate the correlation of the error sources across the H1 data sets. For example, the uncertainty due to the noise subtraction is the same for all data sets leading to one common parameter in the fit (N1), whereas the electron energy uncertainty has two independently varying parameters (E1 and E2) for the H1 NC data sets only.

Parameterisations
The initial parton distributions, xP = xg, xU, xD, xU , xD, are parameterised at Q 2 = Q 2 0 in the following general form The QCD analysis requires choices to be made for the initial scale (Q 2 0 ) and the minimum Q 2 of the data considered in the analysis (Q 2 min ). Variations of both Q 2 0 and Q 2 min are studied. As in [30]  The general ansatz, eq. 16, represents an over-parameterisation of the data considered here. The specific choice of these parameterisations is obtained from saturation of the χ 2 : an additional parameter D, E, F or G is considered only when its introduction significantly improves the χ 2 . The appropriate number of parameters also depends on the data sets included in the fit. The H1 data requires less parameters than the combined H1 and BCDMS data due to the precise BCDMS proton and deuteron data in the large x region, where the cross section variations with x are particularly strong.
The fit to the H1 data uses the following functional forms in which the number of free parameters are further reduced using the constraints and assumptions detailed below.
The number of parameters required by the fit for the different parton distributions follows the expectation. A high x term E g x 2 is not needed in the gluon distribution, since at large x the scaling violations are due to gluon bremsstrahlung, i.e. are independent of the gluon distribution. The xU and xD distributions require more parameters than the anti-quark distributions xU and xD because the former are a superposition of valence and sea quarks, in contrast to the latter. Due to the different electric charges, e 2 u = 4e 2 d , and the y dependence of the charged current cross section, the data are much more sensitive to the up quark than to the down quark distributions. Thus less parameters are needed for xD than for xU.
A number of relations between parameters can be introduced naturally in this ansatz. At low x the valence quark distributions are expected to vanish and the sea quark and the anti-quark distributions can be assumed to be equal. Thus the low x parameters A q and B q are required to be the same for xU, xU and for xD, xD. In the absence of deuteron data from HERA there is no distinction possible of the rise towards low x between xU and xD. Thus the corresponding B parameters are required to be equal, i.e.
Further constraints are the conventional momentum sum rule and the valence quark counting rules.
The ansatz presented above allows the quark distributions xU, xD, xU , xD to be determined. Further disentangling the individual quark flavour contributions to the sea is possible only with additional experimental information and/or assumptions. Assuming that the strange and charm sea quark distributions xs and xc can be expressed as x-independent fractions f s and f c of xD and xU at the starting scale of Q 2 0 = 4 GeV 2 (see table 5), a further constraint is used in the fit: The total number of free parameters of the five parton distributions is thus equal to 10 in the fit to the H1 data. The χ 2 value is hardly improved by including any half integer power of x. The parametric form of eq. 17 is also found starting from an alternative polynomial in x k , which includes half integer powers up to x 5/2 . The addition of the large x BCDMS µp and µD data leads to two additional terms, G U x 4 and F D x 3 , in the polynomials.

NC and CC Cross Sections dσ/dQ 2 , dσ/dx and σ tot CC
The e + p single differential neutral current cross section dσ/dQ 2 measured for y < 0.9 is shown in fig. 5(a). The data are compared with previous H1 e + p measurements made at √ s = 301 GeV. The new cross sections are found to be higher than the measurement from 94 − 97 as expected due to the increase in centre-of-mass energy. Both cross sections, falling by over six orders of magnitude for the measured Q 2 region between 200 GeV 2 and 30 000 GeV 2 , are well described by the H1 PDF 2000 fit. The error band represents the total uncertainty as derived from the QCD analysis by adding in quadrature the experimental and model uncertainty. The experimental uncertainty on the predicted cross sections is significantly larger than the model uncertainty, which is discussed in section 5.3. Fig. 5(b) shows the ratios of the measurements to the corresponding Standard Model expectation determined from the H1 PDF 2000 fit. Note that in this lower figure the H1 data are scaled by the normalisation shift imposed by the QCD fit given in table 3. The new data are given in table 6.
The Q 2 dependence of the charged current cross section from the 99 − 00 data is shown in fig. 6(a). For consistency with the NC cross sections, the data are presented in the range y < 0.9, after correction 4 for the kinematic cuts 0.03 < y < 0.85 and P T,h > 12 GeV (section 3.5). The data are compared with the previous measurement taken at lower centre-of-mass energy. The ratios of data to expectations are shown in fig. 6(b) together with the Standard Model uncertainty. Again in this lower figure the H1 data are scaled by the normalisation shift imposed by the QCD fit, given in table 3. The two data sets agree well with each other, though the new data have a tendency to be higher than the fit result at high Q 2 . The data are listed in table 7. Fig. 7 shows the Q 2 dependences of the NC and CC cross sections representing the total e + p and e − p data sets taken at HERA-I. The e + p data have been combined after scaling the 94 − 97 data to √ s = 319 GeV, using the H1 PDF 2000 fit and the procedure described in [27]. At low Q 2 the NC cross section exceeds the CC cross section by more than two orders of magnitude. The sharp increase of the NC cross section with decreasing Q 2 is due to the dominating photon exchange cross section with the propagator term ∝ 1/Q 4 . In contrast the CC cross section 2 ) approaches a constant at low Q 2 . The CC and NC cross sections are of comparable size at Q 2 ∼ 10 4 GeV 2 , where the photon and Z 0 exchange contributions to the NC cross sections are of similar size to those of W ± exchange. These measurements thus illustrate unification of the electromagnetic and the weak interactions in deep inelastic scattering. While the difference in NC cross sections between e + p and e − p scattering is due to γZ interference, the difference of CC cross sections arises from the difference between the up and down quark distributions and the less favourable helicity factor in the e + p cross section (see eq. 13).
The single differential cross sections dσ/dx are measured for Q 2 > 1 000 GeV 2 for both NC and CC and also for Q 2 > 10 000 GeV 2 in the NC case. The NC data are compared in fig. 8 with the previous H1 e + p measurement at √ s = 301 GeV and the corresponding expectations from the fit. A similar comparison for the CC data is shown in fig. 9. Increases with √ s are observed in both the NC and the CC cross sections, in agreement with the expectations. The fall in the cross sections at low x is due to the restriction y < 0.9. The measurements are summarised in tables 8-10.
The total CC cross section has been measured in the region Q 2 > 1 000 GeV 2 and y < 0.9 after applying a small correction factor of 1.03 for the y and P T,h cuts, determined from the H1 PDF  where the 1.5% normalisation uncertainty is included in the systematic error. This is to be compared with the value from the H1 PDF 2000 fit σ tot CC (e + p) = 16.76±0.32 pb. The difference between the measurement and the fit is 2.3 standard deviations assuming the correlation of uncertainties between measurement and fit is negligible. An unbiased theoretical expectation for σ tot CC (e + p) may be obtained by repeating the H1 PDF 2000 fit but excluding the new 99 − 00 CC data, which leads to 16.66 ± 0.54 pb.
Additionally, the analysis has been performed on the 94−97 data set at the lower centre-of-mass energy, yielding

NC and CC Double Differential Cross Sections
The double differential NC reduced cross section,σ N C (defined in eq. 11), is shown in fig. 10 for both the nominal and high-y 99 − 00 e + p data. In addition the new high-y 98 − 99 e − p data are presented. The data agree well with the expectations of the H1 PDF 2000 fit, which are also shown 5 . The rise of the Standard Model DIS cross section towards low x (high y) departs from the monotonic behaviour of F 2 due to the contribution of the longitudinal structure function F L . This allows F L to be determined in the high y region (section 5.5).
In fig. 11 the e + p NC large x cross section data at √ s = 319 GeV are compared with the data obtained previously [1] at √ s = 301 GeV. The two data sets are found to be in agreement with each other and with the H1 PDF 2000 fit. Fig. 11 also shows the data from the recent H1 measurement at lower Q 2 [30] and the fixed-target data from BCDMS [4]. The fit description of the BCDMS data, which are not used in the fit, is remarkably good except at very large x = 0.65. A similar observation has already been reported in [1,30]. At the highest Q 2 a decrease of the cross section is expected due to the negative γZ interference in e + p scattering.
In fig. 12 the reduced CC cross section,σ CC (defined in eq. 11), is shown for the new data and the data taken at lower energy between 1994 and 1997. These data are found to be compatible with each other considering the weak energy dependence of the reduced CC cross section. An extension of the x range for Q 2 = 3 000 GeV 2 and 5 000 GeV 2 is achieved due to the improved trigger efficiency. The combined 94 − 00 result is compared in fig. 13 with the expectation from the H1 PDF 2000 fit. Also shown is the expected contribution of the xD distribution, which dominates the cross section at large x. The HERA e + p CC data can thus be used to constrain the d quark distribution in the valence region.
All double differential measurements together with the contributions of each of the major systematic uncertainties are listed in tables 11-13.

Fit Results
In this section, the results of the QCD analysis are presented. The χ 2 value for each data set is given in table 3 as well as the optimised relative normalisation as determined from the fit. The total χ 2 value 6 per degree of freedom (χ 2 /ndf) is 540/(621 − 10) = 0.88. The NLO QCD fit requires the lowest Q 2 data (H1 minimum bias 97, Q 2 ≤ 12 GeV 2 ) to be raised by 3.7%, corresponding to 2.3 standard deviations in terms of the normalisation uncertainty excluding the common error of 0.5% (see section 4.1). In contrast all data for Q 2 100 GeV 2 are lowered, by at most 1.9%. It can not yet be decided whether this behaviour is due to inadequacies in the theory (e.g. the missing higher order terms in ln Q 2 ) or experimental effects.   2) leads to an additional contribution [30], which is listed as χ 2 (corr. err.). Also shown is the optimised normalisation of the data set as determined by the fit. The H1 NC 98 − 99 e − p and H1 NC 99 − 00 e + p data include the high-y analyses.
The parameters of the initial parton distributions are given in table 4 (see also [42]) and the distributions are shown in fig. 14. The inner error band describes the experimental uncertainty, while the outer band represents the experimental and model uncertainties added in quadrature.
The experimental accuracy of the initial distributions is typically a few percent in the bulk of the phase space of the H1 data. This accuracy has negligible dependence on Q 2 but a strong dependence on x. The best precision is achieved for the xU quark distribution, which amounts to about 1% for x = 0.01 and reaches 3% and 7% at x = 0.4 and 0.65, respectively. The xD quark distribution is only determined with moderate precision as it is predominantly constrained by the CC e + p cross sections, which are still subject to limited precision. The corresponding uncertainties on xD at the three quoted x values are respectively ∼ 2%, ∼ 10% and ∼ 30%.
These uncertainties reflect the kinematic dependence and size of the measurement errors. However the error size also depends significantly on the fit assumptions. If, for example, the constraint between A U and A D on the low x behaviour of the anti-quark distributions is relaxed, the small uncertainty at low x = 0.01 is much increased to ∼ 6% and ∼ 20% respectively for xU and xD. The measurement of the low x behaviour of up and down quarks and their possible distinction requires electron-deuteron data to be taken at HERA.
The model parameter uncertainties on the PDFs are determined in a similar manner to [30] and the sources of uncertainty are specified in table 5. The model uncertainties are relatively small with respect to those from experimental sources except at small x and low Q 2 , where they reach ∼ 2% and ∼ 6% respectively for xU and xD at x = 0.01 and Q 2 = 10 GeV 2 .
Within the functional form considered (see eq. 16), the parameterisation given in eq. 17 is found uniquely. Possible variations within the ∆χ 2 ≃ 1 region of the parameter space do not lead to noticeably different distributions. Thus in this analysis no account is made of uncertainties due to the choice of parameterisations. A completely different ansatz, however, may well lead to different initial distributions, as seen, for example, in the complicated shape of xg chosen in [33]. The gluon distribution determined in this analysis is consistent with the distribution obtained previously by H1 [30] if the effects of the different heavy flavour treatments are taken into account.  The full curve in fig. 14 is the result [43] of the fit to H1 and BCDMS data, which gives a χ 2 /ndf= 883/(1014 − 12) = 0.88. Excellent agreement of the PDFs between the two fits is observed. For large x, the high Q 2 data of H1 allow distinction between up and down flavours yielding results compatible with those from BCDMS proton and deuteron data. At low x only HERA data are available and thus the two fits are forced to be in agreement, apart from possible small effects due to sum rules.

23
The PDFs from the H1 PDF 2000 fit are further compared in fig. 15 with recent results from the MRST [32] and CTEQ [33] groups for two values of Q 2 at 10 GeV 2 and 1 000 GeV 2 . The H1 PDF 2000 fit is in remarkable agreement with the MRST and in particular the CTEQ analyses, given the many differences in terms of the data sets used, the assumptions made and the treatment of heavy flavours.

Extraction of the Proton Structure Function F 2
The NC structure function term φ N C is obtained from the measured NC double differential cross section according to eq. 1. It is dominated by the structure function F 2 in most of the kinematic range. The structure function F 2 is extracted using Here the correction terms ∆ F 2 and ∆ F 3 account for the effects of Z 0 exchange onF 2 and xF 3 (eqs. 2-4) and ∆ F L originates from the longitudinal structure functionF L . These corrections, shown in table 11, are determined using the H1 PDF 2000 fit (see section 4). At high y and Q 2 < 1 500 GeV 2 , ∆ F L is sizeable and the extraction of F 2 in this Q 2 region is thus restricted to the kinematic range y < 0.6. It is extended to higher y at larger Q 2 (≥ 2 000 GeV 2 ) where the predicted contribution ofF L is small.
The extracted F 2 using the 99 − 00 data is presented in table 11. Fig. 16 shows the F 2 data using the combined 94 − 97 and 99 − 00 high Q 2 e + p data sets. Also shown in the figure are the recent H1 F 2 data at lower Q 2 [30] and structure function data from BCDMS [4] and NMC [44]. The full H1 data cover a range of four orders of magnitude in x and Q 2 . The H1 PDF 2000 fit provides a good description of the data over the whole region except for the BCDMS data at x = 0.65, as seen in fig. 11. The fit also gives a good description of the scaling violations observed in the measurements.

Determination of the Longitudinal Structure Function F L
The structure function term φ N C is used to determine F L at y > 0.63 and Q 2 < 890 GeV 2 . For statistical reasons, the measured cross sections in two neighbouring Q 2 bins are combined, assuming that the systematic uncertainties are fully correlated. The longitudinal structure function is then determined using the formulã for e ± p scattering which, neglecting the small electroweak contributions in the region of this extraction, reduces to the expression The extraction of F L relies upon the extrapolation of the fit result for F 2 into the high y region, that is, to larger Q 2 for given x. In order to avoid a possible influence of the high y data region on this calculation, a dedicated NLO QCD fit (H1 Low y fit) is performed to H1 data with y < 0.35 only and the results are extrapolated using the DGLAP evolution equations. This method was introduced in [45].
Apart from the y range restriction, the H1 Low y fit follows the same procedure as described in section 4. It results in a χ 2 /ndf= 417/(455 − 10) = 0.94 and agrees very well with the H1 PDF 2000 fit over the full y range. The normalisation shifts of the data sets used are found to be within 1% of those from the H1 PDF 2000 fit.
In the extraction of the longitudinal structure function, the experimental cross sections are slightly modified using the results of the H1 Low y fit for the renormalisation and small shifts from the correlated uncertainties common to the low y and the high y region. The combined HERA-I measurements of the structure function term φ N C and the extracted values of F L are listed in table 14. The statistical precision is due directly to the cross section measurement at high y. The systematic uncertainties arise from the measurement errors at high y and the model uncertainties related to the extrapolation of F 2 from the low y to the high y region. The correlations in the systematic uncertainties between low and high y are taken into account.
In fig. 17 the determinations of F L at high Q 2 are shown for both the e + p and the e − p data sets. The results from both data sets are mutually consistent and are in agreement with the H1 Low y fit prediction for F L , based on the gluon distribution derived from the scaling violations of F 2 at larger x. The extreme values allowed for F L (F L = 0 and F L = F 2 ) are clearly excluded by the data. A model independent measurement of F L and the x dependence can, however, only be achieved with reduced beam energies at HERA.

Measurement of the Generalised Structure Function xF 3
At high Q 2 , the NC cross section in e + p scattering has been observed to be significantly smaller than that in e − p scattering [3], confirming the Standard Model expectation of γZ interference. These earlier H1 data were used to obtain a first measurement of the generalised structure function xF 3 in the kinematic range 0.02 ≤ x ≤ 0.65 and 1 500 GeV 2 ≤ Q 2 ≤ 12 000 GeV 2 [3]. A similar measurement has been reported recently by ZEUS [46].
Profiting from the enlarged statistics and the reduced systematic uncertainties, the previous measurement of xF 3 [3] is updated here by using the same published e − p and the full e + p data obtained by H1 at HERA-I. Fig. 18(a) shows the comparison of the e − p and e + p data for three different Q 2 values at 1 500 GeV 2 , 5 000 GeV 2 and 12 000 GeV 2 , together with the expectations determined from the H1 PDF 2000 fit. The generalised structure function xF 3 , given in table 15, is obtained from and is compared in fig. 18(b) with the calculation. Since at high x and low Q 2 the expected sensitivity to xF 3 is smaller than the luminosity uncertainty, the measurement is not performed in this region. The dominant contribution to xF 3 arises from γZ interference, which allows xF γZ 3 to be extracted (eq. 4) according to xF γZ 3 ≃ −xF 3 (Q 2 + M 2 Z )/(a e κQ 2 ) by neglecting the pure Z exchange contribution, which is suppressed by the small vector coupling v e . This structure function is non-singlet and has little dependence on Q 2 . This is illustrated in fig. 18(c). The measured xF γZ 3 at these Q 2 values can thus be averaged taking into account the small Q 2 dependence. The two lowest x bins at x = 0.020 and 0.032 are averaged as well. The averaged xF γZ 3 , determined for a Q 2 value of 1 500 GeV 2 , is shown in table 15 and fig. 18(d) in comparison with the QCD fit result. The structure function xF γZ 3 determines both the shape and magnitude of the valence quark distributions independent of the sea quark distributions. The calculation from the QCD fit, in which the parton densities in the valence region are principally constrained by the NC and CC cross sections rather than the difference between the e ± NC cross sections, gives a good description of the measurement. The averaged structure function is integrated [47] over the measured x range, yielding

The Quark Distributions xu and xd at Large x
The flavour composition of the proton at high x may be disentangled by exploiting the NC and CC cross section measurements. The e + p CC cross section at large x is dominated by the d quark contribution as may be inferred from fig. 13. Similarly the u distribution dominates the e − p CC and e ± p NC cross sections at large x. Using data points for which the xu or xd contribution provides at least 70% of the cross section, as determined from the H1 PDF 2000 fit, the up and down quark distributions are determined locally, using the method described in [1,27]. The extraction relies on weighting the differential cross section measurement with the calculated local flavour contribution and is illustrated in fig. 19, where xu is the combined result from three independent extractions from the NC e ± p and CC e − p data and xd is determined from the CC e + p data only. This method is complementary to performing a QCD fit, since it is based on the local cross section measurements and is less sensitive to the parameterisations and dynamical assumptions used in the fit.
The extracted xu and xd distributions are further compared in fig. 19 with the results of the H1 PDF 2000 fit by subtracting xc and xs from the fitted xU and xD. The two determinations are in good agreement. They also compare well with the recent parameterisations from the MRST [32] and CTEQ [33] groups except for xu at large x = 0.65, where the results of MRST and CTEQ, being constrained mainly by the BCDMS data, yield a larger up quark distribution.

Summary
New measurements are presented of inclusive deep inelastic neutral and charged current scattering cross sections at high momentum transfers Q 2 ≥ 100 GeV 2 from recent e + p data recorded in 1999 and 2000 by the H1 experiment at HERA. This analysis, together with previous analyses of the 1994 − 1997 e + p and 1998 − 1999 e − p data, completes the H1 measurements of the inclusive cross sections at high Q 2 from the first phase of HERA operation.
The accuracy of the neutral current (NC) measurements presented here has reached the level of a few percent in the medium Q 2 range of Q 2 < 3 000 GeV 2 . The very high Q 2 NC and charged current (CC) data are still limited by their statistical precision, which is expected to improve in the high luminosity phase of HERA.
For both e + p and e − p scattering data, the region of very large inelasticity is explored, which allows a determination of the longitudinal structure function F L (x, Q 2 ) for the first time in the large momentum transfer range, 110 GeV 2 ≤ Q 2 ≤ 700 GeV 2 . The observed interference of the photon and Z exchange, differing between e + p and e − p NC scattering at high Q 2 , is used to measure the structure function xF 3 , superseding the earlier measurement.
The NC and CC cross sections in e ± p scattering are sensitive to the sums of up-and down-type quark and anti-quark distributions, xU, xD, xU and xD. Based on these quark distribution decompositions, a novel NLO QCD analysis is performed, resulting in a first determination of the partonic nucleon structure from inclusive DIS scattering data from H1 alone. The experimental precision achieved in this analysis is about 3% and 10% respectively for xU and xD at x = 0.4. The resulting parton distributions are found to be in agreement with those obtained in an analysis also including the BCDMS muon-nucleon data at large x. The QCD analysis leads to a good description of all the fitted NC and CC cross section data and of the derived structure functions over more than four orders of magnitude in x and Q 2 . Charged Current Charged Current   The NC e + p cross section dσ N C /dQ 2 for y < 0.9 after correction (k cor ) according to the Standard Model expectation determined from the H1 PDF 2000 fit for the kinematic cut E ′ e > 6 GeV for Q 2 < 890 GeV 2 . The statistical (δ stat ), uncorrelated systematic (δ unc ), correlated systematic (δ cor ) and total (δ tot ) errors are provided. In addition the correlated systematic error contributions from a positive variation of one standard deviation of the electron energy error (δ E + cor ), of the polar electron angle error (δ θ + cor ), of the hadronic energy error (δ h + cor ), of the error due to noise subtraction (δ N + cor ) and of the error due to background subtraction (δ B + cor ) are given. The normalisation uncertainty of 1.5% is not included in the errors.  The CC e + p cross section dσ CC /dQ 2 for y < 0.9 after correction (k cor ) according to the Standard Model expectation determined from the H1 PDF 2000 fit for the kinematic cuts 0.03 < y < 0.85 and P T,h > 12 GeV. The statistical (δ stat ), uncorrelated systematic (δ unc ), correlated systematic (δ cor ) and total (δ tot ) errors are also given. In addition the correlated systematic error contributions from a positive variation of one standard deviation of the error due to the cuts against photoproduction (δ V + cor ), of the hadronic energy error (δ h + cor ), of the error due to noise subtraction (δ N + cor ) and of the error due to background subtraction (δ B + cor ) are given. The normalisation uncertainty of 1.5% is not included in the errors. The last column gives the correction for QED radiative effects (∆ QED CC ).
x  The NC e + p cross section dσ N C /dx measured for y < 0.9 and Q 2 > 1 000 GeV 2 . The statistical (δ stat ), uncorrelated systematic (δ unc ), correlated systematic (δ cor ) and total (δ tot ) errors are provided. In addition the correlated systematic error contributions from a positive variation of one standard deviation of the electron energy error (δ E + cor ), of the polar electron angle error (δ θ + cor ), of the hadronic energy error (δ h + cor ), of the error due to noise subtraction (δ N + cor ) and of the error due to background subtraction (δ B + cor ) are given. The normalisation uncertainty of 1.5% is not included in the errors.  The NC e + p cross section dσ N C /dx measured for y < 0.9 and Q 2 > 10 000 GeV 2 .
The statistical (δ stat ), uncorrelated systematic (δ unc ), correlated systematic (δ cor ) and total (δ tot ) errors are provided. In addition the correlated systematic error contributions from a positive variation of one standard deviation of the electron energy error (δ E + cor ), of the polar electron angle error (δ θ + cor ), of the hadronic energy error (δ h + cor ), of the error due to noise subtraction (δ N + cor ) and of the error due to background subtraction (δ B + cor ) are given. The normalisation uncertainty of 1.5% is not included in the errors.  The CC e + p cross section dσ CC /dx for y < 0.9 and Q 2 > 1 000 GeV 2 after correction (k cor ) according to the Standard Model expectation determined from the H1 PDF 2000 fit for the kinematic cuts 0.03 < y < 0.85 and P T,h > 12 GeV. The statistical (δ stat ), uncorrelated systematic (δ unc ), correlated systematic (δ cor ) and total (δ tot ) errors are also given. In addition the correlated systematic error contributions from a positive variation of one standard deviation of the error due to the cuts against photoproduction (δ V + cor ), of the hadronic energy error (δ h + cor ), of the error due to noise subtraction (δ N + cor ) and of the error due to background subtraction (δ B + cor ) are given. The normalisation uncertainty of 1.5% is not included in the errors.  The NC e + p reduced cross sectionσ N C (x, Q 2 ), shown with statistical (δ stat ), systematic (δ sys ) and total (δ tot ) errors. Also shown are the total uncorrelated systematic (δ unc ) errors and two of its contributions: the electron energy error (δ E unc ) and the hadronic energy error (δ h unc ). The effect of the other uncorrelated systematic errors is included in (δ unc ). The table also provides the correlated systematic error (δ cor ) and its contributions from a positive variation of one standard deviation of the error on the electron energy (δ E + cor ) and polar angle (δ θ + cor ), of the hadronic energy error (δ h + cor ), of the error due to noise subtraction (δ N + cor ) and background subtraction (δ B + cor ) and of the error due to charge symmetry background subtraction in the high-y analysis (δ S + cor ). The normalisation uncertainty of 1.5% is not included in the errors. The NC structure function term scaled by the helicity factor Y + φ N C /Y + is given as well as the electromagnetic structure function F 2 with the corrections ∆ F 2 , ∆ F 3 and ∆ F L as defined in eq. 18. For Q 2 < 2 000 GeV 2 , the extraction of F 2 is restricted to the region of y < 0.6. The table continues on the next 2 pages.    The CC e + p double differential cross section d 2 σ CC /dxdQ 2 and the structure function term φ CC , shown with statistical (δ stat ), systematic (δ sys ) and total (δ tot ) errors. Also shown are the total uncorrelated systematic error (δ unc ) and its contribution from the hadronic energy error (δ h unc ). The effect of the other uncorrelated systematic errors is included in (δ unc ). In addition the correlated systematic error (δ cor ) and its contributions from a positive variation of one standard deviation of the error due to the cuts against photoproduction (δ V + cor ), of the hadronic energy error (δ h + cor ), of the error due to noise subtraction (δ N + cor ) and of the error due to background subtraction (δ B + cor ) are given. The normalisation uncertainty of 1.5% is not included in the errors. The last column gives the correction for QED radiative effects (∆ QED CC ).  The NC e − p reduced cross sectionσ N C (x, Q 2 ) from the high-y analysis, shown with statistical (δ stat ) and total (δ tot ) errors. Also shown are the total uncorrelated systematic (δ unc ) error and two of its contributions: the electron energy error (δ E unc ) and the hadronic energy error (δ h unc ). The effect of the other uncorrelated systematic errors is included in (δ unc ). In addition the correlated systematic error (δ cor ) and its contributions from a positive variation of one standard deviation of the electron energy error (δ E + cor ), of the polar electron angle error (δ θ + cor ), of the hadronic energy error (δ h + cor ), of the error due to noise subtraction (δ N + cor ), of the error due to background subtraction (δ B + cor ) and of the error due to charge symmetry background subtraction (δ S + cor ) are given. The normalisation uncertainty of 1.8% is not included in the errors. All e − p data not previously reported in [3] are given, including the new high y data and three data points at Q 2 = 100 and 120 GeV 2 from the nominal analysis phase space.  The NC structure function term φ N C (x, Q 2 ) and the structure function F L , shown with its statistical (δ stat ), systematic (δ sys ) and total (δ tot ) absolute error. The total error includes a contribution arising from the model uncertainties in the calculatedF 2 . These are obtained by varying the assumptions of the H1 Low y QCD fit as listed in table 5. The luminosity uncertainties of the e + p and e − p data sets are included in the systematic error.