From f509bbd15dc95eba2118d9794889231686b69cf2 Mon Sep 17 00:00:00 2001 From: caes Date: Mon, 27 Mar 2017 19:52:07 -0400 Subject: [PATCH 1/4] finished lab 4 --- accelerator/data/nuclear_reactions.ods | Bin 13789 -> 14896 bytes accelerator/data/plots.R | 16 +++-- accelerator/report/nuclear.pdf | Bin 5365 -> 5530 bytes accelerator/report/report.tex | 85 ++++++++++++++++++++----- adv_lab.bib | 22 +++++++ 5 files changed, 100 insertions(+), 23 deletions(-) diff --git a/accelerator/data/nuclear_reactions.ods b/accelerator/data/nuclear_reactions.ods index 0a13501292a2c221b6ea87faeeb4c8ed8f77571a..0bff1741dffeb00311e8a91edd060362c033d414 100644 GIT binary patch delta 11229 zcmZvC1ymhD(k>Dp!3iGR-66QUYjAgWCuAVO9WL$?+}&M+yXHc0=LUCp$?m^zcmMbL zoH<`tcTLN5b~gGf+vq<3#}ZZF;Ia!zO-e|-6W!0AidVx`>z>D@g)LEI`7CDJb# zWFhNn2o$`11}_%pXznn?Qp;}>$U7;+mlY1xACe3B47;ou_Il)Q)gGKLBBDmD@K%WC z-TNrHg#%+o>)+ShHN(DveO@CVa3WllT)^bWq)(E3N~xfz<77TL`!4r0NPF4YKlMre-R6IlQq_H)M}h8!GJ$wv;K$3?xG1J03>#{%Uz3L zRl((Q9JusT%4w->xAC?Q-rhmBf}cRJLH-SIKp~2r5RXARRP~cZ3g|O9Wupg`-;i@q zReYJ@dhiX`uiGB557v{cG#mcQfNn&+&OZA4mk2Igxmxka!W{Ms^H<~(#{jOqX$2L= zkH(SFZBkgNfJ1*Cki)SGyb!EBH}iBX*j-MoVRB0ltSy&BZav`uP3%_Fw-SCIoNQvwrA_cXIXlj_2vziC_QiL zvT#l&uTcA|g6rPAVs~ig%{I2vqs-AmYKC|}DEZv_D7L=({VXPD5y>1p;ix*duE0FY z0OED*hxr00Ori0vz>eJSfKObOhua~M3Yx(0czG|}in0;_qg~N7@0K~Xc8yj$be!bk z0%H-$kgf4m^27m!#F%BgHhb3$RgT?3YF>qSVS&U#OB@xt-MhlxeaQf&(%qY~v9c*$ zw=k?YN?--Rf3%Y(@BFP;O^>E)1w4dsf9Y<7%v`=+-(p&Knsj3VGP=_78-#OiGUnFK zfkYc%{{3Yve-|oekmRwQrI=MEi#T0F*eQ$};_9mvxP4pby^!2jzurqT+q}Dk)^h;Oam!TridF+iv=rx0wex^8nDb!os)=o*-0sP+Vfe8$E?o13V$_Oxm?EhszS zueEyjlfjj0r0A0z(Dv9=2(!Q_NvUzXi||03+aa4+&NSVe$bf(!*QMbf)mhkHm;2PFULKA@rpX{K$>0PLtQ;-%_?>YU_d+>I_PUY%wMM)5H3M_ zJ7|agh=l4JK+#TbcE0p%oa6P`sySwK+Sha#aZ;>nN3WHIZb9pY-Z)a5IE#=?B43n; zlaXO+p}v1n;~0&QGe4`q!&rk^yKU8^e1-d$YKaxSZ1FOmW$U)LLR`WTTu{!gl8W{j zETkLDu-haB5JL#`*!sBH7vmsRv|?G~)331QkDlo(vK8qCx5k+)&{TauTBQw(n(;0D z)Oqmv;@8vM@GoPAfH=z%TrT7$+h`^oI>soMOGCseJJwXAxtsjakmrah3o6JBxqQ<6 z_LY%^`K-6w+$Kg5#YJ%f#P7vcU|5~|nR{h;W1i|2P`nvp*QP3#oAzyR+Uei)%Bd;- zkkIa>?s_A3k<^FcF|8YXgd3EVFi9ea!J?LQj1@3iMf0 zgc+05K=&3pmf_zoH)yvlNn<1+eKmidTVZ%e4+JMg?%x-($55@}%%{_{wZH$Y1^2)8FI7nsePF6n4qMxg#HN36A;IDDpwGe?hY%QK>t& zOV`DE-7<*Cwll0N6opX_g*2a@6^W7ZsuO|se!&4pK#Z$QugfY+sy zO;GKUVQ!pa%Kl|u6g$p;$uB3W(t`z7U5$8HYysP{@MlI^1^RknWTaNXjZ&J5!>UzVh~61=YLM8RU;pv z6qr9Nn5e6FuoPeDJz~lFVQ97t>wc-GnqZkl1^n#`A;TSuV8gxc<+E$6*pu4r&B4iK z>?6}c?A_D|`J$5BBtRm|B$1#OJbq&|kezrB^;pGuEH1Fe$>>$;?NCr_G$2&v5t!c; z>g=t@aV!@cs#`O%Hd&5)gWH%E@=D)ukvh08ReV+0Fgyv4E$;lM%Wc<0v`84n!GFBd zL@mFEA)6zRMa1=lWEy0wJHwXCUhpREYjU>oO1yRdI1cPwEv?&oK-3;hJ~_BPbixU= z%^7{#-gt~hGR6bfE!OE}ffnbj1c_}DCQnVYpz~2t)D!hu%RTR0?ku;J2;bVNkxP@} zZ7l#r#1Ii5+1Za)pb7$?=&Qn(uCIE=HJfAb42XkCH{JQ6nLjEi21hR#v&e@K$3$jW zsyztowt#P{vV)0ld`2kwC^#kE;D)fDlGeXD_`qdf=W>5FB|J%sLnTS{;E|l=N_{De zLqkjzc(&@;&*oD0s7dM&*=uE4`&Yy9iMhE+e!h?_PQ=8D;A?QY4$LZ6=x|bR6K@UPYz&8|SCJ3-OaV+5H7zNL>tb zll>sJ?*vYAu~Bo|D3M(Ijd1~j1?q5-0cz@bi^`Ex0F|Y7u0p%Fzk5y6;ffm(g#&Ae zOqigP`29ksd9hAjv~8N2dr+kI{xHGK;ar4=Ty^!gJr@m)*^_mx?#pMFx!Yl%J9d9y zYYvzLwKt|M>CnMNI##exMVO$t@^$#ENCA z{diZACBNl0wf(AoUXiIAH6`*3CM7=INNr7{#zs zMA1mN+Yb=D(oDaFcC`lkHL$|!QgCpN2?4ho!dOO(n=e=&_CAz+rhg%a=Ugx! z7|__1AGblk|2a*?D~gu(lTM>(bc4mb2vjI!0%Ifsjo}j~v>_^zIpU}wqc=JVOyc2A zH>Xq5=xTm0(}@&76GQk)AHgDBv4*#1cA+7;@{*wwD)&H_9;>H5!4mE+8rnQV6qC`-58@`& z7iTjEO~u&pDB(wUP8%LXO7K9s!bmJUN_9=N(=q~!*G5k7dsApfqgAj|xP1)2U&(rj z+Q)7F1>&+!^P$<|GjZYu30>g#(h?oOIF=#`;+^Ik++k&G5HifEV0x$slwmO8_T-B)~^qNG`Bs2c|Di85&KF zuBR`#nzA%%C-d~>Y8jBMwt-#8q3wgm+$5P4s_>G7G_>e#QfNjn_msI&!yCM;oS>%1 z_{}ubq53GIpyhFIlO~ybdwk=%_xEzV!t>9mWE#dU+i>wq zE|+ui1hFK#KOfeXQ)mD<Q+H1lNLT}=)qSA}#QA8}>D zKkDiL=9Zlk`~km*6`7t0PLiGuhAmWnRx5Np!Xz8?Xl~Mg`bQC|-y_x{z-Qp1DP-^y z60gj0B?upx?+?|#1iDtIpVs#6-StFoMN;8xVHk~9{+V?6i)EiSMnZ~I5%Nz5hif~l zGq9`E&4>13J4gTPlRb6$^F`t5PkwY9UD`wd$Bh=Gp|}q@OzWWAm(Xv`uWLr2-qW0p zRh@(vSYJP+KMmv79HQG$z5!x8^Q&Zs2`!7AB7GinWY<^?yeN28jx`dKC(1G#6DZ=f zoY-iC$d9Pc=kCchn~Q$Xc~8<^9zEYTwewdIh#Y4)nhLrLj;QPTX|U39ZMZm_pJo8A zQi_3G+IyVCNEV3&yygb89qf3nqKe`|)52N7M+S<_(p-%bx)a{FwezWYSOf$eRkyAH z8&BE09n~%sNMqhx$Sm!myc^NLt~8m}70NV3P=Dh<=4@bbp% zIZkV5C+B5%c^S;;C)Tc^huOlJ#a2M1b1=6Sclya)PR)mveMVNU+;iUQD9eb6<&^Dx zB=FlhiYvk%?}>{N2pSR^80#u|2>MpSMERCWJ@bXMZb@`i=gh``+gva#nV<=8o;F++ zE)j)_MQ}>*Qdh6j^Xcqxzo3t{pr)uPyLO@yW$K3bQO1tQW@LQT^ZS?d@?)UEOw0qt zo4?lF9J+n2=R*U zJrcBRonR_=8k5OsQt!$+ajI7@SrVwsP<;UF=7&5>VG58}#;`XVvOjx$-f&=hB5v() zF!bPa%eWk4o$>e6>e;B`rv+?fu;tZnLGrHhT0r_gpm_8SM8GvnS1U_V;yWMgvZa|!O$hG8#M*DMshZ1=aw53Z) zhKjfFW8pE$(P8bz4!;-W2HC@m%~?vUhZ{d-v-#PqLGbHAL8b%GB7ov{lls#pfn{w2 z7uAA)DBZTZ?b-NM?bF+&eEX?)M&GrZ+UU}~h82eit9k3~zIAYF?r2jaQ!Rqme1BLr zg|<4wag3p|d!F5!nF$&DMt9Xx5Z}NM>(H%tE$8v?dzk{gCcmbZRio=9YC9q&uK1ja z*w95wt)@1#A-^0$;D8E`#x0q59utcaMfXP#gsij^3bgEP|M6z|l}fQnR8F8=(~%K% zlSQ}71e>yzH>Xz)EyNW5%`)$z$nqiG8#;ASwt2GqoCf_mku72`m07KiRxzEA&P}wI z>3M{?K*M_Pl2xubo@0($S8RIh;E)p7{f~zdbQ=`@ZT_w>0^oy3vo2)#)X=u4m|Afa zyf_VwIDS5)D*#+_cwic)^k1)Qxirj!#y5L|cXQ`NqO+~fe|s2|uj@ofcoEfsC8oz?_1wvMJN z&N2Vm<&ERaY{a`8?MP{2yYpWn521GFPnF8?xNGCf-}{J%hv_K(ZzS@+|49ECrTpo8 zo!ftf3L(J&Inyv>jQ`;zp=2eMB&x)Wzy6cJeEat89|9AUg~bR^@t$qH`-#p>XEWE1 zOl{-05!4w#UAIW_1O04CtH~vV$J5mQxG3w8rY9lWK(l?S>u&23G^I;`o# zfH}*TzPGgoPgSqwR7@2ZF*~<@XI>_5%R|#|q@pcs9!VCOEX-m=t5k6&2>(92H-KQn zurvha?K>{?(_t9sOQ2DXx){!D)TFQB&pP$=s}mFytv$M$oSNc4w{2UH?D}3_m!AL* z%|}1iQE?|#rN&u=cuw*Ca4K)qy6EBbw@X|ldwwj4>V8ufsbwjRjS3QO4lxeM#0!L2T#T5VnP9~+}q_47qjg#MmPEp!& zZ2Of1s}@_!ZriP%(Jao_gx;+t`u%=msEK2L=|m?~@GO5aYK^=`z*118A=Vzi6+xKq zl~{yfEj7hEgX5z-lOmnT6>yaFyGNFmGJ-!*;}GhOVQSTC!SRCXVxx5g+uu3j>x)R$ zTQ?Y~1(E|mbn(Rwi|aFCP<|ENRn$#saC+2?_%q6{SOZCJ^n$6kM%byvMkoRc%_i$| zCRUQ))i9fNUX)Z>ROzk=-n_&dOmWE-5EBIoq2pMpno0VB%RZ7agpFJ)UWrONTSc!K z!A)G;z{EpuF3N9xc{i)`+lc1&q7o+WU%kIPbVn|D z)XMY}R#^~3-bgfx#J4lO1Hrv$N{xG=;~lfa-tM#C8L6gyt-=X{kV%(dn$xmRBNw0_ zT_Y4fzl4gy-rDxitCKe#>b|B020c6)JY?f26B-_VRnw)U ztn4HOhE2^7?Sd5?QMV!+_obFBsAWutgnbuL?94pGan!R3as=|uko z?3iHLUQELqT11**5UwFRszhAQGHZ1dkV~mnM5vPA60!ZnDrO3wijNFW!W@+Sa_7T- zrXSApF8_CKq~iQ(*mP`&sJF}X^YEit|KZEsiL!FR4DD}80Pns&{H=O_Y#-JvD&QO{ zYqU%3ZsoEZ+d(PwkjIT7D}4PpqY>m*^(DS7XQ6PN@-U_Lqo0)ij?&=1MvBgCvUjK3 zb2`Ts`7P^Bl}8#C6JFzpP=tG$FpZfRR^S4SItz;C_d5a9z74dt6dS(NTi>Onr9CdL zg1h`sgQZi8R;0jZIIufTEJs=127qZs#vt#zJkZH|D)BwL>T^O9tr|~8W~BBSr{NLU zdclVAJt?$Dq1h3(#tz;6(n4^Pmw`bATRm_%xcVg-BgG&X_4>~f*@4N@sXy1H1VLu6 z7q%*v?~b5p^N?LlD1G8dwJk3#iF{JUVv#s0#_a(zyA`gI~q;kYiO0d6H*4Paf!93XdjL zP-D^8rB05n(NtfL9OR6c3pSB@=*88IWHfJApxM686W^PZW!c)SNx&M&6MRcC9UfgS zqW{{8rRQ{u(42ySy?s~-Amd$APUOprxv_i$KG)*qNX;-3ma*nWkqIeF5{@uBW!7HH zQ(*~u6FR7v3l(9@(ZW55MD2G1d1@AHs^$vso#CAytXEl_Q--6xFLbJT4i9FW8+5ni zif=(dMV}lVUmTS8@LiXK0?DdNp>SnxR|pA8wBxRDSfm0lqT=@fBs1D4my46qi1WxL zqHD%op@VFAR%7$M(~e8A0~_?P{qG>-r`_B~EvFAu!WGk{DPW~iYgs4#Fmz<-P4ZUM zgzpk;C~71yM@hH{X}8yoxTi8tYmcn$1{E;TF)~I88VXJEhn#K5vy<^H7h(_?@5ME{ z;2P}1FW8$b94t=&(?KTO7eWuih-(5tTAPFXa2kWQTsy}v%(Noh{@N!E8$ML)e| zp&uKZBF8ljKDU4Gxt}?hGKA%u88#RJPogq1jQF+nUJQ!5h=QRTF!-fI9%Y>lkLhHO)jH=* zeKew&4$8=AT1)s@s^tK=AA1*|gF&DU5yQKD%YI?G?JimRRm|3s4LRX(`7vP(ltPVx zK28uN%Ndj~R2!eYf%R$gl+craTk_>c9~kb_76Plln`v0&>mrbq$o_}|tk@k)4LWD+ zhf!)#YYp=~z@A+zh6PvUnziMv+Q?v`K3-zJiW0&8N@M~q?+TR1ZZwBx62j?zCXeFL7?@Jj%}I+>Sy0yWj7}dj8oye~3kQT|e5qWQ zlFcm=f-eMGSGGDE$rz`2Vnd9Q%{a!RI!s*YBpyym0oRY=pDnlTDY@Rq%aLR}Zq|*f zqN7IFdNB`8TT?h?nipwltZZb8+dL{IqkSm%lH~ZotiClSZVqRIV{@tYX^wt&0BVGx zLjx+V&t8dL;}pbYX3?KCnAF*ouY@nEHSq1Ci#oVb*L#R_SZ3XuI_hc+D?)h|uVU?V zyw_5a7{H=*E1IJ2L&~|i2rceh*vhC@8F1-*aIcuvtGz$_0Ii=l^T3(m&Kc9upf}2U z(1w)}Rd3tQ5EH!OwyK%&pz!(NK1XIMW88M39hUo+;F@=K!I7e3fzSKcv$TTH=2Of~ ztFY$%S6?!kB-UI)_X5>dD70+SK3h(CrRuql_<&70(f-fK{_p$BLh&JBPXs}RPOQ>q zf?KMhUhmn)PT}<||n$17f~r;LJPb zTdWh$n9J<1m6KB&%F4Ylhi$0h`pVh7*Ix5ltbam;NG>MM;`QV8i&OD#=e^kArIXpY zkH4xBO*X!Dx5K?ORh*@0+7UDy(@p7Bo&u8Pja0AUbD^*gGJ=0x_esesHhGUy@+EFE78y@|##V>FNPb1N}p#H4DQZ+Z58`}l0Yvh?*}5oJ5POEmkejYUSLc? zQATK4#?{&*)n$3SDXi)zTyo{$m$&cDZqDD`AyLNq{yL%gb#gf$@I;X$Tryf+8)}k) z+RT=^N}gQMo}p@7=f&mV!L=_W($d+w@qBUP?|3;p+7S20at+w73x-35yZusE@UyQ_ z97sXra{6}&=Ru?sOiSAcp{S(3_FlIW|CO39}ZTccde%~G{FwA>JEe$2Z8M_bwm} zQWe2B_*8Png?H^V9kQn^hnbD0^JJwtK~gBuhQ4LXWjL&XtB`*3rBjN57pR5GWHSh<%3#Sy=f0+0U24lP zwVi{?!epIvkurjjP8JmpOb|!F9W~9zeVjP~S`AYh_;T0aPxpKkquA*sBTHHfY)pN; zR3om8p|++AFmpuT%JfA)RROv5IpPb2h0OtFsUopKLN+Tgc6MKL$k%RE&hpgA|DlvE z{;4XL-biGjaMU9+To(^x6r+wY^ccksu83M}L7wg&X2T^I($eSEm`GNHf-zDHRjRwZ zAR#MR6q3l6t^)5&lMtrIY-qf+@RN$sm5c4EQ4Wih-7QTOBcQ{?$EL{sr@}Z%f3VgN1 z`MEnrrxPhXeRK&(O*-8yW!OrU!Azi?ZD%!jc%14+T^CBcr1~nHyeIDV%b%7*!*s{o z=i+QZ8R=;DTe_D_I!(m+4BB`29?2pzjN6s!j*OoRthc8qP#O4w6pv;5r=Im3ZH0;=&{tXBTA+I zR_0IuF*dkq*>{bE5nQpASS_9n$K+FkJiej|@n`?gvpZc$Llws|k%Gp=1jh;_E38%D z{|wz)6FWEQET`^jB<+~+AOo~T%#UHn)uC0hjMSBNE*1hxVzk8SK9WxQQ%|YIN;!tx zGKBHOZ1u?*65jHbEg3bMO`X4o7m8OV03`B}+1-X}k;vYgoCdW0T79_-h{`)^p(^$0 zLX0at?VAtKCbWF(d#o#2^oaM61yLm$Qcb`M|Gn#By0P|g!Rg}5qNW%AYA4xA@gJ}D zWgUW&;ozmR`aJFEMO44*x9;hxHOm<@ORJ#kfWZ z&^+@N-=*Ch*}Z^jGvjw8+STMGP&vB~?Oofaf%W;EvDI~7R>h9)iszvp`F>|vlgQ$a zf=%fXC;v3X!m<)lr`Z6uE~g9qL*qpS0t@lH4&hi8w!urk#mvWf1nzT5U@If~s?-|i z)Br?aJ0I!d7aY0T$`)h$rB%zz4nm9z-p+@7z>O#5Lx&@pKZ}-AA}wWGbp@fz(I#v> zkXyDsY|krBN)46T?#6Jd#iR`sWtuJ3;fYHaH4{XklxRvwO62FKmo1~a=ttQfq4#+c z8T?saAErcJKA-e3eS;G0?{7VTU7 zaK`e+;C@{31xt0pqqEuu+_%2@kXDJf!^Nz`3atiRwda_8q<9yh;5Fzeu~IdukC~?N zl1^*!Gnsn;9qa}xUnR)CUWxGb4t(-F)f`&XSY~(cy5QN(=)Vq}amP-rb%xyD-wa3F zI=lptcimkuIc2T2FY^gwlIS|XM+rfzoP?nf-E$}iwnladOLwJf-8SrnHm%yhjqIe~ zLS#q5Z#Pqu_k2C2S_#OI6%@;SNaUp22-sR|T7B(;@OP{1P000>5hedL#WCu-DY6uU zqe+v?yHV5UJ_5C8{dH#ZENOfp2EL=9g0d7J_7*A@<(O`z1R9!(^=z2HMdKj7#_06V!&i)iu%H;_?cFL%>M@X8?OGY z$gglkg@OWw|980RC-;{I`&apszz5}E;cmk0?O-1bl8;>@#R|TC#>U~$BlyiB6(fC= ze`H#wX)ATa?D|V71V{Xhcjn9ELz-ogj{9N17tcOS74c*xq9M@BXLE%!6P63b$t`xY zjU#v}6ijfX0+ftj*6;KLZ5kl0hBDnz<5jHI`q)JBm8aQcx+`!q@$B((&`g&H8_d}w zrYF@WJS^dKdIP3-xBev8vgRikb9Id84Q=^fjat;L79`%qi1w!bVvVLXT&dlBP5!WQ zEH{)YKOHtzSPZnJ#$iouF!(o!mFU4D+kXuWiFWGKdiP9(n0#K1Dj;?3+h!KsmOko? z6ipK5%qmWVsT)N3ZrfDaHQ7-+S&CQjy_|JN>2DNiDjq+lLaH&^u~AwUPzCKQdFk(n ziBi(S>3G5a=8VK-PGG~**8t>=OV(f+S<}s(h zB>9J3)z7cA9$xFO%R7$(C{+Q^Khkr01dF~?-IDN^O5m}u0bZGy%v1CCd4-n=2%xi! zZ(rjwlpZ$^W*A*nj`Fcr^3(hAnvULVw)Htxz1w)8IuibSu%!IRySdPR4i)y_2m9tP zjamPJItm`<51qNrNDdlheE-iX=f7IMa3F$2nBM*s<(PQ>LJ-rZzp%#i@t;5rD(6H- z{C8p}NS^s0J>Qu5|0vY=FbhKe#a9;dzX-Bm7t1!d;s=x*WY{?{&5L$COi}r3(=oDk;;xrXhIJK<>u~VZ{hZzll(KK6{N&L2m@gS{bHy7yVJiC z6&vV`kr4TxOZg8iJvNX%2lYR;{qKc9K@t9^02R*$+T_FpIWVF8&j9~JJg@&5c#IA7 z#KHA<8b3}#7;<)yDktr~Bwcn;+P|d)_CKTiqYsIT5N4DeM90PSS5updPym7B&(=*G sZ7eO^+?id?E&uoS<%e)q$o^Q36N(B*#2o4Sn4o=7=9!m>b$i+|j%#q(7^@)fY zZ?A|B4-b!tiAhLE7~c*CFcysP;NS?2m0n8gcrP4eNBL$6)dli0c?zO|$B~!!?^j*R zUp6PEM;u&^`r7#GZPK+!6G*n|BC%vR=l0xn$SnsB=2sLx$ctJ z@Px$?qrD#tI$0@XmJZgeB|8|wd2W0|J}3KGY8$nn$LSP7pFxVra;JTu@VkK!cN{XW zKeCZ0JNyrcBk=+T$AixhGx+p}7!ixv89K@Np;pOaoFyJ8OAT_qfrU=E>e34e=;ywP zw0j1pL*O~FD;Gf6)@;41HkYx(Fuk(Z=wRo_QIBH;rnUFP#o1#Bgmtfn@l#if>IP~b zx(-NpQVp$k`w7M6qJ);0p={Bnw5!H_&{g^PpD9C*iun+ z*8jE2Uncq=1o+2#4_%{FT`&B!TA8f&>NWB*TgtCJ=9=j)3XdKzdib`jjJk0ArtST- zSlOA{g0}W++G-dAU_N=V?!)ypP|8axf^#~V2t2n_#)fJpo6!Ro?M{u$tEP!jjrI$= zydeuG9Ad(e`|%&tHV}qh=U3?gs^l&Ms~O|#d4MO49&;CQ_0)qBjC=r7IRkb6HE-$+ z+(kEHrg27xsV4Hh63maNB!w( z`UmLf14&?yd1Q+9+V*iZ!Rf|{hxi|U^HD5vk=$0BmXCawsbMdD*NXL$sT)QgT;Ayi z09~Gk}dxt27pcg)=-U9xpEw-{(LEt*7k8dV@#fAEW z_OX;8<~yOO(8Ll2R_!HR!TyT*Y$z!?KrAQvUR>#?dXFAWx*_wgwLKDEqIyew&B`6f zFS^)MMy-lgsYP>EiUNgOagTd7o%p~A-y#0ShFAm@}O z+5j_YC&VYHuS$yY)Lmk@P1e)@6Bv^_2yL@g2z%Z8n>wzGrT%BGL*FdeM#K@Utx#-< zvep8giRR-Z$7(9s5T+k~4f1eo3bK666nc21&sK|Xok1oB%Jn44<+u^5{7vxZa*mk{ z3q-1SzMrxhA}gVrqL=+hETt*ZKui1)Qs~Xf$&K6g7LzAS1E~R{C*cgWfsG}F2THkM z^ovoSo>JtOokE$l>K3CG_ZTE!A{k6;MM^Yr^~Ja}MkO6kM2)E4@qxw5ooe*#E50k6 zS>sKV=r-ALlR+jQDo(==A`cr?RxYE4<+P$Efs}a7jj!0@3S=<@5amv4U^o1PqR8pi zx$df~MCBn9mdr4hOu;I>ZvrIX7Oz5kQ{E^ZuBd$<$Ka=OQJ`!fp#jUVV2~^lAtkv? z{PejLU4DnqhM7=U;^B}wVTReur#VPw6}X!zfHyV3oK3lCA*0maKy`&N=Q$+PY8Q38 zqBj}GED{$@@kNX5HXq#zKu|JzuaBh-X)h>m8yb*}YhV*mNEn@;s$|fiv|sxmQo46# zZknev_XDe(zh)cTd2yPUVP7xRljQB3`Z7Ik7=66EHX_kiqdq68*n-T~0g6K&Eja_8 z{^3m-wgLm*3wF=du4a)HUyh&iqXDm$+!MRo36PPBnvEt#cU!B~0Np87(ft>XrUbba zCcB*JTTvUkpSfd%4N0dmSnHq+*F_Hnna7dTt#&q{o1z#CO`8|qW@GBhXHa;(2pkKE z3+-Mmys*p1^!q->66OhP3_ZS3ZpkF@#bnH^>)pmsXXs!`gUxvbmD#e_aG{4SPyFJm z7*nKO#$nKbknkB7@Iw@>oVz2Loi(^yRW)X(he}Te?rz2zW^WNjjY7a{PbG^E<27Wx z!&KJWUjt!!t~7(=A|>^9!Yi9J112XeT>Fnk(zkPrH$^oR=A*xWcSf8x;<*o+3_bOp z+hq2|?9T|VjEgU->K~iDz-D`BC1p>N?;2||tg$yyi1;u7#3X}|iW*lK#7$U`eWCcm z9Z8Xm_1eO(F7`LDDA!brVk{mzqW6-N$$Ne+;~YRs5`2qXxpZ3tinQ;1c3j&cyYM}X zJ3}ApqZ|CX#8IqgrquqmA3Xvihh1cTcF%a`-R;_`m2J7#A)oz`FQ=<~;5}TrnhGm{ zIn3#RUDDYu1ix9-T5Z%{&q<&x^5V1msU^Jjf@Zr2em!{{N@iDqJ&l}@V6aDOyAF#Y z3@$_yqhKXFjmPz$t?z9jpD1#q^4%-j4%xO#`aW5Kknta>h6Rcrzjiv+0O97;=Nkq+ zQ0LB=v5cdgA>uDru2b|95#W?q;CgI;eFmkQ`c|IvWs*@)Z;Xa{?4W<({i}>dPAa$O zDd9TpQzr8^^SrZlQhJNBLhV%%{Dr?|`3(vVe)dJ&C+?~C0myFxhRv%xG_SMvE$o`l z9FLOX6YnSI8nu`cDP{(K+>W^#e85GE=v16=H@o^a-(W}VHOH?cQ9dzO4=@Ub!{umw zRn&X59cgn?do$~$TSgJ!OjS*Fu0MW#_4XUWFZtPmMu#?aP5U?XIMq5zcA`Q_;knR0 zyqlDxvYY^ls>8JxIPdsD33PbwC+m{hqt3Zv3ZM1c6kZFQKDSe8azt*EvRU53)+35! z%Lb&cSmKQoutE#tMPHv2rvNdZASOdIvjd{}Sr{OU)3}HvS{*qHeNY%MsgF=nU~^;xEN`CziqNMVssU(YPK4r)=N1SuL}qn+$(o)CO!#6 z`PCT5`&&q#=02KTJc=2t2PuOvAiPpJNJNAq4bfqT%)UJ=ZWAY%9VmzxtIkT2O~cNp zZ}nX3In5lxt;J=^CVhOS0IwFPQY^`{$i-Y&RY-PKLDWMFw?D$>CL%yt_cbjRv+w~W z%_7Vv(*1p^KB^H_j`-}#21<6_ySd4kWTOMOx(us@(j=CbvMqsD0Rd*z5|su;=Qj)X zRKC&-BZrO`UY(Mj*T91T>T39r#`WbvX2*1TswT^``l@Q455Z={T4EhOd-AJqf`lW_ zCrN$Y`etyRc+x{iO(XJ8LTPFTN;evt}v=nBRq@%zP0eA}emwy)7p}ASek_9j)^ZD1&*U5H{$Va$H zPW;b_Ej)B2kb(Fsb_O+2f~3lB?wNPGdA^K&#Kb{Mb7b=fS=Z0nLdbw3dmo&ag6wpW z|3_Gtbw7$0$8@g6>J^UUO-bXp=pHjT5(@zd^LdXBIT_&&m!_^P$Y8@dc4Q*~5gK(mu ztPjC@++t6v*!fgIP}BYKR8r4#%0@LY*92joK2t2Ln_VJZeI0_pR;A6n%nYI*MpL0o z`|b+_N*H8$(hyUQic?p|PsSSVTzR@wpuUv>xp!3r@IXa6F&6k0S<}o*(KB$uehdw6 zS{%=w{PwQGs}q(T8go*#Oy9R0rj1q1o3VrH5XLD=aWsr%eaBaHj_CY&-MfB>bm%Ky zrNJl1{)R7PUr`10MA_oysxmzmo|GibKCRf%OJBEN>Qfx15InTTD7E#tr5D<;;FHpo0#rP)_Ua} zbxLN-^&%{F&9oXL>I6ikq+k9aHw_b{dZO?^N{GQRP)H?a0 z6DM0fFq}ohBxqql6P{C5=+uu}Uh@WXI&EP*W*Ed~)xA5HG!^IG-SvH~@57mnGT(YX zK?E?moTfhO1|xg!Yz>m9A98_Bz56bphW<&WEtVob_EXZKGK&Kbj_ewP*vzs0v)2m! z53tgYg_p^jeji)y{y%}oUq|f6*?GTkY8(s64|UUzwwEr~h8v7`>$K32{T>swdiV+a zd*JQ8LVbO%tPX3p_UDw!8rM>G=kIf1AL7S(m_2W=`bU)1nqT}*-TMfyVJk@2;q_O` zuM2SvmpIO-Fh)4LBO4jLv!7T=%5W5BuEQCsOAou7->~;Cc{z6lW~HA(>mHqJc(HHD z&ksP@gUC~>hCgF{BvkU;5x+}DL6@Z@0oe=OyTc<54hbRM?KJgwdPMi^p2?Cw+s1SY zK65Qz#!)ws-Iz$kio4yxuip1!+b2s3@WKYJG%RhPj+hjSr$6}U+MTTI(>~&kq%(Zh zV;PD{ceQ@93#5ix%T8!-L-YJNmv1mwFJ&8_P3yLRZ$D=H9s{9b>5;E%PQwTvpw*>O zg1ol+_S>9n14^iueO36G+xaZ25{4#iPT*Nce-=K$Zg8*8#rMvr1w43U-@km19f-{#f_PtxAl z2Gk1H$x>M@e{or#L?tz9e|}j7l+U_6cd$1LO1}^|g!}C**m&$*LwcHjl&XFfz;L=u zT0?=*gy4*YGsKoA#oT-N%t2j@+u_jN+j6xk=(i@Ss9;PZ?xx>33=LUeNc51e%fM11 z>nv*iA`3EH*Q}ZbzW!AZ4}xzt0xvJ5K?PU$1F%Qd#G9K$JkLrCZ)PNcA38>(76G*r z{Y#L_@pJ?WnReC3C?Q*VfeSvH=fYqwkfyH+NI*B#qN>HPZ~@(1{zKY ziyGfvNoWdaIYe>IodL0N>{n?`C|@d?bj;`CVU%R{W^TxNg3B3;e}amJORxGDsX{=^ zxQ46lg!Y}g>NyJwdBH$LZEIL|J`~U$EupCqQOWw|oflol7gnM4TptC@SzVk+$m{&2_(^&shKmt652xN7yJs*3#Bx(gr#9O+VVFP*h{oIBxH z9q;YQhfqEew#i?DFjPHC{G-|OVMS9Z zHKUcp+J3Pe{XpaSoRztU3XGy$rhI?OP0HisT8M_J;%Ih*UoA6RlqbWK4}qO;0@Rmh zQfM~9R~7BQQ?9Ul789e9r?5k7=*U-8aqab?XPrVmNKC#UdDf~RYBA|3@!rETV~-zu zzi7v^Iipbn9An>0Bd(yHEB~HJHh0N{zYwM^T{`K%_3u7^Zfv#VKZhw~ zS$K2V2R~xyg?Gp6iMVR6z&}VkvX#y7DfI;4wzWsAo=V)PuJCdZ{3@2xYFJzRuGPAn zcHH~oiYfaUOxz>qOUIa(;=E_5A*JnaJ(eZdTxgf=O{6Z~mNhsSEW*ua6Br(1fqosb zA`{0$%oi^5@=AV$0h69$4UPe9X0bWf#${7)Ui3T?{#JjV^0zmta@g7>!&KeK)lwfIBL8-g?IXVkZ% z6~%LdN}3WPPs#6_-*+F|=>3Z~rjp_y_){vaZF1aiYR{tG-Af6AgL@rk*hc#4S6G+P z(a$ClTxg>7vO9Wnqc=ij-`Phh72iSNs*h{V`SoZg>PnFU3768g8tP718B6DaD_WUW^U?lfz>R|nrugW} zk<-;+8PszN*SYmd7LxC%&8+$DwL7-#EG@)%p(2OfLqT)i42j^;^Z846ZUp=~(=Rj35 zh+}vD!uM*P^dsPV{@MatMr#wX+-;)ic{li+8G=*MHm`e@e_?k$B+N-Ng}f*$cpg+< zgj9G=GXJChbHn6T5HW@mUA)dh*7Xa=tpzf;Rm%D+uf+(U=6sYc|F@ZGWfBHVk1g_3 zVu@F<`TSm1c6MXJ^s*l5O)%$Bq_WphluR#^^huCT`b&f6WeLZ^Sv`yL%^bVd85f)I zVuQR$)3@Hy~@1H_M1brOJ zW$tkP)#I8ubdSm5OX8I%;?mLpt^PY<`20PM^B4Ay0C6OpKgGLGRV?Sxm!g={&+_kX z$-AjPM^@2DeQVV+zYt(#ggLS(Ej@DqcMDXrsR6Z`ViEOUEnK#f3`r*&_5#jzr+o zpeNdvfdT_Ky26wN-G4g#7UtQN64;rcrGI-9F$(-*AKlBK;_?~J3Vm>YnFLk|Pc#kF zrX{xQD8qqgsR<}*K;&UG!_X4dH)~woeK3qatULb9@K}4SE_C%PV^HS89$`$svuwZM zIA4J^8{11CLHPbPtgvyGYVlAI2tm>f=#F1OLymwP(Tmb*C|NW z?$ah{9wm)yF2MnDJ*snn8l?e;+#8aYrVWfCO3VnSwXiy|$4~;Rc%4|^ho1fPG=y(o zKPUTY5CFfIt^Jl;`Q51Rys6Jf+SnMRJK>a^(mNTS^Yhi-jRjk~AcS-a*DhXDY?kxh zZi^E0+-o=DgC~6ctF<(4K|bwi$Imjw!Q7(na|qGzCou;ja7qXyL@*GVJ}17AAe0l@ zn>~;4|Mm+N!P4PZNPlQ(w(5Rg|DxWw&B5C0D-~3FSP``FT{jILQ z9^pJ5vuCS*ii`bb*OGQp!2N}fQtACJdx;= zL?)_SgL$7W1w$7V)>xFeL|w(ukjhy?nw>x+p2aY$5qZd0d8f{Sb*(oPhGM|BdCigx zW8ItiXIu#2yLT0DwCd8bPgF(ur16*1rhUG-;lX>!neTZN>@}S3o0nySNhh(*aK=JZ zMaPuaJj>#01D>ujQOw9Q?6Z9%rXL0m!onuox&o7wDg#Y2#Ma}+i)t*Dd&PDOX{=m)_`<5WMGWw7`~6qO^<4+>4=J42FU*cH<4@J_BcwLTiE0=! z$=CebYxlc8QhW^C@h4#q%^$-g*GwRlu8US>WHie!o*KLs8;NA(CWIb`&VTvfxYNhJ zf<^hMktExC7F=hD^WJr&XYqn7cFkQ<( zbHm4&ss_iATTp$wS>ufWPDZwZ#7I4)I3vgJ-LdXC3>#74Tn z@rk72A$qeiD`}@l2=FfP5P6@Kgfs;?v_M~cOU-oX%t2Ad<|9sk$&ynFZN9~Xhsp80 zfBI|%Vx{p@=b?OVxe?Rk8wOW+{r9BxmIn4^q#uNJD8?e+?}R9chzxljm%sg% zcbzK>a$mM!2zd6PPvl6X=)<<0U6=W&^%K1_32V4AZI`N&-GaLnnAFDwVf8GLb4#Wfl@d^((ls_is|#)D4swst*f zJ#OiBbiO|o>TI&Ou5&PaPiIeBp~Wwj(dczpQE7s+e#2sm-KhrWa{(&&u<^6w~G@Kj@s2dy>DE9@$!UX~j`inz~(GUSmPfKPgC?8hmCfz2r zIA4)U8sU>|j+&JU77!sW3irr|nqO#l)TDP>?W{pVLEIFX9nb|gz zT%jNqGunpFS{kKX-aRP&MJ<*xJH;K5V7IG__gTK- z5NARR#=mp?$}c8v z29{g7V|DPeY3Ao;s&TA4X|6%m8SUIilIf&ePBlN5);VY5q~R(O;kY{4WRv46o#Ew8 zkwOes8y!#tyqnoqxX!oVOG5@A?{lM_ILsVJ*%D zXXLws9h2X219`Ywy$|B`{^}L7hvqW*HD|x2sv7k$5@qp&T8 zP4X1pu$p44xUxFY{2KfHN{_~4iz=dCl@=>u?oHNh1sK9i6amoYbG{&c& z^_oE>uX&hC>3JiydSiW|4Q~Pu;l@Q$#rhFcAEamdfck>2atFLaznukKcTc(;Rp_%z6@46PQw)odYWvD@bJJXQaj9LGha88uRfu2s2fFCZq z{Ah7;9GA;n^NV?iY~8E^;sRvwGT-BX5}-=;F47Q3%JZzccd1fOtiQ8b7*~^&LBiT*pU@sx>PGZ0lhL!ZzpWwMEu$O66JHJwHI%3Oh^iW_O9Ebe`le=bFF&T*x>%?+C;!r!~r z!6nowq{4cdD8S-(mv^eM-Z@9!X~OU95Xzs-wbjpKq;Bcal@hL>{TmJrp+RQA4G)oI z^A|^W*~|f3?^l{KqvaUryO+fCKz0P-squEkKrMPLw8O8|FNshLXlm)Qmay$?0*YhF61qO3PyCYzdFc!V2HZe3E4T1o-y#b)$)IB&P8V^2%p2xf2E zwziCJdF~q~(}1ON-%c&yPWRr0+Qy;g6O`i=Oa7h$6r0sIJ7ovNk|Xa{ z*Kv?$1)yVZ#mC;_odwnGhgy_+FlBfTIor?B%q`G*vxK-Ig>x4EoS#C@c^`K~bMV|fBQ zphbw%(Yy8bn~K4NfL_mt6{wo7y?*W*Pk|T2DAH+!UhpS(icqv@a;AoHm}6LtO4_!@ zs7_(EP`eRD5d+Q+qN=O+QlyB1&$ z7(X~ehdf7fK0}XQ#lp*03;d-hm2_17z5ew0wNz5I{ueX47c;`z{NJ_nBYSU|kJ(jt zx71BiK9nq1nl#q#Oedg75Zic3eM;cb4yQ%Iv!ipAa%^YubjzFin;rTwP% zfnW?#q>ej7plsp@Np_meqZmPS=*Y=RjFj4UI9X%2v6EgGlC+vpHl`;@0KEqOoTGW? zt>g(X;oz(c|2{|iDgO!2|I?m|Y~UQN+|4+>9UVR-f}Iy%kcA#yv9T{0kss@uZ?d3Ciq! zMG<$T{=(l{-ud)vK+e_l%W)O_!50)`O7C)~`GyAR?CIgDN{93+a@+KcVG3z35>KUi zkgb-w^4==>`QB2z+0qqsiYdMpRW<5U_r_i@M!Mj`&jvha5md=zPop0w@B1dAztK0;)5NbM;pgN zQMO46CkOdc%XU}oqh3DGjp5qz`MUx|Rh49&jzmjel8u3mMAgxEY!*aT>pT)r^AGnr zg5xsPMS*i$0%2y48AX9|Ysq&UBGN6aC*k&AW&*hBLB)bBXANNiIt!i|NHX)=D$Ef6bRlNi5f!ue?*KvVPhPL5=$IsWIi!UWHT2b(qk?tO?8Q;LJIH_WbdCp3m?3eZD`DUYSA8DM_p`8l7ta`FDFx z12CrOd@~6c*(qBc18(p1`RE^h8A9fa+9s>$;e9=>Yn69VRp6)pW0grJGz&f(3J24t zuJlsu3d_Pok;NRr)d6GwfxxxB^{#f4Mw&}L*=8WVGa$IrG* z*yzD^Uypk3<}<{qm+Z}}N2bCU)1VlUL#A+4}&vs%*VCgRvCgiik6)_x2+H zSC>ky_dwL$5Xi=@wG4xlN?*&!6TnB^y?;MBJZfk}TvB%$Z*4O9)vkK%;TcL@n^Wz? zKy${u`N0qEX@}ASt();@dhNlhz$#vdL2P?j-hAT=*dB6V?+2)@Z>r8+uyv7E>fXfwRbLWhdzM*}JfxtCm9+`7@+~-z&$9B( zR9=YpO~dL+>wdtAu4u(@KcEBHdb`89e@5Cmv+Mc;aDUFy zUsk%L7aK*vNO_&Fx2>h&iP*?08EC=%i1HGNlM#7@MQgDBD^COFMN|P*kx3Rc`{zG? z`Nl#;?%~+x(L*tQ(t!UEvrt2JekFCC;5BM9gA^WlWKabX}2%EiUmqx zlAl;7plRUh$!9#sSgutk7AJ)ITaxnCYh#4go;(d{>rwuctg_}s9AkjGy^pQJca2fp>!py zEzGRP{hX&<*#$O!yKz$uAy-0iCioDcC6qr2Xh0&ezp0!-%XaZmtvM11G_+)_YJ9#nB$9$vc2~SJ*q;xFeI90~%qJwI z<*s^VKt-!2>%GC?T>`ra^NVDQe9QJHuQf#Ij(^0~_WD48*Uo$7n@8C;DctCd*u&kD zBTJ%HfV%FUxLPTg2D5Q#7N8`gA4$m<=4@&#&tsyUvEf!g(N*>#-Fi=1>I4L;T^CRr}FHXFHON?g22Oe!klzF&ml&-hvjzIS3 z@K*3Jw+qy|s=ZJ%*ASWbAhsKh$h1;M=*+SO6oRHWolhQ)x8DY8Fk_b4!$h?QPzw z#u>J{b0tcV7r7ENWG-J^H5`#BP2bVW#0$yUB z<)cxd&8`Qp?;E!q6yoqdR-D)4ha4wXcMZq=bLe#KL}Nt0aAlzdv&nBgzD#uB$HYb+ z=C|`iy^B7WT&ZhVl`XN&N=^`-ao$Y_xZk`=Xsy!d)aQmreHG+bRysie>C}T zhfJz&!r_ZuRWOi${#qx!Hi*+jbOxD7VV~1#3vVE3+{nywC(Sdbs2pR#<5W+$3+#Cge6dgwC-LoO9b^1qBVwnFv;5TkBMEO-SO$fchIGb+-i`sj5Il+*PNG|eRS zdixiRB`W;kt`)(QS{8fyptVB{BGUimDj;?_Qe)-w!Ad<=xL`Ffq-Pl5?!_G)sVpcD z7ij`buGXcpRcgDZtkfz!F30QSK2XfuLH4I8kMy8U34M3pvcWoIpel?NE@A}TqSsVG zld>TlnA+|r61Y$W(6c+rkQl&x{U%G#tU(F?$jDcjZ|{@1sJnOlLtKwp1}cW}DbCSl zp?~?m(|uB?KE^MKH(|k>Z4pLDPwUO;A)`Nkj6LP_gqpHfPJO6Q^UCV-34LUd?&hgs z^rGXcdn+Q7GnDke{$Rt#6s?W3oac2AZ-wD1Gm;Gx2j{Aq-V6=cawHSJG6RJ#2Spn& z+si?oW)n4-Rsr(hhTE^CrQcGqrpsm5eT{YloHS8@&@H0v!F(%Ko&ZM%S3s=W?Otmv? z%5J>e${n28OzjbPYQ4Fp+o?!IM8?IsWuZy+#Rt7=_mR5;4{eMOY;m>hnMcornR(i* zCX(il8sjuje0Q@#8|fVa69pq~<0oE8w=}zlPxU=H!Y}xDcehaNrd(^3-G^HyU zaVxqP2lvga26go7kaNLT8i72Fnedgpcar794fjj(2#us4T2hQ>eB4+KhvGfuya7w&?y&1+h(S`phL&W>-P$usb8yzo1hJJV)C;5hP`9KY? zl>Zp+)Beqz5&c|KQDM@uJ*BgmyN8n`?y!9#Ykrs%_jUd^&XO?qrnutBN72drrzTe5 zR_LSX^ITnb$#kAXB%iCRa0rTtW3=~026$(q>r0vwd1KfFwW>RhPwyu#q9 z$(Lc~n&O>dP{F0!s|n`kvAyHP4c-`@UR& z*@iMGNPmCZ@_=GO#!h0c0J4cHlUZaR z`PVF3!Y~!k2}4QZF7C(0X7JK9&ydghxOEUs&x1g=2QqpUa{64ZR&7Bn6T+QY+vX*} zelQjo>6e8p=XfajvixAEw>#*f(vD5`s@D&GVc%6ee0_@nWd9lxpu`7=Zo;ADv};q! znu!5u*`$m{56TxG!`%Oj%bkWNNiX#Ud$Ad4n^f} zm7zz%7uHu7(`rVLI{qC~EnGV{)2a*0Crd@h*)(!29&t=87t`S8y(3FT*R>DuDPd)y z-&6_b&Q#tnQn%CgAh>`Q;6Gz=!UFi$fLLdO<9C^hv(x`GoRTb-8@d=1u;1~UTnJ!s SfYCMp#1nvZ>%4bu2mS^>?cc}% diff --git a/accelerator/report/report.tex b/accelerator/report/report.tex index e057cbe..d33742f 100644 --- a/accelerator/report/report.tex +++ b/accelerator/report/report.tex @@ -33,7 +33,7 @@ \maketitle \begin{abstract} -A 2 MeV proton beam is used to annihilate lithium and flourine atoms from a LiF foil. Alpha particle detection is used to verify the nuclear reactions involved. The probabilty cross-sections are computed for the nuclear reactions \FpaO\ and \LipaHe. +A 2 MeV proton beam is used to annihilate lithium and flourine atoms from a LiF foil, transforming them to helium and oxygen through the nuclear reactions \FpaO\ and \LipaHe. Alpha particle detection is used to verify the nuclear reactions involved. Rutherford scattering of a proton against silicon and copper atoms is used to calibrate a kinetic energy scale to the channels of the MCA, but some problems with this are indicated by an imprecise fitting to the nuclear reaction peaks. The nuclear reaction peaks are nonetheless identified and their differential cross-section ratios are computed and compared positively with known values. Absolute cross-sections are also computed, but are dubious due to an unknown foil thickness. \end{abstract} %───────────── @@ -41,7 +41,7 @@ A 2 MeV proton beam is used to annihilate lithium and flourine atoms from a LiF \label{sec:intro} The nuclear properties of lithium are of interest, especially the reaction differential cross-section for \LipaHe. Studies of the Sun's photosphere show the abundance of lithium relative to hydrogen and helium less nearly an order of magnitude. A large cross-section for \LipaHe\ is believed responsible; it predicts that ionized hydrogen will readily collide with lithium, transforming the lithium to helium, and emitting an alpha particle. This is one transition that takes place in stellar nucleosynthesis, and is maintained as fresh lithium is carried toward a star's core by convective currents, but even with this process and reaction in mind, modern astrophysicists have yet to completely explain the lack of abundant lithium in our sun's photosphere. Studies of the nuclear properties of lithium could elucidate a better stellar structure model, but this is outside the scope of our study. \cite{Carroll&Ostlie} - Fluorine can be rearranged in a similar fashion. It is one of the rarest elements observed by astronomers, and thought to be for the same reasons: it is readily rearranged by a proton to produce oxygen and an alpha particle. Figure~\ref{fig:abundances} shows relative abundances of many elements. To judge whether the high-probability explanation is plausible, we will determine the reaction differential cross-sections of the fluorine-proton reaction \FpaO\ and the lithium-proton reaction \LipaHe. In this study, we performed a prompt radiation analysis by observing the alpha particle products of each reaction, and from the kinetic energy spectrum of these products, the differential cross-sections can be computed.\\ + Fluorine can be rearranged in a similar fashion. It is one of the rarest elements observed by astronomers, and thought to be for the same reasons: it is readily rearranged by a proton to produce oxygen and an alpha particle. Figure~\ref{fig:abundances} shows relative abundances of many elements. To judge whether the high-probability explanation is plausible, we will determine the reaction differential cross-sections of the fluorine-proton reaction \FpaO\ and the lithium-proton reaction \LipaHe. In this study, we performed a prompt radiation analysis by observing the alpha particle products of each reaction, and from the kinetic energy spectrum of these products, the differential cross-sections are computed.\\ \begin{figure} \center @@ -57,7 +57,7 @@ A 2 MeV proton beam is used to annihilate lithium and flourine atoms from a LiF \section{Proton Beam and Detector} \label{sec:detector} - The Tandem Van de Graff Accelerator Lab provided a $1.95\pm0.05$ MeV proton beam for three experiments. When incident on a lithium-fluoride foil, we expect the nuclear reactions described in Section~\ref{sec:reactions} to occur. When the beam is incident on a silicon or copper foil, we expect Rutherford scattering. We used a circular normal-faced surface barrier detector to observe the alpha particle products of \FpaO\ and \LipaHe\ and the protons from Rutherford scattering. \cite{ADVLABACCEL} The detector was positioned at $149.95\degree\pm0.05\degree$ from the proton beam, which we define as the lab frame of reference; see figure~\ref{fig:detector}. The detector covers a solid angle $\Omega = 1.895\pm0.001$ steradians. In each experiment, the beam is run until total integrated charged reaches $2\times10^{-4}$ coulombs. + The Tandem Van de Graff Accelerator Lab provided a $1.95\pm0.05$ MeV proton beam for three experiments. When incident on a lithium-fluoride foil, we expect the nuclear reactions described in Section~\ref{sec:reactions} to occur. When the beam is incident on a silicon or copper foil, we expect Rutherford scattering. We used a circular normal-faced surface barrier detector to observe the alpha particle products of \FpaO\ and \LipaHe\ and the protons from Rutherford scattering. \cite{ADVLABACCEL} The detector was positioned at $149.95\degree\pm0.05\degree$ from the proton beam, which we define as the lab frame of reference; see figure~\ref{fig:detector}. The detector covers a solid angle $\Omega = 1.895\pm0.001$ steradians. In each experiment, the beam is run until total integrated charged reaches $2.00000\times10^{-4}$ coulombs, corresponding to a total number of incident beam particles n = 1.24830$\times$10$^{15}$ protons. \begin{figure} \center @@ -138,8 +138,9 @@ A 2 MeV proton beam is used to annihilate lithium and flourine atoms from a LiF \begin{table} \center - \begin{tabular}{rl} - Reaction&Energy (MeV) \\ + \begin{tabular}{rlll} + Reaction&Energy (MeV)&\\ + &Computed&Pub.&Dev. \\ \hline % From Kayani's(?) kinematic factors %$^{28}$Si(p)&0.352$\pm$0.009\\ @@ -148,25 +149,25 @@ A 2 MeV proton beam is used to annihilate lithium and flourine atoms from a LiF $^{28}$Si(p)&1.706$\pm$0.05\\ $^{63}$Cu(p)&1.837$\pm$0.05\\ $^{65}$Cu(p)&1.840$\pm$0.05\\ - \FpaO&7.949\pm0.03\\ - \LipaHe&7.688\pm0.03\\ + \FpaO&7.949$\pm$0.03&6.9&1.049(15\%)\\ + \LipaHe&7.688$\pm$0.03&7.7&0.012(0.2\%)\\ \end{tabular} - \caption{Expected peaks for nuclear reactions and maximum energy of a Rutherford scattered proton. Rutherford scattering energies are computed using equation~\ref{eq:rutherford} with the kinematic factors from table~\ref{tab:nuclei}.} + \caption{Expected kinetic energies for nuclear reactions and maximum energy of a Rutherford scattered proton. Rutherford scattering energies are computed using equation~\ref{eq:rutherford} with the kinematic factors from table~\ref{tab:nuclei}. The expected alpha particle energy for \LipaHe\ is in good agreement with that published in Chu et al., but the expectation for \FpaO\ deviates by 15\%. \cite{Chuetal}} \label{tab:predictions} \end{table} \subsection{Nuclear Reaction Peaks} \label{sec:reactions:peaks} - The kinetic energy of the alpha particle products from \FpaO and \LipaHe are characteristic of those reactions. Binding energies of the involved nuclei are significant compared to their rest masses, so the relativistic equivalence must be considered when computing the energy Q released or absorbed during the reaction, \cite{ADVLABACCEL} so + The kinetic energy of the alpha particle products from \FpaO\ and \LipaHe\ are characteristic of those reactions. Binding energies of the involved nuclei are significant compared to their rest masses, so the relativistic equivalence must be considered when computing the energy Q released or absorbed during the reaction, \cite{ADVLABACCEL} so \begin{equation} Q = \delta(mc^2) = (M1+M2)c^2 - (M3+M4)c^2. \label{eq:Q} \end{equation} - Given this method of computing Q, non-relativistic conservation of total energy and momentum provide a sufficiently useful expression for the kinetic energy of the alpha particle product E3 as a function of the kinematic quantities described in figure~\ref{fig:kinematics}. \cite{Ziegler_1975} With results in table~\ref{tab:predictions}, the expectation values of E3 were computed as + After computing Q this way, non-relativistic conservations of total energy and momentum provide a sufficiently useful expression for the kinetic energy of the alpha particle product E3 as a function of the kinematic quantities described in figure~\ref{fig:kinematics}. \cite{Ziegler_1975} With results in table~\ref{tab:predictions}, the expectation values of E3 were computed as \cite{Meisel1996} \begin{equation} E3^{1/2} = A\pm(A^2+B)^{1/2}, @@ -175,7 +176,7 @@ A 2 MeV proton beam is used to annihilate lithium and flourine atoms from a LiF A = [(M1{\times}M3{\times}E1)^{1/2}/(M3+M4)]cos\theta \end{equation} and \begin{equation} - B = [M4{\times}Q+E1(M4-M1)]/(M3+M4). \cite{Meisel1996} + B = [M4{\times}Q+E1(M4-M1)]/(M3+M4). \end{equation} %\begin{equation} @@ -191,14 +192,12 @@ A 2 MeV proton beam is used to annihilate lithium and flourine atoms from a LiF \section{Results} \label{sec:results} - The Rutherford scattering curves used to calibrate the energy scale are shown in figures \ref{fig:rutherford} and /ref{fig:calibration}. The best fit for this linear scale is + The Rutherford scattering curves used to calibrate the energy scale are shown in figures \ref{fig:rutherford} and \ref{fig:calibration}. The best fit for this linear scale is \begin{equation} x = 0.009\ \textrm{MeV}\times\textrm{channel}+0.00392\ \textrm{MeV}. \label{eq:scale} \end{equation} - An attempt is made to locate the predicted peaks, and instead two peaks are found at . - \begin{figure} \center \includegraphics[width=6in]{rutherford.pdf} @@ -213,21 +212,73 @@ A 2 MeV proton beam is used to annihilate lithium and flourine atoms from a LiF \label{fig:calibration} \end{figure} + An attempt is made to locate the predicted peaks for the alpha particle products of \LipaHe\ and \FpaO\ in figure~\ref{fig:peaks}. Two peaks are visible, but not at the predicted energies under our calibration. Presuming our calibration is off, the counts under these peaks are integrated to be used in computing the differential cross-section. + \begin{figure} \center \includegraphics[width=6in]{nuclear.pdf} - \caption{Two peaks were observed not far from the predicted energies along the calibrated scale. The peaks do not match the expected values, but we progress under the assumption that the scale is wrong and these are in fact the alpha particles predicted as products of the \LipaHe and \FpaO reactions.} + \caption{Two peaks were observed not far from the predicted energies along the calibrated scale. The peaks do not match the computed expected values (blue), but we progress under the assumption that the scale is wrong and the peaks identified (purple) are in fact the alpha particles predicted as products of the \LipaHe\ and \FpaO\ reactions.} \label{fig:peaks} \end{figure} + \subsection{Yield and Differential Cross-Sections} + \label{sec:results:cross-sections} + + We observed the reaction yield Y($\theta$) as the total count integrated under the curve for each reaction. Computing the differential cross-section for a reaction is done using \cite{Ziegler_1975} + \begin{equation} + \sigma(\theta,E) = Y(\theta,E)/(n N \delta x \Omega), + \label{eq:crosssection} + \end{equation} + where \\ + \begin{tabular}[h!]{llll} + ~~~~~~~~~~~~~~~~~~~~~&$\theta$&=&detector angle,\\ + ~~~~~~~~~~~~~~~~~~~~~&E&=&energy referring to a specific reaction,\\ + ~~~~~~~~~~~~~~~~~~~~~&n&=&number of incident beam particles,\\ + ~~~~~~~~~~~~~~~~~~~~~&N&=&volume number density of target atoms,\\ + ~~~~~~~~~~~~~~~~~~~~~&$\delta$x&=&target thickness,\\ + ~~~~~~~~~~~~~~~~~~~~~&$\Omega$&=&detector solid angle.\\ + \end{tabular} + + + Beam and detector parameters can be referenced from section~\ref{sec:detector}. The number densities for Li and F in an LiF foil are N(Li) = N(F) = 6.16 $\times$10$^{22}$ atoms/cm$^3$. \cite{ashcroft1976solid} However, the foil thickness is not known. Our best guess is 200$\pm$20 \AA, but not knowing this thickness may dictate we only compute the ratio of the cross-sections rather than their absolute values. From equation \ref{eq:crosssection}, + \begin{equation} + \frac{\sigma(E_1)}{\sigma(E_2)} = \frac{Y(E_1)}{Y(E_2)}. + \label{eq:ratios} + \end{equation} + + \begin{table} + \center + \begin{tabular}{rllrr} + Reaction&Yield&Ratio&$\sigma$ (cm$^2$)&$\sigma$ (barns)\\ + \hline + \LipaHe&428&3.12&(1.48$\pm$0.15)$\times$10$^{-38}$&(1.48$\pm$0.15)$\times$10$^-14$ \\ + \FpaO&137&0.320&(0.475$\pm$0.047)$\times$10$^{-38}$&(0.475$\pm$0.047)$\times$10$^-14$\\ + + \end{tabular} + \caption{Yields and differential cross-sections computed for the reactions of interest. The ratio of each reaction's cross-section to the other is tabulated. The cross-section is computed using a foil thickness of 200$\pm$20\AA, but this is only a guess.} + \label{tab:results} + \end{table} + + \subsection{Foil Thickness} + \label{sec:results:thickness} + The foil thickness $\delta x$ is not known a priori but it may be determined using the Rutherford scattering of the alpha particle product from \FpaO. The film thickness is less than 1000 \AA, so it can be computed using the surface energy approximation as \cite{Feldman1986} + \begin{equation} + \delta x = \delta E1 \{K [dE/dx]_{E1} + |cos\theta|^{-1} [dE/dx]_{KE1}\}^{-1}. + \end{equation} + + Here, $\delta$ E1 is the full width at half maximum of the Rutherford-scattered peak given a film thickness $\delta$ x. The stopping power terms [dE/dx] are referenced from appendix F of Chu et al., \cite{Chuetal} where E1 refers to the incident alpha particle (the product of the \FpaO\ reaction), and KE1 refers to the Rutherford-scattered alpha particle. The table of stopping cross sections consulted \cite{Chuetal} did not provide values for the known alpha particle energy from this reaction, so an attempt was made to compute the value using the polynomial fit provided in table VII of that text. The output from this attempt was nonsense, so while this is a viable method for computing the foil thickness, thus allowing a better computation of the differential cross-sections of both reactions, the task of computing the foil thickness is not completed here. + + %───────────── \section{Conclusion} \label{sec:conclusion} - The energy scale may be off due to a poor choice in locating the Rutherford scattering energy along the curve, however, look at the different between selecting the 100\% mark over the 20\% reveals a change in the scale of only about .018 MeV, less than 10\% of the scale's deviation from the observed nuclear reaction peaks. This indicates a more significant flaw in the analysis remains to be discovered. + The energy scale may be mis-calibrated due to a poor choice in locating the Rutherford scattering energy along the curve, however, looking at the difference between selecting the 100\% mark rather than the 20\% reveals a scale shift of only ~.018 MeV, less than 10\% of the scale's deviation from the observed nuclear reaction peaks. This indicates a more significant flaw in the analysis remains to be uncovered. However, the peaks in the kinetic energy spectrum related to \FpaO\ and \LipaHe\ were identified near 7.175 MeV and 7.583 MeV respectively, and counts were taken by integrating under those features. + The ratio of the counts was used to compute the ratio of the differential cross-sections between these reactions, as in equation~\ref{eq:ratios}. For \FpaO\ relative to \LipaHe\, the ratio is 0.320. Chu et al. report yields for \FpaO\ and \LipaHe\ at 3 counts/$\mu$C and 9 counts/$\mu$C respectively \cite{Chuetal}, giving a ratio of 1/3. Our result deviates by 4\%, indicating a strong agreement with published values. + + Absolute values for the differential cross-sections of each reaction are computed in cm$^{2}$ and in barns, but these rely on a guess for the LiF foil thickness. An attempt was made to compute the foil thickness, but it ended in no useful values. With further work, this computation may be possible using only the information contained in this study. Without a value, the absolute differential cross-sections have little meaning, and can't inform the availability of these processes for stellar nuclear processes. However, it is clear that there is at least some significant availability, which could be sufficient in the proton-rich environment of stellar plasmas. \printbibliography - \end{document} diff --git a/adv_lab.bib b/adv_lab.bib index 06be0e0..2993ef7 100644 --- a/adv_lab.bib +++ b/adv_lab.bib @@ -21,6 +21,28 @@ + +@book{Feldman1986, + title={Fundamentals of Surface Thin Film Analysis}, + author={Feldman, Leonard C.}, + isbn={9780135005705}, + publisher={Prentice Hall PTR}, + year={1986}, + month=aug +} + + +@book{ashcroft1976solid, + title={Solid state physics}, + author={Ashcroft, N.W. and Mermin, N.D.}, + isbn={9780030493461}, + lccn={lc74009772}, + series={Science: Physics}, + url={https://books.google.com/books?id=FRZRAAAAMAAJ}, + year={1976}, + publisher={Saunders College} +} + @Inbook{Meisel1996, author="Meisel, W.", editor="Long, Gary J. From 9b47887679fb708e17b4e08695ba16ae6b9a04dd Mon Sep 17 00:00:00 2001 From: caes Date: Mon, 27 Mar 2017 19:52:31 -0400 Subject: [PATCH 2/4] finished lab 4 --- accelerator/report/thickness.pdf | Bin 0 -> 6124 bytes 1 file changed, 0 insertions(+), 0 deletions(-) create mode 100644 accelerator/report/thickness.pdf diff --git a/accelerator/report/thickness.pdf b/accelerator/report/thickness.pdf new file mode 100644 index 0000000000000000000000000000000000000000..5ba4cbeb8b653f25979ba59ec30c503867ad39a9 GIT binary patch literal 6124 zcmZ{oc_7ql+sBD4qp~LFd(Qjb|J~pF+V1PG&vgmh(9{-#h)aP4?oQpEzT10uzy&EM zRlnD+lKtM?M46n%xy{LaBD7?<{t&h#vqTH?xe+xrjefSglNmLoxj z?Z?dZ1+4F_99;@*$Z#W~TP(fpE@tNEe@o}WA`@pGPs>kGvFt!W;|NQWM(`L*JNyY6Y$O-1MBU($B|x8WMWyHA zHJJ5Ng&t7XrK-brDy6cBO_yCnNgzetE0J=M<6Mf}n_H9kRApW>mP;QClw>J=+w!hS zQwqooXN*#B75A&&$Fh!`HuF}CX1ft_J@Le?2cQL-1)sr;H+MC2UfxZ6?wPbC&hv5k z1^rXOgkknCH*%6&V|3Gd`)6$Yg9rt`*VZHzS)ThfcEKkOAw^Al_IGROCVET2377 zD{!$V0BMyHAvfCEtkLih&r&F3_sTCC{GmsV2y>zP$-CFwz9`kwjyP_||JW`nQ_u`= z`SOmj^FpTLKnuY~>!=+ceftXm#x~nYEU}iGP1x&hA&f~VbAcDCOg!rH7;btGf6e|J zzV;0LLfOWC1Rp&*z%ebb%xWL|Xy~ySt-ko zxN-40ixuxM%%j;WAeIYkt=x^bLZ8anpYb{xy!jZ+-BVJ1X_I#(X#FyAu|>Zsvv%!% z_gkXA-3IFid+s@^MHRVe0iwQVvddcghGJ92*%OEUuHR11nr6j4nt#w}bYb}RYI5Wu zyP^yW(PckRJAZF#o^BkUa&)e44YcFGiVfi0BfOPcPKXE?w71Mm%MxPTSE8OCIXc~l zPg(>Qz5N{I9BjZqIPmpIt&3mNoE`}};}ysF78cz8lQ%$Vex8^4z_qTZD(4GcEy1t; z94FXQg?jDLhk%#Lp+tSv`t$48S{7{*_&C80@JjNoS{)DH5Y_5_B6pm9#wjOu`n8LfqO!np(@!^U?!gA7fNWNg{F-q8TIY2Ol9 zs^W9rEshKyx%jYM{ab^Ru=SaPd%TefzH*G}0!CoaWb|8n{a2?-kD;byVq2AC&}l`l zXia*XB|ifT`ZQhI*J=8L1lBr;lugX$8EPDr-X%#7AbLi9f#OF#gOU%MYWcLppv9Ap z4=vQ@X@>LE=C-#k(-&8T51Q2JJ<&Dh5#USXqov!wT1#2p!fM7>p&Upf8RikxsE4iu16i5rR3 z-7TR-gp%cN7Nt|v@Q|deB<&}F%-O`uWR@Yh{$ecYBItUG&w#1FZ|JG42rxF|mDj4- zcTSzZsPBd&f6?j@YcqJ)a0Wdj-ZV8!-^Bf$vdJ#NHW@OhN7O zp|a{})8Ss`XvtNP=r@j-=>e_jPT!F0XPBH-YABmbt}UDj3cX*l?IbO)T7+F^4KKQ&|tzp08lUv$S1I=!mBQ5b?` zDyMzT64iwEeq<>f=IFxQ^0qO5TS`SM`T_3;arX6rTb&8VZyJt&ro)2wk2Vs{RnJ9I^Jr)rfseN`tB)`6 z1e8;?Q;NH^jZI_jZ&u0oaTZ+W(fMAsHxSzaeIzs$Wm4)oH)ljJDJ_bZt^Drsz}lua z?YoCy^=ocsF9A}ago?~GiitkwRfZjr0^3+t6G;ZmOXP== z(M2>SVW*z-fM4`KBzZv6N{<9#2e6(V_D5Wu7QxBZHsa^4$uGXUD6)&`1j&Uk*YIhA zl=S%P0_mdPSVaz*utFo+z}3E)!hZZ3D~>+CES+2WSuH|%D9n`UBcq{O%`D0`JpiW% zJ!%vY@#!k{Pfjp2Zc&lb+PeCASmk@J>BBYiRCMMMqcvBp0g0pjM3O&=iVRGD0L^cr z`U9lE`fg5Gz#1R}lnG!<0xh&R378-N*Z}F|23JSp0c#1;o0KG!L_jjOe`78(GyP^H z651h5kvPB_OoB^N1|%M53j%ACcneTe1ZzV92sy|vVgqa8aXO|rGFvJtlEaZwx}!0s zo^Zr31&lS(Q3s@doBx;b$Pzi9HTar_hB_RJbOL0^R<^&Xk$n68M>R60p~y+Vn#enD z2qZZu*pNhm;NLtaEidzb@*KIQzj*HHLmxBKk+YG#^|JLh?lZ+a=)3J|azgiXy>62D zWe7c!-n&E!ei=b0y?PhWJsP3Hx%J%<8%~aPIVIQ=foh}n z$5m1WJra67{8ph*i)8Hj73+ZQR3xMP1~goevlH7hgHLl35AgdV3A#2+ds*XJq3pPC-7Wa$0eJ!%#T)9$2yYdT^HWU;Qj6 zFDKz4=FuQL7;Jlxet+lAL%Z1}EO#FV z(Xf+6Ir|l8erW5N^4vQ$n;E&wOMB#f;VkpPrG_Z*11-(Yg6qKbhn3cCtWhUZ)gx4Q zGt+h%-yh|T)I5y-vdi4}(DFn+-+9%b3lzeanE`r#C(*u@K!z3_s=aQ)$JED9_j^N+ z0%5POAGMvm6daGx(S*KcG4A=9c5w5A6hNV+O0}TMHBZI*j=Jg{^J@ymNvacS!co*- z4v@)EwnR!zfP*chGV9!f&{7AXA{y^3W~ne4&vUrDh)F3qTJgJ`4%B0mCK6{dLPSS%s03n&z-g&UCc5m;x3L@$abE$o)^ZQ zucbJxQJrP%%0guKqrE6WU(k5-1+)9f3I`)Oy7+L8#`aaH9Itq2^}7^LssN_taHIE+ zCYdckDt8N+LztKqatIR4XPdOulkk`5nfXc*S@KS}seVbkkjGmSqnlLNX3a*Ur3JZh z{Q*M|*%Avccxm|`g<uKa2m#tMP$h7pHKZRdV(=iiT05AIv9jdJ1{MJXv3J*0GE;t;8}mZ>)B%A$PR{ z7=!3k;xqY#S!3y+hwG|ehKpt2@6(pkf31&KyTq~;W%$7}C**a6yN10IVULM< zs`pASLocqwi(P~#o9iRh{Etm25l-|bUNsjpH!^Qn5FC;j5-lzpnl17zS{u?E{4%6H zIQ~p>FmP~Su=<(Lvxxy!x5sYbw`Xqs{3utlhCAaC-^soF z0C(S$;r`+EfsX#-eudmEv-Ian*-wk>HTJ`|^Uj#J8t7zIXw2pHAR<8TlaqRioHii^ zIBD-gZ-Lsa+OS$lzZO4Zzk7a5emk4=#4CrF4ta?6L{l0C8a|qnG^I4vr`IB`yk~oF z8*vy>78e`m7%Owmn^OsDA$dlkM;0aBa9dr*OV-M{-vQIjhUkDANPhM38#yIYA*&B5 zfp9@1yLq}!&c5l66ze-{cvhW2L*Tm3GLW45-W?f#vtKMsj6v+pl}sy6%VukjiiK*T zs|Wfi>Z@yDYlI`KtEuY%Vxy!qRqRt$%k~$Wb(@KrlDRv*fhZSLYxO9?XCUJiujAXs zhEGRfVYy*#+1>E|kMAHKUkQ{Z`OYbPmZVrz$TKx?8yvM;&pr58s z=*;zP+1TAUKNnZCP=m2qw9QAIs@ANoN16E8`o#N(?Nc!ZGG^!`w}ux!b19xv>vgdk zw&O25Q})T-Wm|7s1pQ&O_ss|YzQ`%o@1{%G30SY}NcE6UNAEcL-h@c-tKh_AnPV=> zo0MZAvmwt^Ti>PJ4OGj^I;l3K#u(Zis@15(=c((^ww?H9}96QTtL4Mt(4q=MpkrH#8l*U zK6$icBO?G}kw$8a5X#>)OO5mnhJ4U0w}YllKQv9=5R_He(apeQNT3~RM%E2Fh1LT% z-rC*0*=YOGI?MV+l>;t+u7haaSDats2U{!?zp2#Fg6_Us8D8mjJ6QOY_!V$KjkvR| zQIf^Sw49e)L5tP|_86^ZS((@FMWaz)gg?u@c3UneHsICA3vlR{>Wk`4pl4AHyD8nU z6#Kw0_g^7jbY;8uT69@_8(th8tL71giT8qeM=zAuV`&jO?eUMUKKfaA0}eB;CcMxu z*T2bWGh&OpQCYbV=&iS<8wWG+FLhoj`10wE@0+rpYo7-!3e90RcTi;xE8Tu?i9IyA z5ix;}1LgNG$Y{%M2X_-|Yp&%!d-dae(6^u;ccZ8RyW_jN#nQx(rJf~u-s=ucjhz$t z;DN6#?j`8L>PHt&L|%%#5S+fdwDn#s(m&xv{Hs8hU0V;@%G^BG1)FjCp_G1srTE1Q z%IEQ~s$265r~Nxuc6bMrEU-4;@gww2(K)P{hO?PZlvBJHe*BmjK-BB|tvrcV!j^*T z(5vo!BO7x8y-$Da{8aHB-#U27(8CdL=Vbe{cE!0MTd z{`#2xv>mP>sl%7H3xJaHf|A0jLu1%g#MnzA73A?|P~AM=9=ticM|0!is7|Y3pb%Ka zZWn)WZ8$l#uvta)*k-?LX2j<${s1b$epUL&X7ga6EZv^8LpbhVMy}vKZyv%p%be{V z)foq4C&d;$j^8}?VU?^|dATC0@;GRI@7(6va$&4dW)Rg*_{RRfP+R^l3y9EHX?A)Ns*=2wvj68XEk!~dqcdZS$67;g`ExHk>} zqg~J_ zKWuc7?st$lHv}96{@oXV+|j>pd+`{gGe`n}KtK|I9Y97}T2dNt2L95Zq@6bT0I0t- zhzw~C`rjH Date: Mon, 27 Mar 2017 19:52:56 -0400 Subject: [PATCH 3/4] 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>EQtfVXTwg{SM|XI{RD9(;PxeRJVC_15vMFvvXhldy$1ZhIQk#(=z2C|7lAGJDjQ>l0DB`?$q>CicM!d@lsoI z{)d!nCH*Rl+y0yQoJ-m6DvPdxK84z$olBYd;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~35g2sP~O4T zJkAAtWn0~|&r=_uV_x^K6@z{R2Us`G(Rrq2*hWHEsEs%CT~0%FyuQSSmV6x(+XnZo`I3Ft$xjThE-bTXy9NA(FN2=I zIy+{88FwUw&%GMsIb2Y>=}*9Z3NX*#M+k8om4078*LD$~S=a5Tk*jkk5oSm|BTB@)h6P9Nhs`DS($lIsg=$g__ft!D z>&Cqn6lR~BF63b`|aA1H9$l|&DkpDYa8WW8`%pvLFzQ}TQ7 z7$+Y9F~noSiFk5f0p(L@oD=+|#}@TU!rL0EK>mV_=(m)O_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&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*m=ktr5UY-*g znmqE2pSJA3e#+-(+`|p+J(QK1+Y7-bvbKA8CAiu9@GZgf-IzLk#=m`ij)b&E^#m;n-gKy#qhmM6>gfK~l|)35+=L>-Ks{ixKSzMVPHTmm-8X)49+(6>no`~EkFw3gF%gV=&atcrs~Ddkm| z;#jxI?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-`}9=omWvtXqkYJ46g*)vQEneitXvHV{~v zNw`+vALyInxv=wnjC5OGrT~ITxsC{a)FCwaSVeI~`?$@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 zusr7YCb(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~gZX5Glvq`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`jT(%dtdCyOR=rZgZatjsY(sIA1siS@26{@bQ(#xblMIvN?#0Ro zlm)tg8#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+dv17mycDs?F{>dG@G)@NxO-4iNo!lZuOxh#B{z^rt1(M< zVNF`bamts5K;O_9>c}YWn3c(a6Kc%w$?qB_Cc@3l+|+GB>28hy6?2p zQBqqq{Z3DI5&6OwJyeN@7_KPo8>7K!(xrN|>At_~} z4W*CQDJ3(DTTiR^bst_NZLG&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%ym5l%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;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~mrS?rC>Ix7 zpTowrU^q1SKtkhWoluS38t?~+4o;4SZ2`1XBU z1J#}-_!hfio>u94Oiv#UNwye{7qknx&iI1y0PR&+oT&4jY}GInu2+r%^I-1v;E4h4 z`g4(bL@N$}57BO8$&(HHGnq9-re1yivBw1lzW;85;S040gfpkQ! z#m5*-Ju|2xj%Nur^_&zwH?nF12AQJs41L~QREkN*1Pu?T5>73%EV;vHpLH-q%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|uXru)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;&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&*ftD`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@|~Xm0>Lkmi4rE#naE-S&BO;L zQWpr@IQ8j~O0w2IJdrSy0L%JE0%gTx%X7vk%DrV_b%+jyl@umvM zWv0%QA&Ve59VU(pqFp7JU}kiUULbhCW+In~6ebdwsA8g+iENMY*BB9IhAbw=nHkiX zXkeBT%*?#c9Mg(=#}dNc@g*kHy z^9vHpV}dzmFgudLX+c8IU5d><19eb(< zP2C`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~{IRkdJ( 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$f;f~l^y!bU-q!N!Zl0E{ z8%qE|0q=Kh8#a#&jYoiR%E8LQf}>JAC1Au 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=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(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?1hsJ|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(<;`BsQ6=?Ai?{G;L323~z+P-#%&KK==^^~&esjlM@ z*$?VkCu`hpc)Hx6Jv-<5{G4i-xMBEo73_tv#)<`59H7+I|YH3(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>MdqwaWU81(r0;!3G=^-p(NRK`+bpAjhuj)|QM9|6A^*OI@B1v< z^R=V@`}Ce`FW=7w_vX?igN;E2Ja)%Tu-DxBMMb5gkb3*xqh6Hc0i%8U+!`-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>!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 zsSPO^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^k53nI1shbKut@y4G0k?XTSEE^g)a zLz}`0fVq|3pJot9Jm(I9%YxPVtu0)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+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~+>fYmPY6#YD1F@MI(9)>2oqs{h)k@|D{jrOhFOS3n_J5Gq<{aLQS1B%ZNU9Jm{A03;G=h z*BJ@COdzu-P2$KEj^+hW?SeWjSN>4SzOaEJnCaD8n?oBV$GxS;d`}$HnOMgv&PhZ52FrGCA z7( 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>;`O_#<8sKM7gO`tNs*8Y zC*KrF<0}*%HChIP+Gt>(`}OA=HZ1BrgT7HTQSVDba-Oao>R&xbukRnR5hPrCT-Vuj zNldr#NLU;yGjyC^rrkdsn%#XI~(RCROrlR=U4jI@3ryYidDoc zolPYs+R^JLmlnrFX=UqnOU9YCU5wKc@Ordod=$J)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+*RW!P$x20r>5Ta#M4&SAA;nIex zxZ?(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*M9b8VQXX_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=dxV|pB46l;Q5<0dnP;1xh?p145mrv7qN#7! ze{QTXQp~u&n6ydWYF&@v)^F3B#d{A|iUM>#;^F6hJpzv;>)~rOyGQa|vULq!IVHb_ zpApON7T$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}>B5OC2@y+3_MVy*6?JJgeCV^Z8Q_9u10z4=Wo#FB zC-3^yPG)4# 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^Ta5cZ8AWTj2`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-l2AQT12D5CN?cH}Z!&r^<%O#)+gclGZFc~dRDB+>zlRgp+0(zW;?$ z`0Gx5Z4VtriGuJHAPU~O*C3U3vkFC?ax>VyQ(sIP^eD(>5+w7d<0|*OvOs}Pjy;96%p8;G!dync-4;I?xxM-B-xjvcOw$FEDXQ*`6$mhjK zixj9aLvf6P1qFUBR5Z+2_4DC0;Qitn5m3J(+fPGgx>*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>gUcCUFRxNK6X 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(< zE5Ha3663ycfLEwwMijg)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`EU2XRB6u>W`I zi%w)vdSWqRD&@nuj|Mqw8hm9VA0XKQ(ddF$M}3}Q^Ji;t%6gPlf6 zR9W6@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}vjCr2QmnNu)aI zfP4?GA#5xQ3B5lSqWROg3zI z5|CYfP18<%DF(4`-9ZyLR_69dC@i-va|2v@bKjWZ`VUkmj+^0V9-f1REh>`Gt5cVH zjGJ33QH%_UTF6~b9p_hQYF&_#pwnG|4@9N%ALhaz2E_ z6wh4oB2NAU(`(*;TeKNBhkc0@oeJiU{^w0!_K*nT>Z}_ei3)1^-V#0y+%oEu_xMHwPkz6m6%wfn|M_%-bwU3XaI?b*^Itp^ i68@`aj-Z_eTd4lhjYJ`kG&#fjc_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`uP+p4WAoAPys}^mC^BlJRBZ~J`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%9mbe4uhjp6VmqVP49Fm3GygmOF1#&_4Np1;PJ zx&Jng|3~6qZV~7yOd)d&ug8Y0YNLxPCVmZ8^gqphM6${%0tt?hIMy}M-LL$)1pbq% z=0j!yW*%eFI4H&drf};=MLeosN9_~- z?bm%GzsL>k7idF{JP!e9LqvW_q*y*}!==CNz+OpK-1~2qg6hsui}8XG{vuHWY8L_Q9NSsoJ&Y*Ww?_YpzBVKL$&+bbL~ zg;|(KvONfG#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|>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;{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}ufUsPiGLec`Ix$!{c*7MmjNe$RsI9kwvKkT mUuiIV`>!<^j1;Wrtozlr1I+d~sb0@Z+16170`Us;R{0lG_Lf`# From bd0573e5e0673dbf2b00a636096dbff52895b69e Mon Sep 17 00:00:00 2001 From: caes Date: Mon, 27 Mar 2017 19:53:17 -0400 Subject: [PATCH 4/4] finished lab 4 --- accelerator/data/plots.R | 8 ++++---- 1 file changed, 4 insertions(+), 4 deletions(-) diff --git a/accelerator/data/plots.R b/accelerator/data/plots.R index 98eb551..47ef661 100644 --- a/accelerator/data/plots.R +++ b/accelerator/data/plots.R @@ -39,10 +39,10 @@ minor.tick(nx=8,ny=2) dev.print(pdf,"../report/nuclear.pdf",width=6,height=3) -lines(c(0,1023),c(7600,7600),col="blue") -lines(data3,col='red') -lines(c(0,1023),c(18500,18500),col="purple") - +# +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)