  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).

