From 078802df953817a5b420b9b81cbe07d86fc4bd19 Mon Sep 17 00:00:00 2001 From: nam Date: Tue, 28 Oct 2014 22:30:26 +0900 Subject: [PATCH] prog saved --- thesis/chapters/chap6_analysis.tex | 193 +++++++++++++++++++++-------- thesis/figs/al100_resp.pdf | Bin 24204 -> 24277 bytes thesis/thesis.bib | 19 +++ thesis/thesis.tex | 12 +- 4 files changed, 169 insertions(+), 55 deletions(-) diff --git a/thesis/chapters/chap6_analysis.tex b/thesis/chapters/chap6_analysis.tex index 3fa28dd..0bd4bae 100644 --- a/thesis/chapters/chap6_analysis.tex +++ b/thesis/chapters/chap6_analysis.tex @@ -337,13 +337,15 @@ The band of protons is then extracted by cut on likelihood probability calculated as: \begin{equation} p_{i} = \dfrac{1}{\sqrt{2\pi}\sigma_{\Delta E}} - e^{\frac{(\Delta E_{meas.} - \Delta E_i)^2} {2\sigma^2_{\Delta E}}} + \exp{\left[\dfrac{(\Delta E_{meas.} - \Delta E_i)^2} {2\sigma^2_{\Delta + E}}\right]} \end{equation} -where $\Delta E_{\textrm{meas.}}$ is measured energy deposition in the thin +where $\Delta E_{\textrm{meas.}}$ is energy deposition measured by the thin silicon detector by a certain proton at energy $E_i$, $\Delta E_i$ and $\sigma_{\Delta E}$ are the expected and standard deviation of the energy loss -caused by the proton calculated by MC. A cut value of $3\sigma_{\Delta E}$, or -$p_i \ge 0.011$, was used to extract protons (\cref{fig:al100_protons}). +caused by the proton calculated by MC study. A threshold is set to extract +protons at 0.011 (equivalent to $3\sigma_{\Delta E}$), the band of protons is +shown in (\cref{fig:al100_protons}). \begin{figure}[htb] \centering \includegraphics[width=0.47\textwidth]{figs/al100_protons} @@ -354,6 +356,12 @@ $p_i \ge 0.011$, was used to extract protons (\cref{fig:al100_protons}). \label{fig:al100_protons} \end{figure} +The cut efficiency in the energy range from \SIrange{2}{12}{\MeV} is estimated +by MC study. The fraction of protons that do not satisfy the probability cut +is 0.5\%. The number of other charged particles that are misidentified as +protons depends on the ratios between those species and protons. Assuming +a proton:deuteron:triton:alpha:muon ratio of 5:2:1:2:2, the number of +misidentified hits is 0.1\% of the number of protons. %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% \section{Proton emission rate from aluminium} \label{sec:proton_emission_rate_from_aluminium} @@ -365,25 +373,20 @@ of protons is normalised to the number of nuclear muon captures. \subsection{Number of protons emitted} \label{sub:number_of_protons_emitted} -From the particle identification above, number of protons having energy in the -range from \SIrange{2.2}{8.5}{\MeV} hitting the two arms are: +The numbers of protons in the energy range from \SIrange{2.2}{8.5}{\MeV} after +applying the probability cut are: \begin{align} - N_{\textrm{p meas. left}} = 1822 \pm 42.7\\ - N_{\textrm{p meas. right}} = 2373 \pm 48.7 + N_{\textrm{p meas. left}} = 1822\\% \pm 42.7\\ + N_{\textrm{p meas. right}} = 2373% \pm 48.7 \end{align} The right arm received significantly more protons than the left arm did, which -is expected as in \cref{sub:momentum_scan_for_the_100_} it is shown that -muons stopped off centre to the right arm. - -%%TODO -The uncertainties are statistical only. The systematic uncertainties due to -the cut on protons is estimated to be small compared to the statistical ones. +is expected as in \cref{sub:momentum_scan_for_the_100_} where it is shown that +muons stopped off-centred to the right arm. \subsection{Corrections for the number of protons} \label{sub:corrections_for_the_number_of_protons} The protons spectra observed by the silicon detectors have been modified by -the energy loss inside the target so correction (also called unfolding, or -reconstruction) is necessary. +the energy loss inside the target so correction (or unfolding) is necessary. The unfolding, essentially, is finding a response function that relates proton's true energy and measured value. This can be done in MC simulation by generating protons with a spatial distribution as close as possible to the real @@ -413,28 +416,65 @@ method is implemented. \caption{Response functions for the two silicon arms.} \label{fig:al100_resp_matrices} \end{figure} -After training the unfolding code is applied on the measured spectra from the -left and right arms. The unfolded proton spectra in \cref{fig:al100_unfold} -reasonably reflect the distribution of initial protons which is off-centred to -the right arm. The path length to the left arm is longer so less protons at -energy lower than \SI{5}{\MeV} could reach the detectors. The sharp low-energy -cut off on the right arm is consistent with the Coulomb barrier for protons, -which is \SI{4.1}{\MeV} for protons emitted from $^{27}$Mg. +%After training, the unfolding code is applied on the measured spectra from the +%left and right arms. The unfolded proton spectra in \cref{fig:al100_unfold} +%reasonably reflect the distribution of initial protons which is off-centred to +%the right arm. The path length to the left arm is longer so less protons at +%energy lower than \SI{5}{\MeV} could reach the detectors. The sharp low-energy +%cut off on the right arm is consistent with the Coulomb barrier for protons, +%which is \SI{4.1}{\MeV} for protons emitted from $^{27}$Mg. +The unfolded spectra using the two observed spectra at the two arms as input +are shown in \cref{fig:al100_unfold}. The two unfolded spectra generally agree +with each other, except for a few first and last bins. The discrepancy and +large uncertainties at the low energy region are because of only a small +number of protons with those energies could reach the detectors. The jump on +the right arm at around \SI{9}{\MeV} can be explained as the punch-through +protons were counted as the proton veto counters were not used in this +analysis. -Comparing the reconstructed spectra from \SIrange{5}{8}{\MeV}, the protons -yields are consistent with each other: +%Several studies were conducted to assess the performance of the unfolding +%code, including: +%\begin{itemize} + %\item stability against cut-off energy; + %\item comparison between the two arms; + %\item and unfolding of a MC-generated spectrum. +%\end{itemize} +The stability of the unfolding code is tested by varying the lower cut-off +energy of the input spectrum. \cref{fig:al100_cutoff_study} show that the +shapes of the unfolded spectra are stable. The lower cut-off energy of the +output increases as that of the input increases, and the shape is generally +unchanged after a few bins. +\begin{figure}[htb] + \centering + \includegraphics[width=0.85\textwidth]{figs/al100_cutoff_study} + \caption{Unfolded spectra with different cut-off energies.} + \label{fig:al100_cutoff_study} +\end{figure} +The proton yields calculated from observed spectra in two arms are compared in +\cref{fig:al100_integral_comparison} where the upper limit of the integrals +is fixed at \SI{8}{\MeV}, and the lower limit is varied in \SI{400}{\keV} step. +The difference is large at cut-off energies less than \SI{4}{\MeV} due to +large uncertainties at the first bins. Above \SI{4}{\MeV}, the two arms show +consistent numbers of protons. +\begin{figure}[htb] + \centering + \includegraphics[width=0.85\textwidth]{figs/al100_integral_comparison} + \caption{Proton yields calculated from two arms.} + \label{fig:al100_integral_comparison} +\end{figure} +The yields of protons from \SIrange{4}{8}{\MeV} are: \begin{align} - N_{\textrm{p reco. left}} &= (110.9 \pm 2.0)\times 10^3\\ - N_{\textrm{p reco. right}} &= (110.2 \pm 2.3)\times 10^3 + N_{\textrm{p unfold left}} &= (165.4 \pm 2.7)\times 10^3\\ + N_{\textrm{p unfold right}} &= (173.1 \pm 2.9)\times 10^3 \end{align} Therefore, the number of emitted protons is taken as average value: \begin{equation} - N_{\textrm{p reco.}} = (110.6 \pm 2.2) \times 10^3 + N_{\textrm{p unfold}} = (169.3 \pm 2.9) \times 10^3 \end{equation} \begin{figure}[htb] \centering - \includegraphics[width=0.85\textwidth]{figs/al100_unfold} + \includegraphics[width=0.85\textwidth]{figs/al100_unfolded_lr} \caption{Unfolded proton spectra from the 100-\si{\um} aluminium target.} \label{fig:al100_unfold} \end{figure} @@ -458,29 +498,84 @@ the number of nuclear captures are: N_{\mu \textrm{ stopped}} &= (1.57 \pm 0.05)\times 10^7\\ N_{\mu \textrm{ nucl. cap.}} &= (9.57\pm 0.31)\times 10^6 \end{align} -The proton emission rate in the range from \SIrange{5}{8}{\MeV} is therefore: + +\subsection{Proton emission rate} +\label{sub:proton_emission_rate} +The proton emission rate in the range from \SIrange{4}{8}{\MeV} is therefore: \begin{equation} - R_{\textrm{p}} = \frac{110.6\times 10^3}{9.57\times 10^6} = 1.16\times + R_{\textrm{p}} = \frac{169.3\times 10^3}{9.57\times 10^6} = 1.74\times 10^{-2} \label{eq:proton_rate_al} \end{equation} -%\subsection{Uncertainties of the emission rate} -%\label{sub:uncertainties_of_the_emission_rate} -%The uncertainties of the emission rate come from: -%\begin{itemize} - %\item uncertainties in the number of protons: - %\begin{itemize} - %\item statistical uncertainty of the measured spectra; - %\item systematic uncertainty due to misidentification; - %\item systematic uncertainty from the unfolding - %\end{itemize} - %\item uncertainties in the number of nuclear captures: - %\begin{itemize} - %\item statistical uncertainty of the number of X-rays; - %\item uncertainty of the detector acceptance; - %\item uncertainty from the corrections: random summing and transistor - %reset amplifier - %\end{itemize} -%\end{itemize} +The total proton emission rate can be estimated by assuming a spectrum shape +with the same parameterisation as in \eqref{eqn:EH_pdf}. The fit parameters +are shown in . With such parameterisation, the integration in +range from \SIrange{4}{8}{\MeV} is 51\% of the total number of protons. The +total proton emission rate is therefore $3.5\times 10^{-2}$. +\subsection{Uncertainties of the emission rate} +\label{sub:uncertainties_of_the_emission_rate} +The uncertainties of the emission rate come from: +\begin{itemize} + \item uncertainties in the number of nuclear captures: these were discussed + in \cref{sub:number_of_stopped_muons_from_the_number_of_x_rays}; + \item uncertainties in the number of protons: + \begin{itemize} + \item statistical uncertainties of the measured spectra which are + propagated during the unfolding process; + \item systematic uncertainties due to misidentification: this number is + small compared to other uncertainties as discussed in + \cref{sub:event_selection_for_the_passive_targets}; + \item systematic uncertainty from the unfolding + \end{itemize} +\end{itemize} +The last item is studied by MC method using the parameterisation in +\cref{sub:proton_emission_rate}: +\begin{itemize} + \item protons with energy distribution obeying the parameterisation are + generated inside the target. The spatial distribution is the same as that + of in \cref{sub:corrections_for_the_number_of_protons}. MC truth including + initial energies and positions are recorded; + \item the number of protons reaching the silicon detectors are counted, + the proton yield is set to be the same as the measured yield to make the + statistical uncertainties comparable; + \item the unfolding is applied on the observed proton spectra, and the + results are compared with the MC truth. +\end{itemize} +\begin{figure}[htb] + \centering + \includegraphics[width=0.48\textwidth]{figs/al100_MCvsUnfold} + \includegraphics[width=0.48\textwidth]{figs/al100_unfold_truth_ratio} + \caption{Comparison between an unfolded spectrum and MC truth: spectra + (left), and yields (right). The ratio is defined as $\textrm{(Unfold - MC + truth)/(MC truth)}$} + \label{fig:al100_MCvsUnfold} +\end{figure} +\Cref{fig:al100_MCvsUnfold} shows that the yield obtained after unfolding is +in agreement with that from the MC truth. The difference is less than 5\% if +the integration is taken in the range from \SIrange{4}{8}{\MeV}. Therefore +a systematic uncertainty of 5\% is assigned for the unfolding routine. + +A summary of uncertainties in the measurement of proton emission rate is +presented in \cref{tab:al100_uncertainties_all}. +\begin{table}[htb] + \begin{center} + \begin{tabular}{@{}ll@{}} + \toprule + \textbf{Item}& \textbf{Value} \\ + \midrule + Number of muons & 3.2\% \\ + Statistical from measured spectra & 1.6\% \\ + Systematic from unfolding & 5.0\% \\ + Systematic from PID & \textless0.5\% \\ + \midrule + Total & 6.1\%\\ + \bottomrule + \end{tabular} + \end{center} + \caption{Uncertainties of the proton emission rate.} + \label{tab:al100_uncertainties_all} +\end{table} + +The proton emission rate is then $(3.5 \pm 0.2)$\%. diff --git a/thesis/figs/al100_resp.pdf b/thesis/figs/al100_resp.pdf index d4728d2963b8be8c35e581e9652243e89370ee6f..963ea3a9f81834ae5b05954b6b97613c1d7a7be1 100644 GIT binary patch delta 12129 zcmZviby$;M*!UGeVW8wtkrD(3qmdZh9UC!XgrszLJ<&)#(0p6a-rSV71y!Y9Zl%=^5ywhn?P1WU1I z4+Kn5V^|^gOXwOO(F@E@dY@w={%t(`pMuFzjHIJG@0P737T&}eUk7C=nQ+-l*2XM; z9o`AM`r2z#-&U!D5Y;}h-M*PG&Tm$d(nwgmI=wu|-uby8dE@+cV^8XP#`R5h*AZ}b zc(ttGby4jwbG&yomSA=@WwvCdQ)8ExjRAyXm%ktWIX^x8y%*Kxt3Dzwh01nQ!(2@W zn@ByecZ_WEnS3ue|6>NwV7+_r?oXE67Q?)d#7%3l_wI1^mDzB1OmJ6YE6E@6u_2G- z3ZH8qo6q9*7jC@zYkaLX+zC;!meZVn7!j5N$My1D8e%TlvQ8rtMAtM!XTrP( zf_ODP%G1zbrq{v}KT?L{e3Q-gCUgK`Y{1XHtdvqQR7$Siv3`YA0i1CgxkvtyAs!Ea zaGZ7{Uf?PDmGypBV9g<#By=B@$CR%)u`ueDywVcZmH|`a*FP2FhH6SmhNl z9neM$d!)Y13r#d6my;nK6K5C@7?bH2SKxB0w?|E9l74?T(th0Q7{WkH5DdUmw^Osd zowk9_G%4W~2U3RnBt_^KWJ)12v{;g*Z*8Wxa`YRZzMo`9pFF^G39riO>xSA}Kh9z= zjW7OK-&$XV*=XE+c5ZJ9lQP@yhWX9jd&$*&=Q}919$Dc1-97O@=aqLbXnb@`--6`( ze&>r)8X4DFL+F>RxW-Qo^+dpzS>*%&FuoU3&2lI9%|W*=Dsq1{(PHxRTcf*|cDL+b zp@?=Ca~@C1j7Zbbu5A8d)bAs1CGi^XEPN1;cM-uLD+j;F;2?NCS~CKB{#@8;TEoJp z2eE6waB4e_rzAr>QO>~yhZYjVyO*g&mGSwxickRz7a0Ndq1w)6#701ojjUkcd$r?; zx6Yz8_5r>3Y}F6Jsu6c9(i#P>U1FFD2O1xj?*-KpmF{;Qsg`=nR?rfZYU@G}{+dT5 zrK(W%aQHr#eu0%ezG?@S+jCGURMSO+gDN1!tnkM}5&D?YAh{$9{OUh&&c?cUjkZG& zfl2uz=U5{LSj?_sFCkEE#3y@{M%{eJ@x#I_I)py~3go8WoZQbfij!ng}W>6g6{!5JFe? z>5ftmc&XSgEh$XT__Zyw?uBVYNm3Sve4Kn5xty0x;B4y!Eysme4S%x@o!;vAhT<=C z2Lytv)Q^jL=2id~-6!^Jl>@|;aWu^5`~FQwacAO$P(;O>0BSD13K39ROfft~1SI&n zsE8Fo1#~LC`Q;T&3O9^GlTSW|X=~*EW-c6!F68(_(n#U8^1Y!be(vCrm;8B8C1R4? zo~1IJR!#TGB)KFFN94GSb#DQ(>1&o!gsh3-?4fQ$B^1EJdMOmNR}-dAgX3jQenU+4 zKNAL+kl#`E_qG0Q(_Xh;YkPkd8vAGf8cScfD*75FV@(dJp&cyNW&d7m`_8YUMfk}i zu~P(+u7+(6RM~4{RG8-~hJL8ZQncwSL2k*Vc&JcVv`M0Ia=AfJBT5@l`3ZskH^?Cn zULqjiD5AK%6HZonkKpTCyn?<-S9_(hK*RIGO*j;x_u0ULG^CPQqWb3;VFChh_M<&o zd)tQqO29_F@O)>iSS4j6#nc9Yx01#Za{=pXa`~&-EGOHvG_)gYhFx-6TrvD7E`3tm<9nw#EF2s)(BJ;1##|XqRD+nHfI)2Yu?F^7$sLf)T_tkDYh6M z2jQqBS1P1IP=W5>^DF$#$iq$>o-B*FEEM>vzbGQhp5$$7D94_Jw3U$sDhiQBsnNBz z?7x`{QeCN3doeI2+>*i17;QO>Jfd|V6C%d*S12ZKy-M}F)kJ0lIQ1Q-aaA*oGh{!R z1>QNb1fgOMnQ2y}UGY+U)jsM<#Q{Xu`yYWT2Rdj|d1JRd0SIm|6NX5I=|sw4o*@Sl%)Nl8BWdCqv3VmQ6yx z-}ieQGa+-KX!i@Ee1lXkdaew0)S-z8F?vn$7tcZNw5SHc0Oamf>b1s>3{R3=F~yQ7 z{LXq^47U>L|8tK(!TJ|u@>aa;{Q}6fIA*=-&%IQh?eP>Xh?!9ZtQL;^+><~9H^*vR zNn5w9-R9_n(ifQaE7w(09xWk}#d+-?$IxuJTL3X%Xr*Y_fX0AWYTNUn3015+RjL}I z5D_1W$6g5s2nxg&n)f>+nzTL6YTm);^$r_16klGP7`rt1Kw^tzVc@(}s5E2^DK972 zT=#?~1^Ky`NSD>OvLTE90tG#qYIORr@bH`hGMz5nnk*6{^^pBMMqEZ}Ms;K_v(CLj zcA9*bm~QF~@0Me{`Q@%-gnzFVx~gn<*8aBM)FCSX^3j_gf;~S1{V8}hrN80kG~wTe zdBKxcC#1E>eNXAEQT!fVQ7M{?OCdu1@n#}}oL&wj>Q)*Z$8sW5b5B#EU{KAtw*@$;n!q)~i2YsS3+~;j)2^TUyc{;gJ3c|8;fUcfNU$905Ce}mFq|=eXw~fi zn3tA-Ue-r{`qBE}#Kmz1SaYoY8)Db<^Lj@8<9^J7I&{=Mld>ATc3~3vYL7`wR#@)k z6}36RiDa9zPm|(kVsf8YPt2$hPic@u3o%%( zm)^3U13X*0ri8i%<_%u^vAz%@*qPs!%mFBY2gy?2Od^RfWjvXu3SD4wlzmowv@@A= z$;N1MBd*82Pcvr)aSM2Lf{_B`sl2Ibx|N+Gc0)d@&@}&gmXQ5Rs-?n&Rhz`Q^7&uA zTf^;{!@pQb-&nMmxDOw#G~Ui-C2nqePXTvr^2tyAY!k00gp%%p3W_ZQ;#lS$Y5S7Je?LE`c4Ez zT)(+%so&1?Sl1E~{OWTXp!zn{c1{)R&tAsIK%2cd5qp?AOjU@csf4MN&;^(MdiER10Zri5w#Z76Kcvpq8dnG} z*1gJM<;TPmHE4eGwlGEoXqpmlvzf&#A{+k=gPJTV9N|(e_JjbFhaW zdfc;jxLaRjh>r({CYej_4wdp!R9Sr42pQhbU=McGSoE=Pl1Xed>|=_PCb-K{ONd_( z#1bccYbT^qf$>(ew!wAiqxKW5mmFD`k_?!EOs+Q>_(uHs*8$m~QJYP!k_zqEwPP)6 z@obrACFcY)uMdLmvvXSI&m^%BivL0;uFXmYlv(3f94Bt*ySK!TOmf%Nd>F>3#5JbA zTMU^t5I=iLzg48?R+624xq-vSTcVZ&LMu0`#-X=4>U3Xs6VN9)&~s^&hZ*|0*?Mjc zE{3NA_%ar)X@lCDm=ehYPfgd$AXVN*ij6j`TZxuJLdr)p!gj|`?l8E-dIsB~ialHY~)Wrlj{K&N$ok_Zmhxv=QHb8hfB zi+Xe&z!r}OhEy)`5?PUX>xSl#%)L}$XCrO~maKj~sLB=VXZ)nZPGCf6R?RN>Gh%}E z!deQ#1ozWZGplO;I7w1N?-S_I$+ujPbC7Fn7G(L>ka%TWO ztLRwi6Qr2(D|ol^#zYgcoMNXnhH(C0+|6ZkPE|;gc-f z-%>Y$S(_jS^7%(CKRWa!xzXS6yby7}xeh$UY)dvHj`hSIT6x(oTu@{vbNt~LTuags zy`EyfOj4fEbKC4^rZWSRm+UD%ix2D+qYO%CZA4ec5Gk+Ccpiwka4S;VEzSU?y1uR* zH6r{2{#9|kNvbPs*Qk}QW8DuPuG2PUS`rbodY^GeYrT(YqO^o zq7?i{)dsdNhw0~JMh(Ys6ZX@bl%{&D#Pu!nuO$gf?t6{)*=&?fzwjlFyqW%V>~?S* zX!c?I@wYL$^v(_cp>+|Vg@b$s1{2D8OASPR+1jgbA`-zT%gey@qGa%q@P>)j#`oUi z(__*Esf&@Q{F}pXefd)2!IwMfT2hz0%h`R4!09I*;QI7fiXOPyUS?S4=(XrLk~=xw z`>NHnfr_#{TV8J7_2KTA$1FxhnUJ`ws;r`>-)uLD2VJx!NL}uK@&N)bH@<0gT^{c3 z{a%jJTD)0YHUq8?K1p5uLEYnG@pxIalHg*j`5<1W{>{J!SCEFsbKM3jNPF9fw1Qf^ zQMQvNs`gIATM{e-RCI1Xk~YPnp@c^%ywr98ajj5Hf?H|S9QixfPQYQrPt|Gb*%I-V z(%xxRh=tElK`pXpsGC%`VFfa@K9c|gaj{KaYktJ-s38Y_HQ6ASm;cYolL4F+%E6}3 z-v)SaqBj`7)%4~ZcuK9e%CH(;6g{kOD&*{t7F%>F5)+WZqhcz;bq&D@NV2LsSg z1uRNVCf~eZ%f;C68B#HU{B{HaQL1~&6vMBkHkF=jtSbJqSMW|0NveLzVM??#u@Wu8ljFwRa}ZF z2d`&eWY3dexnP5eJALen!G*oeW?Eg2dbC>>;&O#7&4OG_-ZRy_asY0xS=)7NT(S55 z<-~@03})-n-(fVNw&Zl-{c1cJrAkWXY(lSt%ahs%?lH@lx3rBvL6f6?>R*#DBD6FZ zAPp-`?@$7V^F)ofR9)AW;=I00C9=%WLU&+S(z_O5a(+VtNF$l@z<66&@$thvEN&sZ z4*0n$UTF)!&ak@z{fK~}fwvg*`PXB(>|8b*JRU@iLS}C{vURoNY)Exx6}Dro9QbU! z|9ufYldsS(9e~SL>HXsYs0+*w%QgtJcNuUhoi`L)y;FpqKk1XRdlNvQMR|*%j0vckYzE<0I4aO6Oac)E|{|G9jOCq)4TwgtM z*p3b)qrDYaN{0}2_3w0LGieAey`Mx{SCVX;)jxklc2|Xzys1qh|L0 z!0dyj_b$tBFw5Au?YnP`z4az%RW%d(+iSS+58mzsnVbZnmE6tyY4?oJvnxX`&<6RNF|Q=jCs+k7 zUK`!*I2ntn@BAeBCpOcIXE50IAwvvuI}#K-SWKf$b_hYNUv>5C;$ zqq9p%s&m6AJ7$)(`fpdkny(k@h?{K(nxulxE50?D7HYm!6XTd?+KS@Cm6Z}E$Ua+u z%E(>!{A%_4TGP0@>`$+wmk<|NA?^odDy}|I?%sol1?>d?o9>T})?9SJ1kkw6q^P{x zJJW>W?KxzLd|wr1of9LL62JzL{BmIFgtC*H2}<%SA@1>puoU0$ANbalH=JyA4U2Cv zIp(QeZ2n$I+IZOdb6O%0+IXW@spc?489#i#>Sfpzx>Z86<-*fFEjT$d&*jO`$#dxW zRnvPVdHsXp?Y@bFX}NshK_N0^LQG|LQ{0Yw%Q!&v$4;ND?`)L0ttg7R_f6CQH7Ajg z(?}|}iPp<48l4W3_VT{w1ttkRZ1)-@JhNzl+=g$4-sbi`uLl+0ARW)MPfbpugL%ZL zn6@RcQKpAfM-M9geof>G**7Mf%p<(JdlHP7g`(2^+&46b>i%xpZyVDZy_ zqp(WJwV4O^t38&Qgh#X57Q_u%@Ydf^KwielO2FI#jz>w&%RgeaTfb=s^^!jMD3LvO zJZjcXz-+hp>j^_)DEp|0B+cnFrK;a4j|R^^BYvewZwFE^zKoTI@t6*vlk7}|G7)QU z80ud8v%p6m#Q5J615RyUL`iS|Oy_n3!Nf|y-V~3wem@Wk!0-R*;bQdSAp^1ak|kKg zRv+2{0$g34R@V;fa>DSiv?~zaH$@AzODbZ5VAYl_+-FyH(a)J5Jm4f)+ElsPu%lnU zItW9H=PchRCFs4VGvOHHU~PmRIvV;*?kQ)EB1*ji*Sa17-DfSUEH;$(rCBekrkz?U$g<#J{R?nSUY;Y08c>1~e02X5W1s-NN26fkNNwE_xBYF&31Q}}>Jf(egLtH6nd zl-*nX%~8N5jd-3Y^-p|J_|7TeI2}Wwk&gxMbQ|Z{TSzQhHDJ zFH59a;m*5$B=f3Jy${%fNDi~FAH}nRpZ&H6J8k90m;odT#HrA%;sTQp_ib(63ehcp z@SPS@7|4IcS($y?pGvNE0e|r58Ql`CPsd(u0UN;gR=f0%${VbI(m_#WSl*7py}@MU zWLDdwZ-P;Os6I5T>h3k}!O^nYATUJ#-=K#sr-q@QEt?ct+|oxsDm?B?pmChWt!5axcovB;RzqXja z;3;t`v0;JtQ{qImw?(Jlh%NBbD{i%_w>|+2c1Lh-PtuMM=0CN9A=45g(o~;%Bs{eW zasQ`tw6JAT=JTzygiVbl$R4{x8czJwgDRX8&TSPk6TvCXZ51DCMPQ%KNwCMpF?owy zM>2F&EPPNY434#JhqS;;!JSI3Kp_Y!yF?H;C7dwF_~#6w1o;P zNE^25mZ_Nyd#4AF+;Rb$3HPQCSp^4sZzi6Z{N?#-Kfl;kaE_!VMs4Lw-7XBc9TMtI zU_V`btJNDI7;GFuj8{vysqpM#BOl+h*V)Zq?C*9_Km4!;SvM{4o+nQHSA@19?znPv z&@V7;fPD~CwA)KT^0Brj|Z7 z((aUtQOtVfRkd#O%pf-Z>6f4F(r0Js(q~3h@C)3VJ0L!^fMAe-3(jp4?{2OQFBFoF zZ+B5}a+OgP3w*d0y7DdJ?`H#DUoa3<^K3Vi>bmup=^iP6TXuqzZ|A9wn`}RLX>U4% zsE0U}(zGDwFD*8_oKeP?@71@=1jcF6^7;nZY0+^LgNnAS%G!P8bIpsGPW>V`*x%@| z)0)$Jb&Esag@o(f>-|Ge69L=U(Zu2$EZ#~WHDpadA&qL8^Iu!FV91Ga1-L|FOTiz8 zUyb#|{1QkIIW^|{^8SK)`E(IhZ+0;GIEcq>)@I6tTz#|WkErR0#^V}kv1PFR%XZoA zxm8GTm+~P-68P@j!tbpKg?dec1d*z!^sLUUxZ+>!{d1$+)3o!CF?U8i`IYE zKMa2Hqw(qh<8!$$Ur<2*00{LL0d4+j=qJS4!&ATd9^BB7edsr%|IaTnmH@pgGU|)BUuR@m zilWJkrHTr}9cZU3M~O5-#v&fi90m%vLWe#JjmGt`0_c7WL{^NwF=_u={=_;}a=m|x zHaxJa{2uu(exRf9Y%=6rS@xb@bLw3cy?8u57)ZAi&B!8`MtGlNqSNc%lfmdE0p$uz zdnr*g`l4OB(oHP`zNkI;tX6~}O*#iFIY_jm;M#Hl*j#nMh8?i4z@ zD!0VE{H6yfNt_~_X2U-B3?9+5^rODE#k5aqCN^6`!;%7kUmDDoZZHjyXT zoXN6JJa>HAcLA|Ow5G-!svJe5vU z3?79-2$FvOisaZNX=d73YiL~q*%pVY)g+ux2?0^{IdE)HTrl6=)ZNdo0%-~M5xVf$ zLVF0von7pbgK}jF3DRGa;h1l$q(KNIXBU8*nUF^XX&|WH^$Sp*@J6grq6~WEt?>XG z2+xoZ{XUK)@_lvV69?xN_bUQn^jhO&fueop(>3W<`pFi`F>~MGHWPK_apICkA8eq( zz#?>C3qi2MxZUo8T}3`(PlIR$DY%#xiC-|0ghnf0o46$21ccDp%=&4#Yi_0n-Avu? zL~wGDDG``Is+1QhC)(!l0T&QAGc`HNP)o*Dk#aLo6^4ND{*ue95aXdU+-NXOQQG+E zPD4vcjLu%TBw-Ii^MB!_^#Fi(8~3P$fQRs7hkgZtObb}AMmjN|^A%&Wg{b$j4_gp< ztSOnT!|_6@G%bE5ukLFbM#Al?N;2K>N%0){c#H-s|A@?_A7V!@6NdtJW^4Fb4=F^WC zF{NCm1xtg5je^C`wcA4Y({ZGlmqGp)q$IP@emzv)LoAd#$;K(-8v%Q{H;tN1$?VNl zs{MCVIfUd9Pog!^(Rz6vz($r{lT(kkmjo`^Yu zlR2d51iKmD@R~8Y2gVggZ&o6ioJC*HZbJ`aE}$`Cw(B7Q#@g53%p^J}JWI$(>XoeP z4`e5Mp*5c?JuFo01~0BUWUluRXzG<)Ql^!gBKzo3KKwU1YSImwh1XACUkrQCmK@2$ z@8z&ie}sBzUF!V=W;npm_zMPV_yy4ty`$&KHZi)5PJUhD0L>8(n8^gILeE;D??}Itd-rXX))-|-m5Q)w& zrIokGnE6|T(51b1{e1Pw=fmG%eS0VDed06(c?!-)>dJ+NcTouBh8$MQc6BO{gZ}*z zA8ovQBwXD!q)=5H?~bjRq20$@&I)^kBrlLw7(!MEba~6r1G;zen#f=Ajh0S*X-7!^ zn-4|6zQ*2!;I^@N$b4oGU%!w-TE-G=HORVmkCqtf?h(eRBUZP{I5Yuzne(#MLh}u@ z%1+f9?#@#8{JG9W=LsiWYhjQcSM;zew1`^XCL`^97y@mlW)$GyPCavY`S4@x@S8y% zcPoK_HK=mJ@OUSGS3gN}(SGbOFI2Rv|4khrTP1C!I2OI`Lzi{v7^VD1*`>tRDbhzW ziptO2gkcFHF(Dg=o@0W3Ru0mRPKs_lJRZ@=^>XI4;27sX(eF>ad$1cg0Zv*CDt#m| z*%PDu$U+07%;(QvLZ0V67E3~0)0l0NKgQ3;(V?7FdY0?bB|#x~!2wgWXaT30YLNf} z+WC*HBe+YG9!V4mC(nzvl`L3^>S`xtRIp#uTt$oBR&iq5npDhM=#^V1cbYIPC+VhE zwEe=g=}u$iH21{@_n#;2fP8b4Y-jGy?}LjaBM7BX*q}haw9zMJmgb`8Z;bh&b@EF5 zX<$tvE1j`uHw(KNL}~rG|4V_L=H4fOVxIEteR+g#uf(kq)gRzs&#&9(v z=g*p7GpaNDBtXy9%A5JEj)1cx7EJR7?JicD!pOuE>TMDiTH#Q;p~KMwpq71x^KCrE9i2$wrFJJ zd-ic~pt{X_jhE{S`+a>bi@PTN+q(XZQ$L2YTa#haKUr;1d81eJeGI{8zuTlFuMaoM z&4SNQPyg-#pUQV{cK3L0&VJ`#A7R4e^4A-(!t&d4Xd9}=%XT@ZQ0IX|r}F=ZC@WM8 z*8jwP*j8t2nl}6&js4}hnucRK-%=;?XJysm_t&=hsf;wyQ~N>t#5nS3Mw)wu#X&`T zDnrfSLL9bRjL$B{(b%75$7cpVLS+9jVPNn<*SNbP=3>6}Js7705T+E4Cl~v-Wo)>2 z;XvLD(2)MmNUDm(6PytirL*Fb$2d}%3Xa}ed1&OGD%G9(pOZPkPk}hEY?Hy${7SgX zHF_7%$<~E4b-B)cFodIt(BRqspmJBi0Ni!}p0Zcs<34KhCQoEfJf@j zQ}HSpz?n68PJ0Qf>s%IG=OhcQpUCyN@fM@vI@Yq>x}eIwJY|JVx;R*(dX)Uz*SFL= zpPokn7iXW!Z?69aG%~ANj_FJ0dCOC;u&Df0Y)B&g7k{e>Z&@5bb)U>IrF)JI;B9z8 z`t~x*k*LnZH&Ry387jdmi~^6rar}F=5XfBh8iwPb>pOZH2Z8a0!~pBr0&BW5lJK{4 z6Intm|DXI`zdrttB~#Z|=N@Dfv78VyAjNP|bTSm4P4i@%TyU2?8(I^<0X{@iIn64i zy$cTO^4Jf<7u@Rvw-|sy+(dDE_tgN7^Ae?T%;hM7x!Afc(`#V>PV{G8khjm0fRQHT zMnx`i|E++%UuPJr()K@*fjXRHHq_AZrJr~F3`!GyrRj8c=80woCO|ShD+EV}(-?eK z(H$(!Q_Xhz$)wRZ_FX;DHkT@#!?zKI9U^LIUNJ3vyYQ~MFSx%X(bU#qe{=i86?7|V zqV>B?dOd*b?Q=R+y?DxCrE=oBRr&jmBaZr7$_Uq(*~x{=h_)?l#{T}$`czeI^TQTT zmF@j&FClrp!3oI(%FfYvh+=InVwQT%UEt(V>8^Hjr?b+tF~X<;A>Kara4wN1W&`?- z59sEpoIG-WGZe`xm&v&C$<)10=e~>==}mWPKRE*}9|x@n9DF;iTC5@@eDR_6wRLKw zDb$$vFV|lX1f(W|-R|^aer)n-d{f7?!}=@R(H=l(QRp^t1oA^Ocg~EpL&dX7~uY+qRgyUnP>QC5V`MC#L%NDhd>SO=W zUum6lQ};!)!u@63amG%GdZBQY6>nq%Wf>!uDZm|}4^0klg_?0sA}Y09_>KNUzZ)8? z9+FJBEPGhatgKR--BC`HeR<&FqH oUwx2~|7|>&Q2iN-3-6)O{~;_S%%37*sz?MA0kNx_R2VK%2>F6KrCj%=7#uxavM}Km-Y>Nwd zo(8dbO@Ld9v>pvbhqD&xmR&1z++W0{BbC8O@3u34+E2p&E%2GpPAwV=VnUW@@m#^;oMF8I4L7zWiV)%kt z*5J^tFD@w-zeAFllGL6f9N4;K0YhAsWJ64_7nfQUQFECACTWotM{oQ}>(UU=;<=~o zWD3}b8$NL}IT7FeppRdLUC8St6)Rwr1o0A%(EZN_-KsC$aiBDKrHPnig@2M7PkhID zHFw5Ds6kMhB;KR@XO>skPG$5OiXP(Q(ly^s-?SmI8<9f*Xnhty_O2E;uDMx?n-L$x zc1J$!J3==K%}Cvd3`tdpoM|zl*VrB^8xM8C3*|t56PgpsJ z*?JYWUYz5SJ7`!IC;qADQ|RB{C^40!+bYDOsA| ztXa$b@E>Pnxsk8TC()N9RUVJg_d+-MS4^^V7mmkoed?1FJElK++u8E;LpLX%DP(xN zn|M6k%^xd5&(G5c-Ly4`$J3~C9yVADF9*yEsTt&{$;K|+8^=WWcN8C;;vpXnZ+UYS zEkLha_Unjmfq=+;>9C=Ry%*InddN4tA6+myE?U1n8wz#Ve=_MNkD?#_R`^&~a`ce6 zpFkTieaKpVFrYfkDzM8O5vf{TYgbI<-NyIhZ0|Tjp$ugi?@z0gz(#-ZDHciBAVI&@ zjB-yF%Z$aH#rBX?B#{~ho<%aaQeB!@14ohWu3Uz^F7OSCjiAytL|3C3Oob?h0vq9@ z$OL=&%q{u$B$0}gqLqy+B4(t#Yfja2o6;29E#3fwv?ve>JDieOnwHcOf7217?Kj2Q zSC6ty~mD)Ec{w3)sgFZ{9T9c>k82lIL?i#Rm|=%UK0 zC|S$RUNp?V0>`(A%0$&OKVMh_`zO4I7!vZ}Peoz5ESj-0#!}%fn2Z==6!?&Q3J|9F zV8U|F?C!{A#_AH&A2p~s&>n88lam%}vUxa(32SRug*@Yis{P?2mB>6R;mckr7ql(_ zeG$fCDpulNT2vjCpvc~{I@=iHt&BaVe`AllUuX1phA>{KQ3E+ zuR3pxTr&I2R$rP)WHD#PFpw@~8C`t?VULAJdpJmZq`VluvoIS}HiL#kL5*HvYmgHN zz4D2*5V&M$h`7H!vS75zQ=ycnIgcq)BJ+G^-F_*}_dE!`^rg_bz?BJzl%P`HKH-Cw zOJ@EoQc@qn%bFscz^d8ohF#*4On!VBsZeG)2&Lh*<&UHey}(XW5hj=4E?U;bK2ewl zyW)qyMs+TfoI+dDtOdbJ{En}_r%*?Jv1zOfWkF5%*L`PGy&KKo;)mkNbgC|h@JfP^ z0doG1U(ZHhFJ1mxkzW9~+fmJ9?-Z1m)(iZt`5Bgnb8q&FU6wqk8g{iS@TG7C>}ufu*nO#Mh}@dFkZ^?O+25iGwy>NsDf>BdtsYT zK^wXG_3|?lbRJ10-V@mS&2}m#o04bbS~~B=v{0nMR*4pVlcDUDLRvP2RluLZ!s-Ut8*s*D9GAQ2l0TuXmwwK5VpUAq<<_ zGJ`p)T=W%<6#G4RD}n17`lIS=FTgo6!ysBCB|%{9Kqm|O0th#)Lg-|pEIPkN6PCB4 zE5o_|N5#xbKpK_)whOCsMoEQ<`_5REGbL?r!|0yaOK?bD0h{yHzH2>kzkMCCf_v1r zj!a6{Ab107i@Sx1sFFO6LrQyoOHTDS?62upgsMXU9MYrTY{XxV|-H*)zu-fL} z89}UHFF!GoPS*cphX8W5W0a(j_czYmES)Su|HyM3a5EO@*CY2Svh_nkI7UNGQAGiphwkee5J$FeV?IRkDEDBy7@f z`XH6;OaIt>Iro8$83A1JjqQ`!)O2-cdD5CjCM%)sl5bw9?e;fhBGC>2XE;hoH4pLc z4K$c4@g-Y7EEBdcNeh%G8hp6k-Mpm$AO=k94ZeV|B2z`3OZ&mEqo2LM4bTg!iGFn! z5cPoU;)c^*5mh`&X>h1xFOO%>F1C||&rL>>{y7a4dbhk?zM^Vkrlo%jE>S!I^{DeHc&RVy8Jq7m?__FdhKx+%T=10tQV%Mfw9$;XJ2Mb`T?emPI^)Wz{stpYTLA)Y zjSM0dMM_P)RNW6#b-mPKhX;lf8akhYXlaPP=so_1s!FXP=oVV}M#H6DeB_^8x)9Mw zGeZcgT&XX;zGEY`X~~tR54-62^ylS(v&k2Hm}x~JZzEY*rt<@_sulw054uOJnh_W5 zWom_$@oF|N+YtDA5c`p(Pj~~WQ!xNRgl^UOE2?7ApWJqM{mz>Cg<*{ng!t=~tVS`3 zAId3Bt{}jmw5Q&hx0S+UcSC6LP)@k{lHx>%VVDSyUj{2hpb6~9NUbaAcP9A3^ zJTpp-IZ=wW$x!GtBYQw_9pyWaP3D*-b#(xt^1`G#SdUZw1J+k7ezgKClEHH2?7HdQ z!<|2yPI2BX>kk9xx0`8x45}MUnPvJ_OYURB+=bcXasAgAHi;%Tl@7V!q+1*VK0R|YCZ_)d#{Pt|l{L+moU z3tTA(7-X{?pyRULZ@(@L7eCv=%;Nzqr4CkCFA?|VRx`q?Yrvtv-PqUWato8U^9R>G z5IgnH+jwz#WmEbgI)-PaG-VZ`HcuQUPCY*G8W6Y)B^#ML$HA<4e{4sT4WF?*^GtJn zH?u^p;C`4xYCtM$`K*!?x8QSDTCq_0TBHN$D1qKURmf2%&Oq;VV@zgG{g*NNSU5L6 zQ=Dj{fAyo#kB(Vo07e0^=2sd}iR>LLl=?CFN?WsYF=!E|7P7_#IL#-&7@Ic%WM&fS8UiqINsWP&x=ZW6+T zwKiP>8m1+yQCAt_fQA_*D*4Krw7j*5Hay0~Gk&6pWdygJ)|3A9KG{!zpqQb({|U@+ zP47kMZjJ2=#V)s1(o*L4pk!$>3?)ll0TRaK6c@=|`7+RuDA^w&)V4tF?^oE8%Un;C zCW>3tEC~S*?<;OnPV8veRG~CX1^ep9Ahsg*f^)l=T4v^;t48GZRbRXoRXt z_4XTE!|O2O8en5Q_L;%3LfoM>9!ir0q|b6Z_p+p=vk9q6cX+;v=O}4s6nZh zZfzwt5@niDJ73+PotQ;yp67*Y0k)N)c!*_RjRo4M1ZAkW#B?k6Ohy*j!n>`-QZpm4 znrw6hNtlvTZo!-H&O~do-{HgIQoDt{vttyP1=%V{)lT6@964-c>d4l)6^8T zkI@_P`ea!4Bd>)+g%nctwxm;VMPxAIgz0NK=w6yYO{G@1B`jN?;`ufUi*eT>+lHQ>Pe1g z8x21~yB5yl?{BQ&DfB5RcDcncSSnUq7e665F{d*)t$?aSm7ahLwgZe?VK zgw~IOTz4^^N zFzQx~P07BxlDH&CBQz-ZOXAR0Fl)Yk-Kl^gD1$gm%?ekY<|$rV=~G(V%MF?~Hp3uX zHS1p?kTEQIa+c4QT{T8ZTC|z~O@VntFEl#Pf0j(Ud7@vlp~gA^)GnWXv&+oy$FFA1 z9||#LRg=Lqw%>es2RCWc9nn={U0as}b0+YQZU}2(ZKS^|aK?LTz8+{p@F9N0n}fiD z9CNbJ!)WcVkKHS1JnHX;8mC5oK7F=wq9Xfa7E;2zvTijDEHTpFVH7E_Hj~yX9*6VM zq1y+G(Mh@X7*`q2{K$^Hc5jTMw<8L*+%NU*%|M+fQ)3RyDt zJZHq?#!+HX{IDCruXYTtd{mVD8Wm4yo!DUGr@|ETMAq0?-7@s2(eoeV&m4(E=6@Q? zL#FS?Y2>??XaOkiWdRCEQO!=XrT_uW+BZ}nY@L? z&5#<;)Fj&=n+id&MYXl@f_8w@WzmqDC#K7oPM&>Y>m()Nla4c}0$tp_FZ!*;%BaTK z#~+vFjX$H8e^!ITMkvPRQ5cZlZoj91eHgRg);Lb9AXfElfXYGB8)uxB{vqts$-ptj zv!atwWY3(a83_=%)$vC@h^Wx|QOd*Uo>_m?1B4rU%0mIlfkOE66HFNnlOqd(XJ8$!h8Cyf#6yMISNmATYO%ZS&ImD{N4 zUkQ4c4m&u5?mJG`sFM;}9*2bzBAC|5M>dQvn=fZsSDC=<-C=q1%DEc>-{Xkr8>Mnu z2x&y6=^KEzTv1~h^grIwdG1q+N^Kiv5n3t`9dDaz{=Po~^$Ky<9x#L$4pfoOK>l2d+ygxcJ`Y985BDz7K$CM_D@$NnhnfLd+mlwFc zn&M@kX8-(4SGA@6-} zr#gYlFd&fbSG=57Vn-aWH4lr1O>b@q^4rL!xz(Nz6$|Z<4R{tLwK5Q`8fkb+EVZPh z(GDDDY>oNsjq*4ocSO4m>+$t)*NZ3{Mzm+jx$kue@@`2`z4EjITJN0K-B-Y@oqj!MW85kFwn0m#5l2T{jw4WS$=u6 za$vb9|J9lB$oj|fHuFd*XPUCN1;@4-M*v5O2WyUos6@iP>3DSPyZR(mr>TBjeJhV)}pn z9pPn=zTXN1*d(G9c}UlOfkXp=v|s1;sWV>qD8{5$1$(0+oI_{*tQ9*nHs~)~{pPou8=~m!q#?gQIpbHqAzl)9oBpA(RoslJFh%5LHVB%%c)$v zwHg4yMfZ8qsPEIjcD02P$@M(V&nO-Eb z-q-sJ&2z@th<(M5F?t1wr=1xT$lDn75!Lw;!k=%5n9pf+ef8BSLry%Gmp}UOMmZgw z7u`p;sk z?zf=HpX)r_ibRHT;47j zSO3Ee8pkSHoN2$(6Y1;B@ccnc9Tx&7^3xe%70$kz-UWI3%{812cqK_+uH-*ChWcz zUYrmpHgT0e%mWLp_FqfrME)t8>TH(hPBWUkKd`OMnRn8%@~yJ=4L=cf@UAFQiV?#< z`5aRxxLBI4ILYsfKo?s1t%`x!HK3V{uAr2~6BN(Cxm~~J*bmx?SjKNP7G0E@=Bd-!Z|tZxm7dqSIXtp>U*L2(!uvIhMY#vK6nk@}&_B z%eSs1s-sP-dK5#m4%8)=!wgr7N~BWvVH0r9AC7hVwjLT*ts#%tNc#G8Ghch~k2^)y zO%X2lu!^fuykJ)Ed*OmX;AGkyjt^%_3abzMSrr}+;b!J4rzD{?h{B4Pgm)nv#JXTr z8}H8yMlk-{wt#&Bqq{xkDu0=yPLHZpK2k7kC9f>WKgv<8A=1Q zN%ail>v;cSY~6{H7PnP+Y|x4u_D9(Zie({hB;xMFiO5Y?hn7h0&eQ0c)E?O#5a@Mi zHw(bFA|&%GcUp^Ya4Gm>u5Ay~4jbsTVUX-533)P{VKc<0HURsjn)pOV@M>`nY?GR_ zsF)E}?hCd-pi^@tEy4_g*%9h4qv$m*3&WH)k~;LSp~DT@h}<;}{_;R@QJM;unh|Ms ziLLV@3+szTR<$oN1uv8IE}*d|nfaC5PQ zMOAx81#A<`3iFadb6;S*;tcIBl+$bD`?sifW!^}hQI9aQ{MomzA5aaI;h_Gdu^R_> z(QWnR5k$F!>$usnuZP*rN|rSMWB556BA&AfcGUYjvCw##cUW@A0ows64gO9b%J1Fv zGd%{M7+eq0Ie4uz9{)aAQa94p1@AeXaVgp#b1Aq_Btp)(B6PJ^vf6z~FKlvKN78)> zQ5F-w5)8FT-fgO%fe(fxl@Ze(Mz&7Zfq^NT#Vr`2j%lYKa|^v}$fAMN99H^@a{B)C z;y>$Sr0Pk*ykyfU90N1 zmovls>_QaB@|RNjV8FTQw<4~Jf@10Q2GhSs93E)IMjmy%1S540FQi7@-*L?~%_O{D z*%{0dGgO8rQ!0~^|5gz&@dEFDxM~znJTsshLX@+pCX16Q&o2{`8g#gd9)dsFJt-zH zmfp&`;fr1~LzgXu0P`mF3@?Yuv(pzhZmv!4rYF{4Jj99*=qz(Q3?CD2ZnExAH$nq7 z^o*3c#19v3@*plXI&YitVpr9JlMR}SS<>*eYmNZnoxIbc>j7A?tGNFQmhui&?Q*jP zGWQSLv?z;!q{f)*ffL~X6l6*?x-M)?bsEx=Zfv_}w(Lx=e{`O+SYgmCj#-etW^6tx zI)j#AF90iwGjhuVSzQ*p3UPCW$+(ekR_qW2d8}3J;3NuqVeNis)~1uco}~DoI8%YX z#W!DLlsy4a(s-R)7G^%D>A%{E`}kEAN3B4WF)}MyK}u834*N-P2G+%E_D=zp&&seN zc=d@37^v(Qs+w{OBAHlrV;XXZX;wZ?FZJ+;dVo3203&Vn4Yf#Bn#}6NpcS%o78lPw z+7|;+N~IF*&~7Pd-XY#Y=f=yj`5-bVpzyi`GOE0T9qqBl;_<+KWcE0K|Kxk`oeu~2BN2)Grbi*~x$sTyRD#8qgexnju}jPenll`Xi~-(2_5>v^j?uF4vS zD)D=rX#6_MFV8A?9#`6kJFNPuCX&70K(U_wv%kg>sGl6YaMoOHt2qsRn)G~U^d4ybBN97BDyB>L)Rp~MR zaK5AIpF(brm}?9~GlNR|XliW^cxQ;pZ_ ztZBtn?(lolNnZgUTm8F!AB<0){b`T~GqRj@Hg&sJji6vI^}XT@clUB-^BNXU$w<>f zO@9IiZ;P;^NP|>qMe(${RE4xEU}Ea}Dvuny!gNLzmOgAy2HH~CDA%4$EnPZ&0yFJo z-MM;7uB|q*{PdXjp`w76dkT@O|J}k5IXV)vqzdj=PA?9_zBLp;xGemt^-gbAJWa19*d5@Q^`fJTkP z((85ZJmV?htk*#64PMK0KK3Ud>&la(wl_7BIbgl7bof+&?eAms5X&-AfSZ_@`Ki%E4Z0+><`#|5 z1O@GlKh823nTfe8J=)ko1wDsWF#EW^F&0pDOa6gpqt*pDN93B~;D_K2H#~4JMQCvG zU6cVUsTQ!awC5_*(_xP=X++cb!}2>>q{t}x@+rm?yu5|j`C=cI=RUzNAa-qSbF(FE zTIS(gN>!OsL~ed>_z3npn%J-V02`dh74N(lu(}ob1U6lyt|=6LW;M(=gv=L@HkwQ= zlaC@$sLcb8tzhbdT&2FdVrZB48JX0_;v-DTj9TA(W6A^GvT65BSZ!pj;i1Zovo6a} zL01`iO61&0mNP4pc&#GSytlQa74k_1c1VWR1Q?EFh2d~_$2vo?F`o`3DL>rZgq?)# zepV>EkuCWre|s~v`vUQ>8tq=xkH;MGsQ27a09iQfMG+sT`(3QfC#n`w!)-9eeSth8*>l1!rFfkpZmzY(v?6%9LH$LlJ1x3t!*@)v&36FT5OI_0@ zCHodB!Tw8!hTPH=4lcy8CN0oW1f>DcP7UG#VL8<)uvIQiD1II_KhH>-^Z4U;Z|v=e zzr4EoASHC@gzui67J#S?%1W$yUmv8`-z+?9F*)o%B0yXFw?#6KsuEYrn$Sb~-s+h# zZMJQ+F1roR^#$Emw4e*g6WijHWaNKVvrs@iCHdM+Ogzkq1xnfC^D48neZd2$HL&_i z{-qbHxl@BGh^0i4YSIXPfqVyr3@h3fLlk(FgfhsIqVO;fmx7s|qc3oj^4Iyziiw5I zo(8=-{L-u10}Xa4L9fU17}jj4(w77r!P1(vb}brurdH-(v)fgPGMUj7bBF(BMigN+ClAPezDEax|%_cPh8pT(xNHEOO&GZbx9|?<{ zv5-t`)F2oZ%u$|_<1kiTWxA73lVEAyAFA&1?W!6YANPF^;d(NT1JIA~XAL}b8JW|U zN%R$}HcU*pJH7HcP_}-QpMByG-8D;e+}$V`1PfQsp_NRsHj5!yNVk6n4S4m0c5IK`3S`g5QkJezhZ%mE)TORi5Y%4QeF_Cz?Y@@nH@?d z82lAjOx411(Cq&e40Q1yD$YIahW4#!5snY?*Q(*V4Eu1)Am&mDS4~Q6In&IxMhjIX zjBIaxkoXwg)xs|y;0RaMW6FZQs1fx*nP#CMG^uOV!~f1D;`CL(m@U#|B;3^=3}(WK z`+TR89pRn-Fj)ripQ%iXrZtv|h{TqOffFpW=Q5m5hB@93cY zR)Ghwm(jJVU0hMjZ%50yIs1y64xdsqt!tU7q38(OslxU5zN5P^D%V+ zZvJa$)>BhTURFo~8z*>i1kwfBdwZvb*_W?Rqxw$E+|N#P2Xwssu`ruPlk=f$imN84 zT>i9=8b97qdI>o~$Eui*AYdwboQgx+h4S+*zA>B{4=5Z=iqg@l$m5`B+N`c^kOQlY zlMXNmk6|Yi`0n4tYcVty?!7R3if5~Q@&q%(`;0llVXQiyr3t3lV#7FX8$I8QQ_!{Gw(5jq$DcaBJb z*>Ai$pSLXVNu=T~6D6{Q?UmtnC(4i#ymS0RV74Zs>Bw$g2XWNfvV*dN@I8?x3YH@m z8nj{ed%mZBSE)UBP^j%H5lkZWTlZ$Ka^l{nuPwL(Dy(6zifA#eZXOP_HqAA^PfNYO zyeur(W9WqrVP_=&VEr-b5qNutLe+Etx0|lRyJB)Y=J&UEL~$}V`(fsHrza)m(vs4^ z-Qn(5-rdd7GO9`D@8v!4_XK%ye}1|9hXJ?;6{-=b`X!7x=Y=AbGi32zZJSTr9v!tl zDKIJK^S(;x(?UV_l{xk&Jv{#pFr1J7jALjOnZ}<+Lg}@^uvGiS_31vCEx} z_iK(DOq}*5(~a}yS6yPdZ=sI+XSEyr85E`T2P~=aC&j#IUDfb_02gRP?w16fx(Q}( z9q}+ADKh4H({o;`; z*YB+CtwnbQ0*?X&>80x9ICppO^ohJlt~c0;zN21^rprZVG3_4l6ie^7P-{um#tmV9 z-d0EFcKtg9`TK~6ooEBaPE;1wit~5cQV@)ekXiKB`FaxN9KGWAnYfj#Hiq~o(l-bF z+}?V}%!V%RkJ(CGr0)CF6n%-6{x?M@={34n>j>aAgjf36ZaNp1*i3C6!luaHsx+VBte@Clc>BEfzZ(2$5v{>39%U$+ zkubfibicNK5_Kij6&eIbBmc(Lg1ecB>|0wZ;Sn!~&ZvJJJkMKG4SQ=F#il+$YP0;0 zm(d!m`{w^?Fh`N|%K00>t`Sw&1b0Q8juKW-y;L+=EcLUcDsp|Yr0|d7=Y>a2#^y4X zn@`)H{W}wuGHlY_q_K#e)WEa1=s<5_N%0coD=-G5Q)FVWIFZ`H5Wyzd=2BY0B^_<0A?8ldIT@7LQc1&MmPv}NHr^w6S+ zN@R2JG8I zw`d2x(u{y^y4F8D(h8n(xm1#fZ$4#`?yM3F{;vohZ>3KhdwCOCwI|a8xxLfPkAHc? zMvUlY9jo9Ni37(HT%=H*^*}4zLB(dgSl|1DF%)y^fh{^Gstq7oPmpjO@4b-y&s=vjkWU1i5T{G+;z zC*;2e5PH7ye+bplv-6rvLYfy(3QfcU-=kc%e6w)TG<