Measurement of charged particle spectra in deep-inelastic ep scattering at HERA

Charged particle production in deep-inelastic ep scattering is measured with the H1 detector at HERA. The kinematic range of the analysis covers low photon virtualities, 5<Q2<100 GeV2, and small values of Bjorken-x, 10−4<x<10−2. The analysis is performed in the hadronic centre-of-mass system. The charged particle densities are measured as a function of pseudorapidity (η∗) and transverse momentum (\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$p_{T}^{*}$\end{document}) in the range 0<η∗<5 and \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$0<p_{T}^{*} < 10\mbox{~GeV}$\end{document} in bins of x and Q2. The data are compared to predictions from different Monte Carlo generators implementing various options for hadronisation and parton evolutions.

Abstract Charged particle production in deep-inelastic ep scattering is measured with the H1 detector at HERA. The kinematic range of the analysis covers low photon virtualities, 5 < Q 2 < 100 GeV 2 , and small values of Bjorken-x, 10 −4 < x < 10 −2 . The analysis is performed in the hadronic centre-of-mass system. The charged particle densities are measured as a function of pseudorapidity (η * ) and transverse momentum (p * T ) in the range 0 < η * < 5 and 0 < p * T < 10 GeV in bins of x and Q 2 . The data are compared to predictions from different Monte Carlo generators implementing various options for hadronisation and parton evolutions. a e-mail: daum@mail.desy.de b Also at Rechenzentrum, Universität Wuppertal, Wuppertal, Germany.

Introduction
Deep-inelastic scattering (DIS) processes at the ep collider HERA can access small values of Bjorken-x at low four momentum transfers squared Q 2 of a few GeV 2 . In the region of low x, characterised by high densities of gluons and sea quarks in the proton, the parton interaction with the virtual photon may originate from a cascade of partons emitted prior to the interaction as illustrated in Fig. 1. In perturbative Quantum Chromodynamics (QCD) such multi-parton emissions are described only within certain approximations valid in restricted phase space regions. At sufficiently large Q 2 and not too small x the Dokshitzer-Gribov-Lipatov-Altarelli-Parisi (DGLAP) [1][2][3][4][5] evolution equation is expected to be a good approximation. The DGLAP equation corresponds to a strong ordering of the transverse momenta of the propagator partons, k T ,i , with respect to the proton direction, which implies strong ordering of the transverse momenta of the emitted partons, p T   At small x the DGLAP approximation is expected to become inadequate and the Balitsky-Fadin-Kuraev-Lipatov (BFKL) [6][7][8] scheme may be more appropriate, which has no ordering in k T of the partons along the ladder. The Ciafaloni-Catani-Fiorani-Marchesini (CCFM) evolution [9][10][11][12] aims to unify the DGLAP and BFKL approaches. It introduces angular ordering of gluon emissions to implement coherence effects. At small x the CCFM evolution equation is almost equivalent to the BFKL approach, while it reproduces the DGLAP equations for sufficiently large x and Q 2 . Measurements of the proton structure function F 2 (x, Q 2 ) [13,14] are well described by the Next-to-Leading-Order (NLO) or Next-to-Next-to-Leading Order (NNLO) DGLAP evolution [15][16][17][18][19][20][21][22][23][24][25], suggesting that this observable may be too inclusive to exhibit signals for BFKL dynamics. Deviations from the k T ordering at HERA are observed in jet production [26,27], transverse energy flow [28,29], forward jet production [30][31][32] and measurements of forward π 0 production [33]. Studies of the transverse momentum spectrum of charged particles have been proposed in [34] as a more direct probe of the underlying parton dynamics. In that paper it has been shown with the help of QCD models that the high-p T tail is sensitive to parton radiation, while the contribution from hadronisation is small. Previous measurements of the transverse momentum and pseudorapidity dependence of particle densities performed by the H1 collaboration [35] were limited in statistical precision. This paper presents a study of charged particle production in ep collisions for 5 < Q 2 < 100 GeV 2 . The analysis is performed in the hadronic centre-of-mass frame, i.e. in the virtual photon-proton rest frame. The charged particle densities as a function of pseudorapidity and transverse momentum are measured differentially in x and Q 2 .
The data set used for the analysis was collected with the H1 detector in 2006 when positrons and protons collided with energies of 27.6 GeV and 920 GeV, respectively, corresponding to a centre-of-mass energy of √ s = 319 GeV. The integrated luminosity of the data set is 88.6 pb −1 , which is about seventy times larger than what was available for the previous H1 measurement [35]. This allows for a more detailed study of the dynamical features of parton evolution in the proton at small x.

