Next Article in Journal
Deflection Angle of Photons through Dark Matter by Black Holes and Wormholes Using Gauss–Bonnet Theorem
Next Article in Special Issue
Quarkonium Production in the QGP
Previous Article in Journal
Multiverse—Too Much or Not Enough?
Previous Article in Special Issue
Viscous Hydrodynamic Description of the Pseudorapidity Density and Energy Density Estimation for Pb+Pb and Xe+Xe Collisions at the LHC
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Communication

Study of Jet Shape Observables in Au+Au Collisions at s N N = 200 GeV with JEWEL

by
Veronika Agafonova
Nuclear Physics Institute of the Czech Academy of Sciences, CZ-25068 Řež, Czech Republic
Universe 2019, 5(5), 114; https://doi.org/10.3390/universe5050114
Submission received: 31 March 2019 / Revised: 7 May 2019 / Accepted: 8 May 2019 / Published: 11 May 2019
(This article belongs to the Special Issue The Zimányi School and Analytic Hydrodynamics in High Energy Physics)

Abstract

:
Nuclear–nuclear collisions at energies attainable at the large accelerators RHIC and the LHC are an ideal environment to study nuclear matter under extreme conditions of high temperature and energy density. One of the most important probes of such nuclear matter is the study of production of jets. In this article, several jet shape observables in Au+Au collisions at the center of mass energy per nucleon–nucleon pair of s N N = 200 GeV simulated in the Monte Carlo generator JEWEL are presented. Jets were reconstructed using the anti- k T algorithm and their shapes were studied as a function of the jet-resolution parameter R, transverse momentum p T and collision centrality.

1. Introduction

The study of production of jets is one of the most important probes of nuclear matter under extreme conditions of high temperature and energy density. The jet is a collimated spray of hadrons originating from fragmentation of a hard parton created in the initial stage of the nucleus–nucleus collision and can be used for tomography of the nuclear matter (Figure 1). As jets mostly conserve the energy and the direction of the originating parton, they are measured in particle detectors and studied to determine the properties of the original quarks.
To probe the complimentary aspects of the jet fragmentation and constrain theoretical description of jet–medium interactions, different observables related to shapes of jets are studied at the CERN Large Hadron Collider (LHC) [2,3,4]. It is important to perform similar measurements at lower collision energies at the Relativistic Heavy Ion Collider (RHIC) taking advantage of new high statistics data. The jet substructure observables are the perfect tool to understand what is happening when the particles interact with Quark–Gluon Plasma (QGP) medium. This article is focused on two jet shape observables: the girth and momentum dispersion that will be described in more detail below.

2. Jet Shape Observables

The first radial moment (alternatively angularity or girth), g, probes the radial distribution of radiation inside a jet. It is defined as
g = i j e t p T i p T , jet | Δ R i , jet | .
Here p T i represents the momentum of the ith jet constituent and Δ R i , jet is the distance in η × ϕ plane between the constituent i and the jet axis [5], where η is the pseudorapidity and ϕ is the azimuthal angle. This type of shape is sensitive to the radial energy profile or broadening of the jet. In the collinear limit for the polar angle θ 0 the radial moment becomes equivalent to jet broadening.
The next observable is momentum dispersion, p T D . It measures the second moment of the constituent p T distribution in the jet and is connected to hardness or softness of the jet fragmentation. It means that in the case of a large number of constituents and softer momentum the p T D tends to 0, while in the opposite situation the p T D will be close to 1. Its definition is given by the equation:
p T D = i j e t p T , i 2 i j e t p T , i .
These two jet shape observables are infrared and collinear (IRC) safe. It means that if one modifies an event by a collinear splitting or the addition of a soft emission, the set of hard jets that are found in the event should remain unchanged [6].

3. The Anti- k T Algorithm

Jets are commonly reconstructed using the anti- k T clustering algorithm [7]. The anti- k T algorithm is a sequential-clustering algorithm. The algorithm is based on successive pair-wise recombination of particles and it works as follows. Firstly, the distance, d i j , between particles i and j is found as
d i j = min ( k t i 2 , k t j 2 ) Δ i j 2 R 2 ,
where Δ i j 2 = ( y i y j ) 2 + ( ϕ i ϕ j ) 2 and k t i or k t j , y i , ϕ i and R stand for the transverse momenta, rapidity, azimuth, and radius parameter of particle i respectively. Secondly, the algorithm calculates the distance, d i B between the entity i and the beam B as
d i B = k t i 2 .
The next step of the anti- k T jet algorithm is to find the minimum distance, d m i n , between the distances d i j and d i B . In case the smallest distance is d i j , the algorithm performs a recombination of the entities. In other situation i is called to be a jet and is subsequently removed from the list. All these steps are repeated until no particles are left.

