Single-diffractive Drell–Yan pair production at the LHC

We present predictions for single-diffractive low-mass Drell–Yan pair production in pp collisions at the LHC at s=13\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\sqrt{s}=13$$\end{document} TeV. Predictions are obtained adopting a factorised form for the relevant cross sections and are based on a new set of diffractive parton distributions resulting from the QCD analysis of combined HERA leading proton data. We discuss a number of observables useful to characterise the expected factorisation breaking effects.


Introduction
The diffractive physics program pursued at the HERA ep collider in the recent past has substantially improved our knowledge on the dynamics of this class of processes. In the deep inelastic regime, the presence of a hard scale enables the derivation of a dedicated factorisation theorem [1][2][3][4] which allows the investigation of the partonic structure of the colour singlet exchanged in the t-channel. From scaling violations of the diffractive Deep Inelastic Scattering (DDIS) structure functions, quite precise diffractive parton distributions functions (dPDFs) have been extracted by performing QCD analysis [5][6][7][8] of available data.
With this tool available, factorisation tests have been conducted in order to investigate the range of validity of this hypothesis in processes other than DDIS. Factorisation has been shown to hold, as expected theoretically, in diffractive dijets production in DIS, where NLO predictions based on dPDFs well describe experimental cross sections [6,9] both in shape and normalisation. Factorisation breaking effects are expected to appear in diffractive photoproduction of dijets due to the resolved component of the quasi-real photon. In such a case, however, H1 [9] reported a global suppression factor of data over NLO theory around 0.5 while ZEUS [10] found the same ratio compatible with unity. To date, these conflicting results prevent to draw a conclusive statement a e-mail: federico.alberto.ceccopieri@cern.ch about factorisation in this case. We note, however, that the measurements of diffractive dijet photoproduction in ultraperipheral collisions in pp and p A collisions at the LHC [11] may offer an alternative way to settle this issue.
Complementary informations on the nature of diffraction has been provided by hard diffraction measurements in hadronic collisions. As theoretically anticipated in Refs. [1][2][3]12,13] and experimentally observed in pp collisions at Tevatron [14][15][16], factorisation is strongly violated in such a case. In particular, predictions based on a factorised expressions for the relevant cross sections in terms of diffractive parton distributions extracted from HERA data overestimate hard diffraction measurements by a factor O(10) [17]. This conclusion persists even after the inclusion of higher order QCD corrections [18].
A rich program at the LHC is being pursued in diffractive physics by all Collaborations either based on the identification of large rapidity gaps (LRG) [19][20][21][22][23][24] or by using dedicated proton spectrometers [25,26]. Complementing Tevatron ( √ s = 1.96 TeV) results with forthcoming ones from the LHC at higher centre-of-mass energies ( √ s = 8, 13 TeV) will give information on the energy dependence, if any, of the suppression factor, the so called rapidity gap survival (RGS) probability. Hopefully, they will allow to study its kinematic dependences, among which the one on the scale characterising the hard process appears to be particularly relevant. In the simplest scenario, it will be possible to clarify whether factorisation may still hold but revisited in a weak form through a global or local rescaling of diffractive PDFs extracted from DDIS and to study their degree of universality among different hard processes in hadronic collisions.
The purpose of the present paper is to present predictions for the single-diffractive Drell-Yan pair production at the LHC at √ s = 13 TeV, one of the clean and simple measurable process in hadronic collisions. In such a process, the invariant mass of the lepton pair can easily be reconstructed and, depending on experimental capabilities, pushed to rather low values, allowing a detailed characterisation of the hard-scale dependence of the suppression factor. Although estimates of the latter are present in the literature for the specific process at hand [27][28][29], we take a conservative approach and avoid to introduce any suppression factor. We further assume factorisation to hold and adopt factorised expressions for the relevant cross sections. A preliminary set of, newly generated, diffractive parton distributions extracted from combined leading proton HERA data will be used for the calculation. In view of the expected factorisation breaking effects in hard, single-diffractive, measurements in hadronic collisions, the obtained values for the cross sections should be considered as upper bounds.
Given the exploratory nature of the analysis, more intended as a feasibility study, theoretical predictions are calculated to leading order accuracy. We take into account, however, the virtual photon decay into leptons so that cross sections can be studied as a function of, measurable, final state leptons kinematics. This allows us to explore the phase space available for the process and to estimate the impact of typical experimental cuts on the transverse momenta and rapidities of the leptons.
From QCD analyses performed in DDIS and anticipating the results of the next section, we know that the colour singlet exchanged in the t-channel is a gluon-enriched state. Since gluonic contributions to Drell-Yan production starts to O(α s ) in perturbation theory, an accurate estimation of the suppression factor will require the inclusion of higher order corrections. The impact of the latter and a detailed report on the extraction of diffractive parton distributions to NLO accuracy will be presented in a companion publication.
The paper is organised as follows. In Sect. 2 we report in some details the extraction of diffractive PDFs from combined HERA leading proton data. In Sect. 3, making use of such distributions, we present results for single-diffractive Drell-Yan production in pp collisions at the LHC at √ s = 13 TeV. In Sect. 4 we summarise our results.

