Fragmentation in Jets: Cone and Threshold Effects

Fragmentation in Jets: Cone and Threshold Effects

Massimiliano Procura Albert Einstein Center for Fundamental Physics, Institute for Theoretical Physics,University of Bern, CH-3012 Bern, Switzerland    Wouter J. Waalewijn Department of Physics, University of California at San Diego, La Jolla, CA 92093, U.S.A.

We study the fragmentation of (light) quarks and gluons to hadrons inside a jet of cone size . This allows for a more exclusive analysis of fragmentation than is currently the case. The shape of semi-inclusive cross sections in the hadron energy fraction is described by fragmenting jet functions (FJFs), which we calculate in terms of and the jet energy . We introduce a new joint resummation to sum the double logarithms of and in the FJFs, which has a similar application to initial-state radiation at hadron colliders. Our results at next-to-leading logarithmic order indicate that the resummation of the threshold logarithms of is already important for and improves the convergence of perturbation theory. Our framework may be used to analyze LHC and RHIC data and to test and tune Monte Carlo event generators.


High-energy processes with an observed hadron in the final state can be described by factorizing the short-distance (partonic) physics, which is perturbatively calculable in QCD, from (universal) non-perturbative contributions, see e.g. Ref. Collins:1989gx (). In this context, the (unpolarized) fragmentation functions (FFs)  CSfrag () encode the information on the transition from an energetic parton to a hadron , which carries the fraction of its energy, plus a remainder . The knowledge of both perturbative and non-perturbative ingredients in factorization theorems is crucial to have control on theoretical predictions. For example, a better determination of the -quark FF was important for resolving a discrepancy between CDF data and theory for the spectrum of  Cacciari:2003uh (). Furthermore, additional understanding of hadron production at high in collisions is also required to determine more accurately the relative suppression of hadron spectra (jet quenching) seen in heavy-ion collisions Agakishiev:2011dc ().

Parameterizations for the FFs have been constrained by fitting to data for single-inclusive charged hadron production in at next-to-leading order (NLO) in perturbation theory fragee (); Hirai:2007cx (). More recently, global analyses have been performed that also include semi-inclusive deep-inelastic scattering and/or , data from HERA, RHIC and the Tevatron fragglobal (). To illustrate the current level of precision, the dominant is determined with uncertainties at the 10% level for  Hirai:2007cx (). The FFs of the gluon and the non-valence quarks are known even less accurately.

In contrast to the inclusive analyses listed above, the Belle collaboration is studying light-quark fragmentation in their on-resonance data, using a cut on the thrust event shape to remove the large -quark background Seidl:2008xc (). Here, the fragmentation takes place inside a (hemisphere) jet of invariant mass , described by the fragmenting jet functions (FJFs) that we introduced in Ref. Procura:2009vm (). Correspondingly, should not be too small, to avoid contributions from hadrons outside the jets. For , the FJFs can be perturbatively matched onto the FFs Procura:2009vm (), i.e. where the convolution is in the momentum fraction . In performing this OPE, we assume that parton and hadron masses 111The corrections due to the hadron mass are of order , where is the scale of the hard interaction. For light hadrons they are negligible unless is small. are negligible compared to the jet mass. The Wilson coefficients describe the emissions from the parent parton at larger virtualities building up the jet while encodes physics at lower scales where hadronization effects are important.

The one-loop were determined in Refs. Liu:2010ng (); Jain:2011xz (), and allowed us to calculate in Ref. Jain:2011xz () the cross section for with a cut on thrust up to next-to-next-to-leading logarithmic order in Soft-Collinear Effective Theory (SCET)SCET (). There we found correlations between the thrust cut and that are crucial for the analysis of the Belle data. Of course this analysis is subject to non-perturbative power corrections, which can be sizable at the Belle energy, and a comparison with data is desirable to assess their importance.

Here we study fragmentation inside cone jets, defined with a cone jet algorithm. For simplicity we discuss the case of jets, but our method can also be applied to collisions 222For collisions the jet algorithms are defined in terms of and the pseudo-rapidity and thus not azimuthally symmetric around the jet axis, introducing an additional hurdle. We assume the cancellation of Glauber gluons in factorization theorems.. Following Refs. Ellis:2009wj (); Ellis:2010rwa (); Jouttenus:2009ns (), the cone size is denoted by and a cutoff is applied on the energy in the region between jets 333 may also be defined as a jet energy veto, which is equivalent up to .. The jets are required to be energetic and well-separated:


