# Prospects for Open Heavy Flavor Measurements in Heavy Ion andp

Authors: Daniel Kikoła
Journal: Advances in High Energy Physics
 (2015)

---

## Abstract

High luminosity data in a fixed-target experiment allow studying interactions of heavy quarks with nuclear matter in the intermediate energy range with extremely high precision. We present a feasibility study for open charm and bottom production measurements in the energy range of a fixed-target experiment at the LHC (AFTER@LHC). We demonstrate, that high-precision data from AFTER will allow answering two open questions: if there is a collective behavior of charm quarks inp
+
A collisions at RHIC energy and if charm production is suppressed in the energy range of s
N
N = 60–80 GeV. We argue that simultaneous measurement of D
0 suppression as a function of traverse momentum at midrapidity and forward rapidity can help to pin down the mechanism of charm energy loss in the hot and dense nuclear medium.


---

## Body

## 1. Introduction

Relativistic heavy ion collisions provide a unique opportunity for studying quark-gluon plasma (QGP), a new state of nuclear matter with properties determined by quark and gluon degrees of freedom. Its properties can be studied with heavy quarks: charm and bottom. They are produced very early in the collision, in the initial interactions with large momentum transfer. Their production and both total and differential cross sections are well described by perturbative QCD calculations. Moreover, due to large masses, they are expected to interact with the nuclear matter differently compared to light partons.We can infer properties of the QGP by studying modification of heavy quarks production in the heavy ion collisions compared to proton-proton baseline (for details, see [1] and references therein). Energy loss and elliptic flow of open heavy flavor are sensitive to the dynamics of the medium; such measurements could be used to determine the fundamental properties of the QGP, for instance, transport coefficients. Elliptic flow of heavy quarks can give insights into degree of thermalization of the created nuclear matter and can help to discriminate between different models of heavy quark interactions with the QGP.Heavy flavor production at RHIC (measured viaD
0 and electrons from semileptonic decays of heavy flavor hadrons, e
H
F) in Au + Au collisions at s
N
N = 200 GeV is suppressed at high transverse momentum (p
T) and e
H
F have a significant elliptic flow v
2 (elliptic flow is defined as the second harmonic in the Fourier expansion of the particle azimuthal distribution with respect to the reaction plane [2]). These results suggest that charm quarks are strongly coupled with the medium: heavy quarks loose a large amount of energy and acquire a significant elliptic flow during interaction with matter created at top RHIC energy.Beam energy scan (BES) program at RHIC was carried out recently to study the phase diagram of nuclear matter and search for a phase transition and a critical point. BES results show that elliptic flow of inclusive charged hadrons is approximately independent of beam energy (the difference is less than 10% for0.5
<
p
T
<
3 GeV/c) and light hadron production is suppressed at high p
T in the energy range of 39–200 GeV. In Au + Au collisions at s
N
N = 62.4 GeV, production of e
H
F is not suppressed (within large uncertainties). Moreover, measurements of elliptic flow of e
H
F in Au + Au collisions at s
N
N = 62.4 and 39 GeV [3] showed that v
2 of e
H
F is consistent with zero. This is a different behavior than that for light hadrons, where a positive v
2 is observed and the difference between s
N
N = 200 GeV and 39 GeV is small. These results suggest that there is a difference in the interactions of the heavy quarks with the surrounding nuclear matter at 200 GeV compared to the two lower energies and there is already a change of the nuclear matter properties in the energy range of 62.4–200 GeV. High luminosity data of a proposed fixed-target experiment at the LHC (s
N
N = 72 GeV) together with a wide rapidity coverage could shed new light on the energy loss mechanism for charm quarks (which we discuss in Section 5) and give a precise answer if heavy flavor production is suppressed or not at high p
T.For the interpretation of the results from heavy ion collisions, it is important to have a good handle on the so-called cold nuclear matter (CNM) effects. This category includes modifications of the particle yields not related to the QGP formation, such as a shadowing (modification of parton distributions in a nucleus) or a Cronin effect (an enhancement of transverse momentum inp
+
A with respect to p
+
p collisions). Experiments at RHIC and LHC recently reported collective behavior of light hadrons in high multiplicity d + Au [4] and p + Pb collisions [5], where these hadrons have significant elliptic flow with a characteristic mass-splitting pattern [6]. These observations triggered speculations that an enhancement of e
H
F production in central and minimum bias d + Au collisions at midrapidity at RHIC [7] could indicate a collective phenomenon (radial flow) of heavy quarks in d + Au [8]. However, this enhancement could be also due to the CNM effects: it can be reproduced assuming Cronin effect for charm quarks [9]. A direct measurement of elliptic flow v
2 of D
0 will answer the question of heavy quark collectivity inp +A collisions. High luminosityp +A collisions at AFTER will allow quantifying the CNM effects and measuring v
2 with high precision to address these two issues.In this paper we present estimates of the precision expected for open heavy flavor measurements at the proposed fixed target experiment at LHC (AFTER@LHC). We first describe the simulation setup we use and then present and discuss estimates for physical observables for open charm and open bottom inp
+
A collisions at s
N
N = 115 GeV and Pb + Pb collisions at s
N
N = 72 GeV. We argue that these measurements at AFTER will address open questions about heavy quark collectivity in p
+
A collisions and energy dependence of the heavy quark interactions with a nuclear matter.


