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,