QCD models
Parton cascade and hadronisation processes leading to charged particle production in ep collisions are modeled using different Monte Carlo (MC) programs. Brief descriptions of the MC event generators considered in this analysis are given below.
• The RAPGAP [36] MC generator matches first order QCD matrix elements to DGLAP based leading logarithm approximations for parton showers with strongly ordered transverse momenta of subsequently emitted partons. The factorisation and renormalisation scales are set to μ f = μ r = Q 2 +p T 2 , wherep T is the transverse momentum of the outgoing hard parton from the matrix element in the centre-of-mass frame of the hard subsystem. • The DJANGOH [37] MC generator uses the Colour Dipole Model (CDM) as implemented in ARIADNE [38], which models first order QCD processes and creates dipoles between coloured partons. Gluon emission is treated as radiation from these dipoles, and new dipoles are formed from the emitted gluons from which further radiation is possible. The radiation pattern of the dipoles includes interference effects, thus modelling gluon coherence. The transverse momenta of the emitted partons are not ordered, producing a configuration similar to the BFKL treatment of parton evolution [39,40]. • The CASCADE [41,42] MC generator uses off-shell leading order QCD matrix elements, supplemented with parton emissions based on the CCFM evolution equation. The equation requires an unintegrated gluon density (see [41,42]), which takes the transverse momenta of the propagators into account. In contrast to the DGLAP evolution equation, the CCFM equation only contains gluon splitting g → gg.
• The Herwig++ [43] MC program with the POWHEG (POsitive Weight Hardest Emission) option [44], combines the full matrix element including virtual corrections at O(α s ) with a DGLAP-like parton shower simulation. The Herwig++ MC program uses the Coherent Parton Branching algorithm [45,46], which implements coherence via ordering inq which is defined asq = p T /(1 − z) andq = p T /(z(1 − z)) for initial and final state showering, respectively [47]. Here (1 − z) is the light-cone momentum fraction of the emitted parton and p T is the transverse momentum of the splitting products. • Photoproduction background is generated with the PHO-JET [48,49] program, which uses a two-component dual parton model [50] including diffractive processes and vector meson production.
The RAPGAP and DJANGOH predictions are calculated using the CTEQ6L(LO) [51] set of Parton Distribution Functions (PDF), while CASCADE uses the default A0 unintegrated gluon density set [52]. The predictions of Herwig++ were obtained with the default PDF MRST 02 NLO [53]. To simulate hadronisation the Lund string fragmentation model [54] is used, as implemented in JETSET [55] for DJANGOH and PYTHIA [56] for both RAPGAP and CASCADE. The parameters of the Lund string fragmentation model used here are tuned to describe e + e − results [57]. The tuning was performed by the ALEPH collaboration using hadronic Z decay data and the PYTHIA6.1 simulation with Bose-Einstein correlations turned on. In addition, the tune obtained by the Professor tool [58] using LEP data is also tested. Her-wig++ incorporates the cluster model [59,60] of hadronisation, in which colour-singlet clusters of partons form after the perturbative phase and then decay into the observed hadrons.
DJANGOH and RAPGAP are also used together with the H1 detector simulation in order to determine the acceptance and efficiency and to estimate the systematic uncertainties associated with the measurement. The programs are interfaced to HERACLES [61] to simulate the QEDradiative effects. The generated events are passed through a detailed simulation of the H1 detector response based on the GEANT simulation program [62] and are processed using the same reconstruction and analysis program chain as for data. For the determination of the detector effects both the RAPGAP and DJANGOH predictions are studied. Both models describe all relevant control distributions reasonably well [63]. To improve the determination of the detector corrections the transverse momentum, pseudorapidity and multiplicity of charged particles as well as inelasticity y, defined as y = Q 2 /(s · x), are reweighted to the data [63]. The reweighting is applied to the generator quantities.

