From da5d132171b641264723d915a4e6f85e5c22cf06 Mon Sep 17 00:00:00 2001 From: Unknown <unknown> Date: Wed, 2 Oct 2013 10:00:25 +0000 Subject: [PATCH] --- Inputs/DetectorConfiguration/MUST2.detector | 10 +++ .../DetectorConfiguration/Riken_65mm.detector | 20 ++++- Inputs/EventGenerator/10He.reaction | 8 +- Inputs/EventGenerator/10Li-dt.reaction | 21 ----- Inputs/EventGenerator/11Li.beam | 2 +- Inputs/EventGenerator/3He.source | 4 +- Inputs/EventGenerator/alpha.source | 4 +- NPAnalysis/10He_Riken/Analysis | Bin 66708 -> 65030 bytes .../10He_Riken/include/ObjectManager.hh | 1 + NPAnalysis/10He_Riken/src/Analysis.cc | 44 +++++++--- NPAnalysis/10He_Riken/src/GNUmakefile | 2 +- NPAnalysis/macros/GeometricalEfficiency.C | 10 ++- NPLib/GASPARD/TGaspardTrackerPhysics.h | 1 + NPLib/Plastic/TPlasticData.cxx | 1 + NPLib/Plastic/TPlasticData.h | 9 ++ .../EventGeneratorTransfertToResonance.hh | 4 - NPSimulation/src/EventGeneratorBeam.cc | 3 +- .../src/EventGeneratorTransfertToResonance.cc | 77 ++++++++---------- NPSimulation/src/Plastic.cc | 16 ++-- NPSimulation/src/PlasticScorer.cc | 2 +- 20 files changed, 134 insertions(+), 105 deletions(-) delete mode 100644 Inputs/EventGenerator/10Li-dt.reaction diff --git a/Inputs/DetectorConfiguration/MUST2.detector b/Inputs/DetectorConfiguration/MUST2.detector index 4a7ebad47..bf64823e4 100644 --- a/Inputs/DetectorConfiguration/MUST2.detector +++ b/Inputs/DetectorConfiguration/MUST2.detector @@ -22,3 +22,13 @@ M2Telescope CSI= 0 VIS= all +%%%%%%% Telescope 2 %%%%%%% +M2Telescope + THETA= -40 + PHI= 0 + R= 150 + BETA= 0 0 0 + SI= 1 + SILI= 1 + CSI= 1 + VIS= all diff --git a/Inputs/DetectorConfiguration/Riken_65mm.detector b/Inputs/DetectorConfiguration/Riken_65mm.detector index 0c25c3e98..d1a0d9e01 100644 --- a/Inputs/DetectorConfiguration/Riken_65mm.detector +++ b/Inputs/DetectorConfiguration/Riken_65mm.detector @@ -144,5 +144,23 @@ B= 58.66 -56.29 70.79 C= 64.48 -9.85 85.31 D= 17.61 -9.85 104.11 - +%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% +ScintillatorPlastic +%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% +Plastic + THETA= 0 + PHI= 0 + R= 303 + Thickness= 17 + Radius= 25 + LeadThickness= 0 +%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% +Plastic + THETA= 0 + PHI= 0 + R= 317 + Thickness= 8 + Radius= 25 + LeadThickness= 2 +%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% diff --git a/Inputs/EventGenerator/10He.reaction b/Inputs/EventGenerator/10He.reaction index b305e5854..5a8c571e6 100644 --- a/Inputs/EventGenerator/10He.reaction +++ b/Inputs/EventGenerator/10He.reaction @@ -7,18 +7,18 @@ TransfertToResonance Target= 2H Light= 3He Heavy= 10He - ExcitationEnergy= 1.0 + ExcitationEnergy= 2.0 BeamEnergy= 550 BeamEnergySpread= 0 SigmaThetaX= 0.6921330164 SigmaPhiY= 0.963142053 SigmaX= 6.232 SigmaY= 9.069 - ResonanceWidth= 1.3 + ResonanceWidth= 0 ResonanceDecayZ= 2 ResonanceDecayA= 8 CrossSectionPath= 11Li(d,3He)10He.txt - ShootLight= 1 - ShootHeavy= 0 + ShootLight= 0 + ShootHeavy= 1 ShootDecayProduct= 0 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% diff --git a/Inputs/EventGenerator/10Li-dt.reaction b/Inputs/EventGenerator/10Li-dt.reaction deleted file mode 100644 index 59b82da15..000000000 --- a/Inputs/EventGenerator/10Li-dt.reaction +++ /dev/null @@ -1,21 +0,0 @@ -%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% -%%%%%%%%% Reaction file for 11Li(d,3He)10He reaction %%%%%%%%% -%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% -%% Beam energy given in MeV ; Excitation in MeV ; SigmaXY in mm -% SigmaThetaX/PhiY in deg -Transfert - Beam= 11Li - Target= 2H - Light= 3H - Heavy= 10Li - ExcitationEnergy= 1.0 - BeamEnergy= 550 - BeamEnergySpread= 0 - SigmaThetaX= 0.6921330164 - SigmaPhiY= 0.963142053 - SigmaX= 6.232 - SigmaY= 9.069 - CrossSectionPath= flat.txt - ShootLight= 1 - ShootHeavy= 0 -%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% diff --git a/Inputs/EventGenerator/11Li.beam b/Inputs/EventGenerator/11Li.beam index 9bc69caa7..5a06c2798 100644 --- a/Inputs/EventGenerator/11Li.beam +++ b/Inputs/EventGenerator/11Li.beam @@ -7,7 +7,7 @@ Beam ParticleZ= 3 ParticleA= 11 BeamEnergy= 550 - BeamEnergySpread= 0 + BeamEnergySpread= 10 SigmaX= 6.232 SigmaY= 9.069 SigmaThetaX= 0.6921330164 diff --git a/Inputs/EventGenerator/3He.source b/Inputs/EventGenerator/3He.source index 956f62209..c88c8c55a 100644 --- a/Inputs/EventGenerator/3He.source +++ b/Inputs/EventGenerator/3He.source @@ -4,8 +4,8 @@ % Energy are given in MeV , Position in mm % %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% Isotropic - EnergyLow= 0 - EnergyHigh= 25 + EnergyLow= 8 + EnergyHigh= 8 HalfOpenAngleMin= 0 HalfOpenAngleMax= 180 x0= 0 diff --git a/Inputs/EventGenerator/alpha.source b/Inputs/EventGenerator/alpha.source index 1644876d9..e3e16d1f5 100644 --- a/Inputs/EventGenerator/alpha.source +++ b/Inputs/EventGenerator/alpha.source @@ -4,10 +4,10 @@ % Energy are given in MeV , Position in mm % %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% Isotropic - EnergyLow= 0 + EnergyLow= 0 EnergyHigh= 30 HalfOpenAngleMin= 0 - HalfOpenAngleMax= 90 + HalfOpenAngleMax= 45 x0= 0 y0= 0 z0= 0 diff --git a/NPAnalysis/10He_Riken/Analysis b/NPAnalysis/10He_Riken/Analysis index 7a15c0ec1ee2c7c99032279403bdbfc030232674..c071f2cc652182c892867e7892065363517d6008 100755 GIT binary patch literal 65030 zcmeFadwi2c7C-zH5(vnpAYw(-s1;GMHoZ|4wWSSExm8-&ikFb4fd<l~HYpU;waP_m zsiLCd9oNgE%l>3tb#c9wi-^1KT9HLoL0Ls9f?5?3747@|&NGut(v;Wrhwtb8V*^j; z`OcX+bIzGFXD)e?T<;t|Iw>j1P_HgVSA(c+oyC@gwCpjHl4aNpn=#Zl#TaaKLl$XS z^(cc&R6r#w@+`Qzk_Mc01L(Mp1?EK>dGiu4*2T5@F2mr(`bM&(uf)CGYU9-l3=i}* zM7yMS0ckMkgF)lf0*1UO7u(=97I-Y$=2ZyHOVNry4p%I{(J+>GvD$bQf<G@szZ7&7 z@9<Lm&-44r&+}KD>-Pm}>(WBOv<$N@L=(D=nK(t%rJStO9~XHKz(w1T{#N!IZDUS) zYWbyK6wG}(>E5TVy5v9qTmafdxa_!k;woR$|J2)VuD`)B>TfD?_al!d>WzBCNHGRj zQuA#nV%(5F=ZT&@&bAf4Z+X~aB-M}Xxw)!~eaVyo#}(U~C!_tIORdTEXf-)6y=TAl z5$HrqicxPHWmxJ@?UJ>!zTPPAYMI(Q>2Av#chuj!(r#bgqf2>@W9u!KPNyt+xbkro z;&S5RYc#GgxW?kT1lKrRe2q5~V*=8NxQZo)7{*k2KTW2WBfSFGbcva`wR`=x0>CRX zX2`S>=`38cC1&Es6~I+ZfLD!~7_~^};F>G3O6qVk57&H^f%F<&zr}SuE*RcmpDe<) z1Q%ZoW@0qTl*v+Dx8S-J*X_9Oz_lFLow(rE#@)Eo>mGF@(|eKLkLv+k58`?V*GgRK z^{~8o6zT7At&;c&q)+1d1FqG${)mgOr*QoV*VDNEjEk>lajnJmJg&dudI8sqxcGVr z*UPwG!PSgw9WM3Spl*=nE&Ab(M?TA5ocUJK5A}E4xpBaXUIp%J8ZRiQ+L-#6e)W%B z*816gcl+a^i(b4i<+0l)-!@|L{2hOJd-iG1H+|XXQqQPUZaFq}K|{vPmu~CR^@iYe z*Il`#`(1$>^V|RM;ae^F@BQPW2T%U{BOkBu&v@_JX}z~R)Wdgp+}lgPpS0e$v!<)_ z@TX_}KJ>dkE^hhpve!}<Tx-0xIrPA}Ukw`O^`+GO^yfQ|eJ5$n%B1^GIXTPs=ULBY zl)v}c1A)SuZ-6Y>r`~hi(0``f@{f=Dy^|Jt`k(8bT|V!h-7oy-d$(twJJk8F@mKBM z5%}?kGdy><PkZ%&Z5P^)`*v!d{O(0lOTKvI`7a+zUGUGD4|t7|`CFzueAC<41T6(G zpZn_N*VR1p_w84g4gP!eCExD;)8dwM(taO0_@4K@AytR|urF_1%U>S*W?`Vmzn<+| zc(|>1xOiRh?Y8f_|26HbyWdRt;G%~DOHZ41X`i!(?OXKL*IzF@)j0Ua>+bpIA0BJk znlb$GQMoDUJ2QR77v0?B!J95vIA;GVFMLz;r$3JSsQHDDZok%gW?A3M-+O7#zPG>V zV^5y;&bII7?f7W_kkBaSzxMn8_^&(Ny$)3N9`NCVzu)oQHeW+=TkrGlKF0Cy%DGh; zzYpc?^1k=j@jDBrR?c{D+&{*(+`0Iq!Ut?`4w`lSz1Lpda=0z)`OQK5j|X}L&#Zay zLifM!oL4maxihC<a%KDCXTDr;@X|hi`1A+Y$EVsetJZ$wxF@3^cUF&rZU1Qf&7N<w zJZ+kbh}%)NXik6JrZJGo7x6zB5ao;bM5{t(oIJaGEbbjw{?w#c`h5s(>P3F@Eee_Z z-if2HiK{<0j{oIx?Jth2e_oWmVaNT|arzF8(`QwjJjceB{~)e?OI-b%<K!ueqkq^V zDiHPW6@@~+H{<9%<LZ~j+2hQ({<<tq|8N|?MREEpkJG1)lc0Q2zCm&2TjS&z6xV)v zTzia_>P3DxD1ems;y8JakCSI$9Q|)`@|+RZ-qUgNXT;G9;`*mOj^9sl_N|K3H!aRT zCdbKpOI-bb$JNh`qrVg<&+&2i;W+-U#nH`0mT6Slb5%dXxOA_@b`0)MGkP1R3^t5+ zKQ-5M#FkS~2G}9#Tl!!Oer&N_EAc_dYd_s;;~d$~7-WnaWEdHrSZp*C>;HuQ_V2RT zxHcmG{tT@3_C@*kPBDzi?<_W&sh`o;7<{5(YzIFR4>GPqd*?yE(`5aRVUHQ_TWpk_ z{GLX8rQbyPy@G}4^g~g3Zw0;nC9Fg7>TeulEF5SUXQO_qte=NP-Qy@vvy<NzG_;Cz zv_m`_`sYD^h$V3Sso-~p)ix5C`f^@=c4t)I<Bl_og9k0P-m?6f6AWV)^kFx${36KH zjP|Ah_Xq!rVfdMv{htHQ+i&Vax&w-y4|}Tq$vGZ$t^E@qf5icdtzPounr`fWEVh%7 zQvNd7zxb;teE<STGvuK`C{GdWYfp{Zd*|s`qoO{CDeIp@dA0VwLH##Au-Jx4dLb0G zX!2dchQZH{g#6!vyxU<9CC`O$gqAO(^1KPf9t0h3A=<kL9di)=rR;YF_66LlfGq6w zHTq*Y{D*o{z99NzhvqMj!+#f|zo3Soe+&Ecg?!3?0%-pV(3Lz}(O(nM-^#wr(9v7K zU-|D22-0Vd#rCH3hq2%vfPAX_*(f;u6(t{RZHIrr-<3YgpwA$*H$=*FKKx<Ww-(zJ zi5~*r^yX;&X=uOrYm1F*O!ogR=>O%A4{j~=u_Ik3+oMjuhd*xtKjpulL7zLfS!|bp z#`<T1z7+ilGm84RoDBUne|r-4r~zIj=^hm9iS}S-QQnF2x5Hmm`(xl3y<l&Uu~qxB z{anb~qWRZJkT+HLC+NFb^S6!EU-R$Jflr3LRDb^t{dptwDVFX38UEHE@~Qq>i1usH zKdL|f3_<&YpVEI4%8%CUzaM;t!@tTTKic&rj8|2DX&1xT2Arwr@R#>nqy3u+dB$t{ zUxSXVhkk1OUuQRrPf=dU_Zj;8e%Sj}%tQ3gIq-iQ<WuA68pONp&@U|8I|GKb04w|U zMt=>1{1|qk{(8vQ_HxueHp5@)&_6034FJFWC_h=2w?W=wZG0VrhW%Q6Edf4oUbO$d zhMdXGQU2{<oTuscN7!>b>VrM}sgLnn2)YsWgxd)DzGM8>;_ubSUk`gKee%$r?k{&z zF)f~43Bz^Q{P7XAw?@We`mZ1J(=xP&FeTa>h=E-Rd+xEA{#^!twMl=p8z&fdg8r^H zzuXRfEt>yrfTN{o{&xch#t^G*H|BlDt6mViRHJW&eAjCFO~&{djd-$I%JYidFb*Le z!Yzb+7eL=f*IR5~f)C~IiTa;I9_4SpLwoZO4^;ekkc3ElLu`R#KaT!%N`9BYf10)a zNQRttE#5qfavwlG<=;-|KLG7rDeHd^e@fB(>s|<2uEnn>V2^3=Cne9ZXnzL$SH-J8 zV>~Yd|4EYHRL$PWu=m&K&(F2^b_VRd8}X|i<*C1cvy_5hjTgr#hU=0!uF2k+zEId( zQ{Z=pLf(*39j>hLx+`2BcfjNI8!p$(s$jqs3cG8<E|=k&RucAyi##P^hr?6pu5pEH z+`e$A$m1+=<Bzj^^0<=nqEc6Bx~tUTa+c(>%EaNNlgej#J>kGy=PVPKgll|(nFS7K z@i-51D*VNBN@9vu_?^XbKncv9CDM{`wZHRb^)?|vj;nZFNv0|7e^pJ&>&Psf>a7X+ zf`NicuV*&I^ik~G(vonn#y!*PaG3h!`~5*rtV)Fr=Nw^TUm#2q!`hRJCds>rUX?c% zHWxiWu9Mw?ieObHni1BCWC&hpH>-4P#%R-4B@UObq5}UZ&$yCs&K#IESW{F1^H(^X zSs~wiF9ZOw1I<@wvIjRu3^&R+rp7%(3d!5d3I(d&>bA7TE1ZTJMU^WKhQc9Fjjvj1 zHo=Xq<slu}7FW&-`8*-1>9|~&mEvY*PH@jO#nkc&gJCB+&&7)1QW^15I*@52Wg>-@ zwM&K1%0?^vxWN&jWHUiJJK)qpZ&(N*OBPL<91N;{aHJDVsts4y!ntR9!$kpjJqCeu zIQpILT9PhYR4+KZ&>ITZ1m~H`j39@iKnDgTAqq3u>-L1<Q_>ol6N|^^syrwrdMze1 z!vS(dK`=1GH?y|J&AOv~ey`JH=L~pjX3iTQ3@P)nt|L8~4W1<-pHH}MJWUyHLTxCV zQRoi4r9}je?KMaGgefJZ8TmCe?s-B?sFEKD)cV~uB@V6X7<Z`JT~ndu6eFaBeIAtq zr*veLYMF{#QNS1Wx&2VI0yn{cQpk~8+JO|UgP3L7vmh9(sqh8dVXvqN@e0O|brzqS zk(Qp8m6jnj3xz8@=bQtOmSct$F-Sr=?&@l9phCu>i6tR404EHAiWrycnp}<$=k>ZO zyq+3wl{XMBj>!%LE4;39cgRa_s5aaatb!H1&RMg}$QDyN$^=7KW3F>431>l`gt5*! z8MQLen$P3HutS)xGDEEs5l;$`?Ln-rp)<KWb#7OU7vVfOS7?*v3D$;LOoVC`T!hB- zy%>pn0i|x|l1?lM`a>czmV`4#L@ebCguIw55Yy0lx!cbnRNjfzFrp@#ala(wE>VVO zU>jko=?h|Rfn;%$jWg35oDos798*{T6IR$$DfIA%iV{>wU_4JLrl@R87Mz33zQ;Yf zpL3aMlHo3wupm0Q!sTYX7vr=foD0!J&u7s`g<FzG*pVeApCNlFk@){*@obMj7*cUP zBzq1MaYEae6LA7<#jNZs8Sa9D<wYe~E~w5)8j5oMl^z4<i<lDX#Htx_nxChpDy|O< zIk<kVrLuw;o*bR9Y^JLSHlKuyDhxZ(Wj{}j;l7Yd@G<+2Zd*Kha!Ig+%DZOP)w!xM z-3J41zb`z`H78xLrVq$|Lfg?u5FHvZ+4>lUE2+2=i6z4^qu~83Zz6Z@&}CR1=%X85 z7TaTpaaCq0b9;hfa+&GDz-N}ja7_u!^`RFhdqcHVUd}vnO-FJ%u}6l%qDcrA2%WCl zkQXZo5it-(>wIAY(s~2*2^N>jzbdycV7R;$m@Ysl4~5KpWCg@}#F4G$*9q<bmPs{^ z5i)*bRT7yR&4~`4r6NXANeDr&+U=_;f}<2q#9Wm>A~IDwoiKB{tHkTa%<Tw;i{Q1$ z%JO?d80I0aJDKBjI9wB4n0a072q)uBrfNWBiBjxX1~eE{#&sX35w4()D1y^_HG`8G zWskb(xvmn|9Jjv~Vjm&ecCc<?hB9SlmAe|Y=wP}yGewPDk^zf3vqj)EndDB*pInqb zYP>U}%)y8-u|so?;RSxL8&eA>L33UAf7yN}CtfimqpeS4>!#CrLba|LzJOT2V1|iy zrw2=fu<K{`Ty^#TrIZebKRDN0<0=o<Lfpx}mJ>KK;J-D#nU&!x%<ZPkF)Ph6-Y}Qo zk`h1E3L?&Mw8vEr<dSoGZNPUmyc8oLTdX&-$4t&I9-CV=kJ}kXq41ldK00zRI%a$4 znPVgB1RUm(F@4k<b!0(nz27?tC?!sWW(SPnVpLS|Fn*ys0)1z$%gv3XBV%&@#KK7v zoD;_sO>_c@4WKiCWsLXNF=AYZ>BNj3$9z_F%ySuUI^PJ)nSol&7IluyfOoFBERsRA zL)481oetp%A?Vk0TTTal?+k4;9d##(@yy7u6SGZZPw-17o9OB=!u06Ok;5^RXcI9# zMtH{x95bdf%rT8wDN_3AND^~UafG>FtssxO0E%%3msf<N-lFQ&U$p)}T=@0YAKIAp zRpGp>WQ1#cc{-L<nIiOy&3WWpse{2XvPRYV{1qjv>NOVzZogP1sEPN#Uv5=Zsk%qm z??hb3))uQ=bD<cS%YIR}RLhC_s5=W*ypdVA(^VlSA2k(6kDhWd*%VY_sV+?D$SJM# zhTR2jzXvBaJY9(#!kHUy=p9*0d;DSxj@6;Lx%#DRF1R`bOT5g+T3swJV@`q`nUNW? zB;sjDahi!Nm9bgqxG^~L6s79oXl<f%@|erUPbw*K<W$!L!@)qQyJn_0{2T1HOr`W9 z>0`Z_B{(;%&i4!AxTB~~Zk4;vrJwcw_e0EF7tg46Rq}L*JAJh$Iyw<Pe%Z`L(y!4> zHU>(Cc?fHc=uTE@X6wIjr=(){QMXbqq{jvnhdW0QJF+{<_VDcp)73g`0OSdBq;SS% z=0>iwfp|o(#}-;_%wn8A({%b6r^axK*_<Dds~KOocZlxru)!0Dk})ks#!iV(Xj>Bg zg&h~2&5y9iFvC@raL}mvCU&w9&Pb0v)G*&Mkmic-6^A3gqC%aD;5fd+8kla09y-c{ zFN3RJL?+ino==NoORmoo^(mrMTomJp$`QLUYSg1FHn*CQGl8RS6mddCqlgnC;oe6I z$wzP)r<xNY-jc;FoOIQrY!q(Y#JNDH>wsS;utxWYF^fXg`M=b3!)35W!21=C)1<)= zog%#vyUjWJW6qID--;V+Y9Mv87;+^(@pP8+<XJ5ECyy&<_x@6AtNrzn>TAn?^HAl8 zgWs7eRm6K-#pA>%Qc5N~y`>_2n@7GOkA=oBT~(@00lMxN&IY+SM{@BUy3V3?uGq>X z+Qj^Vg`={4BuHdT%rNI}eVLf8uO!E~=l?oEMC;?2(fq3|Y@%nozqq$@%u`36|0g=H zAMWA^87`NVztAJ_P=q?s?YOb+{O$MuixD|OoL<+|dUy(x!*L<Q;Q!mKTi{5vcgf;G zKA5=iPOxyDqaWDfMolN(am1hjH#@@m_NdSF`WqWZojIH_vNbmk#QA>I5rP_ayyzhF z#5`t|n-N)i{c7W$`wQ(|o})MoIMOS;k(r<v4^GW8#*8l-Rp7#77si-^0#_zend3;y zGz!Yf@<$c9cq_8Y%1|T=Po?n_kQV&qzvxRLGc~FF%g6J&iIR9wQq*BB3;x+6^_ipB zGe5y-=BhfRC&_=&JC$>EQ^lX3R!lOwN80@*t(ZK_n)?0E;bqE!Pc3%x?@K((h{tUW zJFaA(4-X4=H!kOchj_8hs~zawV(42ndhw8Xc!bevcq8@EI-bQO|8!iTQ+z0GGpaym zp4nbkjb9hz8$44=8eeZeH5*-TG3Lqg;j?bXGq7EZg_8gNQ?0gccoveMeXx!)jKjd) z@ocVf8d7w<g(+Od#uR*eF$KqDrrq$2CR2PiBZVnA4`vDp?MxweDpP1QlqsI69L5wD zN@t3PVY8T`{amIHdL&akE1Sm@&ukVl#WSd5nXbSyyiD=XZZT6l!&=Jpd3^qtDV`C% zf+?Ogoz4``^j0vf!LznZ@r<XRDV`0kW;)L>!c6fjd>vChBg+)eiq|v6GuaE7+6`j~ zQ#_Nrl<8@BmXzs(c&42xo>g7JbPk@GW{PL0S2D#js*f^-tyeKU6Q6-(Iuy?aGd;^N zo?&`6KKID<9DLr9=`cJi%k*44qs%l7&z>_q51(0J`W&7GW(tSc%G7~pv6*JzGXzXC z@mUR~2ot-RW*bH;(;Pfw&2%`P#b%m|XQ-Kuz_Y(h&qx0=y#W0Wf9Q^f!!1lNG>j<? zAD1lHZ)pX^xcr6<c#_g+yyI|td;6j-VQcG8$Zh>M+@m#&x@DKYzM(z)06eMvj6?X# zI_+oBvnqMJcNeO)pJAY-O49bW2y7vy+wN@^*hY*VH1@6$xEHaVc%{I}#6yXf2|R!p z7BTkL3!Fj>4UN6k0uLt6Bc3j>op>y9vB0Us#l(364<#-mP8WC>@pNLl!0E(PXm7H> zS;W=EhQPVRb;NByA>)jZ#P!6x1kNK~Lfj&7A@MTeW`V~NuOMC{@Oa{t#480ZCSFCn zOyE-DHN^D-ml3Zet`_(T;%4IM0#7I2Ok6B*1#t^;p1_sF+lkW!_7m?SwhLTM+)A7* zaG1D_*bumm_z-d1kJSHKVnjpeFK|6EKTou`Mc{?R$;8b9FCk7LUL){QVmt9lftL{v zC0-`*a^iI2dVyCE^J8aws|CKFIFER`z$=Nz5*G{nC~+}yp1`Y!%ZSqjUQIlm*e>uI z;!5IVfuA9+CN>0KOI%0X#&L4Si^TQBy990~UP9a=@J8Zg#LWV4CSF0jM&LJyR}!xj zxP^EX@iKw860aey7kE4ITH<PfcMvxdPZxL>@n+&;fp-(P5a$WpO1zyoUEux1yNK-q zw-L7zCkuRlxQ*Bl_z>|S;<g{C|6yXxJkVcYMwm()af`r=G?mH3&A^V&u4?$tg57P! zrIRc9If_bt*zD4&mHbRX>nkX5U|Gyd#;5uVSA&gFo8Nfc8r-ZMV7J=NO~#*fhq^SC zwU2l=WN$+B#OLMLCpCy7?E~&4-7G)i-I`A}S;-=4!A8sFS8ixH*z_WzwupZV_NP3Y zf})F>!~Hp~4BYg1&63Kz_LkQ-EQ1V?_%~BaZrHySdM~)nW(lXgR+wt3!v{}LkaKKn zCXB;8epsTxlG<8;!SVvR?R^N-T+_nT+i*|9yYQk$XRig#mPVs#xtU#>YHL7Ulv*lF zEk|{&(z7PKT~hB-JZyr8$d2-0@L@|F|NBm9c)MYf#%VS(8`>J)Wq~(gZ;%vMmIu5; zR(GgIdWl9vcC?X&veZ75+BaA>qBwPvjqF|6U+At_evZOuLDV(%T`t+|{!U|~iob&5 zly&d1CL2}wCR1(&%XAS`W`U~vN31&7G41Fqw7o($*j4sdukS_nN}0XiWGzaW5-Lg; zR6u<NKl&a5pCI*Cng2v41>V<BYw3BGR2EXrKPpS^?L+Y_skvyswH-!i*wC=O0eutV z{r#y79N;C~!jy+o0}DKT2yBmQ*w*k`!{F4`<t%=v;j6tbCyNTA6e=_>6f9U_9V&p8 zs(%{viHM~#dm~d=d@k5<3{b2}O{_Ysc_7wuQ>-PF#T09`6zjwov3>*+Vx4J<B}$oM zkp;yn2S<pd>NlbTB4Vk`Imi@Z?ZT4s->v`NweyP9x_^I;zt!SzrT9Bt{4EoIi^bou z;%}b#n=AgNi@!tZ5W~!#Un2$o3R*Jzd^3Bk%zh8q5KojcHC2=^X%;AO2bQIc$v3PY z356Paf<>}%op8BcPO>UZt!}`uA?r(kQKS7WFjQRbqrH?Lsm#SRoAM(i;MFXt{77Y9 zAxrMvghjWeRkIX-aI99FWwzDSN|Z9SQj`;=Rxb_gOsmUHttypP(>1LWmsh{lv{IRC zkQvcR33wk%Dy>xJVp(!;8@3dhRxMKe(a@UNvahzw?DLQf(~452R*K@Jnb22zu(e3s zSLLQw)0I{;HLVnvFAi#2smyng8PQ4!_$o^(tyJb3S#s~7Sgm$R@vEV=rd6xVz6{yW zN|Z9SQj{Rg1g-jYs+G^ws!VBBrD>(OoS<u^GA+7RO299^(X>*T?;=x-&B9(Webpw# zzZaY|tq#fT1;~b0qLitXqTIu#q16QgI_s;Tsa3JkDx_(pxD-%UHLg@<nv_+<XC>eX zx>hRFqO=NNldE~4A*0qt&^4`WGW(CnhE}4Ksg<H^AVX;7J+?EgYE7-`lvcmhv{GEE z4rp4b%ray~JWvT(z>>-XRc4wjxpxJ&@S0Z1QvA<B*R)EJ+3z45T8UDoR*Lc&8A7XD zF+nFDS2vkj<tnXi)3j1t?*3ZSN@d=N%!pP>z$%thTB*!3S#s}s98hRl*`@d?v04q4 z*-55WqLitXq8u-^dLDCgVy*5mwHm9mTBT{FxUBn1(@JHoMrK4SCE(pGskBm=H_DQG z_uFE$N|)k~g4SAJ<;v`Hkd3|)rA)09C7)(OUwxF^nO4u5TIDIN)@fQPF8jXJv{ISd zkQvcR3Am0Wl~yWqwJf=JV60YoQhYzO*0dTcv!@^%T8UDoR*Ev4W`b6|JJss%rdH`n ztG6_*6qn<4tyCtyJ|L=<5^!Ierj^RvhD_?eHxEZ<ng<q3@$UjBO{+4QeLb?Fl_+Iu zr6_l@X=pVZ>z%|N_%Bndp-QVAnpTR-#gtXWK$Us6lvTt)CE#(oRw~m_TFvSvmS{7j z*eRyi0v4JpoBvRng{jkJ3HpNu-}%DShF_M_urT$0l+vj8pm0P=#bF^0z!i`aB~`;u zArlK))^9YdTd}8XOomx8?+B~gkT0yiAJWKj4Onh%iFl!6G(&RVy9`I3KZoOYaU6e6 zonwB8mEz~r=@(ZAOYmo)RO{W)QSA)w`3$Dnj%9ahl;h?YTQn7>)`2e8oE>Zt*)n|j zIrwnP(1-)RBd|(Utm2IT;nFGG_n0D;aQ^*LIQtPe+7-vaaq33*Go?zn&v8hNHT(Uk zE77~+W|g{Gqi)vXW~Vhb$!bHnW&rq15D7ZJbv76&C{2=zR*_XHOPmGogs@^uv)s%r zO|`IYcUZlV9=qA(scP&2&xpsW%q^6d9(#&v4jg_&=4%U6-#~r7_P)zD%?-H{@LgGN zuLYVmE)`|}+VJe50S%kK{dC>;y&E>KJA_UAS=bl)npZAPG8+DN*7}B)Ma{L21uvSM z!)HQ%wgJ0Z&@19VqW?ZK=9~>}sm$%b(LULT(v2J>qT2d(-=Cre5_vGF$g-%~{#IjQ zs+~OQux!|A&4Z!e-TT3=2m3B>?t9<*Q(G1`hcj4&A&zRY=(#L<MQXarrv;dQoS8p# zU$0K~W7X6H1<T<{WGL78buxP#vYQsBic+BBR#Dt^Jyh(978o<@sbYWRv~ST$odYvO z`c73Hh|I|Rr36S~Ni~0|%>AF5CHGDRceN8;V45*WnR7$KrgcA^5;gW3N|LM@bj$O| zF-B~=w)co_*IPzxsIlU$4!kMxdf!Fg<J^L3k3_eru6$=Nx^B7IRnj%Hki}wSW$!*C zu{wU*2;oVWn_Sw0w&^XirBRp2?1Rt>HWH;wqbf=-RA}8}l2u+S+KQM^WxkEfhzV75 zFY6^$<`b&qKyf~XcGX$Zt+cqBuT<uIWJVY-kzVB0ORCIbz2s6^@&cArhhQpmh%8By z{H+!Ix8JalM|+oF<!ab}<pv)3?a##_9mhhWQw*Q8FthKxa=S!M=&<w3T`U3?!+~~8 z*tO%j?K^Q4a}OBqys~Ylv+Z?q66c0&80UsbeHYz>x%7tXdbPLXd!H691z;~31Fbx5 z*^Mv0Y3;cOJ!eV%nD?K9@mhX;6AIa?dD@7=q_EX~l)*avCyNc|At$vWmsEazE~@i8 zu3E2Wb$&aKUY!(H-vTjFa6D#?)?CQJ2)@%gmO=E2)RdjpV)IT+u&oG9<l&)HqLJ&Q z*9Fi(!}dq)U<p=!vib}&s>!O%WH9XCtzejnYW8|qVW(B*r-*#Pp7~_IP_n;IvcLEf z^j(TN;>u@7VZNg9$H7QwnBU~o<opYB7J$lI8RLHk7ER7Iq!dBOD3>2tQ*wC|{LyYQ zv}60dI4CIx)n#R93bpcR%11vYd%<rnTEaxdC~jtKr(lfylcX%&X?K2W5t*rOBO8h{ z4FWVdU*YYW@>a^cM##LDlsnnRiDnx@y=P<_B&Hlu?6p$tdW!woZZ;-ZW@uuIj8iqS zg*e5?=pZ)Pu+c}*Xp{3b(Y_G-HD&5oh1kucybp>JTd@2Q88i)vZkRBUwGpmi?Pe)< zGsVskVhfhnkfpS?$e1Nr(%M3tdyvsVY_ega`$fOKDcTnrys5-)5n^v9<xI9Atu0uN zmTi#uBxLDa?9EbaerHtc^&injg5^duro<K*e@2FK5+Tlq$mk$8*|5=UA@)0>eR(UL z<UJvF3n|rXLy9d}-YVN5@oUIJv87MP`F4vGJC9;-{*Yn|mc7YRT3cisM(av!A<j98 z#3may`diU&ABgtltrUB|5L=vVKM$5lY{7E7Y=hNC9#L$rdO0NPDRw_0wqQAxETz~Y z!>Nfa#0ev#gSE+qjlL%OZI5VQh`mQSiRCcE$SzX4unlQ#!SYnu28oLxOCoEVzP(F| z-Au8k{fE{TEFVH+N^FsF3o?YYFK==Rab7}32eHY9jV7D6vRKWx5coA=t%F&CU)Tp} zBB6&AR`>u?Dis9#t7Lm5u7_MVY~)0RSy0Y_Dd1w^R0W&Uky?l0bXd%4nE6B*RZz~0 zFF_>p9W!L?%<F|>`1}6r0u2yo!27SizfR!60uLs>(tMX~e+sf1r`Q&(7gKM;fYjD} zG`rJ~@51dLii!)G#iWae8pO9r8x793hA9U=s4}Np6llz|nEb@OEqZS>PU*Go!>&m! zLTa;Q!+NGh!&_&q<4U`$CW)}w53rfp`cyn9k%yiG6R|#?VnZVK30+|uJVept+zr`P zW@RLky)S!w9x@j;*N&8MFkr)mCg(n_-iOc&%o?u165L#YfX)v35d^HXyK}c=?6a=L zU*fA@a~|Mrlk)(*Hd1Fl>zstLP0p=NQ}#nFSaw&O<fns<#u}Yl8=SkEoTAQ;eA`eu zy2#kCy=8d+2NCaep;5ex+2lR6U?Y4+*c2uj(0C<<Xq=MFwx~PbpTep{XtllZ%EPa( zx5(#CWLZ_{+_*w~=T-F1g63qcUy@m#dO|@9stAXSwyd~y?}S@2r87pX$L>Y*Ub9z4 zcb4tMa@(1V6}!{cP;4=KAN*6HS(8(=PGh6>zrh4T*8!$^59lhtMEjs4mCDqHK4fpK z6}^P|er%T-Hb_-U08FbTqw_fcRA)A1rM6B0mE77Q_BO2u8Lu?&@IE#<3EK$QEruTJ zT2G0{FRhI}Qzi<t{p|}j4!B(MpfW;jX>ZC~%nq%k)u{6GSR0Ml1luia?z`|sQ-9`H zffe33AJga@K<-UW5!xrQfladN004|M6h4l37)D{~AQ<PTz{wJL_i$NvJDP1o9U7WF z(2Ma6t<(HUn}a9=6Yc#0>(D0WcC)Ze>(~d-tnl{j4PE#C3#HJ@6o+53O(~`1A2RPw zHW(}7wb+_b>kK#+YdTw^X6cU=G@QybOgK+UeVE5f=@NWF5K$o<+rt<l!=?Bmg<w*a zk0?tIA)lP0pAl@B{U!^>iuiB5e+w&CHJa1;V<@|*S#(zTT;v_ue$Rm)?e_xpn}QBc zZk#ZnVLa@O*a0I)+!4awEX&@CvLf!WV{}_POl=ySyJ2geu$WZobIP8B`ZT$)f47c` z@YxgE+OBN<M%30kvWVJRdB}4j_sAM@k4Fvv4aaG6zCl5nB2Fpxr%^BXf>f#(XJWn= z4HM6lvX4T2W#Y4GvLm`^i7CLo%NwV>(dc}mbNikx3uZJq+d4Y3)T0qIU)OzoDTxfT zk3g}8?Rtn%F*Gka4<Z1V(NxCFedDDr52DYYZR3O#^tu_kF)|n7QH2F<Nzi#9%Q-QP zi(b2(Ug=EU*GJ+e6EC<f+30&G!b9y7azruyb0HKPY0wVBtBf<a6)|4yyoZti<BuUj z5lg-$1d{6j45;V~KETvPMOHRnJ6aAoY_r=ha1IFPP!h0s0GRPM#bhhClB<!&F2;N- ztO$+oz$zy`7A<jpE$d0%&(MRN2Y{N!S{kPu==+@WYnTnz<Cuvo1Fm_OlS!|}ya7$c zNvv87`yCK1A*z2Zrq%@;?U+B!C?1KpX1KkG1)Te3n<E8k+vNN<Y2Ou%n6BX*=A``* z{KXtk!eDPx+bt*^=e*D3eefmDZ|kr7)(H2ORsV*n`;L)dGx<38FSu^M?8JNMWP+y{ zE~CIx)V?0e3Nn{VpowNYK&Do>5v{?FMVp-OB{fcfAMbB+?!Y)k+)t*d;zC9acw=-n zIo|^h%(Od%mEb(bgZbWvqht1T#tYNwo17wA??Gs$`a;IP0L1CN4)2shlF^%-A4MeD z-MJ*I!9hs!ktxY;B}o%)HE5G_AHv5zllF$9T}{&^8hco<YO5D6Hqki++bN{2p$nW{ zuBfqaHT_!5<g(l~@5Gf`AQr)7F@g<7$HzwqnNdy{?sD}cnQBAK94qD^!{Xj_;{#1) z$xU$b=nB>b{=y2#6n%v;A{I!OgMgxxVIXed&yk3;!924NwC~B9=ORMkuil0trsvD? zU4y!<7i>e%;Vh-KKX2qwLC(M0+abM-Ph$8)W_dnE7#m5jen&9r-@qBG?{oNTTT(0| zoNmcWx5imgvxr$vttq(cfowwTFOWvp|H)VjHM}oJV}FzposWUshxgklSF$;j#Z!Y^ zk%>^o9fHY)Et#QDj+{@o&@R0ii^WMJmz&TWhlUIG4{NOD5{VmLE_GTL@vTmr(qWn9 zz^YO#utk47h7h*z<c4<|FxjF!<E<NIIro0gPTB+0FW5Lt?CYX(QCn(hcCKJlE|e8P zSV`_d)&<U9eQ$5>`&=_7I!1Op8D2=K;T0`O%}Fimc6DVld1yvCptiCmiO*V^`fDRW ztS8?CBa^cB$tR!eDe7rug!L}Q`}$GqnHD`KL0;Dcd0)VEW?gX*>xafo&W)7mwzs&I zY3A)>+-`F2GB!E?gX=?FAEDIWNje8!gASOE!_3xp$l5gJ?+jb!xO<T$_M(I=+y{P( zWg0D8i3&|qUR3e%YltDn(i^6^L2FU8^B^`3D8+%IX!A92<WShV2XmfTi|t9bd!AC5 zk?txKR0jwm^i7a8_w|;m#gQx_d5Z%b>LYmjz`JSDqnuHo#DNy!g+d>6BQ>DU2{j{y z&ST3Fl!S8d1wTj2P0$ur5^p@61!FQ#y%jOS<Y#$ID`Q0&&F?nJB7a2@itq`f(ODcs z!+Y})^_jU*<MR^5j?O~Am3%&wd>)MB^MA-^3i!|kz$&Utp(%8)y<{#2hL8Y8Q2qjD za;4XI(KGNP(0lEZ<+;9zEAJNNYx`qiu*+OO#ckdDE)0MR{4R+#4{B<C@@Aa=aVGtQ z*pixsMqB@evWyC}+fT^HopX0ovA=^17m<&NR<SVO^UX!eh!5D0s7)iMZm_ugpY$#< zi!*-U-LtnP<x4PEt|W@?Vrz8MeUL*s=W0rp(0RT>(@josa`p}^!&{W2<^#I$6kjOt z0IIWDP?!Uv44|p>R<n$ca_EPBWl|10z2JY7RzhUVQ6&mZ!&sL=VrqLK2%>NLF4~JZ zKWet+FoSIR31KSXmnXw@qzNL!;QKdBE9^TlTKf~^H919ecmcX9X+_49$e`6Itq4B% zXeC6(5>;ZL^t`#qX!0Hq$1FISGb4qFVOPqc)FlZbo6$WNn+E2b$F1TNSaP=tHRm7) z0fi=Zw(UDZmgA2CHRf3wlNYQ{R_~aa9I_P$o@|II;x@cpEkcr5!%7C<V-ti0{~8hf z@z$Bleu{gmnD+mGoQ~~(#&Tx+17h2MR&j{5-!D=7{MD@l?VqP=NA*We$M$)4A=|gb zwx6vyMB2X%PhWJ@AKsau{il!}r$2H!w*MK+nfh}bBi^w#?f<Oe5NW?(qW0&bSfqWi zx@&TZJyX>`=!oCOjx^yLnWlp`O=(j-djj4R<)Q+sKv9DJq&!z-w1JuM-2RgLDJWs? zU8GNnHhUnK#pr=me~<Uc`!IxJZE`nqI@;vlEN9vzCDtYnD-IEx{0PZA^vMhqi`XRU zlP9X$5uf1hqGS7)v7Fid;Mn$0RU9JiUze!;PY^b=_EY%wZ|!LL)^F{F2;Xgdi}x6W ziVkn_l0x&P!F*y?G2b8lk`+GPELT}rNs^`?&)<%Ppd5p0{dNu_UHsaq5u$aTNB`j> zngYD?1^O}?h37{@;VWp$T<nUzSFsYgs_(CUT`c0|QWooGELkJ_#RIRQpuUC;q4d5% zVxAH+!2+nR-lshiG;?t+*0APAIWgb+<N5v&$M+nK@0&1XEMIwA@*0e07-B6N-Gk5$ zTv~g+7Ty+ft_$^fZfc$@-x+hREan;O-#E6^Zcnz-h~snpTw7SYDnWleq2&$TmnWy% zlfHuV(Y|>frk9jMUxFfzCN%miFya0SA3eIRm)zVx8}9{t3<bKpmWQ|}@3;&!^Lj<@ z@T{fG!!z3Aj?ThKRd}C%B7@@*?km5g_HqhY1Hr{`{Z?%#m!PJMDX)S^HKK}*gT3ST zmm<*IO9RM2V?n>i>|2V2FseTj=+h)OTs($S_j+%L-N~1DkI%Fbbm9wbT)CksJj958 zehiM-Iti`au;Gle7)b`V{(d7Qet{pKY+As-XT7$rz00EZaH`{=d2FFhQq57*zWmQD zwO~`yh@WcSU$Ako*r+uQPDPVw4<2ZF*oF_IEjZL3KK2FP%~-ZzpQXK}u`9A_u`b|8 z+9G9JAqzjtyJ9F-hcCj#FT>|{F+t!h7(_1lnis6K2fW@2dpKwhReL?Y8T0Jr!EmMB zZO5-{h#%&#&%m#1*zrY9$W!0hW#)xyJt!EmpKs5{m&D^YO+xk_Rr6H6G<=_=$)IYU zs+U$Fe`>^DQyVA^qG&}L?hH|XzPEy3<!k~lJI{eHJe%yB?G3mTo!_==;5#+>rH!(* zMED9v!&!G;A^v!LVeyrY=64GlPA7AX(o!V(Mj8H=8JFSjWkzYGFHqvct$BlQ3NgPG za)QX>_d`xF@OAM)e>7W%#^_qVNkaL^z2PkO1LLnzXx=h0@ckzEy%PRViEK?uC%-i< z`YRvCjG_Y9*~Sd{iy*j_cZT@E5yLLNL=^oYOpcy2^3i!=x8EL|VK?kvei6C761`IH z^#=4}{5I;~%o=x9mEBQjzrbFH3%`942pi`2REtCw$WR65*+q6np^-V(>zXhRKW~y2 zt_urPm@(5m-{mNGW#o*DWaj%Lw<SJtYkrfDR##?ps4nv*2TU>eUk^XM9ue&)oKLl1 z4^R(azH}t+>C%d*m0@;0qnDyFgQv(Xc$8}`o>^UqhpU(1;~3TW02jZ{!C=s11rgpb zR`{t^@oKv`ahjTE=3i!}dAG)-^$XMTN%#KaA0n0bcgqJ=jWRY*oLRN*4+W1ryy4;H z8@?|%{KS`ACX^Rne9LFguKae>nv05_dhy<`-kE<<zfaG-;hwv%zVPF#j!(9axaxxK z*S-5q%U3U-pEAAp<%4TZIlsQL!Mkll{fMIRKiq!P1G{oneZ`OM<68lKb@^30s@)L} zO#a>I4{qzwuA-~*iOR*<3tt*>Yo4li=e8lO$E5YU=tz2Ys1N<WdG~S4yN`AsS^Tv9 z$J>AZuL9M6leD+uFY;wS>izeJ{G~hId*-fD|82dd^`Aen;qd{^(>`}JK6XZJp=!Ug z{6B5EpzcQZ+Pw2W`!e~eyl?XrJ}>q1cX!{jKEM5v>ob1bzCK^k*Ug^g8}a6E^4G|E zdoz<KPW(0{PsyqDQhILv^5)jV1AjmA7Ae>1adMq@X!^4c4$CQs(O0v-^2>7z2JibY zXVAqXKf3NWJ#JWZk!nw+s-Fa|AC=c%`te<||Bj!!=fS7CZW+0{V9Q%q)@0xtDfo)| zsgn1a0sm*rf95Z|I9`vGa|=s<IJHNfnVuli?=F7d`SR-{UphVS>``x@vG_k96pp&P zy7iAQ*5vKmn_ZOq<u`ea%WwR;EdBNi-&$Qh@v)&ZFMR!_(CD{M-gn`l$@6nhJNya0 z<-(NXju+l3A2e;(#ZPY^a@@I-oAO@#^PKk<UD1>`_pe`X{QAVEyu&>l59}DR=;FWI z&hmfT^^YSDpI3AG{uj;~WhuG0%{pb?g}>|i(V1I6y6)nO?k@a3`1+Cy$DQ}xU&<T4 zzcBopn_fFm>A&c*d$;Cx{o+sO?|FJ}-+5aux}a|UB{zQh)H@f>;H&@Dd8(h!%6Yzg z(_dy6v}P{2cGj1VUZj3oQ;3;&YqLl(Ld<l4N}EkQJ#S)h>7+^H@f)?`ceq05sR<eL zPV{HA&NI$C(#x>Nls4A@hTR`@<Hu5FV(xL))C6mW*#kk$D;ir-(6G<oFGUz<nA4c} z-K<D+b}_ZtMW$VeAm=;8WBjj&BG)5;|9>?P(7s70cD*1E^TbTFo!s2sJ{|Ky5{^w* zA{~dc1&4js_wZE(NKZvN9c9l#dM(m4q$`k)Li!BSIY?WO-i35O(!V0*j|Xo;%7^8* zAbk<(Zlv3h8au!j=>Vk3IDE<7+1{Rov;t`{(wmX`k*+|x1nJXAA4R$yX*1IONOvIZ zkL|!Aq-974;Ms&Hk!B%%25B+UElB-H4<TKG)b;`RApIvEciM>bU8K8^euDHc(sDf5 zm4bP*8fg~NYmgQreGsW1>ElS3AblF?qe$oP1|OuC;Bmp-NT(v@+!{hU80l}3<|4fh zX(`f|kX9qzjC3i|w~?+w`aaT)NEhRY%vPkUky?7798x>dkCBc<+J>|Y=>R;#%HIed zjdVHEX-L-~tw;I>()W?JAvN~4xA#K-U4(Qf($Po@kxoN89qDwW*CGug<%cDIi}V?! z3y}`P+F=FKu}D`Vtw7q0v>xeJr1v9z@pJSC(vjFe?;u8ch`0^?fsNFQNYjxPW2;t* zv;}E3G4?b&@G+@Jksd<Yj5Hbbwjx~$dMnaZNDm|3h?JY7R-}(0?S+pOA4ZyvwC92L z_VGy5kyautM7j`ZIntF#tC7BlbRp91NLL^|fRvvAO-8#H;p3{qkQO2xi*zkgKhigl zE<ySe(npc@$A@*Bkq$;W4C|?3NXH^gfr<G_nqNHGNV?Y6P5dh8iM=eAB}t^&aXpW> zHqoVqk>AHQYEhToa92@}lfh%45y#WYkaY`aShjVbg^{%ww9EuF?oA&A?UV$x)yVrD zXsHQko00b<Xy}3t?d?X^y$NW8Sr#woK3*u_-TX;ZQ+xOTXl!fM3+?Uv9af9c-)d7f zn4H|zR%`TV8`R^#Y02H>?V;0>yUW|drwJQPM{K+dC+e55AypT*wC%N^^?_XYOl}9- z3eYT|otS{e_cTb0x2+vzx8gdCcWC(-sO9PbQF$+0VBzCaZ1(_Yhf#kR?jtgdS<sa< zE~R_nE$=fD(1w9F5Hz@7v>t9*Z!Bmw&;})-Re;v56Q1>;C4(l;C*tgLzvLMsi)?Q# zX#G*oska9m*r!`T8wuKEoi<j}AM^m|2ooLj;Aaf30PUm%v|*rK1=<juChayBv`Wx0 zq+)ndo(jp6zvUdUw<ga*<afb5lE`;gg7y{WoJ6!2LHiK2YxOpzJbVss6nG};s{^2A zb<$ok=Csp6J5A?_d-m}#(9Q-eQ7jk>+Q0-nh5n#*1C3)eUVqXTg8qL%e-nNgPiMQ2 zf_@d|zeF}@25ll}iTY~?Xp=grcL=nLI;l4R0qtVY64@XfG-oF~$AgyD39S;e(-Y9_ z&~FK3;W}Uy?k!8TF*`2mv#U`5E|f{+dmBMp30k69u}jht*@sU(KAb=%+B^koiDjUj zj{8WT;FkTx^~nm*96C+&=TgwQj){*?<Xa8e??4-(*TuaZ*D}bn8ni2MZ&{<+p=g11 zvE?QapVpxMCn%GsZ{7fHKWJi;6BkQbLF<FHQzGBBVBPaQ>Lrrb4%%0sL9`e<P?wRQ z?dYUl8E9KUOB7?mpsi29(+-}?At#U6|BBS|qGp4Vs12S${T_Hjk$s`}4Q1k)$ttX| z66vxZwBk<W<x>P>JE8HZmb^}Ad7$McpxGhIbhMWa8h>NmvhcX*Tr?@F3lDw{VC`Cu z^6|Rd4_X=4sflQ7L3^<i+E&o)Sl7nav!jUvkZB=kci`S~?7*l@+E^Qab?1Qu^=aR9 z(8_SwT#x&BdB=m6hjnlw+gE~?n}B8qn<bDp9kfZfw+!tl6U(ednafe8jAgR4GIDIN zObg2JqeY@jVaGC*pKEG<<_Sx|4)SxI{Tyg#=rlQ>aP8d;TB3QH>+!XqC5kazqyH7O zMDq~W?W;jcWG}AeAMJ$3_5J;zP0{7VJ@w<BZyRXE326I4<EJ*`$Fdz|srv!c<J#0x zrP)sN3+^WhI~iNtV~hkXCMJRxW#YbLawqk;7r7*XOsvQKi61^*iu-t9<GJXgpvA`k z(LX4jG+(S}lhj(4{8+XXI>zhI{9If+K%0g%($++CoHDDuJ!!us72jVJvc8;T{3Xfy zNEc&!5;)!8#rkO%<D)LtFNPS)x>_IaYCP1H_glr^eO;}uT8u>&>wOkuxrO&HTS#pc ze-Bu!_gjrSt<2eEwPHW~FLB>4{x){A{<)j6x*MrGyIH^PX6zUDA9rWY;qKO^9>!g^ zkssNN9X8%C?!o*Q#NRDFtRM9>_V*b1d{5))p1j}DlhoUe;onuqSfB4@tUqStt-XvT zy?Fn0FXrzNe}C#_y|s_AqW8#;dmHcf=KYdBq`oBnZtY|JxUcbfpOG*0HJ<6q`}g~j zdPg$<KACKNp`WoadE{;VjGO!M{+WKv-z)ys_qX15tZ{Guk$d_ZyTtv?$CCPr_`B^` z>z)Bd+p!~G8elv(fcLuwkb38F{JZ)%E57P$({Us37-%dV$ouC8GJl`=d*kueJ5r4M zj~}`Bc;iEHzchu^b>i>#6zkp-jIUBgzH)-G_5|L4cmk<^8f4vYqR~9a`tcy+eR2QF ziJOx^xaB14n<p86KdB4u?>e~;{{G=)>;99C){~j@{wdZs?Z%-~tasauWp>^_d+MvK zw)AxCkEa`tpFa5Q(~UR8{lU}O&ReO!#ozy&VQn5_{BTCkyM`FI4&nXtLr85MV!h`~ z<MuPH8_qOdJQFUplYiiF#|?f6E5dC_)-RKcPw8tnb-^Eyj8l`2?}OX3tl^|lR#(@C zq@H&q8IQ0CDqJW1)s}7@ZtyGyteTTD=HXBw$-1$N@g>W>)Wv#pSK}QrT+r3}vc-6y zD{eP;wZ7Zccw1!r#bW)&YHYQTfyH{8)wtQpjL)ss)!hyFCSqkb>*j98x^B#Pth*JC zxSn6l{8@Kvdw1iI$be_w+tc{YhU`asSfA@*Jk^64cl5Nrb&Rp5C$it`Y5lCH@rlUz z+cDPVy^TG`Ap7=S)<=685A<S2TQBo?9(vrI)E~r?tT-Z>X8m1~alI8Qn){NhKP4G| zOS0}rGX5p*zfEHPLtXgy$u7@yeZyvaY_($GAF^iS{;%D-yxz@N-_45qJ>5>h-^aUK zzwU0l(%t%2cjIs3{w`Y<@;|eY#UY#Z;U31lJ=pAvJs`k``~z|N4Yk%8=UXo|E<t(= z?TSp>N!Bc6i2GhCNYApS8_vWyaTkM0uM<(>T<c}VC1<lc(dQ>69m_sve=kk4-j-xM z+*yB%3jh1_KMVZN0{^qX|19u73;fRl|FgjVEbxD20lc-_e$P&eE%zG>P8qni<MXWM zi^sO=s=y8JgEHk8DG)D^>9II_C*~9SOwW-iKW9wL505bA*BvpvQl_(HifyOJ-y+iw zWQskg$p22J{3tf-&z0#WnI40)DCVCiQ+_O%c&bb*Wx7nJ8)f<*new+?Sf0Nr!<4@r z!<4^4!<4@jgSYvN`j;#=o_UPNmqC5=Q+sgf4Hj+!jQZOxc#jR|M2jsp6g8F&F$`?w z#iLBvU5gh)HQG+K^1Jzs6#U_y9$y_|;PEG;@&${H_hm@=z7ZdR;`s?X<)`-c;EKOm zZ0MRG+rnFy!Q|)7{sR23UW4)UqN=b7x@^Ps39f^<y5XCu2I3lmD+kwTTvKt)z%>Wg zVqACOdK}ktxHjS1hU*ht2XS@7d-nrz4Z)RzYc#H@xMtv*gKIIayKp^@>p5JTaBai& z39f^<y5Vb62I3lmD+kwTTvKt)z%>WgVqACOdK}ktxHjS1hU*ht2XS@7SN{*hH3U}< zuF<%r;+lbL4z9(x?!xsruIF%V!nF<8C%A-7#uOBsZy!2l;*_)OS!r2m>Gq8D^bzUV z={feHlkvdcSa;aWKR4?vgFh<uiwiqU&cidZV>vs{&FIMa|0jiyY>l+gysEIf+(-+B zYs|D#-T4AxZ%wt476^vDY52{pbHnbL63q<Mrj_&O+Ryb>7->X!fU?p^tC$x6KQj&2 zn3Or*nh+m06~!cQfmY-7v%Gmz?GGDi;x~NL!d^TsnKlCeq+kUfU`_K@x@O=(Os}i5 z0yPx^OOMCpt@C)R!!8JqM@7>-_<h<C1Td3X9u!9-?kb-LrGsFMKWU;hlXH0prSV&i zRd@<ijjx}Ro`^z>FN_IUH&|?00OH(KRwj+F({bU5Ta+7%a$`ZT1FILpmw54<QkAPj zxk{8vl{6-*Tp{l9s6deIGKMWhxo1FReBq1bcumB`IL2~}ZFvZCc?g7zb*LiBFjh^) z#rVc@jD__GVvK?6#j<=|g^O{F<ru?~M;OKu5LLU(=QSG_;~wuB`&OKfb86hPU6!S6 ze9gy&5H0F3c9vaW#!khL<v6$(Xyq7dD^ZRynB^GPRk@|Or>ty;vA+I7oZ;hMy;Qq5 zAw!j0gK}$7E>96b<1Me-wQ_CiEVj0Mybq1viBd1xj;|HC6hF=_<8iDd&VU6qs_%oc z9PLtxa+N4Ij0TV|72F>OCO?*2j?eckzY%BL>>&A4qCNwx%5Aq=ZQB>&H^pcW`BLRx z23F<xxC9kjt%#th+ZI_)G{w1QrP&|I7B5xq9U$^!IsQ&e@=CLx0L4p{dml)Z;}625 ztc<qBf~wp{xM!b{ALqEiEAgHIDl(5RmZNR<<5K+i8*Td^j`9FWlw;X1wQ`)Z+fXir z73GU%czusc@w4G~Xl##k<R_}25ZkAG{FH7ll-sJ=*)9sfA5l(g#$dS%v<!92wE}sm z*O*7~KByuhQ{6J952#oAld$(mSf?R_<*5UGkS?V2%X<ya$mJkYz7#)&^^|dom|uO# z7vmGtuCisuArqS|Fy5G$tz+CEu9(;ywTvGoHhWL4b-IekX2t`PE{cgYlZnlN!0~Nj zuEY#6j!n!+hTn+8cr~%v0UVbm7Vl<>wVjFkFcIU-#NsIsF`i6pj$V!%6SG(F;`qR; ztHG5VUK|H_b(N8w{%&IPT^;(hiOrRwTBCI}1~S{ApPF>86wNiHiTUmi{0jbL;u8d> zf0=lY!0Jy|<HU$xFwcT7SQi$M?I%G#)==^#VTC#rSjm@U9A{($+wogbNYJ~zjd4g< zqdoR~5^7IF%3p|L&k|Q6t=G!ye#G`US@zP}zXAAgKU4lBtd4I5*4w`i*p9e|M6~}W z-pCjEy@Y#yjjNJ}@#td>e*lcn<??NDA>U!7A8fPO_`W>#8HS4dN>`KK*SH9{KfZ*} z<k#02i<@IL%r6ug0Bq{p*Kh-;XqcZcAFSbPfbAOQM*&kcd=KzY4gVf^n1=a`Sh|Lr zfwMIH25_#1`GwabHM|EnPs4m>Y%H+S*NR|Ktl{H;%QSoj@N}*IaNr6JbDvkKVea?* z8vYG%wT9imVGRd>>ohzM_*xAw1g_Wc?Z68)`~dJ0V5o)u?Q1**+@Rt0z&C4n8*rnB z_W>^jR{i-C@G=b_gAle{!zTl;&~O^?{TjXqcqOp1?^NJb8V&-l(eM)BwHjUt+^k`K z{o7^@Zvk%6@NVGk8vYh|mxg;nu~rS+f!j2k3w%h!#Wci^QTaW<HVt0`oUGy7fKxR5 zIIvyAYk`Mq_;uiP4SxWft6}cT^EAx2*2Zf16yRbFrvaB~_(I_68s@WHl^VVRxLU(} zw!2Ql*8taRxDj}XhVKJjrr{@mS7>-G@JbE82E0nc+kn?-_!Ho@8a@cztYItS)MgFy zd9D@>JAk)qxDa@khNl6yYM4JQ(Wc>RfDdW-W?(~$hxY>8H2gGhvW8y;PSNoDz;+FP z3p`ZAeX&_d*YK&pxf&i0oTuRlz+*Lh1#q#3X91UKcph-2hHnP0*6@SCbsByVxL(8C zfR|{P&lxS#aDNQ!6&g+lUa8@Uz^gPo9e9m~`ON)V4POu3tl?XMH*5GI;1&%(1-xCu zuK@4TFrS%f)$qr_Z5sXo_>hKsatYNQ9q%c?HVvN%oUGwo;1msy1Ga0}13Xm2A>eck zHv;Es_(9-24gUprtcF{Fi#6N|T&Ce4fv0P@KbC-%8a^GkTEm&ZbsBa8*K7Dn;3XQa z1zx7%M&K11ei(SAhW`w_O2e-LuhH;3z-u*(X9kRB4f8?6%^L2HC2Wg^PY2$v;Ss>Q zG&~NtRl`>Rw`url;6oa|2^hPO$b9txuua2%0#4TO2H+G8zYT2H@E+iy8s=s@UBfBB zxf(tbI8VbPfX8ZhJaDmwuK+I7@YTT6HM|J8Qp5KFS8I3;aGi!<2d>xfd%#OHd=Pk< zhWVL~6&gMpc%_C%1FzDs2Y8K!=K!zO@M7R*4c`U4S;Kz-Zqe|oz}q#v2Y8o;dtwXI zs^Rm1+ceD2mK@UXbZ&W)qWym@uua2v0w-(uY2XwMzX5F5@E+iy8t#cL)IB>bHlA&8 zy*>9du0d2G?u9F}GMi=Kh51cYYg(xt<$jaDC6O;ZT_C_uU69wo;_z$MyOkoOhUX zbG?3|hPhrJ46NGYdfl#JuGdpF%=P+E4RgIdOv7BSr)!w&^(+l@y`HP#X5f(;=6XF( z!(6WyYMATwu^Q%jy%-o_TiA>1^)d}}y*^#TT(4JZ^|@a6YnbcxY7KL}9@a3|>vbCD zdi`1rbG=@#VXoH~YMATwB^u^>eW`|Lp}+awlkN}8H2N)|4^sQRF1?M%=pS}-e`s_$ z)_5)szYfgv8k91Ux}+E%0IU6UQkQ{7D=>e&*rX3Iz5(VhM3}g*(FI%j;;(S-hm8L4 zr!?T2eHPnkz^4JP0`~8+*tl~i-U~c)uf>KiC-8I}0r0n!`MDz2_X6kbx7c{5OZ;1T zL;xsG2mT8D;6aP+3|YTDjy_l(A*k}_12<!TO@aGie?2aa?v2BZz|Ug5u{`Cy19<n> zQT<*5R_7m-Z)05j_u}wZz@wqRk~dpBtY*>az%815R|E4r)0DrDF&CJ>JB?+nu*Xfn z{b4^GNeTI02Y%xNi*1<1DdrIgs|$Oj0b4$a+WQXRY41nzAAwUedv5~nv&UlNXHckr zCXQI__-zV~f7)-1dE^q|?*Tp!@~Hma9!GBjp8v7Mb}ed<Ul>O;>ii|yIKj99xF_V{ z=iEqth(|d1`v2Wnm$1EN;9jr~8%)M|#T&r$K8?y>fFqs?#It(P*&nljAD8wbzn6gD zZ;h5e)jR?sJ7J&PI2?|{4ZwY2Kc(MKz{}u&2-8B|^Kpc9spij50rOX0l)W|pA3*$N zdz61O@cCMQoS_{#*-<77cscBWZWa75173#lK1t%Ifq7oF6qxpT9k^8MzYl>e-$v_S zz&A|fi~e>6@CV>mhI`uQ4&aT@2i-34*0}Qh5NM{uo(@UB6}aWgsJ)vRh~$g>cK~03 z{!#v!f}0xHXBB8{ZvpT@w1;U|;AY^_u*YVJN8zA#qUN8M^Was!$Ug{t=QfM&65O-D zZv`HP@ddR7zbAnC3qq{I^3TWB{|R{Me=IhHF;V_<INsT4FID2-(h=p0^=|{VzY^v5 zH1HPKm*bK8SnU`;n!WRZb$=ZToQL1P7-^XOH5*vH$o~&^g-m(>0nFcmQRP34qjxzK z;|KkXWv*!dN?`S}<K}8${s0s6$p0>lPW(9VcGzF_&xgR<5brWUqdwmO^L-N2e_<s$ zzI?I%$-wF}{Uqha(Z>SsM}JHPjr^Vh-VXbOB|i6b@Q42?{|*4F7x~u%*P%b?hOGYp z@G#V$BJr{4s9o^?Bw+ggnd~_EVtps@Fl{_s6-R$D4);M|-mJynQy8e_i~L3apA3DJ zzs>>X@0qCnxeGW0{!4>X-WJj^UZw%lKX(BeTKOM|<%|6L4Z--;>@gXbzr@5Mq|X9Y zFXF|({J9k75#J5`ImR1{v3)#EZfu7>7GT=vE#MYye*TFC<%{)G&cytU`RxweQ{IWd z>P31bFyGf>UVpr^835jb`CaL|6F5)v*B|5R4<3ro=0iT!-%j8F<^whVKLV^?Z0{N1 zg(y$Apgt*Q8OCt<tBO|>fImn3R!MIHRxj3nC=NeQ4EgLxXpc95yKC{f;B3Q~0e@p% zws!@vdNJ?bIJ^P)TCKmb&QYb!{+I)NJNjSg`wH+_#2Yo<z5>2NYrm2cw|ud^8{_cp zz^h;n6@R~rqvxCp`-85=%Wr}E!~d0iUItzcdtC(H)bGT!==^KP{5u4gUwlgu*nWUm zzKHJtu805AZH0dUe~S53>C@#r?1!-aQ1o%Y>c#p#;Fl0zlsq+Y^qYW-wfW>};PIOM zK94IuEM0Lo<;?@$q51oz!1daA3jyDV_)&~%lz%s{dQrYi7<LBwM~$CbfFFmxPS9BX zE8ta}9}sYeb22a{U|)WgkK@G)JOK8+QkH)fSiQ)9Qyl&f_+6A&{&99DK7Wh$d!j7) zHv%8T_)zg?C-5NjuafWlEUc%s_*?|6UTn`Bhpz=*5Bsb3md4Rn1K*GK`I&6Ce}1;& zX4*FdT(6Dq+kwlq@%vR=`EfZ|?`Zn1239Z1XLxGDp>XYt8EGEF<tiwhG}$%2sHD{8 z0xFC_mlSBEilA$zKUnVeyTlVFE_ZDmzNxIL+K;bkt4PZpk>M~%aQQ0gh%ycHv6YJ2 zs;YUa9p}VC(ah+{`4gNjh1iU`kG3KnXA!NaZ;VUxK;UWP9hozynh(|#RC+zLkyR2N zpX+k@g5oI}*R+Y*1>?s$i^n^%y?(DJToVkqJeBU6nO-N#W;tB05Wn=z<6?22+wYq% zz5>pLudQ=cc<~5Ot@tWBSGv!K@`_oux7uGDB8TzWvr8vW5qCMWN9T_(5w{~|PjJ^w zDL^|cl1rokUqAszx+qvu8LSDjq^PGOe9xSs2`*EL@OTai4rjgSy)Pi|&4$MZ12ZEH zYrq#!z~K-L6Xs;Gs3}P-Hj7Jo#_UnPaH!Z@GiqMg%Np6UCl*gC;a$$`DS?nz$dxm@ zxCYNlne?K-j3DoFXaB}s6Tnyfn1#i+{g}6-@Kt~8Gbmm%FQgGh<8dcH366}}MR@#b zbg*WwyQV_qWQv?oZqMvUR@UqiSTH~2@do&nLZVo<$Sd^5WM<AT3{DJ&i)(@&FP?Z3 z!WD%I{K0_G$&opGRE;~}shm_k3(g|=q|2=Q06bS3syaef9^pu@^7;L~D)@g5T_BS0 zh~{T>M#zjIq(h5vB2Ll#I8ISQ45#c4(&UJonA&DmOuZ4aeg0sFy6_N{9m7D&kKvM2 z?XC%VJ+Z3c`LmjitVaagfsVPk?!Zh8z+jcmB726nrpl-1WR=6Y^o$Y0;t6umbbM|u zTkeoq?elcV%m~EDkx?59t72egRaK0QfMaZv;fXEoi7W1jE$)aZoe^6&BerbD?3lut zvtvsG$4XKi(?dBm_)ac&j5!^q2X`doRCs6NYd2#|AA<XL;E;(dy;N3JTt-HyMyF?j zv5p<->?Qb8Fc8Bhr(>2BHZcMHpMa3g-tL&6UFod}&)3Be!=htum3OAQLvBVOq*rys zpp2>*nB{eq6zQdh<J&bc0L%3}2gd14x0=G{1pVO{YDTA&9FE$~sWW{?;F`nU?7VI@ zzUj0xYW6JX5aSz}6TDM(nKNQx=Il5;CJvXz;o?}FOFj85rFt7#@fl`ph|kTK9qVr~ zD5EL{X3cenV{D%j@XiefV~0$7rMKEu?GAa|eqEeQL>e7C(%s>zV5l1I6Jx3zcU5&< zmUQgI1UPmAf*3U&^K<D)@eU-c7oV9GpOG1?j4^Na40qf}meXLud`Q_LUqt?nxg3%m za<k`p%VQ_k?65cFck33-p6RX)#f&a7YCGg&bO&nY={0j^d;N2KG2NI|6^wP^5#?Ub z)wSrEm|o5H)YkZ7`j2iG!(2IMOlHIv<8xuCcm~5gSde<-boJNbX>DDHOn+^ij)%jG z3b2^a^BfLjmR1I9L*N)gNC#m87csHUoKYP!-DZd924iz_D&77Wv00hbzL@zhH-AoW zrn{mhs5dy=<MR7vR>nv-!s7~EUF)v##)y_%>B4sz2RyU&>JCRGhEoDUHuGa<HAkk% zO2m4EH&L!5N5wT)l{*x2d1AD2WXW$FPDINNcx%NjBLNZPx;BniCbD#`GOKFSb)4(1 z4k54wVlo{F_f;;;64idJcJ<;J>G;3^R1*7^+8Qqk1wt_v&Tx4APzgI8Q$-i&YjZE7 z)~p$s5$y}Gr$caa`MrUe;Yz(FhhtQ*vL-=o1T}wkrQ22G@yA(Bc#f;elhKjpnC%*c zIc<IwzMeX!mF$wzVppk`@z5QHLNUY)5{qlFl8>ik&K~0ql~#gO8T7~aQih{8;N!<H zysj#*JAh4IAmj^k2#y&)X;l7rm#fHy<$Lbb{K-Z6qsBWk${cJH3xDh?VFT<~DqQ@2 z@;P2t^&DqLddFOR;zNB+>I8SdJ=0qwHmfoG9ocb33NoB?urv1M<mRN!^oCv49#^=s zHZVJ_ypGw}ox0pL2rDjZabh#Df#vQ|vY#}1bcwUnRZ7{hi7&itV*Z4p0@M(@`N+;X z6m(T00QkM!t8@^?F``IqhrN@7!La!9SIReSVo5mHHM!gs#>%g#Bvca4s&@NoipG_M ziYI1HEY2S>Rq5q)mV_PYt`gx(j!?L$B-@3oEI+ynHJv4Hk>hkYToYW?wP9BvSmAXR zqf(}7KxBzhY<V0=!y_ag*Cse(C_17DPVd#u$>U0LV2NKr^jueoYmVDr3$ZCti` zSGhanHLBdcKtdB2ivFfh4(wn!9(^+~Y`yNPA`kK$Jfwh?-Pm*&d7REdhjUI|-k9-4 zqYBE(T$yQErb3xj?rIpbgYDz&9POo&442dC%r5uN#9T**^ZUvXbLm0|b-sXX!<CVi z4x*>d?Gi;@9(UMNnU{yg@d+rGBQ2Bd#}AT;C7}Wb?H!62RdgL3x+XXmcEJGh`-7ei z3I*oQDi%)G(IJNy_`O&>@k{E9EAfi~o)D#%t`94+;ST&JIx<W3(;<`*r#L0yUud^X z3_f!hz}Ls;;%5lN=}?9PcNH;Ta?CI6cb4#kL9@~{h7$3KE@^j9s8(V+LPpf_Mc8)z zOdqeVRzvd_s^xI_gLAz#u5y0UvU9Q=E{=>+`3nO4^#t12^tp(eAEoi1f$?9%g#ebq zVY3lOI-|~{+A!k6@G;&nf8Ie-635M7VQ@QgAeRpZTNoLbGO+HWXJ?B@l09Z}e(~7c zs(I!wAaH_U^usBd`RfehgCRO(Nyuj&8y@wbnS;nP+dI#UJ5je`2#v(ctkSU=qaj7{ zI8TWKC$Say|LgQ%7A&QEP$xmW#Oaz0Nf1}~nZxMdjdipU=x@`C=3sSXOwONJIB9}& z;+Ud|P9RKMI1>sWfN?e8!Y3D_td6>&<3fBVrtTQ8vZC>dqsMgK5m+k(umY^Bb7W%I zXhuOfiFKIuqA|NebVmd?mkYm!+jBbhduM2)u_WjZIeqaV{pZeH3>+Qm#)QvbEwsdp zl1SLnW@oV|`0p1LznZ|qF|28Ou?*@&;F#ExQHUWRR)k`G6sdA_LUUHgH($hE*UTzx z0?bu{TuL2b-N&`cuM%h^V=Qj4_X&4!st%*@R}(}z09HCUtgFIFx0oCfEPJZ`iI#I% z`(jO3QeK3Gdb+FB;W8JgSVd+E^ZufxF2bQ&>c%WaE4-c>FF)7KupIG_=qTu5$D<xZ zwb<svdxS+LBV6Ol)16M_bBdL3ymdCiXJn1yN9^%?O?Z3JT>Rl!jVnR9B*-n5yA`a! zD!>d|Ue!Tj-6^*#KX*BWqx@gTeI1-9(HHm$Dq;QTdg2%LceT|2=yzAT=oYc^>2znq z$y_amMnv|AxtP51=BBx~m{#F0?hyJNIp&XIx&0pg7#BJ#7i)CBznfX2d_S^B^!RC` z5H>B0C&Nq4xjzyH{Mbl#+$v2*9!^@hrNV;J+)<Tc>zVDsN+^M||I%eQd<?=R-jiah zONnDwHjd25(xyZP)@j)C#EImZSy$)6aYrbKAL_!%glkTE{J!fbt|KDbQ0&G#?!1q@ z1~Kb;#TQ6W*>vYb=j1V$jh|Ff;>a<N$V&N6d-yj#dDMs|e$Hk$+vCSOG6|Lqdc2)A zD^YjG2!Wr0%Pa}k_^R`HU^^u+7jn4pOKdndh0zNNu#d(r7c=KkG-z&>yUwNN%}eH> zxtaL=G@Qfo)Xgh4gcz{CU`U(~WHEBlpJScCjd!GEj=2-(Lhrv{r>KqLxX2L8EDVO} zwHQw`s$G?H#WWqtaaUJ!XWPM`)$+Uxr_Ua*KTe|IrL=8e=C3x(#&1ixb@(+_yx9nM z;#rU2JToqLvVSpK|Al*OF(Y9tcmt}hRnH=Z<jfHfu&4qHg9@iJvhMm7w9SRGI61@c zI)avwgPm*-pQ;hDEE?n;4(L&Yb7gvkJml!G_L4JjNq9I4<Cxmp`8r4o*cd7zV2gE; zvjoR|SR|EXxe&a$N{SsN!Wn*6K*<Wizr{wFF|G)MZxS-9uv_lHF`|0LW#Zr$M@#M! zF$;Al+F9^FV${SLr&wLY1gcCkT1I1)dD<CsdU`~z=0wTE3G+aY2TNR7(Pd*=iiE@x z^IaV+G)!>MR12S9=;JC9#2<y2`QkJYVWyL_L^-)<Dkeuf*EQ$fED;c*JvPzou@1hN zkuDdAn02EEOH@|KRO1x~lku}`2R05zenmy(css$VtK^6EU7b&DPcX{GX{qWBI89^? zesp$e8QigBkhvzppZ}lQuC*s&APD-Cs1I5`8dHriQDXhx2-F0FkO;p1y)%2SUhg=i zDPOhd-QL{J?#wngO(6y|@L`vcZSrTK6+oOKeG+)?v<BW4)i<lQ=1*$iYOyMnKs<O5 z#f4xH#f9Jvt#$|XZdD;*OjQzzmbo-_m?2dF9BizvuY+ukaW<HXKk95P@XuKnOOhTu zF4Vt7oZZ_9C|2w+*z&u4l^C`f&gI(lp-h|>E{R68Ed%?cW5QI8j#Lvkq0j4B9o+xL zJNa^33-uVzciv53pyawz0ey?BRqa;0{nR?ct!}9d8l=ZWN&IL9BosYZ76e(wG<NtU zQA4NC!ef{mi7{c^ZE>Bbx61wTR_eg%{g~N63iV*kCN1Lem?>GuEzxqu2;fl0=$>cX zYhN4~0XX&fmhPgQfn+K#3mu=PL~HeV;)(RwF0ye{%z|MQG~Bv=8`qJz`OQxh5GJ^w zq)_~30CDdTL?xE>H_nXG1HARDcnX*>n6Y>>4jh&8>Gfl^slsBFYKW3J;LK#6GeaPR z<wha8S48m=RzC1>HC{Y)J7w_y8Ckw@RlabQC&z_b8almK10P<&dC?XEZl2Xd<q__; z8w_YGjCzfL$izLMJUvim{odde26i#CJ=!R7Rz*Q<{@gJ6$r;PNp*U7vC8sW82`Lwi zIoEXA6YijO>kfB~*t7|rN>i!2Xd+HH-@*{WlgufgGwlLSSaZhqf`bgm?P1ROtBSqK zV(Qm9M<wtdpCQvrbu~OMb@NC3?2>640Z4T6gWohaI&%AfaN~80`eWOOv0Ok9+Pq`l zDa%y3e#t2hbv)IUrBiS+rm$asna8-DySYSt9_#mqE(~(?8J^3}A*F;Q3}qLKQ)2+Y z!$i;SY%ryRo7WxNIcqjAk!^b!^3--^M@x?seeB2&IvWV*k<|j|aHaOMInoNN$%Gt} Yy5$G*RW@2QwjP<HaFtswOAL|p3pAmch5!Hn literal 66708 zcmd3P3w%>m_U<Vp5U@x^bUKO}wIV8%rZ0*j7D|9Bk4g(7iXm+S4YW;dQUt_nnNp;d zDk3O4;@}%!BQqkRqLzmUidqp71+)rE5!5;&qN2Uux6j_mPMY$c@z31fJ+QmZx7ONg zuf6v3oSdB8>KrjNCML#EM?0gvA<&x^TPm2g6{aNBup2g`mvN5K!{`7OWvP|WK@u5| z$qJqYsXb-DsdFGlIu)3SGTLSm2lFD$f7mdXnBRyq<${-NwarILg~3?xSuW-6Kputu zqmVIehCwFESq9T^;Nd8nX$UZrlI46EC`_$`e_dc3X)f$DDf!iqL%r3Z?DsA9=Jzcw z>{IM5t(cMM_a!EqdEuI<+ptk%MPB;JJl&9JyE_tNL;2+weD(RXjuj_0EFF~d<*xT$ zw|u<pt}4i`M6x4wM#^8=?cBR=ubg8TmA7TPPoa$`>+O|>5pQ&#TR6x@V{_I`f3kCz zi)}+bwLEPxVk!r8-h6pG`~0!pPs_2_k3spJ7g%%VKrt>esq-mG{n3b)c%#ymWpuAR zw_WP;%K4Q!?Jaq+G51=w-d%b7a=U$Tr*`?BPOY?DJ&C@IKpKgZgLE|#512J_k@Aqn zAzg#SW4uYk-)r&rI*BKOy5;#~NvD7oAr(tZQlg#nv`pd~LCfX2U(yQD8AvlFCYgma z8|hX89Dgy1Q4Kl|={AWKnU9kiq}wH<7IXp9ok$Cj?n1g7X$jI&Bp&yf#Q3YEBo81x zi1arkbW7tAq(_k+LwW*94k4Tw##7>N=+DzQeHQ6Cq?JhjK>8<=I-Zv&t3h8vdRgLE zK-VI@iu4-N>qzU6)+231dIM<_5|7PDZzF9*dI#xUr1y|`G$3t5`T%J=(uYXu_(+|A zW>)`t_aon=%}aSd``605mu&35G&ak9bM0kWB^wi7Jf-rHYn#3~<ZgM~KXC2k@sHg# z=C1zpX6=6QgK6izRQLVKcbu9~Rg--C)jLjZKgV~=EfcqOyr=XpgIb>c;{Ap}yWZLR z;Ms3K^5xRv$-8Dx=(6p>PTpfPOX4P9^}&K)Mz8noDR1vQ_Vq<i_#b|DUc+zKzL`)p z+jw)ce_5Y@ozdIljW7Rw*m=(w&o4RkqnMS;V;(r?>{RdbQ(s8V-}TKiODC^hRyySN zIq)g%-1|=J^=|wf@9aJ0qeTBd-%WpRZ}Peq7SDXQ<K^!*ZhNF@*P^sOy__G8xMAP! z(%*i)py1w?32!Xhak>4ppT?a$sAKlH-0vQF>HEJWRK0Rz+~jwsEb|z-v$l<W=(Z1T z_F1xSDPQ&W&Kt+~c)RS%pZ5J@Uc)7cPxz1Q^6ctaa`fqgnO8Nu_}GtgOFMn|Lf0Y3 zn!5yY*5xd+{oL{8#Eb5IFaEQEe=A*Z-qfp4zNq)X>i2*6VeYxck!NqY@7<>#t=pdb z-2;<X_j`PBMtstq6mQPJ+dDmY+hucy9eVwhAItw?du+;(XRq2@|H|G)v#l47?|R*? zSN9+M;JcIUaT7k;@$<~xdk^*W4|aZdsQB3rm$+jO7j^0W#e+}m{&|PDCa1Z}rT3oX z__*k%lH@1+>5ZOUkNs)SkgJ}2;l*)9lXqS9&Q%Rd=AAWUneDwZrrvu0>>C@7HK)F` z*=PUla3|k|<qux&{&306>}jhnoOI>HmU*kbuR3z|$%zk5Uj6jfzq-CW*Pc?c=10eU z$ypgwJ7w*7r|FvgKcyBlYh6g39_<w3bi--WN#VH47LMQN1S${eH%7>-Bl2(PsI<*_ z_uQCp{DLJMpMt4T9khS0ScFUt@ydw&zem_VC89pdBkX@4f%~e$K<!A`5%PaVlz$-? z0(sD%Pb1=y8G*M%__KhEkUW@wM}++)5%GLE!rp=i{RblCnGxlk91))fBHAZbIRk1( zs)>--N7&Dci2vsi`aL4zcX5RLvWWOy8ByO4BI45!QJ)7R>a#t<etZPxTA&WbdyWF2 z#O)&doe@#~I}zo59ML~3BkHp&0^c6tPfA4orU?C?BjT|?!k-r+^nZ_NkDDUubAE*X zGa~wrEuy?TBJAB6kv}6s?vH4ncOvXv6M=7u&|e*4zw}hY_~-|VjbS^_=weJb+b|~m zXzp2wuQ=B*5`i63zB&#k%@!N?8I(T*d-Z8n8}|*T7-tx{Jq%+F<aC4bFVNw4d~LCD z??>$Fg?-AAkp0j9Y#2w|TWt*EDMnW#;VkTH(iJ|#7=-d(fPd%7{9hs->zgb#_9fbT z8RhM63E8VY4RdNcYtY|&A#YfReKwA6#z{s8WK2YPY!>G4cLDTaF9CLmA4h#R9k$rc zl=3gn$2iTf+6Dl#e%z~g_J``X>vY5TKQUHY7pcDn87tvG%nAKEl%Iq0CIEMX{qv#t zjTZ0wfj7%|Q2rCzqaN{7?eTFu=KOx4@@o-@BPjnIS>9)-7{;OBEw)N&e<0|)2Q4;) zQPig!;=l7)NPZOMO@e<62<y`a@vYb5J-8?KPsq=1&v>N3|4CZ?KR($o($Sv1rMy>H z>>o9MtKryk*t1LfO=#cgh_CYZG2lsvm-4SY;x}IN?|Jyw6LK_*DE|@oGaL577C002 zHE?Z)zapN?Fka^&9x7gmh|d{_C*wf>x1;@LYxVma{kt>zGuT4@D*D3`*jN2475;pI z_E-Mh|0k?Dn*F5^q-*tg8+bg*qeCpO9`RWGy~XyP>`%wge&aO#DwKctyO6()7@tWS zRejLc;mFsBud4r#sQ+^K10w7XM*muc_>7hKZ|rDqg!1RZpPj#1Y<Gjr{*;G-w;ko9 zTZ#O3&`Gj9*7-{KHv!{CjfXYpFQ;q%K9Bl48!fggp-X$;0CziNv7ws6-WEIJrS<1B z#P2BZ3MsFI{0#JeRelWYuf==;F)n4lZ%JeP;BO@`?eB%+My<WRM16~ZnSuRvE9<X~ zmm|RQ5iiw0IwqifQJ);>vpw!be_aZHP>jF>QT|c1ziRJq&xIWJRQ*$+zfz0;A>@1A zu-eAU{ETY}`h(K{r5&zmfK#OWTJ)zvEgnnJ(YC`MRsY5dF`sDT`3pEwrujDp^79bi zH`<&2UWxH?2mDjxvFZ%8AM6EWdC#HzMOyzni1uE9{-MhI0PS@T>Vs^e{?lN0;f7Fu z;eCT~81D$9koSOMIq(>XOW^;`UqkJah5WIwcafC8fq{2So4>mtzKaKk+V^q<s!{9T zm!khSqCM4odKcoo8~I@l?SHZ{1oQU@_^;NB7&z9f&F?otuN?JN_4@_o>-}jN=Bwpe ze3xS4sngoyHk6l%^4NatZ^<wekMh74^}qQ7l!til$AU!p0L(XQWPi0Ae>Q%Be3I;c z>|fK7u}tfqPhemymvZJGhWyJhf9}ILXaBt!@w*4(tG6s~8!Eb6)9-5My^>J-#iRYU zVt(2z^%Kw@`C9+{5%pj7rp0y_#xcutqP_vxWBuq~5z24D{GsN@JgnEv+I;Z=#!DR9 z+bQ+?aDLX}w+S7uQR^=+KyQjR9-cxx?#FnVDDy9aAIr7=w+{XDL9Ktbpno@M^Isj} z^|dCS0e=p`A2r|2M10~=-e_t6QmsAW&>p{{y#pv;>32eV6l?3DfqPY1zH$$4QyH!+ zr@O{@rg;4UPkC0c+wb@Ijj}*dxyN1TDsYz;c!~{|Yf6c))a4Jj%L6W#;hK;eC=O&7 z<OUp$f+BahD^TwC2K?Ct&RjSCIP=F`m7AZP=gLcR<vCo=+ze(J)h};!{!~vvp!6o^ zR1@b0%DtsivK-Ews|vs=EY6vp8>U)V?97=CN$E{f1<ehV6}MfiUM4(9cja7_n__zV zA7zvNI#Tk+dCL7>Uujm6r(hb~^wRB&yxf4V+&#tPaG3QOR9x&U2(Qu*hjY3Jv9~n9 z5F^@SvPa9aQ69w`j+l!Upw%(%(n4QJ3W^cY2{ME&l$)A2Jb9=Yt6Yc6TUdzy^yjMF zK>Boq)mNTfi0~IWovD8BEDs!ju$9a^!!!?TjxcNVaag%~vUHNCDTNA@xz%Z2xkq#w z)+kiD9G^enFDUnxshW*+qiJ~)lxA~^X8OGa{$NeV32#|bw32A=3`CkWNKuu&D#&q+ zH*2DCu5HDc5Xngy>7HWpf`(9Bp^R9Nmomf`P)?M1$Arp2gy}h|Q4|~|?7?OusWxgI zdPk~jq^sOrI>qBEV5d>t!px^?*Q!;L-2rH#z&tVRWZVjiae}+Z(vgvJyaJ9J;t7a= zNzLrhV|+f<l^jU~qbmYs6&RmWJb~;|3|x$3XFqgZj=0<;5p`X$-w=;KQ0|*)2CY9W zWS6!wkn5*2V?6GH0EUHZ(5q5L<&4Ns%;c0Yo&q0AEUzdCjPVq^1K#QA=H3ECUv$?n zZAXT))KflX<_LFwo-eP+6L4pZl!{?`h`8v_Va#L)3N6g?l}`3fsVGN5zS5!IVvp1G zQX2O8RmT&4I+8;CWEgc9SGWV7>{4&Q>n=7$s++JZT^&B;`n_H;2}H`&fFD`m4<rw9 z2i&q@1P<^2j--)ebMul1m6y9`s%q&xN5(KuAkS0m@fY~YJfnQ&C2&ksK1?eeTBg54 zJY_z=H_%!qITbn>ny8zvboAuxQcUn>%Rv|IIH<I=qS#%Y>(Cq?=JuDl%L_G*j19tE zC=97^Ix_NFi9~&1St}~b=PNJtmLkX^BZOHahC6fmBqt^%rY0uKIQavG1(#d`keH56 zk{d`)a=EYs`EGLgDgrLwWOGcJ%LEn!F}3;A-DPE-(n7g<jLP++c+4n%VZo8%8k3KO z$>VVqdJ4)tC7#ki4i>BsyVO_capk-H9?>mx{rw7jC8)c{Id!U8P?&be2u6JyR&oQW z@F(gD<s5@tDJacb;KF2y#kIs-Ae~sJqyuRMSYgZAj;?|kZdbWSL_pLgwE%ITnpj|i z^9XCqt`X*uw^UWPZBIw#`ilKxJwvSf2P5Sz^?R_wV<kfA`R-!Q8u@KR4U@<ibN0^l zyK`0Gxg7O3tLgP&M*z<vwg}D?k8g6YlKI$plwz9^C@2#3DE4PZ$r2@a1$o$5q+yBY z&cSSZoTK@<12IF=PxYaE>}1?rrNuDH4P?MI(ekP6!D<kq+*IlLWZ6Q|-2V^N(+Y}x zezh9=Wy@jv6g4*H_9;qj#a_yp+s}mx=4a=ox=?j)tx!?!nq-eb=L?1sOMiGZgPrEz ztEN)NcnS^-Ik;f@--=W!=C4wYPUxq&vJvyqV3c6ki6;B^jP&#RUBZspcI>t}L&xO$ za#?xTlo>NzW!R+oO5MfYz)aWlBw?C;;P^uWF10Ze<Ob)J;9SxVizq(T^+fB=t(pvr zoj$tJWZ^A_8L-5h==|PN!(HGLYwDB&41Tg5hHGr;O<uI*7>~cA#KVbN?#3xjA}){< zW%g)H7?@mL6@CvES21f~qMhLl81UFr%6>t0xp6OXqd=FZ5X%z;`F_922dh*p&5ks+ ztr+Pp#Rj+B(O=H+*iQy`1Y*Kr?^3fycCH^2Uzyumo{i3uGYV_tp#H%fiqnZuC%JO5 zpcMoh{y;W*E?B8p-!Rnuyka27>2P4@gbjy_)3TFuPKqi(u!I&HH`==6;7XUBrPjJ@ z&7C@@KZIYd0$#3u+$A5^S*PblPG#nte8MZ4;mUP!%Yh0WXKrrQ4nvYv&{Im>Wr$j< zAV&l=)Ev3Vh^;eCOtYqujB$g;WDgoV!kIkY!TDoUtKuB}vWh)!EJR!~&At2oviuY- z(_#RJN}s^e&F)#?uW(KFmWrJf)}v6vi;bk~-)y<Evj3%*4o9)?CQrF5-&X;5$NaZC zfg_pM2;M10ffDQlOrOKH;^@n~%#xDGL75fo3@3U7ro%6}1E?tV-iTg`k&q^KCuzgR z49XdvQ8JU4td1d~-<<H#k&e+Z%`?*+8=+3X2_ZQ8Pk5o?icN3#Hm06jCnj=S6}nuU zKh?~M(UUx8(5NA!M><Ch%O2$f5|{DLQfzBH|D7;JOd@Rrs`aQ&4UKAU@6GPjA8Tl7 z1y+t3j+9c*P3AU9PLr)>;Lx<vY68L}^WW>3TpEf!leO`3!n;9O-wO^bvGN435&p<> z651^Wjh>kGvL6NuV<MKypzUyn!$xhgIcl-)1hr3$Ct<gnj<XL?JIWL8fHK^w%PYYN zFHyDXAKGDHUik0rFtjo2Ex}#S-2Sc+`AOJbrHIL0+~x@09c&e_<kZ0x-r~YsX7!j$ zKW{!@AE1`j|9nGLQlj#n;C?4&Zd~eO$7^mHgR9jaYL>EK9Z$5MNuhPM&3zx27q!@i z?!NI>P*xGP?IMJZbn_0GySM=NZFo~Bcq7xi5J$a(`|5&XaY>H-pLwzMM|V}|>ToRj zMjHEe8b5(+NN}f%tAf^70LNc{R9@V(8|54`?Aj5db8{W(W#ztruQbnHjw_;Ta9Jfh zj&==^>=4%w*`)>ELUi~+rBgTryym4un9Zc&o|IhNNG%&wEUaC10t1jy;-2Bs@8AFD z1J2wBPcCy6-PG3glDIRNo#ia?1_GW!QB`&GRc^vh^caZPT5}`CPbW@vF=-fTh31`X zb53YuH)odq2d}2o?0>@D&w%&1NX6}<<G3BXf=V+#VGwQH>LNkjFApjwUu$0MwKbHF zJG5{?E-v1}MnQ@>vcd)e#%$Q7{&B6EGZe2cLRW#f2o$$D!%7N{w_LB>QUAe)3+)Gv zvrRE4yj0OaL*X0M#<g;CQur+r^J5nouE<t#I0h9Ks(UQ>Frn2hy3O&L(t2(Z12%Z0 z;{*?+WNx`JW4cE1Cb_tG%8P<%-2jte#O%sjI>&AKP&=Gpw|NB{ywz~Rho`tj#h{2= zRANLNZ{|NvA9Jd;D%inNdEF>`_X!RPuNTGbhBh}J|DCxwbYU8{6IG4>M?=?7&hMBb z|BH9&WWX??1zRKh+6ohCP;z3fqZ}#hTM=VT4Wu@9V6N!5T%Gy6{VlFC#$1)p=KZ75 zRu}ch%dcH?nD=UqJNTU$vWobiE$1q8U5WBhV#F=o@#ajY!V&dGGpo!`rQ}l|*k1h6 zO|QCUK_macr6vP6vRwQ?v8|;bLtJS^yT<whTVK^m!RaJ<RI<6z)%U__`d)aLd)9wv zy3u+=*tq|%?V3V2m;Z2UrJJ`s`B6o*4-ERbIF|a!J?kIr5qR=O-BNCS-QD(+4F8L< z+dudmhL|hGbRees|7m^BazwlOO64aBFyh7s*P`p3_$e@A)U?qY$DLl#&5pBqKH-zL z{>jP-H;&F2aYvpXogHUV!v#4S$h>7AwuMd(?$Z8i<DM5U+Q&>MaKUgS6?%eeMh>3N zn`#Ujkv%xeg=a{NVOd$O6jI4?B&HZy<HrvgobBSNV2>XUl~g<`iSG!v$8)6Z@Xv}A zI>-qfCNeqdoU{W}%zSuU71)$nV9}=KXFf{!tf{G?ED@_@=5A-i$bX?T#rf}0RSv(C z9Ak70mjA~}Gp(q+;Ye|hM@?_=4dtl+OnXn`AzgmG#*P%{_2QBCjz)bi036JBhbGVI zITNUZ@s=jPS(BUjty=ze#w${u1bhB-ym%tnW^6-#mTBs@57VzhH=vBi2hUAIGh2*L zrGDU|MK-qZ0a<==s@2v3&*t*GK-N+C&INEsJX>B#iq0^X6rE{4DW0ufKnjZsNm0&X zQaqcyloXa9Accp^N%0KwqonYE1u32}e})thT1kp$epiv=S>iRMaC$8%{UOCO`5Q^` ztoLS8JX5@t6wjJBkm5TF+ez`P^-fYe!@QgHNj&>Zif5wtk>c6-CQ>|;euxy$9ygP2 z#dkVL@g;+!q<F^t7%9K>f_QYaBXXpFLBvVV!!z`xc$Pbkv<BaeC&jbR@uYaBx(8`b z^aWBpGo3(+Sob2m$S``7UX1UKkgmjc`AIJ^j11CwhB1H?51D6@_QCi6NfYr-0MfpA z#-20@-%TQQ;JfUk$@ne?X^LT7Pl}mj5@{N~(?FVz?_!Yl!!y>T8F==cv_GD~Cf$PX zZjoM!_9wjz?N5q_{pXTij`k-Vi1sJF0_{(liS{QQg!U&LjP@tZLi>{rLHm<B(f*`E z(f*{v(Eg;u(f*{_hB3D0%iO9%mL^2rxNgpd#pQ-kd-t)HmX_*m0c+FQ=#fqLq1!dB zLRsUl+p?i0Er6D4xnKtVGEd6|vvF3$(*upj)^Y)Ra1mt(8U$ufF5<8_P%kihbP;-t zabTsu?Ab+j;^hJ_B<@ALP~gSHh`e#2QsAY;h@^3#OyCEIGl?e&%n?vDoH$2dj)bBd z;!J@#B8tWnCkf1vQ8bCzF7QfXR_H*Sz^jPMhz)_)5YHfP{vC`9))H3|Hww%VR5YKs zL12!gqJ_lu0&_$aEhSzlFh^F=a^mFzHxRENUMTQ(;+4de0`DYVLtG~CZsL04Ndj{O z7HuZZ5tt*fsDU_BV2;S5oy183b7U4Z659oCCT=2*6ZkN3GqEA?QR1V-&A+k!$A~f8 zqW%IK=vhVlKHq@`fi1*w#PtH(h~tS@3LHypCtfaa9C0t=g#vdcP9m-pIG&ha^E^-{ za1Y{4;z<JAiH8&C2%JEiL!2paFXHjUNdosKo<wXHIElE3I8NYH;xb}G;0)p!#LXNh z7Yrb-ByJQqlXyOHgTO<G7ZTSCJe+tb@k)V55HBZQE^rR<3gU$V=Mk?Yt`vAY@fzYX zfv+d7C!QqmB;w7)IRY0FHxOqETtvK+I7#4Q;znY-z-7cu#Bl-#h?|KGfoBjOC2sze z^`A}5{mX$yfh&n^#0>(^C5|Jm7kEB#Jn>3_7ZBTtmkYd*xEJw4ffo}e5mySllsJR9 zOyCEIGl?e&yqtJAagM-`66X+S3cP}NJaLl1&k#={whO$HxQIAT;8nzB#0IeAn;U9A zuiDp~lQ*V_--{~ZmsPJGSH$lWHMRWG(sFoV*dgcF%F9RKOHW3{pxV=NaI)qByoRRF zQv6wWv|Zi!mj0jk?RBd$ri`lfF*QP^rTY@fP5u6#lz+9!N)s_v8!gvO+)#6*ZY@?h zF@IMbihpQ1RIBR)-8e58IO+7JC4sz_hAkTw){RK8z~z=3kHMR3#^qi;=g@p?Eo$Db zy37(td~-;G#b&YEkfR<|ZF&o}B#+;Es<9+AU5rux3e6r!jmLMxs;|W;Apc>|+AK>| zy`|QuTOzo53AP&KtgFtGI(Hy%D91{HCP-Pn%5P(S!42l;$YmL$kmKOlH6PS$3R#&7 zUQKh&CzRg-U+b!sMb3A1)z`_q?Uemk+85lQ{Si{<ThN1Nn)a2I&eHyY3$eWknfn2% zb=4xHsN(g~#3et6OsJyYqzo<|Sc{!fUA0nWe+(hkUzM#iMyZU~AtM{%PBctiH77P! z_dVG(4+w6N<h~@#3LW8@;(spq3O$d$c`)9rs^a|(Jo<SsE>uP>^<+p?xi?9@0~@ge zL_f(uZ)o}_DqgdpW+xrnaTMnh7s`egUX)4iPJ-y*X*D})-mK}7(A0;Sj@JC^z${p* zs}@R7uXQe@nZD4KOrSxSZj>j=S|9L4WG^*6L60<=dTE}-VGma`adBO>=t9Difb?X} zPr4_M&}DkkjNNTrwNMhCET$(`n2}7N7ZFM2iL&+xc*2vF$p3NE$BlcIWul;uzumJe zL((Kkdr4}SG+xp;No}BP?6sm<X8_TKFHkFT-xAzP$=wQWUA52=l~eq~i~({D#U8sh zZqE7vn4oJrBY&K6i|9(RAJWbe`c{K+K+B(igFG$o*HtSsb!;u!OB8P=BPDx@a?V3N z*-I4fYN>bNL0mV4>N8(@-Yu*?3nlj$GSyWJ9Z?^}PmuMQj%$bL^*Jc&!x*djG>7V= z%)H7rm-SJ+C&5$oQO+%)o~)1JvA>8Geizq6q53S9o)3oS<eKfXTylGh`Uo9SAH~mO zaL_)F<C-XXef}lt(@U$*PoerKGv6Eu)kpC@22a&TIrj$jWPKFxA5!nY@3_tj)n|qD zyd2dH)n}#TUJq_vwa^juQT(|K4(hW5*Lu<G^RuXrU8~P;q53E@XKD3OyiQtulyg7+ z7^;ureGZ;zpTw~GtdX8S3^Sqn)JyIHaO<jtj;N2~uVSI7Plq<^(<15<uhpl$H8_5h znUVBU_BX}rFa6{!tDL()tB>NHqUv))oZdd0rRUoq4%MeYa@T@eS1oizeH4EmjiEm2 zm?5LLj}_k6Ey>jCb8@IY%FHc?L-kQSFL<iIDd)ygPxd#(8zl7(+=2THq5ig0dfqOq zK8=!l2$|}tg^s9?;&+$znT$m_dVT&R>SNRDb8e_U%FMGrgzBSs_kpMCqnxXzo~)1J z-6-`AyoUQHq53pQ&y(SKsC}9xw}+^Y&=K`f{2>eu+UFiD^wH~ck*H6cR-cSeeUzC! z{|ePd@wR}c>Z6={fqJq&iuaJzJMb;;3)NM(NRQt_H6<Eit`kmrQhij)Uluh|Xf}<Y z9_z6l*Ek;tsZ~KujRd9aW2L#03<wqXSMX|`57Zjg?bu`2#vz2U+--u|u!f01=fDBz zpq!R%sus%K>!i&Cy-<t)$<Rd2(Ep_C({!=tS;d`38P<BW=NZ$C%IvJEwjeI1t<7N- zsH@Io#AFP^su0``efvE$xCxntI!8A`^U8?sJ~kQ?m-EjG{#nUCYxY=|pJcV|nOI*} z-I2Rkw97G+Ad&Q>LbN8Y<taae_N#R+u3I7)c?lLWTTsnfH0S2JYNfY|9Y{`Winm;L zpo3?G%G(Owo6fCl!2@$}AJW`#DsTQSa~$BU*V>U=Uw*cBS5eL8pT1uAOKi>Nbw|&s z`QRdKDZTY^xVWi#>!S5F4b}A({i|--Y6Mbhom+Wzp{`Q~e`}7nn^7M+-R(OQtLv6> z)Z7f|M(%~sfwtDHU-!#7s)NxMcMi}&_Llb%1)8}HHw)!y6YeM0IyaM>hQbinTG84Z z<*XnDUtn*5hG;D~AOvij>l3=~ajpP~-Q!#ig6kujjQ(rrLv_`?fojKYhBG4i{ZSTf z$<<XiVzXPfMCiiLx+T<8O5MMaN`eVT_JW}rNb%kSPmH%-p(0<QUW(Z)EP$1GLg~ro zQtsU;`3JUkLL^$%y(g<^U4bt5=@xXs)(RJdDXiJ!T)79=05j20dzLLk*?XL8$gmwE zW2s;Wjh_yJQB4L6FIyhc94NA`2x-0p&E0i^Q7tSo*wrEj<1K={7V50!)yP~|Es&Us zdAUgbeTcUiY<3Ck`Xzsk;0vuV{y+p{&%}n1ViEjTp;wvnz*Cd3DkzhBauQa&-cs*C zmyW!QU}q6sCRY4~H%p~AcSm@+R`PeDjii?FQt_AC=;cmSSSqeZsFatA_jmA=m#U!q zsVBWuym?aZz?2SLDjH!@>#%De4R%;Fh8ElzVr(xL-cdiC{|jonXJTVrwNgxAsd9c% zyf|4ZH%`I?3?BPB<elO*fyY`L*n&INb=AF4w&vZ}kP7d}r%xYY`TqFZkGFhs;IqaD zyWYFO%t6WDpgMO~*SSAqBH!bzH>cB}V*T9uu5)j}T(U>Z?3xK7DWJsJ$QImVUCDV^ zjQDzXLd;owNlepg2hg!6)(6$rv{K`{c$)5wP^)jHCajeywVh_6n_H>z%|T7s5o!&s z)a)$uJX1~dLCjRbZ+Trl44qUC^q==EvtdQrQ!V~3TZysNdfbmhgm*tn*`8(G!0d50 zXu*(1w>43m0kIr>Xn+K);d_?Z!+1ReFCM%DZV}?B?Q)(H?L_VLCJS5`jCLDgyoioy zjnal%JEeB~Xz$uDtx2^uqP@ygix%yV(F-G@9W8!~Xdxng`>?D=#&5UaF@6UV_E;O) z_#&4&XSCHhwQ!Rwd;mR(g2VtKpTp!HSp4`~4xb|%on04W04d20P~`wR@ePGnLN!%> zB||F*S(JvN+O<lp?S^81PmWNF(onQ^lBpK0p^`Xa%=J$Vi)amX;LFwx#n+5RHdKn> z{re5YF}f^?t%p*BZuA0!8x-?qQkiVBC1@CO&6^c7B==Tualeg7tR0H~1o&*|pSkSg zW)eT$sO<($%QoD5l9}@n7&T)n-YD?YdZLQHRMS(uo=WeG(48~6G^5yf1WZ|0i?;<l zWm#=on)hjIfn4NwgD2L4Kf~uL=hhf*`!>{UTKD@o=4Ozeq}0aAP2aYcfQ<e-W_Rho z<5o-m4dvFCxUrkNrEB#K^76C}+e@|Od=AEI$8J@7u=>4bQ16pF(1j2)$Z7}Lj9Lbp z!dc0Pc9AU64$k>x{eqg+LCO{|IlZ9N^n+jvjiapRX=2M}D(sh)JlKU9VPl7;AsBBc zjnl-fdwJs;nulOp#a7(~-wB?IZ4SB*EpF8G6mO-br|$4SKs|YfUh(EjJ=!^h7}e~o z^E6j&G_Jkw2G_(o_puH57^3P>My)fJm9BM)?fXcvBiJdC2mWJ!KouCd259%l#@)B< z<aY5TnA|h5d5^Ps3oYRTjX8%Fpw+5wvBh+)z7IQ;Ik&{Nv=~O$>IDF70poNl?_1Bs zi}spc{v54hNvP&|B7(ezuakrlUUZL^hB75A)j=Pun{8N68uy$vA6!cLa|g)IcXl?t z4vjT@JD{GO5@!Dwdxs@q1oje58SsPGjeD$3=HeT_$J%V3i3{E)d>lg?Z?m(ZkS}GA zdJP5C?0lpiCSj_WroO=*$24`+G|*d2O$o?muS68~SS3GR@P&EuY5tHj-@b$C;&+Xr zGst`MZZ;JzM(KRo7ouRkS+dUgB00}P%2Rn2*%>D4oGU4L5>BE^Nh@b+R=#7XaSIx? zWBIWhZ1jWGWoFbA)e@VOJh&q5Wz(MQ?xKJUMuE2tMoV9&fU$;eNu)8=Y&1i6X2JD3 z=j%Lek*Ct<c;WLJN`8cwa+7OTue&UR!i?j(y+*n%HciVuV`0MNBAS$L3&tEqLb@&7 zSqVlfw`qojZa|@R&Nq2FT)6$F3Uwdhc0DBnHMfPyiLwj|=OBd9qK#z*(XN+ni=E?F zpRzDvau*6yZVSd+U?{hRJ3oWb%59oqp|Qg4_jr1jJeAST5N>a#<TjRJUQ(Dfd|Z}6 z;adn{+iq``ZZGGO-B-9ROeWBzjJ9B$DowJJ2zLg9(aLR_VWDG1cl(H^pU6|`c9w9v zfszJPNVzRcvOhC46rL*GmVG*+Z#PJ{8|n6pJ&d+6SwfT2ZNa!+nxxypomwzjxlJ=H z^g+?xKI7^6!tKvgC&?3Tix1=e%ra!Og-P~j<@WUmA)_7L?VZwXzQMU^-EO)qOs+;@ z%5A}T3=Gvtgge{8XyrD|(AWK<-S+deSkz#@a=Tc#-AGA3%aCphlkCsR?Z*(pXl|Q* zyHUE$7cVz;7H$iZEhtR6Eg0W|A>6*M&MDkEE1KIhLywZoSXr$6d;wEFt<Q#qHO3p^ zQx7G~&;l=26@>ZMWO)>R0l(&K<U++|&|CxKQ6~$RDwt=BH1)>EUSd_l$|rP`qFfix zf<b(uk7W;xJ(<u6o_E*IduZpeFa^y|w}>3w1?v9kmQU9S+(Y0V#KAnZnHKHbW}f3I z1E|)HwN<Sb%XUrogr;h^yvH#8(d6-L4$3{z<-?WS8C{;I<QGJj7c057m7Ee8sv8hc zc&u@-e8O2XcC|+B*w}Snw2x_si16EJL$gNX@L*hKowE^&;y&S6o7(vGfoWb5IrkwR zin%DrWRb|nMt=m%x%CwTBpeD@v!Tv;P|L@e1IB8u=dDQLF67L6p~xHP$Y*_wgVr_p z3tNF3-*g`4X`S;h3JT>p#5}{GTj$(fH}(+RLTfcfcs|lBtk$``#@Sfs6nTF8J=ALj z;}95X;Pw85lYpp84GP869(A6h*vnmaLqy#A!}4X$*sk~BW=J3w{dpBO8%?i%tfGlv zcCVdCcWcMS(JMxYj|B8EU8>nxJMq|-^%nU=qSRGN*%3;z*Vn?KIIVxiQJ+ynDD220 zT02y-vUYwO)^J&Xhei6Y@4C2NtDV{Y;@0fIk8lCtjN|o#ttQ7}_7m7oK+Wo$qI5P3 zO8*MM5_RovHpF4Q${#b5XmV9$Rhvvp&3ahIt@rmYTF3k0?ajKeC)#T(L|emlG44Ls zY>?H=1u!EThc3kaqq;;*YC_XUNLf;Iu+_{6qm+r}nK^FFxQR(@0_wP~>6~EQWJ{ps zRrmrde`nRk?$=2htf&ZyYzz9E!^W;)L>Z2ku;c5$3325;w6$hz$S;8@JgzjX*4drb z<r+Ae1#FU8y8|F1P(F=k7>ohgc`(DYM^{CEiyR~UWZs=9mLr{wz}AVye2e>{V45S1 z{qG2L5TXOWVtT-GXDZvYMqh?vp-_vZu>FA#p@kkmcjh4@y_BARLEaJ;81CYv-{@l3 zxyk5s%;{_hg{2#=3((=rzyva-*T*og)h0>s9YL@Pf$$cF`~9T*1B7GJm%a3*lkiV& zJXQ%aWWPtna2H=aCZbqUYi?2=gKl-bXskdV@Qxq9)v(9-y~6s9MT5uHj_h7D0`bNS zgqR1rBb>cm>fQ-mF?+IM^jJI0+SEGtA=X|IF<GT=>3cfzGvp%v9a>l9FdoWbtU2~Y ztha_@y_LLRtW^(LEx5<8;dbb-h8(taxSaq;>Vln;_eQi{@ExVBUPOo|2+smjr0;`~ zUxoN$hU~aqG}m<C;B~n9fcqJ3$M<5Xm|W*<Zrzb(J!+lJHSG^}r6+O<egqZcMY^2W zF);*ZB)R%yzQ6=v&PH-RJvc(v<w3L=YFj%p9<6Rp>lm3{JTq3+9D~7mI`y2`c!*Y8 z#9rx)J9x6hbtbO5CC=!&1k*#slX66HY;t=SKGUEaOucfR5nkb3+)wC50fIk_2E|<R z6Cse?p<qA>Rp`}b%H5KhSy_C=P&wpq^PuH2XZJulJ%Nd3z?^TzzNXG8E|8xAk4=nK zTSO5xz6pD`$hl~~^9Pwv+Fr#T>^uxqH{4P?_HftL&L0pqM2}-8xS6ni$b}{re%9sS z7FKLc4&rxMl!RIR2eCj`ZM0(*HfQnRjB8G}1F7IVB+DEiQ1d3|Pca9tuf?K|&S5V7 zdu8X9E7u3C?KWsfblzo0TUzEjf2zFYCnL~BW_=4;51u3eA2&vMRks|Hjd&kBnP>zt zTn59Y$bBm|Cp0cU>KtXx2ViQ-K9q(yW^Z!tim4rmetf9Txf|mc^L`vd6)|O`!!|}| zopTp#U=`jiqJ+-#Cm28QP-x7a&-uda^mR@#Tkprz%&{wcd=Vg`-jCrSCwMY+lXGv- zlYMP_@(e5pPxhLg>{Fi9F;-`6avsF=anO`)Rk9lyx@cpss$TBy#R|<$0UGC>ZFJIB z(+-_DxRvE)fEWVUp-a6F-;dC%I;6%2H)<@r&MC&nUd;U=OBn8Qe=Lo{qS#`AYmOCj zkl~WW?8b-d#>ds6lZW=fd}>2P0Zh?W7$aiybR7gxWi?L6Df)A8#@S%r+!V53XqvZ8 zM1|MBgNm3vAFeQOIr27DH=yP4Sa8$zJdu~=mG8B*z)v|piQyym)0|1_`M75Io`vfU z?rL+z>be?#ZS!;F45vr@x_2WYsYS%pQ+o<t%fUC{_IIG6tB5%4g=#*Pqp=%wMB`&1 zpUm@}^efIB%A9vitE_5!h)pZxnNyz}IjgoZF0r*a;&O*uQRKv}(5gedYb&@#;+-{a zb()^zGbXtAh;5bwJ5I6r7VYsErZC*P{ip_uE%X_sTJg-pzM5ZHx&t+@+SprMt%m$! zZCOjRafP9LQQ0uh3&7&lz@qxDtLw4Qac0N&iB6{1=oJky^)U_W8r!p&sc2!<0kwWi z5gr1O_1E_JVn6vM44IKx^VCyMbr$)cvoPc*^Iw6_FGBG%yU!(2c<rO`zC-ZLyy6CP zG3s3B+(@h5w}Us&Zoug#XQQ#n`8m=TNPD66Hnpk{Jv4z4i>shvsAApNw>fc{BX2D= zZbWyL%J+!V2FBKg40U7Isu}VJxFJT<Rx{v`H7HqMxGnGHn`MrM84iX6`?2boxmcd; zcro<Coa?@VbLz&3nEXBjLkv26q9E=|ZKQ_zdA!MZUTPfh!ggKtqg+!^iLGymekke# z4N(KMny6+_=_a&RD9g3*1HGY`<K;rs7EuykJh&RcBu{;O@(<_{J2jm?(9wF{D6>e* z1gOx3lOs%&!o-0=nD=R?P_y$T-43lo9i*MVNIRP#)?MjIJ7ZyoT>z#+m04hx7+XW* za%eCoVxaoV$1o=MdR?nmp&vmWTPyXs!->#;Qs`H7!`7hD+(AWL-gKQ?3M=S$G0b_7 zR;`mG{4Pb{siQ)A1F<C`6@@naYl}HAFmCg}$2{cR7pmB)FoiB6@9I}jnGdREpk(aM zaYe%0L1u5>P7~*EX73W~ICt;(tndRuNM^#g>Pmu9#M0PJ55f=GIaBEwyFXar5xx|T z)H%f+!M+GKPoc+JtwPP^4ay6Uoy9^z4uo>GNMLU@Z(@Fr+^W-yo0yGYs7^2JccrKD z?yO+^t~7?2fw3;6d#KYJ5Qw(vT73Ykekg3a5e8ZGlOj~2U!IMwBSR1z6zez1($5U# z{yTVePBA-VY2FG(vhoE(w4Oekt7!;EoYFX5&QvW7o#(K)U4vVW<{~TRF!p~|op!R} zvN>;N;Ie_7zPwl*i%9NMQ_UsdU_xPtoy`X?ka{K5%e2(SRjrRxpAVZB(v$_>_YjMO z|692^61!Mw;1`S&Z24D;`5zw%O7?TSx(X{Fzf~Qj{BNjdmft<ReEw{!un;W&hi{^k zKMQJ6%I~Xkhw2Ybl=@T8EZ-JhewwlnEZ-lk{9n*Jqm=&~*b(&yr*-+?P|qx%I~no$ zycz!&l!ajVKOBlyf2ak^7yCQ8D}3WkcEq2;N1Esxub3UYZfvt^+4sNZTpCw{eW1{& zfrc7If-wUOG0wV4>ra4TUcJaZDR!cFNj3Js7v6~MlYhk!3XjRX;Ixj($J8@p5+5Fu zhm?h2Og=jpB_@-h7K}-#PoAlA2m1uCFIty>E%nUudxV#NuCfp;-xaO=uP|+B<;U}5 z?TVrDBkhXIF?}!ON7sH#6|FwHri9^}j_`4Yl%IRw#SH7$%Y7F1l9cHmJ1oLhP>w-f zL(s$}nCT+-PPK5YZ65m%x6#b~>UU`J&@4PFG!?#%qRj2CXnQqRf~)4Emo|uPyxh!U z=ZsBj@QQKmIw<P9*jJ!^@C>m|iAnt=POI`$Y#(Ii_FC*>%}eFzb_*lzW^D|&`#eIb z?Ea~hUCy7V%9{w9VTj#m=qiM9;O5$m2kX&aTHW<U%*?yaOTz9ti*=?1UZB^@>pfXY z<XvwO@scR*^`yq@bud#dwP*4n>`>c$iqK2R5nd{Yn-iLRDva<7%o$sCORT)Se<89L z@Ua=zKNItwJmX4k^6+fCJi`}T<QbkY7ia7&TvSEx)9-fU1_`e$e`4+B62hm)5kO8N z@{)2sa>_Yn1B9$bsAAdtB7?f8rB5_V834KXjX=AHUAbfnXH<JesZX7}<l>Deb*_H~ zA+ArJ{8C+?h^vy&SAtR2#0_<Uo<`{VyXcTjqfzRd4X5G75><zKG);L8p1%^0Ep}ZM z|6cUwx|Vj;ErA5b5%UIyx~FQ6o|eVWQ>$uIO#k1@Kdsu>LtL`e_DH}55z0ePv^=yB zUk0u^+7dYR6`oCASas0S(oowTtP1Q8Eci-sP`3%b@Ldq|ldwx1fUbTm-ur-s0uR1m z29w83UxmHY<0-TUe0G1Cr@%XTraj*mD6+fl`0Xw6D`57?_{}Xl-d+fvdbg#?^H&r= z(Qm)hJ_tY0g5QAi+dGxaRQVF|=4R7C$xM|mu~7c5nZ3NCG|vas!bF@IBL6H;A%7OZ z1Tfauffs|0@lNxUx|E#nwKwqgYQ7#>YD<KdB^%BeeTU#rfj1yto^1Ylmf>`gYvkpD z3^K;!?|9={{Jqx5EAp1+dU0x=;H6jcCAoqneyhyDkF@xTLu?%yc>NZBnT$V4CN0Pk z<r_mozo=zQ&dzdOY>3|lGwkAh+X$$L8A5LiLi7XfV!Lm$-LQN3Hg0<n!kX{#l<I2c zj{shJDc-}R{Q{lAUu^JADR-BY*d0Ucm)U0^;n$=}147~$qAMgTg_(B2P99>U4EMN3 z&cx4zB?e{$1R9b&#XZaA$af{DUln8yDh{6Jdc~>vvN<iUWVFgI`O*TG5d5!1|GO14 z*6+9j)p9F9C4hP8NSw2)Dxs!Bb|q&gB_o6Pi5u_;+8R9Dx*XqEn~(3Xmf^c*e2<X9 zNsSrA^o23Nw;70|`HJW?A=BhvYtqa+!)WE)#6gsM-g!q*i635kM#*4f^Qb8$>z>Yf z<e?1@E#B};*0CqQ-!?Ko=ZZVNd13ian^q3Yes1mk|N3awz*D~NGv~g0Z@m1=8~zk$ z?|;K(J8${q$A*8sc4_>ioY#)5Jm=ENq8iVR{+0c+NBp|zwq=bOD!;PF^6@sF|GNLG z993@rWn&&5`q^Er%2jfuA5AYJUWS3Vs(ed!^lUmQ@sxqbt7ohHsQ-_jJZ|~q(T)Rh z{%QYh(GwqLsq*V&yp?^ym+h#x-_wH@?B2ENp27cVzPk0_KDpuX?#}bRb<{p~LB$YN zep~#1Y`bj6U)*alFa75GxEnHm8l-UFgzG-pci;L!Emz)}{M*j;gOq&Tw5i_y?_D!! zrObCAC2rKHpW-u>pQ>J}p4-2_z3JHLPYk$2`gMMUU*{d2^umL^)3d_rtHodS%RX5> z4t|k-#uWqh-f~T+In{V`5sy$mQ|%<!m4h?8$$or~Y`;H^+yCHm?Y9khCTrXK6U&qF zCL|uA{L0^(yZ?RI=d<Qs5m}G;K11?;J-1V)N&hVA&sTivd~M5sSI^J9c<=`o%=`Sa zA%pKNYkGEVdFH_bY1tXy|Cm|3_%A<<Pg->O`_JT$daT!!%eTDhANs-B2QNQ5W>&^| z$G*Bu`EmCvALXAhq4A1;?(BJ5pD}frYoDLKtNQx7%$r{RVdD>H)@2^+<XE=5fAtk_ z+b$~psr|D9j`b}+|IjNJ4YuUo-)tQ_^YVw=@4ayQ-dnC1c<+#3d|T#UepTO}U(B!h z<?_Hax4n6|sCeMD_ixW=|J^?>-T%)6U1x3^c-f3uSN`Se=RUf8GLLRIW~z3+DE+1U zO)pN%YD%e^J@xxX2ddvF7k*~mS+B;A*&bJ3Va99m6(+v)iZy3fAG~_a9daBVl{7ii ztmn1$Chab9cZnxuj>^d!J$i&;^c6n>>hG%-XslVGAKUC}^gZ5T*u%8VJ%M2__POzs zWmB+}Im^p^<-P5tKCDNYS)pjyC-YaBj0?<FQT$qFusFL|DeZ!3S1!oq6Z1I!SE6FK zV&eaQwSKT&V$N)TStZsNAIh$Lxus<|ZjPMt0X`7}y$-YqYtCtS&L$QcWe4bT==KL) z3px;#Pce)EJq&sSXdKp~si3_<r-Kd$eG0S?bTw!t=z7oxK)(Q81KMs^OUri9^FW(H z2ZP3<{Cv;{K<@xu1G*e^JLoFVX3(9xTUxr~DTQ`>T3V7o6G2CS=7AQ0=7Y`!oejDi z^kLAopsPW5g6;-A40_(D@EhyvK+xWx3qePK-Va*Dbr*Cl=%=8|LB9oE3%UwVaqb4K z2R#b96|_6n;^BMY59m115uki{nQtAN3py8c0qAnjrJ!p;C+&kB(4?<gT8@DZ0FB4} zbADMQ6|@jE2Xqc-G3X<p^FjX(`Y7m2p!J}9E`B#?F}@;U;N~sAh0_D{@1Pl=Z-eH6 z?f@+V{RMOZXz!+$mS;dSK{ta=2Hgj`5!8nE{|Yn#Gy&f!$pq~UdOc_+=nT-Ipo>8# zfUX2B1l<Z+3_1k;rV?}#=wi^>pesO^g02N!1^Os%zwHD~z=io?V$e8zvnrt(?E(5I z=y1>sTmlw@t^u7-jBBOC*bglLjl=$D1!ynY1$`9q2GDxYeW1HRQ*hCRpQSfyK@&jt zfer`#7ib}9Jie}73ECU<0njT!*MN=(-40p|+6-C=8jEswfc66IiLYY}0qqSs33NH= zT+mgZ+^>BFx)$^o(4C;Mc%Q=%_}c??5-493!y{%^&KM(RwylHsmG(1ZEtdH)l-ZGL zcAzaljX@{d23NQ1g6=Q!aS^J3j5ty@9;_Q6!*;QiECALx$k0r!WDCI>+J@{I$kN-8 zZHBBTWO(IZE1Ubk8rViY&V7k($oK@_*^r$R#U}6fpVx*g4>I}^EngXA7q=l>02zIW z%!gOkhy!gjLNroG(MX2bN_b}@qzOn#FSNAWg+{ZiIm>Dr5*uS=$64D|7(=?YuPEy< zJg$Aaia=b4b`>+?I<~8r4Smt}*s|B)Q=BVVLPY!afb6voF%F{0_@c3QA;UMBTjk4x zY$IfmZEuI_R;0C%@v&aZXuH*>%DlG9!UxCb(_-Xbi#S|_^GKgoLiQG9VLlD3YR`OI zA^QZf3-o+a)&$u;$S@RJ*|ebC9gyJ_G_7QI$hNm(a{y$YLMHCnsCvKzeHkxphWR38 z0mvGW&j}o?pY(GfWS#M(AztqjlHrtn;29|!qsxYe;<gp?0mzH#Gb}It;dR`_kewAp zX2H7G8?v6dO!{MoEFCgT9bq=<&j841^9|Sz)<g5>ddM5F&PE&OGa##PL$(;Qm5{NI z2JK3JRzensH8@&3ZiS5Nd9?DHAp0IR&(m$<oUyTB0y_#>v^mZW*}f=#iuyyg9Wu^6 zK|f{vDZd`_+3@E*^z}$N%bg8*Z|no2`MVUd9*{+AFFvz)RvY;mAd77y-yz7lKo%_q znD~uT+pw7c8TUES{LF;xd+dcGWp>n$`>RICR^Z%Hq>c5fLVdOp`3qsYCv<`_l5PC} zWTlWrjs?bUjg&=;&vwXsZTQ;^*$t4Lua}2&wpT3nQIjFV-NukFT7TyLjNe&_oY&~j zaLCFa>#66(xgDt<ai0a*M4Vd|XmQA{vd**ICg!gB$iE6Y(c0!w$X<p_Tr@_^Gxd;h zFBz@xa(}oPvS|Kt&sYx`Tnp6|r>qP2kx#Xek9$k*Wuwh8+;1+7V$%+r8xi+uZ??3& z3~E`b#UM8ng9XSRgC`Ey7Qw!NQ~I<5vKiQ8Mytz4$n0(S+Xz{F8#3;9<JypM4{VDf zv%{B#h@Ao1e)Mn4q;o@S*yvDQ#&b`M?-x}<HL@-N$a-PF9!<6ovIT9(o`I}6ip-88 z_)JqC_Vjn--142~lQz~4(I@CcE7yhw?S*TBXtD&z;&9_WTI@3+vqh2F;l~>IWI#3= z#9}!w6esBub$IP@33SF&Cr;CmV}m*d{B|iFdutu`4_-Sgh7PvKt>Vq=iFuG+pv&YO zV_@Q12wAi>ya#0SA&WN0WI%R1WYN~4JjiB47A;<7kOkV1Er6^Tva!0KI2ZjNvd1CI zi6Yww*(Z=ij%7P^2cYg7A>*FglC8x~>lc<-jQKXkmL0NK$in7C*rHDZAUm&(eB&WI z1F}e;m@fd?0QBQLoJWq0g^&dxi<|>Q`#?Qrme>KusJ$=kQMW!y{mIWjS_Rnz&|q1k z%yCLqOH0fmYahI%Yo@h6#&{{ly1bq7aSXh^ubp*oJ7ae{>!Kv%w)WNs+8ayS^L(@T z`)PaY^A_Vsduy%5m}}wrGZsqmGD-aX++wY=8pp)>3aj;1tFhWj{%7LvL8}#S#GluJ zoPTz(zS+U37w0c@B<Jmp){kw*-i`zRVKbhx@qC+&{OV5p``1p^Cp#O@cN&OyNc|$t zAMQ;4hvM(S&eoq#GOA7*`2I=8TPN}SUnf!eSS<g(5^LSk#rPn0;6q)E`@8UbeHZe7 z6n|?^w%*s(c;w`P-=A!JEzTEprSuK)cUM>IfjHxbt^?P_88652{EIkBmz=`CPoH9a zshhFkl!143Gj8w3^M7_D|4Z@r*KXFD?#A7x4*dL7<C9Z){+I5QzSG_M-f70J?$%$s z8;8aD+ox5<K=9$|)<u6Z?)XzXoPQwx{x#luEZ$fdZ+$D?*dWe7{PPy7{{2kr(>;t2 z&a@sp)7W>W_3j?Vf*#gS`xqbfSQP{E;MvyIeT>y-Tkkl>c=v4US7#faiSygeA%COz zyZapLD!Xy`9P0wRQElh>vvx{96MuiPTeqESEIQZv)VanZ=komRb1D7p+}k^YEbVQb z+s9bhyT|U{#)snkx84*#*GG8KhxNQKaXS7k>uYW5YdqJt^QOMW>*D<LzLeG`Szk^v zo=>uVk!0*jLdW}+f6y&Y>+v`28+XK5zmGA#W{<qB9R<d@F@HK4Cl^@*F@vqH_BAn` z?~XAZp#s!zk$u{hWbJ3<P`?F_lqS!_g-VQdV>{z}>b=^|dV720BbuvfZ+*>TENhQb zbnH*s8y^V9ix%sTR%1Ii#(2HsT~_0ED;eKft<Q8cYCC|vx`TCd2V-3aG9K$_-D@+} z^Jf^o>1b`~XdD%c-8Sp}osFMuU_aW)y1J9`TqiQ_?rcT(T-h1yU7fApbT+;c4D{E< zU5x!FfxRf!`e>}NES8Ms*eh5i*5me=ZV;Yj-P*~RV0}2oxYdeX&tGG#zsDGF#aMU8 z7$1uBpJK@WTRZ-Js@<ygTWtnj8-R)Es5K4eFL!9SrGv4)gB9odJDh{Rk9V~G(9w9k zqxJoc##`e29$PB--`Hs4sLlFNC*%H3EOu=tIPeAkz@1ZiS!WoRTCX;)1igcC1=Dtx zH5CkT9vcsOku}M1M#pV5a7#?=naI$`daZHg#b(`&vtmwVo3p(a#8~f&F&=8Wy<1wc z_gifBEf%i6@QKf}m<Rd?`-n(|IN`ZZ(wUN0NqQ=-mnr8HdZd>~%5NGG^WjlazQ%!c zqNGzL<@1H)^NU}k4@>&2q`ZeqKCf9xFU4Mpv{BNrs4lTb(stN45}z(<FG;VLbf%>I ztr6<+MS7(CtrJrICJHHky97&&(Xh^9<J!F&uPOiW_ZHuRITm(Qqw#JFFNKUr^DNvF z89(+ixJ3{jablGa2V6GF@OFPx)tG@lT$}J!XagUM8slHL*myo1l=nt(i_PFZfQ|8c zON(#*YZe=t#K*F5Ytgr3t=Ue1|IzU%o>ov9PG4`a^+ZZX8j3UyX)@Auq<Kj9AU%$> z8fg>K4y3P;jv#fwTRKih>WP$&G!$tZ(qyFRNb`{HL3$i%HPR-e9Y|jx9YN}V`}n6L z^+ZZX8j3UyX)@Auq<Kj9AU%$>8fg>K4y3P;jv#e-1No79BBdh@MH+`R8EHDwJfwS& z9!FY@v<Ybk(pN}FkUDHaex#mA=}1G7#vx5cnvOIN=^muVkyayZLfV1!719wT5kveW z(xvuZ!$yt0$ex;*nwVryPD<*Zl$Mlk?==R`6AX6;On#r#iwypL(jQiKn3nq`r-fT~ z^hs`Q`TrB;<3}UWKeHs@&NmYMfpU`;sWWeBz*AmkB$oODo<#h-RG)x*ibPXND-!eh z!|8pzg+?L~p3f*U5({UR!k$S3<)&o1r`*qnQG}WVF38F~#nd-X%8CO<qWCS`#DE8n zT_jEhfXG+K2V4?8MXt$spu^)TDnw3&V6vdV<(X07DGRvZJRUnqEWnSp`r&{{rWQaQ zg}6(+1<>}vF#aTp(oD<werV%o3rp~ji5idp4m}e}oR2t%q+;$#1rXN@GBag7&PU?q zCiR9xZ#V>YoT~%Vr8u}hS9(RzD}r8vl#wXCAviyRd588m=kcC1KmX!<#DjWFqmVf7 zQIB)sBut8xm@qjvB8$-B+%^u0^C9&(cQ&FDoHNxyT^=_eao(dI=e!*QaE}i{RWA8V z(~vk%@|<(!eBAV?$2rTTE`8%M3kjEnW?V4W4#zfu{wRCY<KV7BQhJ=5$3u^EHubm; zD7^(Zr>`uBYeT~z+)u~3I#jv0fuZya%%8?!FjS#9r5@8Fq@W&t1H%}GwFY%o2YuzS z6iMlEE!qu-dEKH8RqlgeP>*pr3caJy8^DUmLoEi611r7ME>>G=ong#n1Ia_#TLrB2 z*2LkLHx`1y1yLSK?=|3{9zIrBvkdW3Ldeu<gVbYD^RYI~f7omfKyfI&kAP^8di))a zrO;EpQ&H)CigTsMAM98TJtfoi_5!fYXpd{+qtM%ojO6iP9@;;Ir0f-A?Jj<_wLR*7 zujz4(FM}S(oI0q(^b3-*w+25CvgWbY_Cy{iv3&Z+Z}Y5$-gb68dDtbjLQgBkpx$K~ zL!EMOL0jth2@TF~%c(=1lD1P!(8g!*8(n8(KZkSbvkvTo>_WP|q_+*Dus@jcQ1%qo zDd#O>zWR#?=O@zkvSiLfCN@jpd}CskZisos#O46v{9$6V_0-;|y_ng|`M{J5HL=Gs zu{jVpzD>-Xn>mh6%$W=a$19HZV&vn5<5FO=<8%C(SlmSt<IKe3(+`8=$;9U9<+w31 zTLt&KFg{EyJ{1$=z{KXv&i-y<^KJ(FwTaE0rP|ZAH}LXv9PFp2oI6W%&uL=b^+0*( zPbU7e!0cZp7N2AZo{7&4>suD%`xbUVj#J8Gu<|g6@-GHE+FroQzZl~*BMmr7#;YAB zimQOd@9u!cAom1N{(=u%mbeIXp{B3*BbLv_a<Nwa9N<+lelgfh-wCXj|5xCZvVB<o zKY-Oid#?iUHFhkA<#Rsb3+ojA44AK~V}&U{26|_s#m4*5wAUMoe5IW!?`jMLo`sh% znfbelIpAgu^VI{hfzA4MHTY)FN)7XyXmd4uGw^&3^D6@jG<+ZMLJdCwyja8he%w+G z*8@MG;jO^SHT((iqZ-~1yh6i=fmZ^n`dTqrtkLjkz_=|HlwSb6S<BxKc&mnaJ=maO zUKeiH@HN0YHS7l7t>IGOMh(vd-lyTYz)c!n1bj%t%Yd7KQ7yD@SK~S0?=`$0_+J{{ z0sMo84+0+sR_*ya@KFt)gemNphR^1d)-j}?2yD^tKw$hZWiWo@fa5gm1CH15d|<nV zmjn0GFkcjtq~UGA85-UPoT=fTfQM_iGb)y&VLR}64QBvP(r^xNk%kL^%QSp5@C*&# z1zf4&$ARZ-n4bkM)bJMIr5gSWc)5mu1YV(G-Y#0H;d6l3XgCqLUc;9IZ`SY#;06uz zv)P>*E(UJY@Xf$Y8m<Ly*6?3}k81cyU_+Zf_{^A1!*2q|X?O>4yoSF5wrltZa4!v8 zF;69FI374d!w%p~4G#exuHgy5IU1e{JYK^$15eWM?Z8DEz8|<u!~X=Hq2V`xD>eKn z@O%yb1iVngU2(CpRKw>2FV}EC;1wF?Yfx8e_<G<q8lDPVui=@%4H~{3c&COR1a8zY zpQ&ro@DAW+4Ic(Rs^M;&I6H;LH$NA(X?PTHoQ5X>$7{F@*skGQfqQBAPT(XB{|z`p z!_NU{YWQ{F;Tql!oTK3{fyZn3SKvt+?u;!|k%syCXPJgC1fHSc4B$!)Uj;m0!v(+# zHS7mos^MDT<r;nvc!h>v1YW7(2H-UsZUU~?@NdAIHQWtbzy=MU54=;ue5rJ!hMmAo z8lDK;tl<jaqZ+OSHadse{~=(ThMxzH)9@R>@f!XJ*sfvR2RC|Ym>;kwX}BA<uo)UY zA2?IP{eg#T_$uHW4fC1F@fyAnc#?*111{3=GT<@|{{wi2hBp9LYWM@-`5N92yimj3 zaxc|zJn(W2UkJQH!~KC*YIp?j8Vz3$T(999fj4Wo8n{8je+AyD;g!IR8r}ljq~Tq_ z%^E%ed{o1+ykx+OhB1r`?l&$5wrO}MaGZt<fa5hh9oVkndBD9id=GGvhMxw`(C{0; znHt^?JY2(_aS4;7;l9A*HOyyECuw*RaFK>*1D9!d3GfUJ{}Z@U!&`ypYj{8KLJfDu zB@{oG;Jpp*w|fFFg&fPgh~G8lB^D5`O}}s9{+-|2HTQF87~J3U>$(ba|6bf;vE_l% z9{1}n<GlxJ|2G^b_`oxWxnExktjgnly<WrIuW!^a_v@QA%>DXS4RgQVpkeOUw`-XD z^_?272i~n=?$;YN%>DX44RgQVq+#yYn}Jn-<9_|9hPhuiaDA%e+^^e!ReiZ%kJT{u z>v0<9e!aVfxnGahF!$>{G|c_FUBle3Cuo@a^<En0e!aJbr=q?28I|52k~BGg8|^7| zz1Oab@ffhWUNzdCYOIdHTY!1L22n6#+Ql260juljn0BWdO~CwdSySHK_z@Tn{z=@` zXopMsoyRP;O0c`3KP3Vm{n27O5BNOb6~L!AS#0M>ybXBA*A^SE<S9QDw*bER!NNz@ z1->1azkSQ+acS>P;LS%YHr@#!UIDBQ-YwV%+%v{%y8!33caeDuL8$O_Xap_<o`mad zR^Sv|&sIdp?~K4N0DCb0s89Q^0vG=ts_%a4%j<Re_iaS}-y(2#?PfJamjJIsJe2>X z5%T+h%QSzV0#^51y5Yh2mB9SnU~GFueCmOhBA#bSf7>C@>9F5h;;R{Wd9eIK;M<!+ z@!tacHR`L{V=wS>E&j)V7k_WD@!3qq!;4!me4RYU3*&i<`mhg}_U{9J0sg4|aI$%e z1`<==lNjaC2BE#DfYtq(IOEU8i@<yx@;>ZSXn!a0Jk7taf!E6WH*uK%jsUkD4*6e# zTR?}j_F4>_DC5QQ_5&Ac{qHLC77^VO{^v*F2P5$F5qK}~64Xc4cRX$}#iM^>S{L<O z4E&K+UMJi@;#=`m{LTVijQN`7vp(km*K6%_E3i83IJpCOJH{89UD$sfINq?@MoS!v z4@`J}YXLChcRoKrkq6u7YT(^k`%VSMN993#8zS`gM_^7wQ&8S`*kU}d03Lz-=q3VJ za-xz4%X<~Lsh!p4z&Z6VLq{F2%^$_=xbmR>eBdw8o~l1@#K}>_Zv|w`-xoJ#{}0N; zGA?j2@K(fYv&7#5f2#GDj+}_(LHlP$;NHN_MvLuA$T*%y0N;c0!?3fw$-w-D2AW|0 z=@I!Kj==kYyB)IFIBhcj5W8VKsQEJvSRJ(IwJT)WzZ1AoD{m$61jLWyo%Y(Fi}9hg z*A>9}co+)Ix4jR5KHJ|5tPYm9I0A2vz(;_`Y4i7m=-87qdqddK<w1Mb1K$IEHQ!7J zz61RW+j8OW`<lKTyf1<IgBds#@$7s)&&-4N&H$c-`l<FG2mBP`mjW5<Ujod}b<FY8 z)d)o7U&Q><K8t=PLcSJwx>mnt;PQhO+Zbv8ngrxWyaN({2;7YISB<yRIdRK_{`LYM zhx{BS^v?;r3jU9k_!ZzvjQ1E|w#R!B`M;sQHXb|lgrYoXzc2!?1}@Uphc_eSp99ZF z|5fAT><dvJ%%`e93<O?-@xlh6zcYaO6Z8{++5hJQFW1`pF<^DjpO=7lYw`F6c#+n> zV|$_hY3-W;tPbYSB*yay<Pl#D9Kim8)u2DKfM;m)V=eIb-$Lu}8enxW|2E*Akl&4S z`fs~PNlp9b0jF#2Jr<asTT-hVKDoUPcs$}y2Frg7n7;&~;`KGKI#^z(ixE%wtNKGS z@Zs-5>(4@9bx{8y;LgxzH>7>*C5G|3Hh*3U96<S2DW3(b4(4AFfgcB64*%?su|K^8 zT&JxsNf?OhQD5d|dD$G8@*r<w1paFT-Uw`gf2w|IeN=8Uo;LxXf&QZ8^}z3<|Eux( z1MnA`Ki<Uf{8bTnQ3QS%I1%wu>qV!&;rWLGZ^HVZ#@ijhOSSmD3w#FfKxAY6FG&jR zFYMS~qyry8f914I`6OVz-HxR&J}ZFLLHrKz*V=yV7vM6~Pt|t-9F2vX-GufEf%7%} zXIT+>&|ZB6{uH<b{;2lfA0a;$fiFl#fB!W!9&&)UW4==Du>e>dEbpxd{3-BkZG3$P zY}e+8Q&P}hVJ`=HS)acEtAqC52R?-MR^zpID)ueFPROXg5V#5bkKbD(-VWT1_TYB_ zI6m8_;XAu%--%K`7g!y%Ul@TafhQrqs^4bdjVP}(bZP&hbkqm)Sr>_K240T-q5Rtk zd`w#}_5-Vf<+bmpklFrc0GFdZRCyOi$gc!G59>R>%f<3{M#xXeKs>eiGYPm!pC5qL z!SW9P&(i95rS_RJ<wikyz#piXJUOwzaJjPbMvrlg$j;4kxqya*p)0dAQK8Q@rP!D6 zE_R70bX@L=8F*_<Nm(&ogHxE8)<4-{P~h?w&LB!L%*TETD@sacs&bs8hKOQ@ju|x4 z=~9TrsPj-M;&C8Rih3hYVgVeUFv5{CWt{n7QdW_tU>aDtfe{%lm)9qrLUK(Qm6kPP zxHD&jBh6FnDF~GNN?ip-?(!)fCv;OCE|;G#)hTdM-Rmy)&JwTKapAQ<u0jtUp{fwC z`f(+Bz0g<2(mZ9w6@FS6kv1)F%vf=jK5gir5xL^D|Fn_r8Dq0h4plOUG~g{&z>y>r zbBlcC0cwhTI>LK=luTF|n~letpxBT3Lg(I6d2SXw!dE&aSg;1Xr3yG4qF}=GRH~Yu zgkw`(%9E!J_6Gbpp7Oyn10Lo`n>H$EbS}@*r;RQ3dxT%<({jr3yqPJ_E}iV-S;n+$ z+~uWsb%LoZ-kxBd4#umL*k-7B?o7WX7>dWmiYag;Ps_&RYeRkIH@V9T1t&#t2D=NU z1zD-nauLBn{sK=aU(q1c(gbgaCybdgZHRA_FOXC2EAZfnG2vXcKdabRD(d7&nKl@& zASx&toj(<wMc7G_tU;ydxiU~?LAmM?j-(QAaj~}q{lA=DAjo%w_{nVvQo;mDs6`-} zr4T>DQb-VHDXo<^>4FoM+hm31>p#s~>}!=5Jw&m?3~2l?E9qtKa=)h_ylQxUue`PC z{-y5H*4zwt=@bk=Ux{uaZL+7l#H({s^U=9<Mt>3UC|ndB$<1KNt(axrf>zAr(l9@g zE5cz(7)&WC35yYMSZtCD!qp2R)C<DZ9bwwZ;mXP3y2;bRlvAdKYXgUSQWn-i>E(FG zi90Nu4zmZh7Ni$?rr@<^VWIb<`?s=?0+z0oS`xuX_LuAO6d2aABZ;kqzT_(nvy<MM zC7q2fK>J4#B(b$y^V5nv<$+nc8)8_r=9YM-xLa|POZ|FQM;J;j34^H~S8leh-4Aa$ z!vM_Jc@B)zDQ>ldP4^WC!lcP<O42!M+m=r89>;1rTeI!FWq6ZZThg?ts6$xaNEzuF zr~8~74pXK@;9(IsF9PR;;|$i5Z;R8*NR4Eeb3-IId0KdX3q#2zVKDV3cOWeG>7|~V z0>1DelT_p>bCtRM1@2<qofOP8I(8(v10_Cx8M;qcsM6gfWf3gdv7-ynv7-pYsA<j5 zU`LAVKq7jP%+yFmimxavylIo&5hGbHgHieLvK3#<{H?hhlC8LDH+k~I7uU3a$6xH$ zBbqkFUEvQKU1HR>;$n1{me164rcd(}Pxpp3V`_;nybJfw_Y~Y%ftCqt)wF_&a&K7s zvD<|iSDiDA8SIOZTm&l8K)(WPNIemCEv~@R?s^?kiYsR5xF32^7B&+)&*1<wugF*7 zhs7{K5(Fc;i3#t_$z@^7ZCc<aUpObd$Xz@+oRw1M4O{;*22J-(aTk{R^aA@8xQe|~ zio(3>U*PiJSm7@Bgt?Yc<jN_=wrrZ7-Qg(0aEc;GBR_0abEF7XG}HY((YTIuHLtl! z+<w2SAgnfyRQcw|XtK0YPldS3h$6(eu86Re0+wE@l#+@h9cOsT{FvBE!<Y_C_a!c@ z5@p5M?ds~uN%+DAswA#kD#|@jD)om&IN4E9j4I)Z$E+e=xmoBnuVvJpH8~|%`z&13 zVRCa7drGGSiu95kj={d7@+i46sTG$Mxn1Q2#Su{xJ;zm2klb44nC2RcHEmW2-X0oO zN?LASjw{c@`OqCeg~Ei%6y}s;Cm$(EnKsPr&ntqc$X6WJmy#V7rCxr0#p5dRxJz-7 zSL*i$I0T1{7(IB<2$w6{h3$LBxItsG2Mr$KOdjuGnb`Q_suD546-%Lu?|+@{ag|MX zCMUJ#;u|ySwR0ogrS2)7a&fU5X5Wz(p^}yCoQ^AFZ+b?0;uKH7RaW2%6jhW?OU$1^ zHm*)x?s7~kE?nY-GjIXRt4C>m^w6QX&OBEheaA)okZVT`8kwDi9O7y|cxCPPxr#6W z6nl8B(#jb}|7>+R>>1<p1;oqr=--4<xq%GVn0!|NJHPB)e{LYP%<V1DzAD$BGwT1= zc4f<LBS~<+$iW79>CGCdN5hqdVi4S2+b<m;2+8n*1Xq9#@7HfubvJHcWCj@F2aAHh zKB}@Zvznm0T3=uN%g#$mFB2o>aV&{6TY3|Pm&P%+3x6p+pI?ze;JL4-mrNa#eJI!X zBiFtGzr<JZDCb)<Dm<QV54Dca7USD6TQ1fK_7FdRm8{b5aFbbnlZaRh%!p<0U}gAU zxw-SEShEXEa;`k}>c=mbAPp{GLY|@8`Dysu!H<run{mUDL4y*|I+VH5BN9CH56b|r z39>V=A3o&yGLZ5{$k)%GZ||4CFTQ+H?x(SykbCe?SoX#))pn%VPtQ?O%E>N#LRF;D zjH4aUnSu?V7#+=fDCg4%zYN}e#h0qU&w}rtKjY2uSu-U*x%B?^bXs}oLeP4r^`Y|C zA@rY<H>?X%Esm3*F`?tzk2Qy5GjwN*I0Vh2rwy;a<JrwYN~70FlGcWeMWt~xZod-l zwrnfm8V)QkJ98ObWUbC0B+308&xdAP$r1S99^qU(o06WoEIqO&N{jZ5Fo-`&XSN|o zU(C-G7fz8ejQ8q;A%8mAoIZ0}Ld3~i_^fvH$V<7=*>i|CcXKQsd3LTRE|XBgp{GR? zeE@;Z^fbhXu<GF@1Dec$3-r)!^P6_Gy+ERiDaIAS4nvlp(10ycluvjOCbygG^`G;@ zPkOSoxS_H_m<jJqKVSPkNhun=G}30{;7~gS3jHVisWGUC83NuMWyjk;oj)*&^<99y zi}ri`(+HA-)<C#F#!8jLrvv&usWu!%iazOS5*aHeoXz#>=HV~7x?Qd$enL5eS9Ani zX_im~q4+{!pzs;d;?D<P+#1ApEN>I9#syx{Q{Z*LUV&8{!4==%g$sSDQP32zMtv+W zyMgXNa7`>O5$aU9<M6SRfnE})GO9!RM=Mt=5?An&uApGNtEY{FU!c6=;rWIl=!}y4 z3k7c_lk&tURdWoIE1986A4Oz4!NH?*hSaCMVuIB%;0iuDgip*hvD||XXPid=V{BGW z2gnP}Domy5VU)yeE&@j8#vnsbRaql|Io^t@?QY16)FNvGRSq!#M2`Z58;!1!;};vY z4Hq%u2TSvFrNtDbPlx|acp_qKa>h$xEFR#>b4h9;UZ=Q?P=cp}ReHN6M8!yLh0-(` zJl5q6JBw@sqAl4*niTQ44}<5BKEqGNpm>=o6tm>U&I}%zC3KC@&s=&}>V7wp60cKX zH%-DWUg3;?r^hMayF1Vd)$|h<AaPN>38F>A*iAeKFVJ?lOhO2%&e#hesupHBi#cn* zPiFqZZ`X*oh51&ewjv5)N!V4cXXnThG-9>7`L>VOO&O5WSd+nGs=qPCj|2Lq7AVgA z3eP4(E^vdG1rBr+fx>L735vyna>_cSVy%XKctkH3ij@Td@2;1{P|QuGx4~<$lt~eS z18YIJg~pyJkhhjhIaekiK(T!7k{X(x3v#{=;v0_{0<AHE@{|Rw{ubUq8YHK1l{_CP zPe9V1-rp6>o;1O?$iSCTucNw_9FQC3(i&bw(|73SRr+s3e*6f5{z$Im=JxCTgXf8< zHVC(L&LI1Dq>WbSG?+zN*`mm(<WYa!P5WkjNb6667_F-o{R&~7gZ+;mZST;AM2^6w zEEmgT5WyD9u8&UyVv*Lql?mbtj9f2!j!xHPh5z&ThB+#HhY%cTZ2MvX-vO=waiFRG zcAw=Qb!0E4<7e^*N6diRI6wIB$})<7yka0vsK@ZW$>|JPO9ObKB;^`1CW=UDSPf7g zT|o1S_B768a8^?n>chQBGD;x8z~f1Mf8%vYr~c_gra)tFPR3(|;~#`^ZA7yzt@`mP zN0F9#!f3Kw$U~IjFb;2``JN|hV&sK&ASG9G0WKx^CxCiBlau_^HK{`}>cBhG(%@>0 ztPPJgK5R4q0&BsXv)Bu$t2ri%DuC?nnPqOPO#(Xo#86xmI3X~PC3q@DlJS%p+|;lS zFPq_Fv=hQ^4p^!;NH4B8%n79jvlRr<8!OZvVu;nuVjTDoo{dopj%+;;8ZVKs4|w4K z_Sd*$K6U3S8}xAUe2)um<I=9$)r;@&qM2D$V9(XKULYT9qpHi@dsoc~$PUtWD*1Rs zDk)?u@00_qr@L<@IQh$e(nCAjjkM@d+*Ewkx2)@GeI>iQ0ILpf*DxE`-nQ3g*JNVM zoT%8tmEOZfbUGuW=OGshbB~DvInK`^m(w~!Xxv7K>-|1&OtiS*#=!>LeUGgCL*5(r zy0uoi&05JYOCW2v6p=SD);eUmAXC+fBoUkAWDB*$DoS*ffFMO|)j>J3#jeyOAQsf4 zzPYua!FAwSk!hMc<g=|lUAK>=_V@>1L~c3JB67=#(%j2_N&I?qj^92C)4PtQqhqT% zaso6!V?Yt*Gn`#SC;6w|+dNmsq$3#V2VWPmFxgY@R8o5PYh^Hc*OV^T06NEMiNLD! zP(fB$aZiX`yxy(u%Z`4vn%SlEM!XY)FGZIPP+;7L3fe4671uONCW>r3H+W`-0)ZEB zvl8ks>#o^5WO~WEF@Zbya+M;deQP;z$Vno;QLWHT-fguquGLy$@bTi<IMTlgpOt|Z zv08b#J5()l;_|lEvubaV%N;EU_Gk)+Vn>T~Zdh2Gzwhr$K}{3+8?={mvp%7Of_(k0 z)_!SG9<^@#7uFkC7(jcgwTO;k!@)3&a;mC9!W^<H;hK@e7xe%yEj0iG%dM@;^E;>3 zcWY__HqLN{{0V1on3LNGK2wj((XhiQo`Y%CmNC{gH>l~d{$t}MUXyu~t);CgH1<TL zu}9HXHMYgvG$Q@oHnEzq!_yXXb~;-;PCD4}t7eLc-b$_&Uhv1uMPY+ebkGd*y6LR+ zaij{(o|KK6TVJ2|y#U%2hf&I&vKb@7@T#~#_l{BH!EW&Z8`5sO4mX@)*=tP{d%Ix- zesOx?Y{@azJ#jh{)CP|;?kdboEk83=z(dRtgk}q-|D;79X!hWVyvo}^T#XsWIq0k; z!Va7KPbYME21LDyYkx<(7oshp#&kRf@93j3qdnNc<F<x99edX+8L$97{>9)UE{-cL zAF$l`FpzplL+;AWVVxz<IE^QS))&`wvPf@_mlCgQ?2J_+$2WT!P3JD)RS()+j-1pE zlX%FL+G!FxLf7Nt!;Z(`^l=lQ=9M`wmf&O<cCUU^%ng&rmT6ZI(#d7y#o&>wuKbu8 rBm=q~Eu+<a$Zo~}SRtz0V?y<1`7cr<G1`YYJ6eNPbC&ZEo00T?ieQXq diff --git a/NPAnalysis/10He_Riken/include/ObjectManager.hh b/NPAnalysis/10He_Riken/include/ObjectManager.hh index 8d2739151..e86658f86 100644 --- a/NPAnalysis/10He_Riken/include/ObjectManager.hh +++ b/NPAnalysis/10He_Riken/include/ObjectManager.hh @@ -28,6 +28,7 @@ // NPL #include "TMust2Data.h" #include "TMust2Physics.h" +#include "TPlasticData.h" #include "NPReaction.h" #include "RootInput.h" #include "RootOutput.h" diff --git a/NPAnalysis/10He_Riken/src/Analysis.cc b/NPAnalysis/10He_Riken/src/Analysis.cc index b51e713f1..2dc14a435 100644 --- a/NPAnalysis/10He_Riken/src/Analysis.cc +++ b/NPAnalysis/10He_Riken/src/Analysis.cc @@ -31,8 +31,8 @@ int main(int argc,char** argv) myDetector -> ReadConfigurationFile(detectorfileName) ; // Attach more branch to the output - double ThinSi=-1 ;double Ex = 0 ; double EE = 0 ; double TT = 0 ; double X = 0 ; double Y = 0 ; int det ; double ResolThetaCM=0; - double ThetaCM=0; + double ThinSi=-1 ;double Ex = 0 ; double EE = 0 ; double TT = 0 ; double X = 0 ; double Y = 0 ; int det ; double ResolThetaCM=0; double FullTheta = 0 ; + double ThetaCM=0; double EPl=0; double DEPl=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") ; @@ -41,6 +41,10 @@ int main(int argc,char** argv) RootOutput::getInstance()->GetTree()->Branch("ThinSi_E",&ThinSi,"ThinSi/D") ; RootOutput::getInstance()->GetTree()->Branch("ThetaCM",&ThetaCM,"ThetaCM/D") ; RootOutput::getInstance()->GetTree()->Branch("ResolThetaCM",&ResolThetaCM,"ResolThetaCM/D") ; + RootOutput::getInstance()->GetTree()->Branch("FullTheta",&FullTheta,"FullTheta/D") ; + RootOutput::getInstance()->GetTree()->Branch("EPl",&EPl,"EPl/D") ; + RootOutput::getInstance()->GetTree()->Branch("DEPl",&DEPl,"DEPl/D") ; + // Get the formed Chained Tree and Treat it TChain* Chain = RootInput:: getInstance() -> GetChain() ; @@ -48,11 +52,14 @@ int main(int argc,char** argv) Chain->SetBranchStatus("ThinSiEnergy",true) ; Chain->SetBranchStatus("InitialConditions",true) ; Chain->SetBranchStatus("fIC_*",true) ; - - TInitialConditions* Init = new TInitialConditions(); + Chain->SetBranchStatus("fPlastic_*",true); + + TInitialConditions* Init = new TInitialConditions(); Chain->SetBranchAddress("ThinSiEnergy" ,&ThinSi ); Chain->SetBranchAddress("InitialConditions" ,&Init ); + TPlasticData* Plastic = new TPlasticData() ; + Chain->SetBranchAddress("Plastic",&Plastic ); double XTarget=0 ; double YTarget=0; double BeamTheta = 0 ; double BeamPhi = 0 ; double E=-1000; @@ -62,11 +69,26 @@ int main(int argc,char** argv) int i; for ( i = 0 ; i < Chain -> GetEntries() ; i ++ ) { - if( i%10000 == 0 && i!=0) cout << i << " Event annalysed " << endl ; + if( i%100000 == 0 && i!=0) {cout << i << " Event annalysed "<<endl; } Chain -> GetEntry(i); - myDetector -> ClearEventPhysics() ; - myDetector -> BuildPhysicalEvent() ; + // Plastic + if(Plastic->GetEnergySize()==2) + { + DEPl=Plastic->GetEnergy(0); + EPl=Plastic->GetEnergy(1); + } + + else + { + EPl =0 ; + DEPl =0 ; + } + + + /* myDetector -> ClearEventPhysics() ; + myDetector -> BuildPhysicalEvent() ; + FullTheta = Init->GetICEmittedAngleThetaLabWorldFrame(0); E = M2 -> GetEnergyDeposit(); XTarget = Init->GetICPositionX(0); YTarget = Init->GetICPositionY(0); @@ -126,9 +148,9 @@ int main(int argc,char** argv) { if(E>18)//CsI are inside a Mylar foil, plus rear alu strip { - // E= He3TargetWind.EvaluateInitialEnergy( E , // Energy of the detected particle - // 3*micrometer , // Target Thickness at 0 degree - // ThetaMM2Surface ); + E= He3TargetWind.EvaluateInitialEnergy( E , // Energy of the detected particle + 3*micrometer , // Target Thickness at 0 degree + ThetaMM2Surface ); E= He3StripAl.EvaluateInitialEnergy( E , // Energy of the detected particle 0.4*micrometer , // Target Thickness at 0 degree ThetaMM2Surface ); @@ -192,7 +214,7 @@ int main(int argc,char** argv) EE = E ; TT = Theta/deg ; - + */ RootOutput::getInstance()->GetTree()->Fill() ; ThinSi = -1 ; } diff --git a/NPAnalysis/10He_Riken/src/GNUmakefile b/NPAnalysis/10He_Riken/src/GNUmakefile index 15a9a7a10..bf3890ce8 100644 --- a/NPAnalysis/10He_Riken/src/GNUmakefile +++ b/NPAnalysis/10He_Riken/src/GNUmakefile @@ -18,7 +18,7 @@ CXXFLAGS += -I$(NPLIB)/include LDFLAGS = `root-config --libs` -lMathMore LDFLAGS+= -L$(NPLIB)/lib -lVDetector -lIORoot -lReaction -lEnergyLoss \ -lMust2Data -lMust2Physics \ - -lAnnularS1Data -lGaspardData \ + -lAnnularS1Data -lGaspardData -lPlasticData\ -lInitialConditions -lInteractionCoordinates LDFLAGS+= -L$(CLHEP_LIB_DIR) -l$(CLHEP_LIB) diff --git a/NPAnalysis/macros/GeometricalEfficiency.C b/NPAnalysis/macros/GeometricalEfficiency.C index 5bbc6f586..fd727bd4e 100644 --- a/NPAnalysis/macros/GeometricalEfficiency.C +++ b/NPAnalysis/macros/GeometricalEfficiency.C @@ -36,8 +36,10 @@ #include "TBranch.h" #include "TH1F.h" -#include "TInitialConditions.h" -#include "TInteractionCoordinates.h" +#include "/home/Adrien/Desktop/NPTool/NPTool.dev/NPLib/include/TInitialConditions.h" +#include "/home/Adrien/Desktop/NPTool/NPTool.dev/NPLib/include/TInteractionCoordinates.h" + +using namespace std ; void GeometricalEfficiency(const char * fname = "Efficiency_10000") { @@ -65,9 +67,9 @@ void GeometricalEfficiency(const char * fname = "Efficiency_10000") // Read the TTree Int_t nentries = tree->GetEntries(); - cout << "TTree contains " << nentries << " events" << endl; + // cout << "TTree contains " << nentries << " events" << endl; for (Int_t i = 0; i < nentries; i++) { - if (i%1000 == 0) cout << "Entry " << i << endl; + //if (i%1000 == 0) cout << "Entry " << i << endl; tree->GetEntry(i); // Fill histos hEmittTheta->Fill(initCond->GetICEmittedAngleThetaLabWorldFrame(0)); diff --git a/NPLib/GASPARD/TGaspardTrackerPhysics.h b/NPLib/GASPARD/TGaspardTrackerPhysics.h index e23742bee..5f8426c05 100644 --- a/NPLib/GASPARD/TGaspardTrackerPhysics.h +++ b/NPLib/GASPARD/TGaspardTrackerPhysics.h @@ -27,6 +27,7 @@ #include <vector> #include "TObject.h" #include "TGaspardTrackerData.h" +#include <cstdlib> using namespace std ; diff --git a/NPLib/Plastic/TPlasticData.cxx b/NPLib/Plastic/TPlasticData.cxx index 1221a2a4f..ac7ae2789 100644 --- a/NPLib/Plastic/TPlasticData.cxx +++ b/NPLib/Plastic/TPlasticData.cxx @@ -22,6 +22,7 @@ void TPlasticData::Clear() { fPlastic_Energy.clear(); fPlastic_Number.clear(); + fPlastic_Time.clear(); } diff --git a/NPLib/Plastic/TPlasticData.h b/NPLib/Plastic/TPlasticData.h index f11a97f57..15b2a5901 100644 --- a/NPLib/Plastic/TPlasticData.h +++ b/NPLib/Plastic/TPlasticData.h @@ -28,6 +28,15 @@ class TPlasticData : public TObject { double GetTime(int i) {return fPlastic_Time[i];} // (N) int GetPlasticNumber(int i) {return fPlastic_Number[i];} + + + double GetEnergySize() {return fPlastic_Energy.size();} + // (T) + double GetTimeSize() {return fPlastic_Time.size();} + // (N) + int GetPlasticNumberSize() {return fPlastic_Number.size();} + + ///////////////////// SETTERS //////////////////////// // (E) void SetEnergy(double E) {fPlastic_Energy.push_back(E);} diff --git a/NPSimulation/include/EventGeneratorTransfertToResonance.hh b/NPSimulation/include/EventGeneratorTransfertToResonance.hh index 51d85928d..762617531 100644 --- a/NPSimulation/include/EventGeneratorTransfertToResonance.hh +++ b/NPSimulation/include/EventGeneratorTransfertToResonance.hh @@ -132,10 +132,6 @@ class EventGeneratorTransfertToResonance : public VEventGenerator G4Event* anEvent , G4ParticleGun* particleGun); - // This method return a random Vector of dimension N and magnitude R - // The return distribution populate uniformely the surface of the N-Sphere of radius R - vector<double> PhaseSpaceUniformGenerator( int N , double R); - void SetEverything(string name1 , //Beam nuclei string name2 , //Target nuclei string name3 , //Product of reaction diff --git a/NPSimulation/src/EventGeneratorBeam.cc b/NPSimulation/src/EventGeneratorBeam.cc index 80550358d..10f13e75f 100644 --- a/NPSimulation/src/EventGeneratorBeam.cc +++ b/NPSimulation/src/EventGeneratorBeam.cc @@ -172,7 +172,7 @@ void EventGeneratorBeam::GenerateEvent(G4Event* anEvent, G4ParticleGun* particle G4double y0 = 1000 * cm ; G4double Beam_thetaX = 0 ; G4double Beam_phiY = 0 ; - + //shoot inside the target with correlated angle if (m_TargetRadius != 0) { @@ -188,7 +188,6 @@ void EventGeneratorBeam::GenerateEvent(G4Event* anEvent, G4ParticleGun* particle RandomGaussian2D( 0 , 0 , 0 , m_SigmaThetaX , x0 , Beam_thetaX ); RandomGaussian2D( 0 , 0 , 0 , m_SigmaPhiY , y0 , Beam_phiY ); } - m_InitConditions->SetICIncidentEmittanceTheta(Beam_thetaX / deg); m_InitConditions->SetICIncidentEmittancePhi(Beam_phiY / deg); diff --git a/NPSimulation/src/EventGeneratorTransfertToResonance.cc b/NPSimulation/src/EventGeneratorTransfertToResonance.cc index 500f81a12..f130ce20e 100644 --- a/NPSimulation/src/EventGeneratorTransfertToResonance.cc +++ b/NPSimulation/src/EventGeneratorTransfertToResonance.cc @@ -588,27 +588,17 @@ void EventGeneratorTransfertToResonance::ResonanceDecay( G4double EnergyHeavy G4double mn = neutron -> GetPDGMass() ; G4double mp = proton -> GetPDGMass() ; - G4double Q = M - md - mn * NumberOfNeutrons - mp * NumberOfProtons ; - - vector<G4double> DecayProductsMomentumCM ; - vector<G4double> DecayProductsMomentumXCM ; - vector<G4double> DecayProductsMomentumYCM ; - vector<G4double> DecayProductsMomentumZCM ; - vector<G4double> DecayProductsThetaCM ; - vector<G4double> DecayProductsPhiCM ; - // Initial Lab Momentum - G4double InitialE = sqrt(EnergyHeavy * EnergyHeavy + M * M) ; - G4double InitialMomentumX = EnergyHeavy * sin(ThetaHeavy) * cos(PhiHeavy) ; - G4double InitialMomentumY = EnergyHeavy * sin(ThetaHeavy) * sin(PhiHeavy) ; - G4double InitialMomentumZ = EnergyHeavy * cos(ThetaHeavy) ; + G4double InitialE = sqrt(EnergyHeavy * EnergyHeavy + M * M) ; + G4double InitialMomentumX = EnergyHeavy * sin(ThetaHeavy) * cos(PhiHeavy) ; + G4double InitialMomentumY = EnergyHeavy * sin(ThetaHeavy) * sin(PhiHeavy) ; + G4double InitialMomentumZ = EnergyHeavy * cos(ThetaHeavy) ; TLorentzVector Initial = TLorentzVector(InitialMomentumX/GeV, InitialMomentumY/GeV, InitialMomentumZ/GeV,InitialE/GeV); - + // Array of masses express in GeV/c2 double* masses = new double[NumberOfDecayProducts+1]; - - // Filling Array + // Filling Array masses[0] = md/GeV ; int ll = 1 ; @@ -617,21 +607,33 @@ void EventGeneratorTransfertToResonance::ResonanceDecay( G4double EnergyHeavy for(int i = 0 ; i < NumberOfProtons ; i++) {masses[ll] = mp/GeV ; ll++;} + // Instentiate a Phase Space Generator, with flat distrution - - TGenPhaseSpace TPhaseSpace ; - if( !TPhaseSpace.SetDecay(Initial, NumberOfDecayProducts+1, masses,"Fermi") ) cout << "Warning: Phase Space Decay forbiden by kinematic"; - TPhaseSpace.Generate() ; + if( !TPhaseSpace.SetDecay(Initial, NumberOfDecayProducts+1, masses) ) cout << "Warning: Phase Space Decay forbiden by kinematic, or more than 18 particles "<<endl; + + double MaxWt=TPhaseSpace.GetWtMax() ; + double Weight = 0 ; + double Rand = 1 ; + + while( Rand > Weight ) + { + Weight = TPhaseSpace.Generate() ; + Rand = CLHEP::RandFlat::shoot()*MaxWt ; + } + TLorentzVector* daugterLV = TPhaseSpace.GetDecay(0); + G4ThreeVector Momentum = G4ThreeVector( daugterLV->X()*GeV , daugterLV->Y()*GeV , daugterLV->Z()*GeV ); + double Energy = Momentum.mag() ; Momentum.unit() ; + //Set the gun to shoot particleGun->SetParticleDefinition(daughter) ; particleGun->SetParticleMomentumDirection(Momentum) ; @@ -639,6 +641,14 @@ void EventGeneratorTransfertToResonance::ResonanceDecay( G4double EnergyHeavy particleGun->SetParticlePosition(G4ThreeVector(x0, y0, z0)) ; // Shoot the Daugter particleGun->GeneratePrimaryVertex(anEvent) ; + + // get theta and phi in the world frame + G4double theta_world = Momentum.theta(); + G4double phi_world = Momentum.phi(); + if (phi_world < 1e-6) phi_world += 2*pi; + // write angles in ROOT file + m_InitConditions->SetICEmittedAngleThetaLabWorldFrame(theta_world / deg); + m_InitConditions->SetICEmittedAnglePhiWorldFrame(phi_world / deg); if (m_ShootDecayProduct) { @@ -653,7 +663,7 @@ void EventGeneratorTransfertToResonance::ResonanceDecay( G4double EnergyHeavy Energy = Momentum.mag() ; Momentum.unit() ; //Set the gun to shoot - particleGun->SetParticleDefinition(neutron) ; + particleGun->SetParticleDefinition(neutron) ; particleGun->SetParticleMomentumDirection(Momentum) ; particleGun->SetParticleEnergy(Energy) ; particleGun->SetParticlePosition(G4ThreeVector(x0, y0, z0)) ; @@ -681,33 +691,12 @@ void EventGeneratorTransfertToResonance::ResonanceDecay( G4double EnergyHeavy jj++; } - delete masses ; + } + delete masses ; } } -//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... -vector<double> EventGeneratorTransfertToResonance::PhaseSpaceUniformGenerator( int N , double R) - { - vector<double> V ; - V.reserve(N) ; - double Norme ; - double Buffer ; - - for(int i = 0 ; i< N ; i++) - { - V.push_back( Buffer = CLHEP::RandFlat::shoot() ); - Norme += Buffer*Buffer ; - } - - Norme = sqrt(Norme) ; - - for(int i = 0 ; i< N ; i++) - V[i] = V[i] / Norme ; - - return V; - } - //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... void EventGeneratorTransfertToResonance::SetEverything( string name1 , string name2 , diff --git a/NPSimulation/src/Plastic.cc b/NPSimulation/src/Plastic.cc index 1f95f4954..3c8142b5e 100644 --- a/NPSimulation/src/Plastic.cc +++ b/NPSimulation/src/Plastic.cc @@ -43,7 +43,7 @@ #include "Plastic.hh" #include "PlasticScorers.hh" #include "RootOutput.h" - +using namespace PLASTIC; // CLHEP header #include "CLHEP/Random/RandGauss.h" @@ -56,7 +56,7 @@ namespace PLASTIC { // Energy and time Resolution const G4double ResoTime = 4.2 ;// = 10ns of Resolution // Unit is MeV/2.35 - const G4double ResoEnergy = 0.42 ;// = 1MeV of Resolution // Unit is MeV/2.35 + const G4double ResoEnergy = 5.0 ;// Resolution in % } @@ -393,7 +393,7 @@ void Plastic::ReadSensitive(const G4Event* event) G4double E = *(Energy_itr->second) ; if (ETrackID == NTrackID) { - m_Event->SetEnergy(RandGauss::shoot(E, ResoEnergy)) ; + m_Event->SetEnergy(RandGauss::shoot(E, E*ResoEnergy/100./2.35)) ; } Energy_itr++; @@ -419,8 +419,10 @@ void Plastic::ReadSensitive(const G4Event* event) } // clear map for next event - DetectorNumberHitMap ->clear(); - EnergyHitMap ->clear() ; + TimeHitMap ->clear() ; + DetectorNumberHitMap ->clear() ; + EnergyHitMap ->clear() ; + } //////////////////////////////////////////////////////////////// @@ -449,10 +451,10 @@ void Plastic::InitializeMaterial() // Plastic - density = 1.243 * g / cm3; + density = 1.032 * g / cm3; m_MaterialPlastic = new G4Material("Plastic", density, ncomponents = 2); m_MaterialPlastic->AddElement(H , natoms = 10); - m_MaterialPlastic->AddElement(C , natoms = 14); + m_MaterialPlastic->AddElement(C , natoms = 9); } diff --git a/NPSimulation/src/PlasticScorer.cc b/NPSimulation/src/PlasticScorer.cc index 4cf0904c2..a6de658bd 100644 --- a/NPSimulation/src/PlasticScorer.cc +++ b/NPSimulation/src/PlasticScorer.cc @@ -12,7 +12,7 @@ * Last update : * *---------------------------------------------------------------------------* * Decription: * - * This class old some of the Platic Scorer, shared by different detector. * + * This class old some of the Platic Scorer. * *---------------------------------------------------------------------------* * Comment: * *****************************************************************************/ -- GitLab