Hybrid beamforming (HBF) architecture provides promising trade-offs between the system performance and computational/hardware complexity in practical implementations of millimeter wave (mm-Wave) massive multiple-input multiple-output (mMIMO) 5th generation (5G) mobile networks compared to its fully digital beamforming (DBF) counterpart. In this thesis, we investigate the future applicability of deploying hybrid beamforming architectures with subarray beam pattern shaping for mm-Wave 5G base stations in spatially heterogeneous user distributions and different propagation scenarios. We propose HBF structures with a cosecant-squared pattern and a flat-top pattern as well as their HBF and DBF benchmarks. In addition to the uniform user distribution, three non-uniform user distributions, i.e., the near-site distribution, the cell-center distribution, and the cell-edge distribution are proposed to represent the traffic flow and mobility of users due to festivals and holidays. We evaluate the performance in a novel 5G new radio (NR) system-level simulation (SLS) model. Numerical results show that the HBF architecture with a cosecant-squared subarray beam pattern is more robust against differences in spatially heterogeneous traffics than the flat-top HBF and benchmark HBF under the line-of-sight (LoS) propagation scenario. Under the non-line-of-sight (NLoS) propagation scenario, more deterministic environment information and radio channel modeling are required to improve the system performance of the shaped HBF beamforming technique.