コード例 #1
0
ファイル: obj_POL_harps.py プロジェクト: ivh/harps-pol
def fix_headers(filename,what):
    hadmrFITS.copy_key(header_base,filename)
    
    hadmrFITS.delete_key(filename,kw_pol_ret25) 
    hadmrFITS.delete_key(filename,kw_pol_ret50)

    # which emodulation?
    if 'diff' in filename:
        demod=kw_pol_demod_diff
    else:
        demod=kw_pol_demod_ratio
    hadmrFITS.write_newkey(filename,[kw_pol_demod, demod, kw_pol_demod_com])

    stokes=[kw_pol_stokes, kw_pol_stokes_v, kw_pol_stokes_com]
    if what=='V':
        hadmrFITS.write_newkey(filename,stokes)
    elif what=='Q':
        stokes[1]=kw_pol_stokes_q
        hadmrFITS.write_newkey(filename,stokes)
    elif what=='U':
        stokes[1]=kw_pol_stokes_u
        hadmrFITS.write_newkey(filename,stokes)
    elif what=='N':
        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)
コード例 #2
0
for i in indexlist: qq[i]=1
ccf_noise=(sum(compress(qq,CCF_slope**2)/compress(qq,CCF_noise_tot[1:-1]**2)))**(-0.5)

WLOG('info',log_opt,'Estimated RV accuracy on stellar spectrum: %.2f[m/s]' %(meanpondref))
WLOG('info',log_opt,'Estimated RV accuracy on CCF             : %.2f[m/s]' %(ccf_noise*1000.))

#
# archive in fits
#
WLOG('',log_opt,'Archiving CCF on file: '+cor_fitsfilename[fiber])
hadmrFITS.Save_ccf(RV_CCF,CCF,average_CCF,cor_fitsfilename[fiber],ccf_mask,\
sum(CCF_max)/sum(pix_passed_all),sum(tot_line[ccf_order_range]),abs(100*CCF_res[0]),CCF_res[1],CCF_res[2]*2.3548,RV,_kw_CCF)

# Add CCF photon noise keyword
kw_CCF_NOISE[1]=ccf_noise
hadmrFITS.write_newkey(cor_fitsfilename[fiber],kw_CCF_NOISE)

#
# copy e2ds file kw into CCF
#
hadmrFITS.copy_key(I_e2ds_fitsfilename,cor_fitsfilename[fiber])

# Update DVRMS keyword
kw_DVRMS[1]=meanpondref
hadmrFITS.update_key(cor_fitsfilename[fiber],kw_DVRMS)

# Add flux correction keywords
kw_FLUX_CORR_MIN[1]=corr.min()
kw_FLUX_CORR_MAX[1]=corr.max()
hadmrFITS.write_newkey(cor_fitsfilename[fiber],kw_FLUX_CORR_MIN)
hadmrFITS.write_newkey(cor_fitsfilename[fiber],kw_FLUX_CORR_MAX)
コード例 #3
0
ファイル: obj_sub_POL_harps.py プロジェクト: ivh/harps-pol
  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)
# Copy keywords of the blaze frame 
  kw_BLAZE_FILE[1]=dic_db['BLAZE_'+fiber][1]
  hadmrFITS.write_newkey(e2ds_fitsfilename[fiber],kw_BLAZE_FILE)