Darth Fader: Using wavelets to obtain accurate redshifts of spectra at very low signal-to-noise

We present the DARTH FADER algorithm, a new wavelet-based method for estimating redshifts of galaxy spectra in spectral surveys that is particularly adept in the very low SNR regime. We use a standard cross-correlation method to estimate the redshifts of galaxies, using a template set built using a PCA analysis on a set of simulated, noise-free spectra. Darth Fader employs wavelet filtering to both estimate the continuum & to extract prominent line features in each galaxy spectrum. A simple selection criterion based on the number of features present in the spectrum is then used to clean the catalogue: galaxies with fewer than six total features are removed as we are unlikely to obtain a reliable redshift estimate. Applying our wavelet-based cleaning algorithm to a simulated testing set, we successfully build a clean catalogue including extremely low signal-to-noise data (SNR=2.0), for which we are able to obtain a 5.1% catastrophic failure rate in the redshift estimates (compared with 34.5% prior to cleaning). We also show that for a catalogue with uniformly mixed SNRs between 1.0 & 20.0, with realistic pixel-dependent noise, it is possible to obtain redshifts with a catastrophic failure rate of 3.3% after cleaning (as compared to 22.7% before cleaning). Whilst we do not test this algorithm exhaustively on real data, we present a proof of concept of the applicability of this method to real data, showing that the wavelet filtering techniques perform well when applied to some typical spectra from the SDSS archive. The Darth Fader algorithm provides a robust method for extracting spectral features from very noisy spectra. The resulting clean catalogue gives an extremely low rate of catastrophic failures, even when the spectra have a very low SNR. For very large sky surveys, this technique may offer a significant boost in the number of faint galaxies with accurately determined redshifts.Comment: 22 pages, 15 figures. Accepted for publication in Astronomy & Astrophysic

Catastrophic Photo-z Errors and the Dark Energy Parameter Estimates with Cosmic Shear

We study the impact of catastrophic errors occurring in the photometric redshifts of galaxies on cosmological parameter estimates with cosmic shear tomography. We consider a fiducial survey with 9-filter set and perform photo-z measurement simulations. It is found that a fraction of 1% galaxies at z_{spec}~0.4 is misidentified to be at z_{phot}~3.5. We then employ both chi^2 fitting method and the extension of Fisher matrix formalism to evaluate the bias on the equation of state parameters of dark energy, w_0 and w_a, induced by those catastrophic outliers. By comparing the results from both methods, we verify that the estimation of w_0 and w_a from the fiducial 5-bin tomographic analyses can be significantly biased. To minimize the impact of this bias, two strategies can be followed: (A) the cosmic shear analysis is restricted to 0.5<z<2.5 where catastrophic redshift errors are expected to be insignificant; (B) a spectroscopic survey is conducted for galaxies with 3<z_{phot}<4. We find that the number of spectroscopic redshifts needed scales as N_{spec} \propto f_{cata}\times A where f_{cata}=1% is the fraction of catastrophic redshift errors (assuming a 9-filter photometric survey) and A is the survey area. For A=1000 {deg}^2, we find that N_{spec}>320 and 860 respectively in order to reduce the joint bias in (w_0,w_a) to be smaller than 2\sigma and 1\sigma. This spectroscopic survey (option B) will improve the Figure of Merit of option A by a factor \times 1.5 thus making such a survey strongly desirable.Comment: 25 pages, 9 figures. Revised version, as accepted for publication in Ap

SuperNova Acceleration Probe (SNAP): Investigating Photometric Redshift Optimization

The aim of this paper is to investigate ways to optimize the accuracy of photometric redshifts for a SNAP like mission. We focus on how the accuracy of the photometric redshifts depends on the magnitude limit and signal-to-noise ratio, wave-length coverage, number of filters and their shapes and observed galaxy type. We use simulated galaxy catalogs constructed to reproduce observed galaxy luminosity functions from GOODS, and derive photometric redshifts using a template fitting method. By using a catalog that resembles real data, we can estimate the expected number density of galaxies for which photometric redshifts can be derived. We find that the accuracy of photometric redshifts is strongly dependent on the signal-to-noise (S/N) (i.e., S/N>10 is needed for accurate photometric redshifts). The accuracy of the photometric redshifts is also dependent on galaxy type, with smaller scatter for earlier type galaxies. Comparing results using different filter sets, we find that including the U-band is important for decreasing the fraction of outliers, i.e., ``catastrophic failures''. Using broad overlapping filters with resolution ~4gives better photometric redshifts compared to narrower filters (resolution >~5) with the same integration time. We find that filters with square response curves result in a slightly higher scatter, mainly due to a higher fraction of outliers at faint magnitudes. We also compare a 9-filter set to a 17-filter set, where we assume that the available exposure time per filter in the latter set is half that of the first set. We find that the 9-filter set gives more accurate redshifts for a larger number of objects and reaches higher redshift, while the 17-filter set is gives better results at bright magnitudes.Comment: 30 pages, 10 figures. Submitted to A

