This paper presents a theoretical framework for investigating the coverage and throughput behavior of sixth generation (6 G) peer-to-peer (P2P) directional slotted Aloha (DirSA) networks managing bursty traffic flows. Proper channel models, accounting for interference, noise, path-loss, random node location, power fluctuation, and beam pointing error, are adopted to derive analytical expressions for the statistic of the received power in ground, air, and space propagation contexts. The resulting coverage probability, obtained in simple integral form for different omnidirectional/directional transmission/reception modes, is exploited to derive multidimensional Markov chains for estimating the throughput in the absence and in the presence of a feedback mechanism, considering also the impact of the initial access procedure and of the beam training overhead. The theoretical results, which are validated by exhaustive Monte Carlo simulations, are used to evaluate the influence of the code-modulation scheme, of the operating signal to interference plus noise ratio (SINR), and of the burst length on the performance of 6 G terrestrial, aerial, and satellite P2P DirSA subnets.