Spinor helicity methods in high-energy factorization: efficient momentum-space calculations in the Color Glass Condensate formalism

We use the spinor helicity formalism to calculate the cross section for production of three partons of a given polarization in Deep Inelastic Scattering (DIS) off proton and nucleus targets at small Bjorken x. The target proton or nucleus is treated as a classical color field (shock wave) from which the produced partons scatter multiple times. We reported our result for the final expression for the production cross section and studied the azimuthal angular correlations of the produced partons in [1]. Here we provide the full details of the calculation of the production cross section using the spinor helicity methods.


Introduction
The Color Glass Condensate (CGC) formalism (see [2] for a review) is an effective field theory approach to QCD at small x where the gluon density in a proton or nucleus is expected to be so large that the standard QCD-improved parton model must break down. In this limit it is more appropriate to treat the hadron or nucleus as a coherent color field rather than a collection of incoherent and individual partons. At the classical level, the CGC generalizes scattering via exchange of a single gluon to multiple gluon exchanges (high energy factorization). By including quantum loop effects, it resums the large logs of Bjorken x [3] (or equivalently energy) as opposed to the leading twist collinear factorization formalism which resums large logs of the hard scale Q 2 . The Color Glass Condensate formalism is the appropriate one to use in the kinematic limit where log 1/x log Q 2 . This is typically the case for particle production in the low to intermediate p t regime and in the forward rapidity region [4] where one probes the small x components of the target wave function.
The Color Glass Condensate formalism has been applied to various high energy processes. A first class of observables is associated with the forward hadron production at RHIC and the LHC (see [5] for a compilation of LHC forward physics observables and [6] for a compilation of the CGC-motivated phenomenology of particle production). While there are strong hints for the presence of gluon saturation effects in the observed suppression of forward hadron production cross section at RHIC and the LHC [7], the effect is far from being universally accepted as coming only from saturation. Models which modify the parton distribution functions of nuclei and add energy loss effects have also been used to fit the data [8]. Azimuthal angular correlations of produced di-hadrons have also been measured (in deuteron-gold collisions) and show a strong decorrelation of the azimuthal angle between the two produced hadrons in the forward rapidity region as predicted in the CGC formalism [9]. The away side peak is observed to be broader which can be understood to be due to multiple scattering. Furthermore, the magnitude of the peak is reduced which is due to the small x evolution of the target wave function which causes "nuclear shadowing". Various angular correlations have also been suggested as possible probes of gluon saturation in the target [10,11]. They all exhibit the same qualitative features in the CGC formalism which will be further probed by the current experiments at the LHC and in the proposed Electron-Ion Collider [12].
While much can be learned from hadron-hadron collisions, observables generally suffer from contaminations due to soft and collinear radiation. A much cleaner environmentboth experimentally and theoretically -is, on the other hand, provided by photon induced processes. While much has been learned from the study of proton structure functions measured at the HERA experiment, ultra-peripheral collisions (UPCs) of protons on protons and nuclei on protons at the LHC, allow to study the currently most energetic photon-proton collisions, where either a large nucleus or a proton serve as the photon source. Among the most prominent observables explored so far is exclusive photo-production process of vector mesons such as the J/Ψ [13,14]. Such processes are of particular interest since they allow to probe the gluon distribution in the proton down to very small values of Bjorken x ∼ 10 −5 , see [15] for CGC-related studies. While such CGC studies provide an excellent description of the energy dependence of the exclusive J/Ψ photo-production cross-section, the same energy dependence has also an excellent description in terms of (dilute) next-to-leading order BFKL evolution (see e.g. [16]) which, unlike the CGC-formalism, does not include corrections due to high densities. From a theoretical point of view the difficulty to distinguish low x evolution without (BFKL) and with (BK, JIMWLK) is related to the type of correlators probed by such observables: both proton structure functions and the amplitude for exclusive photo-production of vector mesons are directly proportional to the two-point correlator of two Wilson lines, which itself can be related to an unintegrated gluon density. The description of such observables is therefore at first identical for low (BFKL) and high (BK, JIMWLK) densities and the presence of high density effects only manifests itself through the particular form of low x evolution. Finding definite evidence for saturation effects requires therefore evolution of the observables far into the low x region, well beyond the kinematical reach of current collider experiments.
In [1] we proposed triple jet/hadron production in DIS as an excellent probe of gluon saturation dynamics in a hadron or nucleus. Unlike the observables discussed above, the theoretical description of observables with multiple final states involves higher order correlators of Wilson lines in the target. As a consequence high density effects manifest themselves directly at the level of the observable without the need to invoke low x evolution. In particular there are substantial differences between a BFKL/low density and a CGC/high density description. While the measurement of such high multiplicity final states is much more cumbersome experimentally, such observables provide important complementary information to inclusive observables, and therefore allow to pin down the valid description of perturbative QCD in the low x region. The case of three jet/hadron production is then of particular interest, since such an observable provides two relative angles rather than one in the case of two-particle correlations, so there is an additional knob to turn. Furthermore there are two away side peaks, each of which will sensitively depend on saturation dynamics and small x evolution. We showed our final expression for the triple differential cross section in [1] and used it to study, quantitatively, the angular dependence of the cross section in certain kinematics. The present paper is dedicated to a presentation of the full details of the calculation, employing spinor helicity formalism (for introductory reviews see [17]) which leads to a major simplification of the calculation. For another recent application of the spinor helicity formalism in the framework of light-cone wave function see [18]; for a related calculation of 3 parton impact factor in DIS in the shock wave framework see also [19,20]; see also the calculation of next-to-leading corrections to inclusive DIS of [21,22]. This paper is organized in the following way: we write down the formal expressions for the production amplitude in Section 2 and proceed to study the space-time structure of the diagrams in Section 2.1 and 2.2. Various cuts of the diagrams are considered and shown to lead to vanishing of some diagrams due to the "wrong cut". This allows us to re-write the expressions for the amplitude in a more compact way. We then give an overview of the spinor helicity formalism in Section 3 and proceed to apply it to the process considered. Major simplifications occur when one considers scattering and production of partons of a fixed helicity due to helicity conservation. We summarize in Section 4 and, using crossing symmetry, point out the connections between our results and Multi-Parton Scattering (MPI) in the projectile, in photon-jet or photon production in forward rapidity proton-nucleus collisions.

