Dustin Lang

emcee: The MCMC Hammer

We introduce a stable, well tested Python implementation of the affine-invariant ensemble sampler for Markov chain Monte Carlo (MCMC) proposed by Goodman & Weare (2010). The code is open source and has already been used in several published projects in the astrophysics literature. The algorithm behind emcee has several advantages over traditional MCMC sampling methods and it has excellent performance as measured by the autocorrelation time (or function calls per independent sample). One major advantage of the algorithm is that it requires hand-tuning of only 1 or 2 parameters compared to ~N2 for a traditional algorithm in an N-dimensional parameter space. In this document, we describe the algorithm and the details of our implementation. Exploiting the parallelism of the ensemble method, emcee permits any user to take advantage of multiple CPU cores without extra effort. The code is available online at http://dan.iel.fm/emcee under the GNU General Public License v2.

THE PANCHROMATIC HUBBLE ANDROMEDA TREASURY

The Panchromatic Hubble Andromeda Treasury is an ongoing Hubble Space Telescope Multi-Cycle Treasury program to image ~1/3 of M31s star-forming disk in six filters, spanning from the ultraviolet (UV) to the near-infrared (NIR). We use the Wide Field Camera 3 (WFC3) and Advanced Camera for Surveys (ACS) to resolve the galaxy into millions of individual stars with projected radii from 0 to 20 kpc. The full survey will cover a contiguous 0.5 deg^(2)area in 828 orbits. Imaging is being obtained in the F275W and F336W filters on the WFC3/UVIS camera, F475W and F814W on ACS/WFC, and F110W and F160W on WFC3/IR. The resulting wavelength coverage gives excellent constraints on stellar temperature, bolometric luminosity, and extinction for most spectral types. The data produce photometry with a signal-to-noise ratio of 4 at m F_(275W) = 25.1, m_(F336W) = 24.9, m_(F475W) = 27.9, m_(F814W) = 27.1, m_(F110W) = 25.5, and m_(F160W) = 24.6 for single pointings in the uncrowded outer disk; in the inner disk, however, the optical and NIR data are crowding limited, and the deepest reliable magnitudes are up to 5 mag brighter. Observations are carried out in two orbits per pointing, split between WFC3/UVIS and WFC3/IR cameras in primary mode, with ACS/WFC run in parallel. All pointings are dithered to produce Nyquist-sampled images in F475W, F814W, and F160W. We describe the observing strategy, photometry, astrometry, and data products available for the survey, along with extensive testing of photometric stability, crowding errors, spatially dependent photometric biases, and telescope pointing control. We also report on initial fits to the structure of M31s disk, derived from the density of red giant branch stars, in a way that is independent of assumed mass-to-light ratios and is robust to variations in dust extinction. These fits also show that the 10 kpc ring is not just a region of enhanced recent star formation, but is instead a dynamical structure containing a significant overdensity of stars with ages >1 Gyr.

Fast particle smoothing: if I had a million particles

We propose efficient particle smoothing methods for generalized state-spaces models. Particle smoothing is an expensive O(N2) algorithm, where N is the number of particles. We overcome this problem by integrating dual tree recursions and fast multipole techniques with forward-backward smoothers, a new generalized two-filter smoother and a maximum a posteriori (MAP) smoother. Our experiments show that these improvements can substantially increase the practicality of particle smoothing.

The Hyper Suprime-Cam software pipeline

In this paper, we describe the optical imaging data processing pipeline developed for the Subaru Telescopes Hyper Suprime-Cam (HSC) instrument. The HSC Pipeline builds on the prototype pipeline being developed by the Large Synoptic Survey Telescopes Data Management system, adding customizations for HSC, large-scale processing capabilities, and novel algorithms that have since been reincorporated into the LSST codebase. While designed primarily to reduce HSC Subaru Strategic Program (SSP) data, it is also the recommended pipeline for reducing general-observer HSC data. The HSC pipeline includes high level processing steps that generate coadded images and science-ready catalogs as well as low-level detrending and image characterizations.

PHAT stellar cluster survey. I. Year 1 catalog and integrated photometry

The Panchromatic Hubble Andromeda Treasury (PHAT) survey is an ongoing Hubble Space Telescope (HST) multi-cycle program to obtain high spatial resolution imaging of one-third of the M31 disk at ultraviolet through near-infrared wavelengths. In this paper, we present the first installment of the PHAT stellar cluster catalog. When completed, the PHAT cluster catalog will be among the largest and most comprehensive surveys of resolved star clusters in any galaxy. The exquisite spatial resolution achieved with HST has allowed us to identify hundreds of new clusters that were previously inaccessible with existing ground-based surveys. We identify 601 clusters in the Year 1 sample, representing more than a factor of four increase over previous catalogs within the current survey area (390 arcmin2). This work presents results derived from the first ~25% of the survey data; we estimate that the final sample will include ~2500 clusters. For the Year 1 objects, we present a catalog with positions, radii, and six-band integrated photometry. Along with a general characterization of the cluster luminosities and colors, we discuss the cluster luminosity function, the cluster size distributions, and highlight a number of individually interesting clusters found in the Year 1 search.

Quasar probabilities and redshifts from WISE mid-IR through GALEX UV photometry

Extreme deconvolution (XD) of broad-band photometric data can both separate stars from quasars and generate probability density functions for quasar redshifts, while incorporating flux uncertainties and missing data. Mid-infrared photometric colors are now widely used to identify hot dust intrinsic to quasars, and the release of all-sky WISE data has led to a dramatic increase in the number of IR-selected quasars. Using forced-photometry on public WISE data at the locations of SDSS point sources, we incorporate this all-sky data into the training of the XDQSOz models originally developed to select quasars from optical photometry. The combination of WISE and SDSS information is far more powerful than SDSS alone, particularly at

THE PANCHROMATIC HUBBLE ANDROMEDA TREASURY. IV. A PROBABILISTIC APPROACH TO INFERRING THE HIGH-MASS STELLAR INITIAL MASS FUNCTION AND OTHER POWER-LAW FUNCTIONS*

z>2

The panchromatic hubble andromeda treasury. VIII. A wide-area, high-resolution map of dust extinction in M31

. The use of SDSS

HIERARCHICAL PROBABILISTIC INFERENCE OF COSMIC SHEAR

+

Cleaning the Usno-B Catalog Through Automatic Detection of Optical Artifacts

WISE photometry is comparable to the use of SDSS