diff --git a/plots/IC5063_x3nl030/IC5063_FOC_POL0_crop_region.png b/plots/IC5063_x3nl030/IC5063_FOC_POL0_crop_region.png index 524ed95..eb28d0a 100644 Binary files a/plots/IC5063_x3nl030/IC5063_FOC_POL0_crop_region.png and b/plots/IC5063_x3nl030/IC5063_FOC_POL0_crop_region.png differ diff --git a/plots/IC5063_x3nl030/IC5063_FOC_errors_POL0_background_location.png b/plots/IC5063_x3nl030/IC5063_FOC_errors_POL0_background_location.png index e7a9984..a98200f 100644 Binary files a/plots/IC5063_x3nl030/IC5063_FOC_errors_POL0_background_location.png and b/plots/IC5063_x3nl030/IC5063_FOC_errors_POL0_background_location.png differ diff --git a/plots/IC5063_x3nl030/IC5063_FOC_errors_background_flux.png b/plots/IC5063_x3nl030/IC5063_FOC_errors_background_flux.png index af6959b..dd40894 100644 Binary files a/plots/IC5063_x3nl030/IC5063_FOC_errors_background_flux.png and b/plots/IC5063_x3nl030/IC5063_FOC_errors_background_flux.png differ diff --git a/src/FOC_reduction.py b/src/FOC_reduction.py index c6d1386..b91afbb 100755 --- a/src/FOC_reduction.py +++ b/src/FOC_reduction.py @@ -20,17 +20,17 @@ from astropy.wcs import WCS 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'] - psf_file = 'NGC1068_f253m00.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'] +# psf_file = 'NGC1068_f253m00.fits' +# globals()['plots_folder'] = "../plots/NGC1068_x274020/" -# globals()['data_folder'] = "../data/IC5063_x3nl030/" -# infiles = ['x3nl0301r_c0f.fits','x3nl0302r_c0f.fits','x3nl0303r_c0f.fits'] -# psf_file = 'IC5063_f502m00.fits' -# globals()['plots_folder'] = "../plots/IC5063_x3nl030/" + globals()['data_folder'] = "../data/IC5063_x3nl030/" + infiles = ['x3nl0301r_c0f.fits','x3nl0302r_c0f.fits','x3nl0303r_c0f.fits'] + psf_file = 'IC5063_f502m00.fits' + globals()['plots_folder'] = "../plots/IC5063_x3nl030/" # globals()['data_folder'] = "../data/NGC1068_x14w010/" # infiles = ['x14w0101t_c0f.fits','x14w0102t_c0f.fits','x14w0103t_c0f.fits', @@ -122,7 +122,7 @@ def main(): crop = False #Crop to desired ROI final_display = False # Polarization map output - figname = 'NGC1068_FOC' #target/intrument name + figname = 'IC5063_FOC' #target/intrument name figtype = '_combine_FWHM020' #additionnal informations SNRp_cut = 5. #P measurments with SNR>3 SNRi_cut = 50. #I measurments with SNR>30, which implies an uncertainty in P of 4.7%.