since the appearance of the first high - energy colliders in the 1980 s , the study of `` jets '' of particles produced in the final state has proved to be a powerful tool for probing the underlying elementary dynamics of the strong force as described by quantum chromodynamics ( qcd ) .
jets have been extensively studied at @xmath2 ( lep , slc ) , @xmath3 ( hera ) and @xmath4/@xmath5 ( tevatron , rhic , lhc ) colliders , with a wide variety of jet algorithms , including infrared and collinear ( irc ) safe algorithms such as those of refs .
@xcite .
a topic of current interest is the use of jets in heavy - ion ( hi ) collisions , where , for example , they can be used to probe the hot , dense medium . in the past few years , the relativistic heavy ion collider ( rhic ) has amassed significant numbers of copper - copper and gold - gold collisions at nucleon - nucleon centre of mass energies up to @xmath6 , and the large hadron collider ( lhc ) should deliver high yields of lead - lead collisions at much higher energies in the near future .
the main obstacle to studying jets in hi collisions is the presence of the huge background given by the underlying event ( ue ) produced simultaneously with the hard nucleon - nucleon collision that initiates the high - transverse momentum jet of interest .
this ue needs to be properly subtracted from the momentum of a given jet in order to reconstruct its `` true '' momentum , _
i.e._the one it would have in the absence of the ue contribution .
this problem is of course well known , and is also present , though to a much smaller extent , in jet studies in proton - proton collisions .
various approaches to address it have been proposed ( see _ e.g._@xcite ) .    in ref .
@xcite two of us proposed a jet area - related technique to determine the transverse momentum density of a sufficiently uniformly distributed background and to subtract it from the jet momenta .
the method of @xcite introduced several novel steps , such as the measurement of jet areas , and procedures to determine the transverse - momentum density of the underlying event and/or pileup .
subsequent work of ours @xcite has sought to provide firmer foundations for these concepts and methods , as well as practical tests in @xmath5 jet reconstruction tasks with simulated events  @xcite .
preliminary experimental jet measurements from the star collaboration at rhic , whose analysis is partially based on the ideas of @xcite , have been presented in ref .
@xcite .    in this article
we give a systematic examination of the performance of such methods for heavy - ion collisions , applying them to monte carlo simulations for rhic and the lhc .
in particular , we determine the accuracy with which jet momenta can be effectively reconstructed , comparing the performance of several different jet algorithms .
jets in hi collisions are produced in an environment that is far from conducive to their detection and accurate measurement .
monte carlo simulations ( and real rhic data ) for gold - gold collisions at @xmath7  gev ( per nucleon - nucleon collision ) show that the transverse momentum density @xmath8 of final - state particles is about 100 gev per unit area ( in the rapidity - azimuth plane ) . for lead - lead collisions at @xmath9  tev at the lhc
this figure is expected to increase by some factor @xmath10 .
this means that jets returned by jet definitions with a radius parameter of , _
e.g. _ , @xmath11 , will contain background contamination of the order of @xmath12 50 and @xmath13 respectively .
a related , and perhaps more challenging obstacle to accurate jet reconstruction is due to the fluctuations both of the background level ( from event to event , but also from point to point in a single event ) and of the jet area itself : knowing the contamination level ` on average ' is not sufficient to accurately reconstruct each individual jet .    a wide range of strategies
have already been advocated to reconstruct jets in hi collisions .
just to mention them briefly , approaches to hi jet reconstruction may involve one or more of the following measures : choosing a jet algorithm returning jets with an area as constant as possible ; eliminating from the clustering all particles below a given transverse momentum threshold , say @xmath14 gev ; measuring the background level in a region of the detector thought to be not affected by the hard event ; or parametrising the average background level ( and fluctuations ) in terms of some other measured properties of the event , such as its centrality .
our framework has the following characteristics :    * we shall restrict ourselves to irc safe jet algorithms , to ensure that measurements can be meaningfully compared to higher - order perturbative qcd calculations . * the jet algorithms we use will
not be limited to those yielding jets of regular shape .
* our analysis will avoid excluding small transverse momentum particles from the clustering .
doing so is collinear - unsafe and inevitably biases the reconstructed jet momenta , which must then be corrected using monte carlo simulations , which can have substantial uncertainties in their modelling of jet quenching and energy loss .
instead , we shall try to achieve a bias - free reconstructed jet , working with all the particles in the event .
cutoffs of their own .
these detector artefacts should not have too large an impact on collinear safety as long as they appear at momenta of the order of the hadronisation scale of qcd , _
i.e._a few hundred mev . ] * the background level will be determined through adaptations of the jet - area / median procedure suggested in @xcite and analysed in detail in  @xcite .
this procedure is designed to give an estimate of the background that is minimally affected by the presence of hard jets .
in order to test the effectiveness of subtracting the underlying event in a heavy - ion collision , and determine the quality of the jet reconstruction , we need access to the idealised hard jets , without the background , as a reference . this can be done by considering a simulated hard @xmath5 event first in isolation , and then embedded in a heavy - ion event . for clarity in the discussions
, we shall adopt the following terminology :    * the _ hard event _ refers to the hard @xmath5 event alone , without a heavy - ion background ; * the _ full event _ refers to combination of the hard event and the hi background ( which possibly includes many semi - hard events of its own ) ; * _ hard jets _ and _ full jets
_ refers to jets from the hard and full events , respectively .
note that though making the distinction between a hard and a full event is not possible with real data , it is feasible in certain monte carlo studies .
one should nevertheless be aware that the extent to which such a distinction is physically meaningful is a question that remains open in view of the numerous issues related to the interaction between an energetic parton and the medium through which it travels .    if the monte carlo program used to simulate heavy - ion events explicitly provides a separation between hard events and a soft part
_ e.g._as does hydjet @xcite  one can extract one of the former as the single hard event and take the complete event as the full one .
alternatively , one can generate a hard event independently and embed it in a heavy - ion event ( obtained from a monte carlo or from real collisions ) to obtain the full event .
this second approach , which we have adopted for the bulk of results presented here because of its greater computational efficiency , is sensible as long as the embedded hard event is much harder than any of the semi - hard events that tend to be present in the background .
for the transverse momenta that we consider , this condition is typically fulfilled .
note however that for studies in which the presence of a hard collision is not guaranteed , _
e.g._the evaluation of fake - jet rates , one should use the first approach .      to measure the performance of our heavy - ion background subtraction procedure
we apply it to both hard and full events .
we then take the two hardest resulting jets in each hard event and compare them to the corresponding subtracted jets in the full event .
to do this in practice , we need a method to match the jets obtained from the full event to those from the hard event alone , to make sure that we are comparing the same object .
in other analyses , this matching was typically performed in terms of the position of the jet in the rapidity - azimuth plane , requiring that the two jet axes not differ by more than a given @xmath15 .
the exact value of @xmath16 is an arbitrary choice , and values between 0.1 and 0.3 have been often used . in order to avoid the arbitrariness of the @xmath16 choice
, we propose a different prescription : a jet reconstructed in the full event is considered matched to a hard jet _ if the constituents common to both the hard and the full jet make up at least 50% of the transverse momentum of the constituents of the hard jet_. of the @xmath1 of the hard jet after ue subtraction as in section  [ sec : subtraction ] . in practice , this detail has negligible impact on matching efficiencies and other results . ]
this definition has the advantage that for a given hard jet , at most one full jet can satisfy this criterion , therefore avoiding having to deal with multiple positive matchings .
another advantage is that it automatically rejects fake matchings given by a soft jet that happens to be close to a hard one .    for a matched pair of jets
, we shall use the notation @xmath17  and @xmath18  for the transverse momentum of the jet in the hard and full event respectively , and @xmath19  and @xmath20  for their subtracted equivalents .
the quantity that we shall mostly concentrate on is the @xmath1 offset , @xmath21 between the momentum of the background - subtracted jet in the full event and its equivalent in the hard event .
large deviations from zero will indicate a poor subtraction .
events without any matched jets are simply discarded and instead contribute to the evaluation of matching inefficiencies ( see section  [ sec : efficiency ] ) .
the average of the @xmath1 shift over many events , @xmath22 , is only one measure that one may examine to establish the quality of jet reconstruction ; its dispersion , @xmath23 can be another important one .
this is especially true in the case of steeply - falling jet spectra , where this dispersion can have a large impact on the measured cross section , necessitating delicate deconvolutions , also known as `` unfolding '' .    in this paper
we shall concentrate on these two quality measures @xmath24 and @xmath25 , keeping only the pair of jets that have been matched to one of the two hardest ( subtracted ) jets in the hard event .
small ( absolute ) values of both @xmath24 and @xmath25 will be the sign of a good subtraction . in practice a trade - off may exist in offset versus dispersion : which one to optimise may depend on the specific observable one wants to measure . and @xmath25 may not , in general , fully characterise the quality of the reconstruction , as the distribution of @xmath26 may be non - gaussian . ]
quantifying the quality of background subtraction using monte carlo simulations has several advantages . besides providing a practical way of generating the hard `` signal '' separately from the soft background
, one can easily check the robustness of one s conclusions by changing the hard jet or the background sample .
one difficulty that arises in gauging the quality of jet reconstruction in heavy - ion collisions comes from the expectation that parton fragmentation in a hot medium will differ from that in a vacuum .
this difference is often referred to as _ jet quenching _
@xcite ( for reviews , see  @xcite ) .
the details of jet quenching are far less well established than those of vacuum fragmentation and can have an effect on the quality of jet reconstruction . here
we shall examine the reconstruction of both unquenched and quenched jets .
for the latter it will be particularly important to be able to test more than one quenching model , in order to help build confidence in our conclusions about any reconstruction bias that may additionally exist in the presence of quenching .    in practice , for this paper we have used both the fortran ( v1.6 ) @xcite and the c++ ( v2.1 ) @xcite versions of hydjet to generate the background .
hard jets have been generated with pythia  6.4 @xcite , either running it standalone , or using the version embedded in hydjet v1.6 .
the quenching effects have been studied using both qpythia @xcite and pyquen @xcite .      as mentioned in the introduction
, the last few years have seen many developments in the field of jet clustering , with the appearance of fast implementations  @xcite of the @xmath0  @xcite and the cambridge / aachen  @xcite sequential recombination algorithms , and the introduction of the siscone  @xcite and anti-@xmath0  @xcite algorithms .    besides these four main infrared - and - collinear - safe jet algorithms , a number of recent papers ( see _ e.g._@xcite ) ,
have introduced and studied jet - cleaning techniques in order to help the reconstruction of jets by reducing the ue contamination .
`` filtering '' @xcite works by reclustering each jet with a radius @xmath27 smaller than the original radius @xmath28 and keeping only the @xmath29 hardest subjets ( background subtraction is applied to each of the subjets before deciding which ones to keep ) .
other similar techniques , `` trimming '' and `` pruning '' , that exploit the substructure of jets have also been introduced recently @xcite and , in @xmath5 environments , have been found to give benefits comparable to those of filtering .
note that the filtering that we use here is unrelated to the gaussian filter approach of @xcite ( not used in this paper as no public code is currently available ) .
in analogy with the systematic analysis of the performance of various jet definitions in the kinematic reconstruction of dijet systems at the lhc @xcite , here we will study how different algorithms behave in the case of hi collisions .
we will use the @xmath0 , cambridge / aachen ( c / a ) , and anti-@xmath0 algorithms , as well as filtering applied to cambridge / aachen ( c / a(filt ) ) with @xmath30 and @xmath31 .
we have not used siscone , as its relatively slower speed compared to the sequential recombination algorithms makes it less suitable for a hi environment .
all the algorithms have been used through their fast implementation available in the package @xcite , version  2.4.2 .
additionally some features of a forthcoming fastjet release have been used to help simplify our analyses .      , the _ circular _ range @xmath32 and the _ doughnut _ range @xmath33 .
the last three are local ranges with a position depending on the jet being subtracted .
see the text for detailed definitions.,scaledwidth=80.0% ]    in order to subtract the hi background from the hard jets , we will mainly follow the method introduced in @xcite .
when clustering the event we determine the 4-vector active area @xmath34 @xcite of each jet @xmath35 , as well as an estimate of the background density @xmath8 .
then for each jet @xmath35 , we subtract from its four - momentum the expected background contamination : @xmath36 the background transverse - momentum density per unit area , @xmath8 , is determined , event - by - event , as proposed in @xcite . in that paper two
ways of determining @xmath8 were explored .
one was to take a global background - estimation range @xmath37 , covering the full rapidity - azimuth plane up to some rapidity @xmath38 ; one then considered the set of ratios of transverse components of the momentum and area 4-vectors , @xmath39 , for the jets within this range .
the median of this set was used as an estimate for @xmath8 : @xmath40 where the subscript @xmath37 explicitly denotes that the background density has been calculated using only jets in the range @xmath37 .
this method assumes that the background level is sufficiently constant within @xmath37 .
this condition is known to be violated when going to sufficiently large rapidity , and the effect is particularly marked in heavy ion collisions .
for this reason , it was alternatively proposed in @xcite to fit the mean @xmath39 as a function of rapidity @xmath41 with a quadratic functional form @xmath42 .    while working on this paper , and partly inspired by discussions with members of the star collaboration ( c.f .
also ref .
@xcite ) , we came to the conclusion that neither of these two methods allows one to extract sufficiently accurate values of @xmath8 in the context of hi collisions .
we propose therefore a variant of the first method , namely making the background - estimation range _ local _ and dependent on a given jet s position .
graphical representations of a global range , as well as of three possible local ranges that we shall use hereafter , are given in fig .
[ fig : ranges ] .
more specifically , for a given jet @xmath35 , our three local ranges are defined as follows ,    * the _ strip _ range , @xmath43 , includes the jets
@xmath44 satisfying @xmath45 , * the _ circular _ range , @xmath46 , includes the jets
@xmath44 satisfying @xmath47 , * the _ doughnut _ range , @xmath33 , includes the jets @xmath44 satisfying @xmath48 .    using @xmath49 to denote any local range around jet @xmath35 ,
the background density will then depend on the jet being subtracted and will be estimated using @xmath50 it is important to be aware that the estimate of the background density @xmath8 is just an input to eq .
( [ eq : subtraction ] ) and that the jet definition that is used for the computation of @xmath8 in eqs .
( [ eq : median ] ) or ( [ eq : local_median ] ) can be different from the jet definition used to obtain the `` physical '' jets and subtract them in eq .
( [ eq : subtraction ] ) .
the only two recommended jet algorithms for the task of determining @xmath8 are the @xmath0 and c / a algorithms ( see @xcite ) ; others are not ideal because they return many jets with very small areas , which distorts the median procedure .
the radius , @xmath51 , used in the jet definition for determining @xmath8 can also differ from that , @xmath28 , used to find the jets .
when using a local range , the underlying idea is to limit the sensitivity to the long - range variations of the background density , by using only the jets in the vicinity of the jet we want to subtract . in practice ,
a compromise needs to be found between choosing a range small enough to get a valid local estimation , but also large enough to contain a sufficiently large number of background ( soft ) jets for the estimation of the median to be reliable .
two effects need to be considered : ( 1 ) statistical fluctuations in the estimation of the background and ( 2 ) biases due to the presence of a hard jet in the region used to estimate the background .    1 .
if we require that the dispersion in the reconstructed jet @xmath1 coming from the statistical fluctuations in the estimation of the background does not amount to more than a fraction @xmath52 of the overall @xmath25 , then as discussed in appendix  [ app : minrange - fluct ] the range needs to cover an area @xmath53 such that @xmath54 taking @xmath55 and @xmath11 this corresponds to an area @xmath56 . for the applications below ,
when clustering with radius @xmath28 , we have used the rapidity - strip ranges @xmath57 and @xmath58 , the circular range @xmath59 and the doughnut range @xmath60 .
one can check that their areas are compatible with the @xmath61 lower - limit estimated above .
the bias in the estimate of @xmath8 due to the presence of @xmath62 hard jets in the range is given roughly by @xmath63 as discussed in @xcite and appendix  [ app : minrange - bias ] .
the ensuing bias on the @xmath1 can be estimated as @xmath64 ( for anti-@xmath0 jets ) . for @xmath11 , @xmath65 , @xmath66 and @xmath67 , the bias in the reconstructed jet @xmath1
is @xmath68 . given @xmath69 in the @xmath70 range ( as we will find in section  [ sec : results ] ) this corresponds to a @xmath71 bias . in order to eliminate this small bias , we will often choose to exclude the two hardest jets in each event when determining of @xmath8 .    a third potential bias discussed in ref .
@xcite is that of underestimating the background when using too small a value for @xmath51 .
given the high density of particles in hi collisions , this will generally not be an issue as long as @xmath72 .
distribution of the background density @xmath8 per unit area ( left ) and its intra - event fluctuations @xmath69 ( right ) .
it has been obtained from 5000 hydjet events with rhic ( auau , @xmath73 ) and lhc ( pbpb , @xmath74 ) kinematics .
the background properties have been estimated using the techniques presented in section [ sec : subtraction ] , using the @xmath0 algorithm with @xmath75 , and keeping only the jets with @xmath76 ( excluding the two hardest).,title="fig:",scaledwidth=50.0% ]   distribution of the background density @xmath8 per unit area ( left ) and its intra - event fluctuations @xmath69 ( right ) .
it has been obtained from 5000 hydjet events with rhic ( auau , @xmath73 ) and lhc ( pbpb , @xmath74 ) kinematics .
the background properties have been estimated using the techniques presented in section [ sec : subtraction ] , using the @xmath0 algorithm with @xmath75 , and keeping only the jets with @xmath76 ( excluding the two hardest).,title="fig:",scaledwidth=50.0% ]    as mentioned already in section  [ sec : mcs ] , for the simulation of the events , we have used pythia 6.4 @xcite to generate the unquenched hard jets ; the background events are generated using hydjet v1.6 @xcite with 0 - 10% centrality , @xmath77 , @xmath78 and ` ptmin`@xmath79 ; for lhc , @xmath80 , @xmath81 , @xmath82 and @xmath83 . in both cases quenching effects are turned on in hydjet , @xmath84 , even when they are not included for the embedded @xmath5 event .
the corresponding pyquen parameters we have used are @xmath85 , @xmath86 , @xmath87 , @xmath88  fm and @xmath89 .
[ foot : pyquen ] ] .
our setup for the background leads , for rhic ( auau , @xmath90 ) , to an average background density per unit area at central rapidity of @xmath91  gev with average fluctuations in a single event of @xmath92  gev and event - to - event fluctuations @xmath93  gev . for the lhc ( pbpb ,
@xmath74 ) the corresponding values are @xmath94  gev , @xmath95  gev and event - to - event fluctuations @xmath96  gev .
[ fig : bkg_distrib ] shows the distributions obtained from the simulations for @xmath8 and @xmath69 .    in the case of the rhic simulation ,
the result for @xmath97 is somewhat higher than the experimental value of @xmath98 quoted by the star collaboration @xcite , however this is probably in part due to limited tracking efficiencies at low @xmath1 at star , and explicit star  @xcite and phenix  @xcite results for @xmath99 correspond to somewhat higher @xmath8 values , about @xmath100 .
the multiplicity of charged particles ( @xmath101 for @xmath102 and 06% centrality ) and the pion @xmath1 spectrum in our simulation are sensible compared to experimental measurements at rhic @xcite . for the lhc our charged particle multiplicity is @xmath103 for @xmath102 and 010% centrality , which is comparable to many of the predictions reviewed in fig .  7 of @xcite .    an independent control analysis has also been performed with hydjet++  2.1  @xcite ( with default parameters ) for the background .
the results at rhic are similar , while for lhc the comparison is difficult because the default tune of hydjet++  2.1 predicts a much higher multiplicity , @xmath104 for @xmath102 and 010% centrality .
most of the results of this section will be obtained without quenching , though in section  [ sec : quenching ] we will also consider the impact on our conclusions of the pyquen  @xcite and qpythia  @xcite simulations of quenching effects .    for the results presented below , we have employed a selection cut of @xmath105 on the jets with @xmath106 for rhic and @xmath107 for the lhc . in that case , in order for information to be derived from jets near the edge of the acceptance and thus bring the available area close to the ideal requirements set out in appendix  [ sec : minimal - range ] .
note that we also use particles beyond @xmath38 in the jet clustering and apply the acceptance cut only to the resulting jets .
] we only consider full jets that are matched to one of the two hardest jets in the hard event .
the computation of the jet areas in , needed both for the subtraction and the background estimation , has been performed using active areas , with ghosts up to @xmath108 , a single repetition and a ghost area of 0.01 .
the determination of the background density @xmath8 has been performed using the @xmath0 algorithm with @xmath75 . though the estimate of the background depends on @xmath109 @xcite ,
we have observed that choices between 0.3 and 0.5 lead to very similar results ( _ e.g._differing by at most a few hundred mev at rhic ) .    with this setup , we have studied the various ranges presented in section [ sec : subtraction ] , with the jet algorithms from section [ sec : jetdef ] . in all cases
, we have taken the radius parameter @xmath11 .
we have adopted this value as it is the largest currently used at rhic .
note that the effect of the background fluctuations of the jet energy resolution increases linearly with @xmath28 , disfavouring significantly larger choices . on the other hand ,
too small a choice of @xmath28 may lead to excessive sensitivity to the details of parton fragmentation , hadronisation and detector granularity .       matching efficiency for reconstructed jets as a function of the jet @xmath1 .
left : rhic , right : lhc .
these results are independent of the choice of background - subtraction range in the heavy - ion events , since background subtraction does not enter into the matching criterion . here and in later figures
, the label `` unquenched '' refers to the embedded @xmath5 event ; the background is always simulated including quenching .
, title="fig:",scaledwidth=50.0% ]   matching efficiency for reconstructed jets as a function of the jet @xmath1 . left : rhic , right : lhc .
these results are independent of the choice of background - subtraction range in the heavy - ion events , since background subtraction does not enter into the matching criterion . here and in later figures , the label `` unquenched '' refers to the embedded @xmath5 event ; the background is always simulated including quenching .
, title="fig:",scaledwidth=50.0% ]    let us start the presentation of our results with a brief discussion of the efficiency of reconstructing jets in the medium .
as explained in section [ sec : quality ] , the jets in the medium are matched to a `` bare '' hard jet when their common particle content accounts for at least 50% of the latter s transverse momentum .
the matching efficiencies we observe depend to some extent on the details of the monte carlo used for the background so our intention is just to illustrate the typical behaviour we observe and highlight that these efficiencies tend to be large .
we observe from fig .
[ fig : efficiency ] that we successfully match at least 95% of the jets above @xmath110 gev at rhic , and at least 99% of the jets above @xmath111 gev at the lhc .
it is also interesting to notice that the anti-@xmath0 algorithm performs best , likely as a consequence of its ` rigidity ' , namely the fact that anti-@xmath0 jets tend to have the same ( circular ) shape , independently of the soft - particles that are present .
effect of the choice of range on the average @xmath1 shift , @xmath26 , as defined in eq .
( [ eq : deltapt ] ) . left : rhic , right : lhc . in this figure and those that follow
, the yellow band corresponds to 1% of the @xmath1 of the hard jet.,title="fig:",scaledwidth=50.0% ]   effect of the choice of range on the average @xmath1 shift , @xmath26 , as defined in eq .
( [ eq : deltapt ] ) . left : rhic , right : lhc . in this figure and those that follow , the yellow band corresponds to 1% of the @xmath1 of the hard jet.,title="fig:",scaledwidth=50.0% ]    we now turn to the results concerning the measurement of the background density and the reconstruction of the jet transverse momentum .
we first concentrate on the impact of the choice of a local range and/or of the exclusion of the two hardest jets when determining @xmath8 . for the determination of @xmath8 in the hard event , without exclusion of any jets .
this ensures that the reference jet @xmath1 is always kept the same .
the impact of subtraction in the hard event is in any case small , so the particular choice of range is not critical . ]    in fig .
[ fig : ptshift_ranges ] we show the average shift @xmath24 for the list of ranges mentioned in section [ sec : subtraction ] .
the results presented here have been obtained with the anti-@xmath0 algorithm with @xmath11 , but the differences among the various range choices have been seen to be similar with other jet definitions .
the label `` 2 excl '' means that the two hardest jets in the event have been excluded from the estimation of the background .
we have found that this improves the precision of the subtractions whenever expected , _ i.e._for all choices of range except the doughnut range , where its central hole already acts similarly to the exclusion of the hardest jets . to keep the figure reasonably readable ,
we have only explicitly shown the effect of removing the two hardest jets for the global range .
the change of 0.40.6  gev ( both for rhic and the lhc ) is in reasonable agreement with the analytic estimate of about 0.6 gev for rhic and the lhc obtained from @xmath112 with @xmath113 calculated using eq .
( [ eq : hard_median_offset ] ) .
note that at lhc the exclusion of the two hardest jets for the global range appears to worsen the subtraction , however what is really happening is that the removal of the two hardest jets exacerbates a deficiency of the global range , namely the fact that its broad rapidity coverage causes it to underestimate @xmath8 , leading to a positive net @xmath113 .
other features that can be understood qualitatively include for example the differences between the two strip and the global ( 2  excl ) range for rhic : while the rapidity width of the global range lies in between that of the two strip ranges , the global range gives a lower @xmath24 than both , corresponding to a larger @xmath8 estimate , which is reasonable because the global range is centred on @xmath114 , whereas the strip ranges are mostly centred at larger rapidities where the background is lower .
the main result of the analysis of fig .
[ fig : ptshift_ranges ] is the observation that all choices of a local range lead to a small residual @xmath26 offset : the background subtraction typically leaves a @xmath115 at both rhic and lhc , _
i.e._better than 1 - 2% accuracy over much of the @xmath1 range of interest .
it is not clear , within this level of accuracy , if one range is to be preferred to another , nor is it always easy to identify the precise origins of the observed differences between various ranges .
another way of viewing this is that the observed differences between the various choices give an estimate of the residual subtraction error due to possible misestimation of @xmath8 . for our particular analysis , at rhic
this comment also applies to the choice of the global range ( with the exclusion of the two hardest jets in the event ) .
this is a consequence of the limited rapidity acceptance , which effectively turns the global range into a local one , a situation that does not hold for larger rapidity acceptances , as we have seen for the lhc results . in what follows
we will use the doughnut(@xmath116 ) choice , since it provides a good compromise between simplicity and effectiveness .      the next potential systematic effect that we consider
is the choice of the jet algorithm used for the clustering .
algorithm with @xmath75 is used for the estimation of the background .
[ fig : deltapt - dist ] shows the distribution of @xmath26 for each of our four choices of jet algorithm , @xmath0 , c / a , anti-@xmath0 and c / a(filt ) , given for rhic collisions and a specific bin of the hard jets transverse momenta , @xmath117 .
one sees significant differences between the different algorithms .
one also observes that gaussians with mean and dispersion set equal to @xmath24 and @xmath25 provide a fair description of the full histograms .
this validates our decision to concentrate on @xmath24 and @xmath25 as quality measures .
one should nevertheless be aware that in the region of high @xmath118 there are deviations from perfect gaussianity , which are more visible if one replicates fig .
[ fig : deltapt - dist ] with a logarithmic vertical scale ( not shown , for brevity ) .
average shift @xmath24 , as a function of @xmath119 , shown for rhic ( left ) and the lhc ( right).,title="fig:",scaledwidth=50.0% ]   average shift @xmath24 , as a function of @xmath119 , shown for rhic ( left ) and the lhc ( right).,title="fig:",scaledwidth=50.0% ]    the first observable we analyse is the average @xmath1 shift .
we show in fig .
[ fig : ptshift ] the @xmath120 results for the four algorithms listed in section [ sec : jetdef ] , as a function of @xmath119 .
we use the doughnut range to estimate the background .
the first observation is that , while the anti-@xmath0 and c / a(filt ) algorithms have a small residual @xmath24 , the c / a and @xmath0 algorithms display significant offsets .
the reason for the large offsets of @xmath0 and c / a is well understood , related to an effect known as _ back - reaction
_  @xcite .
this is the fact that the addition of a soft background can alter the clustering of the particles of the hard event : some of the constituents of a jet in the hard event can be gained by or lost from the jet when clustering the event with the additional background of the full event .
this happens , of course , on top of the simple background contamination that adds background particles to the hard jet .
even if this latter contamination is subtracted exactly , the reconstructed @xmath1 will still differ from that of the original hard jet as a consequence of the back - reaction .
contribution to @xmath24 due to back - reaction .
note that these results are independent of the range used for estimating @xmath8 in the heavy - ion event . left : rhic , right : lhc.,title="fig:",scaledwidth=50.0% ]   contribution to @xmath24 due to back - reaction .
note that these results are independent of the range used for estimating @xmath8 in the heavy - ion event . left : rhic , right : lhc.,title="fig:",scaledwidth=50.0% ]    the effect of the back - reaction can be studied in detail , since in monte carlo simulations it is possible to identify which hard - event constituents are present in a given jet before and after inclusion of the background particles in the clustering . the average @xmath1 shift due to back - reaction can be seen in fig .  [ fig : ptshift_br ] for the different jet algorithms .
as expected @xcite , it is largest for @xmath0 , and smallest ( almost zero , in fact ) for anti-@xmath0 . by comparing fig .
[ fig : ptshift_br ] and fig .
[ fig : ptshift ] one can readily explain the difference between the @xmath24 offsets of the various algorithms in terms of their back - reaction .
the rigidity ( and hence small back - reaction ) of the anti-@xmath0 jets manifestly gives almost bias - free reconstructed jets , while the large back - reaction effects of the @xmath0 algorithm and , to a smaller extent , of the c / a algorithm translates into a worse performance in terms of average shift .
the @xmath1 dependence of the back - reaction is weak .
this is expected based on the interplay between the @xmath121 dependence found in @xcite and the evolution with @xmath1 of the relative fractions of quark and gluon jets .
the case of the c / a(filt ) algorithm is more complex : its small net offset , comparable to that of the anti-@xmath0 algorithm , appears to be due to a fortuitous compensation between an under - subtraction of the background and a negative back - reaction .
the negative back - reaction is very similar to that of c / a without filtering , while the under - subtraction is related to the fact that the selection of the hardest subjets introduces a bias towards positive fluctuations of the background .
this effect is discussed in appendix  [ app : filtbias ] , where we obtain the following estimate for the average @xmath1 shift ( specifically for @xmath122 ) : @xmath123 yielding an average bias of 2  gev for rhic and 4.5  gev at the lhc , which are both in good agreement with the differences observed between c / a with and without filtering in fig .
[ fig : ptshift ] . note that while the bias in eq .
( [ eq : filt - bias - main - text ] ) is proportional to @xmath69 , the back - reaction bias is instead mainly proportional to @xmath8  @xcite .
it is because of these different proportionalities that the cancellation between the two effects should be considered as fortuitous .
since it also depends on the substructure of the jet , one may also expect that the cancellation that we see here could break down in the presence of quenching .
dispersion @xmath25 . left : rhic , right : lhc.,title="fig:",scaledwidth=50.0% ]   dispersion @xmath25 . left : rhic , right : lhc.,title="fig:",scaledwidth=50.0% ]    our results for the @xmath26 dispersion , @xmath25 , are shown in fig .
[ fig : dispersion ] , again using the doughnut range .
( our conclusions are essentially independent of the particular choice of range . )
we first discuss the case of rhic kinematics . for @xmath0 and anti-@xmath0 ,
the observed dispersions are similar to the result of @xmath124 quoted by star  @xcite ( though the number from star includes detector resolution effects , so that the true physical @xmath25 may actually be somewhat lower ) . of note , the advantage enjoyed by anti-@xmath0 in terms of smallest @xmath24 does not hold at the level of the dispersion : c / a and @xmath0 tend to behave slightly better at small transverse momentum .
the algorithm which performs best in terms of dispersion over all the @xmath1 range is now c / a with filtering , for which the result is smaller than that of the other algorithms by a factor of about @xmath125 .
this reduction factor can be explained because the dispersion @xmath25 is expected to be proportional to the square - root of the jet area : the c / a(filt ) algorithm with @xmath122 and @xmath31 produces jets with an area of , roughly , half that obtained with c / a ; hence the observed reduction of @xmath25 .    in the lhc setup ,
the conclusions are quite similar at the lowest transverse momenta shown . as @xmath1 increases
, the dispersion of the anti-@xmath0 algorithm grows slowly , while that of the others grows more rapidly , so that at the highest @xmath1 s shown , the @xmath0 and c / a algorithms have noticeably larger dispersions than anti-@xmath0 , and c / a(filt ) becomes similar to anti-@xmath0 .
the growth of the dispersions can be attributed to an increase of the back - reaction dispersion .
the latter is dominated by rare occurrences , where a large fraction of the jet s @xmath1 is gained or lost to back - reaction , hence the noticeable @xmath1 dependence ( _ cf._appendix  [ app : dispersion - contributions ] ) .
an additional effect , especially for the @xmath0 algorithm , might come from the anomalous dimension of the jet areas , _
i.e._the growth with @xmath1 of the average jet area .
note that the dispersion has some limited dependence on the choice of ghost area  for example , reducing it from 0.01 to 0.0025 lowers the dispersions by about @xmath126 at rhic .
this is discussed further in appendix  [ app : ghost - area - choice ] .
@xmath1 dependence dependence of the @xmath26 dispersion at rhic , for different centrality classes.,scaledwidth=60.0% ]    so far , we have only considered central collisions . since it is known that non - central collisions give rise to elliptic flow @xcite
, one might worry that this leads to an extra source of background fluctuations and/or non - uniformities , potentially spoiling the subtraction picture discussed so far .
one can study this on azimuthally averaged jet samples ( as we have been doing so far ) or as a function of the azimuthal angle , @xmath127 , between the jet and the reaction plane . as above , we use hydjet  v1.6 , whose underlying hydro component includes a simulation of elliptic flow  @xcite .    we have generated heavy - ion background events for rhic in four different centrality bins : 0 - 10% ( as above ) , 10 - 20% , 20 - 40% and 40 - 75% , with @xmath128 values respectively of @xmath129 , @xmath130 , @xmath131 and @xmath132 .
was determined as the average of @xmath133 over all all particles with @xmath134 ( excluding the additional hard @xmath5 event ) . ]
we first examine azimuthally averaged results , repeating the studies of the previous sections for each of the centrality bins .
we find that the results for the average shift , @xmath24 , are largely independent of centrality , as expected if the elliptic flow effects disappear when averaged over @xmath135 .
the results for the dispersion are shown in fig .
[ fig : dispersion_centrality ] .
we observe that the dispersion decreases with increasing non - centrality .
even though one might expect adverse effects from elliptic flow , the heavy - ion background decreases rapidly when one moves from central to peripheral collisions , and this directly translates into a decrease of @xmath25 .
the first conclusion from this centrality - dependence study is therefore that the subtraction methods presented in this paper appear to be applicable also for azimuthally averaged observables in non - central collisions .
@xmath135 dependence of the @xmath26 shift at rhic , for the @xmath136 centrality bin .
left : for three different ranges for the anti-@xmath0 algorithm ; right : for three different jet algorithms for the doughnut(@xmath116 ) range .
the absolute size of the @xmath135 dependence is similar for centralities up to @xmath137 and then decreases beyond .
, title="fig : " ]   @xmath135 dependence of the @xmath26 shift at rhic , for the @xmath136 centrality bin .
left : for three different ranges for the anti-@xmath0 algorithm ; right : for three different jet algorithms for the doughnut(@xmath116 ) range .
the absolute size of the @xmath135 dependence is similar for centralities up to @xmath137 and then decreases beyond .
, title="fig : " ]    we next consider results as a function of @xmath138 , which is relevant if one wishes to examine the correlation between jet quenching and the reaction plane .
an issue in real experimental studies is the determination of the reaction plane , and the extent to which it is affected by the presence of hard jets . in the hydjet simulations ,
this problem does not arise because the reaction plane always corresponds to @xmath139 .
figure  [ fig : ptshift - v - deltaphi ] ( left ) shows the average @xmath26 as a function of @xmath138 for the anti-@xmath0 algorithm and several different background - estimation ranges , for the @xmath136 centrality bin for rhic .
the strip range shows significant @xmath138 dependence , which is because a determination of @xmath8 averaged over all @xmath135 can not possibly account for the local @xmath135-dependence induced by the elliptic flow .
other ranges , such as the doughnut range , instead cover a more limited region in @xmath135 .
they should therefore be able to provide information on the @xmath135-dependence of the background .
however , since their extent in @xmath135 tends to be significantly larger than that of the jet , and @xmath8 varies relevantly over that extent , some residual @xmath135 dependence remains in @xmath24 after subtraction .
the right - hand plot of figure  [ fig : ptshift - v - deltaphi ] shows that the effect is reduced with filtering , as is to be expected since its initial background contamination is smaller .
the conclusion from this part of the study is that residual @xmath135-dependent offsets may need to be corrected for explicitly in any studies of jets and their correlations with the reaction plane .
the investigation of extensions to our background subtraction procedure to address this issue will be the subject of future work .
the last issue we wish to investigate is how the phenomenon of jet quenching ( _ i.e._medium effects on parton fragmentation ) may affect the picture developed so far .
the precise nature of jet quenching beyond its basic analytic properties ( see _ e.g._@xcite ) is certainly hard to estimate in detail , especially at the lhc , where experimental data from , say , flow or particle spectra measurements are not yet available for the tuning of the monte carlo simulations .
additionally , the implementation of monte - carlo generators that incorporate the analytic features of jet quenching models is currently a very active field @xcite .
we may therefore expect a more robust and complete picture of jet quenching in the near future , together with the awaited first pbpbcollisions at the lhc .    in this section ,
we examine the robustness of our hi background subtraction in the presence of ( simulated ) jet quenching . for this purpose
we have used two available models which allow one to simulate quenched hard jets , pyquen @xcite , which is used by hydjet  v1.6 , and qpythia @xcite .
pyquen has been run with the parameters listed in footnote [ foot : pyquen ] for the lhc , and with @xmath140 , @xmath141  fm and @xmath142 for rhic .
for qpythia we have tested two options for the values of the transport coefficient and the medium length ( @xmath143  gev@xmath144/fm , @xmath145  fm and @xmath146  gev@xmath144/fm , @xmath147  fm ) , with similar results .
no serious attempt has been made to tune the two codes with each other or with the experimental data , beyond what is already suggested by the code defaults : in the absence of strong experimental constraints on the details of the quenching effects , this allows us to verify the robustness of our results for a range of conditions .    as in sections
[ sec : choice_range ] and [ sec : choice_alg ] , we have embedded the hard pyquen or qpythia events in a hydjet  v1.6 background and tested the effectiveness of the background subtraction for different choices of algorithm . for the unquenched case .
we will therefore keep employing the doughnut range .
] we shall restrict our attention to the anti-@xmath0 and c / a(filt ) algorithms , as they appear to be the optimal choices from our analysis so far .
average @xmath1 shift for background - subtracted jets with the anti-@xmath0 and c / a(filt ) jet algorithms .
the dashed lines correspond to unquenched hard jets ( pythia ) and the solid ones to quenched hard jets ( pyquen ) .
results are shown for rhic kinematics on the left plot and for the lhc on the right one.,title="fig:",scaledwidth=50.0% ]   average @xmath1 shift for background - subtracted jets with the anti-@xmath0 and c / a(filt ) jet algorithms .
the dashed lines correspond to unquenched hard jets ( pythia ) and the solid ones to quenched hard jets ( pyquen ) .
results are shown for rhic kinematics on the left plot and for the lhc on the right one.,title="fig:",scaledwidth=50.0% ]    we have found that the jet - matching efficiencies are still high , with essentially no changes at rhic , and at lhc a doubling of the ( small ) inefficiencies that we saw in fig .
[ fig : efficiency](right ) .
the dispersions @xmath25 are also not significantly affected within our sample of jet - quenching simulations .
we therefore concentrate on the @xmath24 offset , which is plotted in fig .
[ fig : ptshift_quenched ] for pyquen .
the results are shown for both rhic and the lhc . in the case of rhic , and for the whole @xmath1 range
up to about 50 gev , quenching can be seen not to significantly affect the subtraction offset @xmath148 ( within the usual uncertainty related to the choice of range , which was shown in fig .
[ fig : ptshift_ranges ] ) . in the lhc case , instead
, while the shift obtained using the anti-@xmath0 algorithm is largely similar to the unquenched case , the c / a(filt ) algorithm performance can be seen to deteriorate slightly when quenching is turned on , all the more so at very large transverse momentum .
the @xmath0 and the c / a algorithms are not shown for clarity , but they share the behaviour of c / a(filt ) .
this deterioration of the quality of the subtraction can be traced back to an increased back - reaction compared to the unquenched jets .
anti-@xmath0 jets do not suffer from this effect as a consequence of the usual rigidity of this algorithm . in the case of c / a(filt )
, one should nevertheless emphasise that an error of ( at most ) 10 gev on the reconstruction of a 500 gev jet is still only a 2% effect .
this is modest , both relative to the likely experimental precision and to the expected effect of quenching on the overall jet @xmath1 , predicted by pyquen to be at the level of @xmath149 at this @xmath1 .
though for brevity we have not explicitly shown them , the results with qpythia are very similar .    before closing this section , we reiterate that we have only investigated simple models for quenching and that our results are meant just to give a first estimate of the effects that one might have to deal with in the case of quenched jets .
the expected future availability of new `` quenched '' monte carlo programs , together with specific measurements in the early days of heavy - ion collisions at the lhc , will certainly allow one to address this question more extensively .      to close this section , we examine the relative importance of the average shift and its dispersion , taking the illustrative example of their impact on the inclusive jet cross - section as a function of @xmath1 .
we start from a simple parametrisation of the inclusive - jet @xmath1 spectrum and see how its reconstruction is affected by the average shift and dispersion that remain for @xmath26 after subtraction .
let us assume that the true @xmath1 spectrum decays exponentially _
i.e._@xmath150 while this expression does nt have the @xmath151 form that one expects to see , it is far easier to handle analytically , and not too poor an approximation to observed spectra over quite a broad range of @xmath1 .
after embedding the hard events in a heavy - ion background and applying subtraction , the resulting spectrum will be the convolution of eq .
( [ eq : hard - spectrum ] ) with the @xmath26 distribution .
assuming that the latter is a gaussian of average @xmath24 and dispersion @xmath25 ( _ cf._fig .
[ fig : deltapt - dist ] ) , one obtains a reconstructed spectrum @xmath152 thus the average shift gives a bias by a multiplicative factor @xmath153 and the dispersion by a factor @xmath154 .
the convolution works in such a way that for a given reconstructed @xmath1 , the most likely original true transverse momentum is : @xmath155 where we have neglected the small impact of subtraction on the @xmath5 jets .    to illustrate these effects quantitatively ,
let us first take the example of rhic , where between 10 and 60 gev , the cross - section is well approximated by eq .
( [ eq : hard - spectrum ] ) with @xmath156 gev@xmath157 . both the anti-@xmath0 and c / a(filt ) have @xmath158 , leaving only the dispersion effect . in the case of the anti-@xmath0 ( respectively c / a(filt ) )
algorithm , we see from fig .  [
fig : dispersion ] that @xmath159 ( @xmath160 ) , which gives a multiplicative factor of about 12 ( 3 ) . for a given reconstructed @xmath1 ,
the most likely true @xmath1 is about @xmath161 ( @xmath162 ) smaller . in comparison , for c / a ( @xmath0 ) , with @xmath163 ( @xmath164 ) ( fig .
[ fig : ptshift ] ) and @xmath165 ( similar for @xmath0 ) there is a partial compensation between factors of @xmath166 ( @xmath167 ) and @xmath168 coming respectively from the shift and dispersion , yielding an overall factor of about @xmath169 ( @xmath170 ) , while the most likely true @xmath1 is about @xmath171 ( @xmath172 ) smaller than the reconstructed @xmath1 . at the lhc ( @xmath74 ) ,
( [ eq : hard - spectrum ] ) is a less accurate approximation . nevertheless , for @xmath173
, it is not too unreasonable to take @xmath174 and examine the consequences . for anti-@xmath0 ( respectively c / a(filt ) )
, we have @xmath175 ( also for c / a(filt ) ) and @xmath176 ( @xmath177 ) , giving a multiplicative factor of @xmath178 ( @xmath179 ) , i.e.  far smaller corrections than at rhic . for a given reconstructed @xmath1 ,
the most likely true @xmath1 is about @xmath180 ( @xmath181 ) smaller , rather similar to the values we found at rhic ( though smaller in relative terms , since the @xmath1 s are higher ) , with the increase in @xmath69 being compensated by the decrease in @xmath182 .    in the lhc case , it is also worth commenting on the results for the @xmath0 algorithm , since this is what was used in ref .
@xcite : we have @xmath183 and @xmath176 , giving a multiplicative factor of @xmath184 , which is consistent with the near perfect agreement that was seen there between the @xmath5 and subtracted @xmath185 spectra .
that agreement does not however imply perfect reconstruction , since the most likely @xmath186 is about @xmath162 lower than @xmath187 .
though the above numbers give an idea of the relative difficulties of using different algorithms at rhic and lhc , experimentally what matters most will be the systematic errors on the correction factors ( for example due to poorly understood non - gaussian tails of the @xmath26 distribution ) .
note also that a compensation between shift and dispersion factors , as happens for example with the c / a algorithm , is unlikely to reduce the overall systematic errors .
while the goal of this paper is not to discuss the issue of `` fake - jets '' in detail , it is a question that has been the subject of substantial debate recently ( see for example @xcite ) . here , therefore , we wish to devote a few words to it and discuss how it relates to our background - subtraction results so far .    in a picture in which the soft background and the hard jets are independent of each other
, one way of thinking about a fake jet is that it is a reconstructed jet ( with significant @xmath1 ) that is due not to the presence of an actual hard jet , but rather due to an upwards fluctuation of the soft background .
the difficulty with this definition is that there is no uniquely - defined separation between `` hard '' jets and soft background .
this can be illustrated with the example of how hydjet simulates rhic collisions : one event typically consists of a soft hydro background supplemented with @xmath188 @xmath5 collisions , each simulated with a minimum @xmath1 cut of @xmath189 on the @xmath190 scattering .
to some approximation , the properties of the full heavy - ion events remain relatively unchanged if one modifies the number of @xmath5 collisions and corresponding @xmath1 cut and also retunes the soft background . the fact that this changes the number of hard jets provides one illustration of the issue
that the soft / hard separation is ill - defined . additionally , while there are @xmath188 semi - hard @xmath5 collisions ( @xmath191 mostly central semi - hard jets event separately to obtain a long list of @xmath5 jets from all the separate hard events . ] ) in an event , there is only space within ( say ) the acceptance of rhic for @xmath192 jets .
thus there is essentially no region in an event which does not have a semi - hard jet . from this point of view , every reconstructed jet corresponds to a hard @xmath5 jet and there are no fake jets at all .      for inclusive analyses , such as a measurement of the inclusive jet spectrum ,
this last point is particularly relevant , because every jet in the event contributes to the measurement .
then , the issue of fakes can be viewed as one of unfolding . in that respect
it becomes instructive , for a given bin of the reconstructed heavy - ion @xmath1 , to ask what the corresponding matched @xmath5 jet transverse momenta were .
specifically , we define a quantity @xmath193 , the distribution of the @xmath5 `` origin '' , @xmath194 , of a heavy - ion jet with subtracted transverse momentum @xmath187 .
if the origin @xmath193 is dominated by a region of @xmath194 of the same order as @xmath187 , then that tells us that the jets being reconstructed are truly hard . if , on the other hand , it is dominated by @xmath194 near zero , then that is a sign that apparently hard heavy - ion jets are mostly due to upwards fluctuations of the background superimposed on low-@xmath1 @xmath5 jets , making the unfolding more delicate .    in figure
[ fig : origin ] we show the origins of heavy - ion jets as determined in our hydjet simulations .
jets can match a single heavy - ion jet , i.e.  have at least half their @xmath1 contained in the heavy - ion jet . in evaluating @xmath193
we take only the highest-@xmath1 matched jet .
if there is no matched jet ( this occurs only rarely ) then we fill the bin at @xmath195 .
note also that since we are not explicitly embedding hard jets , all @xmath5 jets in the events have undergone hydjet s quenching . ]
the upper row provides the origin plots for anti-@xmath0 jets at rhic .
each plot corresponds to one bin of @xmath187 , and shows @xmath193 as a function of @xmath194 . at moderate @xmath187 , the @xmath196 bin ,
the origin is dominated by low @xmath186 .
this is perhaps not surprising , given the result in section  [ eq : shift - dispersion - on - spectrum ] that the @xmath186 origin is expected to be @xmath197 lower than @xmath187 for anti-@xmath0 jets  additionally , that result assumed an exponential spectrum for the inclusive jet distribution , whereas the distribution rises substantially faster towards low @xmath186 .
as @xmath187 increases one sees that the contribution of high @xmath186 jets increases , in a manner not too inconsistent with the expected @xmath197 shift , though the @xmath186 distribution remains rather broad and a peak persists at small @xmath186 .
these plots suggest that an inclusive jet distribution measurement with the anti-@xmath0 algorithm at rhic is not completely trivial since , up to rather large @xmath187 , one is still sensitive to the jet distribution at small values of @xmath186 where the separation between `` hard '' jets and the soft medium is less clear .
nevertheless , two points should be kept in mind : firstly , the upper row of fig .
[ fig : origin ] shows that different @xmath187 have complementary sensitivities to different parts of the @xmath186 spectrum .
thus it should still be possible to `` unfold '' the @xmath187 distribution to obtain information about the @xmath186 , unfolding being in any case a standard part of the experimental correction procedure .
secondly star quotes  @xcite a @xmath149 smaller value for @xmath25 than the @xmath198 that we find in hydjet .
such a reduction can make it noticeably easier to perform the unfolding .
the impact of a reduction in @xmath25 is illustrated in the lower row of fig .
[ fig : origin ] , which shows the result for c / a(filt ) . here
, even the @xmath196 bin for @xmath187 shows a moderate-@xmath1 peak in the distribution of @xmath186 , and in the @xmath199 bin the low-@xmath1 `` fake '' peak has disappeared almost entirely .
furthermore , the @xmath186 peak is centred about @xmath162 lower than the centre of the @xmath187 bin , remarkably consistent with the calculations of section  [ eq : shift - dispersion - on - spectrum ] .
overall , therefore , unfolding with c / a(filt ) will be easier than with anti-@xmath0 .
corresponding plots for the lhc are shown in fig .
[ fig : origin_lhc ] .
while c / a(filt ) s lower dispersion still gives it an advantage over anti-@xmath0 , for @xmath200 , anti-@xmath0 does now reach a domain where the original @xmath5 jets are themselves always hard .
procedures to reject fake jets have been proposed , in  @xcite .
they are based on a cut on ( collinear unsafe ) jet shape properties and it is thus unclear how they will be affected by quenching and in particular whether the expected benefit of cutting the low-@xmath186 peak in figs .
[ fig : origin ] and [ fig : origin_lhc ] outweighs the disadvantage of potentially introducing extra sources of systematic uncertainty at moderate @xmath1 .
one final comment is that experimental unfolding should provide enough information to produce origin plots like those shown here .
as part of the broader discussion about fakes it would probably be instructive for such plots to be shown together with the inclusive - jet results .
an example of an exclusive analysis might be a dijet study , in which one selects the two hardest jets in the event , with transverse momenta @xmath201 and @xmath202 , and plots the distribution of @xmath203 . here
one can define `` fakes '' as corresponding to cases where one or other of the jets fails to match to one of the two hardest among all the jets from the individual @xmath5 events .
this definition is insensitive to the soft / hard boundary in a simulation such as hydjet , because it naturally picks out hard @xmath5 jets that are far above that boundary .
and , by concentrating on just two jets , it also evades the problem of high occupancy from the large multiplicity of semi - hard @xmath5 collisions .
this simplification of the definition of fakes is common to many exclusive analyses , because they tend to share the feature of identifying just one or two hard reference jets .
the specific case of the exclusive dijet analysis has the added advantage that it is amenable to a data - driven estimation of fakes .
one divides the events into two groups , those for which the two hardest jets are on the same side ( in azimuth ) of the event and those in which they are on opposite sides ( a related analysis was presented by star in ref .  @xcite ) . for events in which one of the two jets is `` fake ,
'' the two jets are just as likely to be on the same side as on the opposite side .
this is not the case for non - fake jets , given that the two hardest `` true '' jets nearly always come from the same @xmath5 event and so have to be on opposite sides .
, it is nearly always the case that the two hardest jets come from the same @xmath5 event . at the lhc , this happens above @xmath204 .
] thus by counting the number of same - side versus opposite - side dijets in a given @xmath205 bin , one immediately has an estimate of the fake rate .
events , if one has only limited rapidity acceptance then the same - side / opposite - side separation is not infrared safe , because of events in which only one hard jet is within the acceptance and the other `` jet '' is given by a soft gluon emission .
thus to examine the same - side / opposite - side separation in plain @xmath5 events with limited acceptance , one would need to impose a @xmath1 cut on the second jet , say @xmath206 . ]
this is illustrated in fig .
[ fig : dijets ] , which shows the distribution at rhic of the full , subtracted @xmath205 result , together with its separation into opposite - side and same - side components .
one sees that in the peak region , the opposite and same - side distributions are very similar , indicating a predominantly `` fake '' origin for at least one of two hardest jets ( they are not quite identical , because there is less phase - space on the same side for a second jet than there is on the away side ) . however above a certain full , reconstructed @xmath205 value , about @xmath207 for anti-@xmath0 and @xmath208 for c / a(filt )
the same - side distribution starts to fall far more rapidly than the opposite - side one , indicating that the measurement is now dominated by `` true '' pairs of jets .
the lhc results , fig .
[ fig : dijets_lhc ] , are qualitatively similar , with the same - side spectrum starting to fall off more steeply than the opposite - side one around @xmath209 for the anti-@xmath0 algorithm and @xmath210 for c / a(filt )
.    one can also examine origin plots for @xmath211 , in analogy with the monte carlo analysis of section  [ sec : fakes : incl ] .
for brevity , we refrain from showing them here , and restrict ourselves to the comment that in the region of @xmath211 where the result is dominated by opposite - side pairs , the origin plots are consistent with a purely hard origin for the dijets .
in this article we have presented the results of a systematic study of heavy - ion jet reconstruction with jet - area based background subtraction , building on the brief initial proposal of @xcite .
the questions we have examined include those of the choice of range for estimating the background , the choice of algorithm for the jet finding and the robustness of the reconstruction with respect to quenching effects and collision centrality .
we have found that there is little difference between various ranges , as long as they are chosen to be localised to the vicinity of the jet of interest and of sufficient size ( at least 4 units of area for jets with @xmath11 ) . in comparing different algorithms we examined the systematic offset and the dispersion in the reconstructed jet @xmath1 .
the offset can be brought close to zero by using the anti-@xmath0 algorithm , while the @xmath0 algorithm has the largest offset ; the cambridge / aachen ( c / a ) algorithm with filtering also gives a small offset , however this seems to have been due to a fortuitous cancellation between two only partially related effects .
the dispersion is comparable for anti-@xmath0 , c / a and @xmath0 , but significantly smaller for c / a(filt ) ( except at high transverse momenta for lhc ) , as a consequence of its smaller jet area . among the different algorithms
, anti-@xmath0 is the most robust with respect to quenching effects , and c / a(filt ) seems reasonably robust at rhic , though a little less so at the lhc .
the precise numerical results for offset and dispersion can depend a little on the details of the simulation and the analysis , however the general pattern remains .
overall our results indicate that the area - based subtraction method seems well suited for jet reconstruction in heavy - ion collisions .
two jet - algorithm choices were found to perform particularly well : anti-@xmath0 , which has small offsets but larger fluctuations , and c / a with filtering , for which the offsets may be harder to control , but for which the fluctuations are significantly reduced , with consequent advantages for the unfolding of experimentally measured jet spectra .
ultimately , we suspect that carrying out parallel analyses with these two choices may help maximise the reliability of jet results in hi collisions .
we wish to thank n.  armesto , h.  caines , d.  denterria , p.  jacobs and c.  salgado for numerous helpful comments on the manuscript .
we are grateful to them and to b.  cole , l.  cunqueiro , y - s .
lai , m.  ploskon , j.  puthschke , t.  renk , c.  roland , s.  salur , u.  wiedemann and many others , for useful and stimulating conversations on the topic of jets in heavy ion collisions over the past couple of years .
we also thank j.  andersen for helpful suggestions to obtain some of the computing resources used for this paper .
gs would like to acknowledge his affiliation to the brookhaven national laboratory while a large part of this work was performed and jr would like to acknowledge his affiliation to the lpthe in paris during its initial stages .
this work has been supported by the french agence nationale de la recherche , under grants anr-05-jcjc-0046 - 01 and anr-09-blan-0060 and the u.s .
department of energy under contract de - ac02 - 98ch10886 .
in section [ sec : subtraction ] , we gave an estimate of the minimum size of a range one should require for determining @xmath8 , given a requirement that fluctuations in the determination of @xmath8 should be moderate .
we give the details of the computation in this appendix .
we start from the fact that the error made on the estimation of the background density @xmath8 will translate into an increase of the dispersion @xmath25 : on one hand , the dominant contribution to @xmath25 comes from the intra - event fluctuations of the background _ i.e._is of order @xmath212 with @xmath213 the jet area ; on the other hand , the dispersion @xmath214 of the misestimation of @xmath8 leads to an additional dispersion on the reconstructed jet @xmath1 of @xmath215 .
adding these two sources of dispersion in quadrature and using the result from section  3.4 of @xcite , @xmath216 with @xmath217 the area of the range under consideration , we get @xmath218 if we ask _ e.g._that the contribution to the total dispersion coming from the misestimation of the background be no more than a fraction @xmath52 of the total @xmath25 , then we obtain the requirement @xmath219 for anti-@xmath0 jets of radius @xmath28 , with @xmath220 , this translates to @xmath221 where the numerical result has been given for @xmath55 . for @xmath11
, it becomes @xmath222 .
we can also cast this result in terms of the number of jets that must be present in @xmath37 . assuming that the jets used to estimate @xmath8 have a mean area of @xmath223 , or c / a algorithms for the background estimation  @xcite . ]
we find a minimal number of jets , @xmath224 where , as before , we have taken @xmath225 .
taking the numbers quoted above and @xmath65 , as used in the main body of the article , this gives @xmath226 .      from section  3.3 of @xcite
, we know that the presence of hard jets and initial - state radiation leads to a bias in the extraction of @xmath8 of @xmath227 where @xmath228 is a numerical constant and @xmath229 is the average number of `` hard '' jets ( those above the scale of the background fluctuations , including initial - state radiation ) ; @xmath229 is given by @xmath230 where @xmath62 is the number of `` born '' partons from the underlying @xmath231 scattering that enter the region @xmath232 , while @xmath233 is the first coefficient of the qcd @xmath234-function . in the context of @xcite ,
principally directed towards a study of the ue in @xmath5 collisions , @xmath69 was rather small , causing the second term of eq .
( [ eq : nhard ] ) , associated with initial - state radiation above a scale @xmath235 , to be comparable in size to the first term . in our case , @xmath69 is significantly larger and this reduces the impact of the second term sufficiently that we can ignore it .
we thus arrive at the result @xmath236    note that the presence of hard jets and initial - state radiation also affects the fluctuations in the misestimation of @xmath8 and this should in principle have been included in the estimates of appendix  [ app : minrange - fluct ] .
however , while the effect is not completely negligible , to within the accuracy that is relevant for us ( a few tens of percent in the estimation of a minimal @xmath217 ) it does not significantly alter the picture outlined there .
we have seen from fig .  [
fig : ptshift ] in section [ sec : choice_alg ] that the subtraction differs when we use the c / a algorithm with and without filtering .
since this difference is not due to back - reaction ( see fig .
[ fig : ptshift_br ] ) , it has to be due to the subtraction itself .
the difference comes from a bias introduced by the selection of the two hardest subjets during filtering .
the dominant contribution comes when only one subjet , that we shall assume harder than all the others , contains the hard radiation , all the other subjets being pure background . in that case , the selection of the hardest of these pure - background subjets as the second subjet to be kept tends to pick positive fluctuations of the background .
this in turn results in a positive offset compared to pure c / a clustering , as observed in section [ sec : choice_alg ] .    to compute the effect analytically ,
let us thus assume that we have one hard subjet and @xmath237 pure - background subjets of area @xmath238  @xcite .
after subtraction , the momentum of each of the pure - background subjet can be approximated as having a gaussian distribution of average zero and dispersion @xmath239 .
assuming that the `` hard '' subjet s transverse momentum remains larger than that of all the background jets , the 2 subjets that will be kept by the filter are the hard subjet ( subtracted ) and the hardest of all the subtracted background jets .
the momentum distribution of the latter is given by the maximum of the @xmath237 gaussian distributions .
gaussian - distributed random numbers , but with the result replaced by zero if all of them are negative . in the calculations here
we ignore this subtlety , since we will have @xmath240 and only @xmath241 of the time are three gaussian - distributed random numbers all negative . ]
we are only interested here in computing the average bias introduced by the filtering procedure , which is then given by @xmath242 for the typical case @xmath122 and @xmath243 , one finds @xmath244 if we insert in that expression the typical values for the fluctuations quoted in section [ sec : results ] and @xmath11 , we find average biases of 2  gev for rhic and 4.5  gev at the lhc , which are in good agreement with the differences observed between c / a with and without filtering in fig .
[ fig : ptshift ] .
we stated , in section  [ sec : dispersion ] , that the increase of the dispersion at high @xmath1 seen in fig .  [
fig : dispersion ] was mainly due to back - reaction .
this is made explicit in fig .
[ fig : br - v - full - dispersion ] , which decomposes the dispersion into its two components : that associated with the back - reaction , @xmath245 and that associated with background fluctuations and misestimation of @xmath8 ( defined as @xmath246^{\frac12}$ ] ) .
one sees that the background - fluctuation component is essentially independent of @xmath1 , while the back - reaction dispersion has a noticeable @xmath1 dependence .
this is the case because the back - reaction dispersion is dominated by rare events in which two similarly hard subjets are separated by a distance close to @xmath28 ( specifically by @xmath247 with @xmath248 ) .
in such a configuration , the background s contribution to the two subjets can affect whether they recombine and so lead to a large , @xmath249 , change to the jet s momentum . in the limit of a uniform background , @xmath250 , this can be shown to occur with a probability of order @xmath251 .
thus the contribution to the average shift @xmath24 is proportional to @xmath252 ( which in a full analysis is found to be enhanced by a logarithm @xcite for the @xmath0 and c / a algorithms ) , while the contribution to @xmath253 goes as @xmath254 , and so leads to a dispersion that should grow asymptotically as @xmath255 ( fig .
[ fig : br - v - full - dispersion ] is , however , probably not yet in the asymptotic regime ) .
it is worth keeping in mind that even though rare but large back - reaction dominates the overall dispersion , it will probably not be the main contributor in distorting the reconstructed jet spectrum .
such distortions come from upwards @xmath26 fluctuations , whereas large back - reaction tends to be dominated by downwards fluctuations .
the reason is simple : in order to have an upwards fluctuation from back - reaction , there must be extra @xmath1 near the jet in the original @xmath5 event .
this implies the presence of a harder underlying @xmath231 scattering than would be deduced from the jet @xmath1 , with a corresponding significant price to pay in terms of more suppressed matrix elements and pdfs .
figure  [ fig : br - v - full - dispersion ] also shows that the non - back - reaction component is nearly independent of @xmath1 .
this is expected since the anomalous dimension of the jet area is zero for anti-@xmath0 and small for c / a ( with or without filtering ) , and in any case leads to a weak scaling with @xmath1 , as @xmath121 .
furthermore , there is roughly a factor of @xmath125 between anti-@xmath0 and c / a(filt ) , as expected based on a proportionality of the dispersion to the square - root of the jet area .      one further source of @xmath26 fluctuations can come from imperfect estimation of the area of the jets .
we recall that throughout this article we have used soft ghosts , each with area of @xmath256 , in order to establish the jet area .
that implies a corresponding finite resolution on the jet area and related poor estimation of the exact edges of the jets , which can have an impact on the amount of background that one subtracts from each jet , and , consequently , on the final dispersion .
it is therefore interesting to see , figure  [ fig : ghost - area - impact ] , that the dispersion @xmath25 is reduced by about @xmath126 if one lowers the ghost area to @xmath257 .    while this does nt affect any of the conclusions of our paper ,
it does suggest that for a full experimental analysis there are benefits to be had from using a ghost area that is smaller than the default fastjet setting of @xmath256 .
s.  catani , y.  l. dokshitzer , m.  h. seymour , and b.  r. webber , _ nucl .
* b406 * ( 1993 ) 187 ; s.  catani , y.  l. dokshitzer , m.  olsson , g.  turnock , and b.  r. webber , _ phys .
lett . _ * b269 * ( 1991 ) 432 ; s.  d. ellis and d.  e. soper , _ phys . rev . _
* d48 * ( 1993 ) 31603166 , [ http://arxiv.org/abs/hep-ph/9305266[hep-ph/9305266 ] ] .
y.  l. dokshitzer , g.  d. leder , s.  moretti , and b.  r. webber , _ jhep _ * 08 * ( 1997 ) 001 , [ http://arxiv.org/abs/hep-ph/9707323[hep-ph/9707323 ] ] ; m.  wobisch and t.  wengler , http://arxiv.org/abs/hep-ph/9907280[hep-ph/9907280 ] .
m.  cacciari , g.  p.  salam and g.  soyez , jhep * 0804 * ( 2008 ) 063 [ http://arxiv.org/abs/0802.1189[arxiv:0802.1189 ] ] .
g.  p. salam and g.  soyez , _ jhep _ * 05 * ( 2007 ) 086 , [ http://arxiv.org/abs/0704.0292[arxiv:0704.0292 ] ] ; see also , http://projects.hepforge.org/siscone .
o.l . kodolova _
et al _ , `` study of @xmath258+jet channel in heavy ion collisions with cms , '' cms note-1998/063
. v.  gavrilov , a.  oulianov , o.  kodolova and i.  vardanian , `` jet reconstruction with pileup subtraction , '' cms rn 2003/004 .
n.  grau for the atlas collaboration , j.  phys .
g * 35 * ( 2008 ) 104040 [ http://arxiv.org/abs/0805.4656[arxiv:0805.4656 ] ] .
s.  salur [ star collaboration ] , eur .
j.   c * 61 * ( 2009 ) 761 [ http://arxiv.org/abs/0809.1609[arxiv:0809.1609 ] ] .
y.  -s .
lai , b.  a.  cole , .
n.  grau , b.  a.  cole , w.  g.  holzmann , m.  spousta and p.  steinberg [ atlas collaboration ] , eur .
j.   c * 62 * ( 2009 ) 191 [ http://arxiv.org/abs/0810.1219[arxiv:0810.1219 ] ] .
y.  s.  lai for the phenix collaboration , http://arxiv.org/abs/0911.3399[arxiv:0911.3399 ] .
m.  estienne [ alice collaboration ] , http://arxiv.org/abs/0910.2482[arxiv:0910.2482 ] .
j.  kapitan , http://arxiv.org/abs/0911.4754[arxiv:0911.4754 ] .
m.  cacciari and g.  p.  salam , phys .
b * 659 * ( 2008 ) 119 [ http://arxiv.org/abs/0707.1378[arxiv:0707.1378 ] ] .
m.  cacciari , g.  p.  salam and g.  soyez , jhep * 0804 * ( 2008 ) 005 [ http://arxiv.org/abs/0802.1188[arxiv:0802.1188 ] ] .
m.  cacciari , g.  p.  salam and s.  sapeta , jhep * 1004 * , 065 ( 2010 ) [ http://arxiv.org/abs/0912.4926[arxiv:0912.4926 ] ] .
m.  cacciari , j.  rojo , g.  p.  salam and g.  soyez , jhep * 0812 * ( 2008 ) 032 [ http://arxiv.org/abs/0810.1304[arxiv:0810.1304 ] ] ; see also http://quality.fastjet.fr .
s.  salur , nucl .
phys .   a * 830 * ( 2009 ) 139c [ http://arxiv.org/abs/0907.4536[arxiv:0907.4536 ] ] ; m.  ploskon [ star collaboration ] , nucl .
a * 830 * ( 2009 ) 255c [ http://arxiv.org/abs/0908.1799[arxiv:0908.1799 ] ] .
i.  p.  lokhtin and a.  m.  snigirev , eur .
j.   c * 45 * ( 2006 ) 211 [ http://arxiv.org/abs/hep-ph/0506189[arxiv:hep-ph/0506189 ] ] .
i.  p.  lokhtin , l.  v.  malinina , s.  v.  petrushanko , a.  m.  snigirev , i.  arsene and k.  tywoniuk , comput .
commun .   * 180 * ( 2009 ) 779 [ http://arxiv.org/abs/0809.2708[arxiv:0809.2708 ] ] .
r.  baier , y.  l.  dokshitzer , a.  h.  mueller , s.  peigne and d.  schiff , nucl .
b * 484 * ( 1997 ) 265 [ http://arxiv.org/abs/hep-ph/9608322[hep-ph/9608322 ] ] ; b.  g.  zakharov , jetp lett .
* 65 * ( 1997 ) 615 [ http://arxiv.org/abs/hep-ph/9704255[hep-ph/9704255 ] ] .
j.  casalderrey - solana and c.  a.  salgado , acta phys .
b * 38 * ( 2007 ) 3731 [ http://arxiv.org/abs/0712.3443[arxiv:0712.3443 ] ] ; d.  denterria , http://arxiv.org/abs/0902.2011[arxiv:0902.2011 ] ; u.  a.  wiedemann , http://arxiv.org/abs/0908.2306[arxiv:0908.2306 ] ; a.  majumder and m.  van leeuwen , http://arxiv.org/abs/1002.2206[arxiv:1002.2206 ] ; d.  denterria and b.  betz , lect .
notes phys .
* 785 * ( 2010 ) 285 .
t. sjostrand , s. mrenna and p. skands , jhep05 ( 2006 ) 026 [ http://arxiv.org/abs/hep-ph/0603175[hep-ph/0603175 ] ] .
n.  armesto , l.  cunqueiro and c.  a.  salgado , eur .
j.   c * 63 * ( 2009 ) 679 [ http://arxiv.org/abs/0907.1014[arxiv:0907.1014 ] ] ; n.  armesto , l.  cunqueiro , c.  a.  salgado and w.  c.  xiang , jhep * 0802 * ( 2008 ) 048 [ http://arxiv.org/abs/0710.3073[arxiv:0710.3073 ] ] .
m.  cacciari and g.  p.  salam , phys .
b * 641 * ( 2006 ) 57 [ http://arxiv.org/abs/hep-ph/0512210[hep-ph/0512210 ] ] .
j.  m.  butterworth , a.  r.  davison , m.  rubin and g.  p.  salam , phys .
* 100 * ( 2008 ) 242001 [ http://arxiv.org/abs/0802.2470[arxiv:0802.2470 ] ] .
d.  krohn , j.  thaler and l.  t.  wang , jhep * 1002 * ( 2010 ) 084 [ http://arxiv.org/abs/0912.1342[arxiv:0912.1342 ] ] . s.  d.  ellis , c.  k.  vermilion and j.  r.  walsh , phys .
d * 81 * ( 2010 ) 094023 .
m.  cacciari , g.  p.  salam and g.  soyez , http://www.fastjet.fr .
j.  adams _ et al .
_ [ star collaboration ] , phys .
c * 70 * ( 2004 ) 054907 .
s.  s.  adler _ et al .
_ [ phenix collaboration ] , phys .
c * 71 * ( 2005 ) 034908 [ erratum - ibid .
c * 71 * ( 2005 ) 049901 ] .
b.  b.  back _ et al .
_ , phys .
* 91 * ( 2003 ) 052303 [ http://arxiv.org/abs/nucl-ex/0210015[nucl-ex/0210015 ] ] .
n.  armesto , arxiv:0903.1330 [ hep - ph ] .
s.  s.  adler _ et al .
_ [ phenix collaboration ] , phys .
* 91 * ( 2003 ) 182301 [ http://arxiv.org/abs/nucl-ex/0305013[nucl-ex/0305013 ] ] ; j.  adams _ et al . _
[ star collaboration ] , phys .
* 92 * ( 2004 ) 052302 [ http://arxiv.org/abs/nucl-ex/0306007[nucl-ex/0306007 ] ] ; c.  alt _ et al . _
[ na49 collaboration ] , phys .
c * 68 * ( 2003 ) 034903 [ http://arxiv.org/abs/nucl-ex/0303001[nucl-ex/0303001 ] ] ; b.  b.  back _ et al .
_ [ phobos collaboration ] , phys .
c * 72 * ( 2005 ) 051901 [ http://arxiv.org/abs/nucl-ex/0407012[nucl-ex/0407012 ] ] .
j.  y.  ollitrault , phys .
d * 46 * ( 1992 ) 229 ; a.  m.  poskanzer and s.  a.  voloshin , phys .
c * 58 * ( 1998 ) 1671 [ http://arxiv.org/abs/nucl-ex/9805001[nucl-ex/9805001 ] ] .
i.  p.  lokhtin and a.  m.  snigirev , http://arxiv.org/abs/hep-ph/0312204[hep-ph/0312204 ] .
k.  zapp , g.  ingelman , j.  rathsman , j.  stachel and u.  a.  wiedemann , eur .
j.   c * 60 * ( 2009 ) 617 [ http://arxiv.org/abs/0804.3568[arxiv:0804.3568 ] ] .
b.  schenke , c.  gale and s.  jeon , phys .
c * 80 * ( 2009 ) 054913 .