From 610c8de8a6f83a256610933cd91029d42a100451 Mon Sep 17 00:00:00 2001 From: hiro98 Date: Fri, 5 Jun 2020 09:37:25 +0200 Subject: [PATCH] various updates --- notes.org | 8 +- ...2b03d025920448ba653b470ec6492cbdd1e4a7.png | Bin 43641 -> 0 bytes ...e9d3bd60427cf20af835649efbcbaafefbb3e0.png | Bin 9429 -> 0 bytes ...0bffd38987eae852b74eb2483698169c08c4de.png | Bin 12003 -> 0 bytes ...7db0dde9ae59979f271a7cba8dfc46be3f1dd3.png | Bin 42371 -> 0 bytes prog/python/qqgg/analytical_xs.org | 16 +- prog/python/qqgg/figs/pdf/cos_theta_hist.pdf | Bin 24695 -> 24695 bytes prog/python/qqgg/figs/pdf/eta_hist.pdf | Bin 22205 -> 22131 bytes prog/python/qqgg/figs/pdf/eta_hist.pgf | 2300 ++++++++--------- prog/python/qqgg/figs/pdf/o_angle_cs_hist.pdf | Bin 25843 -> 25843 bytes prog/python/qqgg/figs/pdf/pt_hist.pdf | Bin 20809 -> 20809 bytes prog/python/qqgg/monte_carlo.py | 5 +- .../qqgg/parton_density_function_stuff.org | 31 +- prog/python/qqgg/results/pdf/overesimate.tex | 2 +- prog/python/qqgg/tangled/plot_utils.py | 16 +- .../pp_phaeno_299_port/Sherpa_Base.dat | 2 +- prog/runcards/pp_phaeno_299_port/makefile | 2 +- .../qqgg_proton/MC_DIPHOTON_PROTON.plot | 4 +- .../with_pT_and_fragmentation_and_mi/Run.dat | 2 +- sherpa | 2 +- 20 files changed, 1186 insertions(+), 1204 deletions(-) delete mode 100644 prog/python/qqgg/.ob-jupyter/542b03d025920448ba653b470ec6492cbdd1e4a7.png delete mode 100644 prog/python/qqgg/.ob-jupyter/7fe9d3bd60427cf20af835649efbcbaafefbb3e0.png delete mode 100644 prog/python/qqgg/.ob-jupyter/800bffd38987eae852b74eb2483698169c08c4de.png delete mode 100644 prog/python/qqgg/.ob-jupyter/d47db0dde9ae59979f271a7cba8dfc46be3f1dd3.png diff --git a/notes.org b/notes.org index 42ddedc..551bef3 100644 --- a/notes.org +++ b/notes.org @@ -183,10 +183,10 @@ Viele Gruesse, Frank ** TODO lab xs kuerzen ** TODO shower scale anpassen ** TODO effekt shower und kperp -** TODO y-axis a.u.! -** TODO mean, var einzeichnen -** TODO Variance of vegas weighted f! -** TODO look at xs plot -> they seem different +** DONE y-axis a.u.! +** DONE mean, var einzeichnen +** DONE Variance of vegas weighted f! +** DONE look at xs plot -> they seem different ** TODO take new sample: still bias? * Observations diff --git a/prog/python/qqgg/.ob-jupyter/542b03d025920448ba653b470ec6492cbdd1e4a7.png b/prog/python/qqgg/.ob-jupyter/542b03d025920448ba653b470ec6492cbdd1e4a7.png deleted file mode 100644 index 477a1bca36828392aef23b5d79692aebb04bf33f..0000000000000000000000000000000000000000 GIT binary patch literal 0 HcmV?d00001 literal 43641 zcma&Obx>Sg@GXi5&*1I@3_cLt-QC^Y-7UBTCj@s1!GpU8f=dW8xCD3iJKwMFd+(1| zb?Z*m=-G2>pFQ2HSFi3pv8u{47^uXkFfcF}atE%*icA+M*DuBV2JwWp7nyA_O+nWw9xi>ITV zIiMWAyA4~(75@qh3?+=5q?o2}?n$>_vXzbB zzsKjE#`V{K`Finjsqzje5kOJL1T*6K{Lo)YIPhyY+=+w(IQcPe;V=SU8b&Mx%q;^X zePix;`RPc#Re07WGs(YHG!$8^ICF>BN2sYG7u9GwoQ+R*=D!ABkN8a0R98=oO^j8K zDIN-ARhEo(`JMCeO#7V$wwA#O+#Uk2siH}&h9bM!F@3BUE1uFo6wI;LxToF zLqGg?{0;>kj{LuIOZ4c-zW5Pi|>0jahZwUDR<~;wWL5_P7 zB&MM(Ddh}usemXdab)yCQTbRvaW9^^1kDU|s^$Mp)frWuAr(j^hdsxMNEuiBHvJdJ zS^1axn6OwK?cobMJo<%x*pfN*FG@7?P+@NBwdT?+oj#IR$BeWq{Zw+zS@; zR({FxDqLl9NZRZ8y;+T91G(QG>BwcdmZ$PiPRh;7BC%{pnZ&XT}Gb=mH#MG2} z-68`nsz9(X6@psh;7#w-eN(|3)k>4XOQ+uG#5DoF+5E*L^bR zRh@7lA-oU|N|&PP{2m8?={c##0p?HN6uc1Wyd7L4=-@MSZyCQZc?7f{cm9%Ix%FTD z9(O00WL!70w6ug#bpR{&vEB$B`GG4m*VDs@XGD4-H{E>%q3fafeQdn{jKO5>s$u9^%C&A*_H|99>S_u<0EJ+nC)0}W@*ihf z%WZ1uc<@n&4Ip5Y)r2O7LXC_fygTCdSiH& zRM1hf^tiz&3!Zzi_#G=E5HHj;C1NHS*M^6s3>UjAF{4O%=0XLgg~ulf-bk`H0^;V5H()}A zxW&&`X(Q~CJv$Z!*Q?W!YbQ<_V9)9A;zC;otP#khu*&Vs38mj%@`1m(I^Lf4m&f#Y zDGy(VLqk8Jw9v2-Mm3ZL*5h%{hmveOAl1G6b5Y|h+e5ND5;km9j{^3SlF{lK!cm^w z-mopx^?&7K{3udi6~qr}GeYPhN(J`2k(xehdGH+(KU z{P-Ic>MI>*UM3kxs-o)qWtmvy%%o841dbyuP zh|k()LSp+AFb?--L3rJFtWKxb=syk;pjHtisob&V{LQ*?*++Nn1f~5O@x_tJyAAQi z8$S%$v<%E}ikbhXMqV6?2Ld)4=Csn9DkF|LZuxA>yjZp3VZ4{Qzi46NF~I|RHQJqL z`+hv3=g)>+z9~StCp=F8LaAMx86k;0Y3=bDk=-c|?70_Wm&TCW7dvAGqE`L*;AWj! zc(|U)rSEf47kzy6qobqqzFmzn10eL*-qkGF&%R#j{d<2ADWpq;QLWLog)cvbGu38F zE??T&iBQl_l@F*`tY?l_uas;zw03kn$sPg9xho-d zcXX`0zXRojDHiEVeJ!DJ~0YTlE!vj$qT@|8WDUQ1A{f!#oI#5=du z`&%AQDbPQu*(cAia?0J<*vQPnq9#ZEtLOS18VR&XvEaw1(eW6o+#7{WJibeVg`(A} zuz9p^C4MgQ#AR1UQppKYlx&_oJvbu_fil{QgS=mmHYSJ`kqzaEN{q|Qe7LGK2$Sd zi}moj0|SoBS_KP#@!3L^6Uac)RD(trQzfc)frIP*a)Y$@quXfr}yIYWz+Nh<}3)@+S-d{W5EPT;Y@M#t%M z7Zj;|KQ?I1Em?!ncUl0h00d~{;Nanjw#$Z2?lPmbtITA9^L^(Xu_yyPv4u4-{tbLA zF~l@(N;ou$4`X+i*+sDQJwr|H3m1oEQ>xZL(ch;p#R}(p@d46(t*DUMqUi#Wr1Vtt zhNZe<4>ElQuRKg65bhH~H^Gn>A#F^zG_I8H-&-OY=+wr?#~0Tv zPoKi~{S*}yFZ$p=sxmuwQI^9|oC~}ori<%7r33u(CGG^QyOo>jL%H}B)7iTA5?fD_ z6Xa3S2l)8qz}+z7$)5hcZ0BUK6qYnmcG>CUdvA~R z82obBj*#`%#Sh)4-9~W6`uZz3gI9&7qJ5AdY3)tx#@_=2a32Nh-?bruA=)}RkL&k6 zBT!CJp`2H2$m-%Z>-~QH9)5Hbr?u_9z5mcjD+=`o{P>(;K&=VL&_(y~)<2ie!&)M23O>`BcAIzI zbHjzeTOF2-l@$~D0c=tAf;i%KC)sKA>@UzBa}^ntl|_{C?OWyIv2MmNWaI75Nkm03 z0OkF6uQcYNLGzz^ts2zEsc24`UBmOqm9gp+tS*oJlt}SjYHcy(&Ffb zf9%;fH?4Qdp^-^+z1sb4WCXZO>Q!;yLX<3vvu{ma;J0fO+wQ1e5zWN^{ey0y^qZ8D z^$61^YVPBOp0i_wt1L`*ry72xIYcqZwi%MD^ ztEY}wY6?p%wlyJuHM2nwoC*C1zBBh9GS&_8aFqO~hRe>wGXh-{n6$Js!KZ_gfS-Wk z)XW1rCoXRL!J$)p>l)=H$3l6$3j+4&=ve%JDTv=niRc?d}5WLSB!3@P77Sr8aJyj-yC7 z$fA!rewMPzuqZL6*UuawC8F(eM=a5&HPoG{vDLGx7O#1}`?yw4D{xdw4=UqP?aLfQ7t0y~#Vh`ip+|991rYE4CEv=dW^<3G8LKwWwKo>dDg zsz3N-eU$(gPcM<@LFWNow>2(n_bm;i?Za1P-;Al_sDF)&wD0E!a68$ouxK9OguEZe z93i_l&zrqMQE2JMB66K!e|~0Pz#;9)C41N-A{d?D!~d$mf_O;O&o>jFfPpNh_<83W z7Gm^vEPE7%(AaQ1LEX5uZk5!(_uh@x_8c=2lm~D6BB!XAb8i|le)@xEt(ydNCIJlP z6=&xWZ;H+W-v#W8gl)3xp}KBe)zQ-wfmcf%`G$^;G3Zaq#Ke?Y{(EVO&}ZFmzg8u? zCWiX)Qe)Exi%Y4`uL3@+fjZ6kS%TZjdL>zpNq=)GhW8gK1&yAecoud?mUM3SqIS*1 zl88S^-E)OYxgoq`^zZqK=jY^dyLY%-!6Yc^yu7?{A#|ryIiBW?pZ;~QpA9? zRBBg|2wT63n;DC_Sn2Aykau9YLvEmPkftT4(gDJ4MtFt-yPpSu2=QHQ=0Lgbj|D~z z4UJWwjUdCGKtC1bOv7*MCrg@yh1OGtA-gr!dKKYUT(QJ%F>sotR*_%_rDJO}w}n}? z6#}-M#Dtlwyutij%G{jkF9c*!AxBEN4>lrF&4kI?(tt_CbO0dQ)c%#<=d06So$MTZ z!Rcc*e7}5|ss(X|;faX}uAxAtZTdxvUZ%}KD4_OF5^KBv-Z`u&Sa%;7ySn0_p`mGg z9Q96=Q0X(EO>jbhwA@z*ISL@0_!)~l{PHp?FE4L9uI;R)LY?Qlw=8G@ZW8228C&C= zS=CsdYORI+o#VrtJ5Gh^FE&;jH<**pkRuq~_$#^mgIV4LvD5u17Ak% z2?xmLOy8qkVftQl$`a>jt>5OphpH|Qn#==CoboH zZquIPR9M=C^K?asQcm5u!mX~{!gE`uXZa8v}vsuqf8Z7+K6yOTZ! zUxY*a^CaZ@1muVz$Qf3*1nArWZ?<-kwtkVN5SqhVK>24dt#x;lkYE@W$7)DzLd+-h z!Qc8X{K4C^=LJWd^94%u|Z|3DB7=Ra8~zpxHf-Z0@lmPYi+;EsPcuGk36$QUDB z!%Q&{^@fa^tKl@70oab&CFSTpnoU0iRaq&&&BOLO3XawDet%&eSygDV#316c(w+f5zN85ZCyaIIFgh<|| z@KCeBp9xEyr^pcb?~70-?gW33z?D)|R=%2%_~_(cTs`a*TiW$N8UZWLAlV*e&sC@% zZv0`Rz8Jy0cD#jd4PGCt2I)II+gzD%EXG(kxKv${rh7;LmK|r z+*R1m_|d{>clk{pH5j=s2h}G}4=wRE$o?gpxn0HjovKY4xrKj=lKvFrmiyZ@%D-v-j<|7qd@b{Sq^Hiy#;_ zwD)Dach|`!fTR%Z=ks}QflJV_0pv&6M>5#+QWc+1OD#&zj~;=~YrUa-nq|+cEmCN;w?=fM{@-?0~XzEeK^aho1Z| z4;JJvPcvvFIVFD>9v;rg2W!-5`)$Y4KOD%2jQ#n8JCaEDpUp8d zGv8paF4%vC5_YzZv)=bNhu{z@CiRwqA{x`BAnFV&z3={!rrw|1%0{Z632anDL!@%Z z!-!syqR^gN1dJE>ZJk+>Kbha9;S$udpj``W^kwykUflVpKWv6Xd$x!mZfEH9=@`S* zDY0>Z;9|kx&qo)4**WwmD=b)Q`j08imzMoU9PzoFOjfHnKCUApZzzcLuvIX@I%&?C z%Ng@0yibyb;w|`UDMg(?Z!uNB#osl*q!HGRZ`;m$o}Y*)`S_lDrlF3Kzkh2`;~r*Q zC$ZsZJ0%*nL4tQzR&)1C;8|JjE>!Y|itABBI`1a)X3 ztR8x1=I2}Ae61it$~8yfAIDXYN~NuhI^C@!v*!0KY-NP}W3fMLU{HRU4+$kD{`?P) z7x6+s^*eP3M!FeI*g99Qgf<^TC(Ig>Cg$7GjkInRtc#Xd-Ff8XnxzAZe~Hwe6YwPV&8JgZ5&e5*FzujN&rm48UJh=5-=bc}?-S<+2ehl1% zf4Osc=Pul+fU<~0pYCHLh<}ZQLSLs?kE^zD5bFCJ`rV&Qy3JVjfn~&#N|b?NE9m34 zH5xmIIKNQ1nTBUb-}k)fkzqzP;|ULX`H@lpANztaeG)WdDSuw zyp)>C6Bw}(dDep=4>)08zs|qE$eIp}k}m0|cGGG zVB@#%!6zn^3M4$`;l{j6lQmw>X8HN~(hS3~IP_2LSVt2fhxduG)g#|>NjsGQ*b7d# z0@RviQMNqnQBiWlOr!3RM;0zc-tTUbGxKX0?W@U#5ZGe$pwFPz7{}v8!-Ao zHScVR!6AkC?@!3H190BWo?J9kQ_iPoR>eB`71u8BKJpZ=1>h(iVFu@and$1))z#on zpXL|yps7!oU3M<7o481MH(UZ={U?`@nRym7wZgsdJ_*dER8}_A#4lL4{4x4uD@xT{ zYN@|aY2J4!H2Bz&GSIs5d99H7LG#iNWTH1pWYjsiUnXMPif=c>iTqh4yLes~pFcyi zf@*0GR+Tc!?LtWB{5xIFfD=r%o7pJSe0&6`kv(c=6}%v`6UZzyN+O8W?>{oM!bPWs z`ozB=Ft69`7c^Jg4j435W1u=`Zo_RZP2C-gNUz4(>H*&HMbio zz_X7OUwGlikO>~(kMh@5>%2AT5?O$bm7gKlO8B0bfB<{b$Di5+`H(IMFgOO zeV4*iHeWrmPm`FxnQ;#j644mvjxc7HcIVTHj??CKie9|GF-%RreauO{uqN;$jQUTp z5fF&C(jyh}zv1U0GwejJ2qU+DM1_L|fwY@;Rhb>&ywUWPM{gi+mU^Ce4jX>-!=2}C z9WItT8}4g7X(Y7Y-T>?rNZxf4%lPcy-5lViPB?68gH;H- z%4BzOC+TW&t!iE|o$>grXJ&WeXw?EUj`)*M>hp%^)@?<cZ zBK^ZH=74xns9kpJuFwj~B~<<(f3O?n!=#0Ex&PR(d_c}XYcrMmo1f~e{Q4w*sYHTb zHLxww=MDSs0rJ}((MiMXmr(i(V}Bx>Bf|c@ab6snh(-^>#Lu?>ga#gKLbT;F7=cV^ zF~3pfl+ob(n$+TTV&49}|9X<4kie~NVc+uuOj@TwwUz8iItvSn)qwMENMiYf%zPEE zIF~14QW{ph-JVO~`j9@GW*#$1aplOAMy+bleWxd3k{MO8=y7*ZG>2luh7`)@D}F-J z8ND)?bTn2^|8yoT+o#M(?BX39J;dr!iyr?krYFkkc8$4+0d))XJ~?Z*4U#}{EkVvuaSQ|5x3_+6R>+T1<*~nWYq}OQFMSzh zj$-`m$bg7uUS`rhA#_~V!eS@^>F;HUWaE5H5$e`99Fk6~ZtMvclLS!V(`c#!aTFm) zRz1}v7-gsXU1*>Q8e3jVD~gJ!{!?!b=PP^+m2PC_*-U1$oeql7nbI&hww0)jVDI(y zl`lGD-0>Fd%+}oqehI8And@gVit|bxaNTla%AHb2EeiGMWykredD%>jsNb8Z2D9CO*tw z%pGQbY>#7?GG#*j@UP&yWSdO4C1Hm*m)|RHYs;Z5Hm?ti9+cC|kzlsc{uEFXvDtbX zrk73mL*n~Cj-SHvbH}FLr)#RNC^)W~2%%mhm1zD%r3Fwlx4M@2<(3q~b`d0yEweH1>r^?!}O({42-LpW*0 z^Q$2CRQ^G7J6KB*@aNF;N&0JUq0gtJGw?*%xElzC@J5PXL?mRsB*6WXGz-ng{*jjmp@%Kv2)sMnLCQA$7^n*FSi3Y6gJz7P|$@#fl>caurIsb%p z`*21S>DK-nX*%%619yJN+oO?=Pn4OemN4*gVT`90NydSnicw<5aPF&%2hm~2cb;An zF5$=n?x*R+V=o4{QvY7YJ?tjIQT(IYmIzCN=soOqJYt&rO*FLmwHT&3F6_z^z{6pCm{vJdD_(zes%N0!V2RPc%0VD!$vcLSDCu zRD^E-K;;nQXl|u6s1m81$KOo<>|f<(J3UHZ5{NYjcJ}v7a)S8gM+fY)Pxwc=I-`{6 z(`6bBegAe*H({1e`l9ejc5$ea6GI#@vCHDHoIKV=Lrr73WtloWN;o_0k-1~V{7LNp z=>?DgPHLeSLJrl{<3?H4WU0Y@jOi)!Y{k;__m-m-l3OcTw>$-&SV-BF^9b@>2b_RN zjvwuR-eHO~>G=KFhS>2bKZSV#zvDWl9Le*~Jk?!vt{A3}@E3Xn+@TU+*4Eat^XV2NiMJ^H?=a5I`mc731uwhwcAUhdBEh+zY^GcnjeIT@dT$gU!KVbq zzMt6wx6p7=5qlmdU9JpzZw;VkZ9nZ1>Orax7rgoTgZCb3qf?Hx#yITW{C)jXUT=l# zHYksxTXZY2#IqnEWjl7qvmk`%{;c!Dq+Rs#k6YNG-D*k6!~CaA+Ha=*NU%dpSynvf zF{A$0CFXd|oJrp^8aTw8#m_mNn{rqZq49 zMT(oL(5amxJ?p_zAWRAwdW?9Mpg!enywFj+Ka4y+7QFkRM#6b4PYj$=%3%vX{#LpA zS{I1OKuhtwGy+yin-2W1UY7F6fPEzdlhBFgD$^woNp5j>Qwj&i)hmXBlsO)$q)^4u z`QX!l3mg&Z4`Qww2wf!NOZK>P`e?4f&p9N79c!%6klJ4BK1F8TN?aF(-^KUQ zIlG`nvwkW2R8bI2T!tI+5FH|0w`gH$8Da7=Y_ePTc2svK_ zHZ!0V8E^Bicah4kCQPlsV*^vQx)SXSIeUqE-e@lFFL_D{;6+)S!)H|+#($LHD^Lxx*cuhFz{opru?d@8ErPGQgK)acf z^}z>PcFYgFlIyvMWK*qLcxW$p=5vJf(DL>1VX{E;Cx;62KOAs(HckE89f<9I-IN4e zO_o?bo;iByO8+jxkUxOXKjN@XnI)-U8}UYaS*F?s<119lDKKft?X#6c9eoFhgUa?k z7eP{8^e8;CZ3iR9)*Ss3B@bqa145-dNU78ZKs37`*kP?ls@e2hYYNlHYg7iR z1t+a;%3{N%wfCthv}O5Usf$-R+j#1usmSoWBQi~LFXT|`1La0m|5cQ zupsYCP{Ce>`4%~DD(4Tapu>hJk7Gn~(jrPq^ZHqbsZmhzzf*5NU~QKnDvHv|`%q~{ zt4`v9n*@}bNWQyg0p*H1l?ir?Um&vz}N!|G+4%=`X7-LQf0is(t zYGBhdWE8bpw{S&JCgLE7>+}v-_2^I?=yzQ%O~LtNhS9ZSwdk*zh*((!Nnfl)G&b|! zoF?-)M}#NCH{ifU_4o!aQ*DwCFY+|jqy2}s{B6M*UHW7rG_GP5No*-Y~?f%>|lm|C$=)~*g}%C8AKp4?%NDESxmvj&hXyO=K!b(>r*)f`ri)WUKNJD3qfe1k!# z%PLTkrBf^ZhRP75#oT0rASV@dNIlVf)h3$JAgk~Y! zkij{UB(kmUg%FR=Gtx0Q-!OpZ(Q(xcAjrZ4$ zeVl>^1Sp|qW@Wv~>3z)c&2QcR-|nNu>hKm_Bm4tp31oM=Fp6*2rMQf{d_SI+=J}7H z{@1Tqe0afDkw05t1n>#z zyk`(0lJaSP&Nq$xF#}$J)>a1{7Kz2l${t^Vi7~Jw7eYzcoBDK&0RM(*md}f`KM_Kq znNi>;v*P+0dY-7t@OdYK7uFPxP+(Oy6&4-A+6T2tj3ERyRMaq8wjYrx8Dzz}AtLhF zxH@|=;O5LT44R0ccv{*lyB%!n93(uSd$XHpA|BKkHTlb%c@>ZQHJ^JFf2U)uvs6e- z?!CG6npiSrOFUG{G|B=C-xPhkEczI9kbRDa8YwglJLP7J;>(tt!%jL<^c3^=^Pu8} zlXc@JsogegVv$YxExeDrjG-8ghd!ejGn<(bDV6UY?a7V|7gkPc>B65xPAX9bORoq^ zKk~zs@0y4CD@SUQsixl}{l~8OA!sNZF-apcv3x;;W6Q9?7XxbDO(vO6;WKBxM7;?mWThJ+t}L%a2<^c#p53zQK5VNm=+Jm~G(HvdG)aY-?s$w2z8z z^VWMeN2~|%4T!=-^VBla$q66Yyg~)H-7aN18Md109V<1-J?F))AZ70ojL>Ur{?H## zFB2uJQXq~W4@^K32Nh=#fhVen`pe5Ev34}RXZ&^kYBa-NZ((=rj-OA1k;aKCD;amh ziBcSc94mS=QXvs12`SaW#$gN(lU4IpBcZqaI}CU%9VRAN)4=_b4wS>13j(-R9d2;SQU%uo5 zx1R8&el{YZd34a@?^J~WNg{3(6;u~-ZVkwLz4y+P@D}fjW{AVaWMPhYD3t z&U_I>G5)h{CsBYm_=Oh+ElGHE&Ljx_K9M2SIjP~kD{I>_vS@KbtR0u?sKg&m(n^=j z`m6nlk#_Yqk`dS_FM#P#1#ilOu=LnGNwfZpSDDxsxE5{dhlO7!+VFSra|v$TnpXN= zGNat86ZI;#9rjK(eu_G8q7=ZWMo>z+iEj&)LF~{OiDu(NXQ)nYpN9Y{a!!o@TEkTjqWX|Bhv)JG z%n}bwym6C)ZNSb=O2lG7(RXY?WbaNkW>H(LR zic4!Hy|iae=!`)moz4MHbId0L=XGL-mG6b|t_q;ul0m1RcDVWHJDI7a>M}iBdDX^& zpUt*tXbLoGW7ILSr1>@a8E)h@z@^6EaYpm`rsl~d*1|CI*rum{@yd?~9%2{_=IVVd zv}(E3k#geVnFs4zyvi55yq_jHv{w)Y%OTfTZ4&qf$ z1>(j_MF3uVcdXC)#SjB?ppBML7gbbI*?AWJubb#*MgoHBr{^FORUKa_Tx%6QR-b`0 z#%0SUQJL=)-Hd5SJ|}7PL!wT~mR~pcgW`uSW>fmb&uP`_Q(uURuyupK4Y|5H24T0? zDiBHyjhvS-;5I)Jz(fMD_nP||&4|Om=zyrTx(*NVl&GS1y?@^33CQsTF9M?C!*~%5 z;>bA2DA+7IBXV(8*}x@cZ!bI?S|sk|7`Q}PIWUMHl}(WyU5D=~2^eHL+yOzuwt8UgHN)d>HS0lrh5CKJ z>e^ak54HX5y{EdGe_7qyOq&MFAG-2fOZN*x^4;4Naita`Xc)?#)<~fJyA-bnGC@0I z!5=@av0`$;^5iJx`~4$=0$4k=g@oI1s~hVt4j5)I>nvvV}%bgxxSeydg6P3f~JCBo7$cIoz9(qw*#q- z-)&cZ?mTDEBsRF-rLG_Zl$elUv}pdV^Mlaz&dcNVe~|!MvAUQ|3QmWDyzo^@(uG6s zodYD8W%}Cpq+!H&<9YpmiiXk72PoIo&I8za);c8uo7Ox2E{aj;WrAN=6#P!?Tg zM(G}DnH?tkwo$(-gW8P5l2z#(no))mmNeEJ*z+qj#vE`~t5H37UpeegF$gaidnjA< z1zA3WQVvfHH$1MrP5*;k@)#VVKtCE{l%^zuZ_%i~T+ytDp1r{r9USeBFhd)==Z;jP z+tRaVtxdnM$BaoyEu0#wgZN+Z2<4j9o)|~NnUt!fWUoR+Jodl!LlkoB9`5~x!YDVhZ1UR^b6vcQ2ZI#=&k103oVCO`FFlJxcueiT`Kzpt5qvia{QzO6VwZs6I|P&2=4tKTbr5&Ul~ zeF~)NrTGxo-l{e0NhB$FW({#<=Y}swN2Ex5R%ZW#37k}9#mX`>BFm_;%j#);)g&q7 zC2!|72dr`}iN-lmBq$;(oqipBOa;qOsFk2Ry;sRopq0+5L*mC z#6ug=)QvJPvQlNDYfV&Flx_n7jGM1R-dqPz`4;Mv&+qHM1Ku8| zL-KvPcF3Upi!Ht25T^emNqe{+JV1Hmk3?|~SVo$Tjn>@5p08NlHP)|vdu9jRs%2DP zjMMTxK}bf_^JSTmrfs@vvt`Ne_)sg8+2T62l@~MD_ATRzs36@I{Vab;BepCbt1edu>1R!6Kg9;jC;oVl z<3`ax!fDwvKhH6N*BU@{W>;rTCuVqhY*kVIecPwg0* z0$DsS7NPj;9gt=?xxf%x@sD~21W7{z<+_S4v8*0L_e?hNkA$|t%xU>cMayoui;@co ziVH#+uKJjh)-<2k{%&qXp3^QjL|NeP%?s z?i}BrTREW8TM^x6H!r)V{$htg*A|0P->dM3;m4_#EKE`~cpO_lVJ=Fp=9)!DE1WR2 zy@`7EYJz{OKWlT#+?vUh&V3Om8T6NdtjVV(S_7vDEB$o`ujmzZiXoMNwny2e}4U_Ss6|JlLzW#tc{X zWx&RQuxY~xY!r(nIgFL`Zh^IV+frfDI_5-iApFcy!H{9o>g*W|vt}l@iS)>f-@07Z z%?@DgSl`FEz824sAwvxgeYs30LQ;}MuLXl7JwvUoKs`0K!IPi6iI?f`4^F?M*YnA3 zX?`W(r~f_rn!i7L>it}@!h<@mmxGrJ{9R77>B4nq@V4}Ow#$2nlq+RcqbI7(NTvUQ4dnJb}I*SCdnHpH5gDP~W z0iD&yA)=vh7gf6>)nBY(TP`ZjpIvd-sNQK#x0PSEB2FF}r`-_c<*ai$G)iajcr%=F zn5-m{!v7SRN(%Zco0S&p+DoKY3dOUCb*Il=(hinI5Y!zQVnQ^TThQt3*%h49X63^u z7=+|#V*T#A`WsTx5Anh;^~c<0X0CvPd+EGtf5-|b3U~->7$QQYKa@wQv$CNoC(;LA zlAxUPcZcW!I{UTF&Kw3k|Jopa5DVSm&G z;|RXLT(>eqEiN>)aoKedaB-2&t>j5;%K!b0ajZV28~-!djr!wWx-0d^TH6?ww^!@35IKy@UMwtv&xs^{jyZM z^aL1DWByxC)>2?&1zlVF1f4R5@yT?~$geHCLy}{Ah`&|!KEp8JrhScRWb3yj>$zz~ zBp1}*Ue6J+awcw7kM35>oQh}bk*(1lbjy9q?P%7+aG-apeqBLhx`5JPYcy99GdH)C z%F4=|;2O4|Lxul#hEStR{3jw`JnSTG{QOTS%!j7)#}kB>uYH3urkB`4b{fE;Po!~j zF&Y$+3x1Lf#mjZ7voxWoc$NV~aav^M#(9*Sq|_B*(Mpc^o`?*dqW-y|PDD3~V*MGC zW|G5kYGwo}zUdLv2{@z>Dm8XDbJZq11_m)l)t`6_$&9gz{JwUG22T;`t`q5$tm;#; z_1b`@pjX{R^_PpV1ShlsT~#-(BJNIggD8uBU&Q?GP4iJ)20Stou$GFRBT17#_M)a| z0E1nzs_cuZne4Em1P5Pw{hy5VlT43!z^|~%t%ygh@L2ybcZ2#lEK)Mmc8&sdrha#g z%@we*j$tlAm$P_MjftVM98IQP*v?LU(~pr!A)m^jID)IsKDxMT5Dlzz8bM3|<=wXoJ-yp@?1=mKW z)_5tfiMFnl}4<^`*C8 z*}B`d68k}!sFa9N!Q6y6PWIV?F+vqEVyF+UM+Xt$*5q#aBQz92p!UqQ$yITUS ztfhV~;eOFp|EEJsFKzKp8{4?;J7ETAq7By3bq_VLJ83QIMANEn*^6g`rZh#JI??yo}Bi- zc7-Q!ZbkB7@6<#Mnv{$)QBbJ7Izkv|u4KOiB%Njgw>+a1TD=7`g6jz<)sJ-pXCy}_ zp2{c}GyZ6uqRAkwfF->uMT;Uwnv%`sNHgrsqRo%A5LUbw?gY~Bfx&A$B5EriK=XQD z{>Y}QxS-Ji%iQW!P)C76=g@(QNz=yrVNDdQeW!~NH^GHStFsju)*wGHuA}Y`s|a1+ z79d`ck?SX>FdPf4>7LYW{kY3OGfd))WNxHv{m?c*VKfN5xl{wHq_wcUQo1nF>S(>Amr ze9B=ZUeoXUMpSufjO&rJV5_nwDYc@8<2dzimSUdURN247Sq}SfgOs`ZI{A_+NSYDlSf>8HiZORZ;1w zGWAb<++O04EI+VVxnI4y2_Ye)#}r>bEUYHh#~x%>X>lOcWj=OIfQuuM9r~>kla_u+ z#nF+bOjnujT76(vO^glD)Dai<_2*m)DTR`*6K$@a8ao%92ZmgC`u57zyUqEexx?gC z$MlB&G1-0BsZ&5e%XH-rAf9kkL`>y;!2>Z$(f7z?eK=NXx1@+R@le%d_dyIe_t>^m ze2g_X1aW(=BTN&CXq5=h?-8ZHJ}r%0-{s#2?>B3kFhjUX)io%;GWr6?`S~I1ck96g z_j>R9BCdQHC*B&yyx1VbSJqq%+U|^D-F=bOn>CMs{7qAMa+AI{BoODc$c$qh3ZoOz zm#y^-DkX(l={Q#(26&&oW1CUFZ-wNN8}>W(ntUA z{{jRDj<^L|j9up}tG+pqP`{2Kya*Ye-zBh)r(R!sKsPL31QkG)N15_qi>#W8$_$UW zDSQQ@wSK%Fb?3&A`&FKHx7PQLWt8=|XYUm)c~Qh3cMMu|T)VU4pt*rfIVr1Z1ID9N zaT!qK&<|p`etRXb;lr0tVnPwXM$hA@zQ$78recWQB8LLYuRf1&n2e+H{rV8D32Vzy z`C(JWSaAgwuwfEsq;@VQOJ0;!hI}qoy#>4CzB7t1`_Co^uo6OdMj4TVF(0e0gMjt0 zspIm2E-QU%xYwks{v}%yLclTRcsF7HddDxq@mpC71zr-=zbjMP9x1wai=6v9QOj^= zP6|jC5~L{k>kWeZ^QC%y-A=z3ROg-ekw%uS`fcsoNr~PdiyZmO_m8vZAN5XOj!6qP z&QTs-0p_daQ%#Ul9w*0E29mnlRo{}*X8|tBC;X_n)=GoLW}j`?>uybo2c#F}|3%YP z2DSA!-IfBSSaBud(L#Z_e4$`ds;^@MKUdth|_J0_f4H+rOA0|K=3cMdU#a%<@h}&ZsSEERpo%Jig?rR zzLBz~Hiz!m%I~&6D8#9F2HJusb6g^wc6?dowF|&x?dn-hUD*(s+GI z5TCElKi1IS`16KIHL+S)9TMc9P2(J_Xe#x*WdlCGS^Keeh+@E+nwg!=NrJiW-bPMA zF=^e<(T9JhgXv|@gxUM0<;sp{w&FQIK+K>G>+i2nD#0-=q17H7hbMAWVXewY-}rKF zMuYJOoac4I>vDG{Le~!}qTVVnv9fMB%ADUJ1@$;ry^k*3xKQi!*0{V| z9xv~nrBf|55~$d3q#CF6N<(0D=tcqf zU{x1iecbTzGse=%%cg&6=j-!Z|DHiT?glkW@9%yOeAd44SQ{6P!TG&3^~ee zmGQ-C0f?#=g;)DuwO<5!zzDFX+GWbdtNz(K=na*kCltW(2&b*4kuE`-K2cUwC80rs znxbPDz72u$Sak312J!A8-MW)enU)*QOjzI@N;Yj7K!KGzu)Lbbke$QzNdCTYku!##mQA;77hP5bevx7r@Ua|r zcZdNE?)Wsq@?IX16}ch;UQ&YEJHFctZMpQSe8Hmvun7(?XlY2jLwx<%v2+N~_qh+; zP-)tLHL!=${(v`j+XDT)%};*HAZP~Xlv6@mlNcK)q&27 zCZd(XE2pvH6;(QUmo0e@9m`fGFcXLD=7@$6gJE>yryN0CJ z-I#~DwCDhBj8TZT z@Olw)lx~=B0>{|HDzD4D(ETL-jaj>v;b=H1D%Kp%ODVMo6yXZHO3az#DR0<#D(z7j zuf(7*(Db$w;cC>y2F28rdmzFU>m+PF$H8kyI<68Bu0s(&$|Qi*E$NOjV8Nf%Vpe7W zD~)lrium=Td3JMql5{2Dc>FV*+aRE@U{3A2b(ExpiNy*vkPRLU_F1A(Gm}zoCP@~o zJ9`5WSH2ij5^rL>PbJ0|L|JJZXkYnZfh~sZ-Z3M%w^l|q`q=%(zoKC2N-*31;b4x& zVW|zRfJRyh9MX9g-MOyt?^>ayxp`gEe_N7IP;mR|=Rd(Q!K>jV;kn?q@9!$s9$17g zA{88+oqG-snf^o29n9*<=Tw$muMD0KHQx0mJU{y_p8WNJ5_PPE;3K~tyc*{Xh)`2& zNB@HSy6@xDIa5q-K|#;5oSHZ+y*TbA<&PU#q9GKtSYrm3TC5TIAl+o&qA!KCJ6wm& z?AGhaSP~xg3%pG@jWm?cd<(RCeT6*dyNB~uv5KEe%c$ z(`a84#q@kmxYm}wJ-L1=>wG8^{7+YVyX^8meT$@|;0&gl`2qqF@@PdgL8y3!fzVT(G4^dI&cAxpA6j3}BiYE}T@ zJ^XHN2P0K118sChG{0^@Fw`Nv==2KU241&9%C%s$eM5U3>x{kG=KBOg}XS>)xcbh%u( zzDg(9N+8%$du;cu(1;HJ8wjdbBq3e8P+=6qod>CvMFeU^k+nT-EM0YE60>t*CJrSPc4k8c zY)xOD7p$+^OA$y#EuY7a6-oiteEDg|WU3jz=m4|%$(@?839Q%b7L~c zyQ^}_`+8&1zdi+W=NA@`YSrTskj*Ya>Jp9HPm~!U>1M?i{2$Q|jRL81n2ZWhQt|hJ zMvkA$_orFZ%rwi?1aI@qDz?a;dtJIW=yYRP)N(l%9fcZV%ppEY;CvYz5+tPy&zuPf z$7Wh-3`Hht1NCxuWcpMJqgM)cM1sOxwD) zP_^!;9adhh&k*1*30n*tQ!uI`p=p|g~r*xy446&=WGMQdkZ3``>1 zq}m@zh7E0I9Qo1DLukuoh$4*~u)xFBb=H~ytf@jH&Oj}saKG>+_Y9m+qI{N$oXvgJ zB*LIxyMcdMtbi{cM=q@T1OsCUBlq}m;l2=jcf#Jte$DQ>@F>9Cf6tb6be`r490sT) z`S2?w1+`3U(Ff+R8THYbhW0>q%r+a^A@R{G;wW#-w>lIZ4)+heF$~Ga*LvS9q3W3Q z(3O!_v|U&0-D;uqG$%!(AThA@I%tjJ-s(Ak0D)_ zE6ikmRrPlwKEp6kn7a8*2oX} zAs;=*u8OAfFX`C`WdTLgjlDR@&zI}swJR6Di)e%GXoDx}fb+Kjug?E>CenyAIjShJ zY*MS5b0GLe^h-d)>bDzN5aEw#_xEbz^0Hj@!(1X?8f8^PzPRhR-DIGC@Pd>ldk(_Z zLOgG%F$mIXsg+S?@3zx(a3Di444AJY)BG`+eAuYpu`hDPJggh>PulfF$Fr>PnI$8C z59~*JAIEYcAec?kiY3LBnq9vd8omc9ZtS)ka#XGvcXUMz8c&}+*`k~^==Og%H zEO>$q`DrY%`D3rmth_Fi;Iq(zxv;AC1jF+6mr4aCD@b(6l}vAcY0)GV=$_oZ1%7T# zLn4Ps8Nl=k@%#qSjm3IalMyN<(cROHAy6{Za?sd@ZMuFcVd~}z7w&i)tkkz4-o5Jn zt(XxmGqOL-jQ08_WjG3Ch)z7bOE=3AqV-uVY)( z>p7&qjlpf-X6H^RSLZVw)GW^=+t)_^wb6KTdimNB_uG1rp@K8- z^Qf}OR*|%avBeMl+YMV9SUBWLIy^K2aMv7{Bmh>kpbKGJ-VM zM@r~93jvQSF|u|0T&-n_h9V?PBSb~?wFR)G^~rqKR&ab#6dbJdeN}5FAq1ft#%Q-q zw7uz6=v@xh`uww|D9cM1Qp&MFC}-9A=kB2aXjk1+@uTowFu2NBrEW11(FD<|UAKLe zRb=SFWym*MUICgMHe?a>Rj#f%oLEUoF5LB1RDwXx+RpN0%>7z{6i8)G1+bx?r1B`E%r~>zbHVCEQy@S;1d2Bm&L4+g&-uDHiEc_QGq#P6jTxls zc`l`OerUdW9C<>F9M63s^ZyNMbo7kG>VX9Y)D zhKv#v>($BPI}LbRibFubH9MR7N7Pp=(?XdL%8gkSJxgF5}H=}qU2rUETQ^`xG_*8uOO==^kDh&!iS=ttm6U{_pQ3_wOR^$I%BP_ zQ}YL~d6+kVzYVak@3mc6$bh*qadJ+_m)>*!MXBoz^yJQB+@=I@#nyh3(T^muRo zbyf4j-yc4TESUnI;YNj_xL!566 z7BkW)hk|9|B}_?K36A7OxN*yyxKmZcs}9Kfsn@_J6B3)!GV1S)DT&f^^te0`#yHTt zAouNREt5={1>+6YVK_nFxJhftgk8^#+xGyJY+dWwFVH$0O&j7&ECS|{Aa0h(f`Zi$ zhmHcvY2+2iW(qywtp`D(!`*Ig2z!EUg%A19lxv@;SZv*RkUT5%VU1ZsLjlv9xY0MW z>I6$4(96vQE7C@CJ9m?RfHWnxz*+u3IEdbDx0v&UoA<)ryA|?(v}$~O-La=X7Cu1Q zUODhLEI!%MdDk%X(T8fJm6{DRIPu#%xa4W5YD%zIsW^OB8*Zq|p`$5kfQJH^sNdWP z$7q~IX-y{3c6Kte0+Pq+?uID~wA^*O|L*%Fxj9wB9MrlA=iV|_yr_9}i<-EP=6g44y=RkGW}YyBXNZcS$_=q>S$y< zhEet~5>w;;wGGqU{Hc~gjzzn~D~`CLStJp#i);w%WeJ@KdhQU1CvfxDAL7237i^qE zPHLe;1H>}WYpJ`Mq{Uzn=l?KYsa>p14$UR~^1i9NWPNk7x>E^9{PkXN8hp9Bvv9wj zU`x7hYbxjj3z5}C(T4K#W6Z^Bnn-g5>A}quYRC4dt~y|lPIdq-RUCl=Wr)-Ug^T;n z`K-H`5<%F(uO5;SL}84YTwZStyHShZnu_XER*6RuPvnGJPs(0i8)dO45QLg>&Xv1j zjrGgD+S(%{zQ2)IY*FqkuGjXY%)R8yl<}^8{Kwj` z#b&JA(ae*`RUpDb64@y_wS5jcw)dZniWnne`f_}HJbrYu>Iy31e%cbdZ6WHa^>ll% z>w_EssA3KUNxerW`vo*})ey-Tu07=lOe$}>meXxBY9+X%@ThVxP$rU*wENNuKEBpv zg}RUis~@M~DE|OMX!iPZF}>vw4u}=zPFwGv_GRm4@DZJ}V+W~hUDazJ@Wwol!2}ox zM$V>U!G`e9?-VE_Yf#NGM&RQ8e5I&0tpC zf>Av9l(DH4g^0@5YWJQNMPbyXG7RPR7O8#|#w1qwCK)_!EUu+}>lf?3XIJ^7U^CwbcEq;9NBx@(^f?_v*PIHErm)d z$bqkr>4`}-hfdXHv)hWa)x#v@F7Bj3?YlJ{Wb5}q+69(P1kVi-_5&7el+`ou-}zn$ z1j6lHI@*g6l^K4O>2rJ9&gQ%KA9F1$E33(_r*|X4;Gxi9)8XG=hfe0Nb3Wzo*^n>v zQyh}uAxTO2=Dg%<3jUW#4nEr#1jiLUz09Zjlgb#>29CZbiG{*_+=lID&CqwV(6H`h zJw+s+V%s(-qqp41jqb#me1Y!!0;WHMw#Ns&eBK8rjs2T5r(t*v&|!LtLj`R$R-+VTI4Op*~RerIJ(G6u{L zio^-?`NI0hBF(hf)uPIxGNjJ-8I8KKKZrYYB)fn?M!Tk89eQL0g0?UQ)X_`fjVw#; zVhp%4fGje($cPc?ogkEUXn4h@z8zPYz}^WBVy<@?o2twMu)WTlJKhJma7cOXR)?Ko zMDs%}>#Z{vgMMoVejY|G3O+pY5C(zy@GP{fr zg^={Zfn-STca1oOW9Il1e<>(G=A5=nv7Oh3NR)V8hjL;Osg}KzO7~(^?_`xEp;P&C zY`*fZ^Fi=xRPGGsMLUrhij(2|+F#N*W^Fm~_99DOKh{bYaKW}=u{_v%5D62q;{W@? zY7lPK@BS6}^(?4TlD0kNm5HH1Rm5lM><_5DL7mT&Y3frJ4M$395bbi!6YqMc@JZ2t z`d3q%&s~Qc78djf3N|)2_NnJz2Ogn+I;<#b2Y9E!mcRXVD5Fp@&vauF_xKUT9r4sX zWyy^w@EeaEXO9Tg-CelHF>UKiXf*6n>iEap5`zl*9qLg z`&?c#t=2b3*<(KE%;S6^n zFohROGhd00`lYrXqw?JyjP| zu$AR!Z##*JVaGJ4BnqWZsXW965-B0=VOupo5g(`G&Lh^Nm;;QVqfMzHKkcW1Nf5q* zG-Ququ+ggu>>2XrO~&=}q^m7|iI9ds2c^icq}z6Ex#Ja7X4MIe45gDYt&|W)6`O4; zpdqC+VK1`INEd0uGsr&=tCMBAeCnrmrCl$t7vux+%|~Ws!5U>P2v-_28YQO&jV1$g zj_xxJc0U2WWN>Kh(&<)UAqoY>$6`yn{i<++!i|ljm>p-4vI7nZ3r(F#TTk2as>zba z`s|-IQ}LR57YnD70ZOv)OH>ZUeaq?aZ z1Z`bK#j1sNeXn8AaWY%L{m>dU zp?}>&-PcF1xVX5x|=T;X3hznjWvr1^w|b}lc$2vYCqTIau_zaEI*3*fDX zI@T9vXi{QpuS;yh5A#gVdgES8#sy2OK~PS^?nUT2NI>P%ih4oIEhwr1Vpjj-08jlx zm^j_S;kIt9uGmKfrBARUl*9C4k>I6<)Yaka5lFDNsx@0XxI1=F1(qQUZ%u!TA9UEb7 zozoK(f|eFsM@uAyaqHFZlvpgq#3)OCl$87Rf}>UFJacdI6l)i6NMg|GTtqrG-T?X` z2oO$wDuLWl$NVULEf=tT53mpeV$)1d{icl_3f}9rjo8i+SgH&R`Zm zPZvcip;Nm&ZcjF=PP=6lsw|&UlDz9(fQXqFs1c8fY2bH8IB0I`5_G=LVD_)SyP^~@ zWhL&=S6mJVR+28;z(vLzOVZUNCNa@GyO%A-z=5)(H3w*JZES_qUBx`FkeQh=0AnYvOruP=aRczvbsz&I<9&BFxr5H$LK?Kw8`S zJ3gL1NzC@!&>uBHv)=}H2gr7vd8n!Fl*nkWO$@whY-(DVpZC;>8cgkcGV@=5>%Xqo zaz@YpbffVG9gEp*$dVD}g8@+s))!MzcXT$h`W|;s-M;yqaAaOtErDqKtd;M!qA~+% zPiDZg8H)CImSsR=ygabVkem3}fut#Og7)e+dzyDSFFtBr#CTQ2bCeZXZ5>#Jw9gYe zD{OlsvS=+y>RoTm)p144jXH-}X+cYf||v zE#sP}dDb6exVkW+($;PvY-B3U$U#-y+B zI@UK!4bv^%o3ga4%ek$zb7Q;7JgiIRkexo^fH<-ceXzexTN!6NfG?rbg2M6~sMsE# z3!HhU%EyrdO=}#ZJKdK$aUspOxR}gY34$J9Q^81U_=0@$7X;+T)8trVeEuHoMSPvZ z@5K5*bDz-&8Wbn(@5^s|dHHX}+<87Pviz%Lw5cG)4kxFioU(7?=i|GB9G%8CU!6UT zo}I%K{0?=fcRB|S<>=eq%e~2#FBbH^jzCOlA07n0eED)VNc9il-MaPi{OCX?f2Z!| z%5wim>?gXm$s9e3KZq_r!0{ZnV~mu0R=w!ObXv43zmh;4WL}OdaW%$>56xC+xB zGAp#J(3Wu7Zvb>m|9bjeJ?Tt=Esgrtue;Y}$ouVKt_?sfL}A}-zZDoRv@9E8lSXSq8R)m{-%%DLUL;jKHVbMemiR`9A&C+uc!C6twBc#kwF zh#PO@t_P}QBYs4Y4KU2>5iOo29HB}kh%C?YfE)x97-`cRs)6qzvj zc_>RNUq98pX7=22I`iUb=cX){X9ROKHZcCI+fkJxI|Y%oun#8^c1H;q(P-*giWXF- zsY73}hNDD__qDgDE#&ni5t(pVubIF5J;4iiB^I+ewa&5V2> zn`0VZqI$C{W@%s>u7>>c$uY$CV1H{##T3ZFvSRcG+$$0cQO~R!=FWAN?6?r@qSiVh zi^3zTLCU<@N{Ws;aigc{(yk1n&4Ng3k=Dvi*{!f{FroH>AAN9AhX~{&ySTAsa7e0n|!{z<>Q#!I1>#jY}CZl%X>MEKsUS< z;d*Tlc-wEfklez(`bW*bx6?_3GtI#Hv8IJOCfu|2&?pH#s?hC20Sj&x&%+#f2x;?jViN|f$C%nlo(a0 z*s7=WrV>X!emSaWE~ZgPaN|k4N^4dR7R2+0{9w``m{0Hj$st5wC2x++O12rc$5w0C4aq za>q~xY!20Cc-6bmg3treaBuy3=>qR?7d_9QH<6nQ(|WTsc~ktla1f0&~mHQjmy9A z>*^E#&dpuIj&O?EL5Be|GEt9^v{IyLN}=ze(@E8BE=E=g7N$HvN9{?v8{Cxvj13?BFPs z^r&H6VDNTsVc!KFg9_R2P)eXxK5FAo?jp0mt)`eUbPMtNB0k z$)oyZC4DrK7*Kx_7y*(Qon)HCLpj3ahgNi-KJ-hFGg18H3+!10{pd@-j_Idm{urt} zx+bz7H-EBYY34RsEhYg!Us+1hH*feD8zeQf;+wUT=|NEMDrLk_Jezq_gKt!waoUova zkAgQ=_rr)CK&k-Du;Cr5`#?Ucf^L^tZP89VZqM)K1ttTR8yDi|6s#A&Trqw98%bAO zJfLe}fUvahot>RsI<<)5cp6!I!s`6@*g4kUkBaaa>4d+3S%2@ZNOorD*Zkq!qcp6o zx!JHq*ZJSe^X=19FaMK+!DY2tltrrU*azh7LQZhKBCiZ^Pn2lw_`zDF51zu+y|9EM zDA0t6jR+&|@D-p8*ixULgWMJrrgA&q%aYo0oC#Q|&mgp4{HTDAq^1brFT^I5CGhc! z@ABx*G_$%8YHqI7q=Q&s`yxZG%J8BM8}c;)%+`q+r$V$OavxuO90op`u)U?7+lqQd zv^|+?bzOYzt>&;|wL=NmGn?mgW3wgKzi4iCD}fyP;NUs{VgQk_8|CAJo6QDe8d`7t z&Zv3TLaV7owmSq)k^cr3uugte;HbNNuUt3jl;O{q2^zK>T9yB<)jFR`pTpk7aF zm^efm8KGIi==RO!pO0k(&6S7ZcxD?4GLppUQj-2g7O@HxYD`Yg13co@}QBdvzj3Unb{WQm3OX7!e+2^HW{L5u)7L-Iw{f`9xBmV5?b!v z*w@u<$*)Ql+++xf;QZ23{@Ph5i)u;3F5ACLw)-C9j;=by#;2-(4^>wP@$0pa-p{o+ zH*7-w7W~`ra=l2}x!3cZq*EB=2cn^>( zB}Drqt}OW7artY$MATVlXBZ=p+#VFiiiseiHrVlhP#xYyFj)|>Q&bLkJR*!1AG~KU zt)s8@6l`Je`a3QOB5unE9$))GB7>P+%@WaBQR7Z;hK*j$c4KJH>)tB%$b}NWY2wCj z4pqzkkc>Q*L<1{Q!l!}JqN8Z(8!=(^M84`W>JPNgGzb6T^$e6G@T5ua+9Wa+ZN+*V zipLj_i;vAd`xNfgH{*t26}WN{x+F!|S^@>Z!#H7hRG%GyAu${mLqj+{iu z-$C6ZGw{Btky_Vm?sFHeTsv5rh8-F8SmtS+G4`;xrYEp8rOy`=lOo1d?nAxW-LC2! zIs~>3&0$1$NL1HZxD$v^8lcc(&7#*QlH?G1%My=Ud)H?O2Q9&c_g4LbY1+Dd3Qa-~ z46M1{p!NzhDrPW^q)njm>zmdwL(NJ1(l&`N(Bzn6>{4mG0oU{7=H6VkxBP877bDcz zu4b^9B(I;Ddx!x{yZsLy$bGGG2ZNRYR{Aiv5>4diNU*&EM>DJD8A?d5sR_%vWIzdg z!@e8A8uiM?58L?o48ZJdTOZM`kRTuUt`w4dTK_+D{-+Jf7pkv+O}_S5Jug<(bawvB zt~xnOd~oD7bi@e1)kbhq5%#th&j|M-hsXK3(g1sd&bdirRF%7*|Jb?({Zc@U4d@-? zzdXwixDmYSbc?7P^&U&yO<0t9wP|g{ALc%ZxpUJT%ilpaagM5GgLxqv-LZ@hv(RmC zWn3SA@kA~KQKIxs$E7zWN`v462wzB2S=6s;-)r`7n!Q)=T@GuE?Zo-5D$c%>Hq5j9 z$TXH&x%;e|sR!{r_771$b&8kgY-Q>D8)h?)7?U98Gb(SXvRw&E5X}pyTzm!aCt!#a z(}m6jzEN4{31A+vJ$(wuxQPN`6c8k*aIqT{1!H_=-)?haG$U^x^4Lten69lORgCs4C)_T~1z;+xV}_wrgc-XPhACeHV|nP> zbF;h#CtmLjMO{T_o{U6zkBEJUu3R1?k_6LR@UI}q0RwHga7#qZjMe``J%QXo^Q+s3 z81d<{pw%davJIsJojsroi2&bWzysyuu^kzOkaq}2Sa9$g1oSwP-MbR(lmQ}KV`!xB z5OhoC|E#0acAVj-r>tLZUB3ROieGpyb`=`%9;=xl>h^f^ry1#gD)75)1jW*h@~Yp8 zSbxpaK&KA*<@gMxUwcs0g~xT&?Cpp1{bjrd`EGTqF6*csXJ0cdD=k4-3)J)a3uEJ8 zGf8+zM=nuM-+b{@wd)aKomaQ#M$HU<3{h0yf(Mm*|k^SPy_T=K2 zm8M_4)}TQ{W8e~+;#v`EXPeRU5qk9d*F4k5eT2z1jPxv2Jmxy1>bBzFdU(6Mcl*d% zAvMVK~GH- z3<6T8PhYg-=;1eKFXCnt8Gq139cJa5H~$7O$BAD85Q^pDXxDm+ouGV?EH#}(J!-Y< zk-8@t+w&99Mkaf0RPz0`NJ*l+HDBeYbaoCotflGnG^YeS@3!Mzc6%qdBtAHf2ZJpF z2hH~@O*(s8nv)bQrgkl==(pTsPTsRbIT3_)@BU%4XW-$Q{AlwwCtXW$iW5oqbh(&Cn=^HiG6xEX@MpGshB3eDlc`G1QP57_+VX-<(d_yJs;7R0krdGt z)kPW&?)yj@RVGCj*9ULer*G;17`%y9OR1(@u>Dpzz|z4Cs-n1K8f{N-t#IA`0$^E8 zUXA2vR&*e)E+kk)4FJ0t`>Pc4b90yl zMh_vD%7htJ8})|K69@ypG8f|Uu6ll6YI4}SsR|Gx)^o*07=auw)P=i@@`sOP@KSCU z->3c7e(g{FUFo-;$I6VXpk8;uVMROulU|Z z^kQ!odN&>+qYzI0q*3-dvN>*#`Q_j^NP(Kwa_I1TqN%wOmbUz=?dKn&9Hwjak_yGr zl4X^fvHc&PaCoYkCV(Z=q7t|K^$jG~Zd%RRW5b!f&bT}kywxT|N+swTI1;WN^KvpH z)-0b=xXJoCT3eS<59R+nM+|JnD74l$Y1cO~aTS2nXiZ!jm!=P? zieL`~(4v^=B6D)@R%6r2D8k9*k9&a2Az^s`spb^QJ#A{}ufvZ=Pyr4}TrAXJL+{Z- zTAPF+CZ8cu#vIklgF|*L)@fX7>cr}5mSAU~q5EoaySSkiXoEx|iW3rZ1lXosh>_Nb z4lS}GkSSeTnlWPjXC%h`oAa=MpW|d4%9Y`}U()?7Plx$8P6J*fC$Pj)ph^cyDSgO5 z$=&{!nYpn{+=$yWh}fl;Buht1POfXV@!nI17L(eTBN&E~rK>$Ymy7vFvfuAPAjNLeS$Y!e{H`YgnqQ~N zdGIH+-jpELdd^eIp;&~ckTuxm25m02-qh9l=0{qEQblQ&0f97A!w=j3E~{%z=9z>v z7Bo>cKAb!ezu>~4@LV6X+(eKjfN)0jYtp2eb%bKuN*qr{B~Ks9D4WmlFzlVbAYBB_ zwDll&asCWK|Mx|r;CC&s$U=`>?p}1iRe}UWhj@)Uv0}8&gYHs6gcmKIJQ$=Yn-FzG zQWUm_QDiQp92xOf9iz8e%_f%uEd2q%dQY>(9Lz>*mYwc#aMTz%ki6_`#!eFdVR*D{ zWCKE{ItJ`Ko?0XRyy)P888yPeuUFhYVo!FE%xoPovtX`Liu~L5;}8)J3SjqW8)JOQ zEKz!CBf7)m*LpHnk2JK`_d7q*y?(k0(Y4ewZGJX!6Exr)XO0*!_7ol- z_W5layUMyFz3z=0nn$+DUuj#qR`6;p36V%!VkS)Kk4dZWk-;nAgx;@dpjU;O#i}Zs%^b8&?6-uK zA4LLRcC8O7rIW|b>&0xMrEN})%$A)e@~h})m(>ni@oc9*!nb9!vavI4XDv_Tn1Mh4 zX2Om0Ju+h$fqJyIZ8O%L#C+@Yzgpow9`lOBvmPfkq?aMOs|q0(`KPgl9*W6Us%J05-G;H=Yy`o(Vx0 zFaanlscx}qAw5>1XP?YsiTZ9j6q;98m-O=T@&Rc@kSmIwWMB#BltPn=cwfqMtklxgZzo2RA z_tj_QL-3p;6m6GSNMOm;_bW0ivsqVR5_*9!CKVTuq20IVWye&@v2l|w+?`r6QYrI7 z?vR;DGzHToayE(hJtDPple)r7TIX>aLxEra5Ekpi)J~qtb>Sa&Me;d6B6c z52-GiC2Gn;(3pj>l#qKrt7+FV{o+JZ)5%r5w6aqgwP~$vZpjg({2Md#r-)j+7HKI? z#U)qs`ruTR->tO>Gr1s!%#(DDhj9h@D|y=D=~rn?v3Z*4Y1u9ZXVHl< zWjrB!;(Y#6bXJ^AC8%Zx$(t{%gsNwnuxkqcGX(z<>91T{y(!SU6^xdv1Pe`;p*T6= z7)*J!q}s`H|(4Z9Nq5_NKdXQbyKViDF9n%z2eZJb(Ty-l)GSw z?;7Dn06-*~D6q>_uRHo5cW?CQpi&uY(%9a?!GY0AtD6EejjGUeZ!}BSZAL97Hr5Ft zaCPZMAtq0~5DAa4@P92#0$hHvmwjwv@m(v{uHd1spRp#KQ~Lt;6qKh<{De5|11^Of zOo|Zj&$XkAskLTRYw*Z3S>PT!Kd~_xSa;FTSwBJUo%bR)C?%mi>GjHqz(TTen^uwa z8?|a(;OK4c-k=(=li)(GlSnrGThH=fG@UKrsqe8dK`OsS5V;oDe;XN@( z6mD?404zizmmCoQ?DIC%A&27Dv4t3~lU#&_>`)skN<`|SU+;kaxP)!~k*Zb2OwYX7 zKV;vxR>9WWnqn;k{efp}St6C*M8LF&mnKPcEjAvhqGhl#WN`hmkfyl5F)>zcNzj<& z4i?n`X@b)1s+$3k3j!=!VjTFQB&#a4sT?Hx(YVrPkd8EsWhJU<_;)QP36ik>>2evJ zXlb7O{MRng!gW`y8E^4P!Ym&|??fU#PQmy5Y?;5!`QkZc?Bs#8S6dk;s(DqJsy``H z2kzu4u-YZaGQ3!7x6c=>$!_bH&M0^P0n+v z=Qs4_)%l5=5h(-iw`=LL=_j}Y5ve^a(bRX6%G~Dqw%<|6tbfagsm04SPkwmzz?OF+ z1RYt8{r|ZDU05_b&OefI50c56+xLgn<@w>pzP`U0v-b}Vt^MVD*nJ6N%hk)A&f0G; z|D*UJ2Xb0_Z+%$%gD{=nEe1Gxqi&=d{Ge2bd(v!M=_TFVd-T7;YF5rUx)D$wq}%ZQ zhor}pk_Lq$>|!1s9xTuHw>z?4-=4TG01pu(?NRltg}zrqc=cv#U)^m7dvA?YIu{FB z0Vs&5K=qc32|IY)1jpOi&}?rw2sB>RQ&7qTziZ!6EL39fDo6q~+Gri7@ySQn z!2-c{4^pIw5eQT!f*&-SIl&$5kc_1(( z{_`apO|PlO)-0x@N4Z`&FA+c)qIxBx6!L?Jcy43>;P7LNdr z{dCqeo2VaSMUR?s3qe0v5wJ5^yrja2#PCo6icQ?v8XhQ2^6_Kel_2Wdw^lYbUvjVp zE1FgZXEY4t5UjrcfLW0hCnskQtT=n{9zpZe6t)@Q$i7FpOSEq8g`0Yo9%)f!rh;e@ z+o@ryOo*&tw^X1ElO)_*__e*~F}J|Gg?P%e6gAWeiX@G1B+7FAZj~#C*AJp7#YDzu5k-|R4^33rnvx#ZYEBRn>Chj&eMr82WXKHV}pj(Pe(pDlkj4X zC5PFg1IM2}NH6yeyMfBwkYpKaxcJEOJx zN`cp!|E8VF-Uk1zr}U>ne8A+8qrUkwNq#>}&afts$8BPz`HA*3Zh7RXpA04|vmCCW%R$eh6fO1m-Ro!k zEn@sV++W_lvmR|Jnq*+q>y?-Otf2FEQJ;YzTxoK!mb0Kvhbt^qw5}o5tD}z_O`Vh= zvIJ?=P2-UYcGofer)NpRaiGC*2b+pC%S33FZ|SylwKexkqSpf3u39iEjbjT_2S3w;gQSvpzOG)nURX~BqxGPtQ zSCF;BDo;j#@JYn}4*zSWMCZ}UpP~V@BF37O$iqtQ%OQzwo;qF~yDw*UKYfwv+>VRb z8xDe|cscYPEt1=dBq>q?h{0*nygik!9U74I(9o0i6w*=s+A>20Ka@yOGAy3=-%kNZ zFE_W+hsqt&mVpBvu>K<#4@hgBcRHXa-_4cYF5Wa0&M`@2di<@7SEj%{8^rW`dZAl! z%}9t|?-TLjvm)s7y)f3nnnH_7HkX=r0JD-_#<2M++_?Qj?#E9XV|}#}5e^dZTQyUE zBMtL~es^4h5;V&#C~>m5;N!j;CHLu-n>dZ~ez&5nq?5|R+_n)C@#JKGz1m=fQZ(ET zE(}xd(IvV`lnYA~4DotXC%j4gIy!DV%1N@u>7`+AJ)1;<$0v_>1(L+$REqzfKco1q>wfM1DCLqG9*hVu1?R#p=dt3M_wb+1myc!G9_Q63>Rn_4a z`j~874aMDa(OzHn+W`5b9g!Ri59?G|fdt9>G!!^IHTA=|2Z#o7ltYeR#{zz*b1?6j zEj~J3>)@Vnj9E3_6K4vkAF29BaT>f$l0`lIGOu%UDP0%78n<6@-*kX;zPnDdky?We2odJDiUh;GVIV^9Yg*s`a z$eONWOVA^k*fo-_-GF$ut_u%EaU;ltl`$m_5>UhD0?S>yAt&UMlrRu9*nc}Uk_)}C zvaiwMDh+U;uk+B^lamvGgwt86!0pdm$#?)Z9z4UKsZOxGL-dF;imQWW3(@d=mhY~(Fmz*4N-H( zc?1kmp^xA1IUuoaSKXd`GJie#R_4{n50-b2POl^G4I|x#r3LjMl@#$xsY*)lBunW? zEQR5u&$sbQ3awr`XnFDe48WubV_&1s^#6cFcEx@bM7L0{<0YD+nDrcfU7>R(;`}CF z&%im0p~NIc@u_*;)Rc9?7z=@XA+Z`OMqQ&m2Pei284{~hOqxiz(%hH!=^jU18dRyC zx-S~%-a9!1q;t!fTu+@1NJeVEZ@7$+^$o*Kon2#vnvZ}z`X4l{J+5rOn9FZZlwY0I z4a&g;x_~@9-A6~IfvKu%9q+4?DJ19R(PYlTtX7gc)cbBL%vEP06ORR1_>q>LB+QlB zBr@Z~n7~1T5=@&{Hf3`VwD8Je?kqFh-n!0A>(&?lRrFV{?iuCRw=c}|o{kLrN+ZQ- z)x1FYSZRA}Y?-Dt`M&p`SFfEUZRa|@`3olV58M=^kEV0()+(`+OKP0`an(43Tm-ac zO|e334EC#|VZ3&3c`7B8ELweb}QG&*An-BeyZ0&ZR zek%$Ve(S~C9jkBKJGIV~vC=GN^NIfIV z6H{ixn=+DBU9pBpt>V#{M?XvAFXwL@{VN8eH>2ZXXP5m+sq_Wbq=`|0%%1*kt}H5+ zC}uD^BYfCTC4MO43}ni6bacK{RmFC6NERfP>4>o#sS_EKqIZJM!=w}xrY0uA;nY46 zoN50funl^3@yBuNX6Zo^4s{)?`7UJ0-i7fgNoc#$p!IjND>GvL;FShZQBl=R zO_(nn9iJmKqj&c9O!sCg+}z!%Gy~n1E&(xRV&VcT&FF=={w0UZ@@O!8Vnp@lNiw03 z<#;L)@BDh<;}6-;N_hxY<@hhkx<)@H6qcMUAv(tiRV6D6|L}%UhJ$Y)+h3BZLUF%g@;O(cM6Gqw zawL8qN3JY^wCUCXHjlpP2)g!F$ctk$d;MPb=>x&6U*x@l-F;eF(nDM871rXa^!(F} z@~;xv>H1w{nJTdPp6ph>Dnwt%V!KrOBF-28QMmH?@75|6%M3ZyVqP-rp06vM#=IhH z_^~^H%pmMP*TN2`2|e=|oEGKOZnJ-5YL_EZNL*`88)IDhn{FcZDCE_O;gOZK=r>lO z`odsU`PPpws?AM--Bii<1x`ieE!^9c$LXumV) zpS8@nJL`1Lp%1uSqd6Thx{qhCYKsSIli-sy=H>1PIsoq}BIg5IC_k_A-?v46yPpnD z1l4dvMpFpY0%Nl1^gi7j>svc!n2Bde8TtFCF%~pKH!?CxaIcX{#+Gb-MdcAYqsw&- zOhKFM``Mfsg*`0IJLOs~hckpNi0>0U^f@DT7kBOfX05ecQY}Z;F>9=g_S2Zu3Oa^q zi?44{TSV>mhDbn1XtL%!Y#mZ71=iplRNoVwO5W$+xK8a z?zs;^ah3-;5)cwc*i+kig`-3Zj(&7K%B`D;`HKc{mQQ{^9YpE|1_ot%Fz99mAK4&tlNV zMHPgwgk`RUCYIT3*`?tR4XC}r1UZ-$2`(Oe%WA2jHduQ^e_O~FdR;wv=p@K>{s3Wr zR1{)#-v^cxOra(zz8J&Ww))0Oq+R!UI$|Pi!dGdm1@~ewUCeAL<<~B9T3z>9`5Pwl z&m*LmRBTBRv)(DP(EPIAfzy1LNL#JZyod&!1v}20&suEmfcMPZscU3nc}F4Ur+@yQ zxtNOhL=x2>UHYq*sks1aHD~mWx7Aq1*Kz*R`Un#dtIat%4R2mjn5c#Cply?3-)k9G z=<_*MvZnZ!Tv__rx(`KIJ}q{JH^UB38hSx$^CdpFu&~#Qz=x;JN5CB9a@B-)T~=mi zbrCS$`Qh#op_N3FBp@EG!`>y2!yeK?CFW5+C!ngTYN|^d5)x7d3({|}ql;4z-@*8P zD$c~roCZWFDk`D@An2V7H-uiywpNYzCHBU|k>3>U)|n|!!=rbie^d2WP?41Z25@rP z;j7j{?Su;teUe&3gfym`vWU}c!_amp?=dwXvjl%X3CoMMmgY4U^evRe(3N0Ca*46sQ)HJVo; zN+zSAkb%1k+WmeZnBhO!S(X@Y^zaXkpqWWw-A~^^yG6D#U4QixE?tj%R@-MFx;ft- zcy@}F+Io|Nq(ZEK0;@A+cwAf@fL+wo7C_$5rd|>FEwHrh7TcPeZ4dVL)*3tSIMbA1 zdKEKve3T_4{dy8+Zau~(rQ_>SgE4Vzu93qlWl-&19yHSQ_@nXXK6PeEt0T1>!`eD* zXl@N(J%TFB&`bAqx_7qSOZ;Jp?_J&4FoSecxDXII?8jDaLdD{H*bjl`L0>DylkgdFD+d) zEMvWQV}aapuv*DYo_hL#31 zCjH4H)&8|QTdh(8*iX#cUjk4gQyZ?l;;DjS_E$uC^!x18oYs-jvxB^K0|(j+m<$8V zi8)7D`Wj<)O- zo_g;HOK>_xyvjE4AUZ{_Y*o$^?r2S%XN%vSlKoW`oRlQjcrU}KFCi!6D{(z;Wd5!qOcTDC-5K4aKi>{xrKy?bElj5WXnToYI8o3{<$Kxz-2wRO z9RIXzB<)4aN36*R(kq3(ZD5aXUw_H?2{9(9*PQNXt!-HchsmWh%f*6umF)pYe0@}E z=(xDpGA%pZ%*7?%O;njT&CVyPTmepFjn%w5aGCRPodXTLJI{2V(e*nDp)qrEItBY9rqQ#r|DO z_mdP?u!l=xV?j((_g$I|=P*vA9*!>Cl@5~`A8~f^Y9>F#PSX;;$QDCUInvL*QK)#kY4HwZv-FwW|u(!3H&5%wb{?0 z?kRntT+-79J-c3W3Oi;kj$?&$J2n0aUF3d*2YJVh8=V`CMPC9l9?s>Ejfenoc%0hX z&%D?=Tv(sXqYoax6t66&RhD6nFBLcnU%}?oP59t`4EvK>pa;Q+gc;To$EygKB%mpc z36FBNWLAq5N*Y8PNnk%~q4TnO7X~UF+$Dw?hn!{DrtwEK?o5{e`y35xFNpUxy4|G~ zgrp+yoS&qskr7RMdpm~Pqf>?Eqt?Z9F25y9fY(!tyvvVGXu@(HZnykEK4H07>W=l$ zYf$W}3q#ULDI()BAy^&c)DYHZfF@B+yTJB$c<5=-$P9pKLE_MZ9f4=`Hj+LvN8* zdCV!bSYT!vpOqAPxC^vCXvpQ48OY>)skLt#0EV0LYVofKHHf`i(9F9Dt&Jq)-Mb7XC8YRM@IjTYNL1=8Z&JP_nFHjm% z*Jx9auzW*T?;~3Xrqc;(ex0hc{X6+35f-?HS3uuueN8q38BNp|Q))G`VHz*qwSr{7 z9O*F3zW|OqCB9cjYO5Gplw`7}Q+CK7vs8Rtyxq&F<0AkiWMb@=#ZEZG~O}20@5EFdNivY+* zo{W;p#&`7vJ{q#7rY2xm7|z~G>*+v~0jA^=F<>Tm_v-LC^^+p6@`o+fwF?`G_am?Z z#kT511!nQ{$@ zuoh?b03Bs{yFZ)mRlZay1|DS3eP>5vE+%6h0Rv?QhXP}IH&TpLz?yFLJly=mUise2 zDgtQ9X4qCXgG5LJ513ZO9(jKOa%C)c#DjvtwD0Zg?6m)fBi`O=T5QwvJXL=7^GyU( z{CFk8;DUgTkBp22$a@kFzihhOyRDr7Vc?EP%1?iHCQC{V%{9EdBi}kCIU>mgCAt_cO-V@T6aL2At3dJ4MVQP9L*hG!O`TA$!^L)J!1YGkra+ZQ(Rxqgf(J#;-CxV*{f0kC zP5*Yw(lRsG!C*-avx_ku{pISyaPhbYhqc3tto^oSXjZ#1Xx3{-xYoe=+Gj!wZiEX9 z3*ryg6N8uj-~i=pyRXpx{=WRNkZ4DB{NvgOd?&!OfPcV3aFM}tOx1(cutW=CCwdi= zOG|+YFokKo@9TBF>ruVGWutd=pME^+ynOm`kbkC1^KJS|+L9S*@#l;wJi34?HdE5T zsz_7e6|4F;RUrZyNwix)V|Khao+s}(1o+zC-hSPs*@AL?@@@EmF=PB_n$YcL%xA$2d1Av}r#; zv^V<&Rxfaf7aRd{ivF690sIF}IFw-C-SI=%vz*-Mk465L)0!6BA$QKW22S1VEy8Q6 zwWDL>+#({cK@RWk?yhY0D^O}n-|$6 zvAM5NA`9Y2b?DEUSW9K(?3^52tbYyr36{yy!!eiH&M`L)^VZbyjot34E8F}DD%EU z*mA0N_aPaZnh^0wluVYpO%8p3g9pQ^{;qmYDD#Mj3>_X&7lHC2P%GqJHmOb>SqMa1 z0F}CQf6-{S)REQPECO^Di6f7fL-0<5N!pgH?+jHj#zpryfP(?l-=BtSP+U4@o?>v3 z7%&4Ik@CNr8*hML{k%xLqgHe^gYee>-kQZ;gM^B5z21f?P(`2WTowRRl`!P=Z>QK5 zyilD0-BIhXAlk!im%elOyCww4afJV!s|PZ-awavTfSlKx9{&vp!UlVq87ec|@jL-b z=iMQ^oG%g+2>Pn^aPcKzp@7i%Kf;reclavpCDW2$Rt zSOjuQ;LlM1rDbmkZ1jL4JtmslWefx$rmNU+qWeVO@o&271n=Pp?-umu!<+CPcg>!G z|LO1-LxgVLH-DDK5pDWh7Tn>-U8s0Ycjs`aiYnjGqrV5pQj4?KGBp42D=I3W&izmS zFA_T$dA>E3bI;Gfz%cK%7L9a!6fU)nN^^w-o_Y>C#vY{H$!1mFlXu!*(2?-=NNiUt%N5psGNMN_t3c9?Vu#Z$LKIk!Xk8qWDFCN3_u zvy)fg_)H(jeUIJViYJ>3c#g4)ZGj|=5W`n$JV2X(jSz^%figy12!!efSF2&w>6l|c z0)o!U47|Ud_d6YvmoT`BYkIi95m}4kd(6Ultwvi|%Pz4~kbv3ynNU@nfb|)IgaM%z zO=?z_#AeLcO3$|<54J`5ycmcvkJh5cyYX5Ij8B(=j1FUf$p}VL*CT%n9H@ zAT^$chBo!&zzHLWJVh?QwP>MQ{~HfOFfp#c)lq)i=>elEMql;KWnM+%TnjuaA0J<) z%P?zDdN0`V1SokQn`3$b!LH8=jKbjg03@42qC*=piw6r!qI4V7YLAX8CFDATE1w|Z z#N=c|jH9Zr|EJUtkb4<4Q%X`f6Ru#=U8j_Vj*}bP!o&{SPE$YQehZ{|R;)97H=A6-bxMAd)u)H0WJh|Az~OAb!+^3Izyf-v4qz-hiu33b^(hs{Ln3{5KnjZY4mh zETHH6cgzeSN05e#4Scu^+#2PZup{2Fd;*4qpl40p4dHKPy8ds~3-*`wN6)>G5&v3o zwhtR18vs5ix_0$Am>;q(&bmnabO>e^NDwZ$08{^BP2RW^@V~k293D#5g}6|zQXP89 zpWevQcLPi^*#(G(+uY{2Z(lmnXwr+{m#`g zsQ{2{W;kwwLn>H^p#{4`f{N%zc<24O{KlnD?>E%i?kJs?q_Dxn!+y=KZyB%)yFoD4 zgM)(z{miBn7{w3wtdzF(XLa@U$MN7zL`5J3BlfcW(4X-MmPQE?(7lyHWEz}%nY`BV zPlca%#70`amG_0jM#_+(Sn!fH*#Chc%HaldEwlP30`E^58ZHFH1;7FXtz{?pr(&hR zH_H8I6)OZfe~_i#II+ldA#R9>a;)ylkB^=2SD)Xnet{!?fk+mo7R~b4UQecUoJD0s z&^Y*kK9~*xG;%+AUi#|y+LWIM#qyP(4IDIYQ3IB+c{v#8@l|uayciFlmVIfB|fg;|ZxOH6Q>Jt6P zSqWi6aU_Obx$~`_Rc!z(u4gQs6*2MkT;AokwaBl{=sPAMTy#cbRm(2hb;o79c)t~hFk%pWwYD9&SBcGy==BB&Ur3hN3cp9h ztOyB!+rx+Z=U}JL)2lcEco1ib-@BJ$@u)82DL*;GU1OL8lWq|!ldI6r zcFoeC>^FO%_;xi9U0-uoSeH0$VPcL~)F2cWiq$5sc%6TD@WB3_0w2Yc|4;g=k}DJiR$J)^Pb4EEX$~Yj1y)Bw z>082TG{479#8<)itEOER%lpyA4@S!zbbXd;{m*eI-HvTD(JC!F6N2zp%{7MqT=>T`u7r=b1;294?5h`$eP&80n966h*M3 zu<0c4Bet^f<{4-@%-RUA_}*6_Rjsg^Rv^(>S`pzz8xuxF(;whm>z1ZJwpxkk%!bruKyjd9MKk#kI zX@v=hGtGE(4bok-;s&f=o_0Nmu|6q_fF9@v-pf6+R@$#F~|;h!Gx-fy0bkAKHcMgm{2d=QE z5I7Y0fG&~&dm4Om?14Wdp2`|}7lEI^i?-3g|K#rK#-1S1B|5?fJTx8r6WA&3rDEiz z>jv}kee}c*Rsx;?QM%iEQi z0D-tb8qhm>ei>WS{$3deLEArYhI!_rqCH8>B38G0uBb}f@Ueb4`uPEcmfOUu*AX+5 zGTx5^`b4BOCSTX&cuZV*?C4gHBl|E6lZ}|U^(^9Y+U3Pt*Ca(PK0C98s z-ve`TCl&pDQ_!9mjTyCtQ|iX7GjfgyDSR*t9MVS$(1T!yiQOk1)w%{~2OEviTlA$; zGnGb8^rhMdu}HD>4aNpdc)7Ld{l{5)32r9JytvL?PV-^uNu3a*D zzk0^6$NK$)YK0M@=4O{N!;(q9VvE6(y;2K$k(N}PZo|6dm&)UNbI+cUxj8!%TUl9! zrlqC%if~Rck1N3UuMZb#Yj@Y>#`78f?{g*KPH)y^IyCuge)GUg^9O~s#?9Q#z?4uO z&7x`NET>QByK5Ysoo(=kThAk2BjQD`02iPp3kQ)9h2_VGG@f4bOph78;$|%TF=n&~ zzNFP;JuPZ{5IdUfc97=a-?&{=v1R+U2Lrb@7Rv#NS6#DWFYlU0ai6%67!WswqVV$b zDoCP$P-gmZWOOvI=pauhBiaVNkVp*?y1+pqGZpj=!*QePN@!7-7`P~mlc?Bw{Ej$F z~+Bp=Bd%MW<{$vJFy(gvYjy>l;rARz&CfXIxLj9UJ~rhs#A2{qU$DhT1iUD z-oEGpQCQqPy{@e7t5v7;TNg*1okoR5w=Jxo#_ZLGg z-b&<}lkPX*bta<3B|HF`*~rvKw}kLD&jPT$kvOPF_}1;wSwWi>Jl6YHtw zrcAiciP7q!-QA$;Lex9Wn(M~F#yp^vaOQ78e$iX_lW%pXl=B~X)2VmtP1(&$KYG=S zc<6gFh7plIjsBjHK^a;>K}?Ft*6;mkJTg_^$nLwMH0xNiwP6#W%}FDV%I4|pclUPQ za^Moq(!Fc!v8!!LPutDpPYbd^u! z(e9q)*vC;G^yyL`wk_Ocda|=+OW{kCVlWk#-1`!L z5?TbD0d#|cfrTAXj2JU(&m8**b}v7_X^LwO;CSSrNu@ctelvlK+gKayo``vV+j|M~ z9^WpfJlkEhxw>XQ*{pyKdLo6IrS9eNs)@Xof|}TnBn-a`ehKwvY@fR6Mkdvp4vtu+ z_2IZ`1q|nyfrs{&&Z^ifb;&8S|8aKKJdk=$;s$M*E0TEbX&d2AXUF=rEjN}z6zZ?| zs`4kxnmX54zE8?k4szw?eqRJq&r9@=c!3xRJZZQ!l^K4Tvq1|i{K~vN&Qn@0L7a}l zu>bvt+|%g%2qLaCM1@jxUKp56ps}_{0T~5>x02QVxh23z2G48<(eynXw)Va?q`*Hl{b+T4_o)=|v zx6N#rjUv3j)mu2#{DB7-WB(a#;`{qx@%DHOPxV*B=N1VAvVPCD$GI!@~+jLa#CM^W@Xu{D8P@tl8t6 zgGaZe(K%tm!fEF8&5S-GarSC0dydUtg^WcDEh-+Tl)zOs`fK42s^w;I!{2-aHzW;| z!E@M+-p!L6(>bd)n;%()-ytt^NxO9t7~yQ<{9aV(%g#HMU4w9hYl4F>DggKG z&BGfP1};*nR9|Qd8t&bULzB*>FS=gku>BW|4prhhE|OVZ;GdxV@wK?Q@_G77ipSnh z9j*+8E8mZaPYF;URlC14d*;aRx&CF%rFd!LP$s!%43hWZX5m<0JG@=C*l>N?k;c?BfR+=%XRqpI}Y zmX$kfHT~IbBc`iw!qCi=L%`jME5Abb4nPeeO2vpb5iMjU)PW+T_(Wj?hU->8GT4}U z%a&_rL&&i^9fS0Qn}oPEx6LXXc=nOWd*NeuNS)yMh3arVIMcwrIqJ9L7%Umq6Q>KS zVXB+OydXkN1=-Mb25M!5UJNO2jB)KjaTVIeDlipU_s^ERAYY%_nTLUryTZAj$m%V5GT&N^7Fk?zHQV2LB5^M6#a0L?ZY%l6)ak~lnz<6+# zR@pFs-Sq%|>#2GP14zP-ftjg1zaW`U zz9+LZrz*b8sciIfWtTKz8Y2UX;Z6M~cAkvBjmQJ@exsR&+JtSn$-CPF;+iH-YHm+z z-3=J|fSdv|fRgTVXIBK`M{-P!BCaLhrRg!oeHH>R2%(>6$!b1qqG=ZsOvS6XmMM0- z6kwy;KDUL@;|EHs;}`q4X{Lu)w{JCdmxXsa79-yRG+VHB0D$L^(fb&?ZW;6dQT7V_|y;|Cr& z9H~6x55Qb>MMJqttGCAHWwRS(@p!Gs7j%cXs^HE$@R{2hotO^;-j_;qoS6#^lWqOJ zb8iDz(s+p&G_kzn97I~25`lI0vObcNrJ*bb?lLNSx__kw^(t-GnyLRi`81!m78O9w zURAQ)*GWUwDeemi*d9-jCHGUY^AQ%sf1jkElf`N!49ww=<#OYY|6~{cCNHeqPm5~n zV2NS%i=75NFprIafhDE)-S-er- z^mIy9EKZ{^ZxYsc>zk;~Wr^9vWg{p`5GakQ3xl~t~-arr|uGQp`XKZi*EF27;{{^DDr zqZ{S7dhv|ivW$F@!mdag$AIe$Q<}m_G5!F5s7K^LHdq4vCv&Lx1!t3x|&MyuVw&Gkagnx=V z54%R|D%?E*%RT5YWzuxbJ$P?XuZuLQ_4E~trKqbZO?IAiQp4l5Gnd-KJNF}xjgtLrC1A1GSefpy4>TwWh~caT~l|7JsSP$s6aQ$inEk>o*APUJ^vBscU3 zwuBB>ibBz(kMqEM!clG+o9!4Es&cZo1qFQMZS;vTK;G`1^|1M!+fn6(ynI1TwKtTr zxG4iD*M9beevXjHOIi*Jk4D9w=2YApyKeiQ3^K_vi1oo?D}?y2g0Rgqm@>}4BuY6& z>P=8<_$Jxv=H|W0SKt|x;c-Xw-?WTG?p{4CBiYae@%>?Ex#r= znzSnTMz%Q8k?l=^3gV2GH&%~SNlb@zu#5;U&){iEbxxlsfWbJt7b_}g&d|S<26m5< z@6@W?uv%q{#9Q2MEq0?B($i;?xtD6*mODOzNjFyUujo_muhu)|{6{IN14~+khZosJ z%*M?$V#b9)JQ8GOZpj-el@4%K-!xK-d5gu6ccBxYxnEbInTk^%)iV8Pw7e z`lKIK2BVe*f=S4Yrg9Grx?;JtBo6Gwb`}v@9<8#*dI42WAx{Je06V-wDhiHUhDdUG z$E{5IHkZk<{$;5bKgPE==RWWFC-=L~5vM#@IZcsVZHSNtvfuU#P{+5H?xojx2d|z) zoNi})w+D#vyKTy|^vMfEg@YDV-166*W79IL0Bywt1>b0z*>i{t zpbEO3P3zV^u==07kPuFO6Ovb4`T|(+vrV48qi!L!qm$*Jz~d<7ZcKs5&yyi00(Zbu zS#v;_6B~aA@F_12wwoB)+qiqkMRTe~BQ=H7eVjS_d17N_XhtV-?YF`=94vy|{*D%& zj3M>#-_jM%@1zy?9DO71-w8NBl0(tZz~By!09z(D8kqi$*`D>ftsOES4q(&5i;!<$ z>|3hBgF3ZBoIIb#(4%Jh_FpCuPA?GyEPr-Wu;}bjuP8xth8%T? zYGqTw;6~)Spy=q&4)dFA!AAi?(m)Z`dH(p48d<+BT%G)P!m*sn>gbl!gZ4m&JM^_7 zTdWrT&g&LZFGm3zMj3*)foy241gk)kx-mza2cE|ZQOE^d%aCLM@+IW+O-fyt+2C&L(-qiQZ)~{FGn9+D8WWaN{E^20}P?(@o1;F%Iz=x|I;dpH4AQgGkm;@^xB!6o#;E6_2=vzDcBwt zHC>l9kdq2G{coxVEUkQ9n7h;iH`DE74Ra2a1c``4fszFzhH z5s%h^28d5FzljWPup}-Mui_)ennHPHoS1!g>m&!zKG7T6tq<;& zdK|m@z9B6r#gw^b6=>x>e~;g20-uoM6pq=%TD;|k+YoxW;9qu9V&|Ihg)H?QNfl^e ziRP7}AMG^e-OQug3xPYzg;9{4RDfX(h%fWN)O+=+?4ZwKoj1TSa69IHa(TmSv911@=& zDxuJab4hE%b$bmyXPXzq?BUlJ9|9ql-?eQPLQbg0v@H}e)^eZ)Ko)5v0ws!Ny(<_e zcSgbgtlu6ky$vKZtws(Z;HRhHXRx*n;MxI%Y&rB5O-YLIeU+!ol`HeRgxxofz9kvIW?(80?hLz6WA_>WA6_cVbH(E&|0&$7u_~Cd z5y88?20&4{^h%n*7TH6a(+4+MsUOuW@Zv_V&A)T3coo0$FAwo zz&W}a>_*R#W6rM&3W`?6o0>)L0(*d#u}Vbws_f5REFpt~ zeQ^OD0aSl_3+*y!D{*NP7D!xRlkb5_P539P5;4x&zgyB@GU4qF=9TM=rHZozBhznI z+ts@O5M4V9pTNugAG92sI${rUc7#GZL)&@SVFuK47rpB6-xw8waAu}Q#7#gA*qv7F}7)fPq@oM)v^Y#_nv)AX}u)BhB34o#|r%|0{ z47#8?JRE`wzK@f5+c_Q8w$+Ok^liCR=ep^fIxd1e%SUC4SA8Kr9HH_e_VV$8b6t*> z{ap?4osEJU>%J2t{)e%zC4e(-w|>zIDmvnvyTrT~+~6@o*?f{Wc66GKygGgD6P(=> z*Vk%;_oL#~Gf#YBHEA<$s;=&o5mNXlGN-WC7gyYSn8i22E$Fn&7DB%pC%ER9Ob zyl`Qtp+6O{j{eWQd!Y~IwHHa%&pbqhcgT=8l7C|;2q9j~=mkvSiu6h8czHGLLhtMb zIwL6WCCmqb#nqTceHogw-XKC03hS_3e_K$Nk>#sZ;F|qw&)|&Cr!1Q_JexWp;Q2%m z7N0BtW>G%_G^tfcb~6zjZMQ7s!#!yUX}RNAja(r4y6NKI6zo$2H3W25%SJzjkOG|g zDU3I~mdz$Or~!rqhVVs#G0*4sn1?DA38!OvFDziYiaD;{)FfO-MNsn+tDL8(2%eDo zwuxUJvV)_ezZZNN;wFq%NHciY)%}RDIU^t z|L8d6ZgpWsBgTq=>VKINXbdp}pyR%5)~W`gg4M|9PH$nSZ46@tKvO|Y2y?~40M>y; zEKhCOBxT7H94d6Oqxl)4m-rBcN*}MoQ_)mjd0SOAOdN(*;S*E>VtT1FPE#~Empupm zpqG%1g)jcHpOhRue&g?><~RAHw@|?Wh7SK>fh{7)ryODUFg1DtMfV8`2wO@2;0aCT zr2{o7*+Ef)0Oz0-N#)+4&!EsJ79MXyt3GWiQh8k(aDws8A3m9j`5eW|ozxOb5$`$9 z0D{&0Q9JjRV5YUf&SCgo7xwa=X{kc@5^nC>NNJsyHK1tke|kn%cPlUOaVH^TiAXap zt;bjaowB0qI_P7U6AE=GbSUq-RmKP{Ap`~PEnW>JL|0tCC*3bVR*5C*sNW-uD0n*FFUfsr`_jxORgO>~aK3>DK+*rto4} zRVWb3FN{z6h-23hUv^u;e2Ij@xCqUrp^t8~wEjp!P6BH4`Ei>o?4NFI7P$q`iEsr} zEJh{42ntQ&9{s->P&^T>+(T@>iVQPYsnDo@1mLNH@a}-aKU~(1lU6^_|7=;sko-f&2c$BV37p(K^EhEt+EPhSd-CD8w6L8X|VAf(tAm zb({)KtNe)8k!pu(4gx!LTsT@6V4+-?XQhohdgAi`ISwk2k6Bo3H|xzcT=aNbG3AY zJ|^$)+(luJ!6v;Wyi*p`a*0PT0XSs0rT{@q$)~_K)5oV_2Fuu@(2wX-SPk{JCF~Vl zX#-kXja(!Cus^$)0Xp2=N$_9nrcL0V%IA;C__>p!fah1`*%=JHbHGn{7XX-cn?N2} z1a|o@@aLv&A`wJ1`-X@-{$7s-aQSNmDqs$Cl<+@j=wAl@&rJ#KvG^jfrwDdZcx$iI z+)=822S6iNLTv=tO#aRz?PDpug!Hi;5i$vBCe#@}nzK7<=t~WYQQFyFBEH#l=b_l{ z6Yfd$q1TJQI_~l7ZzlgS&$;sC``xC6QsbJhLc6niClsqcSWa?FUr{huWgFEZ9GZS4 zw?9mG2m;6IOR?AbPgL-Q7-N=1oThxBV5J{#tzZrqDepfDq603%U&Ux-|O2@DpA=Dgeh_kK{= z`6?b$a`lL=oKrff1gAGc7K)gr%)snRKkWT4oc$cfhHW0kJPYjy&uad72@zyHAo8ih z(_l*2jxH#H#VY`riJrJm!<)Gfh%6!eC1LjmKF8FdQWz zdO^*vN-`hRgILteTHa!VEu~2mC`sz;OuM!0nj+4`RF<$b1%c3Ukr{YU_5dxc zIaW%X>PT{q0%E#=?Jw{P92Z^>=mnVFQvI;+2_I~Mby@u8>eIoC`_TtxDOF(S^P2w>m r6sd^>PYb&HOUeC$hX1*F8shM{=i@8*%L3pPB#?%R7PL&+>goRj540a6 diff --git a/prog/python/qqgg/.ob-jupyter/800bffd38987eae852b74eb2483698169c08c4de.png b/prog/python/qqgg/.ob-jupyter/800bffd38987eae852b74eb2483698169c08c4de.png deleted file mode 100644 index 67080fd7d9f8a2cc5dafcd0c3d2df5379fd84aae..0000000000000000000000000000000000000000 GIT binary patch literal 0 HcmV?d00001 literal 12003 zcmaKS1yt10x9-qL4v}YrT6}Ynb?*vrp~)?Qie1BXqP>NQoX0fj}TqHC5;f;CEok@Wu^c z@8o80^OnWK(#6%r$x)b3h)jU*++V0x?qqC#kM3A2zhzX-BYV3x@R^#{p#sbXbddGCG1 zmcenMhY0aZGp4mHrbs|Gaw`v4K!`5)%5^%*yw8-WQVi`LoS0k8RLIW1#%^@86qUBl6OhCadgdd0$ql z7DO}nB*cUJM44n}{$+5E>w6a)mA|b)g~@diT2u2Y`l|q4&yI$*%WpWLxGb2;1et)w zhoH>Z>s6LnnTcOva1gOM5xo|s4EQgNBf zXQ(y&?L@Xs@-jQRdUlQS?r1I4B1oy@8fpZmsa=?2;pk`;3^Z3BfwRYyg_087y0TuN zHsct;;WQi#oa(E=HNe)>KR(?*I^LSp`@(O!g9q+0pXKPyzcB3paZ3IcEaHI2i!nPo z_bXF8g&GvtYX4)&;Y&lNfn@NmQc=v{_;F}%ChE6Uw+94~i9#V~aC@4f`=eB`_9)6^ zzzqS7@;UexD%GT_Px23PPW#CPQ&8chZM&z3t%}=O|0qBNoY~VDtM}dpuMfC_<#JrR zf`tCIatU3Nx`E4Ju(+@5w%MqwpX4kwIw_3&Ki1(i`J z9<3JHF8KqMzxy%9umDM(j1c&V2?(3`;lDOH^9f3~sOlUs`ihmnzYPy>==|i8 zQQSAG!v_vF%gox!BBkpK80zHM7GAum$WG^Xb+mipU>iQj9}vVvW=~|V098s6&w!=x z{sE;08!c>lX%@h0Yr) z=OpLMa#MfMt0qBMIwxc}xB(>vgT5xH1I>R`iqt>$GdX|SG~48{pa%07Px>835n5Dg z&E(fy#Fg#Ua)K73jIj^h36Bm zU+(1)(Bdi%n-DEpcz{_iGS!Y%loB;~mgTZ8a;6^e_v=6){h-}~L0%rn!B~)p>17LV z3bucT4m~=3=vnLcp_Ugpy&PCOgJ`LE-4}1X#_QJJJb#k^l<2faP0F~G>tOGq#kiA! z*72qE(q2HW%M+asS`((ea$#B=@PxV>IrEQQHf!18-*%1qL%TVaa)=G{9XPuEnn*bH zwcdBxFlH5K75rjub;WcF(~?D>&~LWN-rp%iBTmf!P*n@tHRJCz)KjtDXS2pN6rz{S zApWA*NkHp(5<)r?z7O7IKTU*72s=1a61d3FmS(V+!X=8Vn*F~z$V{?R<6KrDhtzsO ztYyt>ypwE+4`gU#We}L>5*E0`G}_~hyy-*6H8&X?#f!bU?>b<-0{b;q&Plj zSX;;ydvx8G(*8Fysvr*;^*aM6VdcY7f04@G@8E*`c?CAckVxBiy{7#QUp$F!r0f35 zGR-TTxRup}JYl0!Upz4Y+PTmj_l8!?h;hrnAHI^f6#T`I;>4FG&aF-M?K1b=?aZNP zkp7sag0{F{EEI~-ELN3u#T1ie8>|jkSV$~qIMVVl0!SC{{m_9lW<)X(Q;RTG_IJiS zU#6Jy6`88>SaTPGyY`0WlqK@cbloU`QGdyaPvK~iY`Ho_&>Iwsg@&yy zOtlDaI{=DfQDoe_B~082MbjMp@)ZB_mNwtLnm6Bd=#-)*fEU5sjfyG6y!g@<-5{2f z&KNQcna@!Yxc8NP)%0j*;?i6(8jMYJe|`9MuD=b6>CeAH*Nd_ga*!5ktq6D08XZ0X zC2H6z$C?(7e|#<*(K7ie&Ekl?RrbNS+aAJ&XYN`}*(q|Sfio=>?ty8!tgFQYETM?4 zv`xmhYZ3RJy)a6D3he{E!xlwdR7m)^Ood*0?Oekptlw4JSlr+4xDhu?^^vIg=*yRu zi#vg5vL?{5C zWO8HE;B${ZkJ2+LZZ;pVuft!@MMNV?>3;+>&c#`v>)K`q57>|2xrgbSDf#~V&Fdra z_|uZZx9*~)p%rTSl3@K~0zsM%A2`$l%TE7TR?yH2l!wi4Ufee@cO}k={ThcNep3LO zC|$`S5v3KE$332n7~7J1A<4UHHU1ZW3tano@zN@pbU}dV%l@*aP)wTfoW)zN3LpQM zRxhC9N9?SugKUaJ?*00r5|)B=dgh+ zM1?O3<%qpM|Cgl2L+Wv=mQ-Xilg{GgYpyQK)VJ`+Y9U6IE&RcY*0}syus~_FH9$I8 z2ZdfpilXUUp1}2%;uA;8GUxZ1DY+JL3zadk$jt1k+D1+d7zp3kII54~0%)N9k9+ z=f-~ZOp4=vVX3v|z3}KxlRx6+9QU1de=rES*uIFD5ZZ#5_>$r7`bv4J=5gHD9W8;l zG@2aDKT-xMG$x9Qm!e5nJ1=Gm3oK9 z4Ad?&5(g5rV{tB%dXKE^STj6zH<~043P4TuYW!&bT71*0S|@ZrG>Y2)FdMNR)AFZ# zl`;<4yWjew{*WP6m>o3NIPFTyz-ofTK2}G18%v8C998l5!@wfsRgv(S#w2F`zQ@vq z_(wga$tC6(JJvMW>ANjNkrw(xmpBx@A$zd+q7&x#H%Kr#qcL<;?YDSmjZN6mDJU|i z#aAQ`ZT#{t52c<@l(|CNU>5R8mE>cxWRtcj&d~>J>?|NZ5|U;k5^fsCl0J%ikORF1 zAy~tZMk786UIs~jkqtckkpyh?(|S7#f98+AMgxCjOQ>t_@&nEXfvVnPkqXknTA~_;7nX0i z^z~JG7DaRO#E}$HcEzioqF?VDSQUh(@}*7m{c(-@G`Luc=mQ++xluoxaf8-EFYzkH z!~WjCS~ScWTjLOmw}s`JPP<{MPOl~4wyGti)($-C@-|v+;QVS^`Rr(m<2=3paK`v2 zo*eK?>YkiC_CL{5!t&_}Inal=unsYX)PsHJ+Mi|a4^<~WHG&()?wVd6XuGNWpaPsB z0eL*&e9;rMeGw%8X<&hlTy3ch*MB9!-LhC2z1s3U{Fh$oKbE-FmOsI(E#9i&3}O`` zKKR)|G#iNJlReecFl8x5gv0Ad&XK`P9!wmhkgjhx38-E_6XVrM*S zv!DCRAE$wMp=|9mY})m1LY4<=%40_cPB^3>imafJA+dYi6Tkg02Wfa;+7rnAC@iEO z%2A0!re_>h9TI{a*b&QAQ*AfyNH7q}**X27GgcK`xrdt0pp915$~`Raf~s@qv1eLX zsQ<8?^oN?Pqg|74;LG2K2I}sYG<+{|FOaE`jTMEn3lpPH%xZc zANKc+6|sJoT{#AV3m_&vlXhDxyKe?;swMR4+=Guy@or=7FuBj2p zhlSAzZn@^8qodIeVo@1cJop{Ss>C15&lprS{KPGY|x#D*_k zP#!X{IFWLKX-j9nHcrO+Ek@q(>iwwp`mk0t6dD?;Z(tzgzMW(HE640F&(&0?-`~6sS9&sJQvpq6_CL~n?X%9d zQ@2@K?~eKO8L%a>HdhP~_z{&&k<8d!I{$l;=GrG&zt=*1b{+6K@Y>z6eqXuTah6B{ zP%~XvS67#W*B;9+&)s*sa}9URKwGoFzJ!E?WX)9cNfE=oD~3-M{++IIibJ^^{N-?6 z8+hcgJ*A|p3#Q`G$}KLA{>*LY0<8Km^>-H^hr*=4_v!Z+lrTR(@S0Q$c?wi8)N+HUZI?1m**k^--S4@%{~66;N`!_F zuIlbC+An`A1MHSnZ;ELB*}tfxjnOFi>viqB0fAxmMh#5`U6xdTU+zfXT#PsVbefmD z>)04@9U>)WKhd=GH{twjG0ywud=(7FCui;fepsKb;m`Cxp}#!Y&OxDwHcbQphF&!t z6QSaV*CJUGspo-tU}1_5ygAHFb59y<$kAG>gN+wEyHhn@jQ*gLK$xlWjxLnD@@BV_ z$cw%V`r&?LbsZOYwHVI#s?HCq%;$LHr{ip0Ra7GXBEZ3`rCMpg#`cc<9fOOVp%u@@ zrK9@#7#7bL-jj!jjyJz2SfsV)lvsuzwdoBQ%-TE-{%Y(x6q-P7U?m`^4dYH*bJ=Jy z?$F|s<6H9c_jtosyF_VFY!vU#JUPUmr|itVRp6wH);F*Mt?@!~1Xjia3}Nyv*&=@U zw?^?_^dXYvGZ=W(vpv~y}g~~G`7-X`K-8GI9*oW~G0Rqq7Ua%ol!G5U$E!8z2 z4WmE`Cfo=F4Qux&Tbjs@u%2DtC}#%VoKJ!P+p~QPX&o`DWWh}q%gO~fy3VNO$}TzN zCy}|!g%p0Svw;NQ}b6L0ROLh}S9f@%ugrQvC_@EM#;=j~JP*OK0xqEnSTg(zN zzsnzGA@k82=edtNK6sr##9*{BHM|dmO5=pZ41Bazwy>=PWvcU~rQTtC^-=q${qL2n z4|VEj>^+wFlokD7rembd3XoR@&;auKJk&$z@O^$Fjtwi#ELqd%(qBsXY7onx4G(5L zW4TWY&Vb!s8_XeyqUG00*QYe9_Yy2JD0wbn>g_EFK$^$M)oSIFJj>y=qQxZ(YBB}b zD$ZeRY|PJXZh^^}E$=aGarP_@N9vs!J@t)H;)LkvJNpJqE71Z3_gVk$wmtP(^;Cdr z+|<}u4S>`2xdz`JAz-jP0K2Z%ROFw*k|mtq&TPupo=HehPloo7a>5NL216^-)gNyg zB&#K`-}xHsV31~t$hZ3Xw0l(SfC2OG&A9eEH{;O`)Vf-sP$(e(h{viDYjLtQsjsiU zaCk0pI^%i=py=f!6FMn3Yh)~&dXj+O8RA`U%AdmgjIOMw@2%&5%L2*GJB$RgDwnyf zzewUY?fWWdj=tC^^1j;d@IJ$Y*RE#xp50vT8aD>W_RC%J0A7SVo~cO1ONL2QR5Y4S zFeHiB8MH~Y{hOYz zf-*+~ij1;8iR_vwh_k&n?|QyAUmR~178hF_Z;aZ{)h9oFk1q(!x6XZ2SxE_R-tUlf zd3o7wVM}`TYPE{9B@kU zGS(-`bMUdTjxVn_PA^3%2-g6C`Za-h!3qjE!9g1z-j0i_d*p7DT^w}t+toPE9`%aPj{umn*%m={C#ofs_!u~g5e1g4`GCt3 zkM$w8?ATx&f(_uo-Se4R*W)jZSKop0{Q-PD{-n+~2(;xkC?2pYtz|dLD=X6>XOzeS z=8~_GoDCcxN&NEZWtIEp*v0XB<*zX0QWl*b7Y+f56+1Jvdce%99i|N$ z&;AjU;$r)Z!VK5*GXqe7obK*!cc8`7SSP>*Aa1AZevDSxjVt$vd+xw_`T56z;Q>Eb zadX}uI0k%{_ShQG)0Tx@2Q)fPRjN%_N*Tv7N#_A>V|c8@grNJ!>xK&i3^rN{tlG(_ z%62r#de<08N^8N00k>(TN8dkqS#3{0G|U1acdhWT5w{!Tnb<8cu6V~N<0a6U zqZCCW;$Vxb)-SP0lEh<7naBlcX2$S098rRNPqvz^O2G_j;m8kbM!!6Gr9F2nHbx4O zgs+NFns>}x!IgS+;J&i{+bKbr71q!cDBE`KK|pOA>0|ybJm@fBf|>tPe0bU?UuK!uwZQDG19Pro~c3 zY$zc)CN+-wFqkuN5as3_WY5*r?D3rj)LswV~Xl=#V)1o3B96Y!uuJ6ABMFUk1Vb_WuAP_uHtn2hh&u7pJ@ZHu>rJ zP4TQMUi*J&f%7Bn{=&a{*e{=7SC=yW+Bb&5sW|4x#AQnr;DaYje6KH#n-A{XEC8=! z!2B$4uFmfyi#W}tS!MbsF-d!r1zh{A^dwV3Hm?r*S3euRegl|xn#YulA3*VD0EG1e z03g=wj!x<3DC_F#HUk_(^5$%@_U2;qyxa8V89)$RSNlqaaT+{#=kAyk9x;u8|Akus z;Jbwp3kp6*Ys~->KvsMKjC=0~D0>IM9O%Ry85zXv%@*5Z>H!=;61uMRyxFW~+Pc`& zW<0y0kY}9f2aK9t#QNhRxX}34K<#uaM}-_VI_1O!dsXB3YWp-U!pnxkdxkSXlimTq z`+@Nv)AGc5HV%yD6}+Xh`p7olmRwt}-tipg9F1xV833#S-8&po=NIkV zj7y3QuNTeoA}wTo=DeEj`O|avb-;dtTNn{B@nnS+IY+wNi}bB{cs=wGt^1$& z!YNzGmNRQmoV!=w96LA&!0RI1v-(&TDj?r^!Z|ZET)G%y+Dqtxa~DcXRy&xeGK4fJ zM<#lF>WC3Kf1GpUEIu7fW>{f`wB2ip!7us6!82;aS%Ox%uO~!mh`IUUt)A z`F%&}1sONFf*?cq=!a<27+sVibV;M4$KEUqrG@t-m|&=FRuLyyUrWBYd43+DAXclu;CgL^+X;>LH-%s zg?jUj*jiQTOU{C{B%_UU4B! zfXTQ9Lpi{fHD1pA`{^#y}3wYMX3DS3X=m%nBcQ{YR{)ojyj|7n;M#iAwo2)0|*YTY%?{NofxXck8BaOFCBBt zlnIE~MTOD1XP=f>m-cQx&RfKnvOHr`&g~()oN>Cv-ids1!nW|8up<*c!z2FP&_}l6 z5J!FGN4^^xHx?Yn{`j zGUq{iDSbk?w`UEBCx zq!o*YV*f}~`g2~(M+3h4R~q|X5_ai1}F1D8pYxwkoBC9MJ0<6WIYQ+iL5q!BlbW@!hbJG7g(gEnB3-@U;EnI zCpBLZUh;u))Q%_rXa-SG!l~i+L19o zWhBCun^xu2Z{7eOJTR_XbijZsSqM^pBW54LX|Hj?OOG593 z9*VxEqo8+t)&BHHUxR#oCXkAb_N6h{S*#@=@va5|_d?`d$qi=ys?7NaTLp}80<}ty z)Vbc*)NOY4?aU&z^<_LfOY*X>!$RMoHSsDUPq$_KcKea3s}V{tkz`B# zs~Gl7XpTxBVTR|%oZN2&YYej_D>RHr3q5tZ?VZhZ*t?~=7~U0fi4lVQZSa~-TezW$ zLfK!YW^Wk=w`LG(ky^h;7E{KHD_Mw^^}VUJR?Xuk2Ox%Q4b&6G4v!gaW5>=Eb)5MM zj(MG*mBli#VSn9MTKnj+!*cn7e7UvWh|5@j7c=*K2}58AP12b6guz@?O>g{AqRYmcz8Xwo_nAri*;^wUJ%jbD-Q@@P(YX1qESRylVY zSy^bU2F%VI+s6^NAqoQ?o+~CJdpBS$${HQ{{iO zF*Qslmv+llwl!Yk_jfZSpPu<6@jBkv!>)nzhK=&)z^uwYzBOcq=CX$EU%iB#Mj{vg zRP#*8>y!|D%MlxoNEyX2j|I`M*IRxVWG)+BOIBsjl2NRZx877D9wQ{Tk2GrMYB7ag zW(^%BEiUU8=X$gYgi_YIq3EMc2NC8 z4N(StcWV-|aqMw1$>r#8CHQ^#X6l%dO6jYPg7IG)=x>Q#k3c7O;Mq7ChKe+InT#1D zpRF>tsb$H@xHhi2atn(N?p3Ij%;OfW!Aoiz`3LZb(|r8&`tz~=!RZ)?8dB*L0uHn&RX{2 zA9a}eey69gT=}xWw*H#c^SI2kH`BXDt41P5-Z{&0ZaqGO z#ZJG1N@|`xW(J+Y#~!igTgz;9q-PL#Ek}obOEdTLeJRAlDs%E{b|ku(0aZ;9^cIhH ztDh}iUWNl{qW}?tL&b1dp%w|X4v*aiOe%j~D=n#k3n}97f%|Hx;%_PAbyu9kr)vu# zLVven5SiOv<;cXSPQnVSXDx~kup8`5G8IifOV!UY@VP391x5fV%3kpEAO zj~K%lKg_>1x~<{EVAe%ra0TpGN^gt^Y+Q`We1prJSFrI;P`Oabg$?IREK#-GZHKlL z{;Dc^^=EAJ;@EE_;%0kzAo<6GZIYG7I#?fH>&-4efzZF4!m~)GWh-~oRA2hWVcV64 zsjEE7lPTiXC=+kJhP9tvhj;DGo5D=+DGuxO?`M{1p!=|aBrzGt*Djze`Y-FJXo{Uv zhHcS1Nbc>78SMUgWxL8Lbo}Z9WShhVQVdhoVXpSIc8)vVDwct~r5m&Q4$e&PZ~xZh zx8i4!%Zo!ESi1jFRa%Dv;$oZzp~W(Ba1ed?_DHU5QDcGIHIQlm2JLEzHCYXK8&=IP ze!7dNj57`4BHxXmP<`u%ToxLMo;I&XgDFTtIR;n$;ctyEH;SCHbNxIEo#*%!MIYz6 zFw`dPC*|!y=2kaO^%Pnc*_o+i@eRp?{ZBNveOnjgr?6e~S+wS~ycPmC&@mg)sP5?$g0u2B875Yc>PLVg0h$T0yVcW|`9oZAonDT+Un9ZE3MQxNjU` z{y42uRL9@H{yz_6+rqoEO$~XN<1LI&YX+ARLQW^TteS@+Qf0vcVXw5YlCWBX#@|}S z&E_<*A3m-SvmxHXsxI76q1#qorpqyq|5Ff8_H|iKFi;+JA+==6Z{{ze5=oOclNClr?FO{^&z!uJGQi{&B&ib)5n17Mlffi&Os%1=z%Qx_MbsGrIIc#|pKkn0*Ly zn&M?_FdgY-MALeO>w~{h`1$y5<&!shwpzGD((GKVoB=p8e2`y3KsiJPQiI?GTTgJs zwZ`}g8vR$1=FQnRJ$?5$Z%0Cb7f@7?D^lo3h{=GQd3J0^*#dpC@MDu- zk4yTSD{yi=45HbC00+eQxxBRy>ljGFf500gIS}n~L zU-v;R9(RjE;y8Z-!J%b#^Sbze3F;yq6_AfQ#k$o{Rj0|SRtnps8j_sW7JEn8CHh#v zT)PTI7BxvuwkkR>4sl^?uL*yPHj&)^g>;C~H@QW#P-=3#Xb+G{@Ja8|8C_VLR85=x zhY|D4hnpX)nfqCdat-_HpOsTOJ9{Tc-z#y~;)unHPz_wbk6uef2~tl7{AN-vIDgpG zh!uz0usxj7U*K%t%DE$YCdvCv%F!QBRADhvLZgaq6DlP>s&d(DBSGyc^E?Uz64{U# zS#Sn+7tB2wIrDz-tJEAu)N5}+ijn@~tXI!bh=fx1F|>?9CtDUghW+JDC9$#4>Oxb{ zgM$Nhj8HM3{8)!m<&M1L9wN%OE=fGXu>IbUX!%TG&N&M1Mn}HiG(Br_I>a zgKVeQp~G4SiEsRVu)Sy}1&T5WSAp71M`!IsJCqEBJc=~D{j2lid9j6L{UNe$p8iij z49ue~PEK3P(FT-l?pKFZ%jtNDunkrY%S_(YS3uo?svCOx$ptX2Acw)xM+vA!E$XC4 zjGSW+ofjKsfMbab4!>F9lquzEDGI);k$p|N)T~IfN7DH|oO*#?CsIJV;pQduH&r~J z;DG};{?E#ZOPFMPDx26wN)`PIy_3W2w9=CKw0}MCFk@@{DQO|Sf_QXcbpSHKDFbDY zo2kXg7~A`EtT^2XC*t55rM7(3jlJ^$P!G~f!k6Z_b4q1yz)Eao#$%j(54Yo8 z^7`wh`L@ZFD&z3i=~l|7&e?eb@*-A()=jt!-_0jPy(!mY)2|;F#()>4<-39w#BKy;}P^QRy~98rO<>g|bzsL|vq?lQ5c#0^my99yro z=<>0W(4D3Y9_PZ-J(2{Y!%r> uHRk`)B-{V74Etc=5uO+Rf2-TxEf%Fw<~;ToVH{A94pLLnf|for3;93jN>XD0 diff --git a/prog/python/qqgg/.ob-jupyter/d47db0dde9ae59979f271a7cba8dfc46be3f1dd3.png b/prog/python/qqgg/.ob-jupyter/d47db0dde9ae59979f271a7cba8dfc46be3f1dd3.png deleted file mode 100644 index b80423eb8752ff58da3c9825f8cba82abae298cd..0000000000000000000000000000000000000000 GIT binary patch literal 0 HcmV?d00001 literal 42371 zcma&NWmH^E7c~f^ad-E|-5r9vySux)LnFc6gA?4{-9m7855Wlh0s;}^zhB6O?VvaCMP3g{9S?P9D-Ul|H%kadQx6vhXAcKkGcqqr zH+NfSCk|$IW;RAL8xIc`eioMh_j6`vH*1!XJN^v_2r>v+2~iE7yo+AjIEPKc$0vPz zWiAh7q3^zZB~kvD=KU~`IJ6+0EC^9_DR@XabmJi&q#?E#2^Bzc92($79*qJXff`8e z+|cM#!c)%{VZeB3zk6Uqv9) zrd|uK%_TO#cl-u-=?^i6mjWbDFoz%-A9T$dbk8-eXx+K%zwY^qQlb34T!Y?G_k(O< z=i^DKA3OoTH*`o;wm7y@+o_pcp6r4GDY?!_vi%^ntq^GtszBZKv;F>hT0(fw`{Iw5Q60d^lI0lr5!=fxz=XwEu249gPl& zf*ZqYVkAVXM}u1J{iNKP1zzOCJ(Kam*$(gjxcu1lEjWDr5TEpu*rF6v0^1v|A=@Wi zO@S^}Uv~a+t((lyYR+?{Q9=iqRSN$*MR-eP)2MFth zQvm@g7a2)g@FfvPaEW2)kRs+Pa_Exix*6mV1Fm_7{b>n@g=rht=SO2a5ys;jpTKFuH8 zMGSrnb~zAv*FEy^q=kDHCC{V@9hQ;V>GJ&M*ZXT}X{n%{2O(K75Zr@0N*R`HN>Vx9 zKwoq_>bUEpMDcd}>6Or2&Uyc*2Nhy*Hs|f<>i`vv_)D{0Q-Nf z5o#a4-a?_q-(G=w=}nSr_tU=paIzrHjkwtTAk0&KdxhP&f|HI4YMpjkOWYu#0+Nb= z&6Ore2w)M@KPHP*iNCj4v4A@zShk6IN^E?E};IG`NN(=QUdO-y7)zm4An0pDqay+jo3V$B zgQKIXiw|@aoE;uc+UkP;y8qHVa=kMkfh?n*BJxSf_W8$r_8$I|b2J3?i+`=1Z=O$D zghhts&zTw5mpjWl!QI6+!CC)5JK^tU5(CP3fd`iiMh*B^!@zGra0gf0-CJOLW@qCI z-xr}G>t!@Emb`AaN5_dm9J@nK%0YLvQ4ARbW(L2Mh;@+&Ee!RgZM2gG{1GxCc1=Z=2f%Y6mgq57AFt52u)%_pXFaXCJUOUeP&C^8bB zCdT2vJv6<9iNs(-{ajP%8Owq^d)(*!Rp2!E7-rQxxwS>4XAid9}BNcL~$ zX?B3J6Mx@?OY$y#u%PpR@T+x(lt|*t1|45}jbItcMdFq`v5GUU5B#CGVr1_gxFp)) zO;e|;NlV{Z>M7avb5>ru`>B0yzW`=?0}2DT5EH^$Qh`N_XEz)6^~fLA9JV4drinw8 ziLkmE=zzh#4{JQ(RCGvE^nXFfs#wTzZjX!HJ){=8WATJMHk<@`0^SL@Mn0y~US9Si z6LrRr=4iGO4-0g9+Q|1}7)FQGd3;}HG+)hXZQs^|XMYAiyiP2v^*#QjZ6vg!aL4b(dI={T!uryhJ*qeQZ@xi&^j^JM_XIFGw7ZF z@bEA^rD^+y6^dKR{G?xezYG;E_*xO7R1V&bYRNA|7L!PsB3B-3_er$4eP;kZ0FL~9 zBhw|}`(Uq&$L-$mP^v;7F|;j>>=2=IPb14hP)_<|P^yda<_9Sjc-aFSp6q!UM>lz- z)_dyW9J7!jxE`%;1rdJx&afs!AYdP`_votPEDrLpuQhnFleoMo_X;9oSmhN*o(w9 zDl}jw^$E|gQl4?Y`#fd2 z3Hqar_VLlQMoEN(YA5bR6lK;Z)2Orcw~It#qgLy@G4;{EvD{x`;5B5$PaadUds~Xm z^LHb|xGUUnQe-4_G3~*@=(i=i)-UzfMF>8gt~c$`dRQ#Iq(c3WIw{CgklUq7Mct-P zwR+?51`9tBbz}2#Bx_iTK0bnnUv59~$Uzm{54NMGZR8?!3p+}C!#ImO3{rp|*!y>a z174UbU2(qxZyQ7VKnqaOd~0m>*ojk zdh&8*3YI=M+E<~WUp@r+b6&&4G6D(a2fQbD$J6>5BGE0~+gFDhmz3bmAv61$tr6k#=(($RNE(|!)OC<(Zg*a z!5(a~SATji9vJE=u1(hOn`W;|!;9jXgCpQR*i$JeDBSF$2KzMBhIMjz0jV#!#6gR~)tC=ec~T^uxu6KUmZ4;BY2QR6ZbumFv_lkR{QWJH=}QS3GZX z^Y6TV(T4MR4_0Ji^c!yCPxRWvKJ?e~`hIXw+oL4?KueW`MJOAzFSi~u{d2U}s7-HL zo-mJUG$KZkN}bQ9@3!&_3o4uqUba!r!Vge^{jc@&+sRC|N4Fsu9#%ch)5Ke|0d84| zB0{2$&Hfc=4UA=pZo2H&n~JnX7o1))kSw2hWaZ`#yAm1sol(h^EhN+nC@7nlBiY;A zOQ=eu0Kb$mtz1(SbN^amTgpEtBbL}2OGc^UhphF4Dnhdzk#fvvJi$me3yD^@gN|X& zD>S10+7kLrRBQUk2z|(?Snb$H`!h0k+S)d^wuoTW(_5j&;R7FJ3o9K7j#(r=4)#m0H!PfxU;k;U@{wqTSmkf0nBOqF)eKE1n@1-Vhlro0MZlqODQXkj`V!dG<3(rN?<^mf}%ha_0Wxm zMQ$fn1W}F<1fjY>xlivCN@WtlIA-Agg}En|cGY|9DN0mM77g&KCc7!0wQ=FzoWx2L8<<^^6UQrujA>=>&+`GwajzG43}~GBx^LRuT&u+8=M{lNuUfL>sb?=dP)rk z|71^?!{#L2U^oq!(I!vzb4um8Wq1oCP$}ZVUg$-sV^*FmfFZRtvl z;$^$rx4*O3&m!*#Su;h+K{K}BTV7|M=AE&8&`h65N!u898tC>^k$m95z$!f*n;;$r zm)U^#p22#zdlhf^?v9_q?SQ0xxMtNMeP7+M(ODkQN%0yfwuRXEm}aOP86p^ z9^0S%1|g$EGvVXIKKk?nyGhu6@9C6mtrkfUj*aYBuUUnT*Y_44k)1bPVcLvdLtl}L zA)z?7ykj+L_TQdR#UoX+S$Cu0ch9l0F`tKJZPeox7Z!w$Z6^*M zo~Y|re>ds##B%B5tQq5-j__>0a2bIWR2>EnU!<5e3xw!G8Lhm5&4xboavjlA_$Kxk zn|f0@>if0|m_+Ap`HpSPlMt55IJS6Z9! zflxo2;!uNb!eQO${Lx6t|J!LnK?S@CPa3L?zWI9g688AXpx^)4j>agSB^9 zU0HcvY401Q5%G1F3qSRZk{d9^bGzK&YGOPDqXvUcH_P2&fcLP{m1Qob!|AVlb33~~ zUt@&k^`Hih*vf?w6TbL{$#${&E`X@~f>5her2av_88Lul#qXDFtK* z`-B)VB++NA+^3uRrgj2MeROO_Jhv+jgaQa2_ZwI>0D=Mn(^FGZ8J4fvjgPEgmg3=j z)wWg7pu;6J27VKH1r?WQnQ7_dR-8sDZ;frR4OfPajlNx$%@0P1K2ohd0Etx=S3S4$ z^r=wI!!oC4bA%&r_(p>iYWHZ&TPSp#5;Ijcuu>Zy1F_~$Uoz>!2e2Ex`lE2k*mmz1x>7>t!=E;e)Y%d zIjUmAC3M2!-WXHb;r_Garo}$Ah1{-s_1AJNB}bg~7CAoL!Z!>xDqi0(3>1EOwq^Xh zM+7+;v@ph5S9^{u#UY_t4q9b0KLAn^Ubue#ply9v%Q?ktv{Rp=%Nvx_?OArqTgyR% zesWnF7DA-Z_u;_w-uwn0p2d6kzzVCJnTT!aip6#J>|Gg1kuqMdi?u)myeq20NB4!t zUd*el&S5apCDl>$!YUZZSBhdq(lszR3cU8;?F-N*1<+B6Ev5$G`gwvNW z3@Ll|+WBKy15>i*L9-I}g9oZDIk_T!0_S^HLEGgHT}x6ej+y0Z?Q0I#Yo0|1tyORU zGY3;o%SA$7NZ63t)tVdsMIND9O2tJX^lUazSs9hnPS521{e2ykWz-q`|C3@OL_)Fj zJs%X&5o?F6=3audeiQZ__u&pM&%VV5PE^{sQAZxYLsU*^@*`F|a!=Jc|^Kgt2OrTXJ7M~ba z_noyv@}Uv@<~e|BSFUSpEG;)F&(hms`SO_YFmrh}e8)GGgzoYgzW!l(PIEva1x~^@ zz_3!aR_n|bJFw8hCp z38Oi}(;!3JB#0Tp8CP72sWHNb$X9a;HNqQYnTcAaWZ&6C4^t3fCsaGdU^PoUimcb_ zG?fn0;hGxeW$!E*m%8*Aw3n5fln=M8#$&H!wh`a{9;54r*y z84t&p*hL214I_-3WqHpI;9`=6V^U-scUc^f8sU<3BT~!Y=w#EV*+_zN$V&%o6lnuk zQwER4_&wsKEBDL}P!?E<)2KPY1GE#z>WoK~>KtRxK$pc`6gUB{dSqn&Sg}hRrxzFJ z1p{f(jijhE13Y>tn5>)}Ge<{9`Me*bVQUGjNZ?Q>O~=JryUt>;Dl`|J76hx`Yji43 z`fo0fX5Z8gs+to&El3u!++aMJGWW1qD`?mX?T+uYLWI;P|os@8R5o zpBT)C3_L$Shqkxx0u98Mbf%TC|{~ zOw=4N;EOwns5bnY88&+RAYRKk=DrFNkq1U`={Ge3bd4L0h(nx-^g+5bm4wo@Z~Krw zqikvF5X9m0cx0Hu*|#p5t78v-)@HDYvHddmlxk5)~&^X z8EIruFl!w6$}KPdt;fZUAXL1Vh$y+JEJ`)@C;Uy>H&-DAr)YP=T~CjMd2Xb=f7lWv zl%?U2VpOyl=N7gJ5N@@DgyQn38aC}rh7`vtckUWFd;cI^nI|PO**jJq^Lg37k&k@q7Lmm_ZUSZ3k|h*p9ygHjRWrQ8B`PD-UFS|8DK)m-E#7#$MHn zqPNeNvPCS84{1XfgoYsjvGykdIPHpHo&_XAt8tt9Eja8BK`tXPcOG;eHg7XdqkueG z7^bvgvtJZnl-im0C9_4F&>=+0N&Z_zL%xofTu#hDx+Fqf?L>q58K4M(7MDmO2Q|7v zr;hTCz`7a|a~XnsB;m3vSJCxfzKJW*qFt+3jqLW>sqMWsUYp8PP`nT5Y|Ru7-lrMr zil)vh^i%@o^730HGtLM@yyzAP?8j>K_wt8+<7s@fhboU) zol`Mt$Z%S@@T$d?+E8P$PwnSU<6~nmh=_W@dF z0ex#UIEQ`%=CXEQcb@*siNTv>>D<#O@IiCO@4S|h&(+(`(}|M+O#P}>FMBEcFRms5 zCku>nE%CtwW68FY``JR7s1K!V3So!oBM*Z|n<}9fqxfOteQk0Swh?+jiTPKZ2$bPB%7^qRNsbp(Y zfdV}%9_>)*l#?tm1tqCX9#_72cwlv)!t5hcGGApZ1bHq`*P1eX3y-_Iy$O9qmTxSu zuH*LtwR3}Lm=EL)e%qB3g5-Hjd`K=w$Rgzj-IHw&0PI8_1*%N2 zEf^5u1v1aT0ffz=A*z3+gOnKQe13(v zl!e60K%tmJhXX5?o6~GYRS=*8l#@#N#J;8>aEr~@>ASlbuxoDfyZv#cu79)_=-eKx z`NRZ~6A}{YPxR~(KR(n42THFt1w~FQ+VJ_};c+*SfM*+N{0trRXdRqE`aVzEd(B5a zQdAuk6vmWbUbXzENJI1fwvfGr6{*wjfnH5+=kw!9q#>>rlP<0?HSB;0svfELVKjTt z`aiP#ia<^s2-nQ^X9I zHKlZeGk}eX6O^~1CrygO4v6?=`y{p`)}evefUd_yM6Gfkh=gIT;cSpsn9i?(2m)$8 z7nGtS9Bi4^?3pyU8irFDtUk~^gJ@NsnK>y#Mb++z5ekaYH=g)^echud8%U2ez~5;6 zYTMCdqvJSiA-RVl`=GT8@7SMO_TN+|doM&u22m-mq+_H>EpWCLgV?rN65>KcK3`_t z-~V}*mzQ_A9U}z)s|c3Ek_>ok{(d|Zi3rZTB6r`ncP+iVyl(z!RZEu4@6iMu&^(`g z|2X^Ns13*V8-Mo-Ho5$a64Rw*-UOUb+O_h_G(6bHp<=%OiR~l3`44La{^JDOYwH2H z^?M^w)-9Yt2x@gs37DV_TvoZq+FlYR{0+>2INZ5vH>dT1_PG9*bwqX}bI9C8h6cso zix#3;N?_i=HI^`v%%|%#C$$huC24#Rlc%x1=E8visFUwHyOd+MImsAF)lj{MK2L@K zNGk7jgPiOV`4q@Cnd>N%en6TM0`}NdSG?W~>hEiF5pZw;(Cub8c;F;K)F_2d~JmBZ+WViV;vxZ@oM=!uRoLl&tpLgxs zkQ9(rN{x=!p?sb4UR-2^<@1i^%gr%Ma5m}NASr+`yZ3%U{%u&~Z8%iZr62$BV!b)s z67?>84f$Ob56u03!)3YhQRP*TmkXhd1+%V{3FiL+iL2@yvsOnHg04z&zN8Ey%xO;k zOEqTf-$0Zp8K~71y<(s$f?9_L0)Zo71CE=57M3}IlkCDCz}$+F1-U<{v!`}XwJkvj z0CoabQ`}78Fddg&EUGc9pcR+YoAY({R6R(l~t#Pvr9W`2Bt`*jtO7Mw~ zFhbmh|7yA~>XnKb#s&gwC3>|u&CShU!h&>ccgy56!Q1Cvm3RnqFged}*)DxH_~UQD zzs2t-&$owQ3>XXMBEfotAZrhp;^ibjz7w?e_4WNmYSXG`>FfJ4uz9UKwx%%fsyQH6 zX~b~7kdVsCeb|5n7AgpNoc>Gj!Q?Rfu#53g?#Y}IQE7u*doG;d>Ui>FhQUAVI96V@ z!ZgiHQ-YF~frn`uI-5wo&J-DG6PIWX#gy)>mPkSk7Xuwzc|mEF8?H`vPO%v7loS>? z!5OJl5P(}uw3ZfTAEpet(<|Cdl#bQITznq4#*4U=bxD5gpBa>^m44zEZ8elaOhP$f z##Sn_N^CA@R4suF7alQMKa@Vige6+!bHh;X#}~UKG7psAwi6}Bj?FX!-IToV3wb$n z8YtdBCX-EcFth_T1-QAn@va|OZEB;`ph@8&gI_myPyT-Y6s+(6;IBKs%bpDY?Y`vi znp#h{ z_7LE?ps*)EAe%5q&IL)a`bSMyW+5E>{LnxENvyM^bBFPUdtI;^El5-r5$QC+(!Ri~vY7zjY{tX1 zfOO9=vyJ569vVLRWXs1mk0CaaMP+PgNGj>k;KX}tU`>(Apr;)NH@yX|sz8VwvF2R* zFi)Ifbe@G{1r_M8u40&CNTZR~0trXaI5;?2Sy^}^jESToN<`xUf7HJpw5}R`*#mq0 z@6SsW!q5ZC?@P+a8K=I#p8B4Iggq8NvQsh)8p-!J zkz>$O>;_dB#HtwK>ph8DaF!5YLaC09Gc3w3Lz0%hW!MZ>eJBCp%?EW%eqeK+JH84t zwga!3`2&tOg?9~I;824Uh1ikQX|{sEqVw-BkkX|gA}^E5|MT^;0Vf<_><`BMOdmI> zpGjnJ&ojGkiaUW%|AFJ@_ubDx2;s-Z#9MQbNArLy_wOZO2@5v(1{PLnAk0eHa_ZI6IjqOGHlVn*)r}x=1?W4$83rab~t)+8)*yjNyr9$qg;-{lGakL2?8V zqVo=Y81uB=9;l>O;viv@@+S|%!S-{?5+MEwrB)y>`I^!VXGzOvsFy7C&L1O==Q-&f z%HSTS4!Ld0R!1=x=L_#e>z+*8%x1M4syB~Tu5x42Xt_(d%cYQpb>101q@`2SN?iO? z#Z8Rf96H0V-p~<|dWv8+RIM;Zeb8@w1vmJqZF&Jz`>%tunvbvUdu-|UBEc77xCPEY zWr3nWbo7Qp-hh2Xn%W3TNu*Kja7kwnZU7%qtHp{DVdzcQY{S+bNxT9 z=wh{SooA~n2nKYjdV2yO4R}rxdB><*3=d54y}=<6nG5uCSZ$up*{C^UfoabVT?bvu z&J9lo-z=P?slvDcR5{Mgl*y)iX&SnDm6WZ4tdgW+h|WPc*(zfGWaWAI#tp%e{P@KL z$P!uGy%?x$+M@J>m9B8MFh91EJwlvZVKEv07rc?20k$+K@p6++DmJOPSd;Rq8q;3u z-y>o1b-*crHvvjBmtt9ACuenbjbNwFQNu{czKQ3+Fvyzrv)lXU&!1a71v{>JqBCZ@ z4BoD;uF|NH?`A9rYj_;2aQtBEv|wxg@gTY!`DcCCNOma(t(*A`F?hvD?3AMc0DTrh zjQdGy9vpEly0)EofI~~gw67c)Uc}JxX&@kOi+KzdD5j~XF}cKG(r&y;L3)k95r|vk zQAmJi2QAL=r&8&d3>xgRSVJEYfRoS#=2-;I%Lrp|S^{ByvZ-`~I)!JzlmX1G?FP@W zhZHG>OeKnF4g6WQ=Q`7^Wll=bRsch-EW}ih)E*G%J4p>miMu^H0Q8{%nCr}5+UKID z<)c88i}G7+Sa#A)w7x4Otfu$=fR9$U0lE0k9A8573RPwiaTw%jYLNV8jWR_Py@&!! zwL{JL;<5wap!(zs_0(J^FE78<;TjJ1^+Hr@k73(bF8{Y8iuG+&2=#Sq+7HsHwmUpR zRlpT4FRs-K5QKcFCDS^{#m!xHb5^Rta>BbXrn1Vfm2Xw`9gw;tPY_Dxf)s;sw1IYx zuqdg>fRh$5d2JS*rQoj1#Z^T|p$F2s|ADgg74fq{Y3O;*jaSdw7(een(O6qe6Ox4(lG zM;nX6>F&e>O;i^4qKWDVSSF*KG(&t3VN58g*OAX_4k=KuDr9t(b`FR>M)R8nRsOWUCaC(wZrMKV&?(tEf|D3I ze1}ZJ&C}A`hdRt1*PK`cG@1RYGKHanB?+#nDQUHGaEJrzMWnv3hbL&`o$jYu1qH%$ zlg+u58v^nqRK|<4M0zpO^D!=nk|}F*3B<7l)Uu8o;%Tw~XsC3R8%|eh-SaLDM612S zxme&s@tF4jb`EQ*#Rl_}2j-bJ-947gAjgs%0aAFW7I^*$v1KFV8ey!Scr0{(sKf;_ zju5o;Q=Fg{H3?DloaCB}-u~&#cwO-grTCZ8>ewQb`_*%8G+^+fnVEM~7U?auP{csg z)${r!m^KyP8VuKJc3(Lm4Rv6ZIZB|5xwN|`jwB>9yfu+cA;}>+)vsrD4Y>HprJfYV zY2IA6NYJ$5R;Z`)*&2mRpjkt`**bYzbz%>l-DWkBa_Bp$Dv)d$H0Yl@eQ$Wh!&&aw z&Dk-nfdr9L?FAYCGk#w-cD+3q4fUPz?tk5;^xRw3LoFvJm5M$z4d=Ud;AL!cc&i~3 zhuW>k$C@;o)G2f!j$dgY=mON0>!|=`#W^5Qc}kgo30Al;Sq2ksSf@ zih~-GsxFSQx|C_0;*9G8p!FPzxLP+Orwvvq2CSbN*q1WwI5s#dOff6d1azE|v(z&O zOD_JJYRq4IOI?Y;mdk|8RI-!bDIdGAKn85-^*-U}6KP`^>Mvw64Ax?cw#CG8P|a(J5>h3DU_l?2bgSf_s36i9-j2^9 zau^w2+!sF$K^e$)e)!!!HV}*6QIlAWypNRXBNCG%I0@HWD zpg7;02j-+6VGB4Mh}nj+yO@uHFsa9|jt5*SOyO~^sCvhs#A(nYVoY3zso|%g#qsoK z_mmtmBa=H$Vhi@iV#5bmHJ1T(gYFW_m150^R3f99bjxzdgwt*O;wC@RUHqd&jFnEQ^Ca@JW0r4m`ttLUL zceTE}>-K?tto9u)NL6LMuQ2wtFQ>sD*1rmbx<^HxMH8&vANMW+z3->8Xq!RSgUAe^ zgpk)iPtrk0($6C_!mi&bH%u{I`qRK%*RSs)vtWg;sMg;}$?L-PVD`&ylWO%*%hydI*$VCWu~>zH9!k4;=`Nc6hDcsap2=FCmuwO+ zw&v7QJTW(PKI&IHd%f*B=%HvDSL@W#%M2<13wjm=QM@#!R8gHX$u{T*3tlON>z>Nm zHMtv?Tzt`k1=$+F&h2`s)i$2`w%UK`YDpF#sfv1VxoAyx)qozM0c;(UY$>OLHYfD% z917LtTFD;X3cm`SExqLCqdc$ zdo#49v8A-C$M>O~0I*n&g#j2Wn&}mo;5+1UR|O=AoLUs?_k9F{f)T220X?XWEYavY7M42Vb0=6x#_o!UFBMU-Db}j9_=g7jr^RYkpT(aq!UoL zeG)yN53}sP&-SlleP?l+CT~P}>HZSXfFH7=V($nm+-212B0c)iT+}vBl+M$b zw^T%mP8I_!ZBdxCZ7N@tHj#4b`e}S|VDkWMjh~!N1U~mES1kG-c{-E8E+zhL+q@Ll zUadP%0m3`?y-=&kG56gd>~`uyf2k(k91>8SI+I|Iue%s^*rMB8TeNAIS(iyCra)=L3DO)LE5%9 zUQ~GE&GzTHG}WLm=Z;>4_^q)Lp>dT*>{hnwduPV*1+7bEEr&AMz@KAY*EQ=2Vj<>- z+(;by&9Z|J+%1%{&4@JAh3k#DJ$92`_waisu@d&f{ei@SxJH7GS#?95<}PM;D$4at zXkthTJKnidkf;UzPYurxR+LX418eE&>926FqOUx0ylzlIr}^H`2e|#rmGv$P5=<-(A>~H--M~bbuExm3)Y#QE^}72RWcT3uvAs#P!jA%uV)r$E zy{%WOH=LMMxqz_Xh&1D(W4vNVqZARW@R1@rj2mC*m(tAC5xbGn{4CMOCiREwzGa3R{U3#3V z!G(*R-B4x799+XXgAENuzSl(f$v^FuEq!yMh#AUCzw z$X@T#yp=z{8Oz&VRZj&uol2%S6aJbZP*Ow`e40dC3_9I$d_+`AiRkverifT$FZF01 ze0@Zh0~a0$r1eZKf+(p>%CPz$?fb_v`|hbq)7}ml)kOTdwY9bH-x;2E zMgbcPr@y>~VB$<@4)d86tAP!)ysTMg*gbZ1bY!wM-|qZ;)TwErZNWEu>!|NW>5)Zy46HP1N@-s4Qklv~n_0nfb=pigz1qNNKFfKh z_{p^4bsUz<&Q-2!3q5kyyLoAVAKPb9iA$p6ICR_PDfE|k_^lWj^5+UyC?b!RNiUNt z@vx6BhdyQtB_W5Q?S^?yCBqctmQl%SnrVAF-C(1ucnNZRhH0%M1G|zkJ$qz4RlS^& zqMcLykb>c1cSe$PtzPkuOwM@jyu5~pKi!~eZ$GA`zQIheM{DEKg?rvneC@>gLRhcl zzS%`Q*3jZ4zZeUdSI`o}?^uW}){Jawc0<(RvbCeAuv)mxwWW*QB9^Pu z00x{h63JxB-;w#ewW=6739&^IExyyac1TAZ+qBW9jzLbx-xq-`K zhPRQct5kE&s;pX%!RtTcppq?~8{}ovv3N6)PNgnAbg4Pv7J~-jP|xA)u8^@l&#@YJ zF*9hGZ?*6`PPnuX(9#AlgP-zY#!y|eMytgZifNgx3`bTi;S_X7H_G-7QWuAM-~xM# zSasa$CRwCrlAd^(O=<+#vH!zZlB1L4Dbt~KbH~@p4a?3kEi)CztD|u#v%0E_f58d( zvN2{Mx+gQ!-R`(1YQoEANgKZ@QT+AqyZ%#->#tunX*dd6Z(wKFlR+NBub&}m zph!)aWlU>plh#UJVe%Ux{g=ib)8Jpf1}?D`C*^th%%!E4r?dtg~5fx9THp76srUtg=z4#0{x_#<#FS?8e6tKclWdQ^R2e+V}# zS4Jax*9d-1E7jt19&0PY-;Jj6qby?5Qe*@Dw6?!B#nBv-LtWx|2vbheN~6&ZmS6i% zw9fsxy1{u|ZK;kUAET2DAV*JBl~E=qd00+;4J{g-sQz@M)5P`*f3jvN4NX%4Vzwfy zq88_+(N^<3I<4jNM{(hA!g4!Mn01FAZX7yX5ua+59Ap<=3#E?zy8g>qP9u3JcTMXx z-Jz02s+O%4cxOtn1h`pArJlN&i8U;w`9Re%;tAp^Dd^y~0MWLLLQ@w(otENd7FjiG zS%33Ab#~B}#yKzRiqtfgrozHd+qG)OB?p@ucN@Q|6E=@~ST847xB-{O2hek};qAqC zZHDC*cMN|n6>W!B@+)W|ZI1+b@!%kj(?EM`)t=3->rS{0$H412$?zv{!wj6^!%@z3 z;)`Ywx=}qpY7V#851rY{Z_tzjtgWFd%Bd=>fgFd*eZ_lZnrOq(lvI0mtyAAv`gqCG zNI&*TPp&^du1UYKy`3TsIv)g6@tcnYgCCqKKU& zXH_fL9=Au!ryaOxv>&H?mpt%T7!F2M=xD)#Dzt*GGpDl3s_d*`ElJ~BteuHVWMOJt zkpR=d&x)flpj}+Nl?-K_3NnR#3nlbldFTO^TXE_tl6I~4oqPCAk0cAN+8v|Td^@PY z73iiFj)ctJTz|uJsklYBEkfDdIHBWc03GK?;Vlab`ga1qvTAs zS5lBpuGelcw{Xg|D^ra-OvNe}Amw~x?SJYIlKVTFl?sO0zdW*VdQ;uj;;L9D+-9B8#K_Mv@Y%@fTVcsDQN z($kc$Xhb^C9`BevD%vTKQ?V$v6pqlhQX&_Ao}tjb|G|Dys<(=o(EAEuAfM|jUJz#{yQKAmLvaD zm>is)`{}aYN#DYQAp%|pFJ`P#Z++g%cSR76Up7$h?Z+W|^V#T}4X#_@d2(ku+)sb; zApL$>D3}FSvl7!yg%;`l;{DToB)@%WCuH7*7K1NuH{^{JB~M@cn~*;}E{ZaWH~o)S zdn!7$rB#eo3f7u~)rf9#M(X@`Xb*by3h}5EfvZs*ue=E$lj#f7Gk07 zHOpcLOO}Fv@fBJ%MqNbw z5S#m*mc5m5P~VzjVQJVAoolH=T9;Kj^!9l(HeuQ2I`vv2&>cLn|0^WA+xYTDPB>A? zuyG}+s3ZZj`s5y`No%R=Ob3PAg&cW4FeuB+%+HXBTal4pM}5PX+PNP<+xGi1p_qrg z+Mq>ej&~bE>)w-LZe%g~U%9^ho z6?xNeO@s|AU`-cOWyi*zTOj8FJN^SD9_GavGug7z`GEbUjfwXAzom-*)MoJ6FCrqM z8|BY45h^6{{}hMD?(XSc|EK?Q-P;HRGhQNgB9xbTaGi&z{C3hNTp*(N$tdX)?**DbV`x#PmyUWzqC&%W z>)I2*U?wL*)~sMyHRESjp_L2g=a9!m+UJ3ag?-jCe@0pH!Lld3LeggYY$1$ER8Di! znY!2j%u}{1exFKIMqP^qpG=?3C$%&+7iDp6OrP3mM8>+b?0-zMSJqhQ+RGiDn?dP? zhZC)8+wHw(G>5Jgi~6hpvGcDe*hU&vNKJgL41u0w^f2l=mBTwRa~S%>)OLSJ{!c9c z^hVb zhrAQx$^c`8lw{+2C~F42AC7KX`YRtnowI$hL~JB0yF`qy!&hG`bvAvSDam;(!V+G= zRK=qy(7~X%cz4l{Ol-Yj3zYX5(~5JKQ<>+3_%@20!DaUL`aOLo;t8sLwQ$z9nk6oM z!>D|&89n*tbyr&KzCe|atyAIgDDZLNcFlDHbGWcQ{~@z?@G-1+>Hb5AH)HqS=YwxO zZ+FwmH|^k47kg9b3~6#yV~78p0xr|5b$huxc9zF;irfEY{6?|+6uA53>D<2J=I1xH zzfZooxe1o~dXF=mKgIsOYzo)4-CL*Bsys&S_GZbLJGK^>#u;w7bDhWMv9>)IWKn8W zfp#S3YYB@LO0bZ|Gt3}n5G3f3Dlq!`)VORo%n3RpDrL^wHZeGS*+8uGoKrB;y0M{V zJrJ{@A8Tqire)t&4qT#vLjh4CdYX2pM4ckE*VUQ| zJ9?MWJ0>RRb{Oz}`eN1ZHtJF044GGNci1X>tx>O0UAzgss=uJas{VuMFiP0e-@|t2 zb64g|{%d&^`kxuI7i@Drv6kKoqA0fA`mxTH48CfHlD{xgyM$dRzQI6R>9q(?oej3I z(%_0?zOL=L1&T>9v*H=5wS-6S-}fI1d)N3wAu}93Vsu?YR}eICnRCoLSPJ71e#$q@f%`%Y&qSU%efH(8!33dQ7K=Ev`k- zv-12BRU6~tg#&DNLv{W42)ubdOT3COf2`L zjqNY@B~ttkqkCItdg2gWotnC`2%JqQ-DM`dru-g$;l`4I0*Hi16Fn+)(=^_dC;s$T zgZ-E>w63k}Oi%(sOn%c>?c%uBq3V6u z?GLH5HY4P;FsNq#doRcmc>AYLges*2XW*y7akE3nwl;aDEbHeMV`@ISX+ctCw)QE` zyV%VVItRXlZ$v$dPCRbPBl!m0)|Dlo9A0eWPQqRv_?ER;V)xpAjICoSnccx_go^Ka zqEHJdO^!^y@RLTszf>}m1zmnXPb~{&I{L__=lIi!5pQoQv``ZthF{J{#5YRS?}N6N zU|Ec|vGGa5flNlNUbnv94#Kh^e|d$3X^0+00th$gZD3h3Wb-6pI;7Q_N`x$eTk2(#4q`Y1me?U5 z?qYNG>PF4RoQjCXjT=7PEHF2%2qon|NW!=H!yf32ozKX6J@l>`uzc~&SV)-g?Mymo z2EI!YDZZ=zRsSZtpRn1#ekVF3%*8les@fBwp_XQX<+v3Lt!P*g2JJ6WX>>E=WP_uP ztu#2&NmwfxUCh#0K!b9dBZb^3@F@vC3c^-}{xN@)QJcA(4W@A<=uc1?jf{vLZ?_7+ z3dotmE}qqWMvI5pW>2jKLhsH;(Ibd|zh>z=V$M7g9&Sb6nZc2XcTU=?%k)73L>zvZ zD^`x`iTOsou@!X!&uYcXvQ%qhz@KQPE*g zYn{8I0sZZS%f5Ct)!ByXGG`o%mDp64gJ<7i!;&Vt6~d?#!6Xcg-HgREmU*60ms_rZ zo?zI?7FlpJP3|?gT=n&u2nX}WPz0F$uqz0Wk8EdFHZojDtQclJ9}?id*u?yigXL{xSdBIkQVwDvoT5sY9xZu{7PW2xZndf!qbt65*>b^6pe2FlXsV12ym}*?*Mju3D z?8k|ckWBn}pBUzK3(RJHwDQylt`$$%JYc4i-S`^|-;^1V)Kg=M=rQMP%2Z-GA`90# zDQd?2;wiyv$Rg|1xld@=LWSp=UxwTpeU0~e56BBM_o9m+-CFXN{LZE}>OjtxN{u+I z+qfb!;&m8`K;d<(W#4z|Mtm{Tp8bnEkkfpjNKq`X19|U?N4x4n0nbvp(&*- z%{5bWt|r8-jy*XBapPAJAun^}`O0(!;|b}GdgP=%1A-CX*cu@*Ua#Qv>Zx1Ql0R(y-TE(_}*Vp$}p4?H$2u6lc`I|Rikz{M? ziiSCJEuuO2epP$^Hm)erm3K}F*|xKH&dPkQmnLij`HPhrpXp2F<}mz?oB|=?Zu%eL zDS*&4Hp}?yLf;yg3pC9#L6f2~Efwloc4A|`dDFk_{~XR=IIjJU9{L|~M0WY2ZyNT1 zllRPKJzbm%gRP&Eb>{HOWQENMxjEd{P_3MAI=}wCm(wY|1@`CCR}GQqMl5*TbG3f3 z%gMw-2w$)q3G2)L%bK6)`25lJ!R>?K!=F0G?Z5G1TWax$=XB*Lw|ZUIdhLMAmC%el zm`xh@B#tQ&Cp)^389z$D&EuGi`}+BqfUMJPJr}e%nD2Q)$~K0`B8C&5iQ$JIaaK7d zE9{L^#OX~|1~L>BMR7a|z^iLO@5(On#xLc&KOyX*Oj7t9!` zD}mh{_AnUcO^wl6PU1q!b{{!K-AgFpXl^qHDz}$_>H?LHW7_3tWABvCz)jh!uzZZB zTMQf6A($C9563{#N9R3}p7be1wqxqq+-p%!D}Ksk50jN~_FmN)=`iLSuC^Y3_wRQb z1>9oM@T+j5xx~E_v8k~Wm>u&5Kjx*f;;THfuu2VypO>$GPzUE>@!p2?R%xMMC;rq8 zZh$cJtO&VBsKsc8rPzw&Frc#bT}C_aabZ}D&Pz&x$CmA&wmR5hs#BP>M0(}qSg~AC zjWJ8D5_A&<^8Mv5D>29_{dU3@9GB>RyiZ>@K?4EzUXa*!f?y|ya0^;Bkotr+kOvdj^ z;Z;w+s?PtZ+i?M@)F}TXRkrrCo~Ea}{HgCiGh^Def!q!X#PKo9OZ%s%u7|`uCM@u3 zSV$NFB7A9bx>Sz-jyZ@P=PPVkdDQIn@XR;Uh7UgaGN7dOs-_(MYC%3cb`SM&UH{jL z=}`%fC(~u-)2*4OHNzXz7H(Rr7j>R8I0g@q*#Oz^00Sd9BG+&}x)^B!_ji%bjCd)H!e>M*wt`v@u>SJ36KTKE-ygPaT5qh^ zhc?McNpl)3eSIG?457zt+@YwzQKV47^sz>w1w)oz%dS{E6Xx#{&+P7^! zZ)<_E9wJ>{z1znzOrp-Cr)uIwA%Ai3I{q&i?Kn(&ja;yv z*_5cm)xsb_&Jt^ofWHAJx)^c{NKcf}=nuttrlX^p1k z*E4vxp0u}QV^lV9)J8sN#h_=PP)hmv;wd3lWPC~2eMw~aps~P|#V67xjVoOz@+}`f zuNUkKXbdU_A{m%LHbakU?(>s*_ycR{i;$Z`Op=rRrDMeSdq3MbZB(uTCe5zwfc}(V zAqZbOu@vN+&MRfzk*9@F<642*t{C>G93y z*UVc?bj-B-BH1~>w$L$H6bIyGH-NrwdeRYUlV(8#^bz9*8Gj~iJT_qPnvJSehTOqk zI}wKYodr|uOgt}%z}cR+*f+g` zg&V58ba0j@6!C~6!^T-ueic-!{8xp{V_V78oFS!Bj-->!@_a^K70DRe=}I(1-ix;- z_@b{4*Aq)J`|6L>?`xcHEt4HNf6kg&o6s<~Z^02mBxIPW1oP!YspuL+i1D1oi8)zkKzmy_yB|6ZVFs33c zTWBINtPq#ze_`Z;dWP81Sc07uX}%Lmz1Aa>Xb2I!9bzXWf*m(!-l{cfac5hd>%)V? z5j<8Qq^B9&q5I4kK~hB9>0pF~Uzpl*xg}1?fy-2ZzPA=}64YSS?>~1osHxz1r>#4j z>AP+j#DlLZEMBz9Z!aLd#%VS4F7$6ciX<&S>!zIu(Oro(6beg9lJW9odYnkNTUk|p zClv6PGnu-4H|76kUTT6*(pt>wE^Ime_Hyn3(TDPJzjcJ;Wn$5zL%I<=r}&w{p^>2ZH+1kxHRB4gN}5#cb8Q^ z;Q`bm!MelK_B;B~#_TDW*NSvfI(3=+MWy%>6rT)|`N*V^N4e}EE9xo}mWVHG3COlTPI?jofvU%8cD( zLLA8T8OlzZn>yo8NQFU#Gx4_EyzEBblk6I+E1JM#(47rKbn3!CzaW;J@TG(AVAmBm<4NFNtbyKZ6X$~c9$ISU-Cocck_d+ujfl?#izltku@=}mv5H)Y@_|kQ z6E?i&c_EVLzgI&9wfoBqpxVU9=WXa(y&5@<@ z)x3dy9P69+b7QS~X|dAZ%gt3y+^?WFpl@6G>YYyrGl4cts5I8?VGkfxQpWD?IKR-p1_4zx zSGLUDDPInAra{5gAkoC1H>ASyI(p6}S`yRv{SdyoS*A!mUk2Hqev*B{ zJYuxz_*B`|kM(kf)pdMyo%{~VPWdsh*sx?7xg5tIPRwY>Q=jxh6_%t zL&TXQKLSV##BN}KSQs^w_w-L3Fl-TGatCSZdXze&=wki}BnQ#nUOC8@jAvc|r~JYp%l|WYKr?Ub%1B7`FXTV_@WIBd#nkTO%RS{0)uzs#wQ}8rB^? z6BQlxId*<)%e4dN8(Tt3aKz@0=m4qEZM55i+*=9KeJ|M7z^eKM$b&RG;qf~k;C%0A zYLdWho1023G^S}fy*w7E&|USVum_2BnG6Hs zbCS1KQ~mK1bc7p!{@dy}mbn(5t28n!bBP14{V2@=gq(5uws=akfp?LzGSqf6=@#UB z1*JcA$J|VS<}xsy(6eY9O{4h6@0OxSbV>g_EtJ3TKeLoCKH+o>y)xr;028EnC>?^eMW?h8iCo>u#q){{KuYS3jl-mRG*;kIGa-UFO7exl= zy`=>!t&_zJj$~pxH=S~Li3cHDWhbgTD@hF>Hk8#+_ zyS>x$10SS_cJsC_isT?V|J9dW*X|8ug_!0nYP&;hF>-|vupB6)cbn!tv`atrm4XN5 z9S-BT2GpinPaM3&6O=ik=%sY}*GRe@u}8Ld(?6VXSFJ~V$MLcpq>|}-Dkh7nW$+D0 zxg0Vp+%%3Z)Sd*PMXG)-d$FOXw2S^IgOLTPn?)n39mS!8hqNR`n|?n3vm5A5pAlb! zB*kzOhdaR4+s69tDx4vvOpGJG=)lb~ebg1I(SQ7oVK=d8buu~9)}~%t^K_6QKjMqa z5m~BW>TQ>D;PB|#^aK3yEW9kPDp<8toOH1LVua}*Ps5<~ik6alK=Aea%waOntWSS> z@;CI_eEW5i;On>b$FcRDg|Ej8orev#9j%=}JEv61KT{NZE89U;~s8L>ww9>cfx1inplGGeK3U0Sh^ng`2Yb~?4c5b8RnS+8|oGhw7ggKB-w#@xenzazv9p@VX~{AvKTkVT6Z z|AOGVT?;VR3)GDxlllvC!T2?YYlQVrJEth`EOTdhC3e3b0t3{WBZq8be zM__mfOjyr@b)#5s8rm3aUPA|O73tz z>@_z5cmNy>@(*rf?_9n|@aulIDiQ=CbvR~am{=+0%3rF(&7U$oUW}(R17-3O=R*92qn2SW z&!IRXe~Yu_uSqb&+n@g)>34m7Lg~DXk+W0CRnAeSXD39zYx(iO|8E}u@>AF4C#B;z z(NzZP=Pwu7W3OB?e~ei^#&n(Y{u{*W+D`92o0D82d3{%z{5~$sK_IsG5`gH&_OBc1 zz1|R57LGqK8=wKFrya}6JeQYy zFvjpsSgHe62I)FdyXcuV&$wtCVRL`Z>@m+)Vd7u;8b9g4S5YaQDy+y-7b;)V?9G;|IpwB@`K{_Z zt{4c--qk75gvCt@3>2uK;mJ|v)cl8g%3|Gl{iO?SW|#VW=^4^W5btT|j^z1D1Hxr! zU=;_E3fQcjo;|(y1}aBLG}15SS8W(nl_r3nAdPdRu^CLDq6kI(kCq`*-+Y!4%Ct;= zsVph6+tU2!CeEzLwb!=_?;uWfgeE+QuJ!c1kKKS!*S&u=IL)UYJ4KALxQ^MVs4%N9 znl?;VqHgV4T8*Rp4f!ZR7?!Zn*2L3c)s0m}11UMVi-*Sm;03zFU}%#%O3^ahEdKht z_-^UPlg`%*x`2z)BGOI||Aof%fa|JfydO`0tYi{H{&%_nv1b2%aZq^Rc)6*?^OU^n zE~*Pz`uhXc(zyluz=z^MUJz%VqM=8Z4yXCeqe>k9m>| zCmYNARgr4a?B2utoidgd>mpa@XnWWYL8GRoK03y?;vp`MSA4o0=~g=}rkr&Z>L2w7 zm61nQeCd7QSsE&2S9D=@ozyP&9&ZZtSuU8?XEJYH8np^*hh5=z52hTeIx}~y#{I&7 zY?;aegG}p=Hrl4;KOI{HZV63DU+;f%DPSYOt6EF?eU0tU4&AvL;GatrYjjh@Lu@gg z2nRkIT)ozDgk6XUEL1y2_@c<`i3a`E>mkl}b}T05$d8$A;qSfpD{W)TIsYO`tE$k* z6Oc1bPb~kx7hqt>`p!0Kex%Kj6$3{!@2iFiOaY&-MkqH;Mst9>)DJ1X&PAi>{0NpW~%p8E0_ z3&2qjBoae*Y(|+GInwv^*RF_^xEN^^Hlm(A1c%ec`VXckB2F50HNlMUB_YgpJ1qxA zO01l=7Gt{#ardHI^p)7IV^cPxZnj=0J^(;ews1%Jte=`QYZq83DbT%};1J#ZTn$%b z{ns8zX#2h{5O(O^EF)$)6RS<_8zUp0U91l6&B66oopGqN6PIiA!&AEO|D`CqT3rQu>H4( zptIZc^g;OfdLP}Wg>8OK_Gj%}O>_}+3Oe*i&!ygPR@oQyO&c9_Da{(f*Jz!!c|^!Y zk24S%m7rcx&XPXY=xM7R!ra@m%RP!8Z;sTaehaoEmWm4geW1@eYxfD0HdXe_)Y_Vt zn;$1ApOmd{9a6{!H~665J3zIg zzlPd;P^`{0Ui&via$azh0R$SPWB68OGxzp!`${FP&7dgOA~`bBuI+V{atTSLgn{eY z1finpntRKT4qKt7Gkq#p7w#*7g%Q)b*YJrWCDF3-c=JKkotO%GF%y;ilfGDa(oGkq zT@5*e1d|H>ezyX5UqoCfcnKCgDcgWqA@?Y}(!2+)(5 zza9FUWYzbI!wfiegwH=sToR%q`Q5{fXpw`ONKvjTD_>~IgctHF3`Qy+SQ9s1T6_IDtyMsjB; zr%#b;B9d-H0sWmS4fzi%;E6U>x4|OJdcT{+ED)A)i#sjP!f%`k zPNDKA+@l9sJ!Ou?idCZ8m>$Ns`1EGlwUPJANOt_*;KN_%Ug7^2LJa21Wb05UBXAXc zmsDTpG}j(tOYjf$2#zq^A%p;JU%S=*omv|fQL&bS=4;AaebDb5b1=FBPrUyKu7JraOsSt_r8kt;txKdX7E1>bnu7$@u)i8M=c zmXS53RyUiL+1b{b8*J2qr+S9BiudkVdV~xW_1ji`DrdTkT9r!bJ{umn|0%l2W3c=4 zRrv3?u1i4g67)I?{0E10HhUtejEX{g-ZFZO^1JJ$*h~F-n7ZZx-koXA_IY*U(dGx( z#e5!JAr_h0*?_O`pF?cjg?h57R`DS$;P$r$enjq!wxTY3nAYtMf`L&Dm?<*pt6~5% zxTnp7Qjs|+l=Rkf0WW_zsC#Tu5oJ(Aust4W$X4#vMM0`iwxqMWT#~sa z9s^ugu1Cw#)jRVi_5Lcr^-9Y0k)XDOg)LC@g(QU3^2N3n2-Z;#2PL5?WC~P(K^f&e zp~Q9(I7BD`o+7y_?u$!KGOARV)sxEQ0xQjAC~Qj}P6pDZV#%_Y{lSI!jWeMBK=46` zVK?Iqg(cVyvyFsR#S{hgg2=5bQp+(MD0kDzelTB7$`TuA+wlvIrH|+6)PeMStwGz? zQ69`HCF!A<&UuS}dhYKLCWa^Hm6;)F;m{>PE34J*pCJeg4Qb`a+SBGrz^Pv$Rpy$ZtEQeGwLoLptINdL+WWQXDa@V)-ZW zMF<-*w$`(x4r@`9s5}dP6(uO2WmNS8lCdSjOmmzJ6q;xGwPX_`&bLKoH^fd;?FI;Y z0lQ$@#Qr~&y_cTHCfnEZsK>dkzgqX6H*GPf}Xf*WF;;;#s>}t|@z% zwe2pA)TdDySD8LtaZ{}E7^i@$b-8mPFU#eA))a)KmR$Y7e6CzPKuc^v28wxOX0`Mx z2XFaeRew7RjscgWud>PAZB=)D5%PO6TOrI@m(^O7yph+9MJh9W5#T0y%sOvL$F!LW zen89GhO?G$x%}KrSr99;7`uUnbKoF>=pd@t4r?*aq+>pVg}PI6g|7lFuq`vM0r@ig z*gaK({=Q?0vwm+c>pgS(QSF2v2qnE=mBEWNyKqztJ`ar(69C zjq!d_nyY34aE4TKH)TSStwjI>c3d>L+L2~Cmb;8IZpg6tY-bVP(VIp;;%!;{A1O)D ze?0J@Dncc$L~-Gx!_mYdqU$ob$~+mh-j_!a*jZ0o=zipO8V~~}0i4G_!zV+$uH8Ms z%a+mEQP=a4-$l|7KRp@i0GxR>tZo1#d-dS|_0}^W)wAQspdGLo0xZa)ft9AZ=d-zJ zrOwPLdq#QuB?s|b(&Zy zxC##~{g$>U<5+2C^4q^DJL~9;ShPr*b+KQjiFO@%+b=jue+OMSDU2gA@1#19+DB|r zsw+#)3sA9>FAAH<2G7I3$n>Z0v1iK`Kdyb2Z97!0LZ(C?K-htyl`~<=m@?qu#Uh$m z^9kUJ^*+Gf$!#c)?zz^F@O2^&xxjZ>u+~%jP|P5Yza=jC0*vi+@XY#=H~FC%x27W! zV_eJRDyGV^WqqqGyACpKmUxhUPo9?IKhd)X3?NVJ7yuAp~9b4+ne?D zUAKR81TrBw>ddjzQm{mpC9~J(v&@;rdZr(h&(di{gig98zD@snXNnr=yw41qnhaO^ zh3)3;-8{v+$TS}=$%EKYiQ?Q`$Z-exwLB3dlcH5swr%A$?Do*8RGWjg^4K(-%1*yQ$xt-)J+V;SqprZ)pA)F>1Wi-SN|e_T6#TI0<^zofGfntLTk9P`gY z{PQ0m?z|?>|4)R$ZRO+sDd4reUkHdlH!>VsJc+#F8pA6V6PudQDYYz!Fl7g(3Ni;% zsu51ebnvmkBm8eYza6Edvrzwaw3nzu#>ff_8RR2&I4y*SB?O`2s*0!oR6f&VuOZ-q8WN{TBp~}SKyjZz7VM5T39Vd#z2XLI>`@%Yyr=2+Q%&%v6 z!bEdgeBX0r=*|DatZY3i?=jY!o^^6cR7M7CdFt<>y-pP*u=;g$D_T2x@h2sa8SgT* z!KFf2K#3@pOOu$jvQ><5IBY`zB4K@i<=P`L-j$GgQ66_j%Jp0VPDpZmwmo3`?u7Mr*1 z{w+j!QV|96c}82*Z@nK!WfWp5w0+Bhe}dS!}8-Zuw3xbdS$(%i;D`<_wH8S(=qH)evrcvq%&)k29lu@ z@&Y$KQ%a{~rvxlE7a||+g}Ash8(YE{>Jp0%mJJd*NRr{eu*5(ESt=Meih1a4Y}uD% zCf%f=aOr?AB$q~!IZ8rt!wyrCOH&KrY7kdgu=D7na{?AzAU^y|UF_+u(??ew|BC`M zTQYP_XXC5|Q={&KP-$W(@KnsbPc?q@f9!ScThc@Z-C|J{)ny*LtB9H3T220za&@n@>s9|y(1g_Nb}XRA#UShIpuWz$Q7xK@p2LglCLmvQ7ANP$!kF=v zbZ}|p^&5X+_{JtzoA1+EsGHVX4)}~Glj?eJK*UrdA)?v1crgO!Q zdM_{T_`VjIOZKhxY=6Sp;6Bhv8@1#D6iBRgmb0gcmwS;*T#<9!$E>bP z&%5*=_vxPKol*TkabK_UACm$ulVN`eQ~w zQrBbZ%fW+)~fH!Ht=qW-T&zx|6b+0{$Th#CGuD5^uYbe83vdqJv}{v%iEqO z0t~*&9`Q=;cB}&W_j8q!M062Oil_vTb?|DatqeWJgi5;))5kIDA2})s!_8Y4Bk5xp9V(K= zhm&P?TZ^2~-K#-{Pe-hD1Pl#jLy3Oxe1fB&>GO5)xFi$J`=+9Us+1MeXAC7XciQ*M zCaj^jZ$8GjlJP&s4f0PIUU0~jE>IcP@(uwAsG>YukElsv754@Uoe_FF!hppcgZ#S>UJkMFN zF7p}&kI_RDe5dRF&SmO>7!in7c|SOR?1OFNqFGHlAQ&`s`cOj;9#qn7^SjPoS*rhR zuolpeg968{ah__z5gO-FflU>Z^ZnE}S&ns6_3uFoe`l)N8GifKdY04v{LA3J@GZaa z!_ou*@6X%$XDs&3``x+@?kAt_Cgo7Zb_XUjUpJ zB<_AyUy1wexwQu1XTH409YkNhm6^FBlQX_iauL3!Fslgt*%x0lSyndqnl(vB@k8M9 zyRaA-I^7juL&oBq)>mPmu4uw*DRtriqZb2duECy(V#miMqR6-KW1;lUNe5d-i65KV z=$JN^-ZwZQ3FfTQLki=SLaDMJfKi5_t?2F2L9pAae%l#N(teONFG_atA()I_ z!l6lfFT~HkMwK?iO$q%TmXG`!{osp`oRX5GhJk^>gQNc`g`>uj5Qz{8x|sjO;fw)L zvRp4W2zKrJXP@3axOtY;N;MFJ2{Uk(L4Ow2Bxa27KSr4j4>r!l_D+l4{Hff!Cx+uo zl*}hH^Vu@~HfkvQL#ANxKI$Mk_LExyIwwPnkAX(E1^uBk@t5`dG^{8j`!WrUY<95B z==>0YpP$im%2djOzJgFd&3DN$l>oJzv=7UkS~7GPK5ghD8eJMM3Ipx~b*!K1b!wiEB#c#L&g!6&VY+bczRCZZ1Bfile*)OvMlMOo4r^r?= zXYZ>3{lv{MlJ}wjY6H0P;NkAXvkUINVBo2)Ubpm>!^z1gO^nR71j&~*m5Ml~F|pR{ z9Li7b)sJvCIr+Y6&4yNHzU>${q`emjdc0ZT4fsj#VM*X1$4#w zt)I|2Wxj2Ps@r`CO&kpR8YM*Kp4TY#(Kh}&9lZv^2rwC1HTSkbSDF)s+mNG+RjQ6? z+YnfY#wY=l+!h6Jdq|Rkc(T#O`~4=mM=E1<*;>Nu|UphfHFz_ z`iL^V+2KFW0rSgw#?pe%j!mC8!o9MU(=GJPyb8wUIdZB7sk4+DOvJP#`D3Zlp_=Zi zGpU1Vqpf{rzqd97FqhH=pw&9H-RAq%a*U;B83@bj^-YxOM7Fv3oC&VO1?a7xFh){d z8EA1Ge$iuxxSy>Ca%qaEXR2Ns&AbvLNG=!tN;sW9we(4MkR2>AgLGcPx&3;MIyra1 zYhHph1SM-KZb>^nc%cU?!?xswrM;9x5H=@m`qvZrn=f^z3cJ`V#o}(SE%% z%cW(sWk359TcpmV>l4hgBTA07CFBba9DpFGL}l!4qq!EWajdSc0++Ty#I|aD#z>hA zYhGne?6>@>g*O9-<7@3YON`U6w{7e`%-E%A$UW*7~Q0+xPYuMx4G&hLgvZhsnh5 z%A*@3Murpid=RPJ4!13H{Fql`rWIo@Xr99x_MUh*Y^#i^k^}sQaq6=^6~4yLrW1ON zRC^T30WtwDq{C#FM;8N`vGc=Yu(z{bN^JdSDOsr*;-XL4#71CFSc!>H7Nc*aIwJV9 zm`X{FyA_^6xCNU0h{i4#b9|ruPwV+YLoA&|$JezI!?i4hiWTZ&A?P~b_a^Jj%gRZg zbCqe)KbT4l7bx)9m*O;N7=ys*e##ZWf0C7IDo6)>1TZXRg|YQY*94V1jDnPFU38MO zdfho0@boR{OW#F~EF?@X)lC$hNdKJ?XsqO383%t+L!lbcpkP4MWQpS7+)Io;sub{8 z)gs+kV-2qSR4!dQ!i?NMI>k+6#g;;I+SZ)m`Dj81;4`+O*b)_s`2p|pe*fIlu+rC_ z7~Zv^f42g?vfy&Xb;22dr;nZIdUKyezTL1KD~w=Aj~Lvp`p)L$=y>{i zIr`~|5;L3QnSib{>$us< z&QmJ6&1_^`pN5=WH41Z`5jH_a6cR`eK|R4$lkN2V%bBFBFjXD&A5kjSv2%)ye6gZq zs$?a;a;aQVlrY*~dC59*1Dgi&8RaY3m$l;JUTSicZ;JLRDL?+4mDv85)wxoni54NA zQYuedc$vY%&8${kJ<-z$a)eO+_x*h_>#*Y6f~^mSUPD)I5z46gTJK}R+UpJIunDO1 zKByukmm`QM6NIePD~+ktEwyqoBBtOxvQrbyEUlDuu>R&P{WM}s8=vT zO*lzi|2$p6fP+MzjXG6p=?I#_;o>dP!ikZBt%!%9FO?!b;A$&J5r>=yk*tK-wl|`D z&Z+U7&80GnSCpo$lgOJwZZt7n4$AVi3o#i=e2dM_ZYDAi%XR%BrdCpOpB1PPx{GOnnNI@0+KV%K9*u{}4`MG$F z_s7Sv)n=YC^1pJpo_rlqfrGW7k{3k#-!xFxuj1#&cR1!SzYu{@hBdF1ien;UKujKz*4M1*HJFE(@wjO+R zXg#DA`1I-iAAKLzb|T==yoq$ueyIqsEPwW5iLCfMZpVp?{{!#TxF2brwCt1rd&ets zpMwt$T~8j_eit%f;o&W(oe#0g&^S<4Btaj*pEG#A6nROi^4n!vY#i(Bdy_7Yw0^6n zXtAL0ieZ)MjI>}BAoPmj=|HMwBM7becUgJt-RWw;%hN&YIG&dNRIB5cFX$7O+U;KF zA1DaIOv+V?j|5s_O^s^qQnmhzW~a>)^Q&dd9}0i}>Y8OoJUyB|IAmwg0KGJ-qv+I$ zCOlF>a!$hfJ`%GXJ4p|uVA(O|sxdRg+P4tb`lQ$LMK2NqmvuCNdF;#0(jzB`qk@hY zyd73C*xmaN|BN0`Z;#Ay^ByOSJun-YtHF@^NPZ51-{R`TMYJb2vKpChHkUaNh2R(v zd|B3zlGW}I4B9Vc@p&I)BJF{(P%@w5?P(Mjy3f!izn)4-z3jmtb6gureBD~>apfOj z^CnsmwaopqmH;Jw3bI>+#r)4(ZrUR1Y&xB~ks|#BTe(Vx^WMb}xZ9Lf&R~zb#-43@ zkhZ!V(b#NOpt@=dnkX{% zC|a?3G}(VXo7}2(12?AR{Bsy9eF~z*AM&Y1UPC&ZDH519ZEP{CCeEDTGH;RBsGVXI zhno7H2&B~HASgXAZ`xv7z*kCxfhKD+>T}LGBHTJagQx2Z(etZK16!EDZjOe*CvO5Qzp|p@VZ1 z!p5-IJPNC~Rb~=je*V(Hx6{&TKYApLCAomFN1vCk?#~>k=uw4;fOvIAY=cB3Er%^9`K_R@|6d z#4x-|11UwmDs9f$leUzN%&NAv3DGz71}_$Mvl}xXCBh%viZLu|RncnXxkU$^soSJT z{fufo+(T4#xN~}xO%H^{e1EIk3Pd5vn_?I8&49VfO18qH75S{YZRAqkQ{kx88qR?J zka;@k|G=T5dGjOcEsT^L zm>P_JPLKMdlORpbfi#I*0mkpxi!o#iY1+y<_mb)OZN?Xr8FfqUsqd1krz(}91xduM zMFYHrY8<92Vmcc0h!z9_)^!$AKs=TC2=kx=JcBr@C|L;{sV+1*#qf{VYBATnW5ji*g%$^AWJMx6@w!#^q;JZr=*p&64_vin|n#1PGz}(iz}yl3r@~qA|6`) z9}+c!7;AAVWc7jQz-RI(`?j*XPE$>5(aCh4=8M|Tswgq3hCV1C`aTTWAIAYVsk-*# zX`KPWW;6x~Px}c?zf_bK_*KBm!#XqBI>P^Z0ru0avNxHG?orLm%m7Axu|_X{LT7hJ zM~C3$D3|NWCv}vM=wg2Bj&IlQ+lk=2oO+F4azCY*_l9ek5F)9{4!FD)>t$`bX#v4M z%t%;#J@HAjr7@qVC&u8K%(wk{pC11wS;{Ciw?snKR0lyJL_QT!l1`(j)u$zNyJSqJ zVJ|8AsJ*#EQJ}Q(A|fr3UaRZ_jAVMZc#I)UfA#4I@SkNSVPgw<}}#^i^ka%xSN z6XhpDcC0}5c5|vxkI_oZdHKLy=Y~{T@qDm>-xiF`WjT)zD19>$+zBS~ z^XB9Ex){JwGVStPb#E=mA(Q)M5DXOx7T#AbC|&BH<;5V(cC6H$gSd4elcT7SSuoa)L0kgWq}T%%ywQoaFIY z7d7@t@*|Yk3wnll$T9*pwLs?7J4yJbj$TVWMrVPNTj|b9! z@ygK~uD5rc#)BEF#r#U-4b0BNTK1q?eFV)-O}ByPR}Z|KmhTLE z+4qu3FhePk7aQyd?)|c7ltla$&yA3|)xJ~@Zanfbq7f4JmB?*_$sbnIVgek5%GRO*@QHnKztL4v;6pUxz(K5 z35!4eB-fYszD86wTarG7Ex7a{+QR)&IiD*ZJIcCJhM)+O6SQ;jd#l4mCL}^8Tw8vS zb5{D=lz|dm>~9-}2j;@}MS(1$l+}6SKX29EnH)tc@CBL!_Yo7!nzIqMULYZ9^BHrr zS1B^*CCL(3{$*25zdC3f?{t30K!b)`W#W>^?baRJfZ>oy-XLZUr?-5(Y?4eJD3Ehp6_l3i}@Tym&&S zqe?^+X~HWv!P1ZkW!yS#L%sGjz6w3%E1SmB@c_06+_I!OgJg7ZJoingbMZNAM;JTv_tI0Jg#5>I9wX-FV%68mL}%0CEA91CfES5E2EPVgAM`nA}iuApHG_3B`3#s!2l ztQ3;w4@FLUE}xYB_vr-<;Q>m!6o-dCiWH6_NXjV*@NX@Pc%Q`0h}B9~rOO!I@Zd_m zczEMFUP;fve{gLhSgQA5fNli9N@BVmJ7oAww|q5HFCa&!Bw)t%A_(DnhT6`3!V~gv zZ=u=M)d@R0I}(>7>@3Y?3;hAxLXrI*^mQ%^$}|z_ptZJV(SOwE`?`E=*E_G3QK2VU z(s_Ixd3+N>sAvNe#~Jm&B+bL~UMf9^-$}@nj!I)2bl+uJ=Lz|kk+3DRcIeC*D>$vU zf>0v*y8Am?UQ=PDyXPT?HQ%>fM8;J!ymfTSsVWEYeOV*Q6f5@yKNIZq|N8As%UdKU zi;o7CmQB733;$`J^4ba5wI&x}HZ^!vcT(IPCX--o4B00odNdzILF~hBd49^#N3iq89oWcOuCQNN`2Fvku(Of$>Zy{FD%Nt3Re)(D*AbSPqBRrMRZ)$wN4e-ycm-T zj2Nz+NEW8?166|hzMzWuL}@bIqkEjkRrPg!q@vf}V_SvnsVR4xSqs25bmFth@|?(R zqsr7Teqxug{5qf<9GtA+-Mq1Uh>LiDO~m#hzr$eveIl2<9*QWl0;I`Y4EHe!aGj+r z6d+C;lpjko6|7N6A$#JAJE zwh-dZ0IT3GZ;5|h-kQW+*03HG52sk$tg!09t#h6mZ~EJgMN;8=Fb^x=vHh|47-O$V zrU*H0EB32t8I16Dp&O_ozM8Az@ivF|Jbd`jYNajct5Ibb&@!ZG+VUQ^T!OyyDX}TN zN%gJ;F-GeKUxGwj6J0z%wj-_M3n^J!$A*vg_%BVNwQN6eLs_Y(X6ve7mX`f`ky7(= z=rhrRWwg@+QMXb{Q(EbVacUAwDH==#spHMB%W1VIR$KsBgVYOlQ+C}A%aInb#V;{b z;t&-tPc004?W<;2#^98pWadyf`Rhtdyz+f++#Ys%9Z&yZ5OGSkf$2n=B8TF-0K6|< zvyF5#MSCy;3aNaMPpDc%=u5;M|92fn6cGcZr)`RV^j{3-l!+YAy~^S~|KTT$R9CNX z@W!EoNKN(AMFB85LTd$R+AXFwst#JdTI*veFXJ{F;~* z)MH_sh!*`7f}@KAT4@&jmdt$-~bS#Uln0b#qFSSX1MJ%@DO##|xwl*I9^DnWO0RMQ}$G)9j z|Lt2=B*^Karai0HeaJ%Ila{si+^vkM^-dFC=JuL6afH>UmqmTs zXL?)+tb4Uisq8g$JmJ|z6Zu6GKhuRayxxlGKA06dd@8smD7r(2%}imuDb1vv^sf0r zQV$wGHQ;QfnR5B7sa+VmplU>ExMMVjeM76nI6wZrx4DboCA=AU$N_5Q$%|S254Zy% zp85j+DR4>zJr9rmN{%&~6tL`|hWYWo@fnw=pU=jG5s2q>=|<%OLyaXJ>jeif`lAze z>*CB-9VkJ>=-jhiEe_n@KBoB;v;|e=Wu# z%>6{2q5^WcRi+y8LKpG~W@PyH!_Tx<*}GL;RhpC2W4r}2JP52?IfObe*0KAMdOw&R z&X^`}%Z*RC2-Y;eq*1VlgC9jQVSg?@CFRpsxP8h*M)Kl(M+5SPb)9x6G1XvQ9JB2a zX`EVwN}}UqgnO|8!|u6q$0W6XTf{~PDJFN)?EYe2V(`hAk4Ng%MDu2_IR(!-ZD07u zDv+&PSI)|>SX7f?V*9A&G5SZ z+?P9OI>Zd{*9!Igx7reHAM~|)`y!I-VW&|w^e~bepOr$FFOPyTkAiO#+72RD9esnC zh_WvP?{(xli7S79(o$CKM1&2WJ<*+DtPSYhrP98DKg&rdN#`A@e#WS+AObx>!CviUWh{kRI3B63zxGB7h-$Ln_$$Pk`wE+5tyW698BnGunW zwnEt>wh0OeWCRn_VJ071ZHSQ#p-a7;X$qfh%^XmxLMdK58g$;a-XQDg(Cg4Ejq0|Y zmXBaA;*EPMu}K=YGnfMwX+F zca?!akT@790d_kTWptxO#9@7(NFsn2#^|?oSZea2iTDl?Tk2f$TN&Dx{db*bvQYw` zfne|<)Idkv5%cF4qcbxze|C?K+=AThoKx9k|L$lCKoZ^){8t5?LK)aU%M!DpiSW`; zPYsn!9$0JQE8DB9@|X}Np(`fp5HI{{GK2H3RTTX08?v1R0VikIGc3h_=So>ane#y# zi8S$QL_ss!_;+^3=juw-47>ALN*$@*4CSj)q#3}EJ7zfiWWd3&+(0# zI23BBnZo5}5EgE2tJ*v`oi>t747hXY>QX3*O5rdqj}&QH+mjz03lXUNc*&haXxkJ6 zVlK6ii5yCJfV5hWB>0mE`)G%9@fR!eJD5ZM5wYJ^KyC7-`iDV6K37HU;qQpX*Nl@s- zQ5`(ODtln{i?N7TU5g{&c5aivsCnKNAg;eJE&zh-VO~!al6EvOr;v9PxczK@sX5u= z@2;U4820tm;#=m)$w{6;eFi>1p5@2deJ}yA7nJB1Ft2QTKubtmVthdB@YR)gUS6Jt zhQ?D#$^Jq}n*T1ch=z{N#@06c^Jg+@9;0@k)~;Hh;}v$kACjTLX1WQV8x=nIWPI-Y zg7YLO_+_?hlAL|AZ0fd@S#}-IUQ{)dS%FKaS#{n{=59^IYf{sSBZbQ!DtTDb)7#L2 zvJ?2B)fZ@PIOM3_E#iBK)aAkBF{S)IN+mlrIyIYV!ocRq_ezmEeMo|JovUL)gaZ*I zr3sPz#5B z=5{6Q%pTqTRn~f8VqHlCX!rDEj=zhsj-!t~F&a*~>8=3}@3ou=zG0aeOmFYm!{esx z9ESU%u8Y0E0}m}MUx)&AM;ZnOJ0R>bJOiLkfKwz6*kR<3IdVUJ#8Yy3c&I#cZqsj@ z_5uD^ovHer-&w6NI(O@A1trZHo&N6&c;Kf;R3DSUpR3S zxiF^!Fmg~tGrP=hv!w!y90JEV-!1bIj2O;U?_sl&OPWSPiPpj(f`t^}(8=a9 z`eNFCdZn<(+kzqo0x~-~NKTjwUpchjs#uCK;!@(LgOmYZG%BD9X_Z>+8Z4c!bB>(C zk{o9Pkg|7wc5??4zSD_W9@+O}nMiBS?D~cLBK4`@g5SuQ*PYg6_G-|hB_)*H7@n^SHWLaPlkoR^C*6M> zd1G=u(tW2-v}G=G-qtz8qLm(hkD=n;I-lP>gZ(eBBSNet^x0v6XuY)&ZVz+jgS*op*4sECT}r9&a!`@x;W0SJm>nDbyvH-V-d0 zb^2ZH(UAS{=dWF7ik`6K=&wDt!v#hLn~`cF`35A%Qg#kIF!E?m^98ZB@y~)rRWDq2 zh|HM^P-};DejJq2lFtSsei}V4h_ws_ZUNow z$%dyZRO-Msp>htL1_>`y=X@pYD!lv7)(5b0)kACq^$LHml$t|EJ-(8MtPP)mU)_Y1Yq^Q=r)b=5ecsXz@SDzWP@F~8I4 zGv~KJ9I*;08p4FAMP24QBZ%l4&ZkYU+wo*ZFTBJn9r0FY?&a$ZhA8BWJ{Sc)zw~R9 z@txxXb&Uxr$LzzsJsV5Qkmco9yXTVVqkISaVJKDP%(H1e7y{v1&QW8kPLsI7+y0An z87TYfeXgc!|C2Ps^zFd!ifXF9{lLQx;?6ma>vLgdWZ=HJl!(D9X|~m2OLe7C!w(}h zpO|ediczZAs_}okyQ2Js~*-Qk7SntD3w5LTjis@T0``HSS}26Axzri-)n2= z;%pc8m%cJ28q}Gykq-`+zDZYB)P!S0QTKQipXwj?zHFRL>aOalG$<9fc4E?lx0sOh zvHJKPI$_`Dt4i8g(u9$;)n}D?rvRS2=d`f{SzM3Bma%^MHsCmaI|VuB@vpo%&#e>( z#<80mc%H&Dj4Y_5D3yc5_>IGqqM{-=75`Vd52T$BMsJekaB*>8u#&5!>mS82;S-x= zx!&!=mearfaffszFA(E}=JwhZ=>_>ak6W$$5^h8ao+w46=C$w>n#1u^GqtThA=GU& zUcy+dCHwL}Mp%1Kn%j}ee)i!+U`akB#J$;{&L?DU7gn+vEm8^2dR|FGTIusxTmG@V zRylXM#SBrjl5icdepC;k0Z%lj!q7y}?m8_!NZ*WOlSlkxrEFbIadeJQ*t1x^jDrOA>60D7jgvb^BXPARfs+ z=5~L;rekxea171e*-4EcC102ZGz`iCOa>>o3Bz)rYXph;eAk7`w}SF=x8+x?2?+^i zuua76+M_#WK#AGY9&&nSJo?)eMwG)4vfvK`ju@eaSE?(f&V=`xi!lb}_Z9Ly`}x!} zy>>8Bl^3Ne#VG4+A5@XOaT6aXkMqlU*^cW-rzBEi#D5aX9GXb^0Qw_h`+cW1?lWIQ z;8_jIqc+63s`Yn;;c8t)1_a2dAW2D{3pkQo$t^E}TYEa4%+%)bo(@L?O@!rIc7Hre z-OI05O}a>~9CfyWY8M9tH@|7G>0RV^fK-WZd=T2o_$~v$Xfv-Nk-(E!i=Wa*un~b) zMs|wS@;4LqtE|e;FzP*SpEe_B2Mg0y9Psu6T8ZT(Gcpq(mlqd;PoIuHm{gh2BkcH; z@yt(9NT>spH0G*kq|J0a{(b$U6&et144dsuI-3CSNV3^-p3YfHiZ;9hy3})-V|GHv zviT$@+`L2ITKnA;w>71@jJ+a_^#W?muoM+NhHnSWirV|_Z46JX)*)vj$FFVmv)@yN4UNwKa zc@~l&jbD%DPFF5vFGl5nxBGYQ(R&5?zIX=_s!RS*VYjrxz9e4L-t9Axzf^hejI;b% zRpm6T-BrW^!<@d6uG^)>fc>Z@V##$B&))!SOk>mhhn8reGq#%C`g;`i8q!0JB7wLo zUX#tG!cl8%nb2G~P;!z3^;)*>6LOyzgkV5|B7Np!!Gqql?b^Pus0jA+v&xnRP8R>i z>7jkVsoA)@IRe<#u+`plX}aGsEO4D73u96k4$h)?2yVlwHZNDsJ5(F9v9Y!`>Ijw~h!* zT#qNWg%>*UidC1Cx1~B{!Xf|>^3CtdOE19QSc5TwQWqMe10l)s^{}7k%|Si-ffL{~ zpGrvdfq==Q*9VmNDZZN2#zaTo2hz_+9w3a>Dm5@IR6y6+9kPN%0a9x)V7bSt1a|7l ze9z%hb0YwZ&{cbP`fZRx%+Asv)Nr&QU31@I|4e)|HAzwn8N77IF8_Am?#|gm@VMXX zyG`qr;Aw2Le~hxm?SH@jXRmx4P1ii~?|Js|PjJAWoubCQvSxp9kJ$iT7?9plgq~jk zUdlG$6|94Myu_5@-rin%iOOvldjHr&_ODVzW1j~wFk^$6ZTAe#eW=j0+vI1!6pt&= zBaPdBeEbwp8rH{D9%Xn92LGK2KIs9PCeDiupV8wRC3gQ^VsA`}*~*e_ji@x6LpM4t zcg73Qs@95o7Xx0+(w!ZNNsA2w6t_&V&diHB}KFnF>P_05H@5Ya83@! zQC)~cA|+~}JszYApZaWPfuOy-&>r~WI%udLes9=M%>Iqh+F14#^BLbq?|`NgWg zU8j9p1FJnoByaY#1Ph*BQU7xkL%MIig zJ{0ynS{?%1uPE!SB#00(_+J+b2?{E(C9Hv`Y%v`9<55she29yC)B-AJ2xLxec@aR?M0Ec?(dl@4(OYfjlAP=!Q3Ec>mv7 zEoI;u9@fy;IzY`VE8`oM-21A?$X$?cP%@PB3l4{$>{# z&CbiI-!>y>08a%A_{HMxG?|wZCR8nBq30!L$gEzZN!?H0W9;S|8xo z&W;s8aM0+h*W%#+L$4v^V6R3CwF7@9A$2ipqXB`SE!%zGspw=>k;DU7``OO~O+hIs zsiML_aSct)`l?O)k5d9b^~()2-1kD@%bI3YiGx)e^U@m57K^f(MG>~cTs=|r8s{HO2it9yb5D&y7pa+;fOsEjsC;=rdT>iX*EY{bl zmtzyG8;O6xSot$uU#!O^b^SBJ>)x@;_sgOw3=Hf7RV8`7`lW7{`X%2VGb_q;E`#!E zc(*`{a@7(bohv_A07ZNodlh{#OTXurpGQ}@)SS@%kjFYke(uQ%6Hh8!UpxLB7Z{%A z)|P_4t*u!dw%r_VG#`u3JK8LRs62s^-80u&XruM`2f_7N+*@yY!zJ)W^!MlAj3iJN z?JCVB^eq;^MPIqpyl!U~cv%}=#U9%@u!qCGZarf^MlV7IEk=B zOO10Az(V{Z65Ht%R~9Xzz#aYs$(L>Gt7ILu>$}Tu7I(*yvguxf02AZ|5(69siJlXM zkkeAPpprzL{+>s{W_P#SvdKEAZOMebUT0D6g2YpF>dgK7QM#wN;|G^yrg{IaV`y4) zOQt)|_$|;HZ}l3sqrok7x~Ij3KufO8f-l5gljiQyVlAGz#RrDu7KJuKyA(6d&Z9QE za#MmT--*U<`Ga6}K^)`sI4Cr)d%yefs1Ehk84Q#LLj3D zbK8Nf8%k#T*+aP!=<6+Na`ls4whtUeWD3W?+E0O?{B1*Ka9!g-Y5nrNpdb2Z>;#fk zEBD+%FkI;{!iZ+3p-*?LL36ZYC>Q9jdnG5@*|Kmv_mF3?<^3mpa;`gsKU4z5mJ|xF~%;YY6 z9G-X@YUpFdE)h;DU`gR#m^e^gDnV@9UQbj!ZdDA(78j9t}`<67u?2^K7J7RkOs@(o}Vgb{e5OW zoV)v1`Wt%Y?OE!a|8((BAeQvS_ew%XUU0I1w96=Wr~yG7FFrFhRveZr)?3bk+d@)n zjN_27#PMML|NO}n@IA=fV-rtvyAMV%2nrIjQhWx5;5L~~i+rLVIfo^?v6Vv+I{gos zJ6q)#k3%0EYu#MSo{#d$*9g$XhWdmoRLPj*I4>bl$3~!6c5$iw~UX zyaevEIxM{uIZV40P>3SI^5ES3=d{dT5NKu4eZw;Oz5=dGc!05ps}M!hj%defSGLr` z_dq&N7KM{0V+3PJ-R_p2=66-MlQV}eHT!rVzzynpao zr@!vU8k3_~^SrnT`npDjldemxCdd(9@g5f^dR~?^6^qV`zd!zmN!^*(-2bUJ8Uy^O MK7XlHreGQNKXSI%%K!iX diff --git a/prog/python/qqgg/analytical_xs.org b/prog/python/qqgg/analytical_xs.org index 9edba88..9ad2969 100644 --- a/prog/python/qqgg/analytical_xs.org +++ b/prog/python/qqgg/analytical_xs.org @@ -756,6 +756,18 @@ We define an auxilliary method for convenience. import matplotlib.gridspec as gridspec + def test_compatibility(hist_1, hist_2, confidence=1): + mask = hist_1 > hist_2 + comp_greater = hist_1[mask] - confidence * np.sqrt(hist_1[mask]) <= hist_2[ + mask + ] + confidence * np.sqrt(hist_2[mask]) + comp_lower = hist_1[~mask] + confidence * np.sqrt(hist_1[~mask]) >= hist_2[ + ~mask + ] - confidence * np.sqrt(hist_2[~mask]) + + return (np.count_nonzero(comp_greater) + np.count_nonzero(comp_lower)) / len(hist_1) + + def draw_ratio_plot(histograms, normalize_to=1, **kwargs): fig, (ax_hist, ax_ratio) = set_up_plot( 2, 1, sharex=True, gridspec_kw=dict(height_ratios=[3, 1], hspace=0), **kwargs @@ -774,6 +786,8 @@ We define an auxilliary method for convenience. if "hist" in histogram else np.histogram(histogram["samples"], bins=edges) ) + print(test_compatibility(heights, histograms[0]["hist"][0])) + integral = hist_integral([heights, edges]) errors = np.sqrt(heights) / integral heights = heights / integral @@ -832,7 +846,7 @@ We define an auxilliary method for convenience. errorbars=True, hist_kwargs=dict(color="#1f77b4"), errorbar_kwargs=dict(), - autoau=True, + autoau=True, normalize_to=None, ): """Draws a histogram with optional errorbars using the step style. diff --git a/prog/python/qqgg/figs/pdf/cos_theta_hist.pdf b/prog/python/qqgg/figs/pdf/cos_theta_hist.pdf index ce6724da38f489f8cee0ab6fc5b4890aaff7d32c..964003756e90b37459a191f94861ebf8e4b07342 100644 GIT binary patch delta 20 ccmex1}2-QL_cN$0Av*imH+?% diff --git a/prog/python/qqgg/figs/pdf/eta_hist.pdf b/prog/python/qqgg/figs/pdf/eta_hist.pdf index 229aec734232b74701e5d4d80a4e9fdac6398ecc..a2ff154e118f636a9a45e2213ebf332c779372d5 100644 GIT binary patch delta 10166 zcmZvBbyyVs@-<7x(hW-~EDgI$qjX7kcY`!bDY0~ygtRmwAtfo@(xD(qmmtz0h=A|v z{oecDpWOXxpZUz0nRDj+vCnKrHQ=-gP}>PWD5vr?yDUxCWZKqi#2pUJn5vnitByS5 zFf6Ip*mN?i@V8a2U+-cFCM0t3nG`ZNE=u<-5$OX&&jx_%7s7*X4&=ImPAkm1!hYTC zpLd!6!n(bQ(!aTSaJeGHiyL%xX5|%(KydvU{1u4ae*06MmlqN6<621LjFs*9EoHtC z!v5L`mJ@~OG3V$QxZN1~ZhiBM!GG=6y6(KMlD{x-K_wiq*Y0?Ek;Gv3Yw!H${3OD9 zs`SZ?P&rs9kgxcJ_c5N~E=Wg(VEyvCnc?>EbB)CH(Bb#_8pKmY4`oN4+S`wnM<)ZL zhvo>By(Y{v)rJ>T0z}~3?M02T9FAw5h;O}T?wf@bTk|_R`iZ*6jX%2cv6y`C1@(cF zey++0o6-cn$?<4!!fzq}nGPjxb?>HIj`_S4|DZr@pyE@Af=5(p3vU9;hbw+`7CQH+ zCQ%bo-;eF8&*mHNH{a_LWVOl=51&n+R&Z$MNE~U(Ka>r!j8)*S;p+mLiYyPvJ&JQf%?Hou)KzgfZ+HbVWG!(?q^2*Q zKVD|Ko{*OOV7kjPQ2}wpl*^MK?yQDo_8D_A_GTnA+C+GjK!32B{906GSqaOF^}7)M z$z_sj$wninh1-9W(I2=Y8OvYkyGGb*JY_O$h-7tNLggc2hM&< zPhfwe~#8HmoA$?>zNdy42K1+$)vWR<0FKC zN|+_IyMF9gP**eKX$PFPfgWQ_3@ofavy)t^XV16XZ>DMyWL^Kg&c^8|)mAPVZmElP zL|=i8WZSdE4A0~bb^`uupF)g}&8bQ2@cziUw4EBTr!>{^aRtJ*k@1I23HOuZ?~V5_ zs9N}0Hz<+C9!#$w0+W`;oNSxE1rEnd=t6c{n-}cC9V_+yO2MhDUWqjgNf_sQD`ww$ z%_V;6<99v1IZ5KxbNOViGfJI##7>=>=2K&SyS9~or2p&sYRk9g*U6R=ZuJsb^rdiO z8Dxh(V$tTM#>Yi_b=kGTRa&ygsY1o6??isAK8V_7=(Y?NMJQ0VQ-cV0=%W^GzWhn- zMA(N2vlIXc#`*6C=eyNVK6Y+N3a~4DNRoUx!Zw@fU|IZUj?*1G&fj z+>8`PBoEAapcHW#C)+&}rW)oZXz3XYsclhr6~&6O!V{@o6$VR`HU1w+@&JV=QPVeP zE$fF=!(s%PJrjUtBp@2i(3uH~8f0JQirq5Q4VlNW`2M|^5Z9c{h_WNHl?#i`!~!*| zz!->~ntMziQd6qnOhv2E&q%Gpa17m)?yN*JKhth63#?buV`ui~DXhhvGPLQEPJ%9(|o1EZl~u@3#Ky%P8co?q^lx-+aw zA1PH7Vz|gE6o#&e6w5;25K+nhIhuY-jH2b8T39V|;OJW4m-sprF=UJX(;+#gsv90E|4uEL=#Li=|7`}#S+!TgY5l>j zupSi>(nX}eObS`6|Xxh(zBqM)hOWUt| zr*i+x{ve>31RPn8@$8d>0s8t^EGyr_8`8C?@}uOUb8y$~kXLIw;!_Ie%1{~p3l}l5 zfu0i6uQLJ$ekOj8#EsdKt#~eaiFO1T<$lbm?{-oZjLID)x&g5JC-)O(3lEEe)f!)Q z21NK^E@lkGOE#PczR3>GE@qTfim}-BE6<`ib<-P?w zPowFc6Jmlq5Z@rvXA(R46&qP{`Ssj~Yt*^NMz9R83zYE(Y@f?Ih>E#u3Wb5(@BIzv z&LwJd;H4RM*TRcykQ2pjAML5>pPFIE7y`H8ZHut$S{lNSLnXN{IeIefl=a7_FzeOZ z`nxSIXVJSHl;yt3iu~B_5p{6e>VbZAnKM$An#CHwL|`4|sy>tfr*~@TB@imGko}Iar43rE&(w?dnS#B_J@;72TEOM(Wpu-bTyUVDOJ0+~Vq> zmp=_DGHnC{jvDd(K(`f909Gm^OYfBjL}#v4dvWY>3RexUS9%DjdgtF=%(iUwspD^u z*RG8!Vj>9k#xIF}YI1^zia7HqpT@$$8I-pmwVZLiiwpS{mM`%Ac8V_m?Z^iNZ^s>L zoG&C-?yoaiFsZOeXmf?3@Grg>>qfWH+HD8$!a_XiYWpp^N~z!nQhm(m_(N|r$KJMq zQrsg6HvMM{2sho+REW2|ek_^bFkqf^=FV&Og(3J|tQiJt$doro-PM?8;V2~1medMQ zGTn8mNy3qD zKu{y;@Rdsj4kn8V*UEff8h(SNsUQR*gB9bf`O6o(T_2xh($;~B+KawI)S*RJ^d(ls zME5W~VFiKwo2mRaMDCh>Y#cgy3eQN2ACx0RnQL+{Qz{7ZBMwR1lHDNpKYRKl7J4*@ zP7$v$jONbX_4i@v7lw>ckT>rzT524ak#(Y~0fuNd2u-ebv}y^?R7a)Uhn zFxagTMu`7rBBnwcVx1PEQ!sU%o2^TyR>4F`PmSA<1%2J6o^r+Az4Il!;`C#jm_s%~ z@q?6cR8~!VD{#PGDThm?bIn6upH-acE6Kz6UPZeW97?Gr+QC`N@3KmiFU=Rw*>?*B zRt;x@mZgHM#K$aL+)nORl2+$}ia^2romfi6W_VPdktli`?YoqK(!QsWDm9Ew;1SnKTwA?Ctu z$)*f#SFN$NQllFp4t*C0{HcS&F5~t+;YaVMOCIC36iGOUB!%I{_!KTwaTnL7c^(}} zd1$B;X%VoVh+5)N>B9q2(zPa{bhaurJRw0w@S;r;9?P%30$QsW2E{Yzyy+$?=T(h5 zo)7@rV{(=HkjL0cQHJk&Vi%_nZ8B`FK-sy>h?FpIrLxZZ%sFA{)i5lU43Cm&*%V~C zgO6KY>>`NF3ing16wj## zs#QP9NG57C;T`+<(~;7c;9$1W+`hWWi>PiL+BZ+^(@qt*v8>71yX+fHR(@Wy2!7j1 z74(DEzTin6MGpkc8|GV-yC(SB%HTdFi+}u*D7FGyK$dNtlHrz&^AiER;_kqH@1`~| zJ5pDTk_89t(a!!BdbkGypl2?}H4t2X0SfcUsYazBc543qBFuK&+Qatp7=xY~I6-O9 zkFp0y36Gk{?LzS=Ncvbx%;me9FD+bKVnD#&V}tMV&O8`$)D3)rXFX+CXY|EjS!2hM zxoYLpQ-bwz!IoDsrl|CwuiJ3o;3Ius_M^|f)!AG*QYICBxJ+*mp>xJEyplT(G@`aV zjBM1Ib2U@pc8UxU>-Pj8q-h%zGpZSSYq?F`}%3{j1I@&q6oOY?qQBAm~A2XuJ7>U zM08anl#PEXAfzWaRW~ntcczJx%oj~6iu`@C(wsiy!NSXAj_7-^OHMKj;_C|Pa-k^G zV2ZF=Io~u>@RH3ePOK^Z_VbezIV>WIsT5wxB-79D%@%lHBibZA9-$bGa8uL{E^82& z1-H!UG(9AvsvI5J1oi9VPHe+4FF-_;1tEI7?Ht~rYcyA1SYofh!LQk>WEUUcBH1%! z<~UXE=u?;$K=5U;V4XzKsGs`rG&Y#m-n>E(23hBXO=vjdwh%`-ycPNa3SYygkLAEy zn^6co&tQ9lkV|tdHG%|jLcDRGi#fULiZvP=t4W!p%DRnAxc_X*W9o{JeB!De`-%*2 zLw7+alaa{;3=~=gbz<6Z*hXsL8JVwe;6L+>jNsm=XH{NPTMdb>x94=86zS~@orTsGik~#JVLR!=Ba%8+E}gIWLb;4P*Y8+y^HBe z2l0|_hV3OIa0)s%F3e;yq2X6gAmK$$vqdHgCH7kAvLe~SRc%o_!%FZ(&BnjV=f!ga zYrfXVRF=V(F3t@V&XP6#V#;Jeci(NS#EYddHW`~Kef9Gmx(2#qe(nwznMS*JPwkgNMYZ>B@p zCY8^Ao`Ui6e+|wTJFsPSvZvrw8fz1-e}?J%X14`uPms#>ia$@a15l73eWb<2%~0wI zKtGj_Mjb2tIL{00^TmF3J+FY%$Bc+=a+*lme%UUEC(x6zCdYZ{U2_o}yROv;MvItO zj1yoD@3*X!uAX7owU2c+a6ML@>m00#{W&U?GvnBGaZud>GRrB`)AUQ(D%VPpeA05M zqpQemm_!y^jA->J(Rw!HH*y#!5qNpbYkm$7YOI)SAF0+Iu65V2RnK`WP=N5A#}ks- z%eNT6y>31<@47k5NMaBVyu7V3zd8Rbfmm@`fSJFrozkhMZWc#ziu6eCxIQ}mXla9Z z`~$LZQhcpz9?<}R98J8eN1I;=dHrTi`5vB9TY5W= zKG4=CM3;PU!<>%>bR^ahzpX&L07MXzDk1nra2h-OvEY2UX@0KoWUP-phmq(SI30gx>!oOkiMLkEGr6t<-(2k zLJPxG0>KYUG6HyVm7Vrbi~2OGafjlcO=5UC{?5g!o|m)>`6rk9OW&2w6P0h5CX#L+ zjj45(Y)d$z5he+fm|l;aY#}$@Al;kmCR`^< zCbr*RrtM8Lq`b*^sgkjQUYvTWwK~1$OeS3?LgrqdJ=*U_HqqybP8~9Ky{utvP0g_K ztJvTQ#uR*U+NoF75j|5PAIu`GZ@XF1{yct1P(*xX>>sG*HWCW%{5MoQ@hc;IPwY1mIpFCMa=@nce+|ec%B1!=ChE*r9n;&z z%~IpW)}f}4scV3r>Fi(WGX0~3av3~`3V3G|1v zX|OoMI0N1G&M*VoVeUEP$twHL$;v6p>PuC~pyKo@IWbh&YM)DYr3O~w{rNKFmB6z+ zQO3)Xy>gj4T?G2-P=Q*MN6x0?0dY6b2eoA&&{GC|*s2oyX!Qf< zkP!;g1-1p|1Hj_CcqVVFxwjk0NUTR3Icp|)We&yi9M50>}pB7xerol z?@3Y^%LVs+zE8FXplhMLh0GV*|31mZnvxljw<)}2Zy@JDFkOV&WBQ-{rdN1ri=-EP za0Jc0HRB(n~W%(e?;lB}o0DkR_pI<_6?+!IZ zH{CM)KUl2h4bIi;cHc2MqI=KKRMmU-vHV{qBnL2x-!p3w#b*4sRR0i<4lX$|avRV? zdQc8U=p!}}($A`_s~XPafTx^_Sbyo8`mS909h=fXuNodg;P{MN18JK2X)iNEe$Vd zyykS45j!-R%v=k6CrTtsgU|LYRq3Y9Z6pH_aeU9U3d6E6AoSy5$;+Dc!ty$c;X zvL(5yxJ-@JMPLQB`Xgx5 zQK}0Jv|xyV*7~V%XAl=^fKGv$%Zy6#YtW?_0)a}xLIaRhAv4w)F%Ly^UFSJH8^!!; znfacYA+)yTqSb6josBHX5KD@y4V0DIyjs> zr;W^DgpHt68-Suq!fDd5iw&CVSa|mw>Rd@#c#0;O_Ua}}Y(3+(OokqoigFfOp+3`p zK~qeo9;5Y1lOSH&MoZpqkZsInCG|Fe_t%JZs2XEKZu$_Je7%TLaL7Pvq7kBOpnJWi zVdU(ZvNWbu1HM3Pl%4rb%zyo|8$WKt_7>J<=-4g7!m*rco_4Ti>)J3Q0^F)8>ZdII zB4P@sB8n|7reiTMwNK-vIef_CMb#g;E#mvIPmQZATBvYz7;>7fQ&+uH3tOH#_ z=I(8s7WGbmYN@0J^-lbOY8(P)n#QCvFI$erZS6qC{)TN5BXLV(AjJ`0AS~|h^z&0d zfK~WF9q|VnYC3Lf3H}!L7{e?*=~3lJ4q5s-W$$Uy){0aNN8bCWIb<~y>}pBZ-?FG? z5x7){+{YgHK_nUj-V=K(m*?Pw`!ZzimEhUNqx-E9yrj-fU&13jHC>Aku2KvM1(Myy z0eF1j&o5LpFUXgNokP~|XDIZYnigCg6>xvj5rP?$S#AUv_j6g{QHEVECfJ)evr%HOXF{TVp55CNv-7%g9n>;Jf(Siq)mtsAbjoS>V&cRxw3=tc%3`jO1;5z^- z?Fukc03)yM%$ic z`iSScSM7@kIrXhVzwGjFWp|#jue_hqYb%+51AECOXzn1N_QcVF)18o{u}MLCbbQRR z3J|&8>t!hNAvXtYfE_$!2wyY3;+U6=+73NxNNm`g@r#1S`tr0E}u+bS~o|f;JgFhusfZ_?bg6 zT+sOws%6HAtXDBl@ik5$DW8FZ5%bAjE<2m-1R; zb7X7S(k|MYABJp0Z|VntsQGpd%v*lZAF4n}1<%tqi02MU1JE1OK!kPK&o+qPzcvmC zZ*)e!A)d!>FzNBo2GQ<13e#H&;0@g`Upl=YsuU5XmrO%!&{=eJ96;wh-5z2Ze@z%N zc$>0MXQBG`n99jeh(2upb={D`Xv;A{EqKu{`aCCN1u$nrn^gDQ@i*}6SEjiJ=XNwM zLqU4!{qJQ%7Rq%SO2&j ztx8dV-s&bh!#7x)o+hdl!gk`4m2dhk+H1&IPU zkMBjYM(@Uu0x4`pEJC-K*wh#u6y{nn1;jie*D$H84)2kY*=s1y=$e~&iy=Hh*cno7 z-zljV6qT7Su`-HzgszEx>|Os%AmxGD6I>Kr%zW05qh}@?%F$fHrk0GoM5agkP-d51 z?>j_GP#A{8QpHZ_I(q2zXtUTF6DyCTLG71*?SMSXONZPT z`4Sc}j8WGtJZZIDNgO3kV@`yO)jF_W*ZfIG-}`=Q%*pBSO&jSMIMJ--`ow*GQdTMz zirR4A2d_$4CLksPj(HE&cf2EeRB!`VM)^x#S~xPci3fvO)EJXw$`YkEk43Qx1uJZc zokXqVM{cbvo=mqyXZ8fQ?s=hA#%D4aeEc33v|%(%3Z#8jtswn2TsRDY>MDfJdg7B8 zdO7T|q0JUe3Vfb2q6w}AGpb2l2aw+?KAhZPkEC6zwRn#O@IYN(wkWT8_ed(&@(CV| z{Dc<~5npspU#9FcoM`y)rx%!SMx(TA+gZ@Wh+2Y#`$(28Kcp-?3lPTeC{X6Ewx1`o zrz%^R7{_}T#Rp-uhJ?v$moZG#!V*PwPnlB#gwcwZY_UTIKe_B0`Li-ycR-M%+JNU=RW1UmAB}0(^I3LPB@3LL$O< z@`Rxxsg~Mm2nhdQPf&;eKlDzFkMHkRPzVHe2MZMt`8x&?5rPT)T~7om2>TO5{C|#6 zVUd4={~Z&6-Q^8~{>>2z6@&`@-2?{V`Af6mXzWTox8Nb$Q(X2UeyAQvd(} delta 10288 zcmZvAWn7ef*Dc))Eg;oezeqEXxnw41P zTivEotvSVscl#13RmaZ*?644bZ80Z(Tp~l)a_qi|NmIQbiih)K6tgU~ZuPUUY#);T z{?*VA+_}QP`hBTiYlzqCS(z%ljy!5;1p4;dM$-AQzxzu{r z3Io2f3^*}g1uf~njNZ#;i3g^xud8WXt_ADcfd^ExuSf2l>TPIrXW#Mn{XPA)rGKLX zx;Y_xMOs4mXi2^D_uE;e$I|#TNpLCo>*Mh2i}O;l>%+H~hm-NHd@g*)L{mULf9vJ` zUPJtMWOF&@+soqM?Id|KZt>2on0E+)6| zdg-TBJ~OQviPLp#4eT8fJeFzlp`TL5RqPx1PX?K<+o{eFe9U!KqZvZ6Gp_tE_`B zz#Clwz!X1(Dq!T?P)>VgGO`W zjyB+D`WM*Ix1T%$?wwvghlRiHn|TA?{HF3lrh z(pLOpc{a%S(2Sm{nizk6WwvqQThj0Chp4uhB?#}=G?)CZb&UB!FAk6E*nZJ9it&7< zMT5zVF9DZv0qN6&PEQ!J?TzoLm|mpHK{hW_arHxFS`_q_NPrD?0wH}?a_N2)D=*zJNA4#(ci{zJ6 zDv4j?R1jbOQxIYngNCmZ!>U3ew)-d9GI&NqF%!_=XMDcN96oRTReNS$MXP6}c-fCJ zDBusOML@Kgqk$=^auD9z~E9kP7ci}rvzOY6m`Z0pxvDhx|Ec4_5 ze$9+DNFnaMb?lx=F-L>cG9DQyv^T)*lav&V@yFl=bBBZKUx>K*uR-)qjBe2v%VAOp zy|+BbJKcx*pMwaR`Y}akq?T{}LXx_(MhmU6x^}j?mx+*j#ADl5xpxZf!x>f;ZW;Fn zhtUf;rf0G1eZ*O?R7{t(Pe6;%<5&frx^YnvS}TYxLb(Z=2XI7H)!p7+Y!-{M5er)@ z?Fbp;f6}wdQ~s5{%*{a6^-caIJGUmT0d%E)*GnQ8qKMMqPFJ-InL|65N%#Z}MIj5f z>P~0f!Zpvyl!+<0j2uGS#1N` zC^_!NfxghHrueLgNIhhV`l-es7nS&IS|j>^ z`D9OI4Zu8}`<4Wp)Pvc8rj{nFSyti>nN6O4-4F^|unKwWq)vcdM`XGb&C`LqGx`L` z>J;B_7bsGll4>5kjcxUhXa_ut$<+cKcH=uBe(vm~$ z&n`NVV;yml?_=iJ|8~C^Dp&By2I5!vdim+-U z&)aF?Z-rH34-|3bKP`_*Q8TYIrGLpA_e~n4T}q8EavK^?o70+$6;k0+C($8huN#-o z59-L1ppCn}cdDW<&t=s%5PwZ447|}X_E;|yY&VT)t_=Y7VGu_v58@NY6H5+%y&}M6 zTn?R8Y%cog{Jk|pVjnXwZ>3^x*luXN!2HyLsblcYb z*i-Yzic;o{r@#89e|CaYh-HXXb}aRp3#)hk08=4J4t^%8*k7NgY2hi)Ji_LPzfBtL zLv3U|*xuIM7vNk{dpqsG4K6Hz71s4Zn}6n=)Ivb6oIVD#*B81WZU6k7q_{1*EVlx= z2qFSrGS!?f_zNc`YmHpS4S=Gw{UB@SAs0wlP7L$=oJ3_H19pXpdJxM;wcB1j5=q{A zF|5Wp;>CFvf-Z{Q>p~RaxelnGZ!$U`Dm za=Qx{*{aL$qi!oI^pr>^k84g-Z83&L6rC(SX1nQ4R>F+FEp-aRm9nh5XcG;atF z65lJmcV{N1gFM5a1gIs(YY|F!Vtf(NslHW44kZ!ov_3YU@bm3@Q$IJ7hTtea|$RZWCFfn9{XxTrFukpXWN+uxDVfaIEDyZEmPK8B9dBh z%VfO7`Yb!0MLdFEP*77B)eQ>cNE^V;RO_MWbf_Fup}+Z_dRB=K!e{za`_g? zdk?vVDGRbsaKz&;ftdO+LMF7~U_W(V%2At-t{o`Q;0pO`S+7AQ0XUmKSzF5^J_p}~ zg0|UszRj*zK#=Fb0(Jt1m2#P~tZeo>e-)vam#C;Yc!8DUjC@K|!E<}eLX7zGIyv9g zskC-t6Femi@2584pPqCj9y=gZP)x~^XLY!=SphBx+l0G-d*A}bZz%3W_)W~@fu&oY z&Tt)*|9h%6X}0d)nK7l}+q^#|H|Ne!lrux-ma&frfV))j^@HXy%qm4E#Cu+3nk+vv;h%*@rk(&aB@J})T4 z+u>P7=8K_=KNRX#OI$R+t4hyDV`~j(bDBheKS-qXaWRoJK{7Z^+{`)9F*fmLXaEs8 zeis->@?z~Z*M0R0WF%6{Ibp9ffi%>?bLQJB({(OoKqqw{6_FVd$hkT>T@)hnf8gOD zSU63zQQz=X3l$xTrG-SPOH!)uki7xD;n{@OxKTbFQ#qKA@;t-)`i9>_KL-`*@wV%| zzbGYIY_)~+au%tX}gXKsi@~BbeJF)Nbe&FRhFvueEunDPaTQiGV;V2+Z zUPkR+T1q9EB7bu->y^Y7;UOMB3>TR20-fRNy({QOI7SU-h%r5oBb7^+^*RjOfxu&R z9Hc*U=HwWYD1Aa29VNV0`gAY7tC)*WBV4ZFv8ugvni82&gk-vT#t4rFKIsg9I-lJe zP#R>E?9Hf5`GYf}!>}=lY%wNTW_w_^-EVW1RB~F`%9}?cTvSqhLCBtMvC~_4!&{2P z8H$%KW%ow-u3^lM;wR*%@sZ#JFBI8_ad~)ifxBQ|uQUfsGzUs0^ zJRTsfGKO`=VgEKB0+^ki?)p$kHWmqcA6?Ib(wwm(6CYq9%HuBO+kID4W1%i^n#%fbW#7#C^SF*n9V*(+ zf5R51$LS?~!Qk7WM_4~d0(K>hwZJ6u1O`|z5dX)*+s$7u z8f9KYDN=!v(}`@Wa|Q2A)H_K49nG{Q5U>9RUA?V;Wzb%$$+8UoP&3fqTtmyiOdm5% zYq+&iUAAs+|7hxrm>mCw>o#B0OQtL8ho?${0*v$K(4hJ?5e2&F9E$0C%$jaGX6CGp zDMU7c_Oy+oimhvSc_QD%JYi7(VB`cz*z3H~LvZWk@TYHc8co9;l3k`M$lEY-ON*gB zxufB@;n`N5vJXnmLF7OIzM)H_0J*eGyAOZIl1pW)k_Gfb-{&PV-_JgA>oqaKSEdTG zZ!0RsCrou=R5QZbtr$aUiY9$#bWW0|)g5;gLpjqOOM{Opz7zdCEig!HuuMiQ6WyzO z|MWFS9qwcEOkYeTh7XNme5?c%u9sX<8>JS45ZrQuO_QK|Hp`BHBnDFi*~V`gqc}@r zgJ$fq*w?o$D3gJ^yO z(m7#@77U9NOXq+)x~5`V$_yIXyUz%%wGU~;+AEr|+rWI5tB>-5_SM$RR#B*-ngz1TMsV54xS+ujYX8m+|cSO*1a!QR=SfIF>zHr z%@~U5-Id9`ybESt!CiMT;IAs+54vdE53xQszUr&HkIB*g_&AMsTRY)y@3ed5C}Sua z=(aI`SN+ddK>?dsh^YwIN`9i&PIz|6y+dEUP+F4}C@z&eK=Xa@r1EQ{jB+3|&5v0> zz6w1tm-fVo_e#3W84yM~j&>1>iFFd1Z?eB_$)#KG|Xhk!lw7BI}{t6TMCDg z>BB1ba$U1}hNuzaefRA7b&J&%QSTdhiFJ#H#|g;zr&mG;(kO~rAy&-kC#ha&D7s_rb$KWQJ<&XN48)|-T` z%DKa&e^>1cdQY)96vs+SegMsyM#0xPc4QNxWcH)#!4%Ii)DnE!`3AM|!S{6zbAg{B zfRh|AWKyWED9cQ({DqBphK}N(1D9}WP&6Jx} zSi$_o(`?qrG#54bv6 zY4)@`q_^a+Qz|bewB;|WKk}5hjHpJKjYVBPI3ZQPV7u2PpXzMK7Tl*rDhf8J(3#c_%hEpk7 zHayCGH^H*CME#_}zenWB1<@y=eOI2$hsM|P#~O2UTcQrq{6t#NR3;NjUNYM35`GxZ z&XrO%8JLB!Sl?RQ=I zRd!r#_4lLaIWhj-VV16}Yo~I$f%B#k^5U7s5#Od#@a^RM+Qs=o?iPE1o)UyWIS5fb zYkz>8pYp@OnUB?bPC<+nXZwQ)6L6IueK%1U^QLmIx4mlpvxH0iwN1z~83Q9y%_mSZ zl8wRH=AgI+8WmjYZD43RJN%9X5m7QlTK+ACl2j?v4DpXX#ubYv=DT zl;g8oPEykDoX;P@8CQuXR#yj2ulI{9aT^>ubx=qXdmlgK2o%1^B2!4E=M3yag&b`o|Z^j?B@CgXwg?kqDjR9B7gCGVn<@^JdM_J z(8d}d66WO${Ao?1dZf7a!|%VX72nXUQSJZJTFN-7$i3kAr&R#kL8h%3R4DecGVDJy z7-yLBWF4)x(S5FKLnzjb<(7jA)qT*wHqJH*&iHbZE6FYUlL^9mLL}R#w(}*dpT^th zQ1erR=Cg@dC0%^oz7b7M42~zpgJ&kAK|>INj|w>v($HTlz)wZH;`jpRe@MENylqlA}MR zs*%jaD6CYkdIdYNdU?9Z12utHXr<@%av*)hHo>JC&p^M}rJpf6rH|6IHX5_%s~1mR!Y z1Lelu2>Vea?I{al{_6Vv&FJUwV&I8lb+l5cQoL_tiq93^JaOAw$6QhL)de3ROjJ(9 z^j)Cak?G}b=sa;2z0mvGzY54VvSC@HNql8XGoBQz+|=iD zdT!rxzVlHID;;?sqXE6t2#hCAG?)~6L850`ooxp#O}FGepxxHq%|n%-4mr&~3hG&_ z?BkA+)NeNZ8R>AyGMa$kK;YrM-o*;?f6D{V)SDyr8-|waH>uMIxigz_wr^1T_9?85 z?A-{vh%SbsxU^ONcH`hJ;g@RJoSSBdq}2n*-;Und-dHxdWRkL_E5;4W6&B0fB~G;T z{=yPs$I$AT%l9@2VQRXEA?ZubD!31Sd%v@Ih9B#vt=w>%3E}fKX3Sx8jv!<7TJHu4 z)<e3$Tn5+*p!`+{e`S55C_N_<`P4d>03OT>bY$z%YG#>(&0ls{gyy|5@s? z5JrjV9`HFP!4C(0^2PY){y50kY zDHoWot}nc;#>MssI;1!DG<9L8(caQriNQFo_@YLosoeQo49IR$R@}^w3Y1m*NS!y0 z&x@5hjnB$sR|{~TYOZ)st9(g)hBJtlX+^!I0%k#cC>S|Ow!W@&dd10#<6W z_h|{tS7I8Rs?~78-#tF&uETSvIp_f3fv+0hkW7dQZjEO$Lj@OIaMzc7ea~)+2---C zTS4#KWO!NJ1YV!mJQ_#h?AJQJOnnt?ZsF^(9`*yIWU;=i!-6n&px`arp~)9qBA|{F zK0Eon>FdYq!*<~z?uSn$`*GK3QGQTB*A!!#ugZ$>!w9NhYLGv~CSZ!|@IF33wP>2m zyMUd!5;pO8<^3>OU=2gmMGjvz78f;gLUWRUKFZ+zZKCQzD%80rZTTcOcxa?al> zOu=3ZDNsNAiA3qYve02B9OMe=^_~)cp{wHVm#eWS8ZgEg((8P3>~R&hiqM|2tyu)E zTOQ{4k02+4oMNwAHcYdfR zm`NkA639!VrlLG&;&)EYyPu1g-lb(itfG4OM2NV^6nFUgI9Bj;?HAsjWL2^-2=zh& zmf(Zj-9ucJ0BR_KrF!g59vWu%oA3kob)6vX=3oIC|pPV%p z5+WfSLDOAUD!EcF&H2EpklIHG-c_dZ{$+;G4Hb&r)1UT!GrI9(u#wiklyn0 zCoA+sigcaf$!s=LAv$ehRdK7exZwAmj{R|NO6oJ#%)h~unvv@OwDxi9Ixb&Q`FOo; zGzYbvaSHD>$2u-DlF|0wYr%ZZRF7tj+{+-+hd9wh>&W;8>R1eA^>>*bA;}WBqPoQ| zRYXApuutBJLOc%=B^sT819i3{X9d|91E&t*_WgDLAL%_q=dbXp1u8-f%r+!0yxkG1 z8Xums$GYOow(2+o`g^NHV#L~%SO!n4Ji_Z;_gI@)wGB}zCIHT=xbBSh7|cKMp`ISb z^RoP>gio@m(|$2X3m?5#6K?f4Mx@4;lnT&^=HRE4`!I##V_iFLXg4>9>v{_T z9hca3uzNr9AnhjDc>k^I>3!`s3IFHbG!GZ+-mJgURu{8@Q0=OEAl0#sv?f=o)P%!- zsegK$CR}0p8d!}rJHt4_nrV zsxrFI&VaHeS$TcG)Th4L&QI=VaHjR=m!F`7Ogs=m;hzmeI4i)Hw`F%ZkPNEAzFt+r zIZol66|%!@+WGto!cBs? z2M;M@KR+6%b}VsWqkf}s5;;l$_C^#ZYm~4@X#-C2>AzF;b@U7rylc63no&>FqdD&l z?k-HyJRePWbdXbmj8G1<=igT~K_=6I=YnN1Kop7>M2R@8#R8;17*EsaecTg*#Zd{_ zqqg^%04T%Yi>0lj@&PI6Bmr zenkI#P2>e0(|}r*<+{DB92o0u?)YT`Nc7Pd6)%P zsd=?AufcpdIbOsuHBbg%8B!(6_8Q%K zrG>F0w~Q{sMwf|W4(OGcMEoE$*N{JAl!3;};>W)=&9YLNN<>plByIpR?g0aO2@*$? zN1`8N{rEmdHOegZl8}F*5*1t~n&KJ5=_OdH5H!r8nOywYL()yFD8Di6@qUCgJn}=p z3*_+6o}d^0eyKZYNM&vHoS@1eec{)#s#@w`R2~v;9cQAg&Pts!+IXbizMGLvBUJA+AmKX}G=`?NFzUnOhtdO) z6hG5>#kbyY9ih-uS=(dR2gdDsSQ4$4X68F66&3Js8w!f(d))ehDlqSqJwaJ=OjC)H zGz(U=wNfz507Ipmx_QBjdB5!`=Uh7OV6xyC2RwsX(4@C8ca`-aP&YN{*Rb1DR0 zYAe6A471470XUD@VGvxlmrqrH?3U}1EM+-7`ca#hksgBi5`iy~o|CRXG#yqUkM#)7n3#Bw>)?gF0#w|A!Hgd;}a`M2+zwrICbT-(wE>Z;$l$P>0TxA3ZLL zr2iI8Yf@jHke(hgb95kJKlW*mq6vG8N;i`!twyeT;lZ9V#?atT&D)}fKaR+*KKV@K zT%PP*V}5eA{z7EqV(sbWVPoO^=g-yBfd~R&1+)J7fr&z-J{)s_1aTpvV&dXS1iD24 zOcI^_^Rc?H68-Z71&d3G-H1s*C2z!}#BcUWNkVVd!9@Rxp@}5K|HP6+^dtaDslU5hF!0|bPzdC2KOj&s*nf!s%K@N$fIy`nf13mmg-YB|EGmA(Oo*uDjU__F zz@j&L#h^EwfryELfg9F{{r>=AVmH)@L&1ML0}+?FsaPBazTuVxSp2^(0RQ*tK_o

0] # filter out cubes with zero weight + # print([integrate(f, cube).result for cube, _, _ in cubes]) weights = np.array([weight for _, _, weight in cubes]) integral = weights.sum() weights = np.cumsum(weights / integral) maxima = np.array( [ - find_upper_bound_vector( - f, cube[0], cube[1][1:], tolerance=overestimate_factor - 1 - ) + find_upper_bound_vector(f, cube[0], cube[1][1:], tolerance=0.01) for cube in cubes ] ) diff --git a/prog/python/qqgg/parton_density_function_stuff.org b/prog/python/qqgg/parton_density_function_stuff.org index 7c27f2c..2dd2e36 100644 --- a/prog/python/qqgg/parton_density_function_stuff.org +++ b/prog/python/qqgg/parton_density_function_stuff.org @@ -27,6 +27,8 @@ from tangled import observables #+END_SRC #+RESULTS: +: The autoreload extension is already loaded. To reload it, use: +: %reload_ext autoreload ** Global Config #+begin_src jupyter-python :exports both :results raw drawer @@ -691,7 +693,7 @@ Lets plot how the pdf looks. Overestimating the upper bounds helps with bias. #+begin_src jupyter-python :exports both :results raw drawer - overestimate = 1.0 + overestimate = 1.1 tex_value( (overestimate - 1) * 100, unit=r"\percent", @@ -701,12 +703,7 @@ Overestimating the upper bounds helps with bias. #+end_src #+RESULTS: -:RESULTS: -: \(\SI{0}{\percent}\) -[[file:./.ob-jupyter/542b03d025920448ba653b470ec6492cbdd1e4a7.png]] -[[file:./.ob-jupyter/d47db0dde9ae59979f271a7cba8dfc46be3f1dd3.png]] -[[file:./.ob-jupyter/7fe9d3bd60427cf20af835649efbcbaafefbb3e0.png]] -:END: +: \(\SI{10}{\percent}\) Now we sample some events. Doing this in parallel helps. We let the os figure out the cpu mapping. @@ -718,7 +715,7 @@ figure out the cpu mapping. cubes=xs_int_res.cubes, proc="auto", report_efficiency=True, - cache="cache/pdf/total_xs_10000_000_2_5_take11", + cache="cache/pdf/total_xs_10000_000_2_5_take12", status_path="/tmp/status1", overestimate_factor=overestimate, ) @@ -786,7 +783,11 @@ And now we compare all the observables with sherpa. #+end_src #+RESULTS: -[[file:./.ob-jupyter/800bffd38987eae852b74eb2483698169c08c4de.png]] +:RESULTS: +: 1.0 +: 0.78 +[[file:./.ob-jupyter/1b1562da892b4caff29588dde88a04cc3256ae35.png]] +:END: Hah! there we have it! @@ -817,7 +818,11 @@ both equal. #+end_src #+RESULTS: +:RESULTS: +: 1.0 +: 0.86 [[file:./.ob-jupyter/45cc8fbfa523c8faf704505d716ebc299a1a44fd.png]] +:END: The invariant mass is not constant anymore. #+begin_src jupyter-python :exports both :results raw drawer @@ -860,7 +865,11 @@ The cosθ distribution looks more like the paronic one. #+end_src #+RESULTS: +:RESULTS: +: 1.0 +: 0.82 [[file:./.ob-jupyter/d834978baba81dba0e2f052f1965b62ae736c151.png]] +:END: #+begin_src jupyter-python :exports both :results raw drawer @@ -882,7 +891,11 @@ The cosθ distribution looks more like the paronic one. #+end_src #+RESULTS: +:RESULTS: +: 1.0 +: 0.9 [[file:./.ob-jupyter/22100d74fc1e07aa9ac75e22dd38ce49eaad89d4.png]] +:END: In this case the opening angles are the same because the CS frame is the same as the ordinary rest frame. The z-axis is the beam axis diff --git a/prog/python/qqgg/results/pdf/overesimate.tex b/prog/python/qqgg/results/pdf/overesimate.tex index a37b11f..44722a7 100644 --- a/prog/python/qqgg/results/pdf/overesimate.tex +++ b/prog/python/qqgg/results/pdf/overesimate.tex @@ -1 +1 @@ -\(\SI{0}{\percent}\) \ No newline at end of file +\(\SI{10}{\percent}\) \ No newline at end of file diff --git a/prog/python/qqgg/tangled/plot_utils.py b/prog/python/qqgg/tangled/plot_utils.py index 0873836..b636bcf 100644 --- a/prog/python/qqgg/tangled/plot_utils.py +++ b/prog/python/qqgg/tangled/plot_utils.py @@ -90,6 +90,18 @@ def plot_stratified_rho(ax, points, increment_borders, *args, **kwargs): import matplotlib.gridspec as gridspec +def test_compatibility(hist_1, hist_2, confidence=1): + mask = hist_1 > hist_2 + comp_greater = hist_1[mask] - confidence * np.sqrt(hist_1[mask]) <= hist_2[ + mask + ] + confidence * np.sqrt(hist_2[mask]) + comp_lower = hist_1[~mask] + confidence * np.sqrt(hist_1[~mask]) >= hist_2[ + ~mask + ] - confidence * np.sqrt(hist_2[~mask]) + + return (np.count_nonzero(comp_greater) + np.count_nonzero(comp_lower)) / len(hist_1) + + def draw_ratio_plot(histograms, normalize_to=1, **kwargs): fig, (ax_hist, ax_ratio) = set_up_plot( 2, 1, sharex=True, gridspec_kw=dict(height_ratios=[3, 1], hspace=0), **kwargs @@ -108,6 +120,8 @@ def draw_ratio_plot(histograms, normalize_to=1, **kwargs): if "hist" in histogram else np.histogram(histogram["samples"], bins=edges) ) + print(test_compatibility(heights, histograms[0]["hist"][0])) + integral = hist_integral([heights, edges]) errors = np.sqrt(heights) / integral heights = heights / integral @@ -166,7 +180,7 @@ def draw_histogram( errorbars=True, hist_kwargs=dict(color="#1f77b4"), errorbar_kwargs=dict(), - autoau=True, + autoau=True, normalize_to=None, ): """Draws a histogram with optional errorbars using the step style. diff --git a/prog/runcards/pp_phaeno_299_port/Sherpa_Base.dat b/prog/runcards/pp_phaeno_299_port/Sherpa_Base.dat index 090ec3a..084dc7d 100644 --- a/prog/runcards/pp_phaeno_299_port/Sherpa_Base.dat +++ b/prog/runcards/pp_phaeno_299_port/Sherpa_Base.dat @@ -23,7 +23,7 @@ (selector){ PseudoRapidity 22 -3 3; - PT 22 2 200000000; + PT 22 5 200000000; }(selector) (analysis){ diff --git a/prog/runcards/pp_phaeno_299_port/makefile b/prog/runcards/pp_phaeno_299_port/makefile index 082ccd2..bdf47a1 100644 --- a/prog/runcards/pp_phaeno_299_port/makefile +++ b/prog/runcards/pp_phaeno_299_port/makefile @@ -18,7 +18,7 @@ OUT_DIR = out RES_DIR = results # Event Count -EVENT_COUNT = 10000000 +EVENT_COUNT = 1000 # List of Runcards, see RUNCARD_FOLDER RUNCARDS =\ diff --git a/prog/runcards/pp_phaeno_299_port/qqgg_proton/MC_DIPHOTON_PROTON.plot b/prog/runcards/pp_phaeno_299_port/qqgg_proton/MC_DIPHOTON_PROTON.plot index fd8ae94..7b21d53 100644 --- a/prog/runcards/pp_phaeno_299_port/qqgg_proton/MC_DIPHOTON_PROTON.plot +++ b/prog/runcards/pp_phaeno_299_port/qqgg_proton/MC_DIPHOTON_PROTON.plot @@ -50,7 +50,7 @@ Title=Azimuthal angle between the photons XLabel=$\pi-\Delta\phi_{\gamma\gamma}$ [rad] YLabel=$\mathrm{d}\sigma/\mathrm{d}\Delta\phi_{\gamma\gamma}$ [pb rad$^{-1}$] LogY=0 -LogX=1 +LogX=0 # END PLOT # BEGIN PLOT /MC_DIPHOTON_PROTON/o_angle @@ -68,7 +68,7 @@ LogY=0 # END PLOT # BEGIN PLOT /MC_DIPHOTON_PROTON/total_pT -Title=total $\pT$ +Title=$\pT$ of the $\gamma\gamma$ system XLabel=$\pT$ YLabel=$\mathrm{d}\sigma/\mathrm{d}\pT_{\gamma\gamma}$ [pb GeV$^{-1}$] LogY=1 diff --git a/prog/runcards/pp_phaeno_299_port/runcards/with_pT_and_fragmentation_and_mi/Run.dat b/prog/runcards/pp_phaeno_299_port/runcards/with_pT_and_fragmentation_and_mi/Run.dat index e349df8..9384c75 100644 --- a/prog/runcards/pp_phaeno_299_port/runcards/with_pT_and_fragmentation_and_mi/Run.dat +++ b/prog/runcards/pp_phaeno_299_port/runcards/with_pT_and_fragmentation_and_mi/Run.dat @@ -29,7 +29,7 @@ BEAM_REMNANTS 1; (selector){ PseudoRapidity 22 -3 3; - PT 22 2 200000000; + PT 22 5 200000000; }(selector) (analysis){ diff --git a/sherpa b/sherpa index e4cf1cd..75b7ac4 160000 --- a/sherpa +++ b/sherpa @@ -1 +1 @@ -Subproject commit e4cf1cd24ef048ed99a4447b51c3c76e1116cdf2 +Subproject commit 75b7ac47eb8c0f592660419dad5e87678d5db789