From d11e19a7d88ae92ac57ec9f5ee1e115685e4b372 Mon Sep 17 00:00:00 2001 From: Kamil Sedlak Date: Mon, 31 Jan 2011 10:40:44 +0000 Subject: [PATCH] 31.1.2011 Kamil Sedlak - changes for Geant 4.9.4 Changes in musrPhysicsList.cc that were needed for new version of Geant = Geant4.9.4. The changes remove G4MultipleScattering and MultipleAndCoulombScattering, because multiple scattering has to be treated differently for different particles in Geant4.9.4 and "G4MultipleScattering" is not supported anymore. The documentation has to be updated (i.e. comments about MultipleAndCoulombScattering has to be removed from doc/musrSim.tex). --- doc/musrSim.pdf | Bin 210839 -> 214978 bytes doc/musrSim.tex | 20 +++++++++++- src/musrPhysicsList.cc | 67 ++++++++++++++++++++------------------ src/musrScintSD.cc | 3 +- src/musrSteppingAction.cc | 2 +- 5 files changed, 57 insertions(+), 35 deletions(-) diff --git a/doc/musrSim.pdf b/doc/musrSim.pdf index 572dda6d10820d8542ceff50b3f985d7cfbc1ab9..b432c203cf797c9fc05f424066bac49f6db8bdb7 100644 GIT binary patch delta 48927 zcmV(pK=8kp?+n6^4S=)(v@(A@%a54lse4!bsqVuY`v68@Z2$%{X2W6pFi4=05gGv! znB70WnRRB>mtXgc&|t!WHPPL5W#xU$5yuo=3Sv8RNNkdXet@;@f|-zwW->E4OH^ z?C#ZGE+$xEc0cS7dJ)D1wfla5uuD)vgxzQD`=@)cc%!Vs(Th?}o4;NBH*SuuYLr=w zjO{{pMX=%U^|guhI1BJLjG9#z1NeYm%a zh(X8Qd+@ch%In>q_6L6hdlXLY-r2{nSg+jf0)MDwnZ*RZ`yzr*qF9tx@I$;v z;o2EsFx}cg1(-Xy-J`yt>%NVQBi_2*3;3y8to4R$<2hZjpGU!z&Urfhu)2zB5defs z{2p!~Sh7z7q&|r3E+5#MST4XsFX@(_%7zuy)%C3YCWiLYba{Q?MYNz<^Gv=?$dF^osb2e*sKU18myBIcH$3JfP zwbaj8($!WhPV4gnk}$HG1ep7HTIm22v{(Y(>MNMV3h8Y^a&r7)Z2=gK0q=}`0Eg!T8#4wI7foOjXM^^CXOE?oiDTdug_6aw0k!ur`4BdEkEEb2ppsGi^l5b~m81)i zUy>tT?4^L)Nk~vP|0RF{_8v3Vt;grfGm>5&PzErvP_g+Vyl7{vw{R8%tQPSe)>@n;034Hxi%;KLwtpO_YDgO5pIoD1Et!XlcZ?v zmvzqvtwK5wqN5~qkE z4`hGD&;N3-fI$NH{Tb#IaAfO=j3f;5%eGLL^Dw1Gl?sTsbIr1jRbls5`|$bvwMt`P zxL#aTFnbk(xDVKp28ev|MLWm+LB_(8aBpiZP7g;TV~$8`Er15P+`gpJqfY-gi$K&= zysJ!O9@B3V=(fWeYsWY1NvuvCtbpK=nB{IUb;VjqA3 z+n|nJ@-+Hw)I|)kRk99pdJus~DGk>h)y?`we&mF312e+-P=S1q5g|&-ztynah0oA+G$PX|IkQr|IRTd~(>LO;ek}9oaVT`6y=F#2(Z?iFu9Gek0DoDhg za|btdm*NIhVbE{RosWyVg$VXa0CV)|y2|%g{m-A1G_`Q?9Zp-ihLK6G-t|GXu^?Cd zaL^x{esI0g1Prk5kKx3CSB?Qgm%4xOJ(|Kk0R4kWbQB<<3Z4WA(5*#4+%-T9D2g$< zqX0>2#u0!JHN=_<5R^Kq1LO#yBz-Is3hai?$q>r!fa?gCH8kMH9=Y+vO@Z}WKwVG{ z3B+IQ9Y8X6&9P!7F!S}Ef-$Q-Z4vM8eb|!)PF7h|G}>1B0kL4Hw%OJMJKlc;fgI`B zK00hpGZCuhU+2@p*N^*{{!6^gh=M4*o^7C_$c)K>y{_(%F-4sj6Wo}0H)c#n^k4-W z4?1=|_&szaSs65yva@2RgQ6M_NZz9Z zs;x?fC!YvFBvN@#G(R#aC^COW0;>b$`eN6#-49?^z^)*7rfV-B;wY_wc&Hloq^W#F zrF+yh->wpNh2ve6oS^VZO%!yZx~+~sDUaPL`D#}B@oUa?MRr8yyoVqSiyv&8uLeLF zxF|T)3go!vmsuMI0`k&#NzY7lJ2w=V^@$e}wSNV&^xM3Z7|o!f*S~*=I!mB~oDP^E z;9gd-MewtV4dy)PI)S>WSk1Zr^>)Thz>djxW~4`6`e={#sSB<+8>j8PXy@oYG*+w6 znsEVmqs#)QD6B!}F7QZ`L~C8C>&IOae$;7TvM14lEkiXx$c726_8F;96V!hIERbCt zOK7ry)pTSqVlYKHu(E$&zm4{(15QOw18X}E+9=ose{CPywza>t&%W3CLmU?1DY5ELTB+5~7I;sTFSb+$k&GjHnooR7C}77qCwd6>!nE ziHdfh?AQamq}Ir+>kpzri%R>?C*)4J(w!aXCC~v`GsUH^5$Av4D(t3DJ`Y;E4v?Gu zl!RPBQ<7NE=sVDHF|X~7(4bzZ$a1||&=5{LyBNrtOWPc4WxBXX47(ga7%9NK3{PU+ zUF}hhqN(&VEb4-f)Xv;SDkU;fiRkg)J3}g^0p%mE18lmg$3~!YwOgyi^!zvlp#UP( z{k+CF^jXU?8l-=fwGjlfe*=I~!ezzmjsc?pAfU}32aHtermzXw%TZn8s(IHL!o*08 z?V*~PIw#i`K~{vH1AHA66qIHb^e-c~=+t-t+ww&6TqSRwH;E!)Apw?PC93WbBJ!vz zjRve-UlW#KRUF5k zzQiRRT$m;wR=BFsyVz#BW;K=-51TxF^DA`~szhDT#r7&;{Y7^lz2mnB25PfbEP5{qE z!(v_Mxf*{4MbH(LB(i`*)&AG@USgcv45Evo>KN3b^Ah>;jQy0DTYZ1273^@S{>u%}0k`*fm^%iY@p3@{ zOF>861b=Rk!;4PpXmdx5i7KLp=2J@VG?74=Q%L1N3|66ak7gnO^Yb92O=s}`Xi|0q zFrZ}mgi~JJC160IK|>IiX@DVm^JV6N8a=%LI;ny@5ru(kbUx}yKYmlP<^4S>E;0yD zV<&&aC?cPd;LcV10Vi{kMsL?{k=J!uVSv@Jwo7MEeU+;cQVmRGxJMHC$nfND;xe~J zA75GXUhy7i#Tri831UO2oAWJ)XUwkp9V&nrc2jk8VMyZYoOg9{WJy8X3rkC(xDk z;S%I}b(G@YI*sgL9~&?(X2B%9k5n+7_E6xxAo6Hem3wmOIrq!Mfl|h)ias3iN|=9? z6=aq1ZzgEB<(2R#HgFxVNZyF63*>!t40ay(wz6nE`CT!~k~qD)>ukzK>FoGeY6;19 zNGO={0^P473^0Yb*-JM}9tv)8BHY47fnMMU_zXG`q0TFYLPoI@GzPy0K{@%i8{;2$ z`M3`*;VT~#Qv*)AU;@ZxiH4xRhlziN%qe5Yut{#r^CkoaItFNBz*bLbHgIwksE=ag zbGRd;S3)?(qj4hYy(EDm)W?@m?UZH_dJZN(C{F*LC{8lsmL6F66YNDh0|r(d@jhIz zHW(x+Y{ya%8pqox^oUv|B$3RYw3ls_bD(84dQ5_~%vz_5}HEV>WQHEr-dkRu>4v{rTR za2KXy>6hReR)7^A80_{O&E$V6(SHzy)n6)PDl z`YrjhGCgmPEBta@d4j{KCZ8%#Sj}xV4Ca*GuH|eoa~98wY0f4~jM?LqIi$te1d1dT z0P|gg?ohx&_-XegvB7_r3N=mhVc0k(u&Hd!nv^f!3dXXaH^b4Jbn40UpLA-@nMJ6| zm0$xN?u0o4|7K$r&gz;m8)2iswLwIc+}4zjTcT`O7dI;UsH!R}TYR#o(011jwOjFG zgG`A_FE}$ISkY;>o3cwY`AB-;9}|NWUig~ln`H%vW7v%)X_Pr3&7}kMsd9}9 zm7aAU8U}}LnDx@n8jq^RjUH0KSeDRvSWP9+V=9bH;Y5$d2j*=!qu8pB=QIFX zl_#=zq1zZn#)f}`f^vi-AuCY0rs2I2ggZSeHosFFY8uN{ay5kk^@6j0PEY8;IHeO7YZkpQ=*5v#i_uDr9J}J?YbvbhnJ-VrZgg!jG`1gdepeq44#HG91Dx}>ecUFc}A`W6TcG^K*EFv@0(O@Gyya_rY)DY1V^@niwW5v1p;g{{W8WK$#x z#~lR_tqQH*0-HQwhBk`&|Lf*L|LG7J$5Q7Q8EIQ|!78hUP(LWoF#TX;3^HacS`CsW z?+{!a#ADj0op*AVSUI9Y2`(He*n0w92Rt1+}!s~d4gdx?Mn ziUK8}XQO|Ti7u(;pVXH?*@fAcTAY*h$=$oMQSmgp5`pAy<#8!|fz@*x;Piwks3k{t z5>rTVH6Mf}sPqyZeqAY0TSdz%MogEa36B5lie-`kRp)3hm0K&hmK6ehFP305qcrZy24{!1Beo#%HsF$+5qupXC} z*{6tEa_o|CwLtGG%@n`xWB#*n)-WxU1c?~@HH39_Cl~o}o?(CIC85yWYEjaukY!mwS^YHGjX`0y8b%FFDoZT7 zOB?z5-UU>CaEtSINvvvz6?fmJZG}?92oj0QoUj3d-HqjqsbWT>RwKK>JmE6EjgKN* zD!zFplJUg5P>TtnXU!fVG-gD-Val({eyUIkY{v*f6B56Wb2lphMk#o5LphMJhm0(ra-_Qsyf~E*KL7bEdLs0e9PM(8~=XJ1& z!WSVNj01eDWP;%*n=D7 ztS!$x3Ha5TK{;vN<%fMjo^5Gvbdhb{m1xXIE8~DEJfP4LCe9W~mlznsvmzR{{*{*mhPj zu-H(=DF9Q>D0+dX*6SqQQsS>27ECK@&)+QF0q}^UZ4*H$3i0OLU-&R{{}3AYk8Gh` z!04iMg!?BCJ*G@5dn;9nhI=r=3AlC>r3!Ng?W|LI1KO$`C?S8hUGP)p(K@$}%_RjL zz+wR)%t>Oy7c;U{wFsBzyfLKIWAs?Cb%A4jc}j-!xAQ;am}n_8IvL1Fg%!bsX>;Bs z)p+z&`Ha^w_e+Q=ZEznu54*{L1CTB!W9wRC(%OsQ0G_d*YIt;BnJGI-Yg6 z%24`vuoA#y7=&L3h;aAX|C)Zglq^NGaB}f0v6Ul54aX<0f$)SsIRlC8x2F+OA8JN% z>@=iVyFHhoX%h`)4pXRTB~$T33qnh19{&e@JCT19#x%FmW|3&!+BVc(Awb^s zeFnEkaj>3eNT_qI8Ph37if6^0}z)}h=*B00>9ee2J$W&6uabzDwj zk0Z}cbF2qu`&e$Izc(X-i1@VvKN`m`_FU|mEw-F8Wr6h%b5fg>>2uk77|IRmeaz^m zW6rb|m@yV)PZ?yZq1V?lPgjqm^K`XU>=XZPYswS z@s)b*3>x)GbC^<1-Vq6F-pf8mEsm-IYdtl=9q#fQuJt7O^Tjk%ZcSYytxsCWwzT5O zS8>X#21C6@zsxmckGR-WXD zRJ_W#rsmRAqN;an3gh5&OXnEcn-3Z*_UuxDPb_NSb~`Tnl55D-@?BZhM>Y;>fWv+v zaPs%}ShM0on0A&f0C$mvp+<*DB~=p6V(VL8aDFfX!h5B>%&mN970<9+m!t*KbEgSJ zFUfx-o)!xg+e>2C)(sI^=jXN6m1c7*`*|j}URg&(WCmv=H)2{Mnqw+|=QX!AxKgJ& zIc2b6Au9SUe>##El&E;nukn^-6Pkk7oXvovY2*(~;m!F5OfrDE$UO4pg2C_C(}uKeiY2x=oRGn?@I!{w{7Rdg*`k zca?2hN)g#gj-m20)d~XtBQBeMo?fd*reAWKX}gicJt zj=lxWw9jpudZjw7nsCmT&5Br7r^FLupH=`UsT^WcEqbd)zI0~Nev@zoPklrlKuJ8KNoEW&p%|1LeZ5BlHx=qXwQ6H&o4|Vwhtia;RPgnfv-Fs? zR!YBb5XQ_NYnl{XD=LNp*Xl@RAN?U@N5DX)WVmXfo(|mTYVWoheX$5S+f}|Xt*6V7X_emOCZ?l$Tnr%*z zlR$nbR2lGLj&<%KC)cMsm|1`IXHy)<_q9sB-KaLd28J&)hUZg^s$S3j7;c%clN4B$ ze80TWOONK9HJk3pA~)DXEG&2l8>ug;9Or2(+j-JHn>A@;xH2-#wbnh2n(dj$#}!sd zFr#O3L~)E2`=PQQv+V%pBWa&-wRGmr`(N(iLJqqPZIzB>2%#u@9}W)uL?SJul+ z^7*>?{7U^dLR7D$?kcGV`*d5$Wkmcc&Jv9sPpNORODMJz!4#HDwtIUw`U8+$Im8yTSF@l)u)0kFZg-}hHw_y_v(vpMMnKk|Pp3E9>^V6K1^Z0WRg zvhVkp#qRVl?3?L{6b~0iADY2wea)a{r+vyFvhvdH-Qu*?(`oytZ|3XJ3G@0he`#79 zU4KGm`bF-CY+`yOrE$cv7`6b&;J`$Rye|v{m+=9kyqCAd^WU z&)){#f~Tq^030*5)k@>OkzbD}QJE}s>+hsfa_?B!w8}<8awG~0ije34ro&4qd}b5x zs%DcH*x*kkLqq%6TxWvwdb9C(f-9)mtMc&4LY&LPkrE6piMhh3Iz11umh>);^`l~+ zC>tgxn)`oeoCV$&{|LTTMj6^2heSquix@Uzp>##Id>&Qz?)68KSVgNN+1n*_1W@c$ zS({mpuMnL$SZ}MzhF1v+ClKfu2an)CS`<1z$qC|O1MT~BLf^D8r4g2j%$u5H=ECZF z@$0wy*LZdT`bT?L2$%efq8?9a6@*MO#1XbWz14r2kpuQb&suR%l!v z_o_40m##4Fh&0$-ovK?becR>f=?$qk>>H7KmC}IIKj__#`jv5hZi1qu#pPBTsS~9rjzU7#EE)PkIhVF9TSDryhTpj1_h5sxa!FcC`hyHwC06K5jt;bn{l6#*b-c*C!^EWS zFJN!>gBO>6}7J}m23A0r!>BjKA0**nysEmpSv|)4^u6Z@QE)m?_DX+y6)QL~0{frt2tpSNbNd*h7ESE5o%W zgfJHLwzs(k%CI&Wkc_wjWC@>{xbbO!7xC>bx#xBeWznrxdc3UwoFNm@V9m?P1f0pLUEy<0<^*ow8DFiAhD!U6(FXz*&&30?e*{p zAbGDm+9CH9N7gL_9Z8oZHiXj!$f>kXuye=FKYEWFRgmG$hv%YX=?wwG$=Zco0T9jgoXje2>PcOR1}JyYG%e zWmJ0}2VUa$sHerXfWz^(u!dLs^@xlk-dU%xFVB?ujEBPBVf0~^dVySn{LMUxEGF#h zIds!4Yz5Klh~A!b%(Z`TN-UB;)n)lpPg#Vy&$y!1qVr)_Jp$`Cj&}Zo8@9HmH@ve? z_KQoACU4M^lf)vs>Q*iQ0T>Ce@ zyHsym*I#*Lg5f}|t>U&SdHAmTQ{w(eN5{=Bl?HH`3~%oueR}r)ux>}IWo~41baG{3 zZ3<;>WN%_>3OSb@RRIYDH#IevDOCYSe;l`UerNn6Sb3VM(zY6R;L$GGj#4E>k?7gfj1z`144@A7>VgBEf*AIX92yP1x!xpTx`r)JR9;C;`LnpRffA8|4 zvuYFe;iK0N{_7x)S7!5JH~JvspFRZj-;e$k#!!w9qJ`hI3NX;4r|{fY$1A(-yLbNZ z%~5RHn_VA1KVI#d)VnYy4IvF3;z|X-8Ts%^iZsYyk9yy{@`pFaD_qjy6)wR!W#r-g zao@}?l{JUgcvV`b{o&d1N^W77e>NO`I9|a_FrKTA-?V4QNFR=3$2a$fPv8;hHW?JX zwLJYJ{3>v{>hQ8XeS_2ZEf@zkz;rqUe|XWJeuiV}P4DLLN4V;|a55Zp8dFGl_!@55 zoC+$G)eCJ8&*KfoJ69$%p}a`EiVK9ff-vz-H?d2p4_m!AZE%0iK5W|Ae{w{BI;u^= zV;TogumQ^A3dAO*^y=^qE?T$kLBnRj!nK6AeHs7K!3o%P33obU^F{pmM;KCT7t)Z< z!fV^Zb)p;Z$UN;YGUD?@YXV;BrqT>q(NlOx!RskUCVG>WDI8hln>yq3`d3@zlQ!aY zn-K^?X&^tLf$XG!V^AUwf0Td6=l9(K<(0L?El~R|aFxnRc-q1$%g?wpV>TG^hF&Wp z{)o>6q}IdR`un^1y8sp}w_2c_^m~BUB`bd3e!r%p;?4jWmYv&`M_RZlr++16I{$R1 zb@4T|0ovU?S^afock;&1+r<}0X*T1uAOqzo@zmd?F||?Ce#Xb!e`yW-DXr3Uphhal zo_^m#Fo&N|l4viCGO&>Hd^_V;jjE6VjfYt;Ao#dqrS3@Z+QQ{~CU~uZwM|d&>8wJg z6)9nh$wM5;H-@;$L+YOBkV2%>+up=Od*K>R1PaO1W zoYLw|!Am{{wE?oni$L$4#9@ePjZCUsfQ|6IAcz;Wp(`DwSG?FBdcl$IyEKyU=8j~F zV@eaEMzdy9+EMs$_zuQLSpernUcU{{sO1ITskU$r`T_Xce;6^#Gw)7_GZYx89)KuJ zz?qzKF`L3iNsk@dEhxMp@D@4k_~PLob*I3A>JBAQ-S z+J8gk73j9G`JJBs^$|dz+#O{`q8Z!Hr1wzX(KQv4ag=S_SLGWBg7cTu3W0GJf1Qf^*bIz@xkWyCnC4qj2jfBbLwBd};N49-XJ`~haPu)*%| zd%RV8p|XXsTTe&2IV!tZ=}ns22r1}#{u7@z&L~IifV4{Q3hAnbcZL`#Wd62N4&#rc{yV%|Nd?1wfIpGJ0@c4nV|)h_Yl(&ktEOHJs0Y~x zAer>QfAh}pTBbcyCU<6~I;8WkHcM9&dBVHk>BXC--gcFhN6>JPJ zZdz~Kwy150j{5ePHx}&r7aWR2Qm`I<4+R;xVCZDzPux83U@W~*gHPI@{BqcLtQWHi zN;_WR;C1t<6hP;io}{`gDO3<8y+SRVfbvqi%~``^xCNs{DB?Pe)y}j%(_nx4RYfqL ze{I%q_!=$(0u2?5Xwl)XbfX{y$OA!w06t0bA4D&MY4Kh{Cym%gSW)zL-pw}u3(9?I z69gSSOB>%=g<0AN-P?KvdlSLJS%n3BD%qpx zh#~%}(45m|AhdB76QSjQQk zO-L|7y|it&eUlJv%o%%A!Cca?= z47JnBHBd-}ql&F65&0Eg)0-&6g)9^gmIGzIN$6y^Mp?!u+3wW#KyBKI!ptqqtCil2 z9$Nsr?P~UCq%DlV(gl}x^6hALe@BJ3_IULiHlA(swI#VHK`W!f1SNtBVB(B!3JrsY z`K2`qKHL=LH{C&14rbvq2B@Gxra^Wz_G#v4C)diDsEXlq-9GxJ(E+Sf2Mo#2GuaWF zj+%*Wx8T*_XB+OUxg^@?I-z~FTiLH8ZmjG*4f)sJpHVkBp-Kn4Q@M zeC2h)*L9BF!S`hU%U>UkUV^2R{Ke_vj^VbdtB8+dL+`Bz45NDhn#j!a{=e%^L6!r& zDU_gjtAl6rvZG3_*qr+E%XU}I2S;N!g@Wh(Adm-xtPp_o(}rfwD$l%4?_q-&um^@< zKuGWm^XepDH^nykr-sP~e@AEa{}*QoA@jCfgoL&-RD$|gMz zKh3^Z;RjU~)wg0bZyS=+yw@G)q0ZX;bCZ$`St>*6}*F!bYUmGRtZpEX=9`D z`Do41abI!UKnqLIm%L{&5hF&(FSR!`8yF+d_62c&z#hUg&i0Y3w#Gttde+%slCtmq>vMieJqcJT(Ju1XwvbHCBvvL`3!(=QiC-jRFl`|*Wh~&L{ zOq%Vn-ff2)4MUDeBV*o2e>zBf%AMhq-bWFM*i%Pzdkf+&CEI%5w97X>CYum7n!!A` znjq3gJOA3(c10E&R+~s8R(W<-ms> ziH=7!e~f`cBbEv`_|0lQ!Hz^GxTpIC#TzB>c5EV{?9ncIs9naaN94&K5n)KAN10T+ z71SvRJ%sACLgnCrImU#lxDDRM2!x|SYCQ@hAK ze;~4`(aGv`lB8+NthuI{u}luRGzTv#ozUgkkdQ|~k6K>{vO4e3lU~#02m$9H9#5qQ z-E7+S-;*!)FNlu33mqtVHT?;)_qT8$0vKH}KjWWY9aY#AEVe8NYduJ%1_4MM(pgq| zlzg#|hDVvON>Z(}z8Q#`^%qWg%7fW5e@OJ~N?Aqd&e@RKM}N`1*pPRYV5^u^G$QXV zwu&)>z5q?}7aqG3+ulK_#|*5xVWt_V6Qe-3U> zE)jX!kzr9SHlGs*_L$^OvvaS?5%f09AC9P$ov+#}Cy+(tW;+=4bJ|VeyNZLxAa4hq zmq|0azdgkaCC%RPk#v&Za*oX_#B`H4x=&zY)2!G;lFK8)$yz>A{+@$p?*KY^cTqHg zRITOxIFn6zD?eC90a@TLe=?5Sv6n{{NOJOb4Ss;Opb=9RqbtX>I45;`uSTVPfb&-7 z)kz=r!oIpSeEdhYuShX{oKrjNsjHd%+16C;r+aWOXD<1MQ!E8z0!aDx+eWLOc5GxZ zQB?B4iBL6ph>#I$_cR(_aik0&-mRj^ue_VC{Po+K$2$HRMaWkuf7W*Ic&@pH#X$D} zLp)CN&%Ea}+&^cN0L0@es}aaO)MT+$Xqm9r&uxP^usi~&WSpf;NDvdk&Y)~v=9#J*`}ZypDMy-m+?Z#o#gfNVBF|lUL-`zo}1m9Ux#DprAg-J2_1x0360? zQi58pv)_p9Dm7eP{Q7Oa=;msPjuv>Zm|hYc7$cd*x5(%ie@g}$tXvd%moWM#$J`wy zEQ;zns9VQx9Uc18tVk7?Wjg1Lt(DxOJKlP*UC<6w?S?Wv)yvtQfy0;iXui{b$ z6Ig_82Y9k*p8!*a&AOh9=WGw&IYBHO0t&EdG=Qll>?L(+lZ&qB7y;&hW@Oq4rv&r+ z+7(g2>y=jMe~eYM0BN#Z21X<(0VGROyJ)>swR^dc;fm#(p^MUHx1cj*Z@elDGj-)p zQ(Y2LeT$otzCL5jXK)u)Vl+KD>12`}nQ7Pfu3afs9+d#g#LsP#*F6QBffb;=DU$ z@^7g}mIo1e!i83fE;0h1Y;%&&dtv5CPAKOcqRn_L$wxocBqm`DMV_b=yqL_)!eP-d zF=y3OR>ugip~AlqniV8Ww<`Y51$SroX5w#Tf4h&@o75q;%=8R4ZrHELmVTZee%88I zHk&tRwsOsgj)@Fu?QVh(?%>vWF|o`VS1eSd6w@Sh%%E8xv=8iZj@+;Kk!o|5im> ze`vyarYK4%F`@pcpG^sUF1mE2f*Wa;inLk0y-#pzC}sed4CvOt5!)DOOARbN+bE4K z_^6G+=02L2^-0SWCQ~1oq3mms*uj*AvSrqQ_eyq}2~6jWc93LDZt!wI z#9m_N>i$(%vq+Bj*)TD}HhtGNr&HElf6Z6JeXN|8j}K}2Nwd>gul3YOVBjQT7~Isz zT+b?rm@|gQi+@is5i$kwY8IV9viMV+!^f_3#ku-T9g%26wT>>c^B<8a)6{O({nB0^OwDMM(Ao2e|n#5G?p7ee?E%$RnWe^QnX#LUy9PEPjQaBBq^;m zjnRs`CO-j+_wOsvpCWF_LT-j@+(R^hE-fAC0s2+hzDp{0ieu8M`@m;dnCcE(qRr@# ziQtf_L~~}^Tr_|C%I!2w^SnG`X^qPP;A>BbeLA+p-Th|6G77?YriiC^e;LEU-(8`x zA&WkGM+b}g^o$KjqM5#{O|Ik(iJiQuwyqK8veZ;!37@Y}$}Y7N6-ps+6-vch6ib|Y z6E*uSmX>fd!@H`Y-fK|%%I|>B4de0eyR{LP%RRdktx~S#I);WZCkTisrDu!erIyu0ckEIB>ak-)9U1R3{PH03Hk8YjW}JF@~%SS_>IkzB_7=@rZ@4n|aSn zQF{j<x{3+|8{rx+ zw^IEYJyw&KVeDO)yQsAaZBC}P@_JGM(kG5b2XWwiD!pUXCwF*m)#A?e&!F2k#WV*% zR2E{+DnQa_@)oTye-^n#WqGgDw(kv(jk*$@s7|o}F-7H+IFXHBuHJOPAj}vMz5FW; z=AB`3!~4)HCxqU>o!a5X(%bB8vw5Cz+KyfDA%=*{hDJ2GPwx4a zY!|oJsomvWDb%Q)5zkagjKm>pzP#-;CF=npBte*wq|vjQf0&e?&*g&*eky(yU%KHi zF;AulAn?8C&$p&%X+E`nwoAe;FL$O;eFRf;Su3p4?u!4S3E5iKMk7n^`;nG1kv`)P zt(w?_x*%8AqBQJz-GeZQg(Ujq>*0BEX+H#1p5Za9tlnV)_>+9)i|1_Cw&bFoZBEj_ z_>%mCys=dcf4;41kWUh4>!DUb%9-mp1c*XkO$bnHpW#`EePD2~uuFL%2DD+sfM!<{ z2J9#4F<=~+$AA)L)#UU3Ave+2<`&V+Ud+b+^nU1pisfuxtpb|~LUz@Z>=#YAQJebu4k6dbU$CZAB+1%1A;dM-w6I|tFn^&jrxYi_v9 zK`eC)z;Lg?JSKYvlawZ5pkfqb#( zOlAhtk%E;vEYfG~Rtiy)e38^$emYwNUKNq&udFW@d*{zwMMPeqQ_KfTF>pq{{z)nx zNYL-%hFYo2!i(vA9?OSN)j_qCGz5;$z#x zB_LZkvJRWKAVB#_!kGEms+fz>M^xF9V`+uNiLk-jfJ?R)eS;($Rb<1?x;%3Wpj&(j ze>3;lq0)|gDdLRn_Gd;%JZ?`Dh52)oA{gp0?VBll5(t+udO@W_urV`!nW&e`8TW0tMpI#(Yi<0!l6dV4<6d0D!YYW%K@ zaGOwz((E!*$3eGgavS?Lt3{LBr(~8H7-+0%{K{xsWADb_}MCX41Cje_*!m zJH8hEF{N`f1<+T!6PE_Z&%OGI6~N5rG}d7_V)oQ-&O^K?UexQ(%{jt+zWXwkt9Zj_ z7G7X^$7pn4Zb*E~eB7Ob_^GTJO)6IRF1(8;{RSOGFNMP7YwVvP&NJ`=SaXD(b}@tU z5>~@rMRMv{nPT?m%7UoFG^UkHe|@f}B}NV89;;R!9X-NpPZ9wnX>`Hn7dGvitI_SU zM@*s`b+*^u$ax}v=A5cQ^2xotQPwWVz1Lu<4{1-`1u!~?9OgBB4B8vLL*Ae{<_7POkM;_mgb)mZ;pDw-2tl1jP3|FO4(bhGmnRg=%8y8b}IsBnQE z%@s=Hr`r1Vlznh3H{T;1e_J&!x5el(3tr#e%bzw`dpB*ebdLL01DeiZ=$OhV<>Iwu z>-=~JJ{ObgvNCgH!!zEs)3tl7TZ^-1S4|MmHRmyn3AAAJO>lt&CcNx=<2<^-)Z)tp z#0A*bH0JJ$zD{Y{_U-lJ-dz-cjFoF2*AiXpLfP(BaBnPvp`9nmfBCdrkh!mQ`>+fM zR@JiLRD0az;U#TAOC5S8?#78aHfq;>+>j=1qR0=1Pxv2Y!GKH&yK>yzq1oY(mVAE? zCgds5T&k6qS;D89ZFZU7rB4lCynC+f|D)ABmKdr#tNC0d+K1-z*BPHss9XAo%GMI& zT+-{RdIp;)q=0^je-ZqK5gR1+Cu9FO{MWcrjUS`^Q>Do-`@aBl(_g=U!KKDeMCBGQ z0KV}fY9-8_Oc-ed>DhU{bzxa>X;fr;w(Rh0-nPr!I2nI(4b}+Ru-y>)ZH_Hf93RZ$ zZNXb4a%L+Uv||05Y5$3(OL}oe=Yh;essUu2Lv)}`)NW&RY}>YNb&^iL*tYFAwrwXJ z+qP}n9e2#z_pbkr{%MWQpw_B$YM=e|AW3!bfoxp;`*bV7LZ9{}(!DVP+*}CjP$=nTwNywUxyQ z^a&T(U)FP->u0CHLHv)@HrPgP33MnN#MWpZOsV+0i_LnvVq$fhEj<)AqadF7!moB^ zim6wp;_c6NNBDN+_xQJ&fPvBPtVtv;()C4V>aou6N#wb*bpn# zpk%pIN?)gBb$BMibh&birxMf9JuX-O(~${4Ts6BhwG`g`0FxG^ z#o0VU*Yo)P2IHTLoZatD5&DfR)f*zhnpKGlRqUk9_)t#A>JY8yEnT_#wWX=QL0Sch zmd~|Fe~i@jqqhK`ww1Y6D+)xUQQG1yA>MRTRikeT#6=!RB*9Uk64#BJZ=beQIY1+% zrTo+mw5qQli}h*t3A)+BHfIe8j=|ku{xZ!GfG_!Z_V;a7otgi42s~6LYE<&VsO#Hu zQKe#S&rW&M7qdezn!?8LTh1bC>fqss@5}8)|6rW^uJ7NsS#_i@U#c3y545k@u)Kdf z{qJqj_O9;d#s2l;rY%*~@`0$o9WZ2c@8+j&{IfyTK@bKjTdD9OrL z6#U9<>}Vz1FE|($x9H^oU3(Nd@a4T_aOH@XyIf+4AiyKt?r|vK<6TYO)!4CB1jpM; zT9C7#c)7_q*F{#D+d}GJcB8)+#4InPMj)7M#-eAj;P$m4-TO`1!zjXi4*ce_6ko|e zBkL$XB(_&1*yql5*iGsOlG1d??x}ko+XUi+Wuez%YA;An5Y?7P@SpO0S9ad8p{P}1 zdZ?NFUIh7R!KkIrr2I)~yLH9qR6k_2E&s$I>uA^YX{tsedSEBZ!R0X6N7{$ z9v}Ji_r5>fArCRSAidrzD2$Q~8gh*qI`M-9nG&BQel+Rkk})_muuQuW9)udKp$`vj zm;t(h??Psw8`WI!Gk@6D7wK@uf@@;-7xbgA8S zBiAN|EI9!3$v|WyOn)!j;dmL7_Pj0B9MqIq4g|zlBs2eqXVdwEJ2x_(YmYZODTDun zyDG!@s);Ss(!f&+KxXeqCQ=W|2trEFe^&ofG~-tN&qnD&W_tAp34Dt`?HA91ebXV^ z*FjO-Cb(@c)ypsy-A80##$9!Iuhx6~zqUf-ZdHLy1IyoPL%6Y3;LI#TV?sk?3;YX< zWERq$@MzpULv+1v1Yr7Eke?;-HI@JM%|h)N2Onwq$39*J7?)!$bH?k13rf6U4r;JN z7_()Z%FoApm<7y!39*&64&@(Zv6d0K*bnE^9tDuoZIM^ft@uQ>k%0IPe33kqbnAj& znnH*d&;*de{HI5mT_Gg`N<^pn4(_{gsVL=Hi|wfvM##xUQWSPvMv%Et@CStH8J3UZ zA*24UMcx_+aQ*~iBwOGJyMt*YTf}l=rV8TIj&CH}i)~zh*XuD&QhE;tG9Yn+4WZ?L zOK1!4PfQTz7YExUn&nj=Dh+72#s(m=*MjGT`p(vitjQ@VA@8YZ5*65jd#T(;LE2 zrv^em5%C~U+TH}Z$VA8Vm|(EfoSODeQtO}{e1538wW4gam?S}dXZ53 zJTmY(p*;*L23KrkPWixi zL>kRbD8xamZA&rWnGX&LQ=IY(MJ2YM%TJ#(4`yj&BY{Q2F=PISm_J}*P#!E)YeZMDK4rh2N?3$-eQMKDr^$SP z-f6lvf0f!ayfsQOr#^qNS2&|Ahy_&9OEaa^B8&xaVSi=t7Ma@>|3(yA5q74EYOU3s zk2@`5lvSxCkjS|1J)vBT4It+!w;)@j26?Qk9f`02$n3ERu@)6AWtyA`%sOlmi616>Sh6za(9g=fMRFB>8DO9}xH#(I? zVR1y@V3L33_R?2iDXMVL;W5j{j~kFY90ZFde@D^r%ymdOAy%tgnGT0c?0jnmKY3)6U9t^*p&-pIv7 zz-U<_j=3xpyF^s$qUQuqp_@>$t2c-dH#X1Eae{D~qEMi*vA87=J6Ny*CdwgDD$E>G zxLVNPlD}a>$ql;s%ySABR@roS{mgY5R5;$H>8GqiHBgKPpoppXM{mpExYUa&%6&=Y-lEPb>~tTr;n{YY z!4a_?wzrG%RjqE64gU8Ac(E$q`mia_@+P8mE6kchZO-&y_l_|BS>sDpAD2j$C)ql{ zNN1|+{Q=RT#r)B`mkoVmjNL#7taK zN@;MvdQm~rFK~t}kqDZEL+xbh~aX2_GDHgA6a340hi=@hk9$o*< znjQ%f;^d!t9qg~@X3wr3kB8Xj>(gY)*B@JD0GfWnfho0iVoVsz0l&!BNX*C3ciL{82 z`X`Kg@H}P=KxNV-au}EH`&SvRUsl=ccpGYw1|cA@q6KmpkU7flv~D&cnT0RKLR}gN zhY0O9fDyJ<46jm!y_yYE_2HTDA6kZ2o#1$*EONV0_k~`)b8RBYy7{A2WcWj)_z?hU{JUls;Mb84M(q5QNJ}6 zf)pMq^w>%yV5ua>lCG+MF;U;Q3!nDgoX|J$Ae}vIHjbHoFXcSl@Na%BLUu%o+Jm0W zEc>EENV;5>o>}0;h3Lw4msJq(c#}l4ils|zNF)&hLs7W%M1THcjuu57Cm&c$Lm0)h zCV@tKmT3oNR7=LH)dqPp{PM3UTA>QNkJQT%aay&;1C&g&hS-oz>QzY0wtwi`rDLBfTd8-Wt~|@ zaJq;ASk?Zv8f*ts$_$Y$ir!&Nrw8}=_ymP6O|J4pG5_6++h|iA8}9S%8T%h78nbHr z+Cc+9Ib@Kv76w|NtkEmeS%SL8?_Goo75)?-p%K}z1gGXg8W=b`-oE*G#;qSYHFZg1 zpi2Uf;-=vj(gP7kIh3a*TXNbPDE+XR%?}#SK<+){f^n*@tpdE90((pbC&~q!hsoU7 z{u;M6l^ZMftqXo-wRnv*3$_{M?HUZsqkX0GO&p3XBQQ0TZ^qW8NU)aJU&4xXYTJq#`}=!Pb~Z2?V_XN8Y zK!pm#Jg-%&^hnpr(p>s|Vpw4vFWrhAq>R&_L*9Z?)=@emV!|}F>GQLKpGLI5tW2c} z)Jzm&PTiay$Qit=MQjFgb}&K71y>bH5<`v1EN$6n1b!-SV(A}Q8n95zllAB&Om`|m z>iWv`$?s@Z{vdvEg({w|uc(lUpfOxb0|pWn#jEZ)>%E83%j%m_saeN-NcI8g9i+GzT{`-@>Z0JRZted7MBh!cy-o!8KY)Z^aYXGYJ1cchQ~` z3SccDWwT=XQ6sCGsb3ea!Qu{6MU>ZrgHB;5@W3W+nW&z zV>?O^!S7?NK3G*LYWjlJtR`jVFPrrQxQov$2cJ2L?>4o(zMb&R=9!`#!NxDtee5GV%?+(XPgnn8=EE#9=JPb{CNWVIQMLCR3^aUZDWHa1t&b(lG-y1D0m@o9+kZqi z5N>BArFFmTKP%@68p%~HCm12!co+w1MQTdka9;S8Gsv_ofEujc{yDi-YV9_SVU#=Q zJO~UVnHQUZtGC>BcoF#t;Yba7jPz*xI$ea`R=4bzsFeGG2>CEVzlH|K?Rze+dDK(& zyIi5U3f-!Q%8;B>7@xYsI zEd(bPCPn0O4UXzksX1zUaC(;$k%tsJW}$1&^7F%|N;CyzQrLbhI}`h}X2ecYZw~w4 zOW@ph)@?&n|5aH^q>#g5ID>4f%jH0hD?s`y*5{og^JYsTg#u}UP4l;@L`7pt$iM0MX;V?d4zgu8JLnq9x*R1HwQx!S7 zfwiA4TwN)}Gw~Xu#W8tYqcSD6y(`0PdkD+C2lj`;{${RkNM%h1FE!~o&yt^f`*7pu zA>W=%l~!03K4B3{15oI*^*TFFnO2nF4%K_NZQkW+qF~A*7Qvn@%4az|nkuJHaT`=lE+1f<w74HDTW1J0RNdhtV0^beOaY~ zleg5h;Q*qJbE7?ny|N?ko$zfRa~f8|4s8gL`W*E>1nFuC8&GGD85c<0_A?HS|IvOW zppYyDqoRD+>4}IOkgS-{VBsvjtqK#}arz>Bq%-o%K*X2J?GEIhlh1YcD!3!@a07cP z?4Hf}bV0p%IFrPsaPOCE9v`GEX*HQwK|;vKX~J=mYwtQ%($wBJ2l-HAcXdO+S02gB z_-cHYorcjt2`titVJJefP$#I5M0u!VF>=Bu3B1ne-Q1ac5Q1BZZNR&WU7Yvr{JSZV z|G|tGRF;mG;59Dhol#34{esHRpQRSCE2xk4x>%OTk`x!R_q%-l?k&4$6lzYb6pK0mQP-XJ>(Q%@EV1c`sric zRpqYhoc!X~nw*@SAjNcSTP;1RE6|i0vRo9l_ zd7@Y$13W0q{1xNS(cI?AV^ny|1HU?-q7u9uB=62LImHDil~8H8*2uQaVBnWkI)^>y z`No!r22!IBkz#-T%?9s^PSPMFz`M(6`4y9FmuGanPIcZj3udG#Rd<@b9}KBEeDb=` zK!|2H+=_O(^gs;19y3`*(#nqZ&pe*3Ur!di0TLyyf?4))?Rh`gv$&+r;qVqZ)RVMY zTE z3;}Lh-vZxL`g*h^PrH2i+3knFR-x_tb3KKaS1@oKPpJ|YHBB^PP&lwMZoY$b=$UsJ z0Oe#(%xDe4#GUpvi2&ayzSauxr-@FM@)b6FOjA9B5V%z4&6xD}D9j_W!d*MDGZz-K z;yM#63lMJ*7Rup5AA=5#%~H`TV?hV7B{3xiwcS)MnNM50#CpkNy!uJkvlApdbz; zn(SqEyWHmj*%Nq#<#oGIfd!x4Iir}-NUL)sqtP~ZnS|TlVHc&om0a^~SKRXmMtFpYuGORLREte7h6ovfMUtwTdi|2y|}fQE|gfz^|R+)n-HrCL?vVe&-E zBd@j+leAS0V;x^Kcq|P*jpa`h;PA(c79{aqSLn;E0>ggEVf+4(>%SDdnel3G7fGT5 z^AN6p9cMubrX5Gu-d()$L)d_r%B^>P{`#Ksx$_rh9``hc*J_$ed1wYT2@g%SL9;6z zsuuK|wWMnA2)Hc1TyzvvvZ>fomR;-{|GXZIe_yKI!XfaA_TYOK}xgs3QsD84+&Hg&?Ljz=zbnvySt7kgI(SI;WCXUeQQz_vP zlzhWxT@Z!Q2qa^)4QFbF(7N$gZMD26cFsIgKEBuGOIOtf^m<4=DZR3i!*E*6BRZ#a z<4nBBZC0zVxu01u2IZJ3KtkoV<=x0$2>rL-1+jvE8fvRV4OtgFlI=@WobH~zSUYnJ zfym!VS;iZ1efc+zo+w4?$as6S4#>6v(WQ@xoPLbjE_}b((}A?2k)dVZ8v3|ugSOkA zi;g;hx)EnpYx<}DUv}~C&@1@Yo$xiO&i@g_mRRVoL9BTA({BCH1dKlv?Mbb9ZaGiZ zMTyO!Et?8*mH#lOcDJO&%>OkP(ekg`;T+vv;U_?yC>OdNG=|!r^y-pq#iwvrrPnXp z+;znQ%3Kf_tVHkk6*rU;SX!}wO{^?>T;NnoQ3YQw?em!ay5s3?GclO$5GV~_&?j&kcn=?kR7SHcPs?8Ju zi^W=Ur*ZGrYen4K0=rdFTchbB&!z*;BIYEsWEyDPRIMSKZ9uPR z=+~U>PwgsEqwG-B_wYZa$Emi1jP-$Cjkaxdea4I^;j<$W04_kcC^`5tV8uL29u1si z_+vxKd--+hHo`MyHYIq|OjXxDVI}?>Yl5aWfATG&3^3w%+D#W;+0U%a zxV-+^;_k;Ju7R>mb0+(P&ySVA*w^&#RYhHy|2eNJZhLvY{0_M?Fld!im$W(GW5aiiI-g5@6D@49G!(cO>*vx%Nx|P zysEtG<&lNDNdMkVDAhxPGmi(@#@u3*1ER`SJ;!A;KO~ zG@3XIFbUDSkQwJe2lHLWYR4n4Z0$h84ttm&$?^wO=8QyqBodaSmuHxq6uvhQ9@pFt zxgV$5)9DT)Nb4(THfo8LUz!jIeL&g;HoN4;9@RR;?fney1~TpQ?yqz|Ay7$vz93@D z6$x%jO0d`yN=UE>(E+630kKvGKb6|$Ltfvn6=xSgahKLKpsBHrmuCK+YW;-Y8#DO& z?Fur+ef=2Z>VMmIbtVZ#K%5cudvkE{kYECLq*-CgIGQ2QocyHy?`C2x(}_~5r?yfp z=3i7lLqn(DLm(QSKj3E}Oc8ZH`YEw{V3%7L#wRMwNlFcdDdhr|zKhAbXrYwvT%YF4 zKJ<%=5DViy?t#3D1~M!(o{U;C2n$&NgeBSsO5TM*ef;jRNd`*R|1Isiq8O5>;2Bq%6^ntl*<6i2C@;@2vxpKqgKmBx{n37CJui(aWWC$wql& z-aWSu0eC5_8WC1*ZVQNdm1t9?^q^cM&g0oWoOR#=OHb&6Vn*G|3cTrzw^Dtgwngc| zvme~;Hi-k1J%;%{%^n`nI_L(LO~jnD=`;{Fm@U4MyD#NQ_wm9NSZ zD&q?`peO@}>09@k7n$d=@x4u2?8!Kl6CCf;VLb*QgO3kI(fgq3(9I@PB_l*n)>n~u zfkDOm?2PD5-+LWsBu$RjMwUYHZ3e<QFL>sRKPoKs8wh+$0db{pouYzff+pA#7##Rl~VzIegbH%U~`m#XNgO!fT&pr z1cPJvq`OiEqC6%wR*p!=bhy0X7g{iReX03aKalAti(V$9u6qKAwezw*D(hb=-ba+@ z5c_o@vC=bS5ihHmxLQya9$=7x((1(2Okx2qq_bGm%w8>f8LTh7P$vjJwe!p3z|XLI zn(}5CQX=hY8DxjJf`1#H%k;hrXa3@ic!vmn=p(dJluiRevns!S@*(97R%omIG@>A8OO8m9^JvzvHUX|^-;$RSsUBl3#k7)HRUl>zh zq8NhBV0`VE5#FEXO7zNhCxme#18CNzCWmAE+73ffN!(Exya5;0fmxHxWf~Aax9U7q zV2T?io+6u?rA0r85=a~AIw3uh!y=hjC}|(Ac`QW`mK<&3mjVhy&rYv9Pra{=yy z@GGoyaCZhPYJ7uV2Gl91TjG%~$PW@KCfkNtCBZlEjO48&EFP023u*>2`s@fhDwhfT zJ-iYubb&2Ht-qg`#b$Ihf4rksrrXl3D^6ml%QYHQYI@ z5V>g$hFlk>5n146JdysA$aunhFcki5vAC#5hl<@It9GXcSfQ?vi-4?*)nZe$iq5wy z7=lKr>`cvfZ%}72wBBC=31;sbtS?fZ{bYUaF&|LXHdljr=2Akkx-tPj9EtZ@so-Hu zl@-kR|NQGjF`EC1w=o-Fy;+fpjnQEIB<7x}h4;3q1wY#IR)$)H)Zr-K<0*~)h z90GUpzw+t7IDpLHEzJlaSOHS_%g{k#Lu-o)k)6Xgt+n^H3%bs!Qd|GeRw zsX;vizyH{SyW_>i0$4PV70Kr*&oyWe3gAiZq7D~kA6ok5n*pX+ummCp$h9Hykj}A= z2XJ1?e768uQJA&?)PZvub>T)l%#n0V%taiFjtUGh8G6+m4TJD&pQkbu01c| zu+kfCQsW8BBQqI+yHlCrZF@UzA={Y$Om)ZPQ197fIzgGM$_%}DtD_^Ra>+#Ne{nX> zF^Mv$!-ff|0RSE$uG9eXW$6-u38^7tZFTtO*r0*lOMC=`f2L+TvzVOU&uSP(X0hG~ z9=1eFt8SOUTLHQI9{%$SlVGoiq6a@8G-itAK z!Nx4!+vUqg3oeJm7K30_>bcSC2iD_)CXIGsJ}T zU9dkNR|g7GZLN}9nE5UZXD(Tue{^zO`#-8|5Cjm-c8MDGJtrr$2!P=Qn`s4ABJ@aH z-BB8udZlnsB(Tsm)QV<`&&(+?hKaF?n9lqdNn3kwiw= z&I-{>9q6Pfi#iox{x4LsN08x9Aj>|>t3M@vkK7X27Ro0ZcD=4yG4{)FwDlp9z@?tc z5jV1)s&IQOmouggjc>LZ#g+7pXxGazzE_Rdq+#|HLICDY4XxsMEd@?IGynct-QZF% z}T^*vmibF%Q^9`#^%DC%QqAL&3Ry!akTidB5gY!vX7 zyD*>3_i@;j73q+^IIqU`)siMj9fzrR_4jfQyZb)uFmkcLq!>x%A4!m z31DzUKf}t_-?WYFIj=Hw2TQGKhmc?jj}eUUV<@{0YC>qbZcBDdW%6DZgHcv^s2=2N5oTTKOV5 zP5B3OcpS_q*NUr<5FZ&_qG{0Dxh(5fZFcY?60=A{0Fdm zZ&PCqYb_xO9=T^!Jr%KAVixs=pD!-*TI{@g2{0zmG7{t)=GY&?xEICmSm+7mn_nUmC{=^4@SeHEc# zS@RGJUFh6@Lu25IBCNfdh=WmQJaG97wRlGHv8kza3S1&#@XYaeqPj`th#BWWf7no7 zkiI}~eLlk@wGUl_{izgBUY&0ZZ2`f!)J+s)_Wi*-}j zV=b;CBPGrwUsev({Xk#)?wE)7kjwSlkljrw9@JUW=P>eRBc=LZbQ}vw1;F#Q! zBe%ype-q* zLC7FftUJ4f&TlId#ZhnO!c>%LwAkl=sfFXe=9isy!*uc+0@ENlI=Vw%Z(T;t>P9oy zGZI2GVt(AlBbgGa@aU^}pac?6BH|;357<$CmUtzL!qDR?W66Gh`vcR%{dI|WkG?P* zzTS%1?}83jm9R*EL+h^TG3)hsXnu>aE?p-%36hKHHb#Eqs3~p!m+@+TDno;L(*Eod zqD%tXs5>O#r=&RqO%7;M5gMES80{xyRcAN=`l(_OD-bZ1CN%WwpM8|Te29 zCoeW8lZNIfsxYS|K!73xi@&VNs%flg|v0mTA#C!wx z+K#$E+*-@|7&Pv69W*;zMa5D%&5xDsT{tT}|7*Z^TweQ&4+^|V#Gb|@w_4rRrZtX) z^TGJ>kd4vNYKIK6GB$IR8jA}3C^ZJFK8`|{`N-i6FczD8-YitgC9Kz)WJW4PQ7Jsh zWXj~xS$t%1quk1($A{lESVn8mF<5q=c&#~mfs|}{j_E;xa!>1RIm9wMC}SAsI$QP6 z8dh7frhS6HtpZv}ZZtZ_1CABUAI4cAKm(LN*K3oXbFQVuhsSUbFdcli{lhyaU=t@@6j4`bmdVLg$~6y5sXMAziAPxR<_u^uYcig>q&4h>DPVuew5fmU zKbM^0fax0P>(!OaAfN54es~2aB-yt`JXSbI3Zvc!-pfI~QrH*N=e}hOz+qHj%gc24 zS$3q@gsJ+XuOni9l(eVL#%ABigi5}>x|c{yegY_P66lxX%FvByM{e_qtz6aIYX5;A zl{Rfu+#lOuA(@d_Y=*a+Z0ge-aUj`kUyX4g-_x?Mmn>;B*tVAe!ld$**mWdY$B(>4 z-u;jvDO3C_d*po>y{AJfK|LF753Mevec(yBuGK+sQr z074++zHsAhO>eQe{jBWsVa+ERC15g1FEsb7v6FS$VH}=(1%&Jvp(*nn+Iyly9s%cx zyE=!0j=+P-K^c_t(X6DOT;>P7A05!&l@AJO)TDDW!_UD}S4{Bo z0gr<55>O83 zBX$VZ{?ojCU%XwiSZcbl;AFT07lY;}^2*>LpHu%F<70mTPjJdAsLr^lG~}io%X>&W z2t~@$bf{n}^YdT^E>UWU1rxxj59qDj@D6DFZrMKe9pENp)wqyXh!;vN)GVkVKna3S zJMts=`80f<jdn!#$?8LLzv=}Jq@6%cCla(bz~wKXNCAh9QjD5korPyWIhgb zEIhvR>YT`rLEEC>#o5%K^^LgQ$g7!fUraL^55IhBIVo*V@-ExvUU8OpV2)OhD7YVQ z@_g)illP!3c4q#V)z~>a8)qygzX*7s_br~W0~*}u|+;SP#0&BqUrVI~XKck!IY{Cvdk#;&l zGybX=HO$%RhqRpdq*?VIN%<)oW`djSP!7cF)1=n`YC#6d-|87IR7hVnft-%9Zuw6; zlv;Q!(&vnw6+q?Pyr_tdt+l(=_572w732+7HHG1H3Hw>S()>=5VBL?uqTA!Bccuyo z8p2A;LM^&@gIPQ7uY0Cc;gYkc*}LobUKPw0=JKt5hDuQy(0;taD%J~Sm!FrlYx#)8 zJT=#y9~)74E7?(0g;oOmYwWz1mpsry!P1#RANw9*wg8?J^Gb$f+HgzfCyYq#RH@Y* z>J6t~!M7LFKt?$y4gJH%O&*1^JF<8h=3l$MFjuO8I+zu8X0o=`re}mp5vh!V82p`* zv&vF>DLY_!&ZKP30>*QHw}jO{p2pt4aY^Y3%k4A3^!YRy4`WF9VI`7oojO$inGFSP zIgB8?2-rL6BT8|D)lq`f>x3Rbuo-QFt>weF3qi>{GBCJ;lizQ$tAuxcsd-WMyRc}o zLOZ5xmViNNc`aQPEc&=_B}-j|8J%RVO$L88x)yg*pitorFq6ULu&PKhq{ZgWVDbJ` zhzucnqONGGX=(D$D{u-- z$?Uu@D|iyg6}4maGj(SPWrid!IWUG;Y|j%AH6cMtT(LlpJx~71K{GG5GpGd)IkU}$ zT*Zv)7UnH6BQP$cMZv#*L7+Q%;44XCjB%o2+xzXMybJzI^JdC?zby8b({U*IV)Mkw`w7FJ; z)R+z&Cl12M9YMg|S2O$9%&Wsj&W4>BAhN+X#*fL$7E=O;06~K(S}Wumf)_c-tEP%H zb*cW4wthf(B^s2UM+v^19f3gAex8k~Uy?o_^?dbw$QEb#*V5lq`AiadCP$~bRe0-C+%F#m^N5_wXmHKYW%5d8s^{7;is_Ff0CK+^ z9SX_p`&)VgJ_dqb_1oEfVL39h!5p?yWjm^lss__H$cMyuJR|z$;inCw)HF=_Gp|yI zyNZ!xeD%PGc=h31D?cuNt9 z8F)>)%H%VTZ5}d9(@bYc<)j7_qLR`n^IN42Qv1`LMh5!-8kP&ixk~L; zPKVvP%hD~z5F@Jhvl^AqVtIL^%hO#+@v&%mo;)-1zIp0HcnfM;Gv`wK4Qzh&eCnx7 z7!XZy3DE>LZL!91vkmz@;2*18Z?L2H7ug1deYAZxFwCVGVDv!U zO111OH{$xfKxfB`wDj1adJK?l0;pw_0mc_HM|Y$h%FISHNC4@apjZ3bFr#i`fOMy( z#v6rhda+ZEPM3YJdE=%zAXpADAFSbaY#qtXiWC_TM&_JUf^4$Xp3mYa`T0XRKFjb< zK4?oD!hhjHQ$Em8Ih) z!tTnzkbgt3%~tAC#+9Ug51gUS!{azJ$ni!-VM`Az^GR2pZPvk31&TT4nx5UN9~~5b zgNKMpGvuh2Lb0$kd}4*E!M8gHBzuHMp~h_6=Aplt-JC29BGP1yEnxg6e8QF=cyN($ z%hg$bYDJV9LrbD!j}s#iS;6>?;ZBzv?aU&DGQ5F7MEmUwVbZK-({TYyKWAjw0T|-d+qus`d+&lYt(ep6 z2C9`!A|MR*co*~Y<&g|!#)f$wk!MsPY6qf3!={oE%1A|_pxoIB%7?kkGPZ*=3{^Rz zi$D!tm;DH?%7>0k$-ZN()5yh@Za|K|sc>0oqgP#p;yp5%Sv6OepoP6a(V1&JSMzGY zs+UcxxdVzFfYU7k2+*xC)CZIkg_&MBqyuRw{SY6t+Bb=5CGAOBb4|FH$%ISP81t~; zkHa;2o9^R%CS4riMC1!W+4)nhjplOAvF1=w*#k{j70&-|DCRVSIKuEtnd|A<5u|grc7g-E`iyf-0e=*39#p9P zqV6$zpXomZWI}wIwM=#+X*$O8e~Bp|!`5D$2d{$fPHjiOX~JsQX~)Jd+~W5H`N!$n zLe6pB0Ar)>|GYMhS{>K5Cx6F~lI%S2`a;pT7!FSVTyAW*cet)kxEkz=lMQA>nmKoH z8P|*DV;-$rndC_a%Ziy+8!)yfQSN>+oloEX%R3%jJq0MClSzI`kJGHnYh_N6c zN_LNYCS0H9>2yW)O5pXNz-2$R_HnYXAtC39pU1p)=uE{_$Zf(Zk7QHYqXT^j!}QtF zP47o6fJDds(W;m98UwD%%{&8=v;T1F!=HEKV)#*Wl_P6>)?!`4182|XQwbxYo0ers zfO2t=VPZ+^G(a$O!exBDqWOp5)%OhdCo9rVi~ZgZ@X)B56o(pc$_pZ^Gt)tABvLcg zhG(3ww2f7Qr&n{XZ)hE^bcO$I60&msk4ebQ%*FBlqiI>#SpMgfpP%WUerWDuktRF{ zN&^VUq~S|=3AXaQZz=v}6yn8rKiIU;_{0vC7=hl)HznUku_8g{%oV!1iF%BpsH~7n zvGw4ef2?37W9#va%VS(0VER|mm*DH&!I!$u#N?@ve%)l;{_fwXJVuYUBi2vP`@`d) zFo%rZhp2!2e==;A@!;xX;XZ_hlCe4g4&X{d0=izFB5!IQ5BL`^CW7N;yk9>DfiRBZ zk|^qHObP+|hg+Jh=i%!O`e>d%H9b3b54=(IGt?NBTQrh1jF<1Xle>lIF+k?q}#g%C3tSxW=ez!9W zBbjoQ5d2#t>4ZLBj^B&DI$!rESG$bY-=`G4ZQrn8H~m;Zs8F7t!>b||k)PM=G41U} z>&)EYPkJyA1B@=EKs_(Mn;JjiJx#j(D_|&YK;!sKQZVqUpK6r8`)OK%QNh`)#&4(J zj(bsWNR6R=eIHnY1b-;I=| zTnAd~B@2@8mHPxuKJF3N-00;B1Zua}T|HA7yh=L~P(Wl;Tw2+)QE*SvQ24Yg!u_qg zST4%Z;-Bl9T2n`lL%YOdPz;MT$9!?|KJ$QuSq6VvVkSx zzes|X1!f=A4P@Oa`QMHcLHZxhBSN72(A4^Ymud7mP@&jRFJyOnLFw~0W547`(|E8H z>@&nbf|N=rfi(OOgl)~m^7#*kL4;F|fp25t*NYvrOCLXTea{>8v>p4rMQ-AqwOm4RYjV zj#FW{XsHb^f)!QSME7nQTxJGCf9LVEN8MXZurv9#N-7lPA0{4vGYx90duWw7OzXk9 zqsD)MYu?&wwXI8kE*kR&8!5`1j5k{0^{wQ3Z0dcw^6R|zak;vj@Ajxp?% z7@^RWLu|UFJT*RE;icjxE|~7zAQhK{4ccOp0~r-Hu=HaqCA5Scr2t$Io${=$FO%7u zm*fAom|<{5LdFmP<+%e2n_xlMa&e27fIsw}h_c|-MHLC~tLr~;E!aDRQ-2)2>fz!W zJjLzFl~)yDp`zMgE6TQOjKY55_P_=2v3|JJiBKjM?svj2C4*sZ5n9OxFdjRBaiGf! zp|7LlQ{o&ZI|vTO_0dHW#WPSx{~e{M<>Bxy0riG4;oay42u%iF<@HZ~HQCy~1wtAQ zSIh4Tfo8XJ5o1V2K~ev4@mBLDU%6_vA2MpUZt1m}MLj&h#b>On2A?KVbMPA|q3|$vVMI9~*Rjpy3AmZ6H#svoq@>dxOxx>Du_#dYgrFYSzCy(1M;y7fv0Gq9t(J zAHN`^WbCAX55m)-51bgu>X zsENCf-?WB(O6X{>@b+8@vk=`-^jApp`6w#^WIPGLEZPTQHvr8q(5f^%O9lGcG%+t z-%hHJ7egm5Cfa=I9nnSr{NoUvBX_=N2qYL}h&CP$+aKh6$WhQ$qeu~qUfMVDd44p^ zRwD-(L^Qzu!*cDv!#+X)Usd^;DVJ6AQ_>x+#%<3_%bG3Hk%3vNOBxH$&l1zVWYL2~ zvVmgFf;FrpVXEKC{e^$1!Occi3jWqf_u5N$ry7X{BJeHbtzg*|irC3zsxnO{r!MbK zz1xFg<^uOL1sM;sP4xx|Y+ruJv1)e)1Ze_*REgFG_W>7THYrgbKKE+@FP1Tbh}cwJ zk&}tQ1htyqN`tg>+e-=B3u6cz17{3zhrH|;#a0OmNnc<9`X$S(4(9K<@xsR2)o)bT zyKylmYoX~OP(^rA*x9=B|C-3tZG%fq(6=lCS+mImXS8jBbL3R_`GO;@0GN3lbtIHx$}?9Qc39{U3o3tc zr2f3y<3HE3N2j!55#qxz70V(XLX7_$({q;b2Y7SbY=*$bA;SuDmf7?f_NruU3G9C& zVmd|^X|#4cFJ=>V+a8W#EvY4m1hm&7s$lonX4^Q?1tV(MD>C~;-jh?hR@E~lxHTx4 zdYyHCVCMg6>nfn4?6$T@r*xNeGXs-!cPmIscc(}U-QA5KDIg8fT~ZR#NC=W52ndSE zKi`f2-Vg4q#bT{J=RJGxXYc)-J@d|a4|5*+SB0ldW+&4!7#J2o2rQ*DG+SETs?Og)I~vHuuyUFFivdLjDP+4W_jdWf-oob zi4xDN^5sB8X%+f>w|AJ~WYKbCC+=!6juB4=zSC6U@sR);1{0F&qezaYWNvSHEV)E|E+z2Jddu?F(uuL6X4e6NT+&NcS;I>XJ=>yo z+k2ht!lT>Gpb_D|(HWWdB!y>tXGB!}K+rimO^dvYYwqK3;IKxd1obrT)vR26n2HqU zAJQ2<48xV#+%~@*d1327A%oK=VBsarf>$_O)@@ekv!Z^gVVtizi1ojSIy`{qO6+Tn zdkf*x>W13mvYwO17P=<=QrUXQxW-Z25GWZ07~F~#L{mbbVG~8_ej#ZYZ7jp>>hD$ogm%qmRm5iGX8t zdPFyka8wxaonB@FtdWb6kf-zpdg^tEUOmXzB%T)kd>FlU*@9m9v4aA~PfPr)-<>eE z^W|r^58lq5QgbVE=wPPrg1er7`y-ZxvQztaS zz-@@52~%R%xlD1eAN1VmV4PItkMK=VIqHLyp z{$<5eiYUME&ucM2s2O0gOe?UKG%6{&ux369OwaJgD?ZROcuJH!oZBq;PTg`&1xvw? z;QWz^f7}(D0H(1_QQ1^EHgb0sc4>gpui%dFhg@XNY~OM#ds3Z#5~1ik)*6(0ODo8< z&HGEB_Ng=~|Lbkanxf=cfG8#Gg*OdW297(&n>d7u62PcrCL{|hp7Q7MbnLf-vT^E0 zzui*HcD%Vh@6lqQUZZJ+$!pFN>s{~xiiRRCH8t~Z1FwDT=84QI9cr$*!>*l|`gvqH zYcSL1=-;;FY@>K@^)XVZW6hSdCxVYKmtV=i_y_ZrYh?1?3V(bjyz+)Y;>8D7GPj%f zWHDUCHr^{?i7P=#V2>zir!1342a1%15a9i9T}B0vjVam1+okLrl!1{o%Ftwta$VX& zG!kvyKGjH7NAj`-Tpvv-`%a1_Z9q$hib&%77gVxw_F3<-v|pp=jQ4tmd*!J4Yg5)$ z%*C-z&yRO5IPAgQGudlcduH5P-vtqy-G;I(FjFm=2P+#`o>>1eeVH%qrrO}Tq+MrhZ4rFC0^62F z$cym?xOF{;))8!*grByhQL=CLwhY)%PE6epG@bDDfJv*U*(Q%P?1#sj^Yo_+!*^4~ zT2ObwH$-qrc_p0cPOuG$^eLr(xes;H)NkjTX@9UiAb`}=EiQs{3g&qy56yaf}xY^f?FUv3x=hlCT z1|Z!qqj6ZPU6uQzZC0 zZWOu&Duporc>e4Y_Sct;`Ec8n+;kc7-WQ#LMrQ9g5Q+n8q>%8R=3%N>m#cD%2;Exg z3A&!K@b3?BYT_&>J{)B5)^m2}a@apTVisjM|57h!zPPNhQ5{1~i1T|@O&1_LU%Y_* z-8!a4n6d&-qTZ-~ixNe|=EDvjd*!n{=2cej*O_E?OjZ-^;piD_hJE^tIA5)rvtqwu z<|ymBScJ8C8jSj=3cRW`Np|eCo78K;N(qqtbbI?YSa+^mD}e> z8j+)nr)s1U$uEyd%+*%Hv=dfp2YgXFz#>7j#+TZe!Z(Y=RzDGmv2oq&KdT*IU29PB zr};Y?1|K_W)915&sb^409tWpCHREc)F_&pg8g#p&c&Gg0q`+F6gSFQv?8h|1){?A9 zDF<6WQEYcVNR-h=k}4|I{icAf^aW73O+B6(p$z2m+1p$4lRwebAfAZV5&oW9D?#?1 zR=@9k`x5d1SLGEe-|oEFT7MRfO#F&J)!r;H)E!HVP*sCkQu4e373T;*g>$Ru&NnMx zYR?9z=F9F@tr*Zm30+|prS~=>_)$O(U64aJE-@(q_xVTD&VHy_2W zvd5tx2!$MXEHR=A2PZ^?PpafKW>UKaLb%QBunxuAxuT9p4z^SWc_uUvgiXmcw%4rs z1umB5I1K($w(+=~bk8!099WgJqSGt<@Vxc-=@=%&WpFB7XbZWrkH+V$s`hRlUSi-* zd+>lMBOc&39Vf4attW3q3QUjKwRKM8X@ltX#FMli$+ET7ZwE%-7`+=mpPv=WOp^`glq@c zvGR;pf%(f@ckLC1kM+WA6+*zH()6k48(xK}XLEel7{@3SG^_8gccNIVu6%1_pVIC$ zv|1M9iLs5^xOfH9d)S2bN)Cgt+mo7#9a2{ezCAP7&b9y1PLCojpS|R+tIQ!ZNarwR zuC_=o!A_y7@6XJPz`+!&_}UVC^JFpnqqzq?lc419{QgKc5w8B~7R;`;k9N`1phk~@ zJqQrQxs-RpVaE}M(epTV(%i}pg0P~pXqNM*-|WzBnDa9ByfF$e5!yv0&TUn(y6hZm_X;KcF3yl zz1~fy2q$2 za`*&v28ZO8u6@=Ss*p+e`wunQ+GLy;RgqRnN`d|9p}t*WlD0qmHr zQRRAT3psi<@OYzGH` zSq;*JHn;YOd>Jh~!Ug5k(x(@ZwQ5furkCQXSl3L~&!SG`l)SEV{4}J}P;3gzCWveWuE!Bi3?>ucV<;H^)l^d5@_*IuAO@SA5A{{AU!Hf4TZ!E1GWuU zhj(*VnaP2-ID;IFGlntxA861`n|8hObbUs2_eXPF!gn#ZG~-hT1^@4j5!G8KsyG^F0tdla zdU3{I^<#eegq(4D2%d|%LpRmGe>J~~#H0!={>kshx@pNg7lk^jh9(iPWwudvJePjs z6j4w?f#@y!{(9%jT-w^1P9j;T`)#$K@Dp(^4hDxN^MMK9cb_Q*TDUy30VU-x6$R)= zP+%0HwGz1LWY_XBuSzWyT**EXDY63SN(|5ec;cQ|NcnOAja1D|Z`)XmrI* zewyrh#E;DlW{2NIKbM1EH#a_ib`{#supzY$GbwSS?V+J%2x;cX)214-t+=o%Iv=_& zpl1%EFLNdiTZX0=fjevJK!cx({Pgx8z)cW*;suJ|f{xtZ0Vv-@+a=h&Ye|&NvQl?W zQzdrSm)3ZV5|^-+rEY|^iB%GmBe+F32soxwmrvR15n1}&mso(dFefKV?FSwt!&ggh zD2|Q|Ft)gOn$IO6LV#2FJ?BLh!{=+bQ`ovnw)bw*oX(B^3@tO6ujM;ntrBy;CWIOh zAjcW2+ssxw5d?cuU_d#%vJv*$wJ~dzAYa$_&jTlPUePV)P4xI0DXB4LU(271 znuquql<5#h`Eq>FHl;GJen72nogIeX5rE&U6LqbO{hU*eZS+T%+7FipE`_^QsuNWT z8t6s1H&NDu4O#7MBkrpa*w)pdI;#&?5g7cTWmr?$yH8H$#BGG~@JY#S-VFw=w|;Tp ziv5b}+q}!ElbS&O0{)m)ChGKoZRvttqCxhdwHS)y77BO3+6x`n(#N>24Lb+R9Sc=M zYk#&_9KnhS$)PwB9dkkSL)spvzVNoM{(qv zN&zL0cY9`sM=}yYtRb3)al80nX#}}--XC@BWG3&Rwgpy4Z+_GHV;W*31I)eEk^lt| zSe^GO%#grUb%{DH%VG5e2kOQ7yB+f%uOU0(vxx-@SL<&CBcA*^r)IVnngR6MkVo}< zZ$(Nk+X*N<=CiNiTL|&MQXwiqz4SBFIt$>_OC-=q*R*#NEBO|x@SAeE{{v#uaA{6n zEkZ|!$ROUw$$z&NmrwsOiL`dtwEl3!*8Vt1WT)pX_ZCd!8(u2CDm_hj;`aIFsbQx0 z2DvOTL9 z%P(Ax=Xdx#t~hNG`8v4VhF%}hmiV?z8k-xF1$PDOX|2W$v*USF++^}P z6CwPbQ$|{};LtI~-G5nso5xw43uq;e?(-oxeY#1x zn!TKHVXYJjvvm&Ga8jLo$95k_~J0)XU0h=sc`Yfj)f;iy82mOKL9s=#hvq$ zO}h194iDDLhq_eh*;XFKrqi=m%XH|h-)CB4BTLNi-|j~?Q}%3iD3(6Xq!RiRjPXuk z6jhNlwjv(Yv@tim+5;5SzmY?#9RL2XxBlrA8W(7Ul2dG6uSu=ur1(@=Pt& z7iNdHKg86YKl<3<3lJJvN~>hn{rofEC=g1=0T4|6(k9Gczd*_H#M?@C3)98E;mgDZ z?&;8tZyC$>@=v~fhc3tf(Z^7N-JCg}Z&V@uJ?5#e6JjRL)1BOyPOK55zjyaTv3Dcd zXbjzGA?w;Y>gv(^X9+)PIe6ZmnC=NRr3}3+e9O|BXT$J5a{=AMt}a`-JqaaWT*q0Y@1uoax(NqN#voi?6OwcVgR6Fo%HY{sF8`ULPB?Lc9r#6M{&8r$# z)%ZJxP9tnG#g{fK(1idaN0QqH-^ZbwxuTmBxfVXK`Cjj$PU-G6QVY>vA@B=dNpZ9K zu8%rj$Ddi(k1Yf+GCm1?1w*DI%KSp+>)mzCz=?KE4es zt0u~x7$6=!U-bEWd+RalOe}V`GwGBzYE+QNjs47Tk*A#i>b5Rb?uDarjof230@j-* zT1$O3L_k!D$@m+EnGFe{EVZT!*h@BguZdo+Yrkg%X8gnMrJtS2m@zjdx%iFTYaM9} z@P@qCp$a`~W(dA8ZC6sd7B?IJ_O$`itvWifyrQq#9=fipqraMup#?oHIyRx0fx0il zr-K;9>(U}L3@pCb6a8ZJo85v52kpy>@C{Z$@Hc!YrL>0fUUq!6DqwgrMiJ~O^)xJC zBsg10u%JV?T%*`m?}jSo2@yEN?t~$q@%X2QfBJo zKUTR(+q{|;ZQ|SIAUzV}i*0 zPwJ+9F0Tr`wv|F@Q9Vuk#FQ4~*>Mk57KYDDF8vdfzUFG=vIssLb4&0qV<8s8mxM%1 zY!1*1j~X-~Mni%O+G|;;CGn!(JX>4bC08abRmsc;ne$RidK1U#J}P7JHd}Ma*?9?g zDO4Blji>T-;3fC1wc!A{V(OSU(|`rDM{Pdl?z4*H%<;v~Ws?SBYU zvmRSPMy0P`9ooEiSwkCv@9W(rxSG!8_6ibLuQ@$nN5~&^g2>TCm3Mhpxoe_6gJ(vS z@XZ3Zq%c*OKh=_vPw^Y|>tx{CSb0vS?&XK8^2C(Bb9!U>QtX5}0VZm*h9#;Gt)5O^ zl8A9wli@$)IF$QL7}JT_eP*jDRV2o_D=A)dvWwFhhZ9cEzT~P$tq1sVs@uB!zCL_h z_g#Jp+z9gtz;awNV|dIGwc2?cFBR5=xl^86C2jv6E4m)rIM}=1+K1@ ziIl~OEB4#Y4k6!@G!4_sqw5BrOZd(8(NVh~QRUJ1g0afA;YStOZtVbc?&?pXelOFP zV~hw9PY53$6rVntjj%=xNp!fT zU<=b&iE1`f3i3uiQ*x4F07p zTo#+#QOZ~BX5~iCYWh7mdzko_rSwVpnLtZD2T_W;ks+Y3reKX`vn;8qpcA|8CNlTh zM=L*^bAaVBZN-I%Z4hw*hS;84ug+(`C0W}38;>9$)~_-b{D%#N-54wsd9hdyqck$O z1&B+DbK4r4_ZYZtYKmI!@)l_juVez&fEmVQB@j^^AE$*~*uyQ$Vya7NuP6Meuc8eGrLgYkVGnaLgd^8^QrsK| z)fbOeCT})0_~yv2c9mbp@)>GVhi;k?{v6`xR{74$IxgYtYtsa7#VYbogPLRRBR5H+ z^z1eHm~%o_kBzI|zr~0-Gt>5J(kh9O9sMYA!xyiwPf=O_8qRR>F1^pL(2%wvJA9mB zVe@OzLAyK|(M<*x-Hb_c9%CUEgZ(9I6ch$7{!moDP|vc^rmi2OFNCJAvfVk}!1V_4 z(va6Su1CNl8I5rzRKDw+WmWj2e)?)fdM?`g;#>!P8~@SMCK*OBE4_J^CRrvH@ir|R z^EWMSGq52l@0U^ij`NgLA`*LS&Wd%#yzJ@P0Tz~g(Ld?C%67n&QEcOWl_6K09rzC- zgH&QnD;o4)!rb%nFXMRyC`#mZ1fJ|6hFuj)qjDYK6&yH0-OoB?)Z+E^WR*%zssc%7 zl(cyD3M!uCIuLe*$7E`9&h0MBJ0vNiiCi1KmiQ5B(;_w_s87nL6xU~2_#pw+{7H5q z&G(-ACz>htL4h1(9ln{klCd*fn#7R+}LAE(U%XJ>{EmcSw3sPKFZ`_JBc}HvKgUO^ysVvB?54c8bK5(~Ks$ zT|LTF@otgW^P;f)stxd5u`m)C1QMBT;u+>L@3@B-j(QZQ!1*bW)XS{FEEWsW*) zSgDXkw%-2R&FAOfAKmgJ$itz3K0^K^82jv{?}DNYS2T@OxU+DSggp}@rc810P}l6A z(ap`<(mV#D4g)!1}m?G0YG;14U6G*9aLC4~e}jnSi&W zx6ny*8i(7cD%9IEuV*sx~gAnX`l3F4&S+7E}_^;2IF=kg7BL*v}>=eF}J3t zjE~o8Nu#E45>!R>JEW+IuUjAX%N@y9!!+%`UfxD~=L3e#P!mJA|?jGcc1Ekja&SRuShM3JoffE%6YrJ&6XM?G`e+O5kgS7P(*rkm!jPC0%SJ8DmE zKc%pAT9x7bSEY}~w#e$_V&;-0guRrPEPs~-TNzm=?d$RWVvBZ5S9ws~FJcF$|KsU+ zif6%-hB3c+_Uq;11i0a^js^DmB zotGZ^2O##X(a0uG2kmeD$sC;EI^;PoB&0OM#tB%ZOPT(|n@^1&d1XwMETk}Uuy`^) z$hV~alM&&*vFd=n-p2Nkmx@J=cIC2As$LHujqMQv;m_ch+iv+npgGTuq+S86aT z9#vfn-Ir;j!xQH+PchkF*shFEnI+1Z9#JYVMu{l>_8DF`XW2GIAu*0y7PrL{`O+_6 z_;}~nD{|um;}%G_aIP;_PwNeI{i)j336NpF?Iz&^2E4LW0l`ktsA^_^wxE=!xUHzq z#FkihVj1P9>-#%0;IRp7leynEXO?!BrQMXTx{`^<#zQJXUeM=Io~0bcZic3X49qLb zptj!z)IeKHCaRHsN&vM%^{$nU;Uq)zn^&(ir9DXX*W{(^<8p*A(&Vu=3n=yrQmB2^ zF@S1PreY!4W*u!1z=QI+lO_!kfIiIrlgh~Pe3dfwn?)hj3^bV@(fu_COFGlpL6?yAR72ExZ zEz**V5V@vxJ($D!<@bZ5&-|10ey3$?k80LM;JzmpYKUCXE7U({rx^tw3it!Muge1Dh__^cuA8xK`| zf6@YpvfR=h&8HTN`S(B1g0@h<`L)PpkU!GR;LV?-HTLz3j<%7GP@f7fKj$DOv8X|O zk{H|D;$A>Zw>i|NEfOfXF1CN_nJ6d3T~)zc9}#eWxSQyrfquy|2qu-vi#f^>v5(9M zD`}i;=6i8$C}Q_vM*YE0ujMxc*~up8`(M^qEWWq$$tkmEB1-NfDJ~$!;ZQ15NN-EJtT>@YB-wA^6j(4e6}^ax6JJvV_zgUUgK5 z!{_vRS~bGv+M5}HLWdvIVIArw?q;rxCbYocYC{H&FoHATW`o`GKM( zj)8Aa1h$P8opm+lB)rn{_MVW3ilY&wWcj?`;hc8WUudUifM?HJo+Bv6dKi{J6;H113)e_Zuc}3F8Q^?hpH9(| zNa7B<6FJi7TW`3TtgI>=Vu{T=7ch_7F>B29g);R_>WZkr>OR)HENqpjg$UNL zJUJSHHv_p8YH}s1Pcn&OAyNw?RR&}vfe#lI6Uz+F7#?N27?mizYjna(lG6$mV16*^ zmzkSY5PCW(o9{vmGm0((zYbnt!Eu;MRw)t@h>@C83vG^){qAouyagvX6E;(^-c+eD zkVWp7!YX3RoG+qrB_eUC%G-3eNJlunUn^@)SW)2_b{cTNT=LCT(}f|> zQG|U%(FBeIxtgyVXI1eFH3T>#IKqP)FW7ye=R^{~2;l6%^6Xa`e-|qTk&>+u21y#H?r^E2cd&UQv0g@pbhT~NQFSi3qD-PuUcwwOA z1phhz9S+d577AFOXRO7BEi~lhQelLymkKHhqUMW*o$4g!1LUeqUb1$2z@J`ndVOJCmcQQnJGzt^qG+6hK|Jb zZDi&p>2>4{^crWQA0%7w;k#sT*0f4Tx&F2PTLrmx)%RCNa?Eqt&R>K>TN=mvw^HSB zzBYZCkFmrm7o>ysj#J}3v+cnM-$gqkz0ev_nRd*)ZlL#fa!o@;(3Tno%r^ec>EFzK z_(k!Bj|xx4Xoc6A6;B0OhjKfa@>WyEC6)3Yw_{|Nns5qC4%yl&^_su%_JpRYLsSgD zsZmdAKGt(S2v5DB>r)qHU-2{x_#78yA{3y0=%{e4kpa&$C=^-!5~uqO5tJ1gQj>Sg z;*Q#AX0t{_)RRVzh#}-R(>*bGY)tMIy-E8sNDj)s5rQJc#Pc-YjOM5PC8iqj%iqa= z9-PZiB09t+v~pc?ZJBZ3h_d)Q>uE?64m7-?o?fV0NVlgu;JRv6a@6NDOk$#WXW5G8 zgmKu?V7THdj=M7{9}>zV&m#CO#$1VV#ZFM_WtZ7^tMG!KKM-SX8JJQw26%>_<3~d5 zM(pJWeq&z-l4;SC=GxFAUbHtdVYK#Lz&?@H4$pYYdu2^>=IMg;M?;TPLC(mH>qS{& zqU~~O4F`+gtr}@Q5GoRhWDR-Ui>qU%_ z)wmv-`3^Gw>Oyetmml_iuyQo{?F5Amd>x&LQ<6m;GZ~GJV2rW6I@)k$8S$Jw=PjTS zz{DW_)K-{GiV@zPPjB*79?D-ELP`42;0vR8So6@@+qfJm2h1aOLYF6m{=S5HVZ03a zPmdwYGm~WX0n2J<@+a#vO$`Pr$RjkqKH?Io3{#gQ$NG+_)rfV|>(|=*U?DQ_CZ5G7 zlVsV$S_t&9Z2yN?NcG|H&|^qzFcOYyIk>DalFt3lE1E5$+nPGJJiDGb~2F&5-o z$Z;WdE-pui)4s7dr?YO_ah6=HyiVii28!O=Ob>*}BBC7uz1OtQHO%2DrBnSOWIyM# zKjiv#&wIRAcleq*?CRz$CJBiB~DP_@prD4)~nCq$DKUUSv|bDE~+fQ#@cOTXp&WZ+xw&S zOk$X$5hW+-ufg3$m?wBe*`8r9jWCV$jUjH)U|hr^inov#kDE3+?G2sKT!;TovY$!$fJ4Iih|~ z6!M%lp4V$!`!@`OxmP_L$6N=LlIQDE+MQ{(2W1Oe*-s;({J7X)BJ~4PK>iw4U&xN% zdD5&KB>g#Kxof2P`!S6Ks&`?g?;mS=KR}2k6CbNB)w9-me&ND<&%#{^Mf=lZBKD~f z?aO9Ylos9GZ&(|)0eQYGCk`_mnecOpG6uSU0$&^39If%CB! z3o92SfQVg=RqJfEdTgvp-5z9IiY8`9={sPoOk~hxr|d~~l`Pv%RX4Y=r^iS3BOHPH z`+{0C<NpMSGTc0I~@wWm9wy@aBbu2(kkj!js9&9}tL?~C

{pY{aI@L1LZiaZ?U#MTS^cBn9uC$W8UP}`W5lTLUGVs)ZAtUt9ZmcggtkcOgd?cEir>`XuP%}rqc z7T-wjZ84uo`Upo7zt!0i)luq>BQWb6-}}^2WSX2`SvN>ffCG$%7Avx%VqT%+SEdlX zvDg&jRb6z{Z`3jYRjKRWW%sb)SgSC(_Bmpc`xE8S7c=kO9z|3WVIgMzO-_Laaquf( zRH0+^wA8sk8F+0h&n(BanK(JvLP)9w6QXlxor?aF>#jZgsV*Y4A&MdrNBKwzn`CQD zXE}VG+Hz^rZ}u>C6|{DpYWsA%O%rQgs(G4L@(XFm#-ZA_FWJK)&z?!GI<0ivDZ6dD z?!e@5_2%eYzVS})esxJQ^6N1H#8=f2>R{sRd-^PelJgK2B~M|=aEVXsqt`@4>nX*z$7`=IT(q`3+n zD~N@*emp%CV1H$6ou5h5{-SznQ~oVEd8kee4&Z!LhmoOc+VS5E{jHM>SA1MR!^6(7djl-EPRUq=m*EI?ig&Y#6PhG_-QG6ZP6R|>St zC)|UXu)EsaRBj?Nx*C)N9~D=Te2M4(^k;^Bi6XjChep5fIjik3*-SJlu84r8a0FD7 zYE6kgKj}z_FRaFW>g-p)$uz?L>r9Ac!@e``wJxpgpg4}Fk!OOercp9aoB^Iufvd>S zQ{kNrnfI}>Cz>m*ryWjhYdAjwh#06~)L=2Yb}M+0Sz9ONkkOKMY#xjx3d$edhn{Y%L(RgrByoBIv_o zVnju%ir?Bzo|aX#f!W5QZl8Vr>`Kk>^tiilh#ZkK#%EjC1)XSKdmqJN9>LiUa{l?y zYg{hONz%-m2C6|pgU2bkP;ms1Fe4iO5bM<%E0a=Ix_nqDa4w7{dk7&9D5m$WK9R-V z(QPiG<#zq-{CN#S=tLQT&I=utmiJ=XxX!jreQQNU^^F@w=4TrJz;oOGe7=@ts(GJB ztUXwA1HoiF4isK6iqNRHkz3hq=*wvel$V!xsJsgT6J%gnI0eS-vr4 z;^_dk8}O%=S2^}_6J~fvID}t9qRA^0#T`}Hvg_3ew)>Dz31RNdxwb$J!R^Ln=&H-P ze5Aw{k+)@WoAc4!M0s;qv|NgE3 z01N}c|1JZ2Se}41mS{rTo(Yq;P({^0wDio^k4b`$bnFS{}TRFCIAx@ysH=q z%pGI|00Lk^WH$Wu`ZvA;04N-O&p;pyde1p10&?Fw;_4kSfz+vzV5^hwO1XS=YILIpmVRs^hl!5M7 z0tkq#gqVNK0v5Qh7Yx0p7g=I=9T5Zx3cxc^c~D_mP{@5az>s?$!C``TQypm^^uCu+ zXa*4gl>`-;cp0ye3XtD2K<=9f2i`Xm4!Un91Om^X;YH4ZEYo}NzyR=jM1j1+d#4D1 zf*_b6JR{{Fvk3e%i~f`PKtY%Qk#Wz!<@?K){}Cz}=&#+>f6rOiKeOd; z$pQvK;Ta*wS-4d* zZaO3Vy^9QTxex^2F;D=x`Uu`JMgV{$iaTLKjuRC4CxVFof*&q$&l5qo0Q{cS5WwA~ z;VphlBo(AQFEh{2=HZ3z0JLJ=PF}AbI%?gh(a-^5j2`2q1~>f6RB^ z5jaE;ad$#NII?57qeKvbTrlp)kV*vaz>iE0$o9BN3zY| z2O?2`-OFHP^K!Rm`AaeQu7SudNASLL@H<8$kjp{l?-~OFAXkBV|PTD6Znh1|Iq>dr*x3QyR9G)*$4vf#uF*KuNVo#J;e~{-QEg_ zWc2?K9O9q7`(MI)Wd{er?=l<^DZ7`|aIoOL&H*XAw-Cd@uzN|!g%sYaZ8!w>@A~}z z@JD$`~)a)0=b`55IrD_qWTh56W#? zD|`9!AU6}NFqhvQj(QWu1ai6^`DNa@zd;;eX@i=&DAU z&FBa#KYRiU$uv`x(vm;(FlFhuqmgmbu+5{AdFWqZ6&YRBm+!-0CO8o<-^D$ufE)Y6 z!EPc39WUR3ubow1U;cmja5S(-;pF98hZr{NmAib1Kh(0!VuHW?6V9SxKFb`y17Kv7 z^#EHtHT~;owAyfO0;$grM+sno)i0m5pFT?8hX|+m{-8Hwe58>4;h^AlaCb0T`2geV zJ!nTplx39Fbr}`(<@0toPdk|VpevlxZdS(2>%0v(`xo$Gben%O!w+#DWn@tiK%&!b>U=RBQ$SY1W6 z2>`+){2p!~Sh7z7q&|r3E+5#MST4XsFX@(_H^V;JO=4{qn{zP%E?qb+* z9si`^*H%AcNmpC3Ijw58;<^KA`qKCLh8E?~xR=1ynMMpFXW_ zrIK_Z@=J20hX*O(b`lcQ&3^`9fW60zb?f+ic}CL91IhqK7AiJ>fEVqI^%l-zfYl=2 z!&;kEz)OE|f_ER@`{>>`7uqH90{aOUAV?xm+l2?(!+$|!z5V{Z2k(Dy@%^_iKfm|j ze=g*`2fxSvd;7!R!awi7bMY_t-hKb#qg_ch?y*WTtsT%uS_k%%B^e+#99m1HuSc-? z01`tW#sE-YTj?H+^P92dv&=s5hf&+a$Y8vY08M}RDex?9qC{2#hX+RK%S}XU1D*d_ zLLMiC+JT9YB5dyNVPYh3BkLHLfPjv5&mf=-Ie$ASBma&XW&$2{ag-7^E{ZgM3y@Sv zp$y?-{jV<3o_2xvd{8Z<18FXi(P%9<0n)mu5rkc5Hs-vf;*<#4FFDEj45x=pO9d2G z*DQa7q?SbN6j|hfmiYO`A&M}Yaw~0cD9FvG?B8Pt$0hB)d@Nje+ur}bbWEk0VRW2+QX5^=w zSR}vffO>ccz=ds4r!ILK{U$0T26-!42iZNSvPdZn*Bw>R`bK`>Byj^X!uVK$d_>yr zV0O{2s^zVzX3=a}=?)Agv4JC;Dj5?0uhQK~y&?J5!7tL6N;^wrH<6xgwLTw_jURva z(J#C3?0fzKpe=k5bUSUv(}Gok&;ixElE-&23h)>%_+=I-TIwQZw2~^VWMPb^Qs(i& z0b8>%jvSj2H!4U(n+xYb%ezBygQ_s-H@FTkotukb52~R~47h!N)&Kk{NmC0K-{G{U zYZ#g2>RlgH8w+yP4+s6B=?B+~CSZSnYyS~W40z=jFm$O4-=it)1CTzLL?;0Ps^D3G z0NsLpd-nh_peV-ZP68yU87BZj)DUYbKv3YQ4v-UslJv1mD6ktkCqpQgM_fm^tf2uT z_Q;4Qtp>Qh1=I!QkU;#!!2#2bU308h34DCrQ!r-P(;o3&z6X1D{Of#L`1(m7(|?J#8Bqp>SF;UN6qzwO7ZFYGj4?%> z8IyA$3UADqPUyi3HXd~BqW2wiBv~0Wm9n!=r-P!(4oX6e?v}2;q3r!Et1Zq|-84o` z1@@vo7J+V1CG}Lnd)Zbc!;^nc1RxTryeCQ@nG|#wBZ1Wca(#N)6y4v!tbkoXUYf4B z`~gR44a7q=u%}JiBih{KuJm?ouqzyIqvQl#S8Af06RmAI{-if{r{v37-^Z^x*A>|j zne#q^G;Dq_Yd#+UW#FRVR13&)&CjwT3FzXX9BpQGJ0ox8vzO%knjrLG@$P55!AfytCa54H@|_8=Q3 zwAyE+K26a40kA-Jbu53O$pTi>k-><;3gy7ce*GqzrVcn2ISs7sJZPg}7yPY#XxrBQ z-ah;E(1`q`+HaF5&rDi>RqRD5>V8t4XMNcPdtP6zvI2!flNM5EZ;4acqxAEo#dxYJ zEkMD5eMVeovrUa{2g<%ZKuv0p%<`Uw_w!AELp?jw~F8L5O9ptal`QYj564@d>XX(VCe$C&w2-R)n7ed>wQY6lfMyFeCTq)Bq@>@I?4rC10L5i6UVk z0mfh@s_qXWBB|<)1`J(b2;EE6_t90UPLtgMhug3L4rhNy31J*4^02pW@;(HmXjZ(u?}4eV791NkZ;*=%)0mTmuYIo*!i%sL|saFn}B6nP?2?63$0G zDan7YOSZgsKt#w}))j+CX`XJ-nmnVbWxs$t=}cHu zwc<>;TTLyeNCxuZ8FXcRxdgdho}~CU&LcY*$p$QpSuhFjBNa^N-4u8)BLgOvo^yY| zJRB%xo~r1>3D<mi|F$`5qEh%mrZ;$|=1Fn=hx&53XeKLvV&LAL|cqttoDP{=5Dg2o`& zASfq4cVqnHE}!(lC5+``VrsxiH%xy7c~qhz==5QtB6HRlGHj9_^SlXxfsO&1sHf>E z(*{ni0(Df3d=7VH^df{~JQ**dE=&?ILVbK2)lO+Pp-FbmyT$L{bHzzM+}2|Ye}-*n zXTY|qBi@6%)du5)g*`>3+n2x17b~r{l~JX?-^ag|5z^YrPnwQRKc&C^E3SVpfq9KV zHYe|;i9A8}0 z=6m*dOP~AsF6*pgQnmw7l2OxqywAEyR%NQObjt4%A!y+@R{y;n`eS6t1{Tu?=bAQn zAIJ<47+R~kcesn$$kC0Gx^)2Q@q8vlq7m!%~Q*Z3S(+!n=~~Q4G&*` z_E~nl`7oS_-Dj^mA;Oo-9rFhL<&AgJ1QC*T@~1?40*r^6K9*?^t{8uIV@Y;S&U0|- z0DY>Q`V0~WZfeLfnlm9V;9$d?JClR(#^@$1^^jRWqMV&_yLnf`Q+kZB%H*QA*c;XPj$Du4I4+_c=)dVR<;hKh5Zi$EHSy3UK+9Ca!w9}}Zp9=~u zD@@Pi%0B7#;o!;7u{nP>20@l?++c8(b87U1wvB&GjCoR4X5CniUJnq1f4KMX-xJD$ ziXcQv+!dhWTJK`2yaMc%E=d5m3mpv>bpfGC84AY2D4Q`h{Z(hm@dS$%Y9_@f1RzI{ zo-0DL8r|HD<(A;YsEN{Iqx$Wq2ZUHkz`Fe3E`s_`hiD5{(8Yf^;R4YVt1uhFN5GI_ z`ay6NTn%oNW(ij%4G0rl?H^v3ecE{^cZroFIuu{kv4VZ%69oEH8{st6hOlw^5Kpa> z@~^>Qf?cx-QmpTYJKy&R9U!ZqNTw5+hzu`9S%LBlQ&YUp*1mO(eJ7EAHcx@>Zl%Q+ z=D;S^G{ETzUC@8#jP5MDkitFQu}|>n4X&QACJ`7j&|R_)0>~TL!m=}CN4!VIj+tp% zpF&}bAkSA!3kuMMbUuukJ%lTl95bs*1P~`Wtv{(4F%{=mk}_Wl;z)G8&hxbVd*g)p zPf}T}iM6b44%?F~`Z-t+z!b97t?5X|zl|s(q!NSHmWzLRHIg&=^v45gTX86hBfEdv zXDvw)!IH8%GU!RBI4^9ukU)3OmdZK6+1N_uNS<@bp_Mr6)24HgOx~P)lE_Az}=S3P5z)V;CEuU2e2h0WE%#eXjE zt@HdQ2mODs;Pv@o@;&4nyJW9o1_rz28!djYq%h&veawFr<{GAjk}MH}SBG$~Za|R_ zmqNGTOM=HJs&~AjS!}Si7HyhNM;X{5q}uxrhQC7P^+Mo zZzf7VIg^2=><);(>B}#fDC71QUtcwg@s7N$1#^E{6yXN;3jQ*tn9*xethS!vOFK>Fa7?q&Z?+#PGUWN6CQ zU|9mbClBS0M>XLQd_5RPC^qu~VC5Z2kb$^IlvEP)c^h&62)~8z3ZhJin2<9{g;rUm%i^xkiVL6ttA~@y(r8W)EKtZ6SaC__&}`2}3ant4l&^t?1S!gML3Akyr+0l?DBuK+GOd( zU4-5R$?=Y_^G(NcbNHl~qMc7qAi=DbVV@#wGiFT1t7WPfad_H%cZW*AckryEeO0TD*P!dO%1fb+(Og!pzey((U&oI{pq~#J=K>{#- zaXm44h7KzwJ#{wJOVodZUGa#U;7QX%I$xn1L|1pK=$CQ^0pw2xmL zxG*+b#xo^i0qY;;q-H78M6ET*k~@U^gi%q)Y&zW2d}ppvD@#qmF~aZm#s*xCV3{@wvoqJ0R{&e8?oF0wGx z=oqQ=O2S%f`t%FVA4WiUzf>={#fMh$44cA8T7NfJGb4IQwQH$j7c=Y&xFaI#{JhG# zXg0U9qi27b>xCOcL}qYub|#`ZRyTHrb6cw`HAIqAMjO^Sq5ty76DdN8iU%bdZ%H=E zImlvGktOYzIM48q{-5R+a|u0si20j)1YsQL!+Ixt^=#&wV~mxJY0aCXmJAa9mP?{p zWw}sjkougOi3Ej_T0u1ez91lH@1W`Lda#Dn29$p**#Qtq)=oEW*S4ia4NN;%yMd?^ zp(LN`j^f>Tb+j5%3t!P8T#!q10x+9_SWK@;LF=AAnJOII4V=1|V2w5XH}LAFA{6O& zXmKF$TnZv0Z}zd>M%OMF+U?p~X!^UjtLUZE-&NvmDRN=^X$>_!sR{*HA930A^Yl^> zGW~y+8=2ajY(FB35X8}@0*yDwF|yFZ%33>8VLJLAG}Au!&FMvJSlYbJNX<%DR;R=h zV<$iWD5(fxvn_fvFusA=r2Qu6%=+ z&>jId*=7qa?6v3Nn|G5>Y3)3%FgRAoH$j-gO;~m^1h7(ygbllLvZP(i4LK6w&O(2D z5{~_fd#H;BR!4G5GtOoM{;*}>`9;SX9iw|imnGMGj@cA>lp$c&o+CK?U$MN$gMB_b z*syXy(DVW}YC%#D&+BzH<@+p+t=oLw-PEU(LB1v)#D(kiwD&GkPXR?8b1h9|}h?>khy@QuP@xOXusn z>*YQ!j!y~UtuyB((9q{tV)PTy@wd!#Ck5;TBA zp4kPPdrX9I*(JJtNYl4LHO^S)b^`#2zOr6s63y4m=La9wD4t#f?vmGoonWoxk|KT) zKZ(YMdet}CbjL2C2;#CmU26h!Ftd7?3Ymux-jd%v)mSbWBIC;^qZ}$#WMc#|*Q+s| zX2$(<%$#kZu)hD|xni#q5-opa|MTVCFUVg^XPldd+;7?t_gD4~h5>0WfSr>UEIyyg zpi}noUtkq{BEZnie9qNI=LhxRytJN%6!DZ#ZDQ_q`Z>>degng#?pME>iWw^{J{3s6 z@9z!pPeJ7e8q!-)qyawJ)<0pcfJAKRv~;p>517^N^f2tJsTseAi=%%}8sN0PCeyOh ze#sxQ^3?s^;_b_K}bz=Vpt zFN_73@c{$nPbbLOKv9EKurrzk)*ko5^g{63Boo16)G+Lt2o{5qk%lG+`y=7ry~g2d z+R(`9@eK#XguT;?o0ET1nId|8CC;3CFf!ur@R->8RiAC=-<$I7xr8eRB@L4}Z}rYT zx(_Hp-prc5m1x8!w~2{>EmFQ(tJ|_+Howmy?{s5PL3)QW6%{0bev6_4)Zt21LBfvZ z%b!_RupN#O?C{|Au0P)jOkN3l_NK#}xH=N~<7b!b2z6 z`t&MXMh@5$m16}m$if~jzUk!KTF|&WBU5LnZ>nJW5ov$nww$UPEq&YNk<%TiIP4pd zddX_Q>EF$Tl~QSPWt^*SRE=Kf8=H){ql;zj=+jTAucu|Mr?~X1P1+)@oK-3Zcf;$K!=oznobxbKAcfbSWGbIeYKXR_66MpBgP!t(Y@P@#291y zG^7|CsLNgr-Ap@7a4Xnf3HJBS* zN06<2l{*K|P+_-Su5T40Eb8U~hS5cOQaVMM@ z>0XwsDdx#IC8|yZKA^2ek-JPhi5-k^wW~Q6Hb!?dF<^pDJ z667BkM8;Ak{z*j+Eaf_xj!Y} zi4X=mjly(Cf{*s7O9`xUzZH)|eN=lM4Zg?kQB{kpjKlG_Fo+ixdqT_+2dz`ss%Of7 z#zSEzG1~pV-XYhDc{f`kOA7mHKHYQ+d-3!-qBkuaGcKGGo8(XRsQjs?gu>LPSt?pM zIv;jv8d$e+vK1iQu(dtY;jMo|GGSZ_HhF`_oFqosRX2122*68#n@4Kil3qgmTX@e? z(>p61ne~1%3rmG5SO9wI()lS`iAsY?J;yY?rY7Va#x&6lXoUau#59wDkse=#x*Xoa z$y=3N=fFA~sC=^Q%lFexG-hT6x+noQ@q84AKs1J8L2gp562~T{z14rIF+Dldh*}n_ zTfmUv{V{*Z^Z=M7r6ZZ{^;iYeE|uQa^%u`fFdV40Roqr358t(ACGL+jdE9MWX#khW z@bU(HK2a0z^%2U3u2Uz=atnN`@l>&ec`ofj=#v^`rk_&FOJQ$;itVzI(&<~ zLr=g!cWSf+MVU70233Co(1BLp6M~5EqnTldk7iOLxbB}mXEAzfzvY_r1mN2JQIPWi z##G;H*N0(lCcjd%fWc!swmV!Z0+u%WlfqD|<4-}MK5=kBesbxj>x{&1JuUXuFg-v0sH3y<%Y-BbY^12`}-m-|!! zM}Hkfa(P6X4`{qndQ?5=SHEC$(rFn$`)M#=;^MiuFA`o*Z%HGZ1PIrzw-a5FCTsQxfJ^9`I|@L>d)}s^GDx3lIfPK^1r7q zum12DZi`oOi%vOn_4w;Y`EhyWwMsujjQYPH z|0|54Jsrfzuo)d;pvTYPxzA2FZWAK-aQ($eZYEe^u0KBAh)o$Gt<6IW%ec`oY*t;} zD4EB2mL4=lU*EQ=uiv)Op7)F2oRr%c@v_9O@DrY`z|jTg}leA?rR4E?l`yhvDx?cUo;eN}j<8 zFYhP{x$>4!%#D_6Gs5l(3Ljrnl+&iJ&cG)8glu7gwAzw6JJW9Inu%U!S*7eCtR!p> z3LP6m-k}e13*fk3p0wD4bAKKd0E92#U=;3;*`zQo!g=W}SP3sxpLlw*o02gRN*MxP zM>%f-{ds#bz_>!9WJv={lL=C7!8kg_5AYIjx-qm7)1)SbM63!nO~Mm@!9N0bMq5Ty z+9i&)oA(yh^A2Wou)+TNk9eyJF!zrq1Ft7spd;O#wA-8lsg=(({2VRqOb>YSD%yeCMb_!LqJ7=&sZ4; z=?NF%%r$+XuVxJ~O6mHHiL6G$_jPWQHn0oArvofwkRarJZpn~&!AbJgi%+{i&3oh{ zj-Hi{1v=AJ>Es=bHnbH5^?teyk58%VFR$B#J0cSeXAn&!x_|cS>+Rc?6b(^NU!Kaw zDnBTlb}y2G3y9V#DDbdhWHRy-H!pY?%PchD)BclR4*QPnV$nW%$D8)K?I^TjIp0+z zRbyF~f+!hvY2^fzm&R}2SsuPE7%ia?KWMCeru~`M?OOh;4v8)`DAOq7_2+OAz-MeQ z#K-}9huxK^x)ZF;IO$Y^c7mPa2tX?9gZ?-jCS=U0BAaHAqT zAIy>=(|<+&|9Rh)zhnbv^$6@o?;(sqH_Y`Pl^p&7Q&=E|W*+XGq@VI|z zgMtq?MTgA{P?dvO_>2WA=#XiW9i4r;o>S zP{XGs!#hUcoR8q5vpU=TbxnOG`56U6wj1-28cB~7R5GYtV$=eAn200=FAVTGxx(R{ z3>Yj3?{IYyeIhUMV6;jehWgOQCp5JTfVPkDjew4X1o3)Kp(SJ6( zp(F0|>V|U4$i~l#=T(MDn@Y@Ahs*c}yMEYaLU7XW@ zROq-{iejHu#Lp{; zS*w5buy?;0OyV65eQcK`&0IFpv^NS{2-|Ws}+U1>Vyv3@Z zmC1#QR52l%%RK7%w6P11I-N4F!$xd$?l{D{nDW4@0BOcNh4r;2pf8 z3p?w$%9!%TSeNY2$7_a8`-ixG7D4e@}$9>O!;jn2`; zbkbm$%ujaS+ZMMT+-r3gdVkQMdkHJ;Oefy>c4u64-ABV(k~&q0(`0Q@jF#mx;-+B9 zM(qGFN>t9G=p&N%9${%Q$wto|ZZr(pl{z7V<>?^tX@9^eqnjcTaiG5F{uabzO1Aab zx+^z6A)Amingt`tx++A*zz5nDlC>7_Mr3VC;XzMxc-J8Tac!}-A%Ff|I)KG{gOEGE z7^y+nOh48tz^DSgTMC3aq5cn;29IHh0-*l@YuuA{H)g#>hb@3H>2Mqfy$>)C*`M(Mz}zQ~DWT8eDzm|jVyyk_r2acD zFk6?m@G=j9N*#fUBYz)sB~w$T zImckUlV}Kv_MDd+CoydQ@B-jGDvLI;#{Cq3VUuxBj1;xw0LhSAx9{cD&75Xw2xfBcB8xb?)8;Gp?hcG01rK1%hydEd5GY*5qW8oMCj2qy zxZpGO9IcB#bA1`X2IHcN;s7FPQ##{)p z2Jg_*QQedv0p}y0Ph|w%liJD;|AV<>|Aqi5_&DH_7rUPzgMSZaB7xU6^Jo0|*-6Ju z!(!`#q-t*fn?x}Z2X^v9Hu%A1O0vct3L(|3TJJ))5q}NqGTbF9PnOFh)QdM|6{kBg zBTYyyqX#h~-&=^SB2&qZe1FMR?HP0j7>dL2*tOim0HvNV@am44X3Eb+ap&zamWZza z!f*)Y=CHZgv~8|SPGYRg(Vm#LkBE{xX|@_3L=dsYih(3td3CsfHkG~HD) z$P##S0|@##Z=?)EwL#~N*OSI8RGFOLfg*;QV(O7vCxUG%lFVSo1aAJ|g6gsm}wb4Y~CzzUR$Y)m?5 z7Eqe6kyMW~zY-xiiAF5aW!E8wv`c!FK`AZDnj_nz*`fP*^ppcyb9_e2XqJE-J>{(# zZIfPuQ%%LEU0CrOqXh_uW;&;~Jkp7lAqiuf9z0!-$^N1vP_20HbP=aY@|# zL#Cs_aGnjnihcppH+Xi%AxO2-)Iy)3M=ZtR7Lb#M514-@ORQ{VG&Cf1+v#e?)J<~fab zfjXSY&%UN+LOpHfKx1|TL^IgUW{ zv89f!qRfQNe(oE@A?F!QWg{);MS{4HZjaE`^{!cbpJQgs$HhT^TY;*FOC}|*65>EfxxP5=!`9j5)vr9<}VQqS0DuxYXmQ9_87VS)| z)5{)`^LwmJ!GvW+b6UtlXKePW^^(y}1MWDEOSdG#;a(*!%c#A$$n#-#7;KT&{n82?F)bNRuZBjT$Tb$<=f?APU@ z-;H5 z`dl`4)Vp__-hEFFpssa+f?!ghwfdC6y;93$Mhg=3D=uZSi)GvdV3kEb1(-5!&X2Ub zPXOUP7R2Hq-~ek$BbaW|U4K%owsYByY$L%uFpOY3oWq7LUk}A2Mzn%3xKJ{KX zV;4qVE(;x(z{`?^aL!fs=8n%M?y3nGW)nSW|2!B!!%{iFbM%Z{%2Fw1EgbqgD^*Fq zqW-i@8FstP6yYV~BY*G-*}YJf>1t9On|arg{`1k#4NY3?DC~R(Ez_YQ5zdmS&07eg z9Hdzvd5B(Tro20cAexgnYa%XNw-uvZxlRC(j$)C}s4f}1`)X&Sq#SB$Lo@FkGl6J_c>^Y+3 znYNhzFiA(qi^oXoC#ggV1e70LWi?x&zcE*X1;>FBi6qq3daG9_2jnQz3)Ld4Y_CTQ z{5(COF=mqZ;eX1t{sUt%4g)Y3mUD?@YvN2^bB6m3ym(sb->Y;BO}W67N=Yquz<<~1 zriMfp?K^TAjSeBN#E5~!fiVoCxv$>F>^o@r^*`6TVd<>Iw0Hr&V2TrX#c1pTzz>#P@h*Jfa1 zCt@YswTRuwFPWkXpvQ~W zee(UuVpguU^x|S7qGs;;_^cucHZA->kq{ys6HI^}fn=B@f~n0To86Wk9(q4<#lCo> z2a`ib;&~pks126HI!~85xROW>Dnejs{jco+!hba0m`Wtss-|6z!A`gu3xuoXz3OXg zbC|?+Tu*rDtVEy?l=2!9j_Dac&!G$>wa-4pW_OKV(u0AHC(Ox3wvt_Yip%cuzB zsUl9_WlRVE?kkN=Irhm>IvLmJXY7a)Eq0+dwURd^mh!5Zy;j=GvQ~*vd|bp-Tx@5G zm}2CLn5y?E6+n;HYH?vKjpAgH4~0jg1An0cFD(G4TgLf6wsRvaSDJchTEpBbdMpiV z?w}y1oF10YE2B06XobK$)~&m!?qys|2V(fx)eGZNT(1XizS;Aj?baH|F zU28R1GF$5p)*dOm{{h!8j7{nqWF&|KUn=`U!7xB|3Xu)qu^56O7hj)J2-~97P=BE1 z_a_4~9IJcMG3Lp0bqQ}@LZV{5-{W}neafVR(3LQGardKpk6)?pp< z;LbrCDb&uX69gxF6aXXuqI==9PG=?B6@vk25$Q!-oY|DT0{1cmu@cp zsKo}|!(EJXLPd2U_D~CwZj|@vma)h=s>^$sw_S8lHY-eYs(Q_$#2ls5@=VHly?Qr9 zh%k*rjzX|@p!b$ZARj}qpcRnRE^vFN&K7#Qh<>D1_UZs0YX`K$ZuC6&^nV>=+t!TE ztpv`dUNxng?s4O`i#zbN_Vb<=avW4r%ygkzKD`?>CF=ws)Ij)VJE{hCQzbkOQgE#GHLSJTEQ{h=9ri_QJYr0TaMa^5GPo zsae~Si-Pt!Nr&T0QW46=Hs*WXG$5Zp&elt<;+1=@;}9Uqd%Z(|#)JaTV(JougGHg# zg&5Gq83S6OcNlP-WW<1JU>*Z1j4bbb-)~BnePizt*c>HnYG5CxCV#0|=@#@V#+f1H zP-w}1?S#vdTBkG9%(v0H9=H^j<~gbPK>jZ5912 zibv#HI`T4U7Je_nvVUpYim2H%j4Z+NQyhDZDd)ghqH9k?58qM_l#IMY$ow|U%s@oU zL7x6?N5y@!2^+s^St~!`Sy}8`!R6#X>uDo>#>K}w>goMGgpQUW_2@Vwbnf54zbfBW zdNkZXK52B2q`~y7V5JTV@>#o;LX<3D6m^%MPuzf4MdbM#3xCA$Oq#o}WL zT=iEji1y@Eihq;c2$z6t;m8JT-h%-3GYWGiY_D>zMjugSPfxW`5+}kYM*}X|UG$NX zVpNe0J9K%L20?SubSaD&!y`4lHFe>k(vEye;=p$ME2ATx^rwl!{5eSx40YK4Q58NJ zgv%hkpz$MTWGO(yc7JU|)PKePd^Xc}NyVreGNPn| z(xSuIa{A1Uayv2@Z{++3?Q1{l4AZg*4yhG>j$?6e!&ls=55J{|VdwU-%WiRlmk*P2 zv|IGVJx+XHmNMnCzZ)z4CiS9fyUsLl(0y9qrmoI<(E>L+1knShi5RbKbJ_Y!@)mV% z?y_x6>3`*O?NI$jZ}g#0FvVWSey~-4?9w^90vHS7iA#gy=TZHnqG0BG8|yF}F?(9C z=P6#4E}C^;=MrH)?#+yqE5Y&^h8I}gDH=VNAAb_JG9P#6AbzfNMw5!wy$kQ+S--&m z(QECn@vV62BgAnJSwMC-?Hl`{;Occx9}4z-_m}w#Ye?D5Hy0>0 zK#-2TU9@5$Y0q!{g^f1&k2eosv5hrA#LR($e!pQ27CB6TJ&dEd07jRP!+uS_1nte< zp@97rgP*)ZA=xW@5<=8qlYVl-f*Jh$N`LeEB+Cp;OZKbh+Jo5WejNlfRwr}fQ&KrM z`9DNYOt*rsu-f@dn(GfFjE)!R(Ne55f6HxbUO5K0N?$(0u}$N0TZ}HV7|i{>{CSg2 z@ViY`-g94aK+`!+{ZtvHT)mcT-5>A3=VEe0R%YpWc;H<-U(3h3wK!{b)f53;bAJKT zm_Q3A-vk#pU?R(YG|rP7OfSA`KwN-hjbmxd7}}J!Z=Yxv#c)vqGSYfFs`{TY0uV{N(>d>okH%`>CQM=*mhcsyuMSi4> zq)0N@EDI)NO4ya-<^j!4hqUDT2Y)aj&w=Jzy}Zm4KG%nH$n>tQYxrW`qYePSSj}U# zp?+^QUkXOY(0uuP<5LQCYl~6UTVxz1y{@Wgu%AK>=+_v*?-;Q`vVIEoPs4wXD>eAF z+JA2}`RVuzFgN|{Cos4&_yMWX@&&**eqybLxw8qQtR(F`k6Rno1=rq1ZhvIU0l$`Q zyUdMK@F&-MjgSr77-8Jz)MCZ+!7SbuyhkEuv7$jMw$Gi8A7HwqS7&k_sB)wR;QMr> z!UNFmhdUps0!yMCDYik#M+$2~cS&+oi6Owv@lyY|EW6$cYS%fA_2Xk~?4*Lq)+{xL zHArX$o4iO_GfKdA3V!}1XIt&W>O7jSOkbK+6PHeoyH8Y9?@Idag&GJqUp$lKj4$lM zPDU|V&Iuh3u2NEvcu68i*$<|lI?*VkXn>Fx>Dw|3m&5>MbGLrAC&J08h@eo<)DLUYJtO)NgE;C5_JGEBtTnO zxIkR|=bc$=R=s@PJ%eeOkQE9F9o_HMl9l=8S~*|u#3FYB|5yHf^8D(-ze;I$Prtnq zyWhZnPp`hdlIf9M`RB>=-G`6hvABy%bjkowK(D`<-J>tB^6PTvOVdDCm2I}I*5^BF*?FPk8a?-PY%~^5hD1o|F?r&Ot8Z2 zA0Dp7qKuH%<{^}&LtGoF7pr#HN-i!B^ZcL%%&FY|^+7CFNN?5tba?7ui%t0f9)T%Lbg%1(^oNCqo%G5rF6_KfK#qV1 z_u<2<2M<5qefQ1lFRmW^`%Ydx_$~h5hoAf-{PXam-QQjP>f!Fw+Z`v+YVk3`G;*=! zIKO{uJNxVaq|XAO&<6v?MNjuD%g*=4!5Q?oHKm+IMX9|c%=*!_uu0I zn}EmgtdXV<^yHw_5^M$O(!j__qhA5D2h!8moEG`Cw|cu%OFk0JuA?_c19uuVm% z4qqy}*dPrVyZmfb*ng8=SHPU~+nc;1=k0%q72ULTkhfP8;X0mhzmvs(Q^+n=eH7n% z+#kqHWg?u$pBoX&1VZ;U$NguY^|wBT1u7k5sKSe1z7LPe2m%tv6(EbM;1(H8V}CMM zK6)!9Hx>}Mh>A6_9J7EEeNib?gnYyawMZFdo$MFzmWm;|$jN1W>RS|ZVYEZGtrLGP zbwVCXv{_^ZlUqZ{CkGuDEsR52;g>-E1}gsnB$2RQz5f7S2BndzZog za~RypMMv3Ip2jvYM#+qAgZFi8CF|HYAsPM;2P8`1Fi?iL1RW*}M2^iED4eVq7^*O$ zV9XGhxV1~h1_Wv$MS?&jFA0G_W0#%_9KO)(1 zm%k2t1xS?D1mg;;qMQ&(?vwT>`{0OTpnL-;kN)MCxFV#y-U1|(Xm}PzM!Q77E*75& zppZwqIDdx=aY@{VM-#{bo_yCNlKf0)GXWNkuRVE-TM4oaNG`Cj47+OyvdhQDsL z0&XR=v;EC_`hRCVse81Zn6du4 ztS5Cz>*?IV?_xcvci4JzfHZ#>gurS&Ip-z_{G0XkE?H0Noim)&n;K3d-rm}90w&gY zKF&0l`4%BV0`SEBLzu<{Cn6kU!a`}N-9*5pzD!K15dbB8t#z9ZL}_zAlUTjM<5?31 zIM)n{0VvTXv~5B^U|pFL4r!XRz~<38aR4&y5wKL!APL?USSoAK#dLq#Tof!-SSB`T}Y_5fjT;d#R06WY&!(4C;j8&1E>obecb)eHtW><~oY>i6)D04Mr%~b*W#= zH%k3Yf1t=W2ry?yWZ^cE7w-KZ`$u9ltdT~uD>;cY49uAILTvQ=KOF-cX#iP$WGn%m z$qDc$WddJn|OLX29E7Xjtng2MlogAGe!3TkoGw_kIT`Yo(tqNox1>*17@CpG#+Yvw_YK|TRODj4)^@$PoW0tT;w2T-N-u+5_r zE+Ni|ww{05SkvsYffEg%yU+Hz`D}e9ucw3C^tezZD8J}1cyz>zdNC-1JbKrhkLiOe zm5<&;LQBBiE<<^3I;bAk3}kRJcRoxo%IUH9BmAnXE7p+3F%<%U*;2dpNu+{~ixFW^ zUC9*HjYcpJ$M7Tj)F;0@RNrT1iDt_io&-c4w)TG)RFT04&Kys%R(J=@m|1Zq&6hTm zwU}G5!lU5j$uS95jB!6CEUU4QTP3wqKg5wi_kth#18f^CIxUNLXcOFMw>ZEdm`$QS ziZok*8ab(pZ(#mlV%X4D=$T^d4JcaXgM9L#zAT{k^prNy52K?!%oX&oOR^!SBj<=PyIueHqr5BM#nH8aNP- zJh^qcPVSRyIZYnGR?+QA#`K6N!X#|pF+JN0xr2y2v#J4+=G9E&e+e_NX4kLg`0eRZ~c^$CFqJ&E-O_~h_L_W-%&0NSTVdtMXNhyVq^3=gJ*kFL9xhhqc7 z@ylp+efS1+5$K(4y5bZnN0(>RwP^b@8ibT>lzG8~%Lcv+xGp$Qj|UJ|Ec$=r&xb|Z z^3G?`_R7zkMbnDYOF!b^7{*{(DC~xZ-#PB|is4RmggX&H3wJsOK!56hyyi|nL=i8- zodJxO81D$6u9+)umpidyBdRxm(VsBa7L<{a<)r+|yA7kOqysXl-ey#=q65OH7Feva zhS3W#YPE$V2D)GgV9=gmS*L$9&Yusc_;TU_h186V%MCguvCR^9;PC{hRVKVao)JNa zP-%aS4u=f?^Riw10m?K`B#M!3qj^ zb&?Ylb6QXU(cRmY5_Lf-8GyYp=KmeApNn-KoNi7j+q1AvzeS~7Pf_3L6wg!71nbPt z2j?0HoriP4RPHu77t<^LcA}3=z_5+mWE-Q>t*>ijf+kWWT_LTuQ+!3a#gOR3ix$V9 zi%~FsLej=0v@=!)rxAZ!)dmyr$r%AINK}Y75CAw(_|W1Ce?Qh)jO9h>!gM)r(RfJ< z8d3;f#a7VUs>)2NKRgkz>A?yjJfwE9M)Ltidai&gX(GAp6X^%80O)w`a;*E!3o>Ri zGK4{sK+DyUa2omD(8TP#pP?38oGSn33PEY^DX%iQxuCHuAo73r2V9ues%1mGThVGx ztM#HG^84o{Km~??ZH&ZLl9eTCSX`H?lmXNV%X~qOg2qz(VNs_7spR)LNII>Qs|n4~ zSMbjPj(UD@D_knz5T_cadkTW8m?V*e5Xn28#cNrYadaXA5j0I1u&O<$2vi75&jhYF zhP5Sz#z{N(jFo>nDjKxj*x2g8aVmmogr!Km7bYVj!k^mz2FCOe5TDLYPO^|a z^qCYbNJj1vpofgz+7X_lHX*0_BIcRLZT9SBVs3{SS4}*-GEy%9*t#GAGYKrLU!{Vl zYZ-X3F=D~qsz@)5^i!9j%81TWl|{1St$2pUk(aU*jb(pWsYoj-l1HEEghI+(o=CxW zB&}ZOpQnx#KkAOZ~Q^ zEL;C~q9U4XkM#jw)jW?`7tBvKWh(d3z^+b>_(m8fn-D9_AiZmKx=vV@;_+*Dj8fmzg%va3`d zQ;~x?Q8=bz0u5Ivh61WwwKxp)nmKtEuWweX+`KWjQvHt81*MJl8tPD&UGE4j=H5D6b6krYss?bS^?6nbL(Ri~^x)g}ZZ~m^8 z*uR{YS>?-NJXEprHOM3w$v#y=Zb;(V2Q8zUxOOaN$9Cm~jf#-H?bN*G4Go6Vr%a-% z29NXfh6I}~@78J`lPkTb^cXIfyd2pNe^-Cz^etMVt&UWA1kWoQlu}GD2ES3~OLa-s zl^y`QMPE0UkFp&ct8`gQeyA#_7?M2XLWoWMhRcs~WhG8EpZw*$Ex&o@0w1@=&pgUo zz*D)rwF}O~0?t@7QHjiWlufBpl=+@SkhsuLhCsZF`Rvec)tu0Yr|uBmrMhT!Q8a&R zS5VV(an=m=I<2x~6u(mJS&ege^ugdhpl$99>$t5<^YqTjG@l%iS+TA;Fq`Ic*{K3V z2Fr%*sdLc3w2&=p5_IEw)R!ez^j=YoSo75q52r4bFn5Kvnh6{ULeJn+0^UKSGM;qR z#e}tb#09+HtwHL%r_d)!mYd zAfy1a6uoLs>kLf_4VOy82yAP--zf`{eJ|G?$CcJn!-yICC~+RVU72bP`%M$B;1a{$ zBtiDZm&=q!{dkWXt!%= zaCV3itZahOXa0iqD`yL_$X%M%Q*FUrK<(i6e%WVK(;dpB9{GBFV)f%d&S389R^{<^ z%mtP%VRn6MoAKD)HY|%({5@LEQj5coxk6gteS!%iYgYiEOpvLa4-y-e$zGLZTLCV3(O7?EHVWHKN1r?IWP-UBE2M+#(|<5vSps7VP%vJ#ecava zL5CUP(-lg>`sMIrpkLN*3U)IXY+&VHNy3`M6d3N8q#MfKugTt99Gu?*p^_c%(t|lc zn1&E*i|(lx1gvz#0PID9z9KJzxy?q9=dfql=y8sJaY$~M5N@!8>VtplP!xUy!E+vz zl-|7blsam}xUmal6$}RluX}3l)|3kq$7af{qqOKI_tT~>dJ`uiWHc*%=*IytMB@F(kZVz|zP#VRvdwRr8tSRx)P zR+2f%Fe2Ti8;xRWJa~U#l8_Ll66>?E%hj}M_=L!hu-gr%KPQOG2?DZ0x;dQb#njfLvWzaLYOR8ffn(Q(jiwn z$g;Xc)E!-Z{trslaovre3T?_fJ3bmWgwCn}haS7?=}Wdf)*P50EpHZOUSO&mt03J* z-0nq>77*6&SEzpnSxixFEW4Nnz!`vzNwFCp-0H*jtx)|~g^VPGh;Q9c+9$XW!h1K$leg%n&U#oO&(3kDs0EE(4hONSX3*5Nlwm6kJ5IiNz4=AY zGU^O^d z*Q#aT-59!U{Q}6h{oFU|jj#!b1dJp2&$V7Xbpv&5nt+a-K9vULtvA~16?2;#-3!@T zh4Y(7y=i}N5Z+AvOD-K8#3s5{<)B&b8SoK&a#n#!SGv?FU3WD0y=qtdH#a0PZfe5C zHzcvelMKU#NMD`fO>VGeknNAmyL&^i0cmnm*>Z{wqN!xLsQVy*Ta9c0L++P;d?Ooy z7xDc-zNtB0>48jIbkQAsgiaN94UG5OxvLQGxTJsIZ_naXV7$@nc}(LG#hdFSjxUq* zblI}OmImAO%~Db|sFq!_V(^;SN>-vuw%V6gL+*IfEU7TK?C%ln;IX%YcGfiRCJSYu z?NtX4(yEK}+$zsv-<3095F@21#uwhHG}!4b3_JaK!Qk9-o^9T-RWZK>l4t?lD|LC0 z1l)hq)(^y;K(drE5s)0jnsE<^D0X}xRAjNOh7E@3+Hog!D5_`)^>vTUQ>`#=1!Ad( z@K)j~tCS#VwF%SnQ<0AE1gv_`&$D}`gPvOSz2X-3t$G`{rS%xfN;Mc?7#%a?Hc zaH`z$5jFqZY|HxqSeyS;*Twggcqj~ex&?p00y$c8<<0<$P4cF{*>T}+iYaS#pTyjjJf&T_ObZsKdYi4kuFuqpm$HUYNY zoSE}%F_KG2G}_Rvy9`jv_IhMv^n5*e4-mFiReXQf)kTAn`U5{aWe~n*eEOgEkhBf=p;#sKj}5+bAm7ZA)Glgi&1*gDC-S7S>LH=^&2ZR&?eXg?cxee zhy{PSZ-sW@!mVf*&1KfvZ0})Eb?WH$|ccf6WkM|Prg^O>q6A65*eZ`T|WQhx{Yz8q8Yg&8U~bM!j0!PJJacxml!}2^Qm~`obyA zq`i=^-mKE>F*0fPsaony-LErIPlfVFd;WU$R1&sNul^Sw63e4yZe(+Ga%Ev{3T19& zZ(?c+F)=hCFd%PYY6>vzWXiw_vMqXp9s3;F8%ZT_1&+Y!fndkg(+BR_1)7ip2X+HT@X`$5Fy+- zt0rOZp1ywa>HX`Iw3AZO+<$ixllErT_kTS-sA<;P%lmIm4`%Xmp8fqFm*4;QG*3=T zBk}1;N~`prpZ*R$8~9hX)RR^LzWnqBd`Qrnf)M83EPuwwq~Xs!z(bP~;qF0-$;O9Y zz!%Ld9DMTg(}SGCEVRA<7Di-$bWjdo^P8tn{_y1MJEbR2(|y1{g7TAe?hZyfPtqd| zfA#T`k3RnW-M8O-_xY2L{_{>e`RHHqe}47Jzra5q|N8D9p8V$HyFUymTaR_2teb$a z621c?N|a3ll{Tc|OJxlo{=4PsaQs3FF?oOZ_yw$$F`*j&>yw({+e*TJ&RB+D#(e)N zuA*}a)?__DoON2B;K$Rm*1(=RetH%z3SUj{O&Rl_mzi|7vCy!Cg38cKIPs07jO~gFP*>t^~p@LnXOQ9 z^4C`j=K-Asu%7gaQ+F0xI+NC}(RX-2~&$>8%vW|h&gPriV+ zd!Uoy((O&Us5N`&6(2GDF$MU!46{7_{R!k(n3c%inU$wkTG$kSBEP7}QypMWeVgb7 zZz{8lJ0Zjo(8?tvuKEFybw;z?ZXG*0A~^eR2unn@NedS2F}JXPNKw%qxq(pAEO^c9>nt2I>G zP0ow^-@`+&mOk)*k^Bjskz$s#SM%dje!Q!nKqAnDA0$#BxSdqg9OPslK}F3TzS@dF zAfukdthw0cvix25*2`3Y0WCIG0oq3~b@23sk<$!lq(A-h8UQ4m2sxciS!4=8D7VN7BjYh4jLBOxz**q3^JA&-(QI8MBOTE;85!NJ`AjuW!*0yc*t ziK2AULV(aAVeLgL_Z849Xj5y+Ory^}>j?IhIFHJp!Lo4U2l)0d%|hlS#D}<~VB!t7 zk*X*GGDZOf(jNFNzSYB8A+xK=&5{J!3=hCHfl>$m2KkC2eMVv6LDp&Z9$tol z1kj}E6Ww`#v8@%Mw==dk1upa>W|-mg?BdJvC+u+e0n`8*W(pSp#04R(>Rhxz5_B5K z%NL-3|H`YUk@-B&W3}qYJQ+Y6E+<+5JmPt<_*D9*+{? zJ22nKbYDeqq$306F*Edp3kl08q%43BR*GhSDs)DPQsvp#FS@`@puI&obM)keAqQDC z`HMlKLhUWz8zia_P?t6+QA^n7fTfP`6X;oyFyD>117WTvWn4wu#!sN;Nhm{(0txy> z42^Es9S$I7Z;}iGz2sn}ya*nblV}9579fqo)(6YaOtQi9S__K;_R{~I(BnDVYCz9_ zJ9?E%Shf++#R2lk+RH|l;Kg@>CU9y2%J|Tpug@1N zu^PoW%tb8Xd=rDFpQDIJH-Hh1h<{$hdF4MK^y#_aXA4pEd}Jx*C1bOdu1F3fg)rFJ z8fFqqJc9^9u4|B%*C9guP!KCb7`bkj)9c}zZ-NnF$<#nmDF1Vls|)dQ-YA? z;p$BZ^ob@Zltu?sTHG10@lR=;_G}_&-1M_11$FT71uVQn1;vKn_1ruINfclqQ00K3 zssp*y@~>yzHW#X7>l8X34xcf38GSY5u1dcD5KUJY8nDQ%z9^7Mk1{DaTQF>Y+;;)$sek5ElZrQ4p(WQM%ghx%128Fc$L&^K%;I z_-fCkF>m<_#oAYh;&iPw84Qd-g8K^557tJD&>+@99)SV@n)!Wc+eLxO0y`_+?0k*_ zr3Mp~qD4{zmnbs=WzvAylcMW?{i*^gYgkg5oA4r7Gf?pCYI)NO{2l&}@b~|?Tr4JR z^Ck@lc5xon-y{_a3>rDvV#FeQLD6CKc@HVAA~y}J9T?TgMHYQv7WR=~KCtx8*=TY& zD2(2TEp)Q&deZ2*0PMrv*cpmAZa=?T09Ctp@mkrj2~RnSP;qB5jY)5R04nY(e#zHp z*%yqqZ3GM-LgxeUvf`kI4j(L ziEnh1f;s%qI`c zYK{V-5yuj42lE%;jOPY_1^e^6-k@s-bny_wx`K(N6}N3#Mk+yoI*$e{RIQdEz%Sq_ z8GzS9M-eQ3{FUnnKRD$N4iaEM6-HHx!Dmofk!*?blqiFx}gs|9w%Gv{R zN)T3{G>)E#kQPkCP)%eAIhBS9kRScV@TUpLM8{9cV%C7KCiT;QuB+qY1P-4_kduzv z`O?)fE5r$!k|JC29ZsgSt=|)NzOa*N7jgmDI zC1DYi#0c*jO2TODuGAf5BVjhLtBBi(Lb5CZr6U2Bo}F4*pe0GB)&W;opkKzs88sVj zxO5}(JtOjYdDG2#3H1^ud$Y3X*TNrju;87SSbB+!dxm}#h%RL+^R04o)NvVnXE%%B zFN?#YjmyM;VUTHSET%JLMygyWv@Iu7VcOeZM-qS?nNJtcR(aN z;!-VU@}l%1T9v<}01@KXzgVKrb#$4=Ir|~T!pC*xB>Ayq1Tuza97{V)Pw368{F-^; zR<6{4ut>hj$wHh+WV9E!V;(X+7vvQ)x|OZ`aPOPZy&s&CLOn(>)6l4+9{ke@?E@WF zfj>*ERDg_ey47m2wK^V(;c+Ze@`I{i5q@ zw@uMEoDJ+wfYUXXp!5aF(I?T51=xL%wt+DnhceFci(bx8P7ju~uB109GnB-%W8}VY z3lXN0>t(GIIbX%pWiRY-aFn`iwoP%Kk^?(pvd%=5v{cFr>vTMpByq&lDKzC|Ca8~p zNGCK_c0Xe~BMg>ngBaT4?vv}B*q8wtb7Fm!J=Cu%OIXDZiW>?m7qni$ey@P6=1&)< zC~WzL!i#zcA{3@gmI zlbql{{TX{Y64j0(Sj#I%sZjKEU?`=3?0yy(e8tQT*bpE^q-&zH^r3zpbDR2%e=&Y+ zqEEQVCB}ZKJ_gwRrve7cFyAg=T+U|&Vny~!bn6s2wE1!=>kZ&KD5EwPZiNQzN(Qw~ zhaQh^85;`*BgTbifrn(SP!L9$4luQYo!XV#41qEcp3N?=nw^D`!|Q$&mCp))QM)0Y zSje}5+RsVwxq{-%hGB`KEKZ%@8hD*enFop@j{)S$n_H=>$YF|rB7DPHYxc6PTdU2&fxzXkxFHKkv+i@3@zDvS6Y zsb9{^$Jm+D3owQ3O^f(H@RvZ})DDZNR4lGIdcLT#W@WY)vG1aP7wa0^>r=$gC9GFX zMgX2*bn`lZRTnf402UZ^J4?APz`}-m)~SlLMTwK@M|YTx01JSw@U;(^Ye%OiWBs1C zc62R0u>2}_Go{0A+WFh1tTEmWGG8sY@rJ5a>9KxG z@7&r;fee0()?VZ>y-tP1`-~-CrEu?8MVEXG1O$LK2M~&&42>{i)y(0RrjM~P;n7qP zPLV0ksB(oR->Kreq2m`A)*&>CMp_tz47Yat!YmxUC-nw@?HX{eI8MB#PR|0v21nV( z#!nuMbv~hApD6eTgz8Uzqp0tETOuw)E@-cS<8!7K91^2`#$DyHlqpg;8ao#C8Gm8C zqDf4~cm*|-s(LPT+c4%&!6s&~CkcG3eN5UC!G(ScN?JVy3J`!}8&b&HT^a*GL2`S) zUr+Tc30^;cskUvyY^oV9e)Y}bC(M7{*`NkxuO7pbsvL;6D?O>7z^}qN+N=3KN-x6< z4_ty&2hvJo=`JPcsea*T5(z(3N;D_Y6}dU;26?;&{)vNY;12^sX1P%w!vgP1 zzJc>pQu`cEy~cS;*q#J3rPKxA3S^w9(xE^5jre+hfwD1H?zG5_1JV7;+?FdSn^Tq? zc2MN8QMO9C#M3zbR9!)mKfvC#P9K@iiu4p`F4`OdzR-XtQ66(ElVNHL9p9LFWu3!5?D!vu$$3puhN%`9XyzRU|b{V zl6illoHv4;y8DrDAm=r8+U3;StZlifcR(Uds~fN>TIKV)pmV_$l_aTWz`ltwTRxGt zQSJLir&ex}%j!lVo%)SCMor?WNdb^?&>JCtgKk~C9d(96>LJi}L(8k~KqsZDN5lP4 z)wrcYl{09OR<&)X{UrqjWYMpv?aCNLJ?nIxut7DU+OP9)P^GXe&93%g(Yh+pB|ne- zpa3R}t6;myQi~P0YHHlM6&07t z7S3$Th1_zFWi7BJtZPpSXVfKH&FkCNd39RkW*WN0eqUPD{DodQ@A}2Lg_T&5@qlZx zRXUxLmFd_0k@MaUGbON9Xl%U}-8S2$0dB7tI#senqBBo+T_r%(#J#kvCD`r;-s^52 z^TNUd?Iv{Wd3JZo$6R3KXDitHjO&rKXT1?&EhQT+ov#FRHDNeKviY ziH<7JjLpax8Mr>>W$5w0xGmFyne0Gf4S-VKK%1F=u&12gxmk|_Ir13Idaz@EppsYn ztHLCM-Don&vI5D4CFdY{aDf>ww|lnE=V62d*$*Zrs9I_m#EE%@#bfpeY4^LHVA<-; z29K3uU+N9ck22`jsE9Z0Ym&|u+1ZM~CL?}r5n-yocU_8(q#V|^1D{_Y=He+E>k#p6 zQ^{F)G!xmyNwq~ppOQdgNo>Ch=f@^QgU=PT}i7H1T| z=kOS@288K)i8bV;yzdt|{r#xe_dP3pw@&n@N@4L(UgzPyrQ z)8o{Kwzkxx>6*Hca1%S9Mc{|NuDF>mXxJ62@7xi`3Z%oI0)?b&|JOP#s!p`(mPK}Z z^H)EU{LnmeZ0U?L%qB~Jg_eW#^GJNHhLM|%qQB+pcBvl>#n72+y~%Y0d)@XU^@$ZT zTH50VbToO>%}48Hk3H|s2YWCeA6m1nYNHmu=Zo4(UF72Y+G?m<2gmZvHzqke=aW|0 z-`L)^Fw=)c4xG=SEcL&fr^_?{#9e8z!w`pcE0Z~Og&~_+scbWU&z2QVG7s1K>?{$R zC&-`zXS}Z(H5<1F$}&*xg&(TuV=6%OH3!Ci3;i<5<>)HSdJD95>xFBT?>Bh3&Zm&i z-(qGS+g;I(ja}cH97#F=&aETI`1@A2QQguy)oxzwJgJ%&$I2nWH{_co8KCz{{37Oi zHC3Lm#K-f4URYUw)O*si^pv`|56adlHD&6mb~on};JTfcIFYhaeUTPZ_Y)49s8#x8 z_xIcbrfpNuX*$sKP)dYf__7oGcbvJlL;nxj*h$o7nZvoV()D#>4t(?eHMY zaNSus`LyL*G{nLN)9Jw*SYvWRpG_D<$>}gly-KY+dklDg6}pLYbOVj|Cehl;L4tSL zRb6F|77F5OS36a6zqapFNBA%qqZn(hr9qF4isBx%gVmTQ>d%&muQPOL0cBMK*>-TG zFTLyX!55uS8b4xi`DSwps>ax!v)xLERXt)I{^F4Oix zh&B7-IjP2f49`4`&Kh^gArKT)3-`Y{Jz)Qa9GCUBQE0IRZ$DVjE4ryVwZx}wSSi8Y zhj3=P%@d>k|NV1C>Xw4&+M0tOFk{_PfEK}p3Tu#?&4kckA9R>m-qGQ03d9=DckJ<( z4NZpHjYpGLeqJsJ7@6`0xJV~uzhT!eSXdyZ;AoP61<`9OZ{(#Ii9OCvOjPf6bB=;r z0Bcm*g&3j=zjlv~8<>r0|J`;>+y)PRm_vfKJ2VB=y6rOT{ny2p@?A%xRr@D5tS$0~ zam<&!78>QwpH5g@BfTiYyio+?i@6)gC3R8hAeuE<^r*2;hN)oNWS;ExJSiH6gH+~9 zT8fx|Hq$!tjn7JI*KIr#x62{Ul|{~5Sw&28-vc`=DP&x?uDf!0u$^k^uJX}#6z11J zkRZveJyYsNan6rj@|ZB-%$@Li9aKH=M$eRg1m9UssxHpxjt4qWc~N`T0doC1R_Emf z=^I!@to@Gj?5A(wsa?+VX%=}!WhK;wpK5}CFlEQ{cw+^z?F6P<7`WrUiI2M=n-rJa zq<@vSfPl7N6U4+*b517z!n$CVVLcnWUeifYWQGNMn^}?)&o+e~ z`Cuyfww22^D8tP8W_1ZwBFoCOC(n+z|3*K1;2mHbEOe|HVW;KXkk6b$eOnr<`V(j*V{perviW;fCWSd)}PCq}eQWCH4LIR=IPPXxHrU zOzvCRHE##pf&m)^kGEh)@_Z9rk+bE0F!sB|EFcQY`Sonb-IT2LYWdc>iM?}%k@wVN zwudXCT-M60Y-KJ#Z58Vu*7~EH%QWek*r=*ueY$q7 zMb6y4m_YIZRhS|GO(d39Q;bX3Qp4sW+I+p>ngoMay3&QqVkK^zlgF8>!Z!;eHYZnL zt%m3%EL`7?8S}=R9V^(>yRGtn7k557?i`&;X=?dP$6s#N%v{wYgK9y{%}^5$1Dn_Q z=2{Q9P6_6JWE(1<4$||FT>ibtlZMVW>wBHa5jS><6NS_;!{&Sm`{>*Zb}P-u`nDxB zwN$N}+%77zl)LBWO_iL@GWJ)w|IKL|-!aUv#-Qh`>Ct;;m2Ohp_A{)1(YJSGF<~~W z)jWw*4!4mGoU(~E?6z(;ah|&kow-}~$}?wn8k+RslQ)6V|M}#9=0n7GsiZ`gtKO4* z!&Aawp(^T)D=A~U(Lqt<*omRxp?)L|F+X;HX>+o-nuEqdiB`VjSs(8wW4kgv-Pli# zHF%zqJxp@UVW)%UHzWstpn!bR{JX+DRJ#E3S;X1}6OV1*zjk3QZOG|!jkJLm5xX2s zFi&Hu2zIcvS=VhJAZ>K)oFEo!Cih+6x;{=wq$R4FsoZ!DM>M{0zlK@Q-VDjjIXD(*_V9>P2n_Q>ci*f9L&$v^ZV)MB|EnZ%=ppnjFV+_TUR(=alxuCpQ13+AlPbkBB1{^*p`zt}3E@~_*{|U|wUc|$) zZEZlEM$GUX>#NR^Bxuu97*Af%xIKgP?QZSF<$rkc|1NI`mX|=r0UHE2H!wDrU&aAM ze_dH~9LJH~XDojN%%>TL(X9K>n{`;?hAmSPO+QF}AOQ)|E)R(&rPcb+Z)P1`S=Y<} zn8B_>L1Ch2y1FVW^UJfce>!IiKWF$~`g{56?9rE8%Jb`gp0V>!;otSyPiH**3_f&($YZ`cZ-r926}o0_dVHOn1wNC-R7$;UzF2KZ60qg*di3;TE`)TrbApPE*CA% zF9csq9Oj?f3$<{VmD@buDwv-aV)OG>vBgQZ`F49D7SFgcINT>s&%QYO>0GMCe}x5I z;0JHz!mK;DO2DHX66wRoXOABL<^0{-_urm9`ujOQd-NIp@53h_!Jo$;pa0w0Cy&p+ zn3ArZbEX|f27fyf0QNtVS6F1U%1l!>gdiLwmj?n!tP;tPAr#f_xmV$eJA8Qi=*5+-ub;^6s$6IFxTPO7mmuHy1m*p9lC_q)cb?fr}XR8;W zv_0z9NyhGF>m&m(^m_l+MJTNi9TuUz%3%D>W-Aw^tXHt+fae@WcpE+ne*o$HQD_N2 z03e5p5;*??=ODD*you27?B+YTC;pV61df;(# zF|Cb(T6(<|deKh$G^kdde-P#s4ygc8MIHz@gV43y^!<9%xA-+ESnhzo^TmEIaRbZk zKLo-vW;d5*axd45WyZK0u-AN%T%KQW4tvRXWSw77Wh=)uKtMGeyj3v+i-nm1u>?=g z!;OB27v*9#uYVfoB21d2vDMdvAIa4#7)r>+SWWKg2?LJ@$p&4of22t~49W>48*B;B zWO$I|YLS*)D|}mclx`I4t*|j?G|f~UBL<-ubd82`Jslf(4Mj-JSjk{+e?sd|QlPG#Na4v2dvJUw=6=g~;frwm~Df-V#fDFtNlc zs|k!YkdbILS@5(1OGFz@DVb)2K@-Xm#gZvT#AUI0DvLQ-9={zstT)4yJDwSq+=u@g zJIPp!O6=rze-b<=O3i-_)PIIsA@jmC@NZZc$1S^I1POE%yN8uYxX4{u+0YysL#^Z; z)FSJR1&e&Fp%!(r%W@v>>admrUlAO*qbmz68(i4qk#-2S;CwC|eAtz;jUVek#?*?m z!aHCpU<<%6@LnTOC=cX(9o8tcr&9KA;OkuCpcosYe>521PJsG)74dy~uRQa6I%aq# z+wJ`pWG#>%7Ds%om*r|6Rp^Ya2u2;a;#~-?4s;gikD^K5B}Hbb>K!`&+zIo+3B7Odm?K4>uV6AVuE;E2Cqu?jRKo`f*cxfbvkMdz zr0V4tkq+e-3_zpxKp?~v({Mp>_k*#13B-Vp-+Z$^{c63vU2fxjlsHo)o8p|1tt_!1 zf9EVB8*mBK7%4`&=&10DNcOdObG9$)D$=ei}f94`tuWo}I{^!T zLB#ii_3B6ssi1^pa4Gnzg`RFX2Fxw)1o*o!xT0z%VLk8!h%bP#Kr@!8g{;ILf39B+ zkXpm9y%$+$lVn6$`j~Vpqg8%og@S0NctqL_Yzy`wau(&m%NU^xR&nOLY^%J*$TQH^+loQ2 z)7Qk`ELO+skQ7I8F5^Af0&+_*WXSJ=K&VKkeC-udOc<*i9lMNES!nroVkLbHyJ@ua zSz?qLx{Jyz!g?w&$?K_X8OuQDWR53!Y@Mh6eQtL`FlfAV#DjxzVh z=;+RHjKTn}2eOZNy4-wKJfN9uVD*6Iow3t6Fhf^s81h^FGrHLKS(_31VCjN~ANtA> zTh@WDqXT5nAT^mc791TsbEWsWT)Uqhqkvv7OzPqxMb|wo4A^v;^Z8+!k-MvT@0aQ)DGilF8EF;$8g0?05=9ksQ?^-f6?$GLNIQzH9=p}n`SI;}DNh(55Fr37^@)OU5T#a+r zyefLx7wsyvjL}vgP4=g7%#1WJ{St01aQkkJVJT@dceLJ4_Hlr4i}$9fWSPLH0^CAp z;=9$Y8T%zDa~xRFf41d~%}jTa;RADySlb40T(NJw$dM|-e=b@oj7g&}g2BW$__bYt zrIIsWxGX^1zk@(+T?aVT|M#~tvv>C1+&!%9J+fCq_9`=(H`ya1MaY&dGD5NmA!J2n zN@mKI@&EMw{eJ&nxxdfzcs!2tIs3fMx%WQz{qc;ePIwdQ5S!}xwt`ES6l|YhbJ)Qm zB4&HSspn2d_B^S}(Sm4QPVAM0#zO;bo(O^;kFk6!B{Z0otErw4_`G0-B@s)fJ>Uawm=KGt50m8#lN! zUe*gjXHaO0QtFjL&CrF@B#lq$nQuv{+e!qA~tVS*+QO5I#c=2$7_<(+*YNO;wgCUt}D~1yY@ecH`q7XU#qn=RU)anPe{H{ zmNOFGQ)1oD3%!x7D-_Sy#3rX$J{1BSD~m~WC8TaW1*;A>HFPw4MFq|CPkW~$$hdW) zGdPRI#&~#|nPp2fYV?a_41A63eIKHQ(zVKj4wDAJleWg@?Pk#$0>9cX{m7=}R1Mn>M5w?(mpWsvV{z zUxk(y6#sb_Mqxk9ZqnI1G29|*WpEb9#j}9h*l&riqavo?Yq3K1*s{Eu=lMHnM&ca{ zpRZ54d^|83UFt1M9BV^}CF7R{Neg9a{If!9E?L)GlYg}7@|MzN_%`Z_5^}TXH@DId zzEdmS8$+j%_AKq2N+Q$0n62=Wb&O)@iy-Q3R8Z?#|Bm*t$A^u3kQ&-Do%ce53#djax=bF zRynzmyBAbXw{!&hL~qi`GdDqLyC|0@9KH@Ngx~#gHY>g2?o%qa1&wP@=4$x^-Y{zW#tlot4*O@6Y9BT} za&8uQGWiU}cpZl5-tzzY^6lE%3i%aqzb2^y2xLU0_+98c_ zyf~b}W%*@cfBv2s`O;?Q8qt04w(j!cxsB$>Ddw9)L_BNf?}k3}KQVjevf|&T)?<0I zJxXCTka>eq=e zA<8GQ7)oXIL;gD9ybU{#dc%f$>D{|{#+L&fel4ifTZ|xao-n;7j+scOzej4%V4w^o zd77=NIO3^OJZZ0(xk+B`kc*M<)T-z7KF@?*r5qQJCoI416(q7pSz0#A#(!oRQ8g67 zBSmua#!V8Y8sVSt!vygL57YSwa>x4%{uIU^R~p#*tEd>&kZ(C8yMs=+C}c_-y4Fl@ zuNB(+2~W*#n#n6{1+3_U=!>g8L|%WsH7Mbi-H|h7&IneqbrAdMo_GtBQ^6}!qE)IH z48CN;R3T82KXsYK-eYNWD`t%r+DKNVjE+=Yq9(R_=2W2914$w~WPW0z#V;V;_NV|N zwd)<%D0@?iLIwTNqeE9@y}48G>F5($_CPYIgN}u$)+x0=Yo-O<{L>&z+~BFZ^cx`}LdWoymy0+^TDN;G0zk(uMA zx1?7X%4j2#)3RX65*rnT&x);Qq0cQ$k)Ir{Q$EqfU%#s=pe|BAd3cqEuU;+h`Wd|t zzsPdCKKq7w;n;JXfN*HYTEiu$(18#aX2GF0PIi-bK8r8qNvDY{RcRx{v}DU7{ljd9-tXlb?*! znSMwuY76xuW_;MEtW5H$@nvbIrRgpBt?;b}8Pg7Z`Z>yoTX{vk4?Xj%Z+8n{!W1h6 zlhSeB%~Vw60%xauPUFx!wk^CtbgA=*|0A*4n*K1ye6?4=ef#>SnCcrcMV|ZoRue;n z+Jq95q+?BDE7!AbuMb&w)~Xw1NHfe%k0r!FrG{Qoyev?(NI&urJb7h3_OM(JI2`if z88j{@SP8u^-{dT8s$t9IMbXG>jWne2Fk$L8nS7@Yxhr}dGZd0kZ zd?JmPXotRA)$$6e&0giQ7twC7oq%Mz=kZ*BkjVvcdph_aOBn-K=Et}fd{qaIQWp{P z55_Ch2zDGPR%3tD)pX)87WMg*cqQYQv5JJZjfTyIPM5ppEX025J{2-J;5t!5j+|)A zcd7%As2k;1dAYYkcw>LD`){k&IPkNV$yFbePzh|jBN>Pnpm|?Sl%^ewE+L7}%Ek@h z_3!ZP@#m1yAl<@LrX)p^Z!D*@H`xvDtqSd|f6Fl}`~g)U2wCB^9J`5Lznr! zwE%pD!ou<#A(|&F8~4qOB2V& zkpk}F%5_^ZC5%=iWxKg9bIWf0ArFOcHO^Un#|WmI<9?#O+Cen4eL0!?o{l#2TSk19 zpM63*3N54HyWxZvG+6<0Q!{E&w+yu(qSMKV9??aoP2QB>3={9=wI_NPTR zogx-?o*Gt=0!+Mi5ke6Xc}`C~XHoWwPDMaeiQxMzxby_9(UPwt4BP9~*-l-i?G74vt~20B_y$ zmU!4kpL~{?#nw{=5gd-|d_UWdeUEmH)^{gY%u0lUZMl}amuq9YIL_mBlSj$aZ*$b2 zOquYp=zi7EaTTZBdTqh>Ge_f=-}j!n9?+z6zdqv;2Du=_7hhNVWvTMzX<1#k=y)KB zxpAbQvM6IM~4mYL0tJ8sLCwWh^%IHvw zIwTb_X-qDya&B7~aSzdO#C^Hd(i&bs=xzGLdGDExlYe<(V8O($%S`{+?)!Y|QF{2>AVL!YE=oUIEi-9)=xF`iU*; z1Wt3&1k$J^H>2t^zcgvLr|{2HUg{vaL>t3pS(Kh4ho!8JcfBM*!L{t{(GL)&@5n=@ zs)0P*+Fw2+yN|Dm9!m>kz33#!ulzd0wzMI9vW8KkxRvjMOZ5ilYLjWxrGaJ^k04NS zz|vz|-OH6mbwAlu*FWerj|BN!?TL7^ZB@-kVj06%f6{I<9G0yr=e^bN)cif$AK=53 zTyv+5bAM4$B?Dah!bLnR%H8GiT1%v4r>pk@=`@yXZo@ zSdOF@7z3uV@PpY$thjafW-7^lb{~6=ab_?6y5&r(uqiESIA8uMoiNUN@iSX-X&k(z z0W&y%G??H!5GKXOpoe4>o>hAd9WF6fjNna^t?E=CSJz6Ic|c*2vJvz8q_PT$n0m<> zQ2+~0l_=3S0G`#g#NT4)lqc}%_UO=*9xp2C@RNmFs(JhCY35PBFNND~qa?G91Q3ig=rnsszm18@xjdL8H>km%JLB0L#9~?@W)^B`5H^(KvcEV zUy7GBOg@UJ0Ht&p_!*p6O0q?~b#umq5#a~1hNp$(8@H*mCS;i777`6SxQrXoKQ`Y+ zg@Y$<`W6-wA{?b8YrYF@9*I@DO5eiTF%TQH{!Jju>(PfE$|KQUE!+0VCuK9v<^GtZ zqQ!NAFMPuD9~zN+>R);_JH_pXQ~GfDnb@fWl{hQvNov@r&go=PP1gyc4(xvDvQe`Q zh+GZZmHDA`-JWRGrk(A1RGA-Wm{8?4g0q?=9;dcrH4f*-WeSXy$(wKZLZ5N^a;{xG~*!uTwNyui6$S0bfmG^H=`z zRi|=G3?b>2ln6S0t}7?|0WF&>Ph-`}iT`W+O5Y+*$(JGEDVl?^g*jEXrXGzBV2*OK zGmVYxj8JRy@Su;&RYKG>1Xgj*N|J@4-(v-p*tw11b|J2IMwB$i*{|DKkm+GohITje z?NrK;YYDbL9*Z0kMz(&cd|}9r2q>hY>m7`AyubF`x^X?{xs?e$`}Vbn{r0LU9`4k( z$cC4kslshW;d6{P#2LY+`EqUMaTs~Bpw@|U-_Xcz6S1j;5enu+dBK2NZYRfLGy+29 z(xSyV-#J5#)`nl?GeGGTX=;}5nilZ^gU!WsRht|p*0D|=>Da9w=uuJjq=$g77e8|I z+8Se70i*p0y!(-*NAQ$^TacT(?o~~%>)t9dZqGlDQ(Sfx&dEw?(=MVW8pKSbQ_WI9 z8$x3Isq^h&L-tp}TZa8TTj+mJ$i#zzQWd`dpMzScC$hexC-llZRE2r1QG9|Od`V-f z4_l_1g33|oQFw`55qiojWPX$7Xy=H1(!|KEBtq~WIr07BZq=5E#M}1cC$0h!)9S$U zPjwGWtb&TPre}deGGE``KHp!nG7*;eIG+9Cbn5_V>d+w`vt0M;bw-iNC9BP38uh|j zvN>YR7eC~I31)Vlqn5G9=@lmMT8E0j=Goq&Q+O?6?X;|OB&N79{@1$f`3a|T)wFX5 zgL$Cs`JElrD(%d&abweH5kI7Eo0`3by3i}i)1={+rpDtZ2pUT&@y(!QW1P*ie&K7@ z4_{M|Ak=f&ILTe*vMf*S4{*Mm+;hL$%DvZ!IXQoQyRz0~!Szl%N-XaZnfra(c3R$; z;(J+KcYg{gn|qs7OUtK3M|8t^PS*DGsYBBHRG2Mae&*6-F9cD~7-~4_C)690mJ!pu z>k$w&xk=tPDf*LLmV#_LtfQ45PW_OgH9@U1ZGbrQ1CrX^D}*3VH?`g-=%9lsY_jN+r zJ^03?HkWLxt?%fmEfE;e=+k!7Co-R1A0HQ#B~=a=jSZjXyp!S)U^CD#*S00F74jx@ zR(IfM0RHRjw4FWWt$;V7mZGOApH*mC-oD1P^vV>CTJ`>b1dG|{7^R|?_6YA>^&DC! z-w3hQm=x(f?K-KH@Tvg*zS-x<`pSF?)xXC#2ctiOM1D6JPnHqWYYgejRg{uk&dH&p z-&x};=v7fBzibZQF;`qGO;mTS>wiw*sIbRf==?n3qwwZ4Zsg@Fq{E-vYOWg>Z(qe! zm@CRUJ!XitHYc-Ubth>IdaMC5=uNX3+HDJem}5_)gXslMwJ3pAX^?(*hLAFPYi&?5 zIe&7t7_JiIOL^{DFprSMRSck}vsw%%%T!=-WM|Ho$0LyvTrj!PlOGRNd7M}AP2ZUG zU2x{Z^lNsmHZO(gPUpO%8E{SK+KMk@9*}pFY{IoP55R{3vE;EYh|vx{q`{f~cJZWo zh9AEZURtIdW|!_y1EtfsD?fUvXSx-(?;y&JTn+m)WT&}>qg)CzoHb~=Mn!kyrdvFB zgsdzsPh=3W?WLQsHgv~nrgq_^_cEQ4+eSHtC{zlV3&uPh%Wun9HKjQBB+7;pLr0`+ z0el$Bh6 zEiS8kt~+w>NOP)}q3U(0{A`}nBbc-Oy&{5CaaK^|>rPrj-nLHG3(&_MsZhiKD?}lp z_tT}zc*HI~?9Y^$9piD^L~Hb0rw;vQtamM!Opcx!>{P|MRQN3 zWuYkA#ZW%Eq#}2?BwySwgjqpN-hkZO1f^_F(>z=~Tb3GlhS(U!Nl*Xdvnb0Iq@NYX zH*=~*?yc+(WOa1fKYz{O+r>lHPG`bGkDFt6X()%vYQv@zBMI)fyf{|Cbz?bb&2NlO z`@DLIGleO>CUaLt-d2C^yOqNc4r_R6Y;cAPl^dGg)*cw2!-Be6svD1>r}ybKUu( z+^mfTu3a05nB*CFM8^`BjH@KXBp`s)?@`fS6( zv!(2sQ-p$GI6AsOk8gPs<9+YF04@Df9vU%vmWN4m$hRq%5z+iJg=$Oe+6w8q&bQw- ze|u36eAAsmH8T#>BmYMA(8`eF{pSeOSfURONB1hLdjd|@?!4OU-ecY9tMx>FN1O3R zn7X;HL?630;LZE}fO9hp%T=`rF;~+Q*Cm}F?J)U5J*}bpb_PG>Y3K?CPNHvo)O?5; z@t=^usf8L9EJ)L?-ctD~&>Iub?wTOeB5{{9fo^f*>^ zfN#%J$u7ypH-D4j_2NvQ5Kp?|6`FR{s#lz*MiI_vRvLq8@#&?gA9I{j&E0ZaJC3g| z<;V1H*G@e&%n}Qp5b~90%Jw3lnXqhWzPCQ)&qyyY#u&ff=AbWasQEJB-03}ErRaCb zQMp#OZ|`F@)7lO%r)2lC!Oec;D`q_rPtr-$E>YdZ1?SKWGzBNWx>;dUzwdfjPukmZ zedk52fHOxdW_;C-=6pK*Xdw6_M^Mn>4Xt+{g~}ZRLk31Nct-lw)P20!#}~h8GIfRR z`44YY9qu_Ve=uiLRk9BK8Tmt&edjS!#yA*7@Rfn0+_q&>t)z7!#W2+2K?ANnp$o)D`gwbZe8H)QXmaw8%vHuC?gTfme z{@R7gba4DZ^CWr+xsaJXuq(Nf9MufslKwQ&`#@nOzKnqf0>9#zd|yhW36bLcNkVr# zt9}GELm)N#UNv5(^O>;*ZJ8^59oq=KV|T@%5#L?42OJvi8}x5tXP89MgP?szX+Cyx zHWB^PVa)o?ED;iEv5whGAK#nqhE`2dYW9sKaMM}3G^+JR-k@c9|MCOO-*<7=Yc3~A zF*p$Mz_|MZZ9t&=K&+F+@yoK9!LUcbV8_dQp;=FVBJo}L= z0T{e@KV2)cz6TDm`uj``q9&v23hSfoS*JIu#K4$d`-a@q$7c}UW%KN6V}_eL)}hCsEJr|yHxc{ZLeF+qG#H>CRz=0%kL=6X)P6}Pr9 zt_=SgX8GIAq=7x<#|KE-3dB-t8YW%I94Y{%=jJ`tZs}lI0V*z7lelg-d;o&NYkgYcW(;De&w56IJ)!1)tH%J z&XtDmI9a3|RJjs6*qA-5CKjum=5lMt7iCjZnkMYm)U^JfPYnd0MyB&fg@B)%>1l^t zNxn7`e2vtZIoh7;E*;{9*8w`&XgqCP+VM^eX6gA!d1?y0@z>m<#JgerexvRK{e-9I z&zCr@#Yp1o23j(*@crhZxyivs#QEx31DYG9(GUk)$6SK)iPeY`}Dwepw>I= z5hle)=8)-Bw$*0nux-WUG$~&p|3LMXWm7=t63uttG@m*~1!z}2EUqVRVLS~fhL)=M z@!ZpeOXYQVI(^`n{4}p%2%c1lTUJP zuN0X?e)*&|c8U-dWBn;4ANJ7sWAbCNrieZr7T$?=tqw-io{oFuJhM-+xRD2P&_Y!|(?_9cE>Ssv{?l#dryCg< zkJ1L)A-uKnV^fFDJN|dFD2{qZYgZ*+W_mN$x*BW7;yB)h_r=nGcV<-2E*~yY4!NPq zNPFbIytE%LHeha2CNgQxyQir(di0ZKbeGpGuvK}SHC&^w?2b;W%lOZNHvW9I=dW*6y@r`dCo*o>R5+ zk#kesNq|54TBxIUUhcHr@%{o#yTQSS-)W6SJg2Ew?f1SSc%Iplciy^g@8jt1`k{)lIEY<7^9$D%q zE#5~zcrvL+TzK|0P$Vf72}seG&Cu2K4^E{=)}`py3E;Dk};fg@K^p;J^P1=t6>l*Smk+frcR< zKn(t*0eq0(9{ouJqoHUJ`fnOw6%>s{10lh1fDZ%Q{U5U+AUGO+&^tUl!AXfh)*?@Z&_@EFp|8(~NGy49pf~|*(LlWaQv~4Oo&u+u{l5OB zK`-P2s9``!s$A}(kHq@Y09 zHmClO#-}9!!I98ZC7ItLLu1t)6pF@b0w_@17itl3772rW*m~jE*2CZz0S^oX zL*UR<-mCb07@(9cssNA=tCE1Sz*dI90vjk8cCkf4zyfkH^dK}0`7gZS`m^W(^ukvB zPa-ieK*~is!9dsjP5UhuD<~i!KNsNuC}|j0O~cUWi&YDS1MLpz_!G5(=phjojRBNg zYy&U|iNG=j3I|^tH9)%nx&GS`V5|U>?Z32tqc$+uE?5l#fdCW#%lNx+fU$yvU+4&c z4~)JDsDLq8s14At`M?-7*5CuY!%9CG4y5pJW5CdVV;*4Gf8h-dMPT#6uwf7a2LF>1 z!2kch8L;`{`~gECXza%PpT+>|43>d#>;{JX-J-wa|62)G6d(vV^kVLyz^-vo2^#V5 zg8#!2fDdazK;cNNxd*hU6jloZim`GBY&kIO3Qld8#HYmob$6j8fDHnMjTUg!{}6#f zUc47TCTt)^U^PD&ih#j>1I*tIghXLCIueaFt6(6+@AE!?U;m$dfQB_QV8E>x0}jys zY{-8a2no0X5J1Cf3j`Qdl?# z(6FL~L}KGHK)Z;C09+vcHP8QSAb^N9_kKtJq9*_iOEC(909_me01a}nn^AD+Uzqk! z0|7+rH5P=th-m--p@5A4op&@C*wC@QUrY&r z#8@4HhN7_Q2-q*MmNGQ_KT8khzZDEDiy$l|2q+eCqLCnMqmc+KoI?ZgzgTc+6bklN zM-czUe>AYkq$;W6%UnHl=&`G0rIVtsDTSa|D`kZQu_%_Lva+P!J~m!HcfIT#NP!nM Oz#J##<5SjEA^kreg=H52 diff --git a/doc/musrSim.tex b/doc/musrSim.tex index 9a1bf4a..d8d6375 100644 --- a/doc/musrSim.tex +++ b/doc/musrSim.tex @@ -685,6 +685,7 @@ Three special volumes ``Target, M0, M1 and M2''. %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% \section{Optical photons} +\label{sec:opticalPhotons} Normaly the simulation of ``detected signal'' stops at the level of the deposited energy in a sensitive volume (e.g.\ in a scintilator tile). However, in some special cases, one would like to know how the light is propagated through the scintillators. In such case simulation @@ -945,7 +946,7 @@ The list of variables that can be stored in the Root tree: Note that these variables are by default not calculated (and not stored) and therefore the user has to switch the calculation on by ``/musr/command rootOutput fieldIntegralBx on'' in the macro file. - +\item{\bf nOptPhot} (Int\_t) -- total number of optical photons generated in the current event. \item{\bf det\_n} (Int\_t) -- number of ``detector hits'' in this event. Note that more then 1 detector might be hit, and even the same detector might be hit more than once. The hit might be induced by just one particle, by more then one particle originating from the same particle initially hitting the detector, @@ -986,6 +987,23 @@ The list of variables that can be stored in the Root tree: given volume. This way one can better investigate which (hopefully) single particle caused the hit. Even though even in this case it is not guaranteed that only a single particle gave origin to the hit, it is quite likely, though, that it was in fact just a single particle. If the +\item{\bf odet\_n} (Int\_t) -- number of ``optical photon signals'' in this event. Note that there can be several optical photon signals + per detector. In principle, the optical photons contributing to the same signal may originate from more than + one charged particle. Ideally, there should be the same number of optical photon signals as there are + ``detector hits'', i.e.\ odet\_n = det\_n. However, detector hits are treated independently of the + optical photon signals, and they do not have to match perfectly. +\item{\bf odet\_ID[odet\_n]} (array of Int\_t) -- ID number of the detector where the given optical photon signal occurred. + The detection of optical photons happens at the boundary of surfaces, and it is assigned to the volume ID \emph{to} which + the photon enters (and not \emph{from} which it comes). This e.g.\ allows one to attach more APDs to one scintillator, + in which case odet\_ID[] will correspond to the volume IDs of individual APDs. +\item{\bf odet\_nPhot[odet\_n]} (array of Int\_t) -- number of photons detected in the given optical photon signal. +\item{\bf odet\_timeFirst[odet\_n], odet\_timeLast[odet\_n]} (array of Int\_t) + -- times, when the the first and last photons contributing to the given signal were detected (odet\_timeFirst and odet\_timeLast). +\item{\bf odet\_timeA[odet\_n], odet\_timeB[odet\_n], odet\_timeC[odet\_n], odet\_timeD[odet\_n], odet\_timeE[odet\_n]} (array of Int\_t) + -- time, when the $n^{th}$ photon was detected, where $n$ for \emph{odet\_timeA[i]} is given as \emph{OPSA\_fracA} + multiplied by \emph{odet\_nPhot[i]}. Similary for \emph{odet\_timeB[i]}, ..., \emph{odet\_timeE[i]}. + The variables \emph{OPSA\_fracA}, ..., \emph{OPSA\_fracE}, are defined by ``/musr/command OPSA photonFractions'' command + (see chapter~\ref{sec:opticalPhotons}). \item{\bf save\_n} (Int\_t) -- number of special kind of ``save'' volume that were hit in this event. The ``save volume'' is any volume whose name starts with letters ``save''. Their purpose in the simulation is usually to check positions and momenta of particles at some position of the detector, even if the particle does not deposit any energy in diff --git a/src/musrPhysicsList.cc b/src/musrPhysicsList.cc index 7dee52a..de719e5 100644 --- a/src/musrPhysicsList.cc +++ b/src/musrPhysicsList.cc @@ -214,7 +214,7 @@ void musrPhysicsList::ConstructProcess() #include "G4PhotoElectricEffect.hh" #include "G4RayleighScattering.hh" -#include "G4MultipleScattering.hh" +// #include "G4MultipleScattering.hh" // obsolete in Geant4.9.4 #include "G4eMultipleScattering.hh" #include "G4eIonisation.hh" @@ -346,8 +346,9 @@ void musrPhysicsList::ConstructEM() G4int nr1, nr2, nr3; char charRegion1[100]="", charRegion2[100]="", charRegion3[100]="", charControlString[10]=""; sscanf(&line[0],"%*s %*s %*s %*s %*s %d %d %d %s %s %s %s",&nr1,&nr2,&nr3,charRegion1,charRegion2,charRegion3,charControlString); - if (stringProcessName=="G4MultipleScattering") pManager->AddProcess(new G4MultipleScattering,nr1,nr2,nr3); - else if (stringProcessName=="G4eMultipleScattering") pManager->AddProcess(new G4eMultipleScattering,nr1,nr2,nr3); + // if (stringProcessName=="G4MultipleScattering") pManager->AddProcess(new G4MultipleScattering,nr1,nr2,nr3); + // obsolete in Geant4.9.4 + if (stringProcessName=="G4eMultipleScattering") pManager->AddProcess(new G4eMultipleScattering,nr1,nr2,nr3); else if (stringProcessName=="G4eIonisation") pManager->AddProcess(new G4eIonisation,nr1,nr2,nr3); else if (stringProcessName=="G4eBremsstrahlung") pManager->AddProcess(new G4eBremsstrahlung,nr1,nr2,nr3); else if (stringProcessName=="G4eplusAnnihilation") pManager->AddProcess(new G4eplusAnnihilation,nr1,nr2,nr3); @@ -371,35 +372,37 @@ void musrPhysicsList::ConstructEM() // cks: musrMuScatter could be uncommented here, but testing is needed, because Toni has some strange comments // in his original "musrPhysicsList.cc about implementing musrMuScatter. // else if (stringProcessName=="musrMuScatter") pManager->AddProcess(new musrMuScatter,nr1,nr2,nr3); - else if (stringProcessName=="MultipleAndCoulombScattering") { - G4MultipleScattering* multScat = new G4MultipleScattering(); - // G4CoulombScattering* coulScat = new G4CoulombScattering(); - G4CoulombScatteringModel* coulScatModel = new G4CoulombScatteringModel(); - if (strcmp(charRegion1,"")!=0) { - G4Region* regionForCoulomb = FindG4Region(charRegion1,line); - G4cout<<" Adding Coulomb scattering model to multiple scattering model for region "<AddEmModel(0,coulScatModel,regionForCoulomb); - // multScat->AddEmModel(0,multScat,regionForCoulomb); - } - if (strcmp(charRegion2,"")!=0) { - G4Region* regionForCoulomb = FindG4Region(charRegion2,line); - G4cout<<" Adding Coulomb scattering model to multiple scattering model for region "<AddEmModel(0,coulScatModel,regionForCoulomb); - } - if (strcmp(charRegion3,"")!=0) { - G4Region* regionForCoulomb = FindG4Region(charRegion3,line); - G4cout<<" Adding Coulomb scattering model to multiple scattering model for region "<AddEmModel(0,coulScatModel,regionForCoulomb); - } - if (strcmp(charControlString,"")!=0) { - G4cout<<"More than 3 regions requested for Coulomb Scattering, but presently only up to 3 such regions are supported."<AddProcess(multScat,nr1,nr2,nr3); - } + // + // The G4MultipleScattering is obsolete in Geant4.9.4 + // else if (stringProcessName=="MultipleAndCoulombScattering") { + // G4MultipleScattering* multScat = new G4MultipleScattering(); + // // G4CoulombScattering* coulScat = new G4CoulombScattering(); + // G4CoulombScatteringModel* coulScatModel = new G4CoulombScatteringModel(); + // if (strcmp(charRegion1,"")!=0) { + // G4Region* regionForCoulomb = FindG4Region(charRegion1,line); + // G4cout<<" Adding Coulomb scattering model to multiple scattering model for region "<AddEmModel(0,coulScatModel,regionForCoulomb); + // // multScat->AddEmModel(0,multScat,regionForCoulomb); + // } + // if (strcmp(charRegion2,"")!=0) { + // G4Region* regionForCoulomb = FindG4Region(charRegion2,line); + // G4cout<<" Adding Coulomb scattering model to multiple scattering model for region "<AddEmModel(0,coulScatModel,regionForCoulomb); + // } + // if (strcmp(charRegion3,"")!=0) { + // G4Region* regionForCoulomb = FindG4Region(charRegion3,line); + // G4cout<<" Adding Coulomb scattering model to multiple scattering model for region "<AddEmModel(0,coulScatModel,regionForCoulomb); + // } + // if (strcmp(charControlString,"")!=0) { + // G4cout<<"More than 3 regions requested for Coulomb Scattering, but presently only up to 3 such regions are supported."<AddProcess(multScat,nr1,nr2,nr3); + // } else { sprintf(eMessage,"musrPhysicsList: Process \"%s\" is not implemented in musrPhysicsList.cc for addProcess. It can be easily added.", charProcessName); diff --git a/src/musrScintSD.cc b/src/musrScintSD.cc index ec14fd3..16d0795 100644 --- a/src/musrScintSD.cc +++ b/src/musrScintSD.cc @@ -231,7 +231,8 @@ void musrScintSD::ProcessOpticalPhoton(G4Step* aStep) { if (boundaryStatus!=Detection) { char message[200]; sprintf(message,"musrScintSD.cc::ProcessOpticalPhoton(): Optical photon boundary status is not Detection but %i",boundaryStatus); - musrErrorMessage::GetInstance()->musrError(FATAL,message,false); + // musrErrorMessage::GetInstance()->musrError(FATAL,message,false); + musrErrorMessage::GetInstance()->musrError(WARNING,message,true); } // switch(boundaryStatus){ diff --git a/src/musrSteppingAction.cc b/src/musrSteppingAction.cc index 4ee6742..3b53bd4 100644 --- a/src/musrSteppingAction.cc +++ b/src/musrSteppingAction.cc @@ -353,7 +353,7 @@ void musrSteppingAction::UserSteppingAction(const G4Step* aStep) { // G4cout<<"================================================== BzIntegral = "<GetSecondary(); + const G4TrackVector* secondary = fpSteppingManager->GetSecondary(); G4int n_secondaries= (*secondary).size(); for (G4int i=0; iGetDefinition()->GetParticleName()) == "e+" ) {