FLUID + JETYasuki Tachibana
QCD相転移やQGP生成のモデル化による重イオン衝突の時空発展の理解に向けた理論・実験共同研究会
September 24th, 2021

Jet Evolution in JETSCAPE

Y. Tachibana, QCD相転移やQGP生成のモデル化による重イオン衝突の時空発展の理解に向けた理論・実験共同研究会, September 24th, 2021

JETSCAPE framework
MC event generator package for heavy ion collisions-

General, modular and extensible-

Communication between modules

Cooper-FryeSampling

Viscous Fluid Dynamics for
MediumHadronic Cascade

Hard & Semi-hardHadoronization

Hard Particle
Production

Multi-stage
Jet Shower Evolution
Energy-momentum Deposition

Initial Soft
Density
distribution

Initial geometry of
Nucleus-Nucleus collision

JETSCAPE Event Generator

Y. Tachibana, QCD相転移やQGP生成のモデル化による重イオン衝突の時空発展の理解に向けた理論・実験共同研究会, September 24th, 2021

3

Multi-scale jet evolution

Majumder, Putschke (16), JETSCAPE (17)In-vacuum-

In-vacuum: Virtuality ordered splitting

Small-Q
Large-Q

Q 2: virtuality (off-shellness)
Y. Tachibana, QCD相転移やQGP生成のモデル化による重イオン衝突の時空発展の理解に向けた理論・実験共同研究会, September 24th, 2021

4

Multi-scale jet evolution

Majumder, Putschke (16), JETSCAPE (17)In-medium-In-vacuum: Virtuality ordered splitting

Small-Q
Large-Q

Q 2: virtuality (off-shellness)
Y. Tachibana, QCD相転移やQGP生成のモデル化による重イオン衝突の時空発展の理解に向けた理論・実験共同研究会, September 24th, 2021

4

Multi-scale jet evolution

Majumder, Putschke (16), JETSCAPE (17)In-medium
Small-Q-In-vacuum: Virtuality ordered splitting-

Large-Q: Medium effect on the top of in-vacuum splitting(Medium modified splitting)Large-Q

Q 2: virtuality (off-shellness)
Y. Tachibana, QCD相転移やQGP生成のモデル化による重イオン衝突の時空発展の理解に向けた理論・実験共同研究会, September 24th, 2021

4

Multi-scale jet evolution

Majumder, Putschke (16), JETSCAPE (17)In-medium
Small-Q-In-vacuum: Virtuality ordered splitting-

Large-Q: Medium effect on the top of in-vacuum splitting(Medium modified splitting)Large-Q-

Small-Q: Splitting driven almost purely by medium effect(Medium induced splitting)Q 2: virtuality (off-shellness)
Y. Tachibana, QCD相転移やQGP生成のモデル化による重イオン衝突の時空発展の理解に向けた理論・実験共同研究会, September 24th, 2021

4

Multi-scale jet evolution

Majumder, Putschke (16), JETSCAPE (17)In-medium
Small-Q-In-vacuum: Virtuality ordered splitting-

Large-Q: Medium effect on the top of in-vacuum splitting(Medium modified splitting)Large-Q-

Small-Q: Splitting driven almost purely by medium effect(Medium induced splitting)Cannot be described by a single model→ Combination of multiple models
Q 2: virtuality (off-shellness)
Y. Tachibana, QCD相転移やQGP生成のモデル化による重イオン衝突の時空発展の理解に向けた理論・実験共同研究会, September 24th, 2021

4

Multi-stage description for jet evolution

JETSCAPE (17)Small-Q
Large-Q

Virtuality separation scale: Qsw
Large-Q ( > Qsw)

MATTER
Majumder(13), Kordell, Majumder(17),Cao, Majumder(17)Virtuality ordered splitting
Higher Twist
formalism

Small-Q ( < Qsw)

LBT

MARTINIWang, Zhu(13), Luo, et al.(15,18)
Cao, et al.(16,17), He, et al.(18)

Schenke, Gale, Jeon(09),Park, Jeon, Gale(17,18)

On-shell parton transport
Higher Twist
formalism

AMY
formalism

AdS/CFT
Chesler, Rajagopal(14, 15)
Pablos, et al.(15,16,17)

Diffusion into medium

𝒩 = 4 superYang-Mills

Y. Tachibana, QCD相転移やQGP生成のモデル化による重イオン衝突の時空発展の理解に向けた理論・実験共同研究会, September 24th, 2021

5

Multi-stage description for jet evolution

JETSCAPE (17)Small-Q
Large-Q

Virtuality
separation scale: Qsw
Switching between modules for parton
by parton
Large-Q ( > Qsw)

MATTER
Majumder(13), Kordell, Majumder(17),Cao, Majumder(17)Virtuality ordered splitting
Higher Twist
formalism

Small-Q ( < Qsw)

LBT

MARTINIWang, Zhu(13), Luo, et al.(15,18)
Cao, et al.(16,17), He, et al.(18)

Schenke, Gale, Jeon(09),Park, Jeon, Gale(17,18)

On-shell parton transport
Higher Twist
formalism

AMY
formalism

AdS/CFT
Chesler, Rajagopal(14, 15)
Pablos, et al.(15,16,17)

Diffusion into medium

𝒩 = 4 superYang-Mills

Y. Tachibana, QCD相転移やQGP生成のモデル化による重イオン衝突の時空発展の理解に向けた理論・実験共同研究会, September 24th, 2021

5

Medium response by weakly-coupled recoils

Zapp et al. (’13), Wang, Zhu (13), Luo, et al. (15,18), Park, Jeon, Gale (18), Cao, Majumder (18)

In-medium

Recoils-

Medium partons kicked out by jet partons2

Evolve as small-Q partons in jet shower
Infinite thermalization time (E > Emed)

Hole: Picked up energy and momentum-

Sampled from thermal medium-

Subtracted from final signal
μ

dp
dηdϕ

=
signal

μ

μ

dp
dηdϕ

jet shower

dp
−
dηdϕ

Y. Tachibana, QCD相転移やQGP生成のモデル化による重イオン衝突の時空発展の理解に向けた理論・実験共同研究会, September 24th, 2021

picked-up

6

Multi-stage description for jet evolution

JETSCAPE (17)Small-Q

Graph of parton shower generated by JETSCAPE
Large-Q(0.98,0.95,-0.00,21)68(7.10)(0.62,0.37,0.00,21)45(4.70)(2.18,2.13,0.22,21)(4.50,4.36,0.95,21)(1.81,1.50,0.00,21)(1.46,1.42,0.08,21)

13(1.30)(0.15,0.06,-0.00,-3)39(4.30)(0.93,0.67,0.00,21)(0.11,0.03,0.00,21)40(4.30)(0.60,0.59,-0.00,21)46(4.70)6(0.50)14(1.30)(1.06,1.04,-0.00,-3)72(7.90)(1.69,1.64,-0.00,-1)70(7.10)47(4.70)(8.89,8.64,2.97,21)41(4.30)(9.94,9.70,2.98,21)

21(1.60)(7.96,7.61,2.97,21)42(4.50)24(2.30)(1.73,1.62,0.00,3)(1.68,1.68,0.00,-2)(0.73,0.57,-0.00,-2)15(1.30)(9.68,9.31,2.97,21)

18(1.50)(0.35,0.33,-0.00,21)(0.62,0.16,-0.00,21)56(5.60)

43(4.50)
57(5.60)

