HiGal SED fitter
The source code is hosted on github
Installation
Can be installed through pip:
pip install https://github.com/keflavich/HiGal_SEDfitter/archive/master.zip
Example
import glob
from astropy import units as u
from astropy.io import fits
from higal_sedfitter import smooth,fit,higal_beams
from higal_sedfitter.fit import PixelFitter
pixelfitter = PixelFitter(bfixed=True)
higal_field = '000'
fmt = {'field':higal_field}
for fn in glob.glob("HIGAL{field}*fits".format(**fmt)):
smooth.add_beam_information_to_higal_header(fn, name_to_um=higal_beams.num_to_um)
target_fn = glob.glob('HIGAL{field}*_500_RM.fits'.format(**fmt))[0]
target_header = fits.getheader(target_fn)
smooth.smooth_images_toresolution(45*u.arcsec, skip_existing=False,
globs=['HIGAL{field}*'.format(**fmt)],
regrid=True,
target_header=target_header, clobber=True)
fit.fit_modified_blackbody_tofiles('HIGAL{field}*_{{0}}_RM_smregrid45.fits'.format(**fmt),
pixelfitter=pixelfitter,
name_to_um=higal_beams.num_to_um,
out_prefix='HIGAL_L{field}_'.format(**fmt),
integral=True)
Reference/API
higal_sedfitter.smooth Module
higal_sedfitter.fit Module
Classes
PixelFitter ([tguess, bguess, nguess, ...]) |
Initialize an SED fitter instance with a set of guesses. |