from spectractor import parameters
from spectractor.tools import multigauss_and_bgd
import numpy as np
parameters.CALIB_BGD_NPARAMS = 4
x = np.arange(600., 800., 1)
x_norm = rescale_x_to_legendre(x)
p = [20, 0, 0, 0, 20, 650, 3, 40, 750, 5]
y = multigauss_and_bgd(np.array([x_norm, x]), *p)
plt.plot(x,y,'r-')
plt.show()