arxiv:hep-ph/ Nov 1995

Transcription

1 DESY 95-6 ISSN he-h/95405 November 995 Incluive Dijet Production at HERA: Direct Photon Cro Section in Next-To-Leading Order QCD arxiv:he-h/ Nov 995 M. Klaen, G. Kramer II. Intitut für Theoretiche Phyik Λ Univerität Hamburg D - 76 Hamburg, Germany Abtract We have calculated incluive two-jet cro ection in next-to-leading order QCD for direct hotoroduction in low Q e colliion at HERA. Infrared and collinear ingularitie in real and virtual contribution are cancelled with the hae ace licing method. Analytical formula for the different contribution giving the deendence on the licing arameter are reented. Variou one- and two-jet ditribution have been comuted demontrating the flexibility of the method. Λ Suorted by Bundeminiterium für Forchung und Technologie, Bonn, Germany under Contract 05 6HH93P(5) and EEC Program "Human Caital and Mobility" through Network "Phyic at High Energy Collider" under Contract CHRX-CT (DG COMA)

2 Introduction Electron-roton cattering at HERA roceed dominantly through the exchange of hoton with very mall virtuality (Q ' 0). An imortant fraction of the fl colliion contain high tranvere momentum (E T ) rocee. The reence of thi hard momentum cale allow the alication of erturbative QCD to redict cro ection for the roduction of two or more high-e T jet. In lowet order (LO) QCD, the hard cattering rocee can be claified in two tye [, ]. In the o-called direct roce, the hoton interact in a oint-like fahion with the quark which aear either in the QCD Comton cattering flq! gq or the hoton-gluon fuion roce flg! qμq. In the reolved roce, the hoton act a a ource of arton, which can catter with the arton coming out of the roton. Thee two rocee give two ditinctly different tye of high-e T event toologie. The oint-like interaction reult in event with a characteritic (+)-jet tructure, i.e. two high-e T jet and one jet of low E T fragment of the roton. The reolved hoton event have acharacteritic (+)-jet tructure, ince in addition to the two high-e T jet and the low-e T roton remnant, fragment coming from the hoton will roduce a econd low-e T jet. While in the (+)-jet event the total hoton energy contribute to the hard cattering roce, only a fraction articiate in the (+)-jet event. Several method for detecting the hoton remnant jet have been rooed and teted in Monte Carlo imulation, generally uing a cut on the energy in the hoton direction or other kinematical variable [3, 4]. Uing an E T cluter algorithm, the ZEUS collaboration wa actually able to elect event with two high-e T jet and a third cluter in the aroximate direction of the electron (hoton) beam [5]. They tudied the roertie of the hoton remnant jet and found them to agree with exectation. Thi how that hotoroduction event with a hoton remnant jet can be identified exerimentally. Unfortunately, thi doe not lead to a earation of the direct and reolved rocee, becaue their characterization decribed above i only valid in LO QCD. In NLO, the direct cro ection may alo have contribution with a hoton remnant jet. Therefore it detection doe not uniquely earate the reolved cro ection. In a comlete NLO calculation the direct contribution with a hoton remnant jet could be iolated, o that the remaining contribution could be comared with the (+)-jet configuration. Of coure, in thi cae the election rocedure in the exerimental analyi and in the theoretical calculation mut match. To make further rogre in thi direction, a comlete NLO calculation of the direct cro ection i needed which i flexible enough to incororate the oible election criteria for the high-e T jet and the low-e T remnant jet. NLO calculation of direct hotoroduction have been done in the at for incluive ingle hadron [6], incluive ingle jet [7, 8, 9], and for incluive dijet [7, 0] cro ection. Mot of thee reult are alicable only for ecial obervable in incluive ingle jet and dijet roduction. The leading order calculation of the direct comonent involve only two ubrocee: flq! gq and flg! qμq. By including term u to O(ffff ), new contribution arie. Thee are the one-loo contribution to flq! gq and flg! qμq and the three-body ubrocee flq! qgg, flq! qqμq, flq! qq 0 μq 0, and flg! gqμq. Beide the ultraviolet ingularitie in the one-loo contribution, which can be removed by renormalization, additional infrared and

3 collinear ingularitie occur. They are extracted with the dimenional regularization method. It i well known that for uitably defined incluive obervable the infrared ingularitie in the one-loo contribution will cancel againt thoe in the three-body tree diagram. The collinear ingularitie from initial tate radiation and q μq roduction at the hoton leg can be factorized and aborbed in the reective arton ditribution function. The hard collinear ingularitie from final tate radiation will cancel againt the correonding ingularitie in the one-loo grah. To achieve thi technically and to maintain alo the neceary flexibility for new choice of obervable, jet definition, and exerimental cut, we ue the o-called hae ace licing method. We introduce an invariant ma cut-off y, which allow u to earate the hae ace region containing the ingularitie. If y i choen ufficiently mall, the integration over the oft and collinear region can be erformed analytically in a traightforward way. Then y i jut a technical cut which need not be connected with the exerimental jet definition. Of coure, the uitably defined incluive one- and two-jet cro ection mut be indeendent ofy. Our method i very imilar to that of Baer et al. [7] who ued two ditinct cut-off arameter for oft and collinear ingularitie. The hae ace licing method with invariant ma cut-off ha been ued for many rocee. The calculation for hotoroduction roceed along the ame line a the calculation of higher order QCD correction for multijet cro ection in dee inelatic e cattering firt erformed by Graudenz []. The lan of the aer i a follow: In ection we review the LO calculation and fix our notation. The technique to erform the NLO calculation are exlained in ection 3. In ection 4 we reent ome reult for incluive one- and two-jet cro ection in order to demontrate the flexibility of our method. Section 5 contain a ummary of our concluion. In everal aendice we rovide our analytically calculated reult needed for the Monte Carlo comutation. Leading Order Cro Section. Photon Sectrum The cro ection we have comuted are eentially for kinematical condition a in the HERA exeriment. There, electron of E e =6:7 GeV roduce hoton with mall virtuality Q = q ' 0, which then collide with a roton beam of energy E = 80 GeV. q = e 0 e i the momentum tranfer of the electron to the hoton. The ectrum of the virtual hoton i aroximated by the Weizäcker-William formula x a F fl=e (x a )= ff ß (+( x a) )ln ψ! Q max ( x a ) m e x a ; () where m e i the electron ma and x a = q e ' E fl =E e i the fraction of the initial electron energy tranferred to the hoton. i the four momentum of the incoming roton. In the equivalent hoton aroximation, the cro ection for the roce e! e 0 X with arbitrary final tate X i then given by the convolution dff(e! e 0 X)= Z x a;min dx a F fl=e (x a )dff(fl! X); () where dff(fl! X) denote the cro ection for fl! X with real hoton of energy E fl = x a E e.

4 . Jet Cro Section in LO In leading order, the differential cro ection for two-jet roduction take a very imle form. The artonic ubrocee are flq! gq and flg! qμq, which we denote generically by flb!. b i the arton emitted from the roton with momentum 0 = x b. It ditribution function in the roton F b= (x b ;Mb ) deend on the momentum fraction x b and the factorization cale M b. The final tate arton have momenta and, which can be exreed by their tranvere momentum E T and their raiditie and. The convention i that the z direction i arallel to the roton beam a ued in the exerimental analyi. From energy and momentum conervation one obtain x a = E T e + e ; (3) E e x b = E T E (e + e ) : (4) Thu, the kinematical variable of the two jet are related to the caling variable x a and x b. For examle, the raidity of the econd jet i kinematically fixed by E T,, and the hoton fraction x a through the relation xa E = e ln e : (5) E T In the HERA exeriment, x a i retricted to a fixed interval x a;min x a x a;max <. We hall diregard thi contraint and allow x a to vary in the kinematically allowed range x a;min x a, where x a;min = E E T e E e E E e E T e : (6) From equation (3) and (4), we can exre x b a a function of E T,, and x a : x b = x a E e E T e x a E e E E E T e : (7) Thi function i lotted in Fig. for E T = 0 GeV and four different value of in the range 0 3 a a function of x a. Alo hown i the enveloe of thee curve x b;env = E T x a E e E ; (8) which can be obtained from finding the minimum of x b with reect to for fixed x a. The minimal x a -value for = 0,,, 3 are 0.379, 0.43, , and The raidity of the econd jet i given by equation (5). Deending on x a, different region of x b are exlored. Small value of x b are reached for negative and the maximal oible x a, which ix a = in our cae and x a = x a;max in the HERA exeriment. Thi i of interet for determining the gluon tructure function of the roton in the mall x region []. Since under exerimental condition x a i not contrained very much, it i oible to achieve minimal value of x b by electing the aroriate range in E T,, and according to (3) and (4). 3

5 Figure : Contour of contant in the x a x b lane for E T = 0 GeV. x b;env (full curve) i the enveloe of the curve. Of coure, tudying the gluon tructure of the roton thi way deend very much on the cro ection in the aroriate region. For direct hotoroduction, the two-jet cro ection for e! e 0 + jet + jet + X i obtained from d 3 ff de T d d = X b x a F fl=e (x a )x b F b= (x b ;M b ) dff dt (flb! ): (9) dff/dt tand for the differential cro ection of the arton ubroce flb!. The invariant of thi roce are =( 0 + q) ;t=( 0 ), and u =( 0 ). They can be exreed by the final tate variable E T,, and and the initial tate momentum fraction x a and x b : = 4x a x b E e E ; (0) t = x a E e E T e = x b E E T e ; () u = x a E e E T e = x b E E T e : () So, the deendence of the two-jet cro ection on E T,, and i determined through the hoton ditribution function F fl=e, the arton ditribution function of the roton F b=, and the cro ection of the hard ubrocee, which deend on, t, and u. For the incluive one-jet cro ection, we mut integrate over one of the raiditie in (9). We integrate over and tranform to the variable x a uing (3). The reult i the cro ection for e! e 0 + jet + X, which deend on E T and : d ff de T d = X b Z x a;min Here, x b i given by (7) with =. dx a x a F fl=e (x a )x b F b= (x b ;M b ) 4E e E T dff x a E e E T e dt (flb! ): (3) 4

6 The cro ection for the two ubrocee flq! gq and flg! qμq are well known [3] and are given by dff dt (flq i! gq i ) = ß ffff Q i C F dff dt (flg! q iμq i ) = ffff Q i u u ß u t + t u ; (4) ; (5) where C F =(N C )=(N C ) and N C i the number of color. Thee cro ection are O(ffff ). The index i denote the quark flavor and Q i the quark charge. The magnitude of the hotongluon cro ection determine whether one i enitive to the gluon tructure function of the roton. 3 Next-To-Leading Order Cro Section The higher order correction are calculated a uual with the hel of dimenional regularization. Thi method regularize the ultraviolet ingularitie in the one-loo contribution, which are then ubtracted uing the modified minimal ubtraction (MS) cheme. Dimenional regularization can alo be ued for the infrared and collinear ingularitie. In the following ubection, we hall conider the virtual correction and the variou real correction needed to obtain finite cro ection in the limit dimenion n! Virtual Correction u to O(ffff ) The one-loo diagram for flq! gq and flg! qμq have an additional virtual gluon, which lead to an extra factor ff.interfering with the LO diagram of order O(ffff ) roduce the virtual correction to the! cro ection u to O(ffff ). Thee contribution are well known for many year from e + e! qμqg higher order QCD calculation [4, 5]. For the correonding hotoroduction cro ection, one ubtitute Q = 0 and erform the neceary croing. The reult can be found in [6, 6, 7]. We have alo comared with the reult in [] for dee inelatic cattering eq! e 0 gq and eg! e 0 qμq, which can be exreed by the invariant, t, and u after etting Q = 0. The tructure of the correction deend on the color factor C F, C F N C, and C F N f, where N f i the number of flavor in the qμq loo. For comletene and for later ue, we write the final reult in the form H V (flq i! gq i ) = e g μ 4ffl 8( ffl)q i ff ß ψ 4ßμ! ffl ( ffl) ρ C F V q (; t; u) N CC F V q (; t; u) ( ffl) 4 T q (; t; u) ψ! +C F ffl +ln μ 3 N f 6 N C ψ! H V (flg! q i μq i ) = e g μ 4ffl ff 8( ffl)q 4ßμ ffl ( ffl) i ß ( ffl) ρ C F V g(; t; u) 4 N CV g(; t; u) 4 ) +O(ffl); (6) (7) + ψ! ffl +ln μ 3 N f ) 6 N C T g (; t; u) +O(ffl): 5