## 2. Centrality Estimation

We estimate centrality bins inp + Pb collisions at s
N
N = 115 GeV and P
b
+
P
b at s
N
N = 72 GeV using Glauber Monte Carlo calculations with the PHOBOS Glauber Monte Carlo [10, 11] version 2.1. We assume that the nucleon-nucleon cross section σ
N
N
=
37mb at s
N
N = 72 GeV and the nucleon-nucleon cross section σ
N
N
=
39mb at s
N
N = 115 GeV, and the rest of parameters have the default values.Centrality bins are defined as percentiles of hadronic cross sections. Figure1 shows impact parameter distributions in p-Pb collisions at s
N
N = 115 GeV and (b) Pb-Pb collisions at s
N
N = 72 GeV. Tables 1 and 2 show the average number of participants N
p
a
r
t and binary collisions N
b
i
n and their root mean square (RMS) values for those two systems for the centrality classes.Table 1
Average number of participants〈
N
p
a
r
t
〉 and binary collisions 〈
N
b
i
n
〉 in p
-
P
b collisions at  s
N
N
=
115 GeV for centrality classes defined as percentiles of the hadronic cross section. The mean values and the RMS are obtained with a Glauber Monte Carlo calculation.
Centrality
b
m
i
n [fm]
b
m
a
x [fm]
〈
N
b
i
n
〉
R
M
S
(
N
b
i
n
)
〈
N
p
a
r
t
〉
R
M
S
(
N
p
a
r
t
)
0–10%
0
2.41
8.0
2.7
9.0
2.7
10–20%
2.41
3.42
7.3
2.6
8.3
2.6
20–30%
3.42
4.19
6.6
2.4
7.6
2.4
30–40%
4.19
4.83
5.8
2.3
6.8
2.3
40–50%
4.83
5.41
4.9
2.1
5.9
2.1
50–60%
5.41
5.93
3.9
1.8
4.9
1.8
60–70%
5.93
6.43
2.9
1.5
3.9
1.5
70–80%
6.43
6.95
2.2
1.2
3.2
1.2
80–90%
6.95
7.60
1.6
0.8
2.6
0.8
90–100%
7.60
19.80
1.2
0.4
2.2
0.4Table 2
Average number of participants〈
N
p
a
r
t
〉 and binary collisions 〈
N
b
i
n
〉 in Pb-Pb collisions at s
N
N
=
72 GeV for centrality classes defined as percentiles of the hadronic cross section. The mean values and the RMS are obtained with a Glauber Monte Carlo calculation.
Centrality
b
m
i
n [fm]
b
m
a
x [fm]
〈
N
b
i
n
〉
R
M
S
(
N
b
i
n
)
〈
N
p
a
r
t
〉
R
M
S
(
N
p
a
r
t
)
0–10%
0
4.80
884.0
128.0
338.5
34.0
10–20%
4.80
6.77
557.2
82.2
241.8
25.6
20–30%
6.77
8.30
346.5
58.9
170.8
20.8
30–40%
8.30
9.59
206.2
42.9
117.1
17.2
40–50%
9.59
10.72
115.4
30.6
76.7
14.6
50–60%
10.72
11.75
59.4
20.6
47.0
12.0
60–70%
11.75
12.69
28.0
13.0
26.6
9.5
70–80%
12.69
13.58
12.2
7.5
13.7
6.8
80–90%
13.58
14.55
5.2
4.0
6.8
4.2
90–100%
14.55
20.00
2.3
1.8
3.5
2.1Figure 1
Impact parameter distribution for percentiles of the hadronic cross section obtained from a Glauber Monte Carlo calculation in (a)p-Pb collisions at s
N
N
=
115 GeV and (b) Pb-Pb collisions at s
N
N
=
72 GeV.
(a)
(b)


