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.

Sign up for GitHub

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

Jump to bottom

Calculate integrated intensity #5

Open
low-sky opened this issue Mar 6, 2015 · 1 comment
Open

Calculate integrated intensity #5

low-sky opened this issue Mar 6, 2015 · 1 comment

Comments

@low-sky
Copy link
Contributor

low-sky commented Mar 6, 2015

We'd like to calculate the integrated intensity of W11, W22 and W33 for the model

  • Get a model spectrum out (??)
  • np.sum(ModelSpectrum) * size_of_channel
  • Then sum up significant parts of the spectrum.
@low-sky
Copy link
Contributor Author

low-sky commented Mar 13, 2015

Good so far!

  • Calculate the W11 empirically over the observed data where the model is > than, say 1e-6
  • Calculate the typical root-mean-squared error in a spectrum using np.nanstd or np.std ignoring negative values or the median absolute deviation code here. This gives the RMS in a channel which propagates to W11 empirical.
  • Store output data in an astropy Table.
  • Identify bad fits based differences between W11_theory and W11_empirical.

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Assignees
No one assigned
Labels
None yet
Projects
None yet
Milestone
No milestone
Development

No branches or pull requests
1 participant
@low-sky