Starting with 12 data points: Hipparcos B and V mags, 2MASS J, H, and K mags, and WISE Bands 1-4, I want to fit a stellar photosphere from the Hauschildt et al. 1999 models. The purpose of the SED fitting (in my case) is to determine if there is an excess (flux above the photosphere) at 22 micron (WISE Band 4). This can be indicative of a dusty debris disk around a star.
Downloading the models was step one - I am storing a folder on my computer which contains a set of 94 2-column text files. These text files are just the x and y coordinates for an SED. Each one corresponds to a stellar photosphere of a different temperature and logg. A simple chi-squared fit to the data points should be sufficient to predict the photosphere.
The model SEDs are not sufficient on their own, however. Each data point corresponds to a bandpass with a certain filter response function. This function must be convolved with the model SED over the bandwidth at each data point to properly predict the measurable photospheric flux at that wavelength.
This is proving to be more difficult. I am attempting to keep the relevant files contained in a single folder so that when the script is finished and working, I can send it to anyone so that they can run it on their computer easily.
No comments:
Post a Comment