From e44db9ca6ac5fe7b9404cb3189e22d4583c6a0b2 Mon Sep 17 00:00:00 2001 From: Unknown <unknown> Date: Wed, 2 Oct 2013 10:00:25 +0000 Subject: [PATCH] --- NPAnalysis/10He_Riken/Analysis | Bin 58961 -> 62985 bytes .../10He_Riken/include/ObjectManager.hh | 19 ++- NPAnalysis/10He_Riken/src/Analysis.cc | 127 ++++++++++++++---- NPLib/InitialConditions/TInitialConditions.h | 2 +- NPLib/MUST2/Must2Array.cxx | 6 +- NPLib/Tools/NPEnergyLoss.cxx | 2 +- NPLib/Tools/NPEnergyLoss.h | 2 +- NPSimulation/src/EventGeneratorTransfert.cc | 7 +- NPSimulation/src/ThinSi.cc | 6 +- 9 files changed, 133 insertions(+), 38 deletions(-) diff --git a/NPAnalysis/10He_Riken/Analysis b/NPAnalysis/10He_Riken/Analysis index f5b004c15c19c002bffe3614626ee9f880d82d44..73acf47f2d6e340653baac089ce92f0d1716e695 100755 GIT binary patch literal 62985 zcmd443w%>W+CP2@2?Vx4MRXBajan4El%_X|A{N>Jl>*(iXhGc&+Rz5lq&6uQ7GJAS zrIaeGsK}~|cf9P1cf73?5K-CPil``{tfCe{ML|VG`F+20W|GOJyt}^m|NH~fIp2An znR(`!XP&v7las5QMPrhZk_`3fVstf#maMbbvJhAQ*`#C{cEe^2HcmDM8r_gYT2>=? z5JU!KvLelb(3LdctgAssI2M?JG?vXEKFo`7=w`!UV16UTq>HpcR@)&2I|@tzk@b?^ z1yTQ75Kgz+7&f3l2GUsv!#Ln^sGFexm_gBsKM7SXxY{rlBw1|?^H4s6qW>Oru*ZE= z`9po)=|g>$=lHz-`i8V%AT7hp3(-W|#!j3p@={Lb>5IU!`yo&_r1#AJ)N}rbk?(z9 zb^Z2tpZ@!gD|>!cG9R=t2zG>?2-8>fJ>^e}8X;-pb%pL@SjLbxd!u12Gx{y49A#r+ zSC5+gV9#D>+6u0<oMl1nBYJK)zl(k0<bKB$+uKS|f6v9%ltwV7<fZpKCS^D}(bC^& zwB;N98c*qxwW4uhV{unYS#r{?mW?+xE?QyOm-gr~y~nYQmfugKEaMOg5sDBdAn-F0 zVG=?y0;({|5cn}$ZC+o5sNB3XrXapdUb|&H9dU)cR#BC@h>2$c&ys0A8T%1mj!-8t zqaZ>(LIVLle>5ZGe`I`>j8)W#3n<sPMkY}I3lOeDScuS!(1Os4umpjhKbetnlZ+WH zMSwXPw<0V@;E=f;;SPkqAgn-8pS#tCjPF5wAHo9&D-j+<cnCp#9+nq>L%a&%F^L~X z`~<?22u~qAjqnV@YJ}$yo=12AfuFx4yo|60VJ*TsgjW#wS6!Ulxb5MD=6pN;AQ zabDApH{J6^_O+RB7XH|H)6MJpElbXK|FPw~{F-&C&m7Zu&y<cY_P7t<9~}9@`Tg(v zQ^}ulubsPX>xlPuKl9eCQ=e`9y7%?Rrp|B9SoHgMdw0EhQ}<i^*Ni&+@JDa9kJ|c= zkMBI`je9;>=9{r~-qc>3=G3Ij82{GdgOk>Hx7T%b9{PO91HrrhdTskpQ(jHo5WM4@ zZw8$0@%FF#`H7p4eK%><ilp05J}Jxl#LTBMrf>b?k!3TUzQbRz=xWH3eadae4gP2U z>;LibG4G}YAOC00qaSCy{M6Dp|LlJLKX+`pr(^4l+2;&)zFqW(o!k6B{Wz%N*27a@ zIqv(5dXMT}cv0zB_dNUcUsC5ke|gG`e^%Y$F-qrdnta!FZ~ZY~$-lDh$v3uKUOw=R z+6%wm`PjAXXQe$5{BEmf>*+NI9^RcdzWtf|zFpw&@%B@F3J&e<6)JwY_(t0g-JeSv za_hhPe>n0l{>7)x{C)2sXYX!$bKkxNrx@S;^~&4+`S881n=`(;ea6$n?jN1oKYe?q zw|L~D9(P`M-h#1vUV8r9y2otyRTcbo{Ksw2e|+OS>lx*JF4_9xXS?6}s<%C5>bvj! zFlXDxdrl9IcD}vGoxH!gSHF+$d|=xT?|7Sw_x3vX*55eZtG=Qp<ALC?9iFY@AAIVW zi>haA9siH<?KfZhyMjAx{~9p!s^#-8Z$Gp*>)8zf`%n9O1kR|t^L+Q)H_s`Y_4FCj zF1+mUwNHLM|GVG!e)#hrU7ws{&#YPft>d<g{M?y6^56MK$Hkw0pH;C}GZS$c>=8-n zi%ZvUqVdn&qj3?Pr2LTo5}QJ%?;IJ!e@ap`{aQ;j{<cdrUK2yVJ|_R1nDQ6L)PH$Q z{-KfL4Lh#4$JCb}!*A;v$-wrVgSkw7i0_Wk*OVB2JRief5W_z+rvBf@@TbJICp)J6 z-^SqQWB8k6=*PzNXR>MxVmq&6=;y`cFN@KCc?|#YG5WiYlcD@jUr)!h|CN~j`%8>I z#>eoV5W_z{rv95^>Z31LAJ%`h0*HwZ#q{q*G4+p&sqd{A{vTuHyC|mqb7J@_V(6im z{(2&&{`46BZ87@4Gp2t|h_RpfG4l6}!MQQzeHkOqxEQ=CCjWsLx+8|Z4oZHb!(wAM zooe(l2L29fj?c|?C2^7+>A((2KPSa7_V2XV=1F`8%Cir#+7<&JV+=4dPBe^1KC#$X zE%Q%-0l)sC#m4m+@yIg_!@oOHe!mk8qvL?ZMzcJ|=wp0}{990-i3b>~`(WONe5cC% zYoU*)w_9wq4VE_#{O$WA<t<M&jG3UHEa_WO@wV?H^3MkUy)R*Hj89(-rb~}OIjEoA z!t&ld4gAQTih78v?U?g-S!~EA@Lv&ownf_KI3DYy?<}@nlK#v9!`K1&*=^K^4;k}N z{#4+;C~rCR@o!Dv#~|YiQXiyGz+(PM=uh>>gi{RTbJVBWlMQ)y?6cS!WqBhJui0U- zq1%Lf&B(a%>j*u7^4lN}8%%lbg+A?PM)Z5iAgpVUpJvSV7WX!cGtpmXOZ*q$|Iy?b zdlCr9Z%4%Pw)BU+p#7@;DaT<wtJQxd>K_0)suAU1gNs^}7uFZ-iTfdtg?|2s@pUcQ ztJ=TUZWz6xAGVY7-2y|{rrF0msQ72}54ugzkA-10qkLsAiy+TCpey<IqQ5q2<xd@m z=4kfP3S5Nxl>H8az7~FJvHeTh!%(#MLXF>vf;%uCl)RU~-|T-`wZ|~7K>vLJeX90v z#>G;|bGnpw6zpXs^fOuFRMfxm8;gx=VA{h&u!jw(50Pk}9dYU!vwXI3#c#1^fPGsf zeJJd!rzT$y(95@2Y!`yg^74VZ{o7(gvqX6<Fsyepdt48F?EzjX>7zmK4f|E~H=v;d zp?{FEUHtm+_pt9G$lD0a@`r;lRkN?1Xx|3S{@!4FHGBUScpmhl?BM}4>SySq80E3Q zd&B+~Lmt(CPoe%j@JFhDSJ^>Fd8+*np`d#;eeXdzPkd{!mCN#OLR^FXSN!X{7{;Bz znUek}##gzfkEy6=1N5uf{}MW`5$($bzZg&G-%3qB2k6H`zpr3kq5N0EUao<BYJ6P- ze^-X`LbASA)PJL9kE7tP8qt1L-w^cIBD4o;66JYNZlz{l_roqP#Q0VI=|q%21Gq%; z4~6_2_eT2XeB?J!-VjM2L4BfsmH%A;oP1WKe?I`>UTJUixrbn%d0P8z=%0<q4-)#n zx3L`Z=Rtl%LSDECW2feyFGqSE+N;`kd5k?xg@3tP)Bn%t&vcZ}{-S*?gb>LnAJvNT zPlUhU3H^QsIZ1y4{-nROPrGq~@gC@9(%#q~U!c5p%^pv}Kv*j2%wG@vEW-Fv{%<fE z(yr0pgnS3I_W0pX+Td?CNO?YmJzS*OpAYRTcm?`EgIRue<gZ3~%KrWidH2B|fK9aL zUsx|JU90qoefMcN-#~vlC4Cy~sZG=WC>Ykgnms%PzDqQFzYy)e9^>UQnSUSbX{lyk zAHZI>Y5vOseSDzVV}I1Y1NDQ=ln?DlM){K@eoWJM3iSOV46z#ZDgH+2dzv<1893j` z57c>Zv}3p~ob4*{RC$9TPhGyx9SnMcMs29N&f~6hRk-~X9-rZIRn-Liu3*Sr7jn4_ z*VNLGFH~4j8ge))s@-+2P@UTw3Kmv4OWpY6oL(}%bb4W#t1R7B=5RSnbD3r0u(C<h zXL>3^{wtg_O<WqP^ZKju9nRwM6-cS{70)h>Vy*N!i)Vx4zhb6{OGCB3&a2hygakRR z;_;=KrnLW6HYu+ov+N>IUC<lw=U00wW<g9Z#m+4&4F&4lRUU`KY|kj4FHjNPqymR? zwotLxAEJt(?UKSt@@k?-rHzBmMGs7!m{m3|V~nYZQisc1S&9FwZG34cXEu};s4J|5 zax0zAte|(U2gQRJNAotAa-cLvRB4o8Y@K_C6oHqSl?v3l)n!?oN0<cVjfh?x2!?_c zb>3Pf_XIaOj>m#5wzzsu&|49d?Z}1BC~YR%f<6$v8BP;wVWYw+j&idlTFR8zl*|e4 zDl;vd4TV=q{e;sp3j(3|DsxLcp|WaE$nB&$YrG|qnxJR1Kc&PuU@flpf)TAk6RX1m zgtlAa_E!dKP_eMx_-=7zxhA;k-2N($tAggD48qK(S{~P98SW4`QDK=Fa8jq0K3w3; zRSGgej#R*l3Opg9V98lHsU#3kM&w8*m{cFCt;fi$@`MWg7_u0_&S5ZGj<wQs(Po`- zSb--Pste39RXm&}6#C-|C=F7W5|6tggaIKtdwk}^;-XxYn2}lHsR*FPy84PxiO1&- zd1u4Gy%p#yDmsc8iWDXsm6+i`S(W(#{|s+ceI4=#{A0X6kJGH&>G#xC%_#~5m8~(a zBR!IgGE0MAuNbAVG&Le8)CWTu1@4er>QCTk6LzFem|R+xF{-Z4Jx7VDr#W)QdO~Kc z1)kbK&>J$FqVr^Afd}Kp%Qj4!QRs(nHB|>*C~}nFU+;6*l{z$0#=3*G?z&1XrO=Q3 zZXXm<DGE{zcjT7Ek%XiuS*t2P5U8v4`rRRq$OvM7(Ku)EIT>l`X<2C*QUk$IWyM)% z0i@-imr6q!=`I(%f8YvNpg!aZ%rK4G^w99b?8(w#j=Q$j<FAyC8wsc$!#yZUaOAp5 zrejX<cwCj9iaJk?#~&)jToFn32P!?T>F%J1r6Jp}ia-t8?s3kXX;u`)9VrCEu9H$q zLs^g~;pE^fL9R^H=B;qSlVM(~G3O>HrU)rOb_FK4I`*ThqQUK|^9T)y)?`&c4`dT_ zTX-O|#kAlkiM)Q*+|DJPSQ_vJ#Z&@a4-ad~>koP`xx&k!_UUdPhwb!Ew1!b!V!GDS zpu1Elo^!=;vzguiCSynzGdnvoJ%JhFMo!0^@5c-us;Cz2@C6GKWJ#dBiZaaU*_d59 zX`6kIYji&+e^VvHlnqVCBEroHQjGS}P%cChJ)cDjR>KSRW=Y9s$R0{0{(s4yUEvD^ z)l?ajJ%^P`Lfx1vmjta9%NS?rFc%s)y|6UPg{E^=gNAa!A}t2S7gi<a<>+RHjpjdV zrefCzh8$cd`ajZCRsh43qZ9m@u0rU15)x`K>_nISXN3&&23?{Yv+roO#bZiJ1Ep-d ztE!>FRg3je!0+~XLvviS(?!v=fg=wMh}4NEi1?MTFByhe3^yS?iOe~!%P`yMqZ?fo z-DB{8HKwC;R|Ld_R#kz4&m@Q8n(V*Ai(V}81nX-&bkcGWOmaGLg>;gIli)1iSY7o& z4`w#uG2mnyydeY9di=Brvdi^yjoa%tT%JlyIv`9B2F-LN`Nh2G$X0Xu1h*gS$~wnz z>D{r=3s1+w+0e3-$0#ff!s*qzy>*2!lwwY=qlSl<2Tmu{obD?1_^?`V1Ve?eS|nxp zJV6ZeAQ$RPaXK7WwqUk*(Iq?SZ!%Q{B1v$uV-w9e4$g9g`D&h<VlI&B`k;C_19&;} zart^gVV$1K>BdZdeAFA6>ne54cKhnlz$5tOxPB<eP)g6Nao0j?aY~L+XrwzzGoV{% zws5g#A-NZgDk&T_y2zPP?x63O7+0NRSiaBW#w5h)(p*&kAM4NL6fFjCr1q(--89dN zV7+UG*Duyfn2RDEFII7`|FGw3YyXc@Ivl>h6`nfR^guntE%~*Kz>&cPy0@x2RD)H3 zDRb017`Dv)kEFy7%KWe~9PJU91G(hvU+?!`4lBh-$QCP;?6D=IipS;F%;8SaQ6TK* zsE>{ujE-5JIp)}i7y%tZc=R9jO2u|b@ApnzPpK15oEtS4{ipJr-1;8nZiw*_o#-Kc z&}KyjEm!WQ^$f@K=b(BR1Qd#*57K++4wIoOETi1RP-A-|s8(ewK}GqmAw)Rh`# zTHy{86H9pc|4U|z$U-n&=+T)DhhZ#GCt@ZHmmMu|)S%2T2PLMJF!#}sB<h&w2-Cb; zEFN_MRKam3uLwuIM%AmoXaV6DZu3`KKxkvuTZ6O2(&4V6>FHQKWeVRdcIn|$xj2Po zWR0%(`YKDA)nm@=Zl9Ru)y(<dua|0SRNkX(cfv7a7mEe0xh4!xM!%?AYUNCO)Qts; z+wip6>9UTqhnh<xCtJCg>GP|xDi<nr<cRaZe7COxM>0IG2_Gq%du+5jylk%UiQP7q zcji9omoA!M>JTjPIvLAy7Jd}xknlnli~jfp{*mV&l{eA)Kf@u`|AqbvZzUG3qx@C0 zZ?Cxxz{%W1XUW(pMUzTP9qHpdnWZ>^tR3YOjT?UyTbEnoZgA;`;{W|<GS{#(YF*V= znEiUR2UW}x*F1SX`ePy+W1`YL&NT;bC#x#6k-u<PqulgS7e5!$W3z}On<I!F-pOQ} z$LwrMoX1rnhM6;_n47N7#@rFd05-;A6UEF?Tk0d`GF=7N@sZsYHe2RVV(1s@J~H<n zVM$;*mn`9;A{v<3$+j~iJ^Hl4ynSG{I8Lz+$EeCmbzFhd_P7<R=#B)((;1oZuEM-8 zBetJV?NJ_b8C>30c|xv<JZlz5kX(SvcImz;6TYd;T?gKaasFfOvFIp|pkZgEAC9uv z+zf_y?MGcH;<SfK5vM)EsE_1`kEj8;9}XK>7VmWZIy+SE&BV^V(~17q@!FA1VAQHl zb^b3^-7x8*;pTp&^{S>Y_^hyrM(<;e{+M%Q(zas8ni@!*tae<9PcNO*d5VlRQVv#d z)7ibhRNHC;eWd)_`rbVDIO5=U=E^4GepT^!F^Uqcj$%$n%gtqxQb)owTQ-?DcjR*m zIOX`I3rw{?KqvpgsU{a^LN4By>uf5>6}y^5`<Y*`%v8M;b|)DVGkzUsqV;vus8(x+ zqwaJgr?9`cFLKQDL7wR+I$0m);z$`L_awitzwzvdI(v<q_K$F5{{JvChligQ5S~yB z_RiK1nb~7`tHD|04TU_FVy7XN4@K_j7YFKmm1FAMH6A&_{<q0H-;rn=l*JQ!$l%7! zVRMgi)GZ-qNOjViNAx|g(jzR6kGjPA9gL%1IBYq*g*Q*Wd5h`@9uGTXbdY&89<{y9 z2ruEry665nZiAZ>?RL&loJ}0*m7ef~QjCYsW*TFQ3P<O=@SKJ*Hb38$$yla1(lU+w z^72um3thYv$>rr>$-<L>{AQ~KfB7%+iALli)#NaF`Ik>Zb`z<5G}TP!cRDTjI%T9Z zGpGL7EJ>BA%1biLe_dp(@X?R;8oK`<YvXrYBmbgNBt7w!DF?rw8voBM=~O&iipMk! zJ3@-rizhO>8-Ky~F!*7f@9+$1G`(9C-Kym;K79_5)p#(Rzl)I}>FFphc(NByz}t++ zKxdg|eO<Nkx)?S*<4PJoZ#_B-24FE>l>DI~H{x01F2)8~{*E)Owr+SPmESS3;vx3E zz}@l8^M1xKhXagZIENVH8DL{DVwT7l_1GA9$1|&pQF01nNZ5}t+TEWq<R8cw&(7N! zLqn;I@sROg#;E^n#t=H4F`g06VvJ{za~b2A*b$8R>@#CLLtnrc&jycUjAw3(7~|R9 zV#asjnPkRzR=b=ro{_$kF`k{B#&`{$@n(!?t*aU1;ddWnJUd*=7|-N~82=I9$7PIX zj^{Bx)i4?vUybkDGRCvJ3mM}X;Khsw;ky=$@o?@^#?bLH#)I*lTE=)Nd<El`c&3=~ znRr&8@mcsz4`VzOy^8TU_>KzWG<<i2@lbq+g)t1Ejj;pYgJqn7XXhDb;u&$qa24&0 zvkhZ2;~ad~hw(5xgUt9PJZsDt57+NxJRIMxV0<o~QD=M}`k(O#^grYC(f^D`qW>9R zfc}S_bk9ToGaiNhXFMAH&p03b&$s~n&)A9YRX2Z9I)9I)138ULu3no}jjtu&bm;Km z!%dq)){dW%+A$gC)X@%^$}f3s?cwYWSlZ!1d+?We4iBP5B9~$8+JS6`2hlRCN!!&f zFfFv2=Dw><U|MQ5EWy~dN?-$$SKEnK2y7u9OuR&38!<F%>}nJ^nHbU-yJ`haA<iS7 zCU8IEam2*}_a`nU&J%bbaXE3iz;@zk#CCyGiP@lCDFP2Bt|c}EKAX6Kc<;|h7?e)j zNW4SfEaHX4?E>c#FClIdcm(k>;#C6Y5w9R#A#eflO5!B~k0V}1+$eAn@oM5)fs2XT zh^Gl$M!bQzSm1KvcH%sNFD2eWoG$P*;vK|xfh&nSh*JcvCf-YI2<#(1K)m-Sw!fAb zUKH&YI7H0vf9+})xPdr@xJ}@B#Qlj^3EW6*Cte}&0^-5MO9Wm>oKD;*@M2<qEpS(@ zz)Ohph^Gm>l$ft**i|g>GU8(5Jb`Z~E+<YGcm?q^V!Oci5?2$a2)vTGme>&Z5#k2o zy&NZlRuMN6?-2M&;)TTR0<R`sLfj_s3&hKaR|(ukyn=Xz!0U)t5-$;W1Mw>2Mu9gH zuO_Y)xShC-c$&bQi8l}z3%rH6oj6b6ZNyuM(*@o^yo1;-@J`|m;uL{9i1!j30`DO{ zK)m-ywtp`%S1P-92)v)zM%*s&0pb+mHh~Wj_a|N@Fg;eaop^=7^jy`0iI)gW4_2K{ z+$b<TS#>UPt-$nX)p^9z1g2-J9!FdZ3_A^f9AEsQ`QP(*?kz4WspfZTs`<sl-(OVC z@27O+gL(gwsE_oGjpsjtB8~b{Eyu0G#p)Ei7UY~{{CW96m)7#b!`~0uTc2hn^VcLb z3zozEZYJI2AO3#br|YdOB5D3Q%O#hsZT_zH1u?s${P}zO-_;MSO>LpR9Mc9adc0~$ zW!mBP*VZm+ElRaO<inRA0{6uimHz(fJq>7E^BeQevxL%KEl9N-KqHWY3xtl|=pm-@ z8#K+9)Q*uDvCp&EUC&^yZf(j2X8L+qa!bBtew(GmXuVmamZjR7k+ZdlUV%K<Adcj? zSD+k8yH@45F~3L+=RaHWj7Gfsq~^Dp*GEeE18IAk-)H%2pzzivRbl}stxdyZ-qTh2 z$t+)_hRYu;c|HQq?i0=ORVhCsjpgrp2uqYmai1Z%wMk?YP0W==1b&DVp{jmf6>)(m zLa{Co(kgV3R<}wUtulU)g4bT4o~wSa*jk%L$h!8S0H!V$scbHluw4`@c!Xps-HMKE zZBpp7gOIbo*;JMG7Sh<%-6@eeULv1ViK_Qu$+v4P#j>R4!WufhL&KZbHgBN_J$ZfE z0@?9GiXti9qys4LxaN17Uu_<k+HpA<4>W(X>nE%STbl$E*jt<nSWE?YWg<{E<X4iY z(oR8|=$sdkACmN^BtA`&5Y}tWmo02<5+)=h;S|QwU-@2_<TKPIBnb*h1e1_tA$h$h zrnO1sd>9&5lBl#>ktQU`ME>_W-rKSL4$gh=eYyRP{xVLHu}#LD`$#{4UPgS`687aI zqA$xpV!g*e(@cF%q+Tjhe?$XYn*@(&kxD;}{eoOiVe#COa`l=K@Te_4kw3+_QrJ)O zD=g<`3e*haXUVT41$hp?+1jLv*~l)ErlHbSQ$5l&lyLWwPnw2GTO#>(osETlq&<~V z@==gH(jK2oJxjDl@QC)P^ip;t+Vcq<P2%=!744y7ReN?s+M|j&#EzEjQE44WQ|(d0 zZ6=>=k4kHke7i2i1}4&;kd*v#G&j<oc{249q_#E*9?>3^KA*}#dyeb0Js*ko&_Gms zK8v(R6?3-M9+lQ#YmXAnqP0h*?fuSd_pZg*m_^#NKuUfWvP9amSf;ihwY5p`i1w)T zCs--k<AhsF+&{ZSduZsYJzqrHql&o{GO7MiY2&1vbVW+IEb__zQE7J7o)@s;jI?K| zl>AMQBkj3eraq6<)+WIt+N07xW?^VgJ=|L2_IxGU!vU+>voF#fRm|dVBkfUX4M<b= zri7~`pKOmxE0%n__Fy9&Y0te<@}#KtJR(zfBU5XW;1TUn>Bq|U+y|HR3)-WKd0~H~ zJt}P_(o}nta7)Q2+oRGNCEu=rQSEtBO74K<k^Xr>rk*I;BX~r6RC+#@gZ_CNz9n&c zz7_pL$E*71=ScskV)pEdv`3|FL7Hlh5^f#&WP4QFD#^F22nQamO>anvUqdq`Dlq+O z_k*0N*GYP_XpurKEClUHLOYtBw@27Yz@|omVjjsx)ANpqR5uf8Eza9p4C`hr+*(qg zLJVhNFgEynVK{vu0C-T(;Z3R)O58JLnY*sFA^ZPQ(8*ds{~NDE<HdA*JE!C4A(~nk zT>dp0v!%IdKSq#Q)`qACv^I&plRAuQLU_&a#WyH}Ylhcw29>xF7LS#unBa8rIm@o# zpOyTxihow~PaFSiXl;5+cGp2vA(1R)nCQ-plKvbzu*JEs^=6S!mTF;do{t%#1^wCH z+NAidp#{mAN2U3t1??UXsc$*bUUhC}7w%e*qY-mirlgrCbL`rJ!?%_Ro1gn@%hu}V z4c~wM^1<Zh4KE)!x%sUjSa^BcQm_SRetpQA=Juwx`r-4h+-!t0Tb!G@23M>3;;+mB zcRt!fq0{%8Sgr1-q2;Dd*KzRx3)<Yg=H-JYD+^;;99rO^WxWY4u$Wg)VONOD({Ym4 z;%sMXDJt83*_Ia3-5lm@APISJ*4yt`Mk9!$SwKT~T$<W%dlNZum5hL0#Ad1eM_@#) zO$UHlCbvTvq5X?d7Y?0To7$15^=82fIa_ZgpJF=oOUWb>AxIK(sSZ+UUw}yrxWSPs z-$UB&OtV{9foQYhlij7n`>RafH3U5w-@ZyX!tGYOVYPkm8f-A0fv-!7Fsy#BgSrqJ zNfv+6#z&+Gccb^VJ6D<NQHg*3f{Cl66IUZq<`XV;Ahy|TL^~8<>df4!@jQ#^t3|rt z68tKCGx~7*Wjk7%RKiOrSZPqDJ&ZKvOjSj<l21BQl{R1U?P`G&*zVl2-8vX*dGEFN zcYV0y&OW!cEock<al7yc?|m7S>M>K1<%5@DU@^$K18q-cc`zHfl~}n4N)an@*;=$Q z2hMN%9bBA5($Yj4C-~hn5>{O-D?_b`>d8w`&p+X(V(VEi(pXR8cF$qyrj`7#{f-sj zi#L?1kagQ%NyFYDp)`;0j&F5d0%>x?((F#%ZcRrgn7UHe)RDSYGT*^fca)S|wsMr4 zTbN2|NfXh@Of|PKKAD<})Yc}!BkW71PoqhoZ@!_c!xn~zU@d)-a(Gi~lVUmxN>(<f z(te9HWpk?Pt{R_8`wA)%d;=nzw8((^`twN1YCLU6nkre1{nyARty!f#CHZzaAoG0Z z<|O*gwax2a{`q9nqw>J5B}ID9P0s=u!{3?LYxp}?S%$Bzvp!4jxZt%uO()1r-z4~G z<+m52@RrG&l_P83Wh(U&>BhYvQ_0GW--eb)x*{Z%z7}ndbj3WG`d6f4w?ZVkLZ!cr z^o|vRsFEIMv*moI((aVa=1N^vv`FJqX|omIapLqyo}HkpFVlKSoTnhIQ0vXJVOBJW zY>iK)oucujOTJ#@lZQkq?FTf2eaLd|fbN>Nw0ic=UuR6Y<PWaPTHS}%igUBv7H2Yb z(Bc$xW18@FTO>Lj^276y1%=pvwoTZv?aD3Naf~q&MQ*=r?{??j*H{wH)aLK;&A-x? z)Te22SG;`Z%4D7z_G!8bfE{2QZ{=yl0lajk!~YpN#gh61uh*c|)qK?zm~4A_2nS|T zSmoQS#5(RLiw%bhzccVcAyWDJV`S(1!#YZtoo_PIvy;N?x6ouPsYP(L9fKhUmt@<m zxg6`4ruN@%%`>mWTC&5G)7DB0L?sVMGb=#_&0FrVqev9xV^MsWq*>HBvw*Y3gq4bH z_C^>qB$ic(bWuFhDdOK`@vq3@@A*`impmJ?vY)Wdqwu?+Bvj0|V6{5mWy-~%@=~tP zHlm1DbJg)X2ni#VrM%P;Ddirj8mA#>9qUi#2&5crE;FO0XjUFoX@e-TycbauPONz7 zXjbqD1>=f*l9XjD5?1q#PAp7y8;hZ!CqaN#=ZCz!M_x*q-xo5A12Zq{F#U<yyoIt3 zX4?Z<5{X?2@mrl+q}YuV`}d!)GErnV7AeIR3Exr^QfwhkY9g^&3@fdIu3DX+iTdTG z6#FY7_6|}WN0CZwQRJJl4idAEDE1C1b{oZBzLS-SA}6s(DYi%$sl+~0ic^DxIBl~S zR=SENBo!%yzd)1d9%yI&Iv>ZpIpF(4triYVbUTfq<AOtn4~rp;VJ>(SqjUmyK@~Vb z#~4S#_B`-}uls4&$;b}s2d@dPeggIT;I$847I>h*1Bt_VTJkK~wavW7Nj_LxCfnw( z5q_z;Uus7a#N2L}GHCSvR43`4#PoDUcP6HfQ1n5G=?#i*jiZwy)$<3=a8W)K-T}gz z`f4#+CMUoAQP-sQ7_DA{PBU9X$C(F90HfGtOtz^R?;I4(^}TZ^^q~@O3@1_uX+ifO zaY0-C2nojkHm_}U?$+{g<fE|WOSzE|ss)|qbB4$p>dt#%#%`;)eS%f!<*z#TBLSmr zKPrmk*~2_z!Q1NG+&Xy=#6oZFh>?7PSs7OY&K<2zk>@AORmx07!X8dua)LVaQFt^o zqf*>aY4seyWOd0OV)Ql~C7<U^?sE%Q%E_?d*%+D~t;lX_6Uyw@av8;KnVdqYL?>dt z8fc2tyrt!`L$Bd=!|0y~uVN~UVUm{K0)bLA`%ED}wFp(%kwx@&d}*h3QX1AnxFIuq zO`oN0TI<aI7aO~ic3|b>OyR=G)?91}w=5McYjuj+*)6F3Qz%QcwV&A$`}HQjM@^#3 zRg+a~@+{43P%<{|Zw`5xyXmfG+h~dQmU_|GDBFj%K=WGJ%u)bTqbV>U+K)1c=B(6? z381p3z2ROnHH=y=Gq22XYwAsSt#i@Fmpe`lw@vmMdR}!#i1lxozpmdUvJ5s<s6_S! z<t=7s*HfcZ$FrE#hOdWS7qs<R@Perurq`e-+!7qy;_S!L<s3GN6|9$8`vE{BU_Opl z7>prlJn%4GVXCmV*fBCp=G}s7Invn)?44vzil|*V2Hy7h4m?oNt{>qZFkP9<Hm%d| zK(%0K!Bp0D*W2KN1yG!8k&#kL$v<M+&8#q5#A(<QvDFzcI_7k?N7T}nI}4?lP@a_f z5d2z8x&&VlgqskG?qP^OOo~532qtCun6mT`^2v3=lcE?V|BH;#BEE>323o9XG1nCL zfw!qmbXMpbq#ap*Poq5Q?|HUwGCDk^WkSE^BIp|)2s#g&BZOTfd2ax(@Sf}#UDpn? zH7(Aa(6v`+Og8CD%ASM#RJqW9xA=za$Bi%Qnqyz+dSgV_8<`f?wX%?>Me31TxEVZb zA%|_Nb0Y<54I3plg_;$7MJk&YqvEMTvQU+jeKhhb6`x6!9nnNfO#yaa(lU8ti*sY= z`aM%JX0$r@#v8J1M~icBbJyK{C<%rfhrm5(7<VDkvBNQhJ(8UL;V<9-Om8Ip>Fy%g zmOIgBXl=`c{^)hntz%?*v2mQgHwk8SJo%hhc!*xRkyhzU+1*>>Ruj*^GR5d~Gu%V{ zgK|VsFq}TccpB6L*DL*ukP6qBGUo3YOahcYmIVr5@;xDtT%lk<30Cm6<jED1@~o`B zevBM)*dQG~&)F}OLrGA?9o&w>-->ljt5a+O9zhzr7_+v}B3gU}?hwZMqJ_?VGM_B_ zNm{UTKTzv9OUvZ_eV%sigW8}yj+yXc;=b*iXp$jkYcclj$<4*k-+oaOy!t*dLC;@j z$1H4m@vz4=-7TMsXm;+Ab&e2d?|SF=NxLs?!K4r4FsJ^HrSZy{>n*l+6S!lH_YT}E zT<HA1@yhRwP%oMFb!6TB8wvP04(gkK<sR9Ix6#N%Cy3#~XC*}LtFSm>;qst#qUjHi zs4)jn8}wMX-nliYWdiJYPpfkq#xeYU3RM;3GWcu+MrW&YE6TtuyiI5c#`9YgzU!{Y zm_3dD!Zi9;r|{N%Zh~W1$oLFEOuO&HJw-?|X1(*{up~P>mxPBjLXwY7Np>npTB)l6 z>z%vdK6ab5jf!?TRhMY&&2N({dND(DQGm{=*-9a8&0S!`;iW8>{$dDVg)aGSem{mU zGCsIaW96+*F+M(y>@qMcFx=((SQd&B#S#mt=2$TY8J{~6X1u?(Jf#&z9$5$TrjO79 z5=CEOjEKe4B_M#6%{U&Hu;;MHS!?cS1??b<=GIU&`1W_u5YzG@3iH+>Z-?_e^c;4< z9cS}Gu2qY-96k&=BkR4EJh7gpC*iYMF>4miJ2=kcjMe99{IxACmL5*m_<8TfXj0RN z$*0y7+`mCKA@)~@BfE$ctc99CkfX6Lctq!8Aou3=7Rr@k4rTgXvs5;<E68RkOfy}d z963whre2a;ip6$^OHt$$XFX@P)N_f%(*rJbI-cikE}XAmndQKWQ!Ktkf7}NbhO?=6 zn=#pfpIW*VyqtRvvi2mjeEzz##SS$h7hB6#nw={Onl2g}C3yvsxEr{!t<TeKnCR%) z@x1Mqo2ks?=Juqvr1qD0bY(RcqKB0!Yj#Z18KPj;t*!CJdh!)C-_*?NhaP&Ur^pAM zB@sE9{{q|>jp)xbpR*FAbxn}=6_jV@g;uZL<3pQUo$FZYa5SRT*~ZH>T&{QSFxET& zjqnk|$KZN{q<m-(T~Lezr;gvEiLI00pyM(}-V5YNg1JiO)5T>wbz6=Mt&?9+9&#VV z5Tj|Msc_K3x0!}QY<U7|*7+go<Y3tK8D>2*7weOT_a<~EmoKxys17!S^Yh4@yL-v? zVi`Hi(^Yd($o7b`y(nAW;722Fe@ob*Xb(8p0$QDDruh^An8vc^FggL+fHn9yUamlE zp(XKfp99KdnmRXjgOB)Tji;A1AUO&?kSq^@g(7?jY4jpa0)5vg_<bfmsg?5z#g5EF zx5#of$#U+GDd!zDVAo`nLlZzzk;YtV3Z48O3ztKKl|cu}ULHhYa;?{==}FiT=*e3p zKUX+0{Bs0<eP1jMc9<)unC(rU1%8wQyGvruPqk)!xhbaokefV8rM-zQsd!G`=-32a z`UUFkbfn|TxiiwRk3ois$dk*-WafSDT+|E?fgK5tgG_7Q%p$tJN$V2xI6Vk%lDtI- zN)HsSOi3ghfX;*|?S>rEIM<naOkg}cp$>Uw`IWqcp4lo3et%9LFg``QK((OAQ%E^m zq|#c=L&&Saql{j(x>m|TqZcQJf6+KZ!u5)yz*LNN2_$BlCW9dQrccu@%=!_vJpwhz zsvi`p5_WkKOh>99JPhvsm#Kx_Cq{B#kF-{&@D6XGtx8&vum%a#8l@G^=P`{#B&<*z z$4gIj0}@(2`^9k$jv38KR`{@b$;!5z0Fh0<nTu@$Q-*T0I2oGUsoc<6NP$D4ik*9R z50ZQ}<jb?Pq|9HFqHc_tC1k4-c<v!4i6HJo3zsBTv9f@J7$;cruM+-${vK8TqugCZ z)&DS3;_Lr{d}jUqqU(Q3l@P9fA8zKw)jt<(3F;rJa!1;al=%92dLiq#Mc1FLN(k2< zOjQ3tSZ9LzA4PIZ`;ih~{}<#l+s~DZxG`z!|0z{Mxc+^(6BXZnu!ZXv>pQtBY-*<= ze$Vj~8EL{c#+wG-I(e_^**9^IgUbr60|f`|Nm;H)upvQMZeLmYFeI3}7ip7XC3=cv zqXjnoGu9^mj3E@QlUtDzuao!4XX>PXv`+3)C4_bIA@0P*>0}1j!a9lA<cTVG*e1BY zh_8PN`ONwUM%RCeDj{6ID^dNQ!fj~v_vfv<`Z4mhUH$oR-!|NC>o^UrBJTDaDO6t$ z)W?34_x5@)!;Npqbr#l=r0KUBZp2bhjzR7rSj1WIbg^rv7Kqk)8tsS6Xy$$qk8E&q zkZ^9q6~2V3%;m1=d*v&UN)ut;=8a+*FBh{|Ib+cp-Z9?(IvDj;Yzw$|4-oT|aO!_X z)77223qUiM*J2fGZj}?4TM%1r*MCHpo3EAo8B`fvF8wE(@+y>u-|OUBJhBU+9=N!6 z?-SM*cj^nBnWxS-N1Zy0dFFSJ0#+~gd$N|;Q*WX1ngspzpq4gxcb=STPy7bbNBZUi zs9sXO`x+E+IHA#Jq7d%DoXPXAOqSdGr(&&u&kL~qdGLGkipxMVO*{@LsCbrDruCCo z>~K!1!us^nZXA$sXZbx_FQ<^dL2xl#zgJtzg~%y=%32WFj7Y-@q3_rorf@XNr~qkf z1JUnMJC{NsjOxz>?P--;E*?awYyBRC*q=OzyL+aN#ID5r+Kt%9w1!SMBHxdJF?LKs zjaRR|Z799Tz>e{+Lh9%H<4~h@KK~AR_2t7|nhu9j9p9OU7wS9})kArQmp(zR`RkL0 z|6KRM{B;AxcCBS#DmFu=6<Y4H;LA36+&gsa^Sqj|Wd3f;;r5oUNUF!$z=E&rgn2t4 z3%@6s3~4&H!OW-NJCT?qaAP0dOMd1A>g|4yr_vq@*n_p63h#_L_UVC8wcTyU&n$@F zaj?(8&nwvR-X5f>w>p_=!TJg?2JPqCN8$b6_=S(4y+_R)l`jqN2Q&+)nWOThRmxwc zu-Dc5%K~7nOv9BS^3U~D^4+>7Kw(4qMz|91ERWx%=zM#sfp;YGC47=wBD|8%a5fAr zz@G|lNW7xY{P_jL>13)=R)%PlQI5am#uWUWVw6>T{iR-9niqJZ170?PpOxV6Nys!= zhJ5=~<VO~a8HM?-GY#=05Qbg6rx7X=T7WW*QP5||?Xw4F*bTdfFKf0}Ls8Q`9>30J z{%Fg&=i;>i+OL8be35UU3NNduu{#Rv=h+(&@B;(>kf1mUbcPJYFvl*EGYSm6tI9QD z4t_ZzEz}Sas34=tJ=f)!?#jp+A5I+Q6PIx3_}>V7xC-9mXPgEdz6zib!2IY)T+?(E zQR889BRz$pF@a~c?cMR+V?2Giq7S~~jIWT?;yc#-F1tYo!wkZez$W?4JMr0jLE<<y z&rF|U#(6hH#f=NnMv?CM$3H|Y@!d-Y)QmPZOsuMT`QiL~?pk}-(zOTk4?XzxrU}!F zFS!1Tr&fHwe$~jrM_*X}&AW3)9`pG*SKoH)<>!C$hu@~yhyUTcEmyw(ZTmNWKezw1 z;=g~l>g019tD8OV3~wA>SoGtK*WIxrSLIjbv3`8>`B#@;)uZYien-jOV?O*-T)m2} z_!IHP=*uvmSC#MPcTVs4P1-RdkJQe%{AmBT@855E|K9E+iXXTCbmIeW=d1c#rM^}9 zB3<^Q-hU5|TD)!RlediiZ}YXT`SQWF_xE$2`lX}gzCrZ`s{YQ>AKP?Z!!_>JdFOud zb;=*|z8|IV(9}!b-+9}bQHL+QD&wauYep&h%d=*BhyUy1QLALWU70BpCw|{QPsyp; zrP{gq>qQ-hj(=dp^-`|WV&ppYz_h3CJUb^ps=b>2m0g~bKXCU)IRh>j@$r=x_qe(V z-=gOyVrQzKL~`TkyuQ+oZ;}1?+lxNC^U<!GMm&<g>CMaPGVrZ@ej@pmynpQXx3T}8 zyWoP@cJx1|pzOy}dgPh$2{QiSf)AX3e{IBzr{$eF`mI6N{`<p%(YMxi{Pl&pyxqI9 z3v<8zHm_yrHT%lbZ#@6aN2X7_Z*bN5ue}%?^VUhb&p%KyH}}*-pW-b6rW`jt|L*hw zQ+HhO_?FX;JEx>I?}aC3Z*97?HSdb&_O07@Vr$-^9*#S<4R5;O4cidk_g()w;?U5# z)Al?+WVEGp`CjYfIp^Qq_2V-(e|+TyBX2D@7<g^r`QwNF@XYk)gXf1XzV6li)xMEa zmT%7O`qg9Se)jmTK65sWJg;Hyh1Y!k=)337;HU59d8(g><UBil{WG)jJ2L0boB8#< zBh~LWMva@i8`?~{r^vWpo{9TOJS}fxaoMCvMTRj{{7h1CsG0yVgG7F}Xs9vtNFT!< z#cj^HhTRu%<98*hFtIr6>H>9V+x-E|4qCB-(Xh|ptMZLO<{Tw{s3=^WUCdi{k!V*U z$jJ<T8UGv6psV1>|G%0KsOO{;yPns8`NEI7+twaFoQi{o%{a}w6!A$okX(m3<Q3cu z=|KE%#LK|D3-ObPzeL=QxC?Ho>_OZcF@Gp=0AiltO+=iBcq-yc5zj!}fcRR(OA+6X zcopLHh&Lks6mbXQZri|z`ZEwOMO=V*72<Nl8xhaje)#ZS#P=gk#uEUW5D!Lt0C550 zZdj5}Lu^Mp5Ahhp%MedR{3PN=#O;XJBi@7fBg8h$msua89K<6L7a+bA@ifHOAfAW# zM#Reymw$|{BjQTLn-SL{-ivq-;$+ON2N4fOY~6`+5D!K?4RJc+d5A|KUWV9$Meb_E zoAJoh7Q{Oc???O-;uM@s+wg$e*@y=r9*1}=;!4CNh#L`KhWK{G)reOkPQl}G+Ysj< zK7hCcaX%c81`zXinCBucLVPFUYQ*h`7b4z?_+G?4@eE-b;%dY@5idfFp8z&CARdUg z9dRz=orud2e~!2o@z01CBj$6oD-ow4wqw(hig*O#T*T#wix7tpPeZ%}@i=U9pG3R? zaXT@-0kMY|F@NxM9JW)bh_@jghxk(LJ1U70w_{zCjt$Bl#6^fL$Ttl!e@`)l_;$n# z5wAx4F5+#7hwevv@HMT4h*J^Yf;bN`zY}*U;&#Lhh`&U<6!9U%s}QGr3%f&{ins&u zT*Ukx$k!3u5w{~Afw<s1^fTge#Pbj@Lc9#|ors@A{0L$@)?F_k9)Xx&aONjzZgGi` zG!L)+;V)I5m~62uOd`#WFbZQ3CV`9Iw$V*pdco{PK29d*fJPikD@W4tpkdh;M++gT zH)!alI2!knKfjZh_6TTSf`(Vs#O2$7q;8#*xf8VSIw_Mr?^DoDPEaPFYWuts8qd%v zOQL-IefMuVq1A#$Sz_}IWc~P{ABsdrbT@x8)a(iP1<+Cvmi`^%2}-r>JjQA(NKP^e zQ><O;je<U1>ubA>OX=FBK9tg}OMOF1_b&DGz|TJ1i8^a>Fn1xYWBZ9)!fPSV*$HR^ zLE~K&e7iU<A75^<7_`{_wS#pd!UE9v2&~24&uUY4GiAPo4_r|uz93=&`fCWTV`W+j z+Kr$^$uxF;SLRy<+TEZH((_5$M$jGs4Z|$1%ns1*1P!loilbQ|)9sy<X$S3o(8O7X zY6ps-EF(Z;nNhMxS~+MdL309!+b87=fwl>>6LlJ{X#-0ntwg7di|B0?=oZ{*glmkH zi%ZJ05wvemE?z|yN8^5MKWL}xG%1e-^WD#&!97Ifqda!dSSEifIou9S9?q@HFlQzj z%baf)c0%JE+yEMFG+eHfhx76l&=U0{=jzuxDU<W}v!I=-mx*iYhI4uwXo>uo1rF?y z1Z9f$gLXS;^e<sKW&26b1D$IF{`Pn*o%K!weLvO&iS#oMw4I<O>aS&>ebPxjJ_oqD zlYH%<y$xC-9qa*ZYbRy$eQ&RKLQ4hh*#tB@+P4~dT?yJsTw6wHV|{$YW-F2ZY^-k* z*<K@PIiSUk1?uj0NlT=k)u0VSzC?Q64BF|SorW^PeS%B&*Iv-lK*Q_LA~emOQ?Q20 z13lKSQJ%9w%LVOpJuj~92yM{qXwWXhwPl*7hr;>RYc1Ca-{nL8dEiOZHw!_#1~jq3 zi}5q}f_5)xiEOtGv?ZV=l6M<u3qgZuk+$NJZ8-p1T_^eap`w|fCGs)ppp_>m(~dIN zLAND1Z+rm0!s63(P#n=gE%LvOV|ez3-ZzwKF=%;MM<!~^O3*rBXNhR*K-<v?Z3k#u zI-wl`tvvzF4q28!KN~>%4EAl=thHrQq%9+`zFmxm;~T*r+m>?B_I5%GfmVj~d7^wv zKwF)FW=9o#7HA;W?Kk1tQgK{lUe(50J7oe-qO$i;rcP)!boLhP1rq5u6}0vQG|I%i z!v@eMQKp7?naHyZJo~^?PM!rCj~pB1F(BI+*mDS;rSUwpk4GTKrQpHxH%{LhK${5K zAe|=X6Ye9*K}$4;b8k@$S|T6A{YD9BiRK~hK?*=iq%ZDEMsz~sUL^~($-11l7WNOC zAGG2GH128c1}%0h+rc{m?OO~Q*VUF6HQi};!M*b~@FbF#`xfs1qI@FCqD<UFecnkv z?yEimEmkJx<6g@GJ1)a@?AYM`YXoSqK0x#j*pudp6>yST`?5UpE=W)p({m9TK%0s< zT-QW%oHFU~;iP9Q=ir5(bF6Ji#<NM*6<v(?k|6bMU92B>F}8KF-aN>-uB-L-uExz> zdA&jW{h+J$35)SvS8I#KSYYAxBNkHM7JvV3vCg*|hs5<ttMx^z@wAoH55?czR%>fF z<JxXadAytT)ow<cxPGWRQ`)*)|7kO}bRT)Y&G?Ir*RR@`eo*{f(!+XZPh(||k>B+& zz7p5B_9XRn@png0>zBVV4)h$k{x`-;zv1;~zae#bGXJhhw*I}Bu_1ZnExn8zdhz<1 zUQGX7{B7)Qy}pmJtoO)IdmA5$>+AZE`keUN-p7j9f_3y6`E-i$XbP|2O(AvBG5ovY z80+8q8mo^Pd2L_gu(*DpFVny2Yu$gW(Ri$N#j(b1$MSmbvGbEaIM~nn%yGsO$92K= z0rB_c<E?l9)_CW5>)zvyUE=!3-@ZoH%_msj9ANx#g7r@Wj70;i&kisiA7Fj$RO6=s zzsKJXPRzLLRO6n3)*lBN&kwY|HPF~7uJ;eLB7NWQ4tGH`{}k(2rx<@eW#Fr)7;WOZ z;}nvYpUQ&PoJvs_oHhr4uTQmhq#F09_H0Wvo)y>sP9^pFLDm-s8BY(gelp1TU=Xb1 zXa0fh95?VUSlqpnWc@nH_?%{QT^IbJT^Xk&{kAu*hgd^Nqphy4%}G6PN;2*t4{}^7 zcn&Ao(yhacVx*z0%As?x21~NO+{O5kY|nMEUf0#w%mN#`TA#NV%e&(8m9ExzyBh5x z;Yo}28>{gq*QhYI>#fFhRwlr}7IZfjc0=+L-K?*4GydL<3CsCjUDUHEsV^>mXWiJt zm}<Q{$+*gj)x)2YtUo6iuP0fzB^hsv>+h48{+BNN`%srByKb}@pIEIJ&<Cv9xPGo% zm)E)(Yr0u+{aLq@@%R4j)_vWLm%3Zu>~6d+u5YnrA^i&*i#TAj-qpic-h<V?&;tT| z#6J+{n8DTt<6P_SjSCT9&xRw>_B(4962x_Kf5b!ZrT}MRoVX8ANb-ruaE^6~ap9S= zPkxtlEc=B0u{g>4rzGR9&imu=;lgbe+mrh(n4u^KpQkWC=s)(M2$i_tbwI{*WISKS zunR#yB;y{~YvR)AC*zZ4e5Q=glkqqipM>=i`EQUh{E)z0M=||r8LyXd7S<$8e@w>T z$(Y8#^ebijCmG)(<9B7eU&dVXlb`RFX3Td;Gv>Rf8S|adaMi}WFIjjLW^BY4$iMyh zaA5J(78{B%`rTw<5k|qa7A_x*&1V>#PQ`6Kcw6y-$i_~*O&3iyit&d&9&h(Da4XO# zd(mR!bpc|YG2l4L;QYkS`1$bRz=9VnHgrjVb>TQJ(D9tvPk{f`r$6qPs|=r`Ee8>L zwpnZg5Y9w67vVyLsR%O><{~UYSdQ=z!gC0(BfO9BIl@7No;W2NfN&<lxd;~`OhuT9 zFc)DF!g7R%5S~ML9pQb1&k+tH^u)RM0E9CU&PBKoVJgB*gt-Wd5SAl6gzy}~>j>{7 ze2#Dsp(oy+IRN2IgmV!tM3{;&6JajGB824#4<S5<@H)c#2%jSyMCiE=`4P@UI2YkU zgsBKK5#}NYec}D`=h_F4oj7@jJu59ME#01xo<2N1J3YrfxCD2*$GJmh`Z-xca5DUV zV0M@#56j4oF4=KTMtsTtKVd$yHqwG~YC`VmMjC!~+>EQ$mDi77y{<LV{DF`s4e#nd zC*-b@2)_%PHl05Qd5*WzNF&121Jy=a<s3iCGviR5NtunGSLeeIf=vP!XmuVR`OS-3 zU&u%kKP;Ap-yp}e_(|uqx&VJL*GThJyJqmWF<jM^$f*#D#4n$F8Y(=sAs2+lLlJ2e z_<5`#1TdqS6<|jt?iz0exC1B{f6_#4X36*^b7TtC)Zh^aHTM3KnCAfWJ@gk@S6ggZ z0AkH5Gn2;8X$aV&3BGaQ8wY|NSbg9U#fNKJ#a9i!YVf5>8Y9J5fa~wzV^}Z!+hXwX zY=gdsAM!CwM4*o&AN?Nh_3*v{{UEXk9{RD15a|2JM?c9kTlz=!Aum6FK%kE!AN^a$ zxrVVEL{%@-8D=5SC-Rzpalr_jIpCW0l9#gaGZz6{DN{c9)q?X)zpBb39|!k*jgNk~ z41Dy@<fG45e2Z~SSy>O~fqO5&Jr7)~kE-`NBq+Yk;M)v7p4F)jFBxt`P<%_^XO<LT zFOa3UOj(JRX?&b_Ha3`h0L7__zEkp1FCF0P0N?FlZ*R7hdG6Qv7T}rl1y|!Nh8-k7 zs=OzGDId!-te7{NFn>`Y@}v0v4y^ci6u<LN*e@v}XzH?E^0BHzgRHhgH%Hotq%hyR z8Xw=woP0|}cAf77jgK#2?gu_qm(KSw0Q-#fa{lWNz6}tVY5Xt`)wBnJ?I0h2$6?2< z@#T^CYmJZd`A+b0Y^o1=7!D$+^7;=!y|=}eC-Q)a^|M}nb7Ua+_GsnV1rxdeX<D*D zzVoyMb;)%E%Tk|<m*XCjA|g>;G8TO+8S_>c#tLlyPXZ=C+raV5F{GEb<3+4bxzUs# zRi46nOut3USLyIW|HQbfteJku#AXfjHzsE7re86!Icn)YOl<a^S_gF%p3U?JCS9<J z^^}R(3Ui%hVy@E6acp9GGWaEoR}-6t&T(mCaq=P7eI^zsFJhdTSllr%IG#*Q&Edmw zgHKnn(!=<{_%N}!lO)E0iA~Q=dp9vxhwztSY_1~J`mL*Ryr9!gdC^s@BF*)siTevW z?ass}2u%Aj@c@C<pRUG<VY^_O#rXO#HNbI7oP?Q&Ih1@!SiKGgR`Mko#~IncZBk!d zj9$igVDTPx#7W3K6*1q)&YmT%Mr?n@)GzsUJ7WEuEN!y=tp94@9a4Wu(B}=ndi{R} zUID*{2-elhc#IeF!}4IZ244%sdRRaG5nlkM@Q38Zqvb4s^h1d6-eR$F@4@oU1|we) zWzzc?BY_X%tsrLpKE^m){HS4m&HWIt+1@^en~Xn3Fu(O=(eNLEZ5rm|>B$<t4LC)^ z4*>VmFu#k{U&C#{12w!6*skICfm1d78Sr2Y?*~o?R_(RIS>$T?IN&@D4+0*i<sSxI zq+xz>tysf+T~?WfF9t5xup9VN4f}zoX?PBBrG^&(S8Mo2V4sHX0Imf_v(UeNj7Nbl z*YFzPIt{-A9Mte`;1ICt&!2%CH2fRb$vh2nFWac$G~fjq9tpe<Sn2m7;3XOk058+< zLf{n|UIDyP!+a^pDh+P}UajGsz-=1-9(aR>d!k|O8ny#((Qq#C4h<IrcWAf*c&~>4 z2z)@pe*(tAO}M@H1KTvrHD-#2Ujy#1;SYiB8vYh|u!ehLib&V+$-uc9P6N)<@cF>w zG+YE+tYJQLRjy$l@H7qo5x82zEx@%J{xfich93lO)bMKHg&KYpc!`GJ0bZuzPk~ox z_&ea08n$9eT&3auz^gUv0B+MT*XtWJJQcWI!!v=mX!wu7J2boqxI@Fsf%j_oao__Q z#>pHWri%3c2f#KBe-E6Z;Xc@`^w;nyz;+D}10Jm53Bc(Zz7#lD!!v>NG&~2mSi_5e z%Qbu_@H7p-09>u%cYtd(ydSth!+kNZ8#T=Pg$p%25qOD)rvWe1Fz*kp(C}5jD>Zxr z@G1@e1$ecF9|dmH@JqlOG`tzOUBjOMZ_)6Nz&kYD6HBNL4fhA$tKl<%4`?`-ORk?I z{XHJorr`?U6b%P~`)jxb*skF_fd^~&8Q^pc^L}5hhC6`sH2f3rI1Ts360lgqrvaC1 zI1_l9hMmCG8oms;R>Sqc4H|9%Zq)Eyzza3}1n?3KzXH5W!|wvG(C`<)D>clU(yKJw z7faaH8a@rUO~b>1H)wb~aJz;t1>U0J%Yk=j_&VSY4c`I0SHq71AJFhxVB>J4|K9?( zY4|hX6b*B;-Cx7~f$bVT19-58hXbc;xCl5`!<Pc*Y4~#BaT;y{F4pj$fy*_#3V51^ zUjwez@K)ej4SxsRpy6a}85%WwCh$THj{#ny;R@ho8lDZjLc`YruhcN^YZ<FF{4nrp z4Zi~1rs2<kH)yygwlM7)9tylg!{dN=Xm}cMhlb|?@73_lzy~z^I50Lq;rX7=Y}z#Z z8E}e*dtwV!xXohYeN3*mPY3P~I;MGp>$Qu`EfxzkxSr?wo!<~O*K-35uJ8HPP=&dE zpN9JxWr$fG*Xup8KU3?!ak%)6hPhr(23GZPy`G|BuGjl%nCtca8s>U^poY0#w`-W| z^;8XWy*^mOZNO)1nCta)4RgJorD3kub2ZHMdLA&`w$K;X>*F-c^?I>}xn3{V@^ig@ zsfM{;pQd51*DE#5^?J32xnB2anCtag4RgI7(lFQS4I1WpeV&GAqQ7}RO1FnbjgI@b z#&Wga>(a~MPjRUIs?p_G<LMau8ZggmpdKTsOMl}-V6~r4>T*18_oBbmezQwI<6B_9 zan|JTV|2lme&g2`TO$(s!k*HA_k3=#oeF#^@JisG+by<}CEf)5`iD60l=xU2%>28< zVxwCm{~}=imKo1OS>6r6FYLD1coIOo5?FnBTCfv%z;_ngAY3#55c3E@u<&wB46X!j z!+x6$I0pOK`WX5RG59H955^z)S^kT_)AmK$`x&r0521Wt#^kq}M;Kvw1_IxU_A7Z^ zG4w^i?V5al2Il!D3+Rgn!|wv-@4jJ~EA;Ug@M7o>M^r+-eZXg;{Iew<Z628bksrU6 z3H(2wMD+a}@CRt0(*I`QrJBCK0ABd1#m46s+5V|Ga<Sulfc}8`^P5L7EMAnq2>3zB zqx%1F4E?t_qS?RGVw(qkmiK31bsm#qoM5a3zE+dxIpBfvd?p3{Z9VY5T@iWD$Bn`r zT7OgmpDOjk@}B`N*X-*=^9U#`Z$=FE$Kd(E&Cr)>&jH|M*c;tA%g@3Q(i>X&4+8V8 z(Mmrr0w=@&vO>zg2KY%$|0iolPImC50dK%~K{tx>OM#QWwb&*}{D{WS{#*yl^Q*<U zW`Asrq4&TIlXk5?Pvs31`C<NC;6}7Z)%T|u`cpCZ9pFnKPdW0k{k`Cb@{k{9BJe~y zCi!9cEx_Bpx7ZxGCVvkYYLVtI&ZXhX5BV<weh2-l?D<Aq?16q)g2wzm0oS2EOxpsV zi-XxV=yQX_F92`S?B{JB+{zEj-xY%o0+(;W{Y}s~eh0u&EyDPr8b$vC^9|=Lg82*R z$mEB(HU>Wh-0j~M8{DMG-`6gg5s`loFuy^;RN^tf6QLiDXV!Nc@KUWm-Uim~e;cs+ zFh74cSbd1c#$Ye-g_{3e4$OB~=Yp5^@=Of>zknBkU-^$+z}LWju&fsH6>uQPj~#J2 zFkc|gB-VeWMkj6sE<<}%f8t>q<4*KXCTMK$4q)DwF~>_E<FlCj*3&Saa6cRD;5Uqu zfz^lQ4+WmF!(uDJHOspNxD5IXN!%|La|Gr~W#5y5)ra}3fG<RTnmO|?0A2}sCrf+~ zcn8K;5-|1En~qz4nE!0xmD>1n#?bGH!QTRJ(ENM1(^YAv{_Mc>V1LTKF9UuA_Duy) zzFUF$lb=(8X|E5&<bNKxU2E?l;2X8^(EAM7pVTMiKOI<oh(`lgWBz8btZx?Z7RYM> zX8Y#>7kv|%|I}}Yvkc~63;IUTZxXnd@dL2>kbdl7%ukv<jsfO<L-O>+ozo)VBIuI| ztp9o7Jk4IW0;><pI|y6`c~$?OHU##B`9;lVR|BgL`ELaN8UBT4#PYubenRtSXPt@m zqkgNT&jMB-=D#uq-%1Sm?1-qJhk&ov=7Uqt!t*(3FY~g#bAi=|X%l1c;u!oC@Bytq zPv^jtALg$H?hSiT?SBk-4eVKsx1GT6X!1-wN98vAD-?ql0G|qdDu2B@hCVnA52ZnW zYP{3}FV^(63b;4$NL0i2_8J<Q-|d*+`vZRtd!*YWeKauNTt=zbo<+dwL;NW42N(|= zM&#cHT#5Fm_MMat0_f~!(#HT_r19TCy8N)byJK)0a1G>9{ry@D{R7~Qdn5MV1B%%I z|DyEmq@v`9<<F17w*qg|#?MONpW$DM!Ap5h$WXb>_E!V%K>w-nw-(p}?3DaxWTHIS z6Tdsf{Ih^}Lf?E2o&I1s@b&1w%OwBjz%#Y+HY6*$J|}P)@*gYtTY*=jzMc|)30w>R z(@WyKZ1flUU&(g|@E*;7KLV^i)W?Px{2_21^sVaqDu&)C2j2_C{K@YVv3+*{s}IY2 z2e?u52j2p3)5iPcVT#+-=ih++TKoEN;*%eiXH?XMf}#2uGtw#ymn*+)Qi-dmu(ZtO z0xF0?7v^iE%7ClN7ntt$xx^DME_Zzc-Vst$>%(hID$}xuXE+QJT;9qCqD;eljHa@_ zre=<+$2qY;R5PYz)C8wXAy%WVBejUfYeX&TjVfss5O`{lBeUuv^Fg2dYEQ*1B$bAW za$PQOKs?psnmRE%zi6DZxX6+1@p&pjbpgMtqS{?o<#B>H%i(eb`BIb$7umgTpLec! zMT!frb#YaC@JLd<cvXul-RlLvDlFSm>#Gm4gre+OWhIlvRnDw2ql!w!<?vY(+zpfS zQ4d*ii8SE#E8s{MjHT6qx)3=<J{{pbEQ%&dnOum+g}^wB`6Aa|zq~dpE(-Xo!WC=4 z>sP?x5ET>VWRcaBBpRFSlAbYZv^Nwi_SB7@6Y?-e_N<A;lS+A&Gi$Ox=n-<|%qp(K z^IImp&_5%<tK3-^yX*Y;YQD)VzMXGgj>fBS*k@>X>71ZO7=y>Jd?YwBW)<RburYzU zE8KOJA|+F#jCNPd3MXaFDuo6|1uHy$z5+n7Ws9@|PgG*&tb)MAK&ZGbP~pK7SVFkM zV7@Ql7j1H6&Kivu&Qw%Snm!Z8BFagZNu&I*T&bwqFkM-MBfZA!^LcAv|8+EhaJnOs zp3xa0Gm4OoR)i9j6iJUMDME-UDLYP@9FY>0+f0heH++`Y7l_LX3sK2Y1!(C}rR3DQ z>w=z&=%(R$xVreFhx^_B_|#mtzX}5|P@@--J;PI1<JD8Lro*`OgyBNt2~tsYY-%oR zj!UfdR>UP{_@m^=sE>v<Q82TnCQ3)ZQM$>fh-R;dVXug0cSLb#L^EeZ^JdJ7V$Pfu z%?%taNo`aQ<<#LF8ty1{I!p_WC*)Lms_<H>DAfmH{&6K_B1z}Us)<R+2-fNJOcbnR zM>=~6wiNJ3m6H>nB!x{(K>sHoq_el<)3d8Rb)mVs7-Cq&r`C9?+;OQH{-B=K5rr~p zqF|QCRa&TX55qgHFaW3PX%39jDz}=#W(RzsC~8KhlpK!Q&Z$-2Bb1uM-t0VYttZqe zHG3x75M>*g6Fe8`GG|1?%vmvbYz!`o!Nt)ym+j=c!Sp(^ViQc?5SyAYE85<oP)1D@ z%(}uIiqd_K-*ZJM5Itnlt39=@T6eI*?bF4{gs0K5Bi$XU2?T3lK2fU5ao5zwBuT?g zOn_l0Ac#>DpPoxYiZvjiz1YO8*o4eLb(DIuXSicVvYZAJrbEiObm95qQ#mB#QnRn{ zOpl&ivqPSs&#h}TyUJZ3j2c~H)W)S^bo=Y(=s9y{d3>|IQQeqT6Non9;nO`8m)E0b zqIxyEqQ1@>)qga*sKS+TMkR)AF*X&7iY;JR1s0^9n6~=r@$|UfhD=|5gN}#6it@3T z(9;|aB$ibN>VqgTijWS%1TJEtjX9$>YP!u1T@i>*$*FexW<)1t)_SAnzuZx?16A(I zx`1BcunL#YTU8w;-S7%m@bY?hohM4P+-g^`56iMydUl7S8pA08A)Dz@vzjAQBqb_( zxF=DnBS(2XSB*OubX7#P#*tN97r>BDNXzzn>cuW2K}L-0`k1mZk)*dOv!*^>$GM)` zARL=ND$xOVU*p0oQR~BMS7*;i#}^pTB(ZO)uk(P(AB@s)hNHrVCSk{8HW4patn`|D z8MS83$PBkWAA35s)aUV6g{t+M9FEa}>beBE;naM!)oxc^g)c^H!g5?S6&dj~$1K-q z%xQCL@UF<HTCz*aid|(M`a^dJ4T>UWkXT%Ym3%BEbJkdQu&f%S>VPlGmNFdmelNcc z;&Ih@+<t8G{6TMsLvU=-q|u{_T&_YFmhZV2jVdV|HM+=|QSM-!SomXC2_0a^Qt9IR zJ7;@bwX>ZW>G7%f28w#E*#x)WUFE40o7JfD9oaD~`5DgH*cp3sa&ywEJRw(Yg)3BD z@1K=6y@AQtox0q0a4RlsaiSBjf#vQ|7C&jsm{Mn%tBkT^6JIc8;;0FQ`N$!5^WmL! zFyN|&1MqpcSBVqGF}zT1hdm{MKuEl-jq**MSQ^T8l}vYqu<|P`4VH$oYTe$t!ttfS z;)$6Pi$@K=NVUu9EDbr*U8Ta79KldwX|@YVSw3_Zaym=hBE{)&;4PoE^&yu(Q0Z|N zBU7fTKqLt+A!TlvlnSqZtHf>{rtQcY>j@R+7YERXz*C%Jmp&ea3=0?H<the;5rsz- z*6F$2Su(yf2Ri!|8ky@Vb<KAB>d`<77gnKYp=-K3=rL;CUVp;QC=ixGp&Z!jaO8St zV1RqvHH8&ObMQC={o%&;ys*OQEO0ny=jDwpDjc0(Uhc|F%Q733S>vvS+T(gDrbi<E zS(@Q;I-S|mJyn<$X_`Loba-bP8yuh4FY9n+q@{yc(cpFoR#$~PR8gImhsyEYEtexL zll8|Ar-`M(d<XR%j1^UM9ZbJ2a0T>&q2}`iD&iXCzhY*wFvob)9hUF&V41~t%NAGT z-MbY*N-qZqw8)C%%A4rOEYr`Cz#}ZFH1rGgmWg3%4j7DzqFlVtSDYhdIB-=7`;ENp zkT2DBmhvn@)6!Hr9r0Z>srQOty~H$xjK~lX4&L=2+IVfP8pOY_Sq_IUaD}JNHJ$HJ zbe5RIrN&z_J{RiaAzfyv<JTEMjtssGnCk_sm_ueW9O-mUlj=ioAH(1X_?BczN$g&K zg^}&Zfm~jke_>c)X2J51W}PjRkUh3!RPngnnmK%Nv7><V2%RJj+sxM|7X^Z}#?qkI zJY78MzBLDKXO?G<>2xACLr)uat662^GR8oP;_(&eJ8xwr{+li%?4^#9qtAk-bo1%N z*Oxk7CBlka{HkPR<VI}FQO;g=;X9NO>QUESjLYssd-47!E8>6nR&LWChhvrC$J(!< z!I6pmqUqJ-q!c&bMSO6aw}wYH7Ye_IC37P7d1hz>tTf;YI=!(j`ahkC5<1!_cSw$3 z#nz$>KWxX^Y%3=D|9+A1s|h>|BbB-r3!qK}j&dg%bSGjhICUbH6yYpKXRKxgy>o>t zbyd}1`ESnp(y<?5#mBYBuhO5Y;P5n)TMBnyjRPnrofQ@_*bRl^40U{zwPj?M@C?!1 z9F<@<jYT{hv$w)&dc3%?99L@&YAvqziaIe0j<V8~aXHe(Yly{Lh#fhY*UZ_Xz*Ea* zuBl<HoN!E7gEMw9-y~Q@)%p@G53%aT8nJYGA+`wVt}=(qT=Zg9nkl;X7cHIrLI$yP z#_}z!73{4l>OB00J3W0QnbVJ0^TzepQCF$vTfMkDSy(#URWv=_=|nnbZF$SDvqeKj z)@Z(L0KcSw8=K~05(j!*S(=^TNU$SOoha6^a-;O0F3)fi;nx{H4o<!33;Z~Qxy=eM zQGQW>*UI)E{q9N^-69rKo$i#lK2R&H;o;qBF2)INtD1Y2sg=Is_%(jmdB}0?$PpeS z-|eg5?`B~C80IKdgCV?Kt?-HM3U)<waKp+<uzL>FC=a0z*$nPU)sA%14BVN;{!^^} zMemk|eAxfS@7qd{h67^m@36o%w|8aOF=xB5GD|R0e(4?n1`72h-XLR#j@Dsw6;q$9 zs-eM!BbZ<SzvO|_64&hX*d5+cT;YT_yf`6<KMgqY>d4IN5ii-`yukArb2Ewg!4;Tc z{+>u;js8d5>+A|2?iMB3LFyiP=FCJjBytxS4j4;eexbj@i<79(D1Q|_jn~|LVxKV4 zSu%D?(WKH+{^m<&X{gRyJBp{zll@nC@rXi+Cx}B}^j3a8w))tJb6uEY{(MSdC7gYw z)0yE8IcJhG{wVfQZjHOar4|Vn&W0dWm`ZWn%)>>G*v-N_cmqM<=(Bim!?6-=S8m*R zlPTsg0;T`&7iDT!8`&Ua76d}Ha*Wj(wXW(bgkFL<?%G;z%j5c{Rvzo&kiEj=i;-wp z8MW@OQbYadJH41a^{=&o&W3LlkF*5ln3F&!`!Caoe&GRunEtBp6og;(wdz^;oEXCl zFEW3H%}g#D3%84rdW2SnPo=Ue`0$X}pGRgChXd^>#Bz~VAx~)HmWXn?FAWU?Gfv|@ zoiDmH2alm59K2YBIZJU?i4|FCmJ80FOS0%uVs2IblcT?^0PI_A<momGG597Sp$5yX zxDvxn&zK_4w{b4!E)_Fe9BXHj#StAbj&j9PhME6gZP&KjG!R7jq&Q9!=dnU=RRN*G zl?qQ*QX3R0#ZnTW{5@x8_Tu%fw~e*DrjE6q+c|S)+&Cg|^jiY5B%X-cNU9gHG_GC- zobylHADSpl3{#!Vwq#`A<d#EyMy_Oi6SX6(ZbpS=!P^DX25a?=Slx@K`Alit=)bJE zMMvW!99ssnA^3W?$+GLFF-m+u_Eq--A$-~Ez6qJC<k}!o<vvnL2L`3oI4OB`oRcbn z+^RR!HLw~?x`{Z;0Z9D3SN=eoedk9&BPVlZYCga~LrqT=CxQ`aXIV>^a%j4jH2pP+ zhDIS?J;*v)b0|gBPHse=v5*mY#zKi%oiFMStAgV<Mj)bgR3i+HrN{|@Fj>}-aIvh1 z(dl3Dt;v-IQEhQ&>W-T<KGu{+n+@^)0J!!qLa0ndk{l-!q1>5JL00Jgm&&01e7phf zo5RH)ZMtg4o=6`eR^W6xWnC%p0B?by)&5j73TgPa*)}zJR*JL$$(6)#IFF<34w>Fl zH&oyzH&-um!mj1OA+d!Rf`_TPf5ab!nkPRX7Ne2xZoPC)YJ=1Hcsw?Bsl?@>;c6Wo z0+&0o0GzT!*S;gGOpOCAe8;}u(|hPaUMcToq3qL^XxM(gbLR~q5FE^#zjB`j!=RQ2 z<u;7Syslo6zGdx21%klon1Z3B-2g!I#{*<)R0x<`<h<#lLajylw4KGI-p_KH2%gBr z3QVZ5r}C_3ll`LFs8JM+A?|O!R9)qhW0HDAuE8WLb0%37Z9X+BSkw#1d9;ymsotQ1 z(n%pJd%nHwTg8{0P6s#mH{Zw}&3`Bnc1NNwj;a$NcGmU+nQ#J$>tl=eKBy<!jmWwb zfK@QM_3hG-C_7`HJf>gD<d#-dLN#s0^7u+2?!@>8TTp{?KN#ct0I^CJqg!h;<zUR@ z;>)=Z(5v&STN0qk=lvP4?y&yEtp4M6_t+f|%?<tkTfJ6ux}psbkZ-?qy+)<qXQAX; z(2fbtx5YeyDIPM$@T4?RGAJ)2o{qKwNQj7vko?4U=jQG+%#Al%>R+dphw=r)sP2fB zvMiIH;hNTCqykqb;%c5}ba(OdkG_l^=fe|Q{VAk1Bz6EYYp@|fdtedy8Y(mnCo2?) zhqYRI>agAf?+y2z^0xaPcx<ONoSWA0b}Kvm$(=@mli8I*8vv0%a9SC^jgtZtA)>Dl aDJWJGE@VI1tZ5Aax(8EcD#}OnBk3QEU43`} literal 58961 zcmd754SZC^)jxg%yATv1B34w?Rii~KF`G9MMGPbhL_$cjK_jBeX0xyho86e*K%l6h z5GBS~>PyAfVzrh&)k-TWD%GeI!BQI$73D>$v<5*%MWKr3_x;YD*`3{Ncq~u+fBp-T zx!*bSa^}pLGiT=By?cYBbbLxmilIJzjJ^g@UArYU2l2{VO-hboGg6Ha#u>(NqaTt; z%V|O$1W^EmtVpvU^d${AXA$TKCjc{$#=05Ahh-5;TMUDN<&89xF49i6rj{b4p~5B< zX1k>KLF7YzA7~64P$2{9Y=dDU@I<uDPz=nVXyyNk)Z91MFz)-xlFE>S`WY1cr=TNy zw~wm-!hpZ}!a&W)fIryKm=O+TWSM2bn&{hvvZ<mh`DB?P2&{W30%b$`N5_4>^`|eK zKlZK1?ma$OS^mm5|5@GW25lUI4dEDs>QzI|{AE+qBEx99sl<Ie>o~97)?^sNjG;?v z#-<{Rv1shvrN<07KehNq%b6CmKKhsqBm3BvP91t;xvj7q?H_ZCHLa=dnMPV+<}t^o z6<`o8!;GfXafYSo%sx3QnwpIAzLue>DR)^m-QINb3d45CfIig&PH3|HbOw1TW+V7G z5GEpAMj$_>h$kbIA)qPosgSV~aizSbG_F9n65%SD#^`DU7eX~c4T2Y87J?5!edeeO z8P_7NM~EO`2#t9N^AWB?xE?`$nsBiMp;;y@l`$izj?t>&HeCM#;g>S)4#dAgxC>!9 z!rch>Bk=P8!mq`ZVLT+`hY|l)UNc&W@IMHDKzIz{aRh$;WJboH5wAvg65(ltH3-il z@bjD*8EX;$1tD*jciW}E4&E}h`J3Apge*lv4_-dVJ@9~U_nEfr<1V`E*r$H`<&qCx z?-+Y`R#E<>H*Wd*Hy`KTcx}hw{l=;lDR*^E`)laO^$(2gUvjzi$5n%uOxX8=yX%qg zTmPBA{SO~~<9c-Y%%_~6KJ?7#e|zX(CoUgYG_H1C`kBU8k6rGYxou(n>c8C|EWY{N z+CzW%_R3e%{ZGt(^6)q3c&@JgcyaBTFYS#bbDkRKcynJ%*3CaXcgAH`?=76v@tenP z>=>50Bm3uH{dCZ?U)KN0{*G_n{$HmrKJ(iH0|pL#|AB3P|KI`l&ss0ed1gZ>`}LA< zvZ}Xz{Lveq4*KN2nNK~__Sv9qx8;uf=)B*B`#+m;-txfAZ9gB|^#|Xu`tSTL<@+C- z^g;X9YpaI0KY7ResSn=t#v=6Ns~pe%=RUY?VcARm9Zw(M^w7Xa`HQ0-8JEBNPd9d? zI2UYwxNS>T`E9BH8C-mLf6keApLl=jKTi3_toyxp{bSgBf2+If;HD3be|v=E&C>HT zeitsdal!Vd-<WgOoj0CVJmR%sw?4A{KW}a7yyC;R2lZX_#VHqf??2=8oY@14-tv`? zx%uEW@3xHapI^JO^W%L_+<C&=3#T2odD7n}b<Douo`wJDf5~f``~NDqF?|2X>+XK- z5A}T=hd<bNZusokKOGdv`d#=JC!{Za{fjS_j9GihFxT$w!86mQy}ju6JJ*E<y!m9; zjjKPC!yk=SrA9A>9l4K&L4IggF7B_8X=hrj@pyRx{b$06nt53MOxQ{FA^tlJl>87E zSQIkrzcnHMQwjBLOlWUSLVW`g>bo#jy|8;XC**%PA^&>`^7vIke-|atb7+L+hyB@* zP~YSPe(MwJJCs0QpFl597>|LfF~l~6OA_df3FRvj<auX8{;Lz@@&15V1MJTi3H|Mp zz^^}{|4S3{U!E|2OB34fNNDfASo1<37bW1L1b%G^?Kda1cS=J3{Dk&yPry$l)VCsm z{-=cgUzd=-HX;8b3G&~TFh0u?^rbz4{|N~=KcT*#B*?otfxp2CI5#1GFrj|?APD@Q zrcX7*IL>$pPRIHkmegaBd=~KA$79a!jNx~N8pbr>m6E>a6vN2xw4~k+%<`Lt8^)(v z`IAt8-5as;>jxV~)o`oe%4VEwoOzsKe2x64$@&(a3_U;1n#yh;Z=7Ph52bAU$dWom z(kFty(mj^cCSbOA>#2tE(t9TV<aa*We`}W|m0~9jfZsL$j^*!%`iJhbq|%2V{kby@ z<JLDVsSb(%g7y}IU;3A%PsMmGeI9-#KAa!BPBM%w$X_Ape?b3-fnN?Q`TgrK_%7gA zmH!|eKF+(b{yU&Q!!&u=P%c&0e-^GEgW#uV@~**n&49dak@W{auSEOkmdPLLU-FhE z^{>!-@_R8F9<7b%Y6y10$1#4Y&p`hmpE0uhvyewW)Q4dZ<<CFcFwRANoVKj56?iB3 zSM(jgeIW01CH*3_--Q0>1GE1#AkQ+CSM=*KKbHR`Cci(T;^}hy(~Qr$&OoP^K_2NC zaJEMu<`u}VO4k3W%`j$Tyyy>-{xtMyh9>XPKf-u|e++}5_rv%Y$gku-59RZ<@qZuf zy@dWN`V%PsEb=RT83=hdz7XU0cyxRP<VQO7&5njv0Z)_oAmn$-HcKiMg#5h@{dPk> zYW{>zNB@BBlKv{_+zSDU`c6gtqd!r2kZ}d{BVE&vi_zYO=M?=|<7w#o<(fYIoC^fx z!3rt=+Yxtdj`b%U^bwjq6rYIk+#4IuXE0w1QQjuYdmx{3<o^%~&Hm2kc!57a)IZpG zGI9LC&$E!H>i@UUw_;6yccS9q+IZa#em22gG@(A$|03$!1UewLVxM(UFEgP(${rM* zV;I}D@%jtOXJEXMEaY8@f!v4wsPQbtB(2ffyBy;iL46Rr$Ugz|vmN6dk$4~CGf-aX z+gX?&U!y&x4<|xC0gNBiUX;(U!3J!JwLcv7S9MrYF9V(ZS%&e;#Q0^){A)3OYd(zi z_kK*=V$`Slw;Vm+s?CqzL!VN?zZ%c6C_h@WC!e6B6QO@fpB7;}GoXJPWP9hrULS;h zU>F6z_X0Pfe>C&d&!3~dx6r<tpYK9HM*w3O1^v&^$9<Z9T>^REi1JF_>mW#Fztdp9 z*Q38x8vn~^Up0N}3;EP&`uq;&XO*T;_XD@1ekG3?7{9a7A2ol9K&a9B|9j|5`it0O zgAdC4GIZ=K<o}1XACI9Q2ciFzAN%V;e^z6>E|7Q}+S`xuMKyxH9vAuA_*{#03;0#~ z^cx%Y@|yibl}4I2J|V0}w@UlK_Vz)4Hf!tI-@xBVm=9_^UI5*v>F;A`=uz}v+4Dm5 zs~YlODEa#w{7nPDs{C!ZIE?YSTG9=i9TbJ?y*M^7T$jytRd{Fl!x3+NQNSGzd&5Ru z#8>Zi*SI|HpvN09T&`KQp`a@qao0y&F2gm=83{y6JkE&S?(w<nU6FdXKN2qSIGk?$ zaa31Ka#ojAx+*hWm3EiIna?6+qbjFV&-QvE!Fi6^CU!>Z{lQs9c1QUn4^nCZ<#U~J zS!)80^0}Y{=gk(eGg24mxmmqUaFFLJpXAIox&5z-Nq+6wm6v<#!~RgP$mjLU0h@lZ zonPsUgzDY1ymq_UpRs{J$P?eCV!LCm5V1cPp@<>vijpbvs?4j>CPL<71g4ecR8Gtq zZ%V>xclm2-@Sm+A-5!_Dx+c`SXUfWWnO&nmom*X2)_cX^vgH^P<)LsS?5X$HDGnyP zF$6q3VzuSI`C-2&EIG<Y?>*>gHn@T`#kfS%L<hm%q-cu03SEv*xiT>&d$N0$nHDXE z%4^J?#;VFL4n@o+!4x?3MKzU5RN)TRglbVvDv^nm#4Cl~$_~{8aDihnS<Pgr@PL;W zdn2NEGH1z@icm<6g*}sCN<*Zs0SZ3L8z~7w^Dz+|qcHYVS!bpYhn{g%u{Rv456w6G zSHK!dg54@`hJ`M8-JS^aN=kZCc3F98euYfT%C7KwLTItR!4s+Q2HX+<T<D4413^*1 zahV|oA@#V#EIaC|DGCK=`e!xNqkJeh-XHKf%*Gu-Z~d(KrJ=AIXwkeqGnSr(O5K45 zcf?x~^hf;efJssV#kTb9;K3R8`^Dr)q^V(<+z^gr6}uyDDQ|(}_066+d8)HAYixbJ zd%o(Ho@URR;EkA^6npDJVSmKzx}GO12YI06;5rnXGP5KIlV*w<c`@u`gTaP?yWVNn zSef7s*SYI!v=rHE2(U&}Bm{-hT6IOCP<@R*2oWh_QRzfS`N*t{%#56jEZNI&q{j1; zp8#a!VdP-MIDXD>p1ZEj8?2F54+-clW>^?hSk!!1MKuhn*XyeBdg{Hk-e9CWE;$&g z@w%$rVJ~Y#u~D8-EqdT}%${x9(YV~PN+!jn>$Nx|Ip8O$Wpz}bR5n`kdt5LcSQKha z!|s56l?>#1urkzhG+drWx2xWZ#UnIN^d`p>YKSnKSQwPm5Iv^e#&P5is_yp8X_+$= z2n##xjARu=CFKu>y;zxG8PIyQJHW|W-GkIHcPmV5<_x=?O7IY5q`>T^KZKPGoFy!6 zj%;seX0(&lSTKWF9wQ!~=tm%2lB7ry;dv^tAm(Ce;)-XEJ+3kQTscjVjB>ez1Tn}p zE;m;fF~^;ee6S`)K8Gr*CJEZjk(|$zBb1E&|1x{7ClCsI%_-&#%aOyfpHw#HvY({4 z!Uu3TN4e0!>Jn#;3tgv6fR1vtml^}*i;A)=v71q)`F`D0+4YGbCl{*zUXscQVR~|Q zB7e541Tvq3gj!5HF=XGbl2QJ!OVneI9o4pce1$XQWanM88XH}8a4teYcfcQ+@0yz_ zs-_OeaYEZMOAwnHamo4|*6d2umW;yMg}X!EWa`{)$gsxfvl~MeKVmR|wWcX^dqQFv zo8`g8XOi7;O%2ZTV-zdA;f7i-O|A4GNzNqpNHbV61;zrV)7236Vm%WU17@_*A2GnK zH%OgecIj-^y8S`J<*mWe0YY^+Y^EbAC{{yzu3E__yMu7b>g@&6eES!8ql>Uq<49DZ zCC)I6UY*-tUjjub=b}2cAnFu291wG+%jpfkN3e$@CD2+V<pjK8O!F|kT&6hecGqMV z)_4~O!a;kJ9W_Zz61h0AG-wd0wCjGZO7N^;F^<53!+R}_lWAo!DOB;^VLjjFba9af z+efgr-K1Nbr9_!s>#hSU-9(olrdSW1S&*0`R~SySlKjiZR+Nk#SL(>BveP1zb!*N( zswm)fV`<?cXnLpr$M&<i@QNuJYkeA9H<iv4Zg9=?2Zb+#6(%;E!XtBi&ylOE`#*AN zw+BM=y!Ec?Py^Vl_@R`*o(27__s{Z0YO%JPJjeO8@TRztkd(x!Ruol+URF+f9{7@L zdPC5EEwmIfAy@btxf3eJmQT#DozKmay;$f?ug{J=%#JzU`R3e+DM8$v?)66PIq0=M z?mZY3rvs+h4q>=x6_q`NFYyE0h#ut7ebVN{CM~^cQ+o=q=mi^K?i=mdLGL`%uaQ=; zo7s$6jc)b@HspuXJKFPrccwNSdR?e-sui6#Vg-ph^+&T>#N4}RrM<IKjKW-?OvL&Y ztvjCKxJj91PD(5*(cHbmNlqARy0WZCSUHso*y{|aio2=2BJ_HTYE+NPfd~q-`2#r+ z+MM;*;&{zj;3}=ogin+$Y_~W!hz?gbfo0{4Yw!nZoGj`!S9TmTaG6)*```OZwY92j zFYBE!%-Ea4O*K8e=wftK!%`RRV{h}B8Cy1cbnCc!sI@e9N|le5zQ_k(TZqt}SLyRc z+(qty2gfBmONpM>nLB6nJL(L30%8vh_s-lM9j(&@RR?3q{bV@Xth^W3kf<XISH8R3 ze&qE>m7Ua!{N&fV8(q46`|tJN^g3tOxqLi(>3w$~(J7@1@<TL}3oWiO4^#iwjg>QU zl=~26SbN=K=YxA}?{IX|gZoCIMfV=L9^P6=zWGrsMl)w!X>OBxnk7fnT5M3owugn| zwzfwsWlnB--m(1=wnyfX-qFfFwpt!x5jV|2j?l4~1j>5YC}w5G@A}Qdd$Y&g79YEP zY)y?iHo)0+Hy@Kj1UdDxdCsDDi8BsY84qy9xgPyrxoFW29YGw9n9zEeE$+3Vo8n&Y zv^XMRr^OM8(7+>E%_B-cZa1Uaki%O)Kg<S@$0K56+~Yd>!`Qyq9xKkNQiFfAq8lac zDQwaY<PK^AgMo>vX#6&z_ve&7o4S=am*hn1!I^L+Kk{=_^XL^WNFE%dY7XzwN?Ywu zk5pcBpUqQ+BTjxtzU(4y@{~^!v*?J{f*s*~4tf6@zQECXWy-0>&>iJUn~#GO7jIGZ zWR3I1Rv_6X;3#}ECHtt6$STYFVN8TJf^l<IE!(|rN@GWtM?DsK=Ajyox04-5k8)A} zM@hf#D97917(KA;w(=g~&h-Cbz7<4omkawQrm!&2|J#~cWKZV&=kS~wmAG-EQ>aew zo2`U7(Zg^Y(Ii04j^JSSy7KAg_PF_9=BSu(W$6^BvtXh(yV6|`EBph^vQcnP;*NrO zC*}y&3|mEvoq1N=>r$X|uie1tMFVfotno%|d^w&#nr%!dEg4tj!V?9?grXu>He;D$ z&&W24s;b70D{=8sBv(}-OAa3I!uQhp;yYS>#6Mzwh+>}@&AjR(zVp^k<moGOGrtA@ zQnivyGiy?RiS;zIz7%`|?0@^m@2IJN{c+VNid3o@^V?+I|Mir+d(Q8Z(o=0ZkKd;r zGrtSZ-r?~+!-kON_v6XG{>H<65E38OlHHBoFOEJ=qnDpMACELzjXy@q_c20}o{9Rx zXZZ1?bgJ<b=xoz$udh~LAH!ocX>9+E$L9c9jF)8o$a%Nnk;gv9-+}qbJk^@o56=?v z*>UT1!zcvqk7p!{8Dpp?GKPwjGRDKb<&05dC1X4TUBwvRC78|_HP2uS4r>^LZy#gy zD8Lxch}JQNgd&Xb%wr>Cw7-xsI@H7%&!R42jAuNTGRCuyw=l-TkINY2S=KukKZ0k* z8ROZ}`xxVy%@vIC?C!&iefUlsV?1m6C}TX+x{5I#VqMMn9z1i)7|)KbWsGNB+Zj(Y zjCG9h4CV&LHavUG_)Np-V0;#y?PrXKeYY~kv+dg%L#DeJR~p7{#&~$Plkt!7tSsXZ zc&3!`d4_R-@%e^vi1APGtR`bVe2Q`CKN2IxI0GZi_(H=N$QTNc#@LQ$XBlVVI|7We zQH*gII>;DSA)Rq9p4nxbXBZbS9)<B{oR4Q!85iK0R>l|M*<r>PWBeJ9#`rV71mn*b z4_lWqz7*rnxDeydc&uSeZTXjT@jgo@I5)0Z^zsZBzT$ZM;jXT(=FJgn=XXf$JRa)O zIUKxJUA6J$uH1a|wCkMF_{%a~=M>@!17#R{cA;3;ImN`J?dcGBA~98cPrJaS#28^? z&nkh-iEYFy1g<0=LA*@hDq={^*wZBNbYe)r*i$F)4B|rK83NZ3Pb4lE*hgGWTqtmW zxQaMa;5y<N#5RE=#O%<XG=Uq5>xd137ZNuT@Ba=7=QI&F5$_Us3Gq_m4uO{vFC%Uj z_!i>j#H$2eM!bS}g}`?ZuOwb3@N(i+#7zR<N4$o(PT&>9?Zh(#ewcUzak;=Ni93i3 z1%8xxD{-d4tB7|I+XP-s+)11!@EYR%#D>6Yi4PI)|Car4Cx(?o{{>!0%<uK==@57W zaT;;Ez?+DN5w8-sgV;vALf|dLBZ!v?yp=eUxJlsc#QeI{o;rbd5f>8A5O_B+zs9|% zT;NXPa^gaP_YqeSX9~QZcm}ae-~+@y;xvH|5!VqL0v{%BB;L<?at>#_uZeh<z?}KM zrNkWq(*pRG5w{CWOW<2hyh>nN1m6nc6#~;T_*N1x6POmlw~Dw)U|I^_8sa*EX)%25 z#4`k@<?wAFE*F><#MePwC@?LFZ!2-8z_cj7UBotlX<2-o#AyQ4!ua+R8v@hP_zn^8 z|Azgi#qnW9LjMJ(<?*EwcL+=i<Vz!N7nqjFH;j0dz_ds{8}SN(X_<T@h?fCF&&0?6 z@ii^~T)caKd1Zx<-_!B&%Q-*2+{f<-bl!xF2bRTsq+M&eWFM+C8pgJsxE~j5X5ef8 zBd6ofi--EORdp4-8@9DQ#Re9?l+q%ybPc_ebTfa!yY>HGZ)FuJi`QAMy87jouiDmP zH4t`l@xEaX)*x$hdt?abqk)S7uUOKV*444`<z;Q9=@ziub?sqrbH(M(pDx<RaclY8 z;)^Yjj8}@&Ek1lh0wrdkfzGoaU8eC{NiCN2&T*J+&#~G)Z)2TpYu*UV^jD#Utwol_ z?Uq)f?M{(enV#B$l5NfFWS$!l$4aae=v7I(S(Q&^d660|-!AivN4)p+mN#0~$7-o2 zZGX$V<Y7Go-qx&YoCr!=^WS9I^Hu!=S-(h)*1uNf`8V?HJ=LsV)zU}SzvorB7qRO0 zB5PZ-C@8wPNmlX1e_~arre8o_vbe`3s!&;%3T_pernRlo#;bzg!NuVqe2Rf-YwnP3 z9Yh68JuFhU%G7^~YDFHwnMxmkfoyA5Xz$nHbC}sxmG(B$$mia)SQ{@hpA?B|_i>qT z&n&WKNzaEibRNURms_@S9L7KhEmyb8ffrotlHA>X2=$%V@>a_$EyL40o0;)Y%jbJe zLrrbXA``N=I@(#y?~qp}0u_LN#fhrz9Ha@!?tm`{PKJ?_CRQprG2nx>T)nQXS*VcU zWVhsGKpZFE00~ZRM@?<bB9q``9Xa_DU&E4#K#xPhiW61ay+{+BI8pwcPj=kD3e7O4 zmNQ-<<K;45hM0pl{8J8TClDJ-hZLE5zDV6KQwJec<`G>`>3NhnO3lSgY)xDA(r8TM z)?-jUjW3WuJTRB_+(}+qpzW;rBBY>9*XwP~s+#vG38@Pz?QfKa)CI-d)66GzL8Uz; z^X<6`UODSG`*TQgUWu__s@WfUD&&{?(bg>Ti2kVbnH)LLAHh6H-k(XLKa{EJPkF3A zs+z%Ce^lBx2V?zF%zcP7)gP7ihRnC86}ydCe^MpqOHg&JKlH@NZxE?%%_5KJk4pa) zg#-Fq(C*~@sSy34Zm9lT5$lht=2G%0$4908MC*@Y?j)^0Dy^^T&vU7Ie}+lUSA!hu zk4>gNfYi2Tkw^4LrLSjY=uf{M`*Wq}&vLClSI7FJstJ7+>yJvi3TaB;6mt`pPwJaW z%a-}}?8WXdHa;UH=et3U^(RxNzJ}DcW|2qqN2Pzo%Fv%Y%-Q7Q;}ZQ@ruD}Y>yN7D z*I&l^qtaTDruw6po5y^zKPt^5^X)km{b_5yRC2sfa4b-<Id8q>q&Z*GM~e|sXbdYs zKdLbY*^e!7Q(A{W4j3mPrBoOOA+0*e{pMIz|3p(tyHtgLr(ihyC&pO%57OT0d}r5= z`<Fr-ZOs!TE3cz*iQ-r-l=NS+e%Wi$Au#%wtbflr*z317tE_VmnEkj1bC3Plz<$u4 ziGH-*$%fTDQ9XejQ^ikWOLCs5w7#;Xy<=h=Ux>7}<|&fD;|2c$N&c@BqMoSrPw}t% z{IMD}FE*jAj&*H!irkgy77jyyvTw9FHnlaYN*1DkVojy_WF=%xF>{At<_g&~<<ZtG zP@HJXB|VR|O3|tsDrBvDUd6FWTeHe~>KCSHFNGOsadgB8@1qUHnacbc$AQCtifkR7 z`Wk7kIJQus_DsXcjp;flPPfYvdn(at>*OuZKGwR;*RtW@M=yRouw}!Ght6nu<2)?M z{`NGusV%Rb_fktob9+O<;_J2;k?dB-7W#?G<J`61oWqMvUF}@)nTeGL+K#TYIyzbn z>pHG?sVx4|i(j9i^n`U$OOOYn(e*mUhokZs)^K?_2q!kJj!jH`1C8ytdTU#A4Rot@ z>LwHvV|XeW#R0^QX6oEKpS-j2-97K^dSLKf?MvDtgUr+&j;%XiMA>&bcT7(oiZy<E z`atm)9&eQ#Aq#jl>oz^%g))8m=S&y5ghf*6cR;c`uHMzwtP;ZDPYsVsb0JOH7S+%s z=98+R(niUAds-o<9gbaObB8taOG_$Q1^b97b-qBM^PW$IAed^+N0zuApJ|p>Wu$+0 z6N?B(XYarr)<V`R5bF?r4@+4mB7Nt7Oa~}=?dMAgVK^!=pLX2O^&^&6D$;0!_YO08 zGjphGQfTvq=iykvQQd)_GfB7|j+IP$6-jNNQru?BA&Li-e9mE58*Y|chajbEv&^v` z9fY&W{f1fONhAo#AIGXt_F*lDNVr|)I7PRWR4!)!?I*EOQfa?L8tdG1GTLA4*pfnh zez|4+i{G7LYB@LAt!Yx5H$MYp6uh<YxPrHCuoS#pZ+(VSa>>TQ&CQ3+PTz@nqvlc# zDsP>-MQL^O$7YrX<<xS3OgX8k<&|<XMY)PB--ZswM7f4ZHz29avC0&sNO=q?osA+{ z;5PPC3Qi;}l%3oQ!HH_;XgNf}G?n8-vAdLeE3`HiJgFI=Hfkd|bkaK&o;~uNtmRW_ zRxRH;IR>ABO|^+sX*-a{`q;ZgkWtIlHt+t$>x?U}y2f>NoBQy~*eNgGm*46b$WFI9 z#7?*!s}rWbL|>yAW-ha!l6^qiC-2&R-PRr0=9CHA{vD3}8(9<f+>7^3UwmC^%HZbH zQOTm~2J%E<aP#p1l(q3AD|a(%I4^J83;A2ppW*c*kkJ~xt^t`+XK<5_%%rf%XV{3f z>T^upT{})|M=Gg&{UM6;9UGmaSe$Q5(2JA8;tNn58B4)J=d;i;?x1#9%S~f}c~@m# zrJ?M=PFN|{v6h;KM&6R@l*cjL`p^uvZAlNXs!+758D<3+2y>Z^VzwqI?hdO=A12a8 z^-O2=J7o0-Wc5#fBs4+p9rx@*W7s>A`M*O*Xqa!_XmfnXl#9_FUdnaz8&uKe*hNYT zn1nvbS{~ABxrw62<_W!H`vW=C$VUxg7Di9ettN_cI9QSOJ&2Z+u8RhiAwekT`I46x zkg$et@nL0Z*jNqO=>rC&;Wjz{r7TY3muP13y6MzHQw4{i3JaA5jF7#~$8ba-lnRO` z@~DhbEx$nZIP-zlAz?=$@<gwPYS#l%9Mqi~MXsR&4c)nM=ZgXl7kD^vv`lNEMY~Ql zuW|My+ii`uOJy&FtXiiI6sz*oRAr6&q6b)G+Z?-FV~M_KqG?}*NarJQNqfU+De~cf zEibn@_G;zchcr-G%XBXMVw!?Z8GZ`_7VZ8#<1+SI#gQ2BwXZl1AOQn=fJzcAvyWv? zL*6#WmbR(;!0W))sk;(5KO1y3MoaD3)#eaozU2aCuE1u(J{lZpR4&*P)fe2#FmOuN z<~;-hbJaCsyda{XtyhzY)~RW1i+$$-31+Y8)z;Rl4{v<QBJYdHyeiYkgiLaTThUCK zHZEz*&wiqV78DU`6Kh%3#@lohq-s*th7|}$X9X_}zN1~!F>_Rf%e3a-^vlw?RHe3* zTTHD({psjgn?tltVWag|Ap+6Yp{95b=w1E`<&)4pb;#D*AVvx018_21UY1>P(&H3T zOT*xE<keud<fL~_29>q#kILJWB6z&oyb?OA+w~6g<HgQ1qWnu~W6YF@B5Z%_;&nr> zl69~%qT5p5<hPsy+CZtX%g?~56s(8rmb4FEvexWB(`!)`4rnK|I)<|LHpkmcox%pz z%c4U8u$@5iiM+xzibw?+n%;Q^lq`vMkCJ7#qS;oIp`bYe18IWLI>oPg^A++yM0>u0 zy>4@CH8ZDb1A9N3MTS<#)|S3|-b5~pGTFHi1<9r4{C%d~$p+(DTnH1zUS~nESQ4&# zR4hZV|ASJQf{7GLt`E~A%aq_#f@l{a@gs~WGD@;PS}-Pg`GC9(5d5`?Z9wO0Q4N#- z!Hn@Nz5~|@Qmi$dtj^yeZ*#jCtjI_@VMml7-{D93J;(k{#ek=^P9EA)3VFkJK*&*b z1hY5GyuUzRVfQ#Nx~%PHZ(1F@A#1;onC#Lg<UJ4NDRLqIe%(95dn1&!P04ywOxBy2 z7L~Qqkf%iIk$boWdDz2upg3*jzNsy$lysXlE%=mdvU>?4o{Ixgh<KLdeH_Xw5uZ<y z9Z|$ilYzZgVTXo|ThH=6UuMi|bL{V~$g&@;j{Plt_YNi}(#$@Dj?w3nCK-l7S-nEJ zvcX=!0GQTP+RVMBvM&!{%+TA`$-^+}rs>AaY{ud2;{7S;^N*O%fkjn}+HKTIN7~*& z61SOn@pWm&;5%U+8h$Tl6zx9_bP_BL+JSkMc1Cc8YvG=bAOXUkzzT&eIY<ZuQ*ZXi zbP<rRwNS49%CfTghVgRB4a6?>V#m-(9yvi3_jA_@dppdmR=8Y`B8@|g^;Sp`J-!{s zw~4lBspAV-PS(Ae8tga#)Hc!5I`zQdryO5EY>*!3Ow>0xbvGB2f#9>P96O?cE#;8k z0nrky`WIqpUA)eQ_0zQCQHyJu+cC`G*eBZ@Ezth;j)N(Cr?+BthH{u|^aoOT<?8eX zd%GFA6O{LvLtR}<9S57PJ7`3XlSN-e(Y?n?z#E%rZ}D~e<RIQnB@=_NiFNX>nJ9e& z+zM7Mj~U8Ldw@hOb1Pbd981<awxzUAh92*0b8N>vhTTu2s1i&@9_q&IY;$Zw9aw3% z3n@W)PD15-9*oV|vuQ6(rEhZxYyBZ_Fu|M%9-jtC==X~_W&|hW*E>Fla<aQ;P98-K zf|CzSPIfC!+9<12);spXeC#!8n-uL@iZ0pQTih<aRnx09m%tDHNhVWU`ap@L-wS0n z^;*w&Y+}B}`_%j>V1u#w@c~R`tR_r%(^poNYE#UdE9NBQ9blox2imIA+Mwhy50;0m zLJCL}V}&^)9Hgs2KvwF*k8lZnj#`|T%|ix3`<hjAb1T;TiXEH?Li6SPu1DF<=ibA} z;b5ZkH@uKu=K7soUEp5YCoz4Zt32;=C-^K}cW~gr6>IQQ_?x=4Tv|9?lCP721c?jB zQWzcj`^u-_-UPe}wm(H|S{+I)4ZcvzPB|NgAdeV)Oyohl-b%jG%&AN}WwuT+afe{m z!Zg$L$(eK6cFJX7Yq{8t)7?Z#v9rFQwSg`Y_iA)?I^V-t8+WsCv+VF{g{vsW<F_zj zI0<;W1&b~6)80OZyj**~X6y6N^Tq2f5WBA!U+gV=X%4Qas9JP3j&l!^E_MtYd|UhA zr`oa5(X!*orH8jsm;+lnQrc5GUfk7}&2*qAN&z*`nk3%SH~X*o@xq^+BnpwMHUINJ z|8tBer{xg~!KFBdj!Dnd+n*#!>zgF)Q;5zi3z;t37eJ5O%yZPw-sM)No$UP^m+Kw7 zjP;IxBD|0A0doC~q~{<z41f=t{?1>3*S4vDquDa&-CE|D4`r2^9}t%vlx!Uev`t;B z>?7X?E#}fDQ{13+DB6W!TOJIVZC-(DI2HDMh&9hF#rCAym5>Y5(tQHv)RtYCz8y#q zQx3;d9b|43bC|ah%mnsF)a^&zZOspV2Z(&dkbtazK60=J)H%`3Xr>o1=&>Td1Rv0a z0B6fD&|64JJZSL@s$!Zti+cz8h#%AP+{xx*tH<SLla`l|h1|?dsNz?sV$WDueb%!{ zt7me2J%eODAIW+$66z^qJyTH+RRC4RI<v*hIq)-9E++;XgAA0uoQW*b(;M8p8hQl! zz%OKe`Zo#rKYNGy8-~C!*k$@p30wEUOM<8cdY8hI-)Y^Ni=s?ZBK#Euuj3nuE$MiU z(de9pytE3G+XG0)m1B3TV^@I;6_KYzJGLSNpHs+3%dijFkHGK=+sdUIRXq1LwM(qx zv>!M`enV0gp>m~4B5NI6qnhpoA5wWnle47C^K&rL=6H{7%z`9&iG1wU@UM8LB|0z? z#n~(<OaW1@73tJgGtYleS}DCK)QJS8^rHTg$f-QF6$$+$r+bS{!C3jcAUXLvii)ur z+`Q*oUAfaCDB1Myg{XvH@{m)CAUY|gy(L>eH74}&NNaNl>o7^<RwP`caVw0^d0GyU zFihq6k+f9(k<jKnAdXpZfMzbS>&XAd9DKQG{t_c+TK9Zx8c4m6Tg9o6<ZfkZeu5Mj z6pGlffA2XmUoG<$T3XW<zmyieaUg2QRW<NrLo5+I<rgMNtmLwSuQ5;H@UIfqALq+5 z`El;8;@bZMQo6VQG4q-24~=jCNmWC%{o|9hzW~{iw11&09qT_*y0_1>3)y~ZeEYeo zhG_dd)$HDXXlIi4A4hUR|B=$Y{g0W?>_7b&^Tv&k|C6eQX#2+}YkvW<McWtNU7JJf znZ6Z@crbpZiJD)Awm`w#rtViGdooUq=~TcA6gj9*Y8HuvCM1Y?HbmC`I1<dgi_}Tc z=G`(IHSpUli8^^Nrck_0?m|j;nY_b%rc8##%j7{-LsTZ;g7a=;G85UNGKuNrsj75T zC%C)l-u{)$XSP2)zWp;*4bk?mOV<9sVK%h(hw(N^!+3d%q~Q{n?+rYFzxT}~w?asv z`2Gm-(L&1mC>vSevdz-Vf+tCueu#b>96>n;)%o>LVCfQlC*DTwInDHQO)u*O-cO2I z;RP{M_yU?Ton0~Z%2py(*`t&<gcC1aS@_LxS)=>K%GZ%m_hI`Y_uf;)Iwg$y36P7r zEAlkDY&zG%hc!3K$?HAbkx=i?66$>pB31Rey46ejiLSf?q2YIs=!wVnAd~}LYY%RP zLtlD6SH3guTv@C$cY+IOz1;1|Rua#(#gNw~8L!`KX(RR)%BA+S&(V|E*uZDRgcd&a zA>8vfk0vzwY*fPimt)}K>juir{ga7Wz<aRle<AFiyrLUuriq7F1T`(1HdJ0wTe+wT z?KAy%xyu|V_LT?Od%1+H0^?%39#k93r6?(F%3nd09*NYcNsxEq{!$p4dnf?8_}zhV zkK4DD2xioHCh1R`+;H(2N?q%B8H6N1`xhq^YeAp%k4R~YoNL6MSB64%PC=`SUY<Xa z7G-#6*WbbUbHm_bw=L%1^Imzet50)RB;Ed%d3>QxRL!1t-SGr-Enc5e@Lm1R#p{NP z4O{E*bTn;hqUFI#Jp8=)P*>!H=Xf=9+2XyHu8!8eNNRw8z(=T~c{{-izjLzh0(gm~ zo4dNM#P>_EP~cV!ESUVv4>j0=UT=*p60(KsydM9|`L^m%#AkEc@XloMQfAvsyerv; zHwz$5y(Gd+3paR>F>JfYHWsh##JjS?wgI*CRk;kjAjqtscD^c?Q6pb7ZL4nxR)&zZ zCIeT7D8Im4!<Xxs0F_;6&z$J3@XzrEU5d`vrx|!j65qljb4!GG@)(ZB3ybl`<By1U z^qB8PHXIJ78kLoZ#u`=lTV-5{zgHUgnk?Q3%~wGil|Fya>6azMdzTHHc>NMOF8T%$ z8Dk-(h&x~l&9oUdFJI1N^P&6IUT;v(#uxF0X4SiEYi;&o+r_p<1ibb)7%|LO^odAj z$Wj^R+eC6!v5`H|>zX_tZ)MJiG)6F395AB^(G8g4-{HKb>jr=(0P~|GaSf`Wh*};d zH*q#A8WVV`+R-21^TLyZD+c5HWcV&w9lj^V@8uaZ0xTe=J9LKMG7z8rmnM(X3(fQ^ z&A9LvadFd<jIpG9|NeIoOZ=NVPN^MdY$%&m`{Eyp9(wTQ2k&_K>!QQI|7`Q*>heo( z{rJfh2iLC}Q}X!Qdp>`A!I<Me8oB81yRN<DU)P+JW-GYn;;q-c`(?-H&tEibM)~t! ztvchPCSQy9t%9b4lG1N(yXpR2`Kr9CkL}}IW<NOpsvXsC!TlA#8UNldyS1z6Du1$k z3Gy-w$W@iQ^R06`kIguK%#r%ptvvew<-3nq-hH_L=<+|?zP;^tZx*Tc+oZfz{UTkC zqdtCr7<<e1ZL5DZ?!PV9_R=T6fBBK2j<Y_oxBm8=hGNxzPw9W!d~xH?-D?Uj`uMZ7 zYYGpJRrtd6tKQvx_e*2DF1sP?+pRB+RrD9<%=Q=j<BGAXWVt=rX=P;xhZQP5Rlig} zw|sVU=iwjyZuG5^ud@^QI_uDkCm*;VuPCm+n*5bsjw~9!_x-$6E*<^Bbyo~n)QoQg z@)OfDHBKVAX<XqDsmH&P<9E{KA3pGS-_4^RE!zD0)%98Uh9N((@`~S|4gKSUe=b;Z zX<|QyjV!MG=F9<wW_+@Y|8wb1$MYLU|K;q$^T)k$&W-<kuXx;Db)Ao`tuNfWC$}X3 zvo8x<@A&x_RhhS4^7^CIWxpLU>ynLs36Fo{^u3oHs#uVJ*5QBS8;>R*w?FrG^(oVK zUHa#(=bkvSqOEZ46LYsUPj4%n_v{zzzBsk5@bCcp{o4zgFa2BUd4Yp{9~*u6!uqrK zJ$K$Xi}RlS)~WL^`Ay#se!S&_>n<I0SMk@OjY}_?bm4!Vu5S7IlE@V|y>h@87<1)4 zTk`vU`lpLN{PUi{^EZ#VxN*T{KmX|Qw=bE=&yZ^i)i|G*_e}Nrr{@%PW-nei`?H6~ zsNZ~Ab9(dMpRT#poL_~%XgB%3vORu24J|a&hf4a4!m{$pDN{-f<3jP<1mO$SLWOtT z$9^o}LgT_CeGFS%ZgXukY=Mv)uXUbf!$MMDAF98=77Sq-(5e*~4ckn<)7Ci0h`sGU z+MG?K*~Cu(*c1z9M%rimZ$gJ|fI<KNYF(hbQ%>!BaVFLcFWO%EYFAe|)`_(^h**la zAC4-Ta2R$ePTC$ud==ty<n<s9Aod|%ig+REh?gR6M|=n3?T9xZK7@D+Vm>dh3o(Cd z@Brd=#KYj`?nYdU7+(W7W+NVsxDhd*f4mJbpS@gzxDN3q#J3>+9PvuT`R{dgtwVe> z;@ybvKy2XN{htvJN4ynrKH^Tqm5BfNL04A<@w14RAzqL8QN-hLTW<s6a>Tn4UxgT7 zYd5Y(JRI?ji1QKOhPV>(G&qtG#OHt5)pZBrJjAOII}vX}d=27G#6LsKU#Yquu?_J` z#G?^Efw&6s^N1seXLfdV@i&idLA)CAO2i$A|B84Y;ts^AI9EP^I34l1cy6Z{@kNMd zAfADEA>!u|^SAnTAzq94bUdiF74f-<4<Nn>aT><45b*_wD-lmbd=26n#5IVwBW^@| z2=Oh5htiiooQZff;!?!-VIx_GcsL$pyoDI?N@8qP))6D#g*X#Fy#*T$?h$N=GvWN- zf^P*aLA(<23dHSbXEox{pl?82g?KyS2;!gM`R(P17b9NFbi}(5e}woj;xs%}J`CT= zNk^Q6cpTz##N~(sh-(lpMcjz^VZ^T?Zby7N9<|?&I34jJ#LE#62S1M@&PTixaV6r< z5!WHM97KM^LlGZB%rA8Dld_<^!bn+|+E4tN+o=OBmZd4A*%10-4q@mFW9*>Ram{^> zgL;Z`Tns+Nc_eWntqMt7L4(`YjTS*t2WSv$HyZb9?LE*Q1#J~*_}YB8avPAep@(w2 zL0glAhUeAA2YoLUo$qh{l%3goeAgb-bcD?H7!&xQmNiyuYVp7nqa@ARr@<&5+_#~w z-^8@OeHtQZ{rWUCruFaBun_qvw<_4c3*PMNx(wHevWb9}1AY=^V?)+E5H0|X4<A~7 z-QSw3WPash3m*q!-8|#J0Aq3iuA_BNSlpL1e&KgCXy+uM?E`H*Xi$?_Ib5<Fzqp$Z z8Xny5MoR}RuZKDdK`R7J?6(!Kvd-zU&NyBqtr4_wDCf}ILmxOMcYqcEtwN_w)cOzl z1E6ENb>p=Iv|B(sEeUNOXt#oPu1=G3ONDHggN9*_tCRfjjGuM#HyESx*7zw#dKT{R zB-7^^pq<qNZ6RnwK%>q?>y`X02dy5oWMj1&w3;63>;TOP+F5#?xaJt|18o{;$>uVD z0C0SgIz|6M%L9#bEXt?sKk3Dw?*#oH(62-~+pPipR;*je<kJM&^`Ipiulqo|p@(v7 zKnwOzZVPB3(2~huKWGg-)Hx7_tfmLr2+*8KXg2h3BI+Cm+Dcqo)@XBfQcP!kD8C+g zlIh+O&^kan73HJy!zE?6Lei4SXDw*2C*g_Zwu06U+S$0)>*RPH0Br+kcpXWMrs?xg z_zmxZo@k#~ZzgD)L5u4L(rpOksCPSPSL52UU6Vt}V(X2Tn}mI;L-`-=N<KEXfHoX7 zu}Mj=B`ZN216nfOT?g9vpe5sX7iedJ2G(MI#U=Z47_@;slp6-0#{yb18<PXt=jdl* zoi@}--|YH-bafTL)>t|<8I;FlK!5A@*Rf6`zdknPX&Go+;Tt9E%cGz*^)N0QKx^!Q zwi~p%9%wvA@g<?zz>5v#W`On~^xKl28e5B|#QIW*^6Bs;n~*=TFVjI=+5@c-w0%9$ z?f|VEeraNzHZ-vsJgo-pc3fMY(0J13+7|MJXN(ieQ@;DjQxCL(@GI-!$0n2g2+({< zXg1V25p~W0jfbq34c&QSo&fT+B2N|bY|-+_xxqZkkmm)JXLt8J<bM_TIDkBG-MYzp z6KMYi?HrvZ*ON}r_Jfveowh*9J3&ikV{D*(1X{9nXf$ZsK}#mDD$us{K#PF33ACv? zpSTwK4_ZIye|Zwxqo7>`TC#b^zNdp0K{?BD17f<W>BT;jkMuCN_=BVY(Bf>O<S8Ar zpZ8F%5VQrLCGy0+Ob6{g=y+w4zBGci1++vPAjSvzQx*u1G)4KMtdDtZN!nt1KEhDY zrXh~b!!mQ8GO4R8<^8@R@dB>-*7g+RnH1}aKE^vK;QH=9)(`p^+xu9bJk7YNul2sZ z#+`k6y+Qol+1L7n#rUePwbf!QvGDp)3#o64zyGvY7h8?P;(Ddk`WLJ5l$F%?#NWME zYg<3##(qrsb3f}V{fu^TeQ$rJ{Hee7g;Zl>|1o!^Lh*UMCY9-X#ot8(tiL$MxO>2u z4+j{##P!l+NPSNHef=2gd&e4kjv4dxvBs)ndA;>mQd<Y|?}G!aj~{0|GjPmJ#~Dq> z@p|QPOy4H{es-J{lc#0S7`$DjV-T+o4I=gB!PYO+@bk_&ccvM)rt$idH0zT?jJuD| zd+&H->+!t);}C1d3C6A=)~|;cUx@3CCs@}HH4dC$y<@1+Hk8-DKJn%hRJZp>)<vfp zi%+s{JIT2BB<m9=8IOwVEhjPk&S9HV5Pf;Fb@M64O{WZBb&B!(Q+WNiQ%L^$6mqoo z)VcWEKHS=Rnz4EKF?)s^?~Cgtr;+-{)2usAGu}GQ+I+h4-Dw!e2TsQy4D*S@e+@tE ztrY8LDaJ<};hXx9V4Rt9(jZ)%XN{zcv%31Wq#Sd5it!LLAp3PvM^ZDbqr^Init;X; zZy4KBtS|O4K4G?J`&e)4Yiwa{O?|DuvKaUD#pPf7THo$#bch7#hOFo2lp(l0&5CI` z&H9@Z;|42M?t4?L-=!F@rdYS97;lQ}gDFh^bszrypFXSmZb~)&Wwl~59kS-)`q_Sc zHuf`K>Sx9EhyBjL-$(jezvyqg(BJxcf8$kg{j1a*q<@^sDh{Pu9~@xZGl0#m9RLR2 z=O1+C_z~7d<09)%jmr?<%FZAWKf|1Z1aUoZ7~=D+nT8`d?m?(sDFaVMfsxiLjmyqw z*CFTAQcj?pDdSsGtiMb#9_(4hU0vJoTgEp0{;?@H-ZM8p7!!_*2sOCibx6kivKH}T z8S{G)#Or1JhKzU0c#n(^$k+<s1eeAT8UIekZ^;<;K+r#t@wYNQ4*Lg|`<aY+9!6XX z|A6szGCm;VREUA;r^q-<##3bMkul%1$b5XCB4fUnkul$&h@%nl`{-PM_&e%fe%BSc zWs&*&-F3gT;JgCsX0!Qw&+nXL7%;2iwh%0x_<%(t2XF2|7mQ5&p?$?`3Ju&!GL}DM zN#*rY#6uACb`sYz4#sy~U7-zYEU8_H!KQ&@q0nVdo8tudUwwSIBdQ7vT5Cx?9U%i@ z3_=;gH3&h3>k)26cmUy#2+t$DhOiUi6NK*&2I1SKrz2z_j6o<vxCSAJa6Q7U2oE6q z5#f1+*AR9he1h;D!XTWJo{o@#Fb1Iv;TnV>!u1HZB0PZbM}+4QUPIW4@Cm|q2!rsJ zs?!lN5XK;sAzXtHM7SQ|R)hx-{)q5A!fOaS5k5iq4q;F`$|GbTj6o<vxCTMU2d_`P z$Tng^+0^rFIrw$fOj}lFW<h3dW}a<C1@3cBbVtnekvZoXeAntxXSSO)kIKr8uh~8_ zt9#A=KgoP#X=H@w*GAmcMh1Qy(u{rT3co7tt*<jOf}w~v1FxbT8F9~&2tV|dQO)1? z8R@SvGKg^B*k@$a%nzbIGmg}ol)2_lY{vl?X!YIz^P3lSfrybIes>}RKNX2<@f&>^ z^&#=QkQrW|YbNeLdtJU7lvD^+;s+|djUI1Z#0BPY+dG3lXdMOvW;EM_>}bSY>-Qjc z2o>W`hG@;K89z9QLZMpz3a^@L-zVmto3@7bAqRFP2SE6<vM_1<oQ;63hsZY(`6hy3 z16CiHJn^CbsPg%c&xd^JlEz5oE5<dBWJ7G1_Usnqvw%oj!w>T@lp)ZjF(2*RGGy3= z1lm0m5qW60E=QoPV?Nr&Iv6(EL-k=^ey%~FO=CXVv&-;q(S;zYcA3sF2Z1(_*R+Eh z;9$1nn(Z<#dE;jR0$guVhIVvmfoVroeay$ny;#dfySp6uXpfnXHeThs1=r-2?QD1s z?Ou$0xK<z4?oCKg`3@o9A>^x5M9_H2aGRE|4tA!l5POf2NR%JSj-TZSsy?nWql2cu zEr|G|s(wJ`qg<+xuL}8!qt@N*E6Y5h<+~5h5#Ki#`o#f~A64IKV7ANpI{srx?SRkK z$pM!imG60Am5;~JYnpLRnxwuCnUADWtUIMG=6J~JqWySV%g6T~PDehi8_|3_wS0WR zVGZ&ro>-^4{s7k;Gq%h1&WC&(z%bMJVLsNs4}tw)zD%sknXR!pkR<Xk?`H@qAJ^|3 z<l~vR`Y;c}*9a<K6<()U1qEWd`b4jhiS4spej8vq@&&Z|Y_fnAX<D+ud>3m8>XLo~ z>r$V)+6;r=wNW2+$(UoSK25h{504Go=}2IH@<$yUfS}ixxymrk%tNC5sQMJvW7;iZ zzKVh$+9$?+Wy`cfCN^83y)m&V9oiKWo3obo!^Gz3DZi($uxzG1FzF(j@SjY~UYLHB ziQ&NE!+C6CQ?WT;O+1j1!Fg$7ampb4JQIsk4l&P6EKcVQ&L<O_vzPP6#2gje>%sgm zF}*r-9+=p)?9_J?o4al5wTVryQTc0qjUO>t=&4Dk*J%2aCg!djb^!Wh;*$lYewkR@ zO%Z7(J~d9aEXHSDlmO=`=_y!wSVHla0&jH$u;MSpIMK)jULobx$2iWI1T5ZDiZ}(O zry=GWPdT#0KE%e0rhJ%R*CV#i#j;9ke-ZF%DZdnW%D({C+rJk$4R#L^wCgzIPrQ&H z*7p}&^R+5$hwalI@dYXhzX#k7%nC?9jQBSlmQ<dXvc3zDk*`=W>4S|izytQ<-oD5` z*qDfmV>HY!E)N7Y`#ach1E*=2-((r8;hzBy(=Z=79<JfLfo&T89dNpadGBq6hTDNJ z(C{YUObx#aoTK3nf%7%YJy{{J>M!?Z6E%DyaJh!h0j|=@j{=^qVSX8DhKBjNnHmjW z0qoPT8#tihAaI?A=L1JHyac#W!?yu1)bRblO~B|D#&@vsIPfA3zXbeq4Rg=ASi^gP zmjJ8r{0?}jhL42_yG6sN125BX2Jjsk9s|4_SjqQt;1wDU0k70Bp8;8=;T6DZG|bmb zv}<@X@CFU<2JX=CLEx<#J_a4zrC}Rzr-t)^_iLEXupH8`2iW*5)}NmNr)v0@z-bzO z1bCQ+dEd;Y;f=r}H2fZLriQ-+&e!lUSRx8Fd<O7D4QBwCYxol2Dh=~pV>2{79oVPg z0C1g#e+JyB;a1=#4c`mARKvdqUZ&wSz{@rK3h)XIzXiNf!~X_erQxrD*J#*^C9z$@ z!+<wv*bdyG;bP#e8lDEcOT)8)J2m_>;QboD8TgQf?*TS6`}k+zR1N<XI8DPlfrn}M zAh1otgRxl|q2V)uGc`O4IA6n)feSS}9eARKX9JgOcs}q94c`px)9?epbsAm^+^FHV zfSWXY0C=f}hhSna({Lv6at)UOuh8%e;FTJ#174-!8-Ukn_!q$K8vZr#1`R(B+@awY zfVXOx_o8-b_+P-C8vX`&zlM*23w21t!{}0d5u5Kn22Ry*K5&|bCjk%Bum{+t;V|$B z4f8F$nHqioIA6n00~czz19+l_d0(kq!`}i|X?O@+z!@4o8`!5|p8eHn*a6(A;j4k0 zG~58ZRKu;n%QXBT@Nx}50lY%Pe+6Et;kSWTY4~H{H5%p(;C2lUfeX7q!)F6`Xt)4) ztA-~5@6s^u33Y1tTHyT}z6tn{hVKV94#dX)Pr#`fei=AT!*2i&)9{DDHVt#LJwn67 zfHO7xW8i!Z7XTM(xD<GzhNlCUYxr8=Dh)RS&(QF_z&;JH0<P2WM&L#bZv$@9@K?Y~ zH9QbohGiN)A9%Tj#{;j>um^ahhUWsW((sMIYc%{T;C2oF0eFLk{|elp;SYhgYWNsz zVRmWwLf}pfPXykt;Tgb(G`tYlz&EgACJg#hcLJws_|L#;8r}pvOv4`n+cbO(wott1 zz<Y-Dx6cJ00Xmj>QU4X@77K`b)6cr-zw=vvrayOzLI0j#_fweu`=&jX)JnvxkN)~u z*q<r?Zz3+v)-e6`bYRsU{q+$VroVoHhUu?oYMB0dj)v*4=WCe$`e+T)UoX^fJ8-dv z>90@JF#YvX4bxvQ*D(F{DquAp^w(!-nEtv?!}QndwDR=VBO0c^-l$>v>kBnZf4xb= z^w*bYnEv`w4bxx0MZ@&hmuZ;(`W+gcjq&FFCtV+wYxG+|KSu5M`W$Ec7Fg|9jXozB zPbJ`uz&x*kIv6Q^h8gbxtNnCJpC1{W!2B@|lRnh=l33QC(r2*I2V4654=t%p$S?%@ zlmR?#mnHQq;In{N0$09iN#)L+cn|QDofh7{7kCDa0QlP=d|s2~y}%>>X-VbIpZI$7 zhyaKfH#LE8+-pfa2iL5xD}g@TJVIc0k^drMxqm$#``1YcbZ-J~1wM%R#tyUo+krQA z#`^acV0HdM{?;Xw-<E(s2hKqM6~DRKVKuYP08Z2RyB3(|nIM7xLB>2_{*DUVS|N{{ zfWO{hN#(O2<ZmPJ>MfSk3nU%}LGODvCg1CT>)wdryEqV<d{+T~wat>sX8~E?X*g1u z0e(2&DKERaF$YY3Cj*~^@l)gbQUd*5;A*tL5P4bOH8_${=PhZ*$;KSuVH!Wz0Z*6b zFKMt}Ex_*gV*K*PUn%TWlW6}q<Hx|aN_nvU1;Fbe52a5#%_AS=Hu?QB0o%<ZpeTKO z0-g(e5d5qD?FDXxzQHVu{`?q6Mg^Mw-UrOLtSWh}0^S772H2k`fzN_`6@NptBPJX2 zoC&-T^1v{P`o{w|Vt!AN_}9QZpSlH@<MEVMp5w6_c(gX2-vaLgzpB2=dBa40n14Y6 zz8Clvv{wZh<+&O7QIw~e5T8IpBR?!(4m=U^x8s`h*MQgjJ2pOiK&BYur}Xn;Tug(! zR?7U30PjToSawDE4}dcupA8awanM?$>Bsdvc$FWve_H~62)L>P_eDVCe0&*rDCQf* zME-V=j{B&rg5^I>D1Y)9q?jMluK|v{Wl4oO66K#upuYyZT&wR>;Ai)l^ONn3v!Nic zlK=I<x_&p4E<bGVuQr8D{yt2=UBLTcPvO=Iex{un%g?Mq;H?@zw`+9P{}6CH<gLcz zJ>X2pHybqa`xP+ni<tU5*!VV~{Hd4->bp^#6S@2&Kdj#Yyan>Ez%}c86u2Gzk4Sv} z+0Ys2my*XEVD(}7>wt$M|5RMF{1?DB%%2os>N|h)PkosG9N?jvzK=<u-<N<t0mg$P z!k)9ft_1qYFwFhHztYFcflq-xQenu?GGP9M)ihx0*9u_uAwR2uM{E825;$GcU;Z{B z-*n6}%zp;3`Ve0X{1)a5li6MkE{19Lqyc!1w!Z#`8RduN*8=BbeYzdj<nQwY`oJH< zAAC-=f2=VC_*v+a>fa(@^<jN0fL{VXYP>cAU#_huKOUiSoBhcL4j?~=mHgiad<OJQ z*^}3S-$Hv<Nk9I)`0{5b;7foD(SIAValFcbKh^B}D+&3xCE%0KkFS3maE&%TPbAR4 z1YUsgRs9|F6U-0He>GpO1uoX?{fEHnLwTie;guia<AIk!9?Cv8CeWV*9sxSdFy-|v z@YkArE*XjS0C<d~-vg}HPaD?H-vaY(PHd9xeGIHV#6vRRziIx)rNEmZKh>W_z`Nit zP>tFCtHA2R@;egnH^5!+ht&94E>yYA_D=@R*X(~WaE>O=riA>@C*aM%!!>>X0C+a+ zK{=Wsf0dc>?L7oshVh}All~p>t>~{q;%W$f8Rj4FJ(B(ia5?M+pNZ!DcoX<A<bSoy ze@<3>{bPY&g8Y>}v;d!r_KuPHKLXx~`EZ=XqqDJ(LVb##yMdd~pM21%zpH_#YxebJ zVD%Y-i}w=nSHKqVe*$Razh92ZZOZ#x;3eRP-)$j%bprik;LEi6IypDK{3F2Iwf^xp zOVx+`7@qn_IMOh4W`@Ubxr!>MRJck@oRu#8bh)cI4qaBHk!nJ&S%FZsJKz#efVkWZ zjd&|UZCwDbX{gD_Ey%JPB)I%FjYQdo`B+I!Lv8JR)sCaASTr-fV(esxOCdI+u4Ao; z$4x{l>WvW@9xyzu)Sf--a`VBLBA?eY2T9IIX}-(l4~eH-@MGw?MWquR<)!vqZ@}w` z)Q5sDkB`5z&7axMvAbMhzVyL^pT_n_{O*8%fp`Uk3$Lwk)p+rUPJ?(Ag)7tV=P#kl z%5uGRfrc<^D9xQySus^y<;@vCw$v#u3+7CAH%=`=JIs<#qyc|W0ehy%==6o^Bg`qv z=?E{EP&84?)Dk>CgN&nCE_Uq?%4@UX(ok?#v|$bSg9_N~qG7_k9A-5+iN|JkNza-y z&L0Vvd+W!|k9b)kcTQRP6eq9p=1dKSy@Id2Ipy_up2?(_1ZReLl|SbScYP3FX*M&9 zZ#A2j<M65vju|@coFCQ*<MF<Z015W2IVE^JYJ8}Ep1ZzAq-2Yfac<9?Xj0A`CnPvF z?C}Qq)oGC}SELnt;}Wyy6o<+}k@EVG$BQSD1al?fqChAp`ee_ZGp^np^!TP!&xW#y zdNO6w*dR1l3aTzjR~li@tn~*1{#xjNJyjr@ZjYsB^+d>yBV?i%kz_T+(i3Wm5#nmf z?Z!=>NQo<LCdHL2nBxzGx|M~7sN}c`wDh=I^6K36VXr5?Yj{4azI)XLL3glwYQ8%- z3llI@t5=ab(_3Hb*Hdz;p<H@GfslBTR5YEKn$MQICD!>p-4e5caeQPo#KYP+m|a^N zCnMlE*<^X*vwIS<d*ZX(<8o)kXU>Yxn>8mcbM~D0+`#di)WwZZUOnDD;Eofg-PGXj zguEK>EWCCnPV`}@f43U4k)-F!sZB`83fJrOY*egcdnQK-x)ch=)sxpfNiv(9fbma4 z$mD2uPtWyv>mv(vHpH~(o?7dj<?fc66%6Y|?QtlpHV)=^UCt6c_b9yd2NSSbPqSm5 z&T^|IY;GtJiKAxqNXg@@?U_2ue}r1|IGR0|t@B2Dq~^{>AL4W)d$RX(o#(80m^~)} zPe{O(3Aj8S=d+)Dd6M2nPGW*-8xm8q=EUn;9LlPVgE{lukvQ4s1-<hkq4+72>GRgP z>fB+EJD{_Z4NIeAd!{>58w%G!ed0uw=dP_wNRo=3oB+j6LJ+g2dwM<<DN%uh^b!+u z5)-mRzBuvb&U7ctWVs9`O$V3V(uL*kp2{iNEj4$Zw>o}t&5d}&0k<yE+*$61aNO(? zv$k6*W_PfDzFsnKjyEvZA2*CSwV`+wE~xf;u5G}`#EojMr=i{-H-1#RxXP7s#wA8| zF)<Z_N~~a%2M(z>p|61kJWZ|lAv@5}sN+%4q9QmGdYavi#7bYNA&eU12$>*ErV|sd z%vp7D%WZCCUMM~#&*u)zj8Dq0^T(}!`D5pXX1QzXLwbXwJg$I$mM@OG0*@<vZG*es z8^>C{&s83PTQ)~8Znyg|ostl8nI5;Q*|SAbvZ@Qb$x`ik%Idjl-Qlpy6W1GiPF;Nn zQ$8szH|T8;yNo0SF|QjE>dHov-mC1|hD;sjd+Wk5Y{9rhJIsBp3#&w30N$>iJu4Fr z^rK5+-_lU;MW$dlPQqDsPXJxQj>qhxi|e(ymr*`zR(7=aMcC7^rvYzpR>Y^bWVeqC z`RbFDhEWUD`P{C0Par{RLUUZTo~-UP`yAIetZ55s@z%V!R&t${<*rIE?V&q@4#g3( zNGz|1C!a{ko-@H6uJnQA3kBkIDa+mv^z*wCURSNx9mFOt81_dv1t*kF88^1n<tlN( zeb2vqY(>e~aixx|Dm&YR<Bwe>WPlw@jf?Lmoa=Sf&2?mDc2C7OGSus!CcA^~S>AfF zS&ggTo|}-RD9bSyJ7a%deqP2bZ^TvSaYcL$!8sY#jZDVw)a9;+S#e>D6Q6($EO(Ex z`YGebI~|p-O7f0PeDRfKV<(psp@i7YM|alYkjn=H5b$!Z(v2B=L5bQ9dn-brh<M2r z`I}bejO4p2s$CIyekIPZGm=y1_Sct8a)!&xvdhZH7F@3S<#0G7_Dq*k=#o7gDRJhy zkdza^aG|8b=@uytyWKU}Ro4)41w%DnM>z^*s|G}p$i<c?fix;g^3ylL5ktulS#Wr- zbyQ4p=0Ora0PFcKr)#b|&;Yi{TvRZEZ&$TD>@{lL{$NrO7mM*GQ+DiNI3N8pF>Sr> z+7b`a>^!7^l-<~Lmv|hGV!LB*Vc~?*l5s^<Rj%xe9J52&weC6yvzzP_<Qy9%XO_$1 zaO75dXJM_Q!Ug=*u(?zrm@|J+w&BXk$OO^T=yr*$E{{9n@f8-LaePC`WzWcF`-zjJ z%o#4SQ{Lf3R>jakq3c8QAQw!aKp^Dl)}i3M+2umXx+~<UqJS696W^U#?!y~HJz;V$ zRUcAh!`<pDvu9W8r$fjiPH~)(qm)}VCZ9PCpzEdicnzsI9m=xfswVDbZhSeY!^sl{ zO-j>fO2jv=q})B>28pQ%Suw>IX4~~Wb-b=lP0ge1mfao*&GXi~s`*AaN5u~{UF=zy z>GfRl;q*pKp6!{mFH;&KumGbbcq4p~tE42ikUzjQw&#H_KaQ|42e3H6!=n!83Y(ET zp<-<L#QfU%<||ma6wuD$yv%&(YiTG<)pLgZ<^f@^o5wuZm^t40riF<q2#sabGUilH z%o-0a$|rf8cAUA^;J;}fq9&=AY;_JKrRz!$=Gy6SRe%#%5PthFHccnxV?OdMH``vt zQFh^LfiY^Y%Pzq__aMFQ#wI6bZ200}Qy&Yk&<EkZH8$F_vE?&Om|TInEm<*x+s!n? z{F-b15249iO9S4S+5~fk0%3<g(H?)_WhbdsjkA5sL0y}aI4iP1z}krSElXd*oG^}> zA4oUjCQ4NAwWSzLie*T8FW=X?st)$Lc8<c#=2!`zrw4}P%uiM^%#U#YgfkQ^(mQiI zC+uG!Os{KJEj)hHC6Ip85u7|aCO?R+krl&}!A2+2O@W2J#a7_#3LX!d`<@DHl;NM_ zZQ6bh)Q#MAo8vzi+e%&RnX&h>+Vilond)EctqWnlA?J)UjI+L49F&U{CyCco7f9v} z!gYpA=&UY*pP%Wfw7X0n7fxul7_g)Al!Jl?;VEO9M3oC0Bu~AUU)!dwk0o;f5-x7H zk?Hk`qFGlzZZVcP3tXkunGOfixo*q*c0IWbSvlkQ4SKxp9d{>9-w~&BbR?DH{`b2I zoc8}PN^j=^3wgt8n){)sBXU$@S{LnObl`f&SYTrpb(wnH{Ll+fepEqp8<~%}fcu~3 z_G4O2puD@=9`zFBEVbuV`f%1<<PLcF0)HG(MRPcn!bkU&o`BeaV57r{KdQ0<+uTsC zn(z9w%*p{P<(!VFamq~GLdBL+IP_w4osj@Gs@=C@6-dKLF1KuOW6d2~CAPM?F8E(b zX2a3C^{^xmU-CUKw$kVw_Dc!vxn?ysx^Pw!4&lcKaIE2)o0+(2>qS2%x`)M9zWcWQ z$Ucx+)+^o)VXn;IXB(XB3E*B$lARv*oKaIcdv^DMm>Qht$1VK|Zy1NB5L!_YHqqGP zVkPZGf#lb^8(nH0y=*R;n}t;fhp;?y^NQUdtfM~^7M3%Iju90zeh}Tb9VJuDZTQjo zEoxU7^SQH&LlLSJCdSM<mv0_>?+oX;>*~0z?IyE2dDw-cXOB0Kz|p8mX9RcpXQ_GB z`$j)uPs-iI58*NA!jOpvJwo$Mi`&DN#%%p4H`u}{n1yEpf@-W|BiqZCIUoIj8N(zw zLPw$pH@O}@H6wPlF$-<CgVho^Qyh4Ch|$e~l75aeG76b-4DIdNQPB)0W`i)8!f|mp zahwM)#+l=S@uV9QKM}saS;z@Ng~j%lW~c=7YYGx-;i7b_G1~R4E5*SsW`Nr%R*Y_0 zdosR9G;BD-H1|WWh;gT%M=WJ7T0C+v59D|-M9)STsOVysJ=r}=IUYx;iE6&atb+YY zv&zK%o0(G_0m3{aKL9-9*{ivJ<%oG3lX{sc^=|r<l_|H0*oDN7gvgp;q24PF4-*%& zZq?ZBV{2-nN6<-5MrD0)up0epSAoeWPBfzu?{#aI#f>>8imQxnvpD_a(SaPys11g3 zR;FJ&<>l&ngdLbT8~#CzBhSU6C!V7!&4-?90uJ+#w%3Iu?gUUs;!Xf(ku&mNZGfEl zLfMYApdV3=@<?8-ZBSq=ij-{JDJ32$!wtUq=Ga9QJBNOzFr>XM7an7YW78fzZSokk z66U<4Bw&J~s|(DD+Q-qZ#nOA;$mzX~UlWvsoS8!4z24)Vc6`cup#$H?TO8#NyTr<% zR^8a7I$9kWC4B{~u5zm!i1fVjRbxVHyL73ay)G_IL6Q$G^+OyQaMkX%&Z#{Th_jg! z`Vm!+UfM2Zb6Urny*N%|${Cm3vu-Mji0sMk3A&Qsk#tl$%-NBLooO}udbE0=j&P4u zUejWB<bSc^j(jN}uLJv+xklW{Il2o&ob)gV_iMWCN2Pm}=jt}r7}~h)M^c4Q!cBTt zmjcARhi03HO-Fy<N}ZxW-;UCo$;X`$JmOF*Z-D70P#rV#_`7_3g)`*r?q2fn!8|~~ z9Ta=q9Te(DcjqEs93>??L^_HcuvDQqg^wM!eSaezJ7?rzlJV%GlvQRvg7E`P=-7>T z{YXx70G%j^WiaooVR>}bd-0V4?kcK#I83BsIg@Iws;+koQtU?VQR`lwdB>1XIV7_{ zqr{=?D7ii!rGEPxG1tt87rH6qk#9O5$u1qq2P%l3qKg#{v)pVRrXZKk3BYCTZd3j@ z2erta?0hJPPa>cJH;%D|<J9}>PJ(^tVW^JiTtlOeu#xWdi4}V?gI=ej%uzAn%F-!L zXTd~ob|t^I!iTv^reGxSMF&?y*eh;t<l{6ZRF5Z(>~_;hAIs+}%+=f(rq|SM(3CqL zoTsv}r;|KlKs?EI;=Sz?jKc0B?yfr@4?4nODb75>oOxRws!|}InhD|g7z}LU12fs) z(99&=?)3rVc}w(Z8hG3|Gfu*edB>*bl^Pa=p5hTsLXI?{|4T>1o>}9K`WoeMN`==a JSnRGU{}1x%(sTd- diff --git a/NPAnalysis/10He_Riken/include/ObjectManager.hh b/NPAnalysis/10He_Riken/include/ObjectManager.hh index f42981842..8cef0cdee 100644 --- a/NPAnalysis/10He_Riken/include/ObjectManager.hh +++ b/NPAnalysis/10He_Riken/include/ObjectManager.hh @@ -31,6 +31,7 @@ #include "NPReaction.h" #include "RootInput.h" #include "RootOutput.h" +#include "TInitialConditions.h" // Use CLHEP System of unit and Physical Constant #include "CLHEP/Units/GlobalSystemOfUnits.h" @@ -104,10 +105,20 @@ namespace ENERGYLOSS { // Declare your Energy loss here : - EnergyLoss He3Target = EnergyLoss ( "3He_Mylar.txt" , - 100 , - 1 , - 3 ); + EnergyLoss He3TargetWind = EnergyLoss ( "3He_Mylar.txt" , + 1000 , + 1 , + 3 ); + + EnergyLoss He3TargetGaz = EnergyLoss ( "3He_D2gaz_1b_26K.txt" , + 100 , + 1 , + 3 ); + + EnergyLoss He3StripAl = EnergyLoss ( "3He_Al.txt" , + 100 , + 1 , + 3 ); } diff --git a/NPAnalysis/10He_Riken/src/Analysis.cc b/NPAnalysis/10He_Riken/src/Analysis.cc index e23f787f3..a88cb3a92 100644 --- a/NPAnalysis/10He_Riken/src/Analysis.cc +++ b/NPAnalysis/10He_Riken/src/Analysis.cc @@ -31,50 +31,131 @@ int main(int argc,char** argv) myDetector -> ReadConfigurationFile(detectorfileName) ; // Attach more branch to the output - double Ex = 0 ; double EE = 0 ; double TT = 0 ; double X = 0 ; double Y = 0 ; int det ; + double ThinSi=-1 ;double Ex = 0 ; double EE = 0 ; double TT = 0 ; double X = 0 ; double Y = 0 ; int det ; double ResolTheta=0; RootOutput::getInstance()->GetTree()->Branch("ExcitationEnergy",&Ex,"Ex/D") ; RootOutput::getInstance()->GetTree()->Branch("E",&EE,"EE/D") ; RootOutput::getInstance()->GetTree()->Branch("A",&TT,"TT/D") ; RootOutput::getInstance()->GetTree()->Branch("X",&X,"X/D") ; RootOutput::getInstance()->GetTree()->Branch("Y",&Y,"Y/D") ; + RootOutput::getInstance()->GetTree()->Branch("Y",&Y,"Y/D") ; + RootOutput::getInstance()->GetTree()->Branch("ThinSi_E",&ThinSi,"ThinSi/D") ; + + // Get the formed Chained Tree and Treat it + TChain* Chain = RootInput:: getInstance() -> GetChain() ; // Open the ThinSi Branch - RootInput::getInstance() -> GetChain()->SetBranchStatus("ThinSiEnergy",true) ; - double ThinSi=-1 ; - RootInput::getInstance() -> GetChain()->SetBranchAddress("ThinSiEnergy",&ThinSi) ; + Chain->SetBranchStatus("ThinSiEnergy",true) ; + Chain->SetBranchStatus("InitialConditions",true) ; + Chain->SetBranchStatus("fIC_*",true) ; + + TInitialConditions* Init = new TInitialConditions(); + Chain->SetBranchAddress("ThinSiEnergy" ,&ThinSi ); + Chain->SetBranchAddress("InitialConditions" ,&Init ); + + double TargetX=0 ; double TargetY=0; double BeamTheta = 0 ; double BeamPhi = 0 ; +double TrueE=0 ; double TrueTheta=0 ; + // Get Must2 Pointer: MUST2Array* M2 = (MUST2Array*) myDetector -> m_Detector["MUST2"] ; - // Get the formed Chained Tree and Treat it - TChain* Chain = RootInput:: getInstance() -> GetChain() ; + int i; for ( i = 0 ; i < Chain -> GetEntries() ; i ++ ) { if( i%10000 == 0 && i!=0) cout << i << " Event annalysed " << endl ; Chain -> GetEntry(i); - myDetector -> ClearEventPhysics() ; myDetector -> BuildPhysicalEvent() ; - double E = M2 -> GetEnergyDeposit(); - TVector3 A = M2 -> GetPositionOfInteraction(); + TVector3 HitDirection = M2 -> GetPositionOfInteraction() - TVector3(Init->GetICPositionX(0),Init->GetICPositionY(0),0); - if(ThinSi > 0) E = E + ThinSi ; - - double Theta = ThetaCalculation ( A , TVector3(0,0,1) ) ; - - E= He3Target.EvaluateInitialEnergy( E , // Energy of the detected particle - 15*micrometer , // Target Thickness at 0 degree - Theta ); + BeamTheta = Init->GetICIncidentAngleTheta(0)*deg ; BeamPhi = Init->GetICIncidentAnglePhi(0)*deg ; + +// TVector3 BeamDirection = TVector3(cos(BeamPhi)*sin(BeamTheta) , sin(BeamPhi)*sin(BeamTheta) , cos(BeamTheta)) ; + TVector3 BeamDirection = TVector3(0,0,1) ; BeamDirection.SetTheta(BeamTheta) ; BeamDirection.SetPhi(BeamPhi) ; + double Theta = ThetaCalculation ( HitDirection , BeamDirection ) ; - if(E>-1000) Ex = myReaction -> ReconstructRelativistic( E , Theta ) ; - else Ex = -100 ; + if(E>-1000 && ThinSi>0 ) + { + + E=E+ThinSi; + +// E= He3StripAl.EvaluateInitialEnergy( E , // Energy of the detected particle +// 2*2*0.4*micrometer , // One for ThinSi and one for Must +// 0 ); + +// E= He3StripAl.EvaluateInitialEnergy( E , // Energy of the detected particle +// 2*0.4*micrometer , // Target Thickness at 0 degree +// 0 ); +// +// E= He3TargetWind.EvaluateInitialEnergy( E , // Energy of the detected particle +// 2*15*micrometer , // Target Thickness at 0 degree +// Theta ); + +// E= He3TargetGaz.EvaluateInitialEnergy( E , // Energy of the detected particle +// 3*mm , // Target Thickness at 0 degree +// Theta ); + + + + + E=E+ThinSi; + + Ex = myReaction -> ReconstructRelativistic( E , Theta ) ; + X = HitDirection . X(); + Y = HitDirection . Y(); + } + + else if(ThinSi>0) + { + +// ThinSi= He3StripAl.EvaluateInitialEnergy( ThinSi , // Energy of the detected particle +// 2*0.4*micrometer , // Target Thickness at 0 degree +// 0 ); +// +// ThinSi= He3TargetWind.EvaluateInitialEnergy( ThinSi , // Energy of the detected particle +// 2*15*micrometer , // Target Thickness at 0 degree +// Theta ); +// +// ThinSi= He3TargetGaz.EvaluateInitialEnergy( ThinSi , // Energy of the detected particle +// 3*mm , // Target Thickness at 0 degree +// Theta ); + + E= ThinSi; + + Ex = myReaction -> ReconstructRelativistic( E , Theta ) ; + X = HitDirection . X(); + Y = HitDirection . Y(); + + } + if(E>-1000 ) + { + + + +// E= He3StripAl.EvaluateInitialEnergy( E , // Energy of the detected particle +// 2*0.4*micrometer , +// 0 ); +// +// E= He3TargetWind.EvaluateInitialEnergy( E , // Energy of the detected particle +// 2*15*micrometer , // Target Thickness at 0 degree +// Theta ); +// +// E= He3TargetGaz.EvaluateInitialEnergy( E , // Energy of the detected particle +// 3*mm , // Target Thickness at 0 degree +// Theta ); + + ResolTheta = Theta/deg - Init->GetICEmittedAngleThetaLabWorldFrame(0); + Ex = myReaction -> ReconstructRelativistic( E, Theta ) ; + X = HitDirection . X(); + Y = HitDirection . Y(); + + } + + else {Ex=-100 ; X = -100 ; Y = -100 ;} + EE = E ; TT = Theta/deg ; - if(E>-1000){ - X = A . X(); - Y = A . Y();} - else{X = -1000 ; Y = -1000;} RootOutput::getInstance()->GetTree()->Fill() ; ThinSi = -1 ; @@ -89,6 +170,6 @@ int main(int argc,char** argv) double ThetaCalculation (TVector3 A , TVector3 B) { double Theta = acos( (A.Dot(B)) / (A.Mag()*B.Mag()) ) ; - return Theta ; + return Theta*rad ; } diff --git a/NPLib/InitialConditions/TInitialConditions.h b/NPLib/InitialConditions/TInitialConditions.h index 86260d82f..86db26ad3 100644 --- a/NPLib/InitialConditions/TInitialConditions.h +++ b/NPLib/InitialConditions/TInitialConditions.h @@ -94,7 +94,7 @@ public: // Vertex of interaction Double_t GetICPositionX(Int_t i) {return fIC_Position_X.at(i);} Double_t GetICPositionY(Int_t i) {return fIC_Position_Y.at(i);} - Double_t GetICPositionZ(Int_t i) {return fIC_Position_Z.at(i);} + Double_t GetICPositionZ(Int_t i) {return fIC_Position_Z.at(i);} // Incident particle angles Double_t GetICIncidentAngleTheta(Int_t i){return fIC_Incident_Angle_Theta.at(i);} Double_t GetICIncidentAnglePhi(Int_t i) {return fIC_Incident_Angle_Phi.at(i);} diff --git a/NPLib/MUST2/Must2Array.cxx b/NPLib/MUST2/Must2Array.cxx index 26f51bae4..8041b3ae2 100644 --- a/NPLib/MUST2/Must2Array.cxx +++ b/NPLib/MUST2/Must2Array.cxx @@ -373,10 +373,10 @@ void MUST2Array::AddTelescope( TVector3 C_X1_Y1 , NumberOfTelescope++; // Vector U on Telescope Face (paralelle to Y Strip) (NB: remember that Y strip are allong X axis) - TVector3 U = C_X128_Y1 - C_X1_Y1 ; - U = -U.Unit() ; + TVector3 U = C_X1_Y1 - C_X128_Y1 ; + U = U.Unit() ; // Vector V on Telescope Face (parallele to X Strip) - TVector3 V = C_X1_Y128 - C_X1_Y1 ; + TVector3 V = C_X128_Y128 - C_X128_Y1 ; V = V.Unit() ; // Position Vector of Strip Center diff --git a/NPLib/Tools/NPEnergyLoss.cxx b/NPLib/Tools/NPEnergyLoss.cxx index d478a388f..4e1453037 100644 --- a/NPLib/Tools/NPEnergyLoss.cxx +++ b/NPLib/Tools/NPEnergyLoss.cxx @@ -236,7 +236,7 @@ double EnergyLoss::EvaluateInitialEnergy( double Energy , // Energy of the de Energy = Energy / (double) fNumberOfMass ; if (Angle > halfpi) Angle = pi-Angle ; - TargetThickness = TargetThickness / ( 2*cos(Angle) ) ; + TargetThickness = TargetThickness / ( cos(Angle) ) ; double SliceThickness = TargetThickness / (double)fNumberOfSlice ; Interpolator* s = new Interpolator( fEnergy , fdEdX_Total ) ; diff --git a/NPLib/Tools/NPEnergyLoss.h b/NPLib/Tools/NPEnergyLoss.h index 7ec6bb53f..6388a1222 100644 --- a/NPLib/Tools/NPEnergyLoss.h +++ b/NPLib/Tools/NPEnergyLoss.h @@ -53,7 +53,7 @@ namespace NPL const; // Evaluate Initial Energy of particle before crossing material knowing Angle, final Energy - // and Target Thickness. Interaction is supposed to be in the middle of Target. + // and Target Thickness. double EvaluateInitialEnergy( double energy , // Energy of the detected particle double TargetThickness , // Target Thickness at 0 degree double Angle ) // Particle Angle diff --git a/NPSimulation/src/EventGeneratorTransfert.cc b/NPSimulation/src/EventGeneratorTransfert.cc index 018a695f4..85b6a3fe8 100644 --- a/NPSimulation/src/EventGeneratorTransfert.cc +++ b/NPSimulation/src/EventGeneratorTransfert.cc @@ -334,8 +334,10 @@ void EventGeneratorTransfert::GenerateEvent(G4Event* anEvent , G4ParticleGun* pa G4double Ydir = cos( pi/2. - Beam_phiY ) ; G4double Zdir = sin( pi/2. - Beam_thetaX ) + sin( pi/2. - Beam_phiY) ; + G4double Beam_theta = acos ( Zdir / sqrt( Xdir*Xdir + Ydir*Ydir + Zdir*Zdir ) ); - G4double Beam_phi = atan2( Ydir , Xdir ); + + G4double Beam_phi = atan2( Ydir , Xdir ) ; // write angles to ROOT file m_InitConditions->SetICIncidentAngleTheta(Beam_theta / deg); @@ -394,8 +396,7 @@ void EventGeneratorTransfert::GenerateEvent(G4Event* anEvent , G4ParticleGun* pa // write angles/energy to ROOT file m_InitConditions->SetICEmittedAngleThetaLabIncidentFrame(ThetaLight / deg); - m_InitConditions->SetICEmittedEnergy(EnergyLight); - + m_InitConditions->SetICEmittedEnergy(EnergyLight/MeV); //must shoot inside the target. G4double z0 = (-m_TargetThickness / 2 + RandFlat::shoot() * m_TargetThickness); diff --git a/NPSimulation/src/ThinSi.cc b/NPSimulation/src/ThinSi.cc index 432792a65..88e236f3f 100644 --- a/NPSimulation/src/ThinSi.cc +++ b/NPSimulation/src/ThinSi.cc @@ -46,7 +46,7 @@ namespace THINSI // Geometry const G4double DetectorSize = 68*mm ; const G4double SiliconThickness = 20*micrometer ; - const G4double FrameThickness = 1*mm ; + const G4double FrameThickness = 3*mm ; const G4double SiliconSize = 50*mm ; const G4double AluThickness = 0.4*micrometer ; const G4int NumberOfStrip = 32 ; @@ -589,7 +589,7 @@ void ThinSi::ReadSensitive(const G4Event* event) { G4String DetectorNumber ; bool checkSi = false ; - + m_Energy = 0 ; ////////////////////////////////////////////////////////////////////////////////////// //////////////////////// Used to Read Event Map of detector ////////////////////////// ////////////////////////////////////////////////////////////////////////////////////// @@ -598,6 +598,8 @@ void ThinSi::ReadSensitive(const G4Event* event) std::map<G4int, G4double*>::iterator Energy_itr ; G4THitsMap<G4double>* EnergyHitMap ; + + ////////////////////////////////////////////////////////////////////////////////////// ////////////////////////////////////////////////////////////////////////////////////// G4int HitNumber = 0; -- GitLab