Modularize fit function

[StillerPatrick]

Currently, the fit function is a very big function that includes tracking, fitting, parallelization, etc. Here we should modularize the functionalities to get a more maintainable code.