where is the minimum angular separation between jets and the minimum jet energy. Additional power corrections due to the jet algorithm are suppressed if is not too small DMS (), see also Eq. (Introduction.). SCET can then be applied to separate hard, collinear and soft dynamics, leading to factorization formulae as in Refs. Ellis:2009wj (); Ellis:2010rwa (). Schematically, at leading power,


The hard function describes the hard collision and the soft function the soft radiation (both are matrices in color space). For each of the jets there is a jet function describing the final-state collinear radiation, which is replaced by a FJF for the jet in which the hadron is observed Procura:2009vm (). The phase-space is denoted by . The dependence on the renormalization scale cancels in the cross section up to the order that one is working and is used to estimate the uncertainty from higher-order corrections. Large logarithms of ratios of mass scales are summed by evaluating each of these functions at their natural scale, where they contain no large logarithms, and by running them to a common scale using their renormalization group equations. The FJFs with a cone restriction depend on and rather than , and in the next section we calculate their matching onto FFs at NLO,


To avoid large non-perturbative corrections, should thus not be too small 444From Ref. DMS () non-perturbative effects are expected to scale like for small , which is consistent with the leading non-perturbative correction to the soft function.. Since the other ingredients of the factorization theorem in Eq. (2) do not affect the fragmentation variable , the shape of the cross section in is completely determined by the cone FJF. (This is not the case when invariant masses are measured, as in Ref. Jain:2011xz (), because invariant masses also receive a contribution from soft radiation.) The soft function is sensitive to two different scales Ellis:2010rwa () leading to non-global logarithms ngl () but this only affects the normalization of the cross section. Furthermore, we would like to mention that the generalization of our framework to transverse-momentum-dependent and polarized FFs, can be applied as a tool to obtain information on the proton spin structure by studying azimuthal asymmetries for hadron distributions inside a high- jet in transversely polarized collisions spin ().

Figure 1: The FJF for at LL and NLL order for GeV and . The bands show the perturbative uncertainties from varying the scale in the matching onto up and down by factors of two, and then evolving to . For reference we include the LO result (where is the jet energy ) in the left panel. In the middle panel the same curves and bands relative to this LO result are shown. The dotted lines correspond to the central values without threshold resummation, for which the uncertainties are shown in the right panel.


We now calculate the one-loop matching coefficients for the cone FJFs onto FFs. At one-loop the cone restriction is equivalent to


where is the invariant mass of the jet, the fragmentation variable, the jet energy and the cone size. [Note that and appear in the combination which is invariant under boosts along the jet axis.] We may therefore obtain the matching coefficients from the bare results for the standard FJFs in Ref. Jain:2011xz (). We find, using the scheme,


with the splitting functions in the convention of Eq. (3.7) in Ref. Jain:2011iu (). Anti-quarks have the same coefficients as quarks, and and only start at two-loop order. In Eq. (Calculation.)




We have cross-checked these results with expressions for the jet functions with a cone restriction in Ref. Ellis:2010rwa () (there called “unmeasured jet functions”). The FJFs and jet functions have the same anomalous dimension Procura:2009vm (), consistent with Eq. (2), which provides a check on the UV-divergent terms that we encounter in deriving Eq. (Calculation.). The finite terms were checked with the momentum sum rule Procura:2009vm (); Jain:2011xz ()


and a new sum rule


This new relation follows from quark-number conservation in the perturbative calculation of the FJF .

Figure 2: The FJF for is shown at NLL order for GeV. The curves are for and are shown relative to the result.

Threshold Logarithms.

Closer inspection reveals that and contain large threshold logarithms. For example,


which by Eq. (Introduction.) leads to the large (double) logarithms


We may thus sum both the logarithms of and the threshold logarithms by evaluating the FJF at . (In the plots we evolve the FJF to , to remove the dependence from this scale.) Fig. 1 shows the improved convergence of perturbation theory arising from threshold resummation, which will be discussed more extensively in the next section. We stress that we do not need to perform a threshold expansion; therefore our results are also valid away from the threshold region.

