diff --git a/accelerator/data/.RData b/accelerator/data/.RData
index 8155d66..bd02125 100644
Binary files a/accelerator/data/.RData and b/accelerator/data/.RData differ
diff --git a/accelerator/data/.Rhistory b/accelerator/data/.Rhistory
index 1616b2f..9346fde 100644
--- a/accelerator/data/.Rhistory
+++ b/accelerator/data/.Rhistory
@@ -1,72 +1,512 @@
-data1=read.table("0302171.tab",header=T)
+lines(c(0,1023),c(7600,7600),col="blue")
+abline(h = 1520,col="blue",lwd=1)
+abline(v = 189.12,col="blue",lwd=1)
+axis(2,at=189.12,"189.12")
+axis(3,at=189.12,"0.352±0.009 MeV")
+lines(data3,col='red')
+lines(c(0,1023),c(18500,18500),col="purple")
+abline(h = 3700,col="purple",lwd=1)
+abline(v = 203.63,col="purple",lwd=1)
+axis(3,at=203.63,"0.96±0.04 MeV")
+minor.tick(nx=5,ny=2)
+dev.print(pdf,"../report/calibration.pdf",width=4.5,height=3)
+plot(data2,type='l',xlim=c(185,205))
+lines(c(0,1023),c(7600,7600),col="blue")
+abline(h = 1520,col="blue",lwd=1)
+abline(v = 189.12,col="blue",lwd=1)
+axis(1,at=189.12)
+axis(3,at=189.12,"0.352±0.009 MeV")
+lines(data3,col='red')
+lines(c(0,1023),c(18500,18500),col="purple")
+abline(h = 3700,col="purple",lwd=1)
+abline(v = 203.63,col="purple",lwd=1)
+axis(3,at=203.63,"0.96±0.04 MeV")
+minor.tick(nx=5,ny=2)
+dev.print(pdf,"../report/calibration.pdf",width=4.5,height=3)
+?axis
+plot(data2,type='l',xlim=c(185,205))
+lines(c(0,1023),c(7600,7600),col="blue")
+abline(h = 1520,col="blue",lwd=1)
+abline(v = 189.12,col="blue",lwd=1)
+axis(1,at=189.12,hadj=1)
+axis(3,at=189.12,"0.352±0.009 MeV")
+lines(data3,col='red')
+lines(c(0,1023),c(18500,18500),col="purple")
+abline(h = 3700,col="purple",lwd=1)
+abline(v = 203.63,col="purple",lwd=1)
+axis(3,at=203.63,"0.96±0.04 MeV")
+minor.tick(nx=5,ny=2)
+dev.print(pdf,"../report/calibration.pdf",width=4.5,height=3)
+e")
+> abline(h = 3700,col="purple",lwd=1)
+> abline(v = 203.63,col="purple",lwd=1)
+e")
+> abline(h = 3700,col="purple",lwd=1)
+> abline(v = 203.63,col="purple",lwd=1)
+e")
+> abline(h = 3700,col="purple",lwd=1)
+> abline(v = 203.63,col="purple",lwd=1)
+plot(data2,type='l',xlim=c(185,205))
+lines(c(0,1023),c(7600,7600),col="blue")
+abline(h = 1520,col="blue",lwd=1)
+abline(v = 189.12,col="blue",lwd=1)
+axis(1,at=189.12)
+axis(3,at=189.12,"0.352±0.009 MeV")
+lines(data3,col='red')
+lines(c(0,1023),c(18500,18500),col="purple")
+abline(h = 3700,col="purple",lwd=1)
+abline(v = 203.63,col="purple",lwd=1)
+axis(3,at=203.63,"0.96±0.04 MeV")
+minor.tick(nx=5,ny=2)
+dev.print(pdf,"../report/calibration.pdf",width=4.5,height=3)
+plot(data2,type='l',xlim=c(185,205))
+lines(c(0,1023),c(7600,7600),col="blue")
+abline(h = 1520,col="blue",lwd=1)
+abline(v = 189.12,col="blue",lwd=1)
+axis(1,at=189.12,vadj=1)
+axis(3,at=189.12,"0.352±0.009 MeV")
+lines(data3,col='red')
+lines(c(0,1023),c(18500,18500),col="purple")
+abline(h = 3700,col="purple",lwd=1)
+abline(v = 203.63,col="purple",lwd=1)
+axis(3,at=203.63,"0.96±0.04 MeV")
+minor.tick(nx=5,ny=2)
+dev.print(pdf,"../report/calibration.pdf",width=4.5,height=3)
+?axis
+plot(data2,type='l',xlim=c(185,205))
+lines(c(0,1023),c(7600,7600),col="blue")
+abline(h = 1520,col="blue",lwd=1)
+abline(v = 189.12,col="blue",lwd=1)
+axis(1,at=189.12,padj=1)
+axis(3,at=189.12,"0.352±0.009 MeV")
+lines(data3,col='red')
+lines(c(0,1023),c(18500,18500),col="purple")
+abline(h = 3700,col="purple",lwd=1)
+abline(v = 203.63,col="purple",lwd=1)
+axis(3,at=203.63,"0.96±0.04 MeV")
+minor.tick(nx=5,ny=2)
+dev.print(pdf,"../report/calibration.pdf",width=4.5,height=3)
+plot(data2,type='l',xlim=c(185,205))
+lines(c(0,1023),c(7600,7600),col="blue")
+abline(h = 1520,col="blue",lwd=1)
+abline(v = 189.12,col="blue",lwd=1)
+axis(1,at=189.12,padj=1.3)
+axis(3,at=189.12,"0.352±0.009 MeV")
+lines(data3,col='red')
+lines(c(0,1023),c(18500,18500),col="purple")
+abline(h = 3700,col="purple",lwd=1)
+abline(v = 203.63,col="purple",lwd=1)
+axis(3,at=203.63,"0.96±0.04 MeV")
+minor.tick(nx=5,ny=2)
+dev.print(pdf,"../report/calibration.pdf",width=4.5,height=3)
+plot(data2,type='l',xlim=c(185,205))
+lines(c(0,1023),c(7600,7600),col="blue")
+abline(h = 1520,col="blue",lwd=1)
+abline(v = 189.12,col="blue",lwd=1)
+axis(1,at=189.12,padj=1.3)
+axis(3,at=189.12,"0.352±0.009 MeV")
+lines(data3,col='red')
+lines(c(0,1023),c(18500,18500),col="purple")
+abline(h = 3700,col="purple",lwd=1)
+abline(v = 203.63,col="purple",lwd=1)
+axis(1,at=203.63,padj=1.3)
+axis(3,at=203.63,"0.96±0.04 MeV")
+minor.tick(nx=5,ny=2)
+dev.print(pdf,"../report/calibration.pdf",width=4.5,height=3)
+plot(data2,type='l',xlim=c(185,205))
+lines(c(0,1023),c(7600,7600),col="blue")
+abline(h = 1520,col="blue",lwd=1)
+abline(v = 189.12,col="blue",lwd=1)
+axis(1,at=189.12,padj=1.3,"189.12±0.6")
+axis(3,at=189.12,"0.352±0.009 MeV")
+lines(data3,col='red')
+lines(c(0,1023),c(18500,18500),col="purple")
+abline(h = 3700,col="purple",lwd=1)
+abline(v = 203.63,col="purple",lwd=1)
+axis(1,at=203.63,padj=1.3,"203.63±0.6")
+axis(3,at=203.63,"0.96±0.04 MeV")
+minor.tick(nx=5,ny=2)
+dev.print(pdf,"../report/calibration.pdf",width=4.5,height=3)
+plot(data2,type='l',xlim=c(185,205))
+lines(c(0,1023),c(7600,7600),col="blue")
+abline(h = 1520,col="blue",lwd=1)
+abline(v = 189.12,col="blue",lwd=1)
+axis(1,at=189.12,padj=1.5,"189.12±0.6")
+axis(3,at=189.12,"0.352±0.009 MeV")
+lines(data3,col='red')
+lines(c(0,1023),c(18500,18500),col="purple")
+abline(h = 3700,col="purple",lwd=1)
+abline(v = 203.63,col="purple",lwd=1)
+axis(1,at=203.63,padj=1.5,"203.63±0.6")
+axis(3,at=203.63,"0.96±0.04 MeV")
+minor.tick(nx=5,ny=2)
+dev.print(pdf,"../report/calibration.pdf",width=4.5,height=3)
 plot(data1)
