  We present predictions for the double parton scattering (DPS) four-jet
production cross sections in $pA$ collisions at the LHC. Relying on the
experimental capabilities to correlate centrality with impact parameter $B$ of
the proton-nucleus collision, we discuss a strategy to extract the double
parton scattering contributions in $pA$ collisions, which gives direct access
to double parton distribution in the nucleon. We show that the production cross
sections via DPS of four jets, out of which two may be light- or heavy-quark
jets, are large enough to allow the method to be used already with data
accumulated in 2016 $pA$ run.
