add some comments and try to implement fit function for multiple peaks