Heavy flavor in heavy ion collisions with ATLAS

Qipeng Hu

University of Colorado Boulder

The 2018 RHIC/AGS Annual Users’ Meeting

Illuminating the QCD Landscape

RHIC/AGS Users’ Executive Committee

Jim Thomas, UEC Chair jhthomas@lbl.gov

Rosi Reed, Chair Elect rosijreed@lehigh.edu

Lijuan Ruan, Past Chair ruan@bnl.gov

Conference Coordinator

Kelly Guiffreda guiffreda@bnl.gov

Registration and Program: www.bnl.gov/aum2018

June 12–15, 2018 Brookhaven National LaboratoryTopical workshops: June 12-13 Plenary sessions: June 14–15

The 2018 RHIC/AGS Annual Users’ Meeting

ATLAS detector

�2

Inner Detector ⎢η⎥ < 2.5Calorimeter ⎢η⎥ < 4.9Muon Spectrometer ⎢η⎥ < 2.7

By J. Goodson

ID + MS

ID

HF muon

pT [GeV]14 4010420

CharmoniaBottomonia

D0

Heavy flavor programs in ATLAS

�3

HF muon D meson Charmonia Bottomonia

Large systemPb+Pb

Yield ✔ ✔

Flow ✔ ✔

Small systemp+Pb

Yield ✔ ✔ ✔

Azimuthal anisotropy ✔ ✔

• HF muon yield and flow in 2.76 TeV Pb+Pb• HF muon flow in 8.16 TeV p+Pb• D meson yield and flow in 8.16 TeV p+Pb• Charmonia yield in 5.02 TeV Pb+Pb• J/𝝍 flow in 5.02 TeV Pb+Pb• Charmonia and bottomonia yield in 5.02 TeV p+Pb

prelim.

prelim.

prelim.

Heavy flavor in large system

�4

• Due to large mass, heavy flavors (HF) quarks are most created from hard scattering at early stage of heavy ion collisions, pQCD calculable

• Open heavy flavor gives direct access to the interaction between QGP and heavy quarks, supplements description of quark energy loss in QGP

• Hidden heavy flavor (quarkonia) shows then interplay of color screening and recombination, variety of quarkonium states production provide binding energy dependence information

HF muon RAA in Pb+Pb

�5

• HF muon RAA measured in 5 centrality slices covers 0-60%• Significant suppression ~ 0.4 in most central• Strong centrality dependence. No significant dependence on pT

arXiv: 1805.05220

Not

revi

ewed

,for

inte

rnal

circ

ulat

ion

only

ATLAS DRAFT

[GeV]Tp

4 6 8 10 12 14

]-1 G

eV×

