Universal kinematic scaling as a probe of factorized long-distance effects in high-energy quarkonium production

The European Physical Journal C, Feb 2018

Dimensional analysis reveals general kinematic scaling rules for the momentum, mass, and energy dependence of Drell–Yan and quarkonium cross sections. Their application to mid-rapidity LHC data provides strong experimental evidence supporting the validity of the factorization ansatz, a cornerstone of non-relativistic QCD, still lacking theoretical demonstration. Moreover, data-driven patterns emerge for the factorizable long-distance bound-state formation effects, including a remarkable correlation between the S-wave quarkonium cross sections and their binding energies. Assuming that this scaling can be extended to the P-wave case, we obtain precise predictions for the not yet measured feed-down fractions, thereby providing a complete picture of the charmonium and bottomonium feed-down structure. This is crucial information for quantitative interpretations of quarkonium production data, including studies of the suppression patterns measured in nucleus-nucleus collisions.

A PDF file should load here. If you do not see its contents the file may be temporarily unavailable at the journal website or you do not have a PDF plug-in installed and enabled in your browser.

Alternatively, you can download the file locally and open with any standalone PDF reader:


Universal kinematic scaling as a probe of factorized long-distance effects in high-energy quarkonium production

Eur. Phys. J. C Universal kinematic scaling as a probe of factorized long-distance effects in high-energy quarkonium production Pietro Faccioli Carlos Lourenço Mariana Araújo João Seixas Lisbon Portugal Geneva Switzerland Dimensional analysis reveals general kinematic scaling rules for the momentum, mass, and energy dependence of Drell-Yan and quarkonium cross sections. Their application to mid-rapidity LHC data provides strong experimental evidence supporting the validity of the factorization ansatz, a cornerstone of non-relativistic QCD, still lacking theoretical demonstration. Moreover, data-driven patterns emerge for the factorizable long-distance bound-state formation effects, including a remarkable correlation between the S-wave quarkonium cross sections and their binding energies. Assuming that this scaling can be extended to the P-wave case, we obtain precise predictions for the not yet measured feed-down fractions, thereby providing a complete picture of the charmonium and bottomonium feed-down structure. This is crucial information for quantitative interpretations of quarkonium production data, including studies of the suppression patterns measured in nucleus-nucleus collisions. 1 Introduction Quarkonium production is a paradigm case study for the understanding of hadron formation, but the theoretical description of its basic mechanisms remains a challenge [ 1 ]. Both the attractiveness and the complexity of the problem reside in the quark-antiquark (Q Q) binding process, whose degrees of freedom manifest themselves in the observed spectrum of bound states of different masses and spin– angular-momentum properties, remaining unclear how they are interrelated with the momentum distributions measured for the different states. Non-relativistic QCD (NRQCD) [ 2 ], as well as other approaches to quarkonium production (e.g., colour-singlet and colour-evaporation models [ 3 ]) are based on a pivotal postulate: the long-distance effects related to bound-state formation can be factorized out in the calculation of the cross sections. In this way, short-distance Q Q cross sections (SDCs), containing the dependence on transverse momentum ( pT), rapidity (y), and pp centre-ofmass collision energy (√s), can be calculated perturbatively [ 4–8 ], independently of the long-distance matrix elements (LDMEs), which are proportional to the probability of the transition from the 2s+1lcj pre-resonance Q Q, with c = 1 (colour singlet) or 8 (colour octet), to the final-state observable 2S+1 L J quarkonium. Here, s, l, and j (S, L, and J ) denote the spin, orbital angular momentum, and total angular momentum of the primordial Q Q state (quarkonium), respectively. The LDMEs (in practice, free theory parameters) are considered constant: they depend on the initial Q Q quantum numbers and on the bound-state properties (quarkonium mass M , quark mass m Q , velocity of the heavy quark in the quarkonium rest frame v, etc.), but not on the “laboratory” variables pT, y, and √s, nor on the collision system. The factorization postulate, believed to be valid in the limit of heavy quark mass, for which the time of creation of the Q Q pair, Δt ∼ 1/m Q , is distinctively smaller than the bound-state formation time, ∼ 1 fm, has not yet been demonstrated from QCD principles [ 1 ]. Also because of inconsistencies between the NRQCD predictions and the measurements, especially concerning polarization [ 9,10 ], it is often argued that the validity of factorization might be restricted to specific kinematic domains (high pT) or to the (heavier) bottomonium family. The analysis reported in this paper addresses the problem in a data-driven way, comparing the patterns of momentum, mass, and collision-energy dependences measured for the “reference” case of inclusive Drell–Yan production with the corresponding quarkonium patterns. It uses mid-rapidity pp LHC data and only relies on general kinematic scaling properties of the production cross sections. 2 Dimensional scaling in kinematics of particle production The pT and y double-differential cross section for the inclusive production of a given particle in pp collisions can be written as d2σ d pT dy = d2σˆ d pT dyˆ P(x1) P(x2) dx1 dx2, (1) sˆ/M = Eˆ /M + pˆ/M, where σˆ is the cross section of the parton-initiated process, P(x ) are the parton density functions of the proton, x1 and x2 are the proton momentum fractions carried by the colliding partons, and yˆ is the particle rapidity in the parton-parton centre-of-mass frame, yˆ = y − 21 ln(x1/x2). An implicit delta function sets the particle’s pT to the measured one. The cross section depends on the pp collision energy, √s, through the relation x1x2 = (√sˆ/√s)2, where √sˆ is the parton-parton collision energy, function of pT, yˆ and M : √ Eˆ /M = pˆ/M = 1 + ( pT/M )2 cosh yˆ, sinh2 yˆ + ( pT/M )2 cosh2 yˆ. We start by discussing inclusive Drell–Yan lepton-pair production. After integrating over the lepton emission angles in the dilepton rest frame (and over the degrees of freedom of the recoil system), the kinematics in the rest frame of the colliding partons is fully described by pT, yˆ, and M . Since the dimensionality of the differential partonic cross section is [energy]−3, we can write d2σˆ 1 d pT dyˆ = M 3 f (ξ, yˆ), f being a dimensionless function of the dimensionless variables yˆ and ξ , where we use the definition ξ ≡ pT/M to have more compact equations. When the mixture of production mechanisms does not change with M , the partonic cross section, calculated at any arbitrary reference kinematic point (ξ , yˆ ), scales as M −3. This property can be tested using Drell–Yan cross sections, measured in a range of masses sufficiently smaller than the Z boson mass to stay unaffected by the strongly mass-dependent interference of γ - and Z exchange mechanisms. The M scaling of the observable cross section depends on the pp collision energy via the parton densities, an effect that can be singled out in a purely data-driven way at LHC energies and mid rapidity, given the small average x1,2. For illustration, we parametrize P(x ) as (2) (3) (4) P(x ) = A · x −B · q(x ), explicitly showing the factor describing the low-x behaviour and globally representing the remaining factors by the function q, typical examples of which being q(x ) = (1−x )C (1+ D√x + E x ) or q(x ) = (1 − x )C (1 + Dx E ), all tending to unity in the small-x limit. Therefore, P(x1) P(x2) = A2 = A2 1 (x1 x2)B √s 2B M q(x1) q(x2) √s ˆ M −2B q(x1) q(x2). (5) Moreover, dx1dx2 = 2 √s −2 M √s ˆ M 2 Eˆ /M ξ pˆ/M 1 + ξ 2 dξ dyˆ0 , (6) where yˆ0 = y − yˆ = 21 ln(x1/x2) is the rapidity of the system of colliding partons. Using the previous relations and releasing the ξ integration in Eq. 1, we obtain the pp (“dressed”) version of the partonic cross section, Eq. 3, d2σ 1 d pT dy = M 3 √s b M √s F ξ, y; M where b = 2B − 2 and √s F ξ, y; M = g ξ, y − yˆ0 q(x1) q(x2) dyˆ0. The function g absorbs all the factors depending on ξ and yˆ, including those contained in √sˆ/M , Eˆ /M , and pˆ/M (Eq. 2). F is a fully experimentally determinable dimensionless function of the dimensionless scaling variables ξ and y. It also depends on √s and M , but always through the √ √s/M ratio, via the terms q(x1,2) = q( sˆ/√s e±yˆ0 ) and the integration extremes ± ln(√s/ sˆmin ) = ± ln(√s/M )× 1 + ξ 2) / ln(√s/M )). The sensitivity of F on (√1s−/ Mln(dξec+reases with increasing √s and towards mid rapidity. We normalize F over the (ξ, y) range of the analysis and assimilate its integral I (√s/M ) in a numerical redefinition of the b exponent (always possible to a very good approximation over a small range of √s values, as in our analysis of √s = 7 and 8 TeV data), d2σ 1 d pT dy = M 3 √s β M √s F ξ, y; M with F dξ dy = 1. The ( pT/M, y)-integrated cross section becomes σ = M −2 (√s/M )β and its measurement at two or more √s values provides a determination of β. These considerations lead to the following proposition. In a domain where Drell–Yan production mechanisms do not interfere in varying proportions, the joint distribution of the (7) (8) (9) Drell−Yan CMS 8 TeV CMS 7 TeV x 0.8 ATLAS 7 TeV (2011) x 0.6 ATLAS 7 TeV (2010) x 0.4 20 30 40 50 M (GeV) variables pT/M , y, and √s/M is universal: its shape does not depend on M . At any given kinematic point (( pT/M ) , y ) and √s, the dσ/d pT cross section scales like M −(3+β), where β expresses the increase of the integrated cross section with √s. The scaling of the cross section with M , integrated over pT, was previously discussed in Ref. [ 11 ], but using process-specific arguments rather than completely general dimensional-analysis considerations. We will now consider the Drell–Yan dσ/d M differential cross sections measured by ATLAS [ 12 ] and CMS [ 13,14 ], in pp collisions at √s = 7 and 8 TeV, shown in Fig. 1 for the M < 50 GeV range. It can be recognized from the previous discussion that dσ/d M has the same M −3 (√s/M )β scaling behaviour as dσ/d pT at fixed pT/M and y. Simultaneously fitting all four data sets to a single M −αDY function, only considering point-to-point uncorrelated uncertainties in order to determine the shape parameter with maximum significance, gives a remarkably precise result: αDY = 3.63±0.03. Fitting each data set individually gives: 3.60 ± 0.05 (CMS 7 TeV), 3.63 ± 0.05 (CMS 8 TeV), 3.75 ± 0.07 (ATLAS 7 TeV, 2010) and 3.60 ± 0.07 (ATLAS 7 TeV, 2011). The β exponent, reflecting the √s dependence, can be derived from the ratio between the cross sections measured at 8 and 7 TeV, directly reported by CMS [ 14 ], leading to β = 0.73 ± 0.15, where the uncertainty is obtained assuming a relative uncorrelated luminosity uncertainty of 2%. The resulting α −β difference, 2.90±0.15, is perfectly consistent with 3, as expected. The Drell–Yan cross sections so far reported by the LHC experiments are integrated over pT. The future availability of pT- and pT/M -differential measurements will allow the realisation of more detailed and accurate tests. It is worth noting that the existence of these simple and general kinematic scaling rules has been ignored in all the published analyses of LHC data, an unfortunate situation because of the resulting loss in physics value and also because they can be very useful experimental tools. In fact, detector and trigger acceptances drastically sculpt the reconstructed dilepton distributions, especially in the low momentum and low mass regions. Verifying that the measured pT/M and M spectra, after efficiency corrections, etc., satisfy the expected dimensional scaling constitutes a powerful cross-check of the analysis procedure and a validation of the reported systematic uncertainties. 3 Dimensional scaling in quarkonium production Moving now to our main case study, quarkonium production, we write the differential cross section for the inclusive production of a given S-wave state as d2σ d pT dy = m−Q3 √s β M i Li (m Q , M, ξ, y, √s/M ) m3Q (10) × Fi (m Q , M, ξ, y, √s/M )), where the overall factor m−Q3 matches the global dimensionality of the observable. The only specification of the functions Li is that they have dimension [energy]3, formally compensated by the m3Q denominators, while the Fi are dimensionless shape functions (defined with F dξ dy = 1, as in the Drell–Yan case). The √s/M power-law factor represents the modification from partonic to observable level, as from Eqs. 3 to 9. The coefficient β is the same as for Drell–Yan production. From the previous discussion, we can obtain a precise evaluation of its value at √s 7 or 8 TeV as the difference between the average experimental mass scaling exponent αDY and the expected bare-cross-section scaling exponent: β = αDY − 3 = 0.63 ± 0.03. The expression above, Eq. 10, is built in analogy to the NRQCD factorized expansion. In this limit, Fi and Li encode the information on, respectively, the hard scattering process producing a Q Q pair, of mass 2 m Q , and its ensuing long-distance transition to the final bound state Q, of mass M . In NRQCD the LDMEs Li are independent of laboratory kinematics and only depend on degrees of freedom naturally defined in the Q rest frame (M , m Q , v, quantum numbers, spectroscopic energy levels, etc.), while the kinematic dependence is contained in the SDCs, corresponding to Fi /m6Q (√s/M )β , where the index i indicates one specific Q Q(2s+1lc) → Q(2S+1 L J ) production channel. j In this formulation, a consequence of the factorization hypothesis is that each Fi function should be independent of M , being, in particular, the same for all directly produced 3 S1 charmonia and bottomonia. In fact, it is important to note that the mass difference between the Q and Q Q states plays no role in the pT/M, yˆ dependence of Fi , given that, similarly to the physical transitions between states of the same quarkonium family [ 15 ], the Q Q → Q transition preserves the ratio between laboratory-vector-momentum and mass (and, therefore, both pT/M and yˆ). While these notations have been chosen to accommodate the NRQCD factorization expansion as a limit case, Eq. 10, with generic Li and Fi both redundantly depending of the relevant variables, represents a fully general template of the cross section, with no prejudice on its physical scaling properties. Very interesting physical indications transpire from the unforeseen simplicity of the trends measured by ATLAS and CMS, as discussed in the next paragraphs. The first experimental input to our considerations are the quarkonium cross sections measured at mid-rapidity (|y| 2) by the ATLAS and CMS experiments, shown in the top panels of Fig. 2 for √s = 7 TeV (left) and 13 TeV (right). All the measurements, from the J/ψ to the Υ (3S), scale with pT/M in a state-independent way, at least for not very low pT ( pT/M 2). This “universality” is well illustrated by the middle panels, which clearly show that, after rescaling the state-specific normalizations to those of the J/ψ , the cross section shapes become indistinguishable from each other. The goodness of this universal scaling can be better appreciated by looking at the bottom panels, which show (in a linear scale) the pull distributions, i.e. the differences between each data point and the universal fitted function, divided by the measurement uncertainty. No systematic trends are seen in the pull distributions and the observed deviations are very well compatible with statistical fluctuations. The cross section fits consider correlations between the luminosity uncertainties in each data set and the pulls, evaluated to check the consistency with one common shape, are calculated excluding such uncertainties. A slightly broader pT/M distribution is observed at the higher energy, as shown in Fig. 3-top. The fact that the distinct P- and S-wave states show a compatible kinematic scaling is discussed in Ref. [ 22 ]. Here we focus on the precisely measured cross sections of the closely related 3 S1 states. Their universal pT/M scaling, at a given energy, indicates that the same mixture of processes (or one common dominating process) describes the production of all these states. The second experimental indication for our following discussion is the mass scaling of the cross section from charmonium to bottomonium. By exploiting the pT/M scaling we can determine it without relying on model-dependent extrapolations to low pT: it is sufficient to consider the fitted pT/M distributions at an arbitrary pT/M = ( pT/M ) . We consider for each family the meson closest to the ground state, J/ψ and Υ (1S), at the two energies (Fig. 3-bottom). To obtain the yield of directly produced J/ψ mesons, we subtract from the fitted prompt-J/ψ normalization those of the χc1, χc2 and ψ (2S) states (always considered at the same pT/M value: as mentioned before, pT/M is transferred unchanged from mother to daughter), weighted by the branching fractions [ 23 ] of the respective feed-down processes. The result is a J/ψ direct-production fraction of 31.9%. For the Υ (1S) we assume a direct-production fraction of (50 ± 10)%. Later in this paper we will need the corresponding fractions for the Υ (2S) and Υ (3S), which we assume to be 60 and 70%, respectively, also with ±10% uncertainties. These values were estimated from Υ (nS) yield ratios measured by ATLAS [ 17 ] and CMS [ 20 ], complemented by LHCb data on the (large) χb(mP) → Υ (nS) feed-down fractions [ 24 ]. The resulting mass scaling, measured as (11) dσ/d pT(Υ (1S)) ξ=ξ dσ/d pT(J/ψ ) ξ=ξ = mb mc −α does not show a significant energy dependence: α = 6.6±0.1 and 6.5 ± 0.1 at 7 and 13 TeV, respectively, with 2 mc M (ηc(1S)) = 2.984 GeV and 2 mb M (ηb(1S)) = 9.389 GeV [ 23 ]. Subtracting the β value mentioned above, 0.63 ± 0.03, measured in Drell–Yan production at 7 and 8 TeV, we find that the partonic-level differential cross section changes as m−Q6 between the charmonium and bottomonium families. These experimental facts constrain and specify the functions L and F appearing in Eq. 10. The independence of the pT/M scaling on either M or m Q , at a given energy, indicates that the pT/M and (M, m Q ) dependences do not “mix”: L × F = L(m Q , M, √s/M ) × F (ξ, y, √s/M ) (further studies considering forward-rapidity, low- pT data will assess the combined pT/M, y dependence of F , here effectively a function of only pT/M ). In other words, there is experimental evidence that, at mid rapidity and not too low pT, it is possible to describe quarkonium production “factorizing” short- and long-distance effects, defined as the dependences on, respectively, the laboratory momentum of the detected meson (F ) and the specific bound-state properties (L). The measured m−Q6 scaling of the differential cross section at a fixed pT/M (and y) further specifies L. By precisely equating the explicit m Q dependence of the expression in Eq. 10 (coming from the overall factor and the denominators of the L terms, with F now taken independent of m Q ), such result leaves no margin for a dependence of L on m Q , if not counterbalanced by a dependence on M and/or √s: L = L(m Q /M, m Q /√s). However, given that no significant difference in mass scaling is observed at 7 and 13 TeV, a √ m Q / s dependence would be experimentally indistinguishable from the analogous global energy scaling of all quarkonium cross sections already accounted for by the β power law in Eq. 10. It is therefore always possible to define the longdistance factors as √s-independent, thereby reducing them resent a single universal empirical function, of shape determined by a simultaneous fit to all data points of pT/M > 2 and normalizations specific to each quarkonium state. The bottom panels show the pulls between each data point and the fitted function 7 TeV J/ψ χc1 χc2 ψ(2S) Υ(1S) Υ(2S) Υ(3S) 13 TeV J/ψ ψ(2S) Υ(1S) Υ(2S) Υ(3S) Υ(1S) 13 TeV 7 TeV CMS ATLAS CMS 7 13 J/ψ 13 TeV 7 TeV CMS ATLAS CMS 1 2 10 20 pT/M to functions of the kind L = L(m Q /M ), as is illustrated by the discussion starting in the next paragraph. It should be clear from these considerations that such L functions do not coincide with the LDMEs of NRQCD: instead of being defined by setting an energy scale within the theory, they are built on the basis of dimensional-analysis and data patterns, becoming universal, experimentally definable quantities. Apart from differences in operative definitions, we note that the remarkably simple picture of S-wave quarkonium production emerging from the data mirrors well the primary concepts of NRQCD factorization and of uni103 102 ) eV10 G / b n ( pT 1 d / σ d 10−1 3 Υ(1S) J/ψ in nb/GeV in pb/GeV versality of the long-distance bound-state formation effects. A future verification of the corresponding experimental patterns in different collision systems can further probe these fundamental concepts. By adopting now the guiding idea of factorization of shortand long-distance effects, as supported by data, we can analyse in more detail the mass scaling, using all S-wave states. Figure 4-top shows the values of the direct production cross sections of the S-wave quarkonia, for one specific and arbitrarily chosen point, pT/M = 7, as derived from the F ( pT/M ) functions shown in the top panel of Fig. 2. Given the universality of the function, a different pT/M reference point only leads to an overall vertical rescaling of the points in Fig. 4, not affecting the mass dependence itself. For the J/ψ and Υ (nS), the prompt cross sections were scaled down by the direct-production fractions mentioned earlier. The figure only shows the 7 TeV pattern; completely analogous results are found in the 13 TeV data (apart from a significantly different global normalization scale), as discussed hereafter. 7= 1 M /pT r o f ) Q m 2 ( σ /) ϒ | ψ ( σ10-1 7 TeV From the J/ψ to the Υ (3S) we observe, besides the global drop by three orders of magnitude, a “fine structure”: within each quarkonium family, the cross section decreases faster than from one family to the other. In the factorization perspective, the global drop represents the short-distance scaling, reflecting the change from M (Q Q) 2 mc to 2 mb. The steeper change within each family reflects the M -dependent Q Q → Q transition probability L(m Q , M ). The factorized interpretation is represented by the fit curves shown in the figure. First (red and blue curves) the charmonium and bottomonium cross sections are extrapolated, respectively, to 2 mc and 2 mb (as defined above). From the resulting σ (2mb) / σ (2mc) ratio, the short-distance charmonium-tobottomonium mass scaling (black curve) is determined as ∝ m−QαQQ , with αQ Q = 6.63 ± 0.08, practically identical to the previously obtained result, considering only the J/ψ and the Υ (1S), leading to the m−Q6 partonic-level scaling discussed above. Here we focus on the long-distance mass dependence. The bottom panel of Fig. 4 shows the dependence of the σ (Q)/σ (2 m Q ) ratio on M/(2 m Q ): one common power-law, σ (Q)/σ (2 m Q ) = [M/(2 m Q )]−(9.7±0.3), describes very well both the ψ and the Υ data points. Interestingly, as shown in Fig. 5, this monotonic dependence can be seen as a tight correlation of the longdistance factors to the quarkonium binding energy, defined as 2M (D0)− M (ψ (nS)) or 2M (B0)− M (Υ (nS)) and calculated with mass values from Ref. [ 23 ]. As a function of this variable we see another “universal” trend: the data points are well described by the parametrization σ (Q)/σ (2 m Q ) ∝ Ebδinding, common to charmonium and bottomonium and identical at 7 TeV (δ = 0.63 ± 0.02) and 13 TeV (δ = 0.63 ± 0.04). The central values of the δ and β parameters are identical by coincidence. This result gives further support to the idea that the dependence of the cross sections on the bound-state masses is a “factorizable” long-distance effect, independent of the laboratory-momentum. We have until now considered only S-wave states. Given the absence of detailed χb cross section data at mid rapidity, we can put to test the assumption that the “universal” σ (Q)/σ (2 m Q ) ∝ Ebδinding dependence of the longdistance factors, with δ 0.63, can be extended to the P-wave states, even if with a different normalization constant multiplying Ebδinding to account for a dependence on the orbital angular momentum. In this way, χc data come to constrain the χb(nP) cross sections, through the relation σ (χb J (nP)) = [Ebinding(χb J (nP))/Ebinding(χc J )]δσ (χc J ), corresponding to the assumption that there is no dependence on J , n, or quark flavour. Using branching-ratio measurements from Ref. [ 23 ], the feed-down structure of quarkonium production can be fully predicted. The result is presented in Table 1. It should be kept in mind that these results correspond to the pT/M > 2 region. The predicted feed-down fractions of Υ (nS) production from χb(nP) states are in reasonable agreement with forward-rapidity LHCb measurements merging the J = 1 and J = 2 signals [ 24 ], considered for pT/M > 2. 4 Summary and discussion At the current level of experimental precision, mid-rapidity LHC proton-proton data for inclusive charmonium and bottomonium production are well described by a simple parametrization reflecting a universal (i.e. state-independent) scaling with two variables: the shapes of the pT distributions of all states become one common shape as a function of pT/M and the normalization of this shape scales in a simple monotonic way with Ebinding, at least for the S-wave states. While the cross section shape as a function of pT/M depends on √s, the normalization scaling with Ebinding is found to be identical at 7 and 13 TeV. Having a simple empirical parametrization faithfully describing the proton-proton measurements is certainly very useful for model-independent studies of quarkonium production, especially when considering more complex collision systems such as proton-nucleus and nucleus-nucleus collisions. It is also very convenient for the tuning of Monte Carlo simulations. More importantly, the analysis reported in this paper, exclusively based on non-trivial (albeit potentially misleadingly simple) dimensional analysis arguments, reveals significant experimental evidence supporting that quarkonium production can be understood as being factorized between short- and long-distance phases. This data-driven result mirrors very well the general concept of NRQCD factorization into process-dependent SDCs and universal LDMEs. More generally, seeing evidence of factorization in the LHC data is an important step towards establishing that the Q Q pairs are produced uncorrelated from the rest of the event, allowing us to experimentally probe how they evolve and bind into the observable spectra of quarkonium states. The experimental verification of the factorization ansatz can be extended by including in the analysis the available forward-rapidity LHCb measurements. Such study, currently ongoing, has a higher order of complexity: since the pT/M distributions show a significant y dependence, the determination of scaling patterns is no longer an effectively onedimensional problem as in the case where only mid-rapidity data are considered. Further tests, which may become possible in the future, include the study of the production of quarkonia associated to specific final state particles and the search for simplicity patterns in different kinds of collisions. More data on the production of P-wave quarkonia are needed, especially in the bottomonium family: the χb absolute cross sections have not yet been measured at mid rapidity, and the forward rapidity results, reported relative to Υ (nS) production, do not distinguish between the J = 0, 1, and 2 states. It is reassuring, however, to see that those forward rapidity measurements are in reasonable agreement with the pattern of χb J (nP) → Υ feed-down fractions that we obtained assuming that the simple Ebinding scaling pattern can be extended to the χ states, an assumption that follows very naturally from the work presented in this paper. The detailed feed-down structure of the charmonium and bottomonium families we have determined in this work needs to be tested by new measurements, filling gaps in the experimental picture. It cannot be excluded that future measurements of direct χc J and χb J production will show that the scaling with binding energy is not the same for the P- and S-wave states. Such a result would still be very interesting, naturally. For example, the difference in scaling trends, when known, could be used to see if new resonances, as the X(3872), align better with the S- or the P-wave mass trends, thereby contributing to the understanding of their nature. A complete account of all the feed-down fractions is particularly important in quantitative analyses of the quarkonium production and suppression patterns observed in nucleus-nucleus collisions. It has been argued since long, in particular, that the observed suppression of the J/ψ signal, both from pp to heavy-ion collisions and from peripheral to central nuclear collisions, might be dominantly (or fully) caused by the suppression of the ψ (2S) and χc excited states, more loosely bound and, hence, easier to dissolve in the hot medium created in those interactions. Such hypotheses can only be reliably and quantitatively tested if the hierarchy of feed-down relations between states is taken into account. This computation is now feasible, assuming the global Ebinding scaling that leads to the values collected in Table 1 for pp collisions, complemented by a model for its modification in nucleus-nucleus collisions. Such studies will finally address, and possibly provide evidence for, the concept of sequential quarkonium suppression, an ideal signature of quark-gluon plasma formation [ 25 ] that applies to the directly-produced states, after the subtraction of the feeddown contributions. Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecomm ons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. Funded by SCOAP3. 1. N. Brambilla et al., Eur. Phys. J. C 71 , 1534 ( 2011 ) 2. G. Bodwin , E. Braaten, P. Lepage, Phys. Rev. D 51 , 1125 ( 1995 ) [ Erratum: Phys. Rev. D 55 ( 1997 ) 5853] 3. G. Bodwin , E. Braaten , J. Lee , Phys. Rev. D 72 , 014004 ( 2005 ) 4. M. Butenschön , B.A. Kniehl , Phys. Rev. Lett . 108 , 172002 ( 2012 ) 5. M. Butenschön , B.A. Kniehl , Mod. Phys. Lett. A 28 , 1350027 ( 2013 ) 6. K.-T. Chao , Y.-Q. Ma , H.S. Shao , K. Wang , Y.-J. Zhang , Phys. Rev. Lett . 108 , 242004 ( 2012 ) 7. B. Gong , L.-P. Wan , J.-X. Wang , H.-F. Zhang , Phys. Rev. Lett . 110 , 042002 ( 2013 ) 8. H. -S. Shao , Y.-Q. Ma , K. Wang , K.- T. Chao, Phys. Rev. Lett . 112 , 182003 ( 2014 ) 9. P. Faccioli , C. Lourenço , J. Seixas , H. Wöhri , Eur. Phys. J. C 69 , 657 ( 2010 ) 10. P. Faccioli , V. Knünz , C. Lourenço , J. Seixas , H. Wöhri, Phys. Lett. B 736 , 98 ( 2014 ) 11. N.S. Craigie , Phys. Rep . 47 , 1 ( 1978 ) 12. ATLAS Coll., JHEP 06 , 112 ( 2014 ) 13. CMS Coll., JHEP 12 , 030 ( 2013 ) 14. CMS Coll., Eur. Phys. J. C75 , 147 ( 2015 ) 15. P. Faccioli , M. Araújo , V. Knünz , I. Krätschmer, C. Lourenço , J. Seixas , Phys. Lett. B 773 , 476 ( 2017 ) 16. ATLAS Coll., JHEP 09 , 079 ( 2014 ) 17. ATLAS Coll., Phys. Rev. D 87 , 052004 ( 2013 ) 18. ATLAS Coll., JHEP 07 , 154 ( 2014 ) 19. CMS Coll., Phys. Rev. Lett . 114 , 191802 ( 2015 ) 20. CMS Coll., Phys. Lett. B 749 , 14 ( 2015 ) 21. CMS Coll., Phys. Lett. B. arXiv:1710 .11002 (submitted) 22. P. Faccioli et al., Eur. Phys. J . C arXiv: 1802 . 01106 (submitted) 23. Particle Data Group, Chin. Phys. C 40 , 100001 ( 2016 ) 24. LHCb Coll., Eur. Phys. J. C 74 , 3092 ( 2014 ) 25. S. Digal, P. Petreczky , H. Satz, Phys. Rev. D 64 , 094015 ( 2001 )

This is a preview of a remote PDF: https://link.springer.com/content/pdf/10.1140%2Fepjc%2Fs10052-018-5610-x.pdf

Pietro Faccioli, Carlos Lourenço, Mariana Araújo, João Seixas. Universal kinematic scaling as a probe of factorized long-distance effects in high-energy quarkonium production, The European Physical Journal C, 2018, 118, DOI: 10.1140/epjc/s10052-018-5610-x