# (c) Max Kresch, licensed under BSD, 2008.
#     See LICENSE for further details. 

The fonts used by this program are released under the GPL.  You may download
them thusly:
     $ cd fonts
     $ wget http://ftp.de.debian.org/debian/pool/main/t/ttf-thryomanes/ttf-thryomanes_1.2.4.tar.gz
and untar and unzip them into the directory "./fonts/", and create a symbolic 
link in that directory:
     $ tar xvzf ttf-thryomanes_1.2.4.tar.gz
     $ ln -sf ttf-thryomanes-1.2 ttf-thryomanes 
This should get the code to make prettier plots and complain less about missing 
fonts.

  The purpose of this code is to take S(Q,E) as determined on a time-of-flight 
chopper spectrometer, to correct for multiphonon and multiple-scattering, and 
to return a phonon DOS.  The code is run by editing "user.py" appropriately, 
and then typing at the command prompt ($):
     $ ./getDOS.py
Presumably, the input is a cPickle file with 4 numpy arrays in a tuple:
     ( Q,E,S(Q,E), error in S(Q,E) )
with S(Q,E).shape == ( E.shape, Q.shape).

  The python files are reasonably well documented -- in particular, the 
documentation in "user.py" describes the parameters that the user needs to set.
Also, in an interpreter, using "help(function)" on the various functions (or 
classes) can be quite useful.  The OUTPUT is sort of described in the file 
"Outputs.txt".

  There are some additional snipets of code as well.  If "dos2sqe.py" is edited
appropriately, it should open the 2 column ascii DOS specified by the variable
"dosFile".  When run:
     $ ./dos2sqe.py
it will produce files in "./tmp": "sqe.pkl", "se.?" where ? is a number 
describing the number of phonons involved in the scattering tabulated in the 
file.  Given the "./tmp/sqe.pkl" file:
     $ ./dos2sqePlot.py
produces an eps figure "fig.eps" that shows the intenstiy S(Q,E).
     $ ./dos2sqePlotSe.py
pops open a gnuplot window which shows S(E).

