I've used the
yt.analysis_modules.absorption_spectrum.absorption_spectrum_fit
script to fit an Lya absorption spectrum generated by yt.
But there are some abnormal peaks when I did statistics of line width b and column density N. Two figures are attached. I've tried with different redshifts and all of them had peaks in b-distribution at init_b (which is 40km/s here) and maxb(which is 300km/s here). For the N-distribution I got peaks around 10^12 cm^-2 at every redshift. The parameters I used were:
HI_parameters = {'name':'HI',
'f': [.4164],
'Gamma':[6.265E8],
'wavelength':[1215.67],
'numLines':1,
'maxN': 1E22, 'minN':1E11,
'maxb': 300, 'minb':1,
'maxz': 5, 'minz':0,
'init_b':40,
'init_N':1E14}
Could someone please tell me how to solve this? All the analysis were done with the latest development version of yt and I made no changes to the source code.
Please tell me if you need any datasets.
Thank you!