Line fitting for synthetic spectra using pymc3

You can take a look of Fitting a spectra of gaussians which is a similar question.