more accel lab stuff

This commit is contained in:
caes 2017-03-27 09:11:20 -04:00
parent f5d4608613
commit c471ec5226
3 changed files with 72 additions and 0 deletions

BIN
accelerator/data/.RData Normal file

Binary file not shown.

View File

@ -0,0 +1,72 @@
data1=read.table("0302171.tab",header=T)
plot(data1)
data2=read.table("0302172.tab",header=T)
lines(data2)
lines(data1)
plot()
lines(data1)
sdfklj
?plot
plot(data1,type='c')
?plot
plot(data1,type='l')
lines(data2)
lines(data2,col='red')
plot(data1,type='l',xlim=c(190,210))
lines(data2,col='red')
plot(data1,type='l',xlim=c(195,205))
lines(data2,col='red')
plot(data1,type='l',xlim=c(95,205))
lines(data2,col='red')
data2
ls()
test1=c(1)
ls()
rm(test1)
ls()
max_cu_ruth = 17200
lines
?lines
lines(c(0,1023),c(17200,17200),col="blue")
lines(c(0,1023),c(17400,17400),col="blue")
lines(c(0,1023),c(17600,17600),col="blue")
lines(c(0,1023),c(17500,17500),col="blue")
plot(data1,type='l',xlim=c(95,205))
lines(c(0,1023),c(17500,17500),col="blue")
plot(data1,type='l',xlim=c(95,205))
lines(c(0,1023),c(17600,17600),col="blue")
lines(data2,col='red')
lines(c(0,1023),c(20000,20000),col="green")
lines(c(0,1023),c(20000,20000),col="purple")
.2*20000
?lines
?plot
lines(c(0,1023),c(4000,4000),col="purple")
.2*17600
lines(c(0,1023),c(3520,3520),col="blue")
grid.locator
identify
identify(x,y)
identify(200,16000)
identify(length,coverage)
identify(channel,count)
identify("channel","count")
identify(channel,count)
identify(data1$channel,count)
identify(data1$channel,data1$count)
plot(data1,type='l',xlim=c(95,205))
lines(c(0,1023),c(17600,17600),col="blue")
lines(c(0,1023),c(3520,3520),col="blue")
plot(data1,type='l',xlim=c(95,205))
lines(c(0,1023),c(17600,17600),col="blue")
lines(data2,col='red')
lines(c(0,1023),c(20000,20000),col="purple")
dev.print(pdf,"../report/rutherford.pdf")
?dev.print
dev.print(pdf,"../report/rutherford.pdf"width=6,height=6)
dev.print(pdf,"../report/rutherford.pdf",width=6,height=6)
clear
clear()
dev.print(pdf,"../report/rutherford.pdf",width=4.5,height=6)
dev.print(pdf,"../report/rutherford.pdf",width=4.5,height=6)
quit()

Binary file not shown.