Threshold resummation of the rapidity distribution for Drell-Yan production at NNLO+NNLL

We consider the production of pairs of lepton through the Drell-Yan process at the LHC and present the most accurate prediction on their rapidity distribution. While the fixed order prediction is already known to next-to-next-to-leading order in perturbative QCD, the resummed contribution coming from threshold region of phase space up to next-to-next-to-leading logarithmic (NNLL) accuracy has been computed in this article. The formalism developed in [1-3] has been used to resum large threshold logarithms in the two dimensional Mellin space to all orders in perturbation theory. We have done a detailed numerical comparison against other approaches that resum certain threshold logarithms in Mellin-Fourier space. Our predictions at NNLL level are close to theirs even though at leading logarithmic and next-to-leading logarithmic level we differ. We have also investigated the impact of these threshold logarithms on the stability of perturbation theory against factorisation and renormalisation scales. While the dependence on these scales does not get better with resummed results, the convergence of the perturbative series shows a better trend compared to the fixed order predictions. This is evident from the reduction in the K-factor for the resummed case compared to fixed order. We also present the uncertainties on the predictions resulting from parton distribution functions.

Introduction.-TheStandard Model (SM) has been extremely successful in describing the physics of elementary particles.The production of oppositely charged lepton-pairs, known as the Drell-Yan production [4], is one of the benchmark processes to probe physics at TeV energies at the colliders, namely earlier at Tevatron and now at the Large Hadron Collider (LHC).Because of its large cross section and small systematic uncertainties, Drell-Yan production also serves as luminosity monitor [5] at the LHC.Most importantly, at hadron colliders, Drell-Yan production provides valuable information about the partonic structure of hadrons, its clean electromagnetic probe is best suited for the search of any new physics beyond the SM (BSM).An excess rate over the SM in this channel will potentially indicate the signature of BSM physics.Drell-Yan is the potential background for processes involving Z or W and also for spin-2/graviton searches.With the current LHC-13 and upcoming LHC-14 runs, more events will be available to precisely study the Drell-Yan distributions over a wide kinematic region.
Due to its undeniable importance, Drell-Yan has been studied theoretically to a great extent over many decades [6,7].There has been continuous efforts towards the computation of higher order QCD and electroweak (EW) corrections to unprecedented accuracy.The full inclusive production cross section is known up to nextto-next-to-leading order (NNLO) [8,9] for a very long time.Very recently, predictions at next-to-next-to-nextto-leading-order (N 3 LO) level have become available considering only the dominant soft-virtual (SV) contributions [10,11].Electroweak corrections beyond leadingorder (LO) are also known and at next-to-leading-order (NLO) level, they were computed in [12,13].While inclusive production is important for precise prediction of cross section, differential distributions allow a wider comparison with experiments.Fully differential distributions such as rapidity, transverse momentum for the Drell-Yan are known to up to NNLO level in QCD [14][15][16][17][18]. SV contributions for the rapidity distribution are also known at N 3 LO level [2,19,20].Studies where both QCD and EW corrections are combined can be found in [21].Parton showers matched with NLO QCD results for the Drell-Yan are also available in MC@NLO [22], POWHEG [22,23] and aMC@NLO [24] frameworks.
One of the differential distributions that has been studied extensively is the transverse momentum (p T ) distribution of pair of leptons or the vector bosons such as Z/W ± , see [25][26][27][28][29], often in their large p T region.The rapidity distribution in Drell-Yan was computed in [25] at NLO level in QCD and it was then extended to NNLO level in [14,17] which stabilise the predictions [15] giving only a few percentage sensitivity to renormalisation and factorisation scales, say for example at the Z mass region.However it has to be noted that the result does vary significantly w.r.t. the choices of different parton distribution functions (PDFs).In particular, at large invariant mass or at large rapidity of the final state, the cross sections are sensitive to large Björken x regions of PDFs, where different PDFs show not only differences between them but also exhibit large uncertainties.For a recent review see [30].This sensitivity of PDFs will in turn help to constraint the PDF sets much better.Hence, it is important to study these distributions.
The fixed order predictions are often not reliable in certain regions of phase space where large logarithms of some kinematic variables can appear.For example, at the partonic threshold i.e.where the initial partons have just enough energy to produce the final state such as a pair of leptons or Z/W ± and soft gluons, the phase space available for the gluons become severely constrained which results in large logarithms.These large logarithms however can be systematically resummed to all orders in perturbation theory for reliable predictions.This has made the resummation program an important topic of investigation over many years.For the inclusive production, the resummation of soft gluons in the threshold region was established [1,[31][32][33][34][35][36] in the Mellin space and for the transverse momentum distribution, at small p T , the resulting large logarithms were shown to exponentiate in the impact parameter space [37,38].A modern approach based on soft-collinear effective theory (SCET) demonstrates similar resummation in momentum space, see [39] for inclusive production and [40] for transverse momentum distribution.Resummation for the differential distribution with respect to the Feynman variable x F which describes the longitudinal momentum of the final state was studied in [1] and it was found that there were two thresholds and both could be resummed to all orders.For the resummation w.r.t.x F using a different scheme see [41].For the rapidity distribution, resummation similar to the inclusive one, with a single scaling variable, can be obtained in certain kinematic regions, see [42][43][44][45][46] and an equivalent approach based on SCET can be found in [47,48].In the former one, called the standard direct QCD (dQCD) approach [1,31,32], since the resummation is performed in Mellin space where the phase space of the soft gluons factorises under appropriate Mellin transformation, the threshold limit of the partonic scaling variable z → 1 corresponds to Mellin variable N → ∞, where z = q 2 /ŝ, q 2 = M V , V = Z, W ± and ŝ is the partonic centre of mass energy.In SCET approach [39,47,48], however, resummation can be performed both in Mellin space as well as in z-space using the evolution operators of soft and the hard functions of the coefficient function.
Resummation of large logarithms for rapidity distribution has been an interesting topic and several results are already available to very good accuracy.In [44] the authors have studied resummation of rapidity W ± in Mellin-Fourier (M-F) space following a conjecture (see [42]) and later on same approach was used for Drell-Yan production in [45].A more detailed study in the context of W ± productions as well as production of a pair of leptons was undertaken in [49] emphasising the role of prescriptions that take care of diverging series at a given logarithmic accuracy.
In this article we will follow the dQCD approach [1] to study soft gluon resummation for the rapidity distribution of a pair of leptons produced in hadron colliders.Recently, using the formalism developed in [2,19], we [3] derived a general result, applicable to production of any colorless state in hadron colliders, that resums the soft gluons to all orders in perturbation theory in two dimensional Mellin (M-M) space spanned by N 1 , N 2 .We also investigated their numerical impact on the rapidity distribution of the Higgs boson produced at the LHC.The soft gluon effects show up through delta functions and plus distributions in the partonic cross sections when the partonic scaling variables reach the threshold limits, i.e., z 1 → 1 and z 2 → 1 and these contributions can be resummed to all orders both in z 1 , z 2 space and in N 1 , N 2 space.These resummed results were expanded to the desired accuracy to obtain fixed order predictions for various observables [19,20,50] at the LHC in the SV approximation.This double threshold limit, denoted by a pair of limits, namely (z The corresponding large logarithms are of the form ln n (N i ), where n = 1, • • • and i = 1, 2 and the resummation in M-M space resums terms of the form ω = a s β 0 ln(N 1 N 2 ) through a process independent function g(ω) and a process dependent but N i independent function g 0 .Here β 0 is the leading coefficient of the beta function of the strong coupling constant g s and a s = g 2 s (µ 2 R )/16π 2 with µ R being the renormalisation scale .
The main goal of this article is to study the numerical impact of resummed contributions in the M-M approach on the fixed order predictions for the rapidity distribution of a pair of leptons in the Drell-Yan process at the LHC.At NNLO level, fixed order results show remarkable stability against the factorisation and renormalisation scales.In addition, they demonstrate excellent perturbative convergence.While this is a good news for any phenomenological study with Drell-Yan process, the question remains whether the fixed order predictions will be plagued by presence of large kinematic logarithms resulting from soft gluons in the threshold regions at every order in perturbative expansion.The formalisms that can resum these large logarithms to all orders do exist and it is not a priori clear whether the resulting resummed contributions will not affect the fixed order predictions.Hence, a detailed study taking into account these threshold effects through resummation is warranted.In addition, owing to various ways by which these logarithms can be resummed, a detailed comparison of these approaches is desirable.This article attempts to address all these issues.We start by recapitulating the resummation formalism based on M-M approach and then proceed with a detailed numerical study at the LHC and conclude with our findings.
Theoretical framework.-Inthe QCD improved parton model, for the production of a pair of leptons with invariant mass q 2 and rapidity y, the double differential cross section can be written as where σ q B (x 0 1 , x 0 2 , q 2 ) is the Born prefactor, τ = q 2 /S = x 0 1 x 0 2 with q being the momentum of the final state lepton pairs and S = (p 1 + p 2 ) 2 where p i are the momenta of the incoming hadrons.The hadronic rapidity is defined as y = 1 2 ln p2.q p1.q = 1 2 ln are the PDFs having momentum fractions x 0 2 /z 2 respectively, renormalised at the factorisation scale µ F .∆ q d,ab (a s , z 1 , z 2 , q 2 , µ 2 F ) (shorthand as ∆ q d,ab (z 1 , z 2 )) on the other hand is the Drell-Yan coefficient function for the rapidity distribution mass factorised at µ F .Unlike PDFs, these are calculable order by order in QCD perturbation theory in powers of a s .The coefficients in this power series expansion contain distributions such as δ(1−z i ) and ln m−1 (1−zi) n being the order of perturbation and regular functions of z i .The former ones, namely the distributions, constitute the SV part, denoted by ∆ q,SV d,ab while the latter one, the hard part is ∆ q,H d,ab .In the SV part, these distributions result from certain regions of the phase space in the real emission sub processes and also of the loop integrals in the virtual ones.While these distributions are singular as z i → 1, they are integrable functions.At the level of hadronic cross sections, they often dominate over the hard part when folded with the appropriate PDFs, in the above mentioned kinematic regions at every order in perturbation theory.Hence, they can potentially disturb the reliability of the perturbative predictions.The resolution is to resum these large terms, often the logarithms, to all orders to obtain any sensible prediction.This is indeed the case for the rapidity distribution in Drell-Yan when the scaling variables z 1 → 1 and z 2 → 1. Recall that in the work by Catani and Trentadue [1], a different distribution namely Feynman x F for the Drell-Yan was studied in the context of threshold resummation and it was shown that the potential threshold logarithms can be resummed to all orders in perturbation theory working in M-M space.They had established that these logarithms could be exponentiated and also obtained the resummed result at the next-to-leading-logarithmic (NLL) accuracy.Following this, in [2,19,20] we demonstrated for the rapidity distribution of any colorless particle, resummation of the distributions defined with respect to the scaling variables z i to all orders in perturbation theory in z 1 , z 2 space and later extended it to N 1 , N 2 space in [3] by applying two dimensional Mellin transformations on these distributions to obtain resummed result in the M-M space.The latter one, namely the resummed result in the M-M space, turned out to be more suitable for numerical study and hence, we used this approach to demonstrate the importance of these threshold logarithms for the rapidity distribution of the Higgs boson at the hadron collider to next-to-next-to-leading-logarithmic (NNLL) accuracy [3] over NNLO.In this paper, we explore this approach to study the Drell-Yan process at the LHC.We use the general result obtained for any colorless particle in [3] to study the numerical impact on the rapidity distribution of a pair of leptons produced at the LHC.
Note that the approach followed here [3] differs from earlier ones (see [42,43,45,[47][48][49] in the way the threshold limit(s) is(are) taken).In the latter approach, the threshold contributions from soft gluons in the partonic cross section are defined by considering only those distributions w.r.t. the scaling variable z = z 1 z 2 which appear in the region when z → 1.The remaining contributions contain not only regular terms in z but also distributions and regular functions of partonic rapidity variable (y p ).Here, only distributions in z are resummed to all orders treating the remaining terms as hard part.Interestingly, if one works in M-F space, it can be easily shown that in the limit z → 1, the threshold logarithms resulting from N → ∞ are identical to those of the inclusive cross section.Thus the resummation for the rapidity distributions for the Drell-Yan has been done using a single Mellin variable N corresponding to z and keeping the y p dependent coefficients as it is.
We employ the technique developed in [3] namely the M-M space approach to perform the soft gluon resummation for Drell-Yan rapidity distribution.Thanks to the convolution structure of the hadronic cross section in terms of the PDFs f a,b and the Drell-Yan coefficient functions ∆ q d,ab , the two-dimensional Mellin transformation of the Born normalised hadronic cross section becomes a simple product of fa (N 1 ), fb (N 2 ) and ∆q It was shown in [2,3,19] that the SV part of ∆q d,ab ( ∆SV d,q ) exponentiates the threshold logarithms through the cusp anomalous dimension A q and the collinear functions D q d giving the resummed result where z i = (1 − z i ).The cusp anomalous dimensions for the quark, A q , are known for Drell-Yan up to 3-loops [51][52][53][54][55][56], and the coefficients in A q = ∞ i=1 a i s A q i are given by and D q d s are related [2,3,19] to the D q s of the inclusive cross section for the Drell-Yan.Expanding and n f is the number of active flavours.In [3], following [57], we systematically computed the two dimensional Mellin transformations in the large N i limits and the result takes the following form: where ω is defined as ω = a s β 0 ln(N 1 N 2 ), with The coefficients g d (a s , ω) are process independent and contain purely logarithmically enhanced terms and can be expanded as, Rescaling the constants by appropriate The first three coefficients of the QCD β function, β 0 , β 1 and β 2 are given by [58] The N 1 , N 2 independent terms resulting from integrals have been absorbed in gq d,0 addition to g q d,0 , which however depends on the specific process under study.In principle, these N 1 , N 2 independent terms can also be exponentiated.For Drell-Yan production, expanding , we find up to a 2 s order in strong coupling: To study the numerical impact of our resummed result, we require in addition fixed order results containing only the large logarithms to perform proper matching.They can be obtained by truncating the resummed result and in the following, we present the ∆SV(i) F up to NNLO level: where ω = N 1 N 2 .In the following, we will discuss how these resummed contributions can be systematically included in order to study their phenomenological importance at the LHC.
Phenomenology.-Our next task is to include the resummed contributions consistently in the fixed order predictions and study their numerical impact on the rapidity distribution of lepton pairs produced in the Drell-Yan process at the LHC.We consider the production of both leptons, i.e. + − , where = e, µ through Z and γ * in the collision of two hadrons at the centre of mass energy 14 TeV.Unless otherwise stated, we will mostly focus on the region containing the Z-pole.We take n f = 5 flavors, the MMHT2014(68cl) PDF set [59] and the corresponding a s (M Z ) through the LHAPDF-6 [60] interface at each order in perturbation theory.For the fixed order rapidity distribution, we use the publicly available code Vrap-0.9[14,61].The resummed contribution is obtained from ∆SV d,q (N 1 , N 2 ) in Eq. ( 6) after performing Mellin inversions which are done using an in house Fortran based code.Since the resummed result cannot be simply added to the fixed order one because all the ln(N i ) and N i independent terms present in the resummed exponential g q d and gq d,0 are already present in the fixed order results and hence care is needed to avoid double counting.This can be achieved simply by employing a matching procedure at every order.The matched result is given below where σ q B is given by with α = α(M Z ) = 1/127.925,e q is the quark charge, 03363 and B Z µ = 0.03366.The first term in Eq. ( 12), d 2 σ q,f.o /dq 2 dy , corresponds to contributions resulting from a fixed order perturbative computation.The second term on the other hand contains only threshold logarithms ln(N i ) but to all orders in perturbation theory.The subscript "trunc" in the ∆SV d,q indicates that it is truncated at the same order as the fixed order one after expanding in powers of a s .Hence, at a given order n in a s (i.e., at order a n s ), the non-zero contribution from the second term starts at order a n+1 s and includes terms from all orders.For the fixed n-th order contribution, namely N n LO, the contribution from the second term is called N n LL .Hence, we use the notations LO, NLO and NNLO for the fixed order predictions and correspondingly LO+LL, NLO+NLL and NNLO+NNLL for the resummed ones.It is well known that the resummed expression diverges due to the missing non-perturbative contributions.These divergences show up when ω → 1 in the functions g q d,i ; they are due to the coupling constant a s (µ 2 R ) that diverges near the Landau pole.In order to resolve this, we have adopted the Minimal Prescription (MP) [62].The contours for the integrals corresponding to two Mellin inversions [63] are chosen in such a way that all the poles in the complex plane spanned by N 1 , N 2 remain to the left of the contours except for the Landau pole.
Since the leading order contribution to the Drell-Yan process is due to EW interactions, the dominant theoretical uncertainty comes from the factorisation scale µ F that enters through the parton distribution functions while the dependence on the renormalisation scale µ R starts only from NLO onwards.Unlike the leading order prediction, in the resummed case, the LL contributions do depend on µ R through ω in g q d,1 (w) given in Eq. ( 8).Hence, µ R dependence will show up at even LO+LL level.This will be evident from Fig. 3 that one finds larger scale uncertainty from LO+LL contributions compared to the fixed order one at LO level.It is then important to understand the impact of these two scales at the resummed level and also to determine the optimal choice for the central scale around which the scale uncertainty remains minimal.For the fixed order case it has already been realised in [64] that the optimal choice for the central scale is when both µ R and µ F are set to M Z .In order to obtain the optimised central scale for the resummed case, we have plotted in Fig. 1 the dependence of the rapidity distribution on a The symmetric band is obtained by performing 7-point scale variation [49,57,64] around a given central scale with the constraint (k and by taking maximum absolute deviation from the central scale.From the first and the last panels of Fig. 1, it is clear that the optimal central scale choice is (M Z , M Z ) whereas the middle panel favours ( M Z 2 , M Z ) for the central scale.Comparing all three panels, we find that the choice ( M Z 2 , M Z ) gives minimum uncertainty band.This is demonstrated in Fig. 2 where we compare predictions from fixed order using the central scale (M Z , M Z ) against those from resummed result using the central scale ( M Z 2 , M Z ) for two rapidities y = 0 and y = 2.4.The scale uncertainties from the resummed case at NNLO+NNLL are comparable to what one obtains from NNLO.However, the central values at NLO+NLL and NNLO+NNLL are very close to each other compared to those of fixed order results demonstrating the better perturbative convergence.
As we have discussed in the introduction, in [44,45,49] resummation of threshold logarithms for the rapidity distribution was achieved in the M-F space.Our formalism [3] differs from the other approach in the way the thresh- old contributions are resummed.We resum large logs resulting from the regions where scaling variables z 1 and z 2 approach unity simultaneously while in the case of M-F, only large logarithms from the region where the partonic threshold variable z approaches unity and the partonic rapidity y p is zero, are resummed.In the following we will make the numerical comparison of our predictions, namely the M-M formalism against those of M-F reported in [49].The fixed order contributions are obtained by using Vrap-0.9[14,61]; the resummed contributions up to NNLL for M-F are obtained by using publicly available code ReDY [65] and for M-M, we use our in house Fortran routine.We have set all the parameters including the PDF set (NLO set of NNPDF-2.0 [66] at every order) same as those used in [49].Both our results and those from ReDY are listed in the Table I for various scale choices at the central rapidity.At LL level, both M-F and M-M give positive contributions but the contribution from M-M is about three times larger compared to M-F independent of the scale choice.The additional contribution over LL at NLL for M-F is negative for some scale choices and positive for the rest while for M-M, it is always negative.The magnitude of these additional contributions for M-M is larger than M-F.Interestingly, at NNLL level, the additional contributions over NLL for M-F and M-M are both negative in a such a way that the net NNLL contributions from both approaches become comparable.In the case of M-F, the NLO+NLL is 2% larger compared to LO+LL and NNLO+NNLL is -4.7% larger compared NLO+NLL.For M-M, the corresponding ones are -0.8% and -4.9% respectively at µ R = µ F = 2M Z .
In Fig. 3, using Eq. ( 12), we present the cross section for producing lepton pairs as a function of the rapidity y up to NNLO in the left panel and to NNLO+NNLL in the right panel along with the respective K-factors.The Kfactor at a given perturbative order, say at N n LO (N n LO + N n LL), is defined by the cross section at that order normalised by the same at LO (LO+LL) at the central scale µ R = µ F = M Z .We have made this choice for the scales because the fixed order perturbative prediction is well behaved around this scale [64].The symmetric band at each order is obtained by varying µ R and µ F between , by adding and subtracting to the central scale the highest possible uncertainties originating from all the scale combinations.We find that the magnitude and the sign of the resummed contribution are sensitive to the order of perturbation as well the exact values of y and the scales µ R , µ F .For example, if we choose µ R = M Z /2 and µ F = M Z instead of µ R = µ F = M Z as the central scale, we obtain a negative contribution from NNLL terms for all values of rapidity.Fig. 3 also demonstrates that the inclusion of N n LL contributions increase the cross section at every order for a wide range of rapidity values.In addition, the overlap among various orders is larger for the resummed case compared to the fixed order ones, because the uncertainty band at each order in the resummed case is bigger compared to fixed order.As far as fixed order results are concerned, in particular at NNLO level, several partonic channels open up, effectively reducing the scale uncertainty considerably.On the other hand, resummed contributions come only from quark anti-quark initiated channels to all orders in perturbation theory as other channels do not give threshold logarithms of the type that is resummed.We confirm this through Fig. 4, where we have studied the effects of resummation over the fixed order contributions, by considering a) only q q channel at NNLO and b) all the channels at NNLO.We perform our analysis for y = 0 and set µ F = M Z while varying µ R between M Z /2 to 2M Z .For the q q channel the resum contributions arising from the two extreme scales are of opposite sign and their individual contributions are such that the NNLO+NNLL (q q) curve shows a stable behaviour as compared to NNLO (q q).While the fixed order decreases by 2.36% from M Z /2 to 2M Z , the corresponding decrease for NNLO+NNLL (q q) is 1.53%.This confirms the reduction of scale dependence upon adding resummed terms to the fixed order contributions.However the scenario entirely reverses when we consider all the channels at NNLO.We find that the differential cross section at NNLO (all) increases by 0.29% in the entire range of µ R values; the corresponding increase for NNLO+NNLL (all) is 1.29%.This reduction of the scale dependence at NNLO is due to cancellations among different partonic channels.However the resummation effects come only from q q channel which adds to the fixed order in such a way that the resummed uncertainty increases.This explains the increase of the scale uncertainty at each resummed order depicted in Fig. 3. Furthermore an incomplete cancellation of the factorization scale dependence against the PDFs which do not contain resummed threshold logarithms also increases the band.For the fixed order, the K-factor at NLO varies between 1.3 and 1.2 and at NNLO between 1.37 and 1.3 over the entire rapidity region.On the other hand, the K-factors at both NLO+NLL and NNLO+NNLL significantly overlap with each other over most of the regions of rapidity and stay around 1.2.This demonstrates a better perturbative convergence for resummed case compared to the fixed order.In Table II, we have presented the cross section for benchmark rapidity values along with the percentage scale uncertainties.Note that the differential cross-section at NNLO+NNLL level for the central scale is well approximated by the same at NLO+NLL.In fact, NNLO+NNLL increases approximately by 0.8% with respect to NLO+NLL; the corresponding number for NNLO over NLO is approximately 3%.From the trend that resummed results give, we an- FIG. 4. Drell-Yan rapidity distribution for 14 TeV LHC at y = 0 using MMHT PDFs.The variation of fixed order and resummed results as a function of µR are shown separately for q q channel and also for all the channels added together.
ticipate N 3 LO+N 3 LL cross-section will fall completely within the NNLO+NNLL band.
In Fig. 5, we have plotted both fixed order and resummed results at various orders for the larger invariant mass at q = 1 TeV.Interestingly, the uncertainty bands at NLO+NLL and NNLO+NNLL levels are better compared to those from fixed order.Also the predictions at various orders are closer compared to those from fixed order demonstrating better perturbative convergence of the higher order predictions from the resummed terms.In fact the resummed K-factor for the central rapidity at NNLO+NNLL is 1.25 compared to 1.39 at NNLO.As there are several PDF groups in the literature, each providing sets of PDFs, it is customary to estimate the uncertainty resulting from the choice of PDFs within each set of a given PDF group.Using PDFs from different PDF groups namely MMHT2014nnlo68cl [59],   ABMP [67], NNPDF3.1 [68] and PDF4LHC [69] we have obtained the cross sections along with the corresponding PDF uncertainty.In Fig. 6, we have plotted the uncertainty bands for various PDF sets as function of rapidity in order to demonstrate the correlation of PDF uncertainty with the rapidity values.This will help to better constrain the PDF fits using measurements on rapidity in the Drell-Yan process.In Table III, we have also tabulated the cross sections along with % uncertainties resulting from the choice of different PDFs.
We have studied the q-integrated rapidity distribution at the LHC with 8 TeV centre of mass energy at NNLO+NNLL.The invariant mass is integrated between 60 GeV and 120 GeV and choose µ R = µ F = M Z .Unlike our earlier analysis, we have included both e + e − as well as µ + µ − final states.NLO EW corrections are also included as they are also comparable to QCD corrections at NNLO+NNLL.The EW contributions are obtained by using the publicly available code Horace-3.are taken to be m e = 0.51099 MeV and m µ = 0.10566 GeV respectively.The NNLL contribution increases the cross-section by roughly around 0.5% with respect to   NNLO.However the EW corrections at NLO give negative contribution to the cross-section.The corrections are different for e + e − and µ + µ − pairs.For electrons, the EW contributions are twice that of muons.In total from the electron and muon channels, we see an overall 2.3% decrease in the cross section w.r.t. the NNLO in the central rapidity region.The rapidity distribution in Fig. 7 being inclusive in transverse momenta of the final state leptons, can not be directly compared with the results presented in [74] where a minimum transverse momenta cut is applied in the selection of final state leptons.To really compare this one needs distributions exclusive of transverse momenta which at the moment beyond the scope of the current paper and we leave it to future work.
Both at Tevatron and at the LHC, there are already precise measurements of rapidity distributions for different ranges of invariant mass q.For one of the earliest set of measurements see NuSea [76,77].Since the data from the LHC depends heavily on the kinematic cuts of the final states we cannot directly compare against our predictions.On the other hand CDF [75] has data for the rapidity distributions for wide range of y with invariant mass range 66 < q < 116 GeV.In Fig. 8, we have com- pared our predictions against the data at √ s = 1.8 TeV after integrating q between the above mentioned range for two different choices of PDF sets.The scale uncertainty is obtained as before by using 7-point scale variations around the central value µ R = µ F = M Z .We note that at NNLO+NNLL level, the resummed contributions over the fixed order is very mild, less than 0.5%.We have also observed that the resummed effects become significant for large invariant mass regions.
Discussion and Conclusion.-In this article we have done a detailed study on the role of resummed threshold logarithms for the rapidity distribution of pairs of leptons in the Drell-Yan process at the LHC.Being one of the cleanest channels at the hadron colliders, precise measurements of various observables such as inclusive cross section, transverse momentum and rapidity distributions are already available.Precise predictions from perturbative QCD are known to NNLO level for long and corrections from electroweak theory have become available in recent times.The latter effects being close to the second order effects from QCD so that dedicated efforts to understand the EW effects have been undertaken.Owing to the dominant QCD interactions, soft gluons play vital role in most of the observables.They show up in certain kinematic regions through large logarithms in the perturbative computations.Often they spoil the reliability of the fixed order predictions.In this present article, we have made a detail study on the effect of these soft gluons within the resummation framework.In the literature two different approaches exist.They differ in the kind of logarithms that are resummed to all orders.The approach which uses Mellin-Fourier transformation to achieve the resummation to resum large logarithms of the scaling variable z has been well studied for the rapidity distribution.Threshold logarithms resulting from regions where the scaling variable z i approach unity are successfully resummed using M-M approach.We used the latter approach to get the quantitative predictions at NNLO+NNLL level.Since these formalisms resum different type of logarithms to all orders, they are expected to give different numerical predictions.In this article, we have not only undertaken a detailed study on the numerical impact of M-M approach for the first time for the Drell-Yan process but also made a detailed numerical comparison against the M-F approach.While at LL and NLL level, they differ very much, surprisingly at the NNLL level both the approaches converge to a few percent correction to the fixed order prediction.This could be accidental, however it is desirable to understand this coincidence at NNLL level.Our numerical study on the dependence of renormalisation and factorisation scales shows that the optimal central scales for the resummed result are µ R = M Z /2 and µ F = M Z while it is µ R = µ F = M Z for NNLO.We have also found that, for wide range of rapidity, the scale uncertainties from NNLL contributions at every order are slightly larger than those from fixed order results.We believe that this could be due to an incomplete cancellation of scale dependent terms between resummed result and the PDFs.Note that the PDFs that we use are extracted from data using the fixed order perturbative predictions for the observables and also using evolutions equations controlled by splittings functions computed to desired order in strong coupling constant.Hence, we expect that there will be a better cancellation of scale if appropriate resummed PDF sets are available.We have also presented our predictions for various choices of PDFs from various PDF groups.Each group has several sets and hence we have not only made comparisons with respect to various groups but also estimated the uncertainty from different sets within each PDF group.We have also predicted the q integrated rapidity distribution.Since our resummation formalism can not take into the experimental cuts such as the transverse momentum and/or polar angles of the final state leptons, we can not make any direct comparisons with the existing data on the q integrated rapidity distribution measured at the LHC which are extracted after employing cuts on transverse momentum of final state leptons.On the other hand we have compared our predictions against CDF data at Tevatron for the invariant mass range 66 < q < 116 GeV and found good agreement within both theoretical and experimental uncertainties.We believe that perturbative results that take into account both fixed order as well as the resummed contributions will provide a precise determination of PDFs from the ample data that are already available at the LHC.

4 FIG. 2 .
FIG.2.Fixed order predictions with the central scale µR = µF = MZ and resummed prediction with the central scale µR = MZ /2, µF = MZ for rapidity y=0 and y=2.4 using MMHT2014 PDF at each order.The uncertainties are obtained by using 7-point scale variation (see text for more details) around the central scale.

FIG. 3 .
FIG.3.Drell-Yan rapidity distribution for 14 TeV LHC at q = MZ using MMHT PDFs.The fixed order results are plotted in the left panel and the resummed results in the right panel.Central scale is chosen as µR = µF = MZ for both and the corresponding bands are obtained using 7-point scale variation (see text for more details) around the central scale.The lower panel represents the corresponding K-factors.
Fixed order and the resummed cross sections with % scale uncertainties along with the K-factors at the central scale µR = µF = MZ .

FIG. 6 .
FIG.6.PDF variation at NNLO+NNLL using various sets.The y-axis represents the ratio of extremum variation over the central PDF set.

FIG. 7 .
FIG.7.Rapidity distribution at NNLO+NNLL for 8 TeV LHC in the invariant mass range 60 < q < 120 GeV.The dotted line is the fixed order NNLO contribution, the dashed line represents NNLO+NNLL result and the solid line includes EW corrections.

FIG. 8 .
FIG.8.Comparison between resummed results and the CDF data[75] at √ s = 1.8 TeV in the invariant mass range 66 < q < 116 GeV for two different PDF sets.