SANS2D: I(Q) for sample and background#

In this notebook, we will be reducing a sample and a background measurements to a one-dimensional \(I(Q)\).

It assumes the detector data has been recorded in event mode, while the monitor data has been histogrammed.

The data used in this notebook has been published in Manasi et al. (2021), and we kindly thank the authors for allowing us to use their data.

Outline:

We will begin by loading the data files containing the sample, direct, and background measurements.
We will then apply some corrections to beamline components specific to the SANS2D beamline.
This will be followed by some masking of some saturated or defect detector pixels
Finally, the sample and background measurement will be converted to the \(Q\) dimension

[1]:

import matplotlib.pyplot as plt
import scipp as sc
from ess import loki, sans
from ess.logging import configure_workflow
import scippneutron as scn

[2]:

logger = configure_workflow('sans2d_reduction', filename='sans2d.log')

Define reduction workflow parameters#

We define here whether to include the effects of gravity, as well as common time-of-flight, wavelength and \(Q\) bins for all the measurements.

We also define a range of wavelengths for the monitors that are considered to not be part of the background.

[3]:

# Include effects of gravity?
gravity = True

tof_bins = sc.linspace(dim='tof', start=0, stop=100000, num=2, unit='us')

wavelength_bins = sc.linspace(dim='wavelength', start=2.0, stop=16.0, num=141,
                              unit='angstrom')

q_bins = sc.linspace(dim='Q', start=0.01, stop=0.6, num=141, unit='1/angstrom')

# Define the range of wavelengths for the monitors that are considered
# to not be part of the background
monitor_non_background_range = sc.array(dims=['wavelength'],
                                        values=[0.7, 17.1], unit='angstrom')

Apply corrections to pixel positions#

We apply some corrections to the detector pixel and monitor positions, as the geometry stored in the file is inaccurate.

[5]:

# Custom SANS2D position offsets
sample_pos_z_offset = 0.053 * sc.units.m
bench_pos_y_offset = 0.001 * sc.units.m
# There is some uncertainity here
monitor4_pos_z_offset = -6.719 * sc.units.m

# Geometry transformation
x_offset = -0.09288 * sc.units.m
y_offset = 0.08195 * sc.units.m

[6]:

ds.coords["pixel_width"] = 0.0035 * sc.units.m
ds.coords["pixel_height"] = 0.002033984375 * sc.units.m

# Change sample position
ds.coords["sample_position"].fields.z += sample_pos_z_offset
# Apply bench offset to pixel positions
ds.coords["position"].fields.y += bench_pos_y_offset

for key in ds:
    ds[key].attrs["monitor4"].value.coords["position"].fields.z += monitor4_pos_z_offset

# Now shift pixels positions to get the correct beam center
ds.coords['position'].fields.x += x_offset
ds.coords['position'].fields.y += y_offset

References#

Manasi I., Andalibi M. R., Atri R. S., Hooton J., King S. M., Edler K. J., 2021, Self-assembly of ionic and non-ionic surfactants in type IV cerium nitrate and urea based deep eutectic solvent, J. Chem. Phys. 155, 084902

SANS2D: I(Q) for sample and background

Contents

SANS2D: I(Q) for sample and background#

Define reduction workflow parameters#

Loading data files#

Apply corrections to pixel positions#

Masking#

Mask bad pixels#

Mask Bragg peaks in time-of-flight#

Use to_I_of_Q workflow#

Wavelength bands#

References#