def load_dct(cwd=my_dir, pattern='*.dat'): # FindAllFiles print 'cwd:', cwd print 'pattern:', pattern dct_find = {'cwd': cwd, 'pattern': pattern} x = FindAllFiles(dct_find) # print x.pattern x.get_files() x.print_query(x.dct) print len(x.dct.keys()), 'of', x.total # x.dct (last pos.) # set9 = x.sort_dirname(-1,x.dct) # x.print_ [query,class] # x.query_ [dirname,file,filename](searchstring,pos,dct) # x.remove_[dirname,file,filename](searchstring,pos,dct) # x.print_query(set9) # print len(set9.keys()),'of',x.total # sys.exit() # return set9 return x.dct
def check_hsa_analysis(dct): for i,k in enumerate(dct.keys()): os.chdir(dct[k]['dirname']) # FindAllFiles dct_hsa = {'cwd':dct[k]['dirname'],'pattern':'*.dat'} H = FindAllFiles(dct_hsa) set9 = H.get_files() # H.query_ [dirname,file,filename](searchstring,pos,dct) # H.remove_[dirname,file,filename](searchstring,pos,dct) set9 = H.query_file('hsa',set9) set9 = H.sort_dirname(-1,set9) H.print_query(set9) # set9 print(len(set9.keys()),'of',H.total) # if ((not os.path.exists('hsa_current')) and (len(set9.keys()) != 0)): # if (not os.path.exists('hsa_%s' % new_dir)): if ((dct[traj]['dirname'] == dct[k]['dirname']) and (not os.path.exists('hsa_current'))): print('making hsa_current') os.makedirs('hsa_current') # hcount = 0 for h in set9.keys(): # print h if os.path.basename(set9[h]['dirname']) == os.path.basename(dct[k]['dirname']): new_hsa_cur = os.path.join(dct[k]['dirname'],'hsa_current',set9[h]['filename']) # print 'moving to hsa_current:',new_hsa_cur # hcount += 1 os.rename(set9[h]['file'],new_hsa_cur) lst_hsa = glob.glob('hsa_current/hsa*.dat') # 1 if len(lst_hsa) > 0: print('base:',os.path.basename(dct[k]['dirname'])) print(lst_hsa)
set9 = x.remove_dirname('fail', -2, set9) set9 = x.remove_filename('LOG', set9) set9 = x.remove_filename('contact', set9) set9 = x.remove_filename('mt_angles', set9) set9 = x.remove_filename('indices', set9) set9 = x.remove_filename('tubulin_inter', set9) # x.print_query(set9) set9 = x.query_dirname('3kinesin13', -3, set9) # set9 = x.sort_dirname(-2,set9) # x.print_query(set9) # sys.exit() # projects: 123, 134, 145 set9 = x.query_dirname(dataname, -2, set9) set9 = x.sort_dirname(-2, set9) x.print_query(set9) point_start = 0 # (3000001, 7) -- confirmed point_end = 5000000 ma_value = 2500 step = 10 ts = 20.00 nav = 1000 dcdfreq = 100000 d = PlotSop('gsop', point_start, point_end, ma_value, step, ts, nav, dcdfreq) for i in range(len(set9.keys())): path = set9[i]['file'] print i, path label = dataname + '-' + str(i)