observables in relativistic heavy ion experiments are often reported as functions of measured total inelastic cross section fraction ( centrality ) and are related to the initial collision geometry using monte carlo glauber ( mcg ) or optical glauber models  @xcite .
an example is charged particle multiplicity @xmath4 versus number of participating nucleons @xmath0 .
it is therefore imperative that the level of accuracy of the glauber models be understood and that relevant experimental information be used to estimate collision centrality .
for instance , previous studies  @xcite concluded that the uncertainties in relating peripheral collision data to the corresponding initial collision geometry were very large , thus contributing to the omission of those data from publications  @xcite . in another example , studies of high-@xmath5 suppression of charged particle production for peripheral to mid - central collisions are presently limited by large uncertainties in the estimated number of binary collisions which is used as a scaling variable  @xcite .    in this work
we re - examine the accuracy of monte carlo glauber model descriptions of minimum - bias multiplicity frequency distribution data from rhic using updated density and multipliciity production model parameters and a sensitive power - law inspired representation introduced by trainor and prindle ( tp )  @xcite . after demonstrating the accuracy of our monte
carlo glauber model in describing rhic data we then use the power - law analysis method of ref .
@xcite to estimate the uncertainties in the mapping relationships between initial stage collision geometry quantities and centrality for heavy ion collision systems relevant to the rhic program .
we find that previous uncertainties  @xcite in centrality measures are too pessimistic .
the centrality analysis method of trainor and prindle  @xcite exploits the approximate power - law dependence of the multiplicity frequency distribution data @xmath6 ( @xmath7 is the number of triggered events ) for minimum - bias   collisions and uses proton - proton ( p - p ) multiplicity production data to constrain the peripheral collision end - point of the @xmath6 distribution .
their analysis demonstrates that uncertainties in centrality determination and mcg geometry measures can be significantly reduced , particularly in the peripheral region .    in this
paper mcg results and errors for centrality bin average quantities @xmath8 ( mean impact parameter ) , @xmath9 , @xmath10 , and @xmath11  @xcite are presented for minimum - bias au - au collisions at @xmath3 = 20 , 62 , 130 and 200  gev and for cu - cu at 62 and 200  gev . in this study centrality
is based on charged particle production at midrapidity .
other centrality definitions appear in the literature but will not be considered here . however , the mcg model and analysis methods presented here and in ref .
@xcite can be readily extended to any practical centrality determination method .
the mcg model and updated input parameters are discussed in sec .
[ sec : model ] .
the accuracy of the model for describing @xmath6 data from rhic is demonstrated in sec .
[ sec : datafit ] .
mcg predictions and errors are presented and discussed in sec .
[ sec : results ] .
comparison of the mcg results and analytic parametrizations from ref .
@xcite are presented in sec .
[ sec : analytic ] . in sec .
[ sec : discussion ] we demonstrate via simulations how energy deposition data from transverse particle production , which are generally available at the trigger level from the rhic experiments  @xcite , can be used to further reduce the impact of background contamination and to mitigate the effects of event loss due to trigger and collision vertex finding inefficiencies . a summary and conclusions
are presented in sec .
[ sec : summary ] .
computational details are contained in three appendixes at the end .
the monte carlo glauber collision model used here is based on a standard set of assumptions  @xcite which are appropriate for ultra - relativistic heavy ion collisions .
these assumptions include the characterization of the collision in terms of a classical impact parameter @xmath12 , straight - line propagation of each incident nucleon through the oncoming nucleus , and a fixed transverse interaction range determined by the nucleon - nucleon ( n - n ) total inelastic cross section ( @xmath13 ) in free space .
the impact parameter was selected at random and the positions of the nucleons relative to the geometrical centers of the colliding nuclei were randomly distributed according to a spherical density @xmath14 .
center - of - mass constraints were not imposed on the nucleon positions .
nucleon pairs in the colliding nuclei were assumed to interact hadronically if their relative impact parameter was @xmath15 .    charged hadron multiplicity was assigned using the phenomenological two - component model of kharzeev and nardi  @xcite for `` soft '' plus `` hard '' particle production processes where the mean number of charged hadrons in the acceptance ( @xmath16 ) per unit pseudorapidity ( @xmath17 ) was computed according to & = & ( 1-x)n_pp + x n_pp n_bin .
[ eq1 ] parameters @xmath18 and @xmath19 depend on collision energy .
event - wise multiplicities were obtained by sampling the gaussian distribution  @xcite ( n_ch,|n_ch ) & = & ( - ) [ eq2 ] where @xmath20 is a multiplicity fluctuation width parameter .
alternate distributions ( _ e.g. _ poisson , negative binomial  @xcite ) produce quantitative effects on the multiplicity frequency distribution for very peripheral collisions but do not affect the present estimates of systematic errors . @xmath1 and
@xmath4 were both required to be @xmath21 in order for the simulated collision to be used in the analysis .
centrality bins were defined using the multiplicity frequency distribution @xmath6 .
centrality bin average quantities , @xmath8 , @xmath9 , @xmath10 and @xmath22 , were calculated using the events within each bin .
the acceptance for this study was @xmath23 and full @xmath24 in azimuth .
monte carlo glauber simulations require the distribution of the centers of the nucleons in the nuclear ground state , @xmath25 , the point - matter density . for @xmath26cu and @xmath27au
these were estimated using the measured charge densities  @xcite and the hartree - fock calculations of negele using the density matrix expansion ( dme ) framework  @xcite for the neutron - proton density differences .
the charge densities for @xmath26cu and @xmath27au were represented by a woods - saxon distribution , ( r ) & = & _ 0 \ { 1 + } ^-1 , [ eq3 ] where the radius and diffuseness parameters are listed in table  [ tablei ] .
@cllll density & & + parameter & empirical@xmath28 & dme@xmath29 & empirical@xmath28 & dme@xmath29 + @xmath30 & 4.214@xmath310.026 & 4.232 & 6.38@xmath310.06 & 6.443 + @xmath32 & 0.586@xmath310.018 & & 0.535@xmath310.027 & + @xmath33 & 3.925@xmath310.022 & 3.899 & 5.33@xmath310.05 & 5.423 + @xmath34 & 4.195@xmath310.085 & 4.213 & 6.43@xmath310.10 & 6.495 + @xmath35 & 0.581@xmath310.031 & & 0.568@xmath310.047 & + @xmath36 & 3.901 & 3.875 & 5.41 & 5.502 +    @xmath28 charge density results based on electron scattering analysis  @xcite ; estimates of point matter densities as discussed in the text .
@xmath29 density matrix expansion predictions  @xcite .
the point matter densities assumed for the present analysis are spherically symmetric with a woods - saxon radial distribution where the half - density and rms radii were estimated by adding the hartree - fock dme point matter  charge distribution differences to the measured charge density radii where , c_pt , m & = & c_chrg + _ dme , + r^2_pt , m ^1/2 & = & r^2_chrg ^1/2 + _ dme , [ eq4 ] where @xmath30 and @xmath33 denote the measured radii from electron scattering while the quantities in the square brackets are the dme predictions .
the diffusivity parameter @xmath35 was obtained from the quantities in eq .
( [ eq4 ] ) and the approximate relation , @xmath37 $ ]  @xcite . the nominal radius and diffusivity parameters assumed here for au and cu
are listed in table  [ tablei ] .
the uncertainties in @xmath34 and @xmath35 were obtained by summing the independent errors in the proton  @xcite and neutron point matter densities in quadrature .
the latter is constrained by theoretical and experimental information about the neutron - proton density difference .
theoretical nuclear structure model predictions for the neutron@xmath38proton rms radii differences for isotopes in the cu and au mass range agree to within about @xmath39  fm  @xcite .
in general the theoretical predictions agree with analyses of medium energy proton - nucleus elastic scattering data  @xcite which are typically uncertain by about @xmath40  fm for isotopes in the cu and au mass range .
the uncertainty in the neutron density rms radii relative to the proton density was therefore assumed to be @xmath41  fm and the corresponding uncertainty in the matter rms radii was @xmath42  fm , where @xmath43 are the number of neutrons and nucleons in the isotope .
theoretical contributions to the errors in @xmath34 and @xmath35 were conservatively estimated by requiring each to independently account for the @xmath31(n / a)0.09  fm error .
the latter theoretical errors were combined in quadrature with the corresponding errors for @xmath30 and @xmath32 from analysis of electron scattering data to obtain the final errors listed in table  [ tablei ] .
the n - n total inelastic cross sections used here were based on total cross section measurements for p - p collisions ( @xmath311  mb uncertainty ) and elastic total cross sections for p - p and p-@xmath44 ( @xmath310.5  mb error )  @xcite .
proton - neutron total cross section data are not available in the energy range studied here .
the results for energies @xmath45 = 20 , 62 , 130 and 200  gev are respectively 33 , 35.3 , 38.7 and 41.7  mb , each being uncertain by @xmath31  1.1  mb .
it is possible that the effective interaction cross section between colliding nucleons inside a nucleus differs from that in free space ( density dependence ) or that the strength and range of the effective n - n interaction changes with each successive collision as in the `` used '' nucleon scenario  @xcite or the `` strict '' participant scaling model  @xcite ( limiting case of the used - nucleon model in which nucleons interact only once ) .
the study of density dependent effects is well beyond the scope of the present analysis .
parameters @xmath46 at @xmath45 = 62 , 130 and 200  gev are 2.01 , 2.25 and 2.43 ( @xmath47 error for each ) , respectively , using the energy dependent parametrization of the ua5  @xcite and cdf data given in ref .  @xcite . for the @xmath45 = 20  gev data , which is outside the energy range
parametrized in  @xcite , the average of isr  @xcite and fnal  @xcite measurements summarized in ref .
@xcite was assumed where @xmath48 .
binary scattering parameter @xmath19 in eq .
( [ eq1 ] ) was estimated by fitting @xmath49 data from star  @xcite , phenix  @xcite , and phobos  @xcite versus @xmath2 with eq .
( [ eq1 ] ) rewritten as , & = & n_pp[1+x(-1 ) ] , [ eq5 ] assuming the above values for @xmath50 .
the resulting values of @xmath19 are 0.07 , 0.09 , 0.09 and 0.13 ( @xmath310.03 errors for each ) for the 20 , 62 , 130 and 200  gev data , respectively .
analysis of the 19.6  gev au - au data by the phobos experiment  @xcite assumed @xmath51  @xcite and obtained @xmath52 , leading to claims in the literature that hard scattering contributions to multiplicity do not change with collision energies from 20 to 200  gev  @xcite .
a more recent compilation  @xcite of @xmath50 measurements in this energy range indicates a larger value of @xmath48 .
fitting eq .
( [ eq5 ] ) to the combined phenix  @xcite and phobos  @xcite 19.6  gev data with @xmath53 resulted in a smaller @xmath54 .
phobos 62  gev au - au data  @xcite , when plotted versus @xmath2 , do not linearly extrapolate to @xmath50 in contrast to what is expected from fig .  4 of ref .
@xcite and fig .  1 of ref .
@xcite from the phobos experiment .
however , 62  gev au - au results from the star experiment linearly extrapolate to @xmath50 at @xmath55 , resulting in the fitted value @xmath56 used here .
@xmath19 parameters from the three rhic experiments generally agree for the 130 and 200  gev data ; the phobos values of @xmath56  @xcite and 0.13  @xcite at 130 and 200  gev , respectively , were confirmed and used here .
it is possible that parameter @xmath19 is affected by other processes in addition to hard and semi - hard partonic scattering in the initial collision stage where the latter mechanisms would be expected to follow a @xmath57 dependence .    the variance of the multiplicity distribution @xmath58 in eq .
( [ eq2 ] ) is given by @xmath59 where values of parameter @xmath60 or @xmath61 represent multiplicity fluctuation suppression or excess , respectively , relative to pure statistical fluctuations ( @xmath62 ) .
in general a non - vanishing integral  @xcite of two - particle correlations  @xcite over the acceptance requires @xmath63 . in principle , correlation measures of the type reported in ref .
@xcite could be used to determine parameter @xmath20 using the relationship between fluctuations and correlations developed in refs .
@xcite and discussed in appendix  a. on the other hand , the kharzeev and nardi two - component multiplicity model with distribution @xmath58 constitutes a phenomenology for describing event - wise multiplicity frequency distribution data . the phenomenological approach based on fits to multiplicity distribution data , shown in the next section , was used to estimate parameter @xmath20 in the present analysis .
the uncertainty in the width of @xmath58 was estimated by fitting the data .
the uncertainty in the analytic form of the multiplicity fluctuation distribution was accounted for by comparing monte carlo glauber results assuming the gaussian @xmath58 in eq .
( [ eq2 ] ) ( with @xmath62 ) with results assuming a negative binomial distribution ( nbd )  @xcite in place of eq .
( [ eq2 ] ) as explained in appendix  a.
the 130  gev au - au minimum - bias negatively charged hadron multiplicity frequency distribution data from star  @xcite for 60k events were fit by adjusting parameters @xmath50 ( for negative hadrons ) and @xmath20 where @xmath19 was fixed to 0.09 as discussed above .
the acceptance was defined by transverse momentum @xmath64
gev/@xmath65 , @xmath66 , and @xmath67 .
the fit ( solid histogram ) is shown in the left panel of fig .
[ fig1 ] in comparison with data ( solid dots ) , where the optimum values of @xmath50 and @xmath20 are @xmath68 ( consistent with @xmath69 for charged particle yields ) and @xmath70 , respectively .
similar fits ( not shown ) to the charged particle minimum - bias multiplicity distribution from the same data set resulted in @xmath71 .
parameter @xmath20 was therefore set to @xmath72 for all four energies .        in the tp analysis
@xcite it was shown that minimum - bias multiplicity frequency distributions for relativistic heavy ion collision experiments and monte carlo glauber models approximately follow a power - law distribution .
this is illustrated by plotting the data and mcg fit from the left panel of fig .
[ fig1 ] on log - log axes in the middle panel . except near the end - points the data can be described to within 10% with a slope ( exponent ) of approximately @xmath73 .
the collision event yield is approximately proportional to @xmath74 .
therefore distribution @xmath75 versus @xmath76 is approximately constant , where & = & = 4n_h^-^3/4 as shown in the right - hand panel of fig .
[ fig1 ] for data ( solid dots ) and mcg fit ( solid line ) .
the power - law plotting format in the right - hand panel enables a more sensitive comparison between the model fit and data than in the usual semi - log format ( left panel ) .
the mcg model fit is consistent with the data except for the first two data points at low @xmath76 . except near the end - points the @xmath75 data are constant on @xmath76 to within an overall variation of 20% in comparison to the @xmath77 data which span nearly two orders of magnitude within the same multiplicity range .
the results in fig .
[ fig1 ] demonstrate the efficacy of the monte carlo glauber model with two - component multiplicity production for accurately describing the measured multiplicity frequency distributions at rhic .
based on this outcome we conclude that the present model is reasonable to use in estimating centrality bin average quantities and their systematic errors for rhic data .
ensembles of one - million , minimum - bias ( _ i.e. _ random impact parameter ) monte carlo collisions were generated for each of the six systems studied here .
centrality bin averaged quantities @xmath78 , @xmath8 , @xmath9 , @xmath10 and @xmath2 are listed in tables  [ tableii]-[tablevii ] using the nominal parameter values discussed above .
statistical errors are typically 0.1 - 0.2% and always @xmath79% of the nominal bin averages and in all instances are much less than the systematic errors discussed below . results in table  [ tableii ] for @xmath9 for au - au collisions at 200  gev are systematically larger than published monte carlo glauber predictions in ref .
@xcite by approximately 5% for peripheral centrality bins .
this systematic increase is primarily caused by multiplicity fluctuations in the present model and basing centrality on @xmath4 here rather than on @xmath0 as was done in ref .
@xcite .
average charged particle multiplicity per participant pair for @xmath66 at mid - pseudorapidity as a function of centrality ( @xmath2 ) is shown in fig .
[ fig2 ] for monte carlo simulated au - au collisions at @xmath3 = 200  gev ( solid dots ) .
the p - p limit , @xmath80 , is indicated by the solid square symbol .
the data display a linear dependence on @xmath81 as expected from eq .
( [ eq5 ] ) except for the most - peripheral 90 - 100% centrality bin where multiplicity fluctuations significantly reduce the average @xmath4 .
the slope agrees with @xmath82 ; the linear extrapolation to the p - p limit is evident .
( color online ) monte carlo glauber results for @xmath83 for one - million au - au collisions at @xmath3 = 200  gev using the nominal parameters discussed in the text ( solid dots ) .
linear fit ( dashed line ) to the au - au results ( excluding the 90 - 100% centrality bin ) accurately extrapolates to the p - p multiplicity from ua5  @xcite .
monte carlo data for au - au include statistical errors only which are smaller than the symbols.,width=288 ]      the sources of uncertainty which cause systematic error in the collision geometry quantities can be organized into three categories . the first includes fitting uncertainties that occur when the corrected , multiplicity frequency distribution data are described by adjusting the parameters of the multiplicity production model in eqs .
( [ eq1]-[eq2 ] ) .
the second includes the uncertainties in the matter density , nucleon - nucleon total inelastic cross section , and functional representation of the multiplicity fluctuation distribution @xmath58 .
the third includes the uncertainty in the corrected multiplicity frequency distribution data @xmath6 .
the latter arise from trigger inefficiency , collision vertex finding inefficiency , background contamination , and particle tajectory finding inefficiency .
systematic errors in centrality bin average quantities due to fitting uncertainty were estimated by comparing the nominal results to that obtained with parameters @xmath19 and @xmath20 in eqs .
( [ eq1]-[eq2 ] ) individually varied within their respective fitting errors ( secs .
[ subsec : twocomponent ] and [ sec : datafit ] ) . uncertainty in parameter @xmath50
directly affects @xmath78 but does not affect @xmath8 , @xmath9 , @xmath10 or @xmath2 .
error estimates due to uncertainties in the matter density radius and diffuseness parameters and the n - n inelastic cross section require a refitting of the @xmath6 distribution .
this was accomplished via a phenomenological adjustment of the multiplicity production model described in appendix  b. fit recovery via @xmath84 minimization for mcg simulations with ample statistics ( of order @xmath85 collisions ) is computationally intensive .
the method in appendix  b is fast and accurate .
errors due to the uncertainty in the mathematical representation of @xmath58 were estimated by comparing the nominal mcg results using the gaussian distribution in eq .
( [ eq2 ] ) with results assuming a negative binomial distribution  @xcite as explained in appendix  a.    uncertainties in the @xmath6 data arising from trigger inefficiency , collision vertex finding inefficiency , and background contamination mainly affect the low multiplicity region . minimum - bias trigger efficiencies at rhic are @xmath86%  @xcite , @xmath87%  @xcite , 96%  @xcite and @xmath88%  @xcite . trigger inefficiency causes the lower @xmath4 half - max point of the @xmath89 distribution to shift to larger @xmath90 .
the position of this point for the @xmath89 distribution when corrected for trigger inefficiency ( assuming the power - law behavior ) has a relative uncertainty of about @xmath3110% .  2 to @xmath313% .
the fractional uncertainty in the lower half - max position is @xmath31(0.02 to 0.03)@xmath91% where @xmath92 and @xmath93 are the lower and upper half - max points of the @xmath89 distribution .
] however , knowledge that the lower half - max point is constrained by p - p scattering allows the relative uncertainties of the lower half - max point on @xmath90 to be reduced to 1/4 of the uncertainties in @xmath50 , or to @xmath310.8% , @xmath310.9% , @xmath311% , and @xmath312% for the 200 , 130 , 62 and 20  gev data , respectively .
primary collision vertex reconstruction efficiency is approximately 100% for events with @xmath4 of order a few tens and greater but falls precipitously for @xmath94  @xcite .
phobos  @xcite and star  @xcite both report uncertainties in their overall vertex finding efficiencies of about @xmath95% for minimum - bias collisions .
vertex finding inefficiency increases the slope and half - max position of the low multiplicity edge of the @xmath89 distribution .
an uncertainty of @xmath95% in the overall efficiency results in about @xmath3110% uncertainty in the low @xmath90 half - max position .
however , the p - p data constrain this uncertainty . in this analysis we assumed 100% vertex finding efficiency for collisions producing @xmath96  @xcite , or @xmath97 , and allowed the slope of the lower edge of the @xmath89 distribution to vary such that the half - max position varied by @xmath310.8% , @xmath310.9% , @xmath311% , and @xmath312% for the 200 , 130 , 62 and 20  gev data , respectively , as in the preceding paragraph .
the principle sources of background contamination are from ultra - peripheral two - photon interactions  @xcite and beam - gas collisions .
the former process corresponds to coherent photon - photon interactions which excite both nuclei , followed by neutron decay ( which activates the minimum - bias trigger detectors ) and accompanied by resonance(s ) production which decays into charged particles transverse to the beam direction .
transverse particle multiplicities from upc events are typically @xmath98 ( _ e.g. _ @xmath99-meson decay ) and generally @xmath100  @xcite for @xmath101 at midrapidity .
upc backgrounds are therefore restricted to the most - peripheral ( 90 - 100% ) centrality bin .
upc yields should be approximately proportional to @xmath102  @xcite ( charge numbers for colliding ions 1 and 2 ) whereas beam - gas contamination should scale with beam current
. other background events , _
e.g. _ mutual coulomb dissociation processes such as @xmath103 + a @xmath104 b + n for both nuclei , can be eliminated by requiring minimum transverse particle production .
remaining background events will appear near the lower @xmath4 edge of the @xmath89 distribution .
estimates of background contamination in minimum - bias trigger data at rhic range from @xmath105%  @xcite to 6%  @xcite overall for au - au at @xmath3 = 130  gev corresponding to 0 - 20% and 60% , respectively , of the hadronic collision event yield in the 90 - 100% centrality bin .
most of the upc events occur at @xmath106 and can be eliminated by cuts on the number of transverse charged particles .
the remaining background contributions for @xmath107 are less than the amounts listed above . for the present analysis background contamination
was assumed to be dominated by upc events and to diminish in magnitude with collision energy and @xmath102 . in the present analysis background contamination
was applied to the nominal 90 - 100% centrality bin assuming 3% , 2% and 1% overall contamination levels in the au - au data at @xmath3 = 200 , 130 and 62  gev , respectively .
upc contamination for au - au collisions at 20  gev and for cu - cu at 200 and 62  gev was estimated to be negligible .
however , calculations were done for the latter three systems assuming a 1% overall background contamination ( 10% contamination within the nominal 90 - 100% centrality bin ) in order to provide a reference for further systematic error estimation .
the above sources of systematic uncertainty primarily affect the position , slope and shape of the @xmath89 distribution at the lower @xmath90 end - point .
these changes impact the mcg model which must describe those data and , in turn , the collision geometry quantities like @xmath9 .
the phenomenological method in appendix  b was used to estimate the changes in the collision geometry quantities relative to the nominal values .
the latter differences were taken as the estimated errors .
background levels in collider experiments vary significantly depending on beam quality , beam current and interaction rate .
excessive trigger backgrounds beyond that considered here may result from beam - gas interactions during periods of high integrated beam currents .
collision event pile - up in the particle tracking detectors during periods of high luminosity adversely affect collision vertex finding .
either condition may be so severe as to preclude access to the low multiplicity range of the minimum - bias distribution .
even so , the power - law dependence and p - p end - point constraints enable accurate centrality estimates and collision geometry assignments to be made for the remaining minimum - bias data .
charged particle trajectory reconstruction efficiencies in the large acceptance tracking detectors at rhic are typically 70 - 95%  @xcite and decrease approximately linearly with particle density in the detectors by about 20% from most - peripheral to most - central collisions  @xcite .
uncertainties in the assumed track reconstruction efficiencies were estimated by comparing corrected data for @xmath108 versus @xmath0 between the rhic experiments for au - au collisions at 20  gev  @xcite 130  gev  @xcite and 200  gev  @xcite .
the comparisons indicate uncertainties in both the overall tracking efficiencies and in the dependence on particle density in the tracking detectors .
the latter variation in efficiency from peripheral to central collisions is about @xmath109% .
overall changes in tracking efficiency are compensated in the mcg model by multiplicative adjustments to parameter @xmath50 and therefore have no effect on the centrality measures reported here .
changes in the assumed tracking efficiency dependence on @xmath4 affect the lower and upper end - point positions of the @xmath89 distribution and distort its shape .
the distortions must be accounted for by the mcg model in order to determine the net effect on the collision geometry measures .
the effects of the 8% uncertainty in the multiplicity dependence of the trajectory reconstruction efficiency were estimated by generalizing parameter @xmath50 in eq .
( [ eq1 ] ) to @xmath110 where @xmath111 and @xmath112 for most - central au - au collisions .
the same value for @xmath113 was assumed for cu - cu .
systematic errors were estimated as the differences between the nominal centrality bin averages and those obtained assuming @xmath111 .
the choice to maximize multiplicity variation for most - central collisions is arbitrary . as a result
systematic errors in centrality bin average multiplicities were not included in tables  [ tableii]-[tablevii ] .
however , this ambiguity does not affect the resulting systematic errors in the other centrality measures reported here .
the combined systematic errors ( all components added in quadrature ) in both magnitude and relative percent ( given in parentheses ) are listed in tables  [ tableii ] - [ tablevii ] for the six collision systems studied here .
impact parameter uncertainty is about @xmath114% for au - au and @xmath115% for cu - cu .
uncertainty in @xmath9 is about @xmath116% for au - au and @xmath117% for cu - cu for peripheral centralities , reducing to about @xmath118% or less for central collisions where full geometrical overlap of the two nuclei suppresses the dependence of @xmath0 on variations in the nuclear surface geometry .
systematic errors in @xmath10 vary from about @xmath119% for au - au collisions and @xmath120% for cu - cu for central to peripheral collisions .
errors in @xmath2 ( @xmath121% or less ) are suppressed due to covariation of @xmath9 and @xmath10 .
individual contributions to the systematic errors in collision geometry bin averages @xmath8 , @xmath9 , @xmath10 and @xmath2 for au - au and cu - cu collisions at @xmath3 = 200  gev are summarized in tables  [ tableviii ] and [ tableix ] , respectively .
similar results were obtained for the other four collision systems .
errors are given in percent where the absolute values of positive and negative errors were averaged together .
the three numbers listed for each instance correspond to the average errors within the 60 - 100% , 20 - 60% and 0 - 20% centrality bins , respectively .
the dominant errors are due to uncertainties in the matter density and the n - n total inelastic cross section .
errors due to uncertainties in the analytic form of the phenomenological multiplicity fluctuation model are only significant for peripheral collisions . errors due to uncertainties in the multiplicity dependent particle track reconstruction efficiency and in the trigger and vertex reconstruction inefficiencies are negligible when constrained by p - p data as shown previously by trainor and prindle  @xcite .
percent uncertainties due to possible background contamination are listed for all six collision systems in table  [ tablex ] .
the largest errors occur in the most - peripheral centrality bin as expected .
the absolute magnitudes of the errors decline rapidly with increasing centrality and are negligible for the centrality bins not listed in the table .
reference errors for au - au at 20  gev and cu - cu at @xmath3 = 200 and 62  gev were based on an assumed 10% background contamination in the 90 - 100% centrality bin .
estimates of total systematic error when background contamination differs from that assumed here can be obtained by removing the error contributions in table  [ tablex ] from the total errors in tables  [ tableii ] - [ tableiv ] for au - au collisions at 200 , 130 and 62  gev , respectively , and then adding ( in quadrature ) the appropriately scaled background errors from table  [ tablex ] . for au - au collisions at 20  gev and cu - cu collisions at 200 and 62
gev the scaled errors from table  [ tablex ] should be combined in quadrature with the total systematic errors in tables  [ tablev ] - [ tablevii ] .
if much larger backgrounds than those assumed here are encountered , then the present monte carlo results should not be scaled ; rather the errors should be recalculated .
overall , the systematic errors in collision geometry bin averages for non - peripheral centralities are dominated by uncertainties in the nuclear geometry and @xmath13 .
errors in the more peripheral bins are dominated by background contamination and ambiguities in the analytic form of the multiplicity fluctuation model .
the power - law description of heavy - ion collision centrality developed by trainor and prindle  @xcite prescribes simple , analytic parametrizations of the @xmath0 , @xmath1 , @xmath2 and @xmath4 dependences on total inelastic cross section fraction @xmath122 .
for @xmath0 the running integral relation on @xmath123 is accurately given by  @xcite ( 1 - /_0 ) & = & , [ eq7 ] where ( n_part/2)^ & = & ( n_part , min/2)^ /_0 + ( n_part , max/2)^ ( 1 - /_0 ) .
[ eq8 ] similarly n_bin^ & = & n_bin , min^ /_0 + n_bin , max^ ( 1 - /_0 ) , [ eq9 ] where subscripts @xmath124 and @xmath125 refer respectively to the lower and upper half - max end - point positions of the @xmath126 and @xmath127 distributions .    the above parametrizations are shown as the solid lines in the two upper panels of fig .
[ fig3 ] in comparison with the present mcg model results ( solid dots ) from table  [ tableii ] for au - au collisions at @xmath3 = 200  gev .
the end - point parameters for the power - law parametrizations were @xmath128 , @xmath129 , @xmath130 and @xmath131 .
the parametrization for @xmath132 is compared with the mcg result in the lower - left panel .
the average multiplicities using eq .
( [ eq5 ] ) , the values for @xmath50 and @xmath19 from sec .
[ subsec : twocomponent ] , and the above parametrizations for @xmath0 and @xmath2 ( solid curve ) is compared with the mcg result in the lower - right panel .
the simple power - law parametrizations introduced in ref .
@xcite quantitatively describe the mcg results , thus confirming their utility and precision .
calculation of the analytic power - law parametrizations require the upper half - max end - point values for @xmath0 and @xmath1 from the mcg model .
these quantities are listed in table  [ tablexi ] for the six collision systems studied here .
very peripheral collision data from the four rhic experiments have generally remained unpublished due to concerns with possibly large , and not well understood background contamination , trigger inefficiencies , and primary collision vertex finding inefficiencies .
trainor and prindle  @xcite originally showed , and the present analysis confirms , that the power - law dependence of the multiplicity frequency distribution data together with knowledge of the proton - proton multiplicity enables accurate centrality information to be obtained in spite of these uncertainties . nevertheless , upc events in principle cause the lower end - point portion of the multiplicity frequency distribution for minimum - bias a - a collisions to differ from the p - p limit and these backgrounds will contaminate the a - a data in the most - peripheral centrality bins . in this section
we discuss additional analysis methods using data available to the rhic experiments to minimize backgrounds and to provide contamination level estimates for the accepted centrality bins .
transverse particle production information via scintillators , silicon detectors , and/or calorimetry are available in the rhic experiments at the trigger level  @xcite .
these data record total energy deposition in the sensitive detector material for each ( minimum bias ) triggered collision event including contributions from a - a hadronic collisions , ultra - peripheral collisions , beam - gas and other backgrounds , etc . but do not include particle track finding and primary collision vertex finding inefficiencies . if the integrated yields from these detectors , denoted as @xmath133 , is approximately proportional to @xmath4 , then the frequency distribution will follow an approximate power - law such that the data can be usefully represented by @xmath134 versus @xmath135 , where @xmath136 is the number of triggers .
mcg simulations for the integrated yields from the star central trigger barrel ( ctb )  @xcite plastic scintillator detector ( @xmath137 , @xmath24 azimuth coverage ) were done for p - p and au - au minimum - bias collisions at @xmath3 = 200  gev .
details of this simulation are discussed in appendix  c. the results for @xmath134 versus @xmath135 are shown in fig .
both the p - p and au - au simulations reproduce the general shapes of measured star ctb minimum - bias trigger yields  @xcite .
the simulated au - au ctb yields follow the @xmath138 power - law distribution very well where the lower half - max point coincides with the mode of the p - p distribution .
trigger inefficiency in au - au causes a depletion at low multiplicity as indicated by the green , dashed curve .
over- or under - corrected yields result in the blue dotted curves . from fig .
[ fig4 ] it is clear how the p - p trigger yield can be used to constrain corrections for trigger inefficiency in a - a .
( color online ) monte carlo simulation results for the minimum - bias trigger frequency distribution for transverse particle production integrated trigger detector yield quantity @xmath133 to the 1/4 power .
simulated yields are shown for 1 m 200  gev au - au ( solid black curve , multiplied by 3 ) and 1 m proton - proton ( solid red curve ) collisions .
the yield for au - au collisions when trigger inefficiency is large is illustrated by the lowest ( right - most ) green dashed curve ( hand - drawn sketch ) .
over- and under - corrected yields are similarly illustrated by the upper ( left - most ) and middle dotted blue lines , respectively .
background contamination contributes at lower multiplicity as illustrated by the dashed - dotted magenta curve ( hand - drawn sketch ) .
, width=384 ]    background triggers appear at lower values of @xmath133 where upc events typically produce of order two charged particles at mid - rapidity compared with the average from p - p collisions of about @xmath139 in @xmath137 and therefore appear as a peak or enhancement ( dashed - dotted curve in fig .
[ fig4 ] ) below the mode of the p - p distribution .
comparisons between the measured p - p and au - au @xmath134 distributions and the power - law simulations for au - au provide a reasonable means for defining additional event cuts to reduce background contamination .
this information can also be used to estimate the contamination level for accepted events .
measurements of the dependence of the lower end - point region of the @xmath134 distribution on integrated beam current , luminosity , and ion species would help disentangle background contributions from beam - gas collisions , event pile - up , and upc events
. a similar power - law analysis using reconstructed particle tracks , but without a primary collision vertex requirement , would permit the vertex finding inefficiencies to be estimated and corrected .
a monte carlo glauber and two - component multiplicity production model with fluctuations for high energy heavy ion collisions was used to describe minimum - bias multiplicity frequency distribution data from rhic .
updated woods - saxon radii and diffusivities for the point matter densities of @xmath27au and @xmath26cu were determined using a combination of charge density measurements from electron scattering and theoretical hartree - fock predictions .
the binary scattering parameter @xmath19 was estimated using a compilation of available rhic data from 20 to 200  gev .
the values for @xmath19 obtained here display some energy dependence , although with large uncertainty .
the model was fitted to the 130  gev au - au minimum - bias data  @xcite in both the conventional semi - log format and in the more sensitive power - law inspired format introduced in ref .
quantitative agreement between the present model and data was obtained .
systematic errors in centrality bin averages for the geometrical quantities @xmath8 , @xmath9 , @xmath10 and @xmath11 were estimated based on charged particle multiplicity for @xmath23 and full @xmath24 azimuth acceptance .
the sources of systematic error considered here included uncertainties associated with the model ( nuclear density , nucleon - nucleon inelastic total cross section , multiplicity production model analytic form and parameters ) and with the corrected minimum - bias multiplicity frequency distribution data ( background contamination , trigger and collision vertex finding inefficiencies , and particle trajectory reconstruction inefficiencies ) .
the tp centrality analysis method was applied to the mcg predictions and error analysis for six collision systems relevant to the rhic heavy ion program .
the analysis showed that significantly reduced errors in the collision geometry quantities result when the uncertainties in the minimum - bias multiplicity frequency distribution data are constrained by the empirical power - law behavior and the minimum - bias p - p collision data .
the reduction in errors is in agreement with the original tp power - law analysis ; the resulting errors are in general smaller than previous estimates  @xcite which did not utilize the power - law and p - p constraints . the accuracy of simple parametrizations of centrality dependent quantities developed by trainor and prindle  @xcite was confirmed .
we also discussed how particle production data at the trigger level ( _ e.g. _ scintillator hits and calorimeter energy depositions ) can be used within the power - law and p - p constraint methodology  @xcite to define additional event cuts which minimize background contamination and which enable the detrimental effects of trigger inefficiencies to be reduced .    the systematic errors presented in this paper demonstrate the accuracy of collision geometry quantities which can be achieved with rhic data by exploiting the power - law behavior of the a - a data and using constraints from minimum - bias p - p collision data .
the results presented here assume that backgrounds from upc events , pileup , beam - gas collisions , etc .
are not too large and that the minimum - bias trigger and collision vertex reconstruction efficiencies are not too small for peripheral collisions .
it is intended that the mcg results presented here will serve as a useful resource for the rhic community and that the analysis method developed by trainor and prindle  @xcite and applied in this paper , together with the trigger data analysis method discussed here , will enable accurate description of , and better access to the heavy ion peripheral collision data from rhic .
the authors gratefully acknowledge helpful discussions , detailed comments and the plots in fig .
[ fig3 ] from prof .
t. a. trainor ( univ . of washington ) as well as helpful discussions and comments from drs .
j. c. dunlop and zhangbu xu ( both from brookhaven national lab ) .
this research was supported in part by the u. s. dept . of energy grant
de - fg03 - 94er40845 .
cccccc + centrality ( % ) & @xmath78 & @xmath8 & @xmath9 & @xmath10 & @xmath2 + 90 - 100 & 2.9 & @xmath140 & @xmath141 & @xmath142 & @xmath143 + 80 - 90 & 8.2 & @xmath144 & @xmath145 & @xmath146 & @xmath147 + 70 - 80 & 18.7 & @xmath148 & @xmath149 & @xmath150 & @xmath151 + 60 - 70 & 37.8 & @xmath152 & @xmath153 & @xmath154 & @xmath155 + 50 - 60 & 69.7 & @xmath156 & @xmath157 & @xmath158 & @xmath159 + 40 - 50 & 118.5 & @xmath160 & @xmath161 & @xmath162 & @xmath163 + 30 - 40 & 190.2 & @xmath164 & @xmath165 & @xmath166 & @xmath167 + 20 - 30 & 291.5 & @xmath168 & @xmath169 & @xmath170 & @xmath171 + 10 - 20 & 433.5 & @xmath172 & @xmath173 & @xmath174 & @xmath175 + 5 - 10 & 577.6 & @xmath176 & @xmath177 & @xmath178 & @xmath179 + 0 - 5 & 705.6 & @xmath180 & @xmath181 & @xmath182 & @xmath183 +    cccccc + centrality ( % ) & @xmath78 & @xmath8 & @xmath9 & @xmath10 & @xmath2 + 90 - 100 & 2.7 & @xmath184 & @xmath185 & @xmath186 & @xmath187 + 80 - 90 & 7.4 & @xmath188 & @xmath189 & @xmath190 & @xmath191 + 70 - 80 & 16.6 & @xmath192 & @xmath193 & @xmath194 & @xmath195 + 60 - 70 & 33.2 & @xmath196 & @xmath197 & @xmath198 & @xmath199 + 50 - 60 & 60.1 & @xmath200 & @xmath201 & @xmath202 & @xmath203 + 40 - 50 & 100.4 & @xmath204 & @xmath205 & @xmath206 & @xmath207 + 30 - 40 & 157.8 & @xmath208 & @xmath209 & @xmath210 & @xmath211 + 20 - 30 & 237.7 & @xmath212 & @xmath213 & @xmath214 & @xmath215 + 10 - 20 & 348.7 & @xmath216 & @xmath217 & @xmath218 & @xmath219 + 5 - 10 & 459.5 & @xmath220 & @xmath221 & @xmath222 & @xmath223 + 0 - 5 & 557.6 & @xmath224 & @xmath225 & @xmath226 & @xmath227 +    cccccc + centrality ( % ) & @xmath78 & @xmath8 & @xmath9 & @xmath10 & @xmath2 + 90 - 100 & 2.4 & @xmath228 & @xmath229 & @xmath230 & @xmath231 + 80 - 90 & 6.6 & @xmath232 & @xmath233 & @xmath234 & @xmath235 + 70 - 80 & 14.6 & @xmath236 & @xmath237 & @xmath238 & @xmath239 + 60 - 70 & 29.0 & @xmath240 & @xmath241 & @xmath242 & @xmath243 + 50 - 60 & 52.3 & @xmath244 & @xmath245 & @xmath246 & @xmath247 + 40 - 50 & 87.1 & @xmath248 & @xmath249 & @xmath250 & @xmath251 + 30 - 40 & 136.9 & @xmath252 & @xmath253 & @xmath254 & @xmath255 + 20 - 30 & 205.5 & @xmath256 & @xmath257 & @xmath258 & @xmath259 + 10 - 20 & 300.7 & @xmath260 & @xmath261 & @xmath262 & @xmath263 + 5 - 10 & 395.5 & @xmath264 & @xmath265 & @xmath266 & @xmath267 + 0 - 5 & 480.8 & @xmath268 & @xmath269 & @xmath270 & @xmath271 +    cccccc + centrality ( % ) & @xmath78 & @xmath8 & @xmath9 & @xmath10 & @xmath2 + 90 - 100 & 1.9 & @xmath272 & @xmath273 & @xmath274 & @xmath275 + 80 - 90 & 4.8 & @xmath276 & @xmath277 & @xmath278 & @xmath279 + 70 - 80 & 10.5 & @xmath280 & @xmath281 & @xmath282 & @xmath283 + 60 - 70 & 20.5 & @xmath284 & @xmath285 & @xmath286 & @xmath287 + 50 - 60 & 36.2 & @xmath288 & @xmath289 & @xmath290 & @xmath291 + 40 - 50 & 59.2 & @xmath292 & @xmath293 & @xmath294 & @xmath295 + 30 - 40 & 91.6 & @xmath296 & @xmath297 & @xmath298 & @xmath299 + 20 - 30 & 135.9 & @xmath300 & @xmath301 & @xmath302 & @xmath303 + 10 - 20 & 196.1 & @xmath304 & @xmath305 & @xmath306 & @xmath307 + 5 - 10 & 255.4 & @xmath308 & @xmath309 & @xmath310 & @xmath311 + 0 - 5 & 309.2 & @xmath312 & @xmath313 & @xmath314 & @xmath315 +    cccccc + centrality ( % ) & @xmath78 & @xmath8 & @xmath9 & @xmath10 & @xmath2 + 90 - 100 & 2.2 & @xmath316 & @xmath317 & @xmath318 & @xmath319 + 80 - 90 & 4.7 & @xmath320 & @xmath321 & @xmath322 & @xmath323 + 70 - 80 & 8.4 & @xmath324 & @xmath325 & @xmath326 & @xmath327 + 60 - 70 & 14.3 & @xmath328 & @xmath329 & @xmath330 & @xmath331 + 50 - 60 & 23.3 & @xmath332 & @xmath333 & @xmath334 & @xmath335 + 40 - 50 & 36.6 & @xmath336 & @xmath337 & @xmath338 & @xmath339 + 30 - 40 & 55.4 & @xmath340 & @xmath341 & @xmath342 & @xmath343 + 20 - 30 & 81.5 & @xmath344 & @xmath345 & @xmath346 & @xmath347 + 10 - 20 & 117.4 & @xmath348 & @xmath349 & @xmath350 & @xmath351 + 5 - 10 & 152.2 & @xmath352 & @xmath353 & @xmath354 & @xmath355 + 0 - 5 & 184.8 & @xmath356 & @xmath357 & @xmath358 & @xmath359 +    cccccc + centrality ( % ) & @xmath78 & @xmath8 & @xmath9 & @xmath10 & @xmath2 + 90 - 100 & 1.8 & @xmath360 & @xmath361 & @xmath362 & @xmath363 + 80 - 90 & 3.8 & @xmath364 & @xmath365 & @xmath366 & @xmath367 + 70 - 80 & 6.7 & @xmath368 & @xmath369 & @xmath370 & @xmath371 + 60 - 70 & 11.2 & @xmath372 & @xmath373 & @xmath374 & @xmath375 + 50 - 60 & 17.9 & @xmath376 & @xmath377 & @xmath378 & @xmath379 + 40 - 50 & 27.7 & @xmath380 & @xmath381 & @xmath382 & @xmath383 + 30 - 40 & 41.3 & @xmath384 & @xmath385 & @xmath386 & @xmath387 + 20 - 30 & 60.0 & @xmath388 & @xmath389 & @xmath390 & @xmath391 + 10 - 20 & 84.9 & @xmath392 & @xmath393 & @xmath394 & @xmath395 + 5 - 10 & 108.9 & @xmath396 & @xmath397 & @xmath398 & @xmath399 + 0 - 5 & 132.2 & @xmath400 & @xmath401 & @xmath402 & @xmath403 +    @ccccc error & + source & @xmath404(% ) & @xmath405(% ) & @xmath406(% ) & @xmath407(% ) + @xmath34 & 0.9 , 1.0 , 0.8 & 0.1 , 0.1 , 0.2 & 0.1 , 0.9 , 2.2 & 0.0 , 0.9 , 2.0 + @xmath35 & 1.7 , 1.6 , 1.5 & 3.2 , 1.6 , 0.3 & 5.2 , 5.1 , 2.6 & 2.1 , 3.5 , 2.3 + @xmath13 & 0.1 , 0.1 , 0.0 & 0.4 , 0.5 , 0.3 & 1.0 , 2.0 , 2.5 & 0.6 , 1.5 , 2.2 + @xmath19 & 0.0 , 0.0 , 0.1 & 0.2 , 0.0 , 0.0 & 0.3 , 0.0 , 0.0 & 0.2 , 0.1 , 0.0 + @xmath20 & 0.0 , 0.0 , 0.2 & 0.3 , 0.1 , 0.0 & 0.4 , 0.2 , 0.1 & 0.1 , 0.1 , 0.0 + @xmath408-form & 0.1 , 0.0 , 0.7 & 2.2 , 0.1 , 0.2 & 3.5 , 0.3 , 0.3 & 1.2 , 0.2 , 0.1 + trigger & 0.0 , 0.0 , 0.0 & 0.1 , 0.0 , 0.0 & 0.1 , 0.0 , 0.0 & 0.0 , 0.0 , 0.0 + vertex & 0.0 , 0.0 , 0.1 & 0.1 , 0.0 , 0.0 & 0.1 , 0.0 , 0.0 & 0.0 , 0.0 , 0.0 + tracking & 0.0 , 0.0 , 0.1 & 0.0 , 0.0 , 0.0 & 0.0 , 0.0 , 0.0 & 0.0 , 0.0 , 0.0 +    @ccccc error & + source & @xmath404(% ) & @xmath405(% ) & @xmath406(% ) & @xmath407(% ) + @xmath34 & 1.0 , 1.0 , 1.0 & 0.2 , 0.4 , 0.5 & 0.2 , 1.1 , 2.6 & 0.2 , 0.8 , 2.0 + @xmath35 & 1.9 , 1.8 , 1.7 & 2.1 , 1.8 , 0.8 & 3.2 , 4.3 , 3.2 & 1.2 , 2.5 , 2.3 + @xmath13 & 0.2 , 0.2 , 0.2 & 0.2 , 0.4 , 0.4 & 0.5 , 1.5 , 2.2 & 0.3 , 1.1 , 1.8 + @xmath19 & 0.0 , 0.0 , 0.0 & 0.0 , 0.0 , 0.0 & 0.2 , 0.0 , 0.1 & 0.2 , 0.1 , 0.1 + @xmath20 & 0.1 , 0.1 , 0.3 & 0.8 , 0.2 , 0.1 & 1.2 , 0.4 , 0.2 & 0.4 , 0.2 , 0.1 + @xmath408-form & 0.3 , 0.1 , 1.9 & 2.8 , 0.1 , 0.5 & 4.7 , 0.5 , 0.6 & 1.5 , 0.5 , 0.3 + trigger & 0.0 , 0.0 , 0.1 & 0.2 , 0.1 , 0.0 & 0.2 , 0.1 , 0.0 & 0.1 , 0.0 , 0.0 + vertex & 0.0 , 0.0 , 0.0 & 0.0 , 0.0 , 0.0 & 0.2 , 0.0 , 0.0 & 0.0 , 0.0 , 0.0 + tracking & 0.0 , 0.0 , 0.0 & 0.0 , 0.0 , 0.0 & 0.0 , 0.0 , 0.0 & 0.0 , 0.0 , 0.0 +    .[tablex ] percent changes in collision geometry bin average quantities relative to the reference values in tables  [ tableii ] - [ tablevii ] due to background contamination of the peripheral collision events as explained in sec .
[ sec : results ] .
background contamination errors for centrality bins not listed were negligible . [
cols="^,^,^,^,^,^ " , ]     @xmath28 average percent changes within the combined centrality bin .
@xmath29 reference uncertainties assuming 10% background contamination in the number of collisions in the nominal 90 - 100% centrality bin , equivalent to 1% overall background .    @ccc system & @xmath409 & @xmath410 + au - au 200 gev & 379.3 & 1166 + au - au 130 gev & 377.5 & 1082 + au - au 62 gev & 374.9 & 988 + au - au 20 gev & 373.4 & 926 +
cu - cu 200 gev & 115.8 & 230 + cu - cu 62 gev & 112.7 & 196 +    [ [ section ] ] -    multiplicity frequency distributions for high energy minimum - bias proton - proton collisions are well described by a negative binomial distribution ( nbd )  @xcite given by _ nbd(n ,
n , k ) & = & (    c n+k-1 + k-1 +    ) ( ) ^n ( 1 + n / k)^-k . [ eqa1 ] in eq .
( [ eqa1 ] ) @xmath411 and @xmath412 are random and mean multiplicities , the variance is @xmath413 , the variance excess relative to poisson statistics is @xmath414 , and the multiplicity fluctuation variance excess per final - state particle is @xmath415 .
the latter expression serves as an operational definition of parameter @xmath416 . for p - p collisions both @xmath412 and @xmath417 increase approximately linearly with ln(@xmath418 )  @xcite as expected from pqcd cross sections for parton scattering .
multiplicity fluctuation variance excess per final - state particle within a given acceptance is equal to the integral of two - particle correlations on relative momentum coordinates within the same acceptance  @xcite . for au - au collisions at 130  gev  @xcite and at 62 and 200  gev  @xcite
the integral of charged particle correlations for @xmath23 and @xmath24 azimuth acceptance is dominated by a two - dimensional ( 2d ) correlation peak at small relative opening angles .
the magnitude of this correlation integral increases roughly with @xmath419 .
multiplicity fluctuation variance excess per final - state particle was therefore approximated by & = & [ eqa2 ] where quantity @xmath420 $ ] is the variance excess for p - p collisions , @xmath421 is the average number of binary collisions per incident participant nucleon for the most - central ( 0 - 5% ) a - a data , and @xmath422 is taken to be the ratio of the 2d correlation peak amplitude for most - central a - a collisions to that for p - p collisions .
values of @xmath415 for p - p collisions were obtained from ua5  @xcite using either measured values of @xmath416 or the ua5 energy dependent fit function for @xmath416 .
@xmath412 for @xmath23 is given by the values of @xmath50 assumed in this analysis .
simulations demonstrated that the values of @xmath416 obtained by ua5 with very large @xmath17 acceptance remained approximately constant at the smaller @xmath17 acceptance assumed in this analysis .
values of quantities @xmath420 $ ] were 0.09 , 0.24 , 0.40 and 0.51 for @xmath45 = 20 , 62 , 130 and 200  gev , respectively .
@xmath422 for au - au collisions at 20 , 62 , 130 and 200  gev were estimated to be 4.8 , 5.0 , 5.4 and 5.5 , respectively , from analysis of correlation data at 130  gev  @xcite assuming ln@xmath45 scaling .
the remaining values for cu - cu at 62 and 200  gev were respectively 3.4 and 3.7 assuming @xmath81 scaling .
finally , for each simulated collision with event - wise @xmath2 and @xmath16 from eq .
( [ eq1 ] ) nbd distribution eq .
( [ eqa1 ] ) was sampled to obtain the event - wise multiplicity .
[ [ section-1 ] ] -    fit recovery of the nominal @xmath6 distribution following shifts in the nuclear density parameters and @xmath13 was achieved by constructing a perturbative correction , @xmath423 , which was added to the right - hand side of eq .
( [ eq1 ] ) .
parameter shifts affect the collision distributions on @xmath0 and @xmath1 and in turn the @xmath424 and @xmath6 distributions .
the fit recovery term is defined such that the nominal @xmath424 distribution is maintained when the frequency distribution on @xmath0 changes . running integrals of the distributions @xmath424 and @xmath425 produce a one - to - one correspondence between @xmath426 and @xmath427 as a function of the number of summed events . using the power - law representation this running integral relation
is given by , _
|n_ch , min^1/4^|n_ch^1/4 d|n_ch^  1/4 =
_ 1^(n_part/2)^1/4 d(n^_part/2)^1/4 [ eq13 ] which defines a locus of ordered pairs @xmath428 $ ] , which can be compactly expressed as a discrete function defined by _
ch^1/4 & = & n^1/4_ch , d^ d [ ( n_part/2)^1/4 ] .
[ eq14 ] in eqs .
( [ eq13 ] ) and ( [ eq14 ] ) subscripts @xmath429 and @xmath430 refer to the nominal ( reference ) distribution when @xmath431 and to the shifted distribution when @xmath432 .
lower limit @xmath433 for the reference distribution .
in order to recover the nominal @xmath424 distribution the correction term must be computed by the following difference , ( ) & = & n_ch , s r [ ( ) ^ ] - n_ch , s s [ ( ) ^ ] .
[ eq15 ] throughout this fit recovery procedure the parameters of @xmath58 in eq .
( [ eq2 ] ) remained fixed .    numerically stable results with one - million minimum - bias simulated collisions required coarse binning of the monte carlo distributions @xmath434 and @xmath435 prior to evaluation of eqs .
( [ eq13 ] ) - ( [ eq15 ] ) .
when parameters @xmath34 and @xmath13 were shifted excellent fit recovery was achieved with centrality bins corresponding to 0 - 5% , 5 - 10% , 10 - 40% , 40 - 70% and 70 - 100% total cross section fractions .
when diffusivity parameter @xmath35 was shifted the optimum binning was 0 - 5% , 5 - 10% , 10 - 20% , 20 - 40% , 40 - 60% , 60 - 80% and 80 - 100% . the power - law representation in eq .
( [ eq13 ] ) exploits the approximately uniform statistics throughout the domain of the integrands , making numerical stability of the results less problematic .
variations in the density parameters and @xmath13 produced typical changes in @xmath78 as a function of centrality of order several percent , up to 5% .
fit recovery in terms of @xmath78 was typically within a few tenths of a percent and always less than about 1% where the magnitude of @xmath423 relative to @xmath436 was less than about 1% , 4% and 1.5% for changes in @xmath34 , @xmath35 and @xmath13 , respectively , for all six collision systems .
a similar method was used to account for shifts in @xmath6 representing the systematic uncertainties in trigger inefficiency , background contamination , and collision vertex reconstruction inefficiency .
for these simulations @xmath437 remained fixed while @xmath438 varied .
the required correction is given by , ( ) & = & n_ch , r s [ ( ) ^ ] - n_ch , r r [ ( ) ^ ] .
[ eq16 ]    [ [ section-2 ] ] -    the monte carlo simulation model discussed in sec .  [
sec : discussion ] for transverse particle production and energy deposition in the star ctb for p - p and au - au minimum bias collisions is described here . for p - p collisions the interaction point along the beam - line
was randomly selected within @xmath3125  cm of the geometrical center of the star detector  @xcite .
the minimum - bias charged particle multiplicity was obtained by sampling the negative binomial distribution ( nbd ) fitted to the ua5 @xmath45 = 200  gev p - p data  @xcite with parameters @xmath439 and @xmath440 .
the @xmath5 , @xmath17 and particle species ( @xmath441 , k@xmath442 , proton or antiproton ) for each charged particle produced in the collision was sampled from the measured @xmath443 , @xmath444 , and particle species distributions , respectively .
distribution @xmath443 was assumed to be proportional to @xmath445 , where @xmath446 , @xmath447 is the particle rest mass , and inverse effective temperature @xmath448 = 5.7 , 5.2 and 4.8  gev@xmath449 , for pions , kaons and protons , respectively .
the @xmath448 values were obtained by fitting data  @xcite within the approximate range @xmath450  gev .
@xmath5 was restricted to be less than 2  gev/@xmath65 .
distribution @xmath444 was obtained from the ua5  @xcite measurements for @xmath45 = 200  gev non - singly diffractive p - p collisions for @xmath451 .
the particle species probabilities for pions , kaons and protons at low momentum were assumed to be 0.85 , 0.085 and 0.065 , respectively  @xcite .
particle trajectories were extended outward from the collision point and approximately transverse to the beam direction toward the star central trigger barrel ( ctb ) detector  @xcite which was approximated by a uniform cylinder 220  cm in radius , coaxial with the beam line , and having full @xmath24 azimuthal coverage and total longitudinal length of 484  cm located symmetrically about the geometrical center of the star time projection chamber  @xcite .
charged particles were propagated in a 0.5  t solenoidal magnetic field as helices until they either intersected or missed the ctb .
the simulated mean number of charged particles per unit pseudorapidity at @xmath452 intersecting the ctb was consistent with the measured @xmath50 from ua5  @xcite .
energy deposition in the plastic scintillator material ( ch ) was estimated by sampling the landau distribution  @xcite for the calculated most - probable energy loss  @xcite for 1  cm thick plastic taking into account the range of crossing angles between the helical trajectories and the cylindrical detector . the areal density ( effective thickness ) along the path through the material was estimated to be ( 1.032  gm/@xmath453  @xcite where @xmath113 is the angle of incidence between the particle trajectory helix and a normal to the detector cylinder at the intersection point .
the mean electron excitation ( ionization ) energies for carbon and hydrogen were assumed to be 80  ev and 20  ev , respectively  @xcite .
energy depositions in the carbon and hydrogen components of the detector material were added together .
light attentuation in the plastic scintillator ( 380  cm attenuation length  @xcite ) was included in the simulated ctb output signal . for each p - p collision
the total output from each produced charged particle which deposited energy in the ctb scintillators was summed , resulting in a simulated result for the quantity @xmath133 in sec .
[ sec : discussion ] .
the p - p results in fig .
[ fig4 ] correspond to one - million triggered events with non - zero ctb signal .
the resulting trigger output frequency distribution , @xmath454 , served as input for the au - au trigger simulations described next .    for au - au collisions
the monte carlo glauber model presented here was used to generate an ensemble of @xmath3 = 200  gev simulated minimum - bias collisions . for each au - au collision
the @xmath454 distribution was sampled @xmath455 times corresponding to the average number of p - p collisions required to generate mean multiplicity @xmath16 .
the sampled values of @xmath133 were added to obtain the total ctb trigger output for each simulated au - au collision .
the results in fig .
[ fig4 ] correspond to one - million au - au triggers with non - zero ctb output .
99        glauber r j 1959 _ lectures in theoretical physics _ , edited by brittin w e and dunham l g ( new york : interscience ) 315 glauber r j 1967 _ high energy physics and nuclear structure _ , edited by alexander g ( amsterdam : north - holland ) 311 glauber r j 1970 _ high energy physics and nuclear structure _ , edited by devons s ( new york : plenum ) 207 .
adler c _ et al _ ( star collaboration ) 2001 _ phys .
lett . _ * 87 * 112303 adler c _ et al _ ( star collaboration ) 2002 _ phys .
lett . _ * 89 * 202301 adams j _ et al _ ( star collaboration ) 2006 _ phys .
_ c * 73 * 034906                        ray l , coker w r , hoffmann g w 1978 _ phys .
_ c * 18 * 2641 ray l 1979 _ phys .
1855 ray l , hoffmann g w , barlett m , mcgill j , amann j , adams g , pauletta g , gazzaly m , and blanpied g s 1981 _ phys
_ c * 23 * 828 pauletta g _
et al _ 1981 _ phys .
_ b * 106 * 470 hoffmann g w _ et al _ 1978 _ phys .
_ b * 79 * 376 hoffmann g w _ et al _ 1980 _ phys . rev .
_ c * 21 * 1488                    back b b _ et al _ ( phobos collaboration ) 2001 _ phys .
lett . _ * 87 * 102303 back b b _ et al _ ( phobos collaboration ) 2002 _ phys
_ c * 65 * 061901(r ) back b b _ et al _ ( phobos collaboration ) 2004 _ phys .
_ c * 70 * 021902(r )                                          adams j _ et al _ ( star collaboration ) 2005 _ phys .
_ b * 616 * 8 adams j _ et al _ ( star collaboration ) 2006 _ phys .
_ b * 637 * 161 adler s s _ et al _ ( phenix collaboration ) 2006 _ phys
_ c * 74 * 024904