We present our QCD analysis of the proton structure function $F_2^p(x,Q^2)$ to determine the parton distributions at the next-to-leading order (NLO). The heavy quark contributions to $F_2^i(x,Q^2)$, with $i$ = $c$, $b$ have been included in the frame
work of the `fixed flavour number scheme (FFNS). The results obtained in the FFNS are compared with available results such as the general-mass variable-flavour-number scheme (GM-VFNS) and other prescriptions used in global fits of PDFs. In the present QCD analysis, we use a wide range of the inclusive neutral-current deep-inelastic-scattering (NC DIS) data, including the most recent data for charm $F_2^c$, bottom $F_2^b$, longitudinal $F_L$ structure functions and also the reduced DIS cross sections $sigma_{r,NC}^pm$ from HERA experiments. The most recent HERMES data for proton and deuteron structure functions are also added. We take into account ZEUS neutral current $e^ pm p$ DIS inclusive jet cross section data from HERA together with the recent Tevatron Run-II inclusive jet cross section data from CDF and D{O}. The impact of these recent DIS data on the PDFs extracted from the global fits are studied. We present two families of PDFs, {tt KKT12} and {tt KKT12C}, without and with HERA `combined data sets on $e^{pm}p$ DIS. We find these are in good agreement with the available theoretical models.
We present the results of our QCD analysis for non-singlet unpolarized quark distributions and structure function $F_2(x,Q^2)$ up to N$^3$LO. In this regards 4-loop anomalous dimension can be obtain from the Pade approximations. The analysis is based
on the Jacobi polynomials expansion of the structure function. New parameterizations are derived for the non-singlet quark distributions for the kinematic wide range of $x$ and $Q^2$. Our calculations for non-singlet unpolarized quark distribution functions up to N$^3$LO are in good agreement with available theoretical models. The higher twist contributions of $F_2^{p,d}(x,Q^2)$ are extracted in the large $x$ region in N$^3$LO analysis. The values of $Lambda_{QCD}$ and $alpha_s(M_z^2)$ are determined.
