# Feasibility Studies for Quarkonium Production at a Fixed-Target Experiment Using the LHC Proton and Lead Beams (AFTER@LHC)

Authors: L. Massacrier; B. Trzeciak; F. Fleuret; C. Hadjidakis; D. Kikola; J. P. Lansberg; H.-S. Shao
Journal: Advances in High Energy Physics
 (2015)

---

## Abstract

Being used in the fixed-target mode, the multi-TeV LHC proton and lead beams allow for studies of heavy-flavour hadroproduction with unprecedented precision at backward rapidities, far negative Feynman-x, using conventional detection techniques. At the nominal LHC energies, quarkonia can be studied in detail in p+p, p+d, and p+A collisions at sNN≃115 GeV and in Pb + p and Pb + A collisions at sNN≃72 GeV with luminosities roughly equivalent to that of the collider mode that is up to 20 fb−1 yr−1 in p+p and p+d collisions, up to 0.6 fb−1 yr−1 in p+A collisions, and up to 10 nb−1 yr−1 in Pb + A collisions. In this paper, we assess the feasibility of such studies by performing fast simulations using the performance of a LHCb-like detector.


---

## Body

## 1. Introduction

Since its start-up, the large hadron collider (LHC), the most energetic hadron collider ever built so far, has already made the demonstration of its outstanding capabilities. These can greatly be complemented by the addition of a fixed-target physics program. Its multi-TeV beams indeed allow one to studyp+p, p+d, and p+A collisions at a center-of-mass (c.m.s.) energy sNN≃115 GeV as well as Pb+p and Pb+A collisions at sNN≃72 GeV, with the high precision typical of the fixed-target mode. In this context, the proposal of a fixed-target experiment at the LHC [1], referred to as AFTER@LHC, has been promoted [1] in order to complement the existing collider experiments such as the Relativistic Heavy Ion Collider (RHIC) or the future Electron-Ion Collider (EIC) project in a similar energy range. The idea underlying the AFTER@LHC proposal is a multipurpose detector allowing for the study of a multitude of probes.Various technological ways to perform fixed-target experiment at the LHC exist. On the one hand, the beam can be extracted by means of a bent crystal. This technology [2, 3] is currently developed as a smart beam-collimation solution and is studied by the UA9/LUA9 collaboration, respectively, at SPS and LHC. A bent crystal installed in the halo of the LHC beam would deflect the particles of the halo onto a target, with a flux of 5×108 proton/s without any impact on the LHC performances [3–5].On the other hand, the LHC beam can go through an internal-gas-target system in an existing (or new) LHC experiment. Such a system is already tested at low gas pressure by the LHCb collaboration to monitor the luminosity of the beam [6–8]. Data were taken at a center-of-mass energy of sNN=87 (54) GeV with p+Ne (Pb+Ne) collisions during pilot runs in 2012 and 2013. Although this system, called SMOG, was tested during only few hours in a row, no decrease of the LHC performances was observed.In the bent crystal case, the luminosity achievable with AFTER@LHC would surpass that of RHIC by 3 orders of magnitudes [1]. We have reported in Table 1 the instantaneous and yearly integrated luminosities expected with the proton and Pb beams on various target species of various thicknesses, for the bent crystal as well as internal-gas-target options. Integrated luminosities as large as 20 fb−1 can be delivered during a one-year run of p+H collisions with a bent crystal. Besides, it is worth mentioning that both technologies allow one to polarise the target, which is an important requirement to lead an extensive spin physics programme [1, 11].Table 1
Expected luminosities obtained for a 7 (2.76) TeV proton beam (Pb) extracted by means of a bent crystal or obtained with an internal-gas-target system.
Beam
Target
Thickness
ρ
L
∫
L
(cm)
(g⋅cm−3)
(μb−1⋅s−1)
(pb−1⋅y−1)
p
Liquid H
100
0.068
2000
20000
p
Liquid D
100
0.16
2400
24000
p
Pb
1
11.35
16
160
Pb
Liquid H
100
0.068
0.8
0.8
Pb
Liquid D
100
0.16
1
1
Pb
Pb
1
11.35
0.007
0.007
Beam
Target
Usable gas zone
Pressure
L
∫
L
(cm)
(Bar)
(μb−1⋅s−1)
(pb−1⋅y−1)
p
Perfect gas
100
1
0
-
9
10
100
Pb
Perfect gas
100
1
0
-
9
0.001
0.001Overall, thanks to the large luminosity expected, AFTER@LHC would become a quarkonium [12], prompt photon, and heavy-flavour observatory [1, 13] in p+p and p+A collisions where, by instrumenting the target-rapidity region, gluon and heavy-quark distributions of the proton, the neutron and the nuclei can be accessed at large x and even at x larger than unity in the nuclear case [14]. In addition, the fixed-target mode allows for single-target-spin-asymmetry measurements over the full backward rapidity domain up to xF≃-1 [15, 16]. Also, the versatility in the target choices offer a unique opportunity to study the nuclear matter versus the hot and dense matter formed in heavy ion collisions which can be studied during the one-month lead run. In the latter case, modern detection technology (such as high granularity calorimeter) should allow for extensive studies of quarkonium excited states, from ψ(2S) to χc and χb resonances thanks to the boost of the fixed-target mode [17].In this paper, we report on a feasibility study of quarkonium production at a fixed-target experiment using LHC beams. In Section2, we outline the simulation framework which was used. In Section 3, we describe how a fast simulation of a detector response has been implemented, following a LHCb-like detector setup. In Section 4, we present the charmonium and bottomonium family studies performed with the p+H simulations at s=115 GeV. In Section 5, we present multiplicity studies in p+A and A+p collisions as well as the expected nuclear modification factors for J/ψ and Υ in p+Pb collisions at sNN=115 GeV. Finally in Section 6 some prospects for Pb+A measurement at sNN=72 GeV are given. Section 7 gathers our conclusions.


## 2. Simulation Inputs

In order to get the most realistic minimum bias simulations at AFTER@LHC energy for quarkonium studies in the dimuon decay channels, we have simulated the quarkonium signal and all the background sources separately to have under control the transverse momentum and rapidity input distributions as well as the normalisation of the different sources.The simulation has been performed forp+p collisions at s=115 GeV. On the one hand, the quarkonium signal and the correlated background (Drell-Yan, cc-, bb-) were simulated with HELAC-Onia [18] which produces outputs following the format of Les Houches Event Files [19]. The outputs were then processed with Pythia (Pythia 8.185 [20]) to perform the hadronisation, the initial/final-state radiations, and the decay of the resonances. On the other hand, the uncorrelated background was obtained from minimum bias p+p collisions generated with Pythia.The relative normalisation of the signal and background sources was performed according to the production cross section of the process (taking into account initial phase space cuts, if any). Values of the cross section and the number of simulated eventsNsim, not to be confused with the expected events for a specific luminosity, are reported in Table 2. The cross section values are integrated over rapidity and pT.Table 2
Total cross section for different processes inp+p collisions at s = 115 GeV and number of simulated events Nsim.
σ
t
o
t (mb)
N
s
i
m
J/ψ
1.30× 10−3
1.47× 106
ψ(2S)
1.61× 10−4
1.12× 106
Υ(1S)
4.30× 10−7
1.46× 106
Υ(2S)
1.22× 10−7
1.49× 106
Υ(3S)
5.28× 10−8
1.48× 106
Drell-Yan (M> 2.5 GeV/c2)
2.52× 10−6
4.3× 105
Drell-Yan (M> 7 GeV/c2)
1.49× 10−7
2.0× 106
c
c
-
2.29× 10−1
81.5× 106
b
b
-
4.86× 10−4 (gg→bb-)
32.3× 106 (gg→bb-)
1.49× 10−4 (qq-→bb-)
85.7× 106 (qq-→bb-)
Minimum bias
26.68
11.0× 108


