First measurement of quarkonium polarization in nuclear collisions at the LHC

This is a self-archived version of an original article. This version may differ from the original in pagination and typographic details. Author(s): Title: Year: Version: Copyright: Rights: Rights url: Please cite the original version: CC BY 4.0 https://creativecommons.org/licenses/by/4.0/ First measurement of quarkonium polarization in nuclear collisions at the LHC © 2021 the Authors Published version ALICE Collaboration ALICE Collaboration. (2021). First measurement of quarkonium polarization in nuclear collisions at the LHC. Physics Letters B, 815, Article 136146. https://doi.org/10.1016/j.physletb.2021.136146 2021 Physics Letters B 815 (2021) 136146 Contents lists available at ScienceDirect Physics Letters B www.elsevier.com/locate/physletb First measurement of quarkonium polarization in nuclear collisions at the LHC .ALICE Collaboration a r t i c l e i n f o a b s t r a c t Article history: Received 9 June 2020 Received in revised form 10 February 2021 Accepted 12 February 2021 Available online 22 February 2021 Editor: M. Doser The polarization of inclusive J/ψ and ϒ(1S)produced in Pb–Pb collisions at √sNN =5.02 TeV at the LHC is measured with the ALICE detector. The study is carried out by reconstructing the quarkonium through its decay to muon pairs in the rapidity region 2.5 <y <4and measuring the polar and azimuthal angular distributions of the muons. The polarization parameters λθ, λφand λθφ are measured in the helicity and Collins-Soper reference frames, in the transverse momentum interval 2 <pT<10 GeV/cand pT< 15 GeV/cfor the J/ψ and ϒ(1S), respectively. The polarization parameters for the J/ψ are found to be compatible with zero, within a maximum of about two standard deviations at low pT, for both reference frames and over the whole pTrange. The values are compared with the corresponding results obtained for pp collisions at √s=7and 8 TeV in a similar kinematic region by the ALICE and LHCb experiments. Although with much larger uncertainties, the polarization parameters for ϒ(1S)production in Pb–Pb collisions are also consistent with zero. ©2021 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 SCOAP3. 1. Introduction Quarkonia, bound states of charm (c) and anticharm (c) or bottom (b) and antibottom (b) quarks, represent an important tool to test our understanding of quantum chromodynamics (QCD), since their production process involves both perturbative and nonperturbative aspects. At high energy, the creation of the heavy quark-antiquark pair is a process that can be described using a perturbative QCD approach, due to the large value of the charm and bottom quark masses (mc∼1.3GeV/c2, mb∼4.2GeV/c2[1]). However, the subsequent formation of the bound state is a nonperturbative process that can be described only by empirical models or effective field theory approaches. Among those, models based on Non-Relativistic QCD (NRQCD) [2]give the most successful description of the production cross section, as measured at high-energy hadron colliders (Tevatron, RHIC, LHC) [3–14]. In the NRQCD approach, the non-perturbative aspects are parameterized via long-distance matrix elements (LDME), corresponding to the possible intermediate color, spin and angular momentum states of the evolving quark-antiquark pair. The values of LDMEs need to be fitted on a subset of the available measurements and can be then considered as universal quantities, in the sense that they can be used in the calculation of production cross sections and other observables corresponding, for example, to different collision systems and energies. Other theory approaches, as the Color E-mail address: alice -publications @cern .ch. Singlet Model [15], the Color Evaporation Model [16] and the kTfactorization [17]are also used to describe the quarkonium production process. Among the various charmonium states, the J/ψ meson, with quantum numbers JPC =1−−, was the first to be discovered. It is surely the most studied, also due to the sizeable decay branching ratio to dilepton pairs ((5.961 ±0.033)% for the μ+μ−channel [1]) that represents an excellent experimental signature. While the J/ψ production cross sections are well reproduced by NRQCD-based models, it was soon realized that describing the measured polarization of this state represents a much more difficult problem [18]. The polarization, corresponding to the orientation of the particle spin with respect to a chosen axis, can be accessed via a study of the polar (θ) and azimuthal (φ) production angles, relative to that axis, of the two-body decay products in the quarkonium rest frame. Their angular distribution W(θ, φ) is parameterized as W(θ, φ) ∝1 3+λθ ×1+λθcos2θ+λφsin2θcos2φ+λθφ sin2θcosφ,(1) with the polarization parameters λθ, λφand λθφ corresponding to various combinations of the elements of the spin density matrix of J/ψproduction [19]. In particular, the two cases (λθ=1, λφ=0, λθφ =0) and (λθ=−1, λφ=0, λθφ =0) correspond to the so-called transverse and longitudinal polarizations, respectively. At https://doi.org/10.1016/j.physletb.2021.136146 0370-2693/©2021 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 SCOAP3. ALICE Collaboration Physics Letters B 815 (2021) 136146 leading order, the high-pTproduction is dominated by gluon fragmentation and therefore the J/ψwould be expected to be transversely polarized [18]. However, the results from the CDF experiment at Tevatron showed that the J/ψexhibits a very small polarization [20,21], an observation which was impossible to reconcile with the NRQCD prediction. As of today, on the experimental side, accurate results on inclusive and prompt (i.e., removing contributions from b-quark decays) J/ψpolarization have become available at LHC energies [22–25]. They confirm that this state shows little or no polarization in a wide rapidity (up to y =4.5) and transverse momentum region (from 2 to 70 GeV/c), with the exception of the LHCb measurements at √s=7TeV[24], where the value λθ=−0.145 ±0.027, corresponding to a weak longitudinal polarization, was obtained in the interval 2 <pT<15 GeV/cand 2 <y <4.5, in the helicity frame (its definition will be given later in Sec. 3). On the theory side, a huge effort was pursued in order to move to a complete next-to-leading order (NLO) description of the J/ψproduction process [26,27], and to the calculation of the polarization variables [28,29]. Further important progress includes a quantitative evaluation of the contribution of feed-down processes (J/ψcoming from the decay of χcand ψ(2S)states) on the polarization observables [30]. It was shown that at NLO there are rather large cancellations between contributions corresponding to the different possible combinations of the spin and angular momentum of the intermediate ccstates, reaching a more satisfactory description of the absence of polarization observed in the data [31]. However, those descriptions usually require the inclusion of both cross section and polarization results in the fit of the LDME, leading to a more limited predictive power on the polarization observables and to large variations in the values of the extracted LDME values, depending on the set of data used for their determination. Finally, the description of the J/ψproduction in the NRQCD framework was recently extended to the low-pTregion, and the polarization parameters were studied in a color glass condensate (CGC) + NRQCD formalism, obtaining a fair agreement with LHC data at forward rapidity [32]. Measurements of the polarization parameters are also available for several bottomonium states, and in particular for the ϒ(1S), ϒ(2S)and ϒ(3S)resonances, which were shown to exhibit little or no polarization at LHC energies [33–35]. Approaches similar to that adopted for charmonium, which also need to take into account the rather complex feed-down decay structure for these states, lead to a fair agreement with the experimental results [36]. In this Letter, we move a step forward by presenting the first measurement of J/ψand ϒ(1S)polarization in ultrarelativistic heavy-ion interactions performed by the ALICE Collaboration by studying Pb–Pb collisions at √sNN =5.02 TeV. Such collisions represent an important source of information for the investigation of the phase diagram of QCD [37], and in particular for the study of the properties of the quark–gluon plasma (QGP), a state of matter where quarks and gluons are not confined inside hadrons [38]. Among the experimental observables studied in heavy-ion collisions the suppression of heavy quarkonium production is a fundamental signal, since QGP formation prevents the binding of the heavy-quark pair due to the screening of the color charge [39] and, more generally, has strong effects on the spectral functions [40]. At LHC energies, another mechanism, corresponding to the (re)generation of charmonium states in the QGP and/or when the system hadronizes, becomes relevant [41,42], in particular at low pT, due to the large charm-quark multiplicity (>100 pairs in a central Pb–Pb collision). The presence of a deconfined system may in principle affect also the polarization of quarkonium states. In Ref. [43]the observation of a partial transverse polarization for the J/ψwas predicted in case of QGP formation, due to a modification of the non-perturbative effects in the high energy-density phase. More generally, the observed prompt J/ψare known to be a mixture of direct production and decay products from higher-mass charmonium states (ψ(2S), χc). In nuclear collisions, since suppression effects are expected to affect more strongly the less bound states, the relative contribution of direct and feed-down production would change with respect to that in pp collisions, and the overall measured polarization may be different according to the potentially different polarization of the various states [44,45]. On the other hand, the contribution of the regeneration mechanism in the J/ψformation process by recombination of uncorrelated ccpairs is likely to give rise to unpolarized production at low pT. Finally, the possible presence of polarization is known to strongly affect the acceptance for J/ψdetection in the dilepton decay (up to 20–30% in ALICE [22]), and its measurement is an important requisite for an unbiased evaluation of the absolute yields in nuclear collisions. A first measurement of ϒ(1S)polarization in Pb–Pb collisions is also presented in this Letter, even if the corresponding candidate sample is smaller by a factor ∼30, leading to larger uncertainties. For such a state, considerations similar to those discussed for the J/ψshould hold, except that the contribution of the regeneration mechanism should be negligible due to the much lower multiplicity of bottom quarks with respect to charm. The next sections of the Letter are organized as follows. Section 2contains a short description of the experimental apparatus and some details on the data sample used in this analysis. The analysis procedure and the evaluation of systematic uncertainties are presented in Sec. 3, while the results on the J/ψand ϒ(1S) polarization parameters λθ, λφand λθφ are shown in Sec. 4. The conclusions are presented in Sec. 5. 2. Experimental setup and data sample The measurement described in this Letter is performed with the ALICE detector [46,47], whose main components are a central barrel and a forward muon spectrometer. The latter covers the pseudorapidity region −4 <η<−2.5 and is used to detect muon pairs from quarkonium decays [48]. The muon spectrometer includes a hadron absorber made of concrete, carbon and steel with a thickness of 10 interaction lengths, followed by five tracking stations (cathode-pad chambers), with the central one embedded inside a dipole magnet with a 3 T·m field integral. Downstream of the tracking system, an iron wall filters out the remaining hadrons as well as low-momentum muons originating from pion and kaon decays, and is followed by two trigger stations (resistive plate chambers). Another forward detector, the V0 [49], composed of two scintillator arrays located at opposite sides of the interaction point (IP) and covering the pseudorapidity intervals −3.7 <η<−1.7 and 2.8 <η<5.1, provides the minimum bias (MB) trigger which is given by a coincidence of signals from the two sides. Among the central barrel detectors, the two layers of the Silicon Pixel Detector (SPD), with |η| <2 and |η| <1.4 coverage, and corresponding to the inner part of the ALICE Inner Tracking System (ITS) [50], are used to determine the position of the interaction vertex. Finally, the Zero Degree Calorimeters (ZDC) [51], located on either side of the IP at ±112.5 m along the beam axis, detect spectator nucleons emitted at zero degrees with respect to the LHC beam axis and are used to reject electromagnetic Pb–Pb interactions. The analysis is based on events where, in addition to the MB condition, two opposite-sign tracks are detected in the triggering system of the muon spectrometer (dimuon trigger). The dimuon trigger selects tracks each having a transverse momentum above a threshold nominally set at pμ T=1GeV/c, corresponding to the value for which the single-muon trigger efficiency reaches 50% [52]. The single-muon trigger efficiency reaches a plateau value of 98% at ∼2.5GeV/c. The events are further characterized according to their centrality, i.e., the degree of geometric overlap of the colliding nuclei. It is estimated by means of a Glauber model fit to the V0 signal am2 ALICE Collaboration Physics Letters B 815 (2021) 136146 Fig. 1. Examples of fits to the invariant-mass distributions in the helicity reference frame. The left plot corresponds to the J/ψmass region, while on the right a fit to the ϒ(1S)mass region is shown. The fits are performed using an extended Crystal Ball function for the resonance signals, while the background is parameterized with a variable width Gaussian. plitude distribution [53,54], with more central events leading to a larger signal in the V0. In this analysis, events corresponding to the most central 90% of the inelastic Pb–Pb cross section are selected, as for these events the MB trigger is fully efficient and the residual contamination from electromagnetic processes is negligible. The results of the analysis are obtained using the √sNN = 5.02 TeV Pb–Pb data samples collected by the ALICE experiment during the years 2015 and 2018, corresponding to an integrated luminosity Lint ∼750 μb−1. 3. Data analysis The J/ψand ϒ(1S)candidates are formed by combining opposite-sign muons reconstructed using the tracking algorithm described in Ref. [48]. In order to reject tracks at the edge of the spectrometer acceptance, the condition −4 <ημ<−2.5is required. In addition, tracks must have a radial transverse position at the end of the absorber in the range 17.6 <Rabs <88.9cm. This selection is applied to remove tracks passing through the inner and denser part of the absorber, which are strongly affected by multiple scattering. For each muon candidate, a match between tracks reconstructed in the tracking system and track segments in the muon trigger system is required. The J/ψpolarization parameters λθ, λφand λθφ are studied as a function of transverse momentum in the intervals 2 <pT<4, 4 <pT<6 and 6 <pT<10 GeV/c. For each pTinterval, a two– dimensional (2D) grid of dimuon invariant-mass spectra is created, corresponding to intervals in cosθand φ, where θand φ are the polar and azimuthal emission angles, respectively, of the decay products in the J/ψrest frame, with respect to the reference axis. More in detail, the 2D grid covers the fiducial region −0.8 <cosθ<0.8(17 intervals), 0.5 <φ<π−0.5rad (8 intervals, assuming a symmetric distribution around φ=π), with the choice of the boundaries as well as the width of the intervals dictated by acceptance considerations. The analysis is performed choosing two different reference systems for the determination of the angular variables. In the CollinsSoper (CS) frame the z-axis is defined as the bisector of the angle between the direction of one beam and the opposite of the direction of the other one in the rest frame of the decaying particle, allowing therefore an evaluation of the polarization parameters with respect to the direction of motion of the colliding hadrons. In the helicity (HE) reference frame the z-axis is given by the direction of the decaying particle in the center-of-mass frame of the collision, and therefore the polarization can be evaluated with respect to the momentum direction of the J/ψitself. The φ=0 plane is the one containing the two beams in the J/ψrest frame. For each dimuon invariant-mass spectrum, the J/ψraw yield is obtained by means of a binned maximum likelihood fit in the interval 2.1 <mμμ <4.9GeV/ c2. The background continuum is parameterized with a Gaussian distribution whose width varies linearly with the mass or, alternatively, with a fourth degree polynomial function times an exponential. The J/ψ signal is modeled with a pseudo-Gaussian function or with a Crystal Ball function with asymmetric tails on both sides of the peak [55]. The J/ψ mass is kept free in the fits, while for each interval (i, j) in (cosθ, φ) the width is fixed to σi,j J/ψ =σi,j,MC J/ψ ·(σJ/ψ /σMC J/ψ ), i.e., scaling the resonance width extracted from Monte Carlo (MC) simulations (σi,j,MC J/ψ ) by the ratio between the width obtained by fitting the angle-integrated spectrum in data (σJ/ψ ) and MC (σMC J/ψ ) for the pTinterval under consideration. The parameters of the non-Gaussian tails of the resonance are kept fixed to the MC values. The ψ(2S)contribution, although comparatively negligible, is also taken into account in the fits, with its width and mass fixed in each fit to those of the J/ψaccording to the relations σψ(2S)=σJ/ψ ·σMC ψ(2S)/σMC J/ψ and mψ(2S)=mJ/ψ +mPDG ψ(2S)−mPDG J/ψ , with the Particle Data Group (PDG) masses taken from Ref. [1]. In Fig. 1 (left) an example of a fit to the invariant-mass spectrum in the J/ψ mass region is shown. Due to the stability of the extracted J/ψparameters (mass, width), the fits were carried out directly on the sum of the 2015 and 2018 invariant mass spectra. The J/ψ raw yields as a function of the angular variables are then corrected by the product of the acceptance and detector efficiency (A ×ε), which is evaluated as a function of cosθand φ on a 2D grid via MC simulations. The J/ψ are generated according to pTand ydistributions directly tuned on data [56]via an iterative procedure [57], and their decay muons are propagated inside a realistic description of the ALICE setup, based on GEANT 3.21 [58]. The misalignment of the detection elements and the time-dependent status of each electronic channel during the data taking period are taken into account as well. In the J/ψ generation an isotropic distribution of decay products, corresponding to the assumption of no polarization, is adopted. Due to the choice of relatively small (cosθ, φ) intervals, the A ×εvalues for each interval are quite insensitive to the specific angular distribution assumed in the generation. The three polarization parameters λθ, λφand λθφ are obtained through χ2-minimization fits of the 2D J/ψ distributions, cor3 ALICE Collaboration Physics Letters B 815 (2021) 136146 Fig. 2. Fit to the J/ψ 2D angular distributions in the helicity reference frame projected along cosθ(left) and φ(right) for 2.5 <y <4and 4 <pT<6GeV/c. The displayed uncertainties are statistical. rected for acceptance and efficiency, according to Eq. (1). For each combination of signal and background shape used in the fit to the dimuon invariant-mass spectra, a separate evaluation of the polarization parameters is carried out and their average is taken as the best estimate. The statistical uncertainty is given by the average of the statistical uncertainties of the 2D fits, while the root mean square of the results provides the systematic uncertainty on the signal extraction, with the absolute values ranging between 0.002 and 0.039. The overall procedure described above was checked beforehand with a MC closure test. The 2D fits on the (cosθ, φ) distributions only allow a determination of the absolute value of λθφ, due to the presence of sin2θin the corresponding term that induces an ambiguity in its sign. It is checked that the values of λθ and λφare stable against the choice of the sign of the λθφ term. In the following the λθφ values corresponding to the choice of a positive sign are quoted. Fig. 2illustrates an example of the fit to the angular distributions. For better visibility, both the distribution and the fitted function are projected along one dimension. In addition to the systematic uncertainty related to the choice of the mass shapes for signal and background, several other sources are taken into account. First, an alternative procedure for extracting the J/ψ signal is carried out, by keeping its width as a free parameter in the invariant-mass fits. The corresponding results for the polarization parameters are then obtained and the averages of the values corresponding to fixing the width or not are taken as the central values for λθ, λφand λθφ. Half the difference between the results obtained with free or MC-anchored widths is then considered as a further systematic uncertainty related to the signal extraction. This uncertainty is found to be the leading contribution to the total absolute systematic uncertainty on the polarization parameters, and ranges between 0.001 and 0.063, the latter value corresponding to the uncertainty on λHE θfor 2 <pT<4GeV/c. Another source of systematic uncertainty is related to the evaluation of the trigger efficiency. The muon trigger response function as a function of the single muon transverse momentum pμ Tcan be obtained via MC or with a procedure based on data [59]. Small deviations are found for pμ t<2GeV/cwhich induce an effect on A ×εfor the J/ψ. Therefore, the polarization parameters are recalculated with A ×εvalues weighted in such a way to account for the deviations. The variation of the polarization parameters between the different trigger efficiency estimates is taken as the related systematic uncertainty, with values ranging from 0.001 to 0.043, the highest values being found for λHE θin 2 <pT<4GeV/c. The systematic uncertainty related to the evaluation of the muon tracking efficiency is found to be negligible for this analysis, allowing a significant reduction of the total systematic uncertainty with respect to previous pp analyses [23]. Indeed, although the difference between efficiencies calculated via MC or from data [59]is of the order of 2%, a detailed investigation has shown no dependence on the angular variables and therefore no effect on the polarization parameters. Finally, the systematic uncertainty induced by the choice of the pTand ydistributions used as an input for the calculation of A ×εis evaluated testing alternative pTand yparameterizations, which are obtained by varying within their uncertainties the default distributions directly tuned on Pb–Pb data. The polarization parameters extracted with the modified values of A ×εare compared with those obtained with the default input shapes and the corresponding systematic uncertainty extracted in this way is found to range between 0.001 and 0.030, with the largest value assigned to λHE θfor 2 <pT<4GeV/c. The influence of the choice of the angular distributions of the J/ψdecay products for the A ×ε calculation is also investigated by means of an iterative procedure on these input distributions. The effect is found to be negligible, also due to the fact that the 2D correction procedure on the angular variables is by definition relatively insensitive to the specific choice of the corresponding distributions. A summary of the values of all the absolute systematic uncertainties, which are considered as uncorrelated as a function of pT, is reported in Table 1. The total systematic uncertainties are obtained, for each parameter and pTinterval, as the quadratic sum of the values. A similar procedure is followed for the extraction of the ϒ(1S) polarization parameters. Due to the smaller candidate sample, integrated values over the kinematic interval 2.5 <y <4, pT< 15 GeV/care obtained. The main difference with respect to the 2D approach followed for the J/ψis the use of a simultaneous fit to the 1D angular distributions [23], after integration over the other variables. The requirement pμ T>2GeV/c, which helps reducing the combinatorial background, is included [60]. The ϒ(1S)signal extraction in the various cosθand φintervals is performed by means of invariant-mass fits (see the right panel of Fig. 1for an example). The functions chosen for the resonances are the same as in the J/ψ analysis (pseudo-Gaussian or Crystal Ball), the mass value is fixed to that obtained from a fit to the integrated invariant-mass distribution, while the width for each angular interval is fixed to the MC value scaled by the ratio of the widths between data and MC for the angle-integrated distributions. The tail parameters are fixed to MC values. The small contribution from ϒ(2S)is also included in the fits [60]. The background continuum is parameterized with a Gaussian distribution whose width varies linearly with the mass 4 ALICE Collaboration Physics Letters B 815 (2021) 136146 Table 1 Summary of the absolute systematic uncertainties on the evaluation of the J/ψ polarization parameters. All the uncertainties are considered as uncorrelated as a function of pT. Helicity Collins-Soper Signal J/ψTrigger Input Signal J/ψTrigger Input pT(GeV/c) extr. width eff. MC extr. width eff. MC λθ2<pT<4 0.030 0.063 0.043 0.030 0.026 0.049 0.015 0.009 4<pT<6 0.017 0.046 0.040 0.024 0.002 0.052 0.018 0.007 6<pT<10 0.039 0.005 0.018 0.017 0.022 0.001 0.011 0.006 λφ2<pT<4 0.007 0.030 0.004 0.002 0.024 0.010 0.020 0.003 4<pT<6 0.003 0.035 0.003 0.003 0.002 0.010 0.020 0.003 6<pT<10 0.002 0.009 0.001 0.002 0.005 0.013 0.011 0.002 λθφ 2<pT<4 0.021 0.029 0.024 0.001 0.013 0.010 0.017 0.015 4<pT<6 0.007 0.011 0.017 0.006 0.002 0.042 0.010 0.015 6<pT<10 0.020 0.019 0.007 0.008 0.007 0.042 0.003 0.013 Table 2 J/ψ polarization parameters, measured for Pb–Pb collisions at √sNN =5.02 TeV, in the helicity and Collins-Soper reference frames in the rapidity interval 2.5 <y <4. The first uncertainty is statistical and the second systematic. pT(GeV/c) Helicity Collins-Soper λθ2<pT<40.218 ±0.060 ±0.087 −0.157 ±0.049 ±0.058 4<pT<60.151 ±0.071 ±0.068 −0.057 ±0.059 ±0.055 6<pT<10 −0.070 ±0.068 ±0.047 −0.008 ±0.063 ±0.026 λφ2<pT<4−0.029 ±0.017 ±0.031 0.061 ±0.015 ±0.033 4<pT<6−0.013 ±0.019 ±0.036 0.047 ±0.024 ±0.023 6<pT<10 0.047 ±0.021 ±0.010 0.024 ±0.032 ±0.018 λθφ 2<pT<4−0.124 ±0.028 ±0.043 −0.090 ±0.027 ±0.029 4<pT<6−0.059 ±0.030 ±0.021 −0.040 ±0.034 ±0.046 6<pT<10 −0.025 ±0.031 ±0.030 0.018 ±0.035 ±0.044 or, alternatively, with a second degree polynomial function times an exponential. The systematic uncertainty on the signal extraction is calculated with the same procedure adopted for the J/ψ. An uncertainty related to the choice of the signal width has also been considered, taken as the half-difference between the results obtained with the prescription described above and using as an alternative prescription the pure MC values. The uncertainty on the trigger efficiency is negligible, due to the additional requirement on the single-muon transverse momentum which selects a pT-region where the trigger efficiency is very high and its evaluation via data and MC is consistent. Finally, the procedure for the determination of the uncertainty related to the ϒ(1S)kinematic distributions used in the MC is the same as for the J/ψ. The total systematic uncertainties for the ϒ(1S)analysis are reported in Table 3, together with the results. 4. Results The polarization parameters for J/ψ inclusive production in Pb– Pb collisions at √sNN =5.02 TeV in the helicity and Collins-Soper reference frames are shown in Fig. 3and the corresponding numerical values are reported in Table 2. In Fig. 3, λθ, λφand λθφ are also compared with the LHCb [24] and ALICE [23]measurements in pp collisions at √s=7 and 8 TeV, respectively. For all the pTintervals and in both reference frames the values of the polarization parameters exhibit at most slight deviations from zero. In particular, λHE θindicates a slight transverse polarization at low pT(∼2.1σeffect, calculated using the Gaussian approximation), while λCS θshows a weak longitudinal polarization (∼2.1σ). When increasing pT, the central values of λθbecome close to zero. All values of λφand λθφ are, in absolute value, smaller than 0.1, except for λHE θφ, which is −0.124 at low pTand deviates from zero by ∼2.4σ. When comparing with the pp results, no significant difference is found with respect to ALICE results at √s=8TeV, which are compatible with zero. A significant difference is found with respect to the higher-precision LHCb results at √s=7TeV, reaching Table 3 ϒ(1S)polarization parameters in the helicity and Collins-Soper reference frames measured in Pb–Pb collisions at √sNN =5.02 TeV in the rapidity interval 2.5 <y < 4and for transverse momentum pT<15 GeV/c. The first uncertainty is statistical and the second systematic. Helicity Collins-Soper λθ−0.090 ±0.395 ±0.101 0.418 ±0.526 ±0.178 λφ−0.094 ±0.072 ±0.020 −0.141 ±0.087 ±0.033 λθφ −0.074 ±0.099 ±0.020 0.017 ±0.113 ±0.024 3.3σin the interval 2 <pT<4GeV/cin the helicity reference frame, where pp data [24] indicate a small but significant degree of longitudinal polarization, while the Pb–Pb results favor a slightly transverse polarization. In Pb–Pb collisions at LHC energies, a significant fraction of the detected J/ψoriginates from the recombination of ccpairs in the QGP phase or when the system hadronizes. Moreover, the contribution from higher-mass charmonium states decaying to J/ψcould vary between pp and Pb-Pb due to different suppression effects for each state in nuclear collisions. Therefore, the observed hint for a different polarization in pp and Pb–Pb might be a reflection of the different production and suppression mechanisms in the two systems, but more precise data, along with quantitative theory estimates, are needed for a definite conclusion. It should also be noted that the ALICE results refer to inclusive production, while LHCb has measured prompt J/ψ. However, as discussed in Ref. [22], the size of the non-prompt component is small in the covered pTregion (of the order of 15% at high pT) and its polarization was also measured to be small by CDF (λHE θ∼−0.1[21]), implying that the net effect of this source on inclusive J/ψpolarization should be negligible. In Table 3the values of the ϒ(1S)polarization parameters are shown. The λθvalues are consistent with zero, with large uncertainties that prevent a firm conclusion on the absence of polarization in nuclear collisions. The λφand λθφ values are also consistent with zero. The relatively smaller uncertainties for these parameters are related to a more uniform acceptance distribution as a function of the azimuthal angular variable. 5 ALICE Collaboration Physics Letters B 815 (2021) 136146 Fig. 3. Inclusive J/ψ polarization parameters as a function of transverse momentum for Pb–Pb collisions at √sNN =5.02 TeV, compared with results obtained in pp collisions by ALICE at √s=8TeV [23]and by LHCb for prompt J/ψ at √s=7TeV [24](the LHCb markers were shifted horizontally by +0.3 GeV/cfor better visibility) in the rapidity interval 3 <y <3.5. The error bars represent the total uncertainties for the pp results, while for Pb–Pb statistical and systematic uncertainties are plotted separately as a vertical bar and a shaded box, respectively. In the left part of the plot the polarization parameters in the helicity reference frame are reported, in the right those for the Collins-Soper frame. 5. Conclusions The first measurement of the polarization parameters for J/ψ production in nuclear collisions at LHC energies was carried out by the ALICE Collaboration in Pb–Pb interactions at √sNN =5.02 TeV. The λθ, λφand λθφ parameters were evaluated in the helicity and Collins-Soper reference frames in the rapidity interval 2.5 <y <4 and in the transverse momentum interval 2 <pT<10 GeV/c. All the parameter values are close to zero, with a ∼2.1σindication for a small transverse polarization in the helicity frame at low pT, and a corresponding indication for a small longitudinal polarization in the Collins-Soper frame (∼2.1σeffect). When comparing these results with pp data taken at higher energy at the LHC, an interesting feature is a significant difference in λHE θwith respect to the LHCb results which showed instead a small longitudinal polarization in a similar kinematic domain. This first result obtained for J/ψin nuclear collisions and described in this Letter represents therefore a starting point for future studies connecting such features with the known differences in the production mechanisms between pp and nucleus–nucleus collisions. Results were also obtained for the first time for the ϒ(1S)polarization, integrated over pTand y, showing, within the large uncertainties of the measurement, values compatible with the absence of polarization. Declaration of competing interest The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper. Acknowledgements The ALICE Collaboration would like to thank all its engineers and technicians for their invaluable contributions to the construction of the experiment and the CERN accelerator teams for the outstanding performance of the LHC complex. The ALICE Collaboration gratefully acknowledges the resources and support provided by all Grid centres and the Worldwide LHC Computing Grid (WLCG) collaboration. The ALICE Collaboration acknowledges the following funding agencies for their support in building and running the ALICE detector: A. I. Alikhanyan National Science Laboratory (Yerevan Physics Institute) Foundation (ANSL), State Committee of Science and World Federation of Scientists (WFS), Armenia; Austrian Academy of Sciences, Austrian Science Fund (FWF): [M 2467-N36] and Österreichische Nationalstiftung für Forschung, Technologie und Entwicklung, Austria; Ministry of Communications and High Technologies, National Nuclear Research Center, Azerbaijan; Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq), Financiadora de Estudos e Projetos (Finep), Fun6 ALICE Collaboration Physics Letters B 815 (2021) 136146 dação de Amparo à Pesquisa do Estado de São Paulo (FAPESP) and Universidade Federal do Rio Grande do Sul (UFRGS), Brazil; Ministry of Education of China (MOEC), Ministry of Science & Technology of China (MSTC) and National Natural Science Foundation of China (NSFC), China; Ministry of Science and Education and Croatian Science Foundation, Croatia; Centro de Aplicaciones Tecnológicas y Desarrollo Nuclear (CEADEN), Cubaenergía, Cuba; Ministry of Education, Youth and Sports of the Czech Republic, Czech Republic; The Danish Council for Independent Research | Natural Sciences, the Villum Fonden and Danish National Research Foundation (DNRF), Denmark; Helsinki Institute of Physics (HIP), Finland; Commissariat à l’Énergie Atomique (CEA) and Institut National de Physique Nucléaire et de Physique des Particules (IN2P3) and Centre National de la Recherche Scientifique (CNRS), France; Bundesministerium für Bildung und Forschung (BMBF) and GSI Helmholtzzentrum für Schwerionenforschung GmbH, Germany; General Secretariat for Research and Technology, Ministry of Education, Research and Religions, Greece; National Research, Development and Innovation Office, Hungary; Department of Atomic Energy, Government of India (DAE), Department of Science and Technology, Government of India (DST), University Grants Commission, Government of India (UGC) and Council of Scientific and Industrial Research (CSIR), India; Indonesian Institute of Science, Indonesia; Centro Fermi -Museo Storico della Fisica e Centro Studi e Ricerche Enrico Fermi and Istituto Nazionale di Fisica Nucleare (INFN), Italy; Institute for Innovative Science and Technology, Nagasaki Institute of Applied Science (IIST), Japanese Ministry of Education, Culture, Sports, Science and Technology (MEXT) and Japan Society for the Promotion of Science (JSPS) KAKENHI, Japan; Consejo Nacional de Ciencia (CONACYT) y Tecnología, through Fondo de Cooperación Internacional en Ciencia y Tecnología (FONCICYT) and Dirección General de Asuntos del Personal Academico (DGAPA), Mexico; Nederlandse Organisatie voor Wetenschappelijk Onderzoek (NWO), Netherlands; The Research Council of Norway, Norway; Commission on Science and Technology for Sustainable Development in the South (COMSATS), Pakistan; Pontificia Universidad Católica del Perú, Peru; Ministry of Science and Higher Education, National Science Centre and WUT ID-UB, Poland; Korea Institute of Science and Technology Information and National Research Foundation of Korea (NRF), Republic of Korea; Ministry of Education and Scientific Research, Institute of Atomic Physics and Ministry of Research and Innovation and Institute of Atomic Physics, Romania; Joint Institute for Nuclear Research (JINR), Ministry of Education and Science of the Russian Federation, National Research Centre Kurchatov Institute, Russian Science Foundation and Russian Foundation for Basic Research, Russia; Ministry of Education, Science, Research and Sport of the Slovak Republic, Slovakia; National Research Foundation of South Africa, South Africa; Swedish Research Council (VR) and Knut & Alice Wallenberg Foundation (KAW), Sweden; European Organization for Nuclear Research, Switzerland; Suranaree University of Technology (SUT), National Science and Technology Development Agency (NSDTA) and Office of the Higher Education Commission under NRU project of Thailand, Thailand; Turkish Atomic Energy Agency (TAEK), Turkey; National Academy of Sciences of Ukraine, Ukraine; Science and Technology Facilities Council (STFC), United Kingdom; National Science Foundation of the United States of America (NSF) and United States Department of Energy, Office of Nuclear Physics (DOE NP), United States of America. References [1] Particle Data Group Collaboration, M. Tanabashi, et al., Review of particle physics, Phys. Rev. D 98 (3) (2018) 030001. [2] G.T. Bodwin, E. Braaten, G.P. Lepage, Rigorous QCD analysis of inclusive annihilation and production of heavy quarkonium, Phys. Rev. D 51 (1995) 1125–1171, arXiv:hep -ph /9407339; G.T. Bodwin, E. Braaten, G.P. Lepage, Phys. Rev. D 55 (1997) 5853, Erratum. [3] ALICE Collaboration, S. Acharya, et al., Energy dependence of forward-rapidity J/ψ and ψ(2S)production in pp collisions at the LHC, Eur. Phys. J. C 77 (6) (2017) 392, arXiv:1702 .00557 [hep -ex]. [4] ALICE Collaboration, S. Acharya, et al., Inclusive J/ψproduction at mid-rapidity in pp collisions at √s= 5.02 TeV, J. High Energy Phys. 10 (2019) 084, arXiv: 1905 .07211 [nucl -ex]. [5] LHCb Collaboration, R. Aaij, et al., Measurement of forward J/ψ production cross-sections in pp collisions at √s=13 TeV, J. High Energy Phys. 10 (2015) 172, arXiv:1509 .00771 [hep -ex ]. [Erratum : JHEP05 ,063(2017 )]. [6] CMS Collaboration, A.M. Sirunyan, et al., Measurement of quarkonium production cross sections in pp collisions at √s=13 TeV, Phys. Lett. B 780 (2018) 251–272, arXiv:1710 .11002 [hep -ex]. [7] ATLAS Collaboration, G. Aad, et al., Measurement of the differential crosssections of prompt and non-prompt production of J/ψ and ψ(2S)in pp collisions at √s=7and 8 TeV with the ATLAS detector, Eur. Phys. J. C 76 (5) (2016) 283, arXiv:1512 .03657 [hep -ex]. [8] PHENIX Collaboration, U. Acharya, et al., J/ψ and ψ(2S)production at forward rapidity in p+pcollisions at √s=510 GeV, Phys. Rev. D 101 (5) (2020) 052006, arXiv:1912 .13424 [hep -ex]. [9] STAR Collaboration, J. Adam, et al., Measurements of the transversemomentum-dependent cross sections of J/ψ production at mid-rapidity in proton+proton collisions at √s=510 and 500 GeV with the STAR detector, Phys. Rev. D 100 (5) (2019) 052009, arXiv:1905 .06075 [hep -ex]. [10] CDF Collaboration, F. Abe, et al., J/ψ and ψ(2S)production in p¯ pcollisions at √s=1.8TeV, Phys. Rev. Lett. 79 (1997) 572–577. [11] D0 Collaboration, S. Abachi, et al., J/ψ production in p¯ pcollisions at √s=1.8TeV, Phys. Lett. B 370 (1996) 239–248. [12] LHCb Collaboration, R. Aaij, et al., Measurement of ϒproduction in pp collisions at √s=13 TeV, J. High Energy Phys. 07 (2018) 134, arXiv:1804 .09214 [hep -ex]; R. Aaij, et al., J. High Energy Phys. 05 (2019) 076, Erratum. [13] ALICE Collaboration, J. Adam, et al., Inclusive quarkonium production at forward rapidity in pp collisions at √s=8TeV, Eur. Phys. J. C 76 (4) (2016) 184, arXiv:1509 .08258 [hep -ex]. [14] ATLAS Collaboration, G. Aad, et al., Measurement of Upsilon production in 7 TeV pp collisions at ATLAS, Phys. Rev. D 87 (5) (2013) 052004, arXiv:1211.7255 [hep -ex]. [15] M. Einhorn, S. Ellis, Hadronic production of the new resonances: probing gluon distributions, Phys. Rev. D 12 (1975) 2007. [16] M. Gluck, J.F. Owens, E. Reya, Gluon contribution to hadronic J/ψproduction, Phys. Rev. D 17 (1978) 2324. [17] S. Baranov, A. Lipatov, Prompt charmonia production and polarization at LHC in the NRQCD with kT-factorization. Part III: J/ψ meson, Phys. Rev. D 96 (3) (2017) 034019, arXiv:1611.10141 [hep -ph]. [18] E. Braaten, B.A. Kniehl, J. Lee, Polarization of prompt J/ψ at the Tevatron, Phys. Rev. D 62 (2000) 094005, arXiv:hep -ph /9911436 [hep -ph]. [19] P. Faccioli, C. Lourenco, J. Seixas, H.K. Wohri, Towards the experimental clarification of quarkonium polarization, Eur. Phys. J. C 69 (2010) 657–673, arXiv: 1006 .2738 [hep -ph]. [20] CDF Collaboration, T. Affolder, et al., Measurement of J/ψ and ψ(2S)polarization in p¯ pcollisions at √s=1.8TeV, Phys. Rev. Lett. 85 (2000) 2886–2891, arXiv:hep -ex /0004027 [hep -ex]. [21] CDF Collaboration, A. Abulencia, et al., Polarization of J/ψ and ψ(2S)mesons produced in p¯ pcollisions at √s=1.96-TeV, Phys. Rev. Lett. 99 (2007) 132001, arXiv:0704 .0638 [hep -ex]. [22] ALICE Collaboration, B. Abelev, et al., J/ψ polarization in pp collisions at √s= 7TeV, Phys. Rev. Lett. 108 (2012) 082001, arXiv:1111.1630 [hep -ex]. [23] ALICE Collaboration, S. Acharya, et al., Measurement of the inclusive J/ ψpolarization at forward rapidity in pp collisions at √s=8TeV, Eur. Phys. J. C 78 (7) (2018) 562, arXiv:1805 .04374 [hep -ex]. [24] LHCb Collaboration, R. Aaij, et al., Measurement of J/ψ polarization in pp collisions at √s=7TeV, Eur. Phys. J. C 73 (11) (2013) 2631, arXiv:1307.6379 [hep -ex]. [25] CMS Collaboration, S. Chatrchyan, et al., Measurement of the prompt J/ψ and ψ(2S) polarizations in pp collisions at √s=7TeV, Phys. Lett. B 727 (2013) 381–402, arXiv:1307.6070 [hep -ex]. [26] M. Butenschoen, B.A. Kniehl, Reconciling J/ψ production at HERA, RHIC, Tevatron, and LHC with NRQCD factorization at next-to-leading order, Phys. Rev. Lett. 106 (2011) 022003, arXiv:1009 .5662 [hep -ph]. [27] Y.-Q. Ma, K. Wang, K.-T. Chao, A complete NLO calculation of the J/ψ and ψ production at hadron colliders, Phys. Rev. D 84 (2011) 114001, arXiv:1012 .1030 [hep -ph]. [28] M. Butenschoen, B.A. Kniehl, J/ψpolarization at tevatron and LHC: nonrelativistic-QCD factorization at the crossroads, Phys. Rev. Lett. 108 (2012) 172002, arXiv:1201.1872 [hep -ph]. [29] K.-T. Chao, Y.-Q. Ma, H.-S. Shao, K. Wang, Y.-J. Zhang, J/ψ polarization at hadron colliders in nonrelativistic QCD, Phys. Rev. Lett. 108 (2012) 242004, arXiv:1201.2675 [hep -ph]. 7 ALICE Collaboration Physics Letters B 815 (2021) 136146 [30] B. Gong, L.-P. Wan, J.-X. Wang, H.-F. Zhang, Polarization for prompt j/ψand ψ(2s)production at the tevatron and LHC, Phys. Rev. Lett. 110 (4) (2013) 042002, arXiv:1205 .6682 [hep -ph]. [31] Y. Feng, B. Gong, C.-H. Chang, J.-X. Wang, Remaining parts of the long-standing J/ψ polarization puzzle, Phys. Rev. D 99 (1) (2019) 014044, arXiv:1810 .08989 [hep -ph]. [32] Y.-Q. Ma, T. Stebel, R. Venugopalan, J/ψ polarization in the CGC+NRQCD approach, J. High Energy Phys. 12 (2018) 057, arXiv:1809 .03573 [hep -ph]. [33] CMS Collaboration, S. Chatrchyan, et al., Measurement of the ϒ(1S), ϒ(2S) and ϒ(3S) polarizations in pp collisions at √s=7TeV, Phys. Rev. Lett. 110 (8) (2013) 081802, arXiv:1209 .2922 [hep -ex]. [34] CMS Collaboration, V. Khachatryan, et al., ϒ(nS)polarizations versus particle multiplicity in pp collisions at √s=7TeV, Phys. Lett. B 761 (2016) 31–52, arXiv:1603 .02913 [hep -ex]. [35] LHCb Collaboration, R. Aaij, et al., Measurement of the ϒpolarizations in pp collisions at √s=7and 8 TeV, J. High Energy Phys. 12 (2017) 110, arXiv: 1709 .01301 [hep -ex]. [36] H. Han, Y.-Q. Ma, C. Meng, H.-S. Shao, Y.-J. Zhang, K.-T. Chao, ϒ(nS)and χb(nP) production at hadron colliders in nonrelativistic QCD, Phys. Rev. D 94 (1) (2016) 014028, arXiv:1410 .8537 [hep -ph]. [37] P. Braun-Munzinger, J. Wambach, The phase diagram of strongly-interacting matter, Rev. Mod. Phys. 81 (2009) 1031–1050, arXiv:0801.4256 [hep -ph]. [38] P. Braun-Munzinger, V. Koch, T. Schäfer, J. Stachel, Properties of hot and dense matter from relativistic heavy ion collisions, Phys. Rep. 621 (2016) 76–126, arXiv:1510 .00442 [nucl -th]. [39] T. Matsui, H. Satz, J/ψsuppression by quark-gluon plasma formation, Phys. Lett. B 178 (1986) 416–422. [40] M. Laine, A resummed perturbative estimate for the quarkonium spectral function in hot QCD, J. High Energy Phys. 05 (2007) 028, arXiv:0704 .1720 [hep -ph]. [41] P. Braun-Munzinger, J. Stachel, (Non)thermal aspects of charmonium production and a new look at J/ψsuppression, Phys. Lett. B 490 (2000) 196–202, arXiv:nucl -th /0007059 [nucl -th]. [42] R.L. Thews, M. Schroedter, J. Rafelski, Enhanced J/ψ production in deconfined quark matter, Phys. Rev. C 63 (2001) 054905, arXiv:hep -ph /0007323 [hep -ph]. [43] B.L. Ioffe, D.E. Kharzeev, Quarkonium polarization in heavy ion collisions as a possible signature of the quark gluon plasma, Phys. Rev. C 68 (2003) 061902, arXiv:hep -ph /0306176 [hep -ph]. [44] H.-S. Shao, Y.-Q. Ma, K. Wang, K.-T. Chao, Polarizations of χc1and χc2in prompt production at the LHC, Phys. Rev. Lett. 112 (18) (2014) 182003, arXiv: 1402 .2913 [hep -ph]. [45] CMS Collaboration, A.M. Sirunyan, et al., Measurement of the χc1 and χc2 polarizations in proton-proton collisions at √s=8TeV, arXiv:1912 .07706 [hep -ex]. [46] ALICE Collaboration, K. Aamodt, et al., The ALICE experiment at the CERN LHC, J. Instrum. 3 (2008) S08002. [47] ALICE Collaboration, B. Abelev, et al., Performance of the ALICE Experiment at the CERN LHC, Int. J. Mod. Phys. A 29 (2014) 1430044, arXiv:1402 .4476 [nucl - ex]. [48] ALICE Collaboration, K. Aamodt, et al., Rapidity and transverse momentum dependence of inclusive J/ψ production in pp collisions at √s=7TeV, Phys. Lett. B 704 (2011) 442, arXiv:1105 .0380 [hep -ex]. [49] ALICE Collaboration, E. Abbas, et al., Performance of the ALICE VZERO system, J. Instrum. 8 (2013) P10016, arXiv:1306 .3130 [nucl -ex]. [50] ALICE Collaboration, K. Aamodt, et al., Alignment of the ALICE Inner Tracking System with cosmic-ray tracks, J. Instrum. 5 (2010) P03003, arXiv:1001.0502 [physics .ins -det]. [51] ALICE Collaboration, B. Abelev, et al., Measurement of the cross section for electromagnetic dissociation with neutron emission in Pb-Pb collisions at √sNN = 2.76 TeV, Phys. Rev. Lett. 109 (2012) 252302, arXiv:1203 .2436 [nucl -ex]. [52] ALICE Collaboration, F. Bossu, M. Gagliardi, M. Marchisone, Performance of the RPC-based ALICE muon trigger system at the LHC, J. Instrum. 7 (2012) T12002, arXiv:1211.1948 [physics .ins -det ]. [PoSRPC2012 , 059(2012 )]. [53] ALICE Collaboration, B. Abelev, et al., Centrality determination of Pb-Pb collisions at √sNN =2.76 TeV with ALICE, Phys. Rev. C 88 (4) (2013) 044909, arXiv:1301.4361 [nucl -ex]. [54] ALICE Collaboration, Centrality determination in heavy ion collisions, ALICEPUBLIC-2018-011, http://cds .cern .ch /record /2636623. [55] ALICE Collaboration, Quarkonium signal extraction in ALICE, ALICE-PUBLIC2015-006, https://cds .cern .ch /record /2060096. [56] ALICE Collaboration, S. Acharya, et al., Studies of J/ψproduction at forward rapidity in Pb-Pb collisions at √sNN =5.02 TeV, J. High Energy Phys. 02 (2020) 041, arXiv:1909 .03158 [nucl -ex]. [57] ALICE Collaboration, S. Acharya, et al., Inclusive J/ψproduction at forward and backward rapidity in p-Pb collisions at √sNN =8.16 TeV, J. High Energy Phys. 07 (2018) 160, arXiv:1805 .04381 [nucl -ex]. [58] R. Brun, F. Bruyant, F. Carminati, S. Giani, M. Maire, A. McPherson, G. Patrick, L. Urban, GEANT: Detector Description and Simulation Tool; Oct 1994. CERN Program Library. CERN, Geneva, 1993, http://cds .cern .ch /record /1082634, Long Writeup W5013. [59] ALICE Collaboration, B. Abelev, et al., Measurement of quarkonium production at forward rapidity in pp collisions at √s=7TeV, Eur. Phys. J. C 74 (8) (2014) 2974, arXiv:1403 .3648 [nucl -ex]. [60] ALICE Collaboration, S. Acharya, et al., ϒsuppression at forward rapidity in Pb-Pb collisions at √sNN =5.02 TeV, Phys. Lett. B 790 (2019) 89–101, arXiv: 1805 .04387 [nucl -ex]. ALICE Collaboration S. Acharya141, D. Adamová95, A. Adler74, J. Adolfsson81, M.M. Aggarwal 100, G. Aglieri Rinella 34, M. Agnello30, N. Agrawal 10,54, Z. Ahammed141, S. Ahmad16, S.U. Ahn76, Z. Akbar51, A. Akindinov92, M. Al-Turany107, S.N. Alam40,141, D.S.D. Albuquerque122, D. Aleksandrov 88, B. Alessandro59, H.M. Alfanda6, R. Alfaro Molina71, B. Ali 16, Y. Ali 14, A. Alici10,26,54, N. Alizadehvandchali125, A. Alkin2,34, J. Alme21, T. Alt 68, L. Altenkamper 21, I. Altsybeev113, M.N. Anaam 6, C. Andrei48, D. Andreou34, A. Andronic144, M. Angeletti34, V. Anguelov104, C. Anson 15, T. Antiˇ ci´ c108, F. Antinori57, P. Antonioli 54, N. Apadula80, L. Aphecetche 115, H. Appelshäuser68, S. Arcelli 26, R. Arnaldi 59, M. Arratia80, I.C. Arsene 20, M. Arslandok 104, A. Augustinus34, R. Averbeck 107, S. Aziz78, M.D. Azmi16, A. Badalà56, Y.W. Baek 41, S. Bagnasco 59, X. Bai107, R. Bailhache68, R. Bala 101, A. Balbino30, A. Baldisseri137, M. Ball 43, S. Balouza 105, D. Banerjee3, R. Barbera27, L. Barioglio25, G.G. Barnaföldi145, L.S. Barnby94, V. Barret134, P. Bartalini 6, C. Bartels127, K. Barth 34, E. Bartsch68, F. Baruffaldi 28, N. Bastid134, S. Basu 143, G. Batigne 115, B. Batyunya 75, D. Bauri49, J.L. Bazo Alba 112, I.G. Bearden89, C. Beattie146, C. Bedda63, N.K. Behera61, I. Belikov136, A.D.C. Bell Hechavarria144, F. Bellini34, R. Bellwied 125, V. Belyaev93, G. Bencedi145, S. Beole 25, A. Bercuci48, Y. Berdnikov 98, D. Berenyi145, R.A. Bertens130, D. Berzano59, M.G. Besoiu67, L. Betev34, A. Bhasin101, I.R. Bhat101, M.A. Bhat 3, H. Bhatt49, B. Bhattacharjee42, A. Bianchi25, L. Bianchi25, N. Bianchi52, J. Bielˇ cík37, J. Bielˇ cíková95, A. Bilandzic105, G. Biro145, R. Biswas3, S. Biswas3, J.T. Blair 119, D. Blau88, C. Blume68, G. Boca139, F. Bock 96, A. Bogdanov93, S. Boi23, J. Bok61, L. Boldizsár 145, A. Bolozdynya93, M. Bombara38, G. Bonomi140, H. Borel137, A. Borissov93, H. Bossi146, E. Botta25, L. Bratrud68, P. Braun-Munzinger107, M. Bregant 121, M. Broz37, E. Bruna 59, G.E. Bruno33,106, M.D. Buckland127, D. Budnikov 109, H. Buesching68, S. Bufalino30, O. Bugnon115, P. Buhler 114, P. Buncic 34, Z. Buthelezi72,131, J.B. Butt14, 8