X-Git-Url: https://git.ucc.asn.au/?a=blobdiff_plain;f=research%2FTCS%2Fprocess.py;h=efd5dce86baa111ec8c8c212bc95cb32fa5cc630;hb=b64d0037c97e4d874f1b7e592d8831f0ac72dcde;hp=19834ff504649fcdddd418dccfff8fa666be70d9;hpb=cabb4dd339995baee61abba99bcdd16acc7e447f;p=matches%2Fhonours.git diff --git a/research/TCS/process.py b/research/TCS/process.py index 19834ff5..efd5dce8 100755 --- a/research/TCS/process.py +++ b/research/TCS/process.py @@ -30,6 +30,8 @@ def GetData(filename): def GetTCS(data): result = [] n = 0 + dI = 0 + dE = 0 for i in range(2, len(data)-1): dE = data[i+1][1] - data[i-1][1] if (dE != 0): @@ -42,16 +44,27 @@ def GetTCS(data): result.append([data[i][1], dI / (n * dE)]) return result +def Plot(*args): + gnuplot.plot(args) + +def FitTCS(data): + + def main(): - if (len(sys.argv) != 2): - sys.stderr.write(sys.argv[0] + " - Require 1 argument (filename)\n") + if (len(sys.argv) < 2): + sys.stderr.write(sys.argv[0] + " - Require arguments (filename)\n") return 1 - tcs = GetTCS(GetData(sys.argv[1])) - gnuplot.plot(tcs) - + tcs = [] + gnuplot("set style data lp") + for i in range(1, len(sys.argv)): + tcs.append(GetTCS(GetData(sys.argv[i]))) + if (len(tcs[i-1]) > 0): + gnuplot.replot(tcs[i-1]) + + print("Press enter to exit") sys.stdin.readline()