X-Git-Url: https://git.ucc.asn.au/?a=blobdiff_plain;f=research%2FTCS%2F2012-09-06%2Ffocus%2Fanalysis.py;fp=research%2FTCS%2F2012-09-06%2Ffocus%2Fanalysis.py;h=e900b7fb6fa2a709e81c2063261733043f56b09c;hb=ced554da0fb721d190eed99618e2cd82130f0797;hp=0000000000000000000000000000000000000000;hpb=19cb61f75ddda901224f7ce2bae0d3934f67d9ee;p=matches%2Fhonours.git

diff --git a/research/TCS/2012-09-06/focus/analysis.py b/research/TCS/2012-09-06/focus/analysis.py
new file mode 100755
index 00000000..e900b7fb
--- /dev/null
+++ b/research/TCS/2012-09-06/focus/analysis.py
@@ -0,0 +1,28 @@
+#!/usr/bin/python -u
+
+
+import os, sys
+lib_path = os.path.abspath('../../')
+sys.path.append(lib_path)
+import process
+
+dirs = {}
+
+for f in os.listdir("."):
+	if os.path.isdir(f) == True:
+		data_set = []
+		for datafile in os.listdir("./"+str(f)):
+			data_set.append(process.MaxNormalise(map(lambda e : [e[1], e[2]], process.GetData("./"+str(f)+"/"+str(datafile)))))
+
+		avg = process.Average(data_set)
+		dirs.update({f : avg})
+
+process.gnuplot("set title \"Vary Focus Voltage\"")
+process.gnuplot("set xlabel \"DAC (counts)\"")
+process.gnuplot("set ylabel \"Normalised S(E)\"")
+
+for p in dirs.items():
+	if (len(p[1]) > 0):
+		process.gnuplot.replot(process.Gnuplot.Data(p[1], title=str(p[0])+"", with_="lp"))
+
+sys.stdin.readline()