We identify the primary and secondary Debye relaxation times using the continuous relaxation time distribution method, which we subsequently utilize for discrete model fitting. This approach allows us to obtain accurate Debye relaxation times while minimizing thermal noises and correlations between dielectric relaxations and resonances. We use the Zasetsky-Buchner algorithm implemented in the PyOECP package to obtain g(τ ).