### 2.1. Signal and Correlated Background

#### 2.1.1. Quarkonium Signal

J
/
ψ, ψ(2S), Υ(1S), Υ(2S), and Υ(3S) were simulated in a data-driven way. The amplitude of gg→Q+X (where Q is the quarkonium) is expressed in an empirical functional form [21]: (1)Agg→Q+X2¯=Kexp⁡-κpT2MQ2whenpT≤pTKexp⁡-κpT2MQ21+κnpT2-pT2MQ2-nwhenpT>pT,where K=λ2κs^/MQ2 with s^ the partonic center-of-mass energy and MQ the mass of the quarkonium Q taken from the PDG table [22].The parametersκ, λ, n, and 〈pT〉 were determined by fitting the differential cross section d2σ/dpTdy to the experimental data. The dedicated codes used to perform the fit and to generate unweighted events for quarkonium production have been implemented in HELAC-Onia [18] and we used MSTW2008NLO PDF set [23] provided in LHAPDF5 [24] and the factorisation scale μF=MQ2+pT2. In order to constrain the nontrivial energy dependence of quarkonium production, we used the differential measurements of charmonium production performed by the PHENIX collaboration at RHIC, in p+p collisions at s=200 GeV [9] to predict the corresponding yields at s=115 GeV. Given the lack of such measurements for Υ at RHIC, we performed a combined fit to CDF [25], ATLAS [26], CMS [27], and LHCb [10, 28] data on Υ production. The values of the fitted parameters are listed in Table 3. For illustration, the comparison between fits and the selected experimental data is shown in Figure 1.Table 3
Fit parameters obtained after a combined fit ofd2σ/dpTdy to the PHENIX data [9] for charmonium production and to CDF [25], ATLAS [26], CMS [27], and LHCb [10, 28] data for bottomonium production. We have fixed n=2 and pT= 4.5 (13.5) GeV/c for charmonium (bottomonium) production. The number of fitted data points is also reported.
κ
λ
Number of data points
χ
2
J
/
ψ
0.674
0.380
51
422
ψ
(
2
S
)
0.154
0.351
4
1.12
Υ
(
1
S
)
0.707
0.0837
288
1883
Υ
(
2
S
)
0.604
0.0563
205
856
Υ
(
3
S
)
0.591
0.0411
197
886Figure 1
Some illustrative comparisons between fits and the PHENIX data [9] for charmonium production (a) and LHCb data [10] for bottomonium production (b).
(a)
(b)In order to increase the statistics of the simulated data sample, the decay of the quarkonium in Pythia is forced into the dimuon decay channel. The simulated yields are then weighted by the cross section for this process multiplied by the Branching Ratio (BR).


#### 2.1.2. Open Charm

Open charm production was simulated with the processgg→cc- in HELAC-Onia. In order to avoid the huge theoretical uncertainties in the state-of-the-art perturbative calculations, open charm yields at s=115 GeV are also computed in a data-driven way following the method described in the previous section. Similarly, the matrix element of gg→cc- is determined using (1). The parameters are obtained from a fit to the pT-differential cc- cross section measured by the STAR experiment [29] in p+p collisions at s=200 GeV (see Figure 3). We obtained κ=0.437, λ=3.04, and 〈pT〉=2.86 GeV/c when n=2 by using CTEQ6L1 [30] and by fixing the c quark mass to mc=1.5 GeV/c2 and the factorisation scale to μF=mc2+pT2. The χ2 of the fit is equal to 4.39 with 10 experimental data points. The tuned result is shown in Figure 2. The evolution of the cross section with the energy down to s=115 GeV is then given by HELAC-Onia.Figure 2
A comparison between fit and the STAR data [29] in p+p collisions at s=200 GeV, for cc- production.Figure 3
Drell-Yan cross section as a function of the rapidity in the center-of-mass frame obtained with HELAC-Onia + Pythia at s=38.76 GeV and rescaled by a factor 1.2, together with E866 data extracted from [33]. The invariant mass range considered is 7.2<M<8.7 GeV/c2.After embedding the Les Houches Event File into Pythia, muons from the underlying Pythia event can be produced on top of muons from the initial cc- pair. The combination of those additional muons with a muon from the initial cc- pair is not included in our definition of open charm correlated background. We have however checked that this contribution is negligible. In order to increase the statistics, D0, D0-, D+/-, and Ds+/- were forced to decay into muons and only those decay muons were considered as correlated background. μ+μ- pairs coming from all possible combinations, D0D0-, D+D-, Ds+Ds-, D0D+/-, D0Ds+/-, and D+/-Ds-/+, are considered. The simulated events are weighted by the production cross section times the pair Branching Ratio times the fraction of c quark fragmenting to D0, D0-, D+/-, or Ds+/-. This fraction is obtained from Pythia and found to be 95%.


#### 2.1.3. Open Beauty

The theoretical uncertainty on open beauty production is relatively smaller than the one on open charm production. We therefore calculated open beauty production yields with a Leading Order (LO) matrix element which was normalised to the Next-To-Leading-Order (NLO)K factor. The NLO cross section with the same setup was calculated by MadGraph5_aMC@NLO [31]. We used CTEQ6L1 (CTEQ6M) for the LO (NLO) calculation. The K factor is found to be 1.83. The renormalisation and factorisation scale is μR=μF=mb2+pT2 with the mass of the b quark taken as mb=4.5 GeV/c2. We have adopted a similar definition for the open beauty correlated background as the one of open charm (see the previous section).


#### 2.1.4. Drell-Yan

Drell-Yan (DY) correlated background was simulated with the processqq-→γ*/Z→μ+μ- at LO where qq- is a pair of the same flavour light quarks. The LO calculation was done with the CTEQ6L1 pdf set and the renormalisation and factorisation scale was set to μR=μF=Q/2. In order to have enough statistics in the J/ψ and ψ(2S) mass window, a phase space cut requesting that the invariant mass of the dimuons (M) is greater than 2.5 GeV/c2 was applied. For the simulation of the DY background under the Υ family peaks, a phase space cut M>7 GeV/c2 was applied. The DY cross section obtained with HELAC-Onia at s=38.8 GeV is compared to the existing E866 data at the same energy [32]. A K factor 1.2 is needed to match the data and therefore it was also applied at s=115 GeV. Such a K factor is known to approximately account for the higher-order QCD corrections.


### 2.2. Uncorrelated Background

The uncorrelated background was obtained from a minimum bias Pythiap+p simulation at s=115 GeV using the processSoftQCD:nonDiffractive with the MRSTMCal.LHgrid LHAPDF (6.1.4) set [34]. By comparing our simulation of open charm with a low statistic pure minimum bias Pythia simulation, we have checked that the contribution of dimuons originating from a muon from charm/beauty and a muon from π/K is negligible. The dominant source of uncorrelated opposite-sign muon pairs is the simultaneous semimuonic decay of uncorrelated π and/or K. In order to avoid possible double counting of signal and correlated background processes, the following hard processes have been switched off from the minimum bias simulations:HardQCD:hardccbar,HardQCD:hardbbbar,WeakSingleBoson:ffbar2gmZ (in order to avoid Drell-Yan pair production.),Charmonium:all andBottomonium:all,.


## 2.1. Signal and Correlated Background

### 2.1.1. Quarkonium Signal