-data2=read.table("0302172.tab",header=T)
-lines(data2)
-lines(data1)
-plot()
-lines(data1)
-sdfklj
+plot(data2,type='l',xlim=c(185,205))
+lines(c(0,1023),c(7600,7600),col="blue")
+abline(h = 1520,col="blue",lwd=1)
+abline(v = 189.12,col="blue",lwd=1)
+axis(1,at=189.12,padj=1.5,"189.12±0.6")
+axis(3,at=189.12,"1.706±0.05 MeV")
+lines(data3,col='red')
+lines(c(0,1023),c(18500,18500),col="purple")
+abline(h = 3700,col="purple",lwd=1)
+abline(v = 203.63,col="purple",lwd=1)
+axis(1,at=203.63,padj=1.5,"203.63±0.6")
+axis(3,at=203.63,"1.838±0.05 MeV")
+minor.tick(nx=5,ny=2)
+dev.print(pdf,"../report/calibration.pdf",width=4.5,height=3)
+203.63-189.12
+1.838-1.706
+.132/14.51
+1.706-.009*189.12
+energy_coords=data1$channel*.009+.00392
+energy_coords
+data2$energy=energy_coords
+plot(energy~count,data=data2)
+plot(count~energy,data=data2)
+data3$energy=energy_coords
+data1$energy=energy_coords
+plot(count~energy,data=data1,type='l',xlim=c(7.0,8.0))
+plot(count~energy,data=data1,type='l')
+plot(count~energy,data=data1,type='l',xlim=c(7.0,8.0),ylim=c(0,100))
+plot(count~energy,data=data1,type='l',xlim=c(7.0,8.0),ylim=c(0,100),xlab="energy [MeV]")
+plot(pnorm(50,500)
+plot(pnorm(50,500))
+plot(dnorm(50,50,500))
+plot(count~energy,data=data1,type='l',xlim=c(7.0,8.0),ylim=c(0,100),xlab="energy [MeV]")
+abline(v = 7.949,col="red")
+abline(v = 7.949,col="red")
+axis(3,at=7.949,"7.949 MeV")
+plot(count~energy,data=data1,type='l',xlim=c(7.0,8.0),ylim=c(0,100))
+abline(v = 7.949,col="red")
+axis(3,at=7.949,"7.949 MeV")
+abline(v = 7.688,col="red")
+axis(3,at=7.688,"7.688 MeV")
+plot(count~energy,data=data1,type='l',xlim=c(7.0,8.0),ylim=c(0,100))
+abline(v = 7.949,col="red")
+axis(3,at=7.949,"7.949 MeV")
+abline(v = 7.688,col="red")
+axis(3,at=7.688,"7.688 MeV")
+minor.tick(nx=5,ny=2)
+dev.print(pdf,"../report/rutherford.pdf",width=4.5,height=3)
+plot(data2,type='l',xlim=c(50,205))
+lines(c(0,1023),c(7600,7600),col="blue")
+lines(data3,col='red')
+lines(c(0,1023),c(18500,18500),col="purple")
+minor.tick(nx=5,ny=2)
+dev.print(pdf,"../report/rutherford.pdf",width=4.5,height=3)
+# Rutherford Overlay
+plot(data2,type='l',xlim=c(50,205))
+lines(c(0,1023),c(7600,7600),col="blue")
+lines(data3,col='red')
+lines(c(0,1023),c(18500,18500),col="purple")
+minor.tick(nx=5,ny=2)
+dev.print(pdf,"../report/rutherford.pdf",width=4.5,height=3)
+plot(data2,type='l',xlim=c(50,205))
+lines(c(0,1023),c(7600,7600),col="blue")
+lines(data3,col='red')
+lines(c(0,1023),c(18500,18500),col="purple")
+minor.tick(nx=5,ny=2)
+dev.print(pdf,"../report/rutherford.pdf",width=4.5,height=3)
+plot(data2,type='l',xlim=c(50,205))
+plot(data2,type='line',xlim=c(50,205))
+plot(data2,type=l,xlim=c(50,205))
+plot(data2,type="l",xlim=c(50,205))
 ?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')
+plot(data2,type="l",xlim=c(50,205))
+plot(data2,type="l",xlim=c(50,205))
 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)
