Measurement of the relative prompt production rate of chi(c2) and chi(c1) in pp collisions at sqrt(s) = 7 TeV

A measurement is presented of the relative prompt production rate of chi(c2) and chi(c1) with 4.6 inverse femtobarns of data collected by the CMS experiment at the LHC in pp collisions at sqrt(s) = 7 TeV. The two states are measured via their radiative decays chi(c) to J/psi + gamma, with the photon converting into a dielectron pair for J/psi rapidity abs(y(J/psi))<1.0 and photon transverse momentum pt(gamma)>0.5 GeV. The measurement is given for six intervals of pt(J/psi) between 7 and 25 GeV. The results are compared to theoretical predictions.


Introduction
Understanding charmonium production in hadronic collisions is a challenge for quantum chromodynamics (QCD).The J/ψ production cross section measurements at the Tevatron [1,2] were found to disagree by about a factor of 50 with theoretical color-singlet calculations [3].Soon after, the CDF experiment reported a χ c2 /χ c1 cross section ratio that extended up to p T (J/ψ) 10 GeV/c, where p T is the transverse momentum, and favored χ c1 production over χ c2 [4].The cross section ratio was also studied recently at the Large Hadron Colllider (LHC) in Ref. [5].These measurements independently suggest that charmonium production cannot be explained through relatively simple models.This paper presents a measurement of the prompt relaxχ c2 /χ c1 cross section ratio by the Compact Muon Solenoid (CMS) experiment at the LHC in pp collisions at a center-of-mass energy of 7 TeV.Prompt refers to the production of χ c mesons that originate from the primary pp interaction point, as opposed to the ones from the decay of B hadrons.Prompt production includes both directly produced χ c and also indirectly produced χ c from the decays of shortlived intermediate states, e.g. the radiative decay of the ψ(2S).The measurement is based on the reconstruction of the χ c radiative decays to J/ψ + γ, with the low transverse momentum photons (less than 5 GeV/c) being detected through their conversion into electron-positron pairs.The analysis uses data collected in 2011, corresponding to a total integrated luminosity of 4.6 fb −1 .When estimating acceptance and efficiencies, we assume that the χ c2 and χ c1 are produced unpolarized, and we supply the correction factors needed to modify the results for several different polarization scenarios.
Due to the extended reach in transverse momentum made possible by the LHC energies, the cross section ratio measurement is expected to discriminate between different predictions, such as those provided by the relaxk T -factorization [6] and next-to-leading order nonrelativistic QCD (NRQCD) [7] theoretical approaches.
The strength of the ratio measurement is that most theoretical uncertainties cancel, including the quark masses, the value of the strong coupling constant α s , as well as experimental uncertainties on quantities such as integrated luminosity, trigger efficiencies, and, in part, reconstruction efficiency.Therefore, this ratio can be regarded as an important reference measurement to test the validity of various theoretical quarkonium production models.With this paper, we hope to provide further guidance for future calculations.

CMS detector
A detailed description of the CMS apparatus is given in Ref. [8].Here we provide a short summary of the detectors relevant for this measurement.
The central feature of the CMS apparatus is a superconducting solenoid of 6 m internal diameter.Within the field volume are the silicon pixel and strip tracker, the crystal electromagnetic calorimeter and the brass/scintillator hadron calorimeter.Muons are measured in gasionization detectors embedded in the steel return yoke.In addition to the barrel and endcap detectors, CMS has extensive forward calorimetry.
The inner tracker measures charged particles within the pseudorapidity range |η| < 2.5, where η = − ln[tan(θ/2)], and θ is the polar angle measured from the beam axis.It consists of 1440 silicon pixel and 15 148 silicon strip detector modules.In the central region, modules are arranged in 13 measurement layers.It provides an impact parameter resolution of ∼15 µm.
Muons are measured in the pseudorapidity range |η| < 2.4, with detection planes made using three technologies: drift tubes, cathode strip chambers, and resistive plate chambers.Matching the muons to the tracks measured in the silicon tracker results in a transverse momentum resolution between 1 and 1.5%, for p T values up to 50 GeV/c.The first level (L1) of the CMS trigger system, composed of custom hardware processors, uses information from the calorimeters and muon detectors to select the most interesting events.The high-level trigger (HLT) processor farm further decreases the event rate from around 100 kHz to around 300 Hz, before data storage.The rate of HLT triggers relevant for this analysis was in the range 5-10 Hz.We analyzed about 60 million such triggers.