J
/
ψ, ψ(2S), Υ(1S), Υ(2S), and Υ(3S) were simulated in a data-driven way. The amplitude of gg→Q+X (where Q is the quarkonium) is expressed in an empirical functional form [21]: (1)Agg→Q+X2¯=Kexp⁡-κpT2MQ2whenpT≤pTKexp⁡-κpT2MQ21+κnpT2-pT2MQ2-nwhenpT>pT,where K=λ2κs^/MQ2 with s^ the partonic center-of-mass energy and MQ the mass of the quarkonium Q taken from the PDG table [22].The parametersκ, λ, n, and 〈pT〉 were determined by fitting the differential cross section d2σ/dpTdy to the experimental data. The dedicated codes used to perform the fit and to generate unweighted events for quarkonium production have been implemented in HELAC-Onia [18] and we used MSTW2008NLO PDF set [23] provided in LHAPDF5 [24] and the factorisation scale μF=MQ2+pT2. In order to constrain the nontrivial energy dependence of quarkonium production, we used the differential measurements of charmonium production performed by the PHENIX collaboration at RHIC, in p+p collisions at s=200 GeV [9] to predict the corresponding yields at s=115 GeV. Given the lack of such measurements for Υ at RHIC, we performed a combined fit to CDF [25], ATLAS [26], CMS [27], and LHCb [10, 28] data on Υ production. The values of the fitted parameters are listed in Table 3. For illustration, the comparison between fits and the selected experimental data is shown in Figure 1.Table 3
Fit parameters obtained after a combined fit ofd2σ/dpTdy to the PHENIX data [9] for charmonium production and to CDF [25], ATLAS [26], CMS [27], and LHCb [10, 28] data for bottomonium production. We have fixed n=2 and pT= 4.5 (13.5) GeV/c for charmonium (bottomonium) production. The number of fitted data points is also reported.
κ
λ
Number of data points
χ
2
J
/
ψ
0.674
0.380
51
422
ψ
(
2
S
)
0.154
0.351
4
1.12
Υ
(
1
S
)
0.707
0.0837
288
1883
Υ
(
2
S
)
0.604
0.0563
205
856
Υ
(
3
S
)
0.591
0.0411
197
886Figure 1
Some illustrative comparisons between fits and the PHENIX data [9] for charmonium production (a) and LHCb data [10] for bottomonium production (b).
(a)
(b)In order to increase the statistics of the simulated data sample, the decay of the quarkonium in Pythia is forced into the dimuon decay channel. The simulated yields are then weighted by the cross section for this process multiplied by the Branching Ratio (BR).


### 2.1.2. Open Charm

Open charm production was simulated with the processgg→cc- in HELAC-Onia. In order to avoid the huge theoretical uncertainties in the state-of-the-art perturbative calculations, open charm yields at s=115 GeV are also computed in a data-driven way following the method described in the previous section. Similarly, the matrix element of gg→cc- is determined using (1). The parameters are obtained from a fit to the pT-differential cc- cross section measured by the STAR experiment [29] in p+p collisions at s=200 GeV (see Figure 3). We obtained κ=0.437, λ=3.04, and 〈pT〉=2.86 GeV/c when n=2 by using CTEQ6L1 [30] and by fixing the c quark mass to mc=1.5 GeV/c2 and the factorisation scale to μF=mc2+pT2. The χ2 of the fit is equal to 4.39 with 10 experimental data points. The tuned result is shown in Figure 2. The evolution of the cross section with the energy down to s=115 GeV is then given by HELAC-Onia.Figure 2
A comparison between fit and the STAR data [29] in p+p collisions at s=200 GeV, for cc- production.Figure 3
Drell-Yan cross section as a function of the rapidity in the center-of-mass frame obtained with HELAC-Onia + Pythia at s=38.76 GeV and rescaled by a factor 1.2, together with E866 data extracted from [33]. The invariant mass range considered is 7.2<M<8.7 GeV/c2.After embedding the Les Houches Event File into Pythia, muons from the underlying Pythia event can be produced on top of muons from the initial cc- pair. The combination of those additional muons with a muon from the initial cc- pair is not included in our definition of open charm correlated background. We have however checked that this contribution is negligible. In order to increase the statistics, D0, D0-, D+/-, and Ds+/- were forced to decay into muons and only those decay muons were considered as correlated background. μ+μ- pairs coming from all possible combinations, D0D0-, D+D-, Ds+Ds-, D0D+/-, D0Ds+/-, and D+/-Ds-/+, are considered. The simulated events are weighted by the production cross section times the pair Branching Ratio times the fraction of c quark fragmenting to D0, D0-, D+/-, or Ds+/-. This fraction is obtained from Pythia and found to be 95%.


### 2.1.3. Open Beauty

The theoretical uncertainty on open beauty production is relatively smaller than the one on open charm production. We therefore calculated open beauty production yields with a Leading Order (LO) matrix element which was normalised to the Next-To-Leading-Order (NLO)K factor. The NLO cross section with the same setup was calculated by MadGraph5_aMC@NLO [31]. We used CTEQ6L1 (CTEQ6M) for the LO (NLO) calculation. The K factor is found to be 1.83. The renormalisation and factorisation scale is μR=μF=mb2+pT2 with the mass of the b quark taken as mb=4.5 GeV/c2. We have adopted a similar definition for the open beauty correlated background as the one of open charm (see the previous section).


### 2.1.4. Drell-Yan

Drell-Yan (DY) correlated background was simulated with the processqq-→γ*/Z→μ+μ- at LO where qq- is a pair of the same flavour light quarks. The LO calculation was done with the CTEQ6L1 pdf set and the renormalisation and factorisation scale was set to μR=μF=Q/2. In order to have enough statistics in the J/ψ and ψ(2S) mass window, a phase space cut requesting that the invariant mass of the dimuons (M) is greater than 2.5 GeV/c2 was applied. For the simulation of the DY background under the Υ family peaks, a phase space cut M>7 GeV/c2 was applied. The DY cross section obtained with HELAC-Onia at s=38.8 GeV is compared to the existing E866 data at the same energy [32]. A K factor 1.2 is needed to match the data and therefore it was also applied at s=115 GeV. Such a K factor is known to approximately account for the higher-order QCD corrections.


## 2.1.1. Quarkonium Signal

J
/
ψ, ψ(2S), Υ(1S), Υ(2S), and Υ(3S) were simulated in a data-driven way. The amplitude of gg→Q+X (where Q is the quarkonium) is expressed in an empirical functional form [21]: (1)Agg→Q+X2¯=Kexp⁡-κpT2MQ2whenpT≤pTKexp⁡-κpT2MQ21+κnpT2-pT2MQ2-nwhenpT>pT,where K=λ2κs^/MQ2 with s^ the partonic center-of-mass energy and MQ the mass of the quarkonium Q taken from the PDG table [22].The parametersκ, λ, n, and 〈pT〉 were determined by fitting the differential cross section d2σ/dpTdy to the experimental data. The dedicated codes used to perform the fit and to generate unweighted events for quarkonium production have been implemented in HELAC-Onia [18] and we used MSTW2008NLO PDF set [23] provided in LHAPDF5 [24] and the factorisation scale μF=MQ2+pT2. In order to constrain the nontrivial energy dependence of quarkonium production, we used the differential measurements of charmonium production performed by the PHENIX collaboration at RHIC, in p+p collisions at s=200 GeV [9] to predict the corresponding yields at s=115 GeV. Given the lack of such measurements for Υ at RHIC, we performed a combined fit to CDF [25], ATLAS [26], CMS [27], and LHCb [10, 28] data on Υ production. The values of the fitted parameters are listed in Table 3. For illustration, the comparison between fits and the selected experimental data is shown in Figure 1.Table 3
Fit parameters obtained after a combined fit ofd2σ/dpTdy to the PHENIX data [9] for charmonium production and to CDF [25], ATLAS [26], CMS [27], and LHCb [10, 28] data for bottomonium production. We have fixed n=2 and pT= 4.5 (13.5) GeV/c for charmonium (bottomonium) production. The number of fitted data points is also reported.
κ
λ
Number of data points
χ
2
J
/
ψ
0.674
0.380
51
422
ψ
(
2
S
)
0.154
0.351
4
1.12
Υ
(
1
S
)
0.707
0.0837
288
1883
Υ
(
2
S
)
0.604
0.0563
205
856
Υ
(
3
S
)
0.591
0.0411
197
886Figure 1
Some illustrative comparisons between fits and the PHENIX data [9] for charmonium production (a) and LHCb data [10] for bottomonium production (b).
(a)
(b)In order to increase the statistics of the simulated data sample, the decay of the quarkonium in Pythia is forced into the dimuon decay channel. The simulated yields are then weighted by the cross section for this process multiplied by the Branching Ratio (BR).


