Commit before breaking everything

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

diff --git a/research/analysis and stuff/tcs scripts/au_deposition.py~ b/research/analysis and stuff/tcs scripts/au_deposition.py~
new file mode 100755 (executable)
index 0000000..826df18
--- /dev/null
+++ b/research/analysis and stuff/tcs scripts/au_deposition.py~
@@ -0,0 +1,85 @@
+#!/usr/bin/python
+
+#Deposition for the FIRST sample
+
+from process import *
+import odict
+import time
+
+tcsdir = "../../research/TCS/"
+imagedir = "au_deposition/"
+
+smooth = [4,7]
+
+os.system("mkdir -p " + imagedir.strip("/"))
+
+toplot = odict.odict([
+       ("Si (initial)" , [tcsdir + "2012-08-09/074447.dat",1]),
+       ("+Au (3min)", [tcsdir + "2012-08-09/121616.dat", 1]),
+       ("+Au (8min)", [tcsdir + "2012-08-09/161028.dat", 1]),
+       ("+Au (8min) a 7hrs", [tcsdir + "2012-08-09/233758.dat", 1]),
+       ("+Au (8min) a 15hrs", [tcsdir + "2012-08-10/065228.dat", 1]),
+       ("+Au (18min)", [tcsdir + "2012-08-10/104218.dat", 1]),
+       ("+Au (18min) a 4hrs", [tcsdir + "2012-08-10/142746.dat", 1]),
+       ("+Au (18min) a 8hrs", [tcsdir + "2012-08-10/181247.dat", 1]),
+       ("+Au (18min) a 12hrs", [tcsdir + "2012-08-10/215747.dat", 1]),
+       ("+Au (18min) a 16hrs", [tcsdir + "2012-08-11/014200.dat", 1]),
+       ("+Au (18min) a 20hrs", [tcsdir + "2012-08-11/052314.dat", 1]),
+       ("+Au (18min) a 24hrs", [tcsdir + "2012-08-11/090345.dat", 1]),
+       ("+Au (18min) a 28hrs", [tcsdir + "2012-08-11/124357.dat", 1]),
+       ("+Au (18min) a 32hrs", [tcsdir + "2012-08-11/162419.dat", 1]),
+       ("+Au (18min) a 36hrs", [tcsdir + "2012-08-11/234511.dat", 1]),
+       ("+Au (18min) a 40hrs", [tcsdir + "2012-08-12/032517.dat", 1]),
+       ("+Au (18min) a 44hrs", [tcsdir + "2012-08-12/070514.dat", 1]),
+       ("+Au (18min) a 48hrs", [tcsdir + "2012-08-12/104457.dat", 1]),
+       ("+Au (18min) a 52hrs", [tcsdir + "2012-08-12/142446.dat", 1]),
+       ("+Au (18min) a 56hrs", [tcsdir + "2012-08-12/180443.dat", 1]),
+       ("+Au (18min) a 60hrs", [tcsdir + "2012-08-12/214440.dat", 1]),
+       ("+Au (18min) a 96hrs", [tcsdir + "2012-08-14/094851.dat", 1]),
+       ("+Au (33min)", [tcsdir + "2012-08-16/124728.dat", 1]),
+       #("+Au (33min) a 4hrs", [tcsdir + "2012-08-16/163328.dat", 1]), 
+       #("+Au (33min) a 5hrs", [tcsdir + "2012-08-16/170616.dat", 1]),
+       #("+Au (33min) a 18hrs", [tcsdir + "2012-08-17/062341.dat", 1]),
+       ("+Au (50min)", [tcsdir + "2012-08-17/074929.dat", 1]),
+       ("+Au (50min) a 1hr", [tcsdir + "2012-08-17/084544.dat", 1]),
+       ("+Au (83min) a 57hr", [tcsdir + "2012-08-19/172754.dat", 1]),  
+       ("+Au (100min)", [tcsdir + "2012-08-20/au_on_si(50min)/103405.dat", 1]),
+       ("+Au (100min) a 1min", [tcsdir + "2012-08-20/103534.dat", 1]),
+       ("+Au (100min) a 1hr", [tcsdir + "2012-08-20/113127.dat", 1]),
+       ("+Au (100min) a 2hrs", [tcsdir + "2012-08-20/124236.dat", 1]),
+       #("+Au (100min) a 4hrs", [tcsdir + "2012-08-20/143917.dat", 1]),
+       ("+Au (125min)", [tcsdir + "2012-08-20/153656.dat", 1]),
+       ("+Au (125min) a 20hrs", [tcsdir + "2012-08-21/191221.dat", 1]),
+       ("+Au (125min) a 7days", [tcsdir + "2012-08-27/124551.dat", 1]),
+       ("+Au (145min)", [tcsdir + "2012-08-27/141116.dat", 1]) 
+])
+
+# I(E) (normalised) data
+gnuplot.reset()
+gnuplot("set term postscript colour")
+gnuplot("set output \""+imagedir+"normalised_data.eps\"")
+gnuplot("set title \"Deposition of Au on Si I(E) curves\"")
+gnuplot("set key top left")
+
+for f in toplot.items():
+       d = ShowData(f[1][0], title=f[0], plot=gnuplot.plot, calibrate=True, normalise=True)
+
+
+
+time.sleep(1)
+
+# TCS (normalised) data
+gnuplot.reset()
+gnuplot("set term postscript colour")
+gnuplot("set output \""+imagedir+"tcs_normalised.eps\"")
+gnuplot("set title \"Deposition of Au on Si S(E) curves\"")
+gnuplot("set key top right")
+
+for f in toplot.items():
+       d = ShowTCS(f[1][0], title=f[0], plot=gnuplot.plot, calibrate=True, normalise=True, with_="lp", step=f[1][1], smooth=smooth)
+
+
+time.sleep(1)
+
+#os.system("cd "+imagedir+"; for i in *.eps; do epstopdf --autorotate=All $i; done; rm *.eps")
+