View job log here
scipp._scipp.core.BinEdgeError: The input does not have coordinates with bin-edges.
In provider ess.reduce.unwrap.to_wavelength.monitor_wavelength_data(
monitor_data: ess.reduce.nexus.types.RawMonitor[ess.imaging.types.SampleRun, ess.imaging.types.BeamMonitor1],
lookup: ess.reduce.unwrap.types.ErrorLimitedLookupTable[ess.imaging.types.SampleRun, ess.imaging.types.BeamMonitor1],
ltotal: ess.reduce.unwrap.types.MonitorLtotal[ess.imaging.types.SampleRun, ess.imaging.types.BeamMonitor1],
pulse_stride_offset: ess.reduce.unwrap.types.PulseStrideOffset,
keep_event_time_offset: ess.reduce.unwrap.types.KeepEventTimeOffset,
) -> ess.reduce.unwrap.types.WavelengthMonitor[ess.imaging.types.SampleRun, ess.imaging.types.BeamMonitor1]
.tox/nexusfiles-scipp-odin/lib/python3.12/site-packages/scipp/core/binning.py:992: BinEdgeError
Full output
workflow = <sciline.pipeline.Pipeline object at 0x7fa63a4f77a0>
coda_nexus_file_path = PosixPath('/ess/data/coda/999999/raw/coda_odin_999999_00019227.hdf')
monitor_index = 1
@pytest.mark.parametrize("monitor_index", [1, 2, 3])
def test_can_compute_wavelength__beam_monitor_(
workflow: sciline.Pipeline, coda_nexus_file_path: Path, monitor_index: int
) -> None:
workflow[Filename[SampleRun]] = coda_nexus_file_path
workflow[NeXusName[BeamMonitor1]] = f"beam_monitor_{monitor_index}"
> result = workflow.compute(WavelengthMonitor[SampleRun, BeamMonitor1])
^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
tests/nexusfiles-scipp/odin/odin_reduction_test.py:46:
_ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _
.tox/nexusfiles-scipp-odin/lib/python3.12/site-packages/sciline/pipeline.py:191: in compute
return self.get(tp, **kwargs).compute(reporter=reporter)
^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
.tox/nexusfiles-scipp-odin/lib/python3.12/site-packages/sciline/task_graph.py:122: in compute
return self._scheduler.get(self._graph, [targets], reporter=reporter)[0]
^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
.tox/nexusfiles-scipp-odin/lib/python3.12/site-packages/sciline/scheduler.py:140: in get
return self._dask_get(dsk, list(map(_to_dask_key, keys)))
^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
.tox/nexusfiles-scipp-odin/lib/python3.12/site-packages/dask/threaded.py:115: in get
results = get_async(
.tox/nexusfiles-scipp-odin/lib/python3.12/site-packages/dask/local.py:547: in get_async
raise_exception(exc, tb)
.tox/nexusfiles-scipp-odin/lib/python3.12/site-packages/dask/local.py:351: in reraise
raise exc
.tox/nexusfiles-scipp-odin/lib/python3.12/site-packages/dask/local.py:256: in execute_task
result = task(data)
^^^^^^^^^^
.tox/nexusfiles-scipp-odin/lib/python3.12/site-packages/dask/_task_spec.py:768: in __call__
return self.func(*new_argspec)
^^^^^^^^^^^^^^^^^^^^^^^
.tox/nexusfiles-scipp-odin/lib/python3.12/site-packages/dask/utils.py:81: in apply
return func(*args)
^^^^^^^^^^^
.tox/nexusfiles-scipp-odin/lib/python3.12/site-packages/sciline/_provider.py:145: in __call__
return self._func(*args, **kwargs)
^^^^^^^^^^^^^^^^^^^^^^^^^^^
.tox/nexusfiles-scipp-odin/lib/python3.12/site-packages/ess/reduce/unwrap/to_wavelength.py:565: in monitor_wavelength_data
_compute_wavelength_data(
.tox/nexusfiles-scipp-odin/lib/python3.12/site-packages/ess/reduce/unwrap/to_wavelength.py:482: in _compute_wavelength_data
data = _compute_wavelength_histogram(da=da, lookup=lookup, ltotal=ltotal)
^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
.tox/nexusfiles-scipp-odin/lib/python3.12/site-packages/ess/reduce/unwrap/to_wavelength.py:162: in _compute_wavelength_histogram
).rebin({key: new_bins})
^^^^^^^^^^^^^^^^^^^^^^
.tox/nexusfiles-scipp-odin/lib/python3.12/site-packages/scipp/core/data_group.py:755: in impl
return func(data, *args, **kwargs)
^^^^^^^^^^^^^^^^^^^^^^^^^^^
_ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _
x = <scipp.DataArray>
Dimensions: Sizes[time:10066, frame_time:7000, ]
Coordinates:
* frame_time float64 ...
Data:
float64 [au] (time, frame_time) [0, 0, ..., 0, 0] [0, 0, ..., 0, 0]
arg_dict = {'frame_time': <scipp.Variable> (frame_time: 7002) float64 [µs] [-2.14405e+09, -2.13497e+09, ..., 2.07927e+09, 2.14566e+09]}
kwargs = {}
@data_group_overload
def rebin(
x: Variable | DataArray | Dataset | DataGroup[Any],
arg_dict: IntoStrDict[SupportsIndex | Variable] | None = None,
/,
**kwargs: SupportsIndex | Variable,
) -> Variable | DataArray | Dataset | DataGroup[Any]:
"""Rebin a data array or dataset.
The coordinate of the input for the dimension to be rebinned must contain bin edges,
i.e., the data must be histogrammed.
If the input has masks that contain the dimension being rebinned then those
masks are applied to the data before rebinning. That is, masked values are treated
as zero.
Parameters
----------
x:
Data to rebin.
arg_dict:
Dictionary mapping dimension labels to binning parameters.
**kwargs:
Mapping of dimension label to corresponding binning parameters.
Returns
-------
:
Data rebinned according to the new bin edges.
See Also
--------
scipp.bin:
For changing the binning of binned (as opposed to dense, histogrammed) data.
scipp.hist:
For histogramming data.
Examples
--------
Rebin a data array along one of its dimensions, specifying (1) number of bins, (2)
bin width, or (3) actual binning:
>>> from numpy.random import default_rng
>>> rng = default_rng(seed=1234)
>>> x = sc.array(dims=['row'], unit='m', values=rng.random(100))
>>> y = sc.array(dims=['row'], unit='m', values=rng.random(100))
>>> data = sc.ones(dims=['row'], unit='K', shape=[100])
>>> table = sc.DataArray(data=data, coords={'x': x, 'y': y})
>>> da = table.hist(x=100, y=100)
>>> da.rebin(x=2).sizes
{'x': 2, 'y': 100}
>>> da.rebin(x=sc.scalar(0.2, unit='m')).sizes
{'x': 5, 'y': 100}
>>> da.rebin(x=sc.linspace('x', 0.2, 0.8, num=10, unit='m')).sizes
{'x': 9, 'y': 100}
Rebin a data array along two of its dimensions:
>>> da = table.hist(x=100, y=100)
>>> da.rebin(x=4, y=6).sizes
{'x': 4, 'y': 6}
"""
if isinstance(x, DataGroup):
# Only to make mypy happy because we have `DataGroup` in annotation of `x`
# so that Sphinx shows it.
raise TypeError("Internal error: input should not be a DataGroup")
edges = _make_edges(x, arg_dict, kwargs)
out = x
for dim, edge in edges.items():
> out = _cpp.rebin(out, dim, edge)
^^^^^^^^^^^^^^^^^^^^^^^^^^
E scipp._scipp.core.BinEdgeError: The input does not have coordinates with bin-edges.
E In provider ess.reduce.unwrap.to_wavelength.monitor_wavelength_data(
E monitor_data: ess.reduce.nexus.types.RawMonitor[ess.imaging.types.SampleRun, ess.imaging.types.BeamMonitor1],
E lookup: ess.reduce.unwrap.types.ErrorLimitedLookupTable[ess.imaging.types.SampleRun, ess.imaging.types.BeamMonitor1],
E ltotal: ess.reduce.unwrap.types.MonitorLtotal[ess.imaging.types.SampleRun, ess.imaging.types.BeamMonitor1],
E pulse_stride_offset: ess.reduce.unwrap.types.PulseStrideOffset,
E keep_event_time_offset: ess.reduce.unwrap.types.KeepEventTimeOffset,
E ) -> ess.reduce.unwrap.types.WavelengthMonitor[ess.imaging.types.SampleRun, ess.imaging.types.BeamMonitor1]
.tox/nexusfiles-scipp-odin/lib/python3.12/site-packages/scipp/core/binning.py:992: BinEdgeError
