PDF sensitivity studies using electroweak processes at LHCb arxiv:1011.4260v1 [hep-ex] 18 Nov 2010 University College Dublin E-mail: francesco.de.lorenzi@cern.ch We describe parton density function sensitivity studies, using muon final states produced through the Drell-Yan process via W, Z or γ down to a Q 2 of 10 GeV 2. This makes use of LHCb s unique ability to trigger on low transverse momentum objects. Due to the forward acceptance of LHCb, x values down to 2x10 6 can be probed, where with just 100 pb 1 of data, the gluon PDF can improved up to 70%. XVIII International Workshop on Deep-Inelastic Scattering and Related Subjects, DIS 2010 April 19-23, 2010 Firenze, Italy Speaker. The author would like to acknowledge the support of Irish Research Council for Science, Engineering and Technology. c Copyright owned by the author(s) under the terms of the Creative Commons Attribution-NonCommercial-ShareAlike Licence. http://pos.sissa.it/
1. Introduction LHCb [1] is a detector designed to study b-physics so its pseudo-rapidity coverage extends in the forward region where B mesons are dominantly produced (1.9 < η < 4.9). Part of that range is in common with ATLAS and CMS which have a coverage of η < 2.5; the rest is unique to LHCb. During the first year of data taking, LHC will run at the center of mass energy of 7 TeV and an integrated luminosity of 1 f b 1 will be delivered to LHCb in the first year of data taking. A unique feature of LHCb is its ability to trigger and reconstruct very low transverse momentum muons (p T > 1GeV ). Besides b-physics goals LHCb is able to perform precise measurement of electroweak channels with muon final states at high rapidities. We will show that these measurements can probe parton distribution functions for different values of Bjorken-x. 2. LHCb parton kinematics The current knowledge of the PDFs arises from fitting procedures to HERA, fixed target experiments and Tevatron data (or a subset of those). The foreseen electroweak measurements at LHCb include the analysis of the Z and W ± production and the Drell-Yan pair production of muons in four different mass bins (2.5 5 GeV, 5 10 GeV, 10 20 GeV and 20 40 GeV )[2]. Production at high rapidity involves the interaction of one valence quark carrying a large fraction of proton momentum (high-x) and a sea quark with small fraction of momentum (very small x); according to the relation that link x and Q 2 x 1,2 = Q s e±y Figure 1: Percentage uncertainty on production cross-section of electroweak boson as function of rapidity where s is the energy in the center of mass of the collisions. For the Drell-Yan production with Q in the region 2.5 5 GeV the lowest value of x reachable is 2x10 6. Hence due to its forward rapidity coverage and the ability to reconstruct low invariant mass muon pair systems, LHCb will cover an un-explored region of kinematical space where Parton Distribution Functions are known with less precision. This aspect can be also deduced looking at the theoretical uncertainties on the rapidity distribution for W ±, and Z/γ ; The largest contribution to these uncertainties arises from PDFs uncertainties. In figure 1 the percentage uncertainties on the production cross-sections for W ±, Z and Drell-Yan process at three different mass are reported: the uncertainties due to PDF become large at high rapidity for all the channels studied (corresponding to rapidity LHCb acceptance) and for the Drell-Yan pair production, the smaller the mass, the larger the uncertainty. 3. The fitting algorithm PDFs are obtained via global fits to experimental data by different working groups, and their 2
results have some model-dependence that may differ between each other. Theoretical predictions for the PDF sets MSTW08[3], CTEQ66[4], Alekhin[6] and NNPDF2.0[5] are used in this analysis. The first three groups perform a global fit based on a Hessian method while the fourth uses a neural network approach based on Montecarlo replicas. We now briefly describe our method to constrain the PDFs where the PDFs set have been obtained using these two approaches. 3.1 PDFs built using a Hessian method The Hessian method is used by the groups MSTW, CTEQ and Alekhin. A polynomial parametrization of the PDFs such as x f (x,q 0 ) = (1 x) η (1 + εx 0.5 + γx)x δ is chosen. The global minimum of their fits to the global data produces a central value with the uncertainties given by a N N covariance matrix. Diagonalizing the matrix gives a set of N independent variables. In this way a representation of the PDF in a space where uncertainty can be treated as independent with a set of N eigenvectors is provided. A given physical observable, for example the rapidity distribution, f (y), can be described as a linear combination of the observable obtained using the PDF best-fit, f 0 (y), and the observable f i (y) obtained by moving the i th eigenvector by 1-sigma. f (y) = f 0 (y) + N i=1 [ f 0 (y) f i (y)]λ i (3.1) The one sigma bound is mapped out by sampling λ i from a unit Gaussian with mean zero; this function also provides a form with which the data can be fitted. Given binned data with N j events in the j th bin with error δ j, the normalization and the shape of the distribution can be fitted by minimizing χ 2 = #bins j [ Nj λ 0 ( f0 + N i=1 λ i( f 0 f i ) ) δ j ] 2 + λi 2 λ>0 where λ 0 is the overall normalization or luminosity and λ i are the fitted eigenvalues, constrained by the second term to follow a Gaussian probability. Before the fit the eigenvalues by definition of the model have an uncertainty equal to 1. After our fit, values smaller than 1 imply that we have a better constraint on the PDF. The hessian method provides a diagonalized covariance matrix of the parameters. This fitting operation provides a new covariance matrix, V, in which values on the main diagonal are smaller than 1 and terms off diagonal can be different from zero. This matrix can be used to calculate uncertainties on PDFs using the usual equation for correlated variables: δf = i j F F V i j λ i λ j where λ i are the fitted parameters and F is a PDF (e.g.e u V, d V, sea, gluon etc.). 3.2 PDFs built using a Neural network method NNPDF uses a neural network approach. PDFs are parametrized with a redundant flexible neural-net (i.e. a convenient functional form with 200 parameters) This technique produces a set of equally likely Montecarlo replicas of the PDFs. The best estimation of any observable of the 3
PDFs is obtained by generating this observable using each of the replicas and averaging: f [q] = 1 N N k=1 f [ q (k)] where N is the number of replicas, f is a generic function of the PDF set and ( q (k) are PDF replicas. The uncertainty can be expressed through the standard deviation. σ 2 f = 1 ( [ N 1 N k=1 f q (k) ] f [q] ) ) 2. If more than one observable is generated, it is even possible to compute correlations between quantities. To evaluate the PDF constraint with this PDFset we generated pseudo-experiment as average over the replicas then we generated rapidity distributions with each replica. We compared the prediction of a single replica with the pseudo-experiment minimizing the χ 2 ( χi 2 Nbin λ 0 f i (y) = bin bin where i is the index of the replica, N bin is the pseudo-experiments result in that rapidity-bin and f i (y) is the value of the rapidity distribution generated with the i th replica. This algorithm associates to each replica of the PDFs a χ 2 and so the degree of compatibility of each replica can be evaluated as the chi-square probability. Now it is possible to constrain the PDFs considering the χ 2 -probability as the conditional probability for measuring the set of NNPDF replicas considering the new LHCb data [7]. We don t consider all the PDF replicas as equally likely anymore but we use the χ 2 - probability as a weight in the central value and uncertainties calculations. ) 2 4. Results and comments We present the expected improvement on u-valence, d-valence, sea, s + and gluon distributions when introducing 1 f b 1 of LHCb electroweak data. Here we present results obtained using the PDF set MSTW08. Results for the other PDF sets can be found in [8]. In figure 2 the ratio between the PDF uncertainty after and before the fitting procedure is plotted. As expected the largest improvement is given by the lowest mass bin of the Drell-Yan process. With 1 f b 1 of LHCb data the sea and the gluon distributions at x 10 5 and at high-x can be improved up to 70% while in mid-x regions the improvement is about 50%. Other PDFs can be improved on the full x range of about 60%. For the W ± and Z production the best result is obtained when performing a combined fit of the three rapidity distributions. About 30% improvement can be expected over the whole x range for each of the PDFs. References [1] LHCb Collaboration, "Reoptimised Detector Design and Performance " TDR, JINST3: 508005, 2008; [2] J. M. Keaveney, "Studies of electroweak boson production in the forward region with LHCb" in the same proceedings; [3] A.D. Martin, W.J. Stirling, R.S. Thorne, G. Watt, "Parton distributions for the LHC", arxiv:hep-ex/0901.0002v2; [4] J. Pumplin, D.R. Stump, J. Huston, H.L. Lai, P. Nadolsky, W.K. Tung, "New Generation of Parton Distributions with Uncertainties from Global QCD Analysis", JHEP0207:012,2002; 4
(a) Improvement when fitting W +, W Z rapidity distributions and a combined fit of the three (b) Improvement when fitting Drell-Yan rapidity distributions in 4 different mass bins Figure 2: Ratio between PDF uncertainties before and after adding 1 f b 1 of LHCb data for different distributions as a function of x [5] R. D. Ball, L. Del Debbio, S. Forte, A. Guffanti, J. I. Latorre, A. Piccione, J. Rojo, M. Ubiali,"A first unbiased global NLO determination of parton distributions and their uncertainties" arxiv:hep/-ph:1002.4407; [6] S. I. Alekhin, Phys. Rev. D 63, 094022 (2001); [7] W. T. Giele, S. Keller "Implications of Hadron Collider Observables on Parton Distribution Function Uncertainties" arxiv:hep-ph/9803393v1; [8] F. De Lorenzi, Slides: http://indico.cern.ch/contributiondisplay.py?contribid=184&confid=86184 5