retrieve_examples.ipynb

{
 "cells": [
  {
   "cell_type": "code",
   "execution_count": 1,
   "id": "f50f21d9",
   "metadata": {},
   "outputs": [],
   "source": [
    "%load_ext autoreload\n",
    "%autoreload 2"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 2,
   "id": "909ffefb",
   "metadata": {},
   "outputs": [],
   "source": [
    "%matplotlib inline"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 3,
   "id": "9b03e09f",
   "metadata": {},
   "outputs": [],
   "source": [
    "import sys\n",
    "sys.path.append('code/')"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 4,
   "id": "b3a941c7",
   "metadata": {},
   "outputs": [],
   "source": [
    "from retriever_functions import *"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 5,
   "id": "6ab833f9",
   "metadata": {},
   "outputs": [],
   "source": [
    "citation_data = pd.read_csv('data/FRB_citations_23.csv')\n",
    "abstract_data = pd.read_csv('data/FRB_abstracts.csv')"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 6,
   "id": "7b161df8",
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/html": [
       "<div>\n",
       "<style scoped>\n",
       "    .dataframe tbody tr th:only-of-type {\n",
       "        vertical-align: middle;\n",
       "    }\n",
       "\n",
       "    .dataframe tbody tr th {\n",
       "        vertical-align: top;\n",
       "    }\n",
       "\n",
       "    .dataframe thead th {\n",
       "        text-align: right;\n",
       "    }\n",
       "</style>\n",
       "<table border=\"1\" class=\"dataframe\">\n",
       "  <thead>\n",
       "    <tr style=\"text-align: right;\">\n",
       "      <th></th>\n",
       "      <th>txt_ref</th>\n",
       "      <th>arxiv_id</th>\n",
       "      <th>reasons</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>0</th>\n",
       "      <td>lorimer et al. 2007</td>\n",
       "      <td>arXiv:0709.4301</td>\n",
       "      <td>these references are cited to provide backgrou...</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1</th>\n",
       "      <td>petroff et al. 2022; petroff et al. 2021</td>\n",
       "      <td>arXiv:2107.10113</td>\n",
       "      <td>these references are cited to provide backgrou...</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2</th>\n",
       "      <td>cordes &amp; chatterjee 2019</td>\n",
       "      <td>arXiv:1906.05878</td>\n",
       "      <td>these references are cited to provide backgrou...</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>3</th>\n",
       "      <td>zhang 2022; zhang 2022a; zhang 2022b</td>\n",
       "      <td>arXiv:2212.03972</td>\n",
       "      <td>these references are cited to provide backgrou...</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>4</th>\n",
       "      <td>margalit et al. 2019</td>\n",
       "      <td>arXiv:1907.00016</td>\n",
       "      <td>these references are cited to support the conn...</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>...</th>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2361</th>\n",
       "      <td>hankins et al. 2016</td>\n",
       "      <td>arXiv:1608.08881</td>\n",
       "      <td>referenced to mention the observation of main ...</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2362</th>\n",
       "      <td>mondal et al. 2020</td>\n",
       "      <td>arXiv:2008.11375</td>\n",
       "      <td>cited to mention the possible turnover in the ...</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2363</th>\n",
       "      <td>mcconnell &amp; ma 2013</td>\n",
       "      <td>arXiv:1211.2816</td>\n",
       "      <td>these papers are cited to provide background i...</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2364</th>\n",
       "      <td>reines &amp; volonteri 2015</td>\n",
       "      <td>arXiv:1508.06274</td>\n",
       "      <td>these papers are cited to provide background i...</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2365</th>\n",
       "      <td>cook et al. 2019</td>\n",
       "      <td>arXiv:1710.05016</td>\n",
       "      <td>referenced to describe the palomar observatory...</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "<p>2366 rows × 3 columns</p>\n",
       "</div>"
      ],
      "text/plain": [
       "                                       txt_ref          arxiv_id  \\\n",
       "0                          lorimer et al. 2007   arXiv:0709.4301   \n",
       "1     petroff et al. 2022; petroff et al. 2021  arXiv:2107.10113   \n",
       "2                     cordes & chatterjee 2019  arXiv:1906.05878   \n",
       "3         zhang 2022; zhang 2022a; zhang 2022b  arXiv:2212.03972   \n",
       "4                         margalit et al. 2019  arXiv:1907.00016   \n",
       "...                                        ...               ...   \n",
       "2361                       hankins et al. 2016  arXiv:1608.08881   \n",
       "2362                        mondal et al. 2020  arXiv:2008.11375   \n",
       "2363                       mcconnell & ma 2013   arXiv:1211.2816   \n",
       "2364                   reines & volonteri 2015  arXiv:1508.06274   \n",
       "2365                          cook et al. 2019  arXiv:1710.05016   \n",
       "\n",
       "                                                reasons  \n",
       "0     these references are cited to provide backgrou...  \n",
       "1     these references are cited to provide backgrou...  \n",
       "2     these references are cited to provide backgrou...  \n",
       "3     these references are cited to provide backgrou...  \n",
       "4     these references are cited to support the conn...  \n",
       "...                                                 ...  \n",
       "2361  referenced to mention the observation of main ...  \n",
       "2362  cited to mention the possible turnover in the ...  \n",
       "2363  these papers are cited to provide background i...  \n",
       "2364  these papers are cited to provide background i...  \n",
       "2365  referenced to describe the palomar observatory...  \n",
       "\n",
       "[2366 rows x 3 columns]"
      ]
     },
     "execution_count": 6,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "citation_data"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 7,
   "id": "5dc7963e",
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/html": [
       "<div>\n",
       "<style scoped>\n",
       "    .dataframe tbody tr th:only-of-type {\n",
       "        vertical-align: middle;\n",
       "    }\n",
       "\n",
       "    .dataframe tbody tr th {\n",
       "        vertical-align: top;\n",
       "    }\n",
       "\n",
       "    .dataframe thead th {\n",
       "        text-align: right;\n",
       "    }\n",
       "</style>\n",
       "<table border=\"1\" class=\"dataframe\">\n",
       "  <thead>\n",
       "    <tr style=\"text-align: right;\">\n",
       "      <th></th>\n",
       "      <th>arxiv_id</th>\n",
       "      <th>authors</th>\n",
       "      <th>title</th>\n",
       "      <th>abstract</th>\n",
       "      <th>category</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>0</th>\n",
       "      <td>arXiv:2308.14387</td>\n",
       "      <td>Pavan A. Uttarkar et al.</td>\n",
       "      <td>Searching for the spectral depolarisation of A...</td>\n",
       "      <td>Fast Radio Bursts (FRBs) are extragalactic tra...</td>\n",
       "      <td>astro-ph.HE</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1</th>\n",
       "      <td>arXiv:2308.14325</td>\n",
       "      <td>Zhao-Yang Xia et al.</td>\n",
       "      <td>Polarization Evolution of Fast Radio Burst Sou...</td>\n",
       "      <td>Recently, some fast radio bursts (FRBs) have b...</td>\n",
       "      <td>astro-ph.HE</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2</th>\n",
       "      <td>arXiv:2308.12801</td>\n",
       "      <td>Shivani Bhandari et al.</td>\n",
       "      <td>Constraints on the persistent radio source ass...</td>\n",
       "      <td>We present very-long-baseline interferometry (...</td>\n",
       "      <td>astro-ph.HE</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>3</th>\n",
       "      <td>arXiv:2308.12432</td>\n",
       "      <td>Evan F. Lewis et al.</td>\n",
       "      <td>The Petabyte Project</td>\n",
       "      <td>Transient radio sources, such as fast radio bu...</td>\n",
       "      <td>astro-ph.IM</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>4</th>\n",
       "      <td>arXiv:2308.12118</td>\n",
       "      <td>Danté M. Hewitt et al.</td>\n",
       "      <td>Dense Forests of Microshots in Bursts from FRB...</td>\n",
       "      <td>We report on exceptionally bright bursts (&gt;400...</td>\n",
       "      <td>astro-ph.HE</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>...</th>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1348</th>\n",
       "      <td>arXiv:1611.06481</td>\n",
       "      <td>He Gao, Bing Zhang,</td>\n",
       "      <td>Implications from the upper limit of radio aft...</td>\n",
       "      <td>A $\\gamma$-ray transient, Swift J0644.5-5111, ...</td>\n",
       "      <td>astro-ph.HE</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1349</th>\n",
       "      <td>arXiv:1610.09043</td>\n",
       "      <td>S. Johnston et al.</td>\n",
       "      <td>Radio light curve of the galaxy possibly assoc...</td>\n",
       "      <td>We present observations made with the Australi...</td>\n",
       "      <td>astro-ph.GA</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1350</th>\n",
       "      <td>arXiv:1603.04825</td>\n",
       "      <td>Ye Li, Bing Zhang,</td>\n",
       "      <td>Radio transient following FRB 150418: afterglo...</td>\n",
       "      <td>Recently, Keane et al. reported the discovery ...</td>\n",
       "      <td>astro-ph.HE</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1351</th>\n",
       "      <td>arXiv:1602.08086</td>\n",
       "      <td>Bing Zhang,</td>\n",
       "      <td>On the afterglow and progenitor of FRB 150418</td>\n",
       "      <td>Keane et al. recently detected a fading radio ...</td>\n",
       "      <td>astro-ph.HE</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1352</th>\n",
       "      <td>arXiv:1309.4451</td>\n",
       "      <td>Matthew McQuinn,</td>\n",
       "      <td>Locating the \"missing\" baryons with extragalac...</td>\n",
       "      <td>Recently, Thornton and coworkers (2013) confir...</td>\n",
       "      <td>astro-ph.CO</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "<p>1353 rows × 5 columns</p>\n",
       "</div>"
      ],
      "text/plain": [
       "              arxiv_id                   authors  \\\n",
       "0     arXiv:2308.14387  Pavan A. Uttarkar et al.   \n",
       "1     arXiv:2308.14325      Zhao-Yang Xia et al.   \n",
       "2     arXiv:2308.12801   Shivani Bhandari et al.   \n",
       "3     arXiv:2308.12432      Evan F. Lewis et al.   \n",
       "4     arXiv:2308.12118    Danté M. Hewitt et al.   \n",
       "...                ...                       ...   \n",
       "1348  arXiv:1611.06481      He Gao, Bing Zhang,    \n",
       "1349  arXiv:1610.09043        S. Johnston et al.   \n",
       "1350  arXiv:1603.04825       Ye Li, Bing Zhang,    \n",
       "1351  arXiv:1602.08086              Bing Zhang,    \n",
       "1352   arXiv:1309.4451         Matthew McQuinn,    \n",
       "\n",
       "                                                  title  \\\n",
       "0     Searching for the spectral depolarisation of A...   \n",
       "1     Polarization Evolution of Fast Radio Burst Sou...   \n",
       "2     Constraints on the persistent radio source ass...   \n",
       "3                                  The Petabyte Project   \n",
       "4     Dense Forests of Microshots in Bursts from FRB...   \n",
       "...                                                 ...   \n",
       "1348  Implications from the upper limit of radio aft...   \n",
       "1349  Radio light curve of the galaxy possibly assoc...   \n",
       "1350  Radio transient following FRB 150418: afterglo...   \n",
       "1351      On the afterglow and progenitor of FRB 150418   \n",
       "1352  Locating the \"missing\" baryons with extragalac...   \n",
       "\n",
       "                                               abstract     category  \n",
       "0     Fast Radio Bursts (FRBs) are extragalactic tra...  astro-ph.HE  \n",
       "1     Recently, some fast radio bursts (FRBs) have b...  astro-ph.HE  \n",
       "2     We present very-long-baseline interferometry (...  astro-ph.HE  \n",
       "3     Transient radio sources, such as fast radio bu...  astro-ph.IM  \n",
       "4     We report on exceptionally bright bursts (>400...  astro-ph.HE  \n",
       "...                                                 ...          ...  \n",
       "1348  A $\\gamma$-ray transient, Swift J0644.5-5111, ...  astro-ph.HE  \n",
       "1349  We present observations made with the Australi...  astro-ph.GA  \n",
       "1350  Recently, Keane et al. reported the discovery ...  astro-ph.HE  \n",
       "1351  Keane et al. recently detected a fading radio ...  astro-ph.HE  \n",
       "1352  Recently, Thornton and coworkers (2013) confir...  astro-ph.CO  \n",
       "\n",
       "[1353 rows x 5 columns]"
      ]
     },
     "execution_count": 7,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "abstract_data"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "id": "e62065d3",
   "metadata": {},
   "outputs": [],
   "source": []
  },
  {
   "cell_type": "code",
   "execution_count": 8,
   "id": "4b2103c7",
   "metadata": {},
   "outputs": [],
   "source": [
    "retriever_abstract = gen_retriever('VectorStoreIndex', 'abstract', path_to_db_folder='data/')\n",
    "retriever_citation = gen_retriever('VectorStoreIndex', 'citation', path_to_db_folder='data/')"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "id": "c1aca3f5",
   "metadata": {},
   "outputs": [],
   "source": []
  },
  {
   "cell_type": "markdown",
   "id": "3f3f1b7e",
   "metadata": {},
   "source": [
    "### an example where citation clearly gives more desirable results"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 9,
   "id": "6fe555c3",
   "metadata": {},
   "outputs": [],
   "source": [
    "q = 'what is fast radio burst (frb)'"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 10,
   "id": "000f51bb",
   "metadata": {},
   "outputs": [],
   "source": [
    "rst_all_abstract = retriever_abstract.retrieve(q)\n",
    "rst_all_citation = retriever_citation.retrieve(q)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 11,
   "id": "766668d7",
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAjAAAAGdCAYAAAAMm0nCAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjcuMywgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/OQEPoAAAACXBIWXMAAA9hAAAPYQGoP6dpAAA2DklEQVR4nO3deXxU5d3//3cyyUwWSHIHyEJlkyoQAReWMJVaKikBYsUSW6kUolL4isFfJYoab5TFllD0IS5lae9bgd43lBbrUgFZRMFagmAslUWoIBYsTEC4yUggk2XO74/A0KkZzAyz5GRez8fjPB7XnHOuk8+5uJh8cp3lijEMwxAAAICJxEY6AAAAAH+RwAAAANMhgQEAAKZDAgMAAEyHBAYAAJgOCQwAADAdEhgAAGA6JDAAAMB04iIdQCDcbreOHj2qtm3bKiYmJtLhAACAZjAMQ19++aU6duyo2NjLG0MxZQJz9OhRderUKdJhAACAABw5ckRXXHHFZR3DlAlM27ZtJTU2QEpKSoSjAdBqVVdLHTs2lo8elZKTIxsPokug/a8F91un06lOnTp5fo9fDlMmMBcuG6WkpJDAAAgdi+ViOSWlRf0iQBQItP+ZoN8G4/YPbuIFAACmY8oRGAAIi7g4qajoYhkIp0D7X5T02xjDMIxIB+Evp9Op1NRUVVVVcQkJAACTCObvby4hAQAA02m9Y0sAcLkMQzp7trGclCTx3imEU6D9L0r6LSMwAODL2bNSmzaNy4VfCEC4BNr/oqTfksAAAADTIYEBAACmQwIDAABMhwQGAACYDgkMAAAwHRIYAABgOrwHBgB8sVik22+/WAbCKdD+FyX9lqkEAABAWATz9zcjMAAQJl0fXfO1+3w2tyAMkQDmd1n3wMydO1cxMTF64IEHPOtqampUXFysdu3aqU2bNiosLFRlZaVXvcOHD6ugoEBJSUnKyMjQtGnTVF9ffzmhAACAKBJwArNjxw79+te/Vt++fb3WT506VW+88YZWrVqlLVu26OjRoxo9erRne0NDgwoKClRbW6utW7dq2bJlWrp0qZ544onAzwIAQqG6unEemZiYxjIQToH2vyjptwElMGfOnNHYsWP1X//1X/qP//gPz/qqqiq9+OKLeuaZZ3TzzTerX79+WrJkibZu3apt27ZJkjZs2KC9e/fqf//3f3XddddpxIgRevLJJ7VgwQLV1tYG56wAAECrFlACU1xcrIKCAuXl5Xmtr6ioUF1dndf6nj17qnPnziovL5cklZeXq0+fPsrMzPTsk5+fL6fTqT179jT581wul5xOp9cCAACil9838a5cuVIffvihduzY8ZVtDodDVqtVaWlpXuszMzPlcDg8+/xr8nJh+4VtTSkrK9OsWbP8DRUAALRSfo3AHDlyRD/72c+0fPlyJSQkhCqmrygtLVVVVZVnOXLkSNh+NgAAaHn8SmAqKip0/Phx3XDDDYqLi1NcXJy2bNmi559/XnFxccrMzFRtba1Onz7tVa+yslJZWVmSpKysrK88lXTh84V9/p3NZlNKSorXAgAAopdfCczQoUO1a9cu7dy507P0799fY8eO9ZTj4+O1adMmT539+/fr8OHDstvtkiS73a5du3bp+PHjnn02btyolJQU5eTkBOm0AABAa+bXPTBt27ZV7969vdYlJyerXbt2nvUTJkxQSUmJ0tPTlZKSovvvv192u12DBg2SJA0bNkw5OTkaN26c5s2bJ4fDoenTp6u4uFg2my1IpwUAQWCxSCNHXiwD4RRo/4uSfhv0N/HOnz9fsbGxKiwslMvlUn5+vhYuXOjZbrFYtHr1ak2ePFl2u13JyckqKirS7Nmzgx0KAFyehARpzde/PRcIiUD7X5T0W+ZCAoAwYSoBRLtg/v6+rKkEAAAAIoEEBgB8qa6WkpMbl1b8Sna0UIH2vyjpt8xGDQCXcvZspCNANAu0/0VBv2UEBgAAmA4JDAAAMB0SGAAAYDokMAAAwHRIYAAAgOnwFBIA+BIbK33nOxfLQDgF2v+ipN+SwACAL4mJ0ubNkY4C0SrQ/hcl/bb1pmYAAKDVIoEBAACmQwIDAL5UV0sdOjQurfiV7GihAu1/UdJvuQcGAC7liy8kSb0eX6dz1gSfuzGLNELifP8LWz0TYQQGAACYDgkMAAAwHRIYAABgOiQwAADAdEhgAACA6fAUEgD4Ehsr9e+vv31eJXdMzCV37fromjAFhahxvv95yqGuZzIkMADgS2KitGOHRpGcIBLO97+w1TOZ1puaAQCAVosEBgAAmA4JDAD4cvas1LWr3lt0jxLqaiIdDaLN+f6nrl0by6GuZzLcAwMAvhiG9I9/6ApJMUakg0HUOd//POVQ1zMZRmAAAIDpkMAAAADTIYEBAACmQwIDAABMhwQGAACYDk8hAYAvMTFSTo7+XnlGxqVnEgCC73z/85RDXc9kSGAAwJekJGnPHg1jKgFEwvn+F7Z6JkMCAwAtSHMmhfxsbkEYIgFaNr/ugVm0aJH69u2rlJQUpaSkyG6368033/RsHzJkiGJiYryWe++91+sYhw8fVkFBgZKSkpSRkaFp06apvr4+OGcDAACigl8jMFdccYXmzp2rq666SoZhaNmyZRo1apT++te/6pprrpEkTZw4UbNnz/bUSUpK8pQbGhpUUFCgrKwsbd26VceOHdP48eMVHx+vOXPmBOmUACBIzp6VBgzQhsozurXoGdXEJ0Q6IkST8/1PUuPs0v/y+zQk9UzGrwTm+9//vtfnX/ziF1q0aJG2bdvmSWCSkpKUlZXVZP0NGzZo7969euutt5SZmanrrrtOTz75pB555BHNnDlTVqs1wNMAgBAwDGnvXl0tphJABJzvf55yqOuZTMCPUTc0NGjlypWqrq6W3W73rF++fLnat2+v3r17q7S0VGf/ZSKp8vJy9enTR5mZmZ51+fn5cjqd2nOJG45cLpecTqfXAgAAopffN/Hu2rVLdrtdNTU1atOmjV599VXlnH9c684771SXLl3UsWNHffTRR3rkkUe0f/9+vfLKK5Ikh8PhlbxI8nx2OBw+f2ZZWZlmzZrlb6gAAKCV8juB6dGjh3bu3Kmqqiq9/PLLKioq0pYtW5STk6NJkyZ59uvTp4+ys7M1dOhQHTx4UN27dw84yNLSUpWUlHg+O51OderUKeDjAQAAc/P7EpLVatU3v/lN9evXT2VlZbr22mv13HPPNblvbm6uJOnAgQOSpKysLFVWVnrtc+Gzr/tmJMlms3mefLqwAACA6HXZUwm43W65XK4mt+3cuVOSlJ2dLUmy2+3atWuXjh8/7tln48aNSklJ8VyGAgAA+Dp+XUIqLS3ViBEj1LlzZ3355ZdasWKFNm/erPXr1+vgwYNasWKFRo4cqXbt2umjjz7S1KlTddNNN6lv376SpGHDhiknJ0fjxo3TvHnz5HA4NH36dBUXF8tms4XkBAEgYDExUpcu+vz/zjGVAMLvfP/zlENdz2T8SmCOHz+u8ePH69ixY0pNTVXfvn21fv16fe9739ORI0f01ltv6dlnn1V1dbU6deqkwsJCTZ8+3VPfYrFo9erVmjx5sux2u5KTk1VUVOT13hgAaDGSkqTPPtNgphJAJJzvf2GrZzIxhmG+h8SdTqdSU1NVVVXF/TAAQq45r/cPJ6YSgFkF8/f3Zd8DAwAAEG4kMADgy7lz0oABen3ZVNnqmn5YAQiZ8/1PAwY0lkNdz2SYjRoAfHG7pQ8+0LWSYs13tR1md77/ecqhrmcyjMAAAADTIYEBAACmQwIDAABMhwQGAACYDgkMAAAwHZ5CAoBLad9eJ6trIx0FolX79uGtZyIkMADgS3KydOKE+rWwN/EiSpzvf2GrZzJcQgIAAKZDAgMAAEyHBAYAfDl3ThoyRCtXPMpUAgi/8/1PQ4b4P5VAIPVMhntgAMAXt1vaskWDxFQCiIDz/c9TDnU9k2EEBgAAmA4JDAAAMB0SGAAAYDokMAAAwHRIYAAAgOnwFBIAXEpSks7WNUQ6CkSrpKTw1jMREhgA8CU5WaquVg5TCSASzve/sNUzGS4hAQAA0yGBAQAApkMCAwC+1NRIBQV6adVM2eprIx0Nos35/qeCgsZyqOuZDPfAAIAvDQ3S2rW6WVJsK34lO1qo8/3PUw51PZNhBAYAAJgOCQwAADAdEhgAAGA6JDAAAMB0SGAAAIDpkMAAAADT4TFqAPAlOVkyDHVlKgFEwvn+F7Z6JuPXCMyiRYvUt29fpaSkKCUlRXa7XW+++aZne01NjYqLi9WuXTu1adNGhYWFqqys9DrG4cOHVVBQoKSkJGVkZGjatGmqr68PztkAAICo4FcCc8UVV2ju3LmqqKjQBx98oJtvvlmjRo3Snj17JElTp07VG2+8oVWrVmnLli06evSoRo8e7anf0NCggoIC1dbWauvWrVq2bJmWLl2qJ554IrhnBQAAWrUYw7i8cab09HQ99dRTuv3229WhQwetWLFCt99+uyRp37596tWrl8rLyzVo0CC9+eabuuWWW3T06FFlZmZKkhYvXqxHHnlEJ06ckNVqbdbPdDqdSk1NVVVVlVJSUi4nfADwraZGGjdOa3YdU8ktD8oV17zvqFD7bG5BpENAOJzvf5Kk//kfKSEhtPXCIJi/vwO+ibehoUErV65UdXW17Ha7KioqVFdXp7y8PM8+PXv2VOfOnVVeXi5JKi8vV58+fTzJiyTl5+fL6XR6RnGa4nK55HQ6vRYACLmGBunll1Ww/y9MJYDwO9//9PLL/k8lEEg9k/H7Jt5du3bJbrerpqZGbdq00auvvqqcnBzt3LlTVqtVaWlpXvtnZmbK4XBIkhwOh1fycmH7hW2+lJWVadasWf6GCgBf61I36CbW1ujjMMYCoPn8HoHp0aOHdu7cqffff1+TJ09WUVGR9u7dG4rYPEpLS1VVVeVZjhw5EtKfBwAAWja/R2CsVqu++c1vSpL69eunHTt26LnnntMdd9yh2tpanT592msUprKyUllZWZKkrKwsbd++3et4F55SurBPU2w2m2w2m7+hAgCAVuqyX2TndrvlcrnUr18/xcfHa9OmTZ5t+/fv1+HDh2W32yVJdrtdu3bt0vHjxz37bNy4USkpKcrJybncUAAAQJTwawSmtLRUI0aMUOfOnfXll19qxYoV2rx5s9avX6/U1FRNmDBBJSUlSk9PV0pKiu6//37Z7XYNGjRIkjRs2DDl5ORo3LhxmjdvnhwOh6ZPn67i4mJGWAAAQLP5lcAcP35c48eP17Fjx5Samqq+fftq/fr1+t73vidJmj9/vmJjY1VYWCiXy6X8/HwtXLjQU99isWj16tWaPHmy7Ha7kpOTVVRUpNmzZwf3rAAAQKt22e+BiQTeAwMgWC45TYBhKLHOJUk6F2+TYmLCFNWl8R6YKGEY0tmzjeWkpOb3v0DrhUEwf38zFxIA+BITo3PWlvMSMESZmJjGeY3CVc9kmI0aAACYDgkMAPhgra/T02vm6+k182Wtr4t0OIg2Lpd0112Ni8sV+nomQwIDAD5Y3A26ffcm3b57kyzu1vtKdrRQ9fXSsmWNS3196OuZDAkMAAAwHW7iBQCTueSTU+fxpBJaO0ZgAACA6ZDAAAAA0yGBAQAApkMCAwAATIebeAHAh3PxNt1w/3JPGQirpCTp+PGL5VDXMxkSGADwJSZGp5JSIx0FolVMjNShQ/jqmQyXkAAAgOmQwACAD9b6Os3esEizNyxiKgGEn8slFRc3Lv5OJRBIPZMhgQEAHyzuBo3/6xqN/+saphJA+NXXSwsXNi7+TiUQSD2TIYEBAACmQwIDAABMhwQGAACYDgkMAAAwHRIYAABgOiQwAADAdHgTLwD4UBNv1eB7X/SUgbBKTJQOHbpYDnU9kyGBAQAfjJhYfZ6aGekwEK1iY6WuXcNXz2S4hAQAAEyHBAYAfIhvqFPpOy+p9J2XFN/AVAIIs9paadq0xqW2NvT1TIYEBgB8iGto0P/b/or+3/ZXFNfAVAIIs7o66emnG5c6PxLoQOuZDAkMAAAwHRIYAABgOiQwAADAdEhgAACA6fAeGABohbo+uuZr9/lsbkEYIgFCgxEYAABgOn4lMGVlZRowYIDatm2rjIwM3Xbbbdq/f7/XPkOGDFFMTIzXcu+993rtc/jwYRUUFCgpKUkZGRmaNm2a6uvrL/9sACCIauKt+t49C/S9exYwlQDCLzFR2r27cfF3KoFA6pmMX5eQtmzZouLiYg0YMED19fV67LHHNGzYMO3du1fJycme/SZOnKjZs2d7PiclJXnKDQ0NKigoUFZWlrZu3apjx45p/Pjxio+P15w5c4JwSgAQHEZMrD7p0CXSYSBaxcZK11wTvnom41cCs27dOq/PS5cuVUZGhioqKnTTTTd51iclJSkrK6vJY2zYsEF79+7VW2+9pczMTF133XV68skn9cgjj2jmzJmyWvkrBwAAXNpl3QNTVVUlSUpPT/dav3z5crVv3169e/dWaWmpzp4969lWXl6uPn36KDPz4gRp+fn5cjqd2rNnz+WEAwBBFd9QpwfeW64H3lvOVAIIv9paaebMxsXfqQQCqWcyAT+F5Ha79cADD+jGG29U7969PevvvPNOdenSRR07dtRHH32kRx55RPv379crr7wiSXI4HF7JiyTPZ4fD0eTPcrlccrlcns9OpzPQsAGg2eIaGvTAX34nSfr1wELVWeIjHBGiSl2dNGtWY3naNKm5VygCrWcyAScwxcXF2r17t9577z2v9ZMmTfKU+/Tpo+zsbA0dOlQHDx5U9+7dA/pZZWVlmnXhHwMAAES9gC4hTZkyRatXr9Y777yjK6644pL75ubmSpIOHDggScrKylJlZaXXPhc++7pvprS0VFVVVZ7lyJEjgYQNAABaCb8SGMMwNGXKFL366qt6++231a1bt6+ts3PnTklSdna2JMlut2vXrl06fvy4Z5+NGzcqJSVFOTk5TR7DZrMpJSXFawEAANHLr0tIxcXFWrFihV5//XW1bdvWc89KamqqEhMTdfDgQa1YsUIjR45Uu3bt9NFHH2nq1Km66aab1LdvX0nSsGHDlJOTo3HjxmnevHlyOByaPn26iouLZbPZgn+GAACg1fFrBGbRokWqqqrSkCFDlJ2d7Vl+//vfS5KsVqveeustDRs2TD179tSDDz6owsJCvfHGG55jWCwWrV69WhaLRXa7XT/5yU80fvx4r/fGAAAAXIpfIzCGYVxye6dOnbRly5avPU6XLl20du1af340AACAB5M5AoAPrrh43Tr+GU8ZCKuEBGn79ovlUNczGRIYAPDBHWvRR9lXRzoMRCuLRRowIHz1TIbZqAEAgOkwAgMAPsQ31OnuD/4kSVrS/1bexIvwqq2VnnuusfyznzX/jbqB1jMZEhgA8CGuoUGPbV4iSfqf6wtIYBBedXXSww83lu+7z7+pBAKpZzIkMABara6Prol0CABChHtgAACA6ZDAAAAA0yGBAQAApkMCAwAATIcEBgAAmA5PIQGAD664eI358RxPGQirhATpnXculkNdz2RIYADAB3esRds69410GIhWFos0ZEj46pkMl5AAAIDpMAIDAD7ENdTrx39bJ0n63bXDVW/hKxNhVFcn/eY3jeVJk6T4Zl7GDLSeyfC/EQB8iG+o15MbF0uSXu6dRwKD8KqtlaZMaSzfdVfzE5FA65kMl5AAAIDpkMAAAADTIYEBAACmQwIDAABMhwQGAACYDgkMAAAwHZ4JBAAfauPidfftMzxlIKxsNmn16ovlUNczGRIYAPChIdaid7oPiHQYiFZxcVJBQfjqmQyXkAAAgOkwAgMAPsQ11Ou2vZslSa/lDOFNvAivujpp+fLG8tix/k0lEEg9k+F/IwD4EN9Qr6fXPitJWtNjMAkMwqu2Vrr77sbyD3/o31QCgdQzGS4hAQAA0yGBAQAApkMCAwAATIcEBgAAmA4JDAAAMB0SGAAAYDp+PRNYVlamV155Rfv27VNiYqK+9a1v6Ze//KV69Ojh2aempkYPPvigVq5cKZfLpfz8fC1cuFCZmZmefQ4fPqzJkyfrnXfeUZs2bVRUVKSysjLFxfGIItDadX10zdfu89nclvEW0dq4eN036lFPubUx079FVLLZpD/84WI51PVMxq+MYcuWLSouLtaAAQNUX1+vxx57TMOGDdPevXuVnJwsSZo6darWrFmjVatWKTU1VVOmTNHo0aP1l7/8RZLU0NCggoICZWVlaevWrTp27JjGjx+v+Ph4zZkzJ/hnCAABaoi1aG3PwZEOA9EqLq7xPS7hqmcyMYZhGIFWPnHihDIyMrRlyxbddNNNqqqqUocOHbRixQrdfvvtkqR9+/apV69eKi8v16BBg/Tmm2/qlltu0dGjRz2jMosXL9YjjzyiEydOyGq1fu3PdTqdSk1NVVVVlVJSUgINH0AEhPOv/ub8LFwaIzAIpmD+/r6se2CqqqokSenp6ZKkiooK1dXVKS8vz7NPz5491blzZ5WXl0uSysvL1adPH69LSvn5+XI6ndqzZ0+TP8flcsnpdHotABBqFneDRu57TyP3vSeLuyHS4SDa1NdLq1Y1LvX1oa9nMgHfdOJ2u/XAAw/oxhtvVO/evSVJDodDVqtVaWlpXvtmZmbK4XB49vnX5OXC9gvbmlJWVqZZs2YFGioABMRaX6eFr8+VJPWa+rLOWS0RjghRxeWSfvSjxvKZM42XhkJZz2QCHoEpLi7W7t27tXLlymDG06TS0lJVVVV5liNHjoT8ZwIAgJYroLRsypQpWr16td59911dccUVnvVZWVmqra3V6dOnvUZhKisrlZWV5dln+/btXserrKz0bGuKzWaTrRXfSQ0AAPzj1wiMYRiaMmWKXn31Vb399tvq1q2b1/Z+/fopPj5emzZt8qzbv3+/Dh8+LLvdLkmy2+3atWuXjh8/7tln48aNSklJUU5OzuWcCwAAiBJ+jcAUFxdrxYoVev3119W2bVvPPSupqalKTExUamqqJkyYoJKSEqWnpyslJUX333+/7Ha7Bg0aJEkaNmyYcnJyNG7cOM2bN08Oh0PTp09XcXExoywAAKBZ/EpgFi1aJEkaMmSI1/olS5borrvukiTNnz9fsbGxKiws9HqR3QUWi0WrV6/W5MmTZbfblZycrKKiIs2ePfvyzgRAq8EL1gB8Hb8SmOa8MiYhIUELFizQggULfO7TpUsXrV271p8fDQAA4NE6n60CgCCos8TpoZEPeMpAWFmt0pIlF8uhrmcy/I8EAB/qLXF6uU/e1+8IhEJ8vHT+9oyw1DMZZqMGAACmwwgMAPhgcTfopkMfSpLe7XaDGmJ5Ey/CqL5eWr++sZyf3/w36gZaz2Ra51kBQBBY6+u05OXGaUyYSgBh53JJt9zSWPZ3KoFA6pkMl5AAAIDpkMAAAADTaZ3jSkCI8II1AGgZGIEBAACmQwIDAABMhwQGAACYDvfAAIAPdZY4Pf69ez1lIKysVulXv7pYDnU9k+F/JAD4UG+J0//ccEukw0C0io+XiovDV89kSGAAmFJznggD0HqRwACAD7HuBg38fI8kafsV18jNVAIIp4YG6c9/bix/+9uSpZn9L9B6JkMCA5wXzr/oeZ+MOdjq67Tyd49JYioBREBNjfTd7zaWz5yRkpNDW89keAoJAACYDgkMAAAwHRIYAABgOiQwAADAdEhgAACA6ZDAAAAA0+ExagDwod5i0Zwhd3vKQFjFx0vz5l0sh7qeyZDAAIAPdZZ4/Sa3MNJhIFpZrdK0aeGrZzJcQgIAAKbDCAwA+BDrblDvyoOSpN2Z3ZlKAOHV0CB9+GFj+YYb/JtKIJB6JkMCAwA+2Orr9KfflkhiKgFEQE2NNHBgY9nfqQQCqWcyXEICAACmQwIDAABMhwQGAACYDgkMAAAwHRIYAABgOn4/hfTuu+/qqaeeUkVFhY4dO6ZXX31Vt912m2f7XXfdpWXLlnnVyc/P17p16zyfT506pfvvv19vvPGGYmNjVVhYqOeee05t2rQJ/EyAFqLro2siHQIAtHp+JzDV1dW69tprdc8992j06NFN7jN8+HAtWbLE89lms3ltHzt2rI4dO6aNGzeqrq5Od999tyZNmqQVK1b4Gw4AhEy9xaJnb/yxpwyEVXy8NGPGxXKo65mM3wnMiBEjNGLEiEvuY7PZlJWV1eS2jz/+WOvWrdOOHTvUv39/SdILL7ygkSNH6umnn1bHjh39DQlolZozkvPZ3IIwRBK96izxenbw2EiHgWhltUozZ4avnsmE5B6YzZs3KyMjQz169NDkyZN18uRJz7by8nKlpaV5khdJysvLU2xsrN5///0mj+dyueR0Or0WAAAQvYKewAwfPly//e1vtWnTJv3yl7/Uli1bNGLECDU0NEiSHA6HMjIyvOrExcUpPT1dDoejyWOWlZUpNTXVs3Tq1CnYYQPAV8QYbl114h+66sQ/FGO4Ix0Ooo3bLe3Z07i4/eh/gdYzmaBPJTBmzBhPuU+fPurbt6+6d++uzZs3a+jQoQEds7S0VCUlJZ7PTqeTJAZAyCXU1WrjS8WSLkwlkBDhiBBVzp2TevduLPszJUCg9Uwm5HMhXXnllWrfvr0OHDigoUOHKisrS8ePH/fap76+XqdOnfJ534zNZvvKjcAAuE8GQPQK+XtgPv/8c508eVLZ2dmSJLvdrtOnT6uiosKzz9tvvy23263c3NxQhwMAAFoBv0dgzpw5owMHDng+Hzp0SDt37lR6errS09M1a9YsFRYWKisrSwcPHtTDDz+sb37zm8rPz5ck9erVS8OHD9fEiRO1ePFi1dXVacqUKRozZgxPIAEAgGbxewTmgw8+0PXXX6/rr79eklRSUqLrr79eTzzxhCwWiz766CPdeuutuvrqqzVhwgT169dPf/7zn70uAS1fvlw9e/bU0KFDNXLkSA0ePFi/+c1vgndWAACgVfN7BGbIkCEyDMPn9vXr13/tMdLT03lpHQAACBhzIQEAANMJ+VNIAGBW9RaLfj1wtKcMhFV8vPTQQxfLoa5nMiQwAOBDnSVeZd+9J9JhIFpZrdJTT4WvnslwCQkAAJgOIzAA4EOM4dY3nCckSf9M6SAjhr/5EEZut3T4cGO5c2cptpn9L9B6JkMCAwA+JNTV6r3FEyQxlQAi4Nw5qVu3xrK/UwkEUs9kSGCAVo7pBnA56D9oqVrnuBIAAGjVSGAAAIDpkMAAAADTIYEBAACmw028AIKmOTd8AkAwkMAAgA8NsRb99voCTxkIq7g46b77LpZDXc9kWu+ZAcBlqo2L1xPDJkc6DEQrm01asCB89UyGe2AAAIDpMAIDAL4YhtLPOSVJpxJTpJiYCAeEqGIY0hdfNJbbt29+/wu0nsmQwACAD4l1Ln34wlhJTCWACDh7VsrIaCz7MyVAoPVMhktIAADAdEhgAACA6ZDAAAAA0+EeGADMOAzAdBiBAQAApkMCAwAATIdLSADgQ0OsRS/3HuopA2EVFycVFV0sh7qeybTeMwOAy1QbF6+HCqZGOgxEK5tNWro0fPVMhktIAADAdBiBAQBfDEOJdS5J0rl4W6t9JTtaKMNofKuuJCUl+TeVQCD1TIYRGADwIbHOpY/n366P59/uSWSAsDl7VmrTpnG5kJCEsp7JMAKDFo33kwAAmsIIDAAAMB0SGAAAYDokMAAAwHT8TmDeffddff/731fHjh0VExOj1157zWu7YRh64oknlJ2drcTEROXl5emTTz7x2ufUqVMaO3asUlJSlJaWpgkTJujMmTOXdSIAACB6+H0Tb3V1ta699lrdc889Gj169Fe2z5s3T88//7yWLVumbt266fHHH1d+fr727t2rhIQESdLYsWN17Ngxbdy4UXV1dbr77rs1adIkrVix4vLPCEBINOeGagAIF78TmBEjRmjEiBFNbjMMQ88++6ymT5+uUaNGSZJ++9vfKjMzU6+99prGjBmjjz/+WOvWrdOOHTvUv39/SdILL7ygkSNH6umnn1bHjh0v43QAIHjcsbFa0+NGTxkIK4tFuv32i+VQ1zOZoD5GfejQITkcDuXl5XnWpaamKjc3V+Xl5RozZozKy8uVlpbmSV4kKS8vT7GxsXr//ff1gx/8IJghAUDAXHFWFd9WGukwWjxedxAiCQnSqlXhq2cyQU1gHA6HJCkzM9NrfWZmpmebw+FQRkaGdxBxcUpPT/fs8+9cLpdcrosvkXI6ncEMGwBgEiRLuMAUL7IrKyvTrFmzIh0GWii+0AAg+gQ1gcnKypIkVVZWKjs727O+srJS1113nWef48ePe9Wrr6/XqVOnPPX/XWlpqUpKSjyfnU6nOnXqFMzQ0cpxAyoCkVhbo4/nN95L0GvqyzpnTYhwRIgq1dWN0wFI0pkzUnJyaOuZTFDvSuvWrZuysrK0adMmzzqn06n3339fdrtdkmS323X69GlVVFR49nn77bfldruVm5vb5HFtNptSUlK8FgAAEL38HoE5c+aMDhw44Pl86NAh7dy5U+np6ercubMeeOAB/fznP9dVV13leYy6Y8eOuu222yRJvXr10vDhwzVx4kQtXrxYdXV1mjJlisaMGcMTSAAAoFn8TmA++OADffe73/V8vnBpp6ioSEuXLtXDDz+s6upqTZo0SadPn9bgwYO1bt06zztgJGn58uWaMmWKhg4dqtjYWBUWFur5558PwukAAIBo4HcCM2TIEBmG4XN7TEyMZs+erdmzZ/vcJz09nZfWAQCAgPFmJgAAYDokMAAAwHRM8R4YAIgEd2ys3r6yv6cMhJXFIo0cebEc6nomQwIDAD644qy654czIx0GolVCgrQmgHdYBVrPZPiTAgAAmA4JDAAAMB0SGADwIbG2RnufKdTeZwqVWFsT6XAQbaqrG6cBSE5uLIe6nslwDwwAXEJSnSvSIbQKTLoaoLNnw1vPRBiBAQAApkMCAwAATIcEBgAAmA4JDAAAMB0SGAAAYDo8hQQAPrhjYrStU29PGQir2FjpO9+5WA51PZMhgQEAH1zxNo25c26kw0C0SkyUNm8OXz2Tab2pGQAAaLVIYAAAgOmQwACAD4m1Nap4/k5VPH8nUwkg/KqrpQ4dGhd/pxIIpJ7JcA8MAFxCu3POSIeAaPbFF+GtZyKMwAAAANNhBAYhwcRtAPzVnO+NcP4svqNaNhIYAACaQJLTsnEJCQAAmA4jMAAAhBAjOaFBAgMAPrhjYvS3rKs8ZSCsYmOl/v0vlkNdz2RIYADAB1e8TaOK5kc6DESrxERpx47w1TOZ1puaAQCAVosEBgAAmA6XkBAx4XznAxCIhLoavfXf90mS8n66UDXxCRGOCFHl7FkpJ6exvHevlJQU2nomQwIDv5F4IFrEGNIVzuOeMhBWhiH94x8Xy6GuZzJcQgIAAKZDAgMAAEyHBAYAAJhO0BOYmTNnKiYmxmvp2bOnZ3tNTY2Ki4vVrl07tWnTRoWFhaqsrAx2GAAAoBULyQjMNddco2PHjnmW9957z7Nt6tSpeuONN7Rq1Spt2bJFR48e1ejRo0MRBgAAaKVC8hRSXFycsrKyvrK+qqpKL774olasWKGbb75ZkrRkyRL16tVL27Zt06BBg0IRDgAExIiR/t6us6cMhFVMzMXHof2ZyiLQeiYTkgTmk08+UceOHZWQkCC73a6ysjJ17txZFRUVqqurU15enmffnj17qnPnziovL/eZwLhcLrlcLs9np9MZirABwEtNfIKG/XRhpMNAtEpKkvbsCV89kwl6ApObm6ulS5eqR48eOnbsmGbNmqVvf/vb2r17txwOh6xWq9LS0rzqZGZmyuFw+DxmWVmZZs2aFexQAQCtUDjfVcV7sSIn6AnMiBEjPOW+ffsqNzdXXbp00R/+8AclJiYGdMzS0lKVlJR4PjudTnXq1OmyYwUAoCVoTiL02dyCMERiHiF/jDotLU1XX321Dhw4oKysLNXW1ur06dNe+1RWVjZ5z8wFNptNKSkpXgsAhFpCXY02/Pd92vDf9ymhribS4SDKBNz/zp6VrrmmcTl7NnQBRljIE5gzZ87o4MGDys7OVr9+/RQfH69NmzZ5tu/fv1+HDx+W3W4PdSgA4JcYQ7r65GFdffIwUwkg7ALuf4bROAfS3r2teiqBoF9Ceuihh/T9739fXbp00dGjRzVjxgxZLBb9+Mc/VmpqqiZMmKCSkhKlp6crJSVF999/v+x2O08gAQCAZgt6AvP555/rxz/+sU6ePKkOHTpo8ODB2rZtmzp06CBJmj9/vmJjY1VYWCiXy6X8/HwtXMhd/gAAoPmCnsCsXLnyktsTEhK0YMECLViwINg/GgAARAnmQgIAAKZDAgMAAEwnJG/iBYDWwIiRPk/J8JSBcAq4/8XESF26XCy3UiQwAOBDTXyCBk9+KdJhIEoF3P+SkqTPPgt6PC0Nl5AAAIDpMAIDL8zrAQAwA0ZgAMAHW51Lry+bqteXTZWtzhXpcBBlAu5/585JAwY0LufOhS7ACGMEBgB8iDUMXev4xFMGwing/ud2Sx98cLHcSjECAwAATIcEBgAAmA4JDAAAMB0SGAAAYDokMAAAwHR4CgkALuFkYkqkQ0AUC7j/tW8f3EBaIBIYAPDhnDVB/f6/FZEOA1Eq4P6XnCydOBH8gFoYLiEBAADTYQQGAAATaO5UL5/NLQhxJC0DIzAA4IOtzqWVKx7VyhWPMpUAwi7g/nfunDRkSOPCVAJo6ZqTmUdLVg4ES6xhaNCR3Z4yEE4B9z+3W9qy5WK5lWIEBgAAmA4JDAAAMB0SGAAAYDrcAwMAQJRpDfdNksBEkeY+ggcAQEtHAgMAl3A23hbpEBDFAu5/SUnBDaQFIoEBAB/OWROUU/LHSIeBKBVw/0tOlqqrgx9QC8NNvAAAwHRIYAAAgOlwCckEuPkWiAxbfa0WvTpHkjT5B4/JFWeNcESIJgH3v5oaqbCwsfzHP0oJCSGKMLJIYADAh1i3Wzd/+oGnDIRTwP2voUFau/ZiuZUigQEAoBXp9fg6ffwv5XPW1jkCwz0wAADAdCKawCxYsEBdu3ZVQkKCcnNztX379kiGAwAATCJil5B+//vfq6SkRIsXL1Zubq6effZZ5efna//+/crIyIhUWM3WGl7DDACAWUUsgXnmmWc0ceJE3X333ZKkxYsXa82aNXrppZf06KOPRiosScF76ockBwCA0IhIAlNbW6uKigqVlpZ61sXGxiovL0/l5eVf2d/lcsnlcnk+V1VVSZKcTmdI4nO7zobkuE3pPHVV2H4WAP801NbowrdMg+us3AZPIiF8Au1/weq3ofgde+GYhmFc/sGMCPjnP/9pSDK2bt3qtX7atGnGwIEDv7L/jBkzDEksLCwsLCwsrWA5cuTIZecSpniMurS0VCUlJZ7Pbrdbp06dUrt27RQTExPByMLH6XSqU6dOOnLkiFJSUiIdTqtBuwYfbRoatGto0K6h4atdDcPQl19+qY4dO172z4hIAtO+fXtZLBZVVlZ6ra+srFRWVtZX9rfZbLLZvGfkTEtLC2WILVZKSgr/yUKAdg0+2jQ0aNfQoF1Do6l2TU1NDcqxI/IYtdVqVb9+/bRp0ybPOrfbrU2bNslut0ciJAAAYCIRu4RUUlKioqIi9e/fXwMHDtSzzz6r6upqz1NJAAAAvkQsgbnjjjt04sQJPfHEE3I4HLruuuu0bt06ZWZmRiqkFs1ms2nGjBlfuZSGy0O7Bh9tGhq0a2jQrqERjnaNMYxgPMsEAAAQPsyFBAAATIcEBgAAmA4JDAAAMB0SGAAAYDokMBGyYMECde3aVQkJCcrNzdX27dsvuf/p06dVXFys7Oxs2Ww2XX311Vq7dq1n+8yZMxUTE+O19OzZM9Sn0eL4065Dhgz5SpvFxMSooODiBJuGYeiJJ55Qdna2EhMTlZeXp08++SQcp9KiBLtd77rrrq9sHz58eDhOpUXx93vg2WefVY8ePZSYmKhOnTpp6tSpqqmpuaxjtjbBblO+Wxv50651dXWaPXu2unfvroSEBF177bVat27dZR2zSZc9GQH8tnLlSsNqtRovvfSSsWfPHmPixIlGWlqaUVlZ2eT+LpfL6N+/vzFy5EjjvffeMw4dOmRs3rzZ2Llzp2efGTNmGNdcc41x7Ngxz3LixIlwnVKL4G+7njx50qu9du/ebVgsFmPJkiWefebOnWukpqYar732mvG3v/3NuPXWW41u3boZ586dC9NZRV4o2rWoqMgYPny4136nTp0K0xm1DP626/Llyw2bzWYsX77cOHTokLF+/XojOzvbmDp1asDHbG1C0aZ8t/rfrg8//LDRsWNHY82aNcbBgweNhQsXGgkJCcaHH34Y8DGbQgITAQMHDjSKi4s9nxsaGoyOHTsaZWVlTe6/aNEi48orrzRqa2t9HnPGjBnGtddeG+xQTcXfdv138+fPN9q2bWucOXPGMAzDcLvdRlZWlvHUU0959jl9+rRhs9mM3/3ud8ENvgULdrsaRmMCM2rUqGCHair+tmtxcbFx8803e60rKSkxbrzxxoCP2dqEok35bvW/XbOzs41f/epXXutGjx5tjB07NuBjNoVLSGFWW1uriooK5eXledbFxsYqLy9P5eXlTdb505/+JLvdruLiYmVmZqp3796aM2eOGhoavPb75JNP1LFjR1155ZUaO3asDh8+HNJzaUkCadd/9+KLL2rMmDFKTk6WJB06dEgOh8PrmKmpqcrNzW32Mc0uFO16webNm5WRkaEePXpo8uTJOnnyZFBjb8kCaddvfetbqqio8Ayzf/rpp1q7dq1GjhwZ8DFbk1C06QV8t/rXri6XSwkJCV7rEhMT9d577wV8zKaYYjbq1uSLL75QQ0PDV944nJmZqX379jVZ59NPP9Xbb7+tsWPHau3atTpw4IDuu+8+1dXVacaMGZKk3NxcLV26VD169NCxY8c0a9Ysffvb39bu3bvVtm3bkJ9XpAXSrv9q+/bt2r17t1588UXPOofD4TnGvx/zwrbWLhTtKknDhw/X6NGj1a1bNx08eFCPPfaYRowYofLyclkslqCeQ0sUSLveeeed+uKLLzR48GAZhqH6+nrde++9euyxxwI+ZmsSijaV+G4NpF3z8/P1zDPP6KabblL37t21adMmvfLKK54/uoPVVxmBMQG3262MjAz95je/Ub9+/XTHHXfoP//zP7V48WLPPiNGjNAPf/hD9e3bV/n5+Vq7dq1Onz6tP/zhDxGM3DxefPFF9enTRwMHDox0KK2Kr3YdM2aMbr31VvXp00e33XabVq9erR07dmjz5s2RCdQENm/erDlz5mjhwoX68MMP9corr2jNmjV68sknIx2aaTWnTflu9d9zzz2nq666Sj179pTVatWUKVN09913KzY2uCkHCUyYtW/fXhaLRZWVlV7rKysrlZWV1WSd7OxsXX311V5/mfbq1UsOh0O1tbVN1klLS9PVV1+tAwcOBC/4FiyQdr2gurpaK1eu1IQJE7zWX6gXyDFbi1C0a1OuvPJKtW/fnv56iXZ9/PHHNW7cOP30pz9Vnz599IMf/EBz5sxRWVmZ3G73Zf1btQahaNOm8N3a6FLt2qFDB7322muqrq7WP/7xD+3bt09t2rTRlVdeGfAxm0ICE2ZWq1X9+vXTpk2bPOvcbrc2bdoku93eZJ0bb7xRBw4c8PoP9fe//13Z2dmyWq1N1jlz5owOHjyo7Ozs4J5ACxVIu16watUquVwu/eQnP/Fa361bN2VlZXkd0+l06v333//aY7YWoWjXpnz++ec6efIk/fUS7Xr27Nmv/AV74Y8awzAu69+qNQhFmzaF79bm96uEhAR94xvfUH19vf74xz9q1KhRl31ML82+3RdBs3LlSsNmsxlLly419u7da0yaNMlIS0szHA6HYRiGMW7cOOPRRx/17H/48GGjbdu2xpQpU4z9+/cbq1evNjIyMoyf//znnn0efPBBY/PmzcahQ4eMv/zlL0ZeXp7Rvn174/jx42E/v0jxt10vGDx4sHHHHXc0ecy5c+caaWlpxuuvv2589NFHxqhRo6LyMepgtuuXX35pPPTQQ0Z5eblx6NAh46233jJuuOEG46qrrjJqampCfj4thb/tOmPGDKNt27bG7373O+PTTz81NmzYYHTv3t340Y9+1OxjtnahaFO+W/1v123bthl//OMfjYMHDxrvvvuucfPNNxvdunUz/u///q/Zx2wOEpgIeeGFF4zOnTsbVqvVGDhwoLFt2zbPtu985ztGUVGR1/5bt241cnNzDZvNZlx55ZXGL37xC6O+vt6z/Y477jCys7MNq9VqfOMb3zDuuOMO48CBA+E6nRbD33bdt2+fIcnYsGFDk8dzu93G448/bmRmZho2m80YOnSosX///lCeQosUzHY9e/asMWzYMKNDhw5GfHy80aVLF2PixIlR80v2X/nTrnV1dcbMmTON7t27GwkJCUanTp2M++67z+uXwtcdMxoEu035bm3kT7tu3rzZ6NWrl2Gz2Yx27doZ48aNM/75z3/6dczmiDEMH+NkAAAALRT3wAAAANMhgQEAAKZDAgMAAEyHBAYAAJgOCQwAADAdEhgAAGA6JDAAAMB0SGAAAIDpkMAAAADTIYEBAACmQwIDAABMhwQGAACYzv8PhHox3IxZ/VoAAAAASUVORK5CYII=",
      "text/plain": [
       "<Figure size 640x480 with 1 Axes>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "scores = get_all_scores(rst_all_citation)\n",
    "score_threshold = np.percentile(scores, 99.7)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 12,
   "id": "44d65015",
   "metadata": {},
   "outputs": [],
   "source": [
    "rst_abstract = rearrange_query_results(rst_all_abstract, score_threshold=0.8)\n",
    "rst_citation = rearrange_query_results(rst_all_citation, score_threshold=score_threshold)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 13,
   "id": "8b3ef506",
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "1. lorimer et al. 2007, arXiv:0709.4301\n",
      "\n",
      "Reasons of citation: referenced as the original discovery of fast radio bursts (frbs).; referenced as the first discovery of fast radio bursts (frbs).; referenced as the first discovery of fast radio bursts (frbs).; referenced as the first discovery of fast radio bursts (frbs).; referenced to provide background information on the duration of fast radio bursts (frbs).; referenced to indicate the year of discovery of fast radio bursts (frbs).; referenced to introduce the concept of fast radio bursts (frbs) and their characteristics.; referenced to indicate the discovery of fast radio bursts (frbs) and establish the intense study of frbs in radio astronomy.; referenced to provide a general description of fast radio bursts (frbs) and their characteristics.; referenced to highlight the discovery of fast radio bursts (frbs) through radio astronomy.\n",
      "\n",
      "2. thornton et al. 2013, arXiv:1307.1628\n",
      "\n",
      "Reasons of citation: referenced to provide background information on the duration of fast radio bursts (frbs).; referenced to provide additional support and background information on fast radio bursts (frbs).; referenced to provide additional evidence for the existence of fast radio bursts (frbs).; referenced to indicate the discovery of fast radio bursts (frbs) and establish the intense study of frbs in radio astronomy.\n",
      "\n",
      "3. petroff et al. 2022; petroff et al. 2021, arXiv:2107.10113\n",
      "\n",
      "Reasons of citation: referenced as a recent review of fast radio bursts (frbs) to provide background information on the topic.; these references are cited to indicate that fast radio bursts (frbs) are a new class of radio transient that can be observed from distant galaxies, and to acknowledge that the source of frbs is still unknown.\n",
      "\n",
      "4. petroff et al. 2016, arXiv:1601.03547\n",
      "\n",
      "Reasons of citation: referenced to provide background information on the duration of fast radio bursts (frbs).\n",
      "\n",
      "5. cordes & chatterjee 2019, arXiv:1906.05878\n",
      "\n",
      "Reasons of citation: these references are cited to indicate that fast radio bursts (frbs) are a new class of radio transient that can be observed from distant galaxies, and to acknowledge that the source of frbs is still unknown.\n",
      "\n"
     ]
    }
   ],
   "source": [
    "# citation search\n",
    "print_citation_query_results(citation_data, rst_citation)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 15,
   "id": "cb4f4d51",
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "1. Bo-Jun Wang et al., arXiv:2301.01429\n",
      "\n",
      "Title: Atlas of dynamic spectra of fast radio burst FRB 20201124A\n",
      "\n",
      "Abstract: Fast radio bursts (FRBs) are highly dispersed millisecond-duration radio bursts, of which the physical origin is still not fully understood. FRB 20201124A is one of the most actively repeating FRBs. In this paper, we present the collection of 1863 burst dynamic spectra of FRB 20201124A measured with the Five-hundred-meter Aperture Spherical radio Telescope (FAST). The current collection, taken from the observation during the FRB active phase from April to June 2021, is the largest burst sample detected in any FRB so far. The standard PSRFITs format is adopted, including dynamic spectra of the burst, and the time information of the dynamic spectra, in addition, mask files help readers to identify the pulse positions are also provided.\n",
      "\n",
      "2. Matthew Bailes, , arXiv:2211.06048\n",
      "\n",
      "Title: The discovery and scientific potential of fast radio bursts\n",
      "\n",
      "Abstract: Fast radio bursts (FRBs) are millisecond-timescale bursts of coherent radio emission that are luminous enough to be detectable at cosmological distances. In this review I describe the discovery of FRBs, subsequent advances in our understanding of them, and future prospects. Thousands of potentially observable FRBs reach Earth every day; they probably originate from highly magnetic and/or rapidly rotating neutron stars in the distant Universe. Some FRBs repeat, with this sub-class often occurring in highly magnetic environments. Two repeaters exhibit cyclic activity windows, consistent with an orbital period. One nearby FRB was from a Galactic magnetar during an X-ray outburst. The host galaxies of some FRBs have been located, providing information about the host environments and the total baryonic content of the Universe.\n",
      "\n",
      "3. Griffin Foster et al., arXiv:1808.07809\n",
      "\n",
      "Title: Verifying and Reporting Fast Radio Bursts\n",
      "\n",
      "Abstract: Fast Radio Bursts (FRBs) are a class of short-duration transients at radio wavelengths with inferred astrophysical origin. The prototypical FRB is a broadband signal that occurs over the extent of the receiver frequency range, is narrow in time, and is highly dispersed, following a $\\nu^{-2}$ relation. However, some FRBs appear band-limited, and show apparent scintillation, complex frequency-dependent structure, or multi-component pulse shapes. While there is sufficient evidence that FRBs are indeed astrophysical, their one-off nature necessitates extra scrutiny when reporting a detection as bona fide and not a false positive. Currently, there is no formal validation framework for FRBs, rather a set of community practices. In this article, we discuss potential sources of false positives, and suggest a framework in which FRB-like events can be evaluated as real or otherwise. We present examples of false-positive events in data from the Arecibo, LOFAR, and Nanshan telescopes, which while FRB-like, are found to be due to instrumental variations, noise, and radio-frequency interference. Differentiating these false-positive detections from astrophysical events requires knowledge and tests beyond thresholded single-pulse detection. We discuss post-detection analyses, verification tests, and datasets which should be provided when reporting an FRB detection.\n",
      "\n",
      "4. J. I. Katz, , arXiv:1604.01799\n",
      "\n",
      "Title: Fast Radio Bursts---A Brief Review: Some Questions, Fewer Answers\n",
      "\n",
      "Abstract: Fast Radio Bursts are millisecond bursts of radio radiation at frequencies of about 1 GHz, recently discovered in pulsar surveys. They have not yet been definitively identified with any other astronomical object or phenomenon. The bursts are strongly dispersed, indicating passage through a high column density of low density plasma. The most economical interpretation is that this is the interglactic medium, indicating that FRB are at \"cosmological\" distances with redshifts in the range 0.3--1.3. Their inferred brightness temperatures are as high as $10^{37\\,\\circ}$K, implying coherent emission by \"bunched\" charges, as in radio pulsars. I review the astronomical sites, objects and emission processes that have been proposed as the origin of FRB, with particular attention to Soft Gamma Repeaters and giant pulsar pulses.\n",
      "\n",
      "5. E. Petroff et al., arXiv:1810.10773\n",
      "\n",
      "Title: A fast radio burst with a low dispersion measure\n",
      "\n",
      "Abstract: Fast radio bursts (FRBs) are millisecond pulses of radio emission of seemingly extragalactic origin. More than 50 FRBs have now been detected, with only one seen to repeat. Here we present a new FRB discovery, FRB 110214, which was detected in the high latitude portion of the High Time Resolution Universe South survey at the Parkes telescope. FRB 110214 has one of the lowest dispersion measures of any known FRB (DM = 168.9$\\pm$0.5 pc cm$^{-3}$), and was detected in two beams of the Parkes multi-beam receiver. A triangulation of the burst origin on the sky identified three possible regions in the beam pattern where it may have originated, all in sidelobes of the primary detection beam. Depending on the true location of the burst the intrinsic fluence is estimated to fall in the range of 50 -- 2000 Jy ms, making FRB 110214 one of the highest-fluence FRBs detected with the Parkes telescope. No repeating pulses were seen in almost 100 hours of follow-up observations with the Parkes telescope down to a limiting fluence of 0.3 Jy ms for a 2-ms pulse. Similar low-DM, ultra-bright FRBs may be detected in telescope sidelobes in the future, making careful modeling of multi-beam instrument beam patterns of utmost importance for upcoming FRB surveys.\n",
      "\n"
     ]
    }
   ],
   "source": [
    "# abstract search\n",
    "print_abstract_query_results(abstract_data, rst_abstract)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "id": "503a6f6d",
   "metadata": {},
   "outputs": [],
   "source": []
  },
  {
   "cell_type": "markdown",
   "id": "c2acc4a0",
   "metadata": {},
   "source": [
    "### example where citation is a bit better, but still a lot of noise\n",
    "I know where the petroff papers come from... because I cited those papers in my own paper.  I can imagine why gpt thinks the reason for citation is like that, and it's unavoidable.\n",
    "\n",
    "Perhaps in the future need to let gpt determine what the \"noise\" reasons of citation in each entry.\n",
    "\n",
    "In this case the abstract search is not too bad, and I'm glad that it returned my own paper..."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 22,
   "id": "f3764587",
   "metadata": {},
   "outputs": [],
   "source": [
    "q = 'fast radio burst as a probe of circumgalactic medium (cgm)'"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 23,
   "id": "e8d6c294",
   "metadata": {},
   "outputs": [],
   "source": [
    "rst_all_abstract = retriever_abstract.retrieve(q)\n",
    "rst_all_citation = retriever_citation.retrieve(q)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 24,
   "id": "b6ae7120",
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAigAAAGdCAYAAAA44ojeAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjcuMywgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/OQEPoAAAACXBIWXMAAA9hAAAPYQGoP6dpAAAtlUlEQVR4nO3de3xU9Z3/8ffMJDO5wCTLJQnUcNGqEAVdUWBWaimmRIwrlth6YTEq1RUDW8mKgg8FxK1Q7E/Ulku3RbC7Iiv1toJFEAuuEgTjthVQVhEaLE5QXDISkskkc35/ZDNpAhSGZGa+Z/J6Ph7n8fjOmZPz+Z4zkbz9nnO+47AsyxIAAIBBnInuAAAAQHsEFAAAYBwCCgAAMA4BBQAAGIeAAgAAjENAAQAAxiGgAAAA4xBQAACAcVIS3YEzEQ6HdfDgQXXv3l0OhyPR3QEAAKfBsix9/fXX6tu3r5zOvz5GYsuAcvDgQeXn5ye6GwAA4AwcOHBAZ5111l/dxpYBpXv37pKaD9Dr9Sa4NwDiorZW6tu3uX3woJSZmdj+oFm8Phc+/6QQCASUn58f+Tv+19gyoLRc1vF6vQQUoKtwuVrbXi9/oEwRr8+Fzz+pnM7tGdwkCwAAjGPLERQAXVBKilRa2tqGGeL1ufD5dzkOy7KsRHciWoFAQFlZWaqpqeESDwAANhHN328u8QAAAOMwTgbAHixLOnasuZ2RITEHkhni9bnw+Xc5jKAAsIdjx6Ru3ZqXlj9USLx4fS58/l0OAQUAABiHgAIAAIxDQAEAAMYhoAAAAOMQUAAAgHEIKAAAwDjMgwLAHlwu6frrW9swQ7w+Fz7/Loep7gEAQFxE8/ebERQAxhkwc90pt9m/oDgOPQGQKNyDAgAAjENAAWAL6Q312v+Ta7T/J9covaE+0d1Bi9ra5u/FcTia23avA2MQUAAAgHEIKAAAwDgEFAAAYBwCCgAAMA4BBQAAGIeAAgAAjMNEbQBsIex06s2zL420YQiXS7r66ta23evAGAQUALYQTHHr9u/PTXQ30F5amrTu1DP/2qYOjMH/hgAAAOMQUAAAgHGiCihz586Vw+FoswwaNCjyfn19vcrKytSzZ09169ZNJSUlqq6ubrOPqqoqFRcXKyMjQzk5OZoxY4YaGxs752gAJK30hnrtfrxEux8vYap7k9TWSpmZzUusp7qPRx0YI+p7UC644AK98cYbrTtIad3F9OnTtW7dOq1Zs0ZZWVmaOnWqJkyYoHfeeUeS1NTUpOLiYuXl5Wnr1q36/PPPdcsttyg1NVWPPvpoJxwOgGSWEQomugs4kWPHkqsOjBB1QElJSVFeXt5x62tqarR8+XKtWrVKY8aMkSStWLFCgwcP1rZt2zRy5Eht2LBBu3fv1htvvKHc3FxdfPHFeuSRR3T//fdr7ty5crvdHT8iAABge1Hfg/Lxxx+rb9++OvvsszVx4kRVVVVJkiorKxUKhVRYWBjZdtCgQerXr58qKiokSRUVFRoyZIhyc3Mj2xQVFSkQCGjXrl0nrRkMBhUIBNosAAAgeUUVUEaMGKGVK1dq/fr1Wrp0qfbt26dvfetb+vrrr+X3++V2u5Wdnd3mZ3Jzc+X3+yVJfr+/TThpeb/lvZOZP3++srKyIkt+fn403QYAADYT1SWecePGRdpDhw7ViBEj1L9/fz3//PNKT0/v9M61mDVrlsrLyyOvA4EAIQUAgCTWoceMs7Ozdd555+mTTz5RXl6eGhoadOTIkTbbVFdXR+5ZycvLO+6pnpbXJ7qvpYXH45HX622zAACA5NWhgHL06FHt3btXffr00bBhw5SamqpNmzZF3t+zZ4+qqqrk8/kkST6fTx988IEOHToU2Wbjxo3yer0qKCjoSFcAJLmww6Ft+RdqW/6FCjscie4OWjid0re/3bzE8isI4lUHxojqEs+9996rv//7v1f//v118OBBzZkzRy6XSzfddJOysrI0efJklZeXq0ePHvJ6vZo2bZp8Pp9GjhwpSRo7dqwKCgo0adIkLVy4UH6/Xw8++KDKysrk8XhicoAAkkMw1aMbb16Q6G6gvfR0afPm5KkDY0QVUD777DPddNNNOnz4sHr37q1Ro0Zp27Zt6t27tyRp0aJFcjqdKikpUTAYVFFRkZYsWRL5eZfLpbVr12rKlCny+XzKzMxUaWmp5s2b17lHBQAAbM1hWZaV6E5EKxAIKCsrSzU1NdyPAiShATNP/aVw+xcUx6EnADpTNH+/uZAHwBbSG+pV+dTNqnzqZqa6N0ltrdS7d/MS66nu41EHxoh6JlkASJSedUzSaKQvv0yuOjACIygAAMA4BBQAAGAcAgoAADAOAQUAABiHgAIAAIzDUzwAbCHscOgPeedG2jCE0yldemlr2+51YAwCCgBbCKZ6NL50UaK7gfbS06UdO5KnDoxBDAUAAMYhoAAAAOMQUADYQlqoXm8vvV1vL71daSGmujfGsWPSgAHNy7Fj9q8DY3APCgBbcFjSWYFDkTYMYVnSn/7U2rZ7HRiDERQAAGAcAgoAADAOAQUAABiHgAIAAIxDQAEAAMbhKR4AtmA5pP/p2S/ShiEcDqmgoLVt9zowBgEFgC3Up6Zp7A+XJLobaC8jQ9q1K3nqwBhc4gEAAMZhBAVApxkwc90pt9m/oDgOPQFgd4ygALCFtFC9Nvzqbm341d1MdW+SY8ekCy5oXmI91X086sAYjKAAsAWHJZ13uCrShiEsS9q9u7Vt9zowBiMoAADAOAQUAABgHAIKAAAwDgEFAAAYh5tkASQtHnsG7IuAAsAWLIf0mTcn0oYhHA6pf//Wtt3rwBgEFAC2UJ+aplFTnk50N9BeRoa0f3/y1IExuAcFAAAYh4ACAACMwyUeALbgCQX1/KqZkqQf3LzgtG6ARRzU1UlXXNHcfustKT3d3nVgDAIKAFtwWpYu8n8cacMQ4bD03nutbbvXgTG4xAMAAIxDQAEAAMYhoAAAAOMQUAAAgHEIKAAAwDg8xQMgrjryePDhdG8n9gSdplev5KoDIxBQANhCnTtNw/5pVaK7gfYyM6UvvkieOjAGl3gAAIBxGEEB0KWdziWn/QuK49ATAH+JERQAtuAJBbV61UytXjVTnlAw0d1Bi7o6afTo5qWuzv51YAxGUADYgtOyNPLAzkgbhgiHpS1bWtt2rwNjMIICAACMQ0ABAADGIaAAAADjEFAAAIBxCCgAAMA4PMUDwDaOpXoS3QWcSEZGctWBEQgoAGyhzp2mgvIXEt0NtJeZKdXWJk8dGINLPAAAwDgEFAAAYBwCCgBb8DQ26Ok1c/X0mrnyNDYkujtoUV8vFRc3L/X19q8DY3QooCxYsEAOh0P33HNPZF19fb3KysrUs2dPdevWTSUlJaqurm7zc1VVVSouLlZGRoZycnI0Y8YMNTY2dqQrAJKcMxzWmE/f05hP35OTqc7N0dQkvfZa89LUZP86MMYZB5QdO3boF7/4hYYOHdpm/fTp0/Xqq69qzZo12rJliw4ePKgJEyZE3m9qalJxcbEaGhq0detWPfPMM1q5cqVmz5595kcBAACSyhkFlKNHj2rixIn65S9/qb/5m7+JrK+pqdHy5cv1+OOPa8yYMRo2bJhWrFihrVu3atu2bZKkDRs2aPfu3fr3f/93XXzxxRo3bpweeeQRLV68WA0NDNsCAIAzDChlZWUqLi5WYWFhm/WVlZUKhUJt1g8aNEj9+vVTRUWFJKmiokJDhgxRbm5uZJuioiIFAgHt2rXrhPWCwaACgUCbBQAAJK+o50FZvXq13n//fe3YseO49/x+v9xut7Kzs9usz83Nld/vj2zzl+Gk5f2W905k/vz5evjhh6PtKgAAsKmoRlAOHDigH/3oR3r22WeVlpYWqz4dZ9asWaqpqYksBw4ciFttAAAQf1EFlMrKSh06dEiXXHKJUlJSlJKSoi1btuipp55SSkqKcnNz1dDQoCNHjrT5uerqauXl5UmS8vLyjnuqp+V1yzbteTweeb3eNgsAAEheUV3iufLKK/XBBx+0WXfbbbdp0KBBuv/++5Wfn6/U1FRt2rRJJSUlkqQ9e/aoqqpKPp9PkuTz+fTjH/9Yhw4dUk5OjiRp48aN8nq9Kigo6IxjApCE6txpGnD/2kR3A+1lZkqWlTx1YIyoAkr37t114YUXtlmXmZmpnj17RtZPnjxZ5eXl6tGjh7xer6ZNmyafz6eRI0dKksaOHauCggJNmjRJCxculN/v14MPPqiysjJ5PHwRGAAAiMGXBS5atEhOp1MlJSUKBoMqKirSkiVLIu+7XC6tXbtWU6ZMkc/nU2ZmpkpLSzVv3rzO7goAALAph2XZb8wsEAgoKytLNTU13I8CGGTAzHUx27ensUGPr/1/kqTya/5ZwRR3zGq1t39Bcdxq2U59vTRpUnP73/5NitUDFPGqg5iK5u8338UDwBac4bCK97yj4j3vMNW9SZqapN/8pnmJ9VT38agDYxBQAACAcTr9HhQASDanc+mKy0BA52IEBQAAGIeAAgAAjENAAQAAxiGgAAAA43CTLABbqEv1aPD030TaMERGhnT0aGvb7nVgDAIKAHtwOFTnZnIu4zgczd+Tkyx1YAwu8QAAAOMQUADYgrsxpJ+uW6Sfrlskd2Mo0d1Bi2BQuvXW5iUYtH8dGIOAAsAWXOEmXb9zk67fuUmuMFOdG6OxUXrmmealsdH+dWAMAgoAADAOAQUAABiHgAIAAIxDQAEAAMYhoAAAAOMQUAAAgHGYSRaALdSlenTJtGcjbRgiI0M6dKi1bfc6MAYBBYA9OBz6KiMr0b1Aew6H1Lt38tSBMbjEAwAAjENAAWAL7saQ5m1YqnkbljLVvUmCQamsrHmJ9VT38agDYxBQANiCK9ykW/57nW7573VMdW+SxkZpyZLmJdZT3cejDoxBQAEAAMYhoAAAAOMQUAAAgHEIKAAAwDgEFAAAYBwCCgAAMA4zyQKwhfpUt0bdtTzShiHS06V9+1rbdq8DYxBQANiC5XDqs6zcRHcD7Tmd0oAByVMHxiCgADgtA2auS3QXAHQh3IMCwBZSm0Ka9bunNet3Tyu1ianujdHQIM2Y0bw0NNi/DoxBQAFgCylNTfrH7S/qH7e/qJQmpro3Rigk/fSnzUsohsExXnVgDAIKAAAwDgEFAAAYh4ACAACMQ0ABAADGIaAAAADjEFAAAIBxmKgNgC3Up7r13dsXR9owRHq6tHNna9vudWAMAgoAW7AcTn3cu3+iu4H2nE7pgguSpw6MwSUeAABgHEZQANhCalNIZRXPS5IW+36gkCs1wT2CpOZp5x99tLn9wAOSO0aX3+JVB8ZwWJZlJboT0QoEAsrKylJNTY28Xm+iuwN0CYn+ssD0hnp9uOh6SdLg6b9RnTstof1pb/+C4kR3ITFqa6Vu3ZrbR49KmZn2roOYiubvN5d4AACAcQgoAADAOAQUAABgHAIKAAAwDgEFAAAYh4ACAACMwzwoAGwhmJKqa295PNKGIdLSpO3bW9t2rwNjEFAA2ELY6dIf+5yX6G6gPZdLuuyy5KkDY3CJBwAAGIcRFAC2kNoU0m3v/ackacWl1zLVvSkaGqQnn2xu/+hHsZ3qPh51YAwCCgBbSGlq0gObV0iS/u1viwkopgiFpPvua27ffXfsgkO86sAYUV3iWbp0qYYOHSqv1yuv1yufz6ff/va3kffr6+tVVlamnj17qlu3biopKVF1dXWbfVRVVam4uFgZGRnKycnRjBkz1NjY2DlHAwAAkkJUAeWss87SggULVFlZqffee09jxozR+PHjtWvXLknS9OnT9eqrr2rNmjXasmWLDh48qAkTJkR+vqmpScXFxWpoaNDWrVv1zDPPaOXKlZo9e3bnHhUAALC1Dn+bcY8ePfTYY4/p+uuvV+/evbVq1Spdf33zN45+9NFHGjx4sCoqKjRy5Ej99re/1TXXXKODBw8qNzdXkrRs2TLdf//9+uKLL+Q+zSE7vs0Y6FyJ/qbi08G3GRuKbzNGFOLybcZNTU1avXq1amtr5fP5VFlZqVAopMLCwsg2gwYNUr9+/VRRUSFJqqio0JAhQyLhRJKKiooUCAQiozAnEgwGFQgE2iwAACB5RR1QPvjgA3Xr1k0ej0d33XWXXnrpJRUUFMjv98vtdis7O7vN9rm5ufL7/ZIkv9/fJpy0vN/y3snMnz9fWVlZkSU/Pz/abgMAABuJOqCcf/75+v3vf693331XU6ZMUWlpqXbv3h2LvkXMmjVLNTU1keXAgQMxrQcAABIr6seM3W63vvnNb0qShg0bph07dujJJ5/UDTfcoIaGBh05cqTNKEp1dbXy8vIkSXl5edreMlXxX7zf8t7JeDweeTyeaLsKIIkEU1J1402PRtowRFqa9LvftbbtXgfG6PA8KOFwWMFgUMOGDVNqaqo2bdqkkpISSdKePXtUVVUln88nSfL5fPrxj3+sQ4cOKScnR5K0ceNGeb1eFRQUdLQrAE7ADjfAno6w06Vt/YYmuhtoz+WSRo9OnjowRlQBZdasWRo3bpz69eunr7/+WqtWrdLmzZv1+uuvKysrS5MnT1Z5ebl69Oghr9eradOmyefzaeTIkZKksWPHqqCgQJMmTdLChQvl9/v14IMPqqysjBESAAAQEVVAOXTokG655RZ9/vnnysrK0tChQ/X666/ru9/9riRp0aJFcjqdKikpUTAYVFFRkZYsWRL5eZfLpbVr12rKlCny+XzKzMxUaWmp5s2b17lHBSDppDQ16qY/rJckPXfRVWp0MRG2EUIh6V//tbl9551Saowuv8WrDozR4XlQEoF5UIDTlyyXeJgHxVDMg4IoxGUeFAAAgFghoAAAAOMQUAAAgHEIKAAAwDgEFAAAYBwCCgAAMA4TCQCwhYaUVN12/ZxIG4bweKS1a1vbdq8DYxBQANhCk9Ol351zWaK7gfZSUqTiOMwBE686MAaXeAAAgHEYQQFgCylNjbpu92ZJ0ssFo5nq3hShkPTss83tiRNjO9V9POrAGEx1DyQ5pro3R1JOh89U94gCU90DAABbI6AAAADjEFAAAIBxCCgAAMA4BBQAAGAcAgoAADAOEwkAsIWGlFTdPX5mpA1DeDzS88+3tu1eB8YgoACwhSanS68NGpXobqC9lBTp+99PnjowBpd4AACAcRhBAWALrnCTiv6nQpL0+nk+NTldCe4RJEmNjdJLLzW3v/e95pEOO9eBMfiEAdiCuzGkJa8skNQy1T0BxQjBoPSDHzS3jx6NXXCIVx0Yg0s8AADAOAQUAABgHAIKAAAwDgEFAAAYh4ACAACMQ0ABAADG4TktALYQcqXo3qvvibRhCLdbWrGitW33OjAG/5UDsIVGV4p+M6Qw0d1Ae6mp0q23Jk8dGINLPAAAwDiMoACwBVe4SVfse1+S9NbAS5jq3hSNjdLrrze3i4piO9V9POrAGHzCAGzB3RjSit88LImp7o0SDErXXNPcjvVU9/GoA2NwiQcAABiHgAIAAIxDQAEAAMbhIh4AxMmAmetOuc3+BcVx6AlgPkZQAACAcQgoAADAOFziAWALIVeKHvruXZE2DOF2Sz//eWvb7nVgDIdlWVaiOxGtQCCgrKws1dTUyOv1Jro7QMKczj0NsBfuQUEyi+bvN5d4AACAcRgnBWALznCThn+2S5K0/awLFGaqezM0NUn/9V/N7W99S3LF6HOJVx0Yg4ACwBY8jSGtfu4BSUx1b5T6euk732luHz0qZWbauw6MwSUeAABgHAIKAAAwDgEFAAAYh4ACAACMQ0ABAADGIaAAAADj8JgxAFtodLn06OjbIm0YIjVVWriwtW33OjAGAQWALYRcqfrXESWJ7gbac7ulGTOSpw6MwSUeAABgHEZQANiCM9ykC6v3SpJ25p7DVPemaGqS3n+/uX3JJbGd6j4edWAMAgoAW/A0hvSfvy6XxFT3Rqmvl4YPb27Heqr7eNSBMbjEAwAAjBNVQJk/f74uu+wyde/eXTk5Obruuuu0Z8+eNtvU19errKxMPXv2VLdu3VRSUqLq6uo221RVVam4uFgZGRnKycnRjBkz1NjY2PGjAQAASSGqgLJlyxaVlZVp27Zt2rhxo0KhkMaOHava2trINtOnT9err76qNWvWaMuWLTp48KAmTJgQeb+pqUnFxcVqaGjQ1q1b9cwzz2jlypWaPXt25x0VAACwtajuQVm/fn2b1ytXrlROTo4qKyt1xRVXqKamRsuXL9eqVas0ZswYSdKKFSs0ePBgbdu2TSNHjtSGDRu0e/duvfHGG8rNzdXFF1+sRx55RPfff7/mzp0rt9vdeUcHAABsqUP3oNTU1EiSevToIUmqrKxUKBRSYWFhZJtBgwapX79+qqiokCRVVFRoyJAhys3NjWxTVFSkQCCgXbt2nbBOMBhUIBBoswAAgOR1xgElHA7rnnvu0eWXX64LL7xQkuT3++V2u5Wdnd1m29zcXPn9/sg2fxlOWt5vee9E5s+fr6ysrMiSn59/pt0GAAA2cMaPGZeVlWnnzp16++23O7M/JzRr1iyVl5dHXgcCAUIK0MU0ulx64vKbIu1kNWDmulNus39BcRx6cppSU6U5c1rbdq8DY5xRQJk6darWrl2rt956S2eddVZkfV5enhoaGnTkyJE2oyjV1dXKy8uLbLN9+/Y2+2t5yqdlm/Y8Ho88Hs+ZdBVAkgi5UvXEqImJ7gbac7uluXOTpw6MEdUlHsuyNHXqVL300kt68803NXDgwDbvDxs2TKmpqdq0aVNk3Z49e1RVVSWfzydJ8vl8+uCDD3To0KHINhs3bpTX61VBQUFHjgUAACSJqEZQysrKtGrVKr3yyivq3r175J6RrKwspaenKysrS5MnT1Z5ebl69Oghr9eradOmyefzaeTIkZKksWPHqqCgQJMmTdLChQvl9/v14IMPqqysjFESACflsML65pcHJEmf9MqX5WCeSSOEw9KHHza3Bw+WnDH6XOJVB8aIKqAsXbpUkjR69Og261esWKFbb71VkrRo0SI5nU6VlJQoGAyqqKhIS5YsiWzrcrm0du1aTZkyRT6fT5mZmSotLdW8efM6diQAklpaqEEbny6T1DLVfVqCewRJUl2d9H8PSsR0Cvp41YExogoolmWdcpu0tDQtXrxYixcvPuk2/fv312uvvRZNaQAA0IUwRgYAAIzDtxkDgM3Y7lFk4AwwggIAAIxDQAEAAMYhoAAAAONwDwoAW2h0ufSL4RMibRgiNVW6997Wtt3rwBgEFAC2EHKlav53bk90N9Ce2y099ljy1IExuMQDAACMwwgKAFtwWGF9I/CFJOnP3t5MdW+KcFiqqmpu9+sX26nu41EHxiCgALCFtFCD3l42WRJT3Rulrk5q+eLYWE91H486MAYRFAAAGIeAAgAAjENAAQAAxuEeFMBQp/N9KwCQrBhBAQAAxiGgAAAA43CJB4AtNDld+vXfFkfaMERKinT33a1tu9eBMfiUAdhCQ0qqZo+dkuhuoD2PR1q8OHnqwBhc4gEAAMZhBAWAPViWetQFJElfpXslhyPBHYIkybKkL79sbvfqFbvPJV51YAwCCgBbSA8F9f7PJkpiqnujHDsm5eQ0t2M5BX286sAYXOIBAADGIaAAAADjEFAAAIBxCCgAAMA4BBQAAGAcAgoAADAOjxkDsIUmp0u/ufDKSBuGSEmRSktb23avA2PwKQOwhYaUVN1bPD3R3UB7Ho+0cmXy1IExuMQDAACMwwgKAHuwLKWHgpKkulQPU52bwrKaZ3mVpIyM2E51H486MAYBBYAtpIeC+nDR9ZKY6v50DJi57pTb7F9Q3PFCx45J3bo1t2M91X086sAYXOIBAADGIaAAAADjEFAAAIBxuAcFSIDTuT8AALoyRlAAAIBxCCgAAMA4XOIBYAthp1Przr880oYhXC7p+utb23avA2MQUADYQjDFrbLrZiW6G2gvLU1asyZ56sAY/G8IAAAwDgEFAAAYh4ACwBbSG+q1/yfXaP9PrlF6Q32iu4MWtbXN34vjcDS37V4HxiCgAAAA4xBQAACAcQgoAADAOAQUAABgHAIKAAAwDgEFAAAYh5lkAdhC2OnUm2dfGmnDEC6XdPXVrW2714ExCCgAbCGY4tbt35+b6G6gvbQ0ad265KkDY/C/IQAAwDgEFAAAYBwCCgBbSG+o1+7HS7T78RKmujdJba2Umdm8xHqq+3jUgTG4BwWAbWSEgonuQlIZMPPU93TsX1B86h0dO9YJvTkN8aoDIzCCAgAAjBP1CMpbb72lxx57TJWVlfr888/10ksv6brrrou8b1mW5syZo1/+8pc6cuSILr/8ci1dulTnnntuZJuvvvpK06ZN06uvviqn06mSkhI9+eST6tatW6ccFJBIp/N/pQCAvy7qEZTa2lpddNFFWrx48QnfX7hwoZ566iktW7ZM7777rjIzM1VUVKT6+tZrxhMnTtSuXbu0ceNGrV27Vm+99ZbuvPPOMz8KAACQVKIeQRk3bpzGjRt3wvcsy9ITTzyhBx98UOPHj5ck/frXv1Zubq5efvll3Xjjjfrwww+1fv167dixQ5de2jzp0s9+9jNdffXV+ulPf6q+fft24HAAAEAy6NR7UPbt2ye/36/CwsLIuqysLI0YMUIVFRWSpIqKCmVnZ0fCiSQVFhbK6XTq3XffPeF+g8GgAoFAmwUAACSvTn2Kx+/3S5Jyc3PbrM/NzY285/f7lZOT07YTKSnq0aNHZJv25s+fr4cffrgzuwrAZsIOh7blXxhpwxBOp/Ttb7e27V4HxrDFY8azZs1SeXl55HUgEFB+fn4CewQg3oKpHt1484JEdwPtpadLmzcnTx0Yo1NjaF5eniSpurq6zfrq6urIe3l5eTp06FCb9xsbG/XVV19FtmnP4/HI6/W2WQAAQPLq1IAycOBA5eXladOmTZF1gUBA7777rnw+nyTJ5/PpyJEjqqysjGzz5ptvKhwOa8SIEZ3ZHQAAYFNRX+I5evSoPvnkk8jrffv26fe//7169Oihfv366Z577tG//Mu/6Nxzz9XAgQP10EMPqW/fvpG5UgYPHqyrrrpKd9xxh5YtW6ZQKKSpU6fqxhtv5AkeACeV3lCvt5fdLkkaddfTqnOnJbhHkNQ87fyAAc3t/fubp6K3cx0YI+qA8t577+k73/lO5HXLvSGlpaVauXKl7rvvPtXW1urOO+/UkSNHNGrUKK1fv15paa3/mDz77LOaOnWqrrzyyshEbU899VQnHA6AZNazjif4jPTll8lVB0aIOqCMHj1almWd9H2Hw6F58+Zp3rx5J92mR48eWrVqVbSlgYRjllgAiA+e1QIAAMYhoAAAAOMQUAAAgHEIKAAAwDi2mEkWAMIOh/6Qd26kDUM4nVLLd6vFeqr7eNSBMQgoAGwhmOrR+NJFie4G2ktPl3bsSJ46MAYxFAAAGIeAAgAAjENAAWALaaF6vb30dr299HalheoT3R20OHaseQr6AQOa23avA2NwDwoAW3BY0lmBQ5E2DGFZ0p/+1Nq2ex0YgxEUAABgHAIKAAAwDgEFAAAYh4ACAACMQ0ABAADG4SkeALZgOaT/6dkv0oYhHA6poKC1bfc6MAYBBYAt1KemaewPlyS6G2gvI0PatSt56sAYXOIBAADGIaAAAADjEFAA2EJaqF4bfnW3Nvzqbqa6N8mxY9IFFzQvsZ7qPh51YAzuQQH+z4CZ6xLdBfwVDks673BVpA1DWJa0e3dr2+51YAxGUAAAgHEYQQEAnNSpRhbTG+r1YZz6gq6FERQAAGAcRlDQJXB/CQDYCyMoAADAOIygALAFyyF95s2JtGEGyyGpf//mF7Ge6j4edWAMAgoAW6hPTdOoKU8nuhtopz41Tdq/P/aFMjLiUwfG4BIPAAAwDgEFAAAYh4ACwBY8oaBeeWa6XnlmujyhYKK7g//jCQWlyy5rXurqYleori4+dWAM7kEBYAtOy9JF/o8jbZjBaVnSe+81vwiHY1coHI5PHRiDERQAAGAcAgoAADAOl3hgtNOZAXb/guI49AQAEE+MoAAAAOMQUAAAgHG4xAPb44sAu47D6d5EdwEn0qtXctWBEQgoAGyhzp2mYf+0KtHdQDt17jTpiy9iXygzMz51YAwu8QAAAOMwggIA6BCetkMsMIICwBY8oaBWr5qp1atmMtW9QeL2udTVSaNHNy9Mdd8lMIICwBaclqWRB3ZG2jBD3D6XcFjasqW1jaTHCAoAADAOAQUAABiHgAIAAIxDQAEAAMbhJlkkDDPAAgBOhoACwDaOpXoS3QWcQNw+l4yM+NSBEQgoAGyhzp2mgvIXEt0NtBO3zyUzU6qtjX0dGIN7UAAAgHEYQUHUmNYaABBrBBTEBDfAorN5Ghu09KVHJUlTvveAginuBPcI0ul/Lh39N6GlzphBOdILL0hpaR3aH8xHQAFgC85wWGM+fS/Shhni9blE6nwqqakpZnVgDgIKAKDL4VK1+QgoaINLMwAAEyT0KZ7FixdrwIABSktL04gRI7R9+/ZEdgcAABgiYSMo//Ef/6Hy8nItW7ZMI0aM0BNPPKGioiLt2bNHOTk5iepWUmN0BAAg2eMSV8ICyuOPP6477rhDt912myRp2bJlWrdunZ5++mnNnDkzUd0ykh1+kQAgXgY/tF517pM/xcO/h8khIQGloaFBlZWVmjVrVmSd0+lUYWGhKioqjts+GAwqGAxGXtfU1EiSAoFA7Dt7EhfOef2U2+x8uKhTaoWDx065zemci9PZD2CqpoZ6tfyWNwWPKWzxJI8J4vW5RFOns/49TOTfmFhL1PG37NOyrFNvbCXAn//8Z0uStXXr1jbrZ8yYYQ0fPvy47efMmWNJYmFhYWFhYUmC5cCBA6fMCrZ4imfWrFkqLy+PvA6Hw/rqq6/Us2dPORyOmNUNBALKz8/XgQMH5PV6Y1YHx+PcJw7nPnE494nDuY8Py7L09ddfq2/fvqfcNiEBpVevXnK5XKqurm6zvrq6Wnl5ecdt7/F45PG0/bbM7OzsWHaxDa/Xyy9sgnDuE4dznzic+8Th3MdeVlbWaW2XkMeM3W63hg0bpk2bNkXWhcNhbdq0ST6fLxFdAgAABknYJZ7y8nKVlpbq0ksv1fDhw/XEE0+otrY28lQPAADouhIWUG644QZ98cUXmj17tvx+vy6++GKtX79eubm5ierScTwej+bMmXPc5SXEHuc+cTj3icO5TxzOvXkclnU6z/oAAADET0KnugcAADgRAgoAADAOAQUAABiHgAIAAIzT5QLK4sWLNWDAAKWlpWnEiBHavn37SbcdPXq0HA7HcUtxcesXUVmWpdmzZ6tPnz5KT09XYWGhPv7443gciu109rm/9dZbj3v/qquuiseh2E40516SnnjiCZ1//vlKT09Xfn6+pk+frvr6+g7ts6vq7HM/d+7c437vBw0aFOvDsKVozn0oFNK8efN0zjnnKC0tTRdddJHWr1/foX2igzrly3VsYvXq1Zbb7baefvppa9euXdYdd9xhZWdnW9XV1Sfc/vDhw9bnn38eWXbu3Gm5XC5rxYoVkW0WLFhgZWVlWS+//LL1hz/8wbr22mutgQMHWnV1dXE6KnuIxbkvLS21rrrqqjbbffXVV3E6IvuI9tw/++yzlsfjsZ599llr37591uuvv2716dPHmj59+hnvs6uKxbmfM2eOdcEFF7T5vf/iiy/idUi2Ee25v++++6y+ffta69ats/bu3WstWbLESktLs95///0z3ic6pksFlOHDh1tlZWWR101NTVbfvn2t+fPnn9bPL1q0yOrevbt19OhRy7IsKxwOW3l5edZjjz0W2ebIkSOWx+Oxnnvuuc7tvM119rm3rOaAMn78+M7uatKJ9tyXlZVZY8aMabOuvLzcuvzyy894n11VLM79nDlzrIsuuigm/U0m0Z77Pn36WD//+c/brJswYYI1ceLEM94nOqbLXOJpaGhQZWWlCgsLI+ucTqcKCwtVUVFxWvtYvny5brzxRmVmZkqS9u3bJ7/f32afWVlZGjFixGnvsyuIxblvsXnzZuXk5Oj888/XlClTdPjw4U7tu92dybn/u7/7O1VWVkaGrj/99FO99tpruvrqq894n11RLM59i48//lh9+/bV2WefrYkTJ6qqqip2B2JDZ3Lug8Gg0tLS2qxLT0/X22+/fcb7RMfY4tuMO8OXX36ppqam42aqzc3N1UcffXTKn9++fbt27typ5cuXR9b5/f7IPtrvs+U9xObcS9JVV12lCRMmaODAgdq7d68eeOABjRs3ThUVFXK5XJ16DHZ1Juf+5ptv1pdffqlRo0bJsiw1Njbqrrvu0gMPPHDG++yKYnHuJWnEiBFauXKlzj//fH3++ed6+OGH9a1vfUs7d+5U9+7dY3pMdnEm576oqEiPP/64rrjiCp1zzjnatGmTXnzxRTU1NZ3xPtExXWYEpaOWL1+uIUOGaPjw4YnuSpdzsnN/44036tprr9WQIUN03XXXae3atdqxY4c2b96cmI4mic2bN+vRRx/VkiVL9P777+vFF1/UunXr9MgjjyS6a0nvdM79uHHj9P3vf19Dhw5VUVGRXnvtNR05ckTPP/98Antuf08++aTOPfdcDRo0SG63W1OnTtVtt90mp5M/k4nSZc58r1695HK5VF1d3WZ9dXW18vLy/urP1tbWavXq1Zo8eXKb9S0/dyb77Epice5P5Oyzz1avXr30ySefdKi/yeRMzv1DDz2kSZMm6Yc//KGGDBmi733ve3r00Uc1f/58hcPhDn2eXUkszv2JZGdn67zzzuP3/i+cybnv3bu3Xn75ZdXW1upPf/qTPvroI3Xr1k1nn332Ge8THdNlAorb7dawYcO0adOmyLpwOKxNmzbJ5/P91Z9ds2aNgsGg/uEf/qHN+oEDByovL6/NPgOBgN59991T7rMricW5P5HPPvtMhw8fVp8+fTrc52RxJuf+2LFjx/1fY8slM8uyOvR5diWxOPcncvToUe3du5ff+7/Qkd/RtLQ0feMb31BjY6NeeOEFjR8/vsP7xBlK7D268bV69WrL4/FYK1eutHbv3m3deeedVnZ2tuX3+y3LsqxJkyZZM2fOPO7nRo0aZd1www0n3OeCBQus7Oxs65VXXrH++Mc/WuPHj+cx4xPo7HP/9ddfW/fee69VUVFh7du3z3rjjTesSy65xDr33HOt+vr6mB+PnUR77ufMmWN1797deu6556xPP/3U2rBhg3XOOedYP/jBD057n2gWi3P/z//8z9bmzZutffv2We+8845VWFho9erVyzp06FDcj89k0Z77bdu2WS+88IK1d+9e66233rLGjBljDRw40Prf//3f094nOleXCiiWZVk/+9nPrH79+llut9saPny4tW3btsh73/72t63S0tI223/00UeWJGvDhg0n3F84HLYeeughKzc31/J4PNaVV15p7dmzJ5aHYFudee6PHTtmjR071urdu7eVmppq9e/f37rjjjv4h+Ikojn3oVDImjt3rnXOOedYaWlpVn5+vnX33Xe3+Yf6VPtEq84+9zfccIPVp08fy+12W9/4xjesG264wfrkk0/ieET2Ec2537x5szV48GDL4/FYPXv2tCZNmmT9+c9/jmqf6FwOyzrJuCEAAECCdJl7UAAAgH0QUAAAgHEIKAAAwDgEFAAAYBwCCgAAMA4BBQAAGIeAAgAAjENAAQAAxiGgAAAA4xBQAACAcQgoAADAOAQUAABgnP8Pdvv5Ih8gfAQAAAAASUVORK5CYII=",
      "text/plain": [
       "<Figure size 640x480 with 1 Axes>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "scores = get_all_scores(rst_all_citation)\n",
    "score_threshold = np.percentile(scores, 99.7)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 25,
   "id": "c3e3b612",
   "metadata": {},
   "outputs": [],
   "source": [
    "rst_abstract = rearrange_query_results(rst_all_abstract, score_threshold=0.8)\n",
    "rst_citation = rearrange_query_results(rst_all_citation, score_threshold=score_threshold)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 26,
   "id": "53ae2b32",
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "1. prochaska & zheng 2019, arXiv:1901.11051\n",
      "\n",
      "Reasons of citation: cited to support the statement that fast radio bursts (frbs) are starting to provide an unprecedented probe of the circumgalactic medium (cgm) with the advent of many radio telescopes geared for this science coming online across the globe.; these papers are referenced to highlight the focus on studying the circumgalactic medium (cgm) with frbs, emphasizing its significance in galaxy formation and feedback processes and the challenges in detecting the cgm using traditional means.\n",
      "\n",
      "2. ravi 2019, arXiv:1804.07291\n",
      "\n",
      "Reasons of citation: these papers are referenced to highlight the focus on studying the circumgalactic medium (cgm) with frbs, emphasizing its significance in galaxy formation and feedback processes and the challenges in detecting the cgm using traditional means.; cited to highlight the use of dispersion measure (dm) measurements from fast radio bursts (frbs) to constrain the circumgalactic medium (cgm) gas distribution.\n",
      "\n",
      "3. petroff et al. 2019; petroff et al. 2019a; petroff et al. 2019b, arXiv:1904.07947\n",
      "\n",
      "Reasons of citation: cited to support the statement that fast radio bursts (frbs) are starting to provide an unprecedented probe of the circumgalactic medium (cgm) with the advent of many radio telescopes geared for this science coming online across the globe.\n",
      "\n",
      "4. petroff et al. 2022; petroff et al. 2021, arXiv:2107.10113\n",
      "\n",
      "Reasons of citation: cited to support the statement that fast radio bursts (frbs) are starting to provide an unprecedented probe of the circumgalactic medium (cgm) with the advent of many radio telescopes geared for this science coming online across the globe.\n",
      "\n",
      "5. keating & pen 2020, arXiv:2001.11105\n",
      "\n",
      "Reasons of citation: cited to support the statement that fast radio bursts (frbs) are starting to provide an unprecedented probe of the circumgalactic medium (cgm) with the advent of many radio telescopes geared for this science coming online across the globe.\n",
      "\n"
     ]
    }
   ],
   "source": [
    "# citation search\n",
    "print_citation_query_results(citation_data, rst_citation)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 27,
   "id": "940a48bc",
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "1. Vikram Ravi, , arXiv:1804.07291\n",
      "\n",
      "Title: Measuring the circum- and inter-galactic baryon contents with fast radio bursts\n",
      "\n",
      "Abstract: Over 80% of the cosmic baryon density is likely to be distributed in the diffuse, $\\gtrsim10^{4}$ K circum- and inter-galactic medium (CGM and IGM respectively). We demonstrate that the dispersion measures (DMs) of samples of localized Fast Radio Bursts (FRBs) can be used to measure the distribution of baryons between the CGM and IGM. We propose to separate the CGM and IGM contributions to FRB DMs by including redshift and mass measurements of intervening galaxies. Using simulated samples of FRB sightlines through intervening galaxy halos and an illustrative model for the CGM, and including realistic observational uncertainties, we show that small samples ($O(10^{1})-O(10^{2})$) of localized FRBs are sensitive to the presence of CGM gas. The fractions of baryons in the CGM and IGM can be accurately estimated with 100 FRBs at $z<1$, and the characteristic radial density profiles of CGM halos may also be possible to constrain. The required samples of localized FRBs are expected to be assembled in the coming few years by instruments such as the Australian Square Kilometre Array Pathfinder, the Very Large Array, and the Deep Synoptic Array.\n",
      "\n",
      "2. Z. Zheng et al., arXiv:1409.3244\n",
      "\n",
      "Title: Probing the Intergalactic Medium with Fast Radio Bursts\n",
      "\n",
      "Abstract: The recently discovered fast radio bursts (FRBs), presumably of extra-galactic origin, have the potential to become a powerful probe of the intergalactic medium (IGM). We point out a few such potential applications. We provide expressions for the dispersion measure and rotation measure as a function of redshift, and we discuss the sensitivity of these measures to the HeII reionization and the IGM magnetic field. Finally we calculate the microlensing effect from an isolate, extragalctic stellar-mass compact object on the FRB spectrum. The time delays between the two lensing images will induce constructive and destructive interference, leaving a specific imprint on the spectra of FRBs. With a high all-sky rate, a large statistical sample of FRBs is expected to make these applications feasible.\n",
      "\n",
      "3. Vikram Ravi et al., arXiv:1903.06535\n",
      "\n",
      "Title: Fast Radio Burst Tomography of the Unseen Universe\n",
      "\n",
      "Abstract: The discovery of Fast Radio Bursts (FRBs) at cosmological distances has opened a powerful window on otherwise unseen matter in the Universe. In the 2020s, observations of $>10^{4}$ FRBs will assess the baryon contents and physical conditions in the hot/diffuse circumgalactic, intracluster, and intergalactic medium, and test extant compact-object dark matter models.\n",
      "\n",
      "4. Xiaohan Wu, Matthew McQuinn, , arXiv:2209.04455\n",
      "\n",
      "Title: A measurement of circumgalactic gas around nearby galaxies using fast radio bursts\n",
      "\n",
      "Abstract: The distribution of gas in the circumgalactic medium (CGM) of galaxies of all types is poorly constrained. Foreground CGMs contribute an extra amount to the dispersion measure (DM) of fast radio bursts (FRB). We measure this DM excess for the CGMs of $10^{11}-10^{13}\\ M_\\odot$ halos using the CHIME/FRB first data release, a halo mass range that is challenging to probe in any other way. Because of the uncertainty in the FRBs' angular coordinates, only for nearby galaxies is the localization sufficient to confidently associate them with intersecting any foreground halo. Thus we stack on galaxies within 80 Mpc, optimizing the stacking scheme to approximately minimize the stack's variance and marginalize over uncertainties in FRB locations. The sample has 20-30 FRBs intersecting halos with masses of $10^{11}-10^{12}\\ M_\\odot$ and also of $10^{12}-10^{13}\\ M_\\odot$, and these intersections allow a marginal $1-2\\,\\sigma$ detection of the DM excess in both mass bins. The $10^{11}-10^{12}\\ M_\\odot$ halos bin also shows a DM excess at 1-2 virial radii. By comparing data with different models for the CGM gas profile, we find that all models are favored by the data up to 2-$\\sigma$ level compared to the null hypothesis of no DM excess. With 2000-3000 more bursts from a future CHIME data release, we project a 4-$\\sigma$ detection of the CGM. Distinguishing between viable CGM models by stacking FRBs with CHIME-like localization would require tens of thousands of bursts.\n",
      "\n",
      "5. Takuya Akahori, Dongsu Ryu, B. M. Gaensler, , arXiv:1602.03235\n",
      "\n",
      "Title: Fast Radio Bursts as Probes of Magnetic Fields in the Intergalactic Medium\n",
      "\n",
      "Abstract: We examine the proposal that the dispersion measures (DMs) and Faraday rotation measures (RMs) of extragalactic linearly-polarized fast radio bursts (FRBs) can be used to probe the intergalactic magnetic field (IGMF) in filaments of galaxies. The DM through the cosmic web is dominated by contributions from the warm-hot intergalactic medium (WHIM) in filaments and from the gas in voids. On the other hand, RM is induced mostly by the hot medium in galaxy clusters, and only a fraction of it is produced in the WHIM. We show that if one excludes FRBs whose sightlines pass through galaxy clusters, the line-of-sight strength of the IGMF in filaments, $B_{||}$, is approximately $C(\\langle 1+z \\rangle/f_{DM})(RM/DM)$, where $C$ is a known constant. Here, the redshift of the FRB is not required to be known; $f_{DM}$ is the fraction of total DM due to the WHIM, while $\\langle 1+z \\rangle$ is the redshift of interevening gas weighted by the WHIM gas density, both of which can be evaluated for a given cosmology model solely from the DM of an FRB. Using data on structure formation simulations and a model IGMF, we show that $C(\\langle 1+z \\rangle/f_{DM})(RM/DM)$ closely reproduces the density-weighted line-of-sight strength of the IGMF in filaments of the large-scale structure.\n",
      "\n"
     ]
    }
   ],
   "source": [
    "# abstract search\n",
    "print_abstract_query_results(abstract_data, rst_abstract)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "id": "5aeac330",
   "metadata": {},
   "outputs": [],
   "source": []
  },
  {
   "cell_type": "markdown",
   "id": "883caee2",
   "metadata": {},
   "source": [
    "### example where combining both results should be better\n",
    "Because the abstract search results can return more recent papers when it comes to gravitational lensing, and also the topic better.  But citation search has Lorimer and Cordes & Chatterjee which either are not relevant at all or come from noise like the CGM example above.\n",
    "\n",
    "It seems like openai embedding matches Lorimer and Petroff etc. to everything... need to find out a way of prohibiting that."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 28,
   "id": "bcdc4096",
   "metadata": {},
   "outputs": [],
   "source": [
    "q = 'fast radio burst gravitational lensing'"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 29,
   "id": "9a0d22dd",
   "metadata": {},
   "outputs": [],
   "source": [
    "rst_all_abstract = retriever_abstract.retrieve(q)\n",
    "rst_all_citation = retriever_citation.retrieve(q)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 30,
   "id": "215f4587",
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAigAAAGdCAYAAAA44ojeAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjcuMywgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/OQEPoAAAACXBIWXMAAA9hAAAPYQGoP6dpAAAqH0lEQVR4nO3de3xU9Z3/8ffMJDO5wCTLJQlUELRWQFFXFMhqW4spEWPFgha3LEZKdcXgr5ItCj4UvLTCog/x0iDdXQS7i7LFKm4Bb8WCqwRB/LUroPzUQqHFBNQlI7lMLnN+f0yZNIImk3yH883J6/l4nMfjm5mTD9/D94S8+Z5zvuNzHMcRAACARfxudwAAAODzCCgAAMA6BBQAAGAdAgoAALAOAQUAAFiHgAIAAKxDQAEAANYhoAAAAOukud2BzojFYjp48KB69+4tn8/ndncAAEAHOI6jzz77TAMHDpTf/+VzJN0yoBw8eFCDBg1yuxsAAKATDhw4oFNOOeVL9+mWAaV3796S4gcYDodd7g0AfE5trTRwYLx98KCUne1uf2COqbHtoedIJBLRoEGDEr/Hv0y3DCjHLuuEw2ECCgD7BAKt7XC4x/zy6RFMjW0PP0c6cnsGN8kCAADrdMsZFACwWlqaVFra2oZ3mBpbzpF2+RzHcdzuRLIikYhycnJUU1PDJR4AALqJZH5/c4kHAABYh3klADDNcaS6ung7K0tivSbvMDW2nCPtYgYFAEyrq5N69Ypvx34JwRtMjS3nSLsIKAAAwDoEFAAAYB0CCgAAsA4BBQAAWIeAAgAArENAAQAA1mEdFAAwLRCQrr66tQ3vMDW2nCPtYql7AABwUiTz+5sZFAA4SYbMXd/uPvsWlZyEngD24x4UAABgHQIKAJhWWxv/bBWfL96Gd5gaW86RdhFQAACAdQgoAADAOgQUAABgHQIKAACwDgEFAABYh4ACAACsw0JtAGBaICBdfnlrG95hamw5R9pFQAEA0zIypPXtrxqLbsjU2HKOtItLPAAAwDoEFAAAYB0CCgCYVlsrZWfHN5Yx9xZTY8s50i7uQQGAVKirc7sHSBVTY8s58qWYQQEAANYhoAAAAOsQUAAAgHUIKAAAwDoEFAAAYB2e4gEA0/x+6ZvfbG3DO0yNLedIuwgoAGBaZqa0aZPbvUAqmBpbzpF2dSm2LVq0SD6fT7feemvitYaGBpWVlalv377q1auXJk+erOrq6jbft3//fpWUlCgrK0t5eXmaM2eOmpubu9IVAADgIZ0OKNu3b9fPf/5znXPOOW1enz17tn79619rzZo12rx5sw4ePKhJkyYl3m9paVFJSYkaGxu1ZcsWPfnkk1q5cqXmz5/f+aMAAACe0qmAcvToUU2dOlX/+q//qr/5m79JvF5TU6Ply5froYce0rhx4zRq1CitWLFCW7Zs0datWyVJL7/8snbv3q3/+I//0HnnnacJEybovvvuU0VFhRobG80cFQC4qbZW6t8/vrGMubeYGlvOkXZ1KqCUlZWppKRERUVFbV7fsWOHmpqa2rw+bNgwDR48WJWVlZKkyspKjRw5Uvn5+Yl9iouLFYlEtGvXrhP+edFoVJFIpM0GAFb7+OP4Bu8xNbacI18q6ZtkV69erbffflvbt28/7r2qqioFg0Hl5ua2eT0/P19VVVWJff46nBx7/9h7J7Jw4ULdc889yXYVAAB0U0nNoBw4cEA/+tGPtGrVKmVkZKSqT8eZN2+eampqEtuBAwdO2p8NAABOvqQCyo4dO3To0CGdf/75SktLU1pamjZv3qxHH31UaWlpys/PV2Njo44cOdLm+6qrq1VQUCBJKigoOO6pnmNfH9vn80KhkMLhcJsNAAB4V1KXeC699FK98847bV6bPn26hg0bpttvv12DBg1Senq6Nm7cqMmTJ0uS9uzZo/3796uwsFCSVFhYqJ/+9Kc6dOiQ8vLyJEmvvPKKwuGwRowYYeKYAMDThsxd3+4++xaVnISeAKmTVEDp3bu3zj777DavZWdnq2/fvonXZ8yYofLycvXp00fhcFi33HKLCgsLNXbsWEnS+PHjNWLECE2bNk2LFy9WVVWV7rzzTpWVlSkUChk6LAAA0J0ZX0l2yZIl8vv9mjx5sqLRqIqLi7V06dLE+4FAQOvWrdPMmTNVWFio7OxslZaW6t577zXdFQBwh98vXXBBaxveYWpsOUfa5XMcx3G7E8mKRCLKyclRTU0N96MA6DZMXZrhEg+6q2R+fxPbAACAdQgoAADAOgQUADCtrk4aMiS+1dW53RuYZGpsOUfaZfwmWQDo8RxH+uMfW9vwDlNjyznSLmZQAACAdQgoAADAOgQUAABgHQIKAACwDgEFAABYh6d4AMA0n0869uGnPp+7fYFZpsaWc6RdBBQAMC0rS9q1y+1eIBVMjS3nSLu4xAMAAKzDDAoAtIMP5wNOPmZQAMC0ujrprLPiG8uYe4upseUcaRczKABgmuNIu3e3tuEdpsaWc6RdzKAAAADrEFAAAIB1CCgAAMA6BBQAAGAdAgoAALAOT/EAgGk+n3Tqqa1teIepseUcaRcBBYBnubbAWlaWtG+f+bpwn6mx5RxpF5d4AACAdQgoAADAOgQUADCtvl668ML4Vl/vdm9gkqmx5RxpF/egAIBpsZj01lutbXiHqbHlHGkXMygAAMA6BBQAAGAdLvEAgEU68mg00BMwgwIAAKxDQAEAANbhEg8ApEK/fm73AKliamw5R74UAQUATMvOlg4fdrsXSAVTY8s50i4u8QAAAOsQUAAAgHUIKABgWn29dMkl8Y1lzL3F1NhyjrSLe1AAwLRYTNq8ubUN7zA1tpwj7WIGBQAAWIeAAgAArMMlHgA9GkvLA3ZiBgUAAFiHgAIAAKzDJR4ASIWsLFf/+I5cutq3qOQk9MSDTI2ty+eI7QgoAGBadrZUW+t2L5AKpsaWc6RdXOIBAADWIaAAAADrEFAAwLSGBqmkJL41NLjdG5hkamw5R9rFPSgAYMBf35Sa2digdzdskCQNv2O96oMZbnULprW0SH8ZW7W0uF/Hw5hBAQAA1iGgAAAA6xBQAACAdQgoAADAOgQUAABgHQIKAACwDo8ZA4Bh9cEMDbl9ndvdQCpkZ0uOY08dD2MGBQAAWIeAAgAArENAAQDDQs2Nqli7UBVrFyrU3Oh2d2BSQ4N0zTXxratL3Zuo42EEFAAwzB+LqWTPGyrZ84b8sZjb3YFJLS3SM8/Et64udW+ijocRUAAAgHUIKAAAwDoEFAAAYB0CCgAAsA4BBQAAWIeAAgAArMNS9wBgWH16SMNnP5Now0OysqSjR1vbbtfxMAIKAJjm86k+mOF2L5AKPl/8c3RsqeNhXOIBAADWIaAAgGHB5iY9uH6JHly/RMHmJre7A5OiUen66+NbNOp+HQ8joACAYYFYi67euVFX79yoQIxlzD2luVl68sn41tzsfh0PI6AAAADrJBVQHn/8cZ1zzjkKh8MKh8MqLCzUCy+8kHi/oaFBZWVl6tu3r3r16qXJkyerurq6TY39+/erpKREWVlZysvL05w5c9RMegQAAH8lqYByyimnaNGiRdqxY4feeustjRs3ThMnTtSuXbskSbNnz9avf/1rrVmzRps3b9bBgwc1adKkxPe3tLSopKREjY2N2rJli5588kmtXLlS8+fPN3tUAACgW0vqMePvfOc7bb7+6U9/qscff1xbt27VKaecouXLl+upp57SuHHjJEkrVqzQ8OHDtXXrVo0dO1Yvv/yydu/erd/85jfKz8/Xeeedp/vuu0+333677r77bgWDQXNHBgAAuq1O34PS0tKi1atXq7a2VoWFhdqxY4eamppUVFSU2GfYsGEaPHiwKisrJUmVlZUaOXKk8vPzE/sUFxcrEokkZmFOJBqNKhKJtNkAAIB3JR1Q3nnnHfXq1UuhUEg33XSTnnvuOY0YMUJVVVUKBoPKzc1ts39+fr6qqqokSVVVVW3CybH3j733RRYuXKicnJzENmjQoGS7DQAAupGkV5I988wz9bvf/U41NTV65plnVFpaqs2bN6eibwnz5s1TeXl54utIJEJIAWCt+vSQzr9lVaIND8nKkg4dam27XcfDkg4owWBQX/3qVyVJo0aN0vbt2/XII49oypQpamxs1JEjR9rMolRXV6ugoECSVFBQoG3btrWpd+wpn2P7nEgoFFIoxA85gG7C59OnWTlu9wKp4PNJ/fvbU8fDurwOSiwWUzQa1ahRo5Senq6NGzcm3tuzZ4/279+vwsJCSVJhYaHeeecdHTqWGiW98sorCofDGjFiRFe7AgAAPCKpGZR58+ZpwoQJGjx4sD777DM99dRT2rRpk1566SXl5ORoxowZKi8vV58+fRQOh3XLLbeosLBQY8eOlSSNHz9eI0aM0LRp07R48WJVVVXpzjvvVFlZGTMkAJIyZO56t7vwhYLNTbrz1X+TJP1k3A/VmJbuco9gTDQqHbvl4KGHpM7+7jJVx8OSCiiHDh3Sddddp48++kg5OTk655xz9NJLL+nb3/62JGnJkiXy+/2aPHmyotGoiouLtXTp0sT3BwIBrVu3TjNnzlRhYaGys7NVWlqqe++91+xRAYCLArEWXfd/4wFq4SXTJRFQPKO5WTr2e23x4s4HC1N1PCypgLJ8+fIvfT8jI0MVFRWqqKj4wn1OPfVUbdiwIZk/FgAA9DB8Fg8AALAOAQUAAFiHgAIAAKxDQAEAANYhoAAAAOskvZIsAODLNaQHdfFNyxNteEhmprR3b2vb7ToeRkABAMMcn19/yslvf0d0P36/NGSIPXU8jEs8AADAOgQUADAsvaVJ8377hOb99gmltzS53R2Y1NgozZkT3xob3a/jYQQUADAsraVF/7jtWf3jtmeV1tLidndgUlOT9OCD8a2pC+HTVB0PI6AAAADrEFAAAIB1CCgAAMA6BBQAAGAdAgoAALAOAQUAAFiHlWQBwLCG9KC+/YOKRBsekpkp7dzZ2na7jocRUADAMMfn1/v9T3W7G0gFv1866yx76ngYl3gAAIB1mEEBAMPSW5pUVvlLSVJF4ffUFEh3uUcwprFRuv/+ePuOO6RgJy/hmarjYT7HcRy3O5GsSCSinJwc1dTUKBwOu90dAC4YMne92134QpmNDXp3ydWSpOGzn1F9MMPlHp3YvkUlbneh+6mtlXr1irePHpWys92t080k8/ubSzwAAMA6BBQAAGAdAgoAALAOAQUAAFiHgAIAAKxDQAEAANZhHRQA1rH5EeKOiKal68rrHkq04SEZGdK2ba1tt+t4GAEFAAyL+QP6nwFfc7sbSIVAQLrwQnvqeBiXeAAAgHWYQQEAw9JbmjT9rf+SJK244EqWuveSxkbpkUfi7R/9qGtL3Zuo42EEFAAwLK2lRXdsWiFJ+ve/LSGgeElTk3TbbfH2zTd3PliYquNhXOIBAADWYQYFAHqojjwtxQcKwi3MoAAAAOsQUAAAgHUIKAAAwDoEFAAAYB1ukgUAw6Jp6br27+9PtOEhGRnSb3/b2na7jocRUADAsJg/oK2Dz3G7G0iFQEC65BJ76ngYl3gAAIB1mEEBAMPSWpr1979/UZL09LmXqTnAP7We0dQk/cu/xNs33iild/ISnqk6HsZPDQAYlt7SrPteWSZJeubsIgKKlzQ2SrNmxdvXX9/5YGGqjodxiQcAAFiHWA8A+EIshw+3MIMCAACsQ0ABAADWIaAAAADrEFAAAIB1uEkWAAxrTEvX9KsXJNrwkFBIWreute12HQ8joACAYS3+gH57+oVudwOpkJYmlRh4aslUHQ/jEg8AALAOMygAYFhaS7Ou2r1JkrR2xCWsJOslTU3SqlXx9tSpXVvq3kQdD+OnBgAMS29p1oMbHpYkrT/zYgKKlzQ2StOnx9vXXNO1pe5N1PEwLvEAAADrEFAAAIB1CCgAAMA6BBQAAGAdAgoAALAOAQUAAFiHZ98AwLDGtHTdPHFuog0PCYWkX/6yte12HQ8joACAYS3+gDYMu9jtbiAV0tLi65bYUsfDuMQDAACswwwKABgWiLWo+P9VSpJe+lqhWvwBl3sEY5qbpeeei7e/+934TIibdTyMvxEAMCzY3KSlzy+SJA2f/YzqgwQUz4hGpe99L94+erTzwcJUHQ/jEg8AALAOAQUAAFiHOSUAJ9WQuevd7gKAboAZFAAAYB1mUAAYw+wIAFOSmkFZuHChLrzwQvXu3Vt5eXm66qqrtGfPnjb7NDQ0qKysTH379lWvXr00efJkVVdXt9ln//79KikpUVZWlvLy8jRnzhw1Nzd3/WgAAIAnJDWDsnnzZpWVlenCCy9Uc3Oz7rjjDo0fP167d+9Wdna2JGn27Nlav3691qxZo5ycHM2aNUuTJk3SG2+8IUlqaWlRSUmJCgoKtGXLFn300Ue67rrrlJ6ervvvv9/8EQLASdYUSNOPL7810YaHBIPSihWtbbfreJjPcRyns998+PBh5eXlafPmzfrGN76hmpoa9e/fX0899ZSuvvpqSdJ7772n4cOHq7KyUmPHjtULL7ygK664QgcPHlR+fr4kadmyZbr99tt1+PBhBTswUJFIRDk5OaqpqVE4HO5s9wEYxiWenmnfohK3u4BuIpnf3126SbampkaS1KdPH0nSjh071NTUpKKiosQ+w4YN0+DBg1VZGV9VsbKyUiNHjkyEE0kqLi5WJBLRrl27TvjnRKNRRSKRNhsAAPCuTs89xmIx3Xrrrbrooot09tlnS5KqqqoUDAaVm5vbZt/8/HxVVVUl9vnrcHLs/WPvncjChQt1zz33dLarAHBSBWIt+sbetyVJrw093/NL3Xdk5swzsyzNzdJLL8XbxcVdW+reRB0P6/TfSFlZmXbu3KnXX3/dZH9OaN68eSovL098HYlENGjQoJT/uQDQGcHmJq14Jv6fKpa695hoVLriini7q0vdm6jjYZ36G5k1a5bWrVun1157Taecckri9YKCAjU2NurIkSNtZlGqq6tVUFCQ2Gfbtm1t6h17yufYPp8XCoUUCoU601UAANANJXUPiuM4mjVrlp577jm9+uqrGjp0aJv3R40apfT0dG3cuDHx2p49e7R//34VFhZKkgoLC/XOO+/o0KFDiX1eeeUVhcNhjRgxoivHAgAAPCKpGZSysjI99dRTev7559W7d+/EPSM5OTnKzMxUTk6OZsyYofLycvXp00fhcFi33HKLCgsLNXbsWEnS+PHjNWLECE2bNk2LFy9WVVWV7rzzTpWVlTFLAgAAJCUZUB5//HFJ0iWXXNLm9RUrVuj666+XJC1ZskR+v1+TJ09WNBpVcXGxli5dmtg3EAho3bp1mjlzpgoLC5Wdna3S0lLde++9XTsSAADgGUkFlI4smZKRkaGKigpVVFR84T6nnnqqNmzYkMwfDQAAehA+LBAAAFiH55oAwLCmQJru+vZNiTY8JBiUfvaz1rbbdTyMnxwAMKw5kKZ/P/8Kt7uBVEhPl8rK7KnjYVziAQAA1mEGBQAM88daNPpP8c8W23bKWYp5fKn7HqWlRfrv/463v/51KdDJsTVVx8MIKABgWKi5SaufvkMSS917TkOD9K1vxdtHj0rZ2e7W8TAu8QAAAOsQUAAAgHUIKAAAwDoEFAAAYB0CCgAAsA4BBQAAWIfHjAHAsOZAQPdfMj3Rhoekp0uLF7e23a7jYQQUADCsKZCufxkz2e1uIBWCQWnOHHvqeBgBBQCQckPmrm93n32LSk5CT9BdEFAAdEhHfsEgzh9r0dnVH0qSduafzlL3XtLSIr39drx9/vldW+reRB0PI6AAgGGh5ib91y/KJbHUvec0NEijR8fbXV3q3kQdD+MpHgAAYB0CCgAAsA4BBQAAWIeAAgAArENAAQAA1iGgAAAA6/CYMQAY1hwI6OGL/j7RRseYWmsnpQu+padLCxa0tt2u42E+x3EctzuRrEgkopycHNXU1CgcDrvdHaBHYKE2dBesSGuvZH5/c4kHAABYh0s8AGCYz4npqx8fkCR90G+QHB//F/SMWEx69914e/hwyd/JsTVVx8MIKABgWEZTo155okzSsaXuM1zuEYypr5fOPjve7soS9abqeBiRDQAAWIeAAgAArENAAQAA1iGgAAAA6xBQAACAdQgoAADAOjxmDACGNQcC+vnoSYk2PCQ9Xfrxj1vbbtfxMAIKABjWFEjXwm/9wO1uIBWCQemBB+yp42Fc4gEAANZhBgUADPM5MX0lcliS9Odwf5a695JYTNq/P94ePLhrS92bqONhBBQAMCyjqVGvL5shiaXuPae+Xho6NN7u6lL3Jup4GJENAABYh4ACAACswyUeABoyd73bXQCANphBAQAA1iGgAAAA6xBQAACAdbgHBQAMa/EH9Iu/LUm04SFpadLNN7e23a7jYfytAIBhjWnpmj9+ptvdQCqEQlJFhT11PIxLPAAAwDrMoACAaY6jPvURSdKnmWHJ53O5QzDGcaSPP463+/Xr/NiaquNhBBQAMCyzKaq3H5sqiaXuPaeuTsrLi7e7skS9qToexiUeAABgHQIKAACwDgEFAABYh4ACAACsQ0ABAADW4SkeAICndOTTufctKjkJPUFXEFAAwLAWf0DPnH1pog0PSUuTSktb227X8TD+VgDAsMa0dP24ZLbb3UAqhELSypX21PEw7kEBAADWYQYFAExzHGU2RSVJ9ekhljH3EseJrwIrSVlZXVvq3kQdD2MGBQAMy2yK6t0lV+vdJVcnggo8oq5O6tUrvh0LGG7W8TACCgAAsA4BBQAAWIeAAgAArENAAQAA1iGgAAAA6xBQAACAdVgHBQAMi/n9Wn/mRYk2PCQQkK6+urXtdh0PI6AAgGHRtKDKrprndjeQChkZ0po19tTxMKI9AACwDgEFAABYJ+mA8tprr+k73/mOBg4cKJ/Pp7Vr17Z533EczZ8/XwMGDFBmZqaKior0/vvvt9nn008/1dSpUxUOh5Wbm6sZM2bo6NGjXToQALBFZmOD9v3zFdr3z1cos7HB7e7ApNra+Ofm+Hzxttt1PCzpgFJbW6tzzz1XFRUVJ3x/8eLFevTRR7Vs2TK9+eabys7OVnFxsRoaWn9Ip06dql27dumVV17RunXr9Nprr+nGG2/s/FEAAABPSfom2QkTJmjChAknfM9xHD388MO68847NXHiREnSL37xC+Xn52vt2rW69tpr9e677+rFF1/U9u3bdcEFF0iSHnvsMV1++eV68MEHNXDgwC4cDgAA7Rsyd327++xbVHISeoIvYvQelL1796qqqkpFRUWJ13JycjRmzBhVVlZKkiorK5Wbm5sIJ5JUVFQkv9+vN99884R1o9GoIpFImw0AAHiX0YBSVVUlScrPz2/zen5+fuK9qqoq5eXltXk/LS1Nffr0SezzeQsXLlROTk5iGzRokMluAwAAy3SLp3jmzZunmpqaxHbgwAG3uwQAAFLIaEApKCiQJFVXV7d5vbq6OvFeQUGBDh061Ob95uZmffrpp4l9Pi8UCikcDrfZAACAdxldSXbo0KEqKCjQxo0bdd5550mSIpGI3nzzTc2cOVOSVFhYqCNHjmjHjh0aNWqUJOnVV19VLBbTmDFjTHYHAFwR8/v16mkXJNrwkEBAuvzy1rbbdTws6YBy9OhRffDBB4mv9+7dq9/97nfq06ePBg8erFtvvVU/+clPdMYZZ2jo0KG66667NHDgQF111VWSpOHDh+uyyy7TDTfcoGXLlqmpqUmzZs3StddeyxM8QAp05GkFmBVNC+oH19ztdjeQChkZ0noDP1Om6nhY0gHlrbfe0re+9a3E1+Xl5ZKk0tJSrVy5Urfddptqa2t144036siRI7r44ov14osvKiMjI/E9q1at0qxZs3TppZfK7/dr8uTJevTRRw0cDgAA8AKf4ziO251IViQSUU5OjmpqargfBT0asyNA6rAOinnJ/P7m4igAGJbZ2KDdD03W7ocms9S919TWStnZ8a2rS92bqONhRm+SBQDEZTVF3e4CUqWuzq46HsUMCgAAsA4BBQAAWIeAAgAArENAAQAA1iGgAAAA6/AUDwAYFvP5tHXQ2Yk2PMTvl775zda223U8jIACAIZF00O69vuL3O4GUiEzU9q0yZ46HkZsAwAA1iGgAAAA6xBQAMCwzMYG7Xj0+9rx6PdZ6t5ramul/v3jW1eXujdRx8O4BwUAUqBvfcTtLiBVPv7YrjoexQwKAACwDgEFAABYh0s8gKWGzF3vdhcAwDXMoAAAAOsQUAAAgHW4xAMAhsV8Pv2+4IxEGx7i90sXXNDadruOhxFQAMCwaHpIE0uXuN0NpEJmprR9uz11PIyAAriAG2AB4MsxrwQAAKxDQAEAwzKaGvT64z/Q64//QBlNLHXvKXV10pAh8a2uzv06HsYlHgAwzOdIp0QOJdrwEMeR/vjH1rbbdTyMGRQAAGAdAgoAALAOAQUAAFiHgAIAAKxDQAEAANbhKR7AMBZhg+OT/l/fwYk2PMTnk0aMaG27XcfDCCgAYFhDeobG/3Cp291AKmRlSbt22VPHwwgoAACcQEdmQ/ctKjkJPemZuAcFAABYh4ACAIZlNDXo5X+7WS//280sde8xxsa2rk4666z4xlL3J8QlHgAwzOdIX/tkf6IN7zA2to4j7d7d2sZxmEEBAADWIaAAAADrEFAAAIB1CCgAAMA6BBQAAGAdnuIBAMMcn/SncF6iDe8wNrY+n3Tqqa1tHIeAAvwFq0bClIb0DF088wm3u4EUMDa2WVnSvn1dr+NhXOIBAADWIaAAAADrEFAAwLBQU1TPPzlbzz85W6GmqNvdgUHGxra+XrrwwvhWX2+ugx7CPSgAYJjfcXRu1fuJNrzD2NjGYtJbb7W2cRxmUAAAgHWYQUGP0JEndAAA9mAGBQAAWIeAAgAArENAAQAA1uEeFABIgU8yw253ASlibGz79TNTx6MIKABgWH0wQ6P+z1NudwMpYGxss7Olw4e7XsfDuMQDAACsQ0ABAADW4RIPrGbbJwyzngo6ItQU1ZNrFkiSSq+5R9H0kMs9gimfH9tOq6+XJkyIt194QcrMNNA7byGgoNsjNMA2fsfR2AM7E214h7GxjcWkzZtb2zgOl3gAAIB1CCgAAMA6BBQAAGAd7kEBAKCTbLuR30uYQQEAANZhBgUAUqCOR4s9y9jYZmWZqeNRBBQAMKw+mKER5b9yuxtIAWNjm50t1dZ2vY6HcYkHAABYhxkUpASLpwFAHDfSdg4BBQAMCzU36vHn7pckzfzuHYqmBV3uEUwxNbafr4PjEVCQNGZHgC/nj8U07g9vJdrwDlNjyznSPgIKAADdQEf/c+iVy0XcJAsAAKzj6gxKRUWFHnjgAVVVVencc8/VY489ptGjR7vZJU/j0gwA2Il/n4/nWkD5z//8T5WXl2vZsmUaM2aMHn74YRUXF2vPnj3Ky8tzq1sAAHhed3iyyLWA8tBDD+mGG27Q9OnTJUnLli3T+vXr9cQTT2ju3LludavDusPgAgDQXbkSUBobG7Vjxw7Nmzcv8Zrf71dRUZEqKyuP2z8ajSoajSa+rqmpkSRFIpGU9O/sBS8ZqTN49hojdQB0Ly2NDTr2r1NLtE4xh6c0vMLU2KbyHOnI78ZYtM5InWQdq+k4Tvs7Oy7485//7EhytmzZ0ub1OXPmOKNHjz5u/wULFjiS2NjY2NjY2DywHThwoN2s0C0eM543b57Ky8sTX8diMX366afq27evfD7fSe9PJBLRoEGDdODAAYXD4ZP+5yN5jFn3xLh1P4xZ93Myx8xxHH322WcaOHBgu/u6ElD69eunQCCg6urqNq9XV1eroKDguP1DoZBCobafHpmbm5vKLnZIOBzmB7CbYcy6J8at+2HMup+TNWY5OTkd2s+VdVCCwaBGjRqljRs3Jl6LxWLauHGjCgsL3egSAACwiGuXeMrLy1VaWqoLLrhAo0eP1sMPP6za2trEUz0AAKDnci2gTJkyRYcPH9b8+fNVVVWl8847Ty+++KLy8/Pd6lKHhUIhLViw4LjLTrAXY9Y9MW7dD2PW/dg6Zj7H6cizPgAAACcPn8UDAACsQ0ABAADWIaAAAADrEFAAAIB1CCh/UVFRoSFDhigjI0NjxozRtm3bvnDfSy65RD6f77itpKT1wwEdx9H8+fM1YMAAZWZmqqioSO+///7JOJQew/SYXX/99ce9f9lll52MQ+kxkhkzSXr44Yd15plnKjMzU4MGDdLs2bPV0NDQpZpIjukxu/vuu4/7ORs2bFiqD6NHSWbMmpqadO+99+r0009XRkaGzj33XL344otdqmmMkQ/X6eZWr17tBINB54knnnB27drl3HDDDU5ubq5TXV19wv0/+eQT56OPPkpsO3fudAKBgLNixYrEPosWLXJycnKctWvXOr///e+dK6+80hk6dKhTX19/ko7K21IxZqWlpc5ll13WZr9PP/30JB2R9yU7ZqtWrXJCoZCzatUqZ+/evc5LL73kDBgwwJk9e3anayI5qRizBQsWOGeddVabn7PDhw+frEPyvGTH7LbbbnMGDhzorF+/3vnwww+dpUuXOhkZGc7bb7/d6ZqmEFAcxxk9erRTVlaW+LqlpcUZOHCgs3Dhwg59/5IlS5zevXs7R48edRzHcWKxmFNQUOA88MADiX2OHDnihEIh5+mnnzbb+R7K9Jg5TjygTJw40XRX8RfJjllZWZkzbty4Nq+Vl5c7F110UadrIjmpGLMFCxY45557bkr6i+THbMCAAc7PfvazNq9NmjTJmTp1aqdrmtLjL/E0NjZqx44dKioqSrzm9/tVVFSkysrKDtVYvny5rr32WmVnZ0uS9u7dq6qqqjY1c3JyNGbMmA7XxBdLxZgds2nTJuXl5enMM8/UzJkz9cknnxjte0/VmTH7u7/7O+3YsSMxlfyHP/xBGzZs0OWXX97pmui4VIzZMe+//74GDhyo0047TVOnTtX+/ftTdyA9SGfGLBqNKiMjo81rmZmZev311ztd05Ru8WnGqfTxxx+rpaXluBVs8/Pz9d5777X7/du2bdPOnTu1fPnyxGtVVVWJGp+veew9dF4qxkySLrvsMk2aNElDhw7Vhx9+qDvuuEMTJkxQZWWlAoGA0WPoaTozZt///vf18ccf6+KLL5bjOGpubtZNN92kO+64o9M10XGpGDNJGjNmjFauXKkzzzxTH330ke655x59/etf186dO9W7d++UHpPXdWbMiouL9dBDD+kb3/iGTj/9dG3cuFHPPvusWlpaOl3TlB4/g9JVy5cv18iRIzV69Gi3u4IO+qIxu/baa3XllVdq5MiRuuqqq7Ru3Tpt375dmzZtcqejPdymTZt0//33a+nSpXr77bf17LPPav369brvvvvc7hq+QEfGbMKECbrmmmt0zjnnqLi4WBs2bNCRI0f0y1/+0sWe91yPPPKIzjjjDA0bNkzBYFCzZs3S9OnT5fe7Hw/c74HL+vXrp0AgoOrq6javV1dXq6Cg4Eu/t7a2VqtXr9aMGTPavH7s+zpTE+1LxZidyGmnnaZ+/frpgw8+6FJ/0bkxu+uuuzRt2jT98Ic/1MiRI/Xd735X999/vxYuXKhYLNal8wDtS8WYnUhubq6+9rWv8XNmQGfGrH///lq7dq1qa2v1xz/+Ue+995569eql0047rdM1TenxASUYDGrUqFHauHFj4rVYLKaNGzeqsLDwS793zZo1ikaj+od/+Ic2rw8dOlQFBQVtakYiEb355pvt1kT7UjFmJ/KnP/1Jn3zyiQYMGNDlPvd0nRmzurq64/4Xd+xSm+M4XToP0L5UjNmJHD16VB9++CE/ZwZ05WciIyNDX/nKV9Tc3Kxf/epXmjhxYpdrdllKb8HtJlavXu2EQiFn5cqVzu7du50bb7zRyc3NdaqqqhzHcZxp06Y5c+fOPe77Lr74YmfKlCknrLlo0SInNzfXef75553/+Z//cSZOnMhjxgaZHrPPPvvM+fGPf+xUVlY6e/fudX7zm984559/vnPGGWc4DQ0NKT+eniDZMVuwYIHTu3dv5+mnn3b+8Ic/OC+//LJz+umnO9/73vc6XBNdk4ox+6d/+idn06ZNzt69e5033njDKSoqcvr16+ccOnTopB+fFyU7Zlu3bnV+9atfOR9++KHz2muvOePGjXOGDh3q/O///m+Ha6YKAeUvHnvsMWfw4MFOMBh0Ro8e7WzdujXx3je/+U2ntLS0zf7vvfeeI8l5+eWXT1gvFos5d911l5Ofn++EQiHn0ksvdfbs2ZPKQ+hxTI5ZXV2dM378eKd///5Oenq6c+qppzo33HADv+gMS2bMmpqanLvvvts5/fTTnYyMDGfQoEHOzTff3OYfzvZqoutMj9mUKVOcAQMGOMFg0PnKV77iTJkyxfnggw9O4hF5XzJjtmnTJmf48OFOKBRy+vbt60ybNs3585//nFTNVPE5zhfMuwEAALikx9+DAgAA7ENAAQAA1iGgAAAA6xBQAACAdQgoAADAOgQUAABgHQIKAACwDgEFAABYh4ACAACsQ0ABAADWIaAAAADrEFAAAIB1/j9YVZeqVUUWyQAAAABJRU5ErkJggg==",
      "text/plain": [
       "<Figure size 640x480 with 1 Axes>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "scores = get_all_scores(rst_all_citation)\n",
    "score_threshold = np.percentile(scores, 99.7)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 31,
   "id": "c5aa8552",
   "metadata": {},
   "outputs": [],
   "source": [
    "rst_abstract = rearrange_query_results(rst_all_abstract, score_threshold=0.8)\n",
    "rst_citation = rearrange_query_results(rst_all_citation, score_threshold=score_threshold)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 32,
   "id": "fa61f2a3",
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "1. munoz et al. 2016, arXiv:1605.00008\n",
      "\n",
      "Reasons of citation: these references are cited to support the statement that one method for detecting dark compact objects is through gravitational lensing of radio transients, specifically fast radio bursts (frbs), and to provide a basis for discussing the cosmological distance and abundance of frbs as a candidate for time-domain gravitational lensing searches.; these papers are cited to support the statement that if a pair of gravitational waves (gws) and fast radio bursts (frbs) were lensed by an isolated black hole (bh), there would be multiple temporally resolvable frbs due to strong lensing.\n",
      "\n",
      "2. wucknitz et al. 2021, arXiv:2004.11643\n",
      "\n",
      "Reasons of citation: these references are cited to support the statement that one method for detecting dark compact objects is through gravitational lensing of radio transients, specifically fast radio bursts (frbs), and to provide a basis for discussing the cosmological distance and abundance of frbs as a candidate for time-domain gravitational lensing searches.; referenced to highlight the proposal of using strongly gravitationally lensed frbs for conducting cosmological and astrophysical studies.\n",
      "\n",
      "3. cordes & chatterjee 2019, arXiv:1906.05878\n",
      "\n",
      "Reasons of citation: these papers are cited to support the statement that if a pair of gravitational waves (gws) and fast radio bursts (frbs) were lensed by an isolated black hole (bh), there would be multiple temporally resolvable frbs due to strong lensing.; referenced as an example of a study on extragalactic sources of fast radio bursts (frbs).\n",
      "\n",
      "4. lorimer et al. 2007, arXiv:0709.4301\n",
      "\n",
      "Reasons of citation: these references are cited to highlight the discovery of fast radio bursts (frbs) and their millisecond duration radio signals arising from cosmological distances.; referenced to highlight the discovery of fast radio bursts (frbs) through radio astronomy.\n",
      "\n",
      "5. eichler 2017, arXiv:1711.04764\n",
      "\n",
      "Reasons of citation: these references are cited to support the statement that one method for detecting dark compact objects is through gravitational lensing of radio transients, specifically fast radio bursts (frbs), and to provide a basis for discussing the cosmological distance and abundance of frbs as a candidate for time-domain gravitational lensing searches.\n",
      "\n"
     ]
    }
   ],
   "source": [
    "# citation\n",
    "print_citation_query_results(citation_data, rst_citation)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 33,
   "id": "779a4530",
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "1. Mawson W. Sammons et al., arXiv:2210.09487\n",
      "\n",
      "Title: The Effect of Gravitational Lensing on Fast Transient Event Rates\n",
      "\n",
      "Abstract: Fast cosmological transients such as fast radio bursts (FRBs) and gamma-ray bursts (GRBs) represent a class of sources more compact than any other cosmological object. As such they are sensitive to significant magnification via gravitational lensing from a class of lenses which are not well-constrained by observations today. Low-mass primordial black holes are one such candidate which may constitute a significant fraction of the Universe's dark matter. Current observations only constrain their density in the nearby Universe, giving fast transients from cosmological distances the potential to form complementary constraints. Motivated by this, we calculate the effect that gravitational lensing from a cosmological distribution of compact objects would have on the observed rates of FRBs and GRBs. For static lensing geometries, we rule out the prospect that all FRBs are gravitationally lensed for a range of lens masses and show that lens masses greater than $10^{-5}M_\\odot$ can be constrained with 8000 un-localised high fluence FRBs at 1.4GHz, as might be detected by the next generation of FRB-finding telescopes.\n",
      "\n",
      "2. Noah Pearson, Cynthia Trendafilova, Joel Meyers, , arXiv:2009.11252\n",
      "\n",
      "Title: Searching for Gravitational Waves with Strongly Lensed Repeating Fast Radio Bursts\n",
      "\n",
      "Abstract: Since their serendipitous discovery, Fast Radio Bursts (FRBs) have garnered a great deal of attention from both observers and theorists. A new class of radio telescopes with wide fields of view have enabled a rapid accumulation of FRB observations, confirming that FRBs originate from cosmological distances. The high occurrence rate of FRBs and the development of new instruments to observe them create opportunities for FRBs to be utilized for a host of astrophysical and cosmological studies. We focus on the rare, and as yet undetected, subset of FRBs that undergo repeated bursts and are strongly gravitationally lensed by intervening structure. An extremely precise timing of burst arrival times is possible for strongly lensed repeating FRBs, and we show how this timing precision enables the search for long wavelength gravitational waves, including those sourced by supermassive black hole binary systems. The timing of burst arrival for strongly lensed repeating FRBs is sensitive to gravitational wave sources near the FRB host galaxy, which may lie at cosmological distances and would therefore be extremely challenging to detect by other means. Timing of strongly lensed FRBs can also be combined with pulsar timing array data to search for correlated time delays characteristic of gravitational waves passing through the Earth.\n",
      "\n",
      "3. Tal Adi, Ely D. Kovetz, , arXiv:2109.00403\n",
      "\n",
      "Title: Probing Gravitational Slip with Strongly Lensed Fast Radio Bursts\n",
      "\n",
      "Abstract: The rapid accumulation of observed Fast Radio Bursts (FRBs) originating from cosmological distances makes it likely that some will be strongly lensed by intervening matter along the line of sight. Detection of lensed FRB repeaters, which account for a noteworthy fraction of the total population, will allow not only an accurate measurement of the lensing time delay, but also follow-up high-resolution observations to pinpoint the location of the lensed images. Recent works proposed to use such strongly-lensed FRBs to derive constraints on the current expansion rate $ H_{0} $ as well as on cosmic curvature. Here we study the prospects for placing constraints on departures from general relativity via such systems. Using an ensemble of simulated events, we focus on the gravitational slip parameter $\\gamma_{\\rm PN}$ in screened modified gravity models and show that FRB time-delay measurements can yield constraints as tight as $ \\left| \\gamma_{\\rm PN}-1\\right| \\lesssim 0.04\\times(\\Lambda/100\\rm kpc)\\times[N/10]^{-1/2} $ at $1\\sigma$ with $10$ detections.\n",
      "\n",
      "4. Calvin Leung et al., arXiv:2204.06001\n",
      "\n",
      "Title: Constraining Primordial Black Holes using Fast Radio Burst Gravitational-Lens Interferometry with CHIME/FRB\n",
      "\n",
      "Abstract: Fast radio bursts (FRBs) represent an exciting frontier in the study of gravitational lensing, due to their brightness, extragalactic nature, and the compact, coherent characteristics of their emission. In a companion work [Kader, Leung+2022], we use a novel interferometric method to search for gravitationally lensed FRBs in the time domain using bursts detected by CHIME/FRB. There, we dechannelize and autocorrelate electric field data at a time resolution of 1.25 ns. This enables a search for FRBs whose emission is coherently deflected by gravitational lensing around a foreground compact object such as a primordial black hole (PBH). Here, we use our non-detection of lensed FRBs to place novel constraints on the PBH abundance outside the Local Group. We use a novel two-screen model to take into account decoherence from scattering screens in our constraints. Our constraints are subject to a single astrophysical model parameter -- the effective distance between an FRB source and the scattering screen, for which we adopt a fiducial distance of 1 parsec. We find that coherent FRB lensing is a sensitive probe of sub-solar mass compact objects. Having observed no lenses in $172$ bursts from $114$ independent sightlines through the cosmic web, we constrain the fraction of dark matter made of compact objects, such as PBHs, to be $f \\lesssim 0.8$, if their masses are $\\sim 10^{-3} M_{\\odot}$.\n",
      "\n",
      "5. David Eichler, , arXiv:1711.04764\n",
      "\n",
      "Title: Nanolensed fast Radio Bursts\n",
      "\n",
      "Abstract: It is suggested that fast radio bursts can probe gravitational lensing by clumpy dark matter objects that range in mass from $10^{-3}M_{\\odot}$ to $10^2 M_{\\odot}$. They may provide a more sensitive probe than observations of lensing of objects in the Magellanic Clouds, and could find or rule out clumpy dark matter with an extended mass spectrum.\n",
      "\n"
     ]
    }
   ],
   "source": [
    "# abstract\n",
    "print_abstract_query_results(abstract_data, rst_abstract)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "id": "2f2c2828",
   "metadata": {},
   "outputs": [],
   "source": []
  },
  {
   "cell_type": "markdown",
   "id": "ef4e9d79",
   "metadata": {},
   "source": [
    "### example where citation query fails; abstract is way better\n",
    "I don't understand... there are clearly sentences in my citation database that has scintillation in it and openai embedding doesn't think those are strongly relevant?"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 34,
   "id": "ab1a09c3",
   "metadata": {},
   "outputs": [],
   "source": [
    "q = 'fast radio burst scintillation'"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 35,
   "id": "59018a9b",
   "metadata": {},
   "outputs": [],
   "source": [
    "rst_all_abstract = retriever_abstract.retrieve(q)\n",
    "rst_all_citation = retriever_citation.retrieve(q)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 36,
   "id": "61e5b068",
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAigAAAGdCAYAAAA44ojeAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjcuMywgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/OQEPoAAAACXBIWXMAAA9hAAAPYQGoP6dpAAAwtUlEQVR4nO3de3hU1aH+8XcyyUxCYJIGTAIHAojKRUArSBi1ysFUhJyWKtYbR6NSrRCsQmsBHxXFChztI9YexPYcLvYIpcVTtQWrIArWEkDSH62IUlE0VEiicJJILpNJZv3+SDNxJAFym7UnfD/Psx9X9qy1Zy33TuZlX9a4jDFGAAAADhJnuwMAAABfRUABAACOQ0ABAACOQ0ABAACOQ0ABAACOQ0ABAACOQ0ABAACOQ0ABAACOE2+7A20RCoV06NAh9ejRQy6Xy3Z3AADAKTDG6IsvvlCfPn0UF3ficyQxGVAOHTqkfv362e4GAABog4MHD6pv374nrBOTAaVHjx6SGgbo8/ks9wYA2qCyUurTp6F86JCUnGy3P3Cu9h4rDjrWKioq1K9fv/Dn+InEZEBpvKzj8/kIKABik9vdVPb5CChoWXuPFQcea6dyewY3yQIAAMeJyTMoABDz4uOlvLymMtCS9h4rMXqsuYwxxnYnWquiokIpKSkqLy/nEg8AADGiNZ/fXOIBAACOEzvnegCgKzFGqqpqKHfrJjGnE1rS3mMlRo81zqAAgA1VVVL37g1L44cH0Jz2HisxeqwRUAAAgOMQUAAAgOMQUAAAgOMQUAAAgOMQUAAAgOMQUAAAgOMwDwoA2OB2S9dc01QGWtLeYyVGjzWmugcAAFHRms9vzqAAgIMMmLvhpHU+XpwbhZ4AdnEPCgAAcBwCCgDYUFnZ8J0oLldDGWhJe4+VGD3WCCgAAMBxCCgAAMBxCCgAAMBxCCgAAMBxCCgAAMBxCCgAAMBxmKgNAGxwu6VJk5rKQEvae6zE6LFGQAEAGxITpQ0nnzUWaPexEqPHGpd4AACA4xBQAACA4xBQAMCGykopOblhiaHpx2FBe4+VGD3WuAcFAGypqrLdA8SK9h4rMXiscQYFAAA4DgEFAAA4DgEFAAA4DgEFAAA4DgEFAAA4Dk/xAIANcXHSZZc1lYGWtPdYidFjjYACADYkJUlbttjuBWJBe4+VGD3WYidKAQCA0wYBBQAAOA4BBQBsqKyUzjijYYmh6cdhQXuPlRg91rgHBQCiZMDcpq+8T6qt0Xuffy5JGvrAK6r2JNrqFmLBP48Va+0t4AwKAABwHAIKAABwHAIKAABwHAIKAABwHAIKAABwHJ7iAQALQi6X/pp5drgMtCguTho9uqkc7faWEFAAwIJAgleT85bY7gZiQVKS9Pbb9tpb0q4otXjxYrlcLt1zzz3hdTU1NcrPz1fPnj3VvXt3TZkyRSUlJRHtioqKlJubq27duik9PV333nuv6urq2tMVAADQhbQ5oLz99tv6xS9+oZEjR0asnzVrlv7whz9o3bp12rp1qw4dOqSrr746/Hp9fb1yc3NVW1urbdu26dlnn9WqVav04IMPtn0UAACgS2lTQDl27JimTp2q//qv/9LXvva18Pry8nItX75cTzzxhMaPH69Ro0Zp5cqV2rZtm7Zv3y5J2rhxo/bu3avnnntO559/viZOnKhHHnlES5cuVW1tbceMCgAcLjFYo7eW3aa3lt2mxGCN7e7AyaqqpAEDGpaqqui3t6RNASU/P1+5ubnKycmJWF9YWKhgMBixfsiQIcrKylJBQYEkqaCgQCNGjFBGRka4zoQJE1RRUaF333232fcLBAKqqKiIWAAglrmM1LeiVH0rSuUytnsDRzNG+uSThsW04WBpb3tLWn2T7Nq1a/WXv/xFbzdzw01xcbE8Ho9SU1Mj1mdkZKi4uDhc58vhpPH1xteas2jRIj388MOt7SoAAIhRrTqDcvDgQd19991avXq1EhOj98VW8+bNU3l5eXg5ePBg1N4bAABEX6sCSmFhoUpLS3XBBRcoPj5e8fHx2rp1q5566inFx8crIyNDtbW1Kisri2hXUlKizMxMSVJmZuZxT/U0/txY56u8Xq98Pl/EAgAAuq5WXeK5/PLL9c4770Ssu/XWWzVkyBDNmTNH/fr1U0JCgjZv3qwpU6ZIkvbt26eioiL5/X5Jkt/v16OPPqrS0lKlp6dLkjZt2iSfz6dhw4Z1xJgAIOoGzN1guwtAl9KqgNKjRw8NHz48Yl1ycrJ69uwZXj9t2jTNnj1baWlp8vl8uuuuu+T3+zV27FhJ0hVXXKFhw4bppptu0mOPPabi4mLdf//9ys/Pl9fr7aBhAQCAWNbhM8kuWbJEcXFxmjJligKBgCZMmKCnn346/Lrb7db69es1ffp0+f1+JScnKy8vTwsWLOjorgCAYxmX9PeeWeEy0CKXS2q8wtCWr0Vob3tLXMbE0DNH/1RRUaGUlBSVl5dzPwoAR4jmJZ6PF+dG7b2AjtSaz+/Y+dYgAABw2iCgAAAAxyGgAIAFicEabfzvGdr43zOY6h4nVlUlnXtuw9LWqe7b096SDr9JFgBwci4jnXOkKFwGWmSMtHdvUzna7S3hDAoAAHAcAgoAAHAcAgoAAHAcAgoAAHAcAgoAAHAcnuIBAAuMS/qHLz1cBlrkckn9+zeVo93eEgIKAFhQk5CoS6avsN0NxIJu3aSPP7bX3hIu8QAAAMchoAAAAMchoACABd5gQC89O0svPTtL3mDAdnfgZNXV0oUXNizV1dFvbwn3oACABXHG6LziD8JloEWhkLRrV1M52u0t4QwKAABwHAIKAABwHAIKAABwHO5BAYCTGDB3g+0uAKcdzqAAAADH4QwKAFhyJMlnuwuIFb162W1vAQEFACyo9iRq1A/W2O4GYkFysvTZZ/baW8IlHgAA4DgEFAAA4DgEFACwwBsMaO2auVq7Zi5T3ePEqqulceMalrZOdd+e9pZwDwoAWBBnjMYe3BMuAy0KhaStW5vK0W5vCWdQAACA4xBQAACA4xBQAACA43APCgDEmFOZev/jxblR6AnQeTiDAgAAHIczKABgSVWC13YXECu6dbPb3gICCgBYUO1J1LDZ/2u7G4gFyclSZaW99pZwiQcAADgOAQUAADgOAQUALPDW1WrFuoe0Yt1D8tbV2u4OnKymRsrNbVhqaqLf3hLuQQEAC+JCIY3/aFe4DLSovl56+eWmcrTbW8IZFAAA4DgEFAAA4DgEFAAA4DgEFAAA4DgEFAAA4DgEFAAA4Dg8ZgwAFlR7EjVgznrb3UAsSE6WjLHX3hLOoAAAAMchoAAAAMchoACABd66Wi19cZGWvriIqe5xYjU10ne/27C0dar79rS3hIACABbEhULK3fdn5e77M1Pd48Tq66Xnn29Y2jrVfXvaW0JAAQAAjkNAAQAAjkNAAQAAjkNAAQAAjkNAAQAAjkNAAQAAjsNU9wBgQXWCV0NnPR8uAy3q1k06dqypHO32lhBQAMAGl0vVnkTbvUAscLkavk/HVntLuMQDAAAch4ACABZ46oL66YYl+umGJfLUBW13B04WCEi33NKwBALRb28JAQUALHCH6nXNns26Zs9muUOxM/04LKirk559tmGpq4t+e0u4BwUAuqABczectM7Hi3Oj0BOgbTiDAgAAHIeAAgAAHIeAAgAAHIeAAgAAHIeAAgAAHIeneADAguoEry64a3W4DLSoWzeptLSpHO32lhBQAMAGl0tHu6XY7gVigcslnXGGvfaWcIkHAAA4TqsCyrJlyzRy5Ej5fD75fD75/X798Y9/DL9eU1Oj/Px89ezZU927d9eUKVNUUlISsY2ioiLl5uaqW7duSk9P17333qu6GJrZDgA6gqcuqAUbl2nBxmVMdY8TCwSk/PyGpa1T3benvSWtCih9+/bV4sWLVVhYqF27dmn8+PGaPHmy3n33XUnSrFmz9Ic//EHr1q3T1q1bdejQIV199dXh9vX19crNzVVtba22bdumZ599VqtWrdKDDz7YsaMCAIdzh+p18//boJv/3wamuseJ1dVJTz/dsLR1qvv2tLekVfegfOtb34r4+dFHH9WyZcu0fft29e3bV8uXL9eaNWs0fvx4SdLKlSs1dOhQbd++XWPHjtXGjRu1d+9evfbaa8rIyND555+vRx55RHPmzNFDDz0kj8fTcSMDAAAxq833oNTX12vt2rWqrKyU3+9XYWGhgsGgcnJywnWGDBmirKwsFRQUSJIKCgo0YsQIZWRkhOtMmDBBFRUV4bMwzQkEAqqoqIhYAABA19XqgPLOO++oe/fu8nq9uvPOO/XCCy9o2LBhKi4ulsfjUWpqakT9jIwMFRcXS5KKi4sjwknj642vtWTRokVKSUkJL/369WtttwEAQAxpdUAZPHiwdu/erR07dmj69OnKy8vT3r17O6NvYfPmzVN5eXl4OXjwYKe+HwAAsKvV86B4PB6dddZZkqRRo0bp7bff1s9+9jNdd911qq2tVVlZWcRZlJKSEmVmZkqSMjMztXPnzojtNT7l01inOV6vV14vExkBAHC6aPc8KKFQSIFAQKNGjVJCQoI2b94cfm3fvn0qKiqS3++XJPn9fr3zzjsqbZzRTtKmTZvk8/k0bNiw9nYFAAB0Ea06gzJv3jxNnDhRWVlZ+uKLL7RmzRpt2bJFr776qlJSUjRt2jTNnj1baWlp8vl8uuuuu+T3+zV27FhJ0hVXXKFhw4bppptu0mOPPabi4mLdf//9ys/P5wwJgNNKTYJHl9y5PFwGWpSUJB040FSOdntLWhVQSktLdfPNN+vw4cNKSUnRyJEj9eqrr+qb3/ymJGnJkiWKi4vTlClTFAgENGHCBD399NPh9m63W+vXr9f06dPl9/uVnJysvLw8LViwoGNHBQAOZ1xx+kdKxskrAnFx0oAB9tpb4jLGGNudaK2KigqlpKSovLxcPp/PdncAdHED5m6w3YVO8fHiXNtdwGmmNZ/ffBcPAFiQUB/UvDdWaN4bK5RQz1T3OIHaWuneexuW2trot7eEgAIAFsTX1+v7O3+n7+/8neLrmeoeJxAMSj/9acMSbEOYbW97S1r9mDEAdCVd9fINEOsIKACAFp1KgONeFnQGLvEAAADHIaAAAADHIaAAAADHIaAAAADH4SZZALCgJsGjb962NFy2gSeYYkRSkrRnT1M52u0tIaAAgAXGFacPzuhvuxuIBXFx0rnn2mtvCZd4AACA43AGBQAsSKgPKr/gt5Kkpf5rFXQnWO4RHKu2Vlq4sKF8332Sp5WXBNvb3hK+LBDAac3WfRhJtTV6b8k1kqShs55XtSfRSj86AhO1dbLKSql794bysWNScnJ023cgviwQAADENAIKAABwHAIKAABwHAIKAABwHAIKAABwHAIKAABwHOZBAQALAvEJ+vbNT4TLQIsSE6WdO5vK0W5vCQEFACwIxbn1t97n2O4GYoHbLV14ob32lnCJBwAAOA5nUADAgoT6oG7d9XtJ0srR32aqe7Sstlb62c8aynff3bap7tvT3hICCgBYEF9fr/u2rJQk/c/XcwkoaFkwKP34xw3lGTNaHzDa294SLvEAAADHIaAAAADHIaAAAADHIaAAAADHIaAAAADHIaAAAADH4TFjALAgEJ+g629YGC4DLUpMlN54o6kc7faWEFAAwIJQnFvbs0ba7gZigdstjRtnr70lBBQAMWnA3A22uwCgExFQAMCC+Po63fDXVyRJvz7vStW5+XOMFgSD0i9/2VC+4w4poZWXBNvb3hJ+IwDAgoT6Oj2y6RlJ0vPDcwgoaFltrTRzZkP5lltaHzDa294SnuIBAACOQ0ABAACOQ0ABAACOQ0ABAACOQ0ABAACOQ0ABAACOw3NtAGBBbXyCbr1mfrgMtMjrldavbypHu70lBBQAsKA+zq03Bl1ouxuIBfHxUm6uvfaWcIkHAAA4DmdQAMCC+Po6fWfvFknSi8PGMZMsWhYMSqtXN5SnTm3bVPftaW8JvxEAYEFCfZ1++vKTkqQNgy8hoKBltbXSrbc2lL/73bZNdd+e9pZwiQcAADgOAQUAADgOAQUAADgOAQUAADgOAQUAADgOAQUAADgOz7UBgAW18QmaMXluuAy0yOuVfvvbpnK021tCQAHgOAPmbrDdhU5XH+fWy0Musd0NxIL4+Ib5S2y1t4RLPAAAwHE4gwIAFrhD9Zrw9wJJ0qvn+FUf57bcIzhWXZ30wgsN5auuajgjEs32lsRGLwGgi/HUBfX0S4slSUNnPa9qDwEFLQgEpGuvbSgfO9b6gNHe9pZwiQcAADgOAQUAADgOAQUAADgOAQUAADgOAQUAADgOAQUAADhObDxrBABdTNAdrx9NuidcBlrk8UgrVzaVo93eEn4rAMCCOne8nh+RY7sbiAUJCdItt9hrbwmXeAAAgONwBgUALHCH6nXpgb9Ikt4ceAFT3aNldXXSq682lCdMaNtU9+1pb0ls9BIAuhhPXVArn39YElPd4yQCAenf/q2h3Nap7tvT3pJW9XLRokX63e9+p/fff19JSUm66KKL9B//8R8aPHhwuE5NTY1++MMfau3atQoEApowYYKefvppZWRkhOsUFRVp+vTpeuONN9S9e3fl5eVp0aJFio+R/2kA2m7A3A22uwAgBrQqEWzdulX5+fm68MILVVdXp/vuu09XXHGF9u7dq+TkZEnSrFmztGHDBq1bt04pKSmaOXOmrr76av35z3+WJNXX1ys3N1eZmZnatm2bDh8+rJtvvlkJCQlauHBhx48QANCpTiV0frw4Nwo9QVfSqoDyyiuvRPy8atUqpaenq7CwUJdeeqnKy8u1fPlyrVmzRuPHj5ckrVy5UkOHDtX27ds1duxYbdy4UXv37tVrr72mjIwMnX/++XrkkUc0Z84cPfTQQ/LE0CNQAACgc7TrKZ7y8nJJUlpamiSpsLBQwWBQOTlNj84NGTJEWVlZKigokCQVFBRoxIgREZd8JkyYoIqKCr377rvNvk8gEFBFRUXEAgAAuq42B5RQKKR77rlHF198sYYPHy5JKi4ulsfjUWpqakTdjIwMFRcXh+t8OZw0vt74WnMWLVqklJSU8NKvX7+2dhsAAMSANgeU/Px87dmzR2vXru3I/jRr3rx5Ki8vDy8HDx7s9PcEAAD2tOmxmZkzZ2r9+vV688031bdv3/D6zMxM1dbWqqysLOIsSklJiTIzM8N1du7cGbG9kpKS8GvN8Xq98nq9bekqADhS0B2vB755Z7gMtMjjkf7zP5vK0W5vSat+K4wxuuuuu/TCCy9oy5YtGjhwYMTro0aNUkJCgjZv3qwpU6ZIkvbt26eioiL5/X5Jkt/v16OPPqrS0lKlp6dLkjZt2iSfz6dhw4Z1xJgAwPHq3PH6nwv+zXY3EAsSEqT8fHvtLWlVQMnPz9eaNWv00ksvqUePHuF7RlJSUpSUlKSUlBRNmzZNs2fPVlpamnw+n+666y75/X6NHTtWknTFFVdo2LBhuummm/TYY4+puLhY999/v/Lz8zlLAgAAJLUyoCxbtkySNG7cuIj1K1eu1C3//CKiJUuWKC4uTlOmTImYqK2R2+3W+vXrNX36dPn9fiUnJysvL08LFixo30gAIIbEheo15h8NTy7u7HuuQkx1j5bU10t/+lND+RvfkNytPFba294SlzHG2O5Ea1VUVCglJUXl5eXy+Xy2uwOgFZhJtkFSbY3eW3KNpMap7hMt96hzMVFbO1RWSt27N5SPHZP+OTFq1Np3oNZ8fvNtxgAAwHEIKAAAwHEIKAAAwHEIKAAAwHEIKAAAwHEIKAAAwHGYXxkALKhzu7Vw3K3hMtCihATpsceaytFubwkBBQAsCLoT9MvsKba7gVjg8Uj33muvvSVc4gEAAI7DGRQAsCAuVK/hJR9KkvZkDGKqe7Ssvl76y18ayhdc0Lap7tvT3hICCgBY4K0L6ve/mi2pcar72PjQgAU1NdKYMQ3ltkxV3972lnCJBwAAOA4BBQAAOA4BBQAAOA4BBQAAOA4BBQAAOA4BBQAAOA6PGQOABXVut568+IZwGWhRQoI0f35TOdrtLXEZY4ztTrRWRUWFUlJSVF5eLp/PZ7s7AFphwNwNtrsACz5enGu7C3CA1nx+c4kHAAA4Dpd4AMAClwnprM8PSpL29+on4+Lfi2hBKCS9915DeehQKa6Vx0p721tCQAEACxKDtdq0Il9S41T3iZZ71LlO5dIel4FaUF0tDR/eUG7LVPXtbW8JAQUA0KUQhroGAgqADsMNsAA6SmxciAIAAKcVAgoAAHAcAgoAAHAcAgoAAHAcbpIFAAvq3G79YszV4TLQooQE6Uc/aipHu70lBBQAsCDoTtCif73NdjcQCzwe6fHH7bW3hEs8AADAcTiDAgAWuExI/1LxmSTpU98ZTHWPloVCUlFRQzkrq21T3benvSUEFACwIDFYq7eemSbp9JjqHu1QXS0NHNhQbutU9+1pbwkBBQDgCExRjy+LjfM8AADgtEJAAQAAjkNAAQAAjkNAAQAAjsNNsgCAmHEqN9KiayCgAIAF9XFu/errueEy0KL4eGnGjKZytNtbEjs9BYAupDY+QQ9eMd12NxALvF5p6VJ77S3hHhQAAOA4nEEBABuMUVp1hSTpaJJPcrksdwiOZYz0+ecN5V69Wn+stLe9JQQUALAgKRjQX34+VRJT3eMkqqqk9PSGclumqm9ve0u4xAMAAByHgAIAAByHgAIAAByHgAIAAByHgAIAAByHgAIAAByHx4wBwIL6OLeeH355uAy0KD5eystrKke7vSWx01MA6EJq4xP0o9xZtruBWOD1SqtW2WtvCQEFwCnhW2QBRBMBBQBsMEZJwYAkqTrBGzPTj8MCYxpmg5Wkbt3aNtV9e9pbwk2yAGBBUjCg95Zco/eWXBMOKkCzqqqk7t0blsagEc32lhBQAACA4xBQAACA4xBQAACA4xBQAACA4xBQAACA4xBQAACA4zAPCgBYEIqL04bBF4fLQIvcbumaa5rK0W5vCQEFACwIxHuU/515truBWJCYKK1bZ6+9JcR2AADgOAQUAADgOFziAcAXAVqQVFuj95Y03BcwdNbzqvYkWu4RHKuysmGaekk6dkxKTo5ue0s4gwIAAByn1QHlzTff1Le+9S316dNHLpdLL774YsTrxhg9+OCD6t27t5KSkpSTk6MPPvggos7Ro0c1depU+Xw+paamatq0aTp27Fi7BgIAALqOVgeUyspKnXfeeVq6dGmzrz/22GN66qmn9Mwzz2jHjh1KTk7WhAkTVFNTE64zdepUvfvuu9q0aZPWr1+vN998U3fccUfbRwEAALqUVt+DMnHiRE2cOLHZ14wxevLJJ3X//fdr8uTJkqRf/epXysjI0Isvvqjrr79e7733nl555RW9/fbbGj16tCTp5z//uSZNmqSf/vSn6tOnTzuGA5xeTuXekY8X50ahJwDQsTr0HpQDBw6ouLhYOTk54XUpKSnKzs5WQUGBJKmgoECpqanhcCJJOTk5iouL044dOzqyOwAAIEZ16FM8xcXFkqSMjIyI9RkZGeHXiouLlZ6eHtmJ+HilpaWF63xVIBBQIBAI/1xRUdGR3QYAAA4TE48ZL1q0SA8//LDtbgAxiUeInSkUF6fXzxwdLgMtcrulSZOaytFub0mHBpTMzExJUklJiXr37h1eX1JSovPPPz9cp7S0NKJdXV2djh49Gm7/VfPmzdPs2bPDP1dUVKhfv34d2XUAiKpAvEe3ffch291ALEhMlDa04x8a7W1vSYfG9oEDByozM1ObN28Or6uoqNCOHTvk9/slSX6/X2VlZSosLAzXef311xUKhZSdnd3sdr1er3w+X8QCAAC6rlafQTl27Jj2798f/vnAgQPavXu30tLSlJWVpXvuuUc/+clPdPbZZ2vgwIF64IEH1KdPH33nO9+RJA0dOlRXXnmlbr/9dj3zzDMKBoOaOXOmrr/+ep7gAQBEBU/AOV+rA8quXbv0r//6r+GfGy+95OXladWqVfrxj3+syspK3XHHHSorK9Mll1yiV155RYmJTdM4r169WjNnztTll1+uuLg4TZkyRU899VQHDAcAYkNSbY0K/3OqJGnUzNVMdY+WVVZKjQ+XlJa2bar79rS3pNUBZdy4cTLGtPi6y+XSggULtGDBghbrpKWlac2aNa19awDoUroFAyevBEhSVZXd9hZw6zgAAHAcAgoAAHCcmJgHBQCAWMUNuW3DGRQAAOA4BBQAAOA4XOIBAAtCLpe29xseLgMtiouTLrusqRzt9pYQUADAgkCCV9ffuNh2NxALkpKkLVvstbeEgAIAQDO4udWu2DnXAwAAThsEFACwIKm2RoVP3ajCp25UUm2N7e7AySorpTPOaFgqK6Pf3hIu8QCAJT2rK2x3AbHi88/ttreAMygAAMBxOIMCAEAbncqNtGgbzqAAAADH4QwKAACWnehMTFJtjd6LYl+cgjMoAADAcTiDAgAWhFwu/TXz7HAZaEnjsXJe35S2T3U/enRTOUa4jDHGdidaq6KiQikpKSovL5fP57PdHaBTcPMdgC/rCrPWtubzO3aiFAAAOG0QUAAAgOMQUADAgsRgjd5adpveWnabEoNMdY+WNR4rGjBAqqpq/QaqqhratrW9JdwkCwAWuIzUt6I0XAZaEj5WKiS15bZRY6RPPmkqxwjOoAAAAMchoAAAAMchoAAAAMchoAAAAMfhJlnAAiZhA4ATI6AAgAXGJf29Z1a4DLSk8Vg5J6O71JavRXC5pGHDmsoxgqnuAQs4gwKgtZjqHgAAwDICCgAAcBwCCgBYkBis0cb/nqGN/z2Dqe5xQo3His49t+1T3Z97btvbW8JNskAH4/4SnAqXkc45UhQuAy0JHytH1Pap7vfubSrHCM6gAAAAxyGgAAAAxyGgAAAAxyGgAAAAxyGgAAAAx+EpHgCwwLikf/jSw2WgJY3HSt+vJbV9qvv+/ZvKMYKp7oEOxmPGADoDU90DAABYRkABAACOQ0ABAAu8wYBeenaWXnp2lrzBgO3uwMEajxVdeKFUXd36DVRXN7Rta3tLuEkWACyIM0bnFX8QLgMtCR8rxZJCodZvIBSSdu1qKscIzqAAAADHIaAAAADH4RIP8E88HgwAzsEZFAAA4DicQQEAIEYMfeAVVXsST1inK0zoJhFQAMCaI0nMhI1T0+5jpVevjulIFBFQAMCCak+iRv1gje1uIAa0+1hJTpY++6zjOhQlBBScFrgBFgBiCzfJAgAAxyGgAIAF3mBAa9fM1do1c5nqHifU7mOluloaN65hYap7AMCJxBmjsQf3hMtAS9p9rIRC0tatTeUYwRkUAADgOAQUAADgOFzigaOdytM3XWVSIgBAE86gAAAAxyGgAAAAx+ESDwBYUpXgtd0FxIh2HyvdukX8GAuXzwkoAGBBtSdRw2b/r+1uIAa09lhpNnzc9duG/z6ypUP6FA0EFMQ8prEHgK6HgIJOEQunDwEAzkVAAQALvHW1WvbCQknS9KvuUyDeY7lHcKr2HiuxeqwRUADAgrhQSOM/2hUuAy1p77ESq8caAQXWcO8IAKAlBBS0GsECANDZrE7UtnTpUg0YMECJiYnKzs7Wzp07bXYHAAA4hLUzKL/5zW80e/ZsPfPMM8rOztaTTz6pCRMmaN++fUpPT7fVLUfiiRgAwOnGWkB54okndPvtt+vWW2+VJD3zzDPasGGDVqxYoblz59rq1mmPyzcAACewElBqa2tVWFioefPmhdfFxcUpJydHBQUFx9UPBAIKBALhn8vLyyVJFRUVndK/4fNf7ZTtNmfPwxNOWicUqDppnaxZ6zqiOwCipL62Ro1/weoDVQqZ2Hm6AtHV3mOlre074zO2cZvGmJNXNhZ8+umnRpLZtm1bxPp7773XjBkz5rj68+fPN5JYWFhYWFhYusBy8ODBk2aFmHiKZ968eZo9e3b451AopKNHj6pnz55yuVwttquoqFC/fv108OBB+Xy+aHQ1ahhb7OrK42NssYmxxa5YG58xRl988YX69Olz0rpWAkqvXr3kdrtVUlISsb6kpESZmZnH1fd6vfJ6I7/JMTU19ZTfz+fzxcSOawvGFru68vgYW2xibLErlsaXkpJySvWsPGbs8Xg0atQobd68ObwuFApp8+bN8vv9NroEAAAcxNolntmzZysvL0+jR4/WmDFj9OSTT6qysjL8VA8AADh9WQso1113nT777DM9+OCDKi4u1vnnn69XXnlFGRkZHfYeXq9X8+fPP+7yUFfA2GJXVx4fY4tNjC12deXxuYw5lWd9AAAAosfqVPcAAADNIaAAAADHIaAAAADHIaAAAADHcXRAWbp0qQYMGKDExERlZ2dr586dJ6xfVlam/Px89e7dW16vV+ecc45efvnl8OsDBgyQy+U6bsnPzw/XGTdu3HGv33nnnVbH1lyfXC6XcnObvsHYGKMHH3xQvXv3VlJSknJycvTBBx9EbOfo0aOaOnWqfD6fUlNTNW3aNB07dqzDx9bR4wsGg5ozZ45GjBih5ORk9enTRzfffLMOHToUsZ3m9u/ixYsdPTZJuuWWW457/corr4zYTrT2XUePrbnXXS6XHn/88XAdJ+43SXryySc1ePBgJSUlqV+/fpo1a5Zqampatc2amhrl5+erZ8+e6t69u6ZMmXLcBJVOHNuiRYt04YUXqkePHkpPT9d3vvMd7du3L2Ib0fpb2Rnje+ihh47r+5AhQyK2Eav7zkmfc+3WIV+u0wnWrl1rPB6PWbFihXn33XfN7bffblJTU01JSUmz9QOBgBk9erSZNGmSeeutt8yBAwfMli1bzO7du8N1SktLzeHDh8PLpk2bjCTzxhtvhOtcdtll5vbbb4+oV15ebnVsR44ciejPnj17jNvtNitXrgzXWbx4sUlJSTEvvvii+etf/2q+/e1vm4EDB5rq6upwnSuvvNKcd955Zvv27eZPf/qTOeuss8wNN9zQoWPrjPGVlZWZnJwc85vf/Ma8//77pqCgwIwZM8aMGjUqYjv9+/c3CxYsiNjWsWPHHD02Y4zJy8szV155ZUS9o0ePRmwnGvuuM8b25dcPHz5sVqxYYVwul/nwww/DdZy431avXm28Xq9ZvXq1OXDggHn11VdN7969zaxZs1q1zTvvvNP069fPbN682ezatcuMHTvWXHTRRY4f24QJE8zKlSvNnj17zO7du82kSZNMVlZWxH6Jxt/Kzhrf/PnzzbnnnhvR988++yxiO7G675zyOdcRHBtQxowZY/Lz88M/19fXmz59+phFixY1W3/ZsmXmzDPPNLW1taf8HnfffbcZNGiQCYVC4XWXXXaZufvuu9vc71PR2rF91ZIlS0yPHj3CfyxCoZDJzMw0jz/+eLhOWVmZ8Xq95te//rUxxpi9e/caSebtt98O1/njH/9oXC6X+fTTTztiWGEdPb7m7Ny500gyn3zySXhd//79zZIlS9rc71PRGWPLy8szkydPbrFNtPZdNPbb5MmTzfjx4yPWOXG/5efnH9fP2bNnm4svvviUt1lWVmYSEhLMunXrwnXee+89I8kUFBR0yLhOpR9fdSpj+6rS0lIjyWzdujW8Lhp/K43pnPHNnz/fnHfeeS2+Z1fad7Y+5zqCIy/x1NbWqrCwUDk5OeF1cXFxysnJUUFBQbNtfv/738vv9ys/P18ZGRkaPny4Fi5cqPr6+hbf47nnntNtt9123BcOrl69Wr169dLw4cM1b948VVVVWR3bVy1fvlzXX3+9kpOTJUkHDhxQcXFxxDZTUlKUnZ0d3mZBQYFSU1M1evTocJ2cnBzFxcVpx44dHTE0SZ0zvuaUl5fL5XId951MixcvVs+ePfX1r39djz/+uOrq6to0juZ05ti2bNmi9PR0DR48WNOnT9eRI0fCr0Vj30Vjv5WUlGjDhg2aNm3aca85bb9ddNFFKiwsDJ9u/+ijj/Tyyy9r0qRJp7zNwsJCBYPBiDpDhgxRVlbWKf8/tTG25pSXl0uS0tLSItZ35t9KqXPH98EHH6hPnz4688wzNXXqVBUVFYVf6yr7ztbnXEdx5LcZf/7556qvrz9uVtmMjAy9//77zbb56KOP9Prrr2vq1Kl6+eWXtX//fs2YMUPBYFDz588/rv6LL76osrIy3XLLLRHrb7zxRvXv3199+vTR3/72N82ZM0f79u3T7373O2tj+7KdO3dqz549Wr58eXhdcXFxeBtf3Wbja8XFxUpPT494PT4+XmlpaeE6HaEzxvdVNTU1mjNnjm644YaIL8f6wQ9+oAsuuEBpaWnatm2b5s2bp8OHD+uJJ55o+4C+pLPGduWVV+rqq6/WwIED9eGHH+q+++7TxIkTVVBQILfbHZV9F4399uyzz6pHjx66+uqrI9Y7cb/deOON+vzzz3XJJZfIGKO6ujrdeeeduu+++055m8XFxfJ4PMeF6C//XjpxbF8VCoV0zz336OKLL9bw4cMjttOZfys7c3zZ2dlatWqVBg8erMOHD+vhhx/WN77xDe3Zs0c9evToMvvO1udcR3FkQGmLUCik9PR0/fKXv5Tb7daoUaP06aef6vHHH282oCxfvlwTJ0487iuf77jjjnB5xIgR6t27ty6//HJ9+OGHGjRoUKeP42SWL1+uESNGaMyYMba70ilONr5gMKhrr71WxhgtW7Ys4rXZs2eHyyNHjpTH49H3v/99LVq0yBHTQLc0tuuvvz5cHjFihEaOHKlBgwZpy5Ytuvzyy6PdzTY5leNyxYoVmjp1qhITEyPWO3G/bdmyRQsXLtTTTz+t7Oxs7d+/X3fffbceeeQRPfDAA1b61FFaO7b8/Hzt2bNHb731VsR6p/6tPJXxTZw4MVx/5MiRys7OVv/+/fXb3/622TN8TtHafRern3ONHHmJp1evXnK73cfdMV1SUqLMzMxm2/Tu3VvnnHOO3G53eN3QoUNVXFys2traiLqffPKJXnvtNX3ve987aV+ys7MlSfv372/tMJrVlrE1qqys1Nq1a4/7BWpsd6JtZmZmqrS0NOL1uro6HT169KTv2xqdMb5GjeHkk08+0aZNm0761eLZ2dmqq6vTxx9/3KoxtKQzx/ZlZ555pnr16hU+5qKx7zp7bH/605+0b9++U/6ds73fHnjgAd1000363ve+pxEjRuiqq67SwoULtWjRIoVCoVPaZmZmpmpra1VWVnbK7+uEsX3ZzJkztX79er3xxhvq27fvCfvS0X8rpc4fX6PU1FSdc845Eb9zsb7vbH7OdRRHBhSPx6NRo0Zp8+bN4XWhUEibN2+W3+9vts3FF1+s/fv3R+ykv//97+rdu7c8Hk9E3ZUrVyo9PT3icciW7N69W1JDAOoIbRlbo3Xr1ikQCOjf//3fI9YPHDhQmZmZEdusqKjQjh07wtv0+/0qKytTYWFhuM7rr7+uUCgUPjg7QmeMT2oKJx988IFee+019ezZ86R92b17t+Li4o67PNJWnTW2r/rHP/6hI0eOhI+5aOy7zh7b8uXLNWrUKJ133nkn7YsT9ltVVZXi4iL/PDb+48cYc0rbHDVqlBISEiLq7Nu3T0VFRSf9f2pzbI3/nTlzpl544QW9/vrrGjhw4En70tF/K6XOG99XHTt2TB9++GG477G87xrZ/JzrMLbuzj2ZtWvXGq/Xa1atWmX27t1r7rjjDpOammqKi4uNMcbcdNNNZu7cueH6RUVFpkePHmbmzJlm3759Zv369SY9Pd385Cc/idhufX29ycrKMnPmzDnuPffv328WLFhgdu3aZQ4cOGBeeuklc+aZZ5pLL73U6tgaXXLJJea6665rdpuLFy82qamp5qWXXjJ/+9vfzOTJk5t9zPjrX/+62bFjh3nrrbfM2Wef3WmPGXfk+Gpra823v/1t07dvX7N79+6IR+MCgYAxxpht27aZJUuWmN27d5sPP/zQPPfcc+aMM84wN998s6PH9sUXX5gf/ehHpqCgwBw4cMC89tpr5oILLjBnn322qampCdeLxr7rjOPSGGPKy8tNt27dzLJly457zan7bf78+aZHjx7m17/+tfnoo4/Mxo0bzaBBg8y11157yts0puFR1aysLPP666+bXbt2Gb/fb/x+v+PHNn36dJOSkmK2bNkS8ftWVVVljIne38rOGt8Pf/hDs2XLFnPgwAHz5z//2eTk5JhevXqZ0tLScJ1Y3XfGOONzriM4NqAYY8zPf/5zk5WVZTwejxkzZozZvn17+LXLLrvM5OXlRdTftm2byc7ONl6v15x55pnm0UcfNXV1dRF1Xn31VSPJ7Nu377j3KyoqMpdeeqlJS0szXq/XnHXWWebee+/tlOfDWzu2999/30gyGzdubHZ7oVDIPPDAAyYjI8N4vV5z+eWXHzfGI0eOmBtuuMF0797d+Hw+c+utt5ovvviiw8dmTMeO78CBA0ZSs0vjs/2FhYUmOzvbpKSkmMTERDN06FCzcOHCiA95J46tqqrKXHHFFeaMM84wCQkJpn///ub222+P+JAzJnr7rqOPS2OM+cUvfmGSkpJMWVnZca85db8Fg0Hz0EMPmUGDBpnExETTr18/M2PGDPN///d/p7xNY4yprq42M2bMMF/72tdMt27dzFVXXWUOHz7s+LG19PvWOMdNNP9Wdsb4rrvuOtO7d2/j8XjMv/zLv5jrrrvO7N+/P+I9Y3XfGeOcz7n2chnTwjkvAAAASxx5DwoAADi9EVAAAIDjEFAAAIDjEFAAAIDjEFAAAIDjEFAAAIDjEFAAAIDjEFAAAIDjEFAAAIDjEFAAAIDjEFAAAIDjEFAAAIDj/H/UnFE2+4gcsQAAAABJRU5ErkJggg==",
      "text/plain": [
       "<Figure size 640x480 with 1 Axes>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "scores = get_all_scores(rst_all_citation)\n",
    "score_threshold = np.percentile(scores, 99.7)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 37,
   "id": "8218b7af",
   "metadata": {},
   "outputs": [],
   "source": [
    "rst_abstract = rearrange_query_results(rst_all_abstract, score_threshold=0.8)\n",
    "rst_citation = rearrange_query_results(rst_all_citation, score_threshold=score_threshold)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 38,
   "id": "1d8203f7",
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "1. lorimer et al. 2007, arXiv:0709.4301\n",
      "\n",
      "Reasons of citation: referenced as the first discovery of fast radio bursts (frbs).; referenced as the first discovery of fast radio bursts (frbs).; referenced as the first discovery of fast radio bursts (frbs).; referenced to support the description of fast radio bursts (frbs) as bright, millisecond duration flashes of coherent radio emission.; referenced to support the general description of fast radio bursts (frbs) as short, bright flares of coherent radio emission.; these references are cited to support the statement that fast radio bursts (frbs) are bright, millisecond duration flares of coherent radio emission with large dispersion measures (dms) pointing to an extragalactic origin.; referenced as the original discovery of fast radio bursts (frbs).\n",
      "\n",
      "2. thornton et al. 2013, arXiv:1307.1628\n",
      "\n",
      "Reasons of citation: referenced to support the description of fast radio bursts (frbs) as bright, millisecond duration flashes of coherent radio emission.; referenced to support the general description of fast radio bursts (frbs) as short, bright flares of coherent radio emission.; these references are cited to support the statement that fast radio bursts (frbs) are bright, millisecond duration flares of coherent radio emission with large dispersion measures (dms) pointing to an extragalactic origin.\n",
      "\n",
      "3. cordes & chatterjee 2019, arXiv:1906.05878\n",
      "\n",
      "Reasons of citation: referenced to support the general description of fast radio bursts (frbs) as short, bright flares of coherent radio emission.; these references are cited to establish the significance of fast radio bursts (frbs) as the brightest flashes in the radio sky.\n",
      "\n",
      "4. petroff et al. 2022; petroff et al. 2021, arXiv:2107.10113\n",
      "\n",
      "Reasons of citation: referenced as a review source for fast radio bursts (frbs) and their intense, millisecond-duration radio flashes from extragalactic distances.\n",
      "\n",
      "5. petroff et al. 2016, arXiv:1601.03547\n",
      "\n",
      "Reasons of citation: referenced to support the description of fast radio bursts (frbs) as bright, millisecond duration flashes of coherent radio emission.\n",
      "\n"
     ]
    }
   ],
   "source": [
    "# citation\n",
    "print_citation_query_results(citation_data, rst_citation)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 39,
   "id": "d3530232",
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "1. Qiwu Song et al., arXiv:1707.08979\n",
      "\n",
      "Title: The radiation mechanism of fast radio bursts\n",
      "\n",
      "Abstract: Fast radio bursts are radio transients observed mainly around 1.5 GHz. Their peak frequency decreases at a rate of 100 ~ 500 MHz/s and some of them have a broader pulse with an exponentially decaying tail. Common assumptions for fast radio bursts include a dispersion effect resulting in the peak frequency drifting and a scattering effect resulting in pulse broadening. These assumptions attribute the abnormally large dispersion measure and scattering measure to the environmental medium of the host galaxy. Here we show that the radiation of fast radio bursts can be explained as an undulator radiation and the large dispersion measure can be due to a motion effect mainly from the rotation of the source which is probably variable stars. In our scenario, the pulse broadening is near-field effects and the pulse itself represents a Fresnel diffraction pattern sweeping the observer. Our work is the first analysis of properties of fast radio bursts in the context of a special mechanism of the radiation instead of a special propagation environment of the radiation.\n",
      "\n",
      "2. R. A. Main et al., arXiv:2108.00052\n",
      "\n",
      "Title: Scintillation timescale measurement of the highly active FRB20201124A\n",
      "\n",
      "Abstract: Scintillation of compact radio sources results from the interference between images caused by multipath propagation, and probes the intervening scattering plasma and the velocities of the emitting source and scattering screen. In FRB20201124A, a repeating fast radio burst (FRB) which entered a period of extreme activity, we obtained many burst detections in observations at the upgraded Giant Metrewave Radio Telescope (uGMRT) and the Effelsberg 100\\,m Radio Telescope. Bursts nearby in time show similar scintillation patterns, and we measure a scintillation timescale of $14.3\\pm1.2$\\,min and $7\\pm2$\\,min at Effelsberg (1370\\,MHz) and uGMRT (650\\,MHz), respectively, by correlating burst pair spectra. The scintillation bandwidth scaled to 1\\,GHz is $0.5\\pm0.1$\\,MHz, and the inferred scintillation velocity at Effelsberg is $V_{\\mathrm{ISS}}\\approx (59\\pm7) \\sqrt{d_{l}/2\\,\\rm{kpc}}~{\\rm km~s}^{-1}$, higher than Earth's velocity for any screen beyond a lens distance of $d_{l} \\gtrsim 400\\,$pc. From the measured scintillation bandwidth, FRB20201124A has comparatively lower scattering than nearby pulsars, and is underscattered by a factor of $\\sim 30$ or $\\sim 1200$ compared to the NE2001 and YMW16 model predictions respectively. This underscattering, together with the measured scintillation velocity are consistent with a scattering screen more nearby the Earth at $d_{l} \\sim 400\\,$pc, rather than at the 2\\,kpc spiral arm which NE2001 predicts to be the dominant source of scattering. With future measurements, the distance, geometry, and velocity of the scattering screen could be obtained through modelling of the annual variation in $V_{\\rm ISS}$, or through inter-station time delays or interferometric observations. Scintillation/scattering measurements of FRBs could help improve Galactic electron density models, particularly in the Galactic halo or at high Galactic latitudes.\n",
      "\n",
      "3. Zi-Wei Wu et al., arXiv:2304.14697\n",
      "\n",
      "Title: Scintillation Arc from FRB 20220912A\n",
      "\n",
      "Abstract: We present the interstellar scintillation analysis of fast radio burst (FRB) 20220912A during its extremely active episode in 2022 using data from the Five-hundred-meter Aperture Spherical Radio Telescope (FAST). We detect a scintillation arc in the FRB's secondary spectrum, which describes the power in terms of the scattered FRB signals' time delay and Doppler shift. The arc indicates that the scintillation is caused by a highly localized region of the ionized interstellar medium (IISM). Our analysis favors a Milky Way origin for the localized scattering medium but cannot rule out a host galaxy origin. We present our method for detecting the scintillation arc, which can be applied generally to sources with irregularly spaced bursts or pulses. These methods could help shed light on the complex interstellar environment surrounding the FRBs and in our Galaxy.\n",
      "\n",
      "4. J. I. Katz, , arXiv:1604.01799\n",
      "\n",
      "Title: Fast Radio Bursts---A Brief Review: Some Questions, Fewer Answers\n",
      "\n",
      "Abstract: Fast Radio Bursts are millisecond bursts of radio radiation at frequencies of about 1 GHz, recently discovered in pulsar surveys. They have not yet been definitively identified with any other astronomical object or phenomenon. The bursts are strongly dispersed, indicating passage through a high column density of low density plasma. The most economical interpretation is that this is the interglactic medium, indicating that FRB are at \"cosmological\" distances with redshifts in the range 0.3--1.3. Their inferred brightness temperatures are as high as $10^{37\\,\\circ}$K, implying coherent emission by \"bunched\" charges, as in radio pulsars. I review the astronomical sites, objects and emission processes that have been proposed as the origin of FRB, with particular attention to Soft Gamma Repeaters and giant pulsar pulses.\n",
      "\n",
      "5. Laura G Spitler et al., arXiv:1807.03722\n",
      "\n",
      "Title: Detection of Bursts from FRB 121102 with the Effelsberg 100-m Radio Telescope at 5 GHz and the Role of Scintillation\n",
      "\n",
      "Abstract: FRB 121102, the only repeating fast radio burst (FRB) known to date, was discovered at 1.4 GHz and shortly after the discovery of its repeating nature, detected up to 2.4 GHz. Here we present three bursts detected with the 100-m Effelsberg radio telescope at 4.85 GHz. All three bursts exhibited frequency structure on broad and narrow frequency scales. Using an autocorrelation function analysis, we measured a characteristic bandwidth of the small-scale structure of 6.4$\\pm$1.6 MHz, which is consistent with the diffractive scintillation bandwidth for this line of sight through the Galactic interstellar medium (ISM) predicted by the NE2001 model. These were the only detections in a campaign totaling 22 hours in 10 observing epochs spanning five months. The observed burst detection rate within this observation was inconsistent with a Poisson process with a constant average occurrence rate; three bursts arrived in the final 0.3 hr of a 2 hr observation on 2016 August 20. We therefore observed a change in the rate of detectable bursts during this observation, and we argue that boosting by diffractive interstellar scintillations may have played a role in the detectability. Understanding whether changes in the detection rate of bursts from FRB 121102 observed at other radio frequencies and epochs are also a product of propagation effects, such as scintillation boosting by the Galactic ISM or plasma lensing in the host galaxy, or an intrinsic property of the burst emission will require further observations.\n",
      "\n"
     ]
    }
   ],
   "source": [
    "# abstract\n",
    "print_abstract_query_results(abstract_data, rst_abstract)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "id": "350e622d",
   "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.10.12"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 5
}