From 11646714a8e38a69555ff0d13aa5ce182d856032 Mon Sep 17 00:00:00 2001
From: Oliver Kirsebom <oliver.kirsebom@gmail.com>
Date: Thu, 19 Apr 2018 07:53:14 -0300
Subject: [PATCH] added macros
---
analysis/VKAdd | Bin 469720 -> 469552 bytes
analysis/src/VKBaseData.cpp | 38 +++++++++++-------
mac/allowed.mac | 23 ++++-------
mac/betabgr.mac | 67 ++++++++++++++++++++++++++++++++
mac/f20scan.mac | 28 ++++++++++++-
mac/forbidden.mac | 11 +-----
mac/forbidden_idini.mac | 50 ++++++++++++++++++++++++
mac/forbidden_max.mac | 45 +++++++++++++++++++++
mac/forbidden_min.mac | 45 +++++++++++++++++++++
mac/forbidden_suhonen_IMSRG.mac | 50 ++++++++++++++++++++++++
mac/forbidden_suhonen_USDB.mac | 50 ++++++++++++++++++++++++
mac/test.mac | 42 +++++++++++---------
mac/test2.mac | 38 ++++++++++++++++++
13 files changed, 428 insertions(+), 59 deletions(-)
create mode 100644 mac/betabgr.mac
create mode 100644 mac/forbidden_idini.mac
create mode 100644 mac/forbidden_max.mac
create mode 100644 mac/forbidden_min.mac
create mode 100644 mac/forbidden_suhonen_IMSRG.mac
create mode 100644 mac/forbidden_suhonen_USDB.mac
create mode 100644 mac/test2.mac
diff --git a/analysis/VKAdd b/analysis/VKAdd
index 95c13f6c34272a49ed789a761e9c0b386a807ac5..90db4f2610301d8c6da4a51f799431430a8d68dd 100755
GIT binary patch
delta 39903
zcmc(I33wF67H(J1WS9gb!<v0z-xNYbKtLT9AtGo(5(GsEtH>Hc2y2+IgKUA8G%ApQ
zAe&$?LI4$mdIc3RZU_hx5FrREpdzk$e=VnHn4tH4@4Mgk`d(9KKc`MDUEN(Z9RDJ+
z@bk#rrh0npisml`CBC4EEAfRDClAz&T`aNK#clz_E;hEYleIm&-l?5jb9(H~+GF;8
zz3HB`dlsj>SF~DBtjew$m+t5hUNTRLFO=Hal-inBPjhLB&lo+c4~$P-AlY{)c5;Yh
zyPuJCm!iKFrfK!H>Y8h|ak6?AD>IUUL-e!^V^nZUJ!i7<d~kQYbgFSKcwl_#XsKC)
zN!6~fSQDD*vn74FBx}6n7IQL;3|Cg`v~0<CD)y&!rCPyENgrvGjE}3z0LzT98a+IP
zqa=5p((0kK3RMs_6@4e@+*#agDNkO$fi+Gq$dvRRD(ElllJ3rsbcZIdK0!PxO>-p(
zx^bpP|M1l+!!RYT=8)kwb2E%SHGj|xW*S{;wbbLXjf`3?!wXb3DX?9kvYLgHjs3N<
zcn;%hXN9Fr7Us1mW%4UX@kU#19IPFp7fdux*6yy4e8#9>r@QW+ZH%fjFw8#_PoTzJ
z<+yRSPOEUY%Eu68u(XckrY$IHUbih{g~s5J5FJ5G3uzgju3~JUgv(Uh$If;O3MFH2
zfJ$nHaWZ5e%Q5Z^9i*3J8_$Qf?3unm%2NNaA>%xpg(_J@A6HA#OP-PPsHnKgs;LIZ
zG-`%**R?su<6$jzcee3N*ub!>lcexds<i@CYn2%+JcOM$!opkWX_Jkf;obGbY-443
zOZ3-+;mhL-CrZ6W(uX!6Tsq&Pf+FYdA}HQv-J^}!5d(EC(>NE=Jw8$Og+Qgf6D6-w
zy3(d5U4jwDwOtv;<B`kwFpTUG?w`}e$SGG&m|xU0YC>zhMhoN3N1oMRh%%BNeOzBX
z$9Uz@1+3g?muRpbjSmwe*frx~;<I|axkh$UhmbaN-CDxZDWT<SZdOzzES=l15V0J%
zW*klGsn5OB2=CKcZ!p*B-6t%n1xbGYYemJSJ1Z(GdM!G4;MX~B?X9#q#+*L&^=fmA
zHus6p1IuQ^=Anb>vx`1Vsm^rg9OJ9hGrD`O@kZYv`hht{>|-zK@6IVY{Mgq_Z$HP_
z(to=C@f@Sk<E{0EbBp>tKHjGHn{T}TWVbr4FqZ_!<RBs~Ve#t)C=RTcZ`2>y6|FF0
zV7HJmgoGdx7M-ccjYD7RKQLjj?Z9Q@;J}`>FY|DT(eH82_}PdV)U)o-zKqUKbRHmC
zo<TiBbmZUC>D<4w(5gPqC>zwPR(0^CmqDHTah33S#sg{HL&HI@97B_N#@e(8^(OO*
z&ZkAO!0qU>35RT`R?v_JJRz-y4A<l58F@oG)LSv%t?lnBtXxP~bSPobf6DK@gtl0;
zT~s(>o^fqRLml2ihIZ4F<{3{8ZKe;JSF~_w5DVRgN+&G(v3$d~6%`|L6PC_s7<+lm
zTWE)E!}hVdrMx<G4mmHyp0B7#PFVCu!l5dq^Ngt@Vr&VEelcDc(a010Oeto=5y4GA
zxE}0!B5p+RkLRxk54_|1`kHn<xc`X}!7Wc+53Vtwp)+{?Sqzb0i-$TNT5NBausC-i
zMoPl0-`$w*g69S3n)cA*HyUn%Ov1d1;JL4X0AnSsZ$LezVJQCgYUnCI{Y^y$*M~Ep
z*MmJS$;01Xi!Mr`irkwpO%hz-pS<Sydkq#|n;tyCN3(G3Dj)52nZjKn8W(&?FCMr7
zj!3@;)D`t_`JdqV2ci8>KIi%S>;zxX?bgt#>gp>0?VMH1$}XixW!_UMlYhh9Et!1N
zvVQ=ssa9rX7f!2I!i5ucsS{UE(|!teJ*Z{X2yPkB6Ik$qH((_=#?71OeNS-9!x;fT
z1=r|}j`D8o6Hj&?@J=pQyctW1C!kmpidm%zi>LEUl_V^#Qjyr6tN#@IcyC11@8I>|
z$4g-4LGALk5xFAAp<GfnXzLq1vQI@+=m@;!FWvcfMMZGTTn}R5^?kNsA*z|MxTiB=
zaUWMgvkwvu{%)K1Y4EbQAqSPrhECXOeW?g2`uXDNkgHOWn=bn9?c%r;i}?u-{>}cJ
z-<c<P*^hs54>cZgUiuc$h~SuR{~AyVgc1%RSWUZ>09eS28Y4HV=XtB^y0i&0x&KhX
zIlp`K;ALH96=#SlPIX<nTj0E#iUx-<Bmmt@Wzco02?W%vC>@Ig`w11mZ~3IJd5xNu
zUaJ2m&)bNWl#x_18_c=labil?f*b66Jwkaeu2U@r&2Xqk9IqenC_8XCc`O!x1~q+a
zFDhQ1dsei`9ap0@t@Znm3itug{y3O9ux7p+6Ox$39xCIL*k~s@mh*`x-$_`6|NmG%
zR8Y^9S3QHgjHqEnq^74W#>i0*>X_tT8P!T(ciJc$)!eZjHZga*Mu%6~#zB3a(Ry@q
z$0r;vGX{?icbw$lWn<Rprd5hKh$*{pbXb))IgE+36!0MqV^+R8`tA<j0LE<imV09w
zy8V^V^1TS=p~b0~+9qj&TK*C!`|U^qrsD&D2>C<hE6x~Cj%n@r`zv0tL*;u=zU9Ia
zt^7~q`9MGoE-k8l3fic;tM}sJ+FLGk((;p#=*+XH?Gy98hGg|xMDzUsa3SCREk2`_
z?*z7QC`9)M@{}g<SzOwEm&+M<jBU_hjL;CyI&)`Bvt0>?9@MUKg(N7HZ$E7$k8P(X
zW*AGxHmDy0x~LeM5nZYLE~v{6Bdfdil<~>fR^67I<Z&D-@5al;xJjI;NQ`}pcLt?5
z{j?C^_ozf3fv?*^-0l8TMz?XoZT7g0k1f7YQIRt}R12OvwSq4fZop!BmPH&Uh(B=R
zV2k*uAR;qo%6nPF4j?8hotVa{ZEqM;#>Y11X%v${Q{-$v?t|Karc}%FcBgm+0t)Un
z4vcSHZ_9CDbKVU_)%Xn8`N_@$I}f_l2h<U{7+&R0BOtw%e*8<LZF)<6%a_K8^hZ25
zPe2FW3b^i~AJKcu0|e>E^3#en&&&Ewg5_^3)+pfR^Nw-VcP@X(|D09*zd%VO2n{hL
zmMYRcAm#FzImZ(kr%ej_=^Yxbmq!1LSgo;=`_nsKge&aDfEhX=Cpd>U|J9e$cz-E>
z<Aia2Ldg@`5QSIz(pPZiMRWN?@C#oV%(>z$e~Sl#(VDR6yM#qwKk+22Tw46uE%ICr
zOA1%(8|AycFz%SxR-bypcyQvK`s>Gxi4&jDANj&4ofvjU0hY(X^G>5t%C{UhZcL2s
zdl+_LG-1)(vd10%5(A8nhXc82NJ)DQw1h=JNF)T=*pH#39at%cA8&${kdXbBw~hUn
zadcA1U74Po>Fu@5nBW5wyYeaK5AKSZ`_c|zVv^<>X`dN8Gw-=?YRg#kohjuBvkvh&
zVDaSGIY(LfL>Rx+7(K0g>ZhnQC11{!OGmrJRPpy`yn-HC!M}I~`|+q3yO6Akm3u+h
zAAmmy_CTIhG0_v|eHA?Slwh4H9|i}$*%*dK^3O)1=sqeNMO?s}jk<h3c&EH2sNQrG
ze2#Rqh%&fstNITDv#Njd6I8#R(Jm{}6L{C{QecC@OW`~PfuH!KKuykkc1c+3ag|4k
z#ys=HrMpFPHGlq}qIze_Ym?k%06x7gIiZSJj81vRIFa?BKJXJGcJjSdoO3XN{BYDL
zOYPszCM=+m52B&4-dyw`5_}RfX8ChRdAd`z%R7Q&w9dhM><t&{0$*D`?&G4XlbbOe
zW+SFHX*TU6NG{E6$eV`GSlcdg%^AR#c54`m`F3Q&y@xUb%4-^<r#5TeR7uAv>3&K&
zj!U~Coi3$Ijusu5x`ox!@R%uk>WoruH{jbQ#?<NcBd#3*hT*pr<t2tx@&U@bWxpfF
zp6T^F&HWJB&Bq^v#JgnX4;=ild<^_QwCG)cy%v32K9gHPtni)!@ArP0`C0j>BS!5R
z4fSqEj1DuN#1q}>8Ey5Q9~PC)2w}j=v+K(Y-gnrjG4sBf;|{|DYLAMFzOwF|kv21o
z*@~vkY`|EIvCb2zzxrO$8=l7OfoYKSHV=P{h`Y4!Z^SwHJs#XUJd#+Rhd7>t|CEn9
zR8evHfQa@?`ICo?db6_Y9Se`uQWKWm^*74ohZfy<XwetFntdzgi-UjJ%6q(9^wzBP
zjD2Gao)d#`XW{R=fJ3TgU`+6<;as=z39DlqpL4teJqX|1_u=HH$95Z>Iatc%Y)TUV
zy(mdTkMzA3&NIH7J73>2&&Zhfb;p$6gEL2V=+I;Eq=`+Z+!^~s+>`jz8h<h~CJk?f
z#Y<e@hr0}(IGpSselblLOI+YKF3*pONLSYRjOu&AtBADyg1DZ1%H>zi!?L9JEjjAF
zOL{L!clVHGj(pvt#i>^SB|Rk_Zw6jnkX18Hk;_Oe!|1Ru8jrjE7q)7iLpsV9y}ZyK
zk~VUc?E=_ZnC0ek6u;<QMVV5|Rc1aq!Y|#;Ww7RTPq;7W%NNt{cME#2)@~zxQDT7L
z8E-C%=tJ-E{BxeD-nSI08edgqT((;%yH%N}%Caevo|bJie<s>fsMP(<(rX-g-9fLr
zT(hNmtSa5AEK%m^y&Q|MbzFj_L87K<S6{EFI1ae%ZLCM^prOAXN11!L!c&o3Mbk<_
zh{ans*ZUO}hw;8B5#?39e8@$)1&X;JR8)jw`c1>QU5GLrWdX_@luco%;Nyx44R|5S
zM3k|Z6855WVH)d)wMhZWDJWf-2MbWfq152F<d~<Tq6~nB71vc5a-no0>2WAyQMysK
zN12F{H%1!Dktowq=Ag_$xf*3I%0iR{C=a77L|KNi1m#VXr6^r+T!ykK%Bv{dC^Zij
zu_XXpD37Cbquhdk5;+}ZI?7VKr^-Epz_8Yh#rnNH%3PFbC<{^Mpe)5h!CsVCFF=lG
z64Pg_6Y?nAql`nDh|-O6B+4|DIVe5p09FIYDMtV(b5UMJS%NYYj>}NSp^U{lkwla^
zDAQ4veuq?|bYFo3I7~y?4`nXO<0wl}-b8s7Whk1o<a_9&EW~`f5T)xnnxv{WM2iKG
z3&A*)TTr@DmY_^T*%W;s2W0}v0+b^;{u647^6Cww7>2Z8;fNbV=>k0#WdX`Wl+NmK
zY&2ThLN9yI=(DsXa~ZRiHVIpOq@n`7!@Al48i4=Oje|?uv3|yHOWU<BJycPFo@ZTn
z6{wxZe~10(y+HdD|D_p&myHif1UC%*$iBolx$Kdkeg`TlhNAiHuNd*mA7SantmPrh
zX{^t!70`<tH+Bj9?DARxGX#Eid5Guc=aoIb4DMCXd6%J!NAxb}<+qVP1$x<S=s$yg
z{5JGzNJpupd&s`V9W(^vFc`k+;`B}+lz@)PSr@0LfW8-W|3JrpUU(b&T+p|GjvlA3
zdo(Q);2N%g;k|&9J6Ctj(e-)Nl~hBJGEtLQwBvZd_Q>ar7AvCINaLXuF=6ff!s61A
zXp=Z&)`}*qsj+p%^7xx26%}x#Y2B*3#=tBO{R(7qF)_F}U0pm9JL<UbGh@`s=&HxD
z`00Sa?6J=q#>#ub%JC3}IIN4CnE~-K<Lt_iu;bujgjg5nz6@@uQNsvnz`4Fz>ICky
z2%(hIu`;(m{k-w8aStmnW*H%Dweg(M1n@z?TZ|J%3_EQ6X2j^uG9&!imU`2(M$czk
zI#y%VSle#Qe70_rrg~8I*v{2Mk-9F`on0`XH0@6OXNS+6FDfeLad`Uk#_P{E3MvCW
zgyTz$@1Bj;n;tWQR*9nds<uI{FDoh{5Vbwm7`3W9ZiTy6h4kO@=B@eTd78mn2h&g#
z3fW?hQG>}ARs2R4e~`uJsd)qZT%%rI2wP|*<<$zf0z3y$NZ1rUsc|iJac_~xY@snL
zFNT#G>&4%LdF{{xe#>hY6pOcWNxYeE7#+;E%xz3CLp)AQSg3||ahIJiV``xLwgac9
zfPNDboxjsDpvU4}Zv$>cbjl~KD)69|fFB3G_<$hzy@+qbtO^KxS@Dk>f8~W3cdm|X
zFa=)fb8G5qgYq?0tpry1hcbJM7mU%X?+Wkt`mG7?2Dx5*w!CD#y1F|XS@hfLGd82&
z+MDcrQPR5URqC%2Esk*1YJ<9|_dA&J^^0`_txqAw*o~V#wEU@~qM;QJ<%gCoynphG
zmN=AGTIz@<2ER-#f$#-OUh3oA-3sMzulY4l{MC%RH`lA-{{io+F?@4OO<LE9&IM_F
zsm$Blt=eb(bxk%BzIa{Od?21yJrohGerO!hDu{5bf{tz3So^KoHf4v_JFD%+xQ!R;
zyJovZY-BH1;gpfM*+%zGE%E8D=)G;-7vs}HarL&$Sp(l!ibH#_It-sw%vK4kCBC~9
z4@zL;s^E)D@yCy{R~f#!6wgUwcUNKM#l`*D{d#a39>DmO%UDXF?bl-a6D*G5tB%=r
z0Be8`M#WDIU>?SPG{1V14Q1EN2L`fXdYifCrhzO4=<fqrDn1sONrPC5-ePX?=0WTU
z9Uq{I?L*kXSMgG>_)u+qkIu@>f-pVATot0PW0xiTa46uPCCtL~b(m|-^<nyJdUa9O
zHOGd77_N{85kQ*EGvA5OJHlE;gx&y#8bqS_+3X>U(I~ifZlqqnKL17|uC;7SngrJ|
z1UpeffvPr{SNwUTemDT>GuPGE_whRSZ=g5CCz@hU1AVdX;@@ZZ1^GUKwiWZu@W%R`
zdcr(2p|O4s56sh8kHFU%a}6&C%`4v5SYKvmYl;(E=>1o+;pT=n_1&yRalZq4kJ`*r
zT=1!WNDqIu$7Ia8tV?@Mt0rLW1s1=Ij~R2;7y8}0JHy=b1yYunVU9bdx77P*6t6p`
zBdKW_<{!uP=DcijLXYL;Lnril`mhYM)t7p>^944_L?6vCmz~gC=(!o@>nHU7ysUdt
zkJI;Lm<cEKF8aO<bID0)quhB?4-G%3hUx~H8&r2I>-D0h1xz%{PU=mXeXjW3PG}lr
zk^*h?^aHF9i_^57?Ar~RcF~OcO0O4o)!&@t&YNhC`bv+DT&BjR<!uMcQ=~`Cm%c)r
z1ryEVU+E3&Y*OSls}a%*EFVSTM6<>zy}j$859?L-#YD6JDZL5$#H>?#Gj!DLr_c+u
zOtTcSVWDcUAKN9uPqqETQnXiDlT5STX?>9H&NOGA);sAVGtH9IXjGKHo<<)?%`~IF
z)*FNkSL)S^EcH}vJzJMq-2ZDmg6RdBX4V<Kjs97t`SKY(w)Sx)FQVScs%*si%`7{E
z@UEK9GQED-O(ifurjvI^ydGvHlgzGVdQ<(=N#@uxy-AJJ>JxQMF-h%UPV?*YdbsYK
zWF9Qj8#j+s>_dY@Sm@3Cg&Mzj6S3C35xMVXlgt`t^~UUebHcZJXefWFW(i?B+{qR|
z=>AFO$g_HL*2_HdEvo+LBy%5hczO1$9;x@4WLBI-hZ;D^yz`tMRp+AeeGI;}{j3NU
zhE6gEo<q_XPBQ15(;G(SDv``Cnl`0aW{y{M8{0a`{OAuoDsHNh$gQPmfsg9{flvW!
zzZC;xgr>d1<|uYqXZYB`CQUN0UDs>h{|0+irEfQDrbI>|kqockQ1KK!i~HEEI94Ak
zX0MxA$(lpHK?^s?GN1cKkFC=*%gs}v`LtA<Ec56$di~}PE0#CA{H-3pM0T-+li7VZ
z%M3rSH|~_8<mtv87_Yb9$1C(3#?e<#@fr+JB7V&tH`$zV9?jlmGH><}^Ud?Rt7h-X
z(iye=m@M<mc|A-=`8(V-sN?Hhk-mf7=Ly>aUex61P_$CDLkyM1@qkY%p_7Sl{&
zmif`QXw;%C^E%A%GW-Ji!=5a&{RO>6o41wOvf7$f^=+2$I<IkcG~vWinzoBQrfT8S
zO+L&r%?q%2HOovxbDE_WbQim6p1pv|)tD?g20LqpU(};(TqfhP1Dwth9U$cmrYpv7
zwp2m919W)9w*xd%uB{I6Cab#{ceoKUueI&GDbZ`zjD^}RbJs<^iRaQ7ymV$s5F0y|
zx9<(i7ge8f3&Hd{n5ySYdy`i%@>NZ{YaEAy0ZkH60H7U!_Og4oXxgc#rB-vHt7*<c
zcE3;^2<A6Jp%0*7JVO+)1VDEKa=yrdgswTLT(4Vev|xsTvYbC^XuH|n=E8Elsi&U`
zNmOkHI5N^4MK0ix&VBFzMS4is;gS9&ARg&G=@jY9D$<7DNW+9Gk90J|3K)-cnSgkt
zuX+H=NJ$rsFD~bZ0@eT`-w9zJdG!gLv72pDv8^;TEgt51-@^>GjUBjsU4BwzzQ?HA
z{P4G^F5O}8F`4$601Fte%pO28?Y$~j=DbyzD?))+#+AWUd1cxFl6C3zw8*+=jNtLu
za{=+#Lq@tepz+wNAQBOK9R3atSdnnU*$lOaQzCvX(npcCx~Rp1AFWy(7Cv~3XqXUX
zCE5c*B|LC^=O{P#r}=cAeL|Nv$vC0GYq3&5ycT-^$y$)Ae~bJeWO=P>W^&WJ*;Joe
zouXRBZ)o9Lt(E@ODq=Sjsd|Z~nX@kIwQ4p|Iqv|9tiL(!vL5cKJ<1J7n%1|3=+x>(
zY=q*?@Z-h%WVMJbR#F@Mq<nL~i@l(vreO-e0=D3J6eT5^u3a#HyR5hH?1J_k;Y{=2
z+{V@f9mES1{8zvruk}P?wzD`i)wiIPup^*3VGA#3HkXvJ;})v*cW-5fUedG+%5xfT
z-pDLr=~K|(f~6hos^oUiemE`@U<=}(fs0tkXg9a29fg8OSe8gkppAb*=NE`rVZ|*l
zPhZhndG1o;o@JIeH-C>JKZ0rtD`;U<e?PB&)wl7mi2qiwE05z52lLqq_AqWni>H8j
zknuX35e!??+`u)w36mFF(+x@^@m~SsCcZ^Iwt3j|N(^3bzu;eA<+V;2>I_5vJ0;ng
zB~X9YGW7T)UhjR%&dZ=v#ILfq6?ML<2i9h(e69EQ3wR*`o`*Z5eV}a3c_?3eSJPgy
zs<sctH?vHw+K6aF(;mm&zli;*jCOg)awVhZTfq1QEXkHt%1f;3`91?h?Pm+pvSA@|
zj3mwV_?zHc#gs5ie6CFJ{RxUnCK@Xf&yxiaT2DniMO3;G^DS-40_+}Ec0cj>8{Wev
zDiigP2cAg2oR!t%Ti?h=020}tm#iS)CZXsQ-D+iIs4`+L2>5n4C5t$YdsHBNdxavC
ziM`6iCVvz3)CK1-aY|X>n=)>(u%4Y&7H+LY6Q2kIn(=C+$!hSe9Jg3NhmX(GH2nr5
z2s2w<)$4h}{sRe@bGfD!vojbBSCtF?B8gvY&Ha#g<TKwc^b1<1kj?xV!?@MpFuY2_
zIC-13Mz5>Z>DlM79-YeZi;=)PbNmC~a`pGKvYuNPx~&iY2i>)RHQfk5(==@pivYBP
zg^qP|&zg5B$@ih~izkTn;Kr73(KM1P6UP2Pn48$!pl@Qk@pl(%?qf&!Cz30K`|tCf
z)>`?$6aEW?|9*gWum^l}dA6}iMuGAChv~Vd*YObFEeySXh<lxc=nL3r{9OhcFSF7j
zP4krtoKV2N2aR)IW7=z))>V1QffcOir^8GE3j^a_5l%3mKZWq`sL-QIxCFwzuz39S
zJub}`kM8Nbx?fF@N%C1d_ERcXp#qHpb_Br<g>(UX8_)}YUdJ=wpzC^_J3q0=zX00J
z%G0Ge{ydU~z(f|XdfLVu6U>6^dQ_VXW!=e}1{Z3uOVofrlY~Na7n`O;)H6xOSiU3;
z_YgTY6wT1ac9;jm?*tv8w1dr6O5Pg*-#aAlA!`&bk!!ZQ&(a($G+(N$X(P?UikB0L
zPFJn{WJ`5{Q2n%$YK~5-_Hk}e!9wn|Q=z4~Lul4kTnblr#&JbH>k{j#P#o|+?<Lh_
zr`OP3KZSk)<AbqnHdkr-B37EpnmbI>96f@eQo!1Qz7i4ztTUjm0Aa?BQG$&iqIxK(
zr2ddab;*9q_d1Dd8+*^AdJ3uGnfB^7yr@?et4Tqv+I{)Lec*OX12eh%3P?2sw2kE`
zS#F#^L8|e!NRe|8eYN{93w1ucdq9bX_lx5a-NJ)+eDhsb{DQS422eg@zo?AW5<G4&
z{sSIG0SrzQ2J-;zU`6B1$9~bHJbRRUp>S?|c<Y#44GlRaedTCOa;_Yc`;?Wdu!0`_
zAk4@yc@`pAN@Glh;U>G6@iF-yB^-zlqj{JZWgm$k`M6xhZ^rUz|L8cg#tpqz_gx^#
z!FLJJJ1jk%<9t}g<PfJTpxrD%H4-0x%@G7O);2a`f;sAj9(DJ5<=hD`{Ngi+h?bAP
z5jlMPT~Q)x{GB)V-_RpGzbI+}oLB37MCkD0w+K2rSipEUjA_<7%%?vSt4t_PIimVA
zYQVem&KzC?-kodCqVC*S8H<bH-MIxc3YdrY=W!5`{dprG^jw#c><k&IYY`*;;TBc(
zXTFI}T-(_8+jr;T%3w7qr&a6nx^N%p*^XP1aL@bm_kgys2};=N&->U!MfzUpKlm{|
zfU&2wY7ZKilU=j<@;4pTiWl_ZSz5IOew!Qz##Zs{{EjkE2%&0yp7@AYxYcIBD<D$9
z8UZYxZ4SJN_X7dx=8re^y7%%9gWk`dyaUwyizhm9RZBi<SvYAK_zTeM4C=dn)9dTu
z!c4uE;r_ZQosW|D*g539JD}GYbcAjz6#}1KR0ys7&7{6&nR(hWvqYGIx-ipCnlUf`
zrdMs)O&KWT8R#2|M9B=iVk!SU#|-}+33yzn)iHbjt_OOGRpLZ9kn24$$82L8ZsqVm
z`7675%1>>XDy`a(Rl*madg{&fPCd^lIq%HFcS7<Du98?=-f{yrHq|r|bVVZYA>WnK
zRI*9j-sLp%T)O$!?|P`GKp9PgEQa6~2+IkhxA4p-3^`eBSAu-9;1}kKhC-On74smh
z=8D&pTrLb+ckM0U%ZZ~5(Cf^lI(?k^)OG!iJL}FPb_YPaS*XwKu@p{JD$&zhPY^GS
zh96O0N@Un4MI?OYcoGI)VP7k@m^pT_!wUVeGmpH8Jv-LT0nPb|HCp-5u?1msZxMS}
zaeJ7L{;7v}G;fl{Hvy>%`F?n(QN*q&>6}|_1={$AUVib$3AFJ&z6$H3JBOl@e<(%%
z<g}iYz5S(#$s0EYRk4=~hjf?v=oTu?ir?Hqov12q|5U`0!uZp4lK9JC<QFmOd-arn
z&m&(Eo5Guc59)S`%lmN=Tc}#6w-0qUTdk=4$@MA8;p<LRv_NtA!|N>KOyry*ADtJ7
zBXnN3bl`tCv4ql5AEl3pBb1I=O4zke9HDd7N5}S=)#wP}rjJe(aVGNMY#FFM-z30m
z)Q(tiw1ZWZEbS3<=HGh6oi2saWXIM_hWN|)PtKD$=#z5l*f-CN{avqZe}y&5Fi-!D
zFIh<@tI+F4Jo_FC!*tFLu|_B=DvFpbMYz{$nDG^QBi7RFSD{C;Rpz7$y}rl)GbLLy
z>2nN!{8zwAF$%83V-d@oEKRHe9qo?Sz%G2r$W~NNk<Z`j*g@z#fY&CpovjhCIWD&5
zBk11$1=&rS&nK;ei4dQw%wTaSE=SMCmRKRgFv^rfEzrsRHsr=F<Rt~*t8CNcRU>S1
zt|5?{jsFVR&#-@C{;HR3ZGv9{fgK}an^l)=5!EhPLdF8q#cVNbz8S}CcX=8^khy-r
zUF8F|c&er)4d4{*L-r<CJ4pk(bA-8y1b-o7Od7;VIY-!pcY&mFq+;fWtXmCWLpb(#
z#fPj3l#+%Dp63W#dl=X-p_ViA2-{l+=!kU?^|%grVE7}pwm(Qu3Hgc&XcRXCK2mDr
z<V-ojE>{I%)F(pcbEq6;e6vK-=<m3Dd3jv#gK>l<pdFLOa7GSYu(+!|sN-sW$(`Pb
zo6%8rcr7U7SqpQ3ZVPSJ0vr<0@qEB04+l3R8C0i>-~Ep;$0MLj9B!`FZB1m>K49hn
zyz#uVF?=C=(~uBbQ;x6?Sik=OX(l)yvd%5>lBo@^j%NRmT{sBr9)USOWEoEaYb!9<
zN9=|Az}g86FX%pE=i32qFS)w*K6`k&rVYY3(5^Py8|?2qyR9BuXV$md&ae-P^#EJz
zjr#iw%zL)l>Y*&|x7GGJ!<ksK<#t;fj&+)8+ifjySkqj}%e3s`{o8H#F&s!We|y;$
zUB$MLzg9YDM(?mi;6q%89k!S%LqJ>0elyc{*y`x|6m#+p+x_~%g=X;%TgRHegSnoq
z{T}ZZ+Vn&B>%Yjf@3eK)T~o}?J8jXPy^Az0kBwT74PVc2C<{>a#hSK}EdZpiNZiEg
z02naEjl|&7ZVZz7?MjU0Jk|m9=OCKL?gI3MfMNkfFX68xHob7GmNS**u@69;26i5M
z6A)Bt0W1W_bvCmMK*Oh)*LK=MA`6AsGoV}$P&S~*rF^khJjHDIiY?aj9~F~RM21e#
z3?80<SRaVxu}1+-0)%TT`>_CPKtK5^SE9CaK__C~E>!-8)KNfN*>BH#<?Y-AF6YLd
zM7^Cbd~K<AKU8Z1ybsf;wu*HFn4vsa8k<&&fMhUfToAT*z;?kDlNH*+hc5yxkG%)V
z%bdHFeFku&GD8)sQpOj@dF&i`9}9&mfUa|eZLEtDK0X*$A+`q7U!aX|Mi<vMwswk{
zQ)p{c`%@*je?J)8Sxg@GmzetsZ4Eqim-C0Pr-x~drmJQkX+U?dr(0`|<{bS(qVbxe
z1xN46<=o*}VBfj;Wqe&xIQVzv0$&e&E31yRvOZPDK-qA#m;!ztivsUmNMXip%tg#6
zcG)htOO-U=c*HNZ#<K&D;Cr~nAJDr9jXV|rTK5$kssZSAKwDYY4UnE@&U)1rUib5r
zoYx1GUjc1qE!GL%u2*d#p6&+cbpvJ4vmELGXf+^o*3Z!e9Mjy~uU3LC?#xA3(qZKj
zA$bGrJ63U>D}dqv?P6T)RrajE;;}sTHfU0BkEQqWdM(&TaSMx_X7(+@`0J>=UWGYF
zRnIHj$xJwY8ctqfZ42=F#Ye&!qvx>)LFDXJ>^y4lkYWcSbbdKs%VU3mzFfrn6QCj?
zyYhCj`#gDUtrE#&#ehx-wQYbHI*7@hx{X?N0HVmU)RtIkIZAB<dqhRGo_98W$!6-K
z$XMNq%i*b-S;(tA099@XhqCQ@0FrG-3bC=M*>ZMK<R=ks@`NyN%qLB*#v5}LAk~<U
zDMyK1j~er+kmQXSw3_RnF?(ot0o=v7oNUcoG*xR#Mc$fH5v}Q?w}p+HX13pr#$2kr
z=D?bx>fl}632)31@F5%1SHfyc&X$e&f-(?*itvjXa|tB!7;nr%A<G-HlI*_zvN5@c
zY|LW<;*I$?AlaDyYF1-PHQty~jW;G&+rSE@nJxF&V%pSuj=zp+j)C`>fbIcQAfV2G
znmo_By#WmdB&UG4F)a7k!rPdjDAMU^ZeBCZniYaCi#+jJL45bw%GLvPOgBG<<{hV{
z=HnoBo^IaQW2@I}%yemilH>>-@G7tI(~!|H3&T<OHF*Cz!*mzh!rL_fX*GKslr1b}
zhE%;cSWJ|7bIi}d@C(yslJ2~R_kw$0#16(8=89rlw2pFDu`Sf|gi?&Fho`?s_4kUn
z!>=(d@w?1w_Bo)<Y%x|VBa~n|1gpKxx(LC1FkXP*26j-fZy*44i5x^VT9FE*Mjm?(
z#7`id$Cd*M#DBY+v0Nq9L`2J(T<ih-C&h*W+QwGTFemJ_h3c=)Fc<H&Md_Xy=C-}I
zh!*p~-M|(o8#}td$Oh(7=o=j6SsRb8deOWJ1<zI}m}2_5sHCEW5V|ByL0m-k+G|Rr
zHi@Vy=wnF8Ui*`f<-N9&>^|9PxrFSrVMtcDJjQ!%XMnQT`m0&JR_gIy%k}aY@3mZQ
z1G`}!-Dj)o3CMPfu!2=^j;g(Pb3r~XhY9<9lJ%FeYRN_T^6JJ^HwUz#Q0BK-xe2JT
z>df!%*R*WnI7=blR&mUs`)zd-rufM1zC~`K;(SKbEN(u<wxW=KvwCBO^vdtqYuS+*
z=E428+RZ;#)XzYzu}8&|#p`@darm^0drY>8Q+Ma|XV#kEyl(5FA75)WdBfI_wJYxR
zhHbr#H7x$&E!&g%`q*`5??Lu@C^x4MvfFgmY%^nsJx-6EZEhT5Z=uJ{Ha{6+@296{
zn=OXgqn(?PW0M)#X3|i5xL%NLjs-Q^k&opd+p5xi8iQ!H%G2}g)ogR?Q2YJ*@O<;q
zP<zLk>ma*UtV&s%bsMp#Bj0R0%-#{dnPN^HX0P9GKmmU;!ippyruQ`5=Yg*0E0TbE
z9F15n(0DDNK1Zhul&?qv8uZ(Q*<9h^H+&NKMc^CQEc`uZem~3};S564W2-@h=^Ad2
z#okjhez?7UwYxxF&T8Xt@@zA0xV^4EVuLwlxV^u}|Ce{zng%c7wu&*I%ie*5MbOD*
z9|0-{w3)3zevZy`!wKjO(7BvFvGUGkyFqV*yyUX&fL;N#h0T55D_dP?VD4Ma^5Gzt
zZ3WZ4F@}R{0Zv7xa+v{W7a$DSyI|(jOmqAQdrZwQLE6Ty!rgZEJvuGu1taXyx6s#1
zdMW7ommAHiBkZAi6^|M8l)b)ZB*NdynqfMx;&F>W5*s3vRF=yg0r3N{SF+B4a#<Vv
z<>EWojn(LZO1yxJw_xR%g>zX1L>i+GD_M1+!ti%I@SSX3j<kLkW_@&!6nzFSlFaQ-
z*{kZCW|@1RvPX4%1ERP#vjjvOqa2jML6s{nY~l6?gML=z=rKSJ_}suYAeb1B=^AMd
z54lrB@j57*0ij(F18n9oJCC$SRzJc?Q0NQV+#GW>SNK^dj0Ytg4#1oNa8ZuA2Fy11
z37@M#8OoWju%A@x?jD3k7ql(!={}3@oorT)`Q1o+OHU6KyQ>lJ2JP_tkjuVAt1W~b
zd@seM=PTie(R0~45IK7}yAL;ke#!tX?;Nk?vc9126e0Bj^r?`&j9yVmc3<OMb^{D9
zlFM8Wz9}FZpjd>xf|Uq0^UbI2A(2w*TTrCb2|!$G9edJaUVYjg6WbNoF7_>&ezbtD
z0eTUT99@|nv(+eM>nBh^dIDo%xr)9DDEta=?GPy#1zsp1h<$`{`Lf5%hGw<<B#95z
zZ_F1*+3PoIta^)GsKY>=T$T&=5H1Pr!^q`D^Se>@h^h-fS|L(;(R7WrhqmJcnZ#?e
z+(JJ{{6&9JNtD6xBz_`fc@nRggGSp!Je(j?_N#z+%3|<%B~#{aFd_gA_c64k8qc9r
z<2mGN>(~tyLp1umxSnJ0s9w|<Gl`fhF^9jx;(5+OAkq$!Zero4@*L~pF{_QSM}{6y
zwlaw;6yKX=y2sc<U1t<G1^PRgZMHdbjQt$zVm>_9-c`@tZsv})H`IUJUc6_ly-R;K
zruet_?CDkN#$~`QF57X9>AJ?47RWbndyd-g<Quqm9JT+-_i4ZMvAup6ooldNfSm5Y
zAH3-a<Qufler#V5M#ma#e{psIv9alSR;j&C7@caceXFsadQD<upYtoF_6WX_`ADgK
z0^gAQz$f;tc*}9=6MKoLVXfYeKRB;f!pM<!r)`Jsd|AK0mc5$ovux)&5MJhldpP0q
z?DNi6plsumHk@)S`#kWi9B;?*<JsqJLC*7E+NP~@PO@d{Et#;0eZ#g-v-jJEJ6l4i
zfwLy4FX=j7z;WuZ;m%IB_x0;OB2mr}5P|1hMeGE<SZ6=IX${*XXQpkL$C<lg*W6jN
z`rDc9H?1GIhG!k%?=|N!fppSzyYA>G{yzGg^QB_XXZ9hs`wir?I@EjGgKbUEL-GU7
zIUDdEPfu?)=~b!ieT~KHk(9-p$;B~W+Gjbt(Mb+lFtvf5x2U>3+gzXqawSfsQysP%
zzSOjA^ZXfmU%Y!N?o?)<8Pv@++l>x_E>hjsWbAA=Cu+ix=nLFr+-x@|YQm(GI|=|b
zJ-hht@9bF&H{AtUSM0@ER|4vGVDFmj0Y|90sj=f%*3kSU#SvP(yNP3fgS8io;^r+K
z{TcgMB5PVX(jBbayrZon0>$4QMZpaeKiSstqMiL}GEPImjXv7m;WY2-=qSMf5VNAA
zBb?O}9?Zf{V07r@=*&+`2u3yO8jK2lN(Aan&Exkwn&NB-;73>+so3cOhs&JufTNkD
zhnmYc;R8-+ASKE;A<_*(3u$to+wqIeR?7;E>FRhqFsf{>%ujyUU<q_mWsqholV-Ae
zIcgT)(cAGT!|{+OA3>ciOU75d9j@Z_k2sQ86m12|AE4A{c+LDQrO!MH>2`^b=EqW=
zNpghqa@YNi5c7c~c=$L89{5?4*krI!3Rcx*Xbnw<Ry1vd%6~)IKIQd+qa#9mbP6o}
zk-{zEWJ|BUu!NGciamWDjcgs_=c)SE9VUT(sziM)-ULNk=F2m`?hkA6%2E0NhqE~B
zamP4kw57ILiSN2aO~m!a5XY3LYIBwE+X%`=HzDUNo-owWLXS#R6JbA9GA)(B{y5mu
z(-9e2StP|0F$X>6h%&2>aCESagN>r8Rh#6>AES^=U!vxB<tpY*RU!rgZ&;)O1pcZ>
zb`9Ii@(JtPWTLD#(+2ca;2?{XN?_hZ^V)Ncu>Rz7yJbIxz}*%pkw7xYU8dv{MmCw2
zLLkhiJ>`h>>{B*=Qe}W<IFk*h6|4$Ok-gfo3X1HPC5?>!Yz5etKr+m5kQv;&xGIIo
zw7&iP5~6Ey`Q#KQd<NHYghf+E$cwJE_Nl4q@;Jh~lMJN_H7Af#g$M{FOa56i*$a-a
z`eciw9NP6$Ep>0$gC?4<jB+%t^O;rKL^|?xF~gh??FfsaBTgs0B9w@ri9)0gSuar{
z<b_noB+>mdMfBM!gJ@tx6<wd~Qj$rI^p#vtn9^$?`<*X7c$6bD@S=+As=0HNBgTAV
zoWrRvlRlzoVX;x!z_1{26BZX3nqBD#Wr<gWz#|qZnZVC2QVM}Mg!K=GQpnIpmR4T^
zNt<a&1i}#S`SE2+y4=!ACGZDDvT2jj^9i{j0~n`e%D!Jsm8dL%H7!yif%Pp?GJ%a1
zNmwVOi6Y9mg0Q<Y&4V=@VRw_y5X)K$fg~#XE@4pQeb?hw3J+U~sRRzNNXY~avq-4~
z((yDI7h$U`!DIs0TBKA0*JX+(>rO$WSRyF|l0g|gVPpWKmO$9zt(H*P)?j~I4(%t~
zFEyIwSGG3T$6E7oj_9cWClUXcg2>9VnP&CzjtC2%FB`!JF9yFR^AxB=?6MM)O5k@E
zDVe}WC&><H)3Va@3E8QLvRw$<uc%B*Ch(lEwnXsIc5#wunxM~=yW^SWiZ+h0uH^Gq
zD~v<}ot14LvPg*po|_~LMO{?}NSn|1igZ)@3i3v1CZ#yS13y*P&Y9~cAd7^a!ihwT
zf!gGsP9|c3l54g$9HF)ERf_kk5}cKS6Gfc+ynzykGl{&$L8OGkh=Y4O!hLa(5lVO;
zuL%NYTBN>n6i5w$TNx!Dsphw>)?{d`64kX1`TlC0X6c!oCpscLq()AIv~>K>kr3&q
zpHH;D%48KSOG_UK1`)`YgEO@HmcQjnf#~G(f*&346I(vxs2(*?(TA!M{fkskJ)Wm1
zVv);Rswb*aNd#8B2@qBSEnH{>S*U8PNy`81N#dsIp`R*ZS#pSAE+(*(MM@^HNRh-q
zBBbtQsVWB&Va*j4Arja}k%YEHZ2010b~YVRfen=DrY6gBgx4mIbj}k!GEeqMqKxr|
z=XD@+K!-gC2wOouY(7kax~Mq%OcqNW3-2{qjuWy>hd90CWUcZ^$2F~S5=or(R#b@u
zZje!s4$f%uVG?NJg58y}xyNgOzzY_sFP*uhEMt^VRv()zMoByw`a)^(Zsf1VX})?>
z9pN6*Bd@;0glPMfO>3ACft=BBAcwyy1)`JFzN!g`PQ&DFCI6u+(R)Y*_2a`NnpouW
z*7}))QI{xx8JMAn>Pj;~J~&tL#l@_4T1D_(mKZ2z>@<h-{`W>$qlCZ@M-nNSz>gJ4
z4hX_FjuP<;L&HcL_h$l+Scdu%NUHMQPFUV3^4OQa=PgnSfpEtsO%Iume~Pq*64+Id
zbZtdHGX9Gd#n9?cOO3{=WZF;y+ba^+{6N8fc*QjD<q_alDsSYOX+sF~2_ms_5JM^l
zG2}KuBvuY$h>St{7x1M@K@6iH7O7gJJ_(hB7)D8<V7M5;6b;jc5jf79orKn_6+#hm
zU4p$S8%AK>DD#zRj_}5Gu6dYhqabaald|enHxJKn)Q|R;<&S(mGJMX9@|VG*vyTXW
zP7XS5HKZ`>i*plQh+k$kbP|CSy-k}@h$;Qni}ZecTC7q$gp)F!EC~EjnZB(ZIzYYA
z>d8Y^`z06s>Q>*_R;Q%jrOti_9(KKY;tmkyFNEy)i16&(wdKFN@`?znT?jAqCX7JJ
z{cS5sdL>T_JE??|G7?5$$-gCxPHDIG>D}#>!Z>s1G)F{qWs$giOGLI;l)ns;<0Hb8
z6PR#|)gTL%7P=6g=uHHHG*)h#2+|uQ?ery_l#vJm2mM<j=%BanFpp7FW}+N$sU)-7
zN>+TnSH>;6m9E9-`}Z$@J!HqLCSE`VHNDl9$BX~7S(1sfX7(fk6BX&UDI-0%61q*u
z9*u;^BJxNe>D|_Np)z{e3M!Sr*wL~sxYH6?Uy<Z(kg!IhYgv+wl*B!w#j}*>b7i-t
zrHM1yl`>A}q93(frV^NKk&+3VW7*=b5D8oF6(n$j45BNYmfmZLq!Re1MM@^{phZe0
z@NHEaemt6AH@qSQP8cIuSiTXMp-6UZ!fT#<!sslq%pqZCyy^s&dDRO%`5bd_-)G0j
z%u*~rS+S%Lm~I(MCXl@2VSqq#EtevMeeKmHkOF}&fe1ppDW$W%!DFfRsRYtl6?7j0
zKTDTVc5TK+Pd;H^Sk&w-o_xZ-w5T(YJi<;`RL@>d{t1sl>8Od!DPh$oP*BMPwpAn<
z0b#jbK?3Q#gtSQ5OI|?%@}&r!l8};w9rQ{P_?6|GcU!`~_6ibMc8efk-*^QH{5JU(
z$!`@_?v*5vPFTnw3H!q<NMKd@)xbjyTSYZ2e%VPQ*fU<{RmZt~f+^=V{*Ijx$~qS$
zm~t)z2|1<$;erAVM`9^MLXr@~_sD%%{G^PNG0^cz47+S)Wa$L4Opba`g-NIIATmyg
zWJ{6SL8^9i4jd1E-^kY|@rpyIufh34KCBYweyhtT68OMa>6w)1kTN`bQzAs7No16j
z-wye{GI)Mz<S!G17ZzFSEeFY7G^3tz#01gF?tVBPjk|X=9ozmF^HKgXxXqE(>J#K`
z&IcDyR6e*^UuLN2r;Zg*R6$hX<zvmxiye`z=&bcyO3tnYxonkYe>zuv(Oj|^uMOy|
zbhx=}xg)Z7xN2!?akK=TM9zO*wI}r-UmlGc8Wn;kzVdY?lhMjF-op_{UeILz#zFGh
z6s;K7t8y>A@1J{}>X<>Kl`_ucUlD@n)N~GTJl|QkXPn4tops6!on}T&LKcXc1QBV^
zICJV!M?@SQOr|5qs2PEDkQdEN!)}~x=4jG7r5y3+Zyy=n^pBAYne%;@K(sfNeidzo
zeC0+2GVI_CymBMsJtmlAmN}w>XcXg&v-p~IR2d?@+XNKlFN5YGv#6Iyq%Y222>9YW
z(RBRzvP^Uko!_fE-dw)i5z)J<YD=0)ka#*YoL`{29`zDm9!*x%T?9{dJu(@rOrwDa
zBrm9lzj2VfHbvc{{-vL(*)51f9vm<G869Vy0GwB!@bvMbpFKEF#guEs)|bHN6bbM7
zJo$utqVB$k&O5tH({%pUo-Rzs-BkK?B-yS_Jgilk2)t5e=X>!!P1(o!;!ipvUQ!0m
z-{-{98D%<T45P;@8}&*Ftss)dF=H`@`QY5@-E`u(kG0g8Sf=F@*59gg3W0;t&97Ht
z7?UD!Zm13&>wHeBpcgci!#0{o2Na9<Z{HO1&Pt1>*QB{<I@9|vI#K>ISbO=H6)(Yj
za51I(;DS%bbWh8ci*l(?oDHq7w5bH*_~gL_j;M%A5-FC1*kI6ww*QQ`ULdCuIK?6*
z6F9>nr4pDeUMP2=lab545(E}lq+|j&TBKA0HwimY|6*s2e3xhrdd%^S*cj1;erx|N
zuQ>uwTcl(H&swBZ0?!GjUFo;_Z+ayN{L3OGTi_JAzTq$L6{t@sUT->TF^@|T?(@nJ
zNC&JiBnceql_GG7Me0l7az&C0M8XV9FqObhr-%`51<~0ONg<Ft<4!~1D6bTO%Pi6W
z0?928A|;=&)s|oifp1P}Oz9c$rec5O6(bO_@*zH+5=H^y;g~=Qkh@IECyY)ILXbd;
zMhX(P%5s`SAe?q~IU-{z*;HSt__7tPkc#t{!rNmKp^71m1b&2*c=&?dwX#+cY3;Eh
zN+qoiEm9JJFvPzfMpZVHOokxBCxocsN}K9|IfX>t{6`U=N{A-7jTN7WMJz#dL{e2U
zig1b<jX|IaKZs6AW=#>d8DE^=;qE5>5-S7$(BXZ0jJP$DVumtIhb7@?#9FMk3zWf^
z<!G?*jY^(QMH1f!FKC3VIP?#5{AdN1O5lrAscR<_xLJ|xTKbil`Gg&uDjLYL<**b}
z2qaNCpa|RJ6(kV0_!q4IrL8EJa#PE4mq_4CQ$?>ZOV>M`bsBj21a3q)e1#&8TxvK0
zGgY=Poziqn6<ww~ox}XaGLcH)OVh1161Yu~cs5e=39FYa(%8SUtrW7=dK%g4OW<EK
zBuQqLFfH2>BoOAso07`rlgRv8%lrTW|Mm$rDW9+a<rFglfk8fTB;^wp?iEy^N27x{
zavmb=ORpq>H*YaX7}W>939O-ND=iXM+seTJ0_$bF%`Pr%U-(P46AqM$@kw}9l>xi9
zBG;2oScH|AR05$mej|41ks>t*?so*zX-?TPgi&G<K7o`~89rf@*niEjN5fS)$r_!Q
z`>#2sBVtmJa!lO_z6qpUqk;tTYz}_P5fMvCpu@n(9$pj@=9(tPj3q@`y*0;>;x~c@
zE|u#LRaU|Ko-8MJO%NesT3!o`NRo?~ChkO0UYvTqe_I80QzjpgZ5Kp`fO}6fYYfLi
zm+;&03x?t0kob`2p$1Wt;vg%(yWdUt0Obc2EcE0PHdaxwkRlL@u|?=aq(~j&f73Vq
zue(MoE=hw9z6t!Nz9DKGLqR{L>Sfbbxbq1atlFG}C<(Ww1X7~S2P+d#tCF|*4BqBM
z7%SVH<i@Ht_u^KY`=vH0PnpDNT%PJ6I>?)^+MMv)@I{*wAM#?I@}9C6Zq4d_e^g-5
z3^A9)Q7dk>L<Z0lPh%TNra^mfhDhhwj!L?-r8R=UO^T#zlQn|LBC?+2Ii&>7S-MXV
zNWpPkAD8^z@9Cw)eQYiCc!gw6;b5VdWq2o{Yl708WjP#NUrpA@);vC#q^W#c<xb_$
ze6NP1b`YIxKC1?mLz%o|x*TLf6+S`Yb==B+iYoCit34xLX~Yw0l49|P*(C}t_My@F
zP^6wec}CIcV>(leDb7~}_y1^O-QKh`<1Yg<Uhi@l;fwP}YhRo%Gw5*fe4I$e&g39E
znR;!8=)EmzME+vcW`&!wOC`zNlsbWg|6)eJj0K~G;#e%7#Aq_Fku7d$-sjD(^pBEw
ztovvxfy*pXGJz{CQYwKs>)LviBdUwyZ1GAExWghP6S&JFr4smRwtU{~UZQv(du0fm
zp-zS2p@6^~MZ#kB&zbpzt)3|tt7PkfWh;e15|uY@!sc5ZQV4`Cu~_{d+KRfO;`zaf
zE0I9YOfgkN-B7GQy&?o=BO>0BDimw{Ox2MTej!^<WA53q-*i(YZb}3WnJqV;CM#CI
z*+QjzfvW%8RvuCb9H!12qQ?{Xlp^sQrsfmooGquw%C=I-)|Zy~WCEeZSFVFTR4k&S
ztqHtK9hsC(M((!w1d)>e+$^`ouOd}S+4mH7a+X{klM}B+eiJ8)(=3aeN?49n|5O6K
zuDM`k*PQBgEd%2>N3(g6qh=HZLDloAsY*r+iF~G#B5uEg;25)*KS&6;SZJJ5!v8-j
zXXfv9I3EaC4Ti7_J%-BPT5lfxGbZG`2{YI3c0>qE9mx`MM0Q@Z>?Hf!fs~IG{$3%D
zj$;42AwR8nf31J2h>sqRIcSsB3q>7*=$!Ul9&^VY@O^OJWP=EAXZ3L2bqMRMsCZjL
zAazx_)S`Y!d63kTDhPBYu__`b4P_9m8!EHZPaODIG7+Ck#Q5>~Ggl?rFk^}wQ9*>_
z1UT=0fs1ZsTWF4~YNGNo*dwE$lh^}1ay?G?ZTMnsMtsPNTiFz4Z@R~PXs;tWdXK_m
zRTz8?MEMAtW%4yp^NUL6H&s^GR;>2mUv@~>evSLJSwud`znNpczt_<?C{xKz@tD#3
z&|Gv@J5S*+EB+pZ(^2f^Ia2<d%KVL4vd(lS`j=U<&V=8r%olYgKID0w#qGMd-s_H1
zR|`4`WvJKSrYh+ST%p2o=7|4Gwk_!F5uF{<)VUfui{w4oCFCPWzC_tMp-P;6;`Sl}
zlzxyZy=RsLzvD=C!eERlC#!P4Di^78i7J<?^4Uyj|6lY)ea<WSi>jp45_B}7L~T#U
z-gpr&{R={E6`hVF;LHK%cci`jw5@%rlAEDQI$%)fwJU}O{(&SB5AFY_J@*Th+&{=i
z_;2I)?%|j5`tRZQu_x;DPw_!c@Y|?e$M>qzXB)iWKdk7qMLkY^48SI7!KZ(4=XoVZ
z+tRVwobx-<CU4r^d{)VwS0!yVkM-IW`NF>x#PwUyMsV7EeO2lGgM5Vles1q(ZW-@~
zO8++YM1B4#KKK=U|81?MY8#+ql;rzw>?}}x$O(P_y^U_gPgi{Z-F>lYhgx4{|6PSg
zsC`x59c|J+{Tn1#m0Ylr!wxs@za#wv1;43{RJz*xK)YvX&k605s-^hYgCgRmy(rBU
zy_G7l??mu#Ll^Cd-6EXdf_9D2UW_vJjgfxt_n7){jjcgKp8f>^+A8v2<RiRWP29q-
zpS=kGZ6!YTMEn#N{YL*1Rlnt`eza$Veu@7L#m62Jp-;O-K2`Mds=TP`PuoWV)gFRi
zRn}Hzh$^uOgxl+Q8@i}hdnMONm9#U2KKQ?<@^_EolYAG2_fVyGw}*_MemZ!y;`{9G
z5bZ%5Eoi?6{jNA|V{)tZOIQAB^9Ajd_;2$5`zHR2e8hj@x1dcEzpMIa^r>H!PE}&(
zgUEM4WxB|BGey5kl@F?t_Cnw{+C_fo_uNVTZY9@VmEK(sQl9d*MDhK%G;CM&v~<}%
z^qcF46y3WiLGtM{|08N#(ar#+-;uuU(}#5WmF$m|{0UY1j0;hJv@5S)M1kVdcYYeT
z#P|Qk@6Z1if84%V<G7=Jmquz45vR&{Rk~H#U6mz|yG3yCs4`K}@#6y9KXANc{c6A@
zMbA;?d{xroC{Ha;3RFq|%p{$F{zB2|{P7KiYidz(w<_JLTqjGj@PwmY&jLlHe<G0%
zz|p^jcQsj7qEU)0)6|_XTa~mOVAEvt^$#6wBj}SaeetBvo%l#7TFC#4DDm;ZQt+92
z;YM%b>CHRdl?(oD=pxCHir-P`yH!aq(DD0{LZ5zTlHU5#TVA}E75v-Kh5d1gPrvK9
zQrY*utQGnl3lx2$Dqm6MQuR)-P?hv@kKWq7J45=X7lQx5Z$ZCLNN@G%=K}vhKEnS+
zp}#->H-J9&cs|T)ryOxXw2cGb<64|{?CC&1IdhMCMdMcGL#m{ondz_aG*ynwHY2`v
z<ULHknMTiS&D0YJJ-el=r=LPq(gWxo^(0cNN<0(tARnf81pgy@;vxNRdcdS7^#99X
zTF_JWKk54gA3<`r310L_`rWWzDoNf)P*M37gB88AO4g&Q>{mI7qD9x2p+)IMDW0Y=
z7hP!Xnk|QS9aYAvvau@TRM~#Ex%fjoQ_{;gym^Z_=NMi*_naezvGeA&Zya@OE=@BA
zt<(`sn3i`|<4Wz%YEg=W3;ju5t__}=plj}5;lZ|=Y77{`CDx_P<d7+YwvOOChn
z*v01CmmDowN%6Hyj;9#^*}%sy13zyry6k9$e>-X4Wk(PGLsI$`M^^ko_0u!(FRpz2
z8d*)PdZm9ku1L}GTie^NIPl9GndZqWj_%=fl46L`AL`V!dYVfsnP+zR-qBspnQYER
z8AjbdK*`;NTz&W{GtYkS2njD+;1-N>CizVQB^$r>UGoRH!>=oL`2p@|>g}QAzJwfi
zm-vjify=ov%uj!C9MI!3idX*V7{T<!1?KNp5eI&VvCB0u(iWKdb>O8_&7LU3>8M6c
z71#Lc%5jFd?K-k{y!ga*M-0<TN1MO@1U#p>+0Tx4_)+rN=8#_;EyL**C;#`<`Lz~?
zc)i@&=C)rPKj`VR&F6182FAOeaSJ}>d2DSZr*2&0s0AeEgU&0odY;+*S6=aX=HtIQ
z23GQut;=>Q9A$p`E9!&a4-3Bu`GSeXJ#V7oXy;wOIR^1xC9C<nBMUziw&-_9D?M$p
zxgV%L@)<M!53V-Zoc4#~C9v!L>1bvBb{R!8CqhP3Xg&=#ex&I+6g+ny{3#pd_dj7I
zVZPbrFGtIGnzGQe;$l&<j&%P6Dy~pP$4_Z)_{%Yn|CFZwH|*gj(&GPyJ#CITjpO+7
z$!)w$Utpg6+p#R3cC=Ajt&Ne<ZBczVRM~$Sbl!&xM-@L`;fQX8pM>Anuu7}I^RwMT
z0z<B0m3B3>ITF?82&xuq-rJ~3s~{SE7)Ke6s<aI(Tq?y1&0USEv<qDFoJ5xtUu#t5
zkj|Q#`<hm1$2ytUnpSxraO+OHFxK0=w;8B6%n8jv?N?nUw2`?FR6TpYdA(Vc`}OP5
VW_)v$C1=eHlq}TT)x1iR{{s$HK$HLg
delta 39803
zcmc(|30M`?8$Uj0?sd4J2`I`Yi@0kknyIN3xRP0cfVqon?t&m{q8Hq9Nkc~+D;2cN
zB{hu{wKB{~ODnQ6GZ(Zp{AOj^{;l|Z&hpN^Uajx%`9Hts`A>PzyPx;GXPGl+&KY)n
zTJO|n^$MHnVfB|ce=(r^SxsE!&n`cGuvYP6i4`wS5J2%_a~nGya^UeN+V1`9Wb->x
zMn1XiLTUC}*@u2CS)(UcXF=wbJ9<V`%#-2^rM5ObMAPbPL0b7DvscX_ZH_LG>^l@Y
zHB_?87D;-LqQ4ogX>nRjExy1!T{D+yW=deFo}Xoo3v8(u=bFz3_Rz0QH!lVbX>)Cy
z)C^`)wJR*vf@WcXq>oY-dCDyoXPHxia$D!;Nw!O|Pt=uaWwRuGj6*U$sxBSOFvEj;
zy35B(?tG=yQ)!heCu%ABPSAN~$;t`J<ttd@^s;Per>AoIGpD3yWJ$WA?DPrXUTIoT
zlAxRCgZoG9QxS$KaW#hwx4A3J>{IJ|y=<1*rFKg_InSI@yJbX~DkeF$J4_a{Jl8x>
zJD10>O-OEdey%XDMJbbCLW);fn0Y89R4<!io(}1uS1vN+>h#bv3e0hJhJ^b@VjxP)
zRr;A1>a>c;Q1KY943^iC+*1omn%8Z|SiCtbG*pKZPlUE?Q>c7wpoA+_-6yEvQ84yZ
zQ$fu#Plpaco47Y@s9uq0J{#7u*EQ9ysQow)aqi9*6)d7psx9dWi)D#XP;ph%P!*7E
z)(Y>T$IdYagtyc)^328IL&AfmO5w*<ZB<jXrI{=ulsV1th*o-juGuT1hn|^du83#}
z-ws7Q-KKns)N3U3&;~`w%#SLkDD$`B6mPN_<IFjcL-g2e^I~L=Hkqm|_$%$5NO_S8
zl{QuBQ|Mt_J3h-CP;VLUhV^<z_{MYzV#<}r%rEH`mDO5r-Q0ZRk(K&u^~}^q2k86e
zm@hy2By*YVlT8+A9!rj7!R9x~EA_;==B$*Cp=omyw4|le!z$PQR#la>bY8<<@a15z
z`C&>g{e>20M4#4r*SY4yeZr$2BFXQ5t*W|mXH`{I??o36hJyNL{v2~|pE&)_xh0$X
zMC$(99N0W^sIZ{qSXxb{H=SdCnSNf+m}|b?ces9Tjv4#di~84dN{&8uj_Lj8m|ObK
z)PI;`HX6`c?>4uj-+;*uefoU!y}{k;Jc_CdjLC;bTGEo&ijf@LJ>QHQ(iOEZdPuiW
z4e^2#NsG=`6~-Hy)_+LS634;n=Aj|ILaKvjGx`nSj1V(sXs@~<-i*Y7iGw7|J+xP7
z+)b?hi9)O8JhNhG?^-Rvler8{?8j9)BAaM)P)4(`j$qyDM$dWXx{U66&v_+ZXGAgo
zV}2;a5eG^ZFuVbeNvq+b^gi><)x$g1-#tG;J1|06xtz4<NYbKrD(}65>R9xWDD#+k
z=8wZ0>K*5qp(DEKS@X=nBbw>%c_j-+1hBC8(RGp*{ZP5_QdQOH!lb2n4P(Du`zESl
z+sOT_?loSbxkp@AV!y7cN=;hy-=rgc*XEfsM#nfFT=bLq+~`K`>w&>(4PAkYzt=VG
zdSK9?hOVR~<6VG8#9t2_5;ZCg_Lp=AD|lqXu%so03(-ZA3jWBzI0xFGhCzYzFQE4%
zEqP-LPv%zz&PxF>I`D_G*1uGK{}svsb+}kn_#0P+aY0a)+|pFw{P{jc?gI`lx-|?s
z2uQ2k<ZB7cuJNM-TUI+(AI)Sz5rmZIISqrZ2lmh5?jmQI5qt$%IUYFw2g=Q@eC8v-
z=)k7$!spJuIr;ibE}LXT*|-50l*<UFR(a#1PvOCw(9m@wF)+y45T4%^QD1Xwxp$xp
z_wq74Sz9>q!+U`(k4`Zl7F?rY(%iRW2M$ggbT}{uMOn2=WZ4TL>=TsK?ER!ee>mo!
z2wYbDH*z`xo~B({33zm1OgA2m+sV)Ug9R?z2O)Tq*U*hn1<!C@*&?;(aT8PsZjV>o
zi7O_!a6ne7B()bLDgRFy>4D313=|?Gkx$T->rgB7l6U)*WHI=*{JmFZ%8Zki%tVp>
zstPl4^M)7iAvgnSG`t?z@6h$Y0cAW?&W2YGLs*s^7I?$qwIIEnFYQQr>t7_X35f2>
zl=D2t;BF0FSM<NgK)UEGbFaG4VVXXfpuO8Oj7OzTP7Ev_=@|>d9FLj+9LkS1#m~_A
zH}@fAm6@n^@b2gbUhDS&RqG4rK=OhF?cnR^>Z*%3yTH46D|CmsZ@UH#KAf}&|Nl_g
zRglkDehaq>N4E@&ks(%{x2c&lzPpZ5@agfb^x5akvhmH0L$Ha#?b`SVzxO#9H_r^7
z(A@Zu!<v~mA;S2BgVoJp6Po&c!a<B_g%iU4&T<&T$}Yg)a2R7&`GmVWI^ZeBiZ>Hp
z$iV1vHLP+?MOD><OVV8rE}5zcYUKh@4md@4-#Hj5<d0NlRhaipZ0$aGh8OHe<!YpF
z24N0Uc?M-USS=W#lcU?~&`w%3KIq{kqqH}JP^FazAd$$Ur|lP`VJuAcUPPl|XK*3k
z;mv&1Xyr0s`$s_ZKr=`k9K**rY4-&#XMQ)aL4)2xLky+$d1gyz1tlHnu4SSs4`xB3
za>+UK?n&+S%q(-lqy};AVMr8AHiXJMppJ?|n~=qwb=KTHsa3ZLUvfW=R5s-45(mG(
zB*(rPFWLfcK$Qc}2?0KW$8!(7vvu=qpFL|foGh{(Jo$yk#^Y9+KQl}VoHwJ2-%rsl
z&sX-fi6aH^drnNWiH`~*B6Gg-PMg>f#H6L!8Jt@8SF>McY-1iq(Pf$<W(N%Tcrb}d
zSs8kk7ob{xTXTJ8<N6D~05<pSFcggsHHm{0ha?WY(gxJgg}oPjQ(0+V%50_YJ7w0*
zYN;<gW%kH=#C_^CbdXy<*Io1j8gJz#30A(LSff0w_Z4fSV)X)EIqGw+didL8{O5vd
z{`pFxsn8HTVxl6|0jV%BCPB-0hsA4C18!hUoQHB@%zNl(_-gl+!W)M@2v^vL4%2l?
zexSyy|JsXa)u_a^r_B>n$_6fm7ar*=uOTxJTFA$NpZR?gLqTN`cLbd^Y0-a@7M&Y7
z7*?)~|MVt#E{7X#5N;#OPnzFlx6}KdHXG&KsXzC*nVK_DZ+X((l@oqP8Rn0H^UuQd
z%7tH;Cv&3vZh;*bO<MGpY;jvop@Z@MaBvqYQqop}mbB=5iTuh1&~>zfukf1Y9pG@~
z6cA;*<#qG*Npojz=v_JP{FxoJoS49a*<JYv^Ci!Us{6`PVDfGOjf79lWmE6DuYapp
zw4LddNd-sv7_eko?A#AoWh#tcX^fUu+5Zz1nu0Hf%BACj#87eWQ(i!~EZ|?ffE~Ek
zOM>8G)rx%}9H@pLID2rfT`<uS=6@MD?~Gubuk3~lyra<#mE;?Zc+q@RG)ykwiAG%+
zj+&r)!jb<O!qFlso7=XF-w2pp{8J}T{Q73-w0iC<cik=oU!6n<B6tX{eC!nhH8}Ip
zC26TUsPa$Y{`rAdYK!1%zWg(iU-=Ws9fO=#t4mGt6H~R{7Mc5}b=U7dv10nYeog0K
z{P^%Avm(8J`>#&$_@K)9{KnZ3DdKw?17_unk9f4xwQqOeew`S&*V%BPF7S1gkCc~`
zPjAL_SpDyeCd~$X2+5^<QX~h@#g$w$5BSn<4P!Cfj!C-rNKUoN>*hmGG;8is((y{V
zgOZNt(g~2xl+wi?m8^ea3#+|w5w+Suh!n<}cbc00^Wq{u{s0(y-y)<}7$V90CvTM<
zJ}_73#dRKe9Km(t2V*F2!#Uq`@Q2FY$p67bZv*VT=u+hnZUw#~_cY{w@8>z6RQCG7
z{3)-Y-tYr6V%A_b#+*K@oxbdN$*x(U4A{}xaWaBY@0mZ&zOU9J@4*5}kAjMxGy5wu
z!5z-NDjDEzz*v<zJHMX3;z-Hb{Ko8o0g(074xeMkU)hhjLuKL-?%ZMSNi2^;JdeRU
zmA#HuRegI<czeF`zN03aGpqjI$XG2sY3W_(kfuJk=$8i<eb&3#B{5ta`pZ$-<eieD
zIU5*rn2GaZ@VGW?UW^t0l@VfAFsJ_%^a_4;12fdzH}6zOdQ!Z%??b6$Cw3c_GfYb4
zZ^{q=tt6bni0^%EnP<K~e}R5<p4n%?Id|)ptp>Ik9NRoL?crfLV>)*1K0UK()9H5&
zj33-gFj~aMx4g&0PM$n%)Rd+JTMy>UJ7e24Yv$o~pE7zxQ`pPSnv&DB4QF^PwH_SX
zB37w2YX)~($M=1x%dl+ou_uG8)3OoB=J+S0BI$9TkJ8>3N@9M@>X`eUjPJ#Vaen1L
zq*ArWQLNI;-jbfsQ}N{*i<YdOyUTh>Iv&T{EX=KyugDcjEhEc3urM0Wm|rYx)x4N=
zlr35@&>oaF3YG0L*jkjE#K%2;(E>q*Qro4>d~lp!dIFd6T*{VKH{rq)Av=3s`1=8~
za#6CGusEy_EuZ<;WI6g5zsgmagQiJ5SuMF_s5Dch+H}d!&oj3zj&{eZRT5vbw0c6T
z47AD=Um(?K<t9Vn70Ns<hp|XUr>8LCNQP^fuUA!_0-SLevm7T}MYlMGH0Zr5cU7Su
zgwU^It7}>z(xZ6RFGqR}>xV(_S5<9+V$iXwsxS=d6-XB%y@s?HsfOX9DGbGyS5;}i
z<B=vKEk?QzX(2|(ewgvZVtAd7v=C`A(p^Y3<eTuhyQ-=J0H17JgP}sCu4)iK8jG|H
zX$Pd`NRyFPARUAB8q$2E8s>j%kOm>$g)|oFQKa!mE0882{S9d{(ja7<fwU>oOr!}&
z^W6Z-02CrUg|rOm7C2PS=}503EywZ@pG3MagDykb0ckwuvKdH|k>(@K#KXcqq=lCu
zhe+;18jA(8GNc`lmLp9@T7h&7Qtcb)A$11<SOXvy!}uwr@kp;B%|IH4jPsGkBP~Ok
zj5PL2RaGX^%>N*iNGrZY2FOsmiclksM|uirCeq)K79tHpm1cYgeWb}4#up;p^&_gJ
zx;9)Z15gOTa->_3Rv;}ydJSn)w1L>`7+R4gARWW;8z?E#!e0<#7%E2k8#jnF2=p?f
z#YnFqb=5@1=BA}B^o%3siKQ)Bp~;pt2@k_fgSQ^K+5#Aaf7i^!W$oE1bH=jvt&@*d
zRiWkCR|uf5@h{1T-W#+(@vp+XuxxTbIqn1DXh+T(=76Uk2{?7Is%iwP-`UE1<>^P*
zHIo&FhA;fAstPpw;+A_Ot9-L{VJKT@4lAr(El=RXfKTUm)OO&R&_fHcFRu4C=wpP&
z89*5V`Uy}!0o6bulR4yW2q+2q-lgXHok2*rjX!CiceqV{6F`r@4SgQyO+iP)v9Ef7
z)^hpdRZXj&<mAMfL34F|eoZCS5TqP9TYS8#YBFHw+Gos{mPfHl^S$LU;rlQX@(vT1
zu7}X>GTDkIY>U}!#nWx3AF8Ul2Xw7l&7cV|%ZqdsGC|M{;&gR!SDY~B`qcc-is<V7
zFg5B3&z!~2n6=G&!WZIU4Sv`cH<JhPd~=u?8a^FdG<5so-0k3In#;}52Au0%E1khz
z2^TUs9n*H_m1oRj<~=OIWGh2i2Qz$S6Tpdplg$1rV=y6^u`)*Aa^8G)WlQ~Nh57!<
zmPQzct&R@!@0E3%9Ap7CV-stJA#`18y1JoXXxg3l=R}^d2>5&s|MwZQ)2c=R`M`&B
zyoNb$RkVKeGjp*_wykOxu<K-1RV2K22ATg^)q|n2tq$$q^x)0$<8jJE?G$5Fi2^KZ
z>@lh^O;oAkH?sLdZ9Wgp>)^+m8&-$1rsl`1YgfApJQh%Bco@bWAoj(r)k9>lAZ^94
zOtZfDO|;r$+B3syA5exdErnMTUSwG9Sef~w73$uC@dw4QFP>#*41C3)d)EV}r-8l(
zbl*%TfL?}m&<5O!Xq1EPB5<dk0>2!5vHULhy@_uYt*Yj~UGYzuv#d~a`<i+U!jMZG
zx2CSPNMA+KGO!rSUo$yhde;1Q&0P_vUb{Kq-5}SS4|<hmn`e43SILZL&O6ME=km>d
z>waU-l9cr`{o?i@JAToug_^48O_F);`MUo0%LH@c3!B|Em%&>ZanbaJVnp*0-bVRF
z^A(CM&09n@a{=`_3-5r$6|3;O70#D|`87oRYUbUW>j#%98odGwGDmHWspb173h(+>
zZ|+v(Xn$Rk6@_<3mp315LlYSCh$b@B7imI6I3_eFx4aPYYKTMGp{dInr#b0`%W*Vw
zA$v4~p@<|Gm_0V(g-}7syW4s!!ODDT&34Qcj5ideVLe$LhL;jnt0dMEZ!SuQCb3C=
zcxO@i(WC4YhPM`_b5q#ee#})`+K=6@2j<U75OFAEEY05$TIw9g;u+qHSnUR}26*LC
zI&ctkGZtulIhc)L!PWyq*hoEXuC-|hO9T4H5SEVDEmq1<mZm>6w{-JRHc-baoKoj-
zcIXwX`js9D(f8`QHqY7+uGh6zh3f0|>I!)%3`mGVSU8Y4fxN1>RK$r9Aa+zpgGeAf
z=UIm%^-fH)sv`9Uth&{p9+D8Nr%c8p;o61u^td?wmP1_Y*p>_lu4nK!TSJ1<_MBJx
zX+8aDHAK)_AE)o<#qQrgZ-|#1rS1m$5<Q5&;oukK`}jL{&$l8P>v!s7=2=ON^?SHu
z?#6l~-fviIdFq~5y1%i$%*mFQCbiJ}uV4eLjc@3CSZry(gL=;pHmkJwgnmSic&g_#
z485#N2TiLXVC^~9W*P4`R>5cb-Filrwf8ebEi=oSbW(4r&&VoWe^N(W^Ruk~p3<B1
zw8<BGEKeW&La(pS%d%R1u1C0@W8*9gFj>~JFZ34rt}N@dFZBLAt$SLJ*H30yNvHKL
z`WIQ&Q>USgbmwV3EJ9Zg02^g&P~5Go_w$<8Xo^*FT5r<qC-oq3?-!Z|nG}BqJ<R~?
z!{Rk<C%d#!)BLRXFZKH2LB8fBcheMW+?RT6y`8@79c-^6b+umn5`LCVu}*!dH>~rH
zBDY-wm!4xqNXn;J!DsXiK^MGOudtt|SpCoFP0%I^&gjk1R$n@URuG$Ql|weXgX-HS
zcMJE^9XD8-_6qBnZPh=k57jfWtvP4)&icY^tL!W)73r^M(FUewTT$oq2I2FSdd(7B
zJzd+t4rZ73Kc`1By)4_xJ+HUbZ)98B&+D-vzbSd)^;YKh0>)IU;ym07vRoB<TzFVc
z0?MOlDKebAIbsQ!RpeM*EA*!N^&D$rh2EsrUrMZ&7^ilyra9L6uk~;p>7fd}ar4fK
ze`Kie3(c9Yy77zG5o^tBk!Sv3jum`CZ`^vw&GI~32ttjo#o1~YBzLkWp*AGP8gl_n
zZ9<Oq!6ig<a*nkhW_Wtxf?iLbl4Dg}KqH%-W8HaCkE-LRUQ3^ZJK6zO0uys{tRWYn
zx+TY&dr@y#Z<q2Qr;DadFO|{ch2O^B%d!5vtVhK^qa+GzYnuO~`a2LRW*xSovy9fX
zm)Y}*U6F`9cCb}BR?t^^-OeXe#P_fVmB3JhmSN$J7f{dTIeex#pvZuive(d6f~?_R
zq0YNbwVwG(kFArW3gK0EX;ZBaztZEHk5?>DxcOQjzeKn(x0B%>Kh=u(T5sGrN6FJ2
z-M@|AdOt7D@90@yKF<BmQzAYMAvxE|`x*^lcrI@Uq1GE;>p`_9=E}^dH9R#{wg#ks
zAiL(px43uv*05u$)$0;k!&^7Y|F72Y)>Lc5CA5a}sj@YEGS#{cGdj|U%V-TJr&=8@
z>n++|RAwtMsj2=JOL~o$xF#wvdz_~2X46zDyc)^3Q!VQ<M(`l(waa=V7G_<zj7Hu%
zS2PI5tcY*)=%DIKPgVS1R^V59XrDNh*&en-LA>5OzV2P`-Ib<Y?{BcWoAF3CT85L)
zt@*{)PK4Tv*6wfgChkQOutLpJAl7#xuhw67V{CgoK?r7ImQo}Ci8pw;>b;_A!zXbl
z5YRRORReSp&^~tW7L<C7)S{`0YZtp;s15-$3{MiptPh}W0!jij29WD{79e!3p_O{w
z+Qoty4$9krG;I%yvldqBP2H`OOQL9K<jT{B6uFptTJY#W^7Kz(hkM#Klhe4T697q1
zpI4sV#k0T#Ps4>O_jIFB;GVuKAnxgZ+yJGgq>DNimuqw}YXFhBES?Ydyayoe`FQ2q
z3KLUdn1}0WtTldzRa>6?i>nzr>h#-}YQG}$y*oALHhdEkR#^eFMJBvh7X`$N6)}a1
z)oE;kaH)~_6~+4AE><$cWCbi15HHqFK(bVwtRa<psQa>D@j?V(sN#kA)T<D|a6}X$
z{-#1)xP2l1Qe>?z3Ni0@=*lv5fo0i<UK9{7!#O~x8aFboerq8P3te6}AvsitM+C$R
zF$It;1gZK~%{C#+OZ6cjUaBb70C?40rc$-ph^e2aZkj6@lr^0Xt6B-GP^6eL5eQBe
z@$jzH5_UsLJ>(-LqE`K<W-i#0R|zXtrOhhi79SHi36pKSLR`Z3C^1Wjd8E86VmEtJ
zNlnJUpcS&>XOZlI#+R1!s@}r=BEUN$nC82GjIFH#=f^`lv|QYns%bB=cvRbI(8|~$
z(4^0qo6E}B2R5qpPuR+iya<hP2^zmLcvZ)Y7nV+dJ`<L9upr6psC{)xhRqScUnG~X
z_DbzA6jYRX+5H{-B|N`?Vg8l8z^c5ex9ae)5_e<rsKmMXbVcrpQVTNL5=wu7mtK2I
z)9P4{ey2ycdHJkJhMissoPiZCd9k(kLTE7l6*Dfr0X}bI2b7fbhrg<o*Lue7B*KpG
zzC*UQ0NVR)J9W`yO4t|5&IUIa<lQUmq9V^!K4X@a&S!mpKZn&W#iL!7Z0#4&J@Yn3
zY<6`3<|-Jfr^5JVmZQ}e9gP-TcN)*C-gxO``@?ol(tEz43fou&HN}xz&i%mrp6^QF
z7qEVg0BsgbbXF#2_?qAwuaqn-e5x$)%?yf47J4WPPm=`^g>mB(;5e?2iA>WmzC(;(
zZ~}&pD#LtdgrZ8rdzq;W*F_v4EatUH(O9Zm2|r-y;#(+g(kWs4lukdTV^0G4UNI$$
zTu-T7`OXVPrQ)1Y7M}IBKo48e?=#xQa1m3!VdEwfcw=P!@Pl3_qNi7u)ll0FGN6iD
z_-4KpeN7Lo_APR#X?47&*LT-^2PT4CPitB!JC9ZzG(j5W%S8NQYww0c*H63~#BMZ!
zU2OJG=zpz-g=^Zn8EE^gHQG>(&MTk6!`2fV|E3=BM2;^2E~juI%6cL1yRF}NM{A6Z
z(DVf?7(rwciv+ZTbx=Z{`JiL<k9r+9(JO@VHMq5jy|sn^?!oVF*2_yz<wvBTP@Z4U
zcX_|#*;OLDVv$|QOrGUXFI8UQYOM8l@Xajz0@0nq+VXM?f<DP72)hT5a3`jMUd+bh
zw-;<~XXPa~%lT&&v+qFT+*g_Qs-}%lx#q(P=H!_$gPS@S`4GWf3(!6x{0GW6*{!6@
zAdLy&RZ$o|0sIp(yV-v;rKesK!0Afkn$##}$KhN|9=H7#pvM5c#x$iDFUrYzTxzK;
zwHweL=E_V!MdRX680m0SL<y{ItZ}CG(sezm?FyBfOBq6mUKAzbPZ^&=bT@leiKxlu
ziV6HFBg6{2p&M@E#fG7(+Wsaq`LO;gG<UEUl%{7?=Ua*74dHdgOXiyG+Z^R(;F?jh
zxaPaJYC1@>Qt|S`u&P<3!vtHkt5Cg8s@CosI@Zr8CWxZ#;`wymWDAcK!fh3oJk(u?
z=vKyuQes^biYf2$*22(yi`mcME`xe8<6ZJPpslRw^N6rA#)EG3bV8$;rL+g*k=Yzf
z1Y`liDEfesjDU=4xuBB%{WdiSOpITAx0ATGvCEm(u3youRw@g4eIl+J_pKE9`M;!T
zZ;0%wwqbD+&^A`AWVtahO_AMtt0L#a-5U2-2zB1g9JmosuRg3~JiVH4vl90TSW97m
z<>N64))wAygH7J&g)WA{R>I(FKs(qc%HYR|qU6Jj<@?-xqN}>-NyXkJd|G(y<}R8B
zL$Zsyy#=X{a<=TFUnn!zU<NI^BkagNx(*_;kA4PdALD)Wk`ne;uA`g2C4A@Iv{yc;
z^6~w9B^Muo@<T}Wxrq>vJ#!PF!z>eA$@$XeQ~}z<#;88aJ7WwyqKe<f7ATIQ@Ook$
zHC+2xBz$mAO5li^VTlsCGewbjCyguMowSCMQJs`jevn~K*U=xT|I|^U#qhL7;vr#(
z_q^^fw1YKL+}rdBVpT+7u23U!;_;ii#5m~7F7Y8C*(G`?Yw?l1OT<E>nDHLbd=7_r
zkH`RowjZwq+d+f^Q^H7nfK63B;tswOPh8vBG1VpJD-$&&L#xqc7u1Rw?-C0{e!NR;
z1+<MVSF$|RqD$;&rXqc()H{BpX<OJ>twwkA3B0a`+Dj}G#b_evqjI$xN&E)ed@e65
zo}pW4UnvW_AXTGJ@&`QgR+}~LVTcs7Mu4&b9cH^Ct@6M15GR5?*{b**tHA4i)2p{^
zG&w<JU5ETnJ#3p;ZJXE+=rykR`tN$2-ccC{`W*&3`Wi@U_aXUk#x|hO;|8E83_L0g
z+>1escV)g7@Tjku^sTm;Mo^YTNfKtDF3gM+X6iQ|sdOuN2>RZEP{|OCxAhJHI?O8U
z5Ii9?>$n}+72Jhw>~l}m&?=0mp>6D~n>pNHzM!{B<r%{Dan%^USmevcva`sI#<KlN
z&NG(r&6NB?D~L#8?RZ7<c>)Or)#Jv_TR;AR0oR-8pS4L$3glo@rtD@w9D~iea$aaY
zA@HBiGZzyAG0+@Sg8m{?F&_vQ!hE1fgs>WDK2>ssFi7J}sx}e=a-dlY=rtC9+bO^a
zTk01;dsqjrL1zv!6MZY2$2lcaEGfCKs@%$?`yUC9_`njrfUhzBt=M8<+07~y$_EzS
zZ(mh++rP-ReBN~gV5@Km`&w~3e{A<^(c|*CA9!xV7GVs-xPXfH2?>Pt!-8!It1(d$
zrrxaY@8BDT`6bOX^VDM>-8BNmj#P^LVQV!h%O?unm`WH;3%0vGT#)IW_|{ES8Zl4Z
zL^V|?ewWlV#vip)#4mpdP{OET);%E{#f{v4N-SY(RK>XOR9xxT7DevjMcu>pk?=BX
zcaR*h;IN04DGq<M%_R=%cMto-OQ(=HDEl7vrI*e&H?f4$_g+eG6Gtfh;-z$sI6@~V
zOV$N{=>Cm3LMKddxO25WQB6STG_ZAWhrW|oP};#-C{|Z%(cgOHo$(5#QIV~k0Pzj@
zC&!`AXr*$DnYw@vt9AA6SytuWdV@ASKIOB8ni%tcN|hd7;~n~`cUd@S><DXwq^hcf
zInpHCu<onU8$~yO9%?e(-%)LIh9<68SR*B~*vhZc<L>u;6Ux?(gVuWdD`w^BJ3HZ3
z30s*f4KG#}_3R+LDOs0mRa9ZQfwsCel6l~v?JdaN{TbP~<4Hck9fGwNlpV}O#bxx0
zw529M3f)jixIpK=^PIM^*A-P#p!y1XN0HpiR?ToU4SEdn(}b;`kk7^^SM7GRZTB1q
zzu3fU!i`H{=ZsDe*$0VV%uAY_v!xzgh?NsZ44Y|nVvf7qktY%1pr3Kme4j0ufp<HD
zI7K_g-oS(_Wk@HEu%HsbUx-+x4CSQ!<1FiKAQ>E~nthCQ3kEivV}Dm2V@;rxGD7g&
z$Jshe<x@rqwfxz~*}ghJN1G6J2OaDO!ymA9{Xu$M$m5X6aaP(4_!y~?pFjOL`?fj=
z<K7iIpF-tB#y4lAj6ct_mzO)}Jur^5B#fph6V6)?=#Dyf`hCuGxf8e4582UmU`&n$
z#nz|~#&Py5FKK41l@E>DN#MJJ#6*G_kARUqz}lian#e%C&#Z%3Qol12`66IXAO?=6
z9AWRXe(wNjCOF4fVoNOUwB;qxoX6PZL%{A4nClpuG8kApfdzfQo{Iz4USOJzSbW_c
zcn8VVwfERVGx3}oKgD{);dq_hu&O&9^_gki<#e29Z<O9q&C&V=ef~l#VXLD)($W!I
z9iK9sTD9(d$q|o(m)7`~94&C1(ptmQ{Jhe4UUJ;Wn$SkZHt>&MZ0$V=LMr~PX1{LM
zv^N*>6)fFqw!;yL7aZMpIAWrs7hx%ag-=Tm8Wv>NiuTLjF~?1_vR-$DSgUt9n(MO`
zS#RxdbgH!%3g|@NX`0@)A6jLd#a7@>M<+dgn$>%!Bify_m}fU`1G@SW4&?%RAJ7Zz
zNkIDY<V~y&fIO9Bd<-IgX$6LY)vP1v=}&R!E<no#6btAypcmPu=Wf<=rL)!SeGuC%
z<tlFgf=Um7y8v>X&1?#w`6{Oj<m76>hKsB$FeV7S#eg;nXcnMTfVQyDr&+CDcEswJ
zrdcU3J3_;2VrX5>9tFQOAY5D755;&0on{rh?1-rQxM2Sc%3?rU+3(K^-fr*`cL-iB
zsQ%7XSFvsYS16a9pzzDJX$>0)y3C{A)10=I?SPrGX;zJ0j)>6ipsi-_f-;74x3W(F
zzCF#l56rep1oJ8=J2`V38?Nl18iv6P({Em-{*I<|62Qf^jlD6=TC&U0DCD{lJa7Qy
zmsredB-Y!z91Z%PfTith>`2XMx+)LlF68%yv8^?uIY*aEw2fx8;OLJ69hD3AolCZ3
zj=r2rkIDuAu0r60fp29s(FEhBTM4@z5&rq0V07&<-5RmmaRo=Mtw&yQM8`f4@@iIX
z2O7!>4g~|6252kmx)En|XIM|Y;)rk;nw-}MloCLjS&Q|YmjH`cR_3+*G^@6gizI`g
zujEiSK#{9B)Dch`Ahf7Y(Y6{Zy)v#>h!&;8>{2253)t@o=qjMAfOa!3*7+5-(pU53
zYW5arQgN@Xxcw%@Ev)klYgmaRS|2#WDkyQ(eXJ1z{vvBvtZ4(4AJ@1aG5TuO9YoGv
z#lD8eV-?#Uj`Pd)>T32E=np|+HM;?5vXEUNWUcREd>p5&W~HF463{k4Zve8`<3cJy
zSR#wj)sSf}7*g#iTP<IyZDf;_6B~GI;g_tIDKo4|uc8}Nh+KHJI4r6b@YcWT2zB=X
zMb<+28mfg9K&lp|Da*;RUbV1TNb*{ES3tZLt^nH2xR|Vjn-o<gNIhN&QV*5j_EOx!
zmZ{tZWBB4-!?Apq8?NBm#3sZ29g3eHi>K!5!*+9rXM-O14Aqvmgk4*lEo<vFWuTfc
ziXug+>Pmu$tgZ<{m{-?pGpzDGj=HIwBCBhKfOvHs03@r+*Q8xtQjJ%aRO8jf)i$!S
z8P>gf9Wias0(+G;M^~ziA>vhb51>H;N(A(kfF1@E^&Iumi|A;3F|ee9qB`zhsz@|s
z8G+x5xJ02-=Sr~+02@DHordNe&q&e$kcK^Bu~J9<W=}nlAi_aGF~SDD!b?0BGOtVV
z9&526G1Gde)DhACJV<NU08qBDoS9Pfn_;5DFq7wJWB7&YOVM5T@aAyu^Oz#dwAPlQ
ze<6Ln)DhNVrc#Wrk4L#j^><5nhUYLMRfpX*>{CFS*%C}47R<D&N*$q*kAc3CT~KPj
zAV0Lwd=Qt+v|8_T)b6kp@M`une&4~rJ#3c}Z6e&^yk|iFuPqe=<J;K2nO5FDN7&;h
zl(e|V=5p`VO4DAp9fC!Q9juJO<v)eX?}6C`|2DGCN@PbD$ZTY!@j6F&UGetoXZ?4d
zBeMD^C|F`}@w4jhcZ9X)1ld4?^Ad!7i1_skCp;9geZDMYdHW2uGWR<|$8myeo}USb
zH_yKS$>!;6&~BbmjW<uJ#+xTs+sNwXNgupM6YhMYJnP5(j)<<Y5R~3GQi6fxEDeub
zgwI~;&XATy?2-LhZsO@N!+PX^BTVl*!y0qIQLo+_CG|E&uRSE`Du@1-8CEeQ9yz08
z%kPHk*j0u7O0wC!CVx`YPaqS#SLK(>%cFah2P1u+74VwlqTXY@_5Ev(F8a##R-4xy
z4cU#-$6j}AaInzQ-`{i$F4CJ7S^b7O>m%JfbEwmy#}`<$hCAc+gaYg3;m#I%a)EVb
zxU-*Lm}j*c;f!{TDB@u$$+HHHa7O55dDe7Lqu&Dc8TOv?;4FH=8a0c1mYvA6ULE1Q
zU$0+m{WQYasn&pEUOGO3V{O;JfaR8AD{-W=6FwGV%^vBDbBAr<Pe8cySBrW03EX~x
zZs2$RYV|q#xkTG&)#5nX9yPWBcY3fJ^xLFq6BQ1A!@<Cx2fmRNKuMn^D~Se>YH{7c
zFJy(_y$dPq5QyO-R?H}8J?tg5?i=NdYj_2mr&$PovkRm?4Ij?nLe>EExaY08qn!Qm
zX&39<C}(t+?kIdA`vA6I0JNE{MI-rM<z3i7)2g|8Vc0KZdqBS?WM2Zr-#~3)^Pc6h
z*8QWMk@f!srI4)y<sRXu31}i9^s&33@aHV+iP6rOS|&)_*fkh_iG7E54f;!?ozXYZ
zH%R)Ipx5)rZ)X*tJC<*8Z)TGKHJNSIdE6N{uHlP3k7mzsiY#(+L)?ebSs{A_#06ll
zV2OYVSzG*a@g3|J+!tFb@nSCCf>mN%Eo6}pxhPa>3KfRm7U1k;2i?|&$DOe)-^rJ`
z(4b>9V=LHs*ehh8;CDKNacyQv@bQ7!)(?+6>(*T-Y(ECdTYxsQjSvJca*Q(~ROb5{
zC?oJRf!aF?@WI(u?=jAL?uDELg}$J@;7;I(R!#h3*gAh4+(I@Ayj=pC4Co}+*vRq#
zZdD2aq)_pf5L*J??*du{C<fJwg|pf!8heInS~N(!L-^P9OoDf^4Q}h_G0v9m$CYGI
zBi=gN<N2VFeU2*b2Cc1PsP~pIV)R0`eywCb&F;f(YPvGuulzDzEo6N`9|r@4tT&)#
zfSzXGq6yt1yT5-S`vnXxQpkcJ{HlN)fX)C~&dP417Fi8md}OO#0wmSG0L0bSvstsP
z|BiLW#8v^@%`TyxHbK#LvmXHs1th!W%Gp-Oan6W4UIqoEf#^{$&$fnvboZx{^Z`2P
zaaF<qRl;huIkS+R0q=Jy#{0O#+Bps*NDoy{PLi%8yLmJoMR+Dh5nR851i_^{P&HQ&
zZpJe3T_Lwz1m0K54n7x=!4Fo3$Tj0t>(TK}<Ds#TkilOnWO?v!k=-w2&n0B+O9aGY
z|0N(9`(Uehyfdu36ylMVLOjx3Xg#Z|T!_YuNL<gbuhfXs7=;o;F>Vnrvp%z}zhG*}
za23uW#LebHX<{Dw3_GVfWDaqK?qx;oflQVW2e+j4tV(gxz}d+f6j<34oEO>OR^N%v
zu6pn7)^ih`4fU6{mmZqv?9$&gGzep{Hqy`9z0(;|8gs;%?pHTCD?#|}IHfV{jb}{r
z*E6!L#UDEF<eR?tedzp^@9p0Gkuxrwjzc&uLrypFgGC#EzN!1_N6sgE(rF0CUz}Zy
z*mNj@lMm?OMAouZ#+NkKORq&_Y?VGx?u_Jnp+7HoX7N4Ey*_qs#hS%WA3Mw3&1;YA
z`PhP<V^26;j#u>quFst{9ObhLfxWD|w(AXR1OIqdp{teSiC$%nC$@u8qPtoOxle@L
zc3bY`tU|p(-AfR8K|kO)F>AYHl&huVJ>3<nkj^j`;hN*hhR816mE$<^m=NmBg=&!o
zgzCFu_4t~uY{zuyE_1v3Eq4TLFPK-*-^m<5Y5gFv+Pxm|e>B%gfppe%r*8BUzrEMG
zo-egNaSnGhS%QexglsQoprh&6kbFyX%>jJS-OCft{9G&cb7viV?ya=_=gtCGSaP0<
zWFUpz$ty9>DnIWG3CkpIFu49a%$%E_XZ?BJ*%zN0D^06#&JJiv2R$4$y^YXu4+l3b
zDx51P7p(r)8G?C9Y29y~xeRygCzoAymM*(mt!_u0{IGU+Hn3N?vGJ>p(;plu{h*04
z$iUGL$vS(`aFzCIX{7K&AX4DfRz{`~<(ex?Qq;DO1fEi9faY*niS3L?-K8u}ZD&01
zWT94gClEs=(OTERa9MYB!q4DN#y*xPRT2d^#_H4ANMv_P+MRPt;aFW=4^*lZ_ZuDb
zU=I>x_3UTVw(=h^f~p~R95WH60)DS8_9B?p&US`jbxZ*C+XF^zegH)X0tvF(CK%0d
z6h%^cL8f#;f^mr-N|B@*-Hhncja`k-{x~pjzPEv$Yu55dj4*3(Z^$-y7!iu%pW@<I
zQRiV27^G5oRlO9_GZ@aCR8NL-h)^zlAlVqkqG(rGQJ#`t>(Psoxx}O(mrE%;7us}H
zL{hJNq+3d}<5OXlA2j(e&4^$frHO~qU?Mo3o1jfpMSo;SN5efPdP?a>`$BqiU!#r#
z=US|T{f!8oo<<tghl;s=q}1N6QWUd~(jJ6n(|Uk0$;Guc42Ktr^YDA9@-xGYOnz*J
zGk7Yh0i2-H<s*y+dQ_$w>87icY3T&c$g(~jWz@TkNSc>O>-k1KYeJIIF^0~GEvCU!
zn;KN~ltNZ!S?<wB<b+)c|Im{kfpw?IGBRxtfuV}z)Uf}oh_GCbAc1b1luqE@DOPcs
z5#FD2JZ9&ZM&M^QDVabr$a739B8+S@Esa2!=SA=Eg|bmwc~DI=U8$za*H=aM(mwvU
z5%m~phuHr2C6FvIEtx=ahiQGMD+MyBYw<+{<0KD@Gs2^(<dli7we~8JH7v%6=s_|R
z5|ofY3JDw_kSzH|$GXtQ2(M4JNXpQzpJ}OMIfQ(*){iq{qAFCxZg`?kC!GAISh3MY
zcodyz`rRW!F=#GDQYaI$MwyWcnIyVzz=*z4Ik#6hSD(z1&m>0bQY_$JrPWn7GH?98
zr;K|3eu^I?tom<K=uUxJ@na2F$lHpB9!KCiR{1C+ste5xDm)SdUbRW71pZ``(g^%n
z*y(z!oiwuZovqZDKvHH}3V}Gc#oKsmnw7ZJ>;x_oSd%6!+SJS<Lh4#<yb(1iR`Kri
zs1ex4CM6Tt#U`Z^*h7&V+&CdU6;X}}gr#O%EZ7LYoAPXLTT3I5L}lY842rz*x-*o*
zcv~@@z&x9jO5i-3lujU>Lz8|Hw#ODsCGZWKluqEmY)c<+)ayY`<k%u<1d>7NJz?8D
zf&{`Aue79FZ4L9aWoS1XKL<A}B5ULh(}ocUYqKU9(bgN|3|Ac@{%?L7i7Ib|zW!>q
zb#9UoDUlat;B7>3iYFVjEZ0tByHz#WxYu2#iubWy%ya?+a_DZ6O5kKga%j1kMTC5)
zh_YG;JMGaXP@gId3T=ttdChOCxCsW_P#*l2Z6)?I!n;zQp|+u90-M^TRI>ezO-d$E
zpDGMR1t}J3^Wk2R!g7R209k>`@>nAx;JV@i=VQ2y|CuBDKnOWNClE1N*kuI|Fv3Fm
zE6pJ)1?LsPi6YKm#eo|H;-n!jasVmegyF}j$kq!NA)$y*@t7cRy-n)-yaK5va4)0C
zBi8)B)rJf$Rie7qvB+1A)4cQwmKhOl(xY4hv`qYGB-JzGiino0Z2D=rTILuqh(MVc
zoT0_p8NRF(h)$Urs<{vy3w3WyH)=-BR`PRI%CuBcLHT%~qKQRW-dsM5mnF(q24*}V
zkoDn(bCm!uT-4Vp<-x`raZ`2EXNKx4&oPP<ILIcY68MRCcOayl+N~iw6Jfm+6+RL;
zMUg~jk_bA}v|QPlx+>dAG9V$8C!OI$v)n72B~hO8IOxejhJg-p=1F4#UVPD8hAUU5
z<jUSc_=H^9TgWt>(Dd|{H&mu?snpY3h;ztxE15u?#1zpdC7i$%y~PU`y@f)0(qn=^
z!;WWPI#@}eMsK0eBABALkfEQIA>NvN)i}*dPxKbjqg=gv3(-z1n|5y@0%b<6pE3+p
zt`nUy&Gn&EZ@H*4_|}d>Dyg7+qPGx>vb?!`vbXrkpe2$=?WCCiFI@B%FPvxN|CjP0
zz<P6r5qJOb(R6!BCveypBBc^IN|9uLAng6I!t)V@$_5$0+?v3Vwyi+~lD53T6Lw%M
zWjct!gElGcptWJP5f?yK4~(_S3osC}$E8wzl`?H4f!h^H*On&}OnKGPmM3BFJC6I&
zC&eO*h_}LacqD<B6p5Q6_&TS#HM~mEMhY0<Z(ec)Q(118B~)&@d`sifehz8&bDd_<
zMZSfq0Z&i>_ws)DjV|&nU)ZhXQ<A==^km^%Se1p+yQ$c4TI#vEZ@#6!S<yER|LodL
zp#VY$o#}LOfVSR6GO8F3ZN(Vtr92}p+E<uAC3?y#W)o4qG8kgKMEDRBFj>{0GV6`=
zA+`(grztr+LJ&yt!;N5<gZ!kF^cIW}uDA&&WiLAml<C{rc}O*CI-AT}BxQWWtMyU7
zl=|oeVjZ}3kw;P7FQR;e(6ak_iSYQeEK+&=TV5U+VYUn5n>{fjkj&h+tfW^uM%Zx^
zPRfWGfu;WzGdi$+LB;N$TVN|S-V9dhO_CkOrw1cdf1top>YE8~f3`U*E8mEW@)bdN
zZHcks4%fB)zd7))`P-TzuS7y0ut(4o0;zDfjXCK(BFwl6C0&G^z(<tj+nS-{<#bH@
zra<0Q4jLXP5;sY9x=EBmO{s4ty4I%1W6nyRgUO|@3}Wsj!kctAKUI+1dho9fcr-<3
zT?h~K1f4)Ke%ok~Ud&iw$4xjXBkBak{9DxN?0EHY2_i<A_7a$2lTr!nqDb<NPFVMG
z3AL1RDqsPokCT&GI!WGUoK=6WQLjIp)}Cw|PbYAlO-d#3`J3`1?2t#0z_(;}UFi__
zDO)6+z_T_fmB0%&DV@NJ(vX`DZPy(y%Y(QOxZEc7C2)lzIkl`;-9?1a5oZ}K!m2##
z1ZopxCQvVS7jewsL1Ys|3~gU(+nQ+v7TTm#0x5TROdw?~Cr^a^?a?KW9D%M4P7n(w
zZaTEwO8E*A0_mg-=9L8Am`o+gd%;~q*v~d~))sdWVZYkc*$5tCzuQ#zK6laYZiUhT
z7a3E+T4u`Z-~oYs70Ib_4+z`k5hRe#WJrsIz3UMaAfL6-Nen4T*aeRyf&aDh<_q(L
z{p}GXP*W?o(js9FdtEn`z-p;CNmf&s%Ogo3o#K!t35)ay64*k%z_2T}iK<v#%h3of
znk3_@<7)y0Q_LOw9RMK|buLIS#asvy@{0-r7Zh*`0!tARl7zHSE^|rWD`%qgiw=Wg
zh+d}*@5mIl#i&urbvh^rk!4C`ofHXau1ZXY-th$Mu-G)L1ocrII)e^Qq~JsmXNcVm
zk_j9-LFP<KblMs@d%`3{qDf@29rKPwUNZa<tg)|50G3{I)p`xd`ej->rW?@#la-z6
z6U1AcXgZbtFZQE+W$<E&%(l!FW&*r$@l@!A^Or{fbR7Me31ST;fJ*%G1Z&-5qh0`=
z8b3e5sv3cHe>w~9XH_mXq5|ktSI3D~;1VOUcSlvLRJ({eouw}NQq>)`A8#J@8|oFD
zCtm;IjpOcxMnqFGh??-V1CmEwR8IVaTN~Hle{M{wy#<g~&P0p9;|QQ5-^CM!RU0pw
zC}I_IP~}2LvEg)k(Zte;LWiau8ctpN_5fni*=eju5=e)OQK8i9mdOf@CaphJ=6tol
zONLkSV<bc7yjMPmc2?>8X?b#;kO*Yh;0&x2lJUngt)olPP^rUyFQXDor?g4$HV#Gk
z%Am@~EXpMw!Mt$MP`z-@x9<AHaOu@2iOBoY`M?%d%_T;?-Yrz!(tv@ApySF#bT^<@
z<ISUSikc7SiB{K@EVfdXQ8xrqCMbiiZIC=BMX{p(rP--5%%4O?S)HFkm!ebE%O}ae
z5MC(D-hGpDX_xI<Ujp|l5~iE(B0_#tw_P~r8P=s~I{DjhvIuK@m=dLf(M~P<s8;0L
zZaq?RA&{hT!Z|+yLy8y97oOtjlr$Z4hSgKITJ=att0aws&SIeQ!nxJE>EQDedy<h|
zp%oD}!!C3hf$qsx#}(++qzIfF3fZJY_p20*AW8J%Xri1`9Nx0MLMWOk&6Ylu)}rYs
z@xSOq`O4fh<&mZL&BhDoP6yDb;J;+eMa8Q-Z8v+8kWOHCwaejFku)z6{zk1U9kP7d
zUY|`TaE(n$C2*}xN+%E}E%}a$u5{A!WseAfWi~05z_)EuI)U#<Ls9=?=ykd7S$h8}
zBZ_sQlW^xfwg~*oCZ*coG&%Xe445fUx0X-Et7-a^0AU^p0;6nFDuMNFQaXXyWCJ_m
z?Ve@7>~qH}!NDE{0_o@$`X_;NJyHa2w@G~od|8p?G?B0pTQHr#>(fMcx1AVfi=+`q
zIpc;y;9`#yfje!|AOa~{?nGJ<Vf$>sGy=~~YfLd8bXKvy^N0}$UwL<*OaUVY@i0sv
zImok2DI$!H7ebIg@<s|0w#Uvig+OE){H#$gmV!;?m5Rr0Z-rF6uM}VPRw7g|gpt7S
zris_@U2oM&A+3|PN9m;XtxZZH5Qg{$fT&vyrII0t@X;ZP4nR^lFu;(=*?$!ADuk$l
z+gS1PSojif$rBpB#~_Mu@)?yupmIN8nNnFjO-?wyaDJb^oAB*+1pcAHd$kxbfgr^d
z$}k<dL`DI%jhKKkW$?J{4K|K5m3&K1z}r52FZ@C&HWJYv0|}f??Kqvl!!{|Ez$1#}
z#9w!qT}0T0>7s&cTa9hSGy+Lfb|}J5dISlCE&kr^|D~;{c$LlFc6P}GzB^sEiu)9+
zt4D;ux2MZ#MGwW=J6%RLPT`HgM+@$8P2qKB*pVdg-6x2YO5pp7!~>FEL|9^;h-Cj;
zZKaW|6x)1X0;ANfSQ%WxVigr5Fo7_SJ7Sdk7Ee>iI<q|;L}2tx>AV~^2y3LM7!wF=
zev=?!9X)~sy6NzdG)dU69!Ubjyh@f*L>LtZc@x<BCPBj5+94Q3U?NCX?Pu}cJxT=<
zXIVuTCA@{AJ8?SLT|`(XJM8HMLh;!1Mm;wvQgtBU1ky=USu=!DWZ^%76xRO_K^9M7
z6k3=hkPiR-&mhy8Hz`O#re*|-1X8eXYEcPRoB2j$ECqqiAtQQNF(fR0hUhc46h-ys
z7(<Hpqz3US*6mg5r{!q5E^PWBLZ=zlg%}KOB*<|-P<wFdZ@!fkFjCo^Br7j~4kS;U
zVg2|#S_R>^;dgog&xyo`Johz#>U0^--r){^H{p3IACz#HyNIx*ii#N;flw^ki&jL6
z)FS@hG!6|HHHd$08?Cq`bvsBB_)m>Pls1N(rsK<45lAQ^q(D_V2~iMk4hf`0r5C8s
zEK(`2^gLeaL|CdSodh|~E8T-9@G{^<pKojf_9~k=Nz7v%K!=duRFzKnZTO<ni4S@1
zYt&_BucJL$_6<~nSM!--G>NBbe9smcM1wr_^_%KL#B<_6B|Xg68cpCkill4PG=j+@
z!k+K`lLYKx^l@E*<Q&(f)^K5_#U}tU6O@{lt%Ytcq4isg$b{ueXoE_zldy}joo#0~
zOoZR-=b{SJWf?Z>+D0sPFdU=i9jCD}b(g~DDt!49R>6KFB1D<ZQ7N8Z7cLUp;M)*s
zm16M+*zF4b+>1uXU{O;3P#d*C@s_9*!<d%{zXwDU>-MIl8DANg@#HS23SPK)oc6-S
zq>fHIzc^D)CFsa)aGq$N{?zH}Sn)3zkpYB;-NJt2B}@T{T_;aQmWJupdE)*gA0@jM
zs(WFPeNRj$aHma5CGZuSlujVd>WZh}GQ~RT5h3sco0Ll6$2KXQz!Q0;8+I7AS&s_E
z`@tha;M!RUq9I|9PhhbkVea|o>>|SU&Ehkri0;aSVe6$4NRsj%OxTMaK>}e)Og+gh
zGJq!`0!jP-)>c#vl}$}MyJP~3X342%UB#;B5g~BhEZLsKfnvnt=~=41C{bgU=qvRK
zWUxl66eA#ka|-0!PC9iwy+E{)9%ahHMLW}U0_Ukyk7(BfK4~9><dapxn!*-3P5iC4
z(#Y1Yw)s>7p~Y8ahJLG9L`O{%_^>)wD5Ff-9rbb|t!TeGEcl5^r7RsT{3Bn!8>CD;
z7Wv(nEIw;n<W$0n?eeD+=*gN3-kLS1da{<4!S8a`*1blpC~|^|=T%Y_j2IHRp@Jgr
zqlDlnc6qcB;jef=#f0Hch5!Fqww`pfyQ8WyIKInWq9XW)Cz8G~6LOxAS<!op$gaM|
z5GXQr*fy2wYYI|chWJWBJe?pXuWoP1N2`sm)<2cSOOLx9KnKr1$rr^4ctFMC;n|k|
zUhuteUU5SSf6Q*?y!{Y1Oi?lF5J>G+PMoMQQXC}p57iIQqQn%3GU=wA=;m>%<zCEV
z>Ck$KTgnr@)h&k@I-{=m{x)v8_zWsjnJbtr)#*5R-fTIICj2&hF{2_r<oRr%#~Nkt
zIh7iCN#p-_M3J}&c>SDI`~=k=_>6`Ea&Nwv(O6|AhAS>y3Dne<ui;+?&eVR5|Fv0U
z5ec7lTaES^jRRIHgKJQbeW<LL6~0&DbYlFZ#Qo_^crR=HKEy#+LvCFab~>3|r$EMm
z@UUC?JP!WEuTz@8-}oT71)aqzQ7hi1Dy74A<qF4nC;m5ZThO^DI=B?84hYdXD9_n4
zAs<Qd70S-<D#fuYZZEQt(r>O(&w(ky@3c^zr+7-GmP%h#=~k7}sg{=&{%W?g|1bKY
zJWjRw-A|=-UW3j?R4BSWRobJ!Ow?D=={&@c9O*BeEuc;H^p}v<D*5#)E%(?Jtqp&6
zN%%wO4`|Q-mVcLzoUQnD4gh|7+lwUnZS0Bo|C4{nSMbx6z5Xio+JZ0mz8mk8)msAW
z&K7+7i#`XGKJAFdHh0eNM7zsrTQyVqPL<Lg_XLk!Q9k?yBCg+pcA(R)?;s`j5Au<|
z8_7Mp$fdu&JH);0iSqoDe~=S=-(9j})ZWT9V<g{qmni6>Jmre-yUj6PZCT40EA@T1
z<Y|fzIg$U@3h$&gba{56Nj~jD3sU@6Ds3}P@H^4pObA!|DB@H~`(S9N2<>>fTk)}T
zMC4C9NO~#yqbkL25y8I=U6dbtJUG7v?e(Ay6q?%Y5bv|=0h@e;JpEMz+TZa%<Rd-%
zLuA~4QFfKRNZ-96UiO6l<QM&kfbFV$FRSv=Mi2T!0B02cyyDZwkL!x=r2PXEWqD{P
zNHaxmrP8)4ZLd=75aIrGx(!{#v!9Y1q*B`WLGSM`XUhEO+x;XzT**JKQqN`%As<N}
z8s4Y)UYj>WebAl>+Im4BET_#!87kg|Du3EXLE9ewm%Q&5i2orU`9JtAXbVJydgmXj
z(ncz6s#5G^5b?eZUBvrACHJsON2#<-mG_t`FMS!F<WrSgKb6|M7-W7eC~n&o-*=C~
zF-6axEbE89<^GMLdv+N}KE1Y|q+Z6;n|`I=hF;dwn{oQe_76(_ca?hjg%x<xh^|E+
z6Q>vU@tKnD`;uPKBmc!NkFz!LbE89-?rO4-tkOOz%~0uRm6i>Vo%ygzGZh_QCh&XT
zc+vjLw5t@oSfwwjlqNxY)g+`$rSyj@>45alims_iMqQQ0sx(!lbfos6#H|`%81;M6
z8CCiNi*(G6{+8agR9Oi6VspN_ldV%}u}bkyKlGT<E|Ol)(%VUTHHjB`Jc5zFZ!d|D
zmy3c=AO55zb6O(DnzrEIhAvuqXT=|=7TYpZN(<dn6d$ROr`1_njm7e*uzwr6kYA?w
z^r^^KlpW8Ss?hILrs!|0^dpt-Q0sc-Dy0=VT6+6>rp%vK_5Oixe*%$~=IDck{~#af
zyV~cQpYJlBmpvX2>)|s-d_cLXWV|?g@{F;!BYmAFT`f0dsC2AK>H9P@6rQirg?ZMr
zv&QO&=mT!_H1?o++Ms9nLiN;Au2OomJgJ`UuBjAHd)&!~Xbs^1k-bjz1Wr$n^!WXM
z=}!xKg#IUeAO9ms?l%65MoAwmo1o(4dF&MBZ&9G=!&JN`t91IUaTJyM*LhSbt%TxH
z8N<=#=Jd6od(^WWeJ7}g!jo0nufY1|eLON|C>(2R(=HgJtZQEzM(Mc=Mjw2R$ZGwS
z5z1nDQs?eqt$4BQbfPDzABR>;i(h&>dTp&n>GH1(w?psfuxfo{OyVD-p7V`ylDVw@
zmByQT@nS3DiqVoKls<UHc%1Q19lr7(;7;qpe~ec6-cInhMo<2=sR7>_xox(nuh0ak
zjEns6`FySBEq|A;MA6H0E&Zy|lK+9<Hdl=v5p*JAxY8ft(zN<okXA9zS^+-3*<_tY
z+L}5=H6`~O<l>NzJuEAFZGWksVWoU$g!0dRkNXbU;fs*ZeuwO6K&HPH{<#t_Uu6Bk
z<>Irf=HD9!_2jJ5|9o$ZX8O?uR-bF|qYMVFfpKbqHQ-0!_z3exq!D!FqL#{kGF;+*
z;Y(opbwsUSY3u7o4AZZTv-;csu9X(tFxun$<OSA&pNy6fw64kjd^^9^!4TKS$Eft5
zjqmls0_*zE#*j8;ixLE%;yf`#$*J3!I5Gi=1)%fz?VD%K`GpsJp7qKv#*ka`nWf8m
zDj#b#{}tuIA7`HOE9A?jly3YL1xG#m|85NBUoso_yOE2piCy^JXr<#*TDAT#a`nnZ
z*5W_7TCR2C5939!XZ+V_Wq<FCyqOy*y@|IDfsH>9`y&z_ySP7PrS$m|M#jvyo(0~9
z1~OEwxLA}dBi$8&iYrXf@%hg`{xpX0&wmd13-<6?wZ(tIUhEv}1jq4RO#N@9*A`f9
z{x+7u$wPk|p%F(_lMPemw#QH?O0WNIL^r}m--8?aweok)Ne~iF>%oS8?So?#8aqc&
zwR~$iqynh-q1(LQ(6608^}2Y=zmZ>i{~FIoBQ;7NY~*)DXHiygQ@{4Cwe?_AzX$x+
z?{o@d3D$B>J!iegsR=b@{4Q9*&HUQy^<T4kHS@b)e>TQi%+oO!t@n|#GRwcYUz7g>
D54<`v
diff --git a/analysis/src/VKBaseData.cpp b/analysis/src/VKBaseData.cpp
index f7e1a81..e0f4260 100644
--- a/analysis/src/VKBaseData.cpp
+++ b/analysis/src/VKBaseData.cpp
@@ -12,13 +12,13 @@ VKBaseData::VKBaseData()
// coefficients of the energy-resolution function
//Old parameters:
- //fP0 = 30.44;
- //fP1 = 17.30;
- //fTauRel = 0;
+ fP0 = 30.44;
+ fP1 = 17.30;
+ fTauRel = 0;
- fP0 = 0.0;
- fP1 = 5.5;
- fTauRel = 1.93;
+// fP0 = 0.0;
+// fP1 = 20;
+// fTauRel = 0; //1.93;
double emin = -3000;
double emax = 100;
@@ -122,18 +122,26 @@ void VKBaseData::SampleResolution(Float_t& e)
double fwhm = sqrt(pow(fP0,2) + fP1 * e);
double sigma = fwhm / 2.35;
+
+
double tau = sigma * fTauRel;
- // set parameter values
- fRes->SetParameter(0, 0);
- fRes->SetParameter(1, sigma);
- if (fTauRel > 0) {
- fRes->SetParameter(0, tau); // shift function so mean is at zero
- fRes->SetParameter(2, tau);
+ if (tau > 0) {
+ // set parameter values
+ fRes->SetParameter(0, 0);
+ fRes->SetParameter(1, sigma);
+ if (fTauRel > 0) {
+ fRes->SetParameter(0, tau); // shift function so mean is at zero
+ fRes->SetParameter(2, tau);
+ }
+ if (fwhm > 50 && e > 50) {
+ e += fRes->GetRandom();
+ }
}
-
- if (fwhm > 50 && e > 50) {
- e += fRes->GetRandom();
+ else {
+ if (fwhm > 50 && e > 50) {
+ e += fRNG->Gaus(0, sigma);
+ }
}
if (e<=0) e = 0;
diff --git a/mac/allowed.mac b/mac/allowed.mac
index 9a577b8..c67ebaa 100644
--- a/mac/allowed.mac
+++ b/mac/allowed.mac
@@ -7,11 +7,12 @@
# source dimensions
/VK/source/diskShapedSource true
/VK/source/innerDiameter 0 cm
-/VK/source/outerDiameter 9 mm
+/VK/source/outerDiameter 15 mm
+/VK/source/catcherFoilPosition 0 0 -299.5 mm
# radiation
/gps/particle e-
-/gps/position 0 0 -299.5 mm
+/gps/position 0 5 -299.5 mm
/gps/ang/type iso
/gps/ang/mintheta 90 deg
/gps/ang/maxtheta 180 deg
@@ -34,17 +35,7 @@
# run simulations
-/VK/source/betaMinEnergy 0 MeV
-/VK/field/maxbfield 0.30 tesla
-/VK/output/openFile output/allowed_b30.root
-/run/beamOn 100000
-
-#/VK/source/betaMinEnergy 5.00 MeV
-#/VK/field/maxbfield 0.48 tesla
-#/VK/output/openFile output/allowed_cut_14mm_b48.root
-#/run/beamOn 1000000
-
-#/VK/source/betaMinEnergy 5.00 MeV
-#/VK/field/maxbfield 0.50 tesla
-#/VK/output/openFile output/allowed_cut_14mm_b50.root
-#/run/beamOn 500000
+/VK/source/betaMinEnergy 5.00 MeV
+/VK/field/maxbfield 0.48 tesla
+/VK/output/openFile output/allowed_cut5000_b48.root
+/run/beamOn 1000000
diff --git a/mac/betabgr.mac b/mac/betabgr.mac
new file mode 100644
index 0000000..8087e64
--- /dev/null
+++ b/mac/betabgr.mac
@@ -0,0 +1,67 @@
+
+/control/execute mac/init.mac
+
+# magnetic field strength
+/VK/field/maxbfield 0.48 tesla
+
+# source dimensions
+/VK/source/diskShapedSource true
+/VK/source/innerDiameter 0 cm
+/VK/source/outerDiameter 15 mm
+/VK/source/catcherFoilPosition 0 0 -299.5 mm
+
+# radiation
+/gps/particle e-
+/gps/position 0 5 -299.5 mm
+/gps/ang/type iso
+/gps/ang/mintheta 90 deg
+/gps/ang/maxtheta 180 deg
+
+# energy spectrum
+/VK/source/betaSpectrum true
+/VK/source/betaZ 6
+/VK/source/betaA 12
+
+# output
+/VK/output/zeroSuppression true
+/VK/storeTrajectoryData false
+
+# run simulations
+
+
+/VK/source/betaEndPoint 7.50 MeV
+/VK/source/betaMinEnergy 5.0 MeV
+/VK/source/betaMaxEnergy 7.50 MeV
+
+/VK/field/maxbfield 0.4801 tesla
+/VK/output/openFile output/betabgr7.5MeV_b48.root
+/run/beamOn 1000000
+
+/VK/field/maxbfield 0.5000 tesla
+/VK/output/openFile output/betabgr7.5MeV_b50.root
+/run/beamOn 1000000
+
+/VK/field/maxbfield 0.5603 tesla
+/VK/output/openFile output/betabgr7.5MeV_b56.root
+/run/beamOn 1000000
+
+
+/VK/source/betaEndPoint 6.0 MeV
+/VK/source/betaMinEnergy 5.0 MeV
+/VK/source/betaMaxEnergy 6.0 MeV
+
+/VK/field/maxbfield 0.4801 tesla
+/VK/output/openFile output/betabgr6MeV_b48.root
+/run/beamOn 1000000
+
+/VK/field/maxbfield 0.5000 tesla
+/VK/output/openFile output/betabgr6MeV_b50.root
+/run/beamOn 1000000
+
+/VK/field/maxbfield 0.5603 tesla
+/VK/output/openFile output/betabgr6MeV_b56.root
+/run/beamOn 1000000
+
+
+
+
diff --git a/mac/f20scan.mac b/mac/f20scan.mac
index 2333582..d052533 100644
--- a/mac/f20scan.mac
+++ b/mac/f20scan.mac
@@ -37,7 +37,33 @@
# run simulations
+/VK/field/maxbfield 0.1504 tesla
+/VK/output/openFile output/f20scan/15mm_y5mm/b1500.root
+/run/beamOn 5000000
+
/VK/field/maxbfield 0.2000 tesla
/VK/output/openFile output/f20scan/15mm_y5mm/b2000.root
-/run/beamOn 500000
+/run/beamOn 5000000
+
+/VK/field/maxbfield 0.2504 tesla
+/VK/output/openFile output/f20scan/15mm_y5mm/b2500.root
+/run/beamOn 5000000
+
+/VK/field/maxbfield 0.3000 tesla
+/VK/output/openFile output/f20scan/15mm_y5mm/b3000.root
+/run/beamOn 5000000
+
+/VK/field/maxbfield 0.3504 tesla
+/VK/output/openFile output/f20scan/15mm_y5mm/b3500.root
+/run/beamOn 5000000
+
+/VK/source/betaMinEnergy 2.5 MeV
+/VK/field/maxbfield 0.4000 tesla
+/VK/output/openFile output/f20scan/15mm_y5mm/b4000.root
+/run/beamOn 5000000
+
+/VK/source/betaMinEnergy 2.5 MeV
+/VK/field/maxbfield 0.4020 tesla
+/VK/output/openFile output/f20scan/15mm_y5mm/b4020.root
+/run/beamOn 5000000
diff --git a/mac/forbidden.mac b/mac/forbidden.mac
index 472ecd9..87f9d91 100644
--- a/mac/forbidden.mac
+++ b/mac/forbidden.mac
@@ -25,13 +25,6 @@
/VK/source/betaMinEnergy 0.00 MeV
/VK/source/betaMaxEnergy 7.025 MeV
-##/VK/source/betaShapeFuncA1 1.0
-##/VK/source/betaShapeFuncA2 80.
-##/VK/source/betaShapeFuncA3 -2.0
-##/VK/source/betaShapeFuncA4 -0.28
-##/VK/source/betaShapeFuncA5 -0.009
-##/VK/source/betaShapeFuncA6 0.0080
-
# output
/VK/output/zeroSuppression true
/VK/storeTrajectoryData false
@@ -40,6 +33,6 @@
# run simulation
/VK/source/betaMinEnergy 5.00 MeV
-/VK/field/maxbfield 0.50 tesla
-/VK/output/openFile output/forbidden_cut5000_b50_shape.root
+/VK/field/maxbfield 0.56 tesla
+/VK/output/openFile output/forbidden_cut5000_b56.root
/run/beamOn 1000000
diff --git a/mac/forbidden_idini.mac b/mac/forbidden_idini.mac
new file mode 100644
index 0000000..822c141
--- /dev/null
+++ b/mac/forbidden_idini.mac
@@ -0,0 +1,50 @@
+
+/control/execute mac/init.mac
+
+# magnetic field strength
+/VK/field/maxbfield 0.48 tesla
+
+# source dimensions
+/VK/source/diskShapedSource true
+/VK/source/innerDiameter 0 cm
+/VK/source/outerDiameter 15 mm
+/VK/source/catcherFoilPosition 0 0 -299.5 mm
+
+# radiation
+/gps/particle e-
+/gps/position 0 5 -299.5 mm
+/gps/ang/type iso
+/gps/ang/mintheta 90 deg
+/gps/ang/maxtheta 180 deg
+
+# energy spectrum
+/VK/source/betaSpectrum true
+/VK/source/betaZ 10
+/VK/source/betaA 20
+/VK/source/betaEndPoint 7.025 MeV
+/VK/source/betaMinEnergy 0.00 MeV
+/VK/source/betaMaxEnergy 7.025 MeV
+
+/VK/source/betaShapeFuncA1 1.00
+/VK/source/betaShapeFuncA2 -0.102
+/VK/source/betaShapeFuncA3 -0.0827
+/VK/source/betaShapeFuncA4 -0.0106
+/VK/source/betaShapeFuncA5 0.00274
+/VK/source/betaShapeFuncA6 0.0000416
+
+# output
+/VK/output/zeroSuppression true
+/VK/storeTrajectoryData false
+
+
+# run simulation
+
+/VK/source/betaMinEnergy 5.00 MeV
+/VK/field/maxbfield 0.48 tesla
+/VK/output/openFile output/forbidden_cut5000_b48_shape-idini.root
+/run/beamOn 1000000
+
+/VK/source/betaMinEnergy 5.00 MeV
+/VK/field/maxbfield 0.50 tesla
+/VK/output/openFile output/forbidden_cut5000_b50_shape-idini.root
+/run/beamOn 1000000
diff --git a/mac/forbidden_max.mac b/mac/forbidden_max.mac
new file mode 100644
index 0000000..c935c45
--- /dev/null
+++ b/mac/forbidden_max.mac
@@ -0,0 +1,45 @@
+
+/control/execute mac/init.mac
+
+# magnetic field strength
+/VK/field/maxbfield 0.48 tesla
+
+# source dimensions
+/VK/source/diskShapedSource true
+/VK/source/innerDiameter 0 cm
+/VK/source/outerDiameter 15 mm
+/VK/source/catcherFoilPosition 0 0 -299.5 mm
+
+# radiation
+/gps/particle e-
+/gps/position 0 5 -299.5 mm
+/gps/ang/type iso
+/gps/ang/mintheta 90 deg
+/gps/ang/maxtheta 180 deg
+
+# energy spectrum
+/VK/source/betaSpectrum true
+/VK/source/betaZ 10
+/VK/source/betaA 20
+/VK/source/betaEndPoint 7.025 MeV
+/VK/source/betaMinEnergy 0.00 MeV
+/VK/source/betaMaxEnergy 7.025 MeV
+
+/VK/source/betaShapeFuncA1 1.00
+/VK/source/betaShapeFuncA2 1.24
+/VK/source/betaShapeFuncA3 1.19
+/VK/source/betaShapeFuncA4 -0.0816
+/VK/source/betaShapeFuncA5 -0.0595
+/VK/source/betaShapeFuncA6 0.00998
+
+# output
+/VK/output/zeroSuppression true
+/VK/storeTrajectoryData false
+
+
+# run simulation
+
+/VK/source/betaMinEnergy 5.00 MeV
+/VK/field/maxbfield 0.48 tesla
+/VK/output/openFile output/forbidden_cut5000_b48_shape-max.root
+/run/beamOn 1000000
diff --git a/mac/forbidden_min.mac b/mac/forbidden_min.mac
new file mode 100644
index 0000000..13cfc9d
--- /dev/null
+++ b/mac/forbidden_min.mac
@@ -0,0 +1,45 @@
+
+/control/execute mac/init.mac
+
+# magnetic field strength
+/VK/field/maxbfield 0.48 tesla
+
+# source dimensions
+/VK/source/diskShapedSource true
+/VK/source/innerDiameter 0 cm
+/VK/source/outerDiameter 15 mm
+/VK/source/catcherFoilPosition 0 0 -299.5 mm
+
+# radiation
+/gps/particle e-
+/gps/position 0 5 -299.5 mm
+/gps/ang/type iso
+/gps/ang/mintheta 90 deg
+/gps/ang/maxtheta 180 deg
+
+# energy spectrum
+/VK/source/betaSpectrum true
+/VK/source/betaZ 10
+/VK/source/betaA 20
+/VK/source/betaEndPoint 7.025 MeV
+/VK/source/betaMinEnergy 0.00 MeV
+/VK/source/betaMaxEnergy 7.025 MeV
+
+/VK/source/betaShapeFuncA1 1.00
+/VK/source/betaShapeFuncA2 0.326
+/VK/source/betaShapeFuncA3 0.318
+/VK/source/betaShapeFuncA4 -0.0325
+/VK/source/betaShapeFuncA5 -0.0185
+/VK/source/betaShapeFuncA6 0.00337
+
+# output
+/VK/output/zeroSuppression true
+/VK/storeTrajectoryData false
+
+
+# run simulation
+
+/VK/source/betaMinEnergy 5.00 MeV
+/VK/field/maxbfield 0.48 tesla
+/VK/output/openFile output/forbidden_cut5000_b48_shape-min.root
+/run/beamOn 1000000
diff --git a/mac/forbidden_suhonen_IMSRG.mac b/mac/forbidden_suhonen_IMSRG.mac
new file mode 100644
index 0000000..501ceb3
--- /dev/null
+++ b/mac/forbidden_suhonen_IMSRG.mac
@@ -0,0 +1,50 @@
+
+/control/execute mac/init.mac
+
+# magnetic field strength
+/VK/field/maxbfield 0.48 tesla
+
+# source dimensions
+/VK/source/diskShapedSource true
+/VK/source/innerDiameter 0 cm
+/VK/source/outerDiameter 15 mm
+/VK/source/catcherFoilPosition 0 0 -299.5 mm
+
+# radiation
+/gps/particle e-
+/gps/position 0 5 -299.5 mm
+/gps/ang/type iso
+/gps/ang/mintheta 90 deg
+/gps/ang/maxtheta 180 deg
+
+# energy spectrum
+/VK/source/betaSpectrum true
+/VK/source/betaZ 10
+/VK/source/betaA 20
+/VK/source/betaEndPoint 7.025 MeV
+/VK/source/betaMinEnergy 0.00 MeV
+/VK/source/betaMaxEnergy 7.025 MeV
+
+/VK/source/betaShapeFuncA1 1.00
+/VK/source/betaShapeFuncA2 0.0348
+/VK/source/betaShapeFuncA3 -0.138
+/VK/source/betaShapeFuncA4 -0.00728
+/VK/source/betaShapeFuncA5 0.00429
+/VK/source/betaShapeFuncA6 -0.000163
+
+# output
+/VK/output/zeroSuppression true
+/VK/storeTrajectoryData false
+
+
+# run simulation
+
+/VK/source/betaMinEnergy 5.00 MeV
+/VK/field/maxbfield 0.48 tesla
+/VK/output/openFile output/forbidden_cut5000_b48_shape-suhonen-IMSRG.root
+/run/beamOn 1000000
+
+/VK/source/betaMinEnergy 5.00 MeV
+/VK/field/maxbfield 0.50 tesla
+/VK/output/openFile output/forbidden_cut5000_b50_shape-suhonen-IMSRG.root
+/run/beamOn 1000000
diff --git a/mac/forbidden_suhonen_USDB.mac b/mac/forbidden_suhonen_USDB.mac
new file mode 100644
index 0000000..ee9bcb9
--- /dev/null
+++ b/mac/forbidden_suhonen_USDB.mac
@@ -0,0 +1,50 @@
+
+/control/execute mac/init.mac
+
+# magnetic field strength
+/VK/field/maxbfield 0.48 tesla
+
+# source dimensions
+/VK/source/diskShapedSource true
+/VK/source/innerDiameter 0 cm
+/VK/source/outerDiameter 15 mm
+/VK/source/catcherFoilPosition 0 0 -299.5 mm
+
+# radiation
+/gps/particle e-
+/gps/position 0 5 -299.5 mm
+/gps/ang/type iso
+/gps/ang/mintheta 90 deg
+/gps/ang/maxtheta 180 deg
+
+# energy spectrum
+/VK/source/betaSpectrum true
+/VK/source/betaZ 10
+/VK/source/betaA 20
+/VK/source/betaEndPoint 7.025 MeV
+/VK/source/betaMinEnergy 0.00 MeV
+/VK/source/betaMaxEnergy 7.025 MeV
+
+/VK/source/betaShapeFuncA1 1.00
+/VK/source/betaShapeFuncA2 0.117
+/VK/source/betaShapeFuncA3 -0.164
+/VK/source/betaShapeFuncA4 -0.00411
+/VK/source/betaShapeFuncA5 0.00510
+/VK/source/betaShapeFuncA6 -0.000258
+
+# output
+/VK/output/zeroSuppression true
+/VK/storeTrajectoryData false
+
+
+# run simulation
+
+/VK/source/betaMinEnergy 5.00 MeV
+/VK/field/maxbfield 0.48 tesla
+/VK/output/openFile output/forbidden_cut5000_b48_shape-suhonen-USDB.root
+/run/beamOn 1000000
+
+/VK/source/betaMinEnergy 5.00 MeV
+/VK/field/maxbfield 0.50 tesla
+/VK/output/openFile output/forbidden_cut5000_b50_shape-suhonen-USDB.root
+/run/beamOn 1000000
diff --git a/mac/test.mac b/mac/test.mac
index 704fc81..9cbbfea 100644
--- a/mac/test.mac
+++ b/mac/test.mac
@@ -19,26 +19,32 @@
# energy spectrum
/VK/source/betaSpectrum true
-/VK/source/betaZ 10
-/VK/source/betaA 20
-/VK/source/betaEndPoint 7.025 MeV
-/VK/source/betaMinEnergy 0.00 MeV
-/VK/source/betaMaxEnergy 7.025 MeV
-
-#/VK/source/betaShapeFuncA1 1.0
-#/VK/source/betaShapeFuncA2 80.
-#/VK/source/betaShapeFuncA3 -2.0
-#/VK/source/betaShapeFuncA4 -0.28
-#/VK/source/betaShapeFuncA5 -0.009
-#/VK/source/betaShapeFuncA6 0.0080
+/VK/source/betaZ 6
+/VK/source/betaA 12
# output
-/VK/output/zeroSuppression false
-/VK/storeTrajectoryData true
+/VK/output/zeroSuppression true
+/VK/storeTrajectoryData false
+
+# run simulations
+
+
+/VK/source/betaEndPoint 6.75 MeV
+/VK/source/betaMinEnergy 5.0 MeV
+/VK/source/betaMaxEnergy 6.75 MeV
+
+/VK/field/maxbfield 0.4801 tesla
+/VK/output/openFile output/betabgr6.75MeV_b48.root
+/run/beamOn 1000000
+
+/VK/field/maxbfield 0.5000 tesla
+/VK/output/openFile output/betabgr6.75MeV_b50.root
+/run/beamOn 1000000
+
+/VK/field/maxbfield 0.5603 tesla
+/VK/output/openFile output/betabgr6.75MeV_b56.root
+/run/beamOn 1000000
+
-# run simulation
-/VK/field/maxbfield 0.48 tesla
-/VK/output/openFile output/test1.root
-/run/beamOn 10000
diff --git a/mac/test2.mac b/mac/test2.mac
new file mode 100644
index 0000000..7974a3c
--- /dev/null
+++ b/mac/test2.mac
@@ -0,0 +1,38 @@
+
+/control/execute mac/init.mac
+/VK/detector/volumes/betaDetector true
+/VK/detector/volumes/betaDetectorV2 false
+
+# magnetic field strength
+/VK/field/maxbfield 0.48 tesla
+
+# source dimensions
+/VK/source/diskShapedSource true
+/VK/source/innerDiameter 0 cm
+/VK/source/outerDiameter 15 mm
+/VK/source/catcherFoilPosition 0 0 -299.5 mm
+
+# radiation
+/gps/particle e-
+/gps/position 0 5 -299.5 mm
+/gps/ang/type iso
+/gps/ang/mintheta 90 deg
+/gps/ang/maxtheta 180 deg
+
+# energy spectrum
+/VK/source/betaSpectrum true
+/VK/source/betaZ 6
+/VK/source/betaA 12
+/VK/source/betaEndPoint 13.37 MeV
+/VK/source/betaMinEnergy 0.00 MeV
+/VK/source/betaMaxEnergy 13.37 MeV
+
+# output
+/VK/output/zeroSuppression true
+/VK/storeTrajectoryData false
+
+# run simulations
+
+/VK/field/maxbfield 0.4801 tesla
+/VK/output/openFile output/b12scan/15mm_y5mm_v1/b4800.root
+/run/beamOn 50000000
--
GitLab