Initial version with gaussian estimate
Refactor (does not work)
Integrate source.py into data_generator