Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

SED bulge under 2000Å #185

Open
ZhihaoChen5 opened this issue Aug 7, 2020 · 2 comments
Open

SED bulge under 2000Å #185

ZhihaoChen5 opened this issue Aug 7, 2020 · 2 comments

Comments

@ZhihaoChen5
Copy link

Hi all, I ran Mosfit on some SLSN-I objects using dynesty sampler and my version of Mosfit is 1.1.7. I output the SEDs with command '-x seds sample_wavelengths band_indices'. The SEDs do have a suppression under 3000Å but an unusual bump appears under 2000Å which shouldn't happen on SLSN-I especially early time. It seems the higher the temperature is, the more prominent the bulge is. I thought Mosfit was using linear suppression on SEDs under wavelength_cutoff, but such bulges cannot be explained. Does anyone know which kind of suppression function Mosfit is using?

The blue points are nearby input data, and the grey lines are the SED outputs.
屏幕快照 2020-08-07 14 46 01
屏幕快照 2020-08-07 14 47 51

@mnicholl
Copy link
Collaborator

mnicholl commented Aug 7, 2020

Hi, I am not able to reproduce this issue. I ran the default SLSN model with typical parameters:

mosfit -m slsn -F mejecta 5 vejecta 8000 Pspin 2.5 Bfield 1 texplosion 0 nhhost 1e18 --band-list UVW2 UVM2 UVW1 u g r i z --band-systems AB -x seds sample_wavelengths all_bands

Using the all_bands output to determine the order in which the SEDs are printed, and plotting these versus the sample wavelengths yields the following plot, where the SED looks correct. Maybe there is a problem with the plotting rather than the outputs? Otherwise we may need more detail to diagnose this problem.

sed-test-slsn

@ZhihaoChen5
Copy link
Author

ZhihaoChen5 commented Aug 8, 2020

Hi, thanks for replying me.

I can reproduce the similar figure in your email. Like this, for different samplers at phase=50.0,
image

But I also notice some wired figures at other epoches where wavelength_cutoff is not 3000A but around 6000A, like this, at phase=52.10
image

It seems my code still have some issues when plotting the sed.

May I take a look at your code to plot this figure?

Thank you!
Zhihao

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

No branches or pull requests

2 participants