25(2.30)(0.66,0.42,-0.00,3)22(1.60)

55(5.60)(2.48,2.36,0.29,21)(0.39,0.21,0.00,21)(1.97,1.79,0.08,21)(2.33,2.23,0.25,21)69(7.10)(0.29,0.13,-0.00,21)(0.72,0.69,0.00,-1)(0.90,0.56,-0.00,21)54(5.60)

71(7.90)44(4.50)

Virtuality
separation scale: Qsw
Switching between modules for parton
by parton(1.38,0.28,0.02,21)4(0.20)(13.61,13.18,2.98,21)2(0.10)(143.47,141.75,488.48,21)5(0.20)(141.50,140.13,386.92,21)7(0.50)11(0.90)(4.07,3.98,0.00,-3)20(1.50)12(0.90)26(2.30)

23(1.60)(0.09,0.05,-0.00,-3)(127.90,127.07,197.28,21)(144.85,141.91,836.69,21)

19(1.50)(126.41,126.06,82.46,21)(1.49,1.02,0.02,21)

16(1.40)(2.88,2.87,0.04,21)(123.53,123.30,51.20,21)30(3.00)

31(3.00)

48(5.00)(1.34,1.09,0.00,-2)(0.35,0.35,-0.00,-2)(120.59,120.46,27.13,21)(111.30,111.29,0.86,21)49(5.00)(117.83,117.78,9.82,21)17(1.40)

Large-Q ( > Qsw)0(0.00)(150.00,139.01,3175.39,21)

MATTER(5.15,4.37,0.20,21)

LBT(4.65,4.58,-0.00,-2)(4.12,4.00,-0.00,-2)(0.40,0.23,0.00,-2)9(0.70)36(4.10)(0.45,0.45,-0.00,21)37(4.10)(1.05,0.76,-0.00,21)(0.19,0.09,-0.00,3)62(6.10)MARTINI64(6.20)27(2.40)(2.33,2.15,0.00,2)33(3.20)28(2.40)(1.67,1.62,0.00,21)(0.56,0.56,0.00,-2)66(6.20)

On-shell parton transport
67(6.20)10(0.70)29(2.40)(1.38,1.04,-0.00,2)34(3.20)

Higher Twist
formalism35(3.20)(0.89,0.85,-0.00,21)(3.35,3.05,0.00,-2)(1.10,0.83,-0.00,-2)58(6.00)

AMY
formalism59(6.00)

60(6.00)(1.96,1.92,0.01,21)50(5.00)(6.54,6.51,0.03,21)74(8.80)(2.75,2.75,0.01,21)53(5.50)51(5.00)

63(6.10)Wang, Zhu(13), Luo, et al.(15,18)
Schenke, Gale, Jeon(09),Cao, et al.(16,17), He, et al.(18)(0.63,0.52,0.00,-2)Park,
Jeon, Gale(17,18)65(6.20)(3.13,3.11,0.00,21)(0.74,0.69,-0.00,-1)(1.18,0.57,-0.00,3)38(4.10)(0.59,0.59,-0.00,2)(0.21,0.19,0.00,-1)

Virtuality ordered splitting
Higher Twist
formalism(1.16,0.40,-0.00,-2)(0.57,0.44,-0.00,-2)

Majumder(13), Kordell, Majumder(17),Cao, Majumder(17)8(0.70)(1.48,0.28,0.20,21)3(0.10)61(6.10)

Small-Q ( < Qsw)

32(3.20)1(0.00)(0.51,0.44,0.02,21)52(5.50)73(8.80)

AdS/CFT
Chesler, Rajagopal(14, 15)
Pablos, et al.(15,16,17)

Black arrows: MATTER [Q > Qsw = 2 GeV]

Diffusion into medium

Blue arrows: LBT [Q < Qsw]
Green
Holes
= 4 super𝒩arrows:

Yang-Mills

Y. Tachibana, QCD相転移やQGP生成のモデル化による重イオン衝突の時空発展の理解に向けた理論・実験共同研究会, September 24th, 2021

7

Recent Update

JERSCAPE (in preperation)

Virtuality dependence
in jet quenching parameter q ̂
HTL q

10

-

q (Q2) parametrization: c1=10, c2=100

Kumar, Majumder, Shen (20)Hard Parton

JetScape Logo(s) Ia)q(Q )/ T

3

c)
`

2

1

Based on scale evolution of QGP constituent distribution

Q2`

JETSCAPE JETSCAPEb)Preliminaryd)

Je t Sca pe Je t Sca pe
`

10−1

`

QGP

1

Small-Q 2

Q2

Large-Q 2

E = 100 GeV, T = 0.2 GeV−210

5 10 15 20 25 30 35 40 45 50 55 60
22
Q (GeV )-Less interaction for larger virtuality partons-

Extend the model applicability to further broader scale

Y. Tachibana, QCD相転移やQGP生成のモデル化による重イオン衝突の時空発展の理解に向けた理論・実験共同研究会, September 24th, 2021

Jet and high-pT particle energy loss

JERSCAPE (in preperation)

Jet and high-pT particle energy loss for PbPb@5.02 TeV
Inclusive jet RAA
ATLAS [PLB 790 108 (2019)]

0.6
JetScape Logo(s) I
a)

0.2

c)

1

`

`

anti-k T R=0.4

JETSCAPE JETSCAPE|y |<2.8
Preliminaryjet
Je t Sca pe Je t Sca pe

b)

0.8
0.6

s=5.02 TeV, PbPb (0-10%)

0.4
0.2

d)

`

`

102

1

p jet (GeV)T

103

MATTER+LBT(q̂ without virtuality dependence)q , αfixs =0.25, Q =2GeV JS(MATTER+LBT)

1.2
±
RhAA

jet
RAA

0.8

−0.2

run

jet
run
fix
(q̂ without
virtuality
dependence)
qMATTER+LBT,α
Q0=2GeV
JS(MATTER+LBT)
s =0.25,
HTL
dependent-q)̂
qMATTER+LBT=q(Q2) type (virtuality
5, αfix
s =0.30

1

CMS [JHEP 1704 039 (2017)]

1.4

CMS, |η |<2.0

1.20Charged hadron RAA

pTjet

1.6

1.4

0.4

pT

0
10

0

HTL

a)

Logo(s)Iαfix =0.30
MATTER+LBT
dependent-q)̂
q=JetScapeq(Q2) type5,(virtuality
s
c)
`

`

JETSCAPE JETSCAPEb)Preliminaryd)

Je t Sca pe Je t Sca pe``1

s=5.02 TeV, PbPb (0-10%), |η|<1.0

p

T

102
(GeV)

Y. Tachibana, QCD相転移やQGP生成のモデル化による重イオン衝突の時空発展の理解に向けた理論・実験共同研究会, September 24th, 2021

pT

pT

pT
pT
pT

Jet and high-pT particle energy loss

JERSCAPE (in preperation)

Jet and high-pT particle energy loss for PbPb@5.02 TeV
Inclusive jet RAA

1.6

1.4

ATLAS [PLB 790 108 (2019)]

0.6
JetScape Logo(s) I
a)

0.2

−0.2

c)

1

`

`

anti-k T R=0.4

JETSCAPE JETSCAPE|y |<2.8
Preliminaryjet
Je t Sca pe Je t Sca pe

b)

0.8
0.6

s=5.02 TeV, PbPb (0-10%)

0.4
0.2

d)

`

`

