# 5-sigma depth percentile # standard target map (WFD, NES, SCP, GP) # Slewtime # mask lots of off-meridian space # No pairs # Greedy selection of opsim fields if __name__ == "__main__": survey_length = 365.25 * 10. # days # Define what we want the final visit ratio map to look like target_map = fs.standard_goals()['r'] filtername = 'r' bfs = [] bfs.append(fs.Depth_percentile_basis_function(filtername=filtername)) bfs.append( fs.Target_map_basis_function(target_map=target_map, filtername=filtername, out_of_bounds_val=hp.UNSEEN)) bfs.append(fs.North_south_patch_basis_function(zenith_min_alt=50.)) bfs.append(fs.Slewtime_basis_function(filtername=filtername)) weights = np.array([1., 0.2, 1., 2.]) survey = fs.Greedy_survey_fields(bfs, weights, block_size=1, filtername=filtername) scheduler = fs.Core_scheduler([survey]) observatory = Speed_observatory()
import numpy as np import lsst.sims.featureScheduler as fs from lsst.sims.featureScheduler.observatory import Speed_observatory if __name__ == "__main__": survey_length = 2. # days # Define what we want the final visit ratio map to look like target_map = fs.standard_goals()['r'] bfs = [] bfs.append(fs.Depth_percentile_basis_function()) bfs.append(fs.Target_map_basis_function(target_map=target_map)) weights = np.array([5, 1]) survey = fs.Simple_greedy_survey(bfs, weights) scheduler = fs.Core_scheduler([survey]) observatory = Speed_observatory() observatory, scheduler, observations = fs.sim_runner( observatory, scheduler, survey_length=survey_length, filename='one_survey.db')