From f509bbd15dc95eba2118d9794889231686b69cf2 Mon Sep 17 00:00:00 2001
From: caes <caes>
Date: Mon, 27 Mar 2017 19:52:07 -0400
Subject: [PATCH] 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<L^o6xgoD9a%*l_6-&k)LSSh|E_psH5C6YKe=Zg&ZYZ5jrZ|(GVh?E
zppcM|u(7e@I}CtbjxRw_P*_T`5~3R3tEX9!zH2;H0UT=s$+YWETSSo8WQ5$4hAf*8
zR#xp<nMseOF;r>~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~x<cDFe<OVwG+M^mN`b#$`hD
z3wu>fz<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}<l_dkzN(x;xlmU`uXFwb8ZA*ax7`#oS`4ROkEF
zRZ53;VtXEY8Z1K?Z!3GLieb`<Bep)kd6*e{ncl7s(BFC}+Zk}M)RvEy+on9IEi#Lu
zjBIh&zaAeR!T<t-bmm&`$je@(d>iBX*j-MoVRB0ltSy&BZav`uP3%<o+o7AaWi*e&
z(h1{o(aUn_ZqE7gUj0X$MbRTn3x4&^-&Mc6a8l*U!Zud5DrY6_Jn9kanLpV)vDI9M
zS`Rr!T&j3EmnfFPLfUQux-F{}byw7-vj6{P@yAwP+Lz~O%%|jHD-`_R0lcW73b~w6
z|Fui)5EQMjqI%ND2cP8(!^e6x7v5zQ#+*>_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}D<v`bHn
zF>k)^h;Oam!TridF+iv=rx0wex^8nDb!os)=o*-0sP+Vfe8$E?o13V$_Oxm?EhszS
zueEyjlfjj0r0A0z(Dv9=2(!Q_NvUzXi||03+aa4+&NSVe$bf(!*QMbf)m<TRGbsW$
zADE`ElRW=hSI@vLBqHS-b3If%nx{1Rp|Md9uPbZ=boL`w{c!GC{?$?S%=mj{+vLqF
zA<0i-Ud-oMYa>hkHm;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|8YXg<nX^bE1x}+9&(%amB0c9o!8%rM&8K3t8bXaoD~DomnDu
zk<+S(3QfVK_DL*0715lpA%p2gCDLIL)1jZygnIsg47)3$yvac<b7DB91DOGvqTh9b
z;|ntd;3Ya+b(A2(hdbAlt?E&c<n`fu8Dy|W%*sTaqhhj0_hZM2idJHUp5BA%pbNh{
z=qG*x<*Z`X(`Hjhe?YtpH;O9HPmAwzHCQphKY}%7YH>d3EVFi9ea!J?LQj1@3iMf0
zgc+05K=&3pmf_zoH)yvlNn<1+eKmidTVZ%e4+JMg?%x-($55@}%%{_{wZH$Y1<WXV
zQ!0SXb*S!RzBaK~fQ?onwQgOF0=V+~GpDxE&ynu=9Z&ghv@&w@WN2`uIfH8i67{a@
z2h9p!?@czuJmgn7;1+LMbR!EJloYAuaz-(tNWno=MS?UWGjKytn^@USRFM&+HR$#G
zfJNOfHg*7bDBAyY{5M#{(lRSM&6+!}(9}pZLiaX6()fP4f)u;ET1Y|9q2QvZeR%G4
zBZC7h{5Zfg<p+-HQ%sp~xJ>^2)8FI7nsePF6n4qMxg#HN36A;IDDpwGe?hY%QK>t&
zOV`DE-7<P&JbQr`uIQGH5^PwBE#-#|@D?nFRJy4Ms$3Gdq51c_q<`!lzGb(a)_G8n
zn)0+JP)};uX1vuc+#dE><*Cwll0N6opX_g*2a@6^W7ZsuO|se!&4pK#Z$QugfY+sy
zO;GKUVQ!pa%Kl|u6g<Vtrz0z4Bd{)x3O|~i(L+Rm3Q@Al`8G!lS|=tyfD%yUGLvFE
zFW>$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`&<a3M6Nu1xya}3<eH$K8*9;$9c~EoiHA4&Ir%0X
zv$t>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@^$#EN<nM`eMYY)2z*?;)0q!%6uX#v*rm6T)J>CA
z{diZACBNl0wf(AoUX<AMu13CE)jkVc*X^#GGEC@Hwn?>iIAH6`*3CM7=INNr7{#zs
zMA1mN+Yb=D(oDaFcC`lkHL$|!QgCpN2?4ho!dOO(n=e=&_CAz+rhg%a=<spvj%=0`
zX5DO?j7B3?TJb@8D!DuCd!Mko{aym226UsGy*M?V6CpH@u`9LTIHj4an-@{<hNLT}
zW@bbV)xhO}oMJr&%c%u9qM`a3dc%Sf*tCgtgRM@UCvWD%N!x|ozNbY_%IX}As{+3d
zTwZ=jlSP29qc-L-&fgX~GK1B<1ZIEz4weZJOWsj9BifZ!P<I?w8JkzXf$@JlDy1vn
z;Im2b*`jnCjrf6wcQIF5R#`#`<pH8q=k}m8zWUBEHF?}&k$}aj{;jJ=S7-ZrRkD-#
zKvc|0U=v@O1{Z@*?oADjCYy(($ui)F`%LR4Tg4wDU~I}<LZFjNqFT7`Dv_s4>Ugx!
z7|__1AGblk|2a*?D~gu(lTM>(bc4mb2vjI!0%Ifsjo}j~v>_^zIpU}wqc=JVOyc2A
zH>Xq5=xTm0(}<u0YC?=UH=X7WBWxLnq?a8BEvS;5=!XL`)@@{__F-m9#w#$6UBkK9
zKL6uqg>@&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?o<KuofFf|)30xHMSPJ}r?iB&7p2q-i4Ldgx$_
zIUhSHqC8#=HAdfPPjEAfSsWp<lsATyp!@n$vBoObr=5BFI##EV9}_JaKk}@sOCMhr
z=<QcUIf_skDCFto@PG4pP;r?PojbXDpSuY!yq~ypk;orowVXarfG+}uB6-K-N7yq7
zJWT2o-9o-kpHU1SbkY}wsv)*V4QqJ6?bds$jk0$$2BPL;bsrTj{4@XC;s1}{|K9*)
z`6J&QvEEED>OIF=#`;+^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<PZ+Ndz+KF?49<Ci*!h5Rxv<;zPb
zS4a)NYam^K|Jh-$N*YJz+e+Gp^5swW0T3ara8sk3Q!-(Dwz19^zMDQU!slBK?h5V~
zZS!TFv+lU8=th<rr!jF8DetwT#aaQ9h%~BJ&DN1f-rey(k?3oHe~MiHI`LF~F+!m#
zYJ)45bcGy#-K@QL#3MEwkIn0Qqmv9}Y;tZgStDM+!tyiQ_<7g6@|a^a#&NyRwT=Dj
z6Wf^JV+^Ax*flC2Wp0m1mLdd#L%YwlMeBV1hDfm?nb7<hlN(E2(Oz7Tw>#dU+i>wq
zE|+ui1hFK#KOfeXQ)m<Dx>D<<KeN@NMzAz}h+w$^ZN>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$S<Gu3FX~(R*t^9FIjL>m!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?lF9<ubzq^Q@ifs}5{3Crj?%eGKpSZFiBSsnFxUrz*l851+*LJ7gtb>J+n2=R*U
zJrcBRonR_=8k5OsQt!$+ajI7@SrVwsP<;UF=7&5>VG58}#;`XVvOjx$-f&=hB5v()
zF!bPa%eWk4o$>e6>e;B`rv<EgdaF(s@jYQ9PR^2!iRSFHOlMwNCT`6-VfX0P7R9tN
z>+?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|s<y`$lhnkZ#XCDn!(!
z`u-G}Wh2!k7%Sl?+)VJ;QGCJr=ygo=b8#!ri@hX+`Me27`Fo@H%bc4m+I?3L+%VG^
z+fC0{R1if`EKr-WrYjKb9=E*zo?n}K^n48w{#01V>2|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<L5b+X<>#
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><ri`P{h8AchG$tj<*$;k8I$^
zYMRolznJV4xD=r>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{EpuF3N<XU&Ol+6BkA1@1`z|g=mi7!O>9xc{i)`+lc1&q7o+WU%kIPbVn|D
z)XMY}R#^~3-bgfx#J4lO1Hrv$N{xG=;~lfa-tM#C8L6gyt-=X{kV%(dn$xmRBNw0_
zT_Y4fzl4gy-rDxitCKe#<xN?&WWk)fJDEzzC{U)`Am?t3IOWYg6Y}r_feL(aimlOa
z!wtRwWte=C32m0c48#~#dBO`}bp(gb@5f!FzJIOI;vyr3UVNpq35b$S4V<6+q<wPx
zK*BI<fFsnrRgDA*LW*@PVbtg*`~7*hyiEy4_)Grmii;Oh+#JcE;i~sxt<#rFZW)!`
zeS#m~RN@Z4#^SsG9ASBkkFKZglL=ZUKetS02IA|ppzqks(r-N86@BR+406WRzD)Ea
z@s{XqoQA`?gPz+^oZf?w-f64#LB!6=;&$9a6}3-BjRL$X%P0EYTtfD@U+MMIl^}|2
zt9}UzOS*gXD&H(Fl}oFlnz7Ac(S9sH6XdWaQp*v890KmYGEIIY?K&d4uMe1XxjeJO
z0+5b;dT96%FtsG3CBcWcOq=C#9mLRQq1ZhPsIhEOJGi;){Mh^X5(gD2D9O!H_ZLr8
z6|~(8fm-;R5q?hD)q0zqeg=uc)RI!?+U;A*!`6W>>b|B020c6)JY?f26B-_VRnw)U
z<CA%~5`Di;5?u3P%YtHKz6G3~QH8Nc@W$_P`BUe-dkBQW$A36*F)^W;^SM52Qoqvl
zyFycuM|eB@6<>tn4HOhE2^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<NgE~~fNdLN{%EQniD?L^fL#C+2jfOCD-
zblo6$M2eIK(HB4yok0haRec|0;fv>#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-(<l>5tTAPFXa2kWQTsy}v%(Noh{@N!E8$ML)e|
zp&uKZBF8ljKDU4Gxt}?hGKA%u88#<EnWZ`}92<Pw?!ugVDH>RJ<P*+&@4=e5-7%;T
zRBS^Hh-O~b4IZ0E793;u{!&z>Pogq1jQF+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?z<Yofz
zJLWzyy-x?9<zKOnAFA<U_j%P>CXeFL7?@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?KMhU<c{)c&tO}MS2tH@YM$XlyG!c`!cE=vsHKF+;qo=u_{ew!h_N5i;1WxkI&f7
z`}V8Wn8tOtaybsdE%RknmgB2Ue5;Ng6F0`qCcoPXJw5k1y?LJ3USA$d`uSJ}9rO7@
zS^yi&HPhtX1Grn5&~}0}drT7~GC>hmn)PT}<itV3UgteiQ@eK><||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_esesHhGU<n4epu
ze1w%B9oc<9CE@gQ<NA`t+D}KnVAtk!v@+bnPpd#0nxTNRu+XI1gDOxtQ5Yjj+oT^t
zq>y@+EFE78y@|##V>FNPb1N}p#H4DQZ+Z58`}l0Yvh?*}5o<W>J5POEmkejYUSLc?
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?sOiS<eFdey6+5jSkr0)^2m(3xwuQMn~c0t;6Tg+BJH%G{P+mj`e
z;e8+F>Acp{S(3_FlIW|CO39}ZTccde%~G{FwA>JEe$2<QOJR=U&6^L1k>Z8M_bwm}
zQWe2B_*8Png?H^V9kQn^hnbD0^<SHc{rM`r8??O}^Y01G=8Bq5Su<Lq+JN-H#*z7e
zo#IiwNezgQ>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$sm5c4EQ<Hw=JF0Fltr;wE73AG-FYQD~$I1nx
z*$PwcvE!+@qxKbpylT5gU%_p!$=KGk&aEtSSHS~<nMAmx*3Z7%u_}c#GT7sFl{|{5
zz%N?XOq@{Q<yOp@U6?{1f<t^ZVR`s0+v?=>4Wih-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?2pzj<F{E{wml2gDFD!s*gekB`0sBM|ognj@;o0
zT#%f}hhm}|*S_>N6s!j*OoRthc8qP#<Y)D-^hbADxA?IEv*bnPULw2&9cCB_6_Ifg
zk}~<^l&TmB))Id!r9@-~)GU*=@651tG@K0TM;gm>O<J{`4D{k1qbZjLO!Bu_OBcX|
zwykekQ?=(8o^~om+(|;c2ts{uLP@{~Ta4%&b5y>4w6pv;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<?y<(s?u|BNg52stnutHtA3W%=9
zAECfj|A{v`Jh#Uvt+};Hl!W3Zvtn{rEbRR%Hw_hg^jgKkY^qr5+9D&5B=2oS`v@bF
zUv0<T9G#EG-854_#sk~@r;M0-7j5F56paz-C)Tn^{ORy+wyQFa({MAGPzRokUr(*8
zmw8YaS4O6Ae@nOaxysbf!AR5N10_3e5<ZkqEeHK*J4Fv|&uqaVzzaL7im7y+rS4*>
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<OJ3%;eMS))DBv{ogN)F{adCwFvp^jz9QVVzI>;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<c=B;
zN;YZPqpCjKz5(jVa5TM^U_<pebJc1#1B6~^_lN4jm;nCIgcO=#Bbpx|OwE1zeBmuE
zjwNr!^N9|XL#CB#jbU=xue1t?!U9GQ;N5%q>#_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(<B+NosSy?w_-KTX^iC<(v5cyHz0X>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-Bm<NWWuPf!sXFX`V*(Ep@LkpE?B|G;n9
zpjWmJz+a{RG*G1fQu>7t1!d;s=x*WY{?{&5<e_1(ps=9+N!I$w{pSbxYrp?4k+lD|
z=KuW`+!>L$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^Q3<R1nL`@hyE#R;n6qJb~t`rG<{00z0masU7T

delta 10124
zcmZvC1ymhd&nWKhEmqv+;I753I0X&|m!ie(Fi3$S#oY?U-K97!6p9viC|=;;4h0_l
z@Bh8;_N}+qnw3c=vy;irmSj@nzHE>6N(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|B<IpYL#o@LnXp7?oyO&NiFu|w%a8vTfh=dtG}`I*
zw4I1Pi(7!}rR?_o)js)p77h#M(fjHZOYk8e*e1er1>C%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<Djnag_*QVVl?)l^=;H{>`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!<drHN+H~02mljfT<q`=PfhU&Q1e?;=DtK39-V=-&>w(
z`UmLf14&?yd1Q+9+V*iZ!Rf|{hxi|U^HD5vk=$0BmXCawsbMdD*NXL$sT)QgT;Ayi
z09~Gk<L8gR@5pUC%-RSQ4SiE2!ISwamL+aM*sCk#HmgLGEd(t?6&L)Vyu_K*Dk8()
zI|MBZkw=APzn{Vc63W5QZ|9(fSN<>}dxt27pcg)=-U9xpEw-{(LEt*7k8dV@#fAEW
z_OX;8<~yOO(8Ll2R_!HR!TyT*Y$z!?KrAQvUR>#?dXFAWx*_wgwLKDEqIyew&B`6f
zFS^)MMy-lgsYP>EiUNgO<Ma8)1%^67l)-!)Rpg_#-pH$k!3+YTc1pf<r1$Q1Mg&Ss
zcL*~H<UbGBO0hwHo_n0i)rZfA2Tn@lsc??8br{Y!2@>agTd7o%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{*qH<XhsU<jxRBJLDYq
zVy@s%A7K%3>eNY%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;<PVinK@uPWZ`VjEF<e+IW80SonC~zPbM*VW{NF
z`OWD|=N$<GMQw`a>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(<iIz4FFjbIZM9r&$oW46`X4iDR+E-dM7o
z!!N;?u5;{t573bmDrpYBk3ZQ<@AMKU>@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><YOjnuq?@4LUceX%-3L@zF1!
zyP*;rbO7YlE0VUiwk6$=hd*3!dJqV*8FI}Wg!VdY6AvY6Gpi;um8N!YG$v<$#jd3m
zP3wTL&5~}Z<4&=ibAV%d4x3R9-7PBHo>Qfx15InTTD7E#tr5D<;;FHpo0#rP)_Ua}
zbxLN-^&%{F&9<!)%u{%ILb5c4oxkz3itG4|p96*1%i&_S5jsazT(V;CDLIOn?}8Z#
za9gVkcv;+!$LtQU95b<pK8%kzO0l64>oXL>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<AA4iu#wy=uMwcjcme|yKeJMpSz9tLmH%_Y*tH1
zdAlc$*VN+c&D}!NL)G>{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?WnReC3<g<nSm9yk;>C?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@Sz<SUdOM<7
zGMZwm{YrMRqVL=bqH+OBHn*~<`f1cDNe-L#K4pfKTpUVZ+N{}$jbxsJ(?DPA0-hnu
zZwz)@PO0&QmIi2>Vk+$m{&2_(^&shKmt652xN7yJs*3#Bx(gr#9O+VVFP*h{oIBxH
z9q;YQhfqEew#i?DFjPHC{G-|O<IOD&haN6v1cPB{;v(GY^;1K4%U`qweLDIAzCiQs
zG-CP3#qLe$i#wJSI{Pn@doxJ+bt^z6nzZ={1LVpkz8{;58mC+xy&X;tOw8h>VMS9Z
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^Ii<wBOAC1{UC)ooZz$qHG}g-qD^`4KcBEXM1(*eOVj(VHc_IHk=*(us7Mj24
z<PRGOwD$BGYzi}Ll6M<_XdP{K>p<LQUeN1+&@IZH1kkjRSqV@gX9H^0Br#uVA&LN~
ze*h(Hx}(k7D79A>bn9An>0Bd(yHEB~HJHh0N{zYw<I==A*M+;J(`mL*!V#XaGvDP<
zXx%&_@Ci)Nd{XrVXk=;Ahwa~jd{wFV22v@HFvHukp=ut+2XE^DuE9H*dWTqeNP)g{
zD@&T!o!SmrHK@)$q(LXe3WB)Jkg~k*Yo`G_wphckPC>M^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>!<HY{uCW8)
zuvhf`0rRq#_^~^aA&GSlCYGCAr8BdaAz+p;L24+AYHJ>&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>PnFU3<Ql9U!7SmEePoarOfq+
zX{M~|ahE)vylpEs;dpJjyK?Jquqt<%>768g8tP718B6DaD_WUW^U?lfz>R|nrugW}
zk<-;+<uhNto2fIm?)~vU5yD^F(!UV4|IM2La}UVMb05w4Ljy|So|=-J`pYUA)8K!S
zd1z>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
z<qh1iW`v24jJiMjo9rB%`^&J|IVMcB-{(|VZ)N)hZ09J8*l-g_S(6#G$u**%({2qu
zZyGb=lEy&*I)g%aLVkbneUps3ynnt(9UJx$4qWA`d(<$d?q_w9DS@7m<adl_qDS)b
z#+gMp|6avm&=zaFF;!EWH4ohV{2R`=u;L^cU;a+CZB%>VuQR$)3@Hy~@1H_M1brOJ
zW$tkP)#I8ubdSm5OX8I%;?mLpt^PY<`20PM^B4Ay0C6OpKgGLGRV?Sxm!g={&+_kX
z$-AjPM^@2De<bK+=hkR`eyzTm;bV^~%0Zh<_jZ5lh;k2^H8kKuw<O+ux_s$$zQqr}
zFJIXAQYy+90@-O8XZ6SLh<mb$<-5=`whn}85x-RC5zuVX?JO;6kjf3ASEpiEzB|Fx
zQdT}t&<0}J;%bHQLzUVm?_OyRK1y81#`Ym_x`WPK*j{j=ME-V^PGfgz;V};g?0yOO
zw8*CRJO&RExi>QVV+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(<H{&-}s2uztww57V(3eenn@J
z!4+r%!AY8{PIgxHOdSbP6CNUzf1Ef(h-GdeXbnV=G`;OYU-*rPM6INJ?Kl7?PRP2f
z!HELjx*b=!e}0{2`f%_68)Bt0x!b`h@28TmYdB<#{HDesAa!rh+emuG<>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-~3<Q0MoYv`BDpWYlnQe=Q!cXH<Z(
zquXK&{=*6+PqUg*lm!BZN1qo3byLdP<h(45m~8VD*~L<m=J6vcO-O<S=shX|+N4NI
zL3NvfBwKF&-%)w#J%g6}4nt7s`)pRhBkMiLW2ek&IAH?~;fOi!c7>FSP``FT{jILQ
z9^pJ5vuCS*ii`bb*OGQp!2N}f<tHEncR)}5eI{g$Rz=<}PT__b-CB&Ac#S#xtYGGO
z=|sLOyAK0bm4)MlmccaVwQgQHyQNWUPx`|^VNIF%Xs`;{iUw0!cbciR>QtACJdx;=
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{<z+
zwW-_Yi=OO0Rf>=m)_`<5WMGWw<XG>7`~6qO^<4+>4=J42FU*cH<4@J_BcwLTiE0=!
z$=CebYxlc8QhW^C@h4#q%^$-g*GwRlu8US>WHie!o*KLs8;NA(CWIb`&VTvfxYNhJ
zf<^hMktExC<HS&tJKO+{?p5zMxn8cvxi2$|XF;zE<06ZjSzjy?@TdVwU);HRBr^r}
zC+WfT7U+mPWWV2x)+_QMzbzGbbCGN9lC_dI^hQnF*A^byyU}_}o;u^1SoDoO+%m90
zV~0+zt;5N(o~`5X<)`9}?{8w+<fylWS!PBK`f*cid3)>7F=hD^WJr&XYqn7cFkQ<(
zbHm4&ss_iATTp$wS>ufWPDZwZ#7I4<v0~{AxCf7M)I*S1u>)I3vgJ-LdXC3>#74Tn
z@rk72A$qeiD`}@l2=FfP5P6@Kgfs;?v_M~cOU-oX%t2Ad<|9sk$&ynFZN9~Xhsp80
zfB<FjjwZ<aJ6AA$TJY$$cO}|Nlt<U`<~^<A13tF;CvFu~=WA;KadzO+I=_*HuzbJg
zlcL$%4>I|%Vx{p@=b?OVxe?Rk8wOW+{r9BxmIn4^q#uNJD8?e+?}R9chzxljm%sg%
zcbzK>a$mM!2zd6PPvl6X=)<<0U6=W&^<te`lhy84q?wCcqq0Cm)Q(t24kFSEhg<XG
zMQoGJCB+>%K1_32V4AZI`N&-GaLnnAFDwVf8GLb4#Wfl@d^((ls_is|#)D4swst*f
zJ#OiBbiO|o>TI&Ou5&PaPiIeBp~Wwj(dczpQE7s+e#2sm-KhrWa{(<MvYt~sX=Jj~
z!Q3~2>&&u<^6w~G@Kj@s2dy>DE9@$!UX~j`inz~(GUSmPfKPgC?8<huO&>hmCfz2r
zIA4)U8<!0Vq;<u`rVje>sU>|j+&JU77!sW3irr|nqO#l)TDP>?W{pVLEIFX9nb|gz
zT%jNqGunpFS{k<Z;UA&<d;vncnH+i6EFfbr?UE8xoodeuiw=GpFzWO6T9q-=R3fUy
zSm{F+DQ*9r7+`!U9Fb6qk7^?prj&{F)+EeCIT>KX-aRP&MJ<*xJH;K5V7IG__gTK-
z<wBXhaHAHVVBqO#?oX+LIPIwDt=~zLEuQNO+mU4u-Xf_V#Amr(W3Sh$JWW`qDOlc<
zE-_c_{D7KR;w`V`GG^K@<kF4vNM@b9@Uhlq8mpK-yaTA?Y+g5?xPLOQ?Dh9bge(2A
z@13D75;2)e)?|ni7X4x}Z}hfB<pU!u8P=hGf&Ms;@fccZ+VqR>5NARR#=mp?$}c8v
z29{g7V|DPeY3Ao;s&TA4X|6%m8SUIilIf&ePBlN5);VY5q~R(O;kY{4WRv46o#Ew8
zkwOes8y!#t<AE+@hZwinc9At4C6sp;lAM_o!>yqnoqxX!oVOG5@A?{lM_ILsVJ*%D
zXXLws9h2X219`Ywy$|B`{^}L7hvqW*HD|x2sv7k$5@qp&T8<svY}3UrGmdt0Z3-6>
zP4X1pu$p44xUxFY{2KfHN{_~4iz=dCl<G2N=|kY>@=>u?oHNh1sK9i6amoYbG{&c&
z^<U$x-@VyI1#>_oE>uX&hC>3JiydSiW<ap&&7-&RT^!+lKy^OlvU8FD6BO^79Y1}v
z>|4Q~Pu;l@Q$#rhFcAEamdfck>2atFLaznukKc<e{erGEwe=rJ*Dh6FhMWd_XItT;
z0;pfzc|@*@Z${3X=ZcxO;{7MOX1~36HV-!5enD-fa3#=LHWhNhr9Ao33$#@DE~Dv)
z`{T{PG*4UQM94f>Tc(;Rp_%z6@46PQw)odYWvD@bJJXQaj9LGha88uRfu2s2fFCZq
z{Ah7;9GA;n^NV?iY~8E^;sRvwGT-BX5}-=;F47Q3%JZz<x9IrvD!Co(Xq3Laf{N6P
zd$oIqOI@BO8J`#Wb)Tn!EpG)lU0>ccd1fOtiQ8b7*~^&LBiT*pU@sx>P<j!sr!dj1
zKf$oJp14;=bHCcP=8F2BA3Rth3w?8EaHNz~HX+;Ch3~0fe<R;|E8lu5E2!V=2WSny
zEjgf>GZ0lhL!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|F<ByQj}jKZtaHIh;}y7`
zN#vK*8nZpb#gQD)!kTQH?6h8{Qy0O1a9ed-6ve_Y=Lwx<$}rZ|^tIGURxVR!&@JTG
z1<U(lmPK;w&`%ZY-waREOO`EM>NshLXlm)Qma<lsuqNvy3+Q<4x~5tJw48OIyrzRz
ziD{!5qC$z>y$?0*YhF61qO3PyCYzdFc!V2HZe3E4T1o-y#b)$)IB&P8V^2%p2xf2E
zwziCJdF<Aa!fvpTyIFZ457jG;Zu0)Kiy3^mQqnqQhSvw?g1nCM+;xt_{?wP|8KTYd
zrLiNUwW=wEuUD$K)^wQwL1S%Lx%eV-8^r)+hmlp}gnHUuCD-@({+j|mBcdI#kY`kq
zY2k&C)ug*qs^gB1)&2;Zb<%F7xr1kYIIh2b1yCzp-DBe4)|k|be7Id0^qNZ+?{2QD
zabCzsdKoDZb}#{1y>~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<Y
z$x{crsq}L#UU2L)Dlm_C?OfX=*K;ar`PMm0UoPg<gjg(^6UKMF1T?iT71{9#>|fBQ
zphbw%(Yy8bn~K4NfL_mt6{wo7y?*W*Pk|T2DAH+!UhpS(icqv@a;AoHm}6LtO4_!@
zs7_(EP`eRD5d+Q<fV@udo-)x)+@jzYU5I%ddSt{qbK1N#MP|n5uT<%Y5=Esz@Dobh
zZ=u0_5F^b7&+{p{cG8XgpfsXoVj@EHzU-X&^iB&xa$zJP05o$S->+qN=O+QlyB1&$
z7(X~ehdf7fK0}XQ#lp*03;d-hm2_17z5ew0wNz5I{ueX47c;`z{NJ_nBYSU|kJ(jt
zx71BiK9nq1nl#q#Oedg75Zic3eM;cb4yQ%Iv!ipAa%^Y<pi-tK{>ubjzFin;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?<T0
zm4V$C3JSH3XrB@cfCCwA^ge7A7bvt#eTqtsDMo$ZQ@x7ik)2JA?9Nff-{->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!<GVwk{n5AcTS}D-&h+^5hE+_G9x@
zw}~Flx2zBpe^&8teeVjqZ%uxM5X*qcyX(;}vF-K4oe6HqiH0{I6yDKgTkhe{b%j&b
zNZ#y3{K(9UL<fh`qi;;7Y_fyv**~~u&?o)y2CuWlXt`X@E2Gm*wdDsBy*BLaE9=D+
zC=8Aph4?dz0HTMw_j8O8yWbo%k~>Mj`&jvha5md=zPop0w@B1dAztK0;)5NbM;pgN
zQMO46CkOdc%XU}oqh3DGjp5qz`MUx|Rh49&jzmjel8u3mMAgxEY!*aT>pT)r^AGnr
zg5xsPMS*i$0%2y48AX9|Ysq&UBGN6aC*k&AW&<MjD7pTOECDBBb}szWs3QJ7vh|Lv
z(I^N{X%-NOHjxrU@XsyZe<%_EP_jWne<-sc-9LbnN8leooA@Avh4J6QtHhjtp3Lxw
zJ~7l2@`}U%fm(c)e_&z@Ki}VVS&29N!iIk;wf-HyVg1L={s+cC)o;03xx3pr*|`1j
zw^u6g2t;s1aQ`HS{FMLgu}_ZsuXaF|7%lKj=uh|m{`f!6``;!+|6f=ByX(Yd0SXEn
z#wY3SJ`PrY$^NIxEHO`z0wIe%5fem_s3C~)=L7#xCfF0DK@?d3X!$QrG>*hBLB<z<
z)cLP?I5>)bBXANNiIt!i|NHX)=D$Ef6bRlNi5f!ue?*KvVPhPL5<FCYi}*uF<$MxB
zk{Hc{mWUzzKMMV!N^&Mv{lDBM!W0NzoQeLz9Dn*O7p4$v;CwQKnUkHhm76=KtEKhd
i#_&^C`9?BN^W;BrPf<4UUpk+0B_fJ2q3DDDQ2!54yG)h<

diff --git a/accelerator/data/plots.R b/accelerator/data/plots.R
index 2661f66..98eb551 100644
--- a/accelerator/data/plots.R
+++ b/accelerator/data/plots.R
@@ -26,12 +26,16 @@ 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)
+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)
 
 
diff --git a/accelerator/report/nuclear.pdf b/accelerator/report/nuclear.pdf
index e488d92f96c945627849182d06a9ceb5a10c2e31..516e49805e1366d13f76d0ab2ea5e8536d32aed7 100644
GIT binary patch
delta 1742
zcmZXNdpr{g8^=G+oVgrqE+=!!W!V_YHd-cg$t8Pp8PS|5LRRl_%bi?ib6xZ1ehm#J
zM>=$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^{#f4M<Ma3@lYs4bF!#HFNr>w&}L*=8WVGa$IrG*
z*yzD^U<DkOFn!+}x2v&VR<pIUzJ?^3bXGRRrIYu)?=;{pSGLADI(x121bxAS3M)L7
zR&zp5>ypk3<}<{qm+Z}<L-p0qx3Ib0+svcVFgrRmgL*R_5o_n-5(WD4Jcw?cOs$tf
zyd6D#dRSY<-P}+X-{~QmKly-L>}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>f<EV
zhsw#US^Mqju~e0jf!wNld#=Sw02c`Px>XfwRbLWhdzM*}JfxtCm9+`7@+~-z&$9B(
zR9=YpO~d<Usv@$UOPiuobRjLCGS>L+>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!<t5I(%_-eQ>UEvrt2JekFCC;5BM9gA^WlWKabX}2%EiUmqx
zlAl;7plRUh$!<kjxDSVdPMeV?JZq%qWN3doi~4^EvE-A3ZmRaVqho10sbEZam}FEq
z*6%*i@D>9#sSgutk7AJ)ITaxnCYh#4go;(d{>rwuctg_}s9AkjGy^pQJca2fp>!py
zEzGRP{hX&<*#$O!yKz$uAy-0iCioDcC6qr2Xh0&ezp0!-%XaZmt<oYfY~1@V%aAPQ
zsU<HJFQ$rBxgntdhz8u5r-0ti5TR&5+DG*m_D^cJ$*iDbuZhl877-m7)(7GUkXu3z
zjI`J5(@`4P9y^C!kf5=1iyyD5>vM11G_+)_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
zB<Uow+2n4pYRWP?$!uyjIP}TWi%HDMxnM^-iZ-&W89lp!PLTV{sQB`4xztNF$14i*
z1yU2G>TJ%HfV%FUxLPTg2D5Q#7N8`gA4$m<=4@&#&tsyU<pSmm&Yo)u4s9htX{IOM
z>vEf!g(N*>#-Fi=1>I4L;T^C<P}REuhz*$^ziW2{i0Zjrn)+j<<`=;N-K3(|k#gxO
zy&tW_A7(xPb?*CnTii$z&Anl7Z2rcFIrd%nVj3@E`v~`N$dz#q=Xmz?(vv{dtcduf
zm2#XmUB}VEEWItW^XHU;WBQ5RuA|F2-aE6korZC+Fpl!I%#-Ghp@lyk0v^7d@|4V<
z4hXt+=R!q|ZPF<jqZJ+_gOC9#<wf;r(>Rr}FHXFHON?g22Oe!klzF&ml&-hvjzIS3
z@K*3Jw+qy|s=ZJ%*ASWbA<vjiMzePjF0V~M%V=oM8rmu2x%1vtBrva*b8?@KSCBdi
zzaY+=ehR6OEPRMth`M{MHu9IbQ{;0kucA$k1lhlym}<R~GdR^`Gy8RsLVBxEDzprA
zVx67^F^_*`8~?G|hI%h&ru7w~a_x+I&Q^T@So}i|5<NR@W?h%9`7wB$p^N*71g+6b
zHQI^?y1LoNG!F}tf=vSrg~>hsKh$h1;M=*+SO6oRHWolhQ)x8DY8Fk_b4!$h?QPzw
z#u>J{b0tcV7r7ENWG-J^H5`#B<v#8I{Ntr5=vvgXk5PCWDY0G*)8xJ^*r+^+K+JzX
z0E($fn#Kn}w|Q|54swv1Yw_Cx290h{hqK))ISb46uerik%8R<=ecv>P2bVW#0$yUB
z<)cxd&8`Qp?;E!q6yoqdR-D)4ha4wXcMZq=bLe#KL}Nt0aAlzdv&nBgzD#uB$HYb+
z=C|`iy^B7<BAenT0-%Asr6~rZg*q6m5G~n%4wMnf82g=v#-P6Uu*UfwDyFJflkf5v
Wv@s_CzKX4^k+G319BxNABl|bnszd((

delta 1576
zcmZXQdpOez7{`+dzq#cSJ7{K0F2!$?jc`P1Zq?{?CW~1$l`f)I2tU{3wp`26l|_V#
zgyw`9L#^Yu7OAL>WT&ZhVl`XN&N=^`-ao$Y_xZk`=Xsy!d)aQmreHG+bRysie>C}T
zhfJz&!r_ZuRWOi${#qx!Hi*+jbOxD7VV~1#3vVE3+{nywC(Sdbs2pR#<5W+$3+<yK
zDvmV@I&b)$ubvRcoE?~bWib27U%L3|_jU$_@@(2~Ur-D-<0gzOY*drZ_bLRF3B!(4
z9m{i>#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^Cr<S7h1N@J|PzPx+tZFFDLsJc#wq67EB1+E@ma2CH?Aicx
zoW|o#cv|C8q0c60YdDB>QcGqrpsm5eT{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=<pKIxAu9x@dQn~Uy%63|yQ9FH-
z=XZ0TnI|oMua8o*AhfcZKofxmfbREFs%oBf6Ezx^f%>H!Y}xDcehaNrd(^3-G^HyU
zaVxqP2lvga26<XdRIF--iX%zNWUNU=u%R5pVGh_m>go7kaNLT8i72<Sjra(I6IwA5
z8nE@sAPZ*tU>Fnedgpcar794fjj(2#us4T2hQ>eB4+KhvGfuya7w<o=2zE`FxqD(L
zC#J+7b3321>&?y&1+h(S`<Udz3H3IMGV>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>r0vwd1KfF<QKX<chpk-Ozz^Sje4E7
zI9$4j4@QW<=(V49virnG7eATnBkP*DTsurynt5B#AK9#930r;rI;O;78!*`F2-&+3
zPP|P?qjO?!<c0k5;f2M*cjjI!Q*0F!COwV1tfD1}#L!avx<m3hz6>wW>RhPwyu#q9
z$(L<Q!l`rbpQ&r^irKS??C?##)KBe$&o9)>c~n&O>dP{F0!s|n`kvAyHP4c-`@UR&
z*<K#d-EDS&FR-E(jZ!umXnW7j^$#plYUy&yeTH6!?j_OJi-+{jx=WI4bX*`ldH#Zk
zDXjcy*VuUf*tqSVJMH*8{Kl^xPiiP_ExMP(>@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<tMhwrO10;g4`tb`H)vA$t?Qq>*`$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.