# Neutrino data fit (version 5.1 , modified by Arie)

Used fit procedure (Please correct me!). These are only correct if F2 is corrected for Charm Mass Effects and
gluon diagrams and C-Cbar sea.  So all F2 tables must have correction factors applied.

• 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*                 [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(Nucleon) = [ F2val(Nucleon) + F2sea(Nucleon) ] * fPDF
• F2val-vector(Proton) = 0.5*(1-G2^2)*[                                    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(Proton) = 0.5*                 [                                                      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(Proton) = 0.5*[                              (2*ds+2*str)*Q2/(Q2+C1vs(d)) ],
where C1as(u) and C1as(d) are given from the fit version 4.

• F2sea-axial(Proton) = 0.5*[                                (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(Proton) = [ F2val(Proton) + F2sea(Pproton) ] * fPDF
• F2val-vector(Neutron) = 0.5*(1-G2^2)*[uv*(Q2+C2vv(u))/(Q2+C1vv(u)) +                                 ],
where C1vv(u), C1vv(d), C2vv(u), C2vv(d) are given from the fit version 4.

• F2val-axial(Neutron) = 0.5*                 [uv*(Q2+C2av(u))/(Q2+C1av(u)) +                                  ],
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(Neutron) = 0.5*[ (2*us)*Q2/(Q2+C1vs(u)) +                                                         ],
where C1as(u) and C1as(d) are given from the fit version 4.

• F2sea-axial(Neutron ) = 0.5*[ (2*us)*(Q2+C2as(u))/(Q2+C1as(u)) +                                           ],
where C1as(u), C1as(d), C2as(u), C2as(d) will be fit.

• Get F2(LO) from F2(Neutron) = [ F2val(Neutron) + F2sea(Neutron) ] * 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.

## Without the lowest x bin

20040122-Inkyu.Park@bnl.gov