## 3. Simulation Setup

### 3.1. Detector Setup

We consider a detector with similar apparatus as the LHCb experiment [12]. We assume rapidity coverage in the laboratory frame of 2
<
y
L
a
b
<
5 and similar particle detector capabilities and efficiencies (muon detection capability and precise microvertexing detector). The efficiencies reported by LHCb [13, 14] are used in our estimates (Figure 2). We consider two experimental techniques for open heavy flavor measurements: D
0 measurements via secondary vertex reconstruction and B
→
J
/
ψ determination via a pseudo-lifetime method.Figure 2
Total reconstruction efficiency for (a)D
0 [13] and (b) J
/
ψ [14] in the LHCb experiment.
(a)
(b)Since bothD
0 and J
/
ψ from B-hadron decays are measured at the forward rapidity in a laboratory frame with a precise vertex detector, a large boost will strongly suppress the background. Thus we neglect the background contribution to the statistical uncertainties in our estimations. The large boost also allows for D and B mesons and B
→
J
/
ψ measurements down to zero p
T via displaced vertex reconstruction.


### 3.2. Heavy Quarks Production Cross Sections

We obtain charm and bottom quark production cross sections using FONLL (fixed-order plus next-to-leading logs) program [15, 16], version 1.3.3 [17]. We run the calculation assuming charm quark mass m
c
=
1.5GeV/c2 and bottom mass m
b
=
4.75GeV/c2 using CTEQ6.6 [18] parton distribution function convoluted with EPS09 [19] shadowing parametrization (central value). Other parameters (scales) have the default values.Figures3 and 4 show charm and bottom quarks differential cross section per nucleon-nucleon collision as a function of p
T in p + Pb and Pb + Pb collisions. Only central values from FONLL are shown. Since experimental data for heavy flavor production at RHIC energies are consistent with upper limit on the FONLL calculations [20, 21] we believe that the results in Figures 3 and 4 provide us with conservative estimates of heavy quark cross section for those energies.Figure 3
Charm quark differential cross section per nucleon-nucleon collision as a function ofp
T inp + Pb and Pb + Pb collisions.
(a)
(b)Figure 4
Bottom quark differential cross section per nucleon-nucleon collision as a function ofp
T in p + Pb collisions.
(a)
(b)


### 3.3. Charmed and Bottom Mesons Production

