diff --git a/plots/3C405_x136060/3C405_FOC.png b/plots/3C405_x136060/3C405_FOC.png index c39c05f..dbd9734 100644 Binary files a/plots/3C405_x136060/3C405_FOC.png and b/plots/3C405_x136060/3C405_FOC.png differ diff --git a/plots/3C405_x136060/3C405_FOC_IQU.png b/plots/3C405_x136060/3C405_FOC_IQU.png index 7b226f4..73dc614 100644 Binary files a/plots/3C405_x136060/3C405_FOC_IQU.png and b/plots/3C405_x136060/3C405_FOC_IQU.png differ diff --git a/plots/3C405_x136060/3C405_FOC_P.png b/plots/3C405_x136060/3C405_FOC_P.png index 882e97c..c0b2bb7 100644 Binary files a/plots/3C405_x136060/3C405_FOC_P.png and b/plots/3C405_x136060/3C405_FOC_P.png differ diff --git a/plots/3C405_x136060/3C405_FOC_P_err.png b/plots/3C405_x136060/3C405_FOC_P_err.png index 2aa4b41..3103e98 100644 Binary files a/plots/3C405_x136060/3C405_FOC_P_err.png and b/plots/3C405_x136060/3C405_FOC_P_err.png differ diff --git a/plots/3C405_x136060/3C405_FOC_SNRi.png b/plots/3C405_x136060/3C405_FOC_SNRi.png index 1cc1b86..a159ed1 100644 Binary files a/plots/3C405_x136060/3C405_FOC_SNRi.png and b/plots/3C405_x136060/3C405_FOC_SNRi.png differ diff --git a/plots/3C405_x136060/3C405_FOC_SNRp.png b/plots/3C405_x136060/3C405_FOC_SNRp.png index 9ed6a7c..b720859 100644 Binary files a/plots/3C405_x136060/3C405_FOC_SNRp.png and b/plots/3C405_x136060/3C405_FOC_SNRp.png differ diff --git a/plots/3C405_x136060/3C405_FOC_combine_FWHM050_rot.png b/plots/3C405_x136060/3C405_FOC_combine_FWHM050_rot.png index 08f2d88..db4400f 100644 Binary files a/plots/3C405_x136060/3C405_FOC_combine_FWHM050_rot.png and b/plots/3C405_x136060/3C405_FOC_combine_FWHM050_rot.png differ diff --git a/plots/3C405_x136060/3C405_FOC_combine_FWHM050_rot_IQU.png b/plots/3C405_x136060/3C405_FOC_combine_FWHM050_rot_IQU.png index 6e85cba..7971e2b 100644 Binary files a/plots/3C405_x136060/3C405_FOC_combine_FWHM050_rot_IQU.png and b/plots/3C405_x136060/3C405_FOC_combine_FWHM050_rot_IQU.png differ diff --git a/plots/3C405_x136060/3C405_FOC_combine_FWHM050_rot_P.png b/plots/3C405_x136060/3C405_FOC_combine_FWHM050_rot_P.png index 12f1650..449886d 100644 Binary files a/plots/3C405_x136060/3C405_FOC_combine_FWHM050_rot_P.png and b/plots/3C405_x136060/3C405_FOC_combine_FWHM050_rot_P.png differ diff --git a/plots/3C405_x136060/3C405_FOC_combine_FWHM050_rot_P_err.png b/plots/3C405_x136060/3C405_FOC_combine_FWHM050_rot_P_err.png index 2d2b3c0..bd39e1e 100644 Binary files a/plots/3C405_x136060/3C405_FOC_combine_FWHM050_rot_P_err.png and b/plots/3C405_x136060/3C405_FOC_combine_FWHM050_rot_P_err.png differ diff --git a/plots/3C405_x136060/3C405_FOC_combine_FWHM050_rot_SNRi.png b/plots/3C405_x136060/3C405_FOC_combine_FWHM050_rot_SNRi.png index e2a7dbe..4530fb1 100644 Binary files a/plots/3C405_x136060/3C405_FOC_combine_FWHM050_rot_SNRi.png and b/plots/3C405_x136060/3C405_FOC_combine_FWHM050_rot_SNRi.png differ diff --git a/plots/3C405_x136060/3C405_FOC_combine_FWHM050_rot_SNRp.png b/plots/3C405_x136060/3C405_FOC_combine_FWHM050_rot_SNRp.png index 292ba3d..77c47e8 100644 Binary files a/plots/3C405_x136060/3C405_FOC_combine_FWHM050_rot_SNRp.png and b/plots/3C405_x136060/3C405_FOC_combine_FWHM050_rot_SNRp.png differ diff --git a/plots/3C405_x136060/3C405_FOC_combine_FWHM100_rot.png b/plots/3C405_x136060/3C405_FOC_combine_FWHM100_rot.png index adf77a5..295bddb 100644 Binary files a/plots/3C405_x136060/3C405_FOC_combine_FWHM100_rot.png and b/plots/3C405_x136060/3C405_FOC_combine_FWHM100_rot.png differ diff --git a/plots/3C405_x136060/3C405_FOC_combine_FWHM100_rot_IQU.png b/plots/3C405_x136060/3C405_FOC_combine_FWHM100_rot_IQU.png index fbed9be..8a0c76e 100644 Binary files a/plots/3C405_x136060/3C405_FOC_combine_FWHM100_rot_IQU.png and b/plots/3C405_x136060/3C405_FOC_combine_FWHM100_rot_IQU.png differ diff --git a/plots/3C405_x136060/3C405_FOC_combine_FWHM100_rot_P.png b/plots/3C405_x136060/3C405_FOC_combine_FWHM100_rot_P.png index f5d64e0..69bf678 100644 Binary files a/plots/3C405_x136060/3C405_FOC_combine_FWHM100_rot_P.png and b/plots/3C405_x136060/3C405_FOC_combine_FWHM100_rot_P.png differ diff --git a/plots/3C405_x136060/3C405_FOC_combine_FWHM100_rot_P_err.png b/plots/3C405_x136060/3C405_FOC_combine_FWHM100_rot_P_err.png index f1bf6bb..9fae32c 100644 Binary files a/plots/3C405_x136060/3C405_FOC_combine_FWHM100_rot_P_err.png and b/plots/3C405_x136060/3C405_FOC_combine_FWHM100_rot_P_err.png differ diff --git a/plots/3C405_x136060/3C405_FOC_combine_FWHM100_rot_SNRi.png b/plots/3C405_x136060/3C405_FOC_combine_FWHM100_rot_SNRi.png index a4f66ed..b20fb2a 100644 Binary files a/plots/3C405_x136060/3C405_FOC_combine_FWHM100_rot_SNRi.png and b/plots/3C405_x136060/3C405_FOC_combine_FWHM100_rot_SNRi.png differ diff --git a/plots/3C405_x136060/3C405_FOC_combine_FWHM100_rot_SNRp.png b/plots/3C405_x136060/3C405_FOC_combine_FWHM100_rot_SNRp.png index 69369ed..b9c06fa 100644 Binary files a/plots/3C405_x136060/3C405_FOC_combine_FWHM100_rot_SNRp.png and b/plots/3C405_x136060/3C405_FOC_combine_FWHM100_rot_SNRp.png differ diff --git a/plots/IC5063_x3nl030/IC5063_FOC_combine_FWHM010_rot.png b/plots/IC5063_x3nl030/IC5063_FOC_combine_FWHM010_rot.png index e9e3780..556924c 100644 Binary files a/plots/IC5063_x3nl030/IC5063_FOC_combine_FWHM010_rot.png and b/plots/IC5063_x3nl030/IC5063_FOC_combine_FWHM010_rot.png differ diff --git a/plots/IC5063_x3nl030/IC5063_FOC_combine_FWHM010_rot_P.png b/plots/IC5063_x3nl030/IC5063_FOC_combine_FWHM010_rot_P.png index 71c968f..f8ab48a 100644 Binary files a/plots/IC5063_x3nl030/IC5063_FOC_combine_FWHM010_rot_P.png and b/plots/IC5063_x3nl030/IC5063_FOC_combine_FWHM010_rot_P.png differ diff --git a/plots/IC5063_x3nl030/IC5063_FOC_combine_FWHM010_rot_P_err.png b/plots/IC5063_x3nl030/IC5063_FOC_combine_FWHM010_rot_P_err.png index e3db942..cd89985 100644 Binary files a/plots/IC5063_x3nl030/IC5063_FOC_combine_FWHM010_rot_P_err.png and b/plots/IC5063_x3nl030/IC5063_FOC_combine_FWHM010_rot_P_err.png differ diff --git a/plots/IC5063_x3nl030/IC5063_FOC_combine_FWHM010_rot_SNRi.png b/plots/IC5063_x3nl030/IC5063_FOC_combine_FWHM010_rot_SNRi.png index 4d05609..20548f1 100644 Binary files a/plots/IC5063_x3nl030/IC5063_FOC_combine_FWHM010_rot_SNRi.png and b/plots/IC5063_x3nl030/IC5063_FOC_combine_FWHM010_rot_SNRi.png differ diff --git a/plots/IC5063_x3nl030/IC5063_FOC_combine_FWHM010_rot_SNRp.png b/plots/IC5063_x3nl030/IC5063_FOC_combine_FWHM010_rot_SNRp.png index c64a399..b267ded 100644 Binary files a/plots/IC5063_x3nl030/IC5063_FOC_combine_FWHM010_rot_SNRp.png and b/plots/IC5063_x3nl030/IC5063_FOC_combine_FWHM010_rot_SNRp.png differ diff --git a/plots/IC5063_x3nl030/IC5063_FOC_combine_FWHM020_rot.png b/plots/IC5063_x3nl030/IC5063_FOC_combine_FWHM020_rot.png index f0c21d6..dee3920 100644 Binary files a/plots/IC5063_x3nl030/IC5063_FOC_combine_FWHM020_rot.png and b/plots/IC5063_x3nl030/IC5063_FOC_combine_FWHM020_rot.png differ diff --git a/plots/IC5063_x3nl030/IC5063_FOC_combine_FWHM020_rot_P.png b/plots/IC5063_x3nl030/IC5063_FOC_combine_FWHM020_rot_P.png index e1e571e..d3a78d1 100644 Binary files a/plots/IC5063_x3nl030/IC5063_FOC_combine_FWHM020_rot_P.png and b/plots/IC5063_x3nl030/IC5063_FOC_combine_FWHM020_rot_P.png differ diff --git a/plots/IC5063_x3nl030/IC5063_FOC_combine_FWHM020_rot_P_err.png b/plots/IC5063_x3nl030/IC5063_FOC_combine_FWHM020_rot_P_err.png index b193250..14bf9b8 100644 Binary files a/plots/IC5063_x3nl030/IC5063_FOC_combine_FWHM020_rot_P_err.png and b/plots/IC5063_x3nl030/IC5063_FOC_combine_FWHM020_rot_P_err.png differ diff --git a/plots/IC5063_x3nl030/IC5063_FOC_combine_FWHM020_rot_SNRi.png b/plots/IC5063_x3nl030/IC5063_FOC_combine_FWHM020_rot_SNRi.png index ed9421a..fbf2aa5 100644 Binary files a/plots/IC5063_x3nl030/IC5063_FOC_combine_FWHM020_rot_SNRi.png and b/plots/IC5063_x3nl030/IC5063_FOC_combine_FWHM020_rot_SNRi.png differ diff --git a/plots/IC5063_x3nl030/IC5063_FOC_combine_FWHM020_rot_SNRp.png b/plots/IC5063_x3nl030/IC5063_FOC_combine_FWHM020_rot_SNRp.png index 412455b..b8aa2c3 100644 Binary files a/plots/IC5063_x3nl030/IC5063_FOC_combine_FWHM020_rot_SNRp.png and b/plots/IC5063_x3nl030/IC5063_FOC_combine_FWHM020_rot_SNRp.png differ diff --git a/plots/MKN463_x2rp030/MKN463_FOC_combine_FWHM010_rot.png b/plots/MKN463_x2rp030/MKN463_FOC_combine_FWHM010_rot.png index 62d14ae..f4c1fb7 100644 Binary files a/plots/MKN463_x2rp030/MKN463_FOC_combine_FWHM010_rot.png and b/plots/MKN463_x2rp030/MKN463_FOC_combine_FWHM010_rot.png differ diff --git a/plots/MKN463_x2rp030/MKN463_FOC_combine_FWHM010_rot_P.png b/plots/MKN463_x2rp030/MKN463_FOC_combine_FWHM010_rot_P.png index 514bce7..e1d9baa 100644 Binary files a/plots/MKN463_x2rp030/MKN463_FOC_combine_FWHM010_rot_P.png and b/plots/MKN463_x2rp030/MKN463_FOC_combine_FWHM010_rot_P.png differ diff --git a/plots/MKN463_x2rp030/MKN463_FOC_combine_FWHM010_rot_P_err.png b/plots/MKN463_x2rp030/MKN463_FOC_combine_FWHM010_rot_P_err.png index 8ca56f2..444c506 100644 Binary files a/plots/MKN463_x2rp030/MKN463_FOC_combine_FWHM010_rot_P_err.png and b/plots/MKN463_x2rp030/MKN463_FOC_combine_FWHM010_rot_P_err.png differ diff --git a/plots/MKN463_x2rp030/MKN463_FOC_combine_FWHM010_rot_SNRi.png b/plots/MKN463_x2rp030/MKN463_FOC_combine_FWHM010_rot_SNRi.png index db1d89d..32b23f0 100644 Binary files a/plots/MKN463_x2rp030/MKN463_FOC_combine_FWHM010_rot_SNRi.png and b/plots/MKN463_x2rp030/MKN463_FOC_combine_FWHM010_rot_SNRi.png differ diff --git a/plots/MKN463_x2rp030/MKN463_FOC_combine_FWHM010_rot_SNRp.png b/plots/MKN463_x2rp030/MKN463_FOC_combine_FWHM010_rot_SNRp.png index 9952c3d..17ea8f8 100644 Binary files a/plots/MKN463_x2rp030/MKN463_FOC_combine_FWHM010_rot_SNRp.png and b/plots/MKN463_x2rp030/MKN463_FOC_combine_FWHM010_rot_SNRp.png differ diff --git a/plots/MKN463_x2rp030/MKN463_FOC_combine_FWHM020_rot.png b/plots/MKN463_x2rp030/MKN463_FOC_combine_FWHM020_rot.png index 22064b3..03b6d12 100644 Binary files a/plots/MKN463_x2rp030/MKN463_FOC_combine_FWHM020_rot.png and b/plots/MKN463_x2rp030/MKN463_FOC_combine_FWHM020_rot.png differ diff --git a/plots/MKN463_x2rp030/MKN463_FOC_combine_FWHM020_rot_P.png b/plots/MKN463_x2rp030/MKN463_FOC_combine_FWHM020_rot_P.png index 11b023c..6a1bf3c 100644 Binary files a/plots/MKN463_x2rp030/MKN463_FOC_combine_FWHM020_rot_P.png and b/plots/MKN463_x2rp030/MKN463_FOC_combine_FWHM020_rot_P.png differ diff --git a/plots/MKN463_x2rp030/MKN463_FOC_combine_FWHM020_rot_P_err.png b/plots/MKN463_x2rp030/MKN463_FOC_combine_FWHM020_rot_P_err.png index 70eb946..6440ef1 100644 Binary files a/plots/MKN463_x2rp030/MKN463_FOC_combine_FWHM020_rot_P_err.png and b/plots/MKN463_x2rp030/MKN463_FOC_combine_FWHM020_rot_P_err.png differ diff --git a/plots/MKN463_x2rp030/MKN463_FOC_combine_FWHM020_rot_SNRi.png b/plots/MKN463_x2rp030/MKN463_FOC_combine_FWHM020_rot_SNRi.png index 3aa5c64..a9e81f8 100644 Binary files a/plots/MKN463_x2rp030/MKN463_FOC_combine_FWHM020_rot_SNRi.png and b/plots/MKN463_x2rp030/MKN463_FOC_combine_FWHM020_rot_SNRi.png differ diff --git a/plots/MKN463_x2rp030/MKN463_FOC_combine_FWHM020_rot_SNRp.png b/plots/MKN463_x2rp030/MKN463_FOC_combine_FWHM020_rot_SNRp.png index 4e9a902..989ecad 100644 Binary files a/plots/MKN463_x2rp030/MKN463_FOC_combine_FWHM020_rot_SNRp.png and b/plots/MKN463_x2rp030/MKN463_FOC_combine_FWHM020_rot_SNRp.png differ diff --git a/plots/MKN78_x3nl020/MKN78_FOC_combine_FWHM010_rot.png b/plots/MKN78_x3nl020/MKN78_FOC_combine_FWHM010_rot.png index 40182b7..150af82 100644 Binary files a/plots/MKN78_x3nl020/MKN78_FOC_combine_FWHM010_rot.png and b/plots/MKN78_x3nl020/MKN78_FOC_combine_FWHM010_rot.png differ diff --git a/plots/MKN78_x3nl020/MKN78_FOC_combine_FWHM010_rot_P.png b/plots/MKN78_x3nl020/MKN78_FOC_combine_FWHM010_rot_P.png index 4a04a47..79a1de9 100644 Binary files a/plots/MKN78_x3nl020/MKN78_FOC_combine_FWHM010_rot_P.png and b/plots/MKN78_x3nl020/MKN78_FOC_combine_FWHM010_rot_P.png differ diff --git a/plots/MKN78_x3nl020/MKN78_FOC_combine_FWHM010_rot_P_err.png b/plots/MKN78_x3nl020/MKN78_FOC_combine_FWHM010_rot_P_err.png index aa4e678..fa9a036 100644 Binary files a/plots/MKN78_x3nl020/MKN78_FOC_combine_FWHM010_rot_P_err.png and b/plots/MKN78_x3nl020/MKN78_FOC_combine_FWHM010_rot_P_err.png differ diff --git a/plots/MKN78_x3nl020/MKN78_FOC_combine_FWHM010_rot_SNRi.png b/plots/MKN78_x3nl020/MKN78_FOC_combine_FWHM010_rot_SNRi.png index 5a7360e..ddb5ee5 100644 Binary files a/plots/MKN78_x3nl020/MKN78_FOC_combine_FWHM010_rot_SNRi.png and b/plots/MKN78_x3nl020/MKN78_FOC_combine_FWHM010_rot_SNRi.png differ diff --git a/plots/MKN78_x3nl020/MKN78_FOC_combine_FWHM010_rot_SNRp.png b/plots/MKN78_x3nl020/MKN78_FOC_combine_FWHM010_rot_SNRp.png index f7f8635..eca81b6 100644 Binary files a/plots/MKN78_x3nl020/MKN78_FOC_combine_FWHM010_rot_SNRp.png and b/plots/MKN78_x3nl020/MKN78_FOC_combine_FWHM010_rot_SNRp.png differ diff --git a/plots/MKN78_x3nl020/MKN78_FOC_combine_FWHM020_rot.png b/plots/MKN78_x3nl020/MKN78_FOC_combine_FWHM020_rot.png index 86901b2..076d686 100644 Binary files a/plots/MKN78_x3nl020/MKN78_FOC_combine_FWHM020_rot.png and b/plots/MKN78_x3nl020/MKN78_FOC_combine_FWHM020_rot.png differ diff --git a/plots/MKN78_x3nl020/MKN78_FOC_combine_FWHM020_rot_P.png b/plots/MKN78_x3nl020/MKN78_FOC_combine_FWHM020_rot_P.png index 453d319..06cb4bd 100644 Binary files a/plots/MKN78_x3nl020/MKN78_FOC_combine_FWHM020_rot_P.png and b/plots/MKN78_x3nl020/MKN78_FOC_combine_FWHM020_rot_P.png differ diff --git a/plots/MKN78_x3nl020/MKN78_FOC_combine_FWHM020_rot_P_err.png b/plots/MKN78_x3nl020/MKN78_FOC_combine_FWHM020_rot_P_err.png index 7e904aa..e126a6a 100644 Binary files a/plots/MKN78_x3nl020/MKN78_FOC_combine_FWHM020_rot_P_err.png and b/plots/MKN78_x3nl020/MKN78_FOC_combine_FWHM020_rot_P_err.png differ diff --git a/plots/MKN78_x3nl020/MKN78_FOC_combine_FWHM020_rot_SNRi.png b/plots/MKN78_x3nl020/MKN78_FOC_combine_FWHM020_rot_SNRi.png index d47a2b5..78eee7f 100644 Binary files a/plots/MKN78_x3nl020/MKN78_FOC_combine_FWHM020_rot_SNRi.png and b/plots/MKN78_x3nl020/MKN78_FOC_combine_FWHM020_rot_SNRi.png differ diff --git a/plots/MKN78_x3nl020/MKN78_FOC_combine_FWHM020_rot_SNRp.png b/plots/MKN78_x3nl020/MKN78_FOC_combine_FWHM020_rot_SNRp.png index 6b25635..fdc1174 100644 Binary files a/plots/MKN78_x3nl020/MKN78_FOC_combine_FWHM020_rot_SNRp.png and b/plots/MKN78_x3nl020/MKN78_FOC_combine_FWHM020_rot_SNRp.png differ diff --git a/plots/NGC1068_x14w010/NGC1068_FOC_combine_FWHM010_rot.png b/plots/NGC1068_x14w010/NGC1068_FOC_combine_FWHM010_rot.png index f696fd7..801db7b 100644 Binary files a/plots/NGC1068_x14w010/NGC1068_FOC_combine_FWHM010_rot.png and b/plots/NGC1068_x14w010/NGC1068_FOC_combine_FWHM010_rot.png differ diff --git a/plots/NGC1068_x14w010/NGC1068_FOC_combine_FWHM010_rot_P.png b/plots/NGC1068_x14w010/NGC1068_FOC_combine_FWHM010_rot_P.png index a10f9e8..7ffae4a 100644 Binary files a/plots/NGC1068_x14w010/NGC1068_FOC_combine_FWHM010_rot_P.png and b/plots/NGC1068_x14w010/NGC1068_FOC_combine_FWHM010_rot_P.png differ diff --git a/plots/NGC1068_x14w010/NGC1068_FOC_combine_FWHM010_rot_P_err.png b/plots/NGC1068_x14w010/NGC1068_FOC_combine_FWHM010_rot_P_err.png index aa27775..73df66f 100644 Binary files a/plots/NGC1068_x14w010/NGC1068_FOC_combine_FWHM010_rot_P_err.png and b/plots/NGC1068_x14w010/NGC1068_FOC_combine_FWHM010_rot_P_err.png differ diff --git a/plots/NGC1068_x14w010/NGC1068_FOC_combine_FWHM010_rot_SNRi.png b/plots/NGC1068_x14w010/NGC1068_FOC_combine_FWHM010_rot_SNRi.png index 25721be..915f081 100644 Binary files a/plots/NGC1068_x14w010/NGC1068_FOC_combine_FWHM010_rot_SNRi.png and b/plots/NGC1068_x14w010/NGC1068_FOC_combine_FWHM010_rot_SNRi.png differ diff --git a/plots/NGC1068_x14w010/NGC1068_FOC_combine_FWHM010_rot_SNRp.png b/plots/NGC1068_x14w010/NGC1068_FOC_combine_FWHM010_rot_SNRp.png index 81066de..dc8f6d3 100644 Binary files a/plots/NGC1068_x14w010/NGC1068_FOC_combine_FWHM010_rot_SNRp.png and b/plots/NGC1068_x14w010/NGC1068_FOC_combine_FWHM010_rot_SNRp.png differ diff --git a/plots/NGC1068_x14w010/NGC1068_FOC_combine_FWHM020_rot.png b/plots/NGC1068_x14w010/NGC1068_FOC_combine_FWHM020_rot.png index 46bc56c..37a2c6c 100644 Binary files a/plots/NGC1068_x14w010/NGC1068_FOC_combine_FWHM020_rot.png and b/plots/NGC1068_x14w010/NGC1068_FOC_combine_FWHM020_rot.png differ diff --git a/plots/NGC1068_x14w010/NGC1068_FOC_combine_FWHM020_rot_P.png b/plots/NGC1068_x14w010/NGC1068_FOC_combine_FWHM020_rot_P.png index b72f1c1..3387acf 100644 Binary files a/plots/NGC1068_x14w010/NGC1068_FOC_combine_FWHM020_rot_P.png and b/plots/NGC1068_x14w010/NGC1068_FOC_combine_FWHM020_rot_P.png differ diff --git a/plots/NGC1068_x14w010/NGC1068_FOC_combine_FWHM020_rot_P_err.png b/plots/NGC1068_x14w010/NGC1068_FOC_combine_FWHM020_rot_P_err.png index 484fcaf..d15cabd 100644 Binary files a/plots/NGC1068_x14w010/NGC1068_FOC_combine_FWHM020_rot_P_err.png and b/plots/NGC1068_x14w010/NGC1068_FOC_combine_FWHM020_rot_P_err.png differ diff --git a/plots/NGC1068_x14w010/NGC1068_FOC_combine_FWHM020_rot_SNRi.png b/plots/NGC1068_x14w010/NGC1068_FOC_combine_FWHM020_rot_SNRi.png index 41094c1..2f7e026 100644 Binary files a/plots/NGC1068_x14w010/NGC1068_FOC_combine_FWHM020_rot_SNRi.png and b/plots/NGC1068_x14w010/NGC1068_FOC_combine_FWHM020_rot_SNRi.png differ diff --git a/plots/NGC1068_x14w010/NGC1068_FOC_combine_FWHM020_rot_SNRp.png b/plots/NGC1068_x14w010/NGC1068_FOC_combine_FWHM020_rot_SNRp.png index c3bdfb2..39ba499 100644 Binary files a/plots/NGC1068_x14w010/NGC1068_FOC_combine_FWHM020_rot_SNRp.png and b/plots/NGC1068_x14w010/NGC1068_FOC_combine_FWHM020_rot_SNRp.png differ diff --git a/plots/NGC1068_x274020/NGC1068_FOC_combine_FWHM010_rot.png b/plots/NGC1068_x274020/NGC1068_FOC_combine_FWHM010_rot.png index bded0a4..c000abd 100644 Binary files a/plots/NGC1068_x274020/NGC1068_FOC_combine_FWHM010_rot.png and b/plots/NGC1068_x274020/NGC1068_FOC_combine_FWHM010_rot.png differ diff --git a/plots/NGC1068_x274020/NGC1068_FOC_combine_FWHM010_rot_P.png b/plots/NGC1068_x274020/NGC1068_FOC_combine_FWHM010_rot_P.png index 17b3041..744c37b 100644 Binary files a/plots/NGC1068_x274020/NGC1068_FOC_combine_FWHM010_rot_P.png and b/plots/NGC1068_x274020/NGC1068_FOC_combine_FWHM010_rot_P.png differ diff --git a/plots/NGC1068_x274020/NGC1068_FOC_combine_FWHM010_rot_P_err.png b/plots/NGC1068_x274020/NGC1068_FOC_combine_FWHM010_rot_P_err.png index e3c37f3..250366f 100644 Binary files a/plots/NGC1068_x274020/NGC1068_FOC_combine_FWHM010_rot_P_err.png and b/plots/NGC1068_x274020/NGC1068_FOC_combine_FWHM010_rot_P_err.png differ diff --git a/plots/NGC1068_x274020/NGC1068_FOC_combine_FWHM010_rot_SNRi.png b/plots/NGC1068_x274020/NGC1068_FOC_combine_FWHM010_rot_SNRi.png index 171fb23..a5fff81 100644 Binary files a/plots/NGC1068_x274020/NGC1068_FOC_combine_FWHM010_rot_SNRi.png and b/plots/NGC1068_x274020/NGC1068_FOC_combine_FWHM010_rot_SNRi.png differ diff --git a/plots/NGC1068_x274020/NGC1068_FOC_combine_FWHM010_rot_SNRp.png b/plots/NGC1068_x274020/NGC1068_FOC_combine_FWHM010_rot_SNRp.png index 55f725b..514ed0c 100644 Binary files a/plots/NGC1068_x274020/NGC1068_FOC_combine_FWHM010_rot_SNRp.png and b/plots/NGC1068_x274020/NGC1068_FOC_combine_FWHM010_rot_SNRp.png differ diff --git a/plots/NGC1068_x274020/NGC1068_FOC_combine_FWHM020_rot.png b/plots/NGC1068_x274020/NGC1068_FOC_combine_FWHM020_rot.png index bfcf8ab..922e90e 100644 Binary files a/plots/NGC1068_x274020/NGC1068_FOC_combine_FWHM020_rot.png and b/plots/NGC1068_x274020/NGC1068_FOC_combine_FWHM020_rot.png differ diff --git a/plots/NGC1068_x274020/NGC1068_FOC_combine_FWHM020_rot_P.png b/plots/NGC1068_x274020/NGC1068_FOC_combine_FWHM020_rot_P.png index fabe03d..36c4949 100644 Binary files a/plots/NGC1068_x274020/NGC1068_FOC_combine_FWHM020_rot_P.png and b/plots/NGC1068_x274020/NGC1068_FOC_combine_FWHM020_rot_P.png differ diff --git a/plots/NGC1068_x274020/NGC1068_FOC_combine_FWHM020_rot_P_err.png b/plots/NGC1068_x274020/NGC1068_FOC_combine_FWHM020_rot_P_err.png index 74f2900..43d9122 100644 Binary files a/plots/NGC1068_x274020/NGC1068_FOC_combine_FWHM020_rot_P_err.png and b/plots/NGC1068_x274020/NGC1068_FOC_combine_FWHM020_rot_P_err.png differ diff --git a/plots/NGC1068_x274020/NGC1068_FOC_combine_FWHM020_rot_SNRi.png b/plots/NGC1068_x274020/NGC1068_FOC_combine_FWHM020_rot_SNRi.png index e8ea78f..21ce492 100644 Binary files a/plots/NGC1068_x274020/NGC1068_FOC_combine_FWHM020_rot_SNRi.png and b/plots/NGC1068_x274020/NGC1068_FOC_combine_FWHM020_rot_SNRi.png differ diff --git a/plots/NGC1068_x274020/NGC1068_FOC_combine_FWHM020_rot_SNRp.png b/plots/NGC1068_x274020/NGC1068_FOC_combine_FWHM020_rot_SNRp.png index 3545955..24b5870 100644 Binary files a/plots/NGC1068_x274020/NGC1068_FOC_combine_FWHM020_rot_SNRp.png and b/plots/NGC1068_x274020/NGC1068_FOC_combine_FWHM020_rot_SNRp.png differ diff --git a/src/FOC_reduction.py b/src/FOC_reduction.py index 64a6ad2..59b3612 100755 --- a/src/FOC_reduction.py +++ b/src/FOC_reduction.py @@ -17,11 +17,11 @@ import lib.plots as proj_plots #Functions for plotting data def main(): ##### User inputs ## Input and output locations - globals()['data_folder'] = "../data/NGC1068_x274020/" - infiles = ['x274020at.c0f.fits','x274020bt.c0f.fits','x274020ct.c0f.fits', - 'x274020dt.c0f.fits','x274020et.c0f.fits','x274020ft.c0f.fits', - 'x274020gt.c0f.fits','x274020ht.c0f.fits','x274020it.c0f.fits'] - globals()['plots_folder'] = "../plots/NGC1068_x274020/" +# globals()['data_folder'] = "../data/NGC1068_x274020/" +# infiles = ['x274020at.c0f.fits','x274020bt.c0f.fits','x274020ct.c0f.fits', +# 'x274020dt.c0f.fits','x274020et.c0f.fits','x274020ft.c0f.fits', +# 'x274020gt.c0f.fits','x274020ht.c0f.fits','x274020it.c0f.fits'] +# globals()['plots_folder'] = "../plots/NGC1068_x274020/" # globals()['data_folder'] = "../data/NGC1068_x14w010/" # infiles = ['x14w0101t_c0f.fits','x14w0102t_c0f.fits','x14w0103t_c0f.fits', @@ -44,13 +44,13 @@ def main(): # infiles = ['x3mc0101m_c0f.fits','x3mc0102m_c0f.fits','x3mc0103m_c0f.fits'] # globals()['plots_folder'] = "../plots/3C109_x3mc010/" -# globals()['data_folder'] = "../data/MKN463_x2rp030/" -# infiles = ['x2rp0201t_c0f.fits', 'x2rp0202t_c0f.fits', 'x2rp0203t_c0f.fits', -# 'x2rp0204t_c0f.fits', 'x2rp0205t_c0f.fits', 'x2rp0206t_c0f.fits', -# 'x2rp0207t_c0f.fits', 'x2rp0301t_c0f.fits', 'x2rp0302t_c0f.fits', -# 'x2rp0303t_c0f.fits', 'x2rp0304t_c0f.fits', 'x2rp0305t_c0f.fits', -# 'x2rp0306t_c0f.fits', 'x2rp0307t_c0f.fits'] -# globals()['plots_folder'] = "../plots/MKN463_x2rp030/" + globals()['data_folder'] = "../data/MKN463_x2rp030/" + infiles = ['x2rp0201t_c0f.fits', 'x2rp0202t_c0f.fits', 'x2rp0203t_c0f.fits', + 'x2rp0204t_c0f.fits', 'x2rp0205t_c0f.fits', 'x2rp0206t_c0f.fits', + 'x2rp0207t_c0f.fits', 'x2rp0301t_c0f.fits', 'x2rp0302t_c0f.fits', + 'x2rp0303t_c0f.fits', 'x2rp0304t_c0f.fits', 'x2rp0305t_c0f.fits', + 'x2rp0306t_c0f.fits', 'x2rp0307t_c0f.fits'] + globals()['plots_folder'] = "../plots/MKN463_x2rp030/" # globals()['data_folder'] = "../data/PG1630+377_x39510/" # infiles = ['x3990201m_c0f.fits', 'x3990205m_c0f.fits', 'x3995101r_c0f.fits', @@ -108,10 +108,10 @@ def main(): rotate_stokes = True #rotation to North convention can give erroneous results rotate_data = False #rotation to North convention can give erroneous results # Polarization map output - figname = 'NGC1068_FOC' #target/intrument name + figname = 'MKN463_FOC' #target/intrument name figtype = '_combine_FWHM020_rot' #additionnal informations - SNRp_cut = 20 #P measurments with SNR>3 - SNRi_cut = 130 #I measurments with SNR>30, which implies an uncertainty in P of 4.7%. + SNRp_cut = 3 #P measurments with SNR>3 + SNRi_cut = 30 #I measurments with SNR>30, which implies an uncertainty in P of 4.7%. step_vec = 1 #plot all vectors in the array. if step_vec = 2, then every other vector will be plotted ##### Pipeline start diff --git a/src/lib/reduction.py b/src/lib/reduction.py index 027809b..44cf696 100755 --- a/src/lib/reduction.py +++ b/src/lib/reduction.py @@ -800,7 +800,7 @@ def smooth_data(data_array, error_array, data_mask, headers, FWHM=1., raise ValueError("Not all images in array have same pixel size") FWHM /= pxsize[0].min() - #Define gaussian stdev + # Define gaussian stdev stdev = FWHM/(2.*np.sqrt(2.*np.log(2))) fmax = np.finfo(np.float64).max @@ -821,15 +821,16 @@ def smooth_data(data_array, error_array, data_mask, headers, FWHM=1., dist_rc = np.where(data_mask, fmax, np.sqrt((r-xx)**2+(c-yy)**2)) g_rc = np.array([np.exp(-0.5*(dist_rc/stdev)**2),]*len(data_array)) # Apply weighted combination - if data_mask[r,c]: - smoothed[r,c] = 0. - error[r,c] = 1. - else: - smoothed[r,c] = np.sum(data_array*weight*g_rc)/np.sum(weight*g_rc) - error[r,c] = np.sqrt(np.sum(weight*g_rc**2))/np.sum(weight*g_rc) + smoothed[r,c] = (1.-data_mask[r,c])*np.sum(data_array*weight*g_rc)/np.sum(weight*g_rc) + error[r,c] = np.sqrt(np.sum(weight*g_rc**2))/np.sum(weight*g_rc) + + # Nan handling + error[np.isnan(smoothed)] = 0. + smoothed[np.isnan(smoothed)] = 0. + error[np.isnan(error)] = 0. elif smoothing.lower() in ['gauss','gaussian']: - #Convolution with gaussian function + # Convolution with gaussian function smoothed = np.zeros(data_array.shape) error = np.zeros(error_array.shape) for i,image in enumerate(data_array): @@ -838,12 +839,13 @@ def smooth_data(data_array, error_array, data_mask, headers, FWHM=1., for c in range(image.shape[1]): dist_rc = np.where(data_mask, fmax, np.sqrt((r-xx)**2+(c-yy)**2)) g_rc = np.exp(-0.5*(dist_rc/stdev)**2)/(2.*np.pi*stdev**2) - if data_mask[r,c]: - smoothed[r,c] = 0. - error[r,c] = 1. - else: - smoothed[r,c] = np.sum(data_array*weight*g_rc)/np.sum(weight*g_rc) - error[r,c] = np.sqrt(np.sum(weight*g_rc**2))/np.sum(weight*g_rc) + smoothed[i][r,c] = (1.-data_mask[r,c])*np.sum(data_array*weight*g_rc)/np.sum(weight*g_rc) + error[i][r,c] = np.sqrt(np.sum(weight*g_rc**2))/np.sum(weight*g_rc) + + # Nan handling + error[i][np.isnan(smoothed)] = 0. + smoothed[i][np.isnan(smoothed)] = 0. + error[i][np.isnan(error)] = 0. else: raise ValueError("{} is not a valid smoothing option".format(smoothing))