From 427d4f9cba23493bdd997685635dc72f5c96d0a2 Mon Sep 17 00:00:00 2001 From: Dane Sabo Date: Mon, 10 Mar 2025 13:31:58 -0400 Subject: [PATCH] nuce 2113 lab 6 --- NUCE_2113/lab6/exercise6_4.png | Bin 0 -> 14639 bytes NUCE_2113/lab6/quick_maths.py | 61 +++++++++++++++++++++++++++++++++ 2 files changed, 61 insertions(+) create mode 100644 NUCE_2113/lab6/exercise6_4.png create mode 100644 NUCE_2113/lab6/quick_maths.py diff --git a/NUCE_2113/lab6/exercise6_4.png b/NUCE_2113/lab6/exercise6_4.png new file mode 100644 index 0000000000000000000000000000000000000000..e73b5c40f07bec500f882973b686ca09128b567c GIT binary patch literal 14639 zcmeHuXIN9)+HJs&3U(9}EGP&FDj*%~5~@f?2q*#qLTI5kD>g)tCQ?Lt36b6dih_vr z5&ANRMndcJe^Ip^Nzd+xnI#D`?D)?91OIo|P(G1h~N$`>{=u`!`g zsEug3vzJgPMlTeKp=0f8_(Y`kYd8ESZhubOUJYwv?|8-57^QT@-uedC{)XAry$;5< zc4k;BA%PPDLi~HL+uK{)NeBvB{(gV}*49+8no6X?N!D4*>DZxAELV_!4DqsYW+;?c zDEjOvb*IqbHdiPAo}2V>*W=1s5`N3M)ccLA6ZVQE$3`AG z)}}dk?riaveBBMxy!@%D{F&-5HBBu%yQH@7NU~Do=h$&Oq_p_jid1347WBu0B3WDG z!KWy}qa~la=R8vQw2oIxO%_zUDvb4peSUvMC-5K)ZRY{r26p(rIQJaGN)+nohE<#3 zm!sQG!e=EHQ4A>5hBN>Fi+|^lob`5D;b8gY%_$b?si$A9vIU5tIL|#?Z1hf9pea0Ni))U)VSuINLZ45bPbnk>^OLTL( zDW%Pk+<1;vzrli7)>S*!S)!qdt%#Nl6WTeSATI)Q9))SK1HL?@7MMAOpJQXY-2?JV&YvG1?( zO4)Px9i7SW!Ub=0b94J@>6uxdzMdYlnUOXINo8Z>nBM9LZ2kgzLVEb+)2EyK8O5wS zP9Egq8Vh=V_-TOqHDgYl<_ntB&jQ^0_VYdWbmJRZ-w}DnyP>f$COVp{XdHfsh>EW1 z-bpB;=e2v#EHyMW!|lV%x}KGFjkj4@SUhi7^v+067uxGOP@9fsVh!_#kC-QYwz;>` ztUS*y59<5)Rvk#1o}Tt1>@s^tL>)SG=+l7xbY2UouC9)f*zw{z7E5}rB9qXVq!nR1 zGSO3^ap}^x{I*G8maV}qmVgU+j~uU#=f85tP` zgoPams)%u(?42vs7YORD3Uh6uEKLLz`NKNFXSch2%C)hvdHCoNMTBGP)~y-XZo`x| z@}D2#+m8>24!lCiUqYiF>=4k44-MVj)7zW)`t@laA0J+Re*Xq{W2+_lVoRFwqtRnH z9L_H~>VofJ$J^V9Nl9uqZX_R42_MRKfO)gZ>0GUE{dQH~kQ*S->BIY(3_U!Ex|c6! zx=_t?+pLVQU5igllut-V@XI65=BB3Bjvmv|(IFBn%*``!6wR2}*u==lx1w}9ZPbF4 zqMw+OBE#JWFt&+#|7#eas_L_c4N;Ng*qMQWT8R%>^AEdTp%_tKZ)zL3 z+HAbLsSZndS_D$0AG&bt?OQ|EngI5chgu8Hj8b_8zk_>6bO#CNqk zs6QMN=lx%wc3(?)ecDzc-&6H=y_s2k?1;@u=H=r9ZZReh>A_9O%`x-yM!Y<0sZZuP znf9VkYbDFe>--Pn_2$-ZL4NnT!sCTPq3%V(xKKOqAww|S`9JxK{jn%3fJ&p%$Et2_ z1!?K&nRb0uTM?MlePej@kYezw_;?-ww;AzKfH2)W`-JmPb|=ivIw^$+qct=%PD!kQ zyB^!Ld6ji*mL;!%fTo#Qg1Lo-h+VH@Ww`WZSdQ*hGfq-$GHyBY3JL))UNrZGK6%1i z86u*vapT5$D!~(9=sGi)TEh8^5GlwMlH>w(WYk-$B3Vi3wXv zGEtvh)Jh#zA7ReS%uKh)qmXN0yo__yuCk2NQ!m18;(^0xJ3csgC;AeA)YB?B)n2zd zi`$Fgf{hIgJc5GNSSdk4m3+rhnJ!#}bMG$LH0A?!G4k^AmjEo9QuLEjQc}{}j-yZ~ zwb#@5yu6wBT+-4C^7HdsD1Gkd*YNeLTPek3w$s~X`~xdBli;od&mfUNAl1bv;D){& z1I%}KmnrQU|JXGS@Sd58_5-ACY;JZgbsbQscDaekmfinfFbmLi>&5oD3nv+dI<3VP ztS$#OS`67c6=L+AI`1aN$3LwfPGxB zG0LSqwyTWW==c8&;O8}mhld#qgzE0NGbk!5`qvK(;W7y^F)<{KxAn(tyX64)KMvUc zR!4e!O`sm}oH&t1G71zlNb&RALLSb_Isw=W4^3s=x^<=nK}358hflvD95F)mEcr+L zZJgW?UJkHXT3X60DS25SNCCF|q))z{KxLTZ%wtj*)W$5UHeLcLDL#IL^RC;{;*5xG z_eYtE2(AwwJ{%ShNFh1FP-}?KZs$69p7{Zo`O(<3?ZfV(=W?V4^ctv6HoV|@q zdW=F)c)cC0WHAyHxmR30eSJxpnHou2guvEeO;{O;iHUC5)^85^u7%u&nwo*56d(|m zZFw4ii{{s_Kjnfa6%Vz=u?-g_klVe*heKMoBJfB(vL33c*26m7_6s!RuaM{o$<583 zuaBVHXI+Nt-{;vpzU?TiTSjldj$-k;2Cm`e=kc+z7e`mdMggc;GplRL1N`ZxUcFaW zAElwQYjDroYY0BxN=`uMbupJyhUu$TV+dy)+YhI#+;MnM#Sn%7=-S{?f7+{b0IUjEn>7lO>f=G{m5NgpC6X#>QzAM&Z;UN|8;^TQb zNbxnR4ADoP$gNu!N7s7oLB7#JD+k6Zn)Q6*+;r++Ed&pO8FQ|btoGR@bbfq%{8Lp` zQ&pIxiWq6w!Q$<$RUt>O+8UE|kh$2xhB9MLc_>0L)sl9PMdIy_Cg@9WFb z1U(Y!&fWVwUX&)h+No1F1Nn7cqn{s2C@9baTH2ChV-zu2UsspmvFLI{H$(1>H?wqS z9<%qE1Lq%)(co|UM`80j*xSQuJr8`yqYHkt0JiJk-KY9hfhX{ai|g3+epa=#9M-`@ z{ol0i^_jw(YU1|&PxJEf{QdpAMZ^95qkDUMg)ola-alw5a88d^eDNB92yak?5iuDwm(I&_uJFlVQ~x?~=~3 zd`lfoLTbHDvvJTsT7|%5TR5N+!l^kI`!Bdm)n1^vDgHnfrMe13A6N7+eac(h{GG%u z(>vD3$t|s|D~Kl1>u8L@!NELn_01G>0|f;ILeeeM1H8QQ^|4B7?(U20`uGiqdIt|4 zB)0y>s#P(V@GV=mFlMnP&quK$TmP3YlO8@~g7uO3=FP=RGXc4av>6KkFJc+6bM&A`WZB!>Yc#a9 z;=X**RtlHW@yjr88Gx)9u1p7X0)}y*^Y3mB2783{eZ93Yv``HHvI;pSC(V|_)u&4|0k2LUqi%H z6 zWnc}iPp?3|J&7E?a~brp=i}7mRah!LPyy6vCmlFoIRNaNH` zURoprFe(j`B#6a0V$TrgZNa5#{0~en3H(kay8u?1czJGtV z)t8TtuRc*D$;C%YOUt68_@-$_keW@d-YK1*80Iq#S6A0!(qm!hKcS)Qv6pbwAGRF%+W>$6Wh#64;K6eM=^!A8 zMHPn1nVQDw2?Uu`hQt8$dHvj>eopM$CH4HSJHuli9}OU}a4^)XK6b}ky%0&qkzIIk zn{Q`_o72m%L+!4kjn58TKq?FHsHwS}oE$M=gj6=)2M^Y@h}s_g85Xadn}%-xU_x~C z+4Zaf>dMMbZ0+ncR8=EOOV4lLu|s~4cs$#4$-}Za?ZZrPc*L-z%TzowhZJx$i4Ucv z<`doJ4}E>LfZ36cA<77PP=py8QZlTZ&(+mh78XgG1k91_R{Z0~8vwTsi->4knqima zTnhL4x%utXv--a5`~jEtZ{NP1>P&{4b@kB+#OYFIGC|)w>A7usbOrQ-Ud6A0Tl`0x z!$T6<>qnc@FZfDL+-K&LNj|I=i*!w};@35G0ty?dz5DzCGm2fGI`PD+BsO z5cS=9!IcziIONUY>sl!2b`sA$+JW|qReBlAjxA8t(#Hj~yHN0kd$zzaCsqzKb$m%R z%7n%Y#7s2Xi5tPNy>H{h;L9)DG8p^Eou+__XtV0uKTFB&WdrS?{}Ugffxbw)jgvzT z+LRIJ72z>&hOUc?@JN(hnyV)I&Q11Fif^u%@3_r!3wIrk3O!rAzGI$B$-L2ge7G%N ztE4sJ?#vdXf{-k4^V>tU3xs#m}Ck@0dU&V$R&#Q2T+^AQ>hMxojRKW&PAX z*vHmvj2IN5U_u^!Nkb!mD*=yJ!WJ!NAre!ro&ZuEPTpUYjKbg=HI{`zQ0{0xENsGZ z^>Hesse$?xI38q3cfAUkWhtTs-^4z)t0s^3hUU|igw7}`x-Y|Ztt&QAw z=ANA2aUmgUr#Gj(U#v<*6&3}IHOEo{Z@Q}^!EGg+B(<(=vbq`_cKkTM|TGR-cY*14<%*}nFFU&D~^+vX6kWkfS zY0cHQoDfT*bKicZjx5VoVB5)cb(b2On#@5&4v}`vVq=FwGCJ-g>^5s15Euv|ZhEyi zR6I3MKQQi8jm5cY&n8k49p}DCF(6G3w|*BkptiN;;dps?YWxl=2%FV$fmq82zzV{~ zke!FFZb%csfB<{YtiLAGy($;>t1z}*3MpY-r9So0x7wE<47v_F+gqvauq6<|#6*Vc%!t?Sk~IFp*>8_p@*N2tW4Ad&t^cT|a^V~t692KB%Cbiy zIweH_R^FTs_rZe+xw&)6O>|JrV{4->RDjGeU~i7d-Xa!F>-vHM0?Zl{H400Ul9EKM zJBrg=DjwwMrMyaR3lQ%le0T*PG1uXTycpf0$Tzv&VWp@&X|j_Buz6kxpRGjF1N`o<-F!1`Phua1TfHaxRLC zi+`%D)HrwUE-W1lZEcG!)QZ&9)DS7>wCItDzKZhl+8mp%27s2tL|&+8 zJ$DHFusHh#K~U)Zg-+og62hB6^NH~FR4GnQ<|mV2Wwhksk}j#KHP+PZMb85UXaIqZ zjXkt+=h2f}xs@BXsz^mHEG*bLU|`&lpygfk1I-p#c0-^0OHu5Ih>1=Un&jGXbInxo z3YL{uE#R=t2ZuaY1Myp?q5seMSZ5U(zP`C~yY%s{D9}QUO-(KOs=~NKDGvFWhHn{_ zf{&hpl?7_(h)5rzY(fu?qtEqG7AC7n8d~}>Bm+lHgUi#TD4{=0o_5ma!X6R zNmXb96rm(o1J!dJz86ZoS+oFEprp`3CnqWdDS?-*Y5}LBwG7-fkV-oN8{Tt?ZWs=A z@#M#Ie}Vg6kFfKgPeyYd!1lc0Rh&|&aU{09XLzZr}}I21{0H%#o~UDJ6GQjjOh7? zzRk3?*lv&PBNc-p4zctuS^Alx*$;wK@aXttWaxJ$e=Tq_G%NU4RVf1}K%{Bd_ze2_ zMf-k`YBj2boDQx3MYI16s(n*SrbVVv8MD4^`9bb5_J0*aJTxt2r)tjVrknVLd(xeI zy1PRgbz!MI+doA3?+l;Ch9TftnZ^}6$VSj`aG?ajui<;4g%gL)>!_@(0UnRnEp*8Q>-6~XqNaDSFtKlpDHh`y>R(a*4Hhfws4U02m&W1^Z^~C(0b{(U#88N?=zASBk zE29`6G2XvV`h(Vq^7?p^;q2L@WPcDjGp~Nw05V@v+VKWSM%3N~Rn+EGu!kf^o))x^ zh_0M**^@V^Y;zDqO0`$0{Z|`*XnapbBcJ+k3ZC0D!~JCd&UhxB6NTCj9qwkn-{sB) zo4?7OD3tPJ7>(vHxCg4!7rqJlz4b-i_>WWh|Kpa?U*Jm6DE0Fh%sm7*wA)%t9fRR0 zs08YJ=NU{gm_KSZHmP9oXk5OW^g_tge0r$a{KkzTyFa+!pFe+=>C4E@MtG{r*xNNC zSKr?^w1TsyWfJrI?p^;ZYu_T6Kn(pG#c&E1wYKKtYqR`5C70~pAwJ#hsxlWl1}T{wk-Z6 zlJ=V;zt#a=lK*AsdAoy^mDE?gz7%ZyVRrFgiy6@$yyvw2w6#u~<-qg>d(Ty$(E%E8 z8F!ca4gl|)a+m!MXz{OizPtOWyj<1FDutJqx8ci|_;>G4quW#}N61A6m$Pn+OUD8; z^RovjcaN0QMUWgH`uH?okM=l1W4AnN4O8HNo)KPiH!O>5z3^`hI2?=#I zwJ7-PJ1wF;Cb?4^wOsMe@H3(@1eWWwy%G)VVgbJZlrG%MmoHUZmHsU=!;qk(sTr%0 zq)BBXq@*-~lu~Uj$o=gf(xmV4mgCFu7N+)royV@eHzVZ!TJq~EoQg?LSNk`-3b24H zD=Pykdg}BY+j`i+X>S63MP)YlKdknF#D}hf_7pePaSGt*@QD+e&$(3sxsF`UQe2uR zEG2R-bgTghGR0#ipKMn4WDoE*lDMoaz=qvFh`-Fu&aC|W>C=R}K1@3PayBpMf*_vy zMGw@+QOda^8VYCHleY6|UI+2CCCl=q3US=iG}5 z_uKWD^UyHZTCT47V6n{V+*RBA6Bk{|9@gU@B=C5dc7A+zz%Lpcv`wQxt3na6ICM<( z9@UY#qI=#oumD3qA_UM~hf-a)a1i>Am_WR&%d6;v;*U1F5Z@p3ZC92QYH!`*^D^BB zW>$H5`6=o1r%y9_qQCr*n&RVABP>vlrcWbM%yRGsQZI(VkYTF-6u?JYs5~(CTlnku z(xhQMz8~oxORRr(55DIgk1F~H1eAAL$LHjH5jD{X94T7zK$OcBt2fjjc{j&)XKwqK z8zO0=MFz~?C8I{~?{6m3%0H2T2jeLC?7>^>1!`bB;09uX5@Taa9c7(yP$4MPzE7Rj z*CF=d_i-To$sWmAh!JQ(0;blBVTG{6kTxDusEZ+xM7;~>1MpyxP>Z>F;_3}MUm?@% zBicUJhPBo71lJX}NF`Y`#H)gVF+Ys;q&OCZ-H_t&B~tO&!^e)P-d?@2E~se17&`E< zlJx#H{CB=m2t`(4%<2T`**-w@P4et*7Mbgy0luHF#CjMZdiba6>TC$vm}O>W-nB~K zmuA}|PxONu>%mRoLp+x!0q&WTN?NE)CQ2uq)*B&37yT~7=~Q>#fum2?H4rZBSe}NWxc`dyLF@7gsnhG zUWizI6|107{pxK}@<${n)2-MU_Wbcuc09m;V;VnQ&s`IC54at8(+Kw_tdx*YJ+#;d zld8ptsFw0xh%;#_Dn6!mKUbOW9U>QzLGnn0j%&^?Yeue2t$1IS{lFzaDjd>^E|ki` z(}IK|A7qF`fpW!F2o~ufIU-W)hV@h`l}*IlA5jAkF}s9C*8i+#nQQM0Q~&4Bqh4Z( zI;P$a^MR45&xF%)>hOGbzrr#n=ZZj{sksN8(jpo;|Rwy7(sM*=2gLswF)TEgvQRU6) z5f1@V^M%>*7NUDTOe__a_GNhD&)!Gix6Z*YNuWtsjDPz8`^O0aF4_#Z;kd z{`Be7CFlc?Pte^|9Z=;FMhT(FD2QFg1D{rr8oz@qScn#;NvxA>7KUdE2;CRqNnh93 z(Et>lK@G{N>cR{|R?&5lZ#?~Xh46(TbQMHG9E1YCe$|6fDlbtDmSQRd@)ft}P{&*f zfx*@SoF7DkOXC8fZiY?Qc?c_f6!k4wfJ7DxLa#lvdgmI^&H zqVuXoogc_*wVJFvK(IoTV0+qdwgm_hwi00>61c=2{MxS|u3Inaovb{hn6|H|zGQpx zCP)FPDkwYxe!U&sv;O|;g%h9nkTe1OI!q%~k1vHlU~PaZcwp;Buol!}6;<8c3r8D^ zuQghBe|&~$edI6E=;uJYdC=XB-bD*Nf-R{=-Vnj~#(sd?le{?Mh(nYS`obJI{ST4I zC*+FuY9DY7Qumb25+tm_m|lK|yhPxWZVoe?Ux zSi7mgsvA0Of_FJ3?Bl7R*Xy`W4=pb+qxM4`PUN{N$K)UbkLLNY5*EDd;($t1sJLA+ z)H-db%z?(w8chqq+tjqT&!8@jdXgcmmyL*Cmtj3-AoEY&Jux2#8aL(jBY`wPnGwv) zvwGP!hgQob2qeOw5syIU!Gi})GsGF`<)dM-w!459mJIDl4MZQZK%B!@yN}pu!ZSla z5(f$Kq^zta$lqjg=|G+)9U_rn=4}9ys42_R?b4;x^S-dWm6Dw>R|ioUSU;K z$;a_kH595Ys3;jQ5m!JK)Im6lITQwNd-e^H;WIZuB!)Dy()aBD8Ia<(_ACbJU zNS`S}q7Tq>IA4BB1c^6PNKCgV^wj`Ns29lZw8ynlySG($gM+OV81CLye6}Ri{2ZxjB456lr zAS|qPWRrq0He7b|f%DhFli4SJ!fWgTT-RE<#}r@z!2~*aiRL-qHtG=v28=QY=F^q5 zA&o)^ow`rHBsTzs*v`pG$jA@|!w~HU>IoMbl9A6Ljnqa5(^J@S1jz;=ZF9&8%Me%< zo!LQfG+sYnGn?#@kGO@81oGs-R`#$92?-JQq`PfoKcNWh;u84Y2$+%=ipn z(THUN8PN&8GzVc};g6lx{lGzD5LuaF#p>%V8D=i9VX39;SROFGARgHyLeh(Z|K~Zf zoyM=WS(j|o&3C{e3r8$A z$fFDb^Qv*5Y3rK_2K-k=&d<3tdjV*@@(kIt2%NsC$!b{x#IS`T2%b&6uPlA`@$|~# zlRHo=Yb>+dOd!~rT}_*qoWw&=Zq7;}EgZbFa!}i*T?(7-vSO}Ve$A;(!;S!pvX^G+~EK zWoa*vpzPL>NWun$wJG>fpGPdb7oCP>`_AuLApkBhH7)J9#k)J}JS-4e2wZ`fx^C+s zZ6VXDsWushpoVzq#3 zjl>KQ2fy-#&;`hRaO;l((F1QUX=*6FdJs6Z1e}r%m@|{d=md`BEYU%1Hv$-7tZ|{C zk^m2=dvm2^Ug|v92JKE=8;kfnU{LdAvh>L5X@`xw#cG;zXEKA}y$>`5=+zcVU8bA5 zfaYqdtNW8f!Nn(A#?>dR7sz=<>Oys#g_Kg&i}rBReyu zaq!6cxn<`bjEa@L1O%BLX-k>{b6%Gf^={XyO(%<;r8t1Kfsd-NG~>Ax2L%TMVYm<> z(M@D$ji z^58lZ0+rHC*0yQUdW%pr1Y$(LTHl7*5{5`PVsJw@14d;d?29-6LmsfEk++kOYxK!3 zA*Hu=1&fq*mMX(ZxiMZt6vGU7S{A@zK>4jlfCNT__LZ5BN1E$$czb3p7(58d8rj*K zz%zs`TZBb~h2e#gPc!LQ*@aFPaC7|q{9c1miEfJ-NBl9+1dz6hf zvdB(}w;&X10(+{3t=r$)*QjNHaC=P=(3Y^ WmYB3H&j|3BDD*kyvni(yZ~ZSP_T@SN literal 0 HcmV?d00001 diff --git a/NUCE_2113/lab6/quick_maths.py b/NUCE_2113/lab6/quick_maths.py new file mode 100644 index 0000000..ae499a2 --- /dev/null +++ b/NUCE_2113/lab6/quick_maths.py @@ -0,0 +1,61 @@ +print("Experiment 6.1\n") +measured_activity = 5.2e-6 # curies +elapsed_time = 193.92 # s +sigma_u1 = 35331 # counts +background_time = 193.94 # s +sigma_b = 558 # counts + +R = sigma_u1 / elapsed_time - sigma_b / background_time +print(f"R: {R:.3e} counts / s") +print(f"R: {R/3.7e10*1e6:.3e} microcurie") +eps_ip = 0.101 # from figure 6.2 +f = 0.6617 + +d = 100 # mm +r = 38 / 2 # mm + +G = r**2 / 4 / d**2 + +A = R / eps_ip / G / f + +print(f"G: {G:.3e}") +print(f"A: {A:.3e} decays / s") +print(f"A: {A/3.7e10*1e6:.3e} microcurie") + +A = 4.5 * 3.7e10 / 1e6 +eps_ip_est = R / A / G / f +print(f"eps_ips: {eps_ip_est:.3e}") + +sigma_u1 = [10122, 11763, 14464, 17476, 22073, 28268, 38084] +sigma_b = 130 +d = [100, 90, 80, 70, 60, 50, 40] +t = 60 # s + +import numpy as np + +n = len(sigma_u1) +R = np.empty(n) +G = np.empty(n) +eps_ip = np.empty(n) + +for i, value in enumerate(sigma_u1): + print(i, value) + R[i] = (value - sigma_b) / t + G[i] = r**2 / 4 / d[i] ** 2 + eps_ip[i] = R[i] / A / G[i] / f + +print("R") +print(R.transpose()) +print("G") +print(G.transpose()) +print("eps_ip") +print(eps_ip.transpose()) + +import matplotlib.pyplot as plt + +plt.plot(d, eps_ip, ".b") +plt.xlabel("Distance [mm]") +plt.ylabel("$\epsilon_{ip}$") +plt.grid("both") +plt.savefig("exercise6_4.png") +plt.show()