Neutrino data fit (version 5, modified by Arie)

• Use A, B to get Xw variable, where A and B come from the fit version 4.

• 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.

• F2val-vector(Nucleon) = 0.5*(1-G2^2)*[uv*(Q2+C2vv(u))/(Q2+C1vv(u)) + dv*(Q2+C2vv(d))/(Q2+C1vv(d)) ],
  where C1vv(u), C1vv(d), C2vv(u), C2vv(d) are given from the fit version 4.

• F2val-axial(Nucleon) = 0.5*(1-G2^2)*[uv*(Q2+C2av(u))/(Q2+C1av(u)) + dv*(Q2+C2av(d))/(Q2+C1av(d)) ],
  where C1av(u)=C1vv(d)=0.25, C2av(u)=Z*0.25, C2av(d)=(2Z+0.605)*0.25. Z will be fit.

• F2sea-vector(Nucleon) = 0.5*[ (2*us)*Q2/(Q2+C1vs(u)) + (2*ds+2*str)*Q2/(Q2+C1vs(d)) ],
  where C1as(u) and C1as(d) are given from the fit version 4.

• F2sea-axial(Nucleon) = 0.5*[ (2*us)*(Q2+C2as(u))/(Q2+C1as(u)) + (2*ds+2*str)*(Q2+C2as(d))/(Q2+C1as(d)) ],
  where C1as(u), C1as(d), C2as(u), C2as(d) will be 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.

Fit Result and comparison

Without the lowest x bin