This Twiki page is to document results and ideas of Dmytro Oliinychenko during his summer studentship at CERN. The supervisor is Mateusz Ploskon. The topic of research is studies of jets in heavy ion collisions.

Introduction and main goal

Jets are a nice probe of Quark-Gluon Plasma (QGP). Jet quenching, i.e. attenuation or disappearance of the hadron spray born in hard interaction, was one of the first proposed signals of QGP [1]. Jet quenching was observed experimentally in heavy ion collision experiments at RHIC and LHC and is considered as one of the arguments for QGP existence. The phenomenon of jet quenching is interesting not only qualitatively as a QGP signal, but also quantitatively as a probe that allows to measure properties of the hot and dense quark-gluon medium. There is a number of models for parton energy loss in QGP, which allowed to extract transport coefficient \hat{q}, gluon density and initial energy density loss in QGP \epsilon_{0} from nuclear modification factor and di-hadron observables.

However, there is a hope to extract more information about QGP using jets. Fully recostructed jets seem to be one more promising probe, because they can give an additional information about medium. Taking jet energy profile and its content as observables we can directly (and hopefully model-independently) get the gluon medium-induced bremsstrahlung spectrum (suggested in [3]) and in-medium quark fragmentation functions. Maybe also some other quantities?

Fully reconstructed jets were considered as promising already in 2004[2]. To 2009 it was already carefully studied (e.g. [2,3,4]) at which jet observables one should look and which theoretical quantities one can extract from them. However, there is a reason that prevented people from accomplishing this study already at RHIC, even before LHC start. It is hard to distinguish jets that were born inside fireball from those born outside. The goal of my studies is to do this job, i.e. distinguish "true jets" from background. This can be reformulated as follows: such discriminator should be found that it will cut off almost all "false" jets born outside the fireball.



First step is to study p+p collisions without noise. Technically, procedure is the following: pp event at \sqrt_{S} = 2.76 TeV is modelled using Pythia8 with standard options and HardQCD:all = on. To make the simulation faster and to require less statistics, allowed hard process Q^2 in event is restricted to some bin. Afterwards data for all bins are put together and properly normalized. Final particles from simulation are delivered to FastJet. Jets are formed using anti-kt algorithm with R_jet = 0.4 and default recombination scheme (direct adding). There is a pseudorapidity cut applied to jets: |η|<3.5 Jet constituents are considered as an input for FastJet. Subjets are formed with kt algorithm with R_subjet = 0.2 and default recombination scheme. Results are as follows:

1. Jet pT and η More... Close From them we can conclude that i) number of jets with high transverse momentum is very small, ii) faster jets tend more to midrapidity. JetPt00.png JetEta00.png

2. Jet multiplicity distribution and number of constituents distribution More... Close We see that number of constituents per jet is growing with pT, at low pT bins distribution has clear maximum, while at high pT distribution becomes almost uniform. At the same time average number of jets per event is almost the same for all pT and never exceeds 6. ( I think this comes from lack of statistics.) JetMult00.png JetCon00.png

3. Let us now consider the behaviuor of subjets - their pT, η distributions, number of subjets per jet and fraction of pT carried by subjet - Z. More... Close In subjet pT distribution for jet pT > 20 GeV one can see maximums at the right wing. This means that at jet pT > 20 GeV there is often a situation when only one subjet carries most of the jet momentum. This conclusion is confirmed by Z (subjet pT fraction) distribution. Distribution of subjet pseudorapidity is similar to that of jets. No special effects of recombination scheme are visible. Number of subjets per jet tends to behave similarly for all jet pTs. Strangely, number of subjets distribution rapidly dies at number of subjets = 8 for all pT ( again not enough statistics). Average number of subjets per jet is 4-5. With increasing pT number of subjets per jet remains the same in average, but number of constituents per jet grows. Thus, at higher jet pT subjets have a bigger number of constituents. SubJetPt01.png SubJetEta00.png SubJetMult00.png


Plots for varying subjet radius: More... Close SubJetPt_rsub_var1.pngSubJetMult_rsub_var1.png