H1 detector
A full description of the H1 detector can be found elsewhere [64][65][66] and only the components most relevant for this analysis are briefly mentioned here. The coordinate system of H1 is defined such that the positive z axis is pointing in the direction of the proton beam (forward direction) and the nominal interaction point is located at z = 0. The polar angle θ is then defined with respect to this axis. The pseudorapidity is defined to be η = − ln(tan(θ/2)).
Charged particles are measured within the central tracking detector (CTD) in the polar angle range 20 • < θ < 165 • , which is also used to reconstruct the interaction vertex. The CTD comprises two large cylindrical jet chambers (CJCs), and the silicon vertex detector [67,68]. The CTD is operated inside a 1.16 T solenoidal magnetic field. The CJCs are separated by a cylindrical drift chamber which improves the z coordinate reconstruction. A cylindrical multiwire proportional chamber [69], which is mainly used in the trigger, is situated inside the inner CJC. The trajectories of charged particles are measured with a transverse momentum resolution of σ (p T )/p T ≈ 0.2 %p T /GeV ⊕ 1.5 %. The forward tracking detector (FTD) [70] measures the tracks of charged particles at polar angles 6 • < θ < 25 • . In the region of angular overlap, FTD and short CTD track segments are used to reconstruct combined tracks, extending the detector acceptance for well-reconstructed tracks. Both the CTD and the combined tracks are linked to hits in the vertex detectors: the central silicon tracker (CST) [67,68], the backward silicon tracker (BST) [71] and the forward silicon tracker (FST) [72]. These detectors provide precise spatial coordinate measurements and therefore significantly improve the primary vertex spatial resolution. The CST consists of two layers of double-sided silicon strip detectors surrounding the beam pipe covering an angular range of 30 • < θ < 150 • for tracks passing through both layers. The BST consists of six double wheels of strip detectors measuring the transverse coordinates of charged particles. The FST design is similar to the BST and consists of five double wheels of single-sided strip detectors.
The lead-scintillating fibre calorimeter (SpaCal) [66] covering the region 153 • < θ < 177.5 • , has electromagnetic and hadronic sections. The calorimeter is used to measure the scattered positron and the backward hadronic energy flow. The energy resolution for positrons in the electromagnetic section is σ (E)/E ≈ 7.1 %/ √ E/GeV ⊕ 1 %, as determined in test beam measurements [73]. The SpaCal provides energy and time-of-flight information used for triggering purposes. A backward proportional chamber (BPC) in front of the SpaCal is used to improve the angular measurement of the scattered lepton. The liquid argon (LAr) calorimeter [74] covers the range 4 • < θ < 154 • and is used in this analysis in the reconstruction of the hadronic final state. It has an energy resolution of σ (E)/E ≈ 50 %/ √ E/GeV ⊕ 2 % for hadronic showers, as obtained from test beam measurements [75].

Event reconstruction
The DIS kinematics is reconstructed based on the measurement of the scattered electron and the hadronic final state (HFS) particles. In the so-called eΣ-method [76] the kinematic variables Q 2 , y and x are given by: where s is the square of the centre-of-mass energy, E e and θ e the energy and polar angle of the scattered lepton, respectively, E e being the energy of incoming lepton and where the sum runs over all hadronic final state (HFS) particles i. This method provides an optimum in resolution of the kinematic variables and shows only little sensitivity to QED radiative effects. The HFS particles are reconstructed using an energy flow algorithm [77][78][79]. This algorithm combines charged particle tracks and calorimetric energy clusters, taking into account their respective resolution and geometric overlap, into hadronic objects, while avoiding double counting of energy.