## 2.1.2. Open Charm

Open charm production was simulated with the processgg→cc- in HELAC-Onia. In order to avoid the huge theoretical uncertainties in the state-of-the-art perturbative calculations, open charm yields at s=115 GeV are also computed in a data-driven way following the method described in the previous section. Similarly, the matrix element of gg→cc- is determined using (1). The parameters are obtained from a fit to the pT-differential cc- cross section measured by the STAR experiment [29] in p+p collisions at s=200 GeV (see Figure 3). We obtained κ=0.437, λ=3.04, and 〈pT〉=2.86 GeV/c when n=2 by using CTEQ6L1 [30] and by fixing the c quark mass to mc=1.5 GeV/c2 and the factorisation scale to μF=mc2+pT2. The χ2 of the fit is equal to 4.39 with 10 experimental data points. The tuned result is shown in Figure 2. The evolution of the cross section with the energy down to s=115 GeV is then given by HELAC-Onia.Figure 2
A comparison between fit and the STAR data [29] in p+p collisions at s=200 GeV, for cc- production.Figure 3
Drell-Yan cross section as a function of the rapidity in the center-of-mass frame obtained with HELAC-Onia + Pythia at s=38.76 GeV and rescaled by a factor 1.2, together with E866 data extracted from [33]. The invariant mass range considered is 7.2<M<8.7 GeV/c2.After embedding the Les Houches Event File into Pythia, muons from the underlying Pythia event can be produced on top of muons from the initial cc- pair. The combination of those additional muons with a muon from the initial cc- pair is not included in our definition of open charm correlated background. We have however checked that this contribution is negligible. In order to increase the statistics, D0, D0-, D+/-, and Ds+/- were forced to decay into muons and only those decay muons were considered as correlated background. μ+μ- pairs coming from all possible combinations, D0D0-, D+D-, Ds+Ds-, D0D+/-, D0Ds+/-, and D+/-Ds-/+, are considered. The simulated events are weighted by the production cross section times the pair Branching Ratio times the fraction of c quark fragmenting to D0, D0-, D+/-, or Ds+/-. This fraction is obtained from Pythia and found to be 95%.


## 2.1.3. Open Beauty

The theoretical uncertainty on open beauty production is relatively smaller than the one on open charm production. We therefore calculated open beauty production yields with a Leading Order (LO) matrix element which was normalised to the Next-To-Leading-Order (NLO)K factor. The NLO cross section with the same setup was calculated by MadGraph5_aMC@NLO [31]. We used CTEQ6L1 (CTEQ6M) for the LO (NLO) calculation. The K factor is found to be 1.83. The renormalisation and factorisation scale is μR=μF=mb2+pT2 with the mass of the b quark taken as mb=4.5 GeV/c2. We have adopted a similar definition for the open beauty correlated background as the one of open charm (see the previous section).


## 2.1.4. Drell-Yan

Drell-Yan (DY) correlated background was simulated with the processqq-→γ*/Z→μ+μ- at LO where qq- is a pair of the same flavour light quarks. The LO calculation was done with the CTEQ6L1 pdf set and the renormalisation and factorisation scale was set to μR=μF=Q/2. In order to have enough statistics in the J/ψ and ψ(2S) mass window, a phase space cut requesting that the invariant mass of the dimuons (M) is greater than 2.5 GeV/c2 was applied. For the simulation of the DY background under the Υ family peaks, a phase space cut M>7 GeV/c2 was applied. The DY cross section obtained with HELAC-Onia at s=38.8 GeV is compared to the existing E866 data at the same energy [32]. A K factor 1.2 is needed to match the data and therefore it was also applied at s=115 GeV. Such a K factor is known to approximately account for the higher-order QCD corrections.


## 2.2. Uncorrelated Background

The uncorrelated background was obtained from a minimum bias Pythiap+p simulation at s=115 GeV using the processSoftQCD:nonDiffractive with the MRSTMCal.LHgrid LHAPDF (6.1.4) set [34]. By comparing our simulation of open charm with a low statistic pure minimum bias Pythia simulation, we have checked that the contribution of dimuons originating from a muon from charm/beauty and a muon from π/K is negligible. The dominant source of uncorrelated opposite-sign muon pairs is the simultaneous semimuonic decay of uncorrelated π and/or K. In order to avoid possible double counting of signal and correlated background processes, the following hard processes have been switched off from the minimum bias simulations:HardQCD:hardccbar,HardQCD:hardbbbar,WeakSingleBoson:ffbar2gmZ (in order to avoid Drell-Yan pair production.),Charmonium:all andBottomonium:all,.


## 3. Fast Simulation of the Response of a LHCb-Like Detector

The HELAC-Onia and Pythia generators provide the opposite-sign muon pairs from quarkonium decays, correlated and uncorrelated backgrounds sources, as defined in the previous section. In order to account for the detector resolution and the particle identification capabilities of a given detector and to investigate the feasibility of the quarkonium studies in p+p collisions at s≃115 GeV, the detector response needs to be simulated. For this purpose, we have chosen a detector setup similar to the LHCb detector [35]. A forward detector is very well suited as a fixed-target experiment setup as well, with a good tracking and particle identification capabilities.According to LHCb analysis cuts, muons in our simulations are required to have their transverse momentum satisfyingpT>0.7 GeV/c [36] and their pseudorapidity in the laboratory frame satisfying 2<η<5. The η cut range corresponds to the LHCb detector coverage. Since the momentum resolution reported by LHCb is δp/p~0.4 (0.6)% for a momentum of 3 (100) GeV/c [37], we consider a momentum resolution of δp/p= 0.5%. The single μ identification efficiency is taken to be ϵP= 98%, which is an average efficiency obtained by LHCb for muons coming from J/ψ decays, for p>3 GeV/c and pT>0.8 GeV/c [37]. These cuts and the abovementioned detector response on the muons are applied to simulate the quarkonium states and all the background sources.In the case of uncorrelated background, as discussed in Section2, most of the μ originate from π+/- or K+/- decays. If a π or K decays to a μ before 12 m along the z axis, the μ is rejected by the tracking system and it is not considered in the simulation. 12 m corresponds to the distance where the calorimeters, followed by the muon stations, are placed in the LHCb detector setup. If the μ is produced beyond 12 m or if a π/K is misidentified with μ in the muon stations, a π/K misidentification probability is applied. The misidentification probabilities depend on the total particle momentum and were reported by the LHCb collaboration in [38]. These probabilities are parameterised with the following functions: PMID(π→μ)(p)=(0.5+6.63exp⁡(-0.13p))% and PMID(K→μ)(p)=(0.5+8.6exp⁡(-0.11p))%, and they are shown in Figures 4(a) and 4(b), for π and K, respectively. Based on the single μ identification efficiency ϵμ+/-, the dimuon, μ+μ-, efficiency is calculated as a product of the single efficiencies: ϵμ+μ-=ϵμ+×ϵμ-. For muons coming from π+ or K+/- decays, misidentification probabilities are used: ϵμ+/-=PMID(π→μ)(p) or ϵμ+/-=PMID(K→μ)(p), respectively, for π and K, and for prompt muons ϵμ+/-=ϵP=0.98.Figure 4
Misidentification probability ofπ (a) and K (b) as muon candidates as a function of momentum, PMID(π→μ) and PMID(K→μ), respectively.
(a)
(b)The pair efficiency is extracted in each kinematic phase space point and is shown as a function of the dimuon invariant mass, transverse momentum, and rapidity in Figure5. This efficiency is used to correct dimuon spectra obtained with the uncorrelated background Pythia simulations.Figure 5
Muon pair,μ+μ-, identification efficiency as a function of the pair invariant mass (a), transverse momentum (b), and rapidity (c) for uncorrelated muon background. The efficiency takes into account the identification efficiency of the prompt muons and the π and K misidentification probability, PMID(π→μ) and PMID(K→μ).
(a)
(b)
(c)


