[yt-users] Fitting an absorption spectrum

[Pengfei Chen]

Hi Hilary,

Thank you very much for your update and I'm sorry for the late reply. I've
tested the new code with two datasets. The results are attached. The
abnormal peaks in the statistics of N have disappeared, and the abnormal
peaks for b are much lower now, but still visible.

Thank you very much for your effort in fixing this! Please tell me if you
need any test data from me.

Pengfei


On Fri, Feb 14, 2014 at 3:30 PM, Pengfei Chen <madcpf at gmail.com> wrote:

> Hi Hilary,
>
> Thank you! I'm looking forward to it.
>
> Pengfei
>
>
> On Thu, Feb 13, 2014 at 7:32 PM, Hilary Egan <hilaryye at gmail.com> wrote:
>
>> Hi Pengfei,
>>
>> I have a fix for this issue that just needs a little more testing! I'll
>> push it shortly and reply to this message when it's been accepted.
>>
>> -Hilary
>>
>>
>> On Thu, Feb 13, 2014 at 8:04 PM, Pengfei Chen <madcpf at gmail.com> wrote:
>>
>>> Hi all,
>>>
>>> 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!
>>> Pengfei
>>>
>>> _______________________________________________
>>> yt-users mailing list
>>> yt-users at lists.spacepope.org
>>> http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org
>>>
>>>
>>
>> _______________________________________________
>> yt-users mailing list
>> yt-users at lists.spacepope.org
>> http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org
>>
>>
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.spacepope.org/pipermail/yt-users-spacepope.org/attachments/20140304/8928a13a/attachment.htm>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: sed_f_N_vs_N_100_4.7_40_Mar03.png
Type: image/png
Size: 20463 bytes
Desc: not available
URL: <http://lists.spacepope.org/pipermail/yt-users-spacepope.org/attachments/20140304/8928a13a/attachment.png>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: uvb_f_N_vs_N_100_3.0_40_Feb27.png
Type: image/png
Size: 23099 bytes
Desc: not available
URL: <http://lists.spacepope.org/pipermail/yt-users-spacepope.org/attachments/20140304/8928a13a/attachment-0001.png>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: uvb_f_b_vs_b_100_3.0_40_Feb27.png
Type: image/png
Size: 25789 bytes
Desc: not available
URL: <http://lists.spacepope.org/pipermail/yt-users-spacepope.org/attachments/20140304/8928a13a/attachment-0002.png>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: sed_f_b_vs_b_100_4.7_40_Mar03.png
Type: image/png
Size: 22610 bytes
Desc: not available
URL: <http://lists.spacepope.org/pipermail/yt-users-spacepope.org/attachments/20140304/8928a13a/attachment-0003.png>