- {
- “cells”: [
- {
“cell_type”: “markdown”, “id”: “68864424”, “metadata”: {
- “papermill”: {
“duration”: 0.012029, “end_time”: “2021-08-17T08:11:28.499381”, “exception”: false, “start_time”: “2021-08-17T08:11:28.487352”, “status”: “completed”
}, “tags”: []
}, “source”: [
“# PhAbs”
]
}, {
“cell_type”: “code”, “execution_count”: 1, “id”: “5fba72c6”, “metadata”: {
- “execution”: {
“iopub.execute_input”: “2021-08-17T08:11:28.532115Z”, “iopub.status.busy”: “2021-08-17T08:11:28.530945Z”, “iopub.status.idle”: “2021-08-17T08:11:31.499436Z”, “shell.execute_reply”: “2021-08-17T08:11:31.497920Z”
}, “nbsphinx”: “hidden”, “papermill”: {
“duration”: 2.990727, “end_time”: “2021-08-17T08:11:31.499777”, “exception”: false, “start_time”: “2021-08-17T08:11:28.509050”, “status”: “completed”
}, “tags”: []
}, “outputs”: [], “source”: [
“%%capturen”, “n”, “import numpy as npn”, “n”, “import matplotlib.pyplot as pltn”, “n”, “import warningsn”, “warnings.simplefilter("ignore")n”, “n”, “from astromodels.functions.function import _known_functionsn”, “n”, “n”, “from jupyterthemes import jtplotn”, “jtplot.style(context="talk", fscale=1, ticks=True, grid=False)n”, “%matplotlib inline”
]
}, {
“cell_type”: “code”, “execution_count”: 2, “id”: “c8be52df”, “metadata”: {
- “execution”: {
“iopub.execute_input”: “2021-08-17T08:11:31.527350Z”, “iopub.status.busy”: “2021-08-17T08:11:31.524993Z”, “iopub.status.idle”: “2021-08-17T08:11:31.532682Z”, “shell.execute_reply”: “2021-08-17T08:11:31.533533Z”
}, “nbsphinx”: “hidden”, “papermill”: {
“duration”: 0.023838, “end_time”: “2021-08-17T08:11:31.533855”, “exception”: false, “start_time”: “2021-08-17T08:11:31.510017”, “status”: “completed”
}, “tags”: [
“parameters”
]
}, “outputs”: [], “source”: [
“func_name = "TbAbs"n”, “n”, “x_scale="log"n”, “y_scale="log"n”, “n”, “linear_range = Falsen”, “n”, “wide_energy_range = False”
]
}, {
“cell_type”: “code”, “execution_count”: 3, “id”: “650cfce4”, “metadata”: {
- “execution”: {
“iopub.execute_input”: “2021-08-17T08:11:31.559585Z”, “iopub.status.busy”: “2021-08-17T08:11:31.558481Z”, “iopub.status.idle”: “2021-08-17T08:11:31.566183Z”, “shell.execute_reply”: “2021-08-17T08:11:31.567026Z”
}, “papermill”: {
“duration”: 0.023997, “end_time”: “2021-08-17T08:11:31.567345”, “exception”: false, “start_time”: “2021-08-17T08:11:31.543348”, “status”: “completed”
}, “tags”: [
“injected-parameters”
]
}, “outputs”: [], “source”: [
“# Parametersn”, “func_name = "PhAbs"n”, “wide_energy_range = Falsen”, “x_scale = "log"n”, “y_scale = "log"n”, “linear_range = Falsen”
]
}, {
“cell_type”: “code”, “execution_count”: 4, “id”: “7da50fd5”, “metadata”: {
- “execution”: {
“iopub.execute_input”: “2021-08-17T08:11:31.598602Z”, “iopub.status.busy”: “2021-08-17T08:11:31.597361Z”, “iopub.status.idle”: “2021-08-17T08:11:31.614771Z”, “shell.execute_reply”: “2021-08-17T08:11:31.615620Z”
}, “lines_to_next_cell”: 0, “nbsphinx”: “hidden”, “papermill”: {
“duration”: 0.038698, “end_time”: “2021-08-17T08:11:31.616172”, “exception”: false, “start_time”: “2021-08-17T08:11:31.577474”, “status”: “completed”
}, “tags”: []
}, “outputs”: [], “source”: [
“func = _known_functions[func_name]()n”, “n”, “if wide_energy_range:n”, “n”, ” energy_grid = np.geomspace(1e2,1e4,500)n”, ” n”, “else:n”, ” n”, ” energy_grid = np.geomspace(2e-1,1e1,1000)n”, “n”, “if linear_range:n”, “n”, “tenergy_grid = np.linspace(-5,5,1000)n”, “n”, ” n”, “blue = "#4152E3"n”, “red = "#E3414B"n”, “green = "#41E39E"”
]
}, {
“cell_type”: “markdown”, “id”: “f9b2012d”, “metadata”: {
“lines_to_next_cell”: 0, “papermill”: {
“duration”: 0.009853, “end_time”: “2021-08-17T08:11:31.635962”, “exception”: false, “start_time”: “2021-08-17T08:11:31.626109”, “status”: “completed”
}, “tags”: []
}, “source”: [
“## Description”
]
}, {
“cell_type”: “code”, “execution_count”: 5, “id”: “7764065d”, “metadata”: {
- “execution”: {
“iopub.execute_input”: “2021-08-17T08:11:31.667261Z”, “iopub.status.busy”: “2021-08-17T08:11:31.666028Z”, “iopub.status.idle”: “2021-08-17T08:11:31.670720Z”, “shell.execute_reply”: “2021-08-17T08:11:31.671789Z”
}, “papermill”: {
“duration”: 0.026666, “end_time”: “2021-08-17T08:11:31.672142”, “exception”: false, “start_time”: “2021-08-17T08:11:31.645476”, “status”: “completed”
}, “tags”: []
}, “outputs”: [
- {
- “data”: {
- “text/html”: [
“<ul>n”, “n”, “<li>description: Photometric absorption (phabs implementation), f(E) = exp(- NH * sigma(E)) contributed by Dominique Eckert</li>n”, “n”, “<li>formula: $n.a.$</li>n”, “n”, “<li>parameters: n”, “<ul>n”, “n”, “<li>NH: n”, “<ul>n”, “n”, “<li>value: 1.0</li>n”, “n”, “<li>desc: absorbing column density in units of 1e22 particles per cm^2</li>n”, “n”, “<li>min_value: 0.0001</li>n”, “n”, “<li>max_value: 10000.0</li>n”, “n”, “<li>unit: </li>n”, “n”, “<li>is_normalization: False</li>n”, “n”, “<li>delta: 0.1</li>n”, “n”, “<li>free: True</li>n”, “n”, “</ul>n”, “n”, “</li>n”, “n”, “<li>redshift: n”, “<ul>n”, “n”, “<li>value: 0.0</li>n”, “n”, “<li>desc: the redshift of the source</li>n”, “n”, “<li>min_value: 0.0</li>n”, “n”, “<li>max_value: 15.0</li>n”, “n”, “<li>unit: </li>n”, “n”, “<li>is_normalization: False</li>n”, “n”, “<li>delta: 0.1</li>n”, “n”, “<li>free: False</li>n”, “n”, “</ul>n”, “n”, “</li>n”, “n”, “</ul>n”, “n”, “</li>n”, “n”, “</ul>n”
], “text/plain”: [
” * description: Photometric absorption (phabs implementation), f(E) = exp(- NH * sigma(E))n”, ” * contributed by Dominique Eckertn”, ” * formula: $n.a.$n”, ” * parameters:n”, ” * NH:n”, ” * value: 1.0n”, ” * desc: absorbing column density in units of 1e22 particles per cm^2n”, ” * min_value: 0.0001n”, ” * max_value: 10000.0n”, ” * unit: ‘’n”, ” * is_normalization: falsen”, ” * delta: 0.1n”, ” * free: truen”, ” * redshift:n”, ” * value: 0.0n”, ” * desc: the redshift of the sourcen”, ” * min_value: 0.0n”, ” * max_value: 15.0n”, ” * unit: ‘’n”, ” * is_normalization: falsen”, ” * delta: 0.1n”, ” * free: false”
]
}, “metadata”: {}, “output_type”: “display_data”
}
], “source”: [
“func.display()”
]
}, {
“cell_type”: “markdown”, “id”: “5ac53db5”, “metadata”: {
- “papermill”: {
“duration”: 0.010751, “end_time”: “2021-08-17T08:11:31.693614”, “exception”: false, “start_time”: “2021-08-17T08:11:31.682863”, “status”: “completed”
}, “tags”: []
}, “source”: [
“## Shape n”, “n”, “The shape of the function. n”, “n”, “If this is not a photon model but a prior or linear function then ignore the units as these docs are auto-generated”
]
}, {
“cell_type”: “code”, “execution_count”: 6, “id”: “b8d45292”, “metadata”: {
- “execution”: {
“iopub.execute_input”: “2021-08-17T08:11:31.748678Z”, “iopub.status.busy”: “2021-08-17T08:11:31.747373Z”, “iopub.status.idle”: “2021-08-17T08:11:34.238452Z”, “shell.execute_reply”: “2021-08-17T08:11:34.239985Z”
}, “papermill”: {
“duration”: 2.536191, “end_time”: “2021-08-17T08:11:34.240375”, “exception”: false, “start_time”: “2021-08-17T08:11:31.704184”, “status”: “completed”
}, “tags”: [
“nbsphinx-thumbnail”
]
}, “outputs”: [
- {
- “data”: {
“image/png”: “iVBORw0KGgoAAAANSUhEUgAAAZYAAAEWCAYAAABFSLFOAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjQuMywgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/MnkTPAAAACXBIWXMAAAsTAAALEwEAmpwYAAAoOklEQVR4nO3de3xcdZ3/8ddkkknSptdAgbZAKdIPhJai3NVlAYVVERVcl3JZwcuCi7KoyK4//Snr6q5S23UFF1RWQauACHK/qdhSC0UgSMWW/RRaW5q2FEhvSZrLXM7+cSbpNE3SSTKTycy8n49HHpnznXP5tk3nne/lfE8kCAJERERypaLQFRARkdKiYBERkZxSsIiISE4pWEREJKcULCIiklMKFhERyanKQlegkBobGzsI/w6aC10XEZEiUg8kjjvuuJq+3izrYCH880crKyunFLoiIiLFIpFIDPh+UQaLmVUCpwNnufs1wzhVc2Vl5ZS5c+fmqGYiIqVvxYoVJBKJfnt6Rt0Yi5nNMbMHM7avMrP1ZrbGzOali2cCpwGRQtRRRET6N6qCxcyuAx7K2G4ALgcagJOB+WY2wd1XAzcXppYiIjKQ0dYVtgxoIQwRgHOA2929DWgzsyeAM4B7BntiM6snHHDqsWjRomg0Gh1ejUVEZA+jqsXi7g8Qhku3GcD6jO0mYOoQT38l4Jlfra2t9e3t7UM8nYiI9GW0tVh6SwHxXtsJAHdfB3xhEOe6Abgts6Curu6pWCxW38/+IiIyBKM9WJqA6Rnb04ClQzmRuzfT636VxsbGZCSi8X8RkVwaVV1hfbgPON/Mqs1sGnASQwwWEREZGaO6xeLuq8zsDmAlYTfYZ91dgyIio0QQBCQSkEgGJJMQT4TbyfR2MhmQTEEq1f09LEtlvE72LksF6fd2H5dMBgQBe3ylAiAISPUqD7+CXvuR3i8gSEEAPd9TKYDwmpnlPd8zHobY33MRe5cP9PzEzPf2eD3Ajnudf6Br9XPOvuo0pT7KJy6a0H9lh2jUBYu7LwGWZGwvABYUqj4ixSAIAjo6A3a1B+xqT9HeEdCefr2rPaC9I0V7e/jJUlkZobMroKsroLMroLMz/T1d1tG593vxREAyGZBIQiKxZ5hI8Zp5aFV5BItIuUsmA3bsTLFtR5Kt21Ns255k+84ULS0pWtpStLSmaNuVoqU1oKUtRWtritZdKfaxykZORSJQGYWqmghVlVAZjRCNhq+jlREqoxCNRohGoaICohWR8Hs0QrSCnteZ7/X7fvp7RSS8bvgVIRIJy4hkvFcRIQJEKsK7pyvSnf0VFZH0d/Z4v3v/3fvtWd59vd5/9oG2+3uv93hutsf1ud3vxr72jezxcmxtfsaYFSwiIyCRCHhza5ItbyTZ8kaCLW+Gr7duS9LSmqKyEna2pNi2PcWOltSAXSl9qayESRMqGFNbQW1thNqaSPi6JsKY2gi16depdLdULBahprqCWCxCdSxCTXWk53V1LEJ1dfi9u6yqKgyPymiEysrww16kPwoWkRzp7AzY+FqCpk1xmjYnaNqcYOPmBK+9nqR5WzLdl79v1bEIkydWMGlilEkTwu8TJ1QwflwF48ZWUFdXQd2YCsbVVTCuLkLd2ApqqiN7/UYsUigKFpFBSiQCNmxKsHZ9nDXr4qxdH2fdhjhvNCf7bWlUxyIcsH+UA/aPMmW/yp7X9ZOj1FZHSKZg0oQokyaGLQuFhBQzBYvIAFKpgFc3JljlXax6uYvVa7pYvyFOvI/xjDG1EaZPrWT6QZVMn1rJwVOrmHZQJQfuH2XC+AqFhZQNBYuUvZ0tKcaOCQef2ztSrPIuXnypk1Wru3jp5S7adu3ZDKmsDGfTHD6jquf7YYdUMXmiwkMEFCxSZlKpgFfWxXnuhQ6efaGDl17uoqsrfK/BYvgrXXtNoT10eiUNs2IcNSvGUUfEOGRaFVVVChCR/ihYpOTF4wEv/LmTZc+08+Qz7Wzd3vco+irvIhaDOUdWc0xDjKOtmqNmxagbO9oXqBAZXRQsUhK67wDvbknE4wHP/LGDJU/t4unGjj26sw6dXskJx9ZwwrE1zGmI8fyfOlm3Ic7sI6s58ogYMbVGRIZFwSJFbev2JL95YheP/q6NVzcm+NRHJ7BpS4LFT7azsyVsmUQicLTFeMeJtbzzxBqmT63a4xxvP6GWt59QW4jqi5QkBYsUnSAIu7bueaSV5c917DEmctNPdvS8PtpivPvUMfzVSbVMnqQHuomMFAWLjIjOzoAHftPKho0JPv3xiUPqburqCnh82S7ufrCVtevDx/TEYnD6O8bw1tnV3PDj7UyaUMGZfz2GM08dy7SD9OMtUgj6nyd5FY8HPPx4Gz+7ayfN28KuqTP/egyzj6zO+hy72lPc+0grdz3Qyvad4TkOmhLl3LPreM/pY3sG1//m9DE960iJSOEoWCRv/vB8OzfesoMNm8K7CSsrIZGAjs7sFsJq2xUGyi/vb2Vnaxgoc4+O8eGzx3HK8TV7rVfVvdCgiBSWgkVybsOmODfesoM/PN8BwFsOq+Lyv5/Aw4+3sfjJdrq6Bg6WeDzg3kda+fndLT2BctLbavjoR8Zx1KzsWzoiUhgKFsmZZDLgzvtbuPUXO4nHw9V2P3HhBP7m9DFEoxF++/tdAHTF+w6WIAhY9kwHP/jpdja9Fo7In3xcDX//kfEcdURsxP4cIjI8ChbJiTXruvj2jdtYvSZOJALnnV3HpeeP3+PmwljGPSa9rV0f54YfbWPFyvA2eDu8iis+NpE5R6mFIlJsijJYzGws8Bxwgru3Fro+5SwIAu56sJWbf7aDRAIOmVbJNZ+exNG2dyDEYmGwdMV3l3V0pvjpnS388oEWkknYvz7KJy8az7v+aozGTESKVMGCxczmAN909/ent68CPg8kgC+7+x0DHP5VYG3+aykD2bEzyXXf28bTjR1EInDBueO45O/G9wRIb90tlu4xludWdPCf39/Ga68niUbhovPGcdHfjqOmWkuoiBSzggSLmV0HXAD8Kb3dAFwONABjgEYze8Tdd/Rx7KXAo8ABI1Zh2cv/vtLFV69r5s2tSSZPrOBLn53M2+bUDHhMd+C0tqW44UfbuefhsLE5+8gYn7t8EocdUjXQ4SJSJArVYlkGtAAnp7fPAW539zagzcyeAM4ws5nAmRnHXQycDmwDTgQ+A3wrmwuaWT1Qn1m2aNGiaDSqO7IH6zdPtLHgpm3E43DCsdV88Z8mM2nCvv8eY+nc+MmdO0mlwodffeqSCZxz1lh1e4mUkIIEi7s/YGaZwTIDWJ6xSxMw1d0XAgt7HX4JgJndCnxvEJe9Erg2s6C1tZXaWq0Rla1UKuB/fr6DO+4NWxrzPjSOT1w4Puvnn3d3haVSMOvwKr501WQOmaZWikipGS2D9ykg3mu7j2f07ebulw7yGjcAt2UW1NXVPRWLxer72V8yxOMB131vK79b1k5VFVxzxWTefeqYQZ3jmIZqDjogyhnvHMNHPzJezzQRKVGjJViagOkZ29OApbm8gLs3A82ZZY2NjUkt/7Fv7R0pvragmWf+2Mn4cRV880v1Q7pR8YiZMX5+40F5qKGIjCajJVjuA35qZtcD+wEnAf9Y2CoJhAPtX/z3N1nlXexfH2X+V/fj0OnqvhKR/o2KeZ3uvgq4A1gJLAaudvf2wtZKdrWn+Jevh6Fy8NRKrv/3/RUqIrJPBWuxuPsSYEnG9gJgQaHqI7stf66djs5wva6XXu7ikGmVfOfr+2c180tEZLR0hcko8cjjbXz7xm0929MOrGTBvypURCR7o6IrTEaH5c+1s/D7u0PlwClRFn5tP/abrFARkeypxSIAvLy2i39buJVUCi77+wkcsH+UuQ3VeqSviAyagkXYuj3JV65rprMr4Lyz65j3oXGFrpKIFDF1hZW5rnjAtfObef3NJCccW80/XjKh0FUSkSKnYClz1//PdlZ6F9OnVvKVz9dnvTyLiEh/FCxl7NdL2nj4t22MHRPhG1+s3+OhXCIiQ6VPkjL16sY4/3XzdgCu/tQkLQYpIjmjYClDXV0BX//PrXR0BJxz1lhOe8fgFpMUERmIgqUM/fBnO1izLs7MQ6u44tKJha6OiJQYBUuZWbGyk1891EosBl/9/GSqqzVYLyK5pWApI+0dKeb/91YAPnHBBA7RgpIikgcKljJy8892sHlLktlHxjjv7LpCV0dESpSCpUy88OcO7n2kjepYhH/+9CTdryIieaNgKQPtHameFYs/cdF4pk9VF5iI5I+CpQzccvtONm9JcrTFOPe96gITkfxSsJS4v7wa51cPt1JVCddcoS4wEcm/olvd2MzOAy4CWoEF7v5igas0agVBwE23bieVgvM/OE6zwERkRBSsxWJmc8zswYztq8xsvZmtMbN5Axx6ATAP+BzwpXzXs5g93djBcys6mTSxgos+rKXwRWRkFCRYzOw64KGM7QbgcqABOBmYb2b9rd++ALgB+DygtUj6EY8H3PSTHQB88sIJjKlVr6eIjIxCdYUtA1oIQwTgHOB2d28D2szsCeAMM5sJnJlx3MXAbOBK4AjgsmwvaGb1QH1m2aJFi6LRaGk+IfH+x1pp2pTgLYdV8TenK39FZOQUJFjc/QEzywyWGcDyjF2agKnuvhBYmHmsmTUBtwC1hK2WbF0JXJtZ0NraSm1t7eAqXwR2taf42d0tAHz6YxOpqNCAvYiMnNEyeJ8C4r22E33t6O6PAY8N4Ro3ALdlFtTV1T0Vi8Xq+9m/aN39YCs7dqY48a01zD26utDVEZEyM1qCpQmYnrE9DViaywu4ezPQnFnW2NiYjERK67f5HS1J7rw/bK18/MLxBa6NiJSj0TKiex9wvplVm9k04CRyHCzl4hf3ttC2K+DUU2qZNTNW6OqISBkaFcHi7quAO4CVwGLgandvL2ytis+2HUnuebiNigr42Dy1VkSkMArWFebuS4AlGdsLCKcSyxD96sFWOrsC3n3qGA7VzZAiUiCjosUiw9e2K8W9j7YCcMG5uhlSRApHwVIi7n+slbZdAaccX8Nhh6i1IiKFo2ApAZ2dAXc9GLZWLjxPrRURKSwFSwl4dHEb27anmHt0jKNN962ISGEpWIpcKhXwywfC+1YuPFczwUSk8BQsRe6ZP3aw6bUkMw6u5Phj1VoRkcJTsBS5ex4Ox1bOfV8dpbaKgIgUJwVLEXu1Kc6zL3RSNzbCu0/VCsYiMjooWIrYPY+ErZX3vWsstTX6pxSR0UGfRkWmvT0FQGtbil8v2UUkAh98T12BayUispuCpYj8/O6dnH3xJl58qZNfL2mjvSPg7cfXcNABo2WRahGR0bNsvuxDMhnwo9t2AvC7Zbt4cVUXAB9Qa0VERhm1WIrE8y929rxu3pZk7atxpuwX5bhjNMVYREYXBUuRePR3bT2vn3q2A4D3nDFGjx0WkVFHwVIEWttSPPns7sfTpFIQicB7Th9bwFqJiPRNwVIEfv+Hdrq69ix72zHVHDhFQ2QiMvrsM1jM7LI+ymJmdm1+qrTXtSrN7Ewz+3ZG2Vgze8nMymLk+rdLdwFwxjtre8re9y61VkRkdMqmxXKVmT1pZrMBzOxMwkcIf2gwFzKzOWb2YMb2VWa23szWmNm8AQ6dCZwGZA4mfBVYO5jrF6s3mhO88OdOxtRGOP2d4d314+sqeMeJtfs4UkSkMLIJlrnA3cBiM1tK+Gz67wLHZXsRM7sOeChjuwG4HGgATgbmm9mEvo5199XAzRnHXgo8CryR7fWL2eO/bycI4NRTapnbUI0dXsXHLhhPrEqD9iIyOmXTSZ8EOghbDNOB14A/u3tqENdZBrQQhgjAOcDt7t4GtJnZE8AZZjYTODPjuIvd/c1e5zod2AacCHwG+FY2FTCzeqA+s2zRokXRaDQ6iD/GyOvuBnv3qWOoG1vBTfMPKHCNREQGlk2wvACMAS5091+b2T8Ad5nZUnc/L5uLuPsDZpYZLDOA5Rm7NAFT3X0hsHAf57oEwMxuBb6XzfXTrgT2GBdqbW2ltnb0din95dU4a9fH2W9ylLkNul9FRIpDNl1hDwBz3P3XAO5+M3AUsHMY100B8V7bif52dvd17v6FXmWXunvrIK55A2CZX3V1dc2jOVieWB62Vk57ey3RqLq+RKQ4ZNNiuR4Yb2a9H0/4z8O4bhNht1q3acDSYZxvn9y9GWjOLGtsbEyO5meYLF0e3rty6imjN/xERHrLJlheAwJ2z8oK0t+302vMYhDuA35qZtcD+wEnAf84xHOVpFeb4qzbkKB+cgUNs2KFro6ISNb2GSzu3tNdZmYR4EDgk8CQP+3cfZWZ3UE4bTkFfNbd2/dxWFlZ+nS6tXKylm0RkeIyqFu33T0ANgNfN7M/Al8ZxLFLgCUZ2wuABYO5fjnZHSzqBhOR4rLPYDGzKb2KKoDjAX3i5cnGzQle+UucSRMqmH2kusFEpLgMdYxlA3BFvipV7pY+Hc4Ge+dJmg0mIsVnUGMsMjKefCZcFv+v1A0mIkWo32Axs31OJ3b3+bmtjmzbkeSll7sYOyaimyJFpCgN1GI5asRqIT2eeb6DIIDj59ZQpfXARKQIDRQsS939FjP7uLv/eMRqVOaWPxd2g51yfE2BayIiMjQDBct3zCwG3Jhe52uvX5/d/c681awMxeMBz77QQSQCJ71NwSIixWmgYPk6cD5hoPQ1AywAFCw5tGJVJ+0dAUdbjAnjR/eqyyIi/ek3WLpXGjazW9z9YyNYp7L1tLrBRKQE7HMqsUJlZARBwPLnwrvtTzle04xFpHjpHpVRYsOmBJtfT3LA/lFmHDyolXZEREYVBcso8dyKTgBOOLaG0byUv4jIvmT1q7GZVQH7997f3V/NR6XK0XMvhOMrx8/VTZEiUtyyWYTyMuC7QBRIZrwVED6yWIYpHg9YsbKTigp46xwN3ItIccumxXItMA940N2T+9pZBm/V6i7aOwKOOiLGuDr1TopIccv2U2yxQiV/nluhbjARKR3ZtFj+FfihmX2NXs+Md/fX81GpctP4p3SwHKtuMBEpftkEyw/S3/+uV3lAOO6SV2ZWCZwOnOXu15jZecBFQCuwwN1fzHcd8mlnSwp/Jc6Y2ghHHaGHeolI8Rux57GY2Rzgm+7+/vT2VcDngQTwZXe/o59DZwKnsXutsgsIx3zGAf+d3i5af/xzuJrxsbOrqazUNGMRKX7ZTjd+F/BhYAqwDvipu/8p24uY2XWEAfCn9HYDcDnQQDizrNHMHnH3Hb2PdffVZnYz8Jl00QLgBuBNSmBW2rMvhPevHHeMusFEpDRkO934K8BNwB+BI4Bfm9lX3P3mLK+zDGgBTk5vnwPc7u5tQJuZPQGcYWYzgTMzjrvY3d/sda7ZwJXpelyW5fUxs3qgPrNs0aJF0Wi0cIs9BkHQc//KCcdq4F5ESkM2LZZ/Ac5x9xe6C8zsTuBXQFbB4u4PpJfe7w6WGcDyjF2agKndC1/u43RNwC1ALWFXWrauJJw63aO1tZXa2sKty7VhY4LX30xy0JQo0w7SMi4iUhqy+TSbAKzvVfYKMHEY100B8V7bif52dvd1wBfSrx8DHhvCNW8AbsssqKureyoWi9X3s3/e9Uwz1jIuIlJCsgmW+4EfmNk/u/s6MzsI+Cbw+DCu2wRMz9ieBiwdxvn2yd2b6TVdurGxMVnID/Tu9cE0zVhESkk2M76uBLYDq8wsTth6GQt8chjXvQ8438yqzWwacBJ5DpbRJpkMePGlMFjmNmh8RURKRzbBcrS7X0YYJtOAMe7+EeCQoV7U3VcBdwArgcXA1e7ePtTzFaN1G+K07Qo47JBKxo/TMi4iUjr67QozsxPTL5eb2cnsvo9khpntRzh4n3UfjrsvAZZkbC8gnDpclv60qguAOUeptSIipWWgMZZfpL9H2PvZ9gnC6ccyRN3dYAoWESk1Az3z/jCA9I2L7x25KpW+INg9vnJMg5ZxEZHSks2SLu81syOBD7D7zvtfuPsbea5bydq0JUnzthQHTomyf73uXxGR0rLPUWMz+xDhnfNvATqBU4H/NbMP5LdqpevFVeoGE5HSlc2vy/8OfMTdF3cXmNlZwPcI73GRQdo9vqJuMBEpPdnMc50KPN+rbDnh1GMZghdfCmeEHaMWi4iUoGyCZTFwrZnVQM/zUb4IPJ3PipWqrduTNG1OMGF8BQdP0/iKiJSebILlcsIVhbea2RpgG+EKxP+Qz4qVqpUetlZmW0zrg4lIScpmVtgbwFlmdiDh3favpxeFlCFY5eH4SoNpfEVESlM2z2OpJ2y1zOje38wAcPeP57FuJWnV6rDF0jBL4ysiUpqy6eS/D+gCfgsk81ud0pZIBKxeE6eiAuwtVYWujohIXmQTLMcBB7j7znxXptStWRensyvgiJlV1FRr4UkRKU3ZfLr9Hjg+3xUpByt7usE0viIipWug1Y1vTL/cDjxkZo8Bm4Ggex93vyKvtSsxL60OB+6PNo2viEjpGqgrbEvG95UjUJeStyo91fgotVhEpIQNtLrx1zK3zWwm4SKUm919fb4rVmq2bk+y+fUkE8dXMPWAaKGrIyKSN9lMN55N+DyWw4DXgClmthr4O3d/Oc/1KxkvvZweX9GNkSJS4rIZvL8l/TU2/YyWycBjwE/yWbFuZnaqmd2Z/jotXTbWzF4ys7qRqEMurH4lDBZ7i7rBRKS0ZRMsBlzv7ikAd+8E/g2YO5gLmdkcM3swY/sqM1tvZmvMbN4Ah54KfBT4JHBJuuyrwNrBXL/QfE0cgFkzdf+KiJS2bILlZuBTvcouYu/HFffLzK4DHsrYbiC8m78BOBmYb2YT+jrW3b8BTATmAzeb2aXAo0DRPGgsCAJeXhu2WGYdrhaLiJS2bG6QfBcwx8z+P9AEHEzYHfaKma3q3sndGwY4xzKghTBEAM4Bbnf3NqDNzJ4AzkhPEDgz47iLCVtG5wBfdvdmM7uccCHME4HPAN/K4s/QvTRNfWbZokWLotFo/gfS39yaZNuOFFP2izJpggbuRaS0ZRMsVw33Iu7+gJllBssMwme6dGsCprr7QmBh5rFm9k+E99IsNLOn3P2SdPmthA8by9aVwLWZBa2trdTW1g7iFEPjr4TdYEeoG0xEykA2qxs/kYfrpoB4r+1EP9f/YD/llw7ymjcAt2UW1NXVPRWLxer72T9nVqe7wUzdYCJSBgr1pKkmYHrG9jRgaT4v6O7NQHNmWWNjY3Ikpv6uXqPxFREpH4VaCfE+4HwzqzazacBJ5DlYCiUIwhWNQTPCRKQ8FCRY3H0VcAfhUjGLgavdvb0Qdcm3N5qTbN8ZDtxP1MC9iJSBEesKc/clwJKM7QXAgpG6fqH03L9yuForIlIe9FCQPOseX9HAvYiUCwVLnvXcGDlTwSIi5UHBkmev/CXsCnuLBu5FpEwoWPJo244kW7enmDyxQnfci0jZULDk0dp1YWvl8BlqrYhI+VCw5NGa9QoWESk/CpY8WpNuscw8VAP3IlI+FCx5tEZdYSJShhQseRKPB7y6MU5VJRw8tVBLsomIjDwFS568ujFOIgGHHlxFZaWecS8i5UPBkifqBhORcqVgyRPNCBORcqVgyZOee1gOVbCISHlRsORBEAQZXWGaaiwi5UXBkgdbt6fYvjPF/vVRxo/TX7GIlBd96uXBX17tvjFS3WAiUn4ULHmwbkMYLDMO1v0rIlJ+Rv0nn5mdCnwmvXkjMBm4CGgFFrj7i4WqW3/WbUgAMONgtVhEpPyMWLCY2Rzgm+7+/vT2VcDngQTwZXe/o59DTwU+CsSA7wJ1wDxgHPDfwAV5rvqgdbdYDlWwiEgZGpGuMDO7DngoY7sBuBxoAE4G5pvZhL6OdfdvABOB+cDNwALgBsJQGpPXig9BEASs7w6WaaO+QSgiknMj9cm3DGghDBGAc4Db3b0NaDOzJ4AzzGwmcGbGcRcDc9P7f9ndm83sE8CVwBHAZdlWwMzqgfrMskWLFkWj0dw+gOvNrSnadgUcOCVKba2GsESk/IxIsLj7A2aWGSwzgOUZuzQBU919IbAw81gz+ydgO7DQzJ4C1gO3ALWErZZsXQlcm1nQ2tpKbW3tIE6xb7sH7tUNJiLlqVB9NSkg3ms70deO7v7BPoofG8I1bwBuyyyoq6t7KhaL1fez/5BoRpiIlLtCffo1AdMztqcBS/N5QXdvBpozyxobG5ORSG5XHu4ZX5muFouIlKdCDQLcB5xvZtVmNg04iTwHy0jRVGMRKXcFCRZ3XwXcAawEFgNXu3t7IeqSS0EQsL4pbLEcMl1dYSJSnkbs08/dlwBLMrYXEE4dLhlvNCdp2xVw0JQotTWaESYi5Umffjm0Pt0NphsjRaScKVhySDPCREQULDnVPb6igXsRKWcKlhx6dWPYFXbwNAWLiJQvBUsObdiUDpap6goTkfKlYMmRHS1JduxMMXliBXVj9dcqIuVLn4A5sqGnG0ytFREpbwqWHNndDabxFREpbwqWHOlpsWh8RUTKnIIlRzZsCqcaqytMRMqdgiVHmjZpqrGICChYciKZDNj4WoKqSjhw/9w+kVJEpNgoWHJg85YEiQRMPbCSaDS3z3cRESk2CpYc6JkRpvEVEREFSy5oqrGIyG4KlhzQVGMRkd0ULDnQPdVYT40UERnBJ0gOlZm9C7gMGAP8h7svN7OxwHPACe7eWtAKktliUVeYiMiItVjMbI6ZPZixfZWZrTezNWY2b4BDZwMXAdcC706XfRVYm7/aZq+1LcW2HSkmjq9gXJ0agCIiI/JJaGbXAQ9lbDcAlwMNwMnAfDOb0Nex7v5d4B3A/wDLzOxS4FHgjTxXOyubXgtbK9MOGvWNPxGRETFSn4bLgBbCEAE4B7jd3duANjN7AjjDzGYCZ2YcdzHwAXf/sZmdDPwE6AK2AScCnwG+lU0FzKweqM8sW7RoUTQaHd4NjRsVLCIiexiRT0N3f8DMMoNlBrA8Y5cmYKq7LwQWZh5rZhEz+xlQA/zU3e9Pl98KfG8Q1biSsDutR2trK7W1tYM4xd42bg6DZeoBChYRESjc4H0KiPfaTvS1o7v/CPhRH+WXDvKaNwC3ZRbU1dU9FYvF6vvZPyubtqjFIiKSqVCfhk3A9IztacDSfF7Q3ZuB5syyxsbGZCQyvCVYelosBypYRESgcPex3Aecb2bVZjYNOIk8B0u+9AzeK1hERIACBYu7rwLuAFYCi4Gr3b29EHUZjvaOFM3bUoyrizB+nKYai4jACHaFufsSYEnG9gJgwUhdPx82b0kCGrgXEcmkX7OHoXuqscZXRER2U7AMQ/fAvWaEiYjspmAZBg3ci4jsTcEyDLrrXkRkbwqWYehusWjwXkRkNwXLEHXFA15/M0lNTYRJE/XXKCLSTZ+IQ/TalgRBEI6vDPfufRGRUqJgGaKNGrgXEemTgmWINukeFhGRPilYhqilNUVFBUw9cHjPcxERKTX6dXuILp03gYv/djypVKFrIiIyuihYhqGyUoP2IiK9qStMRERySsEiIiI5pWAREZGcUrCIiEhOKVhERCSnFCwiIpJT5T7duD6RSLBixYpC10NEpGgkEgmA+v7eL/dgSQAkEonmoRzc1dVVHYvFOnNbpcJfP1fnHep5Bntctvtns9++9kmlUpGurq4xsVhsV0VFRZBtHYtBIX+e83ntXJy7GH+W97XfMH+W60l/fvYpCAJ9DfFr1qxZt5Xi9XN13qGeZ7DHZbt/Nvvta59ZoWDWrFmzCvlvP5r/3UfbtXNx7mL8Wd7Xfvn8WdYYi4iI5JSCZXjuLtHr5+q8Qz3PYI/Ldv9s9iv0v2khFfLPns9r5+LcxfizPJTr50QkCEqqm1gkr8xsFuCAufvqQtdHZKjy+bOsFouIiOSUgkVkcJqBr6W/ixSzvP0sqytMRERySi0WERHJKQWLiIjklIJFRERySsEiIiI5pWAREZGcUrCIiEhOKVhERCSnFCwiIpJTChYREcmpcn/Ql8iwmNkc4BtAG/Add3+2wFUSGTIzqwROB85y92uGeh61WET6YWZzzOzBjO2rzGy9ma0xs3np4k8BH09/fbYA1RTJSpY/zzOB04DIcK6lYBHpg5ldBzyUsd0AXA40ACcD881sAlDv7s3u3gHEClJZkX3I9uc5vXz+zcO9nrrCRPq2DGgh/E8HcA5wu7u3AW1m9gRwBtBiZpOALqC1IDUV2bdsf57vycXF1GIR6YO7P0D4n7HbDGB9xnYTMBX4L+CH6a8FI1Q9kUEZxM9zTqjFIpKdFBDvtZ1w95XARwpTJZEh6/PnGcDd1wFfGM7J1WIRyU4TMD1jexqwrjBVERm2vP48K1hEsnMfcL6ZVZvZNOAkYGmB6yQyVHn9eVawiGTB3VcBdwArgcXA1e7eXthaiQxNvn+e9WhiERHJKbVYREQkpxQsIiKSUwoWERHJKQWLiIjklIJFRERySsEiIiI5pSVdRIqUmV0FvAZsAb7v7kcO4thlwOPufm1GWQR4FfhX4B3AN9x9bU4rLWVBLRaREWRm0RydZzzwUeCuIZ7iduC8XmUnAlOAu4HrgflDrqCUNbVYRPpgZu8CvgMcTvgci0+5+1YzuxV4nfA3+rnAr4EL3b3DzGYQPsviFMI7mi9z9xVmdinwYSAg/GXu/Wb2OeBfgHbg58A7gQ8AbwCz3f3ldD2eBH7k7j/uVcV/AO5196SZZdZ7MvB74Ifu/l0zOxb4PjAHeAr4pLuvB34JfNfMjui+FmHQPOLu24EXzOxQMzs6vdCmSNbUYhHpxcwOBm4DrgAOJXzWyrczdrkk/d5Mwg/seWZWAdwP/Ipw+fGfE/7m3+19wE+Bc83sdOAqwjA5DfgQgLvvBH6X3hczqweOT5+3tw+n982s9xjCELw3HSp16e1vES4y+Dzw4/S1XgceT5+n27mELZlui4EP9vPXJNIvBYvI3i4C7nL3Ze7+JvBlwg/dbovcfUX6w/lJwlViTwSi7n6Tu+909+uBqJkdkz7mOXe/y93jwIWELYpX0q2H/8w49z3A2enX7wGeTNehRzrETiBsFXWrJOwW2+7uX06XvR94wd3vTbdCvgKcYmYT0+/3dIeZ2RzgIPYMsRWErS+RQVFXmMjeDgU+ZWZXZBaaWU36ZeYHfQfh/6NDgQYz6734XvfDk7ZllB0CPJGx3ZTx+j7gv8xsLGHL5Vd91K+e8JfCHRllhwPNwFvMbLK7b03X6X191OlAYDthiN2UbqGdB9zXayHCrYQtHZFBUYtFZG9bgG+7e8TdI0AN8Nb0c+0HOubZ7mPSx81lzwDp1kY4SN7tkO4X7r4F+CNwVvrrvj6OrwQi6a9uTYRda08DX8+o0y8z6lNB2LW2On2tHcAjhKFyHnt2g0E4JpQa4M8s0ie1WET2dhfwsJktIpx++x+ELY9zBzjmaWCKmX0I+A3hB/V1wGF97Ps48Bkzu4vw/+AXCKcNd7sHuBZY6+4b+jj+DcJB/4mErQqANnePm9k1wPNmdhNhaHzTzN4JvABcCXzE3d+Wca7b03++SYQTETJNJJyoIDIoarGI9OLufyacsXUvsIl019g+jukgHIT/f4Qf/J8DznP3zj52/wHwGOEYxuL060TG+/cQtnbu7OdaCeBZ4Jg+3lsF3AJ8N936uZhwVtjrwHuBeb0OeZCwa6x7/CfTbMLAFBkUPY9FZISZ2QHQ0+2FmV0OnO7u89LblYTjJQ3uvrGfc1wDjHf3r+Sxnk8STrN+MV/XkNKkFovIyHsv8JCZHWBmRxBOPf4NgJnVAhcQzgbrM1TSbgbOTodQzpnZbGCbQkWGQsEiMvJ+RjhAvxr4A2FX2K3p964Hvkk4Nbhf6enDP2Tvrq1c+Sxhd6DIoKkrTEREckotFhERySkFi4iI5JSCRUREckrBIiIiOaVgERGRnFKwiIhITv0foQLMXn7HzsEAAAAASUVORK5CYII=n”, “text/plain”: [
“<Figure size 432x288 with 1 Axes>”
]
}, “metadata”: {
“needs_background”: “light”
}, “output_type”: “display_data”
}
], “source”: [
“fig, ax = plt.subplots()n”, “n”, “n”, “ax.plot(energy_grid, func(energy_grid), color=blue)n”, “n”, “ax.set_xlabel("energy (keV)")n”, “ax.set_ylabel("photon flux")n”, “ax.set_xscale(x_scale)n”, “ax.set_yscale(y_scale)n”
]
}, {
“cell_type”: “markdown”, “id”: “f9136f8a”, “metadata”: {
“lines_to_next_cell”: 0, “papermill”: {
“duration”: 0.011707, “end_time”: “2021-08-17T08:11:34.263975”, “exception”: false, “start_time”: “2021-08-17T08:11:34.252268”, “status”: “completed”
}, “tags”: []
}, “source”: [
“## F$_{\nu}$n”, “n”, “The F$_{\nu}$ shape of the photon modeln”, “if this is not a photon model, please ignore this auto-generated plot”
]
}, {
“cell_type”: “code”, “execution_count”: 7, “id”: “ce25ec18”, “metadata”: {
- “execution”: {
“iopub.execute_input”: “2021-08-17T08:11:34.332487Z”, “iopub.status.busy”: “2021-08-17T08:11:34.325165Z”, “iopub.status.idle”: “2021-08-17T08:11:35.936644Z”, “shell.execute_reply”: “2021-08-17T08:11:35.935946Z”
}, “papermill”: {
“duration”: 1.661859, “end_time”: “2021-08-17T08:11:35.936993”, “exception”: false, “start_time”: “2021-08-17T08:11:34.275134”, “status”: “completed”
}, “tags”: []
}, “outputs”: [
- {
- “data”: {
“image/png”: “iVBORw0KGgoAAAANSUhEUgAAAZoAAAEUCAYAAAASvPDLAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjQuMywgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/MnkTPAAAACXBIWXMAAAsTAAALEwEAmpwYAAAqvUlEQVR4nO3de3xcdZ3/8ddMJpOkDb0FC7SllAL9QGlhVaBVELks6Aqo4CIgeIFVYF3ZIsjPVX67yMqPCrZy2xWXmytFqCsCFYEVuRQs5VKCFGjxUwq0EC6lpPSS+1zO748zCdOQ0kwykzNJ3s/HI4/M9+RcPoHpvPM953u+JxYEASIiIqUSj7oAEREZ2hQ0IiJSUgoaEREpKQWNiIiUVCLqAspJfX19G+F/k8aoaxERGUTqgPTHP/7x6p5+qKDZWgKoSCQS46MuRERksEin0x/6cwXN1hoTicT4/fffP+o6REQGjeXLl5NOp7d5JkjXaEREpKQUNCIiUlIKGhERKSkFjYiIlJSCRkRESkpBIyIiJaWgERGRktJ9NCIig0yQzUIqRZBKEXSkCDo6wu+pDujILUvllufWoet1B0Eq3fUz8tapqBvH6G9+vej1KmhERAoQZDLhh3R7x1bf2daHeI8f+lu3e/rQ79yGbusEHSnYzp34fZXYfTcFjYhIpyAIoKODbFt77gO/PfwAb2/PfTDnvre358Kg27Ku7+8HBp2vOwOkPa/duU4mE+0vnkgQq6khlqwklkxCZWX4ujJsx3JtOtvJJLFkglhlEpJbr9O1TW79+KhRpSm5JHsdQGY2E7gEaAaucPdlEZckIuQFQWsbQWsrQWsb2dZWgrb23Fdb+EHf1k427/X7y9vCn3Utb99qedDeDgP9hODKBLHqkeGHdFVV14d914d1/od4Tx/63T7g6XFZ9xDotp/44Lu0XtZBkwuRue5+bK49BzgPSAMXuvtC4GzgDMKguRE4NaJyRYaEIJslaGkl29RMtrmJoKk5fN3UTLYp125uDtfpDJG2NrItYZgEre8vJ5stWZ2x6ipi1dXhB3F1Vd4HfpJYVe4r1yaZ165KEqvs1k7mQmOrZd22rawkVlFRst9nKCvboDGzy4BTgOdy7enAWcB0YARQb2b3AXXu3phbJ1nA/usIp7busmDBgooKvZFkCAjSabIbN5HZtInspi1kN2/O+76Z7ObwdZBOQzboCpVsUzNBc0txAiKZJF5TTaymhviImjAUaqrD79VVxKuqul7HqqvDHkLudbzrdd46uWXx6urwwz8W63+NMiDKNmiAJcAWYHaufRxwm7s3A81m9ghwBLDFzMYCHUBTAfs/B7gof0FTUxM1NTX9LlykVLItLWTeeZfM+nfJvPce2Q3vkXlvI5kN75HN+57dtLnPx4hVVxEbOZJ4bS3x2pHEa0cSq+1sh8tiI0cQr6kJg2NEzfuva/Je6482ySnboHH3u80sP2imAI/nrdIATACuBK4jDJpLCzjENcCt+Qtqa2uXJpPJum2sL1IyQUeK9NvrwgBpbIQAMo2NZNatD5e9s570O+sJmpp7tb9YdTXxcWOpGDOK+KhRxEd3ft+B+KhRVORek0gQi8eJ5YdKZWWJf1sZbso2aHqQBVLd2ml3XwGcWOjOcqfbtnp+Qn19fUbdcSmVIJMh8/Y6Ug1vkm54g3TDm6TfeJN0w5tk3lnfq9NV8TGjqfjIjlSM/wgV48YSHzuGinFjqRg7hvjYsV3L4jU9PuhQJBKDKWgagEl57YnAoxHVIrJNQTZL5q11pNasJfXqWlJrXiO9Zi2p1xq2ef9DrLqKxIRdqBj/EeKjRxOriIdhMv4jVIwPgyXxkR2JVVUN8G8j0n+DKWgWATeb2dXAjsAs4B+jLUmGuyAIyLz5Fh2rVtOxajUpX03HSy8TNPdwiqsyQWK3ySQmTSAxcUL4fdIEKidNJF43The3ZcgaNEHj7ivNbCGwgvC02bnu3hpxWTLMBB0ddPx1Fe3Pr6Tj+RW0r/xrj9dNEpMmUDl1Cokpu1G5e/iVmDhBF8hlWCrroHH3xcDivPY8YF5U9cjwE7S30/7Ci7T/ZTntz62gw1dBauvTX4lJE6icthfJaXtQOW1PknvuQbx2ZEQVi5Sfsg4akYGQbW2j/ZnltD3xFO0vvEj6jTeJJRIkpxvtz6+EVN4YlESC5L77UDVzX5L77UvVvvsQ36E2uuJFBgEFjQxL2c1baF36JK2PPkZb/bNbhwnhDY/tzywPg2W/GVR/bH+q9p9Bcu9puiAvUiAFjQwb2S1NtP55KS2PLAlDpHNyxESCqo//DdWzDqRqv33p8JdIv/U2VfvPpGq/fYnrJl6RflHQyJCT2bSJlvsfIvXyq4w+83RSq1+h+X8foPWxx9+/vpJMUv2Jgxjx6YOpnn0Q8ZEjurZPTtszospFhiYFjQwJQRDQvvx5mu++j9Yl7wdKy/0Pvb9SIkH1IZ9gxOGfonr2geqpiAwQBY0MakFHipaHH6Xp9rtIvfxquLAyQXzsGLLvbQybe+7ByM/+LTVHHkrF6NHRFSsyTCloJBLZpmaafreI1OtvMO6Cfy74Anu2qYmmu+6h6a4/kN3wHgAVO+9E7RePYcRnjqRi9Gjalj1DfNwYkntMLcWvICK9pKCRARV0dNB0x+/ZfNvtBFvCybY7vvA5qmbu26vts01NNN1xN1tuv6vrRsnkjOnU/v0XqDl49lY3RFYf+LHi/wIiUjAFjQyIIAho+/PjbPyvG8m8tQ4I5/fqfHLi9mRbW2n67V1bBUz1wbMZ9ZUTSe5jJa1dRPpHQSMll1rzGhuvvpb2Z58HILmPMfrb36Tpzj/Q+tAjBB2pbW4bZDK03P8Qm25aQLZxA5ALmK+dQnKvPQakfhHpHwWNlEyQSrHlttvZfMtvIJ2mYsc6Rp/5DWqO+DSxeJzme+8P1+vouUfT9pfn2HTtDaRWvwJA1Uf3Y/RZZ2j4scggo6CRkujwl9hw+ZWkX10L8Ti1J3+JUV89ZavnpMSS4ZO3u/doMu9tZNO1N9DywGIAErtOYvTZZ1A9+0DNcCwyCClopKiCbJYtv7mDzTctgEyGyj12Z+wFc3rshXQFTXt717bN997Pput+SdDUTGzkCEadfhq1n/8csYTeqiKDlf71StFkGjew4Sc/o73+WYjH2eFrpzDqtJO2GRL5PZr02+/w3uVXdF3HqTnsU4z5p29RUTduoMoXkRJR0EhRtC9/gcaL55LduImK8R9h3A/Pp2q/GR+6TSwZPpu+bekTbP7VrwmaW6jYaTxjzv02NbMOGIiyRWQAKGik35oW3cvG//gvyGSoPng24y6YQ3zUDtvdLlYV9mg6ezEjPnMkY75z1lbzjonI4Dfog8bMZgA/BjYBv3f3OyIuadgI0mk2XvMLmu/+XwBGnfFVdjj1y72+YB8fNarr+5jz/okRhx5cslpFJDplEzRmNhOY6+7H5tpzgPOANHChuy/cxqYHEj7aGeD5khcqQPiwsA0Xz6XtqXpiI2oY94PvUXPwrIL2MeLwT0EsRvUBH9W1GJEhrCyCxswuA04Bnsu1pwNnAdOBEUC9md3n7pt62HwZsAiIAf8JnNzLY9YBdfnLFixYUFGhZ7pvV2bTZhovvJiOlU7FjnXseNnFVO4+peD9xKqqGPmZI4tfoIiUlbIIGmAJsAWYnWsfB9zm7s1As5k9AhxhZlOBo/K2Ow34KLAq187Se+cAF+UvaGpqokZTx3+ozLuNrL/g/5Je+zqJyZPY8bJ/J7HT+KjLEpEyVhZB4+53m1l+0EwBHs9bpQGY4O7zgfn525rZW8ACoBW4uoDDXgPcmr+gtrZ2aTKZrNvG+sNepnED68/7IemGN0juY9RdehEVo0dFXZaIlLmyCJoeZIFUt3a6pxXd/QHggUIP4O6NQGP+svr6+ozuPN9aEAS03PcnMhs30nL/Q2HIzJjOjj/5EfERGh0mIttXrkHTAEzKa08EHo2olmFty22/ZfMNN3e1k9ONHecqZESk9+JRF7ANi4CTzKzKzCYCs1DQDLimP/zvViFTufc0dvzJv+s+FxEpSFn2aNx9pZktBFYQnjY7191bIy5rWGl7qp6NV/4cgHEXfo/E5F2p3G3XrmljRER6q2yCxt0XA4vz2vOAeVHVM5yl1rxG448vg2yW0WefwYgjD4u6JBEZxMr11JlEJLNpE+9eeDFBcwsjP3c0tSceH3VJIjLIKWikS5DJ0PijuWTeWkfV38xkzJx/1PNfRKTfFDTSZfMvb6Fj+QtU7LIz4y76AbHKyqhLEpEhQEEjALQ+sYwtt/4WKiupu+hfdCOmiBSNgkZIv/0OG+aGEy6M+c6ZPT4NU0SkrxQ0w1yQybBh7nyCLU3UHPlpRh772ahLEpEhRkEzzDX9bhEdz6+gYpedGPvdf9LFfxEpOgXNMJZ6dQ2bbrwZYjHGff+7mlZGREpCQTNMBakUG+ZeAak0tV8+nqr9ZkRdkogMUQqaYWrzr/+H1OqXSUyZzOjTT4u6HBEZwhQ0w1DHqtVsueU3UFHBuH85T/OXiUhJKWiGmSCT4b3Lr4Rslh1O/bKGMotIySlohpnm399L6pU1JKZMZtRpJ0VdjogMAwqaYSSzcRObfnkLAGPOOZtYomwm7xaRIUxBM4xsvvFmgqZmaj59CNUf3S/qckRkmBh0f9KaWQI4HDja3S8ws5nAJUAzcIW7L4u0wDLVsWo1zffeT6yqitFnnxF1OSIyjETWozGzmWb2h7z2HDNba2Yvm9nJH7LpVOAwoPMW9rOBM3Jf55ao3EEtCAI2/vwGCAJ2+Mrfk9hpfNQlicgwEknQmNllwD157enAWcB0YDZwuZmN7mlbd18FXJ+3qM7dG929DdA43R60Pfk0Hc+9QMWOdezw5ROiLkdEhpmoTp0tAbYQhgrAccBt7t4MNJvZI8ARZjYVOCpvu9Pc/d1u+9piZmOBDqCptwWYWR1Ql79swYIFFRUVFYX9JmUuyGTYfMOvABj1jVOJVVVFXJGIDDeRBI27321m+UEzBXg8b5UGYIK7zwfmb2d3VwLXEQbNpQWUcQ5wUf6CpqYmampqCthF+Wt5cHE4nHm3XRnxmSOjLkdEhqFyGQyQBVLd2ultrezua4Dv5V6vAE7swzGvAW7NX1BbW7s0mUzWbWP9QSfoSLH5pnA48+hvfp3YEOuticjgUC5B0wBMymtPBB4t5QHdvRFozF9WX1+fGUrT5Dff/yCZd9aTnG5Uf3JW1OWIyDDVp8EAZlZrZsWcU34RcJKZVZnZRGAWJQ6aoS5Ip8NHMwM7fPUUPWdGRCLTqx6NmVUBpwPHA58ARuaWtwJPAHcCN7l7a1+KcPeVZrYQWEF42uzcvu5LQi0PLCbz9joqp+1J9UEfj7ocERnGths0ZnYqMAdYTHjh/XRgAxAAY4AZwBHAEjP7hbtf3+OOunH3xbl9drbnAfMKqF22Ichk2HLr/wAw6rST1JsRkUj1pkdTAxzi7h09/Gxd7utBM7sI+Eoxi5O+aV38Z9INb1I5dYquzYhI5LYbNO5+Q0/LzexTwF7AEndf5e5p4OYi1ycFCoKALbfdDsAOp51ELK7p7EQkWv35FDqEcH6xi3K9GSkD7c8sJ/XKGip22ZmaT30y6nJERPoeNO4+191/4+6nAm+Z2fZurJQB0HT7XQDUfunzum9GRMpCn++jMbMbgF2AnYFKoK1YRUnfpF57nbYnnyY2ciQjP/u3UZcjIgL074bNnwNvAevcPVukeqQfmn63CICRx3yG+Ihi3uYkItJ32z11ZmavdWt/HsDdn3H3txQy0cq2tACQ2bSZlvsfhnic2uOPjbgqEZH39eYazaRu7TtLUYgUbvOChbx57Jdpf+4Fmu/5I0F7OzWHHqznzYhIWenNqbOgW1t3/5WBIJNh8y/DCTNbHv4zbU89DaDejIiUnb6MOusePBKBtief7nqdXvs6mbfWkdhtMskZ0yOsSkTkg3o719lHyOvJdG+7+zvFL00+TPM9f+x63f7scwCMPO6zmm5GRMpOb4ImBrzN1qfM1uW9DgDdsDGAMhve26pHA0AyycijjoimIBGRD9GbKWg0h0mZaXlwMWSzEI+H34ERhx1CfIfaaAsTEelBb4Y39+rplWYWNzNNqjkAWv74IAC1Xzima9nIYz8bVTkiIh+qN6fOdjGzZcBdwEPACnffDGBmY4D9gL8FPgdcV5oypVPH6lfCucx2rKPm0wfTdOfd4SCAffeJujQRkR715tTZ1Wb2S+BbhM+LOdDMYoTXbFqBZcC9wJHuvqmUxcr7vZkRRx9BcsZ0Rn3z61TPOkCDAESkbPVq1Jm7bwF+BvwsFzJ1uW3Xu3umhPV9gJklgMOBo939AjObAfwY2AT83t3vGMh6BlKQTofXZ4ARRx9JLB5n1Fd6dWZTRCQyBc915u4B8G5/D2xmM4G57n5srj0HOA9IAxe6+8JtbDoVOIz3R8EdSPj4Z4Dn+1tXOWt7+i9kN24iuY9RObn7hA0iIuWpP5Nq9pmZXQacAjyXa08HzgKmAyOAejO7r6dTce6+ysyuB76TW7QMWEQYPP8JnNzLGuoIe2ZdFixYUFFRxlPrty7+MwA1RxwacSUiIr0XSdAAS4AtwOxc+zjgNndvBprN7BHgCDObChyVt91p7t69N/VRYFXudSETfJ4DbPXAtqamJmpqagrYxcAJOlK0PvYExGKM+PQhUZcjItJrkQSNu99tZvlBMwV4PG+VBmCCu88HtvdAtbeABYQDE64uoIxrgFvzF9TW1i5NJpN121g/Um3LniFobiG53wwqdizLEkVEetTroDGzM939um7LksAP3P3iftaRBVLd2ultrezua4Dv5V4/ADxQ6AHdvRFozF9WX1+fKdfRWy2LHwXCGzNFRAaTQu76n2Nmj+VGeWFmRwErgC8WoY4Gtn4cwURgTRH2OyQE7e20LX0K4nFqDj046nJERApSSNDsD/wOeNjMHgUWAlcBHy9CHYuAk8ysyswmArOAR4uw3yGh9cmnCVpbqdp/BhXjxkZdjohIQQoJmgzQRji6axLhRJsvFOMJm+6+kjC4VgAPA+e7e2t/9ztUtD6yBICawz4VcSUiIoUrZDDAs4RDj0919z+a2beA283sUXc/odADu/tiYHFeex7hzAOSJ+hIhTM1x2LUHDx7+xuIiJSZQno0dwMz3P2PAO5+PbA34R35UiLty58jaGklOX1vnTYTkUGpkB7Nl4ATzKxUtUgPWh97EoCag2dFXImISN8UEjRn572OAzsBXweeKmpF0iUIAtqWhkFTffAnIq5GRKRveh007v5I92VmdjuwlG532EtxpFatJvNuI4nJk6jcdWLU5YiI9El/n545GhhXjELkg1qXPAFAzSc1CEBEBq9CZgZ4EQjyFsWB3YD/V+yiJNS6NAyaal2fEZFBrK/XaCCcJqbB3V8tYj2Sk357HelX1xIfO4bkPhqAISKDV7+u0UjptD1VD0D1QQcQi/f3DKeISHS2GzQ9nDL7AHefXrSKBMgPmo9FXImISP/0pkfzU+BlYE9gdWnLEYAglaL9meUQj1P98Y9GXY6ISL/0Jmj+g/DhYg8BO/P+I5SlRNpfeJGgrY3k9L2Jj9oh6nJERPqlN0FzH/Bi7vW6vOUBYegEQPk+/3gQanvqaQCqDyrGxNgiItHabtC4+5cAzOxVd9+99CVJ+7JnAAWNiAwNvR7OpJAZGJn175J6ZQ3x0aOonLZn1OWIiPSbxs2WmbbO3syBH9OwZhEZEvRJVmY6r89UHajTZiIyNBQyM0BZMLNDge/kmj8HGoFLgGbgCndfFlVt/RWk07TVPwuxGNUH6v4ZERkaCpnrbDFwB3Cnu7/e3wOb2Uxgrrsfm2vPAc4D0sCF7r5wG5seCnwNSAJXAS3AGYRBcyNwan9ri0rHir8SNLdQaXtRMWZ01OWIiBRFIafObgEOB1aY2VNm9kPr41PQzOwy4J689nTgLGA6MBu43Mx6/KR190uAMcDlwPVAnbs3unsbYfj0toY6M5uW/5XJZCqC4EMnQSiptidzw5pnHRBZDSIixVbIXGc3ADeYWSXwaeAY4A9m1gL8Gljg7m/1cndLgC2EoQJwHHCbuzcDzWb2CHCEmU0Fjsrb7jRg/9z6F7p7o5mdbmZjgQ6gqbe/D3AO3Z6j09TURE1NTQG7KK62+r8A6LSZiAwpfblGUwdMJnxEwI7A88BU4Ckzu87df7y9Hbj73WaWHzRTgMfzVmkAJrj7fGB+/rZm9s/ARmC+mS0FrgSuIwyaSwv4Pa4Bbs1fUFtbuzSZTNYVsI+iyTY1kXr5VWI1NST3nhZFCSIiJVHINZpLgM8B+wAPAncBZ7n7+tzPdwJeAbYbND3IAqlu7XRPK7r7F3pYfGKhB3T3RsKBBF3q6+szsVg0M+y0P78SslmSM/YhVqGJFkRk6CikR7MHMBe4N3eKq7tG4Og+1tEATMprTwQe7eO+BqX25S8AULX/zIgrEREprkKCZjYwC7is2xiANPAu8CfgJ32sYxFws5ldTXg6bhbwj33c16DUvvx5AKr2nxFxJSIixVXIqLPrgDXAPwCfAc4EXgduB34EfIrwvpaCuftKYCGwAngYON/dW/uyr8Eo29JC6qWXiVVXkbS9oi5HRKSoCunRnAPMzF3bAHjJzF4A/uLuPzSzlYRBcXpvdubui4HFee15wLwC6hkyOl54Mbw+s+90YolBdw+tiMiHKqRH0wHs2m3ZuLx9VObWkQK1P6vTZiIydBXy5/O/AH8ys98BrwHjgZOAK8xsEuEpr2uLX+LQ1/6cgkZEhq5CHhOwEPgkYchMBjLAP7j7TwhvlPy6u/9bSaocwrKtbXT4amJVVSRN98+IyNDTqx6NmcWBZ4CD3P0DN0W6+0bgkeKWNjx0vPhXyGRIztyXWLIy6nJERIquVz0ad88CTxOONJMi6nghfEp2cub0iCsRESmNQq7R7AecbmYXAu/l/8Dd9SnZR+0vrASgasY+EVciIlIahQTNBSWrYpgKMhk6Vv4VYjGS++wddTkiIiVRyOzNXddgzCzp7hrK3E+pNWsJWlqpnDqFeO3IqMsRESmJQibVHAP8FDgBGG1mHyOcZv/szok1pTBd12f21WkzERm6Crlh8xdAFfCx3HYvAm8TPtVS+qDz+kxyhi5xicjQVUjQHAOc6+5rgcDdU8APCZ+6KX3Q2aPRQAARGcoKCZpXgO6Pbt4VWFe8coaPzPp3yax7h3jdOCp23inqckRESqaQUWfnAXea2W+BmJn9BDgZ+H5JKhvi2lf+FYCqffchqoetiYgMhEKmoHkQOJDwIWW/AFqBY9z9NyWqbUjr+OsqAJL7dO8kiogMLQXNSe/urwOXlaiWYSXlLwGQ3Fvzm4nI0FbI8OYDgMsJH7m81XbuPrXIdX1YHYcC38k1f074dM8fA5uA37v7HQNVS18FmQwdvhricSqn7RF1OSIiJVVIj+Zm4B7gYqCtvwc2s5nAXHc/NteeQ3gdKA1cmJstuieHAl8DksBVwKNANvez5/tb10BIv/4GQWsriSmTidfURF2OiEhJFRI0k4B/K8Yjls3sMuAU4LlcezpwFjAdGAHUm9l97r6p+7bufomZ7Uz4+Ojrgc3AIiAG/CfhAIWy1nV9RqfNRGQYKCRorgfOBq4ownGXAFuA2bn2ccBt7t4MNJvZI8ARZjYVOCpvu9OA/XPrX+jujWb2VWBV7udZesnM6oC6/GULFiyoqKio6MvvU5COzuszev6MiAwDhQTNIcABZvZ9YENuWYzw5s2Cbm1397vNLD9opgCP563SAExw9/nA/PxtzeyfgY3AfDNbSnh/zwLCUXBXF1DGOYRT6HRpamqiZgBOZb3fo9mr5McSEYlaIUHzf/JeV1OE6zR5skCqWzvd04ru/oUeFj/Qh2NeA9yav6C2tnZpMpms28b6RRF0pEi9/CpUJqicOqWUhxIRKQuFBM1zhJNqHg+MJpzz7CLC02n91UB4DajTRMKL/CXj7o1AY/6y+vr6TKlvnky98iqk01TuPY1YpZ6oKSJDXyFT0FxLONKr+6SaNxShjkXASWZWZWYTgVmUOGiiooEAIjLcFGNSzSP6W4S7rwQWAiuAh4HzizG6rRx1DQTQ9RkRGSYKOXXWOalm/kX7Pk+q6e6LgcV57XnAvL7sazDpWLUagOQ0BY2IDA+aVHMABe3tpNe+TqyqisSuE6MuR0RkQBQ6qeYBaFLNPku9uhayWSqnTiE2APfriIiUg0In1WxAk2r2WcfqVwCo3HPApoYTEYlcIYMBpJ9SnUGzlybSFJHhQ0EzgFKrXwYgqR6NiAwjCpoBEmQypF5ZEz4aYPfdoi5HRGTAKGgGSLrhTYK2dhKTJxGrqoq6HBGRAaOgGSDvnzbT9RkRGV4UNAOka8TZXro+IyLDi4JmgHSNOFOPRkSGGQXNAAiCgNRLYdAk99g94mpERAaWgmYAZNa/S3bzZip2Gk981A5RlyMiMqAUNAMg9coaACrVmxGRYUhBMwBSr64FoHL3KdEWIiISAQXNAEi/ugaAyqm6UVNEhh8FzQDo6tFMUdCIyPBT0OzN5cDMjgTOBEYAlwJNwCVAM3CFuy+LsLwPCNJpUq+9DomEnkEjIsNSZEFjZjOBue5+bK49h/DhamngQndfuI1NZwCnAvsRPl56Z+AMwqC5MfezspFueBNS6fAZNIlBl+siIv0WyakzM7sMuCevPR04C5gOzAYuN7PRPW3r7lcBBwM3AEuAOndvdPc2IFlADXVmNi3/K5PJVARB0PdfrAepNeFps4Qm0hSRYSqqP7GXAFsIQwXgOOA2d28Gms3sEeAIM5sKHJW33WnA5939JjObDfwK2GJmY4EOwtNovXUOcFH+gqamJmpqavr0C21L19BmBY2IDFORBI27321m+UEzBXg8b5UGYIK7zwfm529rZjEzuwWoBm4GXgauIwyaSwso4xrg1vwFtbW1S5PJZF0B+9guDW0WkeGuXC4aZIFUt3a6pxXd/UbCazH5Tiz0gO7eCDTmL6uvr8/EYrFCd/Wh3g8a9WhEZHgql+HNDcCkvPZEYE00pRRPtrWNzFtvExtRQ8VO46MuR0QkEuUSNIuAk8ysyswmArOARyOuqd/Sa1+DIKByym4Uu6ckIjJYlEXQuPtKYCGwAngYON/dW6Otqv9SnTMC6LSZiAxjkV2jcffFwOK89jxgXlT1lIIGAoiIlEmPZqhKrXkNgMSUyRFXIiISHQVNCaVfawCgcrddI65ERCQ6CpoSyba2knlnPbGRI4iPGxt1OSIikVHQlEj69TcAqJy8q0aciciwpqApkdTa1wFITNZpMxEZ3hQ0JZJ+LQyayt0mbWdNEZGhTUFTIp0DAdSjEZHhTkFTIqmuHo2CRkSGNwVNCQTpdPjAs8pKKnbeKepyREQipaApgfSbb0EmQ+WuE4lVVERdjohIpBQ0JZDuHHG2qwYCiIgoaEog9XpuIICuz4iIKGhKobNHUzlZPRoREQVNCaReU49GRKSTgqbIgiAI76GJxaicNDHqckREIhfZ82j6ysyOBM4ERgCXAluAHwObgN+7+x0Rlkf23UaC1lYqdtmJWFVVlKWIiJSFyILGzGYCc9392Fx7DnAekAYudPeF29h0BnAqsB9wDNAAZHM/e76kRfdC6o03AUioNyMiAkR06szMLgPuyWtPB84CpgOzgcvNbHRP27r7VcDBwA3AEmAZ8C3gfMKeTaTSDbmgmTgh4kpERMpDVD2aJYSnvGbn2scBt7l7M9BsZo8AR5jZVOCovO1OAz7v7jeZ2WzgV8C9wKrcz7P0kpnVAXX5yxYsWFBR0c8bLLuCZpKCRkQEIgoad7/bzPKDZgrweN4qDcAEd58PzM/f1sxiZnYLUA3cDLQAC4BW4OoCyjgHuCh/QVNTEzU1NQXs4oPSb6hHIyKSr1wGA2SBVLd2uqcV3f1G4MZuix/owzGvAW7NX1BbW7s0mUzWbWP9Xuns0VSqRyMiApRP0DQA+Xc3TgQeLeUB3b0RaMxfVl9fn+nP0zCDTIb0m29CRYUm0xQRySmX+2gWASeZWZWZTQRmUeKgKYXM+nchlSaxy86aTFNEJKcsgsbdVwILgRXAw8D57t4abVWF00AAEZEPiuzUmbsvBhbntecB86Kqpxg0EEBE5IPKokczVKRffwPQzZoiIvkUNEXU1aPRqTMRkS4KmiLSqTMRkQ9S0BRJkE6TfmsdVFZSMX7HqMsRESkbCpoiSb+9DjIZEhN2IRbXf1YRkU76RCwSDW0WEemZgqZINBBARKRnCpoiyW5pgnhcAwFERLopl7nOBr3R3ziVUad+GYIg6lJERMqKgqaIYpWVUZcgIlJ2dOpMRERKSkEjIiIlpaAREZGSUtCIiEhJKWhERKSkFDQiIlJSGt68tbp0Os3y5cujrkNEZNBIp9MAddv6uYJma2mAdDrd2JeNOzo6qpLJZHtxS4r++MXab1/3U+h2vV2/N+ttb51sNhvr6OgYkUwmW+Lx+JC6WzfK93Mpj12MfQ/G9/L21uvne7mO3Odnj4Ig0FeRvqZNm3brUDx+sfbb1/0Uul1v1+/NettbZ1oomDZt2rQo/9+X8//3cjt2MfY9GN/L21uvlO9lXaMREZGSUtAU1++G6PGLtd++7qfQ7Xq7fm/Wi/r/aZSi/N1Leexi7Hswvpf7cvyiiAXBkDqtLDKgzGwa4IC5+6qo6xHpq1K+l9WjERGRklLQiPRPI3Bx7rvIYFay97JOnYmISEmpRyMiIiWloBERkZJS0IiISEkpaEREpKQUNCIiUlIKGhERKSkFjYiIlJSCRkRESkpBIyIiJaWgERGRktITNkWKyMxmApcAzcAV7r4s4pJE+szMEsDhwNHufkFf96MejUgvmdlMM/tDXnuOma01s5fN7OTc4rOBM3Jf50ZQpkiv9PL9PBU4DIj151gKGpFeMLPLgHvy2tOBs4DpwGzgcjMbDdS5e6O7twHJSIoV2Y7evp9zz6W5vr/H06kzkd5ZAmwh/EcIcBxwm7s3A81m9ghwBLDFzMYCHUBTJJWKbF9v3893FuNg6tGI9IK73034j7PTFGBtXrsBmABcCVyX+5o3QOWJFKSA93NRqEcj0jdZINWtnXb3FcCJ0ZQk0mc9vp8B3H0N8L3+7Fw9GpG+aQAm5bUnAmuiKUWk30r6flbQiPTNIuAkM6sys4nALODRiGsS6auSvp8VNCJ94O4rgYXACuBh4Hx3b422KpG+KfX7ORYEQbH2JSIi8gHq0YiISEkpaEREpKQUNCIiUlIKGhERKSkFjYiIlJSCRkRESkpT0IgMEWY2B3gbWAf8wt33LmDbJcCD7n5R3rIY8BrwI+Bg4BJ3f6WoRcuwoB6NSITMrKJI+xkFfA24vY+7uA04oduyg4DxwO+Aq4HL+1ygDGvq0Yj0gpkdCVwB7EH4HI+z3X2Dmf038A7hX/z7A/cDX3H3NjObQvgsj08Q3nF9prsvN7NvAF8CAsI/9o41s+8C3wdagV8DhwCfB9YDM9z9pVwdjwE3uvtN3Ur8FnCXu2fMLL/uccCfgevc/Soz+xvgF8BMYCnwTXdfC/wWuMrM9uo8FmHw3OfuG4FnzWw3M9s3N3GoSK+pRyOyHWa2K3Ar8G1gN8Jnzfw0b5Wv5342lfAD/GQziwO/B+4gnG7914Q9g06fA24Gjjezw4E5hOFyGPBFAHffDDyUWxczqwMOyO23uy/l1s2vewRhKN6VC5naXPsnhJMmPgPclDvWO8CDuf10Op6wp9PpYeAL2/jPJLJNChqR7TsVuN3dl7j7u8CFhB/CnRa4+/Lch/VjhLPgHgRUuPu17r7Z3a8GKsxsv9w2T7v77e6eAr5C2ONYnetd/Cxv33cCx+RefxZ4LFdDl1yoHUjYa+qUIDyNttHdL8wtOxZ41t3vyvVS/hX4hJmNyf286/SZmc0EdmHrUFtO2DsTKYhOnYls327A2Wb27fyFZlade5n/wd9G+O9qN2C6mXWfTLDzYVLv5S2bDDyS127Ie70IuNLMRhL2bO7oob46wj8aN+Ut2wNoBPY0s3HuviFX0+d6qGlnYCNhqF2b68GdACzqNrHiBsKekEhB1KMR2b51wE/dPebuMaAa+Ki7t21nm2Wd2+S225+tA6VTM+FF906TO1+4+zrgL8DRua9FPWyfAGK5r04NhKfingB+nFfTb/PqiROeiluVO9Ym4D7CkDmBrU+bQXhNKfshv7NIj9SjEdm+24F7zWwB4XDfSwl7Jsd/yDZPAOPN7IvAnwg/uC8Ddu9h3QeB75jZ7YT/Jr9HOEy5053ARcAr7v56D9uvJxxEMIaw1wHQ7O4pM7sAeMbMriUMkblmdgjwLHAOcKK7fyxvX7flfr+xhAMb8o0hHPggUhD1aES2w91fIBwRdhfwJrlTadvZpo3wov4PCIPgu8AJ7t7ew+r/BfyR8BrIw7nX6byf30nYG/qfbRwrDSwD9uvhZyuBXwJX5XpHpxGOOnsH+Dvg5G6b/IHwVFrn9aN8MwgDVKQgeh6NSMTMbCfoOk2GmZ0FHO7uJ+faCcLrLdPd/Y1t7OMCYJS7/2sJ63yMcFj386U6hgxN6tGIRO/vgHvMbCcz24twqPOfAMysBjiFcLRZjyGTcz1wTC6Uis7MZgDvKWSkLxQ0ItG7hfCC/yrgScJTZ/+d+9nVwFzCocjblBuufB0fPBVWLOcSnj4UKZhOnYmISEmpRyMiIiWloBERkZJS0IiISEkpaEREpKQUNCIiUlL/Hz4id/qpkggQAAAAAElFTkSuQmCCn”, “text/plain”: [
“<Figure size 432x288 with 1 Axes>”
]
}, “metadata”: {
“needs_background”: “light”
}, “output_type”: “display_data”
}
], “source”: [
“fig, ax = plt.subplots()n”, “n”, “ax.plot(energy_grid, energy_grid * func(energy_grid), red)n”, “n”, “n”, “ax.set_xlabel("energy (keV)")n”, “ax.set_ylabel(r"energy flux (F$_{\nu}$)")n”, “ax.set_xscale(x_scale)n”, “ax.set_yscale(y_scale)n”, “n”
]
}, {
“cell_type”: “markdown”, “id”: “217a2c98”, “metadata”: {
- “papermill”: {
“duration”: 0.111292, “end_time”: “2021-08-17T08:11:36.066097”, “exception”: false, “start_time”: “2021-08-17T08:11:35.954805”, “status”: “completed”
}, “tags”: []
}, “source”: [
“## $\nu$F$_{\nu}$n”, “n”, “The $\nu$F$_{\nu}$ shape of the photon modeln”, “if this is not a photon model, please ignore this auto-generated plot”
]
}, {
“cell_type”: “code”, “execution_count”: 8, “id”: “5db241d5”, “metadata”: {
- “execution”: {
“iopub.execute_input”: “2021-08-17T08:11:36.764717Z”, “iopub.status.busy”: “2021-08-17T08:11:36.761994Z”, “iopub.status.idle”: “2021-08-17T08:11:52.580642Z”, “shell.execute_reply”: “2021-08-17T08:11:52.581565Z”
}, “papermill”: {
“duration”: 16.393249, “end_time”: “2021-08-17T08:11:52.583101”, “exception”: false, “start_time”: “2021-08-17T08:11:36.189852”, “status”: “completed”
}, “tags”: []
}, “outputs”: [
- {
- “data”: {
“image/png”: “iVBORw0KGgoAAAANSUhEUgAAAZkAAAEUCAYAAAD5i0vIAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjQuMywgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/MnkTPAAAACXBIWXMAAAsTAAALEwEAmpwYAAAlSUlEQVR4nO3deZRc5Xnn8W/tvakldUtsEkgI0GNkBMI2i22MAQ84sSF2nDjgZWxPMjFOjgkOmEkmTEJmzJhApDiY42WM7STIASXxJhNvcWwkLAtsaBuBBTxi09KAhNRS70ttd/6o6lapaUndraq+VV2/zzl1ut5bd3kaSvXr931v3RsJggAREZFKiIZdgIiIzF4KGRERqRiFjIiIVIxCRkREKiYedgHVoqOjY5jCf4+usGsREakh7UD29a9/fcNELypkDooDsXg8flzYhYiI1IpsNnvE1xUyB3XF4/HjzjnnnLDrEBGpGVu2bCGbzR52BEhzMiIiUjEKGRERqRiFjIiIVIxCRkREKkYhIyIiFaOQERGRilHIiIhIxeh7MiIis1A+CBgOsgzk0wwEGQaLPwfyaQaDDAP5DINBmsF8YVkA3NR2UdnrUMiIiFSRfBAwFGQKH/5HCIeB/MHXCuumDwmOwSDDVO4WliCqkBERqVZBEDB0lJ7D6PNyhsPhNETiNEcSNEWTNEUSNEcTNJc8b4okaYomaCn+bIokCIKASCRShqMfVPMhY2Zx4FLgCne/Kex6RKT2lAZEf5AuBMVYABQCoT8/+jxNfzEwRgNiNDzyZYiH8eHQEj0YAqPh0BxN0BwpXZ485LWmSIJYpDqm3Ks6ZMxsJXCbu19ZbF8P3ABkgZvdfR2wDLgEKG/8ikhNCIpzD/1HCIiBfLoQEsXnB8OkfAGRisTGPvibI8lXhUPpa6PLD65TfeFQLlUbMmZ2O/A+4PFiewVwLbACaAI6zOz77r7NzO4GPj6FfbdTuDz1mLVr18ZisVi5yheRKcgHAYNBIQj68mn68iOF50Hh+cD45fkR+sfWHyF3jAGRJDbWGxgNiNHnzdFie6zHUHjeMva8sG4ios+PiVRtyACbgD7gwmL7KuA+dx8ABsxsI3AZ8K1p7Ps64JbSBf39/TQ2Nh5DuSKSDfL05kfozY/Qlx+hNz881j64PE3/WFgUng8cQ08iRoR50YZiECTH5h5Gnx9cXmxHDgbFaGgkFRAVU7Uh4+73m1lpyCwFHipZpRM4aZq7vwu4t3RBS0vL5mQy2X6Y9UXqThAE9AdpDuSGOJAfpjs3zIH8EN25YXpKwqOvJEAGg8y0jpWKxGiJJGmJppgTLf2ZZE4kVfhZbI++Nrq8IRIv+2S1lE/VhswE8kBmXDsL4O7bgU9Odkfu3sW4O2B2dHTk9EaV2SwfBHTnh9mfG6QrN8RAkCYb5NlfDI4DJT8PFIMkS35Kx2iOJGmNJmmNNtAaTTEnmqK15DFnNCCipcGRUk9iFqulkOkEFpe0FwEPhlSLSFVJBzm6coPszQ2wd+znAHuzhfa+3CAH8kNTGpKKE2VBrIn50UbmRxuYF2tgfrSRebEG5kYbmBNNMfeQ8EgRn2WT1nLsailk1gP3mNlngQXABcAfhVuSyMxIB1l2Zwd4OdvH7lwfL2X7eDnbx57cAPtyAxzIDx91HwmiLIw10x5rpC3aSEs0RSISLYRIrIF54362RJIahpJjVjMh4+5Pmtk6YCuFobJPuPtQyGWJlMVQPsOubA+7sr1054bozg+zLzfIS8VQ2ZcbPGIfpDmSYEGsmYWxJhbGmlkYa2bB6PN4EwtizcxRaEgIqjpk3H0DsKGkvRpYHVY9IseqJzfM9mw3OzPd7Mj2sDPTzc5sD6/kBo64XVMkwYnxOZwYayn8jM/hxNgcjo+3sDDWRHM0OUO/gcjUVHXIiNSqIAh4OdfPs5kunk3v59nMfp7LdLE3Nzjh+s2RJEsSczk5Ppf2WBNzow20xxo5qRgmrdGUeiFSkxQyImWwLzfIk+lXeHpkL0+l9/FspouBCU7nnRdtYFliPqfE542FypLEPNqijQoRmZUUMiJTlAvyPJfZz5aR3WxN7+Wp9Cuv6qFEgMXxVk5PtHF6op3TEm2cnmxTmEjdUciIHEU+CHghc4DHRl7msZHdbBnZTX+QPmSd9mgTZyYXsCJ1HGcmF3JGop2maCKkikWqh0JGpGhvboBdmR4WxVtJRmI8OvwSj4x08ujwS3SPO0V4WWI+q1InsDJ5AiuSC1kYbw6papHqppCRupUPAp5K7+VnQzvYPLyLndmew667JD6XVakTWZU6kXNSJzAv1jCDlYrULoWM1JUgCHg6s48Ngy+wcWj7IacOt0ZT9OZHAGiJJHl9w0mc37CINzQsYmFMPRWR6VDISF3Ymenhh4PPsGHwBV7O9Y8tXxKfy8WNS3lz4ymckWinP0izO9vPssT8WXdfD5EwKGRkVunJDfODwWd4bGQ3759zNi9me/n+wDM8kd4zts7ieCuXNJ7KpU2nsjQ+75CzveZEUsxJpsIoXWRWUsjIrPDUyF6+2f8kDw5tJ1O8cvDPhzvHXm+Nprii6XQubzqN0xNtOo1YZIYoZKRm5YI8Pxveydf7tvLr9CtA4crBq1In8NjIbgDOSy3iN5vP4E2Np+hy8iIhUMhIKPrzaZ7P7Oes5PFEp9iryAQ5fjDwDP/S92teyvUB0BZt5N0tZ/LO5uXMjzXycraPOFGdWiwSMoWMzKjhfJZvDzzFur4n6M2P8Kn2y3hz45JJbZsJcvxw4Fm+1rdl7KywZYn5vLfltVzatOyQnsqJ8TkVqV9EpkYhIzMiCAJ+MvQCX+p55JBLsOzLHf1uDbkgz38MPss9vVvYUzwz7MzkQj7cuorzUos0vyJSxRQyUnFPp/fyue5fsLU4b3J28nhaoyk2De8kHeSOuO0vh1/i8z2/4PnMAQAssYCPtJ7L+Q0KF5FaoJCRihnMZ7i751HWDzwNwAmxFq6dex4XNy7hn/seL4ZMdsJtd2V6+GLPIzw0vAuAU+Jz+ejcN/DGhpMVLiI1RCEjFfHzoU4+072ZV3IDpCIxPjDnHH5vzmtJRgpvuVRx/mRkXE8mHWS5t/cJ7u17nCx5WqMpPtJ6Llc2m+4fL1KDFDJSVkP5DHd1/5wfDD4DwKrUCdw4/80sircest5o2IyU9GR+OfwSf9/9EJ3ZXqJE+J2WFXyodRVzovpypEitUshI2WxL7+PW/RvpzPbSHEnwsXnn846mMyYc3hrtyaSDHIP5DJ/v+QXfG9gGFOZdbpj/Js5Its9o/SJSfjUdMmZ2FvApoAf4jrt/M+SS6lIQBHyj/0m+1PMoWfKsSC7k5ra3HvE04tGQeSbTxR/uWc/LuT4aI3H+YO7reVfza3TdMJFZoipCxsxWAre5+5XF9vXADUAWuNnd1x1m0/OgeA0ReKLihcqrDOez/O2BTTww9AIR4ANzzubDrecedf5kdLjsyfReAM5KHseft13MSfp+i8isEnrImNntwPuAx4vtFcC1wAqgCegws++7+0Q3+3gEWE/hbrefA66ZkaIFgN3ZPv6q6yc8m9lPazTFX7VdwusaTprUtq3FeZY4UT7Sei5XzzlLvReRWSj0kAE2AX3AhcX2VcB97j4ADJjZRuAyM1sGXF6y3QeBc4FtxXaeSTKzduCQAf+1a9fGYjFd22qynhrZy190/Yie/AinJdr4VPtlnDCFXsiK5EJunP9mXpNYwGnJtgpWKiJhCj1k3P1+MysNmaXAQyWrdAInufsaYE3ptmb2MrAWGAI+O4XDXgfcUrqgv7+fxsbGqRVfpx4e2sX/2b+B4SDLWxuX8mfz30JDdGpvpVgkyjubl1eoQhGpFqGHzATyQGZce8Jv7Ln7fwL/OY1j3AXcW7qgpaVlczKp05mO5gcDz7D6wM/IE/C7La/lY3PPm/IFLkWkflRjyHQCi0vai4AHy3kAd+8CukqXdXR05PRN8iNb3/8Ud3Y/DMC1c9/A1XNWhlyRiFS7agyZ9cA9ZvZZYAFwAfBH4ZZUv9JBjv78CD8d2sGd3Q8TAW6afxG/0XxG2KWJSA2oupBx9yfNbB2wlcJQ2Sfc/eiX6pWy25cb5BOvfG/sni0R4M/mv4Urmk8PtzARqRlVETLuvgHYUNJeDawOqx6B3vwI/2PvD8cCBhQwIjJ1VREyUl3SQZa/3Pdjtme7WZaYz6rUibw+dRJvbDw57NJEpMYoZOQQ+SDg9v2beCK9h0XxOaxZ8BvMjTWEXZaI1Ch9xVoO8dXeX/LA0Au0RlPc1n65AkZEjolCRsb8x8Cz3Nv3OAli3Nr+NhYn5oZdkojUOIWMAPBMuou/O7AZgD9vu4izUseHXJGIzAYKGaE3P8Jfd/2ENDl+r+UsLm1aFnZJIjJLKGTqXD4I+PT+jbyc62dV6gT+cO7rwy5JRGYRhUyd+3r/Vn4x/CILYk38r7ZLdLl9ESkrfaLUsWfTXXy5p4MIcHPbW2mL6SrUIlJeCpk6NZzPcuv+jWTJ8/45Z3NO6oSwSxKRWUghU6e+2PMIO7M9WGIBH249N+xyRGSWUsjUoZ8PdfKdgadpiMS5ue1i4pqHEZEK0adLnenNj7D6wCYA/mjuefrCpYhUlEKmznyu++d05Yc4P7WIK5st7HJEZJZTyNSRLSO7+dHgczRHEnyy7SJ0J1ARqTSFTJ3IBnnuPPAQAP+t9XUsiDWFXJGI1AOFTJ34Rv9Wtme7OS3RxrtaXhN2OSJSJxQydWBvdoB/6n0MgOvnXahv9YvIjNGnTR34cm8Hw0GWtzedrqsri8iMUsjMctvS+/jR4HM0ROK6+KWIzLiavP2ymcWBS4Er3P2m8e1wq6seQRDwxZ5HAHjfnJW0abJfRGZYqD0ZM1tpZv9e0r7ezHaY2XNmds0RNl0GXAJEDtMW4OHhTh4b2U17tInfbXlt2OWISB0KLWTM7HbguyXtFcC1wArgQuAOM5vw6+juvg24+3DtSRy73cyWlz5yuVwsCIJp/jbVJxfk+VLPowD8/tzX0RhNhFyRiNSjMIfLNgF9FAIF4CrgPncfAAbMbCNwmZktAy4v2e6D7r7vGI99HXBL6YL+/n4aG2fPpe43Dm1nR7abJfF5XNF0WtjliEidCi1k3P1+MysNmaXAQyWrdAInufsaYE2ZD38XcG/pgpaWls3JZLK9zMcJRS7Ic0/xlOUPt67SKcsiEppqmvjPA5lx7ezhVnb37cAnD9c+EnfvArpKl3V0dORmy2VWNgxtZ2e2h1Pj87i4cWnY5YhIHaumP3E7gcUl7UXA9nBKqV25IM/aYi/mQ63nEp0lwSkitamaQmY9cLWZpcxsEXAB8GDINdWcDUMvsDPbw7LEfN7SuCTsckSkzlVNyLj7k8A6YCvwAHCjuw+FW1VtCYKAe3sfB+C/zlmlXoyIhC7UORl33wBsKGmvBlaHVU+te2TkRV7IdrMoPoeLGk8JuxwRkerpycix+5e+XwPw3pazdEaZiFQFfRLNEtvS+/jVyMvMjaZ4e/PpYZcjIgIoZGaNf+3bCsC7W84kFammM9NFpJ4pZGrUUD7Dfb2Psy83yO5sPxuGXiBJjHc3nxl2aSIiY/Qnb436X10/5lcjL/NC9gALYs3kCXh78+nMjTWEXZqIyBiFTA1KBzl+NfIyALsyvTwy/CJQGCoTEakmCpka9ODQ9rHnu7I9DAYZzk4ez6mJ+eEVJSIyAc3J1KDv9m8bez4YFC739q6W14RVjojIYSlkasyL2V62pHcfsmx+tJGLdAkZEalCCpka84OBZwBYlTphbNk7m5eTiMTCKklE5LAUMjUkF+T5wcCzALynZQUAUSJc2bw8zLJERA5LE/815NHhl+jKD3J6oo03NZzCWxuXsiQ+j+PiLWGXJiIyIYVMDfn+YGGo7B3Ny4lGItzSfmnIFYmIHJmGy2pET26YzUM7SRDjbU3Lwi5HRGRSFDI14qdDO8iS542NJzMnmgq7HBGRSVHI1IgfDz4PoF6MiNQUhUwNeCXbz+Pp3TRHklzQsCjsckREJk0hUwMeGHqBALi4cQlJXcZfRGqIQqYGjA6VXaahMhGpMTX1Z7GZxYFLgSvc/SYzuxlYOfq6u18TWnEVsiPTzbOZ/cyPNh7yLX8RkVowrZAxs7cAZwCb3H3b0dafYPuVwG3ufmWxfT1wA5AFbnb3dYfZdBlwCRABcPf/W9z+j4HHplpHLXhg8AUALm1aSiyijqeI1JbpfmpdBAwAt5jZLVPZ0MxuB75b0l4BXAusAC4E7jCzuRNtWwy0u8ft73hgubtvntJvUCM2De8A4JLGU0OuRERk6qbVk3H324pP/8XMPmpma9z9xkluvgnooxAoAFcB97n7ADBgZhuBy8xsGXB5yXYfdPd9E+zvT4EvTKV+M2sH2kuXrV27NhaLVddFJjszPTyfOUBbtJEVyePCLkdEZMqmO1z2ZeBE4AQgAQxPdlt3v9/MSkNmKfBQySqdwEnuvgZYM4ldnunuPtnjF10HHNID6+/vp7GxcYq7qaxNwzsBeHPjKUQjkZCrERGZuulO/H8eeBnY4+75Y6whD2TGtbOHW9ndtwOfLGm/axrHvAu4t3RBS0vL5mQy2X6Y9UPx06HCUNlbdK8YEalR0x0u+2UZa+gEFpe0FwEPlnH/r+LuXUBX6bKOjo5cpIp6C3uzAzyV3ktLJMmq1IlhlyMiMi1HnPg3s53j2r9VgRrWA1ebWcrMFgEXUOGQqQWjE/5vajyZuM4qE5EadbRPr8Xj2t8qdwHu/iSwDtgKPADc6O5D5T5OrdFQmYjMBkcbLgvGtcsynuTuG4ANJe3VwOpy7Hs26MuP8PjIHhoicd6ga5WJSA2b6jjM+NCRCnhk+EXyBJybOpGUrlUmIjXsqJ9gZraQQg8mMq4NgLu/UrHq6tTDw50AXNgwfrRSRKS2HC1kIsDu4s+g+HNPyesBUF3fYKxxuSDPI8WQuUAhIyI17ogh4+5RADNr5uCXJ3/q7ulKF1avPL2PnvwIyxLzOS7eEnY5IiLHZLJzMtuBvwJuA140s/dWrKI6d3Co7OSQKxEROXaTnVU+zd17AczsDcA9Znaau/9N5UqrTw8P7wI0HyMis8OkejKjAVN8/ihwPvA2M/uLShVWj/blBnk2s5/WaIozkwvDLkdE5JhNqidjZm+ncLmX0cfi4uNTwKcrVl2d+UVxqOy81CLdO0ZEZoXJDpd9AXip+HgR2AjcV3wuZfLQUHGorFHzMSIyO0wqZNxdN5evsHSQpWPkJaJEOE/f8heRWUJjMlViy8gehoMsr00eR2s0FXY5IiJloZCpEjqrTERmI4VMlfjl8MsAGioTkVlFIVMF9ueG2JHtpjWaYlmiLexyRETKRiFTBR4bKfRiVqVOIFpFd+cUETlWCpkqcDBkdJtlEZldFDJV4FfFkDlXISMis4xCJmSvZPt5MdtHW7SRU+Jzwy5HRKSsFDIhe2xkN1AYKotoPkZEZhmFTMjGhsoaNFQmIrNPTd1A3swuBj5ebH6ewl05Pwo0AZ9294fCqm06giDQfIyIzGqhhIyZrQRuc/cri+3rgRuALHCzu687zKYXAx8CksCdwGPAB4CzgXcCNRUye3L9vJIb4LhYMyfGdBdMEZl9Zny4zMxuB75b0l4BXAusoHCL5zvMbMIZcHe/FZgH3AHc7e53Am8GvgxsmkIN7Wa2vPSRy+ViQRBM99ealsdH9gBwdup4zceIyKwURk9mE9BHIVAArgLuc/cBYMDMNgKXmdky4PKS7T4InFNc/2Z37zKz33f3r5rZhcA/AQ9MsobrgFtKF/T399PY2DjtX2o6Hk8XQmZl8vgZPa6IyEyZ8ZBx9/vNrDRklnLoMFcncJK7rwHWlG5rZn8CdANrzGwzkDOzrwENwD1TKOMu4N7SBS0tLZuTyWT7FPZxzJ4o9mRWpk6YycOKiMyYapj4zwOZce3sRCu6+7smWPyVqR7Q3buArtJlHR0duZkcsjqQG2JXtoe50RRL9P0YEZmlquEU5k4Kt3IetQjYHk4pM+eJ4lDZWUnNx4jI7FUNIbMeuNrMUma2CLgAeDDkmiruiZJJfxGR2Sr0kHH3J4F1wFYKE/c3uvtQuFVV3uNj8zEKGRGZvUKZk3H3DcCGkvZqYHUYtYRhMJ/hucx+GiJxzkjM6LkGIiIzKvSeTD16Or2XPAErkguJRfS/QERmL33CheDJ9F4AzkwuDLkSEZHKUsiE4CmFjIjUCYXMDAuCQCEjInVDITPDduf66c4Pc2Kshfmxmb2MjYjITFPIzLDR+ZjXqBcjInVAITPDRofKVihkRKQOKGRm2NOajxGROqKQmUHpIMcz6S4SRDl9Zi/4LCISCoXMDHous58MeU5LtJGMxMIuR0Sk4hQyM+ipkeJQWUpDZSJSHxQyM+jpjOZjRKS+KGRm0DPpwn3SlicWhFyJiMjMUMjMkKF8hp3ZHhojcRbHW8MuR0RkRihkZsizmf0EwOmJdqK6E6aI1AmFzAwZGyrTqcsiUkcUMjNkW2YfoPkYEakvCpkZsq3YkzlDPRkRqSMKmRkwnM+yM9tDQyTOyZr0F5E6Eg+7gKkws4uBjxebnwfagA8A/cBqd38irNqO5LnMfvIEnJ5o0+2WRaSuhBIyZrYSuM3dryy2rwduALLAze6+7jCbXgx8CEgCdwItwDXAHOBzwPsqXPq0jM7HaKhMROrNjP9ZbWa3A98taa8ArgVWABcCd5jZ3Im2dfdbgXnAHcDdwGrgLgoB1VTRwo/BM+n9ACxPKGREpL6E0ZPZBPRRCBSAq4D73H0AGDCzjcBlZrYMuLxkuw8C5xTXv9ndu8zsD4DrgDOAj062ADNrBw75xF+7dm0sFqvMRSvHzixL6swyEakvMx4y7n6/mZWGzFLgoZJVOoGT3H0NsKZ0WzP7E6AbWGNmm4EdwD8AjRR6M5N1HXBL6YL+/n4aG8t/O+R0kGV7pptUJMYp8Qk7aCIis1Y1TPzngcy4dnaiFd39XRMs/uE0jnkXcG/pgpaWls3JZPknTbZnuskTcGp8vib9RaTuVEPIdAKLS9qLgAcreUB37wK6Spd1dHTkIhW43MvzmQMALEu2lX3fIiLVrhr+tF4PXG1mKTNbBFxAhUNmJj2XKUz6L0vMD7kSEZGZF3rIuPuTwDpgK/AAcKO7D4VbVfmM9mROS6gnIyL1J5ThMnffAGwoaa+mcDryrBIEgXoyIlLXQu/JzGZd+SF68yMcF2tmTjQVdjkiIjNOIVNBz4/1YjRUJiL1SSFTQQfnYzRUJiL1SSFTQc+lNR8jIvVNIVNBGi4TkXqnkKmQdJBjZ7aHJDEW6x4yIlKnFDIVsjPTTY6ApYl5upyMiNQtffpVyNjlZDRUJiJ1TCFTIaMhc6om/UWkjilkKmRHthuApYl5odYhIhImhUyF7Mh0AwoZEalvCpkKGMpn2J3rpzmSYEG0au8KLSJScQqZCtiZ7QFgSWIelbhHjYhIrVDIVMDYUFl8Xqh1iIiETSFTAduLZ5Yt0XyMiNQ5hUwFbC8Oly3V6csiUucUMhWwY7QnE58bciUiIuFSyJTZ6JllTZEEC2PNYZcjIhIqhUyZ7cr2EKAzy0REQCFTdtuLZ5Yt0ZllIiLEwy5gKszsbcBHgSbg00Af8CmgB/iOu38zxPIAXU5GRKRUKCFjZiuB29z9ymL7euAGIAvc7O7rDrPpWcAHgLOBdwKdQL742hMVLXqS1JMRETloxkPGzG4H3gc8XmyvAK4FVlDooXSY2ffdvWf8tu5+p5m9FfgMcCOwF1gPRIDPAddMsoZ2oL102dq1a2OxWGy6v9YYXbNMROSgMHoymygMc11YbF8F3OfuA8CAmW0ELjOzZcDlJdt9EPgtd/+qmV0I/BPwPWBb8fU8k3cdcEvpgv7+fhobG6f8y5RKBzl25/ppiMR1ZpmICCGEjLvfb2alIbMUeKhklU7gJHdfA6wp3dbMImb2NaABuAcYBNYCQ8Bnp1DGXcC9pQtaWlo2J5PJ9sOsPykvZnvJE7A43kpUZ5aJiFTFxH8eyIxrZyda0d2/Anxl3OL/nOoB3b0L6Cpd1tHRkTvWU447s70AnKwvYYqIANVxCnMnsLikvQjYHk4px2ZXpjCNtDjeGnIlIiLVoRpCZj1wtZmlzGwRcAHwYMg1Tcuu4jXLTk6oJyMiAlUQMu7+JLAO2Ao8ANzo7kPhVjU9YyGj4TIRESCkORl33wBsKGmvBlaHUUs5jc7JaLhMRKQg9J7MbNGTG6Y3P0JbtJHmaDLsckREqoJCpkw0VCYi8moKmTLZNXr6ckJDZSIioxQyZdKpnoyIyKsoZMrk4HdkFDIiIqMUMmUyNlymM8tERMYoZMogF+R5KdtLjAgnxueEXY6ISNVQyJTBntwAGfKcFG8lFtF/UhGRUfpELIPR+RgNlYmIHEohUwbd+SESRHVmmYjIONVwqf+a9/bmM/gvTaeRDnJhlyIiUlUUMmUSi0Rp1HyMiMgh9KkoIiIVo5AREZGKUciIiEjFKGRERKRiFDIiIlIxChkREakYncJ8UHs2m2XLli1h1yEiUjOy2SxA++FeV8gclAXIZrNd091BOp1OJZPJkfKVVB3HLse+p7uPqW432fXLsV4+n4+k0+mmZDI5GI1Gg8nWWAv0Xi7/Pir1Xp7suhV8L7dT/PycUBAEepTpsXz58ntn47HLse/p7mOq2012/XKst7wgWL58+fKw/r9X8//zajz2bHwvT3bdsN7LmpMREZGKUciU1zdm6bHLse/p7mOq2012/XKvN9vovVz+fVTqvTzZdUP5fxoJglk1lCwyY8xsOeCAufu2sOsRma5KvpfVkxERkYpRyIhMXxfwv4s/RWpZxd7LGi4TEZGKUU9GREQqRiEjIiIVo5AREZGKUciIiEjFKGRERKRiFDIiIlIxChkREakYhYyIiFSMQkZERCpGISMiIhWjO2OKlImZrQRuBQaAz7j7IyGXJDJtZhYHLgWucPebprsf9WREJsHMVprZv5e0rzezHWb2nJldU1z8MeD3i49PhFCmyKRM8v28DLgEiBzLsRQyIkdhZrcD3y1prwCuBVYAFwJ3mNlcoN3du9x9GEiGUqzIUUz2/Vy8r8zdx3o8DZeJHN0moI/CP0CAq4D73H0AGDCzjcBlQJ+ZzQfSQH8olYoc3WTfz98qx8HUkxE5Cne/n8I/zFFLgR0l7U7gJODvgS8VH6tnqDyRKZnC+7ks1JMRmbo8kBnXzrr7VuC94ZQkMm0Tvp8B3H078Mlj2bl6MiJT1wksLmkvAraHU4rIMavo+1khIzJ164GrzSxlZouAC4AHQ65JZLoq+n5WyIhMkbs/CawDtgIPADe6+1C4VYlMT6Xfz5EgCMq1LxERkUOoJyMiIhWjkBERkYpRyIiISMUoZEREpGIUMiIiUjEKGRERqRhdVkZkFjCz64HdwB7gi+7+milsuwn4sbvfUrIsAuwE/hp4M3Cruz9f1qKlLqgnIxISM4uVaT+twIeAr09zF/cB7xm37HzgOOAbwGeBO6ZdoNQ19WREjsLM3gZ8BjiNwn04Pubu+83sH4FXKPylfw7wH8D73X3YzJZSuBfHGyl8k/qj7r7FzD4C/A4QUPgj70oz+1Pgz4Ah4J+Bi4DfAvYCZ7n7M8U6fgZ8xd2/Oq7EPwS+7e45Myutuw34KfAld7/TzFYBXwRWApuB/+7uO4B/A+40szNGj0UhdL7v7t3AY2a2xMxeW7wIqMikqScjcgRmdjJwL/DHwBIK94r525JVPlx8bRmFD+9rzCwKfAf4JoVLpv8zhR7BqHcA9wC/bWaXAtdTCJZLgHcDuHsv8JPiuphZO/CG4n7H+53iuqV1N1EIxG8XA6al2P4bChdA/CXw1eKxXgF+XNzPqN+m0MMZ9QDwrsP8ZxI5LIWMyJF9APi6u29y933AzRQ+gEetdfctxQ/qn1G4mu35QMzdv+Duve7+WSBmZmcXt3nU3b/u7hng/RR6Gs8WexV/V7LvbwHvLD7/DeBnxRrGFAPtPAq9pVFxCkNn3e5+c3HZlcBj7v7tYu/kL4E3mtm84utjQ2ZmthI4kUMDbQuFXpnIlGi4TOTIlgAfM7M/Ll1oZg3Fp6Uf+sMU/k0tAVaY2fgLA47eCOpAybJTgI0l7c6S5+uBvzezZgo9mm9OUF87hT8We0qWnQZ0AaebWZu77y/W9I4JajoB6KYQaF8o9tzeA6wfd5HE/RR6QCJTop6MyJHtAf7W3SPuHgEagHPdffgo2zwyuk1xu3M4NExGDVCYYB91yugTd98D/Aq4ovhYP8H2cSBSfIzqpDD89jDwqZKa/q2kniiF4bdtxWP1AN+nEDDv4dChMijMIeWP8DuLTEg9GZEj+zrwPTNbS+GU3k9T6JH89hG2eRg4zszeDfyIwof27cCpE6z7Y+DjZvZ1Cv8eP0nhVORR3wJuAZ53910TbL+XwgkD8yj0NgAG3D1jZjcBvzSzL1AIkNvM7CLgMeA64L3u/rqSfd1X/P3mUziJodQ8Cic5iEyJejIiR+Duv6Zw5te3gZcoDp8dZZthChP4/5NCCPwp8B53H5lg9f8H/JDCnMcDxefZkte/RaEX9K+HOVYWeAQ4e4LXngT+Abiz2Cv6IIWzy14BfhO4Ztwm/05h+Gx0vqjUWRTCU2RKdD8ZkRCZ2fEwNjSGmV0LXOru1xTbcQrzKyvc/cXD7OMmoNXd/7KCdf6MwqnbT1TqGDI7qScjEq7fBL5rZseb2RkUTmf+EYCZNQLvo3BW2YQBU3Q38M5iIJWdmZ0FHFDAyHQoZETC9TUKk/vbgJ9TGC77x+JrnwVuo3C68WEVT0n+Eq8e/iqXT1AYMhSZMg2XiYhIxagnIyIiFaOQERGRilHIiIhIxShkRESkYhQyIiJSMf8fLZt6j+EAiEYAAAAASUVORK5CYII=n”, “text/plain”: [
“<Figure size 432x288 with 1 Axes>”
]
}, “metadata”: {
“needs_background”: “light”
}, “output_type”: “display_data”
}
], “source”: [
“fig, ax = plt.subplots()n”, “n”, “ax.plot(energy_grid, energy_grid**2 * func(energy_grid), color=green)n”, “n”, “n”, “ax.set_xlabel("energy (keV)")n”, “ax.set_ylabel(r"$\nu$F$_{\nu}$")n”, “ax.set_xscale(x_scale)n”, “ax.set_yscale(y_scale)n”
]
}
], “metadata”: {
- “jupytext”: {
“formats”: “ipynb,md”
}, “kernelspec”: {
“display_name”: “Python 3”, “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.7.11”
}, “papermill”: {
“default_parameters”: {}, “duration”: 28.237502, “end_time”: “2021-08-17T08:11:55.707801”, “environment_variables”: {}, “exception”: null, “input_path”: “PhAbs.ipynb”, “output_path”: “../docs/notebooks/PhAbs.ipynb”, “parameters”: {
“func_name”: “PhAbs”, “linear_range”: false, “wide_energy_range”: false, “x_scale”: “log”, “y_scale”: “log”
}, “start_time”: “2021-08-17T08:11:27.470299”, “version”: “2.3.3”
}
}, “nbformat”: 4, “nbformat_minor”: 5
}