In the experiments we measure the final-state particles (D
0, B-meson or J
/
ψ from B-meson decays). In this study we made a few simplifications to estimate the expected yields based on charm and bottom cross sections. For D
0 meson, we assume that it has approximately the same p
T spectrum as charm quarks. Charm hadronization ratio to D
0 is 0.565 [22] and D
0 will be measured via D
0
→
K
-
π
+ (branching ratio BR = 3.8%) via secondary vertex reconstruction. We used TPythia8Decayer from the ROOT framework [23] to simulate B
→
J
/
ψ decays. For the input, we assumed that B
±, B
S, and B
0 have the same p
T spectrum as bottom quarks, rapidity has a uniform distribution in narrow bins we have used (Δ
y
=
0.5), and the hadronization fraction to B meson is f
(
b
→
B
)
=
0.764 (f
(
b
→
B
0
)
=
f
(
b
→
B
±
)
=
(
33.7
±
2.2
)
% and f
(
b
→
B
S
)
=
(
9.0
±
0.9
)
% [24]). The branching ratio for B
→
J
/
ψ
+
X is taken as B
R
=
1.16
%. Figure 4(b) shows an input bottom quark spectrum and J
/
ψ from bottom meson decays for y
<
0.25. At the energy of s
N
N = 115 GeV, we expect bottom quark production with p
T up to 11 GeV/c, which allows for B
→
J
/
ψ measurement with p
T up to 10 GeV/c.To calculate the production rates, we take an integrated luminosityL
=
160
p
b
-
1 for p + Pb collisions at s
N
N = 115 GeV and L
=
7
n
b
-
1 for Pb + Pb collisions at s
N
N = 72 GeV [25]. Finally, we apply the charmed meson and J
/
ψ reconstruction efficiencies reported by LHCb (Figure 2) to the simulated D
0 and B
→
J
/
ψ distributions to estimate the yields and statistical uncertainties expected in the experiment.


### 3.4. Energy Loss of Charm Quarks

The suppression of open heavy flavor production observed at RHIC is comparable to that for light hadrons. This was a surprising result at the beginning because models which described well the light flavor data assuming gluon radiation (radiative energy loss,d
E
/
d
x
r
a
d) predicted that the heavy quark production will be less suppressed than light partons [26]. This led to introduction of collisional energy loss d
E
/
d
x
c
o
l
l due to binary interactions of partons with other objects in the QGP. One difficulty in modeling heavy quark energy loss is that the relative contributions of d
E
/
d
x
c
o
l
l and d
E
/
d
x
r
a
d are unknown. To understand better the interplay of these two processes, precise measurements of suppression of charm and bottom separately are necessary. We argue that measurements in Pb + Pb collisions at s
N
N = 72 GeV in different rapidity ranges can help to estimate what are the shares of d
E
/
d
x
r
a
d and d
E
/
d
x
c
o
l
l in the observed suppression.To model charm quark energy loss, we use predictions from [27] for radiative and collisional energy loss of heavy quark. Figure 5 shows the expected energy loss per unit of distance traveled in a nuclear matter for charm quarks from model I in [27] with a coupling constant α
s
=
0.3, an infrared regulator μ
=
0.4GeV, and temperature T
=
200 MeV (see [27] for details). For low energies of heavy quarks, d
E
/
d
x
c
o
l
l dominates but the radiative energy loss increases fast with increasing heavy quark energy and d
E
/
d
x
r
a
d overshadows d
E
/
d
x
c
o
l
l for E
>
5GeV.Figure 5
Radiative and collisional energy loss for charm quarks from [27] (model I) as a function of charm quark energy. Adopted from [27].A large acceptance of planned AFTER experiment allows probing different regimes of the energy loss, despite limitedD
0
p
T range available experimentally, because charm quark at forward/backward rapidity will have a larger overall energy for the same p
T compared to midrapidity.To model the suppression of heavy quark production, we assume that the average path length in the nuclear matterL is approximately equal to the averaged transverse path length L
x
y (distance traveled in a plane perpendicular to the beam axis) traversed by a quark from a production point to the edge of the created fireball. We adopted a definition of L
=
σ
x
2
+
σ
y
2 from [28], where σ
x
2 and σ
y
2 are variances of x and y positions of the participant nucleons, respectively. σ
x
2 and σ
y
2 are obtained using PHOBOS Glauber Monte Carlo and L
=
4.104 fm for 0–10% most central collisions. We use a nuclear modification factor R
A
A to quantify the modification of the charm quark p
T spectrum due to energy loss:(1)
R
A
A
p
T
=
N
E
Loss
A
A
p
T
N
un.mod
A
A
p
T
,where N
E
L
o
s
s is a yield with energy loss applied and N
u
n
.
m
o
d is a heavy quark yield without any energy loss for a given centrality class and a given rapidity bin.Figure6 shows a nuclear modification factor R
A
A (central value) at midrapidity (Figure 6(a)) and backward rapidity (Figure 6(b)) for 0–10% most central Pb + Pb collisions expected for d
E
/
d
x
c
o
l
l and d
E
/
d
x
r
a
d shown in Figure 5. The difference between d
E
/
d
x
r
a
d and d
E
/
d
x
c
o
l
l is more pronounced at backward rapidity, even at relatively low p
T, because quarks have much larger energy due to larger longitudinal momentum. Measurement at these two rapidity ranges simultaneously could help to pin down interplay of these two mechanisms.Figure 6
Nuclear modification factorR
A
A (central value) for charm quark in Pb + Pb collisions at s
N
N = 72 GeV at midrapidity (y
C
M
S
=
0) and forward rapidity (y
C
M
S
=
-
2.05) for a radiative and collisional energy loss from Figure 5.
(a)
(b)


