#### Abstract

The foreseen capability to cover the far backward region at a Fixed-Target Experiment using the LHC beams allows one to explore the dynamics of the target fragmentation in hadronic collisions. In this report we briefly outline the required theoretical framework and discuss a number of studies of forward and backward particle production. By comparing this knowledge with the one accumulated in Deep Inelastic Scattering on target fragmentation, the basic concept of QCD factorisation could be investigated in detail.

#### 1. Introduction

In hadronic collisions a portion of the produced particle spectrum is characterised by hadrons carrying a sizeable fraction of the available centre-of-mass energy, the so-called leading particle effect. It is phenomenologically observed that for such hadrons their valence-parton composition is almost or totally conserved with respect to the one of initial-state hadrons [1]. In collisions, for example, protons, neutrons, and lambdas show a significant leading particle effect. For such processes, the production cross section peaks at very small transverse momenta with respect to the collision axis, a regime where perturbative techniques cannot be applied, giving insight into nonperturbative aspects of QCD dynamics in high energy collisions.

Quite interestingly, the leading particle effect has been observed in Semi-Inclusive Deep Inelastic Scattering (SIDIS). At variance with the hadronic processes mentioned above, such a process naturally involves a large momentum transfer. The presence of a hard scale enables the derivation of a dedicated factorisation theorem [2–5] which ensures that QCD factorisation holds for backward particle production in DIS. The relevant cross sections can then be factorised into perturbatively calculable short distance cross sections and new distributions, fracture functions, which simultaneously encode information both on the interacting parton and on the spectator fragmentation into the observed hadron. Despite of being nonperturbative in nature, their scale dependence can be calculated within perturbative QCD [6]. The factorisation theorem [2–5] guarantees that fracture functions are universal distributions, at least in the context of SIDIS.

Detailed experimental studies of hard diffraction at HERA have shown to support the hypothesis of QCD factorisation and evolution inherent the fracture function formalism. Furthermore, they led to a quite accurate knowledge of diffractive parton distributions [7–11], a special case of fracture functions in the very backward kinematic limit. For particles other than protons, proton-to-neutron fracture functions have been extracted from a pQCD analysis of forward neutron production in DIS in [12]. A set of proton-to-lambda fracture functions has been obtained by performing a combined pQCD fit to a variety of semi-inclusive DIS lambda production data in [13].

As theoretically anticipated in [2–4, 14, 15] and experimentally observed in hard diffraction in collisions at Tevatron [16, 17], QCD factorisation is violated for fracture functions in hadronic collisions. On general grounds, it might be expected, in fact, that the dynamics of target-remnants hadronisation is affected by the coloured environment resulting from the scattering in a rather different way with respect to the Deep Inelastic Scattering case.

Nonetheless, the tools mentioned above allow us to investigate quantitatively particle production mechanisms in the very backward and forward regions, to test the concept of factorisation at the heart of QCD, and to study the dependencies of factorisation breaking upon the species and the kinematics of the selected final state particle.

This physics program could be successfully carried on at a Fixed-Target Experiment using the LHC beams [18]. Novel experimental techniques are, in fact, available to extract beam-halo protons or heavy ions from LHC beams without affecting LHC performances. Such a resulting beam would be then impinged on a high-density and/or long-length fixed target, guaranteeing high luminosities. Furthermore, most importantly for the physics program to be discussed here, the entire backward hemisphere (in the centre-of-mass system of the collision) would be accessible with standard experimental techniques allowing high precision studies of the target fragmentation. Although measurements of particle production in the very forward region (close to the beam axis) might be challenging experimentally due to the high particle densities and large energy flow, the installation of dedicated detectors, like forward neutron calorimeters and/or proton taggers, could further broaden the physics program outlined above giving access to the beam fragmentation region.

The paper is organised as follows. In Section 2 we first give a brief theoretical introduction on the fracture functions formalism and to higher order corrections to the semi-inclusive Drell-Yan process. In Section 3 we outline different analyses which could be performed at AFTER@LHC with special focus on single hard diffraction. In Section 4 we summarise our results.

#### 2. Collinear Factorisation Formula

Fracture functions, originally introduced in DIS, do depend on a large momentum transfer. Therefore, in order to use them in hadronic collisions, a hard process must be selected. We consider here the semi-inclusive version of the Drell-Yan process:in which one hadron is measured in the final state together with a Drell-Yan pair. In such a process the high invariant mass of the lepton pair, , allows the applicability of perturbative QCD, while the detected hadron can be used, without any phase space restriction, as a local probe to investigate particle production mechanisms.

The associated production of a particle and a Drell-Yan pair in terms of partonic degrees of freedom starts at . One of the contributing diagrams is depicted in Figure 1.