7 H V give the virtual correction to the correonding reaction u to the n-dimenional hae ace factor dps () 4ß ffl = (8) dt ( ffl) tu 8ß and the flux factor =(). The exreion for V q, V q, V g, and V g are collected in aendix A. They contain ingular term / =ffl and =ffl (ffl =4 n), which are alway roortional to the LO cro ection T q = ( ffl) T g = ( ffl) u +ffl; (9) u t u + u ffl (0) t (ee eq. (4) and (5)). g i the quark-gluon couling and μ i the renormalization cale. The term roortional to ( 3 N f 6 N C) reult from the renormalization counter term. The initial gluon in are averaged by alying a factor of =(( ffl)). 3. Real Correction u to O(ffff ) For the calculation of the hard cattering cro ection u to O(ffff ), we mut include all diagram with an additional arton in the final tate, i.e. flq! qgg, flq! qqμq, flq! qq 0 μq 0, and flg! gqμq. The four-vector of thee ubrocee will be labeled by q 0! 3, where q i the momentum of the incoming hoton and 0 i the momentum of the incoming arton. The invariant will be denoted by ij =( i + j ) and the reviouly defined Mandeltam variable, t, and u. For male arton, the! 3 contribution contain ingularitie at ij =0. They can be extracted with the dimenional regularization method and cancelled againt thoe aociated with the one-loo contribution. To achieve thi we go through the following te. Firt, the matrix element for the! 3 ubrocee are calculated in n dimenion. The correonding diagram are hown in Fig. and are claified a I, II,..., VII. They are quared and averaged/ummed over initial/final tate in and color and can be ued for ingoing quark, antiquark, and gluon with the hel of croing. The diagram III' account for ghot to cancel unhyical olarization of the gluon. Thi claification i the ame a in ref. []. Since the hadronic tenor i ymmetric,we need to conider only the roduct I I, II I, III I etc. Furthermore, the roduct V IV, VI IV, VII IV, VI V, VII V, and VII VI are not ingular and vanih with the cut-off decribing the boundary of the ingular region. Such term will be neglected conitently in the following. Next, one can ditinguih three clae of ingularitie in hotoroduction (Q ' 0) deending on which roagator in the quared matrix element vanih. Examle for thee three clae are hown in Fig. 3. The X mark the roagator leading to the divergence. In the firt grah, it i the invariant 3 given by momenta of the final tate. Therefore, thi ingularity will be called final tate ingularity. The econd grah become ingular for q =(q + ) =0, when the hoton and the final quark momentum are arallel. Thi i the cla of hoton initial tate ingularitie. In the third grah, the ingularity occur at 03 =( ), where 0 i the initial arton momentum. Thi tand for arton initial tate ingularity. The firt cla i familiar from imilar calculation for jet roduction in e + e annihilation [6], the third cla from jet roduction in dee inelatic e cattering (Q 6= 0) []. The econd cla occur only 6

