Create phmfit.py