REFINEMENT OF THE PION PDF IMPLEMENTING DRELL- YAN EXPERIMENTAL - - PowerPoint PPT Presentation

refinement of the pion pdf implementing drell yan
SMART_READER_LITE
LIVE PREVIEW

REFINEMENT OF THE PION PDF IMPLEMENTING DRELL- YAN EXPERIMENTAL - - PowerPoint PPT Presentation

REFINEMENT OF THE PION PDF IMPLEMENTING DRELL- YAN EXPERIMENTAL DATA Patrick B Barry 1 , Nobuo Sato 2 , W. Melnitchouk 3 , Chueng-Ryong Ji 1 pcbarry@ncsu.e .edu North Carolina State University 1 University of Connecticut 2 Thomas Jefferson


slide-1
SLIDE 1

REFINEMENT OF THE PION PDF IMPLEMENTING DRELL- YAN EXPERIMENTAL DATA

Patrick B Barry1, Nobuo Sato2, W. Melnitchouk3, Chueng-Ryong Ji1 pcbarry@ncsu.e .edu North Carolina State University1 University of Connecticut2 Thomas Jefferson National Accelerator Facility3

slide-2
SLIDE 2

Drell-Yan process

■ Two hadrons collide

– Do not need to be protons! ■ One donates a quark, other an antiquark ■ Quarks annihilate into a virtual photon ■ Dilepton production ■ Measure differential cross- section of lepton/antilepton pair

slide-3
SLIDE 3

Cross section

■ Cross-section differential in invariant mass of lepton pair, 𝑅", and rapidity, 𝑍 ■ The momentum fractions of the initial hadrons are 𝑦%, 𝑦" ■ Parton distribution functions (PDFs) are 𝑔

( 𝑦, 𝑅"

■ Sum over all partons

slide-4
SLIDE 4

PDFs

■ Parameterize the PDF at 𝑅)

" = 1GeV2 as:

𝑔 𝑦, 𝜈" = 𝑂 ¡𝑦/ 1 − 𝑦1 ■ Definitions: 𝑟3 = 𝑣 53 = 𝑒3; 𝑟7 = 2 ¡𝑣 + 𝑒̅ + 𝑡 ; 𝑕 ■ Use sum rules to fix 𝑂=>, 𝑂

?

■ We fit 𝑏, 𝑐 for the valence, sea, and gluon, and 𝑂 for the sea ■ PDFs are evolved using DGLAP in Mellin space

slide-5
SLIDE 5

Nested Sampling

■ Monte Carlo fitting method ■ Create parameter space to have a uniform prior over a specified range ■ Sample points in parameter space closer and closer to the maximum likelihood ■ Weights produced with each sample based on proximity to maximum likelihood ■ Provides errors without assumption of linear error propagation 𝑊𝑏𝑠 𝒫 ∝ F 𝒫G − 𝐹 𝒫

"

  • G
slide-6
SLIDE 6

Nested Sampling

■ Start with random points on line of 0 < 𝑌 < 1. ■ 𝑌 = 0 is the point

  • f highest

likelihood.

slide-7
SLIDE 7

Nested Sampling

■ Start with random points on line of 0 < 𝑌 < 1. ■ 𝑌 = 0 is the point

  • f highest

likelihood. ■ Delete point of lowest likelihood and make it the upper-bound on new sampling boundary

slide-8
SLIDE 8

Nested Sampling

■ Start with random points on line of 0 < 𝑌 < 1. ■ 𝑌 = 0 is the point

  • f highest

likelihood. ■ Delete point of lowest likelihood and make it the upper-bound on new sampling boundary

slide-9
SLIDE 9

Nested Sampling

■ Start with random points on line of 0 < 𝑌 < 1. ■ 𝑌 = 0 is the point

  • f highest

likelihood. ■ Delete point of lowest likelihood and make it the upper-bound on new sampling boundary ■ Keep randomly sampling until a threshold is reached

slide-10
SLIDE 10

Nested Sampling

■ Start with random points on line of 0 < 𝑌 < 1. ■ 𝑌 = 0 is the point

  • f highest

likelihood. ■ Delete point of lowest likelihood and make it the upper-bound on new sampling boundary ■ Keep randomly sampling until a threshold is reached

slide-11
SLIDE 11

Nested Sampling

■ Start with random points on line of 0 < 𝑌 < 1. ■ 𝑌 = 0 is the point

  • f highest

