DOPPLER TOMOGRAPHY SOFTWARE IN IDL FOR THE RECONSTRUCTION 
OF THE SPECTRA OF SPECTROSCOPIC BINARY STARS

Douglas R. Gies
CHARA, Georgia State University
gies@chara.gsu.edu
04-26-2006

Introduction:
The separation of the spectra of spectroscopic binaries is 
an important tool in studies of the physical properties of 
the component stars.  There are several numerical techniques
available now that use the orbital Doppler shifts to enable
separation of the component spectra.  A brief review of the 
techniques and results is given in the attached PDF file, 
giesdt.pdf, from the article: 
Gies, D. R. 2004, "Doppler Tomography of Binary Stars," 
in Spectroscopically and Spatially Resolving the Components 
of Close Binary Stars (ASP Conf. Ser. Vol. 318), 
ed. R. W. Hilditch, H. Hensberge, and  K. Pavlovski 
(San Francisco: ASP), 61. 
This package contains IDL routines developed at GSU 
for spectral reconstruction by the Iterative Least-Squares
Technique of Doppler tomography. 

Contents:
autorect.pro                Simple rectification procedure
giesdt.pdf                  PDF of review article
ismout.fits                 Spectrum of ISM features
ratio1999m.fits             Ratio of sec/pri flux vs. lambda
readsbcm.pro                Reads in orbital velocities
readstk.pro                 Reads in composite spectra
ry_per.dat                  Wavelengths, times for RY Per
ry_per.fits                 RY Per red spectra 
sbcm.out.pri                Orbital elements for pri. 
sbcm.out.sec                Orbital elements for sec.
tom_demo.pro                Sample run for RY Per
tomography.pro              General tomography procedure
tomography_rlambda.pro      Tomography with variable flux ratio
tomsetup.pro                Generic set up for tomography

Example run for RY Persei:
The basic ILST algorithm is performed in the IDL procedure
tomography.pro.  The application of the procedure will 
depend on details of the binary: for example, the number 
of components, the presence of interstellar lines, the
temperature difference of the components, and whether or 
not the system is eclipsing.  The procedure tom_demo.pro 
gives an example of a reconstruction for the Algol-type 
system RY Per (see Barai et al. 2004, ApJ, 608, 989). 
To run this program in IDL, type the following:
IDL> .run tom_demo
IDL> tom_demo,wave,obsspec,recspec,pshift,ratio,simspec
This will run the tomography algorithm and make plots 
of the reconstructed primary (top) and secondary (bottom). 
See the comments in tom_demo.pro for explanations of 
the various steps. 
