logo

A real-time online data acquisition system for Dragon-I linear induction accelerator

NUCLEAR ELECTRONICS AND INSTRUMENTATION

A real-time online data acquisition system for Dragon-I linear induction accelerator

Jun Wu
Xiang Zhou
Chen Yuan
Ze-Jie Yin
Nuclear Science and TechniquesVol.28, No.3Article number 37Published in print 01 Mar 2017Available online 03 Feb 2017
41900

An extensible high-speed accelerator data acquisition system (ADAS) for the Dragon-I linear induction accelerator (LIA) has been developed. It comprises a VXI crate, a controller, four data acquisition plugins, and a host computer. A digital compensation algorithm is used to compensate for the distortion of high-speed signals arising from long-distance transmission. Compared with the traditional oscilloscope wall, ADAS has significant improvements in system integration, automation and reliability. It achieves unified management of data acquisition and waveform monitoring, and performs excellently with a 107-ps high-accuracy trigger and 32-channel signal monitoring. In this paper, we focus on the system architecture and hardware design of the ADAS, and realization of the trigger and digital compensation algorithm.

ADASLINACBeam position monitorFPGA

1 Introduction

A linear induction accelerator (LIA) accelerates charged particles by utilizing the induced electromotive force produced by magnetic flux change. It has been widely applied in areas of radiography, free electron lasers, high-power microwaves, and heavy ion fusion [1, 2]. The Dragon-I, a new LIA constructed at the Institute of Fluid Physics, Chinese Academy of Engineering Physics, is able to provide 20 MeV2.5 kA electron beams with a pulse width of 70 ns (FWHM) [3, 4].

The accelerator data acquisition system (ADAS), as a key component of Dragon-I, receives output signals from a detector subsystem to provide diagnosis for working conditions of the accelerator and relevant information about beam quality. The detectors include resistive wall, B-Dot and D-Dot [5-7]. The detector unit is shown schematically in Fig. 1, where four or eight detectors, with equal distance and the beam centered, are located around the vacuum chamber [8, 9].

Fig. 1
Schematic showing the position of the equally spaced detectors
pic

Under the circumstances of four detectors (A–D), the average voltage induced is determined by Eq. (1), where, d is the electrode diameter; R is the load resistance; I¯c is the average current; a is the vacuum chamber radius; e is the elementary charge; σz is the bunch length; f is the repetition frequency; (x, y) is the beam position; and δ, θ and F(δ, θ) are defined in Eqs. (2), (3), and (4), respectively. [10].

VA,B,C,Dd2RI¯cF(δ,θA,B,C,D)8a2πeσz2f (1) δ=x2+y2 (2) θA,B,C,D=nπ4tan-1(yx) (3) F(δ,θ)=a2δ2a2+δ22aδcosθ (4)

Excessive number of parameters would lead to complexity in calculation. Thus, we propose a simplified calculation model of normalized voltages. Eqs. (5) and (6) are the definitions of the normalized electrode signals.

U=[(VB+VC)(VA+VD)]/(VA+VB+VC+VD) (5) V=[(VA+VB)(VC+VD)]/(VA+VB+VC+VD) (6)

With U and V, Eqs. (7) and (8) are expressions for the value of beam position (x, y).

x=i=0kj=0iAi-j,jUijVj (7) y=i=0kj=0iBi-j,jUijVj (8)

Dragon-I produces high-speed single pulse beam, and the ADAS shall be able to measure the high-speed signals. In practical operation, long-distance transmission usually causes analog signal distortion, which affects negatively on beam monitoring. In multi-channel operations, the electronic systems need a synchronous trigger with high accuracy for channel alignment. To solve these problems, the ADAS shall possess functions of multi-channel high-speed sampling, high accuracy trigger, and customizable data processing. The key requirements of ADAS are as follows:

• High integration, 32 acquisition channels;

• 1 Gsps sampling rates, signal bandwidth ≥ 100 MHz;

• ENOB ≥ 7 bit @ 100 MHz;

• Digital compensation function;

• Low DC drift, low gain error, self-calibration;

• Multi-channel trigger, trigger accuracy, i.e., root mean square (RMS) < 500 ps;

• Automatic operation, automatic storage, and real-time data processing.

