def test_get_slew_delay(self):
self.model.update_state(0)
# Use old values, to avoid updating final states.
self.model.params.rotator_followsky = True
self.assertEqual(
str(self.model.current_state),
"t=0.0 ra=29.480 dec=-26.744 ang=180.000 "
"filter=r track=False alt=86.500 az=0.000 pa=180.000 rot=0.000 "
"telaz=0.000 telrot=0.000 "
"mounted=['g', 'r', 'i', 'z', 'y'] unmounted=['u']",
)
# This slew will include a CL optics correction.
target = Target()
target.ra_rad = math.radians(60)
target.dec_rad = math.radians(-20)
target.ang_rad = math.radians(0)
target.filter = "r"
delay, status = self.model.get_slew_delay(target)
self.assertAlmostEqual(delay, 85.507, delta=1e-3)
self.model.slew(target)
# This slew simply includes a filter change.
target = Target()
target.ra_rad = math.radians(60)
target.dec_rad = math.radians(-20)
target.ang_rad = math.radians(0)
target.filter = "g"
delay, status = self.model.get_slew_delay(target)
self.assertAlmostEqual(delay, 120, delta=1e-3)
# This slew does not include OL correction, but does involve dome
# crawl.
target = Target()
target.ra_rad = math.radians(50)
target.dec_rad = math.radians(-10)
target.ang_rad = math.radians(10)
target.filter = "r"
delay, status = self.model.get_slew_delay(target)
self.assertAlmostEqual(delay, 17.913, delta=1e-3)
# This slew is only readout.
self.model.slew(target)
delay, status = self.model.get_slew_delay(target)
self.assertAlmostEqual(delay, 2.0, delta=1e-3)
# This slew involves rotator.
target.ang_rad = math.radians(15)
delay, status = self.model.get_slew_delay(target)
self.assertAlmostEqual(delay, 4.472, delta=1e-3)
def test_get_slew_delay_followsky_false(self):
# Test slew time without followsky option. Similar to
# test_get_slew_delay above.
self.model.update_state(0)
self.model.params.rotator_followsky = False
expected_state = (
"t=0.0 ra=29.480 dec=-26.744 ang=180.000 "
"filter=r track=False alt=86.500 az=0.000 pa=180.000 rot=0.000 "
"telaz=0.000 telrot=0.000 "
"mounted=['g', 'r', 'i', 'z', 'y'] unmounted=['u']")
self.assertEqual(str(self.model.current_state), expected_state)
target = Target()
target.ra_rad = math.radians(60)
target.dec_rad = math.radians(-20)
target.ang_rad = math.radians(0)
target.filter = "r"
delay, status = self.model.get_slew_delay(target)
self.assertAlmostEqual(delay, 85.507, delta=1e-3)
self.model.slew(target)
target = Target()
target.ra_rad = math.radians(60)
target.dec_rad = math.radians(-20)
target.ang_rad = math.radians(0)
target.filter = "g"
delay, status = self.model.get_slew_delay(target)
self.assertAlmostEqual(delay, 120, delta=1e-3)
target = Target()
target.ra_rad = math.radians(50)
target.dec_rad = math.radians(-10)
target.ang_rad = math.radians(10)
target.filter = "r"
delay, status = self.model.get_slew_delay(target)
self.assertAlmostEqual(delay, 17.913, delta=1e-3)
self.model.slew(target)
delay, status = self.model.get_slew_delay(target)
self.assertAlmostEqual(delay, 2.0, delta=1e-3)
# Here is the difference when using followsky = False
target.ang_rad = math.radians(15)
delay, status = self.model.get_slew_delay(target)
self.assertAlmostEqual(delay, 2.0, delta=1e-3)
def read_script(self):
scriptfilepath = self.script_file
lines = file(scriptfilepath).readlines()
targetid = 0
self.targetsList = []
for line in lines:
line = line.strip()
if not line: # skip blank line
continue
if line[0] == '#': # skip comment line
continue
targetid += 1
values = line.split()
target = Target()
target.fieldid = eval(values[0])
target.filter = values[1]
target.ra_rad = math.radians(eval(values[2]))
target.dec_rad = math.radians(eval(values[3]))
target.ang_rad = math.radians(eval(values[4]))
target.num_exp = eval(values[5])
target.exp_times = [int(x) for x in values[6].split(',')]
self.targetsList.append(target)
self.log.info("%d targets" % len(self.targetsList))
def test_park(self):
self.model.update_state(0)
self.model.params.rotator_followsky = False
self.model.params.rotator_resume_angle = False
# Start at park, slew to target.
# Use default configuration (dome crawl, CL updates, etc.)
target = Target()
target.ra_rad = math.radians(60)
target.dec_rad = math.radians(-20)
target.ang_rad = math.radians(0)
target.filter = "z"
self.model.slew(target)
expected_state = "t=156.0 ra=60.000 dec=-20.000 ang=243.495 filter=z track=True " \
"alt=61.191 az=76.196 pa=243.224 rot=359.729 telaz=76.196 telrot=-0.271 " \
"mounted=['g', 'r', 'i', 'z', 'y'] unmounted=['u']"
self.assertEqual(str(self.model.current_state), expected_state)
self.check_delay_and_state(
self.model,
self.make_slewact_dict((8.387, 11.966, 0.0, 7.387, 36.0, 18.775,
48.507, 1.0, 120.0, 2.0)),
['telopticsclosedloop', 'filter'], (-3.50, 7.00, 0.0, -1.75, 1.50))
self.model.park()
expected_state = "t=241.1 ra=30.487 dec=-26.744 ang=180.000 filter=z track=False " \
"alt=86.500 az=0.000 pa=180.000 rot=0.000 telaz=0.000 telrot=0.000 " \
"mounted=['g', 'r', 'i', 'z', 'y'] unmounted=['u']"
self.assertEqual(str(self.model.current_state), expected_state)
self.check_delay_and_state(
self.model,
self.make_slewact_dict(
(8.231, 11.885, 1.041, 7.231, 36.0, 18.462, 48.130, 1.0, 0.0,
2.0)), ['telopticsclosedloop', 'domazsettle', 'domaz'],
(3.50, -7.00, 0.520, 1.75, -1.50))
def init_target_with_bad_json(self):
missingfilter = (
'{"targetid": 3, "fieldid": 2573, "ra_rad": 5.24504477561707, '
'"dec_rad": -0.030036505561584215, "ang_rad": 0.7853981633974483, '
'"num_exp": 2, "exp_times": [15.0, 15.0], "_exp_time": null, '
'"time": 0.0, "airmass": 0.0, "sky_brightness": 0.0, "cloud": 0.0, '
'"seeing": 0.0, "propid": 0, "need": 0.0, "bonus": 0.0, "value": 0.0, '
'"goal": 0, "visits": 0, "progress": 0.0, '
'"sequenceid": 0, "subsequencename": "", "groupid": 0, "groupix": 0, '
'"is_deep_drilling": false, "is_dd_firstvisit": false, "remaining_dd_visits": 0, '
'"dd_exposures": 0, "dd_filterchanges": 0, "dd_exptime": 0.0, '
'"alt_rad": 0.7853981633974483, "az_rad": 3.9269908169872414, '
'"rot_rad": 0.5235987755982988, "telalt_rad": 0.7853981633974483, '
'"telaz_rad": 3.9269908169872414, "telrot_rad": 0.5235987755982988, "propboost": 1.0, '
'"slewtime": 0.0, "cost": 0.0, "rank": 0.0, "num_props": 0, "propid_list": [], '
'"need_list": [], "bonus_list": [], "value_list": [], "propboost_list": [], '
'"sequenceid_list": [], "subsequencename_list": [], "groupid_list": [], '
'"groupix_list": [], "is_deep_drilling_list": [], '
'"is_dd_firstvisit_list": [], "remaining_dd_visits_list": [], '
'"dd_exposures_list": [], "dd_filterchanges_list": [], "dd_exptime_list": [], '
'"last_visit_time": 0.0, "note": ""}'
)
t = Target()
t.from_json(missingfilter)
def test_driver_state_copy(self):
target2 = Target()
target2.copy_driver_state(self.target)
self.assertEqual(target2.alt_rad, self.alt_rad)
self.assertEqual(target2.az_rad, self.az_rad)
self.assertEqual(target2.ang_rad, self.ang_rad)
self.assertEqual(target2.rot_rad, self.rot_rad)
self.assertEqual(target2.telalt_rad, self.telalt_rad)
self.assertEqual(target2.telaz_rad, self.telaz_rad)
self.assertEqual(target2.telrot_rad, self.telrot_rad)
def test_slew(self):
self.model.update_state(0)
self.model.params.domaz_free_range = 0
self.model.params.optics_cl_delay = [0, 20.0]
self.model.params.rotator_followsky = True
self.assertEqual(
str(self.model.current_state),
"t=0.0 ra=29.480 dec=-26.744 ang=180.000 "
"filter=r track=False alt=86.500 az=0.000 pa=180.000 rot=0.000 "
"telaz=0.000 telrot=0.000 "
"mounted=['g', 'r', 'i', 'z', 'y'] unmounted=['u']",
)
target = Target()
target.ra_rad = math.radians(60)
target.dec_rad = math.radians(-20)
target.ang_rad = math.radians(0)
target.filter = "r"
self.model.slew(target)
self.assertEqual(
str(self.model.current_state),
"t=74.2 ra=60.000 dec=-20.000 ang=180.000 "
"filter=r track=True alt=60.904 az=76.495 pa=243.368 rot=63.368 "
"telaz=76.495 telrot=63.368 "
"mounted=['g', 'r', 'i', 'z', 'y'] unmounted=['u']",
)
target = Target()
target.ra_rad = math.radians(60)
target.dec_rad = math.radians(-20)
target.ang_rad = math.radians(0)
target.filter = "i"
self.model.slew(target)
self.assertEqual(
str(self.model.current_state),
"t=194.2 ra=60.000 dec=-20.000 ang=180.000 "
"filter=i track=True alt=61.324 az=76.056 pa=243.156 rot=63.156 "
"telaz=76.056 telrot=63.156 "
"mounted=['g', 'r', 'i', 'z', 'y'] unmounted=['u']",
)
def test_get_deep_drilling_time(self):
target = Target()
target.is_deep_drilling = True
target.is_dd_firstvisit = True
target.remaining_dd_visits = 96
target.dd_exposures = 2 * 96
target.dd_filterchanges = 3
target.dd_exptime = 96 * 2 * 15.0
ddtime = self.model.get_deep_drilling_time(target)
self.assertEqual(ddtime, 3808.0)
def test_json_serialization(self):
jsondump = self.target.to_json()
target2 = Target()
target2.from_json(jsondump)
self.assertEqual(self.target.targetid, target2.targetid)
self.assertEqual(self.target.fieldid, target2.fieldid)
self.assertEqual(self.target.filter, target2.filter)
self.assertEqual(self.target.ra_rad, target2.ra_rad)
self.assertEqual(self.target.dec_rad, target2.dec_rad)
self.assertEqual(self.target.num_exp, target2.num_exp)
self.assertListEqual(self.target.exp_times, target2.exp_times)
def __init__(self):
self.log = logging.getLogger("schedulerDriver")
self.params = DriverParameters()
self.location = ObservatoryLocation()
self.observatoryModel = ObservatoryModel(self.location, WORDY)
self.observatoryModel2 = ObservatoryModel(self.location, WORDY)
self.observatoryState = ObservatoryState()
self.sky = AstronomicalSkyModel(self.location)
self.db = FieldsDatabase()
self.build_fields_dict()
self.propid_counter = 0
self.science_proposal_list = []
self.start_time = 0.0
self.time = 0.0
self.targetid = 0
self.survey_started = False
self.isnight = False
self.sunset_timestamp = 0.0
self.sunrise_timestamp = 0.0
self.survey_duration_DAYS = 0.0
self.survey_duration_SECS = self.survey_duration_DAYS * 24 * 60 * 60.0
self.darktime = False
self.mounted_filter = ""
self.unmounted_filter = ""
self.midnight_moonphase = 0.0
self.nulltarget = Target()
self.nulltarget.targetid = -1
self.nulltarget.num_exp = 1
self.nulltarget.exp_times = [0.0]
self.nulltarget.num_props = 1
self.nulltarget.propid_list = [0]
self.nulltarget.need_list = [0.0]
self.nulltarget.bonus_list = [0.0]
self.nulltarget.value_list = [0.0]
self.nulltarget.propboost_list = [1.0]
self.last_winner_target = self.nulltarget.get_copy()
self.deep_drilling_target = None
self.need_filter_swap = False
self.filter_to_unmount = ""
self.filter_to_mount = ""
self.cloud = 0.0
self.seeing = 0.0
def generate_target(self, fb_observation, tag='generate'):
'''Takes an observation array given by the feature based scheduler and generate an appropriate OpSim target.
:param fb_observation: numpy.array
:return: Target
'''
# self.log.debug('%s: %s' % (tag, fb_observation))
self.targetid += 1
filtername = fb_observation['filter']
propid = fb_observation['survey_id']
target = Target()
target.targetid = self.targetid
target.fieldid = fb_observation['field_id']
target.filter = str(filtername)
target.num_exp = fb_observation['nexp']
target.exp_times = [
fb_observation['exptime'] / fb_observation['nexp']
] * fb_observation['nexp']
target.ra_rad = fb_observation['RA']
target.dec_rad = fb_observation['dec']
target.ang_rad = fb_observation['rotSkyPos']
target.propid = propid
target.goal = 100
target.visits = 0
target.progress = 0.0
target.groupid = 1
target.groupix = 0
target.propid_list = [propid]
target.need_list = [target.need]
target.bonus_list = [target.bonus]
target.value_list = [target.value]
target.propboost_list = [target.propboost]
target.sequenceid_list = [target.sequenceid]
target.subsequencename_list = [target.subsequencename]
target.groupid_list = [target.groupid]
target.groupix_list = [target.groupix]
target.is_deep_drilling_list = [target.is_deep_drilling]
target.is_dd_firstvisit_list = [target.is_dd_firstvisit]
target.remaining_dd_visits_list = [target.remaining_dd_visits]
target.dd_exposures_list = [target.dd_exposures]
target.dd_filterchanges_list = [target.dd_filterchanges]
target.dd_exptime_list = [target.dd_exptime]
return target
def setUp(self):
self.targetId = 3
self.fieldId = 2573
self.band_filter = "r"
self.ra = 300.518929
self.dec = -1.720965
self.ang = 45.0
self.num_exposures = 2
self.exposure_times = [15.0, 15.0]
self.alt = 45.0
self.az = 225.0
self.rot = 30.0
self.telalt = 45.0
self.telaz = 225.0
self.telrot = 30.0
self.ra_rad = math.radians(self.ra)
self.dec_rad = math.radians(self.dec)
self.ang_rad = math.radians(self.ang)
self.alt_rad = math.radians(self.alt)
self.az_rad = math.radians(self.az)
self.rot_rad = math.radians(self.rot)
self.telalt_rad = math.radians(self.telalt)
self.telaz_rad = math.radians(self.telaz)
self.telrot_rad = math.radians(self.telrot)
self.target = Target(
self.targetId,
self.fieldId,
self.band_filter,
self.ra_rad,
self.dec_rad,
self.ang_rad,
self.num_exposures,
self.exposure_times,
)
self.target.alt_rad = self.alt_rad
self.target.az_rad = self.az_rad
self.target.rot_rad = self.rot_rad
self.target.telalt_rad = self.telalt_rad
self.target.telaz_rad = self.telaz_rad
self.target.telrot_rad = self.telrot_rad
def test_domecrawl(self):
self.model.update_state(0)
self.assertEqual(
str(self.model.current_state),
"t=0.0 ra=29.480 dec=-26.744 ang=180.000 "
"filter=r track=False alt=86.500 az=0.000 pa=180.000 rot=0.000 "
"telaz=0.000 telrot=0.000 "
"mounted=['g', 'r', 'i', 'z', 'y'] unmounted=['u']")
target = Target()
target.ra_rad = math.radians(35)
target.dec_rad = math.radians(-27)
target.ang_rad = math.radians(0)
target.filter = "r"
# Just test whether dome crawl is faster or not.
# If we test the final slew state, this is including other aspects of slew model (such as CLoptics).
self.model.params.domaz_free_range = 0
delay_nocrawl = self.model.get_slew_delay(target)
self.model.params.domaz_free_range = np.radians(4.0)
delay_crawl = self.model.get_slew_delay(target)
self.assertTrue(delay_crawl < delay_nocrawl)
def test_slewdata(self):
self.model.update_state(0)
# Use old values, to avoid updating final states.
self.model.params.domaz_free_range = 0
self.model.params.optics_cl_delay = [0, 20.0]
self.model.params.rotator_followsky = True
target = Target()
target.ra_rad = math.radians(60)
target.dec_rad = math.radians(-20)
target.ang_rad = math.radians(0)
target.filter = "r"
self.model.slew(target)
self.assertEqual(
str(self.model.current_state),
"t=74.2 ra=60.000 dec=-20.000 ang=180.000 "
"filter=r track=True alt=60.904 az=76.495 pa=243.368 rot=63.368 "
"telaz=76.495 telrot=63.368 "
"mounted=['g', 'r', 'i', 'z', 'y'] unmounted=['u']")
self.check_delay_and_state(
self.model,
self.make_slewact_dict(
(8.387, 11.966, 21.641, 7.387, 20.0, 18.775, 53.174, 1.0, 0.0,
2.0)), ['telopticsclosedloop', 'domazsettle', 'domaz'],
(-3.50, 7.00, 3.50, -1.75, 1.50))
target = Target()
target.ra_rad = math.radians(60)
target.dec_rad = math.radians(-20)
target.ang_rad = math.radians(0)
target.filter = "i"
self.model.slew(target)
self.assertEqual(
str(self.model.current_state),
"t=194.2 ra=60.000 dec=-20.000 ang=180.000 "
"filter=i track=True alt=61.324 az=76.056 pa=243.156 rot=63.156 "
"telaz=76.056 telrot=63.156 "
"mounted=['g', 'r', 'i', 'z', 'y'] unmounted=['u']")
self.check_delay_and_state(
self.model,
self.make_slewact_dict(
(0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 120.0, 2.0)),
['filter'], (0, 0, 0, 0, 0))
target = Target()
target.ra_rad = math.radians(61)
target.dec_rad = math.radians(-21)
target.ang_rad = math.radians(1)
target.filter = "i"
self.model.slew(target)
self.assertEqual(
str(self.model.current_state),
"t=199.0 ra=61.000 dec=-21.000 ang=181.000 "
"filter=i track=True alt=60.931 az=78.751 pa=245.172 rot=64.172 "
"telaz=78.751 telrot=64.172 "
"mounted=['g', 'r', 'i', 'z', 'y'] unmounted=['u']")
self.check_delay_and_state(
self.model,
self.make_slewact_dict((0.683, 1.244, 2.022, 0.117, 0.0, 1.365,
3.801, 1.0, 0.000, 2.000)),
['domazsettle', 'domaz'], (-1.194, 4.354, 1.011, -0.598, 1.425))
def start_night(self, timestamp, filters_mounted_tonight_list, night):
Proposal.start_night(self, timestamp, filters_mounted_tonight_list,
night)
self.fieldsvisitedtonight = {}
self.tonight_filters_list = []
self.tonight_targets = 0
self.tonight_targets_dict = {}
for filter in filters_mounted_tonight_list:
self.tonight_filters_list.append(filter)
self.build_tonight_fields_list(timestamp, night)
self.tonight_fields = len(self.tonight_fields_list)
# compute at start night
for field in self.tonight_fields_list:
fieldid = field.fieldid
# add new fields to fields dictionary
if fieldid not in self.survey_fields_dict:
self.survey_fields_dict[fieldid] = field.get_copy()
self.survey_fields += 1
# add new fields to targets dictionary
if fieldid not in self.survey_targets_dict:
self.survey_targets_dict[fieldid] = {}
self.tonight_targets_dict[fieldid] = {}
for filter in filters_mounted_tonight_list:
if filter in self.params.filter_list:
if self.params.filter_goal_dict[filter] > 0:
if filter not in self.survey_targets_dict[fieldid]:
target = Target()
target.fieldid = fieldid
target.filter = filter
target.num_exp = self.params.filter_num_exp_dict[
filter]
target.exp_times = self.params.filter_exp_times_dict[
filter]
target.ra_rad = field.ra_rad
target.dec_rad = field.dec_rad
target.propid = self.propid
target.goal = self.params.filter_goal_dict[filter]
target.visits = 0
target.progress = 0.0
target.groupid = 1
target.groupix = 1
self.survey_targets_dict[fieldid][filter] = target
self.survey_targets += 1
self.survey_targets_goal += target.goal
self.survey_filters_goal_dict[
filter] += target.goal
target = self.survey_targets_dict[fieldid][filter]
if target.progress < 1.0:
self.tonight_targets_dict[fieldid][filter] = target
self.tonight_targets += 1
for filter in self.params.filter_list:
if self.survey_filters_goal_dict[filter] > 0:
self.survey_filters_progress_dict[filter] = \
float(self.survey_filters_visits_dict[filter]) / self.survey_filters_goal_dict[filter]
else:
self.survey_filters_progress_dict[filter] = 0.0
if self.survey_targets_goal > 0:
self.survey_targets_progress = float(
self.survey_targets_visits) / self.survey_targets_goal
else:
self.survey_targets_progress = 0.0
self.log.debug("start_night tonight fields=%i targets=%i" %
(self.tonight_fields, self.tonight_targets))
self.last_observation = None
self.last_observation_was_for_this_proposal = False
def observe(self, observation):
"""Try to make an observation
Returns
-------
observation : observation object
None if there was no observation taken. Completed observation with meta data filled in.
new_night : bool
Have we started a new night.
"""
start_night = self.night.copy()
# Make sure the kinematic model is set to the correct mjd
t = Time(self.mjd, format='mjd')
self.observatory.update_state(t.unix)
if np.isnan(observation['rotSkyPos']):
observation = self._update_rotSkyPos(observation)
target = Target(band_filter=observation['filter'],
ra_rad=observation['RA'],
dec_rad=observation['dec'],
ang_rad=observation['rotSkyPos'],
num_exp=observation['nexp'],
exp_times=[observation['exptime']])
start_ra = self.observatory.current_state.ra_rad
start_dec = self.observatory.current_state.dec_rad
slewtime, visittime = self.observatory.observe_times(target)
# Check if the mjd after slewtime and visitime is fine:
observation_worked, new_mjd = self.check_mjd(self.mjd +
(slewtime + visittime) /
24. / 3600.)
if observation_worked:
observation['visittime'] = visittime
observation['slewtime'] = slewtime
observation['slewdist'] = _angularSeparation(
start_ra, start_dec, self.observatory.current_state.ra_rad,
self.observatory.current_state.dec_rad)
self.mjd = self.mjd + slewtime / 24. / 3600.
# Reach into the observatory model to pull out the relevant data it has calculated
# Note, this might be after the observation has been completed.
observation['alt'] = self.observatory.current_state.alt_rad
observation['az'] = self.observatory.current_state.az_rad
observation['pa'] = self.observatory.current_state.pa_rad
observation['rotTelPos'] = self.observatory.current_state.rot_rad
observation['rotSkyPos'] = self.observatory.current_state.ang_rad
# Metadata on observation is after slew and settle, so at start of exposure.
result = self.observation_add_data(observation)
self.mjd = self.mjd + visittime / 24. / 3600.
new_night = False
else:
result = None
self.observatory.park()
# Skip to next legitimate mjd
self.mjd = new_mjd
now_night = self.night
if now_night == start_night:
new_night = False
else:
new_night = True
return result, new_night
def __init__(self, propid, field, name, params):
self.propid = propid
self.field = field
self.name = name
self.num_events = params.subsequence_num_events[name]
self.num_max_missed = params.subsequence_num_max_missed[name]
self.time_interval = params.subsequence_time_interval[name]
self.time_window_start = params.subsequence_time_window_start[name]
self.time_window_max = params.subsequence_time_window_max[name]
self.time_window_end = params.subsequence_time_window_end[name]
self.filter_list = params.subsequence_filter_list[name]
self.visits_list = params.subsequence_visits_list[name]
self.filter_num_exp_dict = dict(params.filter_num_exp_dict)
self.filter_exp_times_dict = dict(params.filter_exp_times_dict)
self.num_visits_per_event = 0
self.filters_goal_dict = {}
self.dd_exposures = 0
self.dd_exptime = 0.0
self.num_subevents = len(self.filter_list)
for ix, filter in enumerate(self.filter_list):
if filter not in self.filters_goal_dict:
self.filters_goal_dict[filter] = 0
self.filters_goal_dict[filter] += self.num_events * self.visits_list[ix]
self.num_visits_per_event += self.visits_list[ix]
self.dd_exposures += self.visits_list[ix] * self.filter_num_exp_dict[filter]
self.dd_exptime += self.visits_list[ix] * sum(self.filter_exp_times_dict[filter])
self.goal = self.num_events * self.num_visits_per_event
if (len(self.visits_list) > 1) or (self.visits_list[0] > 1):
self.is_deep_drilling = True
self.dd_visits = self.num_visits_per_event
self.dd_filterchanges = max(len(self.filter_list) - 1, 0)
else:
self.is_deep_drilling = False
self.dd_visits = 0
self.dd_filterchanges = 0
self.subevent_index = 0
self.subevent_num_visits_left = self.visits_list[0]
self.event_visits_list = []
self.all_events_list = []
self.obs_events_list = []
self.mis_events_list = []
self.in_deep_drilling = False
self.target = Target()
self.target.fieldid = field.fieldid
self.target.filter = self.filter_list[self.subevent_index]
self.target.num_exp = self.filter_num_exp_dict[self.target.filter]
self.target.exp_times = self.filter_exp_times_dict[self.target.filter]
self.target.ra_rad = field.ra_rad
self.target.dec_rad = field.dec_rad
self.target.propid = self.propid
self.target.goal = self.goal
self.target.visits = 0
self.target.progress = 0.0
self.target.sequenceid = self.field.fieldid
self.target.subsequencename = self.name
self.target.groupid = 1
self.target.groupix = 1
self.target.is_deep_drilling = self.is_deep_drilling
self.target.is_dd_firstvisit = self.is_deep_drilling
self.target.remaining_dd_visits = self.dd_visits
self.target.dd_exposures = self.dd_exposures
self.target.dd_filterchanges = self.dd_filterchanges
self.target.dd_exptime = self.dd_exptime
self.update_state()