Inclusive J/$\psi$ production in Xe-Xe collisions at $\sqrt{s_{\rm NN}}$ = 5.44 TeV

Inclusive J/$\psi$ production is studied in Xe-Xe interactions at a centre-of-mass energy per nucleon pair of $\sqrt{s_{\rm NN}}= 5.44$ TeV, using the ALICE detector at the CERN LHC. The J/$\psi$ meson is reconstructed via its decay into a muon pair, in the centre-of-mass rapidity interval $2.5<y<4$ and down to zero transverse momentum. In this Letter, the nuclear modification factors $R_{\rm AA}$ for inclusive J/$\psi$, measured in the centrality range 0-90% as well as in the centrality intervals 0-20% and 20-90% are presented. The $R_{\rm AA}$ values are compared to previously published results for Pb-Pb collisions at $\sqrt{s_{\rm NN}}= 5.02$ TeV and to the calculation of a transport model. A good agreement is found between Xe-Xe and Pb-Pb results as well as between data and the model.

A step in this direction is performed in this Letter, where first results on J/ψ production at LHC energies in Xe-Xe, a collision system (A Xe = 129) lighter than Pb-Pb (A Pb = 208), are presented. Data were collected by the ALICE Collaboration at the centre-of-mass energy per nucleon pair √ s NN = 5.44 TeV, during a short run carried out at the end of 2017. Due to the limited integrated luminosity, L int ∼ 0.34 µb −1 , the statistical uncertainties are significantly larger than those of the Pb-Pb results [10], but nevertheless allow a meaningful comparison between the two systems, in terms of the nuclear modification factor R AA . This quantity is obtained as the ratio between the production yields in nucleus-nucleus collisions and the corresponding proton-proton (pp) cross section, normalized to the nuclear thickness function T AA [20]. Values of R AA smaller (larger) than unity indicate suppression (enhancement) effects for the particle under study. The results shown in this Letter correspond to the centre-of-mass rapidity range 2.5 < y < 4, are integrated over transverse momentum (p T ) and were obtained by studying the J/ψ → µ + µ − decay channel. The nuclear modification factor is studied as a function of the centrality of the collision [21], expressed as a percentage of the hadronic Xe-Xe cross section. The results correspond to inclusive J/ψ production, which is the sum of a prompt component (directly produced J/ψ and feed-down from other charmonium states) and a non-prompt component, due to the decay of particles containing a b quark.
ALICE is the LHC experiment dedicated to the study of nuclear collisions, and is described in detail in Refs. [22,23]. The main detector used in this analysis is a muon spectrometer [24], covering the pseudorapidity range −4 < η < −2.5 1 . It includes tracking and trigger chambers, and reconstructs muons with p T larger than a given threshold, which is set at the trigger level. In addition, the V0 [25], a set of scintillator detectors covering 2.8 < η < 5.1 and −3.7 < η < −1.7, is used to define the minimum bias (MB) interaction trigger via a coincidence of signals at positive and negative η values. The V0 is also used for the centrality estimate via a fit of the distribution of the total signal amplitudes in the framework of the Glauber model [21]. The reconstruction of the primary collision vertex is carried out in the two layers of the Silicon Pixel Detector (SPD), the innermost part of the Inner Tracking System of the experiment [26], covering |η| < 2 and |η| < 1.4 respectively. Finally, rejection of non-hadronic Xe-Xe collisions is performed using the Zero Degree Calorimeters (ZDC) [27], which identifies electromagnetic interactions, while the V0 detects beam-gas collisions occurring outside the nominal interaction point region.
The data analyzed in this Letter are taken with a trigger formed by the coincidence of the MB trigger 1 In the ALICE reference frame, the muon spectrometer covers a negative η range and consequently a negative y range. We have chosen to present our results with a positive y notation.
Inclusive J/ψ production in Xe-Xe collisions at √ s NN = 5.44 TeV ALICE Collaboration signal and of at least one muon triggered in the muon spectrometer, with a p T = 0.5 GeV/c threshold. The definition of the trigger is less restrictive than the one usually adopted for Pb-Pb data taking (1 GeV/c threshold and two detected muons), due to the much smaller instantaneous luminosity for Xe-Xe collisions. Standard selection criteria [10] are then applied to such events and to the muon candidates.
In particular, it is required (i) that two opposite-sign tracks reconstructed in the tracking chambers of the muon spectrometer are matched to track segments in the trigger system, (ii) that both muons belonging to the pair (dimuon) have −4 < η µ < −2.5, and (iii) that their transverse position R abs at the end of the hadron absorber of the muon spectrometer satisfies the condition 17.6 < R abs < 89.5 cm. Finally, the reconstructed dimuon should lay in the fiducial rapidity region of the muon spectrometer, 2.5 < y < 4.
The nuclear modification factor R AA for the collision system under study is defined, for the centrality interval i, as where N i J/ψ is the number of detected J/ψ in the i-th centrality interval, BR J/ψ→µ + µ − = (5.96 ± 0.03)% is the branching ratio of the dimuon decay channel [28], N i MB is the number of MB events corresponding to the analyzed triggered event sample, Aε i is the product of the detector acceptance times the reconstruction efficiency, T i AA is the average nuclear thickness function [29], and σ pp J/ψ is the inclusive J/ψ cross section for pp collisions, at the same energy and in the same kinematic range as the Xe-Xe data. Results are given for the centrality interval 0-90% and for the two sub-intervals 0-20% and 20-90%.
Except for the determination of σ pp J/ψ , the other quantities entering the definition of R AA are evaluated following the same procedure used for the analysis of the Pb-Pb data sample and detailed in Ref. [10].
The extraction of N J/ψ is performed with two different approaches. In the first, the raw oppositesign dimuon invariant mass distribution is fitted with a superposition of resonance and background shapes [30], the former being tuned to MonteCarlo (MC) simulations and the latter corresponding to empirical functions. In the second, the background is estimated via a mixed-event invariant mass distribution, obtained from the collected sample of muon-triggered events and subtracted from the raw spectrum [9]. The resulting distribution is then fitted with the sum of a resonance shape and a continuum function accounting for the small residual background component. Due to the low statistical significance of the present data sample, the width of the J/ψ meson, which is usually kept as a free parameter in the invariant mass fits, is fixed to σ J/ψ = 70 MeV/c 2 , corresponding to the value of this quantity obtained in previous analyses [10,31,32]. For each of the two approaches, several fits were performed varying the fit mass range, the signal and background shapes and the J/ψ width by ±1 MeV/c 2 . The obtained value for the centrality interval 0-90% is N J/ψ = 241 ± 47(stat.) ± 26(syst.), where the central value and the statistical uncertainty correspond to the average of the fit results and to the average of the corresponding statistical uncertainties, respectively. The systematic uncertainty is obtained as the root mean square of the distribution of the N J/ψ values obtained with the various fits. The corresponding values for the 0-20% and 20-90% centrality sub-intervals are N J/ψ = 175 ± 42(stat.) ± 23(syst.) and N J/ψ = 77 ± 20(stat.) ± 7(syst.), respectively. Figure 1 shows as an example the results of two fits to the 0-90% Xe-Xe dimuon invariant mass distribution, corresponding to fitting the raw spectrum (left panel) or the mixed-event background subtracted mass distribution (right panel).
The product of the acceptance times the reconstruction efficiency Aε for J/ψ is evaluated via a MC simulation, based on the GEANT3 transport model [33], which takes into account the alignment of the muon spectrometer detectors and their efficiency. The input p T and y distributions for the J/ψ meson cannot be tuned directly to data, due to the low integrated luminosity of the data sample. Therefore, Inclusive J/ψ production in Xe-Xe collisions at √ s NN = 5.44 TeV ALICE Collaboration the differential distributions measured in Pb-Pb collisions at √ s NN = 5.02 TeV [10] for the 20-40% centrality range are used. In this centrality interval, the average number of participant nucleons weighted by the corresponding number of nucleon-nucleon collisions, N w part PbPb,20−40% is equal, within ∼ 2%, to N w part XeXe,0−90% , estimated via a Glauber MC calculation. The reconstruction efficiency slightly depends on the collision centrality, due to the detector occupancy in the muon spectrometer. The effect was evaluated in the analysis of Pb-Pb events [10] by embedding the simulated J/ψ signal into real events corresponding to various centralities. For this analysis, starting from the Pb-Pb results, the decrease in Aε XeXe,0−90% with respect to a simulation containing only J/ψ is estimated to be 4.2% (values for 0-20% and 20-90% centrality ranges are 5.5% and 1.6%, respectively). The resulting value for 0-90% Xe-Xe collisions is Aε XeXe,0−90% = 0.228±0.009(syst.), with a negligible statistical uncertainty. The systematic uncertainty is related to the choice of the J/ψ rapidity and transverse momentum distributions and to the estimate of the trigger and tracking efficiencies. The first term is evaluated by choosing alternative input shapes corresponding to other Pb-Pb centrality ranges, while for the efficiencies the procedure used in Ref.
The normalization factor N MB is evaluated by multiplying the number of opposite-sign dimuon triggers by a factor F norm , corresponding to the inverse of the probability of having a triggered muon in a MB event. This quantity is computed from the event trigger input information and the level-0 trigger mask. The procedure and the evaluation of the systematic uncertainty are described in Ref.
The reference cross section for the calculation of R AA is obtained starting from the measured value of the inclusive J/ψ cross section in pp collisions at √ s = 5.   Finally, a systematic uncertainty on the definition of the centrality intervals is evaluated by varying the Inclusive J/ψ production in Xe-Xe collisions at √ s NN = 5.44 TeV ALICE Collaboration value of the V0 signal amplitude corresponding to 90% centrality by ±0.5% and recalculating correspondingly the centrality intervals. Table 1 shows a summary of the systematic uncertainties for the R AA measurement for the three analyzed centrality ranges. The main contributions come from the estimate of T AA and from the signal extraction. The former is dominated by the uncertainty on the surface thickness of the Xe nucleus. The latter, being estimated in a data-driven way as detailed above, may suffer from the statistical limitations of the data sample. The quoted values can therefore be considered to be a conservative estimate.
The p T -integrated nuclear modification factor for inclusive J/ψ production in Xe-Xe collisions at √ s NN = 5.44 TeV, measured in 2.5 < y < 4 and in the 0-90% centrality range, is R AA = 0.54 ± 0.11(stat.) ± 0.09(syst.). This value can be compared with the corresponding one for Pb-Pb collisions at √ s NN = 5.02 TeV, R PbPb AA = 0.65±0.01(stat.)±0.04(syst.) [10]. Their ratio amounts to 0.84±0.16(stat.)±0.14(syst.), showing that the two values agree within about 0.7σ . Following the approach of Ref. [9], it can be shown that the Xe-Xe nuclear modification factor for prompt J/ψ could be up to 10% higher (lower) than the inclusive R AA if the non-prompt J/ψ component from the decays of hadrons containing a b quark is not (completely) suppressed. In Fig. 2 the R AA values for 0-20% and 20-90% Xe-Xe collisions are plotted, and compared with the centrality dependence of the nuclear modification factor for Pb-Pb collisions [10]. The latter shows, after a decrease up to N part ∼ 100, a saturation at R AA ∼ 0.65−0.7 towards more central events, and the two Xe-Xe points are found to be in agreement, within their larger uncertainties, with the Pb-Pb results. The Xe-Xe and Pb-Pb results are also compared with the calculation of a transport model by Du, Zhao and Rapp [13,14]. A close similarity of the predicted suppression patterns for Pb-Pb and Xe-Xe is observed, which fairly reproduces the experimental results.
In summary, we have measured inclusive J/ψ production in Xe-Xe collisions at √ s NN = 5.44 TeV. Results on the nuclear modification factors were given for various centrality selections and compared to corresponding results for Pb-Pb collisions at √ s NN = 5.02 TeV and to a theoretical model. Within the experimental uncertainties, a good agreement is found between the R AA measured in the two systems and with the calculation. These results show that the relative contribution of suppression and regeneration processes is similar for collisions producing similar N part values from different collision systems.

Acknowledgements
The ALICE Collaboration would like to thank all its engineers and technicians for their invaluable contributions to the construction of the experiment and the CERN accelerator teams for the outstanding performance of the LHC complex. The ALICE Collaboration gratefully acknowledges the resources and support provided by all Grid centres and the Worldwide LHC Computing Grid (WLCG) collaboration. The ALICE Collaboration acknowledges the following funding agencies for their support in building and running the ALICE detector:   [13] X. Zhao and R.   [33] R. Brun, F. Carminati, and S. Giani, "GEANT Detector Description and Simulation Tool," CERN Program Library Long Writeup CERN-W5013 (1994) .
Inclusive J/ψ production in Xe-Xe collisions at √ s NN = 5.44 TeV ALICE Collaboration

A The ALICE Collaboration
Inclusive J/ψ production in Xe-Xe collisions at √ s NN = 5.44 TeV ALICE Collaboration