4. JEWEL

Jet Evolution With Energy Loss (JEWEL) is a Monte Carlo event generator that describes the Quantum Chromodynamics (QCD) evolution of jets in vacuum and in a medium in a perturbative approach [8,9,10]. In this section, the simulation in JEWEL will be described. For this research 50 million events were simulated for the interaction in vacuum and 20 million events for the interaction in medium. The simulation was made for 0–10% central and 60–80% peripheral Au+Au collisions with additional “recoils on/off” option for interaction with medium. “Recoils on” option in JEWEL keeps the thermal partons recoiling against interactions with the jet in the event and let them hadronize together with the jet, while the “recoils off” option ignore the medium response [11]. All events were required to have the center-of-mass (CMS) energy s NN = 200 GeV.
Table 1 contains the parameters used for the vacuum model. Additional parameters for the simulation with the medium can be found in Table 2.
A resolution parameter, R, quantifies the size of the jet. For this study values of the resolution parameter were chosen to be R = 0.2 and R = 0.4, respectively. The charged particles were simulated in pseudorapidity η c e n t = 2.5 and full azimuth. Jets were reconstructed with the anti- k T algorithm included in FastJet software package [12].

5. Results

In this section, only the JEWEL results for central Au+Au collisions at s N N = 200 GeV are presented as they are more appealing from the physical point of view. The jet shape observables are calculated for different values of the resolution parameter R and charged jet p T separately for vacuum and medium with “recoils on/off” option. The distributions will be further compared to the results and the JEWEL simulation from the ALICE collaboration [13].
Figure 2 shows the measured jet shape distributions in 0–10% central Pb–Pb collisions at s NN = 2.76 TeV for anti- k T charged jets at ALICE compared to JEWEL simulation with and without recoils [13]. As the resolution parameter is small, R = 0.2, the effects of medium recoils are also small. It means that the measurement is constrained by purely radiative aspects of the JEWEL shower modification. A good agreement between the data and the model, especially in momentum dispersion, can be observed.
Figure 3 and Figure 4 compare the distributions of angularity for vacuum and medium “recoils on/off” central Au+Au collisions in two different p T ranges 10 < p T < 20 GeV/c and 20 < p T < 30 GeV/c, respectively. As it can be seen, the first radial moment has the same behavior for R = 0.2 as the results from the ALICE experiment (Figure 2). Nevertheless, peaks for the medium “recoils on” and medium “recoils off” simulation of angularity with R = 0.4 are shifted to the right and left, respectively. Distributions for medium “recoils on” collisions with R = 0.4 have a longer tail than others. Also, the spike for g = 0.01 in the case of jets with 10 < p T < 20 GeV/c can be observed for both resolution parameters. That signals the presence of jets with only one constituent. To probe this, the dependence of the number of constituents on the angularity is shown in Figure 5. It can be clearly seen that there is a larger number of particles for R = 0.4 than for R = 0.2 jets.
Figure 6 and Figure 7 compare the results for the momentum dispersion for jets with 10 < p T < 20 GeV/c and 20 < p T < 30 GeV/c, respectively. As for previous observable, there is a better agreement between the models in 20 < p T < 30 GeV/c p T range. However, in contradiction to the ALICE results, the obtained distributions for the momentum dispersion start form p T D = 0 (for R = 0.4 in central and peripheral collisions) and p T D = 0.1 (for R = 0.2 in central collisions) instead of p T D = 0.3. That can be a consequence of the use of different centrality ranges. Also, a shift of the distribution to lower values for the central medium “recoils on” setting for R = 0.4 and 10 < p T < 20 GeV/c can be observed.

6. Conclusions

In this article, the results of study of two jet shape observables, girth and momentum dispersion, in central Au+Au collisions at CMS energy of 200 GeV per nucleon–nucleon pair with the JEWEL Monte Carlo generator were presented. The jet shapes were calculated using the anti-k T jet finding algorithm implemented in the FastJet software package. The chosen observables were studied as a function of the transverse momentum, jet-resolution parameter, and collision centrality. All the obtained results have the same behavior as the results from the ALICE collaboration [13]. In this study, it was shown that the spike in the girth results for g = 0.01 for jets with p T 10–20 GeV/c for both values of resolution parameter, R = 0.2 and R = 0.4, (Figure 3) is due to presence of jets with only one constituent (Figure 5). As for the momentum dispersion, the obtained distributions (Figure 6 and Figure 7) are wider in comparison to the ALICE results.
One of the goals of future work is to perform the background subtraction similarly to the ALICE experiment. It is expected that after the background subtraction the points for medium “recoils on/off” and vacuum models will be closer to each other.

