We summarize the results of the studies of double parton scattering (DPS) off nuclei in the parton
model  and pQCD . A Monte Carlo implementation of the DPS processes is developed which
includes realistic NN correlation in nuclei , transverse geometry of hard and soft NN collisions
and distribution over the number of wounded nucleons. This implementation is an extension of the
Monte Carlo procedure previously developed for processes with one hard collisions [4, 5] and which
was applied within the framework of color fluctuation picture for the centrality dependence of the
jet production in the proton fragmentation region in pA and DA collisions[6, 7]. Numerical results
are presented for the inclusive rate of DPS as well as for the rate of DPS as a function of centrality.
Centrality dependence was calculated using the ATLAS collaboration model of centrality for soft processes based on the \sum E_T distribution at large negative (along the nucleus direction) rapidities.
 M. and Strikman and D. Treleani, Phys. Rev. Lett. 88, 031801 (2002)
 B. Blok, M. Strikman and U. A. Wiedemann, Eur. Phys. J. C 73, no. 6, 2433 (2013)
 M. Alvioli,
H.-J. Drescher and M. Strikman,
Phys. Lett. B 680,
 M. Alvioli and M. Strikman, Phys. Lett. B 722, 347 (2013) doi:10.1016/j.physletb.2013.04.042
 M. Alvioli, L. Frankfurt, V. Guzey and M. Strikman, Phys. Rev. C 90, 034914 (2014)
 M. Alvioli, B. A. Cole, L. Frankfurt, D. V. Perepelitsa and M. Strikman, Phys. Rev. C 93, no.
1, 011902 (2016) doi:10.1103/PhysRevC.93.011902
 M. Alvioli, L. Frankfurt, D. Perepelitsa and M. Strikman, Phys. Rev. D 98, no. 7, 071502