print(get_frame_shifts(frequency_wfm_psc=70))
print("Frame shifts are:", frame_shifts)
# Get frame windows from sample
rebinned = Rebin(ws, rebin_parameters, PreserveEvents=False)
summed = SumSpectra(rebinned)
nz = summed.blocksize()
summed_data = np.zeros([nz, 2])
summed_x = summed.extractX()[0]
summed_y = summed.extractY()[0]
for i in range(nz):
summed_data[i, 0] = 0.5 * (summed_x[i] + summed_x[i + 1])
summed_data[i, 1] = summed_y[i]
frame_parameters1 = get_wfm_windows(data=summed_data,
output_format="tof",
plot=True,
win_threshold=0.30)
print(frame_parameters1)
# Get frame windows from background
rebinned_van = Rebin(van, rebin_parameters, PreserveEvents=False)
summed_van = SumSpectra(rebinned_van)
nz = summed_van.blocksize()
summed_data_van = np.zeros([nz, 2])
summed_x_van = summed_van.extractX()[0]
summed_y_van = summed_van.extractY()[0]
for i in range(nz):
summed_data_van[i, 0] = 0.5 * (summed_x_van[i] + summed_x_van[i + 1])
summed_data_van[i, 1] = summed_y_van[i]
frame_parameters2 = get_wfm_windows(data=summed_data_van,
print(get_frame_shifts(frequency_wfm_psc=70))
print("Frame shifts are:", frame_shifts)
# rebinned = Rebin(ws, binsize, PreserveEvents=False)
rebinned = Rebin(ws, rebin_parameters, PreserveEvents=False)
summed = SumSpectra(rebinned)
nz = summed.blocksize()
summed_data = np.zeros([nz, 2])
summed_x = summed.extractX()[0]
summed_y = summed.extractY()[0]
for i in range(nz):
summed_data[i, 0] = 0.5 * (summed_x[i] + summed_x[i + 1])
summed_data[i, 1] = summed_y[i]
frame_parameters = get_wfm_windows(data=summed_data,
rebin_step_for_string_output=64,
plot=True,
win_threshold=0.30)
print(frame_parameters)
# rebin region reduced to cut off frames that contain no data
# rebin_parameters = '8500,64,43000'
# rebin_parameters = '{},64,{}'.format(frame_parameters[0].split(',')[0],frame_parameters[-1].split(',')[2])
# rebin_parameters = '{},64,{}'.format(tofs_lims[0],tofs_lims[1])
rebin_parameters = '{},64,{}'.format(tof_min, tof_max)
instrument_filename = "V20_geometry.nxs" # theFile # os.path.join(grand_parent_dir, 'IDF', 'V20_Definition_GP2.xml')
# mapping_file_to_write = os.path.join(this_dir, 'mapping_file.txt')
# defining grouping of 2D detector pixels
# grouping_number = 3
# nx_target = grouping_number
# ny_target = grouping_number