in correlation analysis of binned data the quantity of interest is the covariance of observable @xmath0 between arbitrary bins @xmath1 and @xmath2 , given by @xmath3 where @xmath0 is the bin content and over - lines indicate averages over independent measurements . for example , in high - energy collisions between atomic nuclei
@xmath0 could represent the number of sub - atomic particles produced and detected for each collision , or _ event _
, within bins defined by particle 3-momentum .
for example , such bins could be constructed using the transverse momentum ( @xmath4 ) ( component of 3-momentum perpendicular to the direction of the colliding beams ) , the azimuth angle ( @xmath5 ) in the plane transverse to the beam , the pseudorapidity ( @xmath6 ) where @xmath7 and @xmath8 is the polar angle relative to the beam direction , the azimuthal angle difference ( @xmath9 ) for arbitrary particles 1 and 2 , and the pseudorapidity difference ( @xmath10 ) .
quantity @xmath11 is the number of particle - pairs in 2d bin @xmath12 and @xmath13 is calculated by averaging product @xmath11 over all collision events in the event collection .
pairs of particles from the same event are referred to as _ sibling _ pairs .
averages @xmath14 and @xmath15 are calculated using all events in the collection where quantity @xmath16 is the _
reference_. for a given set of measured events the covariance will have a unique numerical value .
the goal of this paper is to calculate the statistical error in this quantity by following and extending the error propagation and error reduction method of reid and trainor  @xcite .    in practical applications measured quantities @xmath17 , @xmath18 and @xmath11
are affected by experimental inefficiency , acceptance and contamination .
inefficiency and acceptance losses in single particle counts and some contamination effects are readily corrected via ratio @xmath19 .
however , two - particle inefficiencies , which occur when signals in the detectors from two particles are unresolved , can not be corrected this way .
such inefficiencies are due to finite detector resolution and , if uncorrected , produce significant artifacts in the correlations for heavy - ion collisions  @xcite .
the conventional correction method  @xcite involves removing particle pairs whose signals ( e.g. induced ionization , secondary particle showers , etc . )
fall within the resolution limits of the detector , and then removing pairs of particles from mixed - events which would have the same , relative locations in the detectors .
corrected results are obtained by increasing the minimum required separation between the detected signals of two nearby particles from zero until ratio @xmath19 stabilizes .
the practical consequence of this correction procedure is that the reference @xmath20 must be calculated by constructing uncorrelated pair counts in 2d bin @xmath12 by averaging over pairs of particles where each particle in a pair is selected from different collision events , referred to as _ mixed - events_.    the statistical error of the covariance in eq .
( [ eq1 ] ) equals the standard deviation of the distribution of covariance values corresponding to independent , statistically equivalent event samples ( event collections ) of underlying parent distributions for quantities @xmath17 , @xmath18 and @xmath11 .
analytical calculations of this error therefore represent @xmath17 , @xmath18 as random , event - wise fluctuating quantities relative to the parent distribution .
the sibling and mixed - event pair - numbers are similarly represented .    in ref .
@xcite reid and trainor derived a practical mixed - event method for calculating the reference in which random , event - wise fluctuations ( noise ) in @xmath17 and @xmath18 , which are common to both sibling and mixed - event pair numbers , cancel in the covariance , significantly reducing the errors . for large data volumes summing the total number of mixed - event pairs
can be computationally demanding .
various event - mixing algorithms have been developed by the heavy - ion community to reduce the necessary computation time while retaining sufficient statistical accuracy .
one such method was discussed in ref .
the choice of the reference and the event - mixing method strongly affects the statistical errors in the final correlation measurement . in this paper
the statistical noise reduction method of reid - trainor will be extended and applied to other , practical event - mixing algorithms .
the consequences of these event - mixing choices , or references , for the statistical uncertainties in the correlations will be quantified .
the present application is for ultra - relativistic heavy - ion collisions such as those measured by the star experiment  @xcite and by the experiments at the large hadron collider ( lhc )  @xcite .
the methods presented here are also directly applicable to correlation analysis of multi - particle production from any type of particle collision .
the present event - mixing technique for constructing an uncorrelated reference distribution has analogues in , for example , cosmology and acoustics .
measurements of the relative distance correlation between galaxies within an angular patch of the sky require an uncorrelated reference distribution .
the latter can be constructed from cross - correlated pairs of galaxies observed in different sky patches or from randomly generated distributions of galaxies  @xcite . in acoustical analysis of multiple , independent time series
the autocorrelation , or time - lag dependence , for each time series must be referenced to a cross correlation between two independent time series having the same lag time  @xcite . in these two examples
angular patches of sky or individual time series correspond to collision events and binned numbers of galaxies or acoustical amplitudes correspond to binned number of particles in the present analysis .
this paper is organized as follows . in sec .
[ sec : general ] the reid - trainor procedure is derived and extended . in sec .
[ sec : other ] their method is applied to other event - mixing algorithms .
monte carlo studies are discussed in sec .
[ sec : montecarlo ] .
conclusions are given in sec .
[ sec : conclusions ] .
in ref .  @xcite the event collection was separated into pairs of events with similar , total number of detected particles , or _
multiplicity_. event - mixing was only applied between the two events in each pair . for arbitrary
2d bin @xmath12 ( @xmath21 ) the total number of sibling pairs of particles in the collection is given by the sum _ mn & = & _ g=1^n_g _ j=1 ^ 2 ( |m+_j)(|n + _ j ) [ eq2 ] and the mixed - event pair sum is given by _ mn & = & _ g=1^n_g _ j^
> j=1 ^ 2 [ ( |m+_j)(|n + _ j^ ) + ( |m+_j^)(|n + _ j ) ] [ eq3 ] where indices @xmath22 and @xmath23 denote events and event groups , respectively . variables @xmath24 and @xmath25 are the parent distribution number of particles in bins @xmath26 ( bins are denoted with subscripts ) and fluctuations are represented with random variables @xmath27 and @xmath28 as in ref.@xcite .
event group index @xmath23 is suppressed in the notation for fluctuations @xmath27 and @xmath28 .
the number of pairs of events in the collection is @xmath29 where @xmath30 is the number of events in the collection .
note that averages of random variables @xmath27 and @xmath28 within an event collection only vanish in the @xmath31 limit . for large event numbers the summations in eqs .
( [ eq2 ] ) and ( [ eq3 ] ) are approximately given by _
mn & & n_events ( |m|n + ) + m_mn & & n_events ( |m|n ) , [ eq4 ] where @xmath32 is non - zero if the fluctuations in bins @xmath1 and @xmath2 are correlated . for relativistic heavy - ion collisions @xmath33
@xcite and for the purpose of calculating statistical errors the small contributions of @xmath34 can be neglected . the large event number limits are defined as @xmath35 .    from the above discussion the correlation quantity of interest is @xmath36 and we calculate the number of correlated pairs per reference pair given by  @xcite & = & - 1 _ mn - 1 .
[ eq5 ] the statistical error in @xmath37 equals the statistical error in @xmath38 , denoted by @xmath39 and is given by ( ) ^2 & = & ( ) ^2 + ( ) ^2 - [ eq6 ] where @xmath40 , etc .
are variances , @xmath41 is a covariance , and @xmath42 . simplifying eq .
( [ eq6 ] ) and averaging over event collections yields ( _ mn)^2 & = & |r_mn^2 + & = & + & = & [ eq7 ] where angle - brackets represent the average over independent event collections , @xmath43 ^ 2 $ ] is the variance of difference @xmath44 , @xmath45 , and @xmath46 .    the key result of ref .  @xcite was to show that event - pair - wise mixing eliminates contributions of single - particle fluctuations , leaving only those contributions from fluctuations in the number of pairs .
using eqs .
( [ eq2 ] ) and ( [ eq3 ] ) the variance in the numerator of eq .
( [ eq7 ] ) simplifies to ( s_mn - m_mn)^2 & = & ^2 + & = & _ j=1^n_events _ j^2 _ j^2 + _ j=1,odd^n_events ( _ j^2 _ j+1 ^ 2 + _ j+1 ^ 2 _ j^2 ) [ eq8 ] where averages over products of bin - wise fluctuations from different events vanish ( see appendix a )
. carrying out the above event averaging gives ( s_mn - m_mn)^2 & = & n_events _ ^2 _ ^2 + n_g _ ^2 _ ^2 + _ ^2 _ ^2 + & = & 2 n_events _ ^2 _ ^2 + & & 2 n_events |m|n 2 n_events |m|n , [ eq9 ] where @xmath47 and @xmath48 are the variances in bins @xmath1 and @xmath2 , respectively .
these variances result from the following averages : _ j=1^n_events _ j^2 _ j^2 & = & _ ^2 _ ^2 , + _ j=1^n_events _ j^2 _ j+1 ^ 2 & = & _ ^2 _ ^2 .
these factorized results are valid if the single - particle fluctuations in bins @xmath1 and @xmath2 are uncorrelated , consistent with the above assumption that correlation corrections to the statistical errors are negligible . in the last two steps in eq .
( [ eq9 ] ) the limit of poisson distributions is assumed , denoted here and in the remainder of this paper with a limit symbol ( arrow ) . in the last line in eq .
( [ eq9 ] ) event collection average @xmath49 is approximated with @xmath50 from a single event - collection .
correlated single - particle fluctuations in the event - wise occupancies of different bins ( other than that caused by collision dynamics ) could also be significant if the range of multiplicities in the event collection is large or if the bin size is large relative to the detector acceptance .
in such cases the event - wise frequency distributions of @xmath51 and @xmath52 would be inconsistent with a poisson distribution . in sec .
[ sec : montecarlo ] the frequency distributions for realistic bin sizes are shown to be accurately represented with poisson distributions . finally , the statistical error in @xmath38 for the event - mixing algorithm in ref .
@xcite is & = & = .
[ eq10 ] note that this error equals @xmath53 , or the square - root of the sum of the relative variances of the sibling and mixed - event pair numbers and is not determined by the errors in the bin - wise single - particle number  @xcite .    for 2d diagonal bins ( @xmath54 )
only unique particle pairs are counted ( no self - pairs ) which requires a separate derivation . in this case
@xmath55 and @xmath56 are given by _
mm & = & _ j=1,odd^n_events [ eq11 ] + m_mm & = & _ j=1,odd^n_events ( |m+_j ) ( |m+_j+1 ) [ eq12 ] and as before the large event number limits are |s_mm & = & n_events ( |m^2 - |m + _ ^2 ) n_events |m^2 [ eq13 ] +    variance quantity yields ( s_mm - m_mm)^2 & = & ^2 + & = & _ j=1,odd^n_events ( _ j^2 + _ j+1 ^ 2 -2 _ j _ j+1 - 2|m - _ j - _
j+1)^2 + & = & _ j=1^n_events ( _ j^4 + _ j^2 + 2|m^2 - 2_j^3 - 4|m_j^2 ) + & & + 6 _ j=1,odd^n_events _ j^2 _ j+1 ^ 2 [ eq15 ] where in going from the first to the second line in eq .
( [ eq15 ] ) cross terms between different events cancel in the poisson limit and in the last line we use the fact that event - collection averages of terms linear in the fluctuation vanish . inserting the moments of poisson distributions ( see appendix b ) , the resulting variance is ( s_mm - m_mm)^2 & & n_events |m^2 [ eq16 ] and the statistical error in @xmath57 , assuming poisson distributions , is & = & = .
[ eq17 ] the multiplicative factor in the numerator for this diagonal bin error is @xmath58 larger than that for off - diagonal bins given in eq .
( [ eq10 ] ) due to the reduced number of unique sibling and mixed - event pairs in diagonal bins .
in general the event collection could be sub - divided into an arbitrary number of groups of events with @xmath59 events per group where @xmath60 and with @xmath61 combinations of mixing event - pairs within each group .
for this general case the sibling and mixed - event particle - pair sums and their limits are given by _
mn & = & _ g=1^n_g _ j=1^ ( |m+_j)(|n+_j ) +    m_mn & = & _ g=1^n_g _ j^ >
j=1^ [ ( |m+_j)(|n+_j^ ) + ( |m+_j^)(|n+_j ) ] +    there are more mixed - event particle pairs than sibling and the ratio @xmath38 must be normalized such that the correlation signal , @xmath44 , equals zero in the absence of correlations . the ratio @xmath38 for the general case
is therefore defined as _ mn & = & = ( - 1 ) [ eq20 ] and _ mn - 1 & = & .
[ eq21 ] the variance is evaluated as before where ^2 & = & \ { _ g=1^n_g [ ( - 1 ) _ j=1^ _ j _ j - _ j^ > j = 1^ ( _ j _
j^ + _ j^ _ j ) ] } ^2 + & = & _ g=1^n_g [ ( - 1 ) _ j=1^
_ j _ j - _
j^ > j = 1^ ( _ j _
j^ + _ j^ _ j ) ] ^2 + & = & _ g=1^n_g [ ( - 1 ) _ j=1^ _
j _ j - _ j^ j _ j _ j^ ] ^2 [ eq22 ] where cross terms among different event groups vanish .
this last expression can be further simplified using the fact that the event - collection average of @xmath62 for different events @xmath63 vanishes , resulting in ^2 & = & _ g=1^n_g [ ( - 1)^2 _ j=1^ _ j^2 _ j^2 + _ j^ j=1^ _ j^2 _ j^^2 ] + & & = _ g=1^n_g [ ( ^2 - 2 + 1 -1 ) _ j=1^ _
j^2 _ j^2 + _ j^ , j=1^ _ j^2 _ j^^2 ] + & & = n_g _ g=1^n_g _ j=1^ _ j^2
_ j^2 + _ g=1^n_g _ j=1^ _ j^2 _ j^=1^ _
[ eq23 ] the first term is an average over all events in the collection and results in the product of variances @xmath64 in the absence of correlated fluctuations between different bins .
the second term is the average over all groups of the product of variances for each group .
( [ eq23 ] ) therefore reduces to ^2 & = & n_g ( ^2 - 2 ) _ ^2 _ ^2 + _ g=1^n_g ( _ ^2 + _ g ) ( _ ^2 + _ g ) + & & = n_g ( - 1 ) _ ^2 _ ^2 ( - 1 ) n_events ( |m|n ) , [ eq24 ] where the variances within each group are expressed as the average variance for all event collections plus random fluctuations @xmath65 and @xmath66 , where averages @xmath67 , @xmath68 , and @xmath69 equal zero .
the error in an off - diagonal bin is therefore & = & = .
[ eq25 ] this error equals that in eq .
( [ eq10 ] ) when @xmath70 and approaches @xmath71 as the number of events in each mixing group increases and approaches the entire number of events in the collection .
again , note that this error equals @xmath53 .
the improvement in statistical accuracy achieved by increasing @xmath59 can be weighed against the increased computational requirements in order to optimize the number of event groups to be used in the analysis of data .
in typical analyses of heavy - ion collision data  @xcite correlations are calculated by reading an unordered list of event data and sorting the particle kinematic information from those events into subsets corresponding to a global event - property such as overall multiplicity , a proxy for the overlap , or _ centrality _ , between the colliding nuclei .
for each centrality all sibling particle - pairs are processed and to reduce computational requirements mixed - event particle - pairs are only processed for the first and second events in the collection , then for the second and third events , and so on until all events in the centrality bin have been processed .
this event - mixing algorithm can be extended such that each event is mixed with the next two events in the list , or the next three events , etc . in this section statistical errors
are derived for the above single , double and multiple - event mixing algorithms .
it will be shown that the noise reduction obtained in ref .
@xcite also occurs for these other types of event - mixing methods .      in this event mixing algorithm particle pairs are constructed for all events by mixing particles in event @xmath22 with those in event @xmath72 . the number of sibling and mixed - event particle pairs and their limits are given by _
mn & = & _ j=1^n_events ( |m+_j)(|n+_j ) _ j=1^n_events s_mn^j +    m_mn & = & _ j=1^n_events - 1 [ ( |m+_j)(|n+_j+1 ) + ( |m+_j+1)(|n+_j ) ] + & & _ j=1^n_events - 1 ( m_mn^j , j+1 + m_mn^j+1,j ) +    definitions @xmath73 and @xmath74 are introduced for brevity .
the approximation @xmath75 occurs because the last event in the event list has no other event to mix with . for @xmath76 those missing pairs
can be ignored in the above summations .    in this single - event mixing algorithm
there are approximately twice as many mixed - event pairs as sibling pairs .
therefore the variance to be calculated is given by ( 2s_mn - m_mn)^2 & = & ^2 .
[ eq28 ] twice the sum over sibling pairs can be written as two , identical summations .
the resulting sibling and mixed - event sums can be separated into sums over even and odd values of event number @xmath22 . as a result eq .
( [ eq28 ] ) can be rearranged into sums with the identical structure as that in ref .
@xcite , thus achieving the same level of noise reduction for this alternate event - mixing algorithm .
the result is ( 2s_mn - m_mn)^2 & = & ^2 .
[ eq29 ] from ref .
@xcite we find that _
mn^j + s_mn^j+1 -m_mn^j , j+1 - m_mn^j+1,j & = & _ j _ j + _ j+1 _
j+1 - _ j _ j+1 - _ j+1 _ j .
+ [ eq30 ] substituting the above result into eq .
( [ eq29 ] ) and combining the even and odd numbered summations gives ( 2s_mn - m_mn)^2 & = & ^2 + & = & _ j=1^n_events ( 2 _ j _ j - _ j _ j+1 - _ j+1 _ j)^2 [ eq31 ] where cross terms for different events vanish as before .
continuing , we find that the above equation in the poisson limit reduces to _ j=1^n_events ( 4 _ j^2 _ j^2 + _ j^2 _ j+1 ^ 2 + _ j+1 ^ 2 _ j^2 ) & = & 6 n_events _ ^2 _ ^2 + & & 6 n_events ( |m|n ) . [ eq32 ]    the statistical error for off - diagonal bins is & = & = .
[ eq33 ] this error is a factor of @xmath77 smaller than that for the reid - trainor method due to the larger number of mixed - event pairs ( twice as many pairs ) generated by this algorithm .
it is important to note that the noise reduction achieved via the mixing algorithm in ref .
@xcite is also achieved with the single - event mixing method and the error is determined by pair number where the error in eq .
( [ eq33 ] ) equals @xmath53 .    for diagonal bins _ mm & = & _ j=1^n_events ( |m+_j)(|m+_j - 1 ) + & = & _ j=1^n_events ( s_mm^j -|m -_j ) +    m_mm & = & _ j=1^n_events - 1 ( |m+_j)(|m+_j+1 ) = _ j=1^n_events - 1 m_mm^j ,
j+1 +    particle pairs are included in the summations .
the variance normalization factor , @xmath78 , for the diagonal bins is also equal to 2 .
the variance quantity is then given by ( 2s_mm - m_mm)^2 & = & ^2 + & & = ^2 [ eq36 ] where @xmath79 .
separating the above into summations over even and odd event numbers , doubling the sibling - pair summation , assuming that @xmath76 ( in order that the contribution of missing mixed - events at the end of the event list is negligible ) , and using the noise reduction result from ref .
@xcite result in ( 2s_mm - m_mm)^2 & = & ^2 + & & = _ j=1^n_events ( _ j^4 + _ j^2 _ j+1 ^ 2 + |m^2 + _ j^2 - 2|m_j^2 - 2_j^3 ) + & & + _ j j^ ( _ j^2 _ j^^2 - |m_j^2 -|m_j^^2 + |m^2 ) . [ eq37 ]
the second term vanishes in the @xmath76 limit where n_events^2 ( _ ^4 - |m _ ^2 - |m _ ^2 + |m^2 ) & & n_events^2 ( |m^2 - |m^2 -    & = & 0 .
[ eq38 ] the first term in eq .
( [ eq37 ] ) , using the moments of a poisson distribution in appendix b , simplifies to ( 2s_mm - m_mm)^2 & = & + & & n_events _ j=1^n_events ( _ j^4 + _ j^2 _ j+1 ^ 2 + |m^2 + _ j^2 - 2|m_j^2 - 2_j^3 ) + & & n_events ( 3|m^2 + |m ) + |m^2 + |m^2 + |m -2|m^2 - 2|m = 3n_events    [ eq39 ]    the statistical error in diagonal bins for the single - event mixing algorithm is & = & = [ eq40 ] where the multiplicative factor in the numerator is @xmath58 larger than that for the off - diagonal errors in eq .
( [ eq33 ] ) .      in this
method particles in each event @xmath22 are mixed with particles in events @xmath72 and @xmath80 .
the sibling and mixed - event pair summations , using the symbols @xmath73 and @xmath74 introduced above , and their limits are given by _
mn & = & _ j=1^n_events s_mn^j +    m_mn & = & _ j=1^n_events - 1 ( m_mn^j , j+1 + m_mn^j+1,j ) + _ j=1^n_events - 2 ( m_mn^j , j+2 + m_mn^j+2,j ) +    & & 4n_events ( |m|n ) . [ eq42 ] the normalization factor for @xmath38 is approximately 4 and the variance to be calculated is ( 4s_mn - m_mn)^2 & = & ^2 .
[ eq43 ]    by writing out the sum over sibling pairs four times it is found that the resulting summations in eq .  ( [ eq43 ] ) can be separated into sums over pairs of events where arbitrary event @xmath22 is mixed with event @xmath81 , and in this example @xmath82 and 2 . for further brevity
the symbol @xmath83 is introduced which is defined by _
mn^j , k & & s_mn^j + s_mn^j+k - m_mn^j , j+k - m_mn^j+k , j .
[ eq43b ] after rearranging the sums the above variance in the @xmath76 limit can be expressed as ( 4s_mn - m_mn)^2 & = & ^2 [ eq44 ] where , as in ref .
mn^j , k & = & _ j _ j + _
j+k - _ j _ j+k - _ j+k _ j .
[ eq45 ] substituting @xmath83 into eq .
( [ eq44 ] ) and recombining the summations gives ( 4s_mn - m_mn)^2 & = & ^2 .
[ eq46 ] eliminating cross terms involving products of fluctuations between different events results in ( 4s_mn - m_mn)^2 & = & _ j=1^n_events ( 16 _ j^2 _
j^2 + _ j^2 _ j+1 ^ 2 + _ j+1 ^ 2 _ j^2 + _ j^2 _ j+2 ^ 2 + _ j+2 ^ 2 _ j^2 ) + & = & 20 n_events _ ^2 _ ^2 20 n_events ( |m|n ) .
[ eq47 ]    the statistical error in off - diagonal bins for the double - event mixing algorithm is given by & = & = [ eq48 ] which is smaller than the above , single - event mixing error and is @xmath84 smaller than the error resulting from the reid - trainor method .
again , this error equals @xmath53 .
the error for diagonal bins is obtained as in the above sections and is given by & = & .
[ eq49 ]      the preceding event - mixing methods can be readily extended to an arbitrary number of mixing events @xmath85 . in the following derivation
however we continue to assume that @xmath86 .
the sums of sibling and mixed - event particle pairs in off - diagonal bins and their limits are _ mn & = & _ j=1^n_events s_mn^j +    m_mn & = & _
j=1^n_events -1(m_mn^j , j+1 + m_mn^j+1,j ) + _ j=1^n_events -2(m_mn^j , j+2 + m_mn^j+2,j ) + & + & + _ j=1^n_events -k(m_mn^j , j+k + m_mn^j+k , j ) +    eq .  ( [ eq43b ] )
the variance is given by ( 2ks_mn - m_mn)^2 & = & ^2 [ eq52 ] where event index @xmath87 , @xmath88 .
recombining the summations over sibling pair quantities @xmath89 results in @xmath90 .
the summations over the last two terms , @xmath91 , include each mixed - event only once and may be simplified to @xmath92 .
the variance can then be expressed as ( 2ks_mn - m_mn)^2 & = & ^2 .
[ eq53 ] as before the average over cross terms between different events vanishes resulting in ( 2ks_mn - m_mn)^2 & = & _ j=1^n_events ( 4k^2 _ j^2 _ j^2 + _ k=1^k ( _ j^2 _ j+k^2 + _ j+k^2 _
j^2 ) ) + & & = n_events ( 4k^2 + 2k)_^2 _ ^2
( 4k^2 + 2k ) n_events(|m|n ) [ eq54 ] and the error in off - diagonal bins is & = & = .
[ eq55 ] evaluating this result for @xmath93 and 2 gives the above errors for single - event and double - event mixing in eqs .
( [ eq33 ] ) and ( [ eq48 ] ) , respectively .
the above error equals @xmath53 .    to summarize , for each event - mixing algorithm presented here the sibling minus mixed - event pair summations
can be accurately rearranged into sums over event - pair - wise differences between sibling and mixed - event pairs as in ref .
@xcite , resulting in significant statistical noise reduction .
the statistical errors for the correlation quantity @xmath94 are determined by the total numbers of sibling and mixed - event pairs of particles in 2d bin @xmath12 , given by @xmath53 .      the correlation quantity @xmath94 in eq .
( [ eq5 ] ) represents the event - averaged number of correlated particle pairs per final - state particle pair produced in the collisions . in some analyses  @xcite
the authors reported the number of correlated particle pairs per final - state particle in which case @xmath94 was scaled by a prefactor  @xcite proportional to total event - multiplicity .
for example , in the analysis of au + au collisions at @xmath95 = 200  gev  @xcite the prefactor was computed using efficiency corrected , charged - particle yields within the detector acceptance defined by pf_ci^ ( , ) & = & ^1/2 = [ eq56 ] for correlations on 2d angular space @xmath96 for arbitrary particle pair ( 1,2 ) and at mid - rapidity where @xmath97 is approximately constant .
subscript `` ci '' ( charge independent ) means that all charge - sign pairs are included , i.e. @xmath98 , @xmath99 , @xmath100 , and @xmath101 .
other kinematic projections are possible . in ref .
@xcite correlations on 2d transverse momentum space were presented with corresponding prefactor pf_ci^(p_t , p_t ) & = & ^1/2 .
[ eq57 ] statistical errors in the prefactor could be included in the final , reported correlation quantities but this factor is usually treated as a simple scaling of the correlation measurement .
insight into the dynamics underlying the observed correlations can be gained by making further selections on the particle pairs .
for example , correlations for positive versus negative charged particle pairs , or like - charge - sign pairs ( @xmath98 and @xmath99 ) versus unlike - sign pairs ( @xmath100 and @xmath101 ) , can help differentiate charge - ordering effects  @xcite in hadronization , resonance decay contributions , quantum correlations for identical particles  @xcite , and photon pair - production contamination .
selection of low-@xmath4 or high-@xmath4 particles can help differentiate soft and hard - scattering processes  @xcite .
similarly , for correlations on 2d transverse momentum - space particle pairs with smaller or larger relative azimuthal angles , e.g. @xmath102 ( near - side ) or @xmath103 ( away - side ) , in conjunction with charge - pair selections , can help differentiate the @xmath4 structures of near - side , jet - like fragmentation from that of away - side , or back - to - back dijet fragmentation .    if the prefactor , which is the square - root of a pair density , is computed using the same charge - sign and/or kinematic selections used in constructing the correlation @xmath94 , then the statistical error in the final quantity , @xmath104 , is unaffected by those selections .
this is a result of the noise reduction inherent in the event - mixing methods presented here and in ref .
@xcite where the statistical error is inversely proportional to the square - root of the number of particle pairs in the bin , @xmath105 .
for example , if only like - sign pairs are selected the total number of pairs is reduced by 1/2 if the numbers of positive and negative particles are the same ( approximately true for relativistic heavy - ion collisions ) .
the prefactor is therefore reduced by @xmath106 while the error in each bin is increased by @xmath58 .
similarly , selecting only near - side or away - side particle pairs reduces the number of pairs by 1/2 which again decreases the prefactor by @xmath106 and increases the bin errors by @xmath58 . in general , such pair selections do not affect the prefactor scaled statistical errors studied here .
this result is summarized as follows : & = & pf_ci [ eq58 ] where the left - hand - side represents the root - mean - square statistical error of quantity @xmath104 for arbitrary charge - sign and/or azimuth angle pair selections .
( [ eq58 ] ) applies to each event - mixing method discussed here .
throughout the preceding derivations it was assumed in the final steps that the frequency distributions for event - wise particle number in bins @xmath1 and @xmath2 were consistent with poisson distributions , i.e. variance @xmath47 equaled mean @xmath24 . in practical data analysis , collision events are usually grouped into subsets based on global event properties such as total multiplicity in the detector acceptance  @xcite .
clearly the multiplicity range chosen for the event collections can affect the frequency distributions in the @xmath1 and @xmath2 bins .
in addition , if the @xmath1,@xmath2 bin sizes are comparable to the acceptance , then the multiplicity distribution of the event collection will also affect that in the bins . in this
section monte carlo simulations are used to determine the differences between actual event - wise particle number frequency distributions in the bins and a poisson distribution .
if the distributions differ significantly from a poisson , then the actual variances should be used in the above derivations .    in au + au collisions at energies of 200  gev per colliding nucleon + nucleon
pair over a thousand charged particles can be produced at mid - rapidity within the tracking acceptance of the star detector  @xcite ( @xmath107  gev/@xmath108 , @xmath109 and @xmath110 in azimuth ) . for events selected with a minimum - bias trigger  @xcite
the multiplicity frequency distribution for total , observed multiplicity , @xmath111 , is approximately proportional to @xmath112 except near the lower and upper multiplicity end - points  @xcite .
the events are then subdivided into multiplicity classes or centrality .
for example , the multiplicity ranges for peripheral , intermediate ( mid - central ) , and head - on ( most - central ) collisions from ref .
@xcite are [ 15,35 ] , [ 152,187 ] and [ 952,1002 ] , respectively .    for the present example , correlations on 2d transverse momentum space were considered where the kinematic variable used is transverse rapidity , defined by @xmath113 $ ] , @xmath114 and for non - identified charged particles @xmath115 , which regulates the singularity in @xmath116 at @xmath117 , was assumed to be the pion mass , 0.14  gev/@xmath118 .
the transverse rapidity spectrum for each multiplicity bin was obtained from the efficiency and contamination corrected @xmath4 spectrum data  @xcite and is accurately represented by a levy distribution given by & = & .
[ eq59 ] parameters @xmath119 and @xmath2 from data fitting for peripheral , mid - central and most - central collisions are ( 0.154  gev , 10.42 ) , ( 0.199  gev , 13.32 ) , and ( 0.226  gev , 17.55 ) , respectively . amplitude @xmath120 is determined by event - wise multiplicity . in the correlation analysis
particle pairs were binned on 2d @xmath121 space where @xmath122 $ ] with @xmath123 2d @xmath12 bins of equal size .    in the simulations , the power - law distribution @xmath112 was randomly sampled within each of the above multiplicity ranges , the single - particle distribution in eq .
( [ eq59 ] ) was sampled for the event - wise multiplicity , the 25 single - particle @xmath116 bins were filled , and the number of particles in each @xmath116 bin for each event was recorded .
the particle number frequency distributions for the 25 @xmath116 bins for the simulated events were compared with poisson distributions .
this was done by calculating the ratio of the computed variance , @xmath124 , to the variance for a poisson distribution , @xmath125 , where @xmath126 is the mean particle number in the @xmath116 bin .
the resulting ratios were compared with unity .
the uncertainty in the ratio of variances for each @xmath116 bin was estimated by repeating the above simulation procedure many times , corresponding to statistically independent collections of events .
the mean and width of the resulting distribution of variance ratios were computed .
results for @xmath127 for six of the twenty - five @xmath116 bins are listed in table  [ table1 ] and shown in fig .
[ figure1 ] for the three event - multiplicity ranges listed above .
for the peripheral , mid - central and most - central collisions @xmath128 , @xmath129 and @xmath129 events were generated , respectively .
one - hundred event collections were simulated for each centrality .
the calculations provided sufficient statistical accuracy to determine a significant difference between @xmath127 and unity for most @xmath116 bins . in all cases @xmath130 agreed with the poisson limit @xmath131 to within 3% or less , thus confirming the poisson assumption for the present application to au + au collisions at 200  gev .
[ cols="^,^,^,^ " , ]      ( color online ) ratios @xmath127 as described in the text and listed in table  [ table1 ] for three event - multiplicity ranges where the solid , dashed - dotted and dashed lines connecting the points correspond to peripheral , mid - central and most - central collisions , respectively.,width=384 ]    monte carlo simulations were also used to verify the noise reduction predicted in ref .
@xcite and for the derivations in sections  [ sec : general ] and [ sec : other ] .
the simulations were also used to check the error scaling factors @xmath58 , @xmath132 and @xmath133 for the reid - trainor , single- and double - event mixing algorithms . in ref .
@xcite and in appendix c the statistical error for correlations obtained with mixed - event particle pairs taken from a different set of events than those used for the sibling particle pairs is shown to be proportional to @xmath134^{1/2}$ ] .
this error depends on particle number and therefore differs significantly from the errors derived in ref.@xcite and in sections  [ sec : general ] and [ sec : other ] which are proportional to @xmath135 .
finally , the simulations were used to verify that the errors in 2d diagonal bins are larger than those in off - diagonal bins by a factor of @xmath58 .    as in the preceding simulation
the @xmath112 frequency distribution was randomly sampled to obtain an event multiplicity , the @xmath116 distribution in eq .
( [ eq59 ] ) was sampled , and the 25 single - particle @xmath116 bins were filled .
sibling and mixed - event pairs were counted for each 2d @xmath12 bin on @xmath121 space using the event - mixing algorithms in ref .
@xcite ( see sec .  [
sec : general ] ) , the single - event mixing described in sec .
[ sec : single ] , or the double - event mixing described in sec .
[ sec : double ] .
the correlation quantity @xmath94 was computed for all @xmath12 bins using all simulated events in the collection .
this entire process was repeated many times thus generating independent , statistically equivalent event collections and corresponding values for @xmath94 . from these distributions of @xmath94 the means and widths in each @xmath12 bin were obtained , including those for on - diagonal and off - diagonal bins .
the resulting monte carlo uncertainties are directly comparable to the analytic results for statistical errors in eqs .
( [ eq10 ] ) , ( [ eq17 ] ) , ( [ eq33 ] ) , ( [ eq40 ] ) , ( [ eq48 ] ) , and ( [ eq49 ] ) .
results are reported below for the mid - central multiplicity bin with @xmath136 $ ] and for two diagonal @xmath137 bins on @xmath121 defined by ranges ( [ 1.56,1.70],[1.56,1.70 ] ) and ( [ 2.96,3.10],[2.96,3.10 ] ) , and for the corresponding off - diagonal @xmath12 bin defined by ( [ 1.56,1.70],[2.96,3.10 ] ) .
the mean , single - particle numbers in these two @xmath116 bins are 14.3 and 3.39 , respectively .    the scaling of the error with particle pair number , rather than with single particle number can be checked by comparing the simulated errors for two diagonal bins . for the above diagonal bins the ratios of the analytical errors in the lower to upper @xmath116 bins are 0.237 for particle - pair scaling ( i.e. with noise reduction ) and 0.462 for single - particle scaling
( statistically independent events for sibling and mixed - event particle - pairs , see appendix c ) .
the simulated ratios for the reid - trainor , single - event mixing , and double event - mixing algorithms are 0.241 , 0.226 and 0.236 , respectively , which agree with the noise reduced , particle - pair scaling
.    the scaling factors for the errors can be checked by comparing the analytic ratio for single - event mixing to reid - trainor mixing , equal to @xmath138 , and the ratio for double - event mixing to reid - trainor mixing , equal to @xmath139 , to the corresponding ratios from the simulations .
the monte carlo error results for the lower and upper @xmath116 diagonal bins are 0.851 and 0.911 ( ratio of single - event mixing to reid - trainor mixing ) and 0.785 and 0.803 ( ratio of double - event mixing to reid - trainor mixing ) , respectively .
both pairs of monte carlo results are in good agreement with the corresponding analytic expressions .    finally , the @xmath58 increase in the diagonal bin errors can be checked via the ratio @xmath140 where @xmath141 is the error in bin @xmath12 .
the analytical value of this ratio equals 1/2 for the three event - mixing algorithms considered in this section .
the simulation results for the reid - trainor , single - event mixing , and double - event mixing algorithms using the above bins are 0.52 , 0.52 and 0.50 , respectively , in good agreement with the analytic errors .
the simulated and analytic absolute errors agree within 2%-6% , where the few percent effects of non - poisson distributions discussed above contribute here as well .
the scaling factors discussed above and the differences between the errors obtained with particle - pair number scaling versus single - particle number scaling are much larger than these few percent differences .
overall , the simulations confirm the analytic expressions .
the poisson distribution assumption should be checked for each application .
the event - mixing and error propagation method of reid - trainor  @xcite for a multiplicity - ordered event collection , grouped into pairs of events for particle - pair mixing , was generalized to include event - mixing groups with arbitrary number of events .
the method was then applied to other event - mixing algorithms which have been used in correlation analysis of relativistic heavy - ion collision data .
the noise reduction arising from the partial cancellation of event - wise single - particle multiplicity fluctuations was shown to occur for those other methods in which particle - pair mixing is done between each event in the collection and the next 1 , 2 or @xmath142 ( where @xmath143 ) events in the list .    the statistical errors for analyses which use the event - mixing group method of ref .
@xcite are given in eqs .
( [ eq10 ] ) and ( [ eq17 ] ) for two events per mixing group and in eq .
( [ eq25 ] ) for @xmath144 events in each event - mixing group . errors for correlation analysis based on single - event mixing and double - event mixing are given in eqs .
( [ eq33 ] ) , ( [ eq40 ] ) and eqs .
( [ eq48 ] ) , ( [ eq49 ] ) , respectively .
correlation errors for all event - mixing algorithms in which each event is mixed with the next @xmath142 events in the list are given in eq .
( [ eq55 ] ) .
statistical errors for correlations reported as the number of correlated pairs per final - state particle , including charge and/or kinematic dependent pair selections  @xcite , are given in eq .
( [ eq58 ] ) .
the statistical errors in correlations for each of the event - mixing algorithms studied here using the reid - trainor method are inversely proportional to the square - root of the number of pairs of particles in a bin .
this result is in contrast to that obtained when the two - particle reference distribution is constructed from a set of events independent from those used for sibling pairs . in the latter case
the errors are inversely proportional to the square - root of the number of particles in the bin and can be much larger than the errors which follow from the event - mixing methods discussed here and in ref .  @xcite .
as long as the sum over sibling particle pairs minus mixed - event particle pairs can be accurately approximated by sibling minus mixed - event particle pair differences summed over pairs of events as in ref .
@xcite and in this paper , then the reid - trainor error reduction will follow and the correlation errors will be inversely proportional to the square - root of the total , bin - wise number of particle pairs .
it should also be noted that error propagation algorithms based on bin - wise single - particle counts which neglect covariance @xmath41 in eq .
( [ eq6 ] ) will produce erroneous results .
monte carlo simulations were employed to test the accuracy of the poisson distribution approximations used in the error calculations .
simulations were also used to verify the predicted reduction in statistical errors for the reid - trainor , single - event and double - event mixing algorithms .
finally , the analytic error results presented here may be helpful in optimizing correlation data analysis when more aggressive event - mixing algorithms are considered and the benefit from reduced statistical errors needs to be weighed against the additional computational cost .    the authors would like to thank professor tom trainor of the univ .
of washington for many informative discussions relevant to this work .
this research was supported in part by the office of science of the u. s. department of energy under grants de - fg02 - 94er40845 and de - sc0013391 .
in deriving eq .  ( [ eq8 ] ) event averages and cross terms involving averages over pairs of different events must be calculated .
the necessary steps are given here .
similar calculations are required for the other event - mixing algorithms described in this paper . starting with the first line in eq .
( [ eq8 ] ) , replacing the group summation in the first part of eq .
( [ eq8 ] ) with event summations , expanding the quadratic term , and retaining the non - vanishing averages yield the following : ( s_mn - m_mn)^2 & = & ^2 + & = & ^2 + & = & _ j=1,odd^n_events ( _ j _ j - _ j _
j+1 ) ^2 + _ j=2,even^n_events ( _ j _ j - _ j _
j-1 ) ^2 + & + & _ j j^,odd ( _ j _ j - _ j _
j+1 ) ( _ j^ _ j^ - _ j^ _ j^+1 ) + & + & _ j j^,even ( _ j _ j - _ j _
j-1 ) ( _ j^ _ j^ - _ j^ _ j^-1 ) + & + & 2 _ j=1,odd^n_events _ j^=2,even^n_events ( _ j _ j - _ j _
j+1 ) ( _ j^ _ j^ - _ j^ _ j^-1 ) + & = & _ j=1^n_events _ j^2 _ j^2 + _ j=1,odd^n_events _ j^2 _ j+1 ^ 2 +
_ j=2,even^n_events _ j^2 _ j-1 ^ 2 .
the last line is equivalent to eq .
( [ eq8 ] ) .
event collection averages over cross terms involving pairs of different events , for example & & _ j=1,odd^n_events _ j^2 _ j _ j+1 ,
_ j j^,odd _ j _ j _ j^ _ j^ , _ j
j^,odd _ j _ j _ j^
_ j^+1 + & & _ j=1,odd^n_events _ j^=2,even^n_events _ j _ j _ j^ _ j^ , vanish because bin - wise multiplicity fluctuations are uncorrelated from event - to - event .
cross term _
j=1,odd^n_events _ j^=2,even^n_events _ j _ j _ j^ _ j^-1 & & + & & = _ j=1,odd^n_events _ ^n_events _ j _ j _ j^ _ j^-1 + _ j=1,odd^n_events _ j _ j+1 _ j^2 also vanishes for the same reason where the @xmath145 term is separated out .
moments of the multiplicity frequency distribution in arbitrary bin @xmath1 involve averages of the event - wise fluctuations _
j=1^n_events _ j^p where exponent @xmath146 and @xmath147 , where @xmath148 is the multiplicity in bin @xmath1 for event @xmath22 . this average can be expressed as _
j=1^n_events ( m_j - |m)^p & = & _ m _
j=1^n(m ) [ ( m_j - |m)^p]|_m_j = m + & = & _ m ( m-|m)^p where @xmath149 is the number of events with @xmath150 particles in bin @xmath1 and @xmath151 sums over all events with @xmath152 . for poisson distributions @xmath153 where @xmath154 .
the moments are calculated via the summations given by _ j=1^n_events ( m_j -    the resulting poisson weighted sums can be related to bell polynomials  @xcite , defined by b_p(|m ) & = & _ m=0^ m^p p(m,|m ) where b_0(|m ) & = & 1 + b_1(|m ) & = & |m + b_2(|m ) & = & |m^2 + |m + b_3(|m ) & = & |m^3 + 3|m^2 + |m + b_4(|m ) & = & |m^4 + 6|m^3 + 7|m^2 + |m and so on .
the first few values of @xmath155 for @xmath156 $ ] are 1 , 0 , @xmath24 , @xmath24 , @xmath157 , respectively .
the statistical error for correlation quantity @xmath94 when independent sets of events are used to calculate the sibling and mixed - event pair sums is derived here for the event - mixing algorithm in ref .  @xcite . as before _
mn & = & _ j=1^n_events ( |m+_j)(|n+_j ) +      @xmath158 and @xmath159 are statistically equivalent .
the error of @xmath94 squared equals @xmath160 , however , the fluctuations in single - particle number no longer cancel as in ref .
@xcite and eq .
( [ eq8 ] ) . substituting the above summations for @xmath55 and @xmath56 into @xmath161 and noting that all cross terms between different events vanish ,
result in ( s_mn - m_mn)^2 & = & _ j=1^n_events ( |m^2_j^2 + |n^2_j^2 + _ j^2 _ j^2 ) + & & + _ j^=1^n_events ( |m^2_j^^2 + |n^2_j^^2 ) + _ j^=1,odd^n_events ( _ j^^2 _ j^+1 ^ 2 + _ j^+1 ^ 2 _ j^^2 ) + & & = 2n_events |m^2 _ ^2 + |n^2 _ ^2 + _ ^2 _ ^2 2n_events ( |m^2|n + |n^2|m + |m|n ) and the error in @xmath38 is & = & ^1/2 . comparing this error to the corresponding statistical error when the same event collection is used for the sibling and mixed - event pairs in eq .
( [ eq10 ] ) yields the ratio & = & . for 200  gev au + au collisions with events selected in the mid - central range
@xmath136 $ ] and for the off - diagonal @xmath12 bin on @xmath121 , defined by the ranges ( [ 1.56,1.70],[2.96,3.10 ] ) where @xmath24 = 14.3 and @xmath25 = 3.39 , the above ratio of statistical errors is 0.23 .
the noise reduction of the reid - trainor method is more than a factor of 4 in this instance .
m. s. daugherity , `` two - particle correlations in ultra relativistic heavy ion collisions , '' ph.d .
thesis , the university of texas at austin , ( 2008 ) , https://drupal.star.bnl.gov/star/files/daugherity_dissertation.pdf .
e. w. oldag , `` two - particle correlations in angular and momentum space in relativistic heavy ion collisions at star , '' ph.d .
thesis , the university of texas at austin , ( 2013 ) , https://drupal.star.bnl.gov/star/files/oldag_dissertation_20132.pdf .