[nb

ηd Tpdµ

HF

N2 d ev

tN1

〉AA

T〈1

1

10

210

310

410

510

610

710

810

910 )5 10× 0 - 10 % ( )4 10×10 - 20 % ( )3 10×20 - 30 % ( )2 10×30 - 40 % ( )1 10×40 - 60 % (

ppppscaled

| < 1η|

-1, 570 nbpp

-1Pb+Pb, 0.14 nb = 2.76 TeVNNsATLAS

Figure 6: The pT dependence of the measured Pb+Pb heavy-flavor muon di�erential per-event yields for di�erentcentrality intervals scaled by the corresponding hTAAi. Also shown is the measured pp heavy-flavor muon di�erentialcross-section. For clarity, the results for the di�erent centralities are multiplied by scale factors that are indicated inthe legend. The pp cross-section is re-plotted multiple times, as dashed lines, multiplied by these scale factors, forcomparison with the results for the di�erent Pb+Pb centralities. The error bars and shaded bands represent statisticaland systematic uncertainties, respectively, and in many cases are too small to be seen.

[GeV] Tp

4 6 8 10 12 14

AA

R

0

0.5

1

-1Pb+Pb, 0.14 nb -1, 570 nbpp

| < 1η|

ATLAS = 2.76 TeVNNs

0-10%20-30%40-60%

-1Pb+Pb, 0.14 nb -1, 570 nbpp

| < 1η|

ATLAS = 2.76 TeVNNs

0-10%20-30%40-60%

-1Pb+Pb, 0.14 nb -1, 570 nbpp

| < 1η|

ATLAS = 2.76 TeVNNs

0-10%20-30%40-60%

-1Pb+Pb, 0.14 nb -1, 570 nbpp

| < 1η|

ATLAS = 2.76 TeVNNs

0-10%20-30%40-60%

-1Pb+Pb, 0.14 nb -1, 570 nbpp

| < 1η|

ATLAS = 2.76 TeVNNs

0-10%20-30%40-60%

[GeV] Tp

4 6 8 10 12 14

AA

R

0

0.5

1

-1Pb+Pb, 0.14 nb -1, 570 nbpp

| < 1η|

ATLAS = 2.76 TeVNNs

10-20%30-40%

-1Pb+Pb, 0.14 nb -1, 570 nbpp

| < 1η|

ATLAS = 2.76 TeVNNs

10-20%30-40%

-1Pb+Pb, 0.14 nb -1, 570 nbpp

| < 1η|

ATLAS = 2.76 TeVNNs

10-20%30-40%

-1Pb+Pb, 0.14 nb -1, 570 nbpp

| < 1η|

ATLAS = 2.76 TeVNNs

10-20%30-40%

Figure 7: The measured Pb+Pb heavy-flavor muon RAA as a function of pT. For clarity, the centrality intervals aresplit between the two panels. The left panel shows results for the 0–10%, 20–30%, and 40–60% centrality intervalswhile the right panel shows results for the 10–20% and 30–40% intervals. The error bars represent statisticaluncertainties. The boxes indicate theoretical uncertainties of hTAAi. The shaded bands represent the experimentalsystematic uncertainties.

12th May 2018 – 19:26 21

HF muon RAA in Pb+Pb

�6 arXiv: 1805.05220

• Similar HF muon and D0 RAA at low pT, both larger than light hadrons

• Less suppression of HF muon at higher pT, hint of less suppression of B-hadrons

• Separation b/c contribution in HF muon is crucial for future measurements

Not

revi

ewed

,for

inte

rnal

circ

ulat

ion

only

ATLAS DRAFT

[GeV] T

p4 6 8 10 12 14

AA

R

0

0.2

0.4

0.6

0.8

1 0-10%

for 0-5%AAR ±h

[GeV] Tp

4 6 8 10 12 14

AA

R

0

0.2

0.4

0.6

0.8

1 20-30% [GeV] T

p4 6 8 10 12 14

AA

R

0

0.2

0.4

0.6

0.8

1 10-20%

[GeV] Tp

4 6 8 10 12 14

AA

R

0

0.2

0.4

0.6

0.8

1 30-40%5 5.2 5.4 5.6 5.8 6

0

0.2

0.4

0.6

0.8

1

| < 1η|

-1, 570 nbpp-1Pb+Pb, 0.14 nb

= 2.76 TeVNNs InternalATLAS

[GeV] Tp

4 6 8 10 12 14

AA

R

0

0.2

0.4

0.6

0.8

1 40-60%

|<1η : |±µ ATLAS HF |<1y : |0 CMS D|<2η : |±h ATLAS

Figure 9: Comparison of the Pb+Pb heavy-flavor muon RAA measured in this analysis to the RAA for inclusivecharged hadrons from ATLAS and the RAA for identified D0 mesons from the CMS Collaboration. The error barsrepresent systematic and statistical uncertainties added in quadrature. The hTAAi errors are identical between thethree measurements and are excluded from the comparison. The inclusive charged hadron RAA values shown in thetop left panel are for the 0–5% centrality interval.

25th April 2018 – 14:47 23

Not

revi

ewed

,for

inte

rnal

circ

ulat

ion

only

ATLAS DRAFT

[GeV] T

p4 6 8 10 12 14

AA

R

0

0.2

0.4

0.6

0.8

1 0-10%

for 0-5%AAR ±h

[GeV] Tp

4 6 8 10 12 14

AA

R

0

0.2

0.4

0.6

0.8

1 20-30% [GeV] T

p4 6 8 10 12 14

AA

R

0

0.2

0.4

0.6

0.8

1 10-20%

[GeV] Tp

4 6 8 10 12 14

AA

R

0

0.2

0.4

0.6

0.8

1 30-40%5 5.2 5.4 5.6 5.8 6

0

0.2

0.4

0.6

0.8

1

| < 1η|

-1, 570 nbpp-1Pb+Pb, 0.14 nb

= 2.76 TeVNNs InternalATLAS

[GeV] Tp

4 6 8 10 12 14

AA

R

0

0.2

0.4

0.6

0.8

1 40-60%

|<1η : |±µ ATLAS HF |<1y : |0 CMS D|<2η : |±h ATLAS

Figure 9: Comparison of the Pb+Pb heavy-flavor muon RAA measured in this analysis to the RAA for inclusivecharged hadrons from ATLAS and the RAA for identified D0 mesons from the CMS Collaboration. The error barsrepresent systematic and statistical uncertainties added in quadrature. The hTAAi errors are identical between thethree measurements and are excluded from the comparison. The inclusive charged hadron RAA values shown in thetop left panel are for the 0–5% centrality interval.

25th April 2018 – 14:47 23

Not

revi

ewed

,for

inte

rnal

circ

ulat

ion

only

ATLAS DRAFT

[GeV] T

p4 6 8 10 12 14

AA

R

0

0.2

0.4

0.6

0.8

1 0-10%

for 0-5%AAR ±h

[GeV] Tp

4 6 8 10 12 14

AA

R

0

0.2

0.4

0.6

0.8

1 20-30% [GeV] T

p4 6 8 10 12 14

AA

R

0

0.2

0.4

0.6

0.8

1 10-20%

[GeV] Tp

4 6 8 10 12 14

AA

R

0

0.2

0.4

0.6

0.8

1 30-40%5 5.2 5.4 5.6 5.8 6

0

0.2

0.4

0.6

0.8

1

| < 1η|

-1, 570 nbpp-1Pb+Pb, 0.14 nb

= 2.76 TeVNNs InternalATLAS

[GeV] Tp

4 6 8 10 12 14

AA

R

0

0.2

0.4

0.6

0.8

1 40-60%

|<1η : |±µ ATLAS HF |<1y : |0 CMS D|<2η : |±h ATLAS

Figure 9: Comparison of the Pb+Pb heavy-flavor muon RAA measured in this analysis to the RAA for inclusivecharged hadrons from ATLAS and the RAA for identified D0 mesons from the CMS Collaboration. The error barsrepresent systematic and statistical uncertainties added in quadrature. The hTAAi errors are identical between thethree measurements and are excluded from the comparison. The inclusive charged hadron RAA values shown in thetop left panel are for the 0–5% centrality interval.

25th April 2018 – 14:47 23

HF muon 𝒗2 in Pb+Pb

�7

• 𝒗2 is extracted based on event plane method

• Checked with scalar product• Combinatorial background is

subtracted by momentum imbalance template fits

• Significant 𝒗2 observed 0 - 60% centrality

arXiv: 1805.05220

Not

revi

ewed

,for

inte

rnal

circ

ulat

ion

only

ATLAS DRAFT

[GeV]T

p4 6 8 10 12 14

2v

0

0.05

0.1Event planeScalar product

0-10 %

[GeV]T

p4 6 8 10 12 14

2v

0

0.05

0.1 10-20 %

[GeV]Tp

4 6 8 10 12 14

2v

0

0.05

0.1 30-40 %

5 5.2 5.4 5.6 5.8 6

0

0.05

0.1

| < 2η|-1Pb+Pb, 0.14 nb

= 2.76 TeVNNsATLAS

[GeV]T

p4 6 8 10 12 14

2v

0

0.05

0.1 20-30 %

[GeV]Tp

4 6 8 10 12 14

2v

0

0.05

0.1 40-60 %

Figure 10: The pT dependence of the Pb+Pb heavy-flavor muon v2. Results are shown for both the EP and SPmethods. Each panel represents a di�erent centrality interval. The error bars and shaded bands represent statisticaland total uncertainties, respectively, and are shown only for the EP v2. The horizontal dashed lines indicate v2 = 0.

12th May 2018 – 19:26 25

HF muon RAA and 𝒗2 vs calculations

�8

Not

revi

ewed

,for

inte

rnal

circ

ulat

ion

only

ATLAS DRAFT

[GeV] T

p4 6 8 10 12 14

AA

R

0

0.2

0.4

0.6

0.8

1 0-10%

[GeV] Tp

4 6 8 10 12 14

AA

R

0

0.2

0.4

0.6

0.8

1 20-30%TAMU values for 20-40%

[GeV] T

p4 6 8 10 12 14

AA

R

0

0.2

0.4

0.6

0.8

1 10-20%

[GeV] Tp

4 6 8 10 12 14

AA

R0

0.2

0.4

0.6

0.8

1 30-40%TAMU values for 20-40%

5 5.2 5.4 5.6 5.8 6

0

0.2

0.4

0.6

0.8

1

| < 1η|

-1, 570 nbpp-1Pb+Pb, 0.14 nb

= 2.76 TeVNNs InternalATLAS

DataDABModTAMU

[GeV] Tp

4 6 8 10 12 14

AA

R

0

0.2

0.4

0.6

0.8

1 40-60%

Figure 14: Comparison of the measured heavy-flavor muon RAA in Pb+Pb collisions with the values predicted fromthe TAMU transport model and the DABMod model. Each panel represents a di�erent centrality interval. For the20–30% and 30–40% centrality intervals, the plotted TAMU values correspond to the 20–40% centrality interval.For the data, the error bars represent statistical uncertainties, the shaded bands represent the experimental systematicuncertainties, and the boxes indicate theoretical uncertainties from hTAAi. For the model calculations the bandsindicate the theoretical systematic uncertainties.

25th April 2018 – 14:47 30

Not

revi

ewed

,for

inte

rnal

circ

ulat

ion

only

ATLAS DRAFT

[GeV] T

p4 6 8 10 12 14

AA

R

0

0.2

0.4

0.6

0.8

1 0-10%

[GeV] Tp

4 6 8 10 12 14

AA

R

0

0.2

0.4

0.6

0.8

1 20-30%TAMU values for 20-40%

[GeV] T

p4 6 8 10 12 14

AA

R

0

0.2

0.4

0.6

0.8

1 10-20%

[GeV] Tp

4 6 8 10 12 14

AA

R

0

0.2

0.4

0.6

0.8

1 30-40%TAMU values for 20-40%

5 5.2 5.4 5.6 5.8 6

0

0.2

0.4

0.6

0.8

1

| < 1η|

-1, 570 nbpp-1Pb+Pb, 0.14 nb

= 2.76 TeVNNs InternalATLAS

DataDABModTAMU

[GeV] Tp

4 6 8 10 12 14

AA

R

0

0.2

0.4

0.6

0.8

1 40-60%

Figure 14: Comparison of the measured heavy-flavor muon RAA in Pb+Pb collisions with the values predicted fromthe TAMU transport model and the DABMod model. Each panel represents a di�erent centrality interval. For the20–30% and 30–40% centrality intervals, the plotted TAMU values correspond to the 20–40% centrality interval.For the data, the error bars represent statistical uncertainties, the shaded bands represent the experimental systematicuncertainties, and the boxes indicate theoretical uncertainties from hTAAi. For the model calculations the bandsindicate the theoretical systematic uncertainties.

25th April 2018 – 14:47 30

Not

revi

ewed

,for

inte

rnal

circ

ulat

ion

only

ATLAS DRAFT

[GeV]T

p

4 6 8 10 12 14

2v

0

0.05

0.10-10 %

[GeV]T

p

4 6 8 10 12 14

2v

0

0.05

0.1 10-20 %

[GeV]T

p4 6 8 10 12 14

2v

0

0.05

0.1 30-40 %TAMU values for 20-40% [GeV]

Tp

4 6 8 10 12 14

2v

0

0.05

0.1 20-30 %TAMU values for 20-40%

[GeV]T

p4 6 8 10 12 14

2v

0

0.05

0.1 40-60 %

5 5.2 5.4 5.6 5.8 6

0

0.05

0.1

| < 2η|-1Pb+Pb, 0.14 nb

= 2.76 TeVNNs InternalATLAS

DATADABModTAMU

Figure 15: Comparison of the Pb+Pb heavy-flavor muon v2 with calculations from the TAMU and DABMod models.Each panel represents a di�erent centrality interval. For the 20–30% and 30–40% centrality intervals, the plottedTAMU values correspond to the 20–40% centrality interval. For the data, the error bars and shaded bands representstatistical and total uncertainties, respectively. For the model calculations, the bands represent theoretical systematicuncertainties.

25th April 2018 – 14:47 31

Not

revi

ewed

,for

inte

rnal

circ

ulat

ion

only

ATLAS DRAFT

[GeV]T

p

4 6 8 10 12 14

2v

0

0.05

0.10-10 %

[GeV]T

p

4 6 8 10 12 14

2v0

0.05

0.1 10-20 %

[GeV]T

p4 6 8 10 12 14

2v

0

0.05

0.1 30-40 %TAMU values for 20-40% [GeV]

Tp

4 6 8 10 12 14

2v

0

0.05

0.1 20-30 %TAMU values for 20-40%

[GeV]T

p4 6 8 10 12 14

2v

0

0.05

0.1 40-60 %

5 5.2 5.4 5.6 5.8 6

0

0.05

0.1

| < 2η|-1Pb+Pb, 0.14 nb

= 2.76 TeVNNs InternalATLAS

DATADABModTAMU

Figure 15: Comparison of the Pb+Pb heavy-flavor muon v2 with calculations from the TAMU and DABMod models.Each panel represents a di�erent centrality interval. For the 20–30% and 30–40% centrality intervals, the plottedTAMU values correspond to the 20–40% centrality interval. For the data, the error bars and shaded bands representstatistical and total uncertainties, respectively. For the model calculations, the bands represent theoretical systematicuncertainties.

25th April 2018 – 14:47 31

RAA 𝒗2

DABMod gives a good description of 𝒗2, underestimates RAA at low pT

TAMU describes the feature of data RAA, but slightly underestimates 𝒗2

Not

revi

ewed

,for

inte

rnal

circ

ulat

ion

only

ATLAS DRAFT

[GeV]T

p

4 6 8 10 12 14

2v

0

0.05

0.10-10 %

[GeV]T

p

4 6 8 10 12 14

2v

0

0.05

0.1 10-20 %

[GeV]T

p4 6 8 10 12 14

2v

0

0.05

0.1 30-40 %TAMU values for 20-40% [GeV]

Tp

4 6 8 10 12 14

2v

0

0.05

0.1 20-30 %TAMU values for 20-40%

[GeV]T

p4 6 8 10 12 14

2v

0

0.05

0.1 40-60 %

5 5.2 5.4 5.6 5.8 6

0

0.05

0.1

| < 2η|-1Pb+Pb, 0.14 nb

= 2.76 TeVNNs InternalATLAS

DATADABModTAMU

Figure 15: Comparison of the Pb+Pb heavy-flavor muon v2 with calculations from the TAMU and DABMod models.Each panel represents a di�erent centrality interval. For the 20–30% and 30–40% centrality intervals, the plottedTAMU values correspond to the 20–40% centrality interval. For the data, the error bars and shaded bands representstatistical and total uncertainties, respectively. For the model calculations, the bands represent theoretical systematicuncertainties.

25th April 2018 – 14:47 31

arXiv: 1805.05220 TAMU: PLB 735 (2014) 445

DABMod: PRC 96 (2017) 064903

Not

revi

ewed

,for

inte

rnal

circ

ulat

ion

only

ATLAS DRAFT

[GeV]T

p

4 6 8 10 12 14

2v

0

0.05

0.10-10 %

[GeV]T

p

4 6 8 10 12 14

2v

0

0.05

0.1 10-20 %

[GeV]T

p4 6 8 10 12 14

2v

0

0.05

0.1 30-40 %TAMU values for 20-40% [GeV]

Tp

4 6 8 10 12 14

2v

0

0.05

0.1 20-30 %TAMU values for 20-40%

[GeV]T

p4 6 8 10 12 14

2v

0

0.05

0.1 40-60 %

5 5.2 5.4 5.6 5.8 6

0

0.05

0.1

| < 2η|-1Pb+Pb, 0.14 nb

= 2.76 TeVNNs InternalATLAS

DATADABModTAMU

Figure 15: Comparison of the Pb+Pb heavy-flavor muon v2 with calculations from the TAMU and DABMod models.Each panel represents a di�erent centrality interval. For the 20–30% and 30–40% centrality intervals, the plottedTAMU values correspond to the 20–40% centrality interval. For the data, the error bars and shaded bands representstatistical and total uncertainties, respectively. For the model calculations, the bands represent theoretical systematicuncertainties.

25th April 2018 – 14:47 31

HF muon flow in Pb+Pb

�9

𝒗3 𝒗4

• Measured HF muon 𝒗3 and 𝒗4 agree with DABMod calculations within uncertainties

[GeV]T

p4 6 8 10 12 14

3v

0

0.01

0.02

0.03

0.04 0-10 %

[GeV]T

p4 6 8 10 12 14

3v

0

0.01

0.02

0.03

0.04 10-20 %

[GeV]Tp

4 6 8 10 12 14

3v

0.05−

0

0.05 30-40 %

5 5.2 5.4 5.6 5.8 6

0

0.01

0.02

0.03

0.04

| < 2η|

-1Pb+Pb, 0.14 nb = 2.76 TeVNNs

ATLAS

DATADABMod

[GeV]T

p4 6 8 10 12 14

3v

0

0.01

0.02

0.03

0.04 20-30 %

[GeV]Tp

4 6 8 10 12 14

3v

0.05−

0

0.05 40-60 %

[GeV]T

p4 6 8 10 12 14

4v

0.01−

0

0.01

0.02

0.03 0-10 %

[GeV]T

p4 6 8 10 12 14

4v

0.04−

0.02−

0

0.02

0.04 10-20 %

[GeV]Tp

4 6 8 10 12 14

4v

0

0.0530-40 % [GeV]

Tp

4 6 8 10 12 14

4v

0.04−

0.02−

0

0.02

0.04 20-30 %

[GeV]Tp

4 6 8 10 12 14

4v

0

0.0540-60 %

5 5.2 5.4 5.6 5.8 60.01−

0

0.01

0.02

0.03

| < 2η|-1Pb+Pb, 0.14 nb

= 2.76 TeVNNsATLAS

DATADABMod

arXiv: 1805.05220

J/𝝍 RAA in Pb+Pb

�10

[GeV]T

p10 15 20 25 30 35 40

AAR

0

0.2

0.4

0.6

0.8

1

1.2

1.4 ATLAS-1 = 5.02 TeV, 0.42 nbNNsPb+Pb,

-1 = 5.02 TeV, 25 pbs, pp, |y| < 2ψPrompt J/

Cent. 0-10 %Cent. 20-40 %Cent. 40-80 %

Correlated systematic uncer.

[GeV]T

p10 15 20 25 30 35 40

AAR

0

0.2

0.4

0.6

0.8

1

1.2

1.4 ATLAS-1 = 5.02 TeV, 0.42 nbNNsPb+Pb,

-1 = 5.02 TeV, 25 pbs, pp, |y| < 2ψNon-prompt J/

Cent. 0-10 %Cent. 20-40 %Cent. 40-80 %

Correlated systematic uncer.

Prompt J/𝝍 Non-Prompt J/𝝍

• J/𝝍 RAA measured in 3 centrality slices covers 0-80%• Prompt and non-prompt J/𝝍 show similar suppression (by accident?)• Non-prompt J/𝝍 RAA is flat in pT, similar to HF muon

J/𝝍 𝒗2 in Pb+Pb

�11 arXiv: 1805.04077 ATLAS-CONF-2018-013

10 15 20 25 30 [GeV]

Tp

0

0.05

0.1

0.15

0.2

2v ATLAS Preliminary-1 = 5.02 TeV, 0.42 nbNNsPb+Pb,

< 2, 0 - 60%y, ψPrompt J/

10 15 20 25 30 [GeV]

Tp

0

0.05

0.1

0.15

0.2

2v ATLAS Preliminary-1 = 5.02 TeV, 0.42 nbNNsPb+Pb,

< 2, 0 - 60%y, ψNon-prompt J/

Prompt J/𝝍 Non-Prompt J/𝝍

• J/𝝍 𝒗2 is extracted from event plane method

• Prompt and non-prompt contributions are separated by 2D yield fit• Significant prompt J/𝝍 𝒗2 at pT ~ 10 GeV

• Non-prompt J/𝝍 results favor a non-zero 𝒗2

J/𝝍 𝒗2 in Pb+Pb

�12

0 5 10 15 20 25 30 [GeV]

Tp

0.05

0.1

0.15

0.2

0.252v ATLAS Preliminary < 2, 0 - 60%y, 5.02 TeV, ψATLAS, Prompt J/

< 2, 0 - 60%y, 5.02 TeV, ψATLAS, Non-prompt J/, 5.02 TeV, 2.5 < y < 4, 20 - 40%ψALICE, Inclusive J/

< 2.4, 10 - 60%y, 2.76 TeV, 1.6 < ψCMS, Prompt J/ < 2.4, 10 - 60%y, 2.76 TeV, ψCMS, Prompt J/

• Prompt J/𝝍 𝒗2 measurements at the LHC show good agreement

• Full Run2 data could allow us to make more differential measurements

Non-prompt J/𝝍 vs. HF muon

�13

Not

revi

ewed

,for

inte

rnal

circ

ulat

ion

only

ATLAS DRAFT

[GeV] T

p4 6 8 10 12 14

AA

R

0

0.2

0.4

0.6

0.8

1 0-10%

[GeV] Tp

4 6 8 10 12 14

AA

R

0

0.2

0.4

0.6

0.8

1 20-30%TAMU values for 20-40%

[GeV] T

p4 6 8 10 12 14

AA

R0

0.2

0.4

0.6

0.8

1 10-20%

[GeV] Tp

4 6 8 10 12 14

AA

R

0

0.2

0.4

0.6

0.8

1 30-40%TAMU values for 20-40%

5 5.2 5.4 5.6 5.8 6

0

0.2

0.4

0.6

0.8

1

| < 1η|

-1, 570 nbpp-1Pb+Pb, 0.14 nb

= 2.76 TeVNNsATLAS

DataDABModTAMU

[GeV] Tp

4 6 8 10 12 14

AA

R

0

0.2

0.4

0.6

0.8

1 40-60%

Figure 14: Comparison of the measured heavy-flavor muon RAA in Pb+Pb collisions with the values predicted fromthe TAMU transport model and the DABMod model. Each panel represents a di�erent centrality interval. For the20–30% and 30–40% centrality intervals, the plotted TAMU values correspond to the 20–40% centrality interval.For the data, the error bars represent statistical uncertainties, the shaded bands represent the experimental systematicuncertainties, and the boxes indicate theoretical uncertainties from hTAAi. For the model calculations the bandsindicate the theoretical systematic uncertainties.

12th May 2018 – 19:26 30

Not

revi

ewed

,for

inte

rnal

circ

ulat

ion

only

ATLAS DRAFT

[GeV] T

p4 6 8 10 12 14

AA

R

0

0.2

0.4

0.6

0.8

1 0-10%

[GeV] Tp

4 6 8 10 12 14

AA

R

0

0.2

0.4

0.6

0.8

1 20-30%TAMU values for 20-40%

[GeV] T

p4 6 8 10 12 14

AA

R

0

0.2

0.4

0.6

0.8

1 10-20%

[GeV] Tp

4 6 8 10 12 14

AAR

0

0.2

0.4

0.6

0.8

1 30-40%TAMU values for 20-40%

5 5.2 5.4 5.6 5.8 6

0

0.2

0.4

0.6

0.8

1

| < 1η|

-1, 570 nbpp-1Pb+Pb, 0.14 nb

= 2.76 TeVNNsATLAS

DataDABModTAMU

[GeV] Tp

4 6 8 10 12 14

AA

R

0

0.2

0.4

0.6

0.8

1 40-60%

Figure 14: Comparison of the measured heavy-flavor muon RAA in Pb+Pb collisions with the values predicted fromthe TAMU transport model and the DABMod model. Each panel represents a di�erent centrality interval. For the20–30% and 30–40% centrality intervals, the plotted TAMU values correspond to the 20–40% centrality interval.For the data, the error bars represent statistical uncertainties, the shaded bands represent the experimental systematicuncertainties, and the boxes indicate theoretical uncertainties from hTAAi. For the model calculations the bandsindicate the theoretical systematic uncertainties.

12th May 2018 – 19:26 30

Non-prompt J/𝝍0-10%

HF muon0-10%

Not

revi

ewed

,for

inte

rnal

circ

ulat

ion

only

ATLAS DRAFT

[GeV]T

p

4 6 8 10 12 14

2v

0

0.05

0.10-10 %

[GeV]T

p

4 6 8 10 12 14

2v

0

0.05

0.1 10-20 %

[GeV]T

p4 6 8 10 12 14

2v

0

0.05

0.1 30-40 %TAMU values for 20-40% [GeV]

Tp

4 6 8 10 12 14

2v

0

0.05

0.1 20-30 %TAMU values for 20-40%

[GeV]T

p4 6 8 10 12 14

2v

0

0.05

0.1 40-60 %

5 5.2 5.4 5.6 5.8 6

0

0.05

0.1

| < 2η|-1Pb+Pb, 0.14 nb

= 2.76 TeVNNs InternalATLAS

DATADABModTAMU

Figure 15: Comparison of the Pb+Pb heavy-flavor muon v2 with calculations from the TAMU and DABMod models.Each panel represents a di�erent centrality interval. For the 20–30% and 30–40% centrality intervals, the plottedTAMU values correspond to the 20–40% centrality interval. For the data, the error bars and shaded bands representstatistical and total uncertainties, respectively. For the model calculations, the bands represent theoretical systematicuncertainties.

25th April 2018 – 14:47 31

Not

revi

ewed

,for

inte

rnal

circ

ulat

ion

only

ATLAS DRAFT

[GeV]T

p4 6 8 10 12 14

2v

0

0.05

0.1Event planeScalar product

0-10 %

[GeV]T

p4 6 8 10 12 14

2v

0

0.05

0.1 10-20 %

[GeV]Tp

4 6 8 10 12 14

2v

0

0.05

0.1 30-40 %

5 5.2 5.4 5.6 5.8 6

0

0.05

0.1

| < 2η|-1Pb+Pb, 0.14 nb

= 2.76 TeVNNsATLAS

[GeV]T

p4 6 8 10 12 14

2v

0

0.05

0.1 20-30 %

[GeV]Tp

4 6 8 10 12 14

2v

0

0.05

0.1 40-60 %

Figure 10: The pT dependence of the Pb+Pb heavy-flavor muon v2. Results are shown for both the EP and SPmethods. Each panel represents a di�erent centrality interval. The error bars and shaded bands represent statisticaland total uncertainties, respectively, and are shown only for the EP v2. The horizontal dashed lines indicate v2 = 0.

12th May 2018 – 19:26 25

Not

revi

ewed

,for

inte

rnal

circ

ulat

ion

only

ATLAS DRAFT

[GeV]T

p4 6 8 10 12 14

2v

0

0.05

0.1Event planeScalar product

0-10 %

[GeV]T

p4 6 8 10 12 14

2v

0

0.05

0.1 10-20 %

[GeV]Tp

4 6 8 10 12 14

2v

0

0.05

0.1 30-40 %

5 5.2 5.4 5.6 5.8 6

0

0.05

0.1

| < 2η|-1Pb+Pb, 0.14 nb

= 2.76 TeVNNsATLAS

[GeV]T

p4 6 8 10 12 14

2v

0

0.05

0.1 20-30 %

[GeV]Tp

4 6 8 10 12 14

2v

0

0.05

0.1 40-60 %

Figure 10: The pT dependence of the Pb+Pb heavy-flavor muon v2. Results are shown for both the EP and SPmethods. Each panel represents a di�erent centrality interval. The error bars and shaded bands represent statisticaland total uncertainties, respectively, and are shown only for the EP v2. The horizontal dashed lines indicate v2 = 0.

12th May 2018 – 19:26 25

Non-prompt J/𝝍0-60%

HF muon10 -20%

RAA 𝒗2

arXiv: 1805.05220 arXiv: 1805.04077 ATLAS-CONF-2018-013

Similar suppression and flow for HF muon and non-prompt J/𝝍

Summary

�14

Large system: • Significant suppression • Significant flow

Summary

�15

‣ How much do cold nuclear matter (CNM) effects (nPDF, final state interaction etc.) contribute?

Large system: • Significant suppression • Significant flow

Summary

�16

‣ How much do cold nuclear matter (CNM) effects (nPDF, final state interaction etc.) contribute?

‣ Is p+Pb a pure system to study CNM effects? “Can we turn the QGP off ?” (quote from Micheal)

Large system: • Significant suppression • Significant flow

Quarkonia production in p+Pb

�17

[GeV]Tp

0 5 10 15 20 25 30 35 40

pPb

R

0.2

0.4

0.6

0.8

1

1.2

1.4

1.6

1.8ATLAS

= 5.02 TeVs, ppATLAS < 2.0y, -1L = 25 pb

= 5.02 TeVNNs+Pb, pATLAS , -2.0 < y* < 1.5-1L = 28 nb

(1S)ϒ

ψPrompt J/

= 5.02 TeVNNs+Pb, pALICE , -1.37 < y* < 0.43-1bµL = 51 ψInclusive J/

[GeV]Tp

0 5 10 15 20 25 30 35 40

pPb

R

0.2

0.4

0.6

0.8

1

1.2

1.4

1.6

1.8

EPS09 NLOR. Nelson et al., Phys. Rev. C 87(2013)014908

• ϒ(1S) is found to be suppressed at low pT

• Difference between ϒ(1S) and J/𝜓 RpPb at 10 < pT < 15 GeV• Results could be described by nPDF / co-mover effects

Quarkonia production in p+Pb

�18

[GeV]Tp

0 5 10 15 20 25 30 35 40

pPb

R

0.2

0.4

0.6

0.8

1

1.2

1.4

1.6

1.8ATLAS

= 5.02 TeVs, ppATLAS < 2.0y, -1L = 25 pb

= 5.02 TeVNNs+Pb, pATLAS , -2.0 < y* < 1.5-1L = 28 nb

(1S)ϒ

ψPrompt J/

= 5.02 TeVNNs+Pb, pALICE , -1.37 < y* < 0.43-1bµL = 51 ψInclusive J/

*y4− 2− 0 2 4

pPb

R

0

0.5

1

1.5

2ATLAS

= 5.02 TeVNNs+Pb, p(1S)ϒ

< 40 GeVTpATLAS,

< 15 GeVTpLHCb,

> 0 GeVTpALICE,

Comover InteractionE. Ferreiro, Phys. Lett. B 731(2014)57

[GeV]Tp

0 5 10 15 20 25 30 35 40

pPb

R

0.2

0.4

0.6

0.8

1

1.2

1.4

1.6

1.8

EPS09 NLOR. Nelson et al., Phys. Rev. C 87(2013)014908

[GeV]Tp

0 5 10 15 20 25 30 35 40

pPb

R

0.2

0.4

0.6

0.8

1

1.2

1.4

1.6

1.8

EPS09 NLOR. Nelson et al., Phys. Rev. C 87(2013)014908

• ϒ(1S) is found to be suppressed at low pT

• Difference between ϒ(1S) and J/𝜓 RpPb at 10 < pT < 15 GeV• Results could be described by nPDF / co-mover effects

Excited quarkonia production

�19

Centrality [%]40 - 90 20 - 40 5 - 20 0 - 5

ψ(2

S) /

J/ψ

Pbpρ

0

0.5

1

1.5

2-1 = 5.02 TeV, L = 28 nbNNs+Pb, p

-1 = 5.02 TeV, L = 25 pbs, pp

ATLAS < 40 GeV

Tp8 <

* < 1.5y-2.0 <

Double Ratioψ(2S) to J/ψPrompt

• Stronger suppression of excited quarkonia states in p+Pb wrt. ground states

• Slight centrality dependence

Charmonia

Excited quarkonia production

�20

Centrality [%]40 - 90 20 - 40 5 - 20 0 - 5

ψ(2

S) /

J/ψ

Pbpρ

0

0.5

1

1.5

2-1 = 5.02 TeV, L = 28 nbNNs+Pb, p

-1 = 5.02 TeV, L = 25 pbs, pp

ATLAS < 40 GeV

Tp8 <

* < 1.5y-2.0 <

Double Ratioψ(2S) to J/ψPrompt

(1S)ϒ(2S)/ϒ (1S)ϒ(3S)/ϒ

(1S)

ϒ(n

S) /

ϒ

Pbpρ

0

0.2

0.4

0.6

0.8

1

1.2ATLAS

< 40 TeVTpATLAS,

< 2.0y, -1 = 5.02 TeV, L = 25 pbs, pp* < 1.5y, -2.0 < -1 = 5.02 TeV, L = 28 nbNNs+Pb, p

CMS < 1.93y, -1 = 2.76 TeV, L = 5.4 pbs, pp

< 1.93*y, -1 = 5.02 TeV, L = 31 nbNNs+Pb, p

Comover Interaction (E.G. Ferreiro, Phys. Lett. B 731(2014)57)

• Stronger suppression of excited quarkonia states in p+Pb wrt. ground states

• Slight centrality dependence• Can be described by final state interactions

Charmonia Bottomonia

Prompt D0 yield in p+Pb

�21 ATLAS-CONF-2017-073

• Prompt D0 with 3 < pT < 30 GeV• FONLL prediction 8 TeV pp collisions scaled by APb

• Data and scaled FONLL agree within uncertainties, small modifications of D meson production

[GeV]T

p5 10 15 20 25 30

b/G

eV]

