Fragmentation to a jet in the large z limit Lin Dai,1,∗ Chul Kim,2,† and Adam K. Leibovich1,‡ 1Pittsburgh Particle Physics Astrophysics and Cosmology Center (PITT PACC) Department of Physics and Astronomy, University of Pittsburgh, Pittsburgh, Pennsylvania 15260, USA 2Institute of Convergence Fundamental Studies and School of Liberal Arts, 7 1 Seoul National University of Science and Technology, Seoul 01811, Korea 0 2 r a Abstract M We consider the fragmentation of a parton into a jet with small radius R in the large z limit, 0 wherez istheratioofthejetenergytothemotherpartonenergy. Inthisregionofphasespace,large 3 logarithms of both R and 1 z can appear, requiring resummation in order to have a well defined − ] perturbative expansion. Using soft-collinear effective theory, we study the fragmentation function h p toa jet(FFJ) inthis endpoint region. We derivea factorizationtheorem forthis object, separating - collinear and collinear-soft modes. This allows for the resummation using renormalization group p e evolution of the logarithms lnR and ln(1 z) simultaneously. We show results valid to next-to- h − leading logarithmic order for the global Sudakov logarithms. We also discuss the possibility of [ non-global logarithms that should appear at two-loops and give an estimate of their size. 2 v 0 6 6 5 0 . 1 0 7 1 : v i X r a E-mail:[email protected] ∗ E-mail:[email protected] † E-mail:[email protected] ‡ 1 I. INTRODUCTION The fragmentation function (FF) [1], which describes an energetic splitting of a parton into a final state, is a very important ingredient in understanding high-energy hadron pro- duction. Using the FF we can systematically separate short- and long-distance interactions related to the production. For instance, inclusive hadron production for e+e− annihilation can be factorized as dσ(e+e− hX) (cid:90) 1 dzdσ (z /z,µ) i h → = D (z,µ), (1) h/i dE z dE h z i h where i denotes the flavor of the produced parton, z = 2E /E , and z = E /E . Here, h h cm h i E is the center of the mass energy of the collision. The partonic scattering cross section cm σ includes the hard interactions for e+e− iX. Long-distance interactions describing the i → fragmenting process from parton i to hadron h are encoded in the FF, D (z). The FF h/i is universal in the sense that it is independent of the hard process and can be applied to other scattering processes. Hence, the FF has long been studied in order to understand its properties. (For details we refer to a recent review [2] and the references therein.) Because we can directly observe a jet using well-defined jet algorithms such as the ones introduced in Refs. [3–7], it is possible to describe the fragmentation function to a jet (FFJ), as long as the the jet radius, R, is enough small [8]. (For a recent review of jet physics see, for example, [9].) Moreover, once the FFJ for the isolated jet is given, we can systematically investigate its substructures (e.g. hadron and subjet fragmentations [10–15], andjetmass[16]andtransversemomentum[17,18]distributions),constructingfactorization theorems in connection with the frgmenting jet functions [19–21]. Analytical results of the FFJ have been calculated up to the next-to-leading order (NLO) in α [12, 14, 22]. Unlike the hadron FF, the FFJ does not have any infrared (IR) divergence s due to the finite size of the jet radius R. However, the presence of large logarithms of R does not give a reliable result in perturbation theory and requires resummation to all order in α . s As shown in Refs. [8, 12, 14, 22], resumming logarithms of R is equivalent to running down to a scale µ QR using Dokshitzer-Gribov-Lipatov-Altarelli-Parisi (DGLAP) evolution ∼ equations, where Q is a hard energy comparable to the jet energy, E . This resummed J result of the FFJ has been successfully applied to inclusive jet [22, 23] and hadron [15] production, where the effects of various values of R have been investigated in detail. If we observe a highly energetic jet, we would expect that most of the energetic splitting processes are captured within the jet radius R since these processes favor small angle ra- diation. This implies that the large z region gives the dominant contribution to the FFJ, where z is the ratio of the jet energy fraction over the mother parton energy. Accordingly, in the perturbative result for the FFJ there are large logarithms of 1 z, which need to be − resummed to all order in α . Already at one loop order there appears a double logarithm s ln(1 z)/(1 z) L2, where L schematic represents a large logarithm. At leading loga- + − − ∼ 2 rithm(LL)accuracy,theresummedcanberepresentedas(cid:80) C (α L2)k exp(Lf (α L)), k=0 k s ∼ 0 s which gives the dominant correction to the perturbative expansion of the FFJ. Thus, for a proper description of the FFJ in the large z limit, we have to systematically handle large logarithms of 1 z as well as large logarithms of R. In general, if some quantity − involves several distinct scales, we try to factorize it so that each factorized part can be well described at one properly chosen scale. Then performing evolutions between these largely separated scales, we resum the large logarithms. For the FFJ, soft-collinear effective theory (SCET) [24–27] provides the appropriate framework for factorization and enable us to resum large logarithms automatically by solving the renormalization group (RG) equations for the factorized parts. Near the endpoint where z 1, the FFJ consists of dynamics with two well-separated → scales. Since an observed jet carries most of energy of the mother parton, radiation outside the jet should be soft with energy E (1 z). Therefore the jet splitting process can be J ∼ − initiated by soft dynamics, while radiation inside the jet is described dominantly by collinear interactions. However, in the effective theory approach wide angle soft interactions are not adequate for explaining the radiation outside the narrow jet because they cannot effectively recognize the jet boundary characterized by the small radius R. Instead, we introduce a more refined soft mode, namely the collinear-soft mode [28, 29], which can resolve the narrow jet boundary and can consistently describe the lower energy, out-of-jet radiations. The collinear-soft mode has previously been used to factorize the cross sections for a narrow jet at a low energy scale [30–33]. In this paper, using SCET we construct a factorization theorem for the FFJ near the endpoint considering collinear and collinear-soft interactions.1 Then we resum the large logarithms of 1 z and R simultaneously. In sec. II we discuss the characteristics of large-z − physics for the FFJ and factorize it into the collinear and the collinear-soft pieces. Then, we confirm our factorized result through NLO by an explicit calculation of each factorized part. In sec. III, based on the factorization, we resum the large logarithms by performing RG evolution. We also discuss large nonglobal logarithms (NGLs) that possibly contribute to NLL accuracy. In sec. IV the numerical results of the FFJ to the accuracy of NLL plus NLO in α are shown. Finally in sec. V we conclude. s 1 In a strict sense our factorization theorem would hold up to NLO in α . Beyond NLO, large nonglobal s logarithms (NGLs) [34, 35] that are sensitive to a restricted jet phase space might appear and require some modification of our factorization theorem presented here. 3 II. THE FFJ IN THE LIMIT z 1 → Using SCET, the FFJ can be defined as [14] (cid:88) zD−3 (cid:16)p+ (cid:17)n/ D (z,µ) = Tr 0 δ J Ψ J (p+,R)X J (p+,R)X Ψ¯ 0 ,(2) Jk/q 2N (cid:104) | z −P+ 2 n| k J ∈/J(cid:105)(cid:104) k J ∈/J| n| (cid:105) c X∈/J,XJ−1 (cid:88) zD−3 D (z,µ) = (3) Jk/g p+(D 2)(N2 1) X∈/J,XJ−1 J − c − (cid:16)p+ (cid:17) Tr 0 δ J ⊥µ,a J (p+,R)X J (p+,R)X ⊥a 0 × (cid:104) | z −P+ Bn | k J ∈/J(cid:105)(cid:104) k J ∈/J|Bnµ| (cid:105). Here Ψ = W†ξ and ⊥µ,a = inρgµνGb ba = inρgµν †,baGb are gauge invariant n n n Bn ⊥ n,ρνWn ⊥ Wn n,ρν collinear quark and gluon field strength respectively. W ( ) is a collinear Wilson line n n W in the fundamental (adjoint) representation [25, 26]. These collinear fields have momentum scaling pµ = (p ,p ,p ) = Q(1,λ,λ2), where λ is a small parameter comparable to small jet n + ⊥ − radius R. p are denoted as p n p = p +nˆ p and p n p = p nˆ p, where nˆ is ± + 0 J − 0 J J ≡ · · ≡ · − · a unit vector in the jet direction and two lightcone vectors nµ = (1,nˆ ) and nµ = (1, nˆ ) J J − have been employed. The expressions for the FFJs in Eqs. (2) and (3) are valid in the jet frame where the transverse momentum of the observed jet, p⊥, is zero. J In this paper, we will consider inclusive k -type algorithms [3–5, 7], where the merging T condition of two light particles is given by θ < R(cid:48). (4) Here θ is the angle between the two particles, and R(cid:48) = R for an e+e− collider and R(cid:48) = R/coshy for a hadron collider, where y (1) is the rapidity for the central region. ∼ O The definitions of the FFJs in Eqs. (2) and (3) hold for z (1), but are not reliable ∼ O near the endpoint where z goes to 1. In the limit z 1, the observed jet takes most of the → energy from the mother parton and hence the jet splitting (out-jet) contributions should be described by soft gluon radiation. If 1 z is power counted as (η) with η 1, the relevant − O (cid:28) soft mode would have momentum scaling k (k ,k ,k ) Q(η,η,η). However, for the + ⊥ − ∼ ∼ proper resummation of lnR, we need a mode that can probe the jet boundary expressed in terms of R. This mode would have a lower resolution than the soft mode while the k + component should still be power counted as (η). Because the jet merging criterion for the O soft gluon radiation is given by [36] R(cid:48) k tan2 > − , (5) 2 k + the proper mode should allow for the hierarchy, k k λ2 k , where λ R. Thus − + + ∼ (cid:28) ∼ this mode should have scaling k Qη(1,λ,λ2). From now on we will call this mode the ∼ collinear-soft mode. 4 We can consistently separate the usual soft mode Q(η,η,η) and the collinear-soft mode ∼ as was first done in the di-jet scattering cross section [30, 31]. Furthermore, the separation of the collinear-soft mode from the collinear fields has been performed in the formulation of SCET [28]. Because the collinear-soft mode can be considered as a subset of the usual soft + mode, we have to subtract the overlapped of the collinear-soft contribution from the soft contribution in loop calculations similar to the usual zero-bin subtractions [37]. If we apply this process to the FFJ with z 1, we see that the soft contributions → can be cancelled by the collinear-soft subtractions. Since the soft mode with a scaling (k ,k ) Q(η,η) cannot resolve the jet boundary in Eq. (5), the real soft gluon radiation + − ∼ does not contribute to the in-jet contribution of the FFJ, while the out-jet contribution from real radiation covers the full phase space of (k ,k ). Thus, independent of R, the + − total soft contributions will be expressed as a function of 1 z, namely S(1 z). For the − − collinear-soft contribution that needs to be subtracted from the soft contribution, we apply the same boundary conditions used for the soft mode. Hence the real collinear-soft radiation have only the out-jet contributions, which are the same as the soft mode. Therefore the net result of the collinear-soft contributions that are to be subtracted are the same as S(1 z), − canceling the soft contribution. Finally we are left with a collinear-soft mode at the lower energy scale. When we apply thistotheFFJ,wehavetokeepthejetboundaryconstraintinEq.(5). Asaresulttheactive collinear-soft contributions can be expressed in terms of 1 z and R simultaneously. As − we will see, the one loop collinear-soft contributions involve double logarithms of lnµ/((1 − z)E R(cid:48)). This fact indicates that the collinear-soft interactions are responsible for large J logarithms of 1 z and its resummation would give the dominant contribution to the FFJ − near the endpoint. A. Factorization of the FFJ when z 1 → With the above reasoning, we can systematically extend the FFJs to the endpoint region including collinear-soft interactions. We first decouple the soft mode Q(η,η,η) from the ∼ collinear mode Q(1,R,R2). Then we introduce the collinear-soft mode Qη(1,R,R2) ∼ ∼ in the collinear sector, classifying collinear and collinear-soft gluons as Aµ Aµ + Aµ . n → n n,cs Accordingly the covariant derivative in the collinear sector decomposes as iDµ = iDµ + c iDµ = µ+gAµ+i∂µ+gAµ , where µ (i∂µ) returns collinear (collinear-soft) momentum. cs P n n,cs P In this decomposition, the commutation relations, [ µ,Aν ] = [∂µ,Aν] = 0, hold. For the P n,cs n factorization of the FFJ, our strategy is simple: after the decomposition into the collinear and collinear-soft modes, we first integrate out collinear interactions with p2 Q2R2. As we c ∼ shall see, this givesan integratedjet functioninside ajet. Then atthe lower scaleµ QηR cs ∼ we will consider the collinear-soft interactions for the jet splitting. As performed in Ref. [28], at low energy we can additionally introduce so called ‘ultra- 5 collinear’ modes after integrating out the collinear interactions with offshellness p2 Q2R2. c ∼ These modes have energy of the same order as the collinear mode, but their fluctuations are much smaller than Q2R2. Then at the low energy scale an external collinear field φ(= ξ,A) n would be matched onto the ultra-collinear fields, φ = φ +φ + , where the lightcone n n1 n2 ··· vectors n reside inside the jet with radius R. Note that collinear interactions between i=1,2,··· different ultra-collinear modes are forbidden since we have already integrated out the large collinear fluctuations Q2R2. Moreover, as these ultra-collinear modes reside within the ∼ collinear interactions, they cannot resolve the jet boundary. Therefore their interactions do not contribute to the FFJs, at least to NLO in α . So for simplicity we will not consider s ultra-collinear interactions in the FFJ. However, in a more refined jet observable identifying subjets, these modes may have to be included. Addingthecollinear-softmode, thequarkinitiatedFFJcanbemoregenericallyexpressed as (cid:88) zD−3 (cid:16)p+ (cid:17)n/ D (z,µ) = Tr 0 δ J n iD ξ J (p+,R)X J (p+,R)X ξ¯ 0 . (6) Jk/q 2N (cid:104) | z − · 2 n| k J ∈/J(cid:105)(cid:104) k J ∈/J| n| (cid:105) c X∈/J,XJ−1 Compared to Eq. (2), W δ(p+/z )W† = δ(p+/z n iD ) has been replaced with n J − P+ n J − · c δ(p+/z n iD) in Eq. (6). J − · In order to satisfy gauge invariances at each order in λ (R) and η, following the ∼ O procedure considered in Ref. [38], we redefine the collinear gluon field, Aµ = Aˆµ +Wˆ [iDµ,Wˆ †] (7) n n n cs n , ˆ ˆ where A are newly defined collinear gluon fields and W is the collinear Wilson line ex- n n ˆ pressed in terms of A . As a consequence the covariant derivative in Eq. (6) can be rewritten n as iDµ = iDµ +W iDµW†, (8) c n cs n where collinear fields on the right-hand side are the redefined fields and we removed the hat for simplicity. Employing Eq. (8), the delta function in Eq. (6) can be rewritten as (cid:16)p+ (cid:17) (cid:16)p+ (cid:17) δ J n iD = W δ J n iD W† (9) z − · n z −P+ − · cs n . Similar to the decoupling of leading ultrasoft interactions from collinear fields [26], we can remove collinear-soft interactions through the term gn A in the Lagrangian of the cs · collinear sector. To accomplish this, the collinear quark and gluon fields can be additionally redefined as ξ Ycsξ , Aµ YcsAµYcs†, (10) n → n n n → n n n where Ycs is the collinear-soft Wilson line that satisfies n iD Ycs = Ycsn i∂ and has the n · cs n n · usual form [26, 39] (cid:34) (cid:35) (cid:90) ∞ Ycs(x) = P exp ig dsn A (sn) . (11) n · cs x 6 Using Eqs. (8) and (10) we rewrite Eq. (6) as (cid:88) zD−3 (cid:16)p+ (cid:17)n/ D (z,µ) = Tr 0 δ J i∂ Ycs†YcsW†ξ J (p+,R)X Jk/q 2N (cid:104) | z −P+ − + 2 n n n n| k J ∈/J(cid:105) c X∈/J,XJ−1 J (p+,R)X ξ¯ W Ycs†Ycs 0 (12) ×(cid:104) k J ∈/J| n n n n | (cid:105), where we used the relation n iD = Ycsi∂ Ycs† and Ycs has the same form as Eq. (11) with · cs n + n n replacement of n n. We also used the crossing symmetry φ X = X φ, where φ φ → ···| (cid:105) (cid:104) |··· φ = W , Ycs. The FFJ in Eq. (12) can describe regions of ordinary z (1) and z 1. n n ∼ O → If z is ordinary and not too close to 1, we can suppress i∂ in the argument of the delta + function,sincep+/z (Q)ispowercountedmuchlargerthani∂ (Qη). Thusthe J −P+ ∼ O + ∼ O collinear-soft Wilson lines cancel by unitarity and we recover the form in Eq. (2). However, when z 1, p+/z becomes the same size as i∂ , and we cannot ignore the term i∂ → J −P+ + + in the delta function, which gives nonzero contributions of collinear-soft interactions. Since returns collinear (label) momentum in Eq. (12), can be fixed as p+ near P+ P+ J the endpoint. Further, it means that collinear interactions are relevant only for jet merging (in-jet) contribution to the FFJ. Therefore the FFJ in the limit z 1 can be expressed as2 → (cid:88) zD−3 n/ D (z 1,µ) = Tr 0 Ycs†Ycs W†ξ J (p+,R)X Jq/q → 2N (cid:104) | n n 2 n n| q J ∈/J(cid:105) c X∈/J,XJ−1 (cid:16)p+ (cid:17) J (p+,R)X ξ¯ W δ J † +i∂ Ycs†Ycs 0 ×(cid:104) q J ∈/J| n n z −P+ + n n | (cid:105) (cid:88) 1 n/ (cid:88) 1 = Tr 0 W†ξ qX J qX J ξ¯ W 0 Tr 0 Ycs†Ycs X 2N (cid:104) |2 n n| c ∈ (cid:105)(cid:104) c ∈ | n n| (cid:105)· N (cid:104) | n n | cs(cid:105) c c Xc∈J Xcs (cid:0) (cid:1) X δ (1 z)p+ +Θ(θ R(cid:48))i∂ Ycs†Ycs 0 (13) ×(cid:104) cs| − J − + n n | (cid:105), where Θ is the step function and we reorganized the final states into collinear states (qX ) c in the jet and collinear-soft states X in order to factorize collinear and collinear-soft in- cs teractions. In the second equality we fixed the collinear label momentum † as p+, and P J then we put the jet splitting constraint in front of i∂ because only the out-jet collinear-soft + radiation gives a nonzero contribution for the region z < 1. From Eq. (5), the jet split- ting constraint Θ(θ R(cid:48)) is equivalent to tan2R(cid:48)/2 < k /k , where k is the collinear-soft − + − momentum. Eq. (13) shows that the quark FFJ in the limit z 1 is factorized as → D (z 1,µ;E R(cid:48),(1 z)E R(cid:48)) = (µ;E R(cid:48),θ < R(cid:48)) S (z,µ;(1 z)E R(cid:48)), (14) Jq/q → J − J Jq J · q − J where is the integrated jet function for the in-jet contribution, defined as q J (cid:88) 1 n/ (µ;E R(cid:48),θ < R(cid:48)) = Tr 0 W†ξ qX J(E ,R(cid:48)) qX J ξ¯ W 0 (15) Jq J 2N p+ (cid:104) |2 n n| c ∈ J (cid:105)(cid:104) c ∈ | n n| (cid:105). Xc∈J c J 2 Note that the splitting q J in the limit z 1 is power suppressed by (1 z) compared to the g → → O − splitting q J . For q J , the splitted parton away from the observed jet is the collinear-soft quark, q g → → which gives a power suppression of (η) compared to the collinear-soft gluon radiation. Similarly, for O gluon splitting, g J dominants for the same reason. g → 7 S isthedimensionlesscollinear-softfunction. WhenwerewriteS = p+S˜ , thedimensionful q q J q ˜ collinear-soft function S can be expressed as q S˜ ((cid:96) ,µ;(cid:96) t) = (cid:88) 1 Tr 0 Ycs†Ycs X X δ(cid:0)(cid:96) +Θ(θ R(cid:48))i∂ (cid:1)Ycs†Ycs 0 (16) q + + N (cid:104) | n n | cs(cid:105)(cid:104) cs| + − + n n | (cid:105), c Xcs where t tanR(cid:48)/2, and (cid:96) t is the scale that will minimize large logarithms in the higher + ≡ order corrections, as we will see later. Using the adjoint representation and taking a similar procedure as we did with the quark case, we obtain the factorization formula for the gluon FFJ, D (z 1,µ) = (µ;E R(cid:48),θ < R(cid:48)) S (z,µ;(1 z)E R(cid:48)), (17) Jg/g → Jg J · g − J where is the gluon integrated jet function, and S is the collinear-soft function defined g g J similar to Eq. (16), with the Wilson lines in the adjoint representation replacing Ycs. n,n B. NLO calculation of the FFJ near the endpoint The integrated jet functions shown in Eqs. (14) and (17) have been explicitly computed at NLO [36, 40, 41] and partially computed at NNLO [31, 32]. The NLO results with the constraint of Eq. (4) read (cid:34) α C 1 1 (cid:16)3 µ2 (cid:17) (µ;E R(cid:48),θ < R(cid:48)) = 1+ s F + +ln Jq J 2π (cid:15)2 (cid:15) 2 p+2t2 UV UV J (cid:35) 3 µ2 1 µ2 13 3π2 + ln + ln2 + , (18) 2 p+2t2 2 p+2t2 2 − 4 J J (cid:34) α C 1 1 (cid:16) β µ2 (cid:17) β µ2 (µ;E R(cid:48),θ < R(cid:48)) = 1+ s A + 0 +ln + 0 ln Jg J 2π (cid:15)2 (cid:15) 2C p+2t2 2C p+2t2 UV UV A J A J (cid:35) 1 µ2 67 23n 3π2 + ln2 + f , (19) 2 p+2t2 9 − 18C − 4 J A where p+t E R(cid:48), β = 11N /3 2n /3, C = N = 3, and n is the number of flavors. J ∼ J 0 c − f A c f For the NLO computation of the collinear-soft function in Eq. (16) we consider virtual and real gluon contributions respectively. Separating ultraviolet (UV) and infrared (IR) divergences carefully, the virtual contributions are given by α C (cid:16) 1 1 (cid:17)2 MS = s F δ((cid:96) ). (20) V − π (cid:15) − (cid:15) + UV IR The real contributions at one loop can be written as α C (µ2eγE)(cid:15) (cid:90) ∞ (cid:104) MS = s F dk dk (k k )−1−(cid:15) δ((cid:96) k )Θ(k t2k ) R π Γ(1 (cid:15)) + − + − + − + − − + 0 − (cid:105) +δ((cid:96) )Θ(t2k k ) MS +MS , (21) + + − − ≡ R1 R2 8 out jet (✓ > R ) 0 � ` ⇤ + + k 2 k = t k + � � k + in jet (✓ < R ) 0 � FIG. 1. Phase space for the real gluon emission in the collinear-soft function. In the (k ,k ) + − plane, the region above the border line k = t2k gives the out-jet contribution and the region − + below gives the in-jet contribution. Λ is the maximum value for the distribution of (cid:96) and can + + be chosen arbitrarily. where k is the momentum of the outgoing collinear-soft gluon and MS (MS ) indicates the R1 R2 contribution from the first (second) term in the square brackets. In Fig. 1 we show the possible phase space for the emitted collinear-soft gluon after the integration over k . MS covers region below the jet border line (k = t2k ). Hence the ⊥ R2 − + result is α C (µ2eγE)(cid:15) (cid:90) ∞ (cid:90) t2k+ MS = s F δ((cid:96) ) dk dk (k k )−1−(cid:15) R2 π Γ(1 (cid:15)) + + − + − 0 0 − (cid:34) (cid:35) α C (cid:18) 1 1 (cid:19)2 (cid:18) 1 1 (cid:19) = s F lnt2 δ((cid:96) ) (22) + . 2π ε − ε − ε − ε UV IR UV IR For MS , k is fixed to be (cid:96) by the delta function, and the possible phase space has R1 + + been denoted as a blue line in the upper plane in Fig 1. However we need to extract the IR divergences as (cid:96) 0. In order to do so, we introduce the so called Λ -distribution, which + + → is defined as (cid:90) L (cid:90) L (cid:90) Λ+ d(cid:96) [g((cid:96) )] f((cid:96) ) = d(cid:96) g((cid:96) )f((cid:96) ) d(cid:96) g((cid:96) )f(0), (23) + + Λ+ + + + + − + + 0 0 0 where f((cid:96) ) is an arbitrary smooth function at (cid:96) = 0. Λ is an arbitrary upper limit for + + + Λ -distribution and is power counted to have the same size as (cid:96) . We can write MS using + + R1 9 this distribution, α C (µ2eγE)(cid:15) (cid:90) ∞ MS = s F (cid:96)−1−(cid:15) dk k−1−(cid:15) R1 π Γ(1 (cid:15)) + − − − t2(cid:96)+ (cid:34) (cid:35) α C (µ2eγE)(cid:15) (cid:90) ∞ = δ((cid:96) )I (Λ ,t)+ s F (cid:96)−1−(cid:15) k−1−(cid:15) , (24) + R1 + π Γ(1 (cid:15)) + − − t2(cid:96)+ Λ+ where the integration region for I corresponds to the green region in Fig. 1. Integrating R1 over this region, we get (cid:34) (cid:35) α C (µ2eγE)(cid:15) (cid:90) ∞ (cid:90) ∞ (cid:90) ∞ (cid:90) ∞ I = s F dk dk (k k )−1−(cid:15) dk dk (k k )−1−(cid:15) R1 + − + − + − + − π Γ(1 (cid:15)) − − 0 t2k+ Λ+ t2k+ (cid:34) α C (cid:18) 1 1 (cid:19)2 (cid:18) 1 1 (cid:19) = s F + lnt2 (25) 2π (cid:15) − (cid:15) (cid:15) − (cid:15) UV IR UV IR (cid:35) (cid:18) 1 1 µ2 1 µ2 π2(cid:19) + ln + ln2 . − (cid:15)2 (cid:15) Λ2t2 2 Λ2t2 − 12 UV UV + + The second term in Eq. (24) is given by (cid:34) (cid:35) (cid:34) (cid:35) α C (µ2eγE)(cid:15) (cid:90) ∞ α C 1 (cid:16) 1 µ2 (cid:17) s F (cid:96)−1−(cid:15) k−1−(cid:15) = s F +ln . (26) π Γ(1 (cid:15)) + − π (cid:96) (cid:15) (cid:96)2t2 − t2(cid:96)+ Λ+ + UV + Λ+ ˜ Finally combining Eqs. (20), (22), (25) and (26) we obtain the bare one loop result of S , q M = MS +MS +MS S V R1 R2 (cid:40) (cid:32) (cid:33) α C 1 1 µ2 1 µ2 π2 = s F δ((cid:96) ) ln ln2 + (27) π + −2(cid:15)2 − 2(cid:15) Λ2t2 − 4 Λ2t2 24 UV UV + + (cid:34) (cid:35) (cid:41) 1 (cid:16) 1 µ2 (cid:17) + +ln . (cid:96) (cid:15) (cid:96)2t2 + UV + Λ+ The one loop result of the collinear-soft function for gluon FFJ is the same if we replace C F with C = N in Eq. (27). A c Since the dimensionless soft-collinear function, S (z) = p+S˜ ((cid:96) ), is a function of z, k=q,g J k + we need to express the Λ -distribution in terms of the standard plus distribution of z. From + Eq. (23) we obtain the relation 1 1 (cid:90) b [g˜((cid:96) )] = [g(z)] + δ(1 z) dz(cid:48)g(z(cid:48)), (28) + Λ+ p+ + p+ − J J 0 where (cid:96) = p+(1 z) and g(z) = p+g˜((cid:96) ). In the Λ -distribution, Λ has been replaced + J − J + + + with p+(1 b), where b is a dimensionless parameter close to 1. J − 10