the phenomenological interest of prompt - photon production in fixed - target experiments  @xcite resides mainly in its use as a gluon probe in structure - function studies .
prompt - photon production is historically our main source of information on the gluon parton density at large  @xmath4 ( e.g. @xmath5 )  @xcite , a region which has very little influence on the evolution of the deep - inelastic - scattering structure functions .
this same region is relevant for hadron colliders in production phenomena at very large transverse momenta , and thus its understanding is crucial in order to disentangle possible new physics signals from the qcd background .    for example , a particularly interesting problem has emerged in the past few years in the production of large - transverse - energy ( ) jets at the tevatron .
an excess over the qcd prediction has been reported by the cdf collaboration  @xcite , for jets with @xmath6  gev . while the excess has not been confirmed by the d@xmath7 data  @xcite ,
it is of interest to study the uncertainty in the high-  tail of the jet distribution due to the gluon density systematics , to see whether there is room for deviations as large as those detected by cdf .
for example , a suitable modification of the gluon density at large @xmath4 has been proposed ( cteq4hj  @xcite ) , which is consistent with the excess observed by cdf .
the study of the recent e706 prompt - photon data  @xcite , however , suggests that a consistent fit of the large-@xmath8 ( @xmath9 ) rate is incompatible with the cteq4hj gluon density  @xcite .
moreover , both jet cross sections and direct - photon cross sections at high transverse energy are affected by soft - gluon effects .
these effects should be understood in both cases in order to be able to claim a discrepancy with qcd predictions .
in particular , these effects can be very important in the direct - photon case , since the typical  values probed are much smaller than in the case of jet production at the tevatron and , therefore , the size of the running coupling @xmath10 at the relevant scales is bigger .
comparisons between theory and prompt - photon experimental results have been carried out recently in refs .
the recent e706 data  @xcite seem to differ most from the next - to - leading order calculation , over the whole @xmath8 range . in refs .
@xcite , an attempt is made to fit the e706 data by introducing an intrinsic transverse momentum of the incoming partons with @xmath11 .
the precise details of how the intrinsic @xmath12 is incorporated in the calculations , however , can significantly affect the impact of these corrections , as shown by the large variations reported in ref .
@xcite .
the use of an intrinsic - transverse - momentum model is sometimes motivated as a way of estimating the effects of soft - gluon emission .
the most prominent effect of soft - gluon emission in drell - yan pair production is the generation of the characteristic transverse - momentum spectrum of the lepton pair .
this can be modeled with an appropriate intrinsic transverse momentum of the incoming partons . as a matter of fact
, the formalism for soft - gluon resummation in drell - yan pair production can be shown to merge , at very small transverse momenta , into some non - perturbative intrinsic transverse momentum of the partons inside the hadron  @xcite .
while this approach is not unreasonable when one considers the transverse momentum of the produced _ pair _ , it can however lead to inconsistencies for the problem of _ single_-photon production .
in fact , for example , it is quite clear that the photon @xmath8 spectrum at large @xmath8 explores the kinematic region of @xmath13 in the parton densities , which is certainly not the case for the transverse - momentum distribution of a drell - yan pair .
thus , as of now , a method for the inclusion of non - perturbative effects in the resummed formulae for the high-@xmath8 limit of the inclusive photon cross section is not available .
furthermore , in the opposite limit of small @xmath8 , it is the multiple emission of hard ( rather than soft ) gluons that leads to a sizeable perturbative broadening of the transverse momenta of the incoming partons  @xcite .    in this work ,
we consider the effect of soft - gluon resummation in prompt - photon production near the threshold limit , that is to say for @xmath14 .
the theoretical evaluation of these effects , at the next - to - leading logarithmic accuracy , has been carried out independently in refs .
@xcite and  @xcite .
we shall review in the next section the necessary formalism , using the language of ref .
@xcite . in the rest of the paper
, we will present its phenomenological applications , and we will thus discuss its numerical implementation as well as its impact on physical cross sections .    as is well known , prompt - photon production takes place both by hard - photon emission from initial- or final - state quarks ( direct component ) , and by collinear radiation from final - state partons .
this last mechanism is not fully calculable in perturbation theory and , in fact , it depends upon the photon fragmentation function .
because of the large suppression of the fragmentation function at large momentum fractions @xmath15 , it is usually believed that this contribution becomes irrelevant when @xmath16 increases .
contrary to common wisdom , we shall instead show that the very - large-@xmath8 behaviour of the direct and of the fragmentation production processes is the same if the incoming hadrons do not contain valence antiquarks , as in the case of @xmath17 collisions . under these circumstances
, resummation should therefore be performed for the fragmentation emission too .
we will show , however , that in the cases of practical interest the corrections due to the fragmentation processes are small , and we shall limit our considerations to the hard - photon part .
the plan of the paper is as follows . in section  [ secnotat ]
we review the formalism for the resummation of threshold effects , and the main formulas valid for the specific case of prompt - photon production .
there we also recall the main issues related to the inversion of the resummed expressions from mellin space back to the physical @xmath4 space . in section  [ results ]
we study numerically the impact of the resummation corrections .
we explore the effects both at the parton and hadron level , considering kinematical configurations and distributions of phenomenological relevance for current experiments .
in particular , we concentrate on the study of the size of the resummation corrections , and of the residual dependence on the choice of renormalization and factorization scales .
section  [ datacomp ] contains a comparison between our results and the data from some recent experiments .
this does not want to be a comprehensive phenomenological study , but a preliminary analysis of the impact of our results on the comparison of theory and data .
our conclusions , and the outlook for future progress , are given in section  [ conclusions ] .
an appendix collects some details of the resummation formulas .
we consider the inclusive production of a single prompt photon in hadron collisions : [ procgamma ] h_1(p_1 ) + h_2(p_2 ) ( p ) + x .
the colliding hadrons @xmath18 and @xmath19 carry momenta @xmath20 and @xmath21 , respectively .
we parametrize the momenta in terms of light - cone coordinates : p^= ( p^+ , p_t , p^- ) , p^ ( p^0 p^3 ) . in their centre - of - mass frame , using massless kinematics , the momenta of the colliding hadrons have the following light - cone coordinates [ mom ] p_1^= ( 1,0,0 ) , p_2^= ( 0,0,1 ) , where @xmath22 is the centre - of - mass energy squared . the photon momentum @xmath23 is thus parametrized as [ momp ] p^= ( e^y , e_t , e^-y ) , where @xmath24 and @xmath25 are the transverse energy and the rapidity , respectively .
we also introduce the customary scaling variable @xmath8 @xmath26 : [ xt ] x_t = .    in the present paper
we are mostly interested in the prompt - photon production cross section integrated over @xmath25 at fixed @xmath24 . according to perturbative qcd
, the cross section is given by the following factorization formula [ 1pxsgamma ] & & = _ a , b _ 0 ^ 1 dx_1 f_a / h_1(x_1,_f^2 ) _
0 ^ 1 dx_2 f_b / h_2(x_2,_f^2 ) + & & _ 0 ^ 1 dx \ { ( x - ) _
ab(x , ( ^2 ) ; e_t^2 , ^2 , _ f^2 , _
f^2 ) . + & & + . _ c _ 0 ^ 1 dz z^2 d_c/(z,_f^2 ) ( x - ) _ abc(x , ( ^2 ) ; e_t^2 , ^2 , _
f^2 ) } . where @xmath27 denotes the parton indices @xmath28 , and
@xmath29 and @xmath30 are the parton densities of the colliding hadrons , evaluated at the factorization scale @xmath31 .
the first and the second term in the curly bracket on the right - hand side of eq .
( [ 1pxsgamma ] ) represent the _ direct _ and the _ fragmentation _ component of the cross section , respectively .
the fragmentation component involves the parton fragmentation function @xmath32 of the observed photon at the factorization scale @xmath33 , which , in general , differs from the scale @xmath31 of the parton densities .
the _ rescaled _ ( @xmath34 ) . ]
partonic cross sections @xmath35 and @xmath36 in eq .
( [ 1pxsgamma ] ) are computable in qcd perturbation theory as power series expansions in the running coupling @xmath37 , @xmath38 being the renormalization scale in the @xmath39 renormalization scheme : [ pxsg ] _
ab(x , ( ^2 ) ; e_t^2 , ^2 , _
f^2 ) = ( ^2 ) , [ pxsc ] _ abc(x , ( ^2 ) ; e_t^2 , ^2 , _
f^2 ) = ^2(^2 ) . note that the ratio between the direct and the fragmentation terms in eqs .
( [ pxsg ] ) and ( [ pxsc ] ) is of the order of @xmath40 , where @xmath41 is the fine structure constant .
this ratio is compensated by the photon - fragmentation function @xmath42 , which ( at least formally ) is of the order of @xmath40 , so that direct and fragmentation components equally contribute to eq .
( [ 1pxsgamma ] ) .    throughout the paper we always use parton densities and parton fragmentation functions as defined in the @xmath39 factorization scheme .
in general , we consider different values for the renormalization and factorization scales @xmath38 , @xmath31 , @xmath33 , although we always assume that all of them are of the order of the photon transverse energy @xmath24 .
the lo terms @xmath43 in eq .
( [ pxsg ] ) are due to the following parton - scattering subprocesses at the tree - level [ loproc ] q + |q g + , q + g q + , + g + . using our normalization , the two independent ( non - vanishing ) partonic cross sections for the direct component are : [ siqqgamma ] _ q|q g^(0)(x ) = e_q^2 ( 2 - x^2 ) [ siqggamma ] _
qg q^(0)(x ) = _
|qg ^(0)(x ) = e_q^2 ( 1 + ) , where @xmath44 is the quark electric charge .
note that , having integrated over the photon pseudorapidity , the expressions  ( [ siqqgamma ] ) and  ( [ siqggamma ] ) are even functions of the photon transverse energy @xmath24 , i.e.  they depend on @xmath45 rather than on @xmath4 .
the nlo terms @xmath46 in eq .
( [ pxsg ] ) were first computed in ref .
@xcite .
the partonic contributions @xmath36 to the fragmentation component of the cross section are exactly equal to those of the single - hadron inclusive distribution .
note that , unlike in the case of the direct component , _ all _ the parton - parton scattering subprocesses @xmath47 ( i.e.  including @xmath48 ) contribute to the fragmentation component already at lo .
the explicit calculation of @xmath36 up to nlo was performed in ref .
@xcite .
the behaviour of the lo and nlo perturbative contributions to the direct component of the prompt - photon cross section is shown in fig .
[ fig : qalonlo ] .
the lo terms @xmath49 and @xmath50 are both singular when @xmath51 : [ lobeh ] _ abd ^(0)(x ) ~ , ( x 1 ) , and they both vanish in the high - energy limit @xmath52 . the integrable singularity in eq .
( [ lobeh ] ) is a typical phase - space effect , while the vanishing behaviour at small-@xmath4 is due to the dominance of fermion ( i.e.  spin @xmath53 ) exchange in the @xmath54-channel .
two new dynamical features appear at nlo . near the threshold region @xmath51 ,
the nlo contributions are double - logarithmically enhanced , [ nlobeh ] ^(1)(x ) ~ ^(0)(x ) ^2 ( 1-x ) , ( x 1 ) , because the radiation of soft and , possibly , collinear partons is strongly inhibited by the kinematics . in the high - energy limit @xmath52 , the partonic cross sections @xmath55 approach constant values  @xcite : this regge plateau follows from the fact that at nlo single - gluon ( i.e.  spin 1 ) @xmath54-channel exchange affects all the partonic subprocesses .
the behaviour of the partonic contributions @xmath56 in the remaining intermediate region of @xmath4 has no straightforward physical interpretation ( e.g.  @xmath57 even becomes negative ) because it strongly depends on the scale - dependent corrections already subtracted in the definition of the parton densities and parton fragmentation functions .
higher - order perturbative qcd corrections in the small-@xmath8 regime can systematically be computed by using the @xmath58-factorization approach  @xcite , which consistently takes into account the _ perturbative _ broadening of the transverse momenta of the incoming partons .
we are interested in this work in the behaviour of the qcd corrections near the partonic - threshold region @xmath51 , i.e.  when the transverse energy @xmath24 of the photon approaches the partonic centre - of - mass energy @xmath59 . in this region , the singularities in eqs .
( [ lobeh ] , [ nlobeh ] ) are enhanced by double - logarithmic corrections due to soft - gluon radiation and the higher - order cross section contributions in eqs .
( [ pxsg ] , [ pxsc ] ) behave as [ hobeh ] ^(n)(x ) ~ ^(0)(x ) .
resummation of these soft - gluon effects to all orders in perturbation theory can be important to improve the reliability of the qcd predictions .
the resummation program of soft - gluon contributions has to be carried out  @xcite in the mellin - transform space , or @xmath60-space .
working in @xmath60-space , we can disentangle the soft - gluon effects in the parton densities from those in the partonic cross section and we can straightforwardly implement and factorize the kinematic constraints of energy and longitudinal - momentum conservation .    the latter point is particularly relevant for soft - gluon resummation in hadron collisions  @xcite .
indeed , all - order soft - momentum recoil can not exactly be taken into account by directly working in @xmath4-space and the ensuing kinematics approximation leads to ( same - sign ) factorially growing coefficients .
this implies  @xcite that no resummed logarithmic hierarchy can consistently be defined in @xmath4-space ( the classes of leading logs @xmath61 , next - to - leading logs @xmath62 and so forth in eq .
( [ hobeh ] ) are not separately summable , because they lead to divergent and not integrable contributions at @xmath63 ) . on the contrary ,
no kinematics approximation ( in the soft limit ) is required in @xmath60-space and the corresponding logarithmic hierarchy of @xmath64-contributions is systematically well defined .
to work in @xmath60-space , it is convenient to consider the mellin transform @xmath65 of the dimensionless hadronic distribution @xmath66 .
the @xmath60-moments with respect to @xmath67 and at fixed @xmath24 are thus defined as follows : [ shn ] _ , n(e_t ) _ 0 ^ 1 dx_t^2 ( x_t^2)^n-1 e_t^3 .
in @xmath60-moment space , eq .
( [ 1pxsgamma ] ) takes a simple factorized form [ 1pxsngamma ] _ ,
n(e_t ) & = & _ a , b f_a / h_1,n+1(_f^2 ) f_b / h_2,n+1(_f^2 ) + & & \ { _ ab , n((^2 ) ; e_t^2 , ^2 , _
f^2 ) + & + & .
_ c _ abc , n((^2 ) ; e_t^2 , ^2 , _
f^2 ) d_c/ , 2n+3(_f^2 ) } , where we have introduced the customary @xmath60-moments @xmath68 and @xmath69 of the parton densities and parton fragmentation functions : [ pdm ] f_a / h , n(^2 ) & & _ 0 ^ 1 dx x^n-1 f_a / h(x,^2 ) , + [ pffm ] d_a/ ,
n(^2 ) & & _ 0 ^ 1 dz z^n-1 d_a/(z,^2 ) . note that the @xmath60-moments of the partonic cross sections in eq .
( [ 1pxsngamma ] ) are again defined with respect to @xmath67 : [ ndef ] _ ab , n((^2 ) ; e_t^2 , ^2 , _
f^2 ) _ 0 ^ 1 dx^2 ( x^2)^n-1 _ ab(x,(^2 ) ; e_t^2 , ^2 , _
f^2 ) . the explicit expressions of the @xmath60-moments @xmath70 , @xmath71 of the lo contributions in eqs .
( [ siqqgamma ] ) and  ( [ siqggamma ] ) were obtained in ref .
@xcite and are recalled in appendix  [ appa ] .
note also the pattern of moment indices in the various factors of eq .
( [ 1pxsngamma ] ) , that is , @xmath72 for the parton densities and @xmath73 for the parton fragmentation functions .
this non - trivial pattern follows from the conservation of the longitudinal and transverse momenta .
the threshold region @xmath74 corresponds to the limit @xmath75 in @xmath60-moment space . in this limit ,
the soft - gluon corrections  ( [ hobeh ] ) to the higher - order contributions of the partonic cross sections become [ hobehn ] ^(n)_n ~ ^(0)_n .
the resummation of the soft - gluon logarithmic corrections to all orders in perturbation theory has been considered in refs .
@xcite and  @xcite . in the following section we recall the main results .      in ref .
@xcite soft - gluon resummation has been performed in detail for the various partonic channels that contribute to the _ direct _ component of the prompt - photon cross section @xmath65 in eq .
( [ 1pxsngamma ] ) .    we discuss first the large-@xmath60 behaviour of the partonic cross sections @xmath76 for the partonic channels @xmath77 that start to contribute at lo .
these cross sections can be written as [ sigreslo ] _ ab , n = _ ab , n^(res ) , ab = q|q , qg , |qg , where @xmath78 denotes terms that contribute beyond lo and are furthermore suppressed by a relative factor @xmath79 at large @xmath60 .
the logarithmically - enhanced soft - gluon corrections are included in the resummed expressions @xmath80 and can be factorized with respect to the corresponding lo cross sections @xmath81 . the _ all - order _ resummation formulae are _
q|q , n^(res)((^2 ) ; e_t^2 , ^2 , _ f^2 , _
f^2 ) & = & ( ^2 ) _ q|qg , n^(0 ) c_q|q ( ( ^2),q^2/^2;q^2/_f^2 ) + [ gammaresqq ] & & _ n+1^q|q g ( ( ^2),q^2/^2;q^2/_f^2 ) , +   + _ qg , n^(res)((^2 ) ; e_t^2 , ^2 , _ f^2 , _
f^2 ) & = & ( ^2 ) _ qgq , n^(0 ) c_qg ( ( ^2),q^2/^2;q^2/_f^2 ) + [ gammaresqg ] & & _ n+1^qg q ( ( ^2),q^2/^2;q^2/_f^2 ) , +   + [ gammaresqbarg ] _
|qg , n^(res)((^2 ) ; e_t^2 , ^2 , _ f^2 , _
f^2 ) & = & _ qg , n^(res)((^2 ) ; e_t^2 , ^2 , _ f^2 , _
f^2 ) , where [ etscale ] q^2 = 2 e_t^2 .    the functions @xmath82 in eqs .  ( [ gammaresqq ] , [ gammaresqg ] )
do not depend on @xmath60 .
thus , the @xmath64-dependence of the resummed cross sections is entirely embodied by the radiative factors @xmath83 .
they depend on the flavour of the qcd partons @xmath84 involved in the lo hard - scattering subprocess @xmath85 and can be expressed in an exponential form : & & _ n^ab d ( ( ^2 ) , ; ) = \ { n g_ab^(1)(b_0(^2)n ) + & & + g_ab^(2)(b_0(^2)n , q^2/^2;q^2/_f^2 ) [ deltanll ] + o((n)^k ) } , where @xmath86 is the first coefficient of the qcd @xmath87-function [ betafirst ] b_0 = . note that the functions @xmath88 , @xmath89 and so forth in the exponent do not depend separately on @xmath10 and @xmath64 .
they are functions of the expansion variable @xmath90 and vanish when @xmath91 .
this means that the exponentiation structure in eq .
( [ deltanll ] ) is not trivial and , in particular , that all the double logarithmic ( dl ) terms @xmath92 in eq .
( [ hobehn ] ) are taken into account by simply exponentiating the lowest - order contribution @xmath93 .
the exponentiation in eq .
( [ deltanll ] ) defines an improved perturbative expansion in the threshold region .
the function @xmath94 resums all the _ leading _ logarithmic ( ll ) contributions @xmath95 in the exponent , @xmath89 contains the _ next - to - leading _ logarithmic ( nll ) terms @xmath96 , and so forth .
once the functions @xmath97 have been computed , we have a systematic perturbative treatment of the region of @xmath60 where @xmath98 , which is much larger than the domain @xmath99 where the fixed - order calculation in @xmath10 is reliable .    the ll and nll functions @xmath88 and @xmath89 in eq .
( [ deltanll ] ) have been explicitly computed in ref .
the ll functions @xmath88 are different for the @xmath100 and @xmath101 partonic channels of eqs .
( [ gammaresqq ] ) and  ( [ gammaresqg ] ) but they can be expressed in terms of parton colour factors and a single ( parton - independent ) function @xmath102 : g_q|q^(1 ) ( ) & = & ( 2c_f - c_a ) h^(1 ) ( ) + c_a h^(1)(/2 ) , + [ g1fun ] g_qg^(1 ) ( ) & = & c_a h^(1 ) ( ) + c_f h^(1)(/2 ) , with [ htll ] h^(1 ) ( ) = .
the explicit expressions of the nll functions @xmath89 are recalled in appendix  [ appa ] .
note that the ll functions @xmath88 do not depend on the factorization scale @xmath31 .
this dependence starts to appear only in the nll functions @xmath89 .
note also the mismatch between the moment index of the radiative factor and that of @xmath81 in eqs .
( [ gammaresqq ] , [ gammaresqg ] ) : the former depends on @xmath103 , like the parton densities in eq .
( [ 1pxsngamma ] ) .
the explicit @xmath31-dependence of @xmath89 exactly matches the scale dependence of the parton densities at large values of @xmath60 .
thus , when ( and only when ) nll resummation is included , we can expect  @xcite better stabilization of the calculation of the cross section at large @xmath8 with respect to variations of the factorization scale @xmath31 ( see sec .
[ hlresults ] ) .
the functions @xmath82 in eqs .  ( [ gammaresqq ] , [ gammaresqg ] ) contain all the terms that are constant in the large-@xmath60 limit .
they are produced by hard _
virtual _ contributions and by subdominant ( non - logarithmic ) soft corrections to the lo hard - scattering subprocesses .
these functions are computable as power series expansions in @xmath10 c_ab((^2),q^2/^2;q^2/_f^2 ) & = & 1 + _ n=1^+ ( ) ^n c_ab^(n)(q^2/^2;q^2/_f^2 ) + [ cgamma ] & = & 1 + c_ab^(1)(q^2/^2;q^2/_f^2 ) + o(^2 ) . at present , we know only the first - order constant coefficients @xmath104 and @xmath105 in eqs .
( [ cgamma ] , [ gammaresqq ] , [ gammaresqg ] ) .
these coefficients can be extracted  @xcite from the complete nlo analytic results of refs .
their values are recalled in appendix  [ appa ] .
the inclusion of the @xmath60-independent function @xmath82 in the resummed formulae does not affect the shape of the cross section near threshold , but improves the soft - gluon resummation by fixing the overall ( perturbative ) normalization of the logarithmic radiative factor .
we can explicitly show  @xcite the theoretical improvement that is obtained by combining the nll radiative factor with the first - order coefficient @xmath106 . expanding the resummation formulae  ( [ gammaresqq ] , [ gammaresqg ] ) in towers of logarithmic contributions as in eq .
( [ hobehn ] ) , we have & & _
n^(res)(;e_t^2,^2,_f^2 ) = _ n^(0 ) \ { 1 + _ n=1^ ^n } , [ sigtow ] where @xmath107 .
the dominant and next - to - dominant coefficients @xmath108 and @xmath109 are controlled by evaluating the radiative factor to nll accuracy . when the nll radiative factor is supplemented with the coefficient @xmath110
, we can correctly control also the coefficients @xmath111 .
in particular , we can predict  @xcite the large-@xmath60 behaviour of the next - to - next - to - leading order ( nnlo ) cross sections @xmath112 in eq .
( [ pxsg ] ) up to @xmath113 .
note also that the coefficients @xmath108 are scale independent and the coefficients @xmath109 depend on the sole factorization scale @xmath31 . in the tower expansion  ( [ sigtow ] ) ,
the first terms that explicitly depend on the renormalization scale @xmath38 ( and on @xmath31 , as well ) are those controlled by @xmath111 .
their dependence on @xmath38 is obtained by combining that of @xmath114 with that of the radiative factor at nll order .
the inclusion of the first - order constant coefficient @xmath110 thus theoretically stabilizes the resummed partonic cross section at large @xmath8 with respect to variations of the renormalization scale .
this scale dependence is numerically studied in sec .
[ hlresults ] .
so far we have only considered the near - threshold behaviour of the partonic cross sections @xmath115 , @xmath116 , @xmath117 in eq .
( [ sigreslo ] ) .
the behaviour of other partonic channels @xmath118 that contribute to the direct component of the prompt - photon cross section was discussed in ref .
it turns out that the partonic channel @xmath119 enters the resummed cross section only at next - to - next - to - leading logarithmic ( nnll ) accuracy and that all the other channels are relatively suppressed in the same way as the correction @xmath78 on the right - hand side of eq .
( [ sigreslo ] ) .
since we are interested in explicitly performing soft - gluon resummation up to nll order , we can limit ourselves to considering the resummed expressions in eqs .
( [ sigreslo])([gammaresqbarg ] ) .
detailed numerical studies of the resummed cross sections are presented in sec .
[ results ] .
however , from the analytical results reviewed in this section , we may already anticipate that soft - gluon resummation increases the perturbative qcd predictions in the large-@xmath8 region .
this conclusion can be argued by a simplified treatment within the dl approximation . to dl accuracy ,
the exponent of the radiative factors in eq .
( [ deltanll ] ) has to be expanded to its first order in @xmath10 , and we obtain [ dlresqg  ] & & \ { ^2 n } = \ { ( c_f + 2c_a ) ^2 n } > 1 , + &  & + [ dlresqq  ] & & \ { ^2 n } = \ { ( 4c_f - c_a ) ^2
n } > 1 , + &  & + [ dlrat ] & & \ { 3(c_a - c_f ) ^2 n } > . for the sake of completeness , in the square bracket on the right - hand side of eqs .
( [ dlresqg  ] ) and  ( [ dlresqq  ] ) we have explicitly separated the positive contributions coming from the initial - state partons and the negative contribution from the final - state recoil . from these equations
we see that the resummed partonic cross sections @xmath120 and @xmath121 are both enhanced with respect to their lo approximations @xmath70 , @xmath71 . moreover , the enhancement in the @xmath101 partonic channel is larger than that in the @xmath100 channel .
we refer the reader to ref .
@xcite for a discussion on the physical origin of this behaviour .
we can now comment on the large-@xmath24 behaviour of the fragmentation component of the prompt - photon cross section , by comparing the direct and fragmentation contributions in eq .
( [ 1pxsngamma ] ) .    the partonic cross sections @xmath76 and
@xmath122 have the same large-@xmath60 behaviour , but , owing to the hard ( although collinear ) emission always involved in any splitting process @xmath123 , the photon fragmentation function @xmath124 is of the order of @xmath125 .
therefore , in the curly bracket on the right - hand side of eq .
( [ 1pxsngamma ] ) the fragmentation component is formally suppressed by a factor of @xmath125 with respect to the direct component .
this suppression is consistent with the fact that the resummed partonic cross sections for the direct processes ( see the right - hand side of eqs .
( [ gammaresqq ] ) and  ( [ gammaresqg ] ) ) turn out to be independent of the photon fragmentation scale @xmath33 .
this argument shows that , in many cases , the fragmentation contributions are subdominant near threshold and , thus , they can be neglected in resummed calculations at large @xmath8 .
the caveat ` in many cases ' in the above conclusion regards the fact that the argument applies to the partonic contributions in the curly bracket of eq .
( [ 1pxsngamma ] ) . in other words
, the argument assumes that all the different initial - state partonic channels @xmath126 give comparable contributions to the hadronic cross section .
this is not always true once the effect of the parton densities is included .
a relevant exception is indeed the case of prompt - photon production in proton - nucleon collisions .
owing to the low antiquark content of the colliding hadrons , the hadronic cross section is mostly due to the partonic channels @xmath127 and @xmath128 : [ xsngammapn ] _ ,
n(e_t ) ~ _ , n^qg(e_t ) + _ , n^qq(e_t ) , ( pn ) .
as for the @xmath101 initial - state contribution @xmath129 , we can use the above argument to conclude that its direct component dominates at large @xmath24 . setting all the scales equal to @xmath24 , for the sake of simplicity
, we can write : [ xsngammapnqg ] _ ,
n^qg(e_t ) ~ _ , n^qg ( dir)(e_t ) ~f_q , n+1(e_t^2 ) f_g , n+1(e_t^2 ) _ qg , n((e_t^2 ) ) .
however , in the case of the @xmath130 initial state , the direct component enters only at nlo and , thus , the cross section is dominated by the fragmentation part and , in particular , by photon fragmentation from a final - state quark of the lo scattering subprocess @xmath131 .
we can write : [ xsngammapnqq ] _ ,
n^qq(e_t ) ~ _ , n^qq ( frag)(e_t ) ~f_q , n+1(e_t^2 ) f_q , n+1(e_t^2 ) _ qqq , n((e_t^2 ) ) d_q/ , 2n+3(e_t^2 ) .    taking the ratio of the two initial - state contributions and replacing @xmath132 and @xmath133 by their lo contributions @xmath134 in eqs .  ( [ pxsg ] , [ pxsc ] )
, we obtain [ qgqqratio ] ~ d_q/ , 2n+3(e_t^2 ) .
the factor @xmath135 on the right - hand side is compensated by the behaviour of the photon fragmentation function @xmath136 . in the large-@xmath60 limit ,
the ratio of the lo partonic contributions @xmath134 is constant and , thus , the fragmentation function produces an @xmath79-suppression factor .
nonetheless , this suppression can be balanced by the parton density contribution @xmath137 since , at large @xmath4 , the gluon density is typically softer than the quark density . as a matter of fact , using the altarelli - parisi evolution equation at lo and under reasonable assumptions on the large-@xmath4 behaviour of the parton densities at the initial evolution scale , it is easy to show that we have the following asymptotic behaviour at very large values of @xmath60 and of the evolution scale : [ eq : largenfg ] f_g , n ~ f_q , n , + [ eq : largendg ] d_q/ , n ~ . combining these results with eq .
( [ qgqqratio ] ) , in the large-@xmath60 limit we therefore get [ eq : fraglargen ] .
this discussion shows that , in the case of large-@xmath8 prompt - photon production in @xmath17 collisions , the contribution of the fragmentation component to the hadronic cross section can become comparable to that of the direct component .
although on the right - hand side of eq .
( [ qgqqratio ] ) we have approximated the partonic contribution @xmath138 by its lo expansion , the inclusion of higher - order terms and , in particular , of resummation effects does not substantially modify the conclusion .    performing soft - gluon resummation in the partonic cross sections @xmath122 of the fragmentation component , we can write an expression that is analogous to eq .
( [ sigreslo ] ) : [ sigreslof ] _ abc ,
n = _ abc , n^(res ) . limiting our treatment to the ll accuracy
, the resummed cross section is given by  @xcite [ gammaresf ] _ abc , n^(res)((^2 ) ; e_t^2 , ^2 , _
f^2 ) ^2(^2 ) _ abdc , n^(0 ) _ n+1^ab
dc((^2),q^2/^2;q^2/_f^2,q^2/_f^2 ) . the radiative factor is [ deltanllf ] _
n^ab d c((^2 ) , ; , ) = \ { n g_abdc^(1)(b_0(^2)n ) + o((n)^k ) } , where the ll function @xmath139 is analogous to those in eq .
( [ g1fun ] ) and it can be expressed in terms of the colour charges @xmath140 ( @xmath141 , if the parton @xmath142 is a quark , and @xmath143 , if @xmath142 is a gluon ) of the partons involved in the lo hard - scattering subprocess : g_abdc^(1 ) ( ) = ( c_a + c_b + c_c - c_d ) h^(1 ) ( ) + c_d h^(1)(/2 ) .
in particular , for the @xmath144 channel we have [ g1qqqq ] g_qqqq^(1 ) ( ) = 2 c_f h^(1 ) ( ) + c_f h^(1)(/2 ) , which is very similar to @xmath145 in eq .
( [ g1fun ] ) because @xmath146 .
more precisely , since @xmath147 , @xmath148 is slightly smaller than @xmath145 as long as they are evaluated in the perturbative region @xmath149 .
we can now come back to the effect of the fragmentation component in @xmath17 collisions . using the resummed partonic cross sections in eqs .
( [ gammaresqg ] , [ gammaresf ] ) rather than their lo approximations , the right - hand side of eq .
( [ qgqqratio ] ) has to be multiplied by an additional contribution , as given by the ratio of the corresponding radiative factors , namely [ deltacorr ] .
because of the relation @xmath150 between the colour charges , this factor does not sizeably differ from unity , as it can be argued by its dl approximation on the right - hand side of eq .
( [ deltacorr ] ) ( see also the comment below eq .
( [ g1qqqq ] ) ) .
we have thus shown that , at least at the ll level , soft - gluon resummation does not enhance the relative importance of the fragmentation component for prompt - photon production at large @xmath8 .
the importance of the fragmentation component in @xmath17 collisions mainly depends on the detailed behaviour of the parton densities at large @xmath4 and on how large are the values of @xmath24 of interest .
this issue , as well as the impact of the nll corrections to the ll results obtained above , require further studies that will be presented in a future work . as for the present study , we limit ourselves to perform soft - gluon resummation in the direct component and we check that the fragmentation component does not sizeably contribute to the hadronic cross section in the actual experimental configurations investigated in the paper ( see sec .
[ datacomp ] ) .
we use soft - gluon resummation to nll accuracy at the parton level to introduce an improved prompt - photon cross section @xmath151 as follows _ , n^(res)(e_t ) & = & _ ab = q|q , qg , |qg f_a / h_1,n+1(_f^2 ) f_b
/ h_2,n+1(_f^2 ) + & & + [ ppcsres ] & + & _ , n^(nlo)(e_t ) , where @xmath152 is the prompt - photon hadronic cross section at nlo , @xmath80 is given in eqs .
( [ gammaresqq])([gammaresqbarg ] ) and @xmath153 represents its perturbative truncation at order @xmath154 ( i.e.  at nlo ) . thus , because of the subtraction in the square bracket on the right - hand side , eq .
( [ ppcsres ] ) exactly reproduces the nlo results and resums soft - gluon effects beyond @xmath155 to nll accuracy . in general , we evaluate @xmath80 using the nll expression  ( [ deltanll ] ) of the radiative factors and including the @xmath156 contribution  ( [ cgamma ] ) of the constant factors @xmath82 .
this defines our nlo+nll predictions .
the resummed formulae presented so far are given in @xmath60-moment space . to obtain cross sections in the physical @xmath8-space
( i.e.  as functions of the centre - of - mass energy ) , one has to perform the inverse mellin transformation : [ mppp ] e_t^3 = _ c_mp - i^c_mp+i dn x_t^-2n _ , n^(res)(e_t ) .    when the @xmath60-moments @xmath157 are evaluated at a fixed perturbative order in @xmath10 , they are analytic functions in a right half - plane of the complex variable @xmath60 . in this case , the constant @xmath158 that defines the integration contour in eq .
( [ mppp ] ) has to be chosen in this half - plane , that is , on the right of all the possible singularities of the @xmath60-moments .
an additional complication occurs when the @xmath60-moments are computed in resummed perturbation theory . in this case , since the resummed functions @xmath159 in eq .
( [ g1fun ] ) ( as well as the nll functions @xmath160 ) are singular at @xmath161 , the soft - gluon factors @xmath162 in eq .
( [ deltanll ] ) have cut singularities that start at the branch - point @xmath163 .
these singularities , which are related to the divergent behaviour of the perturbative running coupling @xmath10 near the landau pole , signal the onset of non - perturbative phenomena at very large values of @xmath60 or , equivalently , in the region very close to threshold .
the issue of how to deal with the landau singularity in soft - gluon resummation formulae for hadronic collisions was discussed in detail in ref .
@xcite . in the evaluation of the inverse mellin transformation  ( [ mppp ] ) we thus use the _ minimal prescription _ introduced in ref .
the constant @xmath158 is chosen in such a way that all singularities in the integrand are to the left of the integration contour , except for the landau singularity at @xmath164 , that should lie to the far right .
this prescription is consistent  @xcite with the perturbative content of the soft - gluon resummation formulae because it converges asymptotically to the perturbative series and it does not introduce ( unjustified ) power corrections of non - perturbative origin .
these corrections are certainly present in physical cross - sections , but their effect is not expected to be sizeable as long as @xmath24 is sufficiently perturbative and @xmath8 is sufficiently far from the hadronic threshold . obviously , approaching the essentially non - perturbative regime @xmath165  gev , @xmath74 , a physically motivated treatment of non - perturbative effects has to be introduced . in the following sections , we limit ourselves to presenting numerical and phenomenological results that do not include any non - perturbative correction
we present in this section some numerical results , to provide an illustration of the size of the effects considered and to show the improvements obtained with respect to scale variations after the inclusion of the nll corrections .
we start by discussing the resummation effects at the level of partonic cross sections .
the resummed partonic cross section can be obtained from eqs .
( [ ppcsres ] ) and ( [ mppp ] ) by assuming parton - density functions of the form @xmath166 , and hence @xmath167 for all complex values of @xmath60 .
we consider first the @xmath168 terms in the expansion of the resummed cross section , in order to estimate to which accuracy this reproduces the exact nlo results . in fig .
[ fig : qqnlo ] ( left ) we plot the function @xmath169 , defined in eq .
( [ pxsg ] ) , as a function of @xmath170 .
the exact @xmath168 result  @xcite is compared with three possible implementations of the resummation procedure , all equivalent at nll accuracy .
the first case ( short - dashed line ) corresponds to our default resummed prediction , as given by eq .
( [ gammaresqq ] ) . in the second case ( dot - dashed line )
, we set the constant @xmath171 introduced in eq .
( [ cgamma ] ) equal to 0 . in the third case
, we keep the contribution of the constant @xmath172 , but we modify it by a term suppressed by a factor of @xmath125 , in order to explore the possible effect of contributions of order @xmath125 which can not be taken into account by the soft - gluon resummation . as a constraint on the form of these corrections
, we must impose that no poles appear on the positive real axis in the @xmath60 plane ( these poles would logarithmically enhance the partonic cross section when @xmath173 ) .
we select a parametrization of the @xmath125 corrections that allows us to bracket the exact result at  : [ eq : cqashift ]
c_q |q ^(1 ) c_q |q ^(1 ) ( 1 + ) , [ eq : cqgshift ] c_q g^(1 ) c_qg^(1 ) ( 1- ) , with @xmath174 . in our applications
we shall consider the two cases with @xmath175 ( namely no correction to the @xmath176 term ) and @xmath177 as a way to establish the size of subleading threshold corrections beyond the nll order .    as one can see from fig .
[ fig : qqnlo ] , the inclusion of the finite term @xmath178 is essential to accurately reproduce shape and normalization of the exact  result not only near threshold , but below it as well .
the agreement deteriorates unavoidably for @xmath179 , as , here , terms subleading in @xmath125 become important .
analogous results for the @xmath101 channel are given in the right panel of fig .
[ fig : qqnlo ] .
note that in both cases the two choices @xmath175 and @xmath177 in eq .
( [ eq : cqgshift ] ) bracket the exact result over a large region of @xmath180 , and thus provide a good estimator of the subleading terms systematics .
the choice @xmath175 , furthermore , provides a very accurate description up to values of @xmath180 of the order of 1/10 .    the fully - resummed parton - level cross sections are shown in fig .  [ fig : qares10 ] for the @xmath181 and @xmath101 channels ( left and right panel , respectively ) , and for @xmath182  gev ( fig .
[ fig : qares100 ] collects the same results for @xmath183  gev ) . here and
in the following we shall define the resummed cross sections as in eq .
( [ ppcsres ] ) , that is , we substitute their @xmath168 terms with the exact nlo result , using the same choice of renormalization and factorization scales . in this way our results are exact up to ( and including ) @xmath168 , and include the nll resummation of terms of @xmath184 and higher .
we compare the fixed - order results ( dashed lines ) with the resummed results . for these
we provide both the @xmath175 and @xmath177 prescriptions .
note that , even at the level of resummed cross sections , the difference between the @xmath175 and @xmath177 results are rather small , in particular for the @xmath185 channel .      in this section
we present some results for the full hadronic cross sections .
the main points we intend to highlight are :    1 .   the size of the nll corrections , relative to the nlo contributions ; 2 .   the scale dependence at nll order .
our goal here is to explore the pure effects of resummation at higher orders .
therefore we shall neglect in this section all production channels which are not improved by the resummation corrections considered in this work .
this includes all processes which first appear at , such as @xmath186 and @xmath187 , as well as all contributions proportional to a parton @xmath188 fragmentation function , as discussed in the previous section .
these terms will however be included , at fixed nlo , in our comparison with experimental data , performed in the next section .    to be more specific ,
we list here the classes of diagrams included at nlo , in addition to the lo processes @xmath189 and @xmath190 ( the possible replacement of quarks with antiquarks in all these cases is understood ) : [ eq : procsin ] q g g , q|q q|q , q|q q|q , q|q gg . for the third set we only include the diagrams proportional to @xmath191 , since the part of the amplitude describing the photon emission from the final - state quarks ( which is by itself gauge invariant ) , can not be considered as a correction to any tree - level process .    as
a default set of parton densities we shall use the cteq4 m set described in ref .  @xcite .
for the purposes of the present study , no significant change is obtained if different sets are used .
figures  [ fig : deltaqa ] and  [ fig : deltatot ] present the ratios : [ deltadef ] , , , where , for simplicity , we indicated here with @xmath192 the differential distribution @xmath193 . for each channel
we present the results using both the @xmath175 ( solid lines ) and @xmath177 ( dashed lines ) prescriptions .
we also show the dependence on the choice of renormalization and factorization scales , which we take equal , and varying within the set @xmath194 . in this section
we shall always keep the fragmentation scale @xmath33 , necessary for the factorization of the singularities from final - state collinear photon emission , equal to .
note that the size of the resummation effects is larger for the larger scales , contrary to the behaviour of the scale dependence of the nlo cross section .
this suggests that the scale dependence of the resummed cross section will be reduced relative to that of the nlo results .
the scale dependence of the resummed cross section ( @xmath195 contributions ) , compared to the nlo one , is given in fig .
[ fig : scdepl0 ] , for the @xmath175 case .
the same result for @xmath177 is given in the left panel of fig .
[ fig : scdepl ] . we plot the distributions as a function of the beam energy ( @xmath196 ) for the fixed value of @xmath197  gev .
different values of @xmath196 , therefore , probe different ranges of , as indicated by the upper labels on the plots .
note the significant reduction in scale dependence , more marked in the @xmath175 case . to display the importance of the inclusion of the constant @xmath198 terms
, we show the same scale - dependence plot with @xmath199 in the right panel of fig .
[ fig : scdepl ] .
while the scale sensitivity is slightly worse than in the cases with @xmath200 , there is still an important improvement over the nlo behaviour .
similar results , for @xmath201  gev , are shown in figs .
[ fig : scdeph0 ] and  [ fig : scdeph ] . the general features of these distributions are similar to those of the plots for @xmath197  gev
. small violations of -scaling can be observed between @xmath197 and @xmath202  gev , due to the evolution of the coupling constant and of parton densities .
we also explored the independent renormalization- and factorization - scale dependence of our calculations .
the large size of this dependence at nlo was stressed already in refs .
the results , for @xmath17 collisions at @xmath203 and 10  gev , are shown in figs .
[ fig : scdep5 ] and [ fig : scdep10 ] , respectively . with the exception of the renormalization - scale dependence at 5  gev ,
a significant improvement in the stability of the results , relative to the dependence at lo and nlo , is observed in all cases .
[ cols="^,^,^,^,^,^,^,^",options="header " , ]     the convergence of the higher - order corrections is displayed in table  [ tab : toppert ] .
the last column includes the sum of all contributions of order @xmath204 and higher , performed using the _ minimal prescription _ of eq .
( [ mppp ] ) .
the fixed - order terms do not have any ambiguity due to the choice of the contour for the mellin transformation in eq .
( [ mppp ] ) .
the contribution from the @xmath181 channel converges very rapidly . in the case of the @xmath101 channel
the convergence is slower , in particular at the larger values of @xmath8 , but even at @xmath205  gev the size of the resummed contributions beyond order @xmath204 is only of the order of 10% of the total .
this supports the validity of the _ minimal prescription _
, since the truncated resummed expansion converges to it very smoothly .
we present in this section a comparison between our calculations and the results of the two most recent measurements of fixed - target direct - photon production .
e706  @xcite studied photons produced in @xmath206 collisions at @xmath207  gev , covering the centre - of - mass rapidity range @xmath208 to indicate the value of photon rapidity in the hadron - hadron centre - of - mass frame . ] .
the data span the  region between 3.5 and 11  gev , approximately .
this corresponds to @xmath209 .
ua6  @xcite studied photons produced in @xmath210 collisions at @xmath211  gev , over the rapidity range @xmath212 .
the  region extends between 4 and 7  gev ( @xmath213 ) , approximately .
what follows is not meant to be a systematic phenomenological study , but will serve as a benchmark to assess the impact of the resummation effects in realistic experimental conditions .
all the calculations in this section have been done assuming @xmath214 . for recent complete studies of all available data , done using the fixed - order nl calculations ,
see refs .
@xcite .    to compare our resummed predictions with actual data , two additional things need to be done : inclusion of the @xmath125-suppressed contributions , and inclusion of realistic experimental cuts .
the class of processes for which we evaluated resummation corrections in section  [ secnotat ] provides the dominant contribution to the production rate in realistic experimental configurations . on the left - hand side of fig .
[ fig : e706frag ] we plot the relative contribution , evaluated at nlo , of the processes with a @xmath101 , @xmath181 , @xmath215 and @xmath216 initial state for the e706 experimental configuration .
as one can see , the sum of @xmath101 and @xmath181 accounts for 90% of the overall rate , independently of @xmath217 .
on the right - hand side of fig .
[ fig : e706frag ] we plot the rate of the direct contributions , relative to the sum of direct and fragmentation , for each given channel . in the case of the @xmath215 and @xmath216 channels we compare the absolute values of the rates , since the direct component is negative after the subtraction of the initial - state mass singularities . the comparison of the two plots in fig .
[ fig : e706frag ] shows that the processes for which we are going to include resummation corrections account at nlo for a fraction of the total rate between 70 and 90% , in the  range 412  gev .
the situation is even better in @xmath210 collisions ( see fig .  [
fig : ua6frag ] , obtained for the ua6 experimental configuration ) , where the to - be - resummed processes account for over 90% of the nlo rate .
we therefore expect that the neglect of the resummation corrections to the @xmath218 and to the fragmentation processes is only a minor correction to the overall picture . for the present study , and in addition to the resummed predictions for the processes listed in eq .
( [ eq : procsin ] ) , we will therefore only include the fixed - order nl determination of these remaining components of the direct - photon production process .
the fragmentation processes are evaluated using the nlo single - inclusive parton  distributions from  @xcite , convoluted with the grv photon fragmentation functions  @xcite .
we found very small sensitivity to the choice of the photon fragmentation functions .
as anticipated above , the comparison of our results with actual data requires the inclusion of realistic detector acceptance cuts .
the resummation formalism discussed so far allows the evaluation of the transverse - energy distributions integrated over the full range of rapidity for the observed photon .
this approximation is technically correct , provided the measurement is performed within a finite range in rapidity including the value of @xmath219 in the collision centre - of - mass frame .
this is because in the large-  limit all production is concentrated at @xmath219 . to include the effect of experimental rapidity cuts , which usually do include the @xmath219 point
, we therefore apply the following acceptance correction to our resummed cross sections : ^(res)(yy ) ^(res)(all   y ) .
the experimental configuration for most experiments of practical interest is wide enough that , for the relevant values of @xmath217 , the rapidity acceptance is very large .
we show two examples in fig .
[ fig : yaccept ] .    one set of curves gives the ratio @xmath220 , evaluated at nlo for the e706 experimental configuration .
the three curves correspond to different choices of scale @xmath38 , and show very small dependence on @xmath38 .
we also checked that the dependence on the pdf set used is at the level of 1 - 2% .
the acceptance loss is of the order of 25% for  values around 4  gev , and becomes totally negligible for @xmath221  gev .
the other set gives the ratio @xmath222 , evaluated at nlo for the ua6 experimental configuration .
the acceptance loss is here more significant , due to the tight cut at negative rapidity .
we present our prediction for the e706 data in fig .
[ fig : e706comp ] has been rescaled by a 1.09 factor , to account for nuclear corrections to the @xmath206 process  @xcite . ] .
notice the significant reduction in scale dependence obtained when going from the fixed - order nl calculation to the resummed result .
this scale reduction is particularly evident at high , where the resummation effects are more important . notice that while at low  the band with the resummed prediction is all contained within the nlo uncertainty band , at high  the nll result becomes larger relative to nlo for all the displayed scale choices .
the plot shows a reasonable agreement between data and theory at large , indicating that no additional significant contribution is required in this region .
the large disagreement between data and theory already present at nlo  @xcite is still present , as no net increase is obtained from the resummation contributions .
their only effect is to reduce the scale dependence .
a similar picture emerges from the comparison of the theory with the @xmath210 ua6 data .
this is shown in fig .
[ fig : ua6comp ] .
the disagreement between data and theory at low  is however much less dramatic here than in the case of e706 .
the extent to which these low-  discrepancies can be removed by the inclusion of non - perturbative effects such as an intrinsic @xmath12 remains to be understood , as the global consistency of the different data sets is not very compelling  @xcite .
we presented in this paper a numerical study of the impact of resummation corrections , at the next - to - leading logarithmic level , on the transverse energy distribution of direct photons produced in hadronic collisions .
we dealt with the resummation of the @xmath223 sudakov logarithms studied theoretically in refs .
@xcite . as a result
, this work is mostly of relevance for typical fixed - target photon production .
the current prompt - photon data from the high - energy hadronic colliders cover in fact the region @xmath224 , where sudakov effects are negligible .
we showed that the inclusion of higher - order sudakov corrections improves significantly the factorization and renormalization scale dependence , relative to what observed in the fixed - order nl calculations .
even when the scales are varied independently , the uncertainty from scale variations is significantly reduced . as a result of the reduced scale dependence ,
the overall size of the resummation contributions depends significantly on the chosen scale . in general , however , the resummed cross sections for different scale choices have values contained within the nlo uncertainty band for   values up to 0.5 , and exceed the upper side of the nlo band by large factors when  approaches 1 .
still , our resummation corrections turn out to be much smaller , at least in the range of existing data , than the effects induced at large  by some implementations of intrinsic @xmath12 effects , as discussed in refs .  @xcite . in these papers
the effect of intrinsic-@xmath12 corrections was evaluated to be as large as factors of 3 and more , for the whole  range .
we believe that this is related to the absence in the intrinsic-@xmath12 models of the appropriate sudakov suppression due to the presence of the hadronic system recoiling against the photon ( represented at ll in our formalism by the negative terms in the exponents of eqs .
( [ dlresqg  ] ) and  ( [ dlresqq  ] ) ) .
in our work we did not include resummation corrections to the fragmentation processes .
we proved that , in @xmath17 collisions , the large-@xmath60 behaviour of the corrections to the @xmath225 processes is formally similar to that of the corrections to the leading one , @xmath226 .
these corrections are therefore not suppressed when @xmath60 increases towards larger values . whether they can be neglected or not , is therefore a pure matter of numerics .
we showed that their contribution is not dominant in the  regions of experimental interest , and limited ourselves to including them at the fixed next - to - leading order .
as discussed in sec .
[ fragcomp ] , we have no reasons to believe that the resummation corrections are any larger for the fragmentation processes than for the @xmath101 channel . a more quantitative study of these statements , and a complete phenomenological assessment of the comparison between theory and the current sets of data in view of the results presented in this paper , will be the subject of future work .
in this appendix we recall ( see ref .
@xcite ) the explicit expressions of the various factors that contribute to the resummed cross sections in eqs .
( [ gammaresqq ] ) and ( [ gammaresqg ] ) .
we use the customary notation for the colour factors in @xmath227 qcd , namely , @xmath228 and @xmath229 .
the @xmath60-moments of the lo partonic cross sections in eqs .
( [ siqqgamma ] , [ siqggamma ] ) are [ qgamman ] _
q|qg , n^(0 ) = e_q^2 ( 2+n ) , [ ggamman ] _ qgq , n^(0 ) = _ |qg , n^(0 ) = e_q^2 ( 7 + 5n ) , where @xmath230 is the euler @xmath231-function .    the formulae for the nll functions @xmath89 in the exponent of the radiative factors in eq .
( [ deltanll ] ) are the following [ g2qq ] g_q|q^(2 ) ( , ; ) & = & ( 2c_f - c_a ) h^(2 ) ( ) + 2 c_a h^(2)(/2 ) + & + & 2 ( 1 - 2 ) + ( 1- ) - + & + & \ { + } , [ g2qg ] g_qg^(2 ) ( , ; ) & = & c_a h^(2 ) ( ) + 2 c_f h^(2)(/2 ) + & + & 2 ( 1 - 2 ) + ( 1- ) - + & + & \ { + } , where @xmath232 is the euler number and @xmath233 are the first two coefficients of the qcd @xmath87-function [ betass ] b_0 = , b_1 = .
the auxiliary function @xmath234 that appears in eqs .
( [ g2qq ] , [ g2qg ] ) has the following expression [ htnll ] h^(2 ) ( ) = - ( 1 - 2 ) - , where the coefficient @xmath235 is given by [ kcoef ] k = c_a ( - ) - t_r n_f .
the first - order coefficients @xmath104 and @xmath105 of the @xmath60-independent functions in eq .
( [ cgamma ] ) are c_q|q^(1)(q^2/^2;q^2/_f^2 ) & = & ^2 ( 2 c_f - c_a ) + - ( 2 c_f - c_a ) 2 + [ cqq1coef ] & + & k - k_q + ( 2 c_f - c_a ) + ( 2 c_f - c_a ) ^2 2 + & - & ( 2 c_f - c_f ) - b_0 , c_qg^(1)(q^2/^2;q^2/_f^2 ) & = & ^2 ( c_f + c_a ) + - ( c_f - 2 c_a ) 2 + [ cqg1coef ] & - & k_q + ( 2 c_f + 19 c_a ) + c_f ^2 2 + & - & ( ( c_f + c_a ) - c_f - b_0 ) - b_0 , where [ kcoeffq ] k_q = ( - ) c_f , and the coefficient @xmath235 is given in eq .
( [ kcoef ] ) .
note that the ll functions @xmath88 are given in eq .
( [ g1fun ] ) .
thus , the formulae presented in this appendix complete all the ingredients that are necessary to evaluate the resummed cross sections with nll accuracy .
99 # 1#2#3sov .
# 1 * ( 19#2 ) # 3 # 1#2#3nucl . phys . * b#1 * ( 19#2 ) # 3 # 1#2#3phys . rep . *
# 1 * ( 19#2 ) # 3 # 1#2#3phys . lett .
* b#1 * ( 19#2 ) # 3 # 1#2#3phys .
# 1 * ( 19#2 ) # 3 # 1#2#3phys . rev .
* d#1 * ( 19#2 ) # 3 # 1#2#3phys . rev .
* d#1 * ( 19#2 ) # 3 # 1#2#3z . phys .
* c#1 * ( 19#2 ) # 3 # 1#2#3eur .
j. * c#1 * ( 19#2 ) # 3 # 1#2#3comm .
# 1 * ( 19#2 ) # 3 # 1#2#3#4#1 * # 2 * ( 19#3 ) # 4 # 1#2#3_ibid .
_ * # 1 * ( 19#2 ) # 3 # 1 hep - ph/#1 # 1 hep - ex/#1 d. adams et al . , e704 coll .
, ; + c. de marzo et al .
, na24 coll .
, ; + a. angelis et al . , r108 coll . , ; + a. angelis et al . , r110 coll . , ; + e. anassontzis et al . , r806 coll . , ; + t. akesson et al . , r807 coll . , ; + g. ballocchi et al . , ua6 coll .
, ; + m. bonesini et al . , wa70 coll .
, ; + g. alverson et al .
e706 coll .
, ; + j. badier et al .
, ; + m. mclaughlin et al . , e629 coll .
. l. apanasevich et al . , e706 collaboration , .
g. ballocchi et al .
, ua6 collaboration , .
p.  aurenche et al . , .
w. vogelsang and a. vogt , .
martin , r.g .
roberts , w.j . stirling and r.s .
thorne , , and references therein .
lai et al . , cteq collaboration , preprint msu
hep903100 , , and references therein .
j. huston et al . , .
f. abe et al . ,
cdf collaboration , .
b. abbott et al . ,
d@xmath7 collaboration , , submitted to phys .
j.  huston et al . , , ( ) .
l. apanasevich et al . , .
p. aurenche et al . , preprint lapth70498 , .
collins and d.e .
soper , ; + j.c .
collins , d.e .
soper and g. sterman , .
s.  catani , m.  ciafaloni and f.  hautmann , , ; + j.c .
collins and r.k .
ellis and d.a .
e. laenen , g. oderda and g. sterman , , ( ) .
s.  catani , m.l .  mangano and p.  nason ,
jhep 9807 ( 1998 ) 024 , ( ) .
s. catani , m.l .
mangano , p. nason and l. trentadue , , , ( ) .
p.  aurenche , r.  baier , a.  douiri , m.  fontannaz and d.  schiff , ; + p.  aurenche , r.  baier , m.  fontannaz and d.  schiff , .
f.  aversa , p.  chiappetta , m.  greco and j.ph .
guillet , .
g.  sterman , . s.  catani and l.  trentadue , .
s.  catani and l.  trentadue , .
r.  bonciani , s.  catani , m.l .  mangano and p.  nason , , ( ) .
contogouris , n.  mebarki and s.  papadopoulos , int . j. mod
a5 * ( 1990 ) 1951 ; + a.p .
contogouris and s.  papadopoulos , mod .
* a5 * ( 1990 ) 901 .
gordon and w.  vogelsang , .
s.  catani , g.  turnock , b.r .
webber and l.  trentadue , ; + s.  catani , l.  trentadue , g.  turnock and b.r .
lai et al . , .
m. glck , e. reya and a. vogt , , erratum .