#! /bin/csh
alias sleep1 'sleep 1'
caput 8ide:asyn_17.TMOD "Write"
sleep1    

caput 8ide:asyn_17.AOUT "*RST"
sleep1    
caput 8ide:asyn_17.AOUT ":SENS:FUNC:CONC OFF"
sleep1
caput 8ide:asyn_17.AOUT ":SOUR:FUNC VOLT"
sleep1
caput 8ide:asyn_17.AOUT ":SENS:FUNC 'CURR:DC'"
sleep1
caput 8ide:asyn_17.AOUT ":SENS:CURR:PROT 100E-6"
sleep1
caput 8ide:asyn_17.AOUT ":SOUR:VOLT:START -4.0"
sleep1
caput 8ide:asyn_17.AOUT ":SOUR:VOLT:STOP 4.0"
sleep1
caput 8ide:asyn_17.AOUT ":SOUR:VOLT:STEP 10E-2"
sleep1
caput 8ide:asyn_17.AOUT ":SOUR:VOLT:MODE SWE"
sleep1
caput 8ide:asyn_17.AOUT ":SOUR:SWE:RANG AUTO"
sleep1
caput 8ide:asyn_17.AOUT ":SOUR:SWE:SPAC LIN"
sleep1
caput 8ide:asyn_17.AOUT ":SOUR:SWE:POIN 81"
sleep1
caput 8ide:asyn_17.AOUT ":TRIG:COUN 81"
sleep1
caput 8ide:asyn_17.AOUT ":SOUR:DEL 0.1"
sleep1
#caput 8ide:asyn_17.AOUT ":OUTP ON"
# The following line should turn the output on/off automatically 
# when a sweep is initiated and is meant to replace the ":OUTP ON" command
caput 8ide:asyn_17.AOUT "SOUR:CLE:AUTO ON"     
sleep1 
caput 8ide:asyn_17.AOUT ":FORM:ELEM VOLT,CURR"
sleep1

#!/APSshare/anaconda3/x86_64/bin/python

import epics
from time import sleep
import sys
arg2 = sys.argv
outpath = "/home/8-id-g/2018-3/tsai201810p/device/"
outfile = outpath + arg2[1]

sweepit = 1
plotit = 0


if sweepit : 
	epics.caput("8ide:2400:ReadArray.PROC", 1)
	# wait for sweep to complete -- hack with fixed delay. Is there a condition to test so we know when sweep is complete?
	print('sweep in progress')
	sleep(20)
	# close fast shutter
	epics.caput("8idg:Unidig1Bo13",1)


# read in the sweep data
indata = epics.caget("8ide:2400:ReadArray");
# read in number of non-zero elements.  
num2pts = epics.caget("8ide:2400:ReadArray.NORD");
print(num2pts);  # development
# Why doesn't this agree?   For 41 pts (82 values), returns 752.
print(sys.getsizeof(indata))  # development

print(indata)  # development

f = open(outfile,"w")
counter = 0
for i in indata :

	#print(i)
	counter=counter+1;
	if (counter % 2 == 1):
		try:
			VOLT = float(i)
		except:
			VOLT = "NaN"
			print("Error converting string to float./n")
	else :
		try:
			CURR = float(i)
		except:
			CURR = "NaN"
			print("Error converting string to float./n")

		print("VOLTAGE=",VOLT," CURRENT=",CURR)
		f.write('{0:>6} {1}\n'.format(VOLT ,CURR))

f.close()
print(arg2[1])

if plotit :
		print("Plot the data here.")

#! /bin/csh
# Set up filter (singal averaging) for sweep scan of Keithely 2400 sourcemeter
# see Keithley 2400 Series Sourcemeter manual, Table 10-4, Staircase Sweep Programming example
# delay will be needed to process each asyn input
alias sleep1 'sleep 1'
# Put asyn comms in Write mode
caput 8ide:asyn_17.TMOD "Write"
sleep1    
# next line chooses filter: REPeating filter averages signal value for each source value
#                           MOVing filter averages together signal for different source values
caput 8ide:asyn_17.AOUT ":SENS:AVER:TCON REP"
sleep1
# next line sets how many values to average together
caput 8ide:asyn_17.AOUT ":SENS:AVER:COUN 10"
sleep1
# next line turns averaging ON or OFF
caput 8ide:asyn_17.AOUT ":SENS:AVER ON"
sleep1

#!/APSshare/anaconda3/x86_64/bin/python

import epics
from time import sleep
import sys
arg2 = sys.argv
outpath = "/home/8-id-g/2018-3/tsai201810p/device/"
outfile = outpath + arg2[1]

sweepit = 1
plotit = 0
verbose = 0
countup = 1

# see documentation for PyEpics release 3.2.4
# https://www.slac.stanford.edu/grp/ssrl/spear/epics/extensions/pyepics/pyepics.pdf
# section 3.1 "The PV class", 3.1.1 for the get() method and 3.1.2 for the timestamp attribute
numelem = epics.PV("8ide:2400:ReadArray.NORD");
data2400 = epics.PV("8ide:2400:ReadArray");

print( data2400.get()[1] )
print( data2400.timestamp)

oldtime = data2400.timestamp;
print("oldtime : ", oldtime)

if sweepit :
	
	data2400.get();
	checktime = data2400.timestamp;
	 
	epics.caput("8ide:2400:ReadArray.PROC", 1)
	# wait for sweep to complete -- hack with fixed delay. Is there a condition to test so we know when sweep is complete?
	print('sweep in progress')
	
	sweepct = 0;
	while (checktime == oldtime) :
		sleep(0.5);
		sweepct=sweepct+1;
		print('waiting for sweep: ', sweepct)
		data2400.get()
		checktime = data2400.timestamp;
		if countup:
			print('checktime: ', checktime);	
	
	# close fast shutter
	epics.caput("8idg:Unidig1Bo13",1)

# read in the sweep data
indata = epics.caget("8ide:2400:ReadArray");
# read in number of non-zero elements.  
num2pts = epics.caget("8ide:2400:ReadArray.NORD");
print(num2pts);  # development
# Why doesn't this agree?   For 41 pts (82 values), returns 752.
print(sys.getsizeof(indata))  # development

if verbose:
	print(indata)  # development

f = open(outfile,"w")
counter = 0
for i in indata :

	#print(i)
	counter=counter+1;
	if (counter % 2 == 1):
		try:
			VOLT = float(i)
		except:
			VOLT = "NaN"
			print("Error converting string to float./n")
	else :
		try:
			CURR = float(i)
		except:
			CURR = "NaN"
			print("Error converting string to float./n")
		if verbose :
			print("VOLTAGE=",VOLT," CURRENT=",CURR)
		f.write('{0:>6} {1}\n'.format(VOLT ,CURR))

f.close()
print(arg2[1])

if plotit :
		print("Plot the data here.")