MATTER+LBT(q̂ without virtuality dependence)q , αfixs =0.25, Q =2GeV JS(MATTER+LBT)

1.2
±
RhAA

jet
RAA

0.8

0

run

jet
run
fix
(q̂ without
virtuality
dependence)
qMATTER+LBT,α
Q0=2GeV
JS(MATTER+LBT)
s =0.25,
HTL
dependent-q)̂
qMATTER+LBT=q(Q2) type (virtuality
5, αfix
s =0.30

1

CMS [JHEP 1704 039 (2017)]

1.4

CMS, |η |<2.0

1.2

0.4

Charged hadron RAA

pTjet

0

HTL

a)

Logo(s)Iαfix =0.30
MATTER+LBT
dependent-q)̂
q=JetScapeq(Q2) type5,(virtuality
s
c)
`

`

JETSCAPE JETSCAPEb)Preliminaryd)

Je t Sca pe Je t Sca pe``1

s=5.02 TeV, PbPb (0-10%), |η|<1.0

0
10

102
- Very first success
p jet
(GeV)
in simultaneous description of data forand hadrons
p jet (GeV)

102

1

-

T

103

Significant effect of virtuality dependence in energy loss

T

Y. Tachibana, QCD相転移やQGP生成のモデル化による重イオン衝突の時空発展の理解に向けた理論・実験共同研究会, September 24th, 2021

pT

pT

pT

pT
pT
pT

Jet and high-pT particle energy loss
Jet and high-pT particle energy loss for PbPb@2.76 TeV

The same parameter set as 5.02 TeV is used
1.4

Inclusive
ATLAS[PLB 790jet
108R(2019)]AA

1.4

ChargedRAA
ATLAS[PLBhadron790 108 (2019)]

1.2

CMS, |η |<2.0

1.6
1.2

CMS, |η |<2.0

jet
run
ATLAS
790
108=2GeV(2019)]
fix =0.25,
̂
MATTER+LBT(
without
virtuality
dependence)q
JetScape Logo(s) I
qHTL, α[PLB
Q
JS(MATTER+LBT)
s
0
a)
c)
CMS, |η
2 |<2.0
fix
jettype (virtuality
dependent-q)̂
qMATTER+LBT
=
q
(Q
)
5,
α
=0.30
s
run

1

1.2

0.8

`

0.4
0.6

`

s=5.02 TeV, PbPb (0-10%)

0.2
0.4

anti-PbPb
k T R(0-10%)=0.4
s=5.02 TeV,|y |<2.8

10

2

10

3

jet

p

jet
p (GeV)
T
jet
T

(GeV)

10
10

3

0.2
−0.2

0 10
10

`

`

JETSCAPE JETSCAPEb)Preliminaryd)

Je t Sca pe Je t Sca pe``

s=5.02 TeV, PbPb (0-10%)
anti-k T R=0.4

0.4
0anti-k|yTjetR|<2.8=0.4
2MATTER+LBT
dependent-q)̂
q= q(Q2) type 5,(virtuality
αfix
s =0.30

0.6
0.2

1

0.2
0

−0.2

`

0

HTL

`

jet
Rh± RAA

jet
RAA

jet
RAA

dependent-q)̂
qMATTER+LBT=q(Q2) type (virtuality
5, αfix
s =0.30

0.6
0.8

−0.20

1.2
0.8
JETSCAPE JETSCAPE 1b)
0.6
Preliminary d)
Je t Sca pe Je t Sca pe
0.8
0.4

fix
(q̂ without
virtuality
dependence)
qMATTER+LBT,α
Q0=2GeV
JS(MATTER+LBT)
s =0.25,
HTL

1

1.4
1

AA

1.4

jet
run
fix =0.25,(q̂ without
virtuality
dependence)
qMATTER+LBT
,
α
Q
=2GeV
JS(MATTER+LBT)
CMS
[JHEP
1704
039
(2017)]s
JetScape Logo(s) I
HTL
0
2fix
a)
run
fix =0.30c)
MATTER+LBT(virtuality
dependent)̂
q
q
=
q
(Q
)
type
5,
α
̂
MATTER+LBT(
without
virtuality
dependence)qq , αs =0.25, Q =2GeV
JS(MATTER+LBT)s

1|y |<2.8jet2s=5.02 TeV, PbPb (0-10%), |η|<1.0

pp jet (GeV)
102
T

T

(GeV)

Y. Tachibana, QCD相転移やQGP生成のモデル化による重イオン衝突の時空発展の理解に向けた理論・実験共同研究会, September 24th, 2021

103Jet and high-pT particle energy loss

JERSCAPE (in preperation)

Jet and high-pT particle energy loss for AuAu@200 GeV

The same parameter set as 5.02 TeV is used
Charged jet RAA
ATLAS [PLB 790 108 (2019)]
a)

CMS,
|
η
|<2.0
STAR [PRC 102, 054913 (2020)]

1.2
0.8
1
0.6
0.8
0.4
0.6
0.2
0.4

ATLAS [PLB 790 108 (2019)]

1.2

CMS, |η |<2.0

jetrun
fix =0.25,(q̂ without
virtuality
dependence)
qMATTER+LBT,α
Q0=2GeV
JS(MATTER+LBT)s
HTL
dependent-q)̂
qMATTER+LBT=q(Q2) type (virtuality
5, αfix
s =0.30

c)

jet
`
run
fix =0.25,(q̂ without
virtuality
dependence)
qMATTER+LBT
,
α
Q
=2GeV
JS(MATTER+LBT)s
HTL0
fix =0.30
qq==qq(Q=0.30
JS(MATTER+LBT)
MATTER+LBT
dependentq)̂ b)Preliminary d)(Q2)2)type
type5,(virtuality5,αfixsαs

1`

JETSCAPE JETSCAPE0.8

Je t Sca pe Je t Sca pe
0.6
`

`

JetScape Logo(s) I

0.4

1

0.2

s=5.02 TeV, PbPb (0-10%)s=200 GeV, AuAu (0-10%)
antik
R
=0.4
T
anti-k T R=0.3
lead,ch|y |<2.8|η |<0.7, p
>5GeV
jet
jet
T

0.20
0
−0.2
−−0.2
0.4

JetScape Logo(s) I

1.4

jet
RAA

jet
jet,ch
RAA
RAA

1.41

Pion RAA

1010

2

15

20
jet
jet,chp p T (GeV)
(GeV)

25

0

a)

c)
`

JETSCAPE JETSCAPEb)Preliminaryd)

Je t Sca pe Je t Sca pe
`

−0.2
3

30 10

s=5.02 TeV, PbPb (0-10%)`10anti-k T R=0.4|y |<2.8

`

jet

3

2
1

p jet (GeV)T

T

Y. Tachibana, QCD相転移やQGP生成のモデル化による重イオン衝突の時空発展の理解に向けた理論・実験共同研究会, September 24th, 2021

10

Jet and high-pT particle energy loss

JERSCAPE (in preperation)

Jet and high-pT particle energy loss for AuAu@200 GeV

The same parameter set as 5.02 TeV is used
Charged jet RAA
ATLAS [PLB 790 108 (2019)]
CMS,
|
η
|<2.0
STAR [PRC 102, 054913 (2020)]

1.2
0.8
1
0.6
0.8
0.4
0.6
0.2
0.4

ATLAS [PLB 790 108 (2019)]

1.2

CMS, |η |<2.0

