#!/usr/bin/env python
import logging
import os.path as op
from ..stpipe import Pipeline
# step imports
from ..ami import ami_analyze_step
from ..ami import ami_average_step
from ..ami import ami_normalize_step
from ..model_blender import blendmeta
__all__ = ['Ami3Pipeline']
# Define logging
log = logging.getLogger(__name__)
log.setLevel(logging.DEBUG)
[docs]
class Ami3Pipeline(Pipeline):
"""
Ami3Pipeline: Apply all level-3 calibration steps to an
association of level-2b AMI exposures. Included steps are:
ami_analyze (fringe detection)
ami_average (average results of fringe detection)
ami_normalize (normalize results by reference target)
"""
class_alias = "calwebb_ami3"
spec = """
save_averages = boolean(default=True)
"""
# Define aliases to steps
step_defs = {'ami_analyze': ami_analyze_step.AmiAnalyzeStep,
'ami_average': ami_average_step.AmiAverageStep,
'ami_normalize': ami_normalize_step.AmiNormalizeStep
}
[docs]
def process(self, input):
log.info('Starting calwebb_ami3')
# Load the input association table
asn = self.load_as_level3_asn(input)
# We assume there's one final product defined by the association
asn_id = asn['asn_id']
prod = asn['products'][0]
self.output_file = prod.get('name', self.output_file)
# Construct lists of all the PSF and science target members
# so that we know what we've been given to work with
psf_files = [
m['expname']
for m in prod['members']
if m['exptype'].upper() == 'PSF'
]
targ_files = [
m['expname']
for m in prod['members']
if m['exptype'].upper() == 'SCIENCE'
]
# Make sure we found some science target members
if len(targ_files) == 0:
log.error('No science target members found in association table')
log.error('Calwebb_ami3 processing will be aborted')
return
# If there aren't any PSF images, we can't do the normalize step
if len(psf_files) == 0:
log.info('No PSF reference members found in association table')
log.info('ami_normalize step will be skipped')
# Run ami_analyze on all the target members
targ_lg = []
for input_file in targ_files:
# Do the LG analysis for this image
log.debug('Do LG processing for member %s', input_file)
result = self.ami_analyze(input_file)
# Save the LG analysis results to a file
result.meta.asn.pool_name = asn['asn_pool']
result.meta.asn.table_name = op.basename(asn.filename)
self.save_model(result, output_file=input_file, suffix='ami', asn_id=asn_id)
# Save the result for use as input to ami_average
targ_lg.append(result)
# Run ami_analyze on all the psf members
psf_lg = []
for input_file in psf_files:
# Do the LG analysis for this image
log.debug('Do LG processing for member %s', input_file)
result = self.ami_analyze(input_file)
# Save the LG analysis results to a file
result.meta.asn.pool_name = asn['asn_pool']
result.meta.asn.table_name = op.basename(asn.filename)
self.save_model(result, output_file=input_file, suffix='ami', asn_id=asn_id)
# Save the result for use as input to ami_average
psf_lg.append(result)
# Average the reference PSF image results
psf_avg = None
if len(psf_files) > 0:
self.log.debug('Calling ami_average for PSF results ...')
psf_avg = self.ami_average(psf_lg)
# Save the results to a file, if requested
if self.save_averages:
psf_avg.meta.asn.pool_name = asn['asn_pool']
psf_avg.meta.asn.table_name = op.basename(asn.filename)
# Perform blending of metadata for all inputs to this
# output file
self.log.info('Blending metadata for averaged psf')
blendmeta.blendmodels(psf_avg, inputs=psf_lg)
self.save_model(psf_avg, suffix='psf-amiavg')
del psf_lg
# Average the science target image results
if len(targ_files) > 0:
self.log.debug('Calling ami_average for target results ...')
targ_avg = self.ami_average(targ_lg)
# Save the results to a file, if requested
if self.save_averages:
targ_avg.meta.asn.pool_name = asn['asn_pool']
targ_avg.meta.asn.table_name = op.basename(asn.filename)
# Perform blending of metadata for all inputs to this
# output file
self.log.info('Blending metadata for averaged target')
blendmeta.blendmodels(targ_avg, inputs=targ_lg)
self.save_model(targ_avg, suffix='amiavg')
del targ_lg
# Now that all LGAVG products have been produced, do
# normalization of the target results by the reference
# results, if reference results exist
if psf_avg is not None:
result = self.ami_normalize(targ_avg, psf_avg)
# Save the result
result.meta.asn.pool_name = asn['asn_pool']
result.meta.asn.table_name = op.basename(asn.filename)
# Perform blending of metadata for all inputs to this output file
self.log.info('Blending metadata for PSF normalized target')
blendmeta.blendmodels(result, inputs=[targ_avg, psf_avg])
self.save_model(result, suffix='aminorm')
result.close()
del psf_avg
del targ_avg
# We're done
log.info('... ending calwebb_ami3')
return