Optimising Spectroscopic and Photometric Galaxy Surveys: Efficient Target Selection and Survey Strategy

The next generation of spectroscopic surveys will have a wealth of photometric data available for use in target selection. Selecting the best targets is likely to be one of the most important hurdles in making these spectroscopic campaigns as successful as possible. Our ability to measure dark energy depends strongly on the types of targets that we are able to select with a given photometric data set. We show in this paper that we will be able to successfully select the targets needed for the next generation of spectroscopic surveys. We also investigate the details of this selection, including optimisation of instrument design and survey strategy in order to measure dark energy. We use color-color selection as well as neural networks to select the best possible emission line galaxies and luminous red galaxies for a cosmological survey. Using the Fisher matrix formalism we forecast the efficiency of each target selection scenario. We show how the dark energy figures of merit change in each target selection regime as a function of target type, survey time, survey density and other survey parameters. We outline the optimal target selection scenarios and survey strategy choices which will be available to the next generation of spectroscopic surveys.Comment: 16 pages, 22 figures, accepted to MNRAS in dec 201

CLASH: Photometric redshifts with 16 HST bands in galaxy cluster fields

Context. The Cluster Lensing And Supernovae survey with Hubble (CLASH) is a Hubble Space Telescope (HST) Multi-Cycle Treasury programme that observes 25 massive galaxy clusters, 20 of which were X-ray-selected to preferably choose dynamically relaxed clusters, and 5 additional “high magnification” clusters, which were selected based on their optical lensing properties. CLASH aims to study the dark matter distribution of the clusters and find magnified high-redshift galaxies behind them. CLASH observations were carried out in 16 bands from UV to NIR to derive accurate and reliable estimates of photometric redshifts. Aims. We present the CLASH photometric redshifts using 16 HST bands and study the photometric redshift accuracy including a detailed comparison between photometric and spectroscopic redshifts for the strong lensing arcs using the measurements from the cluster MACSJ1206.2-0847. Methods. We used the publicly available Le Phare and BPZ photometric redshift estimation codes on 17 CLASH galaxy clusters for which the full photo-z data processing had been completed at the time of this analysis, and derive an estimate of the CLASH photo-z accuracy. Results. Using Le Phare code for objects with a S/N ≥ 10, we reach a precision of 3%(1 + z) for the strong lensing arcs, which is reduced to 2.4%(1 + z) after removing outliers. For galaxies in the cluster field, the corresponding values are 4%(1 + z) and 3%(1 + z). Using mock galaxy catalogues, we show that 3%(1 + z) precision is what is expected using the baseline sky substraction algorithm when taking into account extinction from dust, emission lines, and the finite range of SEDs included in the photo-z template library. An improved method for estimating galaxy colours that yields more accurate photometric redshifts will be explored in a forthcoming paper. We study photo-z results for different aperture photometry techniques and find that the SExtractor isophotal photometry works best. We check the robustness of the arcs photo-z results by rederiving the input photometry in the case of MACS1206. We describe and release a photometric redshift catalogue of the MACS1206 cluster we study. Conclusions. Our photo-z codes give similar results for the strong lensing arcs, as well as for galaxies of the cluster field. Results are improved when optimizing the photometric aperture shape that shows an optimal aperture size around 1′′ radius, giving results that are equivalent to isophotal photometry. Tailored photometry of the arcs improves the photo-z results by showing more consistency between the different arcs of the same strong lensing system

Degradation analysis in the estimation of photometric redshifts from non-representative training sets

