Fitting a crystal ball function

Question

I\'m currently trying to fit the following function using scipy.optimize.curve_fit

def func(x,mu,sigma,A,B,C):
    g = A*np.exp(-(x-mu)**2/(2*sigma**2))
    e