## 4. Quarkonium Production Studies inp+ H Collisions at s=115 GeV

In this section, we show results on the quarkonium production studies in the dimuon decay channels, with the dominant background sources. Simulations have been performed for a 7 TeV proton beam on a hydrogen target (p+p), which gives s=115 GeV. We consider an integrated luminosity of 10 fb−1 which is expected to be obtained after half of a LHC year with the crystal mode, as described in Section 1 and Table 1.


### 4.1. Background Studies

These simulations allow us to quantify the background sources in the quarkonium studies in the dimuon decay channel, which could potentially make the quarkonium signal extraction more difficult or even prevent from obtaining a clear signal. In particular, this may be critical for the excited states. We present here simulations of invariant mass of opposite-sign muon pairs,μ+μ-, from the quarkonia and from the dominant background sources, in two mass ranges; see Figure 6. The first range corresponds to the J/ψ and ψ(2S) invariant mass windows and the second one to the mass range of the Υ(1S), Υ(2S), and Υ(3S). The invariant mass distributions are integrated over the whole transverse momentum and rapidity ranges. The plots show the simulated quarkonium signals and the background, separately from the different sources, and the black solid line is a sum of all contributions. The background sources correspond to an integrated luminosity of 10 fb−1.Figure 6
Dimuon invariant mass distributions forJ/ψ and ψ(2S) (a) and Υ(nS) (b) with different background sources.
(a)
(b)In theJ/ψ and ψ(2S) invariant mass window, the dominant background source is from uncorrelated μ+μ- pairs, mostly from π+/- and K+/- decays. The contributions from Drell-Yan and bb- continuum are very small. In the case of Υ(nS) states, the Drell-Yan contribution is the dominant one. Under the Υ(nS) peak, the contribution from the cc- continuum is negligible, and it is not considered here.The significance (sig=S/(S+B), where S is the number of signal counts and B is the number of background counts, in the invariant mass range MQ±3σQ) and the signal to background ratio (S/B) of each quarkonium state are given in the following:(i)
sigJ/ψ=134.6 102σ, S/BJ/ψ=4.21,(ii)
sigψ(2S)=735.2σ, S/Bψ(2S)=0.16,(iii)
sigΥ(1S)=140.73σ, S/BΥ(1S)=1.75,(iv)
sigΥ(2S)=45.29σ, S/BΥ(2S)=0.48,(v)
sigΥ(3S)=25.75σ, S/BΥ(3S)=0.28,for J/ψ, ψ(2S), Υ(1S), Υ(2S), and Υ(3S), respectively.Transverse momentum and rapidity distributions for the quarkonium signals and for each background source were also studied. As an example, thepT and y distributions in the J/ψ mass range, 3.063<Mμ+μ-<3.129 GeV/c2 (corresponding to MJ/ψ±3σJ/ψ), are shown in Figure 7. It is visible that the distributions for the J/ψ and different backgrounds differ. In more backward or forward rapidity regions, the signal to background ratio increases. This can also be seen in Figure 8, where the dimuon invariant mass distributions in J/ψ and ψ(2S) mass window are shown in three rapidity ranges. In terms of transverse momentum, one can obtain a very clean signal when going to higher pT. Above ~4 GeV/c, the uncorrelated background starts to vanish. Since cc-, bb-, and Drell-Yan simulations are LO simulations, the pT spectra of these correlated background sources are not shown here.Figure 7
p
T (a) and y (b) spectra of J/ψ signal and different background sources in the J/ψ mass range.
(a)
(b)Figure 8
Dimuon invariant mass distributions for three rapidity bins inJ/ψ and ψ(2S) mass window. 2<Yμ+μ-<3, 3<Yμ+μ-<4, and 4<Yμ+μ-<5, shown on panels (a), (b), and (c), respectively.
(a)
(b)
(c)


### 4.2. Quarkonium Simulations

We have also studied thepT and rapidity coverage reach of the quarkonium signals. The transverse momentum distributions are shown in Figure 9(a), for J/ψ, ψ(2S), Υ(1S), Υ(2S), and Υ(3S), from the top to the bottom distribution. Similarly, Figure 9(b) shows the rapidity distribution for each quarkonium state. With an integrated luminosity of 10 fb−1 the quarkonium studies can be carried out in a wide rapidity and pT range. It should be possible to study Υ(nS) signals up to pT≃10 GeV/c, and J/ψ and ψ(2S) could be studied even up to pT≃15 GeV/c. All the quarkonium states can be measured down to pT=0 GeV/c.Figure 9
Transverse momentum (a) and rapidity (b) distributions forJ/ψ, ψ(2S), Υ(1S), Υ(2S), and Υ(3S) from the top to the bottom distribution.
(a)
(b)This study is limited by the rapidity range of2<y<5, in the laboratory frame, due to the pseudorapidity cuts on the decay μ. The red x-axis on the top of Figure 9(b) denotes the rapidity in the center-of-mass frame. The rapidity shift for a 7 TeV proton beam on a fixed-target is −4.8; that is, yCM=0→ylab=4.8. J/ψ and ψ(2S) signals can be studied in the whole mentioned rapidity range, while the lowest rapidity reach for Υ(nS) is ~2.5–3.


## 4.1. Background Studies