2 Realization of the ADAS

2.1 System structure

Figure 2 shows the architecture of ADAS [11]. It consists of four data acquisition plugins, one VXI crate with its controller, and one host computer. The plugins are responsible for data acquisition and processing. Each plugin has eight acquisition channels that offer 8-bit 1 Gsps sampling. The VXI crate provides sufficient space and power supply for over 32 channels. The host computer is responsible for controlling, data receiving and displaying for the entire system, and an IEEE 1394 cable is applied to transmit data between the VXI crate and the host computer.

Fig. 2
Architecture of the ADAS.
pic
2.2 Hardware design

The VXI-based instrument is primarily in charge of signal acquisition and processing [12]. Fig. 3 shows architecture of the plugin. It consists of eight analog-front-ends (AFEs), four analog-to-digital converters (ADCs), two field-programmable gate arrays (FPGAs), the power supply unit, and a clock generator. The function of the VXI bus is to transmit data from the plugin to the host.

Fig. 3
Architecture of the plugin
pic

The AFE focuses on receiving and conditioning the input signals [13]. Since the detectors are located in high-frequency electromagnetic interference (EMI) environment, a five-step LC anti-alias filter (AAF) is needed in front of the ADC. It can block high-frequency noise and limit the pass band at DC-125 MHz to increase signal-to-noise ratio (SNR).

Signals from Dragon-I detectors are of variable dynamic ranges, so the AFE gain should be adjustable to match dynamic range of the ADC input. Based on resistance adjustment, traditional method of gain control fails to meet the requirements of high integration and automation, and a fixed-gain amplifier cooperating with a digital step attenuator (DSA) is used for digital control. The gain parameters can be adjusted by regulating attenuation ratio of the DSA through VXI bus. Taking advantage of the high performance DSA, the AFE gain ranges from −21.5 dB to +10 dB in 0.25 dB steps.

Calibration of the AFE is always necessary, since the AFE would introduce gain error and DC offset, hence the decrease of accuracy. Automatic AFE calibration based on FPGA and ADC is used, rather than traditional manual calibration. As the input bias and the ADC gain are adjustable, the AFE input is first connected to the ground so that the AFE input bias can be removed through adjusting the ADC input bias. Next, a stable reference voltage is used as the AFE input and by regulating the ADC gain, and the AFE gain can be modified. The self-calibration process carries on automatically right after the system is powered on. Meanwhile, manual recalibration is available if necessary.

As function core of the VXI plugin, the two FPGAs are responsible for data receiving and processing, and control of slave devices as well [14, 15]. Fig. 4 shows the architecture of inner logic of the FPGAs. The serializer/deserializer (serdes) and data acquisition units receive data from the ADCs, while the data processing unit accomplishes signal analysis. Other units are responsible for system communication, coordinated control, and operation support.

Fig. 4
Architecture of FPGA logic.
pic
2.3 Data acquisition and processing

To achieve high SNR sampling, a clock with extremely low aperture jitter (tj) is needed. Eq.(9) is the definition of SNR [16].

SNR=20lg[12πftj] (9)

As the jitter of sampling clock of this system must be less than 5 ps, the clock chip of LMK04821 is adopted to provide a 1 GHz synchronous sampling clock for ADCs with clock output jitter less than 300 fs.

ADC08D1000, an 8-bit 1 Gsps ADC is used for analog signal sampling. It has a 1: 2 demultiplexer that deserializes the 8-bit 1 GHz data stream into a 16-bit 500 MHz data stream, exporting the data stream via a low-voltage differential signaling (LVDS) bus. In the FPGA, the serdes unit deserializes the 16-bit 500 MHz high-speed data stream into 64-bit 125 Mbps low-speed data stream, which is then written into ring storage, where a section of the latest data can be captured and stored when triggered. Next, the processing module performs real-time analysis to obtain characteristic parameters of the signals, including amplitude peak, pulse width, rising and falling time, and average voltage of a specified region. Finally, the plugins upload both the raw data and the analysis results to the host.

The ADAS is located considerably far away from the detector, connected by tens or hundreds of meters of coaxial cables, making signal distortion a common issue during high-speed pulse transmission. Therefore, digital compensation in FPGA achieves raw signal reconstruction.

