crp/OpenMC
1
0
Fork 0
mirror of https://github.com/openmc-dev/openmc.git synced 2026-07-23 19:45:34 -04:00
Code Issues Projects Releases Packages Wiki Activity Actions 5
OpenMC/docs/source/pythonapi/examples/pandas-dataframes.ipynb
Will Boyd 5cc884555b Updated IPython Notebooks
2016-03-23 14:51:39 -04:00

2466 lines
135 KiB
Text
Raw Blame History

{
"cells": [
{
"cell_type": "markdown",
"metadata": {},
"source": [
"This notebook demonstrates how systematic analysis of tally scores is possible using Pandas dataframes. A dataframe can be automatically generated using the `Tally.get_pandas_dataframe(...)` method. Furthermore, by linking the tally data in a statepoint file with geometry and material information from a summary file, the dataframe can be shown with user-supplied labels.\n",
"\n",
"**Note:** that this Notebook was created using the latest Pandas v0.16.1. Everything in the Notebook will wun with older versions of Pandas, but the multi-indexing option in >v0.15.0 makes the tables look prettier."
]
},
{
"cell_type": "code",
"execution_count": 1,
"metadata": {
"collapsed": false
},
"outputs": [],
"source": [
"import glob\n",
"from IPython.display import Image\n",
"import matplotlib.pylab as pylab\n",
"import scipy.stats\n",
"import numpy as np\n",
"\n",
"import openmc\n",
"from openmc.statepoint import StatePoint\n",
"from openmc.summary import Summary\n",
"from openmc.source import Source\n",
"from openmc.stats import Box\n",
"\n",
"%matplotlib inline"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"## Generate Input Files"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"First we need to define materials that will be used in the problem. Before defining a material, we must create nuclides that are used in the material."
]
},
{
"cell_type": "code",
"execution_count": 2,
"metadata": {
"collapsed": true
},
"outputs": [],
"source": [
"# Instantiate some Nuclides\n",
"h1 = openmc.Nuclide('H-1')\n",
"b10 = openmc.Nuclide('B-10')\n",
"o16 = openmc.Nuclide('O-16')\n",
"u235 = openmc.Nuclide('U-235')\n",
"u238 = openmc.Nuclide('U-238')\n",
"zr90 = openmc.Nuclide('Zr-90')"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"With the nuclides we defined, we will now create three materials for the fuel, water, and cladding of the fuel pin."
]
},
{
"cell_type": "code",
"execution_count": 3,
"metadata": {
"collapsed": false
},
"outputs": [],
"source": [
"# 1.6 enriched fuel\n",
"fuel = openmc.Material(name='1.6% Fuel')\n",
"fuel.set_density('g/cm3', 10.31341)\n",
"fuel.add_nuclide(u235, 3.7503e-4)\n",
"fuel.add_nuclide(u238, 2.2625e-2)\n",
"fuel.add_nuclide(o16, 4.6007e-2)\n",
"\n",
"# borated water\n",
"water = openmc.Material(name='Borated Water')\n",
"water.set_density('g/cm3', 0.740582)\n",
"water.add_nuclide(h1, 4.9457e-2)\n",
"water.add_nuclide(o16, 2.4732e-2)\n",
"water.add_nuclide(b10, 8.0042e-6)\n",
"\n",
"# zircaloy\n",
"zircaloy = openmc.Material(name='Zircaloy')\n",
"zircaloy.set_density('g/cm3', 6.55)\n",
"zircaloy.add_nuclide(zr90, 7.2758e-3)"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"With our three materials, we can now create a materials file object that can be exported to an actual XML file."
]
},
{
"cell_type": "code",
"execution_count": 4,
"metadata": {
"collapsed": false
},
"outputs": [],
"source": [
"# Instantiate a MaterialsFile, add Materials\n",
"materials_file = openmc.MaterialsFile()\n",
"materials_file.add_material(fuel)\n",
"materials_file.add_material(water)\n",
"materials_file.add_material(zircaloy)\n",
"materials_file.default_xs = '71c'\n",
"\n",
"# Export to \"materials.xml\"\n",
"materials_file.export_to_xml()"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"Now let's move on to the geometry. This problem will be a square array of fuel pins for which we can use OpenMC's lattice/universe feature. The basic universe will have three regions for the fuel, the clad, and the surrounding coolant. The first step is to create the bounding surfaces for fuel and clad, as well as the outer bounding surfaces of the problem."
]
},
{
"cell_type": "code",
"execution_count": 5,
"metadata": {
"collapsed": false
},
"outputs": [],
"source": [
"# Create cylinders for the fuel and clad\n",
"fuel_outer_radius = openmc.ZCylinder(x0=0.0, y0=0.0, R=0.39218)\n",
"clad_outer_radius = openmc.ZCylinder(x0=0.0, y0=0.0, R=0.45720)\n",
"\n",
"# Create boundary planes to surround the geometry\n",
"# Use both reflective and vacuum boundaries to make life interesting\n",
"min_x = openmc.XPlane(x0=-10.71, boundary_type='reflective')\n",
"max_x = openmc.XPlane(x0=+10.71, boundary_type='vacuum')\n",
"min_y = openmc.YPlane(y0=-10.71, boundary_type='vacuum')\n",
"max_y = openmc.YPlane(y0=+10.71, boundary_type='reflective')\n",
"min_z = openmc.ZPlane(z0=-10.71, boundary_type='reflective')\n",
"max_z = openmc.ZPlane(z0=+10.71, boundary_type='reflective')"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"With the surfaces defined, we can now construct a fuel pin cell from cells that are defined by intersections of half-spaces created by the surfaces."
]
},
{
"cell_type": "code",
"execution_count": 6,
"metadata": {
"collapsed": false
},
"outputs": [],
"source": [
"# Create a Universe to encapsulate a fuel pin\n",
"pin_cell_universe = openmc.Universe(name='1.6% Fuel Pin')\n",
"\n",
"# Create fuel Cell\n",
"fuel_cell = openmc.Cell(name='1.6% Fuel')\n",
"fuel_cell.fill = fuel\n",
"fuel_cell.region = -fuel_outer_radius\n",
"pin_cell_universe.add_cell(fuel_cell)\n",
"\n",
"# Create a clad Cell\n",
"clad_cell = openmc.Cell(name='1.6% Clad')\n",
"clad_cell.fill = zircaloy\n",
"clad_cell.region = +fuel_outer_radius & -clad_outer_radius\n",
"pin_cell_universe.add_cell(clad_cell)\n",
"\n",
"# Create a moderator Cell\n",
"moderator_cell = openmc.Cell(name='1.6% Moderator')\n",
"moderator_cell.fill = water\n",
"moderator_cell.region = +clad_outer_radius\n",
"pin_cell_universe.add_cell(moderator_cell)"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"Using the pin cell universe, we can construct a 17x17 rectangular lattice with a 1.26 cm pitch."
]
},
{
"cell_type": "code",
"execution_count": 7,
"metadata": {
"collapsed": false
},
"outputs": [],
"source": [
"# Create fuel assembly Lattice\n",
"assembly = openmc.RectLattice(name='1.6% Fuel - 0BA')\n",
"assembly.dimension = (17, 17)\n",
"assembly.pitch = (1.26, 1.26)\n",
"assembly.lower_left = [-1.26 * 17. / 2.0] * 2\n",
"assembly.universes = [[pin_cell_universe] * 17] * 17"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"OpenMC requires that there is a \"root\" universe. Let us create a root cell that is filled by the pin cell universe and then assign it to the root universe."
]
},
{
"cell_type": "code",
"execution_count": 8,
"metadata": {
"collapsed": false
},
"outputs": [],
"source": [
"# Create root Cell\n",
"root_cell = openmc.Cell(name='root cell')\n",
"root_cell.fill = assembly\n",
"\n",
"# Add boundary planes\n",
"root_cell.region = +min_x & -max_x & +min_y & -max_y & +min_z & -max_z\n",
"\n",
"# Create root Universe\n",
"root_universe = openmc.Universe(universe_id=0, name='root universe')\n",
"root_universe.add_cell(root_cell)"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"We now must create a geometry that is assigned a root universe, put the geometry into a `GeometryFile` object, and export it to XML."
]
},
{
"cell_type": "code",
"execution_count": 9,
"metadata": {
"collapsed": true
},
"outputs": [],
"source": [
"# Create Geometry and set root Universe\n",
"geometry = openmc.Geometry()\n",
"geometry.root_universe = root_universe"
]
},
{
"cell_type": "code",
"execution_count": 10,
"metadata": {
"collapsed": false
},
"outputs": [],
"source": [
"# Instantiate a GeometryFile\n",
"geometry_file = openmc.GeometryFile()\n",
"geometry_file.geometry = geometry\n",
"\n",
"# Export to \"geometry.xml\"\n",
"geometry_file.export_to_xml()"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"With the geometry and materials finished, we now just need to define simulation parameters. In this case, we will use 5 inactive batches and 15 minimum active batches each with 2500 particles. We also tell OpenMC to turn tally triggers on, which means it will keep running until some criterion on the uncertainty of tallies is reached."
]
},
{
"cell_type": "code",
"execution_count": 11,
"metadata": {
"collapsed": true
},
"outputs": [],
"source": [
"# OpenMC simulation parameters\n",
"min_batches = 20\n",
"max_batches = 200\n",
"inactive = 5\n",
"particles = 2500\n",
"\n",
"# Instantiate a SettingsFile\n",
"settings_file = openmc.SettingsFile()\n",
"settings_file.batches = min_batches\n",
"settings_file.inactive = inactive\n",
"settings_file.particles = particles\n",
"settings_file.output = {'tallies': False}\n",
"settings_file.trigger_active = True\n",
"settings_file.trigger_max_batches = max_batches\n",
"source_bounds = [-10.71, -10.71, -10, 10.71, 10.71, 10.]\n",
"settings_file.source = Source(space=Box(\n",
" source_bounds[:3], source_bounds[3:]))\n",
"\n",
"# Export to \"settings.xml\"\n",
"settings_file.export_to_xml()"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"Let us also create a plot file that we can use to verify that our pin cell geometry was created successfully."
]
},
{
"cell_type": "code",
"execution_count": 12,
"metadata": {
"collapsed": false
},
"outputs": [],
"source": [
"# Instantiate a Plot\n",
"plot = openmc.Plot(plot_id=1)\n",
"plot.filename = 'materials-xy'\n",
"plot.origin = [0, 0, 0]\n",
"plot.width = [21.5, 21.5]\n",
"plot.pixels = [250, 250]\n",
"plot.color = 'mat'\n",
"\n",
"# Instantiate a PlotsFile, add Plot, and export to \"plots.xml\"\n",
"plot_file = openmc.PlotsFile()\n",
"plot_file.add_plot(plot)\n",
"plot_file.export_to_xml()"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"With the plots.xml file, we can now generate and view the plot. OpenMC outputs plots in .ppm format, which can be converted into a compressed format like .png with the convert utility."
]
},
{
"cell_type": "code",
"execution_count": 13,
"metadata": {
"collapsed": false
},
"outputs": [
{
"data": {
"text/plain": [
"0"
]
},
"execution_count": 13,
"metadata": {},
"output_type": "execute_result"
}
],
"source": [
"# Run openmc in plotting mode\n",
"executor = openmc.Executor()\n",
"executor.plot_geometry(output=False)"
]
},
{
"cell_type": "code",
"execution_count": 14,
"metadata": {
"collapsed": false
},
"outputs": [
{
"data": {
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAPoAAAD6AgMAAAD1grKuAAAABGdBTUEAALGPC/xhBQAAACBjSFJN\nAAB6JgAAgIQAAPoAAACA6AAAdTAAAOpgAAA6mAAAF3CculE8AAAADFBMVEX///9yEhLpgJFNv8Tq\nQYT7AAAAAWJLR0QAiAUdSAAAAAd0SU1FB+ADFxItHQxw5fwAAAPZSURBVGje7Zs7buMwEIZ9iey5\n0gyNjQpXKTYudIScgkdQYTfut1idwkdQkQNsYQO2Qj0sPiVK+mlQDmwgwIcgg8Cc4fCTSK5W4OeF\nkM8rHv+2I/rgxPZEPZgR7XtQxKdXYuUXJSUnBQ/9WCgo4vOSJ+WFUvF7E08mlia+rn7VcKXP8sRs\nzFX8b2MdX2y6v1Tw6MZUw4H4ojfIjD8mvn/qRL5p4+vvlMqvp2EhR8WBzfiz20hXORmP9fi/bM9E\neUFvV5H/0yRkeSbiGRfFJErxD9ENdz7Mbhig/h89fvtFdMiI/ePUIXV4lXju8K3DKv9NThOZ3q2K\nmUy6grxFES8rjeyic+FFQav+ncg3fXjH+Ts+/iibztFqOiZuZP/Z3OafPX40NGgST2r+uvQkXXp6\ncKvmr+r0e1Eef5um3+JHP3IFF1D/seNZJgaDmvY0Gav1s+2f1fqpIcublfKGt6apotG/NVx3SInW\ntLX+7Vg/Pv1YqOsnun6JSVdOXT/X7vk75f938QP+8OmSBs0fXtymMhJbf8qlPynYmpKCh7OB1fzN\nalOj1sl0ZAruHLiA+RM73pDe/VjMVP89+aTXwjyc/x5n+u991895/utrJTy8/06TXh0r/5JOa2Jm\nYmqi4r/vUm/H4wLmT+z4anhr05X+q6KUXhtzr/9qSff5L5uMT//V/NdU4YuBTPa/8P67l/6r44ds\n+hYuoP5jx9ciy6XTWlibBrmx8V/TdMfjkP+6pOsu/lvM9N90sf7r+f6m/65n+S8p/itN15v0UkW3\n/+48+PRfJX6S9Joo4g+G/1qYG9KroqP/WypcuvyXPf13wH89/hHef7MB6R3Cqn55U4rv4kfH3zaS\ngQuYP7HjVf89tXrbO+hfLdr+Ozv/SP1dgtQ/Ov8C+i/3+q/Zf2D/HWi6bjT6rym9I/v/03/b+LHS\n4cTg/utTsV7/net/Afzz4f0XGX84/2j9xZ4/sePR/of2X7D/o+vPo/sv6h9B/Bfxr9j1Hz2eN/hO\n8/wfff4A848+f/1A/530/I0+/8PvH9D3H9HnT+R49P0b+v4PfP/4E/wXfP8Mvf9G37/D/ovuP8Se\nP7Hj0f0vdP8tqP9O339cyv7p3P1fdP8Z3v9G999j13/seMax8x/o+ZN7+O+E8zdP/8XOf8Hnz9Dz\nb7HnT+x49PxlCp7/BM+fOv13wvnXBfivt2lMvD8TyH/Hnb+Gz3+j589jz5/Y8ej9h4D+W7qQmf57\nefqv239n3T+C7z+h969i13/seMax+3/o/cMcu/8Y2H9n3p+J6r98pv8m4fwXuH+M3n+OO3++AX9c\nlR+4PhbRAAAAJXRFWHRkYXRlOmNyZWF0ZQAyMDE2LTAzLTIzVDE0OjQ1OjI5LTA0OjAw0+qiEQAA\nACV0RVh0ZGF0ZTptb2RpZnkAMjAxNi0wMy0yM1QxNDo0NToyOS0wNDowMKK3Gq0AAAAASUVORK5C\nYII=\n",
"text/plain": [
"<IPython.core.display.Image object>"
]
},
"execution_count": 14,
"metadata": {},
"output_type": "execute_result"
}
],
"source": [
"# Convert OpenMC's funky ppm to png\n",
"!convert materials-xy.ppm materials-xy.png\n",
"\n",
"# Display the materials plot inline\n",
"Image(filename='materials-xy.png')"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"As we can see from the plot, we have a nice array of pin cells with fuel, cladding, and water! Before we run our simulation, we need to tell the code what we want to tally. The following code shows how to create a variety of tallies."
]
},
{
"cell_type": "code",
"execution_count": 15,
"metadata": {
"collapsed": true
},
"outputs": [],
"source": [
"# Instantiate an empty TalliesFile\n",
"tallies_file = openmc.TalliesFile()\n",
"tallies_file._tallies = []"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"Instantiate a fission rate mesh Tally"
]
},
{
"cell_type": "code",
"execution_count": 16,
"metadata": {
"collapsed": false
},
"outputs": [],
"source": [
"# Instantiate a tally Mesh\n",
"mesh = openmc.Mesh(mesh_id=1)\n",
"mesh.type = 'regular'\n",
"mesh.dimension = [17, 17]\n",
"mesh.lower_left = [-10.71, -10.71]\n",
"mesh.width = [1.26, 1.26]\n",
"\n",
"# Instantiate tally Filter\n",
"mesh_filter = openmc.Filter()\n",
"mesh_filter.mesh = mesh\n",
"\n",
"# Instantiate energy Filter\n",
"energy_filter = openmc.Filter()\n",
"energy_filter.type = 'energy'\n",
"energy_filter.bins = np.array([0, 0.625e-6, 20.])\n",
"\n",
"# Instantiate the Tally\n",
"tally = openmc.Tally(name='mesh tally')\n",
"tally.filters = [mesh_filter, energy_filter]\n",
"tally.scores = ['fission', 'nu-fission']\n",
"\n",
"# Add mesh and Tally to TalliesFile\n",
"tallies_file.add_mesh(mesh)\n",
"tallies_file.add_tally(tally)"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"Instantiate a cell Tally with nuclides"
]
},
{
"cell_type": "code",
"execution_count": 17,
"metadata": {
"collapsed": false
},
"outputs": [],
"source": [
"# Instantiate tally Filter\n",
"cell_filter = openmc.Filter(type='cell', bins=[fuel_cell.id])\n",
"\n",
"# Instantiate the tally\n",
"tally = openmc.Tally(name='cell tally')\n",
"tally.filters = [cell_filter]\n",
"tally.scores = ['scatter-y2']\n",
"tally.nuclides = [u235, u238]\n",
"\n",
"# Add mesh and tally to TalliesFile\n",
"tallies_file.add_tally(tally)"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"Create a \"distribcell\" Tally. The distribcell filter allows us to tally multiple repeated instances of the same cell throughout the geometry."
]
},
{
"cell_type": "code",
"execution_count": 18,
"metadata": {
"collapsed": false
},
"outputs": [],
"source": [
"# Instantiate tally Filter\n",
"distribcell_filter = openmc.Filter(type='distribcell', bins=[moderator_cell.id])\n",
"\n",
"# Instantiate tally Trigger for kicks\n",
"trigger = openmc.Trigger(trigger_type='std_dev', threshold=5e-5)\n",
"trigger.scores = ['absorption']\n",
"\n",
"# Instantiate the Tally\n",
"tally = openmc.Tally(name='distribcell tally')\n",
"tally.filters = [distribcell_filter]\n",
"tally.scores = ['absorption', 'scatter']\n",
"tally.triggers = [trigger]\n",
"\n",
"# Add mesh and tally to TalliesFile\n",
"tallies_file.add_tally(tally)"
]
},
{
"cell_type": "code",
"execution_count": 19,
"metadata": {
"collapsed": true
},
"outputs": [],
"source": [
"# Export to \"tallies.xml\"\n",
"tallies_file.export_to_xml()"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"Now we a have a complete set of inputs, so we can go ahead and run our simulation."
]
},
{
"cell_type": "code",
"execution_count": 20,
"metadata": {
"collapsed": false
},
"outputs": [
{
"name": "stdout",
"output_type": "stream",
"text": [
"\n",
" .d88888b. 888b d888 .d8888b.\n",
" d88P\" \"Y88b 8888b d8888 d88P Y88b\n",
" 888 888 88888b.d88888 888 888\n",
" 888 888 88888b. .d88b. 88888b. 888Y88888P888 888 \n",
" 888 888 888 \"88b d8P Y8b 888 \"88b 888 Y888P 888 888 \n",
" 888 888 888 888 88888888 888 888 888 Y8P 888 888 888\n",
" Y88b. .d88P 888 d88P Y8b. 888 888 888 \" 888 Y88b d88P\n",
" \"Y88888P\" 88888P\" \"Y8888 888 888 888 888 \"Y8888P\"\n",
"__________________888______________________________________________________\n",
" 888\n",
" 888\n",
"\n",
" Copyright: 2011-2015 Massachusetts Institute of Technology\n",
" License: http://mit-crpg.github.io/openmc/license.html\n",
" Version: 0.7.1\n",
" Git SHA1: 5f252e2df51930b9175fd41bafa8db01f3eaeb92\n",
" Date/Time: 2016-03-23 14:45:30\n",
" MPI Processes: 1\n",
" OpenMP Threads: 16\n",
"\n",
" ===========================================================================\n",
" ========================> INITIALIZATION <=========================\n",
" ===========================================================================\n",
"\n",
" Reading settings XML file...\n",
" Reading cross sections XML file...\n",
" Reading geometry XML file...\n",
" Reading materials XML file...\n",
" Reading tallies XML file...\n",
" Building neighboring cells lists for each surface...\n",
" Loading ACE cross section table: 92235.71c\n",
" Loading ACE cross section table: 92238.71c\n",
" Loading ACE cross section table: 8016.71c\n",
" Loading ACE cross section table: 1001.71c\n",
" Loading ACE cross section table: 5010.71c\n",
" Loading ACE cross section table: 40090.71c\n",
" Maximum neutron transport energy: 20.0000 MeV for 92235.71c\n",
" Initializing source particles...\n",
"\n",
" ===========================================================================\n",
" ====================> K EIGENVALUE SIMULATION <====================\n",
" ===========================================================================\n",
"\n",
" Bat./Gen. k Average k \n",
" ========= ======== ==================== \n",
" 1/1 0.51036 \n",
" 2/1 0.64436 \n",
" 3/1 0.64874 \n",
" 4/1 0.65998 \n",
" 5/1 0.68369 \n",
" 6/1 0.69058 \n",
" 7/1 0.68288 0.68673 +/- 0.00385\n",
" 8/1 0.69483 0.68943 +/- 0.00350\n",
" 9/1 0.70348 0.69294 +/- 0.00430\n",
" 10/1 0.69969 0.69429 +/- 0.00359\n",
" 11/1 0.67170 0.69052 +/- 0.00477\n",
" 12/1 0.67661 0.68854 +/- 0.00450\n",
" 13/1 0.69571 0.68943 +/- 0.00400\n",
" 14/1 0.67433 0.68776 +/- 0.00390\n",
" 15/1 0.67744 0.68672 +/- 0.00364\n",
" 16/1 0.65256 0.68362 +/- 0.00453\n",
" 17/1 0.66657 0.68220 +/- 0.00437\n",
" 18/1 0.66887 0.68117 +/- 0.00415\n",
" 19/1 0.68238 0.68126 +/- 0.00384\n",
" 20/1 0.64423 0.67879 +/- 0.00435\n",
" Triggers unsatisfied, max unc./thresh. is 1.40549 for absorption in tally 10002\n",
" The estimated number of batches is 35\n",
" Creating state point statepoint.020.h5...\n",
" 21/1 0.66266 0.67778 +/- 0.00419\n",
" 22/1 0.67656 0.67771 +/- 0.00393\n",
" 23/1 0.67643 0.67764 +/- 0.00371\n",
" 24/1 0.66192 0.67681 +/- 0.00361\n",
" 25/1 0.69848 0.67789 +/- 0.00359\n",
" 26/1 0.66274 0.67717 +/- 0.00349\n",
" 27/1 0.69746 0.67810 +/- 0.00345\n",
" 28/1 0.67485 0.67795 +/- 0.00330\n",
" 29/1 0.67427 0.67780 +/- 0.00316\n",
" 30/1 0.66531 0.67730 +/- 0.00308\n",
" 31/1 0.68457 0.67758 +/- 0.00297\n",
" 32/1 0.66592 0.67715 +/- 0.00289\n",
" 33/1 0.65929 0.67651 +/- 0.00286\n",
" 34/1 0.67252 0.67637 +/- 0.00276\n",
" 35/1 0.71827 0.67777 +/- 0.00301\n",
" Triggers satisfied for batch 35\n",
" Creating state point statepoint.035.h5...\n",
"\n",
" ===========================================================================\n",
" ======================> SIMULATION FINISHED <======================\n",
" ===========================================================================\n",
"\n",
"\n",
" =======================> TIMING STATISTICS <=======================\n",
"\n",
" Total time for initialization = 4.5300E-01 seconds\n",
" Reading cross sections = 1.3200E-01 seconds\n",
" Total time in simulation = 1.9780E+00 seconds\n",
" Time in transport only = 1.7780E+00 seconds\n",
" Time in inactive batches = 2.1000E-01 seconds\n",
" Time in active batches = 1.7680E+00 seconds\n",
" Time synchronizing fission bank = 5.0000E-03 seconds\n",
" Sampling source sites = 5.0000E-03 seconds\n",
" SEND/RECV source sites = 0.0000E+00 seconds\n",
" Time accumulating tallies = 1.0000E-03 seconds\n",
" Total time for finalization = 0.0000E+00 seconds\n",
" Total time elapsed = 2.4480E+00 seconds\n",
" Calculation Rate (inactive) = 59523.8 neutrons/second\n",
" Calculation Rate (active) = 21210.4 neutrons/second\n",
"\n",
" ============================> RESULTS <============================\n",
"\n",
" k-effective (Collision) = 0.67866 +/- 0.00337\n",
" k-effective (Track-length) = 0.67777 +/- 0.00301\n",
" k-effective (Absorption) = 0.68234 +/- 0.00332\n",
" Combined k-effective = 0.67987 +/- 0.00255\n",
" Leakage Fraction = 0.34141 +/- 0.00198\n",
"\n"
]
},
{
"data": {
"text/plain": [
"0"
]
},
"execution_count": 20,
"metadata": {},
"output_type": "execute_result"
}
],
"source": [
"# Remove old HDF5 (summary, statepoint) files\n",
"!rm statepoint.*\n",
"\n",
"# Run OpenMC with MPI!\n",
"executor.run_simulation()"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"## Tally Data Processing"
]
},
{
"cell_type": "code",
"execution_count": 21,
"metadata": {
"collapsed": false
},
"outputs": [],
"source": [
"# We do not know how many batches were needed to satisfy the \n",
"# tally trigger(s), so find the statepoint file(s)\n",
"statepoints = glob.glob('statepoint.*.h5')\n",
"\n",
"# Load the last statepoint file\n",
"sp = StatePoint(statepoints[-1])"
]
},
{
"cell_type": "code",
"execution_count": 22,
"metadata": {
"collapsed": false,
"scrolled": true
},
"outputs": [],
"source": [
"# Load the summary file and link with statepoint\n",
"su = Summary('summary.h5')\n",
"sp.link_with_summary(su)"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"**Analyze the mesh fission rate tally**"
]
},
{
"cell_type": "code",
"execution_count": 23,
"metadata": {
"collapsed": false
},
"outputs": [
{
"name": "stdout",
"output_type": "stream",
"text": [
"Tally\n",
"\tID =\t10000\n",
"\tName =\tmesh tally\n",
"\tFilters =\t\n",
" \t\tmesh\t[1]\n",
" \t\tenergy\t[ 0.00000000e+00 6.25000000e-07 2.00000000e+01]\n",
"\tNuclides =\ttotal \n",
"\tScores =\t[u'fission', u'nu-fission']\n",
"\tEstimator =\ttracklength\n",
"\n"
]
}
],
"source": [
"# Find the mesh tally with the StatePoint API\n",
"tally = sp.get_tally(name='mesh tally')\n",
"\n",
"# Print a little info about the mesh tally to the screen\n",
"print(tally)"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"Use the new Tally data retrieval API with pure NumPy"
]
},
{
"cell_type": "code",
"execution_count": 24,
"metadata": {
"collapsed": false
},
"outputs": [
{
"name": "stdout",
"output_type": "stream",
"text": [
"[[[ 0.12916959]]\n",
"\n",
" [[ 0.06336943]]\n",
"\n",
" [[ 0.33288738]]\n",
"\n",
" [[ 0.14666158]]]\n"
]
}
],
"source": [
"# Get the relative error for the thermal fission reaction \n",
"# rates in the four corner pins \n",
"data = tally.get_values(scores=['fission'], filters=['mesh', 'energy'], \\\n",
" filter_bins=[((1,1),(1,17), (17,1), (17,17)), \\\n",
" ((0., 0.625e-6),)], value='rel_err')\n",
"print(data)"
]
},
{
"cell_type": "code",
"execution_count": 25,
"metadata": {
"collapsed": false
},
"outputs": [
{
"data": {
"text/html": [
"<div>\n",
"<table border=\"1\" class=\"dataframe\">\n",
" <thead>\n",
" <tr>\n",
" <th></th>\n",
" <th colspan=\"3\" halign=\"left\">mesh 1</th>\n",
" <th>energy low [MeV]</th>\n",
" <th>energy high [MeV]</th>\n",
" <th>score</th>\n",
" <th>mean</th>\n",
" <th>std. dev.</th>\n",
" </tr>\n",
" <tr>\n",
" <th></th>\n",
" <th>x</th>\n",
" <th>y</th>\n",
" <th>z</th>\n",
" <th></th>\n",
" <th></th>\n",
" <th></th>\n",
" <th></th>\n",
" <th></th>\n",
" </tr>\n",
" </thead>\n",
" <tbody>\n",
" <tr>\n",
" <th>0</th>\n",
" <td>1</td>\n",
" <td>1</td>\n",
" <td>1</td>\n",
" <td>0.00e+00</td>\n",
" <td>6.25e-07</td>\n",
" <td>fission</td>\n",
" <td>2.37e-04</td>\n",
" <td>3.06e-05</td>\n",
" </tr>\n",
" <tr>\n",
" <th>1</th>\n",
" <td>1</td>\n",
" <td>1</td>\n",
" <td>1</td>\n",
" <td>0.00e+00</td>\n",
" <td>6.25e-07</td>\n",
" <td>nu-fission</td>\n",
" <td>5.78e-04</td>\n",
" <td>7.46e-05</td>\n",
" </tr>\n",
" <tr>\n",
" <th>2</th>\n",
" <td>1</td>\n",
" <td>1</td>\n",
" <td>1</td>\n",
" <td>6.25e-07</td>\n",
" <td>2.00e+01</td>\n",
" <td>fission</td>\n",
" <td>7.00e-05</td>\n",
" <td>5.15e-06</td>\n",
" </tr>\n",
" <tr>\n",
" <th>3</th>\n",
" <td>1</td>\n",
" <td>1</td>\n",
" <td>1</td>\n",
" <td>6.25e-07</td>\n",
" <td>2.00e+01</td>\n",
" <td>nu-fission</td>\n",
" <td>1.85e-04</td>\n",
" <td>1.28e-05</td>\n",
" </tr>\n",
" <tr>\n",
" <th>4</th>\n",
" <td>1</td>\n",
" <td>2</td>\n",
" <td>1</td>\n",
" <td>0.00e+00</td>\n",
" <td>6.25e-07</td>\n",
" <td>fission</td>\n",
" <td>4.04e-04</td>\n",
" <td>3.09e-05</td>\n",
" </tr>\n",
" <tr>\n",
" <th>5</th>\n",
" <td>1</td>\n",
" <td>2</td>\n",
" <td>1</td>\n",
" <td>0.00e+00</td>\n",
" <td>6.25e-07</td>\n",
" <td>nu-fission</td>\n",
" <td>9.85e-04</td>\n",
" <td>7.54e-05</td>\n",
" </tr>\n",
" <tr>\n",
" <th>6</th>\n",
" <td>1</td>\n",
" <td>2</td>\n",
" <td>1</td>\n",
" <td>6.25e-07</td>\n",
" <td>2.00e+01</td>\n",
" <td>fission</td>\n",
" <td>1.00e-04</td>\n",
" <td>5.08e-06</td>\n",
" </tr>\n",
" <tr>\n",
" <th>7</th>\n",
" <td>1</td>\n",
" <td>2</td>\n",
" <td>1</td>\n",
" <td>6.25e-07</td>\n",
" <td>2.00e+01</td>\n",
" <td>nu-fission</td>\n",
" <td>2.63e-04</td>\n",
" <td>1.34e-05</td>\n",
" </tr>\n",
" <tr>\n",
" <th>8</th>\n",
" <td>1</td>\n",
" <td>3</td>\n",
" <td>1</td>\n",
" <td>0.00e+00</td>\n",
" <td>6.25e-07</td>\n",
" <td>fission</td>\n",
" <td>5.82e-04</td>\n",
" <td>5.00e-05</td>\n",
" </tr>\n",
" <tr>\n",
" <th>9</th>\n",
" <td>1</td>\n",
" <td>3</td>\n",
" <td>1</td>\n",
" <td>0.00e+00</td>\n",
" <td>6.25e-07</td>\n",
" <td>nu-fission</td>\n",
" <td>1.42e-03</td>\n",
" <td>1.22e-04</td>\n",
" </tr>\n",
" <tr>\n",
" <th>10</th>\n",
" <td>1</td>\n",
" <td>3</td>\n",
" <td>1</td>\n",
" <td>6.25e-07</td>\n",
" <td>2.00e+01</td>\n",
" <td>fission</td>\n",
" <td>1.38e-04</td>\n",
" <td>1.03e-05</td>\n",
" </tr>\n",
" <tr>\n",
" <th>11</th>\n",
" <td>1</td>\n",
" <td>3</td>\n",
" <td>1</td>\n",
" <td>6.25e-07</td>\n",
" <td>2.00e+01</td>\n",
" <td>nu-fission</td>\n",
" <td>3.59e-04</td>\n",
" <td>2.54e-05</td>\n",
" </tr>\n",
" <tr>\n",
" <th>12</th>\n",
" <td>1</td>\n",
" <td>4</td>\n",
" <td>1</td>\n",
" <td>0.00e+00</td>\n",
" <td>6.25e-07</td>\n",
" <td>fission</td>\n",
" <td>6.88e-04</td>\n",
" <td>4.25e-05</td>\n",
" </tr>\n",
" <tr>\n",
" <th>13</th>\n",
" <td>1</td>\n",
" <td>4</td>\n",
" <td>1</td>\n",
" <td>0.00e+00</td>\n",
" <td>6.25e-07</td>\n",
" <td>nu-fission</td>\n",
" <td>1.68e-03</td>\n",
" <td>1.04e-04</td>\n",
" </tr>\n",
" <tr>\n",
" <th>14</th>\n",
" <td>1</td>\n",
" <td>4</td>\n",
" <td>1</td>\n",
" <td>6.25e-07</td>\n",
" <td>2.00e+01</td>\n",
" <td>fission</td>\n",
" <td>1.62e-04</td>\n",
" <td>7.43e-06</td>\n",
" </tr>\n",
" <tr>\n",
" <th>15</th>\n",
" <td>1</td>\n",
" <td>4</td>\n",
" <td>1</td>\n",
" <td>6.25e-07</td>\n",
" <td>2.00e+01</td>\n",
" <td>nu-fission</td>\n",
" <td>4.22e-04</td>\n",
" <td>1.93e-05</td>\n",
" </tr>\n",
" <tr>\n",
" <th>16</th>\n",
" <td>1</td>\n",
" <td>5</td>\n",
" <td>1</td>\n",
" <td>0.00e+00</td>\n",
" <td>6.25e-07</td>\n",
" <td>fission</td>\n",
" <td>7.62e-04</td>\n",
" <td>5.69e-05</td>\n",
" </tr>\n",
" <tr>\n",
" <th>17</th>\n",
" <td>1</td>\n",
" <td>5</td>\n",
" <td>1</td>\n",
" <td>0.00e+00</td>\n",
" <td>6.25e-07</td>\n",
" <td>nu-fission</td>\n",
" <td>1.86e-03</td>\n",
" <td>1.39e-04</td>\n",
" </tr>\n",
" <tr>\n",
" <th>18</th>\n",
" <td>1</td>\n",
" <td>5</td>\n",
" <td>1</td>\n",
" <td>6.25e-07</td>\n",
" <td>2.00e+01</td>\n",
" <td>fission</td>\n",
" <td>1.80e-04</td>\n",
" <td>8.16e-06</td>\n",
" </tr>\n",
" <tr>\n",
" <th>19</th>\n",
" <td>1</td>\n",
" <td>5</td>\n",
" <td>1</td>\n",
" <td>6.25e-07</td>\n",
" <td>2.00e+01</td>\n",
" <td>nu-fission</td>\n",
" <td>4.71e-04</td>\n",
" <td>2.08e-05</td>\n",
" </tr>\n",
" </tbody>\n",
"</table>\n",
"</div>"
],
"text/plain": [
" mesh 1 energy low [MeV] energy high [MeV] score mean \\\n",
" x y z \n",
"0 1 1 1 0.00e+00 6.25e-07 fission 2.37e-04 \n",
"1 1 1 1 0.00e+00 6.25e-07 nu-fission 5.78e-04 \n",
"2 1 1 1 6.25e-07 2.00e+01 fission 7.00e-05 \n",
"3 1 1 1 6.25e-07 2.00e+01 nu-fission 1.85e-04 \n",
"4 1 2 1 0.00e+00 6.25e-07 fission 4.04e-04 \n",
"5 1 2 1 0.00e+00 6.25e-07 nu-fission 9.85e-04 \n",
"6 1 2 1 6.25e-07 2.00e+01 fission 1.00e-04 \n",
"7 1 2 1 6.25e-07 2.00e+01 nu-fission 2.63e-04 \n",
"8 1 3 1 0.00e+00 6.25e-07 fission 5.82e-04 \n",
"9 1 3 1 0.00e+00 6.25e-07 nu-fission 1.42e-03 \n",
"10 1 3 1 6.25e-07 2.00e+01 fission 1.38e-04 \n",
"11 1 3 1 6.25e-07 2.00e+01 nu-fission 3.59e-04 \n",
"12 1 4 1 0.00e+00 6.25e-07 fission 6.88e-04 \n",
"13 1 4 1 0.00e+00 6.25e-07 nu-fission 1.68e-03 \n",
"14 1 4 1 6.25e-07 2.00e+01 fission 1.62e-04 \n",
"15 1 4 1 6.25e-07 2.00e+01 nu-fission 4.22e-04 \n",
"16 1 5 1 0.00e+00 6.25e-07 fission 7.62e-04 \n",
"17 1 5 1 0.00e+00 6.25e-07 nu-fission 1.86e-03 \n",
"18 1 5 1 6.25e-07 2.00e+01 fission 1.80e-04 \n",
"19 1 5 1 6.25e-07 2.00e+01 nu-fission 4.71e-04 \n",
"\n",
" std. dev. \n",
" \n",
"0 3.06e-05 \n",
"1 7.46e-05 \n",
"2 5.15e-06 \n",
"3 1.28e-05 \n",
"4 3.09e-05 \n",
"5 7.54e-05 \n",
"6 5.08e-06 \n",
"7 1.34e-05 \n",
"8 5.00e-05 \n",
"9 1.22e-04 \n",
"10 1.03e-05 \n",
"11 2.54e-05 \n",
"12 4.25e-05 \n",
"13 1.04e-04 \n",
"14 7.43e-06 \n",
"15 1.93e-05 \n",
"16 5.69e-05 \n",
"17 1.39e-04 \n",
"18 8.16e-06 \n",
"19 2.08e-05 "
]
},
"execution_count": 25,
"metadata": {},
"output_type": "execute_result"
}
],
"source": [
"# Get a pandas dataframe for the mesh tally data\n",
"df = tally.get_pandas_dataframe(nuclides=False)\n",
"\n",
"# Set the Pandas float display settings\n",
"import pandas as pd\n",
"pd.set_option('display.float_format', '{:.2e}'.format)\n",
"\n",
"# Print the first twenty rows in the dataframe\n",
"df.head(20)"
]
},
{
"cell_type": "code",
"execution_count": 26,
"metadata": {
"collapsed": false
},
"outputs": [
{
"data": {
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAY8AAAEaCAYAAADpMdsXAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAIABJREFUeJzt3X+c3FV97/HXO8GUViRLRIMmhC0h/BRdow15XBG3coUk\n3BLqD2qslY23vWlJtDwqXEStpLS9/qo1AqVITZW0YEQpPy5GCNQMF60klLABLkvYIBtDiFBJUiH0\nNiH53D++Z5dxdnbmO5udmZ3d9/PxGPme7/ec75xZv5nPnHO+53wVEZiZmdViQrMrYGZmrcfBw8zM\naubgYWZmNXPwMDOzmjl4mJlZzRw8zMysZg4e1lCS9kvaKKlb0r9KmluH93ihyvFjJC0a6fetN0kX\nSLqqzP7LJf1JM+pk45eDhzXanoiYHREdwKeAz9fhPapNXvp14EMH8waSmvVvZ1RPzJI0sdl1sMZw\n8LBGU9H2ZGDnwAHpS5IekbRJ0vlp33mS7knbb5C0WdLr06/wWyWtS/s+W/bNfvmcH0i7PwecnlpA\nf1ySX5KukfSYpLskfU/Se9OxpyR9XtK/Au+X9BZJP06tqJslTU751kmanbZfK+mptD1knSX9rqT1\nqU5/K0lp/+KU937gHRX+rh2S/iXl/e+p7PWSzi16j3+U9Fsln/coSfem931Y0jvS/nmSHpT0kKS7\n074jJN2S/pb/IulNaf/lklZJ+iGwStIESV9Mn6db0h9UqLe1qojwy6+GvYCXgY1AD7ALeGva/17g\nrrT9emArMDWlVwFLgf8NnJ/2XQBsB9qAQ4FHgNnp2C/Sf99X7pzAu4Dbh6jf+4A70vZUsuD23pR+\nCri4KO8m4PS0/WfAX6ftdUV1eS3wk0p1Bk4Ebgcmpnx/A3wYOCrVeQpwCPBD4Moydb4ceAiYlN7v\np6nsGcAtKc/hwJPAhJKyfwJclrYFvBo4Mp1jRtrflv57JfCnafs3gYeK3v8BYFJK/wHwqbQ9KR07\nptnXnl8j+zoEs8Z6KSL6f5XPBf4BeBNwOvAtgIh4TlIB+A3gDuDjwKPAjyPipqJz3R0Ru9O5/imd\nY2PR8XcMcc5KYyKnA99JZZ6VtK7k+LfT+x0OTI6IH6b91wM3UV1xnW9O77cfeBvwQGpxHAo8C5wG\nrIuInSn/t4FZQ5z3tojYCzwv6QfAnIi4XdLfSHot8H7g5og4UFLuAWClpFelc2yS9JvAvRHx0/R3\n2F30t3lv2rdO0hRJh6Vjt6f3BzgLOLWopXd4qvfWHH8faxEOHtY0EXG/pCMlHVnmcHH31tHAAbKW\nwC+dokq60jmHa0+OPC/zSpfwoSXHiuuoovQ3I+LTxRklLSR/nYc67yrg94APAl2DCkXcJ+kM4Bzg\nG5L+Gtg9xPtW+vsW/10EfCwi7s5Zd2tBHvOwRhv4UpJ0Itk1+DxwH/A7qb/8dcA7gQ2SDgFWkn35\n9Uj6RNG53iOpTdKvAueRdesUv0fZc5K1PF4zRP1+BLwvjX1MBTrLZYqIXwC7+scIyL6g703bfcDb\n0/YHSoqW1vlHwA/IxlBel/4uR0iaAawHzkjpV5U5V7GFkialVsa7yFoUkLWILsqqHI+XFkrv81xE\nrCT7O88G7gfeKemY/vqk7PeRdachqRP4eUS8WKYudwEXpv/vkDQrfV4bQ9zysEY7VNJGXvmC/0hE\nBHBL6sbaRNbKuCR1Nf0p8H8i4l8kPUwWUO5IZTcA/wRMA/4hIh5K+wMgIoY6507ggKSHyH7xf7Wo\nfjcD7wb+L7ANeBD49+LzFrkA+Fr6YvwJsDjt/yvgpjRQ/L2SMqV13ggg6TPAWmV3ce0FlkbEBknL\nyb7MdwHdFf6uDwMFsjGPKyLiZ+lv8JykHuCWIcp1ApdI2kcWVD8SET+X9D/I/j8R8BxwNtm4zt9L\n2kTW0vjIEOf8OtAObCwqf16FulsLUvbv1qy1SLoAeFtEfLwO5351ROyRNIXs1/87IuK5EThv3epc\n4T1/jSx4zo6IivNfzGrhlofZYHdIagNeRfYr/qADRzNIOpOsK+rLDhw20tzyMDOzmnnA3CwnZZME\nL06T5F6Q9HfKJiyukfQLSWv1ykTBuZJ+JGlXmmj3rqLzdCmbhPgLSVvS+EL/sXdJ2ibpTyQ9K2m7\npK4mfFyzihw8zGrzXuBM4HjgXGAN8EmyiXUTgY9LeiPZ/JQrIuII4GLg5nQnFGRzOBZExOFkg+xf\nkdRR9B5Hkd0N9kbg94G/6Q9KZqOFg4dZba6KiJ9HxA6yW1fXR8TDaYLcLWS3un4Y+F5E3AUQEf8M\n/CuwIKW/HxF9afs+YC3ZbcT99gJ/HhH7I+L7wIvACQ35dGY5OXiY1ebZou3/KJM+DDgGOF/SzvTa\nRTbb/Q0AkuYrWxPr+XRsPlnLpd/zJTPBX0rnNRs1fLeV2cgKsnWhVkXEktKDkiYB3yVrndwWEQck\n3cLIzH43axi3PMxG3j8C50o6K81uPzQNhL+RbKHASWSzsw9Imk+2FpRZS3HwMMsv11paEbGdbDD9\nU8C/kS0IeDHZirYvki30+J000/2DwG01vq9Z0+Wa5yFpHrCCLNisjIgvlMlzJVnf7R6gKyK685RN\naxV9CTgyInam9XR6gP51eO6PiAuH+fnMzKwOqo55pLV2ria7PfEZsmWjbyteZC01vWdGxCxJpwHX\nAnOrlZU0HXgPg5dq3tK/bLeZmY0+ebqt5gC9EbE1IvYBq4GFJXkWki39TESsByanFUmrlf0KcEmZ\n9/TgoZnZKJYneEwjW12039NpX548Q5ZV9njMbRHxSJn3bFf2WMx1kk7PUUczM2uget2qW7HlkJaw\n/hRZl1VpmWfIHn+5S9lzoG+VdPIQzw0wM7MmyBM8tgMzitLT077SPEeXyTNpiLIzydb735TW+58O\nPChpTlrBdBdARGyU9CTZUhDFjxdFku9AMTOrs4go2xjIEzweAI5Ld0HtILu1cFFJntuBpcC308N3\ndqfnP/+8XNmI6CFbvwfIFpwje97ALmWPJN2Z7oE/FjiO7EE75T5UjupbrZYvX87y5cubXQ2z3HzN\n1kf22768qsEjIvZLWka2/k7/7bY9kpZkh+O6iFgjaYGkLWS36i6uVLbc2/BKt9UZwBWS9pI9/W1J\nROzO+2HNzKz+co15RMSdlCzMFhFfK0kvy1u2TJ5ji7b/iewxndYkfX19za6CWU18zTaeZ5jbIB0d\nHdUzmY0ivmYbr2WfJCgpWrXuZmatQNKQA+ZueZiZWc0cPGyQQqHQ7CqY1cTXbOM5eJiZWc085mFm\nZmV5zMPMzEaUg4cN4v5jazW+ZhvPwcPMzGrmMQ8zMyvLYx5mZjaiHDxsEPcfW6vxNdt4Dh5m1vK6\nu7ubXYVxx8HDBuns7Gx2Fcxqsnu3n9rQaA4eZtby7r///mZXYdzJFTwkzZP0uKQnJF06RJ4rJfVK\n6pbUkbespE9IOiBpStG+y9K5eiSdNZwPZsPn/mNrBYVCYeAJgnfdddfAtq/fxqh6q66kCcATwJnA\nM2SPpf1gRDxelGc+sCwizpF0GvDViJhbrayk6cDXyR4W9baI2CnpJOBG4DfInm1+DzCr9L5c36pb\nP4VCwV1X1lKmTJnCzp07m12NMedgb9WdA/RGxNaI2AesBhaW5FkIrAKIiPXAZElTc5T9CnBJmXOt\njoiXI6IP6E3nsQZx4LBWsGLFCjo7O+ns7GTXrl0D2ytWrGh21caFPI+hnQZsK0o/zeAv83J5plUq\nK+lcYFtEPFLykPVpwI+L0tvTPjOzARdddBEXXXQRAJMmTXJ3VYPleob5MJRt5gwclH4V+BTwnjq9\nvx0Ed1tZKygUCgMBY9++fSxfvhxgoAVi9ZUneGwHZhSlp6d9pXmOLpNn0hBlZwLtwCZlzY7pwEZJ\nc3K+HwBdXV20t7cD0NbWRkdHx8BF039ROe2002Mz/d3vfpdHH32UfrfeeittbW20tbUN7BtN9W2F\ndP92X18f1eQZMJ8IbCYb9N4BbAAWRURPUZ4FwNI0YD4XWJEGzKuWTeWfAmZHxC5JJwM3AKeRdVfd\njQfMzayCCRMmcODAgWZXY8ypNGBeteUREfslLQPWkg2wr4yIHklLssNxXUSskbRA0hZgD7C4Utly\nb0Pq6oqIxyTdBDwG7AMudJQws1IrVqzg1ltvBSAiBn5Fn3feeQNjIVY/ucY8IuJOsttpi/d9rSS9\nLG/ZMnmOLUl/DvhcnrrZyCt4zMNaQEdHx8DM8nvvvXfgmu3o6KhQykZKvQbMzczqqru7e1BfPWTj\nn/7xU39+noeZtTxPEqyPgxrzMDMbDUrmg+U+7h+Z9eGFEW2Q4q4As9EiIoZ8wdIKx6weHDzMbAx4\nf7MrMO54zMPMWp4E/joYeX6GuZmNaZdf3uwajD8OHjaIxzys1XR2FppdhXHHwcPMzGrmMQ8zMyvL\nYx5mZjaiHDxsEI95WKvxNdt4Dh5m1vK++c1m12D88ZiHmbU8z/OoD495mJnZiMoVPCTNk/S4pCck\nXTpEnisl9UrqltRRraykKyRtkvSQpDslHZX2HyPpJUkb0+uag/2QVhv3H1vrKTS7AuNO1eAhaQJw\nNXA2cAqwSNKJJXnmAzMjYhawBLg2R9kvRsRbIuKtwPeA4jmiWyJidnpdeFCf0MzMRlyelsccoDci\ntkbEPmA1sLAkz0JgFUBErAcmS5paqWxEvFhU/tVA8QOIK6+9bHXlB+lY6+lsdgXGnTzBYxqwrSj9\ndNqXJ0/FspL+QtJPgQ8Bny3K1566rNZJOj1HHc1sHPPaVo1XrwHzXC2HiPhMRMwAbgA+lnbvAGZE\nxGzgE8CNkg6rTzWtHI95WKvx2laNl+dJgtuBGUXp6WlfaZ6jy+SZlKMswI3AGmB5ROwF9gJExEZJ\nTwLHAxtLC3V1ddHe3g5kzy3u6OgY6HLp/wJ02mmnx366u7t7VNWnVdP92319fVRTdZ6HpInAZuBM\nslbBBmBRRPQU5VkALI2IcyTNBVZExNxKZSUdFxFbUvmPAe+MiPMlHQnsjIgDko4F7gVOjYjdJfXy\nPA8zszo6qGeYR8R+ScuAtWTdXCvTl/+S7HBcFxFrJC2QtAXYAyyuVDad+vOSjicbKN8K/GHafwZw\nhaS96diS0sBhZmbN5RnmNkihUBhozpq1Al+z9eEZ5mY2pnltq8Zzy8PMWp7XtqoPtzzMzGxEOXjY\nIMW37Zm1hkKzKzDuOHiYmVnNPOZhZi3PYx714TEPMxvTvLZV4zl42CAe87BW47WtGs/Bw8zMauYx\nDzMzK8tjHmZmNqIcPGwQj3lYq/E123gOHmbW8ry2VeN5zMPMWp7nedSHxzzMzGxE5QoekuZJelzS\nE5IuHSLPlZJ6JXVL6qhWVtIVkjZJekjSnZKOKjp2WTpXj6SzDuYDWu3cf2ytp9DsCow7VYOHpAnA\n1cDZwCnAIkknluSZD8yMiFnAEuDaHGW/GBFviYi3At8DLk9lTgbOB04C5gPXSCrbbDIzs+bI0/KY\nA/RGxNaI2AesBhaW5FkIrAKIiPXAZElTK5WNiBeLyr+a7JGzAOcCqyPi5YjoA3rTeaxB/EQ2az2d\nza7AuJMneEwDthWln0778uSpWFbSX0j6KfAh4LNDnGt7mfczMxvgta0ar14D5rm6mSLiMxExA7gB\n+Fid6mI18piHtRqvbdV4h+TIsx2YUZSenvaV5jm6TJ5JOcoC3Eg27rG8wrkG6erqor29HYC2tjY6\nOjoGulz6vwCddtrpsZ/u7u4eVfVp1XT/dl9fH9VUnechaSKwGTgT2AFsABZFRE9RngXA0og4R9Jc\nYEVEzK1UVtJxEbEllf8Y8M6IOD8NmN8AnEbWXXU3MKt0UofneZiZ1VeleR5VWx4RsV/SMmAtWTfX\nyvTlvyQ7HNdFxBpJCyRtAfYAiyuVTaf+vKTjyQbKtwJ/mMo8Jukm4DFgH3Cho4SZ2ejiGeY2SKFQ\nGGjOmrUCX7P14RnmZjameW2rxnPLw8xante2qg+3PMzMbEQ5eNggxbftmbWGQrMrMO44eJiZWc08\n5mFmLc9jHvXhMQ8zG9O8tlXjOXjYIB7zsFbjta0az8HDzMxq5jEPMzMry2MeZmY2ohw8bBCPeVir\n8TXbeA4eZtbyvLZV43nMw8xanud51IfHPMzMbETlCh6S5kl6XNITki4dIs+VknoldUvqqFZW0hcl\n9aT8N0s6PO0/RtJLkjam1zUH+yGtNu4/ttZTaHYFxp2qwUPSBOBq4GzgFGCRpBNL8swHZkbELGAJ\ncG2OsmuBUyKiA+gFLis65ZaImJ1eFx7MBzQzs5GXp+UxB+iNiK0RsQ9YDSwsybMQWAUQEeuByZKm\nViobEfdExIFU/n5getH5yvaxWWP4iWzWejqbXYFxJ0/wmAZsK0o/nfblyZOnLMBHge8XpdtTl9U6\nSafnqKOZjWNe26rx6jVgnrvlIOnTwL6IuDHtegaYERGzgU8AN0o6rA51tCF4zMNajde2arxDcuTZ\nDswoSk9P+0rzHF0mz6RKZSV1AQuAd/fvS91bu9L2RklPAscDG0sr1tXVRXt7OwBtbW10dHQMdLn0\nfwE67bTTYz/d3d09qurTqun+7b6+PqqpOs9D0kRgM3AmsAPYACyKiJ6iPAuApRFxjqS5wIqImFup\nrKR5wJeBMyLi+aJzHQnsjIgDko4F7gVOjYjdJfXyPA8zszqqNM+jassjIvZLWkZ2d9QEYGX68l+S\nHY7rImKNpAWStgB7gMWVyqZTX0XWMrlbEsD96c6qM4ArJO0FDgBLSgOHmZk1l2eY2yCFQmGgOWvW\nCnzN1odnmJvZmOa1rRrPLQ8za3le26o+3PIwM7MR5eBhgxTftmfWGgrNrsC44+BhZmY185iHmbU8\nj3nUh8c8zGxM89pWjefgYYN4zMNajde2ajwHDzMzq5nHPMzMrCyPeZiZ2Yhy8LBBPOZhrcbXbOM5\neJhZy/PaVo3nMQ8za3me51EfHvMwM7MRlSt4SJon6XFJT0i6dIg8V0rqldQtqaNaWUlflNST8t8s\n6fCiY5elc/VIOutgPqDVzv3H1noKza7AuFM1eEiaAFwNnA2cAiySdGJJnvnAzIiYBSwBrs1Rdi1w\nSkR0AL3AZanMycD5wEnAfOAapUcNmpnZ6JCn5TEH6I2IrRGxD1gNLCzJsxBYBRAR64HJkqZWKhsR\n90TEgVT+fmB62j4XWB0RL0dEH1lgmTPcD2i18xPZrPV0NrsC406e4DEN2FaUfjrty5MnT1mAjwJr\nhjjX9iHKmJkBXtuqGeo1YJ67m0nSp4F9EfGtOtXFauQxD2s1Xtuq8Q7JkWc7MKMoPT3tK81zdJk8\nkyqVldQFLADeneNcg3R1ddHe3g5AW1sbHR0dA10u/V+ATjvt9NhPd3d3j6r6tGq6f7uvr49qqs7z\nkDQR2AycCewANgCLIqKnKM8CYGlEnCNpLrAiIuZWKitpHvBl4IyIeL7oXCcDNwCnkXVX3Q3MKp3U\n4XkeZmb1VWmeR9WWR0Tsl7SM7O6oCcDK9OW/JDsc10XEGkkLJG0B9gCLK5VNp76KrGVyd7qZ6v6I\nuDAiHpN0E/AYsA+40FHCzGx08QxzG6RQKAw0Z81aga/Z+vAMczMb07y2VeO55WFmLc9rW9WHWx5m\nZjaiHDxskOLb9sxaQ6HZFRh3HDzMzKxmHvMws5bnMY/68JiHmbWMKVOyYFDLC2rLP2VKcz/jWODg\nYYN4zMOaadeurBVRy2vdukJN+XftavanbH0OHjZI/zpBZmZDcfCwQe68885mV8GsJp5d3ngeMLdB\nXvOa1/DCCy80uxo2TjVi8NsD7Pkc1MKINj4UCoWBsY4XX3yR5cuXA9kvOv+qs9HOa1s1noOHAdk4\nR+ma/pA9J8X/KM2slIOHAbBly5ZfegBM//aWLVuaUyGzGvgHTuN5zGOcksp2Y1blv7nVm8c8Ro+D\nniQoaZ6kxyU9IenSIfJcKalXUrekjmplJb1f0qOS9kuaXbT/GEkvSdqYXtfk/6iWV0QM+YIJFY6Z\njT6em9R4VbutJE0AriZ7lOwzwAOSbouIx4vyzAdmRsQsSacB1wJzq5R9BPht4Gtl3nZLRMwus98a\n4ohmV8DMRrk8LY85QG9EbI2IfcBqYGFJnoXAKoCIWA9MljS1UtmI2BwRvUC5JtHw+lRshHym2RUw\nq4nHPBovT/CYBmwrSj+d9uXJk6dsOe2py2qdpNNz5LcRdVGzK2Bmo1y9ZpgfTMvhGWBG6rb6BHCj\npMNGplqWT6HZFTCricc8Gi/PrbrbgRlF6elpX2meo8vkmZSj7C9J3Vu70vZGSU8CxwMbS/N2dXXR\n3t4OZPMROjo6Bpqv/ReT0047PfbT/eux5c0PBQqF0VP/0ZLu3y6+bX8oVW/VlTQR2Ew26L0D2AAs\nioieojwLgKURcY6kucCKiJibs+w64OKIeDCljwR2RsQBSccC9wKnRsTuknr5Vt06Wb48e5k1g2/V\nHT0OanmSiNgvaRmwlqyba2VE9Ehakh2O6yJijaQFkrYAe4DFlcqmSp0HXAUcCdwhqTsi5gNnAFdI\n2gscAJaUBg6rLwcOM6vGkwRtkILXCbImGk6roNZr1i2PfPwkQTMzG1FueZjZqOIxj9HDLQ8zMxtR\nDh42SFdXodlVMKtJ8a2m1hgOHjbI9dc3uwZmNtp5zMMGcX+wNZPHPEYPj3mYmdmIcvCwMgrNroBZ\nTTzm0XgOHmZmVjMHDxvk8ss7m10Fs5p4RYTG84C5mY0qHjAfPTxgbjVx/7G1Gl+zjefgYWZmNXO3\nlZmNKu62Gj3cbWVmZiMqV/CQNE/S45KekHTpEHmulNQrqVtSR7Wykt4v6VFJ+yXNLjnXZelcPZLO\nGu6Hs+Hx2lbWajzm0XhVg4ekCcDVwNnAKcAiSSeW5JkPzIyIWcAS4NocZR8BfpvsMbPF5zoJOB84\nCZgPXCOpbLPJ6sNrW5lZNXlaHnOA3ojYGhH7gNXAwpI8C4FVABGxHpgsaWqlshGxOSJ6gdLAsBBY\nHREvR0Qf0JvOYw3T2ewKmNXE8zwaL0/wmAZsK0o/nfblyZOnbLX3256jjJmZNVC9BszdzdTSCs2u\ngFlNPObReIfkyLMdmFGUnp72leY5ukyeSTnKlnu/cucapKuri/b2dgDa2tro6OgYaL72X0xOO+30\n2E93d3fXlB8KFAqjp/6jJd2/3dfXRzVV53lImghsBs4EdgAbgEUR0VOUZwGwNCLOkTQXWBERc3OW\nXQdcHBEPpvTJwA3AaWTdVXcDs0ondXieR/0sX569zJrB8zxGj0rzPKq2PCJiv6RlwFqybq6VEdEj\naUl2OK6LiDWSFkjaAuwBFlcqmyp1HnAVcCRwh6TuiJgfEY9Jugl4DNgHXOgo0VgOHGZWjWeY2yCF\nQqGoeW/WWMNpFdR6zbrlkY9nmJuZ2Yhyy8PMRhWPeYwebnmYmdmIcvCwQby2lbWa4ltNrTEcPGwQ\nr21lZtV4zMMGcX+wNZPHPEYPj3mYmdmIcvCwMgrNroBZTTzm0XgOHmZmVjOPeYxxU6bArl31f58j\njoCdO+v/PjYONOrZb/7+qKrSmIeDxxjXqIFBD0DaSPGA+ejhAXOrifuPrdX4mm08Bw8zM6uZu63G\nOHdbWatxt9Xo4W4rMzMbUbmCh6R5kh6X9ISkS4fIc6WkXkndkjqqlZV0hKS1kjZLukvS5LT/GEkv\nSdqYXtcc7Ie02rj/2FqNr9nGqxo8JE0ArgbOBk4BFkk6sSTPfGBmRMwClgDX5ij7SeCeiDgB+AFw\nWdEpt0TE7PS68GA+oJmZjbw8LY85QG9EbI2IfcBqYGFJnoXAKoCIWA9MljS1StmFQP8SfNcD5xWd\nr0E3els5foqgtRpfs42XJ3hMA7YVpZ9O+/LkqVR2akQ8CxARPwNeX5SvPXVZrZN0eo46mplZA9Vr\nwHw4LYf+ex92ADMiYjbwCeBGSYeNWM2sKvcfW6vxNdt4h+TIsx2YUZSenvaV5jm6TJ5JFcr+TNLU\niHhW0lHAcwARsRfYm7Y3SnoSOB7YWFqxrq4u2tvbAWhra6Ojo2Og+dp/MY33NDTm/aBAodD8z+v0\n+Ex3d3fXlN/Xa/l0/3ZfXx/VVJ3nIWkisBk4k6xVsAFYFBE9RXkWAEsj4hxJc4EVETG3UllJXwB2\nRsQX0l1YR0TEJyUdmfYfkHQscC9wakTsLqmX53nk4Hke1mo8z2P0qDTPo2rLIyL2S1oGrCXr5lqZ\nvvyXZIfjuohYI2mBpC3AHmBxpbLp1F8AbpL0UWArcH7afwZwhaS9wAFgSWngMDOz5vIM8zFuOL+w\nCoVCUfO+fu9jVk4jrllfr/l4hrmZmY0otzzGOI95WKtpxOM8/PyZfA5qzMPMrJGG8yPEP14az91W\nNkjxbXtmraHQ7AqMOw4eZmZWM495jHEe87DxwNdffXjMYxwL1JBlJqPof81s7HO31RgnIvtJVsOr\nsG5dzWXkwGFNdMEFhWZXYdxx8DCzltfV1ewajD8e8xjjPOZhZsPlGeZmZjaiHDxsEM/zsFbja7bx\nHDzMzKxmvlV3HKh9raDOmt/jiCNqLmI2YgqFTvwY88bygLkN4sFvazW+ZuvjoAfMJc2T9LikJ9JT\n/8rluVJSr6RuSR3Vyko6QtJaSZsl3SVpctGxy9K5eiSdlf+j2sgoNLsCZjUqNLsC407V4CFpAnA1\ncDZwCrBI0okleeYDMyNiFrAEuDZH2U8C90TECcAPgMtSmZPJnip4EjAfuEZqxCLN9oruZlfArEa+\nZhstT8tjDtAbEVsjYh+wGlhYkmchsAogItYDkyVNrVJ2IXB92r4eOC9tnwusjoiXI6IP6E3nsYbx\nU3+t1fiabbQ8wWMasK0o/XTalydPpbJTI+JZgIj4GfD6Ic61vcz7mdk4I2nIF/xZhWNWD/W6VXc4\n/495uKuB/A/RWk1EDPm64IILhjxm9ZHnVt3twIyi9PS0rzTP0WXyTKpQ9meSpkbEs5KOAp6rcq5B\n/GXWeP6b22h1/fXXV89kIyZP8HgAOE7SMcAO4IPAopI8twNLgW9LmgvsTkHh5xXK3g50AV8ALgBu\nK9p/g6TMktQ8AAAFRUlEQVSvkHVXHQdsKK3UULePmZlZ/VUNHhGxX9IyYC1ZN9fKiOiRtCQ7HNdF\nxBpJCyRtAfYAiyuVTaf+AnCTpI8CW8nusCIiHpN0E/AYsA+40BM6zMxGl5adJGhmZs3jta3GKEkf\nk/SYpOcl/c9hlP9hPeplNhySTpD0kKQHJR07nOtT0p9Jenc96jceueUxRknqAc6MiGeaXRezg5VW\np5gYEf+r2XWxjFseY5CkvwWOBb4v6SJJV6X9H5D0SPoFV0j7Tpa0XtLGtLTMzLT/haLzfSmV2yTp\n/LTvXZLWSfpOWkbmHxr+Qa1lSDomtYSvk/SopDslHZquodkpz2slPVWm7HzgIuCPJP1z2vdC+u9R\nku5N1+/Dkt4haYKkb6T0Jkl/nPJ+Q9J70/aZqcwmSV+X9Kq0/ylJy1MLZ5Ok4xvzF2o9Dh5jUET8\nEdntzZ3ALl6ZQ/OnwFkR8VaymfwAfwisiIjZwNvJJnLSX0bS+4A3R8SpwHuAL6XVAwA6gI8DJwMz\nJf2Xen4ua3nHAVdFxJvIpoS/j8HzuwZ1hUTE98mWPPpKRJxZku9DwJ3p+n0L2TolHcC0iHhzRLwF\n+Ebx+ST9Str3gXT8VcAfFWV5LiLelt7zkuF+2LHOwWNsK72d+YfA9ZJ+n1futPsx8GlJlwDtEfGf\nJWXeAXwLICKeI1uB7jfSsQ0RsSPdDdcNtI/4J7Cx5KmIeCRtb2RkrpcHgMWSPkv2I2cP8BPg1yV9\nVdLZwAslZU4AfhIRT6b09cAZRcdvSf99EDhmBOo4Jjl4jCMRcSHwabJJmA9KOiIivgX8FvD/gDWS\nOqucpjggFQea/fj5MFZZuevlZV75Hjq0/6Ckv0/dq3dUOmFE3Ef2xb8d+KakD0fEbrJWSIGsZf13\nZYpWmifWX09f0xU4eIxdg/5xSDo2Ih6IiMvJZvQfLenXI+KpiLiKbKLmm0vK3wf8TupHfh3wTspM\n2jTLodwXdh9ZdynAB/p3RsRHI+KtEfHfKp1L0gyybqaVwNeB2ZKmkA2u3wJ8BphdUnYzcIykY1P6\n9/Ca7jVzVB27yt1G9yVJs9L2PRHxsKRLJf0e2YTMHcBfFpePiFvSqgGbgAPAJRHxnKSTcryfWbFy\n4xt/BXxH0h8A3xvGuTqBSyTtI+ue+gjZkkbfUPZIiCB7/MNAmYj4T0mLge9KmkjW9fW1IepoQ/Ct\numZmVjN3W5mZWc0cPMzMrGYOHmZmVjMHDzMzq5mDh5mZ1czBw8zMaubgYWZmNXPwMGuyNFHNrKU4\neJgNg6Rfk3RHWn/p4bTc/dsl/SgtbX+/pFdL+pW0TtPDaZnvzlT+Akm3pSXG70n7Lpa0IZW/vJmf\nz6waL09iNjzzgO39ay9JOhx4iGyZ742SDiNbbPKPgQMR8WZJJwBri5aIeStwakT8u6T3ALMiYo4k\nAbdLOj0i/ERHG5Xc8jAbnkeA90j6nKTTgRnAMxGxESAiXoyI/cDpwD+mfZvJFgLsf8DQ3RHx72n7\nrHS+jWTLlZ8A9AcZs1HHLQ+zYYiI3vQEvAXAnwPrchYtXll2T8n+z0VEueXDzUYdtzzMhkHSG4D/\niIgbyVaGPQ14g6S3p+OHpYHw+4DfTfuOJ3uWyuYyp7wL+KikV6e8b0xL4JuNSm55mA3PqWRL3B8A\n9pI9xlTA1ZJ+FXgJ+K/ANcDfSnqYbNn7CyJiXzas8YqIuFvSicCP07EXgA8D/9agz2NWEy/JbmZm\nNXO3lZmZ1czBw8zMaubgYWZmNXPwMDOzmjl4mJlZzRw8zMysZg4eZmZWMwcPMzOr2f8HtPaEGbbk\nhnwAAAAASUVORK5CYII=\n",
"text/plain": [
"<matplotlib.figure.Figure at 0x7f2a35a0db50>"
]
},
"metadata": {},
"output_type": "display_data"
}
],
"source": [
"# Create a boxplot to view the distribution of\n",
"# fission and nu-fission rates in the pins\n",
"bp = df.boxplot(column='mean', by='score')"
]
},
{
"cell_type": "code",
"execution_count": 27,
"metadata": {
"collapsed": false
},
"outputs": [
{
"data": {
"text/plain": [
"<matplotlib.colorbar.Colorbar at 0x7f2a2af63d10>"
]
},
"execution_count": 27,
"metadata": {},
"output_type": "execute_result"
},
{
"data": {
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAVAAAAEZCAYAAADBv319AAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAIABJREFUeJzt3XuwHWWZ7/HvL4Ggcgn3IAlJIAlXHWNkQo632XiBgErQ\n41AwpQhOTcXBeM7UqIUctQiUZ0BmSmuQQWAmZcEoRjwiBEWMqPEIRyAYgyAJ2UETkwiRW5CA5LLz\nnD+6w6ysrFu/e/furLV/n6ouVvd6n+733ew8++3b+yoiMDOz4kZVXQEzs27lBGpmlsgJ1MwskROo\nmVkiJ1Azs0ROoGZmiZxAe5SkYyX9StLzkuZJ+qqkzw5if5dIumEo62jW7eTnQHuTpP8Ano+IT1Zd\nl6Em6XfA30bET6qui41s7oH2rknAb6quRFGSRlddB7NOOYH2IEk/Bk4F/k3SnyRNlfQ1SZfn3x8i\n6Q5Jz0l6RtLPamIvlrQ+j1sh6dR8+6WS/rOm3FmSHpH0rKSfSDq+5rvfSfqkpIfyY3xT0pgmdf2I\npHskfUnS08Clko6R9GNJT0v6o6SvSzogL38TMBG4I6/jp/LtsyTdmx/vV5L+ash/sGZ1nEB7UES8\nE/g58PGIOCAiVtcV+SSwDjgEOBz4X5BdNwU+DrwpIg4ATgfW1O66ptzNwP8ADgN+QJbQ9qop+9fA\nacDRwBuAC1pU+RRgdV6X/w0I+CfgCOAEYAIwP2/b+cDvgffmbfsXSUcC3wMuj4iDgE8B35F0SJsf\nldmgOIGOTNuA1wJHR8RARNybbx8AxgCvk7RXRPw+In7XIP4c4HsR8ZOIGAD+BXg18OaaMv8aERsj\nYhNwBzC9RX02RMS1EbEjIrZExOMR8eOI2B4RzwBfBup7lKr5/CHg+xHxQ4CI+DHwIHBmBz8Ls2RO\noCPTPwOPA4slrZZ0MUBEPA78A1lvb6OkmyUd0SD+SGDtzpXI7kSuA8bXlNlY8/klYL8W9VlXuyLp\n8Py0f72kTcDXgUNbxE8CzskvJzwr6TngLWR/JMxK4wQ6AkXE5oj4VERMAc4C/nHntc6IWBgRbyNL\nSgBfbLCLP9R8v9NRwPrUKtWt/xOwAzgpIg4k62GqRfl1wE0RcXC+HBQR+0fEVYn1MeuIE+gIJOk9\nkqbkqy8A24Ed+bOjp+Y3fLYCfyZLZPVuAd6Tl90rv5HzMvCLIari/sBm4AVJ44FP133/JHBMzfrX\ngfdJOk3SKEmvkvRX+bVRs9I4gfauVg/4TgPulvQCcC/wbxHxM2Af4ErgKbJe5mHAJbvtOGIVWa/w\nmrzse4D3RcT2Do7dicuANwE7r59+p+77K4HP56fr/xgR64E5ZDfDniK7vPAp/PttJfOD9GZmifwX\n2swskROomVkiJ1Azs0ROoGZmqSJij13I7uZ68eKlgmWw/37HFjvemiY5YDawElgFXNykzNVAP7Ac\nmN5pLNkrzTuAg2u2XZLvawVwWrs21r67vIeKJtvn50udjyUc4h+Kh5xw3LKEA8EUHi8c8zb+b8Pt\nP5p/P++ef0rD76YmHOd9z/+gcMze/YVDMv+aELNv483zH4T5Jzf+bv31xQ9zR/EQ/pQQs719kYa2\nNdi2BOhrEXNCwWOcW7B8I88DX+iw7Od2fzEDSaPIHpV7J9ljdUsl3R4RK2vKnAFMiYhpkk4BrgNm\ntYuVNAF4NzVv1Ek6gew15Z3jL9wtaVq0eFSpslN4SbMlrZS0auerhGbWW/bucGliJtAfEWsjYhuw\nkOx531pzgJsAIuJ+YKykcR3EfpndX9CYAyzMx2BYQ9YTndmqfZUk0Jq/DqcDJwHn1Q6HZma9Ya8O\nlybGs+s4CevZdbyFVmWaxko6C1gXEQ+32deGBsfbRVWn8K/8dQCQtPOvw8qWUbvoK6NeXeWYvpb/\nb0eMPr+wyeSqK9DEq4f/kGr5pfRqsjfW3j0UB6sqgTb669Cyq7y7vqGrTZea0jeh6irsEZxA99wE\n2uz0/LF8aWMD2eDZO03It9WXOapBmTFNYqeQ/bgekqR8+zJJMzs83i664CbS/JrPfThxmg293wCP\nlrDfZgnmpHzZqcmNu6XAVEmTgCfI7m2dV1dmEdkg4N+SNAvYFBEb89kNdouNiBVkA3UDr8yvNSMi\nnpO0CPiGpC+RdfKmAg+ktK9sBTL9/PJrYzbC1Se0+tFbUrW4QdRWRAxImgcsJrtfsyAiVkiam30d\nN0TEnZLOlLQaeBG4sFVso8OQn/ZHxKOSbiH7W7INuKjVHXioLoF28pfFzLrcYBNMRNwFHFe37fq6\n9XmdxjYoc0zd+hXAFZ3Wr5IEWuCvg5l1scH0QLtBZddAO/nrYGbdzQnUzCxRBY8xDasuSKAFX5Jb\nf0DxQzxZPOR3R0wuHgT0jV1SOGYr+xSOeT31zwi3t2lsq3nfGjvskc2FY4Bs0uOiEl7LnHBawnEW\nFw8p+qoktHk+poXXD+OxBqsLEsyg9Hr7zKxCPoU3M0vU6wmm19tnZhVyD9TMLFGvJ5heb5+ZVcg9\nUDOzRH6MycwskXugZmaJej3B9Hr7zKxCe3eaYVIniKqYE6iZlWYvJ1AzszR7j666BuVyAjWz0nTc\nA+1SXdC8goODNJkfvKWE8TD2H5s2iMaBbCocM4YthWM2cWDhmGc4pHDMYYc/VDgGgAcTYsamHaqo\ndyXEpAxmmzIACbSZJrKJqu6G7118HJyu0gUJ1My6Vo9nmB5vnplVqsczzKiqK2BmPWyvDpcmJM2W\ntFLSKkkXNylztaR+ScslTW8XK+lySQ9J+pWkuyQdkW+fJOklScvy5dp2zXMCNbPyjO5waUDSKOAa\n4HSySUPPk3R8XZkzgCkRMQ2YC1zXQexVEfGGiHgj8H3g0ppdro6IGflyUbvmOYGaWXkG1wOdCfRH\nxNqI2AYsBObUlZkD3AQQEfcDYyWNaxUbEbV3gPcFdtSsq2jzzMzKMbi78OOBdTXr68kSY7sy49vF\nSvoCcD6wCTi1ptxkScuA54HPR8Q9rSroHqiZlWeQ10ATdNSDjIjPRcRE4BvAJ/LNTwATI2IG8Eng\nZkktJwpzD9TMytMkwyx5Hpa0ny9yAzCxZn0Cu8+PtwE4qkGZMR3EAtwM3AnMj4itwFaAiFgm6XHg\nWGBZswo6gZpZeZrcIOo7OFt2uqzxtKFLgamSJpH1Ds8Fzqsrswj4OPAtSbOATRGxUdLTzWIlTY2I\n1Xn82eTvQUg6FHg2InZIOgaYCvy2VfOcQM2sPIPIMBExIGke2UTTo4AFEbFC0tzs67ghIu6UdKak\n1cCLwIWtYvNdXynpWLKbR2uBj+Xb3w5cLmlr/t3ciGj56qATqJmVZ5AZJiLuAo6r23Z93fq8TmPz\n7R9sUv5W4NYi9XMCNbPy9HiG6fHmmVmlPJhIl3k6IWZC8ZCprG5fqIGtjCkcs2r3s5C29mKgcMwn\nnv9q4Zhd7nMWcWdCzFHti+wmYZikaX+dEPNM8Rj+mBADfP2R4jEfekvBgHuLH6Oh3sswu+jx5plZ\npTygsplZoh7PMD3ePDOrVI9nmB5vnplVyqfwZmaJejzD9HjzzKxSr6q6AuVyAjWz8vgU3swsUY9n\nmB5vnplVqsczTI83z8wq5VN4M7NEPZ5herx5ZlapHs8wXdC8gqM0bD6k+CGWFA/5Be8oHgSMfkPx\nQT4OpOWYrg29jZ8Xjlk39ojCMcf87MnCMQCMTYjZNyGm+I8b+hNiUgYGmZYQA4xLiOkfqsFBivJo\nTGZmiXo8w/R488ysUj2eYXq8eWZWqR6/C+954c2sPIOcF17SbEkrJa2SdHGTMldL6pe0XNL0drGS\nLpf0kKRfSbpL0hE1312S72uFpNPaNc8J1MzKM4gEKmkUcA1wOnAScJ6k4+vKnAFMiYhpwFzgug5i\nr4qIN0TEG4HvA5fmMScC5wAnAGcA10pSq+Y5gZpZeUZ3uDQ2E+iPiLURsQ1YCMypKzMHuAkgIu4H\nxkoa1yo2IjbXxO9LNoUxwFnAwojYHhFryJ7HmNmqeb4GamblGdxoTOOBdTXr69k9oTUqM75drKQv\nAOcDm4BTa/b1i5qYDfm2ppxAzaw8TTLMkt9kSwlannLvFBGfAz6XXxv9BDA/5WBOoGZWnian531/\nkS07XfbthsU2sOu8rxPybfVljmpQZkwHsQA3k10Hnd9iX035GqiZlWdwd+GXAlMlTZI0BjgXWFRX\nZhHZqTiSZgGbImJjq1hJU2vizwZW1uzrXEljJB0NTAUeaNc8M7NyDCLDRMSApHnAYrLO3oKIWCFp\nbvZ13BARd0o6U9Jq4EXgwlax+a6vlHQs2c2jtcDH8phHJd0CPApsAy6KiCipeWZmbQzyQfqIuAs4\nrm7b9XXr8zqNzbd/sMXxrgCu6LR+XZBACw4Osj3hEOsTYoqP7wHAcTxWOGZ/XigccyR/KByzavff\ntbaOeWPiYCIJY74kDdixPCEmZaCTIxNiUn5XgYMTYqadUjDg/oSDNOI5kczMEvX4q5yVJVBJa4Dn\nya5DbIuIlg+smlkX6vEuWpXN2wH0RcRzFdbBzMrkBFoa4ceozHpbjyfQKhNYAD+StFTS31VYDzMr\ny+Dehd/jVfn34S0R8YSkw8gS6YqIuKfC+pjZUOvxHmhlzYuIJ/L/PiXpu2Qv+jdIoPNrPvfli5kN\npSV/ypYh5zmRhp6k1wCjImKzpH2B04DLGpeeP3wVMxuh+g7Ilp0ua/kGeAHugZZiHPBdSZHX4RsR\nsbiiuphZWZxAh15E/A6Y3ragmXU3J1AzszTRxXfYO+EEamalGejxDNMFzdtYrPjKccUPsV/xkL1f\nl3bL8v/x5sIx7+OOwjGPcmLhmMmsKRzzcsrIFsDDR72ucMwUVheOOfiHLxeOSRrkY0tCzMkJMcCb\nio8TU/yC2RANJuIEamaWaMs+YzosubXUepTFCdTMSjMwurcvgjqBmllpBrr5Pc0OOIGaWWm293gC\n9WhIZlaaAfbqaGlG0mxJKyWtyqcgblTmakn9kpZLmt4uVtJVklbk5b8j6YB8+yRJL0lali/Xtmuf\nE6iZlWaA0R0tjUgaBVwDnA6cBJwn6fi6MmcAUyJiGjAXuK6D2MXASRExHegHLqnZ5eqImJEvF7Vr\nnxOomZVmMAmUbICh/ohYGxHbgIXAnLoyc4CbACLifmCspHGtYiPi7ojYkcffRzb/+04q0j4nUDMr\nzRbGdLQ0MR5YV7O+Pt/WSZlOYgE+CvygZn1yfvr+U0lvbdc+30Qys9K0ur5Zko57kJI+SzYf2835\npj8AEyPiOUkzgNsknRgRm5vtwwnUzErT7PR86ZKXeHDJS+3CNwATa9Yn5NvqyxzVoMyYVrGSLgDO\nBN6xc1t+qv9c/nmZpMeBY4FlzSroBGpmpWmWQGf07c+Mvv1fWb/usmcbFVsKTJU0CXgCOBc4r67M\nIuDjwLckzQI2RcRGSU83i5U0G/g08PaIeOUlXEmHAs9GxA5JxwBTgd+2ap8TqJmVZjDPgUbEgKR5\nZHfNRwELImKFpLnZ13FDRNwp6UxJq4EXgQtbxea7/gpZD/VHkgDuy++4vx24XNJWslmD50bEplZ1\nVEQkN7Bs2YDLq4oFTZ9W/EDnFg/hwIQYYMbc4tM+HZ0wyMdoBgrHpAwmckriqBOnb/lh4ZjR23e0\nL1TnVd8uHJI92FJUyqAqb0yIgTZ9oiYOL1ZccyAiCt2R3m0fUvwiOhvF5L9p+aCPVwX3QM2sNH6V\n08ws0dbmjyj1BCdQMytNr78L7wRqZqWp4DnQYdXbrTOzSvkaqJlZIidQM7NEvgZqZpZoK/tUXYVS\nOYGaWWl8Cm9mlsin8GZmifwYk5lZIp/CV+7VxYo/mXCI9QkxiYOJbOKgwjFbeKJwzAvs375QnZTB\nRB7m9YVjADbtU/wHOGafrYVj3n/edwvH7NtffNCSpH9J30+IAWg7TnoDK9oXKYMTqJlZIidQM7NE\nW/wYk5lZGvdAzcwSOYGamSXq9edA284LL+kTkorfOjazEW+AvTpaulXbBAqMA5ZKukXSbOWzMJmZ\ntTPA6I6WZvKcs1LSKkkXNylztaR+ScslTW8XK+kqSSvy8t+RdEDNd5fk+1oh6bR27WubQCPic8A0\nYAFwAdAv6Z8kTWkXa2Yj22ASqKRRwDXA6cBJwHmSjq8rcwYwJSKmAXOB6zqIXQycFBHTyaYQvCSP\nORE4BzgBOAO4tl2HsZMeKJFN3flkvmwHDgL+j6SrOok3s5FpC2M6WpqYCfRHxNqI2AYsBObUlZkD\n3AQQEfcDYyWNaxUbEXdHxM63Je4DJuSfzwIWRsT2iFhDllxntmpfJ9dA/6ekXwJXAfcCr4+Ivwfe\nBPz3dvFmNnIN8hroeGBdzfr6fFsnZTqJBfgocGeTfW1oEvOKTq7eHgx8ICLW1m6MiB2S3ttBvJmN\nUM1Oz9cveZwNSx4v45Ad36OR9FlgW0R8M/VgbRNoRFza4ruK3rA1s27QLIG+tu9YXtt37CvrD1x2\nd6NiG4CJNesT8m31ZY5qUGZMq1hJFwBnAu/oYF9NdXQN1MwsxXZGd7Q0sRSYKmmSpDHAucCiujKL\ngPMBJM0CNkXExlaxkmYDnwbOiogtdfs6V9IYSUcDU4EHWrWvCx7A+lOx4pvKqcVu9ksL++1vTioc\nM/qk7YVjjkwYwelRTiwccwhPF46BtHek38y9hWNW75PwsMjrioek/LwPe2Zz8QMBHF485KlTCv7C\nXphYtzqDecYzIgYkzSO7az4KWBARKyTNzb6OGyLiTklnSloNvAhc2Co23/VXyHqoP8pvst8XERdF\nxKOSbgEeBbYBF+U30JvqggRqZt1qsK9yRsRdwHF1266vW5/XaWy+fVqL410BXNFp/ZxAzaw0W5s/\notQTnEDNrDS9/i68E6iZlaab33PvRG+3zswq5eHszMwSOYGamSXyNVAzs0S+BmpmlsiPMZmZJfIp\nvJlZIp/Cm5kl8l34yhUcSKP4uBvwYELM1IQYSKrfmmeOLhxz6CHPFI5JGRhkH7YWjgE4kOcKxzzD\noYVjxvHHwjEpg6r8mdcUjnn4LcVjAPbnhcIxoxkoGDE0I1U6gZqZJXICHQRJC4D3Ahsj4i/ybQcB\n3wImAWuAcyLi+TLrYWbVSBm2sJuUPaDy18hmxav1GeDuiDgO+An5jHhm1nsGO63xnq7UBBoR98Bu\nF7vmADfmn28Ezi6zDmZWnV5PoFVcAz08H3KfiHhSUsL42mbWDfwcaPlaDpkPX635fDLwl2XWxWxE\nenDJi/xyyUtDvl8/Bzr0NkoaFxEbJR0B7Z4z+fthqZTZSHZy376c3LfvK+s3XJY211W9bj4978Rw\nzMopdp2reRFwQf75I8Dtw1AHM6vAYK+BSpotaaWkVZIublLmakn9kpZLmt4uVtIHJT0iaUDSjJrt\nkyS9JGlZvlzbrn1lP8Z0M9AHHCLp98ClwJXAtyV9FFgLnFNmHcysOlu2pg8mImkUcA3wTuAPwFJJ\nt0fEypoyZwBTImKapFOA64BZbWIfBt4PXM/uVkfEjAbbGyo1gUbE3zT56l1lHtfM9gwD2weVYmYC\n/RGxFkDSQrKneFbWlJkD3AQQEfdLGitpHHB0s9iIeCzfVntmvFOjbU0Nxym8mY1QA9tHd7Q0MR5Y\nV7O+Pt/WSZlOYhuZnJ++/1TSW9sV7u1bZGZWqWbJccc9P2fHvfeUcchCPcg6fwAmRsRz+bXR2ySd\nGBGbmwV0QQLdWKz49m3FD3HE3sVjlhQPSfah4iG/eOgdhWMenHBy4Ziph6wuHAPwUsLgG79MGJx3\nI8PzmPEaJg/LcQCO5A+FY4o/jzk0g4ls39bkuKf0ZctOV13ZqNQGYGLN+oR8W32ZoxqUGdNB7C4i\nYhv5iz8RsUzS48CxwLJmMT6FN7PS7BjYq6OliaXA1Pzu+BjgXLKneGotAs4HkDQL2JS/qNNJLNT0\nWCUdmt98QtIxZGOu/bZV+7qgB2pmXav59c22ImJA0jxgMVlnb0FErJA0N/s6boiIOyWdKWk18CJw\nYatYAElnA18BDgW+J2l5RJwBvB24XNJWYAcwNyI2taqjE6iZleflwaWYiLgLOK5u2/V16/M6jc23\n3wbc1mD7rcCtRernBGpm5UkZ4LyLOIGaWXmcQM3MEjmBmpklSniqsJs4gZpZeYrOZddlnEDNrDw+\nhTczS/Ry1RUolxOomZXHPVAzs0ROoFU7oGD5hIFB7ioewtSEGMimdSpo2zVFfwYkjbi67ZHix1l3\n8lHtCzWw8VXjCsfsP/qFwjEpU0qkDHSSEnMijxaOAfgppxaOKT4AyY3ti3TCCdTMLJEfYzIzS+TH\nmMzMEvkU3swskR9jMjNL5B6omVkiJ1Azs0ROoGZmiXr8MSZPKmdm5RnocGlC0mxJKyWtknRxkzJX\nS+qXtFzS9Haxkj4o6RFJA/n0xbX7uiTf1wpJp7VrnnugZlaeQdyFz2fIvAZ4J9mc7Usl3R4RK2vK\nnAFMiYhpkk4BrgNmtYl9GHg/cH3d8U4AzgFOIJsG+W5J0yIimtXRPVAzK8/2DpfGZgL9EbE2n7N9\nITCnrswc4CaAiLgfGCtpXKvYiHgsIvqpmdK4Zl8LI2J7RKwB+vP9NOUEambl2dbh0th4YF3N+vp8\nWydlOoltd7wN7WK64BT+TwXL/7r4IV4+oXjM5oRBSwCeTohpOTN1E0sSYhJsfvqwtMDJxUOeTbij\nu/bl4wvHjDr+xcIxR44rOlgHLN8yvX2hBsbss7VwzMMDr0861qAN/6uc9b3KUnVBAjWzrtXsj976\nJbBhSbvoDcDEmvUJ+bb6Mkc1KDOmg9hGx2u0r6acQM2sPM0S6BF92bLTA5c1KrUUmCppEvAEcC5w\nXl2ZRcDHgW9JmgVsioiNkp7uIBZ27bEuAr4h6ctkp+5TgQdatM4J1MxKNIjnQCNiQNI8YDHZ/ZoF\nEbFC0tzs67ghIu6UdKak1cCLwIWtYgEknQ18BTgU+J6k5RFxRkQ8KukW4NG85he1ugMPoDbfV0pS\nZO0vovhAvdlTCwVNTrwGenZCTMo10NQBn4sqfokxMzkhJuWtloTHaIbrGugLW/YvHANp10AHBooN\nLP3sXhOIiEFdT5QUfLjD/PKfGvTxquAeqJmVx69ympkl6vFXOZ1Azaw8HpHezCyRT+HNzBI5gZqZ\nJfI1UDOzRFuqrkC5nEDNrDw+ha/acPwf2Fg8ZP2EtEPdlRCTcqiV7YvsZnJCTMrgKADLE2JelRCT\n8LPbsd++hWPWL5xW/ECpv9qzEmKq+pfuU3gzs0R+jMnMLJFP4c3MEjmBmpkl8jVQM7NEfozJzCyR\nT+HNzBL5FN7MLJEfYzIzS+RTeDOzRD2eQEdVXQEz62HbOlyakDRb0kpJqyRd3KTM1ZL6JS2XNL1d\nrKSDJC2W9JikH0oam2+fJOklScvy5dp2zXMCNbPybO9waUDSKOAa4HTgJOA8ScfXlTkDmBIR04C5\nwHUdxH4GuDsijgN+AlxSs8vVETEjXy5q1zwnUDPbU80E+iNibURsAxYCc+rKzAFuAoiI+4Gxksa1\niZ0D3Jh/vpFd58otNDNoF1wDLTpSUspUsRuKh2x/dcJxgDWHFI9JGVkpZQSn+xJiDkyIgbTpkBOm\nKE4awel1CTH7JcSsToiBLBUUlfr/qVrjgXU16+vJEmO7MuPbxI6LiI0AEfGkpMNryk2WtAx4Hvh8\nRNzTqoJdkEDNrPcsyZchlzK3/M7J658AJkbEc5JmALdJOjEiNjcLLPUUXtICSRsl/bpm26WS1tdc\nqJ1dZh3MrErN7hq9BfhszdLQBmBizfoEdj9d3AAc1aBMq9gn89N8JB0B/BEgIrZGxHP552XA48Cx\nrVpX9jXQr5FdxK33pZoLtSlDDJtZVxjEXSRYCkzN746PAc4FFtWVWQScDyBpFrApPz1vFbsIuCD/\n/BHg9jz+0PzmE5KOAaYCv23VulJP4SPiHkmTGnyV0s02s66T/i5nRAxImgcsJuvsLYiIFZLmZl/H\nDRFxp6QzJa0GXgQubBWb7/qLwC2SPgqsBc7Jt78duFzSVmAHMDciNrWqY1XXQOdJ+jDwIPDJiHi+\nonqYWan+PKjo/Az1uLpt19etz+s0Nt/+LPCuBttvBW4tUr8qEui1wOUREZK+AHwJ+NvmxW+r+Xw8\nabdvzaylZ5fAc0tK2HFvjyYy7Ak0Ip6qWf134I7WEWe3/trMBu/gvmzZ6XeXDdGOe/tdzuFIoKLm\nmqekIyLiyXz1A8Ajw1AHM6uEe6DJJN0M9AGHSPo9cClwav6+6g5gDdnrV2bWk9wDTRYRf9Ng89fK\nPKaZ7UncAzUzSzS4u/B7OidQMyuRT+ErtqZg+ZRThnEJMQ8kxAAvp9TvhOIh64sOwgIwuXjIpsQe\nxn0HF4/ZL2Egls0JP+/7Utp0QELMMwkxqceq6p+6T+HNzBK5B2pmlsg9UDOzRO6Bmpklcg/UzCyR\nH2MyM0vkHqiZWSJfAzUzS9TbPdAuntZ4TdUV2APcX3UF9gzbl1Rdgz3Az6quQBODmtJjj+cE2tWc\nQAEYWFJ1DfYAe2oCbTapXP3SnXwKb2Yl6t7eZSecQM2sRL39GJMion2pikjacytn1uMiYlCz50pa\nAzSalbeRtRExeTDHq8IenUDNzPZkXXwTycysWk6gZmaJui6BSpotaaWkVZIurro+VZG0RtJDkn4l\nKXF05+4jaYGkjZJ+XbPtIEmLJT0m6YeSxlZZx7I1+RlcKmm9pGX5MrvKOo4UXZVAJY0CrgFOB04C\nzpN0fLW1qswOoC8i3hgRM6uuzDD6Gtn//1qfAe6OiOOAnwCXDHuthlejnwHAlyJiRr7cNdyVGom6\nKoECM4H+iFgbEduAhcCciutUFdF9//8GLSLuAZ6r2zwHuDH/fCNw9rBWapg1+RlA9jthw6jb/gGO\nB9bVrK/Pt41EAfxI0lJJf1d1ZSp2eERsBIiIJ4HDK65PVeZJWi7pP3r9MsaeotsSqP2Xt0TEDOBM\n4OOS3lp1hfYgI/HZvGuBYyJiOvAk8KWK6zMidFsC3QBMrFmfkG8bcSLiify/TwHfJbu8MVJtlDQO\nQNIRwB9UDeVFAAABiElEQVQrrs+wi4in4r8e6v534C+rrM9I0W0JdCkwVdIkSWOAc4FFFddp2El6\njaT98s/7AqcBj1Rbq2Eldr3etwi4IP/8EeD24a5QBXb5GeR/OHb6ACPr96EyXfUufEQMSJoHLCZL\n/gsiYkXF1arCOOC7+auuewHfiIjFFddpWEi6GegDDpH0e+BS4Erg25I+CqwFzqmuhuVr8jM4VdJ0\nsqcz1gBzK6vgCOJXOc3MEnXbKbyZ2R7DCdTMLJETqJlZIidQM7NETqBmZomcQM3MEjmBmpklcgI1\nM0vkBGpDStLJ+UDPYyTtK+kRSSdWXS+zMvhNJBtyki4HXp0v6yLiixVXyawUTqA25CTtTTbwy5+B\nN4d/yaxH+RTeynAosB+wP/CqiutiVhr3QG3ISbod+CZwNHBkRHyi4iqZlaKrhrOzPZ+kDwNbI2Jh\nPgngvZL6ImJJxVUzG3LugZqZJfI1UDOzRE6gZmaJnEDNzBI5gZqZJXICNTNL5ARqZpbICdTMLJET\nqJlZov8PnPpCLiDrzj8AAAAASUVORK5CYII=\n",
"text/plain": [
"<matplotlib.figure.Figure at 0x7f2a3592a610>"
]
},
"metadata": {},
"output_type": "display_data"
}
],
"source": [
"# Extract thermal nu-fission rates from pandas\n",
"fiss = df[df['score'] == 'nu-fission']\n",
"fiss = fiss[fiss['energy low [MeV]'] == 0.0]\n",
"\n",
"# Extract mean and reshape as 2D NumPy arrays\n",
"mean = fiss['mean'].reshape((17,17))\n",
"\n",
"pylab.imshow(mean, interpolation='nearest')\n",
"pylab.title('fission rate')\n",
"pylab.xlabel('x')\n",
"pylab.ylabel('y')\n",
"pylab.colorbar()"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"**Analyze the cell+nuclides scatter-y2 rate tally**"
]
},
{
"cell_type": "code",
"execution_count": 28,
"metadata": {
"collapsed": false
},
"outputs": [
{
"name": "stdout",
"output_type": "stream",
"text": [
"Tally\n",
"\tID =\t10001\n",
"\tName =\tcell tally\n",
"\tFilters =\t\n",
" \t\tcell\t[10000]\n",
"\tNuclides =\tU-235 U-238 \n",
"\tScores =\t[u'scatter-Y0,0', u'scatter-Y1,-1', u'scatter-Y1,0', u'scatter-Y1,1', u'scatter-Y2,-2', u'scatter-Y2,-1', u'scatter-Y2,0', u'scatter-Y2,1', u'scatter-Y2,2']\n",
"\tEstimator =\tanalog\n",
"\n"
]
}
],
"source": [
"# Find the cell Tally with the StatePoint API\n",
"tally = sp.get_tally(name='cell tally')\n",
"\n",
"# Print a little info about the cell tally to the screen\n",
"print(tally)"
]
},
{
"cell_type": "code",
"execution_count": 29,
"metadata": {
"collapsed": false
},
"outputs": [
{
"data": {
"text/html": [
"<div>\n",
"<table border=\"1\" class=\"dataframe\">\n",
" <thead>\n",
" <tr style=\"text-align: right;\">\n",
" <th></th>\n",
" <th>cell</th>\n",
" <th>nuclide</th>\n",
" <th>score</th>\n",
" <th>mean</th>\n",
" <th>std. dev.</th>\n",
" </tr>\n",
" </thead>\n",
" <tbody>\n",
" <tr>\n",
" <th>0</th>\n",
" <td>10000</td>\n",
" <td>U-235</td>\n",
" <td>scatter-Y0,0</td>\n",
" <td>3.77e-02</td>\n",
" <td>6.49e-04</td>\n",
" </tr>\n",
" <tr>\n",
" <th>1</th>\n",
" <td>10000</td>\n",
" <td>U-235</td>\n",
" <td>scatter-Y1,-1</td>\n",
" <td>2.54e-04</td>\n",
" <td>1.81e-04</td>\n",
" </tr>\n",
" <tr>\n",
" <th>2</th>\n",
" <td>10000</td>\n",
" <td>U-235</td>\n",
" <td>scatter-Y1,0</td>\n",
" <td>3.65e-05</td>\n",
" <td>2.70e-04</td>\n",
" </tr>\n",
" <tr>\n",
" <th>3</th>\n",
" <td>10000</td>\n",
" <td>U-235</td>\n",
" <td>scatter-Y1,1</td>\n",
" <td>-1.70e-04</td>\n",
" <td>2.19e-04</td>\n",
" </tr>\n",
" <tr>\n",
" <th>4</th>\n",
" <td>10000</td>\n",
" <td>U-235</td>\n",
" <td>scatter-Y2,-2</td>\n",
" <td>7.47e-05</td>\n",
" <td>1.54e-04</td>\n",
" </tr>\n",
" <tr>\n",
" <th>5</th>\n",
" <td>10000</td>\n",
" <td>U-235</td>\n",
" <td>scatter-Y2,-1</td>\n",
" <td>-2.35e-04</td>\n",
" <td>1.34e-04</td>\n",
" </tr>\n",
" <tr>\n",
" <th>6</th>\n",
" <td>10000</td>\n",
" <td>U-235</td>\n",
" <td>scatter-Y2,0</td>\n",
" <td>-5.51e-05</td>\n",
" <td>1.79e-04</td>\n",
" </tr>\n",
" <tr>\n",
" <th>7</th>\n",
" <td>10000</td>\n",
" <td>U-235</td>\n",
" <td>scatter-Y2,1</td>\n",
" <td>-1.27e-04</td>\n",
" <td>1.54e-04</td>\n",
" </tr>\n",
" <tr>\n",
" <th>8</th>\n",
" <td>10000</td>\n",
" <td>U-235</td>\n",
" <td>scatter-Y2,2</td>\n",
" <td>1.72e-04</td>\n",
" <td>1.40e-04</td>\n",
" </tr>\n",
" <tr>\n",
" <th>9</th>\n",
" <td>10000</td>\n",
" <td>U-238</td>\n",
" <td>scatter-Y0,0</td>\n",
" <td>2.34e+00</td>\n",
" <td>7.62e-03</td>\n",
" </tr>\n",
" <tr>\n",
" <th>10</th>\n",
" <td>10000</td>\n",
" <td>U-238</td>\n",
" <td>scatter-Y1,-1</td>\n",
" <td>2.46e-02</td>\n",
" <td>1.71e-03</td>\n",
" </tr>\n",
" <tr>\n",
" <th>11</th>\n",
" <td>10000</td>\n",
" <td>U-238</td>\n",
" <td>scatter-Y1,0</td>\n",
" <td>1.15e-03</td>\n",
" <td>2.17e-03</td>\n",
" </tr>\n",
" <tr>\n",
" <th>12</th>\n",
" <td>10000</td>\n",
" <td>U-238</td>\n",
" <td>scatter-Y1,1</td>\n",
" <td>-2.39e-02</td>\n",
" <td>2.15e-03</td>\n",
" </tr>\n",
" <tr>\n",
" <th>13</th>\n",
" <td>10000</td>\n",
" <td>U-238</td>\n",
" <td>scatter-Y2,-2</td>\n",
" <td>-3.92e-03</td>\n",
" <td>1.38e-03</td>\n",
" </tr>\n",
" <tr>\n",
" <th>14</th>\n",
" <td>10000</td>\n",
" <td>U-238</td>\n",
" <td>scatter-Y2,-1</td>\n",
" <td>-1.19e-03</td>\n",
" <td>1.58e-03</td>\n",
" </tr>\n",
" <tr>\n",
" <th>15</th>\n",
" <td>10000</td>\n",
" <td>U-238</td>\n",
" <td>scatter-Y2,0</td>\n",
" <td>3.22e-03</td>\n",
" <td>1.45e-03</td>\n",
" </tr>\n",
" <tr>\n",
" <th>16</th>\n",
" <td>10000</td>\n",
" <td>U-238</td>\n",
" <td>scatter-Y2,1</td>\n",
" <td>1.27e-04</td>\n",
" <td>9.70e-04</td>\n",
" </tr>\n",
" <tr>\n",
" <th>17</th>\n",
" <td>10000</td>\n",
" <td>U-238</td>\n",
" <td>scatter-Y2,2</td>\n",
" <td>-2.70e-03</td>\n",
" <td>1.21e-03</td>\n",
" </tr>\n",
" </tbody>\n",
"</table>\n",
"</div>"
],
"text/plain": [
" cell nuclide score mean std. dev.\n",
"0 10000 U-235 scatter-Y0,0 3.77e-02 6.49e-04\n",
"1 10000 U-235 scatter-Y1,-1 2.54e-04 1.81e-04\n",
"2 10000 U-235 scatter-Y1,0 3.65e-05 2.70e-04\n",
"3 10000 U-235 scatter-Y1,1 -1.70e-04 2.19e-04\n",
"4 10000 U-235 scatter-Y2,-2 7.47e-05 1.54e-04\n",
"5 10000 U-235 scatter-Y2,-1 -2.35e-04 1.34e-04\n",
"6 10000 U-235 scatter-Y2,0 -5.51e-05 1.79e-04\n",
"7 10000 U-235 scatter-Y2,1 -1.27e-04 1.54e-04\n",
"8 10000 U-235 scatter-Y2,2 1.72e-04 1.40e-04\n",
"9 10000 U-238 scatter-Y0,0 2.34e+00 7.62e-03\n",
"10 10000 U-238 scatter-Y1,-1 2.46e-02 1.71e-03\n",
"11 10000 U-238 scatter-Y1,0 1.15e-03 2.17e-03\n",
"12 10000 U-238 scatter-Y1,1 -2.39e-02 2.15e-03\n",
"13 10000 U-238 scatter-Y2,-2 -3.92e-03 1.38e-03\n",
"14 10000 U-238 scatter-Y2,-1 -1.19e-03 1.58e-03\n",
"15 10000 U-238 scatter-Y2,0 3.22e-03 1.45e-03\n",
"16 10000 U-238 scatter-Y2,1 1.27e-04 9.70e-04\n",
"17 10000 U-238 scatter-Y2,2 -2.70e-03 1.21e-03"
]
},
"execution_count": 29,
"metadata": {},
"output_type": "execute_result"
}
],
"source": [
"# Get a pandas dataframe for the cell tally data\n",
"df = tally.get_pandas_dataframe()\n",
"\n",
"# Print the first twenty rows in the dataframe\n",
"df.head(100)"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"Use the new Tally data retrieval API with pure NumPy"
]
},
{
"cell_type": "code",
"execution_count": 30,
"metadata": {
"collapsed": false
},
"outputs": [
{
"name": "stdout",
"output_type": "stream",
"text": [
"[[[ 0.00121338 0.00761835]\n",
" [ 0.00013952 0.00064888]]]\n"
]
}
],
"source": [
"# Get the standard deviations for two of the spherical harmonic\n",
"# scattering reaction rates \n",
"data = tally.get_values(scores=['scatter-Y2,2', 'scatter-Y0,0'], \n",
" nuclides=['U-238', 'U-235'], value='std_dev')\n",
"print(data)"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"**Analyze the distribcell tally**"
]
},
{
"cell_type": "code",
"execution_count": 31,
"metadata": {
"collapsed": false
},
"outputs": [
{
"name": "stdout",
"output_type": "stream",
"text": [
"Tally\n",
"\tID =\t10002\n",
"\tName =\tdistribcell tally\n",
"\tFilters =\t\n",
" \t\tdistribcell\t[10002]\n",
"\tNuclides =\ttotal \n",
"\tScores =\t[u'absorption', u'scatter']\n",
"\tEstimator =\ttracklength\n",
"\n"
]
}
],
"source": [
"# Find the distribcell Tally with the StatePoint API\n",
"tally = sp.get_tally(name='distribcell tally')\n",
"\n",
"# Print a little info about the distribcell tally to the screen\n",
"print(tally)"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"Use the new Tally data retrieval API with pure NumPy"
]
},
{
"cell_type": "code",
"execution_count": 32,
"metadata": {
"collapsed": false
},
"outputs": [
{
"name": "stdout",
"output_type": "stream",
"text": [
"[[[ 0.03284934]]]\n"
]
}
],
"source": [
"# Get the relative error for the scattering reaction rates in\n",
"# the first 30 distribcell instances \n",
"data = tally.get_values(scores=['scatter'], filters=['distribcell'],\n",
" filter_bins=[(i,) for i in range(10)], value='rel_err')\n",
"print(data)"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"Print the distribcell tally dataframe **without** OpenCG info"
]
},
{
"cell_type": "code",
"execution_count": 33,
"metadata": {
"collapsed": false
},
"outputs": [
{
"data": {
"text/html": [
"<div>\n",
"<table border=\"1\" class=\"dataframe\">\n",
" <thead>\n",
" <tr style=\"text-align: right;\">\n",
" <th></th>\n",
" <th>distribcell</th>\n",
" <th>score</th>\n",
" <th>mean</th>\n",
" <th>std. dev.</th>\n",
" </tr>\n",
" </thead>\n",
" <tbody>\n",
" <tr>\n",
" <th>558</th>\n",
" <td>279</td>\n",
" <td>absorption</td>\n",
" <td>7.14e-05</td>\n",
" <td>8.26e-06</td>\n",
" </tr>\n",
" <tr>\n",
" <th>559</th>\n",
" <td>279</td>\n",
" <td>scatter</td>\n",
" <td>1.25e-02</td>\n",
" <td>5.75e-04</td>\n",
" </tr>\n",
" <tr>\n",
" <th>560</th>\n",
" <td>280</td>\n",
" <td>absorption</td>\n",
" <td>8.50e-05</td>\n",
" <td>6.16e-06</td>\n",
" </tr>\n",
" <tr>\n",
" <th>561</th>\n",
" <td>280</td>\n",
" <td>scatter</td>\n",
" <td>1.38e-02</td>\n",
" <td>4.58e-04</td>\n",
" </tr>\n",
" <tr>\n",
" <th>562</th>\n",
" <td>281</td>\n",
" <td>absorption</td>\n",
" <td>1.04e-04</td>\n",
" <td>7.54e-06</td>\n",
" </tr>\n",
" <tr>\n",
" <th>563</th>\n",
" <td>281</td>\n",
" <td>scatter</td>\n",
" <td>1.55e-02</td>\n",
" <td>4.15e-04</td>\n",
" </tr>\n",
" <tr>\n",
" <th>564</th>\n",
" <td>282</td>\n",
" <td>absorption</td>\n",
" <td>1.22e-04</td>\n",
" <td>9.98e-06</td>\n",
" </tr>\n",
" <tr>\n",
" <th>565</th>\n",
" <td>282</td>\n",
" <td>scatter</td>\n",
" <td>1.68e-02</td>\n",
" <td>5.73e-04</td>\n",
" </tr>\n",
" <tr>\n",
" <th>566</th>\n",
" <td>283</td>\n",
" <td>absorption</td>\n",
" <td>1.14e-04</td>\n",
" <td>8.02e-06</td>\n",
" </tr>\n",
" <tr>\n",
" <th>567</th>\n",
" <td>283</td>\n",
" <td>scatter</td>\n",
" <td>1.66e-02</td>\n",
" <td>5.44e-04</td>\n",
" </tr>\n",
" <tr>\n",
" <th>568</th>\n",
" <td>284</td>\n",
" <td>absorption</td>\n",
" <td>1.06e-04</td>\n",
" <td>8.37e-06</td>\n",
" </tr>\n",
" <tr>\n",
" <th>569</th>\n",
" <td>284</td>\n",
" <td>scatter</td>\n",
" <td>1.64e-02</td>\n",
" <td>5.14e-04</td>\n",
" </tr>\n",
" <tr>\n",
" <th>570</th>\n",
" <td>285</td>\n",
" <td>absorption</td>\n",
" <td>1.23e-04</td>\n",
" <td>9.19e-06</td>\n",
" </tr>\n",
" <tr>\n",
" <th>571</th>\n",
" <td>285</td>\n",
" <td>scatter</td>\n",
" <td>1.70e-02</td>\n",
" <td>5.34e-04</td>\n",
" </tr>\n",
" <tr>\n",
" <th>572</th>\n",
" <td>286</td>\n",
" <td>absorption</td>\n",
" <td>1.14e-04</td>\n",
" <td>6.70e-06</td>\n",
" </tr>\n",
" <tr>\n",
" <th>573</th>\n",
" <td>286</td>\n",
" <td>scatter</td>\n",
" <td>1.75e-02</td>\n",
" <td>5.68e-04</td>\n",
" </tr>\n",
" <tr>\n",
" <th>574</th>\n",
" <td>287</td>\n",
" <td>absorption</td>\n",
" <td>1.14e-04</td>\n",
" <td>8.10e-06</td>\n",
" </tr>\n",
" <tr>\n",
" <th>575</th>\n",
" <td>287</td>\n",
" <td>scatter</td>\n",
" <td>1.72e-02</td>\n",
" <td>4.93e-04</td>\n",
" </tr>\n",
" <tr>\n",
" <th>576</th>\n",
" <td>288</td>\n",
" <td>absorption</td>\n",
" <td>1.06e-04</td>\n",
" <td>1.07e-05</td>\n",
" </tr>\n",
" <tr>\n",
" <th>577</th>\n",
" <td>288</td>\n",
" <td>scatter</td>\n",
" <td>1.72e-02</td>\n",
" <td>7.73e-04</td>\n",
" </tr>\n",
" </tbody>\n",
"</table>\n",
"</div>"
],
"text/plain": [
" distribcell score mean std. dev.\n",
"558 279 absorption 7.14e-05 8.26e-06\n",
"559 279 scatter 1.25e-02 5.75e-04\n",
"560 280 absorption 8.50e-05 6.16e-06\n",
"561 280 scatter 1.38e-02 4.58e-04\n",
"562 281 absorption 1.04e-04 7.54e-06\n",
"563 281 scatter 1.55e-02 4.15e-04\n",
"564 282 absorption 1.22e-04 9.98e-06\n",
"565 282 scatter 1.68e-02 5.73e-04\n",
"566 283 absorption 1.14e-04 8.02e-06\n",
"567 283 scatter 1.66e-02 5.44e-04\n",
"568 284 absorption 1.06e-04 8.37e-06\n",
"569 284 scatter 1.64e-02 5.14e-04\n",
"570 285 absorption 1.23e-04 9.19e-06\n",
"571 285 scatter 1.70e-02 5.34e-04\n",
"572 286 absorption 1.14e-04 6.70e-06\n",
"573 286 scatter 1.75e-02 5.68e-04\n",
"574 287 absorption 1.14e-04 8.10e-06\n",
"575 287 scatter 1.72e-02 4.93e-04\n",
"576 288 absorption 1.06e-04 1.07e-05\n",
"577 288 scatter 1.72e-02 7.73e-04"
]
},
"execution_count": 33,
"metadata": {},
"output_type": "execute_result"
}
],
"source": [
"# Get a pandas dataframe for the distribcell tally data\n",
"df = tally.get_pandas_dataframe(nuclides=False)\n",
"\n",
"# Print the last twenty rows in the dataframe\n",
"df.tail(20)"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"Print the distribcell tally dataframe **with** OpenCG info"
]
},
{
"cell_type": "code",
"execution_count": 34,
"metadata": {
"collapsed": false
},
"outputs": [
{
"data": {
"text/html": [
"<div>\n",
"<table border=\"1\" class=\"dataframe\">\n",
" <thead>\n",
" <tr>\n",
" <th></th>\n",
" <th colspan=\"2\" halign=\"left\">level 1</th>\n",
" <th colspan=\"4\" halign=\"left\">level 2</th>\n",
" <th colspan=\"2\" halign=\"left\">level 3</th>\n",
" <th>distribcell</th>\n",
" <th>score</th>\n",
" <th>mean</th>\n",
" <th>std. dev.</th>\n",
" </tr>\n",
" <tr>\n",
" <th></th>\n",
" <th>cell</th>\n",
" <th>univ</th>\n",
" <th colspan=\"4\" halign=\"left\">lat</th>\n",
" <th>cell</th>\n",
" <th>univ</th>\n",
" <th></th>\n",
" <th></th>\n",
" <th></th>\n",
" <th></th>\n",
" </tr>\n",
" <tr>\n",
" <th></th>\n",
" <th>id</th>\n",
" <th>id</th>\n",
" <th>id</th>\n",
" <th>x</th>\n",
" <th>y</th>\n",
" <th>z</th>\n",
" <th>id</th>\n",
" <th>id</th>\n",
" <th></th>\n",
" <th></th>\n",
" <th></th>\n",
" <th></th>\n",
" </tr>\n",
" </thead>\n",
" <tbody>\n",
" <tr>\n",
" <th>0</th>\n",
" <td>10003</td>\n",
" <td>0</td>\n",
" <td>10001</td>\n",
" <td>0</td>\n",
" <td>16</td>\n",
" <td>0</td>\n",
" <td>10002</td>\n",
" <td>10000</td>\n",
" <td>0</td>\n",
" <td>absorption</td>\n",
" <td>1.30e-04</td>\n",
" <td>8.67e-06</td>\n",
" </tr>\n",
" <tr>\n",
" <th>1</th>\n",
" <td>10003</td>\n",
" <td>0</td>\n",
" <td>10001</td>\n",
" <td>0</td>\n",
" <td>16</td>\n",
" <td>0</td>\n",
" <td>10002</td>\n",
" <td>10000</td>\n",
" <td>0</td>\n",
" <td>scatter</td>\n",
" <td>1.98e-02</td>\n",
" <td>6.50e-04</td>\n",
" </tr>\n",
" <tr>\n",
" <th>2</th>\n",
" <td>10003</td>\n",
" <td>0</td>\n",
" <td>10001</td>\n",
" <td>0</td>\n",
" <td>15</td>\n",
" <td>0</td>\n",
" <td>10002</td>\n",
" <td>10000</td>\n",
" <td>1</td>\n",
" <td>absorption</td>\n",
" <td>2.24e-04</td>\n",
" <td>1.44e-05</td>\n",
" </tr>\n",
" <tr>\n",
" <th>3</th>\n",
" <td>10003</td>\n",
" <td>0</td>\n",
" <td>10001</td>\n",
" <td>0</td>\n",
" <td>15</td>\n",
" <td>0</td>\n",
" <td>10002</td>\n",
" <td>10000</td>\n",
" <td>1</td>\n",
" <td>scatter</td>\n",
" <td>3.00e-02</td>\n",
" <td>8.80e-04</td>\n",
" </tr>\n",
" <tr>\n",
" <th>4</th>\n",
" <td>10003</td>\n",
" <td>0</td>\n",
" <td>10001</td>\n",
" <td>0</td>\n",
" <td>14</td>\n",
" <td>0</td>\n",
" <td>10002</td>\n",
" <td>10000</td>\n",
" <td>2</td>\n",
" <td>absorption</td>\n",
" <td>3.16e-04</td>\n",
" <td>2.15e-05</td>\n",
" </tr>\n",
" <tr>\n",
" <th>5</th>\n",
" <td>10003</td>\n",
" <td>0</td>\n",
" <td>10001</td>\n",
" <td>0</td>\n",
" <td>14</td>\n",
" <td>0</td>\n",
" <td>10002</td>\n",
" <td>10000</td>\n",
" <td>2</td>\n",
" <td>scatter</td>\n",
" <td>3.90e-02</td>\n",
" <td>1.25e-03</td>\n",
" </tr>\n",
" <tr>\n",
" <th>6</th>\n",
" <td>10003</td>\n",
" <td>0</td>\n",
" <td>10001</td>\n",
" <td>0</td>\n",
" <td>13</td>\n",
" <td>0</td>\n",
" <td>10002</td>\n",
" <td>10000</td>\n",
" <td>3</td>\n",
" <td>absorption</td>\n",
" <td>3.78e-04</td>\n",
" <td>1.45e-05</td>\n",
" </tr>\n",
" <tr>\n",
" <th>7</th>\n",
" <td>10003</td>\n",
" <td>0</td>\n",
" <td>10001</td>\n",
" <td>0</td>\n",
" <td>13</td>\n",
" <td>0</td>\n",
" <td>10002</td>\n",
" <td>10000</td>\n",
" <td>3</td>\n",
" <td>scatter</td>\n",
" <td>4.86e-02</td>\n",
" <td>1.24e-03</td>\n",
" </tr>\n",
" <tr>\n",
" <th>8</th>\n",
" <td>10003</td>\n",
" <td>0</td>\n",
" <td>10001</td>\n",
" <td>0</td>\n",
" <td>12</td>\n",
" <td>0</td>\n",
" <td>10002</td>\n",
" <td>10000</td>\n",
" <td>4</td>\n",
" <td>absorption</td>\n",
" <td>4.21e-04</td>\n",
" <td>2.14e-05</td>\n",
" </tr>\n",
" <tr>\n",
" <th>9</th>\n",
" <td>10003</td>\n",
" <td>0</td>\n",
" <td>10001</td>\n",
" <td>0</td>\n",
" <td>12</td>\n",
" <td>0</td>\n",
" <td>10002</td>\n",
" <td>10000</td>\n",
" <td>4</td>\n",
" <td>scatter</td>\n",
" <td>5.52e-02</td>\n",
" <td>9.85e-04</td>\n",
" </tr>\n",
" <tr>\n",
" <th>10</th>\n",
" <td>10003</td>\n",
" <td>0</td>\n",
" <td>10001</td>\n",
" <td>0</td>\n",
" <td>11</td>\n",
" <td>0</td>\n",
" <td>10002</td>\n",
" <td>10000</td>\n",
" <td>5</td>\n",
" <td>absorption</td>\n",
" <td>4.86e-04</td>\n",
" <td>2.62e-05</td>\n",
" </tr>\n",
" <tr>\n",
" <th>11</th>\n",
" <td>10003</td>\n",
" <td>0</td>\n",
" <td>10001</td>\n",
" <td>0</td>\n",
" <td>11</td>\n",
" <td>0</td>\n",
" <td>10002</td>\n",
" <td>10000</td>\n",
" <td>5</td>\n",
" <td>scatter</td>\n",
" <td>6.30e-02</td>\n",
" <td>1.35e-03</td>\n",
" </tr>\n",
" <tr>\n",
" <th>12</th>\n",
" <td>10003</td>\n",
" <td>0</td>\n",
" <td>10001</td>\n",
" <td>0</td>\n",
" <td>10</td>\n",
" <td>0</td>\n",
" <td>10002</td>\n",
" <td>10000</td>\n",
" <td>6</td>\n",
" <td>absorption</td>\n",
" <td>5.30e-04</td>\n",
" <td>1.92e-05</td>\n",
" </tr>\n",
" <tr>\n",
" <th>13</th>\n",
" <td>10003</td>\n",
" <td>0</td>\n",
" <td>10001</td>\n",
" <td>0</td>\n",
" <td>10</td>\n",
" <td>0</td>\n",
" <td>10002</td>\n",
" <td>10000</td>\n",
" <td>6</td>\n",
" <td>scatter</td>\n",
" <td>6.93e-02</td>\n",
" <td>1.30e-03</td>\n",
" </tr>\n",
" <tr>\n",
" <th>14</th>\n",
" <td>10003</td>\n",
" <td>0</td>\n",
" <td>10001</td>\n",
" <td>0</td>\n",
" <td>9</td>\n",
" <td>0</td>\n",
" <td>10002</td>\n",
" <td>10000</td>\n",
" <td>7</td>\n",
" <td>absorption</td>\n",
" <td>5.86e-04</td>\n",
" <td>2.02e-05</td>\n",
" </tr>\n",
" <tr>\n",
" <th>15</th>\n",
" <td>10003</td>\n",
" <td>0</td>\n",
" <td>10001</td>\n",
" <td>0</td>\n",
" <td>9</td>\n",
" <td>0</td>\n",
" <td>10002</td>\n",
" <td>10000</td>\n",
" <td>7</td>\n",
" <td>scatter</td>\n",
" <td>7.57e-02</td>\n",
" <td>1.40e-03</td>\n",
" </tr>\n",
" <tr>\n",
" <th>16</th>\n",
" <td>10003</td>\n",
" <td>0</td>\n",
" <td>10001</td>\n",
" <td>0</td>\n",
" <td>8</td>\n",
" <td>0</td>\n",
" <td>10002</td>\n",
" <td>10000</td>\n",
" <td>8</td>\n",
" <td>absorption</td>\n",
" <td>6.30e-04</td>\n",
" <td>2.35e-05</td>\n",
" </tr>\n",
" <tr>\n",
" <th>17</th>\n",
" <td>10003</td>\n",
" <td>0</td>\n",
" <td>10001</td>\n",
" <td>0</td>\n",
" <td>8</td>\n",
" <td>0</td>\n",
" <td>10002</td>\n",
" <td>10000</td>\n",
" <td>8</td>\n",
" <td>scatter</td>\n",
" <td>8.09e-02</td>\n",
" <td>1.49e-03</td>\n",
" </tr>\n",
" <tr>\n",
" <th>18</th>\n",
" <td>10003</td>\n",
" <td>0</td>\n",
" <td>10001</td>\n",
" <td>0</td>\n",
" <td>7</td>\n",
" <td>0</td>\n",
" <td>10002</td>\n",
" <td>10000</td>\n",
" <td>9</td>\n",
" <td>absorption</td>\n",
" <td>7.10e-04</td>\n",
" <td>2.23e-05</td>\n",
" </tr>\n",
" <tr>\n",
" <th>19</th>\n",
" <td>10003</td>\n",
" <td>0</td>\n",
" <td>10001</td>\n",
" <td>0</td>\n",
" <td>7</td>\n",
" <td>0</td>\n",
" <td>10002</td>\n",
" <td>10000</td>\n",
" <td>9</td>\n",
" <td>scatter</td>\n",
" <td>8.94e-02</td>\n",
" <td>1.37e-03</td>\n",
" </tr>\n",
" </tbody>\n",
"</table>\n",
"</div>"
],
"text/plain": [
" level 1 level 2 level 3 distribcell score \\\n",
" cell univ lat cell univ \n",
" id id id x y z id id \n",
"0 10003 0 10001 0 16 0 10002 10000 0 absorption \n",
"1 10003 0 10001 0 16 0 10002 10000 0 scatter \n",
"2 10003 0 10001 0 15 0 10002 10000 1 absorption \n",
"3 10003 0 10001 0 15 0 10002 10000 1 scatter \n",
"4 10003 0 10001 0 14 0 10002 10000 2 absorption \n",
"5 10003 0 10001 0 14 0 10002 10000 2 scatter \n",
"6 10003 0 10001 0 13 0 10002 10000 3 absorption \n",
"7 10003 0 10001 0 13 0 10002 10000 3 scatter \n",
"8 10003 0 10001 0 12 0 10002 10000 4 absorption \n",
"9 10003 0 10001 0 12 0 10002 10000 4 scatter \n",
"10 10003 0 10001 0 11 0 10002 10000 5 absorption \n",
"11 10003 0 10001 0 11 0 10002 10000 5 scatter \n",
"12 10003 0 10001 0 10 0 10002 10000 6 absorption \n",
"13 10003 0 10001 0 10 0 10002 10000 6 scatter \n",
"14 10003 0 10001 0 9 0 10002 10000 7 absorption \n",
"15 10003 0 10001 0 9 0 10002 10000 7 scatter \n",
"16 10003 0 10001 0 8 0 10002 10000 8 absorption \n",
"17 10003 0 10001 0 8 0 10002 10000 8 scatter \n",
"18 10003 0 10001 0 7 0 10002 10000 9 absorption \n",
"19 10003 0 10001 0 7 0 10002 10000 9 scatter \n",
"\n",
" mean std. dev. \n",
" \n",
" \n",
"0 1.30e-04 8.67e-06 \n",
"1 1.98e-02 6.50e-04 \n",
"2 2.24e-04 1.44e-05 \n",
"3 3.00e-02 8.80e-04 \n",
"4 3.16e-04 2.15e-05 \n",
"5 3.90e-02 1.25e-03 \n",
"6 3.78e-04 1.45e-05 \n",
"7 4.86e-02 1.24e-03 \n",
"8 4.21e-04 2.14e-05 \n",
"9 5.52e-02 9.85e-04 \n",
"10 4.86e-04 2.62e-05 \n",
"11 6.30e-02 1.35e-03 \n",
"12 5.30e-04 1.92e-05 \n",
"13 6.93e-02 1.30e-03 \n",
"14 5.86e-04 2.02e-05 \n",
"15 7.57e-02 1.40e-03 \n",
"16 6.30e-04 2.35e-05 \n",
"17 8.09e-02 1.49e-03 \n",
"18 7.10e-04 2.23e-05 \n",
"19 8.94e-02 1.37e-03 "
]
},
"execution_count": 34,
"metadata": {},
"output_type": "execute_result"
}
],
"source": [
"# Get a pandas dataframe for the distribcell tally data\n",
"df = tally.get_pandas_dataframe(summary=su, nuclides=False)\n",
"\n",
"# Print the last twenty rows in the dataframe\n",
"df.head(20)"
]
},
{
"cell_type": "code",
"execution_count": 35,
"metadata": {
"collapsed": false
},
"outputs": [
{
"data": {
"text/html": [
"<div>\n",
"<table border=\"1\" class=\"dataframe\">\n",
" <thead>\n",
" <tr>\n",
" <th></th>\n",
" <th>mean</th>\n",
" <th>std. dev.</th>\n",
" </tr>\n",
" <tr>\n",
" <th></th>\n",
" <th></th>\n",
" <th></th>\n",
" </tr>\n",
" <tr>\n",
" <th></th>\n",
" <th></th>\n",
" <th></th>\n",
" </tr>\n",
" </thead>\n",
" <tbody>\n",
" <tr>\n",
" <th>count</th>\n",
" <td>2.89e+02</td>\n",
" <td>2.89e+02</td>\n",
" </tr>\n",
" <tr>\n",
" <th>mean</th>\n",
" <td>4.15e-04</td>\n",
" <td>1.71e-05</td>\n",
" </tr>\n",
" <tr>\n",
" <th>std</th>\n",
" <td>2.41e-04</td>\n",
" <td>6.82e-06</td>\n",
" </tr>\n",
" <tr>\n",
" <th>min</th>\n",
" <td>1.78e-05</td>\n",
" <td>2.81e-06</td>\n",
" </tr>\n",
" <tr>\n",
" <th>25%</th>\n",
" <td>2.06e-04</td>\n",
" <td>1.16e-05</td>\n",
" </tr>\n",
" <tr>\n",
" <th>50%</th>\n",
" <td>4.03e-04</td>\n",
" <td>1.71e-05</td>\n",
" </tr>\n",
" <tr>\n",
" <th>75%</th>\n",
" <td>6.05e-04</td>\n",
" <td>2.19e-05</td>\n",
" </tr>\n",
" <tr>\n",
" <th>max</th>\n",
" <td>9.35e-04</td>\n",
" <td>4.54e-05</td>\n",
" </tr>\n",
" </tbody>\n",
"</table>\n",
"</div>"
],
"text/plain": [
" mean std. dev.\n",
" \n",
" \n",
"count 2.89e+02 2.89e+02\n",
"mean 4.15e-04 1.71e-05\n",
"std 2.41e-04 6.82e-06\n",
"min 1.78e-05 2.81e-06\n",
"25% 2.06e-04 1.16e-05\n",
"50% 4.03e-04 1.71e-05\n",
"75% 6.05e-04 2.19e-05\n",
"max 9.35e-04 4.54e-05"
]
},
"execution_count": 35,
"metadata": {},
"output_type": "execute_result"
}
],
"source": [
"# Show summary statistics for absorption distribcell tally data\n",
"absorption = df[df['score'] == 'absorption']\n",
"absorption[['mean', 'std. dev.']].dropna().describe()\n",
"\n",
"# Note that the maximum standard deviation does indeed\n",
"# meet the 5e-4 threshold set by the tally trigger"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"Perform a statistical test comparing the tally sample distributions for two categories of fuel pins."
]
},
{
"cell_type": "code",
"execution_count": 36,
"metadata": {
"collapsed": false
},
"outputs": [
{
"name": "stdout",
"output_type": "stream",
"text": [
"Mann-Whitney Test p-value: 1.39844745394e-41\n"
]
}
],
"source": [
"# Extract tally data from pins in the pins divided along y=x diagonal \n",
"multi_index = ('level 2', 'lat',)\n",
"lower = df[df[multi_index + ('x',)] + df[multi_index + ('y',)] < 16]\n",
"upper = df[df[multi_index + ('x',)] + df[multi_index + ('y',)] > 16]\n",
"lower = lower[lower['score'] == 'absorption']\n",
"upper = upper[upper['score'] == 'absorption']\n",
"\n",
"# Perform non-parametric Mann-Whitney U Test to see if the \n",
"# absorption rates (may) come from same sampling distribution\n",
"u, p = scipy.stats.mannwhitneyu(lower['mean'], upper['mean'])\n",
"print('Mann-Whitney Test p-value: {0}'.format(p))"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"Note that the symmetry implied by the y=x diagonal ensures that the two sampling distributions are identical. Indeed, as illustrated by the test above, for any reasonable significance level (*e.g.*, $\\alpha$=0.05) one would **not reject** the null hypothesis that the two sampling distributions are identical.\n",
"\n",
"Next, perform the same test but with two groupings of pins which are not symmetrically identical to one another."
]
},
{
"cell_type": "code",
"execution_count": 37,
"metadata": {
"collapsed": false
},
"outputs": [
{
"name": "stdout",
"output_type": "stream",
"text": [
"Mann-Whitney Test p-value: 0.902458041178\n"
]
}
],
"source": [
"# Extract tally data from pins in the pins divided along y=-x diagonal\n",
"multi_index = ('level 2', 'lat',)\n",
"lower = df[df[multi_index + ('x',)] > df[multi_index + ('y',)]]\n",
"upper = df[df[multi_index + ('x',)] < df[multi_index + ('y',)]]\n",
"lower = lower[lower['score'] == 'absorption']\n",
"upper = upper[upper['score'] == 'absorption']\n",
"\n",
"# Perform non-parametric Mann-Whitney U Test to see if the \n",
"# absorption rates (may) come from same sampling distribution\n",
"u, p = scipy.stats.mannwhitneyu(lower['mean'], upper['mean'])\n",
"print('Mann-Whitney Test p-value: {0}'.format(p))"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"Note that the asymmetry implied by the y=-x diagonal ensures that the two sampling distributions are *not* identical. Indeed, as illustrated by the test above, for any reasonable significance level (*e.g.*, $\\alpha$=0.05) one would **reject** the null hypothesis that the two sampling distributions are identical."
]
},
{
"cell_type": "code",
"execution_count": 38,
"metadata": {
"collapsed": false
},
"outputs": [
{
"name": "stderr",
"output_type": "stream",
"text": [
"/home/wboyd/anaconda2/lib/python2.7/site-packages/ipykernel/__main__.py:4: SettingWithCopyWarning: \n",
"A value is trying to be set on a copy of a slice from a DataFrame.\n",
"Try using .loc[row_indexer,col_indexer] = value instead\n",
"\n",
"See the caveats in the documentation: http://pandas.pydata.org/pandas-docs/stable/indexing.html#indexing-view-versus-copy\n"
]
},
{
"data": {
"text/plain": [
"<matplotlib.axes._subplots.AxesSubplot at 0x7f2a2aeb9990>"
]
},
"execution_count": 38,
"metadata": {},
"output_type": "execute_result"
},
{
"data": {
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAZoAAAEZCAYAAACuIuMVAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAIABJREFUeJzt3XuYXFWd7//3N52uS98Tp0VuSTAQQG5JzoBRx0NQGAFH\nmefnBTPHx1EzGAYjHAeVwBwPaA4ekINIfgwTkskh4hCS0fE6o9PI2HHEUTpyEbWDgpgYEOjGC4KG\nJJDv+WOv3bWreld3dXdVV1f35/U89aRr19p7r6pO72+ttb5rbXN3REREamVWvSsgIiLTmwKNiIjU\nlAKNiIjUlAKNiIjUlAKNiIjUlAKNiIjUlAKNyCQws8vNbEO96yFSDwo00rDM7E/M7Dtm9lsze9rM\nvm1m/2WCx/xLM/t2ybZbzezjEzmuu/9vd3/fRI5RjpkdNLNnzex3ZrbHzK43M6tw39PNbE8t6iUS\nm13vCoiMh5m1A18FVgGfAzLAa4F9Ez00UNVZzGbW5O4vVvOYJRw42d1/bmYvB/4D6Ac2VVI9qvx+\nRUqpRSONahHg7v5PHtnn7ne5+4/iAmZ2gZn1h2/6PzKzxWH7ZWb2SGL7n4ftxwF/D7wqtBB+bWYX\nAP8N+Ego/+VQ9lAz+7yZDZjZz8zsA4nzXmlmnzOzz5rZb4G/DNs+G16fH1oh7zKz3eEYVyT2z5nZ\nZ8L5f2xmHx6l1WHhgbs/CnwHWJw43rsTn8MjZva+sL0F+BpwWKJF9DKLrAllB81sq5l1hX2y4X09\nbWa/MbN7zKx7vL9EmRkUaKRR/RR40cw2m9nZ8YUwZmZvA/4n8E537wDeDPwqvPwI8Jqw/WPAP5rZ\nIe7+EHAh8F13b3f3ue6+Ebgd+KS7d7j7eaFb6qvA/cChwOuBS8zsrEQV3gz8k7t3AVvCttKWw2uA\nY4Azgf9pZseG7VcB84AFwFnAO1P2TRWC5WuBhxObnwLODe/3PcANZrbY3f8AnAP8MrzfDnd/Erg4\n1P+1wGHAb4Cbw7H+EugADgfmhs9rbyV1k5lLgUYakrs/C/wJcBDYAAyY2ZcT365XEgWH+0L5R919\nT/j5n939qfDz54guyqeN4fSnAn/k7le7+4vuvgv4B+AdiTLfdfevhnM8n/YWgKvcfb+7Pwj8ADgl\nvPY24Gp3/527/xJYV0Gd7jOz54i6zHqJWmaE83891BF3/zZwJ1EQKWcV8Lfu/oS7HwA+DrzVzGYB\nB4CXAItCS/J+d3+ugvrJDKZAIw3L3X/i7u9193nAiUTfvj8dXj4S+FnafqHL6v7Q9fMb4ATgj8Zw\n6vnA4aFr69fhGJcDL02UqWSA/anEz38A2sLPhwGPjfFYS9y9DXg78EqgNX7BzM4xs++a2a9CXc9h\n5Pc7H/hi/P6IgtcB4BDgs0APsNXMHjOza8ysqYL6yQymQCPTgrv/FNhMFHAgujgvLC1nZvOIWkAX\nufscd58D/JgwxkF6F1Xptj3Ao6FrbW44Tqe7v2mEfcbiCeCIxPN5FewTj9F8HvgecCWAmWWAzwOf\nBLrD+/06I7/fXwDnlLy/1tDCecHd17r7CcCrgTcB7xr7W5SZRIFGGpKZHWtmf2Nmh4fnRwIrgO+G\nIv8AfMjMlobXF4YyrUTdbU+b2Swzew+F4ARRK+MIM2su2fbyxPM+4Fkz+0gYuG8ysxPM7I/H8hZG\neO2fgMvNrCu8v/eP4bgA1wAXmNlLibLxMsDT7n7QzM4B/jRR9ingJWbWkdh2C/CJEJQxs24ze3P4\nebmZnRi60Z4jaukcHGP9ZIZRoJFG9SxRF9E9ZvYs8J/Ag8CHYOib/dXAFjP7HfBFYK677wSuJ/rW\n/yRRt9ndieN+k6iF86SZDYRtm4ATQlfSF9z9IPBnRJldPwcGgI1Eg+SVKm1JJJ9/HHg8HPtOovTt\nkdK2i44VMu++BXw4jJ9cAnwudIO9A/hyouxPgDuAR8P7exlwYyhzp5k9Q/TZxmNYLyNqIT1D9Dn1\nEnWniZRl9b7xmZmdTdSvPgvY5O7Xlrx+LHArsBS4wt0/Vem+ItOBmV0InO/uZ9S7LiLjUdcWTWh+\n3wS8geib5YqQnpn0K+ADwHXj2Fek4YS5LK8O81mOBS4FvlDveomMV727zk4DHnb33SGNcitwXrKA\nuz/t7vcCL4x1X5EGlSEaJ/kdcBdRt9/fj7iHyBRW7yVoDqc4dfMxKp/PMJF9RaYsd/8FcFK96yFS\nLfVu0YiIyDRX7xbN4xTPETgibKvqvmamRQNFRMbB3StaCXwk9W7R7ACODosMZohSL78yQvnkGx7T\nvu4+5R9XXnll3eugeqqOqqfqGT+qpa4tGnd/0cxWE80ViFOUd5rZquhl32BmhwDfB9qBg2Z2CfAK\nd38ubd86vRURESmj3l1nuPu/AceWbLsl8fNTROtWVbSviIhMLfXuOpOE5cuX17sKFVE9q6cR6giq\nZ7U1Sj2rpe4rA0wGM/OZ8D5FRKrJzPBpkAwgIiLTnAKNiIjUlAKNiIjUlAKNiIjUlAKNiIjUlAKN\niIjUlAKNiIjUlAKNiIjUlAKNiIjUlAKNiIjUlAKNiIjUlAKNiIjUlAKNiIjUlAKNiIjUlAKNiIjU\nlAKNiIjUlAKNiIjUlAKNiIjUlAKNiIjUlAKNiIjUlAKNiIjUlAKNiIjUlAKNiIjUlAKNiIjUlAKN\niIjUlAKNiIjUVN0DjZmdbWYPmdlPzeyyMmXWmdnDZvaAmS1ObP+gmf3IzB40s9vNLDN5NRcRkUrU\nNdCY2SzgJuANwAnACjM7rqTMOcBCdz8GWAWsD9sPAz4ALHX3k4HZwDsmsfoiIlKBerdoTgMedvfd\n7n4A2AqcV1LmPOA2AHe/B+g0s0PCa01Aq5nNBlqAX05OtUVEpFL1DjSHA3sSzx8L20Yq8zhwuLv/\nErge+EXY9lt3v6uGdRURkXGod6AZNzPrImrtzAcOA9rM7C/qWysRESk1u87nfxyYl3h+RNhWWubI\nlDJnAo+6+68BzOwLwKuBLWknuuqqq4Z+Xr58OcuXL59YzSfR4OAgu3btYsGCBXR3d9e7OiIyTW3f\nvp3t27dX/bjm7lU/aMUnN2sCfgK8HngC6ANWuPvORJlzgfe7+xvNbBnwaXdfZmanAZuAU4F9wK3A\nDnf/u5TzeD3f50Tcccc2Vq68iExmAfv372LTpptZseL8eldLRGYAM8PdbcLHqfcF2MzOBm4k6sbb\n5O7XmNkqwN19QyhzE3A28HvgPe5+X9h+JVGm2QHgfuCvQlJB6TkaMtAMDg4yf/5x7N3bC5wMPEg+\nfwa7dz+klo2I1Ny0CTSToVEDzY4dOzjrrAt55pl7h7Z1dCzlrrtu4dRTT61jzURkJqhWoGnYZICZ\nYMGCqLsMHgxbHuTAgd0sWLCgfpUSERkjBZoprLu7m02bbiafP4OOjqXk82ewadPN6jYTkYairrMG\noKwzEakHjdGMQaMHGhGRetAYjYiINAQFGhERqSkFGhERqSkFGhERqSkFGhERqSkFGhERqSkFGhER\nqSkFmgY0ODjIjh07GBwcrHdVRERGpUDTYO64Yxvz5x/HWWddyPz5x3HHHdvqXSURkRFpZYAGotsG\niMhk0soAM9CuXbvIZBYQBRmAk2luns+uXbvqVykRkVEo0DQQ3TZARBqRAk0D0W0DRKQRaYymAem2\nASIyGXSbgDGYboFGRGQyKBlAREQaggKNiIjUlAKNiIjUlAKNiIjUlAKNiIjUlALNNKGFNkVkqlKg\nmQa00KaITGWaR9PgtNCmiNSK5tEIoIU2RWTqU6BpcFpoU0SmuroHGjM728weMrOfmtllZcqsM7OH\nzewBM1uc2N5pZp8zs51m9mMze+Xk1Xxq0EKbIjLV1XWMxsxmAT8FXg/8EtgBvMPdH0qUOQdY7e5v\nDIHkRndfFl7bDHzL3W81s9lAi7v/LuU803aMJqaFNkWk2qbLGM1pwMPuvtvdDwBbgfNKypwH3Abg\n7vcAnWZ2iJl1AK9191vDay+kBZmZoru7m1NPPRVAac4iMqXUO9AcDuxJPH8sbBupzONh21HA02Z2\nq5ndZ2YbzCxf09pOcUpzFpGpaHa9KzABs4GlwPvd/ftm9mlgDXBlWuGrrrpq6Ofly5ezfPnySaji\n5BkcHGTlyovYu7eXvXujNOeVK8/gzDNfp640EanI9u3b2b59e9WPW+9A8zgwL/H8iLCttMyRZcrs\ncffvh58/D6QmE0BxoJlO4rGZ3/zmN2QyC0KQgWSaswKNiFSi9Ev4xz72saoct96BZgdwtJnNB54A\n3gGsKCnzFeD9wDYzWwb81t2fAjCzPWa2yN3jhIL+yat6/d1xxzZWrryITGYB+/Y9ysGDTpTmHLVo\nlOYsIlNBXQONu79oZquBO4nGiza5+04zWxW97Bvc/Wtmdq6ZPQL8HnhP4hAXA7ebWTPwaMlr01pa\nV1lT06vJ5ZaTySzgwIHdSnMWkSlBS9A0qB07dnDWWRfyzDP3JrYuIpN5go9+9HJWrboAQCnPIjJu\n0yW9WcYoXqW5ra1t2IoA8Cv27/8qn/jE9XzhC19SBpqITAlq0TSQ5JjM/v27WLnynWzceBv79nUD\nvwJuBs6nre0kDhx4jH37voUW2hSR8VKLZoZJjsk888y97N3by6ZN/8i///u/kM0+BfwzcD5REsAv\nyGSOQgttishUoEDTIMqt0pzJZLj11g3k828ZWuvsxhv/Dy+8sBsttCkiU0G905ulQsWrNBenL596\n6qmceebrigb+Ozo6WLnyDJqb5ysDTUTqSmM0DSQeo0kGjxUrzi9bfufOnfT19XHaaadx/PHHT2JN\nRWQ6qNYYjQJNg6l0lebSxIGRgpJWfhaRNAo0YzCdAk0lxnJ757EEJBGZWZR1JkPiuTXxrQEqvb1z\nWibbypUX6RYDIlJVCjQNLu3WAJXe3rnSgCQiMhEKNA2sXIsEqOj2zpUGJBGRiVB6cwOLWyRptwZY\nseL8YSnPpbq7u9m06WalQYtITSnQNLCR5tZUKg5I999/PwBLliypQU1FZCZT11kDi1skaV1kY7mt\n8113fZM///MVvP3tl2sBThGpOqU3TwOl82DGkt48lrIiMrNUK71ZXWfTQHd3d1FQGGnsJg5EcWAa\nrayIyESNq+vMzO6rdkWkekbKJivtUrvvvgeUeSYiNTWuQOPuS6tdEamecmM3wLB06A9+cA033HDN\nqKnQIiLjNWLXmZk1AXe5+xmTVB+pkrT05h07djB79uHAPmCQuJts6dLF7N79kNY7E5GaGDHQuPuL\nZnbQzDrd/ZnJqpRUR+nYzX33PcCzzz4CXADsAS4b6iYrLSsiUi2VJAM8B/zQzL4B/D7e6O4X16xW\nUnWDg4N88INrgO8RZ5fBMm644cbUAKMVnUWkWioJNF8ID2lg999/P7NmHUlyXbO2tkUsXbp4WFmt\n6Cwi1TTiPJowRnObu/+3yatS9U33eTSjueOObbz3vRfy/PP7ge+SbNGsX38jq1ZdMFRW82pEJDYp\ntwlw9xeB+WaWmeiJpD7ihTeff/5ywIBlwCnAGcBVfPCDa4puC6AVnUWk2irpOnsU+I6ZfYXiMZpP\n1axWUjW7du0KmWbXAhuAa4B/ABYA3TQ3by2anFmN9dNERJIqmUfzM+BfQtn2xEMaQBQ4dgNHAmcB\njwNZoJu0INLd3c3Kle8kavksApaxcuU71W0mIuNW8VpnZtbi7n+ocX1qYqaP0dxyy0YuvPASooyz\nncBfA3PJ538zbKC/MEbzz8B+4FFyucv5xS9+qmAjMsNM2q2czexVZtYPPBSen2JmN0/0xDJ5Vq26\ngPXrbySbPZ329mvI5Zy1a9/L7t0PDcsmK4zRPAWsADby/PP7ueWWjXWouYhMB6O2aMzsHuCtwFfc\nfUnY9iN3P7EqFTA7G/g0UdDb5O7XppRZB5xDNEb0bnd/IPHaLOD7wGPu/uYy55jRLZpYJXNjBgcH\nmTdvEc8/b8B2lHkmMnNNWosGwN33lGx6caInhqEgcRPwBuAEYIWZHVdS5hxgobsfA6wC1pcc5hKg\nvxr1me66u7s59dRTRwwW3d3dvOlNbwDmoswzEamGSgLNHjN7NeBm1mxmHyLq6K+G04CH3X23ux8A\ntgLnlZQ5D7gNwN3vATrN7BAAMzsCOJcojUqqYHBwkK9+tQf4FVrRWUSqoZJAcyHwfuBwopSlxeF5\nNRxOtOhW7LGwbaQyjyfK3AB8GFC/WMLg4CA7duwomh9TqV27dpHNvpyo4XgGsBR4FVdccam6zURk\nXEadR+PuTwNTbmUAM3sj8JS7P2Bmy4lmI5Z11VVXDf28fPlyli9fXsvq1c1El48pzKM5nij/4xvk\ncu8vWj1ARKan7du3s3379qoft663cjazZcBV7n52eL4G8GRCgJmtB3rdfVt4/hBwOtHYzDuBF4A8\n0dyeL7j7u1LOMyOSAaq1fEwcrJqb53PgwG6tdSYyQ02XWznvAI42s/nAE8A7iHJqk75C1FW3LQSm\n37r7U8AV4YGZnQ5cmhZkZpJq3Zb5zDNfx5e+dAcAS5YsGXb7Z3WhichY1DXQhPvdrAbupJDevNPM\nVkUv+wZ3/5qZnWtmjxClN7+nnnWeyoYvH7Odfft+RltbW8XHSOt6A7Sas4iM27i6zsxsqbvfV4P6\n1MRM6TqDQqBw7+T5558knz8aeJwbbriGpUsXjzqHprTrLZc7HbNZWs1ZZAaa1Hk0Kf56oieW2lix\n4nzuvfdu3H8DfI+9ex9k797LuPDCS3j961cxf/5x3HHHttR901Zubmp66bD72GhOjYiMxbgCjbsr\nBWkKe+6558jlXk4UHAaJVm7+Hs8+ex979/aycuVFqanPxV1vAA/y4osDHDy4B82pEZHxKjtGY2ZL\nR9qxkbrOZprigLGPaOXm0RME7rrrm7zwwn7gVcChZDJP83//7y0ArFx5BrNnz2P//l3ccMMn1W0m\nIhUrO0ZjZr0j7Ofu/rraVKn6ZtIYTSweq2lqOoznnvsZ0crN5cdYisdnmoHPk83+H/bseYTu7m5u\nuWUjl1zyITKZ+bzwwuNKCBCZAao1RlPXeTSTZSYGGoiCx/33309v73ZuvPGWEefF7Nixg7POupBn\nnvkIcBHRjdEeYu3av2XVqgt0e2eRGWjS5tGYWQvwN8A8d3+fmR0DHOvu/zLRk0tt3XXXN4fSkt0P\n8uEPv5VVqy4YCg7JuTELFixg375HifI8thMHlE984gxOO+2Py87PASqaX1PteTjJ41VaBxGpE3cf\n8QFsAz4C/Cg8bwEeGG2/qfSI3ubMMjAw4Pn8XIcfOLjDDzyfn+sDAwPu7r5ly1bP5+d6Z+dSz+fn\n+pYtW33t2qsdjnYYcOhx6PG2thO9p6cn9Vjr128Ydow0aeeaiOTxmpvbPZPprNqxRaQgXDsnfg0e\ntQB8P/x7f2LbD6px8sl6zMRA09fX552dS0NgiB4dHUu8r6+vbBDq7+/35uZ2hzaHlhB08r5mzeX+\n3vde4JB3OMYh7ytX/tWIgSw2WsAbq+LjDTjMKTp2Ntvl/f39ReXj9ywiY1OtQFPJygD7zSxPWCHZ\nzBYSpTLJFDZ8lYBCWnK5pWr27NmD+0GgCfjO0H7XXPMq4CDRQtl/AjzKZz/7YWbNWkAymw0OG5bN\nVq1lcdKPtwM4qqgO+/Z1s2TJMm69dQOgFQ1EpoTRIhFwFvAtogkZtwO7gOXViHKT9WAGtmjcC11M\nHR1LirqVyrUyLrvscodmh4VFLSE42eF2hw6HLoelDrnQwikcA/LDWhPlut1q1aKBuQ69nst1VfW8\nIjMRk9F1RrT0/pHAS4A3An8G/FE1TjyZj5kaaNzLdx3FQai19eSh8ZZcrsuhM+XiPSdc1BeGgOMO\nH3DIhi62hQ6dns0u8M2bN/vAwEDJOEqbZzKdwwLeeCUDaHNzmzc3d4QuvbkOWx3cW1sXeWvrKald\nhyJSmWoFmhG7ztzdzexr7n4S8K/Vb09JrXV3d5ftooq6yfbhfpDBwUGamg4BWolyP84A5gM/Af6W\naHHtJ4gauHHjtg84FPgGcBH79j3JBz5wPRde+N85eNDZv/8/QhfXg2Qyr+Wyy97K6aefTiaTYXBw\ncNwZYitWnM+ZZ75uKNPs6aefZsmSZezb92VgOfGKBtGdwod3HYrI5Bp1Ho2ZfQa4yd13TE6Vqm+m\nzqMpZ3BwkHnzFvH8839HFDieIJc7HYDnnzei9OY4gKwMPz8D7Ccau9kHXAA8kDjq0cA5wP8PbAGu\nBB5OvH4MUa9r89BCn5s23VwUMEoDz86dO+nr6+O0007j+OOPH/E9pd1DB9B9dUQmoFrzaCrpdnqI\n6OZiPyP6evhD4MFqNKcm68EM7jpLE6Uxt4Sxlqi7qaNjia9de7U3N8cZZwt89uxWv/LKj3k22+HQ\nG7ql5jgsSBmfaQnda+7QH54nX88N29bc3OG5XNdQavL69RuGuvlWr74knGORQ95Xr754qP6l3YHx\n8/7+/mHdhP39/b558+aisSMRqQyTmN48P+1RjZNP1kOBpiAtEQDmeC7X5QMDAz4wMOCXXXa5Z7Nd\n3t4ejamsXn3x0JhILtfla9de7evXb/BstisxNtIWjhUHo5eFQHF0+Hd1CGyeeCwM293hWoe8t7cv\nCYEtPdFgy5atnst1eWvrsZ7Ldfnq1Zd4Pj837Nfl69dvGHqv1Z6/IzLTTFqgmQ4PBZqCtPk1sNDX\nrr3a3ctnpJVrLRS3djpLWi29HmWxvSK0doYHuCiLrT/8e3so97YQwJJ1PMbXrVsX5vnMCUGrMzUg\nrV+/oaL5O5pjIzIyBRoFmnEZ7QLc19fn+fxJwwLRmjVXpB6v0MJY5M3NrZ7NnlC0by53XCIYbAjB\nYaEXMsROcniLF3fltaQGkI0bN5YEsts9LRU7m+3wnp6eshNW43qrtSMysmoFmvHe+EwaVHd3N5s2\n3Uw+fwYdHUvJ589g06abefrpp/nMZz7DU089xd69j5C8/ww8wTXXfIpbbtnI4OAgO3bsKLqfTZTd\nlaepqTlx75pBoqSAJ7nuuv9FU9Orgf8OzAGeBK4BjidKGPg68F3gXqA3HPWdwDJgEbCMN7/5DXz7\n23cTJSbEEzSjRIbiuj5Gc/M8gGH31omzzgYHB1m58iL27u3lmWfuHfEePSJSBdWIVlP9gVo0wyS7\njUoH3guPhaGbaquXG7wfPt6T8cL8mqM9k+n06667PoznlCYP5BxWevmxmw6HIx2y3tQUt3JK5/hk\nwvaTQ2vo2qEWWrkJq4XuwwGHPocBzbERSYG6zhRoqqG/vz+1mwo+6XCsFzLJ4gBw+1C5bLbD29uX\nlASJRSXHuzYEnkUpwaTLoySCzqLzz57dHoJQsk7tDktC0Jsbfm7xTOaw0CXX4XCiQ8vQeJN7YXWC\nnp4eHxgY8P7+fl+3bp3Pnt3mhbGeOd7c3Ja6TpvGcGQmq1agqWStM5nG+vr6KL0DJxwOnELU/fUE\n0E3chRZ1V0XlmpvnsX//z0lOioTdwLEU30b634C3lJT7JVE32BNE83L+BHgp8ATuTi53NM8/n6zT\nkcDPgUOAO4BHyeUuB/4AvBJ4hHjeTy6XG3p/yVslPPvsTg4eJJz3RZLruZn916LPJZ6Xo3XSqq/a\nt4yQBlCNaDXVH6hFU1b5Fk2vF1KVF3ouNyfMsUm/VUB7++Kw3xFeyC7rc4iXgYlbIseEcteG1tLt\nofsr4/D3YVtvSp3irrbCCtKrV1/s69dvCNteEVo1H/X0zLPk+0zWK3oku86qveL0dDTe1p6SMBoL\n6jpToKmWP/3Tczw55+WkkxYPmzcz0phHfNG57rrrE0Fkbrj4JwNGrzc1tTocErq75nqUdZZ3ONST\na5VlswvCGmbxONEGLx2fyefnek9Pj2cy80I33JJwjHkpmWebE913w1Otkynco2WsjWQmdLeNN1go\ngDceBRoFmqoo/PF/MVyMv1h23kxcfqQLaTyRs7X1RM9mO/z0019f1Aoxy3r6ygJzQ0sm/jfv+fyC\n0EKKB+2TF/8Bb2lZ5FdccUXKsfLe2npcYuXoXodtXjzuc+1QcM1mu4YmpXZ2LvVcrsubmvLhM+lz\n6K3oghhfgNvbT/JstqNo8uhEf0dTJXhNJFiMdI8kmZoUaBRoqqIWf/zJJWEKF/q+8O8cj+7eWdx1\nFbVG+kJAynqhay1uxQyEVkty9YGFIVi8xIuTFl7h2WxHmYy6jBe67/4/z2Y7/O677w4rV9/ucL1H\nyQmHemk3Xbn3Ga+oEL3XuDV3isddeBMx1bqaJvL/RS2axqNAo0BTFRP54x/tm3b6KgRLQqApTYuO\nAlIm0+mtrccnym/1KLtsvhdSmeNJm8mA0+VxGvZIqwM0N3f47Nkt3tJy9FCro7D22/FeGJ8aeXWE\n0gCwdu3V3t5+0rD9stmu1M+nklbKVLwwT7RO5bpfZWpSoFGgqZrx/PGP9E17eIsm7d42UddVLneC\nQ96z2SM9m+3w6667ftg+2WxHYh5Oj0fjOmlL2rR4JlPoskpb5SCfP9HXrLliaC23QjdZb2hRneDD\nu+mi/bLZjqGutUymOCU7l+vyWbMyHo1zFfZrb1887Nt+JZ9d/O9YWw8T6WardN+JBoup1BUoI1Og\nUaCpqrH88Y/0rTZ5Ec3luvxtbzt/6KKUyXR6c3Pb0AUqXq05ntCZtohn3FooXHDjAHP7sGDQ0nKS\n9/T0DNWzXEZd1E1W2iKa61G3WXqLprDdw7mTAWWrxytel56v9Nt+tD5cV2qZ4s9ujr/rXe9O1DX9\neEkT6WYb676THSxqfT4Fv3TTJtAAZxPdiuCnwGVlyqwjWqvkAWBx2HYE8E3gx0S3Lrh4hHNU4zOX\noNw37eLbNscX8aOLMtfSlvgfbRHP4WWu9bTbDiQvwgMDA75582bPZud5YYLnXM9kDg133kxrEXV5\nIT07DhpHe3Nzh+fzRyXe70Di3MlxpK0OrV6462je3/a284c+ty1btoZFSIsnrw7/7OL6RCndTU0t\no7YeJtoFOtW66JJqPU5VyfFnaiCaFoEGmEU0024+0BwCyXElZc4B/jX8/Erge+HnlyWCThvRrSCP\nK3Oe6nxtZNv/AAAVWUlEQVTq4u7lL0yFtOD09OG0QFNp91Bpd8369Rt8zZr4dgaLiy4QheyvJSFY\nfNTjZIRcrivUfXiLKAoOs73QjfZFh6xv3LgxJamhxaN5O0eF/ZLvOZ4flBtKSih8Zr2pn00UaEoX\nM10SjpP3bdu2jXiR6+npGfOtq+PfxUTSuWut1kGwXAtzpPG4mTSuNF0CzTLg64nna0pbNcB64PzE\n853AISnH+hLw+jLnmfAHLsXS+ukLF4XhF/H4xmqlf7BjuZCkfascvYXU65D11tbjhs65evXFntYi\nilom8ZI4UYJBvDRNcfZa1gtL4rR7lKSQFrhO9Hz+5b5582Zfs+YKL3S3xZNXo3lLq1dfXKabb24I\nWsf45s2bR/xdRF1s5Vt45X5/5cacpkqLppYp0eVamMnxuHx+bupE5anw2UyG6RJo3gJsSDx/J7Cu\npMxXgVcnnt8FLC0ps4DoPsFtZc4z8U9chkm78Je76BVaEuXHJqqRiVR8YYov6Is8k+ksyUTrdfjr\nUM94oudWj9KZM57JHOq5XFdKMEx2lcUBIR8CTmngitZsa2tbHAJbR0kA7PB43lJPT0/o5hu+mGl8\n07f4M0+u3VZct8KY01i72Zqb2yr+HVSaMTeWMb9yZWvVohmphVk8Hhd3YRbS59MSPKYrBZrC8zbg\n+8B5I5zHr7zyyqFHb2/vhH8BUt7AwMBQ6yV9QN+HfTMtd7EZa994JV1UUZfagEcZbMkbrkXpyHff\nfXfRraGLu5aGZ6RFwanH4eoQPBaGi1PGo/Gk+ILVFl6P70q6dehziAJNh0dZb1d7YaWDlqFxni1b\ntoaWx9EOLd7c3JbyuQ44HF72/kHu6a2E9vbF3tPTM+Jnnfy9jtSNNJaupkrKFt73wqH3PVIQreT/\nS/oXkmM8k+ksGY9zL11Mthrzo6aq3t7eomvldAk0y4B/SzyvpOvsobjrDJhNtGLjJaOcZ+K/ARmz\n9AmNlX8zrTQNOG2/kQbdo2Vw4pWbWz2emJnNdg0b54m7lqLlcEZq0fSG19odPu7R2mrJ7q/4gnWk\nR3cd/eKwz6GwhE9hnCd5i+20VPFstnNYZlq0GnX7iAEjrXvxuuuuH/F3UUn33Fi7QispOzAwkJhQ\n21/0uSTLrFlzhWcybd7eftKoLbqenp6Sz613aPJu2tyrqEU6/FYU0910CTRNiWSATEgGOL6kzLmJ\nZIBlcTJAeH4b8KkKzlOFj1wmaixdZJWmUKcdpxBohl8U+/v7E0EjvjC3eSbTVtQ9VXruaK5NvFJA\nPI6zJJyjKfz7Ei++U+jWcHHq80KLZq5HrZa8Z7MLhuofX8ibm6MVCXK5E4reW19f37DB/ujYh/pr\nX3u6D+8CXFiU6l2qsBjp/PBvVKe0b+rFY28neXwPH3BvazvRN2/ePObkjrGU7evrC5Nhrw7vb6kn\nbwexZcvWcIvv6B5I0e/mWs/luoa6GEv/D3Z2Lh2Wbh93k65de3W4a+zJQ63x6PyF9z1VkiVqbVoE\nmuh9cHbIGHsYWBO2rQLelyhzUwhIPwCWhG2vIVrr/QHgfuA+4Owy56jSxy4TNb6ujcJFKC0NuDSt\neaQxi/TVChYW3cNmeJl4nsxRHiUCfNTjFkc22+nbtm3zjRs3pgSwOQ65xO2siwfcm5s7vL+/3wcG\nBsKFMm5ldXpTU3GmWXqLJk5/znv0jbvwjR/yQ4Em7TPftm2bZzIv8+Jxo+ErGfT39/ull17q0bpz\ncYvrlHDB/wuH/ND8p/Ekd5TLXkyWj4Jiejp7f39/aJmUtjI7PFrz7pRR6xafLw5CUfZf3jOZwzyX\n60q9wZ9aNA0WaCbjoUDTeEZPoS4OQPG3y+FBYsBbWxcVXXRHu2hUNvi/wDOZzlFaHAu9ufmlfuml\nl4YAVfx6S8tJQ2NAwxMJWrylpXhQf8uWreGmcMlkgV6HbFistM0L3+qj1kla66+QQfeScKxCndra\nThn6LAvlkmNOye624iy5bLbLr7vu+pClFc0lmj27zdesuWLU7LfkhN70rMT0TMbNmzd7a+uxw14r\nHVfJ5+eGNPUTU//vpAfyaKwvTqevVsJKI1GgUaCZ9kZOoa6kRZMeSCrpwovLtLYu8tJlZVpaTvJ1\n69aNOjYBXZ7LdXl/f7/Pnh2PBQ2v18aNG70wFyfunlnoUYJBcf23bdvmhRWtC4PY0bhT8fFzuTme\nybR5chXq6Nt/PKaU1hKIAlR6unUy+6rPh981tXRB1A94POk1GZRLpY+ZVDY3q3yLpjhTLJM53qOx\nsawnM8oymc6hFl9p91i80GsyGM20SZsKNAo0M0K5FOqRAkUlgaTSFN3RuuqSZdeuvTp0gQ3Pjopa\nI4VVA+ILb5RR1REugHHX2ZyiC2Vpdl7UquktufCWLovjoR5HhIv9UQ5zffbsl4WAEGfPFd8aG1Z7\nJtPp69atC+VKjxe3EnpTAlEy269/WGDI5eZ4T0/PsFUfenp6fN26deFCXzhfa+uJfvXVV4dg2etw\nwVDggrzPmpUb+gyLW1HtIfAm65YNgTW+/9EhDp0+e3ar9/T0JBIx4m7Baz3ZoplJwSVJgUaBZkYb\nLVBU89tnpYEt7vJZs+byYeMMcZ3S58D0eun4TTTGMDAssEXjOW0eZbclWxTJZXHiY8wZOkZ00fxi\nuODGdzMtXcmgy+PWVNTKKg0k8eTUBaHs+aGeR3syXTtKUtjsw28FsdCz2Whpnnz+KG9ubg9B82gf\nng5euF9Q9NpsLyQtxHdSnTN0g7v+/n7ftm2bv/3t53sm0+HZ7IJwnhND8E9rnb3DocXz+WNTXs87\ndHv8hSFel2+mBRwFGgUamUQjzfMZ70BxYTwpbW7OwhBIWoYlKkRjQf1euD9PfHHMhgv+yeFCujVx\nvKNC+fhePIeE8mkTRFu8p6cnrKAQvx7dyyeTie+a2hrqHLfGeksu0kemXLw7vZD23eXDg2uUPBG9\n7zi5odytved6vEpDa+spnsl0+qxZyVtIuMcpy9dff72XjkVFz1tDQBueCh+9/nEvBOr8qGnT05EC\njQKNTAHVuRFYrw+fnT7H4bKhtdKG71OcVQf50NLp9WhsJ3kRL3eh7g3lrkgEoXZvasoNZbz19/eH\n7qvSjLrOcJ5eh4zncnOK1qFbu/bqkKBQ6OqCeIHTreHCXjwwH3XfHe7wZ16cJn61D+/GOzklwMXL\nAhXKxRNRo5ZTaYslrk/a5188xpNMU59JXWkKNAo0MgVMdImUuNstl4u6eqJB63y4iKYPohfPlM97\nU1MhwyxaILJ0rk/Wh4/fLAkXzvgC2uvRYHmhuyoeYyqf7h2NZzQ3HzlsVYHC5/JFj1oMyYAwx6PW\nRFo6eJw9V5xUUXwb7mSgSL6nE720lRena69fv8EzmU7P5U7wbLYrLMja4YXuvUJiRTbbNWx9s+TE\n25kyh8ZdgUaBRqaMiazVlhy36e/vL5t9Ndp4Tywa1I6/6cdjL9mUb/TxN/m8t7Wd6LlcV0qZOUNZ\nc6Ole8eTXWPRatLHetTqSesWbE5Mgo3Tp9scVntaN5dZJpHa3eJp85Kiem3wZCsvmeLd3n7S0F1V\n3ZOTVou72vr7+xP7LA5lCksJqUWjQKNAI3UxnuSDtDku1emKizOmogyr6667PjEZ8USPVh5YMNTN\nFc/lGT4PaIm3ti7yvr6+EdO98/kTi+oXly3M0i8NCHO9peUV3tPTE1a1PiwEpAFPS2qI158rLEMT\nj5u0hRbMQm9u7vCmpsJ9gOJW3mgtzvXrN6TeaiL5O52pc2jcFWgUaKShlbsApt3+ulyLpjSwFQep\naE5OW9uJRanRycVCR58HNKdoTbFK0r3TjxOP1bzCS9cKSyuftpJ0+ooOx3gm01Z0U73SVl4lgbvS\nVHdlnSnQKNBIQxnpAjjWdOr49WqMF1WySvJI9Ut7X62tJ/vFF1/i2WyHt7WdOGyfchNzR7sTazbb\nNazLrtREP5OZToFGgUYaWCUrHIwnnXqi9/YpN/aTVm6s9RupVVBJi2G8762a9zuaaRRoFGikwY3n\nAlitrqBaquWFfbzvrd6fSaOqVqCx6FjTm5n5THif0ngGBwfZtWsXCxYsoLu7u6Ly8+cfx969vcDJ\nwIPk82ewe/dDFe0/Wcb6vmRqMjPc3SZ8nJlwAVagkenkjju2sXLlRTQ3z+fAgd1s2nQzK1acX+9q\nyTSkQDMGCjQy3ajFIJNBgWYMFGhERMauWoFmVjUqIyIiUo4CjYiI1JQCjYiI1JQCjYiI1JQCjYiI\n1JQCjYiI1JQCjYiI1JQCjYiI1JQCjYiI1JQCjYiI1JQCjYiI1JQCjYiI1FTdA42ZnW1mD5nZT83s\nsjJl1pnZw2b2gJktHsu+IiJSX3UNNGY2C7gJeANwArDCzI4rKXMOsNDdjwFWAesr3VdEROqv3i2a\n04CH3X23ux8AtgLnlZQ5D7gNwN3vATrN7JAK9xURkTqrd6A5HNiTeP5Y2FZJmUr2FRGROqt3oBmP\nCd+ER0REJs/sOp//cWBe4vkRYVtpmSNTymQq2HfIVVddNfTz8uXLWb58+XjqKyIybW3fvp3t27dX\n/bh1vZWzmTUBPwFeDzwB9AEr3H1nosy5wPvd/Y1mtgz4tLsvq2TfxDF0K2cRkTGq1q2c69qicfcX\nzWw1cCdRN94md99pZquil32Du3/NzM41s0eA3wPvGWnfOr0VEREpo64tmsmiFo2IyNhVq0XTiMkA\nIiLSQBRoRESkphRoRESkphRoRESkphRoRESkphRoRESkphRoRESkphRoRESkphRoRESkphRoRESk\nphRoRESkphRoRESkphRoRESkphRoRESkphRoRESkphRoRESkphRoRESkphRoRESkphRoRESkphRo\nRESkphRoRESkphRoRESkphRoRESkphRoRESkphRoRESkphRoRESkphRoRESkphRoRESkpuoWaMxs\njpndaWY/MbMeM+ssU+5sM3vIzH5qZpcltn/SzHaa2QNm9s9m1jF5tRcRkUrVs0WzBrjL3Y8Fvglc\nXlrAzGYBNwFvAE4AVpjZceHlO4ET3H0x8HDa/o1m+/bt9a5CRVTP6mmEOoLqWW2NUs9qqWegOQ/4\nTPj5M8Cfp5Q5DXjY3Xe7+wFga9gPd7/L3Q+Gct8DjqhxfWuuUf7zqZ7V0wh1BNWz2hqlntVSz0Dz\nUnd/CsDdnwRemlLmcGBP4vljYVup9wJfr3oNRURkwmbX8uBm9g3gkOQmwIH/kVLcx3mOvwUOuPuW\n8ewvIiK1Ze7jur5P/MRmO4Hl7v6Umb0M6HX340vKLAOucvezw/M1gLv7teH5u4ELgNe5+74RzlWf\nNyki0uDc3SZ6jJq2aEbxFeDdwLXAXwJfTimzAzjazOYDTwDvAFZAlI0GfBj4ryMFGajOByUiIuNT\nzxbNXOCfgCOB3cDb3f23ZnYosNHd/yyUOxu4kWg8aZO7XxO2PwxkgF+FQ37P3S+a5LchIiKjqFug\nERGRmWHarAwwlSeAljtnSZl1ZvZwOP/isexb73qa2RFm9k0z+7GZ/dDMLp6K9Uy8NsvM7jOzr0zV\neppZp5l9Lvyf/LGZvXKK1vODZvYjM3vQzG43s0w96mhmx5rZf5rZ82b2N2PZdyrUc6r9DY30eYbX\nx/Y35O7T4kE01vOR8PNlwDUpZWYBjwDzgWbgAeC48NqZwKzw8zXA/65SvcqeM1HmHOBfw8+vJOoG\nrGjfKn5+E6nny4DF4ec24CdTsZ6J1z8I/CPwlRr+f5xQPYHNwHvCz7OBjqlWT+Aw4FEgE55vA95V\npzr+EfBfgLXA34xl3ylSz6n2N5Raz8TrY/obmjYtGqbuBNCy5yyp+22hHvcAnWZ2SIX7Vsu46+nu\nT7r7A2H7c8BO0uc71bWeEH1zBM4F/qFG9ZtwPUNr+rXufmt47QV3/91Uq2d4rQloNbPZQAvwy3rU\n0d2fdvd7gRfGuu9UqOdU+xsa4fMc19/QdAo0U3UCaCXnLFem0vpWw3jq+XhpGTNbACwG7ql6DdPr\nMNZ63kCUrVjrwcmJ1PMo4GkzuzV0T2wws/xUq6e7/xK4HvhF2PZbd7+rTnWsxb5jVZVzTZG/oZGM\n+W+ooQKNmX0j9AXHjx+Gf9+cUryRJ4A2ZDq2mbUBnwcuCd/KphQzeyPwVPjmaEzdz3k2sBT4O3df\nCvyBaG3AKcXMuoi+Cc8n6kZrM7O/qG+tGtt0/Ruq5zyaMXP3s8q9ZmZPhW6ceALoQEqxx4F5iedH\nhG3xMd5N1CR8XXVqPPo5E2WOTCmTqWDfaplIPQldJ58HPuvuaXOipkI93wq82czOBfJAu5nd5u7v\nmmL1BNjj7t8PP3+eaNyxFiZSzzOBR9391wBm9gXg1UC1v6RVUsda7DtWEzrXFPsbKuc1jOdvqBaD\nTfV4ECUDXBZ+LpcM0ERhECxDNAh2fHjtbODHwEuqXK+y50yUOZfCYOsyCoOto+47FeoZnt8GfGoS\nfs8TqmeizOnUNhlgop/nt4BF4ecrgWunWj2J+vp/COSIvtluBt5fjzomyl4JXDqefetZz7BtyvwN\njVTPxGsV/w3V9A1N5gOYC9xFlK1xJ9AVth8K/Eui3NmhzMPAmsT2h4kmjt4XHjdXsW7DzgmsAt6X\nKHNT+OX/AFg6Wn1r9BmOtZ5LwrbXAC+G/7D3h8/v7ClUz6Upx6hpoKnC7/0UopUxHgC+AHRO0Xpe\nSTRw/SBREk5zPepItKbiHuC3wK+Jxo3ayu1br8+yXD2n2t/QSJ9n4hgV/w1pwqaIiNRUQyUDiIhI\n41GgERGRmlKgERGRmlKgERGRmlKgERGRmlKgERGRmlKgERGRmlKgERGRmlKgEakSM5sfblR2q0U3\n4PtHM3u9md0dnv+xmbWY2SYz+56Z3Wtmb0rs+x9m9v3wWBa2n25mvYmboH22vu9SZOy0MoBIlZjZ\nfKIlPRa7e7+ZfR94wN3/KgSU9wL9wI/dfYtFd4HtI1oS3oGD7r7fzI4G7nD3U83sdOBLwCuAJ4Hv\nAB9y9/+c/HcoMj4NtXqzSAP4ubv3h59/DPx7+PlHwAKilXLfZGYfDtvjFbqfAG4Kt0l+ETgmccw+\nd38CwMweCMdRoJGGoUAjUl37Ej8fTDw/SPT39gLwFnd/OLmTmV0JPOnuJ5tZE7C3zDFfRH+30mA0\nRiNSXaPdCKoHuHiocNSCAegkatUAvItoKXeRaUGBRqS6vMzP8fO1QHN8h1jg4+G1m4F3m9n9wCLg\n9xUcX6QhKBlARERqSi0aERGpKQUaERGpKQUaERGpKQUaERGpKQUaERGpKQUaERGpKQUaERGpKQUa\nERGpqf8HTKaAIu8FxywAAAAASUVORK5CYII=\n",
"text/plain": [
"<matplotlib.figure.Figure at 0x7f2a32c74fd0>"
]
},
"metadata": {},
"output_type": "display_data"
}
],
"source": [
"# Extract the scatter tally data from pandas\n",
"scatter = df[df['score'] == 'scatter']\n",
"\n",
"scatter['rel. err.'] = scatter['std. dev.'] / scatter['mean']\n",
"\n",
"# Show a scatter plot of the mean vs. the std. dev.\n",
"scatter.plot(kind='scatter', x='mean', y='rel. err.', title='Scattering Rates')"
]
},
{
"cell_type": "code",
"execution_count": 39,
"metadata": {
"collapsed": false
},
"outputs": [
{
"data": {
"text/plain": [
"<matplotlib.legend.Legend at 0x7f2a2acef290>"
]
},
"execution_count": 39,
"metadata": {},
"output_type": "execute_result"
},
{
"data": {
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAYgAAAEZCAYAAACNebLAAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAIABJREFUeJzt3Xl4VOX5//H3nbDIYtgKhDURFBGsCqIoVA11qdqiaJXi\nVq1rFZeqVVv9fQturVttXQpuVEFFRdEirmzGlU1BUEGUVWSXLaDs3L8/zgGTOCSTZGbOJPm8rmuu\nzJz1fjLJ3PMs5znm7oiIiBSXEXUAIiKSnpQgREQkJiUIERGJSQlCRERiUoIQEZGYlCBERCQmJQiR\nEpjZX83ssajjEImCEoSknJn9wsw+NLN1Zvadmb1vZodW8Jjnm9n7xZY9aWa3VeS47v4Pd7+0IsfY\nEzPbaWYbzKzAzBab2T/NzOLc9xgzW5yMuER2qRF1AFK9mNnewGjgMuBFoBZwFLCloocGEnrVp5ll\nuvuORB6zGAcOcvcFZtYOeA+YBQyJJzwSXF6R4lSDkFTrALi7j/DAFncf5+6f79rAzC4xs1nhN+vP\nzeyQcPlNZja30PI+4fKOwGDgyPAb+RozuwQ4B7gx3H5UuG0LM3vJzFaa2Twzu6rQeQeY2Ytm9rSZ\nrQPOD5c9Ha7PCb/1/97MFoXHuLnQ/nuZ2dDw/F+Y2Q2lfMu38IG7zwc+BA4pdLwLCv0e5prZpeHy\nusAbQMtCNZBsC/wl3HaVmT1vZg3DfWqH5frOzNaa2WQza1reN1GqByUISbWvgB1m9pSZnbjrA2wX\nMzsT+BtwrrtnAacAq8PVc4Ge4fJbgWfMrLm7fwn8EZjo7nu7e2N3fxx4FrjH3bPc/dSw+WY0MB1o\nARwLXGNmxxcK4RRghLs3BIaHy4p/U+8J7AccB/zNzPYPlw8E2gK5wPHAuTH2jSlMckcBXxdavAI4\nOSzvH4B/mdkh7v4DcBKwNCxvlrsvB64O4z8KaAmsBQaFxzofyAJaAY3D39emeGKT6ksJQlLK3TcA\nvwB2Ao8BK81sVKFvsxcRfKhPC7ef7+6Lw+cj3X1F+PxFgg/Tw8tw+sOAn7n7ne6+w90XAk8A/Qpt\nM9HdR4fn2ByrCMBAd9/q7jOBGcDB4bozgTvdvcDdlwIPxhHTNDPbSNC09A5BTYjw/G+GMeLu7wNj\nCD789+Qy4BZ3X+bu24DbgDPMLAPYBjQBOoQ1t+nuvjGO+KQaU4KQlHP3Oe5+obu3BQ4k+Lb773B1\nG2BerP3Cpp3pYRPJWqAz8LMynDoHaBU2Aa0Jj/FXoFmhbeLp+F1R6PkPQP3weUvg2zIeq4u71wf6\nAt2BertWmNlJZjbRzFaHsZ5EyeXNAV7ZVT6CpLMNaA48DbwNPG9m35rZXWaWGUd8Uo0pQUik3P0r\n4CmCRAHBh2r74tuZWVuCGscV7t7I3RsBXxC24RO7Kaf4ssXA/LAJqnF4nAbu3ruEfcpiGdC60Ou2\nceyzqw/iJWASMADAzGoBLwH3AE3D8r5JyeX9BjipWPnqhTWK7e5+u7t3BnoAvYHfl72IUp0oQUhK\nmdn+ZnadmbUKX7cBzgImhps8AfzZzLqG69uH29QjaJb6zswyzOwP/JhUIPhW39rMahZb1q7Q6ynA\nBjO7MexQzjSzzmbWrSxFKGHdCOCvZtYwLF//MhwX4C7gEjNrRjC6qxbwnbvvNLOTgBMKbbsCaGJm\nWYWWPQr8PUymmFlTMzslfJ5nZgeGzU0bCWoWO8sYn1QzShCSahsImlImm9kG4CNgJvBn2P1N+k5g\nuJkVAK8Ajd19NvBPgm/Zywmalz4odNwJBDWK5Wa2Mlw2BOgcNrm87O47gd8QjBRaAKwEHifovI1X\n8W/uhV/fBiwJjz2GYBhvScN3ixwrHMn1LnBD2D9wDfBi2FzUDxhVaNs5wHPA/LB82cAD4TZjzGw9\nwe92Vx9NNkGNZD3B7+kdgmYnkT2yZN4wyMxaA8MI2kB3Ao+5+0NmNgC4hOAfFOBmd38raYGIRMDM\n/gj8zt17RR2LSHkk+0K57cB17v6pmdUHPjGzseG6+939/iSfXyRlwm/x7QiayzoA1xPfSCaRtJTU\nBBGOzV4ePt9oZrMJxmFDyW25IpVRLYJ+gFxgHUET0OCSdhBJZ0ltYipyIrNcIJ+gY/F64AKC9tCP\ngevdfX1KAhERkbikpJM6bF56Cbgm7HwbBLRz90MIahhqahIRSTNJr0GYWQ3gNeBNd38gxvocYLS7\nHxRjnSYjExEpB3evcDN+KmoQ/wVmFU4OYWfeLqcDn/9kr5C7V9nHgAEDIo9B5VPZVL6q90iUpHZS\nm1lPghk1PzOz6QTjvm8GzrZghs6dwEKCOWRERCSNJHsU04dArPledM2DiEia05XUEcrLy4s6hKSq\nyuWrymUDlU8CKRvmWh5m5ukcn4hIOjIzPAGd1LrlqIj8RG5uLosWLYo6DClFTk4OCxcuTNrxVYMQ\nkZ8Iv4FGHYaUYk/vU6JqEOqDEBGRmJQgREQkJiUIERGJSQlCRERiUoIQkUpln332YcKECbtfP//8\n8zRp0oT33nuPjIwMsrKyyMrKokWLFpxyyimMGzeuyP65ubnUrVuXrKws9t57b7Kysrj66qtTXYxK\nQQlCRCqtoUOHctVVV/HGG2+Qk5ODmbF+/XoKCgqYMWMGxx13HKeddhrDhg3bvY+Z8frrr1NQUMCG\nDRsoKCjgwQd1X6dYlCBEpFJ69NFHueGGGxgzZgzdu3ffvXzXsM9mzZpx9dVXM3DgQG688cYi+2oI\nb3yUIESk0hk0aBADBw5kwoQJdOnSpcRtTz/9dFauXMmcOXNSFF3VoSupRaTM7NbE3DHYB5Tvm/y4\ncePo1asXBx54YKnbtmzZEoA1a9bsXtanTx9q1KiBu2Nm3HvvvVx00UXliqUqU4KQtJedncuKFUWn\nfWjePIflyxdGE5CU+4M9UQYPHswdd9zBRRddxJAhQ0rcdsmSJQA0adJk97JRo0bRq1evpMZYFaiJ\nSdJekBy8yKN4wpDqpXnz5owfP57333+fK664osRtX375ZZo3b06HDh12L1MfRHyUIESkUsrOzmb8\n+PG8/fbbXH/99QBF7qi2cuVKHn74YW6//XbuuuuuKEOttNTEJCKVitmP/R9t2rRh/PjxHHPMMSxf\nvhwzo1GjRrg79erVo1u3brz00kscf/zxRY7Ru3dvMjN/vJfZ8ccfz8iRI1NWhspCs7lK2gs+EIr/\nHWi20WTSbK6Vg2ZzFRGRSChBiIhITEoQIiISkxKEiIjEpAQhIiIxKUGIiEhMShAiIhKTEoSIiMSk\nBCEiVcqBBx7Ie++9F3UYVYIShIjEJTs7FzNL2iM7OzeuOIrfchSCO8sdddRRAHz++eccffTRJR5j\n0aJFZGRksHPnznL9LqoLzcUkInH5cVbdZB2/YjNDFJ6jqTS77gORrOlEduzYUWSup8pKNQgRqVIK\n1zCmTp3KYYcdRoMGDWjRogV//vOfATjmmGMAaNiwIVlZWUyePBl354477iA3N5fs7GwuuOACCgoK\ndh932LBh5Obm0rRpU+64444i57n11ls588wzOe+882jYsCFDhw5l6tSp9OjRg0aNGtGqVSuuuuoq\ntm/fvvt4GRkZDB48mA4dOtCgQQP+9re/MX/+fHr27EnDhg3p169fke2joAQhIpXenmoC11xzDX/6\n059Yv3498+bNo2/fvgC7+ygKCgooKCige/fuPPnkkwwbNox3332X+fPns2HDBq688koAZs2aRf/+\n/XnuuedYtmwZ69evZ+nSpUXO9eqrr9K3b1/WrVvHOeecQ40aNfj3v//NmjVrmDhxIhMmTGDQoEFF\n9hkzZgzTp09n0qRJ3HPPPVx22WUMHz6cxYsX89lnn/Hcc88l+ldVJkoQIlLp9OnTh8aNG+9+9O/f\nP+Z2tWrVYu7cuaxevZq6dety+OGHF1lfOLEMHz6c6667jpycHOrWrcs//vEPXnjhBXbu3MnIkSM5\n5ZRTOPLII6lRowa33XbbT8515JFH0rt3bwBq165Nly5dOPzwwzEz2rZty6WXXsq7775bZJ+bbrqJ\nevXqccABB3DggQdywgknkJOTw957781JJ53E9OnTK/qrqhAlCBGpdEaNGsWaNWt2P4p/M99lyJAh\nzJkzh44dO9K9e3def/31PR5z6dKl5OTk7H6dk5PD9u3bWbFiBUuXLqVNmza719WpU6fILUyBIusB\nvv76a3r37k2LFi1o2LAht9xyC999912RbZo1a1bkmM2bNy/yeuPGjSX8FpJPCUJEKp14O5fbt2/P\n8OHDWbVqFTfeeCNnnHEGmzZtitmh3bJlSxYt+vFWtosWLaJGjRo0b96cFi1a8O233+5et2nTJlav\nXl1k/+LHvPzyyznggAOYN28e69at484776x099hQghCRKuvZZ5/d/a29QYMGmBkZGRk0bdqUjIwM\n5s2bt3vbs846i3/9618sXLiQjRs3csstt9CvXz8yMjI444wzGD16NJMmTWLbtm0MHDiw1HNv2LCB\nrKws6taty5dffsngwYOTVcykUYIQkbg0b54DWNIewfFLV9pw1sLr33rrLTp37kxWVhbXXnstL7zw\nArVr16ZOnTrccsst9OzZk8aNGzNlyhQuvPBCzjvvPI4++mjat29P3bp1efDBBwHo1KkTDz30EL/7\n3e9o2bIlWVlZNGvWjNq1a+8xjvvuu49nn32WrKwsLrvsMvr161diOcoyTDdVknrLUTNrDQwDmgM7\ngcfd/UEzawS8AOQAC4G+7r4+xv665ajolqMR0C1HS/b999/TsGFD5s6dW6TfItUq+y1HtwPXuXtn\n4Eigv5l1BP4CjHP3/YEJwF+THIeISIW89tprbNq0ie+//57rr7+egw46KNLkkApJTRDuvtzdPw2f\nbwRmA62BU4Gh4WZDgT7JjENEpKJGjRpFy5Ytad26NfPmzeP555+POqSkS2oTU5ETmeUC+cCBwGJ3\nb1Ro3Rp3bxxjHzUxiZqYIqAmpsoh2U1MKZmLyczqAy8B17j7RjMrXqI9/iUWHi2Ql5dHXl5eMkIU\nEam08vPzyc/PT/hxk16DMLMawGvAm+7+QLhsNpDn7ivMLBt4x90PiLGvahCiGkQEVIOoHCp7JzXA\nf4FZu5JD6FXggvD5+cCoFMQhIiJlkNQmJjPrCZwDfGZm0wm+Bt4M3A2MMLMLgUVA32TGISJlk5OT\nk5bj8qWoZI+iSlkndXmoiUlATUwiZVWZmphERKQSUoIQEZGYlCBERCQmJQgREYlJCUJERGJSghAR\nkZiUIEREJCYlCBERiUkJQkREYlKCEBGRmJQgREQkJiUIERGJSQlCRERiUoIQEZGYlCBERCQmJQgR\nEYlJCUJERGJSghARkZiUIEREJCYlCBERiUkJQkREYlKCEBGRmJQgpErIzs7FzIo8srNzq9w5RVLJ\n3D3qGPbIzDyd45PUMDOg+N+BUfhvI55toohLJApmhrtbRY+jGoSIiMSkBCEiIjEpQYiISExKECIi\nEpMShIiIxKQEISIiMSlBiIhITEoQIiISkxKEiIjEpAQhIiIxKUGIiEhMSU0QZjbEzFaY2cxCywaY\n2bdmNi18nJjMGEREpHySXYN4EvhVjOX3u3vX8PFWkmMQEZFySGqCcPcPgLUxVlV4lkEREUmuqPog\nrjSzT83sCTNrEFEMIiJSgigSxCCgnbsfAiwH7o8gBhERKUWNVJ/Q3VcVevk4MLqk7QcOHLj7eV5e\nHnl5eUmJS0SkssrPzyc/Pz/hx036HeXMLBcY7e4/D19nu/vy8Pm1wGHufvYe9tUd5UR3lBMpo0Td\nUS6pNQgzGw7kAU3M7BtgANDLzA4BdgILgcuSGYOIiJSP7kktaU81CJGy0T2pRUQkqZQgREQkprgS\nhJm9bGa/NjMlFBGRaiLeD/xBwNnA12Z2l5ntn8SYREQkDcSVINx9nLufA3QlGHk0zsw+MrM/mFnN\nZAYoIiLRiLvJyMyaABcAFwPTgQcIEsbYpEQmUmG1MbMij8zMeuVelp2dW64osrNzE3YskVSKa5ir\nmb0C7A88DTzl7ssKrfvY3bslJTgNcxUqNsw10cvKM7RWw2El1RI1zDXeBHGyu79RbFltd99S0QBK\nOa8ShChBiJRRqq+DuCPGsokVPbmIiKSvEqfaMLNsoBVQx8y68ON9HLKAukmOTUREIlTaXEy/IuiY\nbk3Rabk3ADcnKSYREUkD8fZB/NbdR6YgnuLnVR+EqA9CpIxSMpurmZ3r7s8AuWZ2XfH17q6b/YiI\nVFGlNTHVC3/WT3YgIiKSXjTdt6Q9NTGJlE1Kh7ma2T1mlmVmNc1svJmtMrNzK3pyERFJX/FeB3GC\nuxcAvyGYi2lf4IZkBSUiItGLN0Hs6qv4NfCiu69PUjwiIpIm4r0n9Wtm9iWwCbjczJoCm5MXloiI\nRC3uTmozawysd/cdZlYXyHL35UkNTp3UgjqpRcoqJddBFNOR4HqIwvsMq2gAIiKSnuJKEGb2NNAe\n+BTYES52lCBERKqseGsQ3YBOau8REak+4h3F9DmQncxARCD23deqiyjuPFf8nLrTnRQW72R97wCH\nAFOA3TcJcvdTkheaOqmro7J0Nle1TuooOrN/ek51nlcFqe6kHljRE4mISOUSV4Jw93fNLAfYz93H\nhcNcM5MbmkgZZAH7DING86HGZtiyN6wGliyC9TlRRydSKcU7F9MlwEvAo+GiVsD/khWUSLymL5vO\nyc+eDH8EOrwGthM2N4DaBXAQcOlhcGVHOPav0OSrqMMVqVTi7YP4FDgcmOzuXcJln7n7z5ManPog\nqp2y9CXc/cHd3PfRfdyadytXHHkF7Iixn+2A7Onw8+fgoGfguxXw0avw9a/BM0o8vvogpLJKVB9E\nvAlisrt3N7Pp7t4lvFhumrsfVNEASjmvEkQ1E1+CcPhVBgf3OZjXzn6N1lmt49svYxscUAt+0QUy\nt8L7f4XPzwLPjOOcwTIlCKkMUt1J/a6Z3QzUMbPjgSuA0RU9uUi5/OJuyIX8C/JpuFfD+PfbWRO+\nAL74BNqPhbwB0PNeGAvMc4KkICK7xFuDyAAuAk4g+C96G3gi2V/vVYOofkqtCewzHk47Hx5fghdU\ndJirwwGvwLG/hYJj4e1/woqDS9xXNQipDFLaxBSesCmAu6+q6EnjpQRRtWVn57JixaIYa/bwoV67\nAC7/Obz2KMztQ6FLckrer7RlGQZdB0GvAfDZWfDObbClYYx994rrnIlNEEXP2bx5DsuXLyRRlCCq\nppQkCAv+egYAV/LjiKcdwEPufltFT15qcEoQVVqZv/Wf3D/oOxj9eIztEnBRXN3vgtFOHV6Hsctg\n5k6KNjvFd85E1yCS+QGuBFE1peqWo9cCPYHD3L2xuzcGugM9zezaip5cJG5N5kDnETDuruSd44ef\nBcnnhZfhSOCCPGj6RfLOJ5LmSqtBTAeOd/fvii1vCozZNeQ1acGpBlGllakG0fd0WHI4fHjTHrZL\n8LQaZtDtYcgbCNMuhvf+D7bVi+ucqkFI1FJVg6hZPDnA7n6ImqUd3MyGmNkKM5tZaFkjMxtjZnPM\n7G0za1D2sKVaaQ60ngSTr07dOR2Y2h8Gz4SGi6B/J+iQutOLpIPSEsTWcq7b5UngV8WW/QUY5+77\nAxOAv8ZxHKnOehAkh+11Un/ujS1g5HB49YlgDF+/PtDgm9THIRKB0hLEwWZWEOOxASj1Kmp3/wBY\nW2zxqcDQ8PlQoE+Zo5bqI2sx7Ad8clm0ccw/DgYDSw+Fy7pCz3vinyxfpJIq8U/c3TPdPSvGY293\nL7WJaQ+aufuK8PjLgWblPI5UB4f/B2YAm8twQVyy7CDoi3h8MuS+A5cBbT+IOiqRpEmH70DqEZPY\nMrbBwUPhk6gDKWZte3j2DXgXOKMfnHJRMERWpIqJd6qNRFphZs3dfYWZZQMrS9p44MCBu5/n5eWR\nl5eX3Ogkfez3RvBh/N3yqCOJwWAWMG9WMGXHFZ3h3QFpk8xiXYRY3ovsEnksSY78/Hzy8/MTfty4\nr6Qu9wnMcoHRu2Z+NbO7gTXufreZ3QQ0cve/7GFfDXOtwkod1nnWKTD7NPj0wpK3i/k6FcsKvW4+\nA074MzQYx8tXvEyfjn123y41imGu5T9nekwBIhWT8qk2ynVws+FAHtAEWEFwVfb/gBeBNsAioK+7\nr9vD/koQVViJH4j1lwVDS/+1GLbuveft9vg6FctibNPeOOjPB1G/Vn3uPf5eerTpoQQhKZfq2VzL\nxd3P3sOq45J5XqkCOr0EX/WGrfWjjqRs5sG0S6fxzMxn6PdSP7q17BZ8PVoddWAiZZcOndQiP9X5\nRfjizKijKJfMjEzOP+R85lw5hyNaHwEXAr/5Y1ArEqlElCAk/dRfBs0+g3knRB1JhdSpWYcbe94I\nDxPcI/uKA6HX/wWz0opUAkoQkn46jQyal3bUjjqSxNgEjL0XHp0WXIV91X7QHbbuiGcyApHoKEFI\n+ulUeZuXSrQ+B/43FJ4eC/tCx4c78txnz7HTd0YdmUhMSR/mWhEaxVS1xRwdU8/gyoZw3/JCNYhy\njiiKYhRTWW4ElLsFjg8P8zbBmL5yjGIqy42XNIqpekjVbK4iqbUfMP/4qtO8tEdbYKHD4zvhw+fh\n9NZwGlBvRZmPFCQHL/YQqTglCEkvHYCvfh11FClk8MXv4D+zYSNBR/Zhg8B2RB2YiBKEpJHMrdAO\nmHtS1JGk3tb6MBZ4Kh86vwCXdIeWU6OOSqo5JQhJH20/CC4o+74aT/C7qnOQJCZfDWf3hpOA2uuj\njkqqKSUISR/7vQFfRR1EOjCY8Xv4z6xgroP+nYL7catvQVJMCULSx36vw9dRB5FGNjWG0cCLI+Do\n2+Gck6FR1EFJdaIEIemh0XyosxY0G8VPLe4ZXGS3sBdcAn9//++6yE5SQglC0kP7McHUGmpFiW1n\nTfjwRngMPlr8EYc8cgjvLXov6qikitOFchKZIhdgndk3GN464wISedFa6i+U2wvYwk8l6pzh8Q8A\nTgQWAO8A6+M5VjyxleFCv2LH0k2E0oculJOqw3YG93he8MuoI0mALST3orXw+LMd/lMABTcH98Y+\n+UrYe2kKY/vpsWJfzS2VmRKERK/ZZ7C5ERS0iTqSymXr3jDhzmC22O17BRfZnXJx8PsUSQAlCIne\nPhOqSO0hIj8AY+6Dh+bA2n3gvF/B74H9X9UV2VIh6oOQyOxu2z6rN8w8D77oS6Inzotisr7Iz5m5\nFTrXhsO7Q/3lMPUKmHYRbPpZXMcv361QY+8r0VAfhFQNGdsh531Y0CvqSKqOHbVgJvDEJBjxIjSd\nBVfvC6cCLT6JOjqpRFSDkMiYGbSaBL0vhUdm7FpKlfo2ny7nrLsKujaDbm1hQ0uYcAcsODbm8VWD\nqPwSVYOokYhgRMpN/Q+p8UNT+AD4cH5wx75TLoEVP4fXgQ1RByfpSk1MEq1245UgUskzg76e/8yC\nZV2DIbIHvBx1VJKm1MQkkbGaBjfUh/u/hS0Ndi2lyjf3pNM5Wxmc2RY+uRTevxmoQ0Uu9NP/a3pQ\nE5NUfq0JprfenRwk5ZYAQyYGI8nqL4c3d10AV1iFP2ekklITk0RnH9S8lA42tIShE6DNRDgONCGW\n7KIEIdFRgkgfWxrA028Ht3zt9kjU0UiaUIKQSGzcuhGygcU9og5FdtnUBJ4H8gZCmw+jjkbSgBKE\nROL9Re/DUmBb3ahDkcLWAKOehDN/B3XWRB2NREwJQiIxYcEEmB91FBLT1yfDrN/CSVdHHYlETAlC\nIjFh4YTgXgaSnsb/A1pNgY6vRB2JREgJQlJuzaY1fL3666CJSdLTtrowagicdA3UjDoYiYoShKTc\nuwvfpWfbnqCZqNPbN0fBoqPgqKgDkagoQUjKTVgwgWP3OTbqMCQeY++BbkAjdRhVR0oQknLjF4zn\nl/vo+odKYUMrmAzkDYg6EomAEoSk1LINy1i+cTkHNz846lAkXhOB9mOg6RdRRyIppgQhKfXOwnfI\ny80jMyMz6lAkXluBD2+EX/5f1JFIikWWIMxsoZnNMLPpZjYlqjgktcbNH6f+h8po6hXBsNcW06KO\nRFIoyhrETiDP3bu4++ERxiEp4u6MXzCeY9spQVQ62+vAR9dDz7ujjkRSKMoEYRGfX1Js3tp5bN+5\nnf2b7B91KFIe0y4J7gDYeG7UkUiKRHk/CAfGmtkO4DF3fzzCWCSJduzYwahRo3hz5Zu0t/aMGDEi\n6pCkPLbWh48vhx73wWtRByOpEGWC6Onuy8ysKUGimO3uHxTfaODAgbuf5+XlkZeXl7oIJSHGjh3L\nOef0Z/tpNagxvyWXzHqZrVtnRx2WlMfkq+Cq/SEf2Bh1MLJLfn4++fn5CT9uWtxy1MwGABvc/f5i\ny3XL0Srg9ddf5+xz/kPBZVPgkU+hoDVwD3ATaXsrzmp7zji2OflK2PwfmKBbjqarRN1yNJI+ADOr\na2b1w+f1gBOAz6OIRVJjZ9MC+KFJmBykUpt8FXQFamyOOhJJsqg6iZsDH5jZdGASMNrdx0QUi6TA\n9rarYYFGL1UJq/eHZUDnF6KORJIskj4Id18AHBLFuSUa29ushinHRR2GJMpk4NgHYMbvCZqgpCrS\nMFNJum07t7G95VpYmBd1KJIo84BaG6HNR1FHIkmkBCFJN+f7OWSuqwubGkcdiiSKE/RFdH8w6kgk\niZQgJOk+Xv8xNRY2jToMSbQZ50P7sZC1OOpIJEmUICTpPln/iRJEVbQlC2aeC90ejToSSRIlCEmq\nJQVLWLV1FZnLGkYdiiTDlCuh6+Ma8lpFKUFIUr019y26ZHXBXH9qVdLqDrC8C3TW9ClVkf5rJane\nnPsm3Rp0izoMSaYpV8LhD0cdhSSBEoQkzbYd2xi/YDxdG3SNOhRJpq9PgjqroVXUgUiiKUFI0kz8\ndiLtGrWjUc1GUYciyeSZwQ2FukcdiCSaEoQkzeg5o/n1fr+OOgxJhekXwn6wYuOKqCORBFKCkKRw\nd17+8mVOP+D0qEORVNjcCGbBY588FnUkkkBKEJIUM1fMxN05uPnBUYciqTIFHvnkEbbt2BZ1JJIg\nShCSFK9Nm+nGAAAKi0lEQVR8+QqndTwNM03kVm2sgH0b78srX74SdSSSIEoQkhQvz1bzUnV01eFX\n8dCUh6IOQxJECUISbu6auaz6YRVHtjky6lAkxU7d/1QWrF3Ap8s/jToUSQAlCEm4EV+M4LSOp5Fh\n+vOqbmpm1uTybpfz8BRdOFcV6D9YEsrdeXrm05x70LlRhyIRueTQSxg5eySrf1gddShSQUoQklDT\nlk1j646tHNlazUvVVbN6zejdoTdPTHsi6lCkgpQgJKGemfkM5/78XI1equauP/J6Hpj8AJu3a5bX\nykwJQhJm+87tPP/F82peEg7OPphDWx7Kf6f/N+pQpAKUICRhRn05in0b78t+TfaLOhRJA7ccdQv3\nfHiPLpyrxJQgJGEGfTyI/of1jzoMSRNHtD6C9o3b8+xnz0YdipSTEoQkxOxVs5m1apYujpMibjnq\nFv7xwT/YvnN71KFIOShBSEIM/ngwF3e5mFqZtaIORdJIr9xetNy7JU9OfzLqUKQclCCkwr774Tue\nmfkMf+z2x6hDkTRjZtxz3D0MfHcg32/9PupwpIyUIKTC/j3p35zZ6UxaZemWYvJTh7U6jKNzjub+\nifdHHYqUkRKEVMi6zet45ONHuOkXN0UdiqSxO395Jw9MfoBv1n8TdShSBkoQUiF/f//v9OnYh3aN\n2kUdiqSxdo3acU33a+j/Rn/cPepwJE5KEFJu89fOZ8j0Idze6/aoQ5FK4KZf3MT8tfMZOXtk1KFI\nnJQgpFzcnT+99SeuPeJaWuzdIupwpBKolVmLx37zGFe/eTXLNy6POhyJgxKElMuznz3LgnULuKHH\nDVGHIpVIz7Y9ubjrxZz78rns2Lkj6nCkFEoQUmaL1i3i+jHX89SpT1G7Ru2ow5FKZsAxA9i+czsD\n8gdEHYqUQglCyuSHbT9w2guncWOPGzm05aFRhyOVUGZGJiPOHMHznz/PY588FnU4UoIaUQcglcfW\nHVs5a+RZdG7WmeuOvC7qcKQSa1avGW+d+xZHP3k0dWvW1QzAaUoJQuKyadsmznn5HNydIacM0f0e\npML2bbwv434/jhOfOZGV36/k2iOu1d9VmomsicnMTjSzL83sKzPTVVZp7Jv133D0U0dTK7MWL575\nouZbkoTp1LQTH1z4AUNnDKXvS31Zu2lt1CFJIZEkCDPLAB4GfgV0Bs4ys45RxBKl/Pz8qEMo0Zbt\nW3hg0gMc+tih9O3Ul+d++1yZOqXTvXwVkx91AEmWn7IztW3QlskXT6ZF/RZ0GtSJxz95POmzv1bt\nv83EiaoGcTjwtbsvcvdtwPPAqRHFEpl0/SNdtG4Rd31wF+0fbM9b897ivQve44aeN5S5+p+u5UuM\n/KgDSLL8lJ5trxp78eBJD/LaWa8x/PPhtHugHXe+dyfz185Pyvmq9t9m4kTVB9EKWFzo9bcESUNS\naOuOraz8fiXz1sxj3tp5fLz0Yz5c/CFLCpbw2wN+y//6/Y9uLbtFHaZUI4e2PJR3zn+H6cum88jH\nj9BjSA8a1WlEj9Y96NqiK52adqJ1VmtaZbWibs26UYdb5amTOgJ3vHcHE7+dyFczv2Lys5Nxdxzf\n/RP4ybI9/Szrtpu3b2bd5nWs37yerTu20rReU9o3ak+7Ru3okt2F8w8+n64tulIzs2bCyluzZk02\nb55KVlbv3cu2bJnLli0JO4VUMV1adOHR3o8y+DeDmb5sOlOXTmXasmmMmDWCJQVL+LbgWzIzMqlf\nq/7ux1419iLTMsnMyCTDMn7yvHAN+KvPvmLq8Kl7PL9Rem35Z3V/xlN9nkpEcdOWRTFxlpkdAQx0\n9xPD138B3N3vLradZvUSESkHd6/wkLCoEkQmMAc4FlgGTAHOcvfZKQ9GRERiiqSJyd13mNmVwBiC\njvIhSg4iIuklkhqEiIikv8jnYjKzRmY2xszmmNnbZtZgD9sNMbMVZjazPPtHoQxli3nRoJkNMLNv\nzWxa+DgxddHvWTwXOZrZg2b2tZl9amaHlGXfqJWjfF0KLV9oZjPMbLqZTUld1PErrXxmtr+ZfWRm\nm83surLsG7UKlq0qvHdnh2WYYWYfmNlB8e4bk7tH+gDuBm4Mn98E3LWH7X4BHALMLM/+6Vo2giQ9\nF8gBagKfAh3DdQOA66IuR7zxFtrmJOD18Hl3YFK8+0b9qEj5wtfzgUZRl6OC5fsZcChwe+G/v3R/\n/ypStir03h0BNAifn1jR/73IaxAEF8gNDZ8PBfrE2sjdPwBiXYcf1/4RiSe20i4aTLfJaeK5yPFU\nYBiAu08GGphZ8zj3jVpFygfB+5UO/1d7Umr53P07d/8EKH45c7q/fxUpG1SN926Su68PX04iuOYs\nrn1jSYdfRjN3XwHg7suBZineP5niiS3WRYOtCr2+MmzGeCJNms9Ki7ekbeLZN2rlKd+SQts4MNbM\npprZJUmLsvwq8h6k+/tX0fiq2nt3MfBmOfcFUjSKyczGAs0LLyJ4M/5fjM0r2mue0l73JJdtEHCb\nu7uZ3QHcD1xUrkCjlW61oGTq6e7LzKwpwYfN7LD2K+mvyrx3ZtYL+ANB03y5pSRBuPvxe1oXdjw3\nd/cVZpYNrCzj4Su6f4UkoGxLgLaFXrcOl+HuqwotfxwYnYCQK2qP8Rbbpk2MbWrFsW/UKlI+3H1Z\n+HOVmb1CULVPpw+ZeMqXjH1ToULxVZX3LuyYfgw40d3XlmXf4tKhielV4ILw+fnAqBK2NX76bbQs\n+6daPLFNBfY1sxwzqwX0C/cjTCq7nA58nrxQ47bHeAt5Ffg97L5qfl3Y1BbPvlErd/nMrK6Z1Q+X\n1wNOID3es8LK+h4U/n9L9/ev3GWrKu+dmbUFRgLnufu8suwbUxr0zDcGxhFcWT0GaBgubwG8Vmi7\n4cBSYAvwDfCHkvZPh0cZynZiuM3XwF8KLR8GzCQYcfA/oHnUZdpTvMBlwKWFtnmYYNTEDKBraWVN\np0d5ywfsE75X04HPKmv5CJpMFwPrgDXh/1v9yvD+lbdsVei9exxYDUwLyzKlpH1Le+hCORERiSkd\nmphERCQNKUGIiEhMShAiIhKTEoSIiMSkBCEiIjEpQYiISExKECKAme00s2GFXmea2SozS6cLwURS\nSglCJPA9cKCZ1Q5fH0/Ryc1Eqh0lCJEfvQH8Onx+FvDcrhXhVAxDzGySmX1iZr3D5Tlm9p6ZfRw+\njgiXH2Nm75jZi2Y228yeTnlpRCpICUIk4ARz5J8V1iIOAiYXWn8LMN7djwB+CdxnZnWAFcBx7t6N\nYH6bhwrtcwhwNdAJaG9mPZJfDJHESclsriKVgbt/bma5BLWH1yk6Ud0JQG8zuyF8vWtm2mXAwxbc\nVnUHsF+hfaZ4OEOomX0K5AIfJbEIIgmlBCFS1KvAvUAewe0pdzHgt+7+deGNzWwAsNzdDzKzTGBT\nodVbCj3fgf7fpJJRE5NIYFdt4b/Are7+RbH1bxM0FwUbBzUGgAYEtQgIpgDPTGaQIqmkBCEScAB3\nX+LuD8dYfztQ08xmmtlnwG3h8kHABWY2HehAMBpqj8cXqUw03beIiMSkGoSIiMSkBCEiIjEpQYiI\nSExKECIiEpMShIiIxKQEISIiMSlBiIhITEoQIiIS0/8HBiVkAxzrXDkAAAAASUVORK5CYII=\n",
"text/plain": [
"<matplotlib.figure.Figure at 0x7f2a2ae35f50>"
]
},
"metadata": {},
"output_type": "display_data"
}
],
"source": [
"# Plot a histogram and kernel density estimate for the scattering rates\n",
"scatter['mean'].plot(kind='hist', bins=25)\n",
"scatter['mean'].plot(kind='kde')\n",
"pylab.title('Scattering Rates')\n",
"pylab.xlabel('Mean')\n",
"pylab.legend(['KDE', 'Histogram'])"
]
}
],
"metadata": {
"kernelspec": {
"display_name": "Python 2",
"language": "python",
"name": "python2"
},
"language_info": {
"codemirror_mode": {
"name": "ipython",
"version": 2
},
"file_extension": ".py",
"mimetype": "text/x-python",
"name": "python",
"nbconvert_exporter": "python",
"pygments_lexer": "ipython2",
"version": "2.7.11"
}
},
"nbformat": 4,
"nbformat_minor": 0
}
Reference in a new issue View git blame Copy permalink