def generate_image(filename, sourcename, band, size):
f = open(sourcename, 'r')
fo = open(filename, 'w')
fo.write(f.readline())
#fo.write(f.readline())
for _ in range(size):
z = f.readline().split()
for ik in range(2):
z[ik] = random.normal(z[ik], ERROR)
fo.write('%s %s %s %s\n' % (z[0], z[1], z[2], 0.01))
parsetname = path.basename(filename)
parsetname = parsetname[:parsetname.index('.')] + '.parset'
write_parset(parsetname, filename, FREQUENCY[band], 180.0, 0.0, 5.0)
f.close()
fo.close()
def generate_image(filename, sourcename, band, size):
f = open(sourcename, 'r')
fo = open(filename, 'w')
fo.write(f.readline())
#fo.write(f.readline())
for _ in xrange(size):
z = f.readline().split()
for ik in xrange(2):
z[ik] = random.normal(z[ik], ERROR)
fo.write('%s %s %s %s\n' % (z[0], z[1], z[2], 0.01))
parsetname = path.basename(filename)
parsetname = parsetname[:parsetname.index('.')] + '.parset'
write_parset(parsetname, filename, FREQUENCY[band], 180.0, 0.0, 5.0)
f.close()
fo.close()
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)
def generate_field_parset(filename, ra, decl, radius, size,
frequency=8):
generate_field_file(filename, ra, decl, radius, size)
parsetname = path.basename(filename)
parsetname = parsetname[:parsetname.index('.')] + '.parset'
write_parset(parsetname, filename, FREQUENCY[frequency], ra, decl, radius)
def generate_field_parset(filename, ra, decl, radius, size,
frequency=8):
generate_field_file(filename, ra, decl, radius, size)
parsetname = path.basename(filename)
parsetname = parsetname[:parsetname.index('.')] + '.parset'
write_parset(parsetname, filename, FREQUENCY[frequency], ra, decl, radius)
