gao/analyseScript
1
0
Fork 0
Code Issues Pull Requests repo.project_board Releases Wiki Activity
e6e20a8ae3
analyseScript/Analyser/Joschka/test_rotation.ipynb

1237 lines
2.3 MiB
Plaintext
Raw Normal View History

finished adding rotation
2023-09-14 14:13:49 +02:00
{
"cells": [
{
"cell_type": "code",
"execution_count": 1,
"metadata": {},
"outputs": [],
"source": [
"import lmfit\n",
"import xarray as xr\n",
"import pandas as pd\n",
"import numpy as np\n",
"import copy\n",
"\n",
"import glob\n",
"\n",
"import xrft\n",
"import finufft\n",
"\n",
"from uncertainties import ufloat\n",
"from uncertainties import unumpy as unp\n",
"from uncertainties import umath\n",
"\n",
"from datetime import datetime\n",
"\n",
"import matplotlib.pyplot as plt\n",
"\n",
"#test\n",
"plt.rcParams['font.size'] = 18\n",
"\n",
"from scipy.ndimage import gaussian_filter, rotate\n",
"import matplotlib as mpl\n",
"from scipy.interpolate import CubicSpline\n",
"from scipy.optimize import curve_fit\n",
"mpl.rc('xtick', labelsize=8)\n",
"mpl.rc('ytick', labelsize=8)\n",
"\n",
"import sys\n",
"sys.path.append(\"C:/Users/Jianshun Gao/VisualCodeProjects/analysescript\")\n",
"sys.path.append(\"C:/Users/Jianshun Gao/PycharmProjects/Foglight4/Data\")\n",
"\n",
"from DataContainer.ReadData import read_hdf5_file, read_hdf5_global, read_hdf5_run_time, read_csv_file\n",
"from Analyser.ImagingAnalyser import ImageAnalyser\n",
"from Analyser.FitAnalyser import FitAnalyser\n",
"from Analyser.FitAnalyser import ThomasFermi2dModel, DensityProfileBEC2dModel, Polylog22dModel\n",
"from Analyser.FFTAnalyser import fft, ifft, fft_nutou\n",
"from ToolFunction.ToolFunction import *\n",
"\n",
"import time\n",
"import logging as log\n",
"import scipy.constants as const\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",
"\n",
"import Analyser.FitAnalyser as dp"
]
},
{
"cell_type": "code",
"execution_count": 2,
"metadata": {},
"outputs": [
{
"name": "stdout",
"output_type": "stream",
"text": [
"C:/Users/Jianshun Gao/PycharmProjects/Foglight4/Data/josch/high_cond_1.nc\n"
]
}
],
"source": [
"path = \"C:/Users/Jianshun Gao/PycharmProjects/Foglight4/Data\"\n",
"file = 'high_cond_1'\n",
"file_path = f'{path}/josch/{file}.nc'\n",
"print(file_path)\n",
"with xr.open_dataarray(file_path) as data_x:\n",
" data_x\n",
"\n",
"data = data_x.to_numpy()"
]
},
{
"cell_type": "code",
"execution_count": 3,
"metadata": {},
"outputs": [],
"source": [
"is_debug = True\n",
"fitm = DensityProfileBEC2dModel(is_debug=is_debug, atom_n_conv=144/2)\n",
"\n",
"def calc_cen_bulk(thresh):\n",
" \"\"\"\n",
" returns array in shape of input, containing array with [Y_center,X_center] for each image\n",
" \"\"\"\n",
" shape = np.shape(thresh)\n",
" cen = np.zeros((shape[0], 2))\n",
" for i in range(0, shape[0]):\n",
" cen[i] = fitm.calc_cen_pix(thresh[i])\n",
" return cen"
]
},
{
"cell_type": "code",
"execution_count": 4,
"metadata": {},
"outputs": [],
"source": [
"shape = np.shape(data)\n",
"data_rot = np.empty(shape)\n",
"\n",
"for i in range(0,shape[0]):\n",
" data_rot[i] = rotate(data[i], 36, reshape=False)"
]
},
{
"cell_type": "code",
"execution_count": 5,
"metadata": {},
"outputs": [
{
"name": "stdout",
"output_type": "stream",
"text": [
"[[599.42377261 959.62015504]\n",
" [599.46599496 959.50629723]\n",
" [599.47858942 959.49118388]\n",
" [599.47505938 959.41330166]\n",
" [599.68149883 959.26697892]\n",
" [599.4627907 959.58372093]\n",
" [599.46412556 959.46188341]]\n"
]
}
],
"source": [
"# cropping\n",
"pr_data = data_rot\n",
"\n",
"cut_width_x = 200\n",
"cut_width_y = 150\n",
"thresh = fitm.calc_thresh(pr_data)\n",
"\n",
"center = calc_cen_bulk(thresh)\n",
"\n",
"print(center)\n",
"\n",
"\n",
"if len(shape) == 3:\n",
" cropOD = np.zeros((shape[0], cut_width_y, cut_width_x))\n",
" for i in range(0,shape[0]):\n",
" for j in range(0, shape[1]):\n",
" cropOD[i] = pr_data[i, round(center[i,0]-cut_width_y/2):round(center[i,0]+cut_width_y/2), round(center[i,1]-cut_width_x/2):round(center[i,1]+cut_width_x/2)]\n",
"\n",
"shape = np.shape(cropOD)"
]
},
{
"cell_type": "code",
"execution_count": 6,
"metadata": {},
"outputs": [
{
"data": {
"image/png": "iVBORw0KGgoAAAANSUhEUgAABYkAAADoCAYAAACq22t0AAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjcuMSwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/bCgiHAAAACXBIWXMAAA9hAAAPYQGoP6dpAAEAAElEQVR4nOz9XYxk2XmeCz5NZ5AdJWa0GElV0J0pdliqspQ+anoOPbIICJDmYGBY8I2MgTQYCNBcDAT9zIVuzo10qRsLBgwBhj2wbAMWMMeAB8LYMDAQPGMfA5IsAZYFNORuwyW5ynY0mUkxis6gGUlVdCuS5Fzs/fB7185qUaSqWdFHewGJzIzYe+21vvX9vN/PWvuFL3/5y19mbGMb29jGNraxjW1sYxvb2MY2trGNbWxjG9vYxvansr3veQ9gbGMb29jGNraxjW1sYxvb2MY2trGNbWxjG9vYxvb82hgkHtvYxja2sY1tbGMb29jGNraxjW1sYxvb2MY2tj/FbQwSj21sYxvb2MY2trGNbWxjG9vYxja2sY1tbGMb25/iNgaJxza2sY1tbGMb29jGNraxjW1sYxvb2MY2trGN7U9xG4PEYxvb2MY2trGNbWxjG9vYxja2sY1tbGMb29jG9qe4jUHisY1tbGMb29jGNraxjW1sYxvb2MY2trGNbWxj+1PcxiDx2MY2trGNbWxjG9vYxja2sY1tbGMb29jGNrax/SluY5B4bGMb29jGNraxjW1sYxvb2MY2trGNbWxjG9vY/hS3o+c9gK/WvvSlL/HpT3+a4+NjXnjhhec9nLGNjS9/+ctcX1/z8ssv8773PZ88yygXYzvENsrG2MZ2u41yMbaxPb2NsjG2sd1uo1yMbWxPb6NsjG1st9u7IRcHHyT+9Kc/zbd+67c+72GMbWy32qc+9SnOzs6ey7NHuRjbIbdRNsY2ttttlIuxje3pbZSNsY3tdhvlYmxje3obZWNsY7vdnqVcHHyQ+Pj4uP/rHwEzYAdM6Ib+GPD7HfAy8Ln+77eAbwNu+r8/2H837T/7IPDpvp+3gDmwoSXJUf/5J/vrZ8CLwH8B7gLbuH4W99A/41uA3+8/2/R9vdh/PwWu+2dv+v5e7Mc46a/Z9tfZPtR/75xv+r6P+7kc9zTxOV+IPq6Bj8b9j6PvaT/W/9z3ORvMY9PT1mfsgTv95y/219DP4bh/xk3/jH3/3Sv99a/0/18P6PJR4Hd6erzc/37SP+dDwGX/2Yv9+L6l//96QFuv+dxgnN/eP3fXz+PbKT75s30/zuOt/rdjv46+j4APA/+n4M1vfKtn/zPg/XQ8s6fW9y06uj2mo+fv9dfLt3fo5rLrrzum43P79fen+76/pe8bil83MaJZ//+s72/Xf/4hOhnY95/d9H1P6GTqC3Tr/O20awfFg34u/Y9j3J/tP/tC//z/oX+GdLjur90B/3s62VVvKGcf7D+jH8dnKR0hLXz2Z/uxKE/KoDxzHf/Ld/LWE0pPfJpuXY7759jPp3paX9PJ5b4f2xeA7+iv/Wx/zWcpOXauH6J4YUrJonSj/+xbgDeBU2qtvqV/zot93x8KWr7Yj1ce+VzM54jijw3wfz0Q2fgHFF+qtyb9j7pb2vw+3XzVczs6Hr+meMV7vhD30l8rr77S//1JOv5wHW1v9v25hq7xHbq1lKfe7H+/3F+Ttsk53dCurf086cf4wbj+0/3fH6T4w3m8GH292H8/7cfzYlyXbUfHO9LB8fjMz/X02/efp62hp6V21//p//80ZetyfLledyjb9rgfs3rHZ6hH9v1YfYbNdf9o/8wphSNc73nMz7H+Zzp6Ox5pmvoNSs9c9///zwciF3+botGcbtyfo9VDN3Tr92lqjnNqLtf9/9v+Opv8fNNfI1bb0dFCWtLf90m6dbM/+vv9fB735xh9lvo25XxCxzPqKXnXed3041b/K+dipy2FQTbxmdc7pw2lYz9I2akhvx9F37v+WWmHGTzH67Q76nLlURy8p2RI2/9tlA7/z5Rdla7KR9oIgiZiq09H/+LLD/XfP6a18Yn/tMmXlPyrP16Ovj7ZP+cVYM3h2Iz/Bx1tXu7/F0do08UhYtn/TNmMI7o5py79Pbo1kAev6ez4J/PpdPretVCXSX/XWlt7RK1f6h1lUb657Mf4JrU28oK6EQoHPKbk3Ll+mk4/5vMuKT3heI8pvKMtFY84b3Whtki6pB3NlrL1uehLbD/r75UXP0qH717s5yLvvtxfdxnzeKvv5ztodZN8+rm+j7Qn6hUoPCi+tj3ux5Xy6DiOKP90Qud3fK6nx5yynfbxh8CPHIhc/D+B/97/rb+gTKiXpC10/PYWHU3lG9uTuOaTlN5TZl6h6LSjMLP+gjhsRuks73Vd6L9T9yiX8uVjWp87ZSnxmvhaGydvnVK2RT2+pzDzhsIoyeP2fU3H/z57189b2ZpR/AClc2Y9/dQR8tZjCm+IY4xXXNPx6Gep9Uk7oD22H3GdPoq4SVyojdJHF2tKI22/c9W30mfU3iYep/9fOZQfptHPUfT1zcD/+YBk47uoeeszG79wLse0/rU2VRn6DlrfV99ZrCHuTTwhTt5Ta/tZOl7a02Ib+1ePOs4/S/m82hfX1nXOdknHI9Ctu3r93wOfoOIq+ulQWFBcph+d8mA/yot2TwyesTNlz/scs9d9kvIdxJU3dDbiP9DGP+RDffoJbazDeSQuuw5aigeJa8XBYlBjZ9/WX/9f+nvepOMd125C2R3n5Dz/dxQOVmbEmz7nD3nWWOrgg8RVxj8DXqAjyISOeGfUFE7pCPml/rspnfJa0B29/CXgpf76KfAIOOnvvUM5GS/0famovolO0Cb9NWtgSQscBGkn/T2PgHvA5/tr/E7APKNTcm/3YzuhlMeHKUdMBSoonPbj0fB+E3BFFySUHu+L5yxiXvTf3e/7d84XPV2+1D9bAyQTvtXP96bvTwdwAnycCmBoFL7U962hOutp9v5+nl/sv/tz/fcv9f9/vv9s01+zA76VEraPBE32/fUChD9HJ8AruvV7X//7w5Th+ixlFM+ptXwf8N/6/7/UX+s63sQ9d+gU76OeZjzXLSb17D+gUxQf7v+Xpyb93zf97/uUcT+jFMsbdHz0LXQKf0c3vy/2fSgzX6Kjy03/2ZqO7q4//f+X/fPkR+hod9w/e9LfaxBi19/zpf4+x7/vn0P/3QfoZEZ+eJuOV99P8eUs+nmBTmne6cd9TcnKef/8i54uxz39dv0zPkDJv89/f//9CRUs/v3+mW9TcvoRyohLq5co4K+eej8FXF/safip/nlv92N/qf/7i3Qy+MV+zp/tx/v+nuanPQ2/qR+z+u59/Rps+u8fU7KsnL5EpyumPX32fb8v9HT4g/76b6ICNG9HPxr0P9P30/HlYciGjoROmZ+fUMFUHZAPxz2L/p5Ff89R9DGjo4/6bU6nE/xf2ZB/5eV9/9z7/XPUf1AAQfA0759N/3vdj+P9dGv5/pjXsv/MoCb9c73mCR2/nVK6GwrEGZwSNNmHOvyMji90Oib9/45R3e919M84pRynBGKCv5d6mh/R8aZjzMDxMSVvEzp51Y76nEl/vTb0up/zDZ0N3lAOm02bOO/HMYmfNWVrlNOX+vsvepoYLCPmZfJZh1Xem/Zj67DAYciFTjh0Y/1zFG+rz9PhX9KtrU5F2sQM1k8o3fPF/rukzUsUr2z6awXA8qW44xUq2LKgs+2n0dddSse6dksKG72PTpbeBL63f556+gU6bKH8fZFujR5TjvgH6PS79sPAn9hxTq01lC7VkZC/phRP7fvxfCzm4XhnfR+L+NwAh3Zw0dP6hX5M4lUdq2k/7v9C2dv39/MzoCLecm3VW0nHL/X9foTCWGLGG4ofZpQcqB/fpnSXOOr9/f/vows0GbSASmIeimycUjxs0u1zdFj33/Xfax8/RzfXD/d/n1C6S7uo/2Bw4y06XeDzdnR0/gilj+X7t+lk5ZqOjq9SuvKLdM69ePmE8md+vx/XS3T2+4t0vGTA4X7/nbpL/b+k/Cllc97//Xb/+RfpZCcD1Y7hDh0OEDdqU3d9H9P
"text/plain": [
"<Figure size 1750x250 with 7 Axes>"
]
},
"metadata": {},
"output_type": "display_data"
}
],
"source": [
"fig, ax = plt.subplots(1, shape[0], figsize=(2.5 *shape[0],2.5))\n",
"for i in range(0, shape[0]):\n",
" ax[i].pcolormesh(cropOD[i], cmap='jet')\n",
"\n",
"plt.show()"
]
},
{
"cell_type": "code",
"execution_count": 43,
"metadata": {},
"outputs": [
{
"name": "stdout",
"output_type": "stream",
"text": [
"\n",
" \n",
"file high_cond_1 \n",
" image: 0\n",
"shape: (200, 150)\n"
]
},
{
"data": {
"image/png": "iVBORw0KGgoAAAANSUhEUgAABTkAAAHXCAYAAAB3dl6BAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjcuMSwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/bCgiHAAAACXBIWXMAAA9hAAAPYQGoP6dpAACDoklEQVR4nOzdd3wU1f7/8feGFFoKnUAIoSOCFGmCECIoCCooRVC60iyo1wbqFbgqoNdrRUFQqYqAIBcvIr1IFSlSlE4gREOoCTWkzO+P/Ha+m7C72fRMeD197EOyc/bsmZ1yZj5zis0wDEMAAAAAAAAAYFFe+V0AAAAAAAAAAMgOgpwAAAAAAAAALI0gJwAAAAAAAABLI8gJAAAAAAAAwNIIcgIAAAAAAACwNIKcAAAAAAAAACyNICcAAAAAAAAASyPICQAAAAAAAMDSCHICAAAAAAAAsLQCEeQcO3asbDabbDab1q1bl9/FyTWRkZHmeg4cODC/i5Nj2rVrZ65XdoSFhclmsyksLCxnCoYc98cff+iJJ55QzZo1Vbx4cXO7d+vWLb+LhgIip84HgKOcrj9zur6h/gLy3owZM8zzwowZM7KdX2G9Tnfm8uXLCgkJkc1m03/+85/8Lg48kJP1zMCBA819PTIyssDlB6DgWLdunXl8jx07Nte+Z+PGjbLZbPLx8dH+/fuzlVeWgpz2lXT28vLyUkBAgOrWrau+ffvqf//7X7YKCKBgWL9+ve688059/fXXOnr0qK5du5bfRQLynONDOWcvPz8/VahQQeHh4Xrrrbf0119/ZZin/cYlM6+goCCPynv+/HlNmzZNjz76qOrWrasyZcrIx8dHZcqUUf369TVgwADNmTNHV65cyeYvU7gsXrxYY8eO1dixY3Xx4sX8Lg7ywUcffaSxY8fqo48+yu+iIBvsx3FOBEELk7ffflvR0dGqUKGCRowY4TRNRnVTyZIlVaNGDfXo0UNz585VUlKS2+90vFHOzOv555/3aJ12796tcePGqV27dgoLC1PJkiVVtGhRBQcHq23btnrppZe0ZcuWzP5Ut6Tdu3ebx87u3bvzuzi3jLi4OM2fP18jRoxQixYtzGu2UqVKqWHDhnrqqae0ffv2bH/P0KFD0xxjuRm4srLk5GTt27dPM2bM0LPPPqu77rorTSOfzPxuO3fu1Oeff67BgwerWbNmN52j2rdvr4kTJyomJib3VqiAu/vuu3XPPfcoKSlJI0eOzFZe3jlUJpNhGLp06ZIOHjyogwcP6ptvvlF4eLgWLlyoMmXK5PTXAcgjzz77rK5fvy5J6t+/v9q1a6dSpUpJkoKDg/OzaECBcePGDcXGxio2NlYbNmzQu+++q88++0wDBgzI83K8/fbb+uijj3Tp0qWblp8/f17nz5/X/v37NWvWLJUsWVIvvviiXn75ZZUoUSJPy1oQLV68WDNnzpSU2kLF06AyCo+PPvpIJ06cUNWqVT0OsqDgGTdunCQpPDy80LfO9FR0dLQ+/vhjSdKrr76q4sWLZymfK1eu6NixYzp27JgWLlyoCRMm6L///a+qVauWk8XN0MGDB/XSSy+5bFgTExOjmJgY/fLLL/rPf/6j+vXra+LEierSpUueltNK7AFjKTXY3ahRo/wt0C3gvffe05tvvqmEhISbll28eFEXL17Unj17NHnyZPXt21dffPFFlo7ddevW6csvv8yJIhd6vXr10qJFi3Ikr86dO+v06dNOl9nPUWvWrNE777yj//znPxo6dGiOfK/VjB07VmvWrNGaNWu0fPlydezYMUv5ZDvI+cMPP6T5OyUlRWfPntXWrVs1d+5cXb9+XevXr1fXrl31yy+/OO3CaH9SVNiFhYXJMIz8LkaBRfeGguvUqVPau3evJKljx47mzT+QXmEeciS9Rx99VL17907z3rVr13T06FEtWLBAe/bs0ZUrVzRo0CBVqFBBnTp1yjDPL774QuXLl88wnY+Pj8tlsbGxevjhh7V582bzvfr166tTp06qWbOmSpcurYsXL+rEiRNasWKFduzYocuXL2vcuHEKCwsrkIEA6k8A6Q0cODBHz1e3ynlmwoQJun79ukqVKqXhw4d79Jn0dVNKSoouXryonTt3as6cOYqLi9PevXt17733as+ePRkGX26//Xa9/fbbHn13jRo1XC776aef1KdPH8XHx0uSfH19dc899yg8PFzBwcEqWrSoTp8+rb1792rZsmWKjo7Wvn379MADD1huW3OfVLgdOnTIDHBWr15dHTp0UKNGjVS2bFlduHBBq1ev1sKFC5WcnKw5c+YoNjZWy5Ytk5eX5x1zr127piFDhsgwDJUoUYJePBlITk5O83fp0qVVpkwZHT58OEv5lS1bVi1btlTDhg1VrVo1BQYG6saNGzpy5IgWL16sXbt26fLlyxo2bJi8vb01ePDgnFiNbGvXrl2enS/btGmju+++Wxs3btSYMWOyHOSUkQWSzJc7e/fuNfz9/c20P/30U1a+DgVceHi4R/sDrGvz5s3mNh41alR+FwfIN2PGjDGPhTFjxrhMl5ycbDzxxBNm2ubNm7tMW7VqVTPd8ePHs1W+a9euGU2aNDHzCw0NNZYuXer2M3/++afRq1cvQ5Ixffr0bH2/Vdh/86pVqzpdPmDAgExtk4zyyw/2Mg0YMCC/i2JJBXGbIvPsx3F4eHh+F6VAOHfunFG8eHFDkjF06FC3aT2tm6KioowqVaqYaT///HOn6dauXZuj22PTpk2Gj4+PmWf37t2NEydOuEyfkpJizJs3z6hTp84tf8+SUR03ffp0c7kn1wWZrTNxsyeffNLo0qWLsW7dOpdpNmzYYJQsWdL8rb/++utMfccrr7xiSDIqV65svPDCCx5dz97K3nnnHWPUqFHGggULjGPHjhmGkfbYyMzvtm/fPiMlJcVtmvHjx5t5BwUFGdevX89O8S1r6tSp5u+wadOmLOWRqxMP1a9fX08++aT59/r163Pz6wDkEseuE35+fvlYEsAavLy89O9//1tFihSRJP366695Mo7tSy+9pJ07d0qSatWqpa1bt6pz585uP1O3bl3NmzdPc+fOlb+/f66XEQCQP2bNmqWrV69Kkvr27ZsjeYaEhOill14y/86L+724uDj17NlTiYmJkqRnnnlGCxYsUGhoqMvP2Gw29erVS7t27VL//v1zvYxAZrz33nv63//+p/DwcJdp2rRpowkTJph/Z2as4Z07d5qTjH3yyScKCAjIcllvFa+99pomTJigHj16ZHsYjttvvz3DSVlHjx6tO+64Q1LqEAWbNm3K1ndaVY8ePcx4w+eff56lPHJ9dvW6deua/46Li3OaJqPZ1Z3Ndnj27FmNHTtWDRo0kL+/v/z9/dWkSRNNmDDBrLxduXbtmn744Qc9/fTTaQb1DQwM1O23364RI0bo999/z3Dd0s8inJycrBkzZujee+9V5cqV5e3tnWYGvIxmbczs5BPuZte7ceOGvvrqKz300EOqUqWKihYtqqCgIN1xxx168cUXPe7ycOXKFY0fP15NmjRRQECAAgMDdccdd2jMmDE6e/asR3l4KqNZA53tJ6tXr1b37t3NdaxRo4aGDh160/pdv35dX3zxhVq1aqVy5cqpePHiatCggSZOnOh07BNH8fHx+uabb/TEE0+ocePGCgoKko+Pj0qXLq0mTZroxRdf1NGjRz1ez6NHj+qpp55SzZo1VaxYMZUvX15t27bV1KlTzWbx9vVs165dhvlt2bJFI0aMUL169RQUFKSiRYsqNDRUjz76qJYuXepxuZyx7+MRERHme+PGjbtpX3QmKipKo0aNUpMmTVS6dGn5+fmpcuXKevDBBzVjxoybugC4+m5Pj6/MOn36tF555RXVq1dPJUqUUOnSpdW8eXP95z//Mc8hGe2TmZnpMjNp9+/fr3/84x9q1KhRmt/uoYce0jfffKOUlJQM8/jvf/+rXr16qXr16ipevLiKFi2qypUrq2HDhurZs6c+//xznTt3zulnL168qHfffVfh4eEqX768fH19FRAQoOrVq+uuu+7S008/rWXLljntupDR7OrOZsP9448
"text/plain": [
"<Figure size 640x480 with 1 Axes>"
]
},
"metadata": {},
"output_type": "display_data"
},
{
"name": "stdout",
"output_type": "stream",
"text": [
"y smaller x, 1d fit along y\n",
"\n",
"1d fit initialization\n",
"center = [99.49748744 75.5033557 ]\n",
"BEC widths: [24 13]\n",
"\n",
"1d init fit values\n",
"Name Value Min Max Stderr Vary Expr Brute_Step\n",
"amp_bec 1.805 0 4.694 None True None None\n",
"amp_th 1.805 0 4.694 None True None None\n",
"deltax 39 0 200 None True None None\n",
"sigma_bec 10.66 0 26 None True None None\n",
"sigma_th 26.94 0 inf None False 0.632*sigma_bec + 0.518*deltax None\n",
"x0_bec 75.5 65.5 85.5 None True None None\n",
"x0_th 75.5 65.5 85.5 None True None None\n",
"1d fitted values\n",
"x0_bec: 75.823, (init = 75.503), bounds = [65.50 : 85.50] \n",
"x0_th: 75.762, (init = 75.503), bounds = [65.50 : 85.50] \n",
"amp_bec: 2.929, (init = 1.805), bounds = [0.00 : 4.69] \n",
"amp_th: 0.722, (init = 1.805), bounds = [0.00 : 4.69] \n",
"sigma_bec: 10.851, (init = 10.656), bounds = [0.00 : 26.00] \n",
"sigma_th: 9.869, (init = 26.936), bounds = [0.00 : inf] \n",
"\n"
]
},
{
"data": {
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAh8AAAHXCAYAAAAcDjiVAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjcuMSwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/bCgiHAAAACXBIWXMAAA9hAAAPYQGoP6dpAACJqklEQVR4nO3deXxTVdrA8d9N0r1NukOhLWUXZF8URQEFREVxA5dxA0VABx13ZXQEFEUHZ95xGUVFQUdFUdFh3FkFRJEdVzZb2gKldEvTLc1y3j/SxJYudEmbtH2+n0+UJuee+9wkTZ6ee+55NKWUQgghhBCiheh8HYAQQggh2hdJPoQQQgjRoiT5EEIIIUSLkuRDCCGEEC1Kkg8hhBBCtChJPoQQQgjRoiT5EEIIIUSLkuRDCCGEEC1Kkg8hhBBCtChJPtqhDRs2oGkamqYxb968JveXl5fHI488wuDBg4mIiECn06FpGpGRkZ42Y8aM8eyzNfvll1+49dZb6dGjB6GhoZ5juvzyy5ttn23lufMn7udzzJgxvg5F+DFvf1aKPxh8HUB753A4+PXXX9m+fTs7duxg+/bt7Nmzh9LSUgDmzp3r12/648ePM2LECNLS0prUz+7du/nkk08AuPzyyxk0aFCTY/O2b775hgsvvJCysjJfh9IsCgoK+Ne//gXAoEGDmjWhEkK0b5J8+NjVV1/NypUrfR1Goz355JOexGPkyJHccMMNdOjQAU3TCAgIqHc/u3fvZv78+QCkpKT4ZfJx5513ehKPm266iTFjxhAVFQVAQkKCL0PzioKCAs9rcPPNN0vyIYRoNpJ8+JjD4ajyc3R0NDExMRw4cMBHETXM559/DkBUVBRff/01oaGhNbbbsGFDC0blfZmZmfz4448ATJgwgTfffNPHEQkhmtuYMWOQ2qvNQ5IPHzvjjDPo06cPQ4cOZejQoXTt2pVly5Yxbdo0X4dWLxkZGQD07t271sSjLXAfJ8DgwYN9GIkQQrR+knz42F//+ldfh9Ak5eXlAAQFBfk4kuZltVo9/27rxyqEEM1NrnZpgzZv3sy1115L586dCQ4OJikpicsvv5wvv/zSK/3Pmzev2tUX33zzjec+963yqZbarthYtmwZmqZVGemZNm1atb5SUlKaFHNGRgYPP/wwQ4YMITo6mqCgIDp37syll17KsmXLqp3+Ojnu8847z3Pf/Pnzq8XXWMXFxTz11FMMGTIEo9GIyWRiwIABzJ07l5ycnHr3k56ezosvvsiUKVPo3bs34eHhBAYGEh8fz5gxY3jmmWcwm801bpuWloamaXTt2tVz35tvvlntGDVNqzax+MSJEyxZsoQbbriBfv36YTQaCQgIIDY2lhEjRvDYY4+RlZXVqOfG28fZGEopVqxYweTJk0lOTiY4OJjIyEgGDBjAvffee8rTozVdLZGens59993HaaedRlhYGJGRkZx99tm89NJL2O32esX18ccfM3HiRDp06EBwcDApKSnccMMNbN26Ffjj90rTNJYtW9bg4/7yyy89299111312uauu+7ybPPFF180eJ9uTX2N9+3bR1hYGJqmYTKZSE1NrbVtfn4+ycnJaJqGXq9n48aNVR6v79UuO3fuZNasWfTv39/zOxAfH0/fvn258MILeeKJJ1rNqfQWo4TfWbp0qQIUoObOndugbf/6178qTdM82598mz17tlq/fn2j+1dKqblz59baf+Xb+vXrPduMHj3ac39tx1rXrUuXLg2O023x4sUqJCSkzv779++vUlNTq21bOe66bo2xf/9+1bVr11r7TExMVLt27ar1uXNbv359na+5+xYXF6c2bdpUbfvU1NR6HSNQ5Tk6dOiQ0uv1p9wmLCxMffTRR416jrx5nJW5244ePbrWNllZWeqss86qc18BAQFq4cKFdcZc+Xftiy++UJGRkbX2N378eFVWVlZrf+Xl5WrKlCm1bq/X69Wzzz5b5fdq6dKlp3pqq3E4HJ73ZlRUlCotLa2zfWlpqYqKilKASk5OVg6Ho8H7VMp7r/Grr77qaXvWWWcpu91eY7vJkyd72j3yyCM1xnOqz8q5c+fWK+bLLrusMU9JmyWnXdqQRYsW8dRTTwGudQyuueYaxo8fT0hICHv27GHJkiW8+OKLHDlypEn7ufbaaz1Xo1xxxRUAnH766SxYsKBKu379+p2yr/PPP5+PP/6YdevW8cILLwCuq0rOP//8Ku0aO5/klVdeYdasWZ6fL730UiZOnEhkZCT79+9n6dKlpKam8uOPP3LOOeewa9cu4uLiPO0XLFhATk4OP/30E3/7298AuOaaa7j22msbFY9bXl4e559/PpmZmQB07dqVW265hZ49e5KXl8fKlStZs2YNV155JSaTqc6+ysrKUEpx+umnc95559GnTx9iYmIoKysjIyODTz75hB07dnDixAkuueQSdu/eXWUkKT4+no8//pjs7GxmzpwJwHnnnVfjX7zx8fGef5eXl+NwOOjWrRtjx46lX79+xMXF4XQ6SU9PZ82aNaxbt47i4mKuvfZavv32W4YPH97o56ypx9kQFouFUaNGsX//fsB1NdMtt9zC6aefTklJCatXr+aDDz7AZrMxZ84cnE7nKU+h7t69m0WLFqGUYubMmZx11lkEBQWxfft2Fi9eTHFxMatXr+bJJ5/k8ccfr7GPGTNm8MEHHwAQHBzM1KlTOeuss9Dr9Wzfvp3XX3+d+++/n8mTJzfquN10Oh233XYbf/3rX8nPz+ejjz7i+uuvr7X9hx9+SH5+PgC33norOl3jBtW99RrfdtttfPXVV3z00Ud89913zJ8/v9pz+vrrr/Phhx8CMGLEiEYtafDf//7Xc4VYSEgI1113HSNGjCA6OpqysjIyMzPZvn07q1evbnDfbZ5vcx9Rk8aMfBw8eFAFBQUpQAUGBqrPP/+8Wptjx46p0047rUo23piRj8rc/dT1F6RStY98uDX1L7WapKamekY89Hq9ev/996u1KSkpURMnTvTse/LkyTX21dTRopPdeuutnv4mTJigiouLq7X517/+Va8RlrS0NLV379469/fuu+8qnU6nADV16tQa21QeAbn55ptPeQy5ubnq22+/rbPNunXrVFhYmALUmDFjTtlnXbx1nEqd+n07a9YsT5tzzjlHFRQUVGvz1VdfqeDgYAUog8Ggdu/eXa1N5fcNFaMC+/fvr9Zu69atymAweEYaahr9WLNmjaef2NhY9eOPP1Zrk5qaqrp06VJln439fcrKylIBAQH1+v0eNWqU5/csIyOjUftTyruvcV5enkpMTPTEVXmkZN++fZ73pdFoVIcOHaqxj1P93rs/O/R6fZ2/C6WlpWrr1q11Hld7I8mHH2pM8vGXv/zFs828efNqbbdnz54qQ+VtOfm45557PH0++OCDtbYzm80qISFBAUrTtBq/HLyZfBw/flwFBgYqQMXExKjc3Nxa21YeFm7q3wo33XSTAlRISIgqLy+v9nhDk4/6euyxxzz9pqene63f2pzqOJWq+32bnZ3tSeSNRqM6duxYrftatGiRp6/rr7++2uMnJx8bN26sta/rr7++znaVk+T33nuv1n5O3mdTfp8qn+Kp6fdCKdcXubvNxIkTG72vhqjPa6yUUhs2bPAkKl26dFH5+fmqvLxcDR061BPzf/7zn1q3P9Xvfe/evRWgBgwY4I3Daldkwmkb8fHHHwMQEBDAnXfeWWu7AQMGcMEFF7RUWD7lXrzNYDBw33331drOaDRyxx13AKCU8jyXzeWzzz7zXCU0depUoqOja237wAMPeG2/Z599NgClpaXs3bvXa/3Wd78AP/zwQ4vtr7HH+dlnn3mubrr55pvp2LFjrW3vuOMOIiIiAFi1alWtE5fBdYn2ueeeW+vjlU81/vLLL1UeKysr4+uvvwagU6dOTJkypdZ+xowZw4ABA2p9vCEqn7JcsmRJjW0q3z9jxgyv7PdU6vsajx49mocffhiAw4cPM2vWLB599FF27NgBwPXXX88NN9zQ6Djcp4MzMzO9OtG5PZDkow3Izs4mPT0dcH3A1fVlBjB27NiWCMunsrOzOXz
"text/plain": [
"<Figure size 640x480 with 1 Axes>"
]
},
"metadata": {},
"output_type": "display_data"
},
{
"name": "stdout",
"output_type": "stream",
"text": [
"\n",
"Init Params\n",
"Name Value Min Max Stderr Vary Expr Brute_Step\n",
"amp_bec 3.515 0 6.495 None True None None\n",
"amp_th 0.866 0 6.495 None True None None\n",
"atom_number_bec 942.7 -inf inf None False amp_bec / 5 * 2 * 3.14159265359 * sigmax_bec * sigmay_bec None\n",
"atom_number_th 387.7 -inf inf None False amp_th * 2 * 3.14159265359 * 1.20206 / 1.643 * sigma_th * sigma_th None\n",
"condensate_fraction 0.7086 -inf inf None False atom_number_bec / (atom_number_bec + atom_number_th) None\n",
"rot_angle -36 -66 -6 None False None None\n",
"sigma_th 9.869 0 200 None True None None\n",
"sigmax_bec 19.67 0 48 None True None None\n",
"sigmay_bec 10.85 0 26 None True None None\n",
"x0_bec 99.5 89.5 109.5 None True None None\n",
"x0_th 99.5 89.5 109.5 None True None None\n",
"y0_bec 75.5 65.5 85.5 None True None None\n",
"y0_th 75.5 65.5 85.5 None True None None\n",
"\n",
"Name Value Min Max Stderr Vary Expr Brute_Step\n",
"amp_bec 3.515 0 6.495 None True None None\n",
"amp_th 0.866 0 6.495 None True None None\n",
"atom_number_bec 942.7 -inf inf None False amp_bec / 5 * 2 * 3.14159265359 * sigmax_bec * sigmay_bec None\n",
"atom_number_th 387.7 -inf inf None False amp_th * 2 * 3.14159265359 * 1.20206 / 1.643 * sigma_th * sigma_th None\n",
"condensate_fraction 0.7086 -inf inf None False atom_number_bec / (atom_number_bec + atom_number_th) None\n",
"rot_angle -36 -66 -6 None False None None\n",
"sigma_th 9.869 0 200 None True None None\n",
"sigmax_bec 19.67 0 48 None True None None\n",
"sigmay_bec 10.85 0 26 None True None None\n",
"x0_bec 99.5 89.5 109.5 None True None None\n",
"x0_th 99.5 89.5 109.5 None True None None\n",
"y0_bec 75.5 65.5 85.5 None True None None\n",
"y0_th 75.5 65.5 85.5 None True None None\n",
"bval first fit\n",
"amp_bec: 3.361, (init = 3.515), bounds = [0.00 : 6.50] \n",
"amp_th: 1.033, (init = 0.866), bounds = [0.00 : 6.50] \n",
"x0_bec: 100.041, (init = 99.497), bounds = [89.50 : 109.50] \n",
"y0_bec: 76.006, (init = 75.503), bounds = [65.50 : 85.50] \n",
"x0_th: 99.794, (init = 99.497), bounds = [89.50 : 109.50] \n",
"y0_th: 75.935, (init = 75.503), bounds = [65.50 : 85.50] \n",
"sigmax_bec: 23.817, (init = 19.672), bounds = [0.00 : 48.00] \n",
"sigmay_bec: 11.257, (init = 10.851), bounds = [0.00 : 26.00] \n",
"sigma_th: 9.384, (init = 9.869), bounds = [0.00 : 200.00] \n",
"rot_angle: -36.000, (init = -36.000), bounds = [-66.00 : -6.00] \n",
"\n",
"amp_bec: 3.361, (init = 3.515), bounds = [0.00 : 6.50] \n",
"amp_th: 1.033, (init = 0.866), bounds = [0.00 : 6.50] \n",
"x0_bec: 100.041, (init = 99.497), bounds = [89.50 : 109.50] \n",
"y0_bec: 76.006, (init = 75.503), bounds = [65.50 : 85.50] \n",
"x0_th: 99.794, (init = 99.497), bounds = [89.50 : 109.50] \n",
"y0_th: 75.935, (init = 75.503), bounds = [65.50 : 85.50] \n",
"sigmax_bec: 23.817, (init = 19.672), bounds = [0.00 : 48.00] \n",
"sigmay_bec: 11.257, (init = 10.851), bounds = [0.00 : 26.00] \n",
"sigma_th: 9.384, (init = 9.869), bounds = [0.00 : 200.00] \n",
"rot_angle: -36.000, (init = -36.000), bounds = [-66.00 : -6.00] \n",
"\n",
"\n",
"Atom numbers:\n",
" N_bec: 80596\n",
" N_th: 29742\n",
" N_ges: 110338\n",
" Cond. frac: 73.04 %\n",
"\n",
"Temperature: 31.40 nK\n",
"\n",
" \n",
"file high_cond_1 \n",
" image: 1\n",
"shape: (200, 150)\n"
]
},
{
"data": {
"image/png": "iVBORw0KGgoAAAANSUhEUgAABTkAAAHXCAYAAAB3dl6BAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjcuMSwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/bCgiHAAAACXBIWXMAAA9hAAAPYQGoP6dpAACEUElEQVR4nOzdd3wU1f7/8feGFFoSOgRCCBCqIF0QhBBRQVBBmiIdpVmwK6hX4Ipiu1YUFJWqCAgiXi7Se1dAsNAJhGAINaGGhMzvj3x3fptkd7PpmfB6Ph77eCQ7Z8+e2Sln5jOn2AzDMAQAAAAAAAAAFuWV3wUAAAAAAAAAgOwgyAkAAAAAAADA0ghyAgAAAAAAALA0gpwAAAAAAAAALI0gJwAAAAAAAABLI8gJAAAAAAAAwNIIcgIAAAAAAACwNIKcAAAAAAAAACyNICcAAAAAAAAASysQQc5x48bJZrPJZrNp7dq1+V2cXBMZGWmu56BBg/K7ODmmffv25nplR2hoqGw2m0JDQ3OmYMhxf/31lx599FGFhYWpePHi5nbv1q1bfhcNBUROnQ8ARzldf+Z0fUP9BeS96dOnm+eF6dOnZzu/wnqd7sylS5cUHBwsm82m//znP/ldHHggJ+uZQYMGmft6ZGRkgcsPQMGxdu1a8/geN25crn3Pxo0bZbPZ5OPjoz///DNbeWUpyGlfSWcvLy8vBQQEqG7duurXr5/++9//ZquAAAqGdevWqVmzZvrmm290+PBhXb16Nb+LBOQ5x4dyzl5+fn6qWLGiwsPD9cYbb+jkyZMZ5mm/ccnMq1SpUh6V99y5c5o6daoeeugh1a1bV2XLlpWPj4/Kli2rBg0aaODAgZo9e7YuX76czV+mcFm0aJHGjRuncePG6cKFC/ldHOSDjz76SOPGjdNHH32U30VBNtiP45wIghYmEyZMUHR0tCpWrKiRI0c6TZNR3VSyZEnVrFlTPXv21Jw5c5SUlOT2Ox1vlDPzeuaZZzxap927d2v8+PFq3769QkNDVbJkSRUtWlRBQUFq166dXnjhBW3ZsiWzP9VNaffu3eaxs3v37vwuzk0jLi5O8+bN08iRI9WyZUvzmq106dJq1KiRHn/8ce3YsSNLee/Zs0cvvfSSmjRpovLly8vPz0/BwcFq2bKlnnvuOS1dujSH18b6bty4oT/++EPTp0/XU089pdtvvz1VI5+cCvgNGzYs1TkvNwOJBdkdd9yhO++8U0lJSRo1alS28vLOoTKZDMPQxYsXtX//fu3fv1/ffvutwsPDtWDBApUtWzanvw5AHnnqqad07do1SdKAAQPUvn17lS5dWpIUFBSUn0UDCozr168rNjZWsbGxWr9+vd555x199tlnGjhwYJ6XY8KECfroo4908eLFdMvPnTunc+fO6c8//9TMmTNVsmRJPf/883rxxRdVokSJPC1rQbRo0SLNmDFDUkoLFU+Dyig8PvroIx07dkzVqlXzOMiCgmf8+PGSpPDw8ELfOtNT0dHR+vjjjyVJL7/8sooXL56lfC5fvqwjR47oyJEjWrBggSZOnKiffvpJ1atXz8niZmj//v164YUXXDasiYmJUUxMjDZs2KD//Oc/atCggd5++2116dIlT8tpJfaAsZQS7G7cuHH+Fugm8O677+r1119XQkJCumUXLlzQhQsXtGfPHk2ePFn9+vXTF1984dGxe+XKFT333HOaOnWqkpOTUy2Ljo5WdHS0tm/frm+++YaHumn07t1bCxcuzNXvWLt2rb766qtc/Q4rGTdunFavXq3Vq1dr2bJl6tixY5byyXaQ88cff0z1f3Jyss6cOaOtW7dqzpw5unbtmtatW6euXbtqw4YNTrsw2p8UFXahoaEyDCO/i1Fg0b2h4Dpx4oT27t0rSerYsaN58w+kVZiHHEnroYce0sMPP5zqvatXr+rw4cOaP3++9uzZo8uXL2vw4MGqWLGiOnXqlGGeX3zxhSpUqJBhOh8fH5fLYmNj9eCDD2rz5s3mew0aNFCnTp0UFhamMmXK6MKFCzp27JiWL1+u3377TZcuXdL48eMVGhpaIAMB1J8A0ho0aFCOnq9ulvPMxIkTde3aNZUuXVojRozw6DNp66bk5GRduHBBO3fu1OzZsxUXF6e9e/fq7rvv1p49ezIMvtxyyy2aMGGCR99ds2ZNl8v+97//qU+fPoqPj5ck+fr66s4771R4eLiCgoJUtGhRnTp1Snv37tXSpUsVHR2tP/74Q/fdd5/ltjX3SYXbgQMHzABnjRo1dNddd6lx48YqV66czp8/r1WrVmnBggW6ceOGZs+erdjYWC1dulReXq475l66dEn33Xef1q1bJ0kKCQlRjx491KBBAwUEBCguLk779u3TL7/8ohMnTuTJelrJjRs3Uv1fpkwZlS1bVgcPHsyR/K9evaqhQ4fKMAyVKFGiwPaqat++fZ6dL9u2bas77rhDGzdu1NixY/MvyOlqLL5hw4bpueeeU+vWrXXx4kVt2rRJv/zyi+69997sfiWAPBYVFWX+3aRJk3wsCVBw1K1b12Ud+Morr2jYsGH6+uuvZRiGxo4d61GQ85577snWeFvXrl3Tvffeq507d0pKuaCdPHmyOnfu7DT9hAkTtG/fPo0dO1bz5s3L8vei4AkNDdWxY8c0cOBAuuoCkJTSin/atGmSpF69eqlYsWIefc5V3TRkyBCNHj1arVu3VlRUlA4fPqwZM2a47AJvV65cuWyP575582Z169ZNiYmJkqQePXrogw8+UEhIiNP0hmFo/vz5ev3117V///5sfTeQ02w2m7p06aIXX3xR4eHh6ZYPGzZMGzZsUOfOnXXp0iUtX75cM2bM0ODBg13mOWLECDPA+corr2js2LHy9fVNl+69995Lda+HFLfddpvq1aunZs2aqVmzZqpevbqmT5/u9jfPjHHjxunQoUOqUqWKevfurQ8//DBH8rW6AQMGaOPGjdq2bZs2b96s1q1bZzqPXJ14qEGDBnrsscfM/+0HGQBrcew64efnl48lAazBy8tL7733nooUKSJJ2r59e56MY/vCCy+YAc5atWpp69atLgOcdnXr1tXcuXM1Z84c+fv753oZAQD5Y+bMmbpy5YokqV+/fjmSZ3BwsF544QXz/7y434uLi1OvXr3MAOeTTz6p+fPnuwxwSilBpN69e2vXrl0aMGBArpcRyIx3331X//3vf50GOO3atm2riRMnmv+7e4D5yy+/6Ntvv5UkPf3003rzzTedBjjtqlatmvlCF3KvvPKKJk6cqJ49e+b4MBw7d+40J3375JNPFBAQkKP5W1nPnj3NeMPnn3+epTxyfXb1unXrmn/HxcU5TZPR7OrOZjs8c+aMxo0bp4YNG8rf31/+/v5q2rSpJk6caFberly9elU//vijnnjiiVSD+gYGBuqWW27RyJEj9fvvv2e4bmlnEb5x44amT5+uu+++W1WqVJG3t3eqp54ZzdqY2ckn3LX2uX79ur7++ms98MADqlq1qooWLapSpUrp1ltv1fPPP+9xl4fLly/rrbfeUtOmTRUQEKDAwEDdeuutGjt2rM6cOeNRHp7KaNZAZ/vJqlWr1KNHD3Mda9asqWHDhqVbv2vXrumLL75Q69atVb58eRUvXlwNGzbU22+/7XTsE0fx8fH69ttv9eijj6pJkyYqVaqUfHx8VKZMGTVt2lTPP/+8Dh8+7PF6Hj58WI8//rjCwsJUrFgxVahQQe3atdOXX35pNou3r2f79u0zzG/Lli0aOXKk6tevr1KlSqlo0aIKCQnRQw89pCVLlnhcLmfs+3hERIT53vjx49Pti85ERUVp9OjRatq0qcqUKSM/Pz9VqVJF999/v6ZPn56uC4Cr7/b0+MqsU6dO6aWXXlL9+vVVokQJlSlTRrfddpv+85//mOeQjPbJzMx0mZm0f/75p5577jk1btw41W/3wAMP6Ntvv003po4zP/30k3r37q0aNWqoePHiKlq0qKpUqaJGjRqpV69e+vzzz3X27Fmnn71w4YLeeecdhYeHq0KFCvL19VVAQIBq1Kih22+/XU888YSWLl3qtOtCRrOrO5s
"text/plain": [
"<Figure size 640x480 with 1 Axes>"
]
},
"metadata": {},
"output_type": "display_data"
},
{
"name": "stdout",
"output_type": "stream",
"text": [
"y smaller x, 1d fit along y\n",
"\n",
"1d fit initialization\n",
"center = [99.49748744 75.5033557 ]\n",
"BEC widths: [26 14]\n",
"\n",
"1d init fit values\n",
"Name Value Min Max Stderr Vary Expr Brute_Step\n",
"amp_bec 1.733 0 4.505 None True None None\n",
"amp_th 1.733 0 4.505 None True None None\n",
"deltax 42 0 200 None True None None\n",
"sigma_bec 11.48 0 28 None True None None\n",
"sigma_th 29.01 0 inf None False 0.632*sigma_bec + 0.518*deltax None\n",
"x0_bec 75.5 65.5 85.5 None True None None\n",
"x0_th 75.5 65.5 85.5 None True None None\n",
"1d fitted values\n",
"x0_bec: 75.810, (init = 75.503), bounds = [65.50 : 85.50] \n",
"x0_th: 75.715, (init = 75.503), bounds = [65.50 : 85.50] \n",
"amp_bec: 2.951, (init = 1.733), bounds = [0.00 : 4.51] \n",
"amp_th: 0.574, (init = 1.733), bounds = [0.00 : 4.51] \n",
"sigma_bec: 11.042, (init = 11.475), bounds = [0.00 : 28.00] \n",
"sigma_th: 9.725, (init = 29.008), bounds = [0.00 : inf] \n",
"\n"
]
},
{
"data": {
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAh8AAAHXCAYAAAAcDjiVAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjcuMSwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/bCgiHAAAACXBIWXMAAA9hAAAPYQGoP6dpAACLIUlEQVR4nO3deXhTVfrA8e9N0nRPutFCaUvZF9kFRVBBAXVkBAVxGRwFF2AcdMZ9HQV1RH/qjM44iiMO6KgoKjq4IquAKLKDK1tLW6CFbumeZjm/P9LEli50SZq0fT/PkweanHvve7PcvDn33PNqSimFEEIIIUQr0fk7ACGEEEJ0LJJ8CCGEEKJVSfIhhBBCiFYlyYcQQgghWpUkH0IIIYRoVZJ8CCGEEKJVSfIhhBBCiFYlyYcQQgghWpUkH0IIIYRoVZJ8dEAbNmxA0zQ0TWP+/PktXl9+fj4PPfQQw4YNIzIyEp1Oh6ZpREVFedqMGzfOs8227Mcff+Smm26iV69ehIWFefbp8ssv99k228tzF0jcz+e4ceP8HYoIYN4+VopfGfwdQEfncDj46aef2L59Ozt27GD79u3s2bOH8vJyAB599NGAftPn5OQwatQo0tPTW7Se3bt389FHHwFw+eWXM3To0BbH5m1fffUVl1xyCRUVFf4OxScKCwt5/vnnARg6dKhPEyohRMcmyYefXXXVVaxYscLfYTTbX//6V0/iMWbMGK677joSEhLQNI2goKBGr2f37t0sWLAAgNTU1IBMPm677TZP4nH99dczbtw4oqOjAejSpYs/Q/OKwsJCz2twww03SPIhhPAZST78zOFw1Pg7JiaG2NhYDhw44KeImuazzz4DIDo6mi+//JKwsLA6223YsKEVo/K+rKws9u3bB8DFF1/M66+/7ueIhBC+Nm7cOKT2qm9I8uFnZ511Fv379+fMM8/kzDPPpHv37ixdupRZs2b5O7RGyczMBKBv3771Jh7tgXs/AYYNG+bHSIQQou2T5MPPHnzwQX+H0CKVlZUABAcH+zkS37JarZ7/t/d9FUIIX5OrXdqhzZs3c80119C1a1dCQkJITk7m8ssv54svvvDK+ufPn1/r6ouvvvrKc5/7Vv1US31XbCxduhRN02r09MyaNavWulJTU1sUc2ZmJvfffz/Dhw8nJiaG4OBgunbtymWXXcbSpUtrnf46Ne4LLrjAc9+CBQtqxddcpaWlPPnkkwwfPhyTyYTZbGbw4ME8+uij5ObmNno9GRkZvPjii0yfPp2+ffsSERGB0WgkPj6ecePG8fTTT2OxWOpcNj09HU3T6N69u+e+119/vdY+appWa2DxyZMnWbx4Mddddx0DBw7EZDIRFBREXFwco0aN4pFHHiE7O7tZz42397M5lFIsX76cK6+8kpSUFEJCQoiKimLw4MHceeedpz09WtfVEhkZGdx1113069eP8PBwoqKiGD16NC+99BJ2u71RcX344YdMmjSJhIQEQkJCSE1N5brrrmPr1q3Ar58rTdNYunRpk/f7iy++8Cx/++23N2qZ22+/3bPM559/3uRturX0Nf7ll18IDw9H0zTMZjNpaWn1ti0oKCAlJQVN09Dr9WzcuLHG44292mXnzp3MnTuXQYMGeT4D8fHxDBgwgEsuuYTHH3+8zZxKbzVKBJwlS5YoQAHq0UcfbdKyDz74oNI0zbP8qbd58+ap9evXN3v9Sin16KOP1rv+6rf169d7lhk7dqzn/vr2taFbt27dmhyn26JFi1RoaGiD6x80aJBKS0urtWz1uBu6Ncf+/ftV9+7d611nUlKS2rVrV73Pndv69esbfM3dt06dOqlNmzbVWj4tLa1R+wjUeI4OHTqk9Hr9aZcJDw9XH3zwQbOeI2/uZ3XutmPHjq23TXZ2tjrnnHMa3FZQUJBauHBhgzFX/6x9/vnnKioqqt71TZw4UVVUVNS7vsrKSjV9+vR6l9fr9erZZ5+t8blasmTJ6Z7aWhwOh+e9GR0drcrLyxtsX15erqKjoxWgUlJSlMPhaPI2lfLea/zvf//b0/acc85Rdru9znZXXnmlp91DDz1UZzynO1Y++uijjYp5ypQpzXlK2i057dKOPPPMMzz55JOAax6Dq6++mokTJxIaGsqePXtYvHgxL774IkePHm3Rdq655hrP1ShXXHEFAGeccQZPPPFEjXYDBw487bouvPBCPvzwQ9atW8c///lPwHVVyYUXXlijXXPHk7zyyivMnTvX8/dll13GpEmTiIqKYv/+/SxZsoS0tDT27dvHueeey65du+jUqZOn/RNPPEFubi7ff/89f/nLXwC4+uqrueaaa5oVj1t+fj4XXnghWVlZAHTv3p0bb7yR3r17k5+fz4oVK1izZg1Tp07FbDY3uK6KigqUUpxxxhlccMEF9O/fn9jYWCoqKsjMzOSjjz5ix44dnDx5kt/+9rfs3r27Rk9SfHw8H374ISdOnGDOnDkAXHDBBXX+4o2Pj/f8v7KyEofDQY8ePRg/fjwDBw6kU6dOOJ1OMjIyWLNmDevWraO0tJRrrrmGr7/+mpEjRzb7OWvpfjZFcXEx559/Pvv37wdcVzPdeOONnHHGGZSVlbF69Wree+89bDYbDzzwAE6n87SnUHfv3s0zzzyDUoo5c+ZwzjnnEBwczPbt21m0aBGlpaWsXr2av/71rzz22GN1rmP27Nm89957AISEhDBz5kzOOecc9Ho927dv57XXXuPuu+/myiuvbNZ+u+l0Om655RYefPBBCgoK+OCDD5gxY0a97d9//30KCgoAuOmmm9Dpmtep7q3X+JZbbmHVqlV88MEHfPPNNyxYsKDWc/raa6/x/vvvAzBq1KhmTWnwv//9z3OFWGhoKNdeey2jRo0iJiaGiooKsrKy2L59O6tXr27yuts9/+Y+oi7N6fk4ePCgCg4OVoAyGo3qs88+q9Xm+PHjql+/fjWy8eb0fFTnXk9DvyCVqr/nw62lv9TqkpaW5unx0Ov16t13363VpqysTE2aNMmz7SuvvLLOdbW0t+hUN910k2d9F198sSotLa3V5vnnn29UD0t6errau3dvg9t7++23lU6nU4CaOXNmnW2q94DccMMNp92HvLw89fXXXzfYZt26dSo8PFwBaty4caddZ0O8tZ9Knf59O3fuXE+bc889VxUWFtZqs2rVKhUSEqIAZTAY1O7du2u1qf6+oapXYP/+/bXabd26VRkMBk9PQ129H2vWrPGsJy4uTu3bt69Wm7S0NNWtW7ca22zu5yk7O1sFBQU16vN9/vnnez5nmZmZzdqeUt59jfPz81VSUpInruo9Jb/88ovnfWkymdShQ4fqXMfpPvfuY4der2/ws1BeXq62bt3a4H51NJJ8BKDmJB9/+tOfPMvMnz+/3nZ79uyp0VXenpOPO+64w7POe++9t952FotFdenSRQFK07Q6vxy8mXzk5OQoo9GoABUbG6vy8vLqbVu9W7ilvxWuv/56BajQ0FBVWVlZ6/GmJh+N9cgjj3jWm5GR4bX11ud0+6lUw+/bEydOeBJ5k8mkjh8/Xu+2nnnmGc+6ZsyYUevxU5OPjRs31ruuGTNmNNiuepL8zjvv1LueU7fZks9T9VM8dX0ulHJ9kbvbTJo0qdnbaorGvMZKKbVhwwZPotKtWzdVUFCgKisr1ZlnnumJ+b///W+9y5/uc9+3b18FqMGDB3tjtzoUGXDaTnz44YcABAUFcdttt9XbbvDgwVx00UWtFZZfuSdvMxgM3HXXXfW2M5lM3HrrrQAopTzPpa98+umnnquEZs6cSUxMTL1t77nnHq9td/To0QCUl5ezd+9er623sdsF+O6771pte83dz08//dRzddMNN9xA586d62176623EhkZCcDKlSvrHbgMrku0zzvvvHofr36q8ccff6zxWEVFBV9++SUAiYmJTJ8+vd71jBs3jsGDB9f7eFNUP2W5ePHiOttUv3/27Nle2e7pNPY1Hjt2LPfffz8AR44cYe7cuTz88MPs2LEDgBkzZnDdddc1Ow736eCsrCyvDnTuCCT5aAdOnDhBRkYG4DrANfRlBjB+/PjWCMu
"text/plain": [
"<Figure size 640x480 with 1 Axes>"
]
},
"metadata": {},
"output_type": "display_data"
},
{
"name": "stdout",
"output_type": "stream",
"text": [
"\n",
"Init Params\n",
"Name Value Min Max Stderr Vary Expr Brute_Step\n",
"amp_bec 3.554 0 6.284 None True None None\n",
"amp_th 0.6913 0 6.284 None True None None\n",
"atom_number_bec 1051 -inf inf None False amp_bec / 5 * 2 * 3.14159265359 * sigmax_bec * sigmay_bec None\n",
"atom_number_th 300.6 -inf inf None False amp_th * 2 * 3.14159265359 * 1.20206 / 1.643 * sigma_th * sigma_th None\n",
"condensate_fraction 0.7776 -inf inf None False atom_number_bec / (atom_number_bec + atom_number_th) None\n",
"rot_angle -36 -66 -6 None False None None\n",
"sigma_th 9.725 0 200 None True None None\n",
"sigmax_bec 21.31 0 52 None True None None\n",
"sigmay_bec 11.04 0 28 None True None None\n",
"x0_bec 99.5 89.5 109.5 None True None None\n",
"x0_th 99.5 89.5 109.5 None True None None\n",
"y0_bec 75.5 65.5 85.5 None True None None\n",
"y0_th 75.5 65.5 85.5 None True None None\n",
"\n",
"Name Value Min Max Stderr Vary Expr Brute_Step\n",
"amp_bec 3.554 0 6.284 None True None None\n",
"amp_th 0.6913 0 6.284 None True None None\n",
"atom_number_bec 1051 -inf inf None False amp_bec / 5 * 2 * 3.14159265359 * sigmax_bec * sigmay_bec None\n",
"atom_number_th 300.6 -inf inf None False amp_th * 2 * 3.14159265359 * 1.20206 / 1.643 * sigma_th * sigma_th None\n",
"condensate_fraction 0.7776 -inf inf None False atom_number_bec / (atom_number_bec + atom_number_th) None\n",
"rot_angle -36 -66 -6 None False None None\n",
"sigma_th 9.725 0 200 None True None None\n",
"sigmax_bec 21.31 0 52 None True None None\n",
"sigmay_bec 11.04 0 28 None True None None\n",
"x0_bec 99.5 89.5 109.5 None True None None\n",
"x0_th 99.5 89.5 109.5 None True None None\n",
"y0_bec 75.5 65.5 85.5 None True None None\n",
"y0_th 75.5 65.5 85.5 None True None None\n",
"bval first fit\n",
"amp_bec: 3.445, (init = 3.554), bounds = [0.00 : 6.28] \n",
"amp_th: 0.879, (init = 0.691), bounds = [0.00 : 6.28] \n",
"x0_bec: 100.009, (init = 99.497), bounds = [89.50 : 109.50] \n",
"y0_bec: 75.995, (init = 75.503), bounds = [65.50 : 85.50] \n",
"x0_th: 99.838, (init = 99.497), bounds = [89.50 : 109.50] \n",
"y0_th: 75.953, (init = 75.503), bounds = [65.50 : 85.50] \n",
"sigmax_bec: 24.234, (init = 21.311), bounds = [0.00 : 52.00] \n",
"sigmay_bec: 11.466, (init = 11.042), bounds = [0.00 : 28.00] \n",
"sigma_th: 9.279, (init = 9.725), bounds = [0.00 : 200.00] \n",
"rot_angle: -36.000, (init = -36.000), bounds = [-66.00 : -6.00] \n",
"\n",
"amp_bec: 3.445, (init = 3.554), bounds = [0.00 : 6.28] \n",
"amp_th: 0.879, (init = 0.691), bounds = [0.00 : 6.28] \n",
"x0_bec: 100.009, (init = 99.497), bounds = [89.50 : 109.50] \n",
"y0_bec: 75.995, (init = 75.503), bounds = [65.50 : 85.50] \n",
"x0_th: 99.838, (init = 99.497), bounds = [89.50 : 109.50] \n",
"y0_th: 75.953, (init = 75.503), bounds = [65.50 : 85.50] \n",
"sigmax_bec: 24.234, (init = 21.311), bounds = [0.00 : 52.00] \n",
"sigmay_bec: 11.466, (init = 11.042), bounds = [0.00 : 28.00] \n",
"sigma_th: 9.279, (init = 9.725), bounds = [0.00 : 200.00] \n",
"rot_angle: -36.000, (init = -36.000), bounds = [-66.00 : -6.00] \n",
"\n",
"\n",
"Atom numbers:\n",
" N_bec: 85613\n",
" N_th: 24772\n",
" N_ges: 110385\n",
" Cond. frac: 77.56 %\n",
"\n",
"Temperature: 30.71 nK\n",
"\n",
" \n",
"file high_cond_1 \n",
" image: 2\n",
"shape: (200, 150)\n"
]
},
{
"data": {
"image/png": "iVBORw0KGgoAAAANSUhEUgAABTkAAAHXCAYAAAB3dl6BAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjcuMSwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/bCgiHAAAACXBIWXMAAA9hAAAPYQGoP6dpAACDY0lEQVR4nOzdd3xUVf7/8feENEoKnUAIoSOCFGnSQgQFQQWlKEpXmgV1baCuwoqirutaUBBUuggIsrgs0nsTKYIonUCIhlATaki5vz/ym/udhJnJpOeG1/PxmAdk7pkz584t597PPcVmGIYhAAAAAAAAALAor4IuAAAAAAAAAADkBEFOAAAAAAAAAJZGkBMAAAAAAACApRHkBAAAAAAAAGBpBDkBAAAAAAAAWBpBTgAAAAAAAACWRpATAAAAAAAAgKUR5AQAAAAAAABgaQQ5AQAAAAAAAFhaoQhyjh07VjabTTabTevWrSvo4uSZqKgocz0HDRpU0MXJNR06dDDXKyfCw8Nls9kUHh6eOwVDrvv999/1xBNPqFatWipRooS53Xv06FHQRUMhkVvnA8BRbtefuV3fUH8B+W/69OnmeWH69Ok5zq+oXqc7c/nyZYWGhspms+lf//pXQRcHHsjNembQoEHmvh4VFVXo8gNQeKxbt848vseOHZtn37Np0ybZbDb5+Pho//79OcorW0FO+0o6e3l5eSkwMFD16tVTv3799N///jdHBQRQOKxfv1533nmnvvnmGx09elTXrl0r6CIB+c7xoZyzl5+fnypWrKiIiAi9/fbb+vPPPzPN037jkpVXcHCwR+U9f/68pk6dqkceeUT16tVT2bJl5ePjo7Jly6pBgwYaOHCgZs+erStXruTwlylaFi9erLFjx2rs2LG6ePFiQRcHBeDjjz/W2LFj9fHHHxd0UZAD9uM4N4KgRcn48eMVExOjihUrauTIkU7TZFY3lSpVSjVr1lSvXr00d+5cJScnu/1OxxvlrLyef/55j9Zpz549GjdunDp06KDw8HCVKlVK/v7+CgkJUfv27fXSSy9p69atWf2pbkl79uwxj509e/YUdHFuGfHx8Zo/f75Gjhypli1bmtdspUuXVqNGjfTUU09px44dOf6eYcOGpTvG8jJwZWUpKSn67bffNH36dD377LO666670jXyycrvdu3aNf3444964YUX1LZtW1WoUEG+vr4KDAzUbbfdpsGDB2v16tV5tzIW0LZtW919991KTk7WqFGjcpSXdy6VyWQYhi5duqSDBw/q4MGDmjNnjiIiIrRw4UKVLVs2t78OQD559tlndf36dUnSgAED1KFDB5UuXVqSFBISUpBFAwqNGzduKC4uTnFxcdqwYYPef/99ff755xo4cGC+l2P8+PH6+OOPdenSpZuWnz9/XufPn9f+/fs1c+ZMlSpVSi+++KJefvlllSxZMl/LWhgtXrxYM2bMkJTWQsXToDKKjo8//lgnTpxQtWrVPA6yoPAZN26cJCkiIqLIt870VExMjD755BNJ0quvvqoSJUpkK58rV67o2LFjOnbsmBYuXKgJEyboP//5j6pXr56bxc3UwYMH9dJLL7lsWBMbG6vY2Fht3LhR//rXv9SgQQO999576tatW76W00rsAWMpLdjduHHjgi3QLeCDDz7Qm2++qcTExJuWXbx4URcvXtTevXs1adIk9evXT19++WW2jt1169bpq6++yo0iF3l9+vTRokWLcpzPnDlzNGLECF2+fPmmZUlJSTpw4IAOHDig6dOnq0uXLpo5c6bKly+f4++1orFjx2rNmjVas2aNli9frs6dO2crnxwHOX/44Yd0f6empurs2bPatm2b5s6dq+vXr2v9+vXq3r27Nm7c6LQLo/1JUVEXHh4uwzAKuhiFFt0bCq9Tp05p3759kqTOnTubN/9ARkV5yJGMHnnkET366KPp3rt27ZqOHj2qBQsWaO/evbpy5YoGDx6sihUrqkuXLpnm+eWXX6pChQqZpvPx8XG5LC4uTg899JC2bNlivtegQQN16dJFtWrVUpkyZXTx4kWdOHFCK1as0M6dO3X58mWNGzdO4eHhhTIQQP0JIKNBgwbl6vnqVjnPTJgwQdevX1fp0qU1YsQIjz6TsW5KTU3VxYsXtWvXLs2ePVvx8fHat2+f7rnnHu3duzfT4Mvtt9+u8ePHe/TdNWvWdLnsf//7n/r27auEhARJkq+vr+6++25FREQoJCRE/v7+On36tPbt26dly5YpJiZGv/32m+6//37LbWvuk4q2Q4cOmQHOGjVqqFOnTmrcuLHKlSunCxcuaPXq1Vq4cKFSUlI0e/ZsxcXFadmyZfLy8rxj7rVr1zR06FAZhqGSJUvSiycTKSkp6f4uU6aMypYtq8OHD2cpn+PHj5sBzpCQEN1zzz1q3ry5KlSooCtXrmjjxo1m3Oynn35Sp06dtHXr1mw/gMptHTp0yLfzZbt27dS2bVtt2rRJb731VraDnDKyQZL5cmffvn1GQECAmfZ///tfdr4OhVxERIRH+wOsa8uWLeY2Hj16dEEXBygwb731lnksvPXWWy7TpaSkGE888YSZtkWLFi7TVqtWzUx3/PjxHJXv2rVrRtOmTc38wsLCjKVLl7r9zB9//GH06dPHkGRMmzYtR99vFfbfvFq1ak6XDxw4MEvbJLP8CoK9TAMHDizoolhSYdymyDr7cRwREVHQRSkUzp07Z5QoUcKQZAwbNsxtWk/rpujoaKNq1apm2i+++MJpurVr1+bq9ti8ebPh4+Nj5tmzZ0/jxIkTLtOnpqYa8+bNM+rWrXvL37NkVsdNmzbNXO7JdUFW60zc7MknnzS6detmrFu3zmWaDRs2GKVKlTJ/62+++SZL3/HKK68YkowqVaoYL7zwgkfXs7eyd955xxg9erSxYMEC49ixY4ZhpD82PP3d3n77baNNmzbGjz/+aCQnJztNs3//fiMkJMTM+80338yt1bCcKVOmmL/D5s2bs5VHnk481KBBAz355JPm3+vXr8/LrwOQRxy7Tvj5+RVgSQBr8PLy0j//+U8VK1ZMkvTzzz/nyzi2L730knbt2iVJql27trZt26auXbu6/Uy9evU0b948zZ07VwEBAXleRgBAwZg5c6auXr0qSerXr1+u5BkaGqqXXnrJ/Ds/7vfi4+PVu3dvJSUlSZKeeeYZLViwQGFhYS4/Y7PZ1KdPH+3evVsDBgzI8zICWfHBBx/ov//9ryIiIlymadeunSZMmGD+nZWxhnft2mVOMvbpp58qMDAw22W9Vbz22muaMGGCevXqlaNhOJ5++mlt2rRJ999/v3lfkFH9+vU1ZcoU8+9beRzpXr16mfGGL774Ilt55Pns6vXq1TP/Hx8f7zRNZrOrO5vt8OzZsxo7dqwaNmyogIAABQQEqGnTppowYYJZebty7do1/fDDD3r66afTDeobFBSk22+/XSNHjtSvv/6a6bplnEU4JSVF06dP1z333KMqVarI29s73Qx4mc3amNXJJ9zNrnfjxg19/fXXevDBB1W1alX5+/srODhYd9xxh1588UWPuzxcuXJF7777rpo2barAwEAFBQXpjjvu0FtvvaWzZ896lIenMps10Nl+snr1avXs2dNcx5o1a2rYsGE3rd/169f15ZdfqnXr1ipfvrxKlCihhg0b6r333nM69omjhIQEzZkzR0888YSaNGmi4OBg+fj4qEyZMmratKlefPFFHT161OP1PHr0qJ566inVqlVLxYsXV4UKFdS+fXtNmTLFbBZvX88OHTpkmt/WrVs1cuRI1a9fX8HBwfL391dYWJgeeeQRLV261ONyOWPfxyMjI833xo0bd9O+6Ex0dLRGjx6tpk2bqkyZMvLz81OVKlX0wAMPaPr06Td1AXD13Z4eX1l1+vRpvfLKK6pfv75KliypMmXKqEWLFvrXv/5lnkMy2yezMtNlVtLu379ff/vb39S4ceN0v92DDz6oOXPmKDU1NdM8/vOf/6hPnz6qUaOGSpQoIX9/f1WpUkWNGjVS79699cUXX+jcuXNOP3vx4kW9//77ioiISDcwdo0aNXTXXXfp6ae
"text/plain": [
"<Figure size 640x480 with 1 Axes>"
]
},
"metadata": {},
"output_type": "display_data"
},
{
"name": "stdout",
"output_type": "stream",
"text": [
"y smaller x, 1d fit along y\n",
"\n",
"1d fit initialization\n",
"center = [100.50251256 76.51006711]\n",
"BEC widths: [24 12]\n",
"\n",
"1d init fit values\n",
"Name Value Min Max Stderr Vary Expr Brute_Step\n",
"amp_bec 1.886 0 4.903 None True None None\n",
"amp_th 1.886 0 4.903 None True None None\n",
"deltax 36 0 200 None True None None\n",
"sigma_bec 9.836 0 24 None True None None\n",
"sigma_th 24.86 0 inf None False 0.632*sigma_bec + 0.518*deltax None\n",
"x0_bec 76.51 66.51 86.51 None True None None\n",
"x0_th 76.51 66.51 86.51 None True None None\n",
"1d fitted values\n",
"x0_bec: 76.429, (init = 76.510), bounds = [66.51 : 86.51] \n",
"x0_th: 76.225, (init = 76.510), bounds = [66.51 : 86.51] \n",
"amp_bec: 3.016, (init = 1.886), bounds = [0.00 : 4.90] \n",
"amp_th: 0.716, (init = 1.886), bounds = [0.00 : 4.90] \n",
"sigma_bec: 11.100, (init = 9.836), bounds = [0.00 : 24.00] \n",
"sigma_th: 8.387, (init = 24.864), bounds = [0.00 : inf] \n",
"\n"
]
},
{
"data": {
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAh8AAAHXCAYAAAAcDjiVAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjcuMSwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/bCgiHAAAACXBIWXMAAA9hAAAPYQGoP6dpAACJoUlEQVR4nO3dd3hUVfrA8e+dmcykzaQHCEkIHZFeFEUFC+qKYgPF1VVwEXAXXV17W8GGLuz+1LVgBV1dbKDL2kFAQBSlBFSUJiEJJIS0SZ9MOb8/JjMmppAyJeX9PM88kJlzz33vlDvvnHuKppRSCCGEEEIEiC7YAQghhBCia5HkQwghhBABJcmHEEIIIQJKkg8hhBBCBJQkH0IIIYQIKEk+hBBCCBFQknwIIYQQIqAk+RBCCCFEQEnyIYQQQoiAkuSjC1q/fj2apqFpGvPnz29zfYWFhdx3332MHDkSs9mMTqdD0zSio6O9ZSZOnOjdZ0e2e/du/vjHP9KvXz/Cw8O9x3TJJZf4bZ+d5blrTzzP58SJE4MdimjHfH2uFL8yBDuArs7pdPLTTz+xdetWtm3bxtatW9m5cyeVlZUAPPjgg+36TX/06FHGjRtHRkZGm+pJT0/ngw8+AOCSSy5hxIgRbY7N17788kvOP/98qqqqgh2KXxQXF/Pkk08CMGLECL8mVEKIrk2SjyC74oorWLlyZbDDaLVHH33Um3iMHz+ea665hm7duqFpGiEhIc2uJz09nQULFgCQlpbWLpOPm266yZt4XHvttUycOJGYmBgAevToEczQfKK4uNj7Glx33XWSfAgh/EaSjyBzOp11/o6NjSUuLo59+/YFKaKW+fjjjwGIiYnh888/Jzw8vMFy69evD2BUvpednc33338PwHnnncdrr70W5IiEEP42ceJEZO1V/5DkI8hOOukkTjjhBEaPHs3o0aPp3bs3y5YtY+bMmcEOrVmysrIAGDhwYKOJR2fgOU6AkSNHBjESIYTo+CT5CLJ777032CG0SXV1NQAmkynIkfiXzWbz/r+zH6sQQvibjHbphDZt2sT06dPp2bMnoaGhpKSkcMkll/Dpp5/6pP758+fXG33x5Zdfeu/z3GpfamlsxMayZcvQNK1OS8/MmTPr1ZWWltammLOysrj77rsZNWoUsbGxmEwmevbsyUUXXcSyZcvqXf76bdxnnnmm974FCxbUi6+1ysvLeeyxxxg1ahQWi4WoqCiGDRvGgw8+SH5+frPryczM5JlnnmHatGkMHDiQyMhIjEYjiYmJTJw4kSeeeAKr1drgthkZGWiaRu/evb33vfbaa/WOUdO0eh2Ljx07xssvv8w111zDkCFDsFgshISEEB8fz7hx4/jb3/5Gbm5uq54bXx9nayileOedd5g6dSqpqamEhoYSHR3NsGHD+Otf/3rcy6MNjZbIzMzktttuY9CgQURERBAdHc2pp57Kc889h8PhaFZc77//PpMnT6Zbt26EhoaSlpbGNddcw5YtW4BfP1eaprFs2bIWH/enn37q3f7mm29u1jY333yzd5tPPvmkxfv0aOtrvGfPHiIiItA0jaioKA4ePNho2aKiIlJTU9E0Db1ez4YNG+o83tzRLtu3b2fu3LkMHTrU+xlITExk8ODBnH/++Tz88MMd5lJ6wCjR7ixdulQBClAPPvhgi7a99957laZp3u1/e5s3b55at25dq+tXSqkHH3yw0fpr39atW+fdZsKECd77GzvWpm69evVqcZweS5YsUWFhYU3WP3ToUHXw4MF629aOu6lba+zdu1f17t270TqTk5PVjh07Gn3uPNatW9fka+65JSQkqI0bN9bb/uDBg806RqDOc3TgwAGl1+uPu01ERIRasWJFq54jXx5nbZ6yEyZMaLRMbm6uOuWUU5rcV0hIiFq4cGGTMdf+rH3yyScqOjq60fomTZqkqqqqGq2vurpaTZs2rdHt9Xq9Wrx4cZ3P1dKlS4/31NbjdDq9782YmBhVWVnZZPnKykoVExOjAJWamqqcTmeL96mU717jF1980Vv2lFNOUQ6Ho8FyU6dO9Za77777GozneOfKBx98sFkxX3zxxa15SjotuezSiSxatIjHHnsMcM9jcOWVVzJp0iTCwsLYuXMnL7/8Ms888wyHDx9u036mT5/uHY1y6aWXAnDiiSfyyCOP1Ck3ZMiQ49Z11lln8f7777N27Vr+9a9/Ae5RJWeddVadcq3tT/LCCy8wd+5c798XXXQRkydPJjo6mr1797J06VIOHjzI999/z2mnncaOHTtISEjwln/kkUfIz8/nhx9+4IEHHgDgyiuvZPr06a2Kx6OwsJCzzjqL7OxsAHr37s31119P//79KSwsZOXKlaxZs4bLLruMqKioJuuqqqpCKcWJJ57ImWeeyQknnEBcXBxVVVVkZWXxwQcfsG3bNo4dO8aFF15Ienp6nZakxMRE3n//ffLy8pgzZw4AZ555ZoO/eBMTE73/r66uxul00qdPH84++2yGDBlCQkICLpeLzMxM1qxZw9q1aykvL2f69Ol89dVXjB07ttXPWVuPsyVKS0s544wz2Lt3L+AezXT99ddz4oknUlFRwerVq3n33Xex2+3cc889uFyu415CTU9PZ9GiRSilmDNnDqeccgomk4mtW7eyZMkSysvLWb16NY8++igPPfRQg3XMnj2bd999F4DQ0FBmzJjBKaecgl6vZ+vWrbzyyivcfvvtTJ06tVXH7aHT6bjhhhu49957KSoqYsWKFVx99dWNln/vvfcoKioC4I9//CM6Xesa1X31Gt9www189tlnrFixgq+//poFCxbUe05feeUV3nvvPQDGjRvXqikN/vvf/3pHiIWFhXHVVVcxbtw4YmNjqaqqIjs7m61bt7J69eoW193pBTf3EQ1pTcvH/v37lclkUoAyGo3q448/rlcmJydHDRo0qE423pqWj9o89TT1C1Kpxls+PNr6S60hBw8e9LZ46PV69fbbb9crU1FRoSZPnuzd99SpUxusq62tRb/1xz/+0Vvfeeedp8rLy+uVefLJJ5vVwpKRkaF27drV5P7+85//KJ1OpwA1Y8aMBsvUbgG57rrrjnsMBQUF6quvvmqyzNq1a1VERIQC1MSJE49bZ1N8dZxKHf99O3fuXG+Z0047TRUXF9cr89lnn6nQ0FAFKIPBoNLT0+uVqf2+oaZVYO/evfXKbdmyRRkMBm9LQ0OtH2vWrPHWEx8fr77//vt6ZQ4ePKh69epVZ5+t/Tzl5uaqkJCQZn2+zzjjDO/nLCsrq1X7U8q3r3FhYaFKTk72xlW7pWTPnj3e96XFYlEHDhxosI7jfe495w69Xt/kZ6GyslJt2bKlyePqaiT5aIdak3z85S9/8W4zf/78Rsvt3LmzTlN5Z04+br31Vm+dd955Z6PlrFar6tGjhwKUpmkNfjn4Mvk4evSoMhqNClBxcXGqoKCg0bK1m4Xb+lvh2muvVYAKCwtT1dXV9R5vafLRXH/729+89WZmZvqs3sYc7ziVavp9m5eX503kLRaLysnJaXRfixYt8tZ19dVX13v8t8nHhg0bGq3r6quvbrJc7ST5rbfearSe3+6zLZ+n2pd4GvpcKOX+IveUmTx5cqv31RLNeY2VUmr9+vXeRKVXr16qqKhIVVdXq9GjR3tj/ve//93o9sf73A8cOFABatiwYb44rC5FOpx2Eu+//z4AISEh3HTTTY2WGzZsGOeee26gwgoqz+RtBoOB2267rdFyFouFP/3pTwAopbzPpb989NFH3lFCM2bMIDY2ttGyd9xxh8/2e+qppwJQWVnJrl27fFZvc/cL8O233wZsf609zo8++sg7uum6666je/fujZb905/+hNlsBmDVqlWNdlwG9xDt008/vdHHa19q3L17d53Hqqqq+PzzzwFISkpi2rRpjdYzceJEhg0b1ujjLVH7kuXLL7/cYJna98+ePdsn+z2e5r7GEyZM4O677wbg0KFDzJ07l/vvv59t27YBcPXVV3PNNde0Og7P5eDs7GyfdnTuCiT56ATy8vLIzMwE3Ce4pr7MAM4+++xAhBVUeXl5HDp0CID
"text/plain": [
"<Figure size 640x480 with 1 Axes>"
]
},
"metadata": {},
"output_type": "display_data"
},
{
"name": "stdout",
"output_type": "stream",
"text": [
"\n",
"Init Params\n",
"Name Value Min Max Stderr Vary Expr Brute_Step\n",
"amp_bec 3.665 0 7.277 None True None None\n",
"amp_th 0.8706 0 7.277 None True None None\n",
"atom_number_bec 1006 -inf inf None False amp_bec / 5 * 2 * 3.14159265359 * sigmax_bec * sigmay_bec None\n",
"atom_number_th 281.5 -inf inf None False amp_th * 2 * 3.14159265359 * 1.20206 / 1.643 * sigma_th * sigma_th None\n",
"condensate_fraction 0.7813 -inf inf None False atom_number_bec / (atom_number_bec + atom_number_th) None\n",
"rot_angle -36 -66 -6 None False None None\n",
"sigma_th 8.387 0 200 None True None None\n",
"sigmax_bec 19.67 0 48 None True None None\n",
"sigmay_bec 11.1 0 24 None True None None\n",
"x0_bec 100.5 90.5 110.5 None True None None\n",
"x0_th 100.5 90.5 110.5 None True None None\n",
"y0_bec 76.51 66.51 86.51 None True None None\n",
"y0_th 76.51 66.51 86.51 None True None None\n",
"\n",
"Name Value Min Max Stderr Vary Expr Brute_Step\n",
"amp_bec 3.665 0 7.277 None True None None\n",
"amp_th 0.8706 0 7.277 None True None None\n",
"atom_number_bec 1006 -inf inf None False amp_bec / 5 * 2 * 3.14159265359 * sigmax_bec * sigmay_bec None\n",
"atom_number_th 281.5 -inf inf None False amp_th * 2 * 3.14159265359 * 1.20206 / 1.643 * sigma_th * sigma_th None\n",
"condensate_fraction 0.7813 -inf inf None False atom_number_bec / (atom_number_bec + atom_number_th) None\n",
"rot_angle -36 -66 -6 None False None None\n",
"sigma_th 8.387 0 200 None True None None\n",
"sigmax_bec 19.67 0 48 None True None None\n",
"sigmay_bec 11.1 0 24 None True None None\n",
"x0_bec 100.5 90.5 110.5 None True None None\n",
"x0_th 100.5 90.5 110.5 None True None None\n",
"y0_bec 76.51 66.51 86.51 None True None None\n",
"y0_th 76.51 66.51 86.51 None True None None\n",
"bval first fit\n",
"amp_bec: 3.557, (init = 3.665), bounds = [0.00 : 7.28] \n",
"amp_th: 0.865, (init = 0.871), bounds = [0.00 : 7.28] \n",
"x0_bec: 101.119, (init = 100.503), bounds = [90.50 : 110.50] \n",
"y0_bec: 75.951, (init = 76.510), bounds = [66.51 : 86.51] \n",
"x0_th: 100.508, (init = 100.503), bounds = [90.50 : 110.50] \n",
"y0_th: 76.096, (init = 76.510), bounds = [66.51 : 86.51] \n",
"sigmax_bec: 24.616, (init = 19.672), bounds = [0.00 : 48.00] \n",
"sigmay_bec: 11.485, (init = 11.100), bounds = [0.00 : 24.00] \n",
"sigma_th: 8.771, (init = 8.387), bounds = [0.00 : 200.00] \n",
"rot_angle: -36.000, (init = -36.000), bounds = [-66.00 : -6.00] \n",
"\n",
"amp_bec: 3.557, (init = 3.665), bounds = [0.00 : 7.28] \n",
"amp_th: 0.865, (init = 0.871), bounds = [0.00 : 7.28] \n",
"x0_bec: 101.119, (init = 100.503), bounds = [90.50 : 110.50] \n",
"y0_bec: 75.951, (init = 76.510), bounds = [66.51 : 86.51] \n",
"x0_th: 100.508, (init = 100.503), bounds = [90.50 : 110.50] \n",
"y0_th: 76.096, (init = 76.510), bounds = [66.51 : 86.51] \n",
"sigmax_bec: 24.616, (init = 19.672), bounds = [0.00 : 48.00] \n",
"sigmay_bec: 11.485, (init = 11.100), bounds = [0.00 : 24.00] \n",
"sigma_th: 8.771, (init = 8.387), bounds = [0.00 : 200.00] \n",
"rot_angle: -36.000, (init = -36.000), bounds = [-66.00 : -6.00] \n",
"\n",
"\n",
"Atom numbers:\n",
" N_bec: 89930\n",
" N_th: 21771\n",
" N_ges: 111702\n",
" Cond. frac: 80.51 %\n",
"\n",
"Temperature: 27.44 nK\n",
"\n",
" \n",
"file high_cond_1 \n",
" image: 3\n",
"shape: (200, 150)\n"
]
},
{
"data": {
"image/png": "iVBORw0KGgoAAAANSUhEUgAABTkAAAHXCAYAAAB3dl6BAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjcuMSwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/bCgiHAAAACXBIWXMAAA9hAAAPYQGoP6dpAACDs0lEQVR4nOzdd3wU1f7/8feGFFoKnUAIoSOCFGmCECIoCCooRVC60iyo1wbqFbgqoNdrRUFQqYqAIBcvIr1IFSlSlE4gREOoCTWkzO+PfHd+m7C72fRMeD0fj31Ads6ePbNTzsxnTrEZhmEIAAAAAAAAACzKK78LAAAAAAAAAADZQZATAAAAAAAAgKUR5AQAAAAAAABgaQQ5AQAAAAAAAFgaQU4AAAAAAAAAlkaQEwAAAAAAAIClEeQEAAAAAAAAYGkEOQEAAAAAAABYGkFOAAAAAAAAAJZWIIKcY8eOlc1mk81m07p16/K7OLkmMjLSXM+BAwfmd3FyTLt27cz1yo6wsDDZbDaFhYXlTMGQ4/744w898cQTqlmzpooXL25u927duuV30VBA5NT5AHCU0/VnTtc31F9A3psxY4Z5XpgxY0a28yus1+nOXL58WSEhIbLZbPrPf/6T38WBB3Kynhk4cKC5r0dGRha4/AAUHOvWrTOP77Fjx+ba92zcuFE2m00+Pj7av39/tvLKUpDTvpLOXl5eXgoICFDdunXVt29f/e9//8tWAQEUDOvXr9edd96pr7/+WkePHtW1a9fyu0hAnnN8KOfs5efnpwoVKig8PFxvvfWW/vrrrwzztN+4ZOYVFBTkUXnPnz+vadOm6dFHH1XdunVVpkwZ+fj4qEyZMqpfv74GDBigOXPm6MqVK9n8ZQqXxYsXa+zYsRo7dqwuXryY38VBPvjoo480duxYffTRR/ldFGSD/TjOiSBoYfL2228rOjpaFSpU0IgRI5ymyahuKlmypGrUqKEePXpo7ty5SkpKcvudjjfKmXk9//zzHq3T7t27NW7cOLVr105hYWEqWbKkihYtquDgYLVt21YvvfSStmzZktmf6pa0e/du89jZvXt3fhfnlhEXF6f58+drxIgRatGihXnNVqpUKTVs2FBPPfWUtm/fnu3vGTp0aJpjLDcDV1aWnJysffv2acaMGXr22Wd11113pWnkk5nfbefOnfr88881ePBgNWvW7KZzVPv27TVx4kTFxMTk3goVcHfffbfuueceJSUlaeTIkdnKyzuHymQyDEOXLl3SwYMHdfDgQX3zzTcKDw/XwoULVaZMmZz+OgB55Nlnn9X169clSf3791e7du1UqlQpSVJwcHB+Fg0oMG7cuKHY2FjFxsZqw4YNevfdd/XZZ59pwIABeV6Ot99+Wx999JEuXbp00/Lz58/r/Pnz2r9/v2bNmqWSJUvqxRdf1Msvv6wSJUrkaVkLosWLF2vmzJmSUluoeBpURuHx0Ucf6cSJE6patarHQRYUPOPGjZMkhYeHF/rWmZ6Kjo7Wxx9/LEl69dVXVbx48Szlc+XKFR07dkzHjh3TwoULNWHCBP33v/9VtWrVcrK4GTp48KBeeukllw1rYmJiFBMTo19++UX/+c9/VL9+fU2cOFFdunTJ03JaiT1gLKUGuxs1apS/BboFvPfee3rzzTeVkJBw07KLFy/q4sWL2rNnjyZPnqy+ffvqiy++yNKxu27dOn355Zc5UeRCr1evXlq0aFGO5NW5c2edPn3a6TL7OWrNmjV655139J///EdDhw7Nke+1mrFjx2rNmjVas2aNli9fro4dO2Ypn2wHOX/44Yc0f6ekpOjs2bPaunWr5s6dq+vXr2v9+vXq2rWrfvnlF6ddGO1Pigq7sLAwGYaR38UosOjeUHCdOnVKe/fulSR17NjRvPkH0ivMQ46k9+ijj6p3795p3rt27ZqOHj2qBQsWaM+ePbpy5YoGDRqkChUqqFOnThnm+cUXX6h8+fIZpvPx8XG5LDY2Vg8//LA2b95svle/fn116tRJNWvWVOnSpXXx4kWdOHFCK1as0I4dO3T58mWNGzdOYWFhBTIQQP0JIL2BAwfm6PnqVjnPTJgwQdevX1epUqU0fPhwjz6Tvm5KSUnRxYsXtXPnTs2ZM0dxcXHau3ev7r33Xu3ZsyfD4Mvtt9+ut99+26PvrlGjhstlP/30k/r06aP4+HhJkq+vr+655x6Fh4crODhYRYsW1enTp7V3714tW7ZM0dHR2rdvnx544AHLbWvukwq3Q4cOmQHO6tWrq0OHDmrUqJHKli2rCxcuaPXq1Vq4cKGSk5M1Z84cxcbGatmyZfLy8rxj7rVr1zRkyBAZhqESJUrQiycDycnJaf4uXbq0ypQpo8OHD2cpv7Jly6ply5Zq2LChqlWrpsDAQN24cUNHjhzR4sWLtWvXLl2+fFnDhg2Tt7e3Bg8enBOrkW3t2rXLs/NlmzZtdPfdd2vjxo0aM2ZMloOcMrJAkvlyZ+/evYa/v7+Z9qeffsrK16GACw8P92h/gHVt3rzZ3MajRo3K7+IA+WbMmDHmsTBmzBiX6ZKTk40nnnjCTNu8eXOXaatWrWqmO378eLbKd+3aNaNJkyZmfqGhocbSpUvdfubPP/80evXqZUgypk+fnq3vtwr7b161alWnywcMGJCpbZJRfvnBXqYBAwbkd1EsqSBuU2Se/TgODw/P76IUCOfOnTOKFy9uSDKGDh3qNq2ndVNUVJRRpUoVM+3nn3/uNN3atWtzdHts2rTJ8PHxMfPs3r27ceLECZfpU1JSjHnz5hl16tS55e9ZMqrjpk+fbi735Logs3Umbvbkk08aXbp0MdatW+cyzYYNG4ySJUuav/XXX3+dqe945ZVXDElG5cqVjRdeeMGj69lb2TvvvGOMGjXKWLBggXHs2DHDMNIeG5n53fbt22ekpKS4TTN+/Hgz76CgIOP69evZKb5lTZ061fwdNm3alKU8cnXiofr16+vJJ580/16/fn1ufh2AXOLYdcLPzy8fSwJYg5eXl/7973+rSJEikqRff/01T8axfemll7Rz505JUq1atbR161Z17tzZ7Wfq1q2refPmae7cufL398/1MgIA8sesWbN09epVSVLfvn1zJM+QkBC99NJL5t95cb8XFxennj17KjExUZL0zDPPaMGCBQoNDXX5GZvNpl69emnXrl3q379/rpcRyIz33ntP//vf/xQeHu4yTZs2bTRhwgTz78yMNbxz505zkrFPPvlEAQEBWS7rreK1117ThAkT1KNHj2wPw3H77bdnOCnr6NGjdccdd0hKHaJg06ZN2fpOq+rRo4cZb/j888+zlEeuz65et25d8/9xcXFO02Q0u7qz2Q7Pnj2rsWPHqkGDBvL395e/v7+aNGmiCRMmmJW3K9euXdMPP/ygp59+Os2gvoGBgbr99ts1YsQI/f777xmuW/pZhJOTkzVjxgzde++9qly5sry9vdPMgJfRrI2ZnXzC3ex6N27c0FdffaWHHnpIVapUUdGiRRUUFKQ77rhDL774osddHq5cuaLx48erSZMmCggIUGBgoO644w6NGTNGZ8+e9SgPT2U0a6Cz/WT16tXq3r27uY41atTQ0KFDb1q/69ev64svvlCrVq1Urlw5FS9eXA0aNNDEiROdjn3iKD4+Xt98842eeOIJNW7cWEFBQfLx8VHp0qXVpEkTvfjiizp69KjH63n06FE99dRTqlmzpooVK6by5curbdu2mjp1qtks3r6e7dq1yzC/LVu2aMSIEapXr56CgoJUtGhRhYaG6tFHH9XSpUs9Lpcz9n08IiLCfG/cuHE37YvOREVFadSoUWrSpIlKly4tPz8/Va5cWQ8++KBmzJhxUxcAV9/t6fGVWadPn9Yrr7yievXqqUSJEipdurSaN2+u//znP+Y5JKN9MjMzXWYm7f79+/WPf/xDjRo1SvPbPfTQQ/rmm2+UkpKSYR7//e9/1atXL1WvXl3FixdX0aJFVblyZTVs2FA9e/bU559/rnPnzjn97MWLF/Xuu+8qPDxc5cuXl6+vrwICAlS9enXdddddevrpp7Vs2TKnXRcyml3d2Wy4f/z
"text/plain": [
"<Figure size 640x480 with 1 Axes>"
]
},
"metadata": {},
"output_type": "display_data"
},
{
"name": "stdout",
"output_type": "stream",
"text": [
"y smaller x, 1d fit along y\n",
"\n",
"1d fit initialization\n",
"center = [100.50251256 76.51006711]\n",
"BEC widths: [24 13]\n",
"\n",
"1d init fit values\n",
"Name Value Min Max Stderr Vary Expr Brute_Step\n",
"amp_bec 1.868 0 4.856 None True None None\n",
"amp_th 1.868 0 4.856 None True None None\n",
"deltax 39 0 200 None True None None\n",
"sigma_bec 10.66 0 26 None True None None\n",
"sigma_th 26.94 0 inf None False 0.632*sigma_bec + 0.518*deltax None\n",
"x0_bec 76.51 66.51 86.51 None True None None\n",
"x0_th 76.51 66.51 86.51 None True None None\n",
"1d fitted values\n",
"x0_bec: 76.418, (init = 76.510), bounds = [66.51 : 86.51] \n",
"x0_th: 76.448, (init = 76.510), bounds = [66.51 : 86.51] \n",
"amp_bec: 3.137, (init = 1.868), bounds = [0.00 : 4.86] \n",
"amp_th: 0.548, (init = 1.868), bounds = [0.00 : 4.86] \n",
"sigma_bec: 11.291, (init = 10.656), bounds = [0.00 : 26.00] \n",
"sigma_th: 8.453, (init = 26.936), bounds = [0.00 : inf] \n",
"\n"
]
},
{
"data": {
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAh8AAAHXCAYAAAAcDjiVAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjcuMSwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/bCgiHAAAACXBIWXMAAA9hAAAPYQGoP6dpAACJ30lEQVR4nO3dd3hUVfrA8e+dmkySSQ8tCaFb6IKiqGBBVFwsYFtcBUVgV3RX17q6Imtf3F1xLbg2dFVsoD9siyAgIi7SLaxSJCQBkpA2mbSp5/fHZMaE9GQyk/J+nmeeZOaee+97p9x559xTNKWUQgghhBAiRHThDkAIIYQQ3YskH0IIIYQIKUk+hBBCCBFSknwIIYQQIqQk+RBCCCFESEnyIYQQQoiQkuRDCCGEECElyYcQQgghQkqSDyGEEEKElCQf3dD69evRNA1N03jggQfavL2ioiLuvfdeRo0aRUxMDDqdDk3TiIuLC5SZOHFiYJ+d2e7du7nhhhsYOHAgFoslcEyXXHJJu+2zqzx3HYn/+Zw4cWK4QxEdWLDPleIXhnAH0N15PB7+97//sXXrVrZt28bWrVvZtWsXlZWVACxYsKBDv+nz8vIYN24cmZmZbdrOzp07+eCDDwC45JJLGDlyZJtjC7YvvviC888/n6qqqnCH0i5KSkp48sknARg5cmS7JlRCiO5Nko8wu+KKK1ixYkW4w2i1hx9+OJB4jB8/nmuuuYYePXqgaRpGo7HZ29m5cycLFy4EICMjo0MmHzfffHMg8bj22muZOHEi8fHxAPTq1SucoQVFSUlJ4DW47rrrJPkQQrQbST7CzOPx1LqfkJBAYmIie/fuDVNELfPJJ58AEB8fz2effYbFYqm33Pr160MYVfDl5OTw3XffATB58mReffXVMEckhGhvEydOROZebR+SfITZySefzPHHH89JJ53ESSedRL9+/Vi6dCmzZs0Kd2jNkp2dDcCQIUMaTDy6Av9xAowaNSqMkQghROcnyUeY/elPfwp3CG3idDoBMJvNYY6kfTkcjsD/Xf1YhRCivUlvly5o48aNXHXVVfTp04eIiAjS0tK45JJL+M9//hOU7T/wwAN1el988cUXgcf8t5qXWhrqsbF06VI0TatV0zNr1qw628rIyGhTzNnZ2dx9992MHj2ahIQEzGYzffr04Ve/+hVLly6tc/nr2LjPOuuswGMLFy6sE19rlZeX88gjjzB69GisViuxsbEMHz6cBQsWUFBQ0OztZGVl8fTTT3P55ZczZMgQoqOjMZlMpKSkMHHiRB5//HFsNlu962ZmZqJpGv369Qs89uqrr9Y5Rk3T6jQsPnr0KC+++CLXXHMNQ4cOxWq1YjQaSUpKYty4cdx///3k5ua26rkJ9nG2hlKKd955h+nTp5Oenk5ERARxcXEMHz6c2267rcnLo/X1lsjKyuKPf/wjxx13HFFRUcTFxXHaaafx7LPP4na7mxXX+++/z5QpU+jRowcRERFkZGRwzTXXsHnzZuCXz5WmaSxdurTFx/2f//wnsP4tt9zSrHVuueWWwDqffvppi/fp19bX+KeffiIqKgpN04iNjeXAgQMNli0uLiY9PR1N09Dr9WzYsKHW8ub2dtm+fTvz5s1j2LBhgc9ASkoKJ5xwAueffz4PPvhgp7mUHjJKdDivvPKKAhSgFixY0KJ1//SnPylN0wLrH3ubP3++WrduXau3r5RSCxYsaHD7NW/r1q0LrDNhwoTA4w0da2O3vn37tjhOvyVLlqjIyMhGtz9s2DB14MCBOuvWjLuxW2vs2bNH9evXr8Ftpqamqh07djT43PmtW7eu0dfcf0tOTlZffvllnfUPHDjQrGMEaj1H+/fvV3q9vsl1oqKi1PLly1v1HAXzOGvyl50wYUKDZXJzc9Wpp57a6L6MRqN69NFHG4255mft008/VXFxcQ1ub9KkSaqqqqrB7TmdTnX55Zc3uL5er1dPPPFErc/VK6+80tRTW4fH4wm8N+Pj41VlZWWj5SsrK1V8fLwCVHp6uvJ4PC3ep1LBe43/9a9/Bcqeeuqpyu1211tu+vTpgXL33ntvvfE0da5csGBBs2K++OKLW/OUdFly2aULWbRoEY888gjgG8fgyiuvZNKkSURGRrJr1y5efPFFnn76aQ4dOtSm/Vx11VWB3iiXXnopACeeeCIPPfRQrXJDhw5tcltnn30277//PmvXruWf//wn4OtVcvbZZ9cq19r2JM8//zzz5s0L3P/Vr37FlClTiIuLY8+ePbzyyiscOHCA7777jtNPP50dO3aQnJwcKP/QQw9RUFDA999/z5///GcArrzySq666qpWxeNXVFTE2WefTU5ODgD9+vXj+uuvZ9CgQRQVFbFixQrWrFnDZZddRmxsbKPbqqqqQinFiSeeyFlnncXxxx9PYmIiVVVVZGdn88EHH7Bt2zaOHj3KRRddxM6dO2vVJKWkpPD++++Tn5/P3LlzATjrrLPq/cWbkpIS+N/pdOLxeOjfvz/nnHMOQ4cOJTk5Ga/XS1ZWFmvWrGHt2rWUl5dz1VVX8dVXXzF27NhWP2dtPc6WsNvtnHnmmezZswfw9Wa6/vrrOfHEE6moqGD16tW8++67uFwu7rnnHrxeb5OXUHfu3MmiRYtQSjF37lxOPfVUzGYzW7duZcmSJZSXl7N69Woefvhh/vKXv9S7jTlz5vDuu+8CEBERwcyZMzn11FPR6/Vs3bqVl156idtvv53p06e36rj9dDodN954I3/6058oLi5m+fLlzJgxo8Hy7733HsXFxQDccMMN6HStq1QP1mt84403smrVKpYvX87XX3/NwoUL6zynL730Eu+99x4A48aNa9WQBv/3f/8X6CEWGRnJ1Vdfzbhx40hISKCqqoqcnBy2bt3K6tWrW7ztLi+8uY+oT2tqPvbt26fMZrMClMlkUp988kmdMkeOHFHHHXdcrWy8NTUfNfm309gvSKUarvnwa+svtfocOHAgUOOh1+vV22+/XadMRUWFmjJlSmDf06dPr3dbba0tOtYNN9wQ2N7kyZNVeXl5nTJPPvlks2pYMjMz1bffftvo/t58802l0+kUoGbOnFlvmZo1INddd12Tx1BYWKi++uqrRsusXbtWRUVFKUBNnDixyW02JljHqVTT79t58+YFypx++umqpKSkTplVq1apiIgIBSiDwaB27txZp0zN9w3VtQJ79uypU27z5s3KYDAEahrqq/1Ys2ZNYDtJSUnqu+++q1PmwIEDqm/fvrX22drPU25urjIajc36fJ955pmBz1l2dnar9qdUcF/joqIilZqaGoirZk3JTz/9FHhfWq1WtX///nq30dTn3n/u0Ov1jX4WKisr1ebNmxs9ru5Gko8OqDXJx+9///vAOg888ECD5Xbt2lWrqrwrJx+33nprYJt33nlng+VsNpvq1auXApSmafV+OQQz+cjLy1Mmk0kBKjExURUWFjZYtma1cFt/K1x77bUKUJGRkcrpdNZZ3tLko7nuv//+wHazsrKCtt2GNHWcSjX+vs3Pzw8k8larVR05cqTBfS1atCiwrRkzZtRZfmzysWHDhga3NWPGjEbL1UyS33rrrQa3c+w+2/J5qnmJp77PhVK+L3J/mSlTprR6Xy3RnNdYKaXWr18fSFT69u2riouLldPpVCeddFIg5n//+98Nrt/U537IkCEKUMOHDw/GYXUr0uC0i3j//fcBMBqN3HzzzQ2WGz58OOedd16owgor/+BtBoOBP/7xjw2Ws1qt/O53vwNAKRV4LtvLxx9/HOglNHPmTBISEhose8cddwRtv6eddhoAlZWVfPvtt0HbbnP3C/DNN9+EbH+tPc6PP/440Lvpuuuuo2fPng2W/d3vfkdMTAwAK1eubLDhMvi6aJ9xxhkNLq95qXH37t21llVVVfHZZ58B0Lt3by6//PIGtzNx4kSGDx/e4PKWqHnJ8sUXX6y3TM3H58yZE5T9NqW5r/GECRO4++67ATh48CDz5s3jvvvuY9u2bQDMmDGDa665ptVx+C8H5+TkBLWhc3cgyUcXkJ+fT1ZWFuA7wTX2ZQZwzjnnhCKssMr
"text/plain": [
"<Figure size 640x480 with 1 Axes>"
]
},
"metadata": {},
"output_type": "display_data"
},
{
"name": "stdout",
"output_type": "stream",
"text": [
"\n",
"Init Params\n",
"Name Value Min Max Stderr Vary Expr Brute_Step\n",
"amp_bec 3.835 0 6.902 None True None None\n",
"amp_th 0.67 0 6.902 None True None None\n",
"atom_number_bec 1070 -inf inf None False amp_bec / 5 * 2 * 3.14159265359 * sigmax_bec * sigmay_bec None\n",
"atom_number_th 220 -inf inf None False amp_th * 2 * 3.14159265359 * 1.20206 / 1.643 * sigma_th * sigma_th None\n",
"condensate_fraction 0.8295 -inf inf None False atom_number_bec / (atom_number_bec + atom_number_th) None\n",
"rot_angle -36 -66 -6 None False None None\n",
"sigma_th 8.453 0 200 None True None None\n",
"sigmax_bec 19.67 0 48 None True None None\n",
"sigmay_bec 11.29 0 26 None True None None\n",
"x0_bec 100.5 90.5 110.5 None True None None\n",
"x0_th 100.5 90.5 110.5 None True None None\n",
"y0_bec 76.51 66.51 86.51 None True None None\n",
"y0_th 76.51 66.51 86.51 None True None None\n",
"\n",
"Name Value Min Max Stderr Vary Expr Brute_Step\n",
"amp_bec 3.835 0 6.902 None True None None\n",
"amp_th 0.67 0 6.902 None True None None\n",
"atom_number_bec 1070 -inf inf None False amp_bec / 5 * 2 * 3.14159265359 * sigmax_bec * sigmay_bec None\n",
"atom_number_th 220 -inf inf None False amp_th * 2 * 3.14159265359 * 1.20206 / 1.643 * sigma_th * sigma_th None\n",
"condensate_fraction 0.8295 -inf inf None False atom_number_bec / (atom_number_bec + atom_number_th) None\n",
"rot_angle -36 -66 -6 None False None None\n",
"sigma_th 8.453 0 200 None True None None\n",
"sigmax_bec 19.67 0 48 None True None None\n",
"sigmay_bec 11.29 0 26 None True None None\n",
"x0_bec 100.5 90.5 110.5 None True None None\n",
"x0_th 100.5 90.5 110.5 None True None None\n",
"y0_bec 76.51 66.51 86.51 None True None None\n",
"y0_th 76.51 66.51 86.51 None True None None\n",
"bval first fit\n",
"amp_bec: 3.645, (init = 3.835), bounds = [0.00 : 6.90] \n",
"amp_th: 0.676, (init = 0.670), bounds = [0.00 : 6.90] \n",
"x0_bec: 100.932, (init = 100.503), bounds = [90.50 : 110.50] \n",
"y0_bec: 76.084, (init = 76.510), bounds = [66.51 : 86.51] \n",
"x0_th: 101.313, (init = 100.503), bounds = [90.50 : 110.50] \n",
"y0_th: 75.559, (init = 76.510), bounds = [66.51 : 86.51] \n",
"sigmax_bec: 25.038, (init = 19.672), bounds = [0.00 : 48.00] \n",
"sigmay_bec: 11.686, (init = 11.291), bounds = [0.00 : 26.00] \n",
"sigma_th: 8.691, (init = 8.453), bounds = [0.00 : 200.00] \n",
"rot_angle: -36.000, (init = -36.000), bounds = [-66.00 : -6.00] \n",
"\n",
"amp_bec: 3.645, (init = 3.835), bounds = [0.00 : 6.90] \n",
"amp_th: 0.676, (init = 0.670), bounds = [0.00 : 6.90] \n",
"x0_bec: 100.932, (init = 100.503), bounds = [90.50 : 110.50] \n",
"y0_bec: 76.084, (init = 76.510), bounds = [66.51 : 86.51] \n",
"x0_th: 101.313, (init = 100.503), bounds = [90.50 : 110.50] \n",
"y0_th: 75.559, (init = 76.510), bounds = [66.51 : 86.51] \n",
"sigmax_bec: 25.038, (init = 19.672), bounds = [0.00 : 48.00] \n",
"sigmay_bec: 11.686, (init = 11.291), bounds = [0.00 : 26.00] \n",
"sigma_th: 8.691, (init = 8.453), bounds = [0.00 : 200.00] \n",
"rot_angle: -36.000, (init = -36.000), bounds = [-66.00 : -6.00] \n",
"\n",
"\n",
"Atom numbers:\n",
" N_bec: 95382\n",
" N_th: 16704\n",
" N_ges: 112086\n",
" Cond. frac: 85.10 %\n",
"\n",
"Temperature: 26.94 nK\n",
"\n",
" \n",
"file high_cond_1 \n",
" image: 4\n",
"shape: (200, 150)\n"
]
},
{
"data": {
"image/png": "iVBORw0KGgoAAAANSUhEUgAABTkAAAHXCAYAAAB3dl6BAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjcuMSwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/bCgiHAAAACXBIWXMAAA9hAAAPYQGoP6dpAACDAUlEQVR4nOzdd3wU1f7/8feGNEoKnVBC6IggHRGUEARFQEAFFKUjzYLXa8MOVwXb5VpQBFSqICCIKCJIl44UQZROSABDqAk1JGR+f+Sb+W2S3c0m2ZQJr+fjsQ/IztmzZ3bKmfnMKTbDMAwBAAAAAAAAgEV55XcBAAAAAAAAACAnCHICAAAAAAAAsDSCnAAAAAAAAAAsjSAnAAAAAAAAAEsjyAkAAAAAAADA0ghyAgAAAAAAALA0gpwAAAAAAAAALI0gJwAAAAAAAABLI8gJAAAAAAAAwNIKRJBz9OjRstlsstlsWrNmTX4XJ9dERkaa6zlgwID8Lo7HtG3b1lyvnAgLC5PNZlNYWJhnCgaP++uvvzR48GDVrFlTxYoVM7d79+7d87toKCA8dT4A7Hm6/vR0fUP9BeS9adOmmeeFadOm5Ti/wnqd7silS5dUuXJl2Ww2/fe//83v4sANnqxnBgwYYO7rkZGRBS4/AAVHXtWNR48elY+Pj2w2m5YtW5ajvLIV5ExdSUcvLy8vBQYGqm7duurTp49++umnHBUQQMGwdu1aNW3aVF9//bUOHz6sq1ev5neRgDxn/1DO0cvPz0/ly5dXeHi43nrrLZ08eTLTPFNvXLLyCg4Odqu8586d05QpU/Twww+rbt26Kl26tHx8fFS6dGnVr19f/fv316xZs3T58uUc/jKFy6JFizR69GiNHj1aFy5cyO/iIB989NFHGj16tD766KP8LgpyIPU49kQQtDB5++23deLECZUvX14jRoxwmCazuqlEiRKqUaOGevTooTlz5igpKcnld65ZsybLdZ3NZtO//vUvt9Zp165dGjNmjNq2bauwsDCVKFFC/v7+CgkJUZs2bfT8889r06ZNWf2pbkq7du0yj51du3bld3FuGnFxcZo3b55GjBih22+/3bxmK1mypBo2bKgnnnhC27ZtyzQf+wdAWXm1bds291fSQgzD0IEDBzR79mw999xzatu2rQIDAz0e8Bs7dmya7VDYH7I5U61aNfXr10+S9MwzzygxMTHbeXl7qlCpDMPQxYsXtX//fu3fv1/ffPONwsPDtWDBApUuXdrTXwcgjzz99NO6du2aJKlfv35q27atSpYsKUkKCQnJz6IBBcb169cVGxur2NhYrVu3Tu+9954+++wz9e/fP8/L8fbbb+ujjz7SxYsXMyw/d+6czp07p71792rGjBkqUaKEnnvuOb3wwgsqXrx4npa1IFq0aJGmT58uKaWFirtBZRQeH330kY4dO6aqVau6HWRBwTNmzBhJUnh4+E1745jeiRMn9PHHH0uSXnrpJRUrVixb+Vy+fFlHjhzRkSNHtGDBAo0bN04//PCDqlWr5sniZmr//v16/vnnnTasiYmJUUxMjH777Tf997//Vf369fXuu++qc+fOeVpOK0kNGEspwe5GjRrlb4FuAu+//77eeOMNJSQkZFh24cIFXbhwQbt379bEiRPVp08fTZo0KdvHrjPVq1f3aH5W9/zzz2v8+PG5+h379+/Xf/7zn1z9Dit57bXXNGPGDO3fv19ff/21hg0blq18chzk/P7779P8nZycrDNnzmjz5s2aM2eOrl27prVr16pbt2767bffHHZhTH1SVNiFhYXJMIz8LkaBRfeGguv48ePas2ePJOnee+81b/6B9ArzkCPpPfzww3rkkUfSvHf16lUdPnxY8+fP1+7du3X58mUNHDhQ5cuXV8eOHTPNc9KkSSpXrlym6Xx8fJwui42N1QMPPKCNGzea79WvX18dO3ZUzZo1VapUKV24cEHHjh3T8uXLtX37dl26dEljxoxRWFhYgQwEUH8CSG/AgAEePV/dLOeZcePG6dq1aypZsqSGDx/u1mfS103Jycm6cOGCduzYoVmzZikuLk579uxRhw4dtHv37kyDL7feeqvefvttt767Ro0aTpf9/PPP6t27t+Lj4yVJvr6+ateuncLDwxUSEiJ/f3+dOnVKe/bs0dKlS3XixAn9+eef6tKli+W2NfdJhduBAwfMAGf16tXVvn17NWrUSGXKlNH58+e1cuVKLViwQDdu3NCsWbMUGxurpUuXyssrY8fcdu3aZYjROJKcnKw+ffqYvfMGDhzo2ZWyuBs3bqT5OyAgQFWqVNFff/3lkfwNw9CQIUOUkJCg4sWLF9heVXlZN1arVk2PPPKIZs2apbfffluDBg1yec/jlJENksyXK3v27DECAgLMtD///HN2vg4FXHh4uFv7A6xr48aN5jYeNWpUfhcHyDdvvvmmeSy8+eabTtPduHHDGDx4sJm2RYsWTtNWrVrVTHf06NEcle/q1atGkyZNzPxCQ0ONJUuWuPzM33//bfTq1cuQZEydOjVH328Vqb951apVHS7v379/lrZJZvnlh9Qy9e/fP7+LYkkFcZsi61KP4/Dw8PwuSoFw9uxZo1ixYoYkY+jQoS7Tuls3RUdHG1WqVDHTfv755w7TrV692qPbY8OGDYaPj4+Z50MPPWQcO3bMafrk5GRj7ty5Rp06dW76e5bM6ripU6eay925LshqnYmMHn/8caNz587GmjVrnKZZt26dUaJECfO3/vrrr3P0nUuXLjXzqlWrVo7yKowmTZpkPPvss8Y333xj7Nu3z0hOTk5zHsvp9dXnn39uSDKKFy9ujBkzxmP5Wt3y5cvN3+Kbb77JVh65OvFQ/fr19fjjj5t/r127Nje/DkAuse864efnl48lAazBy8tLH3zwgYoUKSJJ2rp1a56MY/v8889rx44dkqRatWpp8+bN6tSpk8vP1K1bV3PnztWcOXMUEBCQ62UEAOSPGTNm6MqVK5KkPn36eCTPypUr6/nnnzf/zov7vbi4OPXs2dMcs+2pp57S/PnzFRoa6vQzNptNvXr10s6dO81x34CC4v3339dPP/2k8PBwp2nuuusujRs3zvw7p2MNf/311+b/C2Ivnvw2dOhQjR8/Xo8++qjq1Knj0UlVjx8/rlGjRkmS/vOf/7g8d91s2rVrZw6F9/nnn2crj1yfXb1u3brm/+Pi4hymyWx2dUczOp05c0ajR49WgwYNFBAQoICAADVp0kTjxo0zK29nrl69qu+//15PPvlkmkF9g4KCdOutt2rEiBH6448/Ml239LMI37hxQ9OmTVOHDh1UqVIleXt7p5kBL7OZqbI6+YSr2fWuX7+ur776Sl27dlWVKlXk7++v4OBg3XbbbXruuefc7vJw+fJljR07Vk2aNFFgYKCCgoJ022236c0339SZM2fcysNdmc0a6Gg/WblypR566CFzHWvUqKGhQ4dmWL9r165p0qRJatWqlcqWLatixYqpQYMGevfddx2OfWIvPj5e33zzjQYPHqzGjRsrODhYPj4+KlWqlJo0aaLnnntOhw8fdns9Dx8+rCeeeEI1a9ZU0aJFVa5cObVp00aTJ082m8VnZQDoTZs2acSIEapXr56Cg4Pl7++v0NBQPfzww1qyZInb5XIkdR+PiIgw3xszZkyGfdGR6OhojRo1Sk2aNFGpUqXk5+enSpUq6f7779e0adMydAFw9t3uHl9ZderUKb344ouqV6+eihcvrlKlSqlFixb673//a55DMtsnszLTZVbS7t27V//+97/VqFGjNL9d165d9c033yg5OTnTPH744Qf16tVL1atXV7FixeTv769KlSqpYcOG6tmzpz7//HOdPXvW4WcvXLig9957T+Hh4SpXrpx8fX0VGBio6tWr64477tCTTz6ppUuXOuy6kH67pedoNty//vpLw4YNU40aNVS0aFGVLl1ad999t+bMmeNW94jk5GRNnTpVERERKlOmjIoVK6aaNWtq+PDhZpeSzOqZ3FCyZMk03ftSu9PllqioKE2ePFmSVKRIEX3zzTdZGi/3kUce0UMPPZSl7zx+/Lj5uz7
"text/plain": [
"<Figure size 640x480 with 1 Axes>"
]
},
"metadata": {},
"output_type": "display_data"
},
{
"name": "stdout",
"output_type": "stream",
"text": [
"y smaller x, 1d fit along y\n",
"\n",
"1d fit initialization\n",
"center = [100.50251256 75.5033557 ]\n",
"BEC widths: [27 14]\n",
"\n",
"1d init fit values\n",
"Name Value Min Max Stderr Vary Expr Brute_Step\n",
"amp_bec 1.802 0 4.685 None True None None\n",
"amp_th 1.802 0 4.685 None True None None\n",
"deltax 42 0 200 None True None None\n",
"sigma_bec 11.48 0 28 None True None None\n",
"sigma_th 29.01 0 inf None False 0.632*sigma_bec + 0.518*deltax None\n",
"x0_bec 75.5 65.5 85.5 None True None None\n",
"x0_th 75.5 65.5 85.5 None True None None\n",
"1d fitted values\n",
"x0_bec: 75.649, (init = 75.503), bounds = [65.50 : 85.50] \n",
"x0_th: 75.043, (init = 75.503), bounds = [65.50 : 85.50] \n",
"amp_bec: 3.223, (init = 1.802), bounds = [0.00 : 4.68] \n",
"amp_th: 0.355, (init = 1.802), bounds = [0.00 : 4.68] \n",
"sigma_bec: 11.293, (init = 11.475), bounds = [0.00 : 28.00] \n",
"sigma_th: 8.731, (init = 29.008), bounds = [0.00 : inf] \n",
"\n"
]
},
{
"data": {
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAh8AAAHXCAYAAAAcDjiVAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjcuMSwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/bCgiHAAAACXBIWXMAAA9hAAAPYQGoP6dpAACHE0lEQVR4nO3dd3xUVfr48c+dmWRSZ9IDIQmhF+mCoiiggKi44Aq2xQIuArrqrmUta6/oV3d/urqKKyjuqijWZS2LVBFRpKOiNBOSQArpkzbJzJzfH5MZE1NImZbkeb9e8yLMPffe50595txzz6MppRRCCCGEED6i83cAQgghhOheJPkQQgghhE9J8iGEEEIIn5LkQwghhBA+JcmHEEIIIXxKkg8hhBBC+JQkH0IIIYTwKUk+hBBCCOFTknwIIYQQwqck+eiGNm3ahKZpaJrGQw891OHtFRUVce+99zJ69GgiIyPR6XRomkZUVJS7zeTJk9377Mz279/P73//e/r3709YWJj7mC6++GKv7bOrPHaBxPV4Tp482d+hiADm6c9K8QuDvwPo7ux2Oz/++CM7duxg586d7Nixg71791JVVQXAgw8+GNAv+ry8PMaPH09GRkaHtrNnzx4++ugjAC6++GJGjRrV4dg87YsvvuD888+nurra36F4RUlJCc8++ywAo0aN8mpCJYTo3iT58LPLLruMDz74wN9htNvjjz/uTjwmTJjAVVddRWJiIpqmERQU1Ort7Nmzh4cffhiAtLS0gEw+br75Znficc011zB58mSio6MB6Nmzpz9D84iSkhL3c3DttddK8iGE8BpJPvzMbrc3+H9MTAyxsbEcOnTITxG1zaeffgpAdHQ0n3/+OWFhYU2227Rpkw+j8rzs7Gy+++47AKZPn87rr7/u54iEEN42efJkpPaqd0jy4WennXYaQ4YM4dRTT+XUU0+lT58+rFixgvnz5/s7tFbJysoCYNCgQc0mHl2B6zgBRo8e7cdIhBCi85Pkw8/+8pe/+DuEDqmpqQHAaDT6ORLvslqt7r+7+rEKIYS3ydUuXdCWLVu44oor6NWrFyEhIaSkpHDxxRfzv//9zyPbf+ihhxpdffHFF1+473Pd6p9qae6KjRUrVqBpWoOenvnz5zfaVlpaWodizsrK4u6772bMmDHExMRgNBrp1asXv/nNb1ixYkWj01+/jvucc85x3/fwww83iq+9KioqeOKJJxgzZgwmkwmz2cyIESN48MEHKSgoaPV2MjMzeeGFF7j00ksZNGgQERERBAcHk5CQwOTJk3nqqacoLS1tct2MjAw0TaNPnz7u+15//fVGx6hpWqOBxSdOnGDZsmVcddVVDBs2DJPJRFBQEHFxcYwfP54HHniA3Nzcdj02nj7O9lBKsWrVKubMmUNqaiohISFERUUxYsQIbrvttpOeHm3qaonMzExuv/12Bg8eTHh4OFFRUZx55pm8+OKL2Gy2VsX14YcfMmPGDBITEwkJCSEtLY2rrrqKbdu2Ab+8rzRNY8WKFW0+7v/973/u9W+55ZZWrXPLLbe41/nss8/avE+Xjj7HBw4cIDw8HE3TMJvNpKenN9u2uLiY1NRUNE1Dr9ezefPmBstbe7XLrl27WLx4McOHD3e/BxISEhg6dCjnn38+jz76aKc5le4zSgSc1157TQEKUA8++GCb1v3LX/6iNE1zr//r20033aQ2btzY7u0rpdSDDz7Y7Pbr3zZu3OheZ9KkSe77mzvWlm69e/duc5wuS5cuVaGhoS1uf/jw4So9Pb3RuvXjbunWHgcPHlR9+vRpdpvJyclq9+7dzT52Lhs3bmzxOXfd4uPj1Zdfftlo/fT09FYdI9DgMTpy5IjS6/UnXSc8PFy9//777XqMPHmc9bnaTpo0qdk2ubm56owzzmhxX0FBQWrJkiUtxlz/vfbZZ5+pqKioZrc3bdo0VV1d3ez2ampq1KWXXtrs+nq9Xj3zzDMN3levvfbayR7aRux2u/u1GR0draqqqlpsX1VVpaKjoxWgUlNTld1ub/M+lfLcc/zPf/7T3faMM85QNputyXZz5sxxt7v33nubjOdkn5UPPvhgq2KeNWtWex6SLktOu3QhTz/9NE888QTgnMfg8ssvZ9q0aYSGhrJ3716WLVvGCy+8wLFjxzq0nyuuuMJ9Ncpvf/tbAE455RQee+yxBu2GDRt20m2de+65fPjhh2zYsIHnn38ecF5Vcu655zZo197xJC+//DKLFy92//83v/kNM2bMICoqioMHD/Laa6+Rnp7Od999x1lnncXu3buJj493t3/ssccoKCjg+++/5/777wfg8ssv54orrmhXPC5FRUWce+65ZGdnA9CnTx+uu+46BgwYQFFRER988AHr1q3jkksuwWw2t7it6upqlFKccsopnHPOOQwZMoTY2Fiqq6vJysrio48+YufOnZw4cYKLLrqIPXv2NOhJSkhI4MMPPyQ/P59FixYBcM455zT5izchIcH9d01NDXa7nb59+zJlyhSGDRtGfHw8DoeDzMxM1q1bx4YNG6ioqOCKK67gq6++Yty4ce1+zDp6nG1hsViYOHEiBw8eBJxXM1133XWccsopVFZWsnbtWt59911qa2u55557cDgcJz2FumfPHp5++mmUUixatIgzzjgDo9HIjh07WLp0KRUVFaxdu5bHH3+cRx55pMltLFy4kHfffReAkJAQ5s2bxxlnnIFer2fHjh0sX76cO+64gzlz5rTruF10Oh3XX389f/nLXyguLub9999n7ty5zbZ/7733KC4uBuD3v/89Ol37OtU99Rxff/31rFmzhvfff5+vv/6ahx9+uNFjunz5ct577z0Axo8f364pDf7zn/+4rxALDQ3lyiuvZPz48cTExFBdXU12djY7duxg7dq1bd52l+ff3Ec0pT09H4cPH1ZGo1EBKjg4WH366aeN2uTk5KjBgwc3yMbb0/NRn2s7Lf2CVKr5ng+Xjv5Sa0p6erq7x0Ov16t33nmnUZvKyko1Y8YM977nzJnT5LY62lv0a7///e/d25s+fbqqqKho1ObZZ59tVQ9LRkaG2rdvX4v7e+utt5ROp1OAmjdvXpNt6veAXHvttSc9hsLCQvXVV1+12GbDhg0qPDxcAWry5Mkn3WZLPHWcSp38dbt48WJ3m7POOkuVlJQ0arNmzRoVEhKiAGUwGNSePXsatan/uqGuV+DgwYON2m3btk0ZDAZ3T0NTvR/r1q1zbycuLk599913jdqkp6er3r17N9hne99Pubm5KigoqFXv74kTJ7rfZ1lZWe3an1KefY6LiopUcnKyO676PSUHDhxwvy5NJpM6cuRIk9s42fve9dmh1+tbfC9UVVWpbdu2tXhc3Y0kHwGoPcnHH//4R/c6Dz30ULPt9u7d26CrvCsnH7feeqt7m3feeWez7UpLS1XPnj0VoDRNa/LLwZPJR15engoODlaAio2NVYWFhc22rd8t3NHfCtdcc40CVGhoqKqpqWm0vK3JR2s98MAD7u1mZmZ6bLvNOdlxKtXy6zY/P9+dyJtMJpWTk9Psvp5++mn3tubOndto+a+Tj82bNze7rblz57bYrn6S/Pbbbze7nV/vsyPvp/qneJp6Xyjl/CJ3tZkxY0a799UWrXmOlVJq06ZN7kSld+/eqri4WNXU1KhTTz3VHfO///3vZtc/2ft+0KBBClAjRozwxGF1KzLgtIv48MMPAQgKCuLmm29utt2IESM477zzfBWWX7kmbzMYDNx+++3NtjOZTNx4440AKKXcj6W3fPLJJ+6rhObNm0dMTEyzbf/85z97bL9nnnkmAFVVVezbt89j223tfgG+/fZbn+2vvcf5ySefuK9uuvbaa+nRo0ezbW+88UYiIyMBWL16dbMDl8F5ifbZZ5/d7PL6pxr379/fYFl1dTWff/45AElJSVx66aXNbmfy5MmMGDGi2eVtUf+U5bJly5psU//+hQsXemS/J9Pa53jSpEncfffdABw9epTFixdz3333sXPnTgDmzp3LVVdd1e44XKeDs7OzPTrQuTuQ5KMLyM/PJzMzE3B+wLX0ZQYwZcoUX4T
"text/plain": [
"<Figure size 640x480 with 1 Axes>"
]
},
"metadata": {},
"output_type": "display_data"
},
{
"name": "stdout",
"output_type": "stream",
"text": [
"\n",
"Init Params\n",
"Name Value Min Max Stderr Vary Expr Brute_Step\n",
"amp_bec 3.802 0 6.216 None True None None\n",
"amp_th 0.4182 0 6.216 None True None None\n",
"atom_number_bec 1194 -inf inf None False amp_bec / 5 * 2 * 3.14159265359 * sigmax_bec * sigmay_bec None\n",
"atom_number_th 146.6 -inf inf None False amp_th * 2 * 3.14159265359 * 1.20206 / 1.643 * sigma_th * sigma_th None\n",
"condensate_fraction 0.8907 -inf inf None False atom_number_bec / (atom_number_bec + atom_number_th) None\n",
"rot_angle -36 -66 -6 None False None None\n",
"sigma_th 8.731 0 200 None True None None\n",
"sigmax_bec 22.13 0 54 None True None None\n",
"sigmay_bec 11.29 0 28 None True None None\n",
"x0_bec 100.5 90.5 110.5 None True None None\n",
"x0_th 100.5 90.5 110.5 None True None None\n",
"y0_bec 75.5 65.5 85.5 None True None None\n",
"y0_th 75.5 65.5 85.5 None True None None\n",
"\n",
"Name Value Min Max Stderr Vary Expr Brute_Step\n",
"amp_bec 3.802 0 6.216 None True None None\n",
"amp_th 0.4182 0 6.216 None True None None\n",
"atom_number_bec 1194 -inf inf None False amp_bec / 5 * 2 * 3.14159265359 * sigmax_bec * sigmay_bec None\n",
"atom_number_th 146.6 -inf inf None False amp_th * 2 * 3.14159265359 * 1.20206 / 1.643 * sigma_th * sigma_th None\n",
"condensate_fraction 0.8907 -inf inf None False atom_number_bec / (atom_number_bec + atom_number_th) None\n",
"rot_angle -36 -66 -6 None False None None\n",
"sigma_th 8.731 0 200 None True None None\n",
"sigmax_bec 22.13 0 54 None True None None\n",
"sigmay_bec 11.29 0 28 None True None None\n",
"x0_bec 100.5 90.5 110.5 None True None None\n",
"x0_th 100.5 90.5 110.5 None True None None\n",
"y0_bec 75.5 65.5 85.5 None True None None\n",
"y0_th 75.5 65.5 85.5 None True None None\n",
"bval first fit\n",
"amp_bec: 3.919, (init = 3.802), bounds = [0.00 : 6.22] \n",
"amp_th: 0.302, (init = 0.418), bounds = [0.00 : 6.22] \n",
"x0_bec: 100.954, (init = 100.503), bounds = [90.50 : 110.50] \n",
"y0_bec: 75.048, (init = 75.503), bounds = [65.50 : 85.50] \n",
"x0_th: 102.189, (init = 100.503), bounds = [90.50 : 110.50] \n",
"y0_th: 74.351, (init = 75.503), bounds = [65.50 : 85.50] \n",
"sigmax_bec: 24.985, (init = 22.131), bounds = [0.00 : 54.00] \n",
"sigmay_bec: 11.842, (init = 11.293), bounds = [0.00 : 28.00] \n",
"sigma_th: 10.043, (init = 8.731), bounds = [0.00 : 200.00] \n",
"rot_angle: -36.000, (init = -36.000), bounds = [-66.00 : -6.00] \n",
"\n",
"amp_bec: 3.919, (init = 3.802), bounds = [0.00 : 6.22] \n",
"amp_th: 0.302, (init = 0.418), bounds = [0.00 : 6.22] \n",
"x0_bec: 100.954, (init = 100.503), bounds = [90.50 : 110.50] \n",
"y0_bec: 75.048, (init = 75.503), bounds = [65.50 : 85.50] \n",
"x0_th: 102.189, (init = 100.503), bounds = [90.50 : 110.50] \n",
"y0_th: 74.351, (init = 75.503), bounds = [65.50 : 85.50] \n",
"sigmax_bec: 24.985, (init = 22.131), bounds = [0.00 : 54.00] \n",
"sigmay_bec: 11.842, (init = 11.293), bounds = [0.00 : 28.00] \n",
"sigma_th: 10.043, (init = 8.731), bounds = [0.00 : 200.00] \n",
"rot_angle: -36.000, (init = -36.000), bounds = [-66.00 : -6.00] \n",
"\n",
"\n",
"Atom numbers:\n",
" N_bec: 103694\n",
" N_th: 9957\n",
" N_ges: 113652\n",
" Cond. frac: 91.24 %\n",
"\n",
"Temperature: 35.97 nK\n",
"\n",
" \n",
"file high_cond_1 \n",
" image: 5\n",
"shape: (200, 150)\n"
]
},
{
"data": {
"image/png": "iVBORw0KGgoAAAANSUhEUgAABTkAAAHXCAYAAAB3dl6BAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjcuMSwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/bCgiHAAAACXBIWXMAAA9hAAAPYQGoP6dpAACB+ElEQVR4nO3dd3wU1f7/8feGFFoKnVBCgFAFQaQJQgygIiAiTVG6SrF3wQpXBfQq14IiRalKk3LxKtJ7EymCIp1AAEOoCTWkzO+P/DLfDdndbJJNmfB6Ph77gOycPXtmp5yZz5xiMwzDEAAAAAAAAABYlFdeFwAAAAAAAAAAsoMgJwAAAAAAAABLI8gJAAAAAAAAwNIIcgIAAAAAAACwNIKcAAAAAAAAACyNICcAAAAAAAAASyPICQAAAAAAAMDSCHICAAAAAAAAsDSCnAAAAAAAAAAsLV8EOUeMGCGbzSabzaY1a9bkdXFyTGRkpLme/fv3z+vieMw999xjrld2hIaGymazKTQ01DMFg8ft3btXTzzxhMLCwlS0aFFzu3fp0iWvi4Z8wlPnA8Cep+tPT9c31F9A7ps6dap5Xpg6dWq28yuo1+mOXL58WZUqVZLNZtOnn36a18WBGzxZz/Tv39/c1yMjI/NdfgDyD0/Xtc7MnDlTNptNQUFBiomJyVZeWQpypq6ko5eXl5cCAgJUu3Zt9e7dW//73/+yVUAA+cPatWt155136rvvvtPhw4d17dq1vC4SkOvsH8o5evn5+alcuXIKDw/X+++/r1OnTmWYZ+qNS2ZeQUFBbpX3/PnzmjRpkh555BHVrl1bpUqVko+Pj0qVKqV69eqpX79+mjlzpq5cuZLNX6ZgWbRokUaMGKERI0bo4sWLeV0c5IHPPvtMI0aM0GeffZbXRUE2pB7HOXljZkUffPCBTp48qXLlymno0KEO02RUNxUvXlzVq1dX9+7dNWvWLCUmJrr8zjVr1mS6rrPZbHrxxRfdWqddu3Zp5MiRuueeexQaGqrixYurcOHCCg4OVuvWrfXqq69q8+bNmf2pbkm7du0yj51du3bldXFuGbGxsZo7d66GDh2qZs2amddsJUqUUIMGDfT0009r27Ztmcrz0KFDeuONN9SsWTOVLFlSPj4+ZqymT58++uWXX3Jobazvxo0b2r59uyZOnKhBgwbpzjvvlK+vr8cDfvfff3+ac96tWl/16tVLNWvWVGxsrIYPH56tvLw9VCaTYRi6dOmS9u/fr/379+v7779XeHi45s+fr1KlSnn66wDkkueee07Xr1+XJPXt21f33HOPSpQoIUkKDg7Oy6IB+caNGzcUExOjmJgYrVu3Th999JG++uor9evXL9fL8cEHH+izzz7TpUuX0i0/f/68zp8/r7/++kvTp09X8eLF9corr+i1115TsWLFcrWs+dGiRYs0bdo0SSktVNwNKqPg+Oyzz3Ts2DFVqVLF7SAL8p+RI0dKksLDwwt860x3nTx5Up9//rkk6Y033lDRokWzlM+VK1d05MgRHTlyRPPnz9fo0aP13//+V1WrVvVkcTO0f/9+vfrqq04b1kRHRys6Olrr16/Xp59+qnr16mnMmDHq2LFjrpbTSlIDxlJKsLthw4Z5W6BbwMcff6x3331X8fHx6ZZdvHhRFy9e1O7duzV+/Hj17t1bEyZMyPDYHTNmjN59910lJCSked8+VjNz5ky1adNG8+bNU8mSJT26TlZ31113aceOHTn6HdOmTdOyZcty9DusolChQnr77bfVt29fTZ06Va+++qrq1KmTpbyyHeRcuHBhmr+Tk5N19uxZbdmyRbNmzdL169e1du1aPfTQQ1q/fr3DLoypT4oKutDQUBmGkdfFyLfo3pB/nThxQnv27JGU8rQp9eYfuFlBHnLkZo888ogeffTRNO9du3ZNhw8f1rx587R7925duXJFAwYMULly5dS+ffsM85wwYYLKli2bYTofHx+ny2JiYvTwww9r06ZN5nv16tVT+/btFRYWppIlS+rixYs6duyYli1bpu3bt+vy5csaOXKkQkND82UggPoTwM369+/v0fPVrXKeGT16tK5fv64SJUpoyJAhbn3m5ropOTlZFy9e1I4dOzRz5kzFxsZqz549uvfee7V79+4Mgy+33XabPvjgA7e+u3r16k6X/fLLL+rVq5fi4uIkSb6+vmrTpo3Cw8MVHByswoUL6/Tp09qzZ4+WLFmikydP6s8//1SnTp0st625TyrYDhw4YAY4q1Wrpnbt2qlhw4YqXbq0Lly4oJUrV2r+/PlKSkrSzJkzFRMToyVLlsjLy3HH3M8++yxNa7jWrVurY8eOqly5si5cuKCdO3dqxowZio+P16pVq9SxY0dt2LBBhQoVypX1tYKkpKQ0f5cvX15+fn46duyYR/KPiYnRyy+/LEkqVqxYvu1V5em61pXHHntM77zzjo4dO6aRI0dq9uzZWcvIyAJJ5suVPXv2GP7+/mbaX375JStfh3wuPDzcrf0B1rVp0yZzGw8bNiyviwPkmffee888Ft577z2n6ZKSkownnnjCTNu0aVOnaatUqWKmO3r0aLbKd+3aNaNRo0ZmfiEhIcbPP//s8jN///230bNnT0OSMWXKlGx9v1Wk/uZVqlRxuLxfv36Z2iYZ5ZcXUsvUr1+/vC6KJeXHbYrMSz2Ow8PD87oo+cK5c+eMokWLGpKMQYMGuUzrbt0UFRVlVK5c2Uz79ddfO0y3evVqj26PjRs3Gj4+Pmae3bp1M44dO+Y0fXJysjFnzhyjVq1at/w9S0Z13JQpU8zl7lwXZLbORHpPPvmk0bFjR2PNmjVO06xbt84oXry4+Vt/9913DtNduXIlTbpvv/3WYbojR44YFStWNNMtXLjQE6tSYLz00kvGu+++a/z3v/81Tp48aRhG2vuA7F4zp15733HHHUbv3r09lq/Vvfnmm4Yko1ChQi7P6a7k6MRD9erV05NPPmn+vXbt2pz8OgA5xL7rhJ+fXx6WBLAGLy8v/fvf/zafiP/222+5Mo7tq6++anatqVGjhrZs2aIOHTq4/Ezt2rU1Z84czZo1S/7+/jleRgBA3pg+fbquXr0qSerdu7dH8qxUqZJeffVV8+/cuN+LjY1Vjx49zG64zz77rObNm6eQkBCnn7HZbOrZs6d27typvn375ngZgcz4+OOP9b///U/h4eFO07Rq1UqjR482/3Y2duOmTZt0+fJlSVKTJk00cOBAh+mqVq2qYcOGmX+vX78+CyUvuMaOHauRI0eqc+fOqlChgkfzXrx4sebOnSsvLy9NnDiRFrR2Hn/8cUkpLWknTpyYpTxyfHb12rVrm/+PjY11mCaj2dUdzXZ49uxZjRgxQvXr15e/v7/8/f3VqFEjjR492qy8nbl27ZoWLlyoZ555Js2gvoGBgbrttts0dOhQ/fHHHxmu282zCCclJWnq1Km69957VbFiRXl7e6eZAS+jWRszO/mEq9n1bty4oW+//VadO3dW5cqVVbhwYQUFBen222/XK6+84naXhytXrmjUqFFq1KiRAgICFBgYqNtvv13vvfeezp4961Ye7spo1kBH+8nKlSvVrVs3cx2rV6+uQYMGpVu/69eva8KECWrRooXKlCmjokWLqn79+hozZozDsU/sxcXF6fvvv9cTTzyhO+64Q0FBQfLx8VHJkiXVqFEjvfLKKzp8+LDb63n48GE9/fTTCgsLU5EiRVS2bFm1bt1aEydONJvFp67nPffck2F+mzdv1tChQ1W3bl0FBQWpcOHCCgkJ0SOPPKKff/7Z7XI5krqPR0REmO+NHDky3b7oSFRUlIYNG6ZGjRqpZMmS8vPzU8WKFfXggw9q6tSp6boAOPtud4+vzDp9+rRef/111a1bV8WKFVPJkiXVtGlTffrpp+Y5JKN9MjMzXWYm7V9//aWXX35ZDRs2TPPbde7cWd9//72Sk5MzzOO///2vevbsqWrVqqlo0aIqXLiwKlasqAYNGqhHjx76+uuvde7cOYefvXjxoj766COFh4erbNmy8vX1VUBAgKpVq6a77rpLzzzzjJYsWeKwq1dGs6s7mqF
"text/plain": [
"<Figure size 640x480 with 1 Axes>"
]
},
"metadata": {},
"output_type": "display_data"
},
{
"name": "stdout",
"output_type": "stream",
"text": [
"y smaller x, 1d fit along y\n",
"\n",
"1d fit initialization\n",
"center = [99.49748744 75.5033557 ]\n",
"BEC widths: [28 14]\n",
"\n",
"1d init fit values\n",
"Name Value Min Max Stderr Vary Expr Brute_Step\n",
"amp_bec 1.774 0 4.612 None True None None\n",
"amp_th 1.774 0 4.612 None True None None\n",
"deltax 42 0 200 None True None None\n",
"sigma_bec 11.48 0 28 None True None None\n",
"sigma_th 29.01 0 inf None False 0.632*sigma_bec + 0.518*deltax None\n",
"x0_bec 75.5 65.5 85.5 None True None None\n",
"x0_th 75.5 65.5 85.5 None True None None\n",
"1d fitted values\n",
"x0_bec: 75.826, (init = 75.503), bounds = [65.50 : 85.50] \n",
"x0_th: 75.876, (init = 75.503), bounds = [65.50 : 85.50] \n",
"amp_bec: 3.340, (init = 1.774), bounds = [0.00 : 4.61] \n",
"amp_th: 0.193, (init = 1.774), bounds = [0.00 : 4.61] \n",
"sigma_bec: 11.386, (init = 11.475), bounds = [0.00 : 28.00] \n",
"sigma_th: 8.417, (init = 29.008), bounds = [0.00 : inf] \n",
"\n"
]
},
{
"data": {
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAh8AAAHXCAYAAAAcDjiVAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjcuMSwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/bCgiHAAAACXBIWXMAAA9hAAAPYQGoP6dpAACFhUlEQVR4nO3dd3xUVfo/8M+dmUxLMpOeAElI6CC9KIoKCoiKolIsiwVcBNZFv2tZlVVXWV3R1d2f7rqKigu6KjbQxYogICCKdFCUEhKSQHqZTJLJ1PP742auCenJtCSf9+s1SmbOvfe5U+48c+6555GEEAJEREREAaIKdgBERETUvTD5ICIiooBi8kFEREQBxeSDiIiIAorJBxEREQUUkw8iIiIKKCYfREREFFBMPoiIiCigmHwQERFRQDH56Ia2bt0KSZIgSRIef/zxDq+vtLQUDz/8MEaNGoXIyEioVCpIkoSoqCilzaRJk5RtdmZHjhzBb3/7W/Tr1w9Go1HZp2uvvdZv2+wqz10o8T6fkyZNCnYoFMJ8faykX2mCHUB353a78fPPP2PPnj3Yu3cv9uzZg4MHD8JmswEAHnvssZB+0xcUFGD8+PHIysrq0HoOHDiAjz/+GABw7bXXYuTIkR2Ozde++eYbXH755aipqQl2KH5RXl6O559/HgAwcuRIvyZURNS9MfkIsuuvvx7r1q0Ldhjt9te//lVJPCZMmICbb74ZiYmJkCQJYWFhrV7PgQMHsGzZMgBAWlpaSCYfd911l5J43HrrrZg0aRKio6MBAD169AhmaD5RXl6uvAa33XYbkw8i8hsmH0Hmdrvr/R0TE4PY2FgcP348SBG1zeeffw4AiI6OxldffQWj0dhou61btwYwKt/Lzc3F4cOHAQDTpk3DG2+8EeSIiMjfJk2aBNZe9Q8mH0F27rnnYvDgwRgzZgzGjBmD9PR0rF69GvPnzw92aK2Sk5MDABg4cGCTiUdX4N1PABg1alQQIyEi6vyYfATZn/70p2CH0CEOhwMAoNPpghyJf9ntduXfXX1fiYj8jVe7dEE7duzAjTfeiF69ekGv1yMlJQXXXnstvvzyS5+s//HHH29w9cU333yj3Oe91T3V0tQVG6tXr4YkSfV6eubPn99gXWlpaR2KOScnBw899BBGjx6NmJgY6HQ69OrVC1dffTVWr17d4PTX2XFfcsklyn3Lli1rEF97VVVV4amnnsLo0aNhMplgNpsxfPhwPPbYYyguLm71erKzs/Hiiy9izpw5GDhwICIiIqDVapGQkIBJkybhmWeegcViaXTZrKwsSJKE9PR05b433nijwT5KktRgYHFRURFWrlyJm2++GUOHDoXJZEJYWBji4uIwfvx4/PnPf0Z+fn67nhtf72d7CCHw/vvvY/bs2UhNTYVer0dUVBSGDx+Oe++9t8XTo41dLZGdnY377rsPgwYNQnh4OKKionDBBRfgpZdegsvlalVcH330EaZPn47ExETo9XqkpaXh5ptvxq5duwD8+rmSJAmrV69u835/+eWXyvJ33313q5a5++67lWW++OKLNm/Tq6Ov8dGjRxEeHg5JkmA2m5GZmdlk27KyMqSmpkKSJKjVamzbtq3e46292mXfvn1YvHgxhg0bpnwGEhISMGTIEFx++eV44oknOs2p9IARFHJWrVolAAgA4rHHHmvTsn/605+EJEnK8mfflixZIrZs2dLu9QshxGOPPdbk+uvetmzZoiwzceJE5f6m9rW5W+/evdscp9eKFSuEwWBodv3Dhg0TmZmZDZatG3dzt/Y4duyYSE9Pb3KdycnJYv/+/U0+d15btmxp9jX33uLj48X27dsbLJ+ZmdmqfQRQ7znKyMgQarW6xWXCw8PF2rVr2/Uc+XI/6/K2nThxYpNt8vPzxfnnn9/stsLCwsTy5cubjbnuZ+2LL74QUVFRTa5v6tSpoqampsn1ORwOMWfOnCaXV6vV4rnnnqv3uVq1alVLT20DbrdbeW9GR0cLm83WbHubzSaio6MFAJGamircbnebtymE717jV199VWl7/vnnC5fL1Wi72bNnK+0efvjhRuNp6Vj52GOPtSrma665pj1PSZfF0y5dyLPPPounnnoKgDyPwQ033ICpU6fCYDDg4MGDWLlyJV588UWcPn26Q9u58cYblatRrrvuOgDAOeecgyeffLJeu6FDh7a4rksvvRQfffQRNm/ejH/9618A5KtKLr300nrt2jue5JVXXsHixYuVv6+++mpMnz4dUVFROHbsGFatWoXMzEwcPnwYF154Ifbv34/4+Hil/ZNPPoni4mL8+OOPePTRRwEAN9xwA2688cZ2xeNVWlqKSy+9FLm5uQCA9PR03H777ejfvz9KS0uxbt06bNq0CTNnzoTZbG52XTU1NRBC4JxzzsEll1yCwYMHIzY2FjU1NcjJycHHH3+MvXv3oqioCFdddRUOHDhQrycpISEBH330EQoLC7Fo0SIAwCWXXNLoL96EhATl3w6HA263G3369MHkyZMxdOhQxMfHw+PxIDs7G5s2bcLmzZtRVVWFG2+8Ed9++y3GjRvX7ueso/vZFlarFRdffDGOHTsGQL6a6fbbb8c555yD6upqbNy4ER988AGcTieWLl0Kj8fT4inUAwcO4Nlnn4UQAosWLcL5558PnU6HPXv2YMWKFaiqqsLGjRvx17/+FX/5y18aXcfChQvxwQcfAAD0ej3mzZuH888/H2q1Gnv27MHrr7+O+++/H7Nnz27XfnupVCrccccd+NOf/oSysjKsXbsWc+fObbL9hx9+iLKyMgDAb3/7W6hU7etU99VrfMcdd2DDhg1Yu3YtvvvuOyxbtqzBc/r666/jww8/BACMHz++XVMa/O9//1OuEDMYDLjpppswfvx4xMTEoKamBrm5udizZw82btzY5nV3ecHNfagx7en5OHHihNDpdAKA0Gq14vPPP2/QJi8vTwwaNKheNt6eno+6vOtp7hekEE33fHh19JdaYzIzM5UeD7VaLd57770Gbaqrq8X06dOVbc+ePbvRdXW0t+hsv/3tb5X1TZs2TVRVVTVo8/zzz7eqhyUrK0scOnSo2e298847QqVSCQBi3rx5jbap2wNy2223tbgPJSUl4ttvv222zebNm0V4eLgAICZNmtTiOpvjq/0UouX37eLFi5U2F154oSgvL2/QZsOGDUKv1wsAQqPRiAMHDjRoU/d9g9pegWPHjjVot2vXLqHRaJSehsZ6PzZt2qSsJy4uThw+fLhBm8zMTNG7d+9622zv5yk/P1+EhYW16vN98cUXK5+znJycdm1PCN++xqWlpSI5OVmJq25PydGjR5X3pclkEhkZGY2uo6XPvffYoVarm/0s2Gw2sWvXrmb3q7th8hGC2pN8/N///Z+yzOOPP95ku4MHD9brKu/Kycc999yjrPOBBx5osp3FYhE9evQQAIQkSY1+Ofgy+SgoKBBarVYAELGxsaKkpKTJtnW7hTv6W+HWW28VAITBYBAOh6PB421NPlrrz3/+s7Le7Oxsn623KS3tpxDNv28LCwuVRN5kMom8vLwmt/Xss88q65o7d26Dx89OPrZt29bkuubOndtsu7pJ8rvvvtvkes7eZkc+T3VP8TT2uRBC/iL3tpk+fXq7t9UWrXmNhRBi69atSqLSu3dvUVZWJhwOhxgzZowS83//+98ml2/pcz9w4EABQAwfPtwXu9WtcMBpF/HRRx8BAMLCwnDXXXc12W748OG47LLLAhVWUHknb9NoNLjvvvuabGcymXDnnXcCAIQQynPpL5999plyldC8efMQExPTZNs//vGPPtvuBRdcAACw2Ww4dOiQz9bb2u0CwA8//BCw7bV3Pz/77DPl6qbbbrsNSUlJTba98847ERkZCQBYv359kwOXAfkS7YsuuqjJx+ueajxy5Ei9x2pqavDVV18BAHr27Ik5c+Y0uZ5JkyZh+PDhTT7eFnVPWa5cubLRNnXvX7hwoU+225LWvsYTJ07EQw89BAA4deoUFi9ejEceeQR79+4FAMydOxc333xzu+Pwng7Ozc316UDn7oDJRxdQWFiI7Ox
"text/plain": [
"<Figure size 640x480 with 1 Axes>"
]
},
"metadata": {},
"output_type": "display_data"
},
{
"name": "stdout",
"output_type": "stream",
"text": [
"\n",
"Init Params\n",
"Name Value Min Max Stderr Vary Expr Brute_Step\n",
"amp_bec 4.061 0 6.546 None True None None\n",
"amp_th 0.2343 0 6.546 None True None None\n",
"atom_number_bec 1334 -inf inf None False amp_bec / 5 * 2 * 3.14159265359 * sigmax_bec * sigmay_bec None\n",
"atom_number_th 76.31 -inf inf None False amp_th * 2 * 3.14159265359 * 1.20206 / 1.643 * sigma_th * sigma_th None\n",
"condensate_fraction 0.9459 -inf inf None False atom_number_bec / (atom_number_bec + atom_number_th) None\n",
"rot_angle -36 -66 -6 None False None None\n",
"sigma_th 8.417 0 200 None True None None\n",
"sigmax_bec 22.95 0 56 None True None None\n",
"sigmay_bec 11.39 0 28 None True None None\n",
"x0_bec 99.5 89.5 109.5 None True None None\n",
"x0_th 99.5 89.5 109.5 None True None None\n",
"y0_bec 75.5 65.5 85.5 None True None None\n",
"y0_th 75.5 65.5 85.5 None True None None\n",
"\n",
"Name Value Min Max Stderr Vary Expr Brute_Step\n",
"amp_bec 4.061 0 6.546 None True None None\n",
"amp_th 0.2343 0 6.546 None True None None\n",
"atom_number_bec 1334 -inf inf None False amp_bec / 5 * 2 * 3.14159265359 * sigmax_bec * sigmay_bec None\n",
"atom_number_th 76.31 -inf inf None False amp_th * 2 * 3.14159265359 * 1.20206 / 1.643 * sigma_th * sigma_th None\n",
"condensate_fraction 0.9459 -inf inf None False atom_number_bec / (atom_number_bec + atom_number_th) None\n",
"rot_angle -36 -66 -6 None False None None\n",
"sigma_th 8.417 0 200 None True None None\n",
"sigmax_bec 22.95 0 56 None True None None\n",
"sigmay_bec 11.39 0 28 None True None None\n",
"x0_bec 99.5 89.5 109.5 None True None None\n",
"x0_th 99.5 89.5 109.5 None True None None\n",
"y0_bec 75.5 65.5 85.5 None True None None\n",
"y0_th 75.5 65.5 85.5 None True None None\n",
"bval first fit\n",
"amp_bec: 3.959, (init = 4.061), bounds = [0.00 : 6.55] \n",
"amp_th: 0.271, (init = 0.234), bounds = [0.00 : 6.55] \n",
"x0_bec: 99.954, (init = 99.497), bounds = [89.50 : 109.50] \n",
"y0_bec: 76.048, (init = 75.503), bounds = [65.50 : 85.50] \n",
"x0_th: 101.409, (init = 99.497), bounds = [89.50 : 109.50] \n",
"y0_th: 75.044, (init = 75.503), bounds = [65.50 : 85.50] \n",
"sigmax_bec: 25.489, (init = 22.951), bounds = [0.00 : 56.00] \n",
"sigmay_bec: 11.941, (init = 11.386), bounds = [0.00 : 28.00] \n",
"sigma_th: 8.062, (init = 8.417), bounds = [0.00 : 200.00] \n",
"rot_angle: -36.000, (init = -36.000), bounds = [-66.00 : -6.00] \n",
"\n",
"amp_bec: 3.959, (init = 4.061), bounds = [0.00 : 6.55] \n",
"amp_th: 0.271, (init = 0.234), bounds = [0.00 : 6.55] \n",
"x0_bec: 99.954, (init = 99.497), bounds = [89.50 : 109.50] \n",
"y0_bec: 76.048, (init = 75.503), bounds = [65.50 : 85.50] \n",
"x0_th: 101.409, (init = 99.497), bounds = [89.50 : 109.50] \n",
"y0_th: 75.044, (init = 75.503), bounds = [65.50 : 85.50] \n",
"sigmax_bec: 25.489, (init = 22.951), bounds = [0.00 : 56.00] \n",
"sigmay_bec: 11.941, (init = 11.386), bounds = [0.00 : 28.00] \n",
"sigma_th: 8.062, (init = 8.417), bounds = [0.00 : 200.00] \n",
"rot_angle: -36.000, (init = -36.000), bounds = [-66.00 : -6.00] \n",
"\n",
"\n",
"Atom numbers:\n",
" N_bec: 107761\n",
" N_th: 5769\n",
" N_ges: 113530\n",
" Cond. frac: 94.92 %\n",
"\n",
"Temperature: 23.18 nK\n",
"\n",
" \n",
"file high_cond_1 \n",
" image: 6\n",
"shape: (200, 150)\n"
]
},
{
"data": {
"image/png": "iVBORw0KGgoAAAANSUhEUgAABTkAAAHXCAYAAAB3dl6BAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjcuMSwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/bCgiHAAAACXBIWXMAAA9hAAAPYQGoP6dpAACD70lEQVR4nOzdd3wU1f7/8feGFFoKnUAIoSOCFGmCECIoCCgoRVC6SrGgXhuoV+GKYLteCwqCSlWkysWLSC9SRYoUpRMI0RBqQg0p8/sj353fhuxuNj0TXs/HYx+EnbNnz+yUM/OZU2yGYRgCAAAAAAAAAIvyyu8CAAAAAAAAAEB2EOQEAAAAAAAAYGkEOQEAAAAAAABYGkFOAAAAAAAAAJZGkBMAAAAAAACApRHkBAAAAAAAAGBpBDkBAAAAAAAAWBpBTgAAAAAAAACWRpATAAAAAAAAgKUViCDnmDFjZLPZZLPZtG7duvwuTq6JjIw013PQoEH5XZwc065dO3O9siMsLEw2m01hYWE5UzDkuD/++EOPP/64atasqeLFi5vbvXv37vldNBQQOXU+ABzldP2Z0/UN9ReQ96ZPn26eF6ZPn57t/Arrdbozly9fVkhIiGw2m/7973/nd3HggZysZwYNGmTu65GRkQUuPwAFx7p168zje8yYMbn2PRs3bpTNZpOPj4/279+frbyyFOS0r6Szl5eXlwICAlS3bl3169dP//vf/7JVQAAFw/r163XnnXfqm2++0dGjR3Xt2rX8LhKQ5xwfyjl7+fn5qUKFCgoPD9fbb7+tv/76K8M87TcumXkFBQV5VN7z589r6tSpeuSRR1S3bl2VKVNGPj4+KlOmjOrXr6+BAwdq9uzZunLlSjZ/mcJl8eLFGjNmjMaMGaOLFy/md3GQDz7++GONGTNGH3/8cX4XBdlgP45zIghamIwbN07R0dGqUKGCRowY4TRNRnVTyZIlVaNGDfXs2VNz5sxRUlKS2+90vFHOzOv555/3aJ12796tsWPHql27dgoLC1PJkiVVtGhRBQcHq23btnrppZe0ZcuWzP5Ut6Tdu3ebx87u3bvzuzi3jLi4OM2bN08jRoxQixYtzGu2UqVKqWHDhnrqqae0ffv2TOV57NgxvfTSS2rcuLFKlSqlokWLqmrVqnrooYe0YMGCXFqTwiE5OVn79u3T9OnT9eyzz+quu+5K08gnpwJ+Q4cOTXPOy81AYkF2991365577lFSUpJGjhyZrby8c6hMJsMwdOnSJR08eFAHDx7Ut99+q/DwcC1cuFBlypTJ6a8DkEeeffZZXb9+XZI0YMAAtWvXTqVKlZIkBQcH52fRgALjxo0bio2NVWxsrDZs2KD33ntPn3/+uQYOHJjn5Rg3bpw+/vhjXbp0Kd3y8+fP6/z589q/f79mzpypkiVL6sUXX9TLL7+sEiVK5GlZC6LFixdrxowZklJbqHgaVEbh8fHHH+vEiROqWrWqx0EWFDxjx46VJIWHhxf61pmeio6O1ieffCJJevXVV1W8ePEs5XPlyhUdO3ZMx44d08KFCzVhwgT997//VbVq1XKyuBk6ePCgXnrpJZcNa2JiYhQTE6NffvlF//73v1W/fn29++676tKlS56W00rsAWMpNdjdqFGj/C3QLeD999/Xm2++qYSEhHTLLl68qIsXL2rPnj2aNGmS+vXrpy+//DLDY/e9997TG2+8ke4BxMmTJ3Xy5EktXrxY99xzj+bNm0ecxonevXtr0aJFufod69at01dffZWr32ElY8aM0Zo1a7RmzRotX75cHTt2zFI+2Q5y/vDDD2n+n5KSorNnz2rr1q2aM2eOrl+/rvXr16tbt2765ZdfnHZhtD8pKuzCwsJkGEZ+F6PAontDwXXq1Cnt3btXktSxY0fz5h+4WWEecuRmjzzyiPr06ZPmvWvXruno0aOaP3++9uzZoytXrmjw4MGqUKGCOnXqlGGeX375pcqXL59hOh8fH5fLYmNj9dBDD2nz5s3me/Xr11enTp1Us2ZNlS5dWhcvXtSJEye0YsUK7dixQ5cvX9bYsWMVFhZWIAMB1J8AbjZo0KAcPV/dKueZCRMm6Pr16ypVqpSGDx/u0WdurptSUlJ08eJF7dy5U7Nnz1ZcXJz27t2re++9V3v27Mkw+HL77bdr3LhxHn13jRo1XC776aef1LdvX8XHx0uSfH19dc899yg8PFzBwcEqWrSoTp8+rb1792rZsmWKjo7Wvn371LVrV8tta+6TCrdDhw6ZAc7q1aurQ4cOatSokcqWLasLFy5o9erVWrhwoZKTkzV79mzFxsZq2bJl8vJy3jF33Lhx+uc//ykptRdujx49dN999ykwMFBHjx7VjBkzdPDgQa1Zs0ZdunTR2rVrVaxYsTxbXytITk5O8//SpUurTJkyOnz4cI7kf+3aNT355JMyDEMlSpQosL2q2rVrl2fnyzZt2ujuu+/Wxo0b9dZbb2U5yCkjCySZL3f27t1r+Pv7m2l/+umnrHwdCrjw8HCP9gdY1+bNm81tPGrUqPwuDpBv3nrrLfNYeOutt1ymS05ONh5//HEzbfPmzV2mrVq1qpnu+PHj2SrftWvXjCZNmpj5hYaGGkuXLnX7mT///NPo3bu3IcmYNm1atr7fKuy/edWqVZ0uHzhwYKa2SUb55Qd7mQYOHJjfRbGkgrhNkXn24zg8PDy/i1IgnDt3zihevLghyRg6dKjbtJ7WTVFRUUaVKlXMtF988YXTdGvXrs3R7bFp0ybDx8fHzLNHjx7GiRMnXKZPSUkx5s6da9SpU+eWv2fJqI6bNm2audyT64LM1plI74knnjC6dOlirFu3zmWaDRs2GCVLljR/62+++cZpun379hleXl6GJMPb29v48ccf06W5ceOG0atXLzOvsWPH5ti6FBbvvPOOMWrUKGP+/PnGsWPHDMNIe2y4uw/wxCuvvGJIMipXrmy88MILOZav1U2ZMsX8LTZt2pSlPHJ14qH69evriSeeMP+/fv363Pw6ALnEseuEn59fPpYEsAYvLy998MEHKlKkiCTp119/zZNxbF966SXt3LlTklSrVi1t3bpVnTt3dvuZunXrau7cuZozZ478/f1zvYwAgPwxc+ZMXb16VZLUr1+/HMkzJCREL730kvn/vLjfi4uLU69evZSYmChJeuaZZzR//nyFhoa6/IzNZlPv3r21a9cuDRgwINfLCGTG+++/r//9738KDw93maZNmzaaMGGC+X9XYw1PnDhRKSkpkqTnnntOXbt2TZfGx8dH06ZNU7ly5czvt7eIRqrXXntNEyZMUM+ePXN8GI6dO3eak759+umnCggIyNH8raxnz55mvOGLL77IUh65Prt63bp1zb/j4uKcpslodnVnsx2ePXtWY8aMUYMGDeTv7y9/f381adJEEyZMMCtvV65du6YffvhBTz/9dJpBfQMDA3X77bdrxIgR+v333zNct5tnEU5OTtb06dN17733qnLlyvL29k4zA15GszZmdvIJd7Pr3bhxQ19//bUefPBBValSRUWLFlVQUJDuuOMOvfjiix53ebhy5YrGjx+vJk2aKCAgQIGBgbrjjjv01ltv6ezZsx7l4amMZg10tp+sXr1aPXr0MNexRo0aGjp0aLr1u379ur788ku1atVK5cqVU/HixdWgQQO9++67Tsc+cRQfH69vv/1Wjz/+uBo3bqygoCD5+PiodOnSatKkiV588UUdPXrU4/U8evSonnrqKdWsWVPFihVT+fLl1bZtW02ZMsVsFm9fz3bt2mWY35YtWzRixAjVq1dPQUFBKlq0qEJDQ/XII49o6dKlHpfLGfs+HhERYb43duzYdPuiM1FRURo1apSaNGmi0qVLy8/PT5UrV9YDDzyg6dOnp+sC4Oq7PT2+Muv06dN65ZVXVK9ePZUoUUKlS5dW8+bN9e9//9s8h2S0T2ZmpsvMpN2/f7/+8Y9/qFGjRml+uwcffFDffvuteeHizn//+1/17t1b1atXV/HixVW0aFFVrlxZDRs2VK9evfTFF1/o3LlzTj978eJFvffeewoPD1f58uXl6+urgIAAVa9eXXf
"text/plain": [
"<Figure size 640x480 with 1 Axes>"
]
},
"metadata": {},
"output_type": "display_data"
},
{
"name": "stdout",
"output_type": "stream",
"text": [
"y smaller x, 1d fit along y\n",
"\n",
"1d fit initialization\n",
"center = [100.50251256 76.51006711]\n",
"BEC widths: [29 14]\n",
"\n",
"1d init fit values\n",
"Name Value Min Max Stderr Vary Expr Brute_Step\n",
"amp_bec 1.773 0 4.609 None True None None\n",
"amp_th 1.773 0 4.609 None True None None\n",
"deltax 42 0 200 None True None None\n",
"sigma_bec 11.48 0 28 None True None None\n",
"sigma_th 29.01 0 inf None False 0.632*sigma_bec + 0.518*deltax None\n",
"x0_bec 76.51 66.51 86.51 None True None None\n",
"x0_th 76.51 66.51 86.51 None True None None\n",
"1d fitted values\n",
"x0_bec: 76.407, (init = 76.510), bounds = [66.51 : 86.51] \n",
"x0_th: 75.691, (init = 76.510), bounds = [66.51 : 86.51] \n",
"amp_bec: 3.376, (init = 1.773), bounds = [0.00 : 4.61] \n",
"amp_th: 0.106, (init = 1.773), bounds = [0.00 : 4.61] \n",
"sigma_bec: 11.466, (init = 11.475), bounds = [0.00 : 28.00] \n",
"sigma_th: 7.247, (init = 29.008), bounds = [0.00 : inf] \n",
"\n"
]
},
{
"data": {
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAh8AAAHXCAYAAAAcDjiVAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjcuMSwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/bCgiHAAAACXBIWXMAAA9hAAAPYQGoP6dpAACGHklEQVR4nO3deXwU5f0H8M/sbnazmz1yJ4Rc3KDchyIooIKoWDwAj2IFLAJatNVaj3pXK1psf9paxQKCrRZFRUs9ipwCgsgNgiJgQg5ykGt3k2w2ezy/PyY7JuRO9sjxeb9e+0p255mZ7+wx+91nnkMSQggQERERBYkq1AEQERFR98Lkg4iIiIKKyQcREREFFZMPIiIiCiomH0RERBRUTD6IiIgoqJh8EBERUVAx+SAiIqKgYvJBREREQcXkoxvatm0bJEmCJEl4+umn2729kpISPPbYYxgxYgRMJhNUKhUkSUJkZKRSZtKkSco+O7Pjx4/jl7/8Jfr27QuDwaAc0w033BCwfXaV564j8T2fkyZNCnUo1IH5+1xJP9GEOoDuzuPx4LvvvsO+ffuwf/9+7Nu3D4cPH4bD4QAAPPXUUx36TV9QUICxY8ciMzOzXds5dOgQPv74YwDADTfcgOHDh7c7Nn/78ssvcfXVV6OqqirUoQREWVkZXn75ZQDA8OHDA5pQEVH3xuQjxG6++WasW7cu1GG02R//+Ecl8Rg/fjxuv/12JCQkQJIkhIWFtXg7hw4dwjPPPAMASE9P75DJx7333qskHnfccQcmTZqEqKgoAECPHj1CGZpflJWVKa/BnDlzmHwQUcAw+Qgxj8dT5350dDRiYmJw8uTJEEXUOp999hkAICoqCl988QUMBkOD5bZt2xbEqPwvJycHR48eBQBMnToVb731VogjIqJAmzRpEjj3amAw+Qixiy66CIMGDcKoUaMwatQo9OrVC6tXr8a8efNCHVqLZGdnAwAGDBjQaOLRFfiOEwBGjBgRwkiIiDo/Jh8h9vvf/z7UIbRLdXU1AECn04U4ksByOp3K/139WImIAo29XbqgnTt34tZbb0XPnj0RHh6OlJQU3HDDDfjf//7nl+0//fTT9XpffPnll8pjvlvtSy2N9dhYvXo1JEmqU9Mzb968ettKT09vV8zZ2dl45JFHMHLkSERHR0On06Fnz5742c9+htWrV9e7/HV+3Jdffrny2DPPPFMvvraqqKjA888/j5EjR8JsNsNisWDo0KF46qmnUFRU1OLtZGVl4dVXX8WsWbMwYMAAGI1GaLVaxMfHY9KkSXjxxRdhtVobXDczMxOSJKFXr17KY2+99Va9Y5QkqV7D4nPnzmHFihW4/fbbMXjwYJjNZoSFhSE2NhZjx47Fk08+ifz8/DY9N/4+zrYQQmDt2rWYOXMmUlNTER4ejsjISAwdOhQPPPBAs5dHG+otkZWVhd/+9rcYOHAgIiIiEBkZiXHjxuG1116D2+1uUVwfffQRpk2bhoSEBISHhyM9PR2333479uzZA+Cnz5UkSVi9enWrj/t///ufsv59993XonXuu+8+ZZ3PP/+81fv0ae9rfOLECURERECSJFgsFmRkZDRatrS0FKmpqZAkCWq1Gtu3b6+zvKW9XQ4cOIBFixZhyJAhymcgPj4eF1xwAa6++mo8++yzneZSetAI6nBWrVolAAgA4qmnnmrVur///e+FJEnK+uffFi9eLLZu3drm7QshxFNPPdXo9mvftm7dqqwzceJE5fHGjrWpW1paWqvj9Fm2bJnQ6/VNbn/IkCEiIyOj3rq1427q1hY//PCD6NWrV6PbTE5OFgcPHmz0ufPZunVrk6+57xYXFyd27NhRb/2MjIwWHSOAOs/R6dOnhVqtbnadiIgI8eGHH7bpOfLncdbmKztx4sRGy+Tn54tLLrmkyX2FhYWJJUuWNBlz7c/a559/LiIjIxvd3pQpU0RVVVWj26uurhazZs1qdH21Wi1eeumlOp+rVatWNffU1uPxeJT3ZlRUlHA4HE2WdzgcIioqSgAQqampwuPxtHqfQvjvNf7HP/6hlL3kkkuE2+1usNzMmTOVco899liD8TR3rnzqqadaFPP111/flqeky+Jlly5k6dKleP755wHI4xjccsstmDJlCvR6PQ4fPowVK1bg1VdfRW5ubrv2c+uttyq9UW688UYAwIUXXojnnnuuTrnBgwc3u60rrrgCH330EbZs2YK//e1vAOReJVdccUWdcm1tT/LGG29g0aJFyv2f/exnmDZtGiIjI/HDDz9g1apVyMjIwNGjR3HppZfi4MGDiIuLU8o/99xzKCoqwrfffosnnngCAHDLLbfg1ltvbVM8PiUlJbjiiiuQk5MDAOjVqxfuvPNO9OvXDyUlJVi3bh02bdqEm266CRaLpcltVVVVQQiBCy+8EJdffjkGDRqEmJgYVFVVITs7Gx9//DH279+Pc+fO4brrrsOhQ4fq1CTFx8fjo48+QmFhIRYuXAgAuPzyyxv8xRsfH6/8X11dDY/Hg969e+PKK6/E4MGDERcXB6/Xi6ysLGzatAlbtmxBRUUFbr31Vnz11VcYM2ZMm5+z9h5na9jtdkyYMAE//PADALk305133okLL7wQlZWV2LhxI95//324XC48+uij8Hq9zV5CPXToEJYuXQohBBYuXIhLLrkEOp0O+/btw7Jly1BRUYGNGzfij3/8I/7whz80uI0FCxbg/fffBwCEh4dj7ty5uOSSS6BWq7Fv3z6sXLkSDz74IGbOnNmm4/ZRqVS466678Pvf/x6lpaX48MMPMXv27EbLf/DBBygtLQUA/PKXv4RK1bZKdX+9xnfddRc2bNiADz/8ELt378YzzzxT7zlduXIlPvjgAwDA2LFj2zSkwX/+8x+lh5her8dtt92GsWPHIjo6GlVVVcjJycG+ffuwcePGVm+7ywtt7kMNaUvNx6lTp4ROpxMAhFarFZ999lm9Mnl5eWLgwIF1svG21HzU5ttOU78ghWi85sOnvb/UGpKRkaHUeKjVavHee+/VK1NZWSmmTZum7HvmzJkNbqu9tUXn++Uvf6lsb+rUqaKioqJemZdffrlFNSyZmZniyJEjTe7v3//+t1CpVAKAmDt3boNlateAzJkzp9ljKC4uFl999VWTZbZs2SIiIiIEADFp0qRmt9kUfx2nEM2/bxctWqSUufTSS0VZWVm9Mhs2bBDh4eECgNBoNOLQoUP1ytR+36CmVuCHH36oV27Pnj1Co9EoNQ0N1X5s2rRJ2U5sbKw4evRovTIZGRkiLS2tzj7b+nnKz88XYWFhLfp8T5gwQfmcZWdnt2l/Qvj3NS4pKRHJyclKXLVrSk6cOKG8L81mszh9+nSD22juc+87d6jV6iY/Cw6HQ+zZs6fJ4+pumHx0QG1JPn79618r6zz99NONljt8+HCdqvKunHzcf//9yjYfeuihRstZrVbRo0cPAUBIktTgl4M/k4+CggKh1WoFABETEyOKi4sbLVu7Wri9vxXuuOMOAUDo9XpRXV1db3lrk4+WevLJJ5XtZmVl+W27jWnuOIVo+n1bWFioJPJms1nk5eU1uq+lS5cq25o9e3a95ecnH9u3b290W7Nnz26yXO0k+d133210O+fvsz2fp9qXeBr6XAghf5H7ykybNq3N+2qNlrzGQgixbds2JVFJS0sTpaWlorq6WowaNUqJ+V//+lej6zf3uR8wYIAAIIYOHeqPw+pW2OC0i/joo48AAGFhYbj33nsbLTd06FBcddVVwQorpHyDt2k0Gvz2t79ttJzZbMY999wDABBCKM9loHz66adKL6G5c+ciOjq60bK/+93v/LbfcePGAQAcDgeOHDnit+22dL8A8M033wRtf209zk8//VTp3TRnzhwkJiY2Wvaee+6ByWQCAKxfv77RhsuA3EX7sssua3R57UuNx48fr7OsqqoKX3zxBQAgKSkJs2bNanQ7kyZNwtChQxtd3hq1L1muWLGiwTK1H1+wYIFf9tuclr7GEydOxCOPPAIAOHPmDBYtWoTHH38c+/fvBwDMnj0bt99+e5vj8F0OzsnJ8WtD5+6AyUcXUFh
"text/plain": [
"<Figure size 640x480 with 1 Axes>"
]
},
"metadata": {},
"output_type": "display_data"
},
{
"name": "stdout",
"output_type": "stream",
"text": [
"\n",
"Init Params\n",
"Name Value Min Max Stderr Vary Expr Brute_Step\n",
"amp_bec 4.021 0 5.904 None True None None\n",
"amp_th 0.1265 0 5.904 None True None None\n",
"atom_number_bec 1377 -inf inf None False amp_bec / 5 * 2 * 3.14159265359 * sigmax_bec * sigmay_bec None\n",
"atom_number_th 30.54 -inf inf None False amp_th * 2 * 3.14159265359 * 1.20206 / 1.643 * sigma_th * sigma_th None\n",
"condensate_fraction 0.9783 -inf inf None False atom_number_bec / (atom_number_bec + atom_number_th) None\n",
"rot_angle -36 -66 -6 None False None None\n",
"sigma_th 7.247 0 200 None True None None\n",
"sigmax_bec 23.77 0 58 None True None None\n",
"sigmay_bec 11.47 0 28 None True None None\n",
"x0_bec 100.5 90.5 110.5 None True None None\n",
"x0_th 100.5 90.5 110.5 None True None None\n",
"y0_bec 76.51 66.51 86.51 None True None None\n",
"y0_th 76.51 66.51 86.51 None True None None\n",
"\n",
"Name Value Min Max Stderr Vary Expr Brute_Step\n",
"amp_bec 4.021 0 5.904 None True None None\n",
"amp_th 0.1265 0 5.904 None True None None\n",
"atom_number_bec 1377 -inf inf None False amp_bec / 5 * 2 * 3.14159265359 * sigmax_bec * sigmay_bec None\n",
"atom_number_th 30.54 -inf inf None False amp_th * 2 * 3.14159265359 * 1.20206 / 1.643 * sigma_th * sigma_th None\n",
"condensate_fraction 0.9783 -inf inf None False atom_number_bec / (atom_number_bec + atom_number_th) None\n",
"rot_angle -36 -66 -6 None False None None\n",
"sigma_th 7.247 0 200 None True None None\n",
"sigmax_bec 23.77 0 58 None True None None\n",
"sigmay_bec 11.47 0 28 None True None None\n",
"x0_bec 100.5 90.5 110.5 None True None None\n",
"x0_th 100.5 90.5 110.5 None True None None\n",
"y0_bec 76.51 66.51 86.51 None True None None\n",
"y0_th 76.51 66.51 86.51 None True None None\n",
"bval first fit\n",
"amp_bec: 4.093, (init = 4.021), bounds = [0.00 : 5.90] \n",
"amp_th: 0.380, (init = 0.127), bounds = [0.00 : 5.90] \n",
"x0_bec: 101.030, (init = 100.503), bounds = [90.50 : 110.50] \n",
"y0_bec: 75.967, (init = 76.510), bounds = [66.51 : 86.51] \n",
"x0_th: 97.521, (init = 100.503), bounds = [90.50 : 110.50] \n",
"y0_th: 76.807, (init = 76.510), bounds = [66.51 : 86.51] \n",
"sigmax_bec: 25.728, (init = 23.770), bounds = [0.00 : 58.00] \n",
"sigmay_bec: 12.084, (init = 11.466), bounds = [0.00 : 28.00] \n",
"sigma_th: 1.502, (init = 7.247), bounds = [0.00 : 200.00] \n",
"rot_angle: -36.000, (init = -36.000), bounds = [-66.00 : -6.00] \n",
"\n",
"amp_bec: 4.093, (init = 4.021), bounds = [0.00 : 5.90] \n",
"amp_th: 0.380, (init = 0.127), bounds = [0.00 : 5.90] \n",
"x0_bec: 101.030, (init = 100.503), bounds = [90.50 : 110.50] \n",
"y0_bec: 75.967, (init = 76.510), bounds = [66.51 : 86.51] \n",
"x0_th: 97.521, (init = 100.503), bounds = [90.50 : 110.50] \n",
"y0_th: 76.807, (init = 76.510), bounds = [66.51 : 86.51] \n",
"sigmax_bec: 25.728, (init = 23.770), bounds = [0.00 : 58.00] \n",
"sigmay_bec: 12.084, (init = 11.466), bounds = [0.00 : 28.00] \n",
"sigma_th: 1.502, (init = 7.247), bounds = [0.00 : 200.00] \n",
"rot_angle: -36.000, (init = -36.000), bounds = [-66.00 : -6.00] \n",
"\n",
"\n",
"Atom numbers:\n",
" N_bec: 113800\n",
" N_th: 281\n",
" N_ges: 114080\n",
" Cond. frac: 99.75 %\n",
"\n",
"Temperature: 0.80 nK\n"
]
}
],
"source": [
"times = []\n",
"x = np.linspace(0,shape[2],shape[2])\n",
"y = np.linspace(0,shape[1], shape[1])\n",
"\n",
"X,Y = np.meshgrid(x, y)\n",
"X_1d = X.flatten()\n",
"Y_1d = Y.flatten()\n",
"result = []\n",
"for i in range(0, shape[0]):\n",
"# for i in [0]:\n",
" print('\\n ')\n",
" print(f'file {file} \\n image: {i}')\n",
" start = time.time()\n",
" init = fitm.guess(cropOD[i].flatten(), X_1d, Y_1d, rot_angle=36, vary_rot=False)\n",
" init.pretty_print()\n",
"\n",
" res = fitm.fit(cropOD[i].flatten(), x=X_1d, y=Y_1d, params=init)\n",
" fitm.print_bval(res)\n",
" stop = time.time()\n",
" fitm.return_atom_number(res, X_1d, Y_1d)\n",
" fitm.return_temperature(res, omg=1036.5, tof=26e-3, is_print=False)\n",
" fitm.return_temperature(res, tof=26e-3)\n",
"\n",
" # print(f' time = {stop-start:.3f} s')\n",
"\n",
" result.append(res)\n",
" times.append(stop-start)\n",
"\n",
"times = np.array(times)"
]
},
{
"cell_type": "code",
"execution_count": 44,
"metadata": {},
"outputs": [
{
"data": {
"image/png": "iVBORw0KGgoAAAANSUhEUgAABkAAAAkfCAYAAAAqnERsAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjcuMSwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/bCgiHAAAACXBIWXMAAA9hAAAPYQGoP6dpAAEAAElEQVR4nOzdd1gURx/A8e/RQQSxd7D32HtN1CT2EmM32GJNUVNNA9KLb6oliQVLjImaGDXW2GLvvcYC9ooKiEjd949jl7vjDg64Aw5/n+e5h2V3ZnZ2b3d2b2dnRqcoioIQQgghhBBCCCGEEEIIIUQ+4pTbGRBCCCGEEEIIIYQQQgghhLA1qQARQgghhBBCCCGEEEIIIUS+IxUgQgghhBBCCCGEEEIIIYTId6QCRAghhBBCCCGEEEIIIYQQ+Y5UgAghhBBCCCGEEEIIIYQQIt+RChAhhBBCCCGEEEIIIYQQQuQ7UgEihBBCCCGEEEIIIYQQQoh8RypAhBBCCCGEEEIIIYQQQgiR70gFiBBCCCGEEEIIIYQQQggh8h2pABFCCCGEEEIIIYQQQgghRL4jFSBCCCGEEEIIIYQQQgghhMh3pAIkn9LpdNrHVgICArQ0w8PDbZZuRuyxLcL+2rVrp31vW7Zsye3sCAcwYcIEdDodXl5eXLlyJbezI0Secfv2bXx8fNDpdLz44ou5nR0hhLAoIiKCkJAQmjZtip+fH87Oztr94Ny5c3M7e0IIIYQQ4jEkFSBCCIeXlJTE0aNHmT17NmPHjqVRo0a4ublpP7jbtWuX21kUGTh+/DjTpk0D4NVXX6Vs2bJZSmfSpElGlaYBAQEZxpk7d65RHGs+I0eOzFL+ssvex3p8fDwLFiygc+fO+Pv74+HhQalSpWjRogVTpkzhzp07ttmQLLLH9oeHhzNz5kwGDx5M3bp18fPzw9XVlcKFC/PEE08wevRo/v3332zn/eDBg7z99ts0atSIUqVK4e7uTunSpWnQoAHDhw9nwYIF3Lhxw2zcYsWK8dprrwEwZ84c9u/fn+38GDKssA4ODs5U3ODgYIcra8PDwzN9zmf0yex+cwR37txh+fLlvPvuu3Tv3p1atWpp54eXlxdlypThmWee4bPPPuPq1auZTj8vlzdZOSeOHj1KiRIltHiFCxdm37599s1oHnPhwgXq1q1LcHAwe/fu5f79+yQnJ+d2tuwqJ64h2bl+CJEVjlr+5+XrSl5keK3T6XSUL1+euLg4q+Ia3v/179/fzjnNffK8IX22Ovcc6Xe5cECKyJcA7WMr/v7+WpphYWE2Szcj9tgWYX9t27bVvrfNmzfbbT3Lli1TvLy8jI4T00/btm3ttn5hG506dVIApUCBAsqdO3eylMaePXsUJycno+/e398/w3ihoaHpHj/mPiNGjMhSHrPD3sf6qVOnlHr16qWbfvHixZVVq1bZbqMywdbbf/DgQaVJkyZWf+ft2rVTLl68mOl837x5Uxk0aJBV6xg/frzFdO7fv6/4+PhoebElw/I6KCgoU3GDgoIcrqwNCwvL9Dmf0Sez+80RdOnSxertd3d3V4KDg5WkpCSr0s7r5U1mz4n9+/crhQsX1uIUK1ZMOXz4sP0zmiK37tFNdejQQcuHp6en0q1bN2XMmDHK+PHjlfHjxyu7d+/OtbzZWk5cQ2x1/RAisxyx/M/r15W8yPBap36+++47q+Ia3v/169fPzjnNXfK8IX22PPcc5Xe5cEwuCCGEA7t//z4PHz7M7WyIbNixYwdr1qwB4MUXX6RIkSKZTiMhIYGRI0dm+03T6tWr0759+wzDtWjRIlvryQp7HutXrlyhffv2XLt2DdB3PdimTRsqVarE7du32bBhA7Gxsdy6dYuePXuydu1annrqKbvkxRJbb/+ZM2fYu3ev0byqVatSu3ZtihYtyv3799m5c6fWHduWLVto3rw527Zto2LFilat49KlS7Rr146wsDBtXrVq1ahTpw5FihTh4cOHnD9/nsOHD2e4bb6+vowZM4Yvv/ySLVu2sGHDBjp06JDJrRYAPj4+jB8/Pt0we/fu1d7eL126NL169Uo3fJMmTWyWv7yoaNGi1KhRA39/f7y9vXn48CHnzp1j7969JCYmEhcXR3BwMBcuXGDevHnppuUI5U1m7Nq1i06dOhEZGQlAqVKl2LhxIzVq1MjlnOWs69evs2HDBgDc3d05cuQIVapUyeVc2Y+9ryG2vH4IkR2OUP7nt+tKbvr0008ZOXIkXl5euZ2VPEOeN1hmz3MvL/8uFw4qt2tghH1gUCPq6PLTtjxOcqoFiPqWQIkSJZSuXbsqISEhyurVq5VXX31V3shwEF27dlUARafTKf/991+W0vjoo4+073vgwIHadGZbgAQGBmZp/TnBnsd669atjfaZ6ZvLt2/fVtq3b6+FKVy4sHLv3r3sb1Qm2Hr7Fy1apABK5cqVlc8//1y5cuVKmjBJSUnK7Nmzjd76atasmZKcnJxh+vfv31cqVqyoxXvyySeVI0eOmA0bFxenrFmzRlm8eHG6aZ4/f17R6XQKoHTo0MG6DbXC49YCxBr5dbsy46uvvlJ+/PFH5ezZsxbD3LhxQxkwYIDRvdqSJUvSTdcRyhtrz4ktW7Yo3t7eWthy5cqlu7/sJS+0AFm1apWWB1uWT3mVPa8h9rh+CJEZjlb+O8J1JS8y1wIEUD777LMM4z5OLUDkeYNltj73HOV3uXBM8kQ5n8pPlQb5aVseJzlVAXL9+nWzXQrIwyvH8N9//2kPdLP6PZ06dUpxd3dXAGXQoEFGN075qQLEXse64UMrNzc35ejRo2bDPXjwwOiBzOTJkzO9ruyw9fZv2bJFCQ0NVRITEzMM++effxpdi9auXZthnJEjRxr9MLRmPdYw/BFx7Ngxm6QpFSBp5dftsofk5GTlqaeesurht6OUN9acE+vXrzd6sF2xYkUlPDw8R/OpygsVIL/88ouWhyFDhuRKHnKSPa8h9rp+CGFreaH8d5TrSl5keK1r1qyZ0UPqyMjIdOM+ThUg+el5g+Hv3tDQ0GylZY9zz1F+lwvHJIOgCyEcWsmSJSlfvnxuZ0NkUWhoKIqiANCvX79Mx1cUhZEjRxIXF4efnx9ff/21rbOYZ9jrWFcHnwcIDAykTp06ZsMVKFCADz/8UPv/p59+IjEx0eb5scTW29+2bVuGDh2Ks7NzhmF79epl1MXRqlWr0g1/+PBhZs2aBUC5cuWYOXOmVeuxRt++fbXp0NBQm6QpRHbodDqGDRum/X/o0CGLYR2lvMnI33//Tbdu3bQuMapVq8bWrVvx9/fP5ZzlnoSEBG3aySn//8S01zXEntcPIWwtL5T/+eW6ktsGDx5MtWrVALh79y7/+9//cjlHeYc8bzBPzj3haPL/3anQnDlzhgkTJlCjRg28vb3x8fGhbt26TJ48mTt37mQYPyAgAJ1Oh06nIzw8PN2wSUlJzJ49mw4dOlCiRAk8PDwICAigR48eLFu2THvg2a5dOy3NLVu25Ni22FpUVBQ//PAD3bp1IyAgAG9vb9zd3SldujTt27cnJCSEEydOZJjO7t27eemll6hVqxZ+fn54eHhQtmxZnn32WaZOnUpMTIxV+VH3qU6n0+bZap8lJyczb948OnbsSMmSJY2+27/++svqdB4nffv21b6PV1991WK4AwcO4Obmhk6no169eo/FjcHChQu16Z49e2Y6/owZM9ixYwcAX331FcWLF7dV1h4LDx48YOPGjdr/hj9izXnuuefw9vYG9D+Otm7datf85SUtW7bUpjO6Bv7444/a9Pjx4ylYsKDN8tGzZ0+tbP/111+162letWXLFqNrUmY+wnEUK1ZMm46OjjYbxl7ljT3vecz5888/6d27N3FxcQDUqlWLf//9lzJlymQr3cwKDw/XtvvixYva/AoVKpg9n9K7z87q/afh+W34fc6bNy/N+ocOHWqrTXdI1l5D7HH9MHeOHD58mLFjx1KtWjW8vb3x9vamadOmTJ8+3ez95/79+xk6dCg1atSgQIECFClShCeffNLoPi4jly9fJiQkhDZt2lCiRAnc3d1xc3OjSJE
"text/plain": [
"<Figure size 2000x2800 with 42 Axes>"
]
},
"metadata": {},
"output_type": "display_data"
}
],
"source": [
"# fig, axs = plt.subplots(shape[0] , 5, figsize=(14, 4 * shape[0] ),dpi = 100)\n",
"\n",
"fig, axs = plt.subplots(shape[0] , 6, figsize=(20,4 * shape[0] ),dpi = 100)\n",
"fig.suptitle(f'{file}: ν = ({data_x.trap_f[0]}, {data_x.trap_f[1]}, {data_x.trap_f[2]}) Hz, T = {data_x.temp[0]*1e9:.0f}nK, tof = {data_x.tof * 1e3:.0f}ms, a = {data_x.a_fac :.0f}a$_0$, N = {data_x.N:.0e} ', y=-0.005*shape[0]+0.965, fontsize=24)\n",
"# fig.suptitle('test', va='bottom')\n",
"\n",
"for i in range(0,shape[0]):\n",
"\n",
" axs[i,0].set_title(f'image {i}, cond. frac = {fitm.cond_frac(result[i], X, Y) *1e2:.2f} %', loc='left', fontsize=18)\n",
"\n",
" lmfit.fit_report(result[i])\n",
" bval = result[i].best_values\n",
" fit = dp.density_profile_BEC_2d(X,Y, **bval)\n",
" vmax = np.max(fit)\n",
"\n",
" if bval['amp_bec'] > bval['amp_th']:\n",
" cen_str = 'bec'\n",
" else:\n",
" cen_str = 'th'\n",
" cen_x = round(bval[f'x0_{cen_str}'])\n",
" cen_y = round(bval[f'y0_{cen_str}'])\n",
"\n",
" ax = axs[i,0]\n",
" ax.pcolormesh(X, Y, cropOD[i], vmin=0, vmax=vmax, cmap='jet', shading='auto')\n",
" #plt.colorbar(art, ax=ax, label='z')\n",
"\n",
"\n",
" # Plot gaussian 2d Fit + legend including Width parameters\n",
" ax = axs[i,1]\n",
"\n",
" ax.pcolormesh(X, Y, fit, vmin=0, vmax=vmax, cmap='jet', shading='auto')\n",
" #plt.colorbar(art, ax=ax, label='z')\n",
"\n",
" ax = axs[i,2]\n",
"\n",
" ax.pcolormesh(X, Y, fit-cropOD[i], vmin=0, vmax=0.2, cmap='jet', shading='auto')\n",
"\n",
"\n",
" ax = axs[i,3]\n",
"\n",
" ax.plot(x, cropOD[i, cen_y, :])\n",
" ax.plot(x, fit[cen_y, :], label='bimodal')\n",
" ax.plot(x, dp.thermal(x, bval['x0_th'], bval['amp_th'], bval['sigma_th']), label='thermal')\n",
" ax.legend(fontsize=8)\n",
"\n",
" ax = axs[i,4]\n",
"\n",
" ax.plot(y, cropOD[i, :, cen_x])\n",
" ax.plot(y, fit[:, cen_x], label='bimodal')\n",
" ax.plot(y, dp.thermal(y, bval['y0_th'], bval['amp_th'], bval['sigma_th']), label='thermal')\n",
" ax.legend(fontsize=8)\n",
"\n",
" ax = axs[i,5]\n",
" Na = fitm.return_atom_number(result[i], X, Y, is_print=False)\n",
" ax.axis('off')\n",
" col_labels = ['Sim', 'Fit']\n",
" row_labels = ['N', f'N$_{{bec}}$', f'N$_{{th}}$','f', 'T (nk)']\n",
" table_vals = [[f'{data_x.N:.2e}',f\"{Na['N']:.2e}\"],\n",
" [f'{data_x.cond_frac.data[i] * data_x.N:.2e}',f\"{Na['N_bec']:.2e}\"],\n",
" [f'{(1-data_x.cond_frac.data[i]) * data_x.N:.2e}',f\"{Na['N_th']:.2e}\"],\n",
" [f'{data_x.cond_frac.data[i] *1e2:.1f} %', f\"{Na['cond_f'] *1e2:.2f} %\"],\n",
" [f'{data_x.temp[i] *1e9:.0f} nK', f'{fitm.return_temperature(result[i], tof=data_x.tof, is_print=False)*1e9:.1f} nK']\n",
" ]\n",
"\n",
"\n",
" my_table = ax.table(cellText=table_vals,\n",
" colWidths=[0.4]*2,\n",
" rowLabels=row_labels,\n",
" colLabels=col_labels,\n",
" loc='upper right')\n",
" my_table.scale(1,2)\n",
" ax.text(0.1,0.1,f'fitting time = {times[i]:.2f} s', fontsize=14)\n",
"\n",
"t_fonts=20\n",
"axs[0,0].set_title(f'Data \\n ', fontsize=t_fonts)\n",
"axs[0,1].set_title('Fit \\n ', fontsize=t_fonts)\n",
"axs[0,2].set_title('Data - Fit \\n ', fontsize=t_fonts)\n",
"axs[0,3].set_title('cut along x \\n ', fontsize=t_fonts)\n",
"axs[0,4].set_title('cut along y \\n ', fontsize=t_fonts)\n",
"\n",
"\n",
"# plt.tight_layout()\n",
"# plt.savefig(f'fit_output/final/{file}_{saveprefix}.png')\n",
"plt.show()"
]
},
{
"cell_type": "code",
"execution_count": 32,
"metadata": {},
"outputs": [
{
"data": {
"text/plain": [
"array([3.e-08, 3.e-08, 3.e-08, 3.e-08, 3.e-08, 3.e-08, 3.e-08])"
]
},
"execution_count": 32,
"metadata": {},
"output_type": "execute_result"
}
],
"source": [
"data_x.temp"
]
},
{
"cell_type": "code",
"execution_count": null,
"metadata": {},
"outputs": [],
"source": []
}
],
"metadata": {
"kernelspec": {
"display_name": ".venv",
"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.13"
},
"orig_nbformat": 4
},
"nbformat": 4,
"nbformat_minor": 2
}
Reference in New Issue Copy Permalink