X-Git-Url: https://git.ucc.asn.au/?a=blobdiff_plain;ds=sidebyside;f=research%2Fanalysis%20and%20stuff%2Ftcs%20scripts%2Fsi_substrate.py;fp=research%2Fanalysis%20and%20stuff%2Ftcs%20scripts%2Fsi_substrate.py;h=bfa42e9c7ed89e4d1769afc60b6795a34f362be6;hb=70a96cca12cb006506461d26cd112bab179fe74c;hp=0000000000000000000000000000000000000000;hpb=8caf60af39689a3546074f0c68d14c3a2e28191e;p=matches%2Fhonours.git

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
index 00000000..bfa42e9c
--- /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")
+