Start 00:43 22.2.25 https://elog.gsi.de/despec/G-24-00302/57
AIDA data file R9_85
Converting hexadecimal timestamp to decimal: 1740181266259754800
Assuming that this timestamp is in nanoseconds (1 billionth of a second):
GMT: Friday, February 21, 2025 11:41:06.259 PM
Your time zone: Friday, February 21, 2025 11:41:06.259 PM GMT+00:00
End 12:56 22.2.25 https://elog.gsi.de/despec/G-24-00302/85 (start of next degrader setting)
AIDA data file R9_199
Converting hexadecimal timestamp to decimal: 1740224934731157800
Assuming that this timestamp is in nanoseconds (1 billionth of a second):
GMT: Saturday, February 22, 2025 11:48:54.731 AM
Your time zone: Saturday, February 22, 2025 11:48:54.731 AM GMT+00:00
Analysis data files FEB25/R9_85 and FEB25/R9_199
- no timewarps
- deadtime aida05 c. 37% and 39% due to high SC41 scaler input to MACB - issue subsequently fixed
- deadtime aida04 c. 4%, all other FEE64s <<1%
FEE64 configuration
FEE64 a b c
g h
d e f
a b c d e f g h
DSSSD#1 15 3 12 9 1 5 2 4
n+n Ohmic FEE64s 2, 4
Data analysis assumes
- all LEC ADC data channels with valid ADC offset included (507 of 512 channels)
LEC calibration ADC offset only
- no clustering
- no multiplex timestamp correction
- no p+n junction side - n+n Ohmic side correlation time gates
- FEE64 *not* DSSSD strip ordering
- hardware - slow comparator setting p+n junction FEE64s 100keV, n+n Ohmic FEE64s 150keV
- LEC energy difference +/- 168keV
- HEC energy difference +/- 1.68GeV
- valid LEC events
p+n junction side multiplicity = 1 and n+n Ohmic side multiplicity = 1
151keV < LEC energy < 1000keV
to select candidate beta events
standalone analysis of AIDA data, no downstream veto detector
- valid HEC events
p+n junction side multiplicity > 0 and n+n Ohmic side multiplicity > 0
(x,y) strips corresponding to maximum energy
p+n junction and n+n Ohmic side HEC
- HEC veto
not available - only 1x AIDA DSSSD installed
- per pixel implant-decay correlations
- end of event
difference in WR timestamp between successive ADC data items > 2500 and overall event length < 33us
per FEE64 LEC ADC data items 268ms/channel - attachments 3-5
- all
- 150keV < energy < 1500keV
- energy > 1500keV
absence of high instantaneous rates on spill cf. S181, S100, S505 etc
spill structure visible for energies > 1500keV - probably to be expected as majority of such events should be light ions
per FEE64 HEC ADC data items 268ms/channel - attachments 6-8
- all
- 100MeV < energy < 1000MeV
- energy > 1000MeV
Implant & decay event rates 262us/channel - attachment 9
LEC m_p versus m_n - attachment 10
LEC e_p versus e_n - 20keV/channel - attachment 11
HEC x strip versus y strip - attachment 12
- HEC-LEC implant decay time difference <1s
- HEC-LEC implant decay time difference <100s
shows x-y window used to select 82Nb events
HEC m_p versus m_n - attachment 13
HEC e_p versus e_n - 20MeV/channel - attachment 14
HEC-LEC implant-decay time difference (4ms/channel) versus - attachment 15
- HEC energy
- LEC energy
- x & y strip
x strips 130 & 136 disabled
HEC & LEC p strip - n strip time difference (2us/channel, offset=2000 channels) - attachment 16
HEC E (20MeV/channel) versus implant-decay time difference (4.194ms/channel) - attachment 17
LEC E (20keV/channel) versus implant-decay time difference (4.194ms/channel) - attachment 18
LEC e_p - e_n (5.6keV/channel) versus implant-decay time difference (4.194ms/channel) - attachment 19
HEC x & y strip versus implant-decay time difference (4.194ms/channel) - attachment 20
HEC-LEC implant-decay time difference (8.389ms/channel) for candidate 82Nb events - attachment 21
- blue - forward
- cyan - backward
Note
- no FRS PID
- no HEC dE
- no HEC veto
- no HEC energy gates - assume all events within x-y window *stop* in AIDA DSSSD - 82Nb events known to overlap in x-y with (longer-lived, more numerous) 83Nb events
- no front-back time difference
- no clustering
- *all* LEC channels with valid ADC offsets (507 of 512) included *except* 2x strips 130 & 136
HEC-LEC implant-decay time difference (8.389ms/channel) for candidate 82Nb events - attachment 22
Weighted least squares fit channels 1-12 half life 55(9)ms cf. NNDC 50.0(3)ms.
Added 8.2.26 by TD
Attachment 23 - implant timestamp, decay timestamp, decay timestamp - implant timestamp for selected events (attachment 21) with time differences < 255 x 4.194304ms
Attachment 24 - implant timestamp, decay timestamp, implant timestamp - decay timestamp for selected events (attachment 21) with time differences < 255 x 4.194304ms
Attachment 25 - analysis program
Attachment 27 - spectra titles
Attachment 28 - analysis program variables
For data files R9_85 - R9_199 inclusive
Total ADC data items 30117120061
Total ADC events 6551938947
Questions from Giorgio Bruni Campanella Spezza
1) What were the exact number of implants, decays, correlated implants and decays and correlation efficiency?
Implants ( m_p_hec(1) > 0 .and. m_n_hec(1) > 0 ) 1046107
Decays ( m_p_lec(1) = 1 .and. m_n_lec(1) = 1 ) .and. ( ediff.LT.30 .and. ediff.GT.-30 ) 83342477
Selected implants s2112
( 75 < x < 150 .and. 30 < y < 120 ) .and. ( ediff.LT.300 .and. ediff.GT.-300 ) .and . ( 100 < hec_e < 4096 ) 215705 = 32.0/pixel = 0.000732/pixel/s
Selected decays s2040
( 75 < x < 150 .and. 30 < y < 120 ) .and. ( ediff.LT.30 .and. ediff.GT.-30 ) .and. ( decay-implant < 1s ) .and. ( 27 < lec_e < 180 ) 1056 = 0.156/pixel = 0.00000358/pixel/s
Selected decays s2041
( 75 < x < 150 .and. 30 < y < 120 ) .and. ( ediff.LT.30 .and. ediff.GT.-30 ) .and. ( decay-implant < 100s ) .and. ( 27 < lec_e < 180 ) 38483 = 5.70/pixel = 0.000131/pixel/s
2) What is the number of correlated decays for the backward time correlations?
Correlated forward decays s2226 (channels 0-255, 4.19ms/channel) = 1422
Correlated backward decays s2224 (channels 0-255, 4.19ms/channel) = 778
3) How exactly do you define beta detection efficiency for the decay of a specific isotope? I ask this due to the lack of a PID in the analysis, and hence, without any FRS-AIDA implant
correlation, how can you determine the number of implanted ions of a specific isotope?
Yes - to calculate beta efficiency for 82Nb (say) I need an estimate of what fraction of selected implant events are actually 82Nb
4) Do you have any vetoes for the position of implants in AIDA, or do you accept all implants throughout the DSSSD?
See definitions above - there are gates for x and y strips and implant energy
Added TD 12.4.26
Attachment 28 - implant timestamp, decay timestamp, p+n strip channel, n+n strip channel, decay timestamp - implant timestamp for selected events (attachment 21) with time differences < 255 x 4.194304ms, DSSSD#1 p+n and n+n strip #
Attachment 29 - implant timestamp, decay timestamp, p+n strip channel, n+n strip channel, decay timestamp - implant timestamp for selected events (attachment 21) with time differences < 255 x 4.194304ms, DSSSD#1 p+n and n+n strip #
p+n strip channel/n+n strip channel calculated using GREAT format encoded data (channel, module, range)
channel = channel + (64 * FEE64 module #) + (2048 * range)
where range = 0 (LEC data), 1 (HEC data), module = 0-15, channel = 0-63
DSSSD#1 p+n strip # channels 0-383, n+n strip # channels 384-511 |
*trigger
128
*oned
0..15 s 65536
20..35 s 65536
40..55 s 65536
60..75 s 64
80..96 s 64
100..115 s 65536
120..135 s 65536
140..155 s 65536
200..215 s 65536
250..254 s 4096
300..331 s 65536
400..405 s 65536
1400..1403 s 65536
2130..2131 s 65536
2030..2031 s 65536
2200..2201 s 65536
2210..2211 s 65536
2220..2229 s 1024
3030..3033 s 4096
*twod
1000..1001 s 4096 4096
1100..1101 s 512 512
2000..2003 s 512 512
2010..2011 s 512 512
2020..2021 s 1024 1024
2040..2049 s 512 512
2100..2101 s 512 512
2110..2113 s 512 512
2120..2127 s 1024 1024
2300..2301 s 1024 1024
2400..2401 s 1024 1024
2500..2501 s 1024 1024
2600..2601 s 1024 1024
3000..3001 s 512 512
3010..3011 s 1024 1024
3020..3021 s 512 512
3040..3041 s 512 512
*vars
*sort
C 201..232 s 65536
C 250..254 s 4096
C 2500..2505 s 1024 1024
C----67---------------------------------------------------------------72------80
SUBROUTINE init
IMPLICIT none
SAVE
C
C External functions
C
EXTERNAL dtime
C
INTEGER and, lshift, rshift, bb18order
C
LOGICAL btest
C
REAL float, secnds
C
C Parameter variables
C
INTEGER max
PARAMETER (max = 4096)
C
C Local variables
C
INTEGER adc_data, asic, ch, chx, channel, channel_ident
INTEGER count, data(0:4095), disc
INTEGER dsssd_a( 6 ), dsssd_b( 6 ), dsssd_c( 6 ), dsssd_d( 6 )
INTEGER dsssd_e( 6 ), dsssd_f( 6 ), dsssd_g( 6 ), dsssd_h( 6 )
INTEGER*8 dt, tsdata(0:4096), itsdata(0:4096)
INTEGER ediff, fee(32)
INTEGER hec_e_sum_x(2), hec_e_sum_y(2)
INTEGER hec_e_x_max(2), hec_e_y_max(2)
INTEGER hec_e_x_max_ch(2), hec_e_y_max_ch(2)
INTEGER lec_e(2), lec_x(2), lec_y(2)
INTEGER hec_x(2), hec_y(2)
INTEGER m_veto, veto( 1:64 )
INTEGER*8 e_time, e_time_old
INTEGER i
INTEGER*2 i2(2), i2_2(4)
INTEGER i4, i4_2(2)
INTEGER*8 i8, first_ts_value
INTEGER*8 items, items_old, events, events_old
INTEGER idata(0:4095), ierr, information, information_index
INTEGER invalid_id, j, j1, j1_old, j2, j2_old, k
INTEGER j3, j3_old, j4, j4_old
INTEGER l, m, mbs_data(0:3), module, m_disc, ptr, lec(8)
INTEGER m_p_lec(8), m_n_lec(8), m_p_hec(8), m_n_hec(8)
INTEGER sc_channel
INTEGER scaler( 32 ), scaler_old( 32 ), range, time_warp, total
INTEGER*8 timestamp, ts
INTEGER ts28, ts48, ts64
INTEGER*8 old_ts( 8, 2 ), ts_old
INTEGER implant_e( 0:383, 0:127, 1:8 )
INTEGER*8 implant_ts( 0:383, 0:127, 1:8 )
INTEGER*8 decay_ts( 0:383, 0:127, 1:8 )
INTEGER*8 dts( 0:511 ), cts( 0:511 ), min_dts( 0:511 )
INTEGER type, w(0:31), zzz, pulser, pulser_old
INTEGER z_lec, z_lec_veto, z_hec
C
LOGICAL hit(0:4095), ihit(0:4095)
LOGICAL first_event, next_event, first_ts
LOGICAL implant1, implant2, decay1, decay2, other1, other2
LOGICAL dt_tag, warp
LOGICAL first_info_code_4, first_info_code_5
REAL c(0:7), delta_t, gain(0:max-1), offset(0:max-1)
REAL t1, t2, t2_old
DOUBLE PRECISION rates( 32 ), t3, t4, t4_old, t5
C----67---------------------------------------------------------------72------80
C
C Namelists
C
NAMELIST /variables/ gain, offset, c, w,
+ dsssd_a, dsssd_b, dsssd_c, dsssd_d,
+ dsssd_e, dsssd_f, dsssd_g, dsssd_h, veto
C----67---------------------------------------------------------------72------80
C
C Common variables
C
INTEGER*2 gdata(0:15)
C
COMMON /fdata/ gdata
C
C Common variables
C
INTEGER*2 gid(0:15)
C
COMMON /fid/ gid
C
C Common variables
C
INTEGER mult
C
COMMON /fmult/ mult
C
EQUIVALENCE ( i8, i4_2(1) )
EQUIVALENCE ( i4, i2(1) )
C
C----67---------------------------------------------------------------72------80
WRITE ( 6, * ) ' *** FEB25 - implantdecay - March 2025'
WRITE ( 6, * ) ' *** Entry init commences'
C Initialise counters
events = 0
events_old = 0
pulser = 0
pulser_old = 0
time_warp = 0
invalid_id = 0
first_ts = .TRUE.
first_info_code_4 = .false.
first_info_code_5 = .false.
t3 = 0.0D+00
t4 = 0.0D+00
t4_old = 0.0D+00
C Initialise event data
DO i = 0, 4095
data( i ) = 0
tsdata(i) = 0
hit( i ) = .FALSE.
ENDDO
DO i = 1, 32
fee(i) = 0
ENDDO
DO i = 0, 511
min_dts( i ) = 999999999
cts( i ) = 0
ENDDO
first_event = .true.
next_event = .false.
ts_old = 0
e_time_old = 0
count = 0
total = 0
DO i = 0, 4095
gain( i ) = 1.0
offset( i ) = 0.0
ENDDO
DO i = 0, 31
w( i ) = 0
ENDDO
DO i = 0, 7
c( i ) = 1.0
ENDDO
DO i = 1, 64
veto(i) = -1
ENDDO
DO i = 1, 32
scaler( i ) = 0
scaler_old( i ) = 0
rates( i ) = 0.0
ENDDO
DO i = 1, 6
DO j = 1, 2
old_ts( i, j ) = 0
ENDDO
ENDDO
DO i = 0, 383
DO j = 0, 127
DO k = 1, 8
implant_ts( i, j, k ) = 0
decay_ts( i, j , k ) = 0
implant_e( i, j, k ) = 0
ENDDO
ENDDO
ENDDO
C Initialise time
t1 = SECNDS(0.0)
C Read program variables via NAMELIST I/O
OPEN( 1, FILE = '/home/td/FEB25/variables.dat',
+ IOSTAT = ierr )
IF ( ierr.NE.0 ) THEN
WRITE( 6, * ) ' *** OPEN I/O error:', ierr
ENDIF
READ( 1, NML = variables, IOSTAT = ierr )
IF ( ierr.NE.0 ) THEN
WRITE( 6, * ) ' *** READ I/O error:', ierr
ENDIF
CLOSE( 1, IOSTAT = ierr )
IF ( ierr.NE.0 ) THEN
WRITE( 6, * ) ' *** CLOSE I/O error:', ierr
ENDIF
OPEN( 20, FILE = '/home/td/FEB25/20',STATUS='NEW',IOSTAT=ierr )
IF ( ierr.NE.0 ) THEN
WRITE( 6, * ) ' *** OPEN I/O error:', ierr
ENDIF
OPEN( 21, FILE = '/home/td/FEB25/21',STATUS='NEW',IOSTAT=ierr )
IF ( ierr.NE.0 ) THEN
WRITE( 6, * ) ' *** OPEN I/O error:', ierr
ENDIF
C Display program variables
DO i = 0, 511, 8
WRITE( 6, 9001 ) i, i+7, ( gain( j ), j = i, i+7 )
ENDDO
DO i = 0, 511, 8
WRITE( 6, 9002 ) i, i+7, ( offset( j ), j = i, i+7 )
ENDDO
DO i = 0, 7, 8
WRITE( 6, 9003 ) i, i+7, ( c(j), j = 1, 8 )
ENDDO
DO i = 0, 31, 8
WRITE( 6, 9004 ) i, i+7, ( w(j), j = i, i+7 )
ENDDO
DO i = 1, 2
WRITE( 6, 9005 ) i, dsssd_a(i), dsssd_b(i),
+ dsssd_c(i), dsssd_d(i),
+ dsssd_e(i), dsssd_f(i),
+ dsssd_g(i), dsssd_h(i)
ENDDO
WRITE ( 6, * ) ' *** Entry init ends'
RETURN
C----67---------------------------------------------------------------72------80
ENTRY sortin
items = items + 1
IF ( MOD( items, 10000000 ).EQ.0 ) THEN
t2 = SECNDS( t1 )
delta_t = t2 - t2_old
... 1425 more lines ...
|