You will be presented with the full text of a science research paper. Please write an abstract for this paper. Your response should include the abstract and no additional text.

Paper text:
the predictive power of qcd lies in the factorization theorem . in deep inelastic scattering ( dis )
factorization in short and long distance parts allows us to describe the observed hadrons as a convolution of the partonic processes with non - perturbative parton densities and fragmentation functions @xcite .
single meson production in electron proton scattering @xmath0 occurs partonically already in the absence of strong interactions ( @xmath1 ) , where one parton of the proton ( a quark ) interacts with the leptonic current and fragments into a meson ( @xmath2 ) ( naive parton model ) .
+ since we are interested in perturbative qcd effects we require the meson to carry a non - vanishing transversal momentum with respect to the centre - of - mass frame of virtual vector boson coming from the electron and initial proton ( @xmath3 ) .
thus at partonic level at least two final state partons are required to balance the transversal momentum .
the leading order processes with non - vanishing transversal momentum of a fragmenting parton into a meson ( @xmath4 ) are    1 .
@xmath5 2 .
@xmath6 ,    where the virtual photon originates from the electron current .
if the virtuality of the photon @xmath7 ( @xmath8 ) is not too large compared to the squared z - mass , the contribution from z boson exchange is suppressed and may be neglected .
dealing then with c - invariant processes we do not have to calculate separately the processes with interchanged quarks and antiquarks .
factorization , proven generally for dis processes @xcite , describes the events as a convolution of the hard scattering processes with parton densities and fragmentation functions .
the investigation of single hadron production is interesting for several reasons : first of all it is a test of perturbative qcd and factorization .
the predictions depend beside the perturbative partonic calculation essentially on universal parton densities and fragmentation functions .
especially fragmentation functions , fitted to electron positron annihilation data , may be tested , in particular their universality .
further , the predictions allow for a direct comparison with experimental data , in particular there is no need to use any kind of monte carlo procedure .
thus we may expect very meaningful results .
the predictions are due to process ( b ) directly sensitive to the gluon density in the proton and may allow us to draw conclusions concerning the gluon density in the proton . from the experimental side
precise data are available from the hera collider at desy . for instance
@xmath9 mesons were measured in the forward region ( with small angles with respect to the proton remnant ) based on events detected in the h1 detector @xcite and charged hadrons were measured at the zeus experiment @xcite .    in 1978 the process ( [ eqep ] ) with non - vanishing transversal momentum of the hadron ( @xmath10 ) was calculated by mndez @xcite at tree level accuracy ( @xmath4 ) . since qcd corrections are typically large and we are confronted with precise experimental data it is desirable to compare these data with predictions of at least next - to - leading order ( nlo ) ( @xmath11 ) accuracy . also nlo - qcd predictions were computed with the assumption of purely transversal photons , neglecting the longitudinal degrees of freedom of the exchanged virtual photon @xcite .    here
we present a nlo - qcd calculation based on the dipole subtraction formalism @xcite .
in contrast to the more conventional phase space slicing method there is no need to introduce any unphysical parameter to cut the phase space in soft , respectively collinear regions . also
all cancellations of infrared singularities occur before any numerical phase space integration is performed .
thus we may present numerically very stable predictions .
the differential cross section for process ( [ eqep ] ) reads as a convolution with the parton densities and fragmentation functions as @xmath12 where , as usual , the variables @xmath13 , @xmath14 , and @xmath15 are defined as @xmath16 , @xmath17 , @xmath18 with respect to the partonic momenta and the bar quantities @xmath19 , @xmath20 , @xmath21 with respect to the momenta of the hadrons ( with @xmath22 ) .
the angle @xmath23 denotes the azimuthal angle between the planes defined on one hand by the directions of the leptons and on the other hand by the momenta of final state hadron and virtual vector boson in the centre - of - mass frame of vector boson and initial parton .
the sum is over the different initial ( @xmath24 ) and final ( @xmath25 ) state partons and @xmath26 and @xmath27 denote the corresponding parton density respectively fragmentation function .
the hard scattering process may be written as a contraction of a lepton tensor ( @xmath28 ) with a hadron tensor ( @xmath29 ) : @xmath30 if we consider the centre - of - mass system of virtual photon and initial parton , both unpolarized , there can not be any dependence on the azimuthal angle @xmath23 . integrating out this angle dependence we find a decomposition into a transversal and a longitudinal part of the virtual photon :
@xmath31 the computation of the correction to process ( [ eqep ] ) were carried out in the subtraction formalism  @xcite .
the general idea of the subtraction formalism is to subtract from the real correction an artificial counterterm which has the same pointwise singular behaviour in @xmath32 dimensions as the real correction itself .
thus the limit @xmath33 can be performed and the real phase space integral can be evaluated numerically .
the artificial counterterm is constructed in a way that it also can be integrated over the one - parton subspace analytically leading to @xmath34  poles . adding these terms to the virtual part of the correction these poles cancel all singularities in the virtual part analytically and the remaining integration over the phase space
can be carried out numerically .
the advantage compared to phase space slicing methods is that all singularities cancel _ before _ any numerical integration is performed .
also there is no need to introduce an unphysical cut parameter which in phase space slicing methods separates soft , respectively collinear phase space regions from the remaining hard region .    in our case of a convolution of a partonic cross section with distribution functions
additional kinematic constraints have to be taken into account .
the momenta of the partons which enter the convolution have to be kept fixed ( called _ identified _ partons ) .
this leads to modified artificial counterterms .
its analytical integration gives collinear singularities which cancel the singularities of the non - perturbative distribution functions yielding scheme and scale dependent convergent parts .
the whole calculation was done with help of the algebra package form  @xcite . at @xmath11
the real correction is given by the squared matrix elements of the diagrams    1 .
@xmath35 2 .
@xmath36 3 .
@xmath37 ,    where in turn each final state parton serves as an observed hadron . in process ( c )
we have to consider two flavours in the fermion traces .
this yields 4 diagrams for every pair of different flavours in contrast to 8 diagrams in the case of one uniform flavour .
care must be taken to adjust statistical factors properly .
the artificial counterterm was constructed as described by the subtraction formalism .
the phase space integral over the 3 particle final state can be performed yielding a finite real contribution .
the virtual contribution of the correction , i.e. the interference term of the born matrix elements and the one - loop matrix elements are computed . here
we encounter 2-point , 3-point , and 4-point tensor integral contributions which were reduced to scalar integrals via tensor reduction @xcite .
the scalar integrals , containing ultraviolet and infrared singularities , were computed analytically in dimensional regularization .
the analytic expressions were compared with the literature @xcite .
the virtual contribution was renormalized in a mixed scheme , where the wave functions were renormalized on - shell and the strong coupling constant in the @xmath38  scheme , yielding an ultraviolet finite virtual contribution .
the infrared singularities cancel exactly the contributions given by the integrated artificial counterterm in the subtraction formalism .    from the subtraction of remaining singularities into the parton densities and the fragmentation functions we obtain also finite remainders which depend on the factorization scheme and on the factorization scale .
in this context we choose the @xmath38  scheme .
thus we end up with three contributions , the real part , the virtual part , and the part related to identified partons .
all these finite parts have to be integrated over the 2 respectively 3 particle final state phase space . to this purpose
a c - routine was written to perform these integrations numerically .
in the convolution ( [ eqconv ] ) we use the parton densities published by the cteq collaboration @xcite .
herein the parton distribution set , called cteq5 m , where m denotes the @xmath38 scheme , matches our conditions with the assumption of 5 light quarks .
we adopt the kkp fragmentation functions to our calculation @xcite .
since these fragmentation functions are fitted to @xmath39 data with high accuracy and applied here to a dis process the comparison with experimental data will serve as a good check of universality .
we use the value of the strong coupling constant at the @xmath40-scale as given by the particle data group  @xcite @xmath41 and evolve this value with the nlo evolution equation . in the lo approximation
we evolve this value with the lo evolution equation for consistency reasons .    in fig .
[ diffcross ] we show the differential cross section @xmath42 as a function of @xmath10 for fixed values of @xmath43 , @xmath44 , @xmath45 , and @xmath46  gev@xmath47 , where @xmath10 is defined in the centre - of - mass frame of vector boson and initial parton .
we set the renormalization scale ( @xmath48 ) equal to both the initial ( @xmath49 ) and the final ( @xmath50 ) fragmentation scales @xmath51 , where the initial scale is related to the parton densities and the final scale to the fragmentation functions .
the steep fall of the cross section over several orders of magnitude makes it hard to resolve the nlo contribution ( full line ) in contrast to the lo prediction ( dashed line ) graphically .
the lower part of the figure shows the ratio of the nlo correction over the lo result and displays more clearly the higher order effects .
the correction changes from @xmath52  % for small @xmath10 to @xmath53  % for large @xmath10 . for @xmath54
gev we get an even larger correction due to collinear configurations .     for @xmath43 , @xmath44 , and @xmath45 with @xmath46  gev@xmath47 .
the lower part shows the corresponding k - factor @xmath55 of the differential cross sections at nlo and lo.,width=302 ]    the theoretical uncertainty of the predictions of a fixed order calculation is mainly given by the renormalization and fragmentation scale dependence .
the theoretical uncertainties due to the parton densities and fragmentation functions are not considered here .
another source of uncertainty lies in the factorization theorem itself since it predicts correct results with an error of @xmath56 which may be become large for very low @xmath10 .
we check the scale dependence of the cross section for the kinematics already depicted in context with fig .
[ diffcross ] but with a fixed transversal hadron momentum of @xmath57  gev . in fig .
[ scaledependence ] the differential cross section is shown where we vary the common scale @xmath58 over two orders of magnitude with respect to the reference scale @xmath59 .
the scale dependence varies in this range about + 33  % ( + 27  % ) for very low scales to -20  % ( -18  % ) for rather high scales at lo ( nlo ) .
thus there is only an unexpectedly slight scale dependence reduction at nlo compared to lo .    , @xmath44 , and @xmath45 with @xmath46  gev@xmath47 and a fixed @xmath57  gev .
the scale @xmath60 is varied with respect to the reference scale @xmath61  gev@xmath47.,width=302 ]
the calculation of single hadron production in deep inelastic scattering to @xmath11 was presented . owing to the factorization theorem this calculation
was performed as a convolution of the hard scattering process with universal and process independent parton densities and fragmentation functions .
this allows for meaningful predictions , directly comparable to experimental data .    results were shown for specific kinematic values .
a correction of the order of several percent is predicted depending on the transversal momentum of the observed hadron .
the scale dependence of the correction gives slightly reduced theoretical uncertainties compared to lo predictions .
varying the renormalization and factorization scales over two orders of magnitude yield theoretical uncertainties of about @xmath62  % at nlo .
an extensive study of nlo effects will be published elsewhere which in particular will include a detailed comparison of predictions with experimental data from the hera experiments h1 and zeus .
just after finishing the calculation , a preprint appeared on nlo calculations of hadron production with non - vanishing transversal momentum @xcite . in this paper matrix elements
were adopted from the disent program package and the phase space slicing method was applied to handle singularities .
thanks to gustav  kramer and bernd  a.  kniehl for proposing this project and collaboration .
we are very grateful to michael klasen and dominik stckinger for many helpful discussions on the calculation .
we want to thank michael spira for his advises to apply the subtraction formalism in the rather involved case of _ identified _ partons . for comments on the manuscript
we thank ingo schienbein .
j.  c.  collins , d.  e.  soper and g.  sterman , adv .
high energy phys .   * 5 * ( 1988 ) 1 . c.  adloff
_ et al . _
[ h1 collaboration ] , phys .
b * 462 * ( 1999 ) 440 m.  derrick _ et al .
_ [ zeus collaboration ] , z.  phys .
c * 70 * ( 1996 ) 1 a.  mendez , nucl .
b * 145 * ( 1978 ) 199 .
p.  bttner , `` inklusive meson - produktion in tief - inealastischer + elektron - proton - streuung in nchstfhrender ordnung der qcd , '' desy - thesis 1999 - 004 .
s.  catani and m.  h.  seymour , nucl .
b * 485 * ( 1997 ) 291 ( erratum - ibid .
b * 510 * ( 1997 ) 503 ) d.  graudenz , phys .
d * 49 * ( 1994 ) 3291 j.  a.  m.  vermaseren , arxiv : math - ph/0010025 . g.  passarino and m.  j.  g.  veltman , nucl .
b * 160 * ( 1979 ) 151 .
h.  l.  lai _ et al . _
[ cteq collaboration ] , eur .
j.  c * 12 * ( 2000 ) 375 b.  a.  kniehl , g.  kramer and b.  potter , nucl .
b * 582 * ( 2000 ) 514 k.  hagiwara _ et al .
_ [ particle data group collaboration ] , phys .  rev .
d * 66 * ( 2002 ) 010001 .
p.  aurenche , r.  basu , m.  fontannaz and r.  m.  godbole , arxiv : hep - ph/0312359 .
the nlo - qcd correction to single hadron production in deep inelastic scattering is calculated . 
 we require the final state meson to carry a non - vanishing transversal momentum , thus being sensitive to perturbative qcd effects . 
 factorization allows us to convolute the hard scattering process with parton densities and fragmentation functions . 
 the predictions are directly comparable to experimental results at the hera collider at desy . 
 the results are sensitive to the gluon density in the proton and allow us to test universality of fragmentation functions .