Density of States Analysis Example

  • Given sample and empty-can data, compute phonon DOS
  • To use this notebook, first click jupyter menu File->Make a copy
  • Click the title of the copied jupyter notebook and change it to a new title
  • Start executing cells

Preparation


In [ ]:
# where am I now?
!pwd

In [ ]:
# create a new working directory and change into it
workdir = '~/reduction/ARCS/getdos-demo'
!mkdir -p {workdir}
%cd {workdir}

In [ ]:
# Data to reduce. Change the IPTS number and run numbers to suit your need
samplenxs = "/SNS/ARCS/IPTS-17327/data/ARCS_83914_event.nxs"
mtnxs = "/SNS/ARCS/IPTS-17327/data/ARCS_83921_event.nxs"

Run GetDOS


In [ ]:
# import tools
import os, numpy as np
from multiphonon.getdos import notebookUI
import histogram.hdf as hh, histogram as H
%matplotlib notebook

In [ ]:
# create the UI for the first time
notebookUI(samplenxs, mtnxs)

Check output


In [ ]:
ls work/

In [ ]:
dos = hh.load('work/final-dos.h5')
H.plot(dos)

Refine GetDOS


In [ ]:
# if you need to run getdos again with slightly modified options
# you can start from the previous settings
notebookUI(samplenxs, mtnxs, load_options_path="./work/getdos-opts.yaml")