Experimental method
We select χ c1 and χ c2 candidates by searching for their radiative decays into the J/ψ + γ final state, with the J/ψ decaying into two muons.The χ c0 has too small a branching fraction into this final state to perform a useful measurement, but we consider it in the modeling of the signal lineshape.Given the small difference between the J/ψ mass, 3096.916± 0.011 MeV/c 2 , and the χ c1 and χ c2 masses, 3510.66 ± 0.07 MeV/c 2 and 3556.20 ± 0.09 MeV/c 2 , respectively [9], the detector must be able to reconstruct photons of low transverse momentum.In addition, excellent photon momentum resolution is needed to resolve the two states.In the center-of-mass frame of the χ c states, the photon has an energy of 390 MeV when emitted by a χ c1 and 430 MeV when emitted by a χ c2 .This results in most of the photons having a p T in the laboratory frame smaller than 6 GeV/c.The precision of the cross section ratio measurement depends crucially on the experimental photon energy resolution, which must be good enough to separate the two states.A very accurate measurement of the photon energy is obtained by measuring electronpositron pairs originating from a photon conversion in the beampipe or the inner layers of the silicon tracker.The superior resolution of this approach, compared to a calorimetric energy measurement, comes at the cost of a reduced yield due to the small probability for a conversion to occur in the innermost part of the tracker detector and, more importantly, by the small reconstruction efficiency for low transverse momentum tracks whose origin is displaced with respect to the beam axis.Nevertheless, because of the high χ c production cross section at the LHC, the use of conversions leads to the most precise result.
For each χ c1,2 candidate, we evaluate the mass difference ∆m = m µµγ − m µµ between the dimuon-plus-photon invariant mass, m µµγ , and the dimuon invariant mass, m µµ .We use the quantity Q = ∆m + m J/ψ , where m J/ψ is the world-average mass of the J/ψ from Ref. [9], as a convenient variable for plotting the invariant-mass distribution.We perform an unbinned maximum-likelihood fit to the Q spectrum to extract the yield of prompt χ c1 and χ c2 as a function of the transverse momentum of the J/ψ.A correction is applied for the differing acceptances for the two states.Our results are given in terms of the prompt production ratio R p , defined as where σ(pp → χ c + X) are the χ c production cross sections, B(χ c → J/ψ + γ) are the χ c branching fractions, N χ i are the number of candidates of each type obtained from the fit, and ε 1 /ε 2 is the ratio of the efficiencies for the two χ c states.The branching fractions B(χ c1,2 → J/ψ + γ), taken from Ref. [9], are also used to calculate the ratio of production cross sections.