In general, skin effects of cables under high frequency causes signal attenuation and the insulating medium causes high-frequency loss, as described in Eqs. (10) and (11), respectively, where f is signal frequency, E and F are coefficients determined by the cable material and structure, d is the diameter of the conductor, ε is the relative dielectric constant, D is the insulation diameter, and δ is the dielectric loss angle [17].

ar=2.6×103ε(Ed+FD)1lg(D/d)f (10) aG=9.1×105fεtanδ (11)

Considering the influence of DC resistance, the attenuation expression a(f) can be simplified as Eq. (12).

a(f)=C+Bf+Af (12)

Eq. (12) indicates the relationship between signal frequency and signal attenuation ratio through cables. With a digital sweeper, the frequency spectrum of the signals through each coaxial cable can be measured, and the coefficients of A, B and C can be calculated using least square method (LSM). When capturing the valid signal, the data processing unit obtains its discrete frequency spectrum using fast Fourier transform (FFT), and then modifies the spectrum using Eq. (12). Finally, the unit reconstructs raw signal through the inverse fast Fourier transform.

2.4 Realization of the trigger unit

The trigger unit determines the initial time of signal acquisition, playing a key role in data alignment. As the external trigger signal arrives randomly, a 16 ns jitter can occur between any two channels since the inner clock of the FPGA is up to 125 MHz. To reduce jitter between channels, high accuracy trigger position should be determined.

For this purpose, two time-to-digital converters (TDCs) [18-21] are used to obtain time measuring capability using the FPGA’s logic resources. As excessively short measuring time leads to low accuracy, the TDC measures the time interval between rising edges of the trigger and the next third rising edges of data clock (DCLK) to locate trigger position. However, the TDC will miss recording the rising edges if the time interval between the rising edges of the trigger and the nearest DCLK is less than 110 ps, leading to inaccurate calculation. Therefore, once receiving a trigger, the FPGA shall produce a 1250 ps delayed signal called Trigger_d through the IO delay primitive. The correct measured time is then calculated by analyzing the results of the measurements, t and td, as illustrated in Fig. 5.

Fig. 5
Diagram illustrating triggered time measurements t and td.
pic

There are three cases worth considering, each of which is described below.

(1) If ttd is approximately 1250 ps, neither of the TDCs miss recording a rising edge; t is the measured time T.

(2) If ttd is approximately −2750 ps, TDC_A does not miss recording a rising edge; t is the measured time T.

(3) If ttd is approximately 5250 ps, TDC_B does not miss recording a rising edge; td +1250 ps is the measured time T.

Suppose the ADC sampling rate is 1 Gsps, the ADC’s output data are highly synchronous with the DCLK. Thus, suppose that the pre-trigger number is N, the calculated trigger position is (N+12−T) ns. When a plugin receives a trigger signal, it simultaneously broadcasts the signal to the rest of the plugins via the VXI backplane with an emitter-coupled logic (ECL) signal. Since the broadcast signal and clock synchronization pulse use the same transmission path, the time delay caused by different transmission distances cancels each other out. Through this method, the trigger position in the dataflow of every plugin is definite, thus realizing high-accuracy trigger across plugins.

2.5 Software design

The software is in charge of system initialization, data receiving and waveform displaying. Fig. 6 is its workflow. In system initialization, the software sends operation parameters according to the detector’s type and working conditions. It also specifies a plugin that is responsible for external trigger signal receiving. Next, the software commands the specific plugin to send a clock synchronization pulse via the ECL line to synchronize other plugins. On receiving an interrupt, the system reads data from the plugin and stores them on the computer’s hard drive. The host displays the signal waveforms on the monitor. In addition, the software provides an API for other programs to obtain data for further processing.

Fig. 6
Software workflow of our system
pic

3 Test results

3.1 System performance test

For bandwidth test, we used a Fluke 9500B oscilloscope calibration workstation to generate fixed-amplitude sine waves with variable frequencies as input signal. Results showed the bandwidth of the 32 channels ranging from 117 to 132 MHz.

For gain error test, we used an AFG3102C arbitrary function generator to generate linear sweep signal Vs (−4.75 to +4.75 V, 100 kHz). Then the plugin measured the voltage and calculated the gain error by linear fit the measured data. The results indicated that the maximum gain error of each channel was less than 0.75%.

