Project

General

Profile

Support #3420

Import of galaxy SED

Added by Ben Hambrecht over 6 years ago. Updated over 6 years ago.

Status:
Assigned
Priority:
High
Assignee:
Start date:
02/04/2013
Due date:
% Done:

0%

Estimated time:
Duration:

Description

Hi Carlos,

care to comment:

# to be used for converting spectral templates for Sim_Gal_observed

def import_convert_galaxy_sed_spectra_variants():

    #;read in all the variants on the spectra
    tspecfile    =  'k_nmf_derived.default.fits'

    hdr        = headfits(tspecfile)
    nt        = long(sxpar(hdr,'NT'))
    loglam        = np.log10(mrdfits(tspecfile, 11))
    tspec_v0    = mrdfits(tspecfile, 1)
    tspec_v0_nl    = mrdfits(tspecfile, 2)
    tspec_v0_nd    = mrdfits(tspecfile, 3)
    tspec_v0_nd_nl    = mrdfits(tspecfile, 4)
    tspec_v300    = mrdfits(tspecfile, 5)
    tspec_v300_nl    = mrdfits(tspecfile, 6)
    tspec_v300_nd    = mrdfits(tspecfile, 7)
    tspec_v300_nd_nl= mrdfits(tspecfile, 8)
    lspec_v300    = mrdfits(tspecfile, 9)
    tmass        = mrdfits(tspecfile, 17)
    tmetallicity    = mrdfits(tspecfile, 18)
    tmass300    = mrdfits(tspecfile, 19)
    tmass1000    = mrdfits(tspecfile, 20)
    tmremain    = mrdfits(tspecfile, 24)

    with archive.archive(file_data_bank,'a') as data_bank:
        data_bank['/galaxy/flux_spectrum_observed']             = Vfluxobs
        data_bank['/galaxy/flux_spectrum_observed_noise_free']  = Vfluxobs_noise_free
        data_bank['/galaxy/wavelength']                         =  some wavelength
        data_bank['/galaxy/hdr'] = hdr
        data_bank['/'] =  nt
        data_bank['/'] =  loglam
        data_bank['/'] =  tspec_v0              # not smoothed
        data_bank['/'] =  tspec_v0_nl
        data_bank['/'] =  tspec_v0_nd
        data_bank['/'] =  tspec_v0_nd_nl
        data_bank['/'] =  tspec_v300            # smoothed to 300kms(res)
        data_bank['/'] =  tspec_v300_nl
        data_bank['/'] =  tspec_v300_nd
        data_bank['/'] =  tspec_v300_nd_nl
        data_bank['/'] =  lspec_v300
        data_bank['/'] =  tmass                 # stellar stuff
        data_bank['/'] =  tmetallicity
        data_bank['/'] =  tmass300
        data_bank['/'] =  tmass1000
        data_bank['/'] =  tmremain

This seems to import some template spectra, from which module 5 will simulate spectra for all targeted galaxies.
As Measure_Redshift will use the same templates, it would be nice to have written down what all these quantities are.

Cheers
Ben

History

#1 Updated by Carlos Cunha over 6 years ago

I didn't make this, but the labelling seems off. Here's the correct labels, from table 3 of this paper http://adsabs.harvard.edu/abs/2007AJ....133..734B

Number Name Dimensions Description
0 templates Nbasis × Nt Coefficients in basis space for each template
1 spec Nspec × Nt Smoothed spectrum of each template
2 spec nl Nspec × Nt Smoothed spectrum of each template, lines removed
3 spec nd Nspec × Nt Smoothed spectrum of each template, without dust
4 spec nl nd Nspec × Nt Smoothed spectrum of each template, without lines or dust
5 rawspec Nspec × Nt Unsmoothed spectrum of each template
6 rawspec nl Nspec × Nt Unsmoothed spectrum of each template, lines removed
7 rawspec nd Nspec × Nt Unsmoothed spectrum of each template, without dust
8 rawspec nl nd Nspec × Nt Unsmoothed spectrum of each template, without lines or dust
9 lspec Nspec × Nt emission line spectrum
10 extinction Nspec × Nt extinction as a function of wavelength
11 lambda Nspec Wavelength grid (A)
12 sfr Nage × Nt Star-formation rate (M⊙ per year) as a function of time
13 metallicity Nage × Nt Average metallicity as a function of time
14 ages Nage Age grid (yrs)
15 dage Nage Age bin size (yrs); sfr×dage = total star-formation in each bin
16 mass Nt total mass formed in each template
17 mremain Nt total current stellar mass in each template
18 mets Nt average metallicity of current stars in each template
19 m300 Nt total stellar mass formed in last 300 Myrs in each template
20 m1000 Nt total stellar mass formed in last 1 Gyr in each template
21 basis ages Nage × Nmets × Ndust ages of each stellar population basis vector
22 basis mets Nage × Nmets × Ndust metallicity of each stellar population basis vector
23 basis dusts Nage × Nmets × Ndust dust properties of each stellar population basis vector
24 basis mremain Nage × Nmets × Ndust fraction of original stellar mass surviving for each basis vector

#2 Updated by Brian Nord over 6 years ago

  • Tracker changed from Bug to Support
  • Status changed from New to Assigned
  • Assignee changed from Carlos Cunha to Brian Nord
  • Priority changed from Normal to High

Need to check these and put link in the file for fast checking later.



Also available in: Atom PDF