def single_observation(start_search_0ut):
start_search_date = momxml.ephem.Date(start_search_0ut + ' 14:00:00')
target = momxml.TargetSource(name='3C196',
ra_angle=momxml.Angle(shms=('+', 8, 13,
36.0678)),
dec_angle=momxml.Angle(sdms=('+', 48, 13,
2.581)))
auxiliary = momxml.TargetSource(name='3C196-32A',
ra_angle=momxml.Angle(shms=('+', 8, 22,
36.91)),
dec_angle=momxml.Angle(sdms=('+', 51, 12,
29.5)))
transit_date = momxml.next_date_with_lofar_lst(target.ra_angle.as_rad(),
start_search_date)
target.name = '3C196-%4d-%02d-%02d' % start_search_date.tuple()[0:3]
antenna_set = 'HBA_DUAL_INNER'
band = 'HBA_LOW'
stations = momxml.station_list('all')
int_s = 2.0 # 2.0
chan = 64
target_subbands = '62..393' # 380 subbands 115 -- 189 MHz
subbands_main = '70..394'
subbands_aux = '70..232'
target_duration_s = 6 * 3600
sys.stderr.write('MAIN: ' + str(target) + '\n')
sys.stderr.write(' AUX: ' + str(auxiliary) + '\n')
current_date = transit_date - 0.5 * target_duration_s * ephem.second
backend = momxml.BackendProcessing(channels_per_subband=chan,
integration_time_seconds=int_s)
return momxml.Observation(beam_list=[
momxml.Beam(target, subbands_main, storage_cluster='CEP4'),
momxml.Beam(auxiliary, subbands_aux, storage_cluster='CEP4')
],
antenna_set=antenna_set,
frequency_range=band,
start_date=ephem.Date(current_date).tuple(),
duration_seconds=target_duration_s,
stations=stations,
clock_mhz=200,
backend=backend,
bit_mode=8,
initial_status='approved')
def make_obs(self, name, targets, start_date, duration_s):
#stderr(str(targets))
return momxml.Observation(name=name,
antenna_set='HBA_DUAL_INNER',
frequency_range='HBA_LOW',
start_date=ephem.Date(start_date).tuple(),
duration_seconds=duration_s,
stations=self.stations,
clock_mhz=200,
integration_time_seconds=self.int_s,
bit_mode=self.bit_mode,
channels_per_subband=self.channels,
beam_list=[
momxml.Beam(target, self.subbands)
for target in targets
])
sys.stderr.write('MAIN: '+str(target) + '\n')
for cal in cal_fields:
sys.stderr.write(' CAL: '+str(cal)+ '\n')
observations = []
current_date = transit_date - 0.5*target_duration_s*ephem.second
backend = momxml.BackendProcessing(channels_per_subband = chan,
integration_time_seconds = int_s)
observations.append(momxml.Observation(
beam_list = [momxml.Beam(target, target_subbands)] + [momxml.Beam(field, cal_subbands) for field in cal_fields],
antenna_set = antenna_set,
frequency_range = band,
start_date = ephem.Date(current_date).tuple(),
duration_seconds = target_duration_s,
stations = stations,
clock_mhz = 200,
backend = backend,
bit_mode = 8))
obs_folder = momxml.Folder(name = 'EoR-3C196',
description = 'EoR 3C 196 observations',
children = observations,
grouping_parent = True,
mom_id = 449283)
with open('eor-3c196-%d%02d%02d.xml' % (start_search_date.tuple()[0:3]), 'w') as output:
output.write(momxml.xml([obs_folder], 'LC3_028'))
for cal in cal_fields:
sys.stderr.write(' CAL: '+str(cal)+ '\n')
sys.stderr.write('POST: '+ str(post_cal)+'\n')
backend = momxml.BackendProcessing(channels_per_subband = chan,
integration_time_seconds = int_s)
observations = []
current_date = start_date
observations.append(momxml.Observation(
beam_list = [momxml.Beam(pre_cal, target_subbands)],
antenna_set = antenna_set,
frequency_range = band,
start_date = ephem.Date(current_date).tuple(),
duration_seconds = cal_duration_s,
stations = stations,
clock_mhz = 200,
bit_mode = 8,
backend = backend,
initial_status = 'approved'))
current_date += cal_duration_s*ephem.second + 61*ephem.second
observations.append(momxml.Observation(
beam_list = [momxml.Beam(target, target_subbands)] + [momxml.Beam(field, subbands) for field, subbands in zip(cal_fields, cal_subbands)],
antenna_set = antenna_set,
frequency_range = band,
start_date = ephem.Date(current_date).tuple(),
duration_seconds = target_duration_s,
stations = stations,
current_date = transit_date - (0.5 * target_duration_s + cal_duration_s +
61) * ephem.second
source_catalogue = momxml.SourceCatalogue()
pre_cal = source_catalogue.cal_source(current_date, 'HBA')
post_cal = source_catalogue.cal_source(
current_date +
(target_duration_s + 2 * cal_duration_s + 120) * ephem.second, 'HBA')
observations.append(
momxml.Observation(beam_list=[momxml.Beam(pre_cal, target_subbands)],
antenna_set=antenna_set,
frequency_range=band,
start_date=ephem.Date(current_date).tuple(),
duration_seconds=cal_duration_s,
stations=stations,
clock_mhz=200,
backend=momxml.BackendProcessing(
integration_time_seconds=int_s,
channels_per_subband=chan),
bit_mode=8))
current_date += cal_duration_s * ephem.second + 61 * ephem.second
observations.append(
momxml.Observation(
beam_list=[momxml.Beam(target, target_subbands)] +
[momxml.Beam(cal_fields[0], gb87_subbands)] +
[momxml.Beam(field, cal_subbands) for field in cal_fields[1:]],
antenna_set=antenna_set,
frequency_range=band,
