arxiv: v1 [physics.data-an] 12 Sep 2018

Transcription

1 Correcting Differential Correlation functions for Collision Centrality bin width effects Victor Gonzalez,,, Ana Marin, Pedro Ladron de Guevara,, 3 Jinjin Pan,, Suit Basu,, and Claude A. Pruneau Universidad Coplutense de Madrid, Spain GSI Helholtzzentru für Schwerionenforschung, Research Division and ExtreMe Matter Institute EMMI, Darstadt, Gerany 3 Instituto de Fisica, Universidad Nacional Autonoa de Mexico, Mexico City, Mexico Departent of Physics and Astronoy, Wayne State University (Dated: Septeber, 08 arxiv: v [physics.data-an] Sep 08 We present a technique to correct nuber and oentu differential correlation functions R, G and P, for collision centrality bin width averaging. he technique is based on the hypothesis that single- and pair- probability distributions vary slower than single and pair ultiplicity of particles of interest as a function of collision centrality in A A collisions. We test the technique with siulations based on the HIJING and UrQMD odels and conclude it enables a precision better than % for particles in the kineatic range 0. p.0 GeV/c. PACS nubers: 5.75.Gz, 5.75.Ld,.60.Ky,.60.-k victor.gonzalez@uc.es jinjin.pan@cern.ch suits.basu@cern.ch

2 I. INRODUCION Measureents of correlation functions enable in-depth exploration of particle production echaniss in relativistic heavy-ion collisions. ogether with easureents of the nuclear odification factor, R AA, easured nuber correlation functions have provided strong evidence for the foration of a dense and opaque ediu, a Quark Gluon Plasa (QGP, in the idst of high energy heavy ion collisions [ 9]. Measureents of transverse oentu differential correlation functions have also been added to the experiental toolset [0, ]. For these as for differential nuber correlation functions, precision easureents require one accounts for a nuber of experiental conditions and artifacts []. In this paper, we consider corrections required to account for finite collision centrality bin width averaging. Inclusive and conditional particle densities are defined in sec. II. Centrality bin width corrections are first derived for noralized two-particle cuulants in sec. III. he technique is then extended to oentu correlation functions in sec. IV. In sec. V, we present estiates of the precision of the ethod based on siulations carried out with the HIJING [3] and UrQMD [, 5] particle event generators. his work is suarized in sec. VI. II. OBSERVABLE DEFINIIONS Whether analyzing pp, p Pb, or A A collisions, it is legitiate to classify events on the basis of the energy, E ref, or the charged particle ultiplicity,, easured in a reference acceptance Ω ref. Studies [6] have shown that in A A collisions, this ultiplicity aps rather narrowly onto the ipact paraeter b of the collisions, particularly in the case of id- to central-collisions. While such narrow apping does not exist for p A or pp interactions, it reains appropriate to classify collisions based on the ultiplicity as it provides soe easure of the oentu transfer and projectile kinetic energy dissipated in the collisions. In general, one indeed finds that the single and pair density easured in the fiducial acceptance Ω ref are functions of the reference ultiplicity. We thus define single- and two-particle conditional densities (i.e., for a given as ρ (η, ϕ, p = dn p dϕdηdp ( d N ρ (η, ϕ, p,, η, ϕ, p, = p, p, dϕ dη dp, dϕ dη dp, ( ypically, correlation analyses are carried out over fixed p ranges [p,in, p,ax ]. One thus also defines conditional single- and two-particle densities integrated over p as ρ (η, ϕ = ρ (η, ϕ, η, ϕ = ρ (η, ϕ, p dp, (3 ρ (η, ϕ, p,, η, ϕ, p, dp, dp, ( where (here and in the following it is understood that the integrals in p are taken over the fixed range [p,in, p,ax ]. In practice, it is usually not possible (or eaningful to easure the densities ρ (η, ϕ and ρ (η, ϕ, η, ϕ for unit resolution in. One ust then evaluate the densities within finite width bins of ultiplicity [ in,k, ax,k ] (where k =,..., K, represents one of K centrality bins used in the analysis as weighted average of the densities across the bins according to with ρ (k (η, ϕ = Q k ρ (k (η, ϕ, η, ϕ = Q k Q k = ax,k q(ρ (η, ϕ, (5 ax,k q(ρ (η, ϕ, η, ϕ, (6 ax,k q( (7

3 3 and q( representing the probability of events with ultiplicity in the reference acceptance. III. NORMALIZED CUMULANS: R wo-particle noralized cuulants [] are defined according to R (η, ϕ, η, ϕ = ρ (η, ϕ, η, ϕ ρ (η, ϕ ρ (η, ϕ ρ (η, ϕ ρ (η, ϕ (8 and for the basis of any two-particle correlation analyses. Experientally, in the study of A A collisions, it is coon to average R across large collision centrality bins. Given the width of such bins ay depend on the specificities of a given experient (e.g., the size of the reference acceptance, noralized cuulants R are thus bin-width dependent. A correction procedure is then required to account for the finite width of the centrality bins. In principle, one would like to set the collision centrality bins in ters of collision ipact paraeter ranges (alternatively the nuber of wounded nucleons. In practice, only proxies of the ipact paraeter are at best possible and one usually expresses the bins directly in ters of such proxies or in ters of the fractional cross-section easured within a centrality (ultiplicity bin. Let, the ultiplicity easured in the reference acceptance Ω ref represent such a proxy. As for densities, one can thus define conditional noralized cuulants according to R (η, ϕ, η, ϕ = ρ (η, ϕ, η, ϕ. (9 ρ (η, ϕ ρ (η, ϕ However, easureents of R (η, ϕ, η, ϕ for unit resolution are typically not possible due to liitations associated with the finite size of datasets, CPU tie, or storage considerations. It is also ipractical to report differential correlators R (η, ϕ, η, ϕ for a very large nuber of values of. In practice, one thus seeks to report correlators R averaged over reference ultiplicity bins [ in,k, ax,k ]. As for averages of densities defined in sec. II, this is noinally achieved according to R (k (η, ϕ, η, ϕ = Q k ax,k q(r (η, ϕ, η, ϕ, (0 where q( represents the probability of events with ultiplicity in the reference acceptance. deterine R (k according to One ay then R (Bin,k (η, ϕ, η, ϕ = ρ (Bin,k (η, ϕ, η, ϕ, ( ρ (Bin,k (η, ϕ ρ (Bin,k (η, ϕ in which ρ (Bin,k and ρ (Bin,k are single- and pair-densities easured directly using finite width bins in. In the absence of biases or event detection inefficiencies, one has he quantities R (Bin,k and ρ (Bin,k (η, ϕ = ρ (k (η, ϕ, ( ρ (Bin,k (η, ϕ, η, ϕ = ρ (k (η, ϕ, η, ϕ. (3 R (k are thus clearly distinct: R (Bin,k (η, ϕ, η, ϕ = R (k (η, ϕ, η, ϕ = Q k ax,k q(ρ (η, ϕ, η, ϕ [ ] [ ax,k ] ax,k, ( Q k q(ρ (η, ϕ Q k = in,k q( ρ (η, ϕ ax,k ρ (η, ϕ, η, ϕ q(. (5 Q k ρ = (η, ϕ ρ (η, ϕ in,k

4 o quantify the difference, we write the particle densities in ters of single- and two-particle probability densities, P (η, ϕ and P (η, ϕ, η, ϕ according to where, by definition, P and P respectively satisfy ρ (η, ϕ = n P (η, ϕ, (6 ρ (η, ϕ, η, ϕ = n(n P (η, ϕ, η, ϕ, (7 dϕdη P (η, ϕ = (8 dϕ dη dϕ dη P (η, ϕ, η, ϕ =, (9 as probability densities. he quantities n and n(n are the ean nuber of particles and the ean nuber of pairs of particles in the acceptance of the easureent at a given reference ultiplicity ; P (η, ϕ is the probability of finding a particle at η, ϕ when the ultiplicity is, and P (η, ϕ, η, ϕ is the joint probability of easuring particles at η, ϕ and η, ϕ when the ultiplicity is. In general, one expects n and n(n to scale approxiately linearly and quadratically, respectively, with. Let us assue that the shape of P (η, ϕ and P (η, ϕ, η, ϕ change little through a centrality bin. One can then write with R (Bin,k (η, ϕ, η, ϕ = α P (η, ϕ, η, ϕ, (0 P (η, ϕ P (η, ϕ α = ax,k Q k ( q( n(n Q k ax,k q( n ( and P (η, ϕ and P (η, ϕ, η, ϕ are the bin-width averaged values of the single- and two-particle probability densities. Siilarly, one can also write R (k (η, ϕ, η, ϕ = β P (η, ϕ, η, ϕ, ( P (η, ϕ P (η, ϕ with β = ax,k q( n(n Q k n. (3 One thus find that for sufficiently narrow bins (such that P and P can be considered approxiately invariant within, R ( (k (η, ϕ, η, ϕ = βα R (Bin,k (η, ϕ, η, ϕ + ( he above expression assues that the event detection efficiency is either unity or constant across a bin k. If the event detection efficiency, ε(, varies across [ in,k, ax,k ], and the densities P and P are not biased by this dependence, then the coefficients α and β ay be written as α = ax,k q( Q k ( Q k ax,k β = ax,k q( Q k ε( ε( n(n, (5 q( ε( n where q( is now the observed (uncorrected probability distribution of ultiplicity. n(n n, (6

5 5 IV. RANSVERSE MOMENUM CORRELAORS Several fors of transverse oentu correlators have been proposed and reported in the recent literature [0,, 7 ]. In this work, we focus on the G correlator proposed by Gavin et al. [, 8 0] to study transverse oentu current correlations and the P correlator designed to be sensitive to transverse oentu fluctuations [0,, ]. A. G Correlator At fixed reference ultiplicity,, the differential transverse oentu correlator G can be written as [8] G (η, ϕ, η, ϕ = dp, dp, p, p, ρ (η, ϕ, p,, η, ϕ, p, ρ (η, ϕ ρ (η, ϕ p (η, ϕ p (η, ϕ (7 where p (η, ϕ = dp p ρ (η, ϕ, p ρ (η, ϕ is the event inclusive average particle transverse oentu at η, ϕ. Ḡ (k (η, ϕ, η, ϕ proceeds as for R and one writes Introducing Ḡ (k (η, ϕ, η, ϕ = Q k ax,k (8 Evaluation of the centrality-bin averaged q( G (η, ϕ, η, ϕ. (9 S (η, ϕ = p dp ρ (η, ϕ, p, (30 = n P p (η, ϕ, (3 S (η, ϕ, η, ϕ = p, dp, p, dp, ρ (η, ϕ, p,, η, ϕ, p,, (3 with = n(n P pp (η, ϕ, η, ϕ, (33 P p (η, ϕ = P pp (η, ϕ, η, ϕ = p dp P (η, ϕ, p, (3 p, dp, p, dp, P (η, ϕ, p,, η, ϕ, p,, (35 Eq. (7 is written as G (η, ϕ, η, ϕ = S (η, ϕ, η, ϕ ρ (η, ϕ ρ (η, ϕ S (η, ϕ S (η, ϕ ρ (η, ϕ ρ (η, ϕ = n(n n P pp (η, ϕ, η, ϕ P (η, ϕ P (η, ϕ Pp P (η, ϕ (η, ϕ (36 P p (η, ϕ P (η, ϕ. (37 If the ratios P p /P and P pp /P P have a odest dependence on the reference ultiplicity (within a bin k, then it is legitiate to replace the by averages and one gets G (η, ϕ, η, ϕ = n(n n P pp(k (η, ϕ, η, ϕ (η (k, ϕ P (η, ϕ P p(k (η, ϕ (η, ϕ P p(k (η, ϕ (η, ϕ. (38

6 6 he centrality-bin averaged correlator Ḡ(k (η, ϕ, η, ϕ ay then be written as (η, ϕ, η, ϕ = β P pp (k (η, ϕ, η, ϕ (η (k, ϕ P (η, ϕ Ḡ (k P p(k (η, ϕ (η, ϕ P p(k (η, ϕ (η, ϕ. (39 However, if it is not possible to carry out the analysis in fine (unit bins of, the nuerators and denoinators of G ust be separately averaged over the range [ in,k, ax,k ]. Assuing the ratios P p /P and P pp /P P have a odest dependence on the reference ultiplicity, one gets (η, ϕ, η, ϕ = α P pp (k (η, ϕ, η, ϕ (η (k, ϕ P (η, ϕ G (Bin,k Identifying p (η, ϕ (Bin,k = ay be deterined as P p(k P p(k (η, ϕ (η, ϕ P p(k (η, ϕ (η, ϕ. (0 (η, ϕ/ (η, ϕ, one finally gets that the desired correlator Ḡ(k (η, ϕ, η, ϕ Ḡ (k (η, ϕ, η, ϕ = βα ( G (Bin,k (η, ϕ, η, ϕ + p (η, ϕ (Bin,k p (η, ϕ (Bin,k p (η, ϕ (Bin,k p (η, ϕ (Bin,k, ( =βα G (Bin,k (η, ϕ, η, ϕ + ( βα p (η, ϕ (Bin,k p (η, ϕ (Bin,k. he P correlator is defined as B. P Correlator P (η, ϕ, η, ϕ = p, dp, p, dp, ρ (η, ϕ, p,, η, ϕ, p, p, ( ρ (η, ϕ, η, ϕ where p,i p,i p and the p noralization insures P is diensionless. Introducing (3 one next verifies that P is insensitive to ultiplicity fluctuations [0,, ]. P p p (η, ϕ, η, ϕ = p, dp, p, dp, P (η, ϕ, p,, η, ϕ, p, (3 Factorizing the two-particle density as in (7, one gets P (η, ϕ, η, ϕ = ( P P p n(n P p p (η, ϕ, η, ϕ. ( n(n P (η, ϕ, η, ϕ If the ratios P p /P and P p p /P have a odest dependence on the reference ultiplicity (within a bin k, then it is legitiate to replace the by averages and one gets P (η, ϕ, η, ϕ = ( P(k P p(k P p p (k (η, ϕ, η, ϕ (η, (5, ϕ, η, ϕ which is indeed independent of the nuber of pairs and thus ultiplicity fluctuations, since the functions P (η, ϕ, η, ϕ are probability densities. he centrality-bin averaged value of P is calculated according (η, ϕ, η, ϕ = ( P(k P p(k P p p (k (η, ϕ, η, ϕ (η, (6, ϕ, η, ϕ

7 7 R 8 6 HIJING Pb--Pb.76 ev R % 0-% -3% -5% 0. < p <.0 GeV ϕ < average corrected.5 R R FIG.. Longitudinal projection of the two-particle noralized cuulants R (left and R (right for the ost central Pb Pb collisions produced with the HIJING event generator. op panels show uncorrected results for the 0-5% centrality bin together with those fro the 0-%, -3%, and -5% centrality bins, iddle panels show corrected 0-5% centrality bin results copared with the weighted average of those fro the 0-%, -%, -3%, 3-% and -5% centrality bins, while botto panels show their ratio. while the ratio of bin averages is One thus finds that P (Bin,k (η, ϕ, η, ϕ = ( P(k P p(k P p p (k (η, ϕ, η, ϕ (η. (7, ϕ, η, ϕ (η, ϕ, η, ϕ = P (Bin,k (η, ϕ, η, ϕ, (8 featuring a unitary correction factor. he P correlator is thus indeed not affected by the width of the ultiplicity bins so long as the probability densities P and P exhibit only odest shape variations within those bins. V. SIMULAIONS We tested the accuracy of the centrality bin width correction ethods defined by Eqs. ( and ( using HIJING and UrQMD siulations of 5% ost central central Pb Pb collisions at s NN =.76 ev. Overall, 65,000 and 00,000 Pb Pb collisions fro HIJING and UrQMD were used for this study. he collision centrality deterination

8 8 R 8 UrQMD Pb--Pb.76 ev 6 R % 0-% -3% -5% 0. < p <.0 GeV ϕ <.0 0. average corrected 0. R 3.8 R η. 0 η FIG.. Longitudinal projection of the two-particle noralized cuulants R (left and R (right for 5% ost central Pb Pb collisions produced with the UrQMD event generator. op panels show uncorrected results for the 0-5% centrality bin together with those fro the 0-%, -3%, and -5% centrality bins, iddle panels show corrected 0-5% centrality bin results copared with the weighted average of those fro the 0-%, -%, -3%, 3-% and -5% centrality bins, while botto panels show their ratio. for HIJING siulations utilizes the ethod used in ALICE [3, ] which is based on the nuber of charged particles in the pseudo-rapidity range.8 < η < 5. and.7 < η <.7 corresponding to the acceptance of the ALICE V0 detectors. For UrQMD siulations the centrality selection is based on sharp cuts in the ipact paraeter, as given in [3]. We first deterined the R, G, and P correlators for charged particle cobinations (+, +, and ++ as -diensional functions based on generator level charged hadrons in the ranges 0. p.0 GeV/c and η < 0.8 using Eqs. (, (, and (8. Fro these, we calculated the charge independent and charge dependent particle cobinations according to C = [ C (+ + C ( + + C ( + C (++] (9 C = [ C (+ + C ( + C ( C (++]. (50 where C stands for any of the R, G, and P correlators. Projections of these correlators onto η and ϕ were coputed using the integer arithetic technique described in []. Figure presents η projections of the R (left and R (right correlation functions based on 5% ost central HIJING events. he top panels of the figure present uncorrected results for 0-5% (blue, 0-% (green, -3% (red, and -5% (black. he botto panels displays 0-5% centrality results corrected with Eq. ( (blue copared to the weighted ean of the correlators obtained for 0-%, -%, -3%, 3-%, and -5% centralities (red and their ratio. One finds that the results corrected with Eq. ( agree with those obtained with the weighted ean within % for both R and R.

9 9 ((GeV/c G 3 HIJING Pb--Pb.76 ev 0. < p <.0 GeV ((GeV/c G HIJING Pb--Pb.76 ev 0. < p < 0.8 GeV ((GeV/c G 0-5% 0-% -3% -5% < p <.0 GeV ϕ <.0 ((GeV/c G ((GeV/c G 0. 5 ((GeV/c G average corrected FIG. 3. Longitudinal projection of the two-particle transverse oentu correlations G for default transverse oentu range particles (left and for the reduced (see text transverse oentu range particles (center and G for the default transverse oentu range particles (right for the 5% ost central Pb Pb collisions produced with the HIJING event generator. op panels show uncorrected results for the 0-5% centrality bin together with those fro the 0-%, -%, -3%, 3-% and -5% centrality bins, iddle panels show corrected 0-5% centrality bin results copared with the weighted average of those fro the 0-%, -%, -3%, 3-% and -5% centrality bins, while botto panels show their ratio. Siilar results are presented in Fig. for the UrQMD odel. In this case, we find the agreeent between the 0-5% corrected R correlators and the weighted ean of the R correlators obtained in 0-%, -%, -3%, 3-%, and -5% centralities to be also within % for both R and R. We thus conclude, based on the HIJING and UrQMD odel siulations presented, that Eq. ( enables reasonably accurate corrections of the R correlators in the context of these two odels. Given these odels provide relatively realistic representations of single and pair particle spectra, we infer the correction ethod ebodied in Eq. ( should provide reasonably reliable bin-width corrections of R correlation functions easured at any heavy ion collider. Figures 3 and present η projections of the G and G correlation functions based on 5% ost central events fro HIJING and UrQMD odels respectively. Siilarly as in previous figures, the top panels display uncorrected G correlators for 0-5% (blue, 0-% (green, -3% (red, and -5% (black while the botto panels display 0-5% centrality results corrected with Eq. ( (blue copared to the weighted ean of the G correlators obtained for 0- %, -%, -3%, 3-%, and -5% centralities (red and their ratio. Coparing the corrected results with the weighted averages for HIJING events one finds the within 7% for G and within % for G while the coparison for UrQMD events leaves the within % for both G and G. he discrepancy between the correction level achieved for G for collisions siulated with HIJING copared with those produced using UrQMD is hardly satisfactory, especially considering the % correction precision achieved, with both odels, for the R correlator. Given the transverse oentu enters explicitly into the expression of the G correlator, but not in R, we sought to deterine whether the p range considered ay influence the correction precision achievable with Eq. (. We thus repeated the above study for charged particles in the range GeV/c. Central panels of Fig. 3 display η projections of the G correlation function based on 5% ost central events fro HIJING odel for particle within the transverse oentu range GeV/c. As before, the top panel shows a coparison of the correlators calculated within the 0-5% centrality bin with those fro the 0-%, -3% and -5% centrality bins, while the botto panel displays 0-5% centrality results corrected with Eq. ( copared to the weighted ean of the values obtained in the 0-%, -%, -3%, 3-%, and -5% collision centrality bins. We find that the 0-5% centrality corrected result is now within % of the weighted ean of those obtained with finer centrality

10 0 ((GeV/c 3 UrQMD Pb--Pb.76 ev ((GeV/c % 0-% -3% -5% 0. < p <.0 GeV ϕ <.0 G G 0.. ((GeV/c. 0. G average corrected ((GeV/c G η. 0 η 0 η FIG.. Longitudinal projection of the two-particle transverse oentu correlations G (left and G (right for the 5% ost central Pb Pb collisions produced with the UrQMD event generator. op panels show uncorrected results for the 0-5% centrality bin together with those fro the 0-%, -%, -3%, 3-% and -5% centrality bins, iddle panels show corrected 0-5% centrality bin results copared with the weighted average of those fro the 0-%, -%, -3%, 3-% and -5% centrality bins, while botto panels show their ratio. bins. We exained the centrality dependence of the single- and two-particle probability densities (and their p weighted counterparts, Eqs. (3 and (35, and found that those obtained within the GeV/c range feature a quantitatively saller dependence of collision centrality than those integrated in the GeV/c range. he forer range consequently better satisfies the condition, iplicit in Eq. (0, that the densities be approxiately independent of collision centrality, within each centrality bin k. We additionally verified that these densities, calculated within the range GeV/c with UrQMD events, siilarly exhibit very little shape dependence on collision centrality. Our analysis confirs the relevance of the proper behavior of this agnitudes established in section IV A. We thus conclude that the finite centrality bin width correction, included in Eq. (, under the criteria established for its deduction, provides a reasonably accurate correction technique for G correlators to account for collision centrality average when using finite collision centrality bins. It is worth highlighting that, consistently with what was naively expected, the centrality bin width effect is copletely charge independent. We find, indeed, that the aplitude of both R and G are substantially affected by the wide centrality bin averaging while R and G only experience a inor shift owing ( to the very sall differences between the R LS and R US correlator dependence on particle densities and on the transverse oentu of the particles, and ( to the approxiately equal nuber of US and LS pairs observed at large s NN in the central rapidity range. For the sake of copleteness, the behavior of the correction procedures on the ϕ projections of R (for the whole p range and G (for the reduced p range correlators are shown in Fig. 5. We find that the correction procedures yield essentially perfect corrections of these ϕ projections for both correlators. Additionally, the applicability of Eq. (8 for P correlators is deonstrated in Fig. 6 where we show that the P correlator deterined with a 0-5%

11 R 8 HIJING Pb--Pb.76 ev 0. < p <.0 GeV 6 ((GeV/c G % 0-% -3% -5% average 0. < p < 0.8 GeV corrected R ((GeV/c 5 G.5. 0 ϕ (radians. 0 ϕ (radians FIG. 5. Aziuthal projections of the two-particle transverse oentu correlations R for the whole p range (left and G for the reduced p range (right for the 5% ost central Pb Pb collisions produced with the HIJING event generator. op panels show uncorrected results for the 0-5% centrality bin together with those fro the 0-%, -%, -3%, 3-% and -5% centrality bins, iddle panels show corrected 0-5% centrality bin results copared with the weighted average of those fro the 0-%, -%, -3%, 3-% and -5% centrality bins, while botto panels show their ratio. ABLE I. POL0 fit to the ratio 0-5% centrality corrected versus 0-%, -%, -3%, 3-%, and -5% centralities weighted ean. Model proj. R R P P G G HIJING η.000 ± ± ± ± ± ± 0.00 HIJING ϕ.000 ± ± ± ± ± ± 0.00 UrQMD η.00 ± ± ± ± ± ± 0.00 UrQMD ϕ 90 ± ± ± ± ± ± HIJING low p η 95 ± ± HIJING low p ϕ 9 ± ± 0.00 centrality bin is virtually identical to that obtained with a weighted average of the correlators calculated fro the 0-%, -%, -3%, 3-%, and -5% collision centrality bins. Finally, in order to obtain a quantitative assessent of the precision of the centrality bin width correction procedures ebodied by Eqs. (,8,, we fitted the ratios of the centrality bin width corrected correlators to those obtained with weighted averages of correlators obtained in fine width centrality bins, shown in Figs. -6, with a constant polynoial (POL0 in ROO [5]. Fitted values of the ratios, corresponding to ratios averaged across the η and ϕ fiducial ranges, are listed in able I. Overall, we find that the averaged ratios have values consistent with unity within statistical errors (i.e., within a ±σ, 96% confidence interval, thereby iplying that within the statistical precision achieved in this work, we can

12 P 0. HIJING Pb--Pb.76 ev 0. < p <.0 GeV 0. P % 0-% -3% -5% UrQMD Pb--Pb.76 ev ϕ < average corrected 0. P 0. P η FIG. 6. Longitudinal projection of the two-particle transverse oentu correlation P for the 5% ost central Pb Pb collisions produced with the HIJING event generator (left and the UrQMD event generator (right. op panels show uncorrected results for the 0-5% centrality bin together with those fro the 0-%, -%, -3%, 3-% and -5% centrality bins, iddle panels show corrected 0-5% centrality bin results copared with the weighted average of those fro the 0-%, -%, -3%, 3-% and -5% centrality bins, while botto panels show their ratio. conclude that the correction procedures do not introduce a significant bias to the correlation functions. Notable exceptions are the values obtained with both η and φ projections for G when this correlator is calculated for charged particles in the range 0. p.0 GeV/c with HIJING. We note, however, that the ratios of these projections is consistent with unity when the correlator is calculated in a narrower p range in line with the deanded behavior of the probability density distributions. A large discrepancy is also seeingly observed for the P correlator obtained with the UrQMD siulation. We note, however, that the overall agnitude of P predicted by UrQMD is a seventh of that obtained with HIJING events and correlation aplitudes are nearly vanishing, within statistical accuracy, across a wide η range. A proper evaluation of the correction is thus ore challenging in this case. VI. SUMMARY We described techniques to carry out corrections of R, G, and P correlators when these correlators ust be deterined based on wide collision centrality bins. We showed, based on HIJING and UrQMD siulations of Pb Pb collisions at s NN =.76 ev, that the correction techniques enable a precision of the order of % or better when using 5% wide centrality bins. Evidently, better precision should be possible when the centrality bins can be divided in saller intervals or if larger datasets are available.

13 3 ACKNOWLEDGEMENS he authors thank the ALICE collaboration for access to a saple of siulated HIJING events used in this analysis and the GSI Helholtzzentru für Schwerionenforschung for providing the coputational resources needed for producing the UrQMD events used in this analysis. his work was supported in part by the United States Departent of Energy, Office of Nuclear Physics (DOE NP, United States of Aerica under Award Nuber DE- FG0-9ER-073 and by Consejo Nacional de Ciencia y ecnologia (CONACY, through Fondo de Cooperacion Internacional en Ciencia y ecnologia (FONCY and Direccion General de Asuntos del Personal Acadeico (DGAPA, Mexico. [] J. Adas et al. (SAR, Nucl. Phys. A757, 0 (005, arxiv:nucl-ex/ [nucl-ex]. [] K. Adcox et al. (PHENIX, Nucl. Phys. A757, 8 (005, arxiv:nucl-ex/00003 [nucl-ex]. [3] I. Arsene et al. (BRAHMS, Nucl. Phys. A757, (005, arxiv:nucl-ex/0000 [nucl-ex]. [] B. B. Back et al., Nucl. Phys. A757, 8 (005, arxiv:nucl-ex/000 [nucl-ex]. [5] K. Aaodt et al. (ALICE, Phys. Rev. Lett. 05, 530 (00, arxiv:96 [nucl-ex]. [6] K. Aaodt et al. (ALICE, Phys. Rev. Lett. 05, 530 (00, arxiv:9 [nucl-ex]. [7] G. Aad et al. (ALAS Collaboration, Phys. Rev. Lett. 05, 5303 (00. [8] J. Ada et al. (ALICE, Phys. Rev. Lett. 6, 30 (06, arxiv:5.060 [nucl-ex]. [9] J. Ada et al. (ALICE, Phys. Rev. Lett. 6, 330 (06, arxiv:60.09 [nucl-ex]. [0] M. Shara and C. A. Pruneau, Phys. Rev. C79, 0905 (009, arxiv: [nucl-ex]. [] H. Agakishiev et al. (SAR, Phys. Lett. B70, 67 (0, arxiv:06.33 [nucl-ex]. [] S. Ravan, P. Pujahari, S. Prasad, and C. A. Pruneau, Phys. Rev. C 89, 0906 (0. [3] X.-N. Wang and M. Gyulassy, Phys. Rev. D, 350 (99. [] S. A. Bass et al., Prog. Part. Nucl. Phys., 55 (998, [Prog. Part. Nucl. Phys.,5(998], arxiv:nucl-th/ [nucl-th]. [5] M. Bleicher et al., J. Phys. G5, 859 (999, arxiv:hep-ph/ [hep-ph]. [6] R. Rogly, G. Giacalone, and J.-Y. Ollitrault, (08, arxiv: [nucl-th]. [7] J. Adas et al. (SAR, J. Phys. G3, L37 (006, arxiv:nucl-ex/ [nucl-ex]. [8] M. Abdel-Aziz and S. Gavin, Acta Phys. Hung. A5, 55 (006. [9] S. Gavin, G. Moschelli, and C. Zin, Proceedings, 3th Winter Workshop on Nuclear Dynaics (WWND 06: Guadeloupe, French West Indies, J. Phys. Conf. Ser. 736, 000 (06, arxiv: [nucl-th]. [0] S. Gavin, G. Moschelli, and C. Zin, Phys. Rev. C9, 09 (06, arxiv: [nucl-th]. [] S. Acharya et al. (ALICE Collaboration, Phys. Rev. Lett. 8, 630 (07, arxiv: [nucl-ex]. [] S. Acharya et al. (ALICE, (08, arxiv:805.0 [nucl-ex]. [3] B. Abelev et al. (ALICE, Phys. Rev. C88, 0909 (03, arxiv:30.36 [nucl-ex]. [] B. B. Abelev et al. (ALICE, Int. J. Mod. Phys. A9, 300 (0, arxiv:0.76 [nucl-ex]. [5] R. Brun and F. Radeakers, New coputing techniques in physics research V. Proceedings, 5th International Workshop, AIHENP 96, Lausanne, Switzerland, Septeber -6, 996, Nucl. Instru. Meth. A389, 8 (997.