{ "cells": [ { "cell_type": "markdown", "id": "2f1f53e2-8b15-4e78-8428-500ff3554a72", "metadata": {}, "source": [ "# # Time dependence - Gaussian pulse" ] }, { "cell_type": "code", "execution_count": 1, "id": "012b8646-6ff7-40f8-9858-37257b7d639b", "metadata": {}, "outputs": [], "source": [ "import nexus as nx\n", "import numpy as np\n", "import matplotlib.pyplot as plt" ] }, { "cell_type": "code", "execution_count": 2, "id": "df308529-2df4-4ec3-a957-f582f15d4aaa", "metadata": {}, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "Var.value = 20.0, .min = 10.0, .max = 30.0, .fit: True, .id: t0\n", "Var.value = 5.0, .min = 1.0, .max = 10.0, .fit: True, .id: fwhm\n" ] }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAjcAAAGwCAYAAABVdURTAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjguMSwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/SrBM8AAAACXBIWXMAAA9hAAAPYQGoP6dpAAA+pUlEQVR4nO3de3xU1b3///dkkpkk5AoxCZdgwAuIXAuFRutpranU+qUHL/1RSoVSSx+2oGhqK1QuXqqxeqBoD0q1BXz01IPVg54epXhoBJVjCspFpQoqF0EkCRAgN8hlZv/+SPaESMBksmf2zuzX8/HIo2TP3nvW7Ad13qz1WWt5DMMwBAAAECPi7G4AAACAlQg3AAAgphBuAABATCHcAACAmEK4AQAAMYVwAwAAYgrhBgAAxJR4uxsQbcFgUJ999plSU1Pl8Xjsbg4AAOgAwzBUXV2tPn36KC7u3H0zrgs3n332mfLy8uxuBgAACMOBAwfUr1+/c57junCTmpoqqfnhpKWl2dwaAADQEVVVVcrLywt9j5+L68KNORSVlpZGuAEAoJvpSEkJBcUAACCmEG4AAEBMIdwAAICYQrgBAAAxhXADAABiCuEGAADEFMINAACIKYQbAAAQUwg3AAAgphBuAABATLE13Lz++uuaMGGC+vTpI4/HoxdffPELr9mwYYO+9KUvye/368ILL9TKlSsj3k4AANB92BpuamtrNWLECC1durRD5+/du1fXXnutrrzySm3fvl233367fvzjH+uVV16JcEsBAEB3YevGmddcc42uueaaDp+/bNkyDRgwQIsWLZIkXXLJJdq4caN++9vfavz48ZFqJhymoSkoSfLFM6oKADhTt/p2KC0tVWFhYZtj48ePV2lp6Vmvqa+vV1VVVZsfdF8fV9Toiodf1Vd/86o+Kq+2uzkAAAfqVuGmrKxMOTk5bY7l5OSoqqpKJ0+ebPea4uJipaenh37y8vKi0VREyJK/f6jyqnpVVNdrSclHdjcHAOBA3SrchGPu3Lk6ceJE6OfAgQN2NwlhOtUY0Lr3y0O/r3u/XKcaAza2CADgRLbW3HRWbm6uysvL2xwrLy9XWlqakpKS2r3G7/fL7/dHo3mIsHc/PaH6pqCyUnySPDpSU6/tB47rKwN72d00AICDdKuem4KCApWUlLQ5tm7dOhUUFNjUIkTTu58elySNPj9TX87PlCS9c+C4fQ0CADiSreGmpqZG27dv1/bt2yU1T/Xevn279u/fL6l5SGnq1Kmh82+55Rbt2bNHv/zlL7Vz5049/vjj+stf/qI77rjDjuYjynYfrpUkXZSdqkt6p0mSdpVRVAwAaMvWYam3335bV155Zej3oqIiSdK0adO0cuVKHTp0KBR0JGnAgAF6+eWXdccdd+jRRx9Vv3799Ic//IFp4C6x53CNJOmC7B5K9jX/1d1JuAEAfI6t4ebrX/+6DMM46+vtrT789a9/Xdu2bYtgq+BUe44099wMzEpRD79XknSgsk6GYcjj8djZNACAg3SrgmK4V9WpRh2urpckDTyvh+LjmkdUq+ubVHWySenJCXY2DwDgIN2qoBjutael3iY71a/UxAQl+bwts6akA8fq7GwaAMBhCDfoFg5UNgeY83slh471zWz+88Hj7S/gCABwJ8INuoXyqlOSpNz01vWM+mU2//nTY4QbAEArwg26hbITLeEmrXVBxtZww7AUAKAV4QbdQllLz01OWmLoWL+WYakDlfTcAABaEW7QLbQOS7WGmz4tf66oPmVLmwAAzkS4Qbdg9tzkntZzk5XSPERlThEHAEAi3KAbMAxD5VXNAeb0Yams1OZwc6Sm/pyLQQIA3IVwA8c7XteohqagJCn7tIJic52bxoChEycbbWkbAMB5CDdwPHNIqmcPn/zx3tBxf7xXaYnNi2wfqWFoCgDQjHADxzta0yCptafmdObQVAV1NwCAFoQbON7R2ubg0rPHmeHmvBSz7qYhqm0CADgX4QaOV1nbHFx69fCf8VqoqJieGwBAC8INHO9YS7jJ7HHmzt9mz81ham4AAC0IN3C8oy3hpmc7PTfnpbLWDQCgLcINHK91WOrMmhuzDsfs3QEAgHADx2vtuTkz3GQmNw9VHasj3AAAmhFu4Hjn6rnJSG4+dryORfwAAM0IN3A8M9z0bGedm0wz3LBCMQCgBeEGjhYIGqEhp/aGpTJahqWO1zUoGGR/KQAA4QYOd+Jko8w9Mc1emtOZ4SZoSNWnmqLZNACAQxFu4GjHW3ptUv3xSvCe+dfVH+9Vsq95vymKigEAEuEGDmfu9p2WdOYCfiazR4dwAwCQCDdwuKqWoaZzhZvWuhuKigEAhBs4nNlzk54Uf9ZzWmdM0XMDACDcwOFaw80X99wcq6XnBgBAuIHDVZk1N4lfXHNznJobAIAIN3C4qk703LCQHwBAItzA4ToyLGX26rDODQBAItzA4apOffFU8NTE5mLjKnpuAAAi3MDhOtRzk0TPDQCgFeEGjtaRcBPquTlFzw0AgHADh6s6aS7id/Z1bsyaG4alAAAS4QYOx7AUAKCzCDdwrGDQ6FRBcXV9kwJBIyptAwA4F+EGjlVd3ySjJaucaxE/M9xIUk09vTcA4HaEGziWWUPjj49TYoL3rOf5471KTIhrcw0AwL0IN3CsjtTbmFLNomJmTAGA6xFu4Fgd2XrBlGbW3VBUDACuR7iBY3WkmNhknsOwFACAcAPHCmdYip4bAADhBo4VWsAv8ewL+JnSWKUYANCCcAPHqm6Z1p16jmngplBB8Ul6bgDA7Qg3cKyaU2a46UDPTZJZUEzPDQC4HeEGjlVT3xxUUjo0LMVUcABAM8INHMssDk71d6LmhmEpAHA9wg0cy9xKoUM9N+bmmfX03ACA2xFu4Fhmz02KvyOL+FFQDABoRriBY4V6bjowLBXaGZyaGwBwPcINHKszs6Va95ai5wYA3I5wA8fqTM+NWZdjXgMAcC/CDRwpGDQ6VVCc4ms+p6EpqIamYETbBgBwNsINHKm2obUHpiM9Nz383tZr6b0BAFcj3MCRzF6bBK9H/vgv/msa741TYkJcm2sBAO5EuIEj1ZxqrbfxeDwdusbs4Tm91wcA4D6EGziSOeupI/U2ph5muKHnBgBczfZws3TpUuXn5ysxMVHjxo3T5s2bz3n+kiVLNGjQICUlJSkvL0933HGHTp06FaXWIlpaZ0p98QJ+ph4+c60bwg0AuJmt4ebZZ59VUVGRFi5cqK1bt2rEiBEaP368Kioq2j3/mWee0Zw5c7Rw4UJ98MEH+uMf/6hnn31Wv/rVr6LcckRaTSf2lTKZvTy19YGItAkA0D3YGm4WL16sGTNmaPr06RoyZIiWLVum5ORkLV++vN3z33zzTV1++eX6/ve/r/z8fF199dWaPHnyOXt76uvrVVVV1eYHzteZHcFNKQxLAQBkY7hpaGjQli1bVFhY2NqYuDgVFhaqtLS03Wsuu+wybdmyJRRm9uzZozVr1ujb3/72Wd+nuLhY6enpoZ+8vDxrPwgioroTqxObzJqbasINALhax785LHbkyBEFAgHl5OS0OZ6Tk6OdO3e2e833v/99HTlyRF/96ldlGIaampp0yy23nHNYau7cuSoqKgr9XlVVRcDpBjqzOrEppWWtG3puAMDdbC8o7owNGzbowQcf1OOPP66tW7dq9erVevnll3X//fef9Rq/36+0tLQ2P3C+mjBmSzEsBQCQbOy5ycrKktfrVXl5eZvj5eXlys3Nbfea+fPn66abbtKPf/xjSdKwYcNUW1urn/zkJ7r77rsVF9etshrOwey56UxBMcNSAADJxp4bn8+n0aNHq6SkJHQsGAyqpKREBQUF7V5TV1d3RoDxepuHIgzDiFxjEXXVYQ1L0XMDALCx50aSioqKNG3aNI0ZM0Zjx47VkiVLVFtbq+nTp0uSpk6dqr59+6q4uFiSNGHCBC1evFijRo3SuHHj9PHHH2v+/PmaMGFCKOQgNrQOS3V8nRvCDQBAsjncTJo0SYcPH9aCBQtUVlamkSNHau3ataEi4/3797fpqZk3b548Ho/mzZungwcP6rzzztOECRP0wAMP2PURECHVp1qmgoczLMUifgDgaraGG0maNWuWZs2a1e5rGzZsaPN7fHy8Fi5cqIULF0ahZbCTuRBfWMNS7C0FAK5GBS4cyQwoPfwdH25s3VuKFYoBwM0IN3Aks26mRxg9NzXU3ACAqxFu4Ei1Dc29L2GFG2puAMDVCDdwnMZAUA1NQUlSD19nhqWazz3ZGFAgyNIAAOBWhBs4Tt1pNTPJvs7vCi5RVAwAbka4geOYwcTnjZMvvuN/Rf3xXiV4PZIYmgIANyPcwHHMYuLkTsyUMvVgIT8AcD3CDRwnVEzciSEpk3kNM6YAwL0IN3Cc1mngne+5SU0k3ACA2xFu4DihYalwem4YlgIA1yPcwHHMguLObL1g6hFayI9VigHArQg3cBxz+4TkTqxxY0oNLeTXaGmbAADdB+EGjlPXpZ6b5kBkFiUDANyHcAPHMYeUujIVnIJiAHAvwg0cp86cLRVGQXEKBcUA4HqEGziOWVDcmU0zTWyeCQAg3MBxulJQzLAUAIBwA8fpSkFxaFiKjTMBwLUIN3CcmtDeUqxzAwDoPMINHKeuZRp3SjizpVqGsuoYlgIA1yLcwHFq2H4BANAFhBs4Tl19F3YFD9XcMCwFAG5FuIHjtE4FD2e2VMsKxfVNMgzD0nYBALoHwg0cxTCM0JBSOOvcmNc0BQ3VNwUtbRsAoHsg3MBR6puCCrZ0uIQTbpITWnt76hiaAgBXItzAUU5ffC8pofPDUvHeOCUmNP+1pqgYANyJcANHMYuJkxK88sZ5wrpHCqsUA4CrEW7gKF3ZV8pkTiGvY5ViAHAlwg0cpbWYuPNDUiZWKQYAdyPcwFHM9WnCWePGxCrFAOBuhBs4irU9N4QbAHAjwg0cpbYLWy+YUtiCAQBcjXADRzEDSUqXCopbVilmnRsAcCXCDRzFDCRmQAkHm2cCgLsRbuAodRZMBTfrdVihGADciXADR6k1dwSnoBgAECbCDRyFgmIAQFcRbuAo5grFXSsobgk3DEsBgCsRbuAo5rBUVwqKU1qGtOi5AQB3ItzAUeos6LlhthQAuBvhBo5i7geVbMmwFOEGANyIcANHCU0F79KwlNlzQ80NALgR4QaO0rq3lAUrFDMsBQCuRLiBo4TWubFgKnh9U1BNgaAl7QIAdB+EGzhGIGjoZKN1i/hJDE0BgBsRbuAYdacVAHdlWMoXH6cEr0cSRcUA4EaEGziGuReUN84jf3zX/moyHRwA3ItwA8eoCW294JXH4+nSvXqwSjEAuBbhBo5RZ0ExsakHqxQDgGsRbuAYZn1MV4qJTewMDgDuRbiBY1ixxo3J7P2po6AYAFyHcAPHMOtjrByWqmEqOAC4DuEGjtHac2PdsFQdw1IA4DqEGzhGbWi2lHXDUhQUA4D7EG7gGOY6N9YWFDMsBQBuQ7iBY4SGpSzpuWkOSBQUA4D72B5uli5dqvz8fCUmJmrcuHHavHnzOc8/fvy4Zs6cqd69e8vv9+viiy/WmjVrotRaRFLrVHArCoqZCg4AbtX1b5EuePbZZ1VUVKRly5Zp3LhxWrJkicaPH69du3YpOzv7jPMbGhr0zW9+U9nZ2Xr++efVt29fffLJJ8rIyIh+42G50CJ+FgxLpbD9AgC4lq3hZvHixZoxY4amT58uSVq2bJlefvllLV++XHPmzDnj/OXLl6uyslJvvvmmEhISJEn5+fnnfI/6+nrV19eHfq+qqrLuA8BSNRYWFCebKxSz/QIAuI5tw1INDQ3asmWLCgsLWxsTF6fCwkKVlpa2e81f//pXFRQUaObMmcrJydHQoUP14IMPKhA4+xdYcXGx0tPTQz95eXmWfxZYIxIFxfTcAID72BZujhw5okAgoJycnDbHc3JyVFZW1u41e/bs0fPPP69AIKA1a9Zo/vz5WrRokX7961+f9X3mzp2rEydOhH4OHDhg6eeAdUI1NxZOBa+j5wYAXCfsb5Hdu3drxYoV2r17tx599FFlZ2frb3/7m/r3769LL73UyjaGBINBZWdn68knn5TX69Xo0aN18OBBPfLII1q4cGG71/j9fvn9/oi0B9aydPuF0ArF9NwAgNuE1XPz2muvadiwYdq0aZNWr16tmpoaSdI777xz1pDxeVlZWfJ6vSovL29zvLy8XLm5ue1e07t3b1188cXyeluHLS655BKVlZWpoaEhnI8CB6ltKShO9lFQDAAIX1jhZs6cOfr1r3+tdevWyefzhY5/4xvf0D/+8Y8O3cPn82n06NEqKSkJHQsGgyopKVFBQUG711x++eX6+OOPFQwGQ8c+/PBD9e7du0070D2Za9KkWNBzk3zasFQwaHT5fgCA7iOscPPee+/puuuuO+N4dna2jhw50uH7FBUV6amnntLTTz+tDz74QD/96U9VW1sbmj01depUzZ07N3T+T3/6U1VWVmr27Nn68MMP9fLLL+vBBx/UzJkzw/kYcBhzZlOyBeHm9IBU10jdDQC4SVjfIhkZGTp06JAGDBjQ5vi2bdvUt2/fDt9n0qRJOnz4sBYsWKCysjKNHDlSa9euDRUZ79+/X3FxrfkrLy9Pr7zyiu644w4NHz5cffv21ezZs3XXXXeF8zHgII2BoBqamnvkelgwLJWYEKc4jxQ0mjfPtKI3CADQPYT1X/zvfe97uuuuu/Tcc8/J4/EoGAzq//7v/3TnnXdq6tSpnbrXrFmzNGvWrHZf27BhwxnHCgoKOjz0he6j7rQ9oKxY58bj8aiHL17V9U2qqW/SmUtCAgBiVVjDUg8++KAGDx6svLw81dTUaMiQIfqXf/kXXXbZZZo3b57VbYQLmNPAfd44+eKtWaGgda0bhqUAwE3C+ieyz+fTU089pfnz52vHjh2qqanRqFGjdNFFF1ndPriEOasp2YIF/EytqxQzYwoA3KRL/f/9+/dX//79rWoLXMwsJrZiAT8T08EBwJ06/E1SVFTU4ZsuXrw4rMbAvepCC/hZ13NjBiUW8gMAd+lwuNm2bVub37du3aqmpiYNGjRIUvN6M+aqwUBnWblppskMSmzBAADu0uFvkvXr14f+vHjxYqWmpurpp59WZmamJOnYsWOaPn26rrjiCutbiZhn5aaZJjbPBAB3CmtayqJFi1RcXBwKNpKUmZmpX//611q0aJFljYN7WLlppsnsBWK2FAC4S1jhpqqqSocPHz7j+OHDh1VdXd3lRsF9rNw005TCbCkAcKWwws11112n6dOna/Xq1fr000/16aef6r/+679088036/rrr7e6jXABKzfNNJlBiYJiAHCXsP6ZvGzZMt155536/ve/r8bGxuYbxcfr5ptv1iOPPGJpA+EOVm6aaTKHuOoINwDgKmF9kyQnJ+vxxx/XI488ot27d0uSLrjgAvXo0cPSxsE9QptmWjpbyuy5oeYGANykS98kPXr00PDhw61qC1ysNhLr3ISmgtNzAwBuEla4ufLKK+XxeM76+quvvhp2g+BOrTU31g9LMRUcANwlrG+SkSNHtvm9sbFR27dv144dOzRt2jQr2gWXMXtXIrHODQXFAOAuYYWb3/72t+0ev+eee1RTU9OlBsGdQsNSrFAMAOiisKaCn80PfvADLV++3MpbwiVCBcX03AAAusjScFNaWqrExEQrbwmXMKdrR2IqeG19kwzDsOy+AABnC+ub5PML9RmGoUOHDuntt9/W/PnzLWkY3CWSG2cGDam+KajEBOt6hQAAzhXWN0laWlqb2VJxcXEaNGiQ7rvvPl199dWWNQ7uYBhGRDbOPD0o1dQ3EW4AwCXCCjcrV660uBlws4ZAUE3B5mEjK/eW8sZ5lJTg1cnGgOrqA1KKZbcGADhYWDU3AwcO1NGjR884fvz4cQ0cOLDLjYK71J22gnCyxb0rFBUDgPuEFW727dunQODM6bX19fU6ePBglxsFdzGDhz8+TvFeS2vc2RkcAFyoU2MAf/3rX0N/fuWVV5Senh76PRAIqKSkRPn5+ZY1Du7QWm9j3ZCUKZlVigHAdTr1bTJx4kRJksfjOWMl4oSEBOXn52vRokWWNQ7uUBuB1YlN5tTyWjbPBADX6FS4CQaDkqQBAwborbfeUlZWVkQaBXeJxOrEpmSGpQDAdcL6Ntm7d6/V7YCLtW6aaX3PTQ8/w1IA4DYdDjePPfaYfvKTnygxMVGPPfbYOc+97bbbutwwuEfrppnW99ykUHMDAK7T4W+T3/72t5oyZYoSExPPunGm1FyPQ7hBZ0RnWIqaGwBwiw5/m5w+FMWwFKwUiU0zTSkMSwGA61i7qAgQhkhsmmlqnQpOzw0AuEWHv02Kioo6fNPFixeH1Ri4U6jnJgLDUqFF/Oi5AQDX6PC3ybZt2zp03ukbagId0VpzE8HZUkwFBwDX6HC4Wb9+fSTbARdrrblhhWIAQNd1uebmwIEDOnDggBVtgUu11tywQjEAoOvCCjdNTU2aP3++0tPTlZ+fr/z8fKWnp2vevHlqbGy0uo2IcebGmZGouTG3dGBXcABwj7C+TW699VatXr1aDz/8sAoKCiRJpaWluueee3T06FE98cQTljYSsa1148zI1dzUUXMDAK4RVrh55plntGrVKl1zzTWhY8OHD1deXp4mT55MuEGnhDbOjEjPDcNSAOA2YQ1L+f1+5efnn3F8wIAB8vl8XW0TXCY0WyoCBcXmDKyGQFANTUHL7w8AcJ6wws2sWbN0//33q76+PnSsvr5eDzzwgGbNmmVZ4+AOdVHYOFNiaAoA3CKsfypv27ZNJSUl6tevn0aMGCFJeuedd9TQ0KCrrrpK119/fejc1atXW9NSxCTDMELDUpFYoTjBGydffJwamoKqqW9SRjI9iwAQ68L6NsnIyNANN9zQ5lheXp4lDYK71DcFFTSa/xyJdW6k5qGphqZgqHAZABDbwvo2WbFihdXtgEudPkU7KcH6YSmpeWjqWF0j08EBwCXYOBO2MuttkhK88sZFZusOcxZWHTOmAMAVwuq5OXr0qBYsWKD169eroqJCwWDbWSiVlZWWNA6xLzQNPEJDUs33ZiE/AHCTsL5RbrrpJn388ce6+eablZOTw2aZCFvrNPDIDEk135v9pQDATcIKN2+88YY2btwYmikFhKu6PnIzpUyhYSmmggOAK4RVczN48GCdPHnS6rbAhWqjEW5a7l1DzQ0AuEJY4ebxxx/X3Xffrddee01Hjx5VVVVVmx+go2pORSPcNA95MSwFAO4Q9jo3VVVV+sY3vtHmuGEY8ng8CgT4FzI6xizyTUmMfM9NLcNSAOAKYX2jTJkyRQkJCXrmmWcoKEaX1ERwXylTCgXFAOAqYX2j7NixQ9u2bdOgQYOsbg9cxhyWSo1guDH3rKplhWIAcIWwam7GjBmjAwcOWN0WuFAk95UyMRUcANwlrG+UW2+9VbNnz9YvfvELDRs2TAkJCW1eHz58uCWNQ+yrPsWwFADAWmF9o0yaNEmS9KMf/eiM1ygoRmfURqGgODQsxVRwAHCFsL5R9u7da3U74FI1UVjnJoXZUgDgKmF9o5x//vmSpPfff1/79+9XQ0ND6DWPxxN6Hfgi1VFY5ybZZw5L0XMDAG4QVkHxnj17NGLECA0dOlTXXnutJk6cqIkTJ+q6667TxIkTO32/pUuXKj8/X4mJiRo3bpw2b97coetWrVolj8cT1nvCGUIFxREclqLmBgDcJaxwM3v2bA0YMEAVFRVKTk7Wjh079Prrr2vMmDHasGFDp+717LPPqqioSAsXLtTWrVs1YsQIjR8/XhUVFee8bt++fbrzzjt1xRVXhPMR4BDRXKH4ZGNAgaARsfcBADhDWOGmtLRU9913n7KyshQXFyev16uvfvWrKi4u1m233dapey1evFgzZszQ9OnTNWTIEC1btkzJyclavnz5Wa8JBAKaMmWK7r33Xg0cODCcjwCHMIeKojEVXGLzTABwg7DCTSAQUGpqqiQpKytLn332maTmWpxdu3Z1+D4NDQ3asmWLCgsLWxsUF6fCwkKVlpae9br77rtP2dnZuvnmm7/wPerr69n7yqHqmwJqCAQlRXYquD8+Tt645lW0qbsBgNgX1jfK0KFD9c4772jAgAEaN26cHn74Yfl8Pj355JOd6kk5cuSIAoGAcnJy2hzPycnRzp07271m48aN+uMf/6jt27d36D2Ki4t17733drhNiJ7Tg0Yke248Ho+SfV5Vn2pixhQAuEBYPTfz5s1TMNj8L+777rtPe/fu1RVXXKE1a9boscces7SBp6uurtZNN92kp556SllZWR26Zu7cuTpx4kToh5WVncOst0n2eUM9K5FCUTEAuEdY/1weP3586M8XXnihdu7cqcrKSmVmZnZqE82srCx5vV6Vl5e3OV5eXq7c3Nwzzt+9e7f27dunCRMmhI6ZISs+Pl67du3SBRdc0OYav98vv9/f4TYheqrrGyVFdkjKZL5HDeEGAGJeWD037enZs2endwf3+XwaPXq0SkpKQseCwaBKSkpUUFBwxvmDBw/We++9p+3bt4d+vvOd7+jKK6/U9u3blZeX1+XPgegxh6UiuWmmqUfLKsV11NwAQMyL/LfKFygqKtK0adM0ZswYjR07VkuWLFFtba2mT58uSZo6dar69u2r4uJiJSYmaujQoW2uz8jIkKQzjsP5amzouaHmBgBin+3hZtKkSTp8+LAWLFigsrIyjRw5UmvXrg0VGe/fv19xcZZ1MMFBaqIwDdxkrlLMsBQAxD7bw40kzZo1S7NmzWr3tS9aFHDlypXWNwhREVrAL4KrE5tS/AxLAYBb0CUC29RGYdNMEwXFAOAehBvYptqGcMMKxQAQ+wg3sI05LBWVguJQzQ3DUgAQ6wg3sI05LJUahZobc/NMFvEDgNhHuIFtahiWAgBEAOEGtjHDDSsUAwCsRLiBbaLac2OuUNxAzQ0AxDrCDWzDVHAAQCQQbmCb6qgu4seu4ADgFoQb2Caaw1LJbJwJAK5BuIEtDMOI6rBUymkbZxqGEfH3AwDYh3ADW5xqDKop2BwyojEsZdbcBA2KigEg1hFuYIvqU42SpDhP60ymSEr2eeWN87S8N3U3ABDLCDewRdUpc3XiBHk8noi/n8fjUVpLD1FVS7ACAMQmwg1sYQaMaGy9YEpLSmh+75OEGwCIZYQb2MIcGkpLTIjae5pBimEpAIhthBvYwuw9iWrPTUuQYlgKAGIb4Qa2qD6t5iZaQuGGYSkAiGmEG9jC7D1JS4pez01qqKCYYSkAiGWEG9jCnAoezZqbUEExw1IAENMIN7BF67CUDTU3J+m5AYBYRriBLcy6l+j23LDODQC4AeEGtrCj5yaVgmIAcAXCDWzRWlAczdlSrHMDAG5AuIEtbKm5oaAYAFyBcANb2LvODT03ABDLCDewRWtBsR3r3NBzAwCxjHCDqAsGDdU02NBz0zIs1dAU1KnGQNTeFwAQXYQbRF11fZMMo/nPUZ0t5Y+Xx9PSBoqKASBmEW4QdebqxL74OCUmeKP2vnFxHqX4GJoCgFhHuEHUmQW90ay3MZlDU/TcAEDsItwg6syem2jW25hCRcUs5AcAMYtwg6gze03s7LlhWAoAYhfhBlFXZWPPTVqo54ZhKQCIVYQbRN2JliGh9ChuvWAyF/KrpucGAGIW4QZRd7yuJdwk2xBuGJYCgJhHuEHUmT03GTb03KQyLAUAMY9wg6hzwrAUPTcAELsIN4i643UNkqQMW4almntuWOcGAGIX4QZR19pz44v6e7fuDE7PDQDEKsINou64WXNjQ8+NOf38BOEGAGIW4QZRd6LOvpobM1AdJ9wAQMwi3CCqDMNonS1lQ8+N+Z4n6hplmFuTAwBiCuEGUVXbEFBTsDlUZNhQc5OZ3PyeDYGg6hoCUX9/AEDkEW4QVeZMKZ83TokJ0f/rl+zzKsHraW4LQ1MAEJMIN4iq0Eyp5AR5PJ6ov7/H41FGS+/NsdqGqL8/ACDyCDeIKrOY2I7ViU2ZycyYAoBYRrhBVB23cXVik1nrc6yOnhsAiEWEG0SVnTOlTOZ7H6uj5wYAYhHhBlEV2hHchplSptbp4PTcAEAsItwgquzcNNNkTgen5wYAYhPhBlF14qR9m2aazNlSxwk3ABCTCDeIKjNQOKHm5jjDUgAQkwg3iCpzhpK9w1JmQTHhBgBiEeEGUXWstrnnpmcP+wqKzWJmVigGgNhEuEFUHW1ZFdjOcJPZwxyWItwAQCxyRLhZunSp8vPzlZiYqHHjxmnz5s1nPfepp57SFVdcoczMTGVmZqqwsPCc58M5gkEjNBTUq4fftnZkhgqKGxQMsjM4AMQa28PNs88+q6KiIi1cuFBbt27ViBEjNH78eFVUVLR7/oYNGzR58mStX79epaWlysvL09VXX62DBw9GueXorKpTjQq0hAmz98QOZr1P0JCq65tsawcAIDJsDzeLFy/WjBkzNH36dA0ZMkTLli1TcnKyli9f3u75f/7zn/Wzn/1MI0eO1ODBg/WHP/xBwWBQJSUlUW45Ossckkrxx8sf77WtHYkJXiUlNL8/M6YAIPbYGm4aGhq0ZcsWFRYWho7FxcWpsLBQpaWlHbpHXV2dGhsb1bNnz3Zfr6+vV1VVVZsf2OOYA+ptTGYbKtkZHABijq3h5siRIwoEAsrJyWlzPCcnR2VlZR26x1133aU+ffq0CUinKy4uVnp6eugnLy+vy+1GeJxQTGzqlUK4AYBYZfuwVFc89NBDWrVqlV544QUlJia2e87cuXN14sSJ0M+BAwei3EqYzCDRywnhpqUNR2sINwAQa+LtfPOsrCx5vV6Vl5e3OV5eXq7c3NxzXvtv//Zveuihh/T3v/9dw4cPP+t5fr9ffr99M3PQygw3mU4INynNfyeO1Nbb3BIAgNVs7bnx+XwaPXp0m2Jgszi4oKDgrNc9/PDDuv/++7V27VqNGTMmGk2FBcxeEkf03LQMSx2ppucGAGKNrT03klRUVKRp06ZpzJgxGjt2rJYsWaLa2lpNnz5dkjR16lT17dtXxcXFkqTf/OY3WrBggZ555hnl5+eHanNSUlKUkpJi2+fAF6ts6SVxQs1NVss6O0fpuQGAmGN7uJk0aZIOHz6sBQsWqKysTCNHjtTatWtDRcb79+9XXFxrB9MTTzyhhoYG3XjjjW3us3DhQt1zzz3RbDo6qbLO/q0XTGbPDTU3ABB7bA83kjRr1izNmjWr3dc2bNjQ5vd9+/ZFvkGICLPnxgwWdsoya25q6LkBgFjTrWdLoXupbOklMbc/sFOo54ap4AAQcwg3iArDMEJBws59pUxmz01lLftLAUCsIdwgKuoaAqpvCkqSejpgWMrsPQoEDZ04ye7gABBLCDeICrO2JTEhTj189u0rZfLFx4U20GTGFADEFsINoqKiujlAZKcmyuPx2NyaZqG1bpgxBQAxhXCDqKioMsON/fU2ptBaN4QbAIgphBtERXnVKUlSdpqDwk2qOWOKYSkAiCWEG0TF6cNSTmHO2jpSTbgBgFhCuEFUVFQ7r+fmvJYhsgrCDQDEFMINouKwA3tuctOa21LWMmQGAIgNhBtEhRMLinPSW8LNCcINAMQSwg2iwonDUvTcAEBsItwg4uqbAjrWsiN4jgOHpY7XNepUY8Dm1gAArEK4QcSZ9TY+b5wykhNsbk2rtKR4JSY0/1+gnN4bAIgZhBtEnDkb6bxUv2NWJ5Ykj8fTOjRF3Q0AxAzCDSLOLCY+z0HFxKYc6m4AIOYQbhBxh81iYgeGm9yWGVMMSwFA7CDcIOI+PX5SktQnI8nmlpypdViKhfwAIFYQbhBxB481h5t+mc4LN+awFD03ABA7CDeIuIMtPTd9ndhzk07NDQDEGsINIu7Tlp6bvg7uuWG2FADEDsINIupUYyC0zk2/zGSbW3OmvJbAdejESTUGgja3BgBgBcINIupQS49IUoJXmQ5awM90Xqpf/vg4BQ3ps5bhMwBA90a4QUQdPG1IykkL+Jk8Ho/692zuUTpQSbgBgFhAuEFEfXqsTpIzi4lNeS3hZn9lnc0tAQBYgXCDiDJnSjlxGripP+EGAGIK4QYRddDBM6VMeaFhKcINAMQCwg0iKjQN3MHDUqGam2OEGwCIBYQbRNTeo7WSpPxePWxuydnl9WwOXgxLAUBsINwgYqpPNYbWuMnPcnC4aVl/53hdo06cbLS5NQCAriLcIGL2HWnuCclK8Sk9yXlr3Jh6+OOVleKTRN0NAMQCwg0iZs+RGknSAAf32pjMYbPdh2tsbgkAoKsIN4iYvUea6226Q7i5ODdVkvRhebXNLQEAdBXhBhGz57AZblJsbskXuzi7uY0fltNzAwDdHeEGEWP2glyU3Q3CTQ49NwAQKwg3iIjGQDBUvzK4d6rNrfli5rDU/so6nWwI2NwaAEBXEG4QEXsO16oxYCjVH+/oBfxMWSl+9ezhk2FQVAwA3R3hBhGxs6xKkjQoN9WRu4G3xxw+21XG0BQAdGeEG0TEzpaAMCjX+UNSJupuACA2EG4QEe99ekKSNKRPms0t6TizrTs+O2FzSwAAXUG4geWCQUPvfHpckjQyL8PWtnTGiH4ZkqR3D5xQMGjY2xgAQNgIN7DcniO1qj7VpMSEOA3K6U7DUilKSvCqur4ptLoyAKD7IdzAcu8cOC5JGtY3XfHe7vNXLN4bp2F90yVJ2/Yft7cxAICwdZ9vHnQbb+2rlNS9hqRMo/pnSJLe3nfM3oYAAMJGuIHlSvcclSRddkGWzS3pvK9c0EtS62cAAHQ/hBtY6tNjdfrkaJ28cR59eUBPu5vTaV/O76n4OI/2V9bpQGWd3c0BAISBcANL/d/HRyRJw/ulK8Ufb3NrOi/FH68RLcNpG1s+CwCgeyHcwFL/+89ySdKVg7Jtbkn4vn7xeZKkde+X29wSAEA4CDewTG19k95o6e24+tIcm1sTvvFDcyU199zU1DfZ3BoAQGcRbmCZv39QroamoM7vldyt1rf5vIuyU5TfK1kNTUH9nd4bAOh2CDewzLNvHZAkTRzZt9tsltkej8ej60b1k9T6mQAA3QfhBpbYe6RWb+4+Ko9H+u6YfnY3p8tuHNNPHk/zlPDdh1mtGAC6E8INLPH4+o8lSd8YlK1+mck2t6br+mYk6arBzXVDT2zYbXNrAACdQbhBl31UXq0Xth2UJM36xoU2t8Y65md5YdtBfVhebXNrAAAdRbhBlwSChu5+cYeagoYKL8nRqP6ZdjfJMiPzMnT1kJzmz/jCewqwUzgAdAuEG3TJI6/s0ua9lUpK8GrhhCF2N8dyC79zqZJ9Xr2175geXrvT7uYAADrAEeFm6dKlys/PV2JiosaNG6fNmzef8/znnntOgwcPVmJiooYNG6Y1a9ZEqaUwNQWCeuhvO7XsteZ6lIduGKa8nt2/1ubz+mYk6Tc3DJck/f71PXpwzQdqDARtbhUA4FxsDzfPPvusioqKtHDhQm3dulUjRozQ+PHjVVFR0e75b775piZPnqybb75Z27Zt08SJEzVx4kTt2LEjyi13p0DQ0PpdFbr+iTdDweaX3xqkfx3Z1+aWRc6EEX0095rBkqQnX9+j6x9/U6/uLGeYCgAcymMYhq3/hR43bpy+/OUv69///d8lScFgUHl5ebr11ls1Z86cM86fNGmSamtr9dJLL4WOfeUrX9HIkSO1bNmyL3y/qqoqpaen68SJE0pLS7Psc9Q3BXS4ul7tPU3zmCHjc7+ffo7R5ljrfc51zVnu+7njbe937mtPf4+GpqCO1TXoWG2DyqpOaVdZtTbvrdTR2gZJUmpivB68bpgmjOhz5oeOQS+9+5nmrn5P1aeaVy3u2cOnsfk9Nbh3qnqnJ6pnD78ykhPk88Yp3utp+d84ec+x5k84ywF14yWEALiELz5O2amJlt6zM9/ftu5s2NDQoC1btmju3LmhY3FxcSosLFRpaWm715SWlqqoqKjNsfHjx+vFF19s9/z6+nrV19eHfq+qqup6w9vxz8+qdP3jb0bk3k6Tlhiv/29Mnmb8y0DlpFn7l9fJ/t/wPhqb31N/2LhXqzbvV2Vtg9b+s0xr/1lmd9MAwFG+1D9Dq392uW3vb2u4OXLkiAKBgHJy2u5DlJOTo5072y/eLCsra/f8srL2v2CKi4t17733WtPgc4jzeJSY0DzK51HrP63Nf2V7Qr972vx+2qlnnnOWa9sea3vkzGtOv/8X3/f0did445SRnKCeyT717OHTxTmpGto3XWPyM5XgtX1E0xbZaYn61bcv0S/GD9KWT45px8ET+rC8Woer61VZ26ATJxvVGDDUGAiqKWiosSmo4Oe6zdrrKm23x6+dM+3tZwWAjrH7O8LWcBMNc+fObdPTU1VVpby8PMvfZ2Rehnbef43l94UzJXjj9JWBvfSVgb3sbgoA4HNsDTdZWVnyer0qL2+7OWF5eblyc3PbvSY3N7dT5/v9fvn9fmsaDAAAHM/WfiOfz6fRo0erpKQkdCwYDKqkpEQFBQXtXlNQUNDmfElat27dWc8HAADuYvuwVFFRkaZNm6YxY8Zo7NixWrJkiWprazV9+nRJ0tSpU9W3b18VFxdLkmbPnq2vfe1rWrRoka699lqtWrVKb7/9tp588kk7PwYAAHAI28PNpEmTdPjwYS1YsEBlZWUaOXKk1q5dGyoa3r9/v+LiWjuYLrvsMj3zzDOaN2+efvWrX+miiy7Siy++qKFDh9r1EQAAgIPYvs5NtEVqnRsAABA5nfn+dud8XgAAELMINwAAIKYQbgAAQEwh3AAAgJhCuAEAADGFcAMAAGIK4QYAAMQUwg0AAIgphBsAABBTbN9+IdrMBZmrqqpsbgkAAOgo83u7IxsruC7cVFdXS5Ly8vJsbgkAAOis6upqpaenn/Mc1+0tFQwG9dlnnyk1NVUej8fSe1dVVSkvL08HDhxg36oI4jlHB885OnjO0cOzjo5IPWfDMFRdXa0+ffq02VC7Pa7ruYmLi1O/fv0i+h5paWn8HycKeM7RwXOODp5z9PCsoyMSz/mLemxMFBQDAICYQrgBAAAxhXBjIb/fr4ULF8rv99vdlJjGc44OnnN08Jyjh2cdHU54zq4rKAYAALGNnhsAABBTCDcAACCmEG4AAEBMIdwAAICYQrixyNKlS5Wfn6/ExESNGzdOmzdvtrtJ3VpxcbG+/OUvKzU1VdnZ2Zo4caJ27drV5pxTp05p5syZ6tWrl1JSUnTDDTeovLzcphbHhoceekgej0e333576BjP2ToHDx7UD37wA/Xq1UtJSUkaNmyY3n777dDrhmFowYIF6t27t5KSklRYWKiPPvrIxhZ3P4FAQPPnz9eAAQOUlJSkCy64QPfff3+b/Yh4zp33+uuva8KECerTp488Ho9efPHFNq935JlWVlZqypQpSktLU0ZGhm6++WbV1NREpsEGumzVqlWGz+czli9fbvzzn/80ZsyYYWRkZBjl5eV2N63bGj9+vLFixQpjx44dxvbt241vf/vbRv/+/Y2amprQObfccouRl5dnlJSUGG+//bbxla98xbjssstsbHX3tnnzZiM/P98YPny4MXv27NBxnrM1KisrjfPPP9/44Q9/aGzatMnYs2eP8corrxgff/xx6JyHHnrISE9PN1588UXjnXfeMb7zne8YAwYMME6ePGljy7uXBx54wOjVq5fx0ksvGXv37jWee+45IyUlxXj00UdD5/CcO2/NmjXG3XffbaxevdqQZLzwwgttXu/IM/3Wt75ljBgxwvjHP/5hvPHGG8aFF15oTJ48OSLtJdxYYOzYscbMmTNDvwcCAaNPnz5GcXGxja2KLRUVFYYk47XXXjMMwzCOHz9uJCQkGM8991zonA8++MCQZJSWltrVzG6rurrauOiii4x169YZX/va10Lhhudsnbvuusv46le/etbXg8GgkZubazzyyCOhY8ePHzf8fr/xn//5n9FoYky49tprjR/96Edtjl1//fXGlClTDMPgOVvh8+GmI8/0/fffNyQZb731Vuicv/3tb4bH4zEOHjxoeRsZluqihoYGbdmyRYWFhaFjcXFxKiwsVGlpqY0tiy0nTpyQJPXs2VOStGXLFjU2NrZ57oMHD1b//v157mGYOXOmrr322jbPU+I5W+mvf/2rxowZo+9+97vKzs7WqFGj9NRTT4Ve37t3r8rKyto86/T0dI0bN45n3QmXXXaZSkpK9OGHH0qS3nnnHW3cuFHXXHONJJ5zJHTkmZaWliojI0NjxowJnVNYWKi4uDht2rTJ8ja5buNMqx05ckSBQEA5OTltjufk5Gjnzp02tSq2BINB3X777br88ss1dOhQSVJZWZl8Pp8yMjLanJuTk6OysjIbWtl9rVq1Slu3btVbb711xms8Z+vs2bNHTzzxhIqKivSrX/1Kb731lm677Tb5fD5NmzYt9Dzb+28Jz7rj5syZo6qqKg0ePFher1eBQEAPPPCApkyZIkk85wjoyDMtKytTdnZ2m9fj4+PVs2fPiDx3wg0cb+bMmdqxY4c2btxod1NizoEDBzR79mytW7dOiYmJdjcnpgWDQY0ZM0YPPvigJGnUqFHasWOHli1bpmnTptncutjxl7/8RX/+85/1zDPP6NJLL9X27dt1++23q0+fPjxnF2FYqouysrLk9XrPmD1SXl6u3Nxcm1oVO2bNmqWXXnpJ69evV79+/ULHc3Nz1dDQoOPHj7c5n+feOVu2bFFFRYW+9KUvKT4+XvHx8Xrttdf02GOPKT4+Xjk5OTxni/Tu3VtDhgxpc+ySSy7R/v37JSn0PPlvSdf84he/0Jw5c/S9731Pw4YN00033aQ77rhDxcXFknjOkdCRZ5qbm6uKioo2rzc1NamysjIiz51w00U+n0+jR49WSUlJ6FgwGFRJSYkKCgpsbFn3ZhiGZs2apRdeeEGvvvqqBgwY0Ob10aNHKyEhoc1z37Vrl/bv389z74SrrrpK7733nrZv3x76GTNmjKZMmRL6M8/ZGpdffvkZyxl8+OGHOv/88yVJAwYMUG5ubptnXVVVpU2bNvGsO6Gurk5xcW2/2rxer4LBoCSecyR05JkWFBTo+PHj2rJlS+icV199VcFgUOPGjbO+UZaXKLvQqlWrDL/fb6xcudJ4//33jZ/85CdGRkaGUVZWZnfTuq2f/vSnRnp6urFhwwbj0KFDoZ+6urrQObfccovRv39/49VXXzXefvtto6CgwCgoKLCx1bHh9NlShsFztsrmzZuN+Ph444EHHjA++ugj489//rORnJxs/Md//EfonIceesjIyMgw/vu//9t49913jX/9139linInTZs2zejbt29oKvjq1auNrKws45e//GXoHJ5z51VXVxvbtm0ztm3bZkgyFi9ebGzbts345JNPDMPo2DP91re+ZYwaNcrYtGmTsXHjRuOiiy5iKrjT/e53vzP69+9v+Hw+Y+zYscY//vEPu5vUrUlq92fFihWhc06ePGn87Gc/MzIzM43k5GTjuuuuMw4dOmRfo2PE58MNz9k6//M//2MMHTrU8Pv9xuDBg40nn3yyzevBYNCYP3++kZOTY/j9fuOqq64ydu3aZVNru6eqqipj9uzZRv/+/Y3ExERj4MCBxt13323U19eHzuE5d9769evb/W/ytGnTDMPo2DM9evSoMXnyZCMlJcVIS0szpk+fblRXV0ekvR7DOG3ZRgAAgG6OmhsAABBTCDcAACCmEG4AAEBMIdwAAICYQrgBAAAxhXADAABiCuEGAADEFMINAACIKYQbAFGzYcMGeTyeMzbijJaSkhJdcsklCgQCltxvzpw5uvXWWy25FwDrsEIxgIj4+te/rpEjR2rJkiWhYw0NDaqsrFROTo48Hk/U2zR69GgVFRVpypQpltzvyJEjGjhwoLZv366BAwdack8AXUfPDYCo8fl8ys3NtSXYbNy4Ubt379YNN9xg2T2zsrI0fvx4PfHEE5bdE0DXEW4AWO6HP/yhXnvtNT366KPyeDzyeDzat2/fGcNSK1euVEZGhl566SUNGjRIycnJuvHGG1VXV6enn35a+fn5yszM1G233dZmKKm+vl533nmn+vbtqx49emjcuHHasGHDOdu0atUqffOb31RiYmLo2D333KORI0fqT3/6k/Lz85Wenq7vfe97qq6uDp3z/PPPa9iwYUpKSlKvXr1UWFio2tra0OsTJkzQqlWrrHlwACxBuAFguUcffVQFBQWaMWOGDh06pEOHDikvL6/dc+vq6vTYY49p1apVWrt2rTZs2KDrrrtOa9as0Zo1a/SnP/1Jv//97/X888+Hrpk1a5ZKS0u1atUqvfvuu/rud7+rb33rW/roo4/O2qY33nhDY8aMOeP47t279eKLL+qll17SSy+9pNdee00PPfSQJOnQoUOaPHmyfvSjH+mDDz7Qhg0bdP311+v00fyxY8fq008/1b59+8J8WgCsFm93AwDEnvT0dPl8PiUnJys3N/ec5zY2NuqJJ57QBRdcIEm68cYb9ac//Unl5eVKSUnRkCFDdOWVV2r9+vWaNGmS9u/frxUrVmj//v3q06ePJOnOO+/U2rVrtWLFCj344IPtvs8nn3wSOv90wWBQK1euVGpqqiTppptuUklJiR544AEdOnRITU1Nuv7663X++edLkoYNG9bmevOen3zyifLz8zv+kABEDOEGgK2Sk5NDwUaScnJylJ+fr5SUlDbHKioqJEnvvfeeAoGALr744jb3qa+vV69evc76PidPnmwzJGXKz88PBRtJ6t27d+i9RowYoauuukrDhg3T+PHjdfXVV+vGG29UZmZm6PykpCRJzT1QAJyBcAPAVgkJCW1+93g87R4LBoOSpJqaGnm9Xm3ZskVer7fNeacHos/LysrSsWPHOvT+5nt5vV6tW7dOb775pv73f/9Xv/vd73T33Xdr06ZNGjBggCSpsrJSknTeeed15OMCiAJqbgBEhM/ns2w9mdONGjVKgUBAFRUVuvDCC9v8nGsIbNSoUXr//fc7/X4ej0eXX3657r33Xm3btk0+n08vvPBC6PUdO3YoISFBl156aVifB4D16LkBEBH5+fnatGmT9u3bp5SUFPXs2dOS+1588cWaMmWKpk6dqkWLFmnUqFE6fPiwSkpKNHz4cF177bXtXjd+/Hg9/fTTnXqvTZs2qaSkRFdffbWys7O1adMmHT58WJdccknonDfeeENXXHFFaHgKgP3ouQEQEXfeeae8Xq+GDBmi8847T/v377fs3itWrNDUqVP185//XIMGDdLEiRP11ltvqX///me9ZsqUKfrnP/+pXbt2dfh90tLS9Prrr+vb3/62Lr74Ys2bN0+LFi3SNddcEzpn1apVmjFjRpc+DwBrsUIxANf4xS9+oaqqKv3+97+35H5/+9vf9POf/1zvvvuu4uPpCAecgp4bAK5x99136/zzzw8VDHdVbW2tVqxYQbABHIaeGwAAEFPouQEAADGFcAMAAGIK4QYAAMQUwg0AAIgphBsAABBTCDcAACCmEG4AAEBMIdwAAICYQrgBAAAx5f8H79t/llPCsMcAAAAASUVORK5CYII=", "text/plain": [ "
" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "# create a Gaussian pulse\n", "class GaussianPulse(nx.FunctionTime):\n", " def __init__(self, t0, fwhm):\n", " super().__init__(\"Gaussian pulse\")\n", " # create a self. instance of each variable\n", " self.t0 = t0\n", " self.fwhm = fwhm\n", " # register them in the self.fit_variables \n", " self.fit_variables = [t0, fwhm] #\n", "\n", " def Function(self, t):\n", " sigma = self.fwhm.value * 0.424660900144009521361\n", " if sigma == 0:\n", " sigma = 1e-299\n", " ret = np.exp(-1/2 * (t-self.t0.value)**2/ sigma**2)\n", " return ret\n", "\n", "# create Nexus var for a Gaussian pulse\n", "t0 = nx.Var(20, min=10, max=30, fit=True, id=\"t0\")\n", "fwhm = nx.Var(5, min=1, max=10, fit=True, id=\"fwhm\")\n", "\n", "func_gaussian = GaussianPulse(t0, fwhm)\n", "\n", "# print the fit_variables of the function\n", "for elem in list(func_gaussian.fit_variables):\n", " print(elem)\n", "\n", "times = np.linspace(0, 100, 1001)\n", "\n", "pulse = [func_gaussian.Function(t) for t in times]\n", "\n", "plt.plot(times, pulse)\n", "plt.xlabel(\"time (ns)\")\n", "plt.ylabel(\"amplitude\")\n", "plt.savefig(\"gaussian_pulse.png\")\n", "plt.show()" ] } ], "metadata": { "kernelspec": { "display_name": "Python 3 (ipykernel)", "language": "python", "name": "python3" }, "language_info": { "codemirror_mode": { "name": "ipython", "version": 3 }, "file_extension": ".py", "mimetype": "text/x-python", "name": "python", "nbconvert_exporter": "python", "pygments_lexer": "ipython3", "version": "3.12.0" } }, "nbformat": 4, "nbformat_minor": 5 }