CGC amplitudes in momentum space
We consider production of three partons in DIS as depicted in Fig. 1, where the target can be a proton or nucleus and where Q 2 = −l 2 denotes the virtuality of the photon. We will study the process in the limit of high center of mass energy √ s → ∞ with s = (l+P ) 2 . To this end it is convenient to introduce light-cone vectors n,n which are defined through the four momenta of virtual photon and target. With the Sudakov decomposition of we obtain for the momenta of initial particles To include the possibility of arbitrary large gluon densities in the target, we represent the latter by its gluonic field which can reach a maximum strength of A µ ∼ 1/g, with g the gauge coupling. To calculate scattering amplitudes in the high energy limit it is then convenient to treat the gluon field of the target as a background field (shock-wave); in light-cone gauge A · n = 0, the only non-zero component is A − (x + , x t ) = δ(x + )α(x t ), while A t = 0 in the high energy limit. Amplitudes are written in terms of momentum space quark and gluon propagators in the presence of the background field, see e.g. [23], which are directly obtained from Fourier transforming their corresponding counter parts in configuration space. In the above we use the conventional free fermion and gluon propagator, and where with Wilson lines in fundamental (V ) and adjoint (U ) representation. They read with −iT c ab = f acb . To construct amplitudes in the presence of a (strong) background field, it is convenient to extend conventional QCD momentum space Feynman rules by two additional rules: (a) adding the vertices Eqs. (7) and (8) and (b) the requirement that all internal momenta p, i.e. momenta which cannot be expressed in terms of momenta of external particles, are integrated over with the measure d 4 p (2π) 4 , in 1-1 correspondence to conventional loop momenta. In tree diagrams such internal momenta arise if n ≥ 2 vertices from Eqs. (7) and (8), are inserted into a single Feynman diagram; see Fig. 2 for an illustrative example. If the number n of produced colored particles in the final state is small, n ≤ 2, the above method provides an efficient alternative to the calculation of matrix elements in the presence of large gluon densities, see [11,24] for earlier examples. For final states with large multiplicities, n ≥ 3, the method becomes inefficient due to the large number of Feynman diagrams which need to be considered. While the process γ * + target → q +q requires 3 diagrams, one finds already 16 diagrams for the process γ * + target → q +q + g. Moreover, calculations based on configuration space propagators suggest a large redundancy in the result, leading to large cancellations among different diagrams, which further complicate the use of momentum space methods. In the following we demonstrate that there is a direct way to use the emerging simplicity of the configuration space amplitudes for calculations in momentum space.

