We review the nonlinear statistics of Primordial Black Holes that form from the collapse of over-densities in a radiation-dominated Universe. We focus on the scenario in which large over-densities are generated by rare and Gaussian curvature perturbations during inflation. As new results, we show that the mass spectrum follows a power law determined by the critical exponent of the self-similar collapse up to a power spectrum dependent cutoff, and that the abundance related to very narrow power spectra is exponentially suppressed. Related to this, we discuss and explicitly show that both the Press–Schechter approximation and the statistics of mean profiles lead to wrong conclusions for the abundance and mass spectrum. Finally, we clarify that the transfer function in the statistics of initial conditions for Primordial Black Holes formation (the abundance) does not play a significant role.