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 |
|
counts |
| BAFDIFFTDC2 |
0 |
2000 |
1 |
0 |
18 |
1 |
Time difference between PMTs including accelerator / TDC ch |
|
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 |
|
counts |
| BAFLGSSUM |
0 |
4096 |
1 |
0 |
18 |
1 |
Harmonic mean of PMT LGS values |
|
counts |
| BAFSGSUM |
0 |
4096 |
1 |
0 |
18 |
1 |
Geometric mean of PMT SG values |
|
counts |
| BAFSUMTDCS |
-8192 |
8192 |
1 |
0 |
18 |
1 |
Time sum of PMTs / 10 TDC ch |
|
counts |
| BAFSUMTDCS2 |
-8192 |
8192 |
1 |
0 |
18 |
1 |
same as BAFSUMTDCS but with gate on BAFDIFFTDC |
|
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 |
|
counts |
| PLDIFFTDC2 |
0 |
2000 |
1 |
0 |
8 |
1 |
Time difference between PMTs including accelerator / TDC ch |
|
counts |
| PLLG |
0 |
4096 |
1 |
0 |
16 |
1 |
Long Gate raw value |
|
counts |
| PLLGS |
0 |
4096 |
1 |
0 |
16 |
1 |
Long Gate Sensitive raw value |
|
counts |
| PLSG |
0 |
4096 |
1 |
0 |
16 |
1 |
Short Gate raw value |
|
counts |
| PLSGS |
0 |
4096 |
1 |
0 |
16 |
1 |
Short Gate Sensitive raw value |
|
counts |
| PLSGSUM |
0 |
4096 |
1 |
0 |
8 |
1 |
Harmonic mean of PMT SG values |
|
counts |
| PLSUMTDCS |
-8192 |
8192 |
1 |
0 |
8 |
1 |
Time sum of PMTs / 10 TDC ch |
|
counts |
| PLSUMTDCS2 |
-8192 |
8192 |
1 |
0 |
8 |
1 |
same as PLSUMTDCS but with gate on PLDIFFTDC |
|
counts |
| PLSUMTDCS21 |
-8192 |
8192 |
1 |
0 |
8 |
1 |
same as PLSUMTDCS2 but with target 1 |
|
counts |
| PLSUMTDCS23 |
-8192 |
8192 |
1 |
0 |
8 |
1 |
same as PLSUMTDCS2 but with target 3 |
|
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 |
|
counts |
| TDCPLDIFF |
-8192 |
8192 |
2 |
0 |
16 |
1 |
time difference to accelerator / 10 TDC ch |
|
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 |