These simulations allow us to quantify the background sources in the quarkonium studies in the dimuon decay channel, which could potentially make the quarkonium signal extraction more difficult or even prevent from obtaining a clear signal. In particular, this may be critical for the excited states. We present here simulations of invariant mass of opposite-sign muon pairs,μ+μ-, from the quarkonia and from the dominant background sources, in two mass ranges; see Figure 6. The first range corresponds to the J/ψ and ψ(2S) invariant mass windows and the second one to the mass range of the Υ(1S), Υ(2S), and Υ(3S). The invariant mass distributions are integrated over the whole transverse momentum and rapidity ranges. The plots show the simulated quarkonium signals and the background, separately from the different sources, and the black solid line is a sum of all contributions. The background sources correspond to an integrated luminosity of 10 fb−1.Figure 6
Dimuon invariant mass distributions forJ/ψ and ψ(2S) (a) and Υ(nS) (b) with different background sources.
(a)
(b)In theJ/ψ and ψ(2S) invariant mass window, the dominant background source is from uncorrelated μ+μ- pairs, mostly from π+/- and K+/- decays. The contributions from Drell-Yan and bb- continuum are very small. In the case of Υ(nS) states, the Drell-Yan contribution is the dominant one. Under the Υ(nS) peak, the contribution from the cc- continuum is negligible, and it is not considered here.The significance (sig=S/(S+B), where S is the number of signal counts and B is the number of background counts, in the invariant mass range MQ±3σQ) and the signal to background ratio (S/B) of each quarkonium state are given in the following:(i)
sigJ/ψ=134.6 102σ, S/BJ/ψ=4.21,(ii)
sigψ(2S)=735.2σ, S/Bψ(2S)=0.16,(iii)
sigΥ(1S)=140.73σ, S/BΥ(1S)=1.75,(iv)
sigΥ(2S)=45.29σ, S/BΥ(2S)=0.48,(v)
sigΥ(3S)=25.75σ, S/BΥ(3S)=0.28,for J/ψ, ψ(2S), Υ(1S), Υ(2S), and Υ(3S), respectively.Transverse momentum and rapidity distributions for the quarkonium signals and for each background source were also studied. As an example, thepT and y distributions in the J/ψ mass range, 3.063<Mμ+μ-<3.129 GeV/c2 (corresponding to MJ/ψ±3σJ/ψ), are shown in Figure 7. It is visible that the distributions for the J/ψ and different backgrounds differ. In more backward or forward rapidity regions, the signal to background ratio increases. This can also be seen in Figure 8, where the dimuon invariant mass distributions in J/ψ and ψ(2S) mass window are shown in three rapidity ranges. In terms of transverse momentum, one can obtain a very clean signal when going to higher pT. Above ~4 GeV/c, the uncorrelated background starts to vanish. Since cc-, bb-, and Drell-Yan simulations are LO simulations, the pT spectra of these correlated background sources are not shown here.Figure 7
p
T (a) and y (b) spectra of J/ψ signal and different background sources in the J/ψ mass range.
(a)
(b)Figure 8
Dimuon invariant mass distributions for three rapidity bins inJ/ψ and ψ(2S) mass window. 2<Yμ+μ-<3, 3<Yμ+μ-<4, and 4<Yμ+μ-<5, shown on panels (a), (b), and (c), respectively.
(a)
(b)
(c)


## 4.2. Quarkonium Simulations

We have also studied thepT and rapidity coverage reach of the quarkonium signals. The transverse momentum distributions are shown in Figure 9(a), for J/ψ, ψ(2S), Υ(1S), Υ(2S), and Υ(3S), from the top to the bottom distribution. Similarly, Figure 9(b) shows the rapidity distribution for each quarkonium state. With an integrated luminosity of 10 fb−1 the quarkonium studies can be carried out in a wide rapidity and pT range. It should be possible to study Υ(nS) signals up to pT≃10 GeV/c, and J/ψ and ψ(2S) could be studied even up to pT≃15 GeV/c. All the quarkonium states can be measured down to pT=0 GeV/c.Figure 9
Transverse momentum (a) and rapidity (b) distributions forJ/ψ, ψ(2S), Υ(1S), Υ(2S), and Υ(3S) from the top to the bottom distribution.
(a)
(b)This study is limited by the rapidity range of2<y<5, in the laboratory frame, due to the pseudorapidity cuts on the decay μ. The red x-axis on the top of Figure 9(b) denotes the rapidity in the center-of-mass frame. The rapidity shift for a 7 TeV proton beam on a fixed-target is −4.8; that is, yCM=0→ylab=4.8. J/ψ and ψ(2S) signals can be studied in the whole mentioned rapidity range, while the lowest rapidity reach for Υ(nS) is ~2.5–3.


## 5. Quarkonium Measurements inp+A Collisions at s=115 GeV and Pb + H Collisions at s=72 GeV

### 5.1.  Multiplicity in Proton-Nucleus Collisions

In proton-nucleus collisions, the high track multiplicity may induce a high detector occupancy and lead to a reduction of the detector capabilities. Since LHCb has successfully measured theJ/ψ and Υ production in p + Pb collisions at sNN=5 TeV [39, 40], one would expect a good capability of such detector under similar particle multiplicity environment. In the following, the charged particle multiplicity has been generated with the EPOS generator [41, 42] in different configurations: p+Pb collisions at sNN=5 TeV in collider mode (the lead goes in the detector direction), p+Pb collisions at sNN=115 GeV, and Pb+H collisions at sNN=72 GeV in fixed-target mode. The charged particle multiplicity is dominated by the π multiplicity. By comparing these three distributions as a function of the pseudorapidity of the particle in the laboratory frame as shown in Figure 10, one can conclude that the charged particle multiplicity in a fixed-target mode never exceeds the one obtained in p+Pb collisions at sNN=5 TeV in the collider mode in the full pseudorapidity range: a detector with the LHCb capabilities will be able to run in such conditions.Figure 10
Averaged number of charged particles inp+A collisions as a function of the pseudorapidity in the laboratory frame.


### 5.2. Prospects for the Measurements of the Nuclear Modification Factors forJ/ψ and Υ in p+ Pb Collisions at sNN=115 GeV

To illustrate the potential offered by AFTER@LHC inp+Pb collisions at sNN=115 GeV, we have evaluated, in this section, the impact of the nuclear modification of the gluon densities in nucleons within large nucleus, generically referred to as gluon shadowing, and its uncertainty as encoded in the nuclear PDF set EPS09. For that, we have used the probabilistic Glauber Monte-Carlo framework,JIN [43, 44], which allows us to encode different mechanisms for the partonic production and to interface these production processes with different cold nuclear matter effects, such as the aforementioned shadowing, in order to get the production cross sections for proton-nucleus and nucleus-nucleus collisions.JIN also straightforwardly computes any nuclear modification factor for minimum bias collisions or in specific centrality classes. In the case of proton-nucleus (p+A) collisions, it is the ratio of the yield per inelastic collision in p+A collisions to the yield in pp collisions at the same energy multiplied by the average number of binary collisions in a typical p+p collision, 〈Ncoll〉: (2)RpA=dNpANcolldNpp.In the presence of a net nuclear effect, RpA is defined such that it differs fromunity. In the simplest case of minimum bias collisions, one should have (3)RpA=dσpAAdσpp.As in [45], we have used the central curve of EPS09 as well as four specific extreme curves (minimal/maximal shadowing, minimal/maximal EMC effect), which reproduce the envelope of the gluon nPDF uncertainty encoded in EPS09 LO [46].In addition to the modification of the partonic densities, quarkonium production inp+A collisions can be affected by other effects, for instance, by the nuclear absorption which depends much on the nature of the object traversing the nuclear medium. If the meson is already formed, it may be affected more than a smaller preresonant pair. To discuss such an effect, it is useful to introduce the concept of the formation time, tf, based on the Heisenberg uncertainty principle and the time, in the rest frame of the meson, to discriminate between two S states, for instance, the J/ψ and the ψ(2S). In fact, one finds [45, 47] that such a time is similar for the charmonium and bottomonium states and is on the order of 0.3–0.4 fm. Obviously, this time has to be boosted in the frame where the nuclear matter sits. For tf smaller than the nucleus radius, the quarkonium is formed before escaping it. In the fixed-target mode with a proton beam and a nuclear target, the boost factor is simply γ(ylab)=cosh(ylab). We therefore obtain tf as in Table 4.Table 4
Boost and formation time in the (target) Pb rest frame ofJ/ψ and  Υ as a function of its  c.m.s. rapidity at sNN=115 GeV.
y
c
.
m
.
s
y
l
a
b
γ
(
y
l
a
b
)
t
f
J
/
ψ
,
Υ
(
y
)
y
c
.
m
.
s
y
l
a
b
γ
(
y
l
a
b
)
t
f
J
/
ψ
,
Υ
(
y
)
−2.5
2.3
5
1.75 fm
−0.5
4.3
37
13 fm
−1.5
3.3
14
5 fm
0.0
4.8
61
21 fm
−1.0
3.8
22
8 fm
0.5
5.3
100
35 fmOne sees that looking at quarkonium production inp+Pb collisions at different backward rapidities allows one to look at quarkonia traversing the nuclear matter at very different stages of their evolution. This effect could theoretically be studied by giving an ad hoc rapidity dependence to the effective absorption cross section, σabseffective. This is left for a future study since, here, we wish to consider only the nPDF effects and the expected statistics. Other effects to be considered are the coherent energy loss [48] (expected to grow in the forward region) and the rescattering by comovers [49] (expected to grow with the multiplicity along the J/ψ direction).Since we wish to assess the discriminating power of the possible data to be taken with AFTER@LHC, we attribute to the EPS09 central values statistical uncertainties which directly follow from the differential yields, respectively, expected inp+p and p+Pb collisions. For that we take an integrated luminosity of 10 fb−1 for the p+p runs and 100 pb−1 for the p+Pb runs, in accordance with the luminosities discussed above (see Table 1). As this stage, we do not consider additional systematical uncertainties. This simplifying assumption could be lifted in a more detailed study which would also take into account a possible detector acceptance (and related efficiencies) as done in the previous section. In particular, we do not expect that the rapidity region for yCMS>1.5 would be easily accessible.In Figure11, we show the rapidity dependence of Rp+Pb for Υ and its pT dependence near y=0. The only million of Υ to be collected per year allows for the measurement of a Rp+Pb with a much better precision than the gluon nPDF, nearly up to x→1. In addition, one notes that the nuclear modification factor is certainly measurable up to pT≃10 GeV/c.Figure 11
Nuclear modification factor forϒ as a function of (a) yCMS and (b) pT in p+Pb collisions at sNN=115 GeV. The uncertainties attached to the central points are derived from the statistics to be collected with Lp+p=10 fb−1 and Lp+Pb=100 pb−1.
(a)
(b)In Figure12, we also show the rapidity dependence of Rp+Pb for J/ψ and its pT dependence near y=0. In both cases, the luminosity to be taken in a year at AFTER@LHC yields to statistical uncertainties which are largely negligible as compared to the nPDF uncertainties; the statistical uncertainties are not even visible in Figure 12. We expect this to hold also for the ψ(2S) although its yields are down by a factor of 100.Figure 12
Nuclear modification factor forJ/ψ as a function of (a) yCMS and (b) pT in p+Pb collisions at sNN=115 GeV. The uncertainties derived from the statistics to be collected with Lp+p=10 fb−1 and Lp+Pb=100 pb−1 are smaller than the point size.
(a)
(b)As aforementioned, the nPDFs do not account for all the expected nuclear matter effects. However, it is clear that combining the measurements ofΥ, J/ψ, and ψ(2S) for -3<yCMS<0 (as a LHCb-like detector would do) will allow one to pin down the existence of a possible gluon EMC and antishadowing effect. We also stress that the complications induced by a rapidity dependence of σabseffective could be avoided by the parallel measurement of Rp+Pb for nonprompt J/ψ which can only be sensitive to the energy loss since the b quark decay (weakly) into the J/ψ, way outside the nucleus. Figure 13 shows that the trend is similar compared to Υ. Measuring the pT dependence of Rp+Pb for prompt J/ψ and Υ should also avoid the sensitivity on formation time effects.Figure 13
Nuclear modification factor forJ/ψ from b as a function of yCMS in p+Pb collisions at sNN=115 GeV. The uncertainties derived from the statistics to be collected with Lp+p=10 fb−1 and Lp+Pb=100 pb−1 are smaller than the point size.


## 5.1.  Multiplicity in Proton-Nucleus Collisions

In proton-nucleus collisions, the high track multiplicity may induce a high detector occupancy and lead to a reduction of the detector capabilities. Since LHCb has successfully measured theJ/ψ and Υ production in p + Pb collisions at sNN=5 TeV [39, 40], one would expect a good capability of such detector under similar particle multiplicity environment. In the following, the charged particle multiplicity has been generated with the EPOS generator [41, 42] in different configurations: p+Pb collisions at sNN=5 TeV in collider mode (the lead goes in the detector direction), p+Pb collisions at sNN=115 GeV, and Pb+H collisions at sNN=72 GeV in fixed-target mode. The charged particle multiplicity is dominated by the π multiplicity. By comparing these three distributions as a function of the pseudorapidity of the particle in the laboratory frame as shown in Figure 10, one can conclude that the charged particle multiplicity in a fixed-target mode never exceeds the one obtained in p+Pb collisions at sNN=5 TeV in the collider mode in the full pseudorapidity range: a detector with the LHCb capabilities will be able to run in such conditions.Figure 10
Averaged number of charged particles inp+A collisions as a function of the pseudorapidity in the laboratory frame.


## 5.2. Prospects for the Measurements of the Nuclear Modification Factors forJ/ψ and Υ in p+ Pb Collisions at sNN=115 GeV

