research/TCS/2012-09-07/venault/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(".", process.Derivative)
  10
  11
  12 process.gnuplot("set title \"Vary Venault 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
  19 include_range=[0.0, 100.0]
  20
  21 for p in sorted(dirs.items(), key = lambda e : float(e[0].split("=")[1])):
  22         value = float(p[0].split("=")[1])
  23         if (value >= include_range[0] and value <= include_range[1]):
  24                 if (len(p[1]) > 0):
  25                         maxval = max(p[1], key = lambda e : e[1])
  26                         max_curve.append([float(p[0].split("=")[1]), maxval[1], maxval[0]])
  27                         process.gnuplot.replot(process.Gnuplot.Data(p[1], title=str(p[0])+"", with_="lp"))
  28
  29
  30 #process.gnuplot.replot(process.Gnuplot.Data(max_curve, with_="lp", using="1:2", title="Max"))
  31 #process.gnuplot.replot(process.Gnuplot.Data(max_curve, with_="lp", using="1:3", title="Max position"))
  32 sys.stdin.readline()