def get_flist(run):
import glob
fs=shapesim.get_default_fs()
f=shapesim.get_output_url(run, 0, 0, itrial=0, fs=fs)
d=os.path.dirname(f)
if 'ngmix' in run:
flist=glob.glob(d+'/*.fits')
else:
flist=glob.glob(d+'/*.rec')
return flist
def check_is2n(run, is2n, flist, full=False):
fs=get_default_fs()
c = read_config(run)
npair,nsplit = get_npair_nsplit(c, is2n)
ntot=npair*2
for isplit in xrange(nsplit):
f=get_output_url(run, 0, is2n, itrial=isplit, fs=fs)
if f not in flist:
print 'missing:',f
continue
if full:
t=None
try:
t=eu.io.read(f)
except:
print 'failed to read:',f
if t is not None:
if t.size != ntot:
print "expected %d, got %d" % (npair,t.size)