We perform an analysis of photometric redshifts estimated by using a non-representative training sets in magnitude space. We use the ANNz2 and GPz algorithms to estimate the photometric redshift both in simulations as well as in real data from the Sloan Digital Sky Survey (DR12). We show that for the representative case, the results obtained by using both algorithms have the same quality, either using magnitudes or colours as input. In order to reduce the errors when estimating the redshifts with a non-representative training set, we perform the training in colour space. We estimate the quality of our results by using a mock catalogue which is split samples cuts in the $r$-band between $19.4< r< 20.8$. We obtain slightly better results with GPz on single point z-phot estimates in the complete training set case, however the photometric redshifts estimated with ANNz2 algorithm allows us to obtain mildly better results in deeper $r$-band cuts when estimating the full redshift distribution of the sample in the incomplete training set case. By using a cumulative distribution function and a Monte-Carlo process, we manage to define a photometric estimator which fits well the spectroscopic distribution of galaxies in the mock testing set, but with a larger scatter. To complete this work, we perform an analysis of the impact on the detection of clusters via density of galaxies in a field by using the photometric redshifts obtained with a non-representative training set.Comment: 19 pages, 9 figures. Accepted for publication in MNRA

Designing Future Dark Energy Space Missions: II. Photometric Redshift of Space Weak Lensing Optimized Survey

Accurate weak-lensing analysis requires not only accurate measurement of galaxy shapes but also precise and unbiased measurement of galaxy redshifts. The photometric redshift technique appears as the only possibility to determine the redshift of the background galaxies used in the weak-lensing analysis. Using the photometric redshift quality, simple shape measurement requirements, and a proper sky model, we explore what could be an optimal weak-lensing dark energy mission based on FoM calculation. We found that photometric redshifts reach their best accuracy for the bulk of the faint galaxy population when filters have a resolution R~3.2. We show that an optimal mission would survey the sky through 8 filters using 2 cameras (visible and near infrared). Assuming a 5-year mission duration, a mirror size of 1.5m, a 0.5deg2 FOV with a visible pixel scale of 0.15", we found that a homogeneous survey reaching IAB=25.6 (10sigma) with a sky coverage of ~11000deg2 maximizes the Weak Lensing FoM. The effective number density of galaxies then used for WL is ~45gal/arcmin2, at least a factor of two better than ground based survey. This work demonstrates that a full account of the observational strategy is required to properly optimize the instrument parameters to maximize the FoM of the future weak-lensing space dark energy mission.Comment: 25 pages, 39 figures, accepted in A&

Comparison of the properties of two fossil groups of galaxies with the normal group NGC 6034 based on multiband imaging and optical spectroscopy

We collected multiband imaging and spectroscopy for two fossil groups (RX J1119.7+2126 and 1RXS J235814.4+150524) and one normal group (NGC 6034). We computed photometric redshifts in the central zones of each group, combining previous data with the SDSS five-band data. For each group we investigated the red sequence (RS) of the color-magnitude relation and computed the luminosity functions, stellar population ages and distributions of the group members. Spectroscopy allowed us to investigate the large-scale surroundings of these groups and the substructure levels in 1RXS J235814.4+150524 and NGC 6034. The large-scale environment of 1RXS J235814.4+150524 is poor, though its galaxy density map shows a clear signature of the surrounding cosmic web. RX J1119.7+2126 appears to be very isolated, while the cosmic environment of NGC 6034 is very rich. At the group scale, 1RXS J235814.4+150524 shows no substructure. Galaxies with recent stellar populations seem preferentially located in the group outskirts. A RS is discernable for all three groups in a color-magnitude diagram. The luminosity functions based on photometric redshift selection and on statistical background subtraction have comparable shapes, and agree with the few points obtained from spectroscopic redshifts. These luminosity functions show the expected dip between first and second brightest galaxies for the fossil groups only. Their shape is also regular and relatively flat at faint magnitudes down to the completeness level for RX J1119.7+2126 and NGC 6034, while there is a clear lack of faint galaxies for 1RXS J235814.4+150524. RX J1119.7+2126 is definitely classified as a fossil group; 1RXS J235814.4+150524 also has properties very close to those of a fossil group, while we confirm that NGC 6034 is a normal group.Comment: Accepted in A&A, english-improved, 5 jpeg figures, and shortened abstrac