A stellar abundance matching code
Project description
StarFit
Match stellar abundance measurements against a database of stellar models.
Installation
From PyPI
pip install starfit
From git repo
Download data files before installing from the repo
./download-data.sh
pip install .
Usage
starfit.Single
fits an abundance pattern to a single model from the database.
Optional arguments:
z_max
: highest element charge number to fitz_exclude
: element charge numbers to exclude from fitz_lolim
: elements that are model lower limits (effectively the same as observational upper limits)
import starfit
s = starfit.Single(
filename = 'HE1327-2326.dat',
db = 'znuc2012.S4.star.el.y.stardb.gz',
z_max = 30,
z_exclude = [3, 24, 30],
z_lolim = [21, 29],
)
print(s)
The database indices of the best fitting models (sorted from best to worst) are given by:
s.sorted_stars['index']
The corresponding reduced chi^2 values are:
s.sorted_fitness
The physical properties of the models corresponding to these indices can be accessed using the database:
i_bestfit = s.sorted_stars['index'][0]
s.db.fielddata[i_bestfit]
The chemical yield of the models (and the respective element names) are:
list( zip(s.list_db, s.full_abudata[:,i_bestfit]) )
To make the same plots as the web version:
s.plot()
Project details
Release history Release notifications | RSS feed
Download files
Download the file for your platform. If you're not sure which to choose, learn more about installing packages.
Source Distribution
starfit-0.4.3.tar.gz
(14.9 MB
view hashes)