Event reconstruction and selection
In order to select χ c signal events, a dimuon trigger is used to record events containing the decay J/ψ → µµ.The L1 selection requires two muons without an explicit constraint on their transverse momentum.At the HLT, opposite-charge dimuons are reconstructed and the dimuon rapidity y(µµ) is required to satisfy |y(µµ)| < 1.0, while the dimuon p T must exceed a threshold that increased from 6.5 to 10 GeV/c as the trigger configuration evolved to cope with the instantaneous luminosity increase.Events containing dimuon candidates with invariant mass from 2.95 to 3.25 GeV/c 2 are recorded.Our data sample consists of events where multiple pp interactions occur.At each bunch crossing, an average of six primary vertices is reconstructed, one of them related to the interaction that produces the χ c in the final state, the others related to softer collisions (pileup).
In the J/ψ selection, the muon tracks are required to pass the following criteria.They must have at least 11 hits in the tracker, with at least two in the pixel layers, to remove background from decays-in-flight.The χ 2 per degree of freedom of the track fit must be less than 1.8.To remove background from cosmic-ray muons, the tracks must intersect a cylindrical volume of radius 4 cm and total length 70 cm, centered at the nominal interaction point and with its axis parallel to the beam line.Muon candidate tracks are required to have p T > 3.3 GeV/c, |η| ≤ 1.3 and match a well-reconstructed segment in at least one muon detector [10].Muons with opposite charges are paired.The two muon trajectories are fitted with a common vertex constraint, and events are retained if the fit χ 2 probability is larger than 1%.If more than one muon pair is found in an event, only the pair with the largest vertex χ 2 probability is selected.
For the final χ c1 and χ c2 selection, a dimuon candidate must have an invariant mass between 3.0 and 3.2 GeV/c 2 and |y| < 1.0.
In order to restrict the measurement to the prompt J/ψ signal component, the pseudo-proper decay length of the J/ψ ( J/ψ ), defined as J/ψ = L xy • m J/ψ /p T (J/ψ), where L xy is the most probable transverse decay length in the laboratory frame [11], is required to be less than 30 µm.
In the region J/ψ < 30 µm, we estimate, from the observed J/ψ distribution, a contamination of the nonprompt component (originating from the decays of B hadrons) of about 0.7%, which has a negligible impact on the total systematic uncertainty.
To reconstruct the photon from radiative decays, we use the tracker-based conversion reconstruction described in Refs.[12][13][14].We summarize the method here, mentioning the further requirements needed to specialize the conversion reconstruction algorithm to the χ c case.The algorithm relies on the capability of iterative tracking to efficiently reconstruct displaced and low transverse momentum tracks.Photon conversions are characterized by an electron-positron pair originating from a common vertex.The e + e − invariant mass must be consistent with zero within its uncertainties and the two tracks are required to be parallel at the conversion point.
Opposite-sign track pairs are first required to have more than four hits and a normalized χ 2 less than 10.Then the reconstruction algorithm exploits the conversion-pair signature to distinguish between genuine and misidentified background pairs.Information from the calorimeters is not used for conversion reconstruction in our analysis.The primary pp collision vertex associated with the photon conversion, see below, is required to lie outside both track helices.Helices projected onto the transverse plane form circles; we define d m as the distance between the centers of the two circles minus the sum of their radii.The value of d m is negative when the two projected trajectories intersect.We require the condition −0.25 < d m < 1.0 cm to be satisfied.From simulation, we have found that most of the electron-positron candidate pair background comes from misreconstructed track pairs originating from the primary vertex.These typically have negative d m values, thus explaining the asymmetric d m requirements.
In order to reduce the contribution of misidentified conversions from low-momentum displaced tracks that are artificially propagated back to the silicon tracker, the two candidate conversion tracks must have one of their two innermost hits in the same silicon tracker layer.
The distance along the beam line between the extrapolation of each conversion track candidate and the nearest reconstructed event vertex must be less than five times its estimated uncertainty.Moreover, among the two event vertices closest to each track along the beam line, at least one vertex must be in common A reconstructed primary vertex is assigned to the reconstructed conversion by projecting the photon momentum onto the beamline and choosing the closest vertex along the beam direction.If the value of the distance is larger than five times its estimated uncertainty, the photon candidate is rejected.
The primary vertex associated with the conversion is required to be compatible with the reconstructed J/ψ vertex.This requirement is fulfilled when the three-dimensional distance between the two vertices is compatible with zero within five standard deviations.Furthermore, a check is made that neither of the two muon tracks used to define the J/ψ vertex is used as one of the conversion track pair.
The e + e − track pairs surviving the selection are then fitted to a common vertex with a kinematic vertex fitter that constrains the tracks to be parallel at the vertex in both the transverse and longitudinal planes.The pair is retained if the fit χ 2 probability is greater than 0.05%.If a track is shared among two or more reconstructed conversions, only the conversion with the larger vertex χ 2 probability is retained.
Only reconstructed conversions with transverse distance of the vertex from the center of the mean pp collision position larger than 1.5 cm are considered.This requirement suppresses backgrounds caused by track pairs originating from the primary event vertex that might mimic a conversion, such as from π 0 Dalitz decay, while retaining photon conversions occurring within the beampipe.
Finally, each conversion candidate is associated with every other conversion candidate in the event, and with any photon reconstructed using calorimeter information.Any pairs of conversions or conversion plus photon with an invariant mass between 0.11 and 0.15 GeV/c 2 , corresponding to a two-standard-deviation window around the π 0 mass, is rejected.We have verified that the π 0 rejection requirement, while effectively reducing the background, does not affect the R p measurement within its uncertainties.
Converted photon candidates are required to have p T > 0.5 GeV/c, while no requirement is imposed on the pseudorapidity of the photon.
The distribution of the photon conversion radius for χ c candidates is shown in Fig. 1.The first peak corresponds to the beampipe and first pixel barrel layer, the second and third peaks correspond to the two outermost pixel layers, while the remaining features at radii larger than 20 cm are due to the four innermost silicon strip layers.The observed distribution of the photon conversion radius is consistent with the known distribution of material in the tracking volume and with Monte Carlo simulations [14].

