We propose a new parton model and demonstrate that the model describes the relevant experimental data at high energies. The model is based on Pomeron calculus in 1 + 1 space-time dimensions, as suggested in Kovner et al. (JHEP 1608:031, 2016), and on simple assumptions regarding the hadron structure, related to the impact parameter dependence of the scattering amplitude. This parton model evolves from QCD, assuming that the unknown non-perturbative corrections lead to fixing the size of the interacting dipoles. The advantage of this approach is that it satisfies both t-channel and s-channel unitarity, and can be used for summing all diagrams of Pomeron interactions, including Pomeron loops. We can use this approach for all reactions: dilute–dilute (hadron–hadron), dilute–dense (hadron–nucleus) and dense–dense (nucleus–nucleus) for the scattering of parton systems. Unfortunately, we are still far from being able to tackle this problem in the effective QCD theory at high energy (i.e. in the CGC/saturation approach).