Levenberg-Marquardt curve fitting Python with several datasets

Question

I am implementing the Levenberg-Marquardt method to fit a series of gaussian curves in a single dataset. The plot of the data looks like this: Data before curve fitting