Skip to content
Projects
Groups
Snippets
Help
Loading...
Sign in
Toggle navigation
P
Production Test Suite
Project
Project
Details
Activity
Cycle Analytics
Repository
Repository
Files
Commits
Branches
Tags
Contributors
Graph
Compare
Charts
Issues
9
Issues
9
List
Board
Labels
Milestones
Merge Requests
0
Merge Requests
0
Wiki
Wiki
image/svg+xml
Discourse
Discourse
Members
Members
Collapse sidebar
Close sidebar
Activity
Graph
Charts
Create a new issue
Commits
Issue Boards
Open sidebar
Projects
Production Test Suite
Commits
b8b8970c
Commit
b8b8970c
authored
Apr 25, 2014
by
Matthieu Cattin
Browse files
Options
Browse Files
Download
Email Patches
Plain Diff
test43: Add saturation test for no wraparound.
parent
8398db94
Hide whitespace changes
Inline
Side-by-side
Showing
1 changed file
with
406 additions
and
0 deletions
+406
-0
test43.py
test/fmcadc100m14b4cha/python/test43.py
+406
-0
No files found.
test/fmcadc100m14b4cha/python/test43.py
0 → 100755
View file @
b8b8970c
#! /usr/bin/env python
# coding: utf8
# Copyright CERN, 2011
# Author: Matthieu Cattin <matthieu.cattin@cern.ch>
# Licence: GPL v2 or later.
# Website: http://www.ohwr.org
# Last modifications: 30/5/2012
# Import system modules
import
sys
import
time
import
os
# Add common modules and libraries location to path
sys
.
path
.
append
(
'../../../'
)
sys
.
path
.
append
(
'../../../gnurabbit/python/'
)
sys
.
path
.
append
(
'../../../common/'
)
# Import common modules
from
ptsexcept
import
*
import
rr
# Import specific modules
from
fmc_adc_spec
import
*
from
fmc_adc
import
*
from
numpy
import
*
from
pylab
import
*
from
calibr_box
import
*
import
find_usb_tty
from
PAGE.Agilent33250A
import
*
from
PAGE.SineWaveform
import
*
"""
test43: Tests saturation (-> no wraparound)
"""
NB_CHANNELS
=
4
AWG_SET_SLEEP
=
0.3
SSR_SET_SLEEP
=
0.05
BOX_SET_SLEEP
=
0.01
ACQ_TIMEOUT
=
10
PRE_TRIG_SAMPLES
=
10
POST_TRIG_SAMPLES
=
200
NB_SHOTS
=
1
BYTES_PER_SAMPLE
=
2
TRIG_TIMETAG_BYTES
=
16
RANGES
=
[
'10V'
,
'1V'
,
'100mV'
]
def
open_all_channels
(
fmc
):
for
i
in
range
(
1
,
NB_CHANNELS
+
1
):
fmc
.
set_input_range
(
i
,
'OPEN'
)
time
.
sleep
(
SSR_SET_SLEEP
)
def
fmc_adc_init
(
spec
,
fmc
):
print
"Initialise FMC board.
\n
"
# Reset offset DACs
fmc
.
dc_offset_reset
()
# Make sure all switches are OFF
open_all_channels
(
fmc
)
# Set software trigger
fmc
.
set_soft_trig
()
# Set acquisition
fmc
.
set_pre_trig_samples
(
PRE_TRIG_SAMPLES
)
fmc
.
set_post_trig_samples
(
POST_TRIG_SAMPLES
)
fmc
.
set_shots
(
NB_SHOTS
)
# Converts two's complement hex to signed
def
hex2signed
(
value
):
if
(
value
&
0x8000
):
return
-
((
~
value
&
0xFFFF
)
+
1
)
else
:
return
(
value
&
0xFFFF
)
# Converts digital value to volts
def
digital2volt
(
value
,
full_scale
,
nb_bit
):
return
float
(
value
)
*
float
(
full_scale
)
/
2
**
nb_bit
# Converts hex gain value to float
def
gain2float
(
value
):
dec
=
(
value
&
0x8000
)
>>
15
frac
=
value
&
0x7FFF
return
(
float
)
(
dec
+
(
frac
*
1.0
/
2
**
15
))
def
get_corr_values
(
fmc
):
off_corr
=
fmc
.
get_adc_offset_corr
(
1
)
print
(
"Offset corr:0x
%04
X (
%
d)"
%
(
off_corr
,
hex2signed
(
off_corr
)))
gain_corr
=
fmc
.
get_adc_gain_corr
(
1
)
print
(
"Gain corr :0x
%04
X (
%1.6
f)"
%
(
gain_corr
,
gain2float
(
gain_corr
)))
def
acq_channels
(
fmc
,
carrier
,
pause
):
# Make sure no acquisition is running
fmc
.
stop_acq
()
time
.
sleep
(
pause
)
# Start acquisition
fmc
.
start_acq
()
time
.
sleep
(
pause
)
# Trigger
fmc
.
sw_trig
()
# Wait end of acquisition
timeout
=
0
while
(
'IDLE'
!=
fmc
.
get_acq_fsm_state
()):
time
.
sleep
(
.1
)
timeout
+=
1
if
(
ACQ_TIMEOUT
<
timeout
):
print
"Acquisition timeout. Missing trigger?."
print
"Acq FSm state:
%
s"
%
fmc
.
get_acq_fsm_state
()
return
1
# Retrieve data trough DMA
trig_pos
=
fmc
.
get_trig_pos
()
# Enable "DMA done" interrupt
carrier
.
enable_dma_done_irq
()
# Read samples for all channels + trigger timetag
data_length
=
((
PRE_TRIG_SAMPLES
+
1
+
POST_TRIG_SAMPLES
)
*
NB_CHANNELS
*
BYTES_PER_SAMPLE
)
+
TRIG_TIMETAG_BYTES
channels_data
=
carrier
.
get_data
((
trig_pos
-
(
PRE_TRIG_SAMPLES
*
8
)),
data_length
)
trig_timetag
=
[]
data
=
[]
for
i
in
range
(
8
):
data
.
append
(
channels_data
.
pop
(
-
1
))
for
i
in
range
(
0
,
8
,
2
):
trig_timetag
.
append
(((
data
[
i
]
<<
16
)
+
data
[
i
+
1
]))
# Disable "DMA done" interrupt
carrier
.
disable_dma_done_irq
()
#print("raw hex data: 0x%08X"%(channels_data[0]))
channels_data
=
[
hex2signed
(
item
)
for
item
in
channels_data
]
#print("signed data : 0x%08X (%d)"%(channels_data[0], channels_data[0]))
return
channels_data
,
trig_timetag
def
plot_channel
(
ch_data_1
,
ch_data_2
,
ch_data_3
,
ch_data_4
,
sat
,
ylimit
):
sample
=
arange
(
len
(
ch_data_1
))
plot
(
sample
,
ch_data_1
,
'go-'
,
label
=
'reset'
)
plot
(
sample
,
ch_data_2
,
'bo-'
,
label
=
'dc offset to min'
)
plot
(
sample
,
ch_data_3
,
'mo-'
,
label
=
'neg adc offset corr'
)
plot
(
sample
,
ch_data_4
,
'co-'
,
label
=
'with saturation'
)
plot
(
sample
,
[
sat
]
*
len
(
sample
),
'm'
)
plot
(
sample
,
[
-
sat
]
*
len
(
sample
),
'm'
)
plot
(
sample
,
[
32764
]
*
len
(
sample
),
'g'
)
plot
(
sample
,
[
-
32765
]
*
len
(
sample
),
'g'
)
ylim
(
-
ylimit
-
(
ylimit
/
10.0
),
ylimit
+
(
ylimit
/
10.0
))
grid
(
color
=
'k'
,
linestyle
=
':'
,
linewidth
=
1
)
legend
(
loc
=
'upper left'
)
xlabel
(
'Samples'
)
ylabel
(
'Raw ADC values (two
\'
s complement)'
)
title
(
'Variable saturation, limit:
%
d'
%
sat
)
#draw()
show
()
return
0
def
main
(
default_directory
=
'.'
):
# Constants declaration
TEST_NB
=
43
FMC_ADC_BITSTREAM
=
'../firmwares/spec_fmcadc100m14b4cha.bin'
FMC_ADC_BITSTREAM
=
os
.
path
.
join
(
default_directory
,
FMC_ADC_BITSTREAM
)
EXPECTED_BITSTREAM_TYPE
=
0x1
# Calibration box vendor and product IDs
BOX_USB_VENDOR_ID
=
0x10c4
# Cygnal Integrated Products, Inc.
BOX_USB_PRODUCT_ID
=
0xea60
# CP210x Composite Device
# Agilent AWG serial access vendor and product IDs
AWG_USB_VENDOR_ID
=
0x0403
# Future Technology Devices International, Ltd
AWG_USB_PRODUCT_ID
=
0x6001
# FT232 USB-Serial (UART) IC
AWG_BAUD
=
57600
EEPROM_BIN_FILENAME
=
"eeprom_content.out"
EEPROM_BIN_FILENAME
=
os
.
path
.
join
(
default_directory
,
EEPROM_BIN_FILENAME
)
EEPROM_SIZE
=
8192
# in Bytes
CALIBR_BIN_FILENAME
=
"calibration_data.bin"
CALIBR_BIN_FILENAME
=
os
.
path
.
join
(
default_directory
,
CALIBR_BIN_FILENAME
)
start_test_time
=
time
.
time
()
print
"================================================================================"
print
"Test
%02
d start
\n
"
%
TEST_NB
# SPEC object declaration
print
"Loading hardware access library and opening device.
\n
"
spec
=
rr
.
Gennum
()
# Load FMC ADC firmware
print
"Loading FMC ADC firmware:
%
s
\n
"
%
FMC_ADC_BITSTREAM
spec
.
load_firmware
(
FMC_ADC_BITSTREAM
)
time
.
sleep
(
2
)
# Carrier object declaration (SPEC board specific part)
# Used to check that the firmware is loaded.
try
:
carrier
=
CFmcAdc100mSpec
(
spec
,
EXPECTED_BITSTREAM_TYPE
)
except
FmcAdc100mSpecOperationError
as
e
:
raise
PtsCritical
(
"Carrier init failed, test stopped:
%
s"
%
e
)
# Mezzanine object declaration (FmcAdc100m14b4cha board specific part)
try
:
fmc
=
CFmcAdc100m
(
spec
)
except
FmcAdc100mOperationError
as
e
:
raise
PtsCritical
(
"Mezzanine init failed, test stopped:
%
s"
%
e
)
try
:
# Initialise fmc adc
fmc_adc_init
(
spec
,
fmc
)
# Use data pattern instead of ADC data
#fmc.testpat_en(0x1FFF) # max
#fmc.testpat_en(0x0) # mid
#fmc.testpat_en(0x2000) # min
# Set UTC
current_time
=
time
.
time
()
utc_seconds
=
int
(
current_time
)
fmc
.
set_utc_second_cnt
(
utc_seconds
)
#print "UTC core seconds counter initialised to : %d" % fmc.get_utc_second_cnt()
utc_coarse
=
int
((
current_time
-
utc_seconds
)
/
8E-9
)
fmc
.
set_utc_coarse_cnt
(
utc_coarse
)
#print "UTC core coarse counter initialised to : %d" % fmc.get_utc_coarse_cnt()
# Print configuration
#fmc.print_adc_core_config()
# Print ADC config
#fmc.print_adc_config()
# Acquisition parameters
ACQ_PAUSE
=
1
# pause between acq. stop and start, start and trigger
IN_RANGE
=
'CAL'
IN_TERM
=
'OFF'
##################################################
# Reset DC offset
##################################################
fmc
.
dc_offset_reset
()
print
(
"
\n
Reset DC offset."
)
##################################################
# Acquire channel 1
##################################################
print
"Acquiring channel 1"
# Configure analogue input
fmc
.
set_input_range
(
1
,
IN_RANGE
)
fmc
.
set_input_term
(
1
,
IN_TERM
)
time
.
sleep
(
SSR_SET_SLEEP
)
# Perform an acquisition
acq_data
,
trig_timetag
=
acq_channels
(
fmc
,
carrier
,
ACQ_PAUSE
)
channel_data
=
acq_data
[
0
::
4
]
#print("Number of samples: %d"%(len(channel_data)))
# Calculate middle value
ch_max
=
max
(
channel_data
)
ch_min
=
min
(
channel_data
)
ch_mid
=
(
abs
(
ch_max
)
-
abs
(
ch_min
))
ch_mean
=
mean
(
channel_data
)
print
(
"Channel max=
%
d, min=
%
d, mid=
%
d, mean=
%
d"
%
(
ch_max
,
ch_min
,
ch_mid
,
ch_mean
))
ch_data_1
=
channel_data
# Plot channel
#plot_channel(channel_data, ch_mean, 32768)
##################################################
# Apply a DC offset
##################################################
dc_offset
=
0xFFFF
print
(
"
\n
Set dc offset to: 0x
%08
X"
%
(
dc_offset
))
fmc
.
set_dc_offset
(
1
,
dc_offset
)
##################################################
# Acquire channel 1
##################################################
print
"Acquiring channel 1"
# Configure analogue input
fmc
.
set_input_range
(
1
,
IN_RANGE
)
fmc
.
set_input_term
(
1
,
IN_TERM
)
time
.
sleep
(
SSR_SET_SLEEP
)
# Perform an acquisition
acq_data
,
trig_timetag
=
acq_channels
(
fmc
,
carrier
,
ACQ_PAUSE
)
channel_data
=
acq_data
[
0
::
4
]
#print("Number of samples: %d"%(len(channel_data)))
# Calculate middle value
ch_max
=
max
(
channel_data
)
ch_min
=
min
(
channel_data
)
ch_mid
=
(
abs
(
ch_max
)
-
abs
(
ch_min
))
ch_mean
=
mean
(
channel_data
)
print
(
"Channel max=
%
d, min=
%
d, mid=
%
d, mean=
%
d"
%
(
ch_max
,
ch_min
,
ch_mid
,
ch_mean
))
ch_data_2
=
channel_data
# Check for error
if
ch_mid
!=
0
:
raise
PtsError
(
"Unexpected middle value
%
d (must be 0)"
%
ch_mid
)
if
ch_mean
!=
-
32767
:
raise
PtsError
(
"Unexpected mean value
%
d (must be -32767)"
%
ch_mean
)
# Plot channel
#plot_channel(channel_data, ch_mean, 32768)
##################################################
# Apply ADC offset correction
##################################################
offset_corr
=
0xFF8F
print
(
"
\n
Set ADC offset correction to: 0x
%08
X"
%
(
offset_corr
))
fmc
.
set_adc_gain_offset_corr
(
1
,
0x8000
,
offset_corr
)
##################################################
# Acquire channel 1
##################################################
print
"Acquiring channel 1"
# Configure analogue input
fmc
.
set_input_range
(
1
,
IN_RANGE
)
fmc
.
set_input_term
(
1
,
IN_TERM
)
time
.
sleep
(
SSR_SET_SLEEP
)
# Perform an acquisition
acq_data
,
trig_timetag
=
acq_channels
(
fmc
,
carrier
,
ACQ_PAUSE
)
channel_data
=
acq_data
[
0
::
4
]
#print("Number of samples: %d"%(len(channel_data)))
# Calculate middle value
ch_max
=
max
(
channel_data
)
ch_min
=
min
(
channel_data
)
ch_mid
=
(
abs
(
ch_max
)
-
abs
(
ch_min
))
ch_mean
=
mean
(
channel_data
)
print
(
"Channel max=
%
d, min=
%
d, mid=
%
d, mean=
%
d"
%
(
ch_max
,
ch_min
,
ch_mid
,
ch_mean
))
ch_data_3
=
channel_data
# Check for error
if
ch_mid
!=
0
:
raise
PtsError
(
"Unexpected middle value
%
d (must be 0)"
%
ch_mid
)
if
ch_mean
!=
-
32767
:
raise
PtsError
(
"Unexpected mean value
%
d (must be -32767)"
%
ch_mean
)
# Plot channel
#plot_channel(channel_data, ch_mean, 32768)
##################################################
# Set saturation
##################################################
sat_thres
=
10000
print
(
"
\n
Set saturation to: 0x
%08
X"
%
(
sat_thres
))
fmc
.
fmc_adc_csr
.
set_field
(
'CH1_SAT'
,
'VAL'
,
sat_thres
)
fmc
.
fmc_adc_csr
.
set_field
(
'CH2_SAT'
,
'VAL'
,
sat_thres
)
fmc
.
fmc_adc_csr
.
set_field
(
'CH3_SAT'
,
'VAL'
,
sat_thres
)
fmc
.
fmc_adc_csr
.
set_field
(
'CH4_SAT'
,
'VAL'
,
sat_thres
)
##################################################
# Acquire channel 1
##################################################
print
"Acquiring channel 1"
# Perform an acquisition
acq_data
,
trig_timetag
=
acq_channels
(
fmc
,
carrier
,
ACQ_PAUSE
)
channel_data
=
acq_data
[
0
::
4
]
#print("Number of samples: %d"%(len(channel_data)))
# Calculate middle value
ch_max
=
max
(
channel_data
)
ch_min
=
min
(
channel_data
)
ch_mid
=
(
abs
(
ch_max
)
-
abs
(
ch_min
))
ch_mean
=
mean
(
channel_data
)
print
(
"Channel max=
%
d, min=
%
d, mid=
%
d, mean=
%
d"
%
(
ch_max
,
ch_min
,
ch_mid
,
ch_mean
))
ch_data_4
=
channel_data
# Check for error
if
ch_mid
!=
0
:
raise
PtsError
(
"Unexpected middle value
%
d (must be 0)"
%
ch_mid
)
if
ch_mean
!=
-
sat_thres
:
raise
PtsError
(
"Unexpected mean value
%
d (must be
%
d)"
%
(
ch_mean
,
-
sat_thres
))
##################################################
# Plot channel
##################################################
plot_channel
(
ch_data_1
,
ch_data_2
,
ch_data_3
,
ch_data_4
,
sat_thres
,
32768
)
# Make sure all switches are OFF
open_all_channels
(
fmc
)
# Check if an error occured during frequency response test
# if(error != 0):
# raise PtsError('An error occured, check log for details.')
except
(
FmcAdc100mSpecOperationError
,
FmcAdc100mOperationError
,
CalibrBoxOperationError
)
as
e
:
raise
PtsError
(
"Test failed:
%
s"
%
e
)
print
""
print
"==> End of test
%02
d"
%
TEST_NB
print
"================================================================================"
end_test_time
=
time
.
time
()
print
"Test
%02
d elapsed time:
%.2
f seconds
\n
"
%
(
TEST_NB
,
end_test_time
-
start_test_time
)
if
__name__
==
'__main__'
:
main
()
Write
Preview
Markdown
is supported
0%
Try again
or
attach a new file
Attach a file
Cancel
You are about to add
0
people
to the discussion. Proceed with caution.
Finish editing this message first!
Cancel
Please
register
or
sign in
to comment