source_folder_positions = '../processed_data/20180628_Sample_6_bead_1/position_data/'
method = 'fit_ellipse'
tag = 'Sample_6_Bead_1'
particle_diameter = 45 # um
source_folder_positions = '../processed_data/20180806_Sample_9_Bead_2/position_data/'
method = 'fit_ellipse'
tag = 'Sample_9_Bead_2'
particle_diameter = 31 # um
# get all the files and sort them by the run number
# position_file_names = get_position_file_names(source_folder_positions, method=method, runs=list(range(60, 180)))
# position_file_names = get_position_file_names(source_folder_positions, method=method, runs=list(range(83, 85)))
position_file_names = get_position_file_names(source_folder_positions,
method=method,
tag=tag,
runs=list(range(50, 54)))
################################################################################
#### run the script
################################################################################
for i, filename in enumerate(position_file_names):
print(filename)
data, info = load_time_trace(
filename,
source_folder_positions=source_folder_positions,
verbose=False)
return_dict = get_calibration_factor(data,
particle_diameter=particle_diameter,
print_only_names = True
print_only_names = False
################################################################################
#### experiment specific settings
################################################################################
if experiment == 'long term run 2a':
interval_width = 200
interval_width_zoom = 0.9
fo = 101
analysis_method = 1
mode = 'r'
# get all the files and sort them by the run number
position_file_names = get_position_file_names(source_folder_positions,
method=method)
# select the subset of interest
position_file_names = position_file_names[29:] # all
# position_file_names = position_file_names[29:32]
# position_file_names = position_file_names[32:35]
elif experiment == 'long term run r-mode':
interval_width = None
interval_width_zoom = 0.9
fo = None
analysis_method = 2
mode = 'r'
# get all the files and sort them by the run number
method = 'fit_ellipse'
n_smooth = None
################################################################################
#### experiment specific settings
################################################################################
if experiment == 'long term run 2 x':
interval_width = 5
interval_width_zoom = 0.9
fo = 60
mode = 'x'
analysis_method = 1
# get all the files and sort them by the run number
position_file_names = get_position_file_names(source_folder_positions,
method=method)
# select the subset of interest
position_file_names = position_file_names[29:] # all
elif experiment == 'long term run 2 y':
interval_width = 5
interval_width_zoom = 0.9
fo = 70
mode = 'y'
analysis_method = 1
# get all the files and sort them by the run number
position_file_names = get_position_file_names(source_folder_positions,
method=method)
# select the subset of interest
position_file_names = position_file_names[29:] # all
from glob import glob
import matplotlib.pyplot as plt
source_folder_positions = '../processed_data/20180607_Sample_6_bead_1/position_data/'
method = 'fit_ellipse'
image_type = 'rotation vs time'
################################################################################
#### image_type specific settings
################################################################################
if image_type == 'ellipse_spectra_zoom':
position_file_names = get_position_file_names(source_folder_positions, method=method)
image_folder = '../images/'
annotation_dict = {'2r': [13.16877399924, 0.6],
'm': [544.131481993, 0.6],
'r': [6.58438699962, 0.85],
'x': [171.52045677, 0.85],
'y': [177.690573472, 0.6],
'z': [402.859033756, 0.6]}
freq_window = 0.5
method = 'fit_ellipse'
plot_type = 'lin'
elif image_type == 'rotation vs time':
position_file_names = get_position_file_names(source_folder_positions, method=method)
# settings - run 73
run = 73
mode = 'x'
window_length = 5000
frequency_window = 20
source_folder_positions = '../processed_data/20180628_Sample_6_Bead_1/position_data/'
method = 'fit_ellipse'
magnet_diameter = 45
image_folder = '../images/20180628_Sample_6_Bead_1/ring-down/'
################################################################################
#### run the script
################################################################################
filename = get_position_file_names(source_folder_positions,
method=method,
runs=[run])[0]
print(filename)
data, info = load_time_trace(filename,
source_folder_positions=source_folder_positions,
verbose=False)
x_pos = data['ellipse ' + mode]
calib = load_info(filename=filename, folder_positions=source_folder_positions
)['ellipse']['calibration factor: (um/px)']
fo = load_info(filename=filename,
folder_positions=source_folder_positions)['ellipse'][mode]
fs = load_info(filename=filename,
dataset = 'relaxation_run5g'
dataset = 'relaxation_run6'
dataset = 'relaxation_run7a'
dataset = 'relaxation_run7b'
dataset = 'relaxation_run7d'
################################################################################
## end settings ###
################################################################################
################################################################################
#### for real data: 20180607_Sample6_bead_1
################################################################################
if dataset == 'relaxation_run5a':
position_file_names = get_position_file_names(source_folder_positions,
method=method,
runs=list(range(118, 123)))
modes = 'xyr'
elif dataset == 'relaxation_run5b':
position_file_names = get_position_file_names(source_folder_positions,
method=method,
runs=list(range(123, 128)))
modes = 'xyr'
elif dataset == 'relaxation_run5c':
position_file_names = get_position_file_names(source_folder_positions,
method=method,
runs=list(range(128, 132)))
modes = 'xyr'
elif dataset == 'relaxation_run5d':
position_file_names = get_position_file_names(source_folder_positions,
method=method,