2021-10-11 18:31:02 +01:00
{
"cells": [
{
"cell_type": "code",
"execution_count": 1,
"id": "7c1065dd",
"metadata": {},
"outputs": [],
"source": [
"import numpy as np\n",
"import matplotlib.pyplot as plt\n",
"import pyerrors as pe"
]
},
{
"cell_type": "code",
"execution_count": 2,
"id": "20f67709",
"metadata": {},
"outputs": [],
"source": [
"plt.style.use('./base_style.mplstyle')\n",
"plt.rc('text', usetex=True)"
]
},
{
"cell_type": "markdown",
"id": "e5764fd0",
"metadata": {},
"source": [
2022-01-06 17:38:12 +01:00
"We can load data from a preprocessed file which contains a list of `pyerror` `Obs`:"
2021-10-11 18:31:02 +01:00
]
},
{
"cell_type": "code",
"execution_count": 3,
2022-01-06 11:16:00 +01:00
"id": "fbfa65f5",
2021-10-11 18:31:02 +01:00
"metadata": {},
2022-01-06 11:16:00 +01:00
"outputs": [
{
"name": "stdout",
"output_type": "stream",
"text": [
"Data has been written using pyerrors 2.0.0.\n",
"Format version 0.1\n",
"Written by fjosw on 2022-01-06 11:11:19 +0100 on host XPS139305, Linux-5.11.0-44-generic-x86_64-with-glibc2.29\n",
"\n",
"Description: Test data for the correlator example\n"
]
}
],
2021-10-11 18:31:02 +01:00
"source": [
2022-01-06 11:16:00 +01:00
"correlator_data = pe.input.json.load_json(\"./data/correlator_test\")"
2021-10-11 18:31:02 +01:00
]
},
{
"cell_type": "markdown",
"id": "ae93c7c2",
"metadata": {},
"source": [
"With this list a `Corr` object can be initialised"
]
},
{
"cell_type": "code",
"execution_count": 4,
"id": "33a8fdec",
"metadata": {},
"outputs": [],
"source": [
2022-01-06 17:38:12 +01:00
"my_correlator = pe.Corr(correlator_data)"
2021-10-11 18:31:02 +01:00
]
},
{
"cell_type": "code",
"execution_count": 5,
"id": "5f954607",
"metadata": {},
"outputs": [
{
"name": "stdout",
"output_type": "stream",
"text": [
"x0/a\tCorr(x0/a)\n",
"------------------\n",
2022-01-06 11:16:00 +01:00
"8\t 548(13)\n",
"9\t 433(11)\n",
"10\t 343.1(8.6)\n",
"11\t 273.2(6.6)\n",
"12\t 217.5(5.6)\n",
"13\t 172.9(4.9)\n",
"14\t 137.6(4.6)\n",
2021-10-11 18:31:02 +01:00
"\n"
]
}
],
"source": [
"my_correlator.print([8, 14])"
]
},
{
"cell_type": "markdown",
"id": "b00d670b",
"metadata": {},
"source": [
"The `show` method can display the correlator"
]
},
{
"cell_type": "code",
"execution_count": 6,
"id": "b71529d0",
"metadata": {},
"outputs": [
{
"data": {
2022-01-06 17:51:50 +01:00
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAlYAAAGLCAYAAAAF7dxzAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjUuMCwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8/fFQqAAAACXBIWXMAAA9hAAAPYQGoP6dpAAAdmklEQVR4nO3dUYicR4If8H/ZPizu0MysNizWxbuJRgQSnyHc2F44yL7k5M2G4DwEyb4L+7o2JoGQFymCC8bsg2I9XEJIMJJfTXK3EuESX8Dg2UC4wIVbry4BYUISj5M972mz2ZV6xne31m7sykN/LbdaPTPdo+rpGfXvB43m+6q6p3pKmu+vqvqqS601AADcv4fm3QAAgAeFYAUA0IhgBQDQiGAFANCIYAUA0IhgBQDQiGAFANDII/NuwH4qpZQkv5jko3m3BQA4VI4m+eO6ywagCxWs0g9VH867EQDAofR4ku/vVGHRgtVHSfJHf/RHWVpamndbAIBDYGtrK1/84heTCWa8Fi1YJUmWlpYEKwCgOYvXAQAaEawAABoRrAAAGhGsAAAaEawAABoRrAAAGhGsAAAaEawAABoRrAAAGhGsAAAaEawAABoRrAAAGhGsAAAaEawAABoRrAAAGhGsAAAaEawAABoRrAAAGhGsAAAaEawAABp5ZN4NmIdPPq35/fd/nB9+9HG+cPRIvnziWB5+qMy7WQDAIbeQweqr//Q/5v/efvjO8fHlI3nluSfytSePz7FVAMBht5BTgf9n6/Zdxz/Y/Dgvv3ktb1+/MacWAQAPgoUMVqNq9/iN37meG72fzLs5AMAhJVgN+dGf/DTf/N335t0MAOCQEqxG/ORnn8y7CQDAISVYjXj+6S/OuwkAwCG1kHcFjlOSPLZ8JF/9pcfm3RQA4JBayBGr0R2rBsevPPeE/awAgD1byGD1haVH7zp+bPlIXv/6mn2sAID7Umqt827DvimlLCXZvHmrl//24/9n53UAYFdbW1tZXl5OkuVa69ZOdRdyjdXDD5X8ysnPz7sZAMADZiGnAgEAZmFPI1allNeSvN8d3qy1Xh0qO5uk1x2u1Fovjjx3puUAAPMyVbAqpawk+XaSX6219kopa0m+m+7Gui70pNZ6uTs+VUq5VGt9aT/KAQDmaarF66WUS0neHx4lKqWcqrWud1/fSnKi1tobKq+11rIf5RO0fynJ5ubmZpaWliZ+3wDA4ppm8fq0a6xeTHK1lLJaSjmVJEOhajX9qbne6JO6kaWZlo9rbCnl0VLK0uCR5OjkbxUAYDoTB6su2CTJWpKVJBullEtDoWZ17BP766FW9qF8nPNJNoceH25TDwDgvk0zYjUINr1a67Va60aSc0mu7PK8m0mOzan8QpLlocfjO7wOAMB92ctdge8OvugWsK9sNxXX2SkUzbS81no7ye3BcSk2AQUAZmeaEauNbc730h/N2q58pSubdTkAwFxNHKy6qb+N3LvWaSXJu115b2gt1vBz12ddPun7AACYlWnvCjyX5IXBQSnldJL1Wuu17tSFJKdGyi8PPX/W5QAAczP1hzCXUl7MZ3fhfb7Wem6k/Gw+m5p7Zr/Ld2m7fawAgKlMs4/V1MHqMBOsAIBpzXKDUAAAtiFYAQA0IlgBADQiWAEANCJYAQA0IlgBADQiWAEANCJYAQA0IlgBADQiWAEANCJYAQA0IlgBADQiWAEANCJYAQA0IlgBADQiWAEANCJYAQA0IlgBADQiWAEANCJYAQA0IlgBADQiWAEANCJYAQA0IlgBADQiWAEANCJYAQA0IlgBADQiWAEANCJYAQA0IlgBADQiWAEANPLINJVLKaeSvJTknSQbSZ5N8p1a69WhOmeT9LrDlVrrxZHXmGk5AMC8TDtitZLkVJJL3eP9MaEqtdbLtdbLSa6VUi7tVzkAwDyVWuvklUs5nWS91trbpvxWkhPD5aWUWmst+1E+QfuXkmxubm5maWlpovcMACy2ra2tLC8vJ8lyrXVrp7rN1liVUlbTn5rrjSk7Nevybdr0aCllafBIcnTKtwUAMLGp1lh1ni+l3ExyLMnJWuu57vzqNvV76U8hzrp8nPNJXtmmDACgqWlHrK6lPxV4tVvj9H4p5couzxmEsHmUX0iyPPR4fIfXAQC4L1ONWNVaN0ZOfSvJpVLKyg5P2ykUzbS81no7ye3BcSkTLcUCANiTqUasusXrdwytd1pNf/uFcVa6slmXAwDM1cTBqhuVutItIh8+lyQb3WhWb7h8oNa6PuvySd8HAMCsTBysutGpiyPTgS8muTo0cnUh/X2uktwZ4bo8VH/W5QAAczPtPlYr6Yepgc8P3RU4qHM2n03NPbPf5bu03z5WAMBUptnHaqpgddgJVgDAtOayQSgAwKITrAAAGhGsAAAaEawAABoRrAAAGhGsAAAaEawAABoRrAAAGhGsAAAaEawAABoRrAAAGhGsAAAaEawAABoRrAAAGhGsAAAaEawAABoRrAAAGhGsAAAaEawAABoRrAAAGhGsAAAaEawAABoRrAAAGhGsAAAaEawAABoRrAAAGhGsAAAaEawAABoRrAAAGhGsAAAaEawAABp55H6eXEp5p9b67Mi5s0l63eFKrfXifpYDAMzLnkesSimnk5waOXc2SWqtl2utl5NcK6Vc2q9yAIB5KrXW6Z9UykqS55NcqrWWofO3kpyotfaGztVBnVmXT9DupSSbm5ubWVpamvZtAwALaGtrK8vLy0myXGvd2qnuXkesnk/yreETpZTV9KfmeqOVSymnZl0+rpGllEdLKUuDR5Kjk7w5AIC9mDpYdSFmfUzR6jZP6SVZ2Yfycc4n2Rx6fLhNPQCA+7aXEauVWuvGFPVvJjk2p/ILSZaHHo/v8DoAAPdlqrsCSykvdovGp7FTKJppea31dpLbg+NSJlqKBQCwJxOPWJVS1pK8u0OV7UaxVrqyWZcDAMzVNCNWx5KsDS0UP5nc2QJho9Z6tZTSK6Wsjk4V1lrXu7ozLQcAmKeJR6xqreu11ouDR5JL3fmLtdarXbULGdrbqtvranjqcNblAABzs9d9rE4neSHJ6SQXk7wzNKp0Np9NzT1Taz038tyZlu/SbvtYAQBTmWYfqz0Fq8NKsAIAprUfG4QCADBCsAIAaESwAgBoRLACAGhEsAIAaESwAgBoRLACAGhEsAIAaESwAgBoRLACAGhEsAIAaESwAgBoRLACAGhEsAIAaESwAgBoRLACAGhEsAIAaESwAgBoRLACAGhEsAIAaESwAgBoRLACAGhEsAIAaESwAgBoRLACAGhEsAIAaESwAgBoRLACAGhEsAIAaESwAgBoRLACAGjkkWkql1JWkjzfHZ5MsprkG7XW3lCds0kGxyu11osjrzHTcgCAeZl2xOq1JOu11su11nNJbia5MijsQk+68stJrpVSLu1XOQDAPJVa6+SVS3knyTuDUaIu6JyvtX6uO76V5MTICFattZb9KJ+g/UtJNjc3N7O0tDTx+wYAFtfW1laWl5eTZLnWurVT3ammAmutz46ceibJepKUUlbTn5rrjT6vlHIqycYsy2ut69O8FwCA1qYKVsNKKaeTrCQ5051a3aZqr6s36/JxbXw0yaNDp45u8xoAAPdt6rsCSykrpZQX0w86V8aNII24meTYnMrPJ9kceny4w+sAANyXqYNVrbXXLR4frLO61d0tuJ2dQtGsyy8kWR56PL7LawEA7NnEwaobqXptJEStpz8NN1gDNc5KVzbr8nvUWm/XWrcGjyQfbfMaAAD3bZoRq9UkZ3P3CNFK92ev1rqRpNctYr9LrXV91uVTvA8AgJmYOFjVWq8ludgFnIEXklwbCjYX0h+9SnJngfvlofqzLgcAmJtp97FaSfLi0KmTSc6N2Xl9EL6e6TYSzX6V79J++1gBAFOZZh+rqYLVYSdYAQDTmiZY+RBmAIBGBCsAgEYEKwCARgQrAIBGBCsAgEYEKwCARgQrAIBGBCsAgEYEKwCARgQrAIBGBCsAgEYEKwCARgQrAIBGBCsAgEYEKwCARh6ZdwMOkk8+rfmDD27mhx99nC8cPZIvnziWhx8q824
2021-10-11 18:31:02 +01:00
"text/plain": [
"<Figure size 640x395.55 with 1 Axes>"
]
},
"metadata": {
"needs_background": "light"
},
"output_type": "display_data"
}
],
"source": [
"my_correlator.show()"
]
},
{
"cell_type": "markdown",
"id": "c659557e",
"metadata": {},
"source": [
"## Manipulating correlators"
]
},
{
"cell_type": "markdown",
"id": "416cf39a",
"metadata": {},
"source": [
"`Corr` objects can be shifted"
]
},
{
"cell_type": "code",
"execution_count": 7,
"id": "e8d65dd5",
"metadata": {},
"outputs": [],
"source": [
"shifted_correlator = my_correlator.roll(20)\n",
"shifted_correlator.tag = r'Correlator shifted by $x_0/a=20$'"
]
},
{
"cell_type": "markdown",
"id": "634dd613",
"metadata": {},
"source": [
2022-01-06 17:38:12 +01:00
"or symmetrised"
2021-10-11 18:31:02 +01:00
]
},
{
"cell_type": "code",
"execution_count": 8,
"id": "127a661d",
"metadata": {},
"outputs": [],
"source": [
"symmetrised_correlator = my_correlator.symmetric()\n",
"symmetrised_correlator.tag = 'Symmetrised correlator'"
]
},
{
"cell_type": "markdown",
"id": "3d733872",
"metadata": {},
"source": [
2022-01-06 17:38:12 +01:00
"We can compare different `Corr` objects by passing `comp` to the `show` method"
2021-10-11 18:31:02 +01:00
]
},
{
"cell_type": "code",
"execution_count": 9,
"id": "8e264aed",
"metadata": {},
"outputs": [
{
"data": {
2022-01-06 17:51:50 +01:00
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAk0AAAGNCAYAAAAM+kVxAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjUuMCwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8/fFQqAAAACXBIWXMAAA9hAAAPYQGoP6dpAABhB0lEQVR4nO3df3xU1Z0//tcJyo8KZAgUiYDosCIGrDiAP7bSbjVBaP21bfihZeuntSZ1ddut2mTj6pe12uaTtNjdbi2b0B9b1x9I4laU7wpkpFuxVSGJFCFWNOMPxCCFMAEsP2rmfP44cyd3JvfO3DuZmXvvzOv5eMxjMnfuzJyZm+S+55z3eR8hpQQRERERJVfkdAOIiIiIvIBBExEREZEFDJqIiIiILGDQRERERGQBgyYiIiIiCxg0EREREVnAoImIiIjIAlcFTUKIFiFEIHppcLo9RERERJrTnG5AAj+A5wG0A1jicFuIiIiIYoSbKoILISqllK1pPlYAOAvA0cy2ioiIiPLcGAAfyBRBUcZ7moQQPgBLASyRUlYY3F8DIBy96ZNSNurunq9iH5QAgJSy2cZLnwXg/TSaTERERDQFwL5kO2S0p0kIEQAwD4APwDIp5dyE+2sAQAuUhBDlUMFVtcFzdQOYK6UMW3ztsQD69u7di7Fjxw7lbRAREVGBOHLkCKZOnQoAxVLKI8n2zWhPk5SyE0CnEKLSZJc6AOfq9g8KIdoAVEcfM19KWRu9OwyV49Rppw1jx45l0EREREQZl7NEcCGEH2o4LmxwXzmAEAaG7RDd11bARERERJQtuZw95zfZHoYKkIJCiEqtxwnAoHwoPSHECAAjdJvGZKSVRERERAbcUHKgFwOJ39rMOSsz6OoArMxWo4iIiIj03FDcsiTNx9UDKNZdpmSsRUREREQJctnTFDLZ7ktynykp5UkAJ4UQtwO4He4IAImIiChP5SzQkFKGAISjCeGJ9wWH8LwPSynLAFwylPYRERERJZOtoMlsyK0eQLl2I5r0baeAJREREZEjMjo8F+1FqgSwDIC26O52LcFbStkohKjR1XGab1TY0uZrxg/PRfqBt7cCxz4ERp8JTPtroGjYUF6CcqA/IrHt7V4cOHoCE8eMxCXnlmBYkXC6WURERDGuWntuKLSK4H31F2DsCV0VdN/ZwMLvAWXXOdY2Sm7jrh7c/2wXevpOxLaVFo/EymvLsGh2qYMtIyKifHfkyBEUFxcDFiqC51/y9IQLgFuCQN0+dT1xFrDuK0DXM063jAxs3NWD2x7tjAuYAGB/3wnc9mgnNu7qcahlRERE8Tzf05QwPHd+3+HDGOvzDewQiQBrbwIOdAHffJVDdS7SH5G4omHLoIBJIwBMKh6JF2uv5FAdGTpw5AQOHD2Zcr+JY0Zg4tiROWgREXmNnZ4mNxS3HBIp5cMAHtaG51CU0HlWVAQsuBP4eQXw6qPA3JudaCYZaGl/zzRgAgAJoKfvBFra38PyS6blrmHkGY+98h7+7fk3U+73ravOw7crZuSgRUSUzzwfNFky8QJ1/dJPGDS5yJqtb1vej0ETGfnypWejouzM2O23DhzDPz65A/+6bA7+auLo2PaJY0YYPZyIosLhMHz6URoy5PmcJiHE7UKILgDbTHc68Lq6vvyO3DSKLLl1wbkZ3Y8Kz8SxIzF7cnHsogVKfzVxdNx2J4bmamtrUVtbi8bGRjQ3N6O1tTW2Pdc6OztRUVGB6dOn5/y106W1ee7cuWk9PhgMorq6GtXV1bHPfqiqq6sxbtw4BIPJSwsGg0HMnTsXFRVJl1B1ldraWoTDYaeb4XqeD5oGFbeMROJ3iESArQ8BvmnAxSty30AytWTe2SgtHgmzbCUBNYtuybyzc9ks8qj+iMTO98MAgJ3vh9EfcSZfs7OzM3bCbGhoQE1NDaqqqhAIBLBkyZKMncDtCAQCaQdrzc3OlNLT2pzqRF5bW4slS5YM2q59/kuWLIn1oAz1vTQ1NcHvN1t7fkB5eTnq6urQ29s7pNfLNTf0NDU2NqKxsTEW8Brd39zcjObmZjQ2Nua8fZ4PmgZpvQXYuw04eVRdr70J2LMRWPggk8BdZliRwMprywBgUOCk3V55bRmTwCmljbt6cEXDFtzz610AgHt+vQtXNGxxZPblkiVL0NDQgPLy8rjtfr/f8CSQKyUl6S3z2dbWluGWWGelzRUVFVi2bFncts7OTvj9fvh8PpSXl8eORS7fixsCEKtaW1td0StWW1uLmpoa1NTUoKmpCQDi2qUFSVVVVbEvIrn+m8q/nKaDr6ukb41vGrD0EdZpcqlFs0uxekVgUJ2mSazTRBZpZSsS+5W0shWrVwRy9nuk9eYkBkya8vJySz0VbtHc3IxQyPbSoDll9lknBi1eeC9OefLJJ9HS0uJoG8LhMDo7O+Nyq6qrqzF37lyEQiH4/X7U19fj7bcHcmHLy8tRUVERC7BywfNB06CK4N/4HXB4FyuCe8ii2aWoKJvEiuBkW39E4v5nuwYFTICafSkA3P9sFyrKJuXk96m1tdX0JK5JHCZrbGyMBVKhUAg1NTUAVF5MbW1trIdK6yWpqKgw3N7Q0BD3fNqJprKyEmbC4TCam5vh9/vR1taG6upqBAKB2Ou3tbUhFArFvuFrbUun3Vr7jGhtCIfDCIVC8Pl8qKqqit3f2dmJUCiEUCiEQ4cOxZ6rs7MTtbW1CIVC6O7ujm1ramqKtVvrcbLyXow+s8bGRvh8vrR76rTh2N7eXoTD4djrtra2or6+HqFQCC0tLSgvL0coFEJFRQX8fn/KoUDtedva2lBbW4tgMIju7m5UV1fbCszD4bDhewuHw6ivr8f8+fMBqKHJlpaWrPagtbe3IxQKxX4Htfeh/V6YJasHg8GUf3cZI6XMiwuAsQBkX1+fJKLC8Pu3DspptRtSXn7/1sGctAeAbGhosLx/ZWWlbGtri93u7u6W5eXlsdstLS0yEAjItrY22dHRIWtqapJur6yslC0tLbHHl5eXy46ODimllB0dHdLv98e9fk1Njezu7o7d9vv98vDhw7HbbW1tMhAIZKzdRlpaWmRTU1Pcc2m3Ozo6pM/ni3stv98fe09m78toW7L3YvaZ1dTUxLXt8OHDEkBce8y0tbVJAHGfZ1NTk6yqqorbJ7GdVn5/mpqaYs9bVVUlKysrY23Xvxcrmpqa4j5PKdX7DAQCsdfo6OiQKlzIrZaWlthnqH2eiXw+n+33nKivr09Cfc8aK1PEGp7vaSKiwnXgqHmdr3T2y6XOzk4Eg8G4YRG/34/e3t7YN2efz4fOzs7Yt2jtG7jR9lAohNbW1rjnW7JkCZqamkyHL0KhEILBYKxXx+/3IxgMJu2dGkq7zbS0tGDp0qXw+Xzw+/2YN29e7L5wOBzXi6D1CKV6TiuSfWYNDQ1obGzUvpQDUJ+7ndcNBAJxPSNVVVUQQsR64srLy9Hb24vOzs5BvSvJlJSUxJ43FArF8nrSydlqa2uL69UDVG/osmXLYq/R29ub8n1bzS2aO3fuoNczU19fj6ampqS9WyUlJTlNuGfQRESeNXGMtVICVvcbKr/fHxsmMqMNAbW3txueILWhMi1QMDuJJm4PBoPw+Xxx0+G7u7uT5vFowYI2/NHb25vyBDTUdieqrKxEU1MTxo0bh0AggGXLlsUNnSU+j8/ny9hJMtlnpt2XaX6/P5aoDqhA6sknn0QgELA8zKQPatvb2w3zkfTDglqAlkj7XUzU3Nwc93usD4DNZDqvSAvcUgVYuZ6hyKCJiDzrknNLUFo8Evv7ThjmNWlL8Vxybnr5KHZVVlamLCmg9exYrYljduJO3B4OhwedHFOd6Do7O1FfX4+KigosXbo0ZaCj5ZVYYSfgaGtri/VgaSdffeCUDdp7MfvMWltb085jSibxJL9s2TJcddVVaGhoQCgUspWbowXEiZ91KBRCW1tb3Aw0o+dtamoa1EPU2dkJID5Y1fKmcqW1tRXTp0+PC5jMfje1Y5grni85YKm4JRHlJbeVrdASlM2KH+qTbrXE30ShUCiWfGuHNkRn9JpmbbnqqqtQV1eHqqoq+Hy+2L5
2021-10-11 18:31:02 +01:00
"text/plain": [
"<Figure size 640x395.55 with 1 Axes>"
]
},
"metadata": {
"needs_background": "light"
},
"output_type": "display_data"
}
],
"source": [
"shifted_correlator.show(comp=symmetrised_correlator, logscale=True)"
]
},
{
"cell_type": "markdown",
"id": "232e88af",
"metadata": {},
"source": [
"## Effective mass"
]
},
{
"cell_type": "markdown",
"id": "83dc751c",
"metadata": {},
"source": [
"The effective mass of the correlator can be obtained by calling the `m_eff` method"
]
},
{
"cell_type": "code",
"execution_count": 10,
"id": "c686f7e0",
"metadata": {},
"outputs": [],
"source": [
"m_eff = symmetrised_correlator.m_eff()\n",
"m_eff.tag = 'Effective mass'"
]
},
{
"cell_type": "markdown",
"id": "4a9d13b2",
"metadata": {},
"source": [
"We can also use the priodicity of the lattice in order to obtain the cosh effective mass"
]
},
{
"cell_type": "code",
"execution_count": 11,
"id": "5acde8cf",
"metadata": {},
"outputs": [],
"source": [
"periodic_m_eff = symmetrised_correlator.m_eff('periodic')\n",
"periodic_m_eff.tag = 'Cosh effective mass'"
]
},
{
"cell_type": "markdown",
"id": "c658b000",
"metadata": {},
"source": [
"We can compare the two and see how the standard effective mass deviates form the plateau at the center of the lattice"
]
},
{
"cell_type": "code",
"execution_count": 12,
"id": "1d6ea22a",
"metadata": {},
"outputs": [
{
"data": {
2022-01-06 11:16:00 +01:00
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAmkAAAGLCAYAAACcFQXGAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjUuMCwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8/fFQqAAAACXBIWXMAAA9hAAAPYQGoP6dpAABR2klEQVR4nO3dfXhc5WHn/e8tORgClseyLRsT43iEKdgkAWGFQHhJgoRJSdikle2QZEP7ZC2VazchbVIpgnZ5aJcociFPWLpLJTYl5EkBIzXJkqTBkeg2QEPB9piQ2FBsDeHNxsbII/MSDJbu/eOcMzozmrczo5GORr/PdemSzjn3nLmPjmX9dL8dY61FRERERMKlarorICIiIiITKaSJiIiIhJBCmoiIiEgIKaSJiIiIhJBCmoiIiEgIKaSJiIiIhJBCmoiIiEgIzZnuCsxkxhgDLANem+66iIiIyIwyD9hncyxYq5BWmmXAi9NdCREREZmR3gO8lO2gQlppXgN44YUXqKmpme66iIiIyAxw5MgRli9fDnl64hTSJkFNTY1CmoiIiEwqTRwQERERCSGFNBEREZEQClV3pzGmHUi4mxFr7eY85SPABnezHogCm6y1Cfd4H7AFiPvOC4C1Nm6MaQLagAG3TDOwzVrbX/rViIiIiBQvNCHNDWhYa3vd7SZjTI+1ti3Hy7qBbmtt3H1ND9CHE7YAGoCWDK/rB9YDEaDJLRN3z6WAJiIiItMuTN2dnUCvt2GtHQRa87wmSmoIGwLW+rZ7rLXG/wG0WWvX+8qsdI/VewFRREREZLqFoiXNGBPF6d5MZDjW5Aa2Cay1zWm7GgF/2ZRWMbd7c3sJ9ZwLzPXtmlfsuURERERyCUVIw2kRyySB0yWZlzGmxS2bbCXzukHd41EgmiHwbTDGDAO1QL21tiPH23QCNxRSHxEREZFShCWkZeOFp6x8kwciQF+m1jhXR4bxbTEYD3PGmFZjTF9ad6hfF/At3/Y89MQBERERKYOwh7ScAQ3ADWXeZINWY8xhnHFmCa+MMaYhy2vjabvuA3qMMRm7Xq21R4GjvvPmvwIRERGRIoQlpKWHJU8k2zG3Ba0T6PIFqkHGZ2z6x6O14UwqSD9Hi382p7U24QavKG4rm4iITI2DR97i4GtH85armzeXuprjp6BGItMrFCHNXbMsYYyJprduZZs0gBOk2oEefGuruZ8TaWWbgB3+HW7I6zPG1Pu6O73XZwuNIiJSJv/w2PPc+uCevOWuvXQVf9p8+hTUSGR6hSKkubpwwpTXddmCb0kOd+B/i7fArbU2ZozZnBbqNgKxDMEuSlrwclvN0l/fCvTnGNcmIiJl8rnzTqV59ZLk9t6Dr/OVLU/w7Y1nc1rdScn9dfPmZnq5SMUJTUiz1m42xrS74QygMW2gv/d0AP9TCLq8RXBdEeDSDKeP40xCSJf++oU5Jg2IiEgZ1dUcn7Eb87S6kzjrlPnTUCOR6RWakAZOUPNt9qcd68XXsubuS5Aa2rKdtz7L/oJeLyIiU2t0zPLkiwkAnnwxwZkn11BdNX2TtTo6nNWZFi5cSCQSoba2lpaWFjo6Ouju7i76vLFYjI6ODoaHh9mxY0f+FwQQj8fp6emhsbERgG3btgFQX1/PwMAAfX19k/p+MvmMtXa66zBjGWNqgJGRkRFqamqmuzoiIhXhgd/s58Yf72b/yFvJfSfPP54bPrmay886eUrrEovF2LRpE93d3TQ1NSX3x+NxOjo6iMViDA1NmJcWyODgIG1tbSWfJ119fT07duwgEokkt9va2mhvby85XAL09vbS2pr6YKCOjg7i8bgCYB5Hjhxh/vz5APOttUeylQvTY6FERGSWe+A3+7nm+7GUgAbw8shbXPP9GA/8Zv+U1mf9+vUTAhpANBqlrS3Xo6ULV1ubd7WpwAYHB4lEIsmABhCJRJLXUWpAAxgYGJiwr7m5mY0bN5Z8bnGEqrtTRERmr9Exy40/3k2m/h0LGODGH++mefXSKen69Lo40wOap6mpiWg02wNzpl85wp+nt7eXeHziQgjZvldSHIU0EREJhcefHZ7QguZngf0jb/H4s8OcX7+w7PXp7+/PGzq8IOfZvHlzMrjF43Ha28fnpvX29hKNRkkkEsTjcSKRSEp3YSwWIx6PE4/HefXVVwtq7fLeLx6PE41GaWlpYXBwkJ6eHrZv387mzePDrr0xavX19bS2thKJRDK+3l9fP6+ug4ODDAwMEI/Hk+dvb29Pjq+Lx+PJrtv+/n66urqSXaBNTU3E43Gam5uJRqP09PQkv1+56uI3ODiY/L7fcccdxOPx5Ji+np4eent7qa2tZcuWLXR2dtLQML6efSKRSN6HgYEB2traUo5nu0f57l3ZWGv1UeQHUAPYkZERKyIipfnRzhftio6f5P340c4Xp6Q+gO3u7i64fEtLix0YGEhuDw0N2aamJmuttX19fbanpyflmLe9Y8cOG4lEUl4bjUbtjh078r5fX19fcrupqSn5moGBAdvQ0JBSvqGhIeU9cr2+u7vbtre3J4/19fWllM10fu9aotFoyr6BgYEJ+9K/r7nqkol3zvTvWXqd0+vY3t5uh4aGUl5z+PDhZPlM9yjXvSvWyMiIxfm7o8bmyBkakyYiIqFQN6+wpwgUWm4qxWIxBgcHU1reotEow8PDDA46S3f29fWRSCSSx9auXZssm0gkJrw2U3eiJx6P09/fn9LatH79enp6egqqb67XJxIJOjo66OzsTB7bsmVLzvrk0tTUxPDwMLHY+IN8/GPlirmW2tpa4vH4hO+ZX0NDw4Q6x+Px5P3wXuPfznaPct27clJ3p4iIhMIHV9Zy8vzjeXnkrYzj0gywdP7xfHBl+cZa+UWj0bwzLr2uue3bt2ccn+Z1q3V3d9PT08OCBQtoaGhg48aNKV2h6a+NRCIMD2da3tPhTQzwB4yhoaGCg1Su12/fvn3CpINSZ2u2trbS09NDT08Pg4ODbNiwoeRryfQ9q6/PuOJWkncdXrfl8PBw8vvc0tKS8R41NDTkvHflpJAmIiKhUF1luOGTq7nm+zEMpAQ1b5rADZ9cPWXrpbW0tNDf35+zzODgIK2trclWllwGBgaSLW5eK1Gxv+wTiQTRaDSlJSnIoP1cr893zZl4YTWbtrY2zj33XHp6eia0gJV6LUHEYjG6urpobm5mw4YNE+qc7R5N5r0LQt2dIiISGpefdTK3f76BpfNTuzSXzj+e2z/fMKXrpHkD9/0tPH6JRCI5g9IbEJ8uHo/T2NiYHITf0NBAe3s7O3bsYMuWLUXXLVNXnlenUl/f0NCQ8Ty5zu3vyswkGo1SW1tLf3//hFmnpV5LoRKJBJdeeimdnZ3JiRPee8Tj8az3aLLvXRAKaZPg1f3Ps/dXj+T9OLTvuemuqohI6F1+1sk80vExvvHpswD4xqfP4pGOj035QrbgdI91dHRMCGreLEFvHFVDQwNNTU0p5bzg0tLSkizvl6vlKV9AaWpqYu3atRNave67776815Tv9d7MSv/M0EQikXJu/5i5eDyeMkMym7a2NjZt2jShlazUa/HXMZd4PJ4MoR6vqzMWi2W9R0Hv3WTSEwdK4D1xYOC2/0LToe/lLf/o8k2c/8Wby18xEZEK8JuXRvjEbY/wky9dOO3P7kx/LBSQcQmGjo6O5LiooaGhZGuc90vea0WKx+O0trYSj8fp6uqiv7+f7u5u2tvb2bx5M11dXUSjUTo7O7MuReF/P++83hIc3d3dbN++nc7OzpRzrl27lvXr1yfrnun1/nMvXLgwOQEi09MFgOSSHl5Xov9a/LwJCdkmBOSqi1+m9/FfX0dHB7W1tcky7e3tdHZ2EolEknVubm4GnLDV0dHBxo0bk4Et/R55YTF9v3/MXlCFPnFAIa0EXkiLP/1rRt9KJPcffm4XjbF2tjVsZsGKNcn9kcXLWbRsxdRXVERkBjh45C0OvnY0ub334Ot8ZcsTfHvj2ZxWd1Jyf928uRkfxC4yUxQa0jRxYBIsPPlUamrOSm7vBYjBghVrOO0DF05bvUREZpJ/eOx5bn1wz4T9X9nyRMr2tZeu4k+
2021-10-11 18:31:02 +01:00
"text/plain": [
"<Figure size 640x395.55 with 1 Axes>"
]
},
"metadata": {
"needs_background": "light"
},
"output_type": "display_data"
}
],
"source": [
"periodic_m_eff.show([4,47], comp=m_eff, ylabel=r'$am_\\mathrm{eff}$')"
]
},
{
"cell_type": "markdown",
"id": "e3762e68",
"metadata": {},
"source": [
"Arithmetic operations and mathematical functions are also overloaded for the `Corr` class. We can compute the difference between the two variants of the effective mass as follows."
]
},
{
"cell_type": "code",
"execution_count": 13,
"id": "e56d164c",
"metadata": {},
"outputs": [
{
"data": {
2022-01-06 11:16:00 +01:00
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAlIAAAGNCAYAAADaX58UAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjUuMCwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8/fFQqAAAACXBIWXMAAA9hAAAPYQGoP6dpAAArOklEQVR4nO3df2xd5Z3n8c+TlMQwse/NjzExxIE4tGVMujAmToeoaFaDzZLRZItGSYgEXQmtSMioq261VdJUjLKIKplEna7QVGXi/tlk1SHMCBQ0YcG0i6iaBqgHaUM2DNhAaOqMicO9N2mww9rf/eOe61zf3B/nHN8f5977fkkWueecnPswZ5p8eJ7v8z3OzAQAAIDg5tV6AAAAAPWKIAUAABASQQoAACAkghQAAEBIBCkAAICQCFIAAAAhEaQAAABC+kKtB+CHcy4uaZskmdmB2o4GAAAgrV5mpPokLa31IAAAALLVxYyUmT3nnFsiKR709zrnnKSbJF0s97gAAEBDa5X0OyvyGpiqBSlveW6LpM1m1p/n/E5JCe9jvIxLeDdJ+m2Z7gUAAJrLCklnC52sSpByzvVIWqv0jNKSPOd3SpKZDXif+5xzB81sexm+/qIkffzxx2prayvD7QAAQKNLpVLq7OyUSqxoVSVImdmQpCHn3KYCl+yWtCrr+kHn3CuSyhGkJEltbW0EKQAAUFY1LzZ3znUpvZSXyHOuL8T9Fjrn2jI/Sq9vAgAAlF3Ng5SkrgLHE/KKy71A1S+pv8isVsZuScmsH+qjAABARUR5194FefVUZjYoadDn79sn6YdZn1tFmAIAABUQ5SB1TVG6H2Y2KWky8znd/QAAAKD8orC0N1LgeLzIOQAAgJqreZAysxFJCa/oPPec3+U8AACAqqt2kCq0XLdP6dfASJK8gvKBqowIAAAgpKoEKedcl9d0c7ukHufc/uzdd14X87hzbpN3vLdMzTgBAAAqxhV5fUxD8HpJJZPJJA05AQCAL6lUSrFYTJJiZpYqdF3Na6QAAADqFUEKAAAgJIIUAABASAQpAACAkAhSAAAAIRGkAAAAQiJIAQAAhESQAgAACIkgBQAAEBJBCgAAICSCFAAAQEgEKQAAgJAIUgAAACERpAAAAEIiSAEAAIREkAIAAAiJIAUAABASQQoAACAkghQAAEBIBCkAAICQvlDrAQAAANTC1LTpjQ8uaOzihNpbW7Ru1RLNn+cC3YMgBQAAms5LJ0f15NFTGk1OzBzriLVoz8ZuPbCmw/d9nJlVYnyR4Zxrk5RMJpNqa2ur9XAAAEAF+ZlleunkqHYcGlJuAspc9cwjPVq/8g8Ui8UkKWZmqULfx4wUAABoCH5mmUYTn+mJ509eE6IkzRx74vmT+p//6d/5+k6KzQEAQKRNTZuOD4/rhbfP6vjwuKamr41BmVmm7BAlSeeSE9pxaEgvnRyVJD314imdv3Sl6Pedv3RFB/7XaV9jY0YKAABElp9Zpqlp05NHTxWcZXKSnjx6Sv3dy/XZ51O+vtfvdcxIAQCASPI7y3TkrTPXXJPNJI0mJ3TkrTPasrbT13f/5R+v8HUdQQoAAFRdqeW6UrNMUnqWaWra9JPXP/D1nT95/QPdf8dydcRaVKjJgVN6xuvP/uhGX/dkaQ8AAFSVn+W6ILNMj927St/9p5Mlv/exe1dp/jynPRu7tePQkJw0K6hlwtWejd2++0kxIwUAAKrG73JdkFmmzWtX+ppl2rx2pSTpgTUdeuaRHi2Ptcy6bnmsRc880hOojxQzUgAAoCxK9XAKUhRe6VmmB9Z0qL97OZ3NAQBA7ZV7uW7z2pV6+tX3dS45kTd4OaVnkHJnmXLHsLxIt/L585zuWb001L/vzDjobA4AAObCT6fwB9Z06L6//d8a/uT3Je+3+g//QK/+t38/c18p/yxTvmW4crw/T5JSqZSvzubUSAEAgILKubvusXtX+frOzHVhapkys0xfv+tm3bN6aagQFQRLewAAIK9aL9dJ5atlqhRmpAAAwDUqsbsuUxQu6ZoddsVaD1R7likIghQAAJglast1UcbSHgAATaZUQfbL75zztVz38jvnGnK5LgiCFAAATcRP3dOzb33s617PvvWxNnylI1Sn8HK0HogClvYAAGgSfuuerr9uvq/7Za5rtOW6IOgjBQBAAyi1XDea+Ewbf/RLnb90peA9li1aoKPf/Jok+b62I3697zHUE799pOpiac85t0nSEkl3SzpiZoM1HhIAAJHhZ7nuqRdPFQ1GknT+0hU99eIp/fiRu/X9B9cUbYb5/QfXzApRUuMs1wUR+aU951yPJJnZgKRdko7UdkQAAESH3+W6zz6f8nW/zHXNvFwXRD3MSC2R1C/pOTNLOOcuOOd6zGyo1gMDAKCWgrwEeMvaTv3i3U9K3nPL2s6ZXzfS7rpKKUuQcs7FJW2RtNnM+vOc3ykp4X2Mm9kBv/f2lvGyl/KWEKIAAM2gVM1R0K7iHbGWkm0K7r9j+azjzbhcF8Scg5S39LZWUlzp2aPc8zulmaU5Oef6nHMHzWx7iO86KOmxOQ0YAIA64KfuKUhX8a3rbgnVpgDFlW3XnlcQvtvM7s45/qmkVWaWyDpmZua8X++UlC/qjmfPXHn3l5k9F3Bc7NoDANSVTN1T7t/QmYiTqVH62Rsf6bv/dLLk/f7mL9do67pbZu5dKqAhIrv2nHNdSi/lJfKc6zOzQT/LfM65PkkJMxv0ZsASZjZS/hEDAFBZftoUPPH8yaKvZ3ni+ZO6c0W86buKR0Gli827ChxPKL0UWJIXxo54v5bSwazg03bOLZS0MOtQq5/vAQCg0irRpqCZu4pHQa3aH1xQnnqqfMxsxMwWZ/2Uisy7JSWzfn47t6ECADB3tCloTLVqf+ArRIW0T9IPsz63ijAFAKgh2hQ0rkoHqUJ1TPEi5+bEzCYlTWY+e8uBAABUTKm6p5ffOeerTcHL75zT/Xcsp01BHalokDKzEedcwjnXlVsczmteAACNwE/d07NvfezrXs++9bE2fKWDNgV1pJw1UoWW6/ZJ6st88NoYDJTxewEAqAm/dU/XXzff1/0y11H3VD/m3EfK21W3SdJDknokHZD0Zna/J69XVGZGqtfMds3pS4ONjz5SAIDA/LQp2PijXxbdYbds0QId/ebXJMn3tdkvAi41BlSO3z5SZWvIGVUEKQBAUH6W6/7q0G/0zyfPlbzXn69Zrh8/cvfM7JWUf7mOmaZo8RukatX+AACASKJNAYKoVfsDAAAihzYFCIogBQCAhzYFCIqlPQBA05iaNh0fHtcLb5/V8eFxTU3PjkBB2hTMn+e0Z2O3pKt1Thm0KWgezEgBAJqCnwLysG0Kcu+7POe+aFzs2gMANLxMAXnu33i5O+aCtDSgTUFjY9ceAABK93t64vmTBQvITdITz5/UaOIzdcSv1/cfXCOn/Mt1TtL3H1wzK0RJV+uevn7Xzbpn9VJCVBMhSAEA6lqpuqenXjxVdIZJks5fuqKnXjwliTYFCIYaKQBA3fJT9xS035NEmwL4x4wUAKAu+W2cmd3HqZjc61iugx8EKQBA3QlS95Tp91QoBjmlZ7Fy+z0BfhCkAACRU866J/o9oZKokQIAREol6p7o94RKoY8UACAy/PZ7OvZ/RrXj8FDJ+z3zcI82fOVqSKLfE/zy20eKGSkAQFWUCjGl6p6kdN3TnSvivOcOkUGQAgBUnJ/luiB1Tz9+5G7t2ditHYeG5KRZYYq6J1QTxeYAgIry26YgbN0TjTNRS8xIAQBCK+dy3Za1nfrFu5+U/M7sfk80zkStEaQAAKH4Wa77wcvv+lqu+8HL7+rApjupe0LdYWkPABCY3+W61X+4yNf9Vv/hIvo9oS4RpAAAs5Rqhjk1bXry6Kmiy3VPHj2lqWnT8CeXfH1n5jrqnlBvWNoDAMzws1z38jvnrpmJymaSRpMTevmdc/rO/V/Wa//6SdHlvWWLFug793955jN1T6gnzEgBACT5X6579q2Pfd3v2bc+Vkf8en3/wTVyyr9c5yR9/8E16ohfP+scLwxGvSBIAQACvQT4+uvm+7pn5jqW69D
2021-10-11 18:31:02 +01:00
"text/plain": [
"<Figure size 640x395.55 with 1 Axes>"
]
},
"metadata": {
"needs_background": "light"
},
"output_type": "display_data"
}
],
"source": [
"difference_m_eff = np.abs(periodic_m_eff - m_eff)\n",
"difference_m_eff.show([0, 47], logscale=True)"
]
},
{
"cell_type": "markdown",
"id": "472ab97b",
"metadata": {},
"source": [
"## Derivatives"
]
},
{
"cell_type": "markdown",
"id": "d99414fe",
"metadata": {},
"source": [
"We can obtain derivatives of correlators in the following way"
]
},
{
"cell_type": "code",
"execution_count": 14,
"id": "03007f8a",
"metadata": {},
"outputs": [],
"source": [
"first_derivative = symmetrised_correlator.deriv()\n",
"first_derivative.tag = 'First derivative'"
]
},
{
"cell_type": "code",
"execution_count": 15,
"id": "c0311739",
"metadata": {},
"outputs": [],
"source": [
"second_derivative = symmetrised_correlator.second_deriv()\n",
"second_derivative.tag = 'Second derivative'"
]
},
{
"cell_type": "code",
"execution_count": 16,
"id": "165550d9",
"metadata": {},
"outputs": [
{
"data": {
2022-01-06 11:16:00 +01:00
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAlMAAAGLCAYAAADjxBc3AAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjUuMCwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8/fFQqAAAACXBIWXMAAA9hAAAPYQGoP6dpAABHXUlEQVR4nO3df3hcZZ3//9fdsrTWNpmmpbQWrEwtai2C09ZVL5GVThBEXNxP2roK7rULJB+UtbuVnRjd/fhRd62TrVyLen3dSUHdtbq0yUcRVJAZuiJ7IdImdAErS80A0looJZ20BQua3N8/zpzpzGR+ZjIzJzPPx3XNlZxz3zPnfU7Tk3fuX8dYawUAAIDJmVHvAAAAAKYzkikAAIAKkEwBAABUgGQKAACgAiRTAAAAFSCZAgAAqADJFAAAQAVOq3cA1WaMMZJeI+l4vWMBAADTyjxJv7VFFuVs+GRKTiJ1oN5BAACAaeksSQcLVWiGZOq4JD3zzDNqaWmpdywAAGAaOHbsmM4++2yphJ6tZkimJEktLS0kUwAAYMoxAB0AAKACJFMAAAAVIJkCAACoAMkUAABABUimAAAAKkAyBQAAUAGSKQAAgAqQTAEAAFSAZAoAAKACTbMCOgA0osPHTurw8ZeL1ls0b5YWtcyuQURA8yGZyiH95jQ2bvXL347q6Eu/1/w5f6Q3v6ZVM2cYSdycANTfd37xG9187/6i9TatW6G/bT+3BhEBzafiZMoY45O0QdJ6a217jvJQ8tvlkmSt7cpRnkhu+qy1veWUVwM3JwDTxUf++LVqX3lmavvXh0/ob3bs1b9svECvXzQ3tX/RvFn1CA9oChUlU8aYgKQ1knyS2nKUh6213WnbEWNM1E263ETLWtuX3A4aYyJuwlWsvFo+8sev1atPn6kv3vV43jqfvuyNuvKtS6sZBgAUtahlds4W8tcvmqtVS1vrEJHU23vqb94XXnhBXV1dGhgYUCgUKvAuFNLV1aWdO3eqv79fwWCw3uEUNTQ0pO7ubsXjcQ0PD9c7nKqraAC6tXYomejEs8uSLVaB5FdXRFLQGONPbvdI6kv7vJikzrT6xcqrYsHcWfrmA0/lLTeSvvnAU1owl7/0AHjH2LjVIwcSkqRHDiQ0Nm5rHkNXV5eCwaBCoZBCoZDC4bC6u7u1e/fumscyVfr6+opXSuru7tb69eunPIZIJCK/31+8okcEAgF1d3cXr5hDOdfbK6o9m2+NpPR/fTfp8iUTKp+1NpH9pmQLVMHyfAc0xswyxrS4L0nzyg36oSdHdGj0ZN5yK+nQ6Ek99ORIuR8NAFVx92OH9K7wLn36+49Jkj79/cf0rvAu3f3YoZrGsXPnTgUCgYx927Ztq2kMUy0ajZZct729XRs3bqxiNNNHW9uEDquSlHO9vaJqyZS1NmGtnW+tHUrb7SZBcWUmWekScroNi5Xn0yNpNO11oKSA0xw+nj+Rmkw9AKimux87pOu3D034I/DZ0ZO6fvtQTROqRCKheDyzs8Ln82nt2rU1i2Eq9fX1TTifQoLBoDo6OqoYUWMr93p7Ra1n8/VI6rLWJowx+eqMyBl/lShSns8WSTelbc9TmQnVonmlzdArtR4AVMvYuNXn7tynXB16Vs6whM/duU/tKxenZiJXUyAQUHt7uyKRSMbYHne81MDAgLZs2aJ4PJ4a/xOPx9Xe3i6/369IJKJ4PJ7qItq2bZvi8bhGRkY0ODioSCSivr4+tbW1aceOHerp6Um1hMVisUm9T3LGefn9fsXjcfn9fnV0dCgWiykajSoej6fGgYVCodRx/H6/urq6Ui0pGzduzDlOqK+vT36/P5Vo+nw+dXZ25j1uekw+n29SLTzZXWXu8dKPKUnxeDz1b5PvvNrb23PuD4fDRc8hWyKRSF2PaDSqrq6ujH+/XNd7snG78dWEtbbil6QOSYNF6oQlhdK2g87hJ9Q7KmdcVMHyMmJrkWRHR0dtqX579CW7+gv32GXdP8z7Wv2Fe+xvj75U8mcCQDU88OsjBe9V7uuBXx+pSTzDw8PW7/dbObmcDQaDNhqNZtSJRqPW7/dn7AuHwznrpL/X7/fbUCiU2u7v77eBQKDi93V0dNj+/v7UdjAYtIODg6nPyz5G+mdEo1E7ODiY+vzBwcGMc+vv77eRSCTj+rjbhY4bCoUy3nf06FEracK1zCUcDk84X/c4HR0dGZ8xPDxsg8Fg0fPKt7/QOWRfC/e8hoeHU9t+v98ePXo0tZ3vek827kqMjo66P8cttkiuUZMV0I0xHZKGbeayBvna8XzJsmLlVXPb7md05MQrBescOfGKbtv9TDXDAICivDYswe/3a3h4WNFoVKFQSCMjI2pvb9fAwECqTjAY1MjIiIaGhjLel66trU3xeDyjdSu7TiAQmNAlVO774vG4BgYGMlpT1q9fr0gkUvA8fT6fhoaGFAwGFQgECraC9Pf3K5FIpGJZs2ZNweMmEgn19vZmtCb5fL4JY9FySSQS6u7uVk9PT2rfjh07FI/HNTQ0pFgsNuHajIyMKBaLFTyvXPsnc+3i8XjqWO7x07dzqSTuWql6N587WNyeWt7AJ6nNWhs3xiSMMX5rbcb/BuvM2lOx8mpx12154NdH1Hd/PCOxWjj3dHVe6Nc7X7+QdVsA1J1XhyUEg8HUL7/u7m5dd911Gb90Ozs7tWPHDgUCgQm/KF3ZSZDP59Py5cuLHruc98ViMfl8voxf6MPDwyWN2ylldl1HR4cikYjmz5+vQCCgjRs3KhQKqa+vL+9x3ZgmY8+ePfL5fBnv7+/vl3SquzHXeUSj0dS/Qb7zyt4/mWvnxuJ2eY6MjGhkpPBkrj179lQUdy1MVTKVs0M3uQ5VQNJA2nIIHTq13MEWOd15bqKVXlZKeVW467asWtqqay7066EnR3T4+EktmjdbbzunrSbjDgCgFG87p01LWmfr2dGTOcdNGUmLW517V7UlEgnFYrEJY2bC4bB6e3uVSCRSv+Q3btyodevWpVo46rV2UiKRkN/vzzh+oVjccUGSSk54otFoqnUlvdUm33EHBgYmPRPObQErtyxdvvPK3l/utZOcVqYtW7aovb1dGzZsKJoAxePxiuOuhYq6+Ywx/uTCml1y1pQKJxMetwXqXjljpYbTXmGbXO4g2e3nM8Z0JN+31qYtyFmsvBZmzjB6x/IF+tMLluodyxeQSAHwlJkzjD57xUpJTuKUzt3+7BUra3bvyreelN/vz/hlFwgE1NbWVlHiMBVydRVK+ROP9K7JUrgDwQOBgEKhkAYHB1MtcvmOm6+sFIFAIGfsiUQiNdg/Wzwen9Rsy3KvXSKR0Lp169TT06POzk75fL5U3Xzn63bdTWXc1VDpop1xa22vtXa1tdZYa7uttQPJMndpBJP9yvqMXmvtQPI1YYWvYuUA0OwuXbVEX78qoMWtmV15i1tn6+tXBXTpqiU1i6Wvr2/CGJhcrVWSs8BndvdfIaW2UJTzvmAwqDVr1mSM6ZKc9bIkpWapSc4v71LGLWUfO3tmnduak++4fr9fnZ2dGe9LJBIaGhoqeg3c2XTpq9AnEonU+l/BYDDj38dNDieznEOxa5fNbWVKv4ZuF58bR67rPdVxV0WxEerT/aVJzOYDgOnoD2Pj9jsPPmWXdf/QfufBp+wfxsZrevyjR4/aSCRio9GoDYfDGa9choeHbUdHx4T9g4ODtqOjw0pKvTccDlufz5eaHZheJxQK2aNHj076fdaemj2XPvPNFQqFMmbXRaNRGwwGrc/ns+FwODU7LdfxI5FIxueGw+GM2WuFjhsOh21/f7+NRqOpWW7ujLViQqFQ6v3pswLTjxmJRDJmveU7r3z7C51D9nXOvpbuObk/A+nnnn29Jxt3pcqZzWesrf3jBmopuQr66OjoqFpaWuodDgBU1WMHR/X+r/6XfvjX76rbs/lKlUgktGfPnmnxrDk0n2PHjqm1tVWSWq21xwrVrfWinQCAKXT42EkdPv5yavvXh09kfHUtmjcr5wORay19EPrOnTszpv8D0xXJFABMY9/5xW908737J+z/mx17M7Y3rVuhv20/t0ZR5dfd3a3
2021-10-11 18:31:02 +01:00
"text/plain": [
"<Figure size 640x395.55 with 1 Axes>"
]
},
"metadata": {
"needs_background": "light"
},
"output_type": "display_data"
}
],
"source": [
"symmetrised_correlator.show([5, 20], comp=[first_derivative, second_derivative], y_range=[-500, 1300])"
]
},
2022-01-06 17:38:12 +01:00
{
"cell_type": "markdown",
2022-01-06 17:51:50 +01:00
"id": "7fcbcac4",
2022-01-06 17:38:12 +01:00
"metadata": {},
"source": [
"There is a range of addtional methods of the `Corr` class which can be found in the documentation."
]
},
2021-10-11 18:31:02 +01:00
{
"cell_type": "code",
"execution_count": null,
2022-01-06 17:51:50 +01:00
"id": "2fbe1263",
2021-10-11 18:31:02 +01:00
"metadata": {},
"outputs": [],
"source": []
}
],
"metadata": {
"kernelspec": {
2022-01-06 11:16:00 +01:00
"display_name": "Python 3 (ipykernel)",
2021-10-11 18:31:02 +01:00
"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",
2022-01-06 11:16:00 +01:00
"version": "3.8.10"
2021-10-11 18:31:02 +01:00
}
},
"nbformat": 4,
"nbformat_minor": 5
}
