  We present a comparative analysis of the impact of the non-perturbative
intrinsic charm quark content of the proton on differential cross section of
$\gamma +c $-jet in $ pp $ and $ p\bar{p} $ collisions, for the kinematic
regions that are sensitive to this contribution. We discuss the $ Q^2 $
evolution of intrinsic quark distributions and present a code which provide
these distributions as a function of $ x $ and $ Q^2 $ for any arbitrary Fock
state probability. For the $ p\bar{p} $ collisions at the Tevatron, the results
are compared with the recent experimental data of D0 at $\sqrt{s} =1.96$ TeV
and also predictions for $ pp $ collisions at $\sqrt{s} =8$ TeV and $\sqrt{s}
=13$ TeV for the LHC.
