def generate_snapshot(filename):
conn = GSMConnectionManager(database='stress',
use_monet=False).get_connection()
decl_center = -50
while decl_center < -30:
decl_center = degrees(acos(2 * random.random() - 1) - 0.5 * pi)
ra_center = degrees(2 * pi * random.random())
band = random.random_integers(1, 3)
sql = get_field(ra_center, decl_center, 3.0, 2)
cur = conn.get_cursor(sql)
f = open(filename, 'w')
f.write('# RA DEC Total_flux e_Total_flux\n\n')
for data in iter(cur.fetchone, None):
f.write('%s %s %s %s\n' %
(random.normal(data[0], ERROR), random.normal(data[1], ERROR),
random.normal(data[2], FLUX_ERROR), data[3]))
f.close()
parsetname = path.basename(filename)
parsetname = parsetname[:parsetname.index('.')] + '.parset'
write_parset(parsetname, filename, FREQUENCY[band], ra_center, decl_center)
def generate_snapshot(filename):
conn = GSMConnectionManager(database='stress',
use_monet=False).get_connection()
decl_center = -50
while decl_center < -30:
decl_center = degrees(acos(2 * random.random() - 1) - 0.5 * pi)
ra_center = degrees(2 * pi * random.random())
band = random.random_integers(1, 3)
sql = get_field(ra_center, decl_center, 3.0, 2)
cur = conn.get_cursor(sql)
f = open(filename, 'w')
f.write('# RA DEC Total_flux e_Total_flux\n\n')
for data in iter(cur.fetchone, None):
f.write('%s %s %s %s\n' % (random.normal(data[0], ERROR),
random.normal(data[1], ERROR),
random.normal(data[2], FLUX_ERROR),
data[3]))
f.close()
parsetname = path.basename(filename)
parsetname = parsetname[:parsetname.index('.')] + '.parset'
write_parset(parsetname, filename, FREQUENCY[band], ra_center, decl_center)
