Enhance propagation helper with extrapolation

This commit adds a boolean flag to our existing propagation helper used
in the global \chi^2 fitter that enables it to also do an extrapolation,
gathering not only the position where the propagation lands but also a
set of surfaces crossed in between.