Acceptance and efficiencies
In the evaluation of R p , we must take into account the possibility that the geometric acceptance and the photon reconstruction efficiencies are not the same for χ c1 and χ c2 .In order to determine the acceptance correction, a Monte Carlo (MC) simulation sample of equal numbers of χ c1 and χ c2 has been used.This sample was produced using a PYTHIA [15] single-particle simulation in which a χ c1 or χ c2 is generated with a transverse momentum distribution produced from a parameterized fit to the CMS measured ψ(2S) spectrum [16].The use of the ψ(2S) spectrum is motivated by the proximity of the ψ(2S) mass to the states under examination.The impact of this choice is discussed in Section 7.

Conversion radius [cm]
Both χ c states in the simulation are forced to decay to J/ψ + γ isotropically in their rest frame, i.e., assuming they are produced unpolarized.We discuss later the impact of this assumption.The decay products are then processed through the full CMS detector simulation, based on GEANT4 [17,18], and subjected to the trigger emulation and the full event reconstruction.In order to produce the most realistic sample of simulated χ c decays, digitized signals from MCsimulated inelastic pp events are mixed with those from simulated signal tracks.The number of inelastic events to mix with each signal event is sampled from a Poisson distribution to accurately reproduce the amount of pileup in the data.
The efficiency ratio ε 1 /ε 2 for different J/ψ transverse momentum bins is determined using : , where N gen is the number of χ c candidates generated in the MC simulation within the kinematic range |y(J/ψ)| < 1.0, p T (γ) > 0.5 GeV/c, and N rec is the number of candidates reconstructed with the selection above.The resulting values are shown in Table 1, where the uncertainties are statistical only and determined from the MC sample assuming binomial distributions.The increasing trend of ε 1 /ε 2 is expected, because p T (J/ψ) is correlated with the p T of the photon, and at higher photon p T our conversion reconstruction efficiency is approximately constant.Therefore, efficiencies for the χ c1 and the χ c2 are approximately the same at high p T (J/ψ).
This technique also provides an estimate of the absolute χ c reconstruction efficiency, which is given by the product of the photon conversion probability, the χ c selection efficiency, and, most importantly, the conversion reconstruction efficiency, which corresponds to the dominant contribution.This product varies as a function of p T (γ), and goes from 4 × 10 −4 at 0.5 GeV/c to around 10 −2 at 4 GeV/c, where it saturates.

