def focus_run(self, sample_paths, banks, plot_output, instrument, rb_num,
spectrum_numbers):
"""
Focus some data using the current calibration.
:param sample_paths: The paths to the data to be focused.
:param banks: The banks that should be focused.
:param plot_output: True if the output should be plotted.
:param instrument: The instrument that the data came from.
:param rb_num: The experiment number, used to create directories. Can be None
:param spectrum_numbers: The specific spectra that should be focused. Used instead of banks.
"""
if not Ads.doesExist(vanadium_corrections.INTEGRATED_WORKSPACE_NAME
) and not Ads.doesExist(
vanadium_corrections.CURVES_WORKSPACE_NAME):
return
integration_workspace = Ads.retrieve(
vanadium_corrections.INTEGRATED_WORKSPACE_NAME)
curves_workspace = Ads.retrieve(
vanadium_corrections.CURVES_WORKSPACE_NAME)
output_workspaces = [] # List of collated workspaces to plot.
full_calib_path = get_setting(path_handling.INTERFACES_SETTINGS_GROUP,
path_handling.ENGINEERING_PREFIX,
"full_calibration")
if full_calib_path is not None and path.exists(full_calib_path):
full_calib_workspace = LoadAscii(full_calib_path,
OutputWorkspace="det_pos",
Separator="Tab")
else:
full_calib_workspace = None
if spectrum_numbers is None:
for sample_path in sample_paths:
sample_workspace = path_handling.load_workspace(sample_path)
run_no = path_handling.get_run_number_from_path(
sample_path, instrument)
workspaces_for_run = []
for name in banks:
output_workspace_name = str(
run_no) + "_" + FOCUSED_OUTPUT_WORKSPACE_NAME + str(
name)
self._run_focus(sample_workspace, output_workspace_name,
integration_workspace, curves_workspace,
name, full_calib_workspace)
workspaces_for_run.append(output_workspace_name)
# Save the output to the file system.
self._save_output(instrument, sample_path, name,
output_workspace_name, rb_num)
output_workspaces.append(workspaces_for_run)
else:
for sample_path in sample_paths:
sample_workspace = path_handling.load_workspace(sample_path)
run_no = path_handling.get_run_number_from_path(
sample_path, instrument)
output_workspace_name = str(
run_no) + "_" + FOCUSED_OUTPUT_WORKSPACE_NAME + "cropped"
self._run_focus(sample_workspace, output_workspace_name,
integration_workspace, curves_workspace, None,
full_calib_workspace, spectrum_numbers)
output_workspaces.append([output_workspace_name])
self._save_output(instrument, sample_path, "cropped",
output_workspace_name, rb_num)
# Plot the output
if plot_output:
for ws_names in output_workspaces:
self._plot_focused_workspaces(ws_names)
def create_new_calibration(self,
vanadium_path,
sample_path,
plot_output,
instrument,
rb_num=None,
bank=None,
spectrum_numbers=None):
"""
Create a new calibration from a vanadium run and sample run
:param vanadium_path: Path to vanadium data file.
:param sample_path: Path to sample (CeO2) data file
:param plot_output: Whether the output should be plotted.
:param instrument: The instrument the data relates to.
:param rb_num: The RB number for file creation.
:param bank: Optional parameter to crop by bank
:param spectrum_numbers: Optional parameter to crop using spectrum numbers.
"""
van_integration, van_curves = vanadium_corrections.fetch_correction_workspaces(
vanadium_path, instrument, rb_num=rb_num)
sample_workspace = path_handling.load_workspace(sample_path)
full_calib_path = get_setting(path_handling.INTERFACES_SETTINGS_GROUP,
path_handling.ENGINEERING_PREFIX, "full_calibration")
if full_calib_path is not None and path.exists(full_calib_path):
full_calib = LoadAscii(full_calib_path, OutputWorkspace="det_pos", Separator="Tab")
output = self.run_calibration(sample_workspace,
van_integration,
van_curves,
bank,
spectrum_numbers,
full_calib_ws=full_calib)
else:
output = self.run_calibration(sample_workspace, van_integration, van_curves, bank,
spectrum_numbers)
if plot_output:
self._plot_vanadium_curves()
for i in range(len(output)):
if spectrum_numbers:
bank_name = "cropped"
elif bank is None:
bank_name = str(i + 1)
else:
bank_name = bank
difa = output[i].DIFA
difc = output[i].DIFC
tzero = output[i].TZERO
self._generate_tof_fit_workspace(difa, difc, tzero, bank_name)
if bank is None and spectrum_numbers is None:
self._plot_tof_fit()
elif spectrum_numbers is None:
self._plot_tof_fit_single_bank_or_custom(bank)
else:
self._plot_tof_fit_single_bank_or_custom("cropped")
difa = [i.DIFC for i in output]
difc = [i.DIFC for i in output]
tzero = [i.TZERO for i in output]
params_table = []
for i in range(len(difc)):
params_table.append([i, difc[i], difa[i], tzero[i]])
self.update_calibration_params_table(params_table)
calib_dir = path.join(path_handling.get_output_path(), "Calibration", "")
self.create_output_files(calib_dir, difa, difc, tzero, sample_path, vanadium_path, instrument,
bank, spectrum_numbers)
if rb_num:
user_calib_dir = path.join(path_handling.get_output_path(), "User", rb_num,
"Calibration", "")
self.create_output_files(user_calib_dir, difa, difc, tzero, sample_path, vanadium_path,
instrument, bank, spectrum_numbers)