the leading order bfkl kernel  @xcite , derived from feynman diagrams in momentum space , has been investigated in much detail . as a function of the transverse momenta , it is a meromorphic function . in the color singlet exchange channel
, it describes the pomeron contribution in pqcd , consisting of ladder diagrams with reggeized gluons .
this bfkl pomeron couples to the impact factors of colourless particles which , because of gauge invariance , vanish as one of the two reggeized gluons carries a zero transverse momentum .
this property allows us to modify the space of functions to which the pomeron wave function belongs . on the other hand , in the color octet exchange channel
the bfkl equation has the bootstrap solution , which represents a fundamental consistency property derived from the s - channel unitarity .
a manifestation of this bootstrap property also takes place for the color singlet state of two gluons and inside the coupling of three or more gluons to colorless particle impact factors .
in particular , it plays an important role in the odderon solution  @xcite which appears as a bound state of three reggeized gluons .
an important virtue of the color singlet kernel is its invariance under the mbius transformations in the space of transverse coordinates . exploiting the gauge invariance of the colorless impact factor to which the bfkl pomeron couples , the mbius symmetry allows to search solutions of the bfkl eigenvalues equation in the space of mbius functions which vanish as the coordinates of the two reggeized gluons coincide , and to define a spectral representation of the bfkl kernel in terms of conformal eigenfunctions . in this space of functions
, the bfkl kernel also enjoys the property of holomorphic separability .
translating back to the momentum representation , the kernel in the corresponding spectral form acquires @xmath4 function - like pieces which , because of the special properties of the impact factors , do not contribute to physical amplitudes .    as mentioned above , the bfkl kernel , when applied to the color octet wave function in momentum space , satisfies the bootstrap condition , which reflects the s - channel consistency of the bfkl calculation .
when transforming to transverse coordinates one finds that , in the space of mbius functions , this bootstrap condition can not be satisfied , i.e. the colour octet wave function lies outside this space of functions . in this paper
we will define a similarity transform , @xmath5 , which takes us from the mbius space of functions ( m space ) to another space of functions ( named ` analytic feynman diagram ' ( af ) space ) in which the bootstrap holds .
with the same transformation we can also define ` deformed ' mbius transformations and a ` deformed ' spectral representation of the bfkl kernel .
the paper is organized as follows . in the following section we review the properties of the bfkl kernel , both in momentum space and in the space of transverse coordinates . in section 3
we introduce the deformed representation , and we define the similarity transform which takes us from the mbius space ( @xmath6 space ) to the new deformed space ( @xmath7 space ) of functions . in section 4
we explicitly show that , in the @xmath7 representation , the bootstrap properties are fulfilled .
section 5 is devoted to the deformed representation of the conformal algebra which follows from the particular choice of the scalar product ; in particular we compute the transformed generators of the mbius group @xmath0 .
some details of our calculations are collected in two appendices .
in this section we give a brief review of the bfkl approach describing the dynamics of the reggeized gluons in lla of perturbative qcd .
let us start from the schrdinger - like bfkl equation @xcite describing the compound state of two reggeized gluons , @xmath8 here @xmath9 labels the colour representation of the two gluon state and in the singlet and octet channel one has respectively @xmath10 and @xmath11 .
the symbol @xmath12 denotes an integration in the transverse space in the case of the integral operator @xmath13 , while @xmath14 where @xmath15 is the @xmath16-channel angular momentum .
the eigenvalues of @xmath17 give the positions of singularities of the @xmath16 channel partial waves , related to the scattering amplitude by the mellin transformation in the variable @xmath18 , where @xmath19 is the squared total energy of colliding particles .
in lla , in the momentum representation , the gluon trajectory ( scaled by @xmath20 ) is given by the well known expression @xmath21 and the interaction term is defined by its action on the wave functions @xmath22
\phi(\mbox{\boldmath $ k$}_1^{\prime},\mbox{\boldmath $ q$}-\mbox{\boldmath $ k$}_1^{\prime } )   \nonumber\end{aligned}\ ] ] and @xmath23 .
let us note that this form is proper when the kernel acts on amputated ( without gluon propagators ) impact factors , otherwise the propagators should be removed from the kernel , as for example is the case for the expression in eq .
( [ bfkl_mom ] ) .
we stress that this form of the ll bfkl kernel is obtained directly from the feynman diagrams analysis in perturbative qcd and has a simple analytic behavior in the transverse momentum space .
in particular it does not contain @xmath4-functions @xmath24 and @xmath25 .    in the construction of the bfkl kernels
a very important property is the gluon reggeization which can be verified with the bootstrap relation for the amplitude with the octet quantum numbers @xcite .
the bootstrap relation is a consequence of the @xmath19-channel unitarity .
it claims , that the scattering amplitude with octet quantum numbers obtained as a solution of the corresponding bfkl equation should coincide with the born term multiplied by the regge factor @xmath26 .
because the bfkl equation was obtained by summing contributions from the multi - particle production described by the multi - regge amplitudes , the bootstrap relation valid in the leading and next - to - leading orders allows to connect the gluon trajectory and interaction terms . in lla
the bootstrap condition for the bfkl kernel can be written as @xmath27 where the constant @xmath28 is the wave function which can be conveniently obtained after rescaling any function depending only on @xmath29 , and the kernel @xmath30 acts on the amputated amplitudes .
let us note that it is crucial in the bootstrap relation in lla that the two gluons are located at the same point in the transverse coordinate plane .
it has played a crucial role in the discovery of the leading odderon solution in the lla of perturbative qcd  @xcite and also some general relation  @xcite between bound states of @xmath31 and @xmath32 reggeized gluons .    in the following
we shall use its equivalent form  @xcite , which has the virtue of being infrared finite and which uses the standard bfkl singlet kernel : @xmath33    on using an infrared mass regularization ( @xmath34 ) one may write @xmath35 where @xmath36 and @xmath37 , so that a relation equivalent to the bootstrap , but explicitely infraredly finite , is given by @xmath38    the 2-gluon kernel ( [ bfkleq ] ) in the singlet channel has been investigated in great details in the coordinate representation  @xcite . the amplitude for the scattering of colorless objects is factorized , in the high energy limit , in the product of the green s function ( which exponentiates the bfkl kernel ) and two impact factors .
the impact factors vanish as the momentum of one of the attached reggeized gluons goes to zero .
this permits to choose a special representation for the two - gluon propagator and for the full bfkl kernel acting in the space of the so called mbius functions ( for two gluon states these are functions @xmath39 such that @xmath40 )  @xcite .
this is a special choice among infinite others compatible with the gauge freedom of integrating a colorless impact factor with any function belonging to the set defined by the equivalence relation @xmath41 since the corresponding shift in the momentum representation is proportional to @xmath42 .
in such a representation the operators are mbius ( conformal ) invariant   @xcite . using complex coordinates , the bfkl hamiltonian @xmath43 acting on functions with propagators
included can be written in the operator form @xcite @xmath44 where @xmath45 , and we introduced the gluon holomorphic momenta @xmath46 . on the mbius space of functions
the holomorphic separability applies : @xmath47 and it is possible to find more easily the solutions of the homogeneous bfkl pomeron equation which belong to irreducible unitary representations of the mbius group .
in particular the symmetry generators in this representation are @xmath48 for two reggeized gluons one has @xmath49 and the casimir operator is defined as follows @xmath50 where @xmath51 and @xmath52 .
the eigenfunctions of the bfkl kernel are also eigenstates of two casimir operators and given by @xcite @xmath53 where @xmath54 , @xmath55 are conformal weights for the principal series of the unitary representations of the mbius group , @xmath31 is the conformal spin and @xmath56 is the anomalous dimension of the operator @xmath57 describing the compound state of two reggeized gluons ( note , that here and below we use other notations for conformal weights @xmath58 in refs .
the corresponding eigenvalues of the bfkl kernel are given by @xmath59    the action of the bfkl kernel can also be written , again on the space of mbius functions , after a duality transformation  @xcite , in terms of integral operator appearing in the dipole picture @xcite @xmath60    we now consider the pomeron eigenstates in the momentum representation .
different forms of the fourier transform of the function in eq .
( [ pomstatescoord ] ) have been given in the most general case .
one finds a sum of an analytic term @xmath61 , which has been given in an explicit way or in a integral form , and a @xmath4-like one @xmath62 .
more precisely we obtain @xmath63 where a bra - ket compact notation is also introduced . by explicit
computation one finds @xmath64 \frac{i^n}{2\pi } 2^{1-h-\bar{h } } \frac{\gamma(1-\bar{h})}{\gamma(h ) } q^{\bar{h}-1 } q^{*\ , h-1 } \ , , \label{pom_delta}\ ] ] where the complex form @xmath65 of any bidimensional transverse vector @xmath66 ( coordinate or momentum ) is used . in the calculation of a physical cross section one performs integrations with two - gluon colourless impact factors and this @xmath4-like term gives zero contribution , so that all the physics is related to the analytic term .
in fact the @xmath4-terms have been introduced when one has chosen to move to the mbius representation which is characterized by the simple and beautiful form in the coordinate representation .
the analytic part can be written , for example , as given in  @xcite wherein one finds @xmath67 with the coefficient @xmath68 given by @xmath69 the expression for @xmath70 in complex notation was written in terms of the hypergeometric functions @xmath71 another expression for @xmath72 which will be very useful for our purposes has been given in  @xcite in an integral form : @xmath73^{\bar{h}-1 } \left [ \frac{p^ * ( q^*-p^*)}{k_1^*-p^*}\right]^{h-1 } \ , ,   \label{ana_lip}\ ] ] where @xmath74    below we shall consider the completeness relation and the spectral representation for the bfkl kernel ( cf .
@xcite ) .
let us before discuss the scalar product in the space of functions where the kernel is acting on .
for the kernel in the momentum space associated to the feynman diagram derivation there is only one scalar product available .
if we consider non amputated functions of two momenta @xmath75 , i.e. with 2-dimensional propagators @xmath76 included , the scalar product is defined as @xmath77 with the integration measure @xmath78 which kills the extra propagators .    in the mbius space of functions ( in coordinate space ) for the principal series of the mbius group representation ,
whenever functions with conformal weight @xmath79 or @xmath80 are not considered , two possible choices are available  @xcite .
since the mbius functions contain propagators , one possibility is to remove , as before , the propagators in one function , acting with the operator @xmath81 .
the other possibility , related to the form of the casimir operator of the mbius group , is to use instead the measure @xmath82 .
the latter is not equivalent to the former when we consider the conformal wheights @xmath79 or @xmath80 .
note , that , providing that we go from the regge kinematics to the deep - inelastic scattering in the bjorken regime , the additional series of the unitary representations of the mbius group should be also used  @xcite ( cf . ref .
@xcite ) .
we now consider the completeness relation in coordinate space  @xcite , in the space of mbius functions , where the eigenfunctions of the casimir operator of the mbius group , defined in ( [ pomstatescoord ] ) , constitute a suitable spectral basis for the operators acting on the space of the colourless impact factors , which , we stress , never `` feel '' the presence of the terms with a @xmath4-distribution behavior in momentum space , i.e. those of eq .
( [ pom_delta ] ) . for the operators amputated ( spatial propagators removed ) to the left ( with notation @xmath83 and @xmath84 )
one can write @xcite @xmath85 where we have defined the weights @xmath86 , as well as @xmath87 for later use , @xmath88[\nu^2+(n+1)^2/4]}\ ] ] and @xmath89 .
if we try to extrapolate these operators to a wider domain which contains functions which no longer vanish for zero momenta , contrary to the colorless impact factor case , one needs some care .
for example , for the function , which in momentum space depends only on the total momentum , just the case we meet in the bootstrap relation in lla , one obtains the the result @xmath90 in coordinate space ( for the amputated impact factor ) .
but such generalized function is orthogonal to the basis constituted by the functions in eq .
( [ pomstatescoord ] ) since @xmath91 and therefore the application of the operators in eq .
( [ coord_relations ] ) on this function will give zero .
in other words the bootstrap relation can not be fullfilled in such a framework with the spectral representations given in eq .
( [ coord_relations ] ) .    but this should not be surprising , since the choice done thanks to the equivalence relation in eq .
( [ uv_transf ] ) , is based on a restriction of the action of the kernel to the functions similar to colorless impact factors .
moving to the momentum representation , one may write the fourier transform of the relations in ( [ coord_relations ] ) and the same considerations done above may be repeated in this case , leading to the same conclusions . in particular
one can not write a bootstrap relation in a spectral basis built upon the states of eq .
( [ pom_mom ] ) . on the other hand one
notices that the bfkl kernel @xmath92 in momentum space , as constructed from the feynman diagram analysis , has an analytic behavior which must be reflected also in its spectral decomposition .
it can also be noted that the behavior in @xmath93 of the @xmath4-like term ( [ pom_delta ] ) is @xmath94 , which is singular in the forward limit where analytic behavior is expected for fixed @xmath95 .
therefore it is natural to suggest the modified relations for the eigenfunction completeness and for the spectral representation of the bfkl kernel : @xmath96 where only the analytic contribution from each state of the spectral basis is used and the measure for integration is defined in eq .
( [ scal_prod ] ) , according to our choice of considering functions with propagators included .
below we shall verify this anzatz .
firstly we consider a space of functions ( with removed propagators ) which includes functions corresponding to the colorless impact factors and at least one function depending only on the total momentum , which is associated to a particular colored impact factor .
secondly we modify the spectral representation of the ll bfkl kernel , wherein the eigenvalues are assumed to be the same expressions of the conformal weights as for the restriction on the mbius space of function , and the eigenfunctions are deformed in order to be analytic in the momentum space .
a similar idea was already considered in the literature in an attempt to couple the ll bfkl pomeron to a quark  @xcite .
we proceed now with the definition of a transformation between the mbius ( m ) basis and the analytic feynman ( af ) basis .
let us observe that one can go from one basis to the other by a simple transformation , which reads in coordinate representation as follows    @xmath97    these relations are clearly non local and have been obtained by noting that @xmath98 for the case of even conformal spins ( @xmath31 ) the second transformation can be written in a simple local form on observing that @xmath99 while for odd conformal spin such a simple transformation is not possible since @xmath100 .
it is related to the fact , that in the last case both functions @xmath101 and @xmath102 satisfy the colour transparancy property @xmath103 , but only @xmath101 has the simple conformal properties .
it can be easily checked that @xmath104 so that @xmath5 is a 1 - 1 mapping with inverse really given by @xmath105 .
the green s function for the evolution in the rapidity @xmath106 based on the deformed analytic spectral basis was also written above ( see ( [ mom_relations ] ) ) . in the next section
we shall show that with this prescription one is recovering the correct relation compatible with the bootstrap requirement for the gluon reggeization in lla .    since the spectrum in both representations is the same it is natural to expect that the bfkl kernel satisfying the bootstrap relation is conformal invariant , in the new deformed analytic feynman basis obtained by a similarity transformation given by the operator @xmath5 .
namely , the well known generators of the mbius group @xmath107 could be extended to this basis as follows @xmath108 let us summarize the result . if one considers the impact factors derived from feynman diagrams calculations , they are the functions belonging to the af - space and not to the m - space . we remind that we have two different completeness relations based on the two spectral basis @xmath109 and @xmath110 and in these two spaces there are different definitions for scalar products and normalizations of wave functions .
nevertheless the cardinality of the two basis is the same .
the two spaces are related by the 1 - 1 mapping @xmath111 af - space @xmath112 m - space as previously discussed ( see also fig . [ fig_mapping ] ) and
one space is _ shifted _ with respect to the other .
moreover any operator @xmath113 in the af - space which can be decomposed on the spectral basis can be defined on the m - space as @xmath114 and viceversa @xmath115 .
this can be seen as an isometry between the mbius space and the analytical feynman space with the different scalar products in these spaces .
-@xmath28 mapping between the af - space and the m - space.,width=264 ]    the next interesting question is therefore if it is possible to construct explicitely the new representation defined in eq .
( [ derep ] ) .
we shall study this problem in the section @xmath116 , starting from analyzing the fourier transform of the mbius algebra in the m - space and later on moving to the af - space . in the next section
we will show that with the meromorphic prescription we indeed recover the bootstrap properties of the bfkl kernel .
for the sake of using a compact notation , let us introduce for some functions of transverse gluon momenta the corresponding wave functions for a color singlet state .
these are the _ power _
@xmath117 , the _ unity _ @xmath118 and the _ log _ @xmath119 states ( modulus the propagators ) , defined by their representations in the momentum space : @xmath120 since , as previoulsy discussed , one has for the mbius - invariant wave function @xmath121 the equality @xmath122 we obtain @xmath123 let us recall the completeness relation and the bfkl kernel written as @xmath124 on using the previous relations we write the bootstrap relation in the form of eq .
( [ bootstrapirfinite ] ) as @xmath125 ) was used . in order to show that our definitions ( [ gen_relations ] ) are correct
, we should prove that the coefficients of the conformal basis @xmath126 on both the l.h.s . and the r.h.s .
do coincide , i.e. @xmath127 it is therefore enough to calculate the integral given by @xmath128 .
verifying that this expression is finite it is sufficient also to calculate the first two terms of the taylor expansion in @xmath129 to verify the eq .
( [ coeff_bootstrap ] ) , thanks to eq .
( [ power_rel ] ) .
we also note that the left hand side of eq .
( [ coeff_bootstrap ] ) can be calculated easily in an independent way using the relations in eq .
( [ ana_delta ] ) and ( [ pom_delta ] ) , which give : @xmath130 where @xmath131/2 $ ] has support ( and value @xmath28 ) for any even integer number .
let us now find explicitely @xmath132 .
it is convenient to perform the calculations in the momentum space where we have defined the analytic element @xmath133 of the conformal basis .
taking for it the expression given in eq .
( [ ana_lip ] ) we write : @xmath134 using new integration variables induced by the change @xmath135 and @xmath136 with analogous relations for the complex conjugated ones , we obtains @xmath137 where @xmath138 it is convenient to define the function @xmath139 because one can find an explicit expression for it in terms of the generalized hypergeometric functions @xmath140 .
we illustrate some details of the calculation in appendix a.    it is sufficient to represent the result of the integration in power series of @xmath129 for the point @xmath141 : @xmath142    recalling the first two terms in the expansion ( see eq .
( [ d01 ] ) in appendix a ) one can see , on using the definitions of eq .
( [ power_rel ] ) and ( [ hapl ] ) , which imply @xmath143 and @xmath144 , that the relation in eq .
( [ coeff_bootstrap ] ) is verified and therefore also the bootstrap relation . as a check one may evaluate explicitely the quantity @xmath145 using the eq .
( [ cl ] ) .
it is easy to show that it coincides with the expression of eq .
( [ hau ] ) , computed in a completely different way.in order to give an alternative check of the representation proposed , we have also calculated explicitly the third term of the expansion of the bfkl green s function @xmath146 ( the first term is trivial and the second corresponds to the bootstrap relation ) , which contains the double iteration of the kernel acting on the unity ( momentum space ) .
the details are in the appendix b.    using the spectral representation one can compute an arbitrary number of iterations of the bfkl kernel or the action of the bfkl green s function on the impact factor of a quark or gluon line .
in this section we will examine the representation of the @xmath0 algebra in the @xmath7 space of funtions . to this end
we start from the familiar mbius algebra in the @xmath6 space and transform it into momentum space .
we then switch to the @xmath7 space : in order to retain the correct structure of the algebra we have to construct a new representation of the generator @xmath147 .
let us reconsider the generators of the conformal mbius group in @xmath148 dimensions in coordinate space given in eq .
( [ moebius_gen ] ) and the corresponding casimir operator ( [ moebius_cas ] )
. we will be interested in writing similar relations in momentum space .
we start from the conformal eigenfunctions ( [ pomstatescoord ] ) which from now on will carry the superscript @xmath6 , and we define their fourier transforms :    @xmath149(\mbox{\boldmath $ q$};\mbox{\boldmath $ k$}_{1},\mbox{\boldmath $ k$}_{2 } ) & = & \int d\mu ~e^{-i\mbox{\boldmath $ q$}\cdot \mbox{\boldmath $ \rho$}_{0}+i\mbox{\boldmath $ k$}_{1}\cdot \mbox{\boldmath $ \rho$}_{1}+i\mbox{\boldmath $ k$}_{2}\cdot \mbox{\boldmath $ \rho$}_{2}}e_{h{\bar{h}}}^{m}(\mbox{\boldmath $ \rho$}_{0};\mbox{\boldmath $ \rho$}_{1},\mbox{\boldmath $ \rho$}_{2})=   \nonumber \\ & = & ( 2\pi ) ^{2}\delta ^{(2)}(\mbox{\boldmath $ q$}-\mbox{\boldmath $ k$}_{1}-\mbox{\boldmath $ k$}_{2}){\tilde{e}}_{h{\bar{h}}}^{m}(\mbox{\boldmath $ k$}_{1},\mbox{\boldmath $ k$}_{2})\ , , \\
{ \tilde{e}}_{h{\bar{h}}}^{m}(\mbox{\boldmath $ \rho$}_{0};\mbox{\boldmath $ k$}_{1},\mbox{\boldmath $ k$}_{2})=\tilde{\mathcal{f}}[e_{h{\bar{h}}}^{m}](\mbox{\boldmath $ \rho$}_{0};\mbox{\boldmath $ k$}_{1},\mbox{\boldmath $ k$}_{2 } ) & = & \int d\tilde{\mu}~e^{i\mbox{\boldmath $ k$}_{1}\cdot \mbox{\boldmath $ \rho$}_{1}+i\mbox{\boldmath $ k$}_{2}\cdot \mbox{\boldmath $ \rho$}_{2}}e_{h{\bar{h}}}^{m}(\mbox{\boldmath $ \rho$}_{0};\mbox{\boldmath $ \rho$}_{1},\mbox{\boldmath $ \rho$}_{2})= \nonumber \\ & = & e^{i(\mbox{\boldmath $ k$}_{1}+\mbox{\boldmath $ k$}_{2})\cdot \mbox{\boldmath $ \rho$}_{0}}{\tilde{e}}_{h{\bar{h}}}^{m}(\mbox{\boldmath $ k$}_{1},\mbox{\boldmath $ k$}_{2})\ , , \\ { \tilde{e}}_{h{\bar{h}}}^{m}(\mbox{\boldmath $ k$}_{1},\mbox{\boldmath $ k$}_{2})=\tilde{\mathcal{f}}[e_{h{\bar{h}}}^{m}](\mbox{\boldmath $ \rho$}_{0}=0;\mbox{\boldmath $ k$}_{1},\mbox{\boldmath $ k$}_{2 } ) & = & \int d\tilde{\mu}~e^{i\mbox{\boldmath $ k$}_{1}\cdot \mbox{\boldmath $ \rho$}_{1}+i\mbox{\boldmath $ k$}_{2}\cdot \mbox{\boldmath $ \rho$}_{2}}e_{h{\bar{h}}}^{m}(0;\mbox{\boldmath $ \rho$}_{1},\mbox{\boldmath $ \rho$}_{2})\,,\end{aligned}\ ] ]    where @xmath150 , @xmath151 and @xmath152 , @xmath153 are the fourier transform operators in three and two variables respectively . in coordinate space the action of the generators on the eigenfunctions @xmath154 is : @xmath155 note , that the functions @xmath156 can be considered as the clebsch - gordon coefficients in the expansion of the product of the reggeon wave functions @xmath157 in the sum of the irreducible representations    @xmath158    where the conformal weights for the principal series of unitary representations are    @xmath159    and the coefficients @xmath160 are fixed with the use of the normalization conditions for the wave functions @xmath157 and @xmath161 . in the irreducible representations @xmath161
the mbius group generators are given below    @xmath162    the generators of the mbius group can be translated into the momentum space , using the relation @xmath163 : @xmath164 where @xmath165    one can calculate also the action of some bilinear combinations of the generators in the momentum representation    @xmath166 @xmath167 where the casimir operator in the  short  momentum representation is @xmath168 note that in our notations @xmath169 is a generator while @xmath170 is the casimir operator .
the algebra of these generators in momentum space is the same as in coordinate space : @xmath171={\tilde{m}}_{r}^{+}\,,\quad \lbrack { \tilde{m}}_{r}^{+},{\tilde{m}}_{r}^{-}]=-2{\tilde{m}}_{r}^{3}\,,\quad \lbrack { \tilde{m}}_{r}^{-},{\tilde{m}}_{r}^{3}]=-{\tilde{m}}_{r}^{-}\,.\ ] ]    while the action of the generator @xmath172 does not need special comments , it is interesting to look at others .
one finds @xmath173 the last relation follows from the action of the casimir operator @xmath174 and from relation ( [ momm3 m ] ) .
in particular , given that @xmath175 + \left ( { \tilde{m}}^{3}\right ) ^{2}={\tilde{m}}^{+}{\tilde{m}}^{-}+{\tilde{m}}^{3}\left ( { \tilde{m}}^{3}+1\right ) \
, , \label{casimir}\ ] ] where the commutation relation @xmath176=-2{\tilde{m}}^{3}$ ] has been used , one can see that indeed the relation ( [ mommminusm ] ) is satisfied .
these relations mean that in the m - space in momentum representation the mbius algebra is projected to @xmath177 { \tilde e}^m = { \tilde m}^+ { \tilde e}^m\ , , \quad { \tilde m}^-{\tilde m}^+ { \tilde e}^m = 2 { \tilde m}^3 { \tilde e}^m \ , , \quad 0\ , { \tilde e}^m=0\ , { \tilde e}^m\,.\end{aligned}\ ] ]    let us now move to the @xmath7 space .
we start , in coordinate space , from the definition @xmath178 where , according to eq .
( [ mapping_on_e ] ) , @xmath179 and @xmath180 is a projector ( @xmath181 ) defined by its action @xmath182 in order to compute the action of the mbius generators on the analytic part , @xmath183 , of the conformal eigenfunction we first examine the singular piece , @xmath184 .
first of all let us observe that the casimir operator acting on @xmath185 gives zero , as it is trivially checked in the coordinate representation .
this means that @xmath121 in such a case can be obtained with the action of a function of the casimir ( solving the eq .
( [ casimiraction ] ) for @xmath121 ) only before applying the projector @xmath180 . after
the application of @xmath180 one needs to define a new operator to extract the scaling parameter @xmath121 .
it is related to the fact , that @xmath186 is not an eigenfunction of the usual casimir operator .
we should construct a new casimir operator @xmath187 for this representation .    in coordinate representation
we can deduce the following relations for generators @xmath188 which coincides with their action for @xmath189 while the action of @xmath190 is different . according to this fact the action of mbius generators on @xmath191 can be written as @xmath192 and @xmath193 2^{-h-{\bar{h}}}\,.\end{aligned}\ ] ] therefore only the last generator @xmath190 ,
needs to be redefined in the @xmath7 space .
let us consider now again the momentum space for these functions .
first we analyze the action of the @xmath194 generator on @xmath195 , which we have written in terms of hypergeometric functions .
one can immediately see that @xmath196 for any @xmath197 .
this means again that we are not forced to use a non local operator @xmath198 ( solving an equation with the corresponding casimir @xmath199 ) to extract the dimension @xmath121 of an eigenstate , as we had to do in the coordinate representation .
the @xmath200 function is instead very peculiar since we know that the casimir gives zero .
this fact is also verified directly in momentum space , due to the deltas present in @xmath201 .
later we will show it explicitely .
instead we note that the action of the @xmath194 operator is good .
indeed , from the explicit form given in eq .
( [ pom_delta ] ) , we obtain @xmath202 this is compatible with the fact that the @xmath194 operator has the same action on the @xmath203 and @xmath195 functions .
the fact that we may choose to extract the dimension @xmath121 in momentum space using the @xmath194 operator may be useful to simplify the construction of the mappings @xmath5 and @xmath105 . on defining @xmath204 one may write the mapping introduced in eq . ( [ mapping_on_e ] ) as @xmath205 where in general @xmath206 and which imply @xmath207 and we have denoted @xmath208 , @xmath209 and similarly @xmath210 . clearly in the relation for
the mappings one can use different forms for the operator @xmath211 depending on which space it acts .    in momentum space
@xmath212 then one can write @xmath213=0 $ ] ( with some abuse of notation ) and things are further simplified .
proceeding along this line one can give a deformed representation @xmath214 which lives on the af - space . in order to do this ,
let us first study @xmath215 and @xmath216 c_h ( k_1+k_2)^{{\bar{h}}-1 } ( k_1^*+k_2^*)^{h-1 }   \nonumber \\ & = & 2i h ( h-1 ) \frac{1}{k_1^*+k_2^ * } { \tilde e}^\delta_h \,=\ , -h(h-1 ) \left ( { \tilde m}^+ \right)^{-1 } { \tilde e}^\delta_h   \nonumber \\ & = & -\left ( { \tilde m}^+ \right)^{-1 } h(h-1 ) \frac{f}{1-f } p { \tilde e}_h^a \ , .   \label{rel1}\end{aligned}\ ] ] note that @xmath217 has the same form on any of the spaces of functions we consider .
the resulting action of @xmath218 on @xmath219 has to be studied in order to see if it is possible to define a new representation for @xmath218 on af - space which satisfies the conformal algebra . but before let us note that from the previous relation we have @xmath220 so that we reobtain , recalling eq .
( [ casimir ] ) , @xmath221 in general we have defined the mbius generators on the af space ( spanned by the basis functions @xmath219 ) using the isometry with the m - space : @xmath222 let us look at each of the three generators . from the previous discussion , also in the coordinate representation , we observe that @xmath223 which means that in momentum space the actions of these two operators is the same in both m - space and af - space .
we therefore define for the af - space a new generator in momentum representation such that @xmath224 which puts into relation , as before , the action of the casimir operator such that @xmath225 together with the piece of the mbius algebra @xmath226=-2{\tilde m}_{af}^3 $ ] and the action of @xmath227 . according to the previous result ( see eqs .
( [ mminusa ] ) and ( [ rel1 ] ) ) , by construction the choice is @xmath228 where , in the last form , g is a non local operator such that @xmath229    let us finally verify that , with this new representation of @xmath214 , the algebra of the mbius generators is correct .
it is easy to see that @xmath230=-{\tilde m}_{af}^- \,,\ ] ] which is a relation that is projected to zero when applied to the function @xmath219 . in order to check
that the mbius algebra is closed we need to examine the relation @xmath231=+2{\tilde m}_{af}^3\ ] ] understood to act on @xmath219 .
one then has to verify that @xmath232 using eq .
( [ reldelta ] ) we derive also @xmath233 and making the difference between eqs .
( [ v1 ] ) and ( [ v2 ] ) we therefore obtain @xmath234 if @xmath217 commutes with @xmath235 this is equivalent to @xmath236 which coincides with the starting definition of the generator .
the commutation relation is fullfilled since the operator @xmath217 is just a multiplicative operator associated to the total momentum and therefore commutes with the action of the projector @xmath180 and moreover it commutes with the function of the casimir @xmath237 .
finally let us look at the relation between the two casimir operators in m - space and af - space : using @xmath238 + \left ( { \tilde{m}}_{af}^{3}\right ) ^{2}$ ] we find from the previous relations @xmath239 this last relation , applied to @xmath240 , gives @xmath241 which indeed is compatible with @xmath242 .
among the properties of the ll bfkl kernel the bootstrap relation connected to the gluon reggeization is a fundamental consistency condition .
this bootstrap property leads to an important feature of the bfkl kernel also in the color singlet state : this can be demonstrated most easily in momentum space where the bfkl kernel is meromorphic .
many previous investigations , exploiting the conformal symmetry of the kernel , have been carried out in the space of mbius functions ( @xmath6 ) . in this space of functions , however , the bootstrap relation for the bfkl kernel is not satisfied .    in this paper
we have defined a modified space of functions ( @xmath7 ) , in which the bootstrap property is valid . in particular , we have derived a spectral representation of the kernel which is also useful in evaluating the coupling of the bfkl green s function to colored impact factors , and we have verified explicitly that the bootstrap relation is fulfilled . we also have derived the corresponding represention of the mbius algebra ; in order to act in this modified space of functions , one of the mbius generators has to be deformed .
our discussion has been limited to the the case of two - gluon green s functions . since in both spaces , @xmath6 and @xmath7 , the eigenvalues of the casimir operator , @xmath170 and @xmath187 , resp .
, are the same , and since the hamiltonian for the pairwise interaction of two reggeized gluons can be expressed in terms of this casimir operator , one can expect that , for states consisting of more than two gluons all remarkable properties of the multi - colour bfkl dynamics derived in the mbius picture can be generalized to the af picture .
this includes , in particular , the holomorphic separability and integrability @xcite .
we are going to return to these interesting problems in future publications .
part of this work has been done while one of us ( l.n.l ) has been visiting the ii.institut f.theoretische physik , university hamburg , and the university of bologna , bologna .
two of us ( j.b . and
g.p.vacca ) gratefully acknowledge the hospitality of the nuclear physics institute in gatchina , st.petersburg .
the calculation of the integral in eq .
( [ integral ] ) is carried on noting that one can find a third order differential operator in @xmath243 and another in @xmath244 , related to the @xmath140 functions , which , when applied to the integrand , gives a total derivative in an integration variable so that , on applying the stokes theorem , one finds that @xmath245 satisfies the related differential equation in both the holomorphic and antiholomorphic sectors .
this allows to write it as a sum of products of the linearly independent solutions of the two differential equations and on imposing the single valueness and the correct normalization ( at some convenient point ) the full expression is found @xcite .
in particular this has the following , so called , conformal block structure @xmath246 where @xmath247 and @xmath248 are three independent solutions of the generalized hypergeometric differential equations in @xmath243 and @xmath249 .
the coefficients @xmath250 depend through @xmath251 functions on the conformal weights .
the general form of @xmath139 is pretty complicated but it simplifies considerably for @xmath141 using the relations  @xcite @xmath252 on using the above relation , after algebraic simplifications , we find @xmath253 with @xmath254 the expressions of @xmath255 and @xmath256 have a simple pole at @xmath257 and opposite residues for this pole , while @xmath258 is holomorphic at @xmath257 , therefore the full function in eq .
( [ resulti ] ) is not singular at @xmath257 . using the transformation properties of the gamma function ,
each of three terms can be reduced to a common multiplicative factor and expressions involving only trigonometric functions .
these terms are combined together to give the following compact form for @xmath259 @xmath260 now the first few derivatives with respect to @xmath129 can be easily calculated .
in particular we obtain , after some manipulations , @xmath261 \ , c_0 \ , .
\label{d01}\end{aligned}\ ] ]
we recompute here the double iteration of the bfkl kernel in momentum space acting on a constant function  @xcite and compare it with the expressions obtained with the use of the above spectral representation of the kernel ( and therefore of the green s function ) .
let us write before the result of the direct integration in momentum space : @xmath262 \ , .
\label{ksquare1analytic}\end{aligned}\ ] ] before giving some details of the computation let us consider what we obtain from the spectral approach . on using the spectral representation for the kernel , the relation ( [ ksquare1analytic ] )
can be written as @xmath263 a relation which we have checked numerically .
this check has been also performed in the case of the single iteration ( bootstrap relation ) .
the argument of the integration over @xmath264 oscillates with a very slow decay for fixed conformal spin , but summing over some tens of @xmath31 gives a good suppression of the tails and the integral can be computed .
finally we sketch the derivation of eq .
( [ ksquare1analytic ] ) .
there are three kinds of integral involved ; they can be calculated with the dimensional regularization and the feynman parameterization , in the case @xmath265.the needed integrals are : @xmath266 where @xmath267 , @xmath268 and @xmath269 is the gluon trajectory ( [ gluontrajectory ] ) rescaled by @xmath270.the first two integral can be easily calculated exactly , giving : @xmath271 \ , ,   \nonumber \\ j_1 ( \mbox{\boldmath $ k$}_1,\mbox{\boldmath $ k$}_2 ) & = & \pi^{1+\epsilon } \frac{\gamma(1-\epsilon)\gamma^2(\epsilon ) } { \gamma(2 \epsilon ) } \left [ \pi \cot \pi \epsilon + \log \frac{\mbox{\boldmath $ q$}^2}{\mbox{\boldmath $ k$}_1 ^ 2 } + \psi(2 \epsilon ) - \psi(1 ) \right ] ( \mbox{\boldmath $ k$}_1 ^ 2)^\epsilon =   \nonumber \\ & = & \pi^{1+\epsilon } \gamma(1-\epsilon ) ( \mbox{\boldmath $ q$}^2)^\epsilon \left [ \frac{1}{\epsilon^2 } + \frac{1}{\epsilon } \log \frac{\mbox{\boldmath $ q$}^2}{\mbox{\boldmath $ k$}_1 ^ 2 } - \frac{\pi^2}{6 } -\frac{3}{2}\log^2 \frac{\mbox{\boldmath $ q$}^2}{\mbox{\boldmath $ k$}_1 ^ 2 } + \mathcal{o}(\epsilon ) \right ] \ , .\end{aligned}\ ] ] the third integral is calculated to the order @xmath272 : @xmath273 \ , .\end{aligned}\ ] ] putting together all the pieces to construct the complete action of the kernel , all the singulatities cancel and the result , which is finite in the limit @xmath274 , correspond to the expression given in eq .
( [ ksquare1analytic ] ) .
j.  bartels , l.  n.  lipatov and g.  p.  vacca , phys .
b * 477 * ( 2000 ) 178 , [ arxiv : hep - ph/9912423 ] .
g.  p.  vacca , phys .
b * 489 * , 337 ( 2000 ) , [ arxiv : hep - ph/0007067 ] .
l.  n.  lipatov , sov .
jetp * 63 * ( 1986 ) 904 [ zh .
* 90 * ( 1986 ) 1536 ] .
j.  bartels , l.  n.  lipatov and g.  p.  vacca , nucl .
b * 706 * ( 2005 ) 391 , [ arxiv : hep - ph/0404110 ] .
l.  n.  lipatov , phys . lett .
* b309 * ( 1993 ) 394 ; _ high energy asymptotics of multi - colour qcd and exactly solvable lattice models _ , hep - th/9311037 , unpublished .
j.  bartels , m.  a.  braun , d.  colferai and g.  p.  vacca , eur.phys.j.c * 20 * ( 2001 ) 323 , [ arxiv : hep - ph/0102221 ] .
h.  j.  de vega and l.  n.  lipatov , phys .
d * 64 * , 114019 ( 2001 ) , [ arxiv : hep - ph/0107225 ] , phys .
d * 66 * ( 2002 ) 074013 , [ arxiv : hep - ph/0107225 ] .
s.  e.  derkachov , g.  p.  korchemsky , j.  kotanski and a.  n.  manashov , nucl .
b * 645 * ( 2002 ) 237 [ arxiv : hep - th/0204124 ] .
a.  h.  mueller and w.  k.  tang , phys .
b * 284 * ( 1992 ) 123 .
j.  bartels , j.  r.  forshaw , h.  lotter , l.  n.  lipatov , m.  g.  ryskin and m.  wusthoff , phys .
b * 348 * ( 1995 ) 589 [ arxiv : hep - ph/9501204 ] .
j.  s.  geronimo and h.  navelet , j.  math .
* 44 * ( 2003 ) 2293
. bateman project , a. erdelyi , higher trascendental functions , vol .
189 , mcgraw - hill ( 1953 ) .


