Abstract
The polarisation of prompt \(\psi {(2S)} \) mesons is measured by performing an angular analysis of \(\psi {(2S)} \!\rightarrow \mu ^+\mu ^- \) decays using proton-proton collision data, corresponding to an integrated luminosity of 1.0\(\,\text{ fb }^{-1} \), collected by the LHCb detector at a centre-of-mass energy of 7 TeV. The polarisation is measured in bins of transverse momentum \(p_\mathrm{T} \) and rapidity \(y\) in the kinematic region \(3.5< p_\mathrm{T} <15{\mathrm {\,GeV\!/}c} \) and \(2.0<y<4.5\), and is compared to theoretical models. No significant polarisation is observed.
Similar content being viewed by others
Explore related subjects
Discover the latest articles, news and stories from top researchers in related subjects.Avoid common mistakes on your manuscript.
1 Introduction
Measurements of the heavy quarkonium production in hadron collisions can be used to test predictions of quantum chromodynamics (QCD) in the perturbative and non-perturbative regimes. Several theoretical models have been developed within the framework of QCD to describe the quarkonium production cross-section and polarisation as functions of the quarkonium transverse momentum, \(p_\mathrm{T}\), but none can simultaneously describe both of them [1]. Heavy quarkonia can be produced in three ways in \(pp\) collisions: directly in the hard scattering, through feed-down from higher-mass quarkonia states, or via the decay of \(b\) hadrons, with the first two of these being referred to as prompt production. In the case of \(\psi {(2S)} \) mesons, the contribution from feed-down is negligible, allowing a straightforward comparison between measurements of prompt production and predictions for direct contributions.
The \(\psi {(2S)} \) meson has spin, parity and charge-parity quantum numbers, \(J^{\mathrm {PC}}=1^{--}\) and its polarisation can be determined by studying the angular distribution of muons in the \(\psi {(2S)} \!\rightarrow \mu ^+ \mu ^- \) decays [2, 3]. The distribution is described by
where \(\theta \) and \(\phi \) are the polar and azimuthal angles of the \(\mu ^+ \) direction in the rest frame of the \(\psi {(2S)} \) meson, respectively, and \(\lambda _\theta ,\lambda _{\theta \phi }\) and \(\lambda _\phi \) are the polarisation parameters to be determined from the data. The case of \((\lambda _{\theta }, \lambda _{\theta \phi }, \lambda _{\phi } ) = (1,\,0,\,0)\) or \((-1\,,0,\,0)\) corresponds to full transverse or longitudinal polarisation, respectively, while \((\lambda _\theta ,\lambda _{\theta \phi },\lambda _\phi )=(0,\,0,\,0)\) corresponds to the unpolarised state.Footnote 1 In this study of the \(\psi {(2S)} \) polarisation, two choices of polarisation frame are used. These have a common definition of the \(Y\)-axis, taken to be the normal to the production plane, which is formed by the momentum of the \(\psi {(2S)} \) meson and the beam axis in the rest frame of the colliding LHC protons. The helicity frame [4] uses the \(\psi {(2S)} \) momentum as the \(Z\)-axis. In the Collins-Soper frame [5] the \(Z\)-axis is chosen to be the bisector of the angle between the two incoming proton beams in the rest frame of the \(\psi {(2S)} \) meson. In both frames, the \(X\)-axis is defined to complete a right-handed Cartesian coordinate system. The commonly used frame-invariant variable \(\lambda _\mathrm {inv}\) (see [6, 7]) is defined as
Two classes of theoretical models are compared with the measurements in this paper: the colour-singlet model (CSM) [8] and non-relativistic QCD (NRQCD) [9–14], at next-to-leading order (NLO). In the high-\(p_\mathrm{T}\) region, where the quarkonium transverse momentum is much larger than its mass (in natural units), the CSM underestimates significantly the measured prompt \({J /\psi }\) and \(\psi {(2S)}\) production cross-sections [15–17], while the NRQCD model provides a good description of the \(p_\mathrm{T}\)-dependent \({J /\psi }\) and \(\psi {(2S)}\) cross-sections measured by LHCb [16, 17] and CMS [18]. The CSM predicts large longitudinal polarisation for \({J /\psi }\) and \(\psi {(2S)}\) mesons. On the other hand, in the NRQCD model, where quarkonium production is dominated by the gluon fragmentation process in the high-\(p_\mathrm{T}\) region, the gluon is almost on-shell, leading to predictions of large transverse polarisations [11]. Precise measurements of the \({J /\psi }\) polarisation at both the Tevatron [19] and the LHC [20–22], which show no significant longitudinal or transverse polarisations, do not support either the CSM or NRQCD predictions.
The prompt \(\psi {(2S)}\) polarisation has been measured by the CDF experiment [19] in \(p\overline{p}\) collisions at \(\sqrt{s}=1.96\,\mathrm{TeV}\), and by the CMS experiment [21] in \(pp\) collisions at \(\sqrt{s}=7\,\mathrm{TeV}\), using the \(\psi {(2S)} \!\rightarrow \mu ^+ \mu ^- \) decay. The CDF (CMS) measurement used \(\psi {(2S)}\) mesons in the kinematic range \(5<p_\mathrm{T} <30{\mathrm {\,GeV\!/}c} \) (\(14<p_\mathrm{T} <50{\mathrm {\,GeV\!/}c} \)) and rapidity \(|y|<0.6\) (\(|y|<1.5\)). The CDF result for \(p_\mathrm{T} >10{\mathrm {\,GeV\!/}c} \) is in strong disagreement with the NRQCD prediction of large transverse polarisation. At CMS, no evidence of large transverse or longitudinal \(\psi {(2S)}\) polarisation has been observed.
This paper presents the measurement of the prompt \(\psi {(2S)}\) polarisation in \(pp\) collisions at \(\sqrt{s} =7\,\mathrm{TeV}\), using data corresponding to an integrated luminosity of 1\(\,\text{ fb }^{-1} \), from \(\psi {(2S)} \!\rightarrow \mu ^+ \mu ^- \) decays. The \(\psi {(2S)}\) polarisation parameters are determined using unbinned maximum likelihood fits to the two-dimensional angular distribution of the \(\mu ^+ \) in the helicity and Collins-Soper frames. The measurement is performed in the \(\psi {(2S)}\) kinematic range \(3.5<p_\mathrm{T} <15{\mathrm {\,GeV\!/}c} \) and \(2.0<y<4.5\).
2 LHCb detector and data sample
The LHCb detector [23] is a single-arm forward spectrometer covering the pseudorapidity range \(2<\eta <5\), designed for the study of particles containing \(b \) or \(c \) quarks. The detector includes a high-precision tracking system consisting of a silicon-strip vertex detector surrounding the \(pp\) interaction region, a large-area silicon-strip detector located upstream of a dipole magnet with a bending power of about \(4\mathrm{\,Tm}\), and three stations of silicon-strip detectors and straw drift tubes placed downstream. The combined tracking system provides a momentum measurement with relative uncertainty that varies from 0.4 % at 5\({\mathrm {\,GeV\!/}c}\) to 0.6 % at 100\({\mathrm {\,GeV\!/}c}\), and impact parameter resolution of 20\({\,\upmu \mathrm m}\) for tracks with large transverse momentum. Different types of charged hadrons are distinguished by information from two ring-imaging Cherenkov detectors [24]. Photon, electron and hadron candidates are identified by a calorimeter system consisting of scintillating-pad and preshower detectors, an electromagnetic calorimeter and a hadronic calorimeter. Muons are identified by a system composed of alternating layers of iron and multiwire proportional chambers [25].
The trigger [26] consists of a hardware stage, based on information from the calorimeter and muon systems, followed by a software stage, which applies full event reconstruction. The hardware trigger requires the \(p_\mathrm{T}\) of one muon candidate to be larger than \(1.48{\mathrm {\,GeV\!/}c} \), or the product of the transverse momenta of two muon candidates to be larger than \(1.68{\,(\mathrm {GeV\!/}c)^2} \). In a first stage of the software trigger, two oppositely charged muon candidates with \(p_\mathrm{T} >0.5{\mathrm {\,GeV\!/}c} \) and momentum \(p>6{\mathrm {\,GeV\!/}c} \) are selected and their invariant mass is required to be greater than \(2.7{\mathrm {\,GeV\!/}c^2} \). In a second stage of the software trigger, two muon candidates consistent with originating from a \(\psi {(2S)}\) decay are chosen and additional criteria are applied to refine the sample of the \(\psi {(2S)}\) candidates as follows. The invariant mass of the candidate is required to be consistent with the known \(\psi {(2S)}\) mass [27] and, for 0.7\(\,\text{ fb }^{-1} \) of data, the \(p_\mathrm{T}\) of the candidate is required to be greater than 3.5\({\mathrm {\,GeV\!/}c}\).
In the simulation, \(pp\) collisions are generated using Pythia [28] with a specific LHCb configuration [29]. Decays of hadronic particles are described by EvtGen [30], in which final state radiation is generated using Photos [31]. The interaction of the generated particles with the detector and its response are implemented using the Geant4 toolkit [32, 33] as described in Ref. [34]. The prompt charmonium production is simulated in Pythia according to the leading order colour-singlet and colour-octet mechanisms [29, 35], and the charmonium is generated without polarisation.
3 Event selection
The \(\psi {(2S)}\) candidates are reconstructed from pairs of good quality, oppositely charged particles that originate from a common vertex. The \(\chi ^2\) probability of the vertex fit must be larger than 0.5 %. The transverse momentum of each particle is required to be greater than \(1{\mathrm {\,GeV\!/}c} \). Both tracks must also be consistent with the muon hypothesis. As in the measurement of \({J /\psi } \) polarisation [22], the significance \(S_\tau \), which is defined as the reconstructed pseudo-decay time \(\tau \) divided by its uncertainty, is used to distinguish between prompt \(\psi {(2S)}\) mesons and those from \(b\)-hadron decays. The pseudo-decay time \(\tau \) is defined as [17]
where \(z_{\psi {(2S)}}\) (\(z_{\mathrm {PV}})\) is the position of the \(\psi {(2S)} \) decay vertex (the associated primary vertex) in the \(z\)-direction, \(M_{\psi {(2S)}}\) is the known \(\psi {(2S)} \) mass, and \(p_z\) is the measured \(z\) component of the \(\psi {(2S)} \) momentum in the centre-of-mass frame of the \(pp\) collision. The \(z\)-axis of the LHCb coordinate system is defined as the beam direction in the LHCb detector region. The \(\psi {(2S)}\) mesons from \(b\)-hadron decays tend to have large values of \(S_\tau \). The requirement \(S_\tau <4\) reduces the fraction of the selected non-prompt \(\psi {(2S)}\) mesons from about 20 to 3 % [17, 22].
The analysis is performed in five \(p_\mathrm{T}\) and five \(y\) bins of the \(\psi {(2S)}\) meson. As an example, the invariant mass distribution of \(\psi {(2S)}\) candidates for \(5<p_\mathrm{T} <7{\mathrm {\,GeV\!/}c} \) and \(3.0<y<3.5\) is shown in Fig. 1. In each kinematic bin, the mass distribution is fitted with a combination of two Crystal Ball (CB) functions [36] with a common peak position for the signal and a linear function for the combinatorial background. The relative fractions of the narrower and broader CB functions are fixed to 0.9 and 0.1, respectively, determined from simulation.
Using the results of the fit to the mass distribution, the sWeight \(w_i\) for each candidate \(i\) to be signal is computed by means of the sPlot technique [37]. The correlation between the invariant mass of the \(\psi {(2S)}\) candidates and the muon angular variables is found to be negligible, and the sWeights are used to subtract the background from the angular distribution.
4 Polarisation fit
The polarisation parameters are determined from a fit to the (\(\cos \theta ,\phi \)) angular distribution of the \(\psi {(2S)} \!\rightarrow \mu ^+ \mu ^- \) signal candidates in each kinematic bin of the \(\psi {(2S)}\) meson independently. The angular distribution described by Eq. 1 is modified by the detection efficiency \(\epsilon \), which varies as a function of the angular variables (\(\cos \theta ,\phi \)). In each kinematic bin, \(\epsilon \) is obtained from a sample of simulated unpolarised \(\psi {(2S)} \!\rightarrow \mu ^+ \mu ^- \) decays, where \(\cos \theta \) and \(\phi \) are generated according to uniform distributions. As an example, Fig. 2 shows the efficiency in the helicity frame for \(\psi {(2S)}\) candidates in the kinematic bin \(5<p_\mathrm{T} <7{\mathrm {\,GeV\!/}c} \) and \(2.5<y<3.0\). The typical absolute efficiency is about 10 %. For smaller (larger) \(p_\mathrm{T}\) and \(y\) values, the coverage of the reconstructed muon angular variables is narrower (broader). In the regions \(\vert \cos \theta \vert \approx 1\), and \(\vert \phi \vert \approx 0\) or \(\pi \), the efficiency is lower because one of the two muons is likely to escape the LHCb detector acceptance.
Combining the angular distribution given in Eq. 1 with the efficiency, the logarithm of the likelihood function [38], in each \(p_\mathrm{T}\) and \(y\) bin, is defined as
where \(P(\cos \theta _{i},\phi _{i}\vert \lambda _{\theta },\lambda _{\theta \phi }, \lambda _{\phi }) \equiv 1+\lambda _\theta \cos ^2 \theta _{i} + \lambda _{\theta \phi }\sin 2\theta _{i} \cos \phi _{i} + \lambda _\phi \sin ^2 \theta _{i} \cos 2\phi _{i}\), \(w_i\) is the sWeight, and \(N_\mathrm {tot}\) is the number of \(\psi {(2S)} \) candidates in the data. The global factor \(\alpha \equiv \sum _{i=1}^{N_\mathrm {tot}}w_i/\sum _{i=1}^{N_\mathrm {tot}}w_i^2\) is introduced to estimate correctly the statistical uncertainty for the weighted likelihood function. The normalisation \(N(\lambda _{\theta }, \lambda _{\theta \phi }, \lambda _{\phi })\) is defined as
where the sum extends over the \(M_\mathrm {tot}\) candidates in the simulated and reconstructed sample and \(C\) is a constant factor. The last equality holds because the \((\cos \theta ,\phi )\) two-dimensional distribution for the fully simulated unpolarised \(\psi {(2S)}\) mesons is the same as the efficiency \(\epsilon (\cos \theta ,\phi )\) up to a constant global factor.
The angular efficiency is validated in data by using muons from \(B ^+ \!\rightarrow {J /\psi } K ^+ \) decays. Due to angular momentum conservation, the \({J /\psi }\) meson produced in this channel is longitudinally polarised in the \(B ^+ \) meson rest frame. After reweighting the kinematic properties of the simulated \(B ^+ \) and \({J /\psi }\) mesons to reproduce the data, the remaining differences of the angular distributions between the reweighted simulation sample and the data are attributed to imperfections in the modelling of the detector response. Figure 3 compares the \(\cos \theta \) distributions in data for \(B ^+ \!\rightarrow {J /\psi } K ^+ \) candidates in the helicity frame with simulated data after reweighting. The efficiency for simulated events is overestimated for \({J /\psi }\) candidates with \(|\cos \theta |>0.5\), therefore it is corrected further as a function of \(p_\mu \) and \(y_\mu \), the momentum and the rapidity of the muon in the centre-of-mass frame of \(pp\) collisions. A table of weights (corrections) in bins of \(p_\mu \) and \(y_\mu \) are determined by studying the two-dimensional \((p_{\mu },y_{\mu })\) distribution of \(B ^+ \!\rightarrow {J /\psi } K ^+ \) candidates in data and simulation. The normalisation of Eq. 5 is calculated by assigning a weight to each candidate as the product of the weights for \(\mu ^+\) and \(\mu ^-\) depending on their \(p_\mu \) and \(y_\mu \) bins.
5 Systematic uncertainties
Sources of systematic uncertainty are considered for each of the four observables \(\lambda _\theta \), \(\lambda _{\theta \phi } \), \(\lambda _{\phi } \) and \(\lambda _{\mathrm {inv}} \) in both the Collins-Soper and helicity frames. In the Collins-Soper frame, the overall systematic uncertainties are found to be comparable for each of these observables in most kinematic bins, while for the helicity frame the systematic uncertainties assigned to \(\lambda _{\theta \phi } \) and \(\lambda _{\phi } \) are typically a factor of 2–3 smaller than those estimated for \(\lambda _\theta \) and \(\lambda _{\mathrm {inv}} \). For each of the main sources of systematic uncertainty, Table 1 shows the range of values assigned over all kinematic bins, and their average. The total systematic uncertainties for each of the four observables can be found in Tables 2 and 3.
The dominant systematic uncertainty is due to the size of the \(B ^+ \!\rightarrow {J /\psi } K ^+ \) control sample. This leads to non-negligible statistical uncertainties in the correction factors that are applied to simulated events in bins of \(p_\mu \) and \(y_\mu \). The uncertainty on a given correction factor is estimated by varying it by one standard deviation of its statistical uncertainty, while keeping all other factors at their central values. The polarisation parameters are recalculated and the change relative to their default values is considered as the contribution from this factor to the systematic uncertainty. This procedure is repeated for all bins of \(p_\mu \) and \(y_\mu \), and the sum in quadrature of all these independent contributions is taken as the total systematic uncertainty.
The limited size of the sample of simulated events introduces an uncertainty on the normalisation \(N(\lambda _\theta , \lambda _{\theta \phi },\lambda _\phi )\), and this uncertainty is propagated to the polarisation parameters.
The uncertainty of the sWeight of each candidate used for the background subtraction is a source of uncertainty on the polarisation parameters. The effect is studied by comparing the default polarisation parameters with those obtained when varying the definition of the models used to fit the mass distributions and re-evaluating the sWeight for each candidate. Several alternative fitting models are studied, including an exponential function for the background mass distribution, only one CB function for the signal mass distribution, or shapes for signal and background mass distributions fixed to those obtained from fits to the mass distributions in sub-regions of the \((\cos \theta ,\phi )\) distribution space. The largest variation with respect to the default result is assigned as the systematic uncertainty.
In each kinematic bin, discrepancies between data and simulation in the \(\psi {(2S)}\) \(p_\mathrm{T}\) and \(y\) distributions introduce an additional uncertainty. This is evaluated by comparing the default polarisation results with those determined after the \(\psi {(2S)}\) kinematic distribution in the simulation is weighted to that in data. The difference between the two results is quoted as a systematic uncertainty contribution.
The uncertainty due to the contamination of \(\psi {(2S)}\) candidates from \(b\)-hadron decays (\(3~\%\)) is determined by relaxing the \(S_\tau \) selection and studying the variations of the polarisation parameters.
With the exception of the effects due to the differences in the \(\psi {(2S)}\) kinematic spectrum and the size of the sample of simulated events, correlations are expected among \(\psi {(2S)}\) kinematic bins. The correlation between these systematic uncertainties in adjacent bins could be as large as 50 %, as the final state muons may have similar momentum and rapidity. For each kinematic bin, the total systematic uncertainty is calculated as the quadratic sum of the various sources of systematic uncertainties assuming no correlation within each kinematic bin.
6 Results
The results for the polarisation parameters \(\lambda _\theta \), \(\lambda _{\theta \phi }\), \(\lambda _{\phi }\) and \(\lambda _{\mathrm {inv}}\), and their uncertainties, in each \(p_\mathrm{T}\) and \(y\) bin of the prompt \(\psi {(2S)}\) meson sample, are reported in Tables 2 and 3 for the helicity and the Collins-Soper frames, respectively. The systematic uncertainties are similar in size to the statistical uncertainties. The parameters \(\lambda _\theta \) and \(\lambda _{\mathrm {inv}}\) are also shown in Fig. 4 as functions of the \(p_\mathrm{T}\) of the \(\psi {(2S)}\) mesons, for different \(y\) bins.
The frame-invariant polarisation parameter \(\lambda _{\mathrm {inv}}\) is consistent with a negative polarisation with no strong dependence on the \(p_\mathrm{T}\) and \(y\) of the \(\psi {(2S)}\) meson. The values and uncertainties of \(\lambda _{\mathrm {inv}}\) that are measured in the helicity and the Collins-Soper frames are in good agreement with each other, with differences much smaller than the statistical uncertainties. In the Collins-Soper frame, \(\lambda _\theta \) takes small negative values especially in the low-\(p_\mathrm{T}\) region and increases with \(p_\mathrm{T}\). This trend is more significant for the extreme \(y\) bins. In the helicity frame, the polarisation parameter \(\lambda _\theta \) is consistent with zero, with no significant dependence on \(p_\mathrm{T}\) or \(y\) of the \(\psi {(2S)}\) meson. The polarisation parameters \(\lambda _{\theta \phi }\) and \(\lambda _{\phi }\) are consistent with zero in both the helicity and Collins-Soper frames, and their absolute values are below 0.1 for most of the kinematic bins.
In Fig. 5, the measured values of \(\lambda _\theta \) in the helicity frame as a function of \(p_\mathrm{T}\) of the \(\psi {(2S)}\) meson, integrating over the rapidity range \(2.5<y<4.0\), are compared with the predictions of the CSM [39] and NRQCD [39–41] models at NLO. Our results disfavour the CSM calculations, in which the \(\psi {(2S)}\) meson is significantly longitudinally polarised. The three NRQCD calculations in Refs. [39–41], which use different selections of experimental data to determine the non-perturbative matrix elements, provide a good description of our measurements in the low-\(p_\mathrm{T}\) region. However, the prediction of increasing polarisation with \(p_\mathrm{T}\) in these models is not supported by the LHCb data.
7 Conclusion
The polarisation of prompt \(\psi {(2S)}\) mesons is measured as a function of the \(\psi {(2S)}\) \(p_\mathrm{T}\) and \(y\) in the range \(3.5<\) \(p_\mathrm{T} <15{\mathrm {\,GeV\!/}c} \) and \(2.0<y<4.5\), in \(pp\) collisions at \(\sqrt{s}= 7\,\mathrm{TeV}\). The analysis is performed using data corresponding to an integrated luminosity of 1.0 \(\,\text{ fb }^{-1} \), collected by the LHCb experiment in 2011. The polarisation parameters \(\lambda _\theta \), \(\lambda _{\theta \phi }\), \(\lambda _{\phi }\) and \(\lambda _{\mathrm {inv}}\) are determined in the helicity and Collins-Soper frames by studying the angular distribution of the two muons produced in the \(\psi {(2S)} \!\rightarrow \mu ^+ \mu ^- \) decay.
The frame-independent observable \(\lambda _{\mathrm {inv}}\) is consistent with a negative polarisation. The measured values of \(\lambda _{\theta \phi }\) and \(\lambda _{\phi }\) are small over the accessible kinematic range. The \(\lambda _\theta \) distribution in the helicity frame shows that the \(\psi {(2S)}\) meson exhibits neither large transverse nor longitudinal polarisation. Although a direct comparison with previous measurements by CMS and CDF is not possible due to the different kinematic ranges, all results disfavour large polarisation in the high-\(p_\mathrm{T}\) region. The prompt \(\psi {(2S)}\) polarisation measured at LHCb disagrees with the CSM predictions both in the size of the polarisation parameters and the \(p_\mathrm{T} \) dependence. While the NRQCD models provide a good description of the LHCb data in the low-\(p_\mathrm{T}\) region, the predicted transverse polarisation at high-\(p_\mathrm{T}\) is not observed.
Notes
For a \(\psi {(2S)} \) meson in a pure spin state the three polarisation parameters cannot vanish simultaneously.
References
N. Brambilla et al., Heavy quarkonium: progress, puzzles, and opportunities. Eur. Phys. J. C 71, 1534 (2011). arXiv:1010.5827
K. Gottfried, J.D. Jackson, On the connection between production mechanism and decay of resonances at high energies. Nuovo Cim. 33, 309 (1964)
P. Faccioli, C. Lourenço, J. Seixas, H.K. Wöhri, Towards the experimental clarification of quarkonium polarization. Eur. Phys. J. C 69, 657 (2010). arXiv:1006.2738
M. Jacob, G.C. Wick, On the general theory of collisions for particles with spin. Ann. Phys. 7, 404 (1959)
J.C. Collins, D.E. Soper, Angular distribution of dileptons in high-energy hadron collisions. Phys. Rev. D 16, 2219 (1977)
P. Faccioli, C. Lourenço, J. Seixas, H.K. Wöhri, Rotation-invariant observables in parity-violating decays of vector particles to fermion pairs. Phys. Rev. D 82, 96002 (2010). arXiv:1010.1552
P. Faccioli, C. Lourenço, J. Seixas, Rotation-invariant relations in vector meson decays into fermion pairs. Phys. Rev. Lett. 105, 61601 (2010). arXiv:1005.2601
J.P. Lansberg, \(J/\psi \) production at \(\sqrt{s}=1.96\) and 7 Tev: color-singlet model, NNLO* and polarization. J. Phys. G38, 124110 (2011). arXiv:1107.0292
G.T. Bodwin, E. Braaten, G.P. Lepage, Rigorous QCD analysis of inclusive annihilation and production of heavy quarkonium. Phys. Rev. D 51, 1125 (1995). arXiv:hep-ph/9407339
G.T. Bodwin, E. Braaten, G.P. Lepage, Erratum: Rigorous QCD analysis of inclusive annihilation and production of heavy quarkonium. Phys. Rev. D 55, 5853 (1997)
M. Beneke, M. Krämer, Direct \(J/\psi \) and \(\psi ^\prime \) polarization and cross-sections at the Fermilab Tevatron. Phys. Rev. D 55, 5269 (1997). arXiv:hep-ph/9611218
M. Cacciari, M. Greco, M.L. Mangano, A. Petrelli, Charmonium production at the Tevatron. Phys. Lett. B 356, 553 (1995). arXiv:hep-ph/9505379
P. Cho, A.K. Leibovich, Color-octet quarkonia production. Phys. Rev. D 53, 150 (1996). arXiv:hep-ph/9505329
E. Braaten, S. Fleming, Color-octet fragmentation and the \(\psi ^{\prime }\) surplus at the Fermilab Tevatron. Phys. Rev. Lett. 74, 3327 (1995). arXiv:hep-ph/9411365
CDF collaboration, F. Abe et al., Inclusive \(J/\psi \), \(\psi (2S)\), and \(b\)-quark production in \(\bar{p}p\) collisions at \(\sqrt{s} = 1.8\,{\rm TeV}\), Phys. Rev. Lett. 69, 3704 (1992)
LHCb collaboration, R. Aaij et al., Measurement of \(J/\psi \) production in \(pp\) collisions at \(\sqrt{s}=7\,{\rm TeV}\). Eur. Phys. J. C 71, 1645 (2011). arXiv:1103.0423
LHCb collaboration, R. Aaij et al., Measurement of \(\psi (2S)\) meson production in \(pp\) collisions at \(\sqrt{s}=7\,{\rm TeV}\). Eur. Phys. J. C 72 2100 (2012). arXiv:1204.1258
CMS collaboration, S. Chatrchyan et al., \(J/\psi \) and \(\psi (2S)\) production in \(pp\) collisions at \(\sqrt{s}=7\,{\rm TeV}\). JHEP 02 (2012). arXiv:1111.1557
CDF collaboration, A. Abulencia et al., Polarization of \(J/\psi \) and \(\psi (2S)\) mesons produced in \(p \bar{p}\) collisions at \(\sqrt{s} = 1.96\,{\rm TeV}\). Phys. Rev. Lett. 99, 132001 (2007). arXiv:0704.0638
ALICE collaboration, B. Abelev et al, \(J/\psi \) polarization in \(pp\) collisions at \(\sqrt{s}=7\,{\rm TeV}\). Phys. Rev. Lett. 108, 082001 (2012). arXiv:1111.1630
CMS collaboration, S. Chatrchyan et al., Measurement of the prompt \(J/\psi \) and \(\psi (2S)\) polarizations in pp collisions at \(\sqrt{s} = 7\,{\rm TeV}\). Phys. Lett. B 727, 381 (2013). arXiv:1307.6070
LHCb collaboration, R. Aaij et al., Measurement of \(J/\psi \) polarization in \(pp\) collisions at \(\sqrt{s}=7 \,{\rm TeV}\). Eur. Phys. J. C 73, 2631 (2013). arXiv:1307.6379
LHCb collaboration, A.A. Alves Jr. et al, The LHCb detector at the LHC. JINST 3, S08005 (2008)
M. Adinolfi et al., Performance of the LHCb RICH detector at the LHC. Eur. Phys. J. C 73, 2431 (2013). arXiv:1211.6759
A.A. Alves Jr et al., Performance of the LHCb muon system. JINST 8, P02022 (2013). arXiv:1211.1346
R. Aaij et al., The LHCb trigger and its performance in 2011. JINST 8, P04022 (2013). arXiv:1211.3055
Particle Data Group, J. Beringer et al., Review of particle physics. Phys. Rev. D 86, 010001 (2012), and 2013 partial update for the 2014 edition
T. Sjöstrand, S. Mrenna, P. Skands, PYTHIA 6.4 physics and manual, JHEP 05, 026 (2006). arXiv:hep-ph/0603175
I. Belyaev et al., in Handling of the generation of primary events in Gauss, the LHCb simulation framework, Nuclear Science Symposium Conference Record (NSS/MIC) (IEEE, New York, 2010), p. 1155
D.J. Lange, The EvtGen particle decay simulation package. Nucl. Instrum. Meth. A 462, 152 (2001)
P. Golonka, Z. Was, PHOTOS Monte Carlo: a precision tool for QED corrections in Z and W decays. Eur. Phys. J. C 45, 97 (2006). arXiv:hep-ph/0506026
Geant4 collaboration, J. Allison et al., Geant4 developments and applications. IEEE Trans. Nucl. Sci. 53, 270 (2006)
Geant4 collaboration, S. Agostinelli et al., Geant4: a simulation toolkit. Nucl. Instrum. Meth. A 506, 250 (2003)
M. Clemencic et al., The LHCb simulation application, Gauss: design, evolution and experience. J. Phys. Conf. Ser. 331, 032023 (2011)
M. Bargiotti, V. Vagnoni, Heavy quarkonia sector in PYTHIA 6.324: Tuning, validation and perspectives at LHC(b), CERN-LHCB-2007-042
T. Skwarnicki, A study of the radiative cascade transitions between the Upsilon-prime and Upsilon resonances, PhD thesis, Institute of Nuclear Physics, Krakow, 1986, DESY-F31-86-02
M. Pivk, F.R. Le Diberder, sPlot: a statistical tool to unfold data distributions. Nucl. Instrum. Meth. A 555, 356 (2005). arXiv:physics/0402083
Y. Xie, sFit: a method for background subtraction in maximum likelihood fit, arXiv:0905.0724
M. Butenschoen, B.A. Kniehl, \(J/\psi \) polarization at the Tevatron and the LHC: nonrelativistic-QCD factorization at the crossroads. Phys. Rev. Lett. 108, 172002 (2012). arXiv:1201.1872
B. Gong, L.-P. Wan, J.-X. Wang, H.-F. Zhang, Polarization for prompt \(J/\psi \) and \(\psi (2S)\) production at the Tevatron and LHC. Phys. Rev. Lett. 110, 042002 (2013). arXiv:1205.6682
K.-T. Chao et al., \(J/\psi \) polarization at hadron colliders in nonrelativistic QCD. Phys. Rev. Lett. 108, 242004 (2012). arXiv:1201.2675
Acknowledgments
We wish to thank M. Butenschön, B. Gong, H.-S. Shao and Y.-Q. Ma for providing us with the theoretical calculations and helpful discussions. We express our gratitude to our colleagues in the CERN accelerator departments for the excellent performance of the LHC. We thank the technical and administrative staff at the LHCb institutes. We acknowledge support from CERN and from the national agencies: CAPES, CNPq, FAPERJ and FINEP (Brazil); NSFC (China); CNRS/IN2P3 and Region Auvergne (France); BMBF, DFG, HGF and MPG (Germany); SFI (Ireland); INFN (Italy); FOM and NWO (The Netherlands); SCSR (Poland); MEN/IFA (Romania); MinES, Rosatom, RFBR and NRC “Kurchatov Institute” (Russia); MinECo, XuntaGal and GENCAT (Spain); SNSF and SER (Switzerland); NAS Ukraine (Ukraine); STFC (United Kingdom); NSF (USA). We also acknowledge the support received from EPLANET and the ERC under FP7. The Tier1 computing centres are supported by IN2P3 (France), KIT and BMBF (Germany), INFN (Italy), NWO and SURF (The Netherlands), PIC (Spain), GridPP (United Kingdom). We are indebted to the communities behind the multiple open source software packages on which we depend. We are also thankful for the computing resources and the access to software R&D tools provided by Yandex LLC (Russia).
Author information
Authors and Affiliations
Consortia
Corresponding author
Rights and permissions
Open Access This article is distributed under the terms of the Creative Commons Attribution License which permits any use, distribution, and reproduction in any medium, provided the original author(s) and the source are credited.
Funded by SCOAP3 / License Version CC BY 4.0.
About this article
Cite this article
The LHCb Collaboration., Aaij, R., Adeva, B. et al. Measurement of \(\psi {(2S)} \) polarisation in \(pp\) collisions at \(\sqrt{s}\) = 7 TeV. Eur. Phys. J. C 74, 2872 (2014). https://doi.org/10.1140/epjc/s10052-014-2872-9
Received:
Accepted:
Published:
DOI: https://doi.org/10.1140/epjc/s10052-014-2872-9