In this paper, we will expose the work performed for the adaptation and refinement of the GEANT4 simulation toolkit (originally designed for Monte Carlo simulations in High-Energy Physics) in order to simulate Ge detectors in low-level gamma spectrometry. Special emphasis will be given to show and validate our own algorithms implemented in GEANT4 code for variance reduction and data analysis, that have been used for a quicker and precise efficiency calibration at different source-detector configurations of one of the low-level germanium gamma systems available in our laboratory.