With the help of Mateusz I have realised that previous plots were by pTHat, not by pT. In addition they weren't properly normalized. Again, pp collision was simulated via Pythia8, particles with |η| < 3.5 - Q R_jet (Q = 1, R_jet = 0.4) were taken for jet reconstruction. Particles from jets were taken as a new input for FastJet reconstruction with R_subjet = 0.1, 0.2, 0.3. Here are 5 histograms for subjets: subjet Z (pT_subjet/pT_jet), leading subjet Z, subleading subjet Z, subjet normalized projection to jet momentum (p_subjet, p_jet)/p_jet^2, subjet pT: More... Close SubjetPt_pt_bins.pngSubjetPr_pt_bins01.pngSubjetZ_pt_bins01.pngSubjetZL1_pt_bins01.pngSubjetZL2_pt_bins01.png

There is also one interesting quantity refered to as jet girth[5] - linear radial moment ∑ pT(particle)/pT(jet) * r(particle), where r(particle) is the Euclidean distance between particle and jet axis in η - φ plane. I have divided girth by jet radius to obtain value which always lies between 0 and 1. Here is what I got: More... Close JetGirth_pt_bins.png


pp collision was simulated via Pythia8, particles with |η| < 2.8 - Q R_jet (Q = 1, R_jet = 0.4) were taken for jet reconstruction. Particles from jets were taken as a new input for FastJet reconstruction with R_subjet = 0.2. Leading subjet pT fraction distribution is given, bumps in the central region should correspond to gluon jets, right jumps - to quark jets.SubJetZL1_lin00_0.2.png


Using Pythia6 I have generated gluon- and d-jets with narrow Pt spectrum around 45 GeV. To produce such jets I used Py1ent and generated two gluons or d and anti-d quarks with given Pt, that fly into opposite directions.

pythia6->Py1ent(-1, 21, 50., TMath::Pi()/2., 0.);

pythia6->Py1ent( 2, 21, 50., TMath::Pi()/2., TMath::Pi());

R_subjet = 0.2, algorithm for subjet extraction is kt.

More... Close ZL1.png

Analoguos plots for R_subjet = 0.1:

More... Close ZL1_antikt_Rsub0.1.pngZL1_Rsub0.1_qq_gg.png

Kt subjet exctraction algorithm is on the left panel, anti-kt subjet extraction algorithm is on the right panel.

This is similar qualitatively to http://jets.physics.harvard.edu/qvg/ results, but it differs quantitatively. My assumption is that my quark jets are purely d jets, while their quark jets are a mixture. This leads me to an attempt to generate different flavour jets. Kt subjet exctraction algorithm is on the left panel, anti-kt subjet extraction algorithm is on the right panel.

More... Close ZL1_Subjet_kt_flavours.pngZL1_Subjet_antikt_flavours.png


I have gone in for Toy Model to consider background effects. I generate gluon jets and d-jets with Pt around 45 GeV (40 - 50 GeV, rather sharp peak) using Pythia6. To particles generated in this way I add "background" with Boltzmann spectrum, T = 700 MeV, average background multiplicity 1917 particles per event. Eta and phi are uniformly distributed. After that I consider three types of jets: I) hard jets (true jets), II) pure background jets (fake jets) - jets that consist purely from background particles and III) full jets merged with hard (half constituents of the jet match with some hard jet). For these jets Z (subjet Pt fraction) distribution was built. Additionally, constituents lower Pt cut was considered.

More... Close ZL1Hard.pngZL1Fake.pngZL1F.png


Same Toy Model and conditions as earlier but statistics 10 times increased (thank to pdsf.nersc.gov grid) and more quantities are plotted: Pt (hard jets, full jets, full merged with hard), leading subjet Pt fraction Z (hard jets, fake jets, full merged with hard), Z with background exctraction (fake jets and full merged with hard), background density rho.




I tried to generate jets with Pythia8 and check for bumps at Z=1. No bumps. Huge crosses must me due to sides of distribution when summing by PtHat.


Pt versus leading subjet Pt in Pythia 6: PtPtL1subjet_ddbar.png


ToyModel. Pt (hard jets, full jets, full merged with hard), leading subjet Pt fraction Z (hard jets, fake jets, full merged with hard), Z with background exctraction (fake jets and full merged with hard), background density rho. Pt of jets is 1) 15-20GeV and 2) 70-80 GeV.

