{
 "cells": [
  {
   "cell_type": "code",
   "execution_count": 1,
   "id": "e21d310f-957b-4161-9a7a-46f175c92e9a",
   "metadata": {},
   "outputs": [],
   "source": [
    "# imports to write dynamic markdown contents\n",
    "import os\n",
    "from IPython.display import display, Markdown, Latex\n",
    "from IPython.display import HTML"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 2,
   "id": "8d0c4639-5087-4647-a134-af26501a69fd",
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/markdown": [
       "*For the result HTML page:* "
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/html": [
       "<script>\n",
       "code_show=true; \n",
       "function code_toggle() {\n",
       " if (code_show){\n",
       " $('div.input').hide();\n",
       " } else {\n",
       " $('div.input').show();\n",
       " }\n",
       " code_show = !code_show\n",
       "} \n",
       "$( document ).ready(code_toggle);\n",
       "</script>\n",
       "<form action=\"javascript:code_toggle()\"><input type=\"submit\" value=\"Click here to toggle on/off the raw code.\"></form>"
      ],
      "text/plain": [
       "<IPython.core.display.HTML object>"
      ]
     },
     "execution_count": 2,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "# turn off/on code for the result HTML page\n",
    "\n",
    "display(Markdown('*For the result HTML page:* '))\n",
    "    \n",
    "HTML('''<script>\n",
    "code_show=true; \n",
    "function code_toggle() {\n",
    " if (code_show){\n",
    " $('div.input').hide();\n",
    " } else {\n",
    " $('div.input').show();\n",
    " }\n",
    " code_show = !code_show\n",
    "} \n",
    "$( document ).ready(code_toggle);\n",
    "</script>\n",
    "<form action=\"javascript:code_toggle()\"><input type=\"submit\" value=\"Click here to toggle on/off the raw code.\"></form>''')"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "c9f35205-d12d-40ef-a727-e085d92b12a0",
   "metadata": {},
   "source": [
    "ECCE style works with any ECCE software distribution via [Singularity container](https://github.com/ECCE-EIC/Singularity). Or use your own `ROOT` plus simply download the style files to local folder from [`ECCEStyle.C`](https://raw.githubusercontent.com/ECCE-EIC/macros/master/common/ECCEStyle.C) and [`ECCEStyle.h`](https://raw.githubusercontent.com/ECCE-EIC/macros/master/common/ECCEStyle.h). \n",
    "\n",
    "And here is example `ROOT` C++ macro to make an ECCE style plot:\n",
    "\n",
    "```c++\n",
    "\n",
    "{\n",
    "  TCanvas c;\n",
    "  TH1F h1(\"h1\",\"H1;Mass [GeV];Counts [/ GeV];\",20,0,20);\n",
    "  TH1F h2(\"h2\",\"H2;Mass [GeV];Counts [/ GeV];\",20,0,20);\n",
    "  for (int i = 0;i<1000;i++){h1.Fill(gRandom->Gaus(7,2));}\n",
    "  for (int i = 0;i<1000;i++){h1.Fill(gRandom->Exp(.5)*20);}\n",
    "  for (int i = 0;i<1000;i++){h2.Fill(gRandom->Exp(.5)*20);}\n",
    "  h1.Draw(\"e1\");\n",
    "  h2.Draw(\"same\");\n",
    "  h2.SetFillColor(5);\n",
    "  h1.Draw(\"e1same\");\n",
    "\n",
    "  TLegend leg(.55,.70,.85,.93);\n",
    "  leg.AddEntry(\"\",\"#it{#bf{ECCE}} Simulation\",\"\");\n",
    "  leg.AddEntry(\"\",\"e+p 18+275 GeV\",\"\");\n",
    "  leg.AddEntry(&h1,\"Data\",\"LP\");\n",
    "  leg.AddEntry(&h2,\"MC\",\"F\");\n",
    "  leg.Draw();\n",
    "  c.RedrawAxis();\n",
    "  c.Print(\"test_style.pdf\");\n",
    "  c.Print(\"test_style.png\");\n",
    "\n",
    "  TCanvas c2;\n",
    "  c2.SetRightMargin(.2);\n",
    "  TF2 f2(\"f2\",\"exp(-x/5)*sin(y)\",-5,5,-5,5);\n",
    "  //f2.SetContour(100);\n",
    "  f2.SetNpx(100);\n",
    "  f2.SetNpy(100);\n",
    "  f2.Draw(\"surf2z\");\n",
    "  f2.GetXaxis()->SetTitle(\"X\");\n",
    "  f2.GetYaxis()->SetTitle(\"Y\");\n",
    "  f2.GetZaxis()->SetTitle(\"Z\");\n",
    "  c2.Print(\"test_style_f2.pdf\");\n",
    "  c2.Print(\"test_style_f2.png\");\n",
    "}\n",
    "\n",
    "\n",
    "```"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 3,
   "id": "5c347263-b9ba-4cf7-b97c-6e1cc3ff133a",
   "metadata": {},
   "outputs": [],
   "source": [
    "os.environ[\"ROOT_INCLUDE_PATH\"] =  os.getenv(\"ROOT_INCLUDE_PATH\") + \":/cvmfs/eic.opensciencegrid.org/ecce/gcc-8.3/release/release_new/new/rootmacros\""
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 4,
   "id": "29941af0-c799-4bdc-bc46-e0e5ad918b9a",
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Welcome to JupyROOT 6.24/06\n"
     ]
    },
    {
     "data": {
      "text/markdown": [
       "via sPHENIX software distribution at `/cvmfs/sphenix.sdcc.bnl.gov/gcc-8.3/release/release_new/new.4`"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "import ROOT\n",
    "\n",
    "OFFLINE_MAIN = os.getenv(\"OFFLINE_MAIN\")\n",
    "if OFFLINE_MAIN is not None:\n",
    "    display(Markdown(f\"via sPHENIX software distribution at `{OFFLINE_MAIN}`\"))"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 5,
   "id": "c0107825-957d-4863-bbef-30634930f304",
   "metadata": {},
   "outputs": [],
   "source": [
    "%%cpp -d\n",
    "\n",
    "#include <ECCEStyle.C>\n",
    "// #include <sPhenixStyle.C>"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 6,
   "id": "91e6ed86-8173-44cb-b048-3946839f2edd",
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "ECCEStyle: Applying nominal settings.\n",
      "ECCEStyle: ROOT6 mode\n"
     ]
    }
   ],
   "source": [
    "%%cpp\n",
    "\n",
    "// SetsPhenixStyle();\n",
    "SetECCEStyle();\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 7,
   "id": "7c8023bf-9862-474a-949f-52275e8f54f6",
   "metadata": {},
   "outputs": [
    {
     "name": "stderr",
     "output_type": "stream",
     "text": [
      "Info in <TCanvas::Print>: pdf file test_style.pdf has been created\n",
      "Info in <TCanvas::Print>: png file test_style.png has been created\n"
     ]
    },
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAArgAAAHYCAIAAAApvgy/AAAABmJLR0QAAAAAAAD5Q7t/AAAgAElEQVR4nO3du27cSLrA8aKwmSBoBhg5ml1MvOpUp0dKyJ5U77HjC2CH5+wmZEc7qQ3MBfMcdmqyIwvqWA4Xg3McGE52BD0AT/CNasski81uVReL5P8Hw5CoarKaTbK+rmtUlqUCAABoctB3BgAAQLgIFAAAgBWBAgAAsCJQAAAAVgQKAADAikABAABYESgAAAArAgUAAGBFoAAAAKwIFAAAgBWBAgAAsCJQAAAAVgQKAADAikABAABYESgAAAArAgUAAGBFoAAAAKwIFAAAgBWBAgAAsCJQAAAAVgQKAADAikABAABY/cnJXpIkechrsyxzko3eRVHUdxYAAFNRlqWHo7gJFFarlZP9jEDHjy2Koi4pOybrNyWZdJWSTLpKSSaDPXT3lGSyS8ouyR7OWdNDHMfl9lwdPRyRRd/5AgAMj61M8VmsuKlRUA9rfRiTUUY/AIBetJQp3mIFN4FCnueeXwgAADxw0/RQFMVuNQpJklAVAQBAsNwECsvlMoqiMY1fAAAASqmuvSs37OXzlpI4jrMsm2BVQffeqj0ik66QSVfIpCtk0hUyaXJToyBDGNI0jeNYKbVarRaLhdQxFEXh5BAAAMC/vcQjWZYVRWFOrjCROgaCUFfIpCtk0hUy6QqZdMVbJvd7mHrEkKbpiDswcm25QiZdIZOukElXyKQrIwkUtErEEP4HsBuuLVfIpCtk0hUy6QqZdGVgfRQ2GmsVwhCFf/UrMukOmXRl50wWRbFVV61t05u2yuRDDvQQ4/64R2m/gYLMrxBF0WKxkOqEOI6ZZAnAsLRMoyvq34XMp5907m5MpmVZJmnM9PUB5zrZzjP7dsxYl105URRF4zvdlrcMT5CzKZxNEqgul0tzY5qmzLIAYHC6fO2uFLRJkjQulbdarRqri23pl8vlcrk00z+wDqAoisViIT+naSpbVquVLWPBkuZsRtX54ThQyLJs4vGBLaQd0B0IwGSWRjICvM4MFKLavDLq8yV2K+PGt0pfGU3WnpM6iRLqj2WJVMwDBV71W+kmr4LP8M6CqCbZYcnHOolMTdLE4GTnA+LqfAIIhy6Pu6yRaxbelWeg+Se90XxsbpV+hzeii9LGvzosEbYiuUrTdKtXycnZT44Gw9sZcDaFcyU+2Hn1BwAIiv7muvGZZn7NzfO8kr6xw4H+uSzLrdI7p2fM08cyD6d/lRlxos/n05Mt9b4ORVFIG0HlWJWd1+mOFFHT4gBZlsl5Nv/UuE8zt41/lY36cCxE0MxJuBHH8bbx4Ci5Op8YE4c3Gnqhn5Yba0l1StvzML23W/qt6jbqdI1Cl+reyld2fdxKCVKvTjZfJX+tH858C/UahcYmlUr6+uHqN9rGIs/2FnY7vf4pXw8Wnl8uefvYMCD1xxMGxCyW4ia6hDNTbrvnjpmpFGZ1G/dglsFpmrZEDI2BglnqN0YtlTQ7BAqN7SOV3dazV9ZuNElg5q2+pf6m2ltnQjOSQCHPczMcHr2hXF7waUDPHdQ1fmM26eebmbLLnneoHmjPybbH1Rq7lDUGCubDvLFMlZNQOSdbBQqNL6nstp698vMbzVbeV8KC+pvSex5EHztvD5Z9zaOgx+nK8B51P/yX5h8AA7JxAJ5uldcpN8YWonvXh4456agoirIs8zzXEYNexm/jIcwH+J56oWVZVta6a2z73iV9/YPQI0IrRzR/pXdd3V4ChSiK5B4wPyc5+8vlklgBwFDo4txWM6rLlcaJELrYIVCwffPb6qASMZRGXYieYiEE0hFSvnNue27lXNVPrGxh9oVtuZ9wSeKAOI7lw9DdgJMkKcsyiiJiBQCD0/2Lpi2l2T9/407M2YR0z3z51TaXw86k87+M16/M8dCLyvRT9Ykl4Jn7GgUdHzT+tbHmBwACZD6pugcKjc+3JEmkHbb+18b0MnmdOX/dtk0VdTLa0Pb49TZhUfvzX6IEs9uEw8H2tpoGtHMfKLTHfdT8ABiKrZ5U+ot+ZXZaZYz7b9xne/p69cPO5Zzk0Pam9vdY3mrP8sYfGBzYChqKnt24b3qI47glViCgAzAUXWoUpB1dfV5hHkVRmqZSjW/GAWa3rTRN9Z82pjdzslgsbJM3tzdqZFkmHcyTJKm/Hd3rvGUPTuzQE6IeS7XT9TGVtyMfEG3fW3M5hMLoF1OZIaTS18b5QQMx4reGnY37mh+3bR+h7b0H6sMg2/dsdp/sMpKiy0B0c0ymzKNgTshYn3WgktXG/JtbKs9/3ZwhW8yhFrbhkXq2A2l6qMxjUcme+ZYrmdFnTL/N+llqfFO2IZ0Bqmd+Xwfay07vP1R9Wdg+p5FpuYf7zhp6wwUwXBvLZlsxU2creGyzELas+2DTsWyz5bBxOoHKqWg8P/Wdt4c4Uva3zMzYeDb0r/Xd2jLTOI9ke/7L8AKF9g/dTx72ta6obcErPRpilIa1Tiv8kHuBC2NS9JiFxkr+h6d/OHNUxb4PusO7k5eY6W1b1KamhPoLR8NbibPHw8goWKWUdGFN7u3pcCEgUEAdgQKAfRhDoDBBBAqoI1AAsA/eSpx9TeEMAABGYF8TLqn7xR0aOT8oAADYB5fzKOhhxO2VIc7nH21Xnwm18tfGjjCNI5K3moEVAIAxcDV8oj461hxk0svYyMYZSSsZsI0Uqox4blyFrH5Eh+cTI3B9fX1+fi4XTBRF33777fX1dd+ZAjAS3kocN4fR5Wh9ozka1fM637rIz/PcNpGDnrUjr9FpzDlDzP00LrLu441hCK6vr09OTv7617/q4PL09PTk5IRYAYATAwsUGgvOxmkr3FZjtJCjVyoGdJG/VX7q764yYYiZ8mG5xnicn5+bUYKOFb799tu+swZgDLyVOG4GVzQOAJPZtvM8N+dOkH4MlY37IFmqH6iS1S5D1xrT2DY6OZ8YhH10y+X6AdDRkIZHSlfBjl0Uva0eKfmxhSONGWjZWH937euwAQAezuyNjr44G/VQL5KzLGtcIHXbdcB20x4KVOKVyldDsx7CttylVI04XCgdg9Meyx8cHNgSUG2AHcgQrZZSUz9vHz4Hbsuu2ves19LcakBZ437qJcVuKwAURSHrVdruO30sbkwbZ4FCx8+vx9hQXy6NIx1kkumiKFar1WKx8NA4gnGbz+d3d3c3Nzfmxtlsdnh42FeWMGhZlullrOt/MovVhxR7jbsyn4e2PAj93UmepTtkQH2+YLf5ZF6tVlEUbftwNr/1Nb5QSqUu63NOl5OeDt13JSn9L8zVOOShscNjZQRHfTHTltf2df4RIBn1cHp6qj/r2Wz21VdfXV1d9Z01DE+9I3b9T3r0ll6peeej6F3VB3nVB4iJSva6DChrpNt56yn1n7Z9X7bHuAjkUbxD8eEt224O0/G6bBxFuW/mbAodAxQzcX3VVNG4PYSrDeG4vr7+9ttv9eU3n8+JErAVXeq3FAyND7etnniVVzVONtP+eK8/D3d71LfEQ6K9yN9htxuPGLKBBQqlMWmBLUHL9AP7o++xrS4s81okUMADDfcxBGHOJrfb40u+mpdGwR/HcZeHUn3WuHqaxu27TVqz1SEqCRqrV7c6emmJVEy2oelmRNX4MdlKn90ij0AML1AwL2i5B3TNlRkR+/w8dq6nMoMA23XZeB9SHqCOQGG46uX0xq9DjWWzLo3qe+uemfZSvHuNgo5aGnfV+O7asyrvrmVX3eOV3W6WxjF3HetFBn17Di9QKJtqySp81iV0qS5rvL7LWhDQeCXZNu6cYYzVoJ9EE1cpcnTcYHuU1csnUZ/hvuxQC9uYmfr2eq7avyPJs7HjQctNz1I5euVdNwZYG4uA3VoB6m+/sY27cee274FDMchAQVRqEeI43rm+7iG6XHCNXStsvXJyZmbETggUBqqxebG9MNsYKNTTd/8G1XLcesG8sRW4yxHLbrFRfW/1CuaWLoqVVzU2K6Q17RmwNQ2rplUFBtruUA46UAhEYzyrVZLV14OohwWKtR6wEwKFgbJ9cDsHCh27Om2bn9J43HVp5+0eKGws3RurE8puA8psGasHCo1NNpKllioBWwbqlTq2/ISPQOFBbBVf9UChMWX9rqhfqY034aCvOezJ0B9Gk6WL3oqWsrM9UKi/ZKuqb9uFVH8itVcDdAkUzCdeSx3Attd2+w5ttTX558yTqV9i+5haMjz0doeSQMEzXbXVXgdYr/WqoDxAHYHCQLV/2ZBnRXuvLP257y9QaKmuN7e3Z7JxwEJLib5D5sXGqv4uN4t5MhsrG+rnv/G1uw0MCcrGc+WKs5kZB63jjKddZh4FMCZlaylbeW7IfIK26KE+M6BthvjuWuaYV8YsimaW6pk0Xy7z2ccdJkuW5+FWT0VzysX2Pbfs1tyDrAnQJbdmBmQmbNkPM/B24icemQjOJ+q40QbK9sG1VCuq1qYH2yjrh/RRsLXrt3S6bGl6eHh+2v+08Ut8x34MatOotNI+VFW/vPG8DYu3BwvPL5coD1BHoDBQjWVwewHTHihULoNtRwPaEjdubynyWwKF7vlpPw+NR+/4fiVN454be05Ufm3JgPmnxlcNDoHCIFEeoI5AYbh0odUyKqqSfmPhVBk81X1sXnsNgd5zbsxn07gfW6DQvQ942bm3QfdTV3+tMubu0y+XLeZ+KqPSNr59821uzEngBhYo6M/Pyd6Gq+M9hknhAhi0eoeD9lEA7aMeKrva6pm5sew37TCPwlYdAzeeio4Dymwa+3mkxmy5jSPYu7z9srXSIkDdP5T9iTbmo4s8z1erlSxIGsdxxxXHxyeK3JxPjIn0DnN4YTjfIdoVRaH7DLb3fbMtZCx96PI8lxWTJZnzbnQPz6dz3bP08Jc/8FhD5K3EcXwYM2JQSqVpOqmIgUABdQQKMAOFvvOC8fBW4hy43d1isciyTCor0jQtiiKKoiiKsizb2AAGAABC4zhQMMlYVR0xfPfdd1EUJUlCxAAAwFDsMVDQdMTw9u1bpdR33303qfYIABNHiwMGzU0Lx87ziz18YrKg0EcBHtBHAYAaXB+FxWLRcQbN+gupXQAAIFjOmh52CxQAAEDInC0KtVqtpEYUAACMho/OjO1G00EBAIDxcVOjQL8qAABGqf8aBQAIU2QhMzH3nTvAE4bzudTSS4PzDFcYHulNe7+r3aaol8ZW4gx01H4RDnKth4ljHgV4QKDgje1UZ1kmK9rssHwDHx9cGeqiUBNHoAAPKGm8aTnVOlbY9oPg44MrA5twCQAmRTc6VBoRsixLkqSxK0NRFPpVWZaZzRay7rN+FdPQISwl3OF8Yq+ur6/Pz8/lzo2i6Ntvv72+vu47U0GQ1WSca39IxnGslErTtJK+QidI07Tx8VvfzsMZXXi7SKhRAIZhvV5fXl7e3t7Kr2VZ3t3dXV5ertfrfjPWow8fPvzjH//48ssvkyT58ssv//73v3/48MHb0SvdEqUaII5j/XiVlXKlhUISlJ8HH/Jy3d3BLAAUHR4RjH0FCkVRmFd5lmVUqQEP8eLFi5OTk/fv3+stNzc3jx49ev78eY+56tGHDx++//77f//737///rtS6vfff7+9vf3b3/7mLVaodGOUJ5753EuSRGodNkrT1Nyb1DEQKCAQzqZwNiVJslqt4jjWi0OuViul1Gq1kh8IF4C6HSZBv7m52fjCcqT95n766ae//OUvv/zyi97y888/P3ny5Mcff/znP//pPz+N5bo88VokSTLWDwij4T5QyLJM7g0dIMuvUgu3WCyWyyWBAoBGD1wy5ueff1ZK/fDDD91fsnM5bfvGL8+3oig2RgmVvYmtXgV44D5QkJtH33tyz+iKtTRNl8uldPF1fmhg0NpLrIODA1sCvpL2QteY6l8Xi4WZQB53G/ej61yFtFYQLiAc7vsoVK7vyr3ErGTAbubz+enpaWXjbDabz+e95GdPuvfE/uKLLxr3cHx8vFWP7p2zWqk6lSghjmPdLTHLso19FHSUoF9ljqIEQrCXPgqmyr0EYDevXr26vLw8PT2VfglKqdls9vHjx9evX/ebsb48fvz49vZW2hq0p0+fHh0deTi6LsvNLz+6Y1Z38oSsxCt8lUJQ3NcoSARtdgA2Y2q5u4gbgG2dnZ29efPGLAUPDw9fv349shqF7p49e/bbb789efJEb3n69Om//vWvZ8+e7fvQ5hTOLcl27nDQpcEC8GerOrou6neOVKnlea4jBucHDYTP84zJ4orS/u///u9//ud/jo+PlVLHx8f//d///b//+78O99/+8DRnPtCJZXqlPM/NaZTqkzLp18pTUTdYmM9Pc0oGTJbnErw5D/vYqXmt12clq9xdY8LjGx4QKNTt6aliezo3FuH170hpmurnnn4SmjWsjUeRoMFzSYAh8nZ5+FvEqNKrcZRYFAoesKpQyGSIY5Ik5miIyhbdAltJY26svwqoGPbqke2jH0c8NpJAAR4QKABQg149MsuyymDiisViweAfAAAGwdnwSF11Jv/bQgGG/QAAMCDOKi4qk4u1y/N8lK0PND3AA5oeAKgh9lEwaxRWq5VthXX1eReekSFQgAcECgDUEAMFTaYimeZTjEABHhAoAFCDDhSmjEBh6AZRBg8ikwD2bcCjHgAAwGjsJVDIsixqtY+DAgAA59yvHlkUhV7RZOMSqwAAIGTuAwWZQSFNU2ZVArAt/zWO+27l1ZPHjHW0F0bPfdODzKYw2SiB1hZgZ0HdKboJtX1C+pY0SZJEUbS4J8mYdA5bCaER332gMPHmBtvqW33nCxiMsvT0r6OWqeRapqCNokheGMdxmqbyYFytVovFglgB3bUs6ugtD+4DBarXgL0i9PTPVrTbYghZ70aWli6KIsuyoijKspSZ6NpXwwFCs5dFoeI4jqKIqBnA0ElNQGPNgTzi6nWo8mUpjuP6q/SWyTbOYoj2MupBftBRc/1GIoYAMCCNNQdS2NfXy23vp5Wm6XK55BmIAdlLoFC5qbovFrUPehEK1dqguDGN4tsAMEm6aK+0qzY+2XQ1g60RNssyHiAYmJaOEkOX53n9/Uqroale4ZHneSVNfY2r+n7Kshz3+ZyC0d8Ugbs//147M9oyI3d9HMeSMflBk8dLmqbyg/6r+Spgr7w9rMY8hbPUB8ZxnOd5nudyAy+XSzOcz7JMvhaYaSoViXoKKXko6P1QeQhMR6X+QLc7VJLJY4E+3RgVP/GIf41xva5j0FvkV7MKQV5oVhjU01S+Q5gp3b0D9GDcN0X47s9/cDUK8rP5ENCvrTwNpIayscYRcMvbw2ovox7aF3rwM02E1AFU4v1KmN/Ymigv0bNQN75Wfu637wUAbyr1B/LoqLdIqvuHA9WNGJPRNj1IXG+rAJTbWG779kpC2wgo2cLjYDTW6/XFxYX8fHBwcH5+vl6v+80SghLHsf5uYGt36IjZWjEse6lRaKy7yPNcCtfGPobOyfQm9Y3yQ716wGQGAbYWR743jMl6vb68vLy9vZVfy7K8u7u7vLwkVoAmt7wEBxvnamxJwEMDg+OvRkEmOY/juK9ZyYqi0NOlyRbaDiBevHhxcnLy/v17veXm5ubRo0fPnz/vMVcIikQAeih1Y7uDqWV4pJr8VPcYFt9ND/pm839cHSV0rDA0axTa05g29s9gsahetH8K7969M6MEcXNzc3V1xccHTVofNrY7SKWpTmnSw6yYSgEV25YdPp8/vgMF/zX2sjqLdE7M89y8P9uDeslqew+G+l+37U2625sC4F/HXsxJkuhB1FEUyUIPWZYlSaIHWjN+EhU7DEbwlrd+ahS8SZLEXJ2Fm3Oy2u+3ltg8kBsVIdCPr43tDlmW6c5Yy+VysVgsl0uJMLpXagKBiJw/7MzpkOt/klvFzxM2SZL2w2VZtlwu6/etFBvyKunZEMdx5U3JzvM8N4OPKHJ/PuHB+fn53d3dzc2NuXE2mx0eHl5dXfWVqwm6v/W8HU6pPT+OzOdhkiR8V4FD/kqcHao72m2Mtf1MRdJlItWWeZNUbVKm9jR6484ZRo+ur69PTk5OT0/1VTqbzb766qurq6u+szYt97dVEBMuAYHzdvW6XxSqPWr2FlBLW2B7Z4jGFkepXTDDHenBZC4JY5tcAQN1dnb25s0bc4zD4eHh69ev5/N5j7maLDqJAkEZbVV5e49Q/a6l9UHdd1TWyzqYp8UcVymDPPVLKnEPTQ9DZ7Y6wb9eBpLwcWOgvJU4Pg5TX57VwxHbZ2sw37XuyqDVIwAdT2iNPZIIFIaOQAHAUAw+UJDhQJUCuN4lMBx6OueWmEZHBrZOywQKQ0egAGAohh0oVL5/m3Okq1E/hQkUho5AAcBQeCtx3M+jYDbhS4dJWXahvF+AlfnsAAAYir0sCqWamvmVUrLWg2JZFAAABsJ9oCCtDO0LohAoAAAwCO7nUZg4W8MKzd4AgG2F0FjvvkahvXFBDy5wftxA2Ga26jtfAIDhaZkw0Vse9tVHYbFY1McQ6hkLRhwoAEPEqtkAbNwHCpUlVqMoSpJEfpAoQS+qBgAAAreXZaZliVW9FIKeRIG1ngEAGJZxTuHcFyZcGrrJTrg02TcODNewZ2acLAKFoZtseTnZNw4M17BnZqxsSZIkyzLmTgAAYHBcxiN6UENln2Zv6sYZG0eDGgUMFDUKwOAMr0ZBD2rQfRg1s2PjYrGgagEAgKFwE4/Y6hJ2SzZc1ChgoKhRAAZnYJ0Zuz9lxv08IlDAQI37xgRGaXhND/UWh0YyFxOtDwAADIKDQEFK/Y5dFCXZiAOFyKLvfAEAhsdWpvgsVnyvHjniEEFQeQsAcKWlTPEWK/juoyD9Gcc6SJI+Chgo+igAgzO8PgqqQ21BURSsHgkEZb1eX1xcyM8HBwfn5+fr9brfLAEIiptAQRaEbJ8joSiKxWKh7vszAujder2+vLy8vb2VX8uyvLu7u7y8JFYAoDmruMiybLlcys9pmiZJovstCj0d04i7KdD0gGG5uLi4vb19//69ufH09PTo6Ojdu3d95QpAFwObR0HoOgObNE2zLHN1uAARKCA0++juxEUOhGCQgYKQ9Z+k/kDEcSzrQrk9UIAIFBAaAgVgrAYcKEwZgQKG5eDgwHbFciUDgfNW4vieRwFAOObz+d3d3c3NjblxNpsdHh72lSUAoXE5PBLAsLx69erTp0+np6d6y2w2+/jx48uXL3vMFYCgECgA03V2dvbmzZujoyO95fDw8PXr1/P5vMdcAQgKbeou0UcBA8XMjMDgDHJmRgAAMDJ0ZnTMNhqN72oAgG2FsPgwgYJjBAQAAFdCWD2SpgcAAGBFoADAvSiKQqgyBfBwBAoAAMCKQAEAAFgRKAAAACsCBQAAYEWgAAAArAgUAACAFYECAACwIlAAAABWTOEMgKnHAVgRKDjGolAAAFdCmOGUQMExAgKf5BbinAMYKxaFAjA26/X64uJCfj44ODg/P1+v1/1mCcBDECgAcGa9Xl9eXt7e3sqvZVne3d1dXl4SKwDDFVFt61AUcT69oukhNBcXF7e3t+/fvzc3np6eHh0dvXv3rq9cAaPkrcSZSsGWZZn+31QURX2jUipJkvp2vaXxJYpAwTsCBf/21CzKhwhsy1uJM5XOjMvlMo7j+vaiKFar1caXZ1m2XC7NvaVpagsXAAAYjUkECkmS2P5UFIVSKk3T9jQSJUgy+XW5XCZJ0vIqYJTav8EcHBzYElBnAAzUmKvKsywzKwziOJawwNSl7lrS5Hmuw4KiKBaLRX2HND14RtNDaM7Pz+/u7m5ubsyNs9ns8PDw6uqqr1wBo0QfBQeSJDGbFR4YKFTS2DaO+HwGiEAhNDLq4dGjRzpWmM1mHz9+fP369Xw+7zdvwMh4K3HGPDyyKIqyLMuyzPO8Y3rbxnr/BtnS+BJgss7Ozt68eXN0dKS3HB4eEiUAgzbmQGEjXcxHURRF0WKxkB/M4l9+rvdFkC0ECkDF2dmZHglZluXV1RVRAjBokw4UTHEcp2kq9QSLxYIIAAAANfFAQTcrlGUpEyoURZGmqVJqsViYadr3YIq2tIe3NQnMEwxgTLYtO3wWH5MOFLIskxChslF+sDU6mOp/Lbfk6r1MCvMEAxiZbcsOn8XHpAMFGzoqBu7FixcnJyfmPME3NzePHj16/vx5j7kCgFEiUGhg1hPYOi12qW/Aztor3N69e1dZTUApdXNzc3V1FUI1HQCMyaQDhSRJGssPMwiQ/+vTPMsWAgUAwLhNOlAQlcJeT+aot9dbImyTK8CV9pa5luqBENrzAGBMJh0oSHm/Wq1kBQcZ+CDjHcw5mqR742KxkGEROg2LQvVlPp+fnp5WNs5mM8brA4B7O/S0HBwp9WUYZOOfKvI8rySTMZOmNE3re5vI+ezd9fX1ycmJGSvMZrOvvvrq6uqq76zhD9N5vAB98XaLsTaBUkpJdYJSqn1BSF2FYKtLiFjrwZf1ev38+XO9ztB8Pn/58iU1CuGIWIYD2DNvJQ4Fm0sECp5RGgWLjwbYN28lzqT7KAAAgHZ/6jsDAEaIugRgNKhRAAAAVtQoOGYb4s8XLADAtkKYVZZAwTECAgCAKy1lircYgqYHAABgRaAAAACsCBQAAIAVgQIAALAiUAAAAFYECgAAwIrhkRgwBqMCwL5RowAAAKwIFAAAgBWBAgAAsCJQAAAAVgQKAADAilEPjrF6JADAFVaPHCECAgCAK6weCQAAgkagAAAArAgUAACAFYECAACwIlAAAABWBAoAAMCKQAEAAFgRKAAAACsCBQAAYEWgAAAArAgUAACAFWs9OMaiUAAAV1gUaoQICAAArrAoFMIVRVEIkSwAoF8ECgAAwIpAAQAAWBEoAAAAKwIFAABgRaAAAACsCBQAAIAVgQI8YbwlAAwRgQIAALAiUAAAAFYECqhar9cXFxfy88HBwfn5+Xq97jdLAIC+ECg4Fln0na+u1uv15eXl7e2t/FqW5d3d3eXlJbECAPhnK1/96rcAABj+SURBVFN8FisRixg5FEWDP58XFxe3t7fv3783N56enh4dHb179+4he5bLeujnBwAC4a3EGXzBFpRBBAr7iEO7vGsCBQBwyFuJQ9MDAACw+lPfGZiifr9btx/34ODAloDKAACYoKkEClmW6f/riqIoisJM2bKT9jRDN5/P7+7ubm5uzI2z2ezw8HDnfa7X6xcvXsjPBwcH8/n81atXZ2dnD8ooAMCPchqUUnEcN/4pjuPKOcnzvJImTdNKmjRNG4/SMTPBnvnr6+uTk5PT01P9Tmez2VdffXV1dfWQHf71r3/VOzw9PT05Obm+vnabcwCYFG/lyCT6KCRJYvtTlmWr1Uopled5nucSECwWCzNNURTL5VIplaapTrNcLnUlxJicnZ29efPm6OhIbzk8PHz9+vV8Pt9thy9evDg5OTGHUdzc3Dx69Oj58+cPzSsAwAM/8Ugv0jQ1awsaaxTkT2YVgsQBZoVBPU2e54077Hg+B3Hmu2eSKxMA/PP2nBzAcL6dJUkitQUijuNKHUBRFIvFor690tmwse+hbWOX8zmIgYLdM7mneT8CPz8A0C+GRzpQFIVEQ1IBUCd9ElsaJmQnSql6PwbZMsrWh21tCEXtYUT7C32+BQCAzZgDhY7qgYIZBMj/9TSyhUBho/l8bnaNFLPZbOdODwAAnyYdKJgNE9iTV69effr0qTKM4uPHjy9fvuwxVwCAjiYdKLQzaxTa05haFvCorOQhPztZbCnkdaecD6MAgPHpUnb0tSjUVCZcahTHcUulgjQuVHpENqZ5iP/6r/9qTzCC1vqzs7N3794NogsnAPRih2ejt1iBGgUAAGA16UDB1iFRqhB0jUJjGlsnR+e2bcsIqsIKADB0BAobeiFImnrrgxlMAAAwVgQK1SBAJlcwF3eoT5lgm1xBKVWW1n9ffNGcjePjtld1+bfx0JWUXTCZAUJDfRjQi0kHCuo+IIiiSBaQzLJMlnUw14eUnxeLRZZlkkYWg9h2DcnHj9WTJ9WNT5+qx48f8AYAANinMU/hrNmmahb1cQ15nlfaFHQAoaVpWg8UoihqOZ0fPqi//U198436+ec/tjx9qv71L/Xrr+rPf+76XhrJt6wun+R9yh4+dEY94IG4hABT5GsK50kECl3o6Zzbl5qs/FDRHigopT58UD/+qH74QSmljo/V48fq2bMHRQkfPqifflL//KdSSn3xxR87/Ppra3oCBQwXlxBgIlAYpI2Bwn0ypbpVALT78EF9/736y1/UL7/8seXJE/Xbb+rXX62xAoEChotLCDB5CxSm3kdh0H766bMoQSn188/qm2/Ujz/2lycAwLhQo+CS8xqFoXTx5iqCB9QoACZqFAAAQP8mvdZD+NqDxS+/VL//3rD9+Lh5+z4Mpc4DALAbahQci6Lmf/vAxAyYiPV6fXFxIT8fHBycn587WXYVCF8Ik/ETKDj28CkRu3v2TP3222exgkzM8OzZXg4H9GK9Xl9eXt7e3sqvZVne3d1dXl4SK2AKSjtveSBQGLCvv1a//qqOj//49fhYHR05mL4JCMqLFy9OTk7ev3+vt9zc3Dx69Oj58+c95gqYDkY9uOR5HoX97XD7Q3MVYXf7qETlmsToeRv1QGfGHvAEAwAMBYECgJ61fys6ODiwJaDaAPCAQAFA0Obz+d3d3c3NjblxNpsdHh72lSVgUujMCCBor169+vTp0+npqd4ym80+fvz48uXLHnMFTAeBAoCgnZ2dvXnz5ujoSG85PDx8/fr1fD7vMVfAdDDqwaWOox72cFylGPWAsWOtB8DEWg8AAKB/BAoAAMCKUQ9jQF0sAGBPCBQcs00xR1kOANiWz8WfbAgUHCMgAAC40tJd0VsMQR8FAABgRaAAAACsaHoAMAzMoAD0ghoFAABgRaAAAACsCBQAAIAVgQIAALAiUAAAAFYECgAAwIpAAQAAWBEoAAAAKyZccoxFoQAArrAo1AgREAAAXGFRKAAAEDQCBQAAYEWgAAAArAgUAACAFZ0Z4cA++tSwpjAAhIBAAZOwp+7BRDMARo9AAQ+yj4LSeZkewkBkABgoAgVMhduYhtgDwETQmREAAFgRKAAAACsCBQAAYEUfBcdYFAoA4EoIfbEJFBwjIAAAuMKiUAAAIGgECgAAwIpAAQAAWBEoAAAAq6l3ZsyyrCiKxu1JkphbiqLQKbMs23fGEEJfXwDA1AOFoihWq1XjdjNQSJLETLZcLvM8r0QSAACMTzTx5e/ka2ue55XtZhCQZdlyudTJiqKQX+unLoqiaZ/OCZH6jonfPgB6FEWeSnAChUhtetzrYEJHDxI6pGlaaYMgUJgOAoUR6HL7A8EiUPAkiqI4jhu7KYiiKBaLRT1N4yOGQGE6CBRGgEABg+YtUJj0qAcp+3U9ga1Xo/q8JQLYk2g/+n5bAIaNQEEVRSHP08ViEUVRY0xQ3xjHsbLEFsAOKNEBhGnSox6kmJfhDGmaKqWWy+VqtTLrcxrHRAB74rYekdgDwMNNukZBpGlalmWWZVmWlWUpVQUd2xrqNQpRtN0/AP6t1+uLiwv5+eDg4Pz8fL1e95slTFzIrYrUKDRsjKJIVyTEcdxSqVCPJ+gXBQRuvV5fXl6enJzIr2VZ3t3dXV5evnnz5uzsrN+8YbJ26JboLVagRsGK/gfAKL148eLk5OT9+/d6y83NzaNHj54/f95jroBgTbpGoQuZk7EyUaO677vAaAggQDt807q5udn4QgZSYpomXaNgG+Mg5E/yP7ULAIBpmnSgoJSS2gJzi0ycIF0a1X2gUOmmIGlkoASA0JStWqoN2l/o8y0A4Zj0zIwy66JSKk1TCQiyLJOYwDwtrPWAOuczM0Z/zBLoan+yT6WoMK85Pz+/u7uTtgZtNpsdHh5eXV31lStgWxFTOPuhgwBTfWXIyuqRjWkUgcKUECgMl4x6ePTokY4VZrPZx48fX79+PZ/P+81bXcQ807AgUPCqKAo9nXNLrwU9nbMtDYHCdEw2UNjTiCzPD6L1ev38+XNdfzCfz1++fBlglKAIFGBHoDBIBArTMaBAYRB6eRANogweRCbRC2+BAsMjAWxnstEMME0ECsDuAl/JiW+hAB5u6sMjASBMLEiBQBAoOMbiTxNRlnv5h6FztVqPDM24vb2VX/WCFMQKUxPColB0ZnSJzozAtkY2htNV38OLi4vb21tzQQql1Onp6dHR0bt37x64c4wDox4GiUAB2NZkA4VxDDRFjxj1AGBCnJealJeAK/RRAAAHtu17yIIUGAoCBQB9GkeHUOd9D+fz+enpaWXjbDYLc/pIjBt9FFyijwLQu146PTjvezisBSnQC/ooAEBAduhFIWV8eyOC7U9nZ2dv3rx5/vy53nJ4eEiUgF5Qo+ASNQpA7/ZUo7CPQQrdB0fwoEYdNQoAEJD2J/LBwYEtAWU8ho5AAQAeaj6f393d6f4EYjabHR4e9pUlwBVGPQDAQ7169erTp0/mOAXpe/jy5csecwU4QaAAAA8lfQ+Pjo70Fld9D5kdAb2jM6NLrd2bfWYEmK5+54Sm7yHcau9FS2fGQeL5AABwpSUU8LaAJE0PAIAH8bzqMTyjRgHACPVbbnU8Oi0UGARqFAAAO9p2KSwMEYECgFFxvsrUPv6Ng/OlsBAmRj24xBTOALrYX8OIz0e686Ww9mp8A1K8TeFMoOASgQKALoYSKAyli+I0V83wFijQ9AAAvvXYnDHBEQp0pHggAgUAmIQdysuyVfsK2o3u/+ovPKIjxcPR9OASTQ8A+rKnaoKWMuL8/Ny2FNbV1VXjS+7r/51lb+NEnCPuSEEfhUEiUADQF/+BgnxZf/TokY4VZCmslkUu9hQo7EMvhWOYgQJNDwAwBv7Ltf0thTU+UTc7JPaReWoUHGJRKAA92kfZ4XZAgf+mh4ODA9tfH1j8bfuu94FFoQaJgABAX1qeP19+qX7/vWH78XHzdrXPWn2f5vO5rSOF55y0fDr/+If697/VL798tvHJE3V8rP75T+urvH1A1Ci4RB8FAGH6+9/V7a36+efPNj59qo6O1A8/NL+k+4LdIdco7NCRYqP1ev3ixQvpCxlF0Xw+f/Xq1dnZmT2Tw65RoI8CAIzfs2fqt9/Ukyf/2fL0qfrXv9SzZxte6LZlvWNuo8jZ12XnHSkmON6SGgWXqFEAEKwPH9SPP/5Rf3B8rB4/Vs+eqT//2Zp+nwMKrH/68EH99NMf9e1ffPFHJr/+2pp+H9Ue7XYYb7mxHmWHhiGlVBR5qlEgUHCJQAFA4O5L1r7z0eTDB/X99+ovf/lPa/2TJ+q339Svv1pjhT0FCvtoLGg58g4NQ8pjoEDTAwAgCD/99FmUoJT6+Wf1zTfqxx8d7NycGjI0OzcM+UGNgkvUKABAu77GcLqyw3jLLl04t20YUjQ9DBSBAgC0G3qgsNeJq7dqGKLpAQAwQi0rPH3xRfNLjo8fumamQ69evfr06dPp6aneIuMtX7582UNuvCBQAAAE4fHjz9rpxdOn6vHjDS/sOOfxVmzH2uvE1X1FP+1oenCJpgcA2NmHD+pvf1PffPOf/v/Sp+/XX62t9QOaPtJ56UDTAwBgWr7+Wv36qzo+/uPX42N1dNQWJajWhoyd/6GCGgWXqFEAgIcLebKHcHirUWBRKMds9WBc8QDQEQ9MLYS2FQIFx7i+AQCutJQp3mII+ihMTgjx6UZk0hUy6QqZdIVMDg6BAgAAsCJQAAAAVvRR6CrP89VqpZSK43ixWPSdHQAAfCBQ6KQ+SxfDSgEAU0DTw2ZJkiil4jh++/bt27dv4zhWD16tvOOrux9kH11vyKQrZNIVMunEPg5NJv2n9IYJlzbI8/y7775Tn1chSJTw9u3bShtE9wmXoqjrSmJud7iPlGTSVUoy6SolmQz20N1TksluKZnCOQDL5VIplaapuVF+lT8BADBi1Chs0FJ5oGo9FahRCDYlmXSVkky6SjnNQ3dPSSa7pfRRghMobNAYENi2EygEm5JMukpJJl2lnOahu6ckk91S0vQAAAB6RY3CBtvWKHjLGABg4lg9cniIugAAI0PTAwAAsCJQ6CTP876zAABADwgUNpApE2SVBy3LMqWUTNEIAMCI0Zlxs3q/xfrkCkVRFEWhlEqSRKZ8DpDOpLqPdUJWFEWWZTrDQRnEmQz2mpQzZjtvgZzb7pns8fS2Z1KTWynLsl7yuTGT+k/Bnsner8kuGdhvJktsIjUHlbUezFMX/lltbDpJ07TvfLUJ80yW99dDyGey8ePO87zvfP1B7qbGP9XPbV/Zbslk/dzaUu5bx0MHeyYrM96GeSb7fXJ2fHTv+8YJ8UEcoPrHUPlTHMd5nud5Xg8jQqBvQsmkvj9DK+G0+nkOhHnqzDMZTjFc2j/uvvNVlsb9Uv+TeTL7zXZLJs1zWxrPcf8lXEsm68n6ukRbMqlPndxK5q/hZNJ2K3k7mV0e3R5unCCeHYPw9u3bNE3TNH379q3eqC9uM2VoJYdcN5XboDHngTC/Z/Sdl6p6rhpPb48a8yMbe4wL0zQ1o+2Wh7J543jO9sZM2u4an7d8lzOpmd9HfT6Rdvu4S7/1iLt93LLRz/3e8dHt4cYJ7kE8LHKdVT6P0EoO22MizJJY34cBZs92+0kE2UOGmjRekz6fbo0qdXK2h3J9u8/LYGMmbRdA4znvK5MmncZzoLDzx+3zVur4cfd4TXZ5dPu5cRj18CAyGqLSAUc6klQGSvRI7gdbL6HQegtK/9DQciV038DKdukp5j07zSR7lRPY+/ksisJ8rtXJCey30+XGTIpBZFLd57OXod0bMxnCrbQxk5K9ypPc563U5dHt58YhUHAgqC7ldXI/1DfKD0FlvsdHWxc6LiyKQnpoBzguQ8ep+oFbFIUsiR5ONGNTvxrlWRnISc6yrCzLeiYbvzD0K8uy1WqVpmlQudLMQEGGYwR4K+lTp38oikK+yfh5RnV/dO/7xmEK5ynSl3u913GPAn+0aVmWSbmrlFqtVsvlMo7joJ5xZVlGUbRcLnU+lVJ5nod8YsOpgduWnNU4jsM5vRIaxnEceGioH0Tq/lZSgU2EL3HharUy1/Hp8VaqP7r93DjUKExOlmX6UgvnOTKUR5tSSvJpdno3v76HQD/F4jjWDbFB5XBbQcVhWpZlURStVqvQIsWQ2++EFG/1W0kFWTEjP+tbSQc3/jOz7aPb2TXgqrPDNNnOYZjntq8u0F3Uz1iA51CyFPj4kcYeWB2H0nlg63tl63DnvyNe2aHvp3kr9dWPtfuZ7OUcivaudvW7ppenky2T+lOuDyjwfL+3PLr93Dg0PUyFVKCpwCoShJ4Su54x3VUnnO8ZlUyGkzEh9bf1zozy3befPI2LbnsKrSJByKdsztMnW6QrQCC3UhzH0s7YuF36APWRr8/InV5paJC+FD4zGcKjm0DBgUAu6xb6UitDav+rWK1W9ZJMt7KHc4Yb+w2F83RrF2wm5RKtZy+0foI6Sgi8z4fZPUXI/WUGEL2rn0D9pAqB7dqzXav7sPHR7efGoY/Cg0hEXLnxQlsySprZpGKt77w0y7Isr5E/yc+BVIHYOhKHVpi1CDaTjaM6Q6PHj5RNYx8CUb+V9DQP4dxKko36x20bNtkjWyY96PLo9nTjOGnAmLL6aZQt4XQCGOgHHWC2G7sjBDVBcnl/3honXAohky3N/7Zz678TwLYdKXrRfRKtAPsolE0fd19XafvH3fuESzskc37jUKPwUHpODKnTqw+9DURk0Xe+hiRJEvm4oyiSpsokSXRFdN+5+4PkZLlc6mtSd5YOJ5ON5OkWRZHOdmjTP0jV0WKxaLyVAq8OCY3+uOVWMrv09521P8gHulqtzFtJnpk+M7nx0e3jxnEVcUxZvZWh7xz9x8ayoe8MWgWbvfrHHchXTK3xQw8kk+3fgwM5txv76jcKpK9+XZg1CmUwC7G2ZDI31vnznMmtHt37vnGiMtR262HRXYQC6VSMvRrExz2ITDYKcKgL9if8j3sot9L+ziSBAgAAsKKPAgAAsCJQAAAAVgQKAADAikABAABYESgAAAArAgUAAGBFoAAAAKwIFAAAgBWBAgAAsCJQAAAAVgQKAADAikABAABYESgAAAArAgUAAGBFoAAAAKwIFAAAgBWBAgAAsCJQAAAAVgQKAADAikABAABY/anvDAAAMEhFUcgPSZL0mY89o0YBALB3WZZFURRFUUuZWhTFxjTeRAYdEIiiKJIkiaJocU/yXEm21YG6JJCDiizLdjjWbggUAAD+rFYr2598Fn4dpWmapqkZuGRZtlgs5F3EcZymaRzHSqnVarVYLHZ4C/Jy2wtlu6RJkkQfzqsSAIA9S9NUlzt5njem0QniOPaaOXtmKhv1u0jT1PYn27uz0S9syYa5T0lfz8D+UKMAAPCk5duz1Nv38HV5G8vlUimV53n9LWRZJkX4tpUKXdL32xZDoAAA8Kqx9UHKS1sMoZvnkyTpkqbeXWBjgo10K4Ct2JYE9Xe38dC2+Mlsd+iTt7oLAMBkybdtadRXTfXzUiTlea4+b3ow2yxshVdjGrNyfmOCuvpRdA53eOMVlbYVedf1ErnxiDQ9AADGrLHmQL5kN5apurZfl1v179+VNFLuysaOCbrbthWg8dCr1cqsV9D7bKzn6H0MCPMoAAB8q9TP6+ihPhBRKdU47kCn1J0bdBoZHdA9wQPVmwySJJFjyf95npuHzvNcxkeYGYjjWKIHnTKUdgdF0wMAYP9004P+2awk0OVRvemhfVelUW9vaxTYmKBRpYi0tQ7o7SbdLmArZ+vb6/uXEKGeZ5oeAAAjV2l9aGl30IqiyLJMugRWmgySJJEy1Tbx0cYEXdjq/6WGQGusAEhqWvav8yaVLr23OyhFjQIAYP8q1QBSoNZ/rtcoVIpe3R2yUutQn4mo8p17Y4KKehEpW9qrJcyv+42VDS3lr/naxvdYT+YHNQoAAN/ki7JtPKGWZZn8VRfPZit+JWVRFGVZ6q/1y+XS7D2wMcFG7VMoisYuirYCuP4W1H3nR9lPIFNVEigAAHyTIrAoivZ2Bz1koH2FCLNAlZYF+TYvO9+YYKs8VwYsVHLSEvRUdtUSBOj9BNHuQKAAAOiF9PNvmWdJq5SX9aGVy+WyPZJoT9CR2dehcUKnxWJR2dg4V2OWZcvlsjHakPSyn/ZOG155aeAAAExavdHdLAj1xkofBf0qaXqo9BbUr5Jf0zSVHoX1ZRc2JqizFZE6A9JbQvamN9Y7EOjEHQ+t350tgf8+CgQKAIC9a+ydpwtXvaXembHy5Vb+VIkV/MzM2L43nef6nuuJWwIUHXPYEvgPFKKy6T0AABAI3ZVBjy3UW8xa/XqyjftpEUWRUqqliNR7Ux16HW516HbScpGmqbeujgQKAABUbQwU+uI/UKAzIwAAsGKtBwAAmgVVr5AkScfhl24RKAAAUBXEakyfMzs3+JxigT4KAADAij4KAADAikABAABYESgAAAArAgUAAGBFoAAAAKwIFAAAgBWBAgAAsCJQAAAAVgQKAADAikABAABY/T8aJeWOMDBfRAAAAABJRU5ErkJggg==\n",
      "text/plain": [
       "<IPython.core.display.Image object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "%%cpp\n",
    "\n",
    "\n",
    "  TCanvas c;\n",
    "  TH1F h1(\"h1\",\"H1;Mass [GeV];Counts [/ GeV];\",20,0,20);\n",
    "  TH1F h2(\"h2\",\"H2;Mass [GeV];Counts [/ GeV];\",20,0,20);\n",
    "  for (int i = 0;i<1000;i++){h1.Fill(gRandom->Gaus(7,2));}\n",
    "  for (int i = 0;i<1000;i++){h1.Fill(gRandom->Exp(.5)*20);}\n",
    "  for (int i = 0;i<1000;i++){h2.Fill(gRandom->Exp(.5)*20);}\n",
    "  h1.Draw(\"e1\");\n",
    "  h2.Draw(\"same\");\n",
    "  h2.SetFillColor(5);\n",
    "  h1.Draw(\"e1same\");\n",
    "\n",
    "  TLegend leg(.55,.70,.85,.93);\n",
    "  leg.AddEntry(\"\",\"#it{#bf{ECCE}} Simulation\",\"\");\n",
    "  leg.AddEntry(\"\",\"e+p 18+275 GeV\",\"\");\n",
    "  leg.AddEntry(&h1,\"Data\",\"LP\");\n",
    "  leg.AddEntry(&h2,\"MC\",\"F\");\n",
    "  leg.Draw();\n",
    "  c.RedrawAxis();\n",
    "  c.Print(\"test_style.pdf\");\n",
    "  c.Print(\"test_style.png\");\n",
    "  c.Draw();"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 8,
   "id": "d356f857-80c5-4f5a-9924-c595389c0bdb",
   "metadata": {},
   "outputs": [
    {
     "name": "stderr",
     "output_type": "stream",
     "text": [
      "Info in <TCanvas::Print>: pdf file test_style_f2.pdf has been created\n",
      "Info in <TCanvas::Print>: png file test_style_f2.png has been created\n"
     ]
    },
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAArgAAAHYCAIAAAApvgy/AAAABmJLR0QAAAAAAAD5Q7t/AAAgAElEQVR4nO29ebgkR3nm++VZutUS2loLI8kIDMgCJBYDw2I8VJXUghFmAJsBL9gXDxg/slnGcH3ta8PjU+UB+8LlYhaD9WCMYViuQWY1FwHd6soSZjBc0LAJW9eImZHHCIxarYVWd5+qOnX/iHPiZGVGREZGRmR+mfn+nn6kOlm5RFVlZrz5fl98ES0WCwIAAAAAULFSdwMAAAAAwBcIBQAAAABogVAAAAAAgBYIBQAAAABogVAAAAAAgBYIBQAAAABogVAAAAAAgBYIBQAAAABogVAAAAAAgBYIBQAAAABogVAAAAAAgBYIBQAAAABogVAAAAAAgBYIBQAAAABogVAAAAAAgBYIBQAAAABogVAAAAAAgBYIBQAAAABogVAAAAAAgBYIBQAAAABogVAAAAAAgBYIBQAAAABoWau7AQAAsEQUReLFYrGotyWgYjY3zwy387177+F2RolTfTwe9/t985r9fn8ymdTVfjgKAICaiZZZLBbihij+rLt1AIRlMBgMh0PzOpPJpJK2qIm4KSwAQBdIKgDlXUgohuQL0Hq66SgIer1eHMfmNeEoAADajNI2kOaBgcViAWsBtJjxeExEk8kkiiKDVqgRCAUAQBAMyqDog5HYBHIBtJJ+vy+vCJswRPXA0wMAeCM3oFBoV+aQBGgfm8fOCbfzvfe7k9uZkwooiIxFUoUhEHoAADQV54CCM7AWQFuJ43hjY4N2whB1N2cXCAUAQAE8BhSckZGIag4HQGUMh0ORskBEfFIWIBQAADnUrgyUwFoArUSkLPR6PWKTsoBoHwAgTar3reUuYZ+LgKyF1jC96/xwO99z9g+5nSfmzIPhcDgajYio1+uJ3IW62o8LDABQvzLIVk0o1P3Xm+oFfAGhkCKO48FgIP+sq/0o4QxAF6ldGchmiEPLBmRbYqMYpLZQ7gEAtogQgw4Rhqg9vgahAEBX8Dh20bkBKbfAphn2TUUlR9AsLE/UxWJRb6YCrigAWku9tkHJGszZrQplLRCshQbStdBDU4CjAEB74BBQsIkmuAFrof3M0CVxBL8KAM2m3oCCWzTB10HNyOAu5AIAZYBQAKBhtNs2sAHWAgBVAqEAAHdqVwa12AY2wFoAoAIgFADgSO0jFKhu28AGWAsAVACEAgAsgG1QBlgLLWG6WncLgAIIBQDqgdvYRba2gQ2yzbmKAaWZACgKhAIAFVG7Z0BtkQUG7CvYECIRANgBoQBAQDo4dpEJloWfYS0AkAuEAgA+gW3ABFgLTSRCjgJLIBQAKEXtyqDLtoENsBYAKAmEAgCFwdjFBgFrAYCSQCgAkA9sgxYAawEANyAUAFBQuzKQzYBt4Iui1oL9JgC0GwgFALapPaAA26AaLK0FyzWBR6LN9bqbABRAKIDugpJH3aRQ4WdYCwBAKIAOwSGgAFnAB1gLANgAoQDaTO3KANEEzsBaAMAGCAXQNmpPNSDYBk0D1gIABiAUQOOBbQBKAmuBC6jMyBIIBdA8alcGshmwDVpGIWuB8IuDbgChAJpB7QEF2AZdoJC1QIhEgG4AoQCYgrGLoEYsFQAiEaALQCgALnAIKEAWAAGshXqYokviCH4VUCf1BhQQTQC5wFoAAEIBVApsA9AsYC0AAKEAwlK7MoBtALwAawF0FggF4JnalYFsBmwD4BEHa6HQVqAu4jiO41i8Hg6HdTaFK3DJgAcYjl0ETYftT1moYWw/BU+2bnp0uJ2vPu4bqd8ijuPBYJBabWNjA3IhxUrdDQCNJFpmkaCao6dewDZoOtnflC3yHLNpqohE8P9Q3USohF6vNx6Px+PxxsYGEY1GIwiFFFC7wAoOAQU8nLUPw2/aMpeoHZ8iNFU6CsPhcDQa9Xo9GXeghMeAHysJzl2ghdvYRdB0Cv2mDfrdLZuKrIVcqhQK4ucYj8f9fj+5Gn6mLAg9gF3qDSjINogXiCY0He8RIrYGvuXHEZcS20/RNXq9HhGlVIIkaTMAjHroNLUHFDB2sX2EG2/C/9ywnFMKz6xaKqzMqJQCcqFOQHQTCIXOUfsIBcLYxRZRi9RjG5WwtxYIFndpVp90k98dygQFkdUIJEyvN+ARhrYBaDoV/JSNPlvsG9/oj+mdrS89LtzOV590k+GrFrmNhOGRKuAotBDMuwj8wjZCxLaXLVSdCdZCvSSrKWRzGwFBKLSD2j0DgixoI/x/U27tyWKZtWC5ZuuZbq1XfMR+vz+ZTAhGghEIhabCbewi7nFNp9G/KdteFtYCZ6RKwHduBsMjGwPGLgKPtKy6ZaFqibVgWcYR4ycrYzgcTiaTXq/XxBO+YpjKcEAMAgpIQmwfDfopG9RU73TWWjj5xSeF2/kpP/WlbMGlDn7JDiD0wAuMXQQeaXQ0oSRsdQayFlih82/wzSeBUKgZhrYBrpCmA6lHjD8yshYMzBar1RwIhRcLAaFQNbUrA9kM9CXtAFLPDNvn8kLWAuFn9Uq/38f3aQ+EQhXUHlBAX9I+IPUsYfu1FLIWiLHiAa0HQiEIKHkE/AKpVx62Ha1lwzoYiQBMYDc8Uoz9swkg9ft9buOIlGMXK7uq2zHODSTBb+oRtl9dIWsB4ydB9bCT2PIayK2TFU5cy+P2+/1C5Tyzt/VwYOxi+8BvmiTcl8D26y00QwQxlj7O3PWFXridn/3TN7bvG6sGdhdMUiz3ej2DtRDiUpHzgiQxVP/OSvvQN3q29zhQFMgCM/hacmnfVwShwBN2oQfBeDwmoslkYhmG8EIcx3L2sPF4PB6PxWSjcr6QLLoQg7AHyzuEhvJ5oKEgmsAHth6+ZcN83WcAMMNUKCTHrgwGg2rm6hBHESEPEXQYDoe9Xo8SwQh7ykQT0Ze0Bkg9zrD9LRyyFiAXQDiYCgXBYrEQ/fRoNKpg6k8xO0hKE4jjOrsaNtcw+pL2AanXLNj2svbWQjuSHGdba+H+1f3hGgxroUBEcRwL/1+EIYIea2NjQxzLLynJ37LJeABB6jUftr+XbJjlnFKwFkAI2OXCKFMU4ziWiQIytbCCvF95XEM+YxJDbpF4q33JR10Gv2ZQavl62/GbNvdT3PH5K8Pt/LynHm7o11I73B0FgUhZEGGIalIWhDAvpBJSm6deyCI5kPwNBbZBF+D/m8JaANXTDKEgkGGIalIWJCldEumR6xiiCa2JJrYbRIg6Dtsr1PL0w30GeISdQ5UbUEiGIcxreqHf708mE3NFB0khx6+tJVMaTXM92/aB38KM5ffTrPvM9298WridX9A72JTvgRvsLsUoimx65SrtX/srzeHWhrthjaDkEWeY/ChMmlGSpsgFCAWesBsxYm+sec9UEDuspmaDZIGJXipH3voRTQC58D83bKQMBDEoA86bXXQddlBHwcu2wABsg4bC7ffi1h43mD+WwFHgSZOSGashFfUQBoMYcBEUJCr7BUmIwC/8zx/LARFIcgRFYRd6qJHxeDwYDAaDwcbGhhhVMRwOleUaAyEffKkJdyVWZG0DfIEgEGythUKFn3neZza31utuAlDA9IyvC1+zR5aH7c2IA4gmdAH+Py7bFto3jNtHuC3+mXA7f+Dg06w+bIPgdZYwIY5jEYAQU0PZb+j3quMp+WuE200NeCQr/vBzVwOr+wyEAk9wKfokxK2ts7dLdBhdILfquXkdJrBtYSFrgRjIBQgFniBHgTtso4mBwNjFFuOWSsL/BGDbwkJZC8RA8cwW6JI4glEPDaDFicqYQKHFeC+Dzf8SYNtCy4Zh7BVQAqHQGNp0DWPsYosJ9+PyP0nYtrCQtdDWxxLgDIRCKXSTQgWiodcwbIMWU8uPy/8SYNtCWAvAAQiFUiyWqeygnK9hzLvYBer9cfmfSGxbCGsBOIDMkUbCsDQTkhBbDNt6VrUn3+XCtoWWDas4mXq2WK3gKKAocBQaTI2SH9GELsDfE+LWnixsW+hgLcBd6CwQCo2nymuYf88BnGm0+OPfh7FtoX3WAiIRnQVCoQ0EuoYb3XMAMy1LJZFtZtuTsf1WC311sBa6CYRCe/B1Dbej5wBKWv/j8v84bHtZy68O1kIHQTJjKbhdLQ7l1djmqYHydPnHZZtFyLNVSWy+ukBJjlNUZmQJHIVS1DI8Mhcba6H1T5ZdBj8uNeEjc3vMkMBaACkgFNpJ6hpGtkGLwY9rhm1Pxv9nQtYCEEAotJDkRSuu4S4/WbYV2AaW8P9a2PaysBaAAAGh9mAoecSqNBNwoMvZBr5A1oIzlWUtzLbQJXEEjkJTKWQ4Q/I3FNgGHuH/1bG9QotaC2w/CHAD8k1BHMdxHIvXw+GwzqZkKFMpOVCiMvAFbINqgLXgjKW1YLkmaAoQCkvEcTwYDJJLRqPRxsaGTi6EFs7eew5cw3ww/Lj4dYLC/+tle4XisaSbIPSwhFAJvV5vPB6Px+ONjQ0iGo1GOqEQaHhkaMMZ9mCNIJrAB7aXAP9TwnJABCKe7YCpbq2F4XA4Go16vZ6MO1DCY7D5otyeA7JPlpXB9sGlNdT44zYdfGPUli/B3lr48qFfCdeMJ171/hZ8mbXQhrPQF+JsHo/H/X4/u9yXUODWc8AeDAGTH7fRVPkd8v+92LbQvmE2a0Io8AShh116vR4RpVSCJGkzOMDWcIY9WB6UPGo6/H8vti0smrWAW00TgVDYJY7j7Ekv9YFOQChpXM+Ba9gBtuIPOMP/EmDbwkJZC2w/BVCCUQ8mZIKCyGoUWF4PqRf8YRUQYQjGLnYB/r8p2xYWshYIt5pGAaGgReQ2ElFqeKTh5G6BTMagJgHGLnYc/t0Y2xZaNkx5qzlOqwFbBlyBUFCQrKaQzW1sPbKDpO51ipAFQMD/d2fbwqINY6t4gAQ5Cmn6/b4MNywWi66pBEkXkhwbl0oCqof/JcC2hcqGRQn8lp8B4YCjsES/359MJoQOY4dWRiJgGwB7+J8bbFuY9CZTC0GzgFDYZTgcTiaTVMEl0PTMIyQhAi/wvwSYtDCrDOwbdjJCjgJHWJxYTCj/6MzkQg1Es6yFdv8WXQC/YINIioMyOQqfu+HFPpu1zNMPvBtnlBtwFNLoAn7KM4xtdDAEnK0F2AagMnheAkkqaKFbQIH/VweUQChs4xZuSJ30XdANHLIWMHaxfXArbW5ANo9tU0O0ykuqAeevToyB183/13EgFLbp9/vcTly21HVDhyxoHza/KcNORcCzVUlKfnVlAgpmGH51YkbAulvBFAgF4EhoawHRhPbh9pvy/91bI2WqH6GQOuLJmgoudXYYvCUQCsCdEKWZYBu0Bu8Rotb0x9Vj+OrC2QY21FusZTgcxnEshsQDAxAK3WJz80wi2rPnbo/7LBOJgG3QPsJJPf7nRiOkDAobSKASLEFlxg6xeewcmq7RdG3z2Dmbx84RosEXiyKTwmHexdZQS3VLtlnDPE/jaJnFMnW3rk7EjMGLxWI8HtfdFtbAUSgF2xtWluld50ep+N90dZPOJCKaru057Uj5Q+isBdgG7aPeCBH/86d2a6HegAJoGRAKpWjK8MjpHRdqfuuTRETrs81j5xCRF7lAmYdL2AZNh63Uq70/1lF9q4oGFHh+dYWSGZ934NpwLQFJIBQ6TXTPaYszjsk/N4+dE01X18/6F5ddZboQtrIJWMI/sZRbe7KE649Lphrw/+pyue7QNYXWh7BwBkKh/UxvfwDpdXp05IzFGfctrX/X+eKFQTHYPGJyKM0E7GFrG9jA8/mYfH+HIQIKbL86wAcIhZYzvf0BRBRNTYZedOT0rXPuVWx7x4Xixfq539tes+AjZlNq7XWQllW35FzyT+A8LCj5Z4iPxv+rA7UDodByzBLBcuXpHRfSdHX9gn9yu4/AWuBDO2SBAbYfx75hdeUhcvjqNiMMxOMIhEIpmIfh59+5JKJ1mzUVEmGWPjekObF28X8v2pIQpZlALo2OJpSE7fOxclhQ8s/am832qwN1AaFQiqaMelBgdhoyKkGyWJ/PbnuQeB1trq8+9B/tj4lIRDU02ja48KPXixcXfOTTZfbD9iNnaxFyayq39oDagVBoLfO/f5ghhzHFyu37ty64k0gjETSqYrFnOv/OJeK1vWJAJMIv7bANpD7ILvzez11dZs9MhGk2oMCkYQb4txBUA4RCl5iafu6V287fuvDOzCZpiaDLY5h/5xKx8urD/yG3IU25UXKmHbaB/crOcqGuryXXNuD/e/FvIagGCIUWst2FGGVBepOZqvvPUwnRpiIBYv73DyMimq6tPupb5oPCWrCnHbYBFZQIqQ2ZWwvOAQX+irmyFt5XU5fU7/eZ/wT1AqHQHpRdiFoB5KKyDQoNoKD12fwbl8u/dKIB1oKZRtsGSe7/1wd3Xq4R0erKzGEnDK0FLyMU+P+m/FsIggKh0FTMj5hbNz06PbNDCpXfsHrrBfOLFUWWshJB4SUYZcTWTY+Wx1154ldT78JaELTGNpAk9MES863t089BMUhPwlkxiK+399l3nRZtEtG+SNGGU6IZEX3gwKuy2yb/9F5Pifkvzr+FIAT41a0YDofyv0mywxxCf5+WF2qyY7bf+XRrfeVB/7x0OBuJQCqVsHzcJWMj8VZKNHTqNpSVBa1Bpw90FJILaytz9fJoRkTfffZzxJ+Xfuo6ufLqymwfzU6JZqdFm0lZcIpKIiQ5hWZE9O6rflcuadkvxYHk+f+uw/8x3IFecuVb8fO5AUfBitFo1Ov1sstDD490e8Tc+tLjdL/sdMuqrALZSwTKqISMNEmHP9Zncp2tLz1ONmm2WO2CtdCaaEKWovqgKEqJsJbo7NdXp5d+6jqx2t612Y+t3CNlQUoTCAWgI6kMiOhFB1/v3GY3+GtH/i0EHoFQyKff71d8RC99iaUmmC2WevG9yuwEV4lAFkkS6ytT0dS1aH7yi0/60Rd+ilp0G2pfNCFLSX1gYydkJYLUB+urU7HC6srskWu7gTOzFEhhrwxefOgNRPQXB37HfucO8D9JArXwZOSUUwUCA6GgZTgcxnE8mUxCH8hvX3Lyi08iVTXGlCDQbv4/L9r7Y8vRB6dYA2lHUuSccmvRfLZY/dEXfkpYC0f/9qli+VlPCf5DeKTFtgF5dQ7MKkGnD4Q4ePwp/0QWgmBfNFUuf/uB18jXRT2DFx96Q2itIOCvmPm3EJQHQkFLIJUQejIeS0FARNOF+tffS0T2LsL2QneVkHI+ZPuFtXC/p9xIREf/9ql3faFHRLOdJLhz/80N6ubVQRdsAwofWUiSlAhr0ezJp303KQh0fb+ZpDIgopceem3ir6UdHl/ku3EVWwts+2OerQJ+gVDQEsexfDEYDErurcZHTJ0gUDLbWrO1ELaXe5MISoSjcPZP3yherK3MhFa44/NXbjd4R1hsbq1f3P9/cnfokXbbBpJzPnx45+XuL6vLKLTEYCRcfcY33XRAFqMyMCEawEcuUBNOLbZSBpQHQiEINT5i2sgC+Vyu5OjtF599wW2JPdpKBFKqBE2sQakSdHbI0b996tk/vW0tqBtDtGdl+v0bn7aZ2K1f3dAR20CS0AcKZlur5CQXhESQaiB33EEhTizWkuLAXhko2RdNbbQCVRiJIMb9sZdWbVpXnQdVAqFQGJuhDXweMc2aYGnNVCdtqItQLh0hKxF0+iApeqS18MMbr9A2jGjPylRqhdvin9nZ/+5+xLF+4oqPGnZC4SNEbDHrgxS5cuFX9v8XD20y8n9dOUz++b/eMCSiE0VcNC9UphX4n35spQxwBkKhMIZroPbZI+/4/JW5v6lVEoNdIoLA0kXQBRqU7dH5Ij+88YrznnpYvFCukMtaNCei/+/wzymPdfmB6zolC5IUkgiCl51XeBN7jumf5pPiQCiDFMKoKCkXLO0ESZW+AjHuj3m2CpQBQqHNWCY2bhpyBezFgX59LxJBWiNCIpz31MM/vPGK7B5Mn0XD+o77PV2sfevQ826+4fmpI9LOqK0nHHhf0Z3zR6cP/rfzP+vrEEV7XMFpO5kKxxbrStsgKG5trixlgZrQH7OVMqAoEAqtIlcZ2PSjs8Xa//j+pQ8859bsW4X0AfmWCElun1x13lMPihfqJu0eS7G5GIRp3lAmThLR3sWciL5+8JeU47yPLwdWk+s87cq/MB8lNG8fp4sQp3jNeZrvwXiHPxmt2LdBJifa975/eMXrkn9WoAwkbhJBcILWiOgFh96Urf0cDrb9Mc9WAQcgFNqJ/YO1sh/NwkciyJ0IiXBB76B4UchLsPRaklpBIBRDSi7sozkl5IJYR/C5G16s3PNJVdKWcqFgM69j1k27t1fuXFPKZg8pMgxyc8r2Lrbka3vRYE4PTIqDPzj8avnaEIMw4BB3KC8RaoF/f2wvZQoJUFAZEArt4bb4Z5SllgSWgiDbg4aLMpCrREhy++SqC3oHieh/jJ+RWdl2J4aWlEQpLIhoL80powz27vTZWcWwJ9ExK0XDqTQj4yy9upYoEeohRAp6UiukbIOkOCiDW2qCm0rQ6YMXHHoTqeaUCgqsBRAICIVS1J69qMRGE+Q+UhetqxhUIuh2IlwEIREeOPi0fF1IIugaYz9gRLCP5sc99ax7aW5wF/YstnQGg6VcIK+1cvcutuwfBH9v8Ab52kYZFLUTKlMJNhZCxTEI/v0xWykDzEAolCL0pFCFMOuDXGWQ6iz//uilDz/7luU1ihVN8tIrmyVCEikXbj38rPLtMaDrYn2pBIFZK5g5lWYGrVAxSWVARWyD0BLBOdBQKMpQsVYQsO2P+ReaBEq43E1ASb47frZ4YRN9L9Y1htQH5EMiJLn18LMecsUnxQvnJulaZXgK16kE5wd3Z5VgtXOvU++k7ISUMvjj8fYQgEIdcyGVUE06glsWQvUVHagJ1gL/FoIkEAr59Pt9/qe1QR/kygKtx+6jqKKhAb4kQspKERLhIVd88tbDz/IlESiAkWCQAmaV4JzYWAibBAUpEZLiQCqDJJYdc1AXoTJxQMsNe+7Bt37kqle47ack/B/cUy5smXxSEA4IhbbhLgt2kD3iV+59xONP//bunn0MZDA0wM1C0AVcpLVwyw0/W7JhDkaCeSuzDggx9sGmVbtHsZAIr7rijck/leJAYnP35yYRyEkl6FpVl1ZgrhKIaLFY8Mz0AkkgFFpC0af2JIbOw1eKorklJS0E3d6ERLj0yo9JrVD0W3KTCIYNnSWCDeVVgkEiFFIGEu8SoVn6IEmNvgI1wVoAnIFQaA+WWfq5vYWhC/QYXzDvzUEi6PZ2yw0/e+mVHyOibx16nn3bzN+SQ0ZCrgjIXcFgJwSSCElx8KbDv205tMG+e+YjEarJP6hRK0AlgDJAKJSCiWn29YO/pPwpy2iCFKFTFMkYYnDQB8m2CYlw+YHrpFZwsBAE1SctltEHVFAipGyDNx3+7Z2drJzMPVIYiRB6LENRieA6/JLLnZa5teBWXAuEhsvp21BYDY/0IgtSO/GlDwy7EhS1EIqO9iSibx163uUHriOir33uF7Pre7cQqFygwZyO4FEipGyDzE5sCyTUGGWoIEvRefBCSiI843Pv+PTTftNtV17grBIAWyAU2kC2S7C0CnL7ksPHL7ti383JJdVYCOQqEcxFGoREeMzT/m9KyAVuUYaSEiH3N33ZFX+S/DMrDnb2ky8RvFsIoQc61qUPktSuFQTMrQXACgiF9uBgGCy9pd/c4xAGChliUG+YaaGUC186+MvKTXK/yXC5iiESESyVwc5OuOuDoMpgeysXvVJgEw5aASoB2AOh0Hgmh36VMj1Q7mOl/VNvths2Z026jXIkjUrwqA8E8pv50sFffuJV7xcvxBK3gQzb75ZwEUJEGZLi4E8Pv7JQUQQDfuMLQfMPKjAP+GQelAHWAsilDSc68CgLDFSWf5C7T4cMCeVXJCTCE696f3zwhYaWBIoy+C2KkLIN/vTwK7ePQqubOTvxZiGEkAgMkxPd9MGxxZ7kn73Pvmvy9F9z2I93WJVVrqWQJcgFv0rjUXZjHsPkDsEFChBfoBK1IszlEOKDL+xf9V4iSskFtomKomE6ZUCexIGAf4oi82ELKX2QhI9WENSuEgBbIBRKwWR4pI0f4Pxo+/Hp5c9Z/1Zyid8SioYdCnz5B5JsiEFIhP5V740PvtDZnikz0JHs9EEqoJDev7/pIWrRB9wiCwIv/kGz4GAtAFZAKJSi9uGRf3PDr2cTFCiM413l+AW3Co/kOoRBbPjZQy/qH3g3EX320IvS75bLUnT+zl9y5VuTf2bFARWclMGMxxBDoPzE0PqgYlnAzVQQQCWAFBAKbaCMLCCLZ1m/JZYNOwyhDwqlKAqJ8PQD7/7soRfVVUsxKQ6UymB3//mjMyrNPwiUnBiopvJyS/xHFgwkm/fEz7znS//2V92OHhpYC0AAodBsUn2VF1lg6HTdRi6Qa3CBSuQflBnC8MlDL3n6gT8XLxTb+g4x2NgGS/uvagiD98GNfIIL1QcUDM1jqxWqVwmYPZInEArNJuh0grJDYqUPygQXDBvSsgIQEuFZO3LBb3AhpQz+/IZX2DSPKrQQ/PoHgUYuNMg2yIWtVhDAWug4EAoK4jiO45iI+v1+v9+vuTUaoij6wA3pmi0ly/bp+qGPLB7x3OjbySUhggtUQh9QiVkYDCLgukPXENGzDlwrXmSxt3CS4kAoA5sWelEGArb+QcvSDsoM8GPbH/NsFagMCIU0yYTE0WhEzC4SeStZLBYfPPxSMoqDMtMKK3sg55ELVMngBfttCyUnXnfomucduJZ2dIMgVyK84Mp3JP9MigM3uZZuZF4DaimO1Ny0xApsgyzHE1/Coz7zfla3GiVspQwICoTCEsI/6PV6w+GQiIbD4WQyMVwboYc5yEMn9YF8NyURysiCnT1o+56iZRMl3vUBlYgvOA9eEBLheQeuzbo4EoMyyG0Y+dMHxNg8oCL6gK1z4FEc6GDbH4dulVvJbRAa/Cq7xHE8mUzEC7lESIE4jpUxiMg8JhgAACAASURBVEDDI7OyQHl9hlMG5vu+c+YB1aQPyCgC7IsjfeCG3xRqQMqFpDhIKQPLhvkKLnj0D+o1D4LmJFYfULCRBZJHfeb93/i3v0zMXEwlbKUMCAGEwi7CRdjY2Egu3NjYGI1Gw+FQqgfvmG0DHW8fvyq1JJAsEAQSB+RvZKN6HVd9QJrggpAIUh8oxYFN21ilJQrqrYkU1D+oWB8UEgcppFYQsO2PebYKBAJCYRdhJ6Scg+FwOBqNxFteMMiCQtdemZ7GwZp+Jz361+nrqYW16APztpVN8VyXf1BL2YMQaYkMnYMyMYUy4kAwV10vrCZiUMK2YcAjEAppAg1zKCkLbHBWBuZOZffuGRGFqXkgCBRZyF0h6BDT7UNUpQ/IQiLU6BxwS0isINtAh1IWCC779F/d/IxfSC1k2xmzbRjwCIRCENyiCYUo2bXouhPDrTPQnEzO24ablklg1gdMggv8hy20xjko6RkYlEEWpVYQsH2C99IwzB7JE/wqhbHJWAxnG+hw9gxyr8zdPS9/7lx9QE5lE8sMW8hdp+S0TORp8EJ5/8BjWmLtOYkViINqUhGVFNIHSXRagadKIMYNA+WBUCiM4XqobFKojfFriMIrAxW1lD2gAGmJSbzUr8ylAn1Q44DGcM5BW20DHbMtqzOq3dYCYAWEQhsok7lm6DBSPYROIgSKLFB9czbmNoxanZYYKOcgqHNQl2cQVBlc+qnrbnnm83Rbse2M2TYMOAOhkEZXMoEVzk+W9rIgu7c/pCf9Af0dlSt4INDtoWRkQeA2Z6O5Ybs7D68PahnQyEEcVOYcOIuDKg0DgVkrCNg+wbNtGCgKhMIuomRCSiiI4gq9Xq+uVtlguFea+4lccZCilmELodMSqXRkocpqyszFAc9xCvXaBoWUQZZcrcC2M3ZomPPEXSAoUHxLiCSD5HcilozHYxuboRoF/bIbXmd411kZUI7gWCOi4coXlxYGm20hd4WSaYklp2p81RVvfNPh365s2EItOYlMnIPqlUG9skDJdL7+3Wc/J3c1tk/w5oYl333G596hW6081z/9pTy/H/5Avi3R6/Umk0m/35dzPYjlzIMRZQIKmh3mJTwGmKox992SOYl+wwplylxStfqg9joHFdgGNUYTAimD1JIHf+LjuVqBbS/ItmHAEgiFJUTcYTKZDAYDuZDDMIckfqMJy3tWryBv9G7BhaA5iSXFAXEqeOC31IGNPmh0NqKbOGBlGGQFgRcaai0AtkAopIl3IKJ+v2/2EgJNCmVPUFlg3mG4ggflRzNScwoieRQH9dY5YBtQ4KMMSsoCG1NBwLYzzq1I7TwrBwgKhIKCXH1QLyXzDLIYLs7U3n5n1v9P6583N0/X/ddb7YB8KwNdz81zNCOHhINGBBSYGwb2WkHA9gmeZ6uADgiFxqMfCan9ce2VwdIOVWdLyZxEgU4ihM5J3NlJw+okes85YJWK6CAOyngGfAwDA7PE13jxxz9123Oeabkh//64luhtCnsLubNAKDSM5x58q/lXKxlQ2N2POfbPexImMuoDe+fg1U/4M7nkdV/+jeQKFTsH3ocqhEs4KOZe1JSEWEYfVCMI/MLZWqhXKySPPhqNqAnqqnogFJpNNbYBJe7L4QoecJuE6TU/+U6abu9zsT5/9RP+jNZnRPT7f/sqwx7Y1jkIWvsotDioN88gkDJwlgWFTAUBOj8lwj/o9XpymNtkMmErqmoE34hPKjjDDOOMAymDFLOt1TfuOyheN32GRnOH/ZqffKdy+eLUk3v33nPyR/t//0uvSC4vOYnG7n785RwEsg2YBxQYyoIQVkFRrSBg2Asmm/TEz7wn3IG+fPV/SH72OI7F6DbnwjndAY5CKWoPsBUNKEiyN27DfTl15216EWWy6LNf88h300zd1OieNaJ7iOiPnvjWaLr6yq/+xzIH8j5UgUPCQYNkAXlVBuFiB1kcfAWCtZBAuAgbGxvJhaI+73A4FFkLQMBOXTaa0Gq999l3Gd71axsol4tb6hvu99nkwtBFlLePUq6UcqE6SH/8cNP3TER7LvyfJ3+0X7yOdsITQjF4rJBY7wjGQg2gCsUBK8OgSmWgxM1XIE7WQl2Ogs48yNbnBXAUGoxDHqK9bZBEd29twfSM9jMw6Visz4VW+JPHvUWYENd87beyq9U4VCHcCMai4qC5yqB6QWD5qS/86PXf+7mrHfaPjlCAEIMNEApNwiGgIMjeoB1kAe3cLp0jC42ok/h/PPQvZQ6jgZUjp4sXizPu2126NqPZ2rWPeTMR0XT1mptfXr1tQM0PKNQ+JKEaWeBrZOaFH72eiNzkAnGyFgBbIBQajN+AgnoTuzumc2ShmrkZS9Y5iDYVm8uF0R1nEtHWuXcTUSqt4drL3kZEv3rzKwsdLoV3fRAuoABZoNh5gJkglJS3FjgohkLn2zev/pVwLQFJIBSaRFFlQE7OgfK+mdxP7QkHghDlk9/ywPeQShbkMF1duX07ZWHrvHtSb77nsj8RsYmfv8U0rlIQyDagkM5BUX3g3HdylgWVCYJw1K4SivLI699XaH0IC2cgFEpR5aiHR33m/aklhZyDoraB7sY331p79Z3PfN3+T2XfCl1B2UvOgW1/mRt90Kyw8sMzxIuts46l3vrQpW8iomhz/fn/7eWpt2z0AZOAQpXDE0oqgxCagLMgKBmDkHCwFgArIBRKUeWkUEUDCoJC+sCwq+zRm1I+uVh5xHx9sJZ+oWHlh2eKF1vn3Z1cvtgz/fCPv41oW2r8u396RWbTJToVUOCjDDgLAjPl5UKnVIKYMbjuVnAHQqGR+AooGHZl1iX15hyQJ32QOtA7z//Pud2/jmim/eCrt+8nmcSQZH1ORH/zgLcSEU3X/t33f1MsbqJzUEs0wa9h0FxloMQ5ZSFJ9dZCyRKchRAlE1JCQRRX6PV6lTWjEcBi8knQ6+rST12nXG6+1RYKK5D1+Ig/PG839MDBNiAfEyu88/z/rFyeFQGrj/vG/KuPsjpYRnkszllOYkh6GDsrX3G7YoAlsZEFVIdh4EUW1KsG3HrB1ZVZmYOWlwuhSd42L/v0X4U70Ld/5heVHjAqM+YCR6GROIxTUN4ii2U5bO0elEk2InkqnHyC1t63/72kdwUUOHsPR86QrxdnpPMYBIcveLNUJ4Mfvtx2fIT15cw//bBBsiD0Q3By/w6iQUQiqJxi0D0CPeNz79gXzT5yVU74jC29Xm8ymfT7fTnXg1gOlZACjoJPwjkKD/7Ex5XLiwYUBJa2QVIZLK+8RkSvu+CjynfLlEesTBlQpgLS+/a/N72GXgqsPumm+d89VvmW7kszsHbWXcrDLTkZ07UnH82kQAYubFC0o22xLKjSEs/F2WNYW5lTwWKOIof67OjEadEmEe2Ldg99SqRtxgcO5A/wkdToKBBRv9+fTCbJJegTszA6+4E9XmwD3fpmfZCCScKBh/mW7BwC+eUUEgSzhakPmx09R7w45X67eQzZeMcXz35b8rj/+u6X5R43tHPgpgzYygJWasDAfGutqFYQEkFw8ce344ZriZ5+fXUqVltdme2j2SnR7NyV+4joklUijSY4hbRtePGhN4gXf3Hgdwq1s3riHYio3+/DS1DSjAujdlLGlKSy4ZF+UxFT2CuD5OY6iVBXKqJiJ66nd1FXwKwDTAda/pjTe8+Rr8/Yd5dyk/WVKRFNt9b/3zP/VB76yff8BnGVBeRDGfiVBU0RBOVJ6gPSKAMiOnf1R6dFm/uiWUoQGKSAJUIx2MuFWpJIoA9y6coFU5LRaKTMgw09PFIXyyjvHBj6QrM+0ME5rJB7oA/e74O54iClBnLFwdS1dxTf/53HzpVL9u9RiwYiWovmRPTFM/4sedyn3P0bqt1WVMagpCzw2FW0VRDY2Ak6iSD0weV7bxeaoLwUsKGoXADcaOeF5JfqxabUB1IllE9FJI04UAYUDDtZmW+LoWY5B5b5gIW8AQcpoPu2zfzLyW2n4fy9RyjP8FiPZl84889Sh3vK3ddYHstBHJRRBr5kQVs1QRaDSkiJA9rRB08+7btZQbAvmnppj/10qS8+9AZohYbSlavLATEleSrPJQRSFmT1QRaHas32+iC7BykLUkTzlTf8t5/9nR//WPYtFgkHBY61Rnn6QCkIdCqhqBQoJE2kYtjc+U3vv/eIYf21lZlo0hfOvDZ5uN49LxF/VikLGq0JSjY+24W7oVQJyZ1fdb9vixe+dEAuyQPlXrmwFhoKhIKWQCrBIAt0+uD+f30w+0uZai4VyUZU7ierD6K5NmmxykGMJY9CxmLJlg6B+BotBUHR9IXNgukRP9iRDufowxNJRKhicsafpw7Xv/sluk1qNAxCa4LKIuLJAzmIhqw+uPqMb6akgGEMgpKiU4rbsC+a2lzIsBYaB4SCFpEHK14MBoOSe7OXBbmEiCwYnAPFwtn2ygaJ4Ms2IB/6wDyTwgdP+5AXe6CQJrAXBJb99A9O7uY0zBar5xudBsmelaloTHzmn2ePaFAPSsr3uyGUQaNLLj7v7C8XVQCWpHYbQjfoMKgEL/N+Ae9AKAShUDQhl0A5BzpxQBl9IJXB0uaZhV4SDii8cyBJFksu6hDkyoKi3oD9g7uNIpFBiv3rmdLRGfasTJOtXYtmojFSPSSPeOCeFxFLq6ApgmC2tZo0FX5l/38hotPKRQosJ1I3IHSDX7kA26A1QCgUxmZoQ3nbQEchcUAlnIOsOHBQBl5kQZHDWe0qNZmCVlrZOQTevQGHwZaGoMkPNs9J/jnbWrv/3jtydyi1ws6fc9mwQ2e8W3lE8TVefa9iJl+/sqApgkDwsvMOW6YL2HtsWZI6o7xoKMSfXvnq1BJMPtk+IBQKY7gGQpRVOOOvPk/kno0oUUoEN31ARJHejRCUH6oQwjZQYt8xJzWBQR94FwT2wytsfBEZpJAze/2rtTszbcvZz3o0SzVMJE5ef/r7Umum5g87vlP8+/l3vcB8CJ6C4DXnXV9o/dyKZBKlnnBQD6dFU2etcEo00113WUGgAyqhfXROKMgiXAayhZU4YG8bkGvOwdIelhdmlUFyhcpGMHpRBskjfuLUD2aXW5oEfoMF5FsQkN1U4ET0/dn+46rJOx6wetS84Xo0S7V5bWWWatvexTzZkn207bp/+KwPWDb4pKptyoVEtKnvmO+zvt1Zfm/bRzTOe7J3sZXZua10oB31UFQuOGuFN16x4bCVDgdrwe+MoMAXnftVeIoAA4WcA0vbYPetRGefdQ7M4oCq0gduAQX7A7nJAvW0nHUIAsteTSkFzHv7zta5ybceFOXohkIIxZBqVUpVbC/c0RZJcSAXpt7ak+mYaUc9nKopLpQVEKIZuzs3fsN7Ei2xmSxNSIdCcsEXfnVALrAWWkPnhEKuncAZ+8iCs3NAFvpA7vwjX7v6uY9ZcmIrcw4K2QaGd5UqwVfswEYQ+PUGHNSAabXlvd2yOFe3wsMpZ5DF3sXc/jFdt7IQB0ovwfAWEe1ZbBmcBiEgDH6DUr4wZDR4bd1NUIOshabTOaHQLPa8/8ukuvd5jyyQJu1ArqbIZkisH3oEo8RGH1gJEVojoo/vvc5ZE6SWFKrBYMaXJnBTA27rfD06b/cP1boP3brHpjHhMGsFb0ehuY2pUIj/vf96yzXZ9sc8WwXsgVAoRWWTQglsbAPSKAOyS0vUmQe69ZUSwXKQla+EA9vD5Z3tHus32xdgyO3ObUwCG02Q29nbKIbcvlb3XP6Nlf2GRqY618u20smVHChkiuh45eD/9NIYJfz7Y7ZSBpiBUChF6EmhyoxmpCIDFsz6QBGMmLoVUfCWiuhFHHic2ClXGXgRBJa7Kq8JbB6+dZrApjfVPXaLyP1Nq+kYhy4bJpuyZzgxjifCTFmBm90webomT8vkmskTLDUKdLa1KsoHfffZz9E1KRBs+2PZqkLT3YHagVBoAOWdAxtxkDqQepOEPhArmDvsWgIKlmQH9xdKKdApA/uHzvKxgwoEgc1IAV07DSa8IZXPkCFrKQuye7DxvVLnqvKcVEqE+XZJ71VKFBY8Md/7vZ+7WvEZKoGnSkjCv4UgCYRCPv1+v67T2j7ngJzSEpVHSW6Vcg5239p5EY+f1h98LrkOt4CCICsLPrz3o7nKoLxPQJ6SCSoIHORqAvvYwfJW6oMWEgQS1XO/YuVcw0B3lhZ1DqQ+EMkus63VGvVBFrbWgqTi6C1wA0KhAZTJOVCumZt5QCrzYJud14tZdOSPHkZEz/jcYU3D0ziPYEyRKw4cSgEWTTPU9ZoV+AQ265SxCswtNCfr+bIKBB4NA1KJAzdZIBDmgdAHtz3nmcr21w5zlUBEi8UCWoE/EAp8OfWdN0WUbwlYhhV2VjYlH2j9A1rSB0REczry+odpm05EFcqCJLkSQd7rzcrAV0pBaLegfOzAzSfY2dabW6A7GdwMg51tTdEE9SZ6iZCcr2i2WFOaB2yf4Nk2DDQCCAXu6Lv8wmkHpg3z/AOykAi+ahtQiYBCFmV61If3fiy1pLwy4GAVmGVBbgtrjCAY9uZsGJDqnHR2DsTC+dbaD/79VcrGE+MneLYNSxF6YnHgBn6VUgQ1zZwjC6RLbtAnH6T3rNIHRGqJ4MU5sFEG9jeR3NzpagYlts8qoJqUgY1hQEVkgcSgD2TaAeXpgyxsn+DZNgxwBkKhFKGHR5J++iUv5oF6nSIqIYvHgIJHZSAwzHZf1DDIdsOGXj/ocERDk7Z3zkAWGPbjJcPAtKG1bbC9PKEPnDMT2XbGbBsGOAOhwJTT3/a1VIKCzYCFnYUF9UFWGWwv3/6/L4kQOqCQRacMDLLA0idoqCwwTzFgnryjZMrh7vp2hsHOfrzlGQhSZ450DuS7Rc0DM2yf4Nk2DDAEQoE7DmGF7Q0zwQXFaq4SQd5loih6wvV/qVyH+NkGSkpGEELLgkBjEMjVLaAwhkFiV+oyR7rNiwYUdpeH1AcStp0xz4ah4BJPIBSYEjCyQDn6YCH2MI/ufMtPpI+4ow/kXUa8eOJn3iP+rF4ZFJIFom941+kfc+775VsO0xwLGOYWCLwYBrpZP/Q5CsXyDHY3NJZETK2cHMpIibQD8frI869Qts0vbJ/g2TYM8AFCgTXOkYX0aplNlP5BSiIkbYOUPkjhRR/YP0xYigPLCZ+CBhEYygJfKYfbmxQxDBwGJuxumzmQg21AlesDCdvOOLesMgAQChw5643fjChfIuQHF1SSQikRSKiEeURER996Kb1FYRsUwqNt4FcWUInBCElZoJMI4VIOqaoxiuYdGmYKtal+aBNNoIIBBcqcJMmERLntbGt1urV+zy/8G137q4Ftf8yzVYADEAqlCDo8UjveociwBUopA8qIg8Qm27GGt7jcL6oPKNgrA91BYRjk7s1eFpCxerebbWBIRcyeIanIguDkbC8RTbfWiah2iSDg3x+zlTKgLiAUShFoeGSuc5CbeSAwmQeZrbIZCfYYVEItAQXDQd965ie9pBcoJUL1iYc1GgaJvYUKKJA+pkCqsEJyD9uzMGyt3/uLT+XZ7bHtj2ts1XSrWCYyqAYIBb5Y6QOzOKB8fUDlJILg5mf8wmWf/iuqL8+g0KEldVU5pDDjEUIbBjt7sw0omEsaCJQSU5eNuP2uXh+IP9PJB7+w4NklM2xSCp7fG6geCAV2nPPam1cSFRR0AQjtyIXtd3df6vQB+ZAIEi+pBuWjCVlSXZFbboHNlA1usyoL/BoGhh2WH5Ww0zArZbC9E2NhA4Eh58DsHCSX6/IT5cxDPLs9tv0xz1aB6oFQ4ItNfMHgH2yvkFUJs4iI7nz7JR6aaMSvZ2BvGNhkS5TJLRCEkAWUV8CAitsGRZUBFREHltEEshjBSPrIQmrl9FxNdvmJqfE7rGDYpBQ8vzdQGRAK7DDURCKleUBpfbCYavYwi4hoMY+OXvvQsq3UUE2eQZJCs8i8bv+nHHILBFVOpShxiCZQyAxE9R4slAHlBRS211lo9YSMLGxv6zR+AdaCG5W1yma2d1A9+FUUxHEcx7F4PRwODWuGGvWQm3mwvVpiBdVQyaSFQPVJBO8BhULiQD3eoVyGAQUYkkCu0QQqaBuUVwbb+7FIOMh1DnRhheTeUvqAyo1fkNYC8ZML3NqTha2UAUGBUFgijuPBYJBcMhqNNjY2dHIhyKgHc/IBGf2DpMJIbLsQBRKCqYQsHgMK5WWBhE8NA6pKFlDBPAPdfpxtA8pzDsxpB8k0+NnW2n2/9GRlm4vCORJBjBvGs1UgNBAKSwiV0Ov1hDKI43g0Go1GI8qzFvySm3lAhvgC1SMRPEYTyFocFBhkkehvysgC2unCHdIPyWnIokOSAVUYUBCY9UH23DCLA0r8XsKL9qUPUrCNRHBrTxa2UgaEAEJhFyEFer2ejDv0+/1+vz8YDEajUTihkKyUvP/3/p6Sk0bmSgSjPiAimtHRd1VhJNz2nGde/PFPpRb6VQZFB0Bmx2TrJEKgyZOoqiGL5Gl4AnkKKJAx2yC7W6V5QMH0QRJYC27wbBUIBITCLkrnoN/v+z2KYQKFxWJxzu//A5F18oFAqQ+IaEaLeXTXXz7EX9vz8Z6B6OYZmPEiC5QSwW+NIwdZQKXzDMhpBOP2CkZxoJ2aoQ59kIK/tcBWMfhtGGaP5AmEwi69Xm8ymeiUQRzHZUSDbt7FNKn5mZbeWl6SSV/YVgmz7dcVSwQyXuTeQwkCS3Egup8/uPDjJ1UzaDhPq0iuoxIofJ6BYVdkUSlZ4DBOIbuhpT4Qb23+8hN0bQ4Nc2uBZ6uIccOARyAUdpERB+XCQirBft7FFO76gLYlAtWkEpT4DSjYewaGQVa+ZEHQyZcTO6wuoCAxVE0mC2WQ3b/ZP5htrdaoD1KwtRYkbKUM24aB8kAomJCDIDY2NuRCm6ENbvMu7n/VLUt/28UXssvrUgne8wy8KANJ0fRDi168mCzI3a2X4Qmkn/Xbsl4y5SUcKLe1VwZyHT76IAnn8ZPEskkCtg0D5YFQ0DIcDkXWQmp4pOF68D48kqgZEkHwg39/1f3/+mB2uaU48CsLUg34vYs+KV77Sj9M/VlBBmL+3ryOYNxeWV/nILV/RZxia50yP9Zsa3VlHp144b9WHo4PzCMRxLhtbBsGnOmEUEgWUNKRlALJagrj8dh7PqMWjYVAGomwtHAa3fX+BwdsWxE4KIPUEu95BhVkIJp3WCgPkTwlHGQPYWkeENHKPCKiTfYSIQnnSATDJgnKNGxFN7UNqJVOCIVCIxv7/f5kMqGMkRCa/S/9x6WBkURkYSHQNCKipkiEQrmHvg6aorLJlxO7DZVnIHAOKOyuXGQoY3aFXH2wQsTfQtABa8EZtg0DRemEUMi1EyRSJdR7fmfzExVvMZMIWQIZBvayIPl0Uij3kPLyDHQ7JKdogmFvVCIPkayVQXbD7FF0kYWdd9umD1LAWnCAbcNAUTohFCwZDoeTySRZcKkyzr7mO1k7IaUSUhJhdXriyIceUUHbCmEjDgLJArLwLb2PWqRKAgpUbgQjWQ9VyB7IQR+szleOv+hxytY2GlgLzrBtGLAEQmEXkbpYvUrYJeMfUCYRYXV6QrxkqBLM2OgDv7JA8oKLxkqJ4BBQkLIgqw88ZiBSkYAClRMHuQcyKwOB0AdE1EqJkATWggP2DYvmOVO0g1qAUEijG7mgPNe9zR6pzD+Q7FgI4q+12fEffIT77di7LKCCiU7ZO05uNIECz6u0tMMSAQUqbRuQfjTj7iGWZ1tQrrkyj8T3HM2i+379scoDtQ/m4yeJ8RM824YBMxAK27gZCV5mj7z46kN00YMoYx7Ilw2SCB71QUlZkMRXGUSBEAflKyVTQduArPWBg21AmbCRwTwQL6JuWAg6OEciGDZJwL8iNVACobBNv9+v5cR9yJWfuOfiy3cnd9DoA/GCuUogovt+6cmnfvCLySVuuYdmnP1JB2UQtKRBCudUxOy2yoOaEw521snRB0QUzVeiWURE3XERdHCORBDj/phnq4AOCIU6eciVn9j9o+ESQeLdMHCQBdFyFQrFHMoWnkGW/BLO5QIKAmfbQHk4G3FA1vqAOhZlsAHWQhm8RW9BSCAUamNXJagkAjVTJRgoWkrFUh9EmckvdFQwcFEQKKBg2E/uUIXtA2nSErNbKfXByiz60W/+pHLPANaCG/J7E9hfzqBKIBRq4NKnfoiI1ohm66eIJUp9QERr0xP//MmnVNw8LwSSBVT8VtJ/0N95GaGQkgXOEzBK/DoHyjXtIwukEgcEfVAEWAugrUAoVI1QCYK16Yk7LnikUAlJfUBEe+87etv1B6punA9sJEI4WbCiWt/ZMxAUSjgoP4JRtx/bWRuMkYXsJgZ9EM2je1/+GF0LgRJYC21FFOqtslwvHyAUSlE0wJZUCYKUfyBenHLsjltveHbJttXFyV99wt73fDm1sGJlkCRopWTypAxsUhHN6xvGLCjXT+k5mZ9IRJAIJeE8fpJbexrEaDTq9Xp1t6IeIBRKYT888rInv1cGGiR3nfdQ8UJKhNXp8fXNY41TCckPvlgs9r37q7Yb+pYF6f3vdIdFBy5S3iAFe31QyDZQHtchITG7odI5oOWfABLBI5wjEcS4YTypbmpAlkAoVMFlT34vJdSAYLZ+SlIfiBdNUQkpSWR5x6lMFkgeccnXLIcnCLzYBuTJOahSHxARJEII2EYiuLVHUPSSD81wOIzjWEwA1GUgFIIjVEIWoRKSEoGI2KqE0MqgvCYwU1QZCAqlIpZMONAdTjd9hiHtgKAPOAFroblAJQggFAKikwiCe8+6SKqEfceO3HLjz1fSqAKkAgoFNszr9UPLArl/L4UNKEDtI/PK5tEKum2zaQe7rzPxBfEaEqEy+FsLUAxZZMXeOI4Hg0GtbakTCIWASB0gma/vS/4pXIS16QkmKsHNNrAhqDLI3bmvgAIllIGvwsm72TaXegAAHdFJREFUB7UIK+SaB5TnH9DOdwuJUD3MrQWerQIcgFAohSF78VGPvza7UEiH+fq+E6eexUElBFIG9/36Y+/3jv9aoBleZUGS8nMwUkHboPBBjaMVlNsqB6Ba6oOVWXT3Kx+tazmoALbWgqRGKVPoVvCjV6C8R0VAKJRCOepBKRGSrE6Pr01PIaJaJIJzQMFbA4LJgsQhdntNXwEFw65KOgc24oCMzsH2EpU+IKJoHok/IRGYwHn8JLFskpL7vbXA0whBWJQAQsE/3/jKNY99zJuTS2ZrSwMjT556NlWoEsIFFKyOXq0sSHL/h/2DOgMgmHNgEiUWzoFyD5aRBcqIA9pNV1y567cfqWsYqAvmkQji3baSxHGcO2NwN2srKYFQCMJNX/utpFZYm52gHbkwX98XWiLUqwy222ChD4onLrhMGmmTcLC00M42mBsHHeQOVTAcNBtZKGQepJbDRWAO50gEwyb5AiKgEBAKbYCDMkhx78sfc/rbvpZcUo0soERnmZIIRQsfCXS2gcfRCmQnDqiIPoCL0CBgLVRPrp0AkkAohCJlKkhu/uILvey/9lSDXAopA5e5pO2CGmWmYcxdv5BzYJOTSMbgAun1AW2Pe4Q+aCqwFohfwSUggFAIyE1f+y0iknLhG1+5pszeGNoGbpS3Cmw4/dL/bjkTo6BQwgE5jVaQ5PoHuckHqeVic0iEpgNrAfAEQqEUNpNCCblQfufNvT6LKoOi5ZUEK1P1VoWcg0JFDiivfPJu24pHFijzdGWQCNE8Ovq7lytbCBoHrAXADQiFUthPCmUJ/4CCPTb6wE0TkF4WSFIdtrlUc3Z9yks4WH7t6ByQpoplqvhBaqHcFfRBW4G1AFgBoVAzrbENbPBrFeTiUK3ZctYlwx7IouAB6dMSE1vp8hZXtoc+TqMjr7lM1wbQAmAt8KHf73ftIyeBUKia7iiDcG6BzYHmFsWPywxilLiFFcguspDam1gHLkJ3YF6aiXxbC84JTCAoEApV0KaAgj1Hf/fys1//Ld27zj5BIf1R3jnQ7cSXPqAiEoGgEjoJ50gEwyYB70AoBKE7toE9bsrAzZZYPPz7yd49kHMgKJmWmF0nPQJiHomvDoGGjsM5EkFcdQzwAoSCZzhfydVTSBw4hir0WxWyDciuDhLtDjQwDWjc3lxf8yC7jk4fECQC2AHWAqgFCAU/SH1QfuBD6/EuCLL4LXIgsBmwkJuWmF1HGV8QEiGaR3dsPELZGNBlWmwtOCc2gaBAKJQldcWKa5in5K+FCmRBksUsooK2AenFAemdg0LDFkiffJBaE1EGYAOsBVAlEAqlyF6ocglnyV8Z+Sqh3APEQtVDH7/8aPLPQs6BLuk6N7JAec5Bdv8p/0AuhIsALGmxtQBY0WmhMBwOlVODDIfDfr9vs4fsZZC0FqjzcuGOjUecO/r20iK35MTiFeBL6oNcZbAz95Jl9EGhDygjEWge3fGHD1c2GwAlnO8z3NoDnOm0UIjjeDKZKJdbCgUzSXehu9dMEWXgIAiWj0XkNPfS7kIL52Bnc9vMg+z6WX1ARJAIwBnm9xm2DQOWdFooCJUwHo9Ty/2qBGLvENaFuyyY569iWedg+y3jaMbEHmyKIqjNA0rog6XlcBGAJ9jeZ+zb41xeBQSl00JB4EUWZEleG8kxEdyu4Wpw1AQWgiDFsUfcm5QI9raBwN48MFsIJv9AAokAfANrAYSg60Kh1+uFPkTHrQVbiVBcEyiONS0sEcx1Eg1rqvevyU9cWmkeie/kyB89THkgAErC9j4j2wPF0Cy6KxREGqO0E3zlJWRRJjx25zo58kcPO+f3/yG91FUWLPKcSZs6SGSnD3LTErNbZY3T7RV2VoNEANXA3Frg2Sqgo+tCIY7jZImkXq+nHAfhC3ndspX8QQgmC5TYRxbMpZSX17RNTlTseR4RJAKoA/73GVSoawRdFwoin3FjY4OIRqPRZDIJqsGTzhs14TKuADc1sEQyKSGvDhLZBRe0BRVskg8SrZKRF0gEUAucx09SppotKjPypOtzem5sbCwWi+FwOBwOF4uFSFkIFIOQCC0irYXuFH5eTKPsvwLbi5GE2X/LrMwi8U8uieZR8l9qtdT60XxFqoTUhmJboQ+SQxh2726JVi1m27kINIdKADXTqfsM8E7XHYXswiiKlMUVfKEr5shW8pensGfg+lSxeemx3HmYzAMaTRsuf4pscEGwm7w5p8U0uvNNl1o0HIAqaPd9BoSj646CjnCZCrqSzy2W/NrO0s4hMDGLlv4RUcIGEH8qbQPacQ5SZQ+W7AGilWkk/yVXS38EIiLatRCIaE5HXv8wqATAjRbfZ0A4IBTc8TJQIltuQcwpVXK37CgjCFJqIKMMJOluflkZ7Kyzogsu7G64LA6UkQv5EYQ+kFEGmtPiRHTk9Yg1AL7IuevqbkgG3cODl3/Ale6GHqIoMoxxsFEAg8HAb3tSoqFz9mDJ+s3WAxbIkM9ocA6WWaoPMSfaibDARQCNgPn4ScCKTjsKk8kkJRSGwyHZVWHyLsazlytTyV8eO3vAzGIepf7NH3qffFfaBsrIQrIsUuofZQMQy88iqfjCtoWwk5UJlQCaBV9rAXCiu47CeDweDAaDwWBjY0P4B8PhUKQx5iYoCD0RjqTMb4nkd3ILFkXcQnvnwJSWKFHmJ26/tf3f7STNeXTnW37Cvp0AsIL5+EnAge46Cv1+X5ZPEIpBN0dUijiOR6NR0NrP2YTHpkv+O99+ieHdrD0g/pU5YirtgDKZB8p1UuFMpX+w/RZUAmgRSHIEBlrxtFqOOI5lOWeb1IQqayWlrIXQhwvK/pf+Y8m+f5dZZsklJ5Qr5tRMFBjyDzI1JYU+OPonDzv7FbfQLDILIAAaR/X3meRdTlHu3R93/vHDm3v/rJfuhh4klvpArkwWroMvkqd105OPiqmErBQoQn5wwZycqGIxXUpWgJEAWknT7zMgBBAKBRBJDDKnoUqkzG9byZRygkBi5RyQhT7QWAi725YemgEAf9p2nwHlgFCwRaYmhM5kzJJS9zJRuWHXsCdNYHYmyuuDpTqScsPE+ne+/ZLoHTYtBaDBwFoAEpwEtmT1dTWKO1tfgRqbqHz2r33HZrWiqQxrDzqpfsMmvpArEWhXJciMBNw9QXcIfZ9BjgJ/4ChYIesrZO0EscSc6CC3KpQPIUid2W2wFopLgWIs79wcXEjPQ5HRB6KpR699qNcmAtAYKn0scZ2SHgSlbUIhjuPBYDAej0OkEUwmk+x8UaPRSLxQHnE4HMoV5Mq+miclQlOshbL6QDe5lGq3uckHlM0/ECQ2XMwjSAQACJGIbtO2X10IBQrQa2arMIkDiREQyo5fNkbmP4pEhxDNo4bIhbP+w63a94pOMrnD2oNPkm7YQkH/YHu1PBcBt0vQWbzfZ5ZCD78bMvTwBoQeHGmboxAOnQeQG3HY2NhIhh7iOJ5MJsPh0G9SZNJdYH0xuKoBA+bMAzLrA0EiykBENKOj74KRAICCZtxngFcgFAIi4hQpQdDv97NzTJQkedG2dVzT6lRdUokeslNd1BxcoBz/YGc5JAIA+QS6z6Q1PeABhEJARInoClAOi2iQ5NeKABsK6QMySQQiqAQAbIG10B0gFNzJvTyywQWZoxCoGAN/a6GUJtCjeBDRxBe214dEAMAHPO8zwC8dmhRqOBxGUVR9UUWBGMco8x8DNUNnLYQ4lhtHPvSIopuszY4b/skpnokSUzolCyUlprFemm5qtvMPKgGAEmBCqdbTCaEgJIJ4lPeeH+DWnqD7l8UVhG5gO+W8WQHIf1b7mkeKEMOOhaDTB4t5tDixApUAQHnY3mdAedoZepBzQMsBioLkAISKkc/6IplRDH8IcSBdyJCPQ2jb91sw/bHTl3IUloc/pMs2zJaXT6O73v9gXy0BAHgo6xK0FBtwpZ1CgTISIZzbX5Q4jqMoyhZu8kX2+pTSoXHJR2t5CQ1TOp3IqA9mmuVQCQCEoXH3GZBLO4XCZDKRKqFGiSDcCwcPI45j6TeUt0BSiQu1Wws/+MjjLnrWF8rv5/gDztfqAL0+EP+HSgAgKLXfZ4BH2iYUZP8q5mWo10VwGOCQMkLETnxFTJprLehQ1ISGRACAB625z4C2JTPKyRTiOOYTa0j+KeeXUq4sVEKv1xuPx+PxWFRiGI1GXoRCdkwE8+Sj1elxw7+0hTBbykLYfXcaCZWwOj0BlQBAxSDJsQW0SuuJPMFer1f7uAZBdq6H4XAoshOUARExg1Sq/SFmr8hOXV3xaXDx1YfK7+TeBz4o6R+QxkIQlRschmUawEMSAEXRXTXJ5ftf+o/hGnD0HT+By9aN9tzv5DyNrD5RavZIgS5tQoju7LsVzAcvXlT51ZXUCvde9CD5Oh2AWJYIa7PjP/jI48ocKwuEAgAOKG9lEAr8ac/9TpyCNQ6ANCCTE/v9viEgIhwR5bAFCpaVmbxKq0w+yhUK65vH5OvpntNS797zr348vUGiOKOUCETkXSUQhAIAJUjdZyAU+NOS+1242aVrp8qPJsdEVHCsh1z5iTKbH7n48t0/KpQIAggFAEqinPAWQoEnLUlm7Pf74/G4fSdBMssh6IFkMUeqKvno1hue7bztkQseuZ2imEhUXJ2e2Hv86HY9x+mJcCoBAFAe5X1GFlEN8a+uT9oCWiIUiIjJGAePDIdDqRKCxlOS9Z4l/Iu3C3Eg/sliz2vTE3vvO/rPn3xK3a0DAICW0LY6Cu0gWU0hdMGorIueshYoWNQjt/CiDjkFZaoatNjhbdcfKNkwAEAI0v6BKk0BMKQ9jkJr6Pf70khYLBZuKmE4HFqaENnrM7kkWbzdO7fc+PMOW91xwSNpZ0IpIlqbnhD/9t53dHV6HF4CAKyIEiRnqoMyaBBwFHghBj5Q6Yd4UY+hZGOSdgKfaqxJF0FYCKvT43ARAOCA8qFCGdZMPZDAVOAMhAIjRDmm8gWjfIUqsleyd7lgH32YrZ+i3AoqAYB6yYoDQ6UE5Qq7zNSLQb1AxDGiZDc8HA7jOJbzUvqtUJkqt+DxtLnsye+1X/mu8x4qX69Oj69vHiszesIZPP2AjpMUB7mywH6fZ7044PDIu/7iEreG+Z2lr4nAUWCHLiEg9xRPqgTvZO8FNk3yzup0O+4gyjHVohIA6CC6PETKKAPn28JisTj7177jtm0ggs7S1yCQzMiF8k//cRyLFKHxeOyjRWpkRpKv8ZPmmZ+Sa9571kVEtL55bH3z2Nr0BFQCAOGIlpFXPalyDmpqY3CCztLXIOAocKHf7zfiessmLlC5O8U3vnLNox5/re5doRXm6/tox0VYm55wGy4BADCTDSh0OcomZ/qVT3GiBv9gMOiaVoCjABxJPmeUtBa+8ZVrzCusTo+fOPUsgkoAwB9RBtoZuGgQB0FLsbGqzCjm80sJgvZV9rMBQgG4kCrmWL7q89rshPynXmF6AioBgJLoZIEunthZO4GIxAhznTLwmCrOHwgFUBhlMUcq96hx09d+S77WKQaoBACKoks10HmBNsqAf313L4isr+xC8aJT1gJyFDqK5czXSgyjofxWfV6bnZitnUJE8/V9UAkAWKIbu+gmC5RwcBrufs9D81fySmWz9HEDQqGLJO8XIg5X5rJXFoH2O37y5i++0Mt+AGglynhB4/IQ4zim6cX265/5glsL7f/uDzykYIuWGA6H4m7ZweGREAqdQ/gHvV5PnOuiHKSve0pyQikqeKsS0YfHPubNyYW5eY4AdBAbZaCslFxZCx2I45jof6nsWLlJBlINVDlLH08gFLqFLMokL5I4jsVNJ47j8hdAtgxLUWshmawAABBklQErwyCOY1EZtuRO6McqEgr2loCcf6eDRoIEQqFbiBM9FWDb2NgQw4JDlHxuqAsKQO2YByCwMgxS5Qv5Y3mv8zVLX9PBqIduIU76lHMg1IPf8s9Ka6ELmdIAOON9hEI1wxN8HYLbgEM5S1/HVQJhUqiuoQsEhJu7IfXQ0w5roR2fAnDAMEKhEeeYTPEjHzeQc37+26VbpOXOD19WqIXh7oqNA44CCItSlMBaAJ1FaRtQMzukOI5Ho5EoTNRWsvUrO3gHa4ZoBb6o3lHQHaLRar0pT3uAA5Z5iKFPKpmI51A6RUeq0lrJ9kdRtP/5N3tpmJJCjkJ23sgUnboDIJkRVIHyJsgtcxsAL+Qqg4rP/GR0gHZKp5Qf5ic2DzpXbY00ZZa+akDoAQRHWfJZ2q0d9PFA+0iZ0skkROX6lXVCIjpARBsbG8m5kkuOUxCJfhsbGx6LCrRVc7QAPM91C3EXyz5MODiHsl5T+SY1zlpoVmtBCMxjF93OkBDnlRjgl6oBoFxoj7Dlk/Mv+4ok3v+5Xy25BwP/8tHH47J1A45CtxAPE6lhSHLa9UK7Go1G5YczwVoATUE3dpE0VcwdDhGiGxPDnpVzJTtfv8KN4DaaEYQDOQrdQkQrRXkluVA57boZX35jbg0ZAGpEOXax3oSDonifvkg+V2TvGGKJx2RJwAQIhc7R6/Umk0m/30/FDmyubVG90W9pJkHKWqCOJRUDJliOUEjB+VzNducya6FM3HAymWTvAzJlEkKhZfAVwiAcsi6pxPI0SG2YDFL6Rd6veZ6fnJ8ggT3KUJfNwMXQJ0ByyiKP8wskP2+ZUQ/Zq14EI0Q2YhmVcNGzvuC8bS7f+5ufxmXrBu53HUXeidx8wmw2ky+y5RYYnqI8WwUsKZ+HGO4EUA7f9zUdUfKD+714fbmAEAo8Qeiho7CNIyqzFrLLAbDHoAyypoLNmRYuL0GoBJkBIMIE5SMFAtlgYQ32+30kJAIbIBQAX5ISAQ/xwB6dMsieRc4nVYizUeYJyv5bCPrBYJBKQC6JmFw+RLIRaCUQCoAjShcB1gLQYakMmI9QUDoHJZ0/X/VOzLD9SoEXIBTAUuaUjtA3miTm3DHON3pQJdmxi00foSBHJCnfjePYQTT4CltUw+r0eN1NAAogFAC7m4i5fA2yFjqLQz5B9bMqUIkLSqnXk2EIt91SRmS41VgDnQVCAQSpsBZofBf5nrAOcMY8PIFUkYUqmqXH+5zLchCEc92k8Xg8GAwGg4Gcl0FM00D8nhAAWyAUgGey47tGo5HH8V2yM0CSYytRFkOUbzmkIlZmKngfRiRnfSxz+fT7/Y2NjdEOcnn52SNBd4BQAJ6pZnwXwVpoC4Y8xNSarEYoSEKUK02q7fI9+nA4FI0sUzoFdBkIBeCTysZ3wVpoLu0YoSDxrhJk/VNfPpzcLXN9EMfx+uaxulsBFEAoAJ+EGN+lRGlKw1pgC4cRCoEeqaUmVhZVLIpUCR08jfv9Pv2nT9TdCqAAQgG40O/3lTeyEOO7cklOKEWwFnjAbYRCsj1CzjI8SUSaYbgpVABwY6XuBoBWEcdx9v7rZXyXgaR3La0F5Xw/IBzRMosMuUMYwiFOvF6vNx6Px+OxGJjA8AwRCgYqAXADjgIIS/nxXZYosxaI5YNjOyg6tVJdIxRkDkEyRiAaH8jiKolOwXThTF6bnqi7CUABHAUQkOFwKFVC0EHbymKOlBgTAbyg8wyya/IZoSBOvJROFX+yKiQAIwGwBY4CsKJomWeH8V1yc4dcM8PjLJIcy6C0DZqVBSLshNQZJUoUsJoVSZf3A0DtQCgAKwo9exUd3yULywjE6zLDx5VFoJvVvdVFVhkogwjcSiLmwjDEAEBTgFAAVtj7okXHd4mKTEQkS8yKJYPBwEv3kxwTAWtBSWrsIvN6BtVMhwgAkEAoAJ84jO+SIeRk6EEkoImKciWblPs03EGKDkBg9XV5n08B8AHJjDxBMiPwicP4LuX8NNJa8NayjLXQqSRH89hFh4GLdWWJIoIAQPXAUQD+KTS+K/SwSd3RW28tKGdXUn5kty+hyq8uxHwKAABLIBSAN9wMgGxwQWYtBIpDtzVrwWaEQkPlUXmVUEHJBAxbAG0FoQfgDXGjNJC7B2GMi3GV4abBTXWczY1EZAMKtPOJDN92EzszUfFzsViMx+Oi2wrLKqVi5exlftoHQKtp3rMFaDHJ3rqaivdJR8H+Ubuuh3KbCRSUtaeqbG2y5IZ3T0jU5yh6bmR9I7EknBgFbjzq8deG2/k3v/ob6O/cQOgBhKJojSZK3MrFGEsx/CFI4zQFnhlGIpSpBsl3s7kXVTRLRXb6xNFoFLoopw1yrrLU0EqoBABsgFAAoSjTPYhq/OGS18yP3fVG8Q1jEPhMoKBEqIReryd+epFrEjTdxBKRoDCZTJI6hpUcBIAzEAogFJZmQC31c3Iz/6u0FnTKwKNbUMEHkVH/5GSh/X5/MBiMRiPd7+tgO7kR70BONcIB6DLIUQA1o+uS7bvqOI7F8Dlf7cm1Fso/oFuOXWzQCAVd1N/8O9p02Mlf1i1HATSFxz7mzeF2/l+//sqmXE3cgKMAWJAavVYoKT0VFy9JIGvBZg6F3MZ4p8xEXClkHoDyXd3oRPT3APAHQgHUzHg8HgwGg8FAzvUg6kCTneccblhjylqggt22zQiF7AqVWQjeJ+JSdvnJMITbbgEAtYM6CqBm+v2+GOkuJoIaDAZCJdh0WkEzG5TlFsybGMokZzevd3iCnIhrPB6Px2PxE/j1ZuQgiMqKbwIAQgChAOpnOByKWjobGxui61osFrkqQfR21dTMkX1/qjSTYQ4FB0eBqppDITkRlwg6DIdD8U360l7D4VCqhNqHRwIAytCYVCkAUshwQOgRCsocQ/GCyYjKoii/MRGMUOYJFhqbkKym4L2iEZIZgTPNukhZgRwF0EhE9+NQ0NeBbEW/bDJBlekF5SkaC7C3BES5AgpmJGA+BeAMzhxnGnNrA0Ainn1lV1RNzYPUURokC2yQNkAZD0CqhDZ9MwAA5CiAhiFTE6qMfMv8A7mkmmSCLN5nQfQ1EZcYq9Lr9aASAGgZrXoqAo3GMhCum+DHcCbryjGJVD6Xtu4ct/rLx7t94msiLobTZAAAvIAcBcAFG4dAFmLKriyWKAsHxXGsnDai5NN59XkJNh5G0aLIfifi0rUQAgKA5gKhALhg0z+JHm4ymWQ7flk+KNv3y8IMqeXlPfwq+z/LUItzRKbMRFwYgwBAi0HoATSMbJ8k4+uk6ftb4IrLYYFesgV1E3G14IsCAHgHQgE0ntzuLYqiakbeB0qGIN9FI8pPxAUA6A4IPYCWI2cWln+Gm3cgUDJEoKIRZSbiAgB0BwgF0HKEUBABeLkwkMEQIhlCDDuUM2Z5oeREXACAToHQA2g5sgoQ7VQklGmP3k9+79Z9tmKxr0OkZo8UeK+4DABoARAKoOUIoZCqKCwWevcVCiVDHDx48MYbbzSv89rXvpYKFo0ohBxLqRxWCgAAhNADaC6WBQOU65QZCmhoDxVJhhAiIJeiRSMKAX0AAMgFQgE0FS/RdI+5jUWTIXJlynA4vPHGG4sWjQAAAL8g9AA6iujOPUblQyRDOBSNAAAAv0AogJajyxvwnnhYTTIESh0AACoGoQfQfiaTScmaATb5EJUlQwAAQJXAUQAtR4wwJKJszQD7k9/G5NcpCb8xDjgKAICKgVAA7afemgHekyEAAKBKIBRAVwhdM6CyZAgAAKgS5CiArlBBzYDyyRAAAMANOAoA+MFLMgQAAHADQgEAb2ACBQBA+4BQAMAzmEABANAmIBQAAAAAoGWl7gYAAAAAgC8QCgAAAADQAqEAAAAAAC0QCgAAAADQAqEAAAAAAC0QCgAAAADQAqEAAAAAAC0QCgAAAADQAqEAAAAAAC0QCgAAAADQAqEAAAAAAC0QCgAAAADQAqEAAAAAAC0QCgAAAADQAqEAAAAAAC0QCgAAAADQAqEAAAAAAC0QCgAAAADQAqEAAAAAAC0QCgAAAADQAqEAAAAAAC0QCgAAAADQAqEAAAAAAC0QCgAAAADQAqEAAAAAAC0QCgAAAADQAqEAAAAAAC0QCgAAAADQAqEAAAAAAC0QCgAAAADQAqEAAAAAAC0QCgAAAADQAqEAAAAAAC0QCgAAAADQAqEAAAAAAC0QCgAAAADQAqEAAAAAAC0QCgAAAADQAqEAAAAAAC0QCgAAAADQAqEAAAAAAC0QCgAAAADQAqEAAAAAAC0QCgAAAADQAqEAAAAAAC0QCgAAAADQAqEAAAAAAC0QCgAAAADQAqEAAAAAAC0QCgAAAADQ8v8DvqXqm5AcIlgAAAAASUVORK5CYII=\n",
      "text/plain": [
       "<IPython.core.display.Image object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "%%cpp\n",
    "\n",
    "\n",
    "  TCanvas c2;\n",
    "  c2.SetRightMargin(.2);\n",
    "  TF2 f2(\"f2\",\"exp(-x/5)*sin(y)\",-5,5,-5,5);\n",
    "  //f2.SetContour(100);\n",
    "  f2.SetNpx(100);\n",
    "  f2.SetNpy(100);\n",
    "  f2.Draw(\"surf2z\");\n",
    "  f2.GetXaxis()->SetTitle(\"X\");\n",
    "  f2.GetYaxis()->SetTitle(\"Y\");\n",
    "  f2.GetZaxis()->SetTitle(\"Z\");\n",
    "  c2.Print(\"test_style_f2.pdf\");\n",
    "  c2.Print(\"test_style_f2.png\");\n",
    "  c2.Draw();\n"
   ]
  }
 ],
 "metadata": {
  "kernelspec": {
   "display_name": "Python (sPHENIX)",
   "language": "python",
   "name": "sphenix-env"
  },
  "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.8.0"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 5
}