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

   1 #!/usr/bin/python
   2
   3 #Deposition for the FIRST sample
   4
   5 from process import *
   6 import odict
   7 import time
   8
   9 tcsdir = "../../research/TCS/"
  10 imagedir = "au_deposition/"
  11
  12 smooth = [4,5]
  13
  14 os.system("mkdir -p " + imagedir.strip("/"))
  15
  16 toplot = odict.odict([
  17         ("Si (initial)" , [tcsdir + "2012-08-09/074447.dat",1]),
  18         ("+Au (3min)", [tcsdir + "2012-08-09/121616.dat", 1]),
  19         ("+Au (8min)", [tcsdir + "2012-08-09/161028.dat", 1]),
  20         ("+Au (8min) a 7hrs", [tcsdir + "2012-08-09/233758.dat", 1]),
  21         ("+Au (8min) a 15hrs", [tcsdir + "2012-08-10/065228.dat", 1]),
  22         ("+Au (18min)", [tcsdir + "2012-08-10/104218.dat", 1]),
  23         ("+Au (18min) a 4hrs", [tcsdir + "2012-08-10/142746.dat", 1]),
  24         ("+Au (18min) a 8hrs", [tcsdir + "2012-08-10/181247.dat", 1]),
  25         ("+Au (18min) a 12hrs", [tcsdir + "2012-08-10/215747.dat", 1]),
  26         ("+Au (18min) a 16hrs", [tcsdir + "2012-08-11/014200.dat", 1]),
  27         ("+Au (18min) a 20hrs", [tcsdir + "2012-08-11/052314.dat", 1]),
  28         ("+Au (18min) a 24hrs", [tcsdir + "2012-08-11/090345.dat", 1]),
  29         ("+Au (18min) a 28hrs", [tcsdir + "2012-08-11/124357.dat", 1]),
  30         ("+Au (18min) a 32hrs", [tcsdir + "2012-08-11/162419.dat", 1]),
  31         ("+Au (18min) a 36hrs", [tcsdir + "2012-08-11/234511.dat", 1]),
  32         ("+Au (18min) a 40hrs", [tcsdir + "2012-08-12/032517.dat", 1]),
  33         ("+Au (18min) a 44hrs", [tcsdir + "2012-08-12/070514.dat", 1]),
  34         ("+Au (18min) a 48hrs", [tcsdir + "2012-08-12/104457.dat", 1]),
  35         ("+Au (18min) a 52hrs", [tcsdir + "2012-08-12/142446.dat", 1]),
  36         ("+Au (18min) a 56hrs", [tcsdir + "2012-08-12/180443.dat", 1]),
  37         ("+Au (18min) a 60hrs", [tcsdir + "2012-08-12/214440.dat", 1]),
  38         ("+Au (18min) a 96hrs", [tcsdir + "2012-08-14/094851.dat", 1]),
  39         ("+Au (33min)", [tcsdir + "2012-08-16/124728.dat", 1]),
  40         #("+Au (33min) a 4hrs", [tcsdir + "2012-08-16/163328.dat", 1]),
  41         #("+Au (33min) a 5hrs", [tcsdir + "2012-08-16/170616.dat", 1]),
  42         #("+Au (33min) a 18hrs", [tcsdir + "2012-08-17/062341.dat", 1]),
  43         ("+Au (50min)", [tcsdir + "2012-08-17/074929.dat", 1]),
  44         ("+Au (50min) a 1hr", [tcsdir + "2012-08-17/084544.dat", 1]),
  45         ("+Au (83min) a 57hr", [tcsdir + "2012-08-19/172754.dat", 1]),
  46         ("+Au (100min)", [tcsdir + "2012-08-20/au_on_si(50min)/103405.dat", 1]),
  47         ("+Au (100min) a 1min", [tcsdir + "2012-08-20/103534.dat", 1]),
  48         ("+Au (100min) a 1hr", [tcsdir + "2012-08-20/113127.dat", 1]),
  49         ("+Au (100min) a 2hrs", [tcsdir + "2012-08-20/124236.dat", 1]),
  50         #("+Au (100min) a 4hrs", [tcsdir + "2012-08-20/143917.dat", 1]),
  51         ("+Au (125min)", [tcsdir + "2012-08-20/153656.dat", 1]),
  52         ("+Au (125min) a 20hrs", [tcsdir + "2012-08-21/191221.dat", 1]),
  53         ("+Au (125min) a 7days", [tcsdir + "2012-08-27/124551.dat", 1]),
  54         ("+Au (145min)", [tcsdir + "2012-08-27/141116.dat", 1])
  55 ])
  56
  57 # I(E) (normalised) data
  58 gnuplot.reset()
  59 gnuplot("set term postscript colour")
  60 gnuplot("set output \""+imagedir+"normalised_data.eps\"")
  61 gnuplot("set title \"Deposition of Au on Si I(E) curves\"")
  62 gnuplot("set key top left")
  63
  64 for f in toplot.items():
  65         d = ShowData(f[1][0], title=f[0], plot=gnuplot.plot, calibrate=True, normalise=True)
  66
  67
  68
  69 time.sleep(1)
  70
  71 # TCS (normalised) data
  72 gnuplot.reset()
  73 gnuplot("set term postscript colour")
  74 gnuplot("set output \""+imagedir+"tcs_normalised.eps\"")
  75 gnuplot("set title \"Deposition of Au on Si S(E) curves\"")
  76 gnuplot("set key top right")
  77
  78 for f in toplot.items():
  79         d = ShowTCS(f[1][0], title=f[0], plot=gnuplot.plot, calibrate=True, normalise=True, with_="lp", step=f[1][1], smooth=smooth)
  80
  81
  82 time.sleep(1)
  83
  84 #os.system("cd "+imagedir+"; for i in *.eps; do epstopdf --autorotate=All $i; done; rm *.eps")
  85