{ "cells": [ { "cell_type": "markdown", "metadata": {}, "source": [ "# Neutron Powder Diffraction\n", "\n", "In this tutorial demonstrates how neutron-scattering data can be loaded, visualized, and manipulated with generic functionality from `scipp` as well as neutron-specific functionality from `scippneutron`. It focuses on reducing data from the ORNL [POWGEN](https://neutrons.ornl.gov/powgen) neutron diffractometer." ] }, { "cell_type": "code", "execution_count": 1, "metadata": {}, "outputs": [], "source": [ "import numpy as np\n", "import scipp as sc\n", "import scippneutron as scn" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "### Loading Nexus files\n", "\n", "Loading Nexus files requires [Mantid](https://www.mantidproject.org).\n", "See, e.g., [Installation](https://scipp.github.io/getting-started/installation.html) on how to install scipp and Mantid with `conda`.\n", "We are using two files in this tutorial,\n", "[PG3_4844_event.nxs](http://198.74.56.37/ftp/external-data/MD5/d5ae38871d0a09a28ae01f85d969de1e)\n", "and\n", "[PG3_4866_event.nxs](http://198.74.56.37/ftp/external-data/MD5/3d543bc6a646e622b3f4542bc3435e7e).\n", "Both are available as part of Mantid's test data.\n", "\n", "Rename the files upon download.\n", "\n", "We start by loading two files: the sample and the vanadium runs." ] }, { "cell_type": "code", "execution_count": 2, "metadata": {}, "outputs": [ { "name": "stderr", "output_type": "stream", "text": [ "Downloading file 'PG3_4844_event.nxs' from 'https://public.esss.dk/groups/scipp/scippneutron/2/PG3_4844_event.nxs' to '/home/scientist/.cache/scippneutron/2'.\n" ] }, { "name": "stdout", "output_type": "stream", "text": [ "FrameworkManager-[Notice] Welcome to Mantid 6.2.20211217.1715\n", "FrameworkManager-[Notice] Please cite: http://dx.doi.org/10.1016/j.nima.2014.07.029 and this release: http://dx.doi.org/10.5286/Software/Mantid\n", "DownloadInstrument-[Notice] All instrument definitions up to date\n", "Load-[Notice] Load started\n", "Load-[Notice] Load successful, Duration 2.77 seconds\n" ] }, { "name": "stderr", "output_type": "stream", "text": [ "Downloading file 'PG3_4866_event.nxs' from 'https://public.esss.dk/groups/scipp/scippneutron/2/PG3_4866_event.nxs' to '/home/scientist/.cache/scippneutron/2'.\n" ] }, { "name": "stdout", "output_type": "stream", "text": [ "ExtractSpectra-[Notice] ExtractSpectra started\n", "ExtractSpectra-[Notice] ExtractSpectra successful, Duration 0.00 seconds\n", "DeleteWorkspace-[Notice] DeleteWorkspace started\n", "DeleteWorkspace-[Notice] DeleteWorkspace successful, Duration 0.00 seconds\n", "DeleteWorkspace-[Notice] DeleteWorkspace started\n", "DeleteWorkspace-[Notice] DeleteWorkspace successful, Duration 0.00 seconds\n", "DeleteWorkspace-[Notice] DeleteWorkspace started\n", "DeleteWorkspace-[Notice] DeleteWorkspace successful, Duration 0.00 seconds\n", "Load-[Notice] Load started\n", "Load-[Notice] Load successful, Duration 7.34 seconds\n", "DeleteWorkspace-[Notice] DeleteWorkspace started\n", "DeleteWorkspace-[Notice] DeleteWorkspace successful, Duration 0.00 seconds\n" ] } ], "source": [ "sample = scn.load(scn.data.get_path('PG3_4844_event.nxs'),\n", " load_pulse_times=False,\n", " mantid_args={'LoadMonitors': True})\n", "vanadium = scn.load(scn.data.get_path('PG3_4866_event.nxs'),\n", " load_pulse_times=False)" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "The optional `mantid_args` dict is forwarded to the Mantid algorithm used for loading the files – in this case [LoadEventNexus](https://docs.mantidproject.org/nightly/algorithms/LoadEventNexus-v1.html) – and can be used to control, e.g., which part of a file to load.\n", "Here we request loading monitors, which Mantid does not load by default.\n", "The resulting dataset looks as follows:" ] }, { "cell_type": "code", "execution_count": 3, "metadata": {}, "outputs": [ { "data": { "text/html": [ "
\n", "\n", "\n", "\n", "\n", "Show/Hide data repr\n", "\n", "\n", "\n", "\n", "\n", "Show/Hide attributes\n", "\n", "\n", "\n", "\n", "\n", "\n", "
scipp.DataArray (274.59 MB)
" ], "text/plain": [ "\n", "Dimensions: Sizes[spectrum:24794, tof:1, ]\n", "Coordinates:\n", " position vector3 [m] (spectrum) [(1.17451, -1.01106, -2.03797), (1.18148, -0.959466, -2.05334), ..., (1.81429, 0.0956584, 3.84338), (1.81375, 0.149937, 3.8427)]\n", " sample_position vector3 [m] () [(0, 0, 0)]\n", " source_position vector3 [m] () [(0, 0, -60)]\n", " spectrum int32 [dimensionless] (spectrum) [1, 2, ..., 24793, 24794]\n", " tof float64 [µs] (tof [bin-edge]) [19, 16694.3]\n", "Data:\n", " DataArrayView [dimensionless] (spectrum, tof) [\n", "Dimensions: Sizes[event:0, ]\n", "Coordinates:\n", " tof float64 [µs] (event) []\n", "Data:\n", " float32 [counts] (event) [] []\n", "\n", ", \n", "Dimensions: Sizes[event:0, ]\n", "Coordinates:\n", " tof float64 [µs] (event) []\n", "Data:\n", " float32 [counts] (event) [] []\n", "\n", ", ..., \n", "Dimensions: Sizes[event:0, ]\n", "Coordinates:\n", " tof float64 [µs] (event) []\n", "Data:\n", " float32 [counts] (event) [] []\n", "\n", ", \n", "Dimensions: Sizes[event:0, ]\n", "Coordinates:\n", " tof float64 [µs] (event) []\n", "Data:\n", " float32 [counts] (event) [] []\n", "\n", "]\n", "Attributes:\n", " ChopperStatus1 DataArray [dimensionless] () [\n", "Dimensions: Sizes[time:2, ]\n", "Coordinates:\n", " time datetime64 [ns] (time) [2011-08-12T15:50:16.800000000, 2011-08-12T17:22:05.000000000]\n", "Data:\n", " float64 [dimensionless] (time) [4, 4]\n", "\n", "]\n", " ChopperStatus2 DataArray [dimensionless] () [\n", "Dimensions: Sizes[time:2, ]\n", "Coordinates:\n", " time datetime64 [ns] (time) [2011-08-12T15:50:16.816000000, 2011-08-12T17:22:05.000000000]\n", "Data:\n", " float64 [dimensionless] (time) [4, 4]\n", "\n", "]\n", " ChopperStatus3 DataArray [dimensionless] () [\n", "Dimensions: Sizes[time:2, ]\n", "Coordinates:\n", " time datetime64 [ns] (time) [2011-08-12T15:50:16.816000000, 2011-08-12T17:22:05.000000000]\n", "Data:\n", " float64 [dimensionless] (time) [4, 4]\n", "\n", "]\n", " CurrentSP DataArray [dimensionless] () [\n", "Dimensions: Sizes[time:2, ]\n", "Coordinates:\n", " time datetime64 [ns] (time) [2011-08-12T15:50:16.800000000, 2011-08-12T17:22:05.000000000]\n", "Data:\n", " float64 [dimensionless] (time) [300, 300]\n", "\n", "]\n", " EnergyRequest DataArray [dimensionless] () [\n", "Dimensions: Sizes[time:2, ]\n", "Coordinates:\n", " time datetime64 [ns] (time) [2011-08-12T15:50:16.784000000, 2011-08-12T17:22:05.000000000]\n", "Data:\n", " float64 [meV] (time) [287.955, 287.955]\n", "\n", "]\n", " Filename string [dimensionless] () [\"/home/scientist/.cache/scippneutron/2/PG3_4844_event.nxs\"]\n", " LKSRampRate DataArray [dimensionless] () [\n", "Dimensions: Sizes[time:2, ]\n", "Coordinates:\n", " time datetime64 [ns] (time) [2011-08-12T15:50:16.800000000, 2011-08-12T17:22:05.000000000]\n", "Data:\n", " float64 [dimensionless] (time) [0, 0]\n", "\n", "]\n", " LambdaRequest DataArray [dimensionless] () [\n", "Dimensions: Sizes[time:2, ]\n", "Coordinates:\n", " time datetime64 [ns] (time) [2011-08-12T15:50:16.784000000, 2011-08-12T17:22:05.000000000]\n", "Data:\n", " float64 [Å] (time) [0.533, 0.533]\n", "\n", "]\n", " Phase1 DataArray [dimensionless] () [\n", "Dimensions: Sizes[time:1794, ]\n", "Coordinates:\n", " time datetime64 [ns] (time) [2011-08-12T15:50:16.800000000, 2011-08-12T15:50:18.267000007, ..., 2011-08-12T17:22:02.999218750, 2011-08-12T17:22:05.000000000]\n", "Data:\n", " float64 [µs] (time) [8166.72, 8165.16, ..., 8163.85, 8163.85]\n", "\n", "]\n", " Phase2 DataArray [dimensionless] () [\n", "Dimensions: Sizes[time:1793, ]\n", "Coordinates:\n", " time datetime64 [ns] (time) [2011-08-12T15:50:16.800000000, 2011-08-12T15:50:18.267000007, ..., 2011-08-12T17:22:02.999218750, 2011-08-12T17:22:05.000000000]\n", "Data:\n", " float64 [µs] (time) [8335.63, 8334.09, ..., 8332.86, 8332.86]\n", "\n", "]\n", " Phase3 DataArray [dimensionless] () [\n", "Dimensions: Sizes[time:1777, ]\n", "Coordinates:\n", " time datetime64 [ns] (time) [2011-08-12T15:50:16.800000000, 2011-08-12T15:50:19.688999938, ..., 2011-08-12T17:22:04.421093750, 2011-08-12T17:22:05.000000000]\n", "Data:\n", " float64 [µs] (time) [14001.5, 14001.4, ..., 14002, 14002]\n", "\n", "]\n", " PhaseRequest1 DataArray [dimensionless] () [\n", "Dimensions: Sizes[time:2, ]\n", "Coordinates:\n", " time datetime64 [ns] (time) [2011-08-12T15:50:16.800000000, 2011-08-12T17:22:05.000000000]\n", "Data:\n", " float64 [µs] (time) [8164.07, 8164.07]\n", "\n", "]\n", " PhaseRequest2 DataArray [dimensionless] () [\n", "Dimensions: Sizes[time:2, ]\n", "Coordinates:\n", " time datetime64 [ns] (time) [2011-08-12T15:50:16.800000000, 2011-08-12T17:22:05.000000000]\n", "Data:\n", " float64 [µs] (time) [8332.89, 8332.89]\n", "\n", "]\n", " PhaseRequest3 DataArray [dimensionless] () [\n", "Dimensions: Sizes[time:2, ]\n", "Coordinates:\n", " time datetime64 [ns] (time) [2011-08-12T15:50:16.800000000, 2011-08-12T17:22:05.000000000]\n", "Data:\n", " float64 [µs] (time) [14001.8, 14001.8]\n", "\n", "]\n", " S1HCenter DataArray [dimensionless] () [\n", "Dimensions: Sizes[time:2, ]\n", "Coordinates:\n", " time datetime64 [ns] (time) [2011-08-12T15:50:16.800000000, 2011-08-12T17:22:05.000000000]\n", "Data:\n", " float64 [mm] (time) [0, 0]\n", "\n", "]\n", " S1HCenterOffset DataArray [dimensionless] () [\n", "Dimensions: Sizes[time:2, ]\n", "Coordinates:\n", " time datetime64 [ns] (time) [2011-08-12T15:50:16.800000000, 2011-08-12T17:22:05.000000000]\n", "Data:\n", " float64 [dimensionless] (time) [0, 0]\n", "\n", "]\n", " S1HWidth DataArray [dimensionless] () [\n", "Dimensions: Sizes[time:2, ]\n", "Coordinates:\n", " time datetime64 [ns] (time) [2011-08-12T15:50:16.800000000, 2011-08-12T17:22:05.000000000]\n", "Data:\n", " float64 [mm] (time) [10, 10]\n", "\n", "]\n", " S1VCenter DataArray [dimensionless] () [\n", "Dimensions: Sizes[time:2, ]\n", "Coordinates:\n", " time datetime64 [ns] (time) [2011-08-12T15:50:16.800000000, 2011-08-12T17:22:05.000000000]\n", "Data:\n", " float64 [mm] (time) [5, 5]\n", "\n", "]\n", " S1VCenterOffset DataArray [dimensionless] () [\n", "Dimensions: Sizes[time:2, ]\n", "Coordinates:\n", " time datetime64 [ns] (time) [2011-08-12T15:50:16.800000000, 2011-08-12T17:22:05.000000000]\n", "Data:\n", " float64 [dimensionless] (time) [0, 0]\n", "\n", "]\n", " S1VHeight DataArray [dimensionless] () [\n", "Dimensions: Sizes[time:2, ]\n", "Coordinates:\n", " time datetime64 [ns] (time) [2011-08-12T15:50:16.816000000, 2011-08-12T17:22:05.000000000]\n", "Data:\n", " float64 [mm] (time) [30, 30]\n", "\n", "]\n", " SampleTemp DataArray [dimensionless] () [\n", "Dimensions: Sizes[time:467, ]\n", "Coordinates:\n", " time datetime64 [ns] (time) [2011-08-12T15:50:16.784000000, 2011-08-12T15:50:18.001000007, ..., 2011-08-12T17:21:58.623843750, 2011-08-12T17:22:05.000000000]\n", "Data:\n", " float64 [dimensionless] (time) [299.352, 299.446, ..., 300, 300]\n", "\n", "]\n", " Speed1 DataArray [dimensionless] () [\n", "Dimensions: Sizes[time:2, ]\n", "Coordinates:\n", " time datetime64 [ns] (time) [2011-08-12T15:50:16.800000000, 2011-08-12T17:22:05.000000000]\n", "Data:\n", " float64 [Hz] (time) [60, 60]\n", "\n", "]\n", " Speed2 DataArray [dimensionless] () [\n", "Dimensions: Sizes[time:2, ]\n", "Coordinates:\n", " time datetime64 [ns] (time) [2011-08-12T15:50:16.800000000, 2011-08-12T17:22:05.000000000]\n", "Data:\n", " float64 [Hz] (time) [60, 60]\n", "\n", "]\n", " Speed3 DataArray [dimensionless] () [\n", "Dimensions: Sizes[time:2, ]\n", "Coordinates:\n", " time datetime64 [ns] (time) [2011-08-12T15:50:16.800000000, 2011-08-12T17:22:05.000000000]\n", "Data:\n", " float64 [Hz] (time) [60, 60]\n", "\n", "]\n", " SpeedRequest1 DataArray [dimensionless] () [\n", "Dimensions: Sizes[time:2, ]\n", "Coordinates:\n", " time datetime64 [ns] (time) [2011-08-12T15:50:16.784000000, 2011-08-12T17:22:05.000000000]\n", "Data:\n", " float64 [Hz] (time) [60, 60]\n", "\n", "]\n", " SpeedRequest2 DataArray [dimensionless] () [\n", "Dimensions: Sizes[time:2, ]\n", "Coordinates:\n", " time datetime64 [ns] (time) [2011-08-12T15:50:16.784000000, 2011-08-12T17:22:05.000000000]\n", "Data:\n", " float64 [Hz] (time) [60, 60]\n", "\n", "]\n", " SpeedRequest3 DataArray [dimensionless] () [\n", "Dimensions: Sizes[time:2, ]\n", "Coordinates:\n", " time datetime64 [ns] (time) [2011-08-12T15:50:16.784000000, 2011-08-12T17:22:05.000000000]\n", "Data:\n", " float64 [Hz] (time) [60, 60]\n", "\n", "]\n", " TolRequest DataArray [dimensionless] () [\n", "Dimensions: Sizes[time:2, ]\n", "Coordinates:\n", " time datetime64 [ns] (time) [2011-08-12T15:50:16.800000000, 2011-08-12T17:22:05.000000000]\n", "Data:\n", " float64 [dimensionless] (time) [20, 20]\n", "\n", "]\n", " commErrs DataArray [dimensionless] () [\n", "Dimensions: Sizes[time:2, ]\n", "Coordinates:\n", " time datetime64 [ns] (time) [2011-08-12T15:50:16.784000000, 2011-08-12T17:22:05.000000000]\n", "Data:\n", " float64 [dimensionless] (time) [0, 0]\n", "\n", "]\n", " currentsample DataArray [dimensionless] () [\n", "Dimensions: Sizes[time:2, ]\n", "Coordinates:\n", " time datetime64 [ns] (time) [2011-08-12T15:50:16.784000000, 2011-08-12T17:22:05.000000000]\n", "Data:\n", " float64 [dimensionless] (time) [4, 4]\n", "\n", "]\n", " duration float64 [s] () [5508]\n", " end_time string [dimensionless] () [\"2011-08-12T17:22:05\"]\n", " experiment_identifier string [dimensionless] () [\"IPTS-2767\"]\n", " fernsstatus DataArray [dimensionless] () [\n", "Dimensions: Sizes[time:2, ]\n", "Coordinates:\n", " time datetime64 [ns] (time) [2011-08-12T15:50:16.784000000, 2011-08-12T17:22:05.000000000]\n", "Data:\n", " float64 [dimensionless] (time) [3, 3]\n", "\n", "]\n", " file_notes string [dimensionless] () [\"NONE\"]\n", " frequency DataArray [dimensionless] () [\n", "Dimensions: Sizes[time:330473, ]\n", "Coordinates:\n", " time datetime64 [ns] (time) [2011-08-12T15:50:17.000000000, 2011-08-12T15:50:17.016659999, ..., 2011-08-12T17:22:05.085449218, 2011-08-12T17:22:05.102050781]\n", "Data:\n", " float64 [Hz] (time) [0, 60.024, ..., 60.0024, 59.9988]\n", "\n", "]\n", " gd_prtn_chrg float64 [µAh] () [1171.95]\n", " guide DataArray [dimensionless] () [\n", "Dimensions: Sizes[time:2, ]\n", "Coordinates:\n", " time datetime64 [ns] (time) [2011-08-12T15:50:16.784000000, 2011-08-12T17:22:05.000000000]\n", "Data:\n", " float64 [mm] (time) [-55.463, -55.463]\n", "\n", "]\n", " instrument_name string [dimensionless] () [\"POWGEN\"]\n", " monitor1 DataArray [dimensionless] () [\n", "Dimensions: Sizes[tof:200001, ]\n", "Coordinates:\n", " position vector3 [m] () [(0, 0, -1)]\n", " source_position vector3 [m] () [(0, 0, -60)]\n", " tof float64 [µs] (tof [bin-edge]) [0, 1, ..., 200000, 200001]\n", "Data:\n", " float64 [counts] (tof) [25, 10, ..., 0, 0] [25, 10, ..., 0, 0]\n", "Attributes:\n", " instrument_name string [dimensionless] () [\"POWGEN\"]\n", " sample_position vector3 [m] () [(0, 0, 0)]\n", " spectrum int32 [dimensionless] () [1]\n", "\n", "]\n", " proton_charge DataArray [dimensionless] () [\n", "Dimensions: Sizes[time:330473, ]\n", "Coordinates:\n", " time datetime64 [ns] (time) [2011-08-12T15:50:17.000000000, 2011-08-12T15:50:17.016659999, ..., 2011-08-12T17:22:05.085449218, 2011-08-12T17:22:05.102050781]\n", "Data:\n", " float64 [pC] (time) [1.48434e+07, 1.48434e+07, ..., 1.48733e+07, 1.48391e+07]\n", "\n", "]\n", " run_number string [dimensionless] () [\"4844\"]\n", " run_start string [dimensionless] () [\"2011-08-12T15:50:17\"]\n", " run_title string [dimensionless] () [\"diamond cw0.533 4.22e12 60Hz [10x30]\"]\n", " running DataArray [dimensionless] () [\n", "Dimensions: Sizes[time:1, ]\n", "Coordinates:\n", " time datetime64 [ns] (time) [1990-01-01T00:00:00.000000000]\n", "Data:\n", " bool [dimensionless] (time) [True]\n", "\n", "]\n", " s1b DataArray [dimensionless] () [\n", "Dimensions: Sizes[time:2, ]\n", "Coordinates:\n", " time datetime64 [ns] (time) [2011-08-12T15:50:16.784000000, 2011-08-12T17:22:05.000000000]\n", "Data:\n", " float64 [mm] (time) [20, 20]\n", "\n", "]\n", " s1l DataArray [dimensionless] () [\n", "Dimensions: Sizes[time:2, ]\n", "Coordinates:\n", " time datetime64 [ns] (time) [2011-08-12T15:50:16.800000000, 2011-08-12T17:22:05.000000000]\n", "Data:\n", " float64 [mm] (time) [5, 5]\n", "\n", "]\n", " s1r DataArray [dimensionless] () [\n", "Dimensions: Sizes[time:2, ]\n", "Coordinates:\n", " time datetime64 [ns] (time) [2011-08-12T15:50:16.784000000, 2011-08-12T17:22:05.000000000]\n", "Data:\n", " float64 [mm] (time) [-5, -5]\n", "\n", "]\n", " s1t DataArray [dimensionless] () [\n", "Dimensions: Sizes[time:2, ]\n", "Coordinates:\n", " time datetime64 [ns] (time) [2011-08-12T15:50:16.784000000, 2011-08-12T17:22:05.000000000]\n", "Data:\n", " float64 [mm] (time) [-10, -10]\n", "\n", "]\n", " sample PyObject [dimensionless] () []\n", " samplerequest DataArray [dimensionless] () [\n", "Dimensions: Sizes[time:2, ]\n", "Coordinates:\n", " time datetime64 [ns] (time) [2011-08-12T15:50:16.784000000, 2011-08-12T17:22:05.000000000]\n", "Data:\n", " float64 [dimensionless] (time) [4, 4]\n", "\n", "]\n", " start_time string [dimensionless] () [\"2011-08-12T15:50:17\"]\n", " vGuide DataArray [dimensionless] () [\n", "Dimensions: Sizes[time:2, ]\n", "Coordinates:\n", " time datetime64 [ns] (time) [2011-08-12T15:50:16.800000000, 2011-08-12T17:22:05.000000000]\n", "Data:\n", " float64 [dimensionless] (time) [2, 2]\n", "\n", "]\n", " veto_pulse_time DataArray [dimensionless] () [\n", "Dimensions: Sizes[time:1, ]\n", "Coordinates:\n", " time datetime64 [ns] (time) [2011-08-12T15:50:17.000000000]\n", "Data:\n", " float64 [dimensionless] (time) [0]\n", "\n", "]\n" ] }, "execution_count": 3, "metadata": {}, "output_type": "execute_result" } ], "source": [ "sample" ] }, { "cell_type": "code", "execution_count": 4, "metadata": {}, "outputs": [ { "data": { "application/vnd.jupyter.widget-view+json": { "model_id": "b338818e8b2c4c28be816ad5a1ca6dae", "version_major": 2, "version_minor": 0 }, "text/plain": [ "VBox(children=(HBox(children=(VBox(children=(Button(icon='home', layout=Layout(padding='0px 0px 0px 0px', widt…" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "sc.plot(sample)" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "### Instrument view\n", "\n", "Scipp provides a simple 3D instrument view inspired by Mantid's own [instrument view](https://www.mantidproject.org/MantidPlot:_Instrument_View), which can be used to take a quick look at the neutron counts on the detector panels in 3D space or using various cylindrical and spherical projections" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "
\n", "\n", "**Note**\n", "\n", "The support and handling of event data is currently being rewritten.\n", "The interface for manipulating event data (binned) data is not representative for the final API.\n", "\n", "
" ] }, { "cell_type": "code", "execution_count": 5, "metadata": {}, "outputs": [ { "data": { "application/vnd.jupyter.widget-view+json": { "model_id": "cf78004cbcf64995b09696b24147a069", "version_major": 2, "version_minor": 0 }, "text/plain": [ "VBox(children=(HBox(children=(VBox(children=(Button(icon='home', layout=Layout(padding='0px 0px 0px 0px', widt…" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "scn.instrument_view(sample)" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "### Plot against scattering angle $\\theta$ using `groupby`\n", "\n", "*This is not an essential step and can be skipped.*\n", "\n", "Plotting raw data directly yields a hard-to-interpret figure.\n", "We can obtain something more useful by \"binning\" the spectrum axis based on its $\\theta$ value, using the split-apply-combine approach provided by `groupby`:" ] }, { "cell_type": "code", "execution_count": 6, "metadata": {}, "outputs": [], "source": [ "sample.coords['two_theta'] = scn.two_theta(sample)\n", "vanadium.coords['two_theta'] = scn.two_theta(vanadium)\n", "two_theta_bins = sc.Variable(dims=['two_theta'],\n", " unit=sc.units.rad,\n", " values=np.linspace(0.0, np.pi, num=2000))" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "We concatenate events lists from different spectra that fall into a given $2\\theta$ range into longer combined lists:" ] }, { "cell_type": "code", "execution_count": 7, "metadata": {}, "outputs": [ { "ename": "AttributeError", "evalue": "'GroupbyBins' object has no attribute 'concatenate'", "output_type": "error", "traceback": [ "\u001b[0;31m---------------------------------------------------------------------------\u001b[0m", "\u001b[0;31mAttributeError\u001b[0m Traceback (most recent call last)", "\u001b[0;32m/tmp/ipykernel_25962/17277667.py\u001b[0m in \u001b[0;36m\u001b[0;34m\u001b[0m\n\u001b[0;32m----> 1\u001b[0;31m \u001b[0mtheta_sample\u001b[0m \u001b[0;34m=\u001b[0m \u001b[0msc\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mgroupby\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0msample\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0;34m'two_theta'\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0mbins\u001b[0m\u001b[0;34m=\u001b[0m\u001b[0mtwo_theta_bins\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mbins\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mconcatenate\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0;34m'spectrum'\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0m", "\u001b[0;31mAttributeError\u001b[0m: 'GroupbyBins' object has no attribute 'concatenate'" ] } ], "source": [ "theta_sample = sc.groupby(sample, 'two_theta', bins=two_theta_bins).bins.concatenate('spectrum')" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "
\n", "\n", "**Note**\n", " \n", "Use `groupby.sum` instead of `groupby.bins.concatenate` when working with dense (histogrammed) data\n", "
" ] }, { "cell_type": "code", "execution_count": 8, "metadata": {}, "outputs": [ { "ename": "NameError", "evalue": "name 'theta_sample' is not defined", "output_type": "error", "traceback": [ "\u001b[0;31m---------------------------------------------------------------------------\u001b[0m", "\u001b[0;31mNameError\u001b[0m Traceback (most recent call last)", "\u001b[0;32m/tmp/ipykernel_25962/3881907689.py\u001b[0m in \u001b[0;36m\u001b[0;34m\u001b[0m\n\u001b[0;32m----> 1\u001b[0;31m \u001b[0mtheta_sample\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0m", "\u001b[0;31mNameError\u001b[0m: name 'theta_sample' is not defined" ] } ], "source": [ "theta_sample" ] }, { "cell_type": "code", "execution_count": 9, "metadata": {}, "outputs": [ { "ename": "NameError", "evalue": "name 'theta_sample' is not defined", "output_type": "error", "traceback": [ "\u001b[0;31m---------------------------------------------------------------------------\u001b[0m", "\u001b[0;31mNameError\u001b[0m Traceback (most recent call last)", "\u001b[0;32m/tmp/ipykernel_25962/2836602988.py\u001b[0m in \u001b[0;36m\u001b[0;34m\u001b[0m\n\u001b[0;32m----> 1\u001b[0;31m \u001b[0msc\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mplot\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0mtheta_sample\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0m", "\u001b[0;31mNameError\u001b[0m: name 'theta_sample' is not defined" ] } ], "source": [ "sc.plot(theta_sample)" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "### Unit conversion\n", "\n", "*Note: We are back to working with `sample`, not `theta_sample`.*\n", "\n", "`scippneutron` provides means to convert between units (dimensions) related to time-of-flight.\n", "The loaded raw data has dimension `tof`, and we convert to interplanar lattice spacing (dspacing):" ] }, { "cell_type": "code", "execution_count": 10, "metadata": {}, "outputs": [ { "data": { "text/html": [ "
\n", "\n", "\n", "\n", "\n", "Show/Hide data repr\n", "\n", "\n", "\n", "\n", "\n", "Show/Hide attributes\n", "\n", "\n", "\n", "\n", "\n", "\n", "
scipp.DataArray (412.87 MB)
    • spectrum: 24794
    • dspacing: 1
    • dspacing
      (spectrum, dspacing [bin-edge])
      float64
      Å
      0.001, 0.557, ..., 0.003, 2.343
      Values:
      array([[6.33939029e-04, 5.57008885e-01],\n",
       "       [6.32685444e-04, 5.55907426e-01],\n",
       "       [6.31487177e-04, 5.54854572e-01],\n",
       "       ...,\n",
       "       [2.67322756e-03, 2.34882447e+00],\n",
       "       [2.67090407e-03, 2.34678294e+00],\n",
       "       [2.66651403e-03, 2.34292565e+00]])
    • spectrum
      (spectrum)
      int32
      1, 2, ..., 24793, 24794
      Values:
      array([    1,     2,     3, ..., 24792, 24793, 24794], dtype=int32)
    • (spectrum, dspacing)
      DataArrayView
      binned data [len=0, len=0, ..., len=0, len=0]
      Values:
      [<scipp.DataArray>\n",
       "Dimensions: Sizes[event:0, ]\n",
       "Coordinates:\n",
       "  dspacing                  float64             [Å]  (event)  []\n",
       "Data:\n",
       "                            float32         [counts]  (event)  []  []\n",
       "Attributes:\n",
       "  tof                       float64            [µs]  (event)  []\n",
       "\n",
       ", <scipp.DataArray>\n",
       "Dimensions: Sizes[event:0, ]\n",
       "Coordinates:\n",
       "  dspacing                  float64             [Å]  (event)  []\n",
       "Data:\n",
       "                            float32         [counts]  (event)  []  []\n",
       "Attributes:\n",
       "  tof                       float64            [µs]  (event)  []\n",
       "\n",
       ", ..., <scipp.DataArray>\n",
       "Dimensions: Sizes[event:0, ]\n",
       "Coordinates:\n",
       "  dspacing                  float64             [Å]  (event)  []\n",
       "Data:\n",
       "                            float32         [counts]  (event)  []  []\n",
       "Attributes:\n",
       "  tof                       float64            [µs]  (event)  []\n",
       "\n",
       ", <scipp.DataArray>\n",
       "Dimensions: Sizes[event:0, ]\n",
       "Coordinates:\n",
       "  dspacing                  float64             [Å]  (event)  []\n",
       "Data:\n",
       "                            float32         [counts]  (event)  []  []\n",
       "Attributes:\n",
       "  tof                       float64            [µs]  (event)  []\n",
       "\n",
       "]
    • ChopperStatus1
      ()
      DataArray
      {dims=[time: 2], coords=[time]}
      Values:
      <scipp.DataArray>\n",
       "Dimensions: Sizes[time:2, ]\n",
       "Coordinates:\n",
       "  time                    datetime64             [ns]  (time)  [2011-08-12T15:50:16.800000000, 2011-08-12T17:22:05.000000000]\n",
       "Data:\n",
       "                            float64  [dimensionless]  (time)  [4, 4]\n",
       "\n",
       "
    • ChopperStatus2
      ()
      DataArray
      {dims=[time: 2], coords=[time]}
      Values:
      <scipp.DataArray>\n",
       "Dimensions: Sizes[time:2, ]\n",
       "Coordinates:\n",
       "  time                    datetime64             [ns]  (time)  [2011-08-12T15:50:16.816000000, 2011-08-12T17:22:05.000000000]\n",
       "Data:\n",
       "                            float64  [dimensionless]  (time)  [4, 4]\n",
       "\n",
       "
    • ChopperStatus3
      ()
      DataArray
      {dims=[time: 2], coords=[time]}
      Values:
      <scipp.DataArray>\n",
       "Dimensions: Sizes[time:2, ]\n",
       "Coordinates:\n",
       "  time                    datetime64             [ns]  (time)  [2011-08-12T15:50:16.816000000, 2011-08-12T17:22:05.000000000]\n",
       "Data:\n",
       "                            float64  [dimensionless]  (time)  [4, 4]\n",
       "\n",
       "
    • CurrentSP
      ()
      DataArray
      {dims=[time: 2], coords=[time]}
      Values:
      <scipp.DataArray>\n",
       "Dimensions: Sizes[time:2, ]\n",
       "Coordinates:\n",
       "  time                    datetime64             [ns]  (time)  [2011-08-12T15:50:16.800000000, 2011-08-12T17:22:05.000000000]\n",
       "Data:\n",
       "                            float64  [dimensionless]  (time)  [300, 300]\n",
       "\n",
       "
    • EnergyRequest
      ()
      DataArray
      {dims=[time: 2], unit=meV, coords=[time]}
      Values:
      <scipp.DataArray>\n",
       "Dimensions: Sizes[time:2, ]\n",
       "Coordinates:\n",
       "  time                    datetime64             [ns]  (time)  [2011-08-12T15:50:16.784000000, 2011-08-12T17:22:05.000000000]\n",
       "Data:\n",
       "                            float64            [meV]  (time)  [287.955, 287.955]\n",
       "\n",
       "
    • Filename
      ()
      string
      /home/scientist/.cache/scippneutron/2/PG3_4844_event.nxs
      Values:
      '/home/scientist/.cache/scippneutron/2/PG3_4844_event.nxs'
    • L1
      ()
      float64
      m
      60.0
      Values:
      array(60.)
    • L2
      (spectrum)
      float64
      m
      2.560, 2.556, ..., 4.251, 4.252
      Values:
      array([2.56027902, 2.55590928, 2.55268677, ..., 4.25113991, 4.25116338,\n",
       "       4.25188   ])
    • LKSRampRate
      ()
      DataArray
      {dims=[time: 2], coords=[time]}
      Values:
      <scipp.DataArray>\n",
       "Dimensions: Sizes[time:2, ]\n",
       "Coordinates:\n",
       "  time                    datetime64             [ns]  (time)  [2011-08-12T15:50:16.800000000, 2011-08-12T17:22:05.000000000]\n",
       "Data:\n",
       "                            float64  [dimensionless]  (time)  [0, 0]\n",
       "\n",
       "
    • LambdaRequest
      ()
      DataArray
      {dims=[time: 2], unit=Å, coords=[time]}
      Values:
      <scipp.DataArray>\n",
       "Dimensions: Sizes[time:2, ]\n",
       "Coordinates:\n",
       "  time                    datetime64             [ns]  (time)  [2011-08-12T15:50:16.784000000, 2011-08-12T17:22:05.000000000]\n",
       "Data:\n",
       "                            float64             [Å]  (time)  [0.533, 0.533]\n",
       "\n",
       "
    • Ltotal
      (spectrum)
      float64
      m
      62.560, 62.556, ..., 64.251, 64.252
      Values:
      array([62.56027902, 62.55590928, 62.55268677, ..., 64.25113991,\n",
       "       64.25116338, 64.25188   ])
    • Phase1
      ()
      DataArray
      {dims=[time: 1794], unit=µs, coords=[time]}
      Values:
      <scipp.DataArray>\n",
       "Dimensions: Sizes[time:1794, ]\n",
       "Coordinates:\n",
       "  time                    datetime64             [ns]  (time)  [2011-08-12T15:50:16.800000000, 2011-08-12T15:50:18.267000007, ..., 2011-08-12T17:22:02.999218750, 2011-08-12T17:22:05.000000000]\n",
       "Data:\n",
       "                            float64            [µs]  (time)  [8166.72, 8165.16, ..., 8163.85, 8163.85]\n",
       "\n",
       "
    • Phase2
      ()
      DataArray
      {dims=[time: 1793], unit=µs, coords=[time]}
      Values:
      <scipp.DataArray>\n",
       "Dimensions: Sizes[time:1793, ]\n",
       "Coordinates:\n",
       "  time                    datetime64             [ns]  (time)  [2011-08-12T15:50:16.800000000, 2011-08-12T15:50:18.267000007, ..., 2011-08-12T17:22:02.999218750, 2011-08-12T17:22:05.000000000]\n",
       "Data:\n",
       "                            float64            [µs]  (time)  [8335.63, 8334.09, ..., 8332.86, 8332.86]\n",
       "\n",
       "
    • Phase3
      ()
      DataArray
      {dims=[time: 1777], unit=µs, coords=[time]}
      Values:
      <scipp.DataArray>\n",
       "Dimensions: Sizes[time:1777, ]\n",
       "Coordinates:\n",
       "  time                    datetime64             [ns]  (time)  [2011-08-12T15:50:16.800000000, 2011-08-12T15:50:19.688999938, ..., 2011-08-12T17:22:04.421093750, 2011-08-12T17:22:05.000000000]\n",
       "Data:\n",
       "                            float64            [µs]  (time)  [14001.5, 14001.4, ..., 14002, 14002]\n",
       "\n",
       "
    • PhaseRequest1
      ()
      DataArray
      {dims=[time: 2], unit=µs, coords=[time]}
      Values:
      <scipp.DataArray>\n",
       "Dimensions: Sizes[time:2, ]\n",
       "Coordinates:\n",
       "  time                    datetime64             [ns]  (time)  [2011-08-12T15:50:16.800000000, 2011-08-12T17:22:05.000000000]\n",
       "Data:\n",
       "                            float64            [µs]  (time)  [8164.07, 8164.07]\n",
       "\n",
       "
    • PhaseRequest2
      ()
      DataArray
      {dims=[time: 2], unit=µs, coords=[time]}
      Values:
      <scipp.DataArray>\n",
       "Dimensions: Sizes[time:2, ]\n",
       "Coordinates:\n",
       "  time                    datetime64             [ns]  (time)  [2011-08-12T15:50:16.800000000, 2011-08-12T17:22:05.000000000]\n",
       "Data:\n",
       "                            float64            [µs]  (time)  [8332.89, 8332.89]\n",
       "\n",
       "
    • PhaseRequest3
      ()
      DataArray
      {dims=[time: 2], unit=µs, coords=[time]}
      Values:
      <scipp.DataArray>\n",
       "Dimensions: Sizes[time:2, ]\n",
       "Coordinates:\n",
       "  time                    datetime64             [ns]  (time)  [2011-08-12T15:50:16.800000000, 2011-08-12T17:22:05.000000000]\n",
       "Data:\n",
       "                            float64            [µs]  (time)  [14001.8, 14001.8]\n",
       "\n",
       "
    • S1HCenter
      ()
      DataArray
      {dims=[time: 2], unit=mm, coords=[time]}
      Values:
      <scipp.DataArray>\n",
       "Dimensions: Sizes[time:2, ]\n",
       "Coordinates:\n",
       "  time                    datetime64             [ns]  (time)  [2011-08-12T15:50:16.800000000, 2011-08-12T17:22:05.000000000]\n",
       "Data:\n",
       "                            float64             [mm]  (time)  [0, 0]\n",
       "\n",
       "
    • S1HCenterOffset
      ()
      DataArray
      {dims=[time: 2], coords=[time]}
      Values:
      <scipp.DataArray>\n",
       "Dimensions: Sizes[time:2, ]\n",
       "Coordinates:\n",
       "  time                    datetime64             [ns]  (time)  [2011-08-12T15:50:16.800000000, 2011-08-12T17:22:05.000000000]\n",
       "Data:\n",
       "                            float64  [dimensionless]  (time)  [0, 0]\n",
       "\n",
       "
    • S1HWidth
      ()
      DataArray
      {dims=[time: 2], unit=mm, coords=[time]}
      Values:
      <scipp.DataArray>\n",
       "Dimensions: Sizes[time:2, ]\n",
       "Coordinates:\n",
       "  time                    datetime64             [ns]  (time)  [2011-08-12T15:50:16.800000000, 2011-08-12T17:22:05.000000000]\n",
       "Data:\n",
       "                            float64             [mm]  (time)  [10, 10]\n",
       "\n",
       "
    • S1VCenter
      ()
      DataArray
      {dims=[time: 2], unit=mm, coords=[time]}
      Values:
      <scipp.DataArray>\n",
       "Dimensions: Sizes[time:2, ]\n",
       "Coordinates:\n",
       "  time                    datetime64             [ns]  (time)  [2011-08-12T15:50:16.800000000, 2011-08-12T17:22:05.000000000]\n",
       "Data:\n",
       "                            float64             [mm]  (time)  [5, 5]\n",
       "\n",
       "
    • S1VCenterOffset
      ()
      DataArray
      {dims=[time: 2], coords=[time]}
      Values:
      <scipp.DataArray>\n",
       "Dimensions: Sizes[time:2, ]\n",
       "Coordinates:\n",
       "  time                    datetime64             [ns]  (time)  [2011-08-12T15:50:16.800000000, 2011-08-12T17:22:05.000000000]\n",
       "Data:\n",
       "                            float64  [dimensionless]  (time)  [0, 0]\n",
       "\n",
       "
    • S1VHeight
      ()
      DataArray
      {dims=[time: 2], unit=mm, coords=[time]}
      Values:
      <scipp.DataArray>\n",
       "Dimensions: Sizes[time:2, ]\n",
       "Coordinates:\n",
       "  time                    datetime64             [ns]  (time)  [2011-08-12T15:50:16.816000000, 2011-08-12T17:22:05.000000000]\n",
       "Data:\n",
       "                            float64             [mm]  (time)  [30, 30]\n",
       "\n",
       "
    • SampleTemp
      ()
      DataArray
      {dims=[time: 467], coords=[time]}
      Values:
      <scipp.DataArray>\n",
       "Dimensions: Sizes[time:467, ]\n",
       "Coordinates:\n",
       "  time                    datetime64             [ns]  (time)  [2011-08-12T15:50:16.784000000, 2011-08-12T15:50:18.001000007, ..., 2011-08-12T17:21:58.623843750, 2011-08-12T17:22:05.000000000]\n",
       "Data:\n",
       "                            float64  [dimensionless]  (time)  [299.352, 299.446, ..., 300, 300]\n",
       "\n",
       "
    • Speed1
      ()
      DataArray
      {dims=[time: 2], unit=Hz, coords=[time]}
      Values:
      <scipp.DataArray>\n",
       "Dimensions: Sizes[time:2, ]\n",
       "Coordinates:\n",
       "  time                    datetime64             [ns]  (time)  [2011-08-12T15:50:16.800000000, 2011-08-12T17:22:05.000000000]\n",
       "Data:\n",
       "                            float64             [Hz]  (time)  [60, 60]\n",
       "\n",
       "
    • Speed2
      ()
      DataArray
      {dims=[time: 2], unit=Hz, coords=[time]}
      Values:
      <scipp.DataArray>\n",
       "Dimensions: Sizes[time:2, ]\n",
       "Coordinates:\n",
       "  time                    datetime64             [ns]  (time)  [2011-08-12T15:50:16.800000000, 2011-08-12T17:22:05.000000000]\n",
       "Data:\n",
       "                            float64             [Hz]  (time)  [60, 60]\n",
       "\n",
       "
    • Speed3
      ()
      DataArray
      {dims=[time: 2], unit=Hz, coords=[time]}
      Values:
      <scipp.DataArray>\n",
       "Dimensions: Sizes[time:2, ]\n",
       "Coordinates:\n",
       "  time                    datetime64             [ns]  (time)  [2011-08-12T15:50:16.800000000, 2011-08-12T17:22:05.000000000]\n",
       "Data:\n",
       "                            float64             [Hz]  (time)  [60, 60]\n",
       "\n",
       "
    • SpeedRequest1
      ()
      DataArray
      {dims=[time: 2], unit=Hz, coords=[time]}
      Values:
      <scipp.DataArray>\n",
       "Dimensions: Sizes[time:2, ]\n",
       "Coordinates:\n",
       "  time                    datetime64             [ns]  (time)  [2011-08-12T15:50:16.784000000, 2011-08-12T17:22:05.000000000]\n",
       "Data:\n",
       "                            float64             [Hz]  (time)  [60, 60]\n",
       "\n",
       "
    • SpeedRequest2
      ()
      DataArray
      {dims=[time: 2], unit=Hz, coords=[time]}
      Values:
      <scipp.DataArray>\n",
       "Dimensions: Sizes[time:2, ]\n",
       "Coordinates:\n",
       "  time                    datetime64             [ns]  (time)  [2011-08-12T15:50:16.784000000, 2011-08-12T17:22:05.000000000]\n",
       "Data:\n",
       "                            float64             [Hz]  (time)  [60, 60]\n",
       "\n",
       "
    • SpeedRequest3
      ()
      DataArray
      {dims=[time: 2], unit=Hz, coords=[time]}
      Values:
      <scipp.DataArray>\n",
       "Dimensions: Sizes[time:2, ]\n",
       "Coordinates:\n",
       "  time                    datetime64             [ns]  (time)  [2011-08-12T15:50:16.784000000, 2011-08-12T17:22:05.000000000]\n",
       "Data:\n",
       "                            float64             [Hz]  (time)  [60, 60]\n",
       "\n",
       "
    • TolRequest
      ()
      DataArray
      {dims=[time: 2], coords=[time]}
      Values:
      <scipp.DataArray>\n",
       "Dimensions: Sizes[time:2, ]\n",
       "Coordinates:\n",
       "  time                    datetime64             [ns]  (time)  [2011-08-12T15:50:16.800000000, 2011-08-12T17:22:05.000000000]\n",
       "Data:\n",
       "                            float64  [dimensionless]  (time)  [20, 20]\n",
       "\n",
       "
    • commErrs
      ()
      DataArray
      {dims=[time: 2], coords=[time]}
      Values:
      <scipp.DataArray>\n",
       "Dimensions: Sizes[time:2, ]\n",
       "Coordinates:\n",
       "  time                    datetime64             [ns]  (time)  [2011-08-12T15:50:16.784000000, 2011-08-12T17:22:05.000000000]\n",
       "Data:\n",
       "                            float64  [dimensionless]  (time)  [0, 0]\n",
       "\n",
       "
    • currentsample
      ()
      DataArray
      {dims=[time: 2], coords=[time]}
      Values:
      <scipp.DataArray>\n",
       "Dimensions: Sizes[time:2, ]\n",
       "Coordinates:\n",
       "  time                    datetime64             [ns]  (time)  [2011-08-12T15:50:16.784000000, 2011-08-12T17:22:05.000000000]\n",
       "Data:\n",
       "                            float64  [dimensionless]  (time)  [4, 4]\n",
       "\n",
       "
    • duration
      ()
      float64
      s
      5508.0
      Values:
      array(5508.)
    • end_time
      ()
      string
      2011-08-12T17:22:05
      Values:
      '2011-08-12T17:22:05'
    • experiment_identifier
      ()
      string
      IPTS-2767
      Values:
      'IPTS-2767'
    • fernsstatus
      ()
      DataArray
      {dims=[time: 2], coords=[time]}
      Values:
      <scipp.DataArray>\n",
       "Dimensions: Sizes[time:2, ]\n",
       "Coordinates:\n",
       "  time                    datetime64             [ns]  (time)  [2011-08-12T15:50:16.784000000, 2011-08-12T17:22:05.000000000]\n",
       "Data:\n",
       "                            float64  [dimensionless]  (time)  [3, 3]\n",
       "\n",
       "
    • file_notes
      ()
      string
      NONE
      Values:
      'NONE'
    • frequency
      ()
      DataArray
      {dims=[time: 330473], unit=Hz, coords=[time]}
      Values:
      <scipp.DataArray>\n",
       "Dimensions: Sizes[time:330473, ]\n",
       "Coordinates:\n",
       "  time                    datetime64             [ns]  (time)  [2011-08-12T15:50:17.000000000, 2011-08-12T15:50:17.016659999, ..., 2011-08-12T17:22:05.085449218, 2011-08-12T17:22:05.102050781]\n",
       "Data:\n",
       "                            float64             [Hz]  (time)  [0, 60.024, ..., 60.0024, 59.9988]\n",
       "\n",
       "
    • gd_prtn_chrg
      ()
      float64
      µAh
      1171.953902925
      Values:
      array(1171.95390292)
    • guide
      ()
      DataArray
      {dims=[time: 2], unit=mm, coords=[time]}
      Values:
      <scipp.DataArray>\n",
       "Dimensions: Sizes[time:2, ]\n",
       "Coordinates:\n",
       "  time                    datetime64             [ns]  (time)  [2011-08-12T15:50:16.784000000, 2011-08-12T17:22:05.000000000]\n",
       "Data:\n",
       "                            float64             [mm]  (time)  [-55.463, -55.463]\n",
       "\n",
       "
    • incident_beam
      ()
      vector3
      m
      [ 0. 0. 60.]
      Values:
      array([ 0.,  0., 60.])
    • instrument_name
      ()
      string
      POWGEN
      Values:
      'POWGEN'
    • monitor1
      ()
      DataArray
      {dims=[tof: 200001], unit=counts, coords=[tof, position, source_position]}
      Values:
      <scipp.DataArray>\n",
       "Dimensions: Sizes[tof:200001, ]\n",
       "Coordinates:\n",
       "  position                  vector3              [m]  ()  [(0, 0, -1)]\n",
       "  source_position           vector3              [m]  ()  [(0, 0, -60)]\n",
       "  tof                       float64            [µs]  (tof [bin-edge])  [0, 1, ..., 200000, 200001]\n",
       "Data:\n",
       "                            float64         [counts]  (tof)  [25, 10, ..., 0, 0]  [25, 10, ..., 0, 0]\n",
       "Attributes:\n",
       "  instrument_name            string  [dimensionless]  ()  ["POWGEN"]\n",
       "  sample_position           vector3              [m]  ()  [(0, 0, 0)]\n",
       "  spectrum                    int32  [dimensionless]  ()  [1]\n",
       "\n",
       "
    • position
      (spectrum)
      vector3
      m
      [ 1.17451004 -1.01106149 -2.03796699], [ 1.18147634 -0.95946649 -2.05334117], ..., [1.81428985 0.09565841 3.84338287], [1.81375055 0.1499371 3.84269584]
      Values:
      array([[ 1.17451004, -1.01106149, -2.03796699],\n",
       "       [ 1.18147634, -0.95946649, -2.05334117],\n",
       "       [ 1.18844265, -0.90787149, -2.06871534],\n",
       "       ...,\n",
       "       [ 1.81482915,  0.04137972,  3.8440699 ],\n",
       "       [ 1.81428985,  0.09565841,  3.84338287],\n",
       "       [ 1.81375055,  0.1499371 ,  3.84269584]])
    • proton_charge
      ()
      DataArray
      {dims=[time: 330473], unit=pC, coords=[time]}
      Values:
      <scipp.DataArray>\n",
       "Dimensions: Sizes[time:330473, ]\n",
       "Coordinates:\n",
       "  time                    datetime64             [ns]  (time)  [2011-08-12T15:50:17.000000000, 2011-08-12T15:50:17.016659999, ..., 2011-08-12T17:22:05.085449218, 2011-08-12T17:22:05.102050781]\n",
       "Data:\n",
       "                            float64             [pC]  (time)  [1.48434e+07, 1.48434e+07, ..., 1.48733e+07, 1.48391e+07]\n",
       "\n",
       "
    • run_number
      ()
      string
      4844
      Values:
      '4844'
    • run_start
      ()
      string
      2011-08-12T15:50:17
      Values:
      '2011-08-12T15:50:17'
    • run_title
      ()
      string
      diamond cw0.533 4.22e12 60Hz [10x30]
      Values:
      'diamond cw0.533 4.22e12 60Hz [10x30]'
    • running
      ()
      DataArray
      {dims=[time: 1], coords=[time]}
      Values:
      <scipp.DataArray>\n",
       "Dimensions: Sizes[time:1, ]\n",
       "Coordinates:\n",
       "  time                    datetime64             [ns]  (time)  [1990-01-01T00:00:00.000000000]\n",
       "Data:\n",
       "                               bool  [dimensionless]  (time)  [True]\n",
       "\n",
       "
    • s1b
      ()
      DataArray
      {dims=[time: 2], unit=mm, coords=[time]}
      Values:
      <scipp.DataArray>\n",
       "Dimensions: Sizes[time:2, ]\n",
       "Coordinates:\n",
       "  time                    datetime64             [ns]  (time)  [2011-08-12T15:50:16.784000000, 2011-08-12T17:22:05.000000000]\n",
       "Data:\n",
       "                            float64             [mm]  (time)  [20, 20]\n",
       "\n",
       "
    • s1l
      ()
      DataArray
      {dims=[time: 2], unit=mm, coords=[time]}
      Values:
      <scipp.DataArray>\n",
       "Dimensions: Sizes[time:2, ]\n",
       "Coordinates:\n",
       "  time                    datetime64             [ns]  (time)  [2011-08-12T15:50:16.800000000, 2011-08-12T17:22:05.000000000]\n",
       "Data:\n",
       "                            float64             [mm]  (time)  [5, 5]\n",
       "\n",
       "
    • s1r
      ()
      DataArray
      {dims=[time: 2], unit=mm, coords=[time]}
      Values:
      <scipp.DataArray>\n",
       "Dimensions: Sizes[time:2, ]\n",
       "Coordinates:\n",
       "  time                    datetime64             [ns]  (time)  [2011-08-12T15:50:16.784000000, 2011-08-12T17:22:05.000000000]\n",
       "Data:\n",
       "                            float64             [mm]  (time)  [-5, -5]\n",
       "\n",
       "
    • s1t
      ()
      DataArray
      {dims=[time: 2], unit=mm, coords=[time]}
      Values:
      <scipp.DataArray>\n",
       "Dimensions: Sizes[time:2, ]\n",
       "Coordinates:\n",
       "  time                    datetime64             [ns]  (time)  [2011-08-12T15:50:16.784000000, 2011-08-12T17:22:05.000000000]\n",
       "Data:\n",
       "                            float64             [mm]  (time)  [-10, -10]\n",
       "\n",
       "
    • sample
      ()
      PyObject
      <mantid.api._api.Sample object at 0x7f8bea7b65f0>
      Values:
      <mantid.api._api.Sample object at 0x7f8bea7b65f0>
    • sample_position
      ()
      vector3
      m
      [0. 0. 0.]
      Values:
      array([0., 0., 0.])
    • samplerequest
      ()
      DataArray
      {dims=[time: 2], coords=[time]}
      Values:
      <scipp.DataArray>\n",
       "Dimensions: Sizes[time:2, ]\n",
       "Coordinates:\n",
       "  time                    datetime64             [ns]  (time)  [2011-08-12T15:50:16.784000000, 2011-08-12T17:22:05.000000000]\n",
       "Data:\n",
       "                            float64  [dimensionless]  (time)  [4, 4]\n",
       "\n",
       "
    • scattered_beam
      (spectrum)
      vector3
      m
      [ 1.17451004 -1.01106149 -2.03796699], [ 1.18147634 -0.95946649 -2.05334117], ..., [1.81428985 0.09565841 3.84338287], [1.81375055 0.1499371 3.84269584]
      Values:
      array([[ 1.17451004, -1.01106149, -2.03796699],\n",
       "       [ 1.18147634, -0.95946649, -2.05334117],\n",
       "       [ 1.18844265, -0.90787149, -2.06871534],\n",
       "       ...,\n",
       "       [ 1.81482915,  0.04137972,  3.8440699 ],\n",
       "       [ 1.81428985,  0.09565841,  3.84338287],\n",
       "       [ 1.81375055,  0.1499371 ,  3.84269584]])
    • source_position
      ()
      vector3
      m
      [ 0. 0. -60.]
      Values:
      array([  0.,   0., -60.])
    • start_time
      ()
      string
      2011-08-12T15:50:17
      Values:
      '2011-08-12T15:50:17'
    • tof
      (dspacing [bin-edge])
      float64
      µs
      19.0, 1.669e+04
      Values:
      array([   19.        , 16694.30078125])
    • two_theta
      (spectrum)
      float64
      rad
      2.491, 2.504, ..., 0.442, 0.442
      Values:
      array([2.49144445, 2.50372965, 2.51564288, ..., 0.44118918, 0.44157918,\n",
       "       0.44231324])
    • vGuide
      ()
      DataArray
      {dims=[time: 2], coords=[time]}
      Values:
      <scipp.DataArray>\n",
       "Dimensions: Sizes[time:2, ]\n",
       "Coordinates:\n",
       "  time                    datetime64             [ns]  (time)  [2011-08-12T15:50:16.800000000, 2011-08-12T17:22:05.000000000]\n",
       "Data:\n",
       "                            float64  [dimensionless]  (time)  [2, 2]\n",
       "\n",
       "
    • veto_pulse_time
      ()
      DataArray
      {dims=[time: 1], coords=[time]}
      Values:
      <scipp.DataArray>\n",
       "Dimensions: Sizes[time:1, ]\n",
       "Coordinates:\n",
       "  time                    datetime64             [ns]  (time)  [2011-08-12T15:50:17.000000000]\n",
       "Data:\n",
       "                            float64  [dimensionless]  (time)  [0]\n",
       "\n",
       "
" ], "text/plain": [ "\n", "Dimensions: Sizes[spectrum:24794, dspacing:1, ]\n", "Coordinates:\n", " dspacing float64 [Å] (spectrum, dspacing [bin-edge]) [0.000633939, 0.557009, ..., 0.00266651, 2.34293]\n", " spectrum int32 [dimensionless] (spectrum) [1, 2, ..., 24793, 24794]\n", "Data:\n", " DataArrayView [dimensionless] (spectrum, dspacing) [\n", "Dimensions: Sizes[event:0, ]\n", "Coordinates:\n", " dspacing float64 [Å] (event) []\n", "Data:\n", " float32 [counts] (event) [] []\n", "Attributes:\n", " tof float64 [µs] (event) []\n", "\n", ", \n", "Dimensions: Sizes[event:0, ]\n", "Coordinates:\n", " dspacing float64 [Å] (event) []\n", "Data:\n", " float32 [counts] (event) [] []\n", "Attributes:\n", " tof float64 [µs] (event) []\n", "\n", ", ..., \n", "Dimensions: Sizes[event:0, ]\n", "Coordinates:\n", " dspacing float64 [Å] (event) []\n", "Data:\n", " float32 [counts] (event) [] []\n", "Attributes:\n", " tof float64 [µs] (event) []\n", "\n", ", \n", "Dimensions: Sizes[event:0, ]\n", "Coordinates:\n", " dspacing float64 [Å] (event) []\n", "Data:\n", " float32 [counts] (event) [] []\n", "Attributes:\n", " tof float64 [µs] (event) []\n", "\n", "]\n", "Attributes:\n", " ChopperStatus1 DataArray [dimensionless] () [\n", "Dimensions: Sizes[time:2, ]\n", "Coordinates:\n", " time datetime64 [ns] (time) [2011-08-12T15:50:16.800000000, 2011-08-12T17:22:05.000000000]\n", "Data:\n", " float64 [dimensionless] (time) [4, 4]\n", "\n", "]\n", " ChopperStatus2 DataArray [dimensionless] () [\n", "Dimensions: Sizes[time:2, ]\n", "Coordinates:\n", " time datetime64 [ns] (time) [2011-08-12T15:50:16.816000000, 2011-08-12T17:22:05.000000000]\n", "Data:\n", " float64 [dimensionless] (time) [4, 4]\n", "\n", "]\n", " ChopperStatus3 DataArray [dimensionless] () [\n", "Dimensions: Sizes[time:2, ]\n", "Coordinates:\n", " time datetime64 [ns] (time) [2011-08-12T15:50:16.816000000, 2011-08-12T17:22:05.000000000]\n", "Data:\n", " float64 [dimensionless] (time) [4, 4]\n", "\n", "]\n", " CurrentSP DataArray [dimensionless] () [\n", "Dimensions: Sizes[time:2, ]\n", "Coordinates:\n", " time datetime64 [ns] (time) [2011-08-12T15:50:16.800000000, 2011-08-12T17:22:05.000000000]\n", "Data:\n", " float64 [dimensionless] (time) [300, 300]\n", "\n", "]\n", " EnergyRequest DataArray [dimensionless] () [\n", "Dimensions: Sizes[time:2, ]\n", "Coordinates:\n", " time datetime64 [ns] (time) [2011-08-12T15:50:16.784000000, 2011-08-12T17:22:05.000000000]\n", "Data:\n", " float64 [meV] (time) [287.955, 287.955]\n", "\n", "]\n", " Filename string [dimensionless] () [\"/home/scientist/.cache/scippneutron/2/PG3_4844_event.nxs\"]\n", " L1 float64 [m] () [60]\n", " L2 float64 [m] (spectrum) [2.56028, 2.55591, ..., 4.25116, 4.25188]\n", " LKSRampRate DataArray [dimensionless] () [\n", "Dimensions: Sizes[time:2, ]\n", "Coordinates:\n", " time datetime64 [ns] (time) [2011-08-12T15:50:16.800000000, 2011-08-12T17:22:05.000000000]\n", "Data:\n", " float64 [dimensionless] (time) [0, 0]\n", "\n", "]\n", " LambdaRequest DataArray [dimensionless] () [\n", "Dimensions: Sizes[time:2, ]\n", "Coordinates:\n", " time datetime64 [ns] (time) [2011-08-12T15:50:16.784000000, 2011-08-12T17:22:05.000000000]\n", "Data:\n", " float64 [Å] (time) [0.533, 0.533]\n", "\n", "]\n", " Ltotal float64 [m] (spectrum) [62.5603, 62.5559, ..., 64.2512, 64.2519]\n", " Phase1 DataArray [dimensionless] () [\n", "Dimensions: Sizes[time:1794, ]\n", "Coordinates:\n", " time datetime64 [ns] (time) [2011-08-12T15:50:16.800000000, 2011-08-12T15:50:18.267000007, ..., 2011-08-12T17:22:02.999218750, 2011-08-12T17:22:05.000000000]\n", "Data:\n", " float64 [µs] (time) [8166.72, 8165.16, ..., 8163.85, 8163.85]\n", "\n", "]\n", " Phase2 DataArray [dimensionless] () [\n", "Dimensions: Sizes[time:1793, ]\n", "Coordinates:\n", " time datetime64 [ns] (time) [2011-08-12T15:50:16.800000000, 2011-08-12T15:50:18.267000007, ..., 2011-08-12T17:22:02.999218750, 2011-08-12T17:22:05.000000000]\n", "Data:\n", " float64 [µs] (time) [8335.63, 8334.09, ..., 8332.86, 8332.86]\n", "\n", "]\n", " Phase3 DataArray [dimensionless] () [\n", "Dimensions: Sizes[time:1777, ]\n", "Coordinates:\n", " time datetime64 [ns] (time) [2011-08-12T15:50:16.800000000, 2011-08-12T15:50:19.688999938, ..., 2011-08-12T17:22:04.421093750, 2011-08-12T17:22:05.000000000]\n", "Data:\n", " float64 [µs] (time) [14001.5, 14001.4, ..., 14002, 14002]\n", "\n", "]\n", " PhaseRequest1 DataArray [dimensionless] () [\n", "Dimensions: Sizes[time:2, ]\n", "Coordinates:\n", " time datetime64 [ns] (time) [2011-08-12T15:50:16.800000000, 2011-08-12T17:22:05.000000000]\n", "Data:\n", " float64 [µs] (time) [8164.07, 8164.07]\n", "\n", "]\n", " PhaseRequest2 DataArray [dimensionless] () [\n", "Dimensions: Sizes[time:2, ]\n", "Coordinates:\n", " time datetime64 [ns] (time) [2011-08-12T15:50:16.800000000, 2011-08-12T17:22:05.000000000]\n", "Data:\n", " float64 [µs] (time) [8332.89, 8332.89]\n", "\n", "]\n", " PhaseRequest3 DataArray [dimensionless] () [\n", "Dimensions: Sizes[time:2, ]\n", "Coordinates:\n", " time datetime64 [ns] (time) [2011-08-12T15:50:16.800000000, 2011-08-12T17:22:05.000000000]\n", "Data:\n", " float64 [µs] (time) [14001.8, 14001.8]\n", "\n", "]\n", " S1HCenter DataArray [dimensionless] () [\n", "Dimensions: Sizes[time:2, ]\n", "Coordinates:\n", " time datetime64 [ns] (time) [2011-08-12T15:50:16.800000000, 2011-08-12T17:22:05.000000000]\n", "Data:\n", " float64 [mm] (time) [0, 0]\n", "\n", "]\n", " S1HCenterOffset DataArray [dimensionless] () [\n", "Dimensions: Sizes[time:2, ]\n", "Coordinates:\n", " time datetime64 [ns] (time) [2011-08-12T15:50:16.800000000, 2011-08-12T17:22:05.000000000]\n", "Data:\n", " float64 [dimensionless] (time) [0, 0]\n", "\n", "]\n", " S1HWidth DataArray [dimensionless] () [\n", "Dimensions: Sizes[time:2, ]\n", "Coordinates:\n", " time datetime64 [ns] (time) [2011-08-12T15:50:16.800000000, 2011-08-12T17:22:05.000000000]\n", "Data:\n", " float64 [mm] (time) [10, 10]\n", "\n", "]\n", " S1VCenter DataArray [dimensionless] () [\n", "Dimensions: Sizes[time:2, ]\n", "Coordinates:\n", " time datetime64 [ns] (time) [2011-08-12T15:50:16.800000000, 2011-08-12T17:22:05.000000000]\n", "Data:\n", " float64 [mm] (time) [5, 5]\n", "\n", "]\n", " S1VCenterOffset DataArray [dimensionless] () [\n", "Dimensions: Sizes[time:2, ]\n", "Coordinates:\n", " time datetime64 [ns] (time) [2011-08-12T15:50:16.800000000, 2011-08-12T17:22:05.000000000]\n", "Data:\n", " float64 [dimensionless] (time) [0, 0]\n", "\n", "]\n", " S1VHeight DataArray [dimensionless] () [\n", "Dimensions: Sizes[time:2, ]\n", "Coordinates:\n", " time datetime64 [ns] (time) [2011-08-12T15:50:16.816000000, 2011-08-12T17:22:05.000000000]\n", "Data:\n", " float64 [mm] (time) [30, 30]\n", "\n", "]\n", " SampleTemp DataArray [dimensionless] () [\n", "Dimensions: Sizes[time:467, ]\n", "Coordinates:\n", " time datetime64 [ns] (time) [2011-08-12T15:50:16.784000000, 2011-08-12T15:50:18.001000007, ..., 2011-08-12T17:21:58.623843750, 2011-08-12T17:22:05.000000000]\n", "Data:\n", " float64 [dimensionless] (time) [299.352, 299.446, ..., 300, 300]\n", "\n", "]\n", " Speed1 DataArray [dimensionless] () [\n", "Dimensions: Sizes[time:2, ]\n", "Coordinates:\n", " time datetime64 [ns] (time) [2011-08-12T15:50:16.800000000, 2011-08-12T17:22:05.000000000]\n", "Data:\n", " float64 [Hz] (time) [60, 60]\n", "\n", "]\n", " Speed2 DataArray [dimensionless] () [\n", "Dimensions: Sizes[time:2, ]\n", "Coordinates:\n", " time datetime64 [ns] (time) [2011-08-12T15:50:16.800000000, 2011-08-12T17:22:05.000000000]\n", "Data:\n", " float64 [Hz] (time) [60, 60]\n", "\n", "]\n", " Speed3 DataArray [dimensionless] () [\n", "Dimensions: Sizes[time:2, ]\n", "Coordinates:\n", " time datetime64 [ns] (time) [2011-08-12T15:50:16.800000000, 2011-08-12T17:22:05.000000000]\n", "Data:\n", " float64 [Hz] (time) [60, 60]\n", "\n", "]\n", " SpeedRequest1 DataArray [dimensionless] () [\n", "Dimensions: Sizes[time:2, ]\n", "Coordinates:\n", " time datetime64 [ns] (time) [2011-08-12T15:50:16.784000000, 2011-08-12T17:22:05.000000000]\n", "Data:\n", " float64 [Hz] (time) [60, 60]\n", "\n", "]\n", " SpeedRequest2 DataArray [dimensionless] () [\n", "Dimensions: Sizes[time:2, ]\n", "Coordinates:\n", " time datetime64 [ns] (time) [2011-08-12T15:50:16.784000000, 2011-08-12T17:22:05.000000000]\n", "Data:\n", " float64 [Hz] (time) [60, 60]\n", "\n", "]\n", " SpeedRequest3 DataArray [dimensionless] () [\n", "Dimensions: Sizes[time:2, ]\n", "Coordinates:\n", " time datetime64 [ns] (time) [2011-08-12T15:50:16.784000000, 2011-08-12T17:22:05.000000000]\n", "Data:\n", " float64 [Hz] (time) [60, 60]\n", "\n", "]\n", " TolRequest DataArray [dimensionless] () [\n", "Dimensions: Sizes[time:2, ]\n", "Coordinates:\n", " time datetime64 [ns] (time) [2011-08-12T15:50:16.800000000, 2011-08-12T17:22:05.000000000]\n", "Data:\n", " float64 [dimensionless] (time) [20, 20]\n", "\n", "]\n", " commErrs DataArray [dimensionless] () [\n", "Dimensions: Sizes[time:2, ]\n", "Coordinates:\n", " time datetime64 [ns] (time) [2011-08-12T15:50:16.784000000, 2011-08-12T17:22:05.000000000]\n", "Data:\n", " float64 [dimensionless] (time) [0, 0]\n", "\n", "]\n", " currentsample DataArray [dimensionless] () [\n", "Dimensions: Sizes[time:2, ]\n", "Coordinates:\n", " time datetime64 [ns] (time) [2011-08-12T15:50:16.784000000, 2011-08-12T17:22:05.000000000]\n", "Data:\n", " float64 [dimensionless] (time) [4, 4]\n", "\n", "]\n", " duration float64 [s] () [5508]\n", " end_time string [dimensionless] () [\"2011-08-12T17:22:05\"]\n", " experiment_identifier string [dimensionless] () [\"IPTS-2767\"]\n", " fernsstatus DataArray [dimensionless] () [\n", "Dimensions: Sizes[time:2, ]\n", "Coordinates:\n", " time datetime64 [ns] (time) [2011-08-12T15:50:16.784000000, 2011-08-12T17:22:05.000000000]\n", "Data:\n", " float64 [dimensionless] (time) [3, 3]\n", "\n", "]\n", " file_notes string [dimensionless] () [\"NONE\"]\n", " frequency DataArray [dimensionless] () [\n", "Dimensions: Sizes[time:330473, ]\n", "Coordinates:\n", " time datetime64 [ns] (time) [2011-08-12T15:50:17.000000000, 2011-08-12T15:50:17.016659999, ..., 2011-08-12T17:22:05.085449218, 2011-08-12T17:22:05.102050781]\n", "Data:\n", " float64 [Hz] (time) [0, 60.024, ..., 60.0024, 59.9988]\n", "\n", "]\n", " gd_prtn_chrg float64 [µAh] () [1171.95]\n", " guide DataArray [dimensionless] () [\n", "Dimensions: Sizes[time:2, ]\n", "Coordinates:\n", " time datetime64 [ns] (time) [2011-08-12T15:50:16.784000000, 2011-08-12T17:22:05.000000000]\n", "Data:\n", " float64 [mm] (time) [-55.463, -55.463]\n", "\n", "]\n", " incident_beam vector3 [m] () [(0, 0, 60)]\n", " instrument_name string [dimensionless] () [\"POWGEN\"]\n", " monitor1 DataArray [dimensionless] () [\n", "Dimensions: Sizes[tof:200001, ]\n", "Coordinates:\n", " position vector3 [m] () [(0, 0, -1)]\n", " source_position vector3 [m] () [(0, 0, -60)]\n", " tof float64 [µs] (tof [bin-edge]) [0, 1, ..., 200000, 200001]\n", "Data:\n", " float64 [counts] (tof) [25, 10, ..., 0, 0] [25, 10, ..., 0, 0]\n", "Attributes:\n", " instrument_name string [dimensionless] () [\"POWGEN\"]\n", " sample_position vector3 [m] () [(0, 0, 0)]\n", " spectrum int32 [dimensionless] () [1]\n", "\n", "]\n", " position vector3 [m] (spectrum) [(1.17451, -1.01106, -2.03797), (1.18148, -0.959466, -2.05334), ..., (1.81429, 0.0956584, 3.84338), (1.81375, 0.149937, 3.8427)]\n", " proton_charge DataArray [dimensionless] () [\n", "Dimensions: Sizes[time:330473, ]\n", "Coordinates:\n", " time datetime64 [ns] (time) [2011-08-12T15:50:17.000000000, 2011-08-12T15:50:17.016659999, ..., 2011-08-12T17:22:05.085449218, 2011-08-12T17:22:05.102050781]\n", "Data:\n", " float64 [pC] (time) [1.48434e+07, 1.48434e+07, ..., 1.48733e+07, 1.48391e+07]\n", "\n", "]\n", " run_number string [dimensionless] () [\"4844\"]\n", " run_start string [dimensionless] () [\"2011-08-12T15:50:17\"]\n", " run_title string [dimensionless] () [\"diamond cw0.533 4.22e12 60Hz [10x30]\"]\n", " running DataArray [dimensionless] () [\n", "Dimensions: Sizes[time:1, ]\n", "Coordinates:\n", " time datetime64 [ns] (time) [1990-01-01T00:00:00.000000000]\n", "Data:\n", " bool [dimensionless] (time) [True]\n", "\n", "]\n", " s1b DataArray [dimensionless] () [\n", "Dimensions: Sizes[time:2, ]\n", "Coordinates:\n", " time datetime64 [ns] (time) [2011-08-12T15:50:16.784000000, 2011-08-12T17:22:05.000000000]\n", "Data:\n", " float64 [mm] (time) [20, 20]\n", "\n", "]\n", " s1l DataArray [dimensionless] () [\n", "Dimensions: Sizes[time:2, ]\n", "Coordinates:\n", " time datetime64 [ns] (time) [2011-08-12T15:50:16.800000000, 2011-08-12T17:22:05.000000000]\n", "Data:\n", " float64 [mm] (time) [5, 5]\n", "\n", "]\n", " s1r DataArray [dimensionless] () [\n", "Dimensions: Sizes[time:2, ]\n", "Coordinates:\n", " time datetime64 [ns] (time) [2011-08-12T15:50:16.784000000, 2011-08-12T17:22:05.000000000]\n", "Data:\n", " float64 [mm] (time) [-5, -5]\n", "\n", "]\n", " s1t DataArray [dimensionless] () [\n", "Dimensions: Sizes[time:2, ]\n", "Coordinates:\n", " time datetime64 [ns] (time) [2011-08-12T15:50:16.784000000, 2011-08-12T17:22:05.000000000]\n", "Data:\n", " float64 [mm] (time) [-10, -10]\n", "\n", "]\n", " sample PyObject [dimensionless] () []\n", " sample_position vector3 [m] () [(0, 0, 0)]\n", " samplerequest DataArray [dimensionless] () [\n", "Dimensions: Sizes[time:2, ]\n", "Coordinates:\n", " time datetime64 [ns] (time) [2011-08-12T15:50:16.784000000, 2011-08-12T17:22:05.000000000]\n", "Data:\n", " float64 [dimensionless] (time) [4, 4]\n", "\n", "]\n", " scattered_beam vector3 [m] (spectrum) [(1.17451, -1.01106, -2.03797), (1.18148, -0.959466, -2.05334), ..., (1.81429, 0.0956584, 3.84338), (1.81375, 0.149937, 3.8427)]\n", " source_position vector3 [m] () [(0, 0, -60)]\n", " start_time string [dimensionless] () [\"2011-08-12T15:50:17\"]\n", " tof float64 [µs] (dspacing [bin-edge]) [19, 16694.3]\n", " two_theta float64 [rad] (spectrum) [2.49144, 2.50373, ..., 0.441579, 0.442313]\n", " vGuide DataArray [dimensionless] () [\n", "Dimensions: Sizes[time:2, ]\n", "Coordinates:\n", " time datetime64 [ns] (time) [2011-08-12T15:50:16.800000000, 2011-08-12T17:22:05.000000000]\n", "Data:\n", " float64 [dimensionless] (time) [2, 2]\n", "\n", "]\n", " veto_pulse_time DataArray [dimensionless] () [\n", "Dimensions: Sizes[time:1, ]\n", "Coordinates:\n", " time datetime64 [ns] (time) [2011-08-12T15:50:17.000000000]\n", "Data:\n", " float64 [dimensionless] (time) [0]\n", "\n", "]\n" ] }, "execution_count": 10, "metadata": {}, "output_type": "execute_result" } ], "source": [ "dspacing_vanadium = scn.convert(vanadium, 'tof', 'dspacing', scatter=True)\n", "dspacing_sample = scn.convert(sample, 'tof', 'dspacing', scatter=True)\n", "dspacing_sample" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "### Neutron monitors\n", "\n", "*This is an optional section.\n", "The next section does not use the monitor-normalized data produced here.\n", "This section could thus be skipped.*\n", "\n", "If available, neutron monitors are stored as attributes of a data array:" ] }, { "cell_type": "code", "execution_count": 11, "metadata": {}, "outputs": [ { "data": { "text/html": [ "
\n", "\n", "\n", "\n", "\n", "Show/Hide data repr\n", "\n", "\n", "\n", "\n", "\n", "Show/Hide attributes\n", "\n", "\n", "\n", "\n", "\n", "\n", "
scipp.DataArray (4.58 MB)
    • tof: 200001
    • position
      ()
      vector3
      m
      [ 0. 0. -1.]
      Values:
      array([ 0.,  0., -1.])
    • source_position
      ()
      vector3
      m
      [ 0. 0. -60.]
      Values:
      array([  0.,   0., -60.])
    • tof
      (tof [bin-edge])
      float64
      µs
      0.0, 1.0, ..., 2.000e+05, 2.000e+05
      Values:
      array([0.00000e+00, 1.00000e+00, 2.00000e+00, ..., 1.99999e+05,\n",
       "       2.00000e+05, 2.00001e+05])
    • (tof)
      float64
      counts
      25.0, 10.0, ..., 0.0, 0.0
      σ = 5.0, 3.162, ..., 0.0, 0.0
      Values:
      array([25., 10., 12., ...,  0.,  0.,  0.])

      Variances (σ²):
      array([25., 10., 12., ...,  0.,  0.,  0.])
    • instrument_name
      ()
      string
      POWGEN
      Values:
      'POWGEN'
    • sample_position
      ()
      vector3
      m
      [0. 0. 0.]
      Values:
      array([0., 0., 0.])
    • spectrum
      ()
      int32
      1
      Values:
      array(1, dtype=int32)
" ], "text/plain": [ "\n", "Dimensions: Sizes[tof:200001, ]\n", "Coordinates:\n", " position vector3 [m] () [(0, 0, -1)]\n", " source_position vector3 [m] () [(0, 0, -60)]\n", " tof float64 [µs] (tof [bin-edge]) [0, 1, ..., 200000, 200001]\n", "Data:\n", " float64 [counts] (tof) [25, 10, ..., 0, 0] [25, 10, ..., 0, 0]\n", "Attributes:\n", " instrument_name string [dimensionless] () [\"POWGEN\"]\n", " sample_position vector3 [m] () [(0, 0, 0)]\n", " spectrum int32 [dimensionless] () [1]\n" ] }, "execution_count": 11, "metadata": {}, "output_type": "execute_result" } ], "source": [ "mon = sample.attrs['monitor1'].value\n", "mon" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "The monitor could, e.g., be used to normalize the data.\n", "To do so, both data and monitor need to be converted to a unit that accounts for differing flight paths, e.g., wavelength or energy:" ] }, { "cell_type": "code", "execution_count": 12, "metadata": {}, "outputs": [], "source": [ "sample_lambda = scn.convert(sample, 'tof', 'wavelength', scatter=True)\n", "mon = scn.convert(mon, 'tof', 'wavelength', scatter=False)" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "The sample data is in event-mode, i.e., is not histogrammed.\n", "Event data *can* be divided by a histogram (such as `mon` in this case), using a specialized function for scaling (see [Binned data](https://scipp.github.io/user-guide/binned-data.html)).\n", "First we rebin the monitor since the original binning is very fine:" ] }, { "cell_type": "code", "execution_count": 13, "metadata": {}, "outputs": [ { "data": { "text/html": [ "
\n", "\n", "\n", "\n", "\n", "Show/Hide data repr\n", "\n", "\n", "\n", "\n", "\n", "Show/Hide attributes\n", "\n", "\n", "\n", "\n", "\n", "\n", "
scipp.DataArray (23.54 KB)
    • wavelength: 999
    • wavelength
      (wavelength [bin-edge])
      float64
      Å
      0.0, 0.001, ..., 0.999, 1.0
      Values:
      array([0.        , 0.001001  , 0.002002  , 0.003003  , 0.004004  ,\n",
       "       0.00500501, 0.00600601, 0.00700701, 0.00800801, 0.00900901,\n",
       "       0.01001001, 0.01101101, 0.01201201, 0.01301301, 0.01401401,\n",
       "       0.01501502, 0.01601602, 0.01701702, 0.01801802, 0.01901902,\n",
       "       0.02002002, 0.02102102, 0.02202202, 0.02302302, 0.02402402,\n",
       "       0.02502503, 0.02602603, 0.02702703, 0.02802803, 0.02902903,\n",
       "       0.03003003, 0.03103103, 0.03203203, 0.03303303, 0.03403403,\n",
       "       0.03503504, 0.03603604, 0.03703704, 0.03803804, 0.03903904,\n",
       "       0.04004004, 0.04104104, 0.04204204, 0.04304304, 0.04404404,\n",
       "       0.04504505, 0.04604605, 0.04704705, 0.04804805, 0.04904905,\n",
       "       0.05005005, 0.05105105, 0.05205205, 0.05305305, 0.05405405,\n",
       "       0.05505506, 0.05605606, 0.05705706, 0.05805806, 0.05905906,\n",
       "       0.06006006, 0.06106106, 0.06206206, 0.06306306, 0.06406406,\n",
       "       0.06506507, 0.06606607, 0.06706707, 0.06806807, 0.06906907,\n",
       "       0.07007007, 0.07107107, 0.07207207, 0.07307307, 0.07407407,\n",
       "       0.07507508, 0.07607608, 0.07707708, 0.07807808, 0.07907908,\n",
       "       0.08008008, 0.08108108, 0.08208208, 0.08308308, 0.08408408,\n",
       "       0.08508509, 0.08608609, 0.08708709, 0.08808809, 0.08908909,\n",
       "       0.09009009, 0.09109109, 0.09209209, 0.09309309, 0.09409409,\n",
       "       0.0950951 , 0.0960961 , 0.0970971 , 0.0980981 , 0.0990991 ,\n",
       "       0.1001001 , 0.1011011 , 0.1021021 , 0.1031031 , 0.1041041 ,\n",
       "       0.10510511, 0.10610611, 0.10710711, 0.10810811, 0.10910911,\n",
       "       0.11011011, 0.11111111, 0.11211211, 0.11311311, 0.11411411,\n",
       "       0.11511512, 0.11611612, 0.11711712, 0.11811812, 0.11911912,\n",
       "       0.12012012, 0.12112112, 0.12212212, 0.12312312, 0.12412412,\n",
       "       0.12512513, 0.12612613, 0.12712713, 0.12812813, 0.12912913,\n",
       "       0.13013013, 0.13113113, 0.13213213, 0.13313313, 0.13413413,\n",
       "       0.13513514, 0.13613614, 0.13713714, 0.13813814, 0.13913914,\n",
       "       0.14014014, 0.14114114, 0.14214214, 0.14314314, 0.14414414,\n",
       "       0.14514515, 0.14614615, 0.14714715, 0.14814815, 0.14914915,\n",
       "       0.15015015, 0.15115115, 0.15215215, 0.15315315, 0.15415415,\n",
       "       0.15515516, 0.15615616, 0.15715716, 0.15815816, 0.15915916,\n",
       "       0.16016016, 0.16116116, 0.16216216, 0.16316316, 0.16416416,\n",
       "       0.16516517, 0.16616617, 0.16716717, 0.16816817, 0.16916917,\n",
       "       0.17017017, 0.17117117, 0.17217217, 0.17317317, 0.17417417,\n",
       "       0.17517518, 0.17617618, 0.17717718, 0.17817818, 0.17917918,\n",
       "       0.18018018, 0.18118118, 0.18218218, 0.18318318, 0.18418418,\n",
       "       0.18518519, 0.18618619, 0.18718719, 0.18818819, 0.18918919,\n",
       "       0.19019019, 0.19119119, 0.19219219, 0.19319319, 0.19419419,\n",
       "       0.1951952 , 0.1961962 , 0.1971972 , 0.1981982 , 0.1991992 ,\n",
       "       0.2002002 , 0.2012012 , 0.2022022 , 0.2032032 , 0.2042042 ,\n",
       "       0.20520521, 0.20620621, 0.20720721, 0.20820821, 0.20920921,\n",
       "       0.21021021, 0.21121121, 0.21221221, 0.21321321, 0.21421421,\n",
       "       0.21521522, 0.21621622, 0.21721722, 0.21821822, 0.21921922,\n",
       "       0.22022022, 0.22122122, 0.22222222, 0.22322322, 0.22422422,\n",
       "       0.22522523, 0.22622623, 0.22722723, 0.22822823, 0.22922923,\n",
       "       0.23023023, 0.23123123, 0.23223223, 0.23323323, 0.23423423,\n",
       "       0.23523524, 0.23623624, 0.23723724, 0.23823824, 0.23923924,\n",
       "       0.24024024, 0.24124124, 0.24224224, 0.24324324, 0.24424424,\n",
       "       0.24524525, 0.24624625, 0.24724725, 0.24824825, 0.24924925,\n",
       "       0.25025025, 0.25125125, 0.25225225, 0.25325325, 0.25425425,\n",
       "       0.25525526, 0.25625626, 0.25725726, 0.25825826, 0.25925926,\n",
       "       0.26026026, 0.26126126, 0.26226226, 0.26326326, 0.26426426,\n",
       "       0.26526527, 0.26626627, 0.26726727, 0.26826827, 0.26926927,\n",
       "       0.27027027, 0.27127127, 0.27227227, 0.27327327, 0.27427427,\n",
       "       0.27527528, 0.27627628, 0.27727728, 0.27827828, 0.27927928,\n",
       "       0.28028028, 0.28128128, 0.28228228, 0.28328328, 0.28428428,\n",
       "       0.28528529, 0.28628629, 0.28728729, 0.28828829, 0.28928929,\n",
       "       0.29029029, 0.29129129, 0.29229229, 0.29329329, 0.29429429,\n",
       "       0.2952953 , 0.2962963 , 0.2972973 , 0.2982983 , 0.2992993 ,\n",
       "       0.3003003 , 0.3013013 , 0.3023023 , 0.3033033 , 0.3043043 ,\n",
       "       0.30530531, 0.30630631, 0.30730731, 0.30830831, 0.30930931,\n",
       "       0.31031031, 0.31131131, 0.31231231, 0.31331331, 0.31431431,\n",
       "       0.31531532, 0.31631632, 0.31731732, 0.31831832, 0.31931932,\n",
       "       0.32032032, 0.32132132, 0.32232232, 0.32332332, 0.32432432,\n",
       "       0.32532533, 0.32632633, 0.32732733, 0.32832833, 0.32932933,\n",
       "       0.33033033, 0.33133133, 0.33233233, 0.33333333, 0.33433433,\n",
       "       0.33533534, 0.33633634, 0.33733734, 0.33833834, 0.33933934,\n",
       "       0.34034034, 0.34134134, 0.34234234, 0.34334334, 0.34434434,\n",
       "       0.34534535, 0.34634635, 0.34734735, 0.34834835, 0.34934935,\n",
       "       0.35035035, 0.35135135, 0.35235235, 0.35335335, 0.35435435,\n",
       "       0.35535536, 0.35635636, 0.35735736, 0.35835836, 0.35935936,\n",
       "       0.36036036, 0.36136136, 0.36236236, 0.36336336, 0.36436436,\n",
       "       0.36536537, 0.36636637, 0.36736737, 0.36836837, 0.36936937,\n",
       "       0.37037037, 0.37137137, 0.37237237, 0.37337337, 0.37437437,\n",
       "       0.37537538, 0.37637638, 0.37737738, 0.37837838, 0.37937938,\n",
       "       0.38038038, 0.38138138, 0.38238238, 0.38338338, 0.38438438,\n",
       "       0.38538539, 0.38638639, 0.38738739, 0.38838839, 0.38938939,\n",
       "       0.39039039, 0.39139139, 0.39239239, 0.39339339, 0.39439439,\n",
       "       0.3953954 , 0.3963964 , 0.3973974 , 0.3983984 , 0.3993994 ,\n",
       "       0.4004004 , 0.4014014 , 0.4024024 , 0.4034034 , 0.4044044 ,\n",
       "       0.40540541, 0.40640641, 0.40740741, 0.40840841, 0.40940941,\n",
       "       0.41041041, 0.41141141, 0.41241241, 0.41341341, 0.41441441,\n",
       "       0.41541542, 0.41641642, 0.41741742, 0.41841842, 0.41941942,\n",
       "       0.42042042, 0.42142142, 0.42242242, 0.42342342, 0.42442442,\n",
       "       0.42542543, 0.42642643, 0.42742743, 0.42842843, 0.42942943,\n",
       "       0.43043043, 0.43143143, 0.43243243, 0.43343343, 0.43443443,\n",
       "       0.43543544, 0.43643644, 0.43743744, 0.43843844, 0.43943944,\n",
       "       0.44044044, 0.44144144, 0.44244244, 0.44344344, 0.44444444,\n",
       "       0.44544545, 0.44644645, 0.44744745, 0.44844845, 0.44944945,\n",
       "       0.45045045, 0.45145145, 0.45245245, 0.45345345, 0.45445445,\n",
       "       0.45545546, 0.45645646, 0.45745746, 0.45845846, 0.45945946,\n",
       "       0.46046046, 0.46146146, 0.46246246, 0.46346346, 0.46446446,\n",
       "       0.46546547, 0.46646647, 0.46746747, 0.46846847, 0.46946947,\n",
       "       0.47047047, 0.47147147, 0.47247247, 0.47347347, 0.47447447,\n",
       "       0.47547548, 0.47647648, 0.47747748, 0.47847848, 0.47947948,\n",
       "       0.48048048, 0.48148148, 0.48248248, 0.48348348, 0.48448448,\n",
       "       0.48548549, 0.48648649, 0.48748749, 0.48848849, 0.48948949,\n",
       "       0.49049049, 0.49149149, 0.49249249, 0.49349349, 0.49449449,\n",
       "       0.4954955 , 0.4964965 , 0.4974975 , 0.4984985 , 0.4994995 ,\n",
       "       0.5005005 , 0.5015015 , 0.5025025 , 0.5035035 , 0.5045045 ,\n",
       "       0.50550551, 0.50650651, 0.50750751, 0.50850851, 0.50950951,\n",
       "       0.51051051, 0.51151151, 0.51251251, 0.51351351, 0.51451451,\n",
       "       0.51551552, 0.51651652, 0.51751752, 0.51851852, 0.51951952,\n",
       "       0.52052052, 0.52152152, 0.52252252, 0.52352352, 0.52452452,\n",
       "       0.52552553, 0.52652653, 0.52752753, 0.52852853, 0.52952953,\n",
       "       0.53053053, 0.53153153, 0.53253253, 0.53353353, 0.53453453,\n",
       "       0.53553554, 0.53653654, 0.53753754, 0.53853854, 0.53953954,\n",
       "       0.54054054, 0.54154154, 0.54254254, 0.54354354, 0.54454454,\n",
       "       0.54554555, 0.54654655, 0.54754755, 0.54854855, 0.54954955,\n",
       "       0.55055055, 0.55155155, 0.55255255, 0.55355355, 0.55455455,\n",
       "       0.55555556, 0.55655656, 0.55755756, 0.55855856, 0.55955956,\n",
       "       0.56056056, 0.56156156, 0.56256256, 0.56356356, 0.56456456,\n",
       "       0.56556557, 0.56656657, 0.56756757, 0.56856857, 0.56956957,\n",
       "       0.57057057, 0.57157157, 0.57257257, 0.57357357, 0.57457457,\n",
       "       0.57557558, 0.57657658, 0.57757758, 0.57857858, 0.57957958,\n",
       "       0.58058058, 0.58158158, 0.58258258, 0.58358358, 0.58458458,\n",
       "       0.58558559, 0.58658659, 0.58758759, 0.58858859, 0.58958959,\n",
       "       0.59059059, 0.59159159, 0.59259259, 0.59359359, 0.59459459,\n",
       "       0.5955956 , 0.5965966 , 0.5975976 , 0.5985986 , 0.5995996 ,\n",
       "       0.6006006 , 0.6016016 , 0.6026026 , 0.6036036 , 0.6046046 ,\n",
       "       0.60560561, 0.60660661, 0.60760761, 0.60860861, 0.60960961,\n",
       "       0.61061061, 0.61161161, 0.61261261, 0.61361361, 0.61461461,\n",
       "       0.61561562, 0.61661662, 0.61761762, 0.61861862, 0.61961962,\n",
       "       0.62062062, 0.62162162, 0.62262262, 0.62362362, 0.62462462,\n",
       "       0.62562563, 0.62662663, 0.62762763, 0.62862863, 0.62962963,\n",
       "       0.63063063, 0.63163163, 0.63263263, 0.63363363, 0.63463463,\n",
       "       0.63563564, 0.63663664, 0.63763764, 0.63863864, 0.63963964,\n",
       "       0.64064064, 0.64164164, 0.64264264, 0.64364364, 0.64464464,\n",
       "       0.64564565, 0.64664665, 0.64764765, 0.64864865, 0.64964965,\n",
       "       0.65065065, 0.65165165, 0.65265265, 0.65365365, 0.65465465,\n",
       "       0.65565566, 0.65665666, 0.65765766, 0.65865866, 0.65965966,\n",
       "       0.66066066, 0.66166166, 0.66266266, 0.66366366, 0.66466466,\n",
       "       0.66566567, 0.66666667, 0.66766767, 0.66866867, 0.66966967,\n",
       "       0.67067067, 0.67167167, 0.67267267, 0.67367367, 0.67467467,\n",
       "       0.67567568, 0.67667668, 0.67767768, 0.67867868, 0.67967968,\n",
       "       0.68068068, 0.68168168, 0.68268268, 0.68368368, 0.68468468,\n",
       "       0.68568569, 0.68668669, 0.68768769, 0.68868869, 0.68968969,\n",
       "       0.69069069, 0.69169169, 0.69269269, 0.69369369, 0.69469469,\n",
       "       0.6956957 , 0.6966967 , 0.6976977 , 0.6986987 , 0.6996997 ,\n",
       "       0.7007007 , 0.7017017 , 0.7027027 , 0.7037037 , 0.7047047 ,\n",
       "       0.70570571, 0.70670671, 0.70770771, 0.70870871, 0.70970971,\n",
       "       0.71071071, 0.71171171, 0.71271271, 0.71371371, 0.71471471,\n",
       "       0.71571572, 0.71671672, 0.71771772, 0.71871872, 0.71971972,\n",
       "       0.72072072, 0.72172172, 0.72272272, 0.72372372, 0.72472472,\n",
       "       0.72572573, 0.72672673, 0.72772773, 0.72872873, 0.72972973,\n",
       "       0.73073073, 0.73173173, 0.73273273, 0.73373373, 0.73473473,\n",
       "       0.73573574, 0.73673674, 0.73773774, 0.73873874, 0.73973974,\n",
       "       0.74074074, 0.74174174, 0.74274274, 0.74374374, 0.74474474,\n",
       "       0.74574575, 0.74674675, 0.74774775, 0.74874875, 0.74974975,\n",
       "       0.75075075, 0.75175175, 0.75275275, 0.75375375, 0.75475475,\n",
       "       0.75575576, 0.75675676, 0.75775776, 0.75875876, 0.75975976,\n",
       "       0.76076076, 0.76176176, 0.76276276, 0.76376376, 0.76476476,\n",
       "       0.76576577, 0.76676677, 0.76776777, 0.76876877, 0.76976977,\n",
       "       0.77077077, 0.77177177, 0.77277277, 0.77377377, 0.77477477,\n",
       "       0.77577578, 0.77677678, 0.77777778, 0.77877878, 0.77977978,\n",
       "       0.78078078, 0.78178178, 0.78278278, 0.78378378, 0.78478478,\n",
       "       0.78578579, 0.78678679, 0.78778779, 0.78878879, 0.78978979,\n",
       "       0.79079079, 0.79179179, 0.79279279, 0.79379379, 0.79479479,\n",
       "       0.7957958 , 0.7967968 , 0.7977978 , 0.7987988 , 0.7997998 ,\n",
       "       0.8008008 , 0.8018018 , 0.8028028 , 0.8038038 , 0.8048048 ,\n",
       "       0.80580581, 0.80680681, 0.80780781, 0.80880881, 0.80980981,\n",
       "       0.81081081, 0.81181181, 0.81281281, 0.81381381, 0.81481481,\n",
       "       0.81581582, 0.81681682, 0.81781782, 0.81881882, 0.81981982,\n",
       "       0.82082082, 0.82182182, 0.82282282, 0.82382382, 0.82482482,\n",
       "       0.82582583, 0.82682683, 0.82782783, 0.82882883, 0.82982983,\n",
       "       0.83083083, 0.83183183, 0.83283283, 0.83383383, 0.83483483,\n",
       "       0.83583584, 0.83683684, 0.83783784, 0.83883884, 0.83983984,\n",
       "       0.84084084, 0.84184184, 0.84284284, 0.84384384, 0.84484484,\n",
       "       0.84584585, 0.84684685, 0.84784785, 0.84884885, 0.84984985,\n",
       "       0.85085085, 0.85185185, 0.85285285, 0.85385385, 0.85485485,\n",
       "       0.85585586, 0.85685686, 0.85785786, 0.85885886, 0.85985986,\n",
       "       0.86086086, 0.86186186, 0.86286286, 0.86386386, 0.86486486,\n",
       "       0.86586587, 0.86686687, 0.86786787, 0.86886887, 0.86986987,\n",
       "       0.87087087, 0.87187187, 0.87287287, 0.87387387, 0.87487487,\n",
       "       0.87587588, 0.87687688, 0.87787788, 0.87887888, 0.87987988,\n",
       "       0.88088088, 0.88188188, 0.88288288, 0.88388388, 0.88488488,\n",
       "       0.88588589, 0.88688689, 0.88788789, 0.88888889, 0.88988989,\n",
       "       0.89089089, 0.89189189, 0.89289289, 0.89389389, 0.89489489,\n",
       "       0.8958959 , 0.8968969 , 0.8978979 , 0.8988989 , 0.8998999 ,\n",
       "       0.9009009 , 0.9019019 , 0.9029029 , 0.9039039 , 0.9049049 ,\n",
       "       0.90590591, 0.90690691, 0.90790791, 0.90890891, 0.90990991,\n",
       "       0.91091091, 0.91191191, 0.91291291, 0.91391391, 0.91491491,\n",
       "       0.91591592, 0.91691692, 0.91791792, 0.91891892, 0.91991992,\n",
       "       0.92092092, 0.92192192, 0.92292292, 0.92392392, 0.92492492,\n",
       "       0.92592593, 0.92692693, 0.92792793, 0.92892893, 0.92992993,\n",
       "       0.93093093, 0.93193193, 0.93293293, 0.93393393, 0.93493493,\n",
       "       0.93593594, 0.93693694, 0.93793794, 0.93893894, 0.93993994,\n",
       "       0.94094094, 0.94194194, 0.94294294, 0.94394394, 0.94494494,\n",
       "       0.94594595, 0.94694695, 0.94794795, 0.94894895, 0.94994995,\n",
       "       0.95095095, 0.95195195, 0.95295295, 0.95395395, 0.95495495,\n",
       "       0.95595596, 0.95695696, 0.95795796, 0.95895896, 0.95995996,\n",
       "       0.96096096, 0.96196196, 0.96296296, 0.96396396, 0.96496496,\n",
       "       0.96596597, 0.96696697, 0.96796797, 0.96896897, 0.96996997,\n",
       "       0.97097097, 0.97197197, 0.97297297, 0.97397397, 0.97497497,\n",
       "       0.97597598, 0.97697698, 0.97797798, 0.97897898, 0.97997998,\n",
       "       0.98098098, 0.98198198, 0.98298298, 0.98398398, 0.98498498,\n",
       "       0.98598599, 0.98698699, 0.98798799, 0.98898899, 0.98998999,\n",
       "       0.99099099, 0.99199199, 0.99299299, 0.99399399, 0.99499499,\n",
       "       0.995996  , 0.996997  , 0.997998  , 0.998999  , 1.        ])
    • (wavelength)
      float64
      counts
      194.502, 214.933, ..., 206.513, 193.382
      σ = 13.946, 14.661, ..., 14.371, 13.906
      Values:
      array([194.50198662, 214.93282848, 215.43084186, 222.58902438,\n",
       "       194.98808026, 184.28855232, 210.28457907, 205.01589299,\n",
       "       203.80245869, 221.93680173, 224.27265933, 219.71144768,\n",
       "       205.2213808 , 226.06319827, 216.28855232, 224.28855232,\n",
       "       227.82232493, 213.19356806, 192.14626278, 230.03178598,\n",
       "       259.74285997, 197.6641424 , 221.05562547, 203.82195123,\n",
       "       227.33983085, 200.36764358, 200.14626278, 199.13831629,\n",
       "       197.58542483, 215.17767507, 212.09461056, 222.07549171,\n",
       "       197.19716762, 194.43084186, 196.47057434, 241.77899289,\n",
       "       229.24881984, 198.00397325, 211.45468134, 199.01589299,\n",
       "       220.37521639, 209.24162073, 200.16935488, 188.81003149,\n",
       "       211.96858771, 184.72374112, 201.24522029, 209.52507872,\n",
       "       192.48646733, 227.57387878, 196.97541312, 215.29799359,\n",
       "       234.19642023, 205.3203383 , 193.86168371, 230.81400474,\n",
       "       236.18562157, 217.84939027, 192.57313139, 216.96101491,\n",
       "       198.46660109, 227.7710464 , 197.75477971, 192.48609363,\n",
       "       222.37956333, 215.5417191 , 224.25714003, 214.51315898,\n",
       "       234.05525178, 211.35647123, 217.02706535, 194.02383949,\n",
       "       214.72014156, 230.51825331, 203.65110157, 216.78296614,\n",
       "       208.61571604, 204.96536185, 213.27228563, 234.07511802,\n",
       "       197.02421318, 208.62440992, 173.76669946, 198.95666797,\n",
       "       217.53054675, 206.31127072, 219.36404403, 192.19791501,\n",
       "       224.24087334, 194.94362714, 214.59846566, 211.34268301,\n",
       "       204.84219117, 197.37086944, 219.45977567, 191.96784032,\n",
       "       197.53699847, 200.16178208, 211.88304473, 214.69158144,\n",
       "       211.46262784, 215.86093632, 226.58865069, 197.27700627,\n",
       "       229.28383168, 207.04047987, 210.24919354, 225.15780883,\n",
       "       249.19716762, 231.01664038, 222.17295444, 215.37633747,\n",
       "       206.74360736, 193.20076717, 209.73056653, 199.03575923,\n",
       "       224.21740755, 238.24199443, 216.3279111 , 218.93282848,\n",
       "       204.71579462, 216.47020064, 200.61894189, 205.10442553,\n",
       "       214.96858771, 227.61174279, 209.5800942 , 232.60590107,\n",
       "       230.18674265, 203.71467354, 229.72411481, 214.28383168,\n",
       "       212.23975226, 204.70139641, 200.27985844, 216.14626278,\n",
       "       220.893096  , 222.12751763, 221.63992922, 218.84504333,\n",
       "       240.16290316, 222.79885913, 214.34057824, 204.85186874,\n",
       "       215.34057824, 218.30084576, 222.09970489, 214.2365264 ,\n",
       "       211.85883155, 216.52110547, 233.66227393, 176.21852864,\n",
       "       219.70312749, 202.98770656, 205.35460276, 230.27737997,\n",
       "       224.00186847, 212.79053894, 209.72659328, 200.38713613,\n",
       "       214.52942566, 214.56990553, 215.14303692, 219.38713613,\n",
       "       206.05165222, 199.14948864, 203.81400474, 203.05165222,\n",
       "       220.8624311 , 231.00322586, 176.25887111, 192.02818643,\n",
       "       244.5114279 , 201.55363885, 216.21268692, 190.07511802,\n",
       "       230.85298983, 213.95232103, 196.82952403, 198.85957894,\n",
       "       186.88117626, 189.33337913, 226.73888672, 206.88117626,\n",
       "       204.56742706, 201.75837927, 214.56331642, 196.80891041,\n",
       "       203.94946886, 212.07511802, 193.89781664, 193.93282848,\n",
       "       190.12564916, 184.34417779, 195.6948073 , 206.23292685,\n",
       "       182.3131392 , 205.23292685, 216.74398105, 219.85920524,\n",
       "       230.14874125, 198.63273011, 253.37521639, 205.08516929,\n",
       "       189.84616441, 205.08058604, 198.4579072 , 213.74248627,\n",
       "       233.81152627, 187.36156556, 225.61174279, 189.87322976,\n",
       "       192.86168371, 213.41172301, 192.38452027, 192.35969709,\n",
       "       200.48286777, 205.76520468, 203.66961043, 180.91370963,\n",
       "       217.46014937, 188.48422516, 198.95456319, 183.33796238,\n",
       "       194.12354438, 212.09026362, 204.07511802, 206.19108959,\n",
       "       200.49838707, 204.98906395, 189.18749004, 209.99280089,\n",
       "       216.86168371, 219.10106228, 172.374469  , 175.15869361,\n",
       "       224.81250996, 214.321322  , 208.39844587, 192.08951623,\n",
       "       193.80456346, 204.77974028, 213.35287168, 194.92600307,\n",
       "       207.65147527, 205.44734485, 201.57313139, 180.27069194,\n",
       "       208.01005127, 192.50769095, 212.21740755, 204.93259218,\n",
       "       198.50110184, 198.28929971, 228.57897311, 201.57201031,\n",
       "       205.86205741, 193.14116845, 200.72958284, 168.8475218 ,\n",
       "       221.44463007, 167.49254535, 208.0841856 , 175.69728577,\n",
       "       191.16402425, 206.30159315, 169.86168371, 212.92848154,\n",
       "       187.93282848, 189.76917793, 220.37671116, 179.05810394,\n",
       "       182.10479922, 191.83634945, 198.32220678, 178.69840685,\n",
       "       182.41817472, 197.69854424, 205.74459105, 208.39297785,\n",
       "       193.77389856, 184.12130221, 202.34244671, 178.5654212 ,\n",
       "       203.7021438 , 196.32244308, 169.12466546, 199.82045645,\n",
       "       193.58542483, 209.17916985, 172.97106618, 207.30084576,\n",
       "       184.54022433, 172.21740755, 179.22970099, 167.32281677,\n",
       "       201.18973221, 185.7250985 , 188.71368985, 196.74968538,\n",
       "       199.03080229, 187.77773442, 198.54681343, 164.77203009,\n",
       "       162.13845368, 209.33113696, 199.61571604, 167.81288365,\n",
       "       188.40241912, 207.95517319, 189.1299961 , 187.86378849,\n",
       "       158.82307232, 183.55686471, 172.64614463, 182.93096001,\n",
       "       189.28855232, 168.53662477, 161.34529887, 194.51638484,\n",
       "       175.91843027, 178.76571576, 166.60837954, 187.37993702,\n",
       "       182.74944908, 178.58355636, 181.48174669, 175.61758451,\n",
       "       191.48660472, 172.32814741, 171.62403623, 172.76632576,\n",
       "       170.83499206, 186.55041299, 189.43332033, 205.52769458,\n",
       "       167.83747053, 179.16748641, 165.87893409, 168.69817055,\n",
       "       165.8092841 , 173.44724593, 169.60912693, 168.76235251,\n",
       "       177.35969709, 189.43097925, 180.22671143, 173.38788352,\n",
       "       163.88178625, 178.47380019, 165.75702188, 161.27378041,\n",
       "       162.98920134, 178.26036589, 169.65904807, 179.30594009,\n",
       "       159.27116455, 164.81872537, 186.97157727, 162.62688839,\n",
       "       181.37049575, 178.91829287, 177.89758033, 161.62291514,\n",
       "       159.14626278, 171.26321805, 188.13943737, 178.67956278,\n",
       "       152.26976868, 175.70289119, 151.21740755, 148.81984646,\n",
       "       167.76296251, 160.24498398, 183.00397325, 175.25527156,\n",
       "       163.28794232, 172.14687278, 155.18127463, 176.00845758,\n",
       "       161.35633384, 170.14626278, 167.64315508, 170.32977958,\n",
       "       171.07511802, 155.53699847, 160.61323757, 173.47094803,\n",
       "       162.97802898, 172.00159369, 146.66750565, 172.9237609 ,\n",
       "       154.32865849, 183.27948474, 148.84878028, 166.83064512,\n",
       "       169.15519297, 165.91978765, 196.27055455, 173.13322196,\n",
       "       180.24348921, 175.67629845, 165.14527909, 163.11620788,\n",
       "       161.17631769, 171.14626278, 187.01402452, 179.41382778,\n",
       "       164.91272594, 178.2444729 , 143.66526348, 175.25055092,\n",
       "       145.52297395, 147.23839487, 159.10428814, 144.67344628,\n",
       "       170.56892184, 155.9307237 , 169.36180186, 150.62305253,\n",
       "       174.62142036, 158.30977595, 175.59622349, 155.24236812,\n",
       "       140.24820984, 176.99813153, 151.01565669, 149.21156583,\n",
       "       168.07511802, 159.66736826, 169.49268274, 155.73800194,\n",
       "       180.81947276, 139.28855232, 166.34057824, 168.21740755,\n",
       "       166.85635308, 157.25564525, 153.07511802, 170.5047014 ,\n",
       "       148.89459078, 147.82606187, 169.59225024, 172.49937076,\n",
       "       150.15546775, 159.43678249, 151.28658494, 170.09882011,\n",
       "       139.28855232, 145.22796991, 168.86687696, 151.14106954,\n",
       "       173.8326125 , 147.67395737, 162.98882765, 164.19445284,\n",
       "       165.86949282, 158.29588882, 177.44982332, 176.63310381,\n",
       "       145.68176647, 144.89533817, 141.28855232, 159.31089703,\n",
       "       132.43084186, 158.97045618, 152.34044085, 144.68587711,\n",
       "       157.6259047 , 146.09348948, 144.81376843, 150.52199025,\n",
       "       139.78890677, 148.20460303, 172.26298175, 129.70102271,\n",
       "       146.22460666, 149.36689619, 168.11746636, 138.99677414,\n",
       "       161.17828507, 160.13906368, 128.72262003, 145.694571  ,\n",
       "       135.68275016, 165.6097754 , 170.08876884, 142.07834387,\n",
       "       147.85125875, 131.56667967, 153.35001951, 151.94250606,\n",
       "       157.07189216, 159.23924118, 137.56097534, 155.28855232,\n",
       "       149.00645172, 132.43084186, 143.43332033, 160.28359538,\n",
       "       133.64502355, 137.07338693, 143.93306478, 147.21791864,\n",
       "       158.86093632, 158.0002363 , 142.28929971, 132.85472091,\n",
       "       129.93507065, 130.64352877, 173.78581831, 169.64974419,\n",
       "       142.8624311 , 150.44649854, 165.1306061 , 155.57785203,\n",
       "       135.21740755, 144.92338721, 125.43556249, 144.30271423,\n",
       "       135.40723867, 148.29327296, 157.31612875, 168.56369012,\n",
       "       123.252282  , 139.04024356, 151.626141  , 149.23082207,\n",
       "       128.85298983, 136.28855232, 158.06642413, 142.35969709,\n",
       "       130.43953574, 153.32332786, 130.1913259 , 145.95021625,\n",
       "       154.48459885, 155.87271868, 135.25615634, 124.19604653,\n",
       "       166.54470866, 148.854622  , 155.21180213, 140.39681371,\n",
       "       148.27588519, 156.72808806, 144.08866993, 150.92797045,\n",
       "       148.37236422, 150.35969709, 159.15892992, 141.8323762 ,\n",
       "       141.05861503, 167.97466573, 140.34716735, 134.7864283 ,\n",
       "       155.10442553, 154.81984646, 141.84504333, 133.82840295,\n",
       "       149.99083351, 155.50306923, 115.42625861, 153.85710047,\n",
       "       138.86229371, 151.35969709, 128.34244671, 159.14626278,\n",
       "       136.51587375, 140.64427616, 132.27863844, 140.45145548,\n",
       "       141.54321388, 175.62960316, 143.41022823, 134.67619953,\n",
       "       154.01528299, 144.2978562 , 133.11036616, 147.44612485,\n",
       "       143.21740755, 123.59238763, 151.1842257 , 120.6183319 ,\n",
       "       125.71406354, 138.43219924, 147.00125847, 147.14500431,\n",
       "       168.86032633, 155.31052334, 144.59510241, 128.05576286,\n",
       "       147.47739975, 132.13967367, 135.25517265, 143.28196321,\n",
       "       134.34169932, 149.41369039, 144.49455121, 123.90426836,\n",
       "       151.77254118, 135.6337138 , 158.06058241, 127.71593201,\n",
       "       117.21740755, 130.10367814, 140.25550786, 158.31711244,\n",
       "       130.19333176, 132.412333  , 127.83920161, 165.48188408,\n",
       "       138.58318266, 120.13621151, 133.60566476, 142.9389065 ,\n",
       "       153.94923256, 132.01005127, 137.37185313, 147.73369348,\n",
       "       132.4511807 , 151.31290289, 121.51016942, 129.53031045,\n",
       "       147.5427028 , 157.82307232, 133.68078278, 152.21740755,\n",
       "       133.46734848, 138.43107816, 154.0002363 , 116.08048713,\n",
       "       132.24644028, 138.50198662, 158.6173482 , 142.37891485,\n",
       "       136.55204516, 126.85000027, 147.72710437, 171.28855232,\n",
       "       133.28855232, 130.28855232, 124.54246649, 148.58294636,\n",
       "       139.18212093, 133.00859497, 155.03463815, 118.48235669,\n",
       "       145.8785604 , 142.08493298, 144.74608583, 151.11559788,\n",
       "       164.85550678, 132.42014211, 133.21740755, 149.21740755,\n",
       "       151.67008303, 137.95147472, 126.45267548, 146.79539697,\n",
       "       122.35086581, 146.67096781, 136.79937022, 136.00893019,\n",
       "       134.87944518, 164.44761963, 136.36465403, 150.31622766,\n",
       "       132.90613683, 135.09586904, 130.24209334, 138.57313139,\n",
       "       134.57313139, 136.62155776, 125.52470503, 161.47926822,\n",
       "       130.8320025 , 131.17197074, 129.5800942 , 158.15322559,\n",
       "       139.52470503, 150.18772634, 122.69667577, 126.88042887,\n",
       "       133.35188798, 133.8092841 , 133.70448488, 142.44575116,\n",
       "       151.11260833, 158.7567471 , 139.00013739, 138.66315871,\n",
       "       159.96620815, 143.32631741, 119.92871784, 147.84704919,\n",
       "       142.60253782, 131.54769821, 133.45627503, 151.09797381,\n",
       "       148.03351706, 155.76907901, 162.04828897, 165.95965752,\n",
       "       160.2604648 , 162.20389412, 144.74496474, 152.77053531,\n",
       "       145.36492881, 140.7557634 , 143.09512165, 158.14626278,\n",
       "       174.4200432 , 164.85930415, 132.50198662, 171.50198662,\n",
       "       141.4803893 , 131.52358394, 150.88963384, 170.09624273,\n",
       "       143.92932784, 158.52708459, 146.75464232, 143.25330417,\n",
       "       165.32380047, 151.73084131, 156.67952431, 140.66190023,\n",
       "       152.08194343, 184.54582974, 166.00397325, 133.15991361,\n",
       "       149.89183753, 130.62094776, 140.92600307, 149.78371353,\n",
       "       154.50881204, 165.65110157, 170.06544044, 148.66363131,\n",
       "       176.05005853, 145.66933564, 146.19632132, 182.86453588,\n",
       "       163.00682541, 152.74109041, 119.40914562, 147.50769095,\n",
       "       170.94199497, 159.13424413, 148.00570433, 138.28682124,\n",
       "       163.57832464, 156.49842555, 161.21923755, 161.71817418,\n",
       "       168.71715201, 173.71266768, 140.57088923, 162.22525514,\n",
       "       165.64091291, 138.36927575, 156.64203399, 157.50820204,\n",
       "       144.20385564, 132.43705727, 124.72051526, 155.13495304,\n",
       "       154.01018866, 142.28855232, 147.17285552, 145.55672731,\n",
       "       182.15645145, 143.13607412, 142.83223881, 149.7103266 ,\n",
       "       141.45519243, 162.98981134, 137.21309909, 158.12734176,\n",
       "       157.74738278, 170.72991805, 157.09834751, 167.66896195,\n",
       "       162.33898455, 168.84354855, 157.1371952 , 164.15533036,\n",
       "       166.27948474, 162.83050773, 140.87981887, 133.98186484,\n",
       "       169.34183671, 174.05300961, 144.21740755, 151.11111355,\n",
       "       160.83560206, 155.70550705, 156.75141646, 147.85176983,\n",
       "       165.12728132, 157.31347441, 176.52610089, 154.13518934,\n",
       "       155.90277358, 166.10320553, 188.48693993, 143.40078695,\n",
       "       177.45584091, 147.55411145, 183.5190007 , 154.22446927,\n",
       "       159.11223463, 145.27840214, 188.91890287, 147.2514357 ,\n",
       "       168.09213209, 160.73552347, 179.82456709, 162.89482708,\n",
       "       157.13013349, 155.70326488, 152.79783696, 153.43084186,\n",
       "       172.04197465, 163.9918172 , 174.93218   , 177.66526348,\n",
       "       147.58593592, 171.31772244, 169.60230151, 171.50198662,\n",
       "       172.39325263, 150.27177455, 149.62328884, 162.81152627,\n",
       "       166.66923673, 168.76557837, 167.90786791, 196.28855232,\n",
       "       164.45930307, 171.02169623, 163.4345788 , 152.37718377,\n",
       "       180.9153418 , 149.71168398, 143.10007859, 153.28855232,\n",
       "       163.88664428, 153.21740755, 166.76160513, 175.82718296,\n",
       "       158.06343458, 190.94451192, 173.14785647, 156.05946133,\n",
       "       156.48632994, 159.6177219 , 182.04856376, 158.21740755,\n",
       "       168.42789078, 206.00692432, 157.42153797, 175.48235669,\n",
       "       179.26563609, 179.70009945, 175.947463  , 188.74832799,\n",
       "       174.25564525, 186.50198662, 186.02849969, 183.70780965,\n",
       "       171.11335572, 175.43617249, 197.46239154, 169.68115647,\n",
       "       163.50867465, 185.68251386, 175.07104586, 185.29126709,\n",
       "       192.18188462, 174.11189942, 181.32417416, 188.64301769,\n",
       "       183.47033803, 175.21479169, 149.97232466, 183.60477998,\n",
       "       184.28855232, 181.28855232, 176.28855232, 207.59948783,\n",
       "       210.18711635, 183.87220759, 187.79712805, 188.00397325,\n",
       "       178.37025945, 199.1976787 , 175.40312803, 172.9842444 ,\n",
       "       174.36886358, 181.21740755, 170.7550545 , 200.36489033,\n",
       "       174.37423269, 170.71817418, 179.23316315, 174.85893046,\n",
       "       204.87621932, 208.10816247, 209.91937548, 181.96587294,\n",
       "       193.17258074, 203.28855232, 202.19108959, 201.21740755,\n",
       "       175.77889398, 196.79053894, 175.99299872, 172.854012  ,\n",
       "       169.67059412, 179.19478806, 189.57340617, 183.62179406,\n",
       "       197.8618211 , 183.07511802, 197.17271813, 185.13393087,\n",
       "       192.19917348, 181.95928383, 171.14626278, 190.91856766,\n",
       "       217.99588936, 208.72747806, 177.46347414, 198.25989328,\n",
       "       177.84962658, 181.67236368, 187.08774667, 190.21740755,\n",
       "       181.53404739, 217.91913918, 221.9304874 , 215.64661724,\n",
       "       188.28855232, 188.37970072, 189.44921332, 192.59150285,\n",
       "       180.76656207, 168.84331225, 203.72499959, 176.26457544,\n",
       "       193.27891322, 196.80975671, 196.23095946, 213.76258882,\n",
       "       175.44608637, 195.41957059, 235.02149841, 207.63802227,\n",
       "       185.98957503, 182.96475185, 183.58752961, 232.72981914,\n",
       "       187.5268098 , 199.08554298, 188.46093524, 195.2688839 ,\n",
       "       186.23062425, 171.73164913, 207.89112862, 180.69640099,\n",
       "       207.06469305, 198.3974237 , 198.40354021, 180.79020373,\n",
       "       210.2447092 , 200.77366226, 217.64708985, 192.70896921,\n",
       "       173.56667967, 171.42439014, 191.4114867 , 197.03623184,\n",
       "       214.04931114, 211.92637676, 210.12442916, 221.84630181,\n",
       "       194.40336434, 204.37507899, 198.81087779, 194.88817754,\n",
       "       188.99256459, 211.02679056, 209.56817445, 184.58454005,\n",
       "       195.90407053, 206.51339528, 193.38200332])

      Variances (σ²):
      array([194.50198662, 214.93282848, 215.43084186, 222.58902438,\n",
       "       194.98808026, 184.28855232, 210.28457907, 205.01589299,\n",
       "       203.80245869, 221.93680173, 224.27265933, 219.71144768,\n",
       "       205.2213808 , 226.06319827, 216.28855232, 224.28855232,\n",
       "       227.82232493, 213.19356806, 192.14626278, 230.03178598,\n",
       "       259.74285997, 197.6641424 , 221.05562547, 203.82195123,\n",
       "       227.33983085, 200.36764358, 200.14626278, 199.13831629,\n",
       "       197.58542483, 215.17767507, 212.09461056, 222.07549171,\n",
       "       197.19716762, 194.43084186, 196.47057434, 241.77899289,\n",
       "       229.24881984, 198.00397325, 211.45468134, 199.01589299,\n",
       "       220.37521639, 209.24162073, 200.16935488, 188.81003149,\n",
       "       211.96858771, 184.72374112, 201.24522029, 209.52507872,\n",
       "       192.48646733, 227.57387878, 196.97541312, 215.29799359,\n",
       "       234.19642023, 205.3203383 , 193.86168371, 230.81400474,\n",
       "       236.18562157, 217.84939027, 192.57313139, 216.96101491,\n",
       "       198.46660109, 227.7710464 , 197.75477971, 192.48609363,\n",
       "       222.37956333, 215.5417191 , 224.25714003, 214.51315898,\n",
       "       234.05525178, 211.35647123, 217.02706535, 194.02383949,\n",
       "       214.72014156, 230.51825331, 203.65110157, 216.78296614,\n",
       "       208.61571604, 204.96536185, 213.27228563, 234.07511802,\n",
       "       197.02421318, 208.62440992, 173.76669946, 198.95666797,\n",
       "       217.53054675, 206.31127072, 219.36404403, 192.19791501,\n",
       "       224.24087334, 194.94362714, 214.59846566, 211.34268301,\n",
       "       204.84219117, 197.37086944, 219.45977567, 191.96784032,\n",
       "       197.53699847, 200.16178208, 211.88304473, 214.69158144,\n",
       "       211.46262784, 215.86093632, 226.58865069, 197.27700627,\n",
       "       229.28383168, 207.04047987, 210.24919354, 225.15780883,\n",
       "       249.19716762, 231.01664038, 222.17295444, 215.37633747,\n",
       "       206.74360736, 193.20076717, 209.73056653, 199.03575923,\n",
       "       224.21740755, 238.24199443, 216.3279111 , 218.93282848,\n",
       "       204.71579462, 216.47020064, 200.61894189, 205.10442553,\n",
       "       214.96858771, 227.61174279, 209.5800942 , 232.60590107,\n",
       "       230.18674265, 203.71467354, 229.72411481, 214.28383168,\n",
       "       212.23975226, 204.70139641, 200.27985844, 216.14626278,\n",
       "       220.893096  , 222.12751763, 221.63992922, 218.84504333,\n",
       "       240.16290316, 222.79885913, 214.34057824, 204.85186874,\n",
       "       215.34057824, 218.30084576, 222.09970489, 214.2365264 ,\n",
       "       211.85883155, 216.52110547, 233.66227393, 176.21852864,\n",
       "       219.70312749, 202.98770656, 205.35460276, 230.27737997,\n",
       "       224.00186847, 212.79053894, 209.72659328, 200.38713613,\n",
       "       214.52942566, 214.56990553, 215.14303692, 219.38713613,\n",
       "       206.05165222, 199.14948864, 203.81400474, 203.05165222,\n",
       "       220.8624311 , 231.00322586, 176.25887111, 192.02818643,\n",
       "       244.5114279 , 201.55363885, 216.21268692, 190.07511802,\n",
       "       230.85298983, 213.95232103, 196.82952403, 198.85957894,\n",
       "       186.88117626, 189.33337913, 226.73888672, 206.88117626,\n",
       "       204.56742706, 201.75837927, 214.56331642, 196.80891041,\n",
       "       203.94946886, 212.07511802, 193.89781664, 193.93282848,\n",
       "       190.12564916, 184.34417779, 195.6948073 , 206.23292685,\n",
       "       182.3131392 , 205.23292685, 216.74398105, 219.85920524,\n",
       "       230.14874125, 198.63273011, 253.37521639, 205.08516929,\n",
       "       189.84616441, 205.08058604, 198.4579072 , 213.74248627,\n",
       "       233.81152627, 187.36156556, 225.61174279, 189.87322976,\n",
       "       192.86168371, 213.41172301, 192.38452027, 192.35969709,\n",
       "       200.48286777, 205.76520468, 203.66961043, 180.91370963,\n",
       "       217.46014937, 188.48422516, 198.95456319, 183.33796238,\n",
       "       194.12354438, 212.09026362, 204.07511802, 206.19108959,\n",
       "       200.49838707, 204.98906395, 189.18749004, 209.99280089,\n",
       "       216.86168371, 219.10106228, 172.374469  , 175.15869361,\n",
       "       224.81250996, 214.321322  , 208.39844587, 192.08951623,\n",
       "       193.80456346, 204.77974028, 213.35287168, 194.92600307,\n",
       "       207.65147527, 205.44734485, 201.57313139, 180.27069194,\n",
       "       208.01005127, 192.50769095, 212.21740755, 204.93259218,\n",
       "       198.50110184, 198.28929971, 228.57897311, 201.57201031,\n",
       "       205.86205741, 193.14116845, 200.72958284, 168.8475218 ,\n",
       "       221.44463007, 167.49254535, 208.0841856 , 175.69728577,\n",
       "       191.16402425, 206.30159315, 169.86168371, 212.92848154,\n",
       "       187.93282848, 189.76917793, 220.37671116, 179.05810394,\n",
       "       182.10479922, 191.83634945, 198.32220678, 178.69840685,\n",
       "       182.41817472, 197.69854424, 205.74459105, 208.39297785,\n",
       "       193.77389856, 184.12130221, 202.34244671, 178.5654212 ,\n",
       "       203.7021438 , 196.32244308, 169.12466546, 199.82045645,\n",
       "       193.58542483, 209.17916985, 172.97106618, 207.30084576,\n",
       "       184.54022433, 172.21740755, 179.22970099, 167.32281677,\n",
       "       201.18973221, 185.7250985 , 188.71368985, 196.74968538,\n",
       "       199.03080229, 187.77773442, 198.54681343, 164.77203009,\n",
       "       162.13845368, 209.33113696, 199.61571604, 167.81288365,\n",
       "       188.40241912, 207.95517319, 189.1299961 , 187.86378849,\n",
       "       158.82307232, 183.55686471, 172.64614463, 182.93096001,\n",
       "       189.28855232, 168.53662477, 161.34529887, 194.51638484,\n",
       "       175.91843027, 178.76571576, 166.60837954, 187.37993702,\n",
       "       182.74944908, 178.58355636, 181.48174669, 175.61758451,\n",
       "       191.48660472, 172.32814741, 171.62403623, 172.76632576,\n",
       "       170.83499206, 186.55041299, 189.43332033, 205.52769458,\n",
       "       167.83747053, 179.16748641, 165.87893409, 168.69817055,\n",
       "       165.8092841 , 173.44724593, 169.60912693, 168.76235251,\n",
       "       177.35969709, 189.43097925, 180.22671143, 173.38788352,\n",
       "       163.88178625, 178.47380019, 165.75702188, 161.27378041,\n",
       "       162.98920134, 178.26036589, 169.65904807, 179.30594009,\n",
       "       159.27116455, 164.81872537, 186.97157727, 162.62688839,\n",
       "       181.37049575, 178.91829287, 177.89758033, 161.62291514,\n",
       "       159.14626278, 171.26321805, 188.13943737, 178.67956278,\n",
       "       152.26976868, 175.70289119, 151.21740755, 148.81984646,\n",
       "       167.76296251, 160.24498398, 183.00397325, 175.25527156,\n",
       "       163.28794232, 172.14687278, 155.18127463, 176.00845758,\n",
       "       161.35633384, 170.14626278, 167.64315508, 170.32977958,\n",
       "       171.07511802, 155.53699847, 160.61323757, 173.47094803,\n",
       "       162.97802898, 172.00159369, 146.66750565, 172.9237609 ,\n",
       "       154.32865849, 183.27948474, 148.84878028, 166.83064512,\n",
       "       169.15519297, 165.91978765, 196.27055455, 173.13322196,\n",
       "       180.24348921, 175.67629845, 165.14527909, 163.11620788,\n",
       "       161.17631769, 171.14626278, 187.01402452, 179.41382778,\n",
       "       164.91272594, 178.2444729 , 143.66526348, 175.25055092,\n",
       "       145.52297395, 147.23839487, 159.10428814, 144.67344628,\n",
       "       170.56892184, 155.9307237 , 169.36180186, 150.62305253,\n",
       "       174.62142036, 158.30977595, 175.59622349, 155.24236812,\n",
       "       140.24820984, 176.99813153, 151.01565669, 149.21156583,\n",
       "       168.07511802, 159.66736826, 169.49268274, 155.73800194,\n",
       "       180.81947276, 139.28855232, 166.34057824, 168.21740755,\n",
       "       166.85635308, 157.25564525, 153.07511802, 170.5047014 ,\n",
       "       148.89459078, 147.82606187, 169.59225024, 172.49937076,\n",
       "       150.15546775, 159.43678249, 151.28658494, 170.09882011,\n",
       "       139.28855232, 145.22796991, 168.86687696, 151.14106954,\n",
       "       173.8326125 , 147.67395737, 162.98882765, 164.19445284,\n",
       "       165.86949282, 158.29588882, 177.44982332, 176.63310381,\n",
       "       145.68176647, 144.89533817, 141.28855232, 159.31089703,\n",
       "       132.43084186, 158.97045618, 152.34044085, 144.68587711,\n",
       "       157.6259047 , 146.09348948, 144.81376843, 150.52199025,\n",
       "       139.78890677, 148.20460303, 172.26298175, 129.70102271,\n",
       "       146.22460666, 149.36689619, 168.11746636, 138.99677414,\n",
       "       161.17828507, 160.13906368, 128.72262003, 145.694571  ,\n",
       "       135.68275016, 165.6097754 , 170.08876884, 142.07834387,\n",
       "       147.85125875, 131.56667967, 153.35001951, 151.94250606,\n",
       "       157.07189216, 159.23924118, 137.56097534, 155.28855232,\n",
       "       149.00645172, 132.43084186, 143.43332033, 160.28359538,\n",
       "       133.64502355, 137.07338693, 143.93306478, 147.21791864,\n",
       "       158.86093632, 158.0002363 , 142.28929971, 132.85472091,\n",
       "       129.93507065, 130.64352877, 173.78581831, 169.64974419,\n",
       "       142.8624311 , 150.44649854, 165.1306061 , 155.57785203,\n",
       "       135.21740755, 144.92338721, 125.43556249, 144.30271423,\n",
       "       135.40723867, 148.29327296, 157.31612875, 168.56369012,\n",
       "       123.252282  , 139.04024356, 151.626141  , 149.23082207,\n",
       "       128.85298983, 136.28855232, 158.06642413, 142.35969709,\n",
       "       130.43953574, 153.32332786, 130.1913259 , 145.95021625,\n",
       "       154.48459885, 155.87271868, 135.25615634, 124.19604653,\n",
       "       166.54470866, 148.854622  , 155.21180213, 140.39681371,\n",
       "       148.27588519, 156.72808806, 144.08866993, 150.92797045,\n",
       "       148.37236422, 150.35969709, 159.15892992, 141.8323762 ,\n",
       "       141.05861503, 167.97466573, 140.34716735, 134.7864283 ,\n",
       "       155.10442553, 154.81984646, 141.84504333, 133.82840295,\n",
       "       149.99083351, 155.50306923, 115.42625861, 153.85710047,\n",
       "       138.86229371, 151.35969709, 128.34244671, 159.14626278,\n",
       "       136.51587375, 140.64427616, 132.27863844, 140.45145548,\n",
       "       141.54321388, 175.62960316, 143.41022823, 134.67619953,\n",
       "       154.01528299, 144.2978562 , 133.11036616, 147.44612485,\n",
       "       143.21740755, 123.59238763, 151.1842257 , 120.6183319 ,\n",
       "       125.71406354, 138.43219924, 147.00125847, 147.14500431,\n",
       "       168.86032633, 155.31052334, 144.59510241, 128.05576286,\n",
       "       147.47739975, 132.13967367, 135.25517265, 143.28196321,\n",
       "       134.34169932, 149.41369039, 144.49455121, 123.90426836,\n",
       "       151.77254118, 135.6337138 , 158.06058241, 127.71593201,\n",
       "       117.21740755, 130.10367814, 140.25550786, 158.31711244,\n",
       "       130.19333176, 132.412333  , 127.83920161, 165.48188408,\n",
       "       138.58318266, 120.13621151, 133.60566476, 142.9389065 ,\n",
       "       153.94923256, 132.01005127, 137.37185313, 147.73369348,\n",
       "       132.4511807 , 151.31290289, 121.51016942, 129.53031045,\n",
       "       147.5427028 , 157.82307232, 133.68078278, 152.21740755,\n",
       "       133.46734848, 138.43107816, 154.0002363 , 116.08048713,\n",
       "       132.24644028, 138.50198662, 158.6173482 , 142.37891485,\n",
       "       136.55204516, 126.85000027, 147.72710437, 171.28855232,\n",
       "       133.28855232, 130.28855232, 124.54246649, 148.58294636,\n",
       "       139.18212093, 133.00859497, 155.03463815, 118.48235669,\n",
       "       145.8785604 , 142.08493298, 144.74608583, 151.11559788,\n",
       "       164.85550678, 132.42014211, 133.21740755, 149.21740755,\n",
       "       151.67008303, 137.95147472, 126.45267548, 146.79539697,\n",
       "       122.35086581, 146.67096781, 136.79937022, 136.00893019,\n",
       "       134.87944518, 164.44761963, 136.36465403, 150.31622766,\n",
       "       132.90613683, 135.09586904, 130.24209334, 138.57313139,\n",
       "       134.57313139, 136.62155776, 125.52470503, 161.47926822,\n",
       "       130.8320025 , 131.17197074, 129.5800942 , 158.15322559,\n",
       "       139.52470503, 150.18772634, 122.69667577, 126.88042887,\n",
       "       133.35188798, 133.8092841 , 133.70448488, 142.44575116,\n",
       "       151.11260833, 158.7567471 , 139.00013739, 138.66315871,\n",
       "       159.96620815, 143.32631741, 119.92871784, 147.84704919,\n",
       "       142.60253782, 131.54769821, 133.45627503, 151.09797381,\n",
       "       148.03351706, 155.76907901, 162.04828897, 165.95965752,\n",
       "       160.2604648 , 162.20389412, 144.74496474, 152.77053531,\n",
       "       145.36492881, 140.7557634 , 143.09512165, 158.14626278,\n",
       "       174.4200432 , 164.85930415, 132.50198662, 171.50198662,\n",
       "       141.4803893 , 131.52358394, 150.88963384, 170.09624273,\n",
       "       143.92932784, 158.52708459, 146.75464232, 143.25330417,\n",
       "       165.32380047, 151.73084131, 156.67952431, 140.66190023,\n",
       "       152.08194343, 184.54582974, 166.00397325, 133.15991361,\n",
       "       149.89183753, 130.62094776, 140.92600307, 149.78371353,\n",
       "       154.50881204, 165.65110157, 170.06544044, 148.66363131,\n",
       "       176.05005853, 145.66933564, 146.19632132, 182.86453588,\n",
       "       163.00682541, 152.74109041, 119.40914562, 147.50769095,\n",
       "       170.94199497, 159.13424413, 148.00570433, 138.28682124,\n",
       "       163.57832464, 156.49842555, 161.21923755, 161.71817418,\n",
       "       168.71715201, 173.71266768, 140.57088923, 162.22525514,\n",
       "       165.64091291, 138.36927575, 156.64203399, 157.50820204,\n",
       "       144.20385564, 132.43705727, 124.72051526, 155.13495304,\n",
       "       154.01018866, 142.28855232, 147.17285552, 145.55672731,\n",
       "       182.15645145, 143.13607412, 142.83223881, 149.7103266 ,\n",
       "       141.45519243, 162.98981134, 137.21309909, 158.12734176,\n",
       "       157.74738278, 170.72991805, 157.09834751, 167.66896195,\n",
       "       162.33898455, 168.84354855, 157.1371952 , 164.15533036,\n",
       "       166.27948474, 162.83050773, 140.87981887, 133.98186484,\n",
       "       169.34183671, 174.05300961, 144.21740755, 151.11111355,\n",
       "       160.83560206, 155.70550705, 156.75141646, 147.85176983,\n",
       "       165.12728132, 157.31347441, 176.52610089, 154.13518934,\n",
       "       155.90277358, 166.10320553, 188.48693993, 143.40078695,\n",
       "       177.45584091, 147.55411145, 183.5190007 , 154.22446927,\n",
       "       159.11223463, 145.27840214, 188.91890287, 147.2514357 ,\n",
       "       168.09213209, 160.73552347, 179.82456709, 162.89482708,\n",
       "       157.13013349, 155.70326488, 152.79783696, 153.43084186,\n",
       "       172.04197465, 163.9918172 , 174.93218   , 177.66526348,\n",
       "       147.58593592, 171.31772244, 169.60230151, 171.50198662,\n",
       "       172.39325263, 150.27177455, 149.62328884, 162.81152627,\n",
       "       166.66923673, 168.76557837, 167.90786791, 196.28855232,\n",
       "       164.45930307, 171.02169623, 163.4345788 , 152.37718377,\n",
       "       180.9153418 , 149.71168398, 143.10007859, 153.28855232,\n",
       "       163.88664428, 153.21740755, 166.76160513, 175.82718296,\n",
       "       158.06343458, 190.94451192, 173.14785647, 156.05946133,\n",
       "       156.48632994, 159.6177219 , 182.04856376, 158.21740755,\n",
       "       168.42789078, 206.00692432, 157.42153797, 175.48235669,\n",
       "       179.26563609, 179.70009945, 175.947463  , 188.74832799,\n",
       "       174.25564525, 186.50198662, 186.02849969, 183.70780965,\n",
       "       171.11335572, 175.43617249, 197.46239154, 169.68115647,\n",
       "       163.50867465, 185.68251386, 175.07104586, 185.29126709,\n",
       "       192.18188462, 174.11189942, 181.32417416, 188.64301769,\n",
       "       183.47033803, 175.21479169, 149.97232466, 183.60477998,\n",
       "       184.28855232, 181.28855232, 176.28855232, 207.59948783,\n",
       "       210.18711635, 183.87220759, 187.79712805, 188.00397325,\n",
       "       178.37025945, 199.1976787 , 175.40312803, 172.9842444 ,\n",
       "       174.36886358, 181.21740755, 170.7550545 , 200.36489033,\n",
       "       174.37423269, 170.71817418, 179.23316315, 174.85893046,\n",
       "       204.87621932, 208.10816247, 209.91937548, 181.96587294,\n",
       "       193.17258074, 203.28855232, 202.19108959, 201.21740755,\n",
       "       175.77889398, 196.79053894, 175.99299872, 172.854012  ,\n",
       "       169.67059412, 179.19478806, 189.57340617, 183.62179406,\n",
       "       197.8618211 , 183.07511802, 197.17271813, 185.13393087,\n",
       "       192.19917348, 181.95928383, 171.14626278, 190.91856766,\n",
       "       217.99588936, 208.72747806, 177.46347414, 198.25989328,\n",
       "       177.84962658, 181.67236368, 187.08774667, 190.21740755,\n",
       "       181.53404739, 217.91913918, 221.9304874 , 215.64661724,\n",
       "       188.28855232, 188.37970072, 189.44921332, 192.59150285,\n",
       "       180.76656207, 168.84331225, 203.72499959, 176.26457544,\n",
       "       193.27891322, 196.80975671, 196.23095946, 213.76258882,\n",
       "       175.44608637, 195.41957059, 235.02149841, 207.63802227,\n",
       "       185.98957503, 182.96475185, 183.58752961, 232.72981914,\n",
       "       187.5268098 , 199.08554298, 188.46093524, 195.2688839 ,\n",
       "       186.23062425, 171.73164913, 207.89112862, 180.69640099,\n",
       "       207.06469305, 198.3974237 , 198.40354021, 180.79020373,\n",
       "       210.2447092 , 200.77366226, 217.64708985, 192.70896921,\n",
       "       173.56667967, 171.42439014, 191.4114867 , 197.03623184,\n",
       "       214.04931114, 211.92637676, 210.12442916, 221.84630181,\n",
       "       194.40336434, 204.37507899, 198.81087779, 194.88817754,\n",
       "       188.99256459, 211.02679056, 209.56817445, 184.58454005,\n",
       "       195.90407053, 206.51339528, 193.38200332])
    • Ltotal
      ()
      float64
      m
      59.0
      Values:
      array(59.)
    • instrument_name
      ()
      string
      POWGEN
      Values:
      'POWGEN'
    • position
      ()
      vector3
      m
      [ 0. 0. -1.]
      Values:
      array([ 0.,  0., -1.])
    • sample_position
      ()
      vector3
      m
      [0. 0. 0.]
      Values:
      array([0., 0., 0.])
    • source_position
      ()
      vector3
      m
      [ 0. 0. -60.]
      Values:
      array([  0.,   0., -60.])
    • spectrum
      ()
      int32
      1
      Values:
      array(1, dtype=int32)
" ], "text/plain": [ "\n", "Dimensions: Sizes[wavelength:999, ]\n", "Coordinates:\n", " wavelength float64 [Å] (wavelength [bin-edge]) [0, 0.001001, ..., 0.998999, 1]\n", "Data:\n", " float64 [counts] (wavelength) [194.502, 214.933, ..., 206.513, 193.382] [194.502, 214.933, ..., 206.513, 193.382]\n", "Attributes:\n", " Ltotal float64 [m] () [59]\n", " instrument_name string [dimensionless] () [\"POWGEN\"]\n", " position vector3 [m] () [(0, 0, -1)]\n", " sample_position vector3 [m] () [(0, 0, 0)]\n", " source_position vector3 [m] () [(0, 0, -60)]\n", " spectrum int32 [dimensionless] () [1]\n" ] }, "execution_count": 13, "metadata": {}, "output_type": "execute_result" } ], "source": [ "edges = sc.Variable(dims=['wavelength'], unit=sc.units.angstrom, values=np.linspace(0, 1, num=1000))\n", "mon = sc.rebin(mon, 'wavelength', edges)\n", "mon" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "We intend to normalize each event to the relative monitor counts (compared to the total monitor counts).\n", "We use `sum` to compute the total monitor counts and obtain the relative counts using division:" ] }, { "cell_type": "code", "execution_count": 14, "metadata": {}, "outputs": [], "source": [ "mon /= mon.sum('wavelength')" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "The sample data is *event data in bins* and the monitor is a histogram.\n", "Multiplication and division operations for such cases are supported by modifying the weights (values) for each event using the operators of the `bins` property, in combination with the `sc.lookup` helper, a wrapper for a discrete \"function\", given by the monitor:" ] }, { "cell_type": "code", "execution_count": 15, "metadata": {}, "outputs": [ { "data": { "text/html": [ "
\n", "\n", "\n", "\n", "\n", "Show/Hide data repr\n", "\n", "\n", "\n", "\n", "\n", "Show/Hide attributes\n", "\n", "\n", "\n", "\n", "\n", "\n", "
scipp.DataArray (412.87 MB)
    • spectrum: 24794
    • wavelength: 1
    • spectrum
      (spectrum)
      int32
      1, 2, ..., 24793, 24794
      Values:
      array([    1,     2,     3, ..., 24792, 24793, 24794], dtype=int32)
    • two_theta
      (spectrum)
      float64
      rad
      2.491, 2.504, ..., 0.442, 0.442
      Values:
      array([2.49144445, 2.50372965, 2.51564288, ..., 0.44118918, 0.44157918,\n",
       "       0.44231324])
    • wavelength
      (spectrum, wavelength [bin-edge])
      float64
      Å
      0.001, 1.056, ..., 0.001, 1.028
      Values:
      array([[0.00120148, 1.05567339],\n",
       "       [0.00120156, 1.05574713],\n",
       "       [0.00120162, 1.05580152],\n",
       "       ...,\n",
       "       [0.00116986, 1.02789183],\n",
       "       [0.00116986, 1.02789145],\n",
       "       [0.00116984, 1.02787999]])
    • (spectrum, wavelength)
      DataArrayView
      binned data [len=0, len=0, ..., len=0, len=0]
      Values:
      [<scipp.DataArray>\n",
       "Dimensions: Sizes[event:0, ]\n",
       "Coordinates:\n",
       "  wavelength                float64             [Å]  (event)  []\n",
       "Data:\n",
       "                            float32         [counts]  (event)  []  []\n",
       "Attributes:\n",
       "  tof                       float64            [µs]  (event)  []\n",
       "\n",
       ", <scipp.DataArray>\n",
       "Dimensions: Sizes[event:0, ]\n",
       "Coordinates:\n",
       "  wavelength                float64             [Å]  (event)  []\n",
       "Data:\n",
       "                            float32         [counts]  (event)  []  []\n",
       "Attributes:\n",
       "  tof                       float64            [µs]  (event)  []\n",
       "\n",
       ", ..., <scipp.DataArray>\n",
       "Dimensions: Sizes[event:0, ]\n",
       "Coordinates:\n",
       "  wavelength                float64             [Å]  (event)  []\n",
       "Data:\n",
       "                            float32         [counts]  (event)  []  []\n",
       "Attributes:\n",
       "  tof                       float64            [µs]  (event)  []\n",
       "\n",
       ", <scipp.DataArray>\n",
       "Dimensions: Sizes[event:0, ]\n",
       "Coordinates:\n",
       "  wavelength                float64             [Å]  (event)  []\n",
       "Data:\n",
       "                            float32         [counts]  (event)  []  []\n",
       "Attributes:\n",
       "  tof                       float64            [µs]  (event)  []\n",
       "\n",
       "]
    • ChopperStatus1
      ()
      DataArray
      {dims=[time: 2], coords=[time]}
      Values:
      <scipp.DataArray>\n",
       "Dimensions: Sizes[time:2, ]\n",
       "Coordinates:\n",
       "  time                    datetime64             [ns]  (time)  [2011-08-12T15:50:16.800000000, 2011-08-12T17:22:05.000000000]\n",
       "Data:\n",
       "                            float64  [dimensionless]  (time)  [4, 4]\n",
       "\n",
       "
    • ChopperStatus2
      ()
      DataArray
      {dims=[time: 2], coords=[time]}
      Values:
      <scipp.DataArray>\n",
       "Dimensions: Sizes[time:2, ]\n",
       "Coordinates:\n",
       "  time                    datetime64             [ns]  (time)  [2011-08-12T15:50:16.816000000, 2011-08-12T17:22:05.000000000]\n",
       "Data:\n",
       "                            float64  [dimensionless]  (time)  [4, 4]\n",
       "\n",
       "
    • ChopperStatus3
      ()
      DataArray
      {dims=[time: 2], coords=[time]}
      Values:
      <scipp.DataArray>\n",
       "Dimensions: Sizes[time:2, ]\n",
       "Coordinates:\n",
       "  time                    datetime64             [ns]  (time)  [2011-08-12T15:50:16.816000000, 2011-08-12T17:22:05.000000000]\n",
       "Data:\n",
       "                            float64  [dimensionless]  (time)  [4, 4]\n",
       "\n",
       "
    • CurrentSP
      ()
      DataArray
      {dims=[time: 2], coords=[time]}
      Values:
      <scipp.DataArray>\n",
       "Dimensions: Sizes[time:2, ]\n",
       "Coordinates:\n",
       "  time                    datetime64             [ns]  (time)  [2011-08-12T15:50:16.800000000, 2011-08-12T17:22:05.000000000]\n",
       "Data:\n",
       "                            float64  [dimensionless]  (time)  [300, 300]\n",
       "\n",
       "
    • EnergyRequest
      ()
      DataArray
      {dims=[time: 2], unit=meV, coords=[time]}
      Values:
      <scipp.DataArray>\n",
       "Dimensions: Sizes[time:2, ]\n",
       "Coordinates:\n",
       "  time                    datetime64             [ns]  (time)  [2011-08-12T15:50:16.784000000, 2011-08-12T17:22:05.000000000]\n",
       "Data:\n",
       "                            float64            [meV]  (time)  [287.955, 287.955]\n",
       "\n",
       "
    • Filename
      ()
      string
      /home/scientist/.cache/scippneutron/2/PG3_4844_event.nxs
      Values:
      '/home/scientist/.cache/scippneutron/2/PG3_4844_event.nxs'
    • L1
      ()
      float64
      m
      60.0
      Values:
      array(60.)
    • L2
      (spectrum)
      float64
      m
      2.560, 2.556, ..., 4.251, 4.252
      Values:
      array([2.56027902, 2.55590928, 2.55268677, ..., 4.25113991, 4.25116338,\n",
       "       4.25188   ])
    • LKSRampRate
      ()
      DataArray
      {dims=[time: 2], coords=[time]}
      Values:
      <scipp.DataArray>\n",
       "Dimensions: Sizes[time:2, ]\n",
       "Coordinates:\n",
       "  time                    datetime64             [ns]  (time)  [2011-08-12T15:50:16.800000000, 2011-08-12T17:22:05.000000000]\n",
       "Data:\n",
       "                            float64  [dimensionless]  (time)  [0, 0]\n",
       "\n",
       "
    • LambdaRequest
      ()
      DataArray
      {dims=[time: 2], unit=Å, coords=[time]}
      Values:
      <scipp.DataArray>\n",
       "Dimensions: Sizes[time:2, ]\n",
       "Coordinates:\n",
       "  time                    datetime64             [ns]  (time)  [2011-08-12T15:50:16.784000000, 2011-08-12T17:22:05.000000000]\n",
       "Data:\n",
       "                            float64             [Å]  (time)  [0.533, 0.533]\n",
       "\n",
       "
    • Ltotal
      (spectrum)
      float64
      m
      62.560, 62.556, ..., 64.251, 64.252
      Values:
      array([62.56027902, 62.55590928, 62.55268677, ..., 64.25113991,\n",
       "       64.25116338, 64.25188   ])
    • Phase1
      ()
      DataArray
      {dims=[time: 1794], unit=µs, coords=[time]}
      Values:
      <scipp.DataArray>\n",
       "Dimensions: Sizes[time:1794, ]\n",
       "Coordinates:\n",
       "  time                    datetime64             [ns]  (time)  [2011-08-12T15:50:16.800000000, 2011-08-12T15:50:18.267000007, ..., 2011-08-12T17:22:02.999218750, 2011-08-12T17:22:05.000000000]\n",
       "Data:\n",
       "                            float64            [µs]  (time)  [8166.72, 8165.16, ..., 8163.85, 8163.85]\n",
       "\n",
       "
    • Phase2
      ()
      DataArray
      {dims=[time: 1793], unit=µs, coords=[time]}
      Values:
      <scipp.DataArray>\n",
       "Dimensions: Sizes[time:1793, ]\n",
       "Coordinates:\n",
       "  time                    datetime64             [ns]  (time)  [2011-08-12T15:50:16.800000000, 2011-08-12T15:50:18.267000007, ..., 2011-08-12T17:22:02.999218750, 2011-08-12T17:22:05.000000000]\n",
       "Data:\n",
       "                            float64            [µs]  (time)  [8335.63, 8334.09, ..., 8332.86, 8332.86]\n",
       "\n",
       "
    • Phase3
      ()
      DataArray
      {dims=[time: 1777], unit=µs, coords=[time]}
      Values:
      <scipp.DataArray>\n",
       "Dimensions: Sizes[time:1777, ]\n",
       "Coordinates:\n",
       "  time                    datetime64             [ns]  (time)  [2011-08-12T15:50:16.800000000, 2011-08-12T15:50:19.688999938, ..., 2011-08-12T17:22:04.421093750, 2011-08-12T17:22:05.000000000]\n",
       "Data:\n",
       "                            float64            [µs]  (time)  [14001.5, 14001.4, ..., 14002, 14002]\n",
       "\n",
       "
    • PhaseRequest1
      ()
      DataArray
      {dims=[time: 2], unit=µs, coords=[time]}
      Values:
      <scipp.DataArray>\n",
       "Dimensions: Sizes[time:2, ]\n",
       "Coordinates:\n",
       "  time                    datetime64             [ns]  (time)  [2011-08-12T15:50:16.800000000, 2011-08-12T17:22:05.000000000]\n",
       "Data:\n",
       "                            float64            [µs]  (time)  [8164.07, 8164.07]\n",
       "\n",
       "
    • PhaseRequest2
      ()
      DataArray
      {dims=[time: 2], unit=µs, coords=[time]}
      Values:
      <scipp.DataArray>\n",
       "Dimensions: Sizes[time:2, ]\n",
       "Coordinates:\n",
       "  time                    datetime64             [ns]  (time)  [2011-08-12T15:50:16.800000000, 2011-08-12T17:22:05.000000000]\n",
       "Data:\n",
       "                            float64            [µs]  (time)  [8332.89, 8332.89]\n",
       "\n",
       "
    • PhaseRequest3
      ()
      DataArray
      {dims=[time: 2], unit=µs, coords=[time]}
      Values:
      <scipp.DataArray>\n",
       "Dimensions: Sizes[time:2, ]\n",
       "Coordinates:\n",
       "  time                    datetime64             [ns]  (time)  [2011-08-12T15:50:16.800000000, 2011-08-12T17:22:05.000000000]\n",
       "Data:\n",
       "                            float64            [µs]  (time)  [14001.8, 14001.8]\n",
       "\n",
       "
    • S1HCenter
      ()
      DataArray
      {dims=[time: 2], unit=mm, coords=[time]}
      Values:
      <scipp.DataArray>\n",
       "Dimensions: Sizes[time:2, ]\n",
       "Coordinates:\n",
       "  time                    datetime64             [ns]  (time)  [2011-08-12T15:50:16.800000000, 2011-08-12T17:22:05.000000000]\n",
       "Data:\n",
       "                            float64             [mm]  (time)  [0, 0]\n",
       "\n",
       "
    • S1HCenterOffset
      ()
      DataArray
      {dims=[time: 2], coords=[time]}
      Values:
      <scipp.DataArray>\n",
       "Dimensions: Sizes[time:2, ]\n",
       "Coordinates:\n",
       "  time                    datetime64             [ns]  (time)  [2011-08-12T15:50:16.800000000, 2011-08-12T17:22:05.000000000]\n",
       "Data:\n",
       "                            float64  [dimensionless]  (time)  [0, 0]\n",
       "\n",
       "
    • S1HWidth
      ()
      DataArray
      {dims=[time: 2], unit=mm, coords=[time]}
      Values:
      <scipp.DataArray>\n",
       "Dimensions: Sizes[time:2, ]\n",
       "Coordinates:\n",
       "  time                    datetime64             [ns]  (time)  [2011-08-12T15:50:16.800000000, 2011-08-12T17:22:05.000000000]\n",
       "Data:\n",
       "                            float64             [mm]  (time)  [10, 10]\n",
       "\n",
       "
    • S1VCenter
      ()
      DataArray
      {dims=[time: 2], unit=mm, coords=[time]}
      Values:
      <scipp.DataArray>\n",
       "Dimensions: Sizes[time:2, ]\n",
       "Coordinates:\n",
       "  time                    datetime64             [ns]  (time)  [2011-08-12T15:50:16.800000000, 2011-08-12T17:22:05.000000000]\n",
       "Data:\n",
       "                            float64             [mm]  (time)  [5, 5]\n",
       "\n",
       "
    • S1VCenterOffset
      ()
      DataArray
      {dims=[time: 2], coords=[time]}
      Values:
      <scipp.DataArray>\n",
       "Dimensions: Sizes[time:2, ]\n",
       "Coordinates:\n",
       "  time                    datetime64             [ns]  (time)  [2011-08-12T15:50:16.800000000, 2011-08-12T17:22:05.000000000]\n",
       "Data:\n",
       "                            float64  [dimensionless]  (time)  [0, 0]\n",
       "\n",
       "
    • S1VHeight
      ()
      DataArray
      {dims=[time: 2], unit=mm, coords=[time]}
      Values:
      <scipp.DataArray>\n",
       "Dimensions: Sizes[time:2, ]\n",
       "Coordinates:\n",
       "  time                    datetime64             [ns]  (time)  [2011-08-12T15:50:16.816000000, 2011-08-12T17:22:05.000000000]\n",
       "Data:\n",
       "                            float64             [mm]  (time)  [30, 30]\n",
       "\n",
       "
    • SampleTemp
      ()
      DataArray
      {dims=[time: 467], coords=[time]}
      Values:
      <scipp.DataArray>\n",
       "Dimensions: Sizes[time:467, ]\n",
       "Coordinates:\n",
       "  time                    datetime64             [ns]  (time)  [2011-08-12T15:50:16.784000000, 2011-08-12T15:50:18.001000007, ..., 2011-08-12T17:21:58.623843750, 2011-08-12T17:22:05.000000000]\n",
       "Data:\n",
       "                            float64  [dimensionless]  (time)  [299.352, 299.446, ..., 300, 300]\n",
       "\n",
       "
    • Speed1
      ()
      DataArray
      {dims=[time: 2], unit=Hz, coords=[time]}
      Values:
      <scipp.DataArray>\n",
       "Dimensions: Sizes[time:2, ]\n",
       "Coordinates:\n",
       "  time                    datetime64             [ns]  (time)  [2011-08-12T15:50:16.800000000, 2011-08-12T17:22:05.000000000]\n",
       "Data:\n",
       "                            float64             [Hz]  (time)  [60, 60]\n",
       "\n",
       "
    • Speed2
      ()
      DataArray
      {dims=[time: 2], unit=Hz, coords=[time]}
      Values:
      <scipp.DataArray>\n",
       "Dimensions: Sizes[time:2, ]\n",
       "Coordinates:\n",
       "  time                    datetime64             [ns]  (time)  [2011-08-12T15:50:16.800000000, 2011-08-12T17:22:05.000000000]\n",
       "Data:\n",
       "                            float64             [Hz]  (time)  [60, 60]\n",
       "\n",
       "
    • Speed3
      ()
      DataArray
      {dims=[time: 2], unit=Hz, coords=[time]}
      Values:
      <scipp.DataArray>\n",
       "Dimensions: Sizes[time:2, ]\n",
       "Coordinates:\n",
       "  time                    datetime64             [ns]  (time)  [2011-08-12T15:50:16.800000000, 2011-08-12T17:22:05.000000000]\n",
       "Data:\n",
       "                            float64             [Hz]  (time)  [60, 60]\n",
       "\n",
       "
    • SpeedRequest1
      ()
      DataArray
      {dims=[time: 2], unit=Hz, coords=[time]}
      Values:
      <scipp.DataArray>\n",
       "Dimensions: Sizes[time:2, ]\n",
       "Coordinates:\n",
       "  time                    datetime64             [ns]  (time)  [2011-08-12T15:50:16.784000000, 2011-08-12T17:22:05.000000000]\n",
       "Data:\n",
       "                            float64             [Hz]  (time)  [60, 60]\n",
       "\n",
       "
    • SpeedRequest2
      ()
      DataArray
      {dims=[time: 2], unit=Hz, coords=[time]}
      Values:
      <scipp.DataArray>\n",
       "Dimensions: Sizes[time:2, ]\n",
       "Coordinates:\n",
       "  time                    datetime64             [ns]  (time)  [2011-08-12T15:50:16.784000000, 2011-08-12T17:22:05.000000000]\n",
       "Data:\n",
       "                            float64             [Hz]  (time)  [60, 60]\n",
       "\n",
       "
    • SpeedRequest3
      ()
      DataArray
      {dims=[time: 2], unit=Hz, coords=[time]}
      Values:
      <scipp.DataArray>\n",
       "Dimensions: Sizes[time:2, ]\n",
       "Coordinates:\n",
       "  time                    datetime64             [ns]  (time)  [2011-08-12T15:50:16.784000000, 2011-08-12T17:22:05.000000000]\n",
       "Data:\n",
       "                            float64             [Hz]  (time)  [60, 60]\n",
       "\n",
       "
    • TolRequest
      ()
      DataArray
      {dims=[time: 2], coords=[time]}
      Values:
      <scipp.DataArray>\n",
       "Dimensions: Sizes[time:2, ]\n",
       "Coordinates:\n",
       "  time                    datetime64             [ns]  (time)  [2011-08-12T15:50:16.800000000, 2011-08-12T17:22:05.000000000]\n",
       "Data:\n",
       "                            float64  [dimensionless]  (time)  [20, 20]\n",
       "\n",
       "
    • commErrs
      ()
      DataArray
      {dims=[time: 2], coords=[time]}
      Values:
      <scipp.DataArray>\n",
       "Dimensions: Sizes[time:2, ]\n",
       "Coordinates:\n",
       "  time                    datetime64             [ns]  (time)  [2011-08-12T15:50:16.784000000, 2011-08-12T17:22:05.000000000]\n",
       "Data:\n",
       "                            float64  [dimensionless]  (time)  [0, 0]\n",
       "\n",
       "
    • currentsample
      ()
      DataArray
      {dims=[time: 2], coords=[time]}
      Values:
      <scipp.DataArray>\n",
       "Dimensions: Sizes[time:2, ]\n",
       "Coordinates:\n",
       "  time                    datetime64             [ns]  (time)  [2011-08-12T15:50:16.784000000, 2011-08-12T17:22:05.000000000]\n",
       "Data:\n",
       "                            float64  [dimensionless]  (time)  [4, 4]\n",
       "\n",
       "
    • duration
      ()
      float64
      s
      5508.0
      Values:
      array(5508.)
    • end_time
      ()
      string
      2011-08-12T17:22:05
      Values:
      '2011-08-12T17:22:05'
    • experiment_identifier
      ()
      string
      IPTS-2767
      Values:
      'IPTS-2767'
    • fernsstatus
      ()
      DataArray
      {dims=[time: 2], coords=[time]}
      Values:
      <scipp.DataArray>\n",
       "Dimensions: Sizes[time:2, ]\n",
       "Coordinates:\n",
       "  time                    datetime64             [ns]  (time)  [2011-08-12T15:50:16.784000000, 2011-08-12T17:22:05.000000000]\n",
       "Data:\n",
       "                            float64  [dimensionless]  (time)  [3, 3]\n",
       "\n",
       "
    • file_notes
      ()
      string
      NONE
      Values:
      'NONE'
    • frequency
      ()
      DataArray
      {dims=[time: 330473], unit=Hz, coords=[time]}
      Values:
      <scipp.DataArray>\n",
       "Dimensions: Sizes[time:330473, ]\n",
       "Coordinates:\n",
       "  time                    datetime64             [ns]  (time)  [2011-08-12T15:50:17.000000000, 2011-08-12T15:50:17.016659999, ..., 2011-08-12T17:22:05.085449218, 2011-08-12T17:22:05.102050781]\n",
       "Data:\n",
       "                            float64             [Hz]  (time)  [0, 60.024, ..., 60.0024, 59.9988]\n",
       "\n",
       "
    • gd_prtn_chrg
      ()
      float64
      µAh
      1171.953902925
      Values:
      array(1171.95390292)
    • guide
      ()
      DataArray
      {dims=[time: 2], unit=mm, coords=[time]}
      Values:
      <scipp.DataArray>\n",
       "Dimensions: Sizes[time:2, ]\n",
       "Coordinates:\n",
       "  time                    datetime64             [ns]  (time)  [2011-08-12T15:50:16.784000000, 2011-08-12T17:22:05.000000000]\n",
       "Data:\n",
       "                            float64             [mm]  (time)  [-55.463, -55.463]\n",
       "\n",
       "
    • incident_beam
      ()
      vector3
      m
      [ 0. 0. 60.]
      Values:
      array([ 0.,  0., 60.])
    • instrument_name
      ()
      string
      POWGEN
      Values:
      'POWGEN'
    • monitor1
      ()
      DataArray
      {dims=[tof: 200001], unit=counts, coords=[source_position, position, tof]}
      Values:
      <scipp.DataArray>\n",
       "Dimensions: Sizes[tof:200001, ]\n",
       "Coordinates:\n",
       "  position                  vector3              [m]  ()  [(0, 0, -1)]\n",
       "  source_position           vector3              [m]  ()  [(0, 0, -60)]\n",
       "  tof                       float64            [µs]  (tof [bin-edge])  [0, 1, ..., 200000, 200001]\n",
       "Data:\n",
       "                            float64         [counts]  (tof)  [25, 10, ..., 0, 0]  [25, 10, ..., 0, 0]\n",
       "Attributes:\n",
       "  instrument_name            string  [dimensionless]  ()  ["POWGEN"]\n",
       "  sample_position           vector3              [m]  ()  [(0, 0, 0)]\n",
       "  spectrum                    int32  [dimensionless]  ()  [1]\n",
       "\n",
       "
    • position
      (spectrum)
      vector3
      m
      [ 1.17451004 -1.01106149 -2.03796699], [ 1.18147634 -0.95946649 -2.05334117], ..., [1.81428985 0.09565841 3.84338287], [1.81375055 0.1499371 3.84269584]
      Values:
      array([[ 1.17451004, -1.01106149, -2.03796699],\n",
       "       [ 1.18147634, -0.95946649, -2.05334117],\n",
       "       [ 1.18844265, -0.90787149, -2.06871534],\n",
       "       ...,\n",
       "       [ 1.81482915,  0.04137972,  3.8440699 ],\n",
       "       [ 1.81428985,  0.09565841,  3.84338287],\n",
       "       [ 1.81375055,  0.1499371 ,  3.84269584]])
    • proton_charge
      ()
      DataArray
      {dims=[time: 330473], unit=pC, coords=[time]}
      Values:
      <scipp.DataArray>\n",
       "Dimensions: Sizes[time:330473, ]\n",
       "Coordinates:\n",
       "  time                    datetime64             [ns]  (time)  [2011-08-12T15:50:17.000000000, 2011-08-12T15:50:17.016659999, ..., 2011-08-12T17:22:05.085449218, 2011-08-12T17:22:05.102050781]\n",
       "Data:\n",
       "                            float64             [pC]  (time)  [1.48434e+07, 1.48434e+07, ..., 1.48733e+07, 1.48391e+07]\n",
       "\n",
       "
    • run_number
      ()
      string
      4844
      Values:
      '4844'
    • run_start
      ()
      string
      2011-08-12T15:50:17
      Values:
      '2011-08-12T15:50:17'
    • run_title
      ()
      string
      diamond cw0.533 4.22e12 60Hz [10x30]
      Values:
      'diamond cw0.533 4.22e12 60Hz [10x30]'
    • running
      ()
      DataArray
      {dims=[time: 1], coords=[time]}
      Values:
      <scipp.DataArray>\n",
       "Dimensions: Sizes[time:1, ]\n",
       "Coordinates:\n",
       "  time                    datetime64             [ns]  (time)  [1990-01-01T00:00:00.000000000]\n",
       "Data:\n",
       "                               bool  [dimensionless]  (time)  [True]\n",
       "\n",
       "
    • s1b
      ()
      DataArray
      {dims=[time: 2], unit=mm, coords=[time]}
      Values:
      <scipp.DataArray>\n",
       "Dimensions: Sizes[time:2, ]\n",
       "Coordinates:\n",
       "  time                    datetime64             [ns]  (time)  [2011-08-12T15:50:16.784000000, 2011-08-12T17:22:05.000000000]\n",
       "Data:\n",
       "                            float64             [mm]  (time)  [20, 20]\n",
       "\n",
       "
    • s1l
      ()
      DataArray
      {dims=[time: 2], unit=mm, coords=[time]}
      Values:
      <scipp.DataArray>\n",
       "Dimensions: Sizes[time:2, ]\n",
       "Coordinates:\n",
       "  time                    datetime64             [ns]  (time)  [2011-08-12T15:50:16.800000000, 2011-08-12T17:22:05.000000000]\n",
       "Data:\n",
       "                            float64             [mm]  (time)  [5, 5]\n",
       "\n",
       "
    • s1r
      ()
      DataArray
      {dims=[time: 2], unit=mm, coords=[time]}
      Values:
      <scipp.DataArray>\n",
       "Dimensions: Sizes[time:2, ]\n",
       "Coordinates:\n",
       "  time                    datetime64             [ns]  (time)  [2011-08-12T15:50:16.784000000, 2011-08-12T17:22:05.000000000]\n",
       "Data:\n",
       "                            float64             [mm]  (time)  [-5, -5]\n",
       "\n",
       "
    • s1t
      ()
      DataArray
      {dims=[time: 2], unit=mm, coords=[time]}
      Values:
      <scipp.DataArray>\n",
       "Dimensions: Sizes[time:2, ]\n",
       "Coordinates:\n",
       "  time                    datetime64             [ns]  (time)  [2011-08-12T15:50:16.784000000, 2011-08-12T17:22:05.000000000]\n",
       "Data:\n",
       "                            float64             [mm]  (time)  [-10, -10]\n",
       "\n",
       "
    • sample
      ()
      PyObject
      <mantid.api._api.Sample object at 0x7f8b7b451eb0>
      Values:
      <mantid.api._api.Sample object at 0x7f8b7b451eb0>
    • sample_position
      ()
      vector3
      m
      [0. 0. 0.]
      Values:
      array([0., 0., 0.])
    • samplerequest
      ()
      DataArray
      {dims=[time: 2], coords=[time]}
      Values:
      <scipp.DataArray>\n",
       "Dimensions: Sizes[time:2, ]\n",
       "Coordinates:\n",
       "  time                    datetime64             [ns]  (time)  [2011-08-12T15:50:16.784000000, 2011-08-12T17:22:05.000000000]\n",
       "Data:\n",
       "                            float64  [dimensionless]  (time)  [4, 4]\n",
       "\n",
       "
    • scattered_beam
      (spectrum)
      vector3
      m
      [ 1.17451004 -1.01106149 -2.03796699], [ 1.18147634 -0.95946649 -2.05334117], ..., [1.81428985 0.09565841 3.84338287], [1.81375055 0.1499371 3.84269584]
      Values:
      array([[ 1.17451004, -1.01106149, -2.03796699],\n",
       "       [ 1.18147634, -0.95946649, -2.05334117],\n",
       "       [ 1.18844265, -0.90787149, -2.06871534],\n",
       "       ...,\n",
       "       [ 1.81482915,  0.04137972,  3.8440699 ],\n",
       "       [ 1.81428985,  0.09565841,  3.84338287],\n",
       "       [ 1.81375055,  0.1499371 ,  3.84269584]])
    • source_position
      ()
      vector3
      m
      [ 0. 0. -60.]
      Values:
      array([  0.,   0., -60.])
    • start_time
      ()
      string
      2011-08-12T15:50:17
      Values:
      '2011-08-12T15:50:17'
    • tof
      (wavelength [bin-edge])
      float64
      µs
      19.0, 1.669e+04
      Values:
      array([   19.        , 16694.30078125])
    • vGuide
      ()
      DataArray
      {dims=[time: 2], coords=[time]}
      Values:
      <scipp.DataArray>\n",
       "Dimensions: Sizes[time:2, ]\n",
       "Coordinates:\n",
       "  time                    datetime64             [ns]  (time)  [2011-08-12T15:50:16.800000000, 2011-08-12T17:22:05.000000000]\n",
       "Data:\n",
       "                            float64  [dimensionless]  (time)  [2, 2]\n",
       "\n",
       "
    • veto_pulse_time
      ()
      DataArray
      {dims=[time: 1], coords=[time]}
      Values:
      <scipp.DataArray>\n",
       "Dimensions: Sizes[time:1, ]\n",
       "Coordinates:\n",
       "  time                    datetime64             [ns]  (time)  [2011-08-12T15:50:17.000000000]\n",
       "Data:\n",
       "                            float64  [dimensionless]  (time)  [0]\n",
       "\n",
       "
" ], "text/plain": [ "\n", "Dimensions: Sizes[spectrum:24794, wavelength:1, ]\n", "Coordinates:\n", " spectrum int32 [dimensionless] (spectrum) [1, 2, ..., 24793, 24794]\n", " two_theta float64 [rad] (spectrum) [2.49144, 2.50373, ..., 0.441579, 0.442313]\n", " wavelength float64 [Å] (spectrum, wavelength [bin-edge]) [0.00120148, 1.05567, ..., 0.00116984, 1.02788]\n", "Data:\n", " DataArrayView [dimensionless] (spectrum, wavelength) [\n", "Dimensions: Sizes[event:0, ]\n", "Coordinates:\n", " wavelength float64 [Å] (event) []\n", "Data:\n", " float32 [counts] (event) [] []\n", "Attributes:\n", " tof float64 [µs] (event) []\n", "\n", ", \n", "Dimensions: Sizes[event:0, ]\n", "Coordinates:\n", " wavelength float64 [Å] (event) []\n", "Data:\n", " float32 [counts] (event) [] []\n", "Attributes:\n", " tof float64 [µs] (event) []\n", "\n", ", ..., \n", "Dimensions: Sizes[event:0, ]\n", "Coordinates:\n", " wavelength float64 [Å] (event) []\n", "Data:\n", " float32 [counts] (event) [] []\n", "Attributes:\n", " tof float64 [µs] (event) []\n", "\n", ", \n", "Dimensions: Sizes[event:0, ]\n", "Coordinates:\n", " wavelength float64 [Å] (event) []\n", "Data:\n", " float32 [counts] (event) [] []\n", "Attributes:\n", " tof float64 [µs] (event) []\n", "\n", "]\n", "Attributes:\n", " ChopperStatus1 DataArray [dimensionless] () [\n", "Dimensions: Sizes[time:2, ]\n", "Coordinates:\n", " time datetime64 [ns] (time) [2011-08-12T15:50:16.800000000, 2011-08-12T17:22:05.000000000]\n", "Data:\n", " float64 [dimensionless] (time) [4, 4]\n", "\n", "]\n", " ChopperStatus2 DataArray [dimensionless] () [\n", "Dimensions: Sizes[time:2, ]\n", "Coordinates:\n", " time datetime64 [ns] (time) [2011-08-12T15:50:16.816000000, 2011-08-12T17:22:05.000000000]\n", "Data:\n", " float64 [dimensionless] (time) [4, 4]\n", "\n", "]\n", " ChopperStatus3 DataArray [dimensionless] () [\n", "Dimensions: Sizes[time:2, ]\n", "Coordinates:\n", " time datetime64 [ns] (time) [2011-08-12T15:50:16.816000000, 2011-08-12T17:22:05.000000000]\n", "Data:\n", " float64 [dimensionless] (time) [4, 4]\n", "\n", "]\n", " CurrentSP DataArray [dimensionless] () [\n", "Dimensions: Sizes[time:2, ]\n", "Coordinates:\n", " time datetime64 [ns] (time) [2011-08-12T15:50:16.800000000, 2011-08-12T17:22:05.000000000]\n", "Data:\n", " float64 [dimensionless] (time) [300, 300]\n", "\n", "]\n", " EnergyRequest DataArray [dimensionless] () [\n", "Dimensions: Sizes[time:2, ]\n", "Coordinates:\n", " time datetime64 [ns] (time) [2011-08-12T15:50:16.784000000, 2011-08-12T17:22:05.000000000]\n", "Data:\n", " float64 [meV] (time) [287.955, 287.955]\n", "\n", "]\n", " Filename string [dimensionless] () [\"/home/scientist/.cache/scippneutron/2/PG3_4844_event.nxs\"]\n", " L1 float64 [m] () [60]\n", " L2 float64 [m] (spectrum) [2.56028, 2.55591, ..., 4.25116, 4.25188]\n", " LKSRampRate DataArray [dimensionless] () [\n", "Dimensions: Sizes[time:2, ]\n", "Coordinates:\n", " time datetime64 [ns] (time) [2011-08-12T15:50:16.800000000, 2011-08-12T17:22:05.000000000]\n", "Data:\n", " float64 [dimensionless] (time) [0, 0]\n", "\n", "]\n", " LambdaRequest DataArray [dimensionless] () [\n", "Dimensions: Sizes[time:2, ]\n", "Coordinates:\n", " time datetime64 [ns] (time) [2011-08-12T15:50:16.784000000, 2011-08-12T17:22:05.000000000]\n", "Data:\n", " float64 [Å] (time) [0.533, 0.533]\n", "\n", "]\n", " Ltotal float64 [m] (spectrum) [62.5603, 62.5559, ..., 64.2512, 64.2519]\n", " Phase1 DataArray [dimensionless] () [\n", "Dimensions: Sizes[time:1794, ]\n", "Coordinates:\n", " time datetime64 [ns] (time) [2011-08-12T15:50:16.800000000, 2011-08-12T15:50:18.267000007, ..., 2011-08-12T17:22:02.999218750, 2011-08-12T17:22:05.000000000]\n", "Data:\n", " float64 [µs] (time) [8166.72, 8165.16, ..., 8163.85, 8163.85]\n", "\n", "]\n", " Phase2 DataArray [dimensionless] () [\n", "Dimensions: Sizes[time:1793, ]\n", "Coordinates:\n", " time datetime64 [ns] (time) [2011-08-12T15:50:16.800000000, 2011-08-12T15:50:18.267000007, ..., 2011-08-12T17:22:02.999218750, 2011-08-12T17:22:05.000000000]\n", "Data:\n", " float64 [µs] (time) [8335.63, 8334.09, ..., 8332.86, 8332.86]\n", "\n", "]\n", " Phase3 DataArray [dimensionless] () [\n", "Dimensions: Sizes[time:1777, ]\n", "Coordinates:\n", " time datetime64 [ns] (time) [2011-08-12T15:50:16.800000000, 2011-08-12T15:50:19.688999938, ..., 2011-08-12T17:22:04.421093750, 2011-08-12T17:22:05.000000000]\n", "Data:\n", " float64 [µs] (time) [14001.5, 14001.4, ..., 14002, 14002]\n", "\n", "]\n", " PhaseRequest1 DataArray [dimensionless] () [\n", "Dimensions: Sizes[time:2, ]\n", "Coordinates:\n", " time datetime64 [ns] (time) [2011-08-12T15:50:16.800000000, 2011-08-12T17:22:05.000000000]\n", "Data:\n", " float64 [µs] (time) [8164.07, 8164.07]\n", "\n", "]\n", " PhaseRequest2 DataArray [dimensionless] () [\n", "Dimensions: Sizes[time:2, ]\n", "Coordinates:\n", " time datetime64 [ns] (time) [2011-08-12T15:50:16.800000000, 2011-08-12T17:22:05.000000000]\n", "Data:\n", " float64 [µs] (time) [8332.89, 8332.89]\n", "\n", "]\n", " PhaseRequest3 DataArray [dimensionless] () [\n", "Dimensions: Sizes[time:2, ]\n", "Coordinates:\n", " time datetime64 [ns] (time) [2011-08-12T15:50:16.800000000, 2011-08-12T17:22:05.000000000]\n", "Data:\n", " float64 [µs] (time) [14001.8, 14001.8]\n", "\n", "]\n", " S1HCenter DataArray [dimensionless] () [\n", "Dimensions: Sizes[time:2, ]\n", "Coordinates:\n", " time datetime64 [ns] (time) [2011-08-12T15:50:16.800000000, 2011-08-12T17:22:05.000000000]\n", "Data:\n", " float64 [mm] (time) [0, 0]\n", "\n", "]\n", " S1HCenterOffset DataArray [dimensionless] () [\n", "Dimensions: Sizes[time:2, ]\n", "Coordinates:\n", " time datetime64 [ns] (time) [2011-08-12T15:50:16.800000000, 2011-08-12T17:22:05.000000000]\n", "Data:\n", " float64 [dimensionless] (time) [0, 0]\n", "\n", "]\n", " S1HWidth DataArray [dimensionless] () [\n", "Dimensions: Sizes[time:2, ]\n", "Coordinates:\n", " time datetime64 [ns] (time) [2011-08-12T15:50:16.800000000, 2011-08-12T17:22:05.000000000]\n", "Data:\n", " float64 [mm] (time) [10, 10]\n", "\n", "]\n", " S1VCenter DataArray [dimensionless] () [\n", "Dimensions: Sizes[time:2, ]\n", "Coordinates:\n", " time datetime64 [ns] (time) [2011-08-12T15:50:16.800000000, 2011-08-12T17:22:05.000000000]\n", "Data:\n", " float64 [mm] (time) [5, 5]\n", "\n", "]\n", " S1VCenterOffset DataArray [dimensionless] () [\n", "Dimensions: Sizes[time:2, ]\n", "Coordinates:\n", " time datetime64 [ns] (time) [2011-08-12T15:50:16.800000000, 2011-08-12T17:22:05.000000000]\n", "Data:\n", " float64 [dimensionless] (time) [0, 0]\n", "\n", "]\n", " S1VHeight DataArray [dimensionless] () [\n", "Dimensions: Sizes[time:2, ]\n", "Coordinates:\n", " time datetime64 [ns] (time) [2011-08-12T15:50:16.816000000, 2011-08-12T17:22:05.000000000]\n", "Data:\n", " float64 [mm] (time) [30, 30]\n", "\n", "]\n", " SampleTemp DataArray [dimensionless] () [\n", "Dimensions: Sizes[time:467, ]\n", "Coordinates:\n", " time datetime64 [ns] (time) [2011-08-12T15:50:16.784000000, 2011-08-12T15:50:18.001000007, ..., 2011-08-12T17:21:58.623843750, 2011-08-12T17:22:05.000000000]\n", "Data:\n", " float64 [dimensionless] (time) [299.352, 299.446, ..., 300, 300]\n", "\n", "]\n", " Speed1 DataArray [dimensionless] () [\n", "Dimensions: Sizes[time:2, ]\n", "Coordinates:\n", " time datetime64 [ns] (time) [2011-08-12T15:50:16.800000000, 2011-08-12T17:22:05.000000000]\n", "Data:\n", " float64 [Hz] (time) [60, 60]\n", "\n", "]\n", " Speed2 DataArray [dimensionless] () [\n", "Dimensions: Sizes[time:2, ]\n", "Coordinates:\n", " time datetime64 [ns] (time) [2011-08-12T15:50:16.800000000, 2011-08-12T17:22:05.000000000]\n", "Data:\n", " float64 [Hz] (time) [60, 60]\n", "\n", "]\n", " Speed3 DataArray [dimensionless] () [\n", "Dimensions: Sizes[time:2, ]\n", "Coordinates:\n", " time datetime64 [ns] (time) [2011-08-12T15:50:16.800000000, 2011-08-12T17:22:05.000000000]\n", "Data:\n", " float64 [Hz] (time) [60, 60]\n", "\n", "]\n", " SpeedRequest1 DataArray [dimensionless] () [\n", "Dimensions: Sizes[time:2, ]\n", "Coordinates:\n", " time datetime64 [ns] (time) [2011-08-12T15:50:16.784000000, 2011-08-12T17:22:05.000000000]\n", "Data:\n", " float64 [Hz] (time) [60, 60]\n", "\n", "]\n", " SpeedRequest2 DataArray [dimensionless] () [\n", "Dimensions: Sizes[time:2, ]\n", "Coordinates:\n", " time datetime64 [ns] (time) [2011-08-12T15:50:16.784000000, 2011-08-12T17:22:05.000000000]\n", "Data:\n", " float64 [Hz] (time) [60, 60]\n", "\n", "]\n", " SpeedRequest3 DataArray [dimensionless] () [\n", "Dimensions: Sizes[time:2, ]\n", "Coordinates:\n", " time datetime64 [ns] (time) [2011-08-12T15:50:16.784000000, 2011-08-12T17:22:05.000000000]\n", "Data:\n", " float64 [Hz] (time) [60, 60]\n", "\n", "]\n", " TolRequest DataArray [dimensionless] () [\n", "Dimensions: Sizes[time:2, ]\n", "Coordinates:\n", " time datetime64 [ns] (time) [2011-08-12T15:50:16.800000000, 2011-08-12T17:22:05.000000000]\n", "Data:\n", " float64 [dimensionless] (time) [20, 20]\n", "\n", "]\n", " commErrs DataArray [dimensionless] () [\n", "Dimensions: Sizes[time:2, ]\n", "Coordinates:\n", " time datetime64 [ns] (time) [2011-08-12T15:50:16.784000000, 2011-08-12T17:22:05.000000000]\n", "Data:\n", " float64 [dimensionless] (time) [0, 0]\n", "\n", "]\n", " currentsample DataArray [dimensionless] () [\n", "Dimensions: Sizes[time:2, ]\n", "Coordinates:\n", " time datetime64 [ns] (time) [2011-08-12T15:50:16.784000000, 2011-08-12T17:22:05.000000000]\n", "Data:\n", " float64 [dimensionless] (time) [4, 4]\n", "\n", "]\n", " duration float64 [s] () [5508]\n", " end_time string [dimensionless] () [\"2011-08-12T17:22:05\"]\n", " experiment_identifier string [dimensionless] () [\"IPTS-2767\"]\n", " fernsstatus DataArray [dimensionless] () [\n", "Dimensions: Sizes[time:2, ]\n", "Coordinates:\n", " time datetime64 [ns] (time) [2011-08-12T15:50:16.784000000, 2011-08-12T17:22:05.000000000]\n", "Data:\n", " float64 [dimensionless] (time) [3, 3]\n", "\n", "]\n", " file_notes string [dimensionless] () [\"NONE\"]\n", " frequency DataArray [dimensionless] () [\n", "Dimensions: Sizes[time:330473, ]\n", "Coordinates:\n", " time datetime64 [ns] (time) [2011-08-12T15:50:17.000000000, 2011-08-12T15:50:17.016659999, ..., 2011-08-12T17:22:05.085449218, 2011-08-12T17:22:05.102050781]\n", "Data:\n", " float64 [Hz] (time) [0, 60.024, ..., 60.0024, 59.9988]\n", "\n", "]\n", " gd_prtn_chrg float64 [µAh] () [1171.95]\n", " guide DataArray [dimensionless] () [\n", "Dimensions: Sizes[time:2, ]\n", "Coordinates:\n", " time datetime64 [ns] (time) [2011-08-12T15:50:16.784000000, 2011-08-12T17:22:05.000000000]\n", "Data:\n", " float64 [mm] (time) [-55.463, -55.463]\n", "\n", "]\n", " incident_beam vector3 [m] () [(0, 0, 60)]\n", " instrument_name string [dimensionless] () [\"POWGEN\"]\n", " monitor1 DataArray [dimensionless] () [\n", "Dimensions: Sizes[tof:200001, ]\n", "Coordinates:\n", " position vector3 [m] () [(0, 0, -1)]\n", " source_position vector3 [m] () [(0, 0, -60)]\n", " tof float64 [µs] (tof [bin-edge]) [0, 1, ..., 200000, 200001]\n", "Data:\n", " float64 [counts] (tof) [25, 10, ..., 0, 0] [25, 10, ..., 0, 0]\n", "Attributes:\n", " instrument_name string [dimensionless] () [\"POWGEN\"]\n", " sample_position vector3 [m] () [(0, 0, 0)]\n", " spectrum int32 [dimensionless] () [1]\n", "\n", "]\n", " position vector3 [m] (spectrum) [(1.17451, -1.01106, -2.03797), (1.18148, -0.959466, -2.05334), ..., (1.81429, 0.0956584, 3.84338), (1.81375, 0.149937, 3.8427)]\n", " proton_charge DataArray [dimensionless] () [\n", "Dimensions: Sizes[time:330473, ]\n", "Coordinates:\n", " time datetime64 [ns] (time) [2011-08-12T15:50:17.000000000, 2011-08-12T15:50:17.016659999, ..., 2011-08-12T17:22:05.085449218, 2011-08-12T17:22:05.102050781]\n", "Data:\n", " float64 [pC] (time) [1.48434e+07, 1.48434e+07, ..., 1.48733e+07, 1.48391e+07]\n", "\n", "]\n", " run_number string [dimensionless] () [\"4844\"]\n", " run_start string [dimensionless] () [\"2011-08-12T15:50:17\"]\n", " run_title string [dimensionless] () [\"diamond cw0.533 4.22e12 60Hz [10x30]\"]\n", " running DataArray [dimensionless] () [\n", "Dimensions: Sizes[time:1, ]\n", "Coordinates:\n", " time datetime64 [ns] (time) [1990-01-01T00:00:00.000000000]\n", "Data:\n", " bool [dimensionless] (time) [True]\n", "\n", "]\n", " s1b DataArray [dimensionless] () [\n", "Dimensions: Sizes[time:2, ]\n", "Coordinates:\n", " time datetime64 [ns] (time) [2011-08-12T15:50:16.784000000, 2011-08-12T17:22:05.000000000]\n", "Data:\n", " float64 [mm] (time) [20, 20]\n", "\n", "]\n", " s1l DataArray [dimensionless] () [\n", "Dimensions: Sizes[time:2, ]\n", "Coordinates:\n", " time datetime64 [ns] (time) [2011-08-12T15:50:16.800000000, 2011-08-12T17:22:05.000000000]\n", "Data:\n", " float64 [mm] (time) [5, 5]\n", "\n", "]\n", " s1r DataArray [dimensionless] () [\n", "Dimensions: Sizes[time:2, ]\n", "Coordinates:\n", " time datetime64 [ns] (time) [2011-08-12T15:50:16.784000000, 2011-08-12T17:22:05.000000000]\n", "Data:\n", " float64 [mm] (time) [-5, -5]\n", "\n", "]\n", " s1t DataArray [dimensionless] () [\n", "Dimensions: Sizes[time:2, ]\n", "Coordinates:\n", " time datetime64 [ns] (time) [2011-08-12T15:50:16.784000000, 2011-08-12T17:22:05.000000000]\n", "Data:\n", " float64 [mm] (time) [-10, -10]\n", "\n", "]\n", " sample PyObject [dimensionless] () []\n", " sample_position vector3 [m] () [(0, 0, 0)]\n", " samplerequest DataArray [dimensionless] () [\n", "Dimensions: Sizes[time:2, ]\n", "Coordinates:\n", " time datetime64 [ns] (time) [2011-08-12T15:50:16.784000000, 2011-08-12T17:22:05.000000000]\n", "Data:\n", " float64 [dimensionless] (time) [4, 4]\n", "\n", "]\n", " scattered_beam vector3 [m] (spectrum) [(1.17451, -1.01106, -2.03797), (1.18148, -0.959466, -2.05334), ..., (1.81429, 0.0956584, 3.84338), (1.81375, 0.149937, 3.8427)]\n", " source_position vector3 [m] () [(0, 0, -60)]\n", " start_time string [dimensionless] () [\"2011-08-12T15:50:17\"]\n", " tof float64 [µs] (wavelength [bin-edge]) [19, 16694.3]\n", " vGuide DataArray [dimensionless] () [\n", "Dimensions: Sizes[time:2, ]\n", "Coordinates:\n", " time datetime64 [ns] (time) [2011-08-12T15:50:16.800000000, 2011-08-12T17:22:05.000000000]\n", "Data:\n", " float64 [dimensionless] (time) [2, 2]\n", "\n", "]\n", " veto_pulse_time DataArray [dimensionless] () [\n", "Dimensions: Sizes[time:1, ]\n", "Coordinates:\n", " time datetime64 [ns] (time) [2011-08-12T15:50:17.000000000]\n", "Data:\n", " float64 [dimensionless] (time) [0]\n", "\n", "]\n" ] }, "execution_count": 15, "metadata": {}, "output_type": "execute_result" } ], "source": [ "sample_over_mon = sample_lambda.bins / sc.lookup(func=mon, dim='wavelength')\n", "sample_over_mon" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "Finally, we can plot the event data with on-the-fly binning.\n", "By default, the `plot` function uses the coordinates of the binned data to define histogram edges, which, in this case, would give a single bin along the `'tof'` dimension.\n", "For a better representation of the data, we supply the bin edges to obtain a finer binning, yielding a more meaningful figure:" ] }, { "cell_type": "code", "execution_count": 16, "metadata": {}, "outputs": [ { "data": { "application/vnd.jupyter.widget-view+json": { "model_id": "963161521cf94b919001fb5a2b6210dc", "version_major": 2, "version_minor": 0 }, "text/plain": [ "VBox(children=(HBox(children=(VBox(children=(Button(icon='home', layout=Layout(padding='0px 0px 0px 0px', widt…" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "sc.plot(sample_over_mon)" ] }, { "cell_type": "code", "execution_count": 17, "metadata": {}, "outputs": [], "source": [ "del sample_lambda\n", "del sample_over_mon\n", "del sample\n", "del vanadium" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "### From events to histogram\n", "\n", "*Note: We are continuing here with data that has not been normalized to the monitors.*\n", "\n", "We histogram the event data:" ] }, { "cell_type": "code", "execution_count": 18, "metadata": {}, "outputs": [ { "data": { "text/html": [ " vanadiumvanadium(dims=['spectrum', 'dspacing'], shape=[3, 7], unit=counts, variances=True)variances spectrumdspacingvalues spectrumdspacing samplesample(dims=['spectrum', 'dspacing'], shape=[3, 7], unit=counts, variances=True)variances spectrumdspacingvalues spectrumdspacing spect..spectrum(dims=['spectrum'], shape=[3], unit=dimensionless, variances=False)values spectrum dspacingdspacing(dims=['dspacing'], shape=[8], unit=Å, variances=False)values dspacing" ], "text/plain": [ "" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "dspacing_bins = sc.arange(dim='dspacing', start=0.3, stop=2.0, step=0.001, unit=sc.units.angstrom)\n", "hist = sc.Dataset(data={'sample':sc.histogram(dspacing_sample, bins=dspacing_bins),\n", " 'vanadium':sc.histogram(dspacing_vanadium, bins=dspacing_bins)})\n", "sc.show(hist['spectrum', 0:3]['dspacing', 0:7])" ] }, { "cell_type": "code", "execution_count": 19, "metadata": {}, "outputs": [ { "data": { "application/vnd.jupyter.widget-view+json": { "model_id": "635462d6375c4847a7d5b141fb228879", "version_major": 2, "version_minor": 0 }, "text/plain": [ "VBox(children=(VBox(children=(HBox(children=(VBox(children=(Button(icon='home', layout=Layout(padding='0px 0px…" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "sc.plot(hist)" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "### Summing (focussing) and normalizing\n", "\n", "After conversion to `'dspacing'`, generic `sum` and `/` operations can be used to \"focus\" and normalize the diffraction data to the vanadium run:" ] }, { "cell_type": "code", "execution_count": 20, "metadata": {}, "outputs": [ { "data": { "application/vnd.jupyter.widget-view+json": { "model_id": "eb8c0ea4786745329d47570cb1cf31b2", "version_major": 2, "version_minor": 0 }, "text/plain": [ "VBox(children=(HBox(children=(VBox(children=(Button(icon='home', layout=Layout(padding='0px 0px 0px 0px', widt…" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "summed = sc.sum(hist, 'spectrum')\n", "sc.plot(summed)" ] }, { "cell_type": "code", "execution_count": 21, "metadata": {}, "outputs": [ { "data": { "application/vnd.jupyter.widget-view+json": { "model_id": "b5c0d095000749d884992c23f2e3817b", "version_major": 2, "version_minor": 0 }, "text/plain": [ "VBox(children=(HBox(children=(VBox(children=(Button(icon='home', layout=Layout(padding='0px 0px 0px 0px', widt…" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "normalized = summed['sample'] / summed['vanadium']\n", "sc.plot(normalized)" ] }, { "cell_type": "markdown", "metadata": { "tags": [] }, "source": [ "### Focussing with $\\theta$ dependence in event-mode\n", "\n", "The approach used above combines reflections from all crystallographic planes and is therefore of limited use.\n", "We can use `groupby` to focus each of multiple groups of spectra into a distinct output spectrum.\n", "Here we define groups based on a range of scattering angles – a simple $\\theta$-dependent binning.\n", "This also demonstrates how we can postpone histogramming until after the focussing step." ] }, { "cell_type": "code", "execution_count": 22, "metadata": {}, "outputs": [ { "ename": "AttributeError", "evalue": "'GroupbyBins' object has no attribute 'concatenate'", "output_type": "error", "traceback": [ "\u001b[0;31m---------------------------------------------------------------------------\u001b[0m", "\u001b[0;31mAttributeError\u001b[0m Traceback (most recent call last)", "\u001b[0;32m/tmp/ipykernel_25962/1387100095.py\u001b[0m in \u001b[0;36m\u001b[0;34m\u001b[0m\n\u001b[1;32m 1\u001b[0m \u001b[0mtwo_theta\u001b[0m \u001b[0;34m=\u001b[0m \u001b[0msc\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mlinspace\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0mdim\u001b[0m\u001b[0;34m=\u001b[0m\u001b[0;34m'two_theta'\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0munit\u001b[0m\u001b[0;34m=\u001b[0m\u001b[0;34m'rad'\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0mstart\u001b[0m\u001b[0;34m=\u001b[0m\u001b[0;36m0.0\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0mstop\u001b[0m\u001b[0;34m=\u001b[0m\u001b[0mnp\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mpi\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0mnum\u001b[0m\u001b[0;34m=\u001b[0m\u001b[0;36m16\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m 2\u001b[0m \u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0;32m----> 3\u001b[0;31m \u001b[0mfocussed_sample\u001b[0m \u001b[0;34m=\u001b[0m \u001b[0msc\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mgroupby\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0mdspacing_sample\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0;34m'two_theta'\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0mbins\u001b[0m\u001b[0;34m=\u001b[0m\u001b[0mtwo_theta\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mbins\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mconcatenate\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0;34m'spectrum'\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0m\u001b[1;32m 4\u001b[0m \u001b[0mfocussed_vanadium\u001b[0m \u001b[0;34m=\u001b[0m \u001b[0msc\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mgroupby\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0mdspacing_vanadium\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0;34m'two_theta'\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0mbins\u001b[0m\u001b[0;34m=\u001b[0m\u001b[0mtwo_theta\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mbins\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mconcatenate\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0;34m'spectrum'\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m 5\u001b[0m \u001b[0mnorm\u001b[0m \u001b[0;34m=\u001b[0m \u001b[0msc\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mhistogram\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0mfocussed_vanadium\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0mbins\u001b[0m\u001b[0;34m=\u001b[0m\u001b[0mdspacing_bins\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n", "\u001b[0;31mAttributeError\u001b[0m: 'GroupbyBins' object has no attribute 'concatenate'" ] } ], "source": [ "two_theta = sc.linspace(dim='two_theta', unit='rad', start=0.0, stop=np.pi, num=16)\n", "\n", "focussed_sample = sc.groupby(dspacing_sample, 'two_theta', bins=two_theta).bins.concatenate('spectrum')\n", "focussed_vanadium = sc.groupby(dspacing_vanadium, 'two_theta', bins=two_theta).bins.concatenate('spectrum')\n", "norm = sc.histogram(focussed_vanadium, bins=dspacing_bins)\n", "focussed_sample.bins /= sc.lookup(func=norm, dim='dspacing')\n", "normalized = sc.histogram(focussed_sample, bins=dspacing_bins)" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "The normalized output looks as follows:" ] }, { "cell_type": "code", "execution_count": 23, "metadata": {}, "outputs": [ { "data": { "application/vnd.jupyter.widget-view+json": { "model_id": "c4082c9b845e4c43bd214bf8fa8d8f29", "version_major": 2, "version_minor": 0 }, "text/plain": [ "VBox(children=(HBox(children=(VBox(children=(Button(icon='home', layout=Layout(padding='0px 0px 0px 0px', widt…" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "sc.plot(normalized, vmin=sc.scalar(0), vmax=sc.scalar(0.4))" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "As a bonus, we can use slicing and a dict-comprehension to quickly create of plot comparing the spectra for different scattering angle bins:" ] }, { "cell_type": "code", "execution_count": 24, "metadata": {}, "outputs": [ { "ename": "NotFoundError", "evalue": "Expected \n dspacing: (dspacing: 1700) float64 [Å] [0.3, 0.301, ..., 1.998, 1.999] to contain two_theta.", "output_type": "error", "traceback": [ "\u001b[0;31m---------------------------------------------------------------------------\u001b[0m", "\u001b[0;31mNotFoundError\u001b[0m Traceback (most recent call last)", "\u001b[0;32m/tmp/ipykernel_25962/3747691360.py\u001b[0m in \u001b[0;36m\u001b[0;34m\u001b[0m\n\u001b[1;32m 1\u001b[0m \u001b[0;31m# compute centers of theta bins\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0;32m----> 2\u001b[0;31m \u001b[0mangles\u001b[0m \u001b[0;34m=\u001b[0m \u001b[0mnormalized\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mcoords\u001b[0m\u001b[0;34m[\u001b[0m\u001b[0;34m'two_theta'\u001b[0m\u001b[0;34m]\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mvalues\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0m\u001b[1;32m 3\u001b[0m \u001b[0mangles\u001b[0m \u001b[0;34m=\u001b[0m \u001b[0;36m0.5\u001b[0m\u001b[0;34m*\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0mangles\u001b[0m\u001b[0;34m[\u001b[0m\u001b[0;36m1\u001b[0m\u001b[0;34m:\u001b[0m\u001b[0;34m]\u001b[0m \u001b[0;34m+\u001b[0m \u001b[0mangles\u001b[0m\u001b[0;34m[\u001b[0m\u001b[0;34m:\u001b[0m\u001b[0;34m-\u001b[0m\u001b[0;36m1\u001b[0m\u001b[0;34m]\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m 4\u001b[0m sc.plot({f'{round(angles[group], 3)} rad':\n\u001b[1;32m 5\u001b[0m \u001b[0mnormalized\u001b[0m\u001b[0;34m[\u001b[0m\u001b[0;34m'dspacing'\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0;36m300\u001b[0m\u001b[0;34m:\u001b[0m\u001b[0;36m500\u001b[0m\u001b[0;34m]\u001b[0m\u001b[0;34m[\u001b[0m\u001b[0;34m'two_theta'\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0mgroup\u001b[0m\u001b[0;34m]\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n", "\u001b[0;31mNotFoundError\u001b[0m: Expected \n dspacing: (dspacing: 1700) float64 [Å] [0.3, 0.301, ..., 1.998, 1.999] to contain two_theta." ] } ], "source": [ "# compute centers of theta bins\n", "angles = normalized.coords['two_theta'].values\n", "angles = 0.5*(angles[1:] + angles[:-1])\n", "sc.plot({f'{round(angles[group], 3)} rad':\n", " normalized['dspacing', 300:500]['two_theta', group]\n", " for group in range(2,6)})" ] }, { "cell_type": "code", "execution_count": null, "metadata": {}, "outputs": [], "source": [] } ], "metadata": { "kernelspec": { "display_name": "Python 3 (ipykernel)", "language": "python", "name": "python3" }, "language_info": { "codemirror_mode": { "name": "ipython", "version": 3 }, "file_extension": ".py", "mimetype": "text/x-python", "name": "python", "nbconvert_exporter": "python", "pygments_lexer": "ipython3", "version": "3.9.9" } }, "nbformat": 4, "nbformat_minor": 4 }