To illustrate the potential offered by AFTER@LHC inp+Pb collisions at sNN=115 GeV, we have evaluated, in this section, the impact of the nuclear modification of the gluon densities in nucleons within large nucleus, generically referred to as gluon shadowing, and its uncertainty as encoded in the nuclear PDF set EPS09. For that, we have used the probabilistic Glauber Monte-Carlo framework,JIN [43, 44], which allows us to encode different mechanisms for the partonic production and to interface these production processes with different cold nuclear matter effects, such as the aforementioned shadowing, in order to get the production cross sections for proton-nucleus and nucleus-nucleus collisions.JIN also straightforwardly computes any nuclear modification factor for minimum bias collisions or in specific centrality classes. In the case of proton-nucleus (p+A) collisions, it is the ratio of the yield per inelastic collision in p+A collisions to the yield in pp collisions at the same energy multiplied by the average number of binary collisions in a typical p+p collision, 〈Ncoll〉: (2)RpA=dNpANcolldNpp.In the presence of a net nuclear effect, RpA is defined such that it differs fromunity. In the simplest case of minimum bias collisions, one should have (3)RpA=dσpAAdσpp.As in [45], we have used the central curve of EPS09 as well as four specific extreme curves (minimal/maximal shadowing, minimal/maximal EMC effect), which reproduce the envelope of the gluon nPDF uncertainty encoded in EPS09 LO [46].In addition to the modification of the partonic densities, quarkonium production inp+A collisions can be affected by other effects, for instance, by the nuclear absorption which depends much on the nature of the object traversing the nuclear medium. If the meson is already formed, it may be affected more than a smaller preresonant pair. To discuss such an effect, it is useful to introduce the concept of the formation time, tf, based on the Heisenberg uncertainty principle and the time, in the rest frame of the meson, to discriminate between two S states, for instance, the J/ψ and the ψ(2S). In fact, one finds [45, 47] that such a time is similar for the charmonium and bottomonium states and is on the order of 0.3–0.4 fm. Obviously, this time has to be boosted in the frame where the nuclear matter sits. For tf smaller than the nucleus radius, the quarkonium is formed before escaping it. In the fixed-target mode with a proton beam and a nuclear target, the boost factor is simply γ(ylab)=cosh(ylab). We therefore obtain tf as in Table 4.Table 4
Boost and formation time in the (target) Pb rest frame ofJ/ψ and  Υ as a function of its  c.m.s. rapidity at sNN=115 GeV.
y
c
.
m
.
s
y
l
a
b
γ
(
y
l
a
b
)
t
f
J
/
ψ
,
Υ
(
y
)
y
c
.
m
.
s
y
l
a
b
γ
(
y
l
a
b
)
t
f
J
/
ψ
,
Υ
(
y
)
−2.5
2.3
5
1.75 fm
−0.5
4.3
37
13 fm
−1.5
3.3
14
5 fm
0.0
4.8
61
21 fm
−1.0
3.8
22
8 fm
0.5
5.3
100
35 fmOne sees that looking at quarkonium production inp+Pb collisions at different backward rapidities allows one to look at quarkonia traversing the nuclear matter at very different stages of their evolution. This effect could theoretically be studied by giving an ad hoc rapidity dependence to the effective absorption cross section, σabseffective. This is left for a future study since, here, we wish to consider only the nPDF effects and the expected statistics. Other effects to be considered are the coherent energy loss [48] (expected to grow in the forward region) and the rescattering by comovers [49] (expected to grow with the multiplicity along the J/ψ direction).Since we wish to assess the discriminating power of the possible data to be taken with AFTER@LHC, we attribute to the EPS09 central values statistical uncertainties which directly follow from the differential yields, respectively, expected inp+p and p+Pb collisions. For that we take an integrated luminosity of 10 fb−1 for the p+p runs and 100 pb−1 for the p+Pb runs, in accordance with the luminosities discussed above (see Table 1). As this stage, we do not consider additional systematical uncertainties. This simplifying assumption could be lifted in a more detailed study which would also take into account a possible detector acceptance (and related efficiencies) as done in the previous section. In particular, we do not expect that the rapidity region for yCMS>1.5 would be easily accessible.In Figure11, we show the rapidity dependence of Rp+Pb for Υ and its pT dependence near y=0. The only million of Υ to be collected per year allows for the measurement of a Rp+Pb with a much better precision than the gluon nPDF, nearly up to x→1. In addition, one notes that the nuclear modification factor is certainly measurable up to pT≃10 GeV/c.Figure 11
Nuclear modification factor forϒ as a function of (a) yCMS and (b) pT in p+Pb collisions at sNN=115 GeV. The uncertainties attached to the central points are derived from the statistics to be collected with Lp+p=10 fb−1 and Lp+Pb=100 pb−1.
(a)
(b)In Figure12, we also show the rapidity dependence of Rp+Pb for J/ψ and its pT dependence near y=0. In both cases, the luminosity to be taken in a year at AFTER@LHC yields to statistical uncertainties which are largely negligible as compared to the nPDF uncertainties; the statistical uncertainties are not even visible in Figure 12. We expect this to hold also for the ψ(2S) although its yields are down by a factor of 100.Figure 12
Nuclear modification factor forJ/ψ as a function of (a) yCMS and (b) pT in p+Pb collisions at sNN=115 GeV. The uncertainties derived from the statistics to be collected with Lp+p=10 fb−1 and Lp+Pb=100 pb−1 are smaller than the point size.
(a)
(b)As aforementioned, the nPDFs do not account for all the expected nuclear matter effects. However, it is clear that combining the measurements ofΥ, J/ψ, and ψ(2S) for -3<yCMS<0 (as a LHCb-like detector would do) will allow one to pin down the existence of a possible gluon EMC and antishadowing effect. We also stress that the complications induced by a rapidity dependence of σabseffective could be avoided by the parallel measurement of Rp+Pb for nonprompt J/ψ which can only be sensitive to the energy loss since the b quark decay (weakly) into the J/ψ, way outside the nucleus. Figure 13 shows that the trend is similar compared to Υ. Measuring the pT dependence of Rp+Pb for prompt J/ψ and Υ should also avoid the sensitivity on formation time effects.Figure 13
Nuclear modification factor forJ/ψ from b as a function of yCMS in p+Pb collisions at sNN=115 GeV. The uncertainties derived from the statistics to be collected with Lp+p=10 fb−1 and Lp+Pb=100 pb−1 are smaller than the point size.


## 6. Prospects of Pb+A Measurements at s=72 GeV

The charged particle multiplicity has been generated with the EPOS generator [41, 42] in different configurations: Pb+Pb at sNN=5.5 TeV in collider mode, Pb+Ar, Pb+Xe, and Pb+Pb at sNN=72 GeV in fixed-target mode. By comparing these three distributions in the pseudorapidity of the particle in the laboratory frame as shown in Figure 14, one can conclude that the charged particle multiplicity in a fixed-target mode never exceeds the one obtained in Pb+Pb collisions at sNN=5.5 TeV obtained in a collider mode in the full pseudorapidity range: a detector with the ALICE MFT+Muon detector [50] capability will be able to run in such conditions. Detailed studies are needed to evaluate up to which multiplicity a detector such as LHCb would be able to take good quality data.Figure 14
Averaged number of charged particles inA+A collisions as a function of the pseudorapidity in the laboratory frame.


## 7. Conclusion

In summary, we have shown that in a fixed-target mode with an integrated luminosity of 10 fb−1, using 7 TeV LHC proton beam on a hydrogen target, and with a detector setup and performances similar to the LHCb detector, quarkonium studies in the dimuon decay channel can be performed over a wide transverse momentum range and rapidity in the center-of-mass from ~−2.8 for J/ψ and ψ(2S) and ~−2 for Υ states to ~0. We have performed simulations of the dominant background sources contributing to the μ+μ- invariant mass spectrum. The uncorrelated background was obtained using Pythia generator and dimuons from correlated background sources, cc-, bb-, and Drell-Yan, were simulated using both HELAC-Onia and Pythia generators. The estimated background level allows for J/ψ, ψ(2S), Υ(1S), Υ(2S), and Υ(3S) measurements in the dimuon decay channel with good signal to background ratios.These simulations set the stage for further ones including, on the one hand, the detection of photon fromP wave or ηc decay or from the production of a J/ψ+γ pair, whose studies at low transverse momentum can provide important insight on the gluon transverse dynamics [51–54], and, on the other hand, the large combinatorial background typical of p+A and A+A collisions in which the study of excited quarkonium at AFTER@LHC energies is of paramount importance [1, 13]. We note that the Delphes [55] framework seems particularly well suited to account for the photon detectability in such prospective studies.Along our investigations, we have also noted that the main source of dimuons around theΥ(nS) masses is from the Drell-Yan process (see Figure 6(b)). This gives us great confidence that the corresponding cross section can easily be extracted in this mass region in p+p collisions, a fortiori with a vertex detector allowing for tagging the heavy-flavour muons. We therefore consider that the single-spin asymmetries for Drell-Yan pair production can indeed be extracted using a light polarised target. Motivations for such studies are discussed in [15, 16, 56]. Quarkonium polarisation measurements are of course also possible, given the large statistical samples.As regards the case ofp+A collisions, we have had a first look at the charged particle multiplicities as a function of the laboratory pseudorapidity. We have found out that, for all the possible fixed-target modes, p + Pb, Pb + H, and Pb + Pb, these are smaller than the ones reached in the collider modes, where the LHCb was used (p + Pb and Pb + p at 5 TeV). We therefore believe that a detector with similar characteristics as compared to LHCb can very well be used in the fixed-target mode (our observation is obviously supported by the preliminary analysis of the LHCb-SMOG data taken during the pilot run of p+ beam (Pb beam) on a Neon gas target from 2012 (2013) at a c.m.s energy of sNN=87 GeV (54 GeV) [57]).In view of the above, we have evaluated the impact, and its uncertainty, on the nuclear modification of the gluon densities on prompt and nonpromptJ/ψ and Υ in the form of Rp+Pb. We have found that the measurements at backward rapidities allow one to search for the gluon antishadowing, the gluon EMC effect, and even the Fermi motion effect on the gluons with the unheard statistical precisions. The statistics are large enough to perform such measurement with the ψ(2S) and probably also with Υ(2S) and Υ(3S) allowing for thorough investigations of the formation time effect of the meson propagating in the nuclear matter. Overall, our results confirm the great potential of AFTER@LHC for heavy-quark and quarkonium physics.



---

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