ALICE Collaboration / Physics Letters B 798 (2019) 134926 3
Fig. 1. Left: invariant mass distribution for muon pairs satisfying the event selection described in the text. The dashed green line corresponds to the background. The solid
magenta and red lines correspond to Crystal Ball functions representing J/ψ and ψ
signals, respectively. The solid blue line corresponds to the sum of background and signal
functions. Right: transverse momentum distribution for muon pairs in the range 2.85 < m
μμ
< 3.35 GeV/c
2
(around the J/ψ mass).
−4.0 < y < −2.5 and in Fig. 2 in six rapidity subranges. The in-
variant
mass distributions are fitted with a function modeling the
background and two Crystal Ball functions [26]for the J/ψ and the
ψ
peaks. The shape of the background at large invariant masses is
well described by an exponential distribution, as expected if it is
dominated by the process γγ → μ
+
μ
−
. However, at masses be-
low
the J/ψ , the distribution is strongly influenced by the muon
trigger condition. In order to model this, the whole background
distribution is fitted using a template made from reconstructed
STARlight events corresponding to the γγ → μ
+
μ
−
process. The
results of the fit are parametrized using a fourth-order polynomial,
which turns smoothly into an exponential tail as from 4GeV/c
2
.
The coefficients of the polynomial are then kept fixed in the fit to
the experimental data, while the slope of the exponential term and
the normalization are left free. The fitted slope is found to agree
within 2.5 standard deviations with the value obtained from the
generated sample.
The
raw inclusive J/ψ and ψ
yields, N(J/ψ) and N(ψ
), were
obtained by fitting the dimuon invariant mass spectrum in the
range 2.2 < m
μμ
< 6GeV/c
2
. The slope parameters in the Crys-
tal
Ball functions were fixed from fits to the respective Monte
Carlo sets. The width parameter σ
J/ψ
was left free for the J/ψ ,
and was fixed to σ
ψ
= σ
J/ψ
· (σ
MC
ψ
/σ
MC
J
/ψ
) for the ψ
, where the
ratio σ
MC
ψ
/σ
MC
J
/ψ
∼ 1.09 of the ψ
to the J/ψ widths was obtained
from the fits to corresponding Monte Carlo sets. The mass parame-
ter
of the Crystal Ball function was left unconstrained for the J/ψ.
Due to the small ψ
statistics, the ψ
mass was fixed so that the
difference with respect to the J/ψ mass is the same as quoted by
the PDG [27]. The J/ψ mass m
J/ψ
= 3.0993 ± 0.0009 GeV/c
2
, ob-
tained
from the fit in the full rapidity range −4.0 < y < −2.5, is
in agreement with the PDG value within 3 standard deviations.
The
raw inclusive J/ψ yields obtained from invariant mass fits
contain contributions from the coherent and incoherent J/ψ pho-
toproduction,
which can be separated in the analysis of transverse
momentum spectra. The p
T
distributions for dimuons in the range
2.85 < m
μμ
< 3.35 GeV/c
2
are shown in Fig. 1, right, in the full
dimuon rapidity range −4.0 < y < −2.5 and in Fig. 3 in six ra-
pidity
subranges. These distributions were fitted with Monte Carlo
templates produced using STARlight, corresponding to different
production mechanisms: coherent J/ψ , incoherent J/ψ, feed-down
J/ψ from coherent ψ
decays, feed-down J/ψ from incoherent ψ
decays and continuum dimuons from the γγ →μ
+
μ
−
process. In
order to describe the high-p
T
tail, the incoherent J/ψ photopro-
duction
accompanied by nucleon dissociation was also taken into
account in the fits with the template based on the H1 parametriza-
tion
of the dissociative J/ψ photoproduction [28](denoted as dis-
sociative
J/ψ in the following):
dN
dp
T
∼ p
T
1 +
b
pd
n
pd
p
2
T
−n
pd
. (1)
The H1 collaboration provided two sets of measurements corre-
sponding
to different photon–proton center-of-mass energy ranges:
25 GeV < W
γ p
< 80 GeV (low-energy data set) and 40 GeV <
W
γ p
< 110 GeV (high-energy data set). The fit parameters b
pd
=
1.79 ±0.12 (GeV/c)
−2
and n
pd
= 3.58 ±0.15 from the high-energy
data set were used by default, while the corresponding uncer-
tainties
and the low-energy values b
pd
= 1.6 ± 0.2 (GeV/c)
−2
and
n
pd
= 3.58(fixed) were used for systematic checks.
The templates were fitted to the data leaving the normaliza-
tion
free for coherent J/ψ , incoherent J/ψ and dissociative J/ψ
production. The normalization of the γγ → μ
+
μ
−
spectrum was
fixed to the one obtained from the invariant mass fits. The nor-
malization
of the coherent and incoherent feed-down J/ψ tem-
plates
was constrained to the normalization of primary coherent
and incoherent J/ψ templates, according to the feed-down frac-
tions
extracted from the measurement of raw inclusive J/ψ and
ψ
yields, as described below. The extracted incoherent J/ψ frac-
tion
f
I
=
N(incoh J/ψ)
N(coh J/ψ)
for p
T
< 0.25 GeV/c ranges from (4.9 ± 0.6)%
to (6.4 ± 0.8)% depending on the rapidity interval and is consis-
tent
with being constant within the uncertainties of the fits. The
contribution of incoherent J/ψ with nucleon dissociation was also
taken into account in this fraction.
4. Results and discussion
4.1. Ratio of coherent ψ
and J/ψ cross sections
The obtained dimuon invariant mass spectra can be used to ex-
tract
the ratio of coherent ψ
and J/ψ cross sections R =
σ (ψ
)
σ (J/ψ )
and the fraction of feed-down J/ψ from ψ
decays in the raw J/ψ
yields. The fits to the invariant mass distributions for dimuons with
pair p
T
< 0.25 GeV/c in the full rapidity range −4.0 < y < −2.5
result
in the following ratio of the measured raw inclusive ψ
and
J/ψ yields: