Runs on Troxel et al. RomanDESCSims simulated Roman WFI images. Currently developed and tested at Duke Computing Cluster and Perlmutter at NERSC.
phrosty
roman_imsim
get_object_instances.pyto get all the images by filter, pointing, and SCA that contain the specified transient.preprocess.pyto create references, get PSFs, and do cross-convolutions.sfftdiff.pyto run SFFT subtraction on a GPU.postprocess.pyto generate decorrelation kernels, apply those to images, and make stamps.mk_lc.pyto extract a lightcurve from the subtractions.
sfft_and_animate_*.sh as a SLURM array job that has SN ID hardcoded.
This then calls get_object_instances, preprocess, sfftdiff, postprocess, and mk_lc with the SN ID (refered to as oid in the code).
sfft_and_animate_*.py reads the SLURM_ARRAY_JOBID from the environment and uses that to select the band to process. (1-7).
sfft_and_animate_*.py reads input data from input data as defined by phrosty.utils._build_filepath, which uses environment variableSIMS_DIR to look for the data. SIMS_DIR should match where the RomanDESCSims are.
preprocess.py, sfftdiff.py, and postprocess.py write files out to environment variable DIA_OUT_DIR.
mk_lc.py writes files out to environment variable LC_OUT_DIR.
Example usage at the DCC Run on SNID 20172782 (yes, the SN IDs are 8-digit numbers beginning with 20, they are not dates), filter R062, with specified input directory.
python mk_lc.py 20172782 --band R062 --inputdir /work/lna18/