From 17b52791c564d164b631bb8cfab7d8189fb21ade Mon Sep 17 00:00:00 2001
From: caes <caes>
Date: Mon, 27 Mar 2017 19:52:56 -0400
Subject: [PATCH] finished lab 4

---
 accelerator/data/.RData          | Bin 3884 -> 14417 bytes
 accelerator/data/.Rhistory       | 576 +++++++++++++++++++++++++++----
 accelerator/report/thickness.pdf | Bin 6124 -> 6130 bytes
 3 files changed, 508 insertions(+), 68 deletions(-)

diff --git a/accelerator/data/.RData b/accelerator/data/.RData
index 8155d66f0deb865f88c7ddbade8c7fbe2100c2af..bd0212579c857908626cb42ef45387267cb4c3d3 100644
GIT binary patch
literal 14417
zcmeHtcUV)|+V6-36=|X%H3BM4U=#%*6qTkTpwhd5ihu+{fY2F5dWlGr5)cHD8iGIw
zp{Vp4krE(*&_aMffKWmU+{~PrIcL7}o%21neE0rydH7}R_xINIX0N^X+D`y6$A9|m
z@9@&tQD0-L6#13WfVMH4qT9bre3hMiA$C&;di>BU9#8NIu2&dW4=5DM&P5Qd!5ld;
z!+`eox_%7U(WkF>@A~^21VJArzcn%`=yfkuDJ~i+6&$LG)plfq8|Zq+UKQT&pO2N?
z7gZ65E62ON<>66Djaq&;pnwcmel0g81a?!u^gx0m@V*_~f2+$XdB4jlPNi_%rIJd-
zB%5^9p1A9V@s=GCzk-zAuWr?^ot`S3Qt^kmsMSH;v-Z9?Y}`T|;$I_Q!$AZ4KKl84
zep?@GP!fEyZhOOZYIKXzwu8cCJ(^vtjIj$i13M8g?=~tevC)&)ts>Fm=%1~OoY6|C
z3hSlkD`%$wB=IX;s}xgg0s5>k|9UAQVT2D_r?L8Mth@=WD^oAHJ=IiW$sT{od8OK9
zrA{8#jh!3YJrz)+;ndDGtF&hGEq|ORH~adwBH;)4MgCwY6t`Cqyv)({3e8FpasIim
ztw7F7{sMrCDUU!sbze(E*h(DuKf8m!+ofoj-B{F*YT(Jhbt=wm2Ag*|7$M`~BAd1l
z7X{Zb;#}{J{&hDfaZ_X^x4W`K(8p|~vX<;I-7&iw@itmj8Ax{ZjC%Ia$gqDlKZbAn
zO_ieK-1|xSqIU>EQ<MvJ(~K89rqLk1yXB9%cW<*|u{1P;sA2d1ekYCZ;;dkHVES-i
zY|m+;I<_t)y^#PZn4c486C!|OYaceL414#50U@13t^sbVj^vpBq5(zd(s$R0D#oPJ
z{?TCDs&l~(4v2{wX>tfVXTwg{SM|XI{RD9(;Pxe<I|{v-<Oyq9kMRdVxwApGt@-&%
z^}Fb>RJVC_15vMFvvXhldy$1ZhIQk#(=z2C|7lAGJDjQ>l0DB`?$q>CicM!d@lsoI
z{)d!nCH*Rl+y0yQoJ-m6DvPdxK84z$olBY<hMT3s&PHWZ#WIln;lN4o&3WEXH%6v@
zr>d;7^K4SVQ+g87-=0p*8{moYsSKPn3K`5hZJ04uIzQ%oFY0D-Gf-ywzHHQLogZ;C
zWo+zev<=#I8k3)`F-friF$U%E-40H?<#j;uSw9$?2#{4}$rZJ=!{?dX^sNx49uThT
z_UIt7V?T5xB2lJ3c$wxn{QZ+mafW`3r~35<TDVHft{bOa<N4y-I)yC7U8@W0Xs5?r
ztD8tV3vAKf*BPI7l=l)tw@yCAVI;rjP~&I*ehMPHVyMn|sdCJ<!g6s$^>g2sP~O4T
zJkAAtWn0~|&r=_uV_x^K6@z{R2Us`G(Rrq2*<S0b;)B}1Z|VY5fLA@AY|il7N#8GV
z5ApOrV5BC(>hWHEsEs%CT~0%FyuQSSmV6x(+XnZo`I3Ft$xjThE-bTXy9NA(FN2=I
zIy+{88FwUw&%GMsIb2Y>=}*9Z3NX*#M+k<a<m=oQzYov0wU!Yc<`1OHk{|}&SOKV7
z&*HbpLq~xZ&Y7reAP0BgrS818AtS{6sj8!HDew*2xgVQe87p_{MPFDElXcw~d_-_s
z9clh6D58La;q$aFm!seDQYlK%3S&^u9-I#lDJHFqG)tEUJznD*0`c)xmgcZE0v;YP
z3f-`+i#uCEtKs!&OFnAyA|9oa1v?gIAM=C<s{Fg+N~uqkpXgu4lvZ`c<jQniGKh*-
ztC4Fd82@-R^%7M-%QWw%`B#2o8Eh`Guh9#){W2VwK{q9?RqWTZT8P@Ud0B7euf*0D
z%qMpjlI>8om4078*LD$~S=a5Tk*jkk5oSm|BTB@)h6P9Nhs`DS($lIsg=$g__ft!D
z><YZ4-#ozi2IrKRCC)m<-k)ZJ1^K2|2Rz0&L3)C0O#}memTaY3(@)r?OF*O57ng-V
zW`+7G8EvbUW(qKkcY;rrLD%9c=iewkFIoDEk-~<8Oyb~`)9^Inr=Jx*yQPM<rd2M?
zUmw?SyO+GQ{3s}38_1~j(|GSRuWVApUC=!(XjI#+2dMO(DRSaSWUK2L-P*d}8d2O<
z@eI4MqG#=S@v)z!SurYK_efkiWz&9ZQ)1b4MCiEHEtrx)M|{aunJoeky9uNVu#S6G
ztgovG`8m54=iHq(x+;iyYG8<6DYnl0-r06GI*u(KU4e3wIJgQ`AE9k;q=0O!E;S65
zg^fWmbw-OUv9JdJEmt>&Cqn6lR~BF63b`|aA1H9$l|&DkpDYa8WW8`%pvLFzQ}TQ7
z7$+Y9F~noSiFk5f0p(L@oD=+|#}@TU!rL0EK>mV_=(m)O<Fru84VU}q+3~)4m!pa=
z`vv+&3$SRYq%Ee@GN>_l_si->ZwfCpL*-y0Vdpb%tE!+!?<~eay}H^tDt5{aV?rIl
zyz5E|b6>6=^+$iJLD)R2t|d>>AJ5aflb3ivPDfe{dX-J{R%Uk`#mC`uX8Uum)Se11
zS5_<$#U;Rw206e0_kpEHA&<JZo8%X27GxA{E1-Ax_C}4$#(h`%=D-oBiv_eYmb*JF
z9JVJk<0MjQyEKXG-S4s-XDxdI=~bF=1-xn2e$i(~YNdp(qkQi&|A@=bE15ns_xnTH
z=vvj%@!=IGxQkMN_cwz0+zgsD#`!u=A{)NbqoF|3&X?sVp1pR}F4aP6X>r_g2=;Nw
z^-Y@MzRz|jpfUDT^pERV1E~0XJO8z&XJ&%zNXx60?bhT=A{?eEY>)cwUVyI#Q$@G{
zxT4O^79#}({Lq{yWpvh4?v}#GIYujLVXc$PX0o@fgrNHbys-E>+4bHs02R_)UC)7K
z)M}JGSRDIyQN(_kigHnXVTzU36zJC9C&*ov^j;M%P>%r!sq-(r^)4D(N=i)Dyvq0D
zzIkTDcd3cJD~(rEr+(=0f+1*+Z6}weGZCRARayM&f?80z?_gE-2!Hm)u+BV2&9#Q<
zI@)vv+2D%0KYDO|Lc;*=hob2kw}ta3MmM}rnJ|8oX@l>jCeDInEY!XzJ|_{TC{zx$
z@9mVjrrU@uK%*tD&hOGEQ18fz_y(dYKXXZVmwjb)$;lR6G7GHoHZ$068qAGF)bCgh
zls6C;CFV^90cocSm-d>k7%2zrfANN7EZ8~l^qarLv40Z}im0^RWq8fzR7VkC5VK^5
zs^L8ZVbW}t936~*X_Rk}W~LZ4*<bS|I+nrPyB3bkG~RMpe?NLySEs>m=ktr5UY-*g
znmqE2pSJA3e#+-(+`|<KBlCJu6Zx3bISPl6kxxs%qRe=Jm;KoV3B2tl^}7smvoCmz
zYHKfp$d+-apbUSiFaZ+aWe3%arlMOCzFo@(#in6rRx7jKt-G{O#1KD7b4rrzH_I!`
zim^zxj_+v2g6`4?a@M?25j1e&6!9|@!W(1`$?c?gD@i=!Z2Ops2${cu@eiF?02CXD
zkk&^pJ_9>p+J(QK1+Y7-bvbKA8CAiu9@GZgf-IzLk#=<u+$2EZMS865kK*<(&5($U
z0yR+G(imT=_~#D`7vh&()ul)cZ4rahDrXv@xxoeTyGuIkmT#WoTAqLV9E1;0<mab)
zDf(@GY*0jR))xst)2D-EFr)KGC8%0}an(#1_$W{qo%!sBfQ5q=ePs;j!fVy=t-&ro
zL&4>m`ij)b&E^#m;n-gKy#qhmM6>gfK~l|)<jfZ|KYXVMA0KGxJ@L_G`F6L|yor%5
z0-z)t>35+=L>-Ks{ixKSzMVPHTmm-8X)49+(6>no`~EkFw3gF%gV=&atcrs~Ddkm|
z;#jxI?<hi<bU3gdE~g)vw&QF$G;|IBdCl6;Kot_w3j>B(N{nsdLcGyG6mzK)?RGgJ
zT_Cpb=f(xvY{;f9R^^HmGGI4&)veZ>GF~0AbFA-9!3Dfmr!slaoOwGhA#LK8lKIEF
zv^d?&bS=tV-Z&zUxhLg!xY?=M1zKA?eVV|p>Qq|+YH5{CEU0{_$UmgztudSPVqQ^l
znDWM8q#9w+F&&UC1(aoMbO05S$KpN<t5$=|Hsv&W;z5~+#e-{f$d6I(&ac{MHR6ly
z9x8kN_<k`5^pLb-71MUsWmi%cu2EGP>-`}9=omWvtXqkYJ46g*)vQEneitXvHV{~v
zNw`+vALyInxv=wnjC5OGrT~ITxsC{a)FC<G(PQ%z=BXII>waSVeI~`?$@r%ty`;7I
zcsW^VtnpPVMv*dp?cLk$!U*n+&vxjaY4cw%(BYlh*jXu@lbfkrS>prc0q;^zsSA5+
zwftSldS*WEJ`mea)>Lz8DYsD}l4tLaT!e6ewAm%Jn%lLfb)HpUwk}32=}MbHzniv5
zus<x*-ds(02Hr&bkV))Ub5cDSFw=gFxrCHbY}o)>r7YCb(KJuryKAVk$6?nH!QXF<
z_Itm0hi@&sv$R%5(@1yF6r{G=8RRUGL!c`8xrcJ+bPBsBHn7HpCIi`6ljK!Oay?Ot
zczVBuGl?2iP4VoHmZGPR&kaPKF{ai>$$BXd>UK=g56pk|t$Ta@5@U~anD^cgDg0tj
zNZbSJ?r@w@(G^ngZoj-g$9)s3yMV3h;COZ?1JFssC^M2-?TK;lqGvLkmtMH>TzOwW
zoACdUa|=V1;bcT_fwr^qR2TO=Ce}LTlOyo6W1I81TYm5GHTUY_uw-3sj?X$lZl$W}
zZ&#m&v|u_v;KeB+i~aZgtUv9f{IY~gZX5Glv<h^7uR2~;j~}4%B}-i1UhJ?1cP74b
zXbA@OJjm+yRsGmqaGZ2^=#`Ujr-)e=ubZW#E$jiw*{~p=PYGb0MFnE-YaEfm-q)2l
zk~&_hTFm4Beiv{y$q5qd=h$zQr8cTOv5C1!^6Q`~3EwL<%A(6YEPRH_?Of9|bjmBS
z%c(23UU%*n3?XTU&U$C3nCI1*PgHai7sF^%BS-NnuKmfWPg{M2S3C9n(oP3Ib4Sj1
zc1_ocEd%i~h!;?pbfy0*TG3ssM31exG$JMGWQO&VB-*_f4B@zB+}V_c;F9}BwguI(
zoPp`Zb?Qsi2%E_>q`2ve*mWObl9FsIy~}#s%I|aTtoL*qAX#bIIjrbtVlv8P(?Acj
zYv~6pp+`K`14_GNlgonECV5D364~{`qLS1yi+tMU@uZSE*K(xlajE4Un0vs|m@8@Q
zm4T7QeajDZK+CmA@5j2NF~|XDn8M*Cy5a!WkEJMDLwEhC?2A;XWx*w3f|I5oDX4RC
z!)R(Ri5}xmwAj6*=ND4fFc@`$Er~9#XBPCJB8>J`j<c4QL+Da2vvHROq-Uy{wE9=m
z?7vqa`X~ZYfRt0R1oc*nZI23*CQ*c(48WF^GaCqsrg29@&mcYomb_uC+zv0LAxF+e
zn8C4903)S{cb4%+>T(%dtdCyOR=rZgZatjsY(sIA1siS@26{@bQ(#xblMIvN?#0Ro
zlm)tg8<m-cQf>#sn!=n?*`|5R4Q^RQft_=^dtF5?odtCEu_afAu{tKH)MRl)gak3c
zs9)FFFjH@9w9VpraCcjEaXtBvL`5>T*i@kBd87VzAjo$_|D#$m7GQ-^|C}33Ytc{J
z%WDlF#kA`fDoWa&Zm3@q>NDBCI+gAocDp56-nF`QNPYA6_FHAY$K+;{PQHY=SiH=L
z$U_vboJ`&m8STrPuhkw^;M0v8XwReGsx8<@K$}h{Eqb#$iRMd<@R8w?zB6TS;bRAE
z?h1z63xhR}6^8C_Th<3MZ1o7qh$W$LgOlj8K;bo?wFOp%B!u8qDnQTmePDgvs9YEJ
z;r5PO9*yx`o0MUQ2lbvvOmfc|B(uJ~L9$}a{+tj~lI%VmJad6|iiH0t*mBG%?qP$?
z*tnO0%WK+F+dv<F{bMKe$>17mycDs?F{>dG@G)@NxO-4iNo!lZuOxh#B{z^rt1(M<
zVNF`bamts5K;O_9<X3uQ$Bs&A>>c}YWn3c(a6Kc%w$?qB_Cc@3l+|+GB>28hy6?2p
zQBqqq{Z3DI5<m%Ck7OUs-D|g@*vxDK?!_O>&6OwJyeN@7_KPo8>7K!(xrN|>At_~}
z4W*CQDJ3(DTTiR^bst_NZ<a<E%&Hm%zlAS}m{{cOUA`8ygUE>LG&E+V@9%eFXKD;d
z1c0PQL8>Gd#x;iqL1~Z0>rW+TL$f*fcY3V0FBUbG;5!k&Xw!X6(m(7)^S^-1`}2^@
z=(AtlcWR{@pP5o!sQxC~Rk)|1m%{x?b`}7vv8%CR$MxR5@f&#C4%ym<J23vVdmzGa
zE06U32c7{w=g>5l%C|sZ;BLZr$~|guV;dB1D23+qHn2m*O{r`?LvAlzw2oN#c{)G<
zcQg|{0T@ZA6q(Nda=OX;O0n0V?Hsx+nbIBGiMGZCX%$}junZn6H*Cu$3}d|&@ReFo
zleH#%4R3;_tWl}3Y>y&-9321Pdsn?QZj3&TDf+;Q%jq;3-j3s$UY3WY+mM?}r@yPj
zK|LR>*S3u%D=42RL!#Gk<C&4tX2pjjO4$g2E5)i3Y6a_%z}e=GI}Q<#Iyx3XGJ5fK
z?NOgRaf;vyjU7v^$&Sn?)VsB#5^0GM+WHTqmA6%0#{+P0Gt<*j00ZOUMzFe9z0Rw;
zZcZ6bJc<T-ZWYc{_U05?rNBs)ZjZ0*=kpZQlL6{x;;;JqhV(Ofj!#?*`m*a}cwgCu
z?%=gdGgMX^ZvGMuAg@nF!4u0aPa7LkUNzn5vEA{~M^tB8+~6qlx}2ndxRyHI8=RRO
zh-$DEo4jvKDWxkHLpKfH_TB0oJ-V^?FyxR^5Y9Z#QYG>;lB&mi447aVl|SxP3`C!~
z(G&BzDGI)+w)b^})OD?BQNK4Uyj!P(yomq!BP#V7L|lfL2w+<$>8SHq=5Ic4j$?r|
zSON;$<8A{DwK~W;>>yr_=ZrK#?@v`B_}MPPr-S(L94|H{0lO9NXJdtg=U*!)WqF}w
zos2YG45xQZU!O+xCi{sQohh5PGeFR~>gHt9t~<I97Ib($9}9NOnBt&>mrS?rC>Ix7
zpTowrU^q1SKtkhWoluS38t?~+4o;4<q*l1VyNmb?$aIo8TX401FLS3=7>SZ2`1XBU
z1J#}-_!hfio>u94Oiv#UNwye{7qknx&iI1y0PR&+oT&4jY}GInu2+r%^I-1v;E4h4
z`g<uMd)5m-Gw%ngif$`nphpn%L-@5qLS9G2*jnH;I?9(mpg(yjv{!=&{t}UKzfAE0
zVOB!PqM~=np$kuVf+4qVZ3LiCjH}h}B=_org6nY_)f!5EV-2wRy-S1(^rD_hrDSxW
z=CEg@3y5G`7a)+`qtCo1PspBAihDOe@``dJR9cih{P40L?feM0uIm)w**}1PiV46?
ze>4(b<t;g-WK&(dQ4_>L@N$}57BO8$&(HHGnq9-re1yivBw1lzW;85;S040gfpkQ!
z#m5*-Ju|2xj%Nur^_&zwH?nF12AQJs41L~QREkN*1Pu?T5>73%EV;vHpLH-q<Rz=4
zGxF60%>%zkf)wi$`pXoAYeu3SU0wZ$W}Z#iJSN{JTy|V7eXZ6?Sj-B7)+Nt8TQ-t?
zlrBP={eeb%fwIV|FEjMWe+YYWN~yg)`4GrF!CK|d!`mvTc1T00)g3bL4|JFpD1lrj
zgnH8}n+IW=Coq-&aM%q2^~mN!^yUeVR8ZXz^>C|u<dh%id@qm`+-g=i<l8?e$HT4e
zl1+c0lf6L6WYw5Ig^O@2k}~Ji-!#7o#Fo-Li0eGTQU%ouv3+T!M;`ituJQuO!mY^4
z+22a@|5Spy{wC<CpgJJx@%TGr#UJQ6FOWDKKcgJd>Xru)oG1LI%xeEr=z;A2g1<+u
z{(&y_0?EMfv&x88w|t1zJi$~2)dOjW_*0odt~-rVpcgwWyS6L|rFS>_?6IQOGS?k~
zf9Qe8Ar-fCyJSITECGDGh7e5j+1}vis0{q-1-gx5ag9xLMe|p@69iL>QW|`&oQ}T^
z-HWn014=z-Liw@7K`>anc>P89(#~S%w=!OoUXWeg3N4xdx_=X}nCpQlT-B`(22c7x
zkg3;*uko@#+$oGD0TjA=@-q&G9_xM2NiMLP^@;E|&EO79!bjETwOOIZ@>w0sV(??D
zKi5fWB;(Zz`@N^(%qQukPc}i}P_mC5r+}GEUphW<QDt~=YBS{`Ig8&Mjt5>;&v;aR
z4Moo~D_gZtL$%qyb(aT%0MGr7F>~53E(r&9Eottc=KDw@DA;IaPISLA$nIc-b8MVJ
zXvlA9Wj=_bgLZX{E}8dCZ+7rdhlna!cskJb3?`PKzT|{p)v>8z)ObdimbvvVAyK%}
z;4`S~#=N`RnMH))mY)3Kr!$)wi$vk5P3F_9_$!@|iiu6#a#ZSm9NuZMW>H`#6L057
ziLucEg&-b~^NkEWcb#N?yvk}5@iF!F5i&Ildh9U5(aDo_vC1Je0cdXKcXiLyojUOr
z71}EtxaRX6t$A~(P%!lyT9bdVibE1BT-188Qhov}{7t7s2HYA2Y+l@rwy~prUx5sZ
zz+i-1i&Y}Q?ZPV6Xq)1GZR&S5wJ$@HKSR%XAX$Ba&`lm|hapz-Qye@+R@KON-NMf{
z-He|=q@m_D_N$nW9R&NqM5W;k5{Fe_(J5=;HsNxdfT`IeM_a2+<(6$UFcY{CWRK*8
z)@*(d|1})rD&Df;phMX9F~WP(mAV3cxL&*ft<mJzY{+)<(B;AFGK0leUq9?ex=!eG
zdl22h@S;8+_L*9oLlU?Hi2e8I<apBFG?LOK9GX#j<&-mX*K~iIR&wVEG!`iQ%Ep&E
zd?bL#V5xDW+*1lY1tK6Jo}&JpvA|91)waT!7t|XoRQ1Jg5M3laXhj;Dd8YphrU^{F
z5itBs*9VWDo>D`YE?c(>;{r=wGlSKX&XwPL)`Xtkf4$bIxH%Iq@hm{Zm5>5__N*FE
zP~0Z$o(}pdrwR=;m&MilSWvf$pfmM|V8TN$ZeA6SHanPzSE+VQVw%}qW=@G+e?|kg
z4hvUp6$0ix@ho6r?kLsO+Tqb5uG~PdIJAhXrV~BILd~l@@S==sQH@;>*{OPRL5m>X
z)@_gO5(?8cHXh+DJqxNQk2ai&pER&Bs$G&1>W2Hfk0`7Mf~nWG@|~Xm<f+#zy{Xr<
zy{SoMDp1QJCJI!bqd6*2l}}9Mt3Yd+V{NGl)W3uY@l2&DgOMVL6MTU%a@(7_$W(4V
zf&=-C;lNbpX!a}2;6UA?UO1(&7jE^b7mhzcb|#=R76@p>0>Lkmi4rE#naE-S&BO;L
zQWpr@<xIpfQNcu+M?*YN!vqd&G=c*|jNq6DgXoy#1wxC3H+AZR3RH?&juf-3!dw;T
z2@^O*tGL%j;5ZqF)~j0}<cBU0NSuRcy^uk4*UwD&GMz^wr{!G4<1r=1aLfx+IA-x7
z93#Z^%sku64@EO)-th>IQ8j~O0w2IJdrSy0L%JE0%gTx%X7vk%DrV_b%+jyl@umvM
zWv0%QA&Ve59VU(pqFp7JU}kiUULbhCW+In~6ebdwsA8g+iENMY*BB9IhAbw=nHkiX
zXkeBT%*?#c9Mg(=#}<vrAZigrNn;Rgb%hDP$-r~a;*ToO;zA~vGj?T`>dNc@g*kHy
z^9vHpV}dzmFgudL<gJ*zl}x7l1PE_J2C*{3Niat~=7?ef#sq>X+c8IU5d><19eb(<
zP2C`<A+m@w)QT)jmp@gGg2t4P|1Jc)T_D&D{nc3TAH=`t7*T&1t36Spzu6#&BeGaM
zIW4fii<(1_zY8Do7YKhLn8r%eA{AKy?fz6z3Ob^M{CA=4?BHJsCtEeZSGA{V^rQ`f
z(GbLtnd46Zdu|Y2P(sdPhU)aEN>R{EV`B5)^}iagE%>Sa?u&IS`Vanf{=pf)x?ZH9
zQ%lIU%s8aMwJQX@A&9&Mw*HmLZ&b&hdcO@qQo!HEznQ!JVP0L>ThHvq+l@a6+JYbd
zKZIsB9d3sGT!Shhk278Q%-{h`N|u8Dogfdc@e-_tAnwf6%qga95P*m*E6?9W+ai@;
zslV%ci^(7Ui1oKPhPG;mFG0)_5;h2JdfrrV3K~{IR<Fp~Wi~p|42$vrSYUtG>kdJ(
z>-@9+Fn40ERVMX!y#(`AA^q<9dmO)wK%rkln+-wiWCqcQ5^@7GG_!j$6tugon$kar
zzv`Igtn}pH%!&3z5IbgS=B$1v>ilive|P;o&Mb4+`detGUC5{24X@#PNIu1*zA$*}
z<0QAWI~p}yDL9B8FY60yxF`*-?{8D@2|walzpJ(EV#ZQ4UNFbqDCS>A2afmoz8A#t
zf!x*>U&-QP44qG=dj?G`eiYo9^$gmk)k^V?SlcN>{M;E<Wr(4Rwak|~)&Q3s(YYC$
zG6ecb$~Z+Drz<mPPE#)9%LWOl^=`2H!<eshvMCuVNtU>$f;f~l^y!bU-q!N!Zl0E{
z8%qE|0q=Kh8#a#&jYoiR%E8LQf}>JAC1A<xf?EVQ=)fNcSh_#CbIaW-t8$@ZL$ghw
z)Oj$B(o%5E3h%+FX{y3xT%X?QkZrj3+!jVD2+6Q=c07Z_rin~FJ=3|t#*#g&SoUp0
z?`o@txm7sD%ttc7+AN8JLkboz?1C8|kA}9Du21YNnXvhEA;~*qX4rvZx89AZ{*U>u
zDX`I+TBnfdoy6kJ{aAnOhEtQ3ldNy=1|Ky;c>7(;hMCb2;*%eo;*K5!BxZGN(5nNN
zcQ!BJG@`A{CDza5Gy*I>5W0n}5B!^DaH)bd`=P`cK<7zjYtr>?zHDfu6O`AzF5$?*
zBDEz8M?gQPTR$0S<|aJ0-3uf~gw}lBd4i~ku19~(DRoLi-X%=D4QSuC6;p~d%;>D1
zI2NQ`j29R52t8=&VFk5#XeQ*f>fH;2gl?_7Kj%_PvZ$_{uJR;2K>9zqkDMYr*x1B;
z`MOW|vfJxXH`BsCX<+i14G7k34&PN`HTTqfG0+}F=4Zz4*$ej2{7|Xa<R3B<KC8|c
z%?|EhFZJ2g*;bthZ$8#`&!@42U6=^1)(B1zeA|i-+Pu-juEtVrXQ<N69u%)gg!u9)
z<yGdgSu27y7huWGa8q#8N*Xe1mk-RW4QsEcbpM7@fiLDR;pZwI)R<bd;4^-W{pX3D
z>=3@V`q5f>rOPMf>Rkh_DqY3|yJcwyfXv4rk6pTaLBo&D-uKMv5XztS6sBJS6Usj`
z-%{$kdNz$m^bS$Zg1E9Bew-wD5Izht-`%1o*L@O&i`l?U+bXW)J=$cr9)iMAMES;!
z(+$!$2lcQ+_s9`JO0(g+tW7Kq!JfZ1uFnbNo1nJZL&gzSd|b9;xM4MuMjckn_qMa)
zh=Pc2|By7n;A)x7(^+2`lFEV_yOMN;G@SU%h6vTkVx<1ThAaWQNjtPQ?^=_cxb;)e
z2|9!i@V)?|v0TUb9i46Wyon#y=FxZ?3fVC;m=12)kq2FA(4ez<6}Bux$-yMW`c-*i
z)UHW<+A|fEm%B$b5PMSwJXF416*}|wMvr*}-yJk%!VJE9{ULm}%9shJT*Z{Vdm;Ny
zcBji1)4wY<=p(kM_5MC@h@ux&6{d;1)K6gi>#P&`-YT<^8t-;U4?pa=r@gI{uPGwO
zD8zAJI!?&$JtfL}zTt?gRKx3{amD4{7u-I_1707#eHnjZGS=f(jTnK8dYtpM#5cCd
zxRqO5V(Z5bf8%`k(<DqXQu4KA?D3mtZk^FSqj^R*>1L8v((TCQktx6V&hwpOouZxM
z(&wecq(!C0i_RB`z2f>B>2h2<5XbvanETan?bo8;xUX|v=dR-7j=ZCFMofS^G*Y`+
z)J1&Z{DRnm=z_Swc>ys2Q33Jv^XX#gqUqw6=PkvixJ$WCL=HtN9$#%4*U|!N1-$mP
z2a2gMe~Hwb8io;j%QeQ;5vjp7#yu4vrp8^*RnHy7HMQHu;9_te=8}qh8F?%6You{A
zmzIEoq{BIf8xBVtE;yXI9VqN~%Hgtukb|lNi-VZMX@{#0CvOLG`w2LxI2?1h<iP8o
z=)fT+e1-1{=M@0*CnOj0IPx6w2=WZ_6jBJuf;^2pi4;H{L-HaykRr%KNN(f_BtP;f
z@+^`aDU3Xf<Uz6_1(B>sJ|rg+0R9Qg1wIZw2R;Hm13m>70<(xVBm|uV3xJP_dcn0b
zH8XXeh{>Io6O$8_6E8boCRQd|CJs6e5(9~X#Hr`0VpP!zE;=_S7biE23&w58Wyp=?
z!g9-V$#WyP5Zvxu?%c~<%iI7i0CzT5Hn$zuNF*ZCJ#sk`5GmAb*etJAb}xGZz;bY8
z|HK~8mc)wTMD`brmmz;4hDHAq@*k}1R_ncw{7dh@8D8Kz^lx)WirI<%`z%`$Rowq+
zd4I1Y(e(Jgn8F`6cdxSSM{SAz^lz%YtNn-Zzv{@PU;ELxCCc(2`t)D<{a5q*FJ||z
zzMlWqw0FxaH>LKBGPxQ-IV*;f+21t6LQWo>*=KTF1uKRvoXivEKXCtK^nZqso>Ede
z{%g+vGvWSC(ft|&0G5aUS5p1f#Q*jQ{r6_+F(J|vwwzY5E3Ba*#|}>I^Y8s8PGo=8
z_zNL45XAQH^#A1WPegMt`@e|&Pp1E$_-ph3U-si5q@Q)mOm(<;`Bs<xm~qP*A{r-g
z9IwT-cT7~`h`8AidpPS?o$@U~$_If(R!1OPtD^8J!ILkfRvxh*yWyF~qIYr+TtZn`
zF&FMt24Ru<W)ZlB;rQDV%14yTRcN=pOLiBApWS}retShF{4t*W7U*Ocp6i8UPP?9H
z`<0w_`<!;poOXTDv>Q6=?Ai?{G;L32<cdT%@aHUh2>3~z+P-#%&KK==^^~&esjlM@
z*$?VkCu`hpc)Hx6Jv-<5{G4i-xMBEo73_<T5i8!;rX<yJc*y0k#pDg)=%pJ%%0MAM
z=`(nvb0WGovfOXzK8C+`v<VX@zZO>tv#)<`59H7+I|Y<kW<@Jq))o=B3!hZ~Rt^vU
zYH@=WbS{M_-1}@i*_v(R{jn^^8;jxghfg(Ci-ZTBizW!1<g_<;xA(2H4;vs0g$oIu
z1m?6`<dmO%WzF--`gAldpH~EOgeos|UFJ;J$y>H3(I3hV9?Udw%|vm{G(TANf1JVj
z)e#Y|$|q*4FcZx+ki9U!Krt2WDvkHfd0g*Nf^#tkIS+7g5*pE;Uf!bB!+a0;#E1cQ
z3A+lp5f?VG+6SVIhGIK_+|;WNvv}25`#x~hAL;Aosz2V>$i?97W965HY0K-Z`)d((
z*PXPMb=F_$bdaLP4boL*cdjOqoWX*9e5V;XPy7V?ATAOIys0tjiE~C(KMK_@b6=sE
zTN|!tXkm*JPJu0W{4(}LQs=Uxvfj-_45k#IPxV2H&l|?9h0f5LceokX`*ii*pw4p3
z0}|)ZVO=CLeS)^@liMiaA3~Zz&Yg%BvH#*uj@o~G!tWp#_9>MdqwaWU8<lzrY`HG1
zZTG`u-{xlcoW^J;B})C{oK^Jx$2nJ(PDcr}HF(rkuF1sKj$pd3D@)oa@eWytQR$#q
zb#akL<_AQIx2)E!OJ}a$pYgmu^VGpFm(%qEbMw_j!}1i&yORYOKAA*(&+=L{7wa@p
zu7IYxgC}a~?7!UesV8E0tvEs`8~fvw{VK|V+2!fwda4HP#{PKQE{s6c2(TNS)KGs;
zF$M0qtM64U*$oIh%+O-*Ik@~FkVDcXV}w{9+t@Ow;UM2w&9J36PoKVSRfphO4s9m;
zPkq){P0R7+*sh`<csCy4*S>1(r0;!3G=^-p(NRK`+bpAjhuj)|QM9|6A^*OI@B1v<
z^R=V@`}Ce`FW=7w_vX?igN;E2Ja)%Tu-DxBMMb5gkb3*xqh6Hc0i%8U+<iN?{TEke
z$1odGX^oPrHh}c#5Kl;pVZy-@qiX)`ckzRL8n%)vk?meJB_2<q1y}E-LfD(#Ln_h7
zMy(IKmAzEN-OO-(-*f%@tiW>!`-B}+=>P5TuNEr^FUj|_-CX8nQQqfXTAvLgm{219
zRamygSnyAdRHN_AMBb^1`UXDZX`3g^g4ellUwf)LVfbV28>`#|o(~D++!k-7u$DbT
zhnPaP&xq3Ljl5+Lryl9J1#9(T{?7cyvF+=}P;s~7W4Ke@2Id_S#2r?5TF#@G|IO#!
zo1@nf939>$IC?Z*X}9b16?Q~2TEv!{9q>_M5$5JmVKDQ$cx3w(aJwOLSk3|PsqIQT
zC%9c3xxyFfc9mE8f*e-?QmyTqe@Trxq4}Gu5LJf5Pu@qAMOn!dc2h1|b3aNO6jd5G
zS2(VC3~zXOLGx5#ix@yVrY$t^bVR0bs5f7P5RZ|GEh`>!s=EK(bsGBdSleTGs~EBQ
z>Et;-#bXG~YbDK>&D;Tpaed|7?4ehB2G?@!B?OK+uGDIaBed;-@!xhLTi2p=vS3rb
zL60v9c#EFKN}SxjAaGswB<0HKZJE=Q>!<w`Pb2g%c;2}%WOIb7!Gb=|QYCSkqRN85
zD?1$h%cN!64Q7a&vL#{OW^Cx1V+$tQs;ry?%0~9abxtvf%WEycI?)6^my4Q<z?9TB
z2cpR$Uo^uOnJvAk*AiSd4B*{ucgVI}6pyie=Zlcrah+n+4+BufSf1S_KaHGkj$IUq
zx+HQ8Be>6~eOW(pzeU;a@(zFT%yuqcg+-S9j#Mu86%L5NiTPhLXZ2I#TSd){%+Q^*
zvY$4@?(M_o)@HCST=n986n=5T81)B5bq2wRyXz@hjrZ2ewG2;EBoM^w<{~l7i;sPj
z`nx%BRZ<8l0mZGR*cZT~rrf9Ln!9VXTvD%Uw@%htX2-kOx*J%_`TOODNbY#%9!Ttj
z<^m*le#vcPUXF5?C3d27aT5M9XG1U?spKfz4|DR;G@9pCt&m6SGBF+at#(S(c>aNV
zYE0uA`TUu;qGdAvSt2y^b#3hs<8^&)tns?R5@~y?a;|=~qp^jgNUPn3oS_9&4Aq7J
zfwwYw@o&(oy$8aDQB=dvr)lZUG*$5B@pH3S-Sf0}9bL5DmhCmJ-~+W!R8;>#8Zl?T
zs<Fd;l}^~lrmh0U4(x!2Li>SPO^wj4<#+lw_dDksy_FB<+%*q?!KiIq$~K@r_@F;n
zbBDO?dr(js6OtcVO|1&``Oa}5&H7I(^RLVC&jtHGO&>ghu-|W4W?A&AZ460znYO3>
zGXCR2aJKd1tj6X>NAs6kS73@IHw8`Kv&(!BKl{AtWTpMl(Q+0y^Fw>CO((s6JAC<9
zi~VGd+=!P;P5kJdv&uai=UFaU+qas)e$nSE5jT60Sk^8j=euMxE?O>JbAGwTfl+gw
z6a6KRtsC=L5c<`xxBX~uJFxfXrC5^k<aLkR9G>53nI1shbKut@y4G0k?XTSEE^g)a
zLz}`0f<i0rgvxt>Vq|3pJot9Jm(I9%YxPVtu0<G;&DXvA*y`%zfEDbI`IAjL56$IX
zE)|}sOV$Bhj>)v-P)1+QGE;sbPD(K95+x^V_DaX7KIAp%nn0!&mi8ZY2;gXj9M@`)
z3JWwo>DPFnnsIvAK<}}*B=@KOvuzNm=%s>l)dddi8QKY<(ciSWOJKP*jX%B6&|ayF
zId67kqV`Zb@{n(LSOSXWjr`#iYc}z=Gp&vsH{)ehPrR(NI=7L-l?AmY7~UF|j0}8q
z+&lL?El#7_IYd=4r0&?hdC0!G?$I{|Z11zSWSg(JyaBkp0#rL{lxQBNy*z|dng~4X
z*&qVbgjwdbIor*c*v&iysImes6Iqpiu;PDF8Y)m4f^&MlQyMB37eRI0@CHaASjBcU
zeE0q8Mf=20tK~V+RQp0tt4Z`Z{`wJcnvq8gDX}%pvU@acunyP4-dq11-Q{TKx4Q>G
zr6!QiCNl5xX7>ErGj7oxYigwVTw1NapVcZVTh(wKuVvR`WjmX*YZTIBC4!mXABa$B
zf6Ch4rM=v1leH|BN|uGS5OW)y#9_iSG_{$lOI|KmtgW1`wH$K^^F0@p?`pJuoot;F
zp-QBpx|*FS#CX}6jTX<f-2@-c?N*wqpr5djpWi+wttS{0pG-BZ@7Y^h3&@<>+14PS
zl5yPjdAt3v&9#tBc~I(sBV!Lh6YOelG}wo2q*3J$Bzm?rYWEqsjK+iS8V5exxZSlQ
zjAh2we)dM&=Yv=nC42RXGwafO;mTi^+vqLe-hui*oaz7J^!{T3{IAn{sT?wQ4DQca
pThc=s_jXbHzs$+=J<^fW+EbZtoHJ9ewL4%Hw)XugW9g@#{s-~1q|pEX

literal 3884
zcmeHHX;72r8nwe%+A7i(m7TOMAPPh-MU*Yq+T!hq(3Z*~k`ySSLP$wSShBS)L=*^%
zfNY5u6e>$X2wNb8S_&aVNRcIh5TZheEJ-8*gJk)xclps7vG?blJJa{inK|b<&za|W
z=Y2sp-`Vo7W1^TcWq+4?#`O6y`Q5y$Kfhfya^1=-U}x4_CSNhj@v0$q6a69k{cTx$
zc%A$1)L4IbJf}Iv#ARpKftV8?Y>QDHz-*0u@ZJqk-{o!#PP}~B2|VIy8{-dp$AfB6
zRI*U=VKx?!qU01-Nst?Ml~+>fYm<ZgCSorTP(AyI=<c1AQn27GbinoRA|-L2_Ln0v
z>PY6#YD1F@MI(9)>2oqs{h)k@|D{jrOhFOS3n_J5Gq<{aLQS1B%ZNU9Jm{A03;G=h
z*BJ@COdzu-P2$KEj^+hW?SeW<twd(BN$lDRXVrf5d~GplP&W{v=BNqCHs8P**l@mi
z9*J){Xv*6=cSV+uq~WetI2map)v^;ZTza-|s9T<Sp?L{8Z(AIFx$x$9Jq-hO4ZLv2
zW;aXz?A#BD12u`fw!k<?A7qv<)UAkoi(Ei1Cg(d1#35O@23&>jSN>4SzOaEJ<d^=(
zckl{j2U%w+d%9yGU-^IGA3@7cka@G9NqNEfwz7-J!f6-k!Qhhkfm&WIKYqYfU8;^y
z2Oe{wx(W&e`GO)rp`Zku7tcq{vDDYppR4=S`_;bb9QI$K2Ve(bjxaaa$1rEuewZt4
zFSvkSL@%`D3;Cz!Z0Brc<#Xk-$hmv!Mm1dhmzS7y7LJ6g#g*XN*`IPwLY<|Dq)t*W
z<Ri!d$UcY*WDn#Z#1Y~K`559%FQ67si^KSW0l}bvC*Xqzz=L2O_|#HH4m*@^CaDNk
z2rH51NpA^x2p)p38c+?Ycq;yq%uH>nCaD8n?oBV$GxS;d`}$HnOMgv&PhZ52FrGCA
z7(<OgGkb_&;3MDwun%wn_5cR~N5BpE7;pym1Fnc7@%(Bp@d5EUNN)<f4%n-}W5~L>
z@yhE!X8wVud876swR`7g?f#wN_Q9Jq+8_9m-hZ=i$F@xd`_Hrj%QyQrAA8eX{13?v
zI6)L_9IZ57&euBuJffQVC;We{VCD!g2-8~vvg0fNgg)YHAiOd!-@OUHY4PvKwRE%1
z+WoHLwWvR}&9C6U>;<sj*rM$dU(~4X$i+u2H+0|@;tLxK>`O_#<8sKM7gO`tNs*8Y
zC*KrF<0}*%HChIP+Gt>(`}OA=HZ1BrgT7HTQSVDba-Oao>R&xbukRnR5hPrCT-Vuj
zNldr#NLU;yGjyC^rrkdsn<UI;Fgjam`m@KI7#oK}%O{dN&jfQ9hdzrcV}I{&AIV*J
zRo0)fSW8(p^!LwwW=0L?dP8D$X_5rt_!t_Es>%#XI~(RCROrlR=U4jI@3ryYidDoc
zolPYs+R^JLmlnrFX=UqnOU9YCU5wKc@Ordod=$J)<YqCORuviq^<a~yd{}X}>q;zC
z?(RWO6^->;i=|Rg=#p{XB-hkZXDghnwEDp!7GLF%J3irVZXt$0)dUTf0z(#_n1{pR
zLpIFTnu&({`(t&o9|Y_th6#IS|3qtzTVr=L$H2WLE%Rp7IObxBR+<EVU7dIPgf6f#
z9Ho!IkKlS5&#LO+3@Ijp5e06ZIe9FzV!;r+82fnUNmnoth_vNkdy<Y)7aYE(Qd5EG
z5vK*?N$uAZuW1$Kg!VvZ$!~r+%Fk1PDji5-nm&A30S(k<bgD)TrAkQ(_}OIdMrXUG
z;R!(-2R%Hk5~jd@(w*?yPG6vmj`z+IIJ!VdZ71>*RW!P$x20r>5Ta#M4&SAA;nIex
zxZ?(<G7Ak^E_h-b%u}-bFdb5jcE{3^V${&8thTp<idix1ff7_CXrkxY1J@0VM7lI1
zec9WBy)4;_Tc@Q<K)AW*M#)LvgxgYlKa<2Yw|=DUQi<?9k4$d7`GRrQge2=EPvZGX
z1O}#g(7U=C-HiuYqV40i^R)?Gkc`2m>HF1UT~)0)pXz;AuCN}$HW)gvG%aW%`Ijb&
zMW#FhFCqCQbr7FEpN@e_1yLHZ*fVMsV~L*$lIBX+a62?e;jC42M;}dw5%?5~8E&bm
zC4N(C0JPN6#h)n@B&R}De(4gi4$2otNsGmpbZi*M9<J9iYRISY+oSZmMLi}j7zzTn
zDK|D59lY$C3-2Oh4XlmZ66v80?3Vyd*kf2jWE#9&a=QVcj#DqZg6#~R8EKdLY{-2w
z#rELc9m>b8VQXX_HXtY3U-RpUV81j8IxLDAKL6XA7c6zO5|;L?>`HJX`+IlI9B1|M
zbUmVF*vFOdC`Sk}O!r~Y9}M&wyA*B4oMz%3JvHobN(mx;Yp?pfOxdpO{6kj7VTWDl
z--bgUG_f!L2)yE|U+1G=8So337wq#ihW@Dd`xA{*Uw7?bXAp=dx<!W_y|P)SJ!1;X
zT}NF#5h${iZFwL(Ls8@2dao_9K1;kfjR7avB~Eh0bkMxDhxk(eH>V|pB4<F{KRISa
zsiTL`E9;U+3Cdm7BY`=$Cy61g@Vd6Z@tg=jSZvLu`7vG4I32D_KCd4dt1M-B%trSQ
zW|K<|M?+#UW3{Ep?$}YB_Zbjotjd}g5{s<KT&|FFR(xae6MU4%ESif|{DfteTrczS
zIW<;gzZPH9_RuZAe4buu1M%RJoM9eo7rG$ijOOa|!SF77&+{2inzC1XJ%eh&AmeT^
z;ZphVSb3#aSsgP(7hPVoa7s?6BD;tBk6amb8kWPE%%R6)1aYXaEqm5g94nSG8;22%
zBFMu;AGgE}Y};5NoZI3mP7X0lwM2Q1PoCDD4>6l;Q5<0dnP;1xh?p145mrv7qN#7!
ze{QTXQp~u&n6ydWYF&@v)^F3B#d{A|iUM>#;^F6hJpzv;>)~rOyGQa|vULq!IVHb_
zpApON<DpY>7T$kaK7bFB$Ok0bmQ;%p*pf<%_a_YEpwxEI#o~>`lN;|M9Okv@MD5Pq
z-C0T)8*>{2nT_k!E3dJUx21+#N~TD9GQ(Bk!7$5RJ-3mpP2BzJGhnjQ2C9HjJf-i9
zwZvx;DnOBR5>szq7v?ihG~&MQjHGvwa!Y}_D_wgSJ>|zxMx#Lu*mJ=)4j>!os9&41
z8y!r;UML|asVe>AUgfZq-PZLdGf)PcR^3Sy;aN%mO>ENrC^21jfN71LDj~j4xrRop
zT#n*5Y}a%OMTR~@uSWjr)qkfw1afE3|GgI}>B5O<a^1?)8zfZbYO_S+i@J@0wdkrP
zUJA-ZZBVf7yw;%XR$4n~HRe{uQRmzl&v+E?&e<~`;(NHiebJk;a&zR}YBBxB(MS6n
z=GR(!tlVf<uScKlMcpK<JRKxSdi2>C2@y+3_MVy*6?JJgeCV^Z8Q_9u10z4=Wo#FB
zC-3^yPG)4#<P_t3DG2?is223hxDgw#YkfAOXR8v$PdO=(Nje!kJ#t35x1AII_-~MJ
q1s76WEM4u|!(aR)bhd+}SO%o0wdk6!jzmh|&N_JUeEIP$TmAtlZb6v<

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 5ba4cbeb8b653f25979ba59ec30c503867ad39a9..4ee8850cc4ab414d6d39b549f6f73cdd2d54168e 100644
GIT binary patch
delta 2344
zcmZXRdpy(oAIBAn$sxt!xXopfYGbWA8Dr&AJ!%}66|$rga%(Fc*S^eL!YHFSm`f@5
zYc2^34Wmq0(S_KgsMr>{&CWW%-+z7ndcR)J_vih7Js$5bFahk=5!pYGW{1pgW5E9}
zJSMjOVM?>uVgL-JX(N^6e~6vGtEe67AT;S(W4Fqik!5Q^Ta5c<Y^;--AI{|oO1hOU
zWZFXia1I7$>Z8AWTj2`3bztd)5P}v%`*FcREhODo%QU_r;T{u_$ZZDqHuuhpIUl=G
z1Yr3kpWl#ck$<&!Q#jsB=no07+H8yd(LI%>6SBDicy6|}g$G@UiuFx%AJj~@-`SBJ
z7Ao}p=KO3V-gr1iuzNKyG&#Z!Wor@O*?w-l<Jz~9zZovWv8O{nPZFdyv2xfm{<4rq
zB?l>2AQT12D5CN?cH}Z!&r^<%O#)+gclGZFc~dRDB+>zlRgp+0(zW;?<r?I^0z}>$
z`0Gx5Z4VtriGuJHAPU~O*C3U3vkFC?ax>VyQ(sIP^eD(>5+w7d<0|*OvO<NN+D?B_
zwzYF$*umCN+3V7~spIgr{fXaPLHG6C$tQ};`wF#C=XLAyI{EdUd1MV6m8`FZe_Css
z71hP4{2ss#BD*{L>s}Pjy;96%p8;G!dync-4;I?xxM-B-xjvcOw$FEDXQ*`6$mhjK
zixj9aLvf6P1qFUBR5Z+2_4DC0;Qitn5m3J(+fPG<C0z1yHfRA_Qw)Tl!|87)!iJw8
z-7dR2gkvvvT}?tQil!l$khY%A`9pC=00grn_3`cP^j2%ESe>gx>*Z4b7B?)8zc#K8
z_DVS2hqqB|Esj3qHG62?JY1FCD}m3z7ko8WD-bfv>}49>>^l$Xt{+bRtDOVUd&G>#
zFRx5H(Ts?L_IkK$HX^3UiFnDB_g!_`djLAR`XM@er^r|E2-C!F-QT!hmg_p{5Xj!7
z&OT;FL2$@7U@j2e7|rSO=};G?#68>gUcCUFR<nm4@J~6_=&gz<g=Wse_a{|-<ec1k
zXqKS?8(*vFQF!RPj<jm3Iift8=pGi&?Tk7sIM$sXrH`ABT@IouE$zhLtW3A>xNK6X
zH2*qvwrVp-Fh(dqR|Ab|A9IJi+T(?tOJhdtB5`nA9bmP_BeqvCcu0)R29sz-A|Ix6
zZ^#D=_6MUqU^W{iB++^Cw(Z=m73oKO%^v}aJ;J$&aw8bZZ;EJDTR>hpc`>?Hhm7(<
zE5H<c6z-%Ka4SBfFj`t)MLGggQBN7zGw4O>a3663ycfLEwwMijg)6)b0lYk*uzO;M
z31Ey3IjH?sA$`f)@?sHw&EDRMD5%vwYJi&_y@}@d71TPu63)S5m#~%@e3jP$4m^L$
zTLe9vXvb;=V~9|PA+MzPM60Z^dk|u!9P|v>Pmlm?=ij|NE6MB=4Qk|Tu%9?WVU9TM
z4a-snobaW&a5E!}4k@_MX2+0c0%6p3<4hKJ%?X$aPZViCBm8;KPRt2&g-MkgDS_i0
zsfqLXbLH_IAJYc$gI$m(D%X`Yt_Vnuw@yh{U0;&Qiec%iN0&7p%erkguVrm`sVUZA
zxXw35h8;SzW!M5dtn6`Hh&8}hwI(~ukL<`VyzGCCoE$WDdMI;m{^3_!0ej?bPiSO<
zK`EU2X<MI$lJaxC@*EzEI&ElHrg7lbE}75r8tgSwy}VscQTIa~CN^kp{>RB6u>W`I
zi%w)vdS<B_dBDUiY==3>WqRD&@nuj|Mqw8hm9VA0XKQ(ddF$M}3}Q^Ji;t%6gPlf6
zR9W6@<x$jd`PDY^v=mVe@K2MYau5PF#kF-+E#2QBijGIeN2I)0?dudepn|iZDYX%U
zyC^$nwvdrMXqzb6f#vEm<F1XbooLdOAngGuQ+KG`*_+1~<8L2zOt(hvn~~}yUQpMi
z;cu-;)1GaVTi!=fj#xvDl1OPuAz2(Rt(=5lVUu-GjCt{48^3n|G^yr3YeMJLrAXT<
zhx${>cB4=^>bgvkZ|FAY%~e(350T}zMskQT$pA83hPbWc_~8!YbQm}Zin7U+zDV6J
zchay70cY(v`<2S!NfiyUaLUCu$a~}38;i}4=f&1k?FkV`yNq$9R1H3F?<;=qm8S0M
zYmK~;Spv6iVFCBy)@tqKu)@mf5347=Y=|-Dk|k+ie75hJ{M>shQ&#z__~ZqAx_#g|
z_U?XsZ$opJOUVbau2b|(_6r>)nn^X@AUyteBF#km9An71d4)dsJeY}zsO`i6SA7>l
z=!}l&m9lq^X^92AGgb$l>Wp4cB5fwsu@pC^_+dDpnzW#!1PduZ7cPizj6E#9+CaZE
zxha0QbY$PI28>9iiC%s8dAIXeZT6M$i2AQ_P|T4xv!@F51J9k<bQ_z0(1^nz-qbnR
zXI{2ldP32h)DI-32<cYbOb^e~Est7JkC{1NQfuD*=nCYloHxnHOOFuWj6ijhTkxae
zD7OS)ezo$h?aKGrP;`z~pJG&KE#P;9QsQ3B(sV0$aov7TR!iT@>}vjCr2QmnNu)aI
zfP4?GA#5xQ3B5lSqWROg3zI<CQ%j7?r#)s75t{=#HGqy*TLY7&Xmb!ga*S7+kup->
z5|CYfP18<%DF(4`-9ZyLR_69dC@i-va|2v@bKjWZ`VUkmj+^0V9-f1REh>`Gt5cVH
zjGJ33QH%_UTF6~b9<Z$pS#=KPCo5qPrynzC-ZZFA@YL?U?RSfi9AOgCG@biydgAMh
zq<yBI!I|W#F@UhD7SF2YH3=Z;Pt1joQozlxC7>p_hQYF&_#pwnG|4@9N%ALhaz2E_
z6wh4oB2NAU(`(*;TeKNBhkc0@oeJiU{^w0!_K*nT>Z}_ei3)1^-<H)g{PzwIYARPo
zej9ER)x1ufJI5H@%gOgqd)b^mp7!}4r_#toc$XMWnIUPq%|aKnhEv7DTh}9E^NE`M
zk%=wLNvmhhFHZ_mohK}K>V#0y+%oEu_xMHwPkz6m6%wfn|M_%-bwU3XaI?b*^Itp^
i68@`aj-Z_eTd4lhjYJ`kG&<N;9$^MDG_=Lpf&K++HjnE7

delta 2342
zcmZWjeLNG08&7uLg%BET(m08E-|{vS3Zcu(64`NytGw$BYjY?xjFxlgD5=|VxHz2h
z@=I+RHS#iP-r6v2kzumg%jrJ9|9bxVzC7RW^Q=QZLQ9>#fjc<bS(hR7{tv=u_re`W
zTB^FFGjbYuzv0Ms-5Z>_7WfJM&H0@emCF{!-xV$eIcJxQX6uv%mxhMOu78favhzzf
zxK6F$+!Ri7HIX^ARoK_tD;9`E!VgRy$Ko^dVI$S~p~a)5$)t7+CrPr-73~PN(%Xfs
z8MI8iD3~Car8{$)CBu4XiI=oviPsl6Ct;^d;t{oCLA&^~h;}8H2NW+gEYkUzCjENa
z2v?!=XLHs{zs`q%j9b`FC$F9d)paut&9n4gBE4sGwZhkF<c|`mL~c#dSSjkEZr;c*
zPD<~=dAXWgqtHwDe{7G{+UMb-#0PsC|M{v;lkt{~EKoEy$~%%*9(L<WgQrnmn&v1H
zf6C)rWsLEkk2$%e8Xj&+VM@$;)u>`uP+p4WAoAPys}^mC^BlJRBZ~J`<SKQ4_f@jp
z`YEeL1pw8aZ1V+ozj8mRqgLKd);XnX?g}aQNG_G*>FJk8KD|7TsdNJc??3Y9kt<9t
z<7K^%gB-+>QNsiX8b;l3r)p2B2jj0H3w6%iJYBl?@*Nef0yllU=IPI#b&r3pYECL!
zu`+tI`j_%=rbLGByEAq2_F=7RK$vAzt7Z)H5BkCEt&tpyIFPD+Mx7m|#wC949!>i~
zoIFyMMH%<V2{`Xo{KHmvt&2u4tY#9~=vcFv&nGKlHuEX-6W{Dli@~nj;z0N!y@Y&m
z4!^D#H3K%GV`IX~Z^D@aloP#gF#CUZjQl&Z5c+}+7olR6K7xc48qc0e`$0y-eAptl
zIC<t#-~CTx=uoZ87?qIjC8c2sue>9mbe4uhjp6VmqVP49Fm3GygmOF1#&_4Np1;PJ
zx&Jng|3~6qZV~7yOd)d&ug8Y0YNLxPCVmZ8^gqphM6${%0tt?hIMy}M-LL$)1pbq%
z=0j!yW*%eFI4<NY$cM0jF3yj6^bSMehD~n2MThX*ub+V?-{}`F7r@r*T_;mb0y(*v
z#TpwX&HKYZ3{*TN8CwK|hP%DSge9ug?JOiq$It%?HKaFo9>H&drf};=MLeosN9_~-
z?bm%GzsL>k7idF{JP!e9LqvW_q*y*}!==CNz+OpK-1~2qg6<VpEZ=z+c!+Ujjg&9Z
zb#dG;iruOAekf!uE6rlcrTX1`mCspgsaOLcp%1=Fyonu0pR299XQr{^3Sdu3)t_O?
zcu;cRTihi)-d~j<&WNw<F5!E!r*Qj{i#7UCc-lu0)^&Lqq$Q5)CU?}m!yt*7J^BXF
zBm{t8@dEqUEl1?%-g9u>hsui}8XG{vuHWY8L_Q9NSsoJ&Y*Ww?_YpzBVKL$&+bbL~
zg;|<IiXYmd393ZgtO2a16Fa!z#ibGSi6Xq}#v(@5QGsq()S^Gk#Okm!8c~%l%^-nZ
zO**_JciZps_&jQTL3{(0@00=8@PGtDRptx(Fg-ID+v6wtNurk>(KvON<U&v7h?PwL
zQ01yqivLvkjT5xHZR{g%ZSDi4$@?C$iw74jcyRZ8c&O@{uH)XywjK|?JNSYMz=u=?
zECWI=Tt2GY(2+NOw%hmasqc&+##P3ON;}R^<XWlP!NzTF*A-y5DhpH#imtjKmj^X8
z@7jvwRN{_#kCxWQ>fG#D*mHbPu623OGtz@hrG%DqESv8w$rsbIzu#k{R=QCNcumbA
z6B<|Fhc>YQQQ2kp%DnG_9S@Y&lJ?^#;G+$xyUa9B&DdXabY$W;m^pZQA-R6Qb41_c
zFnQ!m^kMRvRGR|b4p>N>u<B<mDfb&L%JoMQGja?sb_^U&D6v6h5vvTr(|*wgVDBM-
z&B2e9m&RnDdk58jyrb+vkIDHqu=F6{_Eq&!=FpuJWcj=?1VcGGE&H#i2l(dB?$Od7
z8pj+NpHrGz4Yp5^&yxTt2HV;8l`F;P%e;a22AdTMigy6_W3Sz`^hs853$}DZ{wi^E
zR}W$gYZcn9nHN5dPGOIEv4=9UPw&&bgy@p%KkLI)!{yv=+`8y++?|}duAN7&1rL5d
z%owE7W}8eE$?HgjdwHJ`9!&Tt%gC;Qo_%?KwX%QD`{jVB>|>cZ;|457$a7%Ept1*q
zy(bDsg3AvO?1VRPxuU-lPq#K)yU&tmR9kd3o1PGk30OVWA?6q@-{3jJHevP9!Y!VW
z&uo%OLJH*pPw}@D!lw_fK8A{iVu9kJbPWaB>LT-ZsPxFmHGvra9Afc&uQ>O~yzNn{
zDe5H$HW`~zu`SNTNj3*!HziC2qmlUPwWwXZ@=4Y{dom_^3ylYa2Tq!7tl3;WcN&g!
zo!asup^a8{9t~?9Fo)owMFYLpcPihMwESe^+g>Mn?{0O;<Z23GD7o`mLJGmI(g5AH
zB0uN7li#G2BE^4@m~v-SwWZ{^J|S~<f^8ku`f@rehV6TUZD3X@cWZNsK4VjQqJ%m<
zUfcb*z?J;T5Bq#P4c&8HFhrC+?UjA3M&Wm)8$|t1PCmfsgm)lRCaCR{9WuhU%G2bm
zVqZ?PQ*Li|93L@wbkykNVvBIBXwbI8j8$;9DUQRz2NKRUkxO9hi}82DBZjLM<4rsN
zG1N+i?7m89%P)UIvup-2Lf1_@w}1F_A-T&;V2K~Jw3@8$T@37^EQYq~&Mq9N7<v(!
zNwR3319q^gr&+eLaG`%fau*>{0PCWFYAr3AL2fANZgIr9PSNfHW#x;}$mAua;EA7m
zb{z~aOR^iWscMp%qQodVUAF=_6pAN>M_3&U|Ken9Cw$2S&Fy9goZiW~-K(-0?k@3r
z*72p2_({Q{7f9PI*OmQ~YLzHE>hPfHE%N5K6S+Z}Z!`jsT}ufUs<COlw$f+j=IO7R
zvPi^$@2H3aa4&)LTLxnBLC`P!kjlhbsYW>PiGLec`Ix$!{c*7MmjNe$RsI9kwvKkT
mUuiIV`>!<^j1;Wrtozlr1I+d~sb0@Z+16170`Us;R{0lG_Lf`#