Description Usage Arguments Value Author(s) References See Also Examples
Calculates Sersic parameters
1 | sersic(mag, re, n, e = 0, r = re)
|
mag |
total magnitude |
re |
half-light radius |
n |
Sersic index |
e |
ellipticity (1 - b/a) |
r |
radius of interest |
mag |
magnitude within r |
magdiff |
difference between total magnitude and magnitude within r |
mu |
surface brightness at r |
muavg |
average surface brightness within r |
inten |
intensity at r |
lum |
luminosity within r |
lumtot |
total luminosity |
lumfrac |
fraction of total luminosity contained within r |
Lee Kelvin, Aaron Robotham
Maintainer: Lee Kelvin <lee.kelvin@uibk.ac.at>
Graham A. W., Driver S. P., 2005, PASA, 22, 118
The astronomy package: astro
.
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 | r = seq(0,8,len=500)
ns = c(0.5,1,2,4,10)
col = hsv(seq(2/3,0,len=length(ns)))
layout(c(1,2))
par("mar"=c(0,4.1,0,2.1))
par("oma"=c(5.1,0,2.1,0))
# surface brightness plot
mu = convmu(r=r, mu=20, re=1, n=ns[1])
aplot(r,abs(mu), ylim=c(28,12), type="l", nxmaj=5, nxmin=1, nymin=4, xlab="",
ylab=bquote(paste(mu," / mag ", arcsec^{-2})), xaxs="i", yaxs="i", las=1,
col=col[1], labels=2)
for(i in 2:length(ns)){
mu = convmu(r=r, mu=20, re=1, n=ns[i])
lines(r, abs(mu), col=col[i])
}
label("top", txt=bquote(paste(mu["e"]," = 20")), bty="n", inset=0.3)
text(0.9,17,labels="n = 10")
text(2,25,labels="n = 0.5")
legend("topright", legend=c(0.5,1,2,4,10), col=col, bty="n", lty=1, inset=0.03)
box()
# magnitude difference plot
md = sersic(mag=0, re=1, n=ns[1], r=r)$magdiff
aplot(r,abs(md), ylim=c(0,3), type="l", nxmaj=5, nxmin=1, nymaj=4, nymin=4,
xlab="", ylab="m (< r) / mag", xaxs="i", yaxs="i", las=1, col=col[1])
for(i in 2:length(ns)){
md = sersic(mag=0, re=1, n=ns[i], r=r)$magdiff
lines(r, abs(md), col=col[i])
}
abline(h=abs(sersic(mag=0, re=1, n=1, r=1)$magdiff), lty=2)
label("top", txt=bquote(paste(m[tot]," = 0")), bty="n", inset=0.3)
text(7,0.4,labels="n = 10")
text(1,2,labels="n = 0.5")
mtext(bquote(paste("r / ",r[e],sep="")), side=1, line=2.5)
box()
|
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.