plot.ProSpectSED: Plot ProSpect SED Outputs

View source: R/ProSpect.R

plot.ProSpectR Documentation

Plot ProSpect SED Outputs

Description

Simple plot of the various luminosity or flux components in a ProSpectSED/ProSpectSEDlike object.

Usage

## S3 method for class 'ProSpectSED'
plot(x, xlim = c(1e3,1e7), ylim = 'auto',
  xlab = 'Wavelength (Ang)', ylab = 'auto', grid = TRUE, type = 'lum', lwd_main = 5,
  lwd_comp = 5, ...)
## S3 method for class 'ProSpectSEDlike'
plot(x, xlim=c(1e3,1e7), ylim='auto',
xlab='Wavelength (Ang)', ylab='auto', grid=TRUE, type='flux', ...)

Arguments

x

Object of class 'ProSpectSED'.

xlim

Numeric vector; length 1/2, x limits.

ylim

Numeric vector; length 1/2, y limits. Default 'auto' generally chooses a sensible range.

xlab

Character scalar; x label.

ylab

Character scalar; y label. Default 'auto' generally chooses a sensible name.

grid

Logical; show background grid?

type

Character scalar; specifies whether to plot luminosity density (type='lum') or flux density (type='flux').

lwd_main

Numeric scalar; width for the main spectral line.

lwd_comp

Numeric scalar; width for the sub-component spectral line.

...

Passed to magplot or plot.

Details

A simple plot for quick visualisation of ProSpectSED or ProSpectSEDlike.

Value

Side effect is to create the plot.

Author(s)

Aaron Robotham

See Also

ProSpectSED, ProSpectSEDlike

Examples

data('BC03lr')
data('Dale_NormTot')
data('AGN_UnOb_Sparse')
data('Dale_M2L_func')
filters=c('FUV_GALEX', 'NUV_GALEX', 'u_SDSS', 'g_SDSS', 'r_SDSS', 'i_SDSS', 'Z_VISTA',
'Y_VISTA', 'J_VISTA', 'H_VISTA', 'K_VISTA', 'W1_WISE' , 'W2_WISE', 'W3_WISE', 'W4_WISE',
'P100_Herschel', 'P160_Herschel', 'S250_Herschel' , 'S350_Herschel', 'S500_Herschel')
filtout={}
for(i in filters){filtout=c(filtout,list(approxfun(getfilt(i))))}

testSED=ProSpectSED(AGNlum=1e43, tau_birth=1, tau_screen=0.3, tau_AGN=2,
alpha_SF_birth=1, alpha_SF_screen=3, alpha_SF_AGN=0, speclib=BC03lr, Dale=Dale_NormTot,
AGN=AGN_UnOb_Sparse, filtout=filtout, Dale_M2L_func=Dale_M2L_func, returnall=TRUE)
plot(testSED)

#An example of a more physical SFH and ZH:

testSED2 = ProSpectSED(massfunc = massfunc_snorm, Z=Zfunc_massmap_box)
plot(testSED2,type = 'lum')

asgr/ProSpect documentation built on Feb. 21, 2025, 1:43 a.m.