generated from erosita/uds
99 lines
2.6 KiB
R
99 lines
2.6 KiB
R
#
|
|
# 2 data sets
|
|
# in Xspec make ipl->wdata
|
|
# then split file to four parts as lecr1/2/3/4
|
|
#
|
|
|
|
|
|
library(sfsmisc)
|
|
|
|
|
|
up=2.0
|
|
#xgrid <- c(30,50,70,100);
|
|
Emin=25.0
|
|
Emax=180.0
|
|
postscript('cutoffpl_solo_gb.eps', horizontal = FALSE, onefile = FALSE, paper = "special",width = 9.0, height = 6.0)
|
|
par(oma=c(1,1,1,1))
|
|
par(lwd=2)
|
|
|
|
layout(matrix(c(1,2), 2, 1, byrow = TRUE), heights = c(1.5,1), TRUE)
|
|
leftx=5.3
|
|
par(mar=c(0.1, leftx, 1, 1))
|
|
par(cex.lab=1.6)
|
|
par(cex.axis=1.6)
|
|
xlim <- c(Emin,Emax)
|
|
ylim <- c(0.2,2)
|
|
|
|
cl="black"
|
|
a <- read.table('./cutoffpl_solo_gb_eeuf.dat', col.names=c("x","dx","y","dy","cutoffpl"))
|
|
plot(a$x, a$y, pch=3, bg="white", col="white", ylim=ylim, xlim=xlim, ylab=expression("keV"^"2"~"(Phot. keV"^"-1"~"cm"^"-2"~"s"^"-1"*")"),type="p",xaxt = 'n',xlab="",log="xy")
|
|
|
|
#axis(side=2, at=c(1e-6,1e-5, 1e-4, 1e-3, 1e-2), labels=expression('-6','-5','-4','-3','-2'))
|
|
|
|
#abline(v=xgrid, col="lightgray", lty="dotted")
|
|
grid()
|
|
|
|
###################################################################################
|
|
##
|
|
## Plot model components
|
|
##
|
|
###################################################################################
|
|
|
|
lines(a$x,a$cutoffpl,pch=1,bg=cl,col="blue", lwd=3, lty="longdash")
|
|
#lines(a$x,a$pow,pch=1,bg=cl,col="red")
|
|
#lines(a$x,a$total,pch=1,bg=cl,col="black")
|
|
|
|
|
|
#
|
|
# Plot upper limits
|
|
#
|
|
upx=a$x[(a$y/a$dy)<up]
|
|
upy=a$y[(a$y/a$dy)<up]
|
|
updx=a$dx[(a$y/a$dy)<up]
|
|
updy=a$dy[(a$y/a$dy)<up]
|
|
|
|
segments(upx-updx,updy*up,upx+updx,updy*up,col=cl)
|
|
segments(upx,updy*up/2,upx,updy*up,col=cl)
|
|
points(upx,updy*up/2,pch=25,bg=cl,col=cl)
|
|
|
|
px=a$x[(a$y/a$dy)>=up]
|
|
py=a$y[(a$y/a$dy)>=up]
|
|
pdx=a$dx[(a$y/a$dy)>=up]
|
|
pdy=a$dy[(a$y/a$dy)>=up]
|
|
|
|
segments(px,py-pdy,px,py+pdy,col=cl)
|
|
segments(px-pdx,py,px+pdx,py,col=cl)
|
|
|
|
text(100, 1.5, "GB",cex=2.5)
|
|
|
|
###################################################################################
|
|
|
|
par(mar=c(4.5, leftx, 0.1, 1))
|
|
ylim <- c(-4.5,4.5)
|
|
|
|
cl="black"
|
|
a <- read.table("./cutoffpl_solo_gb_delchi.dat", col.names=c("x","dx","y","dy"))
|
|
plot(a$x, a$y, pch=3, ylim=ylim, xlim=xlim, ylab=expression(Delta~chi), xlab="",type="p",log="x",xaxt = 'n')
|
|
grid()
|
|
segments(a$x,a$y-a$dy,a$x,a$y+a$dy,col=cl)
|
|
segments(a$x-a$dx,a$y,a$x+a$dx,a$y,col=cl)
|
|
|
|
#abline(v=xgrid, col="lightgray", lty="dotted")
|
|
#grid()
|
|
#abline(h=0, col = "black",lty=2)
|
|
#dev.off()
|
|
|
|
#abline(v=xgrid, col="lightgray", lty="dotted")
|
|
abline(h=0, col = "black",lty=2)
|
|
|
|
mtext(side = 1, text = "Energy, keV", line = 4, cex=1.6)
|
|
|
|
#axis(1, mgp=c(3, 1.5, 0))
|
|
### Log axis using sfsmisc ###
|
|
atx=c(8,30,50,80,100,150)
|
|
eaxis(1, at = atx, labels = pretty10exp(atx, sub10=c(1,100), drop.1=TRUE), las=0)
|
|
abline(v=atx, col="lightgray", lty="dotted")
|
|
###
|
|
|
|
dev.off()
|