We suggest a new theoretical method to describe the baryon clusterization of nuclei in hadron-nucleus reactions. As an example we have explored the nuclei production in π−+C and π−+W collisions at plab=1.7 GeV by using the hybrid approach consisting of the Ultra Relativistic Quantum Dynamics Model (UrQMD) and the Statistical Multifragmentation Model (SMM). The UrQMD describes the production of new baryons, and the propagation toward the subnuclear densities with the fluctuations leading to the formation of excited baryonic clusters. The SMM describes the production of final nuclei and hypernuclei after interaction of baryons inside these clusters. We demonstrate the transverse momenta, rapidity, mass distributions and excitation energies of both primary clusters and final nuclei (including hypernuclei). The results of the UrQMD and UrQMD+SMM model calculations for different clusterization parameters are compared with the available HADES experimental data on baryon production, providing a very promising window for future research on nuclei and hypernuclei formation in these reactions.