{ "cells": [ { "cell_type": "markdown", "id": "6434d0d1-7a2c-424f-a591-dba64e1b8e45", "metadata": {}, "source": [ "# X-Ray Absorption in Cu-W" ] }, { "cell_type": "code", "execution_count": 7, "id": "db62beeb-4162-495e-bb36-2a82ac33a92a", "metadata": { "tags": [] }, "outputs": [], "source": [ "import roentgen\n", "import roentgen.absorption\n", "import astropy.units as u\n", "import numpy as np\n", "import matplotlib.pyplot as plt\n", "from astropy.visualization import quantity_support" ] }, { "cell_type": "code", "execution_count": 3, "id": "3e9857fe-a52b-4610-a854-711415a1e9ec", "metadata": { "tags": [] }, "outputs": [], "source": [ "alloy = roentgen.absorption.Material(\n", " {'Cu': 0.25, 'W': 0.75},\n", " thickness=350*u.micron,\n", " density=14.5*u.g/(u.cm**3),\n", ")" ] }, { "cell_type": "code", "execution_count": 9, "id": "e3aa8a2f-3147-4036-9a69-1b046ce9547e", "metadata": { "tags": [] }, "outputs": [], "source": [ "energy = np.linspace(1,400,1000) * u.keV" ] }, { "cell_type": "code", "execution_count": 10, "id": "3a0e907c-d796-4411-9048-70683778b18e", "metadata": { "tags": [] }, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAiMAAAGxCAYAAACwbLZkAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjUuMywgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/NK7nSAAAACXBIWXMAAA9hAAAPYQGoP6dpAABHPklEQVR4nO3deVzUdf4H8Nd3ZpgZrhnkPuUSrxBUVCLLLCm11uzYzcxNM7PNtX6Z227Z5Va7Ub9+nZtla6ldpumatmmW4ZEaeYDkiYAiIDocIsxww8z39wcwinLLzHeO1/PxmMDvfL8z70+D8erz/RyCKIoiiIiIiCQik7oAIiIicm4MI0RERCQphhEiIiKSFMMIERERSYphhIiIiCTFMEJERESSYhghIiIiSSmkLqA7TCYTzp49C09PTwiCIHU5RERE1A2iKMJgMCA4OBgyWcf9H3YRRs6ePYuwsDCpyyAiIqJeKCwsRGhoaIfP20UY8fT0BNDcGI1GI3E1RERE1B16vR5hYWHm3+MdsYsw0nprRqPRMIwQERHZma6GWHAAKxEREUmKYYSIiIgkxTBCREREkmIYISIiIkkxjBAREZGkGEaIiIhIUgwjREREJCmGESIiIpIUwwgRERFJimGEiIiIJNXjMPLzzz9jypQpCA4OhiAI2LBhQ5fX7NixAyNHjoRKpcKAAQOwcuXKXpRKREREjqjHYaS6uhrx8fFYsmRJt87Py8vD7bffjptuugmZmZlYsGABHn74Yfzwww89LpaIiIgcT483yps8eTImT57c7fOXLl2KyMhIvPnmmwCAIUOGYPfu3Xj77bcxceLEnr59nyox1KGhyWTewEcA0LqXj4CWYwJg3t5HaP94h9e3fm8+1vG5ggDIBAFKBe+cERGRc7H4rr1paWlITk5uc2zixIlYsGBBh9fU19ejvr7e/Ge9Xm+R2v70eToOFlRY5LV7y1OlgJ+nCtH+HogP1WLCkAAMDvTscsdDIiIie2XxMKLT6RAQENDmWEBAAPR6PWpra+Hq6nrFNSkpKXjppZcsXRqUchlULT0RovkfgNjyjSiaD0EUxUu+t1xNhvomGOqbcKqsGluPFeP/fszGsBAtFiTH4ObB/gwlRETkcCweRnpj0aJFWLhwofnPer0eYWFhff4+a/6UdNWvIYpdB5dLww0uOQ40P9f6fZNRxPnqeuj0dTh2Vo+9eeXYmV2Kw0WVmPPpAdweF4RX7xwGrZvLVddNRERkKyweRgIDA1FcXNzmWHFxMTQaTbu9IgCgUqmgUqksXVqfMI8BadNh0fveC62bC6L8PHBdtC8eviEK5dUN+GjnSXyyOw+bDp3DCZ0Bn88ZgyBt+//uiIiI7I3FR0smJSUhNTW1zbGtW7ciKenqeyWcgbe7EotuG4L1f74OgRo1ckuqcO9HaSg11Hd9sYV8faAQW47ozL1CREREV6PHYaSqqgqZmZnIzMwE0Dx1NzMzEwUFBQCab7HMnDnTfP6jjz6KU6dO4W9/+xuysrLwwQcf4Ouvv8aTTz7ZNy1wEnGhXlg3Lwn9vd1QWF6LuZ8dQF2j0ep1FOvr8Ld1h/DoF+l47KuDqKxptHoNRETkWHocRg4cOIARI0ZgxIgRAICFCxdixIgRePHFFwEA586dMwcTAIiMjMSmTZuwdetWxMfH480338THH38s+bReexTazw0rZ4+G1tUFmYUVeH1LltVrqK5vMn+/6dA53PnBHpwqrbJ6HURE5DgE0Q762vV6PbRaLSorK6HRaKQuR3LbT5Rg9or9AIAvH07E2AG+Vnvvk6VVmPDmTgBAiJcriipqoVErsGzmKCRG+VitDiIisn3d/f3NFbbs0E2D/DEjsT8A4PkNR9DQZLLae7dGV62rCzY+NhYJ4f2gr2vCgyv245eTZVarg4iIHAfDiJ1adNsQ+HqokFdWjc/STlvxnZvTiCAAvh4qfPlwIsYP8kNtoxEPrdyP3TkMJERE1DMMI3bKQ6XA3yYOAgC8m5pjtYGkrT0jrZOX1S5yfPRAAm4e7I+6RhMe+fwAjhRVWqUWIiJyDAwjduz3CaEYHOgJQ10TVvySZ5X3bB1gdOlKsCqFHB/+cSTGDvBBTYMRcz7dj7MVtVaph4iI7B/DiB2TyQQ8dvMAAMDy3Xkw1Fm+d+TynpFWKoUcH8xIwMAADxTr6/HQyv2oaWi64noiIqLLMYzYucmxQYj2c4e+rgmr9xVa/P3ES8aMXE7r6oIVs8fAz1OFLJ0Bz284woXRiIioSwwjdk4uEzD3higAwBd782EyWfaX/8Vs0f6S9yFervjX9BGQCcD6jCKsPXDGovUQEZH9YxhxAFOHh8BTrUD++RrszCm16HuZb9N0sv3OtVE++MutzYNrX9h4BCd0BovWRERE9o1hxAG4KuX4Q0LzrsZfpOVb9L3Mt2m6OG/ejdEYP8gP9U0m/GVtJhqN1lsLhYiI7AvDiIN4ICkcALDtRAkKy2ss9j7d6RkBmgfX/u89cdC6uuBIkR4fbD9psZqIiMi+MYw4iEhfd4wd4ANRBDZmFln8/YQu+0YAf40aL0+9BgDwr205XH+EiIjaxTDiQO4cHgIA+OZgkcVmsXS3Z6TVHfHBmBwbiCaTiKfW/sbbNUREdAWGEQcyKTYQKoUMJ0urcfSs3iLv0d0xI60EQcA/7oyFt7sSWToDlu+2zuJsRERkPxhGHIin2gW3DA0A0Nw7YgkXe0a6G0cAHw8VFk0eDAB456ccrs5KRERtMIw4mNZbNd/+dhZGC6w50ttXvGdkKEZH9ENtoxEv//dYn9ZERET2jWHEwYwb6AetqwtKDfXIKLjQ56/fOhalBx0jAJpn17xyZyzkMgFbjuqwPaukz2sjIiL7xDDiYJQKGSYM9gcA/HBE1+evf3GjvJ5fOzhQg4fGRgAAXv7uGAezEhERAIYRh3TrNYEAgC1HdX0+q+biRnm9SCMA/mdCDHzclcgrq8aqvQV9WBkREdkrhhEHdONAP6hdZDhzoRbHzvX1rJre3aZp5al2wYJbBgIA3k3Ngd4KOw0TEZFtYxhxQK5KOcbF+AEAfjha3KevfbFnpPfuGx2GKD93lFc3YOkOrsxKROTsGEYc1MSWWzU/Hu3bcSMXx4z0Po64yGVYNHkIAOCT3Xko4lRfIiKnxjDioG4e7A9BALJ0Bugq6/rsdfuiZwQAkof4IzHSG/VNJvwrNeeq6yIiIvvFMOKg+rkrMTzMCwCw40TfTaMV+yiNCIKAv00aBABYl34GBectt7kfERHZNoYRBzZ+YPMU3+19GUZavl5tzwgAJIR7Y9xAPzSZRPxrG3tHiIicFcOIA7tpcPMg1j2559HQ1DdrevRmOfjOPJkcAwBYf7AIp8uq++Q1iYjIvjCMOLDYYC18PZSoqm/CgfzyPnnNnm6U15UR/fvhpkF+MJpEvMfeESIip8Qw4sBkMgHjBjb3juw8Udo3L2ruGemblwOABcnN645sOFiEk6VVfffCRERkFxhGHNz4QX07buTimJG+SyPxYV5IHuIPkwgs2ZbbZ69LRET2gWHEwY2L8YVMALKLq3C2D9bzEC3QMwI0LxMPABt/O4szFzizhojImTCMODgvNyWGhWgBAGknz1/164no271uWsWFemHsAB8YTSI+3pVnkfcgIiLbxDDiBK4b4AsA2HOy7Kpfq69n01xq3o0DAACr9xfgfFV9n78+ERHZJoYRJ3BdtA+A5p6Rq93Fty/XGbnc2AE+iAvVoq7RhE9/OW2BdyAiIlvEMOIERoV7QymX4VxlHU5f5UqnrWHGAh0jEAQB826MBgB8mpaPqvqmvn8TIiKyOQwjTsBVKceI/l4AgD25V3er5uJGeVdXU0duvSYQUb7uqKxtxFd7CyzzJkREZFMYRpzEddHN40auehCreWsay6QRuUzAn26MAgAs35OHRmPfrBxLRES2i2HESYwd0Dxu5JeTZTCZej9uxLwCq4V6RgBg6vAQ+Hooca6yDt8f0VnujYiIyCYwjDiJuFAvuCnluFDTiCydodev00eb9nZK7SLHH68NBwAs381pvkREjo5hxEkoFTKMifQG0Nw70luipQeNtJiRGA6lXIbMwgqk51+w6HsREZG0GEacSOsU31+uYtyIJaf2XsrPU4Wpw4MBsHeEiMjRMYw4kaSo5kGs+0+Xw9jLcSOWnNp7uYeujwQAfH/kHJeIJyJyYAwjTmRIkCfclXIY6ppwopfjRqzVMwIAQ4I0uC7aByYR+Cwt3wrvSEREUmAYcSIKuQwjw/sBAPbl9e5WjSWXg2/PnJbeka/2FaCai6ARETkkhhEnMyaieRDr/tO9HRTacpumj+rpyk2D/BHp6w5DXRP+k3HGSu9KRETWxDDiZEa3zKjZd7q8V/vUXOwZ6cuqOiaTCZiV1DzN9/O0/KveW4eIiGwPw4iTGR7mBRe5gFJDPfJ7sU/NxTEj1uobAe5OCIWbUo6ckir8eqrcau9LRETWwTDiZNQucsSFegFo7h3pKdGaI1hbaNQuuHNECADgi185kJWIyNEwjDih0a3jRvJ6EUasPGak1QMtK7L+cFSHYn2dld+diIgsiWHECY2JbJ5Rs/8qekasNWak1ZAgDUZH9EOTScRX+7ibLxGRI2EYcUIJ4d4QBOD0+RqUGHrWyyDFmJFWrfvVrNpbwN18iYgcCMOIE9K6umBQgCcAYH9ez6b4WnMF1stNjg2Cr4cKJYZ6bD1WbP0CiIjIIhhGnFTrpnm9uVUDSBNGlAoZpo8JAwB8lnba+gUQEZFFMIw4qdZBrPt6OIjVPGZEgts0ADB9TH/IBODXU+XILu7dkvZERGRbGEacVGvPSJZOD0NdY7evM8+mkSaLINjLFbcMDQDAab5ERI6CYcRJBWjUCPFyhUkEfius7PZ1trAA6gPXRgAA1mcUoaaB+9UQEdk7hhEnltCyaV56fvcHsVp7o7z2XBftgwgfN1TVN+G7385JVgcREfUNhhEn1hpGMgp6EEZavkoXRZr3q7lvTH8AwCquOUJEZPcYRpzYyP4Xw4jJ1L37L1JO7b3U7xNC4SIXkFlYgWNn9dIWQ0REV4VhxIkNDvKEq4schrom5JZWdesaW+gZAQBfDxVuHRoIAFi9n70jRET2jGHEibnIZYgP0wIAMro7bsQGxoy0mt5yq+abjCLUNhglroaIiHqLYcTJtd6q6e4gVqk2ymvPddE+6O/tBkN9E747dFbqcoiIqJcYRpyceUZNNwexSrVRXnuaB7I2r8jKzfOIiOwXw4iTG9HSM3KqtBoXqhu6PP/iMFcbSCNoHsiqkAnIKKhAlo4DWYmI7FGvwsiSJUsQEREBtVqNxMRE7Nu3r9Pz33nnHQwaNAiurq4ICwvDk08+ibq6nu0WS5bh7a5ElK87AOBgYde9I7bUMwIA/p5q84qsq/cVSlwNERH1Ro/DyJo1a7Bw4UIsXrwYGRkZiI+Px8SJE1FSUtLu+atWrcIzzzyDxYsX4/jx4/jkk0+wZs0aPPvss1ddPPWNkT1Y/MyWxoy0ah3Iuj7jDAeyEhHZoR6Hkbfeegtz587F7NmzMXToUCxduhRubm5Yvnx5u+f/8ssvGDt2LO6//35ERETg1ltvxfTp07vsTSHrMS9+ll/R5bm21jMCANcP8EVoP1fo65qw+TBXZCUisjc9CiMNDQ1IT09HcnLyxReQyZCcnIy0tLR2r7nuuuuQnp5uDh+nTp3C5s2bcdttt3X4PvX19dDr9W0eZDmtM2oyCyvQZDR1eu7FdUZsJ43IZIK5d4QDWYmI7E+PwkhZWRmMRiMCAgLaHA8ICIBOp2v3mvvvvx8vv/wyrr/+eri4uCA6Ohrjx4/v9DZNSkoKtFqt+REWFtaTMqmHYvw94KlSoLbRiCydofOTbWQF1sv9ISEUcpmAA/kXkF3cRRuIiMimWHw2zY4dO/Dqq6/igw8+QEZGBtavX49NmzbhlVde6fCaRYsWobKy0vwoLOTAREuSyQSM6Oa4EXPPiI2FEX+NGslD/AGwd4SIyN70KIz4+vpCLpejuLi4zfHi4mIEBga2e80LL7yABx54AA8//DCGDRuGu+66C6+++ipSUlJgMrV/S0ClUkGj0bR5kGUl9O/epnnmMSM2dJumVevmed8cLEJdIweyEhHZix6FEaVSiYSEBKSmppqPmUwmpKamIikpqd1rampqIJO1fRu5XA7g4qZrJL2R4V4AutEzcjGN2JxxMX4I8XJFRU0jfjja/m1DIiKyPT2+TbNw4UIsW7YMn376KY4fP4558+ahuroas2fPBgDMnDkTixYtMp8/ZcoUfPjhh1i9ejXy8vKwdetWvPDCC5gyZYo5lJD0hod5QRCAMxdqUaLveA0YW9korz1ymYA/jAoFwFs1RET2RNHTC6ZNm4bS0lK8+OKL0Ol0GD58OLZs2WIe1FpQUNCmJ+T555+HIAh4/vnnUVRUBD8/P0yZMgX//Oc/+64VdNU81S4YFOCJLJ0BGQUXMCk2qN3zRBvaKK89944Kw7upOfj1VDnyyqoR2bKgGxER2a4ehxEAeOyxx/DYY4+1+9yOHTvavoFCgcWLF2Px4sW9eSuyopHh/ZClMyA9v5Mw0vLVNqMIEOzlihsH+mHHiVKs2V+IZyYPlrokIiLqAvemIbOEbuzgK9ro1N5L3Te6eSDruvQzaOxi3RQiIpIewwiZta7EeqRIj/qmzmej2HAWwYQh/vD1UKGsqh6px4u7voCIiCTFMEJm4T5u8HZXosFowpGi9le9tfUxIwDgIpddMpCVa9QQEdk6hhEyEwTBvDR8Rge3amxxo7z2TBvVvGrvzzmlOHOhRuJqiIioMwwj1IZ507wOFj8TbX0Ea4sIX3ckRflAFIG1B85IXQ4REXWCYYTaGNnfC0DzINb2FqWzxY3yOnLfmObekbUHCmE0cYE9IiJbxTBCbcSFekEhE1BiqEdRRe0Vz18cM2Llwnph4jWB8HJzwdnKOvycXSp1OURE1AGGEWrDVSnH0ODmvYDam+JrL2NGAEDtIsfdI7giKxGRrWMYoSt0NojVnnpGgIu3alKzSjpd5p6IiKTDMEJXuDiItaLDc+xhzAgADAzwxMj+XjCaRKzL4EBWIiJbxDBCVxjZEkaOndOjpqGpzXP2sALr5e4b07wi65r9hTBxICsRkc1hGKErBGvVCNSoYTSJOHSmss1z9nabBgB+FxcET5UC+edr8Oup81KXQ0REl2EYoSsIgoCR4V4ArhzEerFfwX7SiJtSgTuGBwMAvtrPFVmJiGwNwwi1q6NBrPbYMwIA01tu1fxwRIfy6gaJqyEioksxjFC7Ll2J9dLFz+xpau+lYkO0uCZYgwajCd8cLJK6HCIiugTDCLXrmmAtlAoZLtQ0Iq+s2nzcXntGgIsDWVfvK2h3dVkiIpIGwwi1S6mQIS5EC6DtFF97Wg7+clOHB8PVRY6ckqoO994hIiLrYxihDrXeqmkziNUOp/a20qhdcHtcEADgq30cyEpEZCsYRqhDI9oZxGonm/Z2aHrLiqzfHToLfV2jxNUQERHAMEKdaJ3em11iMP/ivjhmxD7jyMj+/RDj74G6RhO+zTwrdTlERASGEeqEv6caYd6uEEUgs2XciAj7HvgpCAKmjW7uHVm9n5vnERHZAoYR6lRC/4tTfAH7nk3T6u6RoVDKZThSpMeRosquLyAiIotiGKFOXT6I1Z5n07TydldiYmwgAOCrfewdISKSGsMIdap1EGtmQQVMJtEhekYAYHrLrZpvM89esRkgERFZF8MIdWpwoCfclHIY6puQU1JltyuwXu7aKB+E+7jBUN+ETYfOSV0OEZFTYxihTinkMsSHegFouVXjID0jMpmAe0e1DmTlmiNERFJiGKEutdmnpuWYvU7tvdQfEkIhlwlIz7+A7GKD1OUQETkthhHqkjmM5F/cNM/+owjgr1FjwmB/AMBqrshKRCQZhhHq0oj+XgCAU2XVKK9uWbXUEdIIgOktm+etP3gG9U1GiashInJODCPUJS83JaL93AFcXG/Enqf2XmrcQD8EadWoqGnED0eLpS6HiMgpMYxQt4xsmeKbV1YNwP4HsLaSywT8oXUgK9ccISKSBMMIdUvruJFWDpJFAAD3jgqFIAC/nDyP/PPVUpdDROR0GEaoW64IIw6URkL7uWFcjB8ATvMlIpICwwh1S7SfBzRqhfnPjjJmpNX0Mc23ataln0Gj0SRxNUREzoVhhLpFJhPMS8MDjtUzAgAThgTA10OJUkM9tmWVSF0OEZFTYRihbht5aRiRsA5LcJHLcE9CKAAOZCUisjaGEeq2NuNGHK1rBMB9o5vXHNmZXYqzFbUSV0NE5DwYRqjb4sO0kLVkEMeLIkCkrzuujfKGSQS+PsCBrERE1sIwQt3mqXbBwABPAA7ZMQLg4oqsaw+cgdEkdnE2ERH1BYYR6pGxA3wBAN7uSokrsYyJ1wRC6+qCoopa7MoplbocIiKnwDBCPbIgOQYfPZCAaaPDpC7FItQuctw1IgQAN88jIrIWhhHqEU+1CyZeEwiVQi51KRbTeqvmp+PFKDXUS1wNEZHjYxghusygQE+M6O+FJpOIdelnpC6HiMjhMYwQtWN6yzTfNfsLIIocyEpEZEkMI0TtuD0uCB4qBU6fr8Gvp8qlLoeIyKExjBC1w12lwJT4YADA6v1ckZWIyJIYRog60Lp53vdHdKioaZC4GiIix8UwQtSBYSFaDA3SoKHJhPUZRVKXQ0TksBhGiDogCIK5d2Q1B7ISEVkMwwhRJ6aOCIGrixzZxVXYf/qC1OUQETkkhhGiTmjULrhzRPNA1k/TTktbDBGRg2IYIerCA9dGAAB+OKJDib5O2mKIiBwQwwhRF4YGazAqvB+aTCJW7eM0XyKivsYwQtQNDySFAwBW7S1Ao9EkcTVERI6FYYSoGybHBsHXQ4USQz1+PFosdTlERA6FYYSoG5QKGe5vmeb7GQeyEhH1KYYRom66PzEccpmAvXnlyNLppS6HiMhhMIwQdVOgVo1bhwYAAD5Py5e4GiIix8EwQtQDrQNZvzlYBH1do8TVEBE5BoYRoh5IivJBjL8HahqMWJ9+RupyiIgcAsMIUQ8IgoCZLb0jn/+az/1qiIj6AMMIUQ/dNTIUHioFTpZW45eT56Uuh4jI7jGMEPWQh0qBu0eGAAA+/eW0tMUQETmAXoWRJUuWICIiAmq1GomJidi3b1+n51dUVGD+/PkICgqCSqXCwIEDsXnz5l4VTGQLWm/V/HS8GIXlNRJXQ0Rk33ocRtasWYOFCxdi8eLFyMjIQHx8PCZOnIiSkpJ2z29oaMAtt9yC06dPY926dThx4gSWLVuGkJCQqy6eSCoD/D0xbqAfTCKwkr0jRERXRRB7OAIvMTERo0ePxvvvvw8AMJlMCAsLw+OPP45nnnnmivOXLl2KN954A1lZWXBxcenWe9TX16O+vt78Z71ej7CwMFRWVkKj0fSkXCKL2XGiBA+u2A8PlQJpi26Gp7p7P99ERM5Cr9dDq9V2+fu7Rz0jDQ0NSE9PR3Jy8sUXkMmQnJyMtLS0dq/59ttvkZSUhPnz5yMgIACxsbF49dVXYTQaO3yflJQUaLVa8yMsLKwnZRJZxbgYP0T7uaOqvgnrOM2XiKjXehRGysrKYDQaERAQ0OZ4QEAAdDpdu9ecOnUK69atg9FoxObNm/HCCy/gzTffxD/+8Y8O32fRokWorKw0PwoLC3tSJpFVyGQCZo+NBNB8q8Zo4jRfIqLesPhsGpPJBH9/f/z73/9GQkICpk2bhueeew5Lly7t8BqVSgWNRtPmQWSL7h4ZAq2rC/LP12BbVvvjpoiIqHM9CiO+vr6Qy+UoLm67hXpxcTECAwPbvSYoKAgDBw6EXC43HxsyZAh0Oh0aGhp6UTKR7XBTKjB9TH8AwPLdeRJXQ0Rkn3oURpRKJRISEpCammo+ZjKZkJqaiqSkpHavGTt2LHJzc2EymczHsrOzERQUBKVS2cuyiWzHzKTm3XzTTp3HsbPczZeIqKd6fJtm4cKFWLZsGT799FMcP34c8+bNQ3V1NWbPng0AmDlzJhYtWmQ+f968eSgvL8cTTzyB7OxsbNq0Ca+++irmz5/fd60gklCwlysmxzb3DC7fw94RIqKeUvT0gmnTpqG0tBQvvvgidDodhg8fji1btpgHtRYUFEAmu5hxwsLC8MMPP+DJJ59EXFwcQkJC8MQTT+Dpp5/uu1YQSeyh6yPx3aFz+DbzLJ6eNBh+niqpSyIishs9XmdECt2dp0wkpTuX7EFmYQUWJMdgQfJAqcshIpKcRdYZIaKOPXR98zTfL37NR31Tx+voEBFRWwwjRH1kcmwggrRqlFU1YOPBs1KXQ0RkNxhGiPqIi1yG2WMjAAD/3nUKJi6CRkTULQwjRH1o+pj+8FQpkFtSxUXQiIi6iWGEqA95ql1w/7XNi6D9++dTEldDRGQfGEaI+thDYyPhIhew73Q5MgouSF0OEZHNYxgh6mMBGjXuHB4CAPj3TvaOEBF1hWGEyAIeGRcFAPjhmA6nSqskroaIyLYxjBBZQEyAJyYM9ocoAst2cYl4IqLOMIwQWUhr78h/Ms6g1FAvcTVERLaLYYTIQsZEemN4mBcamkz49JfTUpdDRGSzGEaILEQQBPyppXfk81/zUV3fJHFFRES2iWGEyIJuvSYQkb7uqKxtxKq9BVKXQ0RkkxhGiCxILhPw6I3NvSP/3nUKdY3cQI+I6HIMI0QWdteIUARr1Sg11GNt+hmpyyEisjkMI0QWplTI8KcbowEAS3ecRKPRJHFFRES2hWGEyAqmjQ6Dr4cKRRW12HCwSOpyiIhsCsMIkRWoXeSYe0MkAOCDHSdhNIkSV0REZDsYRoisZMa14fByc0FeWTU2Hz4ndTlERDaDYYTISjxUCjw0trl3ZMn2XJjYO0JEBIBhhMiqZiVFwEOlQJbOgNSsEqnLISKyCQwjRFakdXPBzKRwAMD723MhiuwdISJiGCGysjnXR0LtIsNvhRXYmV0qdTlERJJjGCGyMh8PFR64trl35O2fctg7QkROj2GESAJ/ujEari5y/FZYge0nOHaEiJwbwwiRBHw9VJh1XQQA4K2t2ewdISKnxjBCJJFHxkXBXSnHkSI9fjxWLHU5RESSYRghkoi3uxKzW9YdeXtrNtcdISKnxTBCJKGHb4iEZ8u6I1uO6qQuh4hIEgwjRBLyclPioesv9o5wzxoickYMI0QSe+j6SGjUCuSUVOG7Q2elLoeIyOoYRogkpnV1wdwbogAA7/6UgyajSeKKiIisi2GEyAY8ODYCXm4uOFVWjfUZRVKXQ0RkVQwjRDbAU+2C+eMHAADe/ikbdY1GiSsiIrIehhEiG/FAUjiCtWqcq6zDZ2mnpS6HiMhqGEaIbITaRY4FtwwEACzZfhKVtY0SV0REZB0MI0Q25J6RoYjx90BlbSM+2nlS6nKIiKyCYYTIhshlAv42aTAAYPmePBTr6ySuiIjI8hhGiGxM8hB/jArvh7pGE95NzZG6HCIii2MYIbIxgiDg6cnNvSNr9hfiVGmVxBUREVkWwwiRDRod4Y0Jg/1hNIl47fssqcshIrIohhEiG/XM5MGQywT8eKwYO7NLpS6HiMhiGEaIbFRMgCcevC4CAPDSt0fR0MRl4onIMTGMENmwJ5Jj4Ouhwqmyaizfkyd1OUREFsEwQmTDNGoXLGoZzPpeag50lZzqS0SOh2GEyMbdNSIECeH9UNNgxD83H5e6HCKiPscwQmTjZDIBL0+9BjIB+O9vZ5F28rzUJRER9SmGESI7cE2wFjMSwwEAi789gkYjB7MSkeNgGCGyE3+5dSD6ubkgu7gKy3dzMCsROQ6GESI74eWmxLO3DQEAvP1TNgrLaySuiIiobzCMENmR3yeE4toob9Q1mvD8hiMQRVHqkoiIrhrDCJEdEQQB/7xrGJRyGXZml+K7Q+ekLomI6KoxjBDZmWg/D8y/aQAA4KX/HkNlTaPEFRERXR2GESI79Oj4KET7uaOsqh6vbeFGekRk3xhGiOyQSiHHq3cNAwB8ta8A+0+XS1wREVHvMYwQ2anEKB9MGxUGAHh2/WHUNxklroiIqHcYRojs2KLbBsPXQ4mckios2ZYrdTlERL3CMEJkx7zclHh5aiwAYMmOkzhSVClxRUREPccwQmTnbhsWhNuHBcFoEvHU2t/Q0MSl4onIvjCMEDmAl6deA293JbJ0Bry/nbdriMi+MIwQOQAfDxVeabld88H2XN6uISK7wjBC5CBujwvCbcMC0cTbNURkZ3oVRpYsWYKIiAio1WokJiZi37593bpu9erVEAQBd955Z2/eloi68PLUWPPtmiW8XUNEdqLHYWTNmjVYuHAhFi9ejIyMDMTHx2PixIkoKSnp9LrTp0/jqaeewg033NDrYomoc74eKrw89RoAwJLtuTh0pkLagoiIuqHHYeStt97C3LlzMXv2bAwdOhRLly6Fm5sbli9f3uE1RqMRM2bMwEsvvYSoqKgu36O+vh56vb7Ng4i65/aW2TVNJhELVmeipqFJ6pKIiDrVozDS0NCA9PR0JCcnX3wBmQzJyclIS0vr8LqXX34Z/v7+mDNnTrfeJyUlBVqt1vwICwvrSZlETq15Z99YBGrUOFVWjX9uOi51SUREnepRGCkrK4PRaERAQECb4wEBAdDpdO1es3v3bnzyySdYtmxZt99n0aJFqKysND8KCwt7UiaR0/NyU+LNe+MBAF/uLUDq8WKJKyIi6phFZ9MYDAY88MADWLZsGXx9fbt9nUqlgkajafMgop4ZO8AXD18fCQD427pDKDXUS1wREVH7FD052dfXF3K5HMXFbf8vq7i4GIGBgVecf/LkSZw+fRpTpkwxHzOZmqcbKhQKnDhxAtHR0b2pm4i64a+TBmF3bhmydAY8/Z9D+GTWKAiCIHVZRERt9KhnRKlUIiEhAampqeZjJpMJqampSEpKuuL8wYMH4/Dhw8jMzDQ/7rjjDtx0003IzMzkWBAiC1Mp5Hj3vhFQKmTYllWCL/YWSF0SEdEVetQzAgALFy7ErFmzMGrUKIwZMwbvvPMOqqurMXv2bADAzJkzERISgpSUFKjVasTGxra53svLCwCuOE5EljEo0BPPTBqMl787hn9uOobESG8MDPCUuiwiIrMeh5Fp06ahtLQUL774InQ6HYYPH44tW7aYB7UWFBRAJuPCrkS25MHrIrAzuxQ7s0vx5y8z8O1jY+Gm7PFffyIiixBEURSlLqIrer0eWq0WlZWVHMxK1Evnq+px23u7UKyvxz0jQ82zbYiILKW7v7/ZhUHkJHw8VHjvvhGQCcB/Ms5g7QFOmSci28AwQuREEqN8sPCWgQCAFzYeQXaxQeKKiIgYRoiczp/HD8ANMb6oazRh/pcZXC6eiCTHMELkZGQyAW9PGw5/TxVySqrwwoajsIOhY0TkwBhGiJyQr4cK714yfmTVPq4/QkTSYRghclJJ0T54auIgAMDfvz2K9PwLEldERM6KYYTIic27MRqTYwPRaBQx74t0lOjrpC6JiJwQwwiRExMEAW/8IR4x/h4oMdTjz19moKHJJHVZRORkGEaInJyHSoGPHkiAp0qBA/kX8I9Nx6QuiYicDMMIESHKzwPv3DccAPBZWj4XRCMiq2IYISIAwIQhAXhiQgwA4LkNRziglYishmGEiMyemBCDW4YGoKHJhD99fgCF5TVSl0REToBhhIjMZDIB70wbjqFBGpRVNeDhTw/AUNcodVlE5OAYRoioDXeVAp88OAp+niqcKDbgf746CKOJK7QSkeUwjBDRFYK0rvh45iioFDJsP1GKf246LnVJROTAGEaIqF3xYV54697hAIDle/Lwxa/50hZERA6LYYSIOnR7XBD+cstAAMDib49iW1axxBURkSNiGCGiTj128wDcMzIURpOIP3+ZgYMFnPJLRH2LYYSIOiUIAl67ZxhuHOiHukYTHlq5H6dKq6Qui4gcCMMIEXXJRS7DBzNGIi5Uiws1jZi5fB9KDNxUj4j6BsMIEXWLu0qB5Q+ORriPG85cqMXsFftRVd8kdVlE5AAYRoio23w9VPjsoTHwcVfi6Fk9Hv08HXWNRqnLIiI7xzBCRD0S7uOOFbNHw00px+7cMjz+1UE0Gk1Sl0VEdoxhhIh6LC7UCx/PHAWlQoatx4rx1NrfuEorEfUawwgR9cp1A3yx9I8joZAJ2Jh5Fs9vOAxRZCAhop5jGCGiXrt5cADeuW84ZALw1b5C/HPTcQYSIuoxhhEiuiq/iwvGa/fEAQA+3p2Ht3/KkbgiIrI3DCNEdNXuHRWGxVOGAgDeS83B21uzJa6IiOwJwwgR9YnZYyPx7G2DAQDvpubgra3ZvGVDRN3CMEJEfeaRcdF47rYhAC72kDCQEFFXGEaIqE/NHReF529vCSTbcvHmjwwkRNQ5hhEi6nMP3xCFF37XPIbk/e25eH3LCQYSIuoQwwgRWcSc6yPxYksgWbrzJF7YeAQmLoxGRO1gGCEii3no+kj8865YCALwxa8FePLrTC4dT0RXYBghIouakRiOd+8bYV6plZvrEdHlGEaIyOLuiA/GspmjoFLIkJpVglnL98FQ1yh1WURkIxhGiMgqbhrsj8/nJMJTpcDevHLcv2wvSg31UpdFRDaAYYSIrGZMpDe+euRaeLsrcbioEnd9sAe5JQapyyIiiTGMEJFVxYZose7RJIT7uOHMhVrc/cEv+OVkmdRlEZGEGEaIyOqi/DzwzZ/HIiG8H/R1TXjgk334cMdJTv0lclIMI0QkCW93Jb58OBF3jwiB0STi9S1ZePizA7hQ3SB1aURkZQwjRCQZtYscb94bj9fuHgalQoZtWSW4/b1dSM+/IHVpRGRFDCNEJClBEHDfmP7Y8OexiPR1x9nKOkz7KA0f7zrFJeSJnATDCBHZhKHBGnz72FjcHheEJpOIf2w6jrmfHcD5Kk7/JXJ0DCNEZDM81S54f/oIvDL1GijlMvx0vAQT39mFHSdKpC6NiCyIYYSIbIogCHggKQIb5o9FjL8Hyqrq8eCK/fj7t0e5jDyRg2IYISKbNDRYg/8+fj1mJYUDAFb+chp3vL8bx8/pJa6MiPoawwgR2Sy1ixwvTY3Fitmj4euhQnZxFaa+vwcf7TwJI9ckIXIYDCNEZPNuGuSPLQtuQPIQfzQYTUj5Pgu/X/oLckuqpC6NiPoAwwgR2QVfDxWWzRyF1+8ZBk+VAgcLKnDbe7uwdOdJNBlNUpdHRFeBYYSI7IYgCJg2uj9+eHIcbhzoh4YmE177Pgv3LE1DTjE33COyVwwjRGR3gr1csXL2aPzv7+PgqVbgt8IK3P7ebizZnotG9pIQ2R2GESKyS4Ig4N5RYfjxyXG4aZAfGowmvPHDCdz+3i4cOF0udXlE1AMMI0Rk14K0rlj+4Gi8dW88vN2VyC6uwu+XpuGZ/xxCRQ033SOyBwwjRGT3BEHA3SNDkbrwRkwbFQYAWL2/EBPe3In1GWe4xw2RjWMYISKH0c9didd/H4ev/5SEGH8PnK9uwMKvf8OMj/dygCuRDRNEO/hfBr1eD61Wi8rKSmg0GqnLISI70NBkwrJdp/Beag7qm0yQywQ8cG04nkweCK2bi9TlETmF7v7+Zs8IETkkpUKG+TcNwNYnb8QtQwNgNIlY+ctp3PTmDny5N58ruBLZEPaMEJFT2JVTipf/eww5Lau2DgnS4O9ThiIxykfiyogcV3d/fzOMEJHTaDSa8MWv+Xh7azb0dU0AgMmxgfjrxEGI8vOQuDoix8MwQkTUgfNV9Xhraza+2lcAkwjIZQLuH9Mf/zMhBn6eKqnLI3IYDCNERF04oTPg9S1Z2JZVAgBwV8oxd1wU5t4QBXeVQuLqiOyfRQewLlmyBBEREVCr1UhMTMS+ffs6PHfZsmW44YYb0K9fP/Tr1w/Jycmdnk9EZC2DAj2x/MHRWP3ItYgP1aK6wYh3fsrBjW/swOe/5nNpeSIr6XEYWbNmDRYuXIjFixcjIyMD8fHxmDhxIkpKSto9f8eOHZg+fTq2b9+OtLQ0hIWF4dZbb0VRUdFVF09E1BeujfLBhvlj8f79IxDu44ayqnq8sOEIbn5zB74+UMhdgYksrMe3aRITEzF69Gi8//77AACTyYSwsDA8/vjjeOaZZ7q83mg0ol+/fnj//fcxc+bMbr0nb9MQkbU0NJnw1b4C/GtbLsqq6gEAET5u+J8JMbgjPhgKOVdEIOoui9ymaWhoQHp6OpKTky++gEyG5ORkpKWldes1ampq0NjYCG9v7w7Pqa+vh16vb/MgIrIGpUKGWddFYNffbsJztw2Bj7sSp8/XYOHXv+HWt3/GxswirlFC1Md6FEbKyspgNBoREBDQ5nhAQAB0Ol23XuPpp59GcHBwm0BzuZSUFGi1WvMjLCysJ2USEV0115bBrD//7SY8PWkwvNxccKqsGk+szsTEd5pDCW/fEPUNq/Y3vvbaa1i9ejW++eYbqNXqDs9btGgRKisrzY/CwkIrVklEdJG7SoF546Ox++mb8deJg6B1dUFuSRWeWJ2Jm95sHuha12iUukwiu9ajMOLr6wu5XI7i4uI2x4uLixEYGNjptf/3f/+H1157DT/++CPi4uI6PVelUkGj0bR5EBFJyUOlwPybBmDX0zfhqVsHwttdicLyWryw4Qiuf307PtxxEoa6RqnLJLJLPQojSqUSCQkJSE1NNR8zmUxITU1FUlJSh9f97//+L1555RVs2bIFo0aN6n21REQS06hd8NjNMdjz9M34+5ShCPFyRVlVPV7fkoXrXtuGN37IQqmhXuoyiexKj2fTrFmzBrNmzcJHH32EMWPG4J133sHXX3+NrKwsBAQEYObMmQgJCUFKSgoA4PXXX8eLL76IVatWYezYsebX8fDwgIdH95Zf5mwaIrJVjUYTvs08iw93nkRuy743SrkMU4cHY/bYSAwN5n+zyHlZdAXW999/H2+88QZ0Oh2GDx+O9957D4mJiQCA8ePHIyIiAitXrgQAREREID8//4rXWLx4Mf7+97/3aWOIiKRiMonYerwYH+44iczCCvPxpCgfzLk+EjcP9odMJkhXIJEEuBw8EZFEMgouYPnuPHx/RGeeBhzh44bZYyPx+4RQLjVPToNhhIhIYkUVtfgs7TS+2ltg3iXYU63APSND8cdr+2OAv6fEFRJZFsMIEZGNqGlown/Sz2DFntM4VVZtPp4Y6Y0Z14Zj0jWBUCq4sis5HoYRIiIbYzKJ+DmnFF/uLUDq8WK0LuTq66HEH0aF4f4x/RHm7SZtkUR9iGGEiMiGnausxVf7CrF6XwFKWqYCCwIwLsYP944KQ/JQf6gUcomrJLo6DCNERHag0WhC6vFifLm3ALtyyszHvdxcMDU+GH8YFYbYEK2EFRL1HsMIEZGdOV1Wja8PFGJ9RhF0+jrz8SFBGvw+IRR3Dg+Gj4dKwgqJeoZhhIjIThlNInbllGJt+hlsPVqMhpYN+VzkAm4e7I+7RoRg/CB/qF14G4dsG8MIEZEDqKhpwLe/ncXaA2dwuKjSfNxTpcCk2EDcMTwYSVE+UMg5G4dsD8MIEZGDydLp8c3BIvw38yzOVl68jeProcLv4oIwdXgwhod5QRC40ivZBoYRIiIHZTKJOJB/ARszi7D58DlcqLm4W3B/bzf8Li4Ik2ODEBuiYTAhSTGMEBE5gUajCbtzyrAxswg/HitGTYPR/FxoP1dMuiYQk4cFYkRYP+6NQ1bHMEJE5GRqGprw0/ESfH/4HHacKEVt48VgEqBRYeI1gZgUG4gxEd4cY0JWwTBCROTEahuM2Jldgu+P6LDteAkM9U3m57zdlUge4o8JQwJw/QBfbtxHFsMwQkREAID6JiP25JZhyxEdfjxWjIpLxpgoFTIkRflgwhB/3DzYH6H9uBw99R2GESIiukKT0YS9eeX46XgxUo+XoKC8ps3zgwM9W4JJAIaHeUHOcSZ0FRhGiIioU6Io4mRpFX46XoJtx0twIL/cvHkf0Hw754YYX4yL8cMNMb7w16ilK5bsEsMIERH1yIXqBuzMLsVPx4uxM7sUhrqmNs8PDvTEuIF+GBfjh1ER/bgCLHWJYYSIiHqt0WhCev4F7Mopxc/ZZThythKX/rZQu8iQGOmDG2J8cV20LwYGeHCGDl2BYYSIiPrM+ap67M4tw66cMuzKKUWxvr7N825KOeJCtRge1g/Dw7QYFuqFYK2ai645OYYRIiKyCFEUkV1chZ+zS/FzTikOFlSgqr7pivN83JUYFqpFXEhzOIkL1SKA406cCsMIERFZhdHUPBD2YMEFHCyowOGiSpzQGdBkuvLXi7+nCnGhWgwLaQ4nsSFa+HmqJKiarIFhhIiIJFPXaESWzoDDZypw6EwlDhdVIrvYgHbyCQI0KgwJ0mBIkAZDW75G+rpzWrEDYBghIiKbUttgxLFzlc3hpCWg5JZWob3fQmoXGQYFNoeToUGeGBqswaBADTy4WqxdYRghIiKbV13fhCydAcfP6XHsnB7Hz+mRdc7QZl+dS4X7uGFokAaDAj0xMMATAwM8EO7jDhfO5LFJDCNERGSXjCYR+eerzeHk2Fk9jp8zQKeva/d8F7mAaD8PxAR4YqB/y9eWkMJbPdJiGCEiIodSXt1gDifZxQZkFxuQU1KFmob2e1FUChmi/TwwMKA5oAwK8MQAfw+E9nPlmihWwjBCREQOz2QSUVRRi5wSA07oqpBTbEB2iQG5JVWoazS1e42LXEC4jzuifN0R7e+BKF93RPl5INrPHV5uSiu3wLExjBARkdMymkQUlteYe0+yiw04oTMgr6wa9U3thxSgeW2UKD93RPl6IMrPHdF+zV/DvN04LqUXGEaIiIguYzKJOFtZi5Ol1ThVWoVTpdU4Vdb89Vxl+2NSAEAhExDm7Yb+3m6I8HFDfx93RPi4IdzHDaH93LhPTwcYRoiIiHqgur4JeWXVOGkOKRcDS0ezewBAEIAgjRrhPu4I93G75Gvz9848HZlhhIiIqA+YTCJ0+jrkn69BQXk1Tp+vQcH5Gpw+X4388zXtLoV/KR93JcJ93BDm7YbQfq4I7Xfxa7CXGiqF4/aqMIwQERFZmCiKKK9uQH55DfJbwknzo/n789UNXb5GgEZ1SUBxrLDCMEJERCQxQ11jS49KDYou1OLMhRqcuVCLwpavHU1LbiUIzfv5hPZzQ4iXK4K81AjWuiJIq0awV/NXb3elze6OzDBCRERkw0RRxIWaRnNAufi1+fvC8tpOx6q0UipkCNKqmwOKtjmwBLUEliCtK4K91NC6ukgSWLr7+9t5R9UQERFJSBAEeLsr4e2uRFyo1xXPXx5Wii7U4mxlLc5V1OFcZS3OVtahrKoeDU0m8+2hjri6yFtCSnNACdSoEaBVI1CjbvleBV93FWQSrVjLMEJERGSDugorANDQZEKxvg7nKlsCSkXbr7rKOpyvbkBto7F5hlBpdYfv98WcRFwf42uh1nSOYYSIiMhOKRUyhHk3z9TpSF2jEbrKOnOvik5fB11l89filu/LquoRoFFZsfK2GEaIiIgcmNpFjghfd0T4und4TpPRBJmEg2AZRoiIiJyc1BsHcqF9IiIikhTDCBEREUmKYYSIiIgkxTBCREREkmIYISIiIkkxjBAREZGkGEaIiIhIUgwjREREJCmGESIiIpIUwwgRERFJimGEiIiIJMUwQkRERJJiGCEiIiJJ2cWuvaIoAgD0er3ElRAREVF3tf7ebv093hG7CCMGgwEAEBYWJnElRERE1FMGgwFarbbD5wWxq7hiA0wmE86ePQtPT08IgtAnr6nX6xEWFobCwkJoNJo+eU1b4+htdPT2AY7fRkdvH+D4bXT09gGO30ZLtk8URRgMBgQHB0Mm63hkiF30jMhkMoSGhlrktTUajUP+cF3K0dvo6O0DHL+Njt4+wPHb6OjtAxy/jZZqX2c9Iq04gJWIiIgkxTBCREREknLaMKJSqbB48WKoVCqpS7EYR2+jo7cPcPw2Onr7AMdvo6O3D3D8NtpC++xiACsRERE5LqftGSEiIiLbwDBCREREkmIYISIiIkkxjBAREZGknDaMLFmyBBEREVCr1UhMTMS+ffukLqlX/v73v0MQhDaPwYMHm5+vq6vD/Pnz4ePjAw8PD9xzzz0oLi6WsOKu/fzzz5gyZQqCg4MhCAI2bNjQ5nlRFPHiiy8iKCgIrq6uSE5ORk5OTptzysvLMWPGDGg0Gnh5eWHOnDmoqqqyYis61lX7HnzwwSs+00mTJrU5x5bbl5KSgtGjR8PT0xP+/v648847ceLEiTbndOfnsqCgALfffjvc3Nzg7++Pv/71r2hqarJmUzrUnTaOHz/+is/x0UcfbXOOrbbxww8/RFxcnHkRrKSkJHz//ffm5+398wO6bqM9f37tee211yAIAhYsWGA+ZlOfo+iEVq9eLSqVSnH58uXi0aNHxblz54peXl5icXGx1KX12OLFi8VrrrlGPHfunPlRWlpqfv7RRx8Vw8LCxNTUVPHAgQPitddeK1533XUSVty1zZs3i88995y4fv16EYD4zTfftHn+tddeE7Varbhhwwbxt99+E++44w4xMjJSrK2tNZ8zadIkMT4+Xvz111/FXbt2iQMGDBCnT59u5Za0r6v2zZo1S5w0aVKbz7S8vLzNObbcvokTJ4orVqwQjxw5ImZmZoq33Xab2L9/f7Gqqsp8Tlc/l01NTWJsbKyYnJwsHjx4UNy8ebPo6+srLlq0SIomXaE7bbzxxhvFuXPntvkcKysrzc/bchu//fZbcdOmTWJ2drZ44sQJ8dlnnxVdXFzEI0eOiKJo/5+fKHbdRnv+/C63b98+MSIiQoyLixOfeOIJ83Fb+hydMoyMGTNGnD9/vvnPRqNRDA4OFlNSUiSsqncWL14sxsfHt/tcRUWF6OLiIq5du9Z87Pjx4yIAMS0tzUoVXp3Lf1mbTCYxMDBQfOONN8zHKioqRJVKJX711VeiKIrisWPHRADi/v37zed8//33oiAIYlFRkdVq746OwsjUqVM7vMae2ieKolhSUiICEHfu3CmKYvd+Ljdv3izKZDJRp9OZz/nwww9FjUYj1tfXW7cB3XB5G0Wx+ZfZpf/hv5y9tbFfv37ixx9/7JCfX6vWNoqi43x+BoNBjImJEbdu3dqmTbb2OTrdbZqGhgakp6cjOTnZfEwmkyE5ORlpaWkSVtZ7OTk5CA4ORlRUFGbMmIGCggIAQHp6OhobG9u0dfDgwejfv7/dtjUvLw86na5Nm7RaLRITE81tSktLg5eXF0aNGmU+Jzk5GTKZDHv37rV6zb2xY8cO+Pv7Y9CgQZg3bx7Onz9vfs7e2ldZWQkA8Pb2BtC9n8u0tDQMGzYMAQEB5nMmTpwIvV6Po0ePWrH67rm8ja2+/PJL+Pr6IjY2FosWLUJNTY35OXtpo9FoxOrVq1FdXY2kpCSH/Pwub2MrR/j85s+fj9tvv73N5wXY3t9Du9gory+VlZXBaDS2+ZcLAAEBAcjKypKoqt5LTEzEypUrMWjQIJw7dw4vvfQSbrjhBhw5cgQ6nQ5KpRJeXl5trgkICIBOp5Om4KvUWnd7n1/rczqdDv7+/m2eVygU8Pb2tot2T5o0CXfffTciIyNx8uRJPPvss5g8eTLS0tIgl8vtqn0mkwkLFizA2LFjERsbCwDd+rnU6XTtfsatz9mS9toIAPfffz/Cw8MRHByMQ4cO4emnn8aJEyewfv16ALbfxsOHDyMpKQl1dXXw8PDAN998g6FDhyIzM9NhPr+O2gjY/+cHAKtXr0ZGRgb2799/xXO29vfQ6cKIo5k8ebL5+7i4OCQmJiI8PBxff/01XF1dJayMeuu+++4zfz9s2DDExcUhOjoaO3bswIQJEySsrOfmz5+PI0eOYPfu3VKXYjEdtfGRRx4xfz9s2DAEBQVhwoQJOHnyJKKjo61dZo8NGjQImZmZqKysxLp16zBr1izs3LlT6rL6VEdtHDp0qN1/foWFhXjiiSewdetWqNVqqcvpktPdpvH19YVcLr9ixHBxcTECAwMlqqrveHl5YeDAgcjNzUVgYCAaGhpQUVHR5hx7bmtr3Z19foGBgSgpKWnzfFNTE8rLy+2y3VFRUfD19UVubi4A+2nfY489hu+++w7bt29HaGio+Xh3fi4DAwPb/Yxbn7MVHbWxPYmJiQDQ5nO05TYqlUoMGDAACQkJSElJQXx8PN59912H+vw6amN77O3zS09PR0lJCUaOHAmFQgGFQoGdO3fivffeg0KhQEBAgE19jk4XRpRKJRISEpCammo+ZjKZkJqa2uZeob2qqqrCyZMnERQUhISEBLi4uLRp64kTJ1BQUGC3bY2MjERgYGCbNun1euzdu9fcpqSkJFRUVCA9Pd18zrZt22Aymcz/QbEnZ86cwfnz5xEUFATA9tsniiIee+wxfPPNN9i2bRsiIyPbPN+dn8ukpCQcPny4TejaunUrNBqNuRtdSl21sT2ZmZkA0OZztOU2Xs5kMqG+vt4hPr+OtLaxPfb2+U2YMAGHDx9GZmam+TFq1CjMmDHD/L1NfY59OhzWTqxevVpUqVTiypUrxWPHjomPPPKI6OXl1WbEsL34y1/+Iu7YsUPMy8sT9+zZIyYnJ4u+vr5iSUmJKIrNU7f69+8vbtu2TTxw4ICYlJQkJiUlSVx15wwGg3jw4EHx4MGDIgDxrbfeEg8ePCjm5+eLotg8tdfLy0vcuHGjeOjQIXHq1KntTu0dMWKEuHfvXnH37t1iTEyMzUx97ax9BoNBfOqpp8S0tDQxLy9P/Omnn8SRI0eKMTExYl1dnfk1bLl98+bNE7Varbhjx4420yJramrM53T1c9k6pfDWW28VMzMzxS1btoh+fn42M22yqzbm5uaKL7/8snjgwAExLy9P3LhxoxgVFSWOGzfO/Bq23MZnnnlG3Llzp5iXlyceOnRIfOaZZ0RBEMQff/xRFEX7//xEsfM22vvn15HLZwjZ0ufolGFEFEXxX//6l9i/f39RqVSKY8aMEX/99VepS+qVadOmiUFBQaJSqRRDQkLEadOmibm5uebna2trxT//+c9iv379RDc3N/Guu+4Sz507J2HFXdu+fbsI4IrHrFmzRFFsnt77wgsviAEBAaJKpRInTJggnjhxos1rnD9/Xpw+fbro4eEhajQacfbs2aLBYJCgNVfqrH01NTXirbfeKvr5+YkuLi5ieHi4OHfu3CuCsi23r722ARBXrFhhPqc7P5enT58WJ0+eLLq6uoq+vr7iX/7yF7GxsdHKrWlfV20sKCgQx40bJ3p7e4sqlUocMGCA+Ne//rXNOhWiaLttfOihh8Tw8HBRqVSKfn5+4oQJE8xBRBTt//MTxc7baO+fX0cuDyO29DkKoiiKfdvXQkRERNR9TjdmhIiIiGwLwwgRERFJimGEiIiIJMUwQkRERJJiGCEiIiJJMYwQERGRpBhGiIiISFIMI0RERCQphhEiIiKSFMMIEfW58ePHY8GCBVKXQUR2gmGEiGzKlClTMGnSpHaf27VrFwRBwKFDh6xcFRFZEsMIEdmUOXPmYOvWrThz5swVz61YsQKjRo1CXFycBJURkaUwjBCRxW3atAlarRZffvklAMBkMiElJQWRkZFwdXVFfHw81q1bBwD43e9+Bz8/P6xcubLNa1RVVWHt2rWYM2eOtcsnIgtjGCEii1q1ahWmT5+OL7/8EjNmzAAApKSk4LPPPsPSpUtx9OhRPPnkk/jjH/+InTt3QqFQYObMmVi5ciUu3VR87dq1MBqNmD59ulRNISILEcRL/7YTEfWB8ePHY/jw4YiJicFzzz2HjRs34sYbbwQA1NfXw9vbGz/99BOSkpLM1zz88MOoqanBqlWrkJWVhSFDhmD79u0YP348AGDcuHEIDw/H559/LkWTiMiCGEaIqM+NHz8eubm5KCkpwZ49ezB69Gjzc0ePHkVsbCzc3d3bXNPQ0IARI0Zg7969AICxY8ciOjoan332GXJzcxETE9MmnBCR41BIXQAROaYRI0YgIyMDy5cvx6hRoyAIAoDmsR9A8ziSkJCQNteoVCrz93PmzMHjjz+OJUuWYMWKFYiOjjb3rhCRY+GYESKyiOjoaGzfvh0bN27E448/bj4+dOhQqFQqFBQUYMCAAW0eYWFh5vPuvfdeyGQyrFq1Cp999hkeeughc6AhIsfCnhEispiBAweab60oFAq888478PT0xFNPPYUnn3wSJpMJ119/PSorK7Fnzx5oNBrMmjULAODh4YFp06Zh0aJF0Ov1ePDBB6VtDBFZDMMIEVnUoEGDsG3bNowfPx5yuRxvvvkmXnnlFfj5+SElJQWnTp2Cl5cXRo4ciWeffbbNtXPmzMEnn3yC2267DcHBwRK1gIgsjQNYiYiISFIcM0JERESSYhghIiIiSTGMEBERkaQYRoiIiEhSDCNEREQkKYYRIiIikhTDCBEREUmKYYSIiIgkxTBCREREkmIYISIiIkkxjBAREZGk/h+QhL3Ezh+lVgAAAABJRU5ErkJggg==\n", "text/plain": [ "
" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "with quantity_support():\n", " plt.plot(energy, alloy.absorption(energy))" ] }, { "cell_type": "code", "execution_count": null, "id": "bc1fe899-4cc9-4b3e-9aad-ca39d88e83da", "metadata": {}, "outputs": [], "source": [] } ], "metadata": { "kernelspec": { "display_name": "Python [conda env:mocksipipeline]", "language": "python", "name": "conda-env-mocksipipeline-py" }, "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" } }, "nbformat": 4, "nbformat_minor": 5 }