Funding

This research was funded by INTER-EXCELLENCE INTER-TRANSFER of Ministry of Education, Youth and Sports of the Czech Republic grant number LTT18002.

Conflicts of Interest

The authors declare no conflict of interest.

Abbreviations

The following abbreviations are used in this manuscript:
LHCLarge Hadron Collider
RHICRelativistic Heavy-Ion Collider
CERNConseil européen pour la recherche nucléaire
QGPQuark Gluon Plasma
QCDQuantum Chromodynamics
CMScenter-of-mass
JEWELJet Evolution With Energy Loss
IRC safeInfrared and Collinear safe
ALICEA Large Ion Collider Experiment

References

  1. Agafonova, V. Study of Jet Shape Observables at RHIC. Available online: https://physics.fjfi.cvut.cz/publications/ejcf/bp_ejcf_17_agafonova.pdf (accessed on 11 May 2019).
  2. Aad, G.; Abbott, B.; Abdallah, J.; Abdelalim, A.A.; Abdesselam, A.; Abdinov, O.; Abi, B.; Abolins, M.; Abramowicz, H.; Abreu, H.; et al. Observation of a Centrality-Dependent Dijet Asymmetry in Lead-Lead Collisions at s N N = 2.77 TeV with the ATLAS Detector at the LHC. Phys. Rev. Lett. 2010, 105, 252303. [Google Scholar] [CrossRef] [PubMed]
  3. Chatrchyan, S. et al. [CMS Collaboration] Observation and studies of jet quenching in PbPb collisions at nucleon-nucleon center-of-mass energy = 2.76 TeV. Phys. Rev. C 2011, 84, 024906. [Google Scholar] [CrossRef]
  4. Vitev, I.; Wicks, S.; Zhang, B.W. A Theory of jet shapes and cross sections: From hadrons to nuclei. J. High Energy Phys. 2008, 11, 93. [Google Scholar] [CrossRef]
  5. Cunqueiro, L. et al. [ALICE Collaboration] Jet shapes in pp and Pb–Pb collisions at ALICE. Nucl. Phys. A 2016, 956, 593–596. [Google Scholar] [CrossRef] [Green Version]
  6. Salam, G.P. Towards Jetography. Eur. Phys. J. C 2010, l67, 637–686. [Google Scholar] [CrossRef]
  7. Cacciari, M.; Salam, G.P.; Soyez, G. The Anti-k(t) jet clustering algorithm. J. High Energy Phys. 2008, 4, 63. [Google Scholar] [CrossRef]
  8. Zapp, K.C. JEWEL 2.0.0: Directions for use. Eur. Phys. J. 2014, 74, 2762. [Google Scholar] [CrossRef]
  9. Zapp, K.C. Geometrical aspects of jet quenching in JEWEL. Phys. Lett. B 2014, 735, 157–163. [Google Scholar] [CrossRef] [Green Version]
  10. Zapp, K.C.; Ingelman, G.; Rathsman, J.; Stachel, J.; Wiedemann, U.A. A Monte Carlo Model for ‘Jet Quenching’. Eur. Phys. J. C 2009, 60, 617–632. [Google Scholar] [CrossRef]
  11. Zapp, K.C.; Kunnawalkam Elayavalli, R. Medium response in JEWEL and its impact on jet shape observables in heavy ion collisions. J. High Energy Phys. 2017, 7, 141. [Google Scholar] [CrossRef]
  12. Cacciari, M.; Salam, G.P.; Soyez, G. FastJet user manual. Eur. Phys. J. C 2012, 72, 1896. [Google Scholar] [CrossRef]
  13. Acharya, S.; Adamová, D.; Adler, A.; Adolfsson, J.; Aggarwal, M.M.; Rinella, G.A.; Agnello, M.; Agrawal, N.; Ahammed, Z.; Ahn, S.U.; et al. Medium modification of the shape of small-radius jets in central Pb-Pb collisions at s NN = 2.76 TeV. J. High Energy Phys. 2018, 10, 139. [Google Scholar] [CrossRef]
