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

   1 #!/usr/bin/python
   2
   3 from process import *
   4 import odict
   5 import time
   6
   7 tcsdir = "../../research/TCS/"
   8 imagedir = "au_thickness/"
   9
  10 os.system("mkdir -p " + imagedir.strip("/"))
  11
  12 allplots = [ ("set1", odict.odict([
  13         #("Si (initial)" , [tcsdir+"2012-09-27/Si/095131.dat", 1])
  14         ("Si", [tcsdir+"2012-09-27/Si (focused)/140416.dat", 5]) # For some reason I altered the focus of the gun, not a good idea
  15         ,("Au (30min 4.0A)", [tcsdir+"2012-10-01/average.dat", 1])
  16         ,("Au (100min 4.0A)", [tcsdir+"2012-10-03/Au (100min) on Si/average.dat", 1])
  17         ,("Au (20min 4.6A + 120min 4.0A)", [tcsdir+"2012-10-03/Au (140min) on Si/average.dat", 1])
  18         ,("Au (80min 4.6A + 120min 4.0A)", [tcsdir+"2012-10-03/Au (200min) on Si/average.dat", 1])
  19         ,("Au (30min 4.8A + 80min 4.6A + 120min 4.0A)", [tcsdir+"2012-10-04/Au (230min) on Si/average.dat", 1])
  20         #,("Au (30min 4.8A + 80min 4.6A + 120min 4.0A) - REFOCUS", [tcsdir+"2012-10-04/Au (230min) on Si - Refocused/average.dat", 1])
  21 ])), ("set2", odict.odict([
  22         ("Au (30min 4.8A + 80min 4.6A + 120min 4.0A)", [tcsdir+"2012-10-09/Au (230min) on Si/average.dat", 1])
  23         ,("Au (30min 5.2A + 30min 4.8A + 80min 4.6A + 120min 4.0A)", [tcsdir+"2012-10-09/Au (260min) on Si/average.dat", 1])]))
  24
  25 ]
  26
  27 for p in allplots:
  28
  29         toplot = p[1]
  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+str(p[0])+"_"+"normalised_data.eps\"")
  38         title="Time Evolution of Au on Si I(E) curves"
  39         if (p[0] == "set1"):
  40                 title += "\\nFocus not altered between curves; different focus to set2"
  41         else:
  42                 title += "\\nFocus not altered between curves; different focus to set1"
  43         gnuplot("set title \""+str(title)+"\"")
  44         gnuplot("set key top left")
  45         gnuplot.replot()
  46
  47         time.sleep(1)
  48
  49         # TCS (normalised) data
  50         gnuplot.reset()
  51         gnuplot("set term wxt")
  52
  53         for f in toplot.items():
  54                 d = ShowTCS(f[1][0], title=f[0], plot=gnuplot.replot, calibrate=True, normalise=True, with_="lp", step=f[1][1])
  55
  56         gnuplot("set term postscript colour")
  57         gnuplot("set output \""+imagedir+str(p[0])+"_"+"tcs_normalised.eps\"")
  58         title="Time Evolution of Au on Si S(E) curves"
  59         if (p[0] == "set1"):
  60                 title += "\\nFocus not altered between curves; different focus to set2"
  61         else:
  62                 title += "\\nFocus not altered between curves; different focus to set1"
  63         gnuplot("set title \""+str(title)+"\"")
  64         gnuplot("set key top right")
  65         gnuplot.replot()
  66
  67         time.sleep(1)
  68
  69
  70 os.system("cd "+imagedir+"; for i in *.eps; do epstopdf --autorotate=All $i; done; rm *.eps")
  71