research/TCS/2012-09-07/deflection/analysis.py

   1 #!/usr/bin/python -u
   2
   3
   4 import os, sys
   5 lib_path = os.path.abspath('../../')
   6 sys.path.append(lib_path)
   7 import process
   8
   9 dirs = process.GetDataSets(".")
  10
  11
  12 #process.gnuplot("set title \"Vary Deflection Voltage\"")
  13 #process.gnuplot("set xlabel \"DAC (counts)\"")
  14 #process.gnuplot("set ylabel \"Normalised S(E)\"")
  15 #process.gnuplot("set key outside right")
  16
  17 max_curve = []
  18 for p in sorted(dirs.items(), key = lambda e : float(e[0].split("=")[1])):
  19         if (len(p[1]) > 0):
  20                 maxval = max(p[1], key = lambda e : e[1])
  21                 max_curve.append([float(p[0].split("=")[1]), maxval[1], maxval[0]])
  22 #               process.gnuplot.replot(process.Gnuplot.Data(p[1], title=str(p[0])+"", with_="lp"))
  23
  24
  25 process.gnuplot.replot(process.Gnuplot.Data(max_curve, with_="lp", using="1:2", title="Max"))
  26 process.gnuplot.replot(process.Gnuplot.Data(max_curve, with_="lp", using="1:3", title="Max position"))
  27 sys.stdin.readline()