Monte Carlo simulation of multilayer delta doped germanium THz laser