This paper evaluates the downlink performance of cellular networks in terms of coverage and electromagnetic field exposure (EMFE), in the framework of stochastic geometry. The model is constructed based on datasets for sub-6~GHz macro cellular networks but it is general enough to be applicable to millimeter-wave networks as well. On the one hand, performance metrics are calculated for $β$-Ginibre point processes which are shown to faithfully model a large number of motion-invariant networks. On the other hand, performance metrics are derived for inhomogeneous Poisson point processes with a radial intensity measure, which are shown to be a good approximation for motion-variant networks. For both cases, joint and marginal distributions of the EMFE and the coverage, and the first moments of the EMFE are provided and validated by Monte Carlo simulations using realistic sets of parameters from two sub-6~GHz macro urban cellular networks, i.e., 5G~NR~2100 (Paris, France) and LTE~1800 (Brussels, Belgium) datasets. In addition, this paper includes the analysis of the impact of the network parameters and discusses the achievable trade-off between coverage and EMFE.