Add generic fit function, generalize prepare_arrays()

The idea is that the functions prepare_arrays, range_mask and func_fit would be the base for fitting of pulse-resolved point data.

@carinanc