dyson
is a package which takes spectra with similar inferred properties, and stacks their fitting residuals together. Spectra can be binned in arbitrary dimension.
The overall workflow is as follows:
- define binning space, consisting of many hyperrectangular bins
B
- define a wavelength grid
L
, having an identical log-wavelength spacing to the data - For each galaxy
G
in a list, assign each spaxelS
within to a binB
, and write out the assignment map to a filePLATEIFU-assign.fits
- for each bin
B
, loop through list of galaxiesG
, and find spaxelsS
inG
which fit inB
. For each spaxelS
: a) find the observed-frame stellar-continuum-fit residual belonging toS
. Deredshift to rest-frame, and drizzle intoL
(along with weight vector) b) add rest-frame, rectified residual (and weight) to running total for bin c) add spaxel redshifts to running list d) write running flux sum, variance, and redshift list to a fileBINNUM-accum.fits