Non-orthogonal multiple access (NOMA), millimeter wave (mmWave), and massive multiple-input-multiple-output (MIMO) have been emerging as key technologies for fifth generation mobile communications. However, less studies have been done on combining the three technologies into the converged systems. In addition, how many capacity improvements can be achieved via this combination remains unclear. In this paper, we provide an in-depth capacity analysis for the integrated NOMA-mmWave-massive-MIMO systems. First, a simplified mmWave channel model is introduced by extending the uniform random single-path model with angle of arrival. Afterward, we divide the capacity analysis into the low signal to noise ratio (SNR) and high-SNR regimes based on the dominant factors of signal to interference plus noise ratio. In the noise-dominated low-SNR regime, the capacity analysis is derived by the deterministic equivalent method with the Stieltjes-Shannon transform. In contrast, the statistic and eigenvalue distribution tools are invoked for the capacity analysis in the interference-dominated high-SNR regime. The exact capacity expression and the low-complexity asymptotic capacity expression are derived based on the probability distribution function of the channel eigenvalue. Finally, simulation results validate the theoretical analysis and demonstrate that significant capacity improvements can be achieved by the integrated NOMA-mmWave-massive-MIMO systems.