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: