Over the last few years, needlets have a emerged as a useful tool for the analysis of Cosmic Microwave Background (CMB) data. Our aim in this paper is first to introduce in the CMB literature a different form of needlets, known as Mexican needlets, first discussed in the mathematical literature by Geller and Mayeli (2009a,b). We then proceed with an extensive study of the properties of both standard and Mexican needlets; these properties depend on some parameters which can be tuned in order to optimize the performance for a given application. Our second aim in this paper is then to give practical advice on how to adjust these parameters in order to achieve the best properties for a given problem in CMB data analysis. In particular we investigate localization properties in real and harmonic spaces and propose a recipe on how to quantify the influence of galactic and point source masks on the needlet coefficients. We also show that for certain parameter values, the Mexican needlets provide a close approximation to the Spherical Mexican Hat Wavelets (whence their name), with some advantages concerning their numerical implementation and the derivation of their statistical properties.

We use the full bispectrum of spherical needlets applied to the WMAP data of the cosmic microwave background as an estimator for the primordial non-Gaussianity parameter f_NL. We use needlet scales up to l_max=1000 and the KQ75 galactic cut and find f_NL=84 +/- 40 corrected for point source bias. We also introduce a set of consistency tests to validate our results against the possible influence of foreground residuals or systematic errors. In particular, fluctuations in the value of f_NL obtained from different frequency channels, different masks and different multipoles are tested against simulated maps. All variations in f_NL estimates are found statistically consistent with simulations.