{ "cells": [ { "cell_type": "code", "execution_count": 1, "id": "7c1065dd", "metadata": {}, "outputs": [], "source": [ "import numpy as np\n", "import matplotlib.pyplot as plt\n", "import pyerrors as pe" ] }, { "cell_type": "code", "execution_count": 2, "id": "20f67709", "metadata": {}, "outputs": [], "source": [ "plt.style.use('./base_style.mplstyle')\n", "plt.rc('text', usetex=True)" ] }, { "cell_type": "markdown", "id": "e5764fd0", "metadata": {}, "source": [ "We can load data from preprocessed pickle files which contain a list of `pyerror` `Obs`:" ] }, { "cell_type": "code", "execution_count": 3, "id": "c49ff771", "metadata": {}, "outputs": [], "source": [ "correlator_data = pe.load_object('./data/correlator_test.p') " ] }, { "cell_type": "markdown", "id": "ae93c7c2", "metadata": {}, "source": [ "With this list a `Corr` object can be initialised" ] }, { "cell_type": "code", "execution_count": 4, "id": "33a8fdec", "metadata": {}, "outputs": [], "source": [ "my_correlator = pe.correlators.Corr(correlator_data)" ] }, { "cell_type": "code", "execution_count": 5, "id": "5f954607", "metadata": {}, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "x0/a\tCorr(x0/a)\n", "------------------\n", "8\t548(13)\n", "9\t433(11)\n", "10\t343.1(8.6)\n", "11\t273.2(6.6)\n", "12\t217.5(5.6)\n", "13\t172.9(4.9)\n", "14\t137.6(4.6)\n", "\n" ] } ], "source": [ "my_correlator.print([8, 14])" ] }, { "cell_type": "markdown", "id": "b00d670b", "metadata": {}, "source": [ "The `show` method can display the correlator" ] }, { "cell_type": "code", "execution_count": 6, "id": "b71529d0", "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAlYAAAGLCAYAAAAF7dxzAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjMuNCwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8QVMy6AAAACXBIWXMAAA9hAAAPYQGoP6dpAAAeLElEQVR4nO3db4hdZ34f8O9jO1gk1cy1XBYr9S7dEaHUWSiRvYGF7ptG2i4U95XklLJvV6rpi9K+sGpIMWZfmPGLpC2UxcpbU0pkSlqn1CBtoaQQSLxDC8aEtpotG+/KTXelmXGSlXdjP30x58pHo3tHc2ee+2dmPh+4aM75PXf0HB3PnK+f85znllprAAA4uEfm3QEAgKNCsAIAaESwAgBoRLACAGhEsAIAaESwAgBoRLACAGjksXl3YJZKKSXJLyb5aN59AQAOlZNJflgfsgDosQpW2Q5VH8y7EwDAofR0kh/s1uC4BauPkuRP/uRPsrS0NO++AACHwNbWVj7/+c8ne7jjddyCVZJkaWlJsAIAmjN5HQCgEcEKAKARwQoAoBHBCgCgEcEKAKARwQoAoBHBCgCgEcEKAKARwQoAoBHBCgCgEcEKAKARwQoAoBHBCgCgEcEKAKARwQoAoBHBCgCgEcEKAKARwQoAoBHBCgCgEcEKAKCRx+bdgXn45NOaP7j54/zpR3fzuZMn8qtfPJVHHynz7hYAcMgdy2D1td/6r/l/Hz96b/v08om88vwz+fqXTs+xVwDAYXcsbwX+362P79v+cPNuXnxzLe+8d2tOPQIAjoJjGax2qt3rN373vdza+Mm8uwMAHFKCVc+P/uyn+dbvvT/vbgAAh5RgtcNPfvbJvLsAABxSgtUOLzz3+Xl3AQA4pI7lU4GjlCRPLZ/I1375qXl3BQA4pI7liNXOFauG2688/4z1rACAfTuWwepzS4/ft/3U8ol8+xtnrWMFABxIqbXOuw8zU0pZSrJ5+85G/vjHf2nldQDgoba2trK8vJwky7XWrd3aHss5Vo8+UvKVM0/OuxsAwBGzr2BVSnkpyUa3ebvW+lavdqnXdFBrfX3He6daBwCYl4nnWJVSrid5q9Z6Ncm7Sa71apeyHXaudvX1UsrqrOoAAPM00RyrLticqbVe6e07W2td676+meR8rXW9V79Ta31iFvU99H8pyebm5maWlpb2fNwAwPE1yRyrSUesVpNc7+/ohapBkpV+6OkMSilnp12f8DgAAJrb8xyrLtgMsh1khvOc+qNXK2PeurFLrWV9bWehlPJ4kv7aCid3+T4AAAcyyYjVMNyc6s1xul5KGc6xOjXmfbe72rTro7ycZLP3+mBMOwCAA5skWA3Dy7vDHbXWG0kulFJ2G1Gap9eSLPdeT8+3OwDAUTbJcgvrO/4c2khydsT+oVPZHlW6PeX6A2qtHyf5eLhdikVAAYDp2fOIVW/S+M7RqUH353pyby7Wzvr6DOoAAHM16VOBaxk9n2mt1rqR7YDzQL3WOvX6JAcBADANkwarK0kuDje6pwPf6o1mrSa5sKN+pff+adcBAOZm4g9hHi4SOtzuLxba1YcfdzNI8uSs6w/puwVCAYCJTLJA6MTB6jATrACASU1z5XUAAMYQrAAAGhGsAAAaEawAABoRrAAAGhGsAAAaEawAABoRrAAAGhGsAAAaEawAABoRrAAAGhGsAAAaEawAABoRrAAAGhGsAAAaEawAABoRrAAAGhGsAAAaEawAABoRrAAAGhGsAAAaEawAABoRrAAAGhGsAAAaEawAABoRrAAAGhGsAAAaEawAABoRrAAAGtl3sCqlDBr2AwDg0HtsksallHNJrve215Ocr7Wu9/Zd6r1lUGt9fcf3mGodAGBeJh2xGiR5tnudqbWeGRGqBrXWq7XWq0nWSymrs6oDAMxTqbXuvXEpF5LcqLVujKnfzIMjWHdqrU/Mor6H/i8l2dzc3MzS0tKejxsAOL62trayvLycJMu11q3d2jabvN7NuVrph57OoJRydtr1NkcBALB/E82x6rxQSrndff3lWuuV7uuVMe03dqm1rK/tLJRSHk/yeG/XyV2+DwDAgUw6YrWe5N1a61u11reS3CylvNHVTo15z+2uNu36KC8n2ey9PhjTDgDgwCYKVrXWtVprf2ToRpJLC7z0wmtJlnuvp+fbHQDgKNvPrcB7aq3rpZRk+1bc7THNTnW1addH9e/jJB8Pt7u+AgBMxZ5HrEopg1LKnVLKSn9fr8n6iH3J9hIN6zOoAwDM1aRzrN7d8VTeSnLvFuFGtgPOA/OdZlGf8DgAAJrbc7Dqgs31HbtfTnKlt72a5MJwo1vQc5Z1AIC5mWiB0CQppbzUfXkmyXe7FdB31jeyfYvuyd5yDDOpP6TvFggFACYyyQKhEwerw0ywAgAmNZeV1wEAjjvBCgCgEcEKAKARwQoAoBHBCgCgEcEKAKARwQoAoBHBCgCgEcEKAKARwQoAoBHBCgCgEcEKAKARwQoAoBHBCgCgEcEKAKARwQoAoBHBCgCgEcEKAKARwQoAoBHBCgCgEcEKAKARwQoAoBHBCgCgEcEKAKARwQoAoBHBCgCgEcEKAKARwQoAoBHBCgCgkccO8uZSyvVa6/kd+y71Nge11tdnWQcAmJd9j1iVUi4kObdj36Vsh52rtdarSdZLKauzqgMAzFOptU7+plIGSS4lWa21lt7+m0nO11rXe/vu1FqfmEV9D/1eSrK5ubmZpaWliY8bADh+tra2sry8nCTLtdat3drud8TqhSRX+zu6sLXSDz2dQSnl7LTr+zwOAIBmJp5j1YWYd0eUVsa8ZWOXWsv62i5tAACmbj8jVs/VWkeFmFNj2t/uatOuP6CU8ngpZWn4SnJyzPcAADiwiYJVKeVCN2n8sHg5yWbv9cF8uwMAHGV7DlbdHKeNXZrcHrP/VFebdn2U15Is915Pj2kHAHBgk8yxeiHJmd5E8TNJUkp5Kcl6khvd9qDWutF736Crr0+5/oBa68dJPh5ul1JGNQMAaGLPwWrnLcBSykqSS/0FOksp69keQdrY8d61WdQBAObpIB9pMxixbzXJheFGt6DnlRnWAQDmZr8LhF5KcjHbK6+/leSNWuvwVuBL2R5RGiR5stZ6Zcd7p1p/SL8tEAoATGSSBUL3FawOK8EKAJjULFZeBwBgB8EKAKARwQoAoBHBCgCgEcEKAKARwQoAoBHBCgCgEcEKAKARwQoAoBHBCgCgEcEKAKARwQoAoBHBCgCgEcEKAKARwQoAoBHBCgCgEcEKAKARwQoAoBHBCgCgEcEKAKARwQoAoBHBCgCgEcEKAKARwQoAoBHBCgCgEcEKAKARwQoAoBHBCgCgEcEKAKARwQoAoBHBCgCgkccmaVxKGSR5ods8k2SQ5EqtdaPX5lLvLYNa6+s7vsdU6wAA8zLpiNVqkndrrVdrrVe6fdeGxS70DLr61STrpZTVWdUBAOap1Fr33riU60muD0eJSikvJVmttZZu+2aS87XW9d577tRan5hFfQ/9X0qyubm5maWlpT0fNwBwfG1tbWV5eTlJlmutW7u1nehWYK31/I5dZ5LcSO7dJlzph57OoJRyNsn6NOu11rVJjgUAoLWJglVfKWUlybkkw7C1Mqbpxi61lnXBCgCYq30Fq26u0+Ukl3sjSKfGNL/d1TamXB/Vz8eTPN7bdXLM9wAAOLB9LbfQTR5/NsmVbp7Vono5yWbv9cF8uwMAHGUHXcdqNclqd1vw9pg2p7ratOujvJZkufd6ekw7AIAD23OwKqUMSinXuknqQ8PbgOeGX++oJ9trXa3PoP6AWuvHtdat4SvJR6PaAQC0MMmI1XCyen8+06D7c71bJHQ9I+Y71VrXpl2f4DgAAKZiz8GqCy9Xdyx38OtJ1mqtN7rt1SQXhsVukvuVXvtp1wEA5mbSBUIH2Z4QPjTIgx9p81K2n+AbJHmyt0L7TOoP6b8FQgGAiUyyQOhEweqwE6wAgElNEqwO+lQgAAAdwQoAoBHBCgCgEcEKAKARwQoAoBHBCgCgEcEKAKARwQoAoBHBCgCgEcEKAKARwQoAoBHBCgCgEcEKAKARwQoAoBHBCgCgkcfm3YFF8smnNX/4vdv504/u5nMnT+RXv3gqjz5S5t0tAGBKWl/7BavOO+/dyqtvv59bm3fv7Tu9fCKvPP9Mvv6l03PsGQAwDdO49rsVmO1/2BffXLvvHzZJPty8mxffXMs7792aU88AgGmY1rX/2AerTz6tefXt91NH1Ib7Xn37/Xzy6agWAMBhM81r/7EPVtfe/f4DabWvJrm1eTfX3v3+7DoFAEzNNK/9xz5Y/fbvf69pOwBgsU3z2n/sg9U3v/rFpu0AgMU2zWv/sQ9WF5/7Qk4vn8i4BytLtp8QuPjcF2bZLQBgSqZ57T/2werRR0peef6ZJHngH3i4/crzz1jPCgCOiGle+499sEqSr3/pdL79jbN5avnEffufWj6Rb3/jrHWsAOCImda1v9R6fJYRKKUsJdnc3NzM0tLSA3UrrwPA8bKXa//W1laWl5eTZLnWurXb97Pyes+jj5R85cyT8+4GADAjra/9bgUCADQiWAEANCJYAQA0MvEcq1LKS92XX06yXmu9sqN+qbc5qLW+Pss6AMC8TDRiVUpZrbW+3r0uJlkppVzr1S9lO+xcrbVeTbJeSlmdVR0AYJ72vNxCKWWQ5DtJfq3WutHtO5vku0nO1FrXSyk3k5yvta733nen1vpE9/VU63s4hl2XWwAA2GmS5RYmnWO10r2GhgFnpQteK/3Q0xmUUs5Ouz7hcQAANLfnOVbdKNXOkaFhyFrP/YGrb2OXWsv62i5tAACm7qBPBV5OcqMbRTo1ps3trjbt+gNKKY+XUpaGryQnx3wPAIAD23ew6m6/nUtysV13mns5yWbv9cF8uwMAHGUHGbFaTfLscCJ7tkeORjnV1aZdH+W1JMu919Nj2gEAHNi+glUp5Y0kl3uhKukmsneTzPsGXW3a9QfUWj+utW4NX0k+GntQAAAHNHGw6taSWh0+nVdKWSmlnO1C1si5VrXWtWnXJz0OAIDWJl0g9EK2R4hWSinnuu0r+WzEaDXJhV77S109M6oDAMzNpAuE3hlVq7WWXruXsr0EwiDJkyM+8maq9YccgwVCAYCJTLJA6J6D1VEgWAEAk5rmyusAAIwhWAEANCJYAQA0IlgBADQiWAEANCJYAQA0IlgBADQiWAEANCJYAQA0IlgBADQiWAEANCJYAQA0IlgBADQiWAEANCJYAQA0IlgBADQiWAEANCJYAQA0IlgBADQiWAEANCJYAQA0IlgBADQiWAEANCJYAQA0IlgBADQiWAEANCJYAQA0IlgBADQiWAEANCJYAQA08tikbyilnEtyudZ6cUTtUm9zUGt9fZZ1AIB52vOIVSnlbCllNcnFJCsj6peyHXau1lqvJlnv2s+kDgAwb6XWOtkbSrmQ5OVa67M79t9Mcr7Wut7bd6fW+sQs6nvs+1KSzc3NzSwtLU103ADA8bS1tZXl5eUkWa61bu3Wtskcq1LKIMlKP/R0Bt1I11TrLY4BAOCgJp5jNcYDtwY7G7vUWtbXdmkDADATrYLVqTH7b3e1jSnXRyqlPJ7k8d6uk+PaAgAc1FFfbuHlJJu91wfz7Q4AcJS1Cla3x+w/1dWmXR/ntSTLvdfTu7QFADiQVsFqPbk3ib1v0NWmXR+p1vpxrXVr+Ery0cMOBABgv5oEq1rrRrYDzgPznWqta9OuH6z3AABt7CdYjZssvprkwnCjW9DzygzrAABztecFQrv1on492+FmJcnVJN/tVkEftnkp20/wDZI8WWu9suN7TLW+h2OwQCgAMJFJFgideOX1w0ywAgAmNfOV1wEAEKwAAJoRrAAAGhGsAAAaafVZgUfWJ5/W/OH3budPP7qbz508kV/94qk8+kiZd7cAgAnN4pouWO3infdu5dW338+tzbv39p1ePpFXnn8mX//S6Tn2DACYxKyu6W4FjvHOe7fy4ptr952AJPlw825efHMt77x3a049AwAmMctrumA1wq2Nn+Q3fve9jFrhq3av3/jd93Jr4ycz7hkAMIlPPq159e33x17Tk+TVt9/PJ5+2WddTsBrhW7/3fn70Zz/dtc2P/uyn+dbvvT+jHgEA+3Ht3e8/MFLVV5Pc2ryba+9+v8nfJ1iN8JOffdK0HQAwH7/9+99r2u5hBKsRXnju803bAQDz8c2vfrFpu4cRrEb42i8/ldPLJzLuAcyS7ScJvvbLT82yWwDAhC4+94U9XdMvPveFJn+fYDXCo4+UvPL8M0nywIkYbr/y/DPWswKABTfra7pgNcbXv3Q63/7G2Ty1fOK+/U8tn8i3v3HWOlYAcEjM8ppeam3zeOFhUEpZSrK5ubmZpaWlPb3HyusAcDTs95q+tbWV5eXlJFmutW7t1tbK6w/x6CMlXznz5Ly7AQAc0Cyu6W4FAgA0IlgBADQiWAEANGKO1T6Y0A4Ai2ue12nBakLvvHcrr779/n2fO3R6+UReef4ZSzAAwJzN+zrtVuAE3nnvVl58c+2BD3P8cPNuXnxzLe+8d2tOPQMAFuE6LVjt0Sef1rz69vsZterXcN+rb7+fTz49PuuCAcCiWJTrtGC1R9fe/f4DCbivJrm1eTfX3v3+7DoFACRZnOu0YLVHv/3739tTu3/9nf+d//Dff5A/uPljo1cAMGWffFrzBzd/nN+8/r/21H6v1/P9Mnl9j7751S/mn//79x7a7oebd/NP/t1/T2JSOwBM06iJ6g/zza9+cYo98lmBe/bJpzV/e/W/5MPNuyPv3478+7I99PhPz/1S/vpf/QVLMwDAPu1cQuHOn/80//jfrk10TX5q+UT+25W/M/F12GcFTsGjj5S88vwzefHNtXuB6WGGbX7rxmfDk8NRrPPPPGUtLADYg1EjU4+UvV2Lk+1QlSSvPP/M1K+1RqwmtJ9hx/v6kO3/EAY//3PZ+Iuf3dt/evlE/sXf+5t54hceF7YAODZGLeaZ5N6+//Ojv8i/vPE/9xyiRjno1JxJRqwEq33o/0ew+p//OD/cZ8h6GKNbABwVowLU9fc/fGCwYvDzP5ck9w0+7NfnTj6ef/UPfuXA1023Aqfs0UdKvnLmySTJT376l3ua1L4fH27ezT96c22i0a2HJf/d9glrAIz7OJi9Xl9Gtb3z5z/Nt/7TgwFqVHhqEaiG/tn5X7p3vZ6VQzliVUq51Nsc1Fpf3+P7moxY9e1nUvs0nF4+kb//t07nP/6PWw9N/qP27TY6luztB2fWbRe1X0fhGBa1X45hMfrlGBajX9M4hlEBaJLry7i2s3aQieqjHOlbgV2ouhemSikXkny51nplD+9tHqySz5bQT/Y+kW7RjJv7NckPzizbLmq/jsIxLGq/HMNi9MsxLEa/pnUMR8EwRn37G2ebLXd01IPVzSTna63rvX13aq1P7OG9UwlWycEntQMAk3ukJP31uKexhuSRDVallEGSO7XWsmN/TfJsrXXtIe+fWrBK7r8vPXyKITm8o1gAsKiGQeDf/MNfmfoT9Ud58vrKmP0bXe2+YFVKeTzJ471dJ6fTrW39Se1J8jee+isjn3bY+Iuf7XktLADgQU8t6KebHLZgdWrM/ttjai8neWV63dnd1790euRk8FGPlwIAox2mTzI5bMFqUq8l+c3e9skkH8yyAztHsZLRgWu3R1GNbgFw1Ezy0NSijk6NctiC1e0x+0+NqtVaP07y8XC7lMVJtqMC19/9ktEtAI6HYViaZJmfw+BQTl5P8kStdaO3fyEmr09L63VGxu0zOgbATi2WgTjsH9t2ZJ8KTMYut1B3Pik45r2HMliNctCVcUft2+tHCxy2dVkWta1+OQbH4BgW7RhafLLHuLaHJUSNctSD1c4FQu/bfsh7j0ywmpZ5rw5sNWT9cgyL3y/HsBj9mtYxHOYANC1HOlglSSnlpWwvsTBI8uReVl3v3idYAQATOcrrWCVJ9vrZgAAAs/TIvDsAAHBUCFYAAI0IVgAAjQhWAACNCFYAAI0IVgAAjQhWAACNCFYAAI0IVgAAjQhWAACNCFYAAI0IVgAAjQhWAACNCFYAAI0IVgAAjQhWAACNCFYAAI0IVgAAjQhWAACNPDbvDszD1tbWvLsAABwSk+SGUmudYlcWSynlryX5YN79AAAOpadrrT/YrcFxC1YlyS8m+SjJyWyHrKe7bRafc3b4OGeHj3N2+Dhns3EyyQ/rQ4LTsboV2P1j/CBJtjNWkuSjWqt7g4eAc3b4OGeHj3N2+DhnM7Onf1uT1wEAGhGsAAAaOc7B6uMkr3Z/cjg4Z4ePc3b4OGeHj3O2QI7V5HUAgGk6ziNWAABNCVYAAI0cq+UWgNkopVyvtZ7fse9Sb3NQa319xt0CmLpjOcfKL/jFV0p5qfvyy0nWa61XdtSdwwVVSrmQ5FqttfT2XUrvPHVtvrzzvDJ73c/aRrd5u9b6Vq/m52yB9M7HIMmTSV6rtW6MqCfO19wcu2DlF/ziK6Ws9s9HKeVaktRaL3bbzuGCKqUMklxKsrojWN1Mcr7Wut7bd6fW+sTse8lQKeV6ksu11vVSytkk3x2eNz9ni6ULwFeHQar7WVuttV7utp2vBXEcg5Vf8Aus+2XxnSS/1vsFcjbJd5Oc6S4AzuGC6n65/06SO70L9KC/3Wtbkzxba12beUcZnqszO/4n5uzwfPg5Wyxjbq/f2+d8LY5jNXm9+wW/0v8PrzPoLt4shpXuNTQ8XyvO4eLq/v3fHVFaGbEv2b79NK7G9K0mud7f0QtVg/g5WzSnelMk7uN8LZZjFaziF/zCq7Vu1Fqf2DGKMTw363EOF9lzY0afTo1pf3uXGlPUXYgH2b7wXupeq70mfs4Wz5Ukq6WU66WUQXe+Lnc152uBHLdg5Rf84XQ5yY3u/8acwwVUSrlQa706736wZ8OL7ala69Xu3F0fzmeMn7OFU2u9keR8knNJ7iT5o94IlfO1QI5bsOKQ6YaxzyW5OO++MFo3+rGxS5PbY/af2qXGdA0vtvdu3XYX7gulFCMcC6g7L2eTPJHkapJrO54CZEEct3Ws/II/fFazPcF5o9t2DhfPC0nO9OZynEnuPcW0nuRGtz3oPxqe7VtRO+eEMBvrO/4c2sj2xXvcefFzNj+rwyejk1zuRhevl1JuxO/FhXLcgtV64hf8YVFKeSPbj4Jv9HY7hwtm5y3A7v+sL/XX0CmlDG/jbux4rycC56B7ujbZviXYPweD7k8/Zwuk+5+W+/7da603SimvZ3tE/3e6ds7XAjhWtwK7/+BGztPxC36xdEPcq8M5BKWUle5R8I04h4tuMGLfapILw43u/FpfZ77WMnr+zZqfs0PjZrYXUN6I87UwjlWw6vgFv+C6he0G2V5e4Vy3fSWf/Z+Xc7ighoG4+/paKeVc8tmoVvf02UvZXj/JqtDzdSW9uYvduXurNyHaz9mC6MLR2W4+Y9+z3dy4xPlaGMdugdDkvo9wGCR50sq0i2O4mOSo2o6VvJ1DOKDhIqHD7REfHeXnbEF0vxtf7jZ/nNEfaeN8LYBjGawAAKbhON4KBACYCsEKAKARwQoAoBHBCgCgEcEKAKARwQoAoBHBCuAhupX/fTgx8FCCFcDDWWgR2BPBCuDhVnof9QIwlmAFsItSytlsf2AxwEMJVgC7u5zkjXl3AjgcBCuA3bkNCOyZYAUwhtuAwKQem3cHAFroQtBK90qSt5JcGNZrra/v49teTrI65u+7kOTLSX6cZD3Jryd5rdYqiMExJlgBh163xtRKrfWtbvtOkjO11sullDeSPJdkP8Fq5G3AUsqlJBdrred72xeSfHO/xwAcDYIVcBScq7Ve7W0Pklzvvt7XGlTjbgN2Ie6NJE/0dq8n2ai1buzn7wKODsEKOAp+Z/hFb4X0G0myM+x09QvZDkMrSd4aMzl93G3AN7r39L/v2eHfBxxvghVw6O0IOeeSrO0yenSt1vpskpRSBkm+k+TZEe3GPQ14Ltuhq+98PhshA44xTwUCR835JO+OKnS39wbD7S58DXZ+DmAp5VzG3wbMiO9/LkasgAhWwBGwIxidS/LdXu1Cr/Zcko0db9/I9q28vovZfVHQeyNZXQhLrXWtlHLWhzXD8SZYAYdaF5xullIG3de3u9fwVt+pXvPBsLbDqR3bI28DdvuGc7OG3/9yPgtr5ywmCsebOVbAYbeW5GqSF7Idms4nuVJKOZUkO54W3MiDISrpha1xtwF7Lia5XEq5me0nAS+WUq6VUl56yPuAY6DUWufdB4CZ6OZYXau1nuntu5nk/HCkqVv3atXIE7AfbgUCx8bOVdG7W3kbO0KUzwYE9s2tQOC4uVhKWU3yR9n+SJqLw0J3G9CyCcC+uRUI0HEbEDgotwIBPnNKqAIOwogVAEAjRqwAABoRrAAAGhGsAAAaEawAABoRrAAAGhGsAAAaEawAABoRrAAAGvn/jZ24vHe/UlsAAAAASUVORK5CYII=\n", "text/plain": [ "
" ] }, "metadata": { "needs_background": "light" }, "output_type": "display_data" } ], "source": [ "my_correlator.show()" ] }, { "cell_type": "markdown", "id": "c659557e", "metadata": {}, "source": [ "## Manipulating correlators" ] }, { "cell_type": "markdown", "id": "416cf39a", "metadata": {}, "source": [ "`Corr` objects can be shifted" ] }, { "cell_type": "code", "execution_count": 7, "id": "e8d65dd5", "metadata": {}, "outputs": [], "source": [ "shifted_correlator = my_correlator.roll(20)\n", "shifted_correlator.tag = r'Correlator shifted by $x_0/a=20$'" ] }, { "cell_type": "markdown", "id": "634dd613", "metadata": {}, "source": [ "Or symmetrised" ] }, { "cell_type": "code", "execution_count": 8, "id": "127a661d", "metadata": {}, "outputs": [], "source": [ "symmetrised_correlator = my_correlator.symmetric()\n", "symmetrised_correlator.tag = 'Symmetrised correlator'" ] }, { "cell_type": "markdown", "id": "3d733872", "metadata": {}, "source": [ "And we can compare different `Corr` objects by passing `comp` to the `show` method" ] }, { "cell_type": "code", "execution_count": 9, "id": "8e264aed", "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAk0AAAGNCAYAAAAM+kVxAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjMuNCwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8QVMy6AAAACXBIWXMAAA9hAAAPYQGoP6dpAABjcUlEQVR4nO3dfXwT150v/s8xCQ8NYGEowQFCIjeEGNJQAXnYhnab2gS2eWprE0hpc/uAXW667W6StS9p8sttk5YrErJ326apzbbddNOEYLfNA78NYIVuIV1SwIIm4CQkVh6AmLhgy5gskAaf+8fRiJE8I83IkmZG+rxfL71kjUbSkWRpvjrne75HSClBRERERKmVON0AIiIiIi9g0ERERERkAYMmIiIiIgsYNBERERFZwKCJiIiIyAIGTUREREQWuCpoEkI0CSH8QgifEKLO6fYQERERac5yugFJ/ADaAYQALHe4LURERERxOQmahBBVAOqllLUG1+l7kHxSytW6y01SyupctImIiIhoKLIaNAkhAgBuBuCD6jVKvr4OukBJCFEjhAhKKRtju/iFEDUAygD0SClbbTy2AHAegP6hPQsiIiIqMmMAvCvTLJMicrGMSizwWSmlnJO0vRNAtZQyotvWK6UcZ3AfnQDmSCmjFh9zMoCDQ2o4ERERFaspUspDqXbIW06TEMIHwK8PmGJ8sR6qMqiASut1ikL1VoUtPkQ/ABw4cABjx44deoOJiIio4B07dgxTp04FLIxU5TMRfNBwXUw0dl0EQJtuuy+2zZaxY8cyaCIiIqKsy2fJgTKT7T0AyqSUYaicpjohRANUInk0b60jIiIiSsFVJQeklM1W9xVCjAAwQrdpTPZbRERERKTks6epx2R7WYrrUlkJoE93YhI4ERER5Uw+g6YIEE8I1/Mhg9wlAKsAlOpOU4bQNiIiIqKU8hY0xfKTIjDIbYrlM9m9v1NSymMAvgzgRQDPD7WNRERERGZyFTSZJX0HAdRoF2LFLhtN9rVESvmwlLISwOVDuR8iIiKiVHJVEbwGaiZcE4B2LcFbStkshGjQKoMDGK+ry0RERETkWjmpCJ5PQojbANwG1Wt2cV9vD8b27gWOvweMPheY9jdAyTCHW0mpnB6Q2PFmD7r7T2LimJG4/MIyDCsRTjeLiIiKwLFjx1BaWgoApbG0H1OeD5o0QoixAPr6Vl2CsSd1VdB95wMLfgBU3uBY28jcxr1d+N6zHejqOxnfVl46EvdeX4mFs8odbBkRERUDO0FTPmfP5ceES4Cvh4CVh9T5xJnA+q8AHc843TJKsnFvF1Y8Fk4ImADgcN9JrHgsjI17uxxqGRER0WCe72kaPDzXi7E+35kdBgaAdbcA3R3At3dzqM4lTg9IXB3cMihg0ggAk0pH4oXGazhUR6a6j51Ed/+ptPtNHDMCE8eOzEOLiMhr7PQ0uaoieCaklA8DeFgbnkNJUudZSQkw/3bg59XA7seAObc60UxK0rLrHdOACQAkgK6+k2jZ9Q6WXD4tfw0jT/n1n97Bvzz/etr9vvPZi/CP1dPz0CIiKmSeD5osmXiJOt/+EwZNLrF225uW92PQRGa+dMX5qK48N375je7j+Icn9+D/3jwbH5s4Or594pgRRjcnIrLF8zlNQojbhBAdAHaY7tT9ijq/6lv5aRSltXz+hVndj4rTxLEjMWtyafykBUofmzg6YTuH5oBoNIpIJJPFF5zjtjY3NjaitrYW4XDqesyhUAj19fVobGRFnULj+Z6mQcNzAwOJOwwMANseAnzTgE8sc6KJZKB27vn4l+ffwOG+kzDKqtNymmrnnp/vppFHnR6QeOlgFADw0sEoLikf62g+nP6AOX78ePj9ftTU1KCxsRHBYDCvbQmHw1i+fDmi0Sg6Ozvz+tiZCoVCaGxsTNvmxsZGRCIRtLS0DLq9tq26uho1NTVGN7dl5cqVuPDCC1FfX59yv6qqKoTDYTz55JNDfsx8Wb16Nerq6uDT5wQ71A4A2LlzJ/x+/6DPSnNzc/zvaDSKhoaGvLbP8z1Ng7R+HTiwAzjVr87X3QLs3wgsuJ9J4C4yrETg3usrAagASU+7fO/1lUwCJ0s27u3C1cEtuOt3ewEAd/1uL64ObnFkBmY4HMacOXNQXV2NYDCIYDCIhoYGBAIB1NbWorW1Ne9tCgQCWLt2bUa31R+k8qmqqspSm+fNm4d58+YN2l5bW4umpibU1tbGA4GhPhefzwe/329p30AgMKTHyredO3c6HjA1NjaioaEBDQ0NaGlpQSQSQW1tbfz65uZmRKNR1NXVoa6uDn6/P++9eYUXNB15RSV9r5qizrs7gMW/Yp0mF1o4qxyPLAtgUmni0Mmk0pF4ZFmAdZrIEreVrqitrUUwGERVVVXCdr/fn7aHIpcyPSC2tbVltyFZVlNTM6i3IRKJoKxMreZVVVUVfy+y8Vy0+y0k4XAY1dXVjrYhGo0iFAohGo3Gt61cuRKtra3xIdpgMJjQY1hTU5P3oN7zw3NJJQeAb/4RYEVwz1g4qxzVlZNYEZwycnpA4nvPdhgO8UqoXsvvPduB6spJefmf0n71JgdMmqqqKss9FU6LRqNYtWqVq3KKMlVIzyUXmpqa8j5kbCQSiSASicR76bTPihYERyKRQZ+faDSKcDict549zwdNg0sODAMunO90s8iGYSUCV1WMd7oZ5EE73uyxVLpix5s9efkfa21tNQ2YNMnDCatXr044OGi9Jq2trWhsbEQgEEB9fX28l2TevHmG27WDXmNjI+bNm4dIJAKfz4e6ujrTtkQikXhS886dO1FdXR1vfygUih/EtDavXLky3mNlt92pDsrNzc3xHhyjdofD4XhbOjs70dTUFN++fPlyAEB7e3t8W1NTE3p6etDY2IiKior4AdfsuaR6zbQcG+259vT0mD4PI9pwbE9PT0IOTnNzM4LBYDwfq6amBpFIBHPmzIHf70dLS0teA2ynh+Z8Ph96e3sTtmlBrt/vNw14fT5fQqCVc1LKgjgBGAtA9vX1SSIqDk/tPiinNW5Ie3pq98G8tAeADAaDlvevq6uTbW1t8cudnZ2yqqoqfrmlpUUGAgHZ3t4u29vbZUNDQ8rtVVVVCfdXU1MTv9zZ2Sn9fn/C41dVVSW01+/3y97e3vjltrY2GQgEstZuIy0tLbKpqSnhvrQ2tbe3S5/Pl/BYfr9ftre3J7Qx+XkZPVez55LqNWtoaEh4fXp7ewe1x0xbW5sEkPB6NjU1ybq6upRtt/r/09LSIltaWmRdXZ1sb2+P37eVtiXfj9FtOjs7ZUNDQ/xxqqqqEp5LPtTV1cX/r7TXM5nf70/4/8lEX1+fhPqNNVamiTU839NERMVr4hhrpQSs7pdPkUgE69evj/eaAOoXdU9PD0KhEKqqqgb9itbOjbaHw2Hs2rUroaeruroaTU1Npr1fTU1NCTk6fr8foVAo5UyzobTbTEtLCxYvXhxPtNa3NxqNJlzWeh2Sh3Aykeo1mzt3LlavXq39KAdgLxFca5u+B6eurg5CCDQ2NiY8T+11s/p8Wltb4+9RW1sbVq1ahZaWFjQ2NtrOuXryyScHzTzUErCff/55+Hy++OuUqjdKn7CdSnV1dcreT004HEYoFIr3ILoFgyYi8qzLLyxDeenItKUrLr8wP8m7fr8/7ZR+LS8jFAoZHiD9fj/a2trSHkSTt+/atQtlZWUJs/Oi0ajhzDL9fWhBUFlZGXp6etIOPw213clqamrQ1NSEcePGIRAI4Oabb05I7E4+UPt8PttDZGZSvWahUGjIQ1ZGt9eCEO31qa+vjwe2+mAoFX0QumvXrvjQp34IVD8smByIaqLRqOH7VF9fj/r6+nj7I5FI2mHn5MBrqBobG9He3h5vg1kw2NPTk9fkfM8HTYMSwYmoaGilK1Y8FoYAEgInJ0pX1NTUpC0pEAqFLP3SHkobrFq9ejV27twZP+ClqyuUq0TqtrY2RCIRhEIhBINBHD16NB4A5OqAqH8uRq9ZrkpDJD+furq6eD0qq8GgPtAJh8ODAppoNIq2trZ4b6A+V02vubkZN99886Db6mtcAer9yefsOi2Q1Aed2nOORqMJ280Cv1zxfKAhpXxYSlkJ4HKn20JE+eem0hXagT4UChleH41G4wfNuXPnGgYh4XA4Ze+QGbP7Mwt0otEoGhsbEw6O2nRvs4rX4XA46+3Wau/4/X7U1dWhvb3d9PXLplTPRRv+009/z4TR7ZOTln0+H6qqqtDY2Ii5c+fauv9QKJRwX9pzWb9+fUJg4fP5DF/TnTt3Dho63bVr16BhRW34MFXQXFtba+mUrkRAc3NzfPhSe07hcDg+NGoUWOazJpbngyYiooWzyvFC4zX44ednAQB++PlZeKHxGkdqfWm5JckHqWg0iubm5nivRiAQiA/JaLSDQybVqwOBwKC6Ndp0bO1vPe3go9+ubTOataQNK2a73drroqfvOUg+SCY/D6PAxGib2XMxe820IE5/nXYAtxpMaTPmNM3NzfGijHra/4uVg384HI7PANRyr/TXAUBnZyfGjz8zW7SsrGxQm82C3OS2abMo/X5/yuVjWlpaLJ1S9bK2trbGl84JhUJobW1FMBiMt6mxsTHh/06bgZhPQp/k5mVayYG+vj6MHTvW6eYQkQP2HurDdT9+ARv+/mrMmlzqaFv0pQUqKioAwPCAoU2LB9TBTt9b1djYiHA4HB/C0XKhjLbr72/8+PHw+XwoKytDTU1N/EAbCoXQ0NAQf4zVq1ejs7MT1dXV8V/y9fX1qK2tjbdVex4VFRUJ7bfbbjNaUKIvOVBXVxcvD6ANZwaDQaxatQqrV69GIBDAypUr4xWhtX2ampoSnmtdXR3q6+vjwUiq55L8mmlWr16dkFezatUqRKPRQYUWkyUXakwuOZBs9erVlpYE0YbOamtr48uMVFdXo6enJ+E9Gz9+fPz+6uvrMWfOnEHPWV92Ibkt2mvh9/vR1NSE6urqeKCZC9FoFOPGjTO8Th+naG2LRqMJw7hDcezYMZSWlgJAqZTyWKp9GTQRUcFwU9BEZJVW3ytbtZKam5sTAtna2lrU19cn5DXV1tZmPXnbq+wETZ5PBCei4tV97CS6+0/FL7/RfTzhXDNxzAhMHOu+sgNUvLSeOq33J5vFJRcvXhwv+gkMLtsQCoUGJYCTNexpIiLP+ue2/fiX519Pu993PnsR/rF6eh5aRGRNKBSK5wjV1dVlvSK3vuRA8pCjNjuNlKLqabJUcmDgNPD2f3E9OoedHpCW15izsy8Vry9dcT6qK89Nu9/EMSPy0Boi6/QLCedCqnwrp5dMyassH/8Lv6ep4xlg83eB6DtntvnOBxb8AKi8Ie/tLFYb93bhe892JKwTVl46EvdeXzlohpOdfYmIiAxZPP7b6Wkq7JIDHc8A678CTJwJfD0ErDykzifOVNs7nnG6hUVh494urHgsPGhh1cN9J7HisTA27u3KaF8iIiJDOTr+F25P08Bp4Eez1Qu05HGgRBcfDgwA624BujuAb+/mUF0OnR6QuDq4xXQlem2ZixcarwEAy/tyqI6IiAzZPP6zpwkAdj+muuTm35H4ggHq8vzbgejbaj/KmZZd75gGQYBa9qKr7yRadr1ja18iIiJDOTz+F27QtP0n6nziJcbXa9u1/Sgn1m570/J+dvYlIiIylMPjf+EGTVd9S513v2J8vbZd249yYvn8Cy3vZ2dfIiIiQzk8/hdu0PSJZSpLftsaNYapNzAAbHsI8E1T+1HO1M49H+WlI2GWgSSgZsbVzj3f1r5EAID+w8C7e9Kf+g871kQiyrMcHv89X6fJVMkwNa1w/VdU0tf821WXXPcr6gXbvxFY/CsmgefYsBKBe6+vxIrHwhBQeUkaLTi69/rKeGK3nX2JsOuXwB/+T/r9Pv2/gM+szH17iMh5OTz+e372XFJxy4ut1WmaBiy4n3Wa8oh1mign+g8n9iId2Q/8djnwhbXABF0F8DGT1Ik8qbGxEZFIBCtXrszZgrHZFo1G0dPTk3KxYsoxi8d/LtibvIwKK4K7AiuCU869uwdo/jRQ9wfgvNmONGH16tUAgKNHjwIAVq5ciebmZksr2HudFtxkeyHYaDSKCy+8EC0tLTmtop0t4XAYy5cvRzQaRWdnp9PNKW4Wjv9FtYyKJSXDgAvnO92KojesROCqivFZ35cIgPpyfHe3+vvd3cCkS/P+46i+vh61tbWDVpP3subm5viisunMmzcP48dn/3Pr8/k81WMTCASwdu3ajN57O683WZDl43/hJoITUfHoeEYVs9vwD+ryhn9Ql/Nc9X/9+vWDekLWrl2b1zZkW1tbm+V9a2pqctajVlZWlpP7zZVM13ez83pT/jFoIiJvc9FySdFoFNFoNGGbz+dDdXV13tqQLdFoND7cRrnH19sbimN4jogK08Bpleg5fWHicglT56nL624BNt8NzPhcXobqqqqqMGfOHDQ1NSX0OGnDLc3NzQgGg/G8n5qaGkQiEcyZMwd+vx8tLS2IRCJobGwEoHqpIpEIenp60N7ejqamJjQ3N6OsrAxPPvlkQmJ0KBTK6HaAykWaN28eIpEIfD4f6urqEAqFEIlEEtqzcuXK+OMEAgHU19fHe0ZuvvlmLF++HADQ3t4ev2/tcQEk3L/Z42q03DBtWK6np8fy+7B69Wr4fL7449bU1CRcp91nJBKJ94y1trYaPq958+YZbg8Gg2mfQ7JIJIJwOAwA2LlzJ6qrq+P/J2avt9ZjZbfdWvsoy6SUBXECMBaA7Ovrk0RUJCJbpbx3rJTv7DC+/p0/qesjW/PSnN7eXhkIBCRUxQxZVVUl29raEvZpa2uTfr8/YVswGDTcR39bv98vGxoa4pdbWlpkIBAY8u2S21hTUxO/3NbWNugx9PfR3t4u29vb4/ef/NxaWlpkU1NT/HJnZ2f8uaZ63IaGhoTXpLe3V/p8vkGvpZG6urqE29bU1MQv19XVJdxHZ2enrKqqSvu8zLaneg6dnZ2D3ueqqqqEtvn9ftnb2xu/bPZ6Z9pusqavr0/7zI6V6WKNdDs4cQJQA6DG5m0YNBEVm5daVFB0st/4+pPH1PUvteS1WdqBSwugWloSHz85sEm+vr29XarftGdUVVUlHAyNDsp2b9fe3i59Pl/C/k1NTbKmpkZKaX4Qb2trG3Q7oza1tLTIqqqqhMBAO7CbPW5vb++g5yCllIFAIG3Q1NnZOei2nZ2dsre3V3Z2dhq2WX+/Zs/LaHu6187o/dHaoqmqqkp4741e76G0m6yxEzS5bnhOCOEDUA+gyeGmEJHbjT5XnXe/oobkkmnLJWj75UkgEIgPfzU2NmL58uUJQ0T19fXxIbzW1taE6zTJicQ+nw8VFRUJ24yGrOzcbteuXSgrK0Nra2v8umg0innzDF7LJFZms9XU1KCpqQnjxo1DIBDAzTffjIaGhviQndHjhkKhjJOow+HwoNtq7Vy/fr1hm/1+P9ra2uLDZGbPK3l7Jq+d3+9HJBLB+vXrUVZWhp6enrTDjqFQaEjtpuzKSdAkhKgCUC+lHDTfUgihH/D1SSlXJ+2yGACnDxBRetP+5sxyCfqcJiBxuYRpf5PzpkSjUYRCoUEBUDAYxOrVqxGNRuMH9Lq6OjQ2NsYLIBrJdLZYJrczCtqMRCIR2wfntrY2RCIRhEIhBINBHD16NB7EGT2uPgjxAquvHaDyknbu3BmvY/Xkk0+m3J9J4e6T1dlzQoiAECIIoBbAoE9WLGDySSmbpZTNACKx/bXrawCsz2abiKiAacsl7N+okr4P7ABO9avzdbeo7Qvuz1u9JrPp4n6/P6EHxOfzoaqqCo2NjZg7d25e2mZk7ty5hgdms4O1lsRsVXNzM6LRKPx+P+rq6tDe3o5QKJTycQOBwKAZiFaZ3TYSiZg+ZjgcttSzlszua6fNjtMX/tTaava6hsPhrLebhiarQZOUMiylbIR5T1EjgFbd/q0A6gBACOEHEJFSRrPZJiIqcJU3qHWkuvcBP68GVk1R590dansel0tav379oJ4Ss+G3xsZGhEIh02VBknugkoMBs8DCzu0CgQBqamrQ3NyccL12ENeGkwBrvUxGj6W/b+0+Uz2uFmDpr9NmnaULprTbajPvtPvVgjFtOFSjDefZ6S3SpHvtktuqvS/67do27TU2er2z3W4aonRJT5mcoBK525O2+dTDDdpXAgjEbqOdWqBymvw2HpOJ4ETF7PSHUu78hUr83vkLdTmPent7ZVNTUzwJPBgMDpoFlszouvb2dllVVSUByLq6Otnb2ysbGhokABkIBGRLS4vhPpneTsozs9WampoGJaU3NDTIhoaG+Cw4LVlZu4/Ozk7DdkupEqO1+2xpaZHBYDAhETrV42rbtesCgYD0+/2D9jOS7vlo96tPkDd7XmbbUz2W/rXQP0YwGJR1dXWypaVFtrW1xWfB6WcYJr/embabrLOTCJ6Ttediw2wrpZRzdNsCsUBKJO3bC2C5VL1O2rYggJ36bRYe03ztOSIqDi5Ye86q1tZWVFVVZZz0TETZ4da158yyE3v018WSyKsABIQQYSklM+GIyFj/YXXSHNmfeK4ZM0mdHNbY2IiKiop4AUQGTETe4rqSA1LKEIA56fYTQowAMEK3aUzOGkVE7rTrl8Af/s/g7b9dnnj50/8L+MzK/LQpherqaoTDYaxevZqLshJ5UD6DJrNiFGUprktlJYB7M28OEXne3K8CFy9Kv58LepkAtcxK8oK+ROQd+QyaIoAqXikTZ8j5tOtsWgXgId3lMQAOZto4IvIglwy7EVFxyGrJgVRigVIEBrlNUkp7xT/UbU7FEra+DOBFAM8PtY1EREREZnIVNJklfQehSgoAiBe7bBzKA0kpH5ZSVgK4fCj3Q0RERJRKVofnYmUFboYKjPxCiCaoMgPNACClbBZCNGiVwQGMl6oYJhEREZGr5aROUz4JIW4DcBtUr9nFtuo0DZwG3v4v4Ph7akHPaX+Tt+UWCt3pAYkdb/agu/8kJo4ZicsvLMOwEpH+hnm+TyIiKgBDOJ7bqdPk+aBJY7u4ZcczwObvAtF3zmzzna/WscrjsguFaOPeLnzv2Q509Z2MbysvHYl7r6/EwlnlrrlPIiIqAEM8ntsJmvKWCO4qHc8A678CTJwJfD0ErDykzifOVNs7nnG6hZ61cW8XVjwWTghuAOBw30mseCyMjXu7XHGfRERUAPJ8PPd8T5Pt4bmB08CPZqsXdMnjQIkubhwYUCujd3cA397NoTqbTg9IXB3cMii40QgAk0pH4oXGaywPq+XiPomIqABk6XheVD1NtmfP7X5MdeHNvyPxBQbU5fm3A9G31X5kS8uud0yDG0CthtjVdxItu94x3Scf90lERAXAgeO554Mm27b/RJ1PvMT4em27th9Ztnbbm1ndL1f3SUREBcCB43nxBU1XfUudd79ifL22XduPLFs+/8Ks7per+yQiogLgwPHc80GTEOI2IUQHgB2WbvCJZSqrftsaNeapNzAAbHsI8E1T+5EttXPPR3npSJhlFgmoGW+1c8939D6JiKgAOHA893zQZDunqWSYmoa4f6NKEjuwAzjVr87X3aK2L7ifSeAZGFYicO/1lQAwKMjRLt97faWthO1c3CcRERUAB47nnp89p8lOnaZp6gVmnaYhYZ0mIiLKmyEez1nckhXBHceK4ERElDesCG7NkJZRISIioqLGOk1EREREWeb5oImIiIgoHxg0EREREVnAoImIiIjIAs8HTbaLWxIRERFlwPNBExPBiYiIKB/OcroBrsPaTba4oXaSG9pARER54uBxmkGTnmFV0fNVmXZWCR/EDVW63dAGIiLKE4eP054fnsuajmeA9V8BJs4Evh4CVh5S5xNnqu0dzzjdQlfZuLcLKx4LJwQrAHC47yRWPBbGxr1dRdEGIiLKExccpz1fEVyT0TIqmoHTwI9mqxd+yeNAiS6WHBhQC/91dwDf3s2hOqjhsKuDWwYFKxoBYFLpSLzQeE3Ohsnc0AYiIsqTHB6ni6oieFZmz+1+THX1zb8j8Y0A1OX5twPRt9V+hM37DpsGKwAgAXT1ncTmfYdz1oaWXe9YakPLrndM9yEiIo9wyXHa80FTVmbPbf+JOp94ifH12nZtvyK3fteBrO6XibXb3szqfkRE5GIuOU57PmjKiqu+pc67XzG+Xtuu7VfkRp1trevT6n6ZWD7/wqzuR0RELuaS4zSDJgD4xDKVfb9tjRob1RsYALY9BPimqf0I91xXiQmjh6fcZ8Lo4bjnusqctaF27vkoLx0Js2wlATWLrnbu+TlrAxER5YlLjtMMmgCVNLbgB8D+jSqZ7MAO4FS/Ol93i9q+4H4mgceU+0bh/ptmQQCDghZt2/03zUK5b1TO2jCsRODe6yvjj5ncBgC49/pKJoETERUClxynOXtOz7D+wzT1RrBO0yBuqJHkhjYQEVGe5OA4bWf2HIOmZKwIbosbqnG7oQ1ERJQnWT5OM2gaStBERERERcNO0MRlVKxiD5SnenS81FYiIjLhsmOv54MmIcRtAG5DLpPauSadp3KHvNRWIiIy4cJjr+dnz2WluGUqLljrxmleWuPNS20lIiITLj32MqcpFa5J56k13rqiJ3D9T17AkeMfmO4zYfRwPPutq3NaDoGIiIYgz8feolp7LqdcstaNk7y0xtt9GzpSBkwAcOT4B7hvQ0eeWkRERLa5+NjLoCkVl6x14yQvrfF24q+ns7ofERE5wMXHXgZNqbhkrRsneWmNt8Vzp2Z1PyIicoCLj70MmlJxyVo3TvLSGm8LZk6y1NYFMyfls1lERGSHi4+9rgqahBA1QogqIUSdEKLO6fakX+vmOWD2l4B9vwPe3KaS1wrA6QGJ7Z1H8fSeQ9jxZg/u+Zw31nizsh7dPZ+7BDve7MHTew5he+dRnB4ojIkQRESeNnBaHUdfblV1marvc3ydOSOumT0nhPABeF5KOSd2WUopLR+Jc1oR3KhWxEc+CggBvN99ZlsB1G4yq3F0w2XleObPXZ6ofVQIz4GIqGiY1WOqvAnoeCrn68E6voyKEKIKQL2UstbgOn0Pkk9KuVp3nU9KGRVCBACsNLp9isfM7TIq+qqkRzuB/1wFTF+osvsnXqLGWLetURHw4l95MnDSahwl/0dokevDt3wC484Z4Ykq28kVwXvf/wC3PW7+3B5ZFmDgRESUb1o9JrPjac2/AeeMz2lFcMeCpliwczMAH4C5Wq+R7vo66AIlIUQNgHlSykbdPjUA6gHUSimjNh47P2vPFWjtJi/VY7KL9ZuIiFzIJcdTx+o0SSnDsQCozWSXRgCtuv1bASTkLsW21QJwfg67ERfXjxiKzfsOW6rHtHnf4fw1KktYv4mIyIU8eDzNWyJ4LGfJL6WMJF3li/VQxcV6mHpiw3zu4uL6EUOxfteBrO7nJqzfRETkQh48nuZz9pzfZHsUgD82Y65Jt70MQHKA5TwX148YilFnW+v6tLqfm1ity3ThhHM4o46IKF88eDzNZ9BUZrK9J3bdegBtsZIDQQDLDXqlnJeufsTWNcA55wJnj/JEGQKtvMAV/jKMHXlWyn0njB6Oe66rzFPLsidd/SbNL/74FpaufRFXB7dwYV8iolzRygucPQoYPQnY9qDr6jGZSX2UzKPYkJyW7xRKt78QYgSAEbpNY3LQrMG02k3rv6KS1Obffibb/z8agPdeBgY+BH67XO3v4jIERlPzjWjBxv03zfJkorRWv2nFY2EIYNAMumSH+05ixWNhzqgjIso2o/IC+zcC65YmzZ576MxsdBdNqspVyYEaqJIBc3TbAgDak2svCSF6oXqVWmGDEOJ/A7g3eXvOZ89pjN54IDZt8k7XlyEwKy9gpFBqGVkNEgFvzxYkInIls/IC//FPZzocNDmox2TGDXWajIImH4BeAOP0pQSEEBLAHCll2OZjGPU0Hcxb0AScqd3U/QrwhyAwZS6w5AmDaZNLga6XgG+EgNLJ+WlbClam4I8deRbuWDAd088d6+p6THZp9ZvaOg7jF398K+3+j3wpgEWXejtYJCJyXN8h4F8/C5RfZnycfGIJcHAX8Lf/SwVTOajHZMaxkgOpxAKlCAxym+wGTLHbnJJSHtNOAPqH3kqbSoYBF84H3v4j8N9HVA+T4bTJO4D+d4FNd+W9iUasTME/dvJD/CnSg6sqxhdMwASoobqrKsbjzSPvW9rfi7MFiYhcZ9NdQH+X+XHyU3cCJ46q4+mF8101JKeXq6DJLOk7CKBGuxArdtlosq8lQojbhBAdAHYM5X6G5IP/Vufppk1q+zmMU/CtzwL871Mfcp06IqKh8thx0kxWgyYhRCA2860RQEAI0aRfNkVK2Rzbr04I0QCgQr+MSiaklA9LKSsBXD6U+xmSQCyzP920ySlz1GKEDsyq0y/CO3uqz9JtrE7V96J7rqvEhNHDU+5TIoA/vdWL76zbw1l1RER26RfhnRLL1kl3nAy4Z6acEdcs2JspIcRtAG6DCgAvzmtOkyZdKfi116hS8KdPndmex1l1RgnQJQIw6zgpliRoLRkeSD+jDuA6dURElhlNlho2Aqj4jEnur3NLkLkypylXXNHTpJUh2L9RvfEHdgCn+tX52muArt3qH+XrIWDlIXU+caaaRdDxTE6bpgUGyTPGUgVMAHDv9ZUFHTABwMJZ5XhkWQCTSkcmbDd72tpL9r1nOzhUR0RkRpslN3Fm0nGv8kx5Af1xct0tavuC+12by6TxfE+TJm8L9qZiGFkPByqucWRWnZVZcsk9ToVSXsAObUZdd/9JvHXkffxz6PW0t+GsOiIiA+lmya39TGzkRXdcymN5ASN2eppcU9wyU0nDc86qvAGY8TlVhuD4e8CuX6iZAKlm1f28Ws0qWPxo1ptjZZbcgASuuGAcbrlyGiaOGVlQ5QWs0mbUAcBXf2ltPsH6XQcYNBERJdNmyS3+d+Pj3t89oI570z4JzP0aMPrcvJYXGCrPB01SyocBPKz1NDndnngZAgD485Pq3Gy2wITp6vxop0qWy9I/jtZz8tZRa9PqPzLiLNw42/n6UW5gdVZd7/sfYHvn0aIMMomIBtHqFh7tVJfTzZI7+xzg0hrjfVzM+d6ZQpZqVl3HM8BPr1B/v/cy8Oh1Kpl8iDlOG/d24ergFixd+yI6uqyVrirkWXJ2WZlVBwB7DvZxRh0REaCOWz+arY5j772stnl8lpwZBk25NOM648V9tSS58tlZTQ43S/o2I6BymBbMnJTR4xWict8o3H/TLAgg7QK/wJl16hg4EVFRSk76bnwHGGNhEd4Z1znT3iHyfNDkiuKWZoxm1Z2IAs/9EzD9WpUkN3UeMGK0Ol/yuFqTZ/Pdtus4dUVP4O6n9lqaOg8U1yw5u8xm1RmRsdPdT+1FV/REzttGROQaA6fV5KfpC9Xxa+o8YFQpsOgBYP8mT8+SM8PZc/lgNKvu6yH1D5bswA6VJLf4V0DljZYf4o71e/Cb8CHL+xfjLDm7tNywh9pew863etPu/8XAZKxZPDv3DSMicoOOp1Uvk9HxrOMZ1UHQf/jMNodnyZkpqtlznqCfVbdpJXD45fRJcrsfsxU0VXx0tKX9rvt4Ob50xTQmMFugzar7zDsTLQVNZw8rwdN7DhXtLEQiKjK7H1PnRsezyhuACz8FBKcBky4Frl3lqVlyZhg05Ys2q25chQqaul8x7mnSkuROva9Kz6eZjqn1hvz+tW5LzRhxVkl8ej1Z0/mX45b2W7fzANbtVAv8siePiAqSNkvu+HvqOAWYH8+O7Ffn4yrOzCr3OM8Pz7liGRU7Mir8ZbzkitHyKKlMGD0cz37rapT7RmXjmRQNK0VCk3HJFSIqOC4r4JwtXEbFzUonA4tWx5LkjJZc2aP+AdPMqrMzU06bCXb/TbMYMGXA7ow6gEuuEFGBMV0aZWaKpVE2AYuCrg6Y7PJ8T5PG1YngRoYQsXfJMls9Hxwqyg67PXsaLrlCRJ7mwaVR7LDT08SgyUn6sWFtyZV0s+oqb8Id+EdLM+Uuv2Ac/rH6YiYlZ5F+nbpfvvAm9hxMX4T+Mxd/FL/8qjc6QomIBll/K9DxVPrjk0eXRimq4TlP05LDL61RJeWBtEuuRA68i4tP7EEJBoz30/nbiyfiqorxDJiySJtRd+PsyTjP4lCntuQKh+mIyFMGTqslvuwujXLhfM8ETHZ5PmhydXFLOywuueLv34G6N7+DF0bdjmtLUj9lq7O+KDNccoWIClYRLY1ih+eDJs8lgptJs+SKTFpyZdLHAnhk+L+YBk4TRg/HnQsuzk/bixSXXCGiglRkS6PY4fmgqWCkWHJFTr8WImnJlZKlTwAXXYt7RzyRMFTHmXL5xSVXiKig9B0CnmuILfVVHEuj2MFEcLfJYMmVJR/cjRcHKgFwppxTtATx4MZXsOdA+uTwv5s1CT9dNicPLSMisiFV0reHlkaxg8uoeFnlDTg9/e/w6p82IfqfP8Yn/7o9bfLd/xz7Au6cVYpR4yZjxhWfxrCz+Lbmm5YgPu4j6XOcAOCto+9je+dRzmwkIufpZ3JHD6pt6ZZGOfdSYGFhLI1iB4+uLnOmFhBwZcln8cnh281L1O9oBkqG4VMn/xPY9Z9q207j6uGUH4vnTsXvX/tL2v06uvqxdO2L7BkkImcZjW4A6ZdG+fQ/FczSKHYwp8lFkqt87xiYgXfFRAxsNUi+2/c05PP3QX6sOm31cMqfBTMnobx0pOXK4UwOJyLHGFX5/uomYPhoYOsDRZ/0bcTzQVOhlBzoip7A3U/thT7DbAAl+N6pW4DXN2HgCV3y3dvbIZ9eARgkiGPJ4yqB77lGldBHeTWsRODe61V+mZXAicnhROQIo4TvEaOBaVcCN/4UeH0z8MSSok76NsJEcJe4Y/0e0yrf15bswL0jHsd5sjvxinTVWS9bCnz+ZzloLaWTyZIrTA4norz53TeBPz9hfhzZtgb4/Q9UvpOmAJK+jTAR3CP0S3KcPcy8X2LTwOVoOzEXl5e8imvPB6pGv4mpb/w6bfVwfHhKVXMtskQ9N1g4qxzVlZOw480e3LdhHzq6+tPe5t2+E3h6zyFMHDOSCeJElBta0veHp9Rls+PI5XXA898HAreq5G+PLY2SKwyaHGK3J2IAJXhxoBKTyybjq2fHeo+MEvW0KaEAsO+36uRjcrgTtBl1f3/NRVjx63Da/fcc6MN31u0BwNIRRJQDRknfZgnfWpXv0x+opVEIQAHkNHlRcsK3VfEq39fcA4wpH1ydVUvqS6oezuRwZ9lNDgeYIE5EWWa7yvcaYMx56nhDcQya8swo4TudQVW+SycDi1bHqrMmVg9XSX1MDncTu8nhABPEiSiLMqryvQlYFFTHG4pj0JRnD25+DUeOf2DrNpNKR+KRZYHEoZrKG4DFvwK696mk7+A0VaV1/p1ASdLbWlICzL8D6H8X2HJfFp4F2WVnuRW9I8c/wIObX8tRq4ioKGy5D+jvGnx80I4jXXvUcWTVFHXe3aG2M6VjEOY05YmW9F0irPU1fO2TF+Cyqb7UScGVNwAzPqeS+nb9QuUvMTnctfTJ4d39J7G98wjW7TyY9nanPhxg9XAiss9K0re+yvfMLwBzv8bjQwoMmvIgk+nn8y4ow6JLLSQBlwxTVVn3/FpdZnK4q2nJ4QCwdX/6yuEAsOGlLmx4qYvJ4URknZ2kb63K91kjirLKtx2eH55ze3HLTJO+O7pSlooYjMnhnnPngosxYbS1teoAJocTkUVM+s4ZFrfMoa7oCVz/kxds5zDdtWgGbvrEZEwcay//Jf5Bmb4QmH+7GpL76RUqYFryROJY9sCASv7regn4RojJfg7RgmoAliYHCKgctxcar+FQHREN1ncI+NfPAuWXJX7vx48P16oc14mXqJ6nbQ+pKt9FnMNkp7il53ua3Mxu0nd56Uj8bFkAdZ+usB8wAUwO9yC7CeISQFffSWzedzi3DSMib2LSd04xpynL9FW+S0edbek2tXOm4AuBKdlJ9LWTHK5tHzYceLmVFV8dkkn18Katnfjg9ACrhxPRmYTv4++p73OASd85wqApizJJ+AaAa2ZMjCcHZ4WV5HAA2NGs9g0/qk4AE8QdoiWIXzD+HEtBE6uHExEA44RvgEnfOcLhuSzJNOEbyCDp2yqz5HAA2Pc08Px9wMeqmSDuIvdcV2krORxggjhR0UpO+F55CPjqJmD4OcDWB5j0nQMMmrIgkyrfmrsWzcCXr5yW9TYBMK4cfqofeHs78PQKVg93oXLfKNx/06x4FXgrtP+77z3bgdMDhTGxg4jSMKryPWI0MO1K4MZHgNc3A08sYaXvLHPV7DkhRF3szzkA2qSUrTZu69jsuTvW78FvwvYCjLwOqZh13349ZNx9e2CHShC8bCnw+Z/lvn00SKZDvY98KWCtvhcRedvvvgn8+Qnz7/Fta4Df/0DlO2l804AF9zP9Iomd2XOuyWkSQgQA9GiBkhBCCiHGSSmjzrbMnJb0ferDgfQ7A1gybyquqhif/+RdfXL48feAN7eqHCZWD3et5Orhv3zhTew52Jf2dj/e8jp8HxnO5HCiQmWlyjcAXF4HPP99IHCrSv7mRJ+syEnQJISoAlAvpaw1uK5Od9EnpVwd+9sPoBqA1rsUiW0L56KNQ5VJT8BfTw/gxtkOdYlqyeEA8EZInbN6uKvpq4dv2nvYUtDU0dWPpWtfZHI4USGyU+W7+xV1fvoD4NKa/LSvCGQ1p0kIERBCBAHUQgU8ydfXQQVKzVLKZgCR2P6I9TA1xvbzASiTUro2YLKb9D1h9HDcueDiHLbKBlYP9xy7CeJMDicqMKzy7Qo5yWkSQtQAWCmlnJO0vRNAtZQyotvWK6Ucl7RfC4AmKWXIxmPmJafJbpVvbYDkkWUBd/3qt109/BZVBO3bu9m96xBWDycqUgOngR/NVgHTksdZ5TvLXFkRPNZ75NcHTDG+WD6Ttl8DbAZM+WS3yvek0pHuC5iADKqH3w5E3wZe3eBMe4nVw4mK1asb1JDc/DtY5dth+UwEHzRcFxONXReO9VCFpZQhLZByyxCdlvRdIqz9Ys9qle9c0SeIb1oJHH45ffXwP/4IOP1XJhU6JJPq4et3HeCMOiKv0Vf5fvERtS1dle9JlwLXruJ3cw7lM2gqM9neA6AsFiS1AIgKFZj4pJSuiDYySfrOepXvXNESxMdVqKApXfXwQ7uA33xdbWOCuCPsVg//71Mf4uk9h7jkCpFXZFrle1wFq3znmGuKW0opw1JKIaUcFzul/GYXQowQQozVTgDG5KJdmVb6zlmV71xZ+ENWD/cYK8nhJQL401u9+M66PVi69kVcHdzC5HAiNxtKle+FP3SmzUUkn0FTj8n2shTXpbISQJ/udDDDdpnKtNJ3Tqt850rK6uH/M0X18IXA5rsTC6hRXlipHp5cIJyz6ohcbOC06mGavpBVvl0qb7PnYongvQASClYKISSAOXZzl4QQIwCM0G0aA+BgNmfP2a30XRC1cTKtHr74V0DljflpIyUwGj4uEYMDJr0Jo4fj2W9djXLfqDy0kIgs6Xha9TKxyndeubIiuJQyKoSIQPUsRZOus53sLaU8BeCUEOI2ALchS71mWsJ3d/9JlI4629JtPJH0bVVy9fAXH1F5TOmqh299ABhVxgREByRXD/9N+0Fsff1IytscOf4BHtz8GtYsnp2fRhKROS3pe+sD6nK6Kt+T5wJXruCEHAfkKmgyS/oOAqgBsBqIF7tsHMoDSSkfBvCwVqdpKPeV6Xpfnkn6tkpfPbzjWRU0pasefvhl4NHrmBzuEH318EO9J9IGTQBw6sMBbO88WhjBPpFXZVLle+wUVvl2SK4qgjcCCAghmvTLpsSqgEMIURerx1ShW0bFUZkmfAMeTPq2wyxBnNXDXavzL8ct7bfhpS4mhxM5KdMq30z4dkxOcpryKWl47uJMcprsVvnWu2vRDNz0icmYONZawUFPYvVwTymYqvVEhYxVvl3DTk6T54MmzVCWUbGb8A0USNK3HUZdyEwOdy27S64ATA4nyqtUSd9a6kO/rpo/E75zxpWJ4G6kJX2f+nAg/c4AlsybiqsqxhdnkUBWD/cUbckVOzl6TA4nyjFW+fY8zwdNmc6eyyTp+6+nB3Dj7CKug8Hq4Z6in1X36z+9jQ0vpc9bOntYCauHE+UCq3wXhKIcntOGLuw8cw5d6PQdAv71s0D5ZYNzmvY9DbTcGhuPv1M3Hr+G4/EO4hA0kYMS8kJjeUqH9wK//gIw7Wpg6TqD3NClQNdLwDdCLFqZY8xpShE0MUk2S5KTw+NfAjXABVczQdxlMpnswP99oiwwS/gGzvzIvGgB8Kk7mfTtkKIKmuzOnivKKt+5wurhnpJJcjjAXlaiIWGVb9ezEzS5ZsHeTEkpH5ZSVgK43Mr+FR8dbel+r/t4OZ5YfiVeaLyGAZOZyhuAb+8Bbt0AfPHnqkotkD5BfPdjeWkeJdKSwyeV2iuPoSWIE1EGtO+7VFW+B06r788v/lx9n357NwMml/J8IrgV+qVRXowctXSbEWeVFFaV71yxWj0cUMN3APB+D/DmNs4GcUDykivbO49g3c70a12zejiRTdpMufdj69GzyndB8PzwnMYsp4mLmeZRugTxp1cAH7x/Zhtn1DmOw9VEOZCculAyDKioYsK3SxXV8JwQ4jYhRAeAHcnXmS2NYhYwidjp/ptmMWDKROlkYNFqYP8mlfR9YAdwql+N2bfcClwwn0uuuMydCy7GhNHDLe9/uO8kVjwW5rIrRGaSl0ZZeQj4zN3A65uBJ5ac+V48sEN9T+7fBCwKMmDyiILtabIyWyi5x4m/orPE6FfWx6pNZtTxV5bTWD2cKEts97Yz4dsNimr2nCY5aPqfj7XjP/YeTnu7Ky4Yh1uunMZiftmmjeeH/rfKc0o3o67yJmDxo/luJcVkUuz1i4HJrB5OpLf+VqDjKfPvu7e3A79cqJK+q/438zpdoqiXUTk9ILG98yjeOvp++p0BfGTEWcVd5TtXtATxUWXqcroZddGDwMutXHLFIZlUD2dyOBESl0aJxiZVmH3fTZqlzkeVscq3RxVc0LTgn/+Av5yyfsBdPHdqDltDCCwD3ticfsmVd7nkitOGlQhcVTEere0HLO2/4aUubHipi8PaVLzsLo2izZQLLMt92ygnCi4R/L1jp6zdDiqHacHMSTlsHWHGdSoI2rZG5TDp7XsaeP4+le/EBHHXYHI4kQVGCd9f3QQMHw1sfWDw993AgKr07ZumvhfJkzwfNNktbgmcWR7i3usrOayQayXDVK/R/o2JM+re3q6SIqdfqxImp84DRoxW50seV9ufa1SJlZRX5b5RuP+mWfHZpOnI2Onup/aiK3oit40jcoO+Q8BzDbHvr8fPfH9NuxK48acpZsptVInfTD/wrIJLBJ/6D+tRMuIjKfflcIIDMl1y5bKlwOd/lp82UgImhxOZ+N03gT8/waVRCkRRJ4Knct3Hy/GlK6YxcdUJlTcAMz53JmHyza1A+FHzhMkJ09X5h6dYPdwhmSSHnz2sBE/vOcTZqFSYtKTvD2NpIKmWRnn++0DgVuDCT3GCi4P0K4Jk43upqIImLo3iMP2SK2+E1LlRwmTHM8Bz/6T+3vdbdWJyuCPsJoev23kA63aqfdmjSwXFqLc8XcL36Q+4NIqDjHrLh/q95PmcJqsmjB6OOxdc7HQzSHPNPcCYcmDbg4kJk1pyZflsJoe7iN3kcIAJ4lRAkpO+G98Bxkwa/P0FxBK+1wBjzlPfc+QIsxVBhvq9VHBBU3KnG5dGcSmjJVdORFUPE5PDXcducjjABHEqEEZJ36NKgUUPxL6/lnJpFJfpip7A3U/tNVzhYKjfS54PmpJLDkwcOyLh+kmlI/HIsgCHCNyo8gZg8a+A7n0q6Ts4Deg/DMy/M3H5AUBdnn8H0P8usOU+Z9pb5BbOKscjywKYVDrS1u2OHP8AD25+LUetIsqxLfcB/V2Dv5e076+uPer7a9UUdd7dobYzlcAxD25+LeUSakDm30uez2mSUj4M4GFt9tzmf/w0Xj36YdaSvijH9Aniu36h8peYHO5a+uTw7v6T2N55BOt2Hkx7O1YPJ8+xkvRdeYNK9A5OA2Z+AZj7NX4vOUSf8H32MGvfMRUfHW37cTwfNCXTElfJQ7QE8T2/VpeZHO5q+s/Y1v1/sXQbVg8nT7GT9H1kvzo/awSXRnFIJuVRAKDzL8dtP5bnh+eogDA53HNYPZwKDpO+PcUs4TudTCeHMWgi92ByuOewejgVFCZ9e0qqhG8zQ50cVnAVwfv6+jB27Finm0NDYdQ1zsrhrpZJ9/inLpqAL86ZwtxDco9Ulb61FIH+w2e2scq3o+5Yvwe/Cdv7wWyUIsCK4ORtdpLDte3DhgMvt7LyrkMyqR6+9fUj2Pr6EQAshEkO0hK+j7+nvkcAJn27nJb0ferDgfQ7A1gybyquqhjPiuBUwKwkhwPAjma1b/hRdQKYIO4Qu9XD9bRcJ5YHobwyWxOTSd+ulUmv9l9PD+DG2dkZPmVOE7mbWXI4AOx7Gnj+PuBj1UwQd5FMqocz14nyLjnhe+Uh4KubgOHnAFsfYNK3C2WS9J3t1UA8HzQlF7ekAmOUHH6qH3h7O/D0CiaIu1Am1cM1LIRJeWGU8D1iNDDtSuDGR4DXNwNPLGHSt4vYTfrO1WogTAQnbzDrRmeCuGtlWjulds4UfCEwhcnhlDupEr4B1aP0+x+ofCcNk74dZTfp206eJBPBqfDok8OPvwe8uVXlMLF6uGslVw9/68j7+OfQ62lv19J+EC3tB5kcTtlnpco3AFxeBzz/fSBwq0r+5gQTx9hN+r7u4+X40hXTcvaji0ETeYeWHA4Ab4TUOauHu5q+evgam8NuTA6nrLJT5bv7FXV++gPg0pr8tI8GyaS3esRZJTldFcTzOU1UpFg93HO+fOU03LVohuX9mRxOWcMq357jhqRvIwyayJtYPdxzJo4dibpPV+BnywIoLx1p+XZMDqchYZVvz3FL0rfhYzERnDyN1cM9SctTsFoIM5vF6ajIsMq35+Qy6duIpxPBhRBVAAJSytVOt4U8gNXDPcluIcx1Ow9g3U61LxPEKS1W+fYc7YdUd/9JnD3M2o+iXCd9G3FV0CSEqANQDWCn020hD2H1cM+6c8HF+MP+v+DI8Q8s34YJ4pQSq3x7TqblSXKd9G0kJzlNQogqIUSLyXV1ulOD/jopZTOAtly0iYoAq4d7TiaFMJkgTqZY5dtzMkn4BvKT9G0kq0GTECIghAgCqAXgN7i+DoBPStkcC5Aisf2Jho7Vwz1p4axyPLIsgEk2ksMBJohTElb59hy7Cd9AfpO+DR8/F4ngQogaACullHOStncCqJZSRnTbeqWU43SXtcDKVk4TE8EpjtXDPUmf07C98wjW7TyY9jZO5DSQS7HKt+fYTfgGcpPT6MpEcCGED4BfHzDF+IQQASllOF9toQLH6uGepC+EuXX/XyzdZsNLXdjwUheTw4sZq3x7jvYDqURY+6HztU9egMum+lwxezafieCDhutiorHrGDRR9rB6uKfZTRBncniRYpVvz8kk6XveBWVYdKk7Ptf5LG5ZZrK9R7suNqxXC6A69jfR0LF6uOfYTRBncngRYpVvz8k06bujK+WIWV65qiK4lLJVSlkdO7Wm2lcIMUIIMVY7ARiTp2aS17B6uCdlkiDO5PAiwSrfnpNJ0jcA3LVoBr585bSctCkT+Rye6zHZXpbiulRWArg38+ZQUam8AVj8K9WV//PqM9sX3wmUJP12KCkB5t+h9ttyH5PDHbRwVjmqKyfZqh5e8dHReWgZOWrLfUB/F7D43xM/v9rn/Ll/Svyc+6ap7Rxyd8yDm1+zVY/NrXmK+QyaIoBKCJdSRnXbfdp1Nq0C8JDu8hgA6afbUPFi9XBPsls9/MXIUUweN8oVSaOURazy7Wmzp/oszZSrnTMFXwhMce1nN29Bk5QyKoSIQPUsRZOus50ELqU8BeCUEOI2ALfBZUON5FKsHu5ZVpLDSwSw9fUj2Pr6EQDu/bVKNrHKtyfpy4j0vG+tl+maGRPzXuXbjlwFGmZJ30EA8QTvWE2mxqE8kJTyYSllJYDLh3I/VGRYPdxzrCSHDyQlTGiz6jbuTT+sRy7FKt+etHFvF64ObsHStS/iO+v24J9Dr1u6nZuSvo1ktbilECIA4GaowMgPoBlAe6z6t7ZPA1RPkw/AeCnlkIIm3f2yuCXZo30ZT18IzL9ddfUf3gv8+ovABfNVgrg+X2JgQCWYdr0EfCPEhFKHGE1ZLhGDAya9CaOH49lvXe1IBWEagr5DwL9+Fii/bPDncd/TQMutwEULgE/dqT6/3a8A2x4C9m9kDpNDuo+dxFO7D+GHz71q+7Z3LZqBmz4xGRPH2lsdYKjsFLfMSUXwfEoanruYQRPZwurhnqTv9n/8xbfxp7d6097m72ZNwk+XzUm7H7kIq3x7zprNr+HHW96wdRunh9FdWRE8V6SUDwN4WOtpcro95DGsHu5J+urhT+22VhbiraPvY3vnUdcmmJIOq3x7VmW5tU6Lf6y6CBdMOMdzEzY8HzQRDRmrh3va4rlT8fvX0i+70tHVj6VrX3T8Vy2lwSrfnqT1/m55tdvS/mXnDMeNs72X4uD5GWdCiNuEEB0AdjjdFioArB7uOQtmTkJ56UhLlcMBJoe7Gqt8e5I+6bul3Vrlnz0HorltVI54PqdJw0RwyprkBPEJ04GfXqECJiaHu5K2PAMAyxWHmRzuMmZJ3/HP47Wq6CwTvl1F++zZiSTc9tmzk9Pk+Z4moqzTqgp371NJ38FpQP9hYH6K6uH976oqxeQILrlSALQq38mfM+3z2LVHfR5XTVHn3R0MmBxmd2kUrVzI/TfNck3AZBdzmoiMsHq452Sy5MqpDweYHO4kVvn2NLtLo0wqgHxCzwdNrAhOOcPq4Z5jd8mVDS91YcNLXUwOdwKrfHuSvtzH2cOs/dC47uPl+NIV0wrixwlzmojSsVJgb/q1alghnm+xhvkWDuqKnsD1P3nB8q9g7Wv8kWUBBk75kJA3eIeusOwXgGlXA0vXMXfQhYwKy1rxxcBkrFk8OzeNyoKiKm6pYdBEOWVaPbwGuOBqkwTxW1Texbd3cxjBAUwOd6mB08CPZqsZckseZ5Vvj8gk4RvwxmeKQRODJsqFTKuHL/4VUHljftpICTL5Zez2X8We1/G0+gHCKt+eYbfnFvBW721RVQRnThPlTXL18BcfAQ7tSl89fOsDwKgyJq46IJPk8LOHleDpPYc8V6nY9bSk760PqMvpqnxPngtcuYITK1zAbsI3UBhJ30Y8HzRxGRXKK3318I5nVdCUrnr44ZeBR69jcrhD7CaHr9t5AOt2qn2ZIJ4lmVT5HjuFVb4dpiV9n/pwIP3OAJbMm4qrKsYX9A8O9s4QZWrhD1k93EPuXHAxJowebus2rB6eBZlW+V74Q2faSwASq3xb6aEFgL+eHsCNsyfjqorxBRkwAQyaiDJXOhlYtBrYv0klfR/YAZyIqh6m6deq5PCp84ARo9X5ksdVIvnmuxPzNSgvyn2jcP9Ns+IF9qyQsdPdT+1FV/RE7hpXqAZOqx6m6QvV///UecCoUmDRA7HPzVL1uTnVr87X3aK2LwpyhpyDtKRvO7mAE0YPx50LLs5hq9yBQRPRUNiuHn47EH0beHWDM+0tcplUDgdYPTxjr25QQ3Lz72CVb48oxirfdng+p4mJ4OQ4fYL4ppUqhyld9fDdj3FGnUP0yeHd/SexvfMI1u1Mv8goq4dnYPdj6jxdle9JlwLXrmLCtwsUY5VvOzwfNDERnFxBSxAfV6GCJrMk18N71fn7PcCb23iQcIiWHA4AW/f/xdJtWD3cBm2m3Ps96nK6Kt/jKljl20HFXuXbDtZpIsqmdNXDn14BfPD+mW2cUec4Vg/PsuSZciXDgIoqVvl2qUKt8m2HnTpNHNIiyiaj5PBT/WpGUMutwAXzOaPOZewmiGs/M7/3bAdODxTGj86sSZ4pt/IQ8Jm7gdc3A08sYdK3i3QfO4nmP3TimzYTvoHiSfo2wp4molww+rX9sWqT5Vb4a9sNMvnF/ciXAlh0KXubAGTQy8oq305as/k1/HjLG7ZuU6i9rOxpInJa5Q3At/cAt25QlY0HTqeYUXcH0P8usOkuR5pKysJZ5Xih8Ro8sfxKVJaPsXSbpq2deHrPIWzvPMpep013Af1dxv/nM28EvvQb9ffkuepz8e3dDJgcVFluv3NhUunIgguY7PJ8Ijhnz5Fracnho8rU5XQz6qIHgZdbuWyEg7QE8QvGn4OOrv60++850IfvrNsDoEirh2sJ38ffU/+/gPn/+aRZ6nxUGZO+HaQlfW95tdvS/l/75AW4bKqvoKt82+H5oImz58j1AsuANzabzyDa0awCpHd3Ab/5utrGBHFH3XNdJXa81WNr6rVWPbxofombLWCdbnmUwLLct40MZTIE3Xfir7hxNtMGNOydIcq1GbF157atGbxsxL6ngefvU/lOTBB3jUyrhwNFkiBulPD91U3A8NFqQV7D5VEeUnlMM65zps1FLNOk72JO+DbDoIko10qGqV6j/RsTZ9S9vV0lx5ouuXIt8FyjSrClvMukergE0NV3Epv3Hc5dw5zWdwh4riH2f/v4mf/baVcCN/40xUy5jSrxm8POeffvL76NHz73quX9i63Ktx2cPUeUL2bDGV8PGQ9nHNihlpaovAlY/GhemkiD6Qv//fKFN7HnYPosgMryMbjnupmFmQOy/lag4ynz/9tta4Df/yBxfUXOlHPUcy93YcWvw5b3L7b8PDuz5zyf00TkGfrlVo6/B2x/ROUxmSXOTpiuzo92snq4g/TVwzftPWwpaOro6sfStS8W1sFHS/o+2qkum/3fXl4HPP994Ly5wFUrOLHBQXaTvmvnTMEXAlMKM9jPEgZNRPmkzagDgGFnq7wQo8TZjmeA5/5J/f3ey8Cj1zE53AXsJogXTHK4US9puoTvq7/N9RUdlEnS94CU8R8IZIw5TUROMUsQ15Jsy2czOdxlirJ6eHLSd+M7wJhJwLYHmfDtUhv3dmEFk75zgkETkVOMEsRPRFUPE5PDXctugrink8ONkr5HlQKLHogtFbSUCd8u0xU9gbuf2gurITqTvu1hIjiR04yGPtIlh1+2FPj8z/LXRhpEyxe5b8M+S4UwZ08txVc/eaG3igT+7pvAn58w/n/UhpD7dcEgE74dd8f6PfhN2PqPqoLKu8tQUSWCsyI4eZ4+QXzXL4B9v01fPXzYcFYPd1jBVg/XV/keNlxtM/p/rLwBuPBTQHAaMPMLwNyv8X/RIfoZnqWjzrZ0GyZ9Z8bzQRMrglNB0BLE9/xaXU5XPTz8qDoBTBB3WEFVD7db5fvIfnV+1ggujeKQTBK+AeCaGROZ9J0B9s4Quck19wBjyo2TbFk93JUKpnq4aZXvc1JU+V4DjDlP/d9S3mWS8K3p6Eo5CkUmGDQRuUnpZGDR6liSLauHe4Xnq4enrPL9SIoq35uARUH1f0t5ZTfhW++uRTPw5SunZb1NxYCJ4ERulGn1cCaIOyqT6uGfufij+OVXL89D61JIlfANsMq3C9lN+AZcnkvnoKJKBCcqSMnVw9/cqnKY0lUP//AUq4c7KJPq4b3vf4DtnUedScjVkr4/PKUup6vyHbhVJX9zAoJjtMC8RFj7X/naJy/AZVN93pq16WKuCpqEEDWxP8sARKSUISfbQ+QoffXwN2IfhXTVw/f9Vp2YHO44qwniew72ObPkSiZVvk9/AFxaM/h6yotMkr7nXVCGRZeyZylbXJPTJITwAaiWUrZKKZsBNDrcJCL3MEsQZ/Vw17KbIK7NqNu4tyvXTcugyjcTvp2WadI3E76zKydBkxCiSgjRYnJdne7UoLtqMYCo7nJUCFGVi/YReY5Rgjirh7uenQRxGTvd/dRedEVP5K5RGVX5ZsK3kzJN+mbCd/ZlNWgSQgSEEEEAtQD8BtfXAfBJKZtjvUmR2P4AUAHgqG73HgC+bLaPyNMqbwAW/wro3qeSvoPTVDXm+XcCJUkf5ZISYP4dQP+7wJb7nGkvAVCB0wuN1+CJ5Vdi9tTStPsfOf4B7tvQkbsGbbkP6O8a/H+j/X917VH/X6umqPPuDrWdQ72OeXDza7bqgJWXjsTPlgVQ9+kKTBxrfUYnpZfVnCYpZRhAOJabNNdgl0YA1br9W4UQa2E+FFeWzfYReR6rh3uSliA+7iPDLe3/bt8JPL3nUPaSd1nl23NY5dud8pYIHstZ8kspI0lX+YQQAQCdUL1NmjIAyfsSEauHe9biuVPx+9f+kna/rC65wirfnsMq3+6Vz0TwQcN1MdHYdeuT9vFx9hxRCqwe7jkLZk5CeelIy5XDgSEmiLPKt+ewyre75TNoMhtq6wFQJqWMAnhSCFETy31qylvLiLyI1cM9Z1iJwL3XVwKwt+RKRgnirPLtOazy7X6uqtMkpWy1uq8QYgSAEbpNY7LfIiKX05J3N39XJe3qpUoQ/3m1Sghm9fC802bU2R1+0RLEf7psjrUbaAnfi/998P/BzBuBnntUle/XN53Z7pvGpG8H2U34BljlO9/yGTT1mGwvS3FdKisB3Jt5c4gKBKuHe87CWeWorpxke8mVt46+n756OKt8ew6rfHtHPoOmCKASwmNDcRofMkv4XgXgId3lMQAOZto4Ik9j9XDP0S+5MnxYCVb8Opz2Nh1d/amrh7PKt+ewyre35C2nKRYoRWCQ2xQrVWD3/k5JKY9pJwD9Q28lUQFg9XDPsZsgbpgczirfnsMq396Tq6DJLOk7CCD+kyaW8D2k5VKEELcJIToA7BjK/RAVDFYP9xy7CeKDksNZ5dtzWOXbm4SUmeTpm9yZqrd0M1Rg5AfQDKA9Vv1b26cBqsyAD8B4KWVW1pgTQowF0NfX14exY8dm4y6JvM1oqObrIeOhmgM7VHL4ZUuZHO6gTIZq/m7WJPz0nLXAn58wfn+1Idn+w2e2+aYBC+7nkKyD7li/B78JW/+RwoTv3Dl27BhKS0sBoDQ2cmUqJxXBkaL3SEq5OpuPKYS4DcBtcNHiw0SuwOrhnqNPEL9vwz50dJlnHZRgAJeXvIppXWEcKC/BVIBVvl2OVb69z1UlBzIhpXwYwMNaT5PT7SFyFVYP9xwtQfzvr7nINDn82pIduHfE4zhPdgPvA3gjdgWrfLsWq3wXBvbOEBUDVg/3HLPk8GtLduCR4f+CSR8LJFT5lsPPgWSVb9fpPnYSzX/oxDdZ5bsgeD5oYiI4kQWsHu45RsnhJRjAvSMeBy66FiVLn0io8i1iVb4lq3y7yr+/+DZ++NyrGd2WSd/uk9VEcCcxEZzIArPFW5kg7lr6YZ0rSzqwbvj95u/XtjWqyvfA6TPbmPTtqOde7rJUg0uPSd/55VgiOBG5HKuHe87CWeWonvFRvPqnTXj7938CPkTaKt+/G34jpn78aowaNxkzrrgWw87iV32+aUnf2yNHLe3/rc9U4KJzx7DKt8vxk0RUbFg93Fs6nsGwzd/FzOg7mKltS1Pl+8njl+LFF9RQXPnWP7DXIs8ySfqeeV4pq3x7AHOaiIoZq4e7m0mVb2lS5Xtg6xq8K87FjoEZ8c2G1cMpJ4aS9M2Eb29gThNRsdMOzNMXAvNvV0NyP71CBUxLngBKdL+tBgZUdemul4BvhJhYnEt9h4B//SxQflni+xB7v+T0ayHm36GG6rpfwcDWNcDrm7Dig+9g08Dlg+5uwujhePZbV6PcNyrPT6R4rNn8Gn685Y30Oya5a9EM3PSJyZg4dmQOWkXp2Mlp8nxPExENUeUNwOJfAd37VNJ3cJqqHj3/zsSACVCX598B9L8LbLnPmfYWiy33Af1dg9+H2Psluvao92vVFODn1Tj8xm7TgAkAjhz/APdt6MhP24tUZbm9H+zlpSPxs2UB1H26ggGTRzCniYjsVQ9ncnhuDZxW78OHp9TlNFW+I+ddj18c/xs83j0VA2l+B7919H1s7zzKROMs0lf5fv2945Zu85WrpmHRrHK+Dx7k+aCJy6gQZYmV6uFMDs8to5IQaap8+z86Gp+86vN4zMK09o6ufixd+yKntGdJplW+r/KPZ5Vvj/J8oCGlfFhKWQnAuE+aiOxhcrgzTJK+Dau4J1X5NqsebobJ4UPDKt/Fy/NBExFlmVH18BNR1cPEyuG50XcIeK4h9vo+rl7XUaXAogdi78PSlFW+jaqHpyJjp7uf2ouu6IlcPrOCxCrfxYuz54jImNFQUbrK4RXXALO/BIw+l7lO6Wi5S8ffU0OinVuMX19tSLT/8JltJlW+MxkuuuKCcbjlymksqmiBlr/U1nEYv/jjW7ZuyyFR97Ize45BExGZ0w7sWnL4ykOqhynZS08CT61IWr6DuU6mzJazMXt9T0TVrMaZXwDmfi1lQKod2O/bsA8dXf22msUDu7lMAlJW+faGoio5wOKWRDmkJYefNUJdjlWcTtDxDPDbeuBj1cx1siI5d2nlIeDvHlTXGb2+QDzpG2eNUO9Hih68YSUCV1WMx99fc5HtpjHXydjGvV1YkUH+0szzSnHj7Mm4qmI8A6YC4fmgiYngRHlglhw+cFr1mDDXyRqj3KURo1Xvke98S0nfVtlNDgdUnhMAfO/ZDpweKIxRiKHqip7A3U/tRSavBhO+C4/ngyYiygOj5PBT/WrYLvoOC2FaZVawsmSYGsrcvwl4YknKpG+r7CaHaySArr6T2LzvcNp9i8GDm1/DkeMf2L4dE74Lk+frNBFRnmiVwzd/VyV967EQZmpWC1Z+oUnlhr2+6cx23zT1umeQG7ZwVjkeWRbIqJbQj7e8Dt9HhhdtLo6WG1Yi7D135oUVNiaCE5E9Q5r1VYTJ4S6YhaivWv34i2/jT2/1Wr5tMQYBmSR9f+2TF6C6clLRBpleVlSJ4ESUZ1py+KU1wA0/YSHMVDItWHnDT9Trmybp2yotOfzG2ZPxf5d8AhNGD7d822JKDh9K0cpzRpzFhO8iwKCJiDLHQpjmhliwMlfKfaNw/02zIMBCmMkyLVrJ/KXi4fnhuaS15y7m8ByRA1wwBOUKOShYmSsshKnohy7fOvI+/jn0uuXbFuPQZSFicUsGTUT5V+yFMHNYsDJXir0QZqYL7tbOmYIvBKYUTOBY7JjTRET5V8yFMHNcsDJXirUQ5lAX3B2QkvlLRYo9TUSUXX2HgH/9LFB+mcpp0uoRDZwGfjQbmFiZuB1QSdDrlgJdLwHfCOU0pyfrCuD5nh6QuDq4BYf7Ttou4jhh9HA8+62rUe4blZO2ZZPWs/Zvf3wTmzrey+g+vPR8yRr2NBGRczIthPnJf1CFMH/zdVXTST9850YDp1U7f/ONvBWszJVMC2ECwJHjH+A763Zje+dRV1cR37i3C1cHt2Dp2hczCpi0xPn7b5rFgKmIsaeJiHLDTo5PxzPApruAvgNntrk5z8noudnK4cpv0rdVmeb4aNya56StHTeUo51bnxsNHRPBGTQRuYOV2WRaPtD0a1VvzcRLVB7QtjXA/o0ZV8POmXh7F6plYk70Ao/XFsxsQf1sst+0H8TW149Yvq3WS/XIsoCjwYX+OZwlBP6/Z/bh6Pv2l0L51EUT8MU5UwpqtiANxqCJQROR+xjl/ljJ+3n3z8DCVYA87UzAoQ/8xDBg00r7z8EluUt2dUVP4PqfvGB77bXx55yN798wCx9KmfeAY6i9ZRrmLhUPBk0MmojcKaGX5vZYL81i816abWuA3//AufIEZkOMKXvL7tD1lj3kzt4yG7ShLQAZD2/la2grG8Nwbukto/wpqqCJxS2JPMZqPlDKYbvngL+9Cxhfkd3eJ32v0tFO4D9XnRmGm3gJsP1h4D9/aCMvy525S3YNtfcml4GINhS3/71jWLN5P46d/HBI98fcpeJTVEGThj1NRB6iBSe//yHwzn8N7rlJNeS172ng6RXAB++f2ZaN3qfkYK5kmKonpX/8N7cBj15n3jP29nbglwuB868CPvNdV+cu2aUFJw+1vYadNhb81Rs78izcsWA6pp871vaQnT5PSRvya+s4nJWhOAC4tvJc/I9PXsjcpSLEoIlBE5E3mNU4MgtO0iWN1/wbcM541VOk9UABZ3qPzLb991Gg5X+kT+4u4PwlqzLNc0pWXjoS93zuEow7Z0RCIDSsRAwKkHrf/wD3/f+JwdGYkWehf4i9Spq7Fs3ATZ+YjIljR2bl/shbGDQxaCLyjuQ8p/gw2KrEYbB4wDJTLYCbHLCs/QzQ3QGc1h3MPzIBECXA+92ptw0brma3aYHQy62qXlSqYcOLFgCfurOg8pesykaek5Hy0pG44bJyPPPnrqz0Hll5PA7FkZ2g6az8NImIyETlDSrQ2PxdNTVfr/uVMz09b/+XGjr74i8GF8d8dQPQ9efEHqgdzcDz30+/7c9PAP9xZ2JxytHnDn58fXs/e49KUH9905ntvmlFETABwMJZ5XhkWSBrQ2Oarr6TaNr6Ztbuz8jnP3Ee/vbiiSwjQBlh0EREzqu8AZjxucFT+7c9eKb353isivPESxJvO3BaBVzTr00sA9D+y/TbAGDUuMH3O+1vVJ7UtjXGvVoHdgAfmehsKQSHLZxVjurKSfFhtD+81o3f7n7X6WalxGE4GioGTUTkDtqCv/rL67+ilhyZfzswYozantz7Y9QDZXUbYNyrpC2Bsv4rKk9pUBmBTUXTq5SKtuAvAFzlH48Zk8bih8+96nCrBuMwHGWL64ImIUQVgICUcrXTbSEiBxkN25UMA7Y+ACxddybwMeqBsroN0PUqPZjYA1V5A1D7qJqpt3/jmf2LaBjOjoljR6Lu0xU4f/xHsj5sl4mhzNQjMuOqoEkIUQegGsBOp9tCRC6gH7b7y2vAe3uB9n9L7P05EZv+ru8pMuo9MstT0vcqPbEkMbn7z4+r0gZzvgacOxP46MVFNwxnl37YLlu1k+zQQqPVNR9nzxJlne3Zc7GeoHopZa3BdXW6i75Meoti92H7tpw9R1QkjIpjJs9+MyoNkK5cgNHsuwIpTumkXM200ySXHuBQHNmVk5IDQogAgJsB+ADMlVLOSbo+IdgRQtQAmCelbLTTeAZNRJSWvnL3oDpLt+tmyt2XWBrAaJu+XIBRnSf2Kg1ZttaD09OCI30yOmfEUSZyWqcpFgytNAiaOgFUSykjum29UspxNu+fQRMR2WfUA/WRjwJCJNVpMtjGHqWcs1Kw0qxOU6pCmERDlfegSQjhA9ArpRRJ+0oAc6SUYSFEg9l96gMkBk1ElLHkHiirFcHZo+QIo6VRjCqCM0CiXHKiuKXfZHs0dl2Ys+GIKOeSyxZorG6jvNKXLLCynchp2Qqayky296S4bpBYL1Zt7O+IlLI1C20jIiIiGjJXlRyIBUmWAiUhxAgAI3SbxuSkUUREREQAStLvYkmPyfayFNcN1UoAfbrTwRw9DhEREVHWgqYIEE8I1/Np1+XAKgClutOUHD0OERERUXaCJillFCo4GpS/JKUMZ+MxDO73VCzL/csAXgTwfC4eh4iIiAjILGgyS+wOAqjRLsRKB9gqbJkJKeXDUspKAJfn+rGIiIioeFlOBNdVBK8B4BdCNAFol1I2A4CUslkI0aDVWQIw3m41cCIiIiK3sl3c0m2EELcBuA2q1+xiFrckIiIiq+wUt8xWIrhjODxHRERE+eCqOk3ZcOxYyiCRiIiIKM5O3OD54TmNEGIyWKuJiIiIMjNFSnko1Q6FFDQJAOcB6IeqDn4QqnZTv5PtIsv4nnkL3y/v4XvmPXzP8mcMgHdlmqCoYIbnYk/0EACo+AkA0J8uqYvcge+Zt/D98h6+Z97D9yyvLL2+nk8EJyIiIsoHBk1EREREFhRq0HQKwPdi5+QNfM+8he+X9/A98x6+Zy5TMIngRERERLlUqD1NRERERFlVMLPniCj3hBBtUsrqpG11uos+KeXqPDeLiCgvCm54jl/g7iaEaIj9OQ9AJHlRZ75/7iWEqAHQIqUUum110L1PsX3mcbFu58U+a9HYxR4pZavuOn7OXEb3nvgAjAewSkoZNbge4HvmmIIKmvgF7m5CiKD+vRBCtACAlLI2dpnvn0sJIXwA6gAEk4KmTgDVUsqIbluvlHJc/ltJGiFEG4B6KWVECBEA0K69b/ycuU8swG3WgqTY5y0opayPXeZ75hKFFjTxC9ylYl8CzwP4rO6LIQCgHUBF7Mud759Lxb601wPo1R18ffrLun0lgDlSynDeG0rae1WR9AMloL0f/Jy5j8mwd3wb3zP3KJhE8NgXuF//TxXjix2cyXn+2EmjvVd+vn/uFXv9dxlc5TfYBqghIbPrKPeCANr0G3QBkw/8nLlRmS51IQHfM3cpmKAJ/AJ3NSllVEo5Lqn3QXtfIuD752ZzTXqNykz270lxHeVQ7ADrgzqg1sVOQd0u/Jy5UyOAoBCiTQjhi71n9bHr+J65SCEFTfwC9556AKHYLyi+fy4khKiRUjY73Q6yTDuIlkkpm2PvXZuWPwh+zlxJShkCUA2gCkAvgJ26niW+Zy5SSEETeUisW7kKQK3TbSFjsV6LaIpdeky2l6W4jnJLO4jGh1NjB+QaIQR7JVwq9t4EAIwD0AygJWm2HLlEIdVp4he4twShkoWjsct8/9xnMYAKXd5EBRCf6RMBEIpd9umnRkMNDyXnX1B+RJLONVGog7LZ+8LPmbOC2ixiAPWxnsE2IUQI/G50lUIKmiIAv8C9QAjRBDUdOqrbzPfPZZKH5WK/huv09WGEENrQajTptpw554DYLFRADdPp3wNf7JyfM5eJ/ShJeO2llCEhxGqo3vj1sf34nrlAwQzPxf6ZDHNj+AXuHrEu56A2Xi+E8MemQ0fB98/tfAbbggBqtAux95e1Y5wVhnGuS5ifM0/phCoAHAXfM9comKAphl/gLhYryOaDKjFQFbvciDO/lvj+uZQW7Mb+bhFCVAFneqNis7QaoOoDsVKxsxqhyxWMvXetusRifs5cJBb4BGI5hHpzYvloAN8z1yio4pZAwtIBPgDjWTHVHbRCiEbXJVWY5vtHNERagUvtssFyRfycuUjs+3Fl7OJRGC+jwvfMBQouaCIiIiLKhUIbniMiIiLKCQZNRERERBYwaCIiIiKygEETERERkQUMmoiIiIgsYNBEREREZAGDJiIqarGq9FzMlojSYtBERMWORQKJyBIGTURU7Py6JUaIiEwxaCKiohVbYZ6LnhKRJQyaiKiY1QNocroRROQNDJqIqJhxaI6ILGPQRERFiUNzRGTXWU43gIgonViA44+dAKAVQI12vZRydQZ3Ww8gaPJ4NQDmATgKIALgZgCrpJQMsoiKGIMmInK1WA0lv5SyNXa5F0CFlLJeCNEEYC6ATIImw6E5IUQdgFopZbXucg2A5Zk+ByIqDAyaiMjtqqSUzbrLPgBtsb8zqrFkNjQXC9CaAIzTbY4AiEopo5k8FhEVDgZNROR267U/dJW7QwCQHMjErq+BCnT8AFpNEr3NhuaaYrfR329AezwiKm4MmojI1ZICmCoA4RS9Pi1SyjkAIITwAXgewByD/cxmzVVBBVR61TjTs0VERYyz54jIS6oB7DK6Ijbk5tMuxwIrX/K6ckKIKpgPzcHg/qvAniYiAoMmInK5pKCnCkC77roa3XVzAUSTbh6FGl7Tq0XqgpbxHqhYgAUpZVgIEeDCvkTFjUETEblWLCjqFEL4Yn/3xE7a8FuZbnefdl2SsqTLhkNzsW1aLpR2//U4E4hVsRAmUXFjThMRuVkYQDOAxVABUTWARiFEGQAkzaqLYnCABOgCKbOhOZ1aAPVCiE6oGXO1QogWIURDmtsRUREQUkqn20BENGSxnKYWKWWFblsngGqthyhW1ynIHiMiygSH54ioICRX644Nr0WTAiSuNUdEGePwHBEVklohRBDATqhlUGq1K2JDcywdQEQZ4/AcERUFDs0R0VBxeI6IikUZAyYiGgr2NBERERFZwJ4mIiIiIgsYNBERERFZwKCJiIiIyAIGTUREREQWMGgiIiIisoBBExEREZEFDJqIiIiILGDQRERERGQBgyYiIiIiC/4f/ECVFZZoGcMAAAAASUVORK5CYII=\n", "text/plain": [ "
" ] }, "metadata": { "needs_background": "light" }, "output_type": "display_data" } ], "source": [ "shifted_correlator.show(comp=symmetrised_correlator, logscale=True)" ] }, { "cell_type": "markdown", "id": "232e88af", "metadata": {}, "source": [ "## Effective mass" ] }, { "cell_type": "markdown", "id": "83dc751c", "metadata": {}, "source": [ "The effective mass of the correlator can be obtained by calling the `m_eff` method" ] }, { "cell_type": "code", "execution_count": 10, "id": "c686f7e0", "metadata": {}, "outputs": [], "source": [ "m_eff = symmetrised_correlator.m_eff()\n", "m_eff.tag = 'Effective mass'" ] }, { "cell_type": "markdown", "id": "4a9d13b2", "metadata": {}, "source": [ "We can also use the priodicity of the lattice in order to obtain the cosh effective mass" ] }, { "cell_type": "code", "execution_count": 11, "id": "5acde8cf", "metadata": {}, "outputs": [], "source": [ "periodic_m_eff = symmetrised_correlator.m_eff('periodic')\n", "periodic_m_eff.tag = 'Cosh effective mass'" ] }, { "cell_type": "markdown", "id": "c658b000", "metadata": {}, "source": [ "We can compare the two and see how the standard effective mass deviates form the plateau at the center of the lattice" ] }, { "cell_type": "code", "execution_count": 12, "id": "1d6ea22a", "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAmkAAAGLCAYAAACcFQXGAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjMuNCwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8QVMy6AAAACXBIWXMAAA9hAAAPYQGoP6dpAABSjUlEQVR4nO3de3hc1X3v//eSXAwB22PZli9gXEa2wy0ByzYEwiUJUkyay0laGQrJCScPP8vl15PSNKkVQds8tA3uqJCW5rRE4rQp+aUQLJ8kJ00bHIlcgELB8jgXbrWtcczFxsLII9sEnFizfn/svcd7RnPbGo20Z/R5Pc880t57zd5rz0ijr9blu4y1FhEREREJl7qproCIiIiIjKUgTURERCSEFKSJiIiIhJCCNBEREZEQUpAmIiIiEkIK0kRERERCSEGaiIiISAjNmOoKVDNjjAGWAEenui4iIiJSVWYB+22BhLUK0sqzBHh5qishIiIiVeks4JV8BxWklecowEsvvcTs2bOnui4iIiJSBY4cOcLSpUuhSE+cgrQJMHv2bAVpIiIiMqE0ccBljIlMdR1EREREPKFqSTPGtPs2I9bariLlI8B17mYTEAE6rLVJ93gv0JbjqXFr7WpjTAvQ5ztfAmi11ibGew8iIiIiEyE0QZoboKUDM2NMmzEmZq3tKPC0GNBtrY27z+kGeoFW9/iw+/2w7znXAw+530eA1e73SQVnIiIiEhahCdKADk4GV1hrtxpj7nP35xMFWoC4uz0I+Fvj+qy1/d6G2/L2uhfUuRJey5uIiIhIWIQiSHODp2iOlqyIMaY5K6hKs9a2Zu1qAvp9x7dmHe8s0jInIiIiEgqhCNJwWsRySbrHcgZpfsYYr1UtO3DzjjcD23Mcus4Y43WHri0UxBljZgIzfbtmFauXiIiIyHiEZXZnQ579wwWOpbnj2XqBjQXGlXXmaFlLAAPW2q3usUF3XFs+ncCI76FEtiIiIlIRYQnSymKt7bHWrgY6jDGbso+7rWy5nhfP6krtB9oLpOPYDMzxPc4qq+IiIiIieYQlSBvOs7+hwLFcYkAsR1DWhtNqVpCvFS5fUHfcWnvEe6A1O0VERKRCwhKkJSBnQtkIeYIrY0zEGNOb9RyvbEtW8etxZn5mP/+wP6BTQlsREREJi1BMHLDWJt1Esg04kwX8x/JNGvAmCvifE3G/Zgd2UXK3yA1kjWGLFrmmiIhUyNCRtxg6erxoucZZM2mcfeok1EhkaoUiSHPFcLolvWS27fhypHmzN621PeAEUsaYnqwg63qc1QT6yRRhbPCXNMb0ZZXrpHBeNhERqZB/eepF7nlkd9Fyt16zgs+0rpyEGolMLWOtneo6pLmD/pM4QdU8fzoML2iz1jb59kVwAitPBN+yUL5yh4FrcrWQ+SYaNAE7vCCwxPrOBkZGRka0wLqISJmyW9L2DB3jDx/6CX97/cUsbzwjvb/WWtKSySTDw8NEo/myUUmtOXLkCHPmzAGY445xzylMLWkUWqvTDZ56svYlKaHly1o7dzzXFBGRydM4+9ScwdfyxjO48Mw5U1Cjkzo6Tv6pmTdvHtFolLa2Njo6OojFYuM+b39/Px0dHSSTSQYHB4s/IYBEIkF3dzdNTU00NDSwfft2kskkTU1NbN++nd7e3gm9HjivUyKRqMi5p6NQBWkiIiIAoynLz15OAvCzl5Oct3g29XVm0usRj8fZsGEDsViMlpaTc9ISiQTr168nHo+XFaS1tLRw3333sX79+omobobVq1ezd+9eIpEI4ARQHR0dtLe3ZwSd49XT00N7e3vGvrVr1zJv3ryyzy2OsMzuFBERAeDhZw5wRewH3PatZwC47VvPcEXsBzz8zIFJr8v69evHBGgA0WiUjRs3Tnp9StXf3080Gk0HaACRSIQ1a9YAlBVYevr6sod1Q1tbG5s2jUlXKuOkIE1ERELj4WcOcMvX4xwYeStj/6sjb3HL1+OTGqh5rU3ZAZqnpaUl1OPIGhqKLtgzLslkMt2tKZWl7k4REQmF0ZTljn99jlzT2SxggDv+9Tlaz180KV2fW7duzRugebK7Dbu6utKBWyKRyGhV6unpSQdOiUSCSCSS0V0Yj8dJJBIkEgkGBwfp7i60SuHJ669duzbjfP39/XR3dzMwMJBRv0QiwebNm4lGo3R2dhKJRHI+338vkUgkXee2tjbAaaXz6umdv7Ozk0QiwYYNGwDYsWNH+p5jsVh6nFpbWxvxeJxrrrmGaDRKb29v+vUqVBc/bxwfwH333UcikWB4eJgdO3bQ3d2dfp0feughOjs7aW5uzngN4nFnDuH27dtpbW3NeI/zvUfF3ruKsdbqMc4HMBuwIyMjVkREyvPEnkN2Wcd3iz6e2HNoUuoD2FgsVnL59vZ229fXl94eHBy0LS0t1lpre3t7bXd3d8Yx79w7duywkUgk47nRaNTu2LGj4PVaWloyntPW1pbe7uvrs83NzRnlm5ubM8oXen57e3vGvbe1tWVs5zq/tz8ajRbdl/26FqpLLt45s1+zTZs2pbd7e3vH1LGlpSXj2tFo1B4+fDhdPtd7VOi9G6+RkRGL87/HbFsgzlB3p4iIhMLQ0beKFwpQbjIlEgm2bNmS0SoTjUYZHh6mv99J3dnb20symUwf85dNJpNjnluoOzEejzMwMJDxnNbW1pJa34o9P5FI0NPTk9EKGIvFSmo5ytX929LSwvDwcLoFC8gYKzeee2loaCCRSIx5zfyam5vTr7enu7s74z6i0Wj6/YH871Gh966S1N0pIiKh0DirtNxnpZYrVzQaLZoWI5FIpP/Q5wpQotEofX19xGIxuru7mTt3Ls3NzVx//fUZQZA/aPG2h4fzL109MDBAQ0MDW7duTe9LJpOsXbu2pHsr9Px4PD6mPuWOvWtvb6e7u5vu7m76+/u57rrryr6XXK9ZU1NTxr7s19ALfrds2UJDQwPDw8PpMm1tbTnfo+bm5oLvXSUpSBMRkVC45JwGFs85lVdH3so5Ls0Ai+acyiXnVGZAfLa2traMwCGX/v7+kscm9fX1kUgk6O/vJxaL8frrr6dnWY53kL83Tmy8cj2/2D3n4gWr+WzcuJHVq1enW+pytUQFvZfxvGZdXV0ZOeIeeuihjOP53qNC710lqbtTRERCob7O8IUPnw84AZmft/2FD58/afnSvD/C/u4wv2QymQ4U1qxZk7N7Mh6Ps3btWnp6ekgmk0SjUdrb29mxY0fe85Yi3/VKnXFZ6Pm5ugmLndvflZlLNBolGo1mDMAvpS4TyZuV6k+0691nPB7P+x5N9HsXhIK0CfD6gRfZ89PHiz4O7d831VUVEQm1ay9czL2faGbRnMwuzUVzTuXeTzRz7YWLJ7U+vb29dHR0jPmjnEwm6enpSbf+NDc309LSktEK5XUbtrW1pcv7+VuesrvlcgVJfs3NzbS1tWWcM5lMFg2WSnm+F4x0dXVlHPMHTf4xc/5WtEL13rhxIx0dHWNazMZ7L8Ves+xtr7x/v7cvkUjkfY+KvXeVFKq1O6uNt3Zn35f/Jy2Hvla0/JNLN3DZzXdVvmIiIlVuNGV5aPuL3PatZ7jzYxdy/dqzp2TFAY8/lYU37ilXN2dHR0f6+ODgYLo1zvsj70/j0N7enk5j4XWbxmIxNm/eTFdXF83NzXR2dhbsBuzo6GDevHnpVBltbW3pFBXxeJxNmzYRi8Xo6OhIn9M/pirX8wudO9dr0tTURHt7O/F4PONesgf+J5NJNm/enLebsNj1PNnXyfWaRaPRMWUikQhdXV0MDg7S2tpKJBJJJyX2r/iQ/R5t2bIl5/7sMXFBlLp2p4K0MnhBWuKFnzP6VjK9/+C+XSwc6OLgmk0sXLYyvT+yYCnzlyyb/IqKiFShZ14Z4UNffpzvfvqKKV+7U2QiVeUC69Vq3uKzmT37wox9y+MHYNlKll90xRTVSkRERKqZgrQJNnriBK/s2kkytYI3du3knAveRf0MvcwiIsUMHXmLoaPH09t7ho5lfPU0zppJ4+zJScMhMpXU3VkGr7tzZGSE2bNns3Pb/Sx56ossTB1MlzlYt5D9l97OqnU3TV1FRUSqwN/07eKeR3YXLXfrNSv4TOvKouVEwkrdnZNs57b7uejJW2HFOrjqc9B4Hgw9z4JH72LBk7eyExSoiYgU8PFLz6b1/IVFyzXOmjkJtRGZempJK4PXkjb8+uv86n+9iwVNq6i74UGo82U2SaVIPXgDQ4M/YcFtz6rrU0REZJortSVNedImwK6BR1iYOkjdVZ/LDNAA6uqou+qzLEq9ygtPbZuaCoqIiEjVUZA2Ad5M7ne+aTwvdwF3/5uHX5mkGomIiEi1U5A2AU6LLHG+GXo+dwF3/2lzz5ykGomIiEi1U5A2AVauuYaDdQtJPXoXpFKZB1MpUo/ezat1izj30nVTU0ERERGpOhrFPgHqZ8xg/6W3s+DJW0k9eAN1V302Pbsz9ejdsHsbBy67h0WaNCAikt/RV51HMbMWOQ+RGqeoYYKsWncTO8HJk7a7Nb1/qG4RBy67R+k3RESKGfgq/Pivipe7+vPw3s7K10dkiikFRxmyk9mCs+LA49+6l9N/fj9vvOMmrvjYLUq7ISJSiuyWtEO74Jsb4Lfvg/m+5LVqSZMqp2S2k+j1Ay8ytDeZ3p4x83QiHOP4zNPZ++x/pvdrgXURkQLyBV/zV8KSiye1KolEglgsRk9PD5FIhPb29vSxZDLJli1biEaj7NixI+M53d3dNDU10dDQwPbt20kmkzQ1NbF9+3Z6e3sntI4dHR0kEokJP6+Eh4K0CTDY103Loa+lt5cD1MHy+GcgfrLck0s3MP/muya9fiIiEkw0GqW7u5uBgQHWrFlDLBbLON7d3U1ra2vGvtWrV7N3714ikQjgBFEdHR20t7fT0dFRVn16enoyAkWAtWvXMm/evLLOK+GmIG0CNLVuZM9b1xctt2LB0kmojYhIDUiNwv6dzvf7d8Kid0Bd/aRXo6GhIe+x9evXp7/v7+8nGo2mAzSASCTCmjVrAMYEeUH19fWNCdLa2trKOqeE37QJ0owxEWttshLnnrf4bGbPvrASpxYRmX6e+w58/3ZIvuhsf/cP4fEvwfu/COd/ZEqrlkwmGR4eJhqNsmbNGpLJZDowKxTQlXO9zZs3k0gkJvzcEn6hCtKMMf5/EyLW2q4i5SPAde5mExABOrxgzBjTAvT5yieAVmttwrcv0DVFRKSCnvsObPkkrLwWfuef0umMeOxuZ/91X5vSQG1gYABwukObm5sBpxXN6xr1d2smEgk2b95MNBqls7OTSCRCR0cHa9euJZFIjBnr1tXVRSQSSQd7bW1t9Pf3k0gkSCQS6XN3dnaSSCTYsGEDQHpcXE9PD7FYLD1Ora2tjXg8zjXXXEM0GqW3t5doNApQsB5+/f396eved999JBIJhoeH2bFjB93d3fT09NDQ0MBDDz1EZ2dn+jXx7j8ed8b8bN++ndbWVlpaWtLHved6Zb165Ns/LVlrQ/EA2oFNvu02IFbkOd1Ac9Z2X9Y5mt1HdCKumfX82YAdGRmxIiJSptET1v7Nhdb+y/XWjo5mHRt19v/NO5xyk6SlpcVGo1G7adMm297ebgHb19c3plxfX59tbm7O2Nfc3JxRtqWlJWO7ra0tvd3e3m5jsVjGMW8717m9/dFotOg+/3mL1SMX75z+Mt5r4unt7R1Tx5aWloxrR6NRe/jw4XT57u7u9LHBwUEbi8Xy7q81IyMjFrDAbFsgzgjTigMdwFZvw1q7FSeIKiQKtPi2B7O2ARLW2rj1tZ6VeU0REamEfU84XZxXfhbqsv481dXBlX8EyX1OuUnU0tJCLBaju7t73DMp4/E4AwMDGS1Jra2tdHd3k0gk6OnpYdOmTeljsVisaOuR1yqWXdfh4eF0CxaQMU6uUD3yaWhoIJFIZDwn+9rNzc0kk8mMfd3d3Rn3EI1G6e/vT2/39vamnxONRtPnz7d/OgpFkOZ2W0ZzBFIRY0xzjqcAYK1ttZndk01Af77yE3FNERGpkGMHna+N5+U+7u33yk2BlpaWjLFn/qCjkIGBARoaGti6dWv6kUwmWbt2LfF4PCOQAsZMQgiivb09HXT19/dz3XXXpY8Vqkch2XWJRCI0NTVl7BseHh5zD8PDw/T09LB161aGh4fTZbxJD3PnzmX16tV0dXXR3Nycd/90FZYxaWP/HXAk3WPxPMfTjDFeq1pr1qHrjDHeT85aa603YKDsa4qIyAQ6Y6Hzdeh5WJojaBh6PrPcFIhEImPGXQWRa0bm1q1bc5TML5FI5GxF82zcuJHVq1enW+lytUQFnRk6nkkRXV1dGfnhHnrooYzjfX19JBIJ+vv7icVivP7668Risbz7p6NQtKQB+d794QLH0tzB/73AxqyWsQQwYK3d6nZlDhpjvDbdwNc0xsw0xsz2HsCsYnUTEZESLbscImc7kwRSqcxjqRQ89iWILHPKhYCXdqMUa9asyRnQJRKJnF2F3rFc/F2ZuUSjUaLRaMYA/FLqMZGSySQdHR0Z3cPePcbjcXp6ekgmk0SjUdrb29mxYwf9/f15909XYQnSymKt7bHWrgY6jDGbfPvj1lr/T3M/0O52dY5HJzDie7w8zvOIiEi2unonzcauh+EbN8JLT8Pxo87Xb9zo7H//X05qvrTsLjy/jo6OkoM0ryuvp6cnvS+ZTBKPx9MBSVdXV8YxL3CKRqPp7/2taLkCO8/GjRvp6OgY02JWqB6FZL8O2dfO3vbK+/d7+xKJBMlkMqMO4Nxnvv3TVSjW7nTHgO2w1pqs/YeBDW4rWCnn8VJuNOWZKIAxxgKr3c1A1zTGzARm+nbNAl72r90pIiJlys6TBk4L2vv/ctLSbxRaFsrriksmk1hr02kq4vE4mzZtIhaL0dHRkR5Pdf3116cnBXR0dDBv3rx0qg1/EFXsGEBTUxPt7e3E43E6Ojro7+/PGIPm8fKr5esmLHQtv+zrxGIxNm/enL63zs5OotHomDKRSISuri4GBwdpbW0lEokQjUbZuHFjRhJgf6qN9vZ2tmzZknP/eMfnhVWpa3eGJUiLAIeBudaXcNYLqLJaw/zPuQ8noEq6+6I4Mzw3AluAve7zE1nXWY3TFRromjnqMGaBdRERmQCpUYh/zUlk+6G/heZPTsmKAyKVUGqQForuTjdISpBjLFiBYMmbKOB/TsT96rWiDWS1qEW9c47zmiIiMhnq6mHJKuf7JasUoMm0FJbZnQAxnGSyXZCeDJBO3ezN3rTW9oATSBljerKCsOuBuLW2331OH5k6/ecsdk0REZlER191Hp5DuzK/emYtch4iNS4U3Z0ed9B/EqdFbJ4vXUY6gLLWNvn2RXACL08E37JQvnOCk0NthxfklXLNEuqr7k4RkYnyw83w478qXu7qz8N7O4uXEwmpqhqTVq0UpImITKDslrR81JImVa7UIC1M3Z0iIjKdKfgSyRCKiQMiIiIikklBmoiIiEgIKUgTERERCSEFaSIiIiIhpCBNREREJIQUpImIiIiEkII0ERERkRBSkCYiIiISQgrSREREREJIQZqIiIhICClIExEREQkhBWkiIiIiIaQgTURERCSEFKSJiIiIhJCCNBEREZEQUpAmIiIiEkIK0kRERERCSEGaiIiISAgpSBMREREJIQVpIiIiIiGkIE1EREQkhBSkiYiIiISQgjQRERGREFKQJiIiIhJCCtJEREREQkhBmoiIiEgIKUgTERERCSEFaS5jTGSq6yAiIiLimTHVFfAzxrT7NiPW2q4i5SPAde5mExABOqy1SV+ZTe63a4GEtbbDd6wF6PNtJ4BWa21i/HchIiIiUr7QBGlugJYOzIwxbcaYmD+oyiEGdFtr4+5zuoFeoNXdjmUFZb3GmF5r7Xp3VwRY7X6fVHAmIiIiYRGm7s4OYKu3Ya3dCrTnLw5AFGjxbQ96224rW0tWN+ZmoM0YE/XtS1hr4wrQREREJExCEaS5gVQ0R6AUMcY053uetbY1q0u0Cej3bUfdhyfh2y8iIiISWmHp7swXNCXdY/FiJ3Bbx1pwuzrdcWlz81zHHwxeZ4wZdr9fW6h71RgzE5jp2zWrWL1ERERExiMULWlAQ579wwWOpbnj2XqBjUW6LTcC/b4yCWDAWrvV7V4ddMe15dMJjPgeLxerm4iIiMh4hCVIK4u1tsdauxro8M3mzOB2m7YA633Pi3uTDlz9QHuBdBybgTm+x1kTUH0RERGRMcISpA3n2d9Q4FguMSCWNTHAf2y1Pz1HNl8LW87uV2vtcWvtEe8BHA1QNxEREZGShSVIS0DOhLIRMsePpRljIm5KDf9zvLItWWW7cbpCk1nPP+wP6JTQVkRERMIiFEGaGzwlyDH+LKs70s+bKOB/TsT9mg7s3PFqMa+VzBgT9c0YHcgawxYtck0RERGRSRGKIM0VA9q8DTe48ieijfpXJHADqZ6sIOt6IG6t7Xef04YTuEWNMS3udgdObrQkvtUGXJ3+a4qIiIhMFWOtneo6pLmD/pM4gdW8rNUC2nGWfGry7YvgBFaeiFsm6R47nOs61lqTdU1wcqztsNb2BKjvbGBkZGSE2bNnl/o0ERERmcaOHDnCnDlzAOa4Y9xzClWQVm0UpImIiEhQpQZpYeruFBERERGXgjQRERGREFKQJiIiIhJCCtJEREREQkhBmoiIiEgIKUgTERERCSEFaSIiIiIhpCBNREREJIQUpImIiIiEkII0ERERkRBSkCYiIiISQgrSREREREJIQZqIiIhICClIExEREQkhBWkiIiIiIaQgTURERCSEFKSJiIiIhJCCNBEREZEQUpAmIiIiEkIK0kRERERCSEGaiIiISAgpSBMREREJIQVpIiIiIiGkIE1EREQkhBSkiYiIiISQgjQRERGREFKQJiIiIhJCCtJEREREQkhBmoiIiEgIzZjqCvgZY9p9mxFrbVeR8hHgOnezCYgAHdbaZKnnDHpNERERkckQmpY0N1iKWGt7rLU9QMIYEyvytBgw4D6nw93XW+o5x3lNERERkYoz1tqprgMAxphBoNVam/DtO2ytnVvgOX1An9f6ZYzZBMSstaaUc47nmlnXnw2MjIyMMHv27CC3KyIiItPUkSNHmDNnDsAca+2RfOVC0ZLmdltG/cGSK2KMac73PGtta1b3ZBPQX8o5x3tNERERkckQljFp0Tz7k+6xeLETGGOiQAvQWuI588l7TWPMTGCmb9esYvUSERERGY9QtKQBDXn2Dxc4luaOLesFNvpaxoqdczzX7ARGfI+Xi9VNREREZDzCEqSVxR34vxrocMelVcpmYI7vcVYFryUiIiLTWFi6O4fz7G8ocCyXGNBnjNlawjkDX9Naexw47m0bYwJUTURERKR0YQnSEuAM9vfnOMPJe5Y9sB+vLHAfsMH3HK9sC7ClyDkDX3OiHNq/j+RrLxUtF1mwlPlLllWyKiIiIhJSoQjSrLVJY0wCpxUrmXUs36QBb6KA/zkR92uilHOO45oTYvf3vsxlL91XtNyTSzcw/+a7KlkVERERCakw5UnzEst25dmOAi1u0lnvOTFfElvcRLQt7vi0Us5ZcLuEOo8rT1p2S9rBfbtYONDFwTWbWLhsZXq/WtJERERqT6l50kITpEE6GW0Sp0VsXlYA1o6z5FOTb18EZ8alJ8LYZaHynrOU40XqOyHJbPf89HGWf+uD7PnYv7H8oivGfR4REREJv1KDtFB0d3oKtWC5LWg9WfuSQMGgqlirmNbqFBERkTCqiRQc1Wz0xAle2bWTgdQKXtm1k9ETJ6a6SiIiIhICE9LdaYz5HBC31v6g/CpVj3K7O3duu58lT32RhamD6X0H6xay/9LbWbXupgmsqYiIiITFZK/dOZ+TKS1+O/ugMebiCbpOzdi57X4uevJWFjStgpv7ofMVuLmfBU2ruOjJW9m57f6prqKIiIhMoYkck+YtWt5qjElmHdsIXD+B16pqoydOsOSpL8KKddTd8CDUubHy0rXU3fAgqQdvYPFTdzJ6zcepnxGqYYMiIiIySSYqAvByliWAZmCe71gEWD1B16kJLzy1jQtSB+Gqz50M0Dx1ddRd9VkW7W7l2ae2ccG7Pzg1lRQREZEpNVFB2kPW2v8DYIzZYK3NyNRqjPmdCbpOTXjz8CvON43n5S7g7k+XExERkWlnosak2TzfewYn6Do14bS5ZzrfDD2fu4C7P11OREREpp2Jaklrcmd4JnDGpGUfXw+sm6BrVb1zL13HwUcWsuDRuzLHpAGkUqQevZuhukWce6leMhERkelqQlrSrLV/DRjgEmAusDzrEZ2I69SK+hkz2H/p7bB7G6kHb4CXnobjR+Glp53t3ds4cOltmjQgIiIyjU1YFOAGahhjrrHWPuI/Zozpm6jr1IpV625iJzh50na3pvcP1S3iwGX3KE+aiIjINDfha3caY+YAMWCPtfYuN29af6FkbdVqItbuHD1xgse/dS+n//x+3njHTVzxsVvUgiYiIlLDJjuZrd/ngV5gL4C19pvAdRW4Tk2onzGDM1euYk3dbs5cuUoBmoiIiACVWWB9wFr7iDHmmgqcuyYc2r+P5GsvpbcP7tsFqcXOV5/IgqXMX7JssqsnIiIiIVCJIO0c96u/H3UN8L8rcK2qtPt7X+ayl06mklsOUAfL45+B+MlyTy7dwPyb75r0+omIiMjUq0SQttMYMwC8bozZiLMCwcYKXKdqrfjAp9nz2keLl1uwtPKVERERkVCa8IkDAMaYKNDubnZba/dO+EVCYCImDoiIiMj0UurEgYqMUrfWJnAmEIiIiIjIOEz47E5jzB+7aTdEREREZJzKCtKMMZ9zH+/z9rlJbXcqUBMREREZv3F3dxpjtgARd/M2N4ltP/AVYCdaCkpERERk3MppSdturX2/+2gA1uIEZ38NDALzJqKCIiIiItPRRK7dGcfJ8qUJAyIiIiJlKqclLW6MuXiiKiIiIiIiJ407SLPWPgK0GGPuNMbMmsA6iYiIiEx75Uwc+ArOwukRoMMYk8CZONAH9BdKziYiIiIihZXT3Zm01jZYa+twlp/8PGCALuCwMWbbRFRQREREZDoqZ+LAHu8bd9mnvcD/AXDTcSgFh4iIiMg4ldOSZty1K8ew1o5Ya3eWce5JZ4yJTHUdRERERDzlTBy4D+g0xlw0UZUxxrT7HptKfM4m99FrjIllHes1xtgcjx3u8RZvH04X7aC7OLyIiIjIlCpn4sBfAZuATRMxacAY0w5ErLVd7nabMSZmre0o8JyM425Q1mutXe/uGgZa3a+e64GH3O8jwGr3+6S7MLyIiIjIlCunu/P1ApMGdhtjNgc8Xwew1duw1m4F2vMVdrsnW7K6KTcDbb7WsD5rbb+1Nu4m20249Y77npNwjytAExERkdAoJ0hLGGP+H5yg5/9Ya3/PWrscaAA6yWy9KsgNtKI5AqWIMaa5wFOjZE5QSPj2e4GeX6fXUiciIiISZuPu7rTWejM5rwEe8e0fwdciVqJ848CS7rF49gFrbRKYm+c8Y1rF3GBve45rXGeM8QLKtYW6V0VEREQmS9lrd7orD6QZY37TWvuLgKdpyLN/uMCxXDbijInL1XXZ6Rur5kngdncCGGMajDHd1tqNuU5ujJkJzPTt0koLIiIiUhHj7u40xlycJwXHXGPMHxtj3ldGvcZTn2agBcgOxMg3Y9M3Vs3TD7QXSMfRCYz4Hi+XU2cRERGRfMoZk3YbziSB140xDxljPmeMudhau9Na+9dAobFk2fKNX2socCxbDFjtdoNmayNHF2g2Xwtcvu7XzcAc3+OsEusmIiIiEkg5QdpDwO/hrN85ALwfiBtjRo0xrwNrA5wrATkTykYoIbgyxnQDG/MEaOCk3RjMek7EGHPY38pWLKGttfa4tfaI9wCOFqubiIiIyHiUE6SdY629z1r7iLX2r62173dTcqwDHrHWXl/qidzgKkGO8WdZ3ZFjuPnVYl4rmDEmmmNGaJTcLXIDWePXvFmhBa8pIiIiUmnlBGnzc+201vYDG9z0HEHEcLolgXTw5U9UG3X34dvXhtPaFnVXD2hznzMmlQfOTFF/PZM4yXf9Ov3XFBEREZkq5czu3G6M2Qa0WWszuv2stSNB18K01va4yzu14wRV87LSYbTgBFA9kO6a7M1zruzZmUlytKRZa7t8y0814SS/7QlSbxEREZFKKCtPmjFmHjBijOnFaZUawGnFiuIEPUHPmTfRrBs89fi2kzirHJRy3ux8aiVdU0RERGSqlJUnzW396ge6cQIo6x5KkCMVhoiIiIiUZiKS2SZwFjHHGHOOu29vuecVERERmc7KDtL8FJyJiIiITIxyZneKiIiISIUoSBMREREJIQVpIiIiIiGkIE1EREQkhBSkiYiIiISQgjQRERGREFKQJiIiIhJCE5onTSrj0P59JF97qWi5yIKlzF+ybBJqJCIiIpWmIK0K7P7el7nspfuKlnty6Qbm33zXJNRIREREKk1BWhVY8YFPs+e1j6a3D+7bxcKBLg6u2cTCZStPlluwdApqJyIiIpWgIK0KzF+ybEw35vL4AVi2kuUXXTFFtRIREZFK0sSBKjN64gSv7NrJQGoFr+zayeiJE1NdJREREakAY62d6jpULWPMbGBkZGSE2bNnV/x6O7fdz5KnvsjC1MH0voN1C9l/6e2sWndTxa8vIiIi5Tty5Ahz5swBmGOtPZKvnLo7q8TObfdz0ZO3wop1cNXnoPE8GHqeBY/exYInb2UnKFATERGpIerurAKjJ06w5Kkvwop11N3wICxdCzPPgKVrne0V61j81J3q+hQREakhCtKqwAtPbWNh6iB1V30O6rLesro66q76LItSr/LCU9umpoIiIiIy4RSkVYE3D7/ifNN4Xu4C7v50OREREal6CtKqwGlzz3S+GXo+dwF3f7qciIiIVD0FaVXg3EvXcbBuIalH74JUKvNgKkXq0bt5tW4R5166bmoqKCIiIhNOQVoVqJ8xg/2X3g67t5F68AZ46Wk4fhReetrZ3r2NA5feRv0MTdYVERGpFfqrXiVWrbuJneDkSdvdmt4/VLeIA5fdo/QbIiIiNUbJbMsw2clswUnH8fi37uX0n9/PG++4iSs+dota0ERERKpIqcls1d1ZZepnzODMlatYU7ebM1euUoAmIiJSoxSkiYiIiISQmmGqwKH9+0i+9lJ6++C+XZBa7Hz1iSxYyvwlyya7eiIiIlIB0yZIM8ZErLXJqa7HeOz+3pe57KX70tvLAepgefwzED9Z7smlG5h/812TXj8RERGZeKEK0owx7b7NiLW2q4TnbHK/XQskrLUdvmMtQJ9vOwG0WmsT5Vxzsq34wKfZ89pHi5dbsLTylREREZFJEZrZnW6wlA6SjDFtwFp/0JXjObGsoKwXwFq73ncOLyBL+oOz8V4z6/mTPrtTREREqlupszvDFKQNMraV67C1dm6e8hHgEeAarxvTGNMM7ACarLUJN+jqz9fNGfSaOZ6vIE1EREQCqaoUHG7AFc1u6QIibuCVT9R9eBK+/ZW6poiIiEjFhWVMWr6gKukei2cfcFvHslu8vPP4A6/rjDHD7vf+rszA1zTGzARm+nbNynMOERERkbKEoiUNaMizf7jAsVw24nRvekFaAhiw1m611m4FBo0x3WVcsxMY8T1eDlA3ERERkZKFpSWtbG4XZQuw2ttnrc1uDesHuo0xJU0MyGEz8CXf9ixCGKhl51XLR3nVREREwissQdpwnv0NBY5liwGrC+VCcycTgNOdGfia1trjwHFv2z1X6GTnVctHedVERETCKyxBWgJyJpyNkDm+LCe3C3Oj/7nuxIC9OIFbwrdvQq4ZZtl51Q7u28XCgS4OrtnEwmUrT5ZTXjUREZHQCkWQZq1NuolmG3AG7vuPjRnA7+fmOov5ArEoJwOtgazZm1H/Ocd7zbCbv2TZmG7M5fEDsGwlyy+6YopqJSIiIkGEZeIAON2Vbd6GG3z5E9VGs1YH8JLPRoCoMabF3e7AWXkgiW+1AVen/5zFrikiIiIyVUKTzBbSSzwlcQKveVmrCbQDHdbaJnc7AhzOdR5rrfE9z1s2qgnYYa3tKfWaJdQ39MlsR0+c4PFv3cvpP7+fN95xE1d87BbqZ4SiAVVERGRaqroVB6pR2IO0ndvuZ8lTX2Rh6mB638G6hey/9HZWrbtpCmsmIiIyfZUapKlJpUbt3HY/Fz15K6xYB1d9DhrPg6HnWfDoXSx48lZ2ggI1ERGREAvTmDSZIKMnTrDkqS/CinXU3fAgLF0LM8+ApWud7RXrWPzUnYyeODHVVRUREZE8FKTVoBee2sbC1EHqrvoc1GW9xXV11F31WRalXuWFp7ZNTQVFRESkKHV31qA3D7/ifNN4Xu4C7n6vnFYoEBERCR8FaTXotLlnOt8MPe90dWYbej6jnFYoEBERCR/N7ixDWGd3jp44waE7z2dB0ypnDJq/yzOVIvXgDQwN/oQFtz1L/YwZY1rS8q1QoJY0ERGR8ml25zRWP2MG+y+9nQVP3krqwRuou+qz6dmdqUfvht3bOHDZPSxy86VphQIREZHwUZBWo1atu4md4ORJ292a3j9Ut4gDl92j9BsiIiIhpyCthq1adxOj13ycH2etOLCowIoDoydO8MqunSRTK3hj107OueBdWqFARERkCuivb42rnzGDM1euYvmzf8KelasKBlzeCgVXpw46yVme/RMOPv/lMSsUaDaoiIhI5SlIEyDYCgWaDSoiIlJ5CtJqUK7ZmqQWO199vJauMSsUeLNB3RUKUg/e4KxQcM3HqZ8xgxUf+DR7XvtoxvlzzQZdsWBpRe9TRESklilIq0HZLV3LAepgefwzED9ZzmvpeuGpbVyQOui0oOVboWB3K88+tY0L3v1BzQYVERGZBArSalB2S1fecm5LV9AVCkRERKTyFKTVoFwtXYUEXaHAT7NBRUREKkMrDpQhrCsOBBV0hQKPNxt0Yepget/BuoWaDSoiIlKAVhyQkgVdoQA0G1RERKTS1JJWhlppSfPkahl7tW4RBy69LaNlTGuDioiIjJ9a0iSwUlco0GxQERGRyqsrXkSmE2+FgjV1uzkzzwoFmg0qIiJSeWpJk8DKmQ0qIoVpoo2IeBSkSeAVCs69dB0HH1nIgkfvyj0m7dG7GapbxLmXrhtzLaXsECnMP9Fm1MILqbN508zkNHucc+tepN445TTRRqT2aeJAGWpl4sCT//i5kmdfXub+UfDP7sw1G/Snl92TMdnAe04pKTtEpjPvn6YXf/YoF+z9KgtTQ+ljB+saefacT3H2O69SS5pIFdPEASlZ0BUKwJlksBOcoGt3a3r/UN0iDuQJ0EpN2aHuHpnO5i9Zxks//xHvGezK+fvynt1d/LTxzPSkm6EjbzF09DjgtFS/9uwPSR19lbpZi1hwwXvTLdWNs2bSOPvUqbqtCeO/30Jq5X5lelNLWhlqpSWtHKMnTvB41mzQ7O7LoCk7xtOyJ1Irgv6+9PzbE/zfx+NcbPbwP2f+G4vtyZa3A6aR/3X8g/zELue/XdFM+wcvn4I7Kixo0OXdL4C18EtO5QT1zGCUt/EWxu0ODuv9ioBa0mSSeLNBlz/7J+zJMxs0aMqO7Ja9fHnV/C17YaKWQPEL+vMQ9Pfl4/WPsHbGt7ioPgHLM1veFj56F3+x+5/56WiUlfWHgfAFLd9+LJ4Ougrxgq6P1z9C+8y7eHh0LX/x60+wtO41GkkyRISXUgv40xlf59r67bxR/zngcrW8SVVTkCYVFzRlhz+vWnqiAWdw4vgbVTHRQCssiF/Qn4egvy+nXvopljz9NWhal9nytnQtdTc8SOrBG1g8uJNTL/1U2fdSCV7QVYwXdJ3+7g08/LYWvv3ww/Se9kWW+FoO95tG/vz4jdD6F1y76hwA/uWpF7nnkd1Fz3/rNSv4TOvKouVEJlO4/9pNImNMxFqbnOp6VIOgs0HHm7LDm2hwdeqgk9Hv2T/h4PNfDv1Eg2pvCZSJFfTnIejvywvP/ZQLUkPFW96e+ykXvPusibuxCXL6uzfAOz8MwGjK8m9PP4fd+f9hVv13PnjJ+dTXOf2Xp89a5JQ5fSE/evQf+YdT7hnTcrjo0bv4h933cNujp9B6xW3UA//9gpl8uHFO+vw/SPySf9m+n4+vXcL7om9Lnz+yYObk37xIEaEK0owx7b7NiLW2q4TnbHK/XQskrLUdpR43xrQAfb7tBNBqrU2M8xamheyWgeUAdbA8/hnw9Vp4LQPjSdkRZKJB2GiFBfEL2jIc9Pel6pNLz1oEsxbx8DMH+Ivv/Jylx/bTyGqGBvYTe2EBf/qRd3DthYvTxZ8efI0/GP1nZ2Z5npbDP9hzP08PbuCyFY3M/68HmP/jv+Lh0bXc8etPcoB5QCOx7Sf42va9fOE3vsa19dvh6s/Dks4peQlE8glNkOYGaOnAzBjTZoyJZQddWc+JZQVdvcaYXmvt+lKOAxFgtft9UsFZaYLOBg26gPvoiRMseeqLBT+EFz91J6PXfHxSuj41xkz8xjubstSW4aC/L7WQXPrhZw7w7Qe+Qu/MB1hyiq/78teN/PkDN8KNv5cO1EZ/8R9OF2eBlsMlu1tJ/OI/YMXHYM2neLj+Pdzy78MYUryr7rn0GLbtqXO55def4d7WhnT3aMb7m7I8u3+Ew7/8NXPf9htcsGROuuVNY9hkMoQmSAM6gHQuB2vtVmPMfe7+MYwxEaAlq5tyM7DDGBMFhgsd9wVkCXVzBpOrpaiYICk7gg6crrRKjzFTEFhdvIHu3mzKi/yzKbfnnk0ZtGU4yO9LOcmlw2A0ZfnRt/+pcPflt0+h9fzbqK8zNJqk88QiLYdeudHTF3LHfzzL++ue5gszH8g5hu2O/7iK1isWUk/wMWyamCCVFIogzQ24ojlasiLGmGZrbb6pP1H34R1P+PYPFzmuVrNJVuoC7kG7byod5JQzxqyUFRYmc6LB6IkTvPDUNt48/AqnzT2Tcy9dN+kTMfzv1+joCV4d/Dknjg0z44wGFjW9g/p6pz6TFZQG/fkJOpsyaMuwV59Zi5pIfqiH7QOPcNqLj/Lm2VexYs01zKqfwZ6fPp6uz+GhV3j2nP/Be3Z35W15e65pE+8cemXSX8/RlOXnr40y/GaKhtPqeMeCet8YMKf+Qbsvm6JN8DhFWw6bok0APL13mHcefZR7CwSBtxyFp/dezGVN8/j4pWdz+in13Pm9F/Le420fOJePrnJaJjUxoTAFseUJRZCGEzTlkiQzyEpzW7/m5jlPothx377rjDHD7vdrC3WvSvlKSdkRtPsmaJAT9I/yeGebltq9VemJBtkZ7C/wZ7B/ZPIz2Hvv187RJpbMGOHt9tDJ+jw3n/0n5rCqfnDSZr8GXYYp6GzKoC3D2T/Pbwfnk3r/T+A7f5fe79Vn9/e+zPu813PPQGbLm5nPgdEo7xuM8eSvhif19cyVIuNvUgv4099wUmR49Q/afVn/m+/mzdPPYmaBlsPjp5/Fab/5buc1OPIGX5j5ACzP/359Yc+DbD9yMzCPefYwX33MmQRVR4pL6l5I1//p1LlY6vjqY7u5+eK3AYv5+KVn03r+wnQV9gwd4w8f+gl/e/3FLG88I72/cdb0nJgw3YLYUoPS02zxMhCeIK0hz/7hAsdy2Qj0Fxhbln08gRPQxQGMMQ3GmG5r7cZcTzbGzAT8v2mzAtRNShS0+8Yf5IyOnmB3VsuD1zLjBTnjbbkKMts0SPdWpSca7P7elzn1Fz/gPfWJ/Bnsd23ltd9835ggthItXSs+8Gl+8EiDk1F/eY7XZ/c2ftDUwTuv+XjZ915qffa89tGTQax5MX3MvwyT9/MTdDZl0JbhoGM+vfKzgOToCV7Ier9m1c9gj6980H9SgpZf8YFP89X4B/nPJ36UM0XGHcdv5MClt/PhZud/5qDdl0PHfs1bl/0pS/t/L2fLodm9jddavsKpx35N4+x6lv/y5yUFgct/+XPgbJ7u7+XAsSbW5ekeveP4jWw7dglP9/dy2e/8AY2zT83ZArS88QwuPHPOmP3TrWVpugWxpQalG9+1uGgZCE+QVjZjTDPQwsmJAEWP5+hG7Qe6jTEdecapdQJfmJAKS15BB057QY4XRH0odTDd8nDwuw8VbLkqJaiDYEFX2CY+RFtvgX/+v0VbfkzrLUB5LV2ldKfObTyTC/b+c8HX5/zB+5nbuInJEHQZpqBBV9CW4aBjPrPLv735PQXLB205DPpPzdxFZ/NfP/l63u7Fe3ffw20/OYVP/pZTz6Ddl86YwENcXncD7XseZoGv5fCQmcd9v76BJ/7tEP9tJE77By/nvFm/dA4Web+8ckNntrBux0MF63/Lr25l6MzrnecffdV5uE49dIwLzF5OPTQHzMkgxJvFWu1j3oLWxx/EjqYsP3s5CcAvf3WC8xbPTnd/14pSg9LT7HFuK+F8YQnShvPsbyhwLFsMWF1gEkCx41hrE8ZZUyRnFyvOxIMv+bZnAS+XWL9pK2heNajs2qBBg7qgQVfYJj4ceumF0lp+XnqBhWc3BW7pCtqdGrbXJ+j7GzToCtvA/qAth0H/qQk6xixo92VG8tusVQ0b7evc/hsPACeT39a5+dWKvV9eucZ580vqHt037/ed/QNfhR//lfP6WMNrqXNpnxHhta1f55y6F6g3biWv/jy8tzNwy1LYugvHW5+HnznAHf/6HAdG3gLgtm89w5d/sIcvfPj8jBQr1a7UltUjR/KuBJUhLEFaAnImlI1QwgB/Y0w3sDFfAJbruDtZYS9O4Jbw7cvLWnscSP8LYUxt/QdQKUHzqnlKnWgwnparIEFd0KCinLxVpUw0CCpofYK2dAXtTq10Xq+g3bVB39+gQVfQluFKC9pyGPifmqApMurqOe2Dm7FbPpm3+/K0674GdfVAZvLbQrzktyy7HCJnYx+7C/O7Y98v+9jdmMgypxxwSf0L1JdQ/4X1LwCNsOZT8PYP8PCeN7nj0SQHjqXSxRefUccXropw7fLTnJY0CNyy5A/qRlOW7z/3Kn//w0F+/71NvP/8RRkpQaDy3anj6b58+JkD3PL1eHZMzasjb3HL1+Pc+4nmmgrUJlIogjRrbdJNJNuAM1nAf6zgom5ufrWYL9CK4uRbixc6jhP8DWSNX4uWck0JJugYm+yWtxkzTyfCMY7PPJ29z/5nev94/8gGDeoq3b3lKXXMW9AxQoEz2Ad8PYN2p1Y6r1fQ7tqg7+94ZlMGbRmupEr/UxN0jBkA538Ec93X4Pu3wz+efH1MZJmz//yPnCzrdhuWaujYr3lrze0s7f897DduwFx58v2yj90Nu7bxom8M27FDLzOnhPofO/Qyc6JOfR7eZ7nl33MEIcdS3PLvw04QsiSzzqW2LHlBXXb5v//hIN+MvzKmfDndqaXkhQsaZI6mLHf863NjXhtwGkINcMe/PkerL+CUk0IRpLliQBvgJbNtx5cjzQ2uWqy1Pb59bTgBV9QXfLV6zyt03A0M06sNuDrJk5dNxi/oGJugLW9B/8gGDUImo3sryB/BoGOK5i89l4N1jUXqs5D5S8/NeJ1KToESsDu10t1/QbtrxzObOMhsyqApNSqt0v/UBB1jlg4S5lwNv/Nj3hh4gCU7/4b9qz7D6WtudFrQXhkZd8uPE7TMYl3drXxh9wMs2XXy/TpgFnLHr25l23dnceubL/KZ1pV8/0XDeq+eBer//RcN6y8ZXxAStGUpSHmvpeuJPYfoeSzBoWO/Speff8YptF8Z5fLl88vuTi01yHx673C6TC4WODDyFk/vHeaypnlKJpwlNEGatbbHGLPJW3kAmJeVDqMFJ4DqgXTXZG+ec20sdtz92uVbNqoJ6PMHgTI1gra8Bf0jGzQIqXT3VtA/gkHHFCX67uXUX89iwe5t+eszGuWtvntZePNdFX89g7ZEBe2+DNpdO97ZxKXOpgyaUqPSKv1PTdAxZl5y4EyfgaeBpx9I7/EnBw7i5NqdrfzSvo9X9v6A+p3/zOiq/8Gb57yPPzb1/DEn1+68uuUj/GrfF/mNAt2jv551Nle3OK17QYOQ0ZED3PHtnxQO6r79U1qXQv2cxYGDwMbZpxJ/8TCbv/fCmOe8fuxXbP7eC9z7ieb0+KigeeEgWNA4dDT/a+PnlQvbGLypFpogDZygqcCxHtwAzd1O4vx85itf8Hgp15SpEbTlLegf2aBByHjGFFVyhYWgY4pWfODTJF/7KD/62aNcMPjVrPos5LmmTZz9zqtY6gYVlX49g7ZEBe2+DPp6Bg0ag86mDPpPR6VVOgg/9OrL/PqCT7Lo6c15x5glL+nkjVdfTicHTk8EKMCbCBCUt3Znhjrgp3/rPDzu2p2N9UfhXe3Q92fwjRvA1z3KY3djdm3jlNY/d8pxekYQkiuvWgrnZ9Ar56X4yMcCB46NplN8BA4CAwZ1886YyVef+EXe8xvgq0/8gpuvdFKmBA4aZ5XW2uWVCxo0hm3260QLVZAmMh5Bg6jxdLeNZ0xRpVZYCNry5gUVyy+6gtETm3g2K0XG+7LqEzhoCdidGrQlKmj3ZdDXs9LJYMezjFolTUYQnm453NMPux9OFzVmBgbL4qfv5MkDrzH/5ruCTwQIyh3YX5R3ft9sTfb0wa6T9aduBmCh70/hV2/AezvTwUXBvGqpS9Llhs5sge17i1Zn6MwW52vAlqigQV2ly19yTgOL55zKqyNv5QzsDLBozqlcco6TEjVo0FjrLW8K0qQmBAmixjvbrtSgy68SKyyUk8KifsaMomktggYtQbtTg7ZEBe2+DPp6Bg0aq12l/6nJaDm0o5x26Fnq3xpm9NQG3px/ARhnlmb69Qw4ESCwoOf3B3WpUXj1Z/DL1+Ft82DRO9OzTL1zXnJOA797xk+480SBvHAz/phLzvktABoXLMRJLFCYU47ALVFBg7pKl6+vM3zhw+dzy9fjGDKzpnhdXV/48PnpsWZBg8CgY/CqjYI0qXrjGZgdJKgLOts0qKB/BCudwmI8LV1BulODCtwdHLBlL2jQWO0O7d/HrEVN/KhpU/73a1ETh/bvY/6SZeNOLn3S1ZN8h2XKDurOypkfPa2elJNX7ZwCedX2P0A9nwfqM1qWTJ5lp/wtS0FbooIGdZUuD3DthYu59xPNGRMNcOudPdEgaBAYdAxetXWPKkiTqjfegdmltowFnW0aNHlv0D+ClU5hETRoCdqdGlTQoDRoy9504/08L8eZHfysf3awfZH3DcZgMJbx+xKmFCKhs+8JTnvjZbjqn/P+E3HaP7bCvifgnCvTLUvffuAr/FmO7tE/P34jH/3w76VbloK2RAUN6ipd3guKzpr7Nr7yidU587w945u9GzQIDDpGrtq6RxWkSdWrdB62oOcfT/LeIH8Ew5bB3q+U7tSgxtN9WcmWvWqX/fM8k5MLEvs74bJ/X8KSQiR0jh10vhb5JyJdDri2bjvrTrkHu2IdXDm2e9TUrQac2aOH9u9j+ehL/NkVb+Pvd7zJoTdPhiPz32b4f5tPY/noIIf2/8r5py9gUPf6seN86vLfzDtQ3wKfuvw3ef3YcRpnnxr4/PmCor//4SB//8PB9LYXFAUNAsfbPeoJ+1qiCtKk6lU6D1vQ8493Nl+pLXthy2BfaUGD0kq37FW7cn9fpjqFSOic4f7BL/JPRLpcahS+fztm5bWY331gTPco37gRvv8ncO4Hoa4+o+Xzk9bw9G+cyxARGklyyegL1A9YGDj5+g8deYuz5r6Nzg+cO2aM1jx3jNZZc9/G0JG3aJx9akktS3d+7wXe+NVoumUpSPdldlCUjxcUBQ0ax9M9WsqyTWExvT+tZFqqdEqEcmbzlTLRAKZX91M5QWklWvamm7ClEAkdd9kpHrsb/EEXQCoFj30JfMtOHX7+R8xNvgi/8085u0e58o/gH1udchdcM2Yixtl7f8A5bp63vefExkzEKBR0HTr2q3Tw47VcZS87VSh5rN+1Fy6m9fxFPLT9RW771jPc+bELuX7t2WNWDcgXFOUTNGgczxi5IKZ6DJuCNJl2wpYSIajp1v0UdKC7TKxq/32puLp6eP8XYcsnnVawK//Il1ftS04KD9/ao0/+9Dl+C4p2jz750+f4rQuuOfn6P/cdZ9ms5Isn87zt+6Zzbd+yWUFbrrKDqIuWRgo+LztoedspM9Jfnz9wctHwiVgbtJSgMWj3qJ9/WaufvZzMuazVVI9hU5AmUmWmW/fTeAa6i0yWoSNvMTTnama3fIVFT/0lp/jWHv3VrLN5teUrHJlzNY1u9+JlF50PuyjaPXrZReef3Pfcd5wgcOW1TgucL7kuWz7pBIFuoBa05SqofEHLHz70k4zt8QYtQYPGoGPkPKUuazXVY9gUpIlUmenW/RR0oLvIZDoZtMyijs2ZKTVeO5fUd+uAx9NBy9zz3lNS9+jc897j7nPGsLHy2szyS9c621lj2CotaEtdpQUdgwfBlrUqZwxbKS11xShIE5liQVN2TLfup+l2v1JdxgQtry2Eb26A374PFrw9vTsdtATsHmXfE04XZ5ExbF6Kj0qrdEtdUEHH4AVN2TFepbbUFaMgTWSKjSdlh4iEw5igxZwCdb+AxlNgSZ6WlvM/4gRi37/dCbA8kWUZXZdA8BQfR191HsVUeqWHSRK0ZS9oyo7xKKWl7vKzTy/pXArSRKbYdOu+FBGcQOzcD0L8a/DdP4QP/S00f3Jsl2XQFB/+tUcLufrz8N7O8dY+NIK27AVN2RFUqS11/37LmpLOpyBNZIqpO0+kimW3XB3alfnV47VcZZc/5fSTX1/9+djyAVN8jFl79L/+HR67y0ma+/bfGrP26HRT6ZQdpbbU7fjF4ZLOpyBNRERkvPK1XH1zQ+a213IVtHzQMWxecOdP2QFOoPbzLWNSdkw35aTsKEWpLXCvHSutnII0ERGR8fK3XBXitVwFLX/0VaclrfXP4T//IXMM26wlzv7I2U457zkBUnZMN+NN2QGlzdYstQVuwRmllTPW5oolpRTGmNnAyMjICLNnz57q6oiISK354eZgY8xSo/B3F0PjBbm7R79xIww9B3+wc1JSdoRV9uxLgMU5lrUKWn40Zbki9oOiLXX/fssaGuZGAOZYa4/kKOqUV5A2fgrSRESkooLO1tz7GNz/Ibi5P/dEg5eedlrjbvrupKTsCLPRlC26rBXkn63plfTnVfOXh9wtdd7szjlz5kCRIE3dnSIiImEVNFVG0JQd01h9neGdZ0UAeOdZkbxdnEHzqpWyAP2RI3njsgx1xYuIiIhIVfCn7MglO2WHFBQkr5rftRcu5vGO93Hnxy4E4M6PXcjjHe8LlMgW1JImIiJSO4Km7JhmsheI3zN0LOOrx1sgvpy8aqW01BWjIE1ERKRWBE3ZMc0EXSC+0nnVilGQJiIiUiuCpuzwT0xIjcKrP4Nfvg5vmweL3jk2/1qVC7qMVKXzqhWjIE1ERKRWFFoW6uh+6PtT5/tcyXXr6p1AzePfnqbLSJWTV20iKEgTERGpFeNJrnvK6dD3Z7C81Vk+Kt09ehfs2ua0vr3zusrWO8RKma1ZKQrSREREakXQbsnTF8D2+5zVCfwTDZauhd990BnXtv1/w2W/X5n6VolrL1xM6/mLSsqrNpGUgkNERGS62veEs77nlZ/NnAkKzvaVfwTJfU65aW4iZmsGpSBNRERkulLy21BTkCYiIjJdKfltqIVqTJoxpt23GbHWdpXwnE3ut2uBhLW2I8g5x3NNERGRmqDkt6EWmpY0N1iKWGt7rLU9QMIYEyvynJi1tst9rAeixpjeUs85nmuKiIjUDC/57a6HnUkCLz0Nx486X79xo7P//X85bZPfTrXQBGlAB7DV27DWbgXa8xU2xkSAFverZzPQZoyJlnjOQNcUERGpKf7ktwd+4iS/3XyW8/XATzOT38qkC0V3pxtoRa21iaxDEWNMs7U2nuepUffhHfeeHzXGDBc6p1t2PNcUERGpDUGT3/pXKCikRlYoCLrW50QLRZCGE2jlkiQzCEuz1iaBuXnOkyjhnPnkvaYxZiYw07drVoHziIiIhFvQ5LeFgjq/GlmhIOhanxMtLEFavkWvhgscy2Uj0G+tTRhjWoqcMzmOa3YCXwhQHxERkfAK2uKVHdQd2gXf3AC/fR/M9wUpNdCKBsHX+iy15e00e5xShCVIK5vbhdkCrK7gZTYDX/JtzwJeruD1REREwiNfUDd/JSy5eNKrU2lB1/osteVt47tKW0oqLEHacJ79DQWOZYsBq91u0FLOGfia1trjQDr8Naby2YZFRERCKTUK+3c63+/fCYveMe1ngZba8naaPc5tJZzPWGuLl6owd+LAYWCuL8jCGGNxAq+Cg/iNMd1AzD8JoNg5ccatjfuabtnZwMjIyAizZ88uep8iIiI14bnvwPdvd5aU8kTOdtJ5nP+RqatXlThy5Ahz5swBmGOtPZKvXChScLhBUoIcY8FKCNDa8QVoxpioOzuz4DnLuaaIiMi09dx3YMsnofECuLkfOl9xvjZe4Ox/7jtTXcOaEYogzRUD2rwNN/jq8G1Hs1YHwBjTBkRwUm60uNsdnEzFUfCcJRwXERERT2rUaUFbea2zQsHStTDzDOfr7z7g7P/+nzjlpGyh6O70uEs8JXECr3n+JZ68AMpa2+RuR3C6K8ew1hrf8/Kes5TjReqr7k4REZk+9j4G93/IaTlbunbs8ZeedhLh3vRdOOfKya9flSi1uzMsEwcAKLRuprtsU49vOwkUHblfbC1OrdUpIiJSomMHna+N5+U+7u33yklZwtTdKSIiImF2hjtzcej53Me9/WcUn+EoxSlIExERkdIsu9yZxfnY3ZBKZR5LpeCxL0FkmVNOyqYgTUREREpTV++k2dj1MHzjRmcM2vGjztdv3Ojsf/9fTvt8aRMlVBMHqo0mDoiIyLSUM0/aMidAU560oqpy4oCIiIhUgfM/Aud+EOJfg+/+IXzob6H5k2pBm2Dq7hQREZHg6uphySrn+yWrFKBVgFrSREREpDRHX3UenkO7Mr96vIXYs8vnk2/h9mlOQZqIiIiUZuCr8OO/Grv/mxsyt6/+PLy3M3/5bF55yaCJA2XQxAEREZlWgraM5Wp5++YG+O37YP7KseWnCU0cEBERkYkVNJjKV37+Slhy8YRVq1Zp4oCIiIhUXmoU9u90vt+/U4uwl0BBmoiIiFTWc9+Bv7vYSdcBzte/u9jZL3kpSBMREZHKee47sOWT0HgB3NwPna84XxsvcPYrUMtLEwfKoIkDIiIiBaRGnRazxgvgdx+AOl/bUCrlLCU19Bz8wc5plWet1IkDakkTERGRytj3hLN01JWfzQzQwNm+8o8guc8pJ2MoSBMREZHKOHbQ+dp4Xu7j3n6vnGRQkCYiIiKVccZC5+vQ87mPe/u9cpJBQZqIiIhUxrLLIXI2PHa3MwbNL5WCx74EkWVOORlDyWwnwJEjecf8iYiITG+X3w7f2gi/XA+X/z4seDu89l/wxN/DYD98rBuOvTHVtZxUpcYNmt1ZBmPMmcDLU10PERERqUpnWWtfyXdQQVoZjDEGWAIcLfNUs3CCvbMm4FzVQPdb23S/tU33W9t0v5N77f22QCCm7s4yuC9s3gi4VE6sB8DRQvlSaoXut7bpfmub7re26X4nVdHraeKAiIiISAgpSBMREREJIQVp4XAcuMP9Oh3ofmub7re26X5rm+43RDRxQERERCSE1JImIiIiEkIK0kQmkTEmMtV1EBGR3ML2Ga0gLQTC9kMxUYwxLcaY3jzH2n2PTZNdt0rId7/ufmuMscBhY8ygMSY6BVWUSVCrv8/iqMX31xjTl2NfzX1Ge7LvN8yf0cqTNkWMMS1An287AbRaaxNTV6uJYYxpBq4HIsCYH3RjTDsQsdZ2udttxpiYtbZjUis6QYrdr7t/tft9shbeYwDfB/daIJH9/rnvsyf9flerQvdba7/PbiBynbvZhPMz3GGtTfrK1Mz7W+x+a+399TPGtAEtWftq6jPaL9f9EuLPaE0cmCLuD4r3gxCqH4qJ4t5jp7V2ddb+QbI+4Iwxh621cye7jhOpwP22Af3+P3DVLvsD22tBtNaud7fHfMgDa6v1Q76E+62p32djTDfQba2N+7aj1tpWd7vW3t9i91tT76/HDU7bgZi11vj21+pndITc9xvaz2h1d06thLU2Xiu/8KVwf0miOe454rZISci572FLVrfPZqDN10XQAWz1Dlprt+J8OFadEu8Xauv3OUpma8Ng1nbNvL+uYvcLtfX+eq4Devw7avwzesz9hp26O2Wy5evnT7rH4pNXlUl1nTFm2P2+alscfKJkvl/eB3rUvc+8H/Jea0WVyXu/vu9rhteC5NME9EPxP+LV+P4Wut9a5QZcAzkO1eRndIH79YTyM1pB2tQK5Q9FhTXk2T9c4Fi1S+D+Fw5gjGkwxnRbazdOcb3Gxe0SyO728D7YE9TYh3wJ9+upyd9nt7WwBfACmZp6f7PluF9Prb2/a6y1PTkmQtTqZ3S++4UQf0aru3PqJIABa+1Wt6tg0B0HITXG7SLx/+HqB9prbJbYRpwxHQlq90Pez3+/UKO/z+7Ys15go+9ea/b9zXO/UGPvrzGmzVpbVd1+5Sh2v2H+jFaQNkXC/ENRYcN59jcUOFZTfB/+oZjiXS63G6EFWD/VdZkMue63Vn+frbU97kSYjlpLw5BLvvutpffXrXOyQJGa+owu4X7HCNNntLo7Q8JamzDGQA10FxSRAOcXJ2smTYQaHNvjfkDsBVZ7v/jV+MFeRAzn/pLudk19yOeQfb9j1ODvcwzoM8ZspfbfX/Ddb66JAlX+/l4HNPkmATRBOsVMAt/Ywxr5jC7lfkP7Ga0gbQpMkz/cOVlrk26OoQay/rupxgHHJRrI+qCPQm3cr9vlszHrw7xmA/Fc91trv89u3e8DNvju03vfWoAtXrlaeH+L3a8xZgs19P5md/u5Y/Da/Xnuaukzutj9uu9laD+j1d05dUL7QzGB8o1PiQFt3oY7DqTaB+FCjvt1P/Szs3l3UgP3675vMd8frqg7uy+J80cu1+tRtT/f+e7XPVxLv8/ewHn/+xdxvyZq8P0teL/u11p6f7NFcuyr1c9oyLrfsH9GqyVtCritSaH9oSiXLwN/G05Khm5gh/cfjTvDZpOXEBOYV80zpUq43y7f+JYmoK/aB+26yR8jOPcbdb9v5eTPsPch7/23WtUf8oXut9Z+n621cWNMT1ZQcj0Qt9Z6aSlq5v0t5X5r6f31c983LyFzL05C3/5a+4z2FLjf0H5Ga8WBKZT1Q7EjLD8UIoW43QOHcx3LyuK9Cae7JEIVf8gHvF+ogd9n9547fbsijF0WqibeXwh0v1AD769UDwVpIiIiIiGkMWkiIiIiIaQgTURERCSEFKSJiIiIhJCCNBEREZEQUpAmIiIiEkIK0kRERERCSEGaiMgkc1crmPLFm0Uk3BSkiYhMvqpN/Coik0dBmojI5ItmLUMkIjKGgjQRkUnkrvVaCwtzi0iFKUgTEZlcG4Huqa6EiISfgjQRkcmlrk4RKYmCNBGRSaKuThEJYsZUV0BEJIzcgCrqPgC2Am3ecWtt1zhOuxGI5bleG7AWeB1IANcDm621CupEpikFaSIiWdwcZlFr7VZ3+zDQZK3daIzpBtYA4wnScnZ1GmPagfXW2lbfdhuwYbz3ICLVT0GaiMhYLdbaHt92BOhzvx9XjrN8XZ1uQNgNzPXtTgBJa21yPNcSkdqgIE1EZKwt3je+lQH6AbIDJ/d4G05gFQW25pkYkK+rs9t9jv+8zd71RGT6UpAmIpIlK2BqAeIFWrV6rbWrAYwxEeARYHWOcvlmdbbgBHB+rZxsuRORaUqzO0VECmsFBnIdcLswI962G8hFstflNMa0kL+rkxznb0EtaSLTnoI0EZEsWUFWC7DDd6zNd2wNkMx6ehKnu9JvPYUT2KZb2NyADmtt3BjTrIXYRaYvBWkiIj5uEDZojIm43w+7D687s8FXPOIdy9KQtZ2zq9Pd541l886/kZOBX4sS34pMXxqTJiKSKQ70ANfhBGCtQIcxpgEga9ZnkrEBGfgCt3xdnT7rgY3GmEGcGZ3rjTG9xphNRZ4nIjXOWGunug4iIlXJHZPWa61t8u0bBFq9FjA3r1pMLWIiEpS6O0VExil7NQC3uzKZFZBprU4RGRd1d4qIlGe9MSYGbMdZ1mm9d8Dt6lQqDREZF3V3iohUiLo6RaQc6u4UEamcBgVoIjJeakkTERERCSG1pImIiIiEkII0ERERkRBSkCYiIiISQgrSREREREJIQZqIiIhICClIExEREQkhBWkiIiIiIaQgTURERCSE/n8b3ZhFxoZ3AQAAAABJRU5ErkJggg==\n", "text/plain": [ "
" ] }, "metadata": { "needs_background": "light" }, "output_type": "display_data" } ], "source": [ "periodic_m_eff.show([4,47], comp=m_eff, ylabel=r'$am_\\mathrm{eff}$')" ] }, { "cell_type": "markdown", "id": "e3762e68", "metadata": {}, "source": [ "Arithmetic operations and mathematical functions are also overloaded for the `Corr` class. We can compute the difference between the two variants of the effective mass as follows." ] }, { "cell_type": "code", "execution_count": 13, "id": "e56d164c", "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAlIAAAGLCAYAAAAMBnwJAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjMuNCwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8QVMy6AAAACXBIWXMAAA9hAAAPYQGoP6dpAAAtoElEQVR4nO3db2xd1Znv8d+TtInJOD4ndhrikqRgTwtj0ltqEhhQ0VxN44hIjYpGNiDBXKka1SlSpaqaKpmMGEURFUxQZ6RqqunEvCxctSRzRRUkuARmbkU1lBI8lSZEMG0MJaXODYl77OSCk47z3Bdnb7N9fP7ss33+bZ/vR7LqvdfyPis98cnDWs96lrm7AAAAUL0VzR4AAABAWhFIAQAAJEQgBQAAkBCBFAAAQEIEUgAAAAkRSAEAACREIAUAAJDQx5o9gHozM5P0SUkXmz0WAACQKmsl/dbLFN1c9oGU8kHUb5o9CAAAkEqbJL1XqrEdAqmLknTmzBl1dXU1eywAACAFZmZmtHnzZqnCilY7BFKSpK6uLgIpAABQUySbAwAAJEQgBQAAkBCBFAAAQEIEUgAAAAkRSAEAACREIAUAAJAQgRQAAEBCBFIAAAAJEUgBAAAkRCAFAACQEIEUAABAQgRSAAAACaXi0GIzGw6+7ZY04e4vNnM8AAAAUgoCKTPLShpy9z3B9XFJBFIAAKDpGhpImdkOSXvcfaRI22jkMuvujwff3yspF2nLmdkOZqUAAECzNSSQMrNBSfdJykrqK9I+qkjwZGbDZnbI3fdJ6pd0IdJ9KngOAABAUzUk2dzdx4Og6HiJLvskHY30PypptERfKZ8rBQAA0FRN37UX5ED1uftEQVM2mMk6Lakncr9bUmFfAACAhmt6IKUiS32BXND2dEGfLPlRAACgFbTCrr1Sy3RTkrrdPWdmPwpKIHRLOlzuYWa2WtLqyK21tRkmAADAQq0QSFUU5EzFtV/SgXqNBQAAINQKS3tTJe53l2kr5zFJmcjXpoTjAgAAKKsVZqQmpHzSubvnIvezSpBU7u6XJV0Or81sicMDAAAorukzUkHwNKEiuVLuPt7wAQEAAMTU6ECqVGL5IUnheXphgc59DRkRAABAQo2ubD4sqc/MDkt63d3HJMndx8xsb1jhXFJPUMATAACgZZm7N3sMdWVmXZKmp6en1dXV1ezhAACAFJiZmVEmk5GkjLvPlOrX9BwpAACAtCKQAgAASKgVyh8AAAA03NxV18/fntK5i7PasLZDt93QrZUrqiubRCAFAADazvMnJ3Xw2ClNTs/O3+vNdOjA7gHdvbU39nNINgcAAMtGnFmm509O6qEnx1UYAYW9vv/goO7c8gexks2ZkQIAAMtCnFmmydyHeviZk4uCKEnz9x5+5qT+5//4b7Fek2RzAADQ0uauul45fUE//sV7euX0Bc1dXRwGhbNM0SBKks5Oz+qhJ8f1/MlJSdIjz57S+UtXyr7e+UtX9Pj/fjPW2JiRAgAALSvOLNPcVdfBY6dKzjKZpIPHTmloYKM+/P1crNeN248ZKQAA0JLizjIdOfHuoj5RLmlyelZHTryre7dtjvXaf/b5TbH6EUgBAICGq7RcV2mWScrPMs1ddT3x8tuxXvOJl9/Wzps3qjfToVJFDkz5Ga8//aNrYz2TpT0AANBQcZbrqpll+updN+iv/tfJiq/71btu0MoVpgO7B/TQk+MyaUGgFgZXB3YPxK4nxYwUAABomLjLddXMMo1s2xJrlmlk2xZJ0t1be/X9Bwe1MdOxoN/GTIe+/+BgVXWkmJECAAA1UamGUzVJ4fWeZbp7a6+GBjZS2RwAADRfrZfrRrZt0Xdf+pXOTs8WDbxM+RmkwlmmwjFsLFOtfOUK0x39PYn+vPPjoLI5AABYijiVwu/e2qsv/t3/0en3/1/F5/V/4g/00l/+9/nnSsVnmYotw9Xi/DxJmpmZiVXZnBwpAABQUi131331rhtivWbYL0kuUzjL9OVbrtMd/T2JgqhqsLQHAACKavZynVS7XKZ6YUYKAAAsUo/ddWFSuKRFO+zKlR5o9CxTNQikAADAAq22XNfKWNoDAKDNVErIfuGNs7GW61544+yyXK6rBoEUAABtJE7e09MnzsR61tMnzmjXZ3sTVQqvRemBVsDSHgAAbSJu3tM1H18Z63lhv+W2XFeNVNSRMrPR4NtbJR1396NV/Cx1pAAAy16cquJfOPQvJZfswiW4n+77U52bmdXu7/1U5y9dKfl66ztX6djXv6De7DWxx5AmcetItfzSnpkNSpoKgyczczNb5+655o4MAIDWUOsyBfff9il9+56tZYthfvuerQuCKGn5LNdVo2ZLe2a2w8yOlGgbjXztrfLRfZKGItcTwT0AANpePcoUSO29XFeNJS/tBTNG90nKStrm7rcWtI9Kyrr748H1sKTt7r6vitfIunvOzLKS3nb3dVX8LEt7AIBlqZrluiMn3o11CPDf/tlW3X/bpxa8xnJZrqtGw5b23H1c0ngQIG0r0mWfIjNK7n7UzJ4I7sd9jVzw7ROSRpKPFgCA9KgUxNS7qrjUnst11ahrjlQwg9Tn7hMFTVkzG3T38XJLfeEsVvCsvZIOu/uL9RktAACtI07eUzXLdfff9qlEZQpQXr3LH5TKZcqFbe7+eKmvsHMw2zXu7i+a2WCwnAgAwLIUN++p3auKt4J679rrLnF/qkzbAkHQdERSzsykfL4V4TIAIJUqLddN5j7Uw8+cLHs8y8PPnNTnNmXbvqp4K2j58gdBDlbsd9fMVktaHbm1tuaDAgAggTjLdY88e6ps/SZJOn/pih559pT+8cFb27qqeCuo99LeVIn73WXalmq/pOnI12/q9DoAAMQWd7nuw9/PxXpe2I/luuaq94zUhPRR+YLI/WzYVgePSfr7yPVaEUwBAJpo7qrr4LFTJZfrTNLBY6c0NLBR927brH996/2Kz7x32+b571mua566BlJB7acJ5WegcgVt43V6zcuSLofXQV4VAAB1Uynv6YU3zsYqU/DCG2e18+aN6s10VMx72nnzxgX3Wa5rjloGUqWSxw9JGpYUFuQcVRU1pAAAaGVx8p6ePnEm1rOePnFGuz7bS5mCFFlyjlRQjuCQ8sHRoJkdjhwyLHcfC/qFx8P0R0sbAACQVnHznq75+MpYzwv7kfeUHks+IqbVcUQMAKAeJnMfavf3flp2h936zlU69vUvSFLsvtGDgNv1eJZWEPeImHrv2gMAIJXmrrpeOX1BP/7Fe3rl9AXNXV048fCdF96KVabgOy+8pd7sNfr2PVtlWlzPJ7z37Xu2LgiipI/ynr58y3W6o7+HIKoFtXwdKQAAGi1O3lP/JzpjPSvsFy7XFT53Y8FzkS4EUgAARIR5T4WJL2HeU5ijdPr9S7GeF+1HmYLlh0AKAIBANcezfGvnjfrJf75fMe/pWztvXHCPMgXLCzlSAIC20Qp5T1hemJECALQF8p5QDwRSAIBlj7wn1AuBFABgWSPvCfVEjhQAINUq5T098uypWHlPjzx7irwnVI0ZKQBAasXJe/rw93OxnhX2I+8J1eCIGABAKpXKewpnksK8p+f+Y1IPPTVe8Xnff2BQuz77UZDE8SztLe4RMcxIAQBSp5q8p503b1RvpkNnp2eL9jflZ5t23rxxwX3ynhAHOVIAgJZTy7ynlStMB3YPSCqe9yRJB3YPMNuERJiRAgC0FPKekCbkSAEAWgZ5T2gV5EgBAFpKpSCGvCekEYEUAKDu4izXVXPO3d/de4sO7B7QQ0+Oy6QFwRR5T2gkks0BAHUVLtdFgyjpo+NZnj85KSn5OXcbMx0L2jdmOuaXAIF6Y0YKAJBYLZfrOOcOaUQgBQBIpNbLdZxzhzRiaQ8AULV6LNdxzh3SiEAKAFCVuauug8dOlV2uO3jslOauetXLdeQ9IW1Y2gMALFAp7+mFN84umomKckmT07N64Y2ziZbryHtCmqQqkDKzYUly96PNHgsALEdx8p6ePnEm1rOePnFGuz7bq2/fs1UPPZkvnlmsTEGx5TrynpAWqVnaM7OspD3NHgcALFdx856u+fjKWM8L+7Fch+WsZkfEmNkOSXvcfaRI22jkMuvujyd4/qikrKSJamakOCIGAOKVKdj9vZ9WXII79vUvSFLsvtGZJo5nQZo07IgYMxuUdJ/yQU5fkfZRRYInMxs2s0Puvq+K1xiW9LSk0Up9AQALxVmue+TZU7HKFDzy7Cn944O3slwHBJa8tOfu40FQdLxEl32Sjkb6H1UVAZGZ9Sk/C5VbyjgBoB3FXa778PdzsZ4X9mO5Dsira7J5kNfU5+4TBU1ZMxt093Ez21vq54NZrMHgWX2StkvqN7PxIs8EAERUKlNgypcpGBrYqHu3bda/vvV+xWfeu23z/PfsrgPqv2tv0VJfIBe0jVfKl4rmQ5nZdkmvEUQBQGXVlCnYefNG9WY6dHZ6tmjgZcrPNu28eeOC+yzXod3Ve9ded4n7U2XaigqS2XdI2hPMTgEAyqimTMHKFaYDuwckFa8qLkkHdg8w2wQUSE35A3d/0d1vdfehcjNSZrbazLrCL0lrGzhMAGiYuauuV05f0I9/8Z5eOX1Bc1cXziVRpgCov3ov7U2VuN9dpm2p9ks6UKdnA0BLiLMT72++NKCfvzNVsUzB33xpYP6avCegOvWekZqQ5pPOo7JhWx08JikT+dpUp9cBgKaIuxMv6SHAYd7Tl2+5Tnf09xBEAWXUNZAKShZMqEg+lLuP1+k1L7v7TPgl6WI9XgcAmmEy96EefuZkyZ14LunhZ05qMvehJJbrgHqr5dJeqeTxQ5KGJYUFOUeVry0FAChQqfp3tYUzJZbrgHqqZWXzYUl9ZnZY0uvuPiZJ7j5mZnsjR7z0VFPVHADaRZy8p2oLZ4YoUwDUR83O2mtVnLUHoNninDEX5j0VfiKHvcJluOf+Y1IPPVU5M+L7Dwxq12dZtgOSathZewCA0uLMMlXKe5LyeU+f25RNXDgTQH2kpo4UAKRN3N1133nhrVh5T9954S0KZwIthkAKAOqgmt11/Z/ojPXMsB878YDWwdIeACRQKe+pmlmmuE6/f2n+e3biAa2BQAoAqhQn7+mWzVn98/h7FZ91y+asdvzRtfrJf75fsQL5t3beuOAeO/GA5mNpDwCqEDfvaX3n6ljPW9+5OnEFcgDNRyAFAIFKhwBXk/d0arLkbukFwn7kPQHpxNIeACjecl01eU/77r5Ja1d/TI8+92bJvn+96ybd8/nr5q/JewLSh0AKQNsrVQwzXK4LZ4Sq2V23oatDo3/Sry09ayoGaFHkPQHpQiAFYFmrtLuummKY0V1z5bC7DmgfBFIAlq1aL9d9a+eN7K4DsADJ5gCWpbi766pZrmN3HYBCBFIAUqeWu+uqXa5jdx2AKJb2AKRKKyzXkfcEIMSMFIDUaKXlujDv6cu3XKc7+nsIooA2RSAFIBVYrgPQiljaA5AKLNcBaEUEUgBaQqV6T0mW6x56clySFsxihU8st1wHAHERSAFoujgJ5EmX6wqfu7FMVXEAqJa5F8s4WD7MrEvS9PT0tLq6upo9HAAFSh3PEs4chTlKk7kPtft7P624XHfs619YMNNUaaYLAIqZmZlRJpORpIy7lzyFnBkpAHVTy+NZWK4D0IoIpADURa3rPf3dvbewXAeg5RBIAai5Ust1Yb2ncLmumgTyELvrALSSVARSZpaVNCppQlKfuz/e3BEBKKWa5bpqE8hDLNcBaBVpKch5RNKYux+VdF+zBwO0s0rn3FVb72l956qyfQvrPQFAK6nZjJSZ7ZC0x91HirSNRi6z1cwomdmgpJy75yTJ3W9d6lgBJBMn76kR9Z4AoFUseUbKzAbN7JCkEUl9RdpHlQ+extx9TNJE0D+ubcFzdpjZaEFQBqBB4p5zx/EsANpJzepImdmwpP2FM0ZmdlrSkLtPRO79zt3XxXzuXkn3hc81s9cljUSfV+HnqSMFLFE1NZwkUe8JQOq1RB2pIEm8r0jQkzWzQXcfDwKlooIlwAlJJyK3pyQNBvcB1EClIKbaMgXUewLQLuq9a2/RUl8gF7SNx8iXelELE8y7JY0vfWgApNrnPUkczwKgfdQ7kOoucX+qTNsC7p4zsx9FcqMei7usB7SzOEtlces9JSlTQL0nAO0gFXWkgrIHsZjZakmrI7fW1n5EQGuLM8s0d9V18NipkvWeTNLBY6c0NLBR39p5o37yn+9XzHsqLFPAch2A5a7edaSmStzvLtO2VPslTUe+flOn1wFaUtzddS+8cXZRnyiXNDk9qxfeODtfpsD0UZ5TKLxHmQIA7ajegdSENJ90HpVV/ZLFH5OUiXxtqtPrAC2nUlVxV76q+GTuQz194kysZ4b9KFMAAIvVdWkvyG+aUH4GKlfQVpeEcXe/LOlyeG1GPgaWj1rurrvm4ytjvWa0H3lPALBQLQOpUsnjhyQNS3pcmi/Qua+Grwu0hTh5T7dszuqfx9+r+KxbNme144+u1c/fmaqY9/Q3XxpYcI+8JwD4SC0rm++TNGhmh6PVx4Nq5gqqku+V1M+hw0B852ZmNfaT0/pakbynyelZfe3JcY395LTOzcxqfefqEk9ZaH3navKeAKAGljwjFSzRjavMLBOBE5DcD372a/3Dv/yqbJ9Hn3tTFy//V+xnnpqc0a7P9lLvCQCWKBXlD4B2NtAb72ijgd4u3fqpdVq7+mN69Lk3S/b761036Z7PXzd/Td4TACRHIAW0uCtzV2P329DVodE/6deWnjUV86miyHsCgGQIpIAmq7QTb8PajjI//ZFoP2aZAKAxCKSAJjg3M6tzFy/r3351XmMvTyzYObe+c5VG7+rTnX+4XhvWrtZtN3SrN9Ohs9OzRetDmfI5TbfdsHDjLLNMAFB/BFJAEzz16rv67ku/LNp2/tKV+Rynb3zx03rg9i36yp3Xl8x7cklfufN6Xbh0WRu64s1eAQBqo96VzYG2NHfV9crpC/rxL97TK6cvaO7qwrmk+7dv1vrOVWWfsb5zle7fvllPvfpu2eRxKb9r76lX313yuAEA1WFGCqixOIUz37nwQawK5O9c+EAP3L5FQwPXSsoHaG/8dlq/++D3Wrfm47r5k5n5vKcNa+PVkAIA1A6BFFAj52Zm9cy/v1d09igsnBmWHjh3sfRhwQueeXFWd/T3LFiy+9zmbK2GDABYIgIpIKZKu+uqKZx5Z//6WK8Zd8ceAKA5CKSAGOIs11VTOPP6njVa37mq4jl31/esWdrAAQB1RbI5UEY159xVUzjzh6+diZUj9cPXziQeOwCg/piRAsqoZrku7ozUqpUr5hPI49SRAgC0LgIptLVKeU/VLNedmpyJ1Tc8MHhDV4e2XpfRX9zVRwVyAEgpAim0rTh5T9Us1/35H3+q6gODJSqQA0CakSOFtlNN3tOqlfF+RVatXDF/YPA/PTio3szC3Xa9mQ7904ODGv2TfqqPA8AywowU2k41eU9xhct1EgcGA0A7IZBC26km7+nWT61juQ4AUBKBFNpONXlP4XLdlp41FfOpAADth0AKy06lnXjV5D2FWK4DABRDIIVlJc5OvGrLFIRYrgMAFDJ3b/YY6srMuiRNT09Pq6srXm4M0qfcgcGhaC5T3L7ssAOA9jQzM6NMJiNJGXcv+V/gzEhhWahmJ95f7ryRvCcAQE0QSGFZqGYnXoi8JwDAUqUikDKzHZKy4bW7H23eaNAMlRLIq9mJF0XeEwBgKVo+kDKzrKQ+dx8Lrg9LIpBqA+dmZnXu4uVYB/tuWBsvlyluPwAA4qhZIBXMGu1x95EibaORy6y7Px73ue6eM7N9ZnbC3ccl5ZY+WqTBU6++q+++9MuibecvXZlPFv/GFz+t+7dv1vrOVQuCrULrO1fp+p41dRkrAKA9LfmsPTMbNLNDkkYk9RVpH1U+eBoLZpUmgv7V2CfpdTN7XdJjSx0z0iEMjspZ37lK92/frB++dqZsECXlg68fvnamlkMEALS5Jc9IBbNE42Y2LGlbkS77JA1F+h81syeC+3ENSdoT/MzrkvqTjxhp8c6FD2IFR+9c+EAP3L5FQwPXxloGBACgVuqaIxXJb5ooaMqa2aC7j5vZ3lI/7+6PBwHa8SDBfMzMjpjZDnd/sY5DRws4d3G2cqeg3x39PdrQ1aGt12X0F3f1sRMPANAQ9U42X7TUF8gFbeMx8qW6JZ2IXP9I0tTSh4ZWlzSBnJ14AIBGqXcg1V3i/lSZtgXcfczM9ppZuGw4FSwnIuUqlTS47YZu9WY6dHZ6VsXq75ukjZn8zwEA0AwtX/5Ayi/xxe1rZqslRRNh1tZ+REiqqpIGXfkq4w89OS6TFgRTYbh1YPcAy3YAgKapdyBVagmuu0zbUu2XdKBOz8YSVVPS4IHbt2jTujXav+umRUFXTxB0bVq3RudmZjkTDwDQFPUOpCakfNK5u+ci97NhWx08JunvI9drJf2mTq+FKt2/fbOeevXXFes95fvFD7q+OfSZuowXAIBy6hpIBcU0J5SfgcoVtNUlz8ndL0u6HF6bsezTSpKUNKiEkgYAgGapZSBVKuP3kKRhSY9L8wU6q6khhWUkSUkDAABa1ZIDKTMblHSf8sFSX3AW3uvh2XiRXXejyi/p9bg7gdQyVWknHmfiAQCWk5pVNleZWaZqdt0hfarZiUdJAwDAcpKK8gdobdXuxPvKndfP3yvkkr5y5/W6cOkyy3oAgJa35EOLgWoOF37q1XdLBlGhR597U0+9+m4thwgAQF0wI4UlS7oTb+6q643fTut3H/xe69Z8XDd/MjOfT8VOPABAGhBIYcmWshPvc5uzdRoVAAD1x9IeloydeACAdsWMFCricGEAAIojkEJZz5+c1MFjpzQ5/dHyXW8mf5jw3Vt7JUkXLl1mJx4AoC2Ze7E5hOXDzLokTU9PT6urq6vZw0mNczOzeubf3yu7w+6vd92kez5/XdnyB1GciQcASIuZmRllMhlJyrj7TKl+zEihqB/87Nf6h3/5Vdk+jz73pi5e/i/9+R9/ip14AIC2RCCFogZ6483eDfR2aUNXBzvxAABtiV17KOrK3NWa9gMAYDkikEJRlDQAAKAyAikUFZY0sBLtpvzuPUoaAADaGYEUilq5wnRg94AkLQqmwusDuwcW1JMCAKDdEEihqHMzs9q0bo3277pJPQUHEvd0rtL+XTdp07o1OjcT73gYAACWI3btoahytaHOX7oyX1+K2lAAgHZGIIWiHrh9y3xtqHKoDQUAaGcEUiiqsDYUAABYjBwpAACAhJiRalNzV10/f3tK5y7OasPafBkDduABAFAdAqk29PzJSR08dkqT0x/tuOvNdOjA7gHdvbW3iSMDACBdWNprI+dmZjX2k9P62pPjC4IoSZqcntXXnhzX2E9OU9IAAICYCKTayA9+9uv5sgWlPPrcm/rBz37doBEBAJBuLbe0Z2Y7JA26++ORe8PBt92SJtz9xaYMLuUGertq2g8AgHbXUoGUmY1KGpL0WuReVtKQu+8Jro9LIpBK4Mrc1Zr2AwCg3VUdSAUzRnvcfaRI22jkMhudVYrD3cfMTJKykdv3SspFrnNmtoNZqeptWBuvLlTcfgAAtLvYOVJmNmhmhySNSOor0j6qfPA05u5jkiaC/kvVL+lC5HpKCwMtxHTbDd3qzXQsOoQ4ZMrv3rvthu5GDgsAgNSKHUi5+7i775N0vESXfZKORvoflTRaou9S8S99AitXmA7sHpCkRcFUeH1g9wD1pAAAiKkmOVJBHlOfu08UNGXNbNDdx81sb6mfr7AEeFr5WalQt6TC10EM52ZmtWndGu3fdZPGXp7Q+UtX5tt6Oldp9K4+bVq3RudmZjkeBgCAGGqVbL5oqS+QC9rGq82Xinha0hOR6yz5Uck89eq7+u5Lvyzadv7SlfnSCN/44qf1zaHPNHJoAACkUq0CqVJLbVNl2hYJyhyMBN9PuPtRd8+Z2Y+Ctm5Jh5c82jb1wO1bNDRwbcV+G9aubsBoAABIv5YqfxDkVR0tcT8WM1stKRoJrK3B0JaFDV0dLNkBAFBDtapsPlXifneZtnrZL2k68vWbBr8+AABoE7UKpCak+aTzqKwanxj+mKRM5GtTg18fAAC0iZos7QV5TBPKz0DlCtrGa/EaVYzlsqTL4XVQ4BMAAKDmksxIlUoePyQpPBMvLNC5L8mgAAAA0iD2jJSZDUq6T/lgqc/MDkt6PahiHh7vsjescC6pJyjgCQAAsCyZuzd7DHVlZl2Spqenp9XV1dXs4QAAgBSYmZlRJpORpIy7z5TqV6tkcwAAgLZDIAUAAJAQgRQAAEBCBFIAAAAJEUgBAAAkRCAFAACQEIEUAABAQgRSAAAACdXkrD20hrmrrp+/PaVzF2e1YW2HbruhWytXcNYgAAD1QiCVcudmZnXu4mX926/Oa+zlCZ2/dGW+bX3nKo3e1ac7/3C9NqxdrQ1dHU0cKQAAyw+BVMo99eq7+u5Lvyzadv7SFT363JuSpG988dP65tBnGjk0AACWPXKkUu7+7Zu1vnNV2T7rO1fp/u2bGzQiAADaB4FUyr1z4YMFy3nFnL90Re9c+KBBIwIAoH0QSKXcuYuzNe0HAADiI5BKuQ1r4yWQx+0HAADiI5BKudtu6FZvpkOlihyYpN5MvhQCAACoLQKplFu5wnRg94AkLQqmwusDuweoJwUAQB0QSKXcuZlZbVq3Rvt33aSegt17PZ2rtH/XTdq0bo3OzZAjBQBArVFHKuWoIwUAQPMQSKXcA7dv0dDAtRX7bVi7ugGjAQCgvRBIpdyGrg6OfgEAoEnIkQIAAEiIQAoAACChllvaM7Mdkgbd/fHIvdHg21slHXf3o00ZHAAAQERLBVJBwDQk6bXIvUFJU2HwZGZuZuvcPdecUQIAAORVHUgFM0Z73H2kSNto5DIbnVWKw93HzEySspHbfcoHV+Es1ERwb7yaZwMAANRa7EAqmBm6T/kgp69I+6giwZOZDZvZIXfft5QBuvtRM3sxeGZWUre7E0QBAICmi51s7u7jQVB0vESXffpo1kjBUtxoib5ViSzjPSFp0UwYAABAM9QkRyqYKepz94mCpqyZDbr7uJntLfXzcZYAg58/7O4vLm20AAAAtVGrZPNFS32BXNA2Xm2+VJSZDQfPeDFYYhTLewAAoNlqFUh1l7g/VaZtkSBgGgm+nwjyowYlHZGUCxPR3d2WOF4AAIAla6nyB0Fe1dGCe+OSYgdOZrZaUvRgubW1GR0AAMBCtapsPlXifneZtnrZL2k68vWbBr8+AABoE7UKpCak+aTzqGzY1kCPScpEvjY1+PUBAECbqMnSnrvnzGxC+RmoXEFbQ5PC3f2ypMvhdZBXBQAAUHNJZqRKJY8fkjQcXgQFOpdUjBMAAKCVJalsPiypz8wOS3rd3cek+eNd9oYVziX1LLWqOQAAQCszd2/2GOrKzLokTU9PT6urq6vZwwEAACkwMzOjTCYjSRl3nynVr1bJ5gAAAG2HQAoAACAhAikAAICECKQAAAASIpACAABIiEAKAAAgIQIpAACAhAikAAAAEiKQAgAASIhACgAAICECKQAAgIQIpAAAABIikAIAAEiIQAoAACAhAikAAICECKQAAAASIpACAABIiEAKAAAgIQIpAACAhAikAAAAEiKQAgAASOhjzR4Aypu76vr521M6d3FWG9Z26LYburVyhTV7WAAAQARSLenczKzOXbysf/vVeY29PKHzl67Mt63vXKXRu/p05x+u14a1q7Whq6OJIwUAoL21XCBlZjskDbr740XahiXJ3Y82fGAN9NSr7+q7L/2yaNv5S1f06HNvSpK+8cVP65tDn2nk0AAAQERL5UiZ2aikPSXasqXalpv7t2/W+s5VZfus71yl+7dvbtCIAABAMVXPSAUzRnvcfaRI22jkMltsVqkcdx8zM0nKFmm+V9Lxap6XVu9c+GDBcl4x5y9d0TsXPlBv9poGjQoAABSKPSNlZoNmdkjSiKS+Iu2jygdPY+4+Jmki6L9kwZLe07V4Vhqcuzhb034AAKA+YgdS7j7u7vtUelZon6Sjkf5HJY2W6BubmfVJmnD33FKflRYb1sZLII/bDwAA1EdNks2D/KU+d58oaMqa2aC7j5vZ3lI/X2EJcDB4jT5J2yX1m9l4kddaNq7vWaP1navKLu+t71yl63vWNHBUAACgUK127S1a6gvkgrbxavOlQtEdema2XdJryzmIkqQfvnYmVo7UD187w649AACaqFaBVHeJ+1Nl2hYJcqFGgu8nCoKoHZJ2SBpc7jNSD9y+RUMD18aqIwUAAJqnpepIBYFT0RpR7v6ipFsrPcPMVkuKRhhrazO6xtnQ1aENXR3ael1Gf3FXH5XNAQBoUbUKpKZK3O8u01Yv+yUdaPBr1s3KFaY7+nuaPQwAAFBErQpyTkjzSedR2bCtgR6TlIl8bWrw6wMAgDZRkxkpd8+Z2YTyM1C5grbxWrxGFWO5LOlyeB0U+AQAAKi5JDNSpZLHD0kaDi+CAp37kgwKAAAgDWLPSJnZoKT7lA+W+szssKTXgyrm4fEue8MK55J6ggKeAAAAy5K5e7PHUFdm1iVpenp6Wl1dXc0eDgAASIGZmRllMhlJyrj7TKl+LVX+oBHmrjrlBAAAQE20VSD1/MlJHTx2SpPTHx3225vp0IHdA7p7a28TRwYAANKoVuUPWt7xU2f10JPjC4IoSTo7PauHnhzX8ycnmzQyAACQVm0TSP3tc2+qWDZYeO/gsVOau7q888UAAEBttU0g9X9nLpdsc0mT07M6cuLdxg0IAACkXtsEUnE88fLbzR4CAABIEQKpiK/edUOzhwAAAFKkbQKpa7tWq1SRA1N+997Iti2NHBIAAEi5tgmk/mrXTZK0KJgKrw/sHqCeFAAAqErbBFJDAxv1/QcHtTHTseD+xkyHvv/gYNE6UnNXXa+cvqAf/+I9vXL6Arv6AADAAm13REzcyuYU7wQAoH3FPSKm7QKpOJ4/OamHnhxfVHcqDLeKzWBVc/QMx9QAANDaOGsvocnch3r4mZNli3c+/MxJfW5TVr3ZayRVN3vFTBcAAMtH2+RIxfXIs6d0/tKVsn3OX7qiR549Jemj2as4R89U0xcAALQ+AqkCH/5+Lna/uauug8dOxTp6ptJMlys/0zWZ+zDZwAEAQMMRSBW4d9vm2P2OnHh30exSVPTome+88Fasma7vvPBWNcMFAABNRCBVYOfNG9Wb6ahYvHPnzRtjHynzxMtvq/8TnbH6xu0HAACaj0CqwMoVpgO7ByRVLt4Z90iZr951g06/fylW37j9AABA8xFIFXH31t5YxTtHtm2JNXs1sm2LvrXzRq3vXFX2ddd3rtK3dt5Ygz8BAABoBOpIlRGn3lO4E0/SgkTyYjWnqukLAACah4KcgaUEUnFRRwoAgOWFQCrQiEBKorI5AADLCZXNG2zlCtMd/T017wsAAFoXyeYAAAAJEUgBAAAkRCAFAACQEIEUAABAQgRSAAAACRFIAQAAJEQgBQAAkBCBFAAAQEJtU5BzZqZkUVIAAIAF4sYN7XBEzHWSftPscQAAgFTa5O7vlWpsh0DKJH1S0sXI7bXKB1ebCu6jtfG+pRPvWzrxvqUT71ttrZX0Wy8TLC37pb3gD78gkszHVpKki+UOIkRr4X1LJ963dOJ9Syfet5qr+P8hyeYAAAAJEUgBAAAk1K6B1GVJB4P/RXrwvqUT71s68b6lE+9bgy37ZHMAAIB6adcZKQAAgCUjkAIAAEho2Zc/QDqZ2Q5Je9x9pEjbaOQy6+6PN25kANB6zOy4uw8V3OOzsgHaLkeKv1itzcwGJd0nKStpm7vfWtA+qsj7ZmbDkra7+75GjxWLmdne4NvtkiYK3xd+/1qPmWUl3Rtc9iv/u7fP3XORPrxvLSz4HDzi7ha5x2dlg7RVIMVfrPQI3pv9RQKp05KG3H0icu937r6u0WPEQmZ2KPq7ZGZHJCmcVeT3rzWZ2WFJh919PHLdF85u8L61tiAQHpV0qCCQ4rOyQdotkOIvVkoUC6SCD4zfRT8sgvsu6dbwHwI0XvDevCTpi+FMRjC7+Lqkfnef4PevNZnZcUnHI4HSXkX+UeZ9a21BoPu0Ip+NfFY2Vtskmwd/sfqiHwaBbPCBj9bXV+J+rkwbGqdPC9+H8Hetj9+/1uXuQwVLdf2SXpT43Gx1wXtwokgTn5UN1DaBlPiLtRx0l7g/VaYNDeDuOXdfV/BfuuHv1YT4/UsFM+uTtEPSnuAW71tr21ZidonPygZqp0CKv1hAY+2R9GIwm8HvX4sLloiOKL9bNpyB4n1rUWY27O5jzR4H2iuQQvpNlbjfXaYNTRAsOeyQtKh8BVqTu48FOYn7Irsv0YKCJddcmS58VjZQO9WR4i9W+k1I+Q+R6NZs5bdrF+ZwoLkOKZ/Umguu+f1Lj0OSjpvZUfG+tap7JfVH8tT6pfmNAhOK5LjxWVl/7RRI8Y9wyrl7zszCZaJcQRu7UFpEsH1+T8HvGb9/LSiY2XhC0lcj70v4fuxQfjcY71uLKVzSC3LbRqObBvisbJy2WdoLPgSK5mrwF6sllcq/OCRpOLwI8jqoZ9MigvfjUJhjY2Z9ZjbI71/LCpPLo+9LNvjfCd631MgWucdnZYO0Wx2pwsJyC67RfJHK5sPKf8iPSXo9+l9gwfR1TvkPjx4KA7aGoPZXn6TwH9ispCEFVbL5/WtNRQqpHpK0I6zhxvvW2oL3Y0T5gPio8sVVw6U9PisboK0CKYm/WEA9hAUAi7UVVFvm96/FBO/d/sitrBYfEcP7BpTQdoEUAABArbRNjhQAAECtEUgBAAAkRCAFAACQEIEUAABAQgRSAAAACRFIAQAAJEQgBQBFBFXZ+5o9DgCtjUAKAIqj6CSAigikAKC4vvDMQAAohUAKAAoEZz5yKC+AigikAGCxPZION3sQAFofgRQALMayHoBYCKQAIIJlPQDV+FizBwAASQVBT1/wJUlHJQ2H7e7+eILH7pF0qMTrDUvaLumCpAlJ90l6zN0JvIA2RSAFIJWCGk997n40uP6dpH5332NmhyVtk5QkkCq6rGdmo5JG3H0ocj0s6atJ/wwA0o9ACkBa7XD3sch1VtLx4PtENaBKLesFQdthSesityck5dw9l+S1ACwPBFIA0urp8JtIBfIXJakwuAnah5UPfvokHS2RTF5qWe9w8DPR5w6GrwegfRFIAUilgqBmh6TxMrNDR9z9Vkkys6yklyTdWqRfqd16O5QPsqKG9NEMGIA2xa49AMvBkKQTxRqC5bpseB0EW9nCc/TMbIdKL+upyPN3iBkpoO0RSAFIpYJAaIek1yNtw5G2bZJyBT+eU35pLmpE5Ytwzs9UBUGX3H3czAY53BhoXwRSAFInCJROm1k2+H4q+AqX7roj3bNhW4Huguuiy3rBvTC3Knz+Hn0UnO2geCfQvsiRApBG45LGJN2rfJA0JGmfmXVLUsFuvpwWB01SJLgqtawXMSJpj5mdVn6n3oiZHTGzvRV+DsAyZ+7e7DEAQN0EOVJH3L0/cu+0pKFwJimoO3WImSUA1WJpD8CyVlh1PFiayxUETZytByARlvYAtIMRMzsk6TXlj3gZCRuCZT3KGABIhKU9AG2NZT0AS8HSHoB2100QBSApZqQAAAASYkYKAAAgIQIpAACAhAikAAAAEiKQAgAASIhACgAAICECKQAAgIQIpAAAABIikAIAAEjo/wOO2DhtvwbXBwAAAABJRU5ErkJggg==\n", "text/plain": [ "
" ] }, "metadata": { "needs_background": "light" }, "output_type": "display_data" } ], "source": [ "difference_m_eff = np.abs(periodic_m_eff - m_eff)\n", "difference_m_eff.show([0, 47], logscale=True)" ] }, { "cell_type": "markdown", "id": "472ab97b", "metadata": {}, "source": [ "## Derivatives" ] }, { "cell_type": "markdown", "id": "d99414fe", "metadata": {}, "source": [ "We can obtain derivatives of correlators in the following way" ] }, { "cell_type": "code", "execution_count": 14, "id": "03007f8a", "metadata": {}, "outputs": [], "source": [ "first_derivative = symmetrised_correlator.deriv()\n", "first_derivative.tag = 'First derivative'" ] }, { "cell_type": "code", "execution_count": 15, "id": "c0311739", "metadata": {}, "outputs": [], "source": [ "second_derivative = symmetrised_correlator.second_deriv()\n", "second_derivative.tag = 'Second derivative'" ] }, { "cell_type": "code", "execution_count": 16, "id": "165550d9", "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAlMAAAGLCAYAAADjxBc3AAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjMuNCwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8QVMy6AAAACXBIWXMAAA9hAAAPYQGoP6dpAABIPklEQVR4nO3df5gcVZkv8O+ZYUmMmZnKTDIkBoh0DD9iEOxk/PUY0dCDIOLibidxVdxnF9JzUTT3RnbacdfLdXWJPRvYC/ostyeg7hr3kvRcZQEF6TYr4oNKMk0WMLKEKY0kBoZkpmcSIEFn3vtHVXWqe/rndFd1Tff38zz9zFSdU1VvVSbdb59z6pQSERARERHRzDTVOgAiIiKi2YzJFBEREVEFmEwRERERVYDJFBEREVEFmEwRERERVYDJFBEREVEFmEwRERERVeCMWgfgNKWUAvAmAMdrHQsRERHNKi0Afi9FJuWs+2QKRiJ1qNZBEBER0ax0NoDDhSo0QjJ1HABeeOEFtLa21joWIiIimgUmJiZwzjnnACX0bDVCMgUAaG1tZTJFREREVccB6EREREQVYDJFREREVAEmU0REREQVYDJFREREVAEmU0REREQVYDJFREREVAEmU0REREQVYDJFREREVAEmU0REREQVqMoM6EqpAIAeEVmfo6zX/LULgC4i4azykG1RE5H+csqJiIiIaqmiZEop5QewEYAGwJejPGJPnpRSMaVUzEq6zEQpnSAppYL2bYqVO2Vk4iRGjp8CAExOCX71+3GMvfoHLJj3J3jrm9rQ3KQAAJ0tc9DZOtfJUIiICrK/XxXC9ysi5ygRqXwnSgUB9InIats6DcCPAVwuIilznR/AEIDlIqIrpYYBdIuIbttuTEQWmL8XLC8xtlYA4+Pj4yU/m++f4s/hjh8fKFpv8+Ur8D+6zy81FCKiquP7FZEzJiYm0NbWBgBtIjJRqK7TDzr2ma+kuWwlRT6l1CgAnz1RMmlm0qUXKheRJBzyiXeeizee2YxbH3o2b50vXnUhrn37UqdCICIqySfeeS66V56VXn5+5AT++859+N8bL8VbOuen13e2zKlFeEQNwbFkymyNym5BsroCdeToFjSlCpTZyx1Lpjrmz8G3Hv9t3nIF4FuP/xbXry0UJhGR8zpb5+bsvntL53ysWtpWg4jICeFwGLquo6+vD36/v9bhlCSVSmF0dBQ+X/1/Vrp9N18PgITZ2tSep86oWVasPCel1BylVKv1AtBSbpBP/GYUR8ZP5i0XAEfGT+KJ34yWu2siIsdMTgmeOpQCADx1KIXJqcqHccxEf38/+vv7EQ6HEQ6HkUql0N/fGPcOhcNhrF8/7V6sivX19SGRSGB0dHZ87iSTSVx++eXo7u6udSiucLqbL83sugsAWF2sboX6ANxSyQ5GjudPpGZSj4jIaQ8/cwRffmB/+ovgF7//DL6++3nccs1KXLlqiWtx9PT0YP369QgEAul1TiQXbhoYGEAoFCpeEUBXVxc6OjqqHoOmabOqhcfv92P79u0z+rcv53p7hZstUxEAq63B6DBamHJpN8uKleezFUCb7XV2uYF2tpR2x0up9YiInPTwM0dw447ktBb1F8dP4sYdSTz8zBHXYtm1a1dGIgUA27dvd+34TojH4yXXDQaD6O3tLV5xBtrb83bKeJKmaTParpzr7RWuJFNKqSiMeahSttW6WaZlVdfMsmLlOYnIKRGZsF4Ajpcb75s75mHh/DML1lk4/0y8uWNeubsmIqqqySnBlx/Yj1wdeta6Lz+w37Uuv1QqhVQqlbFO07RZ2d2TSqXSY5XIebP5ejvezWfOFRWx7spTSvlgzB2VVEpZY6dS9m2sO/WKlTvl3j0v4OiJ1wvWOXriddy75wXeakxENVXOGM93L69+91O2QCCA1atXIxqNZrRQWd02AwMDiEQi0HUdsVgMwWAQuq5j9erV8Pl8iMVi0HUd4bAxneD27duh6zpGR0cxNDSEaDSKgYEBtLe3Y+fOnRkDshOJxIy2A4yxTl1dXdB1HZqmIRQKIZFIQNf1jHissUvhcBh+vx89PT3plpSNGzdi06ZNAIChoaH0vq3jAsjYf77jWqxxZlb3Xjnjpfr7+6FpWvq4wWAwo8zap67r6Za0wcHBnOfV1dWVc30kEil6Dtl0XUcyaXyE79mzB93d3em/k3zX22rhKjduKz5XiEjFLwAhAEM51gcB9MIYKxUwl6Mwkilru96s/ZS8XGJsrQBkfHxcSvXS+Gvy9KGURH/yvKz+yiOyLPxg+rX6K49I9CfPy9OHUvLS+Gsl75OIyAn3PXko4z0q3+u+Jw+5Es/Y2Jj4/X6BkcdJIBCQeDyeUScej4vP58tYF4lEctaxb+vz+aS3tze9HIvFxO/3V7xddozBYDC9HI/Hpx3Dvo+hoSEZGhpK7z/73GKxmESj0fTy8PBw+lwLHbe3tzfjmoyNjYmmadOuZS6hUChj22AwmF4OhUIZ+xgeHpZAIFD0vPKtL3QOw8PD0/6dA4FARmw+n0/GxsbSy/mu90zjrsT4+Lj1d9wqxXKNYhUKbgz4YYyFGjYPGAUQMss06z9T9itrH71WkgSjBSv7GAXLS4ix7GTK7o+TU/L480flvicPyePPH5U/Tk7NaD9ERE54/PmjJSVTjz9/1NW4rA80K7GKxWIZ5dkJT3b50NCQ9XmRFggEMj4kc31Yl7vd0NCQaJqWUT8ajUowGBSR/B/u8Xh82na5YorFYhIIBDISBusDP99xx8bGpp2DiIjf7y+aTA0PD0/bdnh4WMbGxmR4eDhnzPb95juvXOuLXbtc/z5WLJZAIJDxb5/relcSdyXKSaYq6uYTo7stCWDa413EGB+lSthHwftli5U7rblJudI0TkQ0E+84rx1L2ubixfGTOcdNKQCL2+biHee5O3jZ7/enu9HC4TA2bdqU0dXU09OT7gocHBzMKLNkD2DWNA3Lly/PWJer66uc7fbu3Yv29nYMDg6my1KpFLq6uoqeYyl31wWDQUSjUSxYsAB+vx8bN25Eb29vuusv13ETicSMB28nk8lp21px7tq1K2fMPp8P8Xg83d2W77yy18/k2vl8Pui6jl27dqG9vR2jo6NFuy8TiURFcbvBtakRiIio+pqbFG65ZiVu3JGEAjISKuvb7C3XrEw/U9RJqVQKiURiWmIUiUTQ39+PVCqV/qAPhULpOajyfZjO9O61mWyXK5nLRdf1sj+04/E4dF1HIpFAJBLBsWPH0sldruPak5PZoNRrBxjjnvbs2YNYLAYA2LlzZ8H6s2UwutuTdhIRUZVduWoJ7vqkH4vbMqdrWdw2F3d90u/qPFP5bmv3+XwZLSaapiEQCCAcDmPNmjUuRTfdmjVrcn5g5/sQtwZPl2pgYACpVAo+nw+hUAhDQ0NIJBIFj+v3+6fdEVmqfNvqup73mMlksqSWuGzlXjvrbj0rkbLWWTHkkkwmqx63E5hMERHVgStXLcHPwutw60dXAQBu/egq/Cy8ztVECjC6krJbVvJ144XDYSQSibyPR8luscpOEvIlHOVs5/f7EQwGMTAwkFFufbhb3VJAaa1SuY5l37e1z0LHtRIve5l1F1yxJMva1j7jfCqVSidpVreqxeoWLKd1yVLs2mXHav272Ndb66xrnOt6VztuRxQbVDXbX6hwADoR0Wzy9KGULAs/KE8fSrl+7LGxMYlGo+nB55FIZNpdadlylQ0NDUkgEBAAEgqFZGxsTHp7ewWA+P1+icViOevMdDuR03fPRaPRaYPhe3t7pbe3N31XnjVI2trH8PBwzrhFjAHZ1j5jsZhEIpGMAdiFjmutt8r8fr/4fL5p9XIpdj7Wfu0D8/OdV771hY5lvxb2Y0QiEQmFQhKLxSQej6fvyrPf8Zh9vWcad6XKGYCuRGrz7Ca3mM/nGx8fH0dra2utwyEictQzh8fx4a//DA9+9r2ef9Dx4OAgAoHAjAdbEzlpYmICbW1tANAmxiTgeXEAOhHRLDYycRIjx0+ll58fOZHx09LZMgedrbV/BFY4HMby5cvTEzsykaJ6wJYpIqJZ7J/iz+GOHx8oWm/z5Ss88cSGRCKRHlMTCoWYTJFnldMyxWSKiGgWy26ZyscrLVNEswW7+YiIGkRn61wmSUQ1xqkRiIiIiCrAZIqIiIioAkymiIiIiCrAZIqIiIioAhyATkQ0mx1/0XgV07LYeBFR1TGZIiKazfZ+C3j0a8XrXfYF4AN9zsdD1ICYTBERzWZr/gq44KrTy0efA763Cfiz7cBC2ySdbJWaJhwOQ9d19PX15X3YMlEpOGaKiGg2a1kMvOnS0y8rgVp4fuZ6F5IpXdfR09MDpRQWLFiAcDicfq1fvx5KKYTD4XR9a32t9PX1IZFIYHR0tKL9uHketb5mlBtbpoiIqCp8Ph+i0Sj27t2LNWvWIBKJZJSnUils2rQpvdzV1YWOjo6KjzswMJB+1l85NE2Dz+er+PjVOo9suc7LqWNRZZhMERHVi6lJ4PdPGr///klg8cVAU7PrYbS3t+dcr2kaurq60svBYLAqx4vH4zNKpoD8sZajWueRLdd5OXUsqgy7+YiI6sH++4E7LwUe/O/G8oP/3Vjef3/tYjLpuo5UKgUACAQCVdtvKpVKj3uqJ/V6XvWMLVNERLPd/vuBXZ8Czr8S+PNvAp0XASO/Bh67zVi/4V+BlR+pWXiJRAKBQACapqUHeieTyXSX39DQEABgcHAQ4XAYfr8fPT09iMfjAIBIJIKBgYF0K5Ku69A0De3t7dB1Hbqup8di9fX1QdO0vLH09/cDQLp7L9d4qXA4jK6urvRxQqFQ3tg2btw47TwGBgYQiUSg6zpisRiCwSCSySQuv/xy+Hw+xGKx9DUAgD179qC7uzudaCYSiZznpev6jI7l8/lynhNVkYjU9QtAKwAZHx8XIqK6M/lHkX9aJfLdjSKTk1llk8b6f7rYqOeSQCAgPp9Pent7JRgMCgAZHh6eVi8ej4vP58tYF4vFxO/3y9DQkAwNDUlvb6/EYjGJRqPpOsPDwxKJRNL78Pv9JcXV29ub3k5EZGxsTDRNk3g8nhG7fTkYDKaXc8WW7zxyrbMfOxAIZCz7fD4ZGxvL2D7Xec30WPnOifIbHx8XAAKgVYrkGuzmIyKazQ4+DqR+B6z9PNCU9Zbe1ASs3QKkDhr1XBQIBBCJRBCLxaYNRLfkGvytaRp0XYff74ff709vG4vF0l2FPp+v7O7CVCqF/v5+9Pb2ZhzLHkMymcTevXsz9t3d3Y1oNFowtlznEQgEMDo6mm59sra3RKPRjNYhn8+HRCJR9DzKPVaxc6LqqEo3n1IqAKBHRKbdr6mUsrclaiLSX81yIqKGduIl42fnRbnLrfVWvRooN/HJThiCwSCi0SgWLFgAv9+PjRs3ZiRFpUgkEgW7/wBg7969aG9vx+DgYHpdKpXKGDRfzt1/oVAI0WgU0WgUiUQCGzZsyNiPruvYtWsX2tvbMTo6WtEUDfmOVco5UeUqSqaUUn4AGwFoAKb9hZmJUDoBUkoFlVIREQlXo5yIqOHNP8v4OfJr4JwcH5Ajv86sVwPVmBAzHo9D13UkEglEIhEcO3YsZ4uXrusVT3dQrTvmenp6sHr1akSjUei6npFU9vf3Y8+ePenxUzt37iy4r2LnVehYAO8CdFpF3XwikjQTm3ieKmEAg7b6gwBCVSwnImpsy94DaOcag82npjLLpqaAx24HtGVGvVlqYGAAqVQKPp8PoVAIQ0NDebvE7F1ddn6/P91NmM+aNWty3kE307vqfD4ffD5fxuB54PTdelYiZa0D8sefb32xY1X7nCg3x8ZMKaU0AD4Ryf4X05RS/krLHQmaiGi2aWoGrvgH4LmHgXs/DrzwBHDquPHz3o8b66/4quvzTZXSZVUsubHXGxgYyFhntdJY3WVA4dYbKxGz70fXdSSTyXQcfr8fwWAwo04qlSqayBQ6j56eHoTD4YyWIeva2Lez1lnnku+8yj3WTM+JylRshHopLwBBAENZ6/zG7qfVHTPrV1ReRmy8m4+I6t+v/t24q++W1tOvf7rYWO+S4eFhCYVCAkA0TZNQKJT3rrGhoSEJBAICQEKhkIicvoPNWmfdARiNRiUajUosFpNYLCaRSCTjzrfe3l7p7e3NuOMvn0gkkrE/v98vPp9PYrFYxv6setb6fLHlOg+7sbGx9F1/2XGEQiGJxWISj8dleHhYAoFAxjlkn9dMj5XvnKiwcu7mU2IkHBVRSgUB9InIatu6AIC4iKisusMAIgD0SspFJPNrSv7YWgGMj4+Po7W1tfyTIyKaLaYmgeS/GhN2fvh/A/5P1WQGdKJ6MDExgba2NgBoE5GJQnXrbtJOpdQcAHNsq1pqFQsRkauamoE3vd34/U1vZyJF5BInk6l8HebtZlml5fn0AbillACJiGa94y8aL8vR5zJ/WloWGy8iqjonkykdMAaii0jKtl4zyyotz2crgNttyy0ADpUfPhHRLLD3W8CjX5u+/nubMpcv+wLwgT53YiJqMI4lUyKSUkrpMFqSUlllSQCotDzPcU8BOGUtK6XyVSUimv3W/BVwwVXF67FVisgx1Uqm2vOsj8C4M8+adDMEY+6oapUTETU2dt8R1VxFd/PZZkAPwpgBfQDGFAkDtjq9MFqWNAAdkjV7eaXlJcTIu/mIiIioLOXczVeVqRG8jMkUERERlaucZMqxGdCJiIiIGgGTKSIiIqIKMJkiIiIiqkDdzYBORNRIXn71Zbz82stF6y16wyIsmrfIhYiIGg+TKSKiWSz2XAx3/eddRevdeMmN+PSln3YhIqLGw2SKiGgWW3/+erz/nPenl/VxHX2P9WHr2q3wtfnS6xe9ga1SuYTDYei6jr6+Pvj9fte2pfrCZIqIaBZbNC93952vzYeVHStrEBHQ398PADh27BgAoK+vDwMDA+jt7a1JPIX09fXhvPPOQ09Pj6vb2llJWSwWq2g/XjtWI2EyNQtMTgme+M0oRo6fRGfLXLzjvHY0N/ExOUSUaXJqEr86+isAwK+O/goXLLgAzU3NrsbQ09OD9evXIxAIpNetX7/e1RjKoWkafD5f8YpV3tauq6sLHR0dFe8n28DAAEKhkCvHanRMpjxqZOIkRo6fwuPPH8XAYzqOnng9XbZw/pkIrfXhPW9ZiM6WOehsnVvDSInICxIHE9i2dxsOnzgMAPj7X/w97nnmHty85mYElgWKbF09u3btQjQazVi3fft2bNq0Kc8Wtdfenu+JaM5uawkGgxXvI5d4PD4tmXLqWI2OUyN41Hd/+Tt8+Os/w60PPZuRSAHA0ROv49aHnsWHv/4zfPeXv6tRhETkFYmDCWz5yRas0FZgx4d24Jcf/yV2fGgHVmgrsOUnW5A4mHAtllQqhVQqlbFO0zR0d3e7FkOjS6VS6e48cgdbpjzqY13n4Lu/PDgtkbJbOP9MfKzrHBejIiKvmZyaxLa923DZ2ZfhjnV3oEkZ35EvWXQJ7lh3Bzbv3oxte7fhA+d8wJUuv0AggNWrVyMajWZ09WW3kITDYXR1dUHXdWiallHe398PTdPSrT721pT+/v5015qu6+lxWIlEAuGw8ejW7du3Q9d16LqO4eHhaS1l1pguaz+jo6Mln18p2+Y6t8HBQYTDYfj9fvT09CAejwMANm7cmG61GxoaAmB0z0UikfTYpmAwiGQyicsvvxw+ny893imZTAIA9uzZg+7u7vT1TiQS6fO3rklfXx90XZ/RsXw+X8F/LwIgInX9AtAKQMbHx2U2efz5o7Is/GDR1+PPH611qERUQ08ceUJWfXuV7BvZl7P8yZeelFXfXiVPHHnClXjGxsbE7/cLAAEggUBA4vF4Rp3sdcFgML0cCoUkEolklFnLoVAoY7vh4WEJBALp5Xg8Lj6fL6OOz+eToaGh9HJvb2/G/sfGxkTTtGkx5lLKtoXOLRaLid/vl6GhIRkaGpLe3t6MuO1yrbMfOxAIZCz7fD4ZGxvL2N7v9087h5keK9851bPx8XHr77hViuQabJnyqJHjJ6taj4jq08uvGhN2rtBW5CxfsWBFRj2naZqGoaEhJJNJ7Ny5E4lEAt3d3RmtHnv37s1oteru7kY0GoXP58PAwID1RRgAEIlE0N7eDl3Xp43H8vl8GB0dRSKRQCAQSNez79vn80HXdfj9fqRSKfT392fsv9RB5KVsW+jcAoEANE1LxwIg/TPX8QOBAEZHR5FMJtP1NE1Ll0ej0YzxWj6fD4lEouiYqHKPVeycyMBkyqM6W0obVF5qPSKqT9a0CAdSB3DJokumlR8YO5BRzy1+vz/9wRwOh7Fp0yYEg0Hs3bsX7e3tGBwcTNdNpVLo6upCMpnMSBiA0x/+u3btypkI+Hw+xOPx9Ad79vaapqW74hKJxLTyUpWybaFzs8dbqlAohGg0img0ikQigQ0bNmTsx0ow29vbMTo6WlZ3ZanHKuWciMmUZ725Yx4Wzj+z6JipN3fMczEqIvIaf6cfS+cvxd1P3Z0xZgoApmQK9zx9D5bOXwp/p/OTSqZSqZytI5FIBP39/RkD03O1oNg/sGeqGnfXVapad8z19PSkx59lt7j19/djz5496fFTO3fuLLgvXdcLJnKFjgXwLsBieDefR92754WCiRRg3NV3754XXIqIiLyouakZN6+5GY8eehSbd2/GvpF9eOUPr2DfyD5s3r0Zjx56FDevudm1+aasgdXZfD4fNE3DmjVrct5lZu+Ky1WWb7tkMllyK0m+/Vdr20LnNhM+ny/d9WlPEq279ewTb1qxWYPSs+VbX+xY1T6nesVkyqM+8c5z8eBn34svXnUhFs4/M6Ns4fwz8cWrLsSDn30vPvHOc2sUIRF5RWBZALe//3YcSB3AdQ9dh3f927tw3UPX4UDqAG5//+2uzzOV3cI0ODiYbtnw+/0IBoMYGBhIl6dSKSSTSfh8PoRCofQdc1aZlWgFAoGMfVvdgvZWk+yuLnsCZO3ffmxd15FMJosmSqVsW+jcCil07J6eHoTD4ZznaN/OWmclOVY3oLXOapUq91gzPadGo+yD6eqRUqoVwPj4+DhaW1trHc6McAZ0IirF5NQkvnfge/j7X/w9/ue7/if+bMWfuToDeiqVwq5du7BmzRrs3LkTHR0dOHbsGDo6OqY9SiYcDqOjoyM9BYL9A7xY2fLlywEAw8PDiEQiAIzEKhwOI5FIIBQKIRKJYOvWrejv74ff70dfX196P9bUC4DRLbh161akUilEIpGi3VmlbJsrfmvqhmQyiVAohHA4DJ/PNy3u7GkcUqkUtm7dmj5PexzDw8Po7u5OD4S3Zp+3pi2wpkVYvnw5QqHQjI9V7N+kXk1MTKCtrQ0A2kRkolBdJlNERHVk/7H92PjgRuz88M6aPZuPqB6Uk0xxADoR0Sz28qsv4+XXTk97oI/rGT8ti96Q+4HIRFQ5tkwREc1i/7zvn3HXf95VtN6Nl9yIT1/6aRciIqoPbJkiImoQ689fj/ef8/6i9Ra9ga1SRE5hMkVENIstmsfuO6Ja49QIRERERBVwpWVKKWU9XloD0AFgq4ikcpQDgCYi/bblouVEREREteJ4MqWU6gUwYCVPSikNQARAj7kcgi1BUkoFlVIREQmXUk5ERERUS47fzaeUiotId751SqlhAN0iotvKx0RkQSnlJRyfd/MRERFRWcq5m8+NMVPtZuvUNGYrlc+eKJk0pZS/WHn1QyUiIiIqjxvJVBhARCkVV0ppSql0Fx+AfI+wTpllxcqJiIiIasrxMVMiklBKdQOIAxgDsN7W0tSeZ7NRsyxVpHwapdQcAHNsq1rKjZmIiIioVI63TCmlfAD8ABYAGAAQy7o7r9r6AIzbXoccPBYRERE1ODe6+SIi0i8iKRHpAdANIGomWaN5tmk3y4qV57IVQJvtdfaMIyciIiIqwtFkyhwknjF4XEQSAPoBBKwyc6C5nWaWFSufRkROiciE9QJwvJJzICIiIiqkVjOgDwPQzbmndOQY/yQiyWLlDsdIREREVJSjyZSZ8PhztCytNluoAGMCz6BVYI6nsk/IWayciIiIqGbcmLRTgzEoHACOIffjZHph3LmnAejInt28WHmR43PSTiIiIipLOZN2Op5M1RqTKSIiIiqX12ZAJyIiIqpbTKaIiIiIKsBkioiIiKgCTKaIiIiIKsBkioiIiKgCTKaIiIiIKsBkioiIiKgCZ9Q6AJodJqcET/xmFCPHT6KzZS7ecV47mptUrcMiIiKqOSZTlNfIxEmMHD+Fx58/ioHHdBw98Xq6bOH8MxFa68N73rIQnS1z0Nk6t4aREhER1Q6TKcrru7/8He748YGcZUdPvI5bH3oWALD58hX4H93nuxkaERGRZ3DMFOX1sa5zsHD+mQXrLJx/Jj7WdY5LEREREXkPkynK67fHXs3o2svl6InX8dtjr7oUERERkfcwmaK8Ro6frGo9IiKiesRkivLqbCltUHmp9YiIiOoRkynK680d80oaM/XmjnkuRUREROQ9TKYor3v3vFDSmKl797zgUkRERETew6kRKK9PvPNcdK88q6R5poiIiBqVEpFax+AopVQrgPHx8XG0trbWOpxZizOgExFRI5mYmEBbWxsAtInIRKG6bJmikjQ3Kbx7eUetwyAiIvIcjpkiIiIiqgCTKSIiIqIKMJkiIiIiqgCTKSIiIqIKMJkiIiIiqoBrd/MppXoBpMzFUREZtJWFbFU1EenP2rZgOREREVGtuNIypZSKAxgUkQEAewHEbGUhGAnSgFmuK6UipZYTERER1ZLjk3aaydByEQnb1vlFJGn+PgygW0R0W/mYiCwopbyE43PSTiIiIipLOZN2utEyFQEQt6+wJVIaAJ89UTJpSil/sXKH4iUiIiIqmaNjpsxkSIOR/FjjnuytVL48m6YKlNnLkzmOOQeA/WFxLSUHTERERFQmp1umrISo3TbmKa6UssZMtefZbtQsK1aeSx+AcdvrUNlRExEREZXI6WTKSnj2WitEJAEgqJQq1PJUia0A2myvsx06DhEREZHjUyPoWT8tKQD+HOst7TBan0aLlE8jIqcAnLKWlVIlhkpERERUPkdbpmwDx7NboTTzpw6kx1Zll+sllBMRERHVlBt38yWRe3xTUkRSMJKiaeUiUrS8umESERERlc+NZCoMYL21YN7VN2hrtYoACGaVh23bFyunBjE5Jfj58DH8+77D+PnwMUxOOTtHGhERUSkcn7QTOD1xp7Vsn8DTLLceNaMB6Ci3vMixOWlnHXj4mSP48gP7cWT8ZHrdkra5uOWalbhy1ZIaRkZERPWonEk7XUmmaonJ1Ow2MnES9z15GLc+9GzeOl+86kJc+/al6Gyd62JkRERUz7w2AzrRjH3nFwcLJlIAcOtDz+I7vzjoUkRERESZmEyRp61cUlprYqn1iIiIqo3JFHna65NTVa1HRERUbUymyNM6W0obB1VqPSIiompjMkWe9uaOeVg4/8yCdRbOPxNv7pjnUkRERESZmEyRp9275wUcPfF6wTpHT7yOe/e84FJEREREmZx+Nh9RRT7xznPRvfIsPP78UQw8pmckVgvnn4nQWh/e85aF6GyZU8MoiYiokXGeKZo1JqcET/xmFCPHT6KzZS7ecV47mpv4IGsiIqq+cuaZYssUzRrNTQrvXt5R6zCIiIgycMwUERERUQWYTBERERFVgMkUERERUQWYTBERERFVgMkUERERUQWYTBERERFVgMkUERERUQWYTBERERFVgMkUERERUQU4AzpRGfhIGyIiysZkiqhEDz9zBF9+YD+OjJ9Mr1vSNhe3XLMSV65aUsPIiIioltjNR1TEyMRJDDw6jP+2I5mRSAHAkfGT+G87khh4dBgjEyfz7IGIiOoZkymiIr7zi4O49aFnC9a59aFn8Z1fHHQpIiIi8hLXu/mUUnER6c5aF7ItaiLSX045kZNWLmmtaj0iIqovrrZMKaWCAAJZ60IwEqQBERkAoCulIqWWEznt9cmpqtYjIqL64loypZTSAPhyFIUBDFoLIjIIIFRGOZGjzmwu7b9JqfWIiKi+uPnuvwHAgH2FlWCJiJ5VV1NK+YuVOxYpkc3+IxNVrUdERPXFlWTKTHz25ijK1VIFACmzrFg5keOue9cyfPGqCwvW+eJVF+K6dy1zKSIiIvIStwagrxGRAbOlya49T/1RsyxVpHwapdQcAHNsq1pKD5Nous7WuQhdthzndszjPFNERDSN48mUUipoDhx3Sx+AW1w8HjWIK1ctQffKxZwBnYiIMjiaTJktUakCVUbzrG83y4qV57IVwO225RYAhwrEQFSy5iaFdy/vqHUYRETkIU63TG0AsNw2WHw5ACilegHoABLmsiYiKdt2mlmuFymfRkROAThlLSvFVgMiIiJyjqPJVHb3nlLKByBkn3RTKaUjx/goEUmWUk5ERERUS25PjKPlWBcBELQWzEk6w2WUExEREdWMEhF3DmQkQethzIA+CCAqIlY3Xy+MlicNQIeIhLO2LVhe5LitAMbHx8fR2srHfRAREVFxExMTaGtrA4A2ESk4kaBryVStMJmiejM5JbyjkIjIYeUkU64/6JiIZu7hZ45wrisiIo/hw8SIZomHnzmCG3ckMxIpAHhx/CRu3JHEw88cqVFkRESNjckU0SxwJPUa/u6+Z5CrU17M19/d9wyOpF5zOTIiImIyRTQLbHvkv3D0xOsF6xw98Tq2PfJfLkVEREQWJlNEs8Cl52hVrUdERNXDZIpoFlg4f07xSmXUIyKi6mEyRTQL7D9S8K7csusREVH1MJkimgWue9cyfPGqCwvW+eJVF+K6dy1zKSIiIrJwnimiWaCzdS5Cly3HuR3zOM8UEZHHcAZ0olmGM6ATETmPM6AT1bHmJoV3L++odRhERGTimCkiIiKiCrBliojKxq5GIqLTmEwRUVn4sGUiokzs5iOikvFhy0RE0zGZIqKS8GHLRES5MZkiopLwYctERLkxmSKikixfNL+q9YiI6gWTKSIqyfDLJ6paj4ioXjCZIqKS3HzFBVg4/8yCdRbOPxM3X3GBSxEREXkDkykiKskS7Q346rWroABkzyhlrfvqtauwRHuD+8EREdUQkykiKtmVq5bgrk/6sbhtbsb6xW1zcdcn/ZxniogaEh90XMTk1CSSI0m8/OrLWDRvEfydfjQ3NVc/UKJZxAszoHshBiKqX3zQcZUkDiawbe82HD5xOL1u6fyluHnNzQgsC9QwMqLaqvXDljkLOxF5iSstU0qpXvPXLgC6iISzykO2RU1E+sspL3LsGbVMJQ4msOUnW3DZ2ZfhhrfdgBXaChxIHcDdT92NRw89itvffzsTKqIasGZhz37nstqk2N1IRNVQTsuU48mUUipiT56UUjEAEJH15nIItgRJKRUE0GVtU6y8hOOXnUy9+MqL+MQPP4GL2i/CnevuRJM6PbRsSqbwud2fw7Ojz2LHh3Zg8RsXl7RPIqrckdRruOYbPys4eejC+WfigZvey4HwRFSRcpIpRwegK6U0AAHzp2UrgKBSymcuhwEMWoUiMgjA3hJVrLzqvv7k1zHy6gg2vW1TRiIFAE2qCTdcfANeevUlfP3JrzsZBhFl4SzsRORFbtzN5zNfFt1abyZZPhHRs7bRlFL+YuWORAvgvNbzAAArtBU5y1csWJFRj4jcwVnYiciLHE2mRCQlIgtEJGlbbSVWOjKTLLsUpidhucod8ZuJ3wAADqQO5Cw/MHYgox4RuYOzsBORF9VinqkeAAmztak9T51Rs6xY+TRKqTlKqVbrBaCl3AA/+/bPonNeJ7Y/tR1TMpVRNiVTuPvpu3HWvLPw2bd/ttxdE1EFOAs7EXmRq8mU2TUXALDewcP0ARi3vQ6Vu4PFb1yMvnf04aeHforNuzdj38g+vPKHV7BvZB82796Mnx76Kb7wji9w8DmRy7w2C/vklODnw8fw7/sO4+fDxzA5Vd/z9hFRbq5O2qmUigNYLyIpc9kPYEhEVFa9MQCbYHQF5i03B6NnH2MOgDm2VS0ADs1k0k7OM0XkTV6YZ8oLMRCRczw1NUL6QEpFAUTsg8nNAeZjABZYCZa5XgCshpFM5S3PGouV77icAZ2oDtVyBnTOdUVU/zyXTJlzRVnjpGBOi6CJSFIpNQygOyvJEqs1qlh5CceuKJkiIrLjXFdEjcEz80wB6Uk2NRhTIQTM5TBOT5EQARC01Q+Z5SixnIjINZzrioiyOfpsPrMbL5arTER6zJ8DSqlea6ZzAB322c2LlRMRuYlzXRFRNkeTKXOcU9HuuGLP2ivnWXxERE7y2lxXtRw7RkQGR5MpIqJ6c/MVF+DR514uOmbKjbmueEchkTfUYtJOIqJZyytzXVl3FNoTKQB4cfwkbtyRxMPPHHH0+ER0GpMpIqIyXblqCe76pB+L2+ZmrF/cNteVaRGOpF7D3933zLSpGQBAzNff3fcMjqReczQOIjKwm4+IaAauXLUE3SsX12S8Ujl3FN624VLH4yFqdEymiIhmqLlJ4d3LO1w/rtfuKOQgeGp0TKaIiGYZL91RyEHwRBwzNStMTk1iz4t78EP9h9jz4h5MTk3WOiQiqqGbr7gAC+efWbCOG3cUchA8kYEtUx7Hhy0TUTbrjsIbdxiPJ7UPRLc615y+o7DYIHjAGAR/ydkaH6tDdY8tUx6WOJjAlp9swQptBXZ8aAd++fFfYseHdmCFtgJbfrIFiYOJWodIRDVS6zsKvfZYnckpwc+Hj+Hf9x3Gz4ePYXLK+efOElnYMuVRL77yIrY+sRXvO/t9uGPdHWhSRt57yaJLcMe6O/C53Z/D1574GlYtXIXFb1xc42iJqBZqeUehlwbBc9wW1RpbpjzqH/f8I0ZeHcGmt21KJ1KWJtWEGy6+AS+9+hL+cc8/1ihCIvIC647CP710Kd69vMO1u+i8Mgie47bIC5hMedRrfzQm21uhrchZvmLBiox6RERu8sIgeK9NXsquxsbFZMqjPvqWjwIADqQO5Cw/MHYgox4RkZu88FgdL43beviZI3hvZDf+YvsvsPneffiL7b/AeyO72TLWIJhMedS6c9dh6fyluPupuzElUxllUzKFe56+B0vnL8W6c9fVKEIianS1HgTvlXFb7GokDkD3qOamZty85mZs+ckWbN69GddffD1WLFiBA2MHcM/T9+DRQ4/i9vffjuamZlfimZyaRHIkiZdffRmL5i2Cv9Pv2rGJyLtqOQjeC+O2JqcEX35gf96uRgXgyw/sR/fKxa5cE85GXxtKpL77dJVSrQDGx8fH0draWutwyuaFeaa8EAMRUbYjqddwzTd+VrCrb+H8M/HATe91rLvxoaeP4MbvJovWu+sTflx1sbMtdbyrsbomJibQ1tYGAG0iMlGoLpOpWaCWrULWXFeXnX0ZbnjbDVihrcCB1AHc/dTd6dYxJlREVCtWFxuQe/JSp7sb/+pbT+A//uvlovU+cMEifOuv3uFYHNZ1yP5Ed+s62NVL61g5yRS7+WaB5qZmdC3ucv24nOuKiLzOGreV3SKz2KUWmTf8SWlfbEutNxNemo3eK61jbid0bJmivD7/k8/jkYOPYMeHduCSRZdMK983sg/XPXQdrlh2BW57/22Ox8NxW0SUT61aQ7zQ1fjpHUP44TMvFq33oVWL8c+fXO1IDIB3WseqldCxZYqqwktzXXll3BYTOiJvsiYvdZsXnpP42h8mq1pvJrzSOpYvobPurHQqoePUCJSXV+a68sozChMHE7j6+1fjr3/01wg/FsZf/+ivcfX3r3b9GYmTU5PY8+Ie/FD/Ifa8uAeTU869QRJRcbWeImLDmnOqWm8mvvLg/pLm/PrKg/sdi6HYnZWAcWelE5OpsmWK8rLPdWUfMwW4N9eVV8Zt2QfiR94XyRiIv+UnW1wbiM8WOiJvquUUEVe8dTGWtM3Fi+MncyYSCkZid8VbnXuP9ELr2CO/enHaXF92AuDI+Ek88qsXq35nJVumKC9rrqtHDz2Kzbs3Y9/IPrzyh1ewb2QfNu/ejEcPPYqb19zs6IeoF55RODk1iW17t+Gysy/DHevuwCWLLsG8P5mXTuguO/sybNu7zfEWIrbQZfJCC50XYiDvqNVzEpubFG65ZiWA3LPRA8At16x0NB4vtI7t2vtCVeuVY1a0TCmlQrZFTUT6axZMgwksC+D299+ObXu34bqHrkuvXzp/qSutMV4Yt7X7d7tx+MRhRN4XyZnQXX/x9bjuoeuw+3e70f3mbkdiyE7oslvoNu/ejG17t+ED53zA0eSWLXTeigHwRishY6h9DFZX4/964Bm8/IdfQ51xHPLHFiz6k4vwv65Z5XhXoxdax2p5Z6XnkykzkUonUEqpoFIqIiLhGofWMALLAvjAOR+oyZvER9/yUTx2+DEcSB3IeUehG+O2vv/89wEUT+i+//z3HUummNCd5oWEzgsxWHHUOqFjDN6J4YyWX6FtxTa88srv0+va3vgmnNHyNwCcTaas1rEbdyShMIWmeb9JJ3RTr54HoMnx1rEvfXglnvjtqDl2awrNthgmzRgWzj8TX/rwyqofezZ084UBDFoLIjIIIJS/OjnBmuvqQ74PoWtxl2vftrzwjMI3nGHceVJsIL5VzwnlJHROue/AfTh84jBueNsNeRO6wycO474D9zkWgxe6XL0QA+CNbl/G4L0Yzl9wfkYM5y8437UYrly1BJ+++lVoF9yGecu24w1L78W8ZduhXXAbPn31q463jll3Vp7R8gzaVmzLiKFtxTac0fKMc3dWiohnXwA0I8Rp6wWAv8R9tAKQ8fFxodkp/tu4XPzti+WmxE3y5EtPyonXT8iTLz0pNyVukou/fbHEfxt39PhHThyRdbvWyWcSn5HJqcmMssmpSflM4jNy+a7L5ciJI47FsOU/tsiqb6+SfSP7cpY/+dKTsurbq2TLf2xxLIZrvneNrPr2Knnl9Vdylp94/YSs+vYqueZ71zgWw+CzgyVdh8FnB+s6hj9O/lE+OPhBuSlxU86/yZsSN8kHBz8of5z8I2NgDK7EIHL6vfozic/IvpF98srrr8i+kX3ymcRnXHmvzozh05kxxD9ddgzj4+Ni5hutUiTX8HrLlC/P+lSBMqoz1ritA6kDuO6h6/Cuf3sXrnvoOhxIHXClO2XxGxej7x19+Omhn+YciP/TQz/FF97xBUfvJvybrr9B57xObH9qe84WurufvhtnzTsLf9P1N47F8Jcr/xJA8RY6q54T/mX/vwAo3kJn1avXGLzQSsgYvBODNQygWAy7f7fbsRjsd17fue7OjBbbO9fdifed/T587Ymv4cVXik8uOlOTU5PY9vOv4LKz1+LOdV/PjOHyr+OypWux7edfcaTV2OvJVHue9aP5ypRSc5RSrdYLQItj0ZFrAssC+MFHf4BvfvCbiKyN4Jsf/CZ+8NEfuDYWgQkdcO2Ka0vqcr12xbWOxeCFhM4LMXghoWMM3onBC8MAvHDn9X3P/T8cPjWKG94Wyp1Uvm0TDp8axX3P/b+qH9vrydRM9AEYt70O1TYcqpZajduyNHpC54WpMryQ0HkhBi8kdIzBOzHMa55TUgxWPSe89odXAJRw57VZzwnxJ+4oKQarXjV5PZkazbO+vUDZVgBtttfZDsRFDYoJHRM6L8TghYTu2hXXYumcdtz9nwO5Y3hqO5bOaa//GJZfg6XSjLufyhPD09uxVJpx7fJrHIuhb855WDoJ3J1nGMA9T9+NpZNGPaesbzI6i4oldFY9J3zklRMlxWDVqyZPP+hYKaUBGAOwQERStvUCYLWIJEvYBx90TFRltZ7Txwu3oScOJrBtzz/isO02dDdjSE/PsHQtrn/bJqxYsAIHxg7gnqe249HDjzmf3E5NInHXpdgyH7js7Mtw/cU3nI7h6bvx6KFHcfsJIHDjPsCpv43xw0j8awBbWs/IH8PEJAKfigNtS52JYehfkPhxGFvO6jRjuN4Wwz1GDC+NIHB5BFjtUOvUN7qQePWF4jHMOwe4aY8jIUx+Yw2unncKK5ZdhjvW3TntiRmbd38OBw4+ih+8OgfNN+11Joa938LVT91ePIa3bUHzmr8qur9yHnTs6WQKAJRSwwC6RUS3rRMRKWmyCiZTRPWp1gkd9t+PyUf+FsmTL+Hl5mYsmpyEf+5ZaL7iH4CVH3Hl+IkHe7Btydk4PHX6ERpLm+bi5iOHEPhw1Nk49v87sOtTSFwTwbbf3jc9sX3znyLwwBeADf8KrPxTZ2LY9ZfA/vuKx7DyWmCDQ2OWvtEFHH0OiU9+F9v+8xvTY7jkMwjs+CSw8HzHEhkM/QvwwOeKX4dr7nQuofNCUlnlBL/ekqnsSTszlkvYnskUUT2amgQOPg6ceAmYfxaw7D3OtYBk238/sOtTwPlXAms/D3ReBIz8GnjsNuC5h80EwsFEZvwwcPflwJJLMLlhB5JH951OKhdeiuZdnwSOPAXckHCuRea764EDjwB9hzH5J2+Yntj+4VVg69nAiiuAT8SciWHHeuD5EmJ4yxXAJx2KwUxkcH0Ck0v902M4PATc0+1sIjM1Cdx5KdD5Vkxu+E6Ov4frgJH9wOeedO7/iBlDYuE52DZ3cnqL7WtNCBw95GwMQFW/ZNRVMgUASqleGNMhaAA6pIzZz5lMEdWh/fcDj/wtkPrd6XXauYAbrUK2Dy587N+AJtvQ06kp4N6PO//BZbbI4PoEcE7X9PIXnjA+wJ1skdn5l8CvS4jhomuBjQ7FYLaOFY3BydYxL/w9AFkJ/hZbgn+7Owm+LYbJFR9E8uKr8fLcN2LRyVfgf/oHaD7wI3diMOOoRqtx3SVTlWAyReSARm4VsrVEFPwAd7IlwtYigznzp5efOu58i4ytdQwf+785koi/cL51jInM9DimfclYBlzxVXeO75UYgKq8RzGZsmEyRVRljd4qZI6RKZrIODlGxgstMoA3kggvxGDFUSdJRF3EUAVMpmyYTBFVEVuFvBGDF5JKixeSCC/EANRNEkEGJlM2TKao7tTqDdsLH+BeaBXywnUAvNMiA3gjifBCDFRXykmmznAnJCKqilp2sT25wzjun38zM4EAjOW1W4wWmSd3ONci8+6bjFahkV/nbhUa+fXpek5pajau965PGYlTvkTG6Q/ylR8xjvPI3xrX3aItczeRAoxzPW+te8fzagzUsNgyRTRb1LqLja1Cmdi1RFTX2M1nw2SK6oIXkggvjBUC2L1FRK5gMmXDZIqqqlYfnl64e8sLCZ3FK61CRFS3OGaKyAm1Hq8EGK0wuVjrn9zhXDLllbFCgHG9L7yarUJE5AlMpohKYe9a+vNvZo5X2vUp57uWzphn/Cw28Nqq5xQOeiYimobdfETFeKF7ywuzTdtxrBAR1blyuvmaChUSEYBnHzS69tZ+Pv+UAKmDRj2ntC0FruoHnvuRkby98IRx99wLTxjLz/0IuCriTiIFnG4Vujho/GQiRUQNjN18RMV4YbwS4K0uNiIiSmMyRbNHrbqWvDJeCeDAayIiD2IyRbNDLe+ku/JW4NAvgce25R6v9NhtQMubjHpu4MBrIiJP4Zgp8j7rTrrOtxrzLPUdNn52vtVYv/9+Z4/vtfFKRETkKbybj7zNS3excaJIIqKGwUk7qX7s/gpw/Aiw4Tt57qT7vDEYe/dXgI/+H2dj4XglIiLKgckUedvC842fxe6ks+o5jeOViIgoC8dMkbcdfc74ad0xl81ab9UjIiJyGZMp8rZ1XwJalhh30k1NZZbZ76Rb96XaxEdERA2PyRSVZmoS+M1jwNODxs+pSXeOyzvpiIjI43g3HxVXyzmeCsbAO+mIiMgZ5dzNx2SKCrPmeDr/SuPOuc6LjHFKj90GPPewu48x4cN1iYjIJUymbJhMVcBLczwRERG5yFPzTCmles1fuwDoIhLOKg/ZFjUR6S+nnBzkpTmeiIiIPMrRAehKqYiI9Juv9QB8SqmYrTwEI0EaEJEBALpSKlJqOTnMa3M8EREReZBjyZRSSgMQMH9atgIIKqV85nIYwKBVKCKDAOwtUcXKyUmc44mIiKgop6dG8Jkvi26tN5Msn4joWdtoSil/sXJHoqVMnOOJiIioKMeSKRFJicgCEUnaVluJlY7MJMsuhelJWK7ynJRSc5RSrdYLQEtZgdNpnOOJiIioKLefzdcDICEiulIqkKfOKIB2GElTofJ8+gDcMuMIKdPKjxjTHzzyt8Zgc4u2zN1pEYiIiDzKtWTK7JoLAFjt8KG2ArjdttwC4JDDx6xvKz8CXHg153giIiLKoeRkyryzrrtoRSCcY5wTAEQArBaRlLk8mmf7drOsWHlOInIKwClrWSlVLF4qRVMzcN7aWkdBRETkOSUnU+bUBAMzOYhSKgqgx5ZIAeZgdKWUlrVeM8uKlRMRERHVnOMPOjZbtCJWa5VSyqeU8psJko4c459EJFms3NGgiYiIiErk9KSdQRgtST6lVMBcDuN0y1IEQNBWP2SWo8TyxjA1CfzmMeDpQePn1GStIyIiIiKTY8/mM+eJGstVJiLKVq8Xxp17GoCOHI+bKVheQhyz+9l8++837qRL/e70Ou1c4Ip/4J10REREDuGDjm1mdTK1/35g16eA8680noPXeZEx6/hjtwHPPcypCYiIiBzCZMpm1iZT44eBuy8HllwCfOz/Zj5oeGoKuPcvgCNPATckOGkmERFRlZWTTDk+AJ1maPdXgONHgLU3ZyZSgLG89vPA8d8b9YiIiKhmmEx51cLzjZ+dF+Uut9Zb9YiIiKgmmEx51dHnjJ8jv85dbq236hEREVFNMJnyqnVfAlqWAI9tM8ZI2U1NGYPQW95k1CMiIqKaYTLlVW1Lgav6ged+BNz7ceCFJ4BTx42f937cWH9VhIPPiYiIaox383ldznmmlgFXfJXTIhARETmknLv5Sn42H9XIyo8AF14NHHwcOPESMP8sYNl7jAcPExERUc0xmZoNmpqB89bWOgoiIiLKgWOmiIiIiCrAZIqIiIioAkymiIiIiCrAZIqIiIioAkymiIiIiCrAZIqIiIioAkymiIiIiCrAZIqIiIioAkymiIiIiCrAZIqIiIioAkymiIiIiCrAZIqIiIioAkymiIiIiCpwhpsHU0rFRaQ7a13ItqiJSH855URERES15FrLlFIqCCCQtS4EI0EaEJEBALpSKlJqOREREVGtuZJMKaU0AL4cRWEAg9aCiAwCCJVRTkRERFRTbrVMbQAwYF9hJVgiomfV1ZRS/mLljkVKREREVAbHx0yZic/eHEW5WqoAIFWgzF6erCiwUk1NAgcfB068BMw/C1j2HqCp2ZVDExERkfe5MQB9jYgMmC1Ndu156o+aZaki5TkppeYAmGNb1VJamDnsvx945G+B1O9Or9POBa74B2DlR2a8WyIiIqofjnbzKaWC5sBxN/UBGLe9Ds1oL/vvB3Z9Cuh8K3B9Aug7bPzsfKuxfv/9VQyZiIiIZquSW6bMO+u6i1YEwiKimy1RqQL1RvOsbzfLipXnsxXA7bblFpSbUI0fBh7qBc7/IPCxfwOazJzznC5j+d6/AB4KA0tXA21Ly9o1ERER1ZeSkymzhamcVqYNAJbbBosvBwClVC8AHUDCXNZEJGXbTjPL9SLl+eI8BeCUtayUKiNk0+6vAMePABu+czqRsjQ1AWs/D9zTbdT76P8pf/9ERERUNxwbM5XdvaeU8gEI2SfdVErpyDE+SkSSpZQ7ZuH5xs/Oi3KXW+utekRERNSw3HycjJZjXQRA0FowuxLDZZQ74+hzxs+RX+cut9Zb9YiIiKhhuTVpZwhGYgSlVEwpFQBOt14ppUJm999ye8tVsXLHrPsS0LIEeGwbMDWVWTY1BTx2G9DyJqMeERERNTQlIrWOwVFKqVYA4+Pj42htbS19Q+tuvvOvBNZuMbr2Rn4NPHY78NzDwIZ/5fQIREREdWpiYgJtbW0A0CYiE4XqMpkqJOc8U8uAK77KRIqIiKiOlZNMuTFp5+y18iPAhVdzBnQiIiLKi8lUMU3NwHlrax0FEREReZSbd/MRERER1R0mU0REREQVYDJFREREVAEmU0REREQVYDJFREREVAEmU0REREQVYDJFREREVAEmU0REREQVaJhJOycmCs4ET0RERJRWTt7QCM/mWwrgUK3jICIiolnpbBE5XKhCIyRTCsCbAByvdSwVaoGRFJ6N2X8uleB1MPA6GHgdDLwOBl4HA6+DoRrXoQXA76VIslT33XzmBSiYUc4GRk4IADhe7OnV9YzXwcDrYOB1MPA6GHgdDLwOhipdh5K24wB0IiIiogowmSIiIiKqAJOp2eMUgC+bPxsZr4OB18HA62DgdTDwOhh4HQyuXYe6H4BORERE5CS2TBERERFVgMkUERERUQXqfmoEotlMKRUA0CMi63OUhWyLmoj0uxeZu4pch17z1y4AuoiEXQ3ORYWuQ1a9uIh0uxSW64pdB/NvImUujorIoFuxuanE9wcNQAeArSKSci+6xsJkahZolDeGQhrtjUEp5QewEcb5+nKUh2BLoJRSQaVUpN4SiRKuQ8Y5K6ViSqlYsWRjtil2HbLqBgEEXAjLdaVcB6VUHEaCoZv1hwCoXHVnqxL+X/QCGLDeI5VSGoAIgB7XgnRRsS9UrnzxFBG+PPwCEAfgM3/3w5yHtJFeAHrN/wDWsgYgWuu4XDr3IIChHOuHrb8L27qxWsfr5nUw/w6Gsv42/AAk+9rUyyvf30PWNemt9/eJAv8vQgAiWev8tY63BtchXsq6enjl+PeOAYhl/U30Zl2zSLXj4JgpDzOz6aSI6AAgIkkAq2sbVU10i60Vyvy94LfzemZ+y/RZfxc2mvmNtZH4kPm3oNvWN6INAAZqHUQNRWB8AU0z3zcbTbuttaZume+FAfOnZSuAoFLKeg8IA0j35ojRs2NvqaoKJlPexjcGQ0O8MZQhX6KQKlBWd0QkJSILsv5PWOefnWjWPTOR3lvrOGrF/EDVYHypCJmvSG2jqpkwgIhSKq6U0szrUJddfCjwhcrNL55MpjyKbwwZGumNoRTtedaPFihrFD0AEjnePBvBmgb9smWxPlDbRWRARAYAxJVSsVoGVQsikgDQDWPs3BiAPfX4f6KEL1SuffFkMuVdfGMwNcobA1XG/KYZAFBXg89LoZQKmu8Rjcz6IpFunTPfO+xdPg3BPF8/gAUwun1jWYOw65n9C5VrXzyZTHkX3xhMDf7GkMtonvXtBcoaQQTAaqnjuzxzMVuxUzUOwwv0rJ+WFIz3j0YSEZF+s+WmB8aX0Wi9f3bU8gsVp0bwrmJvDI3UMhOR07e695itc3GlVKN25+iA8SGalThoaKy/izSlVBTG7fCpWsdSAxsALLeNAVkOpG8X16VBplIRYyoEwGjVt3f7aDUJqEbMv4OM9wERSSil+mEkGvXcgpn9hcq1L55MpjyKbwyGBn9jyElEUkopqwk7lVXWcGNmzFbKiJVYm9++tUa5Ftnde+b5h6SOJ3EtIInc3TcN8bdQxDDq+MtWni9Urn3xZDeft/GNIb+6fmOwydevH4ExXwqAdEJRVxN2Zsl5HcwJKjUYd+4EzOUw6vdvo5RxHprTQXhAvusQhq2Lx/x/MVjHLdjTroP5JcKfNV0AYLTYJFyJymW5vlAppfxmApVz7FS1v2wpcxIr8iDzUQHrzT5v6w+mW+psdudizBmN19u/WSilotZ1qUe2GY6DMFonB2BMzjdgq2PNjK8B6JA6m/0cKHwdzA+LsVzbiUi9znid9+/BrBeCkUwEYMytE62nD9AS/1+EYHZ1AkCj/b8wyzUAfWb1Y6jjp0aYX6DsPTgajDFiYbMVP/tpERnLVYuDyZS3NcIbQzGN9MZARESlKfULlRtfPJlMEREREVWAY6aIiIiIKsBkioiIiKgCTKaIiIiIKsBkioiIiKgCTKaIiIiIKsBkioiIiKgCTKaIiHIwZ1Gu6wfDElF1MJkiIsqt4SbIJaKZYTJFRJSbr46f6UZEVcRkiogoi/nsMz5QnIhKwmSKiGi6HgDRWgdBRLMDkykiounYxUdEJWMyRURkwy4+IirXGbUOgIhopszEx2e+AGAQQNAqF5H+Gey2B0Akz/GCALoAHAOgA9gIYKuIMPkiamBMpohoVjLngPKJyKC5PAZguYj0KKWiANYAmEkylbOLTykVArBeRLpty0EAm2Z6DkRUH5hMEdFsFRCRAduyBiBu/j6jOaLydfGZiVsUwALbah1ASkRSMzkWEdUPJlNENFvtsn6xzVSeAIDsBMcsD8JIgHwABvMMMM/XxRc1t7Hv128dj4gaG5MpIpqVshKbAIBkgVaimIisBgCllAbgxwBW56iX7y6+AIxEy64bp1vCiKiB8W4+IqoH3QD25iowu+40a9lMuLTs5+4ppQLI38WHHPsPgC1TRAQmU0Q0S2UlQwEAQ7ayoK1sDYBU1uYpGN10dutReKLOdIuVmXhBRJJKKT8fiEzU2JhMEdGsYyZLw0opzfx91HxZ3XjttuqaVZalPWs5Zxefuc4aa2XtvwenE7QAJ/gkamwcM0VEs1ESwACADTASpW4AYaVUOwBk3eWXwvTECbAlWPm6+GzWA+hRSg3DuINvvVIqppTqLbIdETUAJSK1joGIyDHmmKmYiCy3rRsG0G21KJnzUkXYwkREM8FuPiKqa9mzk5vddKmsxInP4iOiGWM3HxE1gvVKqQiAPTAeB7PeKjC7+DjFARHNGLv5iKihsYuPiCrFbj4ianTtTKSIqBJsmSIiIiKqAFumiIiIiCrAZIqIiIioAkymiIiIiCrAZIqIiIioAkymiIiIiCrAZIqIiIioAkymiIiIiCrAZIqIiIioAv8fyZT9PIXTx3EAAAAASUVORK5CYII=\n", "text/plain": [ "
" ] }, "metadata": { "needs_background": "light" }, "output_type": "display_data" } ], "source": [ "symmetrised_correlator.show([5, 20], comp=[first_derivative, second_derivative], y_range=[-500, 1300])" ] }, { "cell_type": "code", "execution_count": null, "id": "ff177781", "metadata": {}, "outputs": [], "source": [] } ], "metadata": { "kernelspec": { "display_name": "Python 3", "language": "python", "name": "python3" }, "language_info": { "codemirror_mode": { "name": "ipython", "version": 3 }, "file_extension": ".py", "mimetype": "text/x-python", "name": "python", "nbconvert_exporter": "python", "pygments_lexer": "ipython3", "version": "3.6.9" } }, "nbformat": 4, "nbformat_minor": 5 }