Cut diagrams according to negative and positive x +
To achieve the required reduction of Feynman diagrams it is necessary to exploit the restriction of the interaction between projectile and target to the light-cone time-slice x + = 0, A − (z) ∼ δ(z + ). Such properties are naturally revealed in the configuration space and/or light-front formalism, which is the main reason for their preferred use over momentum space techniques. To make these properties and the associated reduction of diagrams explicit in momentum space, it is sufficient to associate with each propagator a separate minus-momentum variable. This can be achieved by introducing at each (standard QCD) vertex (where minusmomenta are conserved i.e., excluding the vertices in Eqs. (7) and (8)), a delta function together with a corresponding integration over the newly introduced minus momentum variable; here {p in (out) } denotes the full set of incoming (outgoing) momenta of the vertex. Expressing each of these delta functions as a Fourier integral, we finally obtain for each (standard QCD) vertex an integral over x + i , i = 1, . . . N with N the total number of (standard QCD) vertices. These vertices are then connected by momentum space propagators which are Fourier transformed w.r.t. their plus momentum 1 . At the vertices defined in Eqs. (7), (8) the minus momentum is not conserved and therefore no integration over minus coordinates appears; instead these vertices are associated with the 'light-cone time' x + = 0, i.e. we use e ix + p − x + =0 = 1. Fourier transforming quark and gluon propagators w.r.t. to their plus momenta one finds easily the well known expressions Here Using such propagators together with integrations over light-cone times x + i at each vertex one re-obtains directly the light-cone time ordered Feynman rules of lightfront perturbation theory (old-fashioned perturbation theory), see e.g. [22] which contains Figure 3: Real corrections. We have l + , p + , q + , k + > 0 and l + = p + + k + + q + a recent and compact overview. In particular we find that it is natural to organize QCD interaction vertices into vertices before (x + i < 0) and after the interaction (x + i > 0) with the target field. In a nut-shell this result will allow us to evaluate to zero configurations which are absent in the light-front formalism for momentum space diagrams, simplifying enormously the reduction of diagrams at earlier stages of the calculation.
To be more precise, we first note that potential contact terms in the quark ( ∼ δ(x + ij )/ n) and gluon (∼ δ(x + ij )n µ n ν ) propagators are absent, if the propagators connect with the vertices Eqs. (7) and (8): for the quark due to the identity / n/ n = 0, for the gluon since d µν (p) · n ν = 0 and n 2 = 0 (if the vertex connects to a virtual gluon) and (λ) (p) ν · n ν = 0 (if the vertex connects to a real gluon). Due to the absence of such contact terms it will be sufficient to limit the discussion to scalar skeleton diagrams, independent of the precise nature of particles. We further note that plus momenta are conserved at all vertices (including Eqs. (7) and (8)). On the other hand, since plus coordinates and plus momenta are the only components which receive special constraints from the propagators in the representation of Eq. (11), it is sufficient to study only their effect. The resulting diagrams are depicted in Fig. 3, where for the time being we ignore the presence of potential target interaction vertices Eqs. (7) and (8). In a next step we divide all integrals over light-cone times x + i into a positive x + i > 0 and a negative x + i < 0 sector. While this is motivated by the restriction of the target interaction vertices to light-cone time x + i = 0, this is irrelevant for the following reason. Using now that propagation from a vertex with negative light-cone time x + j < 0 to a positive light-cone time x + i > 0 implies x + ij > 0 and therefore positive light-cone momentum p + > 0 (due to the structure of theta-functions in Eq. (11)), it is straight forward to verify that each skeleton diagram can be organized according to three possible "s-channel" cuts, see Fig. 4. It is easily possible to extend our result to the case of n final state particles: crossing from negative to positive light-cone times requires positive light-cone momentum and -vice versa -crossing from positive to negative light-cone times requires negative light-cone momentum. In the case of a tree skeleton diagram, all light-cone momenta are always positive and therefore crossing from positive to negative light-cone times is impossible. Moreover, once a certain line crosses from the negative to the positive light-cone time sector, all its daughter lines, i.e., lines emerging from a splitting of this line, will have positive light-cone times by default. The resulting picture tells us that each individual diagram is characterized by all possible  The benefit of this result for the study of the interaction with the target field should be apparent by now: since vertices, Eqs. (7) and (8), are limited to light-cone time x + i = 0, such vertices can only be inserted in a "cut" line. Insertions in un-cut lines, are immediately zero, see Fig. 5 for two configurations which cannot occur. Note that this result applies separately for each individual Feynman diagram and holds regardless of whether the actual evaluation takes place in configuration, momentum or mixed i.e. light-front space. While this leads already to a significant reduction in the number of Feynman diagrams to be evaluated, the remaining set of diagrams still contains a sizable fraction of redundant contributions. In particular there are large cancellations between diagrams where we place a target interaction vertex Eqs. (7), (8) at the z + i = 0 cut and diagrams where such interaction is absent, see Fig. 6 for an example.

