def run_dark_current_task(sensor_id): "Single sensor execution of dark current analysis." import lsst.eotest.sensor as sensorTest import siteUtils import eotestUtils file_prefix = '%s_%s' % (sensor_id, siteUtils.getRunNumber()) dark_files = siteUtils.dependency_glob( 'S*/%s_dark_dark_*.fits' % sensor_id, jobname=siteUtils.getProcessName('dark_raft_acq'), description='Dark files:') bias_frame = siteUtils.dependency_glob('%s_sflat*median_bias.fits' % sensor_id, description='Super bias frame:')[0] mask_files = \ eotestUtils.glob_mask_files(pattern='%s_*mask.fits' % sensor_id) gains = eotestUtils.getSensorGains(jobname='fe55_raft_analysis', sensor_id=sensor_id) task = sensorTest.DarkCurrentTask() task.config.temp_set_point = -100. dark_curr_pixels, dark95s \ = task.run(sensor_id, dark_files, mask_files, gains, bias_frame=bias_frame) results_file \ = siteUtils.dependency_glob('%s_eotest_results.fits' % sensor_id, jobname='read_noise_raft')[0] plots = sensorTest.EOTestPlots(sensor_id, results_file=results_file) siteUtils.make_png_file(plots.total_noise, '%s_noise.png' % file_prefix, dark95s=dark95s)
def dark_current_task(run, det_name, dark_files, gains, mask_files=(), temp_set_point=-100., bias_frame=None): """Single sensor execution of the dark current task.""" file_prefix = make_file_prefix(run, det_name) task = sensorTest.DarkCurrentTask() task.config.temp_set_point = temp_set_point return task.run(file_prefix, dark_files, mask_files, gains, bias_frame=bias_frame, linearity_correction=get_nlc_func(det_name))
#!/usr/bin/env python import lsst.eotest.sensor as sensorTest import siteUtils import eotestUtils siteUtils.aggregate_job_ids() sensor_id = siteUtils.getUnitId() dark_files = siteUtils.dependency_glob( '*_dark_dark_*.fits', jobname=siteUtils.getProcessName('dark_acq'), description='Dark files:') mask_files = eotestUtils.glob_mask_files() gains = eotestUtils.getSensorGains() task = sensorTest.DarkCurrentTask() task.run(sensor_id, dark_files, mask_files, gains)