Researcher profile

J. L. Sanz

J. L. Sanz contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
18works
0followers
3topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

18 published item(s)

preprint2012arXiv

Extragalactic point source detection in WMAP 7-year data at 61 and 94 GHz

The detection of point sources in Cosmic Microwave Background maps is usually based on a single-frequency approach, whereby maps at each frequency are filtered separately and the spectral information on the sources is derived combining the results at the different frequencies. On the contrary, in the case of multi-frequency detection methods, source detection and spectral information are tightly interconnected in order to increase the source detection efficiency. In this work we apply the \emph{matched multifilter} method to the detection of point sources in the WMAP 7yr data at 61 and 94 GHz. This linear filtering technique takes into account the spatial and the cross-power spectrum information at the same time using the spectral behaviour of the sources without making any a priori assumption about it. We follow a two-step approach. First, we do a blind detection of the sources over the whole sky. Second, we do a refined local analysis at their positions to improve the signal-to-noise ratio of the detections. At 94 GHz we detect 129 5σ objects at |b|>5° (excluding the Large Magellanic Cloud region); 119 of them are reliable extragalactic sources and 104 of these 119 lie outside the WMAP Point Source Catalog mask. Nine of the total 129 detections are known Galactic sources or lie in regions of intense Galactic emission and one additional (weak) high-Galactic latitude source has no counterpart in low-frequency radio catalogues. Our results constitute a substantial improvement over the NEWPS-3year catalogue.

preprint2011arXiv

A Bayesian technique for the detection of point sources in CMB maps

The detection and flux estimation of point sources in cosmic microwave background (CMB) maps is a very important task in order to clean the maps and also to obtain relevant astrophysical information. In this paper we propose a maximum a posteriori (MAP) approach detection method in a Bayesian scheme which incorporates prior information about the source flux distribution, the locations and the number of sources. We apply this method to CMB simulations with the characteristics of the Planck satellite channels at 30, 44, 70 and 100 GHz. With a similar level of spurious sources, our method yields more complete catalogues than the matched filter with a 5 sigma threshold. Besides, the new technique allows us to fix the number of detected sources in a non-arbitrary way.

preprint2011arXiv

Adaptive Langevin Sampler for Separation of t-Distribution Modelled Astrophysical Maps

We propose to model the image differentials of astrophysical source maps by Student's t-distribution and to use them in the Bayesian source separation method as priors. We introduce an efficient Markov Chain Monte Carlo (MCMC) sampling scheme to unmix the astrophysical sources and describe the derivation details. In this scheme, we use the Langevin stochastic equation for transitions, which enables parallel drawing of random samples from the posterior, and reduces the computation time significantly (by two orders of magnitude). In addition, Student's t-distribution parameters are updated throughout the iterations. The results on astrophysical source separation are assessed with two performance criteria defined in the pixel and the frequency domains.

preprint2011arXiv

Filter Design for the Detection/Estimation of the Modulus of a Vector. Application to Polarization Data

We consider a set of M images, whose pixel intensities at a common point can be treated as the components of a M-dimensional vector. We are interested in the estimation of the modulus of such a vector associated to a compact source. For instance, the detection/estimation of the polarized signal of compact sources immersed in a noisy background is relevant in some fields like Astrophysics. We develop two different techniques, one based on the Maximum Likelihood Estimator (MLE) applied to the modulus distribution, the modulus filter (ModF) and other based on prefiltering the components before fusion, the filtered fusion (FF), to deal with this problem. We present both methods in the general case of M images and apply them to the particular case of three images (linear plus circular polarization). Numerical simulations have been performed to test these filters considering polarized compact sources immersed in stationary noise. The FF performs better than the ModF in terms of errors in the estimated amplitude and position of the source, especially in the low signal-to-noise case. We also compare both methods with the direct application of a matched filter (MF) on the polarization data. This last technique is clearly outperformed by the new methods.

preprint2011arXiv

Joint Bayesian separation and restoration of CMB from convolutional mixtures

We propose a Bayesian approach to joint source separation and restoration for astrophysical diffuse sources. We constitute a prior statistical model for the source images by using their gradient maps. We assume a t-distribution for the gradient maps in different directions, because it is able to fit both smooth and sparse data. A Monte Carlo technique, called Langevin sampler, is used to estimate the source images and all the model parameters are estimated by using deterministic techniques.

