Mariska Kriek
University of California, Berkeley
Mariska Kriek
University of California, Berkeley
FAST
FAST (Fitting and Assessment of Synthetic Templates) is an IDL-based code that fits stellar population synthesis templates to broadband photometry and/or spectra. FAST is compatible with the photometric redshift code EAzY (Brammer et al. 2008) when fitting broadband photometry; it uses the photometric redshifts derived by EAzY, and the input files (photometric catalog, master filter file etc.) are the same. FAST also fits spectra, optionally in combination with broadband photometric data points. Depending on the input parameters, FAST outputs the best-fit redshift, age, dust content, star formation timescale, metallicity, stellar mass, star formation rate (SFR), and their confidence intervals. The main difference with HYPERZ is that (1) FAST fits fluxes instead of magnitudes, (2) you can completely define your own grid of input stellar population parameters, (3) you can easily input photometric redshifts and their confidence intervals, and (4) FAST calculates calibrated confidence intervals for all parameters. However, note that, although it can be used as one, FAST is not a photometric redshift code
How does FAST work?
FAST reads in a parameter file (see example), and makes a cube of model fluxes for the full stellar population grid, and all filters and/or spectral elements. To determine the best-fit parameters, it simply uses χ2 fitting. To avoid skipping over (multiple) minimums, FAST uses no minimum searching algorithm, but fits every point of the model cube. In case spectroscopic or photometric redshifts are provided, the redshift will be fixed to the closest value in the grid.
The confidence levels are calibrated using Monte Carlo simulations. The observed fluxes are modified according to their photometric errors, and these modified fluxes are fitted as well. The 68% (95% or 99%) confidence intervals are defined by the χ2 value in the original grid that encloses 68% (95% or 99%) of these simulations. Thus, the confidence intervals on all properties are the minimum and maximum values allowed by this χ2 threshold. In case photometric redshifts (as provided by EAzY) are assumed, the calculation of the confidence intervals is bit more complicated. In the Appendix of Kriek et al. (2009) you can find more about this issue.
Using FAST
If you are interested in using FAST, please send me an email and I will send you the password. Please do not distribute the code yourself, I would like to keep track of who is using it, so I can update everyone easily in case there is a new version available. Moreover, if you use the code, please cite the following paper (the description of the code can be found in the Appendix): Kriek et al. (2009). Any questions and remarks are welcome, but please first check the documentation in the parameter file and the FAQ page.
Download
The most recent version can be downloaded here. I still hope to make the following modifications in the future:
✦FAST uses the Calzetti et al. (2000) reddening law. I am planning to add more dust laws.
✦The stellar population synthesis libraries used by FAST do not contain emission lines. I plan to add libraries with templates including emission lines.
✦Add the flexible SPS models by Conroy et al. (2009)