def test_that_divide_uses_settings_from_state_if_they_are_set(self):
# Arrange
facility = SANSFacility.ISIS
data_builder = get_data_builder(facility)
data_builder.set_sample_scatter("SANS2D00022024")
data_state = data_builder.build()
scale_builder = get_scale_builder(data_state)
width = 10.
height = 5.
thickness = 2.
scale_builder.set_shape(SampleShape.CylinderAxisAlong)
scale_builder.set_thickness(thickness)
scale_builder.set_width(width)
scale_builder.set_height(height)
scale_state = scale_builder.build()
test_director = TestDirector()
test_director.set_states(scale_state=scale_state,
data_state=data_state)
state = test_director.construct()
divide_factory = DivideByVolumeFactory()
divider = divide_factory.create_divide_by_volume(state)
workspace = ScaleHelperTest._get_workspace()
# Act
output_workspace = divider.divide_by_volume(workspace, scale_state)
# Assert
expected_volume = thickness * math.pi * math.pow(width, 2) / 4.0
expected_value = 0.3 / expected_volume
data_y = output_workspace.dataY(0)
self.assertEqual(data_y[0], expected_value)
def test_that_correct_scale_strategy_is_selected_for_loq_2(self):
# Arrange
facility = SANSFacility.ISIS
data_builder = get_data_builder(facility)
data_builder.set_sample_scatter("LOQ74044")
data_state = data_builder.build()
scale_builder = get_scale_builder(data_state)
scale_builder.set_scale(2.4)
scale_state = scale_builder.build()
test_director = TestDirector()
test_director.set_states(scale_state=scale_state,
data_state=data_state)
state_loq = test_director.construct()
absolute_multiply_factory = MultiplyByAbsoluteScaleFactory()
multiplier = absolute_multiply_factory.create_multiply_by_absolute(
state_loq)
workspace = self._get_workspace()
# Act
output_workspace = multiplier.multiply_by_absolute_scale(
workspace, state_loq.scale)
# Assert
expected_value = 0.3 * 2.4 / math.pi * 100.
data_y = output_workspace.dataY(0)
self.assertEqual(data_y[0], expected_value)
def test_that_divide_uses_settings_from_workspace(self):
# Arrange
facility = SANSFacility.ISIS
data_builder = get_data_builder(facility)
data_builder.set_sample_scatter("SANS2D00022024")
data_state = data_builder.build()
scale_builder = get_scale_builder(data_state)
scale_state = scale_builder.build()
test_director = TestDirector()
test_director.set_states(scale_state=scale_state,
data_state=data_state)
state = test_director.construct()
divide_factory = DivideByVolumeFactory()
divider = divide_factory.create_divide_by_volume(state)
# Apply the settings from the SANS2D00022024 workspace
width = 8.
height = 8.
thickness = 1.
shape = 3
workspace = ScaleHelperTest._get_workspace(width, height, thickness,
shape)
# Act
output_workspace = divider.divide_by_volume(workspace, scale_state)
# Assert
expected_volume = thickness * math.pi * math.pow(width, 2) / 4.0
expected_value = 0.3 / expected_volume
data_y = output_workspace.dataY(0)
self.assertEqual(data_y[0], expected_value)
def test_that_divide_uses_settings_from_workspace(self):
# Arrange
facility = SANSFacility.ISIS
data_builder = get_data_builder(facility)
data_builder.set_sample_scatter("SANS2D00022024")
data_state = data_builder.build()
scale_builder = get_scale_builder(data_state)
scale_state = scale_builder.build()
test_director = TestDirector()
test_director.set_states(scale_state=scale_state, data_state=data_state)
state = test_director.construct()
divide_factory = DivideByVolumeFactory()
divider = divide_factory.create_divide_by_volume(state)
# Apply the settings from the SANS2D00022024 workspace
width = 8.
height = 8.
thickness = 1.
shape = 3
workspace = ScaleHelperTest._get_workspace(width, height, thickness, shape)
# Act
output_workspace = divider.divide_by_volume(workspace, scale_state)
# Assert
expected_volume = thickness * math.pi * math.pow(width, 2) / 4.0
expected_value = 0.3 / expected_volume
data_y = output_workspace.dataY(0)
self.assertEqual(data_y[0], expected_value)
def test_that_divide_uses_settings_from_state_if_they_are_set(self):
# Arrange
facility = SANSFacility.ISIS
data_builder = get_data_builder(facility)
data_builder.set_sample_scatter("SANS2D00022024")
data_state = data_builder.build()
scale_builder = get_scale_builder(data_state)
width = 10.
height = 5.
thickness = 2.
scale_builder.set_shape(SampleShape.CylinderAxisAlong)
scale_builder.set_thickness(thickness)
scale_builder.set_width(width)
scale_builder.set_height(height)
scale_state = scale_builder.build()
test_director = TestDirector()
test_director.set_states(scale_state=scale_state, data_state=data_state)
state = test_director.construct()
divide_factory = DivideByVolumeFactory()
divider = divide_factory.create_divide_by_volume(state)
workspace = ScaleHelperTest._get_workspace()
# Act
output_workspace = divider.divide_by_volume(workspace, scale_state)
# Assert
expected_volume = thickness * math.pi * math.pow(width, 2) / 4.0
expected_value = 0.3 / expected_volume
data_y = output_workspace.dataY(0)
self.assertEqual(data_y[0], expected_value)
def test_that_slice_event_state_can_be_built(self):
# Arrange
facility = SANSFacility.ISIS
file_information = SANSFileInformationMock(run_number=74044)
data_builder = get_data_builder(facility, file_information)
data_builder.set_sample_scatter("LOQ74044")
data_info = data_builder.build()
# Act
builder = get_scale_builder(data_info, file_information)
self.assertTrue(builder)
builder.set_scale(1.0)
builder.set_shape(SampleShape.Cuboid)
builder.set_thickness(3.6)
builder.set_width(3.7)
builder.set_height(5.8)
# Assert
state = builder.build()
self.assertTrue(state.shape is SampleShape.Cuboid)
self.assertTrue(state.scale == 1.0)
self.assertTrue(state.thickness == 3.6)
self.assertTrue(state.width == 3.7)
self.assertTrue(state.height == 5.8)
def test_that_correct_scale_strategy_is_selected_for_loq_2(self):
# Arrange
facility = SANSFacility.ISIS
data_builder = get_data_builder(facility)
data_builder.set_sample_scatter("LOQ74044")
data_state = data_builder.build()
scale_builder = get_scale_builder(data_state)
scale_builder.set_scale(2.4)
scale_state = scale_builder.build()
test_director = TestDirector()
test_director.set_states(scale_state=scale_state, data_state=data_state)
state_loq = test_director.construct()
absolute_multiply_factory = MultiplyByAbsoluteScaleFactory()
multiplier = absolute_multiply_factory.create_multiply_by_absolute(state_loq)
workspace = self._get_workspace()
# Act
output_workspace = multiplier.multiply_by_absolute_scale(workspace, state_loq.scale)
# Assert
expected_value = 0.3 * 2.4 / math.pi * 100.
data_y = output_workspace.dataY(0)
self.assertEqual(data_y[0], expected_value)
def test_that_slice_event_state_can_be_built(self):
# Arrange
facility = SANSFacility.ISIS
file_information = SANSFileInformationMock(run_number=74044)
data_builder = get_data_builder(facility, file_information)
data_builder.set_sample_scatter("LOQ74044")
data_info = data_builder.build()
# Act
builder = get_scale_builder(data_info, file_information)
self.assertTrue(builder)
builder.set_scale(1.0)
builder.set_shape(SampleShape.FlatPlate)
builder.set_thickness(3.6)
builder.set_width(3.7)
builder.set_height(5.8)
# Assert
state = builder.build()
self.assertTrue(state.shape is SampleShape.FlatPlate)
self.assertTrue(state.scale == 1.0)
self.assertTrue(state.thickness == 3.6)
self.assertTrue(state.width == 3.7)
self.assertTrue(state.height == 5.8)
def _get_sample_state(width, height, thickness, shape, scale):
# Arrange
facility = SANSFacility.ISIS
data_builder = get_data_builder(facility)
data_builder.set_sample_scatter("LOQ74044")
data_state = data_builder.build()
scale_builder = get_scale_builder(data_state)
scale_builder.set_scale(scale)
scale_builder.set_thickness(thickness)
scale_builder.set_width(width)
scale_builder.set_height(height)
scale_builder.set_shape(shape)
scale_state = scale_builder.build()
test_director = TestDirector()
test_director.set_states(scale_state=scale_state, data_state=data_state)
return test_director.construct()
def _get_sample_state(width, height, thickness, shape, scale):
# Arrange
facility = SANSFacility.ISIS
file_information = SANSFileInformationMock(instrument=SANSInstrument.LOQ, run_number=74044)
data_builder = get_data_builder(facility, file_information)
data_builder.set_sample_scatter("LOQ74044")
data_state = data_builder.build()
scale_builder = get_scale_builder(data_state, file_information)
scale_builder.set_scale(scale)
scale_builder.set_thickness(thickness)
scale_builder.set_width(width)
scale_builder.set_height(height)
scale_builder.set_shape(shape)
scale_state = scale_builder.build()
test_director = TestDirector()
test_director.set_states(scale_state=scale_state, data_state=data_state)
return test_director.construct()
def _get_sample_state(width, height, thickness, shape, scale):
# Arrange
facility = SANSFacility.ISIS
file_information = SANSFileInformationMock(instrument=SANSInstrument.LOQ, run_number=74044)
data_builder = get_data_builder(facility, file_information)
data_builder.set_sample_scatter("LOQ74044")
data_state = data_builder.build()
scale_builder = get_scale_builder(data_state, file_information)
scale_builder.set_scale(scale)
scale_builder.set_thickness(thickness)
scale_builder.set_width(width)
scale_builder.set_height(height)
scale_builder.set_shape(shape)
scale_state = scale_builder.build()
test_director = TestDirector()
test_director.set_states(scale_state=scale_state, data_state=data_state)
return test_director.construct()
def _get_sample_state(width, height, thickness, shape, scale):
# Arrange
facility = SANSFacility.ISIS
data_builder = get_data_builder(facility)
data_builder.set_sample_scatter("LOQ74044")
data_state = data_builder.build()
scale_builder = get_scale_builder(data_state)
scale_builder.set_scale(scale)
scale_builder.set_thickness(thickness)
scale_builder.set_width(width)
scale_builder.set_height(height)
scale_builder.set_shape(shape)
scale_state = scale_builder.build()
test_director = TestDirector()
test_director.set_states(scale_state=scale_state,
data_state=data_state)
return test_director.construct()
def test_that_correct_scale_strategy_is_selected_for_loq(self):
# Arrange
facility = SANSFacility.ISIS
data_builder = get_data_builder(facility)
data_builder.set_sample_scatter("LOQ74044")
data_state = data_builder.build()
scale_builder = get_scale_builder(data_state)
scale_state = scale_builder.build()
test_director = TestDirector()
test_director.set_states(scale_state=scale_state, data_state=data_state)
state_loq = test_director.construct()
absolute_multiply_factory = MultiplyByAbsoluteScaleFactory()
# Act
multiplier = absolute_multiply_factory.create_multiply_by_absolute(state_loq)
# Assert
self.assertTrue(isinstance(multiplier, MultiplyByAbsoluteScaleLOQ))
def test_that_correct_scale_strategy_is_selected_for_loq(self):
# Arrange
facility = SANSFacility.ISIS
file_information = SANSFileInformationMock(instrument=SANSInstrument.LOQ, run_number=74044)
data_builder = get_data_builder(facility, file_information)
data_builder.set_sample_scatter("LOQ74044")
data_state = data_builder.build()
scale_builder = get_scale_builder(data_state, file_information)
scale_state = scale_builder.build()
test_director = TestDirector()
test_director.set_states(scale_state=scale_state, data_state=data_state)
state_loq = test_director.construct()
absolute_multiply_factory = MultiplyByAbsoluteScaleFactory()
# Act
multiplier = absolute_multiply_factory.create_multiply_by_absolute(state_loq)
# Assert
self.assertTrue(isinstance(multiplier, MultiplyByAbsoluteScaleLOQ))
def test_that_slice_event_state_can_be_built(self):
# Arrange
facility = SANSFacility.ISIS
data_builder = get_data_builder(facility)
data_builder.set_sample_scatter("LOQ74044")
data_info = data_builder.build()
# Act
builder = get_scale_builder(data_info)
self.assertTrue(builder)
builder.set_scale(1.0)
builder.set_shape(SampleShape.Cuboid)
builder.set_thickness(3.6)
builder.set_width(3.7)
builder.set_height(5.8)
# Assert
state = builder.build()
self.assertTrue(state.shape is SampleShape.Cuboid)
self.assertTrue(state.scale == 1.0)
self.assertTrue(state.thickness == 3.6)
self.assertTrue(state.width == 3.7)
self.assertTrue(state.height == 5.8)
def construct(self):
facility = SANSFacility.ISIS
# Build the SANSStateData
if self.data_state is None:
data_builder = get_data_builder(facility)
data_builder.set_sample_scatter("SANS2D00022024")
data_builder.set_can_scatter("SANS2D00022024")
self.data_state = data_builder.build()
# Build the SANSStateMove
if self.move_state is None:
move_builder = get_move_builder(self.data_state)
if hasattr(move_builder, "set_HAB_x_translation_correction"):
move_builder.set_HAB_x_translation_correction(21.2)
move_builder.set_LAB_x_translation_correction(12.1)
self.move_state = move_builder.build()
# Build the SANSStateReduction
if self.reduction_state is None:
reduction_builder = get_reduction_mode_builder(self.data_state)
reduction_builder.set_reduction_mode(ISISReductionMode.Merged)
reduction_builder.set_reduction_dimensionality(ReductionDimensionality.OneDim)
reduction_builder.set_merge_fit_mode(FitModeForMerge.Both)
reduction_builder.set_merge_shift(324.2)
reduction_builder.set_merge_scale(3420.98)
self.reduction_state = reduction_builder.build()
# Build the SANSStateSliceEvent
if self.slice_state is None:
slice_builder = get_slice_event_builder(self.data_state)
slice_builder.set_start_time([0.1, 1.3])
slice_builder.set_end_time([0.2, 1.6])
self.slice_state = slice_builder.build()
# Build the SANSStateMask
if self.mask_state is None:
mask_builder = get_mask_builder(self.data_state)
mask_builder.set_radius_min(10.0)
mask_builder.set_radius_max(20.0)
self.mask_state = mask_builder.build()
# Build the SANSStateWavelength
if self.wavelength_state is None:
wavelength_builder = get_wavelength_builder(self.data_state)
wavelength_builder.set_wavelength_low(1.0)
wavelength_builder.set_wavelength_high(10.0)
wavelength_builder.set_wavelength_step(2.0)
wavelength_builder.set_wavelength_step_type(RangeStepType.Lin)
wavelength_builder.set_rebin_type(RebinType.Rebin)
self.wavelength_state = wavelength_builder.build()
# Build the SANSStateSave
if self.save_state is None:
save_builder = get_save_builder(self.data_state)
save_builder.set_user_specified_output_name("test_file_name")
save_builder.set_file_format([SaveType.Nexus])
self.save_state = save_builder.build()
# Build the SANSStateScale
if self.scale_state is None:
scale_builder = get_scale_builder(self.data_state)
scale_builder.set_shape(SampleShape.Cuboid)
scale_builder.set_width(1.0)
scale_builder.set_height(2.0)
scale_builder.set_thickness(3.0)
scale_builder.set_scale(4.0)
self.scale_state = scale_builder.build()
# Build the SANSAdjustmentState
if self.adjustment_state is None:
# NormalizeToMonitor
normalize_to_monitor_builder = get_normalize_to_monitor_builder(self.data_state)
normalize_to_monitor_builder.set_wavelength_low(1.0)
normalize_to_monitor_builder.set_wavelength_high(10.0)
normalize_to_monitor_builder.set_wavelength_step(2.0)
normalize_to_monitor_builder.set_wavelength_step_type(RangeStepType.Lin)
normalize_to_monitor_builder.set_rebin_type(RebinType.Rebin)
normalize_to_monitor_builder.set_background_TOF_general_start(1000.)
normalize_to_monitor_builder.set_background_TOF_general_stop(2000.)
normalize_to_monitor_builder.set_incident_monitor(1)
normalize_to_monitor = normalize_to_monitor_builder.build()
# CalculateTransmission
calculate_transmission_builder = get_calculate_transmission_builder(self.data_state)
calculate_transmission_builder.set_transmission_monitor(3)
calculate_transmission_builder.set_incident_monitor(2)
calculate_transmission_builder.set_wavelength_low(1.0)
calculate_transmission_builder.set_wavelength_high(10.0)
calculate_transmission_builder.set_wavelength_step(2.0)
calculate_transmission_builder.set_wavelength_step_type(RangeStepType.Lin)
calculate_transmission_builder.set_rebin_type(RebinType.Rebin)
calculate_transmission_builder.set_background_TOF_general_start(1000.)
calculate_transmission_builder.set_background_TOF_general_stop(2000.)
calculate_transmission_builder.set_Sample_fit_type(FitType.Linear)
calculate_transmission_builder.set_Sample_polynomial_order(0)
calculate_transmission_builder.set_Sample_wavelength_low(1.0)
calculate_transmission_builder.set_Sample_wavelength_high(10.0)
calculate_transmission_builder.set_Can_fit_type(FitType.Polynomial)
calculate_transmission_builder.set_Can_polynomial_order(3)
calculate_transmission_builder.set_Can_wavelength_low(10.0)
calculate_transmission_builder.set_Can_wavelength_high(20.0)
calculate_transmission = calculate_transmission_builder.build()
# Wavelength and pixel adjustment
wavelength_and_pixel_builder = get_wavelength_and_pixel_adjustment_builder(self.data_state)
wavelength_and_pixel_builder.set_wavelength_low(1.0)
wavelength_and_pixel_builder.set_wavelength_high(10.0)
wavelength_and_pixel_builder.set_wavelength_step(2.0)
wavelength_and_pixel_builder.set_wavelength_step_type(RangeStepType.Lin)
wavelength_and_pixel = wavelength_and_pixel_builder.build()
# Adjustment
adjustment_builder = get_adjustment_builder(self.data_state)
adjustment_builder.set_normalize_to_monitor(normalize_to_monitor)
adjustment_builder.set_calculate_transmission(calculate_transmission)
adjustment_builder.set_wavelength_and_pixel_adjustment(wavelength_and_pixel)
self.adjustment_state = adjustment_builder.build()
# SANSStateConvertToQ
if self.convert_to_q_state is None:
convert_to_q_builder = get_convert_to_q_builder(self.data_state)
convert_to_q_builder.set_reduction_dimensionality(ReductionDimensionality.OneDim)
convert_to_q_builder.set_use_gravity(False)
convert_to_q_builder.set_radius_cutoff(0.002)
convert_to_q_builder.set_wavelength_cutoff(12.)
convert_to_q_builder.set_q_min(0.1)
convert_to_q_builder.set_q_max(0.8)
convert_to_q_builder.set_q_1d_rebin_string("0.1,0.01,0.8")
convert_to_q_builder.set_use_q_resolution(False)
self.convert_to_q_state = convert_to_q_builder.build()
# Set the sub states on the SANSState
state_builder = get_state_builder(self.data_state)
state_builder.set_data(self.data_state)
state_builder.set_move(self.move_state)
state_builder.set_reduction(self.reduction_state)
state_builder.set_slice(self.slice_state)
state_builder.set_mask(self.mask_state)
state_builder.set_wavelength(self.wavelength_state)
state_builder.set_save(self.save_state)
state_builder.set_scale(self.scale_state)
state_builder.set_adjustment(self.adjustment_state)
state_builder.set_convert_to_q(self.convert_to_q_state)
return state_builder.build()