FIT overview
Diffractive DIS belongs to the Semi-Inclusive lepton-proton DIS process of the type where, along with the outgoing lepton, an additional proton p is detected in the final state. In Eq. (1) X stands for the unobserved part of the hadronic final state and we indicate particles four-momenta in parentheses. In the lp centre-of-mass system, diffractive DIS events are then characterised by outgoing protons with a large momentum fraction of the incident proton and quite small values of the transverse momentum measured with respect to the collision axis, i.e. in the target fragmentation region of the incident proton. The kinematic variables used to describe the DIS process are the conventional Lorentz invariants with q = k − k . Final state protons are instead described by the fractional momentum of the singlet exchange with respect to the proton momentum, x IP , and the invariant momentum transfer t at the proton vertex: where typical DDIS selection requires x IP 0.1 and |t| 1 GeV 2 . In the following we will use the scaled fractional momentum variable β defined by β = x B /x IP . This is interpreted as the fractional momentum of interacting parton with respect to pomeron fractional momentum x IP . The data are often presented in terms of the reduced lp cross section, σ D(4) r , which depends on the diffractive transverse and longitudinal structure functions F D(4) 2 and F D(4) L , respectively. In the one-photon exchange approximation, it reads According to the factorisation theorem [1][2][3][4], structure functions appearing in Eq. (4), are factorised into perturbatively calculable short-distance cross sections and diffractive parton distributions The index i runs on the flavour of the interacting parton. The hard-scattering coefficients C ki (k = 2, L) are perturbatively calculable as a power expansion in the strong coupling α s and depend upon μ 2 F and μ 2 R , the factorisation and renormalisation scales, respectively. The C ki coefficient functions are the same as in fully inclusive DIS. Diffractive PDFs f D i/ p (β, μ 2 F , x IP , t) appearing in Eq. (5) are protonto-proton fracture functions [30] in the very forward kinematical limit and can be interpreted as the number density of interacting partons at a scale μ 2 F and fractional momentum β conditional to the detection of a final state proton with fractional momentum 1 − x IP and invariant momentum transfer t. The t-unintegrated diffractive PDFs appearing in Eq. (5) obey standard DGLAP [31][32][33] evolution equations [34]. The same statement holds when they are integrated over t in a limited range [35]: In this paper we analyse the combined H1 and ZEUS diffractive DIS cross sections measurements [36] of the process in Eq. (1) where leading protons are measured by dedicated forward spectrometers. The centre-of-mass energy for the e + p scattering is √ s = 318 GeV. This data set covers the phase space region 2.5 < Q 2 < 200 GeV 2 and 0.0018 < β < 0.816 and it has widest coverage in the proton fractional energy loss, 0.00035 < x IP < 0.09, subdivided in 10 bins in x IP , with an average of 20 points per-x IP bin for a total of 192 points. At variance with all other DDIS cross sections measurements, the squared four-momentum transfer at the proton vertex, t, is integrated in the restricted range 0.09 < |t| < 0.55 GeV 2 in order to minimise systematic uncertainties originating from t-extrapolation of the various measurements outside their respective measured ranges. The reduced cross sections in Eq. (4) are integrated over t in such a range and the diffractive PDFs in Eq. (6) are defined accordingly. For x IP < 0.03 the data set overlaps with high-statistics LRG data set and for 0.03 < x IP < 0.09 it provides the best experimental information available on diffractive DIS cross sections. The combination procedure, in general, allows a reduction of the systematic uncertainties via cross-calibration of the various measurements. The direct detection of the forward proton allows one to avoid any systematics associated with the large rapidity gap selection. By definition, these data are free from the proton dissociative background which has been found to contribute around 23% of the diffractive DIS cross sections based on LRG selection [5]. Therefore this set of data provides the most precise knowledge about the absolute normalisation of diffractive DIS cross sections. These advantages, however, come at the price of increased uncertainties relative to LRG data given the reduced statistics of the sample. Diffractive parton distributions extracted form this data set will be used in the context of single hard diffraction in hadronic collisions in conjunction with ordinary parton distributions. In order to avoid any mismatch between inclusive and diffractive PDFs we adopt leading order CTEQL1 parton distribution set [37] evolved in the zero-mass variable-flavour-number scheme (ZM-VFNS). The evolution of diffractive PDFs is performed within the same scheme and to the same accuracy by using QCDNUM17 [38] program. The QCD parameters are the ones quoted in Ref. [37]. In particular we set the charm and bottom masses to m c = 1.3 GeV and m b = 4.5 GeV, respectively, and the strong coupling is evaluated at one loop setting In general factorisation theorem [1][2][3][4] for diffractive DIS in the form of Eq. (5) holds at fixed values of x IP and t so that the parton content of the coloursinglet exchange described by f D i is uniquely controlled by the kinematics of the outgoing proton. Therefore, at least in principle, dPDFs may differ at different values of x IP and t. This idea has been successfully tested [7] in the analysis of LRG data from Ref. [5]. In the present context, given the limited number and accuracy of the data points in each x IP bin, we use a simpler approach, namely a fully factorised β − x IP ansatz for the flavour-symmetric singlet and gluon diffractive parton (momentum) distributions defined at the initial scale Q 2 0 : The initial conditions in Eq. (7) are characterised by a common flux factor F(x IP ) controlled by a single power at low x IP . An extra modulation, controlled by parameters f 1 and f 2 , is introduced to accommodate the data at larger values of x IP . In order to guarantee the vanishing of the singlet distribution on the endpoint, we fix the large-β behaviour of the singlet by setting C q = 0.5 but additional freedom at intermediate values of β is allowed leaving D q and E q parameters free in the minimisation. Since the gluon distributions is only indirectly fixed by the slope of the reduced cross section, the gluon parameters B g and C g are highly correlated and we decide to fix C g = 0.5 for a total of 9 free parameters. Such distributions, once evolved, are used to calculate the diffractive structure functions F D 2,L with the help of the QCDNUM17 convolution engine and to reconstruct the diffractive reduced cross sections in Eq. (4) which are then minimised against data [36] with the help of the MINUIT [39] program. The choice of Q 2 0 is optimised performing a scan giving the best χ 2 value for Q 2 0 = 1.5 GeV 2 . By using total errors quoted in the experimental analysis and the standard χ 2 definition, we obtain a total χ 2 /d.o. f. = 0.913. We report in Table 1 the best parameters and the breakdown of the contributions to χ 2 function in each x IP bin. According to these numbers there is no misrepresentation of the data across the x IP bins. The comparison of the best-fit results and the reduced cross sections is presented in Fig. 1 for four representative values of x IP as a function of Q 2 or β. We supplement the best-fit parametrisations with an additional set of parametrisations obtained according to the Hessian method outlined in Refs. [40,41] which allows the experimental uncertainties to propagate to arbitrary observables. We note that our initial conditions assume a common β-shape for the diffractive PDF in all x IP bins.
This theoretical hypothesis, in turn, determines an unrealistic precise determination of the diffractive PDFs if associated with the standard χ 2 = 1 criterion, often exceeding the precision of the data. In order to correct for such an effect and to obtain a more conservative error estimate we choose a tolerance criterion χ 2 = 10 (one unit for x IP bin) and dPDFs alternative parametrisations are obtained with this choice. We have checked by explicit evaluation that each parametrisation gives a consistent value for the χ 2 function, χ 2 best + χ 2 . The error bands presented in the plots are obtained according to this criterion. In Fig. 2 we present the singlet and gluon momentum distributions in two x IP bins at different scales. The singlet shows a bump in the large β region (β 0.5) at the lower scale which is progressively washed away by evolution at higher scales. The rise of the gluon distribution at small β is accelerated by pQCD evolution and it is the dominant contribution for β 0.1.
We have further performed two consistency checks detailed below. The first one concerns the diffractive longitudinal structure function which contributes starting from O(α s ) and it is absent to the accuracy of the present calculation. Since its dominant contributions appear in the large-y region, the fit  x IP = 0.035  has been repeated with the cut y < 0.5 imposed. The second one addresses the issue, reported in previous analyses [5][6][7][8], of the inclusion in the fit of the lowest Q 2 points. For such a reason, the minimisation has been repeated by including only data points for which Q 2 > 6 GeV 2 . In both cases we observe a modest decrease in the χ 2 /d.o. f . However, as shown in Fig. 3, the resulting parametrisations are compatible, within uncertainties, with the ones obtained without imposing the cuts. Given the substantial stability of the results against variation of the phase space boundary of data included in the fit, we consider the "no cut" scenario as our default choice and use the corresponding best-fit parametrisations in the next section. We conclude this section presenting a comparison between dPDFs obtained in this work with the ones already present in the literature. In Fig. 4 we compare our dPDFs with LO Fit B (without error estimation) and NLO Fit B (with error estimation) parametrisations presented in Ref. [5] and obtained from DDIS data selected with the large rapidity gap method. Since the analysis presented in Ref. [5] excludes DDIS data for which Q 2 < 8.5 GeV 2 , the comparison is performed at the higher scale Q 2 = 25 GeV 2 , a value for which no extrapolation is required. The normalisation of dPDFs extracted in Ref. [5] is larger than the one of the present analysis given the larger t-range, |t| < 1 GeV 2 , used in Ref. [5] to define the DDIS cross sections. A further upward normalisation shift is generated by the proton dissociative contribution which amounts to 23% of the DDIS cross section measured with the large rapidity gap method [5] and it is absent in the present analysis. In the left panel of Fig. 4, the comparison is performed at x IP = 0.0025 where both data sets overlap. In this case, the larger data sample in this x IP -bin used in the fit of Ref. [5] induce smaller experimental uncertainties (yellow band) on dPDFs with respect to the ones presented in this work. The virtue of the present parametrisations can be better appreciated for x IP > 0.03, being the latter the maximal x IP -value included in the analysis of Ref. [5]. In the right panel of Fig. 4, the comparison is therefore performed at x IP = 0.05. Diffractive PDFs from Ref. [5] are in a full extrapolation range and we observe a substantial increase of their uncertainties (green error band) driven by the theoretical errors associated with the parametrisations of the flux factor at large x IP . On the contrary, our parametrisations, despite the lower statistical power of the data set used in this analysis, show a much reduced error. To conclude, we note that the singlet distribution from Ref. [5] Q 2 > 6 GeV 2 x I P = 0.075 Best-fit distributions with uncertainties (the band corresponds to χ 2 = 10) are compared with parametrisations returned by the fit with the cut y < 0.5 (solid) or Q 2 > 6 GeV 2 (dashed) imposed is stable in shape and normalisation between the LO and NLO versions and shows the presence of a large-β bump in nearly the same position with respect to the one obtained in the present analysis. On the contrary we observe a steeper behaviour of the LO FitB gluon with respect to the NLO FitB one and to the one from our fits. This is probably due to the fixed flavour number scheme used in Ref. [5] which, by construction, induces a larger gluon distribution with respect to the one obtained with a variable-flavour-number scheme.

