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)