This improvement holds to all orders in perturbation theory, since we find that the logarithms of and are tied together through the anomalous dimension of 555This requires the all-orders structure of the jet function anomalous dimension, which is known for hemisphere invariant mass jets, but not established for jet algorithms.. Using Eq. (Introduction.) this can be obtained from the anomalous dimension of the (unmeasured) jet function Ellis:2009wj (); Ellis:2010rwa () and the FF Korchemsky:1992xv (),


Here is the (quark) cusp anomalous dimension. The non-cusp anomalous dimension only starts at two loops, which is reflected in the absence of single logarithms in Eqs. (Threshold Logarithms.) and (11). In moment space one can see explicitly how the logarithm of and the combine, since the th moment of is .

The above discussion applies to as well. By contrast, the off-diagonal and do not contain threshold logarithms and their contribution should thus be evaluated at the appropriate scale .

Our new joint resummation can also be used in the standard FJFs of Refs. Procura:2009vm (); Liu:2010ng (); Jain:2011xz (). Here the scale choice sums both the logarithms of and . We have checked that this significantly improves the convergence of perturbation theory at large .

Similarly, for initial-state jets described by beam functions  beamf (), we find that the logarithms of and can be simultaneously resummed by evaluating the flavor-diagonal terms for matching beam functions onto PDFs at the scale . At the LHC, the momentum fraction is typically small, suggesting that threshold resummation is less important than in the case studied here. However, it has been argued that the logarithms of are dynamically enhanced through the shape of the PDFs thr (). A detailed study is left for future work.

LO -loop - - -loop
NLO -loop - - -loop
LL -loop - -loop -loop
NLL -loop -loop -loop -loop
Table 1: Order counting for the matching coefficients, the non-cusp and cusp anomalous dimension, and the running.


We now show some numerical results for fragmentation inside a cone jet. Our order counting is shown in Table 1. As input we use the HKNS FFs Hirai:2007cx () at LO or NLO order and their value for .

We first consider a fixed cone and study the convergence of our resummed result. Fig. 1 shows the central values and perturbative uncertainties for the FJF at LL and NLL order. In the middle and right panels these results relative to the LO curve are shown, making the separation between the central curves and the size of the uncertainties more visible. As expected, the cone restriction suppresses the small region because at large the hadron carries off most of the energy and there is not much additional radiation that would be subject to the cone restriction. As is clear from comparing the dotted with the solid curves, the effect of threshold resummation becomes important at . Including threshold resummation improves the convergence and leads to smaller uncertainties. At small the difference between LL and NLL is not fully captured by the uncertainty band. However, in this region the uncertainties on the FFs (which are not shown here) dominate.

We have studied similar plots for and . Here the convergence is not as good when the dominant mixing contribution involving (and ) turns on. For and this happens at NLO and NNLO, respectively.

In Fig. 2 we study the effect of varying the cone size . We stress that the normalization of the curves for different should not be compared because other ingredients in the factorization theorem affect the normalization of the cross section (the shape in is entirely determined by ). For and the small region gets suppressed relative to the large region as one reduces the cone size, as expected. For , the bump-shaped enhancement at intermediate is mainly due to the contribution.


Our framework leads to a reliable description of fragmentation within an identified jet, providing the tools for a more exclusive study of fragmentation. This can be used to reduce the background from other processes, to help identify underlying partonic structures and enables novel tests of the universality of FFs. Our setup provides accurate analytical predictions, where large logarithms are properly summed, for single hadron spectra in a jet which can be used to tune Monte Carlo event generators. This will become more and more relevant as the determination of the fragmentation functions becomes more accurate. We have for the first time established the connection between the resummation of logarithms for exclusive jet production and threshold resummation, and have shown the numerical importance of their interplay in the case of jet fragmentation.

