Device density in cellular networks is expected to increase considerably in the next future. Accordingly, the access point (AP) will equip massive multiple-input multiple-output (mMIMO) antennas, using collimated millimeter-wave (mmW) and sub-THz communications, and increase the bandwidth to accommodate the growing data rate demands. In this scenario, interference plays a critical role and, if not characterized and mitigated properly, might limit the performances of the network. In this context, this paper derives the statistical properties of the aggregated interference power for a cellular network equipping a mMIMO cylindrical array. The proposed statistical model considers the link blockage and other network parameters such as antenna configuration and device density. The findings show that the characteristic function (CF) of the aggregated interference power can be regarded as a weighted mixture of two alpha-stable distributions. Furthermore, by analyzing the service probability, it is found that there is an optimal configuration of the array depending on the AP height and device density. The proposed statistical model can be part of the design of dense networks providing valuable insights for optimal network deployment