Figure 1. A schematic view of jet created in a heavy-ion collision [1].
Figure 1. A schematic view of jet created in a heavy-ion collision [1].
Universe 05 00114 g001
Figure 2. Jet shape distributions g (left) and p T D (right) in 0–10% central Pb–Pb collisions at s NN = 2.76 TeV for R = 0.2 in range of jet p T , jet ch of 40–60 GeV/c compared to JEWEL with and without recoils with different subtraction methods. The colored boxes represent the experimental uncertainty on the jet shapes [13].
Figure 2. Jet shape distributions g (left) and p T D (right) in 0–10% central Pb–Pb collisions at s NN = 2.76 TeV for R = 0.2 in range of jet p T , jet ch of 40–60 GeV/c compared to JEWEL with and without recoils with different subtraction methods. The colored boxes represent the experimental uncertainty on the jet shapes [13].
Universe 05 00114 g002
Figure 3. Girth for jets with p T of 10–20 GeV/c and R = 0.2 (left) and R = 0.4 (right) in central Au+Au collisions at s NN = 200 GeV.
Figure 3. Girth for jets with p T of 10–20 GeV/c and R = 0.2 (left) and R = 0.4 (right) in central Au+Au collisions at s NN = 200 GeV.
Universe 05 00114 g003
Figure 4. Girth for jets with p T of 20–30 GeV/c and R = 0.2 (left) and R = 0.4 (right) in central Au+Au collisions at s NN = 200 GeV.
Figure 4. Girth for jets with p T of 20–30 GeV/c and R = 0.2 (left) and R = 0.4 (right) in central Au+Au collisions at s NN = 200 GeV.
Universe 05 00114 g004
Figure 5. 2D statistics for jets with p T of 10–20 GeV/c and R = 0.2 (left) and R = 0.4 (right) in central “recoils on” Au+Au collisions at s NN = 200 GeV simulated with medium.
Figure 5. 2D statistics for jets with p T of 10–20 GeV/c and R = 0.2 (left) and R = 0.4 (right) in central “recoils on” Au+Au collisions at s NN = 200 GeV simulated with medium.
Universe 05 00114 g005
Figure 6. Momentum dispersion for jets with p T of 10–20 GeV/c and R = 0.2 (left) and R = 0.4 (right) in central Au+Au collisions at s NN = 200 GeV.
Figure 6. Momentum dispersion for jets with p T of 10–20 GeV/c and R = 0.2 (left) and R = 0.4 (right) in central Au+Au collisions at s NN = 200 GeV.
Universe 05 00114 g006
Figure 7. Momentum dispersion for jets with p T of 20–30 GeV/c and R = 0.2 (left) and R = 0.4 (right) in central Au+Au collisions at s NN = 200 GeV.
Figure 7. Momentum dispersion for jets with p T of 20–30 GeV/c and R = 0.2 (left) and R = 0.4 (right) in central Au+Au collisions at s NN = 200 GeV.
Universe 05 00114 g007
Table 1. Parameters of the JEWEL vacuum simulation for central and peripheral collisions [8].
Table 1. Parameters of the JEWEL vacuum simulation for central and peripheral collisions [8].
Name of ParameterName in JEWELValue
Parton Distribution Function setPDFSET10,100
Number of eventsNEVENT100,000
Mass number of Au nucleusMASS197
The CMS energy of the colliding systemSQRTS, [GeV]200
Minimum p T in matrix elementPTMIN, [GeV]3
Maximum p T in matrix elementPTMAX, [GeV]−1
The rapidity rangeETAMAX2.5
Table 2. Parameters of the JEWEL simulation with medium for central and peripheral “recoils on/off” collisions [8].
Table 2. Parameters of the JEWEL simulation with medium for central and peripheral “recoils on/off” collisions [8].
Name of ParameterName in JEWELValue
The initial (mean) temperatureTI, [GeV]0.28
The initial time τ i TAUI, [fm]0.6
An integer mass number of colliding nucleiA197
The lower end of centrality rangeCENTRMIN, [%]060
The upper end of centrality rangeCENTRMAX, [%]1080
The switch of keeping recoilsKEEPRECOLISTF
The nucleus–nucleus cross-sectionSIGMANN, [fm 2 ]4.2

Share and Cite

MDPI and ACS Style

Agafonova, V. Study of Jet Shape Observables in Au+Au Collisions at s N N = 200 GeV with JEWEL. Universe 2019, 5, 114. https://doi.org/10.3390/universe5050114

AMA Style

Agafonova V. Study of Jet Shape Observables in Au+Au Collisions at s N N = 200 GeV with JEWEL. Universe. 2019; 5(5):114. https://doi.org/10.3390/universe5050114

Chicago/Turabian Style

Agafonova, Veronika. 2019. "Study of Jet Shape Observables in Au+Au Collisions at s N N = 200 GeV with JEWEL" Universe 5, no. 5: 114. https://doi.org/10.3390/universe5050114

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop