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)