Contact

Dr. Roland Beyer

Beam line scientist, radiation protection officer FWK
Nuclear Physics
roland.beyer@hzdr.de
Phone: +49 351 260 3281

Data Analysis with LEA

(old version up to July 2013)

The data stream or the list-mode data files written by MBS can easily be analysed using the GSI program LEA. LEA can directly read the GOOSY file format used by MBS.

The LEA project for the whole can be found at the server zirkonia in the following directorie:

zirkonia:/home/ntofdaq/lea/nng_online

This directorie contains the file f_anal_nng.c, which is used to define the data analysis procedures depending on the users wishes. Every change in this file, written in the programming language C, has to be accepted to the total LEA project by running the "make" command.

Note: Before one can use LEA, one has to run the script /home/ntofdaq/.profile!

One can include every standard C command to the f_anal file to manipulate the data. To accumulate histograms, one has to use the commands HISACC1I(1DHIST,X,VALUE) and HISACC2I(2DHIST,X,Y,VALUE). These commands will increase the content of channel X of a one dimensional histogram 1DHIST or channel (X,Y) of a two dimensional histogram 2DHIST, respectively,  by VALUE. The default value of VALUE is 1. Only in special cases, e.g., if you want to do stability investigations or for the count rate in dependence of the measurement time, one should use other values.

The names and dimensions of the histograms are defined in the files hiscre_nng.scom, hisdef_nng.h and hisini_nng.h. To modify these settings one has to run the following command:

lea -mod nng

This will open a window, where all defined histograms are listed. The use of this program is self-explaining. Note that every entry has to be confirmed by pressing the ENTER-key and changed histogram has to be confirmed be clicking the "Apply" button. To save and leave the modifying program one has to click the buttons "Apply", "definition file", "Save", and "Quit" in this order.

To run the analysis one has to do the following steps (Note: every entry in LEA has to be confirmed by pressing the ENTER-key):

0. login to zirkonia and prepare for using LEA ssh -X ntofdaq@zirkonia
password: ***********
source .profile
1. goto working directory cd lea/nng_online
2. start LEA graphical interface lea -pro nng -men
3. choose data source if lmd file:
if MBS stream:
Source = LMD File; Input file/node = /path/filename.lmd
Source = MBS samples; Input file/node = rio42
4. start analysis click "Analyze events"

The following histograms are produced by LEA:

Name xmin xmax xbinsize ymin ymax ybinsize Meaning (x) Meaning (y) Meaning (z)
ADC 0 4096 1       FC ADC raw value counts  
BAFDIFFTDC 0 2000 1 0 18 1 Time difference between PMTs / TDC ch
0...15 : BaF2 1..16
counts
BAFDIFFTDC2 0 2000 1 0 18 1 Time difference between PMTs including accelerator / TDC ch
0...15 : BaF2 1..16
counts
BAFLG 0 4096 1 0 36 1 Long Gate raw value
0...31 : BaF2 DAQch1...32
counts
BAFLGS 0 4096 1 0 36 1 Long Gate Sensitive raw value
0...31 : BaF2 DAQch1...32
counts
BAFSG 0 4096 1 0 36 1 Short Gate raw value
0...31 : BaF2 DAQch1...32
counts
BAFSGS 0 4096 1 0 36 1 Short Gate Sensitive raw value
0...31 : BaF2 DAQch1...32
counts
BAFLGSUM 0 4096 1 0 18 1 Harmonic mean of PMT LG values
0...15 : BaF2 1...16
counts
BAFLGSSUM 0 4096 1 0 18 1 Harmonic mean of PMT LGS values
0...15 : BaF2 1...16
counts
BAFSGSUM 0 4096 1 0 18 1 Geometric mean of PMT SG values
0...15 : BaF2 1...16
counts
BAFSUMTDCS -8192 8192 1 0 18 1 Time sum of PMTs / 10 TDC ch
0...15 : BaF2 1...16
counts
BAFSUMTDCS2 -8192 8192 1 0 18 1 same as BAFSUMTDCS but with gate on BAFDIFFTDC
0...15 : BaF2 1...16
counts
HITBAF 0 32 1 0 36 1 No. of hits per Channel
0 : accelerator
1...32 : BaF2 DAQ ch1...32
33 : trigger
counts
HITFC 0 32 1       No. of hits per Channel counts  
HITPL 0 32 1 0 16 1 No. of hits per Channel
0 : accelerator
1...12 : plastic DAQch1...10
13 : trigger
14 : fission chamber
counts
OPC 0 8192 1 0 32 1 OPC data OPC values (see Data structure)  
PLDIFFTDC 0 2000 1 0 8 1 Time difference between PMTs / TDC ch
0...5 : Plastic 1...6
 counts
