|
|
{ "cells": [ { "attachments": {}, "cell_type": "markdown", "metadata": {}, "source": [ "# Import supporting package" ] }, { "cell_type": "code", "execution_count": 1, "metadata": {}, "outputs": [], "source": [ "import xarray as xr\n", "import numpy as np\n", "import copy\n", "\n", "from uncertainties import ufloat\n", "from uncertainties import unumpy as unp\n", "from uncertainties import umath\n", "import random\n", "import matplotlib.pyplot as plt\n", "plt.rcParams['font.size'] = 12\n", "\n", "from DataContainer.ReadData import read_hdf5_file\n", "from Analyser.ImagingAnalyser import ImageAnalyser\n", "from Analyser.FitAnalyser import FitAnalyser\n", "from Analyser.FitAnalyser import NewFitModel, DensityProfileBEC2dModel\n", "from ToolFunction.ToolFunction import *\n", "\n", "from scipy.optimize import curve_fit\n", "\n", "from ToolFunction.HomeMadeXarrayFunction import errorbar, dataarray_plot_errorbar\n", "xr.plot.dataarray_plot.errorbar = errorbar\n", "xr.plot.accessor.DataArrayPlotAccessor.errorbar = dataarray_plot_errorbar\n", "\n", "imageAnalyser = ImageAnalyser()\n", "\n", "# %matplotlib notebook" ] }, { "attachments": {}, "cell_type": "markdown", "metadata": {}, "source": [ "## Start a client for parallel computing" ] }, { "cell_type": "code", "execution_count": 2, "metadata": {}, "outputs": [ { "data": { "text/html": [ "<div>\n", " <div style=\"width: 24px; height: 24px; background-color: #e1e1e1; border: 3px solid #9D9D9D; border-radius: 5px; position: absolute;\"> </div>\n", " <div style=\"margin-left: 48px;\">\n", " <h3 style=\"margin-bottom: 0px;\">Client</h3>\n", " <p style=\"color: #9D9D9D; margin-bottom: 0px;\">Client-6b96bd98-14cc-11ee-84cc-80e82ce2fa8e</p>\n", " <table style=\"width: 100%; text-align: left;\">\n", "\n", " <tr>\n", " \n", " <td style=\"text-align: left;\"><strong>Connection method:</strong> Cluster object</td>\n", " <td style=\"text-align: left;\"><strong>Cluster type:</strong> distributed.LocalCluster</td>\n", " \n", " </tr>\n", "\n", " \n", " <tr>\n", " <td style=\"text-align: left;\">\n", " <strong>Dashboard: </strong> <a href=\"http://127.0.0.1:8787/status\" target=\"_blank\">http://127.0.0.1:8787/status</a>\n", " </td>\n", " <td style=\"text-align: left;\"></td>\n", " </tr>\n", " \n", "\n", " </table>\n", "\n", " \n", "\n", " \n", " <details>\n", " <summary style=\"margin-bottom: 20px;\"><h3 style=\"display: inline;\">Cluster Info</h3></summary>\n", " <div class=\"jp-RenderedHTMLCommon jp-RenderedHTML jp-mod-trusted jp-OutputArea-output\">\n", " <div style=\"width: 24px; height: 24px; background-color: #e1e1e1; border: 3px solid #9D9D9D; border-radius: 5px; position: absolute;\">\n", " </div>\n", " <div style=\"margin-left: 48px;\">\n", " <h3 style=\"margin-bottom: 0px; margin-top: 0px;\">LocalCluster</h3>\n", " <p style=\"color: #9D9D9D; margin-bottom: 0px;\">4e95eff3</p>\n", " <table style=\"width: 100%; text-align: left;\">\n", " <tr>\n", " <td style=\"text-align: left;\">\n", " <strong>Dashboard:</strong> <a href=\"http://127.0.0.1:8787/status\" target=\"_blank\">http://127.0.0.1:8787/status</a>\n", " </td>\n", " <td style=\"text-align: left;\">\n", " <strong>Workers:</strong> 8\n", " </td>\n", " </tr>\n", " <tr>\n", " <td style=\"text-align: left;\">\n", " <strong>Total threads:</strong> 128\n", " </td>\n", " <td style=\"text-align: left;\">\n", " <strong>Total memory:</strong> 149.01 GiB\n", " </td>\n", " </tr>\n", " \n", " <tr>\n", " <td style=\"text-align: left;\"><strong>Status:</strong> running</td>\n", " <td style=\"text-align: left;\"><strong>Using processes:</strong> True</td>\n", "</tr>\n", "\n", " \n", " </table>\n", "\n", " <details>\n", " <summary style=\"margin-bottom: 20px;\">\n", " <h3 style=\"display: inline;\">Scheduler Info</h3>\n", " </summary>\n", "\n", " <div style=\"\">\n", " <div>\n", " <div style=\"width: 24px; height: 24px; background-color: #FFF7E5; border: 3px solid #FF6132; border-radius: 5px; position: absolute;\"> </div>\n", " <div style=\"margin-left: 48px;\">\n", " <h3 style=\"margin-bottom: 0px;\">Scheduler</h3>\n", " <p style=\"color: #9D9D9D; margin-bottom: 0px;\">Scheduler-06ee832c-dd0b-4c48-a7b6-ae1064b4a1bf</p>\n", " <table style=\"width: 100%; text-align: left;\">\n", " <tr>\n", " <td style=\"text-align: left;\">\n", " <strong>Comm:</strong> tcp://127.0.0.1:62877\n", " </td>\n", " <td style=\"text-align: left;\">\n", " <strong>Workers:</strong> 8\n", " </td>\n", " </tr>\n", " <tr>\n", " <td style=\"text-align: left;\">\n", " <strong>Dashboard:</strong> <a href=\"http://127.0.0.1:8787/status\" target=\"_blank\">http://127.0.0.1:8787/status</a>\n", " </td>\n", " <td style=\"text-align: left;\">\n", " <strong>Total threads:</strong> 128\n", " </td>\n", " </tr>\n", " <tr>\n", " <td style=\"text-align: left;\">\n", " <strong>Started:</strong> Just now\n", " </td>\n", " <td style=\"text-align: left;\">\n", " <strong>Total memory:</strong> 149.01 GiB\n", " </td>\n", " </tr>\n", " </table>\n", " </div>\n", " </div>\n", "\n", " <details style=\"margin-left: 48px;\">\n", " <summary style=\"margin-bottom: 20px;\">\n", " <h3 style=\"display: inline;\">Workers</h3>\n", " </summary>\n", "\n", " \n", " <div style=\"margin-bottom: 20px;\">\n", " <div style=\"width: 24px; height: 24px; background-color: #DBF5FF; border: 3px solid #4CC9FF; border-radius: 5px; position: absolute;\"> </div>\n", " <div style=\"margin-left: 48px;\">\n", " <details>\n", " <summary>\n", " <h4 style=\"margin-bottom: 0px; display: inline;\">Worker: 0</h4>\n", " </summary>\n", " <table style=\"width: 100%; text-align: left;\">\n", " <tr>\n", " <td style=\"text-align: left;\">\n", " <strong>Comm: </strong> tcp://127.0.0.1:62921\n", " </td>\n", " <td style=\"text-align: left;\">\n", " <strong>Total threads: </strong> 16\n", " </td>\n", " </tr>\n", " <tr>\n", " <td style=\"text-align: left;\">\n", " <strong>Dashboard: </strong> <a href=\"http://127.0.0.1:62923/status\" target=\"_blank\">http://127.0.0.1:62923/status</a>\n", " </td>\n", " <td style=\"text-align: left;\">\n", " <strong>Memory: </strong> 18.63 GiB\n", " </td>\n", " </tr>\n", " <tr>\n", " <td style=\"text-align: left;\">\n", " <strong>Nanny: </strong> tcp://127.0.0.1:62880\n", " </td>\n", " <td style=\"text-align: left;\"></td>\n", " </tr>\n", " <tr>\n", " <td colspan=\"2\" style=\"text-align: left;\">\n", " <strong>Local directory: </strong> C:\\Users\\data\\AppData\\Local\\Temp\\dask-worker-space\\worker-he_941fu\n", " </td>\n", " </tr>\n", "\n", " \n", "\n", " \n", "\n", " </table>\n", " </details>\n", " </div>\n", " </div>\n", " \n", " <div style=\"margin-bottom: 20px;\">\n", " <div style=\"width: 24px; height: 24px; background-color: #DBF5FF; border: 3px solid #4CC9FF; border-radius: 5px; position: absolute;\"> </div>\n", " <div style=\"margin-left: 48px;\">\n", " <details>\n", " <summary>\n", " <h4 style=\"margin-bottom: 0px; display: inline;\">Worker: 1</h4>\n", " </summary>\n", " <table style=\"width: 100%; text-align: left;\">\n", " <tr>\n", " <td style=\"text-align: left;\">\n", " <strong>Comm: </strong> tcp://127.0.0.1:62922\n", " </td>\n", " <td style=\"text-align: left;\">\n", " <strong>Total threads: </strong> 16\n", " </td>\n", " </tr>\n", " <tr>\n", " <td style=\"text-align: left;\">\n", " <strong>Dashboard: </strong> <a href=\"http://127.0.0.1:62924/status\" target=\"_blank\">http://127.0.0.1:62924/status</a>\n", " </td>\n", " <td style=\"text-align: left;\">\n", " <strong>Memory: </strong> 18.63 GiB\n", " </td>\n", " </tr>\n", " <tr>\n", " <td style=\"text-align: left;\">\n", " <strong>Nanny: </strong> tcp://127.0.0.1:62881\n", " </td>\n", " <td style=\"text-align: left;\"></td>\n", " </tr>\n", " <tr>\n", " <td colspan=\"2\" style=\"text-align: left;\">\n", " <strong>Local directory: </strong> C:\\Users\\data\\AppData\\Local\\Temp\\dask-worker-space\\worker-3sii11m3\n", " </td>\n", " </tr>\n", "\n", " \n", "\n", " \n", "\n", " </table>\n", " </details>\n", " </div>\n", " </div>\n", " \n", " <div style=\"margin-bottom: 20px;\">\n", " <div style=\"width: 24px; height: 24px; background-color: #DBF5FF; border: 3px solid #4CC9FF; border-radius: 5px; position: absolute;\"> </div>\n", " <div style=\"margin-left: 48px;\">\n", " <details>\n", " <summary>\n", " <h4 style=\"margin-bottom: 0px; display: inline;\">Worker: 2</h4>\n", " </summary>\n", " <table style=\"width: 100%; text-align: left;\">\n", " <tr>\n", " <td style=\"text-align: left;\">\n", " <strong>Comm: </strong> tcp://127.0.0.1:62916\n", " </td>\n", " <td style=\"text-align: left;\">\n", " <strong>Total threads: </strong> 16\n", " </td>\n", " </tr>\n", " <tr>\n", " <td style=\"text-align: left;\">\n", " <strong>Dashboard: </strong> <a href=\"http://127.0.0.1:62919/status\" target=\"_blank\">http://127.0.0.1:62919/status</a>\n", " </td>\n", " <td style=\"text-align: left;\">\n", " <strong>Memory: </strong> 18.63 GiB\n", " </td>\n", " </tr>\n", " <tr>\n", " <td style=\"text-align: left;\">\n", " <strong>Nanny: </strong> tcp://127.0.0.1:62882\n", " </td>\n", " <td style=\"text-align: left;\"></td>\n", " </tr>\n", " <tr>\n", " <td colspan=\"2\" style=\"text-align: left;\">\n", " <strong>Local directory: </strong> C:\\Users\\data\\AppData\\Local\\Temp\\dask-worker-space\\worker-fs_kjhm5\n", " </td>\n", " </tr>\n", "\n", " \n", "\n", " \n", "\n", " </table>\n", " </details>\n", " </div>\n", " </div>\n", " \n", " <div style=\"margin-bottom: 20px;\">\n", " <div style=\"width: 24px; height: 24px; background-color: #DBF5FF; border: 3px solid #4CC9FF; border-radius: 5px; position: absolute;\"> </div>\n", " <div style=\"margin-left: 48px;\">\n", " <details>\n", " <summary>\n", " <h4 style=\"margin-bottom: 0px; display: inline;\">Worker: 3</h4>\n", " </summary>\n", " <table style=\"width: 100%; text-align: left;\">\n", " <tr>\n", " <td style=\"text-align: left;\">\n", " <strong>Comm: </strong> tcp://127.0.0.1:62931\n", " </td>\n", " <td style=\"text-align: left;\">\n", " <strong>Total threads: </strong> 16\n", " </td>\n", " </tr>\n", " <tr>\n", " <td style=\"text-align: left;\">\n", " <strong>Dashboard: </strong> <a href=\"http://127.0.0.1:62934/status\" target=\"_blank\">http://127.0.0.1:62934/status</a>\n", " </td>\n", " <td style=\"text-align: left;\">\n", " <strong>Memory: </strong> 18.63 GiB\n", " </td>\n", " </tr>\n", " <tr>\n", " <td style=\"text-align: left;\">\n", " <strong>Nanny: </strong> tcp://127.0.0.1:62883\n", " </td>\n", " <td style=\"text-align: left;\"></td>\n", " </tr>\n", " <tr>\n", " <td colspan=\"2\" style=\"text-align: left;\">\n", " <strong>Local directory: </strong> C:\\Users\\data\\AppData\\Local\\Temp\\dask-worker-space\\worker-bc9tduav\n", " </td>\n", " </tr>\n", "\n", " \n", "\n", " \n", "\n", " </table>\n", " </details>\n", " </div>\n", " </div>\n", " \n", " <div style=\"margin-bottom: 20px;\">\n", " <div style=\"width: 24px; height: 24px; background-color: #DBF5FF; border: 3px solid #4CC9FF; border-radius: 5px; position: absolute;\"> </div>\n", " <div style=\"margin-left: 48px;\">\n", " <details>\n", " <summary>\n", " <h4 style=\"margin-bottom: 0px; display: inline;\">Worker: 4</h4>\n", " </summary>\n", " <table style=\"width: 100%; text-align: left;\">\n", " <tr>\n", " <td style=\"text-align: left;\">\n", " <strong>Comm: </strong> tcp://127.0.0.1:62930\n", " </td>\n", " <td style=\"text-align: left;\">\n", " <strong>Total threads: </strong> 16\n", " </td>\n", " </tr>\n", " <tr>\n", " <td style=\"text-align: left;\">\n", " <strong>Dashboard: </strong> <a href=\"http://127.0.0.1:62932/status\" target=\"_blank\">http://127.0.0.1:62932/status</a>\n", " </td>\n", " <td style=\"text-align: left;\">\n", " <strong>Memory: </strong> 18.63 GiB\n", " </td>\n", " </tr>\n", " <tr>\n", " <td style=\"text-align: left;\">\n", " <strong>Nanny: </strong> tcp://127.0.0.1:62884\n", " </td>\n", " <td style=\"text-align: left;\"></td>\n", " </tr>\n", " <tr>\n", " <td colspan=\"2\" style=\"text-align: left;\">\n", " <strong>Local directory: </strong> C:\\Users\\data\\AppData\\Local\\Temp\\dask-worker-space\\worker-xjws4_zw\n", " </td>\n", " </tr>\n", "\n", " \n", "\n", " \n", "\n", " </table>\n", " </details>\n", " </div>\n", " </div>\n", " \n", " <div style=\"margin-bottom: 20px;\">\n", " <div style=\"width: 24px; height: 24px; background-color: #DBF5FF; border: 3px solid #4CC9FF; border-radius: 5px; position: absolute;\"> </div>\n", " <div style=\"margin-left: 48px;\">\n", " <details>\n", " <summary>\n", " <h4 style=\"margin-bottom: 0px; display: inline;\">Worker: 5</h4>\n", " </summary>\n", " <table style=\"width: 100%; text-align: left;\">\n", " <tr>\n", " <td style=\"text-align: left;\">\n", " <strong>Comm: </strong> tcp://127.0.0.1:62915\n", " </td>\n", " <td style=\"text-align: left;\">\n", " <strong>Total threads: </strong> 16\n", " </td>\n", " </tr>\n", " <tr>\n", " <td style=\"text-align: left;\">\n", " <strong>Dashboard: </strong> <a href=\"http://127.0.0.1:62917/status\" target=\"_blank\">http://127.0.0.1:62917/status</a>\n", " </td>\n", " <td style=\"text-align: left;\">\n", " <strong>Memory: </strong> 18.63 GiB\n", " </td>\n", " </tr>\n", " <tr>\n", " <td style=\"text-align: left;\">\n", " <strong>Nanny: </strong> tcp://127.0.0.1:62885\n", " </td>\n", " <td style=\"text-align: left;\"></td>\n", " </tr>\n", " <tr>\n", " <td colspan=\"2\" style=\"text-align: left;\">\n", " <strong>Local directory: </strong> C:\\Users\\data\\AppData\\Local\\Temp\\dask-worker-space\\worker-8j774v71\n", " </td>\n", " </tr>\n", "\n", " \n", "\n", " \n", "\n", " </table>\n", " </details>\n", " </div>\n", " </div>\n", " \n", " <div style=\"margin-bottom: 20px;\">\n", " <div style=\"width: 24px; height: 24px; background-color: #DBF5FF; border: 3px solid #4CC9FF; border-radius: 5px; position: absolute;\"> </div>\n", " <div style=\"margin-left: 48px;\">\n", " <details>\n", " <summary>\n", " <h4 style=\"margin-bottom: 0px; display: inline;\">Worker: 6</h4>\n", " </summary>\n", " <table style=\"width: 100%; text-align: left;\">\n", " <tr>\n", " <td style=\"text-align: left;\">\n", " <strong>Comm: </strong> tcp://127.0.0.1:62898\n", " </td>\n", " <td style=\"text-align: left;\">\n", " <strong>Total threads: </strong> 16\n", " </td>\n", " </tr>\n", " <tr>\n", " <td style=\"text-align: left;\">\n", " <strong>Dashboard: </strong> <a href=\"http://127.0.0.1:62913/status\" target=\"_blank\">http://127.0.0.1:62913/status</a>\n", " </td>\n", " <td style=\"text-align: left;\">\n", " <strong>Memory: </strong> 18.63 GiB\n", " </td>\n", " </tr>\n", " <tr>\n", " <td style=\"text-align: left;\">\n", " <strong>Nanny: </strong> tcp://127.0.0.1:62886\n", " </td>\n", " <td style=\"text-align: left;\"></td>\n", " </tr>\n", " <tr>\n", " <td colspan=\"2\" style=\"text-align: left;\">\n", " <strong>Local directory: </strong> C:\\Users\\data\\AppData\\Local\\Temp\\dask-worker-space\\worker-yeltpw2b\n", " </td>\n", " </tr>\n", "\n", " \n", "\n", " \n", "\n", " </table>\n", " </details>\n", " </div>\n", " </div>\n", " \n", " <div style=\"margin-bottom: 20px;\">\n", " <div style=\"width: 24px; height: 24px; background-color: #DBF5FF; border: 3px solid #4CC9FF; border-radius: 5px; position: absolute;\"> </div>\n", " <div style=\"margin-left: 48px;\">\n", " <details>\n", " <summary>\n", " <h4 style=\"margin-bottom: 0px; display: inline;\">Worker: 7</h4>\n", " </summary>\n", " <table style=\"width: 100%; text-align: left;\">\n", " <tr>\n", " <td style=\"text-align: left;\">\n", " <strong>Comm: </strong> tcp://127.0.0.1:62927\n", " </td>\n", " <td style=\"text-align: left;\">\n", " <strong>Total threads: </strong> 16\n", " </td>\n", " </tr>\n", " <tr>\n", " <td style=\"text-align: left;\">\n", " <strong>Dashboard: </strong> <a href=\"http://127.0.0.1:62928/status\" target=\"_blank\">http://127.0.0.1:62928/status</a>\n", " </td>\n", " <td style=\"text-align: left;\">\n", " <strong>Memory: </strong> 18.63 GiB\n", " </td>\n", " </tr>\n", " <tr>\n", " <td style=\"text-align: left;\">\n", " <strong>Nanny: </strong> tcp://127.0.0.1:62887\n", " </td>\n", " <td style=\"text-align: left;\"></td>\n", " </tr>\n", " <tr>\n", " <td colspan=\"2\" style=\"text-align: left;\">\n", " <strong>Local directory: </strong> C:\\Users\\data\\AppData\\Local\\Temp\\dask-worker-space\\worker-z6tlv47i\n", " </td>\n", " </tr>\n", "\n", " \n", "\n", " \n", "\n", " </table>\n", " </details>\n", " </div>\n", " </div>\n", " \n", "\n", " </details>\n", "</div>\n", "\n", " </details>\n", " </div>\n", "</div>\n", " </details>\n", " \n", "\n", " </div>\n", "</div>" ], "text/plain": [ "<Client: 'tcp://127.0.0.1:62877' processes=8 threads=128, memory=149.01 GiB>" ] }, "execution_count": 2, "metadata": {}, "output_type": "execute_result" } ], "source": [ "from dask.distributed import Client\n", "client = Client(n_workers=8, threads_per_worker=16, processes=True, memory_limit='20GB')\n", "client" ] }, { "attachments": {}, "cell_type": "markdown", "metadata": {}, "source": [ "## Start a client for Mongo DB" ] }, { "cell_type": "code", "execution_count": 3, "metadata": {}, "outputs": [], "source": [ "import pymongo\n", "import xarray_mongodb\n", "\n", "from DataContainer.MongoDB import MongoDB\n", "\n", "mongoClient = pymongo.MongoClient('mongodb://control:DyLab2021@127.0.0.1:27017/?authMechanism=DEFAULT')" ] }, { "attachments": {}, "cell_type": "markdown", "metadata": {}, "source": [ "## Set global path for experiment" ] }, { "cell_type": "code", "execution_count": 4, "metadata": {}, "outputs": [], "source": [ "groupList = [\n", " \"images/MOT_3D_Camera/in_situ_absorption\",\n", " \"images/ODT_1_Axis_Camera/in_situ_absorption\",\n", " \"images/ODT_2_Axis_Camera/in_situ_absorption\",\n", "]\n", "\n", "dskey = {\n", " \"images/MOT_3D_Camera/in_situ_absorption\": \"camera_0\",\n", " \"images/ODT_1_Axis_Camera/in_situ_absorption\": \"camera_1\",\n", " \"images/ODT_2_Axis_Camera/in_situ_absorption\": \"camera_2\",\n", "}\n" ] }, { "cell_type": "code", "execution_count": 5, "metadata": {}, "outputs": [], "source": [ "# img_dir = 'C:/Users/control/DyLab/Experiments/DyBEC/'\n", "# SequenceName = \"Repetition_scan\"\n", "# folderPath = img_dir + SequenceName + \"/\" + get_date()\n", "\n", "# mongoDB = mongoClient[SequenceName]\n", "\n", "# DB = MongoDB(mongoClient, mongoDB, date=get_date())" ] }, { "cell_type": "code", "execution_count": 8, "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAm8AAAHPCAYAAAAFwj37AAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjUuMSwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/YYfK9AAAACXBIWXMAAA9hAAAPYQGoP6dpAACo9UlEQVR4nOzdd1yVZRvA8d/DkKEoAgq4DXOQG0dqKm4y90gt90zNPcJcYDlT00xLTSUzSzNH8jpyoGXuzHKkllsCmbIPcM553j+IkwgoRzkekOv7+fB533M/97nP9dye8PK5l6KqqooQQgghhMgXLMwdgBBCCCGEyDlJ3oQQQggh8hFJ3oQQQggh8hFJ3oQQQggh8hFJ3oQQQggh8hFJ3oQQQggh8hFJ3oQQQggh8hFJ3oQQQggh8hFJ3oQQQggh8hGzJ29xcXFMnTqVtm3bUqJECRRFwc/PL8u6586do3Xr1hQpUgRHR0e6devGjRs3sqy7YsUKqlatio2NDRUrVsTf35/U1NRM9cLCwhg4cCAuLi7Y29vTqFEjDh06lGWbBw8epFGjRtjb2+Pi4sLAgQMJCwt76nsXQgghhDCW2ZO3yMhI1qxZQ3JyMl26dMm23pUrV/D29iYlJYWtW7eyfv16rl27RtOmTQkPD89Qd+7cuYwbN45u3bqxf/9+Ro0axbx58xg9enSGesnJybRq1YpDhw6xfPlydu3ahaurKz4+Phw9ejRD3aNHj/L666/j6urKrl27WL58OQcPHqRVq1YkJyfnWn8IIYQQQjyWamZ6vV7V6/WqqqpqeHi4CqizZ8/OVK9nz56qi4uLGhMTYyi7deuWam1trU6dOtVQFhERodra2qrDhw/P8P65c+eqiqKoly5dMpStXLlSBdTjx48bylJTU1VPT0+1QYMGGd5fv3591dPTU01NTTWU/fLLLyqgrlq16uluXgghhBDCSGZ/8qYoCoqiPLaOVqslMDCQ7t27U7RoUUN5+fLladGiBTt27DCU7du3D41Gw6BBgzK0MWjQIFRVZefOnYayHTt2UKVKFRo1amQos7Kyom/fvpw+fZrg4GAAgoODOXPmDP369cPKyspQt3HjxlSuXDnD5wshhBBCmJLVk6uY3/Xr10lKSqJmzZqZrtWsWZMDBw6g0WiwtbXl4sWLANSoUSNDPXd3d1xcXAzXAS5evEjTpk2zbBPg0qVLlC5d2vCe7D7/l19+yTb25OTkDMOqer2eqKgonJ2dn5i0CiGEEKLgUFWVuLg4SpUqhYVF9s/X8kXyFhkZCYCTk1Oma05OTqiqSnR0NO7u7kRGRmJjY0PhwoWzrJveVnq72bX58Oc+6fMfbvNR8+fPx9/f/3G3J4QQQghhcPfuXcqUKZPt9XyRvKV73JOqh6/ltF5u1X1cG9OmTWPixImG1zExMZQrV46bN2/i4OCQ7fueVmpqKkFBQbRo0QJra+tcb7+gkn41DelX05B+NQ3pV9Mwd7+GxIcwPWgMf2nuY6WqjE1U6PL65yiurzz3WOLi4qhYseIT84N8kbw5OzsDZPmEKyoqCkVRcHR0NNTVaDQkJiZib2+fqa6Xl1eGdrNrE/570vakz8/qiVw6GxsbbGxsMpU7OTllmL+XW1JTU7G3t8fZ2Vl+ueQi6VfTkH41DelX05B+NQ1z9uvJ4ONM+WkMD5QUShSCJZZlqDdgC9hn//e6KaXf/5OmVZl9wUJOeHh4YGdnx4ULFzJdu3DhApUqVcLW1hb4b67bo3VDQ0OJiIigevXqhrIaNWpk2yZgqJv+v9nVfbhNIYQQQuRtqqry5fnPGXFwBA/0KbySnMwWNx/q9d1rtsTNGPkiebOysqJjx45s376duLg4Q/mdO3cICgqiW7duhjIfHx9sbW0JCAjI0EZAQACKomTYS65r165cuXKFU6dOGcq0Wi2bNm2iYcOGlCpVCoDSpUvToEEDNm3ahE6nM9Q9efIkV69ezfD5QgghhMi7krRJ+B4YxeLfV6IHOiVo+LKBP27tl4JlvhiQzBvDpnv37iUhIcGQmF2+fJlt27YB0L59e+zt7fH396d+/fp06NABX19fNBoNs2bNwsXFhUmTJhnacnJyYsaMGcycORMnJyfatm3LmTNn8PPzY+jQoXh6ehrqDh48mJUrV9KzZ08WLFhAyZIlWbVqFVevXuXgwYMZYly4cCFt2rShZ8+ejBo1irCwMHx9falevXqmbUmEEEIIkfcExwczfu8griSGYKmqTElUeKvLdyilaps7NKPkieRt5MiR3L592/D6u+++47vvvgPg5s2bVKhQgapVq3LkyBHee+89evTogZWVFS1btmTx4sWUKFEiQ3vTp0/HwcGBlStXsnjxYtzc3PD19WX69OkZ6tnY2HDo0CGmTp3KmDFjSExMpHbt2uzdu5fmzZtnqOvt7c2ePXuYNWsWHTt2xN7eng4dOvDRRx9lOadNCCGEEHnHqeATTD78Lg/0KTjpdCy2KEX9AVugsLO5QzNankjebt26laN6Xl5emZ6IZWfs2LGMHTv2ifVcXV358ssvc9RmmzZtaNOmTY7q5hadTpflmazZSU1NxcrKCo1Gk2GIVzwbc/arpaWlTI4WQoinpKoqG39fzdJ/h0lfSU5mWek3cPP5KN8Mkz4qf0ZdAKiqSmhoKDExMaiqatT73NzcuHv3rmwCnIvM3a82Nja4uLiYZIWyEEK8qJK0SfgHTeJ///wMQKeEJGY28sO2Tl8zR/ZsJHnLo2JiYnjw4AElSpSgcOHCOU4Y9Ho98fHxFClS5LG7MwvjmKtfVVUlNTWVmJgYw3FtksAJIcST/RP/D+P3DuLPxH/+nd+m8lbn71BK1zF3aM9Mkrc8SFVVwsLCKFq0KC4uLka9V6/Xk5KSgq2trSRvucic/WpnZ4eDgwP37t0jIiJCkjchhHiCrOe3fQuFjfs7Na+Sv93zIJ1Oh06nk7+khYGiKBQrVozk5GSj5kAKIURBkja/bQ0jDg7ngT4Fz+Rkvi3Zhvr99r4wiRvIk7c8SavVAmn72wmRLn3Rgk6nkwUMQgjxCI1Wg1/QJP73z08AdIpPYmajWdjW7W/myHKfZAd5mCw4EA+T74MQQmQt6/ltW1DKeD35zfmQJG9CCCGEyLdO/3OSyYfeJVqfTHGdjiWKG/X7fwtFSpo7NJOR5E0IIYQQ+Y6qqmz64wuWnP8EHVAtOYXl7m1wb/8xWL7YU0tkwYJ4rtLPmD179my2dW7duoWiKJnOp30e0uN79MfS0pKZM2dmGdvx48fx8/PjwYMHzz1eIYQoiDRaDe8fGsOifxO3jglJbPTyxb3jpy984gby5E3kQe7u7pw4cQIPDw+zxbBhwwaqVq1qeK3X63FwcMgytuPHj+Pv78/AgQNxdHQ0Q7RCCFFwhMSHMG7fIP5MCMZSVZmcoOftTt+ilK1n7tCeG0neRJ5jY2PDq6++atYYqlevTr16//0i0Ov1xMbG5onYhBCioDrzzykmHRptmN+2GFcaDNjyQs9vy4oMm+YXqgopCTn7SU3Med1n+THi2C5jZDU06efnh6IoXLp0iT59+lCsWDFcXV0ZPHgwMTExj3SVyqpVq6hduzZ2dnYUL16cHj16cOPGjVyPzc/PjylTpgBQsWJFwzDrkSNHnvmzhBBCpEmf3zbswFCi9clUS07hWxdvGvTfX+ASN5Anb/lHaiLMK/XEahaAo8mD+df7/0Chws/r0wDo3r07vXr1YsiQIVy4cIFp06YBsH79ekOdESNGEBAQwNixY1m4cCFRUVHMmTOHxo0b8/vvv+Pq6vrEz9HpdIb99iDtyVtWhg4dSlRUFCtWrGD79u24u7sD4Onp+Sy3KYQQ4l8arYY5R6awO/gIAB0SkpjdYBq29YaYNzAzkuRN5CtDhgwxPOlq3bo1f//9N+vXr2fdunUoisLJkydZu3YtS5YsYeLEiYb3NW3alMqVK7N06VIWLlz4xM/Jamg0PDw8U1mZMmUoV64cAHXq1KFChQpPeWdCCCEelTa/bTB/Jtz7b35bx80o5RqYOzSzkuQtv7C2T3vS9QR6vZ7YuDiKOjiY/gxOa3vTtp+ "text/plain": [ "<Figure size 640x480 with 1 Axes>" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "fx = [751.71168, 736.56597, 885.97705, 1118.95297, 1538.68901, 3085.17389, 4700.81524, 6434.98779, 8223.48942]\n", "dfx = [1.32859, 1.28606, 0.99619, 0.99134, 1.03403, 0.35958, 0.33758, 0.38790, 0.31241]\n", "x_offset_current = [0, 0.03, 0.05, 0.07, 0.1, 0.2, 0.3, 0.4, 0.515]\n", "\n", "fy = [751.71168, 918.67701, 1129.74687, 1429.82911, 1749.75310, 3089.87897, 4314.90128, 5387.64638, 6666.03670]\n", "dfy = [1.32859, 0.90271, 1.33691, 1.77735, 1.42421, 1.58909, 1.73638, 1.40175, 1.00526]\n", "y_offset_current = [0.0, 0.015, 0.03, 0.05, 0.07, 0.15, 0.22, 0.28, 0.35]\n", "\n", "fz = [751.71168, 1223.54627, 1739.13344, 2364.93284, 2798.24971, 3178.49790, 4275.39905, 5352.17283, 6637.80418, 8288.35264, 9573.59333]\n", "dfz = [1.32859, 0.35554, 0.48471, 0.69762, 0.36873, 0.29413, 0.20667, 0.20818, 0.21978, 0.20285, 0.18495]\n", "z_offset_current = [0.0, 0.03, 0.06, 0.095, 0.119, 0.140, 0.2, 0.259, 0.329, 0.419, 0.489]\n", "\n", "f = fy\n", "df = dfy\n", "offset_current = y_offset_current\n", "\n", "f_fit = fy[4:]\n", "df_fit = dfy[4:]\n", "offset_current_fit = y_offset_current[4:]\n", "\n", "x = np.array(offset_current_fit)\n", "y = np.array(f_fit)\n", "\n", "# Degree of the fitting polynomial\n", "deg = 1\n", "# Parameters from the fit of the polynomial\n", "p = np.polyfit(x, y, deg)\n", "m = p[0] # Gradient\n", "c = p[1] # y-intercept\n", "\n", "#print(f'The fitted straight line has equation y = {m:.1f}x {c:=+6.1f}')\n", "\n", "# Model the data using the parameters of the fitted straight line\n", "y_model = np.polyval(p, x)\n", "\n", "# Create the linear (1 degree polynomial) model\n", "model = np.poly1d(p)\n", "# Fit the model\n", "y_model = model(x)\n", "\n", "# Mean\n", "y_bar = np.mean(y)\n", "# Coefficient of determination, R²\n", "R2 = np.sum((y_model - y_bar)**2) / np.sum((y - y_bar)**2)\n", "\n", "#print(f'R² = {R2:.2f}')\n", "\n", "fitted_SlopeInkHz = m\n", "fitted_offsetInkHz = c\n", "muB = 9.274e-24\n", "hbar = 6.626e-34 / (2 * np.pi)\n", "gJ = 1.24\n", "Slope = (((2 * np.pi * fitted_SlopeInkHz * 1e3)*hbar) / (muB*gJ)) * 1e4\n", "Offset = (((2 * np.pi * fitted_offsetInkHz * 1e3)*hbar) / (muB*gJ)) * 1e4\n", "\n", "def calib_fit(x, B):\n", " alpha = ((2 * np.pi * fitted_SlopeInkHz * 1e3)*hbar) / (muB*gJ)\n", " beta = ((2 * np.pi * fitted_offsetInkHz * 1e3)*hbar) / (muB*gJ)\n", " delta_nu = ((muB * gJ) / hbar) * np.sqrt((B**2-beta**2) + ((alpha * x) + beta)**2)\n", " return delta_nu / (2 * np.pi * 1e3)\n", "\n", "\n", "popt, pcov = curve_fit(calib_fit, offset_current, f, np.array([0.1*1e-4]))\n", "Boffset = popt[0] * 1e4\n", "dBoffset = pcov[0][0]**0.5 * 1e4\n", "\n", "fig = plt.figure()\n", "ax = fig.gca()\n", "plt.clf\n", "plt.errorbar(offset_current, f, yerr=df, fmt='o')\n", "xvals = np.linspace(0, 0.55, 500)\n", "plt.plot(np.array(xvals), p[1] + p[0] * np.array(xvals), label=f'Line Fit')\n", "plt.plot(xvals, calib_fit(xvals, *popt), label=f'Curve Fit')\n", "plt.text(0.25, 2200, f'Line Slope = {Slope:.3f} G/A', fontsize=12)\n", "plt.text(0.25, 1500, f'Line Offset = {Offset:=.3f} G', fontsize=12)\n", "plt.text(0.25, 800, f'Bo= {Boffset:=.3f} +/- {dBoffset:=.3f} G', fontsize=12)\n", "plt.xlabel('Y Offset Coil Current (A)', fontsize=12)\n", "plt.ylabel('Resonance Frequency (kHz)', fontsize=12)\n", "plt.xticks(fontsize=12)\n", "plt.yticks(fontsize=12)\n", "plt.legend(fontsize=12)\n", "#plt.xlim(-0.01, 0.04)\n", "plt.ylim(0, 10000)\n", "plt.grid(visible=1)\n", "plt.tight_layout()\n", "plt.show()" ] }, { "attachments": {}, "cell_type": "markdown", "metadata": {}, "source": [ "# X-comp coil" ] }, { "cell_type": "code", "execution_count": 126, "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAjgAAAGyCAYAAADptr7VAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjUuMSwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/YYfK9AAAACXBIWXMAAA9hAAAPYQGoP6dpAAA7sElEQVR4nO3dfVhU553/8c+EwRmRABnGBwhorCZaFGqb+sDuj4rFGJrGRsfYbLauEWPdVtNqm67RxCeaoHbbbkO9sNtsYmlLHmosaGOtrhK0tiZqanMVaUii1mqIREEFEh0EvX9/sJw4AVKGB4HD+3Vdcxm+5zs397nXOp89M+cehzHGCAAAwEZu6OoJAAAAdDQCDgAAsB0CDgAAsB0CDgAAsB0CDgAAsB0CDgAAsB0CDgAAsB1nV0+gK1y9elXvvvuubrzxRjkcjq6eDgAAaAVjjGpqahQbG6sbbvj4azS9MuC8++67io+P7+ppAACANjh16pTi4uI+tqdXBpwbb7xRUsMCRUREdPFsAABAa1RXVys+Pt56Hf84vTLgNL4tFRERQcABAKCHac3HS/iQMQAAsB0CDgAAsB0CDgAAsB0CDgAAsB0CDgAAsB0CDgAAsB0CDgAAsB0CDgAAsB0CDgAAsJ1euZMxAADoHFeuSPv2SadPSzExUkqKFBJy/edBwAEAAB0iP19atEh6550Pa3FxUna25PNd37nwFhUAAGi3/Hzp3nsDw40klZU11PPzr+98CDgAAKBdrlxpuHJjTNNjjbXFixv6rhcCDgAAaJd9+5peubmWMdKpUw191wsBBwAAtMvp0x3b1xEIOAAAoF1iYjq2ryMQcAAAQLukpDTcLeVwNH/c4ZDi4xv6rhcCDgAAaJeQkIZbwaWmIafx5yefvL774RBwAABAu/l80ubN0s03B9bj4hrq3X4fnD//+c+aNm2aYmNjFRYWppEjR+q73/2uLl68GNB3+PBhTZ48WeHh4YqKipLP59Px48ebHXP9+vUaOXKkXC6Xhg4dqszMTNXV1TXpO3PmjObMmSOv16uwsDAlJyersLAw2FMAAACdwOeTTpyQioqk555r+PNvf7v+4UYKMuD89a9/1T/90z/pxIkTevLJJ7Vt2zb9y7/8i7773e/q/vvvt/pKS0uVmpqqy5cva9OmTdq4caPeeustpaSk6OzZswFjZmVladGiRfL5fNq5c6cWLFigNWvWaOHChQF9tbW1SktLU2FhobKzs7V161YNHDhQ6enp2rt3bzuWAAAAdJSQECk1Vbr//oY/u+JrGiRJJgiPPfaYkWSOHj0aUJ8/f76RZM6dO2eMMWbmzJnG6/Waqqoqq+fEiRMmNDTULFmyxKpVVFQYt9tt5s+fHzBeVlaWcTgcpqSkxKrl5OQYSWb//v1Wra6uziQkJJhx48YFcxqmqqrKSAqYHwAA6N6Cef0O6gpOaGioJCkyMjKgHhUVpRtuuEF9+vRRfX29tm3bphkzZigiIsLqGTJkiCZNmqSCggKrtmPHDvn9fmVkZASMl5GRIWOMtmzZYtUKCgo0YsQIJScnWzWn06lZs2bp4MGDKisrC+ZUAACAjQUVcB544AFFRUXp61//uo4fP66amhpt27ZNP/3pT7Vw4UL169dPx44d06VLl5SUlNTk+UlJSTp69Kj8fr8k6ciRI5KkxMTEgL6YmBh5vV7reGNvS2NKUklJSYvzrq2tVXV1dcADAADYV1AB55ZbbtErr7yiI0eOaNiwYYqIiNDUqVP1wAMPKPv/7g+rrKyUJHk8nibP93g8Msbo/PnzVq/L5VK/fv2a7W0cq7G3pTGv/b3NWbt2rSIjI61HfHx8EGcNAAB6mqACzokTJzR16lRFR0dr8+bN2rt3r/7zP/9Tubm5mjdvXkCvo6Xdfj5yrLV9wfZea9myZaqqqrIep06darEXAAD0fM5gmpcuXarq6mq9/vrr1lWXz33uc/J6vZo7d65mz56tQYMGSWr+isq5c+fkcDgUFRUlSYqOjpbf79fFixcVFhbWpPf222+3fo6Ojm5xTKn5K0aNXC6XXC5XMKcKAAB6sKCu4Lz++utKSEho8pbS2LFjJcl666pv374qLi5u8vzi4mINHz5cbrdb0oefvflob3l5uSoqKjR69GirlpiY2OKYkgJ6AQBA7xZUwImNjVVJSYnef//9gPorr7wiSYqLi5PT6dTUqVOVn5+vmpoaq+fkyZMqKiqS75rdftLT0+V2u5WbmxswXm5urhwOh6ZNm2bVpk+frtLSUh04cMCq1dfXKy8vT+PHj1dsbGwwpwIAAOwsmPvPt27dahwOh5kwYYL51a9+ZQoLC01WVpYJDw83CQkJpra21hhjzBtvvGHCw8PN5z73ObN9+3aTn59vRo8ebWJjY82ZM2cCxnziiSeMw+Ewjz76qNmzZ4/5/ve/b1wul/nqV78a0Of3+82oUaNMfHy8efbZZ82uXbvM9OnTjdPpNHv27AnmNNgHBwCAHiiY1++gAo4xxrz88stmypQpZtCgQaZv377mtttuMw8//LCpqKgI6HvttddMWlqaCQsLMxEREWbatGlNNghslJ2dbW677TbTp08fM3jwYLNq1Spz+fLlJn3l5eVm9uzZxuPxGLfbbSZMmGB27doV7CkQcAAA6IGCef12GGNM115Duv6qq6sVGRmpqqqqgM0IAQBA9xXM6zffJg4AAGyHgAMAAGyHgAMAAGyHgAMAAGyHgAMAAGyHgAMAAGyHgAMAAGyHgAMAAGyHgAMAAGyHgAMAAGyHgAMAAGyHgAMAAGyHgAMAAGyHgAMAAGyHgAMAAGyHgAMAAGyHgAMAAGyHgAMAAGyHgAMAAGyHgAMAAGyHgAMAAGyHgAMAAGyHgAMAAGyHgAMAAGyHgAMAAGyHgAMAAGyHgAMAAGyHgAMAAGyHgAMAAGyHgAMAAGyHgAMAAGyHgAMAAGwnqIAzZ84cORyOFh+vvvqq1Xv48GFNnjxZ4eHhioqKks/n0/Hjx5sdd/369Ro5cqRcLpeGDh2qzMxM1dXVNek7c+aM5syZI6/Xq7CwMCUnJ6uwsDDIUwYAAHbnMMaY1jYfO3ZMZ8+ebVKfOnWqXC6X/v73vyskJESlpaUaN26cxowZo6VLl8rv92vlypU6f/68Xn/9dfXv3996blZWllasWKGlS5dqypQpOnTokJYvX64HHnhATz31lNVXW1urz372s7pw4YLWrVunAQMGKCcnR7/97W+1e/duTZw4sdUnXV1drcjISFVVVSkiIqLVzwMAAF0nqNdv00579uwxkszy5cut2syZM43X6zVVVVVW7cSJEyY0NNQsWbLEqlVUVBi3223mz58fMGZWVpZxOBympKTEquXk5BhJZv/+/Vatrq7OJCQkmHHjxgU156qqKiMpYH4AAKB7C+b1u92fwXnmmWfkcDg0d+5cSVJ9fb22bdumGTNmBKSrIUOGaNKkSSooKLBqO3bskN/vV0ZGRsCYGRkZMsZoy5YtVq2goEAjRoxQcnKyVXM6nZo1a5YOHjyosrKy9p4KAACwiXYFnKqqKm3evFlpaWkaOnSopIa3sS5duqSkpKQm/UlJSTp69Kj8fr8k6ciRI5KkxMTEgL6YmBh5vV7reGNvS2NKUklJSYvzrK2tVXV1dcADAADYV7sCzvPPP69Lly7pwQcftGqVlZWSJI/H06Tf4/HIGKPz589bvS6XS/369Wu2t3Gsxt6Wxrz29zZn7dq1ioyMtB7x8fGtPEMAANATtSvgPPPMM4qOjtb06dObHHM4HC0+79pjre0Ltvday5YtU1VVlfU4depUi70AAKDna3PA+ctf/qLXXntNs2bNksvlsurR0dGSmr+icu7cOTkcDkVFRVm9fr9fFy9ebLb32is20dHRLY4pNX/FqJHL5VJERETAAwAA2FebA84zzzwjSZo3b15AfdiwYerbt6+Ki4ubPKe4uFjDhw+X2+2W9OFnbz7aW15eroqKCo0ePdqqJSYmtjimpIBeAADQu7Up4NTW1iovL0/jxo1rEiycTqemTp2q/Px81dTUWPWTJ0+qqKhIPp/PqqWnp8vtdis3Nzd "text/plain": [ "<Figure size 640x480 with 1 Axes>" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "fx = [751.71168, 736.56597, 885.97705, 1118.95297, 1538.68901, 3085.17389, 4700.81524, 6434.98779, 8223.48942]\n", "dfx = [1.32859, 1.28606, 0.99619, 0.99134, 1.03403, 0.35958, 0.33758, 0.38790, 0.31241]\n", "x_offset_current = [0, 0.03, 0.05, 0.07, 0.1, 0.2, 0.3, 0.4, 0.515]\n", "\n", "# fx = [736.56597, 885.97705, 1118.95297, 1538.68901, 3085.17389, 4700.81524, 6434.98779, 8223.48942]\n", "# dfx = [1.28606, 0.99619, 0.99134, 1.03403, 0.35958, 0.33758, 0.38790, 0.31241]\n", "# x_offset_current = [0.03, 0.05, 0.07, 0.1, 0.2, 0.3, 0.4, 0.515]\n", "\n", "data = xr.DataArray(\n", " data=fx,\n", " dims=['x'],\n", " coords={\n", " 'x': x_offset_current\n", " } \n", ")\n", "\n", "data_std = xr.DataArray(\n", " data=dfx,\n", " dims=['x'],\n", " coords={\n", " 'x': x_offset_current\n", " } \n", ")\n", "\n", "data.plot.errorbar(fmt='ob', yerr=data_std)\n", "\n", "plt.show()" ] }, { "cell_type": "code", "execution_count": 127, "metadata": {}, "outputs": [], "source": [ "def magnetic_field_func(x, b0=0, by0=0, alpha=1):\n", " return 1 / (1e3 * 6.626e-34) * (9.273e-24 * 1.24) * np.sqrt( (b0**2 - by0**2) + (alpha * x + by0)**2 )\n", "\n", "data_quadratic = data\n", "\n", "fitModel_quadratic = NewFitModel(magnetic_field_func)\n", "fitAnalyser_quadratic = FitAnalyser(fitModel_quadratic, fitDim=1)\n", "params_quadratic = fitAnalyser_quadratic.fitModel.make_params()\n", "params_quadratic.add(name=\"b0\", value= 0.3, max=np.inf, min=-np.inf, vary=True)\n", "params_quadratic.add(name=\"by0\", value= 0.07, max=np.inf, min=-np.inf, vary=True)\n", "params_quadratic.add(name=\"alpha\", value= 100, max=np.inf, min=-np.inf, vary=True)\n", "fitResult_quadratic = fitAnalyser_quadratic.fit(data_quadratic, params_quadratic).load()\n", "\n", "fitCurve_quadratic = fitAnalyser_quadratic.eval(fitResult_quadratic, x=np.linspace(0, 0.6, 100), dask=\"parallelized\").load()" ] }, { "cell_type": "code", "execution_count": 128, "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAm8AAAHPCAYAAAAFwj37AAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjUuMSwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/YYfK9AAAACXBIWXMAAA9hAAAPYQGoP6dpAAB6WUlEQVR4nO3dd1hUV/4G8HdoQ+8dbFEDoiAK9oLYokajoiYae4uJxpL4iyVWEmPJulmN0ZhiDZYYRY3GhgbLWrGLir0iivQ6MMOc3x8sExFUZphhKO/neXh2595zz3znZNDXe+85VyKEECAiIiKiCsFA3wUQERERUckxvBERERFVIAxvRERERBUIwxsRERFRBcLwRkRERFSBMLwRERERVSAMb0REREQVCMMbERERUQXC8EZERERUgRjpu4D09HR8/fXXuHjxIi5cuICEhATMmTMHc+fOLdL2/PnzmDJlCk6dOgUjIyO0b98eixcvxltvvVWk7bJly7B8+XLcu3cP7u7uGDZsGL788ksYGxsXahcfH48pU6Zg9+7dyMrKQsOGDTFv3jx06NChSJ8HDx7ErFmzcOnSJZibm6N79+749ttv4ezsXOLPq1Qq8eTJE1hZWUEikZT4OCIiIqrchBBIT0+Hu7s7DAxec35N6Nm9e/eEjY2NaNu2rRg1apQAIObMmVOk3fXr14WVlZVo06aN+Ouvv8S2bdtE/fr1hbu7u4iPjy/Udt68eUIikYjp06eLyMhI8e233woTExMxevToQu1kMplo0KCB8PT0FGFhYeLAgQOiZ8+ewsjISBw+fLhQ28OHDwsjIyPRs2dPceDAAREWFiY8PDxEgwYNhEwmK/HnffTokQDAH/7whz/84Q9/+FPsz6NHj16bJSRC6PfZpgVvL5FIkJCQACcnp2LPvL3//vuIjIzEnTt3YG1tDQB48OAB6tati88++wyLFi0CACQmJsLT0xNDhgzBTz/9pDp+/vz5mDlzJqKjo+Hj4wMAWLFiBcaNG4cTJ06gRYsWAACFQoGGDRvC0tISp0+fVh3ftGlTZGZm4tKlSzAyyj9heeLECbRq1QorVqzAJ598UqLPm5qaCltbWzx69Ej1ObRJLpfjwIED6Ny5c5GzjPR6HDvNcew0x7HTHMdOMxw3zel67NLS0lCtWjWkpKTAxsbmle30ftm0JJcOFQoFdu/ejSFDhhQKPDVq1EBwcDC2b9+uCm/79u2DTCbD8OHDC/UxfPhwzJgxAzt27FCFt+3bt8PLy0sV3ADAyMgIgwYNwpdffonY2Fh4eHggNjYWUVFRWLBggSq4AUDLli3x9ttvY/v27SUObwWf19raWmfhzdzcHNbW1vylVBPHTnMcO81x7DTHsdMMx01zZTV2b8pGeg9vJXHnzh1kZ2fDz8+vyD4/Pz9ERERAJpPB1NQU0dHRAABfX99C7dzc3ODo6KjaDwDR0dFo06ZNsX0CwNWrV+Hh4aE65lXvf/z48VfWnpOTg5ycHNXrtLQ0APlfALlc/srjNFXQpy76ruw4dprj2GmOY6c5jp1mOG6a0/XYlbTfChHeEhMTAQD29vZF9tnb20MIgeTkZLi5uSExMRFSqRQWFhbFti3oq6DfV/X54vu+6f1f7PNlCxYsQGhoaJHtBw4cgLm5+SuPK62IiAid9V3Zcew0x7HTHMdOcxw7zXDcNKerscvKyipRuwoR3gq87jTii/tK2k5bbV/Xx/Tp0/H555+rXhdcz+7cubPOLptGRESgU6dOPB2uJo6d5jh2muPYaY5jpxmOm+Z0PXYFV+fepEKENwcHBwAo9gxXUlISJBIJbG1tVW1lMhmysrKKnNlKSkpCQEBAoX5f1Sfwz5m2N71/cWfkCkilUkil0iLbjY2NdfpLo+v+KzOOneY4dprj2GmOY6cZjpvmdDV2Je2zQizSW7t2bZiZmeHKlStF9l25cgV16tSBqakpgH/udXu57dOnT5GQkIAGDRqotvn6+r6yTwCqtgX/+6q2L/ZJREREpEsVIrwZGRmhR48eCA8PR3p6umr7w4cPERkZiZCQENW2Ll26wNTUFGvXri3Ux9q1ayGRSNCrVy/Vtt69eyMmJqbQkiAKhQJhYWFo1qwZ3N3dAQAeHh5o2rQpwsLCkJeXp2p76tQp3Lhxo9D7ExEREelSubhsunfvXmRmZqqC2bVr17B161YAQLdu3WBubo7Q0FA0adIE3bt3x7Rp0yCTyTB79mw4Ojpi8uTJqr7s7e0xc+ZMzJo1C/b29ujcuTOioqIwd+5cjBo1SrVMCACMGDECy5cvR79+/bBw4UI4OztjxYoVuHHjBg4ePFioxkWLFqFTp07o168fxo4di/j4eEybNg0NGjQosiwJERERka6Ui/D2ySef4MGDB6rXf/zxB/744w8AwL1791CzZk14e3vj8OHDmDp1Kvr27Vvo8VhOTk6F+psxYwasrKywfPlyLF68GK6urpg2bRpmzJhRqJ1UKsWhQ4cwZcoUjB8/HllZWfD398fevXsRFBRUqG27du2wZ88ezJ49Gz169FA9Hutf//pXsfe0EREREelCuQhv9+/fL1G7gICAImfEXmXChAmYMGHCG9u5uLhg3bp1JeqzU6dO6NSpU4naEhEREelChbjnjYiIiIjyMbwRERERVSAMb0REREQVCMMbERERUQXC8EZERERUAkIIPC3Z40d1iuGNiIiI6A3keUrM2XUd3142xOl7SXqtheGNiIiI6DVSs+QYtuYMNkU9hlIAt+Mz9FpPuVjnjYiIiKg8up+QiRHronD3eSbMTQzxYa1cDGxWXa818cwbERERUTFO3U1ErxXHcfd5JtxtTLF5VFP42gt9l8XwRkRERPSyLVGPMHjVaaRkydGwmi12fNoK9dys9F0WAF42JSIiIlLJUwos2heDn4/eBQB093PD4n4NYWpsCLlcrufq8jG8EREREQHIzFFg4uaLOHj9GQBgYoe6mNSxLiQSiZ4rK4zhjYiIiKq82JRsjFp3Ftfj0mBiZIDF/RrivYbu+i6rWAxvREREVKVdeJiM0evPISEjB46WUvw8JACNq9vpu6xXYngjIiKiKuvPS0/wxR+XkKNQwtvVCquGNYGHrZm+y3othjciIiKqcoQQWHroFpYcvAUA6FjPGUv6N4KltPxHo/JfIREREZEWyeR5mLL1Mv689AQAMLpNLUzrWg+GBuVrYsKrMLwRERFRlRGfLsNH68/h4qMUGBlIMK9XA/Rvqt8nJqiL4Y2IiIiqhOtxaRi5NgpPUmWwMTPGj4Mao2VtR32XpTaGNyIiIqr0Dl57hgmbLyArNw9vOVpg1bAmqOVooe+yNMLwRkRERJWWEAK/HruH+XuvQwigZW0H/DgwADbmxvouTWMMb0RERFQp5SqUmL0zGpujHgEABjStjq961oexYcV+tDvDGxEREVU6yZm5+GTDOZy6mwQDCTDjXR+MaFWz3D3qShMMb0RERFSp3HmegZFro3A/MQuWUiMsG9AIwd7O+i5LaxjeiIiIqNI4fjsBn4SdQ5pMAQ9bM6waFghvV2t9l6VVDG9ERERUKWw8/RCzdkYjTynQuLotfh4SCEdLqb7L0jqGNyIiIqrQ8pQC8/66hjXH7wMAevq7Y1EfP5gaG+q3MB1heCMiIqIKK10mx4RNFxB54zkAYHKnt/Fp+zqVYmLCqzC8ERERUYX0KCkLo9adxY1n6TA1NsC/+/njXT83fZelcwxvREREVOGce5CMj9afRWJmLpytpPh1aCD8PG31XVaZYHgjIiKiCmXHhVhM2XoZuXlK1He3xq9DA+FmY6bvssoMwxsRERFVCEqlwHcRN/FD5G0AQGcfFyzp7w9zk6oVZ6rWpyUiIqIKKTs3D5P/uIg9V54CAD5pVxtfdPaCgUHlnZjwKgxvREREVK49S5Nh9PqzuPw4FcaGEiwI8UPfAE99l6U3DG9ERERUbkXHpmLUurN4miaDnbkxfhociKa17PVdll4xvBEREVG5tC86Dp/9fgnZ8jzUcbbE6qFNUN3BXN9l6R3DGxEREZUrQgj8eOQOvt13AwDQpq4jlg9sDGtTYz1XVj4wvBE "text/plain": [ "<Figure size 640x480 with 1 Axes>" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "data_linear = data[3:]\n", "\n", "fitAnalyser_linear = FitAnalyser('Linear', fitDim=1)\n", "params_linear = fitAnalyser_linear.guess(data_linear, dask=\"parallelized\")\n", "fitResult_linear = fitAnalyser_linear.fit(data_linear, params_linear).load()\n", "\n", "fitCurve_linear = fitAnalyser_linear.eval(fitResult_linear, x=np.linspace(0, 0.5, 100), dask=\"parallelized\").load()\n", "\n", "fig = plt.figure()\n", "ax = fig.gca()\n", "\n", "data.plot.errorbar(ax=ax, fmt='ob', yerr=data_std)\n", "# fitCurve_linear.plot.errorbar(ax=ax)\n", "fitCurve_quadratic.plot.errorbar(ax=ax)\n", "\n", "plt.xlabel('X Offset Coil Current (A)', fontsize=12)\n", "plt.ylabel('Resonance Frequency (kHz)', fontsize=12)\n", "plt.xticks(fontsize=12)\n", "plt.yticks(fontsize=12)\n", "# plt.legend(fontsize=12)\n", "#plt.xlim(-0.01, 0.04)\n", "# plt.ylim(0, 10000)\n", "plt.grid(visible=1)\n", "plt.tight_layout()\n", "plt.show()" ] }, { "cell_type": "code", "execution_count": 129, "metadata": {}, "outputs": [ { "data": { "text/html": [ "<div><svg style=\"position: absolute; width: 0; height: 0; overflow: hidden\">\n", "<defs>\n", "<symbol id=\"icon-database\" viewBox=\"0 0 32 32\">\n", "<path d=\"M16 0c-8.837 0-16 2.239-16 5v4c0 2.761 7.163 5 16 5s16-2.239 16-5v-4c0-2.761-7.163-5-16-5z\"></path>\n", "<path d=\"M16 17c-8.837 0-16-2.239-16-5v6c0 2.761 7.163 5 16 5s16-2.239 16-5v-6c0 2.761-7.163 5-16 5z\"></path>\n", "<path d=\"M16 26c-8.837 0-16-2.239-16-5v6c0 2.761 7.163 5 16 5s16-2.239 16-5v-6c0 2.761-7.163 5-16 5z\"></path>\n", "</symbol>\n", "<symbol id=\"icon-file-text2\" viewBox=\"0 0 32 32\">\n", "<path d=\"M28.681 7.159c-0.694-0.947-1.662-2.053-2.724-3.116s-2.169-2.030-3.116-2.724c-1.612-1.182-2.393-1.319-2.841-1.319h-15.5c-1.378 0-2.5 1.121-2.5 2.5v27c0 1.378 1.122 2.5 2.5 2.5h23c1.378 0 2.5-1.122 2.5-2.5v-19.5c0-0.448-0.137-1.23-1.319-2.841zM24.543 5.457c0.959 0.959 1.712 1.825 2.268 2.543h-4.811v-4.811c0.718 0.556 1.584 1.309 2.543 2.268zM28 29.5c0 0.271-0.229 0.5-0.5 0.5h-23c-0.271 0-0.5-0.229-0.5-0.5v-27c0-0.271 0.229-0.5 0.5-0.5 0 0 15.499-0 15.5 0v7c0 0.552 0.448 1 1 1h7v19.5z\"></path>\n", "<path d=\"M23 26h-14c-0.552 0-1-0.448-1-1s0.448-1 1-1h14c0.552 0 1 0.448 1 1s-0.448 1-1 1z\"></path>\n", "<path d=\"M23 22h-14c-0.552 0-1-0.448-1-1s0.448-1 1-1h14c0.552 0 1 0.448 1 1s-0.448 1-1 1z\"></path>\n", "<path d=\"M23 18h-14c-0.552 0-1-0.448-1-1s0.448-1 1-1h14c0.552 0 1 0.448 1 1s-0.448 1-1 1z\"></path>\n", "</symbol>\n", "</defs>\n", "</svg>\n", "<style>/* CSS stylesheet for displaying xarray objects in jupyterlab.\n", " *\n", " */\n", "\n", ":root {\n", " --xr-font-color0: var(--jp-content-font-color0, rgba(0, 0, 0, 1));\n", " --xr-font-color2: var(--jp-content-font-color2, rgba(0, 0, 0, 0.54));\n", " --xr-font-color3: var(--jp-content-font-color3, rgba(0, 0, 0, 0.38));\n", " --xr-border-color: var(--jp-border-color2, #e0e0e0);\n", " --xr-disabled-color: var(--jp-layout-color3, #bdbdbd);\n", " --xr-background-color: var(--jp-layout-color0, white);\n", " --xr-background-color-row-even: var(--jp-layout-color1, white);\n", " --xr-background-color-row-odd: var(--jp-layout-color2, #eeeeee);\n", "}\n", "\n", "html[theme=dark],\n", "body[data-theme=dark],\n", "body.vscode-dark {\n", " --xr-font-color0: rgba(255, 255, 255, 1);\n", " --xr-font-color2: rgba(255, 255, 255, 0.54);\n", " --xr-font-color3: rgba(255, 255, 255, 0.38);\n", " --xr-border-color: #1F1F1F;\n", " --xr-disabled-color: #515151;\n", " --xr-background-color: #111111;\n", " --xr-background-color-row-even: #111111;\n", " --xr-background-color-row-odd: #313131;\n", "}\n", "\n", ".xr-wrap {\n", " display: block !important;\n", " min-width: 300px;\n", " max-width: 700px;\n", "}\n", "\n", ".xr-text-repr-fallback {\n", " /* fallback to plain text repr when CSS is not injected (untrusted notebook) */\n", " display: none;\n", "}\n", "\n", ".xr-header {\n", " padding-top: 6px;\n", " padding-bottom: 6px;\n", " margin-bottom: 4px;\n", " border-bottom: solid 1px var(--xr-border-color);\n", "}\n", "\n", ".xr-header > div,\n", ".xr-header > ul {\n", " display: inline;\n", " margin-top: 0;\n", " margin-bottom: 0;\n", "}\n", "\n", ".xr-obj-type,\n", ".xr-array-name {\n", " margin-left: 2px;\n", " margin-right: 10px;\n", "}\n", "\n", ".xr-obj-type {\n", " color: var(--xr-font-color2);\n", "}\n", "\n", ".xr-sections {\n", " padding-left: 0 !important;\n", " display: grid;\n", " grid-template-columns: 150px auto auto 1fr 20px 20px;\n", "}\n", "\n", ".xr-section-item {\n", " display: contents;\n", "}\n", "\n", ".xr-section-item input {\n", " display: none;\n", "}\n", "\n", ".xr-section-item input + label {\n", " color: var(--xr-disabled-color);\n", "}\n", "\n", ".xr-section-item input:enabled + label {\n", " cursor: pointer;\n", " color: var(--xr-font-color2);\n", "}\n", "\n", ".xr-section-item input:enabled + label:hover {\n", " color: var(--xr-font-color0);\n", "}\n", "\n", ".xr-section-summary {\n", " grid-column: 1;\n", " color: var(--xr-font-color2);\n", " font-weight: 500;\n", "}\n", "\n", ".xr-section-summary > span {\n", " display: inline-block;\n", " padding-left: 0.5em;\n", "}\n", "\n", ".xr-section-summary-in:disabled + label {\n", " color: var(--xr-font-color2);\n", "}\n", "\n", ".xr-section-summary-in + label:before {\n", " display: inline-block;\n", " content: '►';\n", " font-size: 11px;\n", " width: 15px;\n", " text-align: center;\n", "}\n", "\n", ".xr-section-summary-in:disabled + label:before {\n", " color: var(--xr-disabled-color);\n", "}\n", "\n", ".xr-section-summary-in:checked + label:before {\n", " content: '▼';\n", "}\n", "\n", ".xr-section-summary-in:checked + label > span {\n", " display: none;\n", "}\n", "\n", ".xr-section-summary,\n", ".xr-section-inline-details {\n", " padding-top: 4px;\n", " padding-bottom: 4px;\n", "}\n", "\n", ".xr-section-inline-details {\n", " grid-column: 2 / -1;\n", "}\n", "\n", ".xr-section-details {\n", " display: none;\n", " grid-column: 1 / -1;\n", " margin-bottom: 5px;\n", "}\n", "\n", ".xr-section-summary-in:checked ~ .xr-section-details {\n", " display: contents;\n", "}\n", "\n", ".xr-array-wrap {\n", " grid-column: 1 / -1;\n", " display: grid;\n", " grid-template-columns: 20px auto;\n", "}\n", "\n", ".xr-array-wrap > label {\n", " grid-column: 1;\n", " vertical-align: top;\n", "}\n", "\n", ".xr-preview {\n", " color: var(--xr-font-color3);\n", "}\n", "\n", ".xr-array-preview,\n", ".xr-array-data {\n", " padding: 0 5px !important;\n", " grid-column: 2;\n", "}\n", "\n", ".xr-array-data,\n", ".xr-array-in:checked ~ .xr-array-preview {\n", " display: none;\n", "}\n", "\n", ".xr-array-in:checked ~ .xr-array-data,\n", ".xr-array-preview {\n", " display: inline-block;\n", "}\n", "\n", ".xr-dim-list {\n", " display: inline-block !important;\n", " list-style: none;\n", " padding: 0 !important;\n", " margin: 0;\n", "}\n", "\n", ".xr-dim-list li {\n", " display: inline-block;\n", " padding: 0;\n", " margin: 0;\n", "}\n", "\n", ".xr-dim-list:before {\n", " content: '(';\n", "}\n", "\n", ".xr-dim-list:after {\n", " content: ')';\n", "}\n", "\n", ".xr-dim-list li:not(:last-child):after {\n", " content: ',';\n", " padding-right: 5px;\n", "}\n", "\n", ".xr-has-index {\n", " font-weight: bold;\n", "}\n", "\n", ".xr-var-list,\n", ".xr-var-item {\n", " display: contents;\n", "}\n", "\n", ".xr-var-item > div,\n", ".xr-var-item label,\n", ".xr-var-item > .xr-var-name span {\n", " background-color: var(--xr-background-color-row-even);\n", " margin-bottom: 0;\n", "}\n", "\n", ".xr-var-item > .xr-var-name:hover span {\n", " padding-right: 5px;\n", "}\n", "\n", ".xr-var-list > li:nth-child(odd) > div,\n", ".xr-var-list > li:nth-child(odd) > label,\n", ".xr-var-list > li:nth-child(odd) > .xr-var-name span {\n", " background-color: var(--xr-background-color-row-odd);\n", "}\n", "\n", ".xr-var-name {\n", " grid-column: 1;\n", "}\n", "\n", ".xr-var-dims {\n", " grid-column: 2;\n", "}\n", "\n", ".xr-var-dtype {\n", " grid-column: 3;\n", " text-align: right;\n", " color: var(--xr-font-color2);\n", "}\n", "\n", ".xr-var-preview {\n", " grid-column: 4;\n", "}\n", "\n", ".xr-index-preview {\n", " grid-column: 2 / 5;\n", " color: var(--xr-font-color2);\n", "}\n", "\n", ".xr-var-name,\n", ".xr-var-dims,\n", ".xr-var-dtype,\n", ".xr-preview,\n", ".xr-attrs dt {\n", " white-space: nowrap;\n", " overflow: hidden;\n", " text-overflow: ellipsis;\n", " padding-right: 10px;\n", "}\n", "\n", ".xr-var-name:hover,\n", ".xr-var-dims:hover,\n", ".xr-var-dtype:hover,\n", ".xr-attrs dt:hover {\n", " overflow: visible;\n", " width: auto;\n", " z-index: 1;\n", "}\n", "\n", ".xr-var-attrs,\n", ".xr-var-data,\n", ".xr-index-data {\n", " display: none;\n", " background-color: var(--xr-background-color) !important;\n", " padding-bottom: 5px !important;\n", "}\n", "\n", ".xr-var-attrs-in:checked ~ .xr-var-attrs,\n", ".xr-var-data-in:checked ~ .xr-var-data,\n", ".xr-index-data-in:checked ~ .xr-index-data {\n", " display: block;\n", "}\n", "\n", ".xr-var-data > table {\n", " float: right;\n", "}\n", "\n", ".xr-var-name span,\n", ".xr-var-data,\n", ".xr-index-name div,\n", ".xr-index-data,\n", ".xr-attrs {\n", " padding-left: 25px !important;\n", "}\n", "\n", ".xr-attrs,\n", ".xr-var-attrs,\n", ".xr-var-data,\n", ".xr-index-data {\n", " grid-column: 1 / -1;\n", "}\n", "\n", "dl.xr-attrs {\n", " padding: 0;\n", " margin: 0;\n", " display: grid;\n", " grid-template-columns: 125px auto;\n", "}\n", "\n", ".xr-attrs dt,\n", ".xr-attrs dd {\n", " padding: 0;\n", " margin: 0;\n", " float: left;\n", " padding-right: 10px;\n", " width: auto;\n", "}\n", "\n", ".xr-attrs dt {\n", " font-weight: normal;\n", " grid-column: 1;\n", "}\n", "\n", ".xr-attrs dt:hover span {\n", " display: inline-block;\n", " background: var(--xr-background-color);\n", " padding-right: 10px;\n", "}\n", "\n", ".xr-attrs dd {\n", " grid-column: 2;\n", " white-space: pre-wrap;\n", " word-break: break-all;\n", "}\n", "\n", ".xr-icon-database,\n", ".xr-icon-file-text2,\n", ".xr-no-icon {\n", " display: inline-block;\n", " vertical-align: middle;\n", " width: 1em;\n", " height: 1.5em !important;\n", " stroke-width: 0;\n", " stroke: currentColor;\n", " fill: currentColor;\n", "}\n", "</style><pre class='xr-text-repr-fallback'><xarray.Dataset>\n", "Dimensions: ()\n", "Data variables:\n", " b0 object 0.435+/-0.016\n", " by0 object -0.164+/-0.017\n", " alpha object 9.52+/-0.05</pre><div class='xr-wrap' style='display:none'><div class='xr-header'><div class='xr-obj-type'>xarray.Dataset</div></div><ul class='xr-sections'><li class='xr-section-item'><input id='section-05cbf367-574f-45e5-98fd-f2532de532ce' class='xr-section-summary-in' type='checkbox' disabled ><label for='section-05cbf367-574f-45e5-98fd-f2532de532ce' class='xr-section-summary' title='Expand/collapse section'>Dimensions:</label><div class='xr-section-inline-details'></div><div class='xr-section-details'></div></li><li class='xr-section-item'><input id='section-2cac1ee7-d44a-409e-a1d6-5fb3b81d7410' class='xr-section-summary-in' type='checkbox' disabled ><label for='section-2cac1ee7-d44a-409e-a1d6-5fb3b81d7410' class='xr-section-summary' title='Expand/collapse section'>Coordinates: <span>(0)</span></label><div class='xr-section-inline-details'></div><div class='xr-section-details'><ul class='xr-var-list'></ul></div></li><li class='xr-section-item'><input id='section-61f8ea8a-7171-490c-aba2-cefa0f0a3af7' class='xr-section-summary-in' type='checkbox' checked><label for='section-61f8ea8a-7171-490c-aba2-cefa0f0a3af7' class='xr-section-summary' >Data variables: <span>(3)</span></label><div class='xr-section-inline-details'></div><div class='xr-section-details'><ul class='xr-var-list'><li class='xr-var-item'><div class='xr-var-name'><span>b0</span></div><div class='xr-var-dims'>()</div><div class='xr-var-dtype'>object</div><div class='xr-var-preview xr-preview'>0.435+/-0.016</div><input id='attrs-f7be1e84-7c68-4419-ad87-ba1f353393d0' class='xr-var-attrs-in' type='checkbox' disabled><label for='attrs-f7be1e84-7c68-4419-ad87-ba1f353393d0' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-e10443e6-b5c8-4635-808f-afc18505a93d' class='xr-var-data-in' type='checkbox'><label for='data-e10443e6-b5c8-4635-808f-afc18505a93d' title='Show/Hide data repr'><svg class='icon xr-icon-database'><use xlink:href='#icon-database'></use></svg></label><div class='xr-var-attrs'><dl class='xr-attrs'></dl></div><div class='xr-var-data'><pre>array(0.4354862392709215+/-0.016241255661794147, dtype=object)</pre></div></li><li class='xr-var-item'><div class='xr-var-name'><span>by0</span></div><div class='xr-var-dims'>()</div><div class='xr-var-dtype'>object</div><div class='xr-var-preview xr-preview'>-0.164+/-0.017</div><input id='attrs-a24f0b66-670b-43b8-994a-d9e43280d6d9' class='xr-var-attrs-in' type='checkbox' disabled><label for='attrs-a24f0b66-670b-43b8-994a-d9e43280d6d9' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-da2b7bea-f86b-4433-b147-cf1694cf90cb' class='xr-var-data-in' type='checkbox'><label for='data-da2b7bea-f86b-4433-b147-cf1694cf90cb' title='Show/Hide data repr'><svg class='icon xr-icon-database'><use xlink:href='#icon-database'></use></svg></label><div class='xr-var-attrs'><dl class='xr-attrs'></dl></div><div class='xr-var-data'><pre>array(-0.16389847359556292+/-0.017148268573828836, dtype=object)</pre></div></li><li class='xr-var-item'><div class='xr-var-name'><span>alpha</span></div><div class='xr-var-dims'>()</div><div class='xr-var-dtype'>object</div><div class='xr-var-preview xr-preview'>9.52+/-0.05</div><input id='attrs-15c77479-0068-4c20-958e-2dc56241d9d5' class='xr-var-attrs-in' type='checkbox' disabled><label for='attrs-15c77479-0068-4c20-958e-2dc56241d9d5' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-551fa910-8da8-4f36-9eb0-5c44fc299dce' class='xr-var-data-in' type='checkbox'><label for='data-551fa910-8da8-4f36-9eb0-5c44fc299dce' title='Show/Hide data repr'><svg class='icon xr-icon-database'><use xlink:href='#icon-database'></use></svg></label><div class='xr-var-attrs'><dl class='xr-attrs'></dl></div><div class='xr-var-data'><pre>array(9.524149639955006+/-0.051341904930731916, dtype=object)</pre></div></li></ul></div></li><li class=' ], "text/plain": [ "<xarray.Dataset>\n", "Dimensions: ()\n", "Data variables:\n", " b0 object 0.435+/-0.016\n", " by0 object -0.164+/-0.017\n", " alpha object 9.52+/-0.05" ] }, "execution_count": 129, "metadata": {}, "output_type": "execute_result" } ], "source": [ "fitAnalyser_quadratic.get_fit_full_result(fitResult_quadratic) * 1e4" ] }, { "cell_type": "code", "execution_count": 130, "metadata": {}, "outputs": [], "source": [ "val_x = fitAnalyser_quadratic.get_fit_value(fitResult_quadratic) * 1e4\n", "std_x = fitAnalyser_quadratic.get_fit_std(fitResult_quadratic) * 1e4\n", "res_x = fitAnalyser_quadratic.get_fit_full_result(fitResult_quadratic) * 1e4" ] }, { "attachments": {}, "cell_type": "markdown", "metadata": {}, "source": [ "# Y-comp coil" ] }, { "cell_type": "code", "execution_count": 131, "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAjgAAAGyCAYAAADptr7VAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjUuMSwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/YYfK9AAAACXBIWXMAAA9hAAAPYQGoP6dpAAA0OElEQVR4nO3df1TU153/8deE0RmRAhlGDQiaVBOtCps29Qe7y4rVKGtj1TE2J621YtJsq+lqm9RoolFOgqbb7anUg22zWcu25EeNBW2s0SjB1K6tmro5i25pqtaqRIOgARodBLzfP+bLhBEwDiA/Ls/HOXPs3M977tz74ZPO63zm87njMMYYAQAAWOSWrh4AAABARyPgAAAA6xBwAACAdQg4AADAOgQcAABgHQIOAACwDgEHAABYx9nVA+gKV69e1XvvvadPfOITcjgcXT0cAABwA4wxqqmpUUJCgm655frnaHplwHnvvfeUlJTU1cMAAABtcPr0aSUmJl63plcGnE984hOSAjsoOjq6i0cDAABuRHV1tZKSkoKf49fTKwNO49dS0dHRBBwAAHqYG7m8hIuMAQCAdQg4AADAOgQcAABgHQIOAACwDgEHAABYh4ADAACsQ8ABAADWIeAAAADrEHAAAIB1euVKxgAA4OZoaJD27ZPOnpXi46W0NCkiovPHQcABAAAdoqBAWrJEOnPmo7bERCknR/L5OncsfEUFAADaraBAuv/+0HAjSWVlgfaCgs4dDwEHAAC0S0ND4MyNMc23NbYtXRqo6ywEHAAA0C779jU/c9OUMdLp04G6zkLAAQAA7XL2bMfWdQQCDgAAaJf4+I6t6wgEHAAA0C5paYG7pRyOlrc7HFJSUqCusxBwAABAu0REBG4Fl5qHnMbn69d37no4BBwAANBuPp+0ZYs0eHBoe2JioL2z18FhoT8AANAhfD5p5kxWMgYAAJaJiJDS07t6FHxFBQAALETAAQAA1iHgAAAA6xBwAACAdQg4AADAOgQcAABgHQIOAACwDgEHAABYh4ADAACsQ8ABAADWIeAAAADrEHAAAIB1CDgAAMA6BBwAAGAdAg4AALAOAQcAAFiHgAMAAKxDwAEAANYh4AAAAOsQcAAAgHUIOAAAwDoEHAAAYB0CDgAAsA4BBwAAWIeAAwAArEPAAQAA1mlTwPntb3+r6dOn69Zbb1W/fv1055136plnngmpOXz4sKZMmaKoqCjFxsbK5/PpxIkTLfa3YcMGjRw5Ui6XS3fccYeysrJUV1fXrK68vFwLFiyQ1+tVZGSkUlNTVVRU1JYpAAAAi4UdcF566SVNnDhRMTEx+tnPfqYdO3boiSeekDEmWFNaWqr09HRduXJFmzdv1qZNm/Tuu+8qLS1N58+fD+kvOztbS5Yskc/n065du7Ro0SKtXbtWixcvDqmrra3V5MmTVVRUpJycHG3btk2DBg1SRkaG3nrrrTZOHwAAWMmE4cyZM6Z///7mG9/4xnXr5s6da7xer6mqqgq2nTx50vTp08csW7Ys2FZRUWHcbrd55JFHQl6fnZ1tHA6HOXr0aLAtNzfXSDL79+8PttXV1ZlRo0aZcePGhTMNU1VVZSSFjA8AAHRv4Xx+h3UG54UXXtCHH36oJ554otWa+vp6bd++XXPmzFF0dHSwfejQoZo0aZIKCwuDbTt37pTf71dmZmZIH5mZmTLGaOvWrcG2wsJCjRgxQqmpqcE2p9OpefPm6eDBgyorKwtnKgAAwGJhBZzf/OY38ng8Ki0t1d133y2n06mBAwfq61//uqqrqyVJx48f1+XLl5WSktLs9SkpKTp27Jj8fr8k6ciRI5Kk5OTkkLr4+Hh5vd7g9sba1vqUpKNHj7Y67traWlVXV4c8AACAvcIKOGVlZbp06ZLmzp2rBx54QHv27NF3vvMd/exnP9P06dNljFFlZaUkyePxNHu9x+ORMUYXL16UJFVWVsrlcql///4t1jb21VjbWp+N21uzbt06xcTEBB9JSUnhTBsAAPQwznCKr169Kr/fr9WrV2v58uWSpPT0dPXt21dLly5VUVGRIiMjJUkOh6PVfppuu9G6cGubWrFihb797W8Hn1dXVxNyAACwWFhncOLi4iRJ06ZNC2n/53/+Z0mBW8Mba1o6o3LhwgU5HA7FxsYG+/P7/bp06VKLtU3P2MTFxbXap9TyGaNGLpdL0dHRIQ8AAGCvsAJOS9fASAreIn7LLbdo2LBh6tevn0pKSprVlZSUaPjw4XK73ZI+uvbm2tpz586poqJCY8aMCbYlJye32qekkFoAANC7hRVw5syZI0l6/fXXQ9p37NghSZowYYKcTqdmzJihgoIC1dTUBGtOnTql4uJi+Xy+YFtGRobcbrfy8vJC+svLy5PD4dCsWbOCbbNnz1ZpaakOHDgQbKuvr1d+fr7Gjx+vhISEcKYCAABsFu496DNmzDAul8s888wzZvfu3WbdunXG7Xab++67L1jzxz/+0URFRZl/+qd/Mjt27DAFBQVmzJgxJiEhwZSXl4f09+yzzxqHw2GefPJJs3fvXvO9733PuFwu87WvfS2kzu/3m9GjR5ukpCTz4osvmt27d5vZs2cbp9Np9u7dG9YcWAcHAICeJ5zP77ADzqVLl8wTTzxhkpKSjNPpNEOGDDErVqwwfr8/pO7tt982kydPNpGRkSY6OtrMmjXLHDt2rMU+c3JyzF133WX69u1rhgwZYlavXm2uXLnSrO7cuXNm/vz5xuPxGLfbbSZMmGB2794d7hQIOAAA9EDhfH47jGnyGwu9RHV1tWJiYlRVVcUFxwAA9BDhfH7za+IAAMA6BBwAAGAdAg4AALAOAQcAAFiHgAMAAKxDwAEAANYh4AAAAOsQcAAAgHUIOAAAwDoEHAAAYB0CDgAAsA4BBwAAWIeAAwAArEPAAQAA1iHgAAAA6xBwAACAdQg4AADAOgQcAABgHQIOAACwDgEHAABYh4ADAACsQ8ABAADWcXb1AAAAuNkaGqR9+6SzZ6X4eCktTYqI6OpR4WYi4AAArFZQIC1ZIp0581FbYqKUkyP5fF03LtxcfEUFALBWQYF0//2h4UaSysoC7QUFXTMu3HwEHACAlRoaAmdujGm+rbFt6dJAHexDwAEAWGnfvuZnbpoyRjp9OlAH+xBwAABWOnu2Y+vQsxBwAABWio/v2Dr0LAQcAICV0tICd0s5HC1vdzikpKRAHexDwAEAWCkiInAruNQ85DQ+X7+e9XBsRcABAFjL55O2bJEGDw5tT0wMtLMOjr1Y6A8AYDWfT5o5k5WMexsCDgDAehERUnp6V48CnYmvqAAAgHUIOAAAwDoEHAAAYB0CDgAAsA4BBwAAWIeAAwAArEPAAQAA1iHgAAAA6xBwAACAdcIKOHv37pXD4Wjx8fvf/z6k9vDhw5oyZYqioqIUGxsrn8+nEydOtNjvhg0bNHLkSLlcLt1xxx3KyspSXV1ds7ry8nItWLBAXq9XkZGRSk1NVVFRUThTAAAAvUCbfqph7dq1mjRpUkjbmDFjgv+7tLRU6enpuvvuu7V582b5/X49/fTTSktL0zvvvKMBAwYEa7Ozs7Vq1SotX75cU6dO1aFDh7Ry5UqVlZXp+eefD9bV1tZq8uTJ+uCDD5STk6OBAwcqNzdXGRkZ2rNnjyZOnNiWqQAAABuZMBQXFxtJ5tVXX71u3dy5c43X6zVVVVXBtpMnT5o+ffqYZcuWBdsqKiqM2+02jzzySMjrs7OzjcPhMEePHg225ebmGklm//79wba6ujozatQoM27cuHCmYaqqqoykkPEBAIDuLZzP7w6/Bqe+vl7bt2/XnDlzFB0dHWwfOnSoJk2apMLCwmDbzp075ff7lZmZGdJHZmamjDHaunVrsK2wsFAjRoxQampqsM3pdGrevHk6ePCgysrKOnoqAACgh2pTwFm8eLGcTqeio6M1bdo0/fa3vw1uO378uC5fvqyUlJRmr0tJSdGxY8fk9/slSUeOHJEkJScnh9TFx8fL6/UGtzfWttanJB09erQtUwEAABYK6xqcmJgYLVmyROnp6YqLi9OxY8f0ve99T+np6fr1r3+ "text/plain": [ "<Figure size 640x480 with 1 Axes>" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "fy = [751.71168, 918.67701, 1129.74687, 1429.82911, 1749.75310, 3089.87897, 4314.90128, 5387.64638, 6666.03670]\n", "dfy = [1.32859, 0.90271, 1.33691, 1.77735, 1.42421, 1.58909, 1.73638, 1.40175, 1.00526]\n", "y_offset_current = [0.0, 0.015, 0.03, 0.05, 0.07, 0.15, 0.22, 0.28, 0.35]\n", "\n", "# fy = [918.67701, 1129.74687, 1429.82911, 1749.75310, 3089.87897, 4314.90128, 5387.64638, 6666.03670]\n", "# dfy = [0.90271, 1.33691, 1.77735, 1.42421, 1.58909, 1.73638, 1.40175, 1.00526]\n", "# y_offset_current = [0.015, 0.03, 0.05, 0.07, 0.15, 0.22, 0.28, 0.35]\n", "\n", "data = xr.DataArray(\n", " data=fy,\n", " dims=['x'],\n", " coords={\n", " 'x': y_offset_current\n", " } \n", ")\n", "\n", "data_std = xr.DataArray(\n", " data=dfy,\n", " dims=['x'],\n", " coords={\n", " 'x': y_offset_current\n", " } \n", ")\n", "\n", "data.plot.errorbar(fmt='ob', yerr=data_std)\n", "\n", "plt.show()" ] }, { "cell_type": "code", "execution_count": 132, "metadata": {}, "outputs": [], "source": [ "def magnetic_field_func(x, b0=0, by0=0, alpha=1):\n", " return 1 / (1e3 * 6.626e-34) * (9.273e-24 * 1.24) * np.sqrt( (b0**2 - by0**2) + (alpha * x + by0)**2 )\n", "\n", "data_quadratic = data\n", "\n", "fitModel_quadratic = NewFitModel(magnetic_field_func)\n", "fitAnalyser_quadratic = FitAnalyser(fitModel_quadratic, fitDim=1)\n", "params_quadratic = fitAnalyser_quadratic.fitModel.make_params()\n", "params_quadratic.add(name=\"b0\", value= 0.3, max=np.inf, min=-np.inf, vary=True)\n", "params_quadratic.add(name=\"by0\", value= 0.07, max=np.inf, min=-np.inf, vary=True)\n", "params_quadratic.add(name=\"alpha\", value= 100, max=np.inf, min=-np.inf, vary=True)\n", "fitResult_quadratic = fitAnalyser_quadratic.fit(data_quadratic, params_quadratic).load()\n", "\n", "fitCurve_quadratic = fitAnalyser_quadratic.eval(fitResult_quadratic, x=np.linspace(0, 0.6, 100), dask=\"parallelized\").load()" ] }, { "cell_type": "code", "execution_count": 133, "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAm8AAAHQCAYAAAD3Qo21AAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjUuMSwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/YYfK9AAAACXBIWXMAAA9hAAAPYQGoP6dpAAB46ElEQVR4nO3dd3hU1dbH8e+kVwhJKAmhSQ8tEARp0gQBAakCClKsgBT1inBBIFdFuBe9IEXlWlAQBamCghQBUao06b33JEACIf28f+TNSEzAzGSGSfl9niePztn77FmzGGBxztl7mwzDMBARERGRPMHJ0QGIiIiISPapeBMRERHJQ1S8iYiIiOQhKt5ERERE8hAVbyIiIiJ5iIo3ERERkTxExZuIiIhIHqLiTURERCQPcXF0AAVNamoqFy9exNfXF5PJ5OhwREREJJcwDIPY2FiCg4Nxcrr39TUVbw/YxYsXKVWqlKPDEBERkVzq3LlzhISE3LNdxdsD5uvrC6T9whQqVMjm4yclJbF69Wpat26Nq6urzcfPz5Q76yl31lPurKfcWUd5s569cxcTE0OpUqXMtcK9qHh7wNJvlRYqVMhuxZuXlxeFChXSb0oLKXfWU+6sp9xZT7mzjvJmvQeVu797rEoTFkRERETyEBVvIiIiInmIijcRERGRPETFm4iIiEgeouJNREREJA9R8SYiIiKSh6h4ExEREclDVLyJiIiI5CEq3kRERETyEBVvIiIiInmIijcRERGRPETFm4iIiEgeouJNREREJA9R8SYiIiKSTdEJjo5AxZuIiIjI30pOSWXiqiO8u9uZfRduOjQWFW8iIiIi9xF1K4E+n23ns9/OkGyY2Hoq2qHxuDj03UVERERysT/O3+DlOTu5eDMeLzdnepRN5IXG5Rwak668iYiIiGRhwe/n6PbxFi7ejKdcoDcLX6pPWIDh6LB05U1ERETkbonJqfxrxQHmbj0LwGNVi/NBj1p4OsMxB8cGKt5EREREzK7ExDPo613sPHMdkwmGt6zEkBYVcHIykZSU5OjwABVvIiIiIgD8fjqagV/v4lpsAr4eLnzYszbNqxRzdFiZqHgTERGRAs0wDOZuPUPE8oMkpxpULu7LJ33CKRvo7ejQsqTiTURERAqs+KQUxizdz8Kd5wF4omYQ/+5aE2/33Fsi5d7IREREROzo/PU4Bs7dxb4LN3Eywci2VXihyUOYTCZHh3ZfKt5ERESkwNl8PJJXvtlN9O1Eini5Mv3pOjSqEOjosLJFxZuIiIgUGIZh8L9NJ5m48jCpBtQoWZiPetchpIiXo0PLNhVvIiIiUiDcTkhmxKI/+OGPSwB0Cw/hnU7V8XB1dnBkllHxJiIiIvne6cjbvDRnJ0euxOLiZGJsh1D6PFIm1z/flhUVbyIiIpKv/Xz4CsO+3UNsfDJFfd356Jk61C3r7+iwrKbiTURERPKl1FSDaT8fZ8q6oxgG1Cntx0e9wyleyMPRoeWIijcRERHJd2Lik3ht/l7WHroCwDP1SzOuQzXcXJwcHFnOqXgTERGRfOXolVhemrOTU5G3cXNx4p0nq/PUw6UcHZbNqHgTERGRfOPHfZf4x3d7iUtMIbiwBx/3CadmiJ+jw7IpFW8iIiKS56WkGvznpyN8vPEEAA0eCmD607UJ8HF3cGS2p+JNRERE8rTrtxMZ8s1ufj0eCcALTcrxZpsquDjn/efbsqLiTURERPKs/Rdu8vLcnZy/fgdPV2f+3a0mHWoFOzosu1LxJiIiInnSkt3nGbloHwnJqZQJ8OKTPuFUKVHI0WHZnYo3ERERyVOSUlJ594dDzN58GoDmlYsypUdtCnu5OjawB0TFm4iIiOQZV2PjeeXr3Ww/HQ3AkBYVePWxSjg55b1trqyl4k1ERETyhF1nrzNw7k6uxCTg4+7CB0/VonW1Eo4O64FT8SYiIiK5mmEYfLP9HOO+309SikGFYj580iec8kV9HB2aQ6h4ExERkVwrPimFccsOMP/3cwC0rV6C/3SvhY97wS1hCu4nFxERkVzt4o07DJy7k73nb+Jkgn88XpmBTctjMhWc59uyouJNREREcp0tJ6J4Zd4uom4n4uflyoc9a/NopaKODitXUPEmIiIiuYZhGHz+22km/HiIlFSD0KBCfNInnFL+Xo4OLddQ8SYiIiK5wp3EFEYu/oNley4C0Ll2SSZ0roGnm7ODI8tdVLyJiIiIw52NiuPFOb9z+HIsLk4mRj9RlX4Nyxb459uyouJNREREHGrDkasM/WY3MfHJBPq4MePpOtR/KMDRYeVaTo4OIDY2lhEjRtC6dWuKFi2KyWRi/PjxWfbdtWsXjz32GD4+Pvj5+dGlSxdOnjyZZd9p06ZRpUoV3N3dKVeuHBERESQlJWXqd/XqVfr160dgYCBeXl40aNCAdevWZTnm2rVradCgAV5eXgQGBtKvXz+uXr1q9WcXEREpyFJTDab/fIz+s3cQE59MWCk/VgxposLtbzi8eIuKimLWrFkkJCTQqVOne/Y7fPgwzZo1IzExkQULFvD5559z9OhRmjRpwrVr1zL0fffddxk2bBhdunThp59+YtCgQUyYMIHBgwdn6JeQkEDLli1Zt24dU6dOZdmyZRQvXpw2bdqwcePGDH03btxI27ZtKV68OMuWLWPq1KmsXbuWli1bkpCQYLN8iIiIFASx8Um8PHcnk1cfxTDg6fqlmf/SI5Qo7OHo0HI9h982LVOmDNevX8dkMhEZGcmnn36aZb+xY8fi7u7OihUrKFSoEADh4eFUrFiRyZMnM2nSJCCtGHznnXd44YUXmDBhAgDNmjUjKSmJMWPGMHz4cEJDQwH47LPP2L9/P5s3b6ZBgwYANG/enFq1ajFixAi2bdtmfv833niDSpUqsXDhQlxc0tJWrlw5GjVqxOeff87AgQPtkyAREZF85vjVW7w053dOXLuNm7MT/3qyGj3rlXZ0WHmGw6+8mUymv30YMTk5mRUrVtC1a1dz4QZphV/z5s1ZsmSJ+diqVauIj4+nf//+Gcbo378/hmGwdOlS87ElS5ZQuXJlc+EG4OLiQu/evdm+fTsXLlwA4MKFC+zYsYM+ffqYCzeAhg0bUqlSpQzvLyIiIvf204HLdJrxGyeu3aZEIQ/mv/SICjcLObx4y44TJ05w584datasmamtZs2aHD9+nPj4eAD2798PQI0aNTL0CwoKIjAw0Nye3vdeYwIcOHAgw5j36nv3mCIiIpJZSqrB5J+O8NKcndxKSKZeOX+WD2lM7dJFHB1anuPw26bZERUVBYC/v3+mNn9/fwzD4Pr16wQFBREVFYW7uzve3t5Z9k0fK33ce4159/v+3fvfPeZfJSQkZHgmLiYmBoCkpKQsJ1DkVPqY9hg7v1PurKfcWU+5s55yZx1H5O1GXBKvL/yDX46l/X3Zt0Fp3ny8Eq7OTnnq18/eucvuuHmieEt3v9urd7dlt5+t+t5vjPfee4+IiIhMx1evXo2Xl/1Wi16zZo3dxs7vlDvrKXfWU+6sp9xZ50Hl7cJt+OyIM1EJJlydDHo+lEodTrLmp6xXi8gL7JW7uLi4bPXLE8VbQEDalOGsrnBFR0djMpnw8/Mz942PjycuLi5TcRQdHU14eHiGce81Jvx5pe3v3j+rK3LpRo0axWuvvWZ+HRMTQ6lSpWjdunWG5/dsJSkpiTVr1tCqVStcXV1tPn5+ptxZT7mznnJnPeXOOg8yb8v/uMSHSw8Qn5RKSBFPZvSqRWiQ7f/ue1Dsnbv0u3N/J08Ub+XLl8fT05N9+/Zlatu3bx8VKlTAwyNtanH6s2779u2jfv365n6XL18mMjKS6tWrm4/VqFHjnmMC5r7p/923bx/t2rXL1PfuMf/K3d0dd3f3TMddXV3t+pvG3uPnZ8qd9ZQ76yl31lPurGPPvCWlpPLej4f5/LdTADxaqSgf9gzDz8vNLu/3oNkrd9kdM09MWHBxcaFDhw4sXryY2NhY8/GzZ8+yfv16unTpYj7Wpk0bPDw8mD17doYxZs+ejclkyrCWXOfOnTl8+HCGJUGSk5OZO3cu9evXJzg4GICSJUtSr14 "text/plain": [ "<Figure size 640x480 with 1 Axes>" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "data_linear = data[3:]\n", "\n", "fitAnalyser_linear = FitAnalyser('Linear', fitDim=1)\n", "params_linear = fitAnalyser_linear.guess(data_linear, dask=\"parallelized\")\n", "fitResult_linear = fitAnalyser_linear.fit(data_linear, params_linear).load()\n", "\n", "fitCurve_linear = fitAnalyser_linear.eval(fitResult_linear, x=np.linspace(0, 0.5, 100), dask=\"parallelized\").load()\n", "\n", "fig = plt.figure()\n", "ax = fig.gca()\n", "\n", "data.plot.errorbar(ax=ax, fmt='ob', yerr=data_std)\n", "# fitCurve_linear.plot.errorbar(ax=ax)\n", "fitCurve_quadratic.plot.errorbar(ax=ax)\n", "\n", "plt.xlabel('Y Offset Coil Current (A)', fontsize=12)\n", "plt.ylabel('Resonance Frequency (kHz)', fontsize=12)\n", "plt.xticks(fontsize=12)\n", "plt.yticks(fontsize=12)\n", "# plt.legend(fontsize=12)\n", "#plt.xlim(-0.01, 0.04)\n", "# plt.ylim(0, 10000)\n", "plt.grid(visible=1)\n", "plt.tight_layout()\n", "plt.show()" ] }, { "cell_type": "code", "execution_count": 134, "metadata": {}, "outputs": [ { "data": { "text/html": [ "<div><svg style=\"position: absolute; width: 0; height: 0; overflow: hidden\">\n", "<defs>\n", "<symbol id=\"icon-database\" viewBox=\"0 0 32 32\">\n", "<path d=\"M16 0c-8.837 0-16 2.239-16 5v4c0 2.761 7.163 5 16 5s16-2.239 16-5v-4c0-2.761-7.163-5-16-5z\"></path>\n", "<path d=\"M16 17c-8.837 0-16-2.239-16-5v6c0 2.761 7.163 5 16 5s16-2.239 16-5v-6c0 2.761-7.163 5-16 5z\"></path>\n", "<path d=\"M16 26c-8.837 0-16-2.239-16-5v6c0 2.761 7.163 5 16 5s16-2.239 16-5v-6c0 2.761-7.163 5-16 5z\"></path>\n", "</symbol>\n", "<symbol id=\"icon-file-text2\" viewBox=\"0 0 32 32\">\n", "<path d=\"M28.681 7.159c-0.694-0.947-1.662-2.053-2.724-3.116s-2.169-2.030-3.116-2.724c-1.612-1.182-2.393-1.319-2.841-1.319h-15.5c-1.378 0-2.5 1.121-2.5 2.5v27c0 1.378 1.122 2.5 2.5 2.5h23c1.378 0 2.5-1.122 2.5-2.5v-19.5c0-0.448-0.137-1.23-1.319-2.841zM24.543 5.457c0.959 0.959 1.712 1.825 2.268 2.543h-4.811v-4.811c0.718 0.556 1.584 1.309 2.543 2.268zM28 29.5c0 0.271-0.229 0.5-0.5 0.5h-23c-0.271 0-0.5-0.229-0.5-0.5v-27c0-0.271 0.229-0.5 0.5-0.5 0 0 15.499-0 15.5 0v7c0 0.552 0.448 1 1 1h7v19.5z\"></path>\n", "<path d=\"M23 26h-14c-0.552 0-1-0.448-1-1s0.448-1 1-1h14c0.552 0 1 0.448 1 1s-0.448 1-1 1z\"></path>\n", "<path d=\"M23 22h-14c-0.552 0-1-0.448-1-1s0.448-1 1-1h14c0.552 0 1 0.448 1 1s-0.448 1-1 1z\"></path>\n", "<path d=\"M23 18h-14c-0.552 0-1-0.448-1-1s0.448-1 1-1h14c0.552 0 1 0.448 1 1s-0.448 1-1 1z\"></path>\n", "</symbol>\n", "</defs>\n", "</svg>\n", "<style>/* CSS stylesheet for displaying xarray objects in jupyterlab.\n", " *\n", " */\n", "\n", ":root {\n", " --xr-font-color0: var(--jp-content-font-color0, rgba(0, 0, 0, 1));\n", " --xr-font-color2: var(--jp-content-font-color2, rgba(0, 0, 0, 0.54));\n", " --xr-font-color3: var(--jp-content-font-color3, rgba(0, 0, 0, 0.38));\n", " --xr-border-color: var(--jp-border-color2, #e0e0e0);\n", " --xr-disabled-color: var(--jp-layout-color3, #bdbdbd);\n", " --xr-background-color: var(--jp-layout-color0, white);\n", " --xr-background-color-row-even: var(--jp-layout-color1, white);\n", " --xr-background-color-row-odd: var(--jp-layout-color2, #eeeeee);\n", "}\n", "\n", "html[theme=dark],\n", "body[data-theme=dark],\n", "body.vscode-dark {\n", " --xr-font-color0: rgba(255, 255, 255, 1);\n", " --xr-font-color2: rgba(255, 255, 255, 0.54);\n", " --xr-font-color3: rgba(255, 255, 255, 0.38);\n", " --xr-border-color: #1F1F1F;\n", " --xr-disabled-color: #515151;\n", " --xr-background-color: #111111;\n", " --xr-background-color-row-even: #111111;\n", " --xr-background-color-row-odd: #313131;\n", "}\n", "\n", ".xr-wrap {\n", " display: block !important;\n", " min-width: 300px;\n", " max-width: 700px;\n", "}\n", "\n", ".xr-text-repr-fallback {\n", " /* fallback to plain text repr when CSS is not injected (untrusted notebook) */\n", " display: none;\n", "}\n", "\n", ".xr-header {\n", " padding-top: 6px;\n", " padding-bottom: 6px;\n", " margin-bottom: 4px;\n", " border-bottom: solid 1px var(--xr-border-color);\n", "}\n", "\n", ".xr-header > div,\n", ".xr-header > ul {\n", " display: inline;\n", " margin-top: 0;\n", " margin-bottom: 0;\n", "}\n", "\n", ".xr-obj-type,\n", ".xr-array-name {\n", " margin-left: 2px;\n", " margin-right: 10px;\n", "}\n", "\n", ".xr-obj-type {\n", " color: var(--xr-font-color2);\n", "}\n", "\n", ".xr-sections {\n", " padding-left: 0 !important;\n", " display: grid;\n", " grid-template-columns: 150px auto auto 1fr 20px 20px;\n", "}\n", "\n", ".xr-section-item {\n", " display: contents;\n", "}\n", "\n", ".xr-section-item input {\n", " display: none;\n", "}\n", "\n", ".xr-section-item input + label {\n", " color: var(--xr-disabled-color);\n", "}\n", "\n", ".xr-section-item input:enabled + label {\n", " cursor: pointer;\n", " color: var(--xr-font-color2);\n", "}\n", "\n", ".xr-section-item input:enabled + label:hover {\n", " color: var(--xr-font-color0);\n", "}\n", "\n", ".xr-section-summary {\n", " grid-column: 1;\n", " color: var(--xr-font-color2);\n", " font-weight: 500;\n", "}\n", "\n", ".xr-section-summary > span {\n", " display: inline-block;\n", " padding-left: 0.5em;\n", "}\n", "\n", ".xr-section-summary-in:disabled + label {\n", " color: var(--xr-font-color2);\n", "}\n", "\n", ".xr-section-summary-in + label:before {\n", " display: inline-block;\n", " content: '►';\n", " font-size: 11px;\n", " width: 15px;\n", " text-align: center;\n", "}\n", "\n", ".xr-section-summary-in:disabled + label:before {\n", " color: var(--xr-disabled-color);\n", "}\n", "\n", ".xr-section-summary-in:checked + label:before {\n", " content: '▼';\n", "}\n", "\n", ".xr-section-summary-in:checked + label > span {\n", " display: none;\n", "}\n", "\n", ".xr-section-summary,\n", ".xr-section-inline-details {\n", " padding-top: 4px;\n", " padding-bottom: 4px;\n", "}\n", "\n", ".xr-section-inline-details {\n", " grid-column: 2 / -1;\n", "}\n", "\n", ".xr-section-details {\n", " display: none;\n", " grid-column: 1 / -1;\n", " margin-bottom: 5px;\n", "}\n", "\n", ".xr-section-summary-in:checked ~ .xr-section-details {\n", " display: contents;\n", "}\n", "\n", ".xr-array-wrap {\n", " grid-column: 1 / -1;\n", " display: grid;\n", " grid-template-columns: 20px auto;\n", "}\n", "\n", ".xr-array-wrap > label {\n", " grid-column: 1;\n", " vertical-align: top;\n", "}\n", "\n", ".xr-preview {\n", " color: var(--xr-font-color3);\n", "}\n", "\n", ".xr-array-preview,\n", ".xr-array-data {\n", " padding: 0 5px !important;\n", " grid-column: 2;\n", "}\n", "\n", ".xr-array-data,\n", ".xr-array-in:checked ~ .xr-array-preview {\n", " display: none;\n", "}\n", "\n", ".xr-array-in:checked ~ .xr-array-data,\n", ".xr-array-preview {\n", " display: inline-block;\n", "}\n", "\n", ".xr-dim-list {\n", " display: inline-block !important;\n", " list-style: none;\n", " padding: 0 !important;\n", " margin: 0;\n", "}\n", "\n", ".xr-dim-list li {\n", " display: inline-block;\n", " padding: 0;\n", " margin: 0;\n", "}\n", "\n", ".xr-dim-list:before {\n", " content: '(';\n", "}\n", "\n", ".xr-dim-list:after {\n", " content: ')';\n", "}\n", "\n", ".xr-dim-list li:not(:last-child):after {\n", " content: ',';\n", " padding-right: 5px;\n", "}\n", "\n", ".xr-has-index {\n", " font-weight: bold;\n", "}\n", "\n", ".xr-var-list,\n", ".xr-var-item {\n", " display: contents;\n", "}\n", "\n", ".xr-var-item > div,\n", ".xr-var-item label,\n", ".xr-var-item > .xr-var-name span {\n", " background-color: var(--xr-background-color-row-even);\n", " margin-bottom: 0;\n", "}\n", "\n", ".xr-var-item > .xr-var-name:hover span {\n", " padding-right: 5px;\n", "}\n", "\n", ".xr-var-list > li:nth-child(odd) > div,\n", ".xr-var-list > li:nth-child(odd) > label,\n", ".xr-var-list > li:nth-child(odd) > .xr-var-name span {\n", " background-color: var(--xr-background-color-row-odd);\n", "}\n", "\n", ".xr-var-name {\n", " grid-column: 1;\n", "}\n", "\n", ".xr-var-dims {\n", " grid-column: 2;\n", "}\n", "\n", ".xr-var-dtype {\n", " grid-column: 3;\n", " text-align: right;\n", " color: var(--xr-font-color2);\n", "}\n", "\n", ".xr-var-preview {\n", " grid-column: 4;\n", "}\n", "\n", ".xr-index-preview {\n", " grid-column: 2 / 5;\n", " color: var(--xr-font-color2);\n", "}\n", "\n", ".xr-var-name,\n", ".xr-var-dims,\n", ".xr-var-dtype,\n", ".xr-preview,\n", ".xr-attrs dt {\n", " white-space: nowrap;\n", " overflow: hidden;\n", " text-overflow: ellipsis;\n", " padding-right: 10px;\n", "}\n", "\n", ".xr-var-name:hover,\n", ".xr-var-dims:hover,\n", ".xr-var-dtype:hover,\n", ".xr-attrs dt:hover {\n", " overflow: visible;\n", " width: auto;\n", " z-index: 1;\n", "}\n", "\n", ".xr-var-attrs,\n", ".xr-var-data,\n", ".xr-index-data {\n", " display: none;\n", " background-color: var(--xr-background-color) !important;\n", " padding-bottom: 5px !important;\n", "}\n", "\n", ".xr-var-attrs-in:checked ~ .xr-var-attrs,\n", ".xr-var-data-in:checked ~ .xr-var-data,\n", ".xr-index-data-in:checked ~ .xr-index-data {\n", " display: block;\n", "}\n", "\n", ".xr-var-data > table {\n", " float: right;\n", "}\n", "\n", ".xr-var-name span,\n", ".xr-var-data,\n", ".xr-index-name div,\n", ".xr-index-data,\n", ".xr-attrs {\n", " padding-left: 25px !important;\n", "}\n", "\n", ".xr-attrs,\n", ".xr-var-attrs,\n", ".xr-var-data,\n", ".xr-index-data {\n", " grid-column: 1 / -1;\n", "}\n", "\n", "dl.xr-attrs {\n", " padding: 0;\n", " margin: 0;\n", " display: grid;\n", " grid-template-columns: 125px auto;\n", "}\n", "\n", ".xr-attrs dt,\n", ".xr-attrs dd {\n", " padding: 0;\n", " margin: 0;\n", " float: left;\n", " padding-right: 10px;\n", " width: auto;\n", "}\n", "\n", ".xr-attrs dt {\n", " font-weight: normal;\n", " grid-column: 1;\n", "}\n", "\n", ".xr-attrs dt:hover span {\n", " display: inline-block;\n", " background: var(--xr-background-color);\n", " padding-right: 10px;\n", "}\n", "\n", ".xr-attrs dd {\n", " grid-column: 2;\n", " white-space: pre-wrap;\n", " word-break: break-all;\n", "}\n", "\n", ".xr-icon-database,\n", ".xr-icon-file-text2,\n", ".xr-no-icon {\n", " display: inline-block;\n", " vertical-align: middle;\n", " width: 1em;\n", " height: 1.5em !important;\n", " stroke-width: 0;\n", " stroke: currentColor;\n", " fill: currentColor;\n", "}\n", "</style><pre class='xr-text-repr-fallback'><xarray.Dataset>\n", "Dimensions: ()\n", "Data variables:\n", " b0 object 0.440+/-0.009\n", " by0 object 0.202+/-0.015\n", " alpha object 10.30+/-0.05</pre><div class='xr-wrap' style='display:none'><div class='xr-header'><div class='xr-obj-type'>xarray.Dataset</div></div><ul class='xr-sections'><li class='xr-section-item'><input id='section-65c281c6-2c36-4496-9e23-a700b90819ce' class='xr-section-summary-in' type='checkbox' disabled ><label for='section-65c281c6-2c36-4496-9e23-a700b90819ce' class='xr-section-summary' title='Expand/collapse section'>Dimensions:</label><div class='xr-section-inline-details'></div><div class='xr-section-details'></div></li><li class='xr-section-item'><input id='section-8318d64d-dd05-491e-a6f6-07f38de2619e' class='xr-section-summary-in' type='checkbox' disabled ><label for='section-8318d64d-dd05-491e-a6f6-07f38de2619e' class='xr-section-summary' title='Expand/collapse section'>Coordinates: <span>(0)</span></label><div class='xr-section-inline-details'></div><div class='xr-section-details'><ul class='xr-var-list'></ul></div></li><li class='xr-section-item'><input id='section-ab018521-2dc5-4f64-a759-8c951922b69c' class='xr-section-summary-in' type='checkbox' checked><label for='section-ab018521-2dc5-4f64-a759-8c951922b69c' class='xr-section-summary' >Data variables: <span>(3)</span></label><div class='xr-section-inline-details'></div><div class='xr-section-details'><ul class='xr-var-list'><li class='xr-var-item'><div class='xr-var-name'><span>b0</span></div><div class='xr-var-dims'>()</div><div class='xr-var-dtype'>object</div><div class='xr-var-preview xr-preview'>0.440+/-0.009</div><input id='attrs-909853d3-0cbb-4b52-a8da-d71ba069a7fb' class='xr-var-attrs-in' type='checkbox' disabled><label for='attrs-909853d3-0cbb-4b52-a8da-d71ba069a7fb' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-4546b020-7e06-402a-9ac2-83db067732a4' class='xr-var-data-in' type='checkbox'><label for='data-4546b020-7e06-402a-9ac2-83db067732a4' title='Show/Hide data repr'><svg class='icon xr-icon-database'><use xlink:href='#icon-database'></use></svg></label><div class='xr-var-attrs'><dl class='xr-attrs'></dl></div><div class='xr-var-data'><pre>array(0.44046053036292887+/-0.008743037526448819, dtype=object)</pre></div></li><li class='xr-var-item'><div class='xr-var-name'><span>by0</span></div><div class='xr-var-dims'>()</div><div class='xr-var-dtype'>object</div><div class='xr-var-preview xr-preview'>0.202+/-0.015</div><input id='attrs-90c38744-62fc-47b9-bb03-795b31d99a9a' class='xr-var-attrs-in' type='checkbox' disabled><label for='attrs-90c38744-62fc-47b9-bb03-795b31d99a9a' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-f6d9ce67-e2b6-4721-8ddc-5840cd33eec8' class='xr-var-data-in' type='checkbox'><label for='data-f6d9ce67-e2b6-4721-8ddc-5840cd33eec8' title='Show/Hide data repr'><svg class='icon xr-icon-database'><use xlink:href='#icon-database'></use></svg></label><div class='xr-var-attrs'><dl class='xr-attrs'></dl></div><div class='xr-var-data'><pre>array(0.2022134511996896+/-0.014697906285531557, dtype=object)</pre></div></li><li class='xr-var-item'><div class='xr-var-name'><span>alpha</span></div><div class='xr-var-dims'>()</div><div class='xr-var-dtype'>object</div><div class='xr-var-preview xr-preview'>10.30+/-0.05</div><input id='attrs-3524648c-1496-4088-bbe9-765d7f05351b' class='xr-var-attrs-in' type='checkbox' disabled><label for='attrs-3524648c-1496-4088-bbe9-765d7f05351b' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-0666538c-5594-4f58-af8b-f4d05c7a1744' class='xr-var-data-in' type='checkbox'><label for='data-0666538c-5594-4f58-af8b-f4d05c7a1744' title='Show/Hide data repr'><svg class='icon xr-icon-database'><use xlink:href='#icon-database'></use></svg></label><div class='xr-var-attrs'><dl class='xr-attrs'></dl></div><div class='xr-var-data'><pre>array(10.299175733301144+/-0.05309408600386392, dtype=object)</pre></div></li></ul></div></li><li class=' ], "text/plain": [ "<xarray.Dataset>\n", "Dimensions: ()\n", "Data variables:\n", " b0 object 0.440+/-0.009\n", " by0 object 0.202+/-0.015\n", " alpha object 10.30+/-0.05" ] }, "execution_count": 134, "metadata": {}, "output_type": "execute_result" } ], "source": [ "fitAnalyser_quadratic.get_fit_full_result(fitResult_quadratic) * 1e4" ] }, { "cell_type": "code", "execution_count": 135, "metadata": {}, "outputs": [], "source": [ "val_y = fitAnalyser_quadratic.get_fit_value(fitResult_quadratic) * 1e4\n", "std_y = fitAnalyser_quadratic.get_fit_std(fitResult_quadratic) * 1e4\n", "res_y = fitAnalyser_quadratic.get_fit_full_result(fitResult_quadratic) * 1e4" ] }, { "attachments": {}, "cell_type": "markdown", "metadata": {}, "source": [ "# Z-comp coil" ] }, { "cell_type": "code", "execution_count": 136, "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAkMAAAG4CAYAAABGhOPcAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjUuMSwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/YYfK9AAAACXBIWXMAAA9hAAAPYQGoP6dpAAAx9UlEQVR4nO3dfXSU5Z3/8c+YwAwhhhDCQ0ICuiBQIOjK8pDuRsKCGG0RGEDXI0sBW08Fa2LbVVAepBrAs6dbcjjQLUsxbfFhKQ6iLEIhBpoV5WGRcwI1RUBKiEEggSQCE5Jw/f6YX8YMkwQCMxMm9/t1zhzM9/7Oleu+DmU+vWfua2zGGCMAAACLuqO1JwAAANCaCEMAAMDSCEMAAMDSCEMAAMDSCEMAAMDSCEMAAMDSCEMAAMDSCEMAAMDSCEMAAMDSWhyGqqqq9MILL2jcuHHq2rWrbDabXnnllUZ7Dxw4oLFjxyo6OlqxsbFyOp06fvx4o70rVqzQgAEDZLfbdffdd2vx4sWqqanx6ztz5oxmzJih+Ph4RUVFKTU1VXl5eY2OuWPHDqWmpioqKkrx8fGaMWOGzpw509JTBgAAbViLw1BZWZlWr16t6upqTZw4scm+oqIipaen68qVK1q/fr3Wrl2rI0eOKC0tTWfPnvXpzc7OVmZmppxOp7Zt26bZs2dryZIlmjNnjk9fdXW1xowZo7y8POXk5GjTpk3q3r27MjIytGvXLp/eXbt26eGHH1b37t21adMm5eTkaMeOHRozZoyqq6tbetoAAKCtMi109epVc/XqVWOMMWfPnjWSzKJFi/z6pk6dauLj401FRYW3duLECdOuXTvzwgsveGvnzp0zDofDPP300z7Pz87ONjabzRw+fNhbW7lypZFkdu/e7a3V1NSYgQMHmuHDh/s8f9iwYWbgwIGmpqbGW/v444+NJLNq1aqWnjYAAGijWnxlyGazyWazNdtTW1urzZs3a/LkyYqJifHWe/furdGjR2vjxo3e2tatW+V2uzVz5kyfMWbOnCljjN577z1vbePGjerfv79SU1O9tcjISE2bNk179+5VSUmJJKmkpET79u3Tv/7rvyoyMtLb+93vflf9+vXz+f0AAMDaIq/f0nLHjh3T5cuXNWTIEL9jQ4YM0fbt2+V2u+VwOHTo0CFJUkpKik9fQkKC4uPjvccl6dChQ0pLS2t0TEk6fPiwevbs6X1OU7//448/vuFzuXr1qr766ivdeeed1w2BAADg9mCMUVVVlRITE3XHHc1f+wlKGCorK5MkxcXF+R2Li4uTMUbnz59XQkKCysrKZLfb1bFjx0Z768eqH7epMRv+3uv9/oZjXqu6utrnM0UlJSUaOHBgk/0AAOD2VVxcrKSkpGZ7ghKG6jV3JaXhsRvtC1Rvc2MsXbpUixcv9qsXFxf7vOUHAABuX5WVlUpOTtadd9553d6ghKEuXbpIUqNXYMrLy2Wz2RQbG+vtdbvdunTpkqKiovx6hw4d6jNuU2NK314Jut7vb+yKUb158+bppz/9qffn+sWMiYkhDAEAEGZu5CMuQdl0sU+fPurQoYMKCwv9jhUWFqpv375yOBySvv2s0LW9p0+f1rlz5zR48GBvLSUlpckxJXl76/9sqrfhmNey2+3e4EMAAgCg7QtKGIqMjNT48ePlcrlUVVXlrZ88eVL5+flyOp3eWkZGhhwOh3Jzc33GyM3Nlc1m89nLaNKkSSoqKtKePXu8tdraWq1bt04jRoxQYmKiJKlnz54aPny41q1bp7q6Om/vp59+qr/+9a8+vx8AAFibzRhjWvqkDz/8UBcvXlRVVZVmzZqlqVOn6rHHHpMkPfLII4qKilJRUZGGDRum+++/X3PnzpXb7dbChQtVXl6ugwcPqmvXrt7xsrOztWDBAs2bN0/jxo3Tvn37NH/+fE2fPl2rV6/29lVXV2vo0KGqrKzUsmXL1K1bN61atUoffPCBduzYoVGjRnl7d+7cqQcffFDjx4/X7NmzdebMGc2dO1edOnXS/v37Zbfbb+hcKysr1alTJ1VUVHCVCACAMNGi1++b2Zyod+/eRlKjjy+//NLbt3//fjNmzBgTFRVlYmJizMSJE83Ro0cbHTMnJ8f069fPtG/f3vTq1cssWrTIXLlyxa/v9OnTZvr06SYuLs44HA4zcuRIs3379kbH/NOf/mRGjhxpHA6HiYuLM9OnTzdff/11i861oqLCSPLZPBIAANzeWvL6fVNXhqyEK0MAAISflrx+8631AADA0ghDAADA0ghDAADA0ghDAADA0oL6dRwAAACNqauTCgqk0lIpIUFKS5MiIlpnLoQhAAAQUi6XlJkpnTr1bS0pScrJkVpjX2TeJgMAACHjcklTpvgGIUkqKfHUXa7Qz4kwBAAAQqKuznNFqLEdDutrWVmevlAiDAEAgJAoKPC/ItSQMVJxsacvlAhDAAAgJEpLA9sXKIQhAAAQEgkJge0LFMIQAAAIibQ0z11jNlvjx202KTnZ0xdKhCEAABASERGe2+cl/0BU//Py5aHfb4gwBAAAQsbplDZskHr29K0nJXnqrbHPEJsuAgCAkHI6pQkT2IEaAABYWESElJ7e2rPw4G0yAABgaYQhAABgaYQhAABgaYQhAABgaYQhAABgaYQhAABgaYQhAABgaYQhAABgaYQhAABgaYQhAABgaYQhAABgaYQhAABgaYQhAABgaYQhAABgaYQhAABgaYQhAABgaYQhAABgaYQhAABgaYQhAABgaYQhAABgaYQhAABgaYQhAABgaYQhAABgaYQhAABgaYQhAABgaYQhAABgaYQhAABgaYQhAABgaYQhAABgaYQhAABgaYQhAABgaYQhAABgaYQhAABgaYQhAABgaYQhAABgaYQhAABgaYQhAABgaZGtPQEAAKyork4qKJBKS6WEBCktTYqIaO1ZWRNhCACAEHO5pMxM6dSpb2tJSVJOjuR0tt68rIq3yQAACCGXS5oyxTcISVJJiafucrXOvKyMMAQAQIjU1XmuCBnjf6y+lpXl6UPoEIYAAAiRggL/K0INGSMVF3v6EDqEIQAAQqS0NLB9CAzCEAAAIZKQENg+BAZhCACAEElL89w1ZrM1ftxmk5KTPX0IHcIQAAAhEhHhuX1e8g9E9T8vX85+Q6FGGAIAIIScTmnDBqlnT996UpKnzj5DoRfUMPTZZ59p4sSJSkxMVFRUlAYMGKBf/OIXunTpkk/fgQMHNHbsWEVHRys2NlZOp1PHjx9vdMwVK1ZowIABstvtuvvuu7V48WLV1NT49Z05c0YzZsxQfHy8oqKilJqaqry8vKCcJwAALeF0SidOSPn50ltvef788kuCUGsJWhj6y1/+ou9+97s6ceKEli9frs2bN+tf/uVf9Itf/EJPPPGEt6+oqEjp6em6cuWK1q9fr7Vr1+rIkSNKS0vT2bNnfcbMzs5WZmamnE6ntm3bptmzZ2vJkiWaM2eOT191dbXGjBmjvLw85eTkaNOmTerevbsyMjK0a9euYJ0yAAA3LCJCSk+XnnjC8ydvjbUiEyQvv/yykWSOHj3qU3/66aeNJFNeXm6MMWbq1KkmPj7eVFRUeHtOnDhh2rVrZ1544QVv7dy5c8bhcJinn37aZ7zs7Gxjs9nM4cOHvbWVK1caSWb37t3eWk1NjRk4cKAZPnx4i86joqLCSPKZHwAAuL215PU7aFeG2rVrJ0nq1KmTTz02NlZ33HGH2rdvr9raWm3evFmTJ09WTEyMt6d3794aPXq0Nm7c6K1t3bpVbrdbM2fO9Blv5syZMsbovffe89Y2btyo/v37KzU11VuLjIzUtGnTtHfvXpWUlATyVAEAQBgLWhj6wQ9+oNjYWD3zzDM6fvy4qqqqtHnzZv3mN7/RnDlz1LFjRx07dkyXL1/WkCFD/J4/ZMgQHT16VG63W5J06NAhSVJKSopPX0JCguLj473H63ubGlOSDh8+HLDzBAAA4S1o31p/11136ZNPPtGkSZPUp08fb/25557T8uXLJUllZWWSpLi4OL/nx8XFyRij8+fPKyEhQWVlZbLb7erYsWOjvfVj1Y/b1Jg "text/plain": [ "<Figure size 640x480 with 1 Axes>" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "fz = [751.71168, 1223.54627, 1739.13344, 2364.93284, 2798.24971, 3178.49790, 4275.39905, 5352.17283, 6637.80418, 8288.35264, 9573.59333]\n", "dfz = [1.32859, 0.35554, 0.48471, 0.69762, 0.36873, 0.29413, 0.20667, 0.20818, 0.21978, 0.20285, 0.18495]\n", "z_offset_current = [0.0, 0.03, 0.06, 0.095, 0.119, 0.140, 0.2, 0.259, 0.329, 0.419, 0.489]\n", "\n", "data = xr.DataArray(\n", " data=fz,\n", " dims=['x'],\n", " coords={\n", " 'x': z_offset_current\n", " } \n", ")\n", "\n", "data_std = xr.DataArray(\n", " data=dfz,\n", " dims=['x'],\n", " coords={\n", " 'x': z_offset_current\n", " } \n", ")\n", "\n", "data.plot.errorbar(fmt='ob', yerr=data_std)\n", "\n", "plt.show()" ] }, { "cell_type": "code", "execution_count": 137, "metadata": {}, "outputs": [], "source": [ "def magnetic_field_func(x, b0=0, by0=0, alpha=1):\n", " return 1 / (1e3 * 6.626e-34) * (9.273e-24 * 1.24) * np.sqrt( (b0**2 - by0**2) + (alpha * x + by0)**2 )\n", "\n", "data_quadratic = data\n", "\n", "fitModel_quadratic = NewFitModel(magnetic_field_func)\n", "fitAnalyser_quadratic = FitAnalyser(fitModel_quadratic, fitDim=1)\n", "params_quadratic = fitAnalyser_quadratic.fitModel.make_params()\n", "params_quadratic.add(name=\"b0\", value= 0.3, max=np.inf, min=-np.inf, vary=True)\n", "params_quadratic.add(name=\"by0\", value= 0.07, max=np.inf, min=-np.inf, vary=True)\n", "params_quadratic.add(name=\"alpha\", value= 100, max=np.inf, min=-np.inf, vary=True)\n", "fitResult_quadratic = fitAnalyser_quadratic.fit(data_quadratic, params_quadratic).load()\n", "\n", "fitCurve_quadratic = fitAnalyser_quadratic.eval(fitResult_quadratic, x=np.linspace(0, 0.6, 100), dask=\"parallelized\").load()" ] }, { "cell_type": "code", "execution_count": 138, "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAm8AAAHPCAYAAAAFwj37AAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjUuMSwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/YYfK9AAAACXBIWXMAAA9hAAAPYQGoP6dpAACAEElEQVR4nO3dd3yN5//H8dfJOhlCJDFCFDVrk6KoXWoWKS2qVidqtL41vnaH0Wp/WkWnakvVqFFqK6o1K6i9Z2olQUJk378/8k3a0wQ5JzlOxvv5eOShue7rvs7n/jT4uO77vi6TYRgGIiIiIpIjODk6ABERERHJOBVvIiIiIjmIijcRERGRHETFm4iIiEgOouJNREREJAdR8SYiIiKSg6h4ExEREclBVLyJiIiI5CAujg4gr0lKSuKvv/7C29sbk8nk6HBEREQkmzAMg6ioKIoVK4aT0z3m1wwHi4yMNN58802jRYsWhr+/vwEY48aNs+iTkJBgfPDBB8aTTz5pFC9e3PDw8DAqVqxoDB8+3Lh+/Xq643788cdGhQoVDDc3N6NUqVLG+PHjjbi4uDT9rly5YvTq1cvw8/MzPDw8jMcee8zYsGFDumOuX7/eeOyxxwwPDw/Dz8/P6NWrl3HlyhWrrvfChQsGoC996Utf+tKXvvSV7teFCxfuWUs4fOYtPDyczz//nOrVq9OxY0e+/PLLNH3u3LnD+PHj6datGy+++CL+/v6EhITwzjvvsGLFCv744w88PDxS+7/77ruMGTOGESNG0LJlS3bv3s3o0aMJDQ3l888/T+0XGxtL8+bNuXHjBh999BGFCxdmxowZtGrVig0bNtC4cePUvlu2bKF169a0bduW5cuXc/XqVYYPH07z5s35448/MJvNGbpeb29vAC5cuED+/PltTdtdxcfHs27dOlq2bImrq2uWj5+bKXe2U+5sp9zZTrmzjfJmO3vnLjIykhIlSqTWCnfj8OKtZMmSXL9+HZPJRFhYWLrFm4eHB2fOnMHPzy+1rUmTJjz00EN06dKFH3/8kR49egDJxeA777zDSy+9xMSJE1P7xsfHM3r0aIYMGUKlSpUA+Oqrrzh48CDbtm2jXr16ADRt2pTq1aszbNgwdu7cmfp5b775JuXLl2fx4sW4uCSnrXTp0jRo0IDZs2fTr1+/DF1vyq3S/Pnz26148/T0JH/+/PpNaSXlznbKne2UO9spd7ZR3mz3oHJ3v8eqHP7Cgslkum+Qzs7OFoVbijp16gDJs1gp1qxZQ0xMDH369LHo26dPHwzDYNmyZaltS5cupUKFCqmFG4CLiws9evRg165dhIaGAhAaGsru3bt5/vnnUws3gPr161O+fHmWLl2a8QsWERERyQSHz7xlxi+//AJA5cqVU9sOHjwIQNWqVS36BgQE4O/vn3o8pW/Dhg3TjFutWjUADh06RPHixVPPSWn/d9/ff//9rjHGxsYSGxub+n1kZCSQXL3Hx8ff+wJtkDKmPcbO7ZQ72yl3tlPubKfc2UZ5s529c5fRcXNs8RYaGsqIESN49NFHadeuXWp7eHg4ZrMZLy+vNOf4+voSHh5u0dfX1zfdfinH//nr3fr+c8x/mzRpEhMmTEjTvm7dOjw9Pe96XmatX7/ebmPndsqd7ZQ72yl3tlPubKO82c5euYuOjs5QvxxZvEVERNCmTRsMw2DBggVpXqe9123Yfx/Lir73GmPkyJG88cYbqd+nPIzYsmVLuz3ztn79elq0aKFnGayk3NlOubOdcmc75c42ypvt7J27lLtz95Pjirfr16/TokULQkND+eWXX3j44Yctjvv5+RETE0N0dHSama2IiAiCgoIs+qY3axYREQH8PdOW8rzd3fqmNyOXwmw2p/smqqurq11/09h7/NxMubOdcmc75c52yp1tlDfb2St3GR3T4S8sWOP69es88cQTnDlzhvXr16f7DFrKs24HDhywaL98+TJhYWFUqVLFou+/+/3z3JS+Kb/ere8/xxQRERGxpxxTvKUUbqdPn2bdunXUrFkz3X6tWrXC3d2dOXPmWLTPmTMHk8lEx44dU9s6derE0aNHLZYESUhIYO7cudStW5dixYoBULx4cerUqcPcuXNJTExM7btjxw6OHTtGcHBw1l2oiIiIyD1ki9umq1ev5vbt20RFRQFw+PBhFi9eDECbNm0wmUw8+eST7N27l2nTppGQkMCOHTtSzy9UqBBlypQBkm91jh49mjFjxuDr65u6SO/48eN58cUXU9d4A+jbty8zZsygS5cuTJ48mcKFCzNz5kyOHTvGhg0bLGKcMmUKLVq0oEuXLvTv35+rV68yYsQIqlSpkmZZEhERERF7yRbFW79+/Th37lzq94sWLWLRokUAnDlzBoDdu3cDMHjw4DTn9+rVy2KmbdSoUXh7ezNjxgymTp1K0aJFGTFiBKNGjbI4z2w2s3HjRoYNG8bAgQOJjo6mRo0arF692mJ3BUhe6HfVqlWMHTuW9u3b4+npSbt27Xj//fczvLuCiIiISGZli+Lt7Nmz9+1jGIZVYw4aNIhBgwbdt1+RIkX45ptvMjRmixYtaNGihVVxiIiIiGSlHPPMm4iIiIioeBMRERHJUVS8iYiIiOQgKt5EREREMuh2NtgSVsWbiIiIyH0kJRnM2nKat/Y6c/LqLYfGki3eNhURERHJrm7eiWfown1sOHIVMLH64BUeKV7QYfGoeBMRERG5i0N/3aTf3BDOR0Tj5uJE8EPxDGxWxqEx6bapiIiISDoW77lI8MxtnI+IJrCgBwterEO9ItatO2sPmnkTERER+YeY+EQmrDjM/F3nAWhaoRD/92wNvFxNnN/v4OBQ8SYiIiKS6kJENP3nhXAg9CYmEwxpXp6Bzcri5GQiPj4bvGqKijcRERERADYfu8qQBfu4ER2Pj6crH3WtSePyhRwdVhoq3kRERCRPS0oy+PiXE3y08QSGAdUCCzDzuVoEFvR0dGjpUvEmIiIiedb123EMWbCPLcevAdC97kOMa18Js4uzgyO7OxVvIiIikif9efEG/eaGEHrjDmYXJ97tVJXOQYGODuu+VLyJiIhInmIYBvN3XWD8T4eIS0yipJ8ns54LolKx/I4OLUNUvImIiEieEROfyJhlB1m05yIATzxShA+eqU4BD1cHR5ZxKt5EREQkTzgXfptX54Zw5FIkTiZ488mKvNLoYZycTI4OzSoq3kRERCTX23D4Cq8v3EdUTAJ+Xm5M71aT+mX9HR2WTVS8iYiISK6VmGTwf+uP88mmkwDUesiHmc8FUbSAu4Mjs52KNxEREcmVwm/FMviHffx2MgyA3vVL8d82j+DmkrO3dlfxJiIiIrnO3vPX6T8vhEs3Y/BwdWby01XpUKO4o8PKEireREREJNcwDIO5O87x1srDxCcaPFzIi097BFG+iLejQ8syKt5EREQkV4iOS+C/Sw6wbN9fALSpWpQpT1fD2z3nLAOSESreREREJMc7fe0Wr87dw/Ert3B2MjGydUVeeLw0JlPOWgYkI1S8iYiISI625uAl/rPoT27FJlDI28wn3WpS92E/R4dlNyreREREJEdKSEzivbXH+PzX0wDUKeXLJ91rUjh/zl0GJCNUvImIiEiOczUqhoHf72XnmQgAXmpYmmGtKuLqnLOXAckIFW8iIiKSo+w+G8GAeSFcjYrFy82Z97tUp03VAEeH9cCoeBMREZEcwTAMvvrtDJNWHyUxyaBc4XzM6hFE2cL5HB3aA6XiTURERLK9W7EJDF/8Jz8fuATAU9WLMSm4Kl7mvFfK5L0rFhERkRzlxJUoXp27h1PXbuPqbGJ020r0rFcyVy4DkhEq3kRERCTbWrH/L4b/+CfRcYkUze/OjOdqEVSyoKPDcigVbyIiIpLtxCUkMWn1Eb7+/SwA9cv48XG3mvjnMzs2sGxAxZuIiIhkK5dvxjDg+xD2nLsOQP8mZRjasgLOTnnzNum/qXgTERGRbGPbqTAGzd9L2K04vN1d+PCZGrSoVMTRYWUrKt5ERETE4QzD4NMtp3l/7VGSDKhY1JtPewRRyt/L0aFlOyreRERExKEiY+IZunA/6w9 "text/plain": [ "<Figure size 640x480 with 1 Axes>" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "data_linear = data[3:]\n", "\n", "fitAnalyser_linear = FitAnalyser('Linear', fitDim=1)\n", "params_linear = fitAnalyser_linear.guess(data_linear, dask=\"parallelized\")\n", "fitResult_linear = fitAnalyser_linear.fit(data_linear, params_linear).load()\n", "\n", "fitCurve_linear = fitAnalyser_linear.eval(fitResult_linear, x=np.linspace(0, 0.5, 100), dask=\"parallelized\").load()\n", "\n", "fig = plt.figure()\n", "ax = fig.gca()\n", "\n", "data.plot.errorbar(ax=ax, fmt='ob', yerr=data_std)\n", "# fitCurve_linear.plot.errorbar(ax=ax)\n", "fitCurve_quadratic.plot.errorbar(ax=ax)\n", "\n", "plt.xlabel('Y Offset Coil Current (A)', fontsize=12)\n", "plt.ylabel('Resonance Frequency (kHz)', fontsize=12)\n", "plt.xticks(fontsize=12)\n", "plt.yticks(fontsize=12)\n", "# plt.legend(fontsize=12)\n", "#plt.xlim(-0.01, 0.04)\n", "# plt.ylim(0, 10000)\n", "plt.grid(visible=1)\n", "plt.tight_layout()\n", "plt.show()" ] }, { "cell_type": "code", "execution_count": 139, "metadata": {}, "outputs": [ { "data": { "text/html": [ "<div><svg style=\"position: absolute; width: 0; height: 0; overflow: hidden\">\n", "<defs>\n", "<symbol id=\"icon-database\" viewBox=\"0 0 32 32\">\n", "<path d=\"M16 0c-8.837 0-16 2.239-16 5v4c0 2.761 7.163 5 16 5s16-2.239 16-5v-4c0-2.761-7.163-5-16-5z\"></path>\n", "<path d=\"M16 17c-8.837 0-16-2.239-16-5v6c0 2.761 7.163 5 16 5s16-2.239 16-5v-6c0 2.761-7.163 5-16 5z\"></path>\n", "<path d=\"M16 26c-8.837 0-16-2.239-16-5v6c0 2.761 7.163 5 16 5s16-2.239 16-5v-6c0 2.761-7.163 5-16 5z\"></path>\n", "</symbol>\n", "<symbol id=\"icon-file-text2\" viewBox=\"0 0 32 32\">\n", "<path d=\"M28.681 7.159c-0.694-0.947-1.662-2.053-2.724-3.116s-2.169-2.030-3.116-2.724c-1.612-1.182-2.393-1.319-2.841-1.319h-15.5c-1.378 0-2.5 1.121-2.5 2.5v27c0 1.378 1.122 2.5 2.5 2.5h23c1.378 0 2.5-1.122 2.5-2.5v-19.5c0-0.448-0.137-1.23-1.319-2.841zM24.543 5.457c0.959 0.959 1.712 1.825 2.268 2.543h-4.811v-4.811c0.718 0.556 1.584 1.309 2.543 2.268zM28 29.5c0 0.271-0.229 0.5-0.5 0.5h-23c-0.271 0-0.5-0.229-0.5-0.5v-27c0-0.271 0.229-0.5 0.5-0.5 0 0 15.499-0 15.5 0v7c0 0.552 0.448 1 1 1h7v19.5z\"></path>\n", "<path d=\"M23 26h-14c-0.552 0-1-0.448-1-1s0.448-1 1-1h14c0.552 0 1 0.448 1 1s-0.448 1-1 1z\"></path>\n", "<path d=\"M23 22h-14c-0.552 0-1-0.448-1-1s0.448-1 1-1h14c0.552 0 1 0.448 1 1s-0.448 1-1 1z\"></path>\n", "<path d=\"M23 18h-14c-0.552 0-1-0.448-1-1s0.448-1 1-1h14c0.552 0 1 0.448 1 1s-0.448 1-1 1z\"></path>\n", "</symbol>\n", "</defs>\n", "</svg>\n", "<style>/* CSS stylesheet for displaying xarray objects in jupyterlab.\n", " *\n", " */\n", "\n", ":root {\n", " --xr-font-color0: var(--jp-content-font-color0, rgba(0, 0, 0, 1));\n", " --xr-font-color2: var(--jp-content-font-color2, rgba(0, 0, 0, 0.54));\n", " --xr-font-color3: var(--jp-content-font-color3, rgba(0, 0, 0, 0.38));\n", " --xr-border-color: var(--jp-border-color2, #e0e0e0);\n", " --xr-disabled-color: var(--jp-layout-color3, #bdbdbd);\n", " --xr-background-color: var(--jp-layout-color0, white);\n", " --xr-background-color-row-even: var(--jp-layout-color1, white);\n", " --xr-background-color-row-odd: var(--jp-layout-color2, #eeeeee);\n", "}\n", "\n", "html[theme=dark],\n", "body[data-theme=dark],\n", "body.vscode-dark {\n", " --xr-font-color0: rgba(255, 255, 255, 1);\n", " --xr-font-color2: rgba(255, 255, 255, 0.54);\n", " --xr-font-color3: rgba(255, 255, 255, 0.38);\n", " --xr-border-color: #1F1F1F;\n", " --xr-disabled-color: #515151;\n", " --xr-background-color: #111111;\n", " --xr-background-color-row-even: #111111;\n", " --xr-background-color-row-odd: #313131;\n", "}\n", "\n", ".xr-wrap {\n", " display: block !important;\n", " min-width: 300px;\n", " max-width: 700px;\n", "}\n", "\n", ".xr-text-repr-fallback {\n", " /* fallback to plain text repr when CSS is not injected (untrusted notebook) */\n", " display: none;\n", "}\n", "\n", ".xr-header {\n", " padding-top: 6px;\n", " padding-bottom: 6px;\n", " margin-bottom: 4px;\n", " border-bottom: solid 1px var(--xr-border-color);\n", "}\n", "\n", ".xr-header > div,\n", ".xr-header > ul {\n", " display: inline;\n", " margin-top: 0;\n", " margin-bottom: 0;\n", "}\n", "\n", ".xr-obj-type,\n", ".xr-array-name {\n", " margin-left: 2px;\n", " margin-right: 10px;\n", "}\n", "\n", ".xr-obj-type {\n", " color: var(--xr-font-color2);\n", "}\n", "\n", ".xr-sections {\n", " padding-left: 0 !important;\n", " display: grid;\n", " grid-template-columns: 150px auto auto 1fr 20px 20px;\n", "}\n", "\n", ".xr-section-item {\n", " display: contents;\n", "}\n", "\n", ".xr-section-item input {\n", " display: none;\n", "}\n", "\n", ".xr-section-item input + label {\n", " color: var(--xr-disabled-color);\n", "}\n", "\n", ".xr-section-item input:enabled + label {\n", " cursor: pointer;\n", " color: var(--xr-font-color2);\n", "}\n", "\n", ".xr-section-item input:enabled + label:hover {\n", " color: var(--xr-font-color0);\n", "}\n", "\n", ".xr-section-summary {\n", " grid-column: 1;\n", " color: var(--xr-font-color2);\n", " font-weight: 500;\n", "}\n", "\n", ".xr-section-summary > span {\n", " display: inline-block;\n", " padding-left: 0.5em;\n", "}\n", "\n", ".xr-section-summary-in:disabled + label {\n", " color: var(--xr-font-color2);\n", "}\n", "\n", ".xr-section-summary-in + label:before {\n", " display: inline-block;\n", " content: '►';\n", " font-size: 11px;\n", " width: 15px;\n", " text-align: center;\n", "}\n", "\n", ".xr-section-summary-in:disabled + label:before {\n", " color: var(--xr-disabled-color);\n", "}\n", "\n", ".xr-section-summary-in:checked + label:before {\n", " content: '▼';\n", "}\n", "\n", ".xr-section-summary-in:checked + label > span {\n", " display: none;\n", "}\n", "\n", ".xr-section-summary,\n", ".xr-section-inline-details {\n", " padding-top: 4px;\n", " padding-bottom: 4px;\n", "}\n", "\n", ".xr-section-inline-details {\n", " grid-column: 2 / -1;\n", "}\n", "\n", ".xr-section-details {\n", " display: none;\n", " grid-column: 1 / -1;\n", " margin-bottom: 5px;\n", "}\n", "\n", ".xr-section-summary-in:checked ~ .xr-section-details {\n", " display: contents;\n", "}\n", "\n", ".xr-array-wrap {\n", " grid-column: 1 / -1;\n", " display: grid;\n", " grid-template-columns: 20px auto;\n", "}\n", "\n", ".xr-array-wrap > label {\n", " grid-column: 1;\n", " vertical-align: top;\n", "}\n", "\n", ".xr-preview {\n", " color: var(--xr-font-color3);\n", "}\n", "\n", ".xr-array-preview,\n", ".xr-array-data {\n", " padding: 0 5px !important;\n", " grid-column: 2;\n", "}\n", "\n", ".xr-array-data,\n", ".xr-array-in:checked ~ .xr-array-preview {\n", " display: none;\n", "}\n", "\n", ".xr-array-in:checked ~ .xr-array-data,\n", ".xr-array-preview {\n", " display: inline-block;\n", "}\n", "\n", ".xr-dim-list {\n", " display: inline-block !important;\n", " list-style: none;\n", " padding: 0 !important;\n", " margin: 0;\n", "}\n", "\n", ".xr-dim-list li {\n", " display: inline-block;\n", " padding: 0;\n", " margin: 0;\n", "}\n", "\n", ".xr-dim-list:before {\n", " content: '(';\n", "}\n", "\n", ".xr-dim-list:after {\n", " content: ')';\n", "}\n", "\n", ".xr-dim-list li:not(:last-child):after {\n", " content: ',';\n", " padding-right: 5px;\n", "}\n", "\n", ".xr-has-index {\n", " font-weight: bold;\n", "}\n", "\n", ".xr-var-list,\n", ".xr-var-item {\n", " display: contents;\n", "}\n", "\n", ".xr-var-item > div,\n", ".xr-var-item label,\n", ".xr-var-item > .xr-var-name span {\n", " background-color: var(--xr-background-color-row-even);\n", " margin-bottom: 0;\n", "}\n", "\n", ".xr-var-item > .xr-var-name:hover span {\n", " padding-right: 5px;\n", "}\n", "\n", ".xr-var-list > li:nth-child(odd) > div,\n", ".xr-var-list > li:nth-child(odd) > label,\n", ".xr-var-list > li:nth-child(odd) > .xr-var-name span {\n", " background-color: var(--xr-background-color-row-odd);\n", "}\n", "\n", ".xr-var-name {\n", " grid-column: 1;\n", "}\n", "\n", ".xr-var-dims {\n", " grid-column: 2;\n", "}\n", "\n", ".xr-var-dtype {\n", " grid-column: 3;\n", " text-align: right;\n", " color: var(--xr-font-color2);\n", "}\n", "\n", ".xr-var-preview {\n", " grid-column: 4;\n", "}\n", "\n", ".xr-index-preview {\n", " grid-column: 2 / 5;\n", " color: var(--xr-font-color2);\n", "}\n", "\n", ".xr-var-name,\n", ".xr-var-dims,\n", ".xr-var-dtype,\n", ".xr-preview,\n", ".xr-attrs dt {\n", " white-space: nowrap;\n", " overflow: hidden;\n", " text-overflow: ellipsis;\n", " padding-right: 10px;\n", "}\n", "\n", ".xr-var-name:hover,\n", ".xr-var-dims:hover,\n", ".xr-var-dtype:hover,\n", ".xr-attrs dt:hover {\n", " overflow: visible;\n", " width: auto;\n", " z-index: 1;\n", "}\n", "\n", ".xr-var-attrs,\n", ".xr-var-data,\n", ".xr-index-data {\n", " display: none;\n", " background-color: var(--xr-background-color) !important;\n", " padding-bottom: 5px !important;\n", "}\n", "\n", ".xr-var-attrs-in:checked ~ .xr-var-attrs,\n", ".xr-var-data-in:checked ~ .xr-var-data,\n", ".xr-index-data-in:checked ~ .xr-index-data {\n", " display: block;\n", "}\n", "\n", ".xr-var-data > table {\n", " float: right;\n", "}\n", "\n", ".xr-var-name span,\n", ".xr-var-data,\n", ".xr-index-name div,\n", ".xr-index-data,\n", ".xr-attrs {\n", " padding-left: 25px !important;\n", "}\n", "\n", ".xr-attrs,\n", ".xr-var-attrs,\n", ".xr-var-data,\n", ".xr-index-data {\n", " grid-column: 1 / -1;\n", "}\n", "\n", "dl.xr-attrs {\n", " padding: 0;\n", " margin: 0;\n", " display: grid;\n", " grid-template-columns: 125px auto;\n", "}\n", "\n", ".xr-attrs dt,\n", ".xr-attrs dd {\n", " padding: 0;\n", " margin: 0;\n", " float: left;\n", " padding-right: 10px;\n", " width: auto;\n", "}\n", "\n", ".xr-attrs dt {\n", " font-weight: normal;\n", " grid-column: 1;\n", "}\n", "\n", ".xr-attrs dt:hover span {\n", " display: inline-block;\n", " background: var(--xr-background-color);\n", " padding-right: 10px;\n", "}\n", "\n", ".xr-attrs dd {\n", " grid-column: 2;\n", " white-space: pre-wrap;\n", " word-break: break-all;\n", "}\n", "\n", ".xr-icon-database,\n", ".xr-icon-file-text2,\n", ".xr-no-icon {\n", " display: inline-block;\n", " vertical-align: middle;\n", " width: 1em;\n", " height: 1.5em !important;\n", " stroke-width: 0;\n", " stroke: currentColor;\n", " fill: currentColor;\n", "}\n", "</style><pre class='xr-text-repr-fallback'><xarray.Dataset>\n", "Dimensions: ()\n", "Data variables:\n", " b0 object 0.4333+/-0.0007\n", " by0 object 0.3254+/-0.0008\n", " alpha object 10.6021+/-0.0022</pre><div class='xr-wrap' style='display:none'><div class='xr-header'><div class='xr-obj-type'>xarray.Dataset</div></div><ul class='xr-sections'><li class='xr-section-item'><input id='section-f77d168d-c7f1-4cbb-be7f-2cef6dfaf365' class='xr-section-summary-in' type='checkbox' disabled ><label for='section-f77d168d-c7f1-4cbb-be7f-2cef6dfaf365' class='xr-section-summary' title='Expand/collapse section'>Dimensions:</label><div class='xr-section-inline-details'></div><div class='xr-section-details'></div></li><li class='xr-section-item'><input id='section-57ddfb08-9489-4068-b143-f0306231c01c' class='xr-section-summary-in' type='checkbox' disabled ><label for='section-57ddfb08-9489-4068-b143-f0306231c01c' class='xr-section-summary' title='Expand/collapse section'>Coordinates: <span>(0)</span></label><div class='xr-section-inline-details'></div><div class='xr-section-details'><ul class='xr-var-list'></ul></div></li><li class='xr-section-item'><input id='section-cc75bd83-4893-4666-8e8c-e9d9cbec9fb9' class='xr-section-summary-in' type='checkbox' checked><label for='section-cc75bd83-4893-4666-8e8c-e9d9cbec9fb9' class='xr-section-summary' >Data variables: <span>(3)</span></label><div class='xr-section-inline-details'></div><div class='xr-section-details'><ul class='xr-var-list'><li class='xr-var-item'><div class='xr-var-name'><span>b0</span></div><div class='xr-var-dims'>()</div><div class='xr-var-dtype'>object</div><div class='xr-var-preview xr-preview'>0.4333+/-0.0007</div><input id='attrs-d69271e3-d2be-4615-b6e9-16971b1d9bee' class='xr-var-attrs-in' type='checkbox' disabled><label for='attrs-d69271e3-d2be-4615-b6e9-16971b1d9bee' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-9859a199-e6b6-4124-869a-40487a3bd516' class='xr-var-data-in' type='checkbox'><label for='data-9859a199-e6b6-4124-869a-40487a3bd516' title='Show/Hide data repr'><svg class='icon xr-icon-database'><use xlink:href='#icon-database'></use></svg></label><div class='xr-var-attrs'><dl class='xr-attrs'></dl></div><div class='xr-var-data'><pre>array(0.43327107997401876+/-0.0007013859839914253, dtype=object)</pre></div></li><li class='xr-var-item'><div class='xr-var-name'><span>by0</span></div><div class='xr-var-dims'>()</div><div class='xr-var-dtype'>object</div><div class='xr-var-preview xr-preview'>0.3254+/-0.0008</div><input id='attrs-001a69b2-d62f-4955-a882-ed50e4b925ea' class='xr-var-attrs-in' type='checkbox' disabled><label for='attrs-001a69b2-d62f-4955-a882-ed50e4b925ea' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-637572eb-5cc6-461d-bf12-6de7cd87d0f5' class='xr-var-data-in' type='checkbox'><label for='data-637572eb-5cc6-461d-bf12-6de7cd87d0f5' title='Show/Hide data repr'><svg class='icon xr-icon-database'><use xlink:href='#icon-database'></use></svg></label><div class='xr-var-attrs'><dl class='xr-attrs'></dl></div><div class='xr-var-data'><pre>array(0.32538643700542086+/-0.000794288942625618, dtype=object)</pre></div></li><li class='xr-var-item'><div class='xr-var-name'><span>alpha</span></div><div class='xr-var-dims'>()</div><div class='xr-var-dtype'>object</div><div class='xr-var-preview xr-preview'>10.6021+/-0.0022</div><input id='attrs-36ac86c4-9b53-49ad-9820-b0daabbf0894' class='xr-var-attrs-in' type='checkbox' disabled><label for='attrs-36ac86c4-9b53-49ad-9820-b0daabbf0894' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-03225bdf-14a6-42bc-9f5d-9f12aa4b7152' class='xr-var-data-in' type='checkbox'><label for='data-03225bdf-14a6-42bc-9f5d-9f12aa4b7152' title='Show/Hide data repr'><svg class='icon xr-icon-database'><use xlink:href='#icon-database'></use></svg></label><div class='xr-var-attrs'><dl class='xr-attrs'></dl></div><div class='xr-var-data'><pre>array(10.602080670806222+/-0.0022465278174072264, dtype=object)</pre></div></li></ul></div> ], "text/plain": [ "<xarray.Dataset>\n", "Dimensions: ()\n", "Data variables:\n", " b0 object 0.4333+/-0.0007\n", " by0 object 0.3254+/-0.0008\n", " alpha object 10.6021+/-0.0022" ] }, "execution_count": 139, "metadata": {}, "output_type": "execute_result" } ], "source": [ "fitAnalyser_quadratic.get_fit_full_result(fitResult_quadratic) * 1e4" ] }, { "cell_type": "code", "execution_count": 140, "metadata": {}, "outputs": [], "source": [ "val_z = fitAnalyser_quadratic.get_fit_value(fitResult_quadratic) * 1e4\n", "std_z = fitAnalyser_quadratic.get_fit_std(fitResult_quadratic) * 1e4\n", "res_z = fitAnalyser_quadratic.get_fit_full_result(fitResult_quadratic) * 1e4" ] }, { "cell_type": "code", "execution_count": null, "metadata": {}, "outputs": [], "source": [] }, { "cell_type": "code", "execution_count": 145, "metadata": {}, "outputs": [ { "data": { "text/plain": [ "0.416688520216393+/-0.009836437571434223" ] }, "execution_count": 145, "metadata": {}, "output_type": "execute_result" } ], "source": [ "umath.sqrt(res_x['by0'].item()**2 + res_y['by0'].item()**2 + res_z['by0'].item()**2)" ] }, { "cell_type": "code", "execution_count": 148, "metadata": {}, "outputs": [ { "data": { "text/plain": [ "0.43640594986928966+/-0.006152786477679057" ] }, "execution_count": 148, "metadata": {}, "output_type": "execute_result" } ], "source": [ "(res_x['b0'].item() + res_y['b0'].item() + res_z['b0'].item())/3" ] }, { "cell_type": "code", "execution_count": null, "metadata": {}, "outputs": [], "source": [] } ], "metadata": { "kernelspec": { "display_name": "base", "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" }, "orig_nbformat": 4 }, "nbformat": 4, "nbformat_minor": 2}
|
