A multi-region collision probability method for determining neutron spectra and reaction rates