8 I II III III IV V VI VII Figure : Three-body diagram for clae I,..., VII a decribed in the text. q q q X X X 3 a) b) c) Figure 3: Three-body diagram with final and initial tate ingularitie. for hotoroduction [9, 0]. When quaring the um of all diagram in Fig., we encounter term, where more than one of the invariant become ingular, e.g. when one of the gluon momenta 3! 0 o that 3 = 0 and 03 = 0. Thee infrared ingularitie are dientangled by a artial fractioning decomoition, o that every term ha only one vanihing denominator. Thi alo allow the earation of the different clae of ingularitie in Fig. 3. Excet for additional term earating the hoton initial tate ingularity, we follow the rocedure in []. It turn out that the reult are alway roortional to the LO cro ection involved in the hard catter like T q and T g in eq. (9) and (0) H F;I;J = KT i : () Here, F; I, and J denote final tate, hoton initial tate and roton initial tate contribution. A the lat te, the decomoed matrix element have tobeintegrated u to ij y. y 7

9 characterize the region, where the two arton i and j cannot be reolved. Then, the ingular kernel K roduce term / =ffl and =ffl, which will cancel againt thoe in the virtual diagram or be aborbed into tructure function, and finite correction roortional to ln y,lny, and y 0.Term of O(y) will be neglected. In the following, we hall give the reult for the different clae of ingularitie earately. 3.. Final State Singularitie In thi ection, we aume that after artial fractioning the! 3 matrix element are ingular only for =0. For the integration over the ingular region of hae ace, we chooe a the coordinate ytem the c.m. ytem of arton and. The angle of the other arton three-momenta 0 and 3 with reect to and are hown in Fig. 4. χ i the angle between 0 and 3, i the angle between 0 and, and ffi i the azimuthal angle between the lane defined by 0 and and 0 and 3, reectively. 0 θ χ φ 3 Figure 4: Diagram for the three-body final tate defining the kinematical variable in the c.m. ytem of arton and. In general, the invariant = ( 0 + q) ; () t = ( 0 ) ; (3) u = ( 0 3 ) (4) differ from the correonding two-body invariant (ee eq. (0)-()), but aroach them for =0or arallel to. Then, t =( 0 μ ) and u =( 0 3 ), where μ = + i the four-momentum of the recombined jet and 3 i the four-momentum of the econd jet. Thi roduce a jet ma μ =( + ) 6= 0, which vanihe in the limit y! 0. For =0,we have μ =, and for arallel to wehave μ /. The definition of t and u are unique only in the limit! 0. Other choice for the three-arton final tate variable differ by finite contribution of O(y), which are neglected throughout. We ubtitute the angle by b = ( co ) (5) 8

10 and introduce the new variable z 0 = 0 q : (6) The three-body hae ace in n dimenion can be factorized into two two-body hae ace term [] dps (3) = dps () dps (r) ; (7) where and with dps () dt 4ß dps (r) = N b and N ffi are normalization factor: = ffl ( ffl) 4ß ffl tu 8ß (8) ( ffl) ( ffl) 6ß ffl dμ F (9) ψ! ffl dμ F =dz 0 z 0 ffl + z0 db b ffl ffl dffi ( b) in ffl ffi: (30) t N b N ffi N b = N ffi = Z 0 Z ß 0 db(b( b)) ffl = ( ffl) ( ffl) ; (3) in ffl ffidffi = ß4ffl ( ffl) : (3) ( ffl) dps () i the ame hae ace factor a in the two-body LO final tate. The full range of integration in dps (r) i given by z 0 [0; t=], b [0; ], and ffi [0;ß]. The ingular region i defined by the requirement that arton and are recombined, i.e. =0or arallel to, o that =0. Weintegrate over thi region u to y, which retrict the range of z 0 to 0 z 0 minf t=; yg y F. The matrix element, which are obtained by quaring the um of diagram in Fig., are exreed by the variable, t, u, z 0, b, and ffi. Deending on the tructure of the interfering diagram and the color factor, we ditinguih even tye of contribution: F, F,..., F 7. Five of them are for an incoming quark or antiquark and two for an incoming gluon. The definition of the F i with the correonding color factor can be read off Tab.. In cla F,wehave QED-tye diagram roortional to C F, where a gluon can become oft or collinear with a quark. F contain QED-tye diagram and interference term with diagram containing a three-gluon vertex roortional to C F N C. Cla F 3 tand for contribution of thi latter tye, where the infrared divergence i caued by the two-gluon final tate. F 4 ha only contribution from non-abelian diagram, and F 5 contain the four-quark diagram. In thi cla, the divergence originate from the collinear qμq air. In F 6 and F 7,wehave contribution with incoming gluon roortional to C F and N C, reectively. We emhaize that thi claification ret uon the ue of the Feynman gauge for the gluon in ummation with all ghot contribution included. Of coure, the um of all term i gauge invariant. After integrating the! 3 matrix element with final tate ingularitie H Fi,we obtain Z ψ! dps (r) H Fi = e g μ 4ffl ff 8( ffl)q 4ßμ ffl ( ffl) i ß ( ffl) 4 F i +O(ffl): (33) 9