jetrun
fix =0.25,(q̂ without
virtuality
dependence)
qMATTER+LBT,α
Q0=2GeV
JS(MATTER+LBT)s
HTL
dependent-q)̂
qMATTER+LBT=q(Q2) type (virtuality
5, αfix
s =0.30c)

1`

JETSCAPE JETSCAPE0.8

Je t Sca pe Je t Sca pe
0.6
`

`

JetScape Logo(s) I

0.4

1

0.2

s=5.02 TeV, PbPb (0-10%)s=200 GeV, AuAu (0-10%)
antik
R
=0.4
T
anti-k T R=0.3
lead,ch|y |<2.8|η |<0.7, p
>5GeV
jet
jet
T

0
−0.2

-

a)

jet
`
run
fix =0.25,(q̂ without
virtuality
dependence)
qMATTER+LBT
,
α
Q
=2GeV
JS(MATTER+LBT)s
HTL0
fix =0.30
qq==qq(Q=0.30
JS(MATTER+LBT)
MATTER+LBT
dependentq)̂ b)Preliminary d)(Q2)2)type
type5,(virtuality5,αfix
sα
s

0.20

−−0.2
0.4

JetScape Logo(s) I

1.4

jet
RAA

jet
jet,ch
RAA
RAA

1.41

Pion RAA

1010

2

0
−0.2

10
20
25
30
jet
jet,chp (GeV)p
T (GeV)
Simultaneous descriptionof differentTc)`

s=5.02 TeV, PbPb (0-10%)`

JETSCAPE JETSCAPEb)Preliminaryd)

Je t Sca pe Je t Sca pe
`

3

15

a)

10

anti-k T R=0.4|y |<2.8

`

jet

3

2
1

p jet (GeV)TsNN with the same parameter set

Y. Tachibana, QCD相転移やQGP生成のモデル化による重イオン衝突の時空発展の理解に向けた理論・実験共同研究会, September 24th, 2021

10

Jet substructure

JERSCAPE (in preperation)

Jet Fragmentation Function-

pT

pT distribution of charged particle inside jetsdNch
1
D( pT) =Njet ∑
dp
Tjet

inI jetJetScape Logo(s)

a)

2.50

jet
126< pT <158 GeV

DPbPb(pT)/Dpp(pT)

2.25
2.00

jet
158<` pT <200 GeV`

JetScape Logo(s) I
a)

c)

jet
200<` pT <251 GeV

`

c)
`

`

JETSCAPE JETSCAPEJETSCAPE JETSCAPEJETSCAPE JETSCAP

b)

b)

b)

Preliminary`

MATTER+LBT(q̂ without virtuality dependence)

1.50

a)

c)

MATTER+LBT(virtuality dependent-q̂)Je t Sca pe

1.75

JetScape Logo(s) I

pTjet

1

d)

Preliminary

ATLAS
Je(2018)]
t Sca pe
Je t Sca pe [PRC 98, no.2, 024908
`

`

d)

Preliminaryd)

Je t Sca pe Je t Sca p

Je t Sca pe

`

`

`

Qswitch = 2 GeV|yjet | < 0.3
anti-kt , R = 0.4
ptrk
T > 1.0 GeV

1

1

Tc = 0.16 GeV

1.25
1.00
0.75

PbPb (0-10%), 5.02 TeV

0.50 010101

pT (GeV)102

101

pT (GeV)102

101

pT (GeV)Y. Tachibana, QCD相転移やQGP生成のモデル化による重イオン衝突の時空発展の理解に向けた理論・実験共同研究会, September 24th, 2021

102

12

Jet substructure

JERSCAPE (in preperation)

Jet Fragmentation Function-

pT

pT distribution of charged particle inside jetsdNch
1
D( pT) =Njet ∑
dp
Tjet

inI jetJetScape Logo(s)

a)

2.50

jet
126< pT <158 GeV

DPbPb(pT)/Dpp(pT)

2.25
2.00

jet
158<` pT <200 GeV`

JetScape Logo(s) I
a)

c)

jet
200<` pT <251 GeV

`

c)
`

`

JETSCAPE JETSCAPEJETSCAPE JETSCAPEJETSCAPE JETSCAP

b)

b)

b)

Preliminary`

MATTER+LBT(q̂ without virtuality dependence)

1.50

a)

c)

MATTER+LBT(virtuality dependent-q̂)Je t Sca pe

1.75

JetScape Logo(s) I

pTjet

1

d)

Preliminary

ATLAS
Je(2018)]
t Sca pe
Je t Sca pe [PRC 98, no.2, 024908
`

`

d)

Preliminaryd)

Je t Sca pe Je t Sca p

Je t Sca pe

`

`

`

Qswitch = 2 GeV|yjet | < 0.3
anti-kt , R = 0.4
ptrk
T > 1.0 GeV

1

1

Tc = 0.16 GeV

1.25
1.00
0.75

PbPb (0-10%), 5.02 TeV

0.50 010101

pT (GeV)102

101

pT (GeV)102

101

pT (GeV)Y. Tachibana, QCD相転移やQGP生成のモデル化による重イオン衝突の時空発展の理解に向けた理論・実験共同研究会, September 24th, 2021

102

12

Jet substructure

JERSCAPE (in preperation)

Jet Fragmentation Function-

pT

pT distribution of charged particle inside jetsdNch
1
D( pT) =Njet ∑
dp
Tjet

inI jetJetScape Logo(s)

a)

2.50

jet
126< pT <158 GeV

DPbPb(pT)/Dpp(pT)

2.25
2.00

jet
158<` pT <200 GeV`

JetScape Logo(s) I
a)

c)

jet
200<` pT <251 GeV

`

c)
`

`

JETSCAPE JETSCAPEJETSCAPE JETSCAPEJETSCAPE JETSCAP

b)

b)

b)

Preliminary`

MATTER+LBT(q̂ without virtuality dependence)

1.50

a)

c)

MATTER+LBT(virtuality dependent-q̂)Je t Sca pe

1.75

JetScape Logo(s) I

pTjet

d)

Preliminary

ATLAS
Je(2018)]
t Sca pe
Je t Sca pe [PRC 98, no.2, 024908
`

`

1

d)

Preliminaryd)

Je t Sca pe Je t Sca p

Je t Sca pe

`

`

`

Qswitch = 2 GeV|yjet | < 0.3
anti-kt , R = 0.4
ptrk
T > 1.0 GeV

1

1

Tc = 0.16 GeV

1.25
1.00
0.75

PbPb (0-10%),
5.02sensitive
TeV- Veryto medium effects during high virtuality phase

0.50 010-Enhancement
at small-pT
pT (GeV)
101

102

101

pT (GeV)102

101

pT (GeV)Y. Tachibana, QCD相転移やQGP生成のモデル化による重イオン衝突の時空発展の理解に向けた理論・実験共同研究会, September 24th, 2021

102

12

Jet substructure

JERSCAPE (in preperation)

Jet splitting function-

pT,2

Momentum fraction in the hardest splitting of jet (zg)zg =5pT,1 + pT,2

PbPb: 0-10%
anti-kt , R = 0.2
Charged jets
60 < pch,jet< 80 GeVT|¥ch,jet | < 0.7

5.02 TeV

4pT,1
PbPb: 0-10%
anti-kt , R = 0.4
Charged jets
80 < pch,jet< 100 GeVT|¥ch,jet | < 0.5

PbPb: 30-50%
anti-kt , R = 0.4
Charged jets
60 < pch,jet< 80 GeVT|¥ch,jet | < 0.5