## 3.1. Detector Setup

We consider a detector with similar apparatus as the LHCb experiment [12]. We assume rapidity coverage in the laboratory frame of 2
<
y
L
a
b
<
5 and similar particle detector capabilities and efficiencies (muon detection capability and precise microvertexing detector). The efficiencies reported by LHCb [13, 14] are used in our estimates (Figure 2). We consider two experimental techniques for open heavy flavor measurements: D
0 measurements via secondary vertex reconstruction and B
→
J
/
ψ determination via a pseudo-lifetime method.Figure 2
Total reconstruction efficiency for (a)D
0 [13] and (b) J
/
ψ [14] in the LHCb experiment.
(a)
(b)Since bothD
0 and J
/
ψ from B-hadron decays are measured at the forward rapidity in a laboratory frame with a precise vertex detector, a large boost will strongly suppress the background. Thus we neglect the background contribution to the statistical uncertainties in our estimations. The large boost also allows for D and B mesons and B
→
J
/
ψ measurements down to zero p
T via displaced vertex reconstruction.


## 3.2. Heavy Quarks Production Cross Sections

We obtain charm and bottom quark production cross sections using FONLL (fixed-order plus next-to-leading logs) program [15, 16], version 1.3.3 [17]. We run the calculation assuming charm quark mass m
c
=
1.5GeV/c2 and bottom mass m
b
=
4.75GeV/c2 using CTEQ6.6 [18] parton distribution function convoluted with EPS09 [19] shadowing parametrization (central value). Other parameters (scales) have the default values.Figures3 and 4 show charm and bottom quarks differential cross section per nucleon-nucleon collision as a function of p
T in p + Pb and Pb + Pb collisions. Only central values from FONLL are shown. Since experimental data for heavy flavor production at RHIC energies are consistent with upper limit on the FONLL calculations [20, 21] we believe that the results in Figures 3 and 4 provide us with conservative estimates of heavy quark cross section for those energies.Figure 3
Charm quark differential cross section per nucleon-nucleon collision as a function ofp
T inp + Pb and Pb + Pb collisions.
(a)
(b)Figure 4
Bottom quark differential cross section per nucleon-nucleon collision as a function ofp
T in p + Pb collisions.
(a)
(b)


## 3.3. Charmed and Bottom Mesons Production