11 Cla Inc. Parton Diagram Color Factor F Quark I I, II II, II I N C C F =N C F Quark II I, III I, III II = N CC F =N C F 3 Quark III I, III II = N CC F =N C F 4 Quark III III N CC F =N C F 5 Quark IV IV, V V, VI VI, VII VII = N C C F =N C F 6 Gluon I I, II II, II I N C C F =N C C F F 7 Gluon II I, III I, III II = N CC F =N C C F Table : Claification of the! 3 matrix element for final tate ingularitie The final reult for the F i are written down in aendix B. They contain the infrared and collinear ingularitie which cancel againt thoe in the virtual correction. At thi oint, it i eential that the ingular term are all roortional to the LO matrix element T q (; t; u) and T g (; t; u) and that the variable, t, and u defined in ()-(4) correond to the two-body variable in the infrared and collinear limit with an aroriate definition of the recombined jet momentum. 3.. Photon Initial State Singularitie To integrate over the ingularity q3 =0,we ue the ame, c.m. coordinate ytem a in the lat ection. Let u aume that we have iolated the term in the three-body cro ection with a ingularity at q3 = 0, i.e. where the outgoing quark momentum 3 i collinear with the hoton momentum q. In thi cae, the quark 3 i art of the hoton remnant and and are the arton momenta of the two-body ubroce. We introduce a new variable z a = 0 q [X a; ]; (34) where X a = 0 e ' E q =E e i now the fraction of the initial electron energy tranferred to the quark in the hard catter. Then, z a arametrize the hae ace of the collinear quark 3 =( z a )q and render the Mandeltam variable in the form With the additional variable the three-body hae ace factorize into = ( 0 + z a q) ; (35) t = ( 0 ) ; (36) u = ( 0 ) : (37) z 00 = q 3 q 0 ; (38) dps (3) = dps () dps (r) : (39) dps () i identical to (8), and 4ß dps (r) = ffl ( ffl) ( ffl) 6ß H a(z 00 )dμ I ; (40) 0

12 where and dμ I =dz 00 z 00 ffl dz a z a za z affl dffi H a (z 00 )= in ffl ( ffl) ffi N ffi ( ffl) (4) ψ+ z00 z a! +ffl ψ z00 z a! ffl =+O(z 00 ) (4) can be aroximated by a it lead only to negligible term of O(y). The full region of integration i given by z 00 [0; u=], z a [X a ; ], and ffi [0;ß], wherea the ingular region i defined by the requirement that arton 3 i collinear to q, o that q3 < y=z a. Thi retrict the range of z 00 to 0 z 00 minf u=; yg y I. The matrix element from Fig. can be ut into even clae I, I,..., I 7 dilayed in Tab.. The claification i again according to incoming arton, color factor, and general tructure. Cla Inc. Parton Diagram Color Factor I Quark I I, II II, II I N C CF =N C I Quark II I, III I, III II = NCC F =N C I 3 Quark IV IV, V V, VI VI, VII VII = N C C F =N C I 4 Quark IV IV, V V, VI VI, VII VII = N C C F =N C I 5 Quark IV IV, V V, VI VI, VII VII C F =N C I 6 Gluon I I, II II, II I I 7 Gluon II I, III I, III II N C CF =N C C F = NCC F =N C C F Table : Claification of the! 3 matrix element for hoton initial tate ingularitie In all clae, the incoming hoton can become collinear to a quark line only, ince there i no direct couling to a gluon line. Furthermore, there are only ingle ole, i.e. only collinear ingularitie. The three-body matrix element are factorized into LO two-body arton-arton cattering matrix element and a ingular kernel K. Since the quark (antiquark) roduced by the hoton become art of the hoton remnant, the remaining two-body roce i induced by an antiquark (quark) in the initial tate, which interact with a quark or antiquark of the ame or differentflavor or a gluon from the roton. The factorized LO arton-arton cattering cro ection are given by the function U i (; t; u) (i =; ; 3) related to the variou rocee a in Tab. 3. The function U i (; t; u) are decomoed with reect to their color tructure a follow: U = C F N C ut ( ffl) ( ffl) t u + u t ffl (43) = C F T a (; t; u) N C T b(; t; u); U = ψ! + u ffl = t T (; t; u); (44) U 3 = N C ( ffl) ψ tu + ffl! = N C T 3 (; t; u): (45)

13 Proce Matrix element jmj =8N C C F g 4 U i qq! gg U (; t; u) qq 0! qq 0 U (; t; u) qq 0! qq 0 U (u; t; ) qq! q 0 q 0 U (t; ; u) qq! qq U (; t; u)+u (; u; t)+u 3 (; t; u) qq! qq U (u; t; )+U (u; ; t)+u 3 (u; t; ) qg! qg U (t; ; u) Table 3: LO arton-arton cattering matrix element The function T a ;T b ;T, and T 3 aear in the final reult Z Z ψ! dz dps (r) a H Ii = e g μ 4ffl 8( ffl) ff 4ßμ ffl ( ffl) X a z a ß ( ffl) 4 I i +O(ffl); (46) which i obtained from the integration of the decomoed matrix element H Ii over dps (r).iti clear that the hoton initial tate ingularitie lead to contribution familiar from LO reolved hotoroduction. to The term I i are collected in aendix C and how exlicitly the ole in =ffl roortional h z a +( z a ) i : (47) P qi ψfl(z a )=N C Q i Thi function aear in the evolution equation of the hoton tructure function a an inhomogeneou or o-called oint-like term. Therefore, the hoton initial tate ingularitie can be aborbed into the hoton tructure function. The neceary te are well known [6, 8]. We define the renormalized ditribution function of a arton a in the electron F a=e (X a ;Ma)a F a=e (X a ;M a )= Z X a dz a z a ffi afl ffi( z a )+ ff ß R aψfl(z a ;M a ) F fl=e Xa z a ; (48) where R ha the general form ψ! R aψfl (z a ;Ma)= ffl P ( ffl) 4ßμ ffl aψfl(z a ) + C ( aψfl(z a ) (49) ffl) with C = 0 in the MS cheme. The renormalized artonic cro ection for flb! jet i then calculated from dff(flb! jet) = dμff(flb! jet) ff ß dμff(flb! jet) i the higher order cro ection before the ubtraction rocedure, and dff(ab! Z M a dz a R aψfl (z a ;Ma )dff(ab! jet): (50) jet) contain the LO arton-arton cattering matrix element U i (; t; u). The factor 4ßμ M a (49) i combined with the factor 4ßμ in (46) and lead to Ma deendent term of the form "ψ 4ßμ ffl P aψfl(z a )! ffl ψ! 4ßμ ffl M a = P aψfl (z a )lnψ M a! in ; (5) o that the cro ection after ubtraction in (50) will deend on the factorization cale M a.

