Level 3 From Aligned Cals

We run F560W with the best options determined from the Improving Alignment example:

run_st_wcs_align_batch.py --input_dir '.' --input_files 'miri_example/*_cal.fits' --outrootdir aligned
        --outsubdir F560W_level2_gaia  --overwr -v --refcat gaia --saveplots -tt -pp --histocut_order dxdy
        --roundness1_lim -0.5 0.5 --objmag_lim 14 21.5  --refmag_lim 16 25 --delta_mag_lim -2 2
        --d2d_max 1.5 --xshift 3 --yshift -4 --filter F560W

Now we run the notebook that runs level3 and creates the mosaic and the corresponding catalog (On GitHub here). That notebook produces a catalog (F560W_snr3_npix10_cat.ecsv), which we choose as our secondary astrometric catalog.

Now we run the rest of the filters. However, we remove –delta_mag_lim, since this cut depends on the filter!!! Note that we define the necessary column names from the catalog, which are different from the defaults.

run_st_wcs_align_batch.py --input_dir '.' --input_files 'miri_example/*_cal.fits' --outrootdir miri_example
        --outsubdir ALL_level2_catF560W  --overwr -v --refcat F560W_snr3_npix10_cat.ecsv --saveplots -tt -pp
        --histocut_order dxdy --roundness1_lim -0.5 0.5 --objmag_lim 14 21.5  --refmag_lim 16 25 --d2d_max 1.5
        --xshift 3 --yshift -4 --iterate_with_xyshifts --refcat_racol sky_centroid.ra
        --refcat_deccol sky_centroid.dec --refcat_magcol aper50_abmag  --refcat_magerrcol aper50_abmag_err
        --filters F560W F1000W F1280W F1130W F1500W F1800W

We find that F1500W F1800W mostly work, but don’t have many stars. Therefore we run it with –d_rotated_Nsigma 0.0: Too few stars to do a 3-sigma cut.

run_st_wcs_align_batch.py --input_dir '.' --input_files 'miri_example/*_cal.fits' --outrootdir miri_example
        --outsubdir ALLRED_level2_catF560W  --overwr -vvv --refcat F560W_snr3_npix10_cat.ecsv --saveplots -tt -pp
        --histocut_order dxdy --roundness1_lim -0.5 0.5 --objmag_lim 14 21.5  --refmag_lim 16 25 --d2d_max 1.5 --xshift 3
        --yshift -4 --iterate_with_xyshifts --refcat_racol sky_centroid.ra --refcat_deccol sky_centroid.dec --refcat_magcol
        aper50_abmag  --refcat_magerrcol aper50_abmag_err  --filters F1500W F1800W F2100W --d_rotated_Nsigma 0.0