import h5py, os, glob, re
import RF, RR, Spikes, misc
from fileconversion import load_cfs
basedir = '/Volumes/BOB_SAGET/Fmr1_voc/voc_ko_nai_20130116'
experiment = 'voc_ko_nai_20130116'
rf_paths = glob.glob(os.path.join(basedir, 'fileconversion', 'RF*.h5'))
rr_paths = glob.glob(os.path.join(basedir, 'fileconversion', 'RR*.h5'))
voc_paths = glob.glob(os.path.join(basedir, 'fileconversion', 'VOC*.h5'))
ix2freq = RF.get_ix2freq()
p = re.compile('RF(\d+).h5')
for rf_path in rf_paths:
penno = p.findall(rf_path)[0]
rr_path = [f for f in rr_paths if penno in f][0]
voc_path = [f for f in voc_paths if penno in f][0]
rf_file = h5py.File(rf_path, 'r')
rf_rast = rf_file['rast'].value
rf_stimparams = rf_file['stimID'].value
rf_file.close()
cfs = load_cfs(experiment)
cf_ix = np.int32(np.round(RF.find_cf(cfs, np.int32(penno))))
cf = ix2freq[20:][cf_ix]
# perform analysis
if len(rr_path) > 0:
