fig1_geometry.ipynb

{
 "cells": [
  {
   "cell_type": "code",
   "execution_count": 1,
   "id": "72b0fe04",
   "metadata": {},
   "outputs": [],
   "source": [
    "import numpy as np\n",
    "\n",
    "# plotting\n",
    "import matplotlib\n",
    "import matplotlib.pyplot as plt"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "fe98ca86",
   "metadata": {},
   "source": [
    "## Figure 1C: barplots."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 3,
   "id": "ff06ff34",
   "metadata": {},
   "outputs": [],
   "source": [
    "def plot_fig1_barplots(height1,height2):\n",
    "  fig, axs = plt.subplots(2,1,figsize=(10,10))\n",
    "\n",
    "  ax = axs[0]\n",
    "  ax.bar(x = ['$X_1$','$X_2$','','$Y_1$','$Y_2$'],\n",
    "       height=height1,\n",
    "      color=['magenta','gray','white','magenta','gray'],\n",
    "      alpha = 0.5)\n",
    "  ax.spines['top'].set_visible(False)\n",
    "  ax.spines['right'].set_visible(False)\n",
    "  ax.tick_params(axis='both', which='both',\n",
    "               bottom=False, top=False, left=False, right=False,\n",
    "               labelbottom=False, labelleft=False)\n",
    "\n",
    "  ax = axs[1]\n",
    "  ax.bar(x = ['$X_1$','$X_2$','','$Y_1$','$Y_2$'],\n",
    "       height=height2,\n",
    "      color=['magenta','gray','white','magenta','gray'],\n",
    "      alpha = 0.5)\n",
    "  ax.spines['top'].set_visible(False)\n",
    "  ax.spines['right'].set_visible(False)\n",
    "  ax.tick_params(axis='both', which='both',\n",
    "               bottom=False, top=False, left=False, right=False,\n",
    "               labelbottom=False, labelleft=False)\n",
    "\n",
    "  plt.subplots_adjust(hspace=0.4)\n",
    "\n",
    "  plt.show()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 4,
   "id": "b7e13a2a",
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAjwAAAIuCAYAAAC7EdIKAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjUuMSwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/YYfK9AAAACXBIWXMAAAsTAAALEwEAmpwYAAAKXElEQVR4nO3dMW7jMABFQWnhg7nUsXQsl75K6ugO3Gq7WEXiFZGHmZbNr4gHqtA6xlgAAMr+zB4AAPC/CR4AIE/wAAB5ggcAyBM8AECe4AEA8m5nh9u2jcfjcdUWAICfWF8dnL7wHMfx/ikAABfzSQsAyBM8AECe4AEA8gQPAJAneACAPMEDAOQJHgAgT/AAAHmCBwDIO/21xGX22QNi9tkDgFP77AEx++wB/AZeeACAPMEDAOQJHgAgT/AAAHmCBwDIEzwAQJ7gAQDyBA8AkCd4AIA8wQMA5AkeACBP8AAAeYIHAMgTPABAnuABAPIEDwCQJ3gAgDzBAwDkCR4AIE/wAAB5ggcAyBM8AECe4AEA8gQPAJAneACAPMEDAOQJHgAgT/AAAHmCBwDIEzwAQJ7gAQDyBA8AkCd4AIA8wQMA5AkeACBP8AAAeYIHAMgTPABAnuABAPIEDwCQJ3gAgDzBAwDkCR4AIE/wAAB5ggcAyBM8AECe4AEA8gQPAJAneACAPMEDAOQJHgAgT/AAAHmCBwDIEzwAQJ7gAQDyBA8AkCd4AIA8wQMA5AkeACBP8AAAeYIHAMgTPABAnuABAPIEDwCQJ3gAgDzBAwDkCR4AIE/wAAB5ggcAyBM8AECe4AEA8gQPAJAneACAPMEDAOQJHgAgT/AAAHmCBwDIEzwAQJ7gAQDyBA8AkCd4AIA8wQMA5AkeACBP8AAAeYIHAMgTPABAnuABAPIEDwCQJ3gAgDzBAwDkCR4AIE/wAAB5ggcAyBM8AECe4AEA8gQPAJAneACAPMEDAOQJHgAgT/AAAHmCBwDIEzwAQJ7gAQDyBA8AkCd4AIA8wQMA5AkeACDvNnsAv8Pz+Zw9Ied+v8+eALzgznu/2XeeFx4AIE/wAAB5ggcAyBM8AECe4AEA8gQPAJAneACAPMEDAOQJHgAgT/AAAHmCBwDIEzwAQJ7gAQDyBA8AkCd4AIA8wQMA5AkeACBP8AAAeYIHAMgTPABAnuABAPIEDwCQJ3gAgDzBAwDkCR4AIE/wAAB5ggcAyBM8AECe4AEA8gQPAJAneACAPMEDAOQJHgAgT/AAAHmCBwDIEzwAQJ7gAQDyBA8AkCd4AIA8wQMA5AkeACBP8AAAeYIHAMgTPABAnuABAPIEDwCQJ3gAgDzBAwDkCR4AIE/wAAB5ggcAyBM8AECe4AEA8gQPAJAneACAPMEDAOQJHgAgT/AAAHmCBwDIEzwAQJ7gAQDyBA8AkCd4AIA8wQMA5AkeACBP8AAAeYIHAMgTPABAnuABAPIEDwCQJ3gAgDzBAwDkCR4AIE/wAAB5ggcAyBM8AECe4AEA8gQPAJAneACAPMEDAOQJHgAgT/AAAHmCBwDIEzwAQJ7gAQDyBA8AkCd4AIA8wQMA5AkeACBP8AAAeYIHAMgTPABAnuABAPIEDwCQJ3gAgDzBAwDkCR4AIE/wAAB5ggcAyBM8AECe4AEA8gQPAJAneACAPMEDAOQJHgAgT/AAAHmCBwDIW8cYrw/X9XNZlo/r5gAAfNsxxti+OjgNHgCAAp+0AIA8wQMA5AkeACBP8AAAeYIHAMgTPABAnuABAPIEDwCQJ3gAgDzBAwDkCR4AIE/wAAB5ggcAyBM8AECe4AEA8gQPAJAneACAPMEDAOQJHgAgT/AAAHmCBwDIEzwAQJ7gAQDyBA8AkCd4AIA8wQMA5AkeACBP8AAAeYIHAMgTPABAnuABAPIEDwCQJ3gAgDzBAwDkCR4AIE/wAAB5ggcAyBM8AECe4AEA8gQPAJAneACAPMEDAOQJHgAgT/AAAHmCBwDIEzwAQJ7gAQDyBA8AkCd4AIA8wQMA5AkeACBP8AAAeYIHAMgTPABAnuABAPIEDwCQJ3gAgDzBAwDkCR4AIE/wAAB5ggcAyBM8AECe4AEA8gQPAJAneACAPMEDAOQJHgAgT/AAAHmCBwDIEzwAQJ7gAQDyBA8AkCd4AIA8wQMA5AkeACBP8AAAeYIHAMgTPABAnuABAPIEDwCQJ3gAgDzBAwDkCR4AIE/wAAB5ggcAyBM8AECe4AEA8gQPAJB3Ozvctm08Ho+rtgAA/MT66uD0hec4jvdPAQC4mE9aAECe4AEA8gQPAJAneACAPMEDAOQJHgAgT/AAAHmCBwDIEzwAQN7pryXgn+fzOXtCzv1+nz0BeMGd936z7zwvPABAnuABAPIEDwCQJ3gAgDzBAwDkCR4AIE/wAAB5ggcAyBM8AECe4AEA8gQPAJAneACAPMEDAOQJHgAgT/AAAHmCBwDIEzwAQJ7gAQDyBA8AkCd4AIA8wQMA5AkeACBP8AAAeYIHAMgTPABAnuABAPIEDwCQJ3gAgDzBAwDkCR4AIE/wAAB5ggcAyBM8AECe4AEA8gQPAJAneACAPMEDAOQJHgAgT/AAAHmCBwDIEzwAQJ7gAQDyBA8AkCd4AIA8wQMA5AkeACBP8AAAeYIHAMgTPABAnuABAPIEDwCQJ3gAgDzBAwDkCR4AIE/wAAB5ggcAyBM8AECe4AEA8gQPAJAneACAPMEDAOQJHgAgT/AAAHmCBwDIEzwAQJ7gAQDyBA8AkCd4AIA8wQMA5AkeACBP8AAAeYIHAMgTPABAnuABAPIEDwCQJ3gAgDzBAwDkCR4AIE/wAAB5ggcAyBM8AECe4AEA8gQPAJAneACAPMEDAOQJHgAgT/AAAHmCBwDIEzwAQJ7gAQDyBA8AkCd4AIA8wQMA5AkeACBP8AAAeYIHAMgTPABAnuABAPIEDwCQJ3gAgDzBAwDkCR4AIE/wAAB5ggcAyBM8AECe4AEA8gQPAJB3mz1gWZZl2WcPiNlnDwBO7bMHxOyzB/AbeOEBAPIEDwCQJ3gAgDzBAwDkCR4AIE/wAAB5ggcAyBM8AECe4AEA8gQPAJAneACAPMEDAOQJHgAgT/AAAHmCBwDIEzwAQJ7gAQDyBA8AkCd4AIA8wQMA5AkeACBP8AAAeYIHAMgTPABAnuABAPIEDwCQJ3gAgDzBAwDkCR4AIE/wAAB5ggcAyBM8AECe4AEA8gQPAJAneACAPMEDAOQJHgAgT/AAAHmCBwDIEzwAQJ7gAQDyBA8AkCd4AIA8wQMA5AkeACBP8AAAeYIHAMgTPABAnuABAPIEDwCQJ3gAgDzBAwDkCR4AIE/wAAB5ggcAyBM8AECe4AEA8gQPAJAneACAPMEDAOQJHgAgT/AAAHmCBwDIEzwAQJ7gAQDyBA8AkCd4AIA8wQMA5AkeACBP8AAAeYIHAMgTPABAnuABAPIEDwCQJ3gAgDzBAwDkCR4AIE/wAAB5ggcAyBM8AECe4AEA8gQPAJAneACAPMEDAOQJHgAgT/AAAHmCBwDIEzwAQJ7gAQDyBA8AkCd4AIA8wQMA5AkeACBP8AAAeYIHAMgTPABAnuABAPIEDwCQJ3gAgDzBAwDkCR4AIE/wAAB5ggcAyBM8AECe4AEA8gQPAJAneACAvHWM8fpwXT+XZfm4bg4AwLcdY4ztq4PT4AEAKPBJCwDIEzwAQJ7gAQDyBA8AkCd4AIC8vwtIPN0rGLuJAAAAAElFTkSuQmCC\n",
      "text/plain": [
       "<Figure size 720x720 with 2 Axes>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "\n",
    "# cis\n",
    "scale = 1\n",
    "p1 = 5\n",
    "p2 = 2.5\n",
    "p3 = 5 * scale\n",
    "p4 = p3*(p2/p1) * scale\n",
    "\n",
    "height1 = [p1,p2,0,p3,p4]\n",
    "height2 = [p2,p1,0,p4,p3]\n",
    "plot_fig1_barplots(height1,height2)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 5,
   "id": "56bca388",
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAjwAAAIuCAYAAAC7EdIKAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjUuMSwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/YYfK9AAAACXBIWXMAAAsTAAALEwEAmpwYAAAKY0lEQVR4nO3dMW7bQABFQTLwwVjyWDqWSl7FtXmHTZXOYmEr2uR5pl2A+pXwsCy4jjEWAICyX7MHAAD8bYIHAMgTPABAnuABAPIEDwCQJ3gAgLy3q8N938f9fn/VFgCA71gfHVze8Jzn+fwpAAAv5pUWAJAneACAPMEDAOQJHgAgT/AAAHmCBwDIEzwAQJ7gAQDyBA8AkHf5aYmXuc0eEHObPQC4dJs9IOb2/Ecex/H8h/5w27ZN/X03PABAnuABAPIEDwCQJ3gAgDzBAwDkCR4AIE/wAAB5ggcAyBM8AECe4AEA8gQPAJAneACAPMEDAOQJHgAgT/AAAHmCBwDIEzwAQJ7gAQDyBA8AkCd4AIA8wQMA5AkeACBP8AAAeYIHAMgTPABAnuABAPIEDwCQJ3gAgDzBAwDkCR4AIE/wAAB5ggcAyBM8AECe4AEA8gQPAJAneACAPMEDAOQJHgAgT/AAAHmCBwDIEzwAQJ7gAQDyBA8AkCd4AIA8wQMA5AkeACBP8AAAeYIHAMgTPABAnuABAPIEDwCQJ3gAgDzBAwDkCR4AIE/wAAB5ggcAyBM8AECe4AEA8gQPAJAneACAPMEDAOQJHgAgT/AAAHmCBwDIEzwAQJ7gAQDyBA8AkCd4AIA8wQMA5AkeACBP8AAAeYIHAMgTPABAnuABAPIEDwCQJ3gAgDzBAwDkCR4AIE/wAAB5ggcAyBM8AECe4AEA8gQPAJAneACAPMEDAOQJHgAgT/AAAHmCBwDIEzwAQJ7gAQDyBA8AkCd4AIA8wQMA5AkeACBP8AAAeYIHAMgTPABAnuABAPIEDwCQJ3gAgDzBAwDkCR4AIE/wAAB5ggcAyBM8AECe4AEA8gQPAJD3NnsA/4fjOGZPyNm2bfYEgB/DDQ8AkCd4AIA8wQMA5AkeACBP8AAAeYIHAMgTPABAnuABAPIEDwCQJ3gAgDzBAwDkCR4AIE/wAAB5ggcAyBM8AECe4AEA8gQPAJAneACAPMEDAOQJHgAgT/AAAHmCBwDIEzwAQJ7gAQDyBA8AkCd4AIA8wQMA5AkeACBP8AAAeYIHAMgTPABAnuABAPIEDwCQJ3gAgDzBAwDkCR4AIE/wAAB5ggcAyBM8AECe4AEA8gQPAJAneACAPMEDAOQJHgAgT/AAAHmCBwDIEzwAQJ7gAQDyBA8AkCd4AIA8wQMA5AkeACBP8AAAeYIHAMgTPABAnuABAPIEDwCQJ3gAgDzBAwDkCR4AIE/wAAB5ggcAyBM8AECe4AEA8gQPAJAneACAPMEDAOQJHgAgT/AAAHmCBwDIEzwAQJ7gAQDyBA8AkCd4AIA8wQMA5AkeACBP8AAAeYIHAMgTPABAnuABAPIEDwCQJ3gAgDzBAwDkCR4AIE/wAAB5ggcAyBM8AECe4AEA8gQPAJAneACAPMEDAOQJHgAgT/AAAHmCBwDIEzwAQJ7gAQDyBA8AkCd4AIA8wQMA5AkeACBP8AAAeYIHAMgTPABAnuABAPIEDwCQJ3gAgDzBAwDkrWOMx4fr+rEsy/vr5gAAfNk5xtg/O7gMHgCAAq+0AIA8wQMA5AkeACBP8AAAeYIHAMgTPABAnuABAPIEDwCQJ3gAgDzBAwDkCR4AIE/wAAB5ggcAyBM8AECe4AEA8gQPAJAneACAPMEDAOQJHgAgT/AAAHmCBwDIEzwAQJ7gAQDyBA8AkCd4AIA8wQMA5AkeACBP8AAAeYIHAMgTPABAnuABAPIEDwCQJ3gAgDzBAwDkCR4AIE/wAAB5ggcAyBM8AECe4AEA8gQPAJAneACAPMEDAOQJHgAgT/AAAHmCBwDIEzwAQJ7gAQDyBA8AkCd4AIA8wQMA5AkeACBP8AAAeYIHAMgTPABAnuABAPIEDwCQJ3gAgDzBAwDkCR4AIE/wAAB5ggcAyBM8AECe4AEA8gQPAJAneACAPMEDAOQJHgAgT/AAAHmCBwDIEzwAQJ7gAQDyBA8AkCd4AIA8wQMA5AkeACBP8AAAeYIHAMgTPABAnuABAPIEDwCQJ3gAgDzBAwDkCR4AIE/wAAB5ggcAyBM8AECe4AEA8gQPAJD3dnW47/u43++v2gIA8B3ro4PLG57zPJ8/BQDgxbzSAgDyBA8AkCd4AIA8wQMA5AkeACBP8AAAeYIHAMgTPABAnuABAPIuPy0BfxzHMXtCzrZtsycwy232gJjb8x/pP+/5Zv/nueEBAPIEDwCQJ3gAgDzBAwDkCR4AIE/wAAB5ggcAyBM8AECe4AEA8gQPAJAneACAPMEDAOQJHgAgT/AAAHmCBwDIEzwAQJ7gAQDyBA8AkCd4AIA8wQMA5AkeACBP8AAAeYIHAMgTPABAnuABAPIEDwCQJ3gAgDzBAwDkCR4AIE/wAAB5ggcAyBM8AECe4AEA8gQPAJAneACAPMEDAOQJHgAgT/AAAHmCBwDIEzwAQJ7gAQDyBA8AkCd4AIA8wQMA5AkeACBP8AAAeYIHAMgTPABAnuABAPIEDwCQJ3gAgDzBAwDkCR4AIE/wAAB5ggcAyBM8AECe4AEA8gQPAJAneACAPMEDAOQJHgAgT/AAAHmCBwDIEzwAQJ7gAQDyBA8AkCd4AIA8wQMA5AkeACBP8AAAeYIHAMgTPABAnuABAPIEDwCQJ3gAgDzBAwDkCR4AIE/wAAB5ggcAyBM8AECe4AEA8gQPAJAneACAPMEDAOQJHgAgT/AAAHmCBwDIEzwAQJ7gAQDyBA8AkCd4AIA8wQMA5AkeACBP8AAAeYIHAMgTPABAnuABAPIEDwCQJ3gAgDzBAwDkCR4AIE/wAAB5ggcAyBM8AECe4AEA8gQPAJAneACAvLfZA5ZlWZbb7AExt9kDAODf4oYHAMgTPABAnuABAPIEDwCQJ3gAgDzBAwDkCR4AIE/wAAB5ggcAyBM8AECe4AEA8gQPAJAneACAPMEDAOQJHgAgT/AAAHmCBwDIEzwAQJ7gAQDyBA8AkCd4AIA8wQMA5AkeACBP8AAAeYIHAMgTPABAnuABAPIEDwCQJ3gAgDzBAwDkCR4AIE/wAAB5ggcAyBM8AECe4AEA8gQPAJAneACAPMEDAOQJHgAgT/AAAHmCBwDIEzwAQJ7gAQDyBA8AkCd4AIA8wQMA5AkeACBP8AAAeYIHAMgTPABAnuABAPIEDwCQJ3gAgDzBAwDkCR4AIE/wAAB5ggcAyBM8AECe4AEA8gQPAJAneACAPMEDAOQJHgAgT/AAAHmCBwDIEzwAQJ7gAQDyBA8AkCd4AIA8wQMA5AkeACBP8AAAeYIHAMgTPABAnuABAPIEDwCQJ3gAgDzBAwDkCR4AIE/wAAB5ggcAyBM8AECe4AEA8gQPAJAneACAPMEDAOQJHgAgT/AAAHmCBwDIEzwAQJ7gAQDyBA8AkCd4AIA8wQMA5AkeACBP8AAAeYIHAMgTPABAnuABAPIEDwCQJ3gAgDzBAwDkCR4AIE/wAAB5ggcAyBM8AECe4AEA8tYxxuPDdf1YluX9dXMAAL7sHGPsnx1cBg8AQIFXWgBAnuABAPIEDwCQJ3gAgDzBAwDk/QadoTzdo5SplgAAAABJRU5ErkJggg==\n",
      "text/plain": [
       "<Figure size 720x720 with 2 Axes>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "# trans\n",
    "\n",
    "\n",
    "p1 = 5\n",
    "p2 = 2.5\n",
    "p3 = 5\n",
    "p4 = 5\n",
    "\n",
    "height1 = [p1,p2,0,p3,p4]\n",
    "height2 = [p2,p1,0,p4,p3]\n",
    "plot_fig1_barplots(height1,height2)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 6,
   "id": "44baacfc",
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAjwAAAIuCAYAAAC7EdIKAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjUuMSwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/YYfK9AAAACXBIWXMAAAsTAAALEwEAmpwYAAAKaUlEQVR4nO3dMXLbMABFQTKjg7HksXQslbyKa/MOSJXOYpE4gv2826LQrzRvgILrGGMBACj7NXsAAMD/JngAgDzBAwDkCR4AIE/wAAB5ggcAyLtdHe77Ph6Px6u2AAD8i/XZweUNz3menz8FAODFPGkBAHmCBwDIEzwAQJ7gAQDyBA8AkCd4AIA8wQMA5AkeACBP8AAAeZeflniZ++wBMffZA4BL99kDYu6zB/AduOEBAPIEDwCQJ3gAgDzBAwDkCR4AIE/wAAB5ggcAyBM8AECe4AEA8gQPAJAneACAPMEDAOQJHgAgT/AAAHmCBwDIEzwAQJ7gAQDyBA8AkCd4AIA8wQMA5AkeACBP8AAAeYIHAMgTPABAnuABAPIEDwCQJ3gAgDzBAwDkCR4AIE/wAAB5ggcAyBM8AECe4AEA8gQPAJAneACAPMEDAOQJHgAgT/AAAHmCBwDIEzwAQJ7gAQDyBA8AkHebPQAAvprjOGZPyNm2bervu+EBAPIEDwCQJ3gAgDzBAwDkCR4AIE/wAAB5ggcAyBM8AECe4AEA8gQPAJAneACAPMEDAOQJHgAgT/AAAHmCBwDIEzwAQJ7gAQDyBA8AkCd4AIA8wQMA5AkeACBP8AAAeYIHAMgTPABAnuABAPIEDwCQJ3gAgDzBAwDkCR4AIE/wAAB5ggcAyBM8AECe4AEA8gQPAJAneACAPMEDAOQJHgAgT/AAAHmCBwDIEzwAQJ7gAQDyBA8AkCd4AIA8wQMA5AkeACBP8AAAeYIHAMgTPABAnuABAPIEDwCQJ3gAgDzBAwDkCR4AIE/wAAB5ggcAyBM8AECe4AEA8gQPAJAneACAPMEDAOQJHgAgT/AAAHmCBwDIEzwAQN5t9gC+h+M4Zk/I2bZt9gSAH8MNDwCQJ3gAgDzBAwDkCR4AIE/wAAB5ggcAyBM8AECe4AEA8gQPAJAneACAPMEDAOQJHgAgT/AAAHmCBwDIEzwAQJ7gAQDyBA8AkCd4AIA8wQMA5AkeACBP8AAAeYIHAMgTPABAnuABAPIEDwCQJ3gAgDzBAwDkCR4AIE/wAAB5ggcAyBM8AECe4AEA8gQPAJAneACAPMEDAOQJHgAgT/AAAHmCBwDIEzwAQJ7gAQDyBA8AkCd4AIA8wQMA5AkeACBP8AAAeYIHAMgTPABAnuABAPIEDwCQJ3gAgDzBAwDkCR4AIE/wAAB5ggcAyBM8AECe4AEA8gQPAJAneACAPMEDAOQJHgAgT/AAAHmCBwDIEzwAQJ7gAQDyBA8AkCd4AIA8wQMA5AkeACBP8AAAeYIHAMgTPABAnuABAPIEDwCQJ3gAgDzBAwDkCR4AIE/wAAB5ggcAyBM8AECe4AEA8gQPAJAneACAPMEDAOQJHgAgT/AAAHmCBwDIEzwAQJ7gAQDyBA8AkCd4AIA8wQMA5AkeACBP8AAAeYIHAMgTPABAnuABAPIEDwCQJ3gAgDzBAwDkCR4AIE/wAAB5ggcAyBM8AECe4AEA8gQPAJAneACAPMEDAOStY4znh+v6vizL2+vmAAD8tXOMsX90cBk8AAAFnrQAgDzBAwDkCR4AIE/wAAB5ggcAyBM8AECe4AEA8gQPAJAneACAPMEDAOQJHgAgT/AAAHmCBwDIEzwAQJ7gAQDyBA8AkCd4AIA8wQMA5AkeACBP8AAAeYIHAMgTPABAnuABAPIEDwCQJ3gAgDzBAwDkCR4AIE/wAAB5ggcAyBM8AECe4AEA8gQPAJAneACAPMEDAOQJHgAgT/AAAHmCBwDIEzwAQJ7gAQDyBA8AkCd4AIA8wQMA5AkeACBP8AAAeYIHAMgTPABAnuABAPIEDwCQJ3gAgDzBAwDkCR4AIE/wAAB5ggcAyBM8AECe4AEA8gQPAJAneACAPMEDAOQJHgAgT/AAAHmCBwDIEzwAQJ7gAQDyBA8AkCd4AIA8wQMA5AkeACBP8AAAeYIHAMgTPABAnuABAPIEDwCQJ3gAgDzBAwDkCR4AIE/wAAB5ggcAyBM8AECe4AEA8gQPAJAneACAPMEDAOQJHgAgT/AAAHmCBwDIEzwAQJ7gAQDyBA8AkHe7Otz3fTwej1dtAQD4F+uzg8sbnvM8P38KAMCLedICAPIEDwCQJ3gAgDzBAwDkCR4AIE/wAAB5ggcAyBM8AECe4AEA8i4/LQF/HMcxe0LOtm2zJwBP+M/7fLP/89zwAAB5ggcAyBM8AECe4AEA8gQPAJAneACAPMEDAOQJHgAgT/AAAHmCBwDIEzwAQJ7gAQDyBA8AkCd4AIA8wQMA5AkeACBP8AAAeYIHAMgTPABAnuABAPIEDwCQJ3gAgDzBAwDkCR4AIE/wAAB5ggcAyBM8AECe4AEA8gQPAJAneACAPMEDAOQJHgAgT/AAAHmCBwDIEzwAQJ7gAQDyBA8AkCd4AIA8wQMA5AkeACBP8AAAeYIHAMi7zR4A8OPcZw+Iuc8ewHfghgcAyBM8AECe4AEA8gQPAJAneACAPMEDAOQJHgAgT/AAAHmCBwDIEzwAQJ7gAQDyBA8AkCd4AIA8wQMA5AkeACBP8AAAeYIHAMgTPABAnuABAPIEDwCQJ3gAgDzBAwDkCR4AIE/wAAB5ggcAyBM8AECe4AEA8gQPAJAneACAPMEDAOQJHgAgT/AAAHmCBwDIEzwAQJ7gAQDyBA8AkCd4AIA8wQMA5AkeACBP8AAAeYIHAMgTPABAnuABAPIEDwCQJ3gAgDzBAwDkCR4AIE/wAAB5ggcAyBM8AECe4AEA8gQPAJAneACAPMEDAOQJHgAgT/AAAHmCBwDIEzwAQJ7gAQDyBA8AkCd4AIA8wQMA5AkeACBP8AAAebfZA5ZlWZb77AEx99kDAOBrccMDAOQJHgAgT/AAAHmCBwDIEzwAQJ7gAQDyBA8AkCd4AIA8wQMA5AkeACBP8AAAeYIHAMgTPABAnuABAPIEDwCQJ3gAgDzBAwDkCR4AIE/wAAB5ggcAyBM8AECe4AEA8gQPAJAneACAPMEDAOQJHgAgT/AAAHmCBwDIEzwAQJ7gAQDyBA8AkCd4AIA8wQMA5AkeACBP8AAAeYIHAMgTPABAnuABAPIEDwCQJ3gAgDzBAwDkCR4AIE/wAAB5ggcAyBM8AECe4AEA8gQPAJAneACAPMEDAOQJHgAgT/AAAHmCBwDIEzwAQJ7gAQDyBA8AkCd4AIA8wQMA5AkeACBP8AAAeYIHAMgTPABAnuABAPIEDwCQJ3gAgDzBAwDkCR4AIE/wAAB5ggcAyBM8AECe4AEA8gQPAJAneACAPMEDAOQJHgAgT/AAAHmCBwDIEzwAQJ7gAQDyBA8AkCd4AIA8wQMA5AkeACBP8AAAeYIHAMgTPABAnuABAPIEDwCQJ3gAgDzBAwDkCR4AIE/wAAB5ggcAyBM8AECe4AEA8gQPAJAneACAPMEDAOQJHgAgT/AAAHmCBwDIEzwAQJ7gAQDyBA8AkCd4AIA8wQMA5AkeACBP8AAAeesY4/nhur4vy/L2ujkAAH/tHGPsHx1cBg8AQIEnLQAgT/AAAHmCBwDIEzwAQJ7gAQDyfgOk7jzhOB1pCwAAAABJRU5ErkJggg==\n",
      "text/plain": [
       "<Figure size 720x720 with 2 Axes>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "# cis + trans\n",
    "p1 = 5\n",
    "p2 = 2.5\n",
    "p3 = 5\n",
    "p4 = 4.0\n",
    "\n",
    "height1 = [p1,p2,0,p3,p4]\n",
    "height2 = [p2,p1,0,p4,p3]\n",
    "plot_fig1_barplots(height1,height2)\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 7,
   "id": "d3c69252",
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAjwAAAIuCAYAAAC7EdIKAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjUuMSwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/YYfK9AAAACXBIWXMAAAsTAAALEwEAmpwYAAAKaklEQVR4nO3du3HjMABFUXBHhTFUWSpLIVtxbPWAbWDFwB9hfX1OikAv0twBAm5zzgEAUPZn9QAAgO8meACAPMEDAOQJHgAgT/AAAHmCBwDIu5wdXq/Xeb/fX7UFAOAztmcHpzc8j8fj66cAALyYJy0AIE/wAAB5ggcAyBM8AECe4AEA8gQPAJAneACAPMEDAOQJHgAg7/TTEgDwGx3HsXpCzr7vS3/fDQ8AkCd4AIA8wQMA5AkeACBP8AAAeYIHAMgTPABAnuABAPIEDwCQJ3gAgDzBAwDkCR4AIE/wAAB5ggcAyBM8AECe4AEA8gQPAJAneACAPMEDAOQJHgAgT/AAAHmCBwDIEzwAQJ7gAQDyBA8AkCd4AIA8wQMA5AkeACBP8AAAeYIHAMgTPABAnuABAPIEDwCQJ3gAgDzBAwDkCR4AIE/wAAB5ggcAyBM8AECe4AEA8gQPAJAneACAPMEDAOQJHgAgT/AAAHmCBwDIEzwAQJ7gAQDyBA8AkCd4AIA8wQMA5AkeACBP8AAAeYIHAMgTPABAnuABAPIEDwCQJ3gAgDzBAwDkCR4AIE/wAAB5ggcAyBM8AECe4AEA8gQPAJAneACAPMEDAOQJHgAgT/AAAHmCBwDIEzwAQJ7gAQDyBA8AkCd4AIA8wQMA5AkeACBP8AAAeYIHAMgTPABAnuABAPIEDwCQJ3gAgDzBAwDkCR4AIE/wAAB5ggcAyBM8AECe4AEA8gQPAJAneACAPMEDAOQJHgAgT/AAAHmCBwDIEzwAQJ7gAQDyBA8AkCd4AIA8wQMA5AkeACBP8AAAeYIHAMgTPABAnuABAPIEDwCQJ3gAgDzBAwDkXVYPGGOMcVs9IOa2egBw6rZ6QMxt9QB+Ajc8AECe4AEA8gQPAJAneACAPMEDAOQJHgAgT/AAAHmCBwDIEzwAQJ7gAQDyBA8AkCd4AIA8wQMA5AkeACBP8AAAeYIHAMgTPABAnuABAPIEDwCQJ3gAgDzBAwDkCR4AIE/wAAB5ggcAyBM8AECe4AEA8gQPAJAneACAPMEDAOQJHgAgT/AAAHmCBwDIEzwAQJ7gAQDyBA8AkCd4AIA8wQMA5AkeACBP8AAAeYIHAMgTPABAnuABAPIEDwCQJ3gAgDzBAwDkCR4AIE/wAAB5ggcAyBM8AECe4AEA8gQPAJAneACAPMEDAORdVg/gZziOY/WEnH3fV08A+DXc8AAAeYIHAMgTPABAnuABAPIEDwCQJ3gAgDzBAwDkCR4AIE/wAAB5ggcAyBM8AECe4AEA8gQPAJAneACAPMEDAOQJHgAgT/AAAHmCBwDIEzwAQJ7gAQDyBA8AkCd4AIA8wQMA5AkeACBP8AAAeYIHAMgTPABAnuABAPIEDwCQJ3gAgDzBAwDkCR4AIE/wAAB5ggcAyBM8AECe4AEA8gQPAJAneACAPMEDAOQJHgAgT/AAAHmCBwDIEzwAQJ7gAQDyBA8AkCd4AIA8wQMA5AkeACBP8AAAeYIHAMgTPABAnuABAPIEDwCQt805nx9u2/sY4+11cwAAPuwx57z+6+A0eAAACjxpAQB5ggcAyBM8AECe4AEA8gQPAJAneACAPMEDAOQJHgAgT/AAAHmCBwDIEzwAQJ7gAQDyBA8AkCd4AIA8wQMA5AkeACBP8AAAeYIHAMgTPABAnuABAPIEDwCQJ3gAgDzBAwDkCR4AIE/wAAB5ggcAyBM8AECe4AEA8gQPAJAneACAPMEDAOQJHgAgT/AAAHmCBwDIEzwAQJ7gAQDyBA8AkCd4AIA8wQMA5AkeACBP8AAAeYIHAMgTPABAnuABAPIEDwCQJ3gAgDzBAwDkCR4AIE/wAAB5ggcAyBM8AECe4AEA8gQPAJAneACAPMEDAOQJHgAgT/AAAHmCBwDIEzwAQJ7gAQDyBA8AkCd4AIA8wQMA5AkeACBP8AAAeYIHAMgTPABAnuABAPIEDwCQJ3gAgDzBAwDkCR4AIE/wAAB5ggcAyBM8AECe4AEA8gQPAJAneACAPMEDAOQJHgAgT/AAAHmCBwDIEzwAQJ7gAQDyBA8AkCd4AIA8wQMA5AkeACDvcnZ4vV7n/X5/1RYAgM/Ynh2c3vA8Ho+vnwIA8GKetACAPMEDAOQJHgAgT/AAAHmCBwDIEzwAQJ7gAQDyBA8AkCd4AIC8009LAPANbqsHxNxWD+AncMMDAOQJHgAgT/AAAHmCBwDIEzwAQJ7gAQDyBA8AkCd4AIA8wQMA5AkeACBP8AAAeYIHAMgTPABAnuABAPIEDwCQJ3gAgDzBAwDkCR4AIE/wAAB5ggcAyBM8AECe4AEA8gQPAJAneACAPMEDAOQJHgAgT/AAAHmCBwDIEzwAQJ7gAQDyBA8AkCd4AIA8wQMA5AkeACBP8AAAeYIHAMgTPABAnuABAPIEDwCQJ3gAgDzBAwDkCR4AIE/wAAB5ggcAyBM8AECe4AEA8gQPAJAneACAPMEDAOQJHgAgT/AAAHmCBwDIEzwAQJ7gAQDyBA8AkCd4AIA8wQMA5AkeACBP8AAAeYIHAMgTPABAnuABAPIEDwCQJ3gAgDzBAwDkCR4AIE/wAAB5ggcAyBM8AECe4AEA8gQPAJAneACAPMEDAOQJHgAgT/AAAHmCBwDIEzwAQJ7gAQDyBA8AkCd4AIA8wQMA5AkeACBP8AAAeYIHAMgTPABAnuABAPIEDwCQJ3gAgDzBAwDkCR4AIE/wAAB5ggcAyBM8AECe4AEA8gQPAJAneACAPMEDAOQJHgAgT/AAAHmCBwDIEzwAQJ7gAQDyBA8AkCd4AIA8wQMA5AkeACBP8AAAeZfVA/gZjuNYPSFn3/fVE4An/Od9vdX/eW54AIA8wQMA5AkeACBP8AAAeYIHAMgTPABAnuABAPIEDwCQJ3gAgDzBAwDkCR4AIE/wAAB5ggcAyBM8AECe4AEA8gQPAJAneACAPMEDAOQJHgAgT/AAAHmCBwDIEzwAQJ7gAQDyBA8AkCd4AIA8wQMA5AkeACBP8AAAeYIHAMgTPABAnuABAPIEDwCQJ3gAgDzBAwDkCR4AIE/wAAB5ggcAyBM8AECe4AEA8gQPAJAneACAPMEDAOQJHgAgT/AAAHmCBwDIEzwAQJ7gAQDyBA8AkCd4AIA8wQMA5AkeACBP8AAAeYIHAMi7rB4wxhjjtnpAzG31AAD4v7jhAQDyBA8AkCd4AIA8wQMA5AkeACBP8AAAeYIHAMgTPABAnuABAPIEDwCQJ3gAgDzBAwDkCR4AIE/wAAB5ggcAyBM8AECe4AEA8gQPAJAneACAPMEDAOQJHgAgT/AAAHmCBwDIEzwAQJ7gAQDyBA8AkCd4AIA8wQMA5AkeACBP8AAAeYIHAMgTPABAnuABAPIEDwCQJ3gAgDzBAwDkCR4AIE/wAAB5ggcAyBM8AECe4AEA8gQPAJAneACAPMEDAOQJHgAgT/AAAHmCBwDIEzwAQJ7gAQDyBA8AkCd4AIA8wQMA5AkeACBvm3M+P9y29zHG2+vmAAB82GPOef3XwWnwAAAUeNICAPIEDwCQJ3gAgDzBAwDkCR4AIO8vXAs84SfkXYgAAAAASUVORK5CYII=\n",
      "text/plain": [
       "<Figure size 720x720 with 2 Axes>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "# cis x trans\n",
    "p1 = 5\n",
    "p2 = 2.5\n",
    "p3 = 5\n",
    "p4 = 10\n",
    "\n",
    "\n",
    "height1 = [p1,p2,0,p3,p4]\n",
    "height2 = [p2,p1,0,p4,p3]\n",
    "plot_fig1_barplots(height1,height2)"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "207f72a5",
   "metadata": {},
   "source": [
    "## Figure 1D: theoretical geometry/diagram. "
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 8,
   "id": "4802e445",
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAABBwAAAFUCAYAAACHqC+jAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjUuMSwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/YYfK9AAAACXBIWXMAAAsTAAALEwEAmpwYAAA8K0lEQVR4nO3deZhc9X3n+89Zqk5Vd1W1drUkBIhVAknYxhhjjMEsQgYbDEIkM5ncuUkmzuRmZhwndq4TOzAo3m7imVxPJjf3Zhk/MxNn/CAQizGW2MwajDG2kQQSILEJSd3au/Zz6pzzu38IUFefQmhpne6ufr/+MvVRoa9b6nrcb/diGWMEAAAAAAAwmuyxPgAAAAAAAHQfggMAAAAAABh1BAcAAAAAADDqCA4AAAAAAGDUERwAAAAAAMCoIzgAAAAAAIBR5x5u/O5z3zV763vTugWTzN//0d+rGlT1he98YaxPwQQ1fdeAPvboOjlRmNhiI33j2i899N1rf+OqMThtTPHajdFyw/f/W9s///YPt8myMrr9po+2Pb7xC6vTPAtd7unBxu9/+YMzvjPWd6Tpt9f89rs/p37b0DYFUdC2X3P2NanfNN787LGPJx6znT6VclPaHvsXv1JK6aLxafN+/4h+3cKp3gm+ZPxb/JcrE4+V/QOKowNtj931q7+Z0kXj1/0v3Z94LOtkNb9vviTp7278O+u9nnvY4LC3vlfzSvOO9z6gI8/1VPbLmlucO9anYAIq7BnU0qeflGu5ktv+UhYZoyDO6KTpF+0ao/PGFK/dGC05N9f2z45lK5KRa7f/74q+rJPmWeh+M8f6gLEUREHifY/XdGnjiLeJJDViI9fi9Wi4ka/P72Wyv52kzm+rKDbK8/6XMPI1SZKaYfOInsuXVACYcPIH9mrJA6vlBsmKHxujIHbVmy2p15vx2hicBwAAAEAEBwATjFcd0tJ1q5VpNhJbLMmPXfVki+rJ8KmCAAAAwFgiOACYMDKNmpauXS2vVklsRpIf2cpniurN5NM/DgAAAEAbggOACcH1m1q6brXy5f2JzUhqRpa8TFGFbPJrzAAAAACkj+AAYNyzW4EWP3inevftTmxGUjOUPLeoYjYv6ci+WRIAAACAE4vgAGBcs6JQ5z58t0q7diS2g5/ZIGXdoopej6xhsaE1pT/FKwEAAACMRHAAMG5ZcaxFj/1QU3e80XFvRkYZp6BSrj02RMVpGrj51rTOBAAAANABwQHA+GSMznxqnWa8/nLHuRkZuU5BJa9X1rCXsjhf0M6Vtyjsm5XWpQAAAAA6IDgAGH+M0Wk/fVT9r2zsODejWLbdq5LXK9saFhuyOe286Ra1ps9P61IAAAAA74HgAGDcOfn5p3XSCz/ruPlRLNsuqM8ryLGcdx83bkaDN/6Jgv7T0zoTAAAAwGEQHACMK3NffE6n/vypjlsQxZLdq5JXkGMfig2ybQ1e90U155+b0pUAAAAA3g/BAcC4MfuVjTrjJ4903FpxrNjqUckryB0eGyxLu675vBqnfzilKwEAAAAcCYIDgHFh+usv66wn13bcwjhWpLxKXkEZ223b9lz1OdUWfTyNEwEAAAAcBYIDgDE3ZfvrWvTofbKMSWyhMWqZnIpeUVkn07bt+8SvqXLesrTOBAAAAHAUCA4AxlRx13ad+8jdsuMosUXGqBVnVfSK8kbEhgMX3qChC29M60wAAAAAR4ngAGDM9O7brSUP3Cmn1UpskTEK4owKXlE5N9u2lT+wTPsv+bW0zgQAAABwDAgOAMZErrxfS9atlhv4iS02RkHsqjdbVN712rbqwou198rPSZaV1qkAAAAAjgHBAUDqsrWylq5drWyjlthiGfmxo3y2qJ5Mrm2rn36+dl/zH4gNAAAAwARAcACQqkyjrqVrVytXHUpsRpIf2cpliipk8m1bc/452nXdFyXHTTwPAAAAwPhDcACQGidoaskDd6hnaF9iM5KaoSXPLaqYbY8Nfv9pGrjhj2VGfC8HAAAAAOMXwQFAKuww1OKH7lZh72BiM5KakZTNFFT0eiQd+pKJ1vSTNLDiT2W8nvSOBQAAAHDcCA4ATjgrirTox/eqb2Bbx70ZSRmnqJLXI2tYbAj7ZmrnylsU95TSOhUAAADAKCE4ADix4lhnP3G/pm/b2nFuRkauU1Cf1yNr2EtS1DtFO1feqqg4Pa1LAQAAAIwiggOAE8cYnfn0Q5r16uaOczMycpyC+rxeWdahl6M416uBlbconDonrUsBAAAAjDKCA4ATZsFzT2jOS8933Pwolm33qpTtlT0sNpiMp4EVX1Uw85S0zgQAAABwAhAcAJwQ89c/o/nrn+m4+VEs2b3q8wpybOfdx43jaOCGL8ufe1ZaZwIAAAA4QQgOAEbdnM2/1IKfPd5xC6JYxupJxAbZtnZ95g/VPGVpSlcCAAAAOJEIDgBG1cytm3TG0w913FqxUWzl1ZcryrXdtm338t9T/cwL0zgRAAAAQAoIDgBGzbRtW7Xw8ftlGZPYwtgolKeSV1RmRGzYe/lvqnruZSldCQAAACANBAcAo6Jv55s655F7ZZk4sUXGqGWyKnlFZZ1M27b/4l9R+fxr0zoTAAAAQEoIDgCOW2H3Ti1+aI3sKExskTEK4oyKXkmek23bhj78aR24aGVaZwIAAABIEcEBwHHpObBXSx64U06rldjit2NDwSsp57bHhsqSy7Xvsv9dsqyULgUAAACQJoIDgGOWqwxp6drblfEbiS2W5MeuerJF5V2vbaud/VHtWfa7xAYAAACgixEcAByTTL2qJetuV7ZeTWxGkh/ZymeL6s3k2rbGgg9o1zW/L9m8/AAAAADdjP/FD+Couc2Glq5brXz5QGIzkpqRpVymqEIm37Y1552tweu+JLmZxPMAAAAAdBeCA4CjYrcCLX5wjXr370lsB2OD5LlFFbPtsSGYdYoGb/yKTDaXeB4AAACA7kNwAHDE7DDUuQ/frdLuHYntndiQdYsqej2SDn1/htbUORq46RbFud70jgUAAAAwpggOAI6IFcda+Nh9mrrjjY57MzLKOAWVvB5Zw2JDWJyugZW3KOqdktKlAAAAAMYDggOA92eMznpyrWa88UrHuRkZuU5BJa9X1rCXlainpIGbb1XYNyutSwEAAACMEwQHAIdnjE5/5hHN3vJCx7kZxbKdXvV5BdnWoZeU2Mtr4KZb1Jo2L61LAQAAAIwjBAcAh3XKL57SvBd/3nHzo1i2XVBftj02GDejgRu/omD2grTOBAAAADDOEBwAvKd5G3+mU375dMctiGLJPviZDY7tvPu4cRwNXv9H8k9alNaZAAAAAMYhggOAjvpf3qDTf/rjjlsrjhVbPSqNiA2yLO2+5vNqnPahlK4EAAAAMF4RHAAkzHj9ZZ351LqOWxgbRcqr5BWVsd22bc+yf6vawovTOBEAAADAOEdwANBm6luvaeGj98kyJrGFsVHLeCp5RWWd9tiw77L/TZWlV6Z1JgAAAIBxjuAA4F2lwe069+G7ZcdRYouMUctkVfSKyjqZtu3ARTdp6ILr0zoTAAAAwARAcAAgSerdO6jFD94pOwoTW2SMgjijgldUzs22beUPfUr7L/7VtM4EAAAAMEEQHAAoP7RPS9fdITfwE1tsjILYVW+2qLzrtW3Vcz6hvZf/lmRZaZ0KAAAAYIIgOACTnFcta+na25Vp1hNbLCM/dtSTLaonk2vb6mdcoN2f+nfEBgAAAAAdERyASSzTqGvJutXyapXEZiT5ka18pqTeTL5ta5y8WLs+8wfS8B+JCQAAAADDEByAScoJmlrywGr1DO1LbEZSM7LkuSUVsu2f2eDPOUODN3xZZsT3cgAAAACA4QgOwCRkhy0teWCNCnt3JbaDsUHKukUVvbykQ18yEcyYr4EVX5XJ5hPPAwAAAIDhCA7AJGNFkc555B6Vdm3vuDcjKesUVfJ6ZA2LDeGU2RpYeYvifDGtUwEAAABMYAQHYDKJYy18/Iea9tZrHedGZJRxConYEBWmaufNtyoqTEvrUgAAAAATHMEBmCyM0Vn//IBmvvZSx7kZGblOQSWvV5Z16KUhzhe0c+UtCvtmp3UpAAAAgC5AcAAmA2N02rOPqf/lDR1nP4pl270qeb2yh8UGk/G0c8WfqjXj5LQuBQAAANAlCA7AJDB//TM6aeOzHTc/iiW7V31eQY516MdcGtfVwI1/omDOGWmdCQAAAKCLEByALjd30y+04LknOm5BFEtWr/q8ohz7UGyQbWvXZ76o5smLU7oSAAAAQLchOABdbNaWF3TG0w913FpxrNjKq5QryB0eGyxLuz/171U/44KUrgQAAADQjQgOQJea/uYWnf3E2o5bGMcKlVfJKypju23bniv+jarnfCKNEwEAAAB0MYID0IWm7HhTi378A1kmTmyhMWqZnEpeUVkn07btu+RfqvLB5WmdCQAAAKCLERyALlPctUPnPrxGdhQmtsgYteKMil5R3ojYMHTB9Rq68Ma0zgQAAADQ5QgOQBfp2b9bSx68U06rldhiYxTEGRW8knJutm2rLL1S+y79dcmy0joVAAAAQJcjOABdIlfer6Xr7pDrNxNbbIz82FVPtqi867VttYUf056rfofYAAAAAGBUERyALpCtVbR03Wpl69XEZiT5saN8tqjeTK5tayz4oHZd83nJ5qUAAAAAwOjiowxggnObDS1dt1q5ylBiMzJqRpZymaIKmXzb1jxpkQav/5LkuInnAQAAAMDxIjgAE5jTCrTkgTvUc2BvYjOSmpEtzy2qmG2PDcHsBRq84Y9lMl7ieQAAAAAwGggOwARlh6HOfWiNinsGEtvB2CBl3YKKXo+kQ9+foTVtnnbe9KeKc73pHQsAAABg0iE4ABOQFUVa9ON7NWXnto57MzLKOEWVvB5Zw2JD2DdTO1feorinL61TAQAAAExSBAdgojFGZz/xI03ftrXj3IiM3Hdjw6F38ai3TztX3qKoNCOtSwEAAABMYgQHYCIxRmc8/ZBmvbqp49yMYjlOQX1er2zr0Lt37PVoYOWtCqfOTetSAAAAAJMcwQGYQE79+ZOau/mXHTc/imXbBfVl22ODyWQ1sOKrCmaektKVAAAAAEBwACaMkzb8VCc//5OOWxDFkt2rPq8gx3befdw4jgY/+2X5885O60wAAAAAkERwACaE/pee12nPPtZxC+JYsdWTiA2yLO369BfUOPW8lK4EAAAAgEMIDsA4N/PVzTrznx/suLXiWLHy6ssV5dpu27Z7+e+pftZFaZwIAAAAAAkEB2Acm7btVS18/IeyjElsYWwUmpxKXlGZEbFh7+W/oeriT6Z1JgAAAAAkEByAcao08JbOeeQeWXGc2CJj1DJZlbyisk6mbdv/sZUqn//ptM4EAAAAgI4IDsA4VNgzqMUPrZEdhYktMkZBnFHBK8lzs23b0PnX6sDHfiWtMwEAAADgPREcgHEmf2CvljywWm7gJ7bYGAWxq95sSfkRsaG6+DLt++RvSJaV1qkAAAAA8J4IDsA44lWHtHTdamWajcQWS/JjVz3ZonoyXttWO/NC7b76/yA2AAAAABg3CA7AOJFp1LR07e3yapXEZiT5ka18pqjeTL5ta5yyVLs+/QVp+I/EBAAAAIAxRnAAxgHXb2rJujuULx9IbEZSM7KUyxRVyObaNn/umRr87P8puZnE8wAAAABgLBEcgDFmtwItfvBOFfbtSmzvxAbPLaqQzUs69CUTwcyTNbDiqzIjIgQAAAAAjAcEB2AMWVGocx++W6VdOxLbwdggZZ2Cil6PrGGxoTWlXwM33aI4V0jxWgAAAAA4cgQHYIxYcaxFj/1QU3e80XFvRkYZp6BSrj02RIWpGrj5VkWFqWmdCgAAAABHjeAAjAVjdNZT6zTj9Zc7zs3IyHUKKnm9soa9m8b5gnbefKvCvllpXQoAAAAAx4TgAKTNGJ3200c1+5WNHedmFMu2e1XyemVbw2JDNqedN92i1vT5aV0KAAAAAMeM4ACk7JRfPq2TXvhZx82PYtl2QX1eQY516MdcGjejwRv/REH/6WmdCQAAAADHheAApGjeC8/plF881XELoliye1XyCnLsQ7FBtq3B67+k5vxzU7oSAAAAAI4fwQFIyexXNur0Zx7puLXiWLHVo5JXkDs8NliWdl3zeTVOOz+lKwEAAABgdBAcgBRMf/1lnfXk2o5bGMeKlFfJKyhju23bnqs+p9qij6dxIgAAAACMKoIDcIJN2f66Fj16nyxjEltojFomp6JXVNbJtG37Lv1Xqpy3LK0zAQAAAGBUERyAE6i4a7vOfeRu2XGU2CJj1IqzKnpFeSNiw4ELb9DQR25I60wAAAAAGHUEB+AE6d23W0seuFNOq5XYImMUxBkVvKJybrZtK39gmfZf8mtpnQkAAAAAJwTBATgB8kP7tWTdarmBn9hiYxTErnqzReVdr22rLvq49l75Ocmy0joVAAAAAE4IggMwyrK1spasW61so5bYYhn5saN8tqieTK5tq59+vnZ/6t8TGwAAAAB0BYIDMIoyjbqWrl2tXHUosRlJfmQrlymqkMm3bc3552jXdV+UHDfxPAAAAACYiAgOwChxgqaWPHCHeob2JTYjqRlZ8tyiitn22OD3n66BG/5YZsT3cgAAAACAiYzgAIwCOwy1+KG7Vdg7mNgOxgYp6xZU9HokHfqSidb0kzSw4qsyXk96xwIAAABACggOwHGyokiLfnyv+ga2ddybkZRxiip5PbKGxYawb6Z2rrxFcU8prVMBAAAAIDUEB+B4xLEWPn6/pm/b2nFuRkauU1Cf1yNr2Ltb1DtFO1feqqg4Pa1LAQAAACBVBAfgWBmjM59+SDNf29xxbkZGjlNQn9cryzr0rhbnejWw8haFU+ekdSkAAAAApI7gAByjBc89oTkvPd9x86NYtt2rktcre1hsMBlPAyu+qmDmKWmdCQAAAABjguAAHIP565/R/PXPdNz8KJbsXvV5BTmW8+7jxnE0cMOX5c89K60zAQAAAGDMEByAozRn0y+04GePd9yCKJaxeg7GBvtQbJBta9dn/lDNU5amdCUAAAAAjC2CA3AUZm7dpDN+8nDHrRUbxVZefbmiXNtt23Yv/z3Vz7wwjRMBAAAAYFwgOABHaNq2rVr4+P2yjElsYWwUylPJKyozIjbsvfw3VT33spSuBAAAAIDxgeAAHIG+nW/qnEfulWXixBYZo5Y5GBuyTqZt2//xX1X5/GvTOhMAAAAAxg2CA/A+Crt3avFDa2RHYWKLjFEQZ1T0ivKcbNs29OHP6MBHb0rrTAAAAAAYVwgOwGH0HNirJQ/cKafVSmzx27Gh4JWUc9tjQ2XJ5dp32b+WLCutUwEAAABgXCE4AO8hVxnS0rW3K+M3ElssyY9d9WSLyrte21Y7+yLtWfa7xAYAAAAAkxrBAeggW69qydrbla1XE5uR5Ee28tmiejO5tq2x4APade3vSzbvWgAAAAAmNz4qAkZwmw0tWbda+cqBxGYkNSNLuUxRhUy+bWvOO1uD131JctzE8wAAAABgsiE4AMPYrUCLH1yj3v17EtvB2CB5blHFbHtsCGadqsEbvyKTzSWeBwAAAACTEcEBeJsdhjr34btV2r0jsb0TG7JuUUWvR9Kh78/QmjZXO1feojjXm96xAAAAADDOERwASVYca+Fj92nqjjc67s3IKOMUVPJ6ZA2LDWFx+sHY0NOX1qkAAAAAMCEQHABjdNaTazXjjVc6zs3IyHUKKnm9soa9y0Q9JQ3cfKui0sy0LgUAAACACYPggMnNGJ3+zCOaveWFjnMzimU7BfV5BdnWoXeX2Mtr4KZb1Jo2L61LAQAAAGBCIThgUjvlF09p3os/77j5USzbLqgv29sWG4yb0cCNX1Ewe0FaZwIAAADAhENwwKQ1b+OzOuWXT3fcgiiW7F71eQU5tvPu48ZxNHj9H8k/aVFaZwIAAADAhERwwKTU//IGnf7TRzturThWbPWoNCI2yLK0+5rPq3Hah9I5EgAAAAAmMIIDJp0Zr72kM59a13FrxbEi5VXyisrYbtu2++rfVW3hxWmcCAAAAAATHsEBk8rUt17Twsd+KMuYxBbGRqHJqeQVlXXaY8O+y/61qkuuSOtMAAAAAJjwCA6YNEqD23Xuw3fLjqPEFhmjlsmq6BWVdTJt24GLbtLQBdeldSYAAAAAdAWCAyaF3r2DWvzgnbKjMLFFxiiIMyp4ReXcbNtW/tCntP/iX03rTAAAAADoGgQHdL380D4tXXeH3MBPbLExCmJXvdmi8q7XtlXP+YT2Xv5bkmWldSoAAAAAdA2CA7qaVy1r6drblWnWE1ssIz921JMtqieTa9vqZ1yg3Z/6d8QGAAAAADhGBAd0rUyjriXrVsurVRKbkeRHtvKZknoz+batcfJi7frMH0jDfyQmAAAAAOCoEBzQlZygqSUPrFbP0L7EZiQ1I0ueW1Ih2/6ZDf6cMzR4w5dlRnwvBwAAAADA0SE4oOvYYUtLHlijwt5die1gbJCyblFFLy/p0JdMBDPma2DFV2Wy+cTzAAAAAABHh+CArmJFkc55+B6Vdm3vuDcjKesUVfJ6ZA2LDeGU2RpYeYvifDGtUwEAAACgqxEc0D3iWAsf/6GmbX+t49yIjDJOIREbosJU7bz5VkWFaWldCgAAAABdj+CA7mCMzvrnBzTztZc6zs3IyHUKKnm9sqxDf+3jfEE7b75VYd/stC4FAAAAgEmB4ICJzxid9uxj6n95Q8fZj2LZdq9KXq/sYbHBZDztXPGnak2fn9alAAAAADBpEBww4c1f/4xO2vhsx82PYskuqM8ryLEO/ZhL47oauPFPFMw5I60zAQAAAGBSIThgQpu76Rda8NwTHbcgiiWr92BssA/FBtm2Bq/7kponL07pSgAAAACYfAgOmLBmbXlBZzz9UMetFceKrbxKuYLc4bHBsrTrmv+gxukfTulKAAAAAJicCA6YkKa/uUVnP7G24xbGsULlVfKKythu27bnin+j2qJL0jgRAAAAACY1ggMmnCk73tSiR+6VZeLEFhqjlsmp5BWVdTJt275L/qUqH1ye1pkAAAAAMKkRHDChFHft0LkPr5EdR4ktMkatOKOiV5Q3IjYMXXC9hi68Ma0zAQAAAGDSIzhgwujZv1tLHrxTTquV2GJjFMQZFbyScm62baucd5X2XfrrkmWldSoAAAAATHoEB0wIufJ+LV27Wq7fTGyxMfJjVz3ZovKu17bVFn5Me678HLEBAAAAAFJGcMC4l61VtHTdamUbtcQWy8iPHeWzRfVmcm1bY8EHteuaz0s2f80BAAAAIG18JIZxzW02tHTdauUqQ4nNyMiPbOUyRRUy+batedIiDV7/JclxE88DAAAAAJx4BAeMW04r0JIH7lDPgb2JzUhqRrY8t6hitj02BLMXaPCGP5bJeInnAQAAAADSQXDAuGSHoc596C4V9wwktoOxQcq6BRW9HkmHvj9Da9o87bzpTxXnetM7FgAAAACQQHDAuGNFkRb9+F5N2flmx70ZGWWcokpej6xhsSHsm6mdK29R3NOX1qkAAAAAgPdAcMD4YozOfuJHmr5ta8e5ERm578aGQ399o94+7Vx5i6LSjLQuBQAAAAAcBsEB44cxOuPphzTr1U0d52YUy3EK6vN6ZVuH/urGXo8GVt6qcOrctC4FAAAAALwPggPGjVN//qTmbv5lx82PYtl2QX3Z9thgMlkNrPiKgpmnpHQlAAAAAOBIEBwwLpy0/qc6+fmfdNyCKJbsXvV5BTm28+7jxnE0+Nkvy5+3MK0zAQAAAABHiOCAMTdn8/M67WePddyCKFZs9SRigyxLuz79BTVOPS+lKwEAAAAAR4PggDFV8H2d8fSDHbdWHCu28urLFeXabtu2e/nvqX7WRWmcCAAAAAA4Bu77/xLgxMg2GypVarKMSWxhbBSanPpyRWVGxIa9l/+Gqos/mdaZAAAAAIBjwGc4YEz0DWxTae9gxy0yRi2TVckrKutk2rb9F9+s8vmfTuNEAAAAAMBxIDggdYU9g1r84JqOn9kQGaMgzqjgleS52bZt6PxrdeCim9M6EwAAAABwHA77JRU/e+zj2ujm0roFk4AdBuobzOnH8Zf0qvmmYivWfwq/8vZqFBtbju3KCWwpOPS8sDRdwY4F0veGxuRuTEwnXTo5v2qM126MlmfffX0+aKv5piRLX6//btvj9X88kN5R6HqT9bUbALoRn+GA1NhRS6VdO2XFccf93dhgtf+1jApTFcw8NYULAQAAAACjheCAVFhRqNKunbKjsOMeG8nuFBt6SvJnnSbJSuFKAAAAAMBoITjghLPiSKXdO2WHrY77O7HBHREb4lyv/P4zJIvYAAAAAAATzWG/SO7xf/xS4v9xBo6W6zc6fhnFwK4tkjG69/5vyRoZFWxHkVfgExtwXP7tpf9rUn4jgyf+8YuyLWesz0AXyDRqbf88sGuLJGnNfavaHo96pqR1EiaB37n0n8b6BADAKDlscLAsR5GxJCV/mgBwJLJBU1ZsdPhyYGSG/cQKY9kK3IxMHBzmOcD7W//G92bog58f6zPSZ2UUdfgpMMDRyrzHa3c84u+XH/F6jdHzys4fLdIHf32szwAAjILDBocbfuf7isXXXeBYGHmDr8qp7u+4hibWXfetkkysz177x+8+HjuuyrPnKnYyaR2KLvbElhVvSZMvOHzmc/8kly9Fwijo2fqztn9ec98qxcbos9d+ue3xffNPT/MsdLlHX1yxSdcQHACgG9AScEJkd7/+nrEhMrHiOEz8/2bGtlWZNYfYAAAAAABdgOCAUZfZu01ueW/HLTKxojiUbbV/Twdj2SrPnKvIzaZxIgAAAADgBCM4YFRl9u9Q5sBgx+2d2GBZsdq+p4NlqTKzX1HWS+dIAAAAAMAJd9jv4fBb/6pPrUjKOnwtMN5f6bn7NP2R70o9yc2PApWbZWXdQM7bX1u+3tqiIG7p6uWbtG++n/K1mAyeeHCsLxgb/+JXSurL8lMqcPwW/MXftP3zeucNBVGgP3S/3vb445/6UppnoctN1tduAOhGfIYDRkVh448PxoYOgqilcrOijHUoNhxkaVepqH3zT0vnSAAAAABAag77GQ7Akeh5+WnNXPvXHbdWHKrsV+RaTbl2e9+qTJ2umtVK40QAAAAAQMr4DAccl/zrz2vWfX8pjfiZ7JIUxqGGmhXZaigzIja8esGlavYW0zoTAAAAAJAyggOOmbf9Jc2++1uyoiixRXGkIb8q29SVHREb3lx6od5a8pG0zgQAAAAAjAGCA45Jdvcb6r/za7JaQWJ7JzYorinrtP8V27HwA3r9/EvSOhMAAAAAMEYIDjhq7v6d6l99m2y/nthiE2soqCmOa/JGxIZdpy3SlouulCx+6gkAAAAAdDuCA46KU96jOatvk1MbSmyxiTXk1xRHVeVG/CjVvfNP10uXfIrYAAAAAACTBMEBR8yuD2nO6lVyh3YnNqNYZb+uMKokYsOBOfO16ZPXyThOWqcCAAAAAMYYwQFHxG7WNOeOP1Nm3/bEZmRU9utqdYgNlRn9euHKGxW7/ARWAAAAAJhMCA54X1bL1+y7vqns4GsdVqOKX1cQVpVzpOG5oT5lujYsW6Eok03rVAAAAADAOEFwwOFFoWbd+23l3trUca4EDflhRTnHtMWGZqFP669eqTDXk86dAAAAAIBxheCA9xbHmnn/f1HPqz/vOFdbDTVbydgQ5Hu1fvlKBb3FdO4EAAAAAIw7BAd0ZoxmPPS3Kmx+quNcazXVCCry7KgtNoReTuuXr1SzNDWdOwEAAAAA4xLBAR1NfeJ7Kj7/YMetEfqqBxVl7VD2sB9zGWUy2nDVCtWnzkzrTAAAAADAOEVwQELfM2s05Zm7Om7NMFDVLytrt+QMiw2x4+qFK25UZdbctM4EAAAAAIxjBAe0Kf5iraY9/r2Omx+1VPEryoyIDcaytemTn9GBuSendSYAAAAAYJwjOOBdvZue0IyH/77jFkQtlf2KMlZTrmW1bS9dslx7Tz4jjRMBAAAAABMEwQGSpJ4tz2rW/f9FMiaxteJQZb8iVw25dvtfmS0fvUK7zjg3rTMBAAAAABMEwQHKvblRs37wbSmOE1sYRyr7VdmmocyI2PD6hz6uHed8KK0zAQAAAAATCMFhksvu3KL+Nd+QFYaJLYojDflVKa4p67T/VXlr8Yf15nkfTetMAAAAAMAEQ3CYxDJ7t2nOnX8mq+Untsi8Exuq8kbEhoGzlujVCy6TRnwvBwAAAAAA3kFwmKTcoUHNuf022Y1qYotNrLJfUxzXErFh96ln6+WPLSM2AAAAAAAOi+AwCTnVfZpz+21yqvsTm3k7NoRRVTmnPSrsO2mBNl96rWTz1wYAAAAAcHh85DjJ2I2K+levkntgMLEZGZX9ulpRVfkRsaE8a55evPx6GcdJ61QAAAAAwARGcJhErKCp/jVfV3bPtsT2TmwIoopyI5pCdfosbVh2o2I3k9KlAAAAAICJjuAwSVhhoNl3fVPejlc6rEaVoKEgPBgbhn9uQ71vmjYsW6kom0vrVAAAAABAFyA4TAZxpFk/+M/Kv7mx41wNmvJbZeUc0xYb/N6iNly9Uq18Tzp3AgAAAAC6BsGh2xmjmT/6r+rZ8mzHudZqqNEqy3PittjQyvVo/fKb5RdK6dwJAAAAAOgqBIduZoymP/IPKrz4eMe53mqqHlTk2ZHsYbkhzHpaf/VNavRNS+tSAAAAAECXITh0salPfV+ln/+o49YIfdWCirJ2KNs6FBtix9XGq1aoNn12WmcCAAAAALoQwaFL9T17r6Y8fUfHrRkGqvoVZe2WnOGxwXb0whXXqzx7XlpnAgAAAAC6FMGhCxU2PKxpj/73jlsQtVTxK8pYQVtsMJalzZdeo/0nnZbWmQAAAACALkZw6DK9m5/SzHV/03ELolBlvyLXasq1rbbtlYuv1p4FC9M4EQAAAAAwCRAcukj+1Z9r5v3fkYxJbK34YGxw1FDGbv9j3/qRyzRw1pK0zgQAAAAATAIEhy7hbd+s2ff8uawoSmxRHKnsV2WbeiI2vPGBi7R98QVpnQkAAAAAmCQIDl0gO/ia+u/8uqywldiiONKQX5XimrJO+x/39kUf0hsfvDitMwEAAAAAkwjBYYLL7Nuu/jtWyfbriS02sYaCquK4Jm9EbBg8/Rxt/ejlkmUlngcAAAAAwPEiOExgTnm3+m+/TU69nNhiE2vIryqOaso57VFh78ln6OVLPkVsAAAAAACcMASHCcquD2nO6lVyK3sTm1Gssl9XGFUTseHAnJO16bLPyNj80QMAAAAAThw+6pyA7GZNc1avUmbfjsRmZFT262pFlURsqMycoxeuvEGx66Z1KgAAAABgkiI4TDBWy9fsNV9XdtfrHVajil9XEFaUc6ThuaE2dYY2XLVCUSab0qUAAAAAgMmM4DCRRKFm3/Pnym1/qeNcCRryO8SGRmmKNly9UmEun86dAAAAAIBJj+AwUcSxZv3wL5V/7Zcd52qroWaropxj2mJD0FPQhqtvVtBTSOVMAAAAAAAkgsPEYIxmPPA36n3pJx3nWqupRlCR58RtsaHl5bV++c1qFvvSuRMAAAAAgLcRHMY7YzTtsf+h4oZHOs6N0Fc9qMizw7Y/zCiT0YZlK1SfMj2dOwEAAAAAGIbgMM5N+ckd6nv23o5bMwxU9cvK2i3Z1qHPbYgdVxuvvFHVmXPSOhMAAAAAgDYEh3Gs9NwPNfXJ73fc/ChQxa8oa7fkDIsNxrb14uXXaWjOyWmdCQAAAABAAsFhnCq88KimP/LfOm5B1FLZryhj+e2xwbK0+ZJrtG/+6WmdCQAAAABARwSHcahny081c+1fd9xacaiyX5ErX65ttW1bPnqFdp++KI0TAQAAAAA4LILDOJN7c4Nm3fttKY4TWxiHGmpWZJuGMiNiw2vnX6Kdiz6Y1pkAAAAAABwWwWEc8Xa8rP4135QVRYktiiMN+VVZpq6s0/7Htm3JR7Rt6YVpnQkAAAAAwPsiOIwT2d1vqP/Or8lq+YktMgdjg+KavBGxYefZ5+m1D39CsqzE8wAAAAAAGCsEh3HA3b9T/atXyW7WEltsYpX9muIOsWH3goV65aIriQ0AAAAAgHGH4DDGnMpezVl9m5zagcRmTKwhv6YoqirntEeFfSct0OZPXCPZ/BECAAAAAMYfPlodQ3a9rDmrV8kd2p3YjGIN+XWFHWLDUP9JevHy62UcJ61TAQAAAAA4KgSHMWL5dfXf+TVl9r6V2IyMyn5draii3IimUJ0+WxuvvEGxm0npUgAAAAAAjh7BYQxYYaD+u74pb2Brh9Wo4tcVhFXlHGn45zbU+6Zpw7KbFGVzaZ0KAAAAAMAxITikLQo1695vK7ftxY5zJWjIDyvKOaYtNjQLfVq/fKVa+Z507gQAAAAA4DgQHNJkjGb+6K/Us/W5jnOt1VCzVZHnxG2xIcj3asPVKxX0ltK5EwAAAACA40RwSIsxmv7Q36mw6cmOc73VVD2oyLMj2cNyQ5j1tOHqlWr0TU3rUgAAAAAAjhvBISVTn/xfKv1yXcetEfqqBRVl7VC2dSg2RG5GG5atUG3azLTOBAAAAABgVBAcUtD307s15Sd3dtyaYaCqX1HWbskZFhti29GLV3xWlVnz0joTAAAAAIBRQ3A4wYrPP6Bpj/3PjpsftVTxK8rYQVtsMJalzZd9WvvnnZrSlQAAAAAAjC6CwwnUu+lJzXjwbztuQRQejA1WU65ltW0vf3y59px6VhonAgAAAABwQhAcTpD8q89p1v3fkYxJbK04VNmvyFFDrt3+R7D1I5/U4JmL0zoTAAAAAIATguBwAuS2vaDZ9/yFFMeJLYwjlf2qbFNXZkRseOODF2v74g+ndSYAAAAAACcMwWGUZQe2avaab8gKW4ktejs2KK4p67S/6befe77e+MBFaZ0JAAAAAMAJRXAYRZm9b6n/jj+THTQTW2xiDQVVxXFV3ojYMHjmYm39yCelEd/LAQAAAACAiYrgMErcoV2as/o2OY1KYotNrCG/qjiqKTciNuw59Sy9fPHVxAYAAAAAQFchOIwCp7pf/bffJqeyL7EZxSr7NYVRVTmnPSrsn3uKNl16rYzNHwMAAAAAoLvwke5xsptV9d+xSpkDA4nNyKjcrKsVVZUfERvKs+bqhSs+K+O4aZ0KAAAAAEBqCA7HwQqa6r/za8rufjOxGRlV/LqCqKqc077Vps3UxqtWKM5kU7oUAAAAAIB0ERyOVdjS7Lv/L3k7XukwGlWChvywopxjNPxzGxqlKVq/7CaFXi6tSwEAAAAASB3B4VjEkWbd95fKv7G+41wNmvJbydjg9xa1fvnNavUU0rkTAAAAAIAxQnA4WsZo5rr/R72vPNNxrrWaarQq8py4LTa0cnmtv3ql/EJfOncCAAAAADCGCA5HwxhN+/F3Vdj4aMe5HvqqB2V5dtj2hg2znjYsW6nGlOmpnAkAAAAAwFgjOByFKU/frr7nfthxa4SBan5ZWTuUbR363IbYcfXClTeoOmN2WmcCAAAAADDmCA5HqPTcfZr61O0dNz8MVPXLytotOcNig7FtvXj59Rrqn5/WmQAAAAAAjAsEhyNQ2PhjTX/kux23IGqp7FeUsYL22GBZ2vyJa7Vv/mlpnQkAAAAAwLhBcHgfPS8/rZlr/7rj1opDlf2KXKsp17batlc+dpV2n7YwjRMBAAAAABh3CA6HkX/9ec267y8lYxJbGIcaalZkq6GM3f5mfPWCSzVw9nlpnQkAAAAAwLhDcHgP3vaXNPvub8mKosQWxZGG/KpsU1d2RGx4c+mFemvJR9I6EwAAAACAcYng0EF29xvqv/NrslpBYnsnNiiuKeu0v/l2LPyAXj//krTOBAAAAABg3CI4jODu36H+1bfJ9uuJLTaxhoKa4rgmb0Rs2HXaIm256ErJshLPAwAAAABgsiE4DOOU92jO6lVyakOJLTaxhvya4qiqnNMeFfbOP10vXfIpYgMAAAAAAG8jOLzNrg9pzupVcod2JzajWGW/rjCqJGLDUP98bfrkdTKOk9apAAAAAACMewQHSXazpjl3/Jky+7YnNiOjsl9Xq0NsqMzo18Yrb1DsummdCgAAAADAhDDpg4PV8jX7rm8qO/hah9Wo4tcVhFXlHGl4bqhPma4Ny1YoynppnQoAAAAAwIQxuYNDFGr2PX+h3FubOs6VoCE/rCjnmLbY0Cz0af3VKxXmetK5EwAAAACACWbyBoc41qz7v6P8a7/oOFdbDTVbydgQ5Hu1fvlKBb3FdO4EAAAAAGACmpzBwRjNePD/U+/mf+4411pNNYKKPDtqiw2hl9P65SvVLE1N504AAAAAACaoyRccjNG0x/6niusf6jg3Ql/1oCLPDmUP+zGXUSajDctuUn3qzLQuBQAAAABgwpp0waHvmTXqe/aejlszDFT1y8rarbbYEDuuXrjiRlVmzknrTAAAAAAAJrRJFRyKv1iraU/8U8fNj1qq+BVl7JacYbHBWLY2ffIzOjD35LTOBAAAAABgwps0waHw4uOa8fDfd9yCqKWyX1HGasq1rLbtpUuWa+/JZ6RxIgAAAAAAXWNSBIeeLc9q5o/+SjImsbXiUGW/IlcNuXb7m2PLR6/QrjPOTetMAAAAAAC6RtcHh9ybGzXrB9+W4jixhXGkcrMq2zSUGREbXjv/Eu0450NpnQkAAAAAQFfp6uCQ3blF/Wu+ISsME1sURxryq5KpKeu0vxneWnyBti29MK0zAQAAAADoOl0bHDJ73tScO/9MVstPbJF5OzbEVXkjYsPAWUv06gWXSiO+lwMAAAAAADhyXRkc3AMDmrN6lexGNbHFJlbZrymOa4nYsHvB2Xr5Y8uIDQAAAAAAHKeuCw5OdZ/mrF4lp7o/sZm3Y0MYVZVz2qPCvnkLtPkT10p2171JAAAAAABIXVd9dG03KupfvUrugcHEZmRU9utqRVXlR8SG8qx5evGK62UcJ61TAQAAAADoal0THKygof47v6bsnm2J7Z3YEEQV5UY0her0Wdqw7EbFbialSwEAAAAA6H5dERysMNDsu74lb+eWDqtRxW8oCA/GhuGf21Dvm6YNy1YqyubSOhUAAAAAgElh4geHKNSsH/xn5d/c2HGuBk35YVk5x7TFBr+3qA1Xr1Qr35POnQAAAAAATCITOzgYo5lr/1o9W57tONdaDTVaZXlO3BYbWrkerV9+s/xCKZ07AQAAAACYZCZucDBG0x/+exVefLzjXG81VQ8q8uxI9rDcEGY9rb/6JjX6pqV1KQAAAAAAk86EDQ5Tn/q+Sr9Y23FrhL5qQUVZO5RtHYoNkZvRxqtWqDZ9dlpnAgAAAAAwKU3I4ND37D2a8vQdHbdmGKjqV5S1W3KGxYbYdvTi5derPHteWmcCAAAAADBpTbjgUFz/kKY9+j86bkHUUsWvKGMFbbHBWJY2X3qt9p+0IK0zAQAAAACY1CZUcOjd/JRmPPD/dtyCKFTZr8i1mnJtq2175eKrtWfB2WmcCAAAAAAANIGCQ/7Vn2vm/d+RjElsrfhgbHDUUMZu/6+09SOXaeCsJWmdCQAAAAAANEGCg/fWJs2+589lRVFii+JIZb8q29QTseGND1yk7YsvSOtMAAAAAADwtnEfHLKDr6l/zTdkha3EFsWRhvyqFNeUddr/q2w/50N644MXp3UmAAAAAAAYZlwHh8y+7eq/Y5Vsv57YYhNrKKgqjqvyRsSGwTPO1dYLL5csK/E8AAAAAABw4o3b4OCUd6v/9tvk1MuJLTaxhvyq4qim3IjYsOeUM/Xyx5cTGwAAAAAAGEPjMjg4tQOas3qV3MrexGYUq+zXFUZV5Zz2qLB/7inafOmnZexx+V8LAAAAAIBJY9x9ZG43a+q/Y5Uy+3YkNiOjsl9XK6okYkN55ly9cMVnFbtuWqcCAAAAAID3MK6CgxU0NXvN15Xd9UaH1aji1xWEFeUcaXhuqE2doY1X3ag4k03rVAAAAAAAcBjjJziELc2+58+V2/5Sx7kSNOR3iA2N4hRtuHqlwlw+nTsBAAAAAMD7Gh/BIY416/7/W/nXn+84V1sNNVsV5RzTFhuCnoI2LL9ZQU8hnTsBAAAAAMARGfvgYIxmPPA36n3pJx3nWqupRlCR58RtsaHl5bV++c1qFvvSuRMAAAAAAByxsQ0Oxmjao/9dxQ2PdJwboa96UJFnh22HRpmsNixbofqU6encCQAAAAAAjsqYBocpP7lDfT/7QcetGQaq+mVl7ZZs69DnNsSOq41X3qDqzDlpnQkAAAAAAI7SmAWH0s/v19Qnv99x86NAFb+irN2SMyw2GNvWi5dfp6E5J6d1JgAAAAAAOAZjEhwKLzym6Q//Q8ctiFoq+xVlLL89NliWNn/iGu2bf3paZwIAAAAAgGOUenDo2fJTzVz7XzturThU2a/IlS/Xttq2LRddqd2nLUrjRAAAAAAAcJxSDQ65N9Zr1g/+kxTHiS2MQw01K7JNQ5kRseG1D39COxd+IKUrAQAAAADA8UotOHg7Xlb/Xd+SFYaJLYojDflVWaaurNN+0ralF2rb0gvTOhMAAAAAAIyCVIJDZvcb6r/za7JafmJ7JzYorskbERt2nn2eXjv/kjROBAAAAAAAo+iEBwd3/07NWb1KdrOW2GITqxzUFHeIDbtOW6hXLrpSsqzE8wAAAAAAwPh2QoODU9mrOatXyakdSGyxiTXk1xRFVeWc9qiw76TT9NIl10j2mP3UTgAAAAAAcBxO2Ef0dr2sOatXyR3aldiMYpX9usIOsWGof75evPw6Gcc5UacBAAAAAIAT7IQEB8uvq//Orymz963EZmRU9utqRRXlRjSF6vTZ2njlDYrdzIk4CwAAAAAApGTUg4MVBuq/65vyBrZ2WI0qfl1BWFXOkYZ/bkO9b5o2LLtJUdYb7ZMAAAAAAEDKRjc4RKFm3ftt5ba92HGuBA35YUU5x7TFhmahT+uXr1Qr3zOq5wAAAAAAgLExesHBGM380V+pZ+tzHedaq6FmqyLPidtiQ5Dv1frlKxX0lkbtFAAAAAAAMLZGJzgYo+kP/a0Km57sONdbTdWDijw7kj0sN4RZTxuuXqlmaeqonAEAAAAAAMaHUQkOU5/4nkq/fKDj1gh91YKKsnYo2zoUGyI3ow3LVqg2beZonAAAAAAAAMaR4w4Ofc+s0ZRn7uq4NcNAVb+irN2SMyw2xLajF674rCqz5h3vbw8AAAAAAMah4woOxecf0LTHv9dx86OWKn5FGTtoiw3GsrTpsk/rwLxTj+e3BgAAAAAA49gxB4feTU9qxoN/23EL3okNVlOuZbVtL398ufaeetax/rYAAAAAAGACOKbgkN/6M826/zuSMYmtFYcq+1U5asi12//1Wy+8XINnLj62SwEAAAAAwIRx1MEht+0Fzb7321IcJ7YwjlT2q7JNXZkRseH1D12s7eeef+yXAgAAAACACeOogkN2YKtmr/mGrLCV2KK3Y4PimrJO+7/2rXM/rDfPu+j4LgUAAAAAABPGEQeHzN5tmnPHKtlBM7FFJtKQX1UcV+WNiA0DZy7Wqx+5TBrxvRwAAAAAAED3OqLg4A7t0pzVq2Q3qoktNrHKfk1xXFNuRGzYc+pZeuXiq4kNAAAAAABMMu8bHJzqPvXffpucyr7EZnQwNoRRVTmnPSrsn3uKNl16rYx9XD95EwAAAAAATEDu4Ua7UdHcNV9T5sBAYjMyKjfrakVV5UfEhvKseXrhis/KOIf91wMAAAAAgC512CLQv+Ybyu5+I/ElEUZGFb+uIKoq57Q/pzZtpjZedaPiTHbUjwUAAAC6zfTXp7/7n/OtvGyr/TOEtz+9Pe2Txp3ZraHEY5Gxla1l2h77yQ8ziV83mQz50RH9up94zvv/oi633jsz8VjgtGRb7T+NcR/vfzqrfFbisdjE6in3vO9zD/v1DrkdL3d41KgSNOSHFeUco+EpolGaovVXr1To5d73NwYAAAAAAN3rqL/BQjVoym8lY4PfW9T65Terle8dxfMAAAAAAMBEdFTBodZqqNGqyHOittjQyuW1/uqV8gt9o3weAAAAAACYiI74uzrWQ1/1oCLPDmUPyw1h1tOGZSvVmDL9MM8GAAAA0Mkd/3DHu/85iiNZI75/Ws7ly5UDP/lhi5ES3+8in7MSv24yacXmiH5dxp7cbydJcurJ7wsSm1gj3zJB/v2/T0G3a4bNxGPGGDn2we8F8q0/+NZ7PveIgkMjDFTzy8raoexhL4Cx42rjlTeqOmP20d4MAAAAoANj2j9ojE38Hr9y8jDq9IG0lXj0yD7cBm+n934bjPy7xvtf8jXpaLxvcPCjlqp+WVm7JWdYbDC2rReuuF7l/pOO+TcHAAAAJrubfuumw+7nzTkvpUvGr02v+0f06xad6p3gS8a3HbXWEf26ub2T+6d5SJI2//TIft3Cj5zYOyaA53c+f8zPPWxwCONQVb8iWw3JWIreDhvGsrT+Y1dqcNZsKagd82+OyY1aCAAAAADdyzrcp0f82j1/9VtbBtZN21PZVE7xJkwuSyXtH+sj0LV2b/3i1u+M9RFp+9Yv9nxe0syxvgPd6ZWdP1r06Iv/cdNY34GuNilfuwGgGx02OAAAAAAAAByLo/qxmAAAAAAAAEeC4AAAAAAAAEYdwQEAAAAAAIw6ggMAAAAAABh1BAcAAAAAADDq/n+ljEA82CWoigAAAABJRU5ErkJggg==\n",
      "text/plain": [
       "<Figure size 1339.2x432 with 3 Axes>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "fig, axs = plt.subplots(1,3, figsize = (18.6,6) )\n",
    "\n",
    "# Data\n",
    "x = np.linspace(-5, 5, 100)\n",
    "\n",
    "axis_width = 10\n",
    "alpha = 0.25\n",
    "line_alpha = 0.7\n",
    "fill_alpha = 0.4\n",
    "\n",
    "\n",
    "ax = axs[0]\n",
    "ax.plot(x, x, color='orangered', linewidth=axis_width,alpha = line_alpha)\n",
    "ax.axhline(0, color='royalblue', linewidth=axis_width, alpha=line_alpha)\n",
    "# ax.plot(x+x, x, color='skyblue', linewidth=axis_width,alpha = line_alpha)\n",
    "# ax.plot(-x-x, x, color='forestgreen', linewidth=axis_width,alpha = line_alpha)\n",
    "ax.axhline(0, color='black', linewidth=2, alpha=0.7, linestyle='-')\n",
    "ax.axvline(0, color='black', linewidth=2, alpha=0.7, linestyle='-')\n",
    "\n",
    "# fill between\n",
    "ax.fill_between(x=[-5,0], y1=0, y2=5, color='forestgreen', alpha = fill_alpha)\n",
    "ax.fill_between(x=[0,5], y1=0, y2=-5, color='forestgreen', alpha = fill_alpha)\n",
    "ax.fill_between(x=[-5,0], y1=-5, y2=[-5,0], color='forestgreen', alpha = fill_alpha)\n",
    "ax.fill_between(x=[0,5], y1=5, y2=[0,5], color='forestgreen', alpha = fill_alpha)\n",
    "ax.fill_between(x=[-5,0], y1=0, y2=[-5,0], color='skyblue', alpha = fill_alpha)\n",
    "ax.fill_between(x=[0,5], y1=0, y2=[0,5], color='skyblue', alpha = fill_alpha)\n",
    "\n",
    "# Set limits and turn off spines\n",
    "ax.set_xlim(-5, 5)\n",
    "ax.set_ylim(-5, 5)\n",
    "ax.spines['top'].set_visible(False)\n",
    "ax.spines['right'].set_visible(False)\n",
    "ax.spines['bottom'].set_visible(False)\n",
    "ax.spines['left'].set_visible(False)\n",
    "ax.tick_params(axis='both', which='both', bottom=False, top=False, left=False, right=False, labelbottom=False, labelleft=False)\n",
    "\n",
    "\n",
    "ax = axs[1]\n",
    "ax.axvline(0, color='orangered', linewidth=axis_width, alpha=line_alpha)\n",
    "ax.axhline(0, color='royalblue', linewidth=axis_width, alpha = line_alpha)\n",
    "# ax.plot(x, x, color='skyblue', linewidth=axis_width,alpha = line_alpha)\n",
    "# ax.plot(-x, x, color='forestgreen', linewidth=axis_width,alpha = line_alpha)\n",
    "ax.axhline(0, color='black', linewidth=2, alpha=0.7, linestyle='-')\n",
    "ax.axvline(0, color='black', linewidth=2, alpha=0.7, linestyle='-')\n",
    "\n",
    "# fill between\n",
    "ax.fill_between(x=[-5,0], y1=0, y2=5, color='forestgreen', alpha = fill_alpha)\n",
    "ax.fill_between(x=[0,5], y1=0, y2=-5, color='forestgreen', alpha = fill_alpha)\n",
    "ax.fill_between(x=[-5,0], y1=0, y2=-5, color='skyblue', alpha = fill_alpha)\n",
    "ax.fill_between(x=[0,5], y1=0, y2=5, color='skyblue', alpha = fill_alpha)\n",
    "\n",
    "\n",
    "# Set limits and turn off spines\n",
    "ax.set_xlim(-5, 5)\n",
    "ax.set_ylim(-5, 5)\n",
    "ax.spines['top'].set_visible(False)\n",
    "ax.spines['right'].set_visible(False)\n",
    "ax.spines['bottom'].set_visible(False)\n",
    "ax.spines['left'].set_visible(False)\n",
    "ax.tick_params(axis='both', which='both', bottom=False, top=False, left=False, right=False, labelbottom=False, labelleft=False)\n",
    "\n",
    "\n",
    "\n",
    "ax = axs[2]\n",
    "ax.axvline(0.5, color='orangered', linewidth=axis_width, alpha=line_alpha)\n",
    "ax.axvline(-0.5, color='royalblue', linewidth=axis_width, alpha=line_alpha)\n",
    "ax.axvline(0.0, color='skyblue', linewidth=axis_width, alpha=line_alpha)\n",
    "ax.axvline(-1, color='forestgreen', linewidth=axis_width, alpha=line_alpha)\n",
    "ax.axvline(1, color='forestgreen', linewidth=axis_width, alpha=line_alpha)\n",
    "ax.axhline(0, color='gray', linewidth=axis_width, alpha = line_alpha)\n",
    "ax.axhline(0, color='black', linewidth=2, alpha=0.7, linestyle='-')\n",
    "\n",
    "# fill between\n",
    "ax.fill_between(x=[-1,-0.5], y1=0, y2=1.0, color='forestgreen', alpha = fill_alpha)\n",
    "ax.fill_between(x=[0.5,1], y1=0, y2=1.0, color='forestgreen', alpha = fill_alpha)\n",
    "ax.fill_between(x=[-0.5,0.5], y1=0, y2=1.0, color='skyblue', alpha = fill_alpha)\n",
    "\n",
    "# Set limits and turn off spines\n",
    "ax.set_xlim(-1.02, 1.02)\n",
    "ax.set_ylim(-0.01, 1)\n",
    "ax.spines['top'].set_visible(False)\n",
    "ax.spines['right'].set_visible(False)\n",
    "ax.spines['bottom'].set_visible(False)\n",
    "ax.spines['left'].set_visible(False)\n",
    "ax.tick_params(axis='both', which='both', bottom=False, top=False, left=False, right=False, labelbottom=False, labelleft=False)\n",
    "\n",
    "plt.subplots_adjust(wspace=0.2)\n",
    "\n",
    "plt.savefig('./figs/fig1C',dpi=450,bbox_inches='tight')\n",
    "plt.show();"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "id": "2f5e8a7c",
   "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.9.7"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 5
}