Assuming that the hadronic cross sections admit a factorisation in terms of long distance nonperturbative distributions and short distance perturbative calculable matrix elements for the partonic process , predictions based on perturbative QCD are obtained convoluting the relevant partonic subprocess cross sections, , with parton distribution functions, and , and fragmentation functions, . The hadronic cross section, at centre of mass energy squared , can be symbolically written as [19, 20]where the convolutions are over the momentum fractions of the incoming and outgoing partons. The partonic indices , , and in the sum run on the available partonic subprocesses. The superscripts label the incoming hadrons and the presence of crossed terms is understood. This type of factorised hadronic cross section is expected to hold for hadrons produced at sufficiently high transverse momentum and it is widely and successfully used to compute cross sections for large momentum transfer processes in hadronic collisions. The Lorentz-invariant variable in (2) is defined byIn the hadronic centre-of-mass frame, where the second identity holds, is just the observed hadron energy, , scaled down by the beam energy . The variable , appearing in (2), is its partonic equivalent. Within this production mechanism, the observed hadron is generated by the fragmentation of the final state parton , and for this reason we address it as* central*. The amplitudes squared [21], however, are singular when the transverse momentum of the final state parton vanishes. In such configurations, the parent parton of the observed hadron is collinear either to the incoming parton or . As these phase space regions are approached, perturbation theory loses its predictivity. This class of collinear singularities escape the usual renormalisation procedure which amounts to reabsorb collinear divergences into a redefinition of bare parton distribution and fragmentation functions. Such singularities are likely to appear in every fixed order calculation in the same kinematical limits spoiling the convergence of the perturbative series. In [19, 20] a generalised procedure for the factorisation of such additional collinear singularities is proposed. The latter is the same as the one proposed in Deep Inelastic Scattering [22] where the same singularities pattern is also found, confirming the universality of collinear radiation between different hard processes. Such a generalised collinear factorisation makes use of fracture functions. These distributions obey DGLAP-type evolution equations which contain an additional inhomogeneous term resulting from the subtraction of collinear singularities in the target fragmentation region [6, 22]. Such equations allow us to resume the corresponding large logarithmic corrections to all orders in perturbation theory. Bare fracture functions, , describe the hadronization of the spectators system in hadron-induced reactions. They express the conditional probability to find a parton initiating the hard scattering while an hadron is produced with fractional momentum in the target fragmentation region of the incoming hadron .

The use of fracture functions allows for particles production already to , since the hadron can be nonpertubatively produced by a fracture function itself. Therefore, the lowest order parton model formula can be symbolically written asand it is sketched in Figure 2. The superscripts in (4) indicate from which incoming hadron, or , the outgoing hadron is produced through a fracture function. In order to complete the calculation to accuracy we should consider higher order corrections to (4). Since in this case the hadron is already produced by fracture functions, final state parton radiation should be integrated and the resulting contribution added to virtual corrections. One of the contributing diagrams is depicted in Figure 3. The general structure of these terms isWe refer to them as to the* target* fragmentation contributions. Their calculation is, apart from minor differences in kinematics, completely analogous to the one of the inclusive Drell-Yan case. The factorisation procedure, first elaborated in [22] in the context of SIDIS, amounts to substitute in (4) the bare fracture and parton distributions functions with their renormalised version [19, 20]. Renormalised parton distributions and fracture functions homogeneous terms do cancel, as in the inclusive Drell-Yan case, all singularities present in (5). The additional singularities in (2) are cancelled by the combination of parton distributions and fracture functions inhomogeneous renormalisation terms. Adding all the various contributions, the resulting hadron- integrated cross section, up to order , is then infrared finite [19, 20] and can be symbolically written aswhere and is the number of colors. We refer to the previous equation as to the collinear factorisation formula for the process under study. The next-to-leading order coefficients and have been calculated in [20], making the whole calculation ready for numerical implementation.

We stress, however, that our ability to consistently subtract collinear singularities in such a semi-inclusive process is a necessary but not sufficient condition for factorisation to hold in hadronic collisions. The one-loop calculation outlined above in fact does involve only the so-called active partons. It completely ignores multiple soft parton exchanges between active and spectators partons, whose effects should be accounted for in any proof of QCD factorisation. Therefore, there is no guarantee that fracture functions extracted from SIDIS can be successfully used to describe forward or backward particle production in hadronic collisions. Reversing the argument, such a comparison may instead offer new insights into nonperturbative aspects of QCD and to the breaking of factorisation.

#### 3. Single Hard Diffraction at AFTER@LHC

