368 lines
69 KiB
Plaintext
368 lines
69 KiB
Plaintext
|
{
|
||
|
|
"cells": [
|
||
|
|
{
|
||
|
|
"cell_type": "markdown",
|
||
|
|
"metadata": {},
|
||
|
|
"source": [
|
||
|
|
"Scripts for Calibration of Cernox Sensors\n",
|
||
|
|
"---------------------------------------\n",
|
||
|
|
"\n",
|
||
|
|
"Make a copy of this notebook for an other run."
|
||
|
|
]
|
||
|
|
},
|
||
|
|
{
|
||
|
|
"cell_type": "code",
|
||
|
|
"execution_count": 1,
|
||
|
|
"metadata": {},
|
||
|
|
"outputs": [],
|
||
|
|
"source": [
|
||
|
|
"import numpy as np\n",
|
||
|
|
"import matplotlib.pyplot as plt\n",
|
||
|
|
"import math\n",
|
||
|
|
"from scipy.interpolate import splrep, splev\n",
|
||
|
|
"from zcalib import read_curve, convert_res, compare_calib, make_calib, logrange, Sensor, CalibRun, nplog, npexp"
|
||
|
|
]
|
||
|
|
},
|
||
|
|
{
|
||
|
|
"cell_type": "code",
|
||
|
|
"execution_count": 5,
|
||
|
|
"metadata": {},
|
||
|
|
"outputs": [
|
||
|
|
{
|
||
|
|
"name": "stdout",
|
||
|
|
"output_type": "stream",
|
||
|
|
"text": [
|
||
|
|
"('lsdat', 'READ /afs/psi.ch/project/SampleEnvironment/SE_internal/Thermometer_calibs/2012/73027 Cernox 5/X75610.dat')\n",
|
||
|
|
"('zdat', 'READ /home/l_samenv/sea/calib_scripts/calib_data/calib2024-11-05_p0_c1.dat')\n",
|
||
|
|
"('zdat', 'READ /home/l_samenv/sea/calib_scripts/calib_data/calib2024-11-05_p1_c1.dat')\n",
|
||
|
|
"('zdat', 'READ /home/l_samenv/sea/calib_scripts/calib_data/calib2024-11-05_p2_c1.dat')\n",
|
||
|
|
"{'selected': 0, 'averaged': 62}\n",
|
||
|
|
"('zdat', 'READ /home/l_samenv/sea/calib_scripts/calib_data/calib2024-11-05_p0_c2.dat')\n",
|
||
|
|
"('zdat', 'READ /home/l_samenv/sea/calib_scripts/calib_data/calib2024-11-05_p1_c2.dat')\n",
|
||
|
|
"('zdat', 'READ /home/l_samenv/sea/calib_scripts/calib_data/calib2024-11-05_p2_c2.dat')\n",
|
||
|
|
"{'selected': 1, 'averaged': 61}\n",
|
||
|
|
"('zdat', 'READ /home/l_samenv/sea/calib_scripts/calib_data/calib2024-11-05_p0_c4.dat')\n",
|
||
|
|
"('zdat', 'READ /home/l_samenv/sea/calib_scripts/calib_data/calib2024-11-05_p1_c4.dat')\n",
|
||
|
|
"('zdat', 'READ /home/l_samenv/sea/calib_scripts/calib_data/calib2024-11-05_p2_c4.dat')\n",
|
||
|
|
"{'selected': 1, 'averaged': 61}\n"
|
||
|
|
]
|
||
|
|
}
|
||
|
|
],
|
||
|
|
"source": [
|
||
|
|
"run = CalibRun([\n",
|
||
|
|
" Sensor(3, 'X75610'), # the reference sensor must be the first\n",
|
||
|
|
" Sensor(1, 'X219219', 'CX-1050-SD-HT'),\n",
|
||
|
|
" Sensor(2, 'X219231', 'CX-1050-SD-HT'),\n",
|
||
|
|
" Sensor(4, 'X219225', 'CX-1050-SD-HT'),\n",
|
||
|
|
" ],\n",
|
||
|
|
" t_points = (1.2, 1.4) + logrange(1.6, 310, n=195) + (330,), # the points to be used in the cal file\n",
|
||
|
|
" caldate = '2024-11-05', # the first measuring day!\n",
|
||
|
|
" logT = False,\n",
|
||
|
|
" logR = True,\n",
|
||
|
|
" calib_data_file = '/home/l_samenv/sea/calib_scripts/calib_data/calib%s_p%d_c%d.dat',\n",
|
||
|
|
" outputpath='%s/%s.340')\n",
|
||
|
|
"# smooth depends on number of measured points (1e-7 for 60, 0.8e-7 for 48 and 0.4e-7 for 24 points)\n",
|
||
|
|
"run.make(diflim=0.001, smoothref=1e-7, smoothtst=0.4e-7)"
|
||
|
|
]
|
||
|
|
},
|
||
|
|
{
|
||
|
|
"cell_type": "code",
|
||
|
|
"execution_count": 3,
|
||
|
|
"metadata": {},
|
||
|
|
"outputs": [
|
||
|
|
{
|
||
|
|
"name": "stdout",
|
||
|
|
"output_type": "stream",
|
||
|
|
"text": [
|
||
|
|
"('z340', 'READ 2018-11-22/X133982.340')\n",
|
||
|
|
"('z340', 'READ 2018-11-20/X133982.340')\n",
|
||
|
|
"('z340', 'READ 2018-11-22/X133928.340')\n",
|
||
|
|
"('z340', 'READ 2018-11-20/X133928.340')\n",
|
||
|
|
"('z340', 'READ 2018-11-22/X131824.340')\n",
|
||
|
|
"('z340', 'READ 2018-11-20/X131824.340')\n",
|
||
|
|
"('z340', 'READ 2018-11-22/X132254.340')\n",
|
||
|
|
"('z340', 'READ 2018-11-20/X132254.340')\n",
|
||
|
|
"('z340', 'READ 2018-11-22/X137461.340')\n",
|
||
|
|
"('z340', 'READ 2018-11-20/X137461.340')\n"
|
||
|
|
]
|
||
|
|
},
|
||
|
|
{
|
||
|
|
"ename": "IOError",
|
||
|
|
"evalue": "[Errno 2] No such file or directory: '2018-11-20/X137461.340'",
|
||
|
|
"output_type": "error",
|
||
|
|
"traceback": [
|
||
|
|
"\u001b[0;31m\u001b[0m",
|
||
|
|
"\u001b[0;31mIOError\u001b[0mTraceback (most recent call last)",
|
||
|
|
"\u001b[0;32m<ipython-input-3-e7e8a820ce5b>\u001b[0m in \u001b[0;36m<module>\u001b[0;34m()\u001b[0m\n\u001b[1;32m 5\u001b[0m \u001b[0;32mfor\u001b[0m \u001b[0msensno\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0mkind\u001b[0m \u001b[0;32min\u001b[0m \u001b[0mtestlist\u001b[0m\u001b[0;34m:\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m 6\u001b[0m \u001b[0mr0\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0mt0\u001b[0m \u001b[0;34m=\u001b[0m \u001b[0mread_curve\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0;34m'2018-11-22/%s.340'\u001b[0m \u001b[0;34m%\u001b[0m \u001b[0msensno\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0;34m'z340'\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0;32m----> 7\u001b[0;31m \u001b[0mr1\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0mt1\u001b[0m \u001b[0;34m=\u001b[0m \u001b[0mread_curve\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0;34m'2018-11-20/%s.340'\u001b[0m \u001b[0;34m%\u001b[0m \u001b[0msensno\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0mkind\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0m\u001b[1;32m 8\u001b[0m \u001b[0mdiff\u001b[0m \u001b[0;34m=\u001b[0m \u001b[0mcompare_calib\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0mr1\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0mt1\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0mr0\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0mt0\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m 9\u001b[0m \u001b[0mplt\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mplot\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0mt0\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0mdiff\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0;34m'-'\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n",
|
||
|
|
"\u001b[0;32m/home/l_samenv/sea/calib_scripts/zcalib.py\u001b[0m in \u001b[0;36mread_curve\u001b[0;34m(filename, kind, instance, **filterargs)\u001b[0m\n\u001b[1;32m 118\u001b[0m \u001b[0;32mprint\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0mkind\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0;34m\"READ %s\"\u001b[0m \u001b[0;34m%\u001b[0m \u001b[0mfilename\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m 119\u001b[0m \u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0;32m--> 120\u001b[0;31m \u001b[0;32mwith\u001b[0m \u001b[0mopen\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0mfilename\u001b[0m\u001b[0;34m)\u001b[0m \u001b[0;32mas\u001b[0m \u001b[0mf\u001b[0m\u001b[0;34m:\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0m\u001b[1;32m 121\u001b[0m \u001b[0mcurves\u001b[0m \u001b[0;34m=\u001b[0m \u001b[0;34m[\u001b[0m\u001b[0;34m[\u001b[0m\u001b[0;34m]\u001b[0m \u001b[0;32mfor\u001b[0m \u001b[0mc\u001b[0m \u001b[0;32min\u001b[0m \u001b[0mfilter\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0moutput\u001b[0m\u001b[0;34m]\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m 122\u001b[0m \u001b[0;32mfor\u001b[0m \u001b[0mline\u001b[0m \u001b[0;32min\u001b[0m \u001b[0mf\u001b[0m\u001b[0;34m:\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n",
|
||
|
|
"\u001b[0;31mIOError\u001b[0m: [Errno 2] No such file or directory: '2018-11-20/X137461.340'"
|
||
|
|
]
|
||
|
|
},
|
||
|
|
{
|
||
|
|
"data": {
|
||
|
|
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAmsAAAFpCAYAAADDbG9yAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDIuMi40LCBodHRwOi8vbWF0cGxvdGxpYi5vcmcv7US4rQAAIABJREFUeJzt3Xm0HVd9p/3nd4Y7ap4n2xIewDYG2wjbYEyTmMGQweRdDIYMfhNoSAe608mbpslMSGd10m8ndOgmZNGEN4YMhjhh4RBoxwyOgeBB2AZjGVu2PEiy5qvpjmeo/f5x6l5dyVeS60qqe42ej9dZVbWrTp2tuudaX+29a1eklJAkSdLsVJnpCkiSJOnYDGuSJEmzmGFNkiRpFjOsSZIkzWKGNUmSpFnMsCZJkjSLGdYkSZJmMcOaJEnSLGZYkyRJmsUMa5IkSbNYbaYrcCotWbIkrV27dqarIUmSdELf+c539qSUlp7ouB+qsLZ27Vo2bNgw09WQJEk6oYh46rkcZzeoJEnSLGZYkyRJmsUMa5IkSbOYYU2SJGkWM6xJkiTNYoY1SZKkWcywJkmSNIsZ1iRJkmYxw5okSdIsZliTJEmaxQxrkiRJs5hhrYh9T8Gm26HdnOmaSJKkM4RhrYgf/BP89VugMTTTNZEkSWcIw9q0pJmugCRJOkMY1oqI6CyTYU2SJJXDsFZIzHQFJEnSGcawJkmSNIsZ1oqwG1SSJJXMsFaI3aCSJKlchrVpsWVNkiSVw7BWhN2gkiSpZIY1SZKkWcywNi22rEmSpHIY1oqwG1SSJJXMsFaId4NKkqRyGdamxZY1SZJUDsNaEXaDSpKkkhnWCrEbVJIklcuwNi22rEmSpHIY1oqwG1SSJJXMsFbIeDeoYU2SJJXDsFbA/VsOADDUaM9wTSRJ0pnCsFbA4FgLgHY7m+GaSJKkM4VhrZBON2hKhjVJklQOw1oRztwhSZJKZlibDu8GlSRJJTGsFRAx3g1qWJMkSeUwrBXSuVyGNUmSVBbDWhHjc+I6z5okSSqJYa2AcFJcSZJUMsNaAcmnTUmSpJIZ1goxrUmSpHIZ1goIH+QuSZJKZlgrxFlxJUlSuQxr0+DdoJIkqSyGtSLGJ8XNDGuSJKkchrUCJp5gYMuaJEkqiWFtWgxrkiSpHIa1QrwbVJIklcuwVsREN6gkSVI5DGsFTDxsKstmtB6SJOnMcUrCWkRcFxGPRMRjEfHBKfZ3R8Rn8/13R8TavHxxRHw9IgYj4n8d9Z6XRcSD+Xs+GhMz0s6cNN6yZtOaJEkqyUmHtYioAh8D3ghcBLwjIi466rB3AftSSucBHwH+KC8fBX4b+LUpTv1x4N8C5+ev6062ricrnBRXkiSV7FS0rF0BPJZS2pxSagA3A9cfdcz1wE35+i3AtRERKaWhlNI36YS2CRGxEpiXUrorpZSATwNvPgV1PSXCUWuSJKkkpyKsrQa2TNrempdNeUxKqQUcABaf4JxbT3BOACLiPRGxISI27N69u2DVC3JSXEmSVLLn/Q0GKaVPpJTWp5TWL1269DR/mpPiSpKkcp2KsLYNOGvS9pq8bMpjIqIGzAf2nuCca05wztIdvsXBsCZJkspxKsLavcD5EbEuIrqAG4BbjzrmVuDGfP0twNfysWhTSiltBw5GxFX5XaA/B3zhFNT1JI3fDWpYkyRJ5aid7AlSSq2IeD9wG1AFPpVSeigiPgxsSCndCvwF8JmIeAwYoBPoAIiIJ4F5QFdEvBl4fUppI/BLwF8CvcCX89fMmvnZQyRJ0hnmpMMaQErpS8CXjir7nUnro8Bbj/Hetcco3wC8+FTU71QZn7rDSXElSVJZnvc3GJTKhjVJklQyw1oRPhtUkiSVzLA2Hd5gIEmSSmJYKyLyy5UcsyZJksphWCtg4gaDGa6HJEk6cxjWihi/wcBuUEmSVBLD2jT4uClJklQWw1oBMX65zGqSJKkkhrUi8m7Q5A0GkiSpJIa1AsJZcSVJUskMawWk8EHukiSpXIa1AmL8Qe5mNUmSVBLD2jQkHLMmSZLKYVgrwm5QSZJUMsNaAd5gIEmSymZYm4awZU2SJJXEsFZE/iB3u0ElSVJZDGsFeDOoJEkqm2FtWoxrkiSpHIa1Isab1jLDmiRJKodhrYDxSXGNapIkqSyGtULGB605Ka4kSSqHYa2I8W5Q29YkSVJJDGuF2A0qSZLKZVgrYKJdzXnWJElSSQxrReTdoD7BQJIklcWwVoB3g0qSpLIZ1qbDljVJklQSw1oR3g0qSZJKZlgrIBh/kPsMV0SSJJ0xDGtFjM+J66S4kiSpJIa1AuKopSRJ0ulmWCvCu0ElSVLJDGvT4KS4kiSpLIa1IiK/XI5ZkyRJJTGsFRDhaDVJklQuw9p02A0qSZJKYlgrxElxJUlSuQxrBUw8G9SsJkmSSmJYK2B8yFqyZU2SJJXEsFZEntbCpjVJklQSw1ohToorSZLKZVgrYGLmDlvWJElSSQxrhXQul08wkCRJZTGsFeGkuJIkqWSGtWmxZU2SJJXDsFaAY9YkSVLZDGsFRP4gd+dZkyRJZTGsTYcta5IkqSSGtQIifDaoJEkql2GtkInnTUmSJJXCsFbExJg1SZKkchjWCjjcC5rNaD0kSdKZ45SEtYi4LiIeiYjHIuKDU+zvjojP5vvvjoi1k/b9el7+SES8YVL5kxHxYEQ8EBEbTkU9T954N6hta5IkqRy1kz1BRFSBjwGvA7YC90bErSmljZMOexewL6V0XkTcAPwR8PaIuAi4AbgYWAV8JSIuSCm18/f9SEppz8nW8VQzqkmSpLKcipa1K4DHUkqbU0oN4Gbg+qOOuR64KV+/Bbg2OrdWXg/cnFIaSyk9ATyWn29WclJcSZJUtlMR1lYDWyZtb83LpjwmpdQCDgCLT/DeBPxzRHwnIt5zCup58mL8chnWJElSOU66G/Q0elVKaVtELANuj4gfpJTuPPqgPMi9B+Dss88+rRUKH+QuSZJKdipa1rYBZ03aXpOXTXlMRNSA+cDe4703pTS+3AV8nmN0j6aUPpFSWp9SWr906dKT/sM8J3aDSpKkkpyKsHYvcH5ErIuILjo3DNx61DG3Ajfm628BvpZSSnn5DfndouuA84F7IqI/IuYCREQ/8Hrg+6egridlvGXNZ4NKkqSynHQ3aEqpFRHvB24DqsCnUkoPRcSHgQ0ppVuBvwA+ExGPAQN0Ah35cZ8DNgIt4H0ppXZELAc+n4ejGvA3KaX/c7J1PWnjY9bMapIkqSSnZMxaSulLwJeOKvudSeujwFuP8d4/AP7gqLLNwEtPRd1OpcNPBjWtSZKkcvgEgwLGu0HDMWuSJKkkhrVpMKpJkqSyGNYKcFJcSZJUNsNaEfkNBsmwJkmSSmJYK8JJcSVJUskMa9MQjlqTJEklMawVEPnkHXaDSpKkshjWCohKNV/LZrQekiTpzGFYK8J51iRJUskMa0VMPG7KljVJklQOw1oBYViTJEklM6wVMRHW7AaVJEnlMKwVEJXxSXFtWZMkSeUwrBVQrdgNKkmSymVYK2B86g5b1iRJUlkMawVUbFmTJEklM6wVMDEprjcYSJKkkhjWCpi4wSBrz3BNJEnSmcKwVkDFedYkSVLJDGsFVKp2g0qSpHIZ1go4fIOB3aCSJKkchrUCKhM3GNgNKkmSymFYK6AycYOB3aCSJKkchrUCwnnWJElSyQxrBfi4KUmSVDbDWgGVSoUsBcm7QSVJUkkMawUMfeUrbPmXRdBozHRVJEnSGcKwVkBrxw6Gd3RDy6k7JElSOQxrRdQ6U3dE27AmSZLKYVgrIKq1zophTZIklcSwVkDUDGuSJKlchrUCwm5QSZJUMsNaEXnLWmSGNUmSVA7DWgHjY9ZS20lxJUlSOQxrBUx0g2aGNUmSVA7DWhETNxgY1iRJUjkMawU4dYckSSqbYa2AqI/fYGDLmiRJKodhrYCodsasYViTJEklMawVkXeDhmPWJElSSQxrBYx3g9qyJkmSymJYK2CiG9QbDCRJUkkMa0XUvMFAkiSVy7BWwMSD3LM0sxWRJEl
|
||
|
|
"text/plain": [
|
||
|
|
"<Figure size 720x432 with 1 Axes>"
|
||
|
|
]
|
||
|
|
},
|
||
|
|
"metadata": {
|
||
|
|
"needs_background": "light"
|
||
|
|
},
|
||
|
|
"output_type": "display_data"
|
||
|
|
}
|
||
|
|
],
|
||
|
|
"source": [
|
||
|
|
"# compare known sensors (skip if not applicable)\n",
|
||
|
|
"plt.figure(figsize=(10, 6))\n",
|
||
|
|
"tmin,tmax,dmax=1.4,310,0.001\n",
|
||
|
|
"testlist = (('X133982', 'z340'),('X133928', 'z340'),('X131824', 'z340'),('X132254', 'z340'),('X137461', 'z340'))\n",
|
||
|
|
"for sensno, kind in testlist:\n",
|
||
|
|
" r0, t0 = read_curve('2023-11-05/%s.340' % sensno, 'z340')\n",
|
||
|
|
" r1, t1 = read_curve('2023-11-05/%s.340' % sensno, kind)\n",
|
||
|
|
" diff = compare_calib(r1, t1, r0, t0)\n",
|
||
|
|
" plt.plot(t0, diff, '-')\n",
|
||
|
|
"plt.plot([tmin,tmax,tmax,tmin,tmin], [-dmax,-dmax,dmax,dmax,-dmax], '-')\n",
|
||
|
|
"plt.legend([sensno + \".\" + kind[-3:] for sensno, kind in testlist] + [\"window\"])\n",
|
||
|
|
"plt.xscale('log')\n",
|
||
|
|
"plt.yscale('symlog',linthreshy=dmax)\n",
|
||
|
|
"plt.grid(True, axis='y')\n",
|
||
|
|
"plt.axis([1.0,350,-1,1])\n",
|
||
|
|
"plt.show()"
|
||
|
|
]
|
||
|
|
},
|
||
|
|
{
|
||
|
|
"cell_type": "code",
|
||
|
|
"execution_count": 6,
|
||
|
|
"metadata": {},
|
||
|
|
"outputs": [
|
||
|
|
{
|
||
|
|
"name": "stdout",
|
||
|
|
"output_type": "stream",
|
||
|
|
"text": [
|
||
|
|
"('z340', 'READ 2024-11-05/X219219.340')\n",
|
||
|
|
"('zdat', 'READ /home/l_samenv/sea/calib_scripts/calib_data/calib2024-11-05_p0_c1.dat')\n",
|
||
|
|
"('zdat', 'READ /home/l_samenv/sea/calib_scripts/calib_data/calib2024-11-05_p1_c1.dat')\n",
|
||
|
|
"('zdat', 'READ /home/l_samenv/sea/calib_scripts/calib_data/calib2024-11-05_p2_c1.dat')\n"
|
||
|
|
]
|
||
|
|
},
|
||
|
|
{
|
||
|
|
"data": {
|
||
|
|
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAYgAAAECCAYAAADkaECYAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDIuMi40LCBodHRwOi8vbWF0cGxvdGxpYi5vcmcv7US4rQAAIABJREFUeJzt3Xl8VeW97/HPb0+Zk01CQiAJg8gooGgEO2jRqgerFbW2iqd6bovl2HuxrT3F2um2t7aH9pzX7W2p9vRQpVStWrQOlKJWqxHr0IITMxKZEhBIIPO0p+f+kaghrkACOyaB7/v1yovsZz3PWr9sl/u711rPXtucc4iIiHTl6+8CRERkYFJAiIiIJwWEiIh4UkCIiIgnBYSIiHhSQIiIiCcFhIiIeFJAiIiIp0B/F/AuM8sAfgVEgDLn3O/7uSQRkZNanx5BmNlSMztgZhu6tM82s61mVm5mt3U0XwU87Jz7EnB5X9YlIiJH19enmJYBszs3mJkfuBO4BJgMzDWzyUAxUNHRLd7HdYmIyFH0aUA451YDh7o0zwDKnXPbnXMR4EFgDlBJe0j0eV0iInJ0/XENooj3jxSgPRhmAouBO8zsUuBP3Q02s/nAQiAcCoWGjh49ug9LFRE58bz11lvVzrn8o/UbMBepnXNNwBd60G8JsASgtLTUrV27tq9LExE5oZjZrp70649TOXuAkk6PizvaRERkAOmPgFgDjDOzMWYWAq4FVvRDHSIicgR9Pc31AeBlYIKZVZrZPOdcDFgAPAVsBpY75zb2ZR0iItJ7fXoNwjk3t5v2VcCqvty2iIgcH00nFRERTwoIERHxpIAQERFPCggREfGkgBAREU8KCBER8aSAEBERTwoIERHxpIAQERFPCggREfGkgBAREU8KCBER8aSAEBERTwoIERHxpIAQERFPCggREfE0oALCzE4xs7vN7OH+rkVE5GSXtIAws6VmdsDMNnRpn21mW82s3MxuO9I6nHPbnXPzklWTiIgcu2R+5egy4A7gnncbzMwP3AlcBFQCa8xsBeAHFnUZ/0Xn3IEk1iMiIschaQHhnFttZqO7NM8Ayp1z2wHM7EFgjnNuEXBZsrYtIiLJl8wjCC9FQEWnx5XAzO46m1ke8GNgupl9qyNIuvaZDywEwuFwmLKysuRWLCIiQN8HRK845w4CNx2lzxJgCUBpaambNWvWh1CZiMjJp69nMe0BSjo9Lu5oExGRAa6vA2INMM7MxphZCLgWWNHH2xQRkSRI5jTXB4CXgQlmVmlm85xzMWAB8BSwGVjunNuYrG2KiEjfSeYsprndtK8CViVrOyIi8uEYUJ+kFhGRgUMBISIinhQQIiLiSQEhIiKeFBAiIuJJASEiIp4UECIi4kkBISIinhQQIiLiSQEhIiKeFBAiIuJJASEiIp4UECIi4kkBISIinhQQIiLiSQEhIiKeFBAiIuIpad8olwxmdgVwKZAN3O2c+0s/lyQictJK5ndSLzWzA2a2oUv7bDPbamblZnbbkdbhnHvMOfcl4CbgmmTVJiIivZfMI4hlwB3APe82mJkfuBO4CKgE1pjZCsAPLOoy/ovOuQMdv3+3Y5yIiPSTpAWEc261mY3u0jwDKHfObQcwsweBOc65RcBlXddhZgb8BHjCOfea13bMbD6wEAiHw2HKysqS9SeIiEgnfX0Nogio6PS4Eph5hP43AxcCOWZ2qnPu1107OOeWAEsASktL3axZs5JXrYiIvGdAXaR2zi0GFvd3HSIi0vfTXPcAJZ0eF3e0iYjIANfXAbEGGGdmY8wsBFwLrOjjbYqISBIkc5rrA8DLwAQzqzSzec65GLAAeArYDCx3zm1M1jZFRKTvJHMW09xu2lcBq5K1HRER+XDoVhsiIuJJASEiIp4UECIi4kkBISIinhQQIiLiSQEhIiKeFBAiIuJJASEiIp4UECIi4kkBISIinhQQIiInkRVv9/x+qQoIEZGTyO83/77HfRUQIiInCeccFfUVR+/YQQEhInKSqG2rpSHa0OP+CggRkZPE7obdveqvgBAROUnsrh/EAWFmk8zs12b2sJl9ub/rERE5kexu2I3Pev6yn8yvHF1qZgfMbEOX9tlmttXMys3stiOtwzm32Tl3E/A54GPJqk1ERNqPIIZnDO9x/2QeQSwDZnduMDM/cCdwCTAZmGtmk81sqpmt7PJT0DHmcuDP6GtKRUSSqqKhgpKskh73T+Z3Uq82s9FdmmcA5c657QBm9iAwxzm3CLism/WsAFaY2Z+B+4+nppo/LKd+5crjWYWIyKCXfdllDLnmc+xu2M3M/At6PC5pAdGNIqDzpNtKYGZ3nc1sFnAVkEI3RxBmNh9YCITD4TBlZWXdbnzIffcRqKwkVlzc68JFRE4EgcpKamtreWloBnVtdby8KdrzsX1YV68558qAsqP0WQIsASgtLXWzZs3qtu+uu5dCOMyoe+9JXpEiIoPIrutvAMA3bSRUQlNrfo/H9vUspj1A5xNexR1tIiLyIXr3MxANDeEej+nrgFgDjDOzMWYWAq4Fen6nKBERSYrdDbsxjEQ0t8djkjnN9QHgZWCCmVWa2TznXAxYADwFbAaWO+c2JmubIiLSM7vrd5MVGAou2OMxyZzFNLeb9lVoyqqISL/a3bCbVIaRGuyHD8qJiMjAVVFfQSKSxylDM3s8RgEhInKCi7k4NW01NDbmcGqBAkJERDq0xVsBqKlXQIiISCdtsTYAEpE8xuYrIEREpENrxxFEIpKrIwgREXlfW6yNdF8uPkKMHpre43EKCBGRE1xrvI1gIp+RuemkBPw9HqeAEBE5wbXFW4m25vXq9BIoIERETmhxFyeaiFLfkMNYBYSIiLyrLd4+gynalturGUyggBAROaG1xt6dwaRTTCIi0sm7RxCJyFAdQYiIyPta4634CJCfmUVOWs/v5AoKCBGRE1pbrA1ckFN7efQACgg5yf2ubBW3/PYH7Khq7O9SBo+2RtjyZ0jE+7sS6YHWeBuJeLDX1x9AASEnsYONTTyw7Zs84/sj314+iyfWrOnvkga+aCs8cC08eB2svAWc6++K5Aiao81EExESiYACQqQ3fvboN9gTggsT2WzOaOXxNbfg9IJ3ZI/+K+x8gfLwx+G138FzP+7viuQIKhoq2n9xwV5foAYFhJykKqr3URZ9galtxs9ueJ5zXT6vp9WwZou+Ebdbe9+ATY9xa+pFXOYL83jax3Ev/RKaDvZ3ZdKN3Q27AXAJnWIS6bH7n/0J9X7jhnE3Yv4An5/xFZp9Pv78wu39XdqA1fLSr/lxbh5PDN9KMGcd38+to8pF4dXf9ndp0o1th3YCkBZIYVh2Sq/HD+6AiLXqHKgckzdrX6QwmuDic78MwIzJVzA65mcd66ltaO7n6gag5kP84+2VPJiTweyRV3L3xXfjS2njxoKRtL3yG4hHj2v18YTjzYpanti4k7f264jkeKyrrOXbj67nqw++zsNvvgHOz9j8bMys1+sK9EF9H5r66q1U3nkWReMvxSZdDkVngW9wZF5LrIXGSCPxeJRoHDJTMhiSltPfZSVPPAY1O2ip3kVjSyvNlklt6nBqCNMcTRD0+8hI8ZOZEmBIeoj8rBRSgz2/y+Tx2FNdweZQCxdFR+ALtM8LNzM+VXgxv6p+ghXP3MENV97a7fia1hpeO/Aaexv3UpJVwun5pzMkdciHUjtAY6SRN6veZEfdDppbawlFmpg4ZALTRn6C9LS+qaPx78v4eW4GmW4I/37e9wj6g3yjdCGL1tzOS3V1nL95BUz5zDGte9fBJhYsf5Zt0T8QyGo/xZcfmMKi87/JOcVTkvlnEE/EWbt/Lc/ufpb89HzOLTqXCbkTuh+QSMDGR6BxP3syT2NrYCItUUdhTiojwqkUZKXi9x39hdc5x876ndS11eEzH2PDY8kIZvSo5sZD+9i2ajFDdj1JbcZYwrO/zeiJ0z1Kdfz0yS389+rtpAX9ZKT4acqtJOhLITPl2F7qB1RAmNls4BeAH7jLOfeTI/WvCAa4JCtKVuUjTN7+IJNckMkFZzBpwuWMnHj
|
||
|
|
"text/plain": [
|
||
|
|
"<Figure size 432x288 with 1 Axes>"
|
||
|
|
]
|
||
|
|
},
|
||
|
|
"metadata": {
|
||
|
|
"needs_background": "light"
|
||
|
|
},
|
||
|
|
"output_type": "display_data"
|
||
|
|
},
|
||
|
|
{
|
||
|
|
"name": "stdout",
|
||
|
|
"output_type": "stream",
|
||
|
|
"text": [
|
||
|
|
"('z340', 'READ 2024-11-05/X219231.340')\n",
|
||
|
|
"('zdat', 'READ /home/l_samenv/sea/calib_scripts/calib_data/calib2024-11-05_p0_c2.dat')\n",
|
||
|
|
"('zdat', 'READ /home/l_samenv/sea/calib_scripts/calib_data/calib2024-11-05_p1_c2.dat')\n",
|
||
|
|
"('zdat', 'READ /home/l_samenv/sea/calib_scripts/calib_data/calib2024-11-05_p2_c2.dat')\n"
|
||
|
|
]
|
||
|
|
},
|
||
|
|
{
|
||
|
|
"data": {
|
||
|
|
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAYgAAAECCAYAAADkaECYAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDIuMi40LCBodHRwOi8vbWF0cGxvdGxpYi5vcmcv7US4rQAAIABJREFUeJzt3Xl8VNX9//HXZ7bsJOxLwqZsgliFCHUlKijWtWoVal1aKvXbar+/aumXVvut2sV+u1oVpdQFV9DihoqiViPagmCpCsgWWcO+ZSHJzGRmzu+PRIzxJiQwMQHez8cjD7nnnnPumXgz77n3nrnXnHOIiIjU52vtAYiISNukgBAREU8KCBER8aSAEBERTwoIERHxpIAQERFPCggREfGkgBAREU+B1h7Ap8wsA7gPiAKFzrknWnlIIiJHtBY9gjCzh8xsu5ktrVc+1sxWmlmRmU2uLb4EmOWcuw64sCXHJSIi+9fSp5imA2PrFpiZH5gCnAsMBsab2WAgD9hYWy3ewuMSEZH9aNGAcM7NA3bXKx4BFDnn1jjnosBM4CKgmJqQaPFxiYjI/rXGNYhcPjtSgJpgGAncDdxrZucBLzbU2MwmApOAnFAo1KlPnz4tOFQRkcPPqlWrdjrnOu+vXpu5SO2cqwC+3YR604BpAPn5+e79999v6aGJiBxWzGx9U+q1xqmcTUDPOst5tWUiItKGtEZALAL6m1lfMwsB44DZrTAOERFpREtPc50BzAcGmlmxmU1wzsWAG4C5wHLgaefcspYch4iINF+LXoNwzo1voHwOMKclty0iIgdH00lFRMSTAkJERDwpIERExJMCQkREPCkgRETEkwJCREQ8KSBERMSTAkJERDwpIERExJMCQkREPCkgRETEkwJCREQ8KSBERMSTAkJERDwpIERExJMCQkREPLWpgDCzo8zsQTOb1dpjERE50iUtIMzsITPbbmZL65WPNbOVZlZkZpMb68M5t8Y5NyFZYxIRkQOXzEeOTgfuBR79tMDM/MAUYAxQDCwys9mAH7izXvvvOOe2J3E8IiJyEJIWEM65eWbWp17xCKDIObcGwMxmAhc55+4Ezk/WtkVEJPmSeQThJRfYWGe5GBjZUGUz6wj8GjjBzH5aGyT160wEJgE5OTk5FBYWJnfEIiICtHxANItzbhdw/X7qTAOmAeTn57uCgoIvYWQiIkeelp7FtAnoWWc5r7ZMRETauJYOiEVAfzPra2YhYBwwu4W3KSIiSZDMaa4zgPnAQDMrNrMJzrkYcAMwF1gOPO2cW5asbYqISMtJ5iym8Q2UzwHmJGs7IiLy5WhT36QWEZG2QwEhIiKeFBAiIuJJASEiIp4UECIi4kkBISIinhQQIiLiSQEhIiKeFBAiIuJJASEiIp4UECIi4kkBISIinhQQIiLiSQEhIiKeFBAiIuJJASEiIp4UECIi4ilpT5RLBjO7GDgPaAc86Jx7rZWHJCJyxErmM6kfMrPtZra0XvlYM1tpZkVmNrmxPpxzzzvnrgOuB65I1thERKT5knkEMR24F3j00wIz8wNTgDFAMbDIzGYDfuDOeu2/45zbXvvvW2vbiYhIK0laQDjn5plZn3rFI4Ai59waADObCVzknLsTOL9+H2ZmwG+BV5xzi722Y2YTgUlATk5ODoWFhcl6CSIiUkdLX4PIBTbWWS4GRjZS/0ZgNJBtZv2cc1PrV3DOTQOmAeTn57uCgoLkjVZERPZpUxepnXN3A3e39jhERKTlp7luAnrWWc6rLRMRkTaupQNiEdDfzPqaWQgYB8xu4W2KiEgSJHOa6wxgPjDQzIrNbIJzLgbcAMwFlgNPO+eWJWubIiLScpI5i2l8A+VzgDnJ2o6IiBy4Dbsqm1xXt9oQETmCfPOF/9fkugoIEZEjhHOOUreiyfUVECIiR4jVO3dAYE+T6ysgRESOEG+v/ahZ9RUQIiJHiMVbP25WfQWEiMgRoqhkFZZIb3J9BYSIyBFiZ3QdWb5eTa6vgBAROQJURKJU+zeTl3F0k9soIEREjgDvrluJ+aoZ3GlQk9soIEREjgDzi5cAcErPoU1uo4AQETkCLNu5Aud8nNxrSJPbtKnnQYiISMvYVPEJIbqSHkxtchsdQYiIHOacc5S7jXRO6dOsdgoIEZHDXFHtLTaOzh7QrHYKCBGRw1zhuppbbAzvPrhZ7RQQIiKHucVbap7TVtD3+Ga1a1MBYWbHmNlUM5tlZv/V2uMRETkcfFK6GotncFT77s1ql8xHjj5kZtvNbGm98rFmttLMisxscmN9OOeWO+euBy4HTknW2EREjmQ7o+vI8vfCzJrVLplHENOBsXULzMwPTAHOBQYD481ssJkNNbOX6v10qW1zIfAyekypiMhBq4hEifo2k5fe9FtsfCqZz6SeZ2Z96hWPAIqcc2sAzGwmcJFz7k7g/Ab6mQ3MNrOXgScPZkx7nnqaspdeOpguREQOaburKhidF6H9RU2/xcanWvqLcrnAxjrLxcDIhiqbWQFwCZBCA0cQZjYRmATk5OTkUFhY2ODG2z/+OIHiYmJ5ec0euIjIYWHDek4tS7CxwBp9v/TSpr5J7ZwrBAr3U2caMA0gPz/fFRQUNFh3/YMPQU4OvR97NHmDFBE5hLx67iiIb2fi175BRiitWW1behbTJqBnneW82jIREfkSxBNxcL5mhwO0fEAsAvqbWV8zCwHjgNktvE0REakVd3EO9K0+mdNcZwDzgYFmVmxmE5xzMeAGYC6wHHjaObcsWdsUEZHGJVwCO8C3+mTOYhrfQPkcNGVVRKRVJFwcn7XyEYSIiLQ9jgSmgBARkfocCfzmP6C2CggRkcNUJBYHEvgUECIiUldJZTVYAr9OMYmISF17KqNAgoDvwOYjKSBERA5TO8qrwBwBX51TTG/c3uT2CggRkcPU9r0lAJ8FhHPw4Ywmt1dAiIgcpnZUlgMQ9NcGxK4iKN/S5PYKCBGRw9SOilIAgp9eg1j7drPaKyBERA5Tu6tqjiACnx5BrH0H2uU2ub0CQkTkMLWnqgyg5otyiQSsewf6jmpyewWEiMhhqjRSAVDzRbntH0PlLuh7epPbKyBERA5T5ZGaU0x+88PaeTWFfU9rcnsFhIjIYaqiuuYIwu/z1QREh6Mgu+mPYFZAiIgcpqrilQD48cH6fzbr9BIoIEREDkvOOcLxSgwfVl0BkTIFhIiIwKaSKpxFap4mV1XzfQj6NP36AyggRKSZVu5czpSFv2N3eHdrD0UasHF3JVc+8B6BQKTmNhvhEugyGDK7NKsfBUQrcM4x852/8NZHz7T2UOQIs2bHXkorqw+orXOOO/7xcy5/6XKmLn+Mr884g4UfPpLkEcrB2lxSxbhpCyiprGZYn/TagChr9tEDJPGZ1EeqcGQvD8/7XzZVbufkQRczpv9FBH3BButXhEv50VOXMZ+tBJ3j9ooSLjhpQouNb+XOT7hj7s/YEFnDsZkn8sNR/8MxXXu32Pa+bJ+sWkyMIAMHDG3tobRpC9cs4bHXf01Z/GNywtmkpl3F/1x9LR0yQk3u4/FFj/P34uc5d28V3SJDeaPdKv578e94LrMH3Y4ec0DjiiccW8vC+M1onxEkJdDIg22cg9KN4BKQ3Qt8be/z7ea9m1m6cykndjuR9qntv/Ttby0N860H36OsqpoZE7/KH5Y8id9R8ztr5vUHUEAclPnzH+G2j3/P5oCRmUjwwoIPKVz2PL+75PEG2/z8qat4z23hvIpOLA/u4PYVf6ZTVg9OOvZcz/qxeDW3PnUlC8IrSSOVk7udzU/H/qJJ93d/Y9lr/GThzQRIcHw0wvzKeXz08jtMPu7/uCD/awf8utuCDxc8yYOLf8OCNEfIOQpe7cbXz76f4YMGNql9LJ5gyrvP8f4nTxGNlzG488lcf9aNdM1K/h91WbSMzXs3Ux4txzlHXlYe3TO6Y2bsroiypbSKlICP7tlpZKQk90+yPFrOna/9L6/sfJ1YqtEtDoszSsiM/4XIA+v55fduIyu14Q80n9pavoP7l/6BwbEoN42dQce+w9k0Yxbzqu/g5//4b6Z
|
||
|
|
"text/plain": [
|
||
|
|
"<Figure size 432x288 with 1 Axes>"
|
||
|
|
]
|
||
|
|
},
|
||
|
|
"metadata": {
|
||
|
|
"needs_background": "light"
|
||
|
|
},
|
||
|
|
"output_type": "display_data"
|
||
|
|
},
|
||
|
|
{
|
||
|
|
"name": "stdout",
|
||
|
|
"output_type": "stream",
|
||
|
|
"text": [
|
||
|
|
"('z340', 'READ 2024-11-05/X219225.340')\n",
|
||
|
|
"('zdat', 'READ /home/l_samenv/sea/calib_scripts/calib_data/calib2024-11-05_p0_c4.dat')\n",
|
||
|
|
"('zdat', 'READ /home/l_samenv/sea/calib_scripts/calib_data/calib2024-11-05_p1_c4.dat')\n",
|
||
|
|
"('zdat', 'READ /home/l_samenv/sea/calib_scripts/calib_data/calib2024-11-05_p2_c4.dat')\n"
|
||
|
|
]
|
||
|
|
},
|
||
|
|
{
|
||
|
|
"data": {
|
||
|
|
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAYgAAAECCAYAAADkaECYAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDIuMi40LCBodHRwOi8vbWF0cGxvdGxpYi5vcmcv7US4rQAAIABJREFUeJzt3Xt8VeWd7/HPb++dnXuyQxICJIFwpyh4aRTrpaVVK9Zbx2NVeqpOa710pHNmprWjbc85nUtrpz3TM8dWTw8dqaNtpdbpKCqOSiv60noBpVUuItdAAiSQhNyzb+s5fyRgCCskgb1NIt/365UXWc961rN+hEW+e6397LXMOYeIiEh/gZEuQERERicFhIiI+FJAiIiILwWEiIj4UkCIiIgvBYSIiPhSQIiIiC8FhIiI+AqNdAGHmFkucD8QA1Y75345wiWJiJzU0noGYWbLzKzBzNb3a19kZpvNbKuZ3dXbfDXwmHPuFuDKdNYlIiKDS/clpgeBRX0bzCwI3AdcCswFFpvZXKAC2N3bLZnmukREZBBpDQjn3EtAU7/ms4GtzrntzrkYsBy4CqilJyTSXpeIiAxuJN6DKOf9MwXoCYYFwL3AT8zsMuDJgTY2s1uBO4FIOBwuqaqqSmOpIiIfPu+9994B51zpYP1GzZvUzrkO4ItD6LcUWApQXV3t1q5dm+7SREQ+VMysZij9RuJSTh1Q2We5ordNRERGkZEIiDXATDObamZh4HpgxQjUISIix5Duaa6PAK8Cs82s1sxuds4lgCXAs8Am4FHn3IZ01iEiIsOX1vcgnHOLB2hfCaxM575FROTEaDqpiIj4UkCIiIgvBYSIiPhSQIiIiC8FhIiI+FJAiIiILwWEiIj4UkCIiIgvBYSIiPhSQIiIiC8FhIiI+FJAiIiILwWEiIj4UkCIiIgvBYSIiPhSQIiIiK9RFRBmNs3MHjCzx0a6FhGRk13KAsLMlplZg5mt79e+yMw2m9lWM7vrWGM457Y7525OVU0iInL8UvnI0QeBnwAPHWowsyBwH3AxUAusMbMVQBC4p9/2X3LONaSwHhEROQEpCwjn3EtmVtWv+Wxgq3NuO4CZLQeucs7dA1yeqn2LiEjqpfIMwk85sLvPci2wYKDOZlYMfBc4w8zu7g2S/n1uBe4EIpFIhNWrV6e2YhERAdIfEMPinGsEbh+kz1JgKUB1dbVbuHDhB1CZiMjJJ92zmOqAyj7LFb1tIiIyyqU7INYAM81sqpmFgeuBFWnep4iIpEAqp7k+ArwKzDazWjO72TmXAJYAzwKbgEedcxtStU8REUmfVM5iWjxA+0pgZar2IyIiH4xR9UlqEREZPRQQIiLiSwEhIiK+FBAiIuJLASEiIr4UECIi4ksBISIivhQQIiLiSwEhIiK+FBAiIuJLASEiIr4UECIi4ksBISIivhQQIiLiSwEhIiK+FBAiIuJLASEiIr5S9kS5VDCzzwKXAQXAA86550a4JBGRk1Yqn0m9zMwazGx9v/ZFZrbZzLaa2V3HGsM597hz7hbgduC6VNUmIiLDl8oziAeBnwAPHWowsyBwH3AxUAusMbMVQBC4p9/2X3LONfR+/+3e7UREZISkLCCccy+ZWVW/5rOBrc657QBmthy4yjl3D3B5/zHMzIDvA884597y24+Z3QrcCUQikQirV69O1V9BRET6SPd7EOXA7j7LtcCCY/T/KnARUGhmM5xzP+3fwTm3FFgKUF1d7RYuXJi6akVE5LBR9Sa1c+5e4N6RrkNERNI/zbUOqOyzXNHbJiIio1y6A2INMNPMpppZGLgeWJHmfYqISAqkcprrI8CrwGwzqzWzm51zCWAJ8CywCXjUObchVfsUEZH0SeUspsUDtK8EVqZqPyIi8sHQrTZERMSXAkJERHwpIERExJcCQkREfCkgRETElwJCRER8KSBERMSXAkJERHwpIERExJcCQkREfCkgRETElwJCROQkct8LW4fcVwEhInISefJPe4bcVwEhInKScM6xq6lzyP0VECIiJ4n97VG6kgeH3F8BISJykqhp7OSUin8ecv9RFRBm9hEz+6mZPWZmXxnpekRExqKklyTpJY9q33mgg/Zw15DHSeUjR5eZWYOZre/XvsjMNpvZVjO761hjOOc2OeduB64FzktVbSIiJ4OEl+DH637Mhb+5kNtW3XbU+h0H9tMYsiGPl8oziAeBRX0bzCwI3AdcCswFFpvZXDObZ2ZP9fsa37vNlcDT6DGlIiLD8tzO51j69lKyQlm8vvd16jvqj1i/t37dsMZL5TOpXzKzqn7NZwNbnXPbAcxsOXCVc+4e4PIBxlkBrDCzp4FfnUhNzb9+lNannjqRIURExoyC9lq+05FkbnEeGxoT7FpxE9054ym4/HKKrruW9rZ3IH/o46UsIAZQDuzus1wLLBios5ktBK4GMhngDMLMbgXuBCKRSITVq1cPuPOiX/yCUG0tiYqKYRcuIjLWtMVbCVuYWHuMDAuzv2M/OdsbOHjwIH8qG09XbMewxkt3QAyLc241sHqQPkuBpQDV1dVu4cKFA/ateWAZRCJMefih1BUpIjJKfe03F1E9oZorLvg+v1v7Ix7e+DCPPTeXkAWZcva5JN/+74xLekMeL92zmOqAyj7LFb1tIiKSQi3RFuo765kZmQnApyZ/ioRL0BptAXqmuHaF26nwhn5ekO6AWAPMNLOpZhYGrgdWpHmfIiInnW0HtwEws6gnIGaPmw1ANBkFYGdjB00ZcSoyCoc8ZiqnuT4CvArMNrNaM7vZOZcAlgDPApuAR51zG1K1TxER6bGleQvA4TOI7FA2OaEc4l4CgO0NDTSHjGn5k4Y8ZipnMS0eoH0lmrIqIpJWWw5uIT8jnwm5Ew63jcsaR8JrA2BP/ZsAVBXNGPKYo+qT1CIicny2NG9hRtEMzN7/INy47HHEvTgArW09n2GeUjpvyGMqIEREPgS2tWxjemT6EW09ZxA9l5i6YjUATC4f8JMGR1FAiIiMcZ3xTlqiLZTnlR/RXpxVTNyLk3SOWHA/JYkkOZEpQx5XASEiMsbVd/bcUqMsp+yI9kNnEN3xJF0Z7VS4DLCRuReTiIiMgEMB0fcNaoDi7GIcjs5YjMaMxLCmuIICQkRkzDt0Uz6/MwiAzlgXLSGjqt8lqMEoIERExrh9HfsAGJ8z/oj2QwHRFe/5NPX04pnDGlcBISIyxtV31lOUWURWKOuI9kMBEfc6AJg8fv6wxlVAiIiMcfWd9ZTllh3VfiggYhYDoHLS0Ke4ggJCRGTM29exjwk5E45qj2RGAMMzj7JEkuzC4T36QAEhIjLGDXQGEQwECVnPHZUmufCwpriCAkJEZEzrSnTREm05agbTIQELAgx7iisoIERExrRDU1z7fwbikEDvPVmHO8W1Z1sRERmzBvoU9WG9D5CbUTJr2GMrIERExrDDAeHzHgSAJR0AU8afNuyxFRAiImPYoQ/J+Z1BJD1HKOkRdFAxjLu4HpKyBwaJiMgHb1frLoqzio/6kBz0PGa0JOmRjyMzf+Kwx9YZhIjIGLajZQfTItN8123a00I+XQTD+cOe4goKCBGRIXHOQXMN7HwFuls+sP3GkjFaov77c86xvWU70wr9A6Ju52YyiRPIHv4UV9AlppNKY3sUz0FpfuZIlyIyZtS117HstX9iRe0LZHtJpsfiXNHeSXn4XPIu/xdOmTb86aNDEffi/Hz9z/nlpl8SS8a491P3ctaEs47os79rP+3xdqYWTvUdI7D7tZ4/s07CgNh7cAdPrf0x86dfSmVk+hHPYh2NnHPsaWrhDxv+k7drV7O7u4YD3kG6iOLhke0CFHuZlAciVGVNYHbBFOaUTKPl4F621r/HlrZddFGIl3seN3/udiYWZg9pv9FEkntXrmRNzQ+JBuJk2nSu/djXueaM09P8Nx7bPM/hOYeLd1K/503+uPsPbG3cRKKziVBXMwVdrVTFPcoTubRmzqC+uJquqRczfXIFcybkkxMe0/+9xrxoIsm7e9vY19pNWUEWp04qIBQc3kWTjY0buWXljXQlu/lEe5I
|
||
|
|
"text/plain": [
|
||
|
|
"<Figure size 432x288 with 1 Axes>"
|
||
|
|
]
|
||
|
|
},
|
||
|
|
"metadata": {
|
||
|
|
"needs_background": "light"
|
||
|
|
},
|
||
|
|
"output_type": "display_data"
|
||
|
|
}
|
||
|
|
],
|
||
|
|
"source": [
|
||
|
|
"# compare calibration files from points 0,1,2\n",
|
||
|
|
"# with the optimized (average or selection of best)\n",
|
||
|
|
"\n",
|
||
|
|
"for sensor in run.sensors:\n",
|
||
|
|
" r0, t0 = read_curve(sensor.outputpath, sensor.outputkind)\n",
|
||
|
|
" plt.figure()\n",
|
||
|
|
" dif = [0,0,0]\n",
|
||
|
|
" for j in range(3):\n",
|
||
|
|
" rr, rt = read_curve(sensor.caldat_file[j], 'zdat')\n",
|
||
|
|
" rc, tc = make_calib(run.rref, run.tref, rr, rt, run.t_points)\n",
|
||
|
|
" dif[j] = compare_calib(r0, t0, rc, tc)\n",
|
||
|
|
" plt.plot(t0, dif[j], '-')\n",
|
||
|
|
" plt.xscale('log')\n",
|
||
|
|
" plt.yscale('symlog', linthreshy=0.001)\n",
|
||
|
|
" plt.grid(True, axis='y')\n",
|
||
|
|
" plt.axis([min(t0),max(t0),-1,1])\n",
|
||
|
|
" tmin,tmax,dmax=1.4,310,0.001\n",
|
||
|
|
" plt.plot([tmin,tmax,tmax,tmin,tmin], [-dmax,-dmax,dmax,dmax,-dmax], '-')\n",
|
||
|
|
" #plt.legend(['dif1','dif2','dif3','est'])\n",
|
||
|
|
" plt.show()"
|
||
|
|
]
|
||
|
|
},
|
||
|
|
{
|
||
|
|
"cell_type": "code",
|
||
|
|
"execution_count": 13,
|
||
|
|
"metadata": {},
|
||
|
|
"outputs": [
|
||
|
|
{
|
||
|
|
"name": "stdout",
|
||
|
|
"output_type": "stream",
|
||
|
|
"text": [
|
||
|
|
"X219219.340 X219225.340 X219231.340\r\n"
|
||
|
|
]
|
||
|
|
}
|
||
|
|
],
|
||
|
|
"source": []
|
||
|
|
},
|
||
|
|
{
|
||
|
|
"cell_type": "code",
|
||
|
|
"execution_count": null,
|
||
|
|
"metadata": {},
|
||
|
|
"outputs": [],
|
||
|
|
"source": []
|
||
|
|
},
|
||
|
|
{
|
||
|
|
"cell_type": "code",
|
||
|
|
"execution_count": null,
|
||
|
|
"metadata": {},
|
||
|
|
"outputs": [],
|
||
|
|
"source": [
|
||
|
|
"run.rref"
|
||
|
|
]
|
||
|
|
},
|
||
|
|
{
|
||
|
|
"cell_type": "code",
|
||
|
|
"execution_count": null,
|
||
|
|
"metadata": {},
|
||
|
|
"outputs": [],
|
||
|
|
"source": [
|
||
|
|
"run.tref"
|
||
|
|
]
|
||
|
|
},
|
||
|
|
{
|
||
|
|
"cell_type": "code",
|
||
|
|
"execution_count": null,
|
||
|
|
"metadata": {},
|
||
|
|
"outputs": [],
|
||
|
|
"source": []
|
||
|
|
},
|
||
|
|
{
|
||
|
|
"cell_type": "code",
|
||
|
|
"execution_count": null,
|
||
|
|
"metadata": {},
|
||
|
|
"outputs": [],
|
||
|
|
"source": []
|
||
|
|
},
|
||
|
|
{
|
||
|
|
"cell_type": "code",
|
||
|
|
"execution_count": null,
|
||
|
|
"metadata": {},
|
||
|
|
"outputs": [],
|
||
|
|
"source": []
|
||
|
|
},
|
||
|
|
{
|
||
|
|
"cell_type": "code",
|
||
|
|
"execution_count": null,
|
||
|
|
"metadata": {},
|
||
|
|
"outputs": [],
|
||
|
|
"source": []
|
||
|
|
},
|
||
|
|
{
|
||
|
|
"cell_type": "code",
|
||
|
|
"execution_count": null,
|
||
|
|
"metadata": {},
|
||
|
|
"outputs": [],
|
||
|
|
"source": []
|
||
|
|
},
|
||
|
|
{
|
||
|
|
"cell_type": "code",
|
||
|
|
"execution_count": null,
|
||
|
|
"metadata": {},
|
||
|
|
"outputs": [],
|
||
|
|
"source": []
|
||
|
|
},
|
||
|
|
{
|
||
|
|
"cell_type": "code",
|
||
|
|
"execution_count": null,
|
||
|
|
"metadata": {},
|
||
|
|
"outputs": [],
|
||
|
|
"source": []
|
||
|
|
},
|
||
|
|
{
|
||
|
|
"cell_type": "code",
|
||
|
|
"execution_count": null,
|
||
|
|
"metadata": {},
|
||
|
|
"outputs": [],
|
||
|
|
"source": []
|
||
|
|
},
|
||
|
|
{
|
||
|
|
"cell_type": "code",
|
||
|
|
"execution_count": null,
|
||
|
|
"metadata": {},
|
||
|
|
"outputs": [],
|
||
|
|
"source": []
|
||
|
|
},
|
||
|
|
{
|
||
|
|
"cell_type": "code",
|
||
|
|
"execution_count": null,
|
||
|
|
"metadata": {},
|
||
|
|
"outputs": [],
|
||
|
|
"source": []
|
||
|
|
},
|
||
|
|
{
|
||
|
|
"cell_type": "code",
|
||
|
|
"execution_count": null,
|
||
|
|
"metadata": {},
|
||
|
|
"outputs": [],
|
||
|
|
"source": []
|
||
|
|
}
|
||
|
|
],
|
||
|
|
"metadata": {
|
||
|
|
"kernelspec": {
|
||
|
|
"display_name": "Python 2.7",
|
||
|
|
"language": "python",
|
||
|
|
"name": "python2"
|
||
|
|
},
|
||
|
|
"language_info": {
|
||
|
|
"codemirror_mode": {
|
||
|
|
"name": "ipython",
|
||
|
|
"version": 2
|
||
|
|
},
|
||
|
|
"file_extension": ".py",
|
||
|
|
"mimetype": "text/x-python",
|
||
|
|
"name": "python",
|
||
|
|
"nbconvert_exporter": "python",
|
||
|
|
"pygments_lexer": "ipython2",
|
||
|
|
"version": "2.7.5"
|
||
|
|
}
|
||
|
|
},
|
||
|
|
"nbformat": 4,
|
||
|
|
"nbformat_minor": 2
|
||
|
|
}
|