BEAST-Fitting / megabeast

Hierarchical Bayesian Model for Ensembles of Dust Extinguished Stellar Populations
http://megabeast.readthedocs.io/
2 stars 12 forks source link

simulations of A_V #31

Closed lea-hagen closed 5 years ago

lea-hagen commented 5 years ago

This creates a simulated A_V dataset, in which the underlying A_V can be a single or double lognormal. Currently, the simulated A_V values have an error of 0.2 and a completeness of 1, which will need to be generalized in the future.

The inputs are the SED grid file, the parameters of the lognormal(s), and a label for the created files. The dimensions of the image and the number of stars per pixel are also optional inputs. The code creates all necessary data files (lnp, nstars, noisemodel) and an input file for the megabeast.

Documentation is included.

Partially addresses #6.

lea-hagen commented 5 years ago

Additional update from this PR: while working on this, we figured out that the normalization of the two lognormals (N1 and N2) are not actually constrained by the current fitting procedure. (Larger N1 orN2 just increased the likelihood of the fit.) Instead, the MegaBEAST now fits for N12_ratio, which is the ratio of the two Ns. Sometime in the future, we'll work through the math to also fit for the total predicted number of stars.

I've also added code (parameter_maps.py) to create maps of the best-fit parameters.

karllark commented 5 years ago

This pull request introduces 7 alerts when merging c66073cb62734a83c9988adf8d0b32098097a0b5 into 7b8914c41a9fcae6b63362aee0f113a498679655 - view on LGTM.com

new alerts:


Comment posted by LGTM.com

karllark commented 5 years ago

This pull request introduces 7 alerts when merging 9a6593b2158487171ce9604d0c1938ea0dd7b3fe into 7b8914c41a9fcae6b63362aee0f113a498679655 - view on LGTM.com

new alerts:


Comment posted by LGTM.com

lea-hagen commented 5 years ago

Now added plotting code (simulate_av_plots.py) to compare the simulated A_V distributions with their assumed model and the best fit.

karllark commented 5 years ago

This pull request introduces 1 alert when merging 766b7bc430da15782835cececffd052630f138c2 into 7b8914c41a9fcae6b63362aee0f113a498679655 - view on LGTM.com

new alerts:


Comment posted by LGTM.com

karllark commented 5 years ago

This pull request introduces 1 alert when merging 02f7b44053b02b3df0e25239a8d43764081968e5 into 7b8914c41a9fcae6b63362aee0f113a498679655 - view on LGTM.com

new alerts:


Comment posted by LGTM.com

lea-hagen commented 5 years ago

@karllark this is finally ready to merge!

karllark commented 5 years ago

Merging.