From 49fd15989dac21d3c1fb702ad8c98b05aa9ad6c4 2017-06-05 21:04:32 From: Tom Bannink Date: 2017-06-05 21:04:32 Subject: [PATCH] Add ETMT plots --- diff --git a/plots/ETMTdistribution.pdf b/plots/ETMTdistribution.pdf new file mode 100644 index 0000000000000000000000000000000000000000..4765394217654d7169c8ce1f55b81c7e31319f49 GIT binary patch literal 34950 zcmV(#K;*wAP((&8F)lR4?5av(28Y+-a|L}g=dWMv9IJ_>Vma%Ev{3V58YeM^fjNs`?@enl?37L^nJb_Yod zff#7eLNm1+W+4zWFk}lPE&A`v&%-<-yQ`X8&=lNC^YDnzJ;GmRcI?>yx$e)>XT9gg z{n_q5M%#akdyQ2-{^=h-{^5V#AOHUUuJoB&{`m8iKmK^#AOCUv%fIeg_s9Rb%4cpL z|NSHHo7czsoV|Sf3DT^T#>dx3PI1kWu8jmKf0ojZcx|~&$xWAMd5Puo_SGe^lnike ziEUcEAlJCOn#4Bd7T296&MdJ*Y;m;8mpFCJO?H80R z?`A+Odf3?2DeoWL;@p2gOm~y0+NlhQt#`_uwb$n@IoE!vpd=N{CD}KWDsk#^;ae#C zXFH!?*_ADqlrI#t)K$<_0Uv_$Oq4st6!(vsC6{|ap1av!7YbTq!OFg2qQbr}qeCTs z-i>mqNpjKk+@|E#+Q+pE%JdC$znLvJ{@@-1jbiUra?$1K%T1Tp^yPi`?qA*`nC+?i z{sfutXN%>g*h-V!ba|91dFnbp{qA4W{pB+aKCgY|H6YhiCAUhs`|{A`mH50pl5n|x zrdcl)@c#N)dy?F0fy^PlMa;z1c-xlu;Hs2MYm)QY!+jTYwP)yl|2_7gNS%{x%X_F5 z${b5u#F$9~BU5rK>q0?UE*PpX8vC|%|NKn3NzSuDwlpB^LU#SCH!tU#cYS?qx$Pg? za^=w2@|ct-n_Kr7fXwxRBCGFDIXIzr!PpZHYX*lx%^Xu(9=S-)xV}folGrju4fQSW z+bYe-c;vo-g5J0-`yAA z%!s?^Z8`5ry=AS9TW(ksP%tGAT%YkFCzm$JTqCxebv>max$O-Fa?zXG1F1a)A9W8T zSI9|n*7aLIxZ(k~hQi_qyz=y#+V+^bD;rps3V0+TZ#I6uewx# zm30>+E?{kss>)$Bxg2g8+YAOB+r5X)UZ`^8*rQzMyf{kr1C1-$9>p4lTV9j$)SK*j zK>h4oAGcih4`a97@8@E{C@d(P7Hz)<$+ucCbVG>-?sr@u4lb~5THH&V+cgUP3TdcT$Pj#;ETi%oX;;O&}dKWBINL**PoH)U@+_^|>x$1gLmPcBU z%l@Xje|`Uw3K9qZmTO;Fl~!Mpc>~Wh6!9Xr$dUBX3uf>hnztA5zCT!Eg~Z9#wX)}1 zf8(x%GYRYN&nB4*@)md13vyF$_s7HGVDS}id4DX)S)JVdc}Yeex4dWaO}W=5xmIimQ!B}qv4nWV-x8+Zmb9Svh1!s@?o-B62w6|Z zfvWLoT5-pINvfc)yX9IXC(C>Qrg@XRPh^m%KHk`0ZK;4Sd&^yN=NsL!Z|MBEBS-zl zEpN-UOm2DldZVBKGPmn3?!z%fZW&u%rE1Hh%zYd?>grrqc0nyH&==lvR|V~HTgG&;k2nAT)58D>*R577vN|HSzX*SZ0iY$M;pmlh0pPI-EqC|jrU>! zv?Ua9Laf9s8Ry?oB^O<9T-~?aecjb@`<1ZGv0zWAyl}m8wcK(UlB?zJ>s%tXJbnG< z{<%urP`Dy(x$g}XGPnLMW3KbO&IZb%p;GHe)d|Q@6v}II)%#XS{PbrY8 zaBR7B%Uis*+W{y6oE#M1dfGB^A+d4hm`!7;k!*_LDs%C)!;K<0|H<>KpIVpONT zmu@M7C}4X~&_jdaCqoHFL&`pWr8~tvH$Nf@a70eo=cj$63KC8qmrF1LP%?%KCI)lT zQnm{UPh%;0FruKiM(!@yw_H@drXEsxKt(~10(>_VwCsZEf~0~Lx9ko4a!FBvTh>0M zF%EvPboWnbj3W=ru78Kz7F*RUz^I0ZcYrCU1?l*~0<0OJ7)reE0Pri>E( z$~w2yC@7h)N2+9k4n0^<)e&gPxw2)ot4flsAbCASLCGqpQKPvAXLO>btAebS@w>Ai zP0MOJrx;~5R9aoaa>M$CXJyEH!Ms!wHE1zT1ky>=N`ADaWr7AXt?fz z;ew&3cJ=}Vy<9Gdds;;vSuV0-fRfecPX&_PByL2>Jq5W zETCMu$3oWlFVZSmpzbwJV-)0dxg@G+{KJCA@u5Rsf<=t7&S2O!A!v1gIjrh+HSK&b za#g;a=GQfObYX#IylcoB`MY4bf?BkAeuT+%br<*vRTBZ0uQF8t8HO**TITC5tvrR2 zh!V`AbjLbhwOp|vxd2KZ3iRP&{&EeC`YIzQ3N*W50CNrP1iPT>M>N@B0^@M?&H#&- zhP(_ND9{bwDq!={0F^No1xwS$n+L33t}*n!c1d*=odIC^a+O|O5D(FGtBqJ-NC7XHXPqzo{b%@^sT+|KZ@t+L7zz*+1!np&- z8$pBVe6VTW37TdR3)qsT?VLt?%a~P=TNF_3Mggrp31%c$(RhI2$OSqW!g5sCyP!zh zzIg&iLA84Qlr-trx9x5ypo0A_pu+q-wMk=UftuD&6%FrQur$1*U{FDO`a4bq$rGMJn?h=W`f^`gw-t9qdQ^vT`m!)|d9Np^s1PEIu{Be3aBG&5m7w zrJV)p>O2KnO<^O_mIr$)3Ko^>bMhT;N7s_EI*oMd zTXc0zd6sx_upoW`9;p+I7p|mes5z?SRBIJt&P{N-Z8YkOZ z_6`2QMxFbSJLK!t<4z^|c|$`w^lgC}_xS=tM)4H>$SVywaQM|&V^aBkzT#M?Km%nL zz+Rp1IAYFQCMpil?vkjY7e8+*fhjzjX{J-4wPt3{5Y4vBB|W`Y%2}jtpIN;Nx*U#P z7huJ4!L(ec{=Ev~4lU19&a-a?{cf*2niTBYI((X)x;KAR%^K3DpLPyvj?YR0Q%{|* z1%1-FalN}grB9!=Wf#Ola@<2@`L;1X%_U`-OUvbw;woA?c0qF$eO=guRKp);ag{l1 z&4R1={UsNj`OiwTQq&(7u_r$Z&5|KPqvdi*RMGfZ^6Rr8xnSBW?gH7wp0#Mvoqv5V ziF;_9N{jxjZ?YeuGhAtVa2;{A|aFyT@GvD%bzvU)b zAcieNv3z9tn3r5?7YkgW>K8g+^s;sb8NTU99_RY9y!FHDRUilA^NW+};ai5L1)10f zm{`h@FBW7~aFqczazHD3h`UU$?&}7l>bVGO6bi4ZYcST4-_Po zQ$aV-3H%nvc>P?X>PGo;0Zh+NgXP$L%TUvv4i65fx4iYb>v@H+9=zUYevJpX#J)+Y zs2|GoV_x41prf9qZ{;6a(OZU6bb6ZQN3?vg0P4$g1N?~JL~YvsIeHP_qW17E2+h{S zt6yZY3>D+)r@V}hEHC33X!;cI_J-^0!xles=@+gQK*rJlCe7t37 zXsJ9HI;%Px+=N`GA`kKB)zQ3mngDS@9auFQhEu-dnXHau% zyS^3l^%^oozkTE?9qH8H`YF_cXZlt-ZO<=Odf&duzCkVCvg=zeLd~2$LtWLFI&#;b z(>tqj*FCb_b$;m=3$$VhELGuY`X*|9Q|lUz2QNu?ZWZH_4CC^$r3*4HJ#vrh$s#Ts zaZri4)Hu+m$WCvVwu57`emGSrFdpzE5{ zrne0JlCRw4yc6W`>*Jp< zIJR!+;jkQ~y8zIb>fYrlaEyNeP)yuqdyi5rI{hH(y@(DMSO>v{7Q>}YgioEV-YUCq~ ztshry{<^TW?l-v?nrWaTx{GyX8mEZUVLi=hJQ;J~zy^z>cU^%JZB%9e!{Y8<Nc$u|M#1y$fJYuHl%+S3w7|wahA8PE9r(R@pyc z*Cj5Z-n0?mG1JgS^ta4_8e)zSe@T_9%Dsnf-$p(H7X;I`ER|JmS?}qBHtr6WW?V~^ zxqi#iKZR|?NA~rCqxgk_vd z!qcGGKdypWY&lP!y|F>=kCBUAFBqHv z#$~Ub7`r!tX$PZd;F*vUjW%Q=)1(&%2JtqS3ufX$To{%6Cg&WmU;IKX$-yWUM-W9U z@6$WU7U{E|b}u*K0f42_{h1xa0@{s4B=iHWPIi!6Dpyz`6^o!e{@bV}`K16tyopQF z{KM|s1tgZe8mL6wAglP68&OE=IklV=$k%2U-+?eTc|i%Xzf?daym@)qG8O&i!P;zo ziP7H-*=Po_G#kiLlTW7!Hk+FL9&H3FNsebg7K>0ys;`qUU`b1$&LOuz43;*btE7XG z1pF;kyY<_WJF?ni+g2dLy;?P!2Z{j=TW0r)(Fa|Y89>MdvXi4E3E57W-Z5;p-5AM$ z5oKc@*}v@5A~27Hq-ZMjBov)a*SmXN7!T~_99msq7Av=cK<*boP*0OgUtqSQsq2S{ zxmgB5(wm$;7q{GWgB(2Fu#jx$O|JohDTz@7{h)e`9IWbLfN*I2CIlvMlf=1Y+r_%Q z7+I?x3$Qa0yTVt{)x%_ea36BqI;fzJ@P^BA);(#A#$bU#7lM-YX$Xq8hbUCHzzB>in?ingY%O*Y*DgCBP|zO?xc`5_@wZ5i?|r|_tc zWS_GRyCZ4EO98tdI_{-#QjcUSn!w`YSob=ulp}}rls?xB1vdOPT_ytUHeO}ugqUug z5CpK4>qA3<2+OHcth!giH>RCwRH8>ASIhQ`3K3cD4wk=eG4lAj>V6em2$JGWSeAv&kWZV8_mY(vByyd;3UdN6-^YM0aV9#(Mb6U&; z0dQ*yk3(;{)p9+S^;u5y^|;Cd&}nTww(DvDW?1Sox4>CX^Z7ZZaHYpoJzMy%hNM1h zWL}c(%Dfv-*t19PWKrI-J=HQo&mK?!KVZ;@a-PrQ=F3t(9cy!%H^Qi^EgmNS;r=Db z;b-VPR88Tqk9R|O&!|8d2gK*|Ea7v{3$r4b!sEMc0mQ@gxCEb1;b>e6dr0BOSwdRb zP=B$Y3FMgXF=8Go<{LK_~8RZb762CjPi-d*B&`-QL7 zqQQD2s~Iq(>P_&ecps~WtzI6L)uVa#%>vx%<|Lyz34gVQy+*WFYQUP@0grcJA(+9w z5?)_no}I>aw(@B6TeSV$fJ%8x!deBHK_6iQL_w38)A&lQ*~{Q<;ym+XOd!lt^&wv@ zkiklLkUNJ#u2@8Go}VkWb0|~}vfCRN(JDk^sC5Kyw)HM(g8v?8Be0u999URcM&uF5 zOgBib(z_E669zuuGO-Rm9URfZJFA0}>ulIw2+CZ?s049Z)*D|ow*(q>@dS6h3v5;s zG6jb^+Dso=R$JFk!Yds!ln{&QX%LRc`51!;PHc}Ek5I%}o-d5R1YIvPkcl$W`HkRs zFN=WVt~08xLIKn2)6q z7{ihtYaV2%wvXK0et+Z2*wyP*ka!+HH`HW*e8^d*Uw?8~t!YPB123Vfo(=#*MeCwx zWv$_JGOnvqUpSw<3Q2pM9jTsBWrLQ4m_B1*%SFq?3aHc+>c<}$$SChyr=YL)`s!5y zEYBlG_)Gi~<9K(*K~nX{@dO~%o!_09Jn%f7oDn0f>hxn@&%~&L#fNzsa@Wk&SFgHq zUOKJO=bK-UyRPdSwPb(IBiFaZ_gBovt?$D)qTly}03L)lV;`*Bub7WZO#1G?bL{`e zc3s|0!|+WpIxorndklMeMmM}@^cvVoU+m^`2z!bB<<7Q5E8^WT?8WcEE(cjMjyL11 z?)Dc^6{2$CXBl-f+z@-0mZRn4B$Pt#5!5ajbLu)rClnaa#G?}IEK*veBjU8^5622JNL#IPO9SRNvx(9(`&GK1> zw!8;?xx>MNvlZPK5}@FgYc7yhp@Ro?O0RnQQ(* z95=2WVc1zc>2{UjFxR_H*ICm9=uUl|#nl_ZzgpYi8{UguUpBu1KCwsl2r~-9912ZI zw{O6)9MQEUe5I22OP;qP0O7;f-x{VA`Pn!IX7W@hQe|g^ zSCv~J;u-O1MYgD!!x<(J5=<+;j5)W(cH0jLKjF;Uh%Wk6D(mY@+Q$OpOSOzyf#}vL z5h^WgG0T^sQcR*$FBRa)vv_MPGW5-2Yp)qzrfjy7JL_w1!RLeRiUsa5MlobNO7WE; zmT;u*Sv_Me6or^h8U92O1D$3TyLJ%KQh0;DcTqeOI!vQQQPkOtsMsR4_rfca`!CUb zh*vtpS9|k@)az_3#k~v!XzavCKh$c$~rkDTpTIrQmCOhOzNqupNqoMtb}h-Wv*IRmZT+u_*L) zbgyMa7TOV>v8yQBb{yAm=ITAFRTRjbatRq0T};;_@2NoGGj|8L3HOjDRN57=4-iUi z8J`N_`p4J3&ygbnQ3}U07Z|eyDb}19F}-)`)Zk)7PQnbrO}(!7Q9` zKU_A23ow?v;Ttg8d$^2>{GYuAl}(Y41&)aDzQd5a@^s9llQ zY7tKDJ>rb9mtql4F;BI(_S=E5ML6ajhH=PPgxj*uqKH_c9e@H)XL}1d{{DR*0xNFc z4xB_ARD>(oAe_U?54EX(Ew1`F9tg-k8nl5)FP4?zVCgOXj+Fa z6%20K?o#7A+#KiyaZt32U#NQq{P_78dsm$?%NkWiIF4JOBK8{u3ZWQ5fuRCN^>Tx; zn^m(bT?|~CbaW40hH)6mc{~p4BOK>}vA~zA%?Cord+b*ZNHax2^?Oh-&jNXW zJRmGsUF$&MF;gPIJHyB*;7=+qdvXt|op?yx&M-nWni>*mws^oT0bs$X<2Uqm_--*I zXyz^>@EK;AK|vQ_qqs%rNI6TePRayuNQCPflQf46_yqABZUFR&65JJ#_lqq-yHgZ` z<%mqlQ+N*{m+_`*5uKVLGZ@At8%-pN7VHcfJR*j{*3yl*!7wP?QVq!=z$C>hoUF;6 zbxOwCkFC1+fmqs0oTvfe$vpgc_D&w|N(8au%V(-H)sNecEnyufQzS)Uu)2z7F;0fn z&6&I|+6st^m%YPI>TP+(4o^*MC|tZKj62UbK;d{CDhNfaU^#T1vJM%SHa{GWslv$R z0U1b?9`knbBSYg{7DW{1ku6M-b#Qb|oF7o~-By}?ax+$_sL#2|>|f-C$Y*0)TtH!8 zkpZ>H2nUPJ?Fumcz(6GAjjV*ZcxG2Qd!YI<3Yyp`@v7;-TjjkmVw^H5nX(&P<@^Nd zQdlJUhJoZuoGeu_uQ_@5SNJSsIeD?jR4VVXPC;Ho+UOM#jWwz5OL?wxX2Z7a>tq9q zn3?Kej5x1BZj{XyW;hqj>Cpfhf0hj%RWA0W6L~eM!?LNUCRu;z>@&J6(U_sR%FthO zF5Tw6GpnA?&ozo2TG=e%$ z{JWVtltv|xQKwjh3!z}BMl>RSndDU2{;u}nZrFbsb6*{j=->nWqo3`#irk{LT zA*dn+OVda2DpdrsMa-QvV|m1(ut8~TZ3l`ECjE6X}@q3B_t(vkYxb!YU zR5?Pl`M?8G$gy?1J0yx623_C;Q&`A1?tm8uW3wFS|Z zj1*Obfybl^AJQJl7vhzgZ(vMmwV9%7yGfxVf}pbSdWj#@@2f+=AwxM%aax5ejGb30 z(K^%sur`_BYkU!hVAhB!1b~zlY3UF|-DrKLQ?%OzJ)6(MEH<~xB5UK?M0dEXEGKZ@{*j56}GoV%gGGZ9`F8> zsG~!`1i2^CC;8`^xskFT-$b5d@VE>;C+?#9q=|y!3iQ&aV0)bcj3%GeAx|*?C1aTl zbixL%4m$F~JGk)g>5eUX32_>jz`6o9~J!5`@sYkTa?KVw!7Y?u85i#F=cYVE%b{j?~NM*Z&7D;cfMZonQ z@P1g=_te{Ncev}#sVcsp6f8mQ4j2^2Un32YzFu?vB^a4WXr*LPX5Z4`2O1KS0G&g+ zBB1Fdp90ErJfowDoa}YK`EkzU4U@C;rKN~-m#L?Tto5}2(51M0DNuxd9ZAGbDRe{7 zsdzLdL+O`@U;Msv&@{;#1bG~uB6@L*Vluuwi`Dad`@^dvdvD8uTQrF;DI+vGZ~vgF zS!4{`vJkVrXq2N!i&}Qjmp0y(#nO~+Q!t{y6$1D_&0plp#t@3Y9WZ|djq&is-sxp) zp>8Dzx1?UQ5)jyG(GW!KWRMkfW}@m=f>24Gsirve6pZqa^)VIO1D}`%pA3(!9R3_x z0$K^`sp|V3kuR`0EQU#d;^YPzLCD8lt(>U0+x;`Kmrk4Qum*?$zm@#=kZ21)^rEmbqzR?wyzh> zPtYCdT9}LH-n&9Y~AQ z;jwZYPyY_8MzD*KF*69cjf|1H1e^UEGfNx3LDp4rpp=8*yRG_SJ!(C$GQyjH@!Oc( z>=qg-V@mOUV9u$yMaIgch{p-$A>Z%iSeExTjM3K%d(vficNmAAhNqo2O#GzS3G7o2U+7zx>XY0#O7+ZsnQ=bUNP2G#^^PSHKl%%5I7PJTGf; z9-LRX46uD-r^@|MgtXPVg~Ss7(9A{5mkAnE#E2P| z>xid24&){d?f0P}E55<84bmz49gLWoXEwVa8CFMSc;6y+v$8 z^YaX5ul~JWU!R;|bL$xrJvmf|gq)NkG~)K}In)S2BjrswC7+HNUmJ#v%EO6&l9(2S|wtV_rt-knjbap zalojK0t{?>&pegCL?h{WXo5J(gVJ`$V4z}q@iKSBRLVinI-Kw&dR^Otf75x?9Wig&+&z+E|So@Y%S!p=`k8Y$jlYfm#}qJNA1)y073`wdJ{{m6od2*B6N zQryo=T!H%KlY;k1YcWKoz4B*LkWhTX6(JS+z{_ZPuf=@@OH^cNDDA-ncf}W-k+n`Z zMw4i=p;r8?W-_7(R7?+Q!lL2C z3%5c&&+Hv{oXR!1%^G;t<8eozkyeLs$Dtep@nS)F(V~qz4&_>Gj9od4(YWIbuGI_j zs0yXauHaL+mcMZyhRfox%?Nd>et=KmNt^=X_A%V49@6{oICA^m6i4nO+jwE`>^P*s zZd{5iW7+k6=+@FZk`LMLyD~9d=iQ{k<~AKM>O;1&P^36h`%|PiQ~Og9!g0MmW(xa+ zJ~m_%yD5nB@E2>aK<W*j4FY2g~P2b`{J`N|5!dEj7 z(_U1uCLC*71SuwFs3IR`bdIO?sM-;Ut$Of3=bahuh!MwEu_jzVAN zV_uQ5%k)Ai2-6V~5`^)h`I7DJAd|(zcBuQaTz8Q^0HV?f7|LI?N+nVMs6e6VCn9CY z4;gZ@gzI+Kk6$dL}VCoU`fh`6F9OYJQ-OmHgI%Annt=li_A-xLQUu@>5;gR z&d?IE1y$C{^5l8s!e0YeU)geawgfcG4mqvS;%ah-7LhC580zzU>X+4H@p63G$#KKX zSWia2e0+(9rdLHow$BiGPG_X515;=#>=7AYUt^EZ1#_lHDmX^I&^v1Ebv0MPut6{j zxS04)E=9@++z>C62Dar4U5uQ;T+hK-HNp>qNP$x{iyA`%EAJr|n?A;UTCfuRB<=&ClwXVG@#X21+1aZ0dXNp8g3 z!18mcU`&NLI6Q__G4nI-$Bvz6#5&M^m6xDojNKvfs41vC%V$eavSVkUq?gDy!_;hy ztY{%F0>7O8la5#Am5kt?MxdB^KjnuNzKtS*XYu!LA<2#0G*SrGz_kz}GzloVU4b-7 zS$4y0Kx2s*19-2W7z#x3z~?;iKwA_$T_ispKQ0c<%j_j zq^T8KikawP(HwcQu%<49nLos=^^&c|H(=XpGiV|k=SUCJZWae06tn0sVQe{ZiiWkYHTb!p`E2-o~4<| zO$hXNP}J7sNO~YUE(-Pxp6xTz%%WRhA-FBWB0gutyeNy<@~Me2O0eNOr4uVOMCi+D z;uIOHMnIbYoU28YL{}=_F3Qxeh=E$eUMR@WLxP~lAyEu-ck|`7ip+Y5Z4#;fiEMK& zhw_qXw5axc0wl+SZS&C&1w6U7L&q1!$8KjE3h130c46c&N6zQX$zM8yUa?CNM?83o zZ$tCVv7+W+jJcg!D-S1kiDWg|E_Nm^Dh{IJu{`AnVzk#iWqX>k2-dL~cE@4S4WG{H zH?X=~ERpjKB>+FA^NvX9l#~a}@=}4^S3%vE>Iob*eRT>rYBq!T3E?w(kvV&+YtXuq zf>VtlVa$KiV&@$)Jf{=WHtYxwny75TjMNPRP}~^g>M&IEXcNmLSxh=cRc8gvLq4yg ztG?M92UeYnC8iD4fMs5w?5+@7fr6+2%`1)A(TCF?5f84YGANRBe2tL`A9vJh@85dU{Jv|+^Ns2r@X-Urs9vj-{| zeuLSeWc%o-Hzs%sg5R<~ON6%(Pi8t|6%e7&tpeuhG|ug$-=4zE9M-NZL@)|`7BwYg zY7w|45>pOTw~P$58p{3|;tW!?vn`W~D^+XxSplKxIUP=CdSaFaw@Z1v%F}fEW0^G1 zy`?&Pa?+syXLMjsZa1_qj{Owzoywre)xft<2Tr%){*gNrcHpV}3^O*kX)8%UH#31B7fXBw2{6LNo}e9BOM1IY{-L2`kxbClu;4m|P24ut|L zrVvwPBzygf;9(Iwi1VxP5!DG$J%(&%homiAB6EJrv)fbcx{Z^Zl5RUHLlIK}6V{0VOtr9{D37 zZ@sxk8onYFPZv!bq)IPatL`0Md-JEkD{d$JQF!t;A{bDtek~4sMfh@f=98lnF}K@^ z?GxS{FpJ`7L^z^a@zJJ;;jTrzg2F8?Vb&aA2@4`&d8YEsWT$m(oNOW0XM>`Mwc5SV zUq!TsM!>2&N2i3Zm_#%fVs5vW)!%8cM3LX^sh9Tbyfalq-tcu?rTJ@|OpR7S)l_h>LWD*1g4 zslZ-~mHWHYGnaEvD?6e_@OcM?Qic)HDe5Gb`cd-Je{$!P$W&ScuP>3Blez3V*MNTF zQAN@RB>{$8!`PhelL9G<$7MEh z7c?s6YogH)62J|VVCa_W5j6fO3UuEfFX>YT^_d-$eDZu(IbKAEF+DV?PQ8(A;Cku8LUBJ;sIR zQ|?w4ic*h5PJTW#*W4n&Q)Z7W_3cu@Dl)q`M0XG8d`PqvlaPbKy@*Xahla?A-b}Mm zjkiZilh7$4&nf$wBVum`ruhifILlzyWBU<6%V`#tS!cw^)wpTMa|N-{m!d5bE#?uB z%Xv^DOr{7+mIGAJT$2h^OCe*Zn^=c3hI)>5k1IGSo@E0d|4Aq1b25oK!wAE;U(9Q> zJoLJ#Qx$6Ybp>zw$!}kjk1aT`Au0^Vk$&A_I8v7&j<5_nhL*f4V(FdTIG@=GhR#O=_{K zDwhpnT&81Q=?P#Q_~im6e0SibC~(qYbszEs)Hq>{@C1097zSnrOL`eW>zo}GY|FnV z?+D9(i~DK#SExL?C#+cwT8AM^e|3nea(bj1$Ce?(evN={2SiYI|C)%TPn0qZYOe?M z^&*~jknrpG5Nn)j@NNv5%IfC{5FRgy7O>C(ixkLIW28kg-0DM&j{uiA|3I!%WvCNmin8J2`+jl{hobh`%|%KNap0?}uQgtP>*pr|1_C1tB<6wo?G| zdg!eYLF^p$#UZ>Lsag4LR51FJ79KIhH-FEu{9=he+1qiFqgPK(!m>n9bP}_dAejD) zkYW0q(r08l;7Cld?yShMExY*T8014(hLGb_B%cuE_DvqSe5>#D+I#1A|E7LH5rswg zm6Re}=?;EwbqMl!C8wv0gO(w?N$MT(Os^PZip3GW4sW_6LYTwPu82!9$)a5IF{C>9>S?E7l{A-giO=Im<{0;G$(nmMQpMF$l2>o6>Iq+Z?AG4AF}qEf(>3%wJJS32dl$-i=OZB zr6jvvwmmHZlygMqCfn>HG%@j*e5GL9smbQGh!0J|Xs6E;ew3-Of_uh#g$axNic*Z{WMF~)&dju_|2L0NpiGa26z zz`MVr0uxe1{0cAxN1qOM|X%Upd<3B@?s}Rp#dWTXz(B zQ#|J_ift7*Dy&4vZ>is9N3m(0=n<#m2X!-|~qyFOtQ#WB6}jX2bzbUVH^DK}_@%0>cVO~K*R zY#ulxFrXog2eZaJIuAa$j)1=jvj(~SZ2(xy6YDwf<#$g&fQ>tqU(d2;re-}Wd9a%p zHyKT8BHZ+O7MY8pAYQiqC+Ao;xnXi9LOUK$Zbf->Q>3cA<-JG|%OJQ%3L@O>+A}D) zEh-Uab+b&1YyV_0toOJo5Zm;Vg>YL&xP9l$HSGo0dgIdW0H7Pk7I}2LAK(yGbLi?c zcp#{ddwsG0XY*tnQr&!Dpu-Q)gjQfzK4FZ`!jf~qK#PHI#er_a^&MhR14!p6#~Sq= za!{SlG{PQlUv{367Pz_y2~L6P^Li1jhmnO*mCOBgZt)SQKqG2RnVO~E-abU;90*$z zljLkiZ?MQ}7i65ZMIJ*GM@8HpqpaqGJK)wsFaDnQRR%FQg$wJI>7+ zVuq5-M^uN93GSG3NUr=EnY-iCzw@W2L>v&^wzZOH44ECd-4p zWXi!BP&O!=mvB)AQw~u#nX^|XmVAIHqEw<=+x4eX4&uxpa&Q{D&y*3Jg>w%|mZ96c zslrm0pnwyYQed(G_hIoV$E#X-8YqW3RgUvAPg+n+h3doXWm@l=VoFvX?rI4lS+RVZ z{X{ZodV4sZWlv?&L1O`-c7a|XuzzCci~Ox@W{Z?-S(Jt~-0-p{TmXa}Crjt$D)54CCBXM7Tz7tjP7P;)Q@Ej=ck|D1x3TD*cF&|<0!=2$d zzXdl?7B-dzD}b7`99Pjcs)V-V>Y()O%s{H*0#~>V5+rbngrx=vCW4RJ!NoXoyOkj7 z7p{$MKeY5&z4MCGh2|&6!2yX5RXHe8WkW0vIX>ZfEytj8QnD8QPUgG>Vyopq=XDXB zwfta5F^I@?VX^7!&U1b0sZ1i zIkHMKT~`Uzk>H>gOAf|Vme>`jY5_v?fkh9o88Acj!GltES~xn&aZ8W2Kn$YD4`4S=}K<&RDr_@eKrjU=_t$AubaY zG$I6NU?xRlkpkUs7OguBR473H*EOG<2DClA;IDX{ z&F)NQSl@flh)5=UPeT`NJd+tIvCp>X~qu1;DJc@0F`1a_Z+c1Vm<3A7S`PrrR?zQziaou#4&F1huev-S=1|A5A)6vQ7ag9B2UL2Z=@4){~2&oVtlWNx3Rh@2yJuH zTjb)#v9{eYK!m8R=O{1E3mhU_T<=ARwhzXZ+FlsPWMZrRc|Xk3(w;;frZl%gr;tq? zeFI~Q_=I|JU9ybXWAo)q9Pmfg!0;kI>2|ic&*;Ni2F&GxD!*qrF z*Wvdn!>_N1v4|iCBWUrNOMei;%%Rf`;ugb*wDQF=QtDJIr^N(QCG08Xc8VRtALL+y zE+2qd_!>)4N*1{!Ban;;TIS%%fkI=hD}wtXO$W>&K*dFS@tsSI-l9}kSBa#^#Yyag z2C*muN(qd2ks+=G4KsG-Op0$VmzO}pT~Oc$G#uDs%CLu&a3K&wT;oki-b3>$u|t+S z>>+o2GE!Bq2c=A^>J?cXxU$-w^2M&f+@Xc$i!B9KEwtYv29i7U?I~p8yZS*G;~=XrLd!jJmq`SktOn)0P*4xCcoM%K zjVQ+L!a9*7R-=eGTJg4Dh5+Xp;3FY`bEtinVZbz!1?)FX8Qsl&1ymkQvnGMy4nc!M zaCdhJ9^5UsdvFMD!QBZO+}+*X-QC?`Uy{%N@7}$8_Ut+P2J4!dmYSYPJ^fVIP)?&V zbM}4RjseFVDjh}oE+j?5ub z!HKCjEwuAWDJUhR)I7PojO4UD{5HaKck7K6d! z=ux(^4dePD=3enbqMcWR2PAy)4>0@?5Fb~Uill53d+goN)J3pNC<-eH#X#R0uWEho zsPdK@c;rr0tV4c>hLpMpyhF@brkbtu3zL8rWTMSHQ*}U&_w#BpMh^1Ym#(?^BsmZW z^kmg;#``nwM5j!?nQ9JQ9%d_AFEa|Mc*e4j!xxr+aJ8^cQSrM2Ef*4Rd~L#srUm4 zCbZ7b;pj2B;@W0zIxCjZ@uRq1r;$+}elnM2pjeZfX{S^7ZVaT188fN~LAa!K?Jk5q zKiKMTHothmc`iFq9A%_Qr5+ADBtF&z+;En-Sg{1UG}vLF!xw{sb`POAc@r=dJR#wcb>Ffnipz=^6AQ8Iq@aG7XZD+D4aF(+ z00T{5*g%nI8YJH~h!umAF55cq{S=sx1-Bsgbj_D)VJk9A4o~S}ec5E>WS~)XhEg!~! zvM*!?;%^@NAvA#!KJZv{D%fm-sHFC;6<_q=?++GbRdOJh;vTit`})K5RN zv_*Z3$rEYI(Vzt7=cIolM!e!b35!L3+;F$1#<+ZkRQO|f)uigXZbG_pY;r;JC8Xyx zi>@`~X?-3zc3*>6tFlN~7k;Y67^N$=l?MVdcdJS=udVRzFxYEf9Ups>_w7jkij;krbAM+`g# zxAs=7B$+vQ<_8`7t;h8D_)WD4rYJ&>Z3)w}QaQMTr&S9cVOM73f_yXJ1Eg4&S6?wg zT)9`fR0+ZTuDw-ty+tGF`}Ep^)(Z{jpc5}ubv^nbkjVp<8df51&Nldw`;1Ye4Oa874r%6ib{nC~n zMXp-npMEZu`!KpX_h6geVJ~6VU2r{@DqvvVmV-HKBOrkNIis=(3)T}B-ezh)*U2C8 zI~j5dkAK+wXYYu6m84Q8vRNHd4()rO25?92yZd(o=k#AR{5u@k16@E+2IUp5D9aEI zI`Y+$61Zmg%1T@fOqD)wGQ=fx@4rPO{d7SIvXnK(5UF?(bCph!UFdQ^71sn;{{Ra@#r+*(FpXy{q=;$O<6)`^l?Z0`mUCm?86zt z^M!)`aOePP+JV#s(hT9jK;3f14VOmEtvW6}M7Mo#$#$Of%gqYU^8pQ7WjyEwX;3H- z{kN$DZXNu?+&1B}fwW?e+nThKw3D#Qvw;?kmvjG1IZ2ONk9()kCz_W{o|pRp+>)2; zfs5rFGy5^jU}@91Z%Won>$cd%H;ZXC*yjachZZBql+DE*l`W}gIOiwMghYY)E+QHp z8!?->E+6s&DE1k9&kOphI4;|mC|gV_k0)V^b?y{0!vbJsVFyYub#G@kFk?Mdg5SKk zr39h&TXCh^0@oURcfI%XJX81SnWU&&^^w(j7s=vqEv9kyG3B6TtI(2MVNL5+y0#Pg z&Mk+=SeOfU92CEY6<4aO9!>G$8w%AeX%Tm(I?)$)K(jsK$+8+cY@ve(i(19TTUFhS z#UwhBea^`!ax}KRd3-g4a(wF#^}8EsP!l=La+v|UmyGS_VNEwl&%tij^TE$rO_t3o zBWX7G*Y!6%cO@%z^XuuO6Z4qp0Tc7&!EJK7v(R=Xj1l%Jm{h$DQX$2?u@f6O5!5`> z8^UIj&lf!$>K+9hKrx3@cZ}2L1!d_knkDc?cN*S6NxAKlcWdftAr!ZJ#ap;X$AYar zy=sP!XV2Fw`zKfoYA4e+ z?mss6SWsnHa{I@h8SDJ^o{Vi!Z2JmA^CUFw6B#>lrJzSQ8H$`vo#6t^d%rz5t(1Z@ z#}~@j`E;Tui0!h&QkWksdyKvCJf5Vv-@!*`?SwpD+SXm*J~K4DJZY?n-F4jW?M#;K zoUS;EKA3SFaPer9?KgnJd{Ji&&{4|K*cnmk9&;eN6H)!Pt&*};@w9K1GXzFU`VJeP z#)voJbZE~aJ1r{lKJEBH>SasGx!R#>Wml!xykreS7hmUXvip#e0HGHX!xAXA1k;Lw zO!EFfhTGa+UwPBo&804+$>YuOM<@32`xN%(f$8Tpc|W)G;W5?`f`kyCKq~5+$~EVhyoI#C(V-7+SQw@EHaX* zl?D$L&{Dms9O`JW;>JeJ=@fT{przjcZAk1pZS$U9l-!QeV$lhu_i>PYocV#bcnS$D?Cr#nbqyEdb~lK;0JaH8Fuxj ze@cHJt{k9M8yyQ#BZ9oO2;D z#u<{D&oWLV=75DAt*nfWmHlCdriw7|pnvea+LyrB>m4*{FCpFtuoum{Q20)?AHByQ zFre@XoIgm{5ovaA-=d`5$e*frW{z(b&q0=4Tr96VZQR_saXUy)7cw!JUKgr7a4?v1 zQy~D23LvnC=M3eUKhB=c5(7=B067m-`d(xgg`Cj5MbHP595xSO*j}8oy^#jG76#%= zpmP~BC-i@MeAm1tbjXnI?rgqP)*h?|q>eOwRnh}>fVE*jsia*e@|1k%fMNh~qU{73 zdtoq=q6lquAxKG~j6<7yT`b%@bPE*Q(p5g`*mK$KvBGcy{Z$?ta}G8dgF>6%w2C$> zC9(q3JisCfnRj)>CVl4I5R&MCv0bGfJHrh>)Qr%r6%m03)gCPbmF8W@KDF29=q~a2 ztD!U1##p;Sh*K-kb-?a5;72f|6-eGy#`{&4@5oIPMB>S2<(o{W1A7Mbl07ZGQLLj7 z(ZR?85s|vi5JK$w1v7DRAcF19Qki&3gHFH&p@vhy5pPw}`OCk1zl)eMnlL&un`#6~ z^mCdlUeDtHOm>(gke^PzdZogMY3}~=4EnTddGvvVg!@tG+V%OkO`MYkmPaupr!A|O zr4jT986V}k$}~%F3?(SpGCOijI)Xu4q?{mu&~P)a(SS$bj`N4tn?wAM zkQr4TE$R(RCdhf+anloNc=Hci*EIp|T8WAbd25@YzC{=QgY0*k;K9cnr7nGuOAT>U zb$LOlM4WYNUHjEmo0(Ruu6iScG*Lb6@M#><`+bGD2E{rn#SmX-GAtKQLl!;4X+PS2 zVl8C*;U;hsE6&lbY@Cd~qForQWu}o_y6oJ0A9GO(=hKTk9X9Dm&sBm5)5WP%&ExMu{+dZ593AlL>PCZ-kNgcFN9Z!NY{y3?Py^GAhqvo-K_$9xWPszESb_3CPGRpcm$U$ljkkEX8VHSl2m_W$mO6u!*7q?!SkHBiZ~g1CHo@JIaAt^u-Doa;!t?-_nMGV{gG*AO zl{Fjg0;NA0&u@yBz8UdYb&6DLPX2iPft?j&mM^o~p6lG($XfOBbIy&YL*muHPpZ zF){j14<^>pl`TCyi92AM))rZ8#h%vWYKDZnU@c6OQ^av5+_E;A?p4MVD2hHT$B(Hh zJsIKNP~3d+^SH;P`a-}h^c~~|6dwi)C1of>y)+)W0 zTaDX{7!yLRT?u;5k|BSLb(+?BU*hdr$t~g^c<#u6hVI=VG#%;ljSl3ca=Hb)>m@Eq z-X53$Yr31KGU$)Z>x+&Dm(N1l+HJl?nET;;e4?=1#JdXYbz&;B6k^rTsH{h(heJWB zzChAScsXdmaBs029R#8sXY!nUa+r2u= zpLs~bm{S*Shyo@LyJ&YI$|KH`G}5jDu{7Tl#C_d5mU*=1D*FUe^D-_u;*-|db@s`@ zJE#qAaJt=k@n%9FoC=r=u-~+<6zdwe1!TeDRGK@(-Tf-*W@s;NWn|j=@z%OpI<(zw z-0xP>7se>Yg_$Tyw#bY^k~+qsOez|)mfbwXBDP8@u|#cEwMcD)J^P{6VvOMHf{Pe; z-$1z8-d3!nuH-w7IDRR`jj#A`&68K95XTn-iM@KG%}Z&TMCIuykPw6-g4F1*UH_bXP?z*2sG z(?KHI2&Dykqn1QyE7-9aDJo~&_Z{PC$)UMs|KfD6XWFA-XX4n4?k2Q(V=u*l#yvF+ zfrr9Ses~p^J3$fo`R-}ESU)S762C!f>1sY(x*)X#hU>X^xoP*y<@MW34f=a(NY`+X z+$}Fxjy-c+9o)5-`CE-pggAlNcG+^Dw2~yl+K$@RN>*bK#;A}-%q6G}@(s+1oGb@Q z5#&PS-A={0O+AUoa*Tn<$H9u@N|AN zV%yqGjyVi^J9VCnj07Jg#-VX(OY*NXaW_w|=_zg;mcTo>uh1GE<(ziG^T22|TJ9(M z-A^4@CY4#N22dd@9?nr&n?oh zC!;qRh&S}6NL$JN>_)}@Z0b%syN1oix5@}<7l)42ve_)kI-(&ldg|Ca;0Bx1C`8~S z=sIHEjc9!uNlnK+W?>Jhno*W?nd7MI3h%vG;-T5eduoV&q)lzJeKjJ7wS!$bR@U*I zl^wg$JM+ox+72m1>(dJ-=GWXw7b`GUP{)hlCs40R7-%nj#8nOKC!{5IL z>os&E<<9I?IZF}VQzhQM`#h~C2l4hJUU-=o_PuAmlUZ2-he_m;!I{eWU=JJDTwg4-PuO9P<~$AS;1PFD z#W3^YTqIoWIYYu6*6vYy`^v^XEDbkk5KDuTy#*;Y3Kq~9f&+#dOw?6!H{!i0$(^*-R( zGB>c6>rwv}cg?choCbUQ{JHk!N`$3+NWoHdloiA#>6| z(Aw~TL29xd2`3m|dDz{UEY&bdTyT<*X6MACmMfVGe1t1)ls_#h`Y2}0$&Hl#$W`wB zlx@?l>W*r7C-vE(sz%z^uOjb$jLFzNBrDk;orM}Yd(CWbA^2*GPtrR`J9kx^U-w^g zy|-MRHN9!N>^P3=DXcGv+TFISKQ&;13>tc`xl-3;FHx5AeEZpPE&8;prO{gJ0OQjP z*|G*LwKJUbsW98KJ@VR_h`lRT>ImiIz;~Z8hQvhH(DJ-684&JJ_M!TmS`UA^tVDUa zV)Qxnax9Sf&D&WskqX!XqeBtepe4-PD=P$z93 znEe|1AClm$C|-2AUJg=2;C+kMR`_vG@K>T(a6(iE%6wskSNDuIEgG4u5W>KXd?fIy zHmS?MElP)tOB1z36T(JoUgii99KwNd^g7CVABpQ8LC$JAWAUS?n`1o>w`%Nb*BDoeW(us%F z3f9B>l4@fuA;MCHW8>0A9Q8|@IEK?feoZ>HDjunG;1m!mPzT9}LzN_>@%6O_pc zdb-g6EdFhwM>vtlJpTJsL3Pf+GHs+qm;qY)N@S5u4#9JiW1?%-kOMkcF>Oje|D@jD zUf_7(+&h&+aXW233iUweyNYrffbM${k$YC20FyuqyfXoJ z5-CN!*hHW{hy{DNE+0T;yC!`K)+pk-slcv83svVuNg3)26gY}T{Jm;KXQZTW#1128 zYf;8sSr4N4a>@$#g831Hn?x#7vOB!5lI8egn61O=e1B)N+jlsFj27Is8V_M$JszfU z%q&z}$d$|vzR6|KaH29$T%=lbjW`470P6a~~aSw*6cIY>qFfkj2^CYaC%4|Y7+e^Bvl*kZ^( zqKOdL&2EXbHL3`#mrr^ab!GPnZL_0Dk=hxfx^$nHk{;7yqcU6_Ob!`X^|Mu2w6ox@ z30UYSTjW^GT|`?HI!4NYDB0n1#XQ*od9-Ry)N17!_UP~+lt5fYZ64ocvSl=7#2*Tq zf1YofFPpcXCqFVcaz6sE#}wQ_oCT-RYS`|>umDx+KJs=k=HX6g7GDwKR1&~u2_lf) zqb?7n{%}<3Zxv^NfGsH8;j|pq%q^Cdj<$zq@ot%!xRr9A<(^rR;DxiC;2Fa1%-_o5 zWIgF*9KDa;SxVPO5x5ee=aC_<1H86ojqsRp>|DCNr!kc%g8Ci{uyS<@#iRTKh~|d0 zv2&<{qmNq^Ed52Y2GN=`ku&U}T`*Goxlib8$3+P?R4vUKj#0XYH??}t!K#jI&&XGo z@3<GS}MSEH}Ck{>= zi$G531p)xku--ovOUL@BI>RTx04(f7tR~E+ggh%}1qlk7|XMS|lzKo$u0d`M| zrss!@nu$P==AKI=81J@rKBMgyDCBh3a=LTNfMd#ourwY|28@M8jQdd%NE1x~*yepNN+r6!9A24DHaXr#UN)$6g*z0S?H{6SA2 zjNbWXe^2^kP5eR($047hE8(!!h5s{fch|OK)L44hgfSEFhzxLoezOi!ywIy|dQ5XDG#t`!xE9T_v5)h%0Jth^*zCEaRIzskOU)*3}Wn@rPE zUonT@PwGNap_N6D4qP4uJs&R9$d6pr9^IHY3DzuJm7XJs#ZNF^6`6MWZr>zVo>blZ z!}|`TR%`K?tS&=rrQlnegVCU>nd%&?A^p$&7IR*| zq+cC9gIg<;6kp(`P5IK?e5=*Wu4QI``_wNSMy@RVZI6QB0Fs2HF(@0yr@T?j>2asJ zF6O|hw=^~Cdzu>pL^TCU5J?#qT6YOTF8nG_yuPz9N!xp4*XNMU=hijyAN7Zr=vHE- z8;jqW<631?5tBg*P}IJ|j531CttKGNRhKdFh01Kt-H^`os;ADH!Y7}?B!HFcOnjf| zl)ZfLg!=`sC$2)1ya50*bLeTU+?+=3w+;pwm~D+{alE#WEsZUvwPyT21;eotbR(S6 zmbFvBDrh4W586HY5~P5jE98*U*PD^b@jeO_eC|4UC->Rr12;^HNZMbe0A9rW z?P!XST&WppJH?6cb&aUK9kM|e`W`eoOpxO~TwAKJ8%y_IWhb}+HRLe9P&n(Uc4`PC zSRSS^!#MIsOogziK_S&pJM_KAc>-I}qMMs2fnMexh33Q&0`sig8`gUwl#gztV> z%V@F_)x+IL|&cOX{Ori!2Wefj)&QO`N_l^QrqZ^l!F)adN`H$xJ zhA}O^iy@EbT=VGSj)0F-+-xT(6LnpnLPo43d^e?UNKeZo1?m_azJgL>_MCd=eoQMSK zKP$YK`JCW!1{JNmoz{f#A@=UT;5$H@YKpXMkjH>PEZ%o_gXU&;j|_c+Hw%fO3AxFR z#CcEZay+GVBO4m|v7KSzc{VY^6K8=M+}Y643!BI>;B8OTsuWVPl-IyVD^N((h}9w@ zbUZ#Z$a>LF#BJU_MO%0g3ndh&{cuQ}!d@;PO9S6N5yeY+aee$m0BRH9tR`|Pr;^Y= zFE(4x@7}Sh&M?0#6y}ULt*L(PPsii=PDY%bT+~=*OOIv4*c_v4%qGpV&Q{zaOQduI z@exL3l@TcncNe047ovF+qPg{(6P+xREXg=q7DOai)x_xOkcmY%=leW#z?5biXCn$R z8D=AbvehL%R)8G{#o+Rb*+rB%R6NMmWK*yUd~+sjQ4t)!r1y12? ztEwkX6dYL{i2~0EeAD=O`QYjxoGzhx`4pu(E(L}dx({698Xa7))Kxmpm21<)lzj#x zjxSWH@6X(_i}P^mAXn;EMWI z%0H?|Wg6mrU=O+>`rt{CyM^o*Zj)_`87&13gfs+Ay+b=NR z>zIxAso1FJ`Fn)x}qz zABTgMyMN}hSH*ZmN@2b?CNL!$nJzkIEEZWj#TFaNzPFbj$_3FQpQ3w>Hz5O7Ku(R8 zU#pQ1W$h)dV2zlYPk)y$_C-Sf2_@!g|07mi{zlM#{IL2oT?$`H>&s;VHxUA!xi@*P zq|QwWrq-NY7N&kjyR$!1GPx+#1i3C6zQ@p{@zu8(DJ62Z3Go=?X4yRb(dZT9)b>R1 zJ25qFvANEmShDlW@oB#Ebb2jwWBrBGGAGcPkgDf?WhY&a9ev%of%94A4ruq>y01*| zhPo)FWl(_f*4n;u`e#=&p!OPswtK0|+2eeNA$4kc52DECF`v=?plQ>68KU-G!{TTaiY@fJFx$g3yg?7v+SIo*9W0fl6$?-G&fu z{sMfdj#aqHH(Jfx6As>JF%2L?95mRt!};b)VV7R+MN0O|^vVdhI465xomQFmv1oWz zr#Ga>pZBQ-J}#OE_vA5|7}&$kk)WWRFwb%gIJJ_10nIf;?N(yw+_)Xd=BpJwwLyk} z4<~{Vk=Gwr)mT3HgbSTD$H>pPFGL-dY?xL+2z|R?jEw(4U5&Qu%s@~8Scy~2*D=~` z!Q}cN2KK?CN!LC>j$+3S;i8k%DxwxbotPQ~7YR22`kMs!RzY+La2j%2V8NtZSE_iw z7BDU{E=U?u+Is6vU;A#`hCj(?`_RJxTtqhrI3w?{mhuYN{95-nTv&mm|6={$xov<= z=~>wRF#ccP6!nY$eWkDk^=)-+jIHb}ZC>fzpRBf&jyZsQ7Eu;eQdOqlvoY2&la&AC zpsK#X!%9M^8)73gB`XSbk%uMeTIV zjCJ`e49)cMXn&r@XRG_laI?`e{Vcxr>a|QoPe=DtAfRLQQQz3m2ta}Zs`)9Bv(qDz({8u*mm7o880{&@(U)9-vw)@Gt%LDlC*JD4W{5rP! zuhsp#J&bK^?F5W;YBBng)9JX<1)4|#G_F%w&1g{HU3@j zvsoo$Jv$@Y-`M?Ey8Pe%(X+e;5->3_|Bt@^t%m;>ve*9nj6(p#6Id!(7{9)B4DaWi z|JMsa_j_Q*U*20r_eb=^zdj278W+GwNB{S*fFmUhGb|;PW{2yZ1CQKi(I~PtsczHp zr3ST{Kn)UWB!_NO61vKRv{ngVJ7QA~Fsyw)s09D7TD)UoTb@vmzSWzMSNQyfWT+Dizz>D5+L`zf zLZJ@5_;F3zk`fYs!H0HOvH3ZO{cPhVil0aunJ}0Fjr!SSYeXbn3pvC7uKk<^H8@Vq zBcS@=7CsiW{caf7twkh2>JqI(oc7QyBNO+Y_3}A8xTMgaIJkv(=Zt=UMuDG>6vw|A zn9Ge$T?0|gTy)Nj+E$j(oOi|@-O;D z8q=!$uB8S6k`bDZJyTVa1h+K7Jj7wd#F(z$cjK7AF1=znPcv5QERU9rchlfXEA3cF z5}_ZHK4Ojt93yL#+67n!Qn678)RDK4wHP+)Hs*p&rn6@ht3A|I&nj6<4)+2b#T^_i z8SyL{d0dKo4!zk_xT|+XUNOBo;kaFnt>LNM5pmn}n0$CidpX*PyKN@c6kmErU{8{O z^4T`;pnqVm&^4PA(~z%XFpyG4`pqkOJ36Lxe8vJ!3VsMf?d$VvP+U*)Ro;XE$&fS)1dbwr~kc6;#dUW}$F;($$lvWnKuu zg;iUZ35Y`eJiaAA)(0IMBGCjZn+0((y(6neW}ySZdO03v$cy9g2^RaeC`Uv{k#PiEup0TA#wE-j5Y=iEDE;~}nF$T^4Vdu6loc%T% zUqrRkWH-;zE&iN;%16Bc&Ca$ClN561xZEKX=6Zl3riV4cZ_LihB1La`AP=6@t^#1c z`RrCBQpI(i*^}uvDlMrBwv!GyO&W9L3l@FsMU}u}C2zg7bg&m)1PPUc&FC3-s*h&5 zA(Y@!R&!|PZDNn4d_XXaZLp8lzN-z5q`cWB?XoEe9jBDUW-4>#7r}^48)wIg%;Dc5{D1=9xTreLO8Fd80$mL&6loISzE1Vn<8) zV@qeTg732sc?8K!UU!JGz&M43*UzxACi0lMj}MFshEcREq`TXC6Ie6vF^Ru9)4D%P@mMvAt1`f39N0m@duMbS0!aTWxM!$6$^skPRGU-U4>ouzR zZ-k!#?<1nvYRZJHkL!jV2zND``E#cDe%W6(t*Vyc(6?|iG?Q8(w1;0jDuTRScC!jx zP;-)B$nTX&;+?|23YO8lUR@1g9^7&n>Nz$%>18W6CT-bBCt-x$DTt4c0$t`FY^#1p z2EVUN*ioH1(5Fv070gOhwU#PS#_Rcvp2vE(($D=R?V$fh$hh^0yPErwh`)DPlY88Z z<=hp}NpQlEPe&7M4n-hmPK*DLHqOMW}MnvA#Dn83|6w^gwf= zKf5|@KyyD+i*0%7D{x^M@L1SoYM>D?@AlR^@^;_r8uDJ(E2&d?-f3vs7C6xa)6AHO zoH^qrKW*kvJ|KiiKy@)a%8-T$D0y>NG*AW#|4?IiN=wud44t9quSgJKGwmQAp<6l@ z8IgdwFWN(Tw8&_ovudzPypw}ou;zmvcE1&4jFn2E86A%s2BzGQpn^=#_9fe^hxOBG z__1ZBzMk5#-{L*<2d%Db-_GJoFYxYCnh$ZQ{9t#(V{q|p;5OgH)#nJlv4SG)K-X$tXjUx^6GN)3>Zs8VbKQ|v zg5AkPJe9@5K-deG=B#uX zUL)7`jFqOwUt4+a8m5=o@c|>)%!vPc%s#Kn3MA2u^g`g+f@KFkPu;qFKGsx&1BF7% zO@OY_c4XF2y3mg>@b;c~S&0~v1GH^COWuF6;;UtTgKPBk41XB$KLb+#?#}-7R9{8E z0D}U*h?J5ng`|$1k-oW(ow2Tt93Y(Z52fFM!#|GyZy*Fc083KQF#(VRJu?F<9z7$& zZ{8)%e~BSj*w}tH`X`2<{{<@k#1OC5{M7lMV2Gct>R)}=zfP|ILlE)b+33Fn6U@IL z(BF-ZNB4Spenrv#4juk4u*83l1ejPE{>7*N71{h}05jdc#QM^|R{TGU^<`jUWByyH zZ;{fHrMd`8TAhPK({r_H$wk@sgds=CdX5I`^;v5<66=#I=cvT-s*PxuDBcr;8L6*6 z@t?k*e%nCN!Lqz}F;1t;)YGcs?Gs6Da z18a=nt5j53)zs2Livt+cIl}8h%=ahKoen$tYTZ&8%b*g3rQmxdG{y*1SWY3!`X*I@wA#DVQ8YE(b?s)(@S+LUZ=g zaVZPB*}@9VuHr-^-c=g#Y9$M2%<|d-RuP^;St|Hsw)r^EH@6mR zlI;tvDQq|nD=m|~u5UCiO_8X_xc5MA7MX9)4k0Hl$JoD-IoPOl_Ac2-9>re15eT(Q z7yBDAy5!rv)|)5BzQz+%BoxcO8)&`vzD}R#$gzSNP{P*dzx#YyZI-4d{WNHFo4v2| z7@6sqRmW;yi%uPpGSIr8?`NXDZN->sKIp>tK||w|>qyt^zR=&35qwg$iq~4e+`=99 z!@KMcZxlE%_ok#`3dIMAOyiT+>wv2;y5PRlS8>o;i=M$SV&6(Y8W4a-859-C)qWs} zn-jBZaN<0}9#YZ6p{@;J{*>)kd`>K=f@C~Ql*gF7co}(ESUk)ax9i%y!T2p#@LtAf z0e7EXhSxL&`%B+ImVNh#KGdW+{{Xp66T4hG4hX=F-)e3T12O8!wU|atw-=!JM3b`@ zESnMTi8-B(!}2MFR|Y;(F`b*CX4kju{u__PlK8!W=l#%PBa1_|v4o*EZouJ4!l%pn zoyO7{v?L|_A1ouFLFM-C)X%%!;C=XqJ}E_x_$2US?679kG)V(SoYjF7L)PwZ?MF~? zBcj(lFs6HI_A%b0L(!#20NnUKCNlj>B@khsc5*7`$pUPJ45@c|;HYp|DB5>n_xD-5 zh@le@?K?s!PX zE}&;aB$p{;)pexetstbZ9ToU4)tIrr{#)kZcdg8owc1e81kBF+ z4I_MkAMTGr?@Syn#6vO^P4YeKh;d*GKHI_t=}48n19gqh-D$EHpCbnZUe3C>`M+03 z3Xv5t;|k8);&ze(wxA!Pwy5LNh&oa{K(Agj!DL+pyVzt6E@u*i@SkEs=Q8EgtbeSWk63ahvYAu zr=qdl9?*Y+!*WIip_eJ7R|5l;w?B|BsqzDH$inv<-Jrg09&Gf^j z!q9r^{NMye$JWQ(H(GS2{EaZQ;CL9WU4?zU=j6P4cjUHMUuWVw(%Z4XE(Ucx@!s>OmkbZ9FQfgKOX=A zC~VW#>g~=dCYlzGit}DUNWP&Gx_Z0yHdM-`ffz|PH$4-@Xv>f zS{lbVt{lopFf`yjoR5?|^qgKg9Ih?oxz(hz!QIi_1jw7tcz7g> zi7=SIl%0Z_Of^$+zgecDW==zll?Z4WvZf;v?!5W>I9w0`O7q~M+&U267y_L!=dZ|* z=rH|76yU=zM@07#V@e2;t}Rv_>#ToSW!@={xA?3!=+kjN%?UR`MN@geum)l35LLk6 zontw5iU|&-IsSB2oORg06anucSmt15JOpuuOO-JS^PkS=PY2<5=aUw&{QLU`t$*PP z0K1(2J7?nWE_|Cw2Q7@>i8~bbf`Ik-5jt>Y!)B>H?tP#EM>)!Kc1wwAzXSZec!S&0CfTx+)0Cl?_L+;Dr7k!H&mL4J}6 zziqYq2MnP9lTuK^Q=L1!`7heye>BcNQOECH zr~3F8;H9SnAYlO3du1_xT3hzl?!-Sy0s3G28UF$x{^Odw?tuISA;pdLUTgXZdj8nw zXb-Rr-CssY|NDBtuU(up0yg?O0A>r1RLRoJz(&U$?=RpJKqOJqQ`1t@;u!!mB#bTW z9q}Alnej-O8L2)S+mYdsirE8*E@n!=$^sn|*)N>>$;ML8URNJbLj*0{Z_1mjXHfw#?G-r)huM-(!O#quOB>%-(+-b^eg~)_BR=z@Jhq{O$N|o z1a#|fGI|z9Ccw=5n~Z^$i4AaD{awbu`XA@gu`#|b#{5l>j*Xe&pZ%j}W?}uO9z8Sb ze~g2kg^n3;fB&sMfQ$|BMEr;BKl?(@#s(lM|Dngm#tN8ce><0fmY(gOwql@VcqR1y zp~uSd&-xhX*;xK*Q-F-^pM7Rv0Mz);`WP5qr^MgeU|^>G4;dpZ%YVxLc75$^bO4;F z4b5+K7wQ|E$slSD}Vgg`uftDE)x)GkQmghXf7?|{4QvhZncqzc z(D?ZHbRtHkw))oZM^k{UzL37Ig`Pe%4-d4Bt+hVD4B9cJR-Hd$w+^<|qrIy?soRi5 zK8qNCRnRM=1B9c7t6mt+$Oja%33v+k@$QLS678gt`wBp&hK~_7{y0DuoTK+xFwFZ}=XyU#+x~jD8ye#Nur^AUZf^f@a=Rm%aeR3d78)@OFNv@F zl_o_>r`=U+0Z>)3{X^hlxIIBhYEY6OCZxMU4_QY&Wt8;iBv*f)TpK@1XtlD?Bq(|Y zWnUw4bOR6Fk3=Z`OMfKLHlSa^3MGRCbZnyENPP#Ebp~{Cp-Gho4enE_rCB?n?8l?*_2F?lK4cL=or*@?z{5jwas+b7&HioT-Um%*P=*vL!xhXmqa8fXxMMo)daoTPLF*hbhX?!*p?X2Z3^#Hd2sb zXu5opAvRbz+yNp~P|_wYs~yoZM#Nd6ei;*m7V8Mw5UWf&x!+CfQC%QQ-kv|9X(yk(1I;?mJF^7z#gFb!7)n zuvvrr+i}kynSemMYZufXTaCl3mUxN7i$iis=7-8@m$>6aPD@}L;#wP19!XvEWha{jG|89Zj?o=gPAUS+RSASMgd%AUH!3J;G!RfT^7`(tT&fBg8eUU;c zp%)p8(dm}w8cm6ylfo^IOVw3P29;}iN&YPj3fQJdNFfI@gls60>GrqG!^@TOP3!g& z!6%||QD)TN@;v%WcvAIc+hOHCNR{G{hfE|=biZZNVi7WZM_kIN*B8fP2gFFy4=bQq zqa1RjaBIgzDQfn{8ajMw#i^)fr_lhpoL?+|kepsvVBwdT@g2_}AAPxDPm7W+hC;FF zG|ph~vZ-0DnD4TiVY6o{&5zED__=K@xu<@MzV<1CMC4+a1O)foCrHwLVYs{Rs zcqyYHo?gqGY9o*-9d>zhh7T_v)n4#ZmUteyuiF?%L200Nkk3brRVG4xTt0YADugUG z-6(2EC&iFNsS@?m&cOTmk2qo4=!wF5kgQ!y2R$X%3+6!#+A&RzAA*Ks?9=9@C15^C9ALo{mUcPaO;YN^#opOnK ze8VJ<877@me@nEI+)Y#m0f1i z!>>GbERI-2OS&1c7>=e!&IRVC0QAF9l3L^&MJ?unz4+zmp!}GpJs|wCRKqpF&nkD=e|&kc>l4JT z9ENAAi_K=n)vU}`mn}$Xy|)k!chr$BE5mN%_>A2I--x50-T8Hq1rzsp%M1Bftr2Cx zy*@%5N%3XXK5`YojT-@xjJ1XzD=PvbE7hG%2ybmmD^??3k4d)m>2e*<$0y>(gO|r+ z$?2DAEAO)D{uW@znH0~{j{UkZV~gt3ME|GBv5WV!&n~JF(g=<={n9tOWDSr(6q6dE z{8L}t0Y94f0NiH|gQRka_EZe?T9y(0)*(0*^_3iZv{!M%1Q_AQW^7)JwW|p*Qt0F` ziml<;iL~#&qm1vq3B#e66f_&IR~&cE-|vu)mJM@xcoG6 zy3gi@#)HabR)VOHm+;=G1&4S@Laz8%_I~6&=tnT&M@}XGxT=UD*151yRDoG+tqcYk znXc|u*0%YtPDDdV!6_H;uRB2=3K{lMi(pHfXG>TS;YQv#q)A_iak{3X*7^<5TP7s~ zs(Co@Ow@>Ky60rCR*0K9Bp{{D^pTaA=EvmhVx{>STXi2@O|(`zY`t2r&znlIB55%q z=lWsA&|$<(w!w2g$1;>oVkXjJOqZ0#3zTdD8-wf*Pmd$S5ma1VrED)k{-53 zg4hLt)czFnY7+8#9DV^lL^4w{`G!&?IivJa?D<{Ca@VOa>0xyuPkuOb^E1;EE ztaTUgYSQy`7R&^xlVmPvy;sS2( znLqUZ(%S^EtvyXF)ns!~YCTPqS6k+8s<9vWroknGI`6+&Wa3kBxrC_UW;RyL&!EeidoM^2|rHsV&lSwWZCI;acnp(FaP3aEmc3LK zK6j#dlk$7G=A1d9&7Y%)I@duA7OsJH{jp{Al_b5b;X!3flj(Mw|Ni?JyUV327~M>g z=-x!2>9dv_^&E%I0I*4h{n)SVxUGA(IVPQQCSOX_n5}qPxJ;N;o9{8pqT;JF zGl)u|%G&di+w&0GBy6W#cMOB6dPmp3t_h>(33jx3M1rB}U7ej`YALiQp=e24zO4RE z-j}$!S@11$WYZGYuIjtXr-ExLYdbe_6}DogZ2PqyCkt8*6@ubjo_-q3hun2<-bIhI zca>ngCWLV`#qsv50+BTR^G_|B@$;)8iGy_U6y0AR!?t2lzFzRQr1HLff~`FvRd?AhT+v)fv`1JDgX=vQo!kWb!L{rwREwsM5-apj#`?n9B{vpmk zj_>lEk%NKd!{MXI_)i)CN4N_++KMXLzKeEzZtf4m`wty}{=I-tC&-V_fKR6jc>jEt z^6wRYebR|onA`rF_fHyn8U}opzfs_$fuA%CH1zna|7oIsH?jR~dT;9e_a9wi|7Vx* z>6Dx-_3`QC00#Phx+e#))_=bm!{6u8N$cww0R$`@@ipE}^!SWy?D!1KjQCoAyanHT z_U>+j|G_GMfr+TKg`MTU_4cp!|5*PXuHySutpVmXmLJ#Bb^6;Zj&Eyi_oo*hLlS)V z64JLf($!ZG75H%ZM=%uhZ7l4pb@grVKLV@pZ=>RPFaHl-V)*Ft{|-za*ZGUA81R`H z*cktVpiWe@OtDo^+Z=BDk38~U#gi#EWco~|R+}~I#WlICkQ{!_%NQHZ+k-HPXI87Z zAO{Ps$|-MTuA|Gn%U*5 zJ{nD5l|V(JJ6>{WvpI7Z|GKwgbC!yz069q-FDz9sw0Jp%_x8MdB{0w`fI&t{*Pa9S zcKuB=Z~?vpvy?WK=P5b2{WQQ`Fz~i_I%KV^P1ChJxQ*(*-STH9_u9)A+oSJ(pCf}c zI4l?Lsf6HnUmnNZW68bcSHf5D)Dt`EHaUvVG--bBP|k;@56{skIOv&gpiSzW^%b@z zhXWRwZ*Y%ye!#rKrW|U1^)!SxGK{o-QNyQ;Kg}dHmP>XL?;RHrObFFp=6iR9l=16s z?!qd2Sk7aL6$Uixl$WM3G_Ksq>B!35X3`tLx+|Wt?iHW;OskCMe>6`mK@>lg+pupO^(KD4NaD$N<>kFHa z(D8+Awy}y3g^GYMox80OC*8Ou3k_W`aTrb1hsc;xpiAHY5FK4GfR2J{5>)y2!FJ}zw$u+$0;NXjebJ=f0b;Gl;CA^8u zeftO(fCC#^`LUcs{AU#XBVU~7=^(u5PpUzxP1}Ac+(Jh)6%E{GeN_ympOSCfalM=N z(D0(9I^QU2ezz)VQ2!toDi1SlB!(*$qsXp3oFzmNxkwaUt;-Ic3HLVb?(xb7k$RJn zjp+Hh&R&b}il3Bgqh^CXH50JlT+MMfpQ*>;$z8ot){QX?VTWPx5DAA&?XXWjI5h+K zJRD$`sPn*qQ&+W82xj+_*M}%QfQSpg*NGU9h7D-^99~hg0XfXi-o*$jstWgR{q9Lb z?a+-FA7RLXJF6(7yOhS)cB>k+Nfqq6RTY{c@|zP&FBdWew|cSkqT<33?wCo6!$O%$ z&hc6I^OL17Ovd8{x_rUxDDBK;ox&E&0%qfC>7^YON`Q0)3Rp`L+0Pbb!QR+LW0+q1 zdpJBx3hY?>o;tG1+7*{O@s0qg7ZkRYAHVYwHPS@uH{zw&J>e5w`eQecXpvlyObm4W z#e2XIM?qYYrF7;as@U!6tkd}jNKUJvL$KvJtwGvz9}(?cXsI8NQ$G2{_s9`J+_Eb$ zVX`PfCyNw{y9ASM&0>Zq^Z_$7QieM@v4l%Td^KmLywxY#8Dy@>cF-=q@mngf#-h@- zc{sIjTsNnd?zlT`B`~Mn0X8efU;o19~QDvpiDp1~_fA@fpS6tTKzqrcKO zY1Qwn;odD85L1cRK765JPXrvs{)+$|!P>*kX7jL*=Xi8lcOj&gW{J#dZm3P$8QW1} zT04*^rIzrrh}25fZ2O>K*E)yDQ}zIrAPqDTJQAj@Kt*rnlEkxM$7q)wXq-NPKZpb7 zXJYq{c#Z&|-Nf#73>yq=12`W#CIgUlsu93CLx9&fHaJ=07D$R7oR2|L_j>UN;GAU_ z8{9FR0|sQ5o}4M5lR{dSd0G?T2M07NYKUz}P((%!k*UNI#=4lV7PMeM*mGPi8dx}P zf@z}A-YFrZ6^fuZE)rO{HdE6s?&cL2*n~om*M>ommo{BoYNfk1I54PS?4{$Ad*zzG zE%cJA_PlwR+ql=%|3-l;Aw1uWm-aC0LcR93Sj;si@uqsh=JMTZY00&_t3l7sEP{Q? zqAZT*F=MUetd-h7;Y43oD)1h)B?kHe|98Y^L@Lr49nqdn5(598de7(E-U2ZAXSSag ze8k_oN#uD$;sv5Cztu@b0V&G?~h}EvYO+Y^0@TtfWx7_?-aYSo_}S zeY>AG5KQJ3h6Vb_XqC>%A)s+*US`8%j|Vj4j@Gj8iyIoC1)d?s*-p>?Pr3I;Tl}r} z=vkTn7FP@(+UbAlkiWFbhu-*qTDaAUY&Y-0hv{=j>)Nto@64qNR?50Dk;D{4fb}N1go(6@HeDcn7d#C|D zMtf8MAIlZD;@a)^IM|b->uFGi$<0oM9?{gVP!d)EX=dCZKiB+c==<_{oSyCDdO#+Z zN~}a2{xwlui|^1Xik-mA;p(r|#E{~;h@bhG7hb1Sk_^g4so$IJto~ap{}##b5n^QI z_-82peF^+8A^D+&h4gK7t&J>gEv!Fu;U9S{12B76vFhqFvJx^>(g0gSeKUZqkuE^- zU1t8%Nd#bKWa@-Z{y*v|{&Fw{7`*FP#(#PX{3Uv67~dzFhJl5N4WE&fi5;JjnThou z>Q~$rU}~hxZ*E|!k5B*SJboM95B1E(!SSc~(b11G4I?AVAA=ykQcT~-z|i(jGk;8q zw)$q@@Y(-5|Dz!Gcj05G?{|wnBmRf7{dcANq3Zv>!~PYPznZhtvHhdymEN`VN8KN* z0Ki874@LanhRDd;##Yb}VEuO#qyYaYvoic?8Jfq=u~y3nk*M`j(Z-c zxzzJTD-WMGPsULMQ48^^?P!u+u`6>|*h{#xBdf!As0tZ+m|=XseOzWgJGA;WDg6kd zmdh0CN{a6G+%j-p$*`{|u*yoXd2XdfVG>^Cv2U$`-SNHU9!4uF%)goFvU8dtUUejgBQ&K$p^||d4`b@qb!^+vNs3p(jcwxKhyLtB8+=RR-DX!;Z1KFWZ zCrtGb5o?@gCtuxMp6^eX>KOT=`8)?XfMsyo-Il6k)-dtf-JI#vOWwrbivr3hl;Ed? zgrzT=u!r4Ka|VkfNV#Mf1o2L)is8d3rj)}OdWDyEj#rPhNJiF%g>B?tf0uk)lzfh8 z@zNE{;jTEvg^0SJ9IJ(%b)I1$c6igg867;w39;HDMGX2{H3=giE6KG zm-S+8lPKXET_LY~1H2~AN(@-p?%XH|6`4(P`PJHvdKnc7JjBfBLomy;(PzchY_LVf zx;7$(Ar7i7Xq6`A6x`D9d%wOdcC(9SqULHYhEF6^>P!A&qrKdh6l@{SaCT*>9hIg@ zc17j>UMM~5C4HFMv?98uh!wUJSze~ad~43ql5dSSZ`4^|B}zFIR-t}vAJwF@^;x>u zd#HX*NA~3zq}`&5?ZzbZ;0_UI7^fiZ*(KfGVVCE)lQuE$M>N*6byzAkYSPzuCGWy?r=$c}cr3wJK0Zc`s>DeSzZKm{p zPDBK?v!4`A9L1@5u+MO6KtwvyD}AMiOFWM&UIglpjBCv{w#5r;G(V|yB3a|yloRcS z`H~+6_oqRea?)};Wb*W`+gZ|h80axg18JPIg64g`Z*6F5?4}Hfh}uF zR|v6<#)*-dWtvZu@MK?5p8xAZuf90Us#PH?=z3UC=)jMk2C29s$5qi#9!2wYFj#Il z%`~?66q#K`fWqV|a+1)ylGmBdP-I6~yNYRtNECGX?nF8c zu>Rt=GFtTa^a4~bq@U#zEZl8ZMv1d+p84NhD16QKgwve`p9vP93V5J{rS{smQ>`psM!)c}H!YGI^RnDqie;db6D6K!dShlzz`!^m zFnA5??XYyeVCkA|aDf(pQ`?x=fU|d|j5&!K-p!xy5Me+<3LgQo78q3E zgc4Ms)Qhj6R1|*(91*g%`2C&!{=JO<{mmkr+?l=_H^I!}6 z6pU=CDr523#nDG%bQK_>;E}RF*D_TIpkE4Xm3cp)srC$f*TRS=ci;$&faMZ#aR3$$20<>Lh0vv^JE;W zbJ{?kdDcCd5eEymlve!w)mSXS&oWJ;u{4=5XnCXFhUvQ*)QDd|YIV&6MvEkqhdaJS ztVeF=yWl2((e95|pf?m}Myc4cv^N8{3mO$}92!5$X-($Gxj?Dm2vg zywQ`5wO7GOIKM^Rm%eepS_Pr;UPKrp*Gp=`Us!BRYJA5)io47>u6VYub) zkO7ie|GbOh@RbgzA&m}+4EmhV3X(Wb@CUjT9U z%~LWA4M6CxD!cI3beIIU4&k2Hs*phEiP+@5g6{;~xwzJd^6<@p1>aeA;yDfx88r_R zLjoUz-(6->Xy;-V*|H7AoJDig)7kmQW9YNbrb$BGE|ag;oExi7@Jy5D-MCkrFCL#} zem7noUTSV;8>P2DCZDodrrAIy5=P?W57hDM<*cvSJ-Huj1V*-ci8aNOhgk%1m)tuy%0g z{;^VDXSmf#Y0D~#%vrk8D=5S( zBx}B+CPaDeVCHI2i`_LQK0A*qbyiR7N`}77P6}>|1p%dHdwE)4TYJtjZroMWxO_8n zesIA9h?6AJ=f8m3hyME;)R@`-1!^Db<$ufQe{8T}_~(N8f7OcbYtaAPVWSl#Y2r(d z;&=Ce>cF7|GszwS7bus9h3+P|j2GEI3vLXK_b*4T}LDR*x@2GvORyu?#KCbE78 zn8Fd+FJ%JWGlIh8-}r$;1%^?)6?XYOVu~rPYn4+JTVr0fq|D2M*eV?0w;Gd#!b1H0Y=5iR8>J{{TKCf&@oJx4x{gQq_$*Q22#J=s zEVuz33-|%Jo?X0wx+ju6oH5DzmCD=o2W}Y`d*#WVLfCn6QeNy2Sq%x?}!l zi~m1@yFbYC9}M&FiH!dZ?*zq#K7#d8{TH(S3-SJM;`Dx&qLa-#WQv>rwD<_d2PP@# z8yLOc)$d@;|DLeY|AWIHP$X-ur*HiMb>!lD?-Xlf>-5hmMLSDNQ+=}!w*3RH|H#e! z=k(z-d`$3P*%uXj4f=NgWTAggNilHXvwj?zng1OBeE)}ymGLhd6C=kvgmU0Bzx!xu z{acRb???Hg-~U@^KJ0(t^#7}v{Ee#kbbk+%M(odh@)y190p2P6J<=ci`rga;0sWQw z`ZoyuEARHN^cv#_!~TyyER6rksQtGUpnr!@26~SFk&4=Ja&lD?t~)Pmy1Q*#IVDb& z$vqQJJnppi)atQT%3~*^oLaA?t8r`G zckU@4yt^&M-XyT4(NcBLCrOH4Q5u1al`_C9S*`i)8*H;gayd_Jd$>~J!lpYBi>%l= zvL#_A4b&xRdXNZjTSZh_fXWcQ;wM9~>>nb_W^zk?-v-hGMheU+H?26P$%wOs(~@yo zl2PoCW@QycwMRtV9Y?3LC;uroF{x%XrXxK+)3-xVY~h`v!nZAzssi;EWni}K%^P_q&f`4;(y}5 z>R%7jWoh1U0Z6OPBJIf7obi?!FYNq5j~QX(%H4mFf>_WwP~7r@7aiwIZh}iwF2nZ?7|_gd>ylr3m7Jya+ewda z-P3AazTC&D@PSk8$gB*F%7hF3#yQ9UURp1c|p(zT-RO8oU5~c&eY&knIwC) zJXY0U4^ZM6K54dIuHMMEBEmZ;QNCsY8P^9{(njA0rm?PHB38(WM-Gk-&Iy{cQlMc| z{;JG0&8hq=Kj}=fRBem9^r7RYd;p=ra+A3BppUOCV{f zsfkGRpcY@2%I3>m)h`5`_6JK{ZZi%iOu~g%Nit87IzJ~_ij0j`MxE9H>Hux!uxt9z zOHx%R@zEsD+@qQ=U;Q8YdV`KKqfFQ54ZrQHe3`?*6LamirXE{Dja%4^#fA}A02e%n zjxd^=)nfD+bLo1pQ8 z7#yx#W2{}z;J#lC_gJ>B*=;6~ z2fm6km@uTv9;bmGRK1KL67CP8BwoI0YbZu!{k$`=L1t7Yz=FBzX`wOF#q^$vbH%D zH_PcH9pmoB9jxl+`=-AO3DWPMF#G18fhsz=tW!D)H}`27h}{R>kR805QSvu0rea(U z!WKh>s+3*F-RE#CucAqRqgYgPg0aXihH_>P{*ag)Zhg`dy$1DojIt-KBx;tmpK$M# zIXdP#=6?ckS?n2}fof0&l_5XSWdQBhuK83uf<;J$mSr)$>9hs5lJZwX+b7AIs*`tc1UZ=&Nk3*R8yEfcRNl^YErn?s)Y(x4}J_`&J43o|3gOcvA z!q8@)?>{sRjp>cAl_s{oJS9pPK?TsNF?F#|})k5He8RvsQKlwwIiwVRbAg~h%In=5@CM>^*7krxW#>sAn6O9g} zC0?^CwJZZUE{Joz@ZUz@;>7e#d-cy1#d2Qy(p8W_T!5#6WE+;IdP|8xn84u50n2QS z=_bQL0fkv{f@Z3#vdNn)&G`+hnSgg;Z;&b1DOSV!yLlP~l?uUKp03$3>2<;yQ$5J8 zPJ>rGziE)0br6Jff@6?Al{kTe-ct0r1aX(wx$oBv#>|&XS$J5S8w`h?`p)O&1e^Qn z1d=+3lSGydsO#Ips?=Sd@LiQIY*jcH&o!?TuKWCYN|DS4xtgw&-%uxHGATdY$b{N< z&Gw6*EhXB7>6b!#`&dNU7Ob1~`5D>zSS3@w@XoTXHITc9RWi71T~>>v;#-ZJ9o>Lu z2dN;JPx7fYkZ%2e2*hdHn%k-Q0t403n#%J@+;KFIxFH4P?9jYsPDR3om=bK}hG!@3 z2Jt+I=WsdJydSH9c=ZN{=EOtA2 z=k@vDLn5C}mNC7l61Z{Rh_xDQhlZ4kb6PE=ERxh89A)Q=;c?;Ko>_{#T8X&3Y%r)i z=o~1ls1K(W^d(oWtVCh*Ve!aJYOL%P;UV%zQ5bJdFyL$oq!?rpeuRZ^nKSk$Z*EHe zBr;YCXFAY?sxk@)`X%*~!;6@hy%YnyU1xWYBF_vaT%YdD8j~h&6Wgk>L4I3!rbKJ9 zB-K-^M)8Ss&^dD3Gl|>#X{x@r=I)O9xJ|8l-YzuU$B+4VZnU5RKM-fL5nNif-E#w< zT;7IxBuT?6yGXv>YX2wu3su(#R|9^sT#P9IBiJo-L%YaG^B169n|i&aQAFj3a4?g+9f;*$vr_AfItU0XdNhQG#W zRDabUpH}5w?(0$;uI(5i)F{x8c?DFI;by0c@lYrplC^FhU`YQpRrc_<9>-Mxv-tO=^Id}+` z2e-DxrMPqsF!I<2li$o_I|99KR5pHN$w?52u#h5jJX^`>2zO?UO zd0cnjlR|ZRwGjv1dEvaQvq7#gVy7VjF1{mmU%TkR^#c7sXO8+Z0UEGwXASIf;`K5C zZn2qM1c@6!a!xNYLzx)Ja`Xd+6!-{_R6L4xz%zf6G_9|h1JOJ{-%)Wn@arXx`lazH z%B_0X?G!3B07ZKqK5`@KfgWn?{7#S-CZIVp9F$N$3W6w^XuPAkBHe{)Drj!klOj<@ z$Q|Jo_W^XOvY)5|b}VvqsdazUA^^ym2l%eTzRa5}>L6-O4g`KBgz=mnL?*I(CL&wb z$09J5dINle+#?i}NHl5~3RUjYa%SpwB*6@kRB&iuN{5O7Cv}f}J^$z9ZWM2o4I&)` zTA{QcO;2JmA3#GIt78K8ZdWQ*7nsy!2Z$H9Y0@VTZ?^TXp@!PoHlRL{ei2E%QSkV< z&Yh02w@$|oEDwIP-8Cs(eemHhJy+G3xH?Pzyxn^}{ET80Nqo)ODjw0z5VCTZws_@q@7cp)#e6rRa+uyh8ikp1LMWh;tcpPEUJm&)uRVbsJtD;j#zfbWy@^Xu` zujB4IT%8pcTg^IjnK}Y#`+g6MyC%sV#K{5wJ`KA~P!*Pfe->q7M1Uh|a^O_vpFG?; z1mD9PU%q#SJ_mWEyyidLCED!)^@cPUT;erp2u~&VD7r4Lx#-Rgg{Q;$={thgeh^fQ z5Vsz8Bm6YuCD5X{ebj0WtbtR6D}NgiZ&doAe(;!Yks;&FL|J0XHLf$ATOyB)td4WU ztIz8YiM197Z$kHu&CUvTY?0TpuMyoH$Jx5LFwhKOK(oJ4XWRu^Z<_ zj}Y5TwwDnXPiXKI54QYZcN!qZ#QE@D%&x-Q1(qiLznWl6*GwWMjlid{I?CBiqT+<}bD*yhVc;!KtUO^uFshD$LZBy5KQ|;b%`P8tb z3$5jzWqnoYQZfe>mdtG)P;ZK@nM46(7V)0Qigbbf4rt;LNzY0!^e_jdAkYPno}l7L zqmH+GlBe^6;5_qdE;qAZ3dDI)XYVP48xrtcS&s4AB(aTC$1|&2#>QNl$_BtYJ&YALY;fibXwbO^|MR`~wuAje}_-LL=HsBc_ ztYv;0zfgCKmD0^;&l>65oZHh@}P`;>-Bz#)!kt2(0i>y}uf7 zsn!BnuB_ z_knj%#(byh>B%`pB8C%Oop4_FYSd_C_-h716{zS6*!ND61>tW2XFi-QF}yw1Uo|nL zWTAbC4r-#@#rmgc*v{yTxlq$_)^we}fOi9J05M;B_*BXvEFwHmso&A;1L`~V8eqpZ zC_Os9_*(~LIye%Wg~RL~TDBVeu-$$LOy18RW&A9$fzMrxmZ5%<@cHwOH*o##-C4ps z8yi`(!?*I0(6j@iSKL*LgI!7K3;P=}a}{)hDx$%(!16%$(6ABGZ2t0T%5l4JBr&R6 zxUy*Pwj%tqD>$5AV?)1YYv6lqA zU!^H;B|4FXwJC&1!)XN{y$CHDQI6&I5usbfGsDDOxEJ^->*i6!`4vb|At2}qjHt^f zyEcB^H1YrDpSC?PnkbnjDiozm+%GjD4f+%)FwwK=K8nIDfm(;_V8jh~>d z`IUv;*jCL-zN%j_Ee(y>C(n@#qLYi>@Ro|2UMz{U024)tYFdnG!T3k(Op1Mjm7tjK zpt5_n(%H3|7Iif4>S# z<0~~CT|YwwoivW-EIUU2Z%+y8FH{?*wsuZxktKZ;By~u2Mf;C=P9bX2vr-fVH*nW; z%wr{8rB=+~fDuC3kZyvre)iw;)GXmtSI*vGgTF5i2@PzxnZwm)ei2KGv(@M1vxLMK z=p;+;VgTrWi%B)l<;g*FWlQvFPE!MroB}fWL31gPBouagS#I8!HuVj(3RtB%7ClT`r&EZ|7ae zhG0Ux!Cz5~wenOL6PMLM$kMh*B_|NL=_H>qvvX>PQEKn+W_5HE^B+`&1Q$sV9!8Rq zF$d=8{*IfYT(WaM8}EBjLSk9Fqkcz8r*55m=_>|l>^krMa}rHObDB=f^iBGz#?WO@ z)nafzNMFH9S6Ql_WXqgVc=<`&Bt*WmfnQ&fF-)F0Nl|^0oJ>piryrwWq-<}x1jAt= zK`nTxLLtG_dWI?s^RI*>?zA^RzT7s*QEtZ3s*Vo$(}M-=U3;NxjHAB@99W-Zh2Ba} zpUA4TnW6C{R6s9+I(I15m6(q=J+ACe%GdMW#RPsQYdMiPSREKckI5p{3(~%A46Sgy z>apvZ&i#E%EC`p}&M@V$*-Wt?RVqlWVNVA~*EfAln?=R`$iX$`vT$&}<-(ha2wtY* zzUF#7KxkdeMoZYdj{n-;w=vmESaqgvxtoh;^H`qd&^pfl0(V?19E*GW0)7a*5k)WN zyLRjA67r&RGtVpUr|kMQRCGYp5T5#f<%l-4$~VUsI?O;Ybmo=o~84jC7j`Z1HUwwNPjT`Wy`Z zM-TSHzo)GM;C$bXNiyz zk~-5l}q zK;4|Wa~lF!20nKPqWn$VJbA*Tb=Q z#X=^dl?$>FopDRn^LQu14L_e4vZz(1QMQD+1G^x=5m;V-mP(%<2XK51Sz^}J_B=;< z?_mt&9-`7Ubske9Pj))Ne14>toTpsSpQjQ2;5X41yOW zDKViKWHv$(!Xbi^SgOb&GbjQp!7(2sVL`>X(XkEY6Sk-$#0Ik-;lR5W^K~m)3aOB6 z-932IXEO|%b}u>=wy*qaO6OeHiQa2@9^cM20<{7orgdp(*-#EEaVjDLRfST#PC(pB zS$Yq{QPlSdlp;PR1>_&oLg^f6Zte~IKz&uZ;kz!@oS!wZrR#c^uj5z`OGK4iw{NYeWXSr-#W)pbaC&9#Yp+ItyFG)FN z^1H<7-C7zhC_=+upY%s5uq7gX^b<9}imtmC_BuMPw~ZclX*M2v zF!IFhcDR;r($E(S+e?a7z}MO!$fj7 zSd>2fa=f*`=52o7RL;6&lSvsU%!b!lC(jE1uCx#SbIlTU~{S#j{r>rwg(X_KU-f91#aE3-f;PZ>GD%m7HTn+-#y0X*UdZc@Vq&O<1M0ms3bg}n z)#HvS@C=cUz{&fT*L1(fuY>~+(oSqq^?1^$6M>jW&`HthNjv7L#sqz}*Ok82fB8Ci zqNh?K-eVl%fPcR?TAmJnyZJ?U^7~hJRTc{#MMvhl`L|~H<=~4|btPW-4B`b{Q}GI* zt>>S$m?p&w(cUj*a|buQ80Ml>E;{crL*LF`$0Gnc-&cX1GvWql z8Q~9y8#9Xau!U9VX>W@R`z?dC;BR&sI-EDjZdCf^di=>@OcTL6a1TG(YyrFi-j%+ zB~X!W_i8s7KbXI{XznJfKpM1QPz58^$Zxz#)c*n}t0F(0Ct06av_ex5xc%ZIc}``T zcH`Pd1`M>=Jii~-+OhjIQ(fH7ap!})1k`A#F~QzvUGC$8Me0i(RG%7dg*;~#`=ak! z3YRH$=G~t2Gb{uPyED6kQYf)4T5DnrmUs}> zY@=Ii7S{w45xj}@Q1z%BLd%l9)|7Yp0Rz47Jp$ont1*PqyLn#G;1Wa;d}%GA(?lU=&H`>65y)$kZTa5uSo z@aMjvsbS~yq2yBG;&l{=okM+EHZ1(A|(L`n4n&_Hw}XEO+B zAe9SYElOA-6$C-`w+Cs^DV@eOSUkBlIgwZ0E4gqc}!e!cZghoYbZd!Gd^Xh`-4(VsK zKtLmtn1kb@fTKa_Id8Gz+t*Ntje!cxYFV7cr~)3R-DRgkg~iv!v&e?`@n4eM}6Z5S3Cwdojx z#!XzEwgmh@F4S$as2R0kMyz;$M6fVjZLzL%NPoADQW~2dr6%l}h@WM_KTE_0FOR5= z87wv~$3_pE8Ji0?^w=2Ops_Ky2~qE!3&F&S&Y74NkLH7zTx^N1br4JAv9^{jgq+=e zunyR$5t|0@*p98{+kjj7^1!uo?#+2oFG83#>v;fdDV3La7UDT1caq$jH)^(DlNCm| N#&IH*n#fNk{sJEF9mN0u literal 0 HcmV?d00001 diff --git a/triangle_etmt_plots.m b/triangle_etmt_plots.m new file mode 100644 index 0000000000000000000000000000000000000000..2c341f6f6d0dd3c1968f34a290475808e17af321 --- /dev/null +++ b/triangle_etmt_plots.m @@ -0,0 +1,55 @@ +(* ::Package:: *) + +Needs["ErrorBarPlots`"] + + +gsraw=Import[NotebookDirectory[]<>"data/graphdata_etmt_partial.m"]; +(* gsraw=SortBy[gsraw,{#[[1,1]]&,#[[1,2]]&}]; (* Sort by n and then by tau. The {} forces a *stable* sort because otherwise Mathematica sorts also on triangle count and other things. *) *) + + +gdata=GatherBy[gsraw,{#[[1,2]]&,#[[1,1]]&}]; +(* Data format: *) +(* gdata[[ tau index, n index, run index , datatype index ]] *) +(* datatype index: +1: {n,tau} +2: etmt +*) +nlabels=Map["n = "<>ToString[#]&,gdata[[1,All,1,1,1]]]; +taulabels=Map["\[Tau] = "<>ToString[#]&,gdata[[All,1,1,1,2]]]; + + +histograms=Map[Histogram[#[[All,2]]]&,gdata,{2}]; + + +histogramsWithLine=Map[Histogram[#[[All,2]],Epilog->Line[{{Mean[N[#[[All,2]]]]+StandardDeviation[N[#[[All,2]]]],0},{Mean[N[#[[All,2]]]]+StandardDeviation[N[#[[All,2]]]],500}}]]&,gdata,{2}]; + + +TableForm[histograms,TableHeadings->{taulabels,nlabels}] + + +gdataSwitched=Transpose[gdata]; +tauChoices={1,3,5,7,9}; +combiHistograms=Map[SmoothHistogram[#[[tauChoices,All,2]],PlotRange->All,PlotLegends->taulabels[[tauChoices]],PlotLabel->"n = "<>ToString[#[[1,1,1,1]]]]&,gdataSwitched,{1}]; + + +combiHistograms[[9]] + + +Export[NotebookDirectory[]<>"plots/ETMTdistribution.pdf",combiHistograms[[9]]] + + +mixingTimesBars=Map[{{#[[1,1,1]],Mean[#[[All,2]]]},ErrorBar[StandardDeviation[#[[All,2]]]]}&,gdata,{2}]; + + +ErrorListPlot[mixingTimesBars[[{1,2,3,5,8}]],Joined->True,PlotMarkers->Automatic,Frame->True,FrameLabel->{"n","ETMT"},PlotLegends->taulabels] + + +mixingTimesDivN=Map[{#[[1,1]],#[[1,2]]/(#[[1,1]])}&,mixingTimesBars,{2}]; +mixingTimesDivNlogN=Map[{#[[1,1]],#[[1,2]]/(#[[1,1]]*Log[#[[1,1]]])}&,mixingTimesBars,{2}]; + + +plotN=ListPlot[mixingTimesDivN,Joined->True,PlotMarkers->Automatic,Frame->True,FrameLabel->{"n","\[LeftAngleBracket]ETMT\[RightAngleBracket]/n"},PlotLegends->taulabels,ImageSize->300] +plotNlogN=ListPlot[mixingTimesDivNlogN,Joined->True,PlotMarkers->Automatic,Frame->True,FrameLabel->{"n","\[LeftAngleBracket]ETMT\[RightAngleBracket]/(n log n)"},PlotLegends->taulabels,ImageSize->300] + + +Export[NotebookDirectory[]<>"plots/ETMTdivN.pdf",plotN]