For ENOB performance test, we generated three sine waves of different frequencies using the Fluke 9500B as analog inputs (9.971 MHz with a 20-MHz low-pass filter, 59.971 MHz with a 100-MHz low-pass filter, and 99.971 MHz with a 200-MHz low-pass filter separately; 9.5 V in amplitude). The ADAS sampled the analog signals and calculated the ENOB. Test results showed that ENOB was always greater than 7.05 bits under the three conditions.

Finally, to test accuracy of the trigger, we used a pulse signal synchronized with the DCLK as trigger input, with repetitive measurement over 1000 times. Test results showed that the RMS accuracy of the trigger was 107 ps.

3.2 Field test

We operated the ADAS on the Dragon-I monitor system. Fig. 7 shows the reconstruction effect report. The power splitter divided the detector output into two signals: the green one, which was sent to the ADAS through a short cable; and the red one, which was transmitted through an 85 m cable; with its reconstructed waveform being shown in blue. Analysis results showed great consistency between the reconstructed signal and the original one, with flat-topped average amplitude error of <1.6% and pulse width error of < 2 ns.

Fig. 7
Results of the waveform reconstruction test
pic

Fig. 8 is a screenshot of the ADAS software during the field test, reflecting the V-t transients of eight types of detected signals with a time base of 50 ns. Among all of the sampled signals, the fastest signal rising edge was 7 ns (see the enlarged section), suggesting that the system is competent in precise signal acquisition with a signal rising edge less than 10 ns.

Fig. 8
The ADAS-detected signals in the field test on Dragon-I.
pic

4 Conclusion

We have developed the extensible high-speed ADAS. Owing to the calibration of the AFE and TDC-based trigger module, the ENOB of the plugin is greater than 7.05 bits and the RMS of the trigger accuracy is better than 107 ps. The data acquisition, control, and trigger algorithms are realized in the FPGAs, by which real-time data acquisition is guaranteed. During long-term operation of the Dragon-I accelerator, performance of the ADAS indicates that it has achieved the anticipated goals of high integration, high automation, and high reliability in data acquisition. The ADAS has potential applications on other large physical experiment facilities, without changing the circuit structure.

References
1. J. E. Leiss,

Induction Linear Accelerators and Their Applications

. IEEE Trans. Nucl. Sci. 3, 3869-3876(1979). doi: 10.1109/TNS.1979.4330636
Baidu ScholarGoogle Scholar
2. J. Cao, R. Xiao, N. Chen, et al.,

Spectral reconstruction of the flash X-ray generated by Dragon-I LIA based on transmission measurements

. Nucl. Sci. Tech. 26, 040403(2015). doi: 10.13538/j.1001-8042/nst.26.040403
Baidu ScholarGoogle Scholar
3. H. Yu, J. Zhu, N. Chen, et al.,

Measurements and effects of backstreaming ions produced at bremsstrahlung converter target in Dragon-I linear induction accelerator

. Rev. Sci. Instrum. 81, 043302(2010). doi: 10.1063/1.3379132
Baidu ScholarGoogle Scholar
4. K.Z. Zhang, L. Wen, H. Li, et al.,

Dragon-I injector based on the induction voltage adder technique

. Phys. Rev. Spec. TOP-AC 9, 080401(2006). doi: 10.1103/PhysRevSTAB.9.080401
Baidu ScholarGoogle Scholar
5. X. He, J. Pang, Q. Li, et al.,

Simulations of axial B-dot monitors inside a groove in the beam pipe

. Nucl. Instrum. Meth. A 729, 14-18(2013). doi: 10.1016/j.nima.2013.06.092
Baidu ScholarGoogle Scholar
6. Q. Luo, B.G. Sun, Q.K. Jia, et al.,

Design and cold test of S-BAND cavity BPM for HLS

. Sci. China Phys. Mech. 54, 292-295(2011). doi: 10.1007/s11433-011-4555-y
Baidu ScholarGoogle Scholar
7. B.E. Carlsten,

Deconvolving long-pulse electron-beam centroid data from resistive wall-current monitors

