How can I fit a good optical Lorentzian on python using scipy.optimize.curve_fit?

Question

I\'m trying to fit a Lorentzian function with more than Imaginary part of Silicon, but the curve_fit function it not working properly, not fitting any peaks with this messag