14 3..3 Parton Initial State Singularitie The comutation of the arton initial tate contribution i very imilar to the lat ection. We again ue the, c.m. coordinate ytem a defined in Fig. 4. Singularitie occur in thi cae for 03 = 0, when the incoming arton 0 radiate a oft gluon ( 3 = 0) or a collinear arton 3. The new variable = 0 q [X b; ] (5) decribe the momentum of arton 3 =( ) 0 and the reduced momentum of the arton involved in the hard catter 0. X b = ' E Pq q=g =E i the energy fraction of the incoming roton involved in the! ubroce. The correonding invariant are Uing the variable we can again factorize the three-body hae ace into = ( 0 + q) ; (53) t = ( 0 ) ; (54) u = ( 0 ) : (55) z 000 = q ; (56) dps (3) = dps () dps (r) ; (57) where dps () tand for (8). 4ß dps (r) = ffl ( ffl) ( ffl) 6ß H b(z 000 )dμ J ; (58) where and dμ J =dz 000 z 000 ffld zb ffl dffi in ffl ( ffl) ffi N ffi ( ffl) (59) H b (z 000 )=( z 000 ) +ffl ψ z000! ffl =+O(z 000 ) (60) can again be aroximated by. The integration over z 000 [0; u=], [X b ; ], and ffi [0;ß] i retricted to the ingular region of z 000 in the range 0 z 000 minf u=; yg y J. The three-body matrix element can be comuted from Fig. and claified according to color factor, incoming arton, and tye of arton articiating in the two-body roce a in Tab. 4. In the cae of final tate ingularitie conidered in ection 3.., the ingoing arton (quark or gluon) wa alway the ingoing arton of the two-body ubroce. A could already been een in the lat ection, thi i not the cae for initial tate ingularitie. For examle, if an ingoing quark i collinear to the outgoing gluon (quark), the factorized -jet matrix element ha an initial quark (gluon). In Tab. 4, J i the cla of QED tye term, where a gluon in the final tate i collinear with the initial tate quark. The ame occur in J, where QED-tye and non-abelian diagram interfere. J 3 contain the four-quark diagram, in which an outgoing quark i collinear with 3

15 Cla Inc. Parton Inc. Parton Diagram Color Factor (3-jet) (-jet) J Quark Quark I I, II II, II I N C CF =N C J Quark Quark II I, III I, III II = NCC F =N C J 3 Quark Gluon IV IV, V V, VI VI, VII VII = N C C F =N C J 4 Gluon Quark I I, II II, II I N C CF =N CC F J 5 Gluon Quark II I, III I, III II = NCC F =N C C F J 6 Gluon Gluon II I, III I, III II = NCC F =N C C F J 7 Gluon Gluon III III NCC F =N C C F Table 4: Claification of the! 3 matrix element for arton initial tate ingularitie an ingoing quark. In J 4 and J 5,wehave an ingoing gluon and an outgoing collinear quark. The clae differ only by their color factor. In J 6 and J 7, an outgoing gluon i collinear with the ingoing gluon, where J 6 contain QED-tye and interference term and J 7 contain ure non-abelian term. The next te i to erform the z 000 integration with the uer limit y J. The reult can be written in the form Z Z dz dps (r) b H Ji = e g μ 4ffl ff 8( ffl)q i X b ß ψ 4ßμ! ffl ( ffl) ( ffl) 4 J i +O(ffl); (6) where the J i are given in aendix D. Beide, t, u, and y J, they till deend on the integration variable. The ingular art take the form J = J = J 3 = J 4 = J 6 = ffl P qψq ( )+ffi( ) C F ffl + ffl N CC F ffi( ) ffl ln u t P gψq ( ) ffl C F CF ln t + 3 C F T q (; t; u)+o(ffl 0 ); (6) T q (; t; u)+o(ffl 0 ); (63) T g(; t; u)+o(ffl 0 ); (64) ffl P qψg( ) C F T q (; t; u)+o(ffl 0 ); (65) P gψg ( )+ffi( ) ffl N C ffl + ln tu ffl N C 6 N C 3 N f (66) N C T g (; t; u)+o(ffl 0 ): 4 J 5 i O(ffl) and doe not contribute. J 6 contain contribution from J 6 a well a J 7. Some of the J i contain infrared ingularitie / =ffl, which mut cancel againt the correonding ingularitie in the virtual contribution. The ingular art are decomoed in uch a way that the Altarelli-Parii kernel in four dimenion roortional to =ffl are lit off. They alo aear in the evolution equation for the arton ditribution function and are defined exlicitly in aendix D. The ingular term roortional to thee kernel are aborbed a uual into the cale deendent tructure function F b= (X b ;M b )= Z X b d ffi bb 0ffi( )+ ff ß R0 bψb 0(;M b ) F b 0 = Xb ; (67) 4

