From 00869a9955abe2956df93fdf6990eb591c58974a Mon Sep 17 00:00:00 2001 From: matta <matta@npt> Date: Fri, 27 Nov 2009 07:37:10 +0000 Subject: [PATCH] * Fixing energy loss Table writting --- NPAnalysis/10He_Riken/Analysis | Bin 55635 -> 56993 bytes .../10He_Riken/include/ObjectManager.hh | 7 + NPLib/Tools/NPEnergyLoss.cxx | 167 ++++++++---------- NPLib/Tools/NPEnergyLoss.h | 9 +- NPSimulation/src/EventGeneratorBeam.cc | 5 + NPSimulation/src/EventGeneratorTransfert.cc | 21 +++ .../src/EventGeneratorTransfertToResonance.cc | 29 +-- NPSimulation/src/Target.cc | 3 +- 8 files changed, 132 insertions(+), 109 deletions(-) diff --git a/NPAnalysis/10He_Riken/Analysis b/NPAnalysis/10He_Riken/Analysis index 592ab93c584150d90ec51a9f190d63e0e5f2bfa5..17a18c713490b00b1f4bb3098d55b3ae7fc7d93b 100755 GIT binary patch literal 56993 zcmeIbdtg+>^#^_fSqT_Sq^OZnT@@9r#%x|BikOfr5C|bnf&z-mX0s#<o86e*K%m-& zhA1IMsZX>%;u{}Tw5a%MgD9e3TBA~p3K#`zBB)qVsbcf{eCE#V&TbxU{o(Jg-|A)V zXU>^3=bSn7+`BiwwHHo_iHR9%$j=aCs1bwv^&4a3(vePEV-nI0s}X06H%>Ih8b=_D zwDbn>;1UH;$b!6BTojfLoW25dTt5TmMH*%E5+9bu^~7q!;KlMrqDkKXykU$bj@L8@ z9E*I`OZpHb+d$t28m~BH@FJac@S287rlD?L`M|sst@z_$#hPV?vF2Z~alE`J%S+LJ z1v=RKe3bliyq@xNycK7AJ-)j7<bXdp)hr9ugl$ucW{a{ADRdfzl(I+TqHjnaV=4Y< z%Iw_7kALLD(m&X*Z~D`nXME!X?E+j@TqAHTSU>8dnAJuDZm%kE9!(i-IaZ)VBX(sD znghaCY+CZ*2+PTF`8#4C2EAeIi1uF$vBn)U`j}#CdkJ_)#K#y7a7JQ|b;!{vNilIo zY?9FsmuHM_IBCe}hK5xQ#Y2t!_?VlGZ9g+sH&j_;hYu+q{xc&saRGJ7!NpH5t~^{* z2%<jIB%Xn^P~Oi(I!oRcBc&>(xJ=KOgo}{O#WfGtCAj9}x)hfKSAU;!WV&!w;;O>6 z2-jj<K3ui9g1FRY32v6+s+SqRmg(h4ufX*iiJ3IuT7jz(7e7t7uENz!fX~%tBL4k7 z{=H7(o3#9!wfkG-{U4CtfvZj8yJdQhOz%Z{AFlf)X7ZrCe;DZ_xGuQUcXe*l_t!7? z$L2ZX4^-_r$(lCN`S%+~o_q6=Pu=(R%H3~t<a+=3_Z#p2EaS?{Iu7l>YtiuM{ioO7 z@$6-Z3v-7RTxdCcPQuEm`~KqWy+82qr{6msS-bG5(l73L=7j&d=i_76#^+70ereoE z#y5Yw-Iu@mwCaNo|1kHJah^Xddh&<wPjyvq{MuGuu=uIT_P_6IPF?-0;dfkh-raw^ zvLh*FXWG2?U--oKZq?G0emF2Ze)LCoy#LqTcQ}9Da$fo~Tm84aeZDpAje_q}%isU( z)3?4j>hs$cKJ`rNmq)#KUB=n>{l4V)pPunRVA!+CXLyoo|NX}qyW2Z2n=`ik$s3l< zzU!*DRzT0j%Z*o_+kM^TMK5}GK7Dk<JtKp~oAx|@Wk*cu@;8slKeT`NeU<s?C*5+) z?QyS<`CH}f?hoI*r#tzsg-_k|deR3I9^0LI?{Pcxx9vXqo$>a+7oM5?K)`>+74vVr z@)!Bz|CV&k^ar1OdRzBJpS*L_&=p^eIm^9hc;4Hd;z_HM15f;IZue*VE}u8@&Ab0F zbk?22J%<+k_SVZ^A9nuV-W+z5@72KVXJ2v4-yYt*@3cU5>H~o{zWQqAr0361s;^F5 z=-9Kv*KqxfFP(AL-X{O>zdzY~<%Tb1Glr5e@u6E|6mBcn!R3cIVYotupl^zzZ#*)R z{$v;m{=aZ=k{`;unuCx05U(7jkSRYOLxTE{et|_Hlb#e+e_mAihob8HJ_?@`7H3#- zUm8WPjIw8Jl>K7BlrpH_wQ&+4A<l{_KP<}LB~khhkCMMPO5W>H^0ej>`ad5<?+&XY z^f@iczM?37E{H0>F-o72T7|Uls3`r1MbU4Is-NSR`mlU_ls;cYwfD}b`gcdQe{B@L zD@y+tqVS_p_AVK!v@zT3nJ9VhMcG##RsPpe@|Q*Fzaq-M!YF%Aj*@?JRQYXT^+fyq zGm8JwDEz0W`0t3S?;TA3_4{JusP@UmNMrrchOt+}qp$#q|136cousc$GK{-4`jH8S zal(h;@++`dXgkpoHxc|t8)J-g$UpcCi-=mQalCQNQJ5c4{ufeS<uQiQeu5>AW|BYc zSi|Vj$`>Ao`SEjeEk^w&p#D`~#Ky5%h!4Tul}}<lh7aQ>Z7k;gJ+X0iN&hqAZ#wil zQR17x6#G$F-Y?;w8=eZ6H=xgzC|@G^mqDM6&<|@Aq2Jvn7{;OQu)RrW|7W|x*b2Tf zDeq+5pD_|^EBH_OcMJUwc~>K$z1`672-L?l7VUW)6=i-L8}~Beka*!3!$^j`lO+9W z=(llqY#hW2dH{~>M0qwl`L_eF=?T+MgT5P3|7num`E%SrzQR{Su~^Oiv!Kr**<Ojp zm%TT`pGol7I5ZIDUxM^`_>*FYPeuD1XfNfTrw}h&pdZv0^s~@@R;@jM3H^ewhwEyV zzXJYv>FKclEr^eH*}kJdKL!5p0OneYxCi!cgMDmpmS2GQ`WfOw`8N&vZP4sldc0vw z2Hht4v%t?i0-%tWfch@kukcaEnP|@~T6><2`bTe2^dpTr*b|5PRegU3|3mNxV~_f@ zp<&w4zN~=ldp*)7$W!sU6!F^yx>eE-LeS-?-}GOCkrw6u7f`<!_A7g1V9z&N|GN(j z@w^tFufjiV)^Pv30P;G{H2WXhlX3JB;zjkpA?O&_XzlSh^zTCZoeDnDzU(hr`*xt= zPDA{#*-1YQ4gH90FXDILkBKOMhNQm*f1LpTpCa)k(Dyn`{!NhgtG8q0ri0G*Ek;LJ zt@YQ5koTt69(BOGfYYSB)$q?EjsF6)Zx`&}gSevq--iC*ApVuTv*7P^Ek1q?f34Hn z>nZ5BK=a=@XqYav_f}coc^Drb1U^~f1oYER)PI)5#n5LJ>{b117uxGR)UW!>{fM_a zwD?KY#>YgAkG-&OyEdL(0soAGd=+p1f`96@c*%x++ckO50k4LCl)Yzz->TVjI^v<? znb^3`(68z5>!JU;Ph;a=m;L*Xu;U7pSMhcc^j(PdQuf^ed3U3|)p&Cm_#Xg2)D-b> z0P(Y4^XC{CzET_Svk^G20xSP*M*N(qwclsJV-TMz{vHMYXxMkTl;4Q)Y=JhO<innL z=%eJl54#?P|L0471AE>)f2|waUBfYbiKE0_=?MhgwRv7=Am9!dHNmP{x3j|Ga{63u zui<c1R{MR9K+stmbT|yhywadISl}uR+H9^WXRRYx>+}Qz1ulE36aU!DOJ<ao7nC{5 zQXFMAhrKk5MT#br%_?8yb_IP)?TbuY8m#sBD)Vgi;u$XFRCtS*lt!>tc<seYK=CbI zB+?>VTG@r}+JMLJ%d2v`7Ne+#qO!_LgZ^4)rQ2pR>&W$b{jSJr@@@7d!fcN(NVDOA zl7d<CuE?$Oroj)wE0kK|^i}w)(@>4jIFun|q2BbeX{l38{Yz~QPeld(r#>@EgPBWU zs=v0N0xqer+tUM{<!&fo7P3VYqAF8soeQNdyiKc6pvI|g%WB=Cy{LHDhGKsp7;x2k zYLo#po$xBpF(|gUYFWVJ3WN;0P$)%%&_u%3S#T&7O;a|_2<6!3m<nq-bNl4Xi_S@# z>8v#KLV{3Tg>+0PFD>66R7zBPO2T!(k<=V!RD~jx!DF_CYKlG%_SBng)YtX4bjM6b zt<zWOcDNWcDlANa%IrP@Q=LI@LX$Gl6{JHeytu&$k<@Ue>_dPT=DUN!laf=xnNA-9 zvsN}iT2XPXEhCBsqeNts2D2O`<>*arcR^{OG?-rF^wbv2C=C=Bv0Qej3cDR$FU3*n z_M%2xAXreEfo>K|_qqcClw>z#j@_0C(>=a`8{M|f=eexTO_yhghVNGsTh`<{kGH}M z8>iR4#O(|IRJx^vb;AGvgF?nTeEtfzqud#Am)ad6g%7`XCWJLGy+bFarl~e?xct6A z(BZ+j=ww%>OGDc9>%LT52}aF4zi*+ZvJQnXPEGN6-S*M|JAS`KY!mXlZfC9VwCt^L zIN~YKiGNHPks)Nu>Ptcdp$-%%m{sEUs{zcGLNKc?SW}0Qy3!pifIS@1?GrFya_lWl z5wV~%PRMr$g0=o-X0vBgLV>T30JAn&>!Oh|xM#q^!Yq}UnpWZ#9;~f%1xwssxMK+h zSC0#WffytsxNTW>pS!klS)sGM%wJaJ4m$H@O2!C2)mI`i2k|?6A_e;cY8(`L+ET*# zspx=CZ=EyfF7SDR9t6J7MU7glOJ`SOMs&Vth0=h>;}PRST_Bj6?+iNSfH0FdazMs( zrst-nPjv?|phE?}Z`Q&BA0tmV*qr8~*lk&0FT)rRaQSQ8MgH1qC@!ilEiKL0Y9`bw zhD(Yu&!Z(Y9+?x#mh1D?d7ZVTHjQzrGf?BKt>~K*Su}t4%$al5IHxg72cwTv2s2fi z*s{v{kVIQEsKbuP^ZRQnJU+CPC<tO+;WT^k*{R7X$?3^y($#@rh3l8U1W3+Aq@hoz zI2>5f_?J5TbwP)Jp*bCxs}8I!7>%U?)P~;8o*FS0Rk&TX?rMyK#Sz(|A!=a$0md{@ zAl>b&@b<TYvX`KmG>m;7mjhEfR$<lVO3GfEiN5M^)jJ)vZeffm$OV&dS0t8RGBP7H z_4t$l11nTi>h}ihqFahPuBza+?9fmsClRcNP<y%4%Tc6!0FPi$DdFNm=vt1Ir4Q@K zpsPyA^kRbQvuy0oXjd6V;|zz{NM`fk9xcQrrVwSDpjx#YOI0VAd}<o?_yc18NM}r` zA)_=9a!<xWZ(X43NBCbXQl)8V!jNG_n4KajRbvDqYJ$sOjV>eRUNI2}7Z*<{DfO2+ zV2Y!%zTQ!TWv$<bF+aG>u_PtbpsEX0o539L_zb7ZFNXR`m&?IUkZdy?vwcfF(5A#4 zsH=8!@R196l2eFXvNIIS!YUpe$Wa$?R~Rn1Y_YTkU8>%LNOF){F7>OOsL$c9z#s-f zc_3isBg<#T%3)uy*f1qyvJoo)JL4-n8?(oWIhfArZ@92=_$<g%1L$0H&}7h02*+_~ zG9PSpWp`i3v49;89qVu#ndK;TaJq31xa?uISr_J(6y#1Ww5QIoG2Wp^TD7wV%ey|S z!KjsA*s{`8ECTImXf1O#56hfKnPyYF0(Fjs9-r7qq0fgM9=-et%k86q-Ck3p)`>q- zOPkH>U+S)Pl>6(TZplyOd|N6atkzRm6|BbaYU&)Z{NbF61-MDcm|BusJT0qw8MlSD zeAy16h3es1rDq}@7Q2_3(GYfjM7$sNAQZbI-QxpvHm<<AC38493aI{%4wzb!Ta-U* zroCutL6IFuY?bXkjP~xI#!XS(=KnKCaY*vI7i#fx*c&4v0z%Oh-Uj_7!_WkT9qkgs z>;RQU#DA(8{}@n1+=r(FXon+%ZhF<63I)_$n_W}uc6oRpWRCVfX0WfWRzufewoCLe zEWR+cniE=R_ZZ$*aSfkWg?U!!V9PYOK2EO-3vX@%LI)=1aSfaw+U2{vV!?_*-i(hQ zJ9L{1-~MMwY<VdhyPqO^$<(=pvr0>CnKiZkpx;;KtgUnhFT%RRWa(!CpK22`a)Hm~ zslW<4*H`Iv7gu@YZbp}xGR>XFEon`zS4f+Yjb)Lq4ok>-Tbj?k)Lh<~+wi`Nf<b2{ zdselx-l5Ll7S=ebmi{FDDS86qps7vc8Rp?rL|*D#b87!lYm^E(M@6(vn%OoHu0)@R zI4c^|TV_!4h$y@|z~Vq`F(PUTwOFY~j`xEcSD7J}9_f+PBGdZ*)^e$7eIrM7jnJa* zF!$9|4)cg-M-k7S#AbrgrWOcV_@^OMqhg=uZG(DM8@bd{^zaOwJm3p0?+*LGY*S9* zTEuk4VEuBC0I=t$n&m@F&~&aCWaJ#?w&Y4dEb#_iU>|mGY1N|3g=~tP#t;AOX-i|+ zMn$X&#{uR;NB^hF_Hw(~Q!{bwTF#;B$2NDKZ0^Cz_u1PG+J)>{(nPm?QSl71(eH2I zh}s0sQOklsVv{EG%wk4q+W*U$RIMb?$Ukxyl!c7|&Yi2}ZYtDUQj1c}SyrDkGxUjg zs&n~I)0@KE#~*D~rg?nCZ!P-Uicb*d>Jy~%f28N|1(4b_4i<ES#DP6KbX<)!F}r}; zg#M@DJ<rBNt>_&|Iyc`C;=}=-Xqv;{_oKpP0RIl!uMqtMZ2f+cRZ@mH5bS4KO$|+} zKUpMaVY1Y!A1qAzH0Gc_wA)fD+@ayA7+-iVGNu+5OwM!QT-unLm*+@h%A7fKaG!1< z&7EA}uqCI37<td9;fLa3$ykF=#)Ln74n}=Q<%j?BK2|)HuK0(dJRS-cPt^-H-18Hb zXqGec6wS<wG5Y^6R_610c=NwuNL8-Ni$Trdd)|dP!X+a){udOM7B2PQ%QZQ7;eio6 zT4GpnC3-w~Dr1=8;zMxwQ0gX)UVPdzpd*ZGjsB=cH_KbJ@<WUrctC{mQ*Z^(xDF`R zs6%<?z4h2)G}2IGMM&Nd<8vtwWFv5*2T$k38COewQ~r>U{Ow~caYs-Po+GhbXc&8d zhvE5@Zl-uXW*<{Lx3Zrp8uS2D)Oe67p07N_l+QilIksWcj45kmiU(KXnNlC7cz$Fw zQ+OzeDV}#3%M{O@SeZhPaZK^t&v>SIzUC~Zc;F?4=?pwi#1zlPWHH6_L=&0fIiMV- zc>XD$>011*gDIZ>DrAc1f{K|s@Vp>XJohq(DV~?Pgy}Vgv4ANa<f&j<hTje_#q%#- zrtn=2Q!D(>bQYcqVtNvuCu8b0j0UE7FlHrFJb<-|={P)N&Gb|}C&&~JD&4^Jbi-K7 zbi85Q#uN{5wJ|-@Fz#jgOFSpU6b}MD!t`uB=fgA^&)YFQ2hSxkML;~y6c60BGp)q) zeN0mgV=Gg1{cTLs@mw0y3_QQYG!xHJF`Zx-JD6tSxhtmGXn&^XqWzhkhxTVW5$(_P ze6&APJlHn7`Qy^YeX-qWdE??0n>SV&M$7ewdV71D-V9p0k75K_5GLLGA!E+PuWs(m z*nvptJ+%w}vP|!(j5L%ojGitO>phhbS4CP+hro=yDn@2cyTFXVDzvE4vtD3EVwIJ+ zO<+c3)p+7H0y8qJ;31=@L0|)#SHV)Fr$*pd;vC`y0>=?gBQ6#=p17DeN8m)_Im9Uf zk0xF~Y!x_(m<IJE3Ots$hS(6;N?cF8|KG?sbsTX6ahJg3iB}PK2z(at8sc_=Q;63R zuNOF-xQ)0?;4I>G#A^heNW7l7LEs$Xjl?wq=M%RRFA#Vd@mAtufeVQ{h;sxkChjCo z5x9)Fi`Xjg9O7=`M1e0M-cM`@yny&1@%|rZe+4l{N7yfL6)~Sp?&%QNOPomDE^rNT z67hP0gTz+iHi7Gj#}ls+_;TVD;s${mh_i@m1YSv;L%cxXRmA*BSWmIQR}&W#=Lo!p zcn)!jz&8*tAhrs;mbi*IQQ+H%Ylsbj+lcFl_p_gzdM|MUahJgBh*uGJ2>b}~8sc_= z*AuTLUN7(l;x^(ofj1JbBVHr$^Tg|k8w74A-bh>{@Jqz)#0v!8O1zc0Sm15M9mF{T zcMx|HrwF{AxQo~-a3^s$aiYLGi1!m40(TJ~B;Nl$?cYPp6-7^%z}>|BiC#~K!25_3 ziQ5I<Ps~Tad)5nlfY?gh28_50eQcjy+WcYTp8dsTB~^UBxrz@p|LVdjK8xAi2<8K8 zB0e0;S}VpKYB>LD%z{SUq~~73F9z>g1F&LmJ8n09KcsPQOm=7B7i(IBrx}gyG53(z zJGzNf!I<4y+r7oY@-Yw<d-43utxJs7=P}ZYu@ZURcWni8Q+sd}JE4J_;jhGwV_t8^ ztD9R3$HhYB-pdX_n~N?i{nd(notSo;|F7}9*kJN2`Qu_&;Ym~!;6&YB1dlS0PjffN zj_WQ$4||SMd)DBY%+{s?L>%*<LfVoS+t?o4Vzk~Ua?8fWHKSx}Q@-T64rRk7?iFad zq%Be9<5*tghRU;7u#OUx*n2|rTg_X-QZ7Ya^Zw@dNdE{9Z*5W%^M(8wvh1-+emv!i z+>rc<lINdD_x{|JucZ8&G|2C1#1!7zq{MxK?A9hxP?%UKMLhPeun5KboDz{FL@3sk zLR*Cvf=Sh_@=8>}r?HG+BP4to8`s*zL4)S*h5+V{6}fX{?q7vi!6P(N`Cp5Cg}(g` z)g_sxs=VFEqp5on!*#q)KIsxw?{ku`$4a$g$8jd;&Sd81=1wX$4^C*l{StJ&)}~6S z)=i>%f%3=2QSFIB?FEv)1aY?anC7>eUuhmYu6rFN9Blq)&kI;kwKfSR6lrOqP*!o( zK?+xe&tPF`jmkR)dBU3UVI^$HQwle+3aRjekfy9u3R`6Po^zo>^X;6uTbuk++pVxx zB5BQZVNI2!uY<K=Yu<x+Va*pwv_>!qYbcb~Jo7EZslpetu(U?yO+}uvMk(P!9;_*( z`=4eN(wd94HdDM4WceO1RA|1vLRcgEr3jK5srUYGbiGejdh_QwtxbYS=*{BPJ6Gtf z3g4sEtMaa8Sg`vXz>2xGDTnol8iTUN-=N0UCy8JRR)0qKJ_r@21jog)lG{WyRSAt% zsVr7?wu+YJEP-g*YnTD8xJIeUm!fipoI^BHy{!5kn8(~|k=r11Tanw^BzT0PD*s73 z1L`b28SBr)6))yserXwj@`?NfFVgsW%DFL#axk<|@-@goncg>Ao0OQ}Kv&hcDsK_> zmmN<jcOm&?$5VNeC11}2*oIKPY0oOD`CA|}*R*Gi%zXj5txbYQ*rW14rZCv^z)1ty zbEmN9Ap3!|=aH~IO3b5-7#XZ8?^fg~dz5lll26*B@@get&pX&&h3#1@H6I#bPn*pB z424>o1dp&s<sT#MdB-}CJsX5Q<0(hl^XITVO3bdW!}h4W*N~^|QOZ3<K5382yIb=0 zaJv|`XPwkM4Vs7BXT8k*g|J8P2zyk%oz6k~a4Xn<yloct#A)`t5w=H(IZCrf<$ZS` zY>!gzAIMYosJu5NU(W<=gTwZ0l$yIBI&4q7%$<qc)+WIt>{0o3bPnv9fGuwS_Pi(T z>DJ=y!>~O{%p~e5<4xtAsoA5HJ65wt<qc8xRAQSSwr8u<{CbeX_H@YHE0EjTBzS~9 zD*rwTgFThly7zC-ZedTSX3wW#dz2X8S7Ccp-aO=~cvH$vA)kylm6s~{dTz$?MQhUs zQsZ9=jRlgE(>`j_)G6sliWXAnhku%OT#qA(<|eNgW)<xXNK(#Lc`wjz4px)HHNK0y z7E!jvu<XG?xg`;fjK^qBH^*V56z={OYM5+rT<@F89Hlq?Psu%7aVqm)h#UUT2yvJ{ zHlU(z2iB!(DlGmIB0HO#wxbb35<4RlY!yUIN@B#0Rxz|9Sx7f=Me<81_wMIAZ|ByY zX}y`wZvdFTYNus9Ug5EG{+iY%E)f|Zx1pp&vXTTS-rCe9>6fDwTbjg-e51%KGv_`2 z*sFz$cS~b;{9cP#$~W8Bh6|=KBxF}m!cLSSv3F8ff%A~p+Vq92Zlq90AgRtVLY;0& z{|t21GRlO1(Rj8@aw1#C)KwIax?TdG_to{b)+QzODoT~Qs=R=dN?nz9*CMaA>07Cb z6U`@))b$ra*Zq>7OR2I`C<!yA)SlnrtWW7G3P@eoV-?<4*S6LsCG~yOs&rL(uP~4p z_32^lb|a7D#W?7M;1k=W6TdL?)OP7a%90u?s~(oJdb)5bNdb!apZm2a-N5xlqSjrS z?d!-~sMYrP=O{3H9i}s}QhM>T$ja}AyzbkfQM3J#7W*U3_VusW?`5a#vEp3x6+63& zvG<Y|`-XSB-|gCY+sotO*Dy+oKyfCL0e;)d9H+q>!W`$qWw~l-v~P%MykbL4!X51^ z+v{#=nR)NCe{6Zbs(I_TpKkgtzIp4WgC{n>bp{r#p7uoSUYlP#<HhEVruKxUkD*rM z755qmE8oW*jmG&NcN($Wk(OgI+WiVFhdr>WW%l~!H_q7fY~!w<=b55MnWB-|ytTP$ zB{X?o{JV4h#`o>_zHeW*>t8$Px9_xX+)>xQ!@Y6m{H;6f?K>?g=(xLH-FZ9DAXukQ zbM!QLc4reQJ1;SPoe2fVbp^QC%TFVReXAlLC&{-+@-W#n{82mEhoM#WUeID+*Vr9% zp8ejs^ynJj-e=1HgA}kuVihY}*z{`1<ylHgT#NnQ<`*}8cOpcIp_dD7u5`jN^oH3V zsK7T+D2$RIy+sU!)C0~o$g>v9W^<tHOS?_bIQaFYwUK7YGW|Kfp)bRo5e%z#+8cIa z)7MBB-ad^^wy$Ew?>=P)&xj<)Wz1+|1{-Qlm{YXR?F+)3Q^C2TRb(^?iI<EUy|YP_ zxMW;BTdB2aJ=kgHVz!NaJFV{Fi%(kY+avmeLhJz;<b@;}j=T7sO`Hczc?UtKyq`l} zYtvX9KeWu=E*fC=TGd6r=s~~Ohsv);<*hf0B50!48_BMio>fdD6A~YUHEKjsc{d_Y zj7Zqus+w1jFU{Npu>xvSt@z~TNb$~=`8{99!Q}bdr9aU8*rrh#l)Y;&EfJDLeud0$ z7WslJBL5nE4YTvMope-dlVaWnYn1XTuM>Gfc^oyU%3dO$94b}bddb%_XBcMPZu(xd zK)m=Dny!)>JrJeoBQpO3SSh)LrYfI5B;Pr|v$aWOOwqJdd8wL~s-hD#EmhtKrR5X& znyR&_19kUp`TL+@sO8N)@Byjhxk5=6i_E`3DB0E8q!{n&4!ct2wIWZsQiwpiEhC?F zrOI<lzMhdn$t3pNKH*vXiBytgM!@-4DZl&qyRUxGv$Hed=FP25if}B%t$aUW<zpC{ z(P3I<uWkIaL5w(~$8~>zQaeNC>dskywq$wnEXDGC1WUrrm?DD1cxG(KrOwLsgq1_M zQk%cMPvPH+!iY36VF`wvE5+=!)7~XIs3nKXFxk~HbF|^d<=t=o16@@XvbQz1kEexp z$ik#^BNdonR;Fa#UWGm+DsO9ialF!eh^Tx=6NyS%zbZ`YM}=9M3Xz36O}&dF3-Kiz z-DCb!p*fL-_)Gfkk3SwjJ<9@(JA9&cBN|i=+W&%CoWMigFz5d5(l776svH|GeY06x z5*G;^OVsG|F&2!F)@@1)=f7*_cPa-%_TO1cSc8`dd+;{6zTwq?C3rZz%-HE0Ve|Uo z>JwpY5fK)S;9rQsgEnuK=0-%AXsMb#12|oUWq&bo-av(JuTkz<+ZgUSzoYAn9Ze(- zs$hS0y}mnKDAWZeL*hVPpkAq<cY$6kt*HWu|Fa5(nhO7Q7dY-e73!}GeDu)(>J88Z zu453J>y{tU1#U2F>(d23!dkn(F7URQ|DU?R-eQbEJ1sSgw%cA7(lsnPmWHLY5WniV zs{BE;ShyGd3M#94LiO15k;nO@XD7zuM*H>{F$%10Y(I;m<#E`uWN%Km4x?p4(_8XD z%RGjhD5ngzFISP~F=W2X%|tG~Fd-5XvC6*``4r3@36D`$aw$6sIlXV<iz8XMfO^W^ zqsq&Wdh*)}RrN_4pUNAl@u?Nx+z-Qct2{gMgx%AjL8JZN7%ZpHv#(3IO_=&Rc9k^h zZf#SLh7v8a?^Ua|{4SH_F}Y$N4Kj<U75ixE7h<(tG-Td3S*A<MP-VKHye^|dV$<V= z(mn6VnoN(E$=sKa+u9^}gvV9>KIC^lEr=?M&qXNzs=U7=PpvCdMVmA}mG?)*cg%BK zLCY`2(4y=QN~rpVtwC#y9E|U5d9TK&^5$rK>f6J7^2u-iR9=eYqnte-Kwf97dw=6g z#@vf9b<A&d9@>nrI2-q6wb<hsbuD(W_p0M6b$+KrBccD!`CVjzkX{JcF|%vO6`kCG zJPx#T{{Efz{jX9IzARm_?|F8ME8=1jn(o8$Y{eDv{6;vTX)OR+^(DNNse1&RyOFOW zLes@H@XJ(&0Vu5DGsq+!a=e8FUS@#@cXHrK{3bRI+wfm(MHV^uI2t(kUa9V-&~hW+ zAf$^Vg(7dINPH25xvbj>J$U-H(^7BFPf0s14d$IVlk3KWN*PCfWvX>5Qp&iJ&DeR5 zfu#dP@%^*i$A(0$G6kF^R^a1M+}Z#y?6k=IB#|$~GoRwmlH&8E_+hXR!63ihJ01%2 zFpU+Ty-QTgw^6k={k9hXnN?trNFMx5hah|{Bk9wiZCKoE6d}d!LPf3it0?3Y3c;<E zP(lgz22zHjIJ6wi<t9>0e<F?NF@&rHg~*I%fIGz89^653qWm+EjIa0jMYNf*4jH1e z{9M-k3hNH3+f4pnn*2gF#cyualvU+QN=eugzd_!{8=_)~`p}NMgt6E1HWNN*wO^}@ z74EXHB4sQqm;L6qy{3Q1Baix${)n*aX5M}&<lU@<R0xl)A>|#2AJAi{N4RT^bk_o| zKke^RAyIfW3M+SsjAd*vo~9+W+C{}{h4w<-`D7k%`by9)m&|6C;K^3yXVGM+HwPJg zyhJrf{IxLfPTnR81MgG@284lYN%;m8ITe~BEvi}&n}sb8kne9o*aJ#Jy^zpG$_7ZF zi_pPDP|4sjgUY_vbWxjhQ5#)!@_TfVC|t?vWkiXL3z4B*Br2{$Mk~airgO<GT_k98 zC9~-w!E>xC4;Km3Gmz28MO1?l=LiEI<LzhiR`wdxMe9g;3#Ka<36UQmLq*X-VawCx zdss+#TDhoExM)2ozej<AU9`@0(R%5kF1o0<lU9qu`K(^LNMxkR>KR3%;u2)E+J&yi zky*M(&@v^n=_0}N74#G?hKq#hqrlV0MO1^jX9xrT!rRy6ZG?+9lCmBmm5YSPzsNdB z%zKyWy~5j7Lc%M`Mb`)ywUgpQ0llNCDB`%?YTszOs9m~<uYl{`@eaj`!Y8qM=^~Ml zD66N7M8)SIgJMwA0qCGyBxpaA%%+P3PX~CEi-hT4Afu0qs0JmD5(d7`+Xv)rgp0P4 z(#mShvkAKpd7G?*M5}1-cX;a-65dfRx>LBQgOnL4FtCfZnl9>)F5+ANx}SNQR*S;_ zhFIkyk?|QaR7VpP$3uLpUFiA-nPn6S+J4Dwx=8T+QI$s&3DdVBqmPTI26g{JG}Z^a zT`6q+KpALuv`$hMv0CXOA#$0j7Q`QfEuWBYvXJnJlJK-}Q5PxWQJ_B;RbXstwRf5> z>XI&Mql<3(JFOOl--lS`B9ZYLGL(x%#XX{8q3fe$mM#*sze;A)MS^FwDt~dSU6}qL zGWxiPYEbuWFc7-OTCfG@@`u`uX92zd6ok-<iWC<30}3d23GvgF3LyRscCC1ci#9F^ z&Glds#dFyP@$6sS)^B=y#R?CrK*1v!6R}2Y$%$>YAL7mI_~v4)ArA@8#(hb00eKsm zR2GignS=6cuv7$XXU=HteynyMk7~fV>s6CWpwYWt-L*;Ju>y}J4wZot+I^gPFI0$c znH|^oqFCNG<K*pnn6uM(Rp_A6li1Xx2l}T^Q}ojQ>1B$ZrqQ8oivh1}`e<lOhwvD4 z;I0=1h2P5Xo0y@xwv+=Zo2?M8hpAf5I;>8joI49`Yg%(59s<ungt8g<rN-<ywedO| zf_SoO-=k$-7s_M>GP3VMCccWADB%>q=FP44y;?c^eF6$=z64d6%7IRGPZwpQ&nSbJ zy<$HApP*MAKt<s)`&ecw%Cy?Ix6a-NUC~lqQJT**E5ioA+1}M^7iE6<cQ~3w#y(`Q z_feqnR;ZUYqf$I=+3G%s{m#Xgs?e6s^`p>v_V|RGc&r(ZI6r_LX7}T7DsRG>qg&=v z-Im#j)QVo>lXuv22>Ux*<{x_X#aQ|Dz2sF)(&NpFNk)7N6iU=uJCXeKA~cVMT+z1S znw2#E2!t+=mWJMkhc&ZbOt_(4v(9XP@qOpy4j#`Y;#*QgRB@~sCy+l5mbKbN?Q9m* z{u-PmY#nVj!~xyp7wE~T`WbE3mO9a%DDTC0V$GYSDWw3W*AfvX3}Dsl&FSO1XM##; z`$H{bdJ#I#H}A}TY5GR=ghjApQ}>A>{iV0jy2?dC*5BFq(&&q&3>qWMmj0%`#ca?z zdW|MOg9StO7Wj?_=^X!Lel<klp@OL`_R*ByY8MmRELN~ZvWy03F(~91-l4+<Wq@D= zIT4Z7_pBaA=VrBhl$fIe9nBVq=ZKBkRaa+>eFGl2sONj^GBF#O%wl|%Uu55oYQfNg znQCay-+4rf)`U7&q9C=Dnt#N+8(CqbikJUD6%$sbBCuG}-Vt`oD17^YpfVj3%#m6j z;<p$n5_~}rG9eh*!cc#LRG%*~4(YO+x(pZk$(6<iA%@wnlQB}ocMl0KR=1d|lKa4m z6IeD@5GUCkgZu9($fN(BqkXf{;E63WM>iM3-xy`!^H4ZK+0~NwTJVYyk`1H#+GbkQ zV&4N_dxXcNdq1b{nJ7<}3;!R{*N{YENWx01^7Xc`ueUKT<ZBfnPl?>YE!++sTKFyk zr`5iVinNA;vH?Ub3cetf=0&-9$@jv=l~VV~D6d?6CS5jYh?bfP?7g^U_O=%Lwt@Y7 zresWQweRm6$kL7$`~K#kdlRUM9A@u<VLYsogEs~SHHPNM*&X8r1^{z3m1E}KLTSq# zXfs&bGBXLSZVuh(nN8o}%e(zC2%%%iXCE!OuVYl&6Zak^ajS_NuShh+*IsoG$|z>v z#MxNXNP~JXyvlJ#XoY)mt~Q<oIDaYyim~KdLLj*UL5B*j3*>4@jjXJ`Zi?)3I7GvA zH9D9{O(5cSo>O4F73-E(yErI%1bJ*?%*?`zu=x6a^&N{=*}szIr0flhVEX}})@iXV zvkxRZW&a9pgZJ2HLaT@;|H{cEz9naLYcbA#vDSnCIH$*QRR2m$!i_IkG1Ho(cxc2m z$Gu5pu<w&~P84YW7W=m`e7<I1GlIjM-gnF3m9x`ZwDwJKM+NUQ2YY*0*}rYL;#(s) zQWkv;MfV;l0U!B5eT`S_lZ|)_gG@An=q{5XQ<VNK79bQZhgSY25cDi$YRoOD4PjHT z#r}Rw%S^=azE=AV^ka<siF8%ekdX=5LdN@$fthxP@DhUOI0)~#E8J&K;do&NeXCuJ z*5ZjW3@1Xzr!h`P*}Vx*@j{a+TkN|-n(P@^lSd#yXtLYXWRKFMmA)FY#l9EA$6k}R zP0=o+>-y_^jqP&vZZ7W4DX{hbsAOF85Cpqg=C^Q}FS-CX0VZG5Yf*fS`&9qPW`*JY zu^U5XSQ5IsxnEF1wJv7%6|<A^cSj<O546rnY(<cV7r=bLQg{KGqOH(J!~*GJ5Wvb9 zI~KQy=g^3=+1#HA+IJMq<2W(r@7MxE%*dDhyB1};Pkj{ujuX%BLtDe^(Hh)CdpSOd z?jsiT97*_Fs;D&!=N){x!x<~#Df}C^s#uP2y2ls466HzFBPO3(Q*a*#-Gtg-APw&n z60sI)-X(kEDDa5JM@Q!GMe#%^)>4UPSLV2DN)=-l$fgwLnM0rKIon^RU*cQ12Q-(P zD2cBG8uy*mQpY6{Pl~zJ>HZTR@WUBAmRUBeYQ+Luw8woI!f+J(PBSK3@Y74zf|qme zcdY#+Sl;;3Sz@Oe){EBCO0#i=pmJesgyt?}am9aKd%{!g>f2mPcQUm`RCL6&$8>D! z8p>)`qlHxfYH>^w4#1`T+B#mWCyx<@sMW@Y9(rhmC<mU%52*M#8h*GI{yiEz)IAGQ zS8V>_EQaaW?q9xyP0Vxeq3gCb^4~LF+0;9vsW&*z_Ki7MsP&5(e!VyRiCm3aVzU2T zyQ}e~vBEtqW5=O#)F*yVch^SzxUKPEZ*bIeyj#<H@uA+1mZ8Y0JJ4Doek3P&(M!6a z3x7vk18us!h~K&R%`m5|82&;~ewO*`tUkB9!W#5j12uR<+rnkma(}SO>a?oYLt7W( z&Dd7F0|<HQJ;7#Ppw0!xfc0E!E?$U=7uW}^!>gC6a>;m~ns_I&wR)K<ms}y=y>6|o z^OgC*T9J%9LzG|cuHXwDO#or%*zoGA63=3{&!OmiKe55=kZ>cr&~tW0&Dpk;LQgoe z)KksZvPppw;U#T`z24;siubjdZ;rk=b^bZ|eR2)?UTc#ofnmohuO&b3OrB7dU}E;1 zbMnDt%rWMgH**cVtOTza=ljKtK5tGo78c|=&NODuE-gzn<g1+xz6KiYV?9TGtfl;e z4tNcC*wKn~j@7Uh)m4|fYpwo;R`GUnwq?27%2yJ@185y%q_w!#Us;Qn&cg@i8rD;- z!>yO&Dsl70v6LJ1RJ+0X8zW#?Z7C_ILTWDL1f5=rHLPyFVBT7V#x8feeTW9p;E*1n zXGp1jZmfp+hZfE~7cckF-tuqojh#@=SzT?l<y+6Q*5ksvl6^rzvE}OwsfuBlRb;2; z8&lKEoaJyW-b&(_xePDCPM(?`4CUmfhG1bqsolUk!5q2X<Y0YJ+?IO8t@(N$S)4SP zk)K-WT<)-yJ5n=e=!{YV$0}S6NPdd}@86iYdVdSh0APM}q~Ggw6u=*b#6QgByDbbX z;PHA1^JNt5tN6t(X~;3S=Ip>Dg8aL5ZFo3k)z9!A1w3SGe)eC`e>yJ5%%5wfIoC#{ z4J(s#Nq7JCuOgNB?i<EbPd2s|RaS3$IPadjHs5u_=I`<jJ^1CDGs}xFxaPAb+rHhh zep10>&)@pbca~2&`qQ&l+;a0}=YM?Zaf#OKOV8`P;=Qjs{`u!~lNJ>J`J448p4(8> z?0!4DA-kaP`|GZ{y(>$VSMpdt-iYv%>#yoj^=98*a`%)Ee&45FMOXa&@kRN|FyL2J z?#8!I>pn91=t+aw*{3}0|N6cAW8b@X*u>%|tUp}$z~A#!{jJj9O1{XK?Wnik!?{=Q zcz?r9lmFLptuKE5;O6^B+fV-7)^gveb@{6Pf$|@J^St`sIXC8<``MR?m*#w%tMEDF zE`D#%EidNwPXBG{51lXOD*C3yi#*w{UzEFEmg`ANEGqgoDM#t4>{51a|8jNrp<^GI zc#YKSlqkJUKDgk?JI>0?i?CPozlzJV^TzJ|D09pO6L(*6(eM>bc*_Pq;W$(6B(fVO z=ZunZe3NXy<1YN<j>m?+Iq{LaH{Y0Fn>tbOY2}r^zaIVQ)DM@hydc_+q_gwOzCUSr zj+q`W(|=vC%l_wAC%$k>&Y6?nI`zsAKggeab4~XjpRdi?+mlg{_2t(&EjRq`t2rsx zo&Uxo<wf_6uRQ<N7Xnk>I$`hm2TPV`oqXsY=P5m|f9{>~G4r}Ec%t*PW6myV&3XP$ zOWtp~q%~*hvtPaR)z4dV4h^^6z9YNog8z#<!~5;fKTbS!PVFiCo;zc5Z0W80Ewh)M zfA`Scr*GeV#RZdY&i~H;>Z<c+ob#`z%bUMDKX}nquN<iIPMUk`_N<{_JbvycPxK@# zdvnrx^~<OK?$gKKIe#HPqb|!)?R-Y&Gv!;JUYysR)_D1%FYlFpHnsdQ|CZq^kWjgt zYunBHxjCvIn(aFM0y91Lf{iCM?fuopYh*t@Xr?*8ZU62!i=Gc&UnTSBN;#u*OnIXv zUXT;>^P%TeVW=$xVf$Ztd#f-8Z`#q@yB_J0J9~Qz@k4_*{(+q)($A2t0`J#I*CG85 z=~kr2;hExYq-P?H!#JIZbUf0@NT(q!LRy9NcSu(uy$0zzq_-m7iu75e-ALa<%0H;~ z9n$eg$8`7hb|5_q>3*cskS5|BEr>J)=_;hfNFV8een{U&+J<xpp0<7o>3F2wNHdYf z;Z(2~>3F0|kQO3cjns>DJ<_X@euH#9(rI{*yAx>%(t}8AkS1a7xe4h+r1vAe1gRH4 z@mPhl9_c!yE0JzRYV7ar-G}rjr15CY<B^_)bRyD1r1?m_NQ;qPjWh{6gbhf$@d#QM z(gR2hJeqSfe%570Iu>aT(rly)kY0?m0qH`dZAdRe+K!aJGwVV+?g0FRW1i_q(~&Mj zT86X{X%Oi(NN+&;1kw#i_aNPY^dM6HUB(me^TDy0XD>%ui1ZGmUZgvaUX64Q(npXU zM7j;>zme`kItpL1#UoC}BIVyBOhLL4X%5m(r1)8n_{1zPE-_**N0aa`82&syHg;9a zSWv9E=I=pUkt_G8xXDdJMq2V;1|2eym4nL;I+`?0&kfTTAp1<v@w2)xeQucEfb7#j z$I>B8pB1LJA$ttyFeOZ%8m6}+`-nm4U7-IH&-)Ejo=covpz|xXkbVBJd@JaKmCupo z2a>-4^w%K&L|y*HVfhXJf!_8X=<WZ3&T7%Ng`#nXnSY+rY+b}WXyb5IJlEUHzhM`< z_6SQH{tUQLkZ2iFXXGagt*bd=TH??lb-~0VhSb$34jWQ;`Os}+N5<rC&fSvxO76>f zp*Wa^v1T~F0iMnZRh$f&A`E4mRDnJje!y>L!}OxCA69{$4*Jh^`bA;-I?z)<KTD@q zh3Q*Cp9p$?al|sPN|fhU&#{HWEpbZ6xs9<`4jpQVp+%-R5!Xo=`^F)z&d}>~n+#Uu zWr3aydW1ex8;3f`ve_`ma|w7d!E-9^Lp}k|9LYn!G=L`yJa`Z{%u`_2LwXzNr-6=_ zeuU{4hV^I%-TEKu*aM!E!82NyQ5>#=efijE8P0GW{B7D;(CvD8*ilcG&1MWu1J5Mz zl;}Lu=-VNam1}8F4fv|ShaoL&k2apH0lftD{>BR{imV5{81&Qh^7QFuv-McVb}0i} zl5ia}Onb&7>j3DJA%lMxHq<Uc*P%}iRmy@ySyl{eJ23YR)W35;-#P%@it-CU-vIjM zxDV9@eTVo=g;~c+@Wf-@94H3v1)cNnKss<<{TTF<^*SyMx9=X%yFnjl43EWp+A)Cq zagaF<^yfj}s>?UqcB66^{mS|K64>!N`dG+T)31Mt?qBubH88giWYb#E`5U(W;($Ke z2zn3b1C{RtefuEg4}$&%=mYtTYm2uAk)H+np9i7yr*jVuK<Aio8`|Us(AVKUwn`h1 zW{5shG(tHx_ee38a4i!Lf1L*YkZ+^nej9j>0?$BWDAzwDKp)7b@d%J%gS5$6p#Kv( zpQ7s^20U{NC<Of*&}}+hw7*kXBl=Vjd~Eyv#%i|h4WJ)@oCsUX&3LBo*tXw+J|Fk7 zS8Kj2XtZ1zd)3e(Gg%#`Ck#)Nw+HgL?i{GCW3fIv6LhgrQ1im_aQ_(x`ZUm~cF6bI zn3WHD${_YsfIbd%G<mqZoNPq@2c7LQkU#DP{ok-RV$3wBR^iW=K<@^<Klv%Bl7E2q z>5kssE#Qk?tNE+c^w-cwhKC&ne(IEj$_nsI9`8doi-`;G_+%bFG0VlOB1Wwz$jhr{ z0A0p|HWAkX(B~ly$t{xYW{*$<=S(q=S<c2=U6xs}%y}lp(l*3+HwISTGQ_eQuT33d z`8><GYN+M5p~j6vdB0WsyKAWBPqD@~LoF?_#>!aUKN3so-^ITl##+$#ZnRi#Il_3! zVtLnMyeaOlJc2ye9bx(O2xIpVmQRKm*A25gFwAHh#{2EVNbMSCxjWAILmcm4i?e(X zXLO4D*5S;#bGT*4aO1Jzme+<Gn?=s55zN11gas?}r$$U#d8F}!xWDg6=68yJKR?oP zDBid-e$w`M<FE0&|2m%3btC!jpGR7@9c8>da?-s=8MhzB`<IVm{z38YH3^p65{>&4 zCLKsHJ`?vhCX%{U{JSgBvhQf)o5V>k9c{dDH19t?n$)$U`0t~mEci9(mQj;#_?glA zGu}V>Gv+spwzM2$Tr=A8_-NxHaevh@^O3*z7|Wx_8V?;i1ozh-XSpxQc;Gn8>&F>e z#r+*gmN$+!-b%83pJW^m_pcpq`D%>Od%Wd=F-F@M-tQRm%rF%0I@!{kX&gG)@~=~j zyH2q@bBgiADV7gUG2T1H^80k-p>h28#c`H9PBq>eXZdEF@rAg*`Bdg_7ys@%)$-D5 z#>&$yx1DC(cpC3tIE~bA#J|5k-O_Wq@!;u}O{W{riu;dFC-o2GH;q8DI-S*jApZS2 z-ST0E(V8*#nGEBJ4Bo$&LF(@^g%>lgiGjYKOt5UqGI}SBxCP;n#rtQoNd6?ta#OZ( zZMNml*~Zh^h{|{P4`TM1v3Ftud^^VSWsLDD!}Y2m_y?a}<D{75j>7#JmSD_ei(_bW z%!un_jC;s~5?2VG-k7))%LIdaXvn&8$eAlZ+7)AYX^8P98DAJ;!QYvFha#^YYWa)B zxPK^a|8J-TUywJ%TCR&VT4H(sTr7zjEtZdm8J!le{9v(seuVMZ5tdCy7|$Zuk1#$t zg2Y|JEUQKs-}9xvZE=?04>#V7vtY#BDehMfCvn9H%fs=EW5d`y!m{B=<LeQYtBy1p zj^zCxk0kNlc+208GB(A7<sb2ucSag(M_L{lY1}iC_Z=fieB~(1mx;z6I-ohh@?C=Q zT!Q7V3C0$2|5XBspC($`er8;CG*})v+H(6SW7pA^{YM);;{K*lBpw<i&QEv-?Zn^8 zHhpHENAn)e3S9@~{ZM2wf22&0lj$`wy+x+?%k*)XJ}=W(WXhLg;(y~~ncjk{wL zM>5TnX|7E9Gg{^^kZBvtB7R$@`IyIuG2)2SiTRoMGMVm>>0X)sAkzd`O}=q5&64SK znL1>;RHhHhbgNAHzG;@{d#Rc7{nbqQo@=IjzjQp((HmmpxF)>>Kc@Qn-@X2;SH#9a zbz}URSRUONbDGS*UA`8ZTNDxRtHx+8J}}wn#2b3yB4a%M;hcd-iw%A~u<!BMINt9- z%KadJ3B<V%Y8iOQ!9N!HXkI_-!f}kh{V}th0sp7ZHhd|m3e1Boi*PN+wHnv0xE{jw zEUwpZy@%^lT;Jgu@nme=7+h!KIv3Y;T=Q@(!nGXNYFxMCdI;CExL(8c9<EPueTQoV zPEf|+IuqBqxTfQphieh8<+xVkx)s+$xSqxJ8m{+neTwTlTqE#X#WA?f#C0yN>A2?M zT7+vkuGP41#q|)bXK}rT>pfhb;`$EPh>a+Z>r7nd;+l?Y9<D{Wmg8ED>sDOCH+awd zxz_Phi)No;O~)TYO|_<`q-5JtQYTo)m*6Y4Y0jXTfA)kk48D~9|7WI`gg24{%c}A1 zo{@~dxN4?V>JEQu-CbK_B>Vh9cQW4De|FGWDG~laRdPB1UiR6Z3L}{a-(*%9$ra0d zkY}dBT9dNG{Im5DzyVsV+e?1)rp6mIlEq(JO2%K~#J%{tjLEfr@fSFg-Bpf-_(Id| zsH#9og&-1t3DjNha@PbM=2w)-{EO}ZC}1XwTwq5f&T5Yf+<pkgKgpstQ!@Uvr;&_5 zJHo$ntNPi0lX8E;afV|<I>wZA0I~L!g-PS*6kHhc1m869O#{IStUlPrix2l3imwWM zRp8^kRecm+J~D9Z>1Vwh%dQ6BHVjl8XZRr>uOeI=&&bCyZVnhaupZ`^ha!T9W7dVZ zIL?ueW8!#Zax7FI^73;jE{<p9<5<>yj$zCLQPs<QUW;*Y{Np{xz_Bo_9rvu4yvmN{ zxW`gnl;Ie<1ANp+$s->-cOx#v$1!&w_&64mkK?=IyBhb@mGyAUADxc<4DQuO)q52( z6kk60^1-)25kccEuj@3v^%!T?XJI;@fK2(JuKcXUrQ~srxunM2KL{fJQKIjVeDq5_ z`0BxzM8o7q4f6K`Qy%%&;K#FTmf`yV0{Kz$HUKNWL*K{79r`u)<6?7VGAh15YkWK! zU3(?Yq7@M|b=x8N7z-OP=WJMQ`WH}q6yG~Qte1Rzv2;85l<s6ye7kV3`1t<lt>9BM zoo_b)+l=x!*KMOb*vmYAScdZV;ZpMW_aL^T!K397f2r|t?(PJil?CO8JiNZcrQ{vN zTe=Un_LV1!fQj`}Kfb5!5coRc1x5U0mA4k~Y1J6yJ5S3{w?~6t$dctIT#NG#MMR#u zWvVs1VZ9Ok0|%%lApdJzw1F|m5YpwPJchLw`lk9Qc?#<($1P(1T$~?{PfUl(nmG=c z*sOu$jfq*iIj)%4?6n*}Ol-EETJsDQBbzxMm~_D=)=nliI|BQ+iAOLI{n*4D$#B+> z{%T^g0oX51%oU!w#xpVAHGvQNnTZnwW`8m<J%_&wf_`ISwhH<(`iF_RqBHw}iOrFn z@owVLObo`giOm(GTB{8;@RnVC7^f!Pe4vE!XkzmlE5@CP@uFyc5MNA&8fa+kGt~Hb zC@z>CYkb*D53rvS$6)4R38h~QR;uHHm3}eCF-8V(m-N>VW27+ySiA`yX$(rwL(2EX zvt@~^kmfvN`iJ~_9I^f&=+A5QuK>P5`Y#5nWQ27@{eJ*Xk?lcwj{~a@<-GvF7ag!3 z*1ruY-*BMt2f%#20SzYo5Yj)q9UI5<Aj&%njC>7(Nl!2)0l$o&HJarUjA^*pqG3L! zzZKZDH^FcMZ_}_3xI@Fg2HvjW7T`_|bN{tN!w&#=X?O$h9u2nxcWZbX@IDQ{2fSay zp8y}w@B!e1z{*|=2a5w?d<<}$hED}f1XlJ=03NO3T;L>Nh=qR>j8fp^HM|gbjE0v3 z|6Ie@0*?h&^6v(=YWOkWaT<ORc)W(+1U^f{9|5Om_%q;i4Ic!~((n)rQ4=+s0Gy-Y zF~Ip6J{5SHhSPuxH9QHp7+CrDLf|<X_5&}_FrPoD(r_DajfVN7s(KB-3EZGze!sO! z!`}k0(eMa3cCCi3z-=1N0$!)#V&L@}b^&kH@UMZ}HT-+vts1@`xI@Dmfjc$)DsY#E zKLGC5@YlfmHOzCDgBm`OQ^wa}|0V;+Y508LL=6`LCu#T+V5^3`z~eRiYv2?Ow*Y5p z_z%E28h#LXnua$57i;(x;5i!R`Nskc{{y&6!`}ecXxM@&v0lSTzzrI<0k6_9&sElF zcpmUt4KD(2)9|l>*J*e)@Olm33cOLneBofbhF=EWs^MM09UA@?xKqOk*sOGE_$1(N z4Nm~xui=@%2Q|#+TnufzUIZMc;bp){8eR=-)$kp_<2C#|aEgZC2F}v(0pJ`Bk3z?u zrr{LeVhtAo&(ZJ#-~}440j|>UZ-Hwxd@XRjhVKM!(C}lxt2F!<;58cF4!l;w9|N~( z_<P`W8XkdpaJ`1Px7euR(}CMHoCUm9!!v+8H0%QI)NlZ}OT#U|-5S0Fc)y0920p0a z4q)TkaDVRxj??fDz=;|jg(YB;hED;uYB&veyoT+-DH@&+oTXvDx+zD)Ex^+>d>3%B zhW`XSN5d}zFVOHiz*QRl47f(a{GzK~!=tc-ZP4&3z^gQz4ZKFfGl17>_!8hY4POSl zPQy6LHr8wScHoT~ejK=6!<&J(YWOYS4h??-+^J!1w!1W(1l+CR(}DMEI2-t&hWWgo zaWLHemjK6U_%h%`4L1QNY4{JoRt>KQ9<Sk7fm1a6K5&+XzX8tCa6GmQ(=>c0aIuD` z0MF5|3wVKsmjG93_)6d!4c`P@ui=M*8#Mef@G1>|0=!1UBd~>8tKoBi+cZ24c%6n9 z0I%2Z<-i*?d?RqXhMxf5s$qO9XLM-z6W~q_kH8kHOT*)VyEQxoc)x}#fDdYz?^HLi z|H80kaDB;V?c+53XW&E)zYm<G;e){Zu7uysa=kwqTR1D|SY`-$cHlGB{$j{T;}YN} z<o?9Ki&v|oa09T~kHrj0GHwD^`?r`O#~OD6^Dk^cyb&{GwDAZqU&v<S1mhXtc>Kw( zlaX;U^hv`Oa{7m{aXdOBo)3J)9-Ko-d<SsMXR&eUR)T&EwwU}|Uwme56!^!OTTW1k z4P!zSo(<fF{TUTF8vCDxQS^o=d_V9N@E7^1&!fO<zs&lcjN*SO3cnBBwL3QMR9XH6 z?QoN<1;FbyedYo4{0k)Lf0VHRn19y{%OGK24e%rTV)>^Dgg#FIFNFNFB#tq+)Ig^F zqk-?}3Hxgq@Cop*^56Bq>pqQ*<8wN+4_n;Ohx(mEj2AF5KIpGJ;8#(<YOjq^{2jof zQQzg@r@S(A3lC)Ka}}`KuO}ME8#e&+H8OiJZ&Ch}z|X?pO21ct56S&~BKEg$0srjt zFnv0X1nQxW;(ydUA_&#@aunVjg%3vI@i<a=5c(_oZU^40$sdIy2EHjz`QtR;4(z8{ zA@%zu@OI7LtK^Xc=~i4fv;0S~acDN-pXW(O`(KTO_P?d^)BiCzkXQ|WsP-PK9chq$ zeiU8~%vaolB;tQ_6#bu3_$VB}e}n#{<WC2V2VRHrY`-UgyI{{&i6^k1Y4KCcfkb|& zPel~21OC<9_}&RL_Lm!hw;*1r2IW5n%(p611k3*=s(dUIJMH7xI1F2Y-wFJt79RmB zE<Y@PB`}{_W)|@Sz<T^V3T#EZO$3epej}><A>h^8_;WT63fs`$ihnu}4&{gP%Yi$= zulnERz;~m+V%Z|}eMjTBB5yA+UrdQxp<kl*9RcN8fp@{aG~Cla{9`>le>39^=P|(Q zHwcW;IZ@@SfUm@OSt8574*WdyLo*6_`+(aKZ#<u%{URuF&kyTA!HV%f8;{Nb-iY?X z$RNslf!j3yT@Abq@j-*A?-pSG3C4NAZ11;$yEOlO1*|@-kAGKpi`JgU1K*(e#|ErE z<ev$=L8D&^yb9yNX`r$GtAIB`f7L(NMwNde3U32mq_uZ1Fu(U;ouly0kA+T_AJ#t> z*bDoVy)D32&Ax|$gV0Cy@9n^G+IWz9ijrpfJ1+{C0iS^K%%(jXakEN}|HKEQ_(zP3 zls68z3-())&GwxDe3v$UJqr8~?4>P~_exax_oMLfoVetNWhMiEg!p2&Vfnf!`aQtk zX!RY>hLInZPm98nqi_lEWY}*5FXcZ2{2s<56~8g3D`_%bfgjQQb185#@FZFOIbb#4 zSTWyh0X|b}@BP5)L;XgM$9m;y)&BT<p1|ru`f}hgu%BjAzXyO<!d~T%9^i@SFK0;l z(YUz|Sn<yzL4GK&ItsS{H=#dRK%@K{qUaAs;Wwjj;+g30TKt~`yg=)(dBEz!`mZC# z_{wHreV+m6z&^XgzxpNCkC?CctUT#A0zZoOnJ?)d1FH|^eFGeX|5f}?!ob)Gc_Tn$ z`8$Evpg*Yb<V)ZUXkS%-(b@Q2i#FcP2UZ`-UlxVe06zeKsruTY=<UFnu#eAov;LxF z_(SWj)xd|k!|R8?0;`X(2lyJzo{P>=R5ShzS8Xs5tXsG+*=0B!d1bRo9QZwNnZp5; zAAzRlX`~9jqtfdycX}P-8+M1Yt{y+fuCDRoXXq8l8QE#+1_=&NMLkibVSa62QCD5P zOx0sA$`{p4DaoB_cPPYa)P1-X@pZnaMg1Z^*#(8?l?Ds591f3PJTl;zS7=MCywH5k zAg{{pS`3b&jJ(2W_ToZYhTH3Q1#A62hpWn2Tj{ogm~L}80{q#!%Rvbqr`NMw{D|Fw zpXEC$+;{|{PW;r~k>c@a!ZO@7-nsxK6lN?gE150sG8a$DEyPcuk(s@Crn7!_9_k^= zBGQ1zr+_U*FqT&NYlES3I`a5JNJyDofX5RcastbR?>#<wZ&qCB_f>`})_}*SfXyZ< zCd^DHtEovOHrXXTb@60RFi`BSoxCjQW{Hf&Ma8p9McLHF1$abciobTLv$jIyq=}r# zPS@g4R{G*nSeYAexqTJ*WxQa^5PA9Th|ILb`Tin*5WlQ<;fM9)%UWCz$n*MrqF!6t z;>q}#yQ^we`64u^kdq>_a(##{X-rL+p1F9oFW?renTv~S@jwf?G8Py37Wz&37ddPB zbA5u$#rE1-zi7Bjh;=&|MQG8|WdV&a1usnSf)EO?lxmOHiy!;L5e$J)zAc=eIuIc( zf{+3a@wEYcB!%;%B!vkPk}~?p$P_sdrOm8}a@mVLUVopmh!B+>AwbKIkdj&BtPQwb zk*4AKo7%pjvwcoq-`p&xuM!>5U#*MCSm>^;_UJk3<!C!SBU^a9pIlTOotwp)`()O5 zTzxWAeGz)3)<weV2$)t~9pNM32;Zc-BH3M0?5;?5TLgD%By(ycZ|dR*=Cs9;+`y5V z)I_vUX05Zr<BV{p&5Yo_gv<&zUVS0PQ5Egk0R(Ly32DgE8z{XxDkC*etJBjUSjV;$ zwi05=?~9O=**8ln+dl#A-wz>$t=%_2qsm<yT&}Ajx<%jIYImixPj0F&pcl18pw#LJ znC^Cz7U<j)@RkmAz;Zp$hJIS<R8!azzc&~`O&ySu$zD4!wbC<))J(SKz-4RP!2zim zi(o@UY^2R}U#RPx8VS=DN8zbaxGV}6N8&8n$+w;8b)-jUm}5h9ZtCL5c#A-()e$g# zsWTYi`%IsEY0w|pWm2l#HI5o*z~%Jn>ZD<$(XlPX8LajPY7jmVuF7;)*F<H>!0w-b z!0v}2dQIQ_ECy0^0159!XQoGIr1`5N+?%n`8P$_*<}}!EKD6wUFGl{px$KgCax<2? z@irdae;Gk{!0U|2sdUx_B6^qTwS98YyM49G^pcs2-QFdhh-OT$_D6<rcDdViSshv? zqE$0ob+w*|_G8#Z2v@-wkr|4`=v+7|TEGMs7NqVdTfKF7T2;3p&0ANm;|YkOJS-;k zJev)fWmW#V03=2bQb3r=MNDKcr`AMFw;92u{>YrnDyMg0WL8>@Cu07~%3b2GbXL^* z^$I7r99~alRfKlgE=S<9I%lmrLba?a2VUXib1l}3+iX?nPW=!vm>)5#+0sN-f1<P9 z{pH#+)u`vFb_N0tSA;dTbom;i{%9FKcb(W}^g~3yu8WeDhAiExwCcJP9cQ^~0vOnQ z5t%j&_tg%}5;a~dT6Ol+6g)5ulf=HIuGS4EUm(K6sWz7vCSk{8nuu3MRd~$3jGDz$ z(?ZtgVNZv_&Ea+XDuY#eO*Y$Pe^qTir7@^^YpR@%T9-G<Ya((S)vnaOG}~gwWXx&H ztMLM?h*~mA%ZeRkZjOh}APkBirjl69*Kb8r(iTs32Fj{Hs`7gyVky;D=kxGcbho40 z?et-j=L>j(vI`c@%F8V*v}Klh9CI9YN0fjFYS|#vxr3IS>lmbLp=YqNpbk_vBQH=e zU{OJgEIGAs*5uqmhoiuO%}UmVxg`a;lMC&sb8Pf5wouqsqs6iPsc`UZb4%QgnkDws zl)kxmz8^dA%&g4hN_Wsv<8lP6>U@in%j=nqO|`>Wi_y@5U0-Ad_TAj{N<C&xnNn(p zhExJ)5BYP8a%UFgp@cYu2ptFn{EjM&TV6M}Sf#-%M@hK@?_+Zplm<$J=`~JIZNZGv zKygu8QE_heh3f5ic6({imf|QCVQmWp3raH_$V&Gja!}G<>J&M4n+>l}tEme*eEtfz zy%>ekR0SeSaIt1vR#}00UXX`1DGm<M%WUaW-NAyqVn4zOc&=R>7|ei>2_YdKE~s!W zQZT5pcK2oWk{P9$@WW4FWR|1UvBc@EgMm~o<U(Pgqud#A8+>^kT<`UiV>o6rVHEWE za&o2?7EI1_q$a0;=&E-*%H5S79H%&guBx1zyg765;}C}}IgM&ZM^sU1AkRjt0#eMJ zIgYgCbg&8sA!us-OHn(b&Fl5MmPB*;mM$t@Vwl<(fhc9rAZca#xeuxm-Y*TNRXc0Y z#C;+ns`bOMQ=00q+wExxS3ByP$4Dy7!uu4(F;J=vcNJ)Jymzn;Es1vK>j&+nJf_fW zp2uD!e!3#P=nB+H%oa@zhm{zg9sgxq)YO<d{b;jnHm`rFyVg;T*8o=7OUw|f_LR)X zg8O*JmR4%}X`;@S%C{48y2pYyXqsV5De?P*v+9EAaud*b_>xLxU??mPQ_s#s^&T7o zA(k=zV-d$#%n-IBc4pQEf~mz-c+;jUkjj3ndhx7<I9NlU6|HZ^RL0bj+~R3j)yw#P zMO(g@kPiDul!;Ed*uBi`qTy)axD*<M(#xi$PC=v>&v2pDJQWrApJ8~|m9eG6d3xXs zFfx|f9VIe;IR4r))c8Bo>2p@PYX|E+4*m>Fc5m2VOD)MQ%AYmUUNp6!$PR=-9|u7` zj5p>IARMu_IgA+xe!vwje%M_QHE0gdB8R<NF(G-~3$;#D>h}ihp6KCG3{$og^WCe3 z{s0GCOo;)Hd2Agqh=C+Vwi(RT%cO#Ui2*M+m6Mw(&R+2o#Z^R*hdR7A5B(IuJOMq6 z9u?F30F_2`@l<wkv3i+Tgb6cL<nRnk=>gAjG3YretFg{DXInW053-)gH`@@i7<K`{ zKG7hWHau&kmGJD$+!T~x7m4K@#x9QwEh`5Mk3S%${_L7sysMaRvsKHN08XQ-ag-|t z%>E|tkVCMPcGbH114@n$;cO0CLU^>ArNNAa_!_0Ezy*FAzm>sS&WUds3S4$fPFQ_m zHG+)_CJd}Ys|&C^Ou@P>(-AF5Ecf~hH@2+F{AsE=tm7btZ`zd;Pe02y)hkf#|GZwq z;jidXhdnyP$P1^UOPXua(4=+vt4ayakOQs^Il`-HGn<olfe&+c7A8!5qlLp0dZEHw z+-LPx6k1=%z_n$XFMoD=U3{N3;w{8cDxx1fEwns$dBwsJD>Fvxgfeq72z3E3*6{yT z{J6>xJsA$nmwf{0$FIo{H*iV+L4#ceHP+^;dBR0l%tmS8sI0GdV2>N{<NeRrjyjg4 zM6Rga;%Ca78B^><_L8Y{3ul#<+A__9h%$bk8NA429_YXpY}3pIH3kQVf1&wW=>F{P zZzGdxLwgqF75H2pY>tDuzDlpVxXNP=2ocreec@@P!CFsEE)T?J`<6ljM~ORtb1B%D zmxmQTwhCO4WSZ}&E~vm(tHN$ibq4K=gln>~q4eQUq`uykhJ#dd^)I(beRf~rReGQG zJl5=4)y{f{8cwG#fg+U{<!}PV^A5LIp`+KRkuIHc9AhZ5C7k#!K<1bm9liO|^8G=! zHsWw$jiYL*Q6qQAIGuC3y;0>Rl+ic7N*Sm>!MZGJHE%B65sXpZ4$U#z7A3>Q@5e=W zhC8y&hA=9?L_u@<Y-;9;qXEu=|N97u?X(ynBf=-mjNWhri3v7xqdsVpn*#$^@#bEY zyDT0#FtQ?Q3I%GZd6pUOk!drXm8NRpJlk9|CwhHkHl*13Vm&axzE_S;Y2j8MWEwXI zk90Avg?G$FbhN0Z)-Aem<WL|em=WqdYT`$v%58VQ{W+~~e-@)Jd$e+wSXNd<_5}`I zhkc$&<rxdY%u&Q;f!OS`yZ09!7+oSqwjv(X3>pz(uZPwjhi9F}{iAs%;KK35;yzv7 zeAuAEYv(D&Pq7D6ZWBjSqK?6a&Ou_%o}X$uTa8Z$LQD;2pV2ntcq#+u*vm}bP;jPm z@)g6zVUG;%vBU|%fXfflvr(g)r{0`OLJJ79Q?XBBqzZ)rB3-RcIf8}r2V9<F2no9* za%fbpD@`%0@do=_!y;O-kS)iAQEskX%WPPHWQpG1?*`d6A<yf^SA<-t%F_tb{XeCz z`*((LHT8!#qo&BdGo)IQf!CbX54*m2t_wp()M%lGf&nHZNB>8`_Hx>UM#;olww!kM z+ZME?$s@?<X?to;AFz>)U5?Td<MEHR=^O#0c=L>)=bkxRI&JghOx&ErdXGy`kJ z875qcQO?_6H2-fqGDJYcN{jC7GrOn;?Dxo?jfDP%hDsAX=EpD5)%L9K%1G?5@Gw~= zz9<vxAr-zq^4K^FXLk6?ty-Sb>O+z`MB(vG@syHMf2qSB?GSFA%<T`(_G}Sndkn?M zUA-+!p7uD)%lz_B>6jb(h#no@l>fKGUTCaOEi%8ri&#o!=v$7dPVs~cnkGB^9fx_! zifO`8>&6qJTuYV@a3a9=E$zFVKkUOt_&EPZADw5Kr^oyzslW4q2@Vb^6J(_R$ne2b zCeP%5^stj1I_SXx7e_erj2|0tbJv`PZY6^JKQ5*6Y}{5ypLC}4+Yv~0Vn!3ofWyB+ zjOs)KIDF7$AdDQ~j9ScVhx-sOLwpg&`J!MJj_~k!tD`R97H4Z&*aZ4(ar$kunM=f6 ze!V4|Ai(8n^dNu(4wKWjRI$;V87|><SD!)SC)ig`2p<KRUzYIO(m}qo!Z`<&li#N3 z=M$-+CA=Id4!^Ql&^6?K@kOBiWsZ$!Q)T$xH{kNuxQqO?)i?+gUkx6{dD)gy;SNpB P#aQ5}4(77=ZQ%b0e7)k2 literal 55635 zcmch=3t&@4_CJ106DU}vV%4tl)hG`^r8enX5nI{-l~S@zVGDvGP17cjCM8Lsu(*~j zRZ6Y0D2OYtRZ(#lwXTS$tXm!;y6kRMRFsERP+CDnMMOpW|D3rq$>h=Qu0Q-oFLOWV z%$zxM=FFK#Zf<Wk7LJRFiBZI-j}og8YT2TROGA2QnaW91bV{5uLb+TStn>wmcxiRW zgNqeFAuVVcT(QIhPFoB<t_y(aA|A=6i+!jpu2l_+LKl@+5>)<b;MEsu<LH_Ifk~iK zy~OW>q!s*D@aT#IfiB`x9du1VB@<9LUHQOt37+X=f4O6xqBu@z;^=arEM0>CYw(dh z;v?i==W>=^=PDoNa=L43lYO3KgIX4<v9^sbn!?ILB-3dSQj(pBi)=&u{V&~nXaC>D ze)Zm+)7Fhy|Jsy86W7nOf_D=x9j*bmX0;BwB4(LVhub>~tml)AmK+^Wf}&ZHgYJN^ z#oOjRF+h7|T>b&gU%;;$JfQ9W`sm^=OuW#fYb!?H0r4?P9gLBXqw8~iN>WUmqDfNf z;_{Tlx-0r5*3~VoGsP<T@iF%+doECx)j4#U{(Z{&U!Z6bW|1yAxagCMD-YK=0%4yC z98W}A$nPg1oy_k|NJ$klF4ZzB<2OjA;hK(X2CkX7ZpCH6741_7q77FiE(fmJxaQz; z<MQJ2;}RcquQDH3Ehqezr;Cu@hU<46Q&NX(F|K-C^nq(Bcj9U!fX`iO!v6gO{$0WG zDv5r-bpHUqe-P<9TrC_w%+p7Bx*qAHxE|w}k|+56UyyFZb<>~RcjY$xaPLBoCeIpw z(sASpUFvA-+bhq#{{C||Kl;s*!*8_bx}JP{<-?z)-!ZrS%!!9)_uuNd%De8xxd|1y z{R*aPuPPn5Wc=~JTe~0g{pHgiEE`u>Y&L)K$P1VK<B^XqTpgb`wrc0lE0j}Du65@x z8(ww#FF#Iub*S^%+0Xs>!!TRbmT&a61#>o!HNJhk(Xi~-{ny=j!^2PB(VmoYFm?L= zza7=T=a_%Rk0<-bCw{bU|343}v;MZ}hO`%Ud)B;jqb~K0f*%ZJ`#<~ituM~|d`-pX z7n;93@BJ0&qaOW3@gF`N`M9s&i^(IMN#38HoOrlx-`vu{ZO`4eV9G;xzO@*7)-O_C zed+LuMMc}42cJK`?vZo-rfo-_xuZSCyztHQ^3R;;|7c}?+7%C6xHj(Jm;9@8t^LC{ zAL&efsABV~e<yvA`SfAK`il?d?>T(_yCaNm7hawGxX*LjZ8KNi@qhUv{*`q1%}+e{ z{GQI=9DVn^*u`I8a;<%K|GamcrZLNseSiJew9e0tFPeVtn-4!2JNZxjoo8m>{>MfC z?swzA-t4!^{hDv>sM{X+*Iy1FAMUF%Jnnnr>#vuL*?MDAZB;^r<;Ved-MuS!j=c6* zgQx%7&voCi=}X>?fg~wDaH|Z$Z5a)4{vmvBe}NRn|9cpJ^EsjX=YmkpPr-?QNZwr( zeE5g(@B0a)%D*NI=OaXj59Pl}E0CK1qA-2_99G{CVfeaWaYcuFa~R(q)}Bkk+Rp+e zmqGe8$8ms!a7I}9eqrtPhw0ldOnz6Gynlwtle&xP|6&-wGpG*JXLwlqio*1{F|7RN zFn!LEDx~(shv}~k<F5&;pW>JJu<~L0d>z)m`@`x#9M<37F#M-5{a*^hPldI2UaZhY z?XMTY<n0S<Us+iB1!3|Rgz3LHtbI3!wMQ2w|H`oPdxGk*{yPzte`6SaCTzU?E3Ce^ zG5MEuY2rxjE0uGVr_NWDw<Vl_1<)zr4V=F;Nl}_4e*b}3!yF2hza5K(x=XZiqmlo7 z<r2jJ`3I7;%xiVZuapbU!<q@@|BuVFUx<0;SK2sgCgs;(q$q7t`I|4s+<Hu1i;;fm zsNZ{B6Gz=b_<OXs?m4W-@FD-v4_1_`Kh(q-Ie#nsFAe%#&hdT7^f@rC>!|#A*k}3X zVELb*e>KV%bN+njvl{wgjl%S6xlB<$`8e3#B=rCKgF&oAJ`<OBCGM{}S5b~h_Sy;i zcfcNZae1FWzY`srI9g+oKF^>c-A9_ZSKx<)%P+wi673zs`O9Fhhmjv*8GkMexgGja zw^RN%ffs%f<PU|u>ps=Q4d?uKE>)CIkze3t)bUdLuZBLPAE;Qaec8Pd_U!Y4CT=J? zh~&>e`WWm7F$@nu|943C*#yJ3LO-a@_*bL<6sbQ)Lcd$l9%29A!X9nU2kjpV|7hX; zI|%h%0sFTC(^`x0QM7+E+D8hKKDWZZenk6(|E5B}R;fKTzfzQI;XitwKLhz`j{wNz zos0TLe<tvG%GK!4R;fR)LH!Av1pgexhxU9W_5aJre;@2YzDN43L&wy?9#jGO^Swyj zkSF}L2L9UyzK-)xLGc-=U-hqniayN#FQNXaXuoK0H_CTN{`3esVvS_~S7D#JOM~O@ zM#yU&7VJ;*qdn+<G5)YIQfj6CcpCb@j{X~leAJ#HXxI@+{=d+13Gh#JJL3;UN3Y`j zMfu-_J%*$FNX~x)_KJi3hj2U{`Yw{>-v@c)-qFO}3_kU55e9-^8n2@vZ>!WFHNd-p zQ@OmQuuqAU|3>t08`^&aens}*3;jFbAELdJVDA*kKNi4VOQrtW1pRU(`(1~Qc^&?? zo7b0(_;@eyl^hSmINgr=ujRN1`gqY^F}@DK|8}B&F<#cg-<G1h@YfV6J|-YOcB6gm zQat-D>~ltH|F^JDspK!2Xy40{ycdA|$S>M^4FwjdJ;UG+B`;{=zQVXBd*2KFR~*&E z{hN>PCy~7Z<%PdZfxag6muTM_$h!;uE#l2w<o^TmLrvxnU%`KtOZH4c!`)K6&xGSV z3oPun4gQlT_1`gI4g6F1-$vyB3FB=Mm%kYCEJunbW6_>du(y!+9@@1UcqY%UV9%T9 z@!GN7RV+8pvlQDaoj$+ao9D9ne0HBw?RR+X)^dx@>bBWkip5e{<#Ahleyi7Su_%`5 zX1~i{U^DymdYi-QwfMbOr{7m#Gn%dV$5>W8(OgzgVkt?nl;|x+a~2gT$}E{&HrsCV zyXPBctJv)KI^C6ddZTHg4U}@1X`VSGYq`s4ng@=1{%n>O=~GLl+Pyxf$DQY}+vcFC zlSE~enEf8FwbHKFt99hMTpn9!HTinuJl1Td+fU7g1&RwM^SdIupiO`sm{mw>vDIDf zsY*pPOydB7$wIwpB@+zeRQ=6*i?h5O|C2rw&HjveXsX9sP!5xn8;xl`=R!LaPz&ip z3Xv+~z19k@3*Dxc3s7wpw<TUX>n~C~*anlw=l9vX&T7$sNmf{u&M`==$+5uawD|%J zn#z=-2BC|Xsk2~EQZ!YxX<~q)FI5$mDAOV+)59sLldP30Eg%Thm2<-cXsP)gzfhve zSsbhbh9u3=3`f92$$8Yi5M2>rU}vq`N0Fx2r&%UhyjFLm-C`rB5pJOh6wQthXt4T` z6PlE;p}-AV?!paDh`5H6_!t75ns4_rOY)ooPO`e;m|ormsYRw-eR>!Vjbff*_Gej& z%P^Yk_5!od>`$w<I=uxG%|25RmCFuPVKidsrC7{%7i!e|`~~K83^RY4%kJ}`Bn?BN z81)%wy3_5mW7yWXopWpKWb$;@@liF=XN|3Ky2@3zvAT@&>~8-rrCUl+Hv|AQDA0I| z+f#10lv#atv(XYz`0QIJGpK=T9WtUJRrG<y=5hP{7ANAOm4-5zG@wmX^BMHTh?;pG zcZIXE289r(#yMSfquECTKWY(uW}eG#^)gHI(F%jZpYp8uN0kxkLi(&o7CZ=bpg_Uo zVvk1zFntPv$u<7!8bs<!yT1VKp@?qGM8KrjYffRlAZN_XxBL8F&jPjEvq?gMJ3@e3 z8(M3lMsnw#h!z%R38Epj*v>5It+Dxw?Jk&O9)hdWhG4*gWJqp(meFnZRxT*CmX&x) z9Cp7oZxYWKl24445MnQOhfTO(k59xwrl&q7NH<^rT3t0(zrDci^gH4BOcxQgs4h7> zv$oU8rk-$~m}Vx``22=^tKZ7Q0mGr;7?YVyHKdKV`w`BefX6+#qQFgl$IPova$(u^ zS;$_32;j4Ms_jJ{Zxz&bk`I{8=6tDUf>OoEV8SeqevtBTN@%uRx4Xt=^_ulk#_?8P zwbff5NeL~QKV{OSX(Gx=nYno}Ldsb)MUUvSN+MXSr^%&*hRE}HyyZ?e`iT_;Gp}%h z(KO1CoRXZDoXSn@^OxJMxdtFP1Ac}Pono<IJ>!{g@znS&o(gp`P}dw-S&%Q9eW(p% zo5pEKBr3Pty!I+YK~o4h5TJU_-$O)W1=8&9a#ytVld%}pq$2V;Z5B-ISc6rm>nNi+ z10&U9tF>CZcGehGkPS`3T@hP&@wqHSQ>R-rpl5}O%pRA|$Oa|5Bgcl>vIBvTPa#+h zq4qMXi(*Jw4;Dd4DW;_b)3priN;g)Kew%~IbYXgmSTsg!w5<fOINhRllG;7EM-R~g zlZn!2ie4?lLe)wOJ~4?pJw7&jq>)dFfMNCp%#&W>s_{8~hW*(Rl{XDt7-(1#=BAKJ z#i{}xm1*--VaTwVmrVoA#HMk@W{=r|rdTR#Yc17S)Oy^A`2Gc!c`1Po6+@s(b!MN_ ztypay7U(N&HVX}cWW8dU;-2q>HpO;dO_iM@0AIcnJB6@~4~BxtSi@rgS!#Uta>WLd z&Ec)Vkg9dUlPr{*FZ8RdsLx_AM-T&}%;!_-Ah}gvIqS<6JtkyKHEjMjl79u~Vj6L5 z1}5V~3m19{p9Og$fKF3`COLg((2oOCd2eegqkS&L0vd1_SZBMDSr)T}CL4Q?iyo>r zYie$BLGIW>qoGt!{ti7-tE|;n+(oPf!`6MlmYEG$0vc1%Tk2dMlsTPbs$FgK)mSQ= zZZ?<oy7&pojnKeotgaTT#Gk39UhneEw|gyRo*JlI{7V^MZ-9q+os|xM6~e2kbI9U{ z=2R@dRZjZ&;#|{&tf~dHDb(lleh4g6&(0b>1O70_zCiVc9{l~RJ0aT>$rkTnu+jRP zHe?nHMFBDXF#rw4xkdStCmD;z7Ze$R*hbmtMzpv8GG+=JHvhjFib9ghULpC*S?`Px z4+wZya1->4grQ8h9kq*v*&Zql@qdHr{|Kmo+-GM1NJk^RZh6I=3I)Vmn_ca-+njV9 zq(=Ln6YQ(1MCdxp?P6mLi!V&A>Vy{9I|esZw1&@fV4h_<=rh!fkJV+v!ke}Mf#VYO zum;8t?D1_bwqQk&SN)@xeQU3wTV41@pCz&3CE?in6d8-hPb-{kHtRF0y&k{EU1IfC z+Wo)5x<kzpWdU!{vl+R-ZF81m1)b}zblFV~C*RG;GE*klQ)xq5o$F%KCT3$<<gURI zvR0q!w$E3Wcj_iQa#7H0^Nd+l)>?}=eXFRpIOhK%<0*Usqo65G<LT<rQwYs4O`Y0* z)*7Xp&ru<LldASjh$%59LQaZ$wU+8sbU+kb9bj?5wiqEb1$xZv<ng|j!z$Iq(n2j_ zC{kM=Z7pX=jr1HgGy;pdvpiM}6z1X2mLfWNVw(x_HnBjE+&>kr8s__S+SaR8rN||g zqGu;`((ztkd3V+WvtAg5*50ZqdK;I$IDj$Vpq39TLDOi(z&+<I_a&_q*b=YT1@>9D zmQ*c<T%b*%)A-q+J@u*NwqYJC+_8uG&=UPj*;r;&M`{KRUCSs`{oL-(<K5j``G~z; zuS3X~#hYk1&NfYC8~td3BWx2`DwYMk_$D=(&MYRHQ~w`lQn8XiC;!Y{P!=`<ICZY# zyQ#ovF%%ipSyrAj)8&bHymjF(Gn#_i$DiG*4E6YkzO#t770+a+>Y3d5Khtvf`bTUU zdvm&8{J@wUIIPB+n1+DZg#M@So~Nfnt?(U58g0HI#EJtv)-`8;?1#Ba5BBYKT*3Q$ z*!uk<tE6;xAQ)v@H3X*BU(A!UFj-2~_vR)Mo!P4mjrx>wdmucS@OAfWWqe`5*gOl) zrIqn{d6ra4iBejM`!oe<?$`p0J~=gzk?!eP_*g}Q|LD0G`9sgas2}3ehyJB|jY5yk ziTtrBk4M5GM|?zXkwYy<+@K_t7K5kVqy4AhRC*d-{jVQVK^3$Z)EvC0yI_uB$<Q4C z3kpgLmipi19>BvTc&tRx;Yx5i@kB;H#Yzvr;X^Y2F7Zvn>50X@${dNmQR1uRwNm*$ z%KLaYgyg5-@{e2rq)~h*Pqeq5o`VjERc;T+>!Td!^1wzuC!WlSQ<m}ks{B3y`CBg5 z#`T3DdIn8970>zr_rtR(dnl!6V<^QlD|;zLckZJUH6EZ84$?s>p6NS6DKzV(6g3{F z6c4SOpcML?q!iDNoTd~OIzuU*bx}}nKRk1yp%i+=QHp1N;wi<mH3^jBVV6WoCn7tg zcqV2rrFfP|M=72G8cHdieHuaOYCMBRDW2U*p%l*qrBP~8lq^c|%*$v>>DePn@5b*w zD8)lO6DTdg?-MA+vo9t}VY?Dab$Awt(#d#+jM6I<Wfr9_MJcBg55+hr#lu%FN{8av zMoNdF-zddHO0|?;g=exT9ib?7l;Yv7C6r!`XO$?u2G9IZiid!fQ#uOIyiuBrXKg6G z4$q)b3WsQ+6c5|2r?gU0Hc)E7ZvrUAz;C5AO;I*cnyx5YD9yk#UX*6y*)2-5@XQpY z+30^tuSfq=dIS2O($VODN^eB}Q;LV$rZj$Ru0O8nBp<tF@s1W79#g&dOm}y8!<&9> z=Xvm;&hOwZoqHj>^p@9lbf>qW$GeB^#J^Ogdl-3|gKoPzP^^0xd7OiIUF{5$=Q+s5 zyV@8g4|JeMm9AEX$rBwq!YvGwM><9jUd}LirUMpIy6PAv4|PBrrK_4@@>E9-;aLom z$2ukuHZe?|>o5_{VVFGFQA#+4Ve(|hEW$d5$)g?Apsoal6<EqqO;}-AL%5dkiJw3i z7Du>_a0kQjgqIR-XE=fIa>8v4ClX#wxRv1~!Yzbb7#>V`1L5Ti>j<|Ju48y8;Vp!# z86H8njqognuO+;ju!-Rm!tI1}7)~R+k8ldZS%f<X>lhwQxRY=K!#RXc5LOt@Cw!Xl zi65!`69^+ZqWugP5~k;pyV@Bx5l$f7#&8MYB*Lu>mlD<yZee%^;Sq$FGdzoM3gJ42 z%L!)@u4dRlIEU~ohFyf|=U`nXhN}si2<I^DCtON6h2dJlvk2=LUPRbIIDz3h!qtQo zhL;eoC47R$$*`q_>j-x+d>7%RgxeWjPIx)tHiqvbyqa(;!>b9m5N=_34dD%hmowZ# zxRr1n!|MreAzaPy2EuKGXED5y@NU8;hFb}@6V73H6XAV?QyAVtxP!2c;jM%_2`4bz zM)(9_h2foqPZK`z1GRrQVOmdgbuhe#F#SBQtDWI?!U=@i7~V^m9{=uYWq2Q99pM&W z{b#o}epr9xgsG(1LC-Tg=pp7`Pj%3<m!0<@!^!0#A0EHzZgk<{Nu?&Y>B4H<Y}tUH z365Hif3}_O(_Gq}{k~7vycs^L)NhYzWLdfsR}x>%pZ&h~lU-U85mUcYbIZ&fji;Ko zBD%BqSbsd}p^eDe(B>aR<4(a%|5r6biPqiz+K%PTg+nz^xqI#zX!Dz?=3g&9z7JDv z<3H+e(D;*I%^#|9;K@^z=s;6Dr^AXwqo=nUHA6e6V^qCFV!P;>@#cmCcp1^RBW=pl z)VFDxl;)L;S~4`Q5ha@&@_C+hC>t!%%FxZ6_j^%3j><D?pgfHPs^eCaICfd%TaCMd zQf7nJc%t!r;(rc<H#Z206i=EPCi1dZ3i<IQpHTzyNAo;CAw71fDqlz$z~y%>#k3t1 z_cgLMH?V@NiG^ImOWy@Wh^k*D9#ri5KPEzCUCy)>$Sc(?Xt#=jTd{PYPWS>H)7(IT zf|~mU1Q7KMqn2{&drU0L!!#507<6QFgFvy!bu3A3s-S)Gy{hN2gkT-}K@%ns^}fpU zbqyo6G(#!+caA0Ej>dhYSOtvGxON7HTXRDtSL<P>7DM?%<4EnnpHpLIasKb2-mwcC z-)Vfcaq!U2O(fxT<5ykVv54ZCP^77Wgi;mvo+ja<a5@#{tr4^<KqIRk8xd6EM$m-9 zf1oP3!p}jPXr)m2d|tlmCaBQ3mS*eb1`pSEFIvkHZ_U-LH4e_-MCJ{&<`al#tvQ)Q ztzns1Ye*=y=9O<roG4r=^^KsJK@+VJN_ap+YYNHyJE#iYno6n9MBY4JzH2U2Xk1&) zTEoUAbCPPV_nA|&-pNAm(M)faiRn$nN$;6VZ&CPZsa`>QklccN|2wRbn;Y0jW;ObG zjmuDD^K*pI6s-7^&eIUen&Ka-p-MI|Z*nk=MX7wMYLxJnJE#P_<(Q&6w8(XtsC*GB zCzrd4by6*_`e`(esFN7Aj#D21wYh=iVTua+HZlg(nTut4Q^MlyIha$L2B3Tb{bCns zd@0FUnM86BT1c`R6qM<Hqq#wdSq5E&s|wmY(w|%v15YSdM)`Ow60{<ouWKW=9wc9F z&r+`W2VfIbZO?K}eGSy+29}4lN6^0{VQ9}o*lt8`&+V)|r)eDU_S_Y0j}Y@Bc?=(V zg7z3_qCG;n6_k&+N6>!9^K~7-mMYku)m-y|A?<15)bCKJxq;<j?Gf}VczfQ!b}D*% z?qls4L2`I|)&<)m#C-ZquswqI4rro1Lb>gfkGDtAp5pns;zQcAfopyfG!OPqE2oZR z?O}OXdj$R0WE}KQOpon(gtaG5YR@yl_6RYTNbM1{zEXRHa$lVcwnxwofyM^!Xl#9h z{j-H@?tz$Kd)hd48mP?;EDvjspf9FM(H{EgUiAKXk+r8&^0#fl_6RXIlb+n)1TBN> zN#jQ-ccs)GK|4>h$AztZusyrE=6?b^*q(MyZ3MNsf#qTC5%gzB7}{gSq#nIJJ6L=6 zN$q((*d8Hf;n%_T2wDYb!rz2)#gvcxo1l&1`MTEN$fCL71Fmr<)0iPXIqe}$8uoGi z#jJ+}>Q6$@j%G~ejSVgqW(DscBne{++H2Hq3RYu-HGT|Q6D!-KXpdlV+>`)A#v_`O znd1;CnYrJE8frFF*!`wxj?kOz&)n^wIEnc$#C3&=!~C&{=8tJKU-PMO#+MMeud!h- zIw2r&Ux<RujERLRi`d~R2Ig?@SCB+2k`yTS-scC`(srKGS~Z<n2QYK#L9KF16L)at z^5%x4+}s;bk|SQp7pyxvIDaWRps9h)$SWDGM4k8gphKEy^{%b&j6bGjYxzdw>gEPf z@FsEzJ`@D)8txLu#sn1@4_b4>7reS(F?AT?>U_i0>E!(Hz!yCunxN%UyHXRVXH;ET z0j}#**gR6#HO&n|>b)eD>ndosgT_1wx(e+cWZHepb*V!4H8(KCb^Vd)dV=#ONxBLN z)40^GJ8`}zbY%s&u7Ad2CQ{dy<_014C~6hD3ff!bNaXctTpx^>FF>Psq4^@XUApEA zl_s`J*N`keY(=Y{=d!v!#z`az5Sf2EA$jR0T3;kc!==%<fihn!)%Lf~QDDjjOlNGR zbm+6t%AbO~&L^Nzqj6)Cabu&g^;P3~8kAkba2EQik%o(MY(|rD(|euobsSvt$_Uss zh>}83oQk-EuX%;yH1Y;8#kpWvS~b)gH^tQ7wkc-dy0#^4H4ijRTL0pcP5T{<yTAQ( z+xPK}ySJUbyz#A(7^%*-1ngcLUmv->vAv;fV8h2ytNyn2%D^T2aYv2D(J$>Zj?j*D z1}4|ex6pF52d!$F(%SgO$ZapyAG$P86+K854PxW&#)c)(WIy}&;LI)ijqCRtHyrxz z;LNsz#w`bG+78&a9GtoPpt0?sHVy;#&}#?R(zyfGStJ>{1eQJ6K%9d!R9laQ0+eee za#8p+P!8j6!5+fd3po1%-Zk_icFKMMdga(nO~wuNoiR5U*Vm+l*Z9$KRsL_efHyf7 zzVaq@y$HE<e$o`zWL)34ecSh!Llg_WwE3o$PSA(mP{#u)(1(>53;Cgt-oyeS=>cOa zlxGE|eKpWU@>Ve(1;0q%ABm?;Q0+OhE;7RdAsLn)G}ax&rtcmy;o1pgWaCmI-2EvL z=v;{BxQz%6M4*l;4d!J1vvyW6rx`g9G&4d2lQ?5&;=u-1V#d&T>ZRs}R%EAU&ZoXH z?xj|DCE`njCga|a@gNY500em^i5iZ(_=62J52*4^gHQ6Vg1qL2!8kZ*nzEO5z?9Wu zi2l%pad8}#--XJXSF$4LqUM#9U1ZuNGBF|~K8Mzbh$LtafyN>d_P3(uX3CeUZi1)+ z(q@6k$2UhJua(ogzK=tbXYS?pK=)&tM#`Y<lgFqfOcJA)bNWg~XSqV?_u{LUgKG|w zp_&^+=F@16P+riEfX0-^icM73PWgDK6tt~8U)Rijm~}hJ_N)iu*}u?qDc9)PFikgd z`Z2VU=VF=)dKdcn;LLr^4T4ZCX(?!9BrQcnS4&z7+C@Uk=ke84b3;4oj_mnopkbir z)pN0vT*+}vNh%gf&u2<@G&hKhPj&`PDQFLXCQQjhpx+uOA2+3-&E@&J1~Da*Xyit? z=ZvFVNn(=+WMaa9@ALOwd%x?0j&%d?-_hJ4I9EX4lKlghJd4<j5z{nfb^WJxEaW5( z?c9fA2ZIG`QWhrjJln4oc{YXS8F)Wt34cF2DpvSvXGz<@C4FeIHgj)8>3OU)oQ=&{ zEW^PiYyvxI>|kS5n?x%yKH4#P(ECa{_kV)X$_p7=>f1(8BS%9P=AD(KK&D!mkhRu< zVZ<tLsoy?AsQw+676+9B4a5}EqN<>aqC%;pLTI6Vs@@|)3(<QuI*0tHLMfqz=*RP& zpMTtgdfF^0ckobc9lDgq?eEYmn!^LuP^bUB+%E6ECJY-a-J;f(MC$|!Cn5kHz*+;+ zI#p>wdb>*BCkzbPZ&NKn4PGJGRgA;w2nU;wrIX}fs`Etpm?>Je2p?++@v)#64`QWz z?cf~Aknl0qSM!hbV0Ph{(Y#`8PleW23j?jL4-TIBWIT8vF?&@oT3zoP4i*Xwf=o#4 zX%N&3734v15{nE{ftdfZ3T8=#|9TKy_MZwx8w6b+^`Kr4gWv&jV|DTJGX}wfYHbmN z;9YfQj2Hy}Q0f0^5F9gM1RT^(Aa7gq3X?8jCRM^*S|Gm|yQ2JQ^jL5lmO^DQS%@*4 z4;oD>T?Y}3>y3M3SR7bg-*zp<%XmClk-cN!3dGBS4gcb&TIyNkXg+UH|I$)YJ&Vle z)Ehy?w<d(xY%J(@&`I!xum+u_ETQ%6AW*vB#J5Mhu!;2KTS!40&Gn>{A5rxPDW9Mv zO8LaHZ`Oywb_-e|Xsq2gLxXzb`WURKZ!m5cxP~?Lb!;uEQNNdV2G^rR)0FjM>2^zp zn&lq8Y`+w2Dk7HcmvXxhM!aBOvGwD$jl9fJE<=<#3gu-P?{chKypWscLtc|=@e)q` z2dK>rEDy7|pnnZ|=gW*KNGDLOu&<yU08K0|L`A!$e1f(~<h$@CT1NA4#?YfQ9)wWw zHQSw1XXGHlQ_uUPe1c|?@`<kyCsIEC6`-JH@_Z!c*awieui1X0ey1|+mRl_|o2_Sd z;Cs&c<5^9{c=Ea?Bin$@#)1+HK#mfj|G}9Zlm$Zi0v(vtap1Onv<uk+kq2j<IA}ca z8cD)8ri+i?MZ@B@xR`+r4`4O5__p}&Zbcc`a2Ei2bv@q1)OkLPyM^9CgszL5Mc=5R za^g@8Cn1w?pI>RIz@t>4^&kbF&aX6a*o*()29PKRJ%xrG^fIZ=-$Ab}^x7a<BymXO zLnIR424ODioC-bY2<)KNq|Q%C2eqZ@9Xpuo#FR=hj*%my6730fSK5pU8uvZY0g(_z zFP-f?Cm?E;D&ShS3?GVOx;j|lpqA5<7@dhHI*E^mSQ0;ki@yvlgfsB(_;l1&I8h_z z(+;tU=@nGX#^;IR1($C5xzE)Q(QIrbPAQZON_n43aw(rfKI%)LB_z_o`6Hb!P+_Qo zX0bWC{s4(*@=l{Bd{;-`DyxK}q!^v`Oy;_fKWteZ)|DiN)ENl{nvH*_+n4z**LgJ4 zc?)sYQyuDpTy5U7ybfY!_Nw+4t~OiI?b%P2F_F2F+KiAQsm;`x3_^t3B!(*WqS9vL zt909ssr{-j^+cw28*x4ZN2tw2o&kYOL(DpukgN?~V%FwXk5jGP#?`K)h3JJ$Z6<OV ziR9L1gauqAS(~Y|8iWY7Neor`66?2r)2)fM;NL>+sZ8zN#L1vKxV4!`Bd>#)uR)id z)!xn3ZX>mCeV=T^M81W}gxZX-9R%SdOq~xwh)|ovP^B8C_PccZ0Ker<Vr6Q#6X$lS zgR9L%-p%VE<{+-N@NHpj8p_SacCI$P1FZA$_egCfaxjVH)@ForxJa@#Q|D?BBGe`^ zROuwvZy(U@NBown?PY4SZ~J$nK|*aN@?8*wliUIm%0n_#?R{KrdKFveiG8Ft6KN-r zTx~`u;UY<Grj8GU2(?KJ>AH&b+fll`hN*p2ILR`mb_a3Jr#iT`nMfV4gP7}K!f0xj z!`GXQ9bD}?QahWe%|t$h%7ogCuoeX7B)2phnL1lRh)|ovP^BNJm1?y3hJivpNxy*# z(3=jSO<)!w4i$I`1w<2=_?f&uVtxp{7Vo5qil#wz4ord+4NX-LPd(9j?RVYXY*xd} z$MVeFv7F78378n0ax{&`Gjua0zR`q<^9;+12aNcv`4sHTOUE6|fdDEC-oczi>3*<u zACG#F^U!N*E`|~hy>@6D!-E+fOgK;mN=Wx{>OE5-zG+Ha{dP8|H{!_QUbN?+g6Er< z{!mfnCsAJ!-xr-fLGaDd`6Ys%D)FIhlLD)3`zSW1o!N~jFxPg*q4Q`ugN~K8B{|?y z*A&9`&`+Xl=;oFT8rO=>Ra-L^kG$r<>!=&(ICn~%SZ$7iAUabr9+8ME0z|5SycJbG z0wNCKM{_t1uyIGT@t9N&Z(=85jWbY%svP)K_9|93{Ak8M_f_Lb*aYL@Bq|D)IZkE9 zqfE1LZ}XJn&=oz=5vKVhwKD8#8jT&zMpov>pMvhp2**JXfn?A-6dhPw8c`{pGi$b= z#v<>QTZM0P<9bl&24nodRamn7<KgIgvB>G%@wTugjG5RplhkdRl0aIKmFSrxEVNkr z_chHt^V)U|f7Y1i6`8oj8$~AW@l8-DLF(-U%1>5=DKwDFMoF+{A&nk+CzFRu!|20f zYT4Td-q$9zPVIm8O;5>NbSjmA?}y+~CXMPRkUtbHYc{glsasI{>o68;Yogi_C*>yp zjVu{fKY5$3sfP6@%DeDwK;sVH6f=NowFI~cIk4#V#x%?p;F7cxfu2#V2pwmtcj~xR zZNo;wBD7;$=j8$YxwX-|!bE<me_#F1#9O!wY7A>Ow>Rl)q7JPgt5K6*z?7K13$~;4 z!>wxjiCzU!_&RlblQEH`^9gV=Rj`W}O$2CCDD#DMhXLp34uS}BIXo-!h!Q7;YPIwb zB}D}?H1$9{g>2L=hC2D!Ddd5Px_-dY8greRnZ;NB72aA@iwsSe2V=Y5rgJIuCe*nD z1xZV;`A0-sNfm~wxadbxao`dI9E(aC+k<AIb6+@>Y8ZbG*ZK?{-==Wz1%W^l{GmMz z^)tEp^t#r7E{92%{!Bl<>e$4@5c%JfF;vC(&M+%hHK{9<N0ArDPSjccQK0p1zs-<G z_IruiHw7J@&@?Hru@Lq~l!46y?g(X<@x1pSFN=`WF|w`oYHON|M__9wvlut;=cIcE z%9F{N{rg5XB!M+#;1Zp%^`4-u_Yf^$YvCcA8MXHo?nNGI;d^kLX5$`Gq&eV}bzn+f z@C9+Hd0{4=_X9I=CD(l{$_o=;O(yHrMa`-L$8KqwvZu+or)T?K%`+OBjVB@<nYW|K zc%m`(*g(>RhuKHaFxsK>;Eli_VrUMZ-4QPk0Muy8W9G3!-j;RfGqkp8QWAPy4c!=- z4d3HyixV+$p^GS=F_Gt9L0)N0ICdV#%_^?HEkR+Y{WVW;FQ&1Hv!}2~gL)8Nd7NQd z;hrroMi2wWA5Q{VEcuoo5MO~{Kn2$Yd^IE@D^*`Jjt@DU9ALUi^k<MJ5V01R;w@XZ zG#lCOcq3@k#h96y7182*v2Y8IMN5re^KxAFCh}n8NucHlnx-iy2W~cg4YR>|G-d*; zh*#TbGKp`>No+RZzzb_V*pKFP4Mp{@*(6-QQ-_&WjpBibtH!-Cl)-qM*EyP@6T6Jx z#?X@m#~a}s>hyk?J1?J|-lEpNiQHk%`@-q&?xn_W>u&p2@t@0!zK){D&f$Q*>_&a{ zw;ku5_y9Q>>jXAj#zH15eLEH)B%Fs<`kMzBSs+T8n@}6<Sg^~uKc;CC{P=jY@c_m# z;(h{|DlBAVKsJ-HA2Kl09$;32^IQz!T@MAv>=23<s?#?cS+r))(jc5L9iK;>4r})| zJiiD{#_cj54rp?uXH7Oj0@LKMs>u<dNi*5%l3m7Q2p`8(-X6i5OQwr9_UhaC>Rnxb zs#D;f{zXd0HTHqCi)DTjE%U_+8;fA-#QXmlmaqP}7$4bGVQ_pLM#v0G!f;pj3qq(g z#MH5(4l?XlR5w1^T$<1fCl4-w>CrM~1rS+ZVT`Z^(k);hE2(!8ZsE^?h_gdI>S4U^ zNi-exu!em88XBT{J|Ev+l<myi0|&?PQ0IAc!%shc|Bvo&XwTyl8$N77PmzRvITE&J zp?L=<+B9Pg+>C$YmYR5klP!MR>tU9ZETZNUYYN)PK{uxM7f6FUg#@gH8V~W&I0$)I z=VKt#ujBA69M)0^>QD~slvpT>1zVLuHL0PGkDO0-lU?GQXb-3^H_>vOgVrCvwyB1e zNOY`5OP$WY(xXi{c*HVGk5w&OV6*;s6d?@fw(mA#vPFKf(%r~QbMN=mNvUXg{myII zRyC*>wU%0`b}kcC#u^*&Kq_kkiB|k8+6HcJ6JKFzJ99{Dctv|mTTJ`5j#$dN20bhs zQ1WA9;f#>CUs}hr_2gBo5NWmL&wu{&09Fop=mXl!9nJpXO7QnLuuJD;$XdK320y5* zKR&o~{Hxfyy!1W>Y;!&RJMz_S-F+Im{X_Mq)DR(7EvnOX-}fx#s^1lp{gd}l{m#M6 zC{2ThqDs^Ut7{%wgCFzMpYHZw@DklsEU!PN>27a|1*zs_J^M*eAa5shq2EgSuSWkY zg4a&N?;>fUilJXY@y`NJjm~Yim+Sl<ov#}2psHA)EA#jrI;&2+{#I9kH~;GJmJ`s# zn^0Anuf~RqKHc@YT)grTFX;B^`d2Lw<&yEnDE8J@UDX0nF1eh)TUzI>ahG_IwLBSj z3M;?RUQVy9Qvrlsr^m~Aik)-pZj0d48(0-02ZSrxh0ams)uZ$&h0Y++?5v^}PjP`9 z;Wbo>vDW7Fvo}(yZ^pgFF!Q?n2&%&0+pBUpP>gt`FVBxVHBW%VG7-7-x_o3(N|kBq z%`^qCW58>a>HWz{#JgdYih?}L)ykwP<`RR#UtX)wi*C_By6eP8XQmg{<2A=YLkre* zIz?AhQ&nd7>O2)X_V#7!%QCx;Ub+ekpm&sWbtbQ;(u<dB!v@zYx?#Hhx<$B(?euC) zlIwR?*^%>iicit$Q&NUOYA)pXtu7L)=<M_YZJh(1U1qnt;SH?AAw5LTfKpLzbc*^S z`RlL8>kFi}m@D*7LMUggs?zE6bvNj0ap7&SZa?Gb^W_W%kzs+3kqzM5Q!SGg;Pt-A z<J0`sG8maC`Gx?*i;gV0u4I3$U!hyGlijNC9N{I2#|im{O6x+4zRY6Cm?&rD5-2Ds z_@@rZ?Fel@VUFs)9iR?C{m6*jgR2Opl!wT56g-JWzuQ-^T*pu6)gRi^7{#v<i3f^8 zE6jF$bF~HEz_#Eal%*Hox8C@1i~2cnQ}lFbj!K`VraAY7q;*S@bBS;N=RaA>@x%9B zQZ-iDT~t}M?Js$cJhbDX`*wVvcjk#N-<(uty6NuEo@@DbSL>L9r?>v`t9KWUIsem9 ziyydu?u{ScdU1j-`_>!w-S+-B?O(loebOw`%coi|zrN1VXn!ZWF1w)chZT3O?Z^`4 zg*>VszuN!B^%wPsdb8ISKRoV(KSb0k_#%I_d|~!d6xdaiTlvoL&U2E_AJbbqBg&)w z-@N~r=Kb~kMw|Ys`*FqNZ|8~no4LJ(d`9Q}DEHrAa_>5@f77b5|7*GC?Vmrf<FQ2J zm7nXI9vxPbFY51!{>+;<)c)SOCFlCjzD&3^=i6L?uN!*H`$rzwp4)x%?S>!sZO;|_ zZF6Qjv;X~@+*V!=50V!ZeVddc^c3w9?cDq2vd%LXJwEzwuGf$-y{<ex>$!E;X5@vm zSF*qG%TalQkA0MJ$xWjV-}amSiyQC;0r~{}O!N~Y*Nx2?#Qk^`@4t(u9$oiz?3<%E z=DqpGOs|1{pN3DcywLZziBFCHaN&}h!rPHFD!=51EBfcC>92VD-A#v#FTXbWZ$omf z9{bj?J3jm%f9(C$olkD{<{axvFUb1xo1CWme*bl8%8DD`*jQHd=!nW2U;CSH+*_9& zyYY1K!mKOLd~$=(<KCCvExTlT$4!6TH~hj;#mzZepPjeAVMcS#{1?C8`SqpEIcNIo z*B;1jxal8pBVFIdJ~{f#b>1PzUm7`9WB%g_?UV&KJ{)`as=bGAyJ^h*`QLk9TYBTf z>%M!wtnvFB{lB^M)sqg_m}!69n-%-TGuI#eYuCU9Z;rX4cHzyx|Mcm1Z>*rtpt(7s zpGRiAP`2y&IeDF_^^0bIxt`lu)$-^3Yr;1kfpR(bw5j*ga>O`P`}O9V)b#qBwp`Y5 z?AKfF=HvLZn&#Zz_Wkc>Z}s2n;Ph!+PGXKKFOlO}IWd>U-rzu}Ed=BE_U`T!#Nb08 zbay+Dj_By_ZpY!kM|cG9B+{RdE=AtHcvf;f(s-n8NJk<)fOIs{(?};EO+>sdL7IYe z7Sckb_aSv4Z9%#O=|-e2NcSM!iu5>AI!lbl(|9M5UW;@I^caJ*1?hC8Tan&@bRW_` zB0Y)pWjuwNgwwQ6q-jVK@l>e^>GeomNGBj&inI*rdZbH{wjo`E^Z?TBNKYe;!-J}c zI4hivG!1DvQWMhOBXuF&fOIL+ElAfR{T&`<--&b?(hj7nke)$02){l^!d#_8nuhcm zq$Z?ANL@&$B3+8K4C#8LS7Sf06Y2MO%H{~tzWC7z&4WXb4n~@cGz;lOq$NnJkX9qD zMS2&~dZZhW-i>r8Qu<X`C(;>6HF%P?7O4*DN~EKawjeD<x)Z4%=_#c5Asv99N3|kN zLAnR&-AGR$ZAD6d_waM1Baog#nvZk<9)OyKGy&-%q(hOeMw*Ot6VfcCBQRf2K$?#f zKR;ogn1!ZdC1z0^Oa*VhG+v`w8bdrCF87h{Zm5nwY<phZ*oHpmqFY!w3P97rBOJ~v z1t}jq{9Z7E=Lg9M9+nFcyyYN`1`n!5@HT>UZ4bQN;0^A9cLcn6@aRjka6M^>Gqi_t zgTYJeftLlI2E5Co$SeUbt_NN<c%)0Ta(9835QT@|D6$XwsE~CM{E<SE+Dmw^54b~d zZF;%8oBqa@raeX*hd)!S6eMW-)F}A_V{5AWPDqIDQ{zwQ+oz^Bp<kbxMaa+k9<jlV zFKTb58p8T_FnBF^91Oqbi{NE}w+g&Vqwq>N@7gH5YVaOIx$wT%A#W$Hb>Pvr!I~u+ zZJf|&TD^w8av*(MQGOk4G!pmW`tAX5BX}Y8j<1iUa-HC91#g&Kj`K8-y$d`%a~dI2 z2i_L&@bk9_-e~YP_aL(ryf*L>qm=V=nIXC$AL+RqyzMAwl*>U+@_~)u9S5&i=1mCN zW)Jw6;ek;M^9WsOtsW0vwDG7xygm=S;c_`{8y$E_;9+?mB9rtP4Ias)zZ)6o3rU|D z;Qtl#LA0@43*N&$@a_X|1$gA60l8eCR`9+8FIr#i0q;l;GCRS03%o03nYbsLYY@QR z2QONT(}A}uicHpi@HT@-u_vG>Z$I&8fPX9IKl(d&;e4uh5%}j}-il_Q)!_YvIWSs( z(KlM%;6*Ff4&LWIlsgXI7vM#+0R{eVdXPC3ydypEa=?2l3QvdjZGp{c9^Zg_O^Y<1 zCkB1C4&}$e2E&mjU?bkwYrr#s7cItY;k;<}*$dvUd$8*X@N&Q#BG-p|>aTdj)d}F? zRboM1B!9jZd|HQu$1dv23E&ljH(V}@dmS!X*USKKChj#mB|8+<Ywyt9$>I~OJys!4 zw7#MB$Qtn2LRQ2?$fA1Nz<U|IXuf*@yg!2%P2bbtJq#XH3$~TZOoXtzdMK9y-V*Sl z#h60yyisK8AeYu;rnkGhUqq^Dmuz4T+JM$gpI|Pcz6kgNZdnXP8nckielcR+qxMr? zS|f(HgXmef&cOesBMrzaQtKlUPB@;{j=~Q(7iimJlow*OEq#>tVj$^(KH9^5lmmUV zUuP+I#%kBZDl22@emDE~P^|V@jdChh+oVx2OW}T_hPZFDe?Qb{*JzcMS}mT%-Ko{S zr&Zo$_jmNAJS+NYKkcg=?yEi8Pg&7V`*=U4r61kz?MK{>e%goQln3MJ{`ENR2XV?i zcHi8eD1Yj&J<wlyy1(}I{>lzUSvr8|>jr2yo}+9YFlNa)%8%^+(Q}BtkNx}kIodPv z${q1z_Qos!jHml=;)%QAT>AIrbG3WUQ{Fju4E{#w+VklCmGg*xn*DqCK<%0Y<*|Wd zP7YK)WA`f)h`XEpdniGB{Cwq9!kC@sD}Otm?ms@CxT^=zzfTR)ZoNR+HE0a}+G_I! zbpPB1M6XNKHeIOPov3{#QTa2wzw^SGpdY(X`_x6spD*fz`+F|dKANOFezEr77t_la zalbA}`^K-7x019!Bq=A^{p-Kdetn73{VVO`mnbcl(EWi+Ug(Fy9an1a&QQ)=sr_z< z^3V|N3qzE@4$*!%M0tORc2%14=b`lP_MzH!!<6@jYEKPSzF_zF4<q_s_V4jw+MUCd zCBwCAhAS(F)BWFu6ZaJR_YYTTyRK56xJtY2D&<9X|It;%eQ?CC0Z8siqw0^af4@)D zev+=NNFV%iy7GKF-G7iy+&eSK693G&I|j+|Ozo~LWl7e6by>>&S#<yREMgze(mt52 z+?TE0k*#dahF2b-f8eth4*nA+ly_pZU&bh(lDpp72mg>iD_6u^d>-ycYW*=|wU*e% zm;v|3D34Gcl(>!M>5hp@(Pk>N7lW*AF;_1J>rjk#XCLKD%J{cF+GVl$le)OAkJY}c zQP#%d7M}lptnwBkJg?DyqgD3Otl6#6{z0oO(-Prxt#)HSrKvB-oBL{a_f@v_CBmcq zw1?xA?exCP&-!V*`zfax;Xs`Bj{}tN=x-FS@2|yQc!ytuaPI)^8|Nsk13=zCK>OJM z<r7AD{T%Im_&dz!fV?7JyFOl78&8B2@oX<m=aDz#@4>3=Q_m*p9t{gzr};g7u|)K9 zd3rHV@8;<PJbjF(&+v3BPhaIJy#x;bD<AXpPgq|P9*MOOr5QZU<*A9Mvv}Hqrcl0j zc$$y6NEp$9r540z!XBQ!&(kiR{=n07(Q3+fB~LSWI)SG%d0NBMCwRJxr}SP`Dvt=l z(!cYR-sej6<4EI?(%++|H5dKe>TiDP_T05t{rl1*mTPF~pp-VKf1h_1HZ~~2-ph)} z%06havK4QJgo%_9_=n;-o+?%7sJQbPO&s0tLrVJ${P<C!`3h<&_<Gxu1UkCc12ZYu zntNV;S~`mX?Z5rB;+rXvcPiRafomSFJ8-SS^%$<rxOU-s2iGUKPT}hNoF?ufT*Gl? z;2MW(Dy|A#^KjjPYZb1?aBar53)efiKEZVgS6^(@F2Xe&R|c+exTfN&z%>uo9k^EE zdJNZQT)S|+gX<Gqr*QSfPjfE9H5^w4u5q}g;;O(k57!;IR^fUK*JfP1aJ_@;6I`co z^~K*EzX;cGTp75=;hKu80@plTci>uui`fQmQ@&m|Vtmn*k-9Ye$-fkxAtfa{B|RlW zH=-C{t4y%^Rr;tjY`^|rGV4{znTGUG$@)=-NXfrk=H45a>|0QUFK?7&{HZ`Sb%;Cs zg<!k4T1j?${PtwLY<iU6TFDXqvPyCp{h8NM&T=K05Wb*tD9Pmu+>ocHey_@zr~b9& z5MTk%Yj;t8^`_e8SCZME2}#DEk;6UvgJQ{E5Bt+`$##dO0^d{FEsk=O6bK^mC+zIC zHhZ<-qJ9IFOn->j2L;q*whh@)iM7gULv9ZQ<DX<!n<^Q9lTJy-`|{}zG>W0~A5&VZ zQ{14jo`x8a2Ef*zyd?4HGXxhx7|S;S`6hs&0~Q}_9oUD~;v$~|`5eeMl=CPN`SNj( zBP|crOR<U0QfNO+af3dTkFFwI6rU&`#jH|fpz~CUVJO1#Pz;)ii{cjLqZl^=M2c<V zLwV_QD=vypl#gQ5%hzF#0H&yy=yc7&Me&U8DW)Z%VQsjldMU4H$3ooWc#oB#n79x5 zNFO1O^3h<e$0hPn4DCcdimjB7;<d<k7w$<{s)u5DLYn#xZ9w)B_1+02^5r034)V<s zOz`NIt`$<g7Q~GfEV9=iulSIz^jVEd$fG%9bd|bRXH5JfM6cud$S(9Z=1Y+;ff~j? zB4|GbO!6q7;}ra09^NlafPaL%O~4}G-qV`6y|vf_2qEB!d@oD+=(wzXG4}b?;rt`= zwF8qpVy(s;vbsUFFV7Z`_pX$W-i^8*`J`$C`3_0>=q0Edkx$e`GR6I2+*6-Xy)<{V zl03ARX!N0cB>y-rq2CtFv0KoQLIzmup757aKAMZSBA<>5@(<;q>w8>6-T}Pe=D^ZO zd8`OBQT?PJ{cf}a`P#|x_(#W6E%Hg#D3p(${ShB=dp-yxOME_WjGmA3OWm&Qy$Ak+ z!@0{q{{|Pefqal$NS1fx39O|sHpNHOE3lkW+#*cBcBBu*CrV>^%@l`JtkyvBM#WUS z8dp@TWtifJiq+l|YnWIT+0=NT@>w>v7E-Y~5NLd>nC-gRI94%5GDV^Bs$#VRXk4n8 zoszLNn~L!Y0eonjshFMWvGJs0b@bA>Q8D!j`W53t#q1=QjRO^{k)8Zq#fg+C<ku=z zSAb&e6|3M)k@%3Gs(kff0rE!`CsCr1->Dcc^1_GwOU3AD=@YA58sb|T<;!le0F6`1 z7K53GN(lX8u<{%MEcA;}E>zNi+qu2^DCa5@f!Vu^k;b6(bfolFWa?SM4x}k~pE447 z!?_Xl`$2D!>R$}Jhubd(t6sQkNc|52C-D9tdCvfg54GoS0Q5Uu@<ghC4^sMpuD~Av z)30+$A>yAwy8InY9G$a}ylat>e#NWu2P$KLSD(=E@&lC#zz;~69=KiuthRTcVg+uI zup4;2gntXXLBdVI8zoHpsa6R;4!lXin}D}SxD9x#g!cfqN%(!>of19@yj#L2f!l#a zd$kA_`y_lJaEFA40e4E}Gl7pwI2ZT?FvP;X0~IsymlCc3{z}3NfxnjUJ-{b{g*_hz zJ}u#=fzL>IJFxO)F#nst8VP>{94Fz=fa4`h-}NR)xDP^9qJ#$mCrS7c;K33e2CS2C zD)3MVj{zP5EbKcKI7Px9;4BF*1<sLh3-AO9)6Zi}5`Gi7RKiDqXG!>5V26YUz_8U4 z)&bW_I19K=!Y1IQ61D*^m+)_aS4;R0z%3Gf40wZt>5RKo!mj~uk?;q=Z4&+lc(;V< zq@!KJmjmyUa58X*gl`1ylyD*N2?^784yPqd-?b~!`2Q_%oP?Wz6D0f~aFT?d0M<!( z3-AaDzY3fpVLH#qlJF<MITAhvJVC-*#3_@6lYmPltOuSYVLB&qNO(GMwS;E_*Gl-e zz;zN{2E0_le*|7G;lBc}mhdaUEfPKiyg|a>0=G(dAT}#oBzy&Mn}jogcT0E@aJz(O z0PmA9y&0lI!V7>;NO&3WX$h|bR=y67&#l055`G6bLBc12lO#L{16wCyI`1DL;UeG^ z3C{w~k}y33kt5;TfhS1#9$=G%{{&nr;irLTN%-%;4hioCu9on}z_k+o0k}@W1F(cz zD&ZvH<r2OMc(sJHfLkOy5qN`yZNRM(_5p8^a1(Hwgx3M@mhkhy?GkPW-Y4Ns;0_7@ z2;3>*L0AHwknj-T(-KalCF03o{~LkhBs>#1LBch_NfK@X)=Bsw;1Lpj7C1%1uK;IB zm|iWEBViojDib73UmTevJP1qJQV9<Mo+aUIV26Y!0#{3T25_x}=K|MB_)g%Z5?%|u zT*A))ua@u*;1&tL1-wDRM}b=<Oq=a35>5hclkio*yCs|r+%Dll;C&LF0o)<sxxk$g zZU8<Z;Rk_FOSlzS`6k%^uK~wNct3E0giiq{NjM%`2Azbj1|A{dalk1OwgG2Jcph+$ zgzo^JAmLTOCJFxqxKzTg0MC-}QDBFJ2Ve_RE#d2cYb87ZxK6^efR{>m5%6*euLNE# z;lBd6NEqLUDH|kw6u4Ew1F(hKBH<ChZ4w>_yj#NM!0i(LJ@7sW({t1v5`G!DQ^Na! zPe}MQFn#wy-*?h_KM`BF)8J#7q0oBO2s}XSFZ!IT%m7}+_a{o93zVubTn8-nV=;Y_ zlvTiD{}$8dBIQrO^ao-fUWw_GsB8qL_X(+Zpz;FnsbiYBD?zvt`lMnDnRZAMM@MIb zX9Az-(8SSBk?=a;V_lj!3@e5@nOUJX6RGzk1l)Ga5tgcT(-3{L^B!~TpEI3N3; ziZFg%7=8@60`{W(RR2@JV!urFJr|aLXBgfO{5td>#>-zO9d1(A0^sG6KGT8e{0l7T zf1WZ6nEo~omO-q2)xfLJetM39^!Y2W3G%PyI7Z!41F7v#1YY_{&|V9G<6vK5zk7jK z9M#0pGbbbuTin2h^t%q2-l#<WLH5cE<6D5Ag+8MH_J;9~1ON1~CT<bRll*FROV8>c z`HujL{d<D)E2R~f-VkyG^BC3l2Jm{=U+8-nSnTf;u;2Xx_{UF!{0bZqlu7>ax_V>~ zsPDrtd@2kN$C1Llu$O4hgTQY~@)B^wKrea{_81D>iv2WIK>CaXep#~j9qN$;t6ztI zR|0<y`=i^KJzfN+_S^-`+AHNJ`*Z{Q(SFh2mq<q%#J?d7*8<aPn!#fJzb%Y^A`Hjl z0KNnE6Y?hjpF(>#pgi^0v%s&TJ-azh$2dD9^+yo}68<55%EGV@cpT0bZU&FW%e}y@ z@K;iU<Ua{aF9jnJRQ_*a<@-RfgpV|F2wN<_1$e7e|6Ed>f2jNt!uTGHNQ55+mi=WT z@!?;ik&wOr9ajE(V80Z9uE9Yey}n80pFjtP{6q3*0dGfsG5!_;--YptWed~yZ7IJF zv@T$JuNrQde&<Qw5s<vgfOn&PskkTmqyy9Wo9b@^mC?ZB`Ev5nsbS?SfU7^$#1-@M z{{(&v`k@<{yib8!;J<V}LGAx3toⅆ)N8CMgc#J{z7D6<!1xeN%mU?yczyM4I+Ja z0@I%nn+{C<y%(6?BPHzjC9wEVd$m^}zto=>1232Ca~-hwQ2vROU*gXI_G11R4j$Ft z0K8fnpR2;kw}#=@flH+R{t=kI_n<ll;hUdWGBW>A{g(hwMf*g18;LKq?+M^rp^q5f zZvuaX_$T65$`Bz<?f={`TnrqC@<gWgY{re3$A7}7!t%!r4V8B#a2wjM1DX0a9k@w~ zUmJn%Lwl($B=425^6!P=i)rHGA1adrya)b8!-mTH!uSsXcS!YJOdZ2NR9+v3bHcC* z_*%4IkGv%RdElLhPr`rI2j*0PByaFlm=7iU%m7XT9>dGO04(Mg9p;yvz==|Se-11@ zB=4LNuphn$CmLZhu=o(J1=gVb)NIo4QD8UPE9`L;csRz(NX{RKn?=APe<?Behvdx( z!;Qdhj0Y`vB!5L1|M4)~9){1m8t1W+|6dNABaPQwVDX{)mji!__E0xaeaC>W1~zg$ z`5LSfKhVU{GuFhv7kD-LXC~(#1{NQZcM|wk*kAbnXavUXkT(E4D!&$ZA>xIICtm=s zL;s5U3rFGEM=9Rj0xUixe|{Lg8~6{fm#A-D82@Ep9p(>u=9ua)OolzA@j3_i6KQ<@ z16X`mdw^@D_Ds1>a8>_MY+k?5UsF+$Y*Q?jypqYqmcjyai3NW^-I5=IZqAcf<sM6= z%Ts1`S=iU>7Hdr{eqviy?ZQus%ahZy4SIzb7H4@ap;Sfv&c3{+s%n9#$5@olsu@?D zJIQDfh^i6y!CKgN`K%W4E9Ybz6rOJO7iL*3P7ixtz%sp1pISLpeaIlsVYkhpKRceD zS2)3FD%7XjU3Qz_>v3Cb4y(7)PJc!{O>eRI=%>gw3rTQVUCxE<hshTF)Y?*R$MX+0 z?B~;#6sJ=XmTs?h)%Zw4Vfvhs;wkJdW6rqTLNmL~o-@f>J0%bGP?ju05^%Z&(5JAB zW{1b?r<|;ujGXQOVp66Q;MoLZ%%pO`d#9V<s}&b|+?9cfCE#=mpx3jC31*~GR#lTw ztY+tY!<?~Bzt3d%j$PolQ;GCBMW)GSR@N}50MCbv^LXc5z2%IO$|z&4wmAV(+8i@l znd`IJ-R1bTGs~9FX!-UKV(OfHPm#xuUrXEYV`<8lHK)Lr=kmB&z53KSWAW2vn`3g> zY;-7-lfp^4Zg>}OOm&c-F=vX~XJ=V6=9s*AnuT(u&na+McvShnv3ltz*aR}>7`<K( z>$nVvwOh%H&@<)*K8Z07KZ19G5pb`RDyPeZA8^AE<N^V@K1et8#7GTcq`*S-Q}75$ zL3)^^AR|OldW4J&MhPjcl0wR5&vCju5oO^af*c}1qK8PysJ422c3WuE@bHZ{QgpW4 z>W-vlS>2TwfSxK@M0$nYTji7~X=Ug;nUKva9)*gk!>L(Ra|E&4X^S8l+#z}xYC>UE z2u!W23b7Gzh;0nE(CoIb?6%PC`jFg)(9DL=yoNa;nN#P4<^~Scq&lRBGQ8Gur!~Z! zdews?85!kvyv~5dQQ_^`J~(ZJgjA5^4oa&EBN%*MnV$;5GS;V1FTs~Q?hrW{ktD8c zbO!oA3L}MjJCdI6uzURrWi{Bah@@88E3FY!gWD$;)rTNMRR~PO`w$A`+?jau0tR52 zOw(hWR$9drHqYbohj0x&ax!St_ROtx_98Wddb8)U)pmc6-1OOKLx^vrPO?vxbvA^; z)Hz{zd>AeX!=_N2MeU?FSjcsxg%i}+5KcAB3H7%SWT*;(Y4ffA5Zh<C?eqPf&>@rJ zuvc5Etv;L8C99K)NF!r?iq&7`@m0fpLQIumt*Q<qamS9%fMZ8tuu&69&mu<(cOYiH zaAI0GA=TpuF>iW>HEbm7)oCy)9a=`vS>%tT(vXaxrq8$I4J)$!(*1Uy%NjzdwAT1S zMi(2k5mbzBw|9YDGGmV2HP0E+jcHY$P#4ZFv)ks@pl3pQHQiR@b%yjGxm}2G;hZ7F zfG>tqVW@BcnKmp)?O|<o)!<Q8xeckVnpzoW!i(~-n2>3DJ%}X^PmK=}Ll`MwOrk|h zs52X?L#Eqw|9np<CBtELRfLjKtDPb9UsmosPo=fo>yazWv{_uvN=Jxx**1%BZjIGz z4^b`4VZm#8+_pJ#alPJw;S`0DPV|sjO`pn0(L`t4qfzx4BI;SHtUjN`7SbAh8h=ks zG+w&fUc+`7QJ5IlHDR(+L6TdQT2+%G<1Bl%4}r}cLewMNS6MJiRJ*Wfm9raC@U%3V z#P%&UUOO_meIXV$=xr`E2|FIOiFgfAxl`TCh*{i_8fbkU_H+nt7MI;!>37IA>Gfkh z4sVpw2x_isht=Y>xx%c*JjYUHGeq+Ab1Y*qr!B0)ON~NmNjI06EG2e|hgLrt6v8wR z(?qXS3g@KG8E^HKIKXmvTp_+>(AT(~^dP$3Qf0TgvB`7$oPIt83n%C078dF=%uY+G z#b^l=5W+3#g*&a+vePWRlr41jRu<fz%BJV}3VJNcn4u-d7fv3VTWGNqSg={inwnc& zkUO@}XeiZ_g|UUgwi-Q-?N7Od-XJ#5ZmFJUG^9jQ@o+zO;2BvN$(44$rP^lkJ8Il> zlFMp|j7_!0>P0lPVAmH)z`mO{y<CsU<Hnhd(2$hC*+c%cqTERZc__gSAp!>iK99wL zxaG3b7R&6<vJ{tD@Qy8ef!SyFr&U{>-hzo{pQ$Ld$dsEsRlM=cXf*rvDHbzxYrW53 zU{1Gyl;(oxprp}kWfY@ckJn;V*Z3`NPr2P_LZMVq0VA<oRI@&-q(D6{$itcx2L~8s z`n2(Oe?gwf19t+RW@HBj6CorsAjC-vDx8ZH^s21UKG#@0(VPK0`~r>4vY0LNtgad~ zkdzCUkhRcKX7$+>dLbE1?{bzQ98)(T3Oe06IpYfp#^zZJ$thskYOR(sd!-Y{DOSJD zk&}~GT8dwWSoF!Mq;|MR6`6f`dTN!AizzL&q$Z~!D>D$B#_O4n+Tm?3m&Z0QJePa^ zY|}hN)kg7!DT5A4Es@WCP!+Sj*`Hcvtwt9|ctlw52YttEuo#WTRJf}V^-U*FD$K$= z^4KwuL65s~^f}%oSA(8JKhrDajAlBfklH+*MiKkn3b&%oSHm&&s3GW9EIwQQOTVbD zR(1N>&C=^#p80mKr3|ljD>oLaE>`6%o|pym(HUE+S^vv;o!)?<<)!H!3tqq441G$m z$K#(|<HwN8#NeS94~hl`-101q><m=z#32xT8Sx*BIP%4G)>in=q#B>!U~=H?i8h~s z#<3X1lPhqrhB3=}U-hZ<@x{5O30YMO=$(1`d^RDGlju{_w;2|CeB|n8pHn?DJnK;^ z0|R)DeStcBgPue&D-fX4N+uY_K?~DF8+z1PUXK6CMbEl2`ZO3;cB&qtrrBsI=DtO- zS)VRq^CYX=T50$8HnJ@A`z${6(FVPtIJYQ&@+4!?_<|xM5GDj161fqJ)a60YoAssS zPZs>LD_H!jn;<Np_Ru3|y;?R|x$G6vfHHerKBF@{bcTj7I0cJq6T^Mh00yQ7ysnZ@ zZK^vP%YPxy3321VaF=GGU&388F?PsiY;y0R(vSgepaIU-F4K!JT?UGrov@VVb1q~d z&r(^1Rkk|U@&Me+Y9?RpL(E~=1NbBQn{{n)&Ppw&lQVToP>ek!7ITPPP8)ibhXkj` z$0q*lYA@bCOYbKYix(fxqN;F^%fe-}sXJg0ETnB-JN<T&B0!K#VTlP3SJUiIufUfm zjshF<>*+fgtmLftj-kM2#MFef7uF)!reM0jDzvHqi^CMG+A=KRg4kj&BB1EA#?r4+ z)u4`J7<zXqpLU`w-o&Utwg2;K4M)FhM4ff-V6hcO#gJ52q=6~x>{pc%oF98!7g7ut zlV&ze-34yU+*z0|@qHGKPRI)7E>pzXttha%;Et=$P+#Y3b=l|*!SJ_04m0z9jI_Yw z+~#7-My$)oTQf`4a32@~F0A1HEB~>Tz<bgym=Pl!=;yD<;5RTybf>`{gEW?Aqv^~< zSjvWJV5zLFwP1(q^WcrX*o<1{rG&1h?d%uIH20+#i;TtNrxi{%oAnv$F+>S{o9X|J zQ$5Z>Td+w}m(vIi7EguxqTgukj<$_4=+U1Ac?E8p6I)|{uDjA@H#wYYKnSTGZ-7oU z`@PQUTsjV$;+_u;EX8&oPNmSkygaP$u~DEkM~3>=;ev8(w91V}gVk@G%}kSx?W7w= zBDJ;pR2-wKYk$5$ir9Mv*XR+ed92v8s;spZ5l(NO2Sq9o<#7H*ryX{-KF6pLkuHtq zIP#&;p0MJZ08UZ2Idb=<=6n3q+wjARYKvpOQqA|rIGeNCU18-iOUO3vO75t?z{)FZ zEw3)z;f!I{4$Lvq1|{7_-;Oi)3=U+y9&S{CiGrFFv89>D4hCor{NE!aHq$IZhPY3v z>b*e^ViRoWHoey_R|5mB;nkfg?Xl><fjldurhucG)sxKNh)kVityEPD()80*O8ELn z?MSxo#d@HJT`!MLsli_FWg1t5M;eQ3!ToX(8Jbn&wX<Oy8VVQ(GeWIXO#JXvzUhuS zo>L>oGmE}7qJ_EG!m&JbEKuk=>+_6(PFUb(mLggfu&q7~_h{~c=n@*)is+c8SC0r< zJ+S^bJL@#sJ*p=HHXKgOi5TkYqXy+JBb`zF61yy6Hg+(@>gX+W_TqEKe1mFi5uf0M zm>SeEBW=g=Oa@M|7pQpy&Y4D&FAE!I-7{zx#m)zMTz;sQ4U2BhT6HQ3EFjcDMPmw) zD&Pk2bg??62o|LGxI9G&37R7`Gz!z1Q`Bl){%C7hcq<mNWtcF^)U|7g9t)5xHrk_Z zf%TbrE<3&;q?IZ^i%`w~OWHcRF_@`Ie^@hYij15g#gYtp)mi<l>x-wl5HiA|g$M;b zOh}gK2f@ZNY7;so18dnbYG>5Gpmr1U=Z>;GKFWG4|GzG9c>C?f*`|r}l&~|AD0`Z) z1L;z}gv-jq8naA2yyS;3Xe~Rzgh>(QT+zJwf7_A41433>WX_1$MRZ`)1AFQuj4yOl zDjP9Be~B)(XOSx-w!gxoWR>`ujID=+`~J*B<1C!Y;S0Aaeo8BcByosB2Q;Q}#b%G$ zVhlG3ZJpHZ4^Hp&A*XlbilMuDeHK4$w3d~4_#a16H}oMRI=CtSZ^K?7)*Fh{ukS*Z zQt9%RW4x6;AA_#R4t~R-p0Q$@uz2lwK9tsyWj&k-Q2&}ExASK`d;|~kfA+z7hI)2P z-y=mkAIP*&NXg`$`ZL`JQyD*z|Jh+DJ8;m011^ei>IpwK;Oee96~l@-`F~tW<>_f# z9e$RXM&FG<sueRDTLzr{6=K*R>cQZ>E(6iX9!{azymqz^@zUAXU^HJ8OvVu&9&WYN z`0VU#Eeo4Kj~D0LdcC?t%%v~4co+C+xf&h>aKNGFbk7&rs5ujD%<8s?pz#arD>H*f zLF!i}^j&E$Us>Ut1IqF5QsnaqLtqKdBgNTQHVe82%+J0Cl)uW+)7exBzV-FlJk|Ch gkGBd3f$WRHvp6s7Q_Aguso8`Do)};@V`K;ae^(!j)c^nh diff --git a/NPAnalysis/10He_Riken/include/ObjectManager.hh b/NPAnalysis/10He_Riken/include/ObjectManager.hh index 1ab104dbc..336fb8c71 100644 --- a/NPAnalysis/10He_Riken/include/ObjectManager.hh +++ b/NPAnalysis/10He_Riken/include/ObjectManager.hh @@ -109,17 +109,21 @@ namespace ENERGYLOSS // 3He Energy Loss EnergyLoss He3TargetWind = EnergyLoss ( "He3_Mylar.G4table" , + "G4Table", 10000 ); EnergyLoss He3TargetGaz = EnergyLoss ( "He3_D2.G4table" , + "G4Table", 10000 ); EnergyLoss He3StripAl = EnergyLoss ( "3He_Al.txt" , + "LISE", 10000 , 1 , 3 ); EnergyLoss He3StripSi = EnergyLoss ( "3He_Si.txt" , + "LISE", 10000 , 1 , 3 ); @@ -150,16 +154,19 @@ namespace ENERGYLOSS // proton Energy Loss EnergyLoss protonTargetWind = EnergyLoss ( "proton_Mylar.txt" , + "LISE", 1000 , 1 , 1 ); EnergyLoss protonTargetGaz = EnergyLoss ( "proton_D2gaz_1b_26K.txt" , + "LISE", 1000 , 1 , 1 ); EnergyLoss protonStripAl = EnergyLoss ( "proton_Al.txt" , + "LISE", 100 , 1 , 1 ); diff --git a/NPLib/Tools/NPEnergyLoss.cxx b/NPLib/Tools/NPEnergyLoss.cxx index 6285a0247..f41ac75f9 100644 --- a/NPLib/Tools/NPEnergyLoss.cxx +++ b/NPLib/Tools/NPEnergyLoss.cxx @@ -56,7 +56,7 @@ EnergyLoss::~EnergyLoss() {} //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... -EnergyLoss::EnergyLoss(string Path , int NumberOfSlice=100 , int LiseColumn , int NumberOfMass) +EnergyLoss::EnergyLoss(string Path , string Source, int NumberOfSlice=100 , int LiseColumn , int NumberOfMass) { fNumberOfSlice = NumberOfSlice ; @@ -67,105 +67,86 @@ EnergyLoss::EnergyLoss(string Path , int NumberOfSlice=100 , int LiseColumn , i cout << "///////////////////////////////// " << endl ; cout << "Initialising an EnergyLoss object " << endl ; - - //If LiseColumn is set to 0 File type is expected to be from SRIM or Geant4 - if (LiseColumn == 0) - { - // Opening dE/dX file + + ifstream TableFile ; + TableFile.open(Path.c_str()) ; - ifstream TableFile ; - TableFile.open(Path.c_str()) ; + // Opening dE/dX file + if(!TableFile) cout << "ERROR: TABLE FILE NOT FOUND" << endl; + + if (Source == "G4Table") + { + cout << "Reading Energy Loss File: " << Path << endl ; + // Reading Data + double energy, total; + string dummy; + //skipped first line + getline(TableFile,dummy); + while ( TableFile >> energy) + { + fEnergy.push_back ( energy*MeV ) ; + TableFile >> total; + fdEdX_Total.push_back ( total*MeV/micrometer ) ; + } + + // Close File + TableFile.close(); + + } + + else if (Source == "SRIM") + { + // Reading Data + double energy, nuclear, electronic; + string unit, dummy; + + while ( TableFile >> energy >> unit + >> electronic >> nuclear + >> dummy >> dummy >> dummy + >> dummy >> dummy >> dummy ) + { + if ( unit == "keV" ) energy = energy*keV ; + if ( unit == "MeV" ) energy = energy*MeV ; + if ( unit == "GeV" ) energy = energy*GeV ; + fEnergy .push_back ( energy ) ; + fdEdX_Nuclear .push_back ( nuclear ) ; + fdEdX_Electronic .push_back ( electronic ) ; + fdEdX_Total .push_back ( nuclear + electronic ) ; + } + + // Close File + TableFile.close(); + } - if ( !TableFile ) - { - cout << "Failed to open file " << Path << endl; - return; - } - else - { - cout << "Reading Energy Loss File: " << Path << endl ; - // Reading Data - double energy, total; - string dummy; - //skipped first line - getline(TableFile,dummy); - while ( TableFile >> energy) - { - fEnergy.push_back ( energy*MeV ) ; - TableFile >> total; - fdEdX_Total.push_back ( total*MeV/micrometer ) ; + else if(Source == "LISE") + { + cout << "Reading Energy Loss File: " << Path << endl ; + // Reading Data + double energy=0, energyloss=0; + string dummy; + // skipping comment first line + getline(TableFile,dummy); + + while ( TableFile >> energy ) + { + for (int k = 0 ; k < 11 ; k++ ) + { + TableFile >> dummy ; + if (k+1==LiseColumn) energyloss = atof(dummy.c_str()) ; } - - // Close File - TableFile.close(); - } - /* else - { - // Reading Data - double energy, nuclear, electronic; - string unit, dummy; - - while ( TableFile >> energy >> unit - >> electronic >> nuclear - >> dummy >> dummy >> dummy - >> dummy >> dummy >> dummy ) - { - if ( unit == "keV" ) energy = energy*keV ; - if ( unit == "MeV" ) energy = energy*MeV ; - if ( unit == "GeV" ) energy = energy*GeV ; - fEnergy .push_back ( energy ) ; - fdEdX_Nuclear .push_back ( nuclear ) ; - fdEdX_Electronic .push_back ( electronic ) ; - fdEdX_Total .push_back ( nuclear + electronic ) ; - } - - // Close File - TableFile.close(); - }*/ - + fEnergy.push_back (energy*MeV) ; + fdEdX_Total.push_back(energyloss*MeV/micrometer); + } + + // Close File + TableFile.close(); + } - } - - //Else File is expected to be from Lise, and LiseColumn gives which model to take else { - // Opening dE/dX file - - ifstream TableFile ; - TableFile.open(Path.c_str()) ; - - if ( !TableFile ) - { - cout << "Failed to open file " << Path << endl; - return; - } - - else - { - cout << "Reading Energy Loss File: " << Path << endl ; - // Reading Data - double energy=0, energyloss=0; - string dummy; - // skipping comment first line - getline(TableFile,dummy); - - while ( TableFile >> energy ) - { - for (int k = 0 ; k < 11 ; k++ ) - { - TableFile >> dummy ; - if (k+1==LiseColumn) energyloss = atof(dummy.c_str()) ; - } - fEnergy.push_back (energy*MeV) ; - fdEdX_Total.push_back(energyloss*MeV/micrometer); - } - - // Close File - TableFile.close(); - } - - + cout << "ERROR : Wrong Source Type" << endl ; } + fInter = new Interpolator( fEnergy , fdEdX_Total ) ; cout << "///////////////////////////////// " << endl ; } diff --git a/NPLib/Tools/NPEnergyLoss.h b/NPLib/Tools/NPEnergyLoss.h index 005b0b2b8..080c3298a 100644 --- a/NPLib/Tools/NPEnergyLoss.h +++ b/NPLib/Tools/NPEnergyLoss.h @@ -60,10 +60,11 @@ namespace NPL public : // Constructor EnergyLoss(); - EnergyLoss( string Path , // Path of dE/dX table file - int NumberOfSlice , // Low number = Faster, High Number = more accurate / typical: 100 to 1000 - int LiseColumns=0 , // Indicate which model to read in a lise File, set to 0 (Default value) for a SRIM file - int NumberOfMass=1 ); // Number of mass A of the nucleus (used only for Lise file) + EnergyLoss( string Path , // Path of dE/dX table file + string Source , // Type of file : Geant4,Lise,SRIM + int NumberOfSlice , // Low number = Faster, High Number = more accurate / typical: 100 to 1000 + int LiseColumns=0 , // Indicate which model to read in a lise File, set to 0 (Default value) for a SRIM file + int NumberOfMass=1 );// Number of mass A of the nucleus (used only for Lise file) ~EnergyLoss(); private : // dE/dX, slice parameter diff --git a/NPSimulation/src/EventGeneratorBeam.cc b/NPSimulation/src/EventGeneratorBeam.cc index b3da7f51c..a8c8db278 100644 --- a/NPSimulation/src/EventGeneratorBeam.cc +++ b/NPSimulation/src/EventGeneratorBeam.cc @@ -185,6 +185,11 @@ void EventGeneratorBeam::ReadConfiguration(string Path) //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... void EventGeneratorBeam::GenerateEvent(G4Event* anEvent, G4ParticleGun* particleGun) { + //--------------write the DeDx Table ------------------- + if(m_Target!=0) + m_Target->WriteDEDXTable(m_particle ,0, m_BeamEnergy+4*m_BeamEnergySpread); + + m_InitConditions->Clear(); /////////////////////////////////////////////////////////////////////// diff --git a/NPSimulation/src/EventGeneratorTransfert.cc b/NPSimulation/src/EventGeneratorTransfert.cc index db689adb8..894f882a5 100644 --- a/NPSimulation/src/EventGeneratorTransfert.cc +++ b/NPSimulation/src/EventGeneratorTransfert.cc @@ -326,6 +326,27 @@ while(ReadingStatus){ //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... void EventGeneratorTransfert::GenerateEvent(G4Event* anEvent , G4ParticleGun* particleGun) { + // If first time, write the DeDx table + if(anEvent->GetEventID()==0) + { + //-------------- Before living, wrtie the DeDx Table ------------------- + + G4int LightZx = m_Reaction->GetNucleus3()->GetZ() ; + G4int LightAx = m_Reaction->GetNucleus3()->GetA() ; + + G4int BeamZx = m_Reaction->GetNucleus1()->GetZ() ; + G4int BeamAx = m_Reaction->GetNucleus1()->GetA() ; + + if(m_Target!=0) + { + m_Target->WriteDEDXTable(G4ParticleTable::GetParticleTable()->GetIon(LightZx,LightAx, 0.) ,0, m_BeamEnergy+4*m_BeamEnergySpread); + m_Target->WriteDEDXTable(G4ParticleTable::GetParticleTable()->GetIon(BeamZx,BeamAx, 0.) ,0, m_BeamEnergy+4*m_BeamEnergySpread); + } + + } + + + // Clear contents of Precedent event (now stored in ROOTOutput) m_InitConditions->Clear(); diff --git a/NPSimulation/src/EventGeneratorTransfertToResonance.cc b/NPSimulation/src/EventGeneratorTransfertToResonance.cc index 06e8d7de2..270f2eada 100644 --- a/NPSimulation/src/EventGeneratorTransfertToResonance.cc +++ b/NPSimulation/src/EventGeneratorTransfertToResonance.cc @@ -69,7 +69,7 @@ EventGeneratorTransfertToResonance::EventGeneratorTransfertToResonance() //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... EventGeneratorTransfertToResonance::~EventGeneratorTransfertToResonance() { - //------------- Default Destructor ------------ + //------------- Default Destructor ------------ delete m_InitConditions; delete m_Reaction ; } @@ -77,13 +77,7 @@ EventGeneratorTransfertToResonance::~EventGeneratorTransfertToResonance() void EventGeneratorTransfertToResonance::SetTarget(Target* Target) { if(Target!=0) - { m_Target = Target; - G4int LightZ = m_Reaction->GetNucleus3()->GetZ() ; - G4int LightA = m_Reaction->GetNucleus3()->GetA() ; - m_Target->WriteDEDXTable(G4ParticleTable::GetParticleTable()->GetIon(LightZ,LightA, 0.) ,0, m_BeamEnergy); - } - } //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... EventGeneratorTransfertToResonance::EventGeneratorTransfertToResonance( string name1 , @@ -377,9 +371,24 @@ while(ReadingStatus){ void EventGeneratorTransfertToResonance::GenerateEvent(G4Event* anEvent , G4ParticleGun* particleGun) { - G4int LightZx = m_Reaction->GetNucleus3()->GetZ() ; - G4int LightAx = m_Reaction->GetNucleus3()->GetA() ; - m_Target->WriteDEDXTable(G4ParticleTable::GetParticleTable()->GetIon(LightZx,LightAx, 0.) ,0, 550); + // If first time, write the DeDx table + if(anEvent->GetEventID()==0) + { + //-------------- Before living, wrtie the DeDx Table ------------------- + + G4int LightZx = m_Reaction->GetNucleus3()->GetZ() ; + G4int LightAx = m_Reaction->GetNucleus3()->GetA() ; + + G4int BeamZx = m_Reaction->GetNucleus1()->GetZ() ; + G4int BeamAx = m_Reaction->GetNucleus1()->GetA() ; + + if(m_Target!=0) + { + m_Target->WriteDEDXTable(G4ParticleTable::GetParticleTable()->GetIon(LightZx,LightAx, 0.) ,0, m_BeamEnergy+4*m_BeamEnergySpread); + m_Target->WriteDEDXTable(G4ParticleTable::GetParticleTable()->GetIon(BeamZx,BeamAx, 0.) ,0, m_BeamEnergy+4*m_BeamEnergySpread); + } + + } // Clear contents of Precedent event (now stored in ROOTOutput) m_InitConditions->Clear(); diff --git a/NPSimulation/src/Target.cc b/NPSimulation/src/Target.cc index 122bbe8bf..51eba3244 100644 --- a/NPSimulation/src/Target.cc +++ b/NPSimulation/src/Target.cc @@ -653,8 +653,7 @@ void Target::WriteDEDXTable(G4ParticleDefinition* Particle ,G4double Emin,G4doub // G4double dedx = emCalculator.ComputeTotalDEDX(E, Particle, m_WindowsMaterial); G4double dedx = emCalculator.ComputeDEDX( E, Particle , "ionIoni", m_WindowsMaterial); -cout << dedx<<endl ; - File << E/MeV << "\t" << dedx/(MeV/micrometer) << endl ; + File << E/MeV << "\t" << dedx/(MeV/micrometer) << endl ; } } File.close(); -- GitLab