+q
+plot(count~channel,data=data2,type='l',xlim=c(50,205))
+lines(c(0,1023),c(7600,7600),col="blue")
+lines(count~channel,data=data3,col='red')
+lines(c(0,1023),c(18500,18500),col="purple")
+minor.tick(nx=5,ny=2)
+dev.print(pdf,"../report/rutherford.pdf",width=4.5,height=3)
+plot(count~energy,data=data1,type='l',xlim=c(7.0,8.0),ylim=c(0,100))
+abline(v = 7.949,col="red")
+axis(3,at=7.949,"7.949 MeV")
+abline(v = 7.688,col="red")
+axis(3,at=7.688,"7.688 MeV")
+minor.tick(nx=5,ny=2)
+dev.print(pdf,"../report/nuclear.pdf",width=4.5,height=3)
+# Rutherford Overlay
+plot(count~channel,data=data2,type='l',xlim=c(50,205))
+lines(c(0,1023),c(7600,7600),col="blue")
+lines(count~channel,data=data3,col='red')
+lines(c(0,1023),c(18500,18500),col="purple")
+minor.tick(nx=5,ny=2)
+dev.print(pdf,"../report/rutherford.pdf",width=6,height=3)
+# calibration close-up
+plot(data2,type='l',xlim=c(185,205))
+lines(c(0,1023),c(7600,7600),col="blue")
+abline(h = 1520,col="blue",lwd=1)
+abline(v = 189.12,col="blue",lwd=1)
+axis(1,at=189.12,padj=1.5,"189.12±0.6")
+axis(3,at=189.12,"1.706±0.05 MeV")
+lines(data3,col='red')
+lines(c(0,1023),c(18500,18500),col="pudatarple")
+abline(h = 3700,col="purple",lwd=1)
+abline(v = 203.63,col="purple",lwd=1)
+axis(1,at=203.63,padj=1.5,"203.63±0.6")
+axis(3,at=203.63,"1.838±0.05 MeV")
+minor.tick(nx=5,ny=2)
+dev.print(pdf,"../report/calibration.pdf",width=6,height=3)
+# nuclear reactions
+plot(count~energy,data=data1,type='l',xlim=c(7.0,8.0),ylim=c(0,100))
+abline(v = 7.949,col="red")
+axis(3,at=7.949,"7.949 MeV")
+abline(v = 7.688,col="red")
+axis(3,at=7.688,"7.688 MeV")
+minor.tick(nx=5,ny=2)
+dev.print(pdf,"../report/nuclear.pdf",width=6,height=3)
+plot(data2,type='l',xlim=c(185,205))
+lines(c(0,1023),c(7600,7600),col="blue")
+abline(h = 1520,col="blue",lwd=1)
+abline(v = 189.12,col="blue",lwd=1)
+axis(1,at=189.12,padj=1.5,"189.12±0.6")
+axis(3,at=189.12,"1.706±0.05 MeV")
+lines(data3,col='red')
+lines(c(0,1023),c(18500,18500),col="pudatarple")
+abline(h = 3700,col="purple",lwd=1)
+abline(v = 203.63,col="purple",lwd=1)
+axis(1,at=203.63,padj=1.5,"203.63±0.6")
+axis(3,at=203.63,"1.838±0.05 MeV")
+minor.tick(nx=5,ny=2)
+dev.print(pdf,"../report/calibration.pdf",width=6,height=3)
+plot(count~channel,data=data2,type='l',xlim=c(185,205))
+lines(c(0,1023),c(7600,7600),col="blue")
+abline(h = 1520,col="blue",lwd=1)
+abline(v = 189.12,col="blue",lwd=1)
+axis(1,at=189.12,padj=1.5,"189.12±0.6")
+axis(3,at=189.12,"1.706±0.05 MeV")
+lines(count~channel,data=data3,col='red')
+lines(c(0,1023),c(18500,18500),col="purple")
+abline(h = 3700,col="purple",lwd=1)
+abline(v = 203.63,col="purple",lwd=1)
+axis(1,at=203.63,padj=1.5,"203.63±0.6")
+axis(3,at=203.63,"1.838±0.05 MeV")
+minor.tick(nx=5,ny=2)
+dev.print(pdf,"../report/calibration.pdf",width=6,height=3)
+data1[7.8-7.9]
+data1[energy=7.8-7.9]
+data1[energy>7.8]
+data1[$energy>7.8]
+data1[data1$energy>7.8]
+plot(count~channel,data1)
+plot(count~channel,data1,xlim=c(840,870)))
+plot(count~channel,data1,xlim=c(840,870))
+plot(count~channel,data1,xlim=c(840,870),ylim=c(0,100))
+plot(count~channel,data1,xlim=c(800,870),ylim=c(0,100))
+plot(count~channel,data1,xlim=c(750,870),ylim=c(0,100))
+plot(count~channel,data1,xlim=c(700,870),ylim=c(0,100))
+plot(count~channel,data1,xlim=c(700,970),ylim=c(0,100))
+plot(count~channel,data1,xlim=c(7800,850),ylim=c(0,100))
+plot(count~channel,data1,xlim=c(780,850),ylim=c(0,100))
+plot(count~channel,data1,xlim=c(780,850),ylim=c(0,100),type=l)
+plot(count~channel,data1,xlim=c(780,850),ylim=c(0,100),type='l')
+plot(count~channel,data1,xlim=c(780,850),ylim=c(0,100),pch='|')
+plot(count~energy,data=data1,type='l',xlim=c(7.0,8.0),ylim=c(0,100))
+abline(v = 7.949,col="red")
+axis(3,at=7.949,"7.949 MeV")
+abline(v = 7.688,col="red")
+axis(3,at=7.688,"7.688 MeV")
+minor.tick(nx=6,ny=2)
+dev.print(pdf,"../report/nuclear.pdf",width=6,height=3)
+plot(count~energy,data=data1,type='l',xlim=c(7.0,8.0),ylim=c(0,100))
+abline(v = 7.949,col="red")
+axis(3,at=7.949,"7.949 MeV")
+abline(v = 7.688,col="red")
+axis(3,at=7.688,"7.688 MeV")
+minor.tick(nx=8,ny=2)
+dev.print(pdf,"../report/nuclear.pdf",width=6,height=3)
+plot(count~channel,data1,xlim=c(780,850),ylim=c(0,100),pch='|')
+data1
+peak1_interval = c(793,800)
+peak2_interval = c(838,847)
+interval
+seq
+seq(peak2_interval)
+seq(838,847(
+seq(838,847)
+peak2 = seq(838,847)
+peak1 = seq(793,800)
+peak1
+data1[peak1]
+data1[index==peak1]
+?subset
+data1[,peak1]
+data1[1,peak1]
+data1[1,]
+data1[,]
+data1[,1]
+data1[,peak1]
+data1[,index]
+data1[,$index]
+data1[,2]
+data1[,1]
+data1[,1==peak1]
+sum
+?sum
+data1[$count]
+data1[count]
+data1["count"]
+data1["count"][peak1]
+data1["count",peak1
+data1["count",peak1]
+data1["count"]
+peak1
+data1
+names(data1)
+data1["v1"]
+data1["channel"]
+data1["channel","count"]
+data1["channel":"count",]
+data1["channel":"count"]
+data1["channel","count"]
+data1["channel"]
+data1[c("channel","count")]
+data1[5:10]
+5:10
+data1[c(5:10)]
+data1[,c(5:10)]
+data1[c(5:10),]
+data1[peak1,]
+data1[peak2,]
+data1[0,]
+data1[1,]
+data1[801,]
+data1[802,]
+data1[801,]
+peak1
+peak1=793:801
+data1[801,]
+data1[peak1,]
+peak1=793:802
+data1[peak1,]
+data1[peak2,]
+peak2=837:849
+data1[peak2,]
+peak2=837:850
+data1[peak2,]
+peak2=836:850
+data1[peak2,]
+peak2=835:850
+data1[peak2,]
+peak2=834:850
+data1[peak2,]
+peak2=833:850
+data1[peak2,]
+?sum
+data1[peak2,"energy"]
+data1[peak2,"count"]
+sum(data1[peak2,"count"])
+sum(data1[peak1,"count"])
+plot(count~energy,data=data1,type='l',xlim=c(7.0,8.0),ylim=c(0,100))
+abline(v = 7.949,col="red")
+abline(v = 7.7,col="blue")
+abline(v = 6.9,col="blue")
+axis(3,at=7.949,"7.949 MeV")
+abline(v = 7.688,col="red")
+axis(3,at=7.688,"7.688 MeV")
+minor.tick(nx=8,ny=2)
+dev.print(pdf,"../report/nuclear.pdf",width=6,height=3)
+7.949-6.9
+1.049/6.9
+7.688-7.7
+0.012/7.7
+sum(data1[peak1,"count"])
+sum(data1[peak2,"count"])
+3.12/3
+plot(count~energy,data=data1,type='l',xlim=c(7.0,8.0),ylim=c(0,100))
+abline(v = 7.949,col="blue")
+abline(v = 7.175,col="purple")
+abline(v = 7.588,col="purple")
+axis(3,at=7.949,"7.949 MeV")
+abline(v = 7.688,col="blue")
+axis(3,at=7.688,"7.688 MeV")
+minor.tick(nx=8,ny=2)
+dev.print(pdf,"../report/nuclear.pdf",width=6,height=3)
+plot(count~energy,data=data1,type='l',xlim=c(7.0,8.0),ylim=c(0,100))
+abline(v = 7.949,col="blue")
+abline(v = 7.175,col="purple")
+abline(v = 7.58,col="purple")
+axis(3,at=7.949,"7.949 MeV")
+abline(v = 7.688,col="blue")
+axis(3,at=7.688,"7.688 MeV")
+minor.tick(nx=8,ny=2)
+dev.print(pdf,"../report/nuclear.pdf",width=6,height=3)
+plot(count~energy,data=data1,type='l',xlim=c(7.0,8.0),ylim=c(0,100))
+abline(v = 7.949,col="blue")
+abline(v = 7.175,col="purple")
+abline(v = 7.583,col="purple")
+axis(3,at=7.949,"7.949 MeV")
+abline(v = 7.688,col="blue")
+axis(3,at=7.688,"7.688 MeV")
+minor.tick(nx=8,ny=2)
+dev.print(pdf,"../report/nuclear.pdf",width=6,height=3)
+plot(count~energy,data=data1,type='l',xlim=c(7.0,8.0),ylim=c(0,100))
+abline(v = 7.949,col="blue")
+abline(v = 7.175,col="purple")
+axis(3,at=7.175,"7.175 MeV")
+abline(v = 7.583,col="purple")
+axis(3,at=7.583,"7.583 MeV")
+axis(3,at=7.949,"7.949 MeV")
+abline(v = 7.688,col="blue")
+axis(3,at=7.688,"7.688 MeV")
+minor.tick(nx=8,ny=2)
+dev.print(pdf,"../report/nuclear.pdf",width=6,height=3)
+plot(count~energy,data=data1,type='l',xlim=c(7.0,8.0),ylim=c(0,100))
+abline(v = 7.949,col="blue")
+abline(v = 7.175,col="purple")
+axis(3,at=7.175,"7.175 MeV",padj=1)
+abline(v = 7.583,col="purple")
+axis(3,at=7.583,"7.583 MeV",padj=1)
+axis(3,at=7.949,"7.949 MeV")
+abline(v = 7.688,col="blue")
+axis(3,at=7.688,"7.688 MeV")
+minor.tick(nx=8,ny=2)
+dev.print(pdf,"../report/nuclear.pdf",width=6,height=3)
+plot(count~energy,data=data1,type='l',xlim=c(7.0,8.0),ylim=c(0,100))
+abline(v = 7.949,col="blue")
+abline(v = 7.175,col="purple")
+axis(3,at=7.175,"7.175 MeV",padj=1.2)
+abline(v = 7.583,col="purple")
+axis(3,at=7.583,"7.583 MeV",padj=1.2)
+axis(3,at=7.949,"7.949 MeV")
+abline(v = 7.688,col="blue")
+axis(3,at=7.688,"7.688 MeV")
+minor.tick(nx=8,ny=2)
+dev.print(pdf,"../report/nuclear.pdf",width=6,height=3)
+plot(count~energy,data=data1,type='l',xlim=c(7.0,8.0),ylim=c(0,100))
+abline(v = 7.949,col="blue")
+abline(v = 7.175,col="purple")
+axis(3,at=7.175,"7.175 MeV",padj=1.2)
+abline(v = 7.583,col="purple")
+axis(3,at=7.583,"7.583 MeV",padj=1.2)
+axis(3,at=7.949,"7.949 MeV",padj=-.3)
+abline(v = 7.688,col="blue")
+axis(3,at=7.688,"7.688 MeV",padj=-.3)
+minor.tick(nx=8,ny=2)
+dev.print(pdf,"../report/nuclear.pdf",width=6,height=3)
+plot(count~energy,data=data1,type='l',xlim=c(7.0,8.0),ylim=c(0,100),xlab="energy (MeV)")
+abline(v = 7.949,col="blue")
+abline(v = 7.175,col="purple")
+axis(3,at=7.175,"7.175")
+abline(v = 7.583,col="purple")
+axis(3,at=7.583,"7.583")
+axis(3,at=7.949,"7.949")
+abline(v = 7.688,col="blue")
+axis(3,at=7.688,"7.688")
+minor.tick(nx=8,ny=2)
+dev.print(pdf,"../report/nuclear.pdf",width=6,height=3)
+plot(count~energy,data=data1,type='l',xlim=c(0.5,2.5),ylim=c(0,100),xlab="energy (MeV)")
+minor.tick(nx=8,ny=2)
+plot(count~energy,data=data1,type='l',xlim=c(0.5,2.5),ylim=c(0,20000),xlab="energy (MeV)")
+minor.tick(nx=8,ny=2)
+dev.print(pdf,"../report/thickness.pdf",width=6,height=3)
+plot(count~energy,data=data1,type='l',xlim=c(0.5,2.5),ylim=c(0,25000),xlab="energy (MeV)")
+minor.tick(nx=8,ny=2)
+dev.print(pdf,"../report/thickness.pdf",width=6,height=3)
 quit()
diff --git a/accelerator/report/thickness.pdf b/accelerator/report/thickness.pdf
index 5ba4cbe..4ee8850 100644
Binary files a/accelerator/report/thickness.pdf and b/accelerator/report/thickness.pdf differ