µ) [yd Tp

/ (d

σ2 d

5−10

4−10

3−10

2−10

1−101

10

210

310

410

510 Preliminary ATLAS-1bµ = 8.16 TeV, 76.3 NNs+Pb p

Production0DPrompt

FONLL* < 0.5y, 0.0 < 3 10×Data * < 0.0y, -0.5 < 2 10×Data * < -0.5y, -1.0 < 1 10×Data * < -1.0y, -1.5 < 0 10×Data

5 10 15 20 25 30

Dat

a / T

heor

y

123 * < 0.0y-0.5 <

* < 0.5y0.0 <

5 10 15 20 25 30

Dat

a / T

heor

y

123 * < -0.5y-1.0 <

[GeV]Tp

5 10 15 20 25 30

Dat

a / T

heor

y

0123 * < -1.0y-1.5 <

[GeV]T

p5 10 15 20 25 30

b/G

eV]

µ) [yd Tp

/ (d

σ2 d

5−10

4−10

3−10

2−10

1−101

10

210

310

410

510 Preliminary ATLAS-1bµ = 8.16 TeV, 76.3 NNs+Pb p

Production0DPrompt

FONLL* < 0.5y, 0.0 < 3 10×Data * < 0.0y, -0.5 < 2 10×Data * < -0.5y, -1.0 < 1 10×Data * < -1.0y, -1.5 < 0 10×Data

5 10 15 20 25 30

Dat

a / T

heor

y

123 * < 0.0y-0.5 <

* < 0.5y0.0 <

5 10 15 20 25 30

Dat

a / T

heor

y

123 * < -0.5y-1.0 <

[GeV]Tp

5 10 15 20 25 30

Dat

a / T

heor

y

0123 * < -1.0y-1.5 <

[GeV]T

p5 10 15 20 25 30

b/G

eV]

µ) [yd Tp

/ (d

σ2 d

5−10

4−10

3−10

2−10

1−101

10

210

310

410

510 Preliminary ATLAS-1bµ = 8.16 TeV, 76.3 NNs+Pb p

Production0DPrompt

FONLL* < 0.5y, 0.0 < 3 10×Data * < 0.0y, -0.5 < 2 10×Data * < -0.5y, -1.0 < 1 10×Data * < -1.0y, -1.5 < 0 10×Data

5 10 15 20 25 30

Dat

a / T

heor

y

123 * < 0.0y-0.5 <

* < 0.5y0.0 <

5 10 15 20 25 30

Dat

a / T

heor

y

123 * < -0.5y-1.0 <

[GeV]Tp

5 10 15 20 25 30

Dat

a / T

heor

y

0123 * < -1.0y-1.5 <

Prompt D0 yield in p+Pb

�22

[GeV]Tp

5 10 15 20 25 30

FBR

00.20.40.60.8

11.21.41.61.8

2ATLAS Preliminary

= 8.16 TeVNNs+Pb pπK→0DPrompt

Forward, 0 < y* < 0.5Backward, -0.5 < y* < 0

ATLAS-CONF-2017-073

Prompt D0

• Forward: 0 < y* < 0.5, Backward: -0.5 < y* < 0• No obvious modification in forward wrt. backward for prompt D

[GeV]Tp

5 10 15 20 25 30

FBR

00.20.40.60.8

11.21.41.61.8

2ATLAS Preliminary

= 8.16 TeVNNs+Pb p

sππK→sπ0D→*DPrompt

Forward, 0 < y* < 0.5Backward, -0.5 < y* < 0

Prompt D*

HF muon v2 in p+Pb

�23

• Event collected by muon + high multiplicity trigger• Two particle correlation between HF muon and charged hadrons• Non-flow subtracted with template fit and fake HF muon background

suppressed with momentum imbalance cut

Low Multiplicity High Multiplicity

ATLAS-CONF-2017-006

φ∆

0

2

4

η∆

-4-2

02

4

)φ

∆,η

∆C

(

0.9

1

1.1

1.2

ATLAS Preliminary +Pbp-1=8.16 TeV, 171 nbNNs

Correlationsµ-h

<5 GeVh

Tp0.5<

<4.5 GeVµ

Tp4<

<40 rec chN≤0

φ∆

0

2

4

η∆

-4-2

02

4

)φ

∆,η

∆C

(0.98

1

1.02

1.04

ATLAS Preliminary +Pbp-1=8.16 TeV, 171 nbNNs

Correlationsµ-h

<5 GeVh

Tp0.5<

<4.5 GeVµ

Tp4<

200≥ rec chN

φ∆

-1 0 1 2 3 4

)φ

∆C

(

0.98

1

1.02

1.04

1.06

Correlations)µ-h) (φ∆C(

) + Gφ∆(periph

FC

)φ∆(templ

C

(0)periph

) +FCφ∆(ridge

C

(0)periph

G + FC

ATLAS Preliminary-1=8.16 TeV, 171 nbNNs+Pb p

<6 GeVµ

Tp<5 GeV, 4<a

Tp0.5<

<160 rec

chN≤|<5, 140η∆1<|

φ∆

-1 0 1 2 3 4

)φ

∆C

(

0.97

0.98

0.99

1

1.01

1.02

1.03

1.04

1.05

Correlations)µ-h) (φ∆C(

) + Gφ∆(periph

FC

)φ∆(templ

C

(0)periph

) +FCφ∆(ridge

C

(0)periph

G + FC

ATLAS Preliminary-1=8.16 TeV, 171 nbNNs+Pb p

<6 GeVµ

Tp<5 GeV, 4<a

Tp0.5<

200≥ rec

chN|<5, η∆1<|

HF muon v2 in p+Pb

�24

recchN

0 50 100 150 200 250 300

)b T

p(2v

0

0.05

0.1

Correlations)h-h<5 GeV (b

Tp0.5<

Correlations)µ-h<6 GeV (µ

Tp4<

ATLAS Preliminary-1=8.16 TeV, 171 nbNNs+Pb p

<5 GeVa

Tp0.5<

• No obvious multiplicity dependence of HF muon 𝒗2, no obvious bias by combining events collected by different high multiplicity triggers

• 𝒗2 (HF muon) ~ 0.6 × 𝒗2 (hadron)

ATLAS-CONF-2017-006

[GeV]µ

Tp

4 4.5 5 5.5 6 6.5 7 7.5 8

)µ T

p(2

v

0

0.02

0.04

0.06

0.08ATLAS Preliminary

-1=8.16 TeV, 171 nbNNs+Pb p

<5 GeVa

Tp0.5<

100≥ rec chN

Correlationsµ-h

D*—hadron correlation in p+Pb

�25

• Template fit performed to background subtracted correlation function, result favors 𝒗2,2 > 0

• 1~2σ hint for cos(2Δ𝜙) azimuthal modulation of D*

φΔ

0 0.5 1 1.5 2 2.5 3

)φΔ (

C

0.94

0.96

0.98

1

1.02

1.04

1.06

1.08

1.1

1.12ATLAS Preliminary

= 8.16 TeVNNs+Pb p < 160chN120 <

Simultaneous fit

)φΔ (C)φΔ( periphCf + g

)φΔ( templC(0) periphCf) + φΔ( ridgeC

0.4(syst)± 3.3(stat) ± = 8.5 310×2,2v

φΔ

0 0.5 1 1.5 2 2.5 3

)φΔ (

C

0.96

0.98

1

1.02

1.04

1.06 ATLAS Preliminary = 8.16 TeVNNs+Pb p

< 240chN160 < Simultaneous fit

)φΔ (C)φΔ( periphCf + g

)φΔ( templC(0) periphCf) + φΔ( ridgeC

0.3(syst)± 2.6(stat) ± = 3.1 310×2,2v

ATLAS-CONF-2017-073

Summary

�26

Large system: • Significant suppression • Significant flow

Small system: • No evidence of suppression due to QGP interaction, modifications

could be described by CNM effects • Significant HF muon azimuthal anisotropy is observed, collective

flow?

Backup slides

Heavy flavor probes in ATLAS

�28

Momentum imbalance Δ𝑝 = 𝑝ID - ( 𝑝MS + 𝑝Calo )template fit

HF muons

Not

revi

ewed

,for

inte

rnal

circ

ulat

ion

only

ATLAS DRAFT

IDp/pΔ

0.4− 0.3− 0.2− 0.1− 0 0.1 0.2 0.3 0.4 0.5 0.6

IDp/p

Δdµ

Nd

µN1

0

2

4

6

8 SimulationATLAS = 2.76 TeVNNs

< 6 GeVTp5 <

| < 1η|

Pb+Pb 0-60%signalbackground

ppsignalbackground

Figure 1: Signal and background template distributions in pp collisions (square points) and Pb+Pb collisions (circularpoints) in the 0–60% centrality interval for muons having 5 < pT < 6 GeV and |⌘ | < 1. The signal and backgrounddistributions are separately normalized such that their integral is unity. For clarity, the background distribution isbinned more coarsely.

in the calorimeter have, on average, smaller pT compared to the parent particle. As a result, background204

muons are expected to have �p/pID > 0.205

Distributions for �p/pID are obtained from the simulated samples separately for signal muons and for206

background muons. The signal muons include muons directly produced in electromagnetic decays of207

hadrons, in decays of ⌧-leptons, in decays of W and Z bosons, in decays of top quarks, and in semi-208

leptonic decays of heavy-flavor hadrons; this last contribution dominates the signal sample, contributing209

about 99% of the muons over the pT range measured in this analysis (Ref. [55] and references therein). The210

di�erent contributions to the background – pion decays in flight, kaon decays in flight, muons produced211

by secondary interactions of prompt particles, and mis-associations – are evaluated separately. Figure 1212

shows MC distributions of �p/pID for signal and background muons having 5 < pT < 6 GeV for Pb+Pb213

collisions in the centrality range 0–60% and for pp collisions. The �p/pID distribution for signal muons214

is centered at zero while the distribution for background muons is shifted to positive values. The signal215

distributions show only modest di�erences between pp and Pb+Pb collisions. Similarly, when making216

separate templates for di�erent Pb+Pb collision centralities, a weak dependence of the signal templates217

on centrality is observed. The background �p/pID distributions are much broader and are insensitive to218

the centrality-dependent e�ects seen in the signal distributions.219

A template-fitting procedure is used to estimate statistically the signal fraction for each kinematic and220

centrality selection used in the analysis. The measured �p/pID distribution is assumed to result from a221

12th May 2018 – 19:26 9

𝑝ID 𝑝Calo

Signal 𝜇

𝜇

Punch-throughLate K/𝜋 decay

𝑝MS

Inner Detector Calorimeter Muon Spectrometer

arXiv: 1805.05220HF muon

pT [GeV]14 4010420

Heavy flavor probes in ATLAS

�29

I Ca M

[GeV]µµm2.6 2.8 3 3.2 3.4 3.6 3.8 4 4.2

Wei

ghte

d Ev

ents

/ 20

MeV

210

310

410

510 ATLAS-1 = 5.02 TeV, L = 28 nbNNs+Pb p

< 11.0 GeVµµ

Tp10.0 < * < 1.5y-2.0 <

/ndof = 1.342χ DataFit

(nS)ψPrompt (nS)ψNon-prompt

Bkg

[ps]µµτ

4− 2− 0 2 4 6 8 10

Wei

ghte

d Ev

ents

/ 0.

1 ps

1

10

210

310

410

510 ATLAS-1 = 5.02 TeV, L = 28 nbNNs+Pb p

< 11.0 GeVµµ

Tp10.0 < * < 1.5y-2.0 <

/ndof = 1.662χ DataFit

(nS)ψPrompt (nS)ψNon-prompt

Bkg

• Non-Prompt J/𝝍 refers to J/𝝍 from B hadron decays• Measured via dimuon decay channel• Extracted from 2D fit in invariant mass and pseudo

proper lifetime

Fit projection on mass Fit projection on pseudo proper lifetime

EPJC 78 (2018) 171

Heavy flavor probes in ATLAS

�30

) [GeV]π(Km1.75 1.8 1.85 1.9 1.95 2

Wei

ghte

d Ev

ents

/ 5

MeV

0

1

2

3

4

5

6

7

8310×

ATLAS Preliminary = 8.16 TeVNNs+Pb p

DataFit

Signal0D+0DComb. Bkg

Swappedπk

* < 0.5y-0.5 < < 6.0 GeV

Tp5.0 <

/ndof = 0.92χ

D0→K𝜋

D0→K𝜋 Decay

ATLAS-CONF-2017-073

x

y

decay lengthLxy

point

ing an

gle α xy

PV

SV

𝜋

K

refitted K𝜋 vertex pT

D0

• Event collected by MinBias trigger

• Dominated by open charm contribution

D0→K𝜋 decay vertex selection: Vertex probabilityPointing angleDecay length significance Lxy / σ(Lxy)

) [MeV]π(Km) - ππ(Km140 145 150 155 160 165 170

Wei

ghte

d Ev

ents

/ 0.

5 M

eV

0

0.5

1

1.5

2

2.5

3

3.5

310×

ATLAS Preliminary = 8.16 TeVNNs+Pb p

DataFit

* Signal)±(DBkg

* < 0.5y-0.5 < < 6.0 GeV

Tp5.0 <

) < 1.96 GeVπK(m1.75 <

/ndof = 1.22χ

D*→D0𝜋→K𝜋𝜋

HF muon 𝒗2 in Pb+Pb

�31

• Good agreement between ATLAS and ALICE

• Smaller uncertainties of ATLAS results would provide tight constraints on models

[GeV]T

p4 6 8 10 12 14

2v

0

0.1

0.2

0.30-10 %

[GeV]T

p4 6 8 10 12 14

2v

0

0.1

0.2

0.310-20 %

[GeV]Tp

4 6 8 10 12 14 2v

0

0.1

0.2

0.3 30-40 %ALICE values for 20-40%

5 5.2 5.4 5.6 5.8 6

0

0.1

0.2

0.3

| < 2η|-1Pb+Pb, 0.14 nb

= 2.76 TeVNNsATLAS

|<2η |±µATLAS HF |y|<0.7±eALICE HF 2.5<y<4±µALICE HF

[GeV]T

p4 6 8 10 12 14

2v

0

0.1

0.2

0.320-30 %ALICE values for 20-40%

[GeV]Tp

4 6 8 10 12 14

2v

0

0.1

0.2

0.3 40-60 %

arXiv: 1805.05220

HF muon in pp vs FONLL

�32

Not

revi

ewed

,for

inte

rnal

circ

ulat

ion

only

ATLAS DRAFT

[GeV]T

p

4 6 8 10 12 14

] -1

GeV

× [n

b ηd Tpdµ

HF

σ2 d

1

10

210

310ATLAS

= 2.76 TeVs -1, 570 nbpp

| < 1η|

, 2.76 TeVppFONLL(CTEQ6.6)FONLL(CTEQ6.6) bottom onlyFONLL(CTEQ6.6) charm only

[GeV]T

p4 6 8 10 12 14

FONL

Lσ/

ppσ 1

2

[GeV]Tp

4 6 8 10 12 14

cσ/ b

σ

0

1From FONLL

Figure 5: Top panel: the pT dependence of the measured heavy-flavor muon cross-section inp

s = 2.76 TeV ppcollisions. The data points are plotted at the average muon pT within a given pT interval. The vertical bars and bandson the data points indicate statistical and systematic uncertainties, respectively. The cross-section for heavy-flavordecays from FONLL calculations is also shown, along with the individual contributions from bottom and charmquarks. For the FONLL calculations, the vertical width of the band represents theoretical systematic uncertainties.Middle panel: the ratio of the measured and FONLL cross-sections integrated over each pT interval. Statistical andsystematic uncertainties in the data are indicated by error bars and gray shaded boxes, respectively. The systematicuncertainty of the ratio from FONLL is indicated by the shaded band centered on unity. Bottom panel: the ratio ofthe bottom contribution to the charm contribution in the FONLL calculations. All results are averaged over |⌘ | < 1.

12th May 2018 – 19:26 20

arXiv: 1805.05220

HF muon yields wrt. event plane

�33

|2Ψ - φ2|0 0.5 1 1.5 2 2.5 3

|) 2Ψ-φ

d(2|

µNd

µ

N1

0.9

1

1.1

10-20%

| < 2η| < 4.5 GeV

Tp4 <

= 2.76 TeVNNs-1Pb+Pb, 0.14 nb

ATLAS

|2Ψ - φ2|0 0.5 1 1.5 2 2.5 3

|) 2Ψ-φ

d(2|

µNd

µ

N1

0.9

1

1.1

40-60%

| < 2η| < 4.5 GeV

Tp4 <

= 2.76 TeVNNs-1Pb+Pb, 0.14 nb

ATLAS

|2Ψ - φ2|0 0.5 1 1.5 2 2.5 3

|) 2Ψ-φ

d(2|

µNd

µ

N1 1

1.1

10-20%

| < 2η| < 10 GeV

Tp8 <

= 2.76 TeVNNs-1Pb+Pb, 0.14 nb

ATLAS

|2Ψ - φ2|0 0.5 1 1.5 2 2.5 3

|) 2Ψ-φ

d(2|

µNd

µ

N10.9

1

1.1

40-60%

| < 2η| < 10 GeV

Tp8 <

= 2.76 TeVNNs-1Pb+Pb, 0.14 nb

ATLAS

arXiv: 1805.05220

J/𝝍 RAA in Pb+Pb

�34

partN0 50 100 150 200 250 300 350 400

AAR

0.2

0.4

0.6

0.8

1

1.2

1.4 ATLAS-1 = 5.02 TeV, 0.42 nbNNsPb+Pb,

-1 = 5.02 TeV, 25 pbs, pp, |y| < 2ψPrompt J/

< 40 GeVT

9 < pCorrelated systematic uncer.

0-5%5-10%

10-20%

20-30%

30-40%

40-50%50-60%

60-80%

partN0 50 100 150 200 250 300 350 400

AAR

0.2

0.4

0.6

0.8

1

1.2

1.4 ATLAS-1 = 5.02 TeV, 0.42 nbNNsPb+Pb,

-1 = 5.02 TeV, 25 pbs, pp, |y| < 2ψNon-Prompt J/

< 40 GeVT

9 < pCorrelated systematic uncer.

0-5%5-10%10-20%

20-30%

30-40%40-50%

50-60%

60-80%

D reconstruction

�35

Trigger selection MinBias (+HMT for correlation)

D0 selection

Two tracks, pT > 1GeV 𝜋 and K masses assigned in turn, 1.7 < m(K𝜋) < 2.0 GeV

Vertex probability Pointing angle cos𝛼xy Lxy / σ(Lxy)

D* selection

A selected D0 vertex An additional track (𝜋 mass), same charge with the 𝜋 in D0 with soft pion pT > 400 MeV (for yield) or 250 MeV (for correlation) 1.75 < m(K𝜋) < 1.96 GeV as SR m(K𝜋) < 1.76 GeV or m(K𝜋) > 1.96 GeV as CR

Cross sections

�36

*y1.5− 1− 0.5− 0 0.5

b]µ [ ydσd

0

20

40

60

80

100

120

140 ATLAS Preliminary-1bµ = 8.16 TeV, 76.3 NNs+Pb p

Production0DPrompt < 30 GeV

Tp3 <

FONLL

*y1.5− 1− 0.5− 0 0.5

b]µ [ ydσd

2

4

6

8

10

12

14 ATLAS Preliminary-1bµ = 8.16 TeV, 76.3 NNs+Pb p

* ProductionDPrompt < 30 GeV

Tp5 <

FONLL

‣ Prompt D0 (3 < pT < 30 GeV) and prompt D* (5 < pT < 30 GeV)‣ |ylab(D)| < 1.0 for better mass resolution ⟶ -1.5 < y* < 0.5‣ FONLL (fixed-order next-leading-logarithm) prediction

extrapolated from 7 and 8 TeV calculates, and scaled by 208‣ Relatively small modification in p+Pb

ATLAS-CONF-2017-073

Prompt D*Prompt D0 FONLL uncertainties‣ renormalisation scale‣ factorization scale ‣ charm quark mass‣ pdf

Recoil jet bias estimation

�37

Not

revi

ewed

,for

inte

rnal

circ

ulat

ion

only

ATLAS DRAFT

[GeV]T

p4 6 8 10 12 14

2v

0

0.05

0.12vEvent plane

Jet-bias

0-10 %

[GeV]T

p4 6 8 10 12 14

2v

0

0.05

0.1 10-20 %

[GeV]Tp

4 6 8 10 12 14

2v

0

0.05

0.1 30-40 %

5 5.2 5.4 5.6 5.8 6

0

0.05

0.1

| < 2η|-1Pb+Pb, 0.14 nb

= 2.76 TeVNNsATLAS

[GeV]T

p4 6 8 10 12 14

2v

0

0.05

0.1 20-30 %

[GeV]Tp

4 6 8 10 12 14

2v

0

0.05

0.1 40-60 %

Figure 18: The pT dependence of the Pb+Pb heavy-flavor muon v2. Each panel represents a di�erent centralityinterval. The error bars and shaded bands represent statistical and total uncertainties, respectively. The horizontalcontinuous lines indicate the jet bias on the v2 estimated using the data-overlay MC. Half of the estimated bias isused as a correction to the measured v2 and half is used as a systematic uncertainty (Section 3.5). The horizontaldashed lines indicate v2 = 0.

12th May 2018 – 19:26 42

Template fits — 𝛍-h in p+Pb

�38

φ∆

-1 0 1 2 3 4

)φ

∆C

(

0.96

0.98

1

1.02

1.04

1.06

1.08

1.1

1.12

Correlations)µ-h) (φ∆C(

) + Gφ∆(periph

FC

)φ∆(templ

C

(0)periph

) +FCφ∆(ridge

C

(0)periph

G + FC

ATLAS Preliminary-1=8.16 TeV, 171 nbNNs+Pb p

<6 GeVµ

Tp<5 GeV, 4<a

Tp0.5<

<80 rec

chN≤|<5, 60η∆1<|

φ∆

-1 0 1 2 3 4

)φ

∆C

(

0.96

0.98

1

1.02

1.04

1.06

1.08

Correlations)µ-h) (φ∆C(

) + Gφ∆(periph

FC

)φ∆(templ

C

(0)periph

) +FCφ∆(ridge

C

(0)periph

G + FC

ATLAS Preliminary-1=8.16 TeV, 171 nbNNs+Pb p

<6 GeVµ

Tp<5 GeV, 4<a

Tp0.5<

<120 rec

chN≤|<5, 100η∆1<|

φ∆

-1 0 1 2 3 4

)φ

∆C

(

0.98

1

1.02

1.04

1.06

Correlations)µ-h) (φ∆C(

) + Gφ∆(periph

FC

)φ∆(templ

C

(0)periph

) +FCφ∆(ridge

C

(0)periph

G + FC

ATLAS Preliminary-1=8.16 TeV, 171 nbNNs+Pb p

<6 GeVµ

Tp<5 GeV, 4<a

Tp0.5<

<160 rec

chN≤|<5, 140η∆1<|

φ∆

-1 0 1 2 3 4

)φ

∆C

(

0.97

0.98

0.99

1

1.01

1.02

1.03

1.04

1.05

Correlations)µ-h) (φ∆C(

) + Gφ∆(periph

FC

)φ∆(templ

C

(0)periph

) +FCφ∆(ridge

C

(0)periph

G + FC

ATLAS Preliminary-1=8.16 TeV, 171 nbNNs+Pb p

<6 GeVµ

Tp<5 GeV, 4<a

Tp0.5<

200≥ rec

chN|<5, η∆1<|

Systematics for mu-h correlation

�39

• Choice of peripheral bin 0-20, 20-40, 10-20, 20-30, 30-40

• Background muons • Efficiency correction • Track/muon selection • Pileup • Acceptance

systematics for muon-hadron correlation

|η∆|1 1.5 2 2.5 3 3.5 4 4.5 5

2,2

v

0

0.005

0.01

0.015

Template Fits

Fourier Transform

ATLAS Preliminary-1=8.16 TeV, 171 nbNNs+Pb p

Correlationsµ-h

100≥ rec chN

<5 GeVh

Tp0.5<

<6 GeVµ

Tp4<

ATLAS Preliminary-1=8.16 TeV, 171 nbNNs+Pb p

Correlationsµ-h

100≥ rec chN

<5 GeVh

Tp0.5<

<6 GeVµ

Tp4<

|η∆|1 1.5 2 2.5 3 3.5 4 4.5 5

2v

0.05

0.1

0.15

0.2

Template Fits

Fourier Transform

ATLAS Preliminary-1=8.16 TeV, 171 nbNNs+Pb p

Correlationsµ-h

100≥ rec chN

<5 GeVh

Tp0.5<

<6 GeVµ

Tp4<

D*-hadron correlation

�40

• Event collected by MinBias and high multiplicity triggers • Third soft pion from D* decay with pT > 250 MeV for larger D* fiducial

volume ( 3 < pT < 30 GeV and -1.5 < y* < 0.5 ) • Charged particles pT > 0.25 GeV, Δ𝜂 > 1 for more statistics for 2PC • Using sideband region 150 < Δm < 170 MeV to estimation the background

correlation function

) [MeV]π(Km) - ππ(Km140 145 150 155 160 165 170

Wei

ghte

d Ev

ents

/ 0.

5 M

eV

0

0.5

1

1.5

2

2.5

3

3.5

310×

ATLAS Preliminary = 8.16 TeVNNs+Pb p

DataFit

* Signal)±(DBkg

* < 0.5y-0.5 < < 6.0 GeV

Tp5.0 <

) < 1.96 GeVπK(m1.75 <

/ndof = 1.22χ

φΔ

0 0.5 1 1.5 2 2.5 3

)φΔ(

C

0.9

0.95

1

1.05

1.1

1.15ATLAS Preliminary

= 8.16 TeVNNs+Pb pSide BandsBefore subtractionAfter subtraction

> 250 MeVtrkTp

< 80chN10 < > 1.0ηΔ

h-h correlation

D*-h correlation

ATLAS-CONF-2017-073Signal Region

Sideband Region

φΔ

0 0.5 1 1.5 2 2.5 3

)φΔ (

C0.95

1

1.05

1.1

1.15ATLAS Preliminary

= 8.16 TeVNNs+Pb p* < 0.5y-1.5 <

< 30 GeVTp3 <

Correlationh*-D < 80chN10 <

)φΔ (CFit

> 1.0ηΔ > 250 MeV, trkTp

D*-h correlation fit

�41

• Assuming weak multiplicity dependence of near-side D*-h long range correlation

• Peripheral reference bin 10 < Nch < 80

• Simultaneous template fit method applied to low and high Nch single correlation functions φΔ

0 0.5 1 1.5 2 2.5 3

)φΔ (

C

0.920.940.960.98

11.021.041.061.08

1.11.12

ATLAS Preliminary = 8.16 TeVNNs+Pb p

< 120chN80 < Simultaneous fit

)φΔ (C)φΔ( periphCf + g

)φΔ( templC(0) periphCf) + φΔ( ridgeC

0.9(syst)± 8.5(stat) ± = 15.4 310×2,2v