Data selection
DIS events were recorded using triggers based on electromagnetic energy deposits in the SpaCal calorimeter. The trigger efficiency is determined using independently triggered data. For DIS events the trigger inefficiency is negligible in the kinematic region of the analysis. The scattered lepton, defined by the most energetic SpaCal cluster, is required to have an energy E e larger than 12 GeV. The kinematical phase space is defined by 5 < Q 2 < 100 GeV 2 and 0.05 < y < 0.6, corresponding to the geometric acceptance of the SpaCal. The upper cut on y reduces background from photoproduction. In addition, x is required to be in the range of 0.0001 < x < 0.01.
Additional selections are made to reduce QED radiation effects and to suppress background events. The z coordinate of the event vertex is required to be within 35 cm of the nominal interaction point. Events with high energy initial state photon radiation are rejected by requiring 35 < i (E i − p z,i ) < 75 GeV. Here, the sum extends over all HFS particles and the scattered electron. This cut further suppresses photoproduction background events to a level of about 0.5 %.
The tracks used in the analysis are measured in the CTD alone (central tracks) or result from combinations of CTD and FTD information (combined tracks). Central tracks are required to have transverse momenta in the laboratory frame p T > 150 MeV. The momentum of a combined track is required to be larger than 0.5 GeV to ensure that the track has enough momentum to cross the endwall of the CJC. Both central and combined tracks are required to originate from the primary event vertex and to be in the pseudorapidity range −2 < η < 2.5 measured in the laboratory frame. Using only tracks assigned to the event vertex, the contributions from in-flight decays of K 0 S , Λ and from photon conversions and from other secondary decays are reduced. Further track quality cuts [63] are applied to ensure a high purity of the track reconstruction.

Definition of experimental observables
The results of this analysis are presented in the hadronic centre-of-mass frame (HCM), to minimise the effect of the transverse boost from the virtual photon. The transformation to the HCM frame is reconstructed with the knowledge of the kinematic variables Q 2 and y [63]. The transverse momentum and pseudorapidity of charged particles in the HCM frame are labelled as p * T and η * . Since in this frame the positive z * axis is defined by the direction of the virtual photon, HFS particles with η * > 0 belong to the current hemisphere and particles with η * < 0 originate from the target (proton remnant) hemisphere.
Charged particle densities as a function of transverse momentum and pseudorapidity are defined as (1/N )(dn/dp * T ) and (1/N )(dn/dη * ), respectively. Here, dn is the total number of charged particles with transverse momentum (pseudorapidity) in the dp * T (dη * ) bin and N denotes the number of selected DIS events. For distributions measured differentially in x and Q 2 , dn and N are the numbers for the respective (x, Q 2 ) bin.
Hadronisation is expected to be more relevant at small transverse momenta, while the hard parton radiation is expected to contribute more significantly at high p * T (p * T > 1 GeV) [34]. To distinguish hadronisations effects from parton evolution signatures, the charged particle density is measured as a function of η * for 0 < p * T < 1 GeV and for 1 < p * T < 10 GeV. The p * T dependence of the charged particle densities is studied in two different pseudorapidity intervals, 0 < η * < 1.5 and 1.5 < η * < 5, referred to as the "central region" and "current region", respectively, as illustrated in Fig. 2. Such division approximately defines the regions where the sensitivity to the hard scatter is largest (current region), and where the parton shower models can be tested (central region). The target region, η * < 0, is not accessible in this analysis. Fig. 2 The two pseudorapidity regions analysed in this paper. The region 0 < η * < 1.5 and 1.5 < η * < 5, are denoted as "central" and "current" regions, respectively