preprint2010arXiv

Constraints on f_nl and g_nl from the analysis of the N-pdf of the CMB large scale anisotropies

[Abridged] In this paper we explore a local non-linear perturbative model up to third order as a general characterization of the CMB anisotropies. We focus our analysis in large scale anisotropies. At these angular scales, the non-Gaussian description proposed in this work defaults (under certain conditions) to an approximated local form of the weak non-linear coupling inflationary model. In particular, quadratic and cubic terms are governed by the non-linear coupling parameters f_nl and g_nl, respectively. The extension proposed in this paper allows us to directly constrain these non-linear parameters. Applying the proposed methodology to WMAP 5-yr data, we obtain -5.6 x 10^5 < g_nl < 6.4 x 10^5, at 95% CL. This result is in agreement with previous findings obtained for equivalent non-Gaussian models and with different non-Gaussian estimators. A model selection test is performed, indicating that a Gaussian model is preferred to the non-Gaussian scenario. When comparing different non-Gaussian models, we observe that a pure f_nl model is the most favoured case, and that a pure g_nl model is more likely than a general non-Gaussian scenario. Finally, we have analyzed the WMAP data in two independent hemispheres, in particular the ones defined by the dipolar pattern found by Hoftuft et al. 2009. We show that, whereas g_nl is still compatible with zero for both hemispheres, it is not the case for f_nl (with a p-value 0.04). However, if anisotropy of the data is assumed, the distance between the likelihood distributions for each hemisphere is larger than expected from Gaussian and anisotropic simulations, also for g_nl (with a p-value of 0.001 in the case of this parameter). This result is an extra evidence for the CMB asymmetries previously reported in WMAP data.

preprint2009arXiv

A multifrequency method based on the Matched Multifilter for the detection of point sources in CMB maps

In this work we deal with the problem of simultaneous multifrequency detection of extragalactic point sources in maps of the Cosmic Microwave Background. We apply a linear filtering technique that uses spatial information and the cross-power spectrum. To make this, we simulate realistic and non-realistic flat patches of the sky at two frequencies of Planck: 44 and 100 GHz. We filter to detect and estimate the point sources and compare this technique with the monofrequency matched filter in terms of completeness, reliability, flux and spectral index accuracy. The multifrequency method outperforms the matched filter at the two frequencies and in all the studied cases in the work.

preprint2009arXiv

Analysis of non-Gaussian CMB maps based on the N-pdf. Application to WMAP data

We present a new method based on the N-point probability distribution (pdf) to study non-Gaussianity in cosmic microwave background (CMB) maps. Likelihood and Bayesian estimation are applied to a local non-linear perturbed model up to third order, characterized by a linear term which is described by a Gaussian N-pdf, and a second and third order terms which are proportional to the square and the cube of the linear one. We also explore a set of model selection techniques (the Akaike and the Bayesian Information Criteria, the minimum description length, the Bayesian Evidence and the Generalized Likelihood Ratio Test) and their application to decide whether a given data set is better described by the proposed local non-Gaussian model, rather than by the standard Gaussian temperature distribution. As an application, we consider the analysis of the WMAP 5-year data at a resolution of around 2 degrees. At this angular scale (the Sachs-Wolfe regime), the non-Gaussian description proposed in this work defaults (under certain conditions) to an approximative local form of the weak non-linear coupling inflationary model (e.g. Komatsu & Spergel 2001) previously addressed in the literature. For this particular case, we obtain an estimation for the non-linear coupling parameter of -94 < F_nl < 154 at 95% CL. Equally, model selection criteria also indicate that the Gaussian hypothesis is favored against the particular local non-Gaussian model proposed in this work. This result is in agreement with previous findings obtained for equivalent non-Gaussian models and with different non-Gaussian estimators. However, our estimator based on the N-pdf is more efficient than previous estimators and, therefore, provides tighter constraints on the coupling parameter at degree angular resolution.

preprint2009arXiv

Detection/estimation of the modulus of a vector. Application to point source detection in polarization data