the study of open heavy flavor production in relativistic @xmath18 collisions is sensitive to different kinds of strong interaction physics . because the leading order ( lo ) production mechanism is gluon fusion  @xcite , open heavy flavor production rates are directly related to modification of the gluon parton distribution function ( pdf ) , i.e.  shadowing or saturation  @xcite .
also , the initial and/or final state partons can scatter and lose energy in the cold nuclear medium  @xcite , thereby modifying and producing a nuclear modification of open heavy flavor production . recently
, the possibility of flow even in small collision systems such as @xmath18 has raised the question of modified charm momentum distributions @xcite .
modification to heavy quark production rates and kinematics in  collisions at the relativistic heavy ion collider ( rhic ) is well established .
electron production from open heavy flavor decay is enhanced  @xcite , while @xmath19 production is suppressed  @xcite at midrapidity . at positive rapidity , defined with the positive @xmath20 axis as the direction of the deuteron
, there is a suppression of heavy - flavor decay muons  @xcite and a larger suppression of @xmath19  @xcite .
while @xmath21@xmath13 correlations from open heavy flavor decays have not been published at rhic to date , correlations involving light flavor hadrons have shown modification in collisions at rhic .
a suppression has been observed of positive rapidity @xmath22 mesons associated with midrapidity trigger hadrons , especially in the back - to - back peak at @xmath23 , indicating @xmath24 scatterings  @xcite .
this suppression increases as @xmath25 , the fraction of the nucleon momentum carried by the gluon , decreases .
these results are in quantitative agreement with energy loss models  @xcite and saturation models  @xcite .
this paper presents measurements of azimuthal correlations of electron - muon pairs produced from heavy flavor decays , primarily @xmath17 , in and collisions using the phenix detector at rhic .
the heavy - flavor @xmath21@xmath13 correlations are free of backgrounds from other sources that contribute to other dilepton analyses ( @xmath26 or @xmath27 ) , such as resonance decay and drell - yan . while analysis of dilepton mass and @xmath28 provides a way to separate charm and bottom contributions , the azimuthal correlations have an important advantage for studying the charm production process .
the leading order production , @xmath29 and @xmath30 , will produce back - to - back open heavy flavor pairs that can semileptonically decay and produce azimuthally correlated @xmath21@xmath13 pairs .
next - to - leading order ( nlo ) processes like flavor excitation and gluon splitting produce much less correlated @xmath31 and thus much less correlated @xmath21@xmath13 pairs .
therefore , modification to different portions of the azimuthal correlations can be attributed to modifications of @xmath17 pairs from different production mechanisms . in energy loss models such as ref .
@xcite , a broadening of the back - to - back azimuthal correlation should accompany a suppression of the peak due to the multiple scattering that the incoming gluons and/or the outgoing @xmath17 undergo in the cold nuclear medium .
this paper is organized as follows .
the phenix detector is outlined in section  [ sec : phenix ] .
section  [ sec : method ] describes the details of the method used to measure the correlations , the background subtraction method , and the tests of the method .
section  [ sec : results ] presents the results in and compares them to monte carlo models .
the results are presented and compared to the results in section  [ sec : dauresults ] .
conclusions are given in section  [ sec : end ] .
the phenix detector at rhic is multi - purposed and optimized for precision measurements of electromagnetic probes for relativistic hadronic and heavy ion collisions .
a complete overview of the detector can be found in ref .
the data presented here are from 2006 and 2008 data taking at rhic .
figure  [ fig : phenix ] shows a schematic of the detector during those years .
this analysis uses the central spectrometer arms for electron detection and the forward rapidity muon spectrometer arms , labeled north and south in fig .
[ fig : phenix ] , for muon identification . for the 2008 collisions ,
the deuteron beam moves toward the north arm , which defines positive rapidity for both and .
the forward produced muons come from a high-@xmath25 parton in the deuteron interacting with a low-@xmath25 parton in the gold .
pythia@xcite indicates that the average @xmath25 of a parton producing a heavy flavor muon from 1@xmath32  gev/@xmath3 in the forward muon spectrometer is about @xmath33 .
this analysis focuses only on the muons measured in the north arm utilizing the deuteron beam as a probe of low-@xmath25 partons in the gold nucleus .
the central spectrometer comprises two arms subtending @xmath34 in azimuth and covering @xmath1 .
charged tracks are measured using a drift chamber ( dc ) and a set of multi - wire proportional chambers with pad readout ( pc1 and pc3 ) .
the dc measures the bend angle in the @xmath35 plane due to a central magnetic field directed along the beam axis .
pc1 is used to measure the longitudinal coordinate of the track .
these tracks are then projected into pc3 , where a hit is required to ensure high track quality .
the momentum resolution of the tracks in this data is @xmath36%@xmath371.16%@xmath7 , where @xmath7 is the total momentum measured in gev/@xmath3 .
electrons can be identified from associated hits in the ring imaging erenkov ( rich ) detector and the electromagnetic calorimeters ( emcal ) .
electrons above 17 mev/@xmath3 passing through the co@xmath38-filled rich will emit erenkov radiation .
the emcal comprises eight sectors , six of lead - scintillator and two of lead - glass , used to collect the energy from electron and photon showers .
the nominal energy resolution for the lead - scintillator and lead - glass is 8.1%@xmath39}}\oplus$]2.1% and 6.0%@xmath40}}\oplus$]0.9%@xcite , respectively .
the north muon spectrometer is located at @xmath41 and covers 2@xmath42 in azimuth .
the spectrometer measures tracks in the muon tracker ( mutr ) and the muon identifier ( muid ) .
prior to entering the muon arm , particles pass through approximately 20  cm of copper and 60  cm of iron
. particles that are not absorbed pass through the mutr , which comprise three stations of cathode strip chambers with multiple ionization regions and located inside a radial magnetic field . after the mutr , particles pass through the muid , which comprises five alternating steel absorbers and muid detector planes , called gaps , with iarocci tubes .
muid roads reconstructed from muid hits are projected back to mutr tracks and to the measured vertex to provide the complete information for a track through the spectrometer .
trigger and global event characterization in and are provided by the beam - beam counter ( bbc ) .
the bbc is a set of 64 hexagonal erenkov counters located from @xmath43 and covering full azimuth .
the vertex of the collision along the beam line ( @xmath44 ) is determined by the time difference between the bbcs on either side of the collision region .
the minimum bias ( mb ) trigger requires that there is at least one hit in each of the bbcs . from vernier scans and verified by monte carlo studies , the bbc mb trigger is sensitive to 55@xmath125% of the inelastic cross section and 88@xmath124% of the inelastic cross section@xcite .
the trigger used for this analysis is a combination of the bbc trigger and a deep muon trigger .
the deep muon trigger requires three or more muid gaps with a signal in both the @xmath25 and @xmath45 direction tubes and that the last pair of hits be in the last ( 5th gap ) or next to last gap ( 4th gap ) .
after quality cuts and requiring a vertex within 25  cm of the @xmath20=0 vertex , an integrated luminosity of 2.1 pb@xmath46 in and a -equivalent of 7.7 pb@xmath46 in was sampled .
the primary goal of this analysis is to identify @xmath47 where the opposite sign electron - muon pair is from the @xmath17 pair decay .        only muon candidates with @xmath48  gev/@xmath3 are used in the analysis because real muons with total momentum less than about 2.7  gev/@xmath3 are stopped in the muon arm before reaching the 5th ( and last ) gap .
single muon candidates are constructed from muid roads projected and matched to mutr tracks .
cuts on muid roads and mutr tracks are designed to reject hadrons that mimic a muon signal and to reject tracks that did not originate from the collision vertex . for the muid roads ,
at least three of five gaps with @xmath49 hit information are required , including a pair of hits in the 5th gap .
these muid roads must project back near the nominal vertex position .
those muons that do not typically come from beam - related backgrounds . for the mutr tracks , cuts that reject hadrons
are detailed in ref .
the muid roads are then projected and matched to mutr tracks at the 1st muid gap .
an identified muon candidate is the closest mutr track that matches a muid road within at least 10@xmath50 in slope and 10  cm in distance .
muon candidates are further restricted to @xmath4 . during both the and data taking periods
, there were backgrounds primarily from beam - related particles interacting with material in the accelerator upstream of phenix .
collimators were used in the accelerator to reduce this background but it was not totally eliminated . restricting the @xmath51 range of the muon candidates helped minimize this background .
electrons with @xmath2  gev/@xmath3 are identified by matching a track in dc , pc1 , and pc3 to a signal in the rich and a cluster in the emcal .
the relevant details on measuring electrons in phenix are given in ref .
@xcite . for this analysis ,
the projected track must match within 3@xmath52 in position to a cluster in the emcal .
clusters are also required to have a matching profile , when compared to an electromagnetic shower shape profile at the measured energy .
once a track matches both the rich and the emcal , an @xmath53 cut is applied , where it is required that the energy measured in the the emcal @xmath54 be approximately equal to the reconstructed track momentum @xmath7 .
this is sufficient to remove most combinatorial matches and background from real electrons resulting from long - lived particle decays occurring near the dc , which have mismeasured momentum .
a cut of -2@xmath52 to + 3@xmath52 from the mean @xmath53 in the data and -1.5@xmath52 to + 3@xmath52 from the mean in the data is applied .
the asymmetry of the cuts is due to the dominance of backgrounds below 2 or 1.5@xmath52 of the mean .
the tighter cut in the data was necessary because of the increased background from the hadron blind detector ( hbd ) support material not present during 2006 data taking .
after particle identification cuts have been applied to an event , all pairs of identified electrons and muons are formed in each of the four charge - sign combinations .
the fully corrected invariant - pair yield , calculated for each sign combination , is @xcite @xmath55    where @xmath56 is the number of sampled bbc triggered events , @xmath3 is the mb trigger bias accounting for events missed by the bbc trigger@xcite , @xmath57 and @xmath58 are the rapidity ranges of the electrons and muons , respectively , @xmath59 is the inclusive electron - muon pair yield , and @xmath60 is the mixed - event electron - muon pair distribution .
the two - particle acceptance and efficiency is corrected by the mixed - event technique , where electrons from one event are paired with muons from a different event .
pools of inclusive electrons and muons are kept in bins 2.5  cm - wide @xmath20 vertex bins , and , in the case of , 10%-wide centrality bins . when mixing events , the pair distribution is weighted by the @xmath45- and @xmath61-averaged efficiency of each particle , @xmath62 and @xmath63 .
both @xmath62 and @xmath63 were determined by generating single electrons and single muons with a flat distribution in @xmath28 , @xmath61 , @xmath64 or @xmath65 and collisions @xmath20-vertex location and running them through a geant-3 simulation of the phenix detector .
the output was subjected to the same analysis cuts applied to the data .
the efficiency is defined as the ratio of particles reconstructed through the analysis to the number simulated .
these simulations demonstrated that @xmath62 and @xmath63 are independent of the @xmath20-position of the event vertex , @xmath66 is independent of @xmath51 and @xmath67 has a slight @xmath51-dependence .
pair yields are reported with the average pseudorapidity @xmath68 , which include the @xmath51-dependence of both single inclusive muons and the single particle efficiency .
inclusive muon and electron candidates come from both heavy- and light - flavor decays and from misidentified hadrons .
the fully - corrected inclusive electron - muon pair yield for each sign combinations can be written as    @xmath69    here @xmath70 indicates the fully - corrected inclusive pair yield defined in eq .
[ eq : pairyield ] ; @xmath71 is the fully - corrected pair yield produced from a heavy flavor pair decay ; @xmath72 is the fully - corrected pair yield from correlating a heavy flavor decay product with a light flavor decay product ; and @xmath73 is the fully - corrected pair yield from correlating pairs of light - flavor decay products or misidentified hadrons .
pairs from the semileptonic decay of a @xmath17 pair have opposite signs .
[ eq : basecorr ] can be decomposed into its like- and unlike - sign pieces :    @xmath74    while semileptonic decays of @xmath75 can also produce both like- and unlike - sign @xmath21@xmath13 signals , in this analysis , pythia indicates that only about 1% of the final heavy - flavor @xmath21@xmath13 pair yield is from @xmath75 and is neglected .
if we assume muon ( electron ) candidates from light flavors are not charge - correlated with electron ( muon ) candidates from light flavors , then @xmath76 if only one of the pair is from heavy flavor , then again , we assume they are not charge correlated , and @xmath77 therefore , the heavy flavor @xmath21@xmath13 signal distributions is the difference between the unlike - sign and the like - sign inclusive correlations : @xmath78    figure  [ fig : sb ] shows the fully - corrected inclusive like - sign ( @xmath79 ) and unlike - sign ( @xmath80 ) @xmath21@xmath13 pair distributions in and .
the inset figures show the signal - to - background distributions given the assumptions above .    ) and unlike - sign ( @xmath81 ) distributions for ( a ) and ( b ) , as a function of @xmath6 .
the inset shows the unlike - like difference divided by the like - sign distribution , which is the heavy flavor signal - to - background in the inclusive unlike - sign distribution .
, title="fig : " ] ) and unlike - sign ( @xmath81 ) distributions for ( a ) and ( b ) , as a function of @xmath6 .
the inset shows the unlike - like difference divided by the like - sign distribution , which is the heavy flavor signal - to - background in the inclusive unlike - sign distribution .
, title="fig : " ]    we have checked the like - sign subtraction method using pythia leading order quantum chromodynamics ( qcd ) events . with all events containing a heavy quark in the final state removed ,
the pair yields as a function of @xmath6 for like - sign and unlike - sign electron - muon pairs were the same within 3% over all @xmath6 .
while this corroborates the basic idea of the subtraction , the assumption was further tested with data . in the following sections we detail the results of different methods to tag electrons and muons from light flavor decay to examine the validity of eq .
[ eq : subtrcorr ] and to quantify the systematic uncertainty of the method .
the general method is to use a sample of single electrons paired with single muons , where one or both are likely from light - hadron decays .
if the method is correct , the like - sign subtraction should produce no correlation at all .
if there are statistically significant correlations after like - sign subtraction , these are subtracted from the final @xmath21@xmath13 pair yield and uncertainties on the residual correlation strength are propagated as a systematic uncertainty on the final pair yield . if no statistically significant yield is found after like - sign subtraction , the statistical uncertainty on the zero yield is propagated as the systematic uncertainty .
one source of background to the single muons is from hadrons that penetrate to the 5th gap , called punch - through hadrons .
after single particle cuts there is some small fraction ( roughly 1 out of every 250  @xcite ) of candidate tracks with @xmath48  gev/@xmath3 that are hadrons that punch through . while this represents an irreducible background to the single muons , we can obtain a clean sample of hadrons that punch through and stop in the 4th gap of the muid .
[ fig : pz ] shows the @xmath82 distribution of muon candidates that stop in the 4th gap .
the peak at 2.3  gev is composed of muons that have insufficient energy to penetrate further .
the broader portion of the distribution comprises light hadrons that are not stopped by the upstream absorber materials but are subsequently absorbed in the steel just after the 4th gap , thus not leaving a hit in the 5th gap .
we identify punch - through hadrons as having stopped in the 4th gap with @xmath82 larger than 3  gev .     for tracks that stop in the next - to - last muid gap ( 4th gap ) .
the peak at lower @xmath82 is due to muons , while the broad distribution is from hadrons that punch through the absorber to the 4th gap .
the solid line is a two - gaussian fit to this distribution with the solid line indicating the hadronic background in the muon peak region . ]
/ndf = 22.7/24 . in it is a flat line and a gaussian centered at @xmath42 with @xmath83/ndf = 26.3/22.,title="fig : " ] /ndf = 22.7/24 . in it
is a flat line and a gaussian centered at @xmath42 with @xmath83/ndf = 26.3/22.,title="fig : " ]    fig .
[ fig : punch - throughcorr ] shows the fully corrected like - sign subtracted pair yield of central - arm electrons and the punch - through hadrons in the muon arms for both and collisions .
if both the like- and unlike - sign pair yields were dominantly from light hadron decays , the like - sign subtraction should produce zero pair yield . to determine the magnitude of the residual correlation strength after like - sign subtraction
, the data were fitted with a flat line .
this is shown as the solid line in fig .
[ fig : punch - throughcorr]a .
the fit uncertainty is shown as the shaded band around the solid line .
the flat fit in had a @xmath83/ndf of 22.7/24 and gave a value that was nonzero with greater than 1@xmath52 significance .
this means there is yield in the final @xmath21@xmath13 correlations from these punch - through hadrons .
the fitted yield was subtracted from the final pair yield and its uncertainty was propagated as a systematic uncertainty on the final pair yield . for the case , we fitted the residual correlation to a flat line and found reasonable agreement with a @xmath83/ndf of 30.9/24 or a @xmath7-value of 14% .
however , there is a possible excess of counts near @xmath14 , which when included as a gaussian component fixed at @xmath14 and the width and yield as free parameters , a slightly better @xmath83/ndf of 26.3/22 or a @xmath7-value of 26% was found .
if there is any correlated yield beyond a pedestal , it would show up in the back - to - back peak .
therefore , we subtract the gaussian fit , shown as the solid line in fig .  [
fig : punch - throughcorr ] from the final pair yield and propagate the uncertainty , shown as the shaded region around the solid line , on the fit to the systematic uncertainty in the final pair yield .
two additional corrections to this data are applied before subtraction from the final pair yield . because the punch - through hadrons are measured in the 4th gap , the yields need to be scaled to match the rate of hadrons at the last gap .
the rate of hadrons at the 5th gap was determined by using pion and kaon nlo perturbative qcd spectra  @xcite and passing them through a geant-3 model of the phenix muon arms .
the muid absorber steel cross section was modified until there was agreement between data and the simulation for the rate of punch - through hadrons in the 3rd and 4th gap .
we extrapolated to the 5th gap and find the rate of hadrons is 2.81@xmath120.30 times the rate of punch - through hadrons in the 4th gap  @xcite .
the 3  gev/@xmath3 @xmath82 cut removes some fraction of the punch - through hadrons .
based on the two - component fit to the @xmath82 distribution shown in fig .
[ fig : pz ] , the yield is scaled up to account for those hadrons rejected by the @xmath82 cut . in the end ,
the pair yield uncertainty is 2.17@xmath8410@xmath85 ( rad)@xmath46 in . in there is a @xmath6-independent uncertainty on the final pair yield that is 1.42@xmath8410@xmath86 ( rad)@xmath46 and the gaussian uncertainty that ranges from 0 to 6.30@xmath8410@xmath86 ( rad)@xmath46 .
one source of real muons is from decays of light hadrons , predominantly charged pions and kaons , before the absorber material
. the observed rate of muons into the north arm is higher , when the collision vertex is farther from the spectrometer arm . because heavy flavor decays ( including drell - yan , heavy quarkonia , etc . )
have a much shorter @xmath87 than light flavor decays , heavy flavor decay muons have a much weaker vertex dependence .
therefore , we assume there are two components to the muon rate : a component that follows the primary vertex distribution , attributable to heavy flavor decays , and a component that folds the linear component due to light hadron decays with the primary vertex distribution .
pair yield in ( a ) and ( b ) collisions .
both are consistent with no residual correlation after like - sign subtraction .
the solid lines and shaded bands indicate the flat line fits and their uncertainty with @xmath83/ndf of 27.1/24 and 18.0/24 in and  , respectively .
, title="fig : " ]   pair yield in ( a ) and ( b ) collisions
. both are consistent with no residual correlation after like - sign subtraction .
the solid lines and shaded bands indicate the flat line fits and their uncertainty with @xmath83/ndf of 27.1/24 and 18.0/24 in and  , respectively .
, title="fig : " ]    muons that are near the detector ( @xmath88  cm ) and far from the detector ( @xmath89  cm ) , where @xmath44 is the measured collisions vertex , are separately correlated with central arm electrons .
because the signal heavy flavor muons follow the primary collision vertex distribution , subtracting the near - vertex pair yield from the far - vertex pair yield , should remove these and only residual correlations from decay muons should be present .
the pair yields in and after subtracting near- and far - vertex muons and after like - sign subtraction are shown in fig .
[ fig : decaycorr ] .
the correlations are consistent with a flat line with zero yield with a @xmath83/ndf of 18.0/24 .
the data seems to have a residual shape . however , this shape is asymmetric about @xmath6 of zero and is not physical .
therefore , we fit with a flat line that results in zero correlation yield and a @xmath83/ndf of 27.1/24 .
the fits are shown in fig .
[ fig : decaycorr ] as solid lines and shaded bands indicating the statistical uncertainties .
these uncertainties were propagated into the systematic uncertainties of the final pair yields .    to propagate the uncertainties ,
additional corrections are needed .
first , in the far - near subtraction , some fraction of the decay muons are removed .
second , light hadron decays outside the @xmath1230  cm vertex cut are not counted in the subtraction . to account for both effects , a fit to the vertex dependence of the muon yield
is extrapolated to a point one interaction length inside the absorber , a distance of about 56  cm from the nominal @xmath20 vertex and about 16  cm into the absorber .
it is assumed that the decay contribution to the muons is negligible at that point , which fixes the fraction of muons that are from light decays within the measured vertex window of the analysis . under this assumption , only 22% of the decay muons are measured within the vertex window after the like - sign subtraction .
the fit uncertainties are increased to account for those muons not measured .
the final systematic uncertainties on the final pair yield are 1.13@xmath8410@xmath86 ( rad)@xmath46 and 5.05@xmath8410@xmath86 , independent of @xmath6 for and , respectively .
electrons can result from light hadron decays through internal and external photon conversions .
the dominant photonic source of electrons are from @xmath22 decays .
we assume that , if we measure the @xmath22-decay electrons correlations with muons , this will represent the other photonic sources ( such as @xmath51 and @xmath90 decay ) in shape and yield . to tag decay or converted electrons , we construct the invariant mass distribution of all pairs of electrons and photons in an event .
electrons paired with photons within the @xmath22 mass peak are then correlated with muon candidates
. the signal - to - background of pairs in the @xmath22 mass range is about one . to remove correlations from combinatorial electron - photon pairs that fall within the @xmath22 mass window ,
muon candidates were also correlated with the @xmath21@xmath91 pairs in a `` sideband '' @xmath22 mass region from 0.20.4  gev/@xmath92 . after scaling by the appropriate signal - to - background under the @xmath22 mass region ,
the `` sideband '' correlations were subtracted from the in - mass electron - muon correlations for each of the @xmath21@xmath13 charge types .
[ fig : photoniccorr ] shows the `` sideband''-subtracted and like - sign subtracted correlation between electrons tagged in the @xmath22 mass region with muons from and data .
flat fits to these correlations produced a yield consistent with zero with @xmath83/ndf of 33.2/24 and 20.2/24 in and data , respectively .
the statistical uncertainty from the fitted yield to the @xmath22-tagged correlations is a factor of 10 smaller than the other background correlations after accounting for reconstruction efficiency and additional sources of photonic electrons .
this uncertainty is negligible compared to those from the muon backgrounds .
pair yield in ( a ) and ( b ) collisions .
both are consistent with no residual correlation after like - sign subtraction .
the solid lines and shaded bands indicate the flat line fits and their uncertainties with @xmath83/ndf of 33.2/24 and 20.2/24 in and  , respectively .
, title="fig : " ]   pair yield in ( a ) and ( b ) collisions
. both are consistent with no residual correlation after like - sign subtraction .
the solid lines and shaded bands indicate the flat line fits and their uncertainties with @xmath83/ndf of 33.2/24 and 20.2/24 in and  , respectively .
, title="fig : " ]      in this analysis there are three general types of uncertainty that we identify as type a , point - to - point uncorrelated , type b , point - to - point but correlated , and type c , total normalization uncertainty . except for statistical uncertainties
there are no type a uncertainties in this analysis .
the type b uncertainties are from the subtraction of known backgrounds discussed in section  [ sec : subtr ] .
the fully corrected pair yield uncertainties in are @xmath93 ( rad)@xmath46 and @xmath94 ( rad)@xmath46 from punch - through hadron and decay hadron subtraction uncertainties , respectively .
these values are independent of @xmath6 . in the flat - line fit contributions to the systematic uncertainty are @xmath95 ( rad)@xmath46 and @xmath96 ( rad)@xmath46 from punch - through hadron and decay hadron subtraction uncertainties , respectively .
the additional uncertainty from the gaussian fit to the punch - through hadron correlations resulted in a @xmath6-dependent uncertainty ranging in absolute value of 0 at @xmath97 2 rad to @xmath98 ( rad)@xmath46 at @xmath99 .
the type - b systematics are summarized in table  [ tab : sys ] .
cccc type & description & & + b & @xmath6 dependent &  & 0%6.30@xmath8410@xmath86(rad)@xmath46 + b & punch - through & @xmath93 ( rad)@xmath46 & @xmath95(rad)@xmath46 + b & decay muons & @xmath94 ( rad)@xmath46 & @xmath96(rad)@xmath46 +   + c & muon cuts & 7.8% & 8.3% + c & electron cuts & 8.3% & 9.3% + c & muon efficiency & 2.2% & 2.2% + c & electron efficiency & 1.0% & 1.0% + c & trigger efficiency & 11.1% & 4.2% +   +
c & total & 16.1% & 13.4% +    the type c uncertainties are attributable to several sources and are given in table  [ tab : sys ] .
one source of systematic uncertainty is evaluated by tightening the single particle cuts for this analysis .
each single particle cut was tightened independently and the analysis , including reevaluation of the single particle efficiency , was performed .
the uncertainty from each of the individual single particle cuts was combined using the correlation amongst the cuts .
the values of these are different in and data , because of the higher backgrounds in collisions .
however , these uncertainties are highly correlated between and , because the same cuts are applied to both data sets .
another source of uncertainty is in the evaluation of the single particle efficiencies .
the single particles were generated flat in @xmath28 and then weighted to match the measured phenix heavy flavor lepton spectra  @xcite . for the uncertainty determination ,
the single particle efficiency was re - evaluated without the weighting applied .
this was estimated to be 1.0% for the electrons and 0.8% for the muons . for muons
there is an additional 2.0% uncertainty due to the run - by - run variation in muon acceptance .
the final portion of the type c systematic uncertainty is due to the trigger efficiency . to evaluate this uncertainty ,
the data were analyzed for several data - taking periods defined by the muon trigger performance .
the difference in fully - corrected yields between data sets was taken to be the uncertainty in the muon trigger efficiency .
this is combined with the uncertainties in the bias factor @xmath3 in eq .
[ eq : pairyield ] .
the total uncertainty for the trigger is 11.1% for and 4.2% for .
as indicated in table  [ tab : sys ] , combining all type c uncertainties gives 16.1% for data and
13.4% for data .
the fully - corrected like - sign subtracted @xmath21@xmath13 pair yield as a function of @xmath6 for electrons with @xmath2  gev/@xmath3 and @xmath100 , with opposite - signed forward muons with @xmath5  gev/@xmath3 and @xmath4 , in is shown in fig .
[ fig : ppcorr ] .
the average muon @xmath51 in these correlations is 1.75 .
the error bars are statistical uncertainties only , while the boxes are the type b systematic uncertainties .
we note that the distribution has two components : a nonzero continuum as well as a back - to - back peak near @xmath101 .
@xmath13 pair yield in .
the error bars are statistical only .
the boxes show the type b systematic uncertainty from the punch - through hadron and light hadron decay muon background subtraction .
the 16.1% type c systematic uncertainty is not shown . ]
points ) with heavy flavor production in powheg ( [ blue ] dashed line ) , pythia ( [ black ] solid line ) and mc@nlo ( [ green ] long dashed line ) . the @xmath21@xmath13 pair yield from the subset of pythia events , when the @xmath17 is not produced at the event vertex
is plotted as the dotted [ black ] line .
each monte carlo curve was scaled by a single parameter to match the observed yield .
the resulting cross sections are consistent with the previously measured phenix results ( see table  [ tab : xsect ] ) . ]    to interpret these data , we compare the results to several different monte carlo generators , pythia , powheg@xcite , and mc@nlo@xcite .
the pythia mb qcd events were generated to model the lo gluon fusion process and also model next - to - leading order processes , like flavor excitation and gluon splitting .
events with a @xmath17 pair and an electron and a muon in the measured kinematic range as the corrected data ( @xmath102  gev/@xmath3 and @xmath103 , @xmath104  gev/@xmath3 and @xmath105 ) were correlated and a like - sign subtraction was performed .
an overall scale factor was used to fit the pythia curve to the data . in the fit
, the @xmath83 was calculated for different scale parameters using the statistical error on the data .
we report the cross section for the scale factor that minimizes that @xmath83 and report a statistical error on the cross section as the value that changes the @xmath83 by one unit . to evaluate the systematic uncertainty on the cross section , the data were increased and decreased by their combined type b and type c systematic uncertainty and the process to determine the scale factor by finding a minimum @xmath83 using the statistical uncertainty in the data was repeated .
we find the pythia correlation is consistent with the data with a @xmath17 cross section of @xmath106 = 340@xmath1229(stat)@xmath12116(syst )  @xmath13b with a @xmath83/ndf of 20.5/24 .
this is shown as the solid curve in fig .
[ fig : compare ] .
the other model comparisons are from nlo generators , powheg and mc@nlo .
events were generated to produce the hard scattering heavy flavor event vertex and then interfaced to pythia , which performed the fragmentation and underlying event generation .
the qualitative features of the data are present in these correlations : the continuum and the back - to - back peak .
as described for the pythia fit , a single scale parameter was used to calculate a @xmath83 between the generated @xmath21@xmath13 correlations and the data using the data s statistical uncertainty .
the resulting best fits for powheg and mc@nlo are shown in fig .
[ fig : compare ] as the short dashed and the long dashed lines , respectively .
the extracted cross sections are @xmath106 = 511 @xmath12 44 ( stat )  @xmath12 198 ( syst )  @xmath13b with @xmath83/ndf of 23.5/24 for powheg and @xmath106 = 764 @xmath12 64 ( stat )  @xmath12 284 ( syst )  @xmath13b with @xmath83/ndf of 19.2/24 for mc@nlo .
we combine the cross sections from the three models and report a measured cross section of @xmath11 538 @xmath12 46 ( stat ) @xmath12 197 ( data  syst ) @xmath12 174 ( model  syst ) .
the central value of the cross section is the average of the three model cross sections , while the model systematic uncertainty is the standard deviation of the three model cross sections .
this value can be compared with previous phenix measurements . from the heavy flavor electron spectra at midrapidity ,
phenix found @xmath11 567 @xmath12 57 ( stat ) @xmath12 224 ( syst )  @xcite and from the dielectron mass spectrum at midrapidity , phenix extracted @xmath11 554 @xmath12 39(stat ) @xmath12 142  ( data  syst )  @xmath12 200  ( model  syst )  @xcite . within the data systematics
the value extracted here is consistent with previously published phenix results .
cccc & description & @xmath106 ( @xmath13b ) & + & pythia @xmath21@xmath13 & 340@xmath1229(stat)@xmath12116(syst ) & + & powheg @xmath21@xmath13 & 511@xmath1244(stat)@xmath12198(syst ) & + & mc@nlo @xmath21@xmath13 & 764@xmath1264(stat)@xmath12284(syst ) & + & combined @xmath21@xmath13 & 538@xmath1246(stat)@xmath12197(data syst)@xmath12174(model syst ) & +   + & phenix single @xmath107  @xcite & 567@xmath1257(stat)@xmath12224(syst ) & +   + & phenix dilepton ( @xmath26 )  @xcite & 554@xmath1239(stat)@xmath12142(data syst)@xmath12200(model syst ) & +    [ tab : xsect ]    using the pythia event record , it is possible to separate the @xmath17 production into an lo component , where the @xmath108 and a component from the pythia model of nlo mechanisms of flavor excitation and gluon splitting , where the @xmath17 pair is produced in the initial or final - state shower .
the `` pythia ( no lo ) '' dashed line in fig .
[ fig : compare ] shows the correlations from the sample of produced pythia events , where the @xmath17 were not generated in the primary event vertex of pythia .
the back - to - back peak at @xmath14 is dominated by the lo gluon fusion process while the continuum is due to the correlations from the higher order processes . from an accounting from pythia
, we find that 32% of the @xmath21@xmath13 pair yield results from gluon fusion , consistent with the expectations from charm production  @xcite .    throughout the analysis
it has been assumed that semileptonic @xmath17 decay is the dominant contribution to the correlations .
however , @xmath75 semileptonic decays would produce a signal in both the like- and the unlike - sign pair distributions .
up to four semileptonic decays can occur where @xmath109-quarks semileptonically decay to @xmath3-quarks , which subsequently semileptonically decay .
we have used pythia and powheg to check these contribution from bottom . in both cases , for electrons and muons in the kinematic region that we measure ,
the bottom contribution is about a factor of 100 below the charm yield .
this is further corroborated by the phenix heavy flavor electron measurements that show that bottom becomes significant only at @xmath28 above 3 gev/@xmath3  @xcite . in this analysis only 3% of the sampled electrons
have a @xmath28 above 3  gev/@xmath3 , so we expect that the bottom contribution is negligible in this measurement especially compared to the background subtraction systematic uncertainties
.      @xmath13 pair yield in .
the error bars are statistical only .
the boxes show the type b systematic uncertainty from the punch - through hadron and light hadron decay muon background subtraction .
the 13.4% type c systematic uncertainty is not shown . ]
the fully - corrected like - sign subtracted pair yield as a function of @xmath6 for electrons with @xmath102  gev/@xmath3 and @xmath103 with forward muons with @xmath111  gev/@xmath3 and @xmath105 in 0%100% , corresponding to the total inelastic cross section , is shown in fig .
[ fig : daucorr ] . a nonzero correlations strength is observed .
however , unlike the data , there is a much less distinct back - to - back peak near @xmath6 of @xmath42 .
[ fig : paircorr ] shows the overlay of the and pair correlations .
the pair correlations are scaled by the @xmath112  @xcite .
the peak in is suppressed compared to , indicating a medium modification to the yield per collision in .
@xmath13 pair yield in ( [ red ] circles ) @xmath113-scaled ( [ blue ] boxes ) , shifted in @xmath6 for clarity .
the bars are statistical uncertainty .
the boxes are the type b systematic uncertainty from the decay and punch - through background subtraction .
the overall normalization uncertainties of 16.1% and 13.4% in and , respectively , and 5.7% uncertainty from @xmath114 are not included . ]
is plotted as a function of @xmath6 .
the vertical bars are statistical uncertainties , the black boxes are the type b systematic uncertainties , and the gray band around 1.0 on the left is the type c systematic uncertainty .
the type b systematics are symmetric around the central value but in some cases are outside the range of the plot .
( b ) alternative rebinning of the data in ( a ) .
, title="fig : " ]   is plotted as a function of @xmath6 .
the vertical bars are statistical uncertainties , the black boxes are the type b systematic uncertainties , and the gray band around 1.0 on the left is the type c systematic uncertainty .
the type b systematics are symmetric around the central value but in some cases are outside the range of the plot .
( b ) alternative rebinning of the data in ( a ) .
, title="fig : " ]    to quantify the difference between and yields , we calculate the ratio @xmath16 defined as the ratio of a pair yield in to the @xmath114-scaled pair yield in . @xmath115
any deviation from unity of this ratio would indicate modification to the yield .
when taking this ratio several systematic uncertainties common to and cancel .
these are dominantly from identical cuts used in the analyses with the same systematic uncertainties .
the noncanceling type c systematic uncertainties in the and yields are 7.7% and 8.9% , respectively .
[ fig : jdadphi]a shows a plot of @xmath16 as a function of @xmath6 for all bins in @xmath6 .
the bars are statistical uncertainties and the type b systematic uncertainties are plotted as boxes . the noncanceling type - c uncertainty is 14.1% and is indicated by the shaded box around one on the left . while the points near @xmath116 are consistent with unity with large error bars , the points near @xmath14 , where @xmath16 is about 0.4 . for clarity ,
[ fig : jdadphi]b shows a rebinning of fig .
[ fig : jdadphi]a .
we find @xmath117 for the bin near @xmath14 .
this value is 3.5@xmath52 different from unity after combining the statistical and systematic uncertainties .
these results show that , in the measured kinematics , charm pairs are modified in the cold nuclear medium .
these results are in a different kinematic region than either the single electrons , which are enhanced at midrapidity  @xcite , or the single muons , which are suppressed at forward rapidity  @xcite . from the pythia simulation ,
the @xmath21@xmath13 correlations arise from partons in the gold nucleus with @xmath118 at @xmath119 , on the edge of the shadowing region . as discussed in section  [ sec : yields ] , the back - to - back peak is dominated by leading order gluon fusion , while the continuum is dominated by other processes like flavor excitation and gluon splitting .
the observed back - to - back peak and pedestal in and should help lead to an understanding of the mechanism or mechanisms responsible for the modification .
for example , the back - to - back peak is dominated by low-@xmath25 gluons participating in the hard scattering , whereas the continuum has a larger contribution of quarks participating in the hard scattering .
quarks are probably less shadowed than gluons at the @xmath25 and @xmath120 where this analysis is measured .
it is possible that there are kinematic differences between the final state charm quarks in the peak and the continuum .
these differences could affect the amount of final state energy loss and multiple scattering that modify the measured pair yields .
it may be possible to combine these results with other cold nuclear matter charm measurements to disentangle the effects of shadowing , saturation , and energy loss .
we presented phenix results for heavy flavor production of azimuthally - correlated unlike sign @xmath21@xmath13 pairs in and collisions at @xmath121 of 200  gev .
the yield shows a nonzero continuum as well as a back - to - back peak structure centered at @xmath101 .
when compared with several models , we find the charm cross section @xmath11 538 @xmath12 46 ( stat ) @xmath12 197 ( data  syst ) @xmath12 174 ( model  syst ) @xmath13b .
this is also consistent with previously measured @xmath17 cross sections at this center of mass energy . in collisions a yield reduction in the back - to - back peak is observed , where we measure @xmath122 = 0.433 @xmath12 0.087 ( stat ) @xmath12 0.135 ( syst ) .
this indicates that the nuclear medium modifies the @xmath17 correlations .
such a suppression could arise due to nuclear pdf shadowing , saturation of the gluon wavefunction in the au nucleus , or initial / final state energy loss and multiple scattering .
we thank the staff of the collider - accelerator and physics departments at brookhaven national laboratory and the staff of the other phenix participating institutions for their vital contributions .
we acknowledge support from the office of nuclear physics in the office of science of the department of energy , the national science foundation , a sponsored research grant from renaissance technologies llc , abilene christian university research council , research foundation of suny , and dean of the college of arts and sciences , vanderbilt university ( u.s.a ) , ministry of education , culture , sports , science , and technology and the japan society for the promotion of science ( japan ) , conselho nacional de desenvolvimento cientfico e tecnolgico and fundao de amparo  pesquisa do estado de so paulo ( brazil ) , natural science foundation of china ( p.  r.  china ) , ministry of education , youth and sports ( czech republic ) , centre national de la recherche scientifique , commissariat  lnergie atomique , and institut national de physique nuclaire et de physique des particules ( france ) , bundesministerium fr bildung und forschung , deutscher akademischer austausch dienst , and alexander von humboldt stiftung ( germany ) , hungarian national science fund , otka ( hungary ) , department of atomic energy and department of science and technology ( india ) , israel science foundation ( israel ) , national research foundation and wcu program of the ministry education science and technology ( korea ) , physics department , lahore university of management sciences ( pakistan ) , ministry of education and science , russian academy of sciences , federal agency of atomic energy ( russia ) , vr and wallenberg foundation ( sweden ) , the u.s .
civilian research and development foundation for the independent states of the former soviet union , the us - hungarian fulbright foundation for educational exchange , and the us - israel binational science foundation .