Data corrections
The data are corrected to the number of stable charged particles including charged hyperons, with proper lifetime cτ > 10 mm, in the phase space given in Table 1. Correction factors are calculated for each analysis bin from the ratio of the number of generated stable charged particles to the number of reconstructed tracks. The bin widths are chosen such that a purity 1 of more than 75 % is ensured in all bins. The correction takes into account detector effects like limited resolution and losses near the phase space boundaries, as well as a small residual contamination from weak decays of neutral particles (e.g. K 0 and Λ).
In addition to migrations between bins inside the measurement phase space, there are migrations from outside of the analysis phase space and there is background from photoproduction. These contributions are subtracted prior to applying the correction factors according to the procedure outlined in [32].
The DJANGOH MC was used to correct the data. The differences to the correction factors obtained from RAPGAP Transverse momentum in the laboratory frame p T > 150 MeV Pseudorapidity in the HCM frame 0 < η * < 5 Transverse momentum in the HCM frame 0 < p * T < 10 GeV 1 The purity is defined as the ratio of the number of charged particles generated and reconstructed in a given bin to the total number of charged particles in the phase space of the analysis which are reconstructed in this bin.
are taken as systematic uncertainties. The correction factors strongly depend on η * and to a lesser extent on p * T . In the 1.5 < η * < 5 region they vary between 1 and 1.8 with the largest values seen at high p * T and large η * . In the 0 < η * < 1.5 region, the correction factors rise up to 2.6 at high p * T , due to the limited detector acceptance in this region. The two MC models predict very similar correction factors for most of the phase space region, but differences up to 5.5 % are observed at small η * and high p * T .

Systematic uncertainties
The following sources of systematic uncertainties are considered for all measured quantities.
• The systematic uncertainty on the SpaCal energy scale is 1 % [80], which results in a systematic uncertainty of typically 0.4 % for the measured single differential distributions. • The SpaCal angular resolution of 1 mrad leads to a systematic uncertainty of about 0.1 % for the measured distributions. • The hadronic energy scale uncertainty is known to a precision of 2 % [81]. The scale uncertainty enters into the uncertainty of the phase space calculations, which depend on E − P z of the HFS, and also affects the boost to the HCM frame. The systematic effect on the present measurements is about 0.3 %. • The systematic uncertainty arising from the model dependence of the data correction is taken as the difference of the correction factors calculated using RAPGAP and DJANGOH MC. The resulting uncertainty on the measurements varies between 0.2 % and 5.5 %. • The systematic uncertainty associated with the track reconstruction (e.g. track reconstruction efficiency, vertex reconstruction efficiency, weak decays and nuclear interaction uncertainties) is estimated to be: -1 % for central tracks. The track reconstruction efficiency was determined from the analysis of curling tracks and of the decay angle distribution of charged pions from the decay K 0 S → π + π − . Data agree with the simulation within 0.5 % for p T > 100 MeV [82]. Nuclear interactions of pions and kaons have been investigated by the analysis of secondary vertices of charged particles located in the material between the two CJCs and originating from interactions with the detector material. The nuclear interaction cross sections for these hadrons are found to be smaller in the simulation than in data. After correcting for these deficits, the agreement between data and MC in the overall track efficiency is found to be better than 1 %. -10 % for combined tracks [63,70]. This was checked using all selected central tracks, as well as by using pions from K 0 S decays, as a function of transverse momentum and pseudorapidity. Consistent results are obtained from both samples showing agreement of data and MC within 10 %. The systematic uncertainty associated with the track reconstruction is applied as an independent uncorrelated uncertainty on every data point. The resulting uncertainty on the measurements varies between 1 % and 5.4 % and is 1.6 % on average. An additional systematic uncertainty of 0.2 % is assigned due to the different K 0 S contamination seen in data and MC for both central and combined tracks. The corresponding effect arising from Λ contamination is expected to be negligible.
• The systematic uncertainty on the remaining photoproduction background is estimated to be 30 %. This results in an uncertainty on the measured densities up to 0.9 % at small x and Q 2 . At large x and Q 2 the contribution from photoproduction is small and its uncertainty is negligible.
The dominant systematic uncertainties are the uncertainties arising from the model dependence of the data correction, which are largely correlated, and uncertainties associated with the track reconstruction, which are uncorrelated. The systematic uncertainties shown in the figures and tables are calculated by adding all contributions in quadrature. The total systematic uncertainty for the single differential measurements is below 2.5 % for most analysis bins.

