没有合适的资源?快使用搜索试试~ 我知道了~
温馨提示
ϒ(1S),ϒ(2S)和ϒ(3S)量子态的产量通过在CMS检测器中通过衰变成μ对来测量,在每个核子对质心能量的PbPb和pp碰撞中 2.76 TeV。 该数据分别对应于PbPb和pp碰撞的166μb-1和5.4pb-1的综合亮度。 差异生产横截面据报道是速度y高达2.4的函数,横向动量pT高达20GeV / c的函数。 相对于pp碰撞,在PbPb中观察到了强烈的依赖于中心点的抑制作用,分别对于1(1S)和ϒ(2S)状态,其受阻系数高达≈2和8。 没有观察到这种抑制作为y或pT的函数的显着依赖性。 在PbPb碰撞中未观察到ϒ(3S)状态,这对应于在95%置信度下对中心积分数据至少抑制了≈7倍。 所观察到的抑制与对夸克胶子等离子体中的夸克态连续熔化建模的理论场景相一致。
资源推荐
资源详情
资源评论
Physics Letters B 770 (2017) 357–379
Contents lists available at ScienceDirect
Physics Letters B
www.elsevier.com/locate/physletb
Suppression of ϒ(1S), ϒ(2S), and ϒ(3S) quarkonium states in PbPb
collisions at
√
s
NN
= 2.76 TeV
.The CMS Collaboration
CERN, Switzerland
a r t i c l e i n f o a b s t r a c t
Article history:
Received
4 November 2016
Received
in revised form 26 March 2017
Accepted
13 April 2017
Available
online 19 April 2017
Editor:
M. Doser
Keywords:
CMS
Physics
Heavy
ions
Dimuons
Bottomonia
The production yields of ϒ(1S), ϒ(2S), and ϒ(3S) quarkonium states are measured through their decays
into muon pairs in the CMS detector, in PbPb and pp collisions at the centre-of-mass energy per nucleon
pair of 2.76 TeV. The data correspond to integrated luminosities of 166 μb
−1
and 5.4 pb
−1
for PbPb and
pp collisions, respectively. Differential production cross sections are reported as functions of ϒ rapidity
y up to 2.4, and transverse momentum p
T
up to 20 GeV/c. A strong centrality-dependent suppression is
observed in PbPb relative to pp collisions, by factors of up to ≈2and 8, for the ϒ(1S) and ϒ (2S) states,
respectively. No significant dependence of this suppression is observed as a function of y or p
T
. The ϒ (3S)
state is not observed in PbPb collisions, which corresponds to a suppression for the centrality-integrated
data by at least a factor of ≈7at a 95% confidence level. The observed suppression is in agreement with
theoretical scenarios modeling the sequential melting of quarkonium states in a quark gluon plasma.
© 2017 The Author(s). Published by Elsevier B.V. This is an open access article under the CC BY license
(http://creativecommons.org/licenses/by/4.0/). Funded by SCOAP
3
.
1. Introduction
At large energy density and high temperature, strongly interact-
ing
matter is predicted by lattice QCD calculations to consist of a
deconfined system of quarks and gluons [1]. This state, often re-
ferred
to as “quark gluon plasma” (QGP) [2], constitutes the main
object of studies using high energy heavy ion collisions.
The formation of QGP in nuclear collisions is studied in a va-
riety
of ways. One of its most striking signatures is the sequential
suppression of quarkonium states, both in the charmonium (J/ψ,
ψ
, χ
c
, etc.) and the bottomonium (ϒ(1S, 2S, 3S), χ
b
, etc.) fami-
lies.
Historically, this phenomenon was proposed as direct evidence
of deconfinement because, in the deconfined medium, the binding
potential between the constituents of a quarkonium state, a heavy
quark and its antiquark (QQ), should be screened by the colour
charges of the surrounding light quarks and gluons [3,4]. The sup-
pression
of quarkonium production is predicted to occur above the
critical temperature of the medium (T
c
) and to depend on the
QQbinding energy. Since the ϒ(1S) is the most tightly bound state
among all quarkonia, it is expected to have the highest dissocia-
tion
temperature. Estimates of dissociation temperatures are given
in Ref. [5]: T
dissoc
≈ 2 T
c
, 1.2 T
c
, and 1 T
c
for the ϒ(1S), ϒ(2S),
and ϒ(3S) states, respectively. Other medium effects, such as re-
generation
from initially uncorrelated quark–antiquark pairs [6,7]
E-mail address: cms-publication-committee-chair@cern.ch.
or absorption by comoving particles [8,9] can modify quarkonium
production in heavy ion collisions. Furthermore, nuclear effects
such as modifications of parton distributions inside nuclei [10]
or
energy loss processes in nuclear matter [11] are expected to
affect the production of quarkonia independently of any QGP for-
mation.
An admixture of several of the above-mentioned effects in
the context of bottomonium production is investigated in Refs. [12,
13] and
a recent review on quarkonium production can be found
in Ref. [14].
The suppression of ϒ(1S) production in heavy ion collisions rel-
ative
to pp yields scaled by the number of binary nucleon–nucleon
(NN) collisions was first measured by CMS [15] in the midrapidity
range |y| < 2.4, then by ALICE at forward rapidities 2.5 < y < 4
[16].
Both measurements were done at the CERN LHC in PbPb
collisions at a centre-of-mass energy per nucleon pair,
√
s
NN
, of
2.76 TeV. Alarger suppression of the ϒ(2S) and ϒ(3S) was first
suggested [17] then observed [18] by CMS. In pPb collisions at
√
s
NN
= 5.02 TeV, ALICE [19] and LHCb [20] reported ϒ(1S) yields
that are slightly suppressed along the p-going forward direction,
possibly indicating the importance of nuclear effects. Lacking pp
reference data at
√
s
NN
= 5.02 TeV, the pp yields were estimated
by interpolating results at 2.76, 7, and 8TeV [19], or by scaling
data at 8TeV [20]. The ϒ(2S) and ϒ(3S) were reported by CMS
to be slightly more suppressed than the ϒ(1S) ground state in
pPb collisions [21]. At the BNL RHIC, STAR reported no significant
suppression of the overlapping ϒ(1S +2S +3S) states in dAu colli-
sions
at
√
s
NN
= 200 GeV, while observing a suppression in central
http://dx.doi.org/10.1016/j.physletb.2017.04.031
0370-2693/
© 2017 The Author(s). Published by Elsevier B.V. This is an open access article under the CC BY license (http://creativecommons.org/licenses/by/4.0/). Funded by
SCOAP
3
.
358 The CMS Collaboration / Physics Letters B 770 (2017) 357–379
AuAu collisions at the same energy [22]. Altogether, these results
are interpreted as a sequential suppression of the three states in
nucleus–nucleus collisions [12,13], with the tighter bound states
disappearing less in the QGP.
This Letter reports the production yields of ϒ(1S), ϒ(2S), and
ϒ(3S) for PbPb and pp data at the same
√
s
NN
= 2.76 TeV, using
integrated luminosities of 166 μb
−1
and 5.4 pb
−1
, respectively. The
two sets of data correspond to approximately the same number of
NN collisions. The pp sample collected in 2013 contains 20 times
more events than the 2011 data used previously [15,17,18], allow-
ing
further differential studies with respect to the ϒ meson rapid-
ity
and transverse momentum. Muon reconstruction is improved in
PbPb collisions relative to Ref. [18], yielding a 35% increase in the
number of measured ϒ candidates. In total, the improved recon-
struction
and a relaxed muon-p
T
selection provide almost twice
the number of ϒ(1S) candidates used in Ref. [18]. The yields in
PbPb and pp events are used to extract nuclear modification fac-
tors,
R
AA
.
2. The CMS detector
A detailed description of the CMS detector, together with a def-
inition
of the coordinate system used and the relevant kinematic
variables, can be found in Ref. [23]. The central feature of the CMS
apparatus is a superconducting solenoid of 6m internal diameter.
A silicon tracker, a crystal electromagnetic calorimeter, and a brass
and scintillator hadron calorimeter reside within the magnetic field
volume.
Muons
are detected in the pseudorapidity interval |η| < 2.4
using
gas-ionization detectors made of three technologies: drift
tubes, cathode strip chambers, and resistive-plate chambers, em-
bedded
in the steel flux-return yoke of the solenoid. The silicon
tracker is composed of pixel detectors (three barrel layers and two
forward disks on either side of the detector, made of 66 million
100×150 μm
2
pixels) followed by microstrip detectors (ten barrel
layers, and three inner and nine forward disks on either side of the
detector, with strips of pitch between 80 and 180 μm). The trans-
verse
momentum of muons matched to tracks reconstructed in the
silicon detector is measured with a resolution better than 1.5% for
p
T
values smaller than 100 GeV/c [24]. This high resolution is the
result of the 3.8 T magnetic field and the high granularity of the
silicon tracker.
In
addition, CMS has extensive forward calorimetry, includ-
ing
two steel and quartz-fibre Cherenkov hadron forward (HF)
calorimeters, that cover the range 2.9 < |η| < 5.2. These detectors
are used in the present analysis to select events and to determine
the centrality of PbPb collisions, as described in the next section.
3. Data selections
3.1. Event selection and centrality
To select purely inelastic hadronic PbPb collisions, contributions
from ultraperipheral collisions and noncollision beam backgrounds
are removed, as described in Ref. [25]. Events are preselected if
they contain a primary vertex built from at least two tracks, and
at least three signals (one in the case of pp collisions) in HF tow-
ers
on each side of the interaction point with deposited energies of
at least 3GeV in each tower. To further suppress beam-gas events,
the distribution of hits in the pixel detector along the beam direc-
tion
is required to be compatible with particles originating from
the event vertex. These criteria select (97 ± 3)%of the inelas-
tic
hadronic PbPb collisions [25], yielding an efficiency-corrected
number of minimum bias (MB) events N
MB
= (1.16 ± 0.04) × 10
9
for the MB sample corresponding to this analysis. The pp data
Table 1
Average
values of the number of participating nucleons (N
part
, with the root-mean-
square
of its distribution in each bin), and nuclear overlap function (T
AA
, with its
systematic uncertainty) for the centrality bins used in the ϒ(1S) (upper rows) and
ϒ(2S) (middle) analyses. The centrality-integrated values are given in the last row.
Centrality (%) N
part
(RMS) T
AA
(mb
−1
)
0–5 381 (19) 25.9 ±1.1
5–10 329 (22) 20.5
±0.9
10–20 261 (30) 14.5
±0.8
20–30 187 (23) 8.80
±0.58
30–40 130 (18) 5.09
±0.43
40–50 86.2 (13.6) 2.75
±0.30
50–70 42.0 (14.4) 0.98
±0.14
70–100 8.8 (6.0) 0.125
±0.023
0–10 355 (33) 23.2
±1.0
10–30 224 (46) 11.6
±0.7
30–50 108 (27) 3.92
±0.37
50–100 22.1 (19.3) 0.47
±0.07
0–100 113 (115) 5.67
±0.32
correspond to an integrated luminosity of 5.4 pb
−1
, known to an
accuracy of 3.7% coming from the uncertainty in the calibration
based on a van der Meer scan [26].
The
measurements are based on events that were first selected
by the Level-1 trigger, a hardware-based system that uses informa-
tion
from the muon detectors and calorimeters. The presence of at
least two muons was required, with no selection applied on their
momenta. The events were then further filtered using a software-
based
high-level trigger, and rejected if muons were poorly recon-
structed,
hence likely to be misidentified. The pp and PbPb data
were collected using the same trigger logic.
The
centrality of PbPb collisions is defined as the fraction of
the total number of inelastic hadronic collisions, with 0% repre-
senting
collisions with the largest overlap of the two nuclei. This
fraction is determined from the distribution of total energy in both
HF calorimeters. Variables related to the centrality, such as the
number of nucleons participating in the collision (N
part
) and the
nuclear overlap function (T
AA
) [27], are estimated using a Glauber
model simulation described in Ref. [25]. The value of T
AA
at a
given centrality is equal to the number of binary NN collisions
divided by the NN cross section and can be interpreted as the NN-
equivalent
integrated luminosity per heavy ion collision.
It
is to be noted that the PbPb hadronic cross section (7.65 ±
0.42 b) computed with this Glauber simulation corresponds to an
integrated luminosity of 152 ±9 μb
−1
, compatible within 1.2 sigma
with the experimental value of 166 ± 8 μb
−1
based on the van der
Meer scan. The mean values of T
AA
and N
part
are presented in Ta-
ble 1 for
the narrow centrality bins used in the ϒ(1S) analysis,
the wider bins used in the ϒ(2S) analysis, and the centrality-
integrated
estimate. The most peripheral bins are rather wide
and, since quarkonium yields scale with the number of nucleon–
nucleon
collisions, most bottomonia are produced close to the
most central edge of the bins, namely 70% and 50%. The N
part
values shown in the following figures and reported in Table 1 are
computed by averaging over all MB events in a given centrality bin,
and are therefore not corrected for any bias introduced by requir-
ing
the presence of the ϒ . Also presented is the root-mean-square
(RMS) of the N
part
distribution in each bin. The uncertainty on T
AA
is computed by varying the Glauber parameters and the event se-
lection
inefficiency, as described in Ref. [25]. In this Letter, N
part
is used to show the centrality dependence of the measurements,
while T
AA
directly enters into the nuclear modification factor cal-
culation:
R
AA
= N
PbPb
/
T
AA
σ
pp
where N
PbPb
is the number of ϒ
produced per PbPb collision in a given kinematic range and σ
pp
the corresponding ϒ cross section in pp collisions.
The CMS Collaboration / Physics Letters B 770 (2017) 357–379 359
3.2. Muon selection
Muons are reconstructed using a global fit to a track in the
muon detectors that is matched to a track in the silicon tracker.
The offline muon reconstruction algorithm used for the PbPb data
has been improved relative to that used previously [18]. The
efficiency has been increased by running multiple iterations in
the pattern recognition step, raising the number of reconstructed
ϒ(1S) candidates by approximately 35%. Background muons from
cosmic rays and heavy-quark semileptonic decays are rejected by
imposing a set of selection criteria on each muon track. These
criteria are based on previous studies of the performance of the
muon reconstruction algorithm [28]. The track is required to have
a hit in at least one pixel detector layer, and a respective transverse
(longitudinal) distance of closest approach of less than 3 (15) cm
from the measured primary vertex, primarily to reject cosmic ray
muons and muons from hadron decays in flight. To ensure a good
p
T
measurement, more than 10 hits are requested in the tracker,
and the χ
2
per number of degrees of freedom of the trajectory
fits is limited to be smaller than 10 when using the silicon tracker
and the muon detectors, and smaller than 4 when using only the
tracker. Pairs of oppositely charged muons are considered when
the χ
2
fit probability of the tracks originating from a common ver-
tex
exceeds 1%.
For
the ϒ(2S) and ϒ(3S) analyses, the transverse momentum
of each muon (p
μ
T
) is required to be above 4 GeV/c, as in previ-
ous
publications [15,17,18], while one of them is relaxed down to
3.5 GeV/c for the ϒ(1S) analysis. Reducing this p
T
threshold raises
the ϒ(1S) yield by approximately 40%, and its statistical signifi-
cance
by up to 50%, depending on the p
T
and y of the dimuon
system. Relaxing the criterion on the second muon was also con-
sidered
then discarded, since it did not significantly raise the ac-
ceptance
for the ϒ states. The resulting invariant mass distribu-
tions
are shown on Fig. 1 for the entire pp and PbPb data samples.
4. Analysis
4.1. Signal extraction
To extract the ϒ(1S), ϒ(2S), and ϒ(3S) meson yields, unbinned
maximum likelihood fits to the μ
+
μ
−
invariant mass spectra are
performed between 7.5 and 14 GeV/c
2
. The results for the p
T
-, y-
and
centrality-integrated case are displayed as solid lines on Fig. 1.
Each ϒ resonance is modelled by the sum of two Crystal Ball (CB)
functions [29] with common mean but different widths to account
for the pseudorapidity dependence of the muon momentum res-
olution.
The CB functions are Gaussian resolution functions with
the low-side tail replaced by a power law describing final-state ra-
diation.
This choice was guided by simulation studies, as well as
analyses of large pp event samples collected at
√
s = 7 TeV [30].
Given the relatively large statistical uncertainties, the only signal
model parameters that are left free in the fit are the mean of the
ϒ(1S) peak, and the ϒ(1S), ϒ(2S) and ϒ(3S) meson yields. The
other parameters, such as the width of the ϒ(1S) peak are fixed
in every bin to the corresponding value obtained from simulations.
The mean and width of the CB functions describing the ϒ(2S) and
ϒ(3S) peaks are set by the fitted ϒ(1S) peak mean and the fixed
ϒ(1S) width, respectively, multiplied by the world-average mass
ratio [31]. The parameters describing the tail of the CB function
are fixed to values obtained from simulations, kept common in
the three ϒ states, then allowed to vary when computing the as-
sociated
systematic uncertainties. The background distribution is
modelled by an exponential function multiplied by an error func-
tion
(the integral of a Gaussian) describing the low-mass turn-on,
with all parameters left free in the fit.
Fig. 1. Dimuon invariant mass distributions in pp (top) and centrality-integrated
PbPb (bottom) data at
√
s
NN
= 2.76 TeV, for muon pairs having one p
T
greater than
4 GeV/c and the other greater than 3.5 GeV/c. The solid (signal + background) and
dashed (background only) lines show the result of fits described in the text.
With one muon having p
T
greater than 4 GeV/c and the other
greater than 3.5 GeV/c, this fitting procedure results in ϒ(1S) me-
son
yields and statistical uncertainties of 2534 ±76 and 5014 ±87
in
centrality-integrated PbPb and pp collisions, respectively. With
both muons’ transverse momenta above 4 GeV/c, it yields 173 ±41
for
ϒ(2S) and 7 ±38 for ϒ(3S) (hence unobserved) in PbPb colli-
sions,
and 1214 ± 51 for ϒ(2S) and 618 ± 44 for ϒ(3S) states in
pp collisions.
4.2. Acceptance and efficiency
To correct yields for acceptance and efficiency in the two
data samples, the three ϒ states have been simulated using the
pythia 6.412
generator [32] and embedded in PbPb events sim-
ulated
with hydjet 1.8 [33], producing Monte Carlo (MC) events
with the same settings as in Ref. [18], including radiative tails
handled by photos [34]. Acceptance is defined as the fraction of
ϒ in the |y| < 2.4range that decay into two muons, each with
|η
μ
| < 2.4, and p
μ
2
T
> 4 GeV/c and p
μ
1
T
> 3.5or 4 GeV/c for the
ϒ(1S) and ϒ(2S)/ϒ(3S) states, respectively. For the ϒ(1S) state,
the acceptance over the analyzed phase space averages to 35%. For
all three ϒ states, the acceptance is constant over most of the ra-
pidity
range, with a drop at large |y|. When the ϒ meson has
360 The CMS Collaboration / Physics Letters B 770 (2017) 357–379
p
T
≈ 5 GeV/c, the lower p
T
decay muon often falls below the re-
quired
momentum to reach the muon detector, resulting in a drop
in acceptance for intermediate p
T
. For ϒ(2S) and ϒ(3S) states,
where p
T
for both muons is required to be above 4 GeV/c, the ac-
ceptance
is 28 and 33%, respectively. Within this acceptance, the
average reconstruction and trigger efficiencies are 68, 74 and 75%
for the ϒ(1S), ϒ(2S), and ϒ(3S) states, respectively. The slightly
lower efficiency for the ϒ(1S) state arises from including lower-p
T
muons, which have smaller reconstruction efficiencies, in particu-
lar
at midrapidity.
The
individual components of the efficiency are crosschecked
using collision data and muons from J/ψ meson decays, with a
technique called tag-and-probe, similar to the one described in
Ref. [30]. The method consists of fitting the J/ψ candidates in
data and MC samples, with and without applying the probed se-
lection
criterion on one of the muons. The muon reconstruction,
identification, and trigger efficiencies in the muon detectors are
probed by testing the selection response in a sample collected
with single-muon triggers. The small discrepancies observed be-
tween
the results for data and simulation are used to determine
p
T
- and η-dependent single-muon correction factors that are ap-
plied
to muons in the simulation. The net correction factors to the
ϒ meson yields range from 3 to 18%, the largest being located
at low p
T
or at large |y|. The tracker efficiency, larger than 99%,
is also evaluated with this method by checking the presence of a
track for muons that are primarily reconstructed in the muon de-
tectors.
The corresponding uncertainty is evaluated to be 0.3 and
0.6% for each muon, for the pp and PbPb data, respectively.
4.3. Systematic uncertainties
The uncertainty from the fitting procedure is estimated by per-
forming
seven changes in the fitting functions. Five of them consist
of releasing one by one the originally fixed signal-shape parame-
ters,
to accommodate for possible imperfections in the simulation.
The other two changes consist of adding to the default background
function a first- or second-order Chebychev polynomial. The max-
ima
of the five signal and of the two background variations are
summed in quadrature, yielding systematic uncertainties from 4 to
25% in the PbPb data and from 1 to 10% in the pp data, for the
ϒ(1S) meson yield. For the less significant ϒ(2S) signal, the un-
certainties
range from 13 to 71% in PbPb, and from 1 to 15% in pp
data.
The
systematic uncertainty from the acceptance and efficiency
estimation includes changes of the generated p
T
and y spectra, as
well as variations of the distribution of ϒ candidates across event
centrality, within limits imposed by the data themselves. These are
propagated into bin-by-bin systematic uncertainties of 0.7 and 1.1%
on average, in pp and PbPb collisions, respectively.
Single-muon
efficiencies obtained from the tag-and-probe
method
are assigned a systematic uncertainty from varying re-
quirements
for the tag selection, the dimuon mass range, and
the distributions of the invariant mass peak and the underlying
backgrounds. The maximum deviation in each p
μ
T
and η
μ
inter-
val
is retained as the systematic uncertainty on the single-muon
correction factors. Next, the single-muon correction factors are
changed within their statistical uncertainties derived from data. To
do so, one hundred variations of the single-muon efficiencies are
computed, resulting in one hundred dimuon efficiency correction
factors in each analysis bin. The RMS of the resulting efficien-
cies,
summed in quadrature with the systematic uncertainty in the
efficiency correction factors, represent the overall uncertainty in
muon efficiency. The resulting systematic uncertainties range from
3.2 to 7.7% from midrapidity in pp collisions to the most forward
bins in PbPb collisions. In addition, the uncertainty in the tracking
Fig. 2. Differential cross section for ϒ states as a function of their transverse mo-
mentum
and per unit of rapidity in pp (top) and PbPb (bottom) collisions. The PbPb
results are integrated over centrality and divided by the number of elementary NN
collisions. Statistical (systematic) uncertainties are displayed as error bars (boxes).
Global relative uncertainties of 3.7% (pp) and 6.5% (PbPb) are not displayed.
efficiency of 0.3 and 0.6% for each track is considered as fully cor-
related
and thus doubled for dimuon candidates, and taken as a
global uncertainty (common to all points).
The
relative uncertainties in the integrated luminosity of pp
data (3.7%) or the number of PbPb MB events (3%) are also con-
sidered
as global uncertainties. The uncertainties in the T
AA
values
are given in Table 1.
5. Results
5.1. Cross sections
Figs. 2 and 3 show the differential cross sections as functions
of p
T
(per unit of rapidity) and |y|, respectively, in pp (top) and
PbPb (bottom) collisions. Measured yields are corrected for the ac-
ceptance
and efficiency, then divided by the width of the bin in
consideration. To put the pp and PbPb data on a comparable scale,
the corrected yields are normalized by the measured integrated
luminosity in pp collisions, and by the product of the number of
corresponding MB events and the centrality-integrated T
AA
value
for PbPb collisions. The statistical uncertainties in pp collisions al-
low
a measurement for the three states using the same binning:
five bins in p
T
with edges at 0, 2.5, 5.0, 8.0, 12.0, and 20.0 GeV/c,
剩余22页未读,继续阅读
资源评论
weixin_38737630
- 粉丝: 1
- 资源: 929
上传资源 快速赚钱
- 我的内容管理 展开
- 我的资源 快来上传第一个资源
- 我的收益 登录查看自己的收益
- 我的积分 登录查看自己的积分
- 我的C币 登录后查看C币余额
- 我的收藏
- 我的下载
- 下载帮助
最新资源
资源上传下载、课程学习等过程中有任何疑问或建议,欢迎提出宝贵意见哦~我们会及时处理!
点击此处反馈
安全验证
文档复制为VIP权益,开通VIP直接复制
信息提交成功