Signal extraction
We extract the numbers of χ c1 and χ c2 events, N χ c1 and N χ c2 , respectively, from the data by performing an unbinned maximum-likelihood fit to the Q spectrum in various ranges of J/ψ transverse momentum.
Because of the small intrinsic width of the χ c states we are investigating, the observed signal shape is dominated by the experimental resolution.The signal probability density function (PDF) is derived from the MC simulation described in Section 5, and is modeled by the superposition of two double-sided Crystal Ball functions [19] for the χ c1 and χ c2 and a single-sided Crystal Ball function for the χ c0 .Each double-sided Crystal Ball function consists of a Gaussian core with exponential tails on both the high-and low-mass sides.We find this shape to provide an accurate parameterization of the Q spectra derived from MC simulation.When fitting the data, we fix all the parameters of the Crystal Ball function to the values that best fit our MC simulation and use a maximum-likelihood approach to derive N χ c1 and N χ c2 , which are the integrals of the PDFs for the two resonances.Because the Q resolution depends on the p T of the J/ψ, a set of shape parameters is determined for each bin of p T (J/ψ).Simulation shows that the most important feature of the χ c0 signal shape is the low-mass tail due to radiation from the electrons, while the high-mass tail is overwhelmed by the combinatorial background and the low-mass tail of the other resonances.Hence the choice to use a single-sided Crystal Ball function to fit the χ c0 mass distribution.Different choices of the χ c0 signal parameterization are found to cause variations in the measured R p values that are well within the quoted systematic uncertainties given below.
The background is modeled by a probability distribution function defined as where α 1 and β 1 are free parameters in the fit, and q 0 is set to 3.2 GeV/c 2 .
In Fig. 2 we show the Q distribution for two different ranges, 11 < p T (J/ψ) < 13 GeV/c (left) and 16 < p T (J/ψ) < 20 GeV/c (right).This procedure is repeated for several ranges in the transverse momentum of the J/ψ in order to extract N χ c1 and N χ c2 in the corresponding bin.
The results are shown in Table 2, where the reported uncertainties are statistical only.

Systematic uncertainties
Several types of systematic uncertainties are addressed.In particular, we investigate possible effects that could influence the measurement of the numbers of χ c1 and χ c2 from data, the Table 2: Numbers of χ c1 and χ c2 events extracted from the maximum-likelihood fit, and the ratio of the two values.Uncertainties are statistical only.evaluation of ε 1 /ε 2 from the MC simulation, and the derivation of the R p ratio.In Table 3 the various sources of systematic uncertainties and their contributions to the total uncertainty are summarized.The following subsections describe how the various contributions are evaluated.

Uncertainty from the mass fit and χ c1 and χ c2 counting
The measurement of the ratio N χ c2 /N χ c1 could be affected by the choice of the functional form used for the maximum-likelihood fit.The use of an alternative background parameterization, a fourth-order polynomial, results in systematically higher values of the ratio N χ c2 /N χ c1 , while keeping the overall fit quality as high as in the default procedure.From the difference in the numbers of signal events using the two background parameterizations, we assign the systematic uncertainty from the background modeling shown in Table 3.
We evaluate the systematic uncertainty related to the parameterization of the signal shape by varying the parameters derived from the MC simulation within their uncertainties.The results fluctuate within 1-3% in the various transverse momentum ranges.We assign the systematic uncertainties from this source, as shown in Table 3.
The method to disentangle and count the χ c1 and χ c2 states is validated by using a PYTHIA MC simulation sample of inclusive J/ψ events, including those from χ c decay, produced in pp collisions and propagated through the full simulation of the detector.The ratio N χ c2 /N χ c1 derived from the fit to the Q distribution of the reconstructed candidates in the simulation is consistent with the actual number of χ c events contributing to the distribution, within the statistical uncertainty, for all J/ψ momentum ranges.Therefore, we do not assign any further systematic uncertainty in the determination of N χ c2 /N χ c1 .
The stability of our analysis as a function of the number of primary vertices in the event has been investigated.The number of χ c candidates per unit of integrated luminosity, once trigger conditions are taken into account, is found to be independent of the instantaneous luminosity, within the statistical uncertainties.In addition, the measured ratio N χ c2 /N χ c1 is found to be constant as a function of the number of primary vertices in the event, within the statistical uncertainties.Thus, no systematic uncertainty due to pileup is included in the final results.