. Rev. Sci. Instrum. 73, 3772-3777(2002). doi: 10.1063/1.1515763
Baidu ScholarGoogle Scholar
8. L. Qin, L. Hong, C. Nan, et al.,

B-dot monitor for intense electron beam measurement

. High Power Laser and Particle Beams 9, 025(2009).(in Chinese)
Baidu ScholarGoogle Scholar
9. L. Zhao, X. Gao, X. Hu, et al.,

Beam Position and Phase Measurement System for the Proton Accelerator in ADS

. IEEE Trans. Nucl. Sci. 61, 538-545(2014). doi: 10.1109/TNS.2013.2291779
Baidu ScholarGoogle Scholar
10. T. Shintake, M. Tejima, H. Ishii, et al.,

Sensitivity calculation of beam position monitor using boundary element method

. Nucl. Instrum. Meth. A 254, 146-150(1987). doi: 10.1016/0168-9002(87)90496-7
Baidu ScholarGoogle Scholar
11. J. Gál, G. Hegyesi, J. Molnár, et al.,

The VXI electronics of the DIAMANT particle detector array

. Nucl. Instrum. Meth. A 516, 502-510 (2004). doi: 10.1016/jnima.2003.08.158
Baidu ScholarGoogle Scholar
12. P. Wang, R.Y. Zhang, Y.Y. Yan, et al.,

An acquisition system of digital nuclear signal processing for the algorithm development

. Nucl. Sci. Tech. 24, 60408-060408(2013). doi: 10.13538/j.1001-8042/nst.2013.06.012
Baidu ScholarGoogle Scholar
13. S. Ivanenko, A. Khilchenko, E. Puryga, et al.,

Prototype of Fast Data Acquisition Systems for ITER Divertor Thomson Scattering Diagnostic

. IEEE Trans. Nucl. Sci. 62, 1181-1186(2015). doi: 10.1109/TNS.2015.2428195
Baidu ScholarGoogle Scholar
14. S.P. Li, X.F. Xu, H.R. Cao, et al.,

A real-time n/γ digital pulse shape discriminator based on FPGA

. Appl. Radiat. Isotopes 72, 30-34(2013). doi: 10.1016/j.apradiso.2012.10.004
Baidu ScholarGoogle Scholar
15. J.H. Zhang, Y.Y. Wang, G.Y. Nan, et al.,

Trigger signal pre-processing in nuclear physics experiment data acquisition systems

. High Power Laser and Particle Beams 24, 2727-2730 (2012). doi: 10.3788/HPLPB20122411.2727
Baidu ScholarGoogle Scholar
16. M. Shinagawa, Y. Akazawa, T. Wakimoto.

Jitter analysis of high-speed sampling systems

. IEEE J. SOLID-ST. Circ. 25, 220-224(1990). doi: 10.1109/4.50307
Baidu ScholarGoogle Scholar
17. Y. Zhang.

Factors affecting attenuation of RF coaxial cables

. Electric wire and cable 6, 13-15(2008).(in Chinese)
Baidu ScholarGoogle Scholar
18. J. Song, Q. An, S.B. Liu.

High-Resolution Time-to-Digital Converter Implemented in Field-Programmable-Gate-Arrays

. IEEE Trans. Nucl. Sci. 53, 236-241(2006). doi: 10.1109/TNS.2006.869820
Baidu ScholarGoogle Scholar
19. X.J. Hao, S.B. Liu, L. Zhao, et al.,

A digitalizing board for the prototype array of LHAASO WCDA

. Nucl. Sci. Tech. 22, 178-184(2011). doi: 10.13538/j.1001-8042/nst.22.178-184
Baidu ScholarGoogle Scholar
20. L. Dong, J.F. Yang, K.Z. Song.

Carry-chain propagation delay impacts on resolution of FPGA-based TDC

. Nucl. Sci. Tech. 25, 030401-030401 (2014). doi: 10.13538/j.1001-8042/nst.25.030401
Baidu ScholarGoogle Scholar
21. T. Suwada, F. Miyahara, K. Furukawa, et al.,

Wide dynamic range FPGA-based TDC for monitoring a trigger timing distribution system in linear accelerators

. Nucl. Instrum. Meth. A 786, 83-90(2015). doi: 10.1016/j.nima.2015.03.019
Baidu ScholarGoogle Scholar