20 GeV: More... Close PtHard20GeV.pngPtFull20GeV.pngPtF20GeV.pngZL1Hard20GeV.pngZL1Fake20GeV.pngZL1F20GeV.pngZL1FakeB20GeV.pngZL1FB20GeV.pngrho20GeV.png

80 GeV: More... Close PtHard80GeV.pngPtFull80GeV.pngPtF80GeV.pngZL1Hard80GeV.pngZL1Fake80GeV.pngZL1F80GeV.pngZL1FakeB80GeV.pngZL1FB80GeV.pngrho80GeV.png


I have looked into jet contents. How many mesons and baryons are there in gluon/light quark jets? Here is the plot to answer it.



With the help of Mateusz and Xiaoming I tried to analyze data, keeping in mind my previous MC studies. It took me a week to get some minimal understanding how data analysis is done at ALICE.


[1] J. D. Bjorken, FERMILAB-PUB-82-059-THY (1982)

[2]C. A. Salgado and U. A. Wiedemann, Phys. Rev. Lett. 93 , 042301 (2004)

[3]I. Vitev, S. Wicks and B. W. Zhang, arXiv:0810.2807 [hep-ph]

[4]F. Arleo, arXiv:0810.1193v1 [hep-ph], 2008

[5]J. Gallicchio and M. D. Schwartz; Quark and Gluon Tagging at the LHC; 2011

Topic attachments
I Attachment History Action Size Date Who Comment
Postscriptps JetCon.ps r1 manage 38.4 K 2013-07-06 - 19:04 OliinychenkoDmytro  
Postscriptps JetEta.ps r1 manage 49.4 K 2013-07-06 - 19:04 OliinychenkoDmytro  
Postscriptps JetMult.ps r1 manage 23.9 K 2013-07-06 - 19:04 OliinychenkoDmytro  
JPEGjpg JetPt.jpg r2 r1 manage 56.9 K 2013-07-06 - 19:29 OliinychenkoDmytro  
Postscriptps JetPt.ps r1 manage 29.8 K 2013-07-06 - 19:03 OliinychenkoDmytro  
Postscriptps SubJetEta.ps r1 manage 51.2 K 2013-07-06 - 19:05 OliinychenkoDmytro  
Postscriptps SubJetMult.ps r1 manage 23.0 K 2013-07-06 - 19:05 OliinychenkoDmytro  
PNGpng SubJetMult_rsub_var1.png r1 manage 71.2 K 2013-07-08 - 17:46 OliinychenkoDmytro  
Postscriptps SubJetPt.ps r1 manage 78.9 K 2013-07-06 - 19:04 OliinychenkoDmytro  
PNGpng SubJetPt_rsub_var1.png r1 manage 112.1 K 2013-07-08 - 17:46 OliinychenkoDmytro  
Postscriptps SubJetZ.ps r1 manage 52.3 K 2013-07-06 - 19:05 OliinychenkoDmytro  
PNGpng SubJetZL1_lin00_0.2.png r1 manage 63.2 K 2013-07-18 - 15:33 OliinychenkoDmytro  
PNGpng ZL1F.png r1 manage 87.3 K 2013-07-29 - 17:09 OliinychenkoDmytro  
PNGpng ZL1Fake.png r1 manage 64.0 K 2013-07-29 - 17:09 OliinychenkoDmytro  
PNGpng ZL1Hard.png r1 manage 98.9 K 2013-07-29 - 17:09 OliinychenkoDmytro  
PNGpng ZL1_Rsub0.1_qq_gg.png r1 manage 28.7 K 2013-07-24 - 13:58 OliinychenkoDmytro  
Edit | Attach | Watch | Print version | History: r19 < r18 < r17 < r16 < r15 | Backlinks | Raw View | WYSIWYG | More topic actions
Topic revision: r19 - 2013-08-23 - OliinychenkoDmytro
    • Cern Search Icon Cern Search
    • TWiki Search Icon TWiki Search
    • Google Search Icon Google Search

    Sandbox All webs login

This site is powered by the TWiki collaboration platform Powered by PerlCopyright & 2008-2020 by the contributing authors. All material on this collaboration platform is the property of the contributing authors.
or Ideas, requests, problems regarding TWiki? use Discourse or Send feedback