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_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 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()