In this paper, a simple analytical model for the Gaussian’s peak response part of the timing jitter of single photon avalanche diodes (SPADs) is proposed using Fourier series in the multiplication time calculation. The multiplication time characterizes avalanche multiplication process speed in which low multiplication time suggests a swifter response time and a higher avalanche speed. This paper presents an analytical solution which results in a more accurate multiplication time. The model is verified for SPADs implemented in 0.15 and 0.18 μm standard CMOS process, and the accuracy of the proposed analytical method in full-width at half-maximum (FWHM) calculation is improved by 25% and 5% with respect to the numerical model, respectively.