As an application of the formalism presented in the previous sections we will consider single hard diffractive production of a Drell-Yan pair: where we have indicated in parentheses the four momenta of the relevant particles. We present in the following cross sections differential in the virtual photon variables. The subsequent decay of the virtual photon into a lepton pair can be easily included so that realistic cuts on leptons rapidity and transverse momentum can then be applied. We consider the AFTER@LHC kinematic setting in which a 7 TeV proton beam collides on a fixed target proton leading to a centre-of-mass energy of GeV. We consider the projectile proton moving in the positive direction and at rest in the laboratory. The diffractively produced proton has in general almost the incoming projectile proton energy and very small transverse momentum as measured with respect to the collision axis. The detection of such fast protons will in general require the installation of forward proton taggers. The lepton pair instead will be measured by the main AFTER@LHC detector. This kinematical configuration is pictorially represented in the right plot of Figure 2.

Diffractive processes have been intensively analysed in DIS at HERA collider, revealing their leading twist nature. From scaling violations of the diffractive structure functions [7, 10, 11] and dijet production in the final state [8, 9] quite precise diffractive parton distributions functions (dPDFs) have been extracted from HERA data, which parametrise the parton content of the color singlet exchanged in the -channel. The comparison of QCD predictions for single diffractive hard processes based on diffractive parton distributions measured at HERA against data measured at Tevatron [16, 17] ( TeV), adopting a factorised ansatz as in (4), has indeed revealed that these processes are, not unexpectedly [14, 15], significantly suppressed in hadronic collisions. This conclusion persists even after the inclusion of higher order QCD corrections [23]. Complementing these results with the forthcoming ones from LHC at higher centre-of-mass energies ( TeV) and the ones from AFTER@LHC at GeV will give insight into the energy dependence of the so-called rapidity gap survival (RGS) probability in a wide range in . Since the theoretical computation of the RGS factor is highly model dependent, we decided not to include it in our predictions. Therefore, the latter must be considered as upper bounds.

Diffractive parton distributions are in general proton-to-proton fracture functions . They depend upon the final state proton fractional energy loss, with given in (3), the fractional momentum of the interacting parton with respect to the pomeron momentum, and the virtuality . In general fracture functions may depend also upon the invariant momentum transfer at the proton vertex [24]. In all diffractive structure functions measurements at HERA, out of which dPDFs are determined, is integrated up to some . In this case dPDFs obey ordinary DGLAP evolution equations [25] as their extended, -dependent version [24]. In the present paper we use dPDFs form [7] which are defined by GeV^{2}. Since they are extracted from large rapidity gap data where the proton is not directly measured, they contain a contribution (23%) from the so-called proton dissociation contribution. In order to use dPDFs in the present context we first note thatThe extra factor comes from the Jacobian of the change . We then rearrange (4) in terms of new variables obtainingwith . For simplicity we consider here leading order formulas but the extension to higher order is straightforward. In (9) we use parton distribution functions from [26]. We show explicitly the dependence of fracture and parton distributions functions upon the factorisation scale, . Predictions are obtained with this scale set to . Theoretical errors associated with higher order corrections are instead estimated varying such scale in the range and .

In Figure 4 we present predictions for the distribution. In left plot we consider a Drell-Yan pair of mass GeV^{2}. The distribution shrinks as lower -values are approached whereas, from hard diffraction at HERA, it is well known that diffractive cross sections rise as an inverse power of . Such an effect therefore is then attributed to phase space threshold effects. The Drell-Yan invariant mass constraint can be rewritten in the diffractive case as , which can be cast (for and ) in upper bound on the invariant mass at fixed and . This hypothesis is further supported in the right plot of Figure 4, where differential distributions are presented for three values of . The lowest values of are then accessed only by lowering the invariant mass of the pair. We note that, even considering the maximum value of , single diffractive production of and is beyond the kinematic reach at AFTER@LHC. In the left panel of Figure 5 we present the prediction for the distribution at a fixed value of . The cross section, as expected, is fast falling as an inverse power of . The distribution is particularly instructive since it allows us to study the possible dependence of the RGS factor on and therefore to get an insight into the underlying dynamics. In the right panel of Figure 5 we present the -differential cross section multiplied by a factor , which compensates the fast fall off of the electromagnetic cross section so that all the dependence is accounted for by that of fracture and parton distributions. These curves and the corresponding slopes, however, can not be readily interpreted as genuine results of QCD evolution of fracture and parton distributions functions because of the threshold effect mentioned above appearing at such moderate values of .

**(a)**

**(b)**

**(a)**

**(b)**

