Commit before breaking everything

[matches/honours.git] / research / analysis and stuff / tcs scripts / au_thickness.py

diff --git a/research/analysis and stuff/tcs scripts/au_thickness.py b/research/analysis and stuff/tcs scripts/au_thickness.py
new file mode 100755 (executable)
index 0000000..3245151
--- /dev/null
+++ b/research/analysis and stuff/tcs scripts/au_thickness.py
@@ -0,0 +1,71 @@
+#!/usr/bin/python
+
+from process import *
+import odict
+import time
+
+tcsdir = "../../research/TCS/"
+imagedir = "au_thickness/"
+
+os.system("mkdir -p " + imagedir.strip("/"))
+
+allplots = [ ("set1", odict.odict([
+       #("Si (initial)" , [tcsdir+"2012-09-27/Si/095131.dat", 1])
+       ("Si", [tcsdir+"2012-09-27/Si (focused)/140416.dat", 5]) # For some reason I altered the focus of the gun, not a good idea
+       ,("Au (30min 4.0A)", [tcsdir+"2012-10-01/average.dat", 1])
+       ,("Au (100min 4.0A)", [tcsdir+"2012-10-03/Au (100min) on Si/average.dat", 1])
+       ,("Au (20min 4.6A + 120min 4.0A)", [tcsdir+"2012-10-03/Au (140min) on Si/average.dat", 1])
+       ,("Au (80min 4.6A + 120min 4.0A)", [tcsdir+"2012-10-03/Au (200min) on Si/average.dat", 1])
+       ,("Au (30min 4.8A + 80min 4.6A + 120min 4.0A)", [tcsdir+"2012-10-04/Au (230min) on Si/average.dat", 1])
+       #,("Au (30min 4.8A + 80min 4.6A + 120min 4.0A) - REFOCUS", [tcsdir+"2012-10-04/Au (230min) on Si - Refocused/average.dat", 1])
+])), ("set2", odict.odict([
+       ("Au (30min 4.8A + 80min 4.6A + 120min 4.0A)", [tcsdir+"2012-10-09/Au (230min) on Si/average.dat", 1])
+       ,("Au (30min 5.2A + 30min 4.8A + 80min 4.6A + 120min 4.0A)", [tcsdir+"2012-10-09/Au (260min) on Si/average.dat", 1])]))
+
+]
+
+for p in allplots:
+
+       toplot = p[1]
+       # I(E) (normalised) data
+       gnuplot.reset()
+
+       for f in toplot.items():
+               d = ShowData(f[1][0], title=f[0], plot=gnuplot.replot, calibrate=True, normalise=True)
+
+       gnuplot("set term postscript colour")
+       gnuplot("set output \""+imagedir+str(p[0])+"_"+"normalised_data.eps\"")
+       title="Time Evolution of Au on Si I(E) curves"
+       if (p[0] == "set1"):
+               title += "\\nFocus not altered between curves; different focus to set2"
+       else:
+               title += "\\nFocus not altered between curves; different focus to set1"
+       gnuplot("set title \""+str(title)+"\"")
+       gnuplot("set key top left")
+       gnuplot.replot()
+
+       time.sleep(1)
+
+       # TCS (normalised) data
+       gnuplot.reset()
+       gnuplot("set term wxt")
+
+       for f in toplot.items():
+               d = ShowTCS(f[1][0], title=f[0], plot=gnuplot.replot, calibrate=True, normalise=True, with_="lp", step=f[1][1])
+
+       gnuplot("set term postscript colour")
+       gnuplot("set output \""+imagedir+str(p[0])+"_"+"tcs_normalised.eps\"")
+       title="Time Evolution of Au on Si S(E) curves"
+       if (p[0] == "set1"):
+               title += "\\nFocus not altered between curves; different focus to set2"
+       else:
+               title += "\\nFocus not altered between curves; different focus to set1"
+       gnuplot("set title \""+str(title)+"\"")
+       gnuplot("set key top right")
+       gnuplot.replot()
+
+       time.sleep(1)
+
+
+os.system("cd "+imagedir+"; for i in *.eps; do epstopdf --autorotate=All $i; done; rm *.eps")
+