Direct, on-the-fly calculation of unresolved resonance region cross sections in Monte Carlo simulations