likelihood. ■ Delete point of lowest likelihood and make it the upper-bound on new sampling boundary ■ Keep randomly sampling until a threshold is reached

slide-12
SLIDE 12

Datasets & Constrictions

■ For Drell-Yan, we use E615 and NA10 datasets – 𝜌N beam incident on a Tungsten target – Consider only 0 < 𝑦O < 0.9 and 4.16 < 𝑅 < 8.34 to avoid 𝐾/Ψ and Υ production

slide-13
SLIDE 13

Drell-Yan fits

slide-14
SLIDE 14

Drell-Yan fits

  • 𝑟3 is well-constrained by Drell-

Yan

  • 𝑟7has large spread in

parameters

  • 𝑕 has almost not constrain
slide-15
SLIDE 15

Drell-Yan fits

slide-16
SLIDE 16

Leading Neutron

■ Add in data from HERA (ZEUS & H1) to perform global fit ■ Detect neutrons in coincidence with outgoing electrons: ■ Neutron has most of the energy of the proton ■ Incoming electron barely strikes the surface of the proton, knocking out a pion from the pion cloud ■ Focuses on small 𝑦Z, whereas Drell-Yan focuses on large 𝑦Z

slide-17
SLIDE 17

Leading Neutron

■ Observable in H1 data is 𝐺

" \] ^ 𝑦, 𝑅", 𝑧 = 𝑔Z`a(𝑧) ¡𝐺" Z(𝑦Z, 𝑅")

– Where 𝑔Z`a(𝑧) ¡is the splitting function from the proton, and ¡𝐺"

Z(𝑦Z, 𝑅") is the pion

structure function (depends on pion PDF) ■ Observable in ZEUS is 𝑠 𝑦Z, 𝑅", 𝑧 = 𝑔Z`a 𝑧 𝐺"

Z 𝑦Z, 𝑅"

𝐺

" d 𝑦, 𝑅"

Δ𝑧 – where 𝐺

" d 𝑦, 𝑅" is the proton structure function

slide-18
SLIDE 18

Datasets & Constrictions

■ For Drell-Yan, we use E615 and NA10 datasets – 𝜌N beam incident on a Tungsten target – Consider only 0 < 𝑦O < 0.9 and 4.16 < 𝑅 < 8.34 to avoid 𝐾/Ψ and Υ production ■ For Leading Neutron, we use H1 and ZEUS datasets – We consider cuts on data based on maximum 𝑧 = 𝑦Z/𝑦 values

slide-19
SLIDE 19

LN results

𝑧-cut of 0.2

slide-20
SLIDE 20

LN results

𝑧-cut of 0.2

slide-21
SLIDE 21

LN results

𝑧-cut of 0.2

slide-22
SLIDE 22

Conclusion

■ First attempted fit to both high-𝑦Z and low- 𝑦Z regions using Drell-Yan and Leading Neutron data ■ Use of nested sampling algorithm to improve errors ■ Next steps: to include threshold resummation in our calculation

slide-23
SLIDE 23

BACKUP SLIDES

slide-24
SLIDE 24

Prediction of E866

■ Can make a prediction of E866 data for 𝑒̅ − 𝑣 5 = (𝑔Z`a − "

^ 𝑔Zfg``) ⊗ 𝑟

53

Z using our

valence 𝜌 PDF, where 𝑔Z`a and 𝑔Zfg`` are the splitting functions from the proton

slide-25
SLIDE 25

Kinematics - DY

0.2 0.4 0.6 0.8 xπ 101 102

2 × 101 3 × 101 4 × 101 6 × 101

Q2

E615 NA10-194GeV NA10-286GeV

  • Hard cut-offs for

4.16" < 𝑅" < 8.34"

  • More available data for

large-𝑦i

slide-26
SLIDE 26

Kinematics - LN

10−3 10−2 10−1 xπ 101 102 103 Q2

H1 ZEUS

  • 𝑧-cut of 0.2
slide-27
SLIDE 27

Normalization Parameterization

■ For all datasets with overall normalization uncertainty, we fit to within the reported percentage around 1.

LN DY

slide-28
SLIDE 28

Mellin Transformation

■ Analogous to the Fourier transform ■ Transform from x-space to Mellin space (exponents of x)

WH WHY??

■ We know how PDFs evolve in scale based on DGLAP:

slide-29
SLIDE 29

Mellin Inversion

■ After evolution, invert back into x-space

  • For each value on the contour,

we do the DGLAP evolution

  • At large enough contour

radius, integrand converges