TBL: Orca image normalization workflow#

This notebook shows how to use the workflow to compute normalized images recorded by the Orca detector on the TBL instrument.

[1]:

import ess.tbl.data  # noqa: F401
from ess import tbl
from ess.imaging.types import *
import scipp as sc
import plopp as pp

%matplotlib widget

Workflow setup#

[2]:

wf = tbl.OrcaNormalizedImagesWorkflow()

wf[Filename[SampleRun]] = tbl.data.tbl_lego_sample_run()
wf[Filename[DarkBackgroundRun]] = tbl.data.tbl_lego_dark_run()
wf[Filename[OpenBeamRun]] = tbl.data.tbl_lego_openbeam_run()
wf[NeXusDetectorName] = 'orca_detector'

wf[MaskingRules] = {}  # No masks to begin with
wf[UncertaintyBroadcastMode] = UncertaintyBroadcastMode.upper_bound

Downloading file 'ymir_lego_sample_run.hdf' from 'https://public.esss.dk/groups/scipp/ess/tbl/3/ymir_lego_sample_run.hdf' to '/home/runner/.cache/ess/tbl'.
Downloading file 'ymir_lego_dark_run.hdf' from 'https://public.esss.dk/groups/scipp/ess/tbl/3/ymir_lego_dark_run.hdf' to '/home/runner/.cache/ess/tbl'.
Downloading file 'ymir_lego_openbeam_run.hdf' from 'https://public.esss.dk/groups/scipp/ess/tbl/3/ymir_lego_openbeam_run.hdf' to '/home/runner/.cache/ess/tbl'.

[3]:

wf.visualize(NormalizedImage, compact=True, graph_attr={"rankdir": "LR"})

[3]:

../_images/tbl_orca-image-normalization_4_0.svg

Adding masks#

If we want to mask some part of the image, we update the masking rules. For example, here we mask the upper part of the image:

[6]:

wf[MaskingRules] = {'y_pixel_offset': lambda x: x > sc.scalar(0.082, unit='m')}

pp.slicer(wf.compute(NormalizedImage), autoscale=False)

[6]:

Intermediate results#

We can also inspect intermediate results, which is useful for debugging:

[7]:

results = wf.compute([
    RawDetector[SampleRun],
    CorrectedDetector[SampleRun],
    BackgroundSubtractedDetector[SampleRun]
])

fig = pp.tiled(2, 2, hspace=0.3, wspace=0.3)
fig[0, 0] = results[RawDetector[SampleRun]]['time', 0].plot(title='Raw data')
fig[0, 1] = results[CorrectedDetector[SampleRun]]['time', 0].plot(title='Masks applied')
fig[1, 0] = results[BackgroundSubtractedDetector[SampleRun]]['time', 0].plot(title='Background subtracted')
fig[1, 1] = image['time', 0].plot(title='Final image')
fig

[7]: