SIMSTACK for Python
An open-source, Python version of SIMSTACK.
For those familiar with the IDL version, viero_quick_stack is a similar wrapper to the stack_in_redshift_slices function.
The python version of the code is intended to be object oriented, with map objects (containing maps, noisemaps, beams, color corrections, etc.) and catalog objects (with RA/DEC lists for different selection criteria) collected into libraries, and SIMSTACK performed with one line of code.
As of 1/28/2017 SIMSTACK just got easier to use! Edit the parameter file to set the paths to your maps and catalogs, as well as to choose the binning and method of stacking (use example.cfg as a guide) and run from the command line with: ../run_simstack_cmd_line.py example.cfg
Future projects include adding priors from SIMSTACK, and moving from a least-squares solution to marginalized probabilities using a full MCMC.
SEDSTACK is an attempt to fit entire SEDs to the full set of maps --- rather than fitting for flux densities one wavelength at a time --- in order to dig still deeper into the confusion-dominated data. Outstanding issues are dealing with systematic offsets, covariance between maps, careful consideration of beam and pixel size differences and how the affect the chi-squared, and using actual beams instead of Gaussian approximations.
Hopefully, with community input (which I encourage!) the code will grow in utility.
- Marco Viero
- Lorenzo Moncelsi
This is work in progress! If you use any of the code or ideas here in your research, please cite us as (Viero et al. in preparation).
All content Copyright 2015 the authors. The code is available for use under the MIT license.