Used fit procedure (Please correct me!)
- Use A, B to get Xw variable, where A and B come from the fit.
- Call grv98_lo_du(Xw,Fs,uv,dv,us,ds,str,chm,btm,top,glu), and get parton distribution values of each quarks, uv, dv, us, ds, str, chm, btm, top, glu at a given (Xw,Q2) point.
- No ccbar bbar productions are considered.
- Calculate valence quark contribution as follow:
F2val(P) = (1-G2^2) * [ uv*(Q2+Cv2(u))/(Q2+Cv1(u)) + dv*(Q2+Cv2(d))/(Q2+Cv1(d)) ],
where Cv1(u), Cv1(d), C2(u), C2(d) come from the fit.
- Calculate sea quark contribution as follow:
F2sea(P)= [ (2*us)*Q2/(Q2+Cs1(u)) + (2*ds+2*str)*Q2/(Q2+Cs1(d)) ],
where Cs1(u) and Cs1(d) come from the fit.
- Get F2(LO) from F2(P) = [ F2val(P) + F2sea(P) ] * fPDF
- For Deutreum target data (WA25), we multiply by
0.9902 + 0.3233*Xb - 2.164*Xb^2 + 6.070*Xb^3 - 8.539*Xb^4 + 4.613*Xb^5
- For Neon target data (WA59), we further multiply by
1.0963 - 0.36427*Xw - 0.27805*exp(-21.936*Xw)+2.7715*Xw^14.417
- For CCFR data, F2(Heavy)/F2(light) correction table was used to scale the CCFR data as Un-Ki directed. I.e., F2(CCFR)=F2(CCFR)/Correction
- red solid line is with ANUCL(Xb) correction (F2(iron)/F2(Deutereum), while black dotted line is with only Deutereum correction.
WA25
Reference: Allasia et al. Zeit.Phys.C28(91)321
o F2(nucleon) from nu-deuterium interactions
F2(nucleon) from Nu-Deutereum scattering
WA59
Reference: Varvell et al. Zeit.Phys.C36(91)1
o F2(nucleon) and xF3(nucleon) from neutrino-neon scattering
F2(nucleon) from Nu-Neon scattering
CCFR
Reference: Un-Ki's thesis
o F2(nucleon) from neutrino-iron scattering
F2(nucleon) from Nu-Iron scattering
