Medium modified leading hadrons, jets and sub-jets in a single formalism

Medium modified leading hadrons, jets and sub-jets in a single formalism

   Shanshan Cao and Abhijit Majumder
Department of Physics and Astronomy
   Wayne State University    Detroit    Michigan 48201    USA

Vacuum and medium modified fragmentation functions are used to study the properties of hadrons produced in hard interactions, in - and - collisions, respectively. We study the modification of reconstructed jets and high transverse momentum (high-) hadrons as they propagate through a dense medium. Jets and jet modification are studied within the framework of a factorized jet function. The vacuum jet function is found consistent with PYTHIA simulations. Both the modification of the fragmentation function and the jet function are carried out within the identical higher-twist energy loss formalism. This approach is then extended to di-hadron production within a jet. Its further extension to a “di-prong” function and its application to the soft-drop measurements are discussed. The methodology of extending these semi-analytical calculations to full Monte-Carlo simulations is outlined. All calculations are carried out within a (2+1)-D viscous hydrodynamic simulation and compared to experimental data from a variety of collision energies.


Medium modified leading hadrons, jets and sub-jets in a single formalism \FullConferenceInternational Conference on Hard and Electromagnetic Probes of High-Energy Nuclear Collisions
30 September - 5 October 2018
Aix-Les-Bains, Savoie, France

1 Introduction

Jets, collimated sprays of particles, are produced in the early stage of relativistic heavy-ion collisions. They gather information about the Quark-Gluon Plasma (QGP) as they traverse it [1]. Usually, jets are defined according to their energy and radius. Small radius jets are widely studied in heavy ion collisions which are expected to have little contamination from the QGP background. A jet function is usually defined to study small radius jets in theoretical calculations [2, 3]. Dasgupta et al. (DDSS) [2] assumed the jet function at to be a delta function and performed an angular evolution down to smaller radii. Kang et al. (KRV) [3] calculated the jet function at one scale using the “Soft Collinear Effective Field Theory” (SCET) [4] and evolved it to different scales by using “Dokshitzer-Gribov-Lipatov-Altarelli-Parisi" (DGLAP) [5] evolution. Also, there has been considerable interest in studying jet substructure in the last decade [6, 7], especially “Soft Drop”, a widely used jet substructure technique [8]. Even though there have been many theoretical studies on both small radius jets and “Soft Drop”, none have provided a common formalism to describe medium modified jets and “Soft Drop” together with leading hadrons.

In earlier work, leading hadrons with medium modified fragmentation functions have been studied in the “higher-twist” formalism [9, 10, 11, 12]. In this current effort, we extended the evolution of the single hadron fragmentation function [9, 10] to the single jet function and succeeded in describing the small radius medium modified jets. The jet function was first calculated at lower energy and mass scale and then boosted to a higher energy and evolved to a higher scale. Di-hadron fragmentation function evolution [11, 12] was reformulated to evolve the “di-prong function” in order to describe the soft dropped jet substructure. This calculation allows using a single formalism to describe medium modified hadrons, jets as well as jet substructure in a soft drop algorithm.

2 Single hadron fragmentation function and single jet function

Although the vacuum evolution of single hadron fragmentation function through DGLAP evolution is similar in most studies, medium evolution is different in different formalisms. The quark jet function is defined as the probability to find a jet with momentum fraction , energy and radius within the fragmentation products of a quark with energy , and the operator definition is given by,


Here, and sum over is assumed. The quark and gluon vacuum jet functions at lower energy ( GeV) and large radius were calculated from PYTHIA simulation in lieu of experimental data. Then they were boosted up to a higher energy ( GeV). Jets with high transverse momentum become collinear as they are boosted to higher energy. Since the virtuality, , and the momentum fraction, , are boost invariant, the boost shrinks the jet radius as the energy increases. For example, if we boost a jet function from  GeV to  GeV, the jet radius reduces by a factor of . Finally, the boosted jet functions are evolved in scale up to a higher scale ( GeV to  GeV) using the DGLAP equations,

Figure 1: Diagrammatic representation of DDSS approach, KRV approach, and our approach to jet function evolution

The evolution method that we use is different from the DDSS and KRV approaches. Fig. 1 illustrates these three evolution approaches in detail.

Our evolution approach was then extended to the medium modified evolution of jet functions. We started with the vacuum jet functions at lower energy calculated from PYTHIA simulations and boosted them up to a higher energy so that the jets become significantly narrower, and the formation time of the jets are larger so that they form outside the QGP medium. Then we evolved them up to a higher scale using the following medium modified evolution equation,




3 Di-hadron fragmentation function and Soft Drop function

The di-hadron fragmentation function and its evolution are calculated in an extended DGLAP formalism [11, 12]. We define a new non-perturbative function, namely "Di-Prong Function" (DPF), as the distribution of the combination of two prongs with momentum fractions and produced in the hadronization of the outgoing hard parton. The operator definition of DPF is given by,


Even though the operator definition of the DPF is similar to the definition of di-hadron fragmentation function, the evolution equation does not contain the terms responsible for a wider seperation of the two prongs than the radius of the jet, which contains the two prongs [11, 12]. However, the angular evolution changes the radius of the jet, and hence partonic splittings can cause a significant effect.

Similar to the previous procedure for the single jet function, we first calculated the DPFs at lower energy ( GeV) from PYTHIA simulations and then boosted them to a higher energy ( GeV) and finally evolved them to a higher scale ( GeV) by using the following evolution equation,


4 Results and discussion

The jet functions calculated from PYTHIA simulations at  GeV and  GeV as well as the evolved jet function at  GeV are shown in Fig. 2. The evolved jet functions have a reasonable agreement with a direct calculation from PYTHIA at  GeV. This validates our evolution approach for the jet function.

Figure 2: Calculated quark and gluon jet functions at  GeV and  GeV by using PYTHIA simulations and the evolved jet functions at  GeV

In addition to the vacuum evolution, we performed a medium modified evolution within the higher-twist formalism. The nuclear modification factor calculated for the extremely narrow jets was in good agreement with the single hadron results from CMS [13] as well as our direct calculation of single hadron by using the higher-twist formalism in Fig. 4. Even though this result was expected, this is the first time it has been explicitly demonstrated. When the radius of a jet becomes very small, most jets contain a single hadron, and hence the jet spectrum follows the single hadron spectrum.

Figure 3: Medium evolution of extremely narrow jets and leading hadrons in higher-twist formalism compared to data[13].
Figure 4: Calculated DPFs at  GeV and  GeV (with ) from PYTHIA and our evolved DPF at  GeV

Figure 4 shows the calculated DPFs from PYTHIA at  GeV and  GeV as well as our evolved DPFs at  GeV. The evolved DPF results have good agreement with the calculated DPF results in PYTHIA.

5 Summary and outlook

The vacuum evolution of the single jet function was in good agreement with the PYTHIA simulations. Medium evolution was performed for very narrow jets since the jet should be created outside the QGP medium for the evolution to work. We are working on the angular evolution approach to evolve the jet functions "in angle" to different jet radii inside a medium. Once the angular evolution is calculated, we can use the results from Relativistic Heavy Ion Collider (RHIC) to predict the results from Large Hadron Collider (LHC).

6 Acknowledgment

This work was supported in part by the U.S. Department of Energy (DOE) under grant number DE-SC0013460 and in part by the National Science Foundation (NSF) under grant number ACI-1550300 within the framework of the JETSCAPE collaboration. CS would like to thank Sara Tipton for the proof reading of the manuscript.


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