Further reduction of diagrams
To reduce the number of diagrams further, it is necessary to study the Fourier transform of the complete propagators (containing both interacting and non-interacting) parts, Eq. (4). Using these propagators for all lines in diagrams (for external lines this implies the use of the LSZ-reduction procedure, see e.g. [11]) the full process can be represented in terms of two diagrams only. On the level of scalar skeleton graphs these are precisely those of Fig. 3. It is well known that such complete propagators can be written in configuration space in three parts: one term is associated with crossing from negative to positive light-cone time x + and is directly proportional to either a Wilson line (positive light-cone momentum fraction) or a hermitian conjugate Wilson line (negative light-cone momentum fraction). The other two terms describe free propagation between two points with either negative or positive light-cone time. For the following discussion the following form is sufficient: are given by Eq. (11). The second term necessarily involves crossing the slice x + = 0 and is given bỹ where the condition p + > 0 (p + < 0) selects automatically the configuration x + > 0 > y + (x + < 0 < y + ). Similarly one has in the case of gluons G ab,(U,U † ) µν (x + , y + ; p + , p t ; q + , q t ) =G (0),ac µα (x + , p + , p t )τ cd G (p, q)G (0),db where we defined Applying now these results to the (skeleton) diagrams, Fig. 4, we now find that each "cut" line must necessarily come with a vertex, Eqs. (15) and (16). The "un-cut" lines can only come with a free propagator. At first these propagators are limited to positive or negative light-cone time only and therefore cannot be directly related to their momentum space counter-parts.
Due to the results of Sec. 2.1, adding a free propagator which crosses from negative to positive light-cone time will however give only a zero contribution. Adding such a zero contribution it is then straight forward to Fourier transform our result back to momentum space. Therefore the complete amplitude can be calculated from the six diagrams of Fig. 4 with a vertex, Eqs. (15) and (16) at each cut, minus the diagram with the background field A + set to zero; if the initial particle is is not colored (as in our case) the total number of diagrams reduces finally to four.

The minimal set of amplitudes
The final set of diagrams which we need to evaluate for the process described by Eq. (1) is depicted in Fig. 7. The four matrix elements corresponding to Fig. 7 read where ν (l), * µ (k) denote polarization vectors of the incoming virtual photon and the outgoing gluon respectively and With flux factor F = 2l + , photon momentum fractions {z 1 , z 2 , z 3 } = {p + /l + , q + /l + , k + /l + }, and the three-particle phase space, the differential 3-parton production cross-section reads where . . . A + denotes the average over background field configurations, F = 2l + and

Spinor helicity methods
While spinor helicity methods are well established for calculations within conventional perturbative QCD calculations, its application to calculations in the QCD high energy/Regge limit are rather limited, see [25] for some examples. We therefore start this section by recalling basic definitions which will further serve to fix our notation.

