dopcode/ctio_init.pro at master · mattgiguere/dopcode · GitHub

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
function ctio_init, obsnm, iss_nm, iss_bc, iss_obnm=iss_obnm,tag=tag

; fischer oct 2009
; dopenv=ctio_init('rqa10.7353','iod',0.0,iss_obnm='rqa10.7353')

; SETUP YOUR DIRECTORIES AND RENAME THE PATHS BELOW
  root_path='/home/fischer/dop2/'            ; the code lives here
  obs_dir = '/mir7/fitspec/'                         ; the fits files (obs) live here

 ; for later - the initial guesses are stored in the "ipcf" file
  ipcf_dir = root_path+'ipcf_files/' ; the ipcf directory
  ipcf_st = 'ipcf_ctio'+tag+'.dat'
  ipcf_file = ipcf_dir+ipcf_st

 ; the FTS spectrum is here
  fts_dir = root_path+'atlas/'
  fts_atlas = 'ftseso50.bin'    ; fts_ctio.dat gives poor fits??

 ; the barycentric correction file is here
  bary_file = '/mir7/bary/qbcvel.dat/'
  files_dir = '/mir7/files/'    ; store the output (vd) files here

 ; 1851 pixels per order at Lick, start pixel = 50, end pixel = 1809
 ; iodine starts in order 54, at 5020 A, pixel 1090 to 1809 (9 chunks)
 ; good iodine lines in orders 35 - 53 (22 chunks per order = 418 chunks)
 ; useable iodine in order 34 from pixel 50 to 1170 = 14 chunks
 ; total number of chunks: 441

 ; 2048 pixels per order at CTIO, start pixel=520, end pixel=1479
 ; iodine starts in order 5013 A, 12 chunks per order
 ; useable orders 5 - 23, 12 chunks/ord * 18 orders = 216 chunks

  st_pix = 280   ; at CTIO, trim the first 520 pixels because low SNR
  st_ord = 5     ; at CTIO, iodine starts in order 5
  n_pix = 80
  n_ord = 18
  nch_ord = 18               ; at CTIO, number of chunks per order
  n_chunks= n_ord*nch_ord    ; 710 for keck
  n_wcof = 2                 ; number of coefficients for fitting wavelength of one chunk

  ; I store CTIO stuff in /mir7, hence variable names for the PSF parameters
  ; NOTE: in pass 0 of the doppler code, the FWHM of the central gaussian is
  ; a free parameter.  Therefore, m7_psfsig will be reset and then fixed.
  ;           [ 1      2   3     4     5     6      7     8    9    10    11   12    13    14   15   16  17]
  m7_ampl =  [0.25,  0.0,0.001,0.003,0.006, 0.03, 0.11, 0.21, 0.21, 0.11,0.03,0.007,0.002,0.001,0.0, 0.0, 0.0]
;  m7_psfpix = [0.0, -4.0, -3.6, -3.0, -2.4, -1.8, -1.2, -0.8, -0.4, 0.4, 0.8,  1.2,  1.8,  2.4, 3.0, 3.6, 4.0]
;  m7_psfsig = [0.7, 0.0,  0.8,  0.7,  0.7,  0.7,  0.7,  0.6,  0.5, 0.5, 0.6,  0.7,  0.7,  0.7, 0.7, 0.8, 0.0]
  m7_psfpix = [0.0, -4.0, -3.6, -3.0, -2.4, -1.8, -1.2, -0.8, -0.4, 0.4, 0.8,  1.2,  1.8,  2.4, 3.0, 3.6, 4.0]
  m7_psfsig = [0.9, 0.0,  0.8,  0.7,  0.7,  0.7,  0.7,  0.6,  0.5, 0.5, 0.6,  0.7,  0.7,  0.7, 0.7, 0.8, 0.0]

  psfpix = m7_psfpix
  psfsig = m7_psfsig
  psf_ampl = m7_ampl
  osamp = 4

 ; del_ord and del_pix are used for smoothing the PSF in the second
 ; pass - this suggests that the PSF is constant over this spatial
 ; size on the CCD.  Here, I'm guessing it is smooth over 3
 ; orders (cross-disp) and 120 pixels (disp direction)
  del_ord = 3
  del_pix = 200
  ord_dist = 30  ;typical separation between orders in pixels

; for SBs it can be useful to cross correlate instead of using
; the barycentric shift
  xcorl_lambda = 4976.
  xcorl_pix1 = 900
  xcorl_npix = 500
  xcorl_ord = 3     ;[findgen(500)+900,3]
  xcorl_shiftrange=5     ; make this larger for SB's - otherwise used to refine the wav guess
  decker = '100m_fiber'  ; "B5" or "B1" for Keck, "64" or "80" for Lick
  file_ext = 0

;stop
  dopenv = dop_init(obsnm, iss_nm, iss_bc, iss_obnm = iss_obnm, bary_file = bary_file, del_ord = del_ord, $
                del_pix = del_pix, ord_dist = ord_dist, npix = n_pix, files_dir = files_dir, obs_dir=obs_dir, $
                fts_dir = fts_dir, fts_atlas = fts_atlas, ipcf_file = ipcf_file,$
                psfpix = psfpix, psfsig = psfsig, psf_ampl = psf_ampl, osamp=osamp, $
                st_pix=st_pix, st_ord=st_ord, n_ord=n_ord, nch_ord = nch_ord, n_chunks=n_chunks, $
                decker = decker, dpad = dpad, n_wcof = n_wcof, xcorl_lambda = xcorl_lambda, $
                xcorl_pix1=xcorl_pix1, xcorl_ord=xcorl_ord, xcorl_shiftrange=xcorl_shiftrange,xcorl_npix=xcorl_npix,$
                file_ext = file_ext)


  return,dopenv
end