#!/usr/bin/python
import spidev
import RPi.GPIO as GPIO
import time

import numpy as np
import matplotlib
matplotlib.use('Agg')
import matplotlib.pyplot as plt

def StartUp():
    GPIO.setmode(GPIO.BCM)
    PRESET = 25
    IO4 = 26
    GPIO.setup(PRESET,GPIO.OUT)
    GPIO.setup(IO4,GPIO.IN, pull_up_down=GPIO.PUD_DOWN)
    print "Reset 25 - Low 1s"
    GPIO.output(PRESET,GPIO.LOW)
    time.sleep(3)
    print "Reset 25 - High 0.2s"
    GPIO.output(PRESET,GPIO.HIGH)
    time.sleep(0.2)

    spi = spidev.SpiDev()
    spi.open(0,1) # CS2 - FPGA, on CE1 = IO4
    spi.mode = 0b01
    print "spi.cshigh is " + str(spi.cshigh)
    print "spi mode is " + str(spi.mode)
    spi.max_speed_hz = 2000000
    print "spi maxspeed is "+str(spi.max_speed_hz)+"hz"
    return spi

def TestSPI(spi,ncycles):
    i = 0
    while i < ncycles:
        WriteFPGA(spi,0xEB,0x01) # 0: single mode 1 continious mode
        time.sleep(0.5)
        WriteFPGA(spi,0xEB,0x00) # 0: single mode 1 continious mode
        time.sleep(0.5)  
        i = i+1

def LoopSPI(spi):
    while 1:
        WriteFPGA(spi,0xEB,0x01) # 0: single mode 1 continious mode
        WriteFPGA(spi,0xEB,0x00) # 0: single mode 1 continious mode

def LoopAcq(spi):
    while 1:
        WriteFPGA(spi,0xEB,0x00) # Doing 1 shot 
        WriteFPGA(spi,0xEF,0x01) # Cleaning memory pointer
        WriteFPGA(spi,0xEA,0x01) # Software Trig : As to be clear by software
        time.sleep(0.1) # sleep 10ms

def ClearMem(spi):
    WriteFPGA(spi,0xEF,0x01) # To access memory

def ConfigSPI(spi):
    WriteFPGA(spi,0xE0,0x14) # set sEEPon
    WriteFPGA(spi,0xE1,0x00) # set sEEPoff MSB
    WriteFPGA(spi,0xE2,0x35) # set sEEPoff LSB
    WriteFPGA(spi,0xE3,0x00) # set sEEDelayACQ sEEDelayACQ MSB
    WriteFPGA(spi,0xE4,0x05) # set sEEDelayACQ EEDelayACQ LSB
    WriteFPGA(spi,0xE5,0x32) # set sEEPeriod MSB // durée de l'acquisition ? 
    WriteFPGA(spi,0xE6,0xC8) # set sEEPeriod LSB 0x32C8 acq (13000t * 10ns = 130us)
    WriteFPGA(spi,0xE7,0x00) # Period of one cycle MSB
    WriteFPGA(spi,0xE8,0xC3) # Period of one cycle 15 to 8
    WriteFPGA(spi,0xE9,0x50) # Period of one cycle LSB
    #WriteFPGA(spi,0xEA,0x00) # Software Trig : As to be clear by software
    WriteFPGA(spi,0xEB,0x00) # 0: single mode 1 continious mode
    WriteFPGA(spi,0xEC,0x11) # Voltage gain control: 0V to 1V
    WriteFPGA(spi,0xED,0x03) # Frequency of ADC acquisition / sEEADC_freq (3 = 16Msps, 1 = 32, 0 = 64, 2 = 21Msps)
    WriteFPGA(spi,0xEE,0xA0) # How many cycles in countinious mode
    print "Config FPGA done!"

def WriteFPGA(spi,adress,value):
    spi.xfer([0xAA] )
    spi.xfer([adress] )
    spi.xfer([value] )

spi = StartUp()
ConfigSPI(spi)
N = 0

Reset 25 - Low 1s
Reset 25 - High 0.2s
spi.cshigh is False
spi mode is 1
spi maxspeed is 2000000hz
Config FPGA done!

#LoopSPI(spi)

#LoopAcq(spi)

#WriteFPGA(spi,0xEA,0x01) # Software Trig : As to be clear by software

#WriteFPGA(spi,0xEA,0x01) # Cleaning memory pointer

TestSPI(spi,3) # LED2 clignote 3x

# Set PulseOn
POn = 200 #ns
HPon = POn/10
print hex(HPon)
WriteFPGA(spi,0xE0,HPon) # set sEEPon
# Set Lengh between Pon and Poff
DeltaPP = 100 #ns
HPP =DeltaPP/10
print  hex(HPP)
WriteFPGA(spi,0xD0,HPP) # set sEEPon
# Setting Poff
sEEPoff = 2000 #ns
POff = sEEPoff/10 
POffMSB, POffLSB = 0xFF00&POff/0xFF,0x00FF&POff 
print "Poff:", hex(POffMSB),hex(POffLSB)
WriteFPGA(spi,0xE1,POffMSB) # set sEEPon MSB
WriteFPGA(spi,0xE2,POffLSB) # set sEEPon LSB

