Commit before breaking everything

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

diff --git a/research/analysis and stuff/tcs scripts/si_substrate.py~ b/research/analysis and stuff/tcs scripts/si_substrate.py~
new file mode 100755 (executable)
index 0000000..b88e055
--- /dev/null
+++ b/research/analysis and stuff/tcs scripts/si_substrate.py~
@@ -0,0 +1,79 @@
+#!/usr/bin/python
+
+from process import *
+import odict
+import time
+
+tcsdir = "../../research/TCS/"
+
+
+smooth=[2,5]
+#os.system("rm -rf si_substrate")
+os.system("mkdir -p si_substrate")
+
+toplot = odict.odict([
+       ("Si Substrate" , tcsdir+"2012-09-20/171922.dat"),
+       ("+5min" , tcsdir+"2012-09-20/172442.dat"),
+       ("+1hour" , tcsdir+"2012-09-20/181935.dat"),
+       ("+2hour" , tcsdir+"2012-09-20/192017.dat"),
+       ("+7day" , tcsdir+"2012-09-27/Si/095131.dat")
+])
+
+# I(E) data
+gnuplot.reset()
+
+for f in toplot.items():
+       d = ShowData(f[1], title=f[0], plot=gnuplot.replot, calibrate=True)
+
+gnuplot("set term postscript colour")
+gnuplot("set output \"si_substrate/raw_data.eps\"")
+gnuplot("set title \"Evolution of Si Substrate I(E) curves\"")
+gnuplot("set key top left")
+gnuplot.replot()
+
+time.sleep(1)
+
+#sys.exit(0)
+
+# NORMALISED I(E) data
+gnuplot.reset()
+gnuplot("set term wxt")
+
+for f in toplot.items():
+       d = ShowData(f[1], title=f[0], plot=gnuplot.replot, calibrate=True, normalise=True)
+gnuplot("set term postscript colour")
+gnuplot("set output \"si_substrate/normalised_data.eps\"")
+gnuplot("set title \"Evolution of Si Substrate Normalised I(E) curves\"")
+gnuplot("set key top left")
+gnuplot.replot()
+time.sleep(1)
+
+# TCS data
+gnuplot.reset()
+gnuplot("set term wxt")
+gnuplot("set term postscript colour")
+gnuplot("set output \"si_substrate/tcs.eps\"")
+gnuplot("set title \"Evolution of Si Substrate TCS\"")
+gnuplot("set key top right")
+
+for f in toplot.items():
+       #gnuplot.reset()
+       d = ShowTCS(f[1], title=f[0], plot=gnuplot.replot, calibrate=True, smooth=smooth)
+
+
+time.sleep(1)
+# NORMALISED TCS data
+gnuplot.reset()
+gnuplot("set term wxt")
+gnuplot("set term postscript colour")
+gnuplot("set output \"si_substrate/tcs_normalised.eps\"")
+gnuplot("set title \"Evolution of Si Substrate Normalised TCS\"")
+gnuplot("set key top right")
+
+for f in toplot.items():
+       gnuplot.reset()
+       d = ShowTCS(f[1], title=f[0], plot=gnuplot.replot, normalise=True, calibrate=True, smooth=smooth)
+
+time.sleep(1)
+#os.system("cd si_substrate; for i in *.eps; do epstopdf --autorotate=All $i; done; rm *.eps")
+