Neutronic calculations for CANDU thorium systems using Monte Carlo techniques