# Setting Poff to Acq delay sEEDelayACQ

DeltaAcq = 7000 #ns
hDA = DeltaAcq/10 
hDAMSB, hDALSB = hDA/255 , 0x00FF&hDA 
print "hDA:", hex(hDAMSB),hex(hDALSB)
WriteFPGA(spi,0xE3,hDAMSB) # set sEEPon MSB
WriteFPGA(spi,0xE4,hDALSB) # set sEEPon LSB

# sEEACQ
LAcq = 130 #us
LAcqNs = LAcq*1000/10 #ns
LAcqMSB, LAcqLSB = LAcqNs/256,0x00FF&LAcqNs 
print "Acq:",  hex(LAcqMSB),hex(LAcqLSB)
WriteFPGA(spi,0xE5,hDAMSB) # set sEEPon MSB
WriteFPGA(spi,0xE6,hDALSB) # set sEEPon LSB

# sEEPeriod
lEPeriod = 1000 #us
lEPNs = lEPeriod*1000/10 #ns
EPNsMSB, EPNs, EPNsLSB = 0x00FF&lEPNs/(256*256),0x00FF&lEPNs/256,0x0000FF&lEPNs 
print "Period:",  hex(EPNsMSB),hex(EPNs),hex(EPNsLSB) 
WriteFPGA(spi,0xE7,EPNsMSB) # Period of one cycle MSB
WriteFPGA(spi,0xE8,EPNs) # Period of one cycle 15 to 8
WriteFPGA(spi,0xE9,EPNsLSB) # Period of one cycle LSB

# sEETrigCounter
NAcTrig = 10
WriteFPGA(spi,0xEE,NAcTrig) # Period of one cycle LSB

0x14
0xa
Poff: 0x0 0xc8
hDA: 0x2 0xbc
Acq: 0x32 0xc8
Period: 0x1 0x86 0xa0

f = 0x01
LAcq = 130 #us
WriteFPGA(spi,0xED,f) # Frequency of ADC acquisition / sEEADC_freq (3 = 16Msps, 1 = 32, 0 = 64, 2 = 21Msps)
VGA = 0x44
NAcTrig = 200 # 200 acquisitions
WriteFPGA(spi,0xEC,VGA)
WriteFPGA(spi,0xEB,0x01) # Doing several lines
WriteFPGA(spi,0xEE,NAcTrig) # Acquire 100 lines
WriteFPGA(spi,0xEF,0x01) # Cleaning memory pointer
WriteFPGA(spi,0xEA,0x01) # Software Trig : As to be clear by software

time.sleep(0.2) # sleep 10ms

Fech = int(64/((1+f)))
print Fech
Nacq = LAcq * Fech * NAcTrig
print "-> "+str(Nacq) + ' samples'

10
-> 260000 samples

LAcq * Fech

1300

spi.max_speed_hz = 9000000
A = []
limit = 2*Nacq-10
x = 20000
print limit
print limit/x
for i in range(limit):
    A.append ( spi.xfer([0x00] )[0] )  
    if not(i%x):
        print len(A)/x
a = np.asarray(A).astype(int)
print 'acq done'
np.savetxt( "FI-"+str(NAcTrig)+"-"+str(N)+"-VGA@"+str(hex(VGA))+"-spimode"+str(spi.mode)+"-"+str(Fech)+"msps.csv", a.astype(int), fmt='%i', delimiter=";")
N = N+1

519990
25
0
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
acq done

519990/2

259995

spi.close()

f = 0x02
WriteFPGA(spi,0xED,f) # Frequency of ADC acquisition / sEEADC_freq (3 = 16Msps, 1 = 32, 0 = 64, 2 = 21Msps)
VGA = 0x22
WriteFPGA(spi,0xEC,VGA)
WriteFPGA(spi,0xEB,0x00) # Doing several lines 
WriteFPGA(spi,0xEF,0x01) # Cleaning memory pointer
WriteFPGA(spi,0xEA,0x01) # Software Trig : As to be clear by software

time.sleep(0.2) # sleep 10ms

Fech = int(64/((1+f)))
print Fech
Nacq = LAcq * Fech
print "-> "+str(Nacq) + ' samples'

21
-> 2730 samples

A = []
for i in range(2*Nacq-10):
    A.append ( spi.xfer([0x00] )[0] )  
a = np.asarray(A).astype(int)

np.savetxt( "One-"+str(N)+"-VGA@"+str(hex(VGA))+"-spimode"+str(spi.mode)+"-"+str(Fech)+"msps.csv", a, delimiter=";")
N = N+1

LAcq * 64.0/3

2773.3333333333335