Results
The measurements of the charged particle densities as a function of pseudorapidity and transverse momentum in the phase space summarised in Table 1 are listed in Tables 2 to 9 and shown in Figs. 3 to 10. 6.1 Charged particle densities as a function of pseudorapidity The charged particle densities as a function of η * were measured separately for p * T < 1 GeV and for 1 < p * T < 10 GeV, as shown in Fig. 3. In the soft p * T region, the pseudorapidity distribution is almost flat in the 1.5 η * 3 range with about 1.7 charged particles per unit of pseudorapidity. The distribution falls at small η * due to the cut on pseudorapidity in the laboratory frame. In the hard p * T region the distribution becomes more peaked near η * = 2.5, with a maximum of 0.23 charged particles per unit of pseudorapidity. For 1 < p * T < 10 GeV the density increases rather strongly up to η * ≈ 2.5, a behaviour expected from the strong ordering of transverse momentum towards the hard scattering vertex. Table 2 Charged particle densities as a function of η * for 0 < p * T < 1 GeV with relative statistical (stat.) and systematic (sys.) uncertainties given in per cent. The phase space is defined in Table 1 0 < p  Table 3 Charged particle densities as a function of η * for 1 < p * T < 10 GeV with relative statistical (stat.) and systematic (sys.) uncertainties given in per cent. The phase space is defined in Table 1 1 < p * T < 10 GeV η * range 1/N · dn/dη * stat. show similar results and predict less particles compared to calculations using the default LO PDF set CTEQ6L(LO). All predictions are close to the data. At large p * T , differences between the NLO PDF sets are observed, with CTEQ6L(LO) being closest to the data, although the differences to the data are still larger than the differences between the various PDF predictions. Similar PDF uncertainties are observed when using the CDM model as implemented in DJANGOH.
To check the sensitivity to hadronisation effects, the RAPGAP predictions obtained with three sets of fragmentation parameters are compared to the data in Fig. 4: parameters tuned by ALEPH [57], by the Professor tuning tool [58] and default PYTHIA6.424 fragmentation parameters. Significant differences between these three samples are seen in the soft p * T region, where the data are best described by the ALEPH tune. At large transverse momenta they give similar results but none of them describes the data.
Predictions from models with different approaches for QCD radiation (see Sect. 2) are shown in Fig. 5. The data Table 4 Charged particle densities as a function of η * for 0 < p * T < 1 GeV for different Q 2 and x intervals with relative statistical (stat.) and systematic (sys.) uncertainties given in per cent. The phase space is defined in Table 1 0 < p   Table 5 Charged particle densities as a function of η * for 1 < p * T < 10 GeV for different Q 2 and x intervals with relative statistical (stat.) and systematic (sys.) uncertainties given in per cent. The phase space is defined in Table 1 1 < p * T < 10 GeV    Table 6 Charged particle densities as a function as a function of p * T in the region 0 < η * < 1.5 shown with relative statistical (stat.) and systematic (sys.) uncertainties given in per cent. The phase space is defined in Table 1 0 < η * < 1.5  Table 7 Charged particle densities as a function as a function of p * T in the region 1.5 < η * < 5 shown with relative statistical (stat.) and systematic (sys.) uncertainties given in per cent. The phase space is defined in Table 1 1.5 < η * < 5 p * T range [GeV] 1/N · dn/dp * T [GeV  Table 8 Charged particle densities as a function of p * T in the region 0 < η * < 1.5 for different Q 2 and x intervals shown with relative statistical (stat.) and systematic (sys.) uncertainties given in per cent. The phase space is defined in Table 1 0 < η * < 1.5 1/N · dn/dp * In the soft p * T region, DJANGOH and RAPGAP describe the data within the PDF uncertainties (Fig. 3). Her-wig++, which uses the cluster fragmentation model, provides a reasonable description of the data in the central region. The effect of not using the POWHEG option in Her-wig++ also has been investigated. Only small differences were observed which are not considered further in this paper. CASCADE predicts too high multiplicities in most of the measured η * range. In the region of 1 < p * T < 10 GeV the best description of the data is achieved by DJANGOH. RAP-GAP strongly undershoots the data in the central region. Her-wig++ predicts a spectrum which is even below the prediction of RAPGAP. CASCADE is significantly above the data in a wide range of η * .
The charged particle densities as a function of η * are shown in Fig. 6 for p * T < 1 GeV in eight (x, Q 2 ) bins. The data are compared to predictions of the DJANGOH, RAP-GAP, Herwig++ and CASCADE generators. DJANGOH provides a good description of the data over the full kinematic Table 9 Charged particle densities as a function of p * T in the region 1.5 < η * < 5 for different Q 2 and x intervals shown with relative statistical (stat.) and systematic (sys.) uncertainties given in per cent. The phase space is defined in Table 1 1.5 < η * < 5 1/N · dn/dp * T  range. In general the description of the data by RAPGAP is somewhat worse, with overshooting the data by about 10 % at low x. Herwig++ predicts smaller charged particle densities than observed in data in most of the phase space with differences of the order of 10 % at the highest Q 2 . CAS-CADE is significantly above the data for η * < 3 in all (x, Q 2 ) bins.
In Fig. 7 the charged particle densities as a function of η * are shown in (x, Q 2 ) bins for 1 < p * T < 10 GeV. The shape of the distributions changes with x and Q 2 more strongly than what is observed for p * T < 1 GeV (Fig. 6). At small values of x and Q 2 the measured distribution is less dependent on η * compared to the region at high x and Q 2 . None of the models describes all aspects of the data. In general DJANGOH is closest to the data. However it fails to describe the data at low and medium x in the central pseudorapidity region, with downwards deviations of the order of 20 %. The RAPGAP prediction is below the data, with the strongest deviation observed at small x and small η * . Her-wig++ significantly undershoots the data. The prediction of CASCADE agrees reasonably well with the measurement at low x and Q 2 , but overshoots the data significantly as x or Q 2 increases. 6.2 Charged particle densities as a function of transverse momentum In Fig. 8 the charged particle densities as a function of p * T are shown for two pseudorapidity intervals, 0 < η * < 1.5 (central) and 1.5 < η * < 5 (current). The shapes of the measured p * T distributions in the two pseudorapidity ranges are similar. The spectrum falls by more than four orders of magnitude from p * T < 1 GeV to p * T ≈ 8 GeV. The measurements are compared to the predictions of the DJANGOH, RAPGAP, Herwig++ and CASCADE generators. DJANGOH provides in general a good description of the data, while only at high p * T in the current region deviations from the measurement are observed. The other models fail to describe the data, with the strongest deviations being observed in the central region. The ratio of RAPGAP to data shows a sharp drop at p * T ≈ 1 GeV. The p * T spectra predicted by Herwig++ are even softer than those predicted by RAPGAP. CASCADE in general produces higher particle densities than measured.
In Figs. 9 and 10 the charged particle densities as a function of p * T are shown for eight (x, Q 2 ) bins for the central and the current region, respectively. In the central region the measurement shows a dependence on x, such that the number of soft particles is decreasing with increasing x for fixed Q 2 . In the current region this effect is less pronounced. The DJANGOH model provides in general a good description of the data over the full kinematic range in both pseudorapidity Fig. 3 Charged particle density as a function of η * for (a) p * T < 1 GeV and for (b) 1 < p * T < 10 GeV compared to RAPGAP predictions with different proton PDFs. The predictions are obtained using the ALEPH tune Fig. 4 Charged particle density as a function of η * for (a) p * T < 1 GeV for (b) 1 < p * T < 10 GeV compared to RAPGAP predictions for three different sets of fragmentation parameters. The predictions are obtained using CTEQ6L(LO) PDF regions, degrading at high p * T in the lowest (x, Q 2 ) bin. Significant deviations of the RAPGAP predictions from data are observed in the central region at low x and Q 2 . The description becomes somewhat better at larger values of x and Q 2 . The same trend is observed for the current region, but the overall data description is better. Herwig++ fails to describe the measurements at high p * T in the whole phase space. At lowest x and Q 2 the spectrum is much softer than the one obtained with RAPGAP, while at high x and Q 2 both predictions are similar. CASCADE describes the data in the lowest x and Q 2 bin at high p * T only.

Conclusion
This paper presents a study of charged particle production in ep collisions at low Q 2 measured with the H1 detector. The kinematic range of the analysis covers low photon virtualities, 5 < Q 2 < 100 GeV 2 , and small values of x, 10 −4 < x < 10 −2 . The analysis is performed in the hadronic centre-of-mass system. The charged particle densities as a function of pseudorapidity (η * ) and transverse momentum (p * T ) are measured in bins of x and Q 2 . The charged particle densities as a function of pseudorapidity show different   6 Charged particle density as a function of η * for p * T < 1 GeV for eight intervals of Q 2 and x compared to DJANGOH, RAPGAP, Herwig++ and CASCADE Monte Carlo predictions Fig. 7 Charged particle density as a function of η * for 1 < p * T < 10 GeV for eight intervals of Q 2 and x compared to DJANGOH, RAPGAP, Herwig++ and CASCADE Monte Carlo predictions Fig. 8 Charged particle density as a function of p * T in the ranges (a) 0 < η * < 1.5 and (b) 1.5 < η * < 5 compared to DJANGOH, RAPGAP, Herwig++ and CASCADE Monte Carlo predictions. The ratios of MC predictions to the measurements are shown on the bottom of the figure  Fig. 9 Charged particle density as a function of p * T in the range 0 < η * < 1.5 for eight intervals of Q 2 and x compared to DJANGOH, RAPGAP, Herwig++ and CASCADE Monte Carlo predictions shapes, depending on the p * T range. For 0 < p * T < 1 GeV the density of particles is approximately constant for 1 < η * < 3.5, while for 1 < p * T < 10 GeV the density increases with increasing η * up to η * ≈ 2.5, a behaviour expected from the strong ordering of transverse momentum towards the hard scattering vertex. The charged particle densities as a function of transverse momentum show an x dependence at small η * (0 < η * < 1.5), such that the number of soft particles is decreasing with increasing x, while in the 1.5 < η * < 5 range this effect is less visible.
In order to relate the charged hadron spectra to the parton dynamics at small x, the data are compared to QCD models with different evolution approaches for simulating the parton cascade and with different hadronisation schemes. The data allow the validity of different models to be tested. At small p * T , the data are reasonably well described by DJAN-GOH (based on the Colour Dipole Model), as well as by RAPGAP (based on the DGLAP shower evolution). At high p * T and at low η * , RAPGAP severely undershoots the data. The differences are most pronounced at lowest x and Q 2 , and decrease with increasing x and Q 2 values. Herwig++ which is also based on DGLAP but uses a cluster fragmentation model reasonably describes the low p * T region but undershoots the data even more than RAPGAP at high p * T . CASCADE (based on CCFM) gives a reasonable description only at the lowest x and Q 2 , but overall predicts higher charged particle densities than observed in data. The Colour Dipole Model implemented in DJANGOH is the best among the considered models and provides a reasonable description of the data. Fig. 10 Charged particle density as a function of p * T in the range 1.5 < η * < 5 for eight intervals of Q 2 and x compared to DJANGOH, RAPGAP, Herwig++ and CASCADE Monte Carlo predictions