3

1

dæ
æjet,inc dzg

min( pT,1, pT,2)2
MATTER (vacuum)

1a)

MATTER+LBT(virtuality dependent-q̂)

0

c)

a)

c)
`

`

Preliminaryb)d)Je t Sca pe Je t Sca pe`

1.2

JetScape Logo(s) I

JetScape Logo(s) I

PbPb[ALICE-PUBLIC-2020-006] JETSCAPE JETSCAPEJETSCAPE JETSCAPESoft Drop, zcut = 0.2, Ø = 0`b)

PbPb/pp

pp[ALICE-PUBLIC-2020-006] a)

JetScape Logo(s) I

d)

`

1

`

JETSCAPE JETSCAPEb)Preliminaryd)

Je t Sca pe Je t Sca pe

Je t Sca pe Je t Sca pe

`

`

`

Preliminary

c)

`

`

`

1

1

1.0
0.8
0.2

0.3

0.4

zg

0.2

0.3

0.4

zg

0.2

0.3

0.4

0.5

zg

Y. Tachibana, QCD相転移やQGP生成のモデル化による重イオン衝突の時空発展の理解に向けた理論・実験共同研究会, September 24th, 202113

Jet substructure

JERSCAPE (in preperation)

Jet splitting function-

pT,2

Momentum fraction in the hardest splitting of jet (zg)zg =5pT,1 + pT,2

PbPb: 0-10%
anti-kt , R = 0.2
Charged jets
60 < pch,jet< 80 GeVT|¥ch,jet | < 0.7

5.02 TeV

4pT,1
PbPb: 0-10%
anti-kt , R = 0.4
Charged jets
80 < pch,jet< 100 GeVT|¥ch,jet | < 0.5

PbPb: 30-50%
anti-kt , R = 0.4
Charged jets
60 < pch,jet< 80 GeVT|¥ch,jet | < 0.5

3

1

dæ
æjet,inc dzg

min( pT,1, pT,2)2
MATTER (vacuum)1

MATTER+LBT(virtuality dependent-q̂)a)

PbPb/pp

c)

a)

c)
`

`

Preliminaryb)d)Je t Sca pe Je t Sca pe`

1.2

JetScape Logo(s) I

JetScape Logo(s) I

PbPb[ALICE-PUBLIC-2020-006] JETSCAPE JETSCAPEJETSCAPE JETSCAPESoft Drop, zcut = 0.2, Ø = 0
`

b)

0

pp
[ALICE-PUBLIC-2020-006] a)

JetScape Logo(s) I

d)

`

1

`

JETSCAPE JETSCAPEb)Preliminaryd)

Je t Sca pe Je t Sca pe

Je t Sca pe Je t Sca pe

`

`

`

Preliminary

c)

`

`

`

1

1

1.0

Good agreement with experimental data
0.2
0.3no medium
0.4
0.2modification
0.3
0.4 in hardest
0.2
0.3splittings
0.4
0.5- Almost
-

0.8

zg

zg

zg

Y. Tachibana, QCD相転移やQGP生成のモデル化による重イオン衝突の時空発展の理解に向けた理論・実験共同研究会, September 24th, 2021

13

Fluid + Jet Approach for Hydro Response

Y. Tachibana, QCD相転移やQGP生成のモデル化による重イオン衝突の時空発展の理解に向けた理論・実験共同研究会, September 24th, 2021

Hydrodynamic Medium Response to Hard Partons

Hard
PartonIn-medium thermalization-

Drop to typical scale of thermalized medium constituents (E ∼ Emed)

Transition to hydrodynamic transport → Hydrodynamic medium response

Y. Tachibana, QCD相転移やQGP生成のモデル化による重イオン衝突の時空発展の理解に向けた理論・実験共同研究会, September 24th, 2021

15

