snews_feedbacks.ipynb

{
 "cells": [
  {
   "cell_type": "markdown",
   "id": "410f81c1",
   "metadata": {},
   "source": [
    "# Heartbeats and Feedbacks on the Heartbeats\n",
    "\n",
    "20.01.2023"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 24,
   "id": "e91e743f",
   "metadata": {},
   "outputs": [],
   "source": [
    "import pandas as pd\n",
    "import numpy as np\n",
    "from datetime import datetime, timedelta\n",
    "import time\n",
    "import matplotlib.pyplot as plt\n",
    "import os, json\n",
    "from snews_pt.snews_pub import SNEWSTiersPublisher"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "8aa7778b",
   "metadata": {},
   "source": [
    "On the server both the coincidence tracking script and the feedback generation script should be running.<br>\n",
    "\n",
    "```bash\n",
    "snews_cs run-coincidence\n",
    "snews_cs run-feedback\n",
    "```"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "3b7350fb",
   "metadata": {},
   "source": [
    "Now we invoke 7 heartbeat messages, sleep 5 sec in between and wait. <br>\n",
    "There are two main functionalities we can investigate. Once the server acquires enough data, it starts measuring the heartbeat frequency. Every minute, server checks the registered heartbeats from the last 24 hours for each detector. <br>\n",
    "\n",
    "As soon as there is an inconsistency e.g.<br>\n",
    "Detector was sending beats every 10minutes+/-1min, and at the time of next check, if it has been more than 10+3*1 minutes, it issues a warning message to the users listed for that detector."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 25,
   "id": "d14ab115",
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Message Generated for Heartbeat\n",
      "\u001b[94m----------------------------------------------------------------\u001b[0m\n",
      "\u001b[91mSending message to Heartbeat on kafka://kafka.scimma.org/snews.experiments-firedrill\u001b[0m\n",
      "_id                :19_Heartbeat_2023-01-20T11:59:35.646130\n",
      "detector_name      :XENONnT\n",
      "machine_time       :2023-01-20T11:59:35.646130\n",
      "detector_status    :ON\n",
      "meta               :\n",
      "schema_version     :1.2.0\n",
      "sent_time          :2023-01-20T11:59:35.646130\n"
     ]
    }
   ],
   "source": [
    "msg = SNEWSTiersPublisher(detector_name='XENONnT', detector_status='ON',) \n",
    "msg.send_to_snews()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 26,
   "id": "783fc3c0",
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Message Generated for Heartbeat\n",
      "\u001b[94m----------------------------------------------------------------\u001b[0m\n",
      "\u001b[91mSending message to Heartbeat on kafka://kafka.scimma.org/snews.experiments-firedrill\u001b[0m\n",
      "_id                :19_Heartbeat_2023-01-20T12:02:56.624173\n",
      "detector_name      :XENONnT\n",
      "machine_time       :2023-01-20T12:02:56.624173\n",
      "detector_status    :ON\n",
      "meta               :\n",
      "schema_version     :1.2.0\n",
      "sent_time          :2023-01-20T12:02:56.624173\n",
      "Message Generated for Heartbeat\n",
      "\u001b[94m----------------------------------------------------------------\u001b[0m\n",
      "\u001b[91mSending message to Heartbeat on kafka://kafka.scimma.org/snews.experiments-firedrill\u001b[0m\n",
      "_id                :19_Heartbeat_2023-01-20T12:03:05.124962\n",
      "detector_name      :XENONnT\n",
      "machine_time       :2023-01-20T12:03:05.124962\n",
      "detector_status    :ON\n",
      "meta               :\n",
      "schema_version     :1.2.0\n",
      "sent_time          :2023-01-20T12:03:05.124962\n",
      "Message Generated for Heartbeat\n",
      "\u001b[94m----------------------------------------------------------------\u001b[0m\n",
      "\u001b[91mSending message to Heartbeat on kafka://kafka.scimma.org/snews.experiments-firedrill\u001b[0m\n",
      "_id                :19_Heartbeat_2023-01-20T12:03:13.604559\n",
      "detector_name      :XENONnT\n",
      "machine_time       :2023-01-20T12:03:13.604559\n",
      "detector_status    :ON\n",
      "meta               :\n",
      "schema_version     :1.2.0\n",
      "sent_time          :2023-01-20T12:03:13.604559\n",
      "Message Generated for Heartbeat\n",
      "\u001b[94m----------------------------------------------------------------\u001b[0m\n",
      "\u001b[91mSending message to Heartbeat on kafka://kafka.scimma.org/snews.experiments-firedrill\u001b[0m\n",
      "_id                :19_Heartbeat_2023-01-20T12:03:22.076504\n",
      "detector_name      :XENONnT\n",
      "machine_time       :2023-01-20T12:03:22.076504\n",
      "detector_status    :ON\n",
      "meta               :\n",
      "schema_version     :1.2.0\n",
      "sent_time          :2023-01-20T12:03:22.076504\n",
      "Message Generated for Heartbeat\n",
      "\u001b[94m----------------------------------------------------------------\u001b[0m\n",
      "\u001b[91mSending message to Heartbeat on kafka://kafka.scimma.org/snews.experiments-firedrill\u001b[0m\n",
      "_id                :19_Heartbeat_2023-01-20T12:03:30.568242\n",
      "detector_name      :XENONnT\n",
      "machine_time       :2023-01-20T12:03:30.568242\n",
      "detector_status    :ON\n",
      "meta               :\n",
      "schema_version     :1.2.0\n",
      "sent_time          :2023-01-20T12:03:30.568242\n",
      "Message Generated for Heartbeat\n",
      "\u001b[94m----------------------------------------------------------------\u001b[0m\n",
      "\u001b[91mSending message to Heartbeat on kafka://kafka.scimma.org/snews.experiments-firedrill\u001b[0m\n",
      "_id                :19_Heartbeat_2023-01-20T12:03:39.046397\n",
      "detector_name      :XENONnT\n",
      "machine_time       :2023-01-20T12:03:39.046397\n",
      "detector_status    :ON\n",
      "meta               :\n",
      "schema_version     :1.2.0\n",
      "sent_time          :2023-01-20T12:03:39.046397\n",
      "Message Generated for Heartbeat\n",
      "\u001b[94m----------------------------------------------------------------\u001b[0m\n",
      "\u001b[91mSending message to Heartbeat on kafka://kafka.scimma.org/snews.experiments-firedrill\u001b[0m\n",
      "_id                :19_Heartbeat_2023-01-20T12:03:47.494588\n",
      "detector_name      :XENONnT\n",
      "machine_time       :2023-01-20T12:03:47.494588\n",
      "detector_status    :ON\n",
      "meta               :\n",
      "schema_version     :1.2.0\n",
      "sent_time          :2023-01-20T12:03:47.494588\n"
     ]
    }
   ],
   "source": [
    "for i in range(7):\n",
    "    time.sleep(5)\n",
    "    msg = SNEWSTiersPublisher(detector_name='XENONnT', detector_status='ON',) \n",
    "    msg.send_to_snews()"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "51fc24aa",
   "metadata": {},
   "source": [
    "<!-- ![](explosion.png) -->"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "5c4acb2b",
   "metadata": {},
   "source": [
    "Now that we send beats every 5 seconds, the server should realize that after 7th message there were no messages and sends a warning. After sending this warning, it resets itself until there are messages again."
   ]
  },
  {
   "cell_type": "markdown",
   "id": "9e9874a0",
   "metadata": {},
   "source": [
    "This is how it looks like in the server;\n",
    "```bash\n",
    "[DEBUG] >>>>> Checking if beat skipped                                                                                  \n",
    "[DEBUG] >>>>> mean:5.61, std:3.46, trigger at 307.00                                                                  \n",
    "[DEBUG] >>>>> Delay since last: 23.98                                                                                   \n",
    "[DEBUG] >>>>> Running minute: 5 \n",
    "```\n",
    "Then when it notices a skipped beat after 307 it issues the following\n",
    "```bash\n",
    "[DEBUG] >>>>>                                                                                                             \n",
    "Your (XENONnT) heartbeat frequency is every 8.07+/-5.61 sec                                                            \n",
    "Expected a heartbeat at 2023-01-20T11:59:47.754112 +/- 5.61 sec                                                         \n",
    "Since last heartbeat there has been 35.48 sec                                                                           \n",
    "Is everything alright? Do you wanna talk about it?                                                                                                                                                                                             'smtps://user:password@smtp.gmail.com' mutt   \n",
    "-F /dev/null   \n",
    "-e 'set from=snews_heartbeats@snews.org'   \n",
    "-e 'set smtp_url=$SMTP_URL'   \n",
    "-s 'SNEWS Server Heartbeat for XENONnT is skipped!'   kara@kit.edu << EOM  \n",
    "Dear XENONnT,   Your (XENONnT) heartbeat frequency is every 8.07+/-5.61 sec                                                                             \n",
    "Expected a heartbeat at 2023-01-20T11:59:47.754112 +/- 5.61 sec                                                         \n",
    "Since last heartbeat there has been 35.48 sec                                                                           \n",
    "Is everything alright? Do you wanna talk about it? \n",
    "Provided by the SNEWS server. \n",
    "--Cheers EOM\n",
    "```"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "cfc3ba55",
   "metadata": {},
   "source": [
    "## Feedback Requests\n",
    "\n",
    "Now let's have a look at the second functionality. <br>\n",
    "\n",
    "In this case, the user should be allowed to \"request\" a feedback on their heartbeats within last 24 hours. Once the server receives such a request, it compiles a png file with some information, (and a csv?) and sends it to the requesting user via email."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 28,
   "id": "9046b26f",
   "metadata": {},
   "outputs": [],
   "source": [
    "from snews_pt.remote_commands import get_feedback"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 29,
   "id": "66d185e0",
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Heartbeat Feedback is requested! Expect an email from us!\u001b[0m\n"
     ]
    }
   ],
   "source": [
    "get_feedback(detector_name='XENONnT', email_address='kara@kit.edu')"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "9e993a42",
   "metadata": {},
   "source": [
    "Which appears on the server as;\n",
    "\n",
    "```bash\n",
    "[DEBUG] >>>>>                                                                                                    \n",
    "{'_id': '0_Get-Feedback', 'email': 'kara@kit.edu', 'detector_name': 'XENONnT', 'meta': {}}                                                                                                                                                     \n",
    "> [DEBUG] These mails are passed kara@kit.edu for detector: XENONnT \n",
    "'smtps://$USER:$PASSWORD@smtp.gmail.com' mutt   \n",
    "-F /dev/null   -e 'set from=snews_heartbeats@snews.org'   \n",
    "-e 'set smtp_url=$SMTP_URL'   -s 'SNEWS Server Feedback for %detector'   \n",
    "-a /snews_cs/beats/XENONnT_2023-01-20_12H08M.png \n",
    "--    kara@kit.edu << EOM  \n",
    "\n",
    "Dear XENONnT, \n",
    "\n",
    "Attached /snews_cs/beats/XENONnT_2023-01-20_12H08M.png please find the feedback information,  Your heartbeat frequency is every 46.44+/-378.659624 sec                                    \n",
    "The last heartbeat received at 2023-01-20 12:03:51.407667                                                               \n",
    "The received heartbeat frequency, together with the computed latency is plotted, and sent in the attachment.            \n",
    "Provided by the SNEWS server. --Cheers EOM \n",
    "```"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 30,
   "id": "be59855b",
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Heartbeat Feedback is requested! Expect an email from us!\u001b[0m\n"
     ]
    }
   ],
   "source": [
    "get_feedback(detector_name='XENONnT', email_address='kara@kit.edu;test@hotmail.com')"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "ea2f8390",
   "metadata": {},
   "source": [
    "You can also pass multiple emails separated with a semicolon, if some of these email addresses are not registered for the given detector, the server will ignore it \n",
    "\n",
    "```bash\n",
    "> [DEBUG] These mails are passed kara@kit.edu; test@hotmail.com for detector: XENONnT\n",
    "> [DEBUG] The given email: test@hotmail.com is not registered for XENONnT!\n",
    "```"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "b06245b3",
   "metadata": {},
   "source": [
    "Below is an example how the server looks the data within last 24 hours to create the feedback."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 33,
   "id": "47f17526",
   "metadata": {},
   "outputs": [],
   "source": [
    "from snews_cs.heartbeat_feedbacks import plot_beats"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 32,
   "id": "d9a2aba9",
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAwAAAADgCAYAAAC0LQGRAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjUuMSwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/YYfK9AAAACXBIWXMAAAsTAAALEwEAmpwYAAAvwUlEQVR4nO3dd5hdVb3G8e+bJmFSIBBCCV0E8RoCBJCiRqpIl6ISqkgR0Chw1VhgAiIKRhFBLwHpKBAUDUU6CeCNQBoQyqVHEhJCCaQSUn73j7UnOZmcc2bPzJk5M5n38zznmb3XXu0Uwl57NUUEZmZmZmbWMXSqdgXMzMzMzKz1uAFgZmZmZtaBuAFgZmZmZtaBuAFgZmZmZtaBuAFgZmZmZtaBuAFgZmZmZtaBuAFgZqsNSbWSbmqlsiTpWkmzJT3ZGmV2BJKuk/TzatejnPZQRzOzctwAMLNWIekNSXvXCztB0uP14iyUNC+7sb5b0sYtVJ/m3sTtAewD9I+InStQn+0lzZH0yYKwHSV9IGmz7Lzw86l7XZ5dO0FSSPpBvXynSRpccL6tpNGSPpQ0V9IjknYruL5Zls899fK5SVJtjvfxz4K6LZb0ccH5/zTt02mc+r+rEnF+Lenl7DN4UdJxJeIdl30e32qZ2pqZtT43AMysrTkoInoAGwBvA7+vcn1K2RR4IyLmNzahpC71wyJiEnA5cFXWu9AVuAY4NyLeKIh6UET0KHidWXDtfeAHknqWKHdL4F/As8DmwIbAHcD9knatF32XwoZBXhGxf13dgJuBiwvqelpj82tB84GDgN7A8cDv6r9fSWsDPwaea2ohxb5rM7NqcwPAzNqkiPgIuB3YtlQcSZtLGps9xX0AWLfe9VGSZmZPux+V9Jks/BRgCOlmeZ6kO7PwH0l6NcvveUmHlSj3JOBqYNcs/fAs/GRJr0h6P3vKvmFBmpB0hqSXgZdLvKXhpIbPKaQbz3mkRkFeLwDjgLNKXK8FxkXETyLi/YiYGxGXATcCv6oX92LgwmKZSBqc9SycLWmWpBmSTmxEPXORtLakuyS9k/UI3SWpf8H1EyS9ln1fr0saIunTwP+w4rv5oFjeEXFeRLwYEcsi4gngMaB+I+gi4DLg3SJZrJ31UM2V9ETWuKqrV57v2sysatwAMLM2SdKawNeAf5eJ9mdgAunG/wLSk9xC/wS2AtYDJpKeSBMRI1n56fRBWfxXgc+TngoPB26StEH9QiPiT8BppJvpHhFxnqQ9STeMR5Fu4qcCt9RLeiiwCyUaNRGxCDiJdDN+NnBSRCwr8/6L+RnwPUl9ilzbBxhVJPw2YHdJ3QvC/gB8SvWGbRVYn/Q5bZTV+YrsiXkldQKuJfW2bAIsJGsQSaoh3ZzvHxE9gd2AyRHxAit/N2s1VEj2vnei4Em/pJ2BQaTGRDFfJ/1G1gZeYdXG0qGU+a7NzKrJDQAza01/z8a0f5A9mf1DqTjAh6Qb1kuKZSRpE9JN288iYlFEPArcWRgnIq7JnnIvIj393k5S71KVi4hREfFW9lT4VtLT27zj+4cA10TExKy8YaSn0JsVxLkoe/K+sEw+U4AlwLMR8WKR6yt9hpJOrvceJgMPAD8sknZdYEaR8Bmk/x8UNhoWkm5qS82TWAycHxGLI+IeUm/F1qXfVuNFxHsR8deIWBARc7P6fLEgyjLgvyR1j4gZEdHUoTr/AzwN3AcgqTPpt3lmmQbYHRHxZEQsITUmB9a7nue7NjOrCjcAzKw1HRoRa9W9gNNLxQHWAM4Exkpav0i8DYHZ9cbgT607kNRZ0i+zIT1zgDeySysNEyqUTficXNBA+a9y8YvUZ3n5ETEPeI/0hLzOmznyGQGMBfpL+nqR6yt9hhFxVZE45wLfltSvXvi7pN6J+jYg3UzPrhd+NdBP0kGrJuG97Oa3zgKgR5F4TSZpTUlXSpqafYePAmtJ6px9718jPe2fkQ3H2aYJZVxC+p6PiojIgk8HnomIcr1PMwuOi733PN+1mVlVuAFgZm1SRCyNiL8BS0kr7tQ3gzQOu6YgbJOC46OBQ4C9SUNVNsvCVVdEYWaSNgWuIjU61skaIVMK4jfkLdJQlbr8aoB1gOmFb6tcBtlwm4OBU4FvkyamFhvKU1bWc/A34Cf1Lj0IHFkkyVGkITML6uXzMWmYywXk/xwq6WxSr8IuEdEL+EIWrqx+90XEPqQGzIuk7w8a+JzrZHM39gf2jYg5BZf2Ag7L5o/MJA0vGqFsxaWcctXBzKwa3AAwszZJySGkMdYv1L8eEVOB8cBwSd0k7UFa1aVOT2AR6Sn8msAv6mXxNrBFwXkN6abtnaz8E0lPhvP6C3CipIGSPpGV90S9FXxKyhoMI4HvR8S72bCaB4DfNqIOhYYDJwJr1QvbTdKFkvpI6inpO8BxFB8yBGmC8BrAl5tYj1Vkk2QH54jakzQU6YOsIXReQR79JB2SfW6LSEOQ6obrvE3qQelWpg7DSI3EvSPivXqXTwA+TRrWM5Dsd8aqDSozs3bJDQAza2vulDQPmEMa8318mbHdR5MmWr5Pujm8oeDaDaQhOdOB51l1MvGfgG2z4T5/j4jnScNvxpFuID9LWjIzl4h4kDQB96+k3oktSRNF8/oF8GJE3FwQ9j1gf0n7FITdqZX3AbijRH1eJ9281xSEvUzqTdmONCRqBnA4sF9EFH2vEbGUNKSo0T0RxSjt6zCXtBRpQy4FupOGLv0buLfgWifSakdvkb7/L5J6TQAeJk3onSmp2Ao+kD7vTYBXCj7LHwNExAcRMbPuBXwMzImID/O/UzOztksrhjyamZm1LEnHAJ+JiGHVrouZWbsjfk+5IYbBd/Nk4w1KzMys1UTETdWug5lZOza+Epm4B8DMzMzMrANxD4CZmZmZWXsi+pIWb9iWtFBDEuyZJ7knAZuZmZmZtS83k1bI25y0StkbwFN5E3sIkJmZmZlZeyImEOyIeIZgQBb2FMFOeZK3qyFAnTp1iu7du1e7GmZmZma2mluwYEFERFsdLbM4+zsDcQBpSeTcyzW3qwZA9+7dmT9/frWrYWZmZmarOUkLq12HMn6O6E3aMf33QC/g+3kTt6shQDU1NeEGgJmZmZm1NEkLIqKm4ZjtT1vt1mgblgFLgGnAPVWui5mZmZkZgLgesVbB+dqIa/ImdwOglOdInSm/BF4CDgK2BD6qZqXMzMzMzBhA8MHys2A2sH3exG4AlLIPaWXVHwN7Au8As4Ejq1kpMzMzMzM6IdZefib60Ii5ve1qEnCr+Q8wA5jMiiZSH+Aq4JYq1cnMzMzMLBkBjEOMAgQcAVyYN7EnARfzPnAjMLQuoDb9WSLo0n4+LzMzM7P2r7Yqpbb5ScBiW9I4lQAeIXg+b1L3ABTTB9ibtMfakLrAWhhE2m/tjupUy8zMzMws05X09L/uODfPAShlHHAsMAA4+quwEfA88POq1srMzMzMOjoxlPSoel1gPeAmxHdyJ/cQoDKeAn4AzHwHPtsXLgU2bL3izczMzKw62vQQIPEMsCvB/Oy8BhhHMCBP8hbvAZB0jaRZkqYUhNVKmi5pcvb6SkvXo0l2Ah4BXrgCbsM3/2ZmZmbWFghYWnC+lBXDgRrUGnMArgMuB26oF/7biPh1K5RvZmZmZrY6uRZ4Ai2fmXoo8Ke8iVu8ARARj0rarKXLMTMzMzPrEILfIMYAe2QhJxJMypu8mqsAnSnpOGA8cHZEzK5iXczMzMzM2ra04VedN7LXimvB+3myqdYqQH8EtgQGkrbcGlEqoqRTJI2XNH7JkiWtVD0zMzMzs/IkdZY0SdJdRa6dIOmdgjmv36pAkRNID88nFByPLzjOpSo9ABHxdt2xpKuAVT60grgjgZGQVgFq+dqZmZmZmeUyFHgB6FXi+q0RcWbFSgs2r0Q2VekBkLRBwelhwJRScc3MzMzM2hpJ/YEDgKurXZfGavEeAEl/AQYD60qaBpwHDJY0kLR18RvAqS1dDzMzMzOzRugiqXBYzchsZEqdS0k7RvUsk8fhkr4AvAR8PyLerHw1G681VgH6RpHg3MsUmZmZmZlVwZKIGFTsgqQDgVkRMUHS4BLp7wT+EhGLJJ0KXA/s2SI1baRqrgJkZmZmZtYe7Q4cnG1muwbQS9JNEXFMXYSIeK8g/tXAxRWtgVgb2JjC+/lgYp6kbgCYmZmZmTVCRAwDhgFkPQDnFN78Z+EbRMSM7PRg0mThyhAXACcAr5KG1JP9zdXD4AaAmZmZmVkFSDofGB8Ro4HvSjoYWAK8T7phr5SjgC0JPm5KYkW0n5U1a2pqYv78+VUouTZ7mZmZmVlHIGlBRNRUux5Fib8C3yaY1ZTk7gEwMzMzM2tfLgImIaYAi5aHBgfnSewGgJmZmZlZ+3I98CvgWWBZYxO7AWBmZmZm1r4sILisqYndADAzMzMza18eQ1wEjGblIUBeBtTMzMzMbDW0ffb3cwVhXgbUzMzMzGy1FHypOck75YkkqZukAZI+K6lbcwo0MzMzM7NmEL0Rv0GMz14jEL3zJm+wASDpANIuY5cBlwOvSNq/6TU2MzMzM7NmuAaYS9oQ7ChgDnBt3sR5hgCNAL4UEa8ASNoSuBv4Z6OramZmZmZmzbUlweEF58MRk/MmzjMEaG7dzX/mNVKLw8zMzMzMWt9CxB7Lz8TuwMK8iUv2AEj6anY4XtI9wG2k2cVHAk81qapmZmZmZtZcpwE3FIz7nw0cnzdxuSFABxUcvw18MTt+B+jemBqamZmZmVkFiM7AsQTbIXoBEMxpTBYlGwARcWLzamdmZmZmZhUVLF0+/KeRN/51GpwELKkvcDKwWWH8iPhmUwo0MzMzM7NmmYQYDYwC5i8PDf6WJ3GeVYD+ATwGPAgsbUIFzczMzMysctYA3mPlnX8DKtcAWDMiftiEipmZmZmZWaWIXxH8ELiHYFRTs8mzDOhdkr7S1ALMzMzMzKwivoIQMKw5mZRbBnQuqStBwI8lLQIWZ+cREb2aU7CZmZmZmTXKvaQlP3uglSYACwiCXPfnJXsAIqJnRPTK/naKiO4F57lv/iVdI2mWpCkFYX0kPSDp5ezv2nnzMzMzMzPrkIL/JlgLuJugV8GrZ96bf8g3BKi5rgO+XC/sR8BDEbEV8FB2bmZmZmZmDQkOaU7yFm8ARMSjwPv1gg8Brs+OrwcObel6mJmZmZlZ6/QAFNMvImZkxzOBfqUiSjpF0nhJ45csWdI6tTMzMzMza4CkzpImSbqrTJzDJYWkQa1Zt3IabABIujFPWFNFRJAmG5e6PjIiBkXEoC5d8qxaamZmZmbWKoYCL5S6KKlnFueJipYqhuYKKyFPD8BnVspb6gzsmLeAEt6WtEGW3wbArGbmZ2ZmZmbWaiT1Bw4Ari4T7QLgV8BHFS7++CJhJ+RNXLIBIGlYthToAElzstdc0s36PxpdzZWNZkXFj69AfmZmZmZmrelS4AfAsmIXJe0AbBwRd1esRPENxJ3A5ojRBa9HWHXObUklx9RExEXARZIuiogmbzYg6S/AYGBdSdOA84BfArdJOgmYChzV1PzNzMzMzFpAF0njC85HRsRIAEkHArMiYoKkwfUTSuoE/IaGnsqLs3LUYz7Bldnx/wIzgHWBEQVx5gLP5MgrFZuG4JeJIO0OTI6I+ZKOAXYAfhcRU/MWUik1NTUxf/781i4WqM1eZmZmZtYRSFoQETUlrl0EHAssAdYAegF/i4hjsuu9gVeBeVmS9UlP6A+OiBWNCjED+CNpI69ShhB8qnnvZmV5ZtX+EdhO0nbA2aRxTjcAX6xkRczMzMzM2oNsdMwwgKwH4Jy6m//s+oekp/RkccZkccazshsJzi9bmFi1ESK+SppbsB6p8VCZnYALLMlW6jkEuDwirgB65snczMzMzKyjkHS+pINzJwh+0MQ4FwMHE/RuqZ2A50oaBhwD3J2NaeqatwAzM2t9b731FkcefSTde3Snx1o9+OYp3+SDDz6odrXMzBq0aNEizq09lw0325A+6/fhhG+dwMyZM6tdrZIiYkxEHJgdnxsRo4vEGVzk6f8K4heItQrO10b8vEyxbxOllx9tSJ45AOsDRwNPRcRjkjYBBkfEDU0ttKk8B8DMrGELFy5kq223YuYmM1m681JYCt3+1Y2tl27N5Ccn06lTtfaANDNr2P4H78+YV8fw0R4fwSegy8QurDdtPV589kV69my9QSjl5gBUvjAmEWxfL2wiwQ4l4v+ONK/g78Ci5eHB3/IU1+D/BSJiZkT8JiIey87/U42bfzMzy+e2227jwzU/ZOmeS6EH0Bs+3v9jXn/ndR5++OFqV8/MrKSnn36aR//3UT766kewIbAOLNlnCR/2/pAbb6zYPrRtUWfEJ5afie5QcL6qXsACYF/goOx1YN7CGpwELOlzwO+BTwPdgM7AvIjonbcQMzNrPU8/8zTzNpi3cqBgcf/FTJkyhb333rs6FTMza8CkSZPQ5lrlDnX+xvMZN34cp3N6dSrW8m4GHkJcm52fCFxfMnZwYnMKy7MK0OXA14FRwCDgOKjsUkRtXe2YMXgIkJm1F692fpWuM7qymMUrAgPizWD8x+OpHVNbtbqZmZUz9cOpLJ6+GIKVFsZc49012Ga3bapWrxYX/ArxNFD3hOYCgvtKxhefIq3U2Y/gvxADSJOCy80bWJE8xxyA8RExSNIzETEgC5sUEduXTdgCqjUHoHbMYGoHj2n1cs3MmmL+/Plsuc2WvPupd5fPAej6eFe2mLcFz09+3nMAzKzNigg+u8NneanHSyzeY3FaduZp6PWvXrz03Ev069ev1erSqnMAAMSmwFYEDyLWBDoTzC0Rdyzw38CVy+cOiCkE/5WnqDz/F1ggqRswWdLFkr6fM52ZmVVBTU0NT/7rSfbruR+dR3Sm62VdOWyLw3j84cd9829mbZokHrn/Efbrux/dLu1G14u7MvDtgYx9cGyr3vy3OnEycDss3/F3I9IE31LWJHiyXtiSvMXlGQJ0LGnc/5nA94GNgcPzFmBmZq1vk0024e6/301dL69UbpNJM7O2o2/fvtz5tztZuHAhixcvplev3Mvbt2dnADsDTwAQvIxYr0z8dxFbkgZLgTgCmJG3sAYbABExNTtcCAzPm7GZmVWfb/zNrL3q3r073bt3r3Y1Wssigo+Xz3sQXai7uS/uDGAksA1iOvA6MCRvYSUbAJKeLVdw3XwAMzMzMzNrlrGIHwPdEfsApwN3lowdvAbsjagBOpWcK1BCuR6A3GuJmpmZmZlZk/0IOAl4FjgVuAe4epVY4iDgGYK6ETpnA4cjpgJDCV7PU1jJBkDB0B8zMzMzM2spwTLgKuAqRB+gP1F0JM6FwOcAEAcCxwDfALYH/gfYL09xXg7CzMzMzKyaxBhEr+zmfwKpIfDbIjGDYEF2/FXgTwQTCK4G+uYtzg0AMzMzM7Pq6k0wh3RTfwPBLsBeReIJ0QPRKbv+UMG1NfIWVrYBIKmzpJvzZmZmZmZmZo3WBbEBcBRwV5l4lwKTgfHACwTjARDbU6llQCNiqaRNJXWLiI/zZmpmZmZmZrmdD9wHPE7wFGIL4OVVYgXXIO4D1gOeLrgyEzgxb2F5NgJ7DfiXpNHA/OXlR/wmbyFmZmZmZlZCMAoYVXD+GqU23g2mA9PrheV++g/55gC8SuqK6AT0LHiZmZmZmVlTiVMqEqeR8uwEPBxA0poRsaCh+GZmZmZmlsuPEO+WuS5gKGnX34ppsAEgaVfgT0APYBNJ2wGnRsTpzS1c0hvAXGApsCQiBjU3TzMzMzOzdmIscFADcR4oe1WcQjSugZBnDsClpE0FRgNExNOSvtCYQhrwpYgo1/IxMzMzM1v9RP6Ju2WcRiN7CHLtAxARb9YLWtqYQszMzMzMrEWosQnyNADelLQbEJK6SjoHeKHRVSsugPslTZBU8QkOZmZmZmaruYaGEK0iTwPgNOAMYCPSkkMDgWaP/8/sERE7APsDZxQbWiTpFEnjJY1fsmRJhYo1MzMzM2uebNPcSZJW2bxL0mmSnpU0WdLjkrYtnRGdm1yJYFpjk+RpAGwdEUMiol9ErBcRxwCfbnztVhUR07O/s4A7gJ2LxBkZEYMiYlCXLnmmLJiZmZmZtYqhlB4Z8+eI+GxEDAQuBsrtofUy4hJE6UZCBeVpAPw+Z1ijSKqR1LPuGNgXmNLcfM3MzMzMWpqk/sABwNXFrkfEnILTGtLQ91K2A14Crkb8G3EKolfxgumE2K1ptU5KPlLPlv/cDegr6ayCS72gGd0UK/QD7pBUV48/R8S9FcjXzMzMzKy5ukgaX3A+MiIKV9u5FPgBZTbIlXQGcBbQDdizZEnBXOAq4CrEF4E/A79F3A5cQPBKQdxliCuA7Rv7huqUG1PTlbT2fxdWfmNzgCOaWmCdiHiN1NrJbdGSZXztynHNLbrR3vjga7zwf61frpmZmZlVTck9qiQdCMyKiAmSBpfKICKuAK6QdDTwU+D4ohHTHIADgBOBzYARwM3A54F7gE/VS/EQ4nDgb0TZnoWiyjUAzouIvSR9pm43YDMzMzMzY3fgYElfAdYAekm6KZsrW8wtwB/L5Pcy8AhwCcH/FoTfjii2/9appJ6FpYiFpKVAgygxbKgeRRRvNEh6HvgWaRfgo6m3xmhETMxTQCXV1NTE/PnzW7tYascMpnbwmFYv18zMzMyqQ9KCiKjJEW8wcE5EHFgvfKuIeDk7Poj0cL1ojwKiB8G8Zlc6p3I9AOcCPwP6s+qs5aDcOCYzMzMzsw5G0vnA+IgYDZwpaW9gMTCbUsN/kisQQwk+SBmxNjCC4JvFC0LAEGBzggsQGwMbEDyZp54lGwARcTtwu6SfRcQFeTIzMzMzM+tIImIMMCY7PrcgfGgjshmw/OYfIJiNyk7y/QOwjPRA/gJgHnAFsFOewhpcWD8iLpC0NrAVaYxTXfijeQowMzMzM7OyOiHWJpgNgOhD+fv0XQh2QEwC6hoM3fIW1mADQNK3SJsc9AcmA58DxuEhQGZmZmZmlTACGIcYRZp3ewRwYZn4i7OVg9JkXtGX1COQS56NwIaSuhOmRsSXSGuOfpC3ADMzMzMzKyO4ATgceBuYCXyV4MYyKS4D7gDWQ1wIPA78Im9xDfYAAB9FxEeSkPSJiHhR0tZ5CzAzMzMzswa9SJosnO7PxSYE/ykaM7gZMQHYi9RjcCjBC3kLytMAmCZpLeDvwAOSZgNT8xZgZmZmZmZliO8A55F6AJZSt64/DCgR/zLgFoIrmlJcnknAh2WHtZIeAXoD9zalMDMzMzMzW8VQYGuC93LGnwD8FLE1aSjQLQTj8xaWpwdguYgY25j4ZmZmZmbWoDeBD3PHDq4Hrs9WCzoc+FU2ZGirPMkb1QAwMzMzM7OKew0Yg7gbWLQ8NFbZjLe+TwLbAJtCZecAmJmZmZlZy/lP9uqWvcoTFwOHAa8CtwAXrLSRWAPy7ANQAyyMiGWSPkVqZfwzIhbnLcTMzMzMzEoIhgMg1iRYUDauEDAX2JXg3aYUl2cfgEeBNSRtBNwPHAtc15TCzMzMzMysHrEr4nnSUqAgtkP8oWjcIICjmnrzD/kaAIqIBcBXgT9ExJHAZ5paoJmZmZmZreRSYD/IVgEKnga+UCb+RMROTS0szxwASdoVGAKclIV1bmqBZmZmZmZWT/AmWilkaZnYuwBDEFOB+dTtGxAl9g2oJ08D4HvAMOCOiHhO0hbAI3kyNzMzMzOzBr2J2A0IRFfSvgDlVvXZrzmF5dkIbCwwtuD8NeC7zSnUzMzMzMyWOw34HbARMJ007/b0krGDqYg9gK0IrkX0BXrkLaxkA0DSnaQtiIuXG3Fw3kLMzMzMzKykrQmGrBQidgf+VTS2OA8YBGwNXAt0BW4Cds9TWLlJwL8GRgCvAwuBq7LXPNKao2ZmZmZm1ny/zxlW5zDgYNL4fwjeAnrmLaxkD0A29AdJIyJiUMGlOyWNz1tAOZK+TOru6AxcHRG/rES+lfLWW2/xox/9iIcnTmThAT9k+PDhrLHGGtWulpmZmZmtDsSuwG5AX8RZBVd6UX7RnY8JAmWjdURNY4rNswxoTTbxN+UvbQ6NK6QYSZ2BK4D9gW2Bb0jatrn5Vsro0aPpv3l/bvz3jUzvO5eLb7iYtdZbi2nTplW7amZmZma2euhGGrvfhfQEv+41BziiTLrbEFcCayFOBh4kjdTJRRElh/mnCOkp/UjgNdISQ5sCp0bEfXkLKZHvrkBtROyXnQ8DiIiLSqWpqamJ+fPnN6fY3NbssyYLd10IO2cBAdwOg3oO4qlxT7VKHczMzMysOiQtiIhmP/TOVxibEkxtZJp9gH1J9+f3ETyQN2meVYDulbQVsE0W9GJELGpUBYvbCHiz4HwaaU3Tqnv55ZdZOHch7FgQKGAPmHjDxGpVy8zMzMxWTwsQl5A2210x3jzYs2jsNOTnYYIHEFsDWyO6EizOU1iefQAg3QpvlsXfThIRcUPOtM0i6RTgFIBu3bq1RpF06ZJ9LPU7R5alHRZqx9S2Sj3MzMzMOrrawbXVrkJruBm4FTiQtCTo8cA7ZeI/CnwesTZwLzAe+BrUW0mohDxDgG4EtgQms2JHsoiIZu0F0NaHAPVctyfztp8He2QBy4BbYI8N9uCxMY+1Sh3MzMzMrDryDAHK5rSOB6ZHxIH1rp0FfAtYQrqZ/2ZEFB/mIyYQ7Ih4ZvluvuIpgp1KxJ9IsAPiO0B3gosRkwkG5nlveXoABgHbRkMthcZ7Ctgqm1Q8Hfg6cHSFy2iyO2+/k72+vBfLXlyWBiv9H/To3IN/jPtHtatmZmZmZm1D3Y69vYpcmwQMiogFkr4NXEx6Sl9M3dCdGYgDgLeAPmXKVbaC0BDgpCys3KpBK8mzCtAUYP28GeYVEUuAM4H7SB/cbRHxXKXLaarBgwcze9Zsvn/k9/nKOp9kxM9G8OGsD+nTp9x3YWZmZmYdgaT+wAHA1cWuR8QjEbEgO/030L9Mdj9H9AbOBs7J8vxemfhDgWHAHQTPIbYAHsld9xxDgB4BBgJPAssn/1ZjJ+DWHAK0strsZWZmZmYdgaSPgWcLgkZGxMiC67cDF5GW7Tyn/hCgenldDsyMiJ/nrwDfI7i0sfXOI88QoNqWKNjMzMzMrA1bUm8z3OUkHQjMiogJkgaXy0TSMaQh9V9sZPlnQZUaABExVlI/WD4J4cmImNUSlTEzMzMzawd2Bw6W9BXSsp29JN0UEccURpK0N/AT4ItNWEZflanqqhqcAyDpKNLwnyOBo4AnJJXbmczMzMzMbLUVEcMion9EbEZayObhIjf/2wNXAgc38eF56XH6YvdcYSXkGQL0E2CnuopL6kvabvj2vIWYmZmZma3uJJ0PjI+I0cAlQA9glCSA/6wyh1bMpfiNvoDuZYr6PbBDjrCi8jQAOtVrtbxHvtWDzMzMzMxWaxExBhiTHZ9bEL53w4np2ajC0tKfuwF9EWcVXOlFI5YBzdMAuFfSfcBfsvOvAf/MW4CZmZmZmVVEN1KvQhdYqfEwB8g9RL/BZUABJH2VFXviPhYRd+SvZ+V4GVAzMzMzaw15dgKuGrEpQfFdhXNosAcg26n3noj4W3beXdJmEfFGUws1MzMzM7MmW4C4BPgMaRWiJNgzT+I8Y/lHAcsKzpdmYWZmZmZm1vpuBl4ENgeGA28AT+VNnKcB0CUiPq47yY67Na6OZmZmZmZWIesQ/AlYTDCW4JuQ7+k/5GsAvCNp+ZJFkg4B3m18Pc3MzMzMrAIWZ39nIA5AbA/0yZs4zypApwE3S7qCtE7pNOC4RlfTzMzMzMwq4eeI3sDZpPX/ewHfz5u4wQZARLwKfE5Sj+x8XhMramZmZmZmzRXclR19CHypsckbHAIkqZ+kPwGjImKepG0lndTYgszMzMzMrPryzAG4DrgP2DA7fwn4XgvVx8zMzMzMWlCeBsC6EXEb2VKgEbGEtBSomZmZmZm1M3kmAc+XtA5pAjCSPkcab2RmZmZmZq1FnFX2evCbPNnkaQCcBYwGtpT0L6AvcESezM3MrIreBvqRBnFuCmxT3eqYmVmz9cz+bg3sRLpHBzgIeDJvJnlWAZoo6YtZQQL+LyIWN5DMzMyq6afA3cBE4C3gWOBc4MxqVsrMzJolGA6AeBTYgWBudl5L+lc/l5JzACTtJGl9WD7uf0fgQmCEpNwbDZiZWSubAFwPPEh6bHMiaYP4c0k7uZiZWXvXD/i44PzjLCyXcpOAr6zLWNIXgF8CN5DG/49sdDXNzKx1/J20XeM6BWGbArewYu9IMzNrz24AnkTUZk//nyA9+sml3BCgzhHxfnb8NWBkRPwV+KukyU2srJmZtbSzgW5FwgcA67dyXczMrPKCCxH/BD6fhZxIMClv8nI9AJ0l1TUQ9gIeLriWZ/JwSZJqJU2XNDl7faU5+ZmZWYHepOdA/ykImwJsC8yuSo3MzKzy1gTmEPwOmIbYPG/CcjfyfwHGSnoXWAg8BiDpk1RmGdDfRsSvK5CPmZkVEmmoz0DgANIuLvcAfwTWrl61zMysQsR5wCDSIj3XAl2Bm4Dd8yQv2QCIiAslPQRsANwfEZFd6gR8pzl1NjOzFvZd4FDg3uz8t8B6VauNmZlV1mHA9qS13iB4Cy1fIrRBZYfyRMS/i4S91MgKlnKmpOOA8cDZEVG0Y1rSKcApAN26FRvUamZmRW1C9q+nmZmtZj4mCJQ26kXUNCZxuTkAzSLpQUlTirwOIXVEb0nqoJ4BjCiVT0SMjIhBETGoS5dmTT0wMzMzM6sYSZ0lTZJ0V5FrX5A0UdISSZXeRPc2xJXAWoiTSQs/X5U3cYvdUUfE3nniSboKWOVDMzMzMzNr44YCLwC9ilz7D3ACcE7FSw1+jdgHmEOaB3AuwQN5k1flkbqkDSJiRnZ6GGl9CjMzMzOzdkFSf9JSCxcCZ9W/HhFvZPGWVb5wzgJubcxNf6Fqjam5WNJAIIA3gFOrVA8zMzMzs6a4FPgB5J98W0E9gfsR7wO3AqMI3s6buCoNgIg4thrlmpmZmZnl1EXS+ILzkRExEkDSgcCsiJggaXCr1ywYDgxHDCBt2DsWMY0g1xB8z6o1MzMzM1vVkogYVOLa7sDB2Wa2awC9JN0UEce0XvUAmAXMBN6jEYs9t9gqQGZmZmZmq6OIGBYR/SNiM+DrwMOtevMvTkeMAR4C1gFOJhiQN7l7AMzMzMzMKkDS+cD4iBgtaSfgDtIe7AdJGh4Rn6lQURsD3yOY3KR6rtjgt+2rqamJ+fPnV6Hk2uxlZmZmZh2BpAUR0agNtlqc6EUwB9Gn6PXg/TzZuAfAzMzMzKx9+DNwIDCBtJqmCq4FsEWeTNwAMDMzMzNrD4IDs7+bNycbTwI2MzMzM2tPhBDHIH6WnW+C2DlvcjcAzMzMzMzalz8AuwJHZ+dzgSvyJvYQIDMzMzOz9mUXgh0QkwAIZiO65U3sHgAzMzMza5s+rnYF2qzFiM6kib8g+gLL8iZ2A8DMzMzM2o6ZwPGkm//LgV2Ap6pao7boMtIeA+shLgQeB36RN7GHAJmZmZlZ27AMOADYj/RseyjQLwt7GtigelVrU4KbEROAvUhLgR5K8ELe5G4AmJmZmVnb8CiwBLiQFSvcDwGeBx7Ojg3E54DniGzir+iF2IXgiVzJvRNwHrVVKNPMzMysg3lyI/hPbzji+YLA2qpUpU3uBFwnTf7dgVg+B6ATMJ5ghzzJ3QOQS221K2BmZma2+tsZ2J40/r9wTZsfAN8EtqlGpdokLb/5BwiWofz39Z4EbGZmZmZtx9nAgcD/As8AZwB3ARtXs1JtzmuI7yK6Zq+hwGt5E7sBYGZmZmZtxwhgb+A04AigKzAWaJuDcarlNGA3YDowjbRW0il5E3sOgJmZmZlZPW16DkAzuQfAzMzMzKw9EZ9CPISYkp0PQPw0b3I3AMzMzMzM2pergGHAYgCCZ4Cv503sBoCZmZmZWfuyJsGT9cKW5E3crpYBXbBgQUhaWO162GqhC434D8WsHfNv3ToK/9at0rpXuwJlvIvYEpbvA3AEMCNv4nY1CdisUiSNj4hB1a6HWUvzb906Cv/WrUMRWwAjSSsBzQZeB4YQTM2TvF31AJiZmZmZdXjBa8DeiBrSkP4FpDkAuRoAngNgZmZmZtYeiF6IYYjLEfuQbvyPB14BjsqbjXsArKMaWe0KmLUS/9ato/Bv3TqCG0lDfsYBJwM/AQQcRjA5byaeA2BmZmZm1h6IZwk+mx13Jk383YTgo8Zk4yFAZmZmZmbtw+LlR8FSYFpjb/7BDQBrByRdI2mWpCkFYRdIekbSZEn3S9qwRNp7JX0g6a564WdKekVSSFq3TNnHS3o5ex1fEL6jpGezPC6TpEq8V+vY2uhvfYyk/8vKnyxpvUq8V+vYqvxbL5V+c0lPZHncKqlbc9+nWQvYDjEne80FBiw/FnPyZuIGgLUH1wFfrhd2SUQMiIiBwF3AuSXSXgIcWyT8X8DelJktL6kPcB6wC7AzcJ6ktbPLfySNvdsqe9Wvn1lTXEfb+60DDImIgdlrVp43YtaA66jCb72B9L8CfhsRnySNsT6pgXzMWl/QmaBX9upJ0KXguFfebNwAsDYvIh4F3q8XVtjKraFuI4xV0z4EzC0SPiki3mig6P2AByLi/YiYDTwAfFnSBkCviPh3pEk0NwCH5nw7ZiW1td96I6pu1ihV/K0XTZ/14u4J3J4FXY//XbfVmBsA1m5JulDSm8AQsidFkgZJuroZeRam3wh4s+DytCxso+y4frhZi6jib73OtdmwjJ95uJu1pFb4rZeyDvBBRNTtJOx/12215gaAtVsR8ZOI2Bi4GTgzCxsfEd9qRp7NSm/WEqr8Wx8SEZ8FPp+9ig2dMKsI/7tu1jrcALDVwc3A4S2Q73Rg44Lz/lnY9Oy4frhZS2vt3zoRUfd3LvBn0hwBs5bWUr/1Ut4D1pJUtz+S/1231ZobANYuSdqq4PQQ4MUWKOY+YF9Ja2cTIvcF7ouIGcAcSZ/LhkMcB/yjBco3q+pvXVKXutVUJHUFDgSmlMnHrMla6bdeVDaf6xHgiCzoePzvuq3G3ACwNk/SX0g73m0taZqkk4BfSpoi6RnSzcrQLO5KYz0lPQaMAvbK0u6XhX9X0jTSU55n6tIUpo+I94ELgKey1/lZGMDpwNWkrbdfBf7Zsp+CdQRt8Lf+CVJD4BlgMumJ6FUt/kHYaq9av/Vy6YEfAmdJeoU0J+BPLfohmFWRdwI2MzMzM+tA3ANgZmZmZtaBuAFgZmZmZtaBuAFgZmZmZtaBuAFgZmZmZtaBuAFgZmZmZtaBuAFgZlYhkpZKmpwtZXinpLUqnP89lchTUq2kc+qF/SSr++SC9zE5W1rxaknbNrdcMzNrG7wMqJlZhUiaFxE9suPrgZci4sIqV2sVkmqBeRHx6xLXl78PMzNb/bgHwMysZYwDNgKQtKWkeyVNkPSYpG2y8H6S7pD0dPbaLQs/RtKT2RP4KyV1zsLfkLSupF9KOqOuoMIn+pL+W9JTkp6RNLwgzk8kvSTpcWDrxrwRSWMkDcqO50m6RNJzkh6UtHN2/TVJB2dxOmdx6upxaha+gaRHC3pJPt/0j9fMzJrKDQAzswrLbtj3AkZnQSOB70TEjsA5wB+y8MuAsRGxHbAD8JykTwNfA3aPiIHAUmBIvSJuBY4qOD8KuFXSvsBWwM7AQGBHSV+QtCPw9SzsK8BOzXh7NcDDEfEZYC7wc2Af4DDg/CzOScCHEbFTVtbJkjYHjgbuy97XdqTdhc3MrJV1qXYFzMxWI90lTSY9+X8BeEBSD2A3YJSkunifyP7uCRwHEBFLgQ8lHQvsCDyVxe8OzCosJCImSVpP0oZAX2B2RLwpaSiwLzApi9qD1CDoCdwREQsAJI2m6T4G7s2OnwUWRcRiSc8Cm2Xh+wIDJB2RnffO6vEUcI2krsDfI2JyM+phZmZN5AaAmVnlLIyIgZLWBO4DzgCuAz7InnrnIeD6iBjWQLxRwBHA+qQegbq0F0XElStlKH0vZ9l5LI4Vk8eWAYsAImKZpLr/p4jU43Ff/cSSvgAcAFwn6TcRcUMF62ZmZjl4CJCZWYVlT9q/C5wNLABel3QkgJLtsqgPAd/OwjtL6p2FHSFpvSy8j6RNixRzK2lYzxGkxgCkRsc3s14HJG2U5fMocKik7pJ6AgdV/E2v7D7g29mTfiR9SlJN9j7ejoirgKtJw57MzKyVuQfAzKwFZMN0ngG+QRrD/0dJPwW6ArcATwNDgZGSTiKN9f92RIzL4t0vqROwmNSTMLVe/s9lN/PTI2JGFnZ/NodgXDZ8aB5wTERMlHRrVuYs0lCclnQ1aTjQRKWKvAMcCgwG/lvS4qxux7VwPczMrAgvA2pmZmZm1oF4CJCZmZmZWQfiBoCZmZmZWQfiBoCZmZmZWQfiBoCZmZmZWQfiBoCZmZmZWQfiBoCZmZmZWQfiBoCZmZmZWQfiBoCZmZmZWQfy/43I9C67qEl2AAAAAElFTkSuQmCC\n",
      "text/plain": [
       "<Figure size 864x216 with 2 Axes>"
      ]
     },
     "metadata": {
      "needs_background": "light"
     },
     "output_type": "display_data"
    }
   ],
   "source": [
    "plot_beats(df, \"XENONnT\", \"test.png\")"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "id": "04fb6678",
   "metadata": {},
   "outputs": [],
   "source": []
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "id": "3e1a1a90",
   "metadata": {},
   "outputs": [],
   "source": []
  }
 ],
 "metadata": {
  "kernelspec": {
   "display_name": "snews",
   "language": "python",
   "name": "snews"
  },
  "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
}