By changing variable from to the virtual photon centre-of-mass rapidity, ,equation (9) can be further manipulated to give the three-differential cross section:The rapidity range for diffractive Drell-Yan production readswhich, as expected, turns out to be asymmetric given the kinematic constraint . The rapidity range for the inclusive Drell-Yan case is recovered simply setting . The rapidity distribution is particularly sensitive to the shape of the diffractive parton distributions. This distribution will be useful to investigate any possible kinematic dependence of the RGS factor. In the left panel of Figure 6 we present the centre-of-mass rapidity distribution at fixed GeV^{2} and . In this frame the distribution is shifted at negative values of . Therefore, on average, the parton originating from the target proton carries more momentum than the one originating from the pomeron. Since the rapidity is additive under boost along the collision axis we may easily boost the to the laboratory frame by usingwith the proton mass. In the AFTER@LHC kinematics this implies a rapidity shift . The rapidity distributions in the laboratory frame for a Drell-Yan pair of mass GeV^{2} and for three different values are presented in the right panel of Figure 6. One may notice from the plot that for increasing the Drell-Yan pair spans a wider rapidity range and the corresponding spectrum is increasingly more forward. It might be useful to discuss the single diffractive Drell-Yan pair production in conjunction with the analogous inclusive process. Such a comparison is presented for the centre-of-mass rapidity distributions in the left plot of Figure 7 for a common Drell-Yan pair of mass squared GeV^{2} and, for the diffractive case, integrated in the range . The rapidity distribution in the single diffractive case is strongly asymmetric whereas in the inclusive case it is symmetric around . This effect is primarily due to the different kinematics of the two processes and to the different fractional momentum distributions between parton and fracture distributions. In the right plot of Figure 7 the ratio between the two previous distributions is presented. The ratio gives direct information on the suppression factor between the single diffractive to the inclusive process, assuming a factorised expression for the former (11). Such a ratio might be convenient from the experimental side since many lepton detection systematics will cancel. On the theoretical one it is expected to be more stable against the inclusion of higher order corrections. In the present case, for example, the factorisation scale is simultaneously varied both on the numerator and denominator resulting in a reduced theoretical error band with respect to the one obtained for absolute cross sections.

**(a)**

**(b)**

**(a)**

**(b)**

We wish to end this section with a brief overview of other possible applications of the proposed formalism. A completely analogous program can be performed for the associated production of forward neutron and a Drell-Yan pair, . The production of forward neutron in DIS at HERA has shown a leading twist nature. From scaling violations of the semi-inclusive neutron structure functions a set of proton-to-neutron fracture functions has been extracted from data in [12] which can be used to predict forward neutron rate in hadronic collisions. As in the case of hard diffraction, both physics programs would highly benefit from the installation of a dedicated instrumentation for the measurements of fast neutrons and protons quite close to the beam axis. Measurements in the forward region, although problematic experimentally, give in fact direct access to the study of the beam fragmentation region.

As a third application we consider hyperon production associated with a Drell-Yan pair, , where generically indicates either a or hyperon. At very low transverse momentum, the longitudinal momentum spectrum should show a significant leading particle effect, which can be predicted, assuming factorisation as in , by the proton-to-lambda fracture functions set obtained from a fit to SIDIS data in [13]. On the other hand, the spectrum in the same kinematical conditions should instead show almost no leading particle effect, giving access to the proton-to- fracture functions. We note, in general, that the particle-to-antiparticle fracture function is indeed an interesting and almost unknown distribution. On the other hand, if one considers or at sufficiently large transverse momentum, their combined analysis, described by , should allow for the investigation of parton hadronisation into hyperons in the QCD vacuum as parametrised by fragmentation functions.

As a last application we consider the associated production of one particle and a Drell-Yan pair in the context of multiparton interactions. The latter process has already been used to investigate the contamination of the so-called underlying event [27] to jet observable and has been successfully used to study underlying event properties [28]. If the detected hadron is measured at sufficiently large transverse momentum, the latter constitutes a natural infrared regulator for the partonic matrix elements. In this kinematical conditions we also expect a rather small contribution from fracture functions. Therefore, the central term, , can be used to estimate the single parton scattering contribution to the process. The latter might be considered as the baseline to study the double (or multiple) parton scattering contributions to the same final state, where, for example, the primary scatterer produces a Drell-Yan pair while the secondary one produces the detected hadron .

#### 4. Conclusions

We have briefly reviewed a perturbative approach to single particle production associated with a Drell-Yan pair in hadronic collisions. On the theoretical side we have shown that the introduction of fracture functions allows a consistent factorisation of new class of collinear singularities arising in this type of processes. The factorisation procedure coincides with the one used in DIS confirming, as expected, the universal structure of collinear singularities and supports the proposed collinear factorisation formula. On the phenomenological side we have outlined some areas in which the formalism can be fully tested. In particular, focusing on the AFTER@LHC kinematical range, we have presented numerical predictions for the single diffractive production of virtual photons. The study of such a process might improve our understanding of nonperturbative aspects of QCD and it allows one to explore in detail the nature of factorisation breaking at intermediate energies.

#### Conflict of Interests

The author declares that there is no conflict of interests regarding the publication of this paper.