  We calculate $c\bar{c}c\bar{c}$ and $b\bar{b}b\bar{b}$ production cross
section in Double Parton Scattering (DPS) processes at the energies of LHC. In
our calculation we include saturation effects in the gluon distribution. For
charm production, our results confirm the prediction of a previous work in
which saturation effects were neglected, namely that at the energies of LHC
charm production in DPS processes becomes comparable with charm production in
SPS. We also estimate the production of $c\bar{c}b\bar{b}$ in DPS. Our results
indicate that at $\sqrt{s}=14$ TeV half of the total amount of bottom produced
in proton-proton collisions comes from this process.
