  Using recent and updated world data on polarized structure functions $g_1$
and $g_2$ we perform QCD analysis at the next-next-to-leading-order (NNLO)
accuracy. We include also target mass correction and higher twist effect to get
more precise results in our fitting procedure. To confirm the validity of our
fitting results several sum rules are examined and we do a comparison for them
with results from other models. In our analysis we employ Jacobi polynomials
approach to obtain analytical solutions of the DGLAP evolution equations for
parton distribution functions (PDFs). Using the extracted PDFs from our data
analysis as input we also compute the $x$- and $\mathbf{p}_{T}$-dependence of
some transverse momentum dependence (TMD) PDFs in polarized case, based on
covariant parton model. These functions are naively even time-reversal (T-even)
at twist-2 approximation. The results for TMDs are indicating proper and
acceptable behaviour with respect to what are presented in other literatures.