PLDIFFTDC2 0 2000 1 0 8 1 Time difference between PMTs including accelerator / TDC ch
0...5 : Plastic 1...6
 counts
PLLG 0 4096 1 0 16 1 Long Gate raw value
0...11 : Pl DAQch1...12
counts
PLLGS 0 4096 1 0 16 1 Long Gate Sensitive raw value
0...11 : Pl DAQch1...12
counts
PLSG 0 4096 1 0 16 1 Short Gate raw value
0...11 : Pl DAQch1...12
counts
PLSGS 0 4096 1 0 16 1 Short Gate Sensitive raw value
0...11 : Pl DAQch1...12
counts
PLSGSUM 0 4096 1 0 8 1 Harmonic mean of PMT SG values
0...5 : Plastic 1...6
counts
PLSUMTDCS -8192 8192 1 0 8 1 Time sum of PMTs / 10 TDC ch
0...5 : Plastic 1...6
counts
PLSUMTDCS2 -8192 8192 1 0 8 1 same as PLSUMTDCS but with gate on PLDIFFTDC
0...5 : Plastic 1...6
counts
PLSUMTDCS21 -8192 8192 1 0 8 1 same as PLSUMTDCS2 but with target 1
0...5 : Plastic 1...6
counts
PLSUMTDCS23 -8192 8192 1 0 8 1 same as PLSUMTDCS2 but with target 3
0...5 : Plastic 1...6
counts
RATE1...3 0 8192 1 0 32 1 number of scaler readouts input channel of scaler 1...3 count rate = counts divided by time flag
TDCBAF 0 81920 10 0 36 1 TDC raw value
0 : accelerator
1...32 : BaF2 DAQch1...32
33 : trigger
counts
TDCBAFACC 0 81920 10 0 36 1 TDC raw value of accelerator used for TDCBAFDIFF
1...32 : BaF2 DAQch1...32
counts
TDCBAFDIFF -8192 8192 2 0 36 1 Time difference to accelerator / 10 TDC ch
1...32 : BaF2 DAQch1...32
counts
TDCBAFH0 0 81920 10 0 36 1 TDC raw value of first hit in matchwindow
0 : accelerator
1...32 : BaF2 DAQch1...32
counts
TDCFCACC 0 81920 10       TDC raw value of accelerator used for TDCFCDIFF counts  
TDCFCDIFF -8192 8192 2       Fission chamber time difference to accelerator / 10 TDC ch counts  
TDCFCH0 0 81920 10       Fission chamber TDC raw value of first hit in matchwindow counts  
TDCPL 0 81920 10 0 12 1 TDC raw value
0 : accelerator
1...12 : Pl DAQch1...12
13 : trigger
14 : fission chamber
15 : veto end
counts
TDCPLACC 0 81920 10 0 16 1 TDC raw value of accelerator used for TDCPLDIFF
1...12 : Pl DAQch1...12
counts
TDCPLDIFF -8192 8192 2 0 16 1 time difference to accelerator / 10 TDC ch
1...12 : Pl DAQch1...12
counts
TDCPLH0 0 81920 10 0 16 1 TDC raw value of first hit in matchwindow
0 : accelerator
1...12 : Pl DAQch1...12
counts
TRGTIME 0 1048580 4       Time difference to previous trigger time tag / 800 ns  counts  
TWODMULT 0 16 1 0 8 1 BaF2 that fired Plastic that triggered counts
TWODLGSGSUM01...16 0 4096 16 0 4096 16 LGSUM value of BaF2 1..16 SGSUM value of BaF2 1...16 counts
VETO 0 653600 25       Global dead time / ns counts  
VETOLONG 0 655360 100 0 2 1 Readout dead time length / ns
0 : Plastic TDT
1 : BaF2 TDT
counts
VETOSHORT 0 25600 25 0 16 1 Busy length / ns
0...7 : BaF2 TAPS busy
8...10 : Pl TAPS busy
11 : ADC busy
12 : Trigger
counts