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)
