+- name: str, optional
source name, Default is None. If None, it will use plateid+mjd+fiberid as the name. If there are no
+such parameters, it will be empty.
+- nsmooth: integer number, optional
do n-pixel smoothing to the raw input flux and err spectra. The default is set to 1 (no smooth).
+It will return the same array size. We note that smooth the raw data is not suggested, this function is in case of some fail-fitted low S/N spectra.
+- and_mask: bool, optional
If True, and and_mask or or_mask is not None, it will delete the masked pixels, and only return the remained pixels. Default: False
+- or_mask: bool, optional
If True, and and_mask or or_mask is not None, it will delete the masked pixels, and only return the remained pixels. Default: False
+- reject_badpix: bool, optional
reject 10 most possible outliers by the test of pointDistGESD. One important Caveat here is that this process will also delete narrow emission lines
+in some high SN ratio object e.g., [OIII]. Only use it when you are definitely clear about what you are doing. It will return the remained pixels.
+- deredden: bool, optional
correct the Galactic extinction only if the RA and Dec are available. It will return the corrected flux with the same array size. Default: True.
+- wave_range: 2-element array, optional
trim input wavelength (lam) according to the min and max range of the input 2-element array, e.g.,
+np.array([4000.,6000.]) in Rest frame range. Default: None
+- wave_mask: 2-D array
mask some absorption lines or sky lines in spectrum, e.g., np.array([[2200.,2300.]]), np.array([[5650.,5750.],[5850.,5900.]])
+- decompose_host: bool, optional
If True, the host galaxy-QSO decomposition will be applied. If no more than 100 pixels are negative, the result will be applied. The Decomposition is
+based on the PCA method of Yip et al. 2004 (AJ, 128, 585) & (128, 2603). Now the template is only available for redshift < 1.16 in specific absolute
+magnitude bins. For galaxy, the global model has 10 PCA components and first 5 will enough to reproduce 98.37% galaxy spectra. For QSO, the global model
+has 50, and the first 20 will reproduce 96.89% QSOs. If have i-band absolute magnitude, the Luminosity-redshift binned PCA components are available.
+Then the first 10 PCA in each bin is enough to reproduce most QSO spectrum. Default: False
+- host_prior: bool, optional
This parameter is only useful when the decompose_host is True and BC03 is False. If True, the code will
+adopt the prior parameters given in the pca file to perform host decomposition.
+- host_prior_scale: float, optional
If the prior decomposition is performed, the code will use this parameter to scale the prior penalty to the
+original chi2. Default: 0.2
+- host_line_mask: bool, optional
If True, the line region of galaxy will be masked when subtracted from original spectra. Default: True
+- BC03: bool, optional
if True, it will use Bruzual1 & Charlot 2003 host model to fit spectrum, high shift host will be low resolution R ~ 300, the rest is R ~ 2000. Default: False
+- Mi: float, optional
i-band absolute magnitude. It only works when decompose_host is True. If not None, the Luminosity redshift binned PCA will be used to decompose
+the spectrum. Default: None
+- npca_gal: int, optional
the number of galaxy PCA components applied. It only works when decompose_host is True. The default is 5,
+which is already account for 98.37% galaxies.
+- npca_qso: int, optional
the number of QSO PCA components applied. It only works when decompose_host is True. The default is 20,
+No matter the global or luminosity-redshift binned PCA is used, it can reproduce > 92% QSOs. The binned PCA
+is better if have Mi information.
+- Fe_uv_op: bool, optional
if True, fit continuum with UV and optical FeII template. Default: True
+- poly: bool, optional
if True, fit continuum with the polynomial component to account for the dust reddening. Default: False
+- BC: bool, optional
if True, fit continuum with Balmer continua from 1000 to 3646A. Default: False
+- rej_abs_conti: bool, optional
if True, it will iterate the continuum fitting once, rejecting 3 sigma outlier absorption pixels in the continuum
+(< 3500A), which might fall into the broad absorption lines. Default: False
+- rej_abs_line: bool, optional
if True, it will iterate the emission line fitting twice, rejecting 3 sigma outlier absorption pixels
+which might fall into the broad absorption lines. Default: False
+- n_pix_min_conti: float, optional
minimum number of negative pixels for host continuuum fit to be rejected. Default: 100
+- param_file_name: str, optional
name of the qso fitting parameter FITS file. Default: ‘qsopar.fits’
+- MC: bool, optional
if True, do Monte Carlo resampling of the spectrum based on the input error array to produce the MC error array.
+if False, the code will not save the MLE minimization error produced by lmfit since it is biased and can not be trusted.
+But it can be still output by using the lmfit attribute. Default: False
+- MCMC: bool, optional
if True, do Markov Chain Monte Carlo sampling of the posterior probability densities after MLE fitting to produce the error array.
+Note: An error will be thrown if both MC and MCMC are True. Default: False
+- nburn: int, optional
the number of burn-in samples to run MCMC chain if MCMC=True. It only works when MCMC is True. Default: 20
+- nsamp: int, optional
the number of trials of the MC process to produce the error array (if MC=True) or number samples to run MCMC chain (if MCMC=True). Should be larger than 20. It only works when either MC or MCMC is True. Default: 200
+- linefit: bool, optional
if True, the emission line will be fitted. Default: True
+- save_result: bool, optional
if True, all the fitting results will be saved to a fits file, Default: True
+- plot_fig: bool, optional
if True, the fitting results will be plotted. Default: True
+- save_fig: bool, optional
if True, the figure will be saved, and the path can be set by “save_fig_pathâ€. Default: True
+- plot_corner: bool, optinoal
whether or not to plot the corner plot results if MCMC=True. Default: True
+- save_fig_path: str, optional
the output path of the figure. If None, the default “save_fig_path†is set to “pathâ€
+- save_fits_path: str, optional
the output path of the result fits. If None, the default “save_fits_path†is set to “pathâ€
+- save_fits_name: str, optional
the output name of the result fits. Default: “result.fitsâ€
+- verbose: bool, optional
turn on (True) or off (False) debugging output. Default: False
+- kwargs_plot: dict, optional
extra aguments for plot_fig for plotting results. See LINK TO PLOT_FIG_DOC. Default: {}
+- kwargs_conti_emcee: dict, optional
extra aguments for emcee Sampler for continuum fitting. Default: {}
+- kwargs_line_emcee: dict, optional
extra arguments for emcee Sampler for line fitting. Default: {}
+
+