stokes[1]=kw_pol_stokes_n
hadmrFITS.write_newkey(filename,stokes)
elif what=='I':
stokes[1]=kw_pol_stokes_i
hadmrFITS.write_newkey(filename,stokes)
## RUN THE DEMODULATION
WLOG('info',log_opt,'Starting the demodulation calculation.')
if DOUB: # We have double set of files
X,I,N=ratio_double(data)
hadmrFITS.write_data(X_filename,X)
fix_headers(X_filename,STOKES)
hadmrFITS.write_data(I_filename,I)
fix_headers(I_filename,'I')
hadmrFITS.write_data(N_filename,N)
fix_headers(N_filename,'N')
X,I,N=diff_double(data)
hadmrFITS.write_data(X_filename.replace('.fits','_diff.fits'),X)
fix_headers(X_filename.replace('.fits','_diff.fits'),STOKES)
hadmrFITS.write_data(I_filename.replace('.fits','_diff.fits'),I)
fix_headers(I_filename.replace('.fits','_diff.fits'),'I')
hadmrFITS.write_data(N_filename.replace('.fits','_diff.fits'),N)
fix_headers(N_filename.replace('.fits','_diff.fits'),'N')
if dic_db.has_key('TH_'+fiber)==0:
WLOG('error',log_opt,'No wavelength Calibration for fiber '+fiber+' FATAL!')
exit(3)
th_file=os.path.join(dirfits,dic_db['TH_'+fiber][1])
wave[fiber],param_ll[fiber],param_x=hadmrFITS.get_e2ds_ll(th_file,_kw_TH_CAL_)
WLOG('',log_opt,'Wavelength calibration is set using '+os.path.split(th_file)[1])
#############################
# Save E2DSFF on fits file #
#############################
WLOG('',log_opt,'Saving E2DS spectrum of Fiber '+fiber+' in '+os.path.split(e2ds_fitsfilename[fiber])[1])
hadmrFITS.write_data(e2ds_fitsfilename[fiber],e2dsff[fiber])
# Copy keywords of the raw frame
hadmrFITS.copy_key(fitsfilename,e2ds_fitsfilename[fiber])
hadmrFITS.write_newkey(e2ds_fitsfilename[fiber],kw_version)
hadmrFITS.write_newkey(e2ds_fitsfilename[fiber],kw_CCD_SIGDET)
hadmrFITS.write_newkey(e2ds_fitsfilename[fiber],kw_CCD_CONAD)
# Copy keywords of the localization frame
kw_LOCO_FILE[1]=dic_db['LOC_'+fiber][1]
hadmrFITS.write_newkey(e2ds_fitsfilename[fiber],kw_LOCO_FILE)
hadmrFITS.copy_keys_root(loco_file,e2ds_fitsfilename[fiber],kw_root_drs_loc)
# Copy keywords of the flat field frame
kw_FLAT_FILE[1]=dic_db['FLAT_'+fiber][1]
hadmrFITS.write_newkey(e2ds_fitsfilename[fiber],kw_FLAT_FILE)
hadmrFITS.copy_keys_root(flat_file,e2ds_fitsfilename[fiber],kw_root_drs_flat)
