import glob
import bead_util as bu
import tkinter
import tkinter.filedialog
import os, sys
from scipy.optimize import curve_fit
import bead_util as bu
from scipy.optimize import minimize_scalar as minimize
import pickle as pickle
import time
dirs = [
24,
]
ddict = bu.load_dir_file("/dirfiles/dir_file_july2017.txt")
#print ddict
calibrate = True
cal_drive_freq = 41.
respdir = 'X'
cant_axis = 2
lpf = 150 # Hz
bin_size = 1 # um
load_charge_cal = True
maxfiles = 100
fig_title = 'Force vs. Cantilever Position: Dipole Response'
import matplotlib.pyplot as plt import glob import bead_util as bu import tkinter import tkinter.filedialog import os, sys from scipy.optimize import curve_fit import bead_util as bu from scipy.optimize import minimize_scalar as minimize import pickle as pickle import time dirs = [45, 46] ddict = bu.load_dir_file( "/home/charles/opt_lev_analysis/scripts/dirfiles/dir_file_june2017.txt" ) files = [0,1] axes_to_plot = [2, 4] maxfiles = 2 load_charge_cal = True calibrate = True tf_path = './trans_funcs/Hout_20170627.p' step_cal_path = './calibrations/step_cal_20170627.p' thermal_cal_file_path = '/data/20170627/bead4/1_6mbar_nocool.h5' cal_drive_freq = 41. init_data = [0., 0., 20.]
limitdata_path2 = '/home/charles/opt_lev_analysis/scripts/gravity_sim/data/limitdata_20160914_datathief.txt'
limitdata2 = np.loadtxt(limitdata_path2, delimiter=',')
limitlab2 = 'With Decca 2'
init_data = [0, 0, 0]
derp_ind = 0
minval = -200.
stepval = 40.
### Load data directories
dirs = [2, 3, 4, 5, 6, 7, 8, 9]
ddict = bu.load_dir_file("/dirfiles/dir_file_aug2017.txt")
maxfiles = 3 # Maximum number of files to load from a directory
SWEEP_AX = 2 # Cantilever sweep axis, 1 for Y, 2 for Z
RESP_AX = 1 # Response Axis
bin_size = 1 # um, Binning for final force v. pos
lpf = 150 # Hz, acausal top-hat filter at this freq
cantfilt = True #True # Notch filter at cantilever drive
fig_title = 'Force vs. Cantilever Position:'
xlab = 'Distance along Cantilever [um]'
# Locate Calibration files
tf_path = '/calibrations/transfer_funcs/Hout_20170707.p'
step_cal_path = '/calibrations/step_cals/step_cal_20170707.p'
import bead_util as bu
from scipy.optimize import minimize_scalar as minimize
import pickle as pickle
plot_vs_time = True
plot_vs_sep = False
filstring = ''
filnames = ['/data/20170903/bead1/1_5mbar_nocool.h5', \
'/data/20170903/bead1/turbombar_xyzcool_discharged.h5']
labs = ['1.5 mbar', '1e-6 mbar']
use_labs = True
ddict = bu.load_dir_file('/dirfiles/dir_file_sept2017.txt')
dirs = []
chan_to_plot = [0, 1, 2]
chan_labs = ['X', 'Y', 'Z']
NFFT = 2**14
xlim = [1, 2500]
ylim = [6e-19, 1.5e-15]
maxfiles = 100
calibrate = True
tf_path = '/calibrations/transfer_funcs/Hout_20170903.p'
step_cal_path = '/calibrations/step_cals/step_cal_20170903.p'
import os, sys
from scipy.optimize import curve_fit
import bead_util as bu
from scipy.optimize import minimize_scalar as minimize
import pickle as pickle
import time
#dirs = [42,38,39,40,41]
#dirs = [99,100,101,102,103,104,107,108,109,110] # diff_speeds
#dirs = [113,114,115,116,117,118,119,120,121,122,123,124] # no bead
#dirs = [125,126,127] # Cant far background
#dirs = [128,129,130] # No red in
dirs = [212, 213]
ddict = bu.load_dir_file(
"/home/charles/opt_lev_classy/scripts/cant_force/dir_file.txt")
#print ddict
load_from_file = False
show_each_file = False
show_avg_force = False
fft = False
calibrate = True
respdir = 'Y'
cant_axis = 2
load_charge_cal = True
maxfiles = 1000
fig_title = 'Force vs. Cantilever Position: NO BEAD Background with Cantilever Far Away'