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.
Fit Result and comparison
Without the lowest x bin
20040122-Inkyu.Park@bnl.gov