PDF Studies at LHCb! Simone Bifani University of Birmingham (UK) On behalf of the LHCb collaboration MENU 213 Rome 3th September - 4th October 213
Outline Introduction Analyses» Z/γ* ll, l =, e, τ#» W ν#» Z/γ* +jet Summary and Outlook Simone Bifani MENU 213 2
LHCb A Forward Spectrometer Designed to look at CP violation in B decays @ LHC Fully instrumented within 2. < η < 5.» 2.-2.5: overlap with ATLAS & CMS» 2.5-5.: complementary range 25 mrad 1 mrad Beam Beam Simone Bifani MENU 213 3
Analyses based on 21/211 data samples at 7 TeV» L 21 ~37 pb -1» L 211 ~1 fb -1 Datasets Due to luminosity levelling, same running conditions throughout fills and throughout year EW trigger thresholds unchanged between 21, 211 and 212 Simone Bifani MENU 213 4
W & Z Production and PDFs LHCb s forward acceptance provides very interesting possibilities to study the proton Parton Density Functions Take large-x from one proton and a small-x from the other probe two distinct regions in the (x, Q 2 ) space 1 9 1 8 1 7 1 6 ATLAS/CMS LHCb CDF/D W&Z HERA Fixed Target LHC 7 TeV Kinematics Can probe the low-x, high-q 2 region inaccessible to other experiments» W and Z (x of 1-4 and 1-1 )» Low-mass Drell-Yan (x down to 1-6 ) Q 2 [GeV 2 ] 1 5 1 4 1 3 1 2 1 1 W&Z γ* γ* y =6 y =4 y =2 y = y =2 y =4 y =6 1 1 7 1 6 1 5 1 4 1 3 1 2 1 1 1 x x 1,2 = M p s e ±y Q 2 = M 2 Simone Bifani MENU 213 5
W & Z Production and PDFs Theoretical predictions» Partonic cross-sections known at NNLO to few %» PDF uncertainty dominates at large rapidities (1% for y<2, 6-8% at y~5) Experimental measurements» Muon modes provide clean signature and easily reconstructible final states» Electron and tau modes provide complementary measurements and test lepton universality Cross-section measurements @ LHCb can constrain PDFs Simone Bifani MENU 213 6
Cross Section The cross-section for boson production can be expressed as BR= 1 R L N Candidates A " T rigger " Tracking " P articleid " Selection f FSR Measurements performed in the forward region (2.<η<4.5) for leptons with p T >2 GeV/c: Acceptance = 1 for Z and W ν, but obtained from MC for Z ee and Z ττ Efficiencies determined mostly from data and cross checked with simulation FSR correction evaluated using HERWIG++ and PHOTOS+PYTHIA Simone Bifani MENU 213 7
Z/γ*! LHCb-CONF-213-7 Preliminary L ~1 fb -1 (211) Single muon trigger: p T > 1 GeV/c 2 reconstructed muons» p T > 2 GeV/c» 2. < η < 4.5» 6 < M < 12 GeV/c 2 Backgrounds» Heavy flavour (Data)» K/π mis-id (Data)» Z ττ (MC) -» tt (MC)» WW (MC) N Candidates = 52626 Purity > 99% ) 2 Events / (1 GeV/c 7 6 5 4 3 2 1 LHCb preliminary Data Z/γ* + 6 8 1 12 Dimuon invariant mass [GeV/c ] 2 Simone Bifani MENU 213 8
Z/γ* ee! LHCb-PAPER-212-36 arxiv:1212.462 L ~945 pb -1 (211) Single electron trigger: p T > 15 GeV/c 2 reconstructed electrons» p T > 2 GeV/c» 2. < η < 4.5» E ECal / pc >.1» E HCal / pc <.5» M ee > 4 GeV/c 2 Backgrounds» K/π mis-id (Data)» Heavy flavour (Data)» Z ττ (MC)» tt - (MC) N Candidates = 2142 Purity ~ 95% Events / (2 GeV/c 2 ) 12 1 8 6 4 2 LHCb 4 6 8 1 12 e - ) [GeV/c 2 ] Simone Bifani MENU 213 9 M (e Data Signal Background Mass peak distorted by bremsstrahlung +
Z/γ* ττ! LHCb-PAPER-212-29 arxiv:121.6289 L ~1 fb -1 (211) Single lepton trigger:» p T,(e) > 1(15) GeV/c 5 final states:, e, e, h, eh» p T,1 > 2 GeV/c, p T,2 > 1 GeV/c» 2. < η l < 4.5, 2.25 < η h < 3.75» M 12 > 2 GeV/c 2 Backgrounds» QCD (Data)» EW (Data) -» tt (MC)» WW (MC)» Z ll (MC) Simone Bifani MENU 213 1
Z/γ* ττ! LHCb-PAPER-212-29 arxiv:121.6289 Events / (5 GeV/c 2 ) 4 Data 35 3 25 2 15 1 LHCb Ô s =7TeV Z æ QCD EWK t t æ WW æ Z æ Events / (5 GeV/c 2 ) 8 Data 7 Z æ QCD 6 EWK t t æ e 5 WW æ e 4 3 2 LHCb Ô s =7TeV Events / (5 GeV/c 2 ) 4 Data 35 Z æ QCD 3 EWK t t æ e 25 WW æ e 2 15 1 LHCb Ô s =7TeV 5 1 5 2 4 6 8 1 12 M [GeV/c 2 ] 2 4 6 8 1 12 M e [GeV/c 2 ] 2 4 6 8 1 12 M e [GeV/c 2 ] N Candidates = 99 Events / (5 GeV/c 2 ) 35 3 25 2 15 1 5 Data Z æ QCD EWK t t æ h WW æ h Z æ LHCb Ô s =7TeV Events / (5 GeV/c 2 ) 3 25 2 15 1 5 Data Z æ QCD EWK t t æ eh WW æ eh Z æ ee LHCb Ô s =7TeV Purity ~6-7% 2 4 6 8 1 12 M h [GeV/c 2 ] 2 4 6 8 1 12 M eh [GeV/c 2 ] Simone Bifani MENU 213 11
Results Z/γ* ll LHCb-CONF-213-7 Preliminary Fiducial volume: p T,l > 2 GeV/c, 2. < η l < 4.5, 6 < M ll < 12 GeV/c 2 LHCb Preliminary Z LHCb preliminary (75.4 ±.3 stat ± 1.9 syst ± 2.6 lumi ) pb Z ee (76. ±.8 stat ± 2. syst ± 2.6 lumi ) pb Z ττ NNLO(MSTW8) (71.4 ± 3.5 stat ± 2.8 syst ± 2.5 lumi ) pb (74.7 +1.6 1.4 PDF ±.4 T heory) pb 5 6 7 8 9 1 σ(z ll) [pb] Results agree with NNLO Measurements limited by statistics (ττ) or luminosity (ee, ) Main systematics: reconstruction (), tracking efficiency (ee), backgrounds and efficiency (ττ) Simone Bifani MENU 213 12
Results Z/γ* ll LHCb-CONF-213-7 Preliminary dσ/dy [pb] 7 LHCb preliminary Data (stat) 6 Data Data (tot) (tot) MSTW8 NNLO - MSTW8 (NNLO) 5 ResBos CT1 Powheg (NNLO) + PS 4 NNPDF23 (NNLO) ABM11 (NNLO) 3 2 1 2 3 4 Z boson rapidity dσ/dy [pb] 7 LHCb preliminary 6 5 4 3 2 1 Graph Data (stat) Data (tot) Z ee 2 2.5 3 3.5 4 Z boson rapidity Good agreement with NNLO Good agreement between and ee final states Good agreement with ATLAS (after adapting to ATLAS cuts) dσ/dy [pb] 14 12 1 8 6 4 2 LHCb 211 preliminary, Z extrapolated LHCb 211, Z ee extrapolated ATLAS 21, Z, Z ee.5 1 1.5 2 2.5 3 3.5 4 4.5 y(z) Simone Bifani MENU 213 13
Results Z/γ* ll LHCb-CONF-213-7 Preliminary dσ/dφ* [pb] 3 1 LHCb preliminary Data (stat) Data (tot) NNLO - MSTW8 ResBos Powheg + PS 2 1-2 1-1 1 1 φ* prediction/data 1.7 1.6 1.5 1.4 1.3 1.2 1.1 1.9.8.7 LHCb preliminary Data (stat) Data (tot) NNLO - MSTW8 ResBos Powheg + PS -2 1-1 1 1 φ* = tan (( )/2) cosh ( /2) p T m Z NNLO fails to describe φ* Bett er agreement when compared to ResBos and POWHEG dσ/dφ* [pb] 3 1 LHCb preliminary 1 2 1 1-2 1 Graph -1 1 1 Data (stat) Data (tot) Z ee φ* Simone Bifani MENU 213 14
W L ~37 pb-1 (21) LHCb-PAPER-212-8 arxiv:124.162 ν! Single muon trigger: pt > 1 GeV/c ν# 1 reconstructed & isolated muon» pt > 2 GeV/c» 2. < η < 4.5» Cone pt (R=.5) < 2 GeV/c W+ + (charged & neutral information) Backgrounds» Z/γ*» W τν and Z ττ» Heavy flavour» K/π punchthrough» K/π decay in flight Simone Bifani R (MC) (MC) (Data) (Data) (Data) R= MENU 213 p 2 + 2 15
W ν! LHCb-PAPER-212-8 arxiv:124.162 Specific cuts implemented to reduce each background component Z/γ*» Veto on 2 nd muon with p T > 2 GeV/c W τν, Z ττ and Heavy flavour» Impact parameter < 4 m K/π punchthrough» E / pc <.4 Largest residual backgrounds» Decays in flight of K/π» Z events with one muon outside the acceptance.1.2.3.4.5 IP [mm] Simone Bifani MENU 213 16 Event probability Event probability -1 1 1 1-2 -3 1-4 1 1 1-1 1-2 -3 1-4 1 LHCb LHCb Pseudo-W (Data) W ν (Simulation) bb X (Simulation) cc X (Simulation) Pseudo-W (Data) W ν (Simulation) Hadrons (Data).5 1 1.5 E / pc
W ν LHCb-PAPER-212-8 arxiv:124.162 Fit positive and negative muon p T spectra in data to expected shapes for signal and backgrounds in 5 η bins Events / (1 GeV/c) 8 7 6 5 4 3 2 1 LHCb + (a) 2 3 4 5 6 7 2 3 4 5 6 7 2. < η < 4.5 p [GeV/c] Data K/π ν Fit Z W ν T - W τν & Z ττ Heavy Flavour Signal and decay in flight templates free to vary; remaining backgrounds fixed N + Candidates = 1466 N - Candidates = 11618 Purity + ~ 79% Purity - ~ 78% W and Z muon p T spectrum at NNLO Simone Bifani MENU 213 17
Results W ν LHCb-PAPER-212-8 arxiv:124.162 Fiducial volume: p T,l > 2 GeV/c, 2. < η l < 4.5 (831 ± 9 stat ± 27 syst ± 29 lumi ) pb (656 ± 8 stat ± 19 syst ± 23 lumi ) pb 1.27 ±.2 stat ±.1 syst Results in general agreement with NNLO predictions Main systematics: background templates, efficiencies (statistical in nature) Simone Bifani MENU 213 18
l Results W ν LHCb-PAPER-212-8 LHCb-CONF-213-5 [pb] dσ W ν/dη 8 LHCb, s = 7 TeV 7 6 5 4 3 2 1 p T > 2 GeV/c + W + ν stat + W + ν tot - - W ν stat - - W ν tot 2 2.5 3 3.5 4 4.5 Data Data Data Data MSTW8 ABKM9 JR9 NNPDF21 HERA15 CTEQ6M (NLO) η [pb] dσ/dη 7 6 5 4 3 2 1 + LHCb 21, W ν - LHCb 21, W ν + ATLAS 21, W ν - ATLAS 21, W ν extrapolated extrapolated.5 1 1.5 2 2.5 3 3.5 4 4.5 η l Good agreement with NNLO Good agreement with ATLAS (after adapting to ATLAS cuts) Simone Bifani MENU 213 19
Results W ν LHCb-PAPER-212-8 LHCb-CONF-213-5 A.6.4.2 -.2 -.4 LHCb, Datastat Datatot s = 7 TeV MSTW8 ABKM9 JR9 NNPDF21 HERA15 CTEQ6M (NLO) -.6 1 2 3 4 p T > 2 GeV/c η Lepton charge asymmetry.4.2 -.2 -.4 -.6 l p T LHCb 21, W ν CMS 21, W ν CMS 21, W eν > 25 GeV/c.5 1 1.5 2 2.5 3 3.5 4 4.5 l eta Good agreement with NNLO Good agreement with CMS (after adapting to CMS cuts) Simone Bifani MENU 213 2
Z/γ* +jet! LHCb-CONF-212-16 Preliminary L ~1 fb -1 (211) Same trigger & selection as Z Jet reconstruction» Particle Flow: tracks, calorimeter clusters, Vs» Anti-k t algorithm (R=.5)» Jet ID: veto on isolated particles, require tracks from primary vertex» Jet Energy Correction from MC Jet requirements» p T > 1 GeV/c» 2. < η < 4.5» ΔR(jet, ) >.4 Event Rate [GeV/c] T p 8 7 6 5 4 3 2 1.12.1.8.6.4.2 LHCb Preliminary s = 7 TeV Data 3 LHCb preliminary s = 7 TeV 2 Inclusive Z Inclusive Z + Jet 6 8 1 12 Reconstructed Z Decay Muons Jet 1 φ -1-2 -3 M() [GeV] 2 3 η 4 5 Simone Bifani MENU 213 21
Results Z/γ* +jet! LHCb-CONF-212-16 Preliminary Number of Candidates 3 1 2 1 LHCb preliminary s = 7 TeV Data Simulation Number of Candidates 16 14 12 1 8 6 4 LHCb preliminary s = 7 TeV Data Simulation 1 2 4 6 8 1 p (Jet) [GeV] T 2 2 2.5 3 3.5 4 4.5 η(jet) Simulation well describes the data JEC is validate on data 1 st step for jet measurements at LHCb Event Rate.8.7.6.5.4.3.2.1 LHCb Preliminary s = 7 TeV Data.5 1 1.5 2 p 211 Data Simulation T (Jet) / p (Z) T Simone Bifani MENU 213 22
Results Z/γ* +jet! LHCb-CONF-212-16 Preliminary + n Jets) σ(z ) σ (Z 1-1 1 LHCb preliminary s = 7 TeV ] -1 ) [GeV / σ(z.7.6.5.4 LHCb preliminary s = 7 TeV -2 1 dσ (Z +Jet) d pt (Z ).3.2.1 1 2 >2 Jet Multiplicity 2 4 6 8 1 p (Z T o ) [GeV] Z!+jet Z! Data :.229 ±.6 stat ±.9 syst Theory :.219 +.6.9 PDF ±.12 scale Simone Bifani MENU 213 23
Summary and Outlook LHCb probes a unique region in η and low p T reach at LHC Z and W production cross-sections measured at 7 TeV in the forward region by using different final states Results consistent with theoretical predictions calculated at NNLO using recent parameterisation of the PDFs Upcoming updates with full 211 and 212 data sets for all analyses Z/γ* +jet paper coming out soon Low mass Drell-Yan: LHCb-CONF-212-12 Simone Bifani MENU 213 24
Backup
LHCb A Forward Spectrometer Complementary η range to ATLAS & CMS» Overlap for cross check 2. < η < 2.5» Unique to LHCb 2.5 < η < 5. η# in LHCb 17%(16%) of W + (W - ) within LHCb 8% of Z within LHCb Simone Bifani MENU 213 26
W & Z Production and PDFs Cancel or highlight PDF uncertainties with ratios» A = (dσ W+ - dσ W- ) / (dσ W+ + dσ W- ) tests u V and d V difference A # R W#» R W = dσ W+ / dσ W- tests d V /u V ratio» R WZ = dσ W / dσ Z almost insensitive to PDFs precise test of SM Many systematic errors cancel Simone Bifani MENU 213 27
Reconstruction Efficiencies Trigger» Tag: triggered muon#» Probe: offline identified muon Tracking Tracking» Tag: identified muon track» Probe: trajectory from muon stub and minimal tracking information Particle ID Particle ID» Tag: identified muon#» Probe: reconstructed track Simone Bifani MENU 213 28
Z/γ* ll Simone Bifani MENU 213 29
Z/γ*! Simone Bifani MENU 213 3
Z/γ* ee! Simone Bifani MENU 213 31
Z/γ* ττ! Simone Bifani MENU 213 32
Z/γ* ττ! Simone Bifani MENU 213 33
W ν Event probability -1 1 1 1-2 LHCb Pseudo-W (Data) W ν (Simulation) IP > 8 m (Data) Event probability 1 1-1 1-2 LHCb Pseudo-W (Data) W ν (Simulation) IP > 8 m (Data) -3 1-3 1-4 1 1 2 3 [GeV/c] pcone T -4 1 1 2 3 [GeV] cone E T Simone Bifani MENU 213 34
W ν Fit positive and negative muon p T spectra in data to expected shapes for signal and backgrounds in 5 η bins Events / (1 GeV/c) 25 2 15 1 5 2 3 4 5 6 7 2 3 4 5 6 7 LHCb + - (b) 2 3 4 5 6 7 2 3 4 5 6 7 2. < η < 2.5 p [GeV/c] T Events / (1 GeV/c) 2 3 4 5 6 7 2 3 4 5 6 7 3 LHCb + - 25 (c) 2 15 1 5 2 3 4 5 6 7 2 3 4 5 6 7 2.5 < η < 3. p [GeV/c] T 22 2 LHCb + - 18 (d) 16 14 12 1 8 6 4 2 2 3 4 5 6 7 2 3 4 5 6 7 3. < η < 3.5 p [GeV/c] Events / (1 GeV/c) T Events / (1 GeV/c) 1 8 6 LHCb + (e) - Events / (1 GeV/c) 25 2 15 LHCb + (f) - 4 1 2 5 2 3 4 5 6 7 2 3 4 5 6 7 3.5 < η < 4. p [GeV/c] T 2 3 4 5 6 7 2 3 4 5 6 7 4. < η < 4.5 p [GeV/c] W and Z muon p T spectrum at NNLO Simone Bifani MENU 213 35 T
W ν Simone Bifani MENU 213 36
Z/γ* +jet! Simone Bifani MENU 213 37
Z/γ* +jet! Jet Efficiency 1.8.6 LHCb simulation s = 7 TeV Jet Efficiency 1.8.6 LHCb simulation s = 7 TeV.4.2 p p jet T jet T > 5 GeV > 15 GeV.4.2 1.5<η <5. jet 2.<η <4.5 jet 2 3 4 5 jet η 1 2 3 4 jet p T (GeV) Simone Bifani MENU 213 38
References Simone Bifani MENU 213 39
Low Mass DY! LHCb-CONF-212-12 Preliminary L ~37 pb -1 (21) Dimuon trigger: p T > 2.5 GeV/c 2 reconstructed muons» p > 1 GeV/c» p T > 3(15) GeV/c (M > 4 GeV/c 2 )» 2. < η < 4.5» 5 < M < 12 GeV/c 2 Backgrounds» Heavy flavour (Data)» K/π mis-id (Data)» Υ (MC) Fit muon isolation spectrum in data to expected shapes for signal and backgrounds in 9 M bins isolation = p T /p-jet T Purity varies from ~7% (at low mass) to ~1% Simone Bifani MENU 213 4
Results Low Mass DY LHCb-CONF-212-12 Preliminary )] 2 [pb/(gev/c dσ γ* /dm Theory/Data 2 1 LHCb Preliminary, s = 7 TeV 1 1 1.5 1.5 1 p > 1 GeV/c p > 3 (15) GeV/c T 2. < η < 4.5 Data stat Data tot LO PYTHIA (CTEQ5L) NLO FEWZ (MSTW8) NLO DYNNLO (MSTW8) 2 2 ] 1 Dimuon invariant mass [GeV/c [pb] γ* dσ γ* /dy 4 35 3 25 2 15 1 5 p p T > 1 GeV/c > 3 GeV/c 1.5 < M < 2 GeV/c 2 LHCb Preliminary, Data stat Data tot s = 7 TeV LO PYTHIA (CTEQ5L) NLO FEWZ (MSTW8) NLO FEWZ (NNPDF2) NLO FEWZ (CTEQ66) 2 2.5 3 3.5 4 4.5 Dimuon rapidity Main systematic at low mass comes from the fit PYTHIA predictions agree in shape but too low in normalisation Reasonable agreement with NLO predictions where calculations are available Simone Bifani MENU 213 41
Low Mass DY! Events Events Events 22 LHCb preliminary 2 21 data 18 16 14 12 1 8 6 4 2 5 4 3 2 1 12 1 2 5. < M < 7. GeV/c Drell-Yan Upsilon Heavy Flavour Hadron mis-id 2. < y < 4.5.5 1 Minimum muon isolation.5 1 Minimum muon isolation 8 6 4 2 2 12.5 < M < 15. GeV/c 2 3. < M < 4. GeV/c 2. < y < 4.5 2. < y < 4.5.5 1 Minimum muon isolation Events Events Events 2 7. < M < 9. GeV/c 2. < y < 4.5 24 22 2 18 16 14 12 1 8 6 4 2.5 1 Minimum muon isolation 2 15. < M < 2. GeV/c 2. < y < 4.5 45 4 35 3 25 2 15 1 5.5 1 Minimum muon isolation 8 7 6 5 4 3 2 1 2 4. < M < 6. GeV/c 2. < y < 4.5.5 1 Minimum muon isolation Events Events Events 8 7 6 5 4 3 2 1 3 25 2 15 1 2 1.5 < M < 12.5 GeV/c.5 1 Minimum muon isolation 5 2 2. < M < 3. GeV/c 2 6. < M < 12 GeV/c 2. < y < 4.5 2. < y < 4.5.5 1 Minimum muon isolation 2. < y < 4.5 9 8 7 6 5 4 3 2 1.5 1 Minimum muon isolation Simone Bifani MENU 213 42
Low Mass DY! Simone Bifani MENU 213 43