250 likes | 501 Views
Jettiness and Reaction Plane determination at LHC. University of Tsukuba Dousatsu Sakata. Workshop at Wuhan 2008/12/4. Outline. HYDJET++(HYDJET) Reaction Plane resolution Jet Reconstruction to be done. HYDJET++(HYDJET).
E N D
Jettiness and Reaction Plane determination at LHC University of Tsukuba Dousatsu Sakata Workshop at Wuhan 2008/12/4
Outline • HYDJET++(HYDJET) • Reaction Plane resolution • Jet Reconstruction • to be done... Work shop at Wuhan 2008/12/4
HYDJET++(HYDJET) HYDJET++ is event generator for high energy heavy ion collision. HYDJET is previous one. http://lokhtin.web.cern.ch/lokhtin/hydjet++/ Hard part: PYTHIA or PYQUEN (modifies PYTHIA6.4 for jet quench effect) Soft part: HYDJET++ (represents the ”thermal” hadronic state FASTMC) Set Parameter (LHC default) A :207(Pb) CME : 5500GeV Set :hydro + pyquen hard(&quenching) soft pthatmin: 7.GeV/c nucl. shadowing : on T0 : 0.8GeV Tf :130MeV t0 :0.1fm/c Tc :170MeV # of flavor:0 ttf :10fm/c max ly flow:4. max ty flow:1.1 unsure very peripheral |h|>3 Work shop at Wuhan 2008/12/4
Particle emission of HYDJET++ These figure are output of sample code of HYDJET++ to compare RHIC data and HYDJET++ Data RHIC AuAu 200GeV left: pT distribution right: rapidity distribution point : data from STAR point: data from PHOBOS black line: HYDJET++ black line: HYDJET++ Work shop at Wuhan 2008/12/4
Reaction Plane resolution HYDJET [hydro+pythia] Py Pz pT weight Reaction plane Px pT cut for RP <2GeV/c ITS+TPC A: h>0 B: h<0 -0.9<h<0.9 -5.1<h<-1. 1.7<h< 3.4 FMD ITS+TPC(pT wight) have a good Reaction Plane resolution. -5.0<h<-4.5 2.9<h< 3.3 V0 -3.2<h<-2.9 5.0<h< 6.0 T0 Work shop at Wuhan 2008/12/4
Compare mesured-v2 and true-v2 ptcut for RP : 3GeV/c ratio= mesuredv2 /truev2 [truev2: v2 with respect to true RP] No pT dependence centrality 20-40% we have to minimize non-flow contribution |h|<1 Work shop at Wuhan 2008/12/4
Methodfor the Reconstruction seed particle Purpose estimate R & parton pT dependence of eTjet/eTpart & resolution Data set PYTHIA - hardQCD:all on softQCD:all off - primary parton Pt :20-200GeV/c(20GeV/c step) HYDJET++ - LHC default • Find seed particle find highest eT track (only a particle par event) • Sum up eT in Cone (eTjet) • Compare parton eT (eTpart) calculate eTjet/eTpart Cone Radius Primary parton Figure PYTHIA (hard only) pThat 100-120GeV/c Work shop at Wuhan 2008/12/4
Jet Reconstruction (PYTHIA) distribution of eTjet/eTpart (R=0.4) R dependence of eTjet/eTpart&resolution pThat dependence of eTjet/eTpart&resolution(R=0.4) R=0.4 color : differential pT of primary parton [20-200GeV/c(20GeV/c step)] eTjet/eTpart & resolution are better when R is larger, and the larger pT of primary parton also makes resolution better. at pthat=40 eTjet/eTpart~0.9 resolution[RMS/Mean] < 20% Work shop at Wuhan 2008/12/4
Soft and Hard particle in HYDJET++ In eT>3GeV , particles come from hard more than soft. Almost particles comes from hard collision in eT>8GeV Work shop at Wuhan 2008/12/4
Jet Reconstruction (HYDJET++) we like to get information of jet. So we estimate R dependence of hard eT ratio. setting HYDJET++ [HYDJET++ & PYQUEN] eT cut [seedeTcut=0] Centrality seed eT cut [eTcut=0] Yaxis: SeThard/SeTall Xaxis: cone radius R Smaller R is better to get information of Jets. Work shop at Wuhan 2008/12/4
Combine PYTHIA & HYDJET++ To analyze jets with information of primary parton. Mix hard event of PYTHIA and soft event of HYDJET++ HYDJET++ (soft) PYTHIA (hard) HYDJET++&PYTHIA h-f distribution Radius direction: h [-10<h<10] Anguler direction: f [-p<f<p] Hight direction: pT[GeV/c] Work shop at Wuhan 2008/12/4
To be done... • Find & reconstructjets in HYDJET++ with di-jets. • Estimate # of jets at # of Ncol. And mix HYDJET++. • Fix cut parameter to reconstruct & find jets • Include Detector effects. • Estimate trigger rate ,efficiency, ・・・ Work shop at Wuhan 2008/12/4
Conclusion • EstimateRP resolution & compare v2 and true-v2 • Testof simple reconstruction - PYTHIA eTjet/eTpart~0.9 resolution<20% at pThat=40GeV/c larger R is better to reconstruct jets - HYDJET++ smaller R is better to find jets and seed particle in eT>8GeV/c is almost come from hard collision in eT > 3GeV/c particles in cone come from hard more than soft And we have to do much more work! Work shop at Wuhan 2008/12/4
Backup Work shop at Wuhan 2008/12/4
ALICE detector ITS -0.9<η<0.9 -5.1<η<-1.7 1.7<η< 3.4 FMD V0 -5.0<η<-4.5 2.9<η< 3.3 T0 -3.2<η<-2.9 5.0<η< 6.0 Work shop at Wuhan 2008/12/4
Pt & Et of primaly parton & seed particle violet : parton pT or eT other color : trigger particle pT or eT Work shop at Wuhan 2008/12/4
Particle distribution of jet axis HYDJET++ hard ratio vs R PYTHIA particle distribution vs dφ dφ= φ-φparton Work shop at Wuhan 2008/12/4
Detail of jet reconstruct R=0.4 R=0.4 Work shop at Wuhan 2008/12/4
Njet correlation Work shop at Wuhan 2008/12/4
Reaction Plane Introduction Z Py Pz Y X Px Reaction plane Work shop at Wuhan 2008/12/4 3
Method y y φ Ψ x x Work shop at Wuhan 2008/12/4 6