Given a set of images, whose pixel values can be considered as the components of a vector, it is interesting to estimate the modulus of such a vector in some localised areas corresponding to a compact signal. For instance, the detection/estimation of a polarized signal in compact sources immersed in a background is relevant in some fields like astrophysics. We develop two different techniques, one based on the Neyman-Pearson lemma, the Neyman-Pearson filter (NPF), and another based on prefiltering-before-fusion, the filtered fusion (FF), to deal with the problem of detection of the source and estimation of the polarization given two or three images corresponding to the different components of polarization (two for linear polarization, three including circular polarization). For the case of linear polarization, we have performed numerical simulations on two-dimensional patches to test these filters following two different approaches (a blind and a non-blind detection), considering extragalactic point sources immersed in cosmic microwave background (CMB) and non-stationary noise with the conditions of the 70 GHz \emph{Planck} channel. The FF outperforms the NPF, especially for low fluxes. We can detect with the FF extragalactic sources in a high noise zone with fluxes >= (0.42,0.36) Jy for (blind/non-blind) detection and in a low noise zone with fluxes >= (0.22,0.18) Jy for (blind/non-blind) detection with low errors in the estimated flux and position.

preprint2009arXiv

Polarization of the WMAP Point Sources

The detection of polarized sources in the WMAP 5-year data is a very difficult task. The maps are dominated by instrumental noise and only a handful of sources show up as clear peaks in the Q and U maps. Optimal linear filters applied at the position of known bright sources detect with a high level of significance a polarized flux P from many more sources, but estimates of P are liable to biases. Using a new technique, named the &#34;filtered fusion technique&#34;, we have detected in polarization, with a significance level greater than 99.99% in at least one WMAP channel, 22 objects, 5 of which, however, do not have a plausible low radio frequency counterpart and are therefore doubtful. Estimated polarized fluxes P < 400 mJy at 23 GHz were found to be severely affected by the Eddington bias. The corresponding polarized flux limit for Planck/LFI at 30 GHz, obtained via realistic simulations, is 300 mJy. We have also obtained statistical estimates of, or upper limits to the mean polarization degrees of bright WMAP sources at 23, 33, 41, and 61 GHz, finding that they are of a few percent.

preprint2003arXiv

Goodness-of-fit tests of Gaussianity: constraints on the cumulants of the MAXIMA data

In this work, goodness-of-fit tests are adapted and applied to CMB maps to detect possible non-Gaussianity. We use Shapiro-Francia test and two Smooth goodness-of-fit tests: one developed by Rayner and Best and another one developed by Thomas and Pierce. The Smooth tests test small and smooth deviations of a prefixed probability function (in our case this is the univariate Gaussian). Also, the Rayner and Best test informs us of the kind of non-Gaussianity we have: excess of skewness, of kurtosis, and so on. These tests are optimal when the data are independent. We simulate and analyse non-Gaussian signals in order to study the power of these tests. These non-Gaussian simulations are constructed using the Edgeworth expansion, and assuming pixel-to-pixel independence. As an application, we test the Gaussianity of the MAXIMA data. Results indicate that the MAXIMA data are compatible with Gaussianity. Finally, the values of the skewness and kurtosis of MAXIMA data are constrained by |S| \le 0.035 and |K| \le 0.036 at the 99% confidence level.

preprint2003arXiv

Tests of Gaussianity

We review two powerful methods to test the Gaussianity of the cosmic microwave background (CMB): one based on the distribution of spherical wavelet coefficients and the other on smooth tests of goodness-of-fit. The spherical wavelet families proposed to analyse the CMB are the Haar and the Mexican Hat ones. The latter is preferred for detecting non-Gaussian homogeneous and isotropic primordial models containing some amount of skewness or kurtosis. Smooth tests of goodness-of-fit have recently been introduced in the field showing some interesting properties. We will discuss the smooth tests of goodness-of-fit developed by Rayner and Best for the univariate as well as for the multivariate analysis.

preprint2001arXiv

Spherical Mexican Hat wavelet: an application to detect non-Gaussianity in the COBE-DMR maps

The spherical Mexican Hat wavelet is introduced in this paper, with the aim of testing the Gaussianity of the Cosmic Microwave Background temperature fluctuations. Using the information given by the wavelet coefficients at several scales, we have performed several statistical tests on the COBE-DMR maps to search for evidence of non-Gaussianity. Skewness, kurtosis, scale-scale correlations (for two and three scales) and Kolmogorov-Smirnov tests indicate that the COBE-DMR data are consistent with a Gaussian distribution. We have extended the analysis to compare temperature values provided by COBE-DMR data with distributions (obtained from Gaussian simulations) at each pixel and at each scale. The number of pixels with temperature values outside the 95% and the 99% is consistent with that obtained from Gaussian simulations, at all scales. Moreover, the extrema values for COBE-DMR data (maximum and minimum temperatures in the map) are also consistent with those obtained from Gaussian simulations.