16 where F Xb b0 = i the LO tructure function before abortion of the collinear ingularitie and R 0 bψb 0(;Mb )= ψ ffl P bψb 0(z ( ffl) 4ßμ b) ( ffl) M b! ffl + C 0 bψb 0() (68) with C 0 = 0 in the MS cheme. Then, the renormalized higher order hard cattering cro ection dff(flb!jet) i calculated from dff(flb! jet) =dμff(flb! jet) ff ß Z d R 0 bψb 0(;M b )dff(flb0! jet): (69) dμff(flb! jet) i the higher order cro ection before the ubtraction rocedure, and dff(flb 0! jet) contain the lowet order matrix element T q (; t; u) and T g (; t; u) inn dimenion. Thi well known factorization recrition [8] remove finally all remaining collinear ingularitie. It i univeral and lead for all rocee to the ame definition of tructure function if the choice concerning the regular function C 0 in (68) i ket fixed. Similar to the cae of hoton initial tate ingularitie, the higher order cro ection in (69) will deend on the factorization cale M b due to term of the form P bψb 0( )ln. 3.3 Jet Cro Section in NLO M b To obtain a finite cro ection, we mut add the four art conidered in ection 3. and Then, the ole in =ffl mut cancel, and we can take the limit ffl! 0. The reult i a ecial kind of two-jet cro ection, where the recombination of two arton into one jet or the recombination of a arton with the hoton or roton remnant jet i done with an invariant ma cut-off y. In LO, thi cro ection wa written down in (9). Including the NLO correction, it ha the form d 3 ff de T d d = X b " x a F fl=e (x a ) x b F q= (x b ;Mb ) dff dt (flq! )+ d~ff dt (flq! ) (70) +x b F g= (x b ;Mb )dff dt (flg! )+ d~ff dt (flg! ) +x b F b= (x b ;M b ) d~ff dt (ab! ) In (70), dff dt (flq! ) and dff dt (flg! ) tand for two-body contribution in LO and NLO together with analytically integrated contribution of the oft and collinear divergent region of the three-arton final tate. The contribution from the hoton initial tate ingularitie are denoted d~ff dt (ab! ). d~ff dt (flq! ) and d~ff dt (flg! ) originate from the integration of the oft and collinear ingularitie of the arton initial tate. The quark induced cro ection can be written a : where dff dt (flq! )=ff (μ )CT q + ff (μ ) ß C(T qa q + B q ); (7) A q = 6 (N C N f )ln μ + N f 3 5 ln y F 5 N C 3 6 ln y F (7)

17 t +(C F N C ) ln u t " 4 ln + N C ln y F y F ln u + u 4 ln + C F " ß u B q = 3C F ln u + + u 3 7 ln y F 3 ln y F (C F N C ) ln t ; +lny F ln t Li yf t yf Li ; ln u + t u + u ß ß +ln t u T q = u ßff ; and C = u Q i C F : (74) For the gluon induced cro ection, we have where (73) dff dt (flg! ) = 3 8 ff (μ )CT g + 3 ff (μ ) 8 ß C(T ga g + B g ); (75) A g = 6 (N C N f )ln μ + C F " +N C tu ln ß yf yf Li Li ; t u t u B g = 3C F ln u + u t ln t + + t t ln ; u T g = t u + u t : ln t +ln u ln y F 3lny F 6 y F ln t y F ln u +ln y F + t 4 ln + u 4 ln + (C F N C ) ln tu + + u u ln t The cut deendence of the final tate NLO correction in A q and A g i contained in the y F deendent term. For y F! 0, thee term behave like ( ln y F ), which lead to unhyical negative cro ection for very mall y F. Thu if y i ued a a hyical cut, it mut be ufficiently large. In mot of the alication, we hall ue thee reult for comuting incluive cro ection, in which the y deendence of the two-jet cro ection cancel againt the y deendence of the numerically calculated three-jet cro ection. The quadratic term in (ln y F ) are connected with the oft divergence of the three-arton final tate. The contribution from the hoton initial tate ingularitie necearily have a hoton remnant. They deend on the cut-off y I only through term roortional to (ln y I ). Thi contribution can be conidered earately only if y I iahyical cut, i.e. if it i ufficiently large. In cae one introduce other kinematical condition to define the hoton remnant term in the NLO direct cro ection, one mut add the aroriate three-body contribution. The two-jet cro ection i calculated from d~ff dt (ab! ) = d^ff dt (ab! ) ff " ß dz a P qi ψfl(z a ) X a z a Z 6 ψ ln ψ! ( za )y I z a M a! Q i +N C : (76) (77) (78) (79)

