214 ALICE Collaboration / Physics Letters B 766 (2017) 212–224
Fig. 1. (Colour online.) Invariant mass distributions of US dimuons with 2.5 < y < 4and p
T
< 12 GeV/c. The top (bottom) row shows the distribution before (after) background
subtraction with the event-mixing technique. The left panels correspond to the most central events (0–10%) while the right panels to a peripheral (70–80%) centrality range.
The fit curves shown in blue represent the sum of the signal and background shapes, while the red lines correspond to the J/ψ signal and the grey ones to the background.
dimuon continuum component is still present and was fitted us-
ing
the sum of two exponentials. Several fitting sub-ranges, within
the interval 2 < m
μμ
< 5GeV/c
2
, were used for both signal ex-
traction
procedures.
Fig. 1 sho
ws examples of fits to the US dimuon invariant mass
distributions with and without background subtraction using the
event-mixing technique, for different selections in centrality. The
raw J/ψ yield in each centrality or p
T
interval was determined
as the average of the results obtained with the two fitting ap-
proaches,
the various parameterisations of signal and background
and the different fitting ranges, while the corresponding systematic
uncertainties were defined as the RMS of these results. A further
contribution to the systematic uncertainty was estimated by using
a different set of resonance tails obtained using in the MC simu-
lation
a different particle transport model (GEANT4 [37] instead of
GEANT3 [38]). The total number of J/ψ , integrated over centrality,
p
T
and y, is N
J/ψ
= 2.77 ±0.02(stat) ±0.05(syst) ·10
5
. The system-
atic
uncertainty ranges from 1.6% to 2.8% as a function of centrality
and from 1.2% to 3.1% as a function of p
T
.
The
nuclear modification factor, as a function of the centrality
class i of the collision and for the J/ψ transverse-momentum in-
terval
p
T
, is calculated as
R
i
AA
(p
T
) =
N
i
J
/ψ
(p
T
)
BR
J/ψ→μ
+
μ
−
N
i
MB
Aε
i
(p
T
)T
i
AA
σ
pp
J
/ψ
(p
T
)
,
(1)
where N
i
J
/ψ
(p
T
) is the number of extracted J/ψ in a given cen-
trality
and p
T
range, BR
J/ψ→μ
+
μ
−
= 5.96 ± 0.03% is the branching
ratio of the dimuon decay channel [39], N
i
MB
is the number of
equivalent minimum-bias events, Aε
i
(p
T
) is the product of the
detector acceptance times the reconstruction efficiency, T
i
AA
is
the average of the nuclear overlap function, and σ
pp
J
/ψ
(p
T
) is the
inclusive J/ψ cross section for pp collisions at the same energy and
in the same kinematic range as the Pb–Pb data.
The
Aε values were determined from MC simulations, with the
generated p
T
and y distributions for the J/ψ adjusted on data,
and separately tuned for each centrality class using an iterative
approach. Unpolarised J/ψ production was assumed [21]. For the
tracking chambers, the time-dependent status of each electronic
channel during the data taking period was taken into account as
well as the misalignment of the detection elements. The efficien-
cies
of the muon trigger chambers were determined from data and
were then applied in the simulations. Finally, the dependence of
the efficiency on the detector occupancy was taken into account
by embedding MC-generated J/ψ into real minimum-bias Pb–Pb
events.
For
J/ψ produced within 2.5 < y < 4 and p
T
< 12 GeV/c, in
0–90% most central collisions, the Aε value is 0.136 ± 0.007(syst).
A relative decrease of the efficiency by 14% was observed when
going from peripheral to central collisions. As a function of p
T
,
Aε has a minimum value of about 0.12 at p
T
≈ 1.5GeV/c, and
then steadily increases up to about 0.4 at the upper end of the
considered range. The following sources of systematic uncertainty
on Aε were considered. A first contribution of 2% due to the in-
put
MC p
T
and y distributions was estimated by (i) varying the
input shapes that were tuned on data within their statistical un-
certainties
and (ii) taking into account the effect of possible p
T
− y
correlations
by comparing, as a function of centrality, the Aε val-
ues
with the corresponding result of a 2-D acceptance calculation
in classes of p
T
and y. A second contribution comes from the
tracking efficiency and it was estimated by comparing the single-
muon
tracking efficiency values obtained, in MC and data, with
a procedure that exploits the redundancy of the tracking-chamber
information [21]. A 3% systematic uncertainty on the dimuon track-
ing
efficiency is obtained and is approximately constant as a func-
tion
of centrality and kinematics. The systematic uncertainty on
the dimuon trigger efficiency represents the third contribution and
it has two origins: the intrinsic efficiencies of the muon trigger
chambers and the response of the trigger algorithm. The first one
was determined from the uncertainties on the trigger chamber ef-
ficiencies
measured from data and applied to simulations and it
amounts to 1.5%. The second one was estimated by comparing the
p
T
dependence, at the single-muon level, of the trigger response
function between data and MC and it varies between 0.2% and 4.6%
as a function of p
T
. Combining the two sources, a systematic un-
certainty
ranging from 1.5% to 4.8% as a function of the J/ψ p
T
is
obtained. Finally, there is a 1% contribution related to the choice