Basic definitions
The presentation in this paragraph follows closely those of the reviews [17]. For massless fermions, helicity is a good, i.e. conserved, quantum number. One defines (on-shell) helicity eigenstates (spinors) as It is further convenient to introduce the following short-hands Which allows to define the basic spinor products by Further note that To evaluate such spinor products it is necessary to pick a certain representation of the Dirac γ matrices. In the Dirac representation with σ i , i = 1, . . . , 3 the Pauli-matrices. Using the convention of Eq. (2) to define light-cone momenta, the massless spinors can be written as follows, Using these expressions it is possible to obtain explicit formulae for the spinor brackets All of the above versions are equivalent to each other, but prove useful in specific circumstances. We further note that where the last relation is particularly useful to invert helicities in an amplitude. While the above relations are valid in general, they are particularly useful for the description of high energy factorized amplitudes where light-cone momenta k + i are conserved. In particular this enables us to deal with explicit expressions of brackets in the evaluation and leads to a further simplification of diagrams. For instance, for brackets involving light-cone vectors n,n the following relations hold Since the (physical) Hilbert space of a massless vector is isomorphic to the Hilbert space of a massless spinor, it is further possible to express gluon and photon polarization vectors in terms of the above spinors. One has where k denotes the on-shell four-momentum of the boson and n the gauge vector. Obviously one has · k = 0 = · n. These polarization vectors obey the following polarization sum Using the Fierz identity, polarization vectors contracted with Dirac matrices can be written as While the treatment of the real gluon is straight-forward within this formalism, the case of the virtual gluon is at first special. Since l 2 = −Q 2 = 0 in general, the photon has two transverse polarizations which are both annihilated by the light-cone vectors n,n, i.e. one has µ λ=T · n = 0 = µ λ=T ·n. It is therefore straight forward to define transverse polarization vectors within the spinor helicity as where g t µν has only transverse components. The longitudinal polarization vector is given by In the actual calculation we will make use of QED gauge invariance of the scattering amplitude and define¯ While adding a term proportional to l µ will give a zero contribution upon contraction with the amplitude, the above form will facilitate actual calculations.

Calculation of amplitudes
To evaluate the set of amplitudes Eq. (17) we first note that due to the usual spin sums of spinors we have for any on-shell momentum p 2 = 0 which allows us to expand contractions of gamma matrices with momenta p, q, k of the final state particles as well as with the light-cone vectors n,n in terms of spinors. For the momentum of the initial state off-shell photon we use the representation in Eq. (3) and apply Eq. (39) to each of the light-cone vectors. For the internal momentum k 1,2 , which are integrated over and which are therefore intrinsically off-shell, we use the following decomposition sincek 2 1,2 = 0, both vectors used to represent k 1,2 have a spinor representation. While there are many possibilities to present a massive four vector in terms of two mass-less four-vectors, the above presentation has two advantages that make it special: a) since plus-momenta are conserved during interaction with the background field, the denominators introduced in Eq. (40) can be directly expressed in terms of plus momenta of external particles b) making use of Eq. (25) and the fact that Eq. (15) is directly proportional to / n, terms proportional to n µ will be set to zero in most cases. To keep our notation more compact, we further use In the case of gluons a similar simplifications can be achieved using the decomposition of Eq. (40), since For the current calculation the gluon polarization tensor will be always contracted with an external polarization vector. In this case one finds where we made use of the polarization sum, Eq. (33), and Equipped with this information, evaluation of the numerators of the diagrams Eq. (17) is straight forward. While a direct evaluation is possible -most efficiently done using Computer algebra systems -the spinor helicity formalism allows to break down the calculation into independent sub-amplitudes, which allows for a very economic evaluation. In the following we explain this in full detail using diagram 1 and 3 as examples; the case of diagram 2 and 4 follows then directly.

Evaluation of diagrams 1 and 3
To put all the relations above to work and get some practice with the spinor helicity methods, we consider at first the first and third diagrams given by iA 1 and iA 3 . We have where N 1 is the numerator which further includes the denominator of the external quark line Here λ γ = (L, T = ±), λ q = (±), λq = (±), λ g = (±) encode the possible helicity states of initial and final state particles. In the case of amplitude 3 we find with N 3 defined as N λγ ;λqλqλg 3 It is now possible to factorize both N 1 and N 3 which in turn will allow for their efficient evaluation. To this end we split N 1,3 into two, using the spinor representation of / k 1 , which itself relies on the decomposition in Eq. (40). Note that in this decomposition the term proportional to / n gives a zero result for N 1 . We therefore have 2 where Q λγ λq We start with the splitting γ * → qq, with configuration where the incoming photon has longitudinal polarization (λ = L). We find where me made use of Eq. (30) to invert the helicities in the last expression. For transverse photon polarizations T = ± we find For the splitting q → qg it is obviously necessary to distinguish between amplitude one (interaction before the splitting) and three (interaction after the splitting). For the first case we obtain where we used (p + k) 2 = [pk] kp . Note that these splittings allow for a straight forward identification (and if needed isolation) of collinear and soft singularities. Indeed, these splittings encode (together with the correspondingq →qg splittings in amplitude 2) the complete soft (k → 0) and collinear (k ∝ p) singularities of the process, corresponding to the vanishing of the brackets p ± |k ∓ . The corresponding splitting which occurs in the third amplitude is on the other hand free of such potentially singular configurations. One finds 3

