Chaotic waveforms with Gaussian distributions are significant for laser-chaos-based applications such as random number generation. By exploring the injection parameter space of the optical injection semiconductor lasers, we numerically investigate the associated probability density functions of the generated chaotic waveforms when different high-pass filters with different cutoff frequencies are used. Our results demonstrate that the chaotic waveforms with Gaussian probability density functions can be obtained once the cutoff frequency of the high-pass filter is larger than the laser relaxation resonance frequency. Especially, we find that the Gaussian probability density function can reach a superhigh coefficient of determination R2 ≥ 99.5% and an ultralow skewness |S|<0.1 in a large parameter space by jointly controlling the injection parameter and cutoff frequency.