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,