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')
rot = pos_angle_inc*(i) + pos_angle_offset
l, m = rotate_lm_CCW(l0, m0, rot)
ra_angle, dec_angle = ra_dec_from_lm(l, m, target.ra_angle, target.dec_angle)
cal_fields.append(momxml.TargetSource('3C196-' + fieldname,
ra_angle = ra_angle,
dec_angle = dec_angle))
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,
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,
start_date=ephem.Date(current_date).tuple(),
duration_seconds=target_duration_s,
stations=stations,
