Software for fitting and classifying IRS galaxy spectra with the NMF components described in Hurley et al. 2013
An example idl script (example_script.pro) is provided to demonstrate the fitting and classifying procedures.
To run the code:
-
Download into a suitable directory within your IDL path
-
Download the bvls.pro script from http://www-astro.physics.ox.ac.uk/~mxc/idl/bvls.pro
-
Start idl within the NMF_software directory and run example_script.pro
-
IDL should output a postscript plot (example_IRAS_10378+1108fit.ps) of the NMF fit and give the probability of being in the Sy1, Sbrst, Sy2, Sy1.5 or Dusty cluster
-
The NMF components and wavelength array are contained within NMF_components.sav (for use within idl), as well as NMF_components.ascii for use outside idl.
Feel free to edit/use code and NMF_components and classifier, but please cite Hurley et al. 2013.
GMM_NMF_classifier.pro: Calculates the probability of belonging to one of the 5 clusters
GMMclusters.sav: IDL structure containing normalisation, centres and covariance matrices to decribe each cluster
IRAS_10378+1108.dat: Example IRS spectrum
NMF7_fit.pro: Main code to fit the 7 NMF componentns to a galaxy spectrum (requires bvls.pro)
NMF_components.ascii: Wavelength and flux density of NMF components in ascii form for use outside IDL
NMF_components.sav: Wavelength and flux density of NMF components in IDL's native '.sav' format.
ajs_gaussnd.pro: Calculates the pdf value for each of the five clusters, at any point in NMF space. (uses GMMclusters.sav)
closeps.pro: Useful bit of IDL code to to tidy up postscript device within IDL when plotting
example_script.pro: Example IDL script to demonstrate use of NMF components and classifier
my_chisq.pro: Calculates chi square nearest.pro: Does nearest neighbour interpolation