1 / 34

Low-x at the LHC and with UHECR’s Low-x 2006 Workshop, Lisbon A.M Cooper-Sarkar (Oxford)

Low-x at the LHC and with UHECR’s Low-x 2006 Workshop, Lisbon A.M Cooper-Sarkar (Oxford). The LHC is a low-x machine- at least for the first few years What can W/Z rapidity spectra tell us about the low-x gluon PDF? What can the W decay lepton asymmetry tell us about low-x valence PDFs?

zizi
Download Presentation

Low-x at the LHC and with UHECR’s Low-x 2006 Workshop, Lisbon A.M Cooper-Sarkar (Oxford)

An Image/Link below is provided (as is) to download presentation Download Policy: Content on the Website is provided to you AS IS for your information and personal use and may not be sold / licensed / shared on other websites without getting consent from its author. Content is provided to you AS IS for your information and personal use only. Download presentation by click this link. While downloading, if for some reason you are not able to download a presentation, the publisher may have deleted the file from their server. During download, if you can't get a presentation, the file might be deleted by the publisher.

E N D

Presentation Transcript


  1. Low-x at the LHC and with UHECR’sLow-x 2006 Workshop, Lisbon A.M Cooper-Sarkar (Oxford) • The LHC is a low-x machine- at least for the first few years • What can W/Z rapidity spectra tell us about the low-x gluon PDF? • What can the W decay lepton asymmetry tell us about low-x valence PDFs? • Can we use the conventional DGLAP NLO QCD formalism? • Future information from Ultra-High Energy Cosmic Rays

  2. LHC is a low-x machine (at least for the early years of running) Low-x information comes from evolving the HERA data W/Z production have been considered as good standard candle processes possibly even monitors of the luminosity But how do the uncertainties on PDFs affect these cross-sections?

  3. fa pA x1 • where X=W, Z, D-Y, H, high-ET jets, prompt-γ • and  is known • to some fixed order in pQCD and EW • in some leading logarithm approximation (LL, NLL, …) to all orders via resummation ^ x2 pB fb X Transporting PDFs to hadron-hadron cross-sections QCD factorization theorem for short-distance inclusive processes Need PDFs to calculate any cross-section Need estimates of the uncertainties on PDFs to estimate the uncertainty on any cross-section

  4. What has HERA data ever done for us? W+ W- Z Pre-HERA W+/W-/Z and W→ lepton rapidity spectra ~ ± 15% uncertainties ! Why? Because the central rapidity range AT LHC is at low-x (5 10-4 to 5 10-2) NO WAY to use these cross-sections as a good luminosity monitor Lepton+ Lepton-

  5. W+ W- Z Lepton+ Lepton- Post-HERA W+/W-/Z and W → lepton rapidity spectra ~ ± 5% uncertainties ! Why? Because there has been a tremendous improvement in our knowledge of the low-x glue and thus of the low-x sea

  6. Pre-HERA sea and glue distributions Post HERA sea and glue distributions

  7. 9 The uncertainty on the W/ Z rapidity distributions is dominated by –- gluon PDF dominated eigenvectors We can decompose the uncertainties on the PDF parameters into eigenvector combinations -the largest uncertainty is along eigenvector 9 –which is dominated by the low- x gluon uncertainty . Differences are visible within the measurable range of rapidity It may at first sight be surprising that W/Z distns are sensitive to gluon parameters BUT our experience is based on the Tevatron where Drell-Yan processes can involve valence-valence parton interactions. At the LHC we will have dominantly sea-sea parton interactions at low-x And at Q2~MZ2 the sea is driven by the gluon- which is far less precisely determined for all x values

  8. So – let’s put some experimental realism into the predictions for W/ZW/Z production have been considered as good standard candle processes insensitive to PDF uncertainties……? This is true WITHIN a PDFset generator level But how about comparing PDFsets? We actually measure the decay lepton spectra Generate with HERWIG+k-factors (checked against MC@NLO) using CTEQ6.1MZEUS_S MRST2001 PDFs with full uncertainties from LHAPDF eigenvectors At y=0 the total uncertainty is ~ ±6% from ZEUS ~ ±4% from MRST01E ~ ±8% from CTEQ6.1 ZEUS to MRST01 central value difference ~5% Could we improve the situation with ALTAS data? We’d NEED to control systematic errors to within ~4% electron positron ATLFAST electron positron

  9. Study the effect of including the LHC W Rapidity distributions in global PDF fits by how much can we reduce the PDF errors with early LHC data? Generate data with 4% error using CTEQ6.1 PDF, pass through ATLFAST detector simulation and then include this pseudo-data in the global ZEUS PDF fit Central value of prediction shifts and uncertainty is reduced BEFORE including W data AFTER including W data AMCS, A. Tricoli (Hep-ex/0509002) Improvement in PDF errors of ~30% Lepton+ rapidity spectrum data generated with CTEQ6.1 PDF compared to predictions from ZEUS PDF AFTER these data are included in the fit Lepton+ rapidity spectrum data generated with CTEQ6.1 PDF compared to predictions from ZEUS PDF Specifically the low-x gluon shape parameter λ, xg(x) = x –λ , was λ = -.199 ± .046 for the ZEUS PDF before including this pseudo-data It becomes λ = -.181 ± .030 after including the pseudodata

  10. The uncertainty on the W/Z rapidity distributions is dominated by –- gluon PDF dominated eigenvectors andthere is cancellation in the ratios AW = (W+ - W-)/(W+ + W-) ZW = Z/(W+ + W-) Remaining uncertainty comes from valence PDF related eigenvectors Well Known?Gold plated? We will measure the lepton asymmetry generator level Within each PDF set uncertainty in the lepton asymmetry IS LESS than in the lepton rapidity spectra, e.g about 2% for the asymmetry at y=0, as opposed to about 4% for the lepton rapidity spectra themselves (using MRST2001 PDFS) However the PDF sets differ from each other more strikingly- MRST01and CTEQ6.1 differ by about 13% at y=0! But this is an opportunity to use ATLAS measurements to increase knowledge of the valence PDFs at x~0.005 ATLFAST

  11. Relationship of valence PDFs to Asymmetry measurement • Dominantly, at LO Aw= (u d – d u) • (u d + d u) • And u = d = q at small x ( I will return to this..it can be challenged) • So Aw~ (u – d) = (uv – dv) • (u + d) (uv + dv + 2 q ) • Actually this pretty good even quantitatively – Make this simple LO calculation using the CTEQ and the MRST PDFs at Q2=MW2 and x~0.006 (corresponding to zero rapidity) and you’ll find that it DOES explain the Aw difference between MRST and CTEQ. - let’s see the evidence

  12. CTEQ6.1 uv MRST02 Q2=Mw2 uv – dv Q2=Mw2 dv x- range affecting W asymmetry in the measurable rapidity range MRST and CTEQuv – dv distributions are significantly different : at Q2=MW2 and x~0.006 CTEQ: uv=0.17, dv=0.11, uv-dv=0.06 MRST: uv=0.155, dv=0.11,uv-dv=0.045 Q2=Mw2

  13. CTEQ6.1 dbar MRST02 dbar-ubar ubar Evidence that dbar = ubar for both PDFs at small-x at Q2=MW2 and x~0.006 MRST and CTEQ dbar=ubar=0.7 So Aw=0.06/(0.28+1.4) = 0.036 CTEQ Aw=0.045/(0.265+1.4) = 0.027 MRST…pretty close!

  14. Generate data with 4% error using CTEQ6.1 PDF, compare to ZEUS-PDF predictions and then include this pseudo-data in the global ZEUS PDF fit - shape remains very similar and uncertainty is reduced from ~8% to ~5% at y=0 (this improvement is spread amongst the valence parameters, not in any one parameter) BEFORE including AW data AFTER including AW data Cteq6.1 pseudodata ZEUS-S prediction CTEQ6.1 pseudodata input to the ZEUS fit MRST02 pseudodata ZEUS-S prediction

  15. Generate data with 4% error using MRST02 PDF, compare to ZEUS-PDF predictions and then include this pseudo-data in the global ZEUS PDF fit - cannot fit with standard ZEUS-Global fit parametrisation – need to extend the valence parametrisations The valence parametrisations at Q20 become xuv = Au xau (1-x)bu(1+du √x + cu x) , xdv = Ad xad (1-x)bd (1+dd√x + cd x)- 16 rather than 14 parameters are needed to fit MRST pseudodata MRST02 pseudodata ZEUS-S prediction MRST02 pseudodata input to the ZEUS fit with 16 parameters Clearly there will be low-x valence information from the ATLAS lepton asymmetry data

  16. But all of this was within the standard DGLAP NLOQCD theoretical paradigm Is NLO (or even NNLO) DGLAP good enough? The QCD formalism may need extending at small-x

  17. MRST have produced a set of PDFs derived from a fit without low-x data –ie do not use the DGLAP formalism at low-x- called MRST03 ‘conservative partons’. These give VERY different predictions for W/Z production to those of the ‘standard’ PDFs. Z W+ W- MRST02 W+ Z W- MRST03

  18. Differences persist in the decay lepton spectra and even in theirratio and asymmetry distributions MRST02 MRST02 MRST02 MRST02 MRST03 MRST03 MRST03 MRST03 Reconstructed Electron Pseudo-Rapidity Distributions (ATLAS fast simulation) 200k events of W+- -> e+-generated with HERWIG 6.505 + NLO K factors Reconstructed e- Reconstructed e+ 6 hours running h h Reconstructed e+- e- Asymmetry Reconstructed e- / e+ Ratio h

  19. Note of caution. MRST03 conservative partons DO NOT describe the HERA data for x< 5 10-3 which is not included in the fit which produces them. So there is no reason why they should correctly predict LHC data at non-central y, which probe such low x regions. What is really required is an alternative theoretical treatment of low-x evolution which would describe HERA data at low-x, and could then predict LHC W/Z rapidity distributions reliably – also has consequences for pt distributions. The point of the MRST03 partons is to illustrate that this prediction COULD be very different from the current ‘standard’ PDF predictions. When older standard predictions for HERA data were made in the early 90’s they did not predict the striking rise of HERA data at low-x. This is a warning against believing that a current theoretical paradigm for the behaviour of QCD at low-x can be extrapolated across decades in Q2 with full confidence. → The LHC measurements may also tell us something new about QCD

  20. UHECR section The Pierre Auger Array is sensitive to neutrinos of energy > 108 GeV Eν ~108 -1010 GeV → x ~10-4 – 10-6 at Q2 ~104 GeV2 For HERA such x values have low Q2 And LHC cannot access such low-x except at high rapidity- a hadron collider of cm energy ~103 TeV would be necessary for central rapidity and Q2~104 So the UHE νN cross-section is dependent on low-x physics We (hep-0605086) have evaluated the expected energy dependence of the cross-section using the ZEUS global PDF fit analysis within the framework of DGLAP NLO QCD using a general mass heavy quark treatment andincluding all HERA-I data and accountig for the PDF uncertainties σNLO = 6.04 ±0.40 (Eν/GeV) 0.358±0.005 pb • But we should extend the theoretical formalism – two possibilities have been used • BFKL/DGLAP + screening –Kutak Kwiecinski • CGC-Henley-JalilianMarian

  21. The problem with νN cross-section measurements is that there are two different benchmark cosmic fluxes- the Cosmogenic and the Waxman-Bahcall Butνμ and ντare maximally mixed so cosmic ν’s will contain a substantial component of ντ and earth skimming (ES) ντ’s give rise to upward going showers By contrast ν’s of all flavours give deeply penetrating quasi-horizontal (QH) showers with a distinctive em component. The ratio of QH/ES events gives a robust estimate of the νN cross-section- almost independent of the flux The figure shows the numbers of QH and ES events for pseudo-measurements made with a Super-Auger array of 10000 sqm over 10 years! Such measurements could distinguish theoretical predictions for low-x –particularly if the true flux is closer to the WB flux (Patience is a virtue!)

  22. Summary Measuring W/Z rapidity spectra at the LHC can reduce the uncertainty on the low-x gluon PDF Measuring the W decay lepton asymmetry can give us information on the valence PDFs at x~0.005 on a proton target The LHC may show us that we need to extend the conventional DGLAP NLO QCD formalism In future the Pierre Auger array may measure Ultra-High Energy neutrino cross-sections and be able to distinguish alternative theoretical formalisms at low-x

  23. Extra after ere

  24. But perhaps we need to look more closely at the tiny difference in dbar-ubar– look at MRST and CTEQ dbar-ubar on a scale blown up by x10 Could this play a role? Without approximations..but still LO, Aw = d (uv-dv) +dvΔ d (uv+dv+2d) -2dΔ –dvΔ Where Δ=dbar-ubar=0.016 for CTEQ Aw = 0.7x(0.06) + 0.11x0.016 0.7X(1.68) – 2x0.7x0.016-0.11x0.016 And Δ=0.010 for MRST AW=0.7x(0.045) + 0.11x0.011 0.7x(1.665) – 2x0.7x0.011-0.11x0.011 The terms involving the difference of ubar and dbar are simply too small to matter compared to the terms involving the valence difference.

  25. Now look at PDF uncertainties more generally- this section mostly about high-x PDF uncertainties at the LHC! The general trend of PDF uncertainties is that The u quark is much better known than the d quark The valence quarks are much better known than the sea and the gluon at high-x The valence quarks are poorly known at small-x - but they are not important for physics in this region since the sea quarks are dominant The sea and the gluon are well known at low-x The sea is poorly known at high-x, but the valence quarks are more important in this region The gluon is poorly known at high-x And it can still be very important for physics e.g.– high ET jet xsecn need to tie down the high-x gluon

  26. Example of how PDF uncertainties matter– Tevatron jet data were originally taken as evidence for new physics-- i These figures show inclusive jet cross-sections compared to predictions based on theory in the form (data - theory)/ theory Something seemed to be going on at the highest E_T And special PDFs like CTEQ4/5HJ were tuned to describe it better- note the quality of the fits to the rest of the data deteriorated. But this was before uncertainties on the PDFs were seriously considered

  27. Today Tevatron jet data are considered to lie within PDF uncertainties. (Example from CTEQ hep-ph/0303013) We can decompose the uncertainties into eigenvector combinations of the fit parameters-the largest uncertainty is along eigenvector 15 –which is dominated by the high x gluon uncertainty

  28. And we can translate the current level of PDF uncertainty into the uncertainty on LHC jet cross-sections. This will impact on any BSM physics signal which can be expressed as a contact term

  29. SM 2XD 4XD 6XD E.G. Dijet cross section potential sensitivity to compactification scale of extra dimensions (Mc) reduced from ~6 TeV to 2 TeV. (Ferrag et al) Mc = 2 TeV, no PDF error Mc = 6 TeV, no PDF error Mc = 2 TeV, with PDF error

  30. g q H W/Z t H W/Z g q W/Z And how do PDF uncertainties affect the Higgs discovery potential? S Ferrag

  31. Inputting data to a PDF fit needs a prediction for the cross-section which can be easily obtained analytically –true for DIS inclusive cross-section. But many NLO cross-sections can only be computed by MC and can take 1-2 CPU days to compute. This cannot be done for every iteration of a PDF fit. Recently grid techniques have been developed to include DIS jet cross-sections in PDF fits (ZEUS-JETs fit)

  32. This technique has been extended to LHC high-ET jet cross-sections Use data at lower PT and higher η-where new physics is not expected T. Carli Dan Clements, Glasgow Claire Gwenlan, Oxford

  33. Such PDF uncertainties on the jet cross sections compromise the potential for discovery. E.G. Dijet cross section potential sensitivity to compactification scale of extra dimensions (Mc) reduced from ~6 TeV to 2 TeV. (Ferrag et al) Mc = 6 TeV, no PDF error The reduced uncertainties can then be used in background calculations for new physics signals

  34. Prompt photon cross-sections may also be amenable to this treatment –compare photon pt and η distributions for Cteq61 up and down eigenvector 15 -emphasizing the uncertainties in the high-x gluon (PT > 330 GeV)Ivan Hollins Birmingham

More Related