galaxies accurately, we add UV and NIR photometric data points to the optical data. As Ma et al. (2009b) discussed, Kaviraj et al. (2007) showed that the combination of FUV and NUV photometry with optical observations in the standard broad bands enables one to efficiently break the age-metallicity degeneracy. The optical broadband colors have the more obvious problem of age-metallicity degeneracy (Worthey 1994; MacArthur et al. 2004). Again, de Jong (1996) showed that such degeneracy can be partially broken by adding NIR photometry to optical colors, which has also been stated by Bell et al. (2000), Wu et al. (2005) and Li et al. (2007).

Our observational data consist of integrated luminosities through a given set of filters (FUV, NUV, ugriz, JHKs), thus we convolved the theoretical SEDs with these filter response curves to obtain synthetic ultraviolet, optical and NIR photometry for comparison. The synthetic magnitude in the AB magnitude system for the $i$th filter can be computed as $$m_i=-2.5{\rm log}\frac{\int_{\lambda} F_\lambda\varphi_i(\lambda)d\lambda}{\int_{\lambda}\varphi_i(\lambda)d\lambda}-48.60,$$ where $F_{\lambda}$ is the theoretical SED and $\varphi_i$ is the response curve of the $i$th filter of the sets of filter we used.

We use a $\chi^2$ minimization test to examine which model SEDs are most comparable with 