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

renormalization in calculation of DT flux #104

Open
jsitarek opened this issue Oct 12, 2021 · 0 comments
Open

renormalization in calculation of DT flux #104

jsitarek opened this issue Oct 12, 2021 · 0 comments

Comments

@jsitarek
Copy link
Collaborator

jsitarek commented Oct 12, 2021

Investigating the DT SED computation I realized that this part is quite problematic. If you compute DT SED for nu that are not sorted (e.g. if you generate lambdas and convert them to nu) the value will be wrong (I got negative). Also if the calculations are done only at course frequencies (or not covering the whole emission) this renormalization can screw up the integral emission. BB emission can be integrated analitically, so it would be better to use such analytical integration instead.

L = (np.trapz(sed_dt / nu, nu) * 4 * np.pi * d_L ** 2).to("erg s-1")

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

1 participant