Diagrams 2 and 4 from symmetry
The corresponding expressions for diagrams 2 and 4 read: where N 2 is the numerator which further includes the denominator of the external anti-quark line and with N 4 defined as Expanding N 2,4 in terms of spinors one realizes immediately that the resulting expression is identical to N 1,3 up to a) interchange of momenta p ↔ q b)the order in which spinors are written and c) a minus sign for each propagator numerator. In particular one finds the following factorization of N 2 and N 4 : whereQ λγ λq Evaluating these expressions one realizes, that they are identical to the corresponding expressions of diagrams 1 and 3, up to q ↔ p and an overall minus sign (due to the different sign in the numerator of the internal propagators). In case of theq →qg splitting of diagram 2 we find for instance Using now anti-symmetry of the brackets in the numerator as well as the bracket [kq] in the denominator, this expression can be transformed into 4 We therefore find that the corresponding expressions of diagrams 2 and 4 can be obtained from the corresponding ones of diagram 1 and 3 through a) inverting the quark helicity b) interchanging momenta p and q and c) an overall minus sign. In detail we havē

Wilson lines and color algebra
To this end we recall that the amplitudes in Eq. (17)  where # = 1, . . . , 4. To write down the differential cross-section we are in general dealing with expressions of the form with indices a, i, j summed over. It is now convenient to rewrite this as which allows us to define amplitudes with four fundamental color indices only. Extracting factors of Wilson lines and SU(N c ) generators, we find in this way where the first two lines corresponds to diagram one and two respectively while the last line gives the corresponding factor of diagrams three and four. To determine the operators of Wilson lines at cross-section level, we restrict at first to the leading N c terms and subtract contributions without target interaction (as indicated in Eq. (20)). Extracting an overall factor N 2 c /2 and using the conventional definitions of dipole and quadrupole Diffractive reactions (where the color exchange is at amplitude level restricted to the color singlet) require to project the qqg system onto an overall color-singlet. This leads to the following replacement of Eq. (64) At the level of Eq. (66) this corresponds to replacing t a kl by t a ji on the right-hand side which translates into contracting the left-hand side with δ jk δ li . As a results one finds that all quadrupole operators in Eq. (68) factorize into a product of dipole operators. One finds for the leading N c terms: To obtain sub-leading terms in N c , it is again necessary to replace the operators N

Integrals
Integrations over plus-momenta are carried out trivially using the delta-functions associated with the vertices Eqs. (15) and (16). Integrations over minus-momenta can be performed using contour integrations. In the case of diagram one and two one finds with photon momentum fractions z 1 = p + /l + , z 2 = q + /l + and z 3 = k + /l + . In the case of diagrams three and four two (closely related) integrals are needed for each diagram. One has If one wishes to evaluate Eq. (20) with gluon correlators in the target given in transverse momentum space, the above form of the integrals is already sufficient. If gluon correlators in the target are provided in transverse coordinate space, the corresponding Fourier transforms can be be carried out using the integrals of appendix B of [19], which then allow to write our result in the large N c limit in the form already presented in [1]. For completeness we present the large N c result below 5 .

The (large N c ) result
With α em and α s the electromagnetic and strong coupling constants, and e f the electromagnetic charge of the quark with flavor f we obtain the following leading N c result 6 : where ψ i ≡ ψ i , i = 1, . . . , 3 and z 3 = 1 − z 1 − z 2 , while c L = 1, c T = 1/2. To obtain sub-leading terms in N c all operators N (4) , N (22) , N (24) , N (44) are to be replaced by 1/N c · N (4) (x 1 , x 1 , x 2 , x 2 ). With φ ij the azimuthal angle of x ij , i, j = 1 . . . , 3 and X 2 j = x 2 12 (z j + z 3 ) (1 − z j − z 3 ) , j = 1, 2, X 2 3 = z 1 z 2 x 2 12 + z 1 z 3 x 2 13 + z 2 z 3 x 2 23 , (74) we obtain amplitudes one can relate this process to Multi Parton Interactions (MPI) in a proton at large x in processes where one produces a (real or virtual) photon and a jet (in case of DPI) or a (real or virtual) photon in case of TPI in the forward rapidity kinematics. If one assumes that the target is accurately described by the CGC formalism, one can then extract valuable information on intrinsic parton correlations at large x in a proton [26].