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 ])
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, grouping_parent = True, mom_id = 449283)
dec_angle=dec_3c61_1)) sys.stderr.write('PRE : ' + str(pre_cal) + '\n') sys.stderr.write('MAIN: ' + str(target) + '\n') 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)) current_date += cal_duration_s * ephem.second + 61 * ephem.second observations.append( momxml.Observation( beam_list=[momxml.Beam(target, target_subbands)] + [momxml.Beam(field, cal_subbands) for field in cal_fields],
sys.stderr.write('MAIN: '+str(target) + '\n') for field in aux_pointings: sys.stderr.write(' GRID: '+str(field)+ '\n') all_fields = [target] + aux_pointings 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(field, subbands) for field in all_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, mom_id = 322103) print momxml.xml([obs_folder], 'LC1_039')
sys.stderr.write(' CAL: '+str(cal)+ '\n') observations = [] 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') backend = momxml.BackendProcessing(channels_per_subband = chan, integration_time_seconds = int_s) 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)) 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)] +
target_duration_s = 8 * 3600 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, aux_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, clock_mhz=200, backend=backend, bit_mode=8, initial_status='approved')) obs_folder = momxml.Folder(name='EoR-3C196', description='EoR 3C 196 observations', children=observations, grouping_parent=True)
#cal_fields.append(scintillation_3c220_3) sys.stderr.write('PRE : '+str(pre_cal)+ '\n') sys.stderr.write('MAIN: '+str(target) + '\n') 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,