preprint2000arXiv

Geometrical estimators as a test of Gaussianity in the CMB

We investigate the power of geometrical estimators on detecting non-Gaussianity in the cosmic microwave background. In particular the number, eccentricity and Gaussian curvature of excursion sets above (and below) a threshold are studied. We compare their different performance when applied to non-Gaussian simulated maps of small patches of the sky, which take into account the angular resolution and instrumental noise of the Planck satellite. These non-Gaussian simulations are obtained as perturbations of a Gaussian field in two different ways which introduce a small level of skewness or kurtosis in the distribution. A comparison with a classical estimator, the genus, is also shown. We find that the Gaussian curvature is the best of our estimators in all the considered cases. Therefore we propose the use of this quantity as a particularly useful test to look for non-Gaussianity in the CMB.

preprint2000arXiv

Tests of Gaussianity of CMB maps

We study two different methods to test Gaussianity in CMB maps. One of them is based on the partition function and the other on the morphology of hot and cold spots. The partition function contains information on all the moments and scales, being a useful quantity to compress the large data sets expected from future space missions like Planck. In particular, it contains much richer information than the one available through the radiation power spectrum. The second method utilizes morphological properties of hot and cold spots such as the eccentricity and number of spots in CMB maps. We study the performance of both methods in detecting non-Gaussian features in small scale CMB simulated maps as those which will be provided by the Planck mission.

preprint1997arXiv

Weak lensing correlations in open and flat universes

Correlations between the magnification or polarization of background sources, induced by gravitational lensing due to the large-scale structure, and the positions of foreground galaxies are investigated. We found that their amplitude is enhanced with respect to correlations for a single population. We analize the dependence of the correlations with the density parameter Ωconsidering a nonlinear evolution of the matter power spectrum. The contribution of the linear evolution is negligible at scales below several arcminutes. Detailed results on the dependence of the correlations on the redshift of the foreground and background populations for different cosmological models are presented. The effect reaches its maximum amplitude for foreground populations with relatively small redshifts due to the fast increase of the nonlinear matter power spectrum at recent times.

preprint1995arXiv

Genus and spot density in the COBE DMR first year anisotropy maps

A statistical analysis of texture on the {\it COBE}-DMR first year sky maps based on the genus and spot number is presented. A generalized $χ^2$ statistic is defined in terms of ``observable&#39;&#39; quantities: the genus and spot density that would be measured by different cosmic observers. This strategy together with the use of Monte Carlo simulations of the temperature fluctuations, including all the relevant experimental parameters, represent the main difference with previous analyses. Based on the genus analysis we find a strong anticorrelation between the quadrupole amplitude $Q_{rms-PS}$ and the spectral index $n$ of the density fluctuation power spectrum at recombination of the form $Q_{rms-PS}= 22.2 \pm 1.7 - (4.7 \pm 1.3) \times n\ μ$K for fixed $n$, consistent with previous works. The result obtained based on the spot density is consistent with this $Q_{rms-PS} (n)$ relation. In addition to the previous results we have determined, using Monte Carlo simulations, the minimum uncertainty due to cosmic variance for the determination of the spectral index with the genus analysis. This uncertainty is $δn\approx 0.2$.

preprint1995arXiv

Omega from the COBE-DMR anisotropy maps

We have made a likelihood statistical analysis of the angular correlations in the {\it COBE}-DMR two-year sky maps by Monte Carlo simulation of the temperature fluctuations. We assume an open universe and consider as primordial power spectrum the Harrison-Zeldovich one, $P(k)=Ak$. We find that the flatness of the universe is not implied by the data. The quadrupole normalization amplitude, $Q_{rms-PS}$, is related to the density parameter, $Ω$, by $Q_{rms-PS} = 10.67 + 55.81 Ω- 128.59 Ω^2 + 81.26 Ω^3\ μ$K. We have determined the p.d.f. of $Ω$ due to cosmic plus sampling (i.e. $20^\circ$ galactic cut) variance which generically shows a bimodal shape. The uncertainty as given by the r.m.s. is $\approx 0.35$, therefore to better constrain $Ω$ experiments sensitive to higher multipoles ($l>20$) should be considered.