Hydrodynamic Medium Response Effect on Jets{

Jet

Jet-induced flow in medium-

Transport momentum deposited by jet-

Modify particle emission around jet

Hydro medium response contribution-

Soft, spread out from jet-

Jet-correlated, cannot be subtracted-

Affect structures inside/around jet

Y. Tachibana, QCD相転移やQGP生成のモデル化による重イオン衝突の時空発展の理解に向けた理論・実験共同研究会, September 24th, 2021

16

Hydrodynamic Medium Response Effect on Jets{

Jet

Jet-induced flow in medium-

Transport momentum deposited by jet-

Modify particle emission around jet

Hydro medium response contribution-

Soft, spread out from jet-

Jet-correlated, cannot be subtracted-

Affect structures inside/around jet

Y. Tachibana, QCD相転移やQGP生成のモデル化による重イオン衝突の時空発展の理解に向けた理論・実験共同研究会, September 24th, 2021

16

Description of Hydro Medium Response to Jets
Medium fluid evolution with energy-momentum deposition
Coupled Jet-Fluid (YT, N.-B. Chang, G.-Y. Qin), CoLBT-hydro (W. Chen, X.-N. Wang, W. Zhao,…)
EPOS3-HQ (I. Karpenko,…), JETSCAPE (JETSCAPE), Hybrid+Linearized Hydro (X. Yao, D. Pablos,…),LEXUS+MUSIC (C.Shen, B. Schenke,…), DCCI2 (Y. Kanakubo, YT, T. Hirano,…), JAM (Y.Nara),…

Hydrodynamic transport of jet energy-momentum via thermal partons- Evolution together with the bulk QGP fluid-

Hydrodynamic equation with source termμνν∇μ Tmed(x) = Jjet(x)Energy-momentum tensor
of the QGP fluid

Energy and momentum
deposited into the fluidν- Source term Jjet constructed from jet-shower evolution calculation-

Bulk part particle with hydro response obtained via the Cooper-Frye

Y. Tachibana, QCD相転移やQGP生成のモデル化による重イオン衝突の時空発展の理解に向けた理論・実験共同研究会, September 24th, 2021

17

Causal Formulation for Source Term
YT, C. Shen and A. Majumder, arXiv:2001.08321

Energy-momentum deposition

Causal source profile

CoLBT-hydro (W. Chen, X.-N. Wang et al.)JETSCAPE (JETSCAPE)-

Soft partons

-

Holes

-

Relativistic diffusion equation∂∂+ τdiff 2 − Ddiff ∇2 j ν(x) = 0
∂t
[ ∂t
]
2

dep
E < Ecutdep
E < Ecutwith initial condition
ν
ν
(3)
⃗
⃗ )
j (t = tdep, x ) = pdep.δ ( x ⃗ − x dep
jiν(t ∼ ti, x)jiν(t = ti + Δt, x)

< cΔt
dep
Ecut : Energy scale for in-medium thermalization
Y. Tachibana, QCD相転移やQGP生成のモデル化による重イオン衝突の時空発展の理解に向けた理論・実験共同研究会, September 24th, 2021

18

Causal Formulation for Source Term
YT, C. Shen and A. Majumder, arXiv:2001.08321

Energy-momentum deposition

Causal source profile

CoLBT-hydro (W. Chen, X.-N. Wang et al.)JETSCAPE (JETSCAPE)-

Soft partons

-

Holes

-

Relativistic diffusion equation∂∂+ τdiff 2 − Ddiff ∇2 j ν(x) = 0
∂t
[ ∂t
]
2

p

μ

dep
E < Ecutpμdep
E < Ecutwith initial condition
ν
ν
(3)
⃗
⃗ )
j (t = tdep, x ) = pdep.δ ( x ⃗ − x dep
jiν(t ∼ ti, x)jiν(t = ti + Δt, x)

< cΔt
dep
Ecut : Energy scale for in-medium thermalization
Y. Tachibana, QCD相転移やQGP生成のモデル化による重イオン衝突の時空発展の理解に向けた理論・実験共同研究会, September 24th, 2021

18

Causal Formulation for Source Term
YT, C. Shen and A. Majumder, arXiv:2001.08321

Energy-momentum deposition

Causal source profile

CoLBT-hydro (W. Chen, X.-N. Wang et al.)JETSCAPE (JETSCAPE)-

Soft partons

-

Holes

−p

μ

−p

-

Relativistic diffusion equation∂∂+ τdiff 2 − Ddiff ∇2 j ν(x) = 0
∂t
[ ∂t
]
2

μ

p

μ

dep
E < Ecutpμdep
E < Ecutwith initial condition
ν
ν
(3)
⃗
⃗ )
j (t = tdep, x ) = pdep.δ ( x ⃗ − x dep
jiν(t ∼ ti, x)jiν(t = ti + Δt, x)

< cΔt
dep
Ecut : Energy scale for in-medium thermalization
Y. Tachibana, QCD相転移やQGP生成のモデル化による重イオン衝突の時空発展の理解に向けた理論・実験共同研究会, September 24th, 2021

18

Evolution of Hydrodynamic Medium Response
MATTER +MATTERLBT + Causal+ IdealHydro[Static
Brick,
Tbrick= 250 MeV]+ LBTDiff.
+ CausalDiff.
+ IdealHydro[Static
Brick]
YT, C. Shen, A.YT,
Majumder,C. Shen, arXiv:2001.08321
A. Majumder, arXiv:2001.08321-

Jet-Induced flow induced by a parton shower propagating in the x direction

z

x

y

z
Orange:Region with T > 250 MeV

x

Blue: Region
with T < 250 MeV
y
Red: Energetic PartonsY. Tachibana, QCD相転移やQGP生成のモデル化による重イオン衝突の時空発展の理解に向けた理論・実験共同研究会, September 24th,
24th, 2021
2021

19

Evolution of Hydrodynamic Medium Response
MATTER +MATTERLBT + Causal+ IdealHydro[Static
Brick,
Tbrick= 250 MeV]+ LBTDiff.
+ CausalDiff.
+ IdealHydro[Static
Brick]
YT, C. Shen, A.YT,
Majumder,C. Shen, arXiv:2001.08321
A. Majumder, arXiv:2001.08321-

Jet-Induced flow induced by a parton shower propagating in the x direction

z

x

y

z
Orange:Region with T > 250 MeV

x

Blue: Region
with T < 250 MeV
y
Red: Energetic PartonsY. Tachibana, QCD相転移やQGP生成のモデル化による重イオン衝突の時空発展の理解に向けた理論・実験共同研究会, September 24th,
24th, 2021
2021

19

Recoil vs Recoil+Hydro response: Static brick case
YT, C. Shen and A. Majumder, arXiv:2001.08321

x
Einit = 140 GeV

T = 0.25 GeV

Jet propagation in medium with static uniform initial condition- Energy loss by MATTER+LBT with recoils--

Medium evolves upto t = 10 fm

Y. Tachibana, QCD相転移やQGP生成のモデル化による重イオン衝突の時空発展の理解に向けた理論・実験共同研究会, September 24th, 2021

20

0

100

100

50

10000

5000

5000

0
0

°5000
0
10

2

5

0.0

0.5

(c-2) Expanding, Inward
µ
0
0.0

0.5

(a-2) Brick

0.0

5

0.5

0
0.5
(c-3) Expanding, Inward
µ
µ
15000
0
0.0

1

1

µ
µ

µ

2

2

1.0

µ

1.5

1.0

1.0

1.0

E(µ) (GeV)

415

1.5

µ

415

10
2

1.5

1.5

0
(a-3) Recoil+Hydro
Brick particlizedµStrong,
10000

Recoil+Hydro(w/o particlization)

Strong
RecoilWeak
VacuumVacuum

0

3

µ
3

100

2000

1000

0
0

-

50
0.0

￼

￼

dP
µ
μ
P
(θ)
=
dθ′
,
0
(b-2) Expanding, Outward
µ∫
dθ′
0

1

ΩE (µ)/Ωvac
E (µ)
3

100

-

50

3000

Angular structure of jet in brickJet energy, momentum and mass

2

0.5

θ

0.0

1

µ

1.0

1
dE
ρE(θ) =
,
µ
E(θ
=
1)
dθ
0
(b-3) Expanding, Outward
0.5

1.0

2

μ 1.5

3

Detailed substructure

2

1

1.5

µ

3

E(µ) (GeV)

CMS [JHEP 05, 006 (2018)]
PbPb 5.02 TeV (0-10%)

Px(µ) (GeV)

50
100
50

M (θ)0 =
P
(θ)P
(θ)
μ
(c-2) Expanding, Inward

2
μ

100

dPx
1
ρPx(θ) =
µ
P
(θ
=
1)
dθ
0
x
(c-3) Expanding, Inward

50

15000

0
0

sha1_base64="wFMugZqUWMZFtRU5ugq9AHBJ7EQ=">AAAB/HicdVDLSgMxFL1TX7W+qi7dBIvgapgpgl0W3LisYB/QDiWTZtrYJDMkGbEM9Rvc6tqduPVfXPonZtoK1seBkMM593LvPWHCmTae9+4UVlbX1jeKm6Wt7Z3dvfL+QUvHqSK0SWIeq06INeVM0qZhhtNOoigWIaftcHyR++1bqjSL5bWZJDQQeChZxAg2Vmr1QpHdTfvliudWvRzoN/Hd2e9VYIFGv/zRG8QkFVQawrHWXd9LTJBhZRjhdFrqpZommIzxkHYtlVhQHWSzbafoxCoDFMXKPmnQTP3ekWGh9USEtlJgM9I/vVz8y+umJqoFGZNJaqgk80FRypGJUX46GjBFieETSzBRzO6KyAgrTIwNaGlKKPJMvg5H/5NW1fU91786q9Rri3SKcATHcAo+nEMdLqEBTSBwAw/wCE/OvfPsvDiv89KCs+g5hCU4b5+tkJXj</latexit>

sha1_base64="wFMugZqUWMZFtRU5ugq9AHBJ7EQ=">AAAB/HicdVDLSgMxFL1TX7W+qi7dBIvgapgpgl0W3LisYB/QDiWTZtrYJDMkGbEM9Rvc6tqduPVfXPonZtoK1seBkMM593LvPWHCmTae9+4UVlbX1jeKm6Wt7Z3dvfL+QUvHqSK0SWIeq06INeVM0qZhhtNOoigWIaftcHyR++1bqjSL5bWZJDQQeChZxAg2Vmr1QpHdTfvliudWvRzoN/Hd2e9VYIFGv/zRG8QkFVQawrHWXd9LTJBhZRjhdFrqpZommIzxkHYtlVhQHWSzbafoxCoDFMXKPmnQTP3ekWGh9USEtlJgM9I/vVz8y+umJqoFGZNJaqgk80FRypGJUX46GjBFieETSzBRzO6KyAgrTIwNaGlKKPJMvg5H/5NW1fU91786q9Rri3SKcATHcAo+nEMdLqEBTSBwAw/wCE/OvfPsvDiv89KCs+g5hCU4b5+tkJXj</latexit><latexit

sha1_base64="wFMugZqUWMZFtRU5ugq9AHBJ7EQ=">AAAB/HicdVDLSgMxFL1TX7W+qi7dBIvgapgpgl0W3LisYB/QDiWTZtrYJDMkGbEM9Rvc6tqduPVfXPonZtoK1seBkMM593LvPWHCmTae9+4UVlbX1jeKm6Wt7Z3dvfL+QUvHqSK0SWIeq06INeVM0qZhhtNOoigWIaftcHyR++1bqjSL5bWZJDQQeChZxAg2Vmr1QpHdTfvliudWvRzoN/Hd2e9VYIFGv/zRG8QkFVQawrHWXd9LTJBhZRjhdFrqpZommIzxkHYtlVhQHWSzbafoxCoDFMXKPmnQTP3ekWGh9USEtlJgM9I/vVz8y+umJqoFGZNJaqgk80FRypGJUX46GjBFieETSzBRzO6KyAgrTIwNaGlKKPJMvg5H/5NW1fU91786q9Rri3SKcATHcAo+nEMdLqEBTSBwAw/wCE/OvfPsvDiv89KCs+g5hCU4b5+tkJXj</latexit><latexit

sha1_base64="wFMugZqUWMZFtRU5ugq9AHBJ7EQ=">AAAB/HicdVDLSgMxFL1TX7W+qi7dBIvgapgpgl0W3LisYB/QDiWTZtrYJDMkGbEM9Rvc6tqduPVfXPonZtoK1seBkMM593LvPWHCmTae9+4UVlbX1jeKm6Wt7Z3dvfL+QUvHqSK0SWIeq06INeVM0qZhhtNOoigWIaftcHyR++1bqjSL5bWZJDQQeChZxAg2Vmr1QpHdTfvliudWvRzoN/Hd2e9VYIFGv/zRG8QkFVQawrHWXd9LTJBhZRjhdFrqpZommIzxkHYtlVhQHWSzbafoxCoDFMXKPmnQTP3ekWGh9USEtlJgM9I/vVz8y+umJqoFGZNJaqgk80FRypGJUX46GjBFieETSzBRzO6KyAgrTIwNaGlKKPJMvg5H/5NW1fU91786q9Rri3SKcATHcAo+nEMdLqEBTSBwAw/wCE/OvfPsvDiv89KCs+g5hCU4b5+tkJXj</latexit><latexit

latexit<150

4ΩE (µ)/Ωvac
E (µ)

150
(b-1) Expanding, Outward4ΩPx (µ)/Ωvac
Px (µ)

0(a-1) BrickΩPx (µ)/Ωvac
Px (µ)

50

Px(µ) (GeV)E(c-1) Expanding, Inward-

Backward suppression due
10000 to particlization (Recoil+Hydro)-

Backward suppression due
toholes(Recoil)5000
M 2(µ) (GeV2)

50
Px

100

x

100

vacvacΩ(µ)/Ω
E
ΩE (µ)/Ω E(µ)(µ)

150vac vacΩ(µ)/Ω
PxPx (µ)
ΩP (µ)/Ω (µ)E(µ)(GeV)
(GeV)
E(µ)
150

M 2(µ) (GeV2)

Px(µ)
(GeV)
Px(µ)
(GeV)

5

2
M 2(GeV(µ) (GeVM 2(µ)) 2)

5

Recoil vs Recoil+Hydro response: Static brick case
YT, C. Shen and A. Majumder, arXiv:2001.08321(c-1) Expanding, Inward

2

x

0.0
0.5

0.0

1

µ

µ

0.5

2

1.0

1.0

1.5

θ
Einit = 140 GeV

2

1.5

µ

Strong, particlized
Strong
Weak
Vacuum3µ

Y. Tachibana, QCD相転移やQGP生成のモデル化による重イオン衝突の時空発展の理解に向けた理論・実験共同研究会, September 24th, 2021

21

Recoil vs Recoil+Hydro response: Expanding medium case
YT, C. Shen and A. Majumder, arXiv:2001.08321

xy

Tcenter = 0.5 GeV-

Jet propagation in medium with oblate 3-D Gaussian (σT = 1.5 fm, σz = 0.75 fm)

Medium size is chosen to reproduce the full simulation results- Radial flow following the jet propagation-

Y. Tachibana, QCD相転移やQGP生成のモデル化による重イオン衝突の時空発展の理解に向けた理論・実験共同研究会, September 24th, 2021

22

0
0.5

x

µ

0

1.0

50

5

1.0

1000

0
0
0.0 100

3g, Outwardµ 2000

1
0

1

0.0

5000

0
2
3
0
1
2
3
Y. Tachibana,
QCD相転移やQGP生成のモデル化による重イオン衝突の時空発展の理解に向けた理論・実験共同研究会,
September 24th, 2021

µ
0.5

µ

0.5

0
50
(b-3)
Expanding,Outward
µ
3000

µ

10000

2
1.0

(b-2) Expanding, Outward
µ

50
2

100
1.0

CMS [JHEP 05, 006 (2018)]
PbPb 5.02 TeV (0-10%)(c-1) Expanding, Inward

4

1.5

4

1.5

E(µ) (GeV)ΩE (µ)/Ωvac
E (µ)

0.0
0.5

2

0.0
0.5

0
(c-3) Expanding, Inwardµ15000

Recoil+Hydro(w/o particlization)

Strong
RecoilWeak
VacuumVacuum

3

100

-

50

1.5
0.0

￼

￼

1.5

2

µ

dP
μ
P
(θ)
=
dθ′
,
0
(c-2) Expanding, Inward
µ∫
dθ′
0

µ 100

1.0

(c-2) Expanding, Inward
µ

-

50

15000

1.0
1.5

10000

Recoil+Hydro
Strong, particlized

5000

0
0

0.5

θ

0.0

-

Collimation due to push by the radial flow-

Small effect of particlization (Recoil+Hydro)-

Clear difference at very large angle region (θ > 1)
1

µ

1
dE
ρE(θ) =
,
µ
E(θ
=
1)
dθ
0
(c-3) Expanding, Inward
0.5

2

1.0

1.0

2

1.5

µ

sha1_base64="wFMugZqUWMZFtRU5ugq9AHBJ7EQ=">AAAB/HicdVDLSgMxFL1TX7W+qi7dBIvgapgpgl0W3LisYB/QDiWTZtrYJDMkGbEM9Rvc6tqduPVfXPonZtoK1seBkMM593LvPWHCmTae9+4UVlbX1jeKm6Wt7Z3dvfL+QUvHqSK0SWIeq06INeVM0qZhhtNOoigWIaftcHyR++1bqjSL5bWZJDQQeChZxAg2Vmr1QpHdTfvliudWvRzoN/Hd2e9VYIFGv/zRG8QkFVQawrHWXd9LTJBhZRjhdFrqpZommIzxkHYtlVhQHWSzbafoxCoDFMXKPmnQTP3ekWGh9USEtlJgM9I/vVz8y+umJqoFGZNJaqgk80FRypGJUX46GjBFieETSzBRzO6KyAgrTIwNaGlKKPJMvg5H/5NW1fU91786q9Rri3SKcATHcAo+nEMdLqEBTSBwAw/wCE/OvfPsvDiv89KCs+g5hCU4b5+tkJXj</latexit>

sha1_base64="wFMugZqUWMZFtRU5ugq9AHBJ7EQ=">AAAB/HicdVDLSgMxFL1TX7W+qi7dBIvgapgpgl0W3LisYB/QDiWTZtrYJDMkGbEM9Rvc6tqduPVfXPonZtoK1seBkMM593LvPWHCmTae9+4UVlbX1jeKm6Wt7Z3dvfL+QUvHqSK0SWIeq06INeVM0qZhhtNOoigWIaftcHyR++1bqjSL5bWZJDQQeChZxAg2Vmr1QpHdTfvliudWvRzoN/Hd2e9VYIFGv/zRG8QkFVQawrHWXd9LTJBhZRjhdFrqpZommIzxkHYtlVhQHWSzbafoxCoDFMXKPmnQTP3ekWGh9USEtlJgM9I/vVz8y+umJqoFGZNJaqgk80FRypGJUX46GjBFieETSzBRzO6KyAgrTIwNaGlKKPJMvg5H/5NW1fU91786q9Rri3SKcATHcAo+nEMdLqEBTSBwAw/wCE/OvfPsvDiv89KCs+g5hCU4b5+tkJXj</latexit><latexit

sha1_base64="wFMugZqUWMZFtRU5ugq9AHBJ7EQ=">AAAB/HicdVDLSgMxFL1TX7W+qi7dBIvgapgpgl0W3LisYB/QDiWTZtrYJDMkGbEM9Rvc6tqduPVfXPonZtoK1seBkMM593LvPWHCmTae9+4UVlbX1jeKm6Wt7Z3dvfL+QUvHqSK0SWIeq06INeVM0qZhhtNOoigWIaftcHyR++1bqjSL5bWZJDQQeChZxAg2Vmr1QpHdTfvliudWvRzoN/Hd2e9VYIFGv/zRG8QkFVQawrHWXd9LTJBhZRjhdFrqpZommIzxkHYtlVhQHWSzbafoxCoDFMXKPmnQTP3ekWGh9USEtlJgM9I/vVz8y+umJqoFGZNJaqgk80FRypGJUX46GjBFieETSzBRzO6KyAgrTIwNaGlKKPJMvg5H/5NW1fU91786q9Rri3SKcATHcAo+nEMdLqEBTSBwAw/wCE/OvfPsvDiv89KCs+g5hCU4b5+tkJXj</latexit><latexit

sha1_base64="wFMugZqUWMZFtRU5ugq9AHBJ7EQ=">AAAB/HicdVDLSgMxFL1TX7W+qi7dBIvgapgpgl0W3LisYB/QDiWTZtrYJDMkGbEM9Rvc6tqduPVfXPonZtoK1seBkMM593LvPWHCmTae9+4UVlbX1jeKm6Wt7Z3dvfL+QUvHqSK0SWIeq06INeVM0qZhhtNOoigWIaftcHyR++1bqjSL5bWZJDQQeChZxAg2Vmr1QpHdTfvliudWvRzoN/Hd2e9VYIFGv/zRG8QkFVQawrHWXd9LTJBhZRjhdFrqpZommIzxkHYtlVhQHWSzbafoxCoDFMXKPmnQTP3ekWGh9USEtlJgM9I/vVz8y+umJqoFGZNJaqgk80FRypGJUX46GjBFieETSzBRzO6KyAgrTIwNaGlKKPJMvg5H/5NW1fU91786q9Rri3SKcATHcAo+nEMdLqEBTSBwAw/wCE/OvfPsvDiv89KCs+g5hCU4b5+tkJXj</latexit><latexit

latexit

μ 1.5
<

150

4

2

ΩE (µ)/Ωvac
E (µ)(b-1) Expanding, OutwardΩPx (µ)/Ωvac
Px (µ)CMS [JHEP 05, 006 (2018)]
PbPb 5.02 TeV (0-10%)1Px(µ) (GeV)50150

ΩPx (µ)/Ωvac
Px (µ)

2

vac
Ω
(µ)/Ω
E
E(µ) (GeV) E (µ)
3

M 2(µ) (GeV2)

µ
100

vac
Ω
(µ)/Ω
P
x
P (µ) (GeV) Px (µ)

100

M (µ) (GeV )0.5

2g, Outward

E(µ) (GeV)150

2g, Outwardµ

Px(µ) (GeV)0M 2(µ) (GeV2)

5

Recoil vs Recoil+Hydro response: Expanding medium case
YT, C. Shen and A. Majumder, arXiv:2001.08321

Angular structure of jet in oblate medium(c-1) Expanding, Inward

Jet energy, momentum and massx

M 2(θ) = P μ(θ)Pμ(θ)θ

1.5

Detailed substructure

Einit = 140 GeV

4

dPx
1
ρPx(θ) =Px(θ = 1) dθ

Strong, particlized
Strong
Weak
Vacuum3µ

23Summary

Jet evolution in JETSCAPE-

Multi-stage description by switching modules with virtuality
2

Further extension by Q -dependence in jet quenching parameter q̂

-Simultaneous description of jet and high-pT particle energy loss in various-

Further systematic studies for jet substructuresNN

Hydrodynamic medium response
Transport of thermalized jet energy and momentum- Described by hydrodynamic equation with source terms- Difference from recoils led by background flow-

Y. Tachibana, QCD相転移やQGP生成のモデル化による重イオン衝突の時空発展の理解に向けた理論・実験共同研究会, September 24th, 2021

24

BackupY. Tachibana, QCD相転移やQGP生成のモデル化による重イオン衝突の時空発展の理解に向けた理論・実験共同研究会, September 24th, 2021

Causal Formulation for Source Term
YT, C. Shen and A. Majumder, arXiv:2001.08321

0(a-3) Brick
10000

5

µ

1.0

M 2(µ) (GeV2)

µ

µ

1.0

3

2

µ

ΩE (µ)/Ωvac
E (µ)E(µ) (GeV)

1000

3

1.0

1.5

1.0

1.5

4

50

2
0.0

10000

1

µ

0.5

µ

0
(c-3) Expanding, Inward
µ
15000

2000

0
0

100

0.5

(c-2) Expanding, Inward
µ

1.5

0
µ
3000 (b-3) Expanding, Outward

0

2

1
0.5

2
0.0

0

2

50

4

1.5

3

0.0

µ

1

µ

1.0

(b-2) Expanding, Outward
µ

1.5

5000

°5000
0

100

0.5

50

(c-1) Expanding, InwardΩPx (µ)/Ωvac
Px (µ)

0

15

0.5

1
0.0

µ

0
0.0

50

1.5

10

50

2

100

Px(µ) (GeV)µ

1.0

CMS [JHEP 05, 006 (2018)]
PbPb 5.02 TeV (0-10%)

M 2(µ) (GeV2)

0.5

3

150

ΩPx (µ)/Ωvac
Px (µ)

0
0.0

Px(µ) (GeV)

Px(µ) (GeV)1005(a-2) BrickΩPx (µ)/Ωvac
Px (µ)

50

100(b-1) Expanding, Outward
ΩE (µ)/Ωvac
E (µ)

15
10

0

M 2(µ) (GeV2)

E(µ) (GeV)

100

150

ΩE (µ)/Ωvac
E (µ)E(µ) (GeV)150(a-1) Brick

5000
0
0

Strong, particlized
Strong
Weak
Vacuum1

2

3µ

Y. Tachibana, QCD相転移やQGP生成のモデル化による重イオン衝突の時空発展の理解に向けた理論・実験共同研究会, September 24th, 2021

26