diff --git a/notebooks/strexp_fitting.ipynb b/notebooks/strexp_fitting.ipynb
index 459c2cf..4e601dd 100644
--- a/notebooks/strexp_fitting.ipynb
+++ b/notebooks/strexp_fitting.ipynb
@@ -48,18 +48,9 @@
},
{
"cell_type": "code",
- "execution_count": 1,
+ "execution_count": 2,
"metadata": {},
- "outputs": [
- {
- "name": "stderr",
- "output_type": "stream",
- "text": [
- "/sw/miniconda3/lib/python3.6/site-packages/h5py/__init__.py:36: FutureWarning: Conversion of the second argument of issubdtype from `float` to `np.floating` is deprecated. In future, it will be treated as `np.float64 == np.dtype(float).type`.\n",
- " from ._conv import register_converters as _register_converters\n"
- ]
- }
- ],
+ "outputs": [],
"source": [
"from __future__ import (absolute_import, division, print_function)\n",
"\n",
@@ -71,7 +62,7 @@
"from lmfit.models import LinearModel, LorentzianModel, ConstantModel, LinearModel\n",
"\n",
"import qef\n",
- "from qef.io.loaders import load_nexus\n",
+ "from qef.io.loaders import load_dave\n",
"from qef.models.deltadirac import DeltaDiracModel\n",
"from qef.models.tabulatedmodel import TabulatedModel\n",
"from qef.models.resolution import TabulatedResolutionModel\n",
@@ -88,7 +79,7 @@
},
{
"cell_type": "code",
- "execution_count": 2,
+ "execution_count": 3,
"metadata": {},
"outputs": [],
"source": [
@@ -115,23 +106,14 @@
},
{
"cell_type": "code",
- "execution_count": 5,
+ "execution_count": 4,
"metadata": {},
"outputs": [
{
"name": "stdout",
"output_type": "stream",
"text": [
- "First, rewinding head to replay your work on top of it...\n",
- "Fast-forwarded master to ad9365c5454e952489223a29cacd15dc003afe2a.\n"
- ]
- },
- {
- "name": "stderr",
- "output_type": "stream",
- "text": [
- "From https://github.com/jmborr/qef_data\n",
- " ee29426..ad9365c master -> origin/master\n"
+ "Current branch master is up to date.\n"
]
}
],
@@ -148,7 +130,7 @@
},
{
"cell_type": "code",
- "execution_count": 24,
+ "execution_count": 5,
"metadata": {},
"outputs": [
{
@@ -176,9 +158,812 @@
},
{
"cell_type": "code",
- "execution_count": 26,
+ "execution_count": 27,
"metadata": {},
- "outputs": [],
+ "outputs": [
+ {
+ "data": {
+ "application/javascript": [
+ "/* Put everything inside the global mpl namespace */\n",
+ "window.mpl = {};\n",
+ "\n",
+ "\n",
+ "mpl.get_websocket_type = function() {\n",
+ " if (typeof(WebSocket) !== 'undefined') {\n",
+ " return WebSocket;\n",
+ " } else if (typeof(MozWebSocket) !== 'undefined') {\n",
+ " return MozWebSocket;\n",
+ " } else {\n",
+ " alert('Your browser does not have WebSocket support.' +\n",
+ " 'Please try Chrome, Safari or Firefox ≥ 6. ' +\n",
+ " 'Firefox 4 and 5 are also supported but you ' +\n",
+ " 'have to enable WebSockets in about:config.');\n",
+ " };\n",
+ "}\n",
+ "\n",
+ "mpl.figure = function(figure_id, websocket, ondownload, parent_element) {\n",
+ " this.id = figure_id;\n",
+ "\n",
+ " this.ws = websocket;\n",
+ "\n",
+ " this.supports_binary = (this.ws.binaryType != undefined);\n",
+ "\n",
+ " if (!this.supports_binary) {\n",
+ " var warnings = document.getElementById(\"mpl-warnings\");\n",
+ " if (warnings) {\n",
+ " warnings.style.display = 'block';\n",
+ " warnings.textContent = (\n",
+ " \"This browser does not support binary websocket messages. \" +\n",
+ " \"Performance may be slow.\");\n",
+ " }\n",
+ " }\n",
+ "\n",
+ " this.imageObj = new Image();\n",
+ "\n",
+ " this.context = undefined;\n",
+ " this.message = undefined;\n",
+ " this.canvas = undefined;\n",
+ " this.rubberband_canvas = undefined;\n",
+ " this.rubberband_context = undefined;\n",
+ " this.format_dropdown = undefined;\n",
+ "\n",
+ " this.image_mode = 'full';\n",
+ "\n",
+ " this.root = $('
');\n",
+ " this._root_extra_style(this.root)\n",
+ " this.root.attr('style', 'display: inline-block');\n",
+ "\n",
+ " $(parent_element).append(this.root);\n",
+ "\n",
+ " this._init_header(this);\n",
+ " this._init_canvas(this);\n",
+ " this._init_toolbar(this);\n",
+ "\n",
+ " var fig = this;\n",
+ "\n",
+ " this.waiting = false;\n",
+ "\n",
+ " this.ws.onopen = function () {\n",
+ " fig.send_message(\"supports_binary\", {value: fig.supports_binary});\n",
+ " fig.send_message(\"send_image_mode\", {});\n",
+ " if (mpl.ratio != 1) {\n",
+ " fig.send_message(\"set_dpi_ratio\", {'dpi_ratio': mpl.ratio});\n",
+ " }\n",
+ " fig.send_message(\"refresh\", {});\n",
+ " }\n",
+ "\n",
+ " this.imageObj.onload = function() {\n",
+ " if (fig.image_mode == 'full') {\n",
+ " // Full images could contain transparency (where diff images\n",
+ " // almost always do), so we need to clear the canvas so that\n",
+ " // there is no ghosting.\n",
+ " fig.context.clearRect(0, 0, fig.canvas.width, fig.canvas.height);\n",
+ " }\n",
+ " fig.context.drawImage(fig.imageObj, 0, 0);\n",
+ " };\n",
+ "\n",
+ " this.imageObj.onunload = function() {\n",
+ " fig.ws.close();\n",
+ " }\n",
+ "\n",
+ " this.ws.onmessage = this._make_on_message_function(this);\n",
+ "\n",
+ " this.ondownload = ondownload;\n",
+ "}\n",
+ "\n",
+ "mpl.figure.prototype._init_header = function() {\n",
+ " var titlebar = $(\n",
+ " '');\n",
+ " var titletext = $(\n",
+ " '');\n",
+ " titlebar.append(titletext)\n",
+ " this.root.append(titlebar);\n",
+ " this.header = titletext[0];\n",
+ "}\n",
+ "\n",
+ "\n",
+ "\n",
+ "mpl.figure.prototype._canvas_extra_style = function(canvas_div) {\n",
+ "\n",
+ "}\n",
+ "\n",
+ "\n",
+ "mpl.figure.prototype._root_extra_style = function(canvas_div) {\n",
+ "\n",
+ "}\n",
+ "\n",
+ "mpl.figure.prototype._init_canvas = function() {\n",
+ " var fig = this;\n",
+ "\n",
+ " var canvas_div = $('');\n",
+ "\n",
+ " canvas_div.attr('style', 'position: relative; clear: both; outline: 0');\n",
+ "\n",
+ " function canvas_keyboard_event(event) {\n",
+ " return fig.key_event(event, event['data']);\n",
+ " }\n",
+ "\n",
+ " canvas_div.keydown('key_press', canvas_keyboard_event);\n",
+ " canvas_div.keyup('key_release', canvas_keyboard_event);\n",
+ " this.canvas_div = canvas_div\n",
+ " this._canvas_extra_style(canvas_div)\n",
+ " this.root.append(canvas_div);\n",
+ "\n",
+ " var canvas = $('');\n",
+ " canvas.addClass('mpl-canvas');\n",
+ " canvas.attr('style', \"left: 0; top: 0; z-index: 0; outline: 0\")\n",
+ "\n",
+ " this.canvas = canvas[0];\n",
+ " this.context = canvas[0].getContext(\"2d\");\n",
+ "\n",
+ " var backingStore = this.context.backingStorePixelRatio ||\n",
+ "\tthis.context.webkitBackingStorePixelRatio ||\n",
+ "\tthis.context.mozBackingStorePixelRatio ||\n",
+ "\tthis.context.msBackingStorePixelRatio ||\n",
+ "\tthis.context.oBackingStorePixelRatio ||\n",
+ "\tthis.context.backingStorePixelRatio || 1;\n",
+ "\n",
+ " mpl.ratio = (window.devicePixelRatio || 1) / backingStore;\n",
+ "\n",
+ " var rubberband = $('');\n",
+ " rubberband.attr('style', \"position: absolute; left: 0; top: 0; z-index: 1;\")\n",
+ "\n",
+ " var pass_mouse_events = true;\n",
+ "\n",
+ " canvas_div.resizable({\n",
+ " start: function(event, ui) {\n",
+ " pass_mouse_events = false;\n",
+ " },\n",
+ " resize: function(event, ui) {\n",
+ " fig.request_resize(ui.size.width, ui.size.height);\n",
+ " },\n",
+ " stop: function(event, ui) {\n",
+ " pass_mouse_events = true;\n",
+ " fig.request_resize(ui.size.width, ui.size.height);\n",
+ " },\n",
+ " });\n",
+ "\n",
+ " function mouse_event_fn(event) {\n",
+ " if (pass_mouse_events)\n",
+ " return fig.mouse_event(event, event['data']);\n",
+ " }\n",
+ "\n",
+ " rubberband.mousedown('button_press', mouse_event_fn);\n",
+ " rubberband.mouseup('button_release', mouse_event_fn);\n",
+ " // Throttle sequential mouse events to 1 every 20ms.\n",
+ " rubberband.mousemove('motion_notify', mouse_event_fn);\n",
+ "\n",
+ " rubberband.mouseenter('figure_enter', mouse_event_fn);\n",
+ " rubberband.mouseleave('figure_leave', mouse_event_fn);\n",
+ "\n",
+ " canvas_div.on(\"wheel\", function (event) {\n",
+ " event = event.originalEvent;\n",
+ " event['data'] = 'scroll'\n",
+ " if (event.deltaY < 0) {\n",
+ " event.step = 1;\n",
+ " } else {\n",
+ " event.step = -1;\n",
+ " }\n",
+ " mouse_event_fn(event);\n",
+ " });\n",
+ "\n",
+ " canvas_div.append(canvas);\n",
+ " canvas_div.append(rubberband);\n",
+ "\n",
+ " this.rubberband = rubberband;\n",
+ " this.rubberband_canvas = rubberband[0];\n",
+ " this.rubberband_context = rubberband[0].getContext(\"2d\");\n",
+ " this.rubberband_context.strokeStyle = \"#000000\";\n",
+ "\n",
+ " this._resize_canvas = function(width, height) {\n",
+ " // Keep the size of the canvas, canvas container, and rubber band\n",
+ " // canvas in synch.\n",
+ " canvas_div.css('width', width)\n",
+ " canvas_div.css('height', height)\n",
+ "\n",
+ " canvas.attr('width', width * mpl.ratio);\n",
+ " canvas.attr('height', height * mpl.ratio);\n",
+ " canvas.attr('style', 'width: ' + width + 'px; height: ' + height + 'px;');\n",
+ "\n",
+ " rubberband.attr('width', width);\n",
+ " rubberband.attr('height', height);\n",
+ " }\n",
+ "\n",
+ " // Set the figure to an initial 600x600px, this will subsequently be updated\n",
+ " // upon first draw.\n",
+ " this._resize_canvas(600, 600);\n",
+ "\n",
+ " // Disable right mouse context menu.\n",
+ " $(this.rubberband_canvas).bind(\"contextmenu\",function(e){\n",
+ " return false;\n",
+ " });\n",
+ "\n",
+ " function set_focus () {\n",
+ " canvas.focus();\n",
+ " canvas_div.focus();\n",
+ " }\n",
+ "\n",
+ " window.setTimeout(set_focus, 100);\n",
+ "}\n",
+ "\n",
+ "mpl.figure.prototype._init_toolbar = function() {\n",
+ " var fig = this;\n",
+ "\n",
+ " var nav_element = $('')\n",
+ " nav_element.attr('style', 'width: 100%');\n",
+ " this.root.append(nav_element);\n",
+ "\n",
+ " // Define a callback function for later on.\n",
+ " function toolbar_event(event) {\n",
+ " return fig.toolbar_button_onclick(event['data']);\n",
+ " }\n",
+ " function toolbar_mouse_event(event) {\n",
+ " return fig.toolbar_button_onmouseover(event['data']);\n",
+ " }\n",
+ "\n",
+ " for(var toolbar_ind in mpl.toolbar_items) {\n",
+ " var name = mpl.toolbar_items[toolbar_ind][0];\n",
+ " var tooltip = mpl.toolbar_items[toolbar_ind][1];\n",
+ " var image = mpl.toolbar_items[toolbar_ind][2];\n",
+ " var method_name = mpl.toolbar_items[toolbar_ind][3];\n",
+ "\n",
+ " if (!name) {\n",
+ " // put a spacer in here.\n",
+ " continue;\n",
+ " }\n",
+ " var button = $('');\n",
+ " button.addClass('ui-button ui-widget ui-state-default ui-corner-all ' +\n",
+ " 'ui-button-icon-only');\n",
+ " button.attr('role', 'button');\n",
+ " button.attr('aria-disabled', 'false');\n",
+ " button.click(method_name, toolbar_event);\n",
+ " button.mouseover(tooltip, toolbar_mouse_event);\n",
+ "\n",
+ " var icon_img = $('');\n",
+ " icon_img.addClass('ui-button-icon-primary ui-icon');\n",
+ " icon_img.addClass(image);\n",
+ " icon_img.addClass('ui-corner-all');\n",
+ "\n",
+ " var tooltip_span = $('');\n",
+ " tooltip_span.addClass('ui-button-text');\n",
+ " tooltip_span.html(tooltip);\n",
+ "\n",
+ " button.append(icon_img);\n",
+ " button.append(tooltip_span);\n",
+ "\n",
+ " nav_element.append(button);\n",
+ " }\n",
+ "\n",
+ " var fmt_picker_span = $('');\n",
+ "\n",
+ " var fmt_picker = $('');\n",
+ " fmt_picker.addClass('mpl-toolbar-option ui-widget ui-widget-content');\n",
+ " fmt_picker_span.append(fmt_picker);\n",
+ " nav_element.append(fmt_picker_span);\n",
+ " this.format_dropdown = fmt_picker[0];\n",
+ "\n",
+ " for (var ind in mpl.extensions) {\n",
+ " var fmt = mpl.extensions[ind];\n",
+ " var option = $(\n",
+ " '', {selected: fmt === mpl.default_extension}).html(fmt);\n",
+ " fmt_picker.append(option)\n",
+ " }\n",
+ "\n",
+ " // Add hover states to the ui-buttons\n",
+ " $( \".ui-button\" ).hover(\n",
+ " function() { $(this).addClass(\"ui-state-hover\");},\n",
+ " function() { $(this).removeClass(\"ui-state-hover\");}\n",
+ " );\n",
+ "\n",
+ " var status_bar = $('');\n",
+ " nav_element.append(status_bar);\n",
+ " this.message = status_bar[0];\n",
+ "}\n",
+ "\n",
+ "mpl.figure.prototype.request_resize = function(x_pixels, y_pixels) {\n",
+ " // Request matplotlib to resize the figure. Matplotlib will then trigger a resize in the client,\n",
+ " // which will in turn request a refresh of the image.\n",
+ " this.send_message('resize', {'width': x_pixels, 'height': y_pixels});\n",
+ "}\n",
+ "\n",
+ "mpl.figure.prototype.send_message = function(type, properties) {\n",
+ " properties['type'] = type;\n",
+ " properties['figure_id'] = this.id;\n",
+ " this.ws.send(JSON.stringify(properties));\n",
+ "}\n",
+ "\n",
+ "mpl.figure.prototype.send_draw_message = function() {\n",
+ " if (!this.waiting) {\n",
+ " this.waiting = true;\n",
+ " this.ws.send(JSON.stringify({type: \"draw\", figure_id: this.id}));\n",
+ " }\n",
+ "}\n",
+ "\n",
+ "\n",
+ "mpl.figure.prototype.handle_save = function(fig, msg) {\n",
+ " var format_dropdown = fig.format_dropdown;\n",
+ " var format = format_dropdown.options[format_dropdown.selectedIndex].value;\n",
+ " fig.ondownload(fig, format);\n",
+ "}\n",
+ "\n",
+ "\n",
+ "mpl.figure.prototype.handle_resize = function(fig, msg) {\n",
+ " var size = msg['size'];\n",
+ " if (size[0] != fig.canvas.width || size[1] != fig.canvas.height) {\n",
+ " fig._resize_canvas(size[0], size[1]);\n",
+ " fig.send_message(\"refresh\", {});\n",
+ " };\n",
+ "}\n",
+ "\n",
+ "mpl.figure.prototype.handle_rubberband = function(fig, msg) {\n",
+ " var x0 = msg['x0'] / mpl.ratio;\n",
+ " var y0 = (fig.canvas.height - msg['y0']) / mpl.ratio;\n",
+ " var x1 = msg['x1'] / mpl.ratio;\n",
+ " var y1 = (fig.canvas.height - msg['y1']) / mpl.ratio;\n",
+ " x0 = Math.floor(x0) + 0.5;\n",
+ " y0 = Math.floor(y0) + 0.5;\n",
+ " x1 = Math.floor(x1) + 0.5;\n",
+ " y1 = Math.floor(y1) + 0.5;\n",
+ " var min_x = Math.min(x0, x1);\n",
+ " var min_y = Math.min(y0, y1);\n",
+ " var width = Math.abs(x1 - x0);\n",
+ " var height = Math.abs(y1 - y0);\n",
+ "\n",
+ " fig.rubberband_context.clearRect(\n",
+ " 0, 0, fig.canvas.width, fig.canvas.height);\n",
+ "\n",
+ " fig.rubberband_context.strokeRect(min_x, min_y, width, height);\n",
+ "}\n",
+ "\n",
+ "mpl.figure.prototype.handle_figure_label = function(fig, msg) {\n",
+ " // Updates the figure title.\n",
+ " fig.header.textContent = msg['label'];\n",
+ "}\n",
+ "\n",
+ "mpl.figure.prototype.handle_cursor = function(fig, msg) {\n",
+ " var cursor = msg['cursor'];\n",
+ " switch(cursor)\n",
+ " {\n",
+ " case 0:\n",
+ " cursor = 'pointer';\n",
+ " break;\n",
+ " case 1:\n",
+ " cursor = 'default';\n",
+ " break;\n",
+ " case 2:\n",
+ " cursor = 'crosshair';\n",
+ " break;\n",
+ " case 3:\n",
+ " cursor = 'move';\n",
+ " break;\n",
+ " }\n",
+ " fig.rubberband_canvas.style.cursor = cursor;\n",
+ "}\n",
+ "\n",
+ "mpl.figure.prototype.handle_message = function(fig, msg) {\n",
+ " fig.message.textContent = msg['message'];\n",
+ "}\n",
+ "\n",
+ "mpl.figure.prototype.handle_draw = function(fig, msg) {\n",
+ " // Request the server to send over a new figure.\n",
+ " fig.send_draw_message();\n",
+ "}\n",
+ "\n",
+ "mpl.figure.prototype.handle_image_mode = function(fig, msg) {\n",
+ " fig.image_mode = msg['mode'];\n",
+ "}\n",
+ "\n",
+ "mpl.figure.prototype.updated_canvas_event = function() {\n",
+ " // Called whenever the canvas gets updated.\n",
+ " this.send_message(\"ack\", {});\n",
+ "}\n",
+ "\n",
+ "// A function to construct a web socket function for onmessage handling.\n",
+ "// Called in the figure constructor.\n",
+ "mpl.figure.prototype._make_on_message_function = function(fig) {\n",
+ " return function socket_on_message(evt) {\n",
+ " if (evt.data instanceof Blob) {\n",
+ " /* FIXME: We get \"Resource interpreted as Image but\n",
+ " * transferred with MIME type text/plain:\" errors on\n",
+ " * Chrome. But how to set the MIME type? It doesn't seem\n",
+ " * to be part of the websocket stream */\n",
+ " evt.data.type = \"image/png\";\n",
+ "\n",
+ " /* Free the memory for the previous frames */\n",
+ " if (fig.imageObj.src) {\n",
+ " (window.URL || window.webkitURL).revokeObjectURL(\n",
+ " fig.imageObj.src);\n",
+ " }\n",
+ "\n",
+ " fig.imageObj.src = (window.URL || window.webkitURL).createObjectURL(\n",
+ " evt.data);\n",
+ " fig.updated_canvas_event();\n",
+ " fig.waiting = false;\n",
+ " return;\n",
+ " }\n",
+ " else if (typeof evt.data === 'string' && evt.data.slice(0, 21) == \"data:image/png;base64\") {\n",
+ " fig.imageObj.src = evt.data;\n",
+ " fig.updated_canvas_event();\n",
+ " fig.waiting = false;\n",
+ " return;\n",
+ " }\n",
+ "\n",
+ " var msg = JSON.parse(evt.data);\n",
+ " var msg_type = msg['type'];\n",
+ "\n",
+ " // Call the \"handle_{type}\" callback, which takes\n",
+ " // the figure and JSON message as its only arguments.\n",
+ " try {\n",
+ " var callback = fig[\"handle_\" + msg_type];\n",
+ " } catch (e) {\n",
+ " console.log(\"No handler for the '\" + msg_type + \"' message type: \", msg);\n",
+ " return;\n",
+ " }\n",
+ "\n",
+ " if (callback) {\n",
+ " try {\n",
+ " // console.log(\"Handling '\" + msg_type + \"' message: \", msg);\n",
+ " callback(fig, msg);\n",
+ " } catch (e) {\n",
+ " console.log(\"Exception inside the 'handler_\" + msg_type + \"' callback:\", e, e.stack, msg);\n",
+ " }\n",
+ " }\n",
+ " };\n",
+ "}\n",
+ "\n",
+ "// from http://stackoverflow.com/questions/1114465/getting-mouse-location-in-canvas\n",
+ "mpl.findpos = function(e) {\n",
+ " //this section is from http://www.quirksmode.org/js/events_properties.html\n",
+ " var targ;\n",
+ " if (!e)\n",
+ " e = window.event;\n",
+ " if (e.target)\n",
+ " targ = e.target;\n",
+ " else if (e.srcElement)\n",
+ " targ = e.srcElement;\n",
+ " if (targ.nodeType == 3) // defeat Safari bug\n",
+ " targ = targ.parentNode;\n",
+ "\n",
+ " // jQuery normalizes the pageX and pageY\n",
+ " // pageX,Y are the mouse positions relative to the document\n",
+ " // offset() returns the position of the element relative to the document\n",
+ " var x = e.pageX - $(targ).offset().left;\n",
+ " var y = e.pageY - $(targ).offset().top;\n",
+ "\n",
+ " return {\"x\": x, \"y\": y};\n",
+ "};\n",
+ "\n",
+ "/*\n",
+ " * return a copy of an object with only non-object keys\n",
+ " * we need this to avoid circular references\n",
+ " * http://stackoverflow.com/a/24161582/3208463\n",
+ " */\n",
+ "function simpleKeys (original) {\n",
+ " return Object.keys(original).reduce(function (obj, key) {\n",
+ " if (typeof original[key] !== 'object')\n",
+ " obj[key] = original[key]\n",
+ " return obj;\n",
+ " }, {});\n",
+ "}\n",
+ "\n",
+ "mpl.figure.prototype.mouse_event = function(event, name) {\n",
+ " var canvas_pos = mpl.findpos(event)\n",
+ "\n",
+ " if (name === 'button_press')\n",
+ " {\n",
+ " this.canvas.focus();\n",
+ " this.canvas_div.focus();\n",
+ " }\n",
+ "\n",
+ " var x = canvas_pos.x * mpl.ratio;\n",
+ " var y = canvas_pos.y * mpl.ratio;\n",
+ "\n",
+ " this.send_message(name, {x: x, y: y, button: event.button,\n",
+ " step: event.step,\n",
+ " guiEvent: simpleKeys(event)});\n",
+ "\n",
+ " /* This prevents the web browser from automatically changing to\n",
+ " * the text insertion cursor when the button is pressed. We want\n",
+ " * to control all of the cursor setting manually through the\n",
+ " * 'cursor' event from matplotlib */\n",
+ " event.preventDefault();\n",
+ " return false;\n",
+ "}\n",
+ "\n",
+ "mpl.figure.prototype._key_event_extra = function(event, name) {\n",
+ " // Handle any extra behaviour associated with a key event\n",
+ "}\n",
+ "\n",
+ "mpl.figure.prototype.key_event = function(event, name) {\n",
+ "\n",
+ " // Prevent repeat events\n",
+ " if (name == 'key_press')\n",
+ " {\n",
+ " if (event.which === this._key)\n",
+ " return;\n",
+ " else\n",
+ " this._key = event.which;\n",
+ " }\n",
+ " if (name == 'key_release')\n",
+ " this._key = null;\n",
+ "\n",
+ " var value = '';\n",
+ " if (event.ctrlKey && event.which != 17)\n",
+ " value += \"ctrl+\";\n",
+ " if (event.altKey && event.which != 18)\n",
+ " value += \"alt+\";\n",
+ " if (event.shiftKey && event.which != 16)\n",
+ " value += \"shift+\";\n",
+ "\n",
+ " value += 'k';\n",
+ " value += event.which.toString();\n",
+ "\n",
+ " this._key_event_extra(event, name);\n",
+ "\n",
+ " this.send_message(name, {key: value,\n",
+ " guiEvent: simpleKeys(event)});\n",
+ " return false;\n",
+ "}\n",
+ "\n",
+ "mpl.figure.prototype.toolbar_button_onclick = function(name) {\n",
+ " if (name == 'download') {\n",
+ " this.handle_save(this, null);\n",
+ " } else {\n",
+ " this.send_message(\"toolbar_button\", {name: name});\n",
+ " }\n",
+ "};\n",
+ "\n",
+ "mpl.figure.prototype.toolbar_button_onmouseover = function(tooltip) {\n",
+ " this.message.textContent = tooltip;\n",
+ "};\n",
+ "mpl.toolbar_items = [[\"Home\", \"Reset original view\", \"fa fa-home icon-home\", \"home\"], [\"Back\", \"Back to previous view\", \"fa fa-arrow-left icon-arrow-left\", \"back\"], [\"Forward\", \"Forward to next view\", \"fa fa-arrow-right icon-arrow-right\", \"forward\"], [\"\", \"\", \"\", \"\"], [\"Pan\", \"Pan axes with left mouse, zoom with right\", \"fa fa-arrows icon-move\", \"pan\"], [\"Zoom\", \"Zoom to rectangle\", \"fa fa-square-o icon-check-empty\", \"zoom\"], [\"\", \"\", \"\", \"\"], [\"Download\", \"Download plot\", \"fa fa-floppy-o icon-save\", \"download\"]];\n",
+ "\n",
+ "mpl.extensions = [\"eps\", \"pdf\", \"png\", \"ps\", \"raw\", \"svg\"];\n",
+ "\n",
+ "mpl.default_extension = \"png\";var comm_websocket_adapter = function(comm) {\n",
+ " // Create a \"websocket\"-like object which calls the given IPython comm\n",
+ " // object with the appropriate methods. Currently this is a non binary\n",
+ " // socket, so there is still some room for performance tuning.\n",
+ " var ws = {};\n",
+ "\n",
+ " ws.close = function() {\n",
+ " comm.close()\n",
+ " };\n",
+ " ws.send = function(m) {\n",
+ " //console.log('sending', m);\n",
+ " comm.send(m);\n",
+ " };\n",
+ " // Register the callback with on_msg.\n",
+ " comm.on_msg(function(msg) {\n",
+ " //console.log('receiving', msg['content']['data'], msg);\n",
+ " // Pass the mpl event to the overridden (by mpl) onmessage function.\n",
+ " ws.onmessage(msg['content']['data'])\n",
+ " });\n",
+ " return ws;\n",
+ "}\n",
+ "\n",
+ "mpl.mpl_figure_comm = function(comm, msg) {\n",
+ " // This is the function which gets called when the mpl process\n",
+ " // starts-up an IPython Comm through the \"matplotlib\" channel.\n",
+ "\n",
+ " var id = msg.content.data.id;\n",
+ " // Get hold of the div created by the display call when the Comm\n",
+ " // socket was opened in Python.\n",
+ " var element = $(\"#\" + id);\n",
+ " var ws_proxy = comm_websocket_adapter(comm)\n",
+ "\n",
+ " function ondownload(figure, format) {\n",
+ " window.open(figure.imageObj.src);\n",
+ " }\n",
+ "\n",
+ " var fig = new mpl.figure(id, ws_proxy,\n",
+ " ondownload,\n",
+ " element.get(0));\n",
+ "\n",
+ " // Call onopen now - mpl needs it, as it is assuming we've passed it a real\n",
+ " // web socket which is closed, not our websocket->open comm proxy.\n",
+ " ws_proxy.onopen();\n",
+ "\n",
+ " fig.parent_element = element.get(0);\n",
+ " fig.cell_info = mpl.find_output_cell(\"\");\n",
+ " if (!fig.cell_info) {\n",
+ " console.error(\"Failed to find cell for figure\", id, fig);\n",
+ " return;\n",
+ " }\n",
+ "\n",
+ " var output_index = fig.cell_info[2]\n",
+ " var cell = fig.cell_info[0];\n",
+ "\n",
+ "};\n",
+ "\n",
+ "mpl.figure.prototype.handle_close = function(fig, msg) {\n",
+ " var width = fig.canvas.width/mpl.ratio\n",
+ " fig.root.unbind('remove')\n",
+ "\n",
+ " // Update the output cell to use the data from the current canvas.\n",
+ " fig.push_to_output();\n",
+ " var dataURL = fig.canvas.toDataURL();\n",
+ " // Re-enable the keyboard manager in IPython - without this line, in FF,\n",
+ " // the notebook keyboard shortcuts fail.\n",
+ " IPython.keyboard_manager.enable()\n",
+ " $(fig.parent_element).html('
');\n",
+ " fig.close_ws(fig, msg);\n",
+ "}\n",
+ "\n",
+ "mpl.figure.prototype.close_ws = function(fig, msg){\n",
+ " fig.send_message('closing', msg);\n",
+ " // fig.ws.close()\n",
+ "}\n",
+ "\n",
+ "mpl.figure.prototype.push_to_output = function(remove_interactive) {\n",
+ " // Turn the data on the canvas into data in the output cell.\n",
+ " var width = this.canvas.width/mpl.ratio\n",
+ " var dataURL = this.canvas.toDataURL();\n",
+ " this.cell_info[1]['text/html'] = '
';\n",
+ "}\n",
+ "\n",
+ "mpl.figure.prototype.updated_canvas_event = function() {\n",
+ " // Tell IPython that the notebook contents must change.\n",
+ " IPython.notebook.set_dirty(true);\n",
+ " this.send_message(\"ack\", {});\n",
+ " var fig = this;\n",
+ " // Wait a second, then push the new image to the DOM so\n",
+ " // that it is saved nicely (might be nice to debounce this).\n",
+ " setTimeout(function () { fig.push_to_output() }, 1000);\n",
+ "}\n",
+ "\n",
+ "mpl.figure.prototype._init_toolbar = function() {\n",
+ " var fig = this;\n",
+ "\n",
+ " var nav_element = $('')\n",
+ " nav_element.attr('style', 'width: 100%');\n",
+ " this.root.append(nav_element);\n",
+ "\n",
+ " // Define a callback function for later on.\n",
+ " function toolbar_event(event) {\n",
+ " return fig.toolbar_button_onclick(event['data']);\n",
+ " }\n",
+ " function toolbar_mouse_event(event) {\n",
+ " return fig.toolbar_button_onmouseover(event['data']);\n",
+ " }\n",
+ "\n",
+ " for(var toolbar_ind in mpl.toolbar_items){\n",
+ " var name = mpl.toolbar_items[toolbar_ind][0];\n",
+ " var tooltip = mpl.toolbar_items[toolbar_ind][1];\n",
+ " var image = mpl.toolbar_items[toolbar_ind][2];\n",
+ " var method_name = mpl.toolbar_items[toolbar_ind][3];\n",
+ "\n",
+ " if (!name) { continue; };\n",
+ "\n",
+ " var button = $('');\n",
+ " button.click(method_name, toolbar_event);\n",
+ " button.mouseover(tooltip, toolbar_mouse_event);\n",
+ " nav_element.append(button);\n",
+ " }\n",
+ "\n",
+ " // Add the status bar.\n",
+ " var status_bar = $('');\n",
+ " nav_element.append(status_bar);\n",
+ " this.message = status_bar[0];\n",
+ "\n",
+ " // Add the close button to the window.\n",
+ " var buttongrp = $('');\n",
+ " var button = $('');\n",
+ " button.click(function (evt) { fig.handle_close(fig, {}); } );\n",
+ " button.mouseover('Stop Interaction', toolbar_mouse_event);\n",
+ " buttongrp.append(button);\n",
+ " var titlebar = this.root.find($('.ui-dialog-titlebar'));\n",
+ " titlebar.prepend(buttongrp);\n",
+ "}\n",
+ "\n",
+ "mpl.figure.prototype._root_extra_style = function(el){\n",
+ " var fig = this\n",
+ " el.on(\"remove\", function(){\n",
+ "\tfig.close_ws(fig, {});\n",
+ " });\n",
+ "}\n",
+ "\n",
+ "mpl.figure.prototype._canvas_extra_style = function(el){\n",
+ " // this is important to make the div 'focusable\n",
+ " el.attr('tabindex', 0)\n",
+ " // reach out to IPython and tell the keyboard manager to turn it's self\n",
+ " // off when our div gets focus\n",
+ "\n",
+ " // location in version 3\n",
+ " if (IPython.notebook.keyboard_manager) {\n",
+ " IPython.notebook.keyboard_manager.register_events(el);\n",
+ " }\n",
+ " else {\n",
+ " // location in version 2\n",
+ " IPython.keyboard_manager.register_events(el);\n",
+ " }\n",
+ "\n",
+ "}\n",
+ "\n",
+ "mpl.figure.prototype._key_event_extra = function(event, name) {\n",
+ " var manager = IPython.notebook.keyboard_manager;\n",
+ " if (!manager)\n",
+ " manager = IPython.keyboard_manager;\n",
+ "\n",
+ " // Check for shift+enter\n",
+ " if (event.shiftKey && event.which == 13) {\n",
+ " this.canvas_div.blur();\n",
+ " event.shiftKey = false;\n",
+ " // Send a \"J\" for go to next cell\n",
+ " event.which = 74;\n",
+ " event.keyCode = 74;\n",
+ " manager.command_mode();\n",
+ " manager.handle_keydown(event);\n",
+ " }\n",
+ "}\n",
+ "\n",
+ "mpl.figure.prototype.handle_save = function(fig, msg) {\n",
+ " fig.ondownload(fig, null);\n",
+ "}\n",
+ "\n",
+ "\n",
+ "mpl.find_output_cell = function(html_output) {\n",
+ " // Return the cell and output element which can be found *uniquely* in the notebook.\n",
+ " // Note - this is a bit hacky, but it is done because the \"notebook_saving.Notebook\"\n",
+ " // IPython event is triggered only after the cells have been serialised, which for\n",
+ " // our purposes (turning an active figure into a static one), is too late.\n",
+ " var cells = IPython.notebook.get_cells();\n",
+ " var ncells = cells.length;\n",
+ " for (var i=0; i= 3 moved mimebundle to data attribute of output\n",
+ " data = data.data;\n",
+ " }\n",
+ " if (data['text/html'] == html_output) {\n",
+ " return [cell, data, j];\n",
+ " }\n",
+ " }\n",
+ " }\n",
+ " }\n",
+ "}\n",
+ "\n",
+ "// Register the function which deals with the matplotlib target/channel.\n",
+ "// The kernel may be null if the page has been refreshed.\n",
+ "if (IPython.notebook.kernel != null) {\n",
+ " IPython.notebook.kernel.comm_manager.register_target('matplotlib', mpl.mpl_figure_comm);\n",
+ "}\n"
+ ],
+ "text/plain": [
+ ""
+ ]
+ },
+ "metadata": {},
+ "output_type": "display_data"
+ },
+ {
+ "data": {
+ "text/html": [
+ "
"
+ ],
+ "text/plain": [
+ ""
+ ]
+ },
+ "metadata": {},
+ "output_type": "display_data"
+ },
+ {
+ "data": {
+ "text/plain": [
+ ""
+ ]
+ },
+ "execution_count": 27,
+ "metadata": {},
+ "output_type": "execute_result"
+ }
+ ],
"source": [
"filenames = dict(dat='data.grp', bkg='background.grp', res='resolution.grp')\n",
"# Check files exists\n",
@@ -191,10 +976,154 @@
"dat = load_dave(pjn(data_dir, 'data.grp'))\n",
"bkg = load_dave(pjn(data_dir, 'resolution.grp'))\n",
"\n",
- "f, (ax1, ax2, ax3) = plt.subplots(1, 3)\n",
- "[ax1.semilogy(res['x'], spectrum) for spectrum in res['y']]; ax1.set_title('resolution')\n",
- "[ax2.semilogy(dat['x'], spectrum) for spectrum in dat['y']]; ax2.set_title('sample')\n",
- "[ax3.semilogy(bkg['x'], spectrum) for spectrum in bkg['y']]; ax3.set_title('empty can')"
+ "plot_args=dict(color='black', marker='o', markersize=1, markerfacecolor='none')\n",
+ "f, ax = plt.subplots()\n",
+ "ax.scatter(res['x'], res['y'][7])\n",
+ "#ax.semilogy(res['x'], res['y'][7], **plot_args)\n",
+ "#[ax.semilogy(res['x'], spectrum) for spectrum in res['y']]; ax.set_title('resolution')\n",
+ "#f, (ax1, ax2, ax3) = plt.subplots(1, 3)\n",
+ "#[ax1.semilogy(res['x'], spectrum) for spectrum in res['y']]; ax1.set_title('resolution')\n",
+ "#[ax2.semilogy(dat['x'], spectrum) for spectrum in dat['y']]; ax2.set_title('sample')\n",
+ "#[ax3.semilogy(bkg['x'], spectrum) for spectrum in bkg['y']]; ax3.set_title('empty can')"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "execution_count": 9,
+ "metadata": {},
+ "outputs": [
+ {
+ "data": {
+ "text/plain": [
+ "array([2.07352e-03, 1.98801e-03, 2.40648e-03, 2.11926e-03, 1.94248e-03,\n",
+ " 2.65906e-03, 2.74229e-03, 1.97901e-03, 2.49010e-03, 1.90315e-03,\n",
+ " 1.90426e-03, 2.77858e-03, 2.14998e-03, 2.10387e-03, 2.36193e-03,\n",
+ " 2.12497e-03, 2.28415e-03, 2.49503e-03, 2.36365e-03, 2.22517e-03,\n",
+ " 2.27831e-03, 2.07177e-03, 2.15702e-03, 2.49062e-03, 1.97939e-03,\n",
+ " 2.28361e-03, 2.44326e-03, 2.23229e-03, 1.99061e-03, 1.68617e-03,\n",
+ " 2.28203e-03, 1.64718e-03, 2.16069e-03, 2.65759e-03, 1.88809e-03,\n",
+ " 2.31075e-03, 2.11627e-03, 2.24369e-03, 1.73030e-03, 1.92596e-03,\n",
+ " 2.78282e-03, 1.47648e-03, 2.16421e-03, 2.32116e-03, 1.98224e-03,\n",
+ " 2.44584e-03, 2.49590e-03, 2.73845e-03, 1.85247e-03, 1.97539e-03,\n",
+ " 2.03021e-03, 2.91694e-03, 2.18518e-03, 2.09944e-03, 2.32502e-03,\n",
+ " 2.06023e-03, 2.44385e-03, 2.66492e-03, 2.36131e-03, 2.27538e-03,\n",
+ " 2.02805e-03, 2.11112e-03, 3.04458e-03, 2.27854e-03, 2.19993e-03,\n",
+ " 2.48600e-03, 2.34966e-03, 2.78541e-03, 2.45633e-03, 2.02416e-03,\n",
+ " 2.06155e-03, 2.27499e-03, 2.53365e-03, 2.23751e-03, 2.44426e-03,\n",
+ " 2.09845e-03, 1.98645e-03, 1.76600e-03, 1.93790e-03, 1.42635e-03,\n",
+ " 2.02145e-03, 2.79043e-03, 2.45066e-03, 2.57754e-03, 2.39610e-03,\n",
+ " 2.95587e-03, 2.05948e-03, 2.23060e-03, 1.81066e-03, 2.47072e-03,\n",
+ " 2.39268e-03, 2.48634e-03, 2.27539e-03, 2.69679e-03, 2.92821e-03,\n",
+ " 2.64264e-03, 2.05084e-03, 2.78656e-03, 2.27340e-03, 2.26677e-03,\n",
+ " 2.88654e-03, 2.13273e-03, 2.68241e-03, 2.24374e-03, 2.14203e-03,\n",
+ " 1.86914e-03, 2.60970e-03, 2.55469e-03, 2.67880e-03, 2.77404e-03,\n",
+ " 2.10167e-03, 2.35183e-03, 2.75417e-03, 2.43942e-03, 1.86918e-03,\n",
+ " 2.17290e-03, 2.68981e-03, 2.58961e-03, 2.52443e-03, 1.98795e-03,\n",
+ " 2.83062e-03, 2.71118e-03, 2.38238e-03, 2.06025e-03, 2.61622e-03,\n",
+ " 2.57453e-03, 2.69858e-03, 2.77589e-03, 2.97541e-03, 2.95205e-03,\n",
+ " 2.45388e-03, 2.45735e-03, 2.39754e-03, 2.51626e-03, 2.79324e-03,\n",
+ " 2.65575e-03, 2.66312e-03, 2.53157e-03, 2.58195e-03, 2.53595e-03,\n",
+ " 2.39898e-03, 2.70431e-03, 2.50084e-03, 2.94866e-03, 3.17008e-03,\n",
+ " 1.90523e-03, 2.62332e-03, 2.85942e-03, 2.57729e-03, 2.32550e-03,\n",
+ " 2.99415e-03, 2.72313e-03, 2.98580e-03, 2.99370e-03, 2.78571e-03,\n",
+ " 3.01404e-03, 3.80445e-03, 2.24861e-03, 2.89206e-03, 2.24179e-03,\n",
+ " 2.66175e-03, 3.12115e-03, 2.19479e-03, 2.53968e-03, 2.81251e-03,\n",
+ " 2.95061e-03, 3.04304e-03, 3.18596e-03, 3.05680e-03, 2.38256e-03,\n",
+ " 3.02446e-03, 3.57089e-03, 2.59584e-03, 2.59105e-03, 2.75340e-03,\n",
+ " 2.69895e-03, 2.46750e-03, 2.66874e-03, 2.92053e-03, 2.45609e-03,\n",
+ " 2.07409e-03, 3.00083e-03, 2.57886e-03, 2.51004e-03, 2.44170e-03,\n",
+ " 2.85621e-03, 2.96456e-03, 3.05431e-03, 2.70712e-03, 3.14910e-03,\n",
+ " 2.68484e-03, 2.36652e-03, 3.45014e-03, 3.10367e-03, 3.03759e-03,\n",
+ " 3.23444e-03, 3.37453e-03, 2.77932e-03, 3.10750e-03, 3.23938e-03,\n",
+ " 2.97690e-03, 3.20192e-03, 3.45549e-03, 3.37409e-03, 3.75603e-03,\n",
+ " 3.19663e-03, 4.42901e-03, 3.23041e-03, 4.15401e-03, 5.82597e-03,\n",
+ " 5.05707e-03, 5.28717e-03, 5.42936e-03, 5.77280e-03, 5.22316e-03,\n",
+ " 6.01651e-03, 7.01220e-03, 7.47903e-03, 7.31737e-03, 7.78484e-03,\n",
+ " 7.11050e-03, 9.29424e-03, 9.83936e-03, 9.67883e-03, 9.55568e-03,\n",
+ " 9.63449e-03, 1.11162e-02, 1.32040e-02, 1.32394e-02, 1.40030e-02,\n",
+ " 1.50712e-02, 1.48919e-02, 1.53750e-02, 1.62214e-02, 1.68546e-02,\n",
+ " 1.88030e-02, 1.97669e-02, 2.21957e-02, 2.13952e-02, 2.19078e-02,\n",
+ " 2.33465e-02, 2.41731e-02, 2.72144e-02, 2.61240e-02, 2.96669e-02,\n",
+ " 3.14108e-02, 3.37740e-02, 3.35188e-02, 3.49078e-02, 3.78900e-02,\n",
+ " 3.95885e-02, 4.15876e-02, 4.41046e-02, 4.68040e-02, 4.59056e-02,\n",
+ " 5.27376e-02, 5.58091e-02, 5.58518e-02, 6.10066e-02, 6.42855e-02,\n",
+ " 6.47374e-02, 7.09028e-02, 7.48037e-02, 8.15661e-02, 8.40971e-02,\n",
+ " 8.85314e-02, 1.01524e-01, 1.06127e-01, 1.16434e-01, 1.18878e-01,\n",
+ " 1.32713e-01, 1.46167e-01, 1.57076e-01, 1.82480e-01, 2.03719e-01,\n",
+ " 2.28850e-01, 2.56709e-01, 2.86639e-01, 3.35071e-01, 3.93904e-01,\n",
+ " 4.39578e-01, 5.23814e-01, 6.25825e-01, 7.47127e-01, 9.12108e-01,\n",
+ " 1.11099e+00, 1.37497e+00, 1.66811e+00, 2.08018e+00, 2.65054e+00,\n",
+ " 3.35505e+00, 4.22075e+00, 5.37934e+00, 6.91791e+00, 8.83160e+00,\n",
+ " 1.12629e+01, 1.41890e+01, 1.71804e+01, 1.95246e+01, 2.02338e+01,\n",
+ " 1.86709e+01, 1.51660e+01, 1.05711e+01, 6.32952e+00, 3.22560e+00,\n",
+ " 1.38432e+00, 5.08688e-01, 1.71896e-01, 6.10896e-02, 3.29144e-02,\n",
+ " 2.49313e-02, 2.31227e-02, 1.89193e-02, 1.87901e-02, 1.85200e-02,\n",
+ " 1.82864e-02, 1.66453e-02, 1.59132e-02, 1.44735e-02, 1.41853e-02,\n",
+ " 1.36843e-02, 1.37253e-02, 1.36083e-02, 1.27431e-02, 1.20175e-02,\n",
+ " 1.07525e-02, 1.07982e-02, 1.16278e-02, 1.08813e-02, 1.09815e-02,\n",
+ " 1.03583e-02, 1.04614e-02, 1.04650e-02, 9.86538e-03, 9.04935e-03,\n",
+ " 8.65516e-03, 9.15849e-03, 8.48274e-03, 8.88216e-03, 8.28317e-03,\n",
+ " 8.47669e-03, 7.24619e-03, 7.85788e-03, 6.94695e-03, 7.05628e-03,\n",
+ " 6.77768e-03, 5.74367e-03, 5.90941e-03, 5.30114e-03, 5.40820e-03,\n",
+ " 5.82286e-03, 5.23451e-03, 5.41669e-03, 4.95516e-03, 4.94773e-03,\n",
+ " 3.90986e-03, 4.77555e-03, 4.75833e-03, 5.70458e-03, 4.48677e-03,\n",
+ " 4.56655e-03, 4.68920e-03, 4.39629e-03, 3.43943e-03, 4.25240e-03,\n",
+ " 3.36397e-03, 3.77910e-03, 3.94559e-03, 3.73069e-03, 3.25141e-03,\n",
+ " 3.81802e-03, 3.33750e-03, 4.65665e-03, 3.53463e-03, 3.17860e-03,\n",
+ " 3.33898e-03, 3.39096e-03, 3.35002e-03, 3.49346e-03, 3.17282e-03,\n",
+ " 3.46696e-03, 3.74702e-03, 3.64913e-03, 2.99306e-03, 3.22060e-03,\n",
+ " 3.22626e-03, 3.04175e-03, 3.05020e-03, 2.60630e-03, 2.73593e-03,\n",
+ " 3.25625e-03, 2.72771e-03, 2.74346e-03, 2.66109e-03, 3.31551e-03,\n",
+ " 2.63189e-03, 2.91678e-03, 2.67406e-03, 2.69288e-03, 2.99564e-03,\n",
+ " 2.99774e-03, 2.86850e-03, 3.35939e-03, 3.22409e-03, 2.69796e-03,\n",
+ " 2.96511e-03, 2.34836e-03, 2.48178e-03, 2.61465e-03, 2.97272e-03,\n",
+ " 3.36323e-03, 2.56860e-03, 2.74271e-03, 3.53368e-03, 3.03068e-03,\n",
+ " 2.56234e-03, 1.86588e-03, 3.39876e-03, 3.08911e-03, 3.05531e-03,\n",
+ " 2.74532e-03, 2.83829e-03, 2.57814e-03, 2.61691e-03, 2.79071e-03,\n",
+ " 2.65320e-03, 2.78391e-03, 2.35336e-03, 2.05396e-03, 2.95993e-03,\n",
+ " 3.18196e-03, 2.64338e-03, 3.01607e-03, 2.92460e-03, 2.27586e-03,\n",
+ " 2.79255e-03, 3.09374e-03, 2.66650e-03, 2.52324e-03, 2.74674e-03,\n",
+ " 2.71370e-03, 2.62634e-03, 2.06339e-03, 2.97489e-03, 2.40118e-03,\n",
+ " 3.50194e-03, 3.59058e-03, 2.77072e-03, 2.91274e-03, 2.47781e-03,\n",
+ " 2.89326e-03, 2.73156e-03, 2.56152e-03, 2.53920e-03, 2.86390e-03,\n",
+ " 2.54181e-03, 2.86398e-03, 2.81895e-03, 2.64348e-03, 2.88787e-03,\n",
+ " 3.05978e-03, 2.68863e-03, 2.51372e-03, 2.19842e-03, 3.00046e-03,\n",
+ " 2.72275e-03, 2.68651e-03, 2.49482e-03, 2.54159e-03, 3.24878e-03,\n",
+ " 3.38886e-03, 3.37550e-03, 2.63820e-03, 2.45954e-03, 2.64090e-03,\n",
+ " 2.54267e-03, 2.67391e-03, 2.50333e-03, 2.57664e-03, 2.77382e-03,\n",
+ " 2.44347e-03, 2.15113e-03, 2.82968e-03, 2.51458e-03, 2.19794e-03,\n",
+ " 3.86735e-03, 2.33218e-03, 2.73011e-03, 2.66817e-03, 2.70032e-03,\n",
+ " 2.82107e-03, 2.95210e-03, 2.75802e-03, 2.08020e-03, 2.48008e-03,\n",
+ " 2.18425e-03, 2.89944e-03, 2.71299e-03, 2.66804e-03, 2.92648e-03,\n",
+ " 2.72186e-03, 2.49946e-03, 1.85661e-03, 2.33070e-03, 2.58615e-03,\n",
+ " 2.54713e-03, 2.54148e-03, 2.66770e-03, 2.67363e-03, 3.21592e-03,\n",
+ " 2.22920e-03, 2.07532e-03, 3.07635e-03, 2.62717e-03, 2.58792e-03,\n",
+ " 2.94388e-03, 2.40489e-03, 3.12705e-03, 2.84948e-03, 2.17925e-03,\n",
+ " 2.80748e-03, 2.63197e-03, 2.53013e-03, 2.84187e-03, 2.56897e-03,\n",
+ " 2.41366e-03, 2.53796e-03, 2.40467e-03, 2.22654e-03, 2.35061e-03,\n",
+ " 2.08604e-03, 3.06315e-03, 2.26903e-03, 2.50106e-03, 2.90944e-03,\n",
+ " 2.84619e-03, 2.31831e-03, 2.13081e-03, 2.26269e-03, 2.71843e-03,\n",
+ " 2.72076e-03, 2.49025e-03, 2.88677e-03, 2.75622e-03, 2.80286e-03,\n",
+ " 2.27806e-03, 2.62367e-03, 2.85965e-03, 2.75378e-03, 2.45441e-03,\n",
+ " 2.35238e-03, 3.24105e-03, 2.48270e-03, 2.64116e-03, 2.63295e-03,\n",
+ " 2.45561e-03, 2.72581e-03, 2.33879e-03, 2.27319e-03, 2.59789e-03,\n",
+ " 2.72411e-03, 2.74688e-03, 2.66762e-03, 2.95966e-03, 2.02116e-03,\n",
+ " 2.84203e-03, 2.29798e-03, 1.53611e-03, 2.40691e-03, 2.38609e-03,\n",
+ " 2.61860e-03, 3.39788e-03, 2.58306e-03, 3.05087e-03, 2.48685e-03,\n",
+ " 2.27617e-03, 2.95029e-03, 3.07029e-03, 2.43194e-03, 2.59492e-03,\n",
+ " 2.58083e-03, 2.82092e-03, 3.25174e-03, 2.83641e-03, 2.76988e-03,\n",
+ " 2.29295e-03, 2.44398e-03, 2.82941e-03, 2.32943e-03, 2.82407e-03,\n",
+ " 2.31433e-03, 2.98565e-03, 2.00406e-03, 1.97199e-03, 3.12078e-03,\n",
+ " 2.62855e-03, 2.83925e-03, 2.62244e-03, 2.38216e-03, 2.25435e-03])"
+ ]
+ },
+ "execution_count": 9,
+ "metadata": {},
+ "output_type": "execute_result"
+ }
+ ],
+ "source": [
+ "res['y'][0]"
]
},
{