diff --git a/data/3C273_x0u20/README.md b/data/3C273_x0u20/README.md new file mode 100755 index 0000000..c264204 --- /dev/null +++ b/data/3C273_x0u20/README.md @@ -0,0 +1,5 @@ +Observation data +Instrument : FOC +Proposal ID : 2451 +Target name : 3C273 +Dataset : X0U20 diff --git a/plots/3C273_x0u20/3C2173_FOC.png b/plots/3C273_x0u20/3C2173_FOC.png new file mode 100644 index 0000000..da588b2 Binary files /dev/null and b/plots/3C273_x0u20/3C2173_FOC.png differ diff --git a/plots/3C273_x0u20/3C2173_FOC_IQU.png b/plots/3C273_x0u20/3C2173_FOC_IQU.png new file mode 100644 index 0000000..744855a Binary files /dev/null and b/plots/3C273_x0u20/3C2173_FOC_IQU.png differ diff --git a/plots/3C273_x0u20/3C2173_FOC_POL0_crop_region.png b/plots/3C273_x0u20/3C2173_FOC_POL0_crop_region.png new file mode 100644 index 0000000..54a20a4 Binary files /dev/null and b/plots/3C273_x0u20/3C2173_FOC_POL0_crop_region.png differ diff --git a/plots/3C273_x0u20/3C2173_FOC_center_image.png b/plots/3C273_x0u20/3C2173_FOC_center_image.png new file mode 100644 index 0000000..93c4d88 Binary files /dev/null and b/plots/3C273_x0u20/3C2173_FOC_center_image.png differ diff --git a/plots/3C273_x0u20/3C2173_FOC_combine_FWHM010.png b/plots/3C273_x0u20/3C2173_FOC_combine_FWHM010.png new file mode 100644 index 0000000..fe359da Binary files /dev/null and b/plots/3C273_x0u20/3C2173_FOC_combine_FWHM010.png differ diff --git a/plots/3C273_x0u20/3C2173_FOC_combine_FWHM010_IQU.png b/plots/3C273_x0u20/3C2173_FOC_combine_FWHM010_IQU.png new file mode 100644 index 0000000..a4c3c6a Binary files /dev/null and b/plots/3C273_x0u20/3C2173_FOC_combine_FWHM010_IQU.png differ diff --git a/plots/3C273_x0u20/3C2173_FOC_combine_FWHM010_I_err.png b/plots/3C273_x0u20/3C2173_FOC_combine_FWHM010_I_err.png new file mode 100644 index 0000000..f2bfed7 Binary files /dev/null and b/plots/3C273_x0u20/3C2173_FOC_combine_FWHM010_I_err.png differ diff --git a/plots/3C273_x0u20/3C2173_FOC_combine_FWHM010_P.png b/plots/3C273_x0u20/3C2173_FOC_combine_FWHM010_P.png new file mode 100644 index 0000000..4bd0e88 Binary files /dev/null and b/plots/3C273_x0u20/3C2173_FOC_combine_FWHM010_P.png differ diff --git a/plots/3C273_x0u20/3C2173_FOC_combine_FWHM010_POL0_crop_region.png b/plots/3C273_x0u20/3C2173_FOC_combine_FWHM010_POL0_crop_region.png new file mode 100644 index 0000000..aef990a Binary files /dev/null and b/plots/3C273_x0u20/3C2173_FOC_combine_FWHM010_POL0_crop_region.png differ diff --git a/plots/3C273_x0u20/3C2173_FOC_combine_FWHM010_P_err.png b/plots/3C273_x0u20/3C2173_FOC_combine_FWHM010_P_err.png new file mode 100644 index 0000000..4b2a69d Binary files /dev/null and b/plots/3C273_x0u20/3C2173_FOC_combine_FWHM010_P_err.png differ diff --git a/plots/3C273_x0u20/3C2173_FOC_combine_FWHM010_P_flux.png b/plots/3C273_x0u20/3C2173_FOC_combine_FWHM010_P_flux.png new file mode 100644 index 0000000..d974d5d Binary files /dev/null and b/plots/3C273_x0u20/3C2173_FOC_combine_FWHM010_P_flux.png differ diff --git a/plots/3C273_x0u20/3C2173_FOC_combine_FWHM010_SNRi.png b/plots/3C273_x0u20/3C2173_FOC_combine_FWHM010_SNRi.png new file mode 100644 index 0000000..9db04c7 Binary files /dev/null and b/plots/3C273_x0u20/3C2173_FOC_combine_FWHM010_SNRi.png differ diff --git a/plots/3C273_x0u20/3C2173_FOC_combine_FWHM010_SNRp.png b/plots/3C273_x0u20/3C2173_FOC_combine_FWHM010_SNRp.png new file mode 100644 index 0000000..157461d Binary files /dev/null and b/plots/3C273_x0u20/3C2173_FOC_combine_FWHM010_SNRp.png differ diff --git a/plots/3C273_x0u20/3C2173_FOC_combine_FWHM010_center_image.png b/plots/3C273_x0u20/3C2173_FOC_combine_FWHM010_center_image.png new file mode 100644 index 0000000..cb1191d Binary files /dev/null and b/plots/3C273_x0u20/3C2173_FOC_combine_FWHM010_center_image.png differ diff --git a/plots/3C273_x0u20/3C2173_FOC_combine_FWHM010_crop_region.png b/plots/3C273_x0u20/3C2173_FOC_combine_FWHM010_crop_region.png new file mode 100644 index 0000000..992324f Binary files /dev/null and b/plots/3C273_x0u20/3C2173_FOC_combine_FWHM010_crop_region.png differ diff --git a/plots/3C273_x0u20/3C2173_FOC_combine_FWHM010_deconvolved.png b/plots/3C273_x0u20/3C2173_FOC_combine_FWHM010_deconvolved.png new file mode 100644 index 0000000..ababdea Binary files /dev/null and b/plots/3C273_x0u20/3C2173_FOC_combine_FWHM010_deconvolved.png differ diff --git a/plots/3C273_x0u20/3C2173_FOC_combine_FWHM010_deconvolved_IQU.png b/plots/3C273_x0u20/3C2173_FOC_combine_FWHM010_deconvolved_IQU.png new file mode 100644 index 0000000..18869e8 Binary files /dev/null and b/plots/3C273_x0u20/3C2173_FOC_combine_FWHM010_deconvolved_IQU.png differ diff --git a/plots/3C273_x0u20/3C2173_FOC_combine_FWHM010_deconvolved_I_err.png b/plots/3C273_x0u20/3C2173_FOC_combine_FWHM010_deconvolved_I_err.png new file mode 100644 index 0000000..3bb859f Binary files /dev/null and b/plots/3C273_x0u20/3C2173_FOC_combine_FWHM010_deconvolved_I_err.png differ diff --git a/plots/3C273_x0u20/3C2173_FOC_combine_FWHM010_deconvolved_P.png b/plots/3C273_x0u20/3C2173_FOC_combine_FWHM010_deconvolved_P.png new file mode 100644 index 0000000..130eec4 Binary files /dev/null and b/plots/3C273_x0u20/3C2173_FOC_combine_FWHM010_deconvolved_P.png differ diff --git a/plots/3C273_x0u20/3C2173_FOC_combine_FWHM010_deconvolved_P_err.png b/plots/3C273_x0u20/3C2173_FOC_combine_FWHM010_deconvolved_P_err.png new file mode 100644 index 0000000..1243188 Binary files /dev/null and b/plots/3C273_x0u20/3C2173_FOC_combine_FWHM010_deconvolved_P_err.png differ diff --git a/plots/3C273_x0u20/3C2173_FOC_combine_FWHM010_deconvolved_P_flux.png b/plots/3C273_x0u20/3C2173_FOC_combine_FWHM010_deconvolved_P_flux.png new file mode 100644 index 0000000..7cce111 Binary files /dev/null and b/plots/3C273_x0u20/3C2173_FOC_combine_FWHM010_deconvolved_P_flux.png differ diff --git a/plots/3C273_x0u20/3C2173_FOC_combine_FWHM010_deconvolved_SNRi.png b/plots/3C273_x0u20/3C2173_FOC_combine_FWHM010_deconvolved_SNRi.png new file mode 100644 index 0000000..ef62a2b Binary files /dev/null and b/plots/3C273_x0u20/3C2173_FOC_combine_FWHM010_deconvolved_SNRi.png differ diff --git a/plots/3C273_x0u20/3C2173_FOC_combine_FWHM010_deconvolved_SNRp.png b/plots/3C273_x0u20/3C2173_FOC_combine_FWHM010_deconvolved_SNRp.png new file mode 100644 index 0000000..d9d6c47 Binary files /dev/null and b/plots/3C273_x0u20/3C2173_FOC_combine_FWHM010_deconvolved_SNRp.png differ diff --git a/plots/3C273_x0u20/3C2173_FOC_combine_FWHM010_errors_background_flux.png b/plots/3C273_x0u20/3C2173_FOC_combine_FWHM010_errors_background_flux.png new file mode 100644 index 0000000..4fe987a Binary files /dev/null and b/plots/3C273_x0u20/3C2173_FOC_combine_FWHM010_errors_background_flux.png differ diff --git a/plots/3C273_x0u20/3C2173_FOC_combine_FWHM010_errors_background_location.png b/plots/3C273_x0u20/3C2173_FOC_combine_FWHM010_errors_background_location.png new file mode 100644 index 0000000..42e8fd0 Binary files /dev/null and b/plots/3C273_x0u20/3C2173_FOC_combine_FWHM010_errors_background_location.png differ diff --git a/plots/3C273_x0u20/3C2173_FOC_crop_region.png b/plots/3C273_x0u20/3C2173_FOC_crop_region.png new file mode 100644 index 0000000..848b41b Binary files /dev/null and b/plots/3C273_x0u20/3C2173_FOC_crop_region.png differ diff --git a/plots/3C273_x0u20/3C2173_FOC_errors_background_flux.png b/plots/3C273_x0u20/3C2173_FOC_errors_background_flux.png new file mode 100644 index 0000000..ad21397 Binary files /dev/null and b/plots/3C273_x0u20/3C2173_FOC_errors_background_flux.png differ diff --git a/plots/3C273_x0u20/3C2173_FOC_errors_background_location.png b/plots/3C273_x0u20/3C2173_FOC_errors_background_location.png new file mode 100644 index 0000000..e0a9e78 Binary files /dev/null and b/plots/3C273_x0u20/3C2173_FOC_errors_background_location.png differ diff --git a/src/FOC_reduction.py b/src/FOC_reduction.py index 0534b98..75981f0 100755 --- a/src/FOC_reduction.py +++ b/src/FOC_reduction.py @@ -17,11 +17,11 @@ from lib.convex_hull import image_hull 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', @@ -81,9 +81,13 @@ def main(): # infiles = ['x25d0401t_c0f.fits','x25d0402t_c0f.fits','x25d0403t_c0f.fits'] # globals()['plots_folder'] = "../plots/PictorA_x25d040/" + globals()['data_folder'] = "../data/3C273_x0u20/" + infiles = ['x0u20101t_c0f.fits','x0u20102t_c0f.fits','x0u20103t_c0f.fits','x0u20104t_c0f.fits','x0u20105t_c0f.fits','x0u20106t_c0f.fits','x0u20201t_c0f.fits','x0u20202t_c0f.fits','x0u20203t_c0f.fits','x0u20204t_c0f.fits','x0u20205t_c0f.fits','x0u20206t_c0f.fits','x0u20301t_c0f.fits','x0u20302t_c0f.fits','x0u20303t_c0f.fits','x0u20304t_c0f.fits','x0u20305t_c0f.fits','x0u20306t_c0f.fits'] + globals()['plots_folder'] = "../plots/3C273_x0u20/" + ## Reduction parameters # Deconvolution - deconvolve = False + deconvolve = True if deconvolve: psf = 'gaussian' #Can be user-defined as well psf_FWHM = 0.10 @@ -106,16 +110,16 @@ def main(): display_data = False # Smoothing smoothing_function = 'combine' #gaussian_after, gaussian or combine - smoothing_FWHM = 0.20 #If None, no smoothing is done + smoothing_FWHM = 0.10 #If None, no smoothing is done smoothing_scale = 'arcsec' #pixel or arcsec # Rotation 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 - figtype = '_combine_FWHM020_rot_withaxiserror' #additionnal informations - SNRp_cut = 15. #P measurments with SNR>3 - SNRi_cut = 80. #I measurments with SNR>30, which implies an uncertainty in P of 4.7%. + figname = '3C2173_FOC' #target/intrument name + figtype = '_combine_FWHM010_deconvolved' #additionnal informations + SNRp_cut = 5. #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