  The effect of full $7$ sets of LHC ATLAS jet cross sections data at $
\sqrt{s} = 7$ TeV on the proton parton distribution functions (PDFs) up to
next-to-next-to-next-to-leading order (NNNLO or N3LO) corrections is
investigated for the first time. Phenomenologically, the proton central PDFs in
this perturbative Quantum Chromo Dynamics (pQCD) analysis are defined based on
the full seven data sets from HERA I and II combined. It is shown that, the LHC
ATLAS jet cross sections data at $ \sqrt{s} = 7$ TeV on the HERA I and II
combined data reduces the error band of proton PDFs. Particularly, the
uncertainties of the gluon $xg(x,Q^2)$ and charm $xc(x,Q^2)$ distributions
decrease dramatically. Adding the LHC ATLAS jet cross sections data at $
\sqrt{s} = 7$ TeV on the central proton PDFs improves the quality of the fit up
to $\sim 1.53$ \%, $\sim 2.72$ \% and $\sim 2.80$ \% corresponding to
next-to-leading order (NLO), next-to-next-to leading order (N2LO) and N3LO,
respectively.