Single-diffractive Drell-Yan production
The signature of hard diffraction in hadronic collisions is the presence of hard-scattering process associated with the pro-duction of a leading proton. Among many others, we consider here the simplest hard-scattering process, namely the Drell-Yan pair production. Therefore we consider the reaction The invariant mass of the pair q 2 = ( p 3 + p 4 ) 2 = Q 2 is chosen to be large enough so that perturbative QCD can be applied. In hadronic collisions, the Lorentz-invariant variable z is used to characterise final state hadrons and is defined by In the hadronic centre-of-mass frame, where the second identity holds, z is just the observed proton energy, E * p , scaled down by the beam energy, In Eq. (8), we assume that the proton with momentum P 1 is moving in the +z direction and the leading proton with momentum P is produced quasi-collinearly to P 1 at large and positive rapidities. At the cross section level, diffractive parton distributions for the proton with momentum P 1 will be used. The same process, of course, may occur also in the opposite emisphere and, since the hadronic initial state is symmetric, will be not considered here.
At the partonic level and to lowest order in the strong coupling, the process proceeds via the annihilation of a quarkantiquark pair into a massive virtual photon which subsequently decays into a opposite-sign lepton pair. To be definite we consider here the decay into muons: Before discussing our results we found useful to sketch some details of the calculation. The incoming parton momenta in the hadronic centre-of-mass-system are given by with s = (P 1 + P 2 ) 2 . We choose as final state variables the lepton rapidities y 3 , y 4 , and lepton transverse momentum, p t . In terms of the latter, the four-momenta of the leptons are given by with p t = | p t | and q = p 3 + p 4 . Assuming factorisation to hold, the differential cross section, to leading order accuracy, involves appropriate products of diffractive and ordinary parton distributions functions. It reads where the sum runs over quark and antiquarks. In actual calculations we have set the factorisation scale to μ F = M μμ . We further introduce the leptons rapidity sum, Y , and differenceȳ: The former defines the rapidity of the virtual photon. The partonic Mandelstam invariants appearing in Eq. (15) are then given by In terms of these variables the momentum fractions are given by Since the two momentum fractions cannot exceed unity, the following bounds can be derived: with τ = M μμ / √ s. Given the kinematic constraint x 1 x IP , the pair rapidity spans an increasingly asymmetric range as x IP decreases. For x IP < √ τ , the pair is entirely in the Y < 0 rapidity range. Formally, the rapidity range for the inclusive Drell-Yan case is recovered simply setting x IP = 1 in Eq. (22).
In the present analysis we focus on diffractive processes tagged with dedicated instrumentation [25,26]. We choose dσ/dp μ t [pb/GeV] the proton fractional momentum loss to be in the range 10 −4 < x IP < 10 −1 , with maximal overlap with the range measured at HERA [36]. Predictions presented in the following are integrated over the t-range of the data [36] out of which dPDFs are estracted, i.e. 0.09 < |t| < 0.55 GeV 2 . We set the centre-of-mass energy of the pp collisions to √ s = 13 TeV. The invariant mass of the muon pair is required to be in the range 2 < M μμ < 20 GeV, a range of virtualities in line with those measured at HERA. We assume that the J/ and ϒ contributions, which both lie within this mass range, can be properly subtracted from the data sample. We require both muons to have rapidity |y μ | < 2.45 but we do not apply cuts either on the muons transverse or three momenta.
The resulting fiducial cross sections for single-diffractive Drell-Yan pair production is reported in Table 2. In the case that proton spectrometers are installed on both sides of the interaction point, such a result for the fiducial cross section should be doubled. As already anticipated, the quoted result does not include any rapidity gap suppression factor and predictions refer to virtual photon decay in the muon channel. The first error represents the propagation of experimental uncertainties as obtained in the diffractive PDF fit. The second one, of theoretical nature, is obtained varying the factorisation scale μ 2 F appearing in both distributions in Eq. (15) in the range 1/2M 2 μμ < μ 2 F < 2M 2 μμ . In this regime of relatively low Q 2 values where diffractive and inclusive parton distributions evolve faster, we find that the dominant error source, of theoretical nature, is associated with missing higher order corrections. The latter are known to high accuracy for a number of distributions and will be considered in a separate publication. In the present contest, predictions can be stabilised against factorisation scale variation by considering appropriate ratios of diffractive over inclusive cross sections. This issue will be discussed in some details at the end of this section. We begin our overview of predictions showing in Fig. 5 single-differential cross sections in the fiducial phase space region specified in Table 2. The pair invariantmass distribution is shown in the top left panel and rapidly falls as an inverse power ofŝ = M 2 μμ typical of annihilation processes into massive states, as it can be read out from Eq. (15). In the top right panel the muon transverse momentum distribution is presented. Its kinematically allowed range extends up to p t = M max μμ /2. Given the fast falling nature of the M μμ -distribution, dominated by low values of the invariant, the muon transverse momentum distribution shows a maximum (the Jacobian peak) at p t = M min μμ /2. The muonpair-rapidity distribution, presented in the bottom left panel, is slightly asymmetric and indicates a preference for the virtual photon to populate the negative rapidity hemisphere (the one containing the dissociated proton, in the chosen reference frame). We note that, despite phase space limitations introduced by Eq. (22) and the difference between diffractive and ordinary parton distributions, the muon pair populates the available rapidity range, as defined by the muon rapidity cuts and by Eq. (16). In the bottom right panel we present the x IP distribution. In general, it is well known that such distribution behaves approximately as an inverse power of x IP at small x IP . In the present case, the flattening of the distribution at small x IP can be ascribed to the shrinkage of phase space for the production of massive pair, since the maximum partonic centre-of-mass energy is reduced to √ x IP s. The kinematic of the scattered proton induces peculiar features on muon-pair production whose effects are presented in Fig. 6 in terms of the muon-pair rapidity, Y, in various bins of x IP . The distri-butions is strongly asymmetric at the lowest values of x IP , where the muon pair populates the negative rapidity range (dissociated proton direction) due to the kinematic constraint x 1 < x IP . In the intermediate x IP range the pair starts to populate the positive emisphere (diffractive proton direction) with a tendency to show a maximum in this range. At even higher values of x IP , the available centre-of-mass for the reaction increases and the distribution progressively turns into a symmetric one. Given the relatively light masses produced, this regime is sensitive to parton distributions evaluated at relatively small values of β and x 2 , the symmetry of the rapidity distribution indicates that the shapes of the sea component of the diffractive and ordinary distributions are similar, both being driven by QCD evolution. This complicated pattern is further illustrated in Fig. 7 where the single-differential cross section as a function of Y is shown in four different ranges of the pair invariant mass and integrated over x IP . In all mass bins, the distributions show a maximum in the negative rapidity range, a signal that the interacting parton from the dissociated proton carries, on average, slightly more momentum with respect to the one originating from the scattered proton. In Fig. 8 we present single-differential distributions as a function of x IP in four different invariant-mass ranges. As the invariant mass increases, we observe a progressive flattening  of the distributions at small x IP . This effect is due to the phase space reduction induced by the constraint M 2 μμ = βx IP x 2 s, which at low x IP disfavours the production of increasingly massive pair. In Fig. 9 we present single-differential cross section as a function of β, the fractional momentum of the interacting parton with respect to the one of the colour singlet exchanged in the t-channel, integrated in various bins of M μμ and x IP . Such distributions offer insight in the sensitivity of the cross section to diffractive parton distributions, modulo kinematics effects. In the lowest x IP bin the distributions span all the allowed β range and progressively shrinks at large β as x IP increases, a natural consequence of momentum conservation. As already shown in Figs. 6 and 9, the distributions in the pair rapidity Y are asymmetric around Y = 0. The asymmetry decreases both as the mass of the pair increases and as x IP increases. Such an effect is absent in the inclusive Drell-Yan case initiated by a symmetric initial state. This effect is better appreciated considering the absolute and relative asymmetries, A a and A r , respectively, defined by They differ by a different normalisation and are shown in the left panel of Fig. 10 along with the associated theoretical errors estimated via factorisation scale variations. The absolute asymmetry, integrated over all masses and proton energy loss, reaches its maximum 3% at Y 1 and is always negative, implying that the muon pair is produced mostly in the hemisphere opposite to the one containing the scattered proton. The relative asymmetry shows an even more pronounced behaviour. Since for Y → 0 the rapidity distribution becomes symmetric and the cross section is maximal, see the bottom left panel of Fig. 5, the asymmetry vanishes and the corresponding error band shrinks in that limit. Both asymmetries, being normalised to single-diffractive cross sections, are not affected by uncertainties due to the rapidity gap survival factor. They are primarily sensitive to the shape of diffractive PDFs and show reasonably stability against scale variations. Depending on the accumulated integrated luminosity, this predicted behaviour, absent in the inclusive case, could be exploited to correlate the forward proton detection with the central Drell-Yan production. As discussed at the beginning of the section, predictions are affected by large theoretical errors associated with scale variations. Such uncertainties can be conveniently reduced by considering the ratio R of diffractive to inclusive cross sections which also offers the advantage reducing the experimental systematics errors. In Eq. (24) Y stands for the selected hardscattering process (Drell-Yan in this case) and X for the unobserved part of the final state. At Tevatron the ratio R has been measured in a variety of final state [14][15][16] and it shows a quite stable behaviour with a value close to 1%. For the single-diffractive Drell-Yan production in pp collisions at √ s = 13 TeV, the ratio R is presented in the right panel of Fig. 10. Given our leading order estimate of the inclusive Drell-Yan cross section, R varies between 6 and 8% and decreases mildly as a function of the invariant mass of the pair, M μμ . This prediction, however, does not take into account the RGS suppression factor. In this respect it would be interesting to check whether the data follow at least the shape of the ratio as a function of M μμ .

Conclusions
In this paper we have considered the single-diffractive production of low-mass Drell-Yan pair in pp collisions at the LHC at √ s = 13 TeV. Predictions are based on a fully factorised approach for the cross section which makes use of a set of diffractive parton distributions obtained from a QCD fit to combined leading proton DIS data from HERA. A number of distributions are presented both in terms of Drell-Yan pair and scattered proton variables. Examples of asymmetries and ratio are constructed in order to minimise theoretical and experimental uncertainties. In view of the foreseen measurements of this type of process at the LHC Run-II, these results constitute a baseline for the characterisation of the expected factorisation breaking effects.