Automatic commit. Wed Sep 5 00:00:03 WST 2012
[matches/honours.git] / research / TCS / interface.py
1 #!/usr/bin/python
2
3 # @file "interface.py"
4 # @author Sam Moore
5 # @purpose AVR Butterfly Datalogger control for Total Current Spectroscopy Measurements
6 #       The AVR responds to a limited number of RS232 commands to read ADC channels and set the DAC
7 #       This script can be used to control the DAC and monitor initial energy & sample current over RS232
8
9 import sys
10 import os
11 import time
12 import serial
13 import datetime
14
15 import odict
16
17
18 # TODO: Insert variables for calibration purposes here.
19 calibrate = {
20         "ADC_Counts" : [2**10,2**10,2**10,2**10,2**10,2**10,2**10], # Maxed counts on ADC channels (10 bits = 2**10 = 1024)
21         "DAC_Counts" : 2**12, # Maxed counts on DAC channel (12 bits = 2**12 = 4096)
22
23         # Calibration data for DAC and ADC
24         "DAC" : None, 
25         "ADC" : [None, None, None, None, None, None, None, None],
26
27         # Data used for rough calibration if above data is not present
28         "Vref" : 3.3, # The output of the voltage regulator (Vref = Vcc) relative to signal ground
29         "ADC_Rin" : [None, None, None, None, 15000, 1100, 1100, 1100],
30         "ADC_Rvar" : [None, None, None, None, 1000, 5000, 10000, 50000],
31         "DAC_Gain" : 1
32 }
33
34 # TODO: Adjust aqcuisition parameters here
35 aquire = { "DAC_Sweep" : "0.0 + 10.0*int(step)", # DAC Sweep value (t is in STEPS, not seconds!)
36         "ADC_Averages" : 200,
37         #"ADC_Vi" : 5, # ADC channel to read back Vi (set by DAC) through
38         #"ADC_Is" : 4, # ADC channel to read back Is through
39         #"ADC_Ie" : 4, # ADC channel to read back Ie through
40         "DAC_Settle" : 2.0, # Time in seconds to wait for DAC to stabilise
41         #"response_wait" : 0.2, # Time to wait in seconds between sending data and reading back
42         "start_date" : None,
43         "open_files" : []
44 }
45
46 #Setup the serial connection parameters
47 ser = serial.Serial(
48         port="/dev/ttyUSB1", # Modify as needed (note: in linux need to run `sudo chmod a+rw /dev/ttyUSBX' to set permissions)
49
50         # Do not change the values below here (unless AVR butterfly is reprogrammed to use different values)
51         baudrate=4800,
52         parity="N",
53         stopbits=1,
54         bytesize=8,
55         timeout=None,
56         xonxoff=0,
57         rtscts=0
58 )
59
60 parameters = odict.odict([
61         ("Accelerating Voltage" , None),
62         ("Focus Voltage" , None),
63         ("Deflection Voltage" , None),
64         ("Venault Voltage" , None),
65         ("Initial Voltage" , None),
66         ("Heating Current" , None),
67         ("Heating Voltage" , None),
68         ("Chamber Pressure" , None),
69         ("610B Zero" , None),
70         ("602 Zero" , None),
71         ("610B Scale" , None),
72         ("602 Scale" , None),
73         ("602 0.1 Battery" , None),
74         ("602 0.03 Battery" , None),
75         ("602 0.01 Battery" , None),
76         ("602 0.003 Battery" , None),
77         ("602 0.001 Battery" , None), 
78         ("ADC Regulator" , None),
79         ("Sample", None),
80         ("Sample Angle", None),
81         ("Title" , None),
82         ("Comment" , None),
83         ("Data" , None),
84         ("Parameters last checked", None)
85 ])
86
87 def getTime():
88         return str(datetime.datetime.now()).split(" ")[1].split(".")[0].replace(":","")
89
90 def getDate():
91         return str(datetime.datetime.now()).split(" ")[0]
92
93 # Used for when I press Control-C to stop things
94 def set_exit_handler(func):
95     if os.name == "nt":
96         try:
97             import win32api
98             win32api.SetConsoleCtrlHandler(func, True)
99         except ImportError:
100             version = ".".join(map(str, sys.version_info[:2]))
101             raise Exception("pywin32 not installed for Python " + version)
102     else:
103         import signal
104         signal.signal(signal.SIGTERM, func)
105         signal.signal(signal.SIGINT, func)
106
107 def killed_handler(signal, frame):
108         reason = ""
109         sys.stdout.write("\n# Reason for killing program? ")
110         reason = sys.stdin.readline().strip("\r\n ")
111         for out in aquire["open_files"]:
112                 sys.stdout.write("# Closing file " + str(out) + "\n")
113                 out.write("# Recieved KILL signal.\n# Reason: " + str(reason) + "\n")
114                 log_close(out)
115
116         
117
118 def cleanup():
119         for out in aquire["open_files"]:
120                 out.write("# Program exits.\n")
121                 log_close(out)
122
123 def log_open(a, b):
124         result = open(a, b,0)
125         if (b == "w"):
126                 result.write("# File opened at " + str(datetime.datetime.now()) + "\n")
127                 aquire["open_files"].append(result)
128         return result
129
130 def log_close(afile):
131         if (afile in aquire["open_files"]):
132                 afile.write("# File closed at " + str(datetime.datetime.now()) + "\n")
133                 aquire["open_files"].remove(afile)
134         
135         afile.close()
136
137
138 def init():
139         #import atexit
140         #atexit.register(cleanup)
141         set_exit_handler(killed_handler)
142         
143         aquire["start_date"] = getDate()
144
145
146
147         # Connect serial
148         ser.open()
149         ser.isOpen()
150 #       print("Waiting for \"# hello\" from device...")
151 #       while (ser.readline().strip("\r\n") != "# hello"):
152 #               pass
153         #while (ser.readline().strip("\r\n ") != "#"):
154         #       pass
155         time.sleep(1.0)
156
157         ser.write("a "+str(aquire["ADC_Averages"]) + "\r\n")
158         print(ser.readline().strip("\r\n"))
159         print(ser.readline().strip("\r\n"))
160         print(ser.readline().strip("\r\n"))
161
162         #print("Writing config information to config.dat...")
163         #output = log_open("config.dat", "w", 1)
164
165         #output.write("# Initialise " + str(datetime.datetime.now()) + "\n")
166
167         #for field in calibrate:
168         #       output.write("# calibrate["+str(field)+"] = "+str(calibrate[field]) + "\n")
169         #output.write("\n")
170         #for field in aquire:
171         #       output.write("# aquire["+str(field)+"] = "+str(aquire[field]) + "\n")
172
173         #output.write("# Ready " + str(datetime.datetime.now()) + "\n# EOF\n")
174         #output.close()
175
176 def main():
177
178         init()
179         
180         # I haven't ever used calibrated results, and yet this code is still here, why???
181         #if (loadCalibration_DAC() == False):
182         #       if (calibrateDAC() == False):
183         #               return -1
184         #if (loadCalibration_ADC(aquire["ADC_Is"]) == False):
185         #       if (calibrateADC_usingDAC(aquire["ADC_Is"], False) == False):
186         #               if (calibrateADC(aquire["ADC_Is"]) == False):
187         #                       return -1
188
189         #if (loadCalibration_ADC(aquire["ADC_Vi"]) == False):
190         #       if (calibrateADC_usingDAC(aquire["ADC_Vi"], True) == False):
191         #               if (calibrateADC(aquire["ADC_Vi"]) == False):
192         #                       return -1
193         
194
195         # Make directory for today, backup calibration files
196         os.system("mkdir -p " + getDate())
197         #os.system("cp *.dat " + getDate() +"/")
198
199         checkList()
200
201         
202         
203                 
204
205         # Experiment
206         # TODO: Modify data to record here
207         sweep = 1
208         while True:
209                 os.system("mkdir -p " + getDate())
210                 record_data([5, 0], getDate()+"/"+str(getTime())+".dat", None, 4000)
211                 sweep += 1
212         
213
214 def checkList():
215         try:
216                 input_file = log_open(getDate()+"/checklist", "r")
217         except:
218                 input_file = None
219
220         if (input_file != None):
221                 for line in input_file:
222                         k = line.split("=")
223                         item = None
224                         if (len(k) >= 2):
225                                 item = k[0].strip("# \r\n")
226                                 value = k[1].strip("# \r\n")
227
228                         if (item in parameters):
229                                 parameters[item] = value
230         
231                 print("Checklist found. Overwrite? [Y/n]")
232                 response = sys.stdin.readline().strip(" \r\n")
233                 if (response == "" or response == "y" or response == "Y"):
234                         input_file = log_open(getDate()+"/checklist.old", "w")
235                         for item in parameters:
236                                 input_file.write("# " + str(item) + " = " + str(parameters[item]) + "\n")
237                         input_file.write("\n")
238                         log_close(input_file)
239                         input_file = None
240         
241         if (input_file == None):
242                 for item in parameters:
243                         if item == "Parameters last checked":
244                                 continue
245                         sys.stdout.write("\""+str(item)+"\" = " + str(parameters[item]) + " New value?: ")
246                         response = sys.stdin.readline().strip("\r\n ")
247                         if (response != ""):
248                                 parameters[item] = response
249                         sys.stdout.write("\n")
250                 parameters["Parameters last checked"] = str(datetime.datetime.now())
251                         
252
253         checklist = log_open(getDate()+"/checklist", "w")
254         for item in parameters:
255                 checklist.write("# "+str(item) + " = " + str(parameters[item]) + "\n")
256                 #output_file.write("# "+str(item) + " = " + str(parameters[item]) + "\n")
257         log_close(checklist)
258         
259
260 def record_data(ADC_channels, output, pollTime = None, dac_max = None):
261         
262         if (output != None):
263                 output = [log_open(output, "w"), sys.stdout]
264         else:
265                 output = [sys.stdout]
266
267         for field in aquire:
268                 for out in output:
269                         out.write("# aquire["+str(field)+"] = "+str(aquire[field]) + "\n")
270         
271         for out in output:
272                 out.write("# Parameters:\n")
273
274         for field in parameters:
275                 for out in output:
276                         out.write("# "+str(field)+" = " + str(parameters[field]) + "\n")
277
278
279         start_time = time.time()
280         
281         for out in output:
282                 out.write("\n")
283                 out.write("# Experiment " + str(datetime.datetime.now()) + "\n")
284                 out.write("# Polling for " + str(pollTime) + "s.\n")
285                 out.write("\n")
286                 out.write("# Data:\n")
287                 out.write("# time\tDAC")
288                 for channel in ADC_channels:
289                         out.write("\tADC"+str(channel))
290                 out.write("\n")
291
292         step = 0
293         dacValue = int(eval(aquire["DAC_Sweep"]))
294         if (setDAC(dacValue) == False):
295                 setDAC(dacValue)
296         while (pollTime == None or time.time() < start_time + pollTime):
297                 if (aquire["DAC_Sweep"] != None):
298                         nextDacValue = int(eval(aquire["DAC_Sweep"]))
299                         if (nextDacValue != dacValue):
300                                 dacValue = nextDacValue
301                                 if (dacValue < 0):
302                                         break
303                                 setDAC(dacValue)
304                         step += 1
305                 
306
307                 if (dac_max != None and dacValue >= dac_max):
308                         break
309
310                 measure_start = time.time()
311
312                 raw_adc = []
313
314                 for channel in ADC_channels:
315                         read = readADC(channel)
316                         if read == False:
317                                 for out in output:              
318                                         print("# Abort data collection due to failed ADC read")                                 
319                                         if out != sys.stdout:
320                                                 log_close(out)
321                                         return False
322                         raw_adc.append((channel, read[0], read[1]))
323
324                 end_time = time.time()
325
326                 for out in output:
327                         out.write(str((measure_start + (end_time - measure_start)/2.0 - start_time)))
328                         out.write("\t"+str(dacValue))
329                         for adc in raw_adc:
330                                 out.write("\t" + str(adc[1]) + "\t" + str(adc[2]))
331                         out.write("\n") 
332         
333                 
334         for out in output:              
335                 if out != sys.stdout:
336                         log_close(out)
337         return True
338
339 def loadCalibration_ADC(channel):
340         try:
341                 input_file = log_open("calibrateADC"+str(channel)+".dat")
342         except:
343                 print("Couldn't find calibration file for ADC " + str(channel))
344                 return False
345         
346         calibrate["ADC"][channel] = []
347         for l in input_file:
348                 l = l.split("#")[0].strip("\r\n ")
349                 if (l == ""):
350                         continue
351                 else:
352                         split_line = l.split("\t")
353                         calibrate["ADC"][channel].append((float(split_line[0]), float(split_line[1])))
354         log_close(input_file)
355
356         if (len(calibrate["ADC"][channel]) <= 0):
357                 print("Empty calibration file for ADC " + str(channel))
358                 return False
359         return True
360
361 def loadCalibration_DAC():
362         try:
363                 input_file = log_open("calibrateDAC.dat")
364         except:
365                 print("Couldn't find calibration file for DAC")
366                 return False
367         
368         calibrate["DAC"] = []
369         for l in input_file:
370                 #print("Line is: "+str(l))
371                 l = l.split("#")[0].strip("\r\n ")
372                 if (l == ""):
373                         continue
374                 else:
375                         split_line = l.split("\t")
376                         if (len(split_line) >= 3):
377                                 calibrate["DAC"].append((int(split_line[0]), float(split_line[1]), float(split_line[2])))
378                         else:
379                                 calibrate["DAC"].append((int(split_line[0]), float(split_line[1])))
380
381
382         log_close(input_file)
383
384         if (len(calibrate["DAC"]) <= 0):
385                 print("Empty calibration file for DAC")
386                 return False
387         return True
388
389 def getADC_Voltage(channel, counts):
390         if (calibrate["ADC"][channel] == None or len(calibrate["ADC"][channel]) <= 0):
391                 if (calibrate["ADC_Rin"][channel] != None and calibrate["ADC_Rvar"][channel] != None):
392                         print("Warning: Using rough calibration for ADC"+str(channel) + " = " + str(counts))
393                         ratio = float(calibrate["ADC_Rin"][channel]) / (float(calibrate["ADC_Rin"][channel]) + float(calibrate["ADC_Rvar"][channel]))
394                         return ratio * (float(counts) / float(calibrate["ADC_Counts"][channel]) * Vref)
395                 else:
396                         print("Error: No calibration for ADC"+str(channel))
397                         return False
398
399         c = calibrate["ADC"][channel]
400         valueIndex = 1
401         for i in range(0, len(c)-1):
402                 if (c[i][0] <= counts and i + 1 < len(c)):
403                         grad = (float(c[i+1][valueIndex]) - float(c[i][valueIndex])) / (float(c[i+1][0]) - float(c[i][0]))
404                         value = float(c[i][valueIndex]) + grad * (float(counts) - float(c[i][0]))
405                         return value
406
407         
408         print("Warning: Extrapolating outside calibration range for DAC = " + str(counts))
409
410         grad = (float(c[len(c)-1][valueIndex]) - float(c[len(c)-2][valueIndex])) / (float(c[len(c)-1][0]) -  float(c[len(c)-2][0]))
411         value = float(c[len(c)-1][valueIndex]) + grad * (float(counts) - float(c[len(c)-1][0]))
412                 
413 def readADC(channel):
414         #for i in range(0, aquire["ADC_Averages"]):
415         #       ser.write("r "+str(channel)+"\r") # Send command to datalogger
416                 #time.sleep(aquire["response_wait"])
417         #       response = ser.readline().strip("#\r\n ")
418         #       if (response != "r "+str(channel)):
419         #               print("Received wierd response reading ADC ("+response+")")
420         #               return False                    
421                 #time.sleep(aquire["response_wait"])
422         #       values.append(int(ser.readline().strip("\r\n ")))
423         #       adc_sum += float(values[len(values)-1])
424         ser.write("r "+str(channel)+"\r")
425         response = ser.readline().strip("#\r\n")
426         if (response != "r "+str(channel)):
427                 print("Received wierd response reading ADC ("+response+")")
428                 return False                    
429         return ser.readline().strip("\r\n").split(" ")
430
431 def getDAC_Voltage(counts, gain = True):
432         if (calibrate["DAC"] == None or len(calibrate["DAC"]) <= 0):
433                 if (calibrate["DAC_Gain"] != None):
434                         print("Warning: Using rough calibration for DAC = " + str(counts))
435                         return float(counts) / float(calibrate["DAC_Counts"]) * float(calibrate["Vref"]) * float(calibrate["DAC_Gain"])
436                 else:
437                         print("Error: No calibrate for DAC")
438                         return False
439
440         valueIndex = 1
441         if (gain == False):
442                 valueIndex = 2
443                 if (len(calibrate["DAC"][0]) < 3):
444                         print("Error: No data for unamplified DAC")
445                         return False
446         
447         c = calibrate["DAC"]
448         if (len(c) == 1):
449                 print("Warning: Only one point in calibration data!")
450                 return float(c[0][valueIndex])
451         
452         for i in range(0, len(c)-1):
453                 if (c[i][0] <= counts and i + 1 < len(c)):
454                         grad = (float(c[i+1][valueIndex]) - float(c[i][valueIndex])) / (float(c[i+1][0]) - float(c[i][0]))
455                         value = float(c[i][valueIndex]) + grad * (float(counts) - float(c[i][0]))
456                         return value
457
458         
459         print("Warning: Extrapolating outside calibration range for DAC = " + str(counts))
460
461         grad = (float(c[len(c)-1][valueIndex]) - float(c[len(c)-2][valueIndex])) / (float(c[len(c)-1][0]) -  float(c[len(c)-2][0]))
462         value = float(c[len(c)-1][valueIndex]) + grad * (float(counts) - float(c[len(c)-1][0]))
463         return value
464
465
466 def setDAC(level):
467         
468         ser.write("d "+str(level)+"\r") 
469         
470         response = ser.readline().strip("#\r\n ")
471         if (response != "d "+str(level)):
472                 print("Recieved wierd response setting DAC to "+str(level) + " ("+response+")")
473                 return False
474         #time.sleep(aquire["response_wait"])
475         #time.sleep(aquire["DAC_Settle"])
476         #time.sleep(aquire["DAC_Settle"])
477         response = ser.readline().strip("#\r\n ")
478         if (response != "DAC "+str(level)):
479                 print("Recieved wierd response setting DAC to "+str(level) + " ("+response+")")
480                 return False
481         #time.sleep(aquire["response_wait"])
482         time.sleep(aquire["DAC_Settle"])
483                 
484
485 def calibrateADC_usingDAC(channel, gain = True):
486         if (calibrate["DAC"] == None):
487                 print("ERROR: DAC is not calibrated, can't use to calibrate ADC!")
488                 return False
489
490         calibrate["ADC"][channel] = []
491         outfile = log_open("calibrateADC"+str(channel)+".dat", "w")
492         outfile.write("# Calibrate ADC " + str(channel) + "\n")
493         outfile.write("# Start " + str(datetime.datetime.now()) + "\n")
494
495         print("Connect DAC output to ADC " + str(channel) + " and press enter\n")
496         sys.stdin.readline()
497         
498         for dac in calibrate["DAC"]:
499                 if (setDAC(dac[0]) == False):
500                         return False
501                 
502                 value = readADC(channel)
503                 if (value == False):
504                         return False
505                 canadd = (len(calibrate["ADC"][channel]) <= 0)
506                 if (canadd == False):
507                         canadd = True
508                         for adc in calibrate["ADC"][channel]:
509                                 if adc[0] == value[0]:
510                                         canadd = False
511                                         break
512
513                 if (canadd):            
514
515                         if gain:
516                                 input_value = dac[1]
517                         else:
518                                 input_value = dac[2]
519
520                         #input_value = getDAC_Voltage(dac[0], gain)
521                         if (input_value == False):
522                                 return False
523         
524                         outfile.write(str(value[0]) + "\t" + str(input_value) + "\t" + str(value[1]) + "\n")
525                         print(str(value[0]) + "\t" + str(input_value) + "\t" + str(value[1]))
526
527                         calibrate["ADC"][channel].append((value[0], input_value, value[1]))
528
529         outfile.write("# Stop " + str(datetime.datetime.now()) + "\n# EOF\n")
530         log_close(outfile)
531         if (setDAC(0) == False):
532                 return False
533         if (len(calibrate["ADC"][channel]) <= 0):
534                 print("Error: No calibration points taken for ADC " + str(channel))
535                 return False
536         return True
537
538 def calibrateADC(channel):      
539         calibrate["ADC"][channel] = []
540         outfile = log_open("calibrateADC"+str(channel)+".dat", "w")
541         outfile.write("# Calibrate ADC " + str(channel) + "\n")
542         outfile.write("# Start " + str(datetime.datetime.now()) + "\n")
543
544         print("Calibrating ADC...\n")
545         print("Enter measured voltage, empty line stops.\n")
546
547         while True:
548                 read = sys.stdin.readline().strip("\r\n ")
549                 if (read == ""):
550                         break
551                 input_value = float(read)
552                 output_value = readADC(channel)
553                 if (output_value == False):
554                         return False
555                 
556                 calibrate["ADC"][channel].append((output_value[0], input_value))
557                 print(str(output_value[0]) + "\t" + str(input_value))
558                 
559         calibrate["ADC"][channel].sort(lambda e : e[1], reverse = False)
560         
561         
562         
563
564         outfile.write("# Stop " + str(datetime.datetime.now()) + "\n# EOF\n")
565         log_close(outfile)
566         if (len(calibrate["ADC"][channel]) <= 0):
567                 print("Error: No calibration points taken for ADC " + str(channel))
568                 return False
569
570         return True
571
572 def calibrateDAC():
573         calibrate["DAC"] = []
574
575         outfile = log_open("calibrateDAC.dat", "w")
576         outfile.write("# Calibrate DAC\n")
577         outfile.write("# Start " + str(datetime.datetime.now()) + "\n")
578
579         print("Calibrating DAC...")
580         sys.stdout.write("Number of counts to increment per step: ")
581         read = sys.stdin.readline().strip("\r\n")
582         if (read == ""):
583                 return False
584         stepSize = max(1, int(read))
585         sys.stdout.write("\n")
586
587         outfile.write("# Increment by " + str(stepSize) + "\n")
588         print("Input measured DAC voltage for each level; empty input ends calibration.")
589         print("You may optionally input the DAC voltage before it is amplified too.")
590
591         level = 0
592         while (level < calibrate["DAC_Counts"]):
593
594                 setDAC(level)
595
596                 sys.stdout.write(str(level) + " ?")
597                 read = sys.stdin.readline().strip("\r\n ")
598                 if (read == ""):
599                         break
600                 else:
601                         read = read.split(" ")
602                         if (len(calibrate["DAC"]) > 0 and len(calibrate["DAC"][len(calibrate["DAC"])-1]) != len(read)):
603                                 print("Either give one value for EVERY point, or two values for EVERY point. Don't mess around.")
604                                 return False
605                         if (len(read) == 2):
606                                 calibrate["DAC"].append((level, float(read[0]), float(read[1])))
607                                 outfile.write(str(level) + "\t" + str(float(read[0])) + "\t" + str(float(read[1])) + "\n")
608                         else:
609                                 calibrate["DAC"].append((level, float(read[0])))
610                                 outfile.write(str(level) + "\t" + str(float(read[0])) + "\n")
611                 
612
613
614                 level += stepSize
615
616         outfile.write("# Stop " + str(datetime.datetime.now()) + "\n# EOF\n")
617         log_close(outfile)
618         if (len(calibrate["DAC"]) <= 0):
619                 print("Error: No calibration points taken for DAC")
620                 return False
621         return setDAC(0)
622
623 # Run the main function
624 if __name__ == "__main__":
625         sys.exit(main())

UCC git Repository :: git.ucc.asn.au