{
 "cells": [
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "# Closest Point Lookup"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 1,
   "metadata": {},
   "outputs": [],
   "source": [
    "# Import of necessary libraries\n",
    "import pathlib\n",
    "import detl\n",
    "import numpy as np"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Parsing and Accessing the Raw Data\n",
    "\n",
    "Alternatively, instead of using pathlib.Path(), a filepath can be given as well (as a rawstring)."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": [
    "# Parsing of data file\n",
    "path_to_file = pathlib.Path('..', 'tests', 'testfiles', 'v4_NT-WMB-2.Control.csv')\n",
    "ddata = detl.parse(path_to_file) "
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Application Examples"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Example 1: Retrieve data closest to a set of given timepoints"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 4,
   "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>timestamp</th>\n",
       "      <th>duration</th>\n",
       "      <th>process_time</th>\n",
       "      <th>volume_pv</th>\n",
       "      <th>temperature_sp</th>\n",
       "      <th>temperature_pv</th>\n",
       "      <th>temperature_out</th>\n",
       "      <th>stirrer_speed_sp</th>\n",
       "      <th>stirrer_speed_pv</th>\n",
       "      <th>stirrer_torque_pv</th>\n",
       "      <th>...</th>\n",
       "      <th>loop_b_sp</th>\n",
       "      <th>loop_b_pv</th>\n",
       "      <th>loop_b_out</th>\n",
       "      <th>loop_c_sp</th>\n",
       "      <th>loop_c_pv</th>\n",
       "      <th>loop_c_out</th>\n",
       "      <th>loop_d_sp</th>\n",
       "      <th>loop_d_pv</th>\n",
       "      <th>loop_d_out</th>\n",
       "      <th>internal_a_start_oszillation</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>6181</th>\n",
       "      <td>2022-03-31 08:56:04+00:00</td>\n",
       "      <td>17.197553</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>1001.642</td>\n",
       "      <td>30.0</td>\n",
       "      <td>29.916</td>\n",
       "      <td>3.613</td>\n",
       "      <td>400.000</td>\n",
       "      <td>398.806</td>\n",
       "      <td>24.956</td>\n",
       "      <td>...</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>6900</th>\n",
       "      <td>2022-03-31 10:56:05+00:00</td>\n",
       "      <td>19.198009</td>\n",
       "      <td>1.999722</td>\n",
       "      <td>1002.739</td>\n",
       "      <td>30.0</td>\n",
       "      <td>30.017</td>\n",
       "      <td>2.803</td>\n",
       "      <td>428.657</td>\n",
       "      <td>426.134</td>\n",
       "      <td>26.928</td>\n",
       "      <td>...</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>7619</th>\n",
       "      <td>2022-03-31 12:56:10+00:00</td>\n",
       "      <td>21.199410</td>\n",
       "      <td>4.001111</td>\n",
       "      <td>1004.966</td>\n",
       "      <td>30.0</td>\n",
       "      <td>30.004</td>\n",
       "      <td>2.008</td>\n",
       "      <td>576.682</td>\n",
       "      <td>577.423</td>\n",
       "      <td>28.522</td>\n",
       "      <td>...</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>8337</th>\n",
       "      <td>2022-03-31 14:56:04+00:00</td>\n",
       "      <td>23.197583</td>\n",
       "      <td>5.999444</td>\n",
       "      <td>1008.937</td>\n",
       "      <td>30.0</td>\n",
       "      <td>29.981</td>\n",
       "      <td>1.866</td>\n",
       "      <td>616.291</td>\n",
       "      <td>618.838</td>\n",
       "      <td>28.102</td>\n",
       "      <td>...</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>8696</th>\n",
       "      <td>2022-03-31 15:56:02+00:00</td>\n",
       "      <td>24.197055</td>\n",
       "      <td>6.998889</td>\n",
       "      <td>1010.868</td>\n",
       "      <td>30.0</td>\n",
       "      <td>30.019</td>\n",
       "      <td>1.910</td>\n",
       "      <td>612.672</td>\n",
       "      <td>613.727</td>\n",
       "      <td>28.898</td>\n",
       "      <td>...</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>9056</th>\n",
       "      <td>2022-03-31 16:56:10+00:00</td>\n",
       "      <td>25.199314</td>\n",
       "      <td>8.001111</td>\n",
       "      <td>1013.138</td>\n",
       "      <td>30.0</td>\n",
       "      <td>29.965</td>\n",
       "      <td>1.976</td>\n",
       "      <td>596.205</td>\n",
       "      <td>607.313</td>\n",
       "      <td>27.884</td>\n",
       "      <td>...</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>9774</th>\n",
       "      <td>2022-03-31 18:56:02+00:00</td>\n",
       "      <td>27.197161</td>\n",
       "      <td>9.998889</td>\n",
       "      <td>1018.850</td>\n",
       "      <td>30.0</td>\n",
       "      <td>29.996</td>\n",
       "      <td>1.936</td>\n",
       "      <td>610.679</td>\n",
       "      <td>611.356</td>\n",
       "      <td>28.209</td>\n",
       "      <td>...</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "<p>7 rows × 59 columns</p>\n",
       "</div>"
      ],
      "text/plain": [
       "                     timestamp   duration  process_time  volume_pv  \\\n",
       "6181 2022-03-31 08:56:04+00:00  17.197553      0.000000   1001.642   \n",
       "6900 2022-03-31 10:56:05+00:00  19.198009      1.999722   1002.739   \n",
       "7619 2022-03-31 12:56:10+00:00  21.199410      4.001111   1004.966   \n",
       "8337 2022-03-31 14:56:04+00:00  23.197583      5.999444   1008.937   \n",
       "8696 2022-03-31 15:56:02+00:00  24.197055      6.998889   1010.868   \n",
       "9056 2022-03-31 16:56:10+00:00  25.199314      8.001111   1013.138   \n",
       "9774 2022-03-31 18:56:02+00:00  27.197161      9.998889   1018.850   \n",
       "\n",
       "      temperature_sp  temperature_pv  temperature_out  stirrer_speed_sp  \\\n",
       "6181            30.0          29.916            3.613           400.000   \n",
       "6900            30.0          30.017            2.803           428.657   \n",
       "7619            30.0          30.004            2.008           576.682   \n",
       "8337            30.0          29.981            1.866           616.291   \n",
       "8696            30.0          30.019            1.910           612.672   \n",
       "9056            30.0          29.965            1.976           596.205   \n",
       "9774            30.0          29.996            1.936           610.679   \n",
       "\n",
       "      stirrer_speed_pv  stirrer_torque_pv  ...  loop_b_sp  loop_b_pv  \\\n",
       "6181           398.806             24.956  ...        0.0        0.0   \n",
       "6900           426.134             26.928  ...        0.0        0.0   \n",
       "7619           577.423             28.522  ...        0.0        0.0   \n",
       "8337           618.838             28.102  ...        0.0        0.0   \n",
       "8696           613.727             28.898  ...        0.0        0.0   \n",
       "9056           607.313             27.884  ...        0.0        0.0   \n",
       "9774           611.356             28.209  ...        0.0        0.0   \n",
       "\n",
       "      loop_b_out  loop_c_sp  loop_c_pv  loop_c_out  loop_d_sp  loop_d_pv  \\\n",
       "6181         0.0        0.0        0.0         0.0        0.0        0.0   \n",
       "6900         0.0        0.0        0.0         0.0        0.0        0.0   \n",
       "7619         0.0        0.0        0.0         0.0        0.0        0.0   \n",
       "8337         0.0        0.0        0.0         0.0        0.0        0.0   \n",
       "8696         0.0        0.0        0.0         0.0        0.0        0.0   \n",
       "9056         0.0        0.0        0.0         0.0        0.0        0.0   \n",
       "9774         0.0        0.0        0.0         0.0        0.0        0.0   \n",
       "\n",
       "      loop_d_out  internal_a_start_oszillation  \n",
       "6181         0.0                           0.0  \n",
       "6900         0.0                           0.0  \n",
       "7619         0.0                           0.0  \n",
       "8337         0.0                           0.0  \n",
       "8696         0.0                           0.0  \n",
       "9056         0.0                           0.0  \n",
       "9774         0.0                           0.0  \n",
       "\n",
       "[7 rows x 59 columns]"
      ]
     },
     "execution_count": 4,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "# Define some timepoints to look up (parameter: 'process_time')\n",
    "sampling_times = np.array([0, 2, 4, 6, 7, 8, 10])\n",
    "\n",
    "# Retrieve data from reactor 2 at the defined points\n",
    "ddata[2].get_closest_data(sampling_times)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Example 2: Retrieve data closest to a set of defined total pump volumes"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 6,
   "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>timestamp</th>\n",
       "      <th>duration</th>\n",
       "      <th>process_time</th>\n",
       "      <th>volume_pv</th>\n",
       "      <th>temperature_sp</th>\n",
       "      <th>temperature_pv</th>\n",
       "      <th>temperature_out</th>\n",
       "      <th>stirrer_speed_sp</th>\n",
       "      <th>stirrer_speed_pv</th>\n",
       "      <th>stirrer_torque_pv</th>\n",
       "      <th>...</th>\n",
       "      <th>loop_b_sp</th>\n",
       "      <th>loop_b_pv</th>\n",
       "      <th>loop_b_out</th>\n",
       "      <th>loop_c_sp</th>\n",
       "      <th>loop_c_pv</th>\n",
       "      <th>loop_c_out</th>\n",
       "      <th>loop_d_sp</th>\n",
       "      <th>loop_d_pv</th>\n",
       "      <th>loop_d_out</th>\n",
       "      <th>internal_a_start_oszillation</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>6723</th>\n",
       "      <td>2022-03-31 10:26:22+00:00</td>\n",
       "      <td>18.702739</td>\n",
       "      <td>1.462500</td>\n",
       "      <td>1002.014</td>\n",
       "      <td>30.0</td>\n",
       "      <td>30.023</td>\n",
       "      <td>3.157</td>\n",
       "      <td>401.879</td>\n",
       "      <td>402.688</td>\n",
       "      <td>27.091</td>\n",
       "      <td>...</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>7708</th>\n",
       "      <td>2022-03-31 13:10:52+00:00</td>\n",
       "      <td>21.444261</td>\n",
       "      <td>4.203889</td>\n",
       "      <td>1004.968</td>\n",
       "      <td>30.0</td>\n",
       "      <td>29.989</td>\n",
       "      <td>2.909</td>\n",
       "      <td>592.199</td>\n",
       "      <td>593.634</td>\n",
       "      <td>34.763</td>\n",
       "      <td>...</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>8590</th>\n",
       "      <td>2022-03-31 15:38:10+00:00</td>\n",
       "      <td>23.899420</td>\n",
       "      <td>6.659167</td>\n",
       "      <td>1010.046</td>\n",
       "      <td>30.0</td>\n",
       "      <td>29.999</td>\n",
       "      <td>2.086</td>\n",
       "      <td>630.550</td>\n",
       "      <td>630.419</td>\n",
       "      <td>35.504</td>\n",
       "      <td>...</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "<p>3 rows × 59 columns</p>\n",
       "</div>"
      ],
      "text/plain": [
       "                     timestamp   duration  process_time  volume_pv  \\\n",
       "6723 2022-03-31 10:26:22+00:00  18.702739      1.462500   1002.014   \n",
       "7708 2022-03-31 13:10:52+00:00  21.444261      4.203889   1004.968   \n",
       "8590 2022-03-31 15:38:10+00:00  23.899420      6.659167   1010.046   \n",
       "\n",
       "      temperature_sp  temperature_pv  temperature_out  stirrer_speed_sp  \\\n",
       "6723            30.0          30.023            3.157           401.879   \n",
       "7708            30.0          29.989            2.909           592.199   \n",
       "8590            30.0          29.999            2.086           630.550   \n",
       "\n",
       "      stirrer_speed_pv  stirrer_torque_pv  ...  loop_b_sp  loop_b_pv  \\\n",
       "6723           402.688             27.091  ...        0.0        0.0   \n",
       "7708           593.634             34.763  ...        0.0        0.0   \n",
       "8590           630.419             35.504  ...        0.0        0.0   \n",
       "\n",
       "      loop_b_out  loop_c_sp  loop_c_pv  loop_c_out  loop_d_sp  loop_d_pv  \\\n",
       "6723         0.0        0.0        0.0         0.0        0.0        0.0   \n",
       "7708         0.0        0.0        0.0         0.0        0.0        0.0   \n",
       "8590         0.0        0.0        0.0         0.0        0.0        0.0   \n",
       "\n",
       "      loop_d_out  internal_a_start_oszillation  \n",
       "6723         0.0                           0.0  \n",
       "7708         0.0                           0.0  \n",
       "8590         0.0                           0.0  \n",
       "\n",
       "[3 rows x 59 columns]"
      ]
     },
     "execution_count": 6,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "# Define pump volumes for data retrieval\n",
    "pump_volumes = np.array([2, 5, 10])\n",
    "\n",
    "# Retrieve data from reactor 4 at the defined points\n",
    "ddata[4].get_closest_data(pump_volumes, reference='pump_b_volume_pv')"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": []
  }
 ],
 "metadata": {
  "kernelspec": {
   "display_name": "Python 3 (ipykernel)",
   "language": "python",
   "name": "python3"
  },
  "language_info": {
   "codemirror_mode": {
    "name": "ipython",
    "version": 3
   },
   "file_extension": ".py",
   "mimetype": "text/x-python",
   "name": "python",
   "nbconvert_exporter": "python",
   "pygments_lexer": "ipython3",
   "version": "3.9.12"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 2
}