research/analysis and stuff/tcs scripts/au_deposition2.py

   1 #!/usr/bin/python
   2
   3 # Deposition for the SECOND sample
   4
   5 from process import *
   6 import odict
   7 import time
   8
   9 tcsdir = "../../research/TCS/"
  10 imagedir = "au_deposition2/"
  11
  12 smooth = [4,3]
  13
  14 os.system("mkdir -p " + imagedir.strip("/"))
  15
  16 toplot = odict.odict([
  17         ("Si (initial)" , [tcsdir+"2012-09-27/Si/095131.dat", 1]),
  18         ("Si (refocused)", [tcsdir+"2012-09-27/Si (focused)/140416.dat", 5]), # For some reason I altered the focus of the gun, not a good idea
  19         ("Au (30min 4.0A e-8mbar)", [tcsdir+"2012-09-27/Au on Si/151444.dat", 5]),
  20         #("Au +5min", tcsdir+"2012-09-27/151951.dat"),
  21         #("Au +1hour", tcsdir+"2012-09-27/161548.dat"),
  22         #("Au +2hour", tcsdir+"2012-09-27/171618.dat"),
  23         ("Au +8hour", [tcsdir+"2012-09-27/Au on Si/235714.dat", 5]),
  24         ("Au (avg over +8hour to +32hour)", [tcsdir+"2012-09-28/average.dat", 1]),
  25         ("Au (avg over +32hour to +56hour)", [tcsdir+"2012-09-29/average.dat", 1]),
  26         ("Au (avg over +56hour to +80hour)", [tcsdir+"2012-09-30/average.dat", 1]),
  27         ("Au (avg over +80hour to +104hour)", [tcsdir+"2012-10-01/average.dat", 1])
  28 ])
  29
  30 # I(E) (normalised) data
  31 gnuplot.reset()
  32
  33 for f in toplot.items():
  34         d = ShowData(f[1][0], title=f[0], plot=gnuplot.replot, calibrate=True, normalise=True)
  35
  36 gnuplot("set term postscript colour")
  37 gnuplot("set output \""+imagedir+"normalised_data.eps\"")
  38 gnuplot("set title \"Deposition of Au on Si I(E) curves\"")
  39 gnuplot("set key top left")
  40 gnuplot.replot()
  41
  42 time.sleep(1)
  43
  44 # TCS (normalised) data
  45 gnuplot.reset()
  46 gnuplot("set term wxt")
  47 gnuplot("set term postscript colour")
  48 gnuplot("set output \""+imagedir+"tcs_normalised.eps\"")
  49 gnuplot("set title \"Deposition of Au on Si S(E) curves\"")
  50 gnuplot("set key top right")
  51
  52 for f in toplot.items():
  53         d = ShowTCS(f[1][0], title=f[0], plot=gnuplot.plot, calibrate=True, normalise=True, with_="lp", step=1, smooth=smooth)
  54 time.sleep(1)
  55
  56 #os.system("cd "+imagedir+"; for i in *.eps; do epstopdf --autorotate=All $i; done; rm *.eps")
  57