Terahertz conductivity of doped silicon calculated using the ensemble Monte Carlo/finite-difference time-domain simulation technique