In the experiments we measure the final-state particles (D
0, B-meson or J
/
ψ from B-meson decays). In this study we made a few simplifications to estimate the expected yields based on charm and bottom cross sections. For D
0 meson, we assume that it has approximately the same p
T spectrum as charm quarks. Charm hadronization ratio to D
0 is 0.565 [22] and D
0 will be measured via D
0
→
K
-
π
+ (branching ratio BR = 3.8%) via secondary vertex reconstruction. We used TPythia8Decayer from the ROOT framework [23] to simulate B
→
J
/
ψ decays. For the input, we assumed that B
±, B
S, and B
0 have the same p
T spectrum as bottom quarks, rapidity has a uniform distribution in narrow bins we have used (Δ
y
=
0.5), and the hadronization fraction to B meson is f
(
b
→
B
)
=
0.764 (f
(
b
→
B
0
)
=
f
(
b
→
B
±
)
=
(
33.7
±
2.2
)
% and f
(
b
→
B
S
)
=
(
9.0
±
0.9
)
% [24]). The branching ratio for B
→
J
/
ψ
+
X is taken as B
R
=
1.16
%. Figure 4(b) shows an input bottom quark spectrum and J
/
ψ from bottom meson decays for y
<
0.25. At the energy of s
N
N = 115 GeV, we expect bottom quark production with p
T up to 11 GeV/c, which allows for B
→
J
/
ψ measurement with p
T up to 10 GeV/c.To calculate the production rates, we take an integrated luminosityL
=
160
p
b
-
1 for p + Pb collisions at s
N
N = 115 GeV and L
=
7
n
b
-
1 for Pb + Pb collisions at s
N
N = 72 GeV [25]. Finally, we apply the charmed meson and J
/
ψ reconstruction efficiencies reported by LHCb (Figure 2) to the simulated D
0 and B
→
J
/
ψ distributions to estimate the yields and statistical uncertainties expected in the experiment.


## 3.4. Energy Loss of Charm Quarks

The suppression of open heavy flavor production observed at RHIC is comparable to that for light hadrons. This was a surprising result at the beginning because models which described well the light flavor data assuming gluon radiation (radiative energy loss,d
E
/
d
x
r
a
d) predicted that the heavy quark production will be less suppressed than light partons [26]. This led to introduction of collisional energy loss d
E
/
d
x
c
o
l
l due to binary interactions of partons with other objects in the QGP. One difficulty in modeling heavy quark energy loss is that the relative contributions of d
E
/
d
x
c
o
l
l and d
E
/
d
x
r
a
d are unknown. To understand better the interplay of these two processes, precise measurements of suppression of charm and bottom separately are necessary. We argue that measurements in Pb + Pb collisions at s
N
N = 72 GeV in different rapidity ranges can help to estimate what are the shares of d
E
/
d
x
r
a
d and d
E
/
d
x
c
o
l
l in the observed suppression.To model charm quark energy loss, we use predictions from [27] for radiative and collisional energy loss of heavy quark. Figure 5 shows the expected energy loss per unit of distance traveled in a nuclear matter for charm quarks from model I in [27] with a coupling constant α
s
=
0.3, an infrared regulator μ
=
0.4GeV, and temperature T
=
200 MeV (see [27] for details). For low energies of heavy quarks, d
E
/
d
x
c
o
l
l dominates but the radiative energy loss increases fast with increasing heavy quark energy and d
E
/
d
x
r
a
d overshadows d
E
/
d
x
c
o
l
l for E
>
5GeV.Figure 5
Radiative and collisional energy loss for charm quarks from [27] (model I) as a function of charm quark energy. Adopted from [27].A large acceptance of planned AFTER experiment allows probing different regimes of the energy loss, despite limitedD
0
p
T range available experimentally, because charm quark at forward/backward rapidity will have a larger overall energy for the same p
T compared to midrapidity.To model the suppression of heavy quark production, we assume that the average path length in the nuclear matterL is approximately equal to the averaged transverse path length L
x
y (distance traveled in a plane perpendicular to the beam axis) traversed by a quark from a production point to the edge of the created fireball. We adopted a definition of L
=
σ
x
2
+
σ
y
2 from [28], where σ
x
2 and σ
y
2 are variances of x and y positions of the participant nucleons, respectively. σ
x
2 and σ
y
2 are obtained using PHOBOS Glauber Monte Carlo and L
=
4.104 fm for 0–10% most central collisions. We use a nuclear modification factor R
A
A to quantify the modification of the charm quark p
T spectrum due to energy loss:(1)
R
A
A
p
T
=
N
E
Loss
A
A
p
T
N
un.mod
A
A
p
T
,where N
E
L
o
s
s is a yield with energy loss applied and N
u
n
.
m
o
d is a heavy quark yield without any energy loss for a given centrality class and a given rapidity bin.Figure6 shows a nuclear modification factor R
A
A (central value) at midrapidity (Figure 6(a)) and backward rapidity (Figure 6(b)) for 0–10% most central Pb + Pb collisions expected for d
E
/
d
x
c
o
l
l and d
E
/
d
x
r
a
d shown in Figure 5. The difference between d
E
/
d
x
r
a
d and d
E
/
d
x
c
o
l
l is more pronounced at backward rapidity, even at relatively low p
T, because quarks have much larger energy due to larger longitudinal momentum. Measurement at these two rapidity ranges simultaneously could help to pin down interplay of these two mechanisms.Figure 6
Nuclear modification factorR
A
A (central value) for charm quark in Pb + Pb collisions at s
N
N = 72 GeV at midrapidity (y
C
M
S
=
0) and forward rapidity (y
C
M
S
=
-
2.05) for a radiative and collisional energy loss from Figure 5.
(a)
(b)


