Full dimension (FD) multiple-input multiple-output (MIMO) technology has attracted substantial research attention from both wireless industry and academia in the last few years as a promising technique for next-generation wireless communication networks. FD-MIMO scenarios utilize a planar two-dimensional (2D) active antenna system (AAS) that not only allows a large number of antenna elements to be placed within feasible base station (BS) form factors, but also provides the ability of adaptive electronic beam control over both the elevation and the traditional azimuth dimensions. This paper presents a tutorial on elevation beamforming analysis for cellular networks utilizing FD Massive MIMO antenna arrays. In contrast to existing works that focus on the standardization of FD-MIMO in the 3rd Generation Partnership Project (3GPP), this tutorial is distinguished by its depth with respect to the theoretical aspects of antenna array and 3D channel modeling. In an attempt to bridge the gap between industry and academia, this preliminary tutorial introduces the relevant array and transceiver architecture designs proposed in the 3GPP Release 13 that enable elevation beamforming. Then it presents and compares two different 3D channel modeling approaches that can be utilized for the performance analysis of elevation beamforming techniques. The spatial correlation in FD-MIMO arrays is characterized and compared based on both channel modeling approaches and some insights into the impact of different channel and array parameters on the correlation are drawn. All these aspects are put together to provide a mathematical framework for the design of elevation beamforming schemes in single-cell and multi-cell scenarios. Simulation examples associated with comparisons and discussions are also presented. To this end, this paper highlights the state-of-the-art research and points out future research directions.