EPJ Web of Conferences 6, 55 () DOI:.5/epjconf/655 Owned by the authors, published by EDP Sciences, Quarkonium Production at Monica Pepe Altarelli,a On behalf of the collaboration CERN Abstract. I will review a selection of results on the production of heavy quarkonium states in pp collisions, including recent results on J/ψ and Υ(nS ) (n =,, ) production at s = 8 ev, as well as preliminary results on J/ψ production in proton-lead collisions at s NN = 5 ev. Introduction Heavy quarkonium states are particularly interesting systems to test our understanding of strong interactions, both at perturbative and non-perturbative level, and are therefore the subject of intense theoretical and experimental studies. An effective field theory, non-relativistic QCD (NRQCD) [, ], provides the foundation for much of the current theoretical work. According to NRQCD, the production of heavy quarkonium factorises into two steps: a heavy quark-antiquark pair is first created at short distances and subsequently evolves non-perturbatively into quarkonium at long distances. he NRQCD calculations include colour-singlet (CS) and colour-octet (CO) amplitudes (see [] and references therein), which account for the probability of a heavy quark-antiquark pair in a particular colour state to evolve into a heavy quarkonium state. he CS model (CSM) was initially used to describe experimental data. However, it underestimates the observed cross-section for single J/ψ production at high transverse momentum (p ). o resolve this discrepancy, the CO mechanism was introduced. More recent higherorder calculations to the CS predictions raise substantially the cross-sections at large p bringing them closer to the experimental data. However, none of these approaches can reproduce in a consistent way the available experimental results on both cross-section and polarisation (see [] and references therein, [4]). Heavy charmonium is also produced from b-hadron decays. his production mechanism can be distinguished from prompt quarkonium production by exploiting the b- hadron finite decay time. In this case QCD predictions are based on the Fixed-Order-Next-to-Leading-Log (FONLL) formalism [5, 6], which combines a Fixed-Order calculation with an all-order resummation of Leading and Nextto-Leading logarithms. In the following, results from heavy quarkonium production at will be compared with various theoretical a e-mail: monica.pepe.altarelli@cern.ch models to provide direct tests of the underlying production mechanism. Furthermore, preliminary results on J/ψ production in proton-lead collisions at s NN = 5 ev will be reported. Production of J/ψ and Υ(nS ) mesons he differential production cross-sections of prompt J/ψ and Υ(nS ) mesons produced at the pp collision point either directly or via feed-down from higher mass charmonium or bottomonium states, are measured at s = 8 ev [7] with the detector [8]. he measurements are performed in the range of rapidity. < y < 4.5 and p < 4 GeV/c for the J/ψ meson and p < 5 GeV/c for the Υ(nS ), in bins of p and y with bin sizes p = GeV/c and y =.5. he fraction of J/ψ from b (abbreviated as from b in the following) is also measured in the same fiducial region. he Υ(nS ) meson analysis is based on a data sample, corresponding to an integrated luminosity of about 5 pb of pp interactions, while the analysis for the more abundant J/ψ mesons is based on data, corresponding to an integrated luminosity of about 8 pb. he J/ψ and Υ(nS ) mesons are reconstructed through their decay into a pair of muons following the strategy described in Refs. [7, 9, ]. he resulting dimuon invariant mass distributions are displayed in Figs. and for the J/ψ and Υ(nS ) mesons. he invariant mass resolution is on average.5 MeV/c for the J/ψ and 4 MeV/c for the Υ(S ) meson. he production is studied under the assumption of zero polarisation, and no corresponding systematic uncertainty is assigned on the cross-section for this effect. his assumption is justified by the small polarisation measured for the J/ψ meson at s = 7 ev by [] and AL- ICE [], in a kinematic range similar to that used in this analysis, and for the Υ(nS ) by CMS [] at large p and central rapidity. he measured differential cross-sections for the production of prompt J/ψ mesons as a function of p are his is an Open Access article distributed under the terms of the Creative Commons Attribution License., which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited. Article available at http://www.epj-conferences.org or http://dx.doi.org/.5/epjconf/655
EPJ Web of Conferences Candidates / (5 MeV/c ) 5 5 (a) s = 8 ev.5<y <. <p <4 GeV/c [nb/(gev/c] d σ(j/ψ) d p 5 prompt J/ψ,. < y < 4.5 4 Direct NLO CSM,. < y < 4.5 Direct NLO NRQCD,. < y < 4.5 Direct NNLO*,. < y < 4.5 5 5 m( µ + µ ) [MeV/c ] Figure. Dimuon invariant mass of the selected J/ψ µ + µ candidates in a selected bin in y and p (.5 < y <. and < p < 4 GeV/c). Candidates / (5 MeV/c ) s = 8 ev 8 6 4 9 m( µ + µ ) [MeV/c ] Figure. Invariant mass distribution of the selected Υ(nS ) µ + µ candidates in the range p < 5 GeV/c and. < y < 4.5. he three peaks correspond to the Υ(S ), Υ(S ) and Υ(S ) meson signals (from left to right). compared in Fig. to three theoretical models: an NRQCD model at next-to-leading order (NLO) [4, 5], an NNLO* CSM [6, 7] where the notation NNLO* indicates that the calculation at next-to-next leading order is not complete and neglects part of the logarithmic terms and an NLO CSM [8]. In these comparisons the predictions are for direct J/ψ meson production, i.e., they do not include feed-down from higher charmonium states, whereas the experimental measurements do include feed-down. In particular, the contribution from J/ψ mesons produced in radiative χ c decays in the considered fiducial range was measured to be at the level of % at s = 7 ev [9]. Allowing for this contribution, both the NNLO* CSM and the NLO NRQCD models provide reasonable descriptions of the experimental data. In contrast, the CSM at NLO underestimates the cross-section by an order of magnitude. he fraction of J/ψ from b is shown is Fig.4 as a function of the J/ψ transverse momentum in bins of rapidity. his fraction is measured to be smaller at low p and for forward data. In Fig.5 the data for the differential production cross-section for J/ψ from b as a function of p at s = 8 ev are compared to the FONLL predictions. Very good agreement is ob- - s = 8 ev 5 Figure. Comparison of the differential cross-section for the production of prompt J/ψ meson (under the assumption of zero polarisation) as a function of p with direct production in an NLO NRQCD model [4, 5] (orange diagonal shading), an NNLO* CSM [7] (solid yellow) and an NLO CSM [8] (blue vertical shading). Fraction of J/ψ from b.4... s = 8 ev.5<y<..<y<.5.<y<.5.5<y<4. 4.<y<4.5 5 Figure 4. Fraction of J/ψ from b as a function of the J/ψ transverse momentum, in bins of rapidity. served. he integrated cross-section for J/ψ from b production in the defined fiducial region is determined to be σ (J/ψ from b, p < 4 GeV/c,. < y < 4.5) =.8 ±. ±. µb [7], where the first uncertainty is statistical and the second is systematic, in very good agreement with the FONLL prediction of.4 +.6.49 µb [6]. he results on J/ψ production at s = 8 ev provide an extra measurement, which complements those already published by at s =.76 [] and 7 ev[9]. Figure 6 shows the results for the cross-section for J/ψ from b at s =.76, 7 and 8 ev as well as the FONLL prediction as a function of centre-of-mass energy. Excellent agreement with the theoretical calculation is observed. From the measurement of the production cross-section for J/ψ from b in the fiducial region, using the Monte Carlo simulation based on Pythia to extrapolate to the full polar angle range (with an extrapolation factor 55-p.
LHCP [nb/(gev/c] d σ(j/ψ) d p J/ψ from b,. < y < 4.5 (a) FONLL,. < y < 4.5 [nb/(gev/c)] dσ S dp S B (a) ϒ(S) data,.<y <4.5 Direct NNLO* CSM,.<y <4.5 Direct NLO CSM,.<y <4.5 - - s = 8 ev 5 5 Figure 5. Differential production cross-section for J/ψ from b as a function of p in the fiducial range. < y < 4.5. he FONLL prediction [5, 6] is shown in yellow. [nb/(gev/c)] dσ S dp - -4 s = 8 ev 5 5 (b) ϒ(S) data,.<y <4.5 Direct NNLO* CSM,.<y <4.5 Direct NLO CSM,.<y <4.5 - [µb] 4.5.5.5.5 J/ from b,. < y < 4.5, p FONLL,. < y < 4.5, p < 4 GeV/c < 4 GeV/c 5 5 s [ev] S B [nb/(gev/c)] dσ S dp S B - -4 s = 8 ev 5 5 (c) ϒ(S) data,.<y <4.5 Direct NNLO* CSM,.<y <4.5 Direct NLO CSM,.<y <4.5 - Figure 6. Predictions based on the FONLL formalism [5, 6] for the production cross-section for J/ψ from b in the fiducial range < p < 4 GeV/c and. < y < 4.5 (yellow band). he uncertainty includes contributions from the renormalisation scale, quark masses and the choice of PDF set. he black dotted line shows the central value of the prediction. he points show the measurements at s =.76 [], 7 [9], and 8 ev [7]. of 5.4 from the measured cross-section to the full kinematic region) and the inclusive b J/ψ X branching fraction, one can derive the total bb production cross-section at s = 8 ev, σ(pp bbx) = 98 ± ± 6 µb, where the first uncertainty is statistical and the second is systematic. In Fig.7 the cross-sections times dimuon branching fractions for the three Υ meson states are compared to the NNLO CSM [6] and an NLO CSM [8] predictions as a function of p. he NNLO* CSM provides a reasonable description of the experimental data, particularly for the Υ(S ) meson, which is expected to be less affected by feed-down. As for the prompt J/ψ meson production, the CSM at NLO underestimates the cross-section by an order of magnitude. Production of ψ(s ) mesons he production of ψ(s ) mesons was studied by in pp collisions at 7 ev []. he ψ(s ) mesons were - -4 s = 8 ev 5 5 Figure 7. Comparison of the differential production crosssections times dimuon branching fractions for (a) Υ(S ), (b) Υ(S ) and (c) Υ(S ) mesons as a function of p with direct production in an NNLO CSM [6] (solid yellow) and an NLO CSM [8] (blue vertical shading) model. reconstructed in the decay channel ψ(s ) µ + µ and ψ(s ) J/ψπ + π and the results from the two channels were combined. he ψ(s ) production from b-hadron decays is distinguished from promptly produced charmonium by exploiting the b-hadron finite decay time. he interest of measuring ψ(s ) production arises from the fact that no significant feed-down from excited charmonium states is expected, which facilitates the comparison with theoretical model. Figures 8 and 9 show the differential production cross-section as a function of p for prompt ψ(s ) and ψ(s ) from b-hadron-decays. he data are compared with various theoretical models: MWC [] and KB [5] are NLO NRQCD calculations, while AL [6, 7] is a CSM including the dominant NNLO terms. he NLO NRQCD and FONLL [5, 6] calculations provide a good description of the ψ(s ) prompt cross-section and of the cross-section of the ψ(s ) from b-hadron decays. 55-p.
EPJ Web of Conferences [htb] Figure 8. Differential production cross-section as a function of p for prompt ψ(s ). he data are compared with various theoretical models: MWC [] and KB [5] are NLO NRQCD calculations, while AL [6, 7] is a CSM including the dominant NNLO terms Figure. Differential production cross-section for J/ψ pairs as a function of the invariant mass of the J/ψ pair system. he shaded area corresponds to the prediction by the model described in Ref. [4]. tion of the invariant mass of the J/ψ pair system. he data are compared with a theoretical prediction [4], which includes both direct production and feed-down from ψ(s ) decays and no contribution from DPS. Within the available statistics good agreement is observed. For double J/ψ production, the predictions from single and double parton scattering are expected to be fairly close in magnitude [5], but with somewhat different kinematic distributions for the J/ψ pairs. he analysis of the higher statistics already collected by will allow an in-depth study of the kinematic properties of double J/ψ production and of the different production models. Figure 9. Differential production cross-section as a function of p for ψ(s ) from b-hadron decays. he data are compared with FONLL calculations [5, 6]. 4 Double J/ψ production Given the large charmonium production cross-section at the LHC, the question of multiple production naturally arises. his was studied by using an integrated luminosity of about 7 pb of pp interactions at 7 ev []. Simultaneous production of two J/ψ mesons can be expected in single parton collisions via higher order O(α 4 s) gluon-gluon diagrams. his production could also be enhanced by double parton scattering (DPS) in which two gluons are extracted from a proton and the two J/ψ mesons are produced in two independent sub-processes. he cross-section for double J/ψ production in the fiducial range < y J/ψ < 4.5 and p J/ψ < GeV/c is measured to be σ J/ψJ/ψ = 5. ±. ±.nb [], where the first uncertainty is statistical and the second is systematic, in good agreement with predictions from leading-order QCD calculations [4]. he suppression factor with respect to single prompt J/ψ production is σ J/ψJ/ψ /σ J/ψ = (5. ±. ±.6 +.. ) 4, where the first uncertainty is statistical, the second systematic and the third accounts for the unknown J/ψ polarisation. Figure shows the differential production cross-section for J/ψ pairs as a func- 5 Exclusive production of J/ψ and ψ(s ) mesons Exclusive production of J/ψ and ψ(s ) mesons through processes such as the one displayed in Fig. were observed by [6] in the dimuon channel. is very well suited for this type of studies as it has access to high rapidities, with some sensitivity to backward tracks using the Vertex Locator (VELO), it operates at relatively low pileup and has sensitivity to low momentum and transverse momentum particles. Exclusively produced J/ψ and ψ(s ) mesons into dimuons are characterised by a very distinct topology: their selection requires the presence of two reconstructed muons in the forward region and no other tracks and photons in the detector. An additional rapidity gap is obtained by also excluding tracks in the backward region using the VELO detector. Based on a data sample, corresponding to an integrated luminosity of about 6 pb of pp interactions, the cross-sections times branching fractions to two muons for exclusive J/ψ and ψ(s ) with pseudorapidities between. and 4.5 is measured to be σ pp J/ψ( µ + µ ) = 7 ± ± 6 pb and σ pp ψ(s )( µ + µ ) = 7.8 ±. ±. pb, where the first uncertainty is statistical and the second is systematic, in good agreement with theoretical predictions [7 ]. During the run, a new trigger was implemented to increase the rate of exclusive events, which makes use of upstream silicon sensors to veto any backward activity 55-p.4
LHCP Figure. Feynman diagram contributing to exclusive J/ψ photoproduction. σ prompt J/ψ /A [µb] 9 8 7 6 5 4 Preliminary p < 4 GeV/c.5 < y < 4. (pp,pa) 4. < y <.5 (Ap) pp (Rescaled) 4 6 8 s NN [ev] Ap pa and of soft p cuts to gain sensitivity for hadronic final states. Data corresponding to an integrated luminosity of about.4 fb were already collected with such a trigger, which will allow to study the production of charmonium decaying into hadronic final states, the production of open charm and to perform a search for higher mass charmonium states. 6 Production of J/ψ mesons in proton-lead collisions at s NN = 5 ev has recently presented preliminary results on production of J/ψ mesons in proton-lead collisions at s NN = 5 ev [, ]. Data were collected with a proton beam energy of 4 ev and a lead beam energy of.58 ev per nucleon, resulting in a centre-of-mass energy of the proton-nucleon system of 5 ev. he directions of the proton and lead beams were swapped during data taking to produce both pa and Ap collisions. Data collected with the inverted beam directions (i.e., in Ap collision) allow to measure backward production. he protonnucleon centre-of-mass system has a rapidity in the laboratory frame of +.47 for pa and -.47 for Ap collisions, resulting in a rapidity coverage in the laboratory frame ranging from about.5 to 4. for pa collisions and from -5. to.5 for Ap collisions. he analysis is based on data samples corresponding to integrated luminosities of.75 nb of pa and. nb of Ap collisions. he double-differential production is measured for prompt J/ψ and J/ψ from b in bins of the kinematic variables y and p. he measured cross-sections for prompt J/ψ production integrated over y and p, scaled by a factor of /A and rescaled to the common rapidity range of the protonnucleon system.5 < y < 4. for pp and pa collisions and 4. < y <.5 for Ap collisions are shown is Fig. as a function of s NN. he results indicate that the J/ψ production cross-section is suppressed in protonlead collisions, as expected from previous experiments; however such suppression is less enhanced in the backward region. One can also derive the attenuation factor R pa (y, s NN ) = A dσ pa dy (y, s NN ) dσpp dy (y, s NN ) by using for the J/ψ production cross-section in pp collisions at 5 ev a linear interpolation of the published measurements at s =.76 [], 7 [9], and 8 ev [7]. he result is shown in Fig. and Figure. Integrated cross-sections for prompt J/ψ production scaled by a factor of /A and rescaled to the common rapidity range of the proton-nucleon system.5 < y < 4. for pp and pa collisions and 4. < y <.5 for Ap collisions as a function of snn R pa (y).8.6.4..8.6.4. Preliminary pa/ap = 5 ev s NN., prompt J/ψ E. loss E. loss + saturation -5 5 Figure. Attenuation factor R pa compared to theoretical predictions []. he black triangles are measurements, the red solid line is the theoretical prediction based on parton energy loss effects, the blue dashed line takes additional saturation effects into account. compared to theoretical predictions []. Within the available statistics good agreement is observed. 7 Conclusion has performed a wealth of measurements on quarkonium production, and many more are coming! A few selected ones were briefly illustrated in this report. Measurements of charmonium [9,, ] and bottomonium [7, ] production were performed at various centre-of-mass energies. hey allow an in depth comparison with theoretical models. A simple CSM is disfavoured by the data, while a combination of CS and CO, as implemented in the NRQCD formalism, or CS improved by QCD corrections, provide a good description of prompt y 55-p.5
EPJ Web of Conferences quarkonium production. Charmonium production from b- hadron decays is very well reproduced by FONLL calculations. Some other studies were illustrated in this report, such as charmonium exclusive production [6], double J/ψ production [], as well as first preliminary results on J/ψ production in proton-lead collisions at snn = 5 ev []. Many of these results are based on the analysis of a partial dataset and will gain in precision with the analysis of the full data sample available to. References [] W.E. Caswell, G.P. Lepage, Phys.Lett. B67, 47 (986) [] G.. Bodwin, E. Braaten, G.P. Lepage, Phys. Rev. D 5, 5 (995) [] N. Brambilla et al., Eur. Phys. J. C7, 54 (),.587 [4] C. Lourenco (CMS collaboration), these proceedings (LHCp ) [5] M. Cacciari, M. Greco, P. Nason, JHEP 5, 7 (998), hep-ph/984 [6] M. Cacciari, S. Frixione, N. Houdeau, M.L. Mangano, P. Nason et al., JHEP, 7 (), 5.644 [7] R. Aaij et al. ( collaboration), JHEP 6, 64 (), 4.6977 [8] A.A. Alves Jr. et al. ( collaboration), JINS, S85 (8) [9] R. Aaij et al. ( collaboration), Eur. Phys. J. C7, 645 (),.4 [] R. Aaij et al. ( collaboration), Eur.Phys.J. C7, 5 (),.6579 [] R. Aaij et al. ( collaboration) (), in preparation [] B. Abelev et al. (ALICE collaboration), Phys.Rev.Lett. 8, 8 (),.6 [] S. Chatrchyan et al. (CMS collaboration), Phys.Rev.Lett. (), 9.9 [4] M. Butenschön, B.A. Kniehl, Phys.Rev. D84, 55 (), 5.8 [5] M. Butenschön, B.A. Kniehl, Phys.Rev.Lett. 6, (), 9.566 [6] P. Artoisenet, J.M. Campbell, J. Lansberg, F. Maltoni, F. ramontano, Phys.Rev.Lett., 5 (8), 86.8 [7] J.P. Lansberg, Eur. Phy. J. C 6, 69 (8), 8.45 [8] J.M. Campbell, F. Maltoni, F. ramontano, Phys.Rev.Lett. 98, 5 (7), hep-ph/7 [9] R. Aaij et al. ( collaboration), Phys.Lett. B78, 4 (), 4.46 [] R. Aaij et al. ( collaboration), JHEP, 4 (),.45 [] R. Aaij et al. ( collaboration), Eur.Phys.J. C7, (), 4.58 [] Y.Q. Ma, K. Wang, K.. Chao, Phys.Rev. D84, 4 (),. [] R. Aaij et al. ( Collaboration), Phys.Lett. B77, 5 (), 9.96 [4] A. Berezhnoy, A. Likhoded, A. Luchinsky, A. Novoselov, Phys.Rev. D84, 94 (),.588 [5] A.V. Berezhnoy, A.K. Likhoded, A.V. Luchinsky, A.A. Novoselov, Phys. Rev. D 86, 47 () [6] R. Aaij et al. ( collaboration), J.Phys. G4, 45 (),.784 [7] L. Harland-Lang, V. Khoze, M. Ryskin, W. Stirling, Eur.Phys.J. C65, 4 (), 99.4748 [8] S.R. Klein, J. Nystrand, Phys. Rev. Lett. 9, 4 (4) [9] L. Motyka, G. Watt, Phys.Rev. D78, 4 (8), 85. [] V. Goncalves, M. Machado, Phys.Rev. C84, 9 (), 6.6 [] -CONF-8 () [] B. Schmidt ( collaboration), these proceedings (LHCp ) [] F. Arleo, S. Peigne, JHEP, (),.44 55-p.6