Move P_int, PA_int computation in pipeline, replot NGC1068, MKN463, IC5063

This commit is contained in:
Thibault Barnouin
2022-03-28 10:42:38 +02:00
parent 453afca32c
commit 8d4c33603a
156 changed files with 86 additions and 138 deletions

View File

@@ -126,6 +126,10 @@ def save_Stokes(I_stokes, Q_stokes, U_stokes, Stokes_cov, P, debiased_P, s_P,
header['targname'] = (ref_header['targname'], 'Target name')
header['orientat'] = (ref_header['orientat'], 'Angle between North and the y-axis of the image')
header['filename'] = (filename, 'Original filename')
header['P_int'] = (ref_header['P_int'], 'Integrated polarization degree')
header['P_int_err'] = (ref_header['P_int_err'], 'Integrated polarization degree error')
header['PA_int'] = (ref_header['PA_int'], 'Integrated polarization angle')
header['PA_int_err'] = (ref_header['PA_int_err'], 'Integrated polarization angle error')
#Create HDUList object
hdul = fits.HDUList([])