Massive connections and high data rate services are key players in 5G ecosystem and beyond. To satisfy the requirements of these types of services, non orthogonal multiple Access (NOMA) and 3-dimensional beam-forming (3DBF) can be exploited. In this paper, we devise a novel joint radio resource allocation and 3D multiple input single output (MISO) BF algorithm in NOMA-based heterogeneous networks (HetNets) at which our main aim is to maximize the profit of mobile virtual network operators (MVNOs). To this end, we consider multiple infrastructure providers (InPs) and MVNOs serving multiple users. Each InP has multiple access points as base stations (BSs) with specified spectrum and multi-beam antenna array in each transmitter that share its spectrum with users by employing NOMA. To realize this, we formulate a novel optimization problem at which the main aim is to maximize the revenue of MVNOs, subject to resource limitations and quality of service (QoS) constraints. Since our proposed optimization problem is non-convex, NP-hard and mathematically intractable, we transform it into a convex one by introducing a new optimization variable and converting the variables with adopting successive convex approximation. More importantly, the proposed solution is assessed and compared with the alternative search method and the optimal solution that is obtained with adopting the exhaustive search method. In addition, the proposed algorithm is studied from the computational complexity, convergence, and performance perspective. Our simulation results demonstrate that NOMA-3DBF has better performance and increases system throughput and revenue of MVNOs compared to orthogonal multiple access with 2DBF by approximately 64%. Especially, by exploiting 3DBF the revenue of MVNOs is improved by nearly 27% in contrast to 2DBF in a high order of antennas.