{
 "cells": [
  {
   "cell_type": "markdown",
   "id": "8ad2e16a",
   "metadata": {},
   "source": [
    "# LHE"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 2,
   "id": "8564ec16",
   "metadata": {},
   "outputs": [],
   "source": [
    "from pyfeyn2.render.all import AllRender\n",
    "import pylhe\n",
    "from feynml.interface.lhe import lhe_event_to_feynman\n",
    "from pyfeyn2.auto.position import feynman_adjust_points\n",
    "from feynml.feynml import FeynML\n",
    "import smpl_io.io as io\n",
    "from pyfeyn2.interface.dot import feynman_to_dot,dot_to_positions"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 5,
   "id": "2fcfa7c3",
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "<LesHouchesEvents version=\"3.0\">\n",
      "<!--\n",
      "file generated with POWHEG-BOX-V2\n",
      "Input file powheg.input contained:\n",
      "vdecaymode  1     ! Z decay products (default 2): 1 for electronic\n",
      "\n",
      "numevts 20000\t! number of events to be generated\n",
      "ih1   1           ! hadron 1 (1 for protons, -1 for antiprotons)\n",
      "ih2   1           ! hadron 2 (1 for protons, -1 for antiprotons)\n",
      "ebeam1 2510d0      ! energy of beam 1\n",
      "ebeam2 2510d0      ! energy of beam 2\n",
      "\n",
      "\n",
      "bornktmin 0.25     !  (default 0d0) generation cut. Minimum kt in underlying Born\n",
      "##bornsuppfact 10d0 ! (default 0d0) mass param for Born suppression factor. If < 0 suppfact = 1\n",
      "\n",
      "\n",
      "! To be set only if using LHA pdfs\n",
      "lhans1 27100 ! pdf set for hadron 1 (LHA numbering)\n",
      "lhans2 27100 ! pdf set for hadron 2 (LHA numbering)\t\n",
      "\n",
      "! Parameters to allow or not the use of stored data\n",
      "use-old-grid    1 ! if 1 use old grid if file pwggrids.dat is present (<> 1 regenerate)\n",
      "use-old-ubound  1 ! if 1 use norm of upper bounding function stored in pwgubound.dat, if present; <>\n",
      "\n",
      "ncall1 10000    ! number of calls for initializing the integration grid\n",
      "itmx1    2      ! number of iterations for initializing the integration grid\n",
      "ncall2 10000     ! number of calls for computing the integral and finding upper bounding envelope\n",
      "itmx2    2      ! number of iterations for computing the integral and finding upper bouding envelope\n",
      "foldcsi   1     ! number of folds on csi integration\n",
      "foldy    1     ! number of folds on  y  integration\n",
      "foldphi   1     ! number of folds on phi integration\n",
      "nubound 30000  ! number of calls to set up the upper bounding norms for radiation\n",
      "\n",
      "! OPTIONAL PARAMETERS\n",
      "\n",
      "min_Z_mass 0.5\n",
      "max_Z_mass 10\n",
      "##runningscale 0\n",
      "\n",
      "alphas_from_pdf 1\n",
      "pdf_cutoff_fact 0.01d0\n",
      "#bornonly 1\n",
      "#LOevents 1\n",
      "\n",
      "storemintupb 1\n",
      "\n",
      "!! parallel\n",
      "iseed 10103        ! initialize random number sequence\n",
      "manyseeds 1\n",
      "parallelstage 4\n",
      "xgriditeration 2\n",
      "maxseeds 2000\n",
      "clobberlhe 1\n",
      "\n",
      "!! shower/analysis param\n",
      "NoPythia 0\n",
      "NoRivet 0\n",
      "QEDShower 0\n",
      "QCDShower 1\n",
      "hadronization 1\n",
      "mpi 0\n",
      "pythiaHepMC 1\n",
      "HepMCEventOutput 0\n",
      "HepMCEventInput 0\n",
      "rivetAnalysis \"ALICE_2020_I1797621\"\n",
      "rivetAnalysis1 \"ALICE_2020_I1797621:PHOTON_Z=1\"\n",
      "rivetAnalysis2 \"POWHEG_2023_DILEPTONS\"\n",
      "rivetWeight 0\n",
      "End of powheg.input content\n",
      "PDF package: lha\n",
      " Random number generator initialized with:         1000             0             0\n",
      "-->\n",
      "<init>\n",
      "     2212     2212  2.51000E+03  2.51000E+03     -1     -1     -1     -1     -4      1\n",
      "  4.37712E+06  1.47821E+04  1.00000E+00  10001\n",
      "</init>\n",
      "<event>\n",
      "      6  10001  1.68713E+07  8.94427E-01 -1.00000E+00  6.65838E-01\n",
      "       2    -1     0     0   501     0  0.000000000E+00  0.000000000E+00  1.749324861E+02  1.749324861E+02  0.000000000E+00  0.00000E+00  9.000E+00\n",
      "      21    -1     0     0   502   501  0.000000000E+00  0.000000000E+00 -1.592747830E-03  1.592747830E-03  0.000000000E+00  0.00000E+00  9.000E+00\n",
      "      23     2     1     2     0     0 -2.408277128E-01  1.052626414E-01  1.565325189E+02  1.565338272E+02  5.835294129E-01  0.00000E+00  9.000E+00\n",
      "      11     1     3     3     0     0  7.521133062E-02  1.659522570E-01  1.774845055E+01  1.774938573E+01  5.110000000E-04  0.00000E+00  9.000E+00\n",
      "     -11     1     3     3     0     0 -3.160390434E-01 -6.068961564E-02  1.387840684E+02  1.387844415E+02  5.110000000E-04  0.00000E+00  9.000E+00\n",
      "       2     1     1     2   502     0  2.408277128E-01 -1.052626414E-01  1.839837441E+01  1.840025160E+01  0.000000000E+00  0.00000E+00  9.000E+00\n",
      "#rwgt            1          24   6209355.4982931120             1000         419           0\n",
      "</event>\n",
      "<event>\n",
      "      6  10001 -1.68713E+07  8.94427E-01 -1.00000E+00  6.63949E-01\n",
      "      -4    -1     0     0     0   502  0.000000000E+00  0.000000000E+00  9.983585819E-01  9.983585819E-01  0.000000000E+00  0.00000E+00  9.000E+00\n",
      "      21    -1     0     0   502   501  0.000000000E+00  0.000000000E+00 -7.764707284E-01  7.764707284E-01  0.000000000E+00  0.00000E+00  9.000E+00\n",
      "      23     2     1     2     0     0 -2.586998664E-01 -1.510185871E-01  6.026979797E-01  1.290320707E+00  1.100886206E+00  0.00000E+00  9.000E+00\n",
      "      11     1     3     3     0     0 -6.211468495E-01 -2.712056290E-01  1.429668701E-01  6.926872952E-01  5.110000000E-04  0.00000E+00  9.000E+00\n",
      "     -11     1     3     3     0     0  3.624469831E-01  1.201870419E-01  4.597311097E-01  5.976334118E-01  5.110000000E-04  0.00000E+00  9.000E+00\n",
      "      -4     1     1     2     0   501  2.586998664E-01  1.510185871E-01 -3.808101262E-01  4.845086034E-01  0.000000000E+00  0.00000E+00  9.000E+00\n",
      "#rwgt            1           3  -1231479.7875413459             1000        4036           0\n",
      "</event>\n",
      "</LesHouchesEvents>\n",
      " #Random number generator exit values:         1000      58113976             0\n",
      "\n"
     ]
    }
   ],
   "source": [
    "print(io.read(\"example.lhe\"))"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 8,
   "id": "793907e4",
   "metadata": {
    "scrolled": false
   },
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "<?xml version=\"1.0\" encoding=\"UTF-8\"?>\n",
      "<feynml version=\"0.1.7.22\">\n",
      "  <head>\n",
      "    <meta name=\"feynml\" content=\"0.1.7.22\"/>\n",
      "    <style/>\n",
      "  </head>\n",
      "  <diagram id=\"FeynmanDiagram9\" style=\"layout: dot\" class=\"test\" sheet=\"\">\n",
      "    <propagator id=\"Leg3\" pdgid=\"23\" type=\"boson\" style=\"\">\n",
      "      <name>Z0</name>\n",
      "      <source>Vertex0</source>\n",
      "      <target>Vertex3</target>\n",
      "    </propagator>\n",
      "    <vertex id=\"Vertex0\" style=\"\"/>\n",
      "    <vertex id=\"Vertex3\" style=\"\"/>\n",
      "    <leg id=\"Leg1\" pdgid=\"2\" type=\"fermion\" style=\"color: red\">\n",
      "      <name>u</name>\n",
      "      <target>Vertex0</target>\n",
      "      <sense>incoming</sense>\n",
      "    </leg>\n",
      "    <leg id=\"Leg2\" pdgid=\"21\" type=\"gluon\" style=\"color: green\">\n",
      "      <name>g</name>\n",
      "      <target>Vertex0</target>\n",
      "      <sense>incoming</sense>\n",
      "    </leg>\n",
      "    <leg id=\"Leg6\" pdgid=\"2\" type=\"fermion\" style=\"\">\n",
      "      <name>u</name>\n",
      "      <target>Vertex0</target>\n",
      "      <sense>outgoing</sense>\n",
      "    </leg>\n",
      "    <leg id=\"Leg4\" pdgid=\"11\" type=\"fermion\" style=\"\">\n",
      "      <name>e-</name>\n",
      "      <target>Vertex3</target>\n",
      "      <sense>outgoing</sense>\n",
      "    </leg>\n",
      "    <leg id=\"Leg5\" pdgid=\"-11\" type=\"anti fermion\" style=\"\">\n",
      "      <name>e+</name>\n",
      "      <target>Vertex3</target>\n",
      "      <sense>outgoing</sense>\n",
      "    </leg>\n",
      "  </diagram>\n",
      "</feynml>\n",
      "\n",
      "tikz:\n"
     ]
    },
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAALUAAABBCAMAAACgjhWUAAAAIGNIUk0AAHomAACAhAAA+gAAAIDoAAB1MAAA6mAAADqYAAAXcJy6UTwAAADeUExURf///wAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAP8AAP8AAP8AAP8AAP8AAAAAAAAAAAAAAP8AAP8AAP8AAP8AAP8AAP8AAP8AAP8AAAAAAP8AAP8AAP8AAP8AAP8AAP8AAP8AAAD/AAD/AAD/AAD/AJhoAGScADJWAAAAAAAAAAD/AAD/AAD/AAD/AAD/AAD/AAD/AAD/AAD/AAD/AAAAAAAAAAAAAAAAAAAAAAAAAAAAAP8AAAD/AP///ysOM5AAAABGdFJOUwAimRGI3WZEu6pVYKeEd+7NzzNpxLvMiFUR48jx3UQiZpmqdzPM7sfWyEC/wBGq3Zmn3/HsdYhVuzNEZiLMd+7Rn/rHklxQ/zUvAAAAAWJLR0QAiAUdSAAAAAlwSFlzAAAAZAAAAGQAD5bF3QAAAAd0SU1FB+cFGxMaFtL0jNcAAAPwSURBVGje7Zp7e6JGFMYHZoCBttjSGt3dxlztHWd62eEmuLfifv9PtOeM2k0bY0BRnOfp+weJiZKfb95zzjBIiHmy7L4J2osyx+2bYQ9m1jfDPsy0b4i2co1ktq2+GdrKM5GZcxOZvb4Z2sryDWS2naBvhj2YzRuEnmtYNjzviy+/8swZKiHngxr1tTnM4HK90jcmQRPCNfQg6pujnb7V1GbVYRQE3wG0UQ0PmCNMtkkXWZbn6jTbvG+SVszrMFNj2od5g9BcZmNCsdH/zCdj3nJFaIV9c+1SFGxjDvy67ptsF3PweIlEQweXIX2zPZYbYpCtbcyRvV6n7j7D8GI0fnFaaGYzJyL05Sv+SK++X0PXfJcIeTEeXQ5HkxPQ0pBhiKlDSeQwtn1v1LObeK11dX1zeXt3f3SPXdwl0LiU8afaHWVO81zf390C+dXnH4y6eheU2Rhg6kCWpyxwnu/OAW9Tjfej4eXN9Zr8h+VdN9TMDnBXg+keHLJGV1Ztd/UmowtdoJPlcnl7UJ1S28cVkQUeU+67DTw+SFCgFxc/AvXy8pCU2CyaQoaZvphy2Sk2Ra9+WmrtkRLqc9ya0x7bnB3b439pudbPv7R9pR1a3KbE1h57J70pMV4xD4e//va5QHfLmg6wqaHHhDlsekLaja5vRuPN6JlAa2kwQX0X6m/jsXUOy09doKPx1t9ZjvOPx6Fj9+HxLnKYoNtmP9Qem1pgtX4L53j593D2R84Uio5YuHkROPx89zBiPODsxwLlHgQiWntMz9HjlWZCrr+b/P7Hn38x5r12pn7fVE8rVgkcBUkzfJSnQhZzy3OZzQecA/3ZOa1TIVWZI3Wh9KPk4RMiL2T+mdHPhQB/0zifQULyHKljsfWZjemP28kTCT5XKk/B3YLAMRezgqDxO1/2LP3gWB/xmBO0tCoJZgKOWZkBNck+qnwuZdzoHE/R07o+zl1EJZA0JTKHhKDLpFiAwXGpKplWLU/2X3p9u6jbu8yxRKgSIwHFp+CBkhhjiHcqRXHIqTf0b1bXsx19ciLGf3ylMLpQfBLyWynEnIlknoiExHmzZDynt+vNg3fvO+g5sUArIbcQDlXEKQZiAZyJKFQqs06AtdgK2v/w92EdM0GoXGFHqyoIBuRBYZ+oIBazct4dsBav6yl7gLlvvxeFiAlkYpaD47obZwsJzaKMSZZ3zAyrWXd7w25LnxIMsIhVjs1CYk1KaBaqwUuPoKb0IsaZp1IxxzCLUoqegFvRK6nXF0oUWV6C63nXUe6aXveGWC7gayag7c2SQ/7AiehJqQdHssAZ0mF3Oy49xBiaRQXdbd59szimYIaobobeyfQJO417BfcmDG0AAAAldEVYdGRhdGU6Y3JlYXRlADIwMjMtMDUtMjdUMTk6MjY6MjIrMDA6MDC4+Rr2AAAAJXRFWHRkYXRlOm1vZGlmeQAyMDIzLTA1LTI3VDE5OjI2OjIyKzAwOjAwyaSiSgAAACh0RVh0ZGF0ZTp0aW1lc3RhbXAAMjAyMy0wNS0yN1QxOToyNjoyMiswMDowMJ6xg5UAAAATdEVYdHBkZjpWZXJzaW9uAFBERi0xLjVTNlrBAAAASnRFWHRzaWduYXR1cmUAYjRkMThmZDg1OGJlY2NjMTEzMTdiZmE4YjNhOTlmNGE5YTY1M2E2YzFjODNmYjZhYmY2YjlhODhkZTJmMTMxYrUfnS8AAAAASUVORK5CYII=",
      "text/plain": [
       "<wand.image.Image: b4d18fd 'PDF' (181x65)>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "pyx:\n"
     ]
    },
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAWkAAACBCAQAAACoEoOeAAAAIGNIUk0AAHomAACAhAAA+gAAAIDoAAB1MAAA6mAAADqYAAAXcJy6UTwAAAACYktHRAD/h4/MvwAAAAlwSFlzAAAAyAAAAMgAY/rnrQAAAAd0SU1FB+cFGxMaFtL0jNcAABHESURBVHja7Z2/q/xcXsdfn+ULC6JgttpS8gV1OyXfwu5pzgVxGxFytbAQhVzczmZzC7ezyGAhiogz/8DKDGihhTgp3EqLGZ4VLVxkAouw2uyN8MD6WB2LnMnNr8mPmWQyc+95pbkzyWQyc99z8jmfX0c0FsvbQEJCnA9zX4bFMgYSEOICKytpy92TyzlmoWMractdI4oIDyNnACtpy90iihAF7HnO5AxW0pY7JZdzwkKvSnusx8Nyb4iLa+ScvI7OR+wobbkrcjmn7OtyBitpyx0hDt5Rzux12nyUlbTlLhAHDw+nXc5gJW25A3I5Q0yik/ajraQtN46o/nIGK2nLTSMKFxfYs+8jZ7CSttwsBTknet//dVbSlhtEPFw8BssZrKQtN0cu54T9UDmDlbTlphAXz8g5aQ6kdGMlbbkR2sPc/bGSttwAxTB3eyClGytpy8z0C3P3x0raMiOluOAIcgYractsDAtz98dK2jILeSBlVDmDlbRlBs4Jc/fHStpyVcTDOyfM3R8racvVuCTM3R8ractVKIS5LwqkdGMlbZmcMcLc/bGStkxKd/nr2FhJWyZjzDB3f6ykLZMgv8QDP8+PrytnsJK2XIxkPuYsDugBDp7Z9c9857pyBitpywVISNSy+3v87viBlB5XZRuIWc5HAiKcxl3/ymfXHp/NNVlJWy5BPJa5ofFKyqc5RmiAr8z8jVjuHL3ngVXt6Ye5BG0lbbkcl5i/Lz3zNF24uxs7PbRcQB7m/ht+NX+y0u/56tdkbWnLeZTLX3MhbfTjzNdlJW0ZTj3MbYS052EeP0fh2qykLcNoLn8Vzax+jlesLW0ZQEeX5xn9HK9YSVt60ln+Oquf4xUraUsvOrs8z+znKFyptaUtXUxb/jo2dpS2tHJul+f5sJK2nOQ65a9jYyVtaeTSLs/zYSVtqXHd8texsZK+YUQR4RGTsLjWtGysLs/zYT0eN4sELPMHKQ/T3/6vX809yaewkr5NxOUAPOuFOEQEwNOUnt+xuzzPhzU8bpUAWOgF6JQnSYhYCtOIeoouz/NhR+kbRQ64fHy1oI0Z8mls82OqLs/zYUfpm0R8XOKiwPRKHCK2MqpNPV2X5/mwkr5NPGBTfkovxCVgKSNlJM8Z5hYHn4DnKSah1vC4SWSHVzQ78ue3KFb66eLzT97lueW9XQICHGL9MMn5raRvD3F4IdEfG/cccHjUm+Fnzc8xY5hbfAKUeXDRpziNNTxuEQ9ovCXrVJ5Ys5QzTYX5wtziEuDj5k8k0wjaSno0JMAnJWE1gl2aRe8a0RtZELJk8E17vjC3KHyCypOLyd6NkPi+0lJuD3FY57fTEUIisiTg4ZTsxGGHy7MeJAr5O74JfMk/8i9X/Gq+yi/yK/xs7fn/KURGR+YDDr54JPB23DhXJ0SR8EyKImQprn6+6HwuJ0dp0Kk8siOUTf//lwR8gx/wOf/N/13xe/kqv8XXG/f8+3RvKijIR5iEBO4/fnRdzJTta9m3ZhofXuSVkBccLa1HRIRDPAYSEfLH/MMMX88v8xnfrD37Bb8x1RuKxkSQin2B96Rwv4kr10VCoqKExWOLc4n5IZq9/tR6xEDjo92Umfwbqk4NYcKMlZITz0jbzd88BivtLmSHVw5Ui2LL2f808dh1j8Ci2PbvmyFbVJOf+6rfU9GBR7OTcpT3afJLl/q6Z7lZkNpJZBPi8FIfUyVgeW5CqCi2LLqt8SHGh+zw2k2ZK31bxzALTDdO65YND4UiJCIiIkShcDX9NlTfI+95I0ATNTwfodnhnHHGsPmMteMcDmiCXufUvMz9TRWuO2CH5jDN+XtFDyW7XQycRMqBzYUz/ztAlgTNcTDZovqMtrXXRYT9ImvG+PjYPZ0XPVX4+VzEI2A/xTg9KCAuitok8rS0ZY1PzOPb9p/IDu/o7ajscdmdE7yWNX7fqZxEhN29QkWxnb+j6LU4I8ej7yRSIkIg5fG2JpjiEOCQMsCv23I2fXqiIz7rfqNo6VVbVPOPpPGz7HC7fjbisz7nfnGnXGATuSgUgbG0IxQKr7Bf5X/2sAyvZskpDuNdFwrNtmX/Gs164DkP6IFX8NJusxPR0+Z+C9sFOR46s6mR44p3CkAUe7JI5OtsPxTvNgwQcVnjELPBJSAU58JETQdo82o84eGLP8j4cBlw99CxbPBZ0mZWODDknHfOiONf5h85jtkBih8Vdh+KI/hcGxH52IzHC5rdxecLW4/oMYpWXtE67jcc7/BC6yjMFt3fU3Xv2xSnVCh8IiL+q7IrHP/dBl7bS/Gfi8sOzfKC80XoLmclEQOMjy5T5uRrXk6Llpchpsy9b9Od2OG7DU+vz/HUjnZNARVx4fLCBTZ1n/HP+I/9nmdU51wPESd/CLgM/pHc8zbdiYMTO37E789lgrCkJi1jfpw5eeo3lcOnt/HRN9BSeZXDjhN3wVPBoLe6TXVa1br7uyYSeWVpc0DXhWUEd9a10DMG1t/z0ceUaXydh0Y3fYqmH/Jb3iapPTS5aEViUhJS9vwUP6EaibxKpra4HJoz3CQiPGctKHF46ReV618zOCTQUnllSNSUVyIH3L5+7rfABIVa4rAGYiPhtDG+GOdBdhdArpGprThV0fcsHoqQocEID+h1xXnNYNz5Cbvcgqffo7Epgijc95UBP4GkdUqPtEGzWt4+z9R2ASVTZmpnRf7NPLEjlPiMd+55f9Eb2eD3+Nl4cK4A9ZO4KLYlUfvAjayich1upOnBNTK123OGJWBJwqchcpKIsH+g2RgfHUbFZQlG4rDFezWixOFleEj+vrmRCnGdZiLOM7WzSKQ3aqa2B6dtdr0S1RmFa6L3zy43Plp+NqL6nu3kezywxeMgjzoWlzWweU+CvplRunZZxyB7sQihYxKZF/Wf6PDWVf7UbxQtvWLwVE7W+G11iWMkGOXNezMG3nneAHO7XDpcU6oWZFdN/l2TVF7cdmVfMx6dAQd8NIf+wSC26GGhIxO8PulS6w6w93yfkG0WRJoztDXPdqOjdJW2TG3xWJqWKxtS9ng4pngz4fnoNutX/iRr/AHW8QF3aPlTe8Lp+S48y5E7kbS52KZJ5K/xOzjELMpCEEWIguNtXiLC7orqYcZHW650y6uWBKemgN3tDixdfGXuCxiEi4PC5T/5ITEJit/jD3D4KxZV/7CO9QOPpASyk2PYp3OaqVOegaU4XUeCnJuy+UyCkrDhjB5O/+mmpZkb8Xh0IQ5hXlmckRDzfb4F/DU/RJUytQ16I3vWeGzlITdaOtAr8XuGXXoHWirvkHVLiqTe6FBxZpjF8spdGB4SGTnH+T9c5RJd6SfjH2kMshtP7QoX1e+WbmoGO42PS7wTzU0RZIuylvTFzD0/7TF7X6LRLMtJnDjm+UJpQSFTOyLCz/wjxsvwRf+cYUJ6eD4u804QUclxxmWyRgDvabsxw+M11HAcq2RJ0JSMo1PJsiE8tvKU3cJNkD0utUNT7PlzvsNP8/3eP/OFqB7Gx9nZGAD6WRwCtvKYf7KA2mIWljOY+zeVj1HKjLqv25rAjGVNKZNmTCNC05iOiUvAslA8+2NT8qu6kzdNaUBrwunl5U+s0bxkXmpUueLGbuduN2FL5/2ZU2LjQ/DyyOGJhdEkYJlZsl3NusTF5et8m5/wRWnH8fi4sWVDSNQeeRvD4SZLAjI/u4J31JhgQm5A0ia7OmFR7r2TLzjc+I8u9uLM8p3be31CJYOink1RkLWOu7sldfcX7fXpfZY4QMrzlKvPvh+uLGnJI4C5rZwJesVzdTwUlwMpnNhb6nIkO7zhY1yHxP+Nv+Bn+DafN47io3YwEo/knWViTMYVJS0BQck7vGfDhi1ucyJPFu8jziRfPUI0qf5a/shjd3kSZU3iv8Cv8yV/aTrnl0bxftFIy/W5kqTFJzJZF3sSwDVJ/8DJ4LBZctg0yaqIuppXLEuCsW1RUfwhn/E9/oj6KP6b/Bzf4j+yB9aXfDtcRdKmO14lD0M8Qnxgz1N9alesFDTGSalyryZpl0Nx3B7pyrOfk/mplEbxLV/yp6WDK7a4ZR5G90vn1nJe7yYRYdPkR+8lBZIsYF0Tdbay9fHIzso9nchmcKuuTnQqj2wJJdGrolDFB/62lI+hiuO4lOscbbv5KzKipKsrckhCzAaH8KSLzSflExFBg6hLN/rGyr1q1kaMjzd2sELv5ZklSymvFaiAcqViXLPFixL3TRlDxogSlxCHhZ1avjKS4SEOIVlu2Z49KU5eYQInltYVn3VmIZsIYWndEYkIi/kOxv/xseTjKJXy91vj5KxPl3m+86bCfSv6Wj0qI0ncfC/WAZgzyihtUoNSVsW1WcXBJzIFV01jZ7aONaCfBALWtHh5dSIrAoJ8VdM9Hl7Jft0LlVWbRkKvxCNgK8fVRSJg1T0yXmMU14ns8VhKyJO14OGSlul5EmdelfxYrfczI4gDTW66csV21bNcHaXNKJwHNyQkqno4pkygN4GfDQkKr/9aVq3nHGUUl5DI/BnzZJdjHbqwhY9fyFreEOuVbFHNvYoy2bE54VkuCbDqWW5qKFDsG9R08xcN09WE5HG+Ez6aEd7hLImLy6Fw3OK9W9a9JS0ekfmajyNn9ijFOdV8S7YoPum9ccJVBVoJKJc9y021guUFKWVJUA52nFc4NeDLMoVi17JbCxJ3apPhosT/jG+U9izedQioZ85YwAuaA0ExixifLRrNrjmDjDz/1/QH9St7SxXbuBSWMmtqIVvOUK52VcZn8CIS97ThHPMIS1nhERH/VDv48D6W6Gv8nk487ZWS07PGuqdW99M0JoCiKEis3lCAumTXFGTPgUqyZbWtbHltQba8oxVJShL/QeMBe37b7L+B9Reu+L1UvqTAjLvZtiXCNW1eG8XCAY3HkgZR41MW4LryWFfrTMqSLS5C0bRfo2GH5oUAxZp3WhOSDyvV7X/5nD/JR/KwMMa/aYkXbGkJjMvtaKW99jpymhN0Xn0QJz3LhdfVPMu1xdDE5fDqWRaXQzljuZ4qJA5LfPPgzAWO7xvx2NWejImJSQr2d7st/qaim8YvnTedWrEphDc8fFPI2jyH9jGlRcazvKQa6Ch8VTXPcgKoosdaJ1KICepE4kq5lEelNEqnPIoy5QM9PMVvkGXh74QN+0JYq+in3ohTkfhVopvXx4zSssNrWnTTtGqB055lU3NivNPFUbm2cHDFsxywrOXXlZZvN7XaJsAhHrtpPRr3h/FJp9m4PMwnXZN4mbuVuGiM2dDsWQ5YsmLDutGzvC02Eqh5lpuCIV2e5W25NYEJcKzY4LA8ZQC9V8QlImEzjuDeisRFt3ZpkzU+j3pz0rNcCm00eJaro3Dds/xU9PKWE/uhFOCgreenZWzuVeKis5HzxCoirxE+I+pyiLoiwHLOctPaKOWYoCi2pXG9MfFIHAI8UvY2NWdO7kXiQsDyZF1JSWL1PvkNY6oZ1c0jU5dS2F+JCZYX/qnHAy23y61K/AMBcEpEpeYrDd3iklrmWzlneUNI0NrgZYHCYyfPJPgEpO9rXZF75rhyg+FmPCpCS15EdYLX6Fkuj8INnuXi2etTxryHB0DK07hVKZb5mGsU/0Bb06pKS9qGnGWFV2xIqxNJX3+dmWdZooKEC+VX5piUBwnwcElY2NTIt8MFo/hF62B+YMACOsCKAP90sIRM5q9ki68lxrPsEpQu/fjhrbHxDhgk8WP/2YxBEv9AW6vC2on0XhI8cc1bbFjiy3NrkWtWueeyAtb0qgaxvAc6JV6Ypw2RuKBPJ803epZLMcEGz3Jt/ZK8ERhwzrLGlvdJhy1+UuLtknZ4qXmWSxO8mme54RUgLqFZG/adV1xYzqevxD+0WdI6LZkZGeW0oVg2+GzzQtisBU31PAk26me5kN6GSo9c3HKGsio/g8MOzQ5luufbHsl2m2XDMdU+HYfVFluoJvbnoj5u76iuxG63uHUsEqcTYlyJCk9li2aWbgj6E0/ExGx4sA45y7x0Voibqoljc/Isf/qjDYlYbpXOpTz1nidgK6GpfXGKHZUsllujVx8P00w3Y5Q+QxbLVPw/O+ukpF5GvIEAAAAldEVYdGRhdGU6Y3JlYXRlADIwMjMtMDUtMjdUMTk6MjY6MjIrMDA6MDC4+Rr2AAAAJXRFWHRkYXRlOm1vZGlmeQAyMDIzLTA1LTI3VDE5OjI2OjIyKzAwOjAwyaSiSgAAACh0RVh0ZGF0ZTp0aW1lc3RhbXAAMjAyMy0wNS0yN1QxOToyNjoyMiswMDowMJ6xg5UAAAA3dEVYdHBkZjpIaVJlc0JvdW5kaW5nQm94ADEzMC4wNjR4NDYuMzc1OSsxMS42Njg3NzYrNS4wNjIxMjVPeZBQAAAAE3RFWHRwZGY6VmVyc2lvbgBQREYtMS40JDFqVwAAAEp0RVh0c2lnbmF0dXJlADIzZTZmYTg0NGY0ODhmMjMzMDNlYWFiODAzNDI5NWY1YWVlNzAyNTdjYWNjM2JiN2M4NWIxYzg1YmQxYzU0YjbaaZAGAAAAAElFTkSuQmCC",
      "text/plain": [
       "<wand.image.Image: 23e6fa8 'PDF' (361x129)>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "feynmp:\n"
     ]
    },
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAALQAAABvCAQAAADvPSJVAAAAIGNIUk0AAHomAACAhAAA+gAAAIDoAAB1MAAA6mAAADqYAAAXcJy6UTwAAAACYktHRAD/h4/MvwAAAAlwSFlzAAAAZAAAAGQAD5bF3QAAAAd0SU1FB+cFGxMaGUJLkUYAAAnZSURBVHja5Z2/j/RGGcc/gxRRIIoJ/4EpELUpEEWqeRtEkcZ0SGnwtVAZUdEx19BGPoo0kSLtFamgwCslHc2uEE2KSD4hgQgIyadIkVJFQ+Hxz933bu0dj8f3zhS39q7tma8fP/N9nvmOD8O0isagpx71EisphvTiX088eYyhRK7dyRAqB8zlSHyLaUUDN+Zx4lEvsAhFzN0EJCbeRcNubUsKo1JgiC7//SSLFgp46G1LocVBGFGKdG0L81tEhGJvHiYcMukuZhhUuyU5YKgoqHjDPDd5H4mLjpj048E4S4HlH0g0hmzt7nuDWWI4TDxm0o8jDHkP9F3v0hXl2gB4A1ozgdjNANpASYW0wFZ9Z0HOBLKz5TrPqKbSuzskOaCQ3A7IjeLhDaF9CZL7yUdNvp8FBo1mMBiQ8sbEi81TPfGoyQdElBgqTG/PDkM1hVVut5Iwy6TmXEqyw2A4oMkpMBgKIgOSjB0l+csFfWqgcgXQBhvnVxhqHq0MBmK7p8TMu+vhVxQt7/IDtMYQE/f2SEoMGdJASjnvvodaUaTIOYHKtUArRgHKMH360gZHNIaK9zEUM88w87ARj67Tp4NfVHObFGK1o5LBjkbT61Qe3ZDCR244krU7Yvb8vPteRBw5rk14HZao/aQoRSHiyWe44i4XtO6DoqN77aM2MUgNuY42Dx54dO/QmlHvkENHYenfjMaEWokGm1WfBHgAuoW07AZCYjRVAzOSlJzdvHE6nIoabCazznF1I7KWPZd2V1k7jZZXGwzFHCsIpaJ7GzOTwW4aklBiKNiha75pIKLCkBERkTObFoVQyduPs4IVZ0Dbez6gPezoTxLkmO3atE0zXDXuuGrKYBrApl56NoxiwwGM/TMjZ9fVmTz6hCTecYvsscuYe24GP9lvlVeLhkO/uirf7uyuR1RdqpSRT2bHZnMfvIvh6qjAZYNSDBWJDci7vIckZ/JkZjiVjzB8cPVZnDYpocKQk3UW0OfVBhS6zvBtpSL5hM8dnMdxsyIO9uOBwjLr6oRXV9sZGIn5G+8GB7TBWm39sWDX2K+FWaNstnojmRAyfufiCVyqeeXYJ3OgTZlvRwVCjJ4Xco+rI3p3Uo7EQvUokibm1uwt03nknqj/fbAlBvYuTrQU0LdALmS7nfLAbe/7TShARETE0Y1aZSGgzZEbHpqJARFzHMltJPfGiaUsWmJwFWYtZdGYOySZ2AkJxKjOhoUUOWn4Ni0kMQ+TpLlPArLcQNJkq5N+ysmmVXdtji8JNWIkQbtLhC3d2Kz9qNGWWVd1TtcmVw0mRMENksylEHn5BkfoNogxFGhrywmGipyMggAFZSiX9uwBaIPNfQxS5lauYDtCRnBJVD5ymyrw1ezR0gtyBpNClFRrQzto768xfI1295z5avhQxxQzTqOO5AprV3IM/8RgyN1MLftsejtBS9III+22pFjPdZCRU1B0roIUw8FA2mplB5kZ5PTW+utOjM1Wj5dh2NVdqyw1QloomzyjRlu1d2MUsZ2aLZthvJm4m+bBfXYqsfk7SdF5ZNI6h91ue9WAUHSQtcKJE30dkc2qWyJqb8thCivxaz/9bHWCRltere2NaKzLUwr1zFgRv/5Gk9rWH9td1eWZPa9AW7spBjsaXp1aW6ptx4sKBM1E3RGqpyu17Q8U6LaDGWok+7V+0fpNDwHMnCnjkQ6vNo4LvPU6QJ+oPCjoL36evYDhiWvGpGjSDlgUM56cnmqpq+Xz3noVoO1w0jaOdAzsHAguhkej7HRbNXXoRfb0hP36Fb8KE2hFT4DAjsGq26FcwcHVCurhN24HNMNc+a0i6a1G69e/PuWGVgLa+umGWB36M4iWV7tLUCqGcjWFRvdHiNlnlihSNB/zhd31weuekfWAbrLV8fB1K3aO3KGH9qWS4reW+BXn6OlqQBvabHVJ5zfr9YsNr86tIOcqmDj4y6OciyMDANpAfMJL8wGvPtgo7PLAIBo6hSUYzLMt6OLIbsBfF2jbiAMGTdJBVA+Wllenl6pA2g62QbQN8FdQZttht0s1rA2yBWgsYx/Sv4vWq5JTP7J6lChyIoCZ1S/Vk+KvDbJt0FC2PnqDwCUydlJ6MRrKpj4dpu6v7OPaDWhh0oPIUI949bPzd77C9q0DXc8hJh1og+9Knn2NEIEvR1pMQDNRXfLIK2AntN3RylZEyoFopHM6FbskDN7IF2BZ+073bLLOVpd2MJOoNkM9nEhK2FFSsWufgFoSHLDjMAiz9p0eFJGikINd98M3JgpNBhx5RAFH9sQo4Mbcrd36J8vad/qMZTdC9kFS036XYjhYjhyxs7x5A+tyA7PopoiUhL25PdlfIvn+Jl/7tvadfq1dH2jD8XbfnwlQPHZhf9ZuwGsbVk9olWTEBmLe40sM32Aw44TNFmqgrqMuIqF5XfJ3eZu3+Au/53vERDxwNJtaiRs00CAkCsUP+Db/5mPzod2rUMDe1bIHLz0JG2j7SuyIhxHJi4iJeWS/FbsOHGgRo5Ac2Z/arlDESI4cnS1/WLIn4QItJOppq23t+hj+wqNggT7vMs78rrH5wO06UKCtW9hfYqlC2iD8ol+v1qPwgH7eZZw5Jg6d9AUHtIXs7PD35HGtXYdJ+gIDumHI85yA9etBkr6AgLYu49nh75mzBEr6ggH6KcY88UxBkr5AgBYK5fKRD8+uAwD6UsY88ay1IwqG9K0OtHUZi8AREulbdRZcSJHzIW8v9e4Oc+SePRGJUN1LWtqZdr9ltcR+zPt8hYc3GxDZRRW1kk/jWfS4WuJfSBJS6tdr/sn8zMs1m8Hxf3wCwJ25ufKUU4tnO1ZW81zXzzxeOSJB83m7w/OSaH8XkmTti1Dq+uh75o8/Dja9vnnPz0VO5eYzl+pc1YroZEWVR6h9AX260/ubHM/c7Jn/JWhO9ULvzJHfjHbdriDgejgjg7yf8R9VZoLgx5okn/Y2VxPYEtuF/l390o8L80LvRITih/yS71jL+tG6GWMRkZDQ2PJ/+bH5x+LXXB7oJsgGCuCRV+sHxGDhfoefAv/iF+bThS+3uMtI0LWsywT5j52Q/IG/8wXvLRuhLmrRNi/XyzGLNEQVs4iI+Qn/WTKDvSDQ101L+S7WwTlO1vbOvwzQbqal/JZlp3cXAdrdtJTvslwGewGg7bRUUDN2k9vv3K4dA73MtJTvsoRW1SnQ23UZZ/rieHrXGdBzhFxhF7eyBUdAvwyXcaZfzrSqToCeov3cWnGlVb0a6JfnMs700QHpuxLoudrPrZXrg5mrgN5WkH1tuU6rOhvoLQbZ15eW9E1+gmcC/ZIY88SezyR9s4B2rf3cWpkTzMwAWmhiHrh7U2EGEJKMGLi5FOr/A71QhzntW8JmAAAAJXRFWHRkYXRlOmNyZWF0ZQAyMDIzLTA1LTI3VDE5OjI2OjI1KzAwOjAwfV4keAAAACV0RVh0ZGF0ZTptb2RpZnkAMjAyMy0wNS0yN1QxOToyNjoyNSswMDowMAwDnMQAAAAodEVYdGRhdGU6dGltZXN0YW1wADIwMjMtMDUtMjdUMTk6MjY6MjUrMDA6MDBbFr0bAAAAE3RFWHRwZGY6VmVyc2lvbgBQREYtMS41UzZawQAAAEp0RVh0c2lnbmF0dXJlADMzOGViODQ4Y2E0N2VhNTIyYzUwYzU4NmVjMGRjY2EzZTE1ZWRmYzJhZTExNDYyZjU0ZGQ2ZTBmNDg1MWQyZTR0Pq5GAAAAAElFTkSuQmCC",
      "text/plain": [
       "<wand.image.Image: 338eb84 'PDF' (180x111)>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "feynman:\n"
     ]
    },
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAgcAAAIHCAYAAAALof87AAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjcuMSwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/bCgiHAAAACXBIWXMAAA9hAAAPYQGoP6dpAABiNUlEQVR4nO3dd1hU17oG8HfoiIoVbIgtKkqzV+wJVuwCKrZYo0nsJjFGTXLUaGyJGntXsItdsGAFCwKCFRQLoCAWQOkw94+5ji5mgAGBGeD9Pc99bvbae80s7s05vOy19/dJpFKpFERERET/T0vdCyAiIiLNwnBAREREAoYDIiIiEjAcEBERkYDhgIiIiAQMB0RERCRgOCAiIiIBwwEREREJGA6IiIhIwHBAREREAoYDIiIiEjAcEBERkYDhgIiIiAQMB0RERCRgOCAiIiIBwwEREREJGA6IiIhIwHBAREREAoYDIiIiEjAcEBERkYDhgIiIiAQMB0RERCRgOCAiIiIBwwEREREJGA6IiIhIwHBAREREAoYDIiIiEjAcEBERkYDhgIiIiAQMB0RERCRgOCAiIiIBwwEREREJGA6IiIhIwHBAREREAoYDIiIiEjAcEBERkYDhgIiIiAQMB0RERCRgOCAiIiIBwwEREREJGA6IiIhIwHBAREREAoYDIiIiEjAcEBERkYDhgIiIiAQMB0RERCRgOCAiIiIBwwEREREJGA6IiIhIwHBAREREAoYDIiIiEjAcEBERkYDhgIiIiAQMB0RERCRgOCAiIiIBwwEREREJGA6IiIiyEB8fr+4lFDiGAyIiokxcv34dNWrUgKurq7qXUqAkUqlUqu5FEBERaZr4+Hg0atQIDx8+BAA4OjpizZo1KFeunJpXlv9454CIiEiJWbNmyYMBAOzZswcbN25U44oKDu8cEBERZeDh4QF7e3thrFmzZrhy5Qp0dXXVtKqCw3BARET0mTdv3sDKygoRERHyMUNDQ/j5+aFevXpqXFnB4bYCERHRZyZOnCgEAwBYvHhxsQkGAO8cEBERybm5ucHZ2VkY++abb3Dy5EloaRWfv6cZDoiIiACEh4fD0tIS7969k4+VLVsWgYGBqFq1qvoWpgbFJwYRERFlQiqVYtSoUUIwAIA1a9YUu2AAMBwQERFhzZo18PDwEMacnJzg5OSkphWpF7cViIioWHvw4AEaNWqEhIQE+ViVKlUQGBhYLAoeKcM7B0REVGylpKTAxcVFCAYAsGXLlmIbDACGAyIiKsYWLFiAGzduCGMTJ07EN998o6YVaQZuKxARUbF048YNtGrVCmlpafKxunXrws/PDyVKlFDjytSPdw6IiKjYiY+Ph4uLixAMtLW1sXPnzmIfDACGAyIiKoZ++uknPHjwQBj79ddf0axZMzWtSLNwW4GIiIoVT09PhWcKmjZtiqtXrxaLpkqqYDggIqJi4+3bt7CyskJ4eLh8zMDAAH5+fqhfv74aV6ZZuK1ARETFxqRJk4RgAMiaKjEYiHjngIiIioU9e/YoVDzs0qULTp8+XayaKqmC4YCIiIq88PBwWFlZ4e3bt/KxMmXKIDAwENWqVVPjyjQToxIRERVpUqkU3377rRAMAFk/BQYD5RgOiIioSPvvv/9w+vRpYczR0RHOzs5qWpHm47YCEREVWQ8fPoStrS2bKuUQ7xwQEVGRlJqaqrSp0ubNmxkMssFwQERERdLChQtx/fp1Yey7776Dvb29mlZUeHBbgYiIihxfX1+0bNkSqamp8rGvvvoKfn5+MDIyUuPKCgfeOSAioiIlISEBQ4cOFYKBtrY2duzYwWCgIoYDIiIqUn766Sfcv39fGPvll1/QokULNa2o8OG2AhERFRlnzpzB119/LYw1adIE3t7ebKqUAwwHRERUJLx79w5WVlYICwuTj7GpUu5wW4GIiIqESZMmCcEAAP766y8Gg1zgnQMiIir09u3bh0GDBgljnTt3hoeHB5sq5QLDARERFWovXryApaUl3rx5Ix8zNjZGYGAgzMzM1LiywotxioiICq2PTZU+DwYAsHr1agaDL8BwQEREhda6detw8uRJYWzgwIEYPHiwmlZUNHBbgYiICqXg4GDY2toiPj5ePla5cmUEBgaifPnyalxZ4cc7B0REVOh8bKr0eTAAZE2VGAy+HMMBEREVOosWLcK1a9eEsQkTJqBr165qWlHRwm0FIiIqVNhUKf/xzgERERUaCQkJcHFxEYKBlpYWtm/fzmCQhxgOiIio0Pjll19w7949hbGWLVuqaUVFE7cViIioUDh37hw6d+4sjDVu3Bg+Pj5sqpTHGA6IiEjjvXv3DtbW1nj+/Ll8TF9fH7du3UKDBg3UuLKiidsKRESk8b7//nshGACyNxYYDPIH7xwQEZFG279/PwYOHCiMderUCZ6enmyqlE8YDoiISGOxqZJ6MHIREZFGyqyp0qpVqxgM8hnDARERaaT169crNFUaMGAAhgwZoqYVFR/cViAiIo0TEhICGxsboXdCpUqVEBQUxN4JBYB3DoiISKOkpqZi2LBhbKqkRgwHRESkURYvXgxvb29hbNy4cejWrZuaVlT8cFuBiIg0xq1bt9CiRQuhd0KdOnXg5+eHkiVLqnFlxQvvHBARkUZITExU2lRpx44dDAYFjOGAiIg0wi+//IK7d+8KYz///DObKqkBtxWIiEjtzp8/j06dOgljjRo1go+PD/T09NS0quKL4YCIiNQqJiYGVlZWbKqkQbitQEREavXDDz8oNFVauHAhg4Ea8c4BERGpzcGDB9G/f39hrGPHjjhz5gybKqkRwwEREanFy5cvYWlpidevX8vHSpcujcDAQFSvXl2NKyPGMiIiKnAfmyp9HgwAWVMlBgP1YzggIqICt2HDBpw4cUIY69+/P4YOHaqmFdHnuK1AREQF6tGjR7CxscGHDx/kY6ampggKCkKFChXUuDL6iHcOiIiowKSlpWHYsGFCMACATZs2MRhoEIYDIiIqMIsXL8bVq1eFsbFjx6JHjx5qWhEpw20FIiIqEP7+/mjevDlSUlLkY7Vr14a/vz97J2gY3jkgIqJ8l5iYiKFDhwrBQEtLC9u3b2cw0EAMB0RElO9+/fVX3LlzRxibNWsWWrduraYVUVa4rUBERPnKy8sLHTt2FMZsbW1x7do1NlXSUAwHRESUb16/fo1atWohNjZWPqanp4dbt26hYcOGalwZZYXbCkRElC8iIiLQoEEDIRgAsqZKDAaajXcOiIgoz3l4eKB///54//69MG5hYYGgoCA2VdJw/P8OERHlmdTUVPz666/o2rWrQjAAAFdXVwaDQoB3DoiIKE9ERERg0KBBuHr1KpT9aqlduzZCQkLUsDLKKcY3IiL6Yh4eHrC0tMS1a9eUBgMAGDVqVAGvinKL4YCIiHItNTUVc+bMQdeuXRETE4PU1NRMr+3Zs2cBroy+BLcViIgoVyIiIuDo6IgrV65kerfgo0qVKiEiIgISiaSAVkdfgncOiIgoxz5uI/j4+GQbDHR0dODg4MBgUIgwHBARkcpyso3w+Zzu3bsXwOoor3BbgYiIVJKTbYTP6ejo4O3bt2ywVIjwzgEREWXL09MTlpaW8Pb2zlEw0NLSQocOHRgMChmGAyIiytKKFStgb2+Pt2/fIi0tLcfz+ZZC4cNwQEREWYqKilK4W6CtrQ1dXd1s56anp6NHjx75tTTKJ3zmgIiIspWUlISIiAg8f/4cYWFhCAsLw9atW3Hv3r0s59WoUQOhoaEFtErKKzrqXgAREWk+fX191KxZEzVr1gQAHD58ONtgoKOjgz59+hTA6iivcVuBiIhyJDIyEmPGjBHG9PX1FeoYpKamckuhkGI4ICIilUmlUowZMwbR0dHC+MqVK1G+fHlhzMDAAHZ2dgW5PMoj3FYgIiKVbd68GUePHhXG+vTpg7FjxyIpKQmTJ0+GVCqFtrY27O3toa+vr6aV0pfgA4lERKSSx48fw8bGBu/fv5ePmZiYICgoCBUrVkRiYiLMzc0RFRUFANiwYQNGjx6truXSF+C2AhERZSstLQ3Dhw8XggEAbNy4ERUrVgQg20aYO3eu/Fy3bt0KdI2Ud3jngIiIsrV48WLMmjVLGPv222+xceNGYSwpKQk1a9ZE+fLlERgYWJBLpDzEcEBERFkKCAhAs2bNkJKSIh+rWbMmAgICUKpUKYXrExISYGBgwC6MhRgfSCQiokwlJSXBxcVFCAYSiQTbt29XGgwAwNDQsKCWR/mEzxwQEVGm5syZo7A9MHPmTLRt21ZNK6KCwG0FIiJS6uLFi+jQoYPQV8Ha2hrXr1/nK4pFHMMBEREpiI2NhY2NDZ48eSIf09PTw82bN2FlZaW+hVGB4LYCEREpmDJlihAMAODPP/9kMCgmeOeAiIgE7u7uCg2T7OzscP78eWhra6tnUVSgGA6IiEguKioKlpaWePXqlXysZMmSuH37trwjIxV93FYgIiIAn5oqfR4MAOCff/5hMChmGA6IiAgAsGXLFhw5ckQY6927N0aMGKGeBZHacFuBiIgQGhoKa2trhaZKgYGBMDExUePKSB1454CIqJjLrKnShg0bGAyKKYYDIqJibtmyZbh06ZIwNmrUKDg4OKhpRaRu3FYgIirGbt++jWbNmiE5OVk+VqNGDQQEBKB06dJqXBmpE+8cEBEVUx+bKn0eDD42VWIwKN4YDoiIiqm5c+fi9u3bwtj06dNhZ2enphWRpuC2AhFRMXTp0iW0b99eaKpkZWWFGzdusKkSMRwQERU3cXFxsLGxQWhoqHxMT08PN27cgLW1tRpXRpqC2wpERMXMlClThGAAAH/88QeDAcnxzgERUTFy5MgR9O7dWxhjUyXKiOGAiKiYePXqFSwtLREVFSUfY1MlUobbCkRExYBUKsXYsWOFYAAAK1asYDAgBQwHRETFwLZt23D48GFhrFevXhg1apR6FkQajdsKRERF3JMnT2BtbY24uDj5WMWKFREYGAhTU1M1row0Fe8cEBEVYR+bKn0eDABg/fr1DAaUKYYDIqIibPny5bh48aIwNnLkSPTp00c9C6JCgdsKRERFVGBgIJo2bcqmSpRjvHNARFQEZdZUadu2bQwGlC2GAyKiImjevHkICAgQxqZNm4Z27dqpaUVUmHBbgYioiLly5QratWuH9PR0+ZilpSVu3LgBAwMDNa6MCguGAyKiIiQuLg62trZ4/PixfExXVxc3btyAjY2NGldGhQm3FYiIipBp06YJwQAAfv/9dwYDyhHeOSAiKiKOHj0KBwcHYaxNmza4cOECmypRjjAcEBEVAZk1VQoICECtWrXUuDIqjLitQERUyEmlUowbN06hqdLy5csZDChXeOeAiKiQ27ZtG0aMGCGM9erVC+7u7pBIJOpZFBVqDAdERIXY06dPYWVlJfROqFChAoKCgtg7gXKN2wpERIVUeno6RowYwaZKlOcYDoiICqkVK1bAy8tLGBs+fDj69u2rngVRkcFtBSKiQigoKAhNmzZFUlKSfKx69eq4ffs2jI2N1bgyKgp454CIqJBJTk7G0KFDhWAgkUiwfft2BgPKEwwHRESFjLKmSlOnTkX79u3VtCIqaritQERUiFy9ehV2dnZsqkT5iuGAiKiQeP/+PWxsbBSaKl2/fh22trbqWxgVOdxWICIqJJQ1VZo/fz6DAeU53jkgIioEjh8/jp49ewpjrVu3xsWLF9lUifIcwwERkYaLjo6GpaUlIiMj5WNGRkYICAhA7dq11bgyKqq4rUBEpME+NlX6PBgAsqZKDAaUX3jngIhIg23fvh3Dhw8Xxnr27IkjR46wqRLlG4YDIiIN9ezZM1hZWSE2NlY+Vr58eQQFBaFSpUpqXBkVddxWICLSQB+bKn0eDABZUyUGA8pvDAdERBpo5cqVOH/+vDA2bNgw9OvXT00rouKE2wpERBrmzp07aNKkCZsqkdrwzgERkQZJTk6Gi4uLEAwAYOvWrQwGVGAYDoiINMjvv/8OPz8/YWzKlCno2LGjmlZExRG3FYiINIS3tzfatm0rNFVq0KABfH192VSJChTDARGRBnj//j1sbW3x6NEj+Ziuri6uXbuGRo0aqXFlVBxxW4GISANMnz5dCAYAMG/ePAYDUgveOSAiUrMTJ06gR48ewlirVq1w8eJF6OjoqGlVVJwxHBARqVF0dDSsrKzw8uVL+ZiRkRH8/f1Rp04dNa6MijNuKxARqYlUKsWECROEYAAAS5cuZTAgteKdAyIiNdm5cydcXFyEse7du+PYsWNsqkRqxXBARKQGz58/h5WVFWJiYuRj5cuXR2BgICpXrqzGlRFxW4GIqMB9bKr0eTAAgHXr1jEYkEZgOCAiKmD//PMPzp07J4y5uLigf//+aloRkYjbCkREBeju3bto3Lix0DvBzMwMgYGB7J1AGoN3DoiICgibKlFhwXBARFRA/vjjD9y6dUsYmzx5Mjp16qSmFREpx20FIqIC4OPjgzZt2ghNlSwsLODr6wtDQ0M1roxIEcMBEVE++/DhA2xtbRESEiIf09HRwbVr19C4cWM1roxIOW4rEBHlsxkzZgjBAADmzp3LYEAai3cOiIjy0cmTJ9G9e3dhrGXLlrh06RKbKpHGYjggIsonr1+/hqWlpdA7oUSJEggICGDvBNJo3FYgIsoHRbmpUmxsLHbu3InvvvsOCQkJ6l4O5QPeOSAiyge7du3C0KFDhbFu3brh+PHjhbKpUmxsLI4cOYI9e/bg1KlTSE1NBQCcOnUK9vb2al4d5TWGAyKiPKasqVK5cuUQFBRUqHonKAsE2traSEtLg0QiQePGjXHjxo1CGXYoa3wahogoD6Wnp2PkyJEKTZXWrl1bKIJBVoEAgPx/S6VSLF26lMGgiGI4ICLKQ6tWrcLZs2eFsSFDhmDgwIFqWlH2VA0EH+no6ODrr79G+/bt1bFcKgDcViAiyiP37t1Do0aNhN4JlStXxsWLFzWud0JcXBxOnToFd3d3nDt3TiEQZEUikcDLywsWFhYFsFLNUbFiRXUvocAwHKiDVAq8eQOUL6/ulRBRHklJSUGrVq3g6+ur7qVQPilOvy75KqM6bN4M1KsHHD2q7pUQUR6Jjo5W9xKI8gzvHBS0+/eBJk2A+HjZ8cSJwJIlABuvEOWbgIAAHD16FNHR0ahatSratm2LFi1aQEsr7/4+Cg4OxpkzZ/Ddd9/l2WeSZilOvy4ZDgpSUhLQsiXg7y+OW1oCbm5Aw4ZqWRZRUSWVSjFr1iwsWbJE4VyNGjUwbdo0fPvtt3nSFXH27NlYsGDBF38Oaa7i9OuS4aAgxcYCI0cCBw8qnjMwAJYtA8aPB/hqEFGe2LhxI8aMGZPlNZUqVcL//vc/jBgx4ovuJBw6dAj9+vWTH5cqVQrBwcF5enfiS0mlUvj5+cHd3R2HDh3Cy5cvVX4I8SNtbW106NABrq6u+bhSzcQHEin/SKXA+vXAlCmAsrKjffoAGzfyYUWiLySVSlG7dm2EhoaqdH3Lli2xZs0aNGrUKFffFxYWBjMzM2EsODhYY0slS6VS3LhxA3v37oWrqysiIiKgo6Mjr3yYnQsXLqBdu3b5vEpSF4YDdbl7F3B2Bm7fVjxXtSqwYwfQsWPBr4uoiHj8+DFq164tjHXu3Bnh4eG4f/++0jk6OjqYPXs2Zs+eDV1d3Rx9n1QqRcWKFfH69Wv52MGDB9G3b9+cL76A5TQoaGtro3Hjxrh27RqLIBVRmnO/q7hp0AC4dg344QfFc+HhQOfOwOzZQEpKwa+NqAh49OiRcFymTBl4enri7t27uHDhArp06aIwJzU1FfPnz0eLFi0QGBiYo++TSCSwtrYWxnL6GeoikUjQvHlz/P333wgLC8O1a9fw448/okqVKgCg0Fo6LS0NN27cwKFDh9SxXCoADAfqZGAArFwpe6WxQgXxnFQKLFgA2NkBjx+rZ31EhVjGboFly5aFRCKBRCJBu3bt4OnpiZMnTyq97e/n54cmTZpg4cKFKt9mB6AQDm4ruzOo4XISFGbMmJGj//tQ4cFwoAl69pRtLyj5SwbXrgG2tsDu3QW+LKLCLD09PdtrunbtiqCgIPz666/Q1tYWzqWkpOCXX36BnZ0dHj58qNJ3WllZCccPHjxQfcEaKLug8PjxY3h4eKh5lZQf+MyBJklPB5YuBX75BVCWxocNA1atAkqVKvi1ERUyly9fhp2dnfy4ZMmSiIuLy/T6mzdvYvjw4bh7967COUNDQ/z111+YOHFilm8fnDt3Dp07d5YfGxsb4927d7n7ATTYx2cUrl69ih9++EGj3sigvMFwoIlu3pQ9rBgSoniuTh3ZXYRmzQp+XUSFSEhICL766ithLCYmBqVLl850TmJiIubNm4clS5YovfPQqVMnbNmyBdWrV1c6Pzg4GHXr1s3RdxJpIsY9TdS0KXDrFjB8uOK5kBCgdWtg8WLZnQYiUqpatWoKWwXZPSBoYGCARYsW4fLlywrBApDdGbCyssKmTZuUhofySl5Bfv/+fQ5XTqR+DAeaqlQpYOtWYNcuxW2E1FRg1izA3h548UItyyPSdAYGBgpdA/38/FSa26pVK/j5+WHSpEkK52JjYzF69Gi0adMGN2/eFM5lfKofUGx3TFQYMBxousGDZeWWW7RQPHfmDGBtDRw7VuDLIioMMhY0UjUcAICRkRH+/fdfeHp6KhQ3AgAfHx80b94cw4cPR1BQEADg7du3Ctfp6+vncNVE6sdwUBjUqgVcuiSre5Cx4Eh0NNCrl6xeQmKietZHpKG+JBx81KVLFwQGBmK4km0+qVSK7du3w8rKCi1btoSDg4NwvlSpUsWq5C4VHQwHhYWuLvDnn8DZs7IKihn9+6/s7oKSJ62JiquM4SAoKAjJyck5/hxjY2Ns3boVJ06cUPosAgBcu3ZNoa6Bvb09KwhSocRwUNh07AgEBAC9eyueu31b9jDjunWyIkpExZytra1wnJKSovRVRVV169YNgYGBWLRoEYyMjLK8VktLCzNnzsz1dxGpE8NBYVS+PHDoELBmjazK4ucSEmSdHQcMAN68Uc/6iDREmTJlULNmTWEsN1sLn9PX18esWbPw+PFjzJ07FxUyVjf9/2s2btyIZnzlmAop1jko7IKCZDUR/v+BKEG1asDOnUD79gW/LiIVvHnzBmfOnEFUVBSqVKmC1q1bo1KlSnn6HQMGDMCBAwfkx99//z3++eefPPv8pKQkXLx4Ed7e3njz5g3Mzc3Rv3//TGshEBUGDAdFQUICMHOmrHpiRhKJ7EHGuXMBJa9ZEanL1q1bMWnSJHz48EE+JpFI0L59e8yYMQPdunXLk/36//3vf/j111/lx23btsWlS5e++HOJijKGg6LkyBFg1Cjgs5axcq1aySor1qhR4MsiysjDwwP29vZZXmNra4tly5ah4xe2Lj9x4gR69OghPy5ZsiRiYmJY8pcoC/xPR1Hi4CB7WLFTJ8Vz3t6AjQ3g5lbw6yLKYO7cudle4+/vj06dOsHR0RHh4eG5/q6Mbyy8f/8eIcpKkxORHMNBUVO1KuDhASxcqLiNEBsrez5h5EiAJV1JTV6/fg0fHx9hzNzcHKUyaSi2d+9eNGzYENu2bUNubnRWrlwZJiYmwlh2ZZSJijuGg6JIWxv46SfgyhVZAaWMtm4FGjcGfH0LfGlE9+7dE4719fURHByMqKgobNmyRaFxESBrXjRixAg4ODggIiIix9+ZsZUywwFR1hgOirLmzQE/P2DoUMVzwcGy5xD+/psNnKhAxcTECMcmJibQ1dWFgYEBRowYgbt37+K///5DuXLlFOYeO3YMlpaW2LVrV47uIjAcEOUMw0FRV7o0sGOH7H9KlhTPpaQAM2YA3boBL1+qZ31U7GRsTpSSkiIca2trY/z48Xj48CGGDRumMP/t27cYOnQoBgwYgKioKJW+M2M4uHPnTg5XTVS8MBwUF0OHyho4NW+ueM7DQ9bA6cSJAl8WFT9ly5YVjpU1KwJk7Y+3bduGw4cPw9TUVOH8wYMHYWlpiYMHD2b7nRlLHj9//jxXzy8QFRcMB8VJ7drA5cvAzz8rNnB69Qro0QOYPBlISlLL8qh4yFhRMCkpCZGRkZle37t3b9y5cwfOzs4K5169eoX+/ftj6NChmYYMAKhWrZpwHB8fj3fv3uVs4UTFCMNBcaOrCyxYIGv3XLmy4vmVK2UNnDI8NEaUV8zNzVGiRAlhzN/fP8s55cuXx+7du7F//36l5Yp37doFS0tLnDp1Sul8ZXceXiurB0JEABgOiq9OnWSNmnr1UjwXEAA0aQJs2MAGTpTntLW1YW1tLYxlFw4+6t+/P+7cuYM+ffoonIuIiEC3bt0wbNgwhTcatLW1c7tcomKJ4aA4q1ABcHeXlV3W1xfPJSQAY8cCgwYBWdyuJcqNjIWJVA0HgOzthoMHD2LHjh0wNjZWOL9jxw7UrVsXc+fORXR0NAAIJZo/YmAgyhzLJ5NMYKCsQJKyp7jNzIBduwA7u4JfFxVJ69evx7hx4+TH9evXV6h/oIqwsDCMHj0ap0+fVnre0NAQffv2RVJSktB8SUdHB/Hx8dDV1c354omKAYYD+iQhAZg2DfjvP8VzWlrAr78Cc+awgRN9sevXr6NFixbyY4lEgri4OBgZGeX4s6RSKTZs2IAZM2YgNjZWpTnt27eHl5dXjr+LqLjgtgJ9YmgIrFkDHDoEZCxAk54O/P470KED8PSpWpZHRYelpaXQ+Egqlea6MJFEIsHYsWMRHByM0aNHq9TJ8ccff8zVdxEVFwwHpKhPH9lDiR06KJ67ckXWwGnv3oJeFRUhJUqUQL169YSxnDx3oIyJiQk2bNiAgIAAuLi4KBRb+mjmzJno27fvF30XUVHHbQXKXFoa8NdfwG+/yf45o1GjgH/+AXJxK5g0W3p6OgIDA/HixQtUqlQJDRs2zPP9+SFDhmD37t3y43HjxmHt2rV59vkvXryAu7s7vL29ER0djerVq8PZ2Rnt2rXLs+8gKqoYDih7Pj7A4MFAaKjiubp1AVdXWSMnKhIuXbqEUaNGCW2NjY2N0a9fP0ydOhWWlpZ58j1LlizBzJkz5cctWrRQ6NZIROrBbQXKXsuWsgZOgwcrnnv4UHZ+2TLVGjhJpbKeDqSRgoKCYG9vLwQDQNYsacuWLbCyskKfPn3w4MGDL/4uW1tb4fj27dtIU3aHiogKHMMBqcbYGNi5E9i2TXkDp2nTZOWXsyiDC6lUVsK5XTsgNTV/10u58uuvvyIhISHLa9zd3WFlZYWZM2ciLi4u199lY2MjHCckJODhw4e5/jwiyjsMB6Q6iQQYNkx2F6FpU8Xzp07JGjhlUsIWb9/KtiZ8fID//S9/10o5lpiYiBMqNt9KSUnBkiVLYG1tjfPnz+fq+0xMTFClShVhzM/PL1efRUR5i+GAcq5OHdlbC5/tF8tFRclaQE+dqtjA6fNnFubPl30GaYwHDx4otE9+8+YNnj9/jt9++02hmyIAPHnyBJ06dcKkSZPw/v37HH9nxq2Fu3fv5vgziCjvMRxQ7ujpyd5k8PQEKlVSPL98OdCqFfD53vTjx5/+WUsLcHQE2BlPY2TcIjA2NkbZsmVRrVo1zJ8/H0+fPsWsWbOUvrWwevVqWFtb48KFCzn6TgsLC+H4/v37OV84EeU5hgP6Ml26yBo49eiheM7PT/YWw6ZNsucNQkOBj/Xs09KAly+BcePY3ElDpGd4oPTzIkUAUKpUKSxatAhBQUFKXwcMDQ1Fx44dMXXq1GyfW/goYzjITQllIsp7DAf05SpWBI4eldU8yNjAKT4eGD0acHKStYH+/BdOWpqsmNL27QW7XlIqYxvluLg4hcAAAHXr1sX58+excuVKGBoaCuekUimWL1+ORo0a4fr169l+Z/369YXj4OBgvrFApAFY54Dy1u3bn4JARgYGQGKi8vHAQNmzDKQ2ERERqFq1qjD2sQhSZkJCQjBy5EhcvnxZ4Zy2tjZ+/vlnzJkzB3p6eip/Z1RUFCpWrJiLn4CI8grvHFDesrYGbt6UbRdkpCwYALLXGgcNApKT83dtlCVTU1OFNsaPP39ORIk6derAy8sLS5cuhX6Gu0ZpaWn4888/0bx5c9y+fVvp/PLlyyuMvX79OocrJ6K8xnBAea9ECWDtWuDAAUDJE+4KUlMBf39g7tx8XxplTltbO1f9DrS1tTF16lT4+fmhqZJXXAMCAtC0aVPMmTMHHz58EM7p6+sr3FXIzVsPRJS3GA4o/3TpAowZo1qLZ6lU9vbDuXP5vy7KVOMMZbBv3bql8lwLCwt4e3vjjz/+UGh6lJKSgj///BP169fH1q1bkfz/d4mSk5Pl//yRgYFBLldPRHmFzxxQ3nv9Gli5ElixAvjwQbWyyoDsYcUKFYC7dwElt5sp/y1btgzTpk2THzdu3Bi+vr45/hw/Pz8MGzYMQUFBSs9Xq1YNI0aMQEpKCv766y/h3Lt372BsbJzj7ySivMNwQHnn7Vtg4ULg339lzw+oGgo+p60tey3y8GFZRUYqUF5eXujYsaP8WE9PD3FxcZk+UJiVpKQkzJs3D0uWLFH5DQQLCwsWQiLSANxWoLzj7w8sWSJ78DA3wQCQvd545Aiwfn2eLo1Uk7FiYXJycq5/Wevr62PhwoW4ffs2OnfurNKcCRMm5Oq7iChvMRxQ3mnZEti1C+jfX/Z6IqDa8wbK/PCDbHuBClSZMmVQs2ZNYexL+x00aNAAnp6eOH78OOzs7DK9rlevXgwHRBqC2wqUPxITgbNngUOHgIMHZVsOOjo568Zobg7cv/8paBAA2dP8kZGRqFixIkqXLp3nnz9gwAAcOHBAfjxp0iT8+++/efb5169fx759++Dt7Y2oqCiYmZlh0KBB+PbbbxUeZCQi9WA4oPyXlgZ4e8uCwr59wPPnsmcLVNmHbtxYNjcXe95FTWhoKL7//nucOHECH/9ja2VlBUdHR4wfP15pzYDc+N///odff/1VftymTRulRY6IqOhiOKCCJZUCQUGyoLB/v6wyopaWbDyzfxVr1wZOnADq1i3YtWqQly9fonHjxnjx4oXS84aGhhg7dizmzJnzxSHh5MmT6N69u/zYyMgIMTExCgWSiKjoYjgg9Xr6VPZmwoEDwOXLmQeEEiWA1auB4cOL5VsMEyZMwNq1a7O9rmzZsvj9998xYcKEXP8yf/nyJSpXriyMPXjwAHWLcTgjKm4YDkhzREfLGjgtXQrcuaP8GicnWfXFYvQefFpaGipUqIB3OWhv3bx5c2zZsgUNGjTI1XdWqlQJkZGR8uO9e/di4MCBufosIip8+LYCaY4KFYCRI2XbDlevAhka8gAA3NwAW1vZcwjFxLNnzxSCwfXr13Hx4kUMHTpU6UN8169fR6NGjbBgwQKk5uQh0P9nY2MjHKtSRpmIig6GA9JMrVoBDx7Iyi9n9OQJYGcH/PGHag81FnLR0dHCsZ6eHpo2bQo7Ozvs2LEDISEhcHR0VJiXnJyM2bNno2XLljmuVZAxHAQEBOR84URUaDEckOYyMpIVQ9q3DyhTRjyXlgb89hvQubPs7YdiRFtbG5LPnrswNzeHm5sbvLy8lD4X4Ovri8aNG+eoUmHGcMCqhUTFC8MBab4BA4CAANndgowuXABsbGRvPxRRGUsXJyUlKd0qaN++Pfz9/TFjxgxoaWkpzJk5cybatWuH4ODgbL/zq6++Eo6fP3+ucrAgosKP4YAKh+rVZR0b58+Xvfr4ubdvgX79gPHjgfh44VRcUhxC34YiLDYMaemF85dblSpVhOP09HS8fPlS6bWGhoZYvHgxrl69ivr16yucv3r1KmxsbLBq1SqkZ1HiukaNGsJxamoqIiIicr54IiqU+LYCFT5XrgCDBwPPnimes7BA4s6tWJfsjc3+mxEYGQgpZP+K62vro7VZa/S36A8XGxeU1s/76oL5QSqVwsDAQGhtfPnyZbRp0ybLeQkJCZgzZw6WLVsGZf8x79ixIzZv3qwQBD5+p46OjhAg/P39FbYbiKho4p0DKnzatJFtMwwapHAqJPIemm5uicmnJ+N25G15MACApLQknH9yHpNOTkK1ZdWw4NICJKQkFOTKc0UikaB27drC2O3bt7OdZ2hoiL///hsXL15UmA8A58+fh5WVFdavX69wF0EikcAgQ9nqpKSkXKyeiAojhgMqnMqUkb3WuGmTrEASgEgjoPMw4E7F7G+GxSXHYfa52Wiyvgn8XnxZY6GC0LhxY+H41q1bKs9t27YtAgICMHHiRIVz79+/x7hx49CsWTN4eXnJx6VSqcIzBrzJSFR8MBxQ4SWRAKNGAbduAY0aYUpX4FmZnH3Eveh7aLWpFfYE7cmXJeaVLwkHgKwE8qpVq3DmzBlUr15d4fytW7fQsWNHtGjRAps2bcLOnTsV7hRkfPaBiIouPnNARcLjyHuovVasBmgZCaw8BTQLB+J//A4+wzpi253dOHz/sLDd8NGa7mswoZlmtgz28vJCx44d5cd6enqIi4tTeJNBFbGxsZg6dSo2bdqk8pzKlSsjPDxceIWSiIou3jmgImF/yFHh2DgROLMd6BQKlEoGTJesQe+Rf+Fgk8W4OfYmmlRuovAZ3534Dtv8txXUknPE1tZWOE5OTs517YHSpUtj48aNOHXqlNI3GpRxcXFhMCAqRhgOqEi4Hn5dOB4UUw2mHzJcdPMm0KgRGnsG4eqoK5jQVPEuwZijY3D5mea1Jy5Tpgxq1aoljOV0ayEje3t73L59G6tXr4aZmVmm11laWgotnImo6GM4oCIhPC5cOG404idg2TJAV1e88P17YPhw6A0fhdVtF2Je+3nC6ZT0FAzYOwCvPrzK1Tryc5fuS587UEZXVxffffcdHj9+jAMHDqBfv34wMTGBlpYWqlatih9//BGXL19GqVKlvvi7iKjwYDigIiEpVXx4zlDPCJgyBbh2DahXT3HC7t2QNGqE3/S/wdSWU4VTkR8iMe7YOJV/0cfFxWHGjBkwMzODtrY2TExM0LdvX7i6uiIlJSXXP1NGGcOBn1/evWWho6ODfv364cCBA4iMjERaWhrCwsKwYsUKGBejDphEJMNwQEVCWcOywrH8L/9GjQBfX+DbbxUnhYZCYmeHxTfLwr7WN8KpQ/cP4djDY9l+b0JCAjp06IC///4bYWFhkEqlePXqFQ4fPozBgwejZs2aWLFihVDAKLcaNWokHPv7+7OkMRHlC4YDKhJqlqkpHPtH+n86MDICNm4E9uwBMv4VnJYG7V/nYNuO96hoUF44NeX0FKSkZf2X/+LFi7O8vR8eHo4pU6bA0tISR48ezfQ6VWQMB/Hx8Xj48OEXfSYRkTIMB1QkZHz74MzjM4q9FAYNklVWVFJ22NTzKpYeSRTGHr19hJ23d2b6nVKpFNu2qfZ2Q3BwMBwcHDBgwABERkaqNEdhjaamCrUGVKmUSESUUwwHpHHuRN3BlFNT0HxDc1RfXh21VtZCx20dMefcHNx6ofyv9G9qi9sCUR+icDb0rOKF5uaAlxcwd65CA6ehVz+gRZh4+YLLC5AuVd6g6M2bNwgNDRXG1q5di61bt2ba9+DAgQNo0KABdu3alauHF62trYXjgICAHH8GEVF2GA5IYySlJuH7E9/D8j9LrLi2AjcibuB57HOEvguF1xMv/HnpTzRZ3wQtN7bE2cfiL/7a5WrD2lT8xfnPtX+Uf5GODjBvniwkfPYKnwTAPC/x0pA3ITgXek7pxyjrUvjtt99i+PDhuHz5Mi5fvowWLVooXPPmzRsMHToUffr0ybS7YmYyNj7inQMiyg8MB6QR0tLT4LjfEaturMr22mvh19BlRxc47nfEm4Q38vFxTcYJ1x0PPg6fMJ/MP8jOTrbNMGCAfMg+RFZZ8XObbqleSVBHR0f+z23atMHVq1exbds2lC9fXuHaI0eOwNLSEvv27VP58zPeOQgMDFR5LhGRqhgOSCP8c+0fuD9wz9GcvXf2wmatDfxf+gMAhtkMQ4USFYRrJp2YhNT01Mw/pGxZYO9eYMMGwNAQEgBjfcVLjt/ej6TIcIWpuhlrKECxc6GWlhaGDRuGu3fvYuDAgQrXv379GoMGDYKTkxNev36d+Tr/X8aKhmFhYXn6uiQREcBwQBogPiUe8y/MF8YMdQzxU5ufcMTpCPYN3IcZrWegWulqCnPDYsPQdnNbnH18FiX1SuKnNj8J531f+GLxlcVZL0AiAUaPljVwsrXFwLuA5LPHAeK0U3GhlzVwTtxeqFSpksJHPX/+XOlXmJiYYO/evdi/fz9MTEwUzu/Zs0elNxoyNk1KT09HeLhicCEi+hIMB6R2xx4eQ0xSjDg2+BgWdlmIXvV6YUCDAVj89WKE/hiK9T3Xo5xhOeHaDykf0Mu1Fy4+vYhJzSfBooKFcP6387/B64lX9gupXx/w8UGl0ZMVHky8VPIN0KUL8PPPwP//pV6mTBmFAkFPnjzJ8iv69++PoKAg9O/fX+Hcy5cv4eDggFGjRiEmJkbJbKB8+fIwMDAQxnL79gMRUWYYDkjtrjy7Ihx3qdUFnWp2UrhOR0sHY5qMQeCEQNhVtxPOJaQmoI9bHzyPfY4NvTZAS/LpX+00aRr67umLe6/uZb8YfX1g+XK0ayb+8vY2AyCVAosWAW3bAo8eAQBq1hTrKzx+/Djbr6hYsSL27duH3bt3o2zZsgrnt2zZAisrK3h4eCick0gkMDQ0FMbyosASEdHnGA5I7ULehgjHbc3aZnl9lVJV4OniiYENxD38t4lv0dutNxpXbqzQM+Fd4jt03t4Z96Pvq7Sm5nbOwvG9zx9luH5dVnlx507Url1buC4oKEilz5dIJHB2dkZQUBB69OihcP758+ewt7dH3759ERISouQTPklNzeKZCiKiXGA4ILWLT4kXjjOWQlZGX0cfu/rtQq+6vYTxu6/uYtaZWfjF7hc41HMQzr14/wLttrSD93PvbD+/XgWxH0NEaeDD588fxsUBLi6YGxKCz1sS5bQZUpUqVXD06FFs2rRJaXOjw4cPo379+nBycsL58+dx584dvH37VrjG1NQ0R99JRJQdhgNSu5J6JYXjmETl++0Z6WrrYnf/3Qr1Df69/i9uRNyAa39XtKgq1hl4Ff8KHbd1xHrf9VkWIapuXF1hLLphDYUxq4AA+AFo/v/Huel3IJFIMGrUKAQGBqJTJ8XtlLS0NOzZswedOnWCpaWlcM7IyEjh7gUR0ZdiOCC1q1qqqnB8L1qFZwP+X0m9knDr7wYDHfEhve9Pfg8DHQOcGnpKISAkpSVh3LFx6LunL57HKH+7oKReSeG5BQB4d2A3MHKkwrW1AVwGMAtA/IcPCA4OVnn9nzM3N4enpyfWrVuHChUqZD8BQL9+/aCvr5+r7yMiygzDAaldo0piQ6ELTy9kWrJYGYuKFvij4x/C2M2Im3C/744yBmXg4eKBrnW6Ksxzf+AOi9UWmOc1D+8S3wnntCRaCuEg3dAA2LwZcHUFSpeWj78sCRxuACwAYAVg6dKlePPmDXJDS0sLY8eORXBwMKZPn44SJUpkem3FihWxYMGCXH0PEVFWJNLcFHgnykRqeipOh5zGpWeXEBEnKy9sbmyO1mat0aFGBxjqGirMuR99HxarxdcPL428hLbVs34w8XPJacmw/s8aD14/kI81qtQIvmN9IZFIkJqeihkeM7Di2gql8431jTHcZjhGNRoFa1NrpEvTofOHjnBN0IQgNDRpKDt48gRvhg/EGtzE0lbAez1g5Spg4v8/DqCjowN7e3s4OzvDwcFB6fMEqnj79i02bdqEXbt2wd/fH4BsG+Lrr7/GqlWr8NVXX+Xqc4mIssJwQHnm+MPjmHhiIp7GPFV6vrR+aYyyHYWZbWaicqnKwrn6q+oLv9gHWw3Grn67cvT9bkFucD4gvmVwZdQVtDZrLT8+9vAYRrmPwqv4V5l+To0yNdCwYkMcDz4ujHu6eCJdmo6AlwG48PQCPB55ICX9U3XCiqdL45V3rMLnGRgYoGfPnnB2dka3bt0UXkVUVVxcHKKjo1GhQoVchw0iIlUwHFCeWH19NSadnKTStUa6RpjXYR6mtJwCbS1tAMCSK0sw88xM+TVaEi3c+e4O6leon9nHKEhLT0PdVXXx+O2nWgMjbEdgS+8twnWv41/jl7O/YMOtDZAij/71TwOwDMCHrC8rVaoU+vTpA2dnZ3Tp0kVpCWYiInVjOKAvdvHpRXTY2iHHv2g71uiIPQP2oKJRRbyOfw3zFeb4kPLpt2u3Ot1wfPBxSCQSlT/zr8t/4aezn0oolzEog6jpUdDVVvwlHBQVhN8v/I79d/d/eUhIA+wS7GAcbAzPU54KPRaUKVeuHAYMGABnZ2fY2dlBW1v7y9ZARJRHGA7oi0ilUjRZ3wR+L/2E8aZVmqJd9XZITU+F7wtfXH1+Vekv4Npla+PU0FOoU64OZp+djQWXxQfsdvTdgaHWQ1VeT1hsGMyWmwljni6e6FKrS6ZzQt+GYuOtjdgZuBPPYp6p/F0AoKeth34W/TCrzSzYVrIFAMTGxsLd3R1ubm7w8PBQqUhR5cqV5Q2YWrRokaNARESU1xgO6Iv4Rvii6Yamwtivdr/i946/C7/gHr99jAWXFmCz32aFkFCtdDVcGHEB5Q3Lw2K1BV68fyE/V1KvJG6OualQlCgrzTc0x42IG8J6/uj0RxYzZKRSKW5H3sb+u/vx56U/lV5ToUQFVDeujmZVmsGuuh161u0JYwNjpdcCQHR0NA4ePAg3Nzd4eXllWVvhoxo1asDJyQlOTk6wtrZmUCCiAsdwQF8k4238GmVqIOT7EPmzBBldfHoRQw4OQVis2NmoTrk6uDb6Gs6HnseAfQMUznl/663QjjkzszxnYfHVT50YO9boiHPDz2UxQ/T31b8xw3OG/FgCCV5Me4FyhuWUbk+oKiIiAvv27YObmxt8fHxUmmNhYSEPCnXr1s31dxMR5QTrHNAXuRt9Vzi2r22faTAAgHbm7XB99HWF2gYhb0LgtN8JfS36YqTtSIVzXXd2xZsE1WoHtDJrJRwHRanW7wAAYpNi8ffVv4WxAQ0GwLSk6RcFA0BWKvnHH3+Et7c3Hj9+jIULF8LGxibLOffu3cPcuXNRr149NGnSBEuWLMGzZznb+iAiyimGA/oiGUsdZ6x2qEzlUpVxZtgZhbLHno89scx7GVZ1XyXfv//I94UvOm/vjPDY8Gw/P2PL5lfxrxSKHGVmzrk5iPwgtkCe0nKKSnNzombNmvjpp5/g7++Pu3fv4rfffsu2ZsGtW7cwc+ZMmJubo23btli1ahXbNRNRvmA4oC+ip60nHCelZf+UPgCUMyyH44OPo1LJSsL47HOz8TzmOY4PPg6z0uKDhf4v/dF8Y3P4hGV9S75m2ZoKYy/iXii5UnTo3iH8c/0fYaz7V90V7kTkNQsLC8yfPx8PHjzArVu3MGPGDJiZmWU558qVK/j+++9RpUoVfP3119i0aZNCQyYiotxiOKAvYmokdgQMeZN1e+HPVStdDXsG7IEEnx64S05LxuTTk1G5ZGWcHXZW4U5ERFwE2m5uiz8u/IHktGSln6unrQcjXSNhLLs7B2cen1EooKSnrYfl9stV/nm+lEQiQaNGjbB48WI8efIEV65cwaRJk2BiYpLpnPT0dJw5cwajR4+GqakpHBwcsHv3brx//77A1k1ERQ/DAX2RjFsDV59fVemJ/I/ambdTuG1/KuQUzj85j6/Kf4ULIy6gdlmx62CaNA2/ef0G27W2OBVySun3ZSzTnJiaqPT7pVIp1txYg267uinc9VhhvwJ1y6vnIUAtLS20bt0a//77L8LDw3HmzBl8++23KFOmTKZzUlJScPToUQwZMgQmJiZwdHTEoUOHkJio/GcnIsoM31YgQUxiDHzCfPA89jkkkMDM2AzNqjRDWcOySq+/E3UHlv+JbYRvjrmJJlWaqPydcUlxqLeqnvAKY3vz9vAa4QVAVtGw/97+uPD0gtL5Laq2wI8tfkQ/i37Q15F1KDReZIzYpE+ljDOWUQaA6+HX8dOZn3D+yXmFzxxpOxKbHDZp3GuEycnJ8PDwgKurK9zd3fHhQzYlGQGULl0affv2hZOTEzp37syqjESULYYDAgC8fP8Sv5z9BbsCdyncrteSaMGuuh0mNJ2A/g36Q0frU0MiqVSK2v/URui7UPnYuCbjsLbn2hx9/9qbazHh+ARhzG+cn/zBxNT0VPx58U/8cfGPTDs2ljUoi+5fdUd78/YYe2yscC5gfADqla+HoKggXHx6Efvu7oN3mLfSzxnUcBB29dsl/JyaKD4+HsePH4erqytOnDihUlXGChUqYMCAAXBycoKdnR20tHjzkIgUMRwQbr24hW67uiHqQ1S211qaWOKfrv+gY82O8rG55+fi94u/y48NdAwQ8n0IqpbO/s2Fj5JSk1BzZU3h7sGEphOwpsca4brr4dcx4fgE3HpxS+XPBgBtiTbSpGnZXjerzSz8r9P/snwdUxPFxMTA3d0drq6u8PT0RFpa9j9rlSpV4OjoCCcnJzRr1kzj7pIQkfowHBRzL9+/hM1aG5WCweemtpyKRV0WQVdbF2GxYai5siZS0z+VCR5pOxKbe2/O0Wf+dv43/HHxUyXD8obl8XL6S4W/4NPS07A9YDv+uPiHcMfiS1Q3ro7/evyH7l91z5PPU6fo6GgcOHAAbm5uuHDhgkrPgNSqVUtebMnS0pJBgaiYYzgo5kYcHoFtAduEMV0tXViayJ4juBd9L9OH+b6u9TX2DdwHYwNjTDg2AWt9xa2Ec8POCXcYshP6NhS1/qkljHkN90L7Gu2VXp+SloID9w5gve96pc8NqKJSyUqY1moaxjcdj5J6JXP1GZosPDxcXpXx2rVrKs1p0KCBPChkV3uBiIomhoNi7NWHV6iyrIrwF3978/ZwG+Amrz+QkJKAIw+OYMHlBbgdeVvhM1pUbQEPFw+8T36P+qvqIy45Tn6uaqmq8B/vr3LZYwBotK4R/F/6y49V7YvwPOY5jj08htOPTsP9gXuW135V7iu0rd4W/S364+vaXyvUaiiqHj9+jD179sDNzQ23byv+/1KZJk2awMnJCY6OjtnWXiCiooPhoBjbdXsXhh761PHQQMcAT358AtOSpgrXpqWnYdX1VZh1ZpbCK3/tzdvDw8UD626uww+nfhDOdajRAaeGnJK/RZCdL+2LcCL4BHrs7iGMLbdfDisTK5iWNIVZabMsGyUVF3fv3oWbmxtcXV0REqJabYq2bdvCyckJAwYMgKmp4r8jRFR0MBwUYzM9Z2LJ1SXy4551e+Ko89Es51wLuwYHNweFZxTGNB6D/3r8B/ud9jgbelY4l5On/488OILebr3lx8b6xng7661Ke+BSqRTttrbD5WeX5WNWJlYIGB/APfRMSKVS3Lp1C25ubnBzc0NYWFi2c7S0tNC5c2c4OTmhb9++KFtW+WuuRFR48T2mYuxV/CvhuE7ZOtnOaVGtBS6MuAATI7Fq34ZbG7D/7n7s6rcLlUtWFs7tvbMXzgecM3124XNWJlbCcUxSDKLjo7OdBwBuQW5CMACAyS0nMxhkQSKRyBs6PX36FJcuXcLEiROzrcro6emJb7/9FqampujduzdcXV1ZlZGoCGE4oByrX6E+Tg45CUMdsQrh+OPjIYUUxwcfVyhfvP/ufnTe3hmR77NuFFTduLrCHYYn755ku6bQt6H47sR3wliNMjUw2GpwtnNJRktLS97QKTw8HB4eHhg1ahSMjTPfhklJScGRI0cwePBgmJqawsnJCYcPH2ZVRqJCjuGgGKtgKD4o+Dz2ucpzG1dujPW91gtj7xLf4aczP6FR5UY47HQYBjoGwvmrz6/CZq0NToWcyvRztbW0Uc6wnDAWkxSTydUyke8j0XVXV4X+CUu/WaqwBlKNjo6OvKFTZGQk3N3d4ezsjBIlSmQ6Jz4+Hnv27EHfvn1hamqKkSNH4vTp00hJSSnAlRNRXmA4KMbqV6gvHN+MuJmjvghDrYfCsaGjMLYtYBv8XvihS60uODH4hMLrgZEfItFtVzcMPjA40/bLGee8T878dvWdqDtos7kNHr5+KIw7WTqhb/2+Kv8slDl9fX15Q6eoqCi4ubmhT58+0NPL/C2P2NhYbN26FV27dkWVKlXw3Xff4eLFi0hPV17dkog0C8NBEROfEo+QNyEIfh2M+JT4LK9tXrW5cPw05inuvLqTo+9b0XUFSumVEsY+FjLqWLMjroy6gurG1RXmuQa5os6/dTD51GQ8fvtYOPf5q5WAYltoQPaK5aLLi9B0Q1M8evtIOGdRwQLreq7jswb5wMjISN7QKSoqClu2bIG9vT20tTOvKBkdHY3//vsP7du3R/Xq1TFt2jTcuHEjR0GUiAoW31YoAtKl6dh7Zy9WXV8FnzAfoUywpYkletXthTGNx6Bm2ZrCPKlUCvMV5sJ2wtSWU7HUfmmOvn/BpQWYfW62MPZg0gN5R8NXH15hpPtIHA8+nulndK7ZGf0s+sG+tj2arG8ibCV8LKaULk3Hnag72Hd3H9b7rkfkB8XnF2qVrYWLIy7mqHQzfblXr15h//79cHNzw6VLl1T6xV+7dm2hKiMRaQ6Gg0IuNikWTvudcDLkZJbXSSDBMJthWNB5AaqUqiIfn+ExA397/y0/NtY3RuiPoZl2YcxsDTVW1MDbxLfysemtpmPJN59ek5RKpVjvux4/nf1J4dmA7FiZWKGUfincibqT5fMHzao0wyHHQwwGahYWFoa9e/fCzc0NN27cUGlOw4YN4ezsDEdHR9Spk/1bM0SUvxgOCrHU9FR8s+ObHJUONtY3xtqea+Fk6QQACHkTgq/+FUvkzmozC4u6LMrRWqZ7TMdS7093HEyNTBExLQJaEnHnKupDFOacm4PN/psVtg9ySwIJvm/+PRZ1WQRDXcPsJ1CBefToEfbs2QNXV1cEBQWpNKdp06byqozVqlXL5xUSkTIMB4XYSp+VmHx6cq7mTms1DYu/XgwtiRYG7RuEfXf3yc/paunCb5wfGpo0VPnz7r26hwZrGghjV0ZdQWuz1kqvf/LuCZZ7L8f229tzfCfhc1/X+hp/dvpT4fkJ0jxBQUHyoPDo0aPsJwCws7ODk5MTBg4ciIoVK+bzConoI4aDQio1PRXmK8wRERchH6tYoiJWd1+NrnW6QkuihYDIAOy9sxeb/DYpfeJ/iNUQbOuzDcFvgmH1n5Xwl7yNqQ18Rvvk6FXAhmsa4u6ru/Lj2Xaz8WenP7OcE58Sj2MPj+How6M4HXJaoTCTMubG5uhv0R9DrIegceXGKq+PNINUKoWvr6+8KmN4uPK3Vj6nra0tVGUsU6ZM/i+UqBhjOCikLj69iPZbxW6FF0ZcQDvzdgrXRr6PxDSPadgVuEvh3CjbUdjosBGzz83GwssLhXMjbEdgs8NmlZ/6/9K+CI/ePEKdf8X9ZisTK/Su1xuVSlaCmbEZmlRuwmcKipD09HRcuXIFbm5u2LdvH169yj4c6unpoVu3bnByckKvXr1gZGSU7RwiyhmGg0JqmfcyTPOYJj+2MrHC7QlZd9rb4rcFY4+NVdjrX9xlMSY1n4RmG5opvMo4t/1czOswT6U1ZeyLYKRrhLif41QOF5NOTMLqG6vlx6X1SyNyeiQLGRUTqampOHfuHNzc3HDw4EHExGRd/AoASpQoAQcHBzg5OaFr167Q11etwRcRZY11DgqpjAWEbCvZZjtnZKOROOp8VKFuwKwzs+D7whd7B+5FCV2xAt78C/Mx32u+Sq+mZbzF/yHlA168f5HtPED2zMI633XC2JjGYxgMihEdHR1888032Lx5MyIjI3H48GE4OTllW5XxY1EmU1NTjBo1Ch4eHkhNzZuHXYmKK4aDQipj22RVf4l2rdMVbv3dIMGnv+alkMLlkAvMSpthd7/dwjkAmHdhHsYfG4+UtKzL4FYpVUWh38KjN9k/eJaUmgSXQy7CHQ09bT380OKHLGZRUaavry9v6BQVFQVXV1f07t07y6qMMTEx8qJMVapUwcSJE3Hp0iVWZSTKBYaDQqqsgViHIGML5az0teiL3zv+Low9efcEf178E73r98aaHmsU5qy/tR7tt7ZH6NvQTD9XS6IF05KmwtjntQ+USZem49sj38L3ha8wPr3VdKWVFan4MTIykjd0ioyMxObNm/HNN99kWZXx1atXWLNmDdq1awdzc3NMnz4dN2/mrDw4UXHGcFBI1SpbSzj2f+mfo/k/t/1Z4eHF5T7LEfImBOObjsd/Pf5TmOMd5g2btTZY6bMy0xoFxvpiB7+YxMz3jRNTE+FyyEXhQcl65evhF7tfVP1RqBgpU6aMvKFTREQEVq9eDTs7uyznhIWFYenSpWjWrBnq1q2LOXPm4M6dnJUJJypuGA4KqSZVmgjHT2Oe4kH0A5Xna2tpY0OvDdDV0pWPpaSnYMGlBQCA8U3HY9/AfQrbBHHJcZh8ejKs/rOCa6CrQkiQQvzLTFtL+V93vhG+aLWpFXYH7hbGDXUMsXfgXhjp8Ql0ypqJiYm8odOzZ8/w999/o2nTplnOCQkJwZ9//glLS0tYW1tjwYIFKtdcyMrZs2cZOKhI4dsKGiAsNgwbfDfg3JNzCH0bitT0VFQtXRVNKzdFz7o90e2rbtDR0hHmpEvTUW1ZNeGBv9/a/Yb5Hefn6LunnZ6GZT7L5Mc6Wjp4/MNjmBmbAZD9Enfc76jQ3OijaqWrYaTtSPSz6AcbUxvUXVUXIW9C5Of3D9yP/g36AwDS0tNw+dllrLqxCgfuHlAIErpaujjifARd63TN0c9A9Lng4GB5saW7d+9mPwFA8+bN4eTkhEGDBqFq1Zy9Kvvq1SuYmJjAwMAAmzZtwuDBg3OzbCKNwnCgRlKpFP9e/xczPWcqPGD4ObPSZvi57c8Y02SMEBImHp+INTc/PR9QqWQlPP7hcY5KCEe+j0TNlTWRkJogH8v4+mJcUhymnJ6CTX6bsvwsEyMThWcfHBs6onbZ2rjz6g68w7wzfTaipF5J7BmwB92/6q7y2omyExgYKC+29Pjx42yvl0gkaNeuHZycnNC/f3+VqjJKpVIYGxsjLi4OADB+/HisWLGCr1VSocZwoEZ/XPgDv3n9pvL11qbW2NVvFyxNZB3sAl4GwHadrXDN31//jWmtpymZnbkJxyZgre9a+bFZaTM8nfxUoT7B5WeX8cPJH+D30i9Hn5+dhhUbYs+APTkq10yUE1KpFDdv3oSrqyv27NmDiIiIbOdoa2ujS5cucHZ2Rp8+fWBsbJzpta1atYKPjw8AQEtLCzY2Njh48CBq1KiRVz8CUYFiOFCTy88uw25L1g9SKaOvrY/NvTdjsJXs1uXXO77Gmcdn5OdL65fG3e/u5qiKoP9LfzRa10gYuzHmBppWUdy/lUqlcH/gjoWXF+J6+PUcr/9zpfVLY1abWZjeerpC7QWi/JKeno7Lly/D1dUV+/btw+vXr7Odo6enh+7du8PZ2Rk9e/ZUqL0wbtw4bN68WV5fQUdHByVKlMDu3bvRo0ePfPk5iPITw4GatN/aHhefXhTGvm30LRzqOUBXSxeBUYE4cO9Apr+AF3dZjBltZuB6+HW02NhCONe1TlccH3xcoSNiZqRSKeqvro+Hrx/Kx1R5fsHvhR+2B2yH+wN3hL7L/BXHjKxNrTHEagjGNRkHY4PM/xojym8pKSk4d+4cXF1dcejQIcTGxmY7x8jICA4ODnB2dsY333wDfX19rFy5ElOnThVqKmhpaSE9PR2//PIL5s+fDx0dnSw+lUizMByowdN3T1FjZQ1hbFHnRZjVdpbCtV5PvPD9ye8RFKXY7nbZN8swpdUUjHQfia3+W4Vz/+v0vxy9Dpix5XJO+iJIpVLcjrytsMUBAHbV7WBa0hTVSlVD48qN0dqsNWqXq63yuogKSmJiIk6dOgVXV1ccPXoUCQkJ2c4pU6YM+vXrh7p16+Knn35Seo1EIoGdnR327t0LU1NTpdcQaRqGAzVwDXTF4IOfnmgua1AWL6e/zPTWemJqIqadniY8fPjRwUEH0c68HRquaYjID5HCObf+bnC0dFRpTRn7IpTQLYHYn2IzfRUxo7U312LC8QnC2JMfn8C8jLlK84k0yfv373H06FG4urri1KlTSEnJujpodrS1tVGuXDkcOHAg27oMRJqAdQ7U4PFb8anp5lWbZ7nnbqBjgNU9VmPpN0sVzg09NBTR8dHY1W+XQtnjoYeG4sDdAyqtqXnV5sJxfEo8wmLDVJr7Pvk9/rj4hzDWtU5XBgMqtEqWLAlnZ2ccOXIEkZGR2LRpE7p06QItrdz9V2ZaWhpev36NDh064O+//2alRtJ4DAdq8CHlg3Bc1rBsJleKpraaikWdFwlj8SnxGHpoKNqZt8NfXf4SzqWmp8JxvyPW3lyL7JgamaK0fmlh7PNnELIyy3MWIuLEp7+ntcrZGxNEmqps2bIYNWoUPD09ERERgVWrVqFNmzY5/pz09HSkp6djxowZ6NOnD969e5f3iyXKIwwHalBKr5RwnFWJ4YxmtpmJEbYjhLGbETfxz7V/ML31dExsNlE4lyZNw4TjEzDh2AQkpGS+hyqRSGBW2kwYi46PznY92wO2K2x32Ne2R+eanbOdS1TYmJqaYuLEibh8+TKePn2KJUuWoFatWtlPzODYsWOwsbGBn1/evhZMlFcYDtSgWulqwvGdV6qXXZVIJFjdfTXqla8njM+/MB+RHyLxT7d/ML7JeIV5a33XoumGprjy7Eqmn13GoIxw/C7xXZZr2Xl7J0a6jxTGDHUMsbLrSoUaCURFTfXq1VG+fHmEhobm+N/39PR0PHv2DM2aNcPSpUu5zUAah+FADRpVFmsKPIt5pvAcQlZK6JbAlt5bhLG45DgsubIEWhItrOmxBj+3/Vlh3t1Xd9F2S1u4HHJB8OtghfMZHz7MWN74o/iUeEw9PRUuh1yQLhXb4a7qvgr1KtRTOo+oKFmyZAlGjRoFqVSa61/uaWlpmD59OipXrox///0X0dHZ360jKggMB19IKpXi0tNL+O38bxhycAgG7RuEyacmY3vAdkS+j1Q6x6KCBSqWEMuy7ruzL0ff28qsFYbbDBfG/rv5H159eAWJRIIFnRdge5/tCo2TANlf/PVX10e/Pf1wKuQUUtJkT2InpiYK133elAkAElISsMF3AxquaYjlPssVPndaq2kYaTtSYZyoKJFKpZg5cyZmzpyZZ58ZGRmJH374AZUqVUL37t2xfft2lWouEOUXvsr4BW5H3sYo91HwfeGr9LyWRAt96vfBz21/Vqg2OPrIaKFXQc0yNfHw+4cKDZay8uTdE3z171dCZ8S/uvyFmW0+/ZfW/ej7GHxgcJYlj431jdG5VmccvHdQGJ/fYT5aVWuFe9H3cOX5FZwIPoH3ye+VfsaPLX7Ecvvl3E6gIm/KlClYsWJFvn+Pvr6+vCpjjx49FKoyEuUnhoNc8nrihe67ugsNi7Iy0nYkltsvl1cEvPr8KtpsFp943tJ7i8LDhtlxOeSCnbd3yo/rlKuDh5MeCr+kU9NTser6Kvx2/jfEJcfl6POzo6eth5VdV2J8U8XnHIiKorlz52LNmjVIS0uTv4GQnp6OtLQ0SKVS4Z/z6r9ejYyM0Lt3b3lVRj09lhun/MVwkAuR7yNhsdoCbxPf5mieubE5Dgw6gCZVmkAqlaL5xua4GXFTfr5KqSq4+93dHJUUvvLsCtpuaSuM+Y/zh00lG4VrX314haXeS7Hq+iqF1ylzw666Hdb1XAeLihZf/FlERZFUKkVaWpo8SGT85w8fPuDMmTM4dOgQzpw5g6SkzLuzflS2bFn069cPzs7O6NChA7S1VStURpQTDAe5kLGLIQDULlsbver2goGOAe6/vg/PR55KfwGX0C2BvQP2okfdHjgZfBLdd4stisc2Hot1vdapvBapVIra/9QWehvMaz8PczvMzXROTGIMdgfuxhb/LbgRcUPl7wIACSToXKszZraeiS61unAbgSiPxMXF4ciRI3Bzc8Pp06dVqspoamqKgQMHwtnZGS1btsx1kSaijBgOcighJQEmf5sIe++96vbCvoH7oK/zqX97bFIs1txYgwWXFijcytfV0sVhp8PoVqcbuu3qhtOPTgvnc1L2GACmnJqCFddWyI871OiA88PPqzQ3Ii4Cp0NOY9SRUUrPlzcsj6qlq6Jx5cZoVa0VHOo5oFLJSiqvjYhy7s2bNzh48CDc3Nxw/vx5oaFTZqpXrw5HR0c4OzvD1taWwZ2+CMNBDl18ehHtt7aXH2tJtPDkxycwMzZTev2zmGcYenAoLj27JIwb6Bjg0shLKG9YHpb/WSI+JV5+zlDHEOeHn0eLai0yfpxSGfsiGOoYIvbnWJUfblTWPvrssLNob95e5d4KRJQ/Xr58iX379sHNzQ1Xr15VaU7dunXh5OQEZ2dn1K9fP59XSEUR70Hl0L1X94RjKxOrTIMBAFQ3ro4zw87AxdpFGE9MTURvt94w0jPC6u6rhXMJqQnosbsH/F/6q7SmVtVaKcwPfataC2WpVIq5XuIWRM0yNRkMiDREpUqV8P333+PKlSt48uQJ/vrrLzRq1CjLOQ8fPsTvv/8OCwsL2NraYtGiRXjy5EnBLJiKBIaDHIpJEksdVy5VOds5etp62Npnq0JAiIiLwJijYzDcZjhGNxotnHud8BqdtnXC5WeXs/38ikYVUd6wvDCmal+E3YG7cS5UbM38Q4sfGAyINJC5uTlmzpyJW7du4f79+5g3bx7q1cu66FhAQAB+/vln1KxZE61atcLKlSvx4sWLAloxFVYMBzmUsXtiVv0KPqcl0cImh01ob95eGD/y4Ah2Be7Cmh5rYF/bXjj3NvEtOm/vjI23Nmb7SlTNsjWF45fvX2a7pnuv7mH8cfEVRLPSZhjbZGy2c4lIverVq4e5c+fi3r178Pf3x6xZs2BunnUnVB8fH0yePBlVq1ZFp06dsH79erx+/bqAVkyFCcNBDmV8GO/R20cqz9XV1sXegXthYmQijM/wnIGE1AQcGHQAHWt0FM4lpyVjzNExcDrghKgPUZl+djnDcsJxdn0RQt6E4OsdXysUNVpmvwwldFlshaiwkEgksLGxwaJFixAaGgpvb295tcXMSKVSnD9/HuPGjUOlSpXQo0cP7Nixg1UZSY7hIIesTKyE47DYMDyLeabyfBMjE6zrKb6q+PL9S6z0WQkjPSMcG3wMPb7qoTBv7529qLeqHv699i+SUhXfhTbQMRCOM5ZC/pzHIw+02NgC4XHhwvgI2xEY0GCAyj8LEWkWiUSCli1bYuXKlQgLC8O5c+cwduxYlCtXLtM5qampOHHiBIYNGwZTU1MMGDAA+/fvR0KCandFqWhiOMih+hXqK/yVfujeoRx9Ru96vdGtTjdhbLnPcsQmxaKEbgm4O7ljcovJCvPeJb7DD6d+QJ1/62DJlSXCnYTktGTh2s9fq/wo5E0Ihh0aBvud9niT8EY416xKM6zqtipHPwcRaS5tbW107NgR69atw4sXL3D8+HG4uLigZMmSmc5JTEzEgQMHMHDgQJiYmMDFxQXHjx9HcnJypnOoaOKrjLkw7NAw7Li9Q37csGJDBE4IzNF7xbcjb8NmrVjFcFW3VZjYfKL8+OC9gxh9ZHSmlRh1tHTQoUYHdKvTDdM8pgnn1vZYi5GNRuLh64e4+vwqDt0/BI9HHgpdFAHAtpItzg47qxB6iKjoSUhIwIkTJ+Dm5oZjx44hMTHzu4wflStXDv3794eTkxPat2/PqozFAMNBLng+8sQ3O78Rxtyd3OFQzyFHn9PHrQ/cH7jLj61MrBAwPkAIGRFxEZjuMR2uQa5ftuhM9K7XGzv67kAp/VL58vlEpLni4uLg7u4ur8qYmpqa7ZxKlSph0KBBcHJyQsuWLVlsqYhiOMiFdGk6GqxugAevH8jHLCpYwH+8v8LbDFnxeOQB+53iGwp3vruDBhUbKFx78elF/Hb+N1x4eiH3C/9MKb1SWNRlEcY3HQ8tCXeXiIq7169fC1UZVfnVYG5uDicnJzg5OcHGxoZBoQhhOMil3YG7MeTgEGHs9w6/Y077OSp/Rro0HdWXVxceDPyz45+Y3W52pnOuPLuCdb7rsPfOXiSlZd+kJaOSeiUxoekETG01lWWQiUipFy9eyKsyent7qzSnfv368qCQXe0F0nwMB7mUlp6G5hub49aLW/IxbYk2Lo68iNZmrVX+nEknJmH1jU8VEjvW6Ihzw89lMUMmLikOHo88cDLkJDb5bcryWnNjc7Ss1hJ96vdBz7o9UVIv8weSiIg+9+TJE+zZswdubm7w9/dXaY6trS2cnZ3h6OiYbe0F0kwMB1/AN8IXzTc2Fx7yq1yyMnxG+6C6cXWVPiNjX4QSuiXwbtY76GrrqjT/RdwLmK8wR0r6pw5uQ62HYpTtKFQ0qoiqpaqirGFZFX8iIqLM3b9/H25ubnB1dcXDh6pVYW3VqhWcnZ0xcODALGsvkGZhOPhC87zmYf6F+cKYRQULnB9+HqYlTbOdHx0fjYpLKgpjmT13oMx0j+lY6r1UfqyvrY/wqeEoX6J8FrOIiHJPKpUiICAArq6ucHNzw7Nn2dd60dLSQocOHeDs7Ix+/fplWXuB1I/h4AulpafBfqc9zoaeFcbrV6gPTxdPVCtdLdvPqLy0slDu+LDjYfSu3zuLGTL3o+/D+j9r4a7BmMZjsL7X+hz8BEREuSeVSuHj4wNXV1fs3bsXkZGR2c7R0dGBvb09nJ2d4eDggFKl+LaUpuFj6l9IW0sb+wbug0UFC2H8fvR9NN/QHNfDr2f7GXXK1RGOn8Y8zXZOQkoCHPc7CsFAV0sXM9vMVHHlRERfTiKRoFWrVvjnn38QHh6Os2fPYvTo0ShbNvPtzNTUVBw/fhxDhw6FiYkJBg4ciIMHD7IqowZhOMgDZQ3L4tTQU6hVtpYw/uL9C7Td3BZLrixBWnpapvMrlKggHGfXFyE5LRkD9w3E7cjbwviUllMUggYRUUHR1tZGp06dsGHDBrx8+RLHjh3DkCFDsq3KuH//fvTv3x+mpqYYNmwYTp48iZSUlEznUP5jOMgj1Y2r48KIC6hbvq4wnpKegplnZqLlppa4FnZN6dyMbw9kbIb0uej4aNjvtMfx4OPCuKWJJeZ2mJvL1RMR5S09PT306NEDO3fuRGRkJPbt24d+/fpBX1+xtPtHcXFx2LFjB7p3745KlSph3Lhx8PLyQlpa5n9cUf5gOMhD1UpXg/e33gqdFQHgZsRNtNzUEg6uDrj87LJQYCQ1XaxKpqyQUlp6Gnbd3oWGaxrC64mXcM5Y3xj7Bu5jN0Ui0kglSpTAgAEDcODAAURFRWH79u3o3r07dHR0Mp3z5s0brF+/Hh07doSZmRkmT54MHx8flYoz0ZfjA4n5ICUtBfO85mHh5YWQQvn/eS0qWGBgg4Ho/lV3zD43W3ig8WMhJKlUikdvH+Hw/cPYeGujUJHxIyNdI3i6eKKVWat8+3mIiPJDdHS0vCqjl5eXSr/4a9SoIS+2ZG1tzaqM+YThIB9dfHoR446Nw/3o+zme28asDR6+fohX8a8yvaZqqao47HQYTas0/ZJlEhGpXUREhLwqo4+Pj0pzLCws5EGhbt262U8glTEc5LPktGT8e+1fLLqyCNHx0Xn2uQ71HLCu5zqWQCaiIic0NFRelTEgIEClOY0bN4aTkxMcHR1RvbpqRegocwwHBeR98nus912PtTfXIvhNcK4/x8bUBvM6zEPver15O42Iiry7d+9iz549cHV1RXCwav/d2aZNGzg5OWHgwIEwNc2+GB0pYjgoYFKpFJefXYb7A3ccDz6u0pZDxRIV4VDPAc6WzuhUsxNDAREVO1KpFH5+fnBzc4ObmxueP3+e7RwtLS106tQJTk5O6NevX5a1F0jEcKBGKWkpKLe4nMKri3Pbz0WFEhVQpVQV2FayRc0yNRkIiIj+X3p6Ory9veHm5oa9e/ciKioq2zm6urro2rUrnJyc4ODgkGXtBWI4UKsdATsw7PAwYezW2FtoVLmRmlZERFS4pKamwsvLC25ubjhw4ADevXuX7RxDQ0P06tULTk5O6NatGwwMDPJ/oYUMw4GafEj+gPqr6yMsNkw+1sasDS6PuqzGVRERFV7Jyck4ffo03Nzc4O7ujg8fPmQ7p3Tp0ujbty+cnJzQuXNn6Oqq1hG3qGM4UJPRR0Zjk98mYeyI0xH0qtdLTSsiIio64uPjcezYMbi5ueHEiRNISkrKdk6FChUwYMAAODk5wc7ODlpaxbdOIMOBGqzwWYEpp6cIY11qdYHHUA8+W0BElMdiYmJw+PBhuLm5wdPTU6VyzFWqVMGgQYPg7OyMZs2aFbv/bmY4KEBSqRQLLy/E7HOzhfGSeiXhO9ZXoS8DERHlrejoaOzfvx9ubm64ePGiSlUZa9asCScnJwwbNgz169cvgFWqH8NBAQp+HQzrtdZITE0Uxnf124XBVoPVtCoiouIpPDwc+/btg6urK65fv57t9YsWLcKsWbMKYGXqx3BQwHYH7saQg0Pkx391+Qsz28xU44qIiOjx48fyYkuBgYFKrwkNDUWNGjUKdmFqwnCgBj+f+Rl/XfkLS75egmmtp6l7OURE9Jk7d+7Ig0JISAgAoGXLlvD29lbzygoOw4EapKWn4UbEDbSs1lLdSyEiokxIpVLcunULbm5usLKywrBhw7KfVEQwHBAREZGg+L7ESUREREoxHBAREZGA4YCIiIgEDAdEREQkYDggIiIiAcMBERERCRgOiIiISMBwQERERAKGAyIiIhIwHBAREZGA4YCIiIgEDAdEREQkYDggIiIiAcMBERERCRgOiIiISMBwQERERAKGAyIiIhIwHBAREZGA4YCIiIgEDAdEREQkYDggIiIiAcMBERERCRgOiIiISMBwQERERAKGAyIiIhIwHBAREZGA4YCIiIgEDAdEREQkYDggIiIiAcMBERERCRgOiIiISMBwQERERAKGAyIiIhIwHBAREZGA4YCIiIgEDAdEREQkYDggIiIiAcMBERERCRgOiIiISMBwQERERAKGAyIiIhIwHBAREZGA4YCIiIgEDAdEREQkYDggIiIiAcMBERERCRgOiIiISMBwQERERAKGAyIiIhIwHBAREZGA4YCIiIgEDAdEREQkYDggIiIiAcMBERERCRgOiIiISMBwQERERAKGAyIiIhIwHBAREZGA4YCIiIgEDAdEREQkYDggIiIiAcMBERERCRgOiIiISMBwQERERAKGAyIiIhIwHBAREZGA4YCIiIgEDAdEREQkYDggIiIiAcMBERERCRgOiIiISMBwQERERAKGAyIiIhIwHBAREZGA4YCIiIgEDAdEREQkYDggIiIiAcMBERERCRgOiIiISMBwQERERAKGAyIiIhIwHBAREZGA4YCIiIgEDAdEREQkYDggIiIiAcMBERERCRgOiIiISMBwQERERAKGAyIiIhIwHBAREZGA4YCIiIgEDAdEREQkYDggIiIiAcMBERERCRgOiIiISMBwQERERAKGAyIiIhIwHBAREZGA4YCIiIgEDAdEREQkYDggIiIiAcMBERERCRgOiIiISMBwQERERAKGAyIiIhIwHBAREZGA4YCIiIgEDAdEREQkYDggIiIiAcMBERERCRgOiIiISMBwQERERAKGAyIiIhIwHBAREZHg/wCmU00Scc/m/AAAAABJRU5ErkJggg==",
      "text/plain": [
       "<Figure size 500x500 with 1 Axes>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "dot:\n"
     ]
    },
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAgoAAADJCAQAAADoK9LtAAAAIGNIUk0AAHomAACAhAAA+gAAAIDoAAB1MAAA6mAAADqYAAAXcJy6UTwAAAACYktHRAD/h4/MvwAAAAlwSFlzAAAAZAAAAGQAD5bF3QAAAAd0SU1FB+cFGxMaHUUmVV8AABS2SURBVHja7d1/bBvnfQbwxyVXqjuLJ1MltVGmHClTNlIWIM0GEsMGJxSrHTezoQUJqnpwIFQFWjYLgrngXKGbA3soDFeoi8DpmAzVoNWoY89dJzht1aSoIRAu7A42rEEWhYWdVIsWC5El46N1G7mSyf54dSGPklPZJnX88Xz+iR0oypFxHr187nvvu+mD/8MmEBGtMON93DT6IoioYjR+zOgrIKLKwlAgIh2GAhHpMBSISIehQEQ6DAUi0mEoEJEOQ4GIdBgKRKTDUCAiHYYCEekwFIhIh6FARDoMBSLSYSgQkQ5DgYh0GApEpMNQICIdhgIR6TAUiEiHoUBEOgwFItJhKBCRjrm03y5kOePYqR5U7DmjXxhRrRu3Am/LQFPuZLSU37fEoWDP7lSvS6edgFfZqQ4lN/ItIqplcdPPpYT5ugRMbQYiFldmW3pLtiPTmX4iU9p/06YP0uU4ISpuuiT/VL4iuzI9y/uU/tTGvG1EtaMwBFImxVwYAp4Sx4BOY5lCQROy/Ei+IV2Ru9Uedb/iVcv4UoiqmmEhUKzcoaAJShPylDQtdas96uHEhr5EoopUMSFQbKNCQROU3rRNbY5Y9ig71GcUhgPVj4oNgWIbHQqacevbclBWzAcSvFdBtalqQqCYUaEgxE2X5OtSUAa8CutIqm5VGwLFjA0FTdw01nxD4r0KqiY1EwLFKiMUNCHLj+TLsqgjea+CKkvNhkCxygoFjXavgnUkGUmbGKzxEChWmaGgGbe+LYt7FZ9WWEdSueVDoHhisKZDoFhlh4Jm1KbVkbxXQaVTGAIzUpdapyFQrDpCQeDoND2a4hBoyj6Wbso9uWzLsb8qUE2hoBGj0z9u3q2wjqSPwhB4KNUYCpr86PSnFNaRxBAokWoOBU3h6PRgwp4DBtp3qP6Y0ddF5cYQKItaCAWNGJ0GvMq/2pMf/+zSmQVWkrWFIbAhaikUhLjpnz75d+0A0Jh9afHrvzb6emj9Qpa97gOJfInMEDBE7YUCMOL4246nlMwm6f34771n/tNl1pHVYtz6lx7xKymXQ88yQ8AQtRgKIUth6cg6snIFpaTpF5vvmn7VcNccsbgyTdmfy/fMwJ8nX73N/1KGqcVQWNvqOpI22uoQKF4J7HIDo3MMBEPVTyho8nUkpyPLrTAEZqQuFehZ/uiPA0GJHxQMV3+hIIidHN5q5nRkKa0VAsA+BeA7XEXqNRQ02sayrkyfwjryQTEEalK9h4KGG8uuB0OgLjAU9LixrP7dYAjUIYbC2uqzjmQIEBgKH63WN5ZlCNAaGArrUSt1JEOA1oGh8CBClrPNWh1ZDeHAEKCHwFB4GEFpQp6UK6+ODFnetYQb5izvmW83MAToITEUHoXR51wVhkDE4sowBKgEGAqPbuPqSIYAbQCGQumU/pwrhgAZgKFQesXnXAWlCfnIUuGHi5DljGOnOpRc/U8WhoCcteYYArThGArl8+1Pjjl+Pxf+RMrcu/zLT/wspBWSPtdbzS9HxxyDsWeUtVYCO9Xm7G61XgamqOIwFEovbjrm3Kf0p7q2H0xesnmVBUtb5mdbFiz9v7H/NvOxi/a29OP/+2vLbUvy493LO+9xJUAVhqFQSsPOP0oPJQfat2SD8sHkJdvMraD0Z91fWgxEgP2d/9H4202NuSdTL/xGhIDP1ZFZved0UPqFxJ2oyUAMhVIYtSXN/tiI47J817wtfUW+MxU39XbvUc7PAz4XEIgAcROg/0gQsjzfOXOr8O8EpaNtQFNWzp6fN/o1Ud1q/JjRV1DNRm3DzrgpKJ12zlm6tv+geTh6dfZ2gysD2HPHF243AMA+5XstAGDPFbcEnow1F5Tyvw9Khzq/sHR1diKsmEccRr82ql8MhYcQlHyukCVueqXtrqm3+03bkWgg4lWuyV4VGJ2LWADgoHJNDlmAJzL3zOPWtb/Ps4l//PB//nHroc5zYXFH4psLY461/71Gv3KqBwyFBxA3+VzjVuBQZ1Nur/uY80AiEHk5+nprcxY4EW1NhyyAJ5Myj9oAe64xe7YZ8GRa02/La38/f+yKHDcBwIjja23nwtrTFMVrCAAYtXVt97V3bRdfT1Q+DIV1GXaO2oCX2oCvtflcrszJ6Lnw660A4I89nUiYAXvOlTnjAIAu9acyAHwmMW8BgNfmPpe83/d9OdrnHrUNtI85JmcLH7B6NvGmrfDrRhynnRfDM7cGYwefMPq9oFpnNvoCKtuIA/DHhp3zlkn5unRFvjMVjz7e+1dLgFf90qL4msfS328eSgI9y0EZEaBPmZQB4MuxsAX46NuN/pgt+/3m5xLFxeJgorc7ENF+N9B+u2Fy1p4D/LExh/5cC6JS40phTaJCHLeOOW5IXdsn5WPRq7NBWVSIf7/wq5UK8a1mANip/mTlr6HNAHBk6dkEAHjXmFlcbSg5EV79dfbcHkUrGwfagauzWk05GDuzqm8YcXRt39rDxoFKg6Ggo68Q35aPRM/Pe5VrsicDXAyLCvEZ5SfNokJcbAhKwC5VnHo4lPzqAgDYc6WYMzgWfdUZlMatXduBwnWEPyaiSBM3DbRfli+Gz4UPdYYsRr9/VAsYCgBEhRiU9BXi4VhxhbjYMGoDPJnG7I9kUSFOyIAn84253SoAnIyW7oo8mXNhX/sbLV9fKP5osUcZ/bBviJv63FuyE2FPxqseXmMNQfTg6j4Uhp3jVlEhHurMV4g7VeBk9ClFqxDPNgPAU4qoEHcrcxYAOL6wXwEAf6wcTyp41ZlbE+HVjcSXY99pEb8KWfrcgzGteziypF9DCCOOXe5RG4jWrW5DYdg54gB8rinpVKvPdUUORG5Of6+lZ1lUiNclAOhZFiHQs3zJBgB9ihhIGo7uUwBgKGnElmxeVRSg49a97sFY/qOKPVe4hgCAoNS1fc5ydPG0k7FA61d3oaBViJPyDWmXe2rzd+cKK8T3zACwUxU/czsyF1rE74srRGMfYLr07phja88bLefC+vZiIPH9grVCUDrU+fWFQKQ/dTH8SpuR10vVpW5CISj5XHFTyKJViF9YOj+/LX1NtueAwLyoEJ9UL7TETUDnSoX4jFKOCvHR2XMzt+5MTYSL1yn9qWlJKxtHbYc6z618+PBkitcQwohj2Gn0a6HKU/OhkK8Qgd7uM47CCvHMQms6bgK86mLDuFUszMeaAa/amP2FVK4KsZwOx044ATHsdK4gNPRrCPEVW3vmLFOSz2X0NVOlqdlQiJtEhfhCh6gQu9VA5LU5rUL0LGsV4ukWAHhKEaPITydEhfituWfKWCGWz8kosL9za89lWT8h2Z+algoHpIedY46b04HIRDgo3+/JDKpXNRgKokI85hQV4rQUiLwz+72Wx9JAf+qzS0kzAHhX7iNsS0+uVIlTm4HCCrFapwbPz39x6Z3ZiXBxmB1InG7Rfj3Qfsk2uTIOdSR6vnmt7zNu5VMW9aqGQkFUiKO2KemyvL+zsEL0ZLQKcYf6qhMAmnI/bha/vyYDwInoF5YA4yvEUuhPrRVoL8XOOsT/5gPtwOSHE5JDSe2hrDyfa2vPGy19bsZCfaqBUCisEPvcP5WfS0yEAa1CTJkAoDO9UiGmFxtCFuBw4p45KAH+mFYhrmckuZp5Mi9H+9zDzq09W7Ln5wvXEQcSYwVrhbhpoP09883pibBXOcYasi5VcSjETT5XyAL42oHe7hPOA4lAxKtcaDmoAN+dE3sdedXQ5qAkHkq6JAP9qcbsuxbAk/nqgjsNVE+F+Oj8scEY8M5s/kEr4aVYfvcGMSEpQiMQeauZa4V6VIWhoO1qcPAJYK9bqxAvtOxTgEDEszzbANhzTyni559nWTyE/HRCDCR9a+6JDACcjFZXhVgK/tjJ6OqPFp6MNSfKRv2EZPEaAhDv/dYePnxV26oqFLQKETjV6nOlTKJC3JIVFWK4AQC6VbGb0bb0ZRkAvEqtVIjlc2rhrzuC0ohDPyGpX0MAQFDqczflbk4fX/C1G33NVD5VsZ/CqO2XDSejI44pCbghTUuTs0Bvt6gQX4yKrUzaM2MOfwxoz3zbiXlgh/oPbQDwUuyqCogZBFqbVz2+cLStZ/mdWX1gih2gtHcuKB3qfG2uPwUMJa9LI47KGOSi0qvolUJQGnZqFWLX9svyF5cmwopZMdtz9txrc6JCfHJZ3Ed4cjm0OWQB9iv3zCELMJh4MQoAnkytV4ilMJS8OhuIrF5D5XeRLJyQXL2GAICQxefa5eaWs9WvIkNB7GoAHOq8a+rtPuE8HAtEvMpPmvtTwHfnUua4CehPiQpxt3rPPGrLV4he9asLngxgz9VPhVg+/ti0NGqLm4ad+glJTwbQ794w4tjr7siMzr3qZN9Q7SoqFLQK8WibqBAPJAKR4wsXWp5cBgIRsauBPedKi/qrNf2mDbDnPMvFFaLRr6OWXAz/TUdv912TfkKyeAeoEceY4+a0P+bJvDZ3lA9fVbkKCQVRIYqNUYedKVMgci4szksYSj6d0CpE8cewLSMqxD0rFWJgXmyMygqxHDyZ1PU7U4FI8b2awUR+9wafa8yhjUP1p1Km1TtAhSw+Fx++qhYGh8KobdgJjDimpMvyQPsVORCZnD3rsOYAr/riyk/8x9I/aBZ/FX8MP6VMSwBw7MMpRJaIG8+eO5DwubRhp/yE5Fq7SPpcz3d2ZCZl9g3VwaBQyJ+tJCrE4ehE+HaDnBVnK2kVolgr7FOuyXET8LnkYgMADCYOxwBWiEYLRILyLndvd/GEZOEaAoib9ne+Z5655Y9devdVJ3eRrAYbHAr6jVH/0TG4UiGunK3UAABDSXEfYbcqzlbarQI/l0SFCLBCrBwzt44u3pkqnpAs3AEqbupzP5YWu0zac2vvIhmy8OGryrJBoaBViL72/MaoxxcutHSm8xWiJ2PNigVma/psM2DPibOV7LnvvMsKsTKt/fjYQOK0M25aPSG5ehfJuGmX+/nON1pe6DD6lVBe2UNBXyEC+bOV8hWidrZSt3pDAoA9ypQE5M9WYoVYXfpTXuWFjoH24glJe26PUrh7Q1Dqcz+bmLk1EZazrCErR9lCYdg57ASGnfOWy/JA+4+bRYXoVfRnK11e2c1A/ATpUcXGqF+OPZcAgP4UK8TqFIh8Stmh3pkqnnocSOR3bxB7SIqvOLNwds0Skh2EEUoeCtrGqJdsd03ibKWJ8O2GLlVUiOImYuHZSoUV4oOdrUSVzB9bawy6P3VFFsNN+gnJ1TtRA0Fpl/v5TutOBsNGK3EoxE2vtAG93eebtePZPZm1jmcXG6MeVIorRM7T17rX5nzt49aBdv2EJPDplb2wNCOOv3AfXZy59WL0BD9YbLASh8Il+UDiZPT4woWW33U8+4QM2HPVtjEqPar+1JHoGy071Jlb+o+G+h2gRhxjjv++2Z8CTkavyJWxVhi11cs9khKHwnWpIyMqxLWPZ9+taMezl/NsJapkQ8mJ8ForwgMJbaenwglJ4OXoVypicPo7LY/3imdyal2JH53uyNiywP2PZx+OJk3ARx/PTvXpRPTx3s8lgZNOOVs4IemP3ZD+sLcp25Q19vpC0j3z662vt0o5a3Zb2uh3q9CzidJ+7C5xKGgX15F5vRVhYKcqbj8eWRJhwbsJdD/23A9nD3W6Mqv/iJ+fH7eeagWOLhr5w2SX+5rcmP1M4tNKrdfgmz5I42bpv23I0tX776H+FMDNOKg0xq2nWpuy31wwamZlf2ePemSpDj7uNpYpFIARx2CiDt5A2lCjtlfa9ijHohxmK6PyhQJReYw4XnUyGMqosUL2UyBaL3/sztQO9flOn6tebhFuNIYCVSF/TGzey2AoB4YCVSV7LhC5OS2CwehrqTUMBapaWjBs7eGeTqXEUKCqZs8FIu/M3pAYDKXDUKCq58mcnxfBUPykJT0M3pKkmhGyDHUYPfdYAzinQLWlEgaiqxxDgWqP0QPRVY7DS1R7+lNXZ59L7HUPtNfDg86lx1CgmjSUvDO1Q2UwPAyGAtUsDkQ/HIYC1TQORD84hgLVOA5EPyiGAtUBDkQ/CIYC1QkORK8XQ4HqiDYQ3bW98Pg60uPwEtWhoHS0jXOP98GJRqpXHIi+D4YC1bNx69faulXu96jDMWeqZ/2pmVuceyzGUKA6J+Ye97o53qRhKBDBH7szxblHDUOBCACQn3us92BgKBCt0OYe+9z1Pd7EUCAqYM8FIhfD9T33yFAgKpLfCLY+5x45p0B0H3W6ESyHl4g+Sh3OPTIUiH6Xceup1m3pupl75EQj0e/Sn7o6u0Pd697h+WGj0deyEbhSIFqnz2/7lz94f9Mfq1uywNVZo6+mbLhSIFqvf779b7PAf0nX5D/5H6OvpZwYCkTr1p/6xhzw8fff3lLLc48MBaIH4I99afHzv/7Pmp57ZKdA9EDiJsCeA0KWE84r8stRf8zoKyoxdgpED8aes+eAWp575EqB6JHU3Nwjh5eIHl1NzT0yFIhKo2bmHtkpEJWGNvdY/fs9MhSISqY29ntkKBCVlD92Z6opV83bujEUiEruZLSa93tkKBCVgbbfY2939c09MhSIyqRaz7lmKBCVUX7ucdRm9LWsF+cUiDZAyPKVtrvmqhhv4vAS0UapkrlHhgLRRhq3nmptyn5zoYLnHjnRSLSR+lNXZ59LVPbcI0OBaIMNJcXcY6UGA0OByACVPBDNUCAyiD9WmXOPDAUiw+TnHn0uo68lj6FAZCgtGCpn7pGhQGS4yhqIZigQVYTK2QiWw0tEFcXwjWA50UhUeQwdiGYoEFUmwwaiOeZMVJmMG4hmKBBVLGMGohkKRBVt4weiGQpEFW9jB6IZCkRVoHAgutzBwFAgqhJaMPS5yzv3yFAgqiL2XCByMVzegWiGAlGVyQ9ElycYGApEVUgLhl3u0j8pwYlGoio2bj3V2rMciJTwW3LMmYh0OOZMRHoMBSLSYSgQkQ5DgYh0GApEpMNQICIdhgIR6TAUiEiHoUBEOgwFItJhKBCRDkOBiHQYCkSkw1AgIh2GAhHpMBSISIehQEQ6DAUi0mEoEJEOQ4GIdBgKRKTDUCAiHYYCEels+iCLD4y+CCKqGJv+HwG4PoNMMj/9AAAAJXRFWHRkYXRlOmNyZWF0ZQAyMDIzLTA1LTI3VDE5OjI2OjI5KzAwOjAwuv5ODAAAACV0RVh0ZGF0ZTptb2RpZnkAMjAyMy0wNS0yN1QxOToyNjoyOSswMDowMMuj9rAAAAAodEVYdGRhdGU6dGltZXN0YW1wADIwMjMtMDUtMjdUMTk6MjY6MjkrMDA6MDCcttdvAAAAE3RFWHRwZGY6VmVyc2lvbgBQREYtMS41UzZawQAAAEp0RVh0c2lnbmF0dXJlADViYWZhMWVjZjY0NGQyMDUxNWU5Y2Q5ZGZiNWRlNzhjZWZjNzNkOGEyMzYyMTRhZTQ5YjQ3YmRlOWUwZmIwMTnmWUEeAAAAAElFTkSuQmCC",
      "text/plain": [
       "<wand.image.Image: 5bafa1e 'PDF' (522x201)>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "asciipdf:\n"
     ]
    },
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAZ8AAAErCAMAAAA2Z37FAAAAIGNIUk0AAHomAACAhAAA+gAAAIDoAAB1MAAA6mAAADqYAAAXcJy6UTwAAAC0UExURf///wAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAP8AAP8AAP8AAP8AAP8AAP8AAP8AAP8AAP8AAP8AAP8AAAAAAAAAAP8AAP8AAP8AAAD/AAD/AAD/AAD/AAD/AAD/AAD/AAD/AAD/AAD/AAD/AAD/AAD/AAD/AAD/AAD/AAD/AAD/AAD/AAD/AAD/AAA2AAD/AABOAAD/AAD/AABVAP///+AterAAAAA5dFJOUwBVqhGIdyJmzJlEl2vSr8dEEVUiiKp3ZszSr2kzmZfHRHcimWa7Ee6J3fEziOOfp1XSl8zf10AwqoHM4fAAAAABYktHRACIBR1IAAAACXBIWXMAAABkAAAAZAAPlsXdAAAAB3RJTUUH5wUbExoe3C8E5QAABIhJREFUeNrt3Vtb2zYABuBA0gBhh3Jqd6BdCJSVJutY2QH2/3/YAo1DbGOQmR0p5H0vIlXRhRo9jV1/kdTpAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAACwLjZiD4BHbXa6sYdApc2N3qvN2IMgr9/LdDpb21uxh0O1zVe93k7sQVCpP73+9GMPgke4fwMAAAA6nYEAIWn9Xs8MJW06Qx6CJs0MpWIevu0upHBT3e1v8g0kpLu1M4g9Bqp0fbslzLUnbRvurgEAgJUm1UmbVCd1ntukzgwlp1+Z6uyKd1Ij1UmZVCdlrj1pk+oAAADQIrFb2sRuqfNgLXVmKHXdb78rpG6730vfUiF2S5nYLWWuPWkTuwEAACTl9V7sEcTSr/Porb3OD9hfqB8cHq7DDJUOp6gVzzXTOfiAjKPO4pRMZ+hgGR9RVIXDxQZ1PvCmOgcecHa0f/jmKNfyombo7ncF3XxROFxssFPjKKumOocfcPb2h7e3UzJz17T340+5P6+sQW9rszfo/Zwrjm8PF7tfD5edY1XeALFcaapzjQPOjt4cHr7Ltey9ffc69gfbkNvI4HhQKIqHi9VK5ZrqHHzA2cH0+rP4dbb3kr7dHlY6XGyjzpdWM51rHHBWuH976bNzq3TvFOEa9KwDzvbX4e76Qcd1/pfSXmcAACBB65sUrYY1SYpW2Jo8C1ph7c3QLK+Zl6WG4A4vWxYMvf8lX2bB0Dwpmr/x/vD/R0ez2OY4K0sNwR1if35xtZUUZUFNVpYagjusszVIilaYu4O0rW9SBAAAUQxPTk5Gi5V5Awk4PZu9ZJV5Awk4+XD7ev7rKKt8zBpy3epsXGiTwwZdfP2HcnqaVT5lDbmfttZZ7GSTw+aMb4Z35eS3WeVz1nBRWJpUZ2si2xg1ZZhNx82s8nvW8G9paZIZWr7RzeVdOfljVvmSNVzllkd97f3QeWFVWxh2t5/103/yPkzuiqtJVvlzVv5VWppUZ+NCmxw2ZHxxPn0dXp1nlb9n5T+FpUl1Ni60yWFzxtfDztn1+X1l3lC4f3PtiWQ8HOUr84Z7dTYutMkhAECAUnhTbBDzxFQKb4oNYp6YSuFNseFjWMyzXOFBURs9l6kU3hQbPlXHPPGEB0Vt9FyicTG8KaY5nx+JeWIKfxTURs+lGRbDm2KaU/nG5Cby0NdihkbF8KaY5nypjnnaGlO/MiDazadK1adhlSKlZ/RM40CtYnhTSnMq35hEHnl4UNRGz2UphjelNKfyjfOo4w4PitrouUSl8KbYUB3zxLMW155MKbwpNoTEPMsVHhS10RMAAKRsaZOyJS1wMRUhGk7yLjuPLqairrAkL3gRxnXnfJxoyraqQp7VBZ6fdX05ObkeJpuyraqnZij8/Kyzq7PHFlPF/oumqjolrF5yttgl9Pys65PJ5LR6MVXslG1VPZXkBZ+fNZpef0bVi6nipmyr6ukkr8b5WZe3LymmbKsqLMmru4g2vZRtVUnyAEjaf1ybhwaML1WSAAAAJXRFWHRkYXRlOmNyZWF0ZQAyMDIzLTA1LTI3VDE5OjI2OjMwKzAwOjAw48wLQQAAACV0RVh0ZGF0ZTptb2RpZnkAMjAyMy0wNS0yN1QxOToyNjozMCswMDowMJKRs/0AAAAodEVYdGRhdGU6dGltZXN0YW1wADIwMjMtMDUtMjdUMTk6MjY6MzArMDA6MDDFhJIiAAAAE3RFWHRwZGY6VmVyc2lvbgBQREYtMS41UzZawQAAAEp0RVh0c2lnbmF0dXJlADg2ODZlZGI5MmNjZjYzODkyNWQ4NzcyYWM4YTYxMWI5MmUzOWUwMjcxMTVhNjI5YWI0MmQxYWNjNWUxYjRlNDLgR1H9AAAAAElFTkSuQmCC",
      "text/plain": [
       "<wand.image.Image: 8686edb 'PDF' (415x299)>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "unicodepdf:\n"
     ]
    },
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAd0AAACSCAQAAACqTdjoAAAAIGNIUk0AAHomAACAhAAA+gAAAIDoAAB1MAAA6mAAADqYAAAXcJy6UTwAAAACYktHRAD/h4/MvwAAAAlwSFlzAAAAZAAAAGQAD5bF3QAAAAd0SU1FB+cFGxMaIvNAeGIAAAs5SURBVHja7Z3teeI4FIWP9tkGPCWwJXhLcEpgSmBKCCWQEpgSSAmmhFACKWEoQfsDJ4stgyz0ee3z8mcCHnSu5Qs24Y2UBiHEBbUD9DZ3ir9yByBEEmqtWlRYqVY1mZPwXZcQF1SDHYBf+pQ3x9+5dwQhwrjgAuCSOwbfdQlxQu1wBNDkvtpl6xIiEn5MRcgE1ErtVJU7xS281iXkLmp388MGQPZfCd1k4wkzIXZUhVofc6foJWLrEiIRXusSIhK2LiEiYesSIhK2LiEiYesSIhK2LiEiYesSIhK2LiEiYesSIhK2LiEiYesSIhK2LiEiYesSIhK2Lpk9piY/vMf2c4lQtSczxdTkh/coPP65JLF+pD76umTumJr88B7bzyXC1iVEJLzWJUQkbF1CRMLWJUQkbF1CRMLWJUQkbF1CRMLWJUQkbF1CRMLWJUQkbF1CRMLWJUQkbF1CRMLWJUQkbF1SOHYNXqIo7w9Ve1IkNi3eLs6XLcr7Q1+XFI5dg5coygfYL2xdQiTCa11CRMLWJUQkbF1CRMLWJUQkbF1CRMLWJUQkbF1CRMLWJUQkbF1CRMLWJUQkbF1CRMLWJUQkbF0SFdq2saCvSyJA2zY+lP5IVGjbxoKtS4hIeK1LiEjYuoSIhK1LiEjYuoSIhK1LiEjYuoSIhK1LiEgCta7aqA/VqlbVjv9vp9qHj+6mP5d9DDPlc7kfpw7D7RjDlHH2dujc80hdLkFaV21Q40W/YAuPb5/2p0itVYsKK9WqZnwL/5Qxcrs/7po7TerQe1tm6pIJ0LqqwgZbfQH0CSeswwTT73hDjRW2Yb4UZ6aMkzs0w5RMnS517jyPCfGu2+CISn2oD7VDhRVwfe1TtWqVvr4GqkodlFZaHbr3u1qdvx67y6W7jTIcwTqGmdK4J0hqC85jDFM6pw6R23kM59Qx9naAY6RktPcNOzRo0WhAY41WQwMNNA7YoEGtoYEWe1So0GKvoYEz9hpY4XzdvnuenfG8ze19t1sMR7CNYaYcueeJ1Gbukf1zW4PjGMOU7qmn7e3HuV3HcE/93DHiltr9GCn5FuIpWmzwoYEKf7DqdskKGuvvLWpoNBoaWEP3ft49PS29EexjmClH7nkitXPrOo4xTOmaOsbeto/hmrrUY6TkWxhfd4ULgBrv3b+gPxXw+f14BaD9MjJVde8k2OlsoT/ClDGMlMN7wqVWK2zwpkf+xxNjDHM7pdaX6bmfTT06hlvqKHs7yDFSLGFa91rkGhc0OPV21JUTgJebj5suauRJtFWuNra4/QzQPoaZciy3Y+rbVEN9/KqLj9TlNkY/Ze2Y+noAO6d+uLenjOGW+oljZOLe9j1GiiXEx1QXfOKiWpywRqXfzQ30BZ/YAIDaqX33fxpAVeE+xbOOYaa05HZPrbdfN7zhp/2laNIYw5TuqS25g6S27e0Ix0iEvW1NWRQBrnXXY1cFX1cR3U8VDtA4o+0+5qhxhsahdx3T2j4Y6G/RH8E2hplyLLd7anvu4eNuYwxTuqeOsbdtY7inLvMYKfkW5knO3xf/uxDPF6lUI6XM3Eyd8hgp9xbkD9yoFfaoABzHP+QoAzOlzNxMnS51yfBvUxEiEppDhIiErUuISNi6hIikmNaVanOWm5vGbxmpY+UupHVDeZJx3dlcxi/d2XCp53OMFNG6MjxJGr9M7Zo65jFSROtGM34tOI4RwPhN4c76G78p3Fl/4zfFMZLDCp+cO/d3QjTCGb9x3dkQxm8Kd9bf+E3hzvobvymOkRxW+NS9nb1tu/BBjN+47qy/8ZvCnfU3flO4syGM3xTHSHorfPreLmV93YjGb0B31tP4TeHOBjB+U7izAYzfFMdIeit8+jFSSusGMX6ju7O+xm8Kd9bf+E3hzoYyflMcI2mt8OnHSO6TZQ0NHLDGAS0217P87t6+ZnbGoTtxuJ4s/cFOA1X/Omb02au++NV7zGEMM+VY7sfPOD13qNRmSnvqOHvbbc9MSZ3iGHHbMymPkextq6GDGb86qjsbwvjVCdxZf+N32t7Ob/ymOEbc9kzKY+SJRotxk+FJ0vhl6nKOkUKkPxmeJI1fpnZNHS93Ia1LCHGjjG9TEUIcYesSIhK2LiEiYesSIpKFt66chZDD1rnUumdVee7fg2X9HdwGe1QaqNGiijZK/wsBFfb4A30deWyL+HWmqXvw9fkC6p7XjC/4XTeTvH1Ag5/4Bw12/k/2TJ1LrXtuM77g1h2RoDvFWb12ErRdih6K161aq1ZpdR4/IVM1Gmz1UX/i9/ihM8wwQSi3ZbCr9t5j2iq31x1h70/R3gXP+JJbt8YRe2z1v3jFETWAPY5aYXtdUgrAEQ1e8Y6XzuA54IIf+IHq+/XTvGeNX/iBE15vRjrdjAkcu/sq1YxsMcxgT2HLMKzTXvczYz6u3F53+L1v1jmvGY95rVH2zZC5VxhK6tOl6K972u77tnck6a/tNG6/hP7/FmYGB6H8TgabtB5iTFvltrrHKvfd+3btXfaMl+Lr5qEvQbvr4eY9Nr63UKvHj09NMSnDY2k9zphudSPK3rdp76JnfMknzAMJ+s40DKVo1d0uo/eMoLc38vYFUNcrnhrfJ03/b3H3ULifYkqGoew9pW7fMfuVW+q+W7nf3h/WPasZX3Lrji2W3XR/7GQEc2Hl0SWiH6JPOGKtKrXCK37fOdweZDDHnJDBvkx1+DHd6w6/96cshi54xpfcuu9o9E/9on/rf/QvAMAL1uoPdvh953/8Cyitzqi/tzDvMVBt789y/kSFPzjjhO3oFrYM5pi2DMM63et2H9Osy1Z3+L1v1jmrGc/+YVHWD6ruSNCxvyowKVvADFOl9RLqDpliuvZeQuWuGZb8rgu84FV9qA/UeAMAtVFnVeH6S4RMRMkwqLPEuqOkGNZdZOVPZ8j9WlPSDRX20NBxvx5YXoYS6l5u5c9m4F/JIEQkyz5hJkQsbF1CRMLWJUQkxbTujBRox0qXUvlS645Web5PFHufsmWQv8vQv5eqvXPGfSsv4l0324rl2fXvecnfxdc9q8qLaN3RFcsLlKAnrHFegPbuX3cC7Z0z7j3jZbSuqUCXKUEPM5SqvY/VXZb2zhn3n/GYZ/mTrwbGV/4uToJ2WOM8o/Y+Vndp2jtn3H/GS1Hth1J0kRK0dY1zau/3U9jWqOeMT6y8NNV+IEUXLEGba5xTe39Ge+eMe854Ga1rKtCiJWhq70/UzRkXqdqbCrRoCZra+xN1c8Zlqvb3V/6WKEG7VFqy/B1Te+eM+854Ge+6I1K0aAnapdIC5e8k2jtn3HfGc7++3X2tEixBSxuzjBQlVC6pbqr2hIiklBNmQogTbF1CRBKodfO4l2qnEvkf9+tcinW6VNu24BkPcqGdwzldo8UeB7Rf30UdbhFh/ExLLdO2vVf3cmc8yC+H8riX+h1vqLHCVif6VcLcl1qeXvlS6y5rxkOcMAdwL12XWgYAXLrbvR0f2r0MsNRyCts2wkLXEWzbHItrz2vGQ7RuCPfSdanl6+NveDO++RnPvQyx1HIK2zb8QtcxbNtQi2svd8YDnJV7u5fuSy0/uEKI5l76L7WcwraNsdB1eNs21OLay53xUL6up3v5hHlpJ4Z76bvUcgrbNo7xG9i2zWIZz2zGw/xyyNe9dHdOrVtEcS+H5qXrUssJbNsci2snmfHRupc742Fa19u9dDcvJ+YK6156L7WcyLbN4ZxGn3F/y3huMx6idUO4l+7OqX2L0O5liKWWU9i2OZzTFDM+UvdyZxyBvpKxEPdy3kstu1TOGc8/42GudZfiXs57qWWXyjnj+Wc80use3ctF1b3cyvPV/R9PBvROQbH3lwAAACV0RVh0ZGF0ZTpjcmVhdGUAMjAyMy0wNS0yN1QxOToyNjozMyswMDowMNIkEdwAAAAldEVYdGRhdGU6bW9kaWZ5ADIwMjMtMDUtMjdUMTk6MjY6MzMrMDA6MDCjealgAAAAKHRFWHRkYXRlOnRpbWVzdGFtcAAyMDIzLTA1LTI3VDE5OjI2OjMzKzAwOjAw9GyIvwAAABN0RVh0cGRmOlZlcnNpb24AUERGLTEuNVM2WsEAAABKdEVYdHNpZ25hdHVyZQBkZGE0NTdlYjQ2YmE4ZWMzNWQ5MmQ2OTAxZTkyNWE1NWY3OTViNTk5YzgzZjViYmQ3YTI3YTQ4OWIwY2M4NWE17eQPiwAAAABJRU5ErkJggg==",
      "text/plain": [
       "<wand.image.Image: dda457e 'PDF' (477x146)>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "events = pylhe.read_lhe(\"example.lhe\")\n",
    "for event in events:\n",
    "    fd = lhe_event_to_feynman(event).with_class(\"test\").with_style(\"layout : dot\")#.with_style_property(\"layout\", \"dot\")\n",
    "    fd.get_point(\"Leg1\").with_color(\"red\")\n",
    "    fd.get_point(\"Leg2\").with_color(\"green\")\n",
    "    print(FeynML(diagrams=[fd]).to_xml())\n",
    "    fd = feynman_adjust_points(fd,clear_vertices=False)\n",
    "    AllRender(fd).render()\n",
    "    break"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "id": "69e6d4f4",
   "metadata": {},
   "outputs": [],
   "source": []
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "id": "0b990896",
   "metadata": {},
   "outputs": [],
   "source": []
  }
 ],
 "metadata": {
  "kernelspec": {
   "display_name": "Python 3 (ipykernel)",
   "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.11.3"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 5
}