18 In (79), M a i the factorization cale at the hoton leg. Thi deendence mut cancel againt the M a deendence of the LO reolved contribution (ee eq. (48)). Thi ha been demontrated exlicitly in [9] for the incluive ingle jet cro ection. The d^ff dt (ab! ) are the cro ection for the! ubrocee which were given in ection 3.. in term of the function U i (; t; u). Thee are the well known LO cro ection alo ued for the calculation of hadron-hadron cro ection. From the arton initial tate ingularitie, we deduce d~ff dt (flq! ) = ff (μ ) ß CT q (80) ( Z d Xb " X b F q= ;Mb C F J ~ ( )+P qψq ( )ln X b Z X b d X b F g= Xb ;M b d X b F q= Xb ;M b " P qψg ( ) ψ ln M b ψ! ( zb )y J! ) + ; + X b Mb d~ff dt (flg! ) = 3 ff (μ ) 8 ß CT g (8) ( Z d Xb " X b F g= ;Mb N C J X b ~ 6 ( )+P gψg ( )ln Mb Z " ψ ψ!! ( zb )y J + P gψq ( ) ln + C F ff : M b The function ~ J i are the regular art of the function J i defined in aendix D. They contain the cut-off y J and are ditribution in. The dominant term for y J! 0 are again roortional to (ln y J ) and are related to the oft divergence in the arton initial tate. J doe not contribute after the aroriate croing of (t $ u). The reult in (70)-(8) are equivalent to thoe of Baer et al. [7]. They cannot be comared ince thee author introduce a different cut for integrating over the oft ingularitie. Intead of invariant mae, they ue a cut on the energy of the gluon. Term that are indeendent of oft ingularitie agree. For examle, thi i the cae for the cro ection in (79) coming from the collinear hoton initial tate ingularity. 4 Incluive One- and Two-Jet Cro Section In thi ection, we reent ome characteritic numerical reult for one- and two-jet cro ection which have been obtained with our method of licing the hae ace with invariant ma cut. In a hort communication we have ued thi method to calculate the differential two-jet cro ection a a function of average raidity μ = ( + ) of the two jet of highet tranvere energy with the cut on j j < 0:5 [0] and comared it to recent exerimental data of the ZEUS collaboration at HERA []. For thi aer, we have calculated variou one- and two-jet ditribution without alying ecial cut on kinematical variable of the initial or final tate dictated by the exerimental analyi, although our aroach i articularly uitable for thi 7

19 a ha been hown in [0]. The inut for our calculation i lited in the following. We have choen the CTEQ3M roton tructure function [], which i a NLO arametrization with MS factorization and Λ (4) = 39 MeV. Thi Λ value i alo ued to calculate the two-loo ff value at the cale μ = E T. The factorization cale i alo M = E T.We do not aly a ecial cut on the hoton energy fraction x a in () but integrate from x a;min to. Q max = 4 GeV i the ame a ued in the jet analyi of the ZEUS exeriment []. The further calculation i baed on two earate contribution a et of two-body contribution and a et of three-body contribution. Each et i comletely finite, a all ingularitie have been cancelled or aborbed into tructure function. Each art deend earately on the cut-off y. If y i choen large enough, the two art determine hyically well defined two-jet and three-jet cro ection. Our analytic calculation are valid only for very mall y, ince term O(y) have been neglected in the analytic integration. For very mall y, the two iece have no hyical meaning. In thi cae, the (ln y) term force the two-body contribution to become negative, wherea the three-body cro ection are large and oitive. In [0], we have lotted uch cro ection for y =0 3.However, when the two- and three-body contribution are uerimoed to yield a uitable incluive cro ection, a for examle the incluive inglejet cro ection, the deendence on the cut-off y will cancel. Then, the earation of the two contribution with the cut-off y i only a technical device. The cut-off only erve to ditinguih the hae ace region, where the integration are done analytically, from thoe where they are done numerically. Furthermore, y mut be choen ufficiently mall o that exerimental cut imoed on kinematical variable do not interfere with the cancellation of the y deendence. Firt, we conider the incluive ingle-jet cro ection. To achieve thi, we mut chooe a jet definition, which recombine two nearly collinear arton. A uual, we adot the jet definition of the Snowma meeting [3]. According to thi definition, two arton i and j are recombined if R i;j <R, where R i = q ( i J ) +(ffi i ffi J ) and J ;ffi J are the raidity and the azimuthal angle of the recombined jet. We chooe R = in all of the following reult. Thi mean that two arton are conidered a two earate jet or a a ingle jet deending whether they lie outide or inide the cone with radiu R around the jet momentum. In ome cae, it may occur that two arton i and j qualify both a two individual jet i and j and a a combined jet ij. In thi cae, we make no further election a it i done in [4], where only the combined jet i counted. In NLO, the final tate may conit of two or three jet. The three-jet amle conit of all three-body contribution, which do not fulfill the cone condition. In Fig. 5, the reult for the incluive one-jet cro ection d ff=de T d are hown for variou raiditie =0;, and a a function of E T. The NLO and the LO rediction are lotted. The LO curve i calculated with the ame roton tructure function and the ame ff a the NLO curve. Only the hard cattering cro ection i calculated in LO. For R =, thee two cro ection differ only very little. Of coure, the NLO cro ection deend on R, wherea the LO curve doe not. To earate the curve for the three value, the reult for =0( =) ha been multilied by 0: (0:5). The maximum of the one-jet cro ection i near =. For E T = 0 GeV, thi can be een from Fig. 6, where d ff=de T d i lotted a a function of between = and ' 4:. At thi 8

20 Figure 5: Incluive ingle-jet cro ection d ff=de T d a a function of E T for variou raiditie =0; ; in LO and NLO. The cro ection for =0( =)imultilied by 0: (0:5). Figure 6: Incluive ingle-jet cro ection d ff=de T d a a function of for E T = 0 GeV in LO and NLO. 9

21 value of E T, the kinematical boundarie, where the cro ection goe to zero, differ only little in LO and NLO. The k-factor for the NLO rediction i k ' : for the central region of 0 3. Our ingle-jet ditribution have been comared with the earlier reult of Bödeker [9], who ued the ubtraction method to cancel oft and collinear ingularitie. We found erfect agreement. To achieve thi, we had to ue a rather mall value of y =0 3. Furthermore, we had to add ome finite contribution of O() for y! 0, which correond to two-jet contribution where an outgoing arton i collinear to another outgoing or an ingoing arton without the correonding ole term in the! 3 matrix element after artial fractioning. We alo checked that our reult are indeendent of y, if y i ufficiently mall. The more eential reult are on incluive two-jet cro ection hown in the following figure. In Fig. 7, we reent d 3 ff=de T d d a a function of E T for = and variou choice of =0; ;. Here, E T and are the tranvere energy and the raidity of the o-called trigger jet. i the raidity of the econd jet, o that E T and E T are the jet with the highet tranvere momenta for the three-jet contribution. For exactly two jet in the final tate, we have E T = E T. In Fig. 7, we can ee how the cro ection decreae when i choen away from the maximum region at =. =ialway ket fixed. To dientangle the three curve, we have multilied the =( = 0) cro ection by a factor of 0:5 (0:) a in Fig. 5. The LO and NLO rediction are very much arallel. The k-factor i k ' :, imilar a in the one-jet cro ection. The only large difference between the LO and NLO cro ection occur at the boundary of the hae ace. For examle, the cro ection for = 0 ha an additional tail in the E T ditribution coming from the three-jet cro ection, where all three arton lie outide the cone with radiu R =. Such lot can be roduced alo for other choice of outide the region where the cro ection i maximal. We have tudied the incluive two-jet cro ection alo a a function of and for fixed E T. A an examle, we how in Fig. 8 the two-dimenional ditribution d 3 ff=de T d d for E T = 0 GeV in form of a lego-lot in the interval ; [ ; 4]. Over the whole range of and, the cro ection i aroximately ymmetric in and.we how thi cro ection in LO and NLO, and we ee clear difference between the LO and the NLO rediction. Thi become even clearer when we lot the rojection for fixed a a function of.in Fig. 9, thee rojection are hown for =0; ;, and 3 for LO and NLO. We ee that in NLO the cro ection i alway larger than the LO cro ection which i conitent with the k-factor larger than given above. It i clear that many more ditribution or artially integrated cro ection uing other two-jet variable can be calculated with the hae ace licing method. A a final alication, we have tudied the R and the y deendence of the two-jet cro ection, integrated over ; and E T with a minimal E T ;min = 0 GeV. The reult i hown in Fig. 0, where the LO cro ection i indeendent of R. The NLO cro ection i the um of oitive three-body contribution and negative two-body contribution, both calculated with the invariant ma cut y =0 3, and thi i not indeendent ofr. However, the deendence i rather weak. The cro ection increae with increaing R from R =0: tor =by a factor of and lie above the LO reult for R 0:6. For R<0:, the y cut become effective. Therefore we have lotted the curve only for R 0:. The um of the two-body and three-body contribution deend 0

22 Figure 7: Incluive dijet cro ection d 3 ff=de T d d a a function of E T for = and three value of =0; ;. The cro ection for =0( =)imultilied by0: (0:5). Figure 8: The LO and NLO trile differential incluive dijet cro ection for E T = 0 GeV a a function of and.

23 Figure 9: Projection of the LO and NLO trile differential incluive dijet cro ection for E T = 0 GeV a a function of for fixed =0; ;, and 3. on R like a + b ln R + cr due to the recombination of two arton momenta i and j inide the cone with radiu R. Thi way, the ingled out tranvere energy of a jet can differ from the individual arton E T '. Thi R deendence ha been tudied before for the incluive one-jet cro ection with imilar reult [8, 9]. We have alo calculated the ame cro ection a a function of the invariant ma cut y without the cone recombination. In thi cae, the cut y oerate a a hyical cut. Of coure, the reult for the two-body contribution (being negative for y<3 0 ) and the three-body contribution (being above the um for y<3 0 ) only make ene for ufficiently large value of the y cut. For y! 0, the two-body (three-body) contribution decreae (increae) like ln y (ln y). The reult are reented in Fig.. Here we ee that the um of both contribution deend only moderately on y (/ ln y) due to the recombination of arton in the final tate. Thi deendence on y of the ummed contribution i equivalent to the R deendence of the integrated cro ection hown in Fig. 0. We remark that in Fig. the reective cro ection mut till be corrected by finite term of O(y). The two-body contribution i calculated from the analytical formula given above, where thee term are neglected. Thi hould influence our reult only for y> 0. It i, however, no roblem to calculate thee term numerically.

24 Figure 0: Integrated incluive dijet cro ection for E T ;min = 0 GeV a a function of R in NLO (dahed-dotted) and LO (full). The body (dahed) and 3 body (dotted) contribution with y =0 3 are lotted earately. Figure : Integrated incluive dijet cro ection for E T ;min = 0 GeV a a function of y in NLO (dahed-dotted) and LO (full). The body (dahed) and 3 body (dotted) contribution are lotted earately. 3

25 5 Summary Differential cro ection d 3 ff=de T d d have been calculated in NLO for direct hotoroduction. Infrared and collinear ingularitie are cancelled with the hae ace licing method uing an invariant ma cut-off. Thi method allow to incororate variou cut on the final tate dictated by the analyi of exerimental data and to erform calculation for different choice of jet algorithm. Analytical formula for the different contribution giving the deendence on the licing arameter are derived. Thi i needed if one want to develo a comlete decrition of hotoroduction event including hadronization and other correction. Numerical reult for the two-jet incluive cro ection at HERA have been reented, and LO and NLO rediction have been comared. For a cone radiu of R =, the NLO correction lead to an increae of the order of 0% comared to the LO rediction which i calculated with the ame roton tructure function and the ame ff a the NLO rediction. 4

26 A Virtual Correction In thi aendix, we give the interference term of the virtual correction with the LO QCD Comton cattering flq! gq (= T q (; t; u)) and the LO hoton-gluon fuion flg! qμq (= T g (; t; u)) matrix element a defined in equation (6) and (7). The reult deend only on the Mandeltam variable, t, and u: V q = V q = V g = V g = " ffl + ffl ln u " ffl + ffl ln u " ln t 3 + ß 3 t 7+ln u +4ln t 3 u ln u + ß ln t ln u +4ln t + u + u T q (; t; u) ß +ln t u t 3 +ln T q (; t; u) u t ß +ln u ffl 3ffl ß + 3 t u 7+ln +ln +ln t +ln u +3u t " ffl + ffl ln t ln u +ln t +ln u + + u t + u T g (; t; u) t ln +3t u ln + u + u t + ß tu +ln T 3 g (; t; u) ln u + + t u ln t : + t ln u ; ln t ; (A.) (A.) (A.3) u ln + + t t ln u ; (A.4) B Final State Correction Thi aendix contain the real correction coming from the integration of the! 3 matrix element over the final tate ingular region of hae ace. In addition to the Mandeltam variable, the reult deend on the integration cut-off y F : yf F = F = F 3 = F 4 = " ffl + ffl ln t ln y F 3 ln y F t ln ffl " ffl + ffl yf + Li u + 3 ln t +lny F ln t Li C F T q(; t; u); t ln +lny F ln t Li yf t ln u 5 3ffl ln y F 3 9 y 4+lnyF +ln y F +ln F N CC F T q (; t; u); N CC F T q (; t; u); 5 t ß (B.) N CC F T q (; t; u); (B.) u u ln + ß 3 (B.3) (B.4)

27 F 5 = 3ffl ln y F C F T q (; t; u); (B.5) F 6 = F 7 = " ffl + 3 ffl ß 3 +7 ln y F 3lny F CF T g(; t; u); (B.6) tu y F y F ln ffl +ln t +ln yf + Li N C u 4 u ln y F ln t T g (; t; u): u ln + Li yf t (B.7) C Photon Initial State Correction Like gluon, real hoton can lit into q μq air. After integration over the collinear region of hae ace, the Born matrix element for qμq! gg (= T a;b (; t; u)), qq 0! qq 0 (= T (; t; u)), and qq! qq (= T ;3 (; t; u)) can be factorized out. The reult deend on the Mandeltam variable, on the integration cut-off y I, and on the integration variable z a : I = I = I 3 = I 4 = I 5 = I 6 = I 7 = " " " " ffl ffl ffl ffl " ffl P qi ψfl(z a )+ P qi ψfl(z a )ln N C N C P qi ψfl(z a )+ N C P qi ψfl(z a )+ N C P qi ψfl(z a )+ N C N C P qi ψfl(z a )ln N C P qi ψfl(z a )ln N C P qi ψfl(z a )ln P qi ψfl(z a )+ P qi ψfl(z a )ln N C N C + zycl. ermutation of,t, and u; P qi ψfl(z a )+ P qi ψfl(z a )ln N C N C " ffl +(t $ u); P qi ψfl(z a )+ N C +(t $ u): " ffl za y z I a za y z I a za y z I a za y z I a za y z I a za z a y I za P qi ψfl(z a )ln N C z a y I + Q i + Q i + Q i + Q i + Q i + Q i + Q i C F T a(; t; u); (C.) ( N C C F )T b (; t; u);(c.) C F T (; t; u); C F T (t; ; u); C F T 3 (; t; u) N C ( C F )T a (t; ; u) N C T b (t; ; u) The Altarelli-Parii litting function P qi ψfl i equal to N C Q i P qψg (ee Aendix D). (C.3) (C.4) (C.5) (C.6) (C.7) D Parton Initial State Correction In thi aendix, the arton initial tate correction are given a a function of the Mandeltam variable, t, u, of the cut-off y J, and of the additional integration variable : J = ffl C F P qψq ( )+ffi( ) ffl + ffl 6 ln t ln t + ß (D.)