phy-4660/xrd/plots.R

x = seq(from = 5.00547782,to = 119.99571965,length.out = 13763)

grease_count = read.table("grease.tab")
plastic_count = read.table("plastic.tab")
wood_count = read.table("wood.tab")

plot(x=x,y=plastic_count[,1],pch='-',log="x",ylab=expression(count),xlab=expression(2~theta~degree),xlim=c(12,55))
abline(v=42.5)
axis(3,at=42.5)
abline(v=30.5)
axis(3,at=30.5)
abline(v=14)
axis(3,at=14)

dev.print(pdf,"../report/plastic_peaks.pdf",width=3,height=3)




plot(x=x,y=grease_count[,1],pch='-',log="x",ylab=expression(count),xlab=expression(2~theta~degree))
abline(v=12)
axis(3,at=12)
abline(v=22)
axis(3,at=22)
dev.print(pdf,"../report/grease_peaks.pdf",width=3,height=3)


plot(x=x,y=wood_count[,1],pch='-',log="x",ylab=expression(count),xlab=expression(2~theta~degree),xlim=c(15,55))
abline(v=16)
axis(3,at=16)
abline(v=22.5)
axis(3,at=22.5)
abline(v=34.5)
axis(3,at=34.5)
abline(v=46)
axis(3,at=46)
dev.print(pdf,"../report/wood_peaks.pdf",width=3,height=3)