## 4. Nuclear Modification Factors and Elliptic Flow inp + Pb Collisions at (
s
N
N
) = 115 GeV

The sampled luminosity calculation inp
+
p collisions is a significant source of the uncertainty when an energy loss is studied comparing heavy ion and p
+
p collisions. Alternatively, we can use a central-to-peripheral ratio R
C
P to quantify modification of the yields in p
+
A or A
+
A collisions:(2)
R
C
P
p
T
=
N
bin
P
e
r
i
p
N
bin
Central
N
C
e
n
t
r
a
l
p
T
N
P
e
r
i
p
p
T
,where N
Central and N
P
e
r
i
p are yields in central and peripheral A
+
A collisions, respectively, for a given rapidity range and N
b
i
n
C
e
n
t
r
a
l and N
bin
P
e
r
i
p are the average numbers of binary collisions for central and peripheral collisions. We use 0–10% most central collisions as a central bin and 60–80% bin as a baseline in R
C
P calculations.Figure8(a) shows the expected statistical uncertainties for nuclear modification factor R
C
P for D
0 at midrapidity and backward rapidity. Figure 7 presents the results for B
→
J
/
ψ. The expected precision will allow a precision study of the cold nuclear matter effects for both charm and bottom for p
T
<
6GeV/c. Figure 8(b) shows statistical uncertainties for elliptic flow measurement. High luminosity data expected from AFTER would give a precise answer if there are any signs of collective behavior of charm quarks in p
+
A collisions at the intermediate energy at RHIC.Figure 7
Expected statistical uncertainty ofR
C
P for B
→
J
/
ψ in p + Pb collisions at s
N
N = 115 GeV.Figure 8
Expected statistical uncertainties of (a)D
0
R
C
P and (b) D
0 elliptic flow in p + Pb collisions at s
N
N = 115 GeV.
(a)
(b)


## 5. Nuclear Modification Factor for Charmed Meson in Pb + Pb Collisions ats
N
N = 72 GeV

Figure9 shows the expected precision of R
C
P for charmed meson as a function of p
T at midrapidity and backward rapidity. We show estimates for pure collisional and radiative energy loss scenarios based on Figure 5. The expected integrated luminosity for a single year will allow addressing the question if charm production is suppressed at high transverse momentum at s
N
N
=
72 GeV. The comparison of results at different rapidity ranges will help to pin down the interplay of energy loss mechanisms at this energy range. The measurements of D
0 suppression for various collisions systems (which is feasible in a multiple-year physics program at AFTER) are even more interesting since they will allow studying the energy loss as a function of path length and energy density of the created system.Figure 9
Expected statistical uncertainties of nuclear modification factorR
C
P for D
0 meson in Pb + Pb collisions at s
N
N = 72 GeV.
(a)
(b)


## 6. Summary

We presented prospects for open charm and bottom production measurements in a fixed target experiment at LHC. We argue that such measurements will address important open issues about interactions of heavy quarks with the nuclear matter. High luminosityp
+
A and Pb + Pb data will help answer questions if there is a collective behavior of heavy quark in p
+
A collisions at RHIC and what is the mechanism of energy loss of heavy quarks. Measurements at s
N
N
=
72 GeV can shed new light on energy dependence of heavy quark interaction with the nuclear matter and thus on the phase diagram of nuclear matter.



---

*Source: 783134-2015-10-04.xml*