Implementation of a double Gaussian source model for the BEAMnrc Monte Carlo code and its influence on small fields dose distributions