We thank T. Becher, S. Ellis, A. Jain, C. Lee, A. Manohar, J. Walsh and F. Würthwein for discussions and comments. M.P. acknowledges support by the “Innovations- und Kooperationsprojekt C-13” of the Schweizerische Universitätskonferenz SUK/CRUS and by the Swiss National Science Foundation. W.W. is supported by DOE grant DE-FG02-90ER40546.


  • (1) J. C. Collins, D. E. Soper and G. F. Sterman, Adv. Ser. Direct. High Energy Phys. 5, 1 (1988).
  • (2) J. C. Collins, and D. E. Soper, Nucl. Phys. B 193, 381 (1981); ibid. 194, 445 (1982).
  • (3) M. Cacciari et al., JHEP 07, 033 (2004).
  • (4) G. Agakishiev et al., STAR Coll., Phys. Rev. Lett. 108, 072302 (2012).
  • (5) S. Albino, B. A. Kniehl and G. Kramer, Nucl. Phys. B 725, 181 (2005); 734, 50 (2006).
  • (6) M. Hirai et al., Phys. Rev. D 75, 094009 (2007).
  • (7) D. de Florian, R. Sassot and M.  Stratmann, Phys. Rev. D 75, 114010 (2007); ibid. 76, 074033 (2007); S. Albino, B. A. Kniehl, G. Kramer, Nucl. Phys. B 803, 42 (2008).
  • (8) R. Seidl et al., Phys. Rev. D 78, 032011 (2008).
  • (9) M. Procura and I. W. Stewart, Phys. Rev. D 81, 074009 (2010).
  • (10) X. Liu, Phys. Lett. B 699, 87 (2011).
  • (11) A. Jain, M. Procura and W. J. Waalewijn, JHEP 05, 035 (2011).
  • (12) C. W. Bauer, S. Fleming and M. E. Luke, Phys. Rev. D 63, 014006 (2000); C. W. Bauer et al., Phys. Rev. D 63, 114020 (2001); C. W. Bauer and I. W. Stewart, Phys. Lett. B 516, 134 (2001); C. W. Bauer, D. Pirjol and I. W. Stewart, Phys. Rev. D 65, 054022 (2002).
  • (13) S. D. Ellis et al., Phys. Lett. B 689, 82 (2010).
  • (14) S. D. Ellis et al., JHEP 11, 101 (2010).
  • (15) T. T. Jouttenus, Phys. Rev. D 81, 094017 (2010).
  • (16) M. Dasgupta, L. Magnea and G. P. Salam, JHEP 02, 055 (2008).
  • (17) M. Dasgupta and G. P. Salam, Phys. Lett. B 512, 323 (2001); JHEP 08, 032 (2002).
  • (18) F. Yuan, Phys. Rev. Lett. 100, 032003 (2008); Phys. Rev. D 77, 074019 (2008); R. Fersch, J. Phys. Conf. Ser. 295, 012048 (2011).
  • (19) A. Jain et al., JHEP 04, 132 (2012).
  • (20) G. P. Korchemsky and G. Marchesini, Nucl. Phys. B 406 225 (1993).
  • (21) S. Fleming, A. K. Leibovich and T. Mehen, Phys. Rev. D 74, 114004 (2006); I. W. Stewart, F. J. Tackmann and W. J. Waalewijn, Phys. Rev. D 81, 094035 (2010); S. Mantry and F. Petriello, Phys. Rev. D 81, 093007 (2010).
  • (22) D. Appell, G. F. Sterman and P. B. Mackenzie, Nucl. Phys. B 309 259 (1988); S. Catani, and D. de Florian and M. Grazzini JHEP 05, 025 (2001); T. Becher, M. Neubert and G. Xu, JHEP 07, 030 (2008).
Comments 0
Request Comment
You are adding the first comment!
How to quickly get a good reply:
  • Give credit where it’s due by listing out the positive aspects of a paper before getting into which changes should be made.
  • Be specific in your critique, and provide supporting evidence with appropriate references to substantiate general statements.
  • Your comment should inspire ideas to flow and help the author improves the paper.

The better we are at sharing our knowledge with each other, the faster we move forward.
The feedback must be of minimum 40 characters and the title a minimum of 5 characters
Add comment
Loading ...
This is a comment super asjknd jkasnjk adsnkj
The feedback must be of minumum 40 characters
The feedback must be of minumum 40 characters

You are asking your first question!
How to quickly get a good answer:
  • Keep your question short and to the point
  • Check for grammar or spelling errors.
  • Phrase it like a question
Test description