Uncertainty in the ratio of efficiencies
The statistical uncertainty in the measurement of ε 1 /ε 2 from the simulation, owing to the finite size of the MC sample, is taken as a systematic uncertainty, as shown in Table 3.
Since the analysis relies on photon conversions, the effect of a possible incorrect simulation of the tracker detector material is estimated.Two modified material scenarios, i.e., special detector geometries prepared for this purpose, in which the total mass of the silicon tracker varies by up to 5% from the reference geometry, are used to produce new MC simulation samples [20].With these models, local variations of the radiation length with respect to the reference simulation can be as large as +8% and −3%.No significant difference in the ratio of efficiencies is observed and the corresponding systematic uncertainty is taken to be negligible.
Several choices of the generated p T (χ c ) spectrum are investigated.In particular, the use of the measured J/ψ spectrum [11] gives values that are compatible with the default ψ(2S) spectrum used for the final result.The choice of the spectrum affects the values of ε 1 /ε 2 only inasmuch as we perform an average measurement in each bin of p T (J/ψ), and the size of these bins is finite.We choose to assign a conservative systematic uncertainty by comparing the values of ε 1 /ε 2 obtained with the ψ(2S) spectrum with those obtained in the case where the p T (χ c ) spectrum is taken to be constant in each p T bin.The corresponding systematic uncertainties are given in Table 3.

χ c polarization
The polarizations of the χ c1 and χ c2 are unknown.Efficiencies are estimated under the assumption that the two states are unpolarized.If the χ c states are polarized, the resulting photon angular distribution and transverse momentum distributions will be affected.This can produce a change in the photon efficiency ratio ε 1 /ε 2 .
In order to investigate the impact of different polarization scenarios on the ratio of the effi-ciencies, we reweight the unpolarized MC distributions to reproduce the theoretical χ c angular distributions [21,22] for different χ c polarizations.We measure the efficiency ε 1 /ε 2 for the χ c1 being unpolarized or with helicity m χ c1 = 0, ±1, in combination with the χ c2 being unpolarized or having helicity m χ c2 = 0, ±2 in both the helicity and Collins-Soper [23] frames.The ratio of efficiencies for the cases involving m χ c2 = ±1 is between the cases with m χ c2 = 0 and m χ c2 = ±2.Tables 4 and 5 give the resulting ε 1 /ε 2 values for each polarization scenario in different J/ψ transverse momentum bins for the two frames, relative to the value of the ratio for the unpolarized case.These tables, therefore, provide the correction that should be applied to the default value of ε 1 /ε 2 in each polarization scenario and each range of transverse momentum.
Table 4: The efficiency ratio ε 1 /ε 2 for different polarization scenarios in which the χ c1 is either unpolarized or has helicity m χ c1 = 0, ±1 and the χ c2 is either unpolarized or has helicity m χ c2 = 0, ±2 in the helicity frame, relative to the unpolarized case.

Branching fractions
The measurement of the prompt χ c2 to χ c1 production cross section ratio is affected by the uncertainties in the branching fractions of the two states into J/ψ + γ.The quantity that is directly accessible in this analysis is R p , the product of the ratio of the χ c2 to χ c1 cross sections and the ratio of the branching fractions.
In order to extract the ratio of the prompt production cross sections, we use the value of 1.76 ± 0.10 for B(χ c1 → J/ψ + γ)/B(χ c2 → J/ψ + γ) as derived from the branching fractions and associated uncertainties reported in Ref. [9].

Results and discussion
The results of the measurement of the ratio R p and of the ratio of the χ c2 to χ c1 prompt production cross sections for the kinematic range p T (γ) > 0.5 GeV/c and |y(J/ψ)| < 1.0 are reported in Tables 6 and 7, respectively, for different ranges of p T (J/ψ).The first uncertainty is statistical, the second is systematic, and the third comes from the uncertainty in the branching fractions in the measurement of the cross section ratio.Separate columns are dedicated to the uncertainty derived from the extreme polarization scenarios in the helicity and Collins-Soper frames, by choosing from Tables 4 and 5 the scenarios that give the largest variations relative to the unpolarized case.These correspond to (m χ c1 , m χ c2 ) = (±1, ±2) and (m χ c1 , m χ c2 ) = (0, 0) for both the helicity and Collins-Soper frames.Figure 3 displays the results as a function of the J/ψ transverse momentum for the hypothesis of unpolarized production.The error bars represent the statistical uncertainties and the green bands the systematic uncertainties.
Our measurement of the ratio of the prompt χ c2 to χ c1 cross sections includes both directly produced χ c mesons and indirectly produced ones from the decays of intermediate states.To convert our result to the ratio of directly produced χ c2 to χ c1 mesons requires knowledge of the amount of feed-down from all possible short-lived intermediate states that have a decay mode into χ c2 or χ c1 .The largest known such feed-down contribution comes from the ψ(2S).
Using the measured prompt J/ψ and ψ(2S) cross sections in pp collisions at 7 TeV [16], the branching fractions for the decays ψ(2S) → χ c1,2 + γ [9], and assuming the same fractional χ c contribution to the total prompt J/ψ production cross section as measured in pp collisions at 1.96 TeV [24], we estimate that roughly 5% of both our prompt χ c1 and χ c2 sample comes from ψ(2S) decays.The correction in going from the prompt ratio to the direct ratio is about 1%.In comparing our results with the theoretical predictions described below, we have not attempted to correct for this effect since the uncertainties on the fractions are difficult to estimate, the correction is much smaller than the statistical and systematic uncertainties, and our conclusions on the comparisons with the theoretical predictions would not be altered by a correction of this magnitude.
We compare our results with theoretical predictions derived from the k T -factorization [6] and NRQCD [7] calculations in Fig. 4. The k T -factorization approach predicts that both χ c1 and χ c2 are produced in an almost pure helicity-zero state in the helicity frame.Therefore, in our comparison, we apply the corresponding correction on the ratio of efficiencies from Table 4, amounting to a factor of 0.73, almost independent of p T .The theoretical calculation is given in the same kinematic range (p T (γ) > 0.5 GeV/c, |y(J/ψ)| < 1.0) as our measurement.There is no information about the χ c polarization from the NRQCD calculations, so we use the ratio of efficiencies estimated in the unpolarized case for our comparison.The prediction is given in the kinematic range p T (γ) > 0 GeV/c, |y(J/ψ)| < 1.0.We use the same MC simulation described in Section 5 to derive the small correction factor (ranging from 0.98 to 1.02 depending on p T , with uncertainties from 1 to 4%) needed to extrapolate the phase space of our measurement to the one used for the theoretical calculation.The uncertainty in the correction factor stemming from the assumption of the χ c transverse momentum distribution is added as a systematic uncertainty.The values of R p after extrapolation are shown in Table 8.The comparison of our measurements with the k T -factorization and NRQCD predictions are shown in the leftand rightplots of Fig. 4, respectively.The k T -factorization prediction agrees well with the trend of R p versus transverse momentum of the J/ψ, but with a global normalization that is higher by about a factor two with respect to our measurement.It is worth noting that this calculation assumes the same wave function for the χ c1 and the χ c2 .On the other hand, the NRQCD prediction is compatible with our results within the experimental and theoretical uncertainties, though, since predictions for χ c1 or χ c2 polarizations were not provided, the level of agreement can vary considerably.
A direct comparison of our results with previous measurements, in particular from [4] and [5], is not straightforward, because of the different conditions under which they were carried out.Specifically, there are differences in the kinematical phase space considered and, in the case of [4], in the initial-state colliding beams and center-of-mass energy used.However, with these caveats, a direct comparison shows that the three results are compatible within their uncertainties.In particular, all three results confirm the trend of a decreasing ratio of χ c2 to χ c1 production cross sections as a function of p T (J/ψ), under the assumption that the χ c2 and χ c1 polarizations do not depend on p T (J/ψ).

Summary
Measurements have been presented of the ratio as a function of the J/ψ transverse momentum up to relaxp T (J/ψ) = 25 GeV/c for the kinematic range p T (γ) > 0.5 GeV/c and |y(J/ψ)| < 1.0 in pp collisions at √ s = 7 TeV with a data sample corresponding to an integrated luminosity of 4.6 fb −1 .The corresponding values for the ratio of the χ c2 to χ c1 production cross sections have been determined.
The results have also been shown after extrapolating the photon acceptance down to zero p T .The effect of several different χ c polarization scenarios on the photon reconstruction efficiency has been investigated and taken into account in the comparison of the experimental results with two recent theoretical predictions.This is among the most precise measurements of the χ c production cross section ratio made in hadron collisions, and extends the explored J/ψ p T range of previous results.These measurements will provide important input to and constraints on future theoretical calculations of quarkonium production, as recently discussed in [25] for the bottomonium family.σ(χ c1 )B(χ c1 ) for the given p T (J/ψ) ranges in the fiducial kinematic region p T (γ) > 0.5 GeV/c, |y(J/ψ)| < 1.0, assuming unpolarized χ c production.The first uncertainty is statistical and the second is systematic.The last two columns report the additional uncertainties derived from the extreme polarization scenarios in the helicity (HX) and Collins-Soper (CS) frames.Figure 3: Ratio of the χ c2 to χ c1 production cross sections (circles) and ratio of the cross sections times the branching fractions to J/ψ + γ (squares) as a function of the J/ψ transverse momentum with the hypothesis of unpolarized production.The error bars correspond to the statistical uncertainties and the green band corresponds to the systematic uncertainties.For the cross section ratios, the 5.6% uncertainty from the branching fractions is not included.σ(χ c1 )B(χ c1 ) values with theoretical predictions from the k T -factorization [6] (left) and NRQCD [7] (right) calculations (solid red lines).The error bars and green bands show the experimental statistical and systematic uncertainties, respectively.The measurements in the leftplot use an acceptance correction assuming zero helicity for the χ c , as predicted by the k T -factorization model.The measurements in the rightplot are corrected to match the kinematic range used in the NRQCD calculation and assume the χ c are produced unpolarized.The measurements assuming two different extreme polarization scenarios are shown by the long-dashed blue and short-dashed green lines in the plot on the right.The 1standard-deviation uncertainties in the NRQCD prediction, originating from uncertainties in the color-octet matrix elements, are displayed as the dotted red lines.Table 7: Measurements of σ(χ c2 )/σ(χ c1 ) for the given p T (J/ψ) ranges derived using the branching fractions from Ref. [9], assuming unpolarized χ c production.The first uncertainty is statistical, the second is systematic, and the third from the branching fraction uncertainties.The last two columns report the uncertainties derived from the extreme polarization scenarios in the helicity (HX) and Collins-Soper (CS) frames.σ(χ c1 )B(χ c1 ) for the given p T (J/ψ) ranges after extrapolating the measurement to the kinematic region p T (γ) > 0 and assuming unpolarized χ c production.The first uncertainty is statistical and the second is systematic.The last column reports the largest variations due changes in the assumed χ c polarizations.

Figure 1 :
Figure 1: Distribution of the conversion radius for the χ c photon candidates.

Figure 2 :
Figure 2: The distribution of the variable Q = m µµγ − m µµ + m J/ψ for χ c candidates with p T (J/ψ) ranges shown in the figures.The line shows the fit to the data.

Table 1 :
Ratio of efficiencies ε 1 /ε 2 as a function of the J/ψ transverse momentum from MC simulation.The uncertainties are statistical only.

Table 3 :
Relative systematic uncertainties on R p for different ranges of J/ψ transverse momentum from different sources and the total uncertainty.

Table 6 :
Measurements of