AIDA GELINA BRIKEN nToF CRIB ISOLDE CIRCE nTOFCapture DESPEC DTAS EDI_PSA 179Ta CARME StellarModelling DCF K40
  AIDA, Page 27 of 46  ELOG logo
New entries since:Thu Jan 1 01:00:00 1970
ID Date Author Subjectdown
  758   Fri Nov 2 01:19:49 2018 TDOffline analysis R9 - alpha background
 
Offline analysis of Oct18/R9_0 to 17 (c. 46h data)

alpha background, BNC PB-4 pulser OFF, LEC/MEC slow comparator 0x64 (~1MeV)

*** ENTRY finish

 *** data items:  404616704 (  10823.34 Hz)
 *** ADC events: ********** ( -30856.57 Hz)
 *** time warps:          0 (      0.00 Hz)
 *** DSSSD # 1 count:      1004 old count:      1004 dt: 160265.10 s  LEC rate:      0.01 Hz
 *** DSSSD # 2 count:       423 old count:       423 dt: 160265.10 s  LEC rate:      0.00 Hz
 *** DSSSD # 3 count:       364 old count:       364 dt: 160265.10 s  LEC rate:      0.00 Hz
 *** DSSSD # 4 count:       363 old count:       363 dt: 160265.10 s  LEC rate:      0.00 Hz
 *** DSSSD # 5 count:       455 old count:       455 dt: 160265.10 s  LEC rate:      0.00 Hz
 *** DSSSD # 6 count:       590 old count:       590 dt: 160265.10 s  LEC rate:      0.00 Hz
 *** DSSSD # 1 count:         0 old count:         0 dt: 160265.10 s  HEC rate:      0.00 Hz
 *** DSSSD # 2 count:         0 old count:         0 dt: 160265.10 s  HEC rate:      0.00 Hz
 *** DSSSD # 3 count:         0 old count:         0 dt: 160265.10 s  HEC rate:      0.00 Hz
 *** DSSSD # 4 count:         0 old count:         0 dt: 160265.10 s  HEC rate:      0.00 Hz
 *** DSSSD # 5 count:         0 old count:         0 dt: 160265.10 s  HEC rate:      0.00 Hz
 *** DSSSD # 6 count:         0 old count:         0 dt: 160265.10 s  HEC rate:      0.00 Hz

 *** ENTRY finish ends

attachment 1 - per DSSSD m_p versus m_n
           2 - per DSSSD x versus y
           3 - per DSSSD E_p versus E_n
           4 - DSSSD channel versus E_p

           5 - implantdecay4.f
           6 - titles
           7 - variables 
               N.B. ADC offsets not set, 0 by default
                    | E_p - E_n | < 2000


DSSSD configuration looks OK. No obvious issues with missing ASICs  
Attachment 1: 20.png
20.png
Attachment 2: 21.png
21.png
Attachment 3: 22.png
22.png
Attachment 4: 23.png
23.png
Attachment 5: implantdecay4.f
*trigger
128  
*oned
1..32       s 65536
101..132    s 65536
201..232    s 65536
250..254    s 4096
301..332    s 65536
400         s 64
500..511    s 4096
2130..2135  s 65536
2030..2035  s 65536
2200..2205  s 65536
2210..2215  s 65536
*twod
1000..1001  s 4096 4096
1100..1105  s 512 512
2000..2005  s 128 128
2010..2015  s 128 128
2020..2025  s 512 512
2100..2105  s 128 128
2110..2115  s 128 128
2120..2125  s 512 512
2300..2305  s 1024 1024
2310..2315  s 1024 1024
2400..2405  s 1024 1024
2410..2415  s 1024 1024
2500..2505  s 1024 1024
2510..2515  s 1024 1024
2600..2605  s 256 256
2700..2705  s 1024 1024
2710..2715  s 1024 1024
2800..2805  s 1024 1024
2810..2815  s 1024 1024
2900..2905  s 1024 1024
2910..2915  s 1024 1024
3000..3005  s 1024 1024
3010..3015  s 1024 1024
3100..3105  s 1024 1024
4000        s 64 64
*vars
*sort
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, channel, channel_ident
      INTEGER count, ctr, data(0:4095), disc
      INTEGER dsssd_a( 6 ), dsssd_b( 6 ), dsssd_c( 6 ), dsssd_d( 6 )
      INTEGER*8 dt, dt1, dt2
      INTEGER ediff, events, events_old, fee(32)
      INTEGER hec_e_sum_x, hec_e_sum_y
      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 idata(0:4095), ierr, information, information_index
      INTEGER invalid_id, items, items_old, j, j1, j1_old, j2, j2_old, k
      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 next_ch, next_adc_data, next_module
      INTEGER sc_channel
      INTEGER scaler( 32 ), scaler_old( 32 ), range, time_warp, total
      INTEGER*8 ts
      INTEGER ts28
      INTEGER ts48
      INTEGER*8 old_ts( 8, 2 ), ts_old, next_ts
      INTEGER implant_e( 0:127, 0:127, 1:8 )
      INTEGER*8 implant_ts( 0:127, 0:127, 1:8 )
      INTEGER*8 decay_ts( 0:127, 0:127, 1:8 )
      INTEGER*8 old_decay_ts( 0:127, 0:127, 1:8 )
      INTEGER*8 older_decay_ts( 0:127, 0:127, 1:8 )
      INTEGER type, w(0:31), zzz
      INTEGER z_lec, z_lec_veto, z_hec
      INTEGER*8 hit_ts(0:4095), ihit_ts(0:4095), tdiff1, tdiff2
       
C
      LOGICAL first, first_ts, hit(0:4095), ihit(0:4095)
      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

C----67---------------------------------------------------------------72------80
C
C     Namelists
C
      NAMELIST /variables/ gain, offset, c, w,
     +                     dsssd_a, dsssd_b, dsssd_c, dsssd_d

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, * ) ' *** AIDA @ RIKEN - implant-decay - June 2017'
      WRITE ( 6, * ) ' *** Entry init commences'

C     Initialise counters

      events = 0
      events_old = 0

      time_warp = 0
      invalid_id = 0

      first_ts = .TRUE.
      t3 = 0.0D+00
      t4 = 0.0D+00
      t4_old = 0.0D+00

C     Initialise event data

      DO i = 0, 4095
       data( i ) = 0
       hit( i ) = .FALSE.
       hit_ts( i ) = 0
       ihit_ts( i ) = 0
      ENDDO
      DO i = 1, 32
       fee(i) = 0
      ENDDO

      first = .true.
      ts_old = 0
      e_time_old = 0
      count = 0
      total = 0
      ctr = 0

      DO i = 0, 127
       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, 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

C     Initialise time

      t1 = SECNDS(0.0)

C     Read program variables via NAMELIST I/O

      OPEN( 1, FILE = '/homes/npg/td/GREAT/variables.dat',
     +         IOSTAT = ierr )
      IF ( ierr.NE.0 ) THEN
       WRITE( 6, * ) ' *** OPEN I/O error:', ierr
       RETURN
      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
       RETURN
      ENDIF  

C     Display program variables

      DO i = 0, 4095, 8
       WRITE( 6, 9001 ) i, i+7, ( gain( j ), j = i, i+7 )
      ENDDO

      DO i = 0, 4095, 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, 6
       WRITE( 6, 9005 ) i, dsssd_a(i), dsssd_b(i),
     +                     dsssd_c(i), dsssd_d(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
       WRITE( 6, 9010 ) items, FLOAT( items - items_old ) / delta_t,
     +                  events, FLOAT( events - events_old ) / delta_t
       items_old = items
       events_old = events
       t2_old = t2

      ENDIF

      IF ( MOD( items, 50000 ).EQ.0 ) THEN

       t4 = DFLOAT( ts ) * 10.0D-9 - t3
       delta_t = t4 - t4_old

       IF ( delta_t.GE.1.0D+00 ) THEN

        DO i = 1, 32
         rates( i ) = DFLOAT( scaler( i ) - scaler_old( i ) ) / delta_t
        ENDDO
C        WRITE( 6, * ) ' '
        DO i = 1, 6
C         WRITE( 6, 9020 ) i, scaler( i ), scaler_old( i ),
C     +                    delta_t, rates( i )
         sc_channel = INT ( MOD( t4, 65536.0 ) + 0.5 )
         CALL set1d( i+300, sc_channel, INT( rates( i ) + 0.5 ) )
        ENDDO
        DO i = 7, 12
C         WRITE( 6, 9030 ) i-6, scaler( i ), scaler_old( i ),
C     +                    delta_t, rates( i )
         sc_channel = INT ( MOD( t4, 65536.0 ) + 0.5 )
         CALL set1d( i+300, sc_channel, INT( rates( i ) + 0.5 ) )
        ENDDO
        DO i = 1, 32
         scaler_old( i ) = scaler( i )
        ENDDO
        t4_old = t4
       ENDIF

      ENDIF

C----67---------------------------------------------------------------72------80
C     Event data

C     For Linux/x64 & Windows/x86 least significant 16-bit word i2(1)
C     For Solaris/SPARC least significant 16-bit word i2(2)

C     GREAT 3.2.1 format
... 1616 more lines ...
Attachment 6: titles.dat
s1 nnaida1: ADC data items (20us/ch)
s2 nnaida2: ADC data items (20us/ch)
s3 nnaida3: ADC data items (20us/ch)
s4 nnaida4: ADC data items (20us/ch)
s5 nnaida5: ADC data items (20us/ch)
s6 nnaida6: ADC data items (20us/ch)
s7 nnaida7: ADC data items (20us/ch)
s8 nnaida8: ADC data items (20us/ch)
s9 nnaida9: ADC data items (20us/ch)
s10 nnaida10: ADC data items (20us/ch)
s11 nnaida11: ADC data items (20us/ch)
s12 nnaida12: ADC data items (20us/ch)
s13 nnaida13: ADC data items (20us/ch)
s14 nnaida14: ADC data items (20us/ch)
s15 nnaida15: ADC data items (20us/ch)
s16 nnaida16: ADC data items (20us/ch)
s17 nnaida17: ADC data items (20us/ch)
s18 nnaida18: ADC data items (20us/ch)
s19 nnaida19: ADC data items (20us/ch)
s20 nnaida20: ADC data items (20us/ch)
s21 nnaida21: ADC data items (20us/ch)
s22 nnaida22: ADC data items (20us/ch)
s23 nnaida23: ADC data items (20us/ch)
s24 nnaida24: ADC data items (20us/ch)
s25 nnaida25: ADC data items (20us/ch)
s26 nnaida26: ADC data items (20us/ch)
s27 nnaida27: ADC data items (20us/ch)
s28 nnaida28: ADC data items (20us/ch)
s29 nnaida29: ADC data items (20us/ch)
s30 nnaida30: ADC data items (20us/ch)
s31 nnaida31: ADC data items (20us/ch)
s32 nnaida32: ADC data items (20us/ch)

s101 nnaida1: HEC ADC data items (20us/ch)
s124 nnaida1: HEC ADC data items (20us/ch)

s201 nnaida1: Disc data items (20us/ch)
s202 nnaida2: Disc data items (20us/ch)
s203 nnaida3: Disc data items (20us/ch)
s204 nnaida4: Disc data items (20us/ch)
s205 nnaida5: Disc data items (20us/ch)
s206 nnaida6: Disc data items (20us/ch)
s207 nnaida7: Disc data items (20us/ch)
s208 nnaida8: Disc data items (20us/ch)
s209 nnaida9: Disc data items (20us/ch)
s210 nnaida10: Disc data items (20us/ch)
s211 nnaida11: Disc data items (20us/ch)
s212 nnaida12: Disc data items (20us/ch)
s213 nnaida13: Disc data items (20us/ch)
s214 nnaida14: Disc data items (20us/ch)
s215 nnaida15: Disc data items (20us/ch)
s216 nnaida16: Disc data items (20us/ch)
s217 nnaida17: Disc data items (20us/ch)
s218 nnaida18: Disc data items (20us/ch)
s219 nnaida19: Disc data items (20us/ch)
s220 nnaida20: Disc data items (20us/ch)
s221 nnaida21: Disc data items (20us/ch)
s222 nnaida22: Disc data items (20us/ch)
s223 nnaida23: Disc data items (20us/ch)
s224 nnaida24: Disc data items (20us/ch)
s225 nnaida25: Disc data items (20us/ch)
s226 nnaida26: Disc data items (20us/ch)
s227 nnaida27: Disc data items (20us/ch)
s228 nnaida28: Disc data items (20us/ch)
s229 nnaida29: Disc data items (20us/ch)
s230 nnaida30: Disc data items (20us/ch)
s231 nnaida31: Disc data items (20us/ch)
s232 nnaida32: Disc data items (20us/ch)

s250 LEC & HEC channels (DSSSD order)
s251 LEC & HEC channels (ASIC order)

s301 DSSSD #1 LEC rate (Hz) versus time 
s306 DSSSD #6 LEC rate (Hz) versus time
s307 DSSSD #1 HEC rate (Hz) versus time 
s312 DSSSD #6 HEC rate (Hz) versus time

s1000 4*(module # - 1) + asic # (0-95 ADC; 100-195 Disc) v. timestamp (20us/ch)
s1001 HEC channel (2048-3584) versus energy

s2000 m_p versus m_n DSSSD #1 (LEC)
s2005 m_p versus m_n DSSSD #6 (LEC)
s2010 x versus y DSSSD #1 (LEC)
s2015 x versus y DSSSD #6 (LEC)
s2020 E_p versus E_n DSSSD #1 (LEC)
s2025 E_p versus E_n DSSSD #6 (LEC)
s2030 LEC-LEC time DSSSD #1 (2.56us/ch)
s2035 LEC-LEC time DSSSD #6 (2.56us/ch)

s2100 m_p versus m_n DSSSD #1 (HEC)
s2105 m_p versus m_n DSSSD #6 (HEC)
s2110 x versus y DSSSD #1 (HEC)
s2115 x versus y DSSSD #6 (HEC)
s2120 E_p versus E_n DSSSD #1 (HEC)
s2125 E_p versus E_n DSSSD #6 (HEC)

s2130 HEC-HEC time DSSSD #1 (2.56us/ch)
s2135 HEC-HEC time DSSSD #6 (2.56us/ch)

s2200 per pixel HEC-LEC time DSSSD #1 (2.56us/ch)
s2205 per pixel HEC-LEC time DSSSD #6 (2.56us/ch)

s2210 per pixel HEC-LEC time DSSSD #1 (163.840us/ch)
s2215 per pixel HEC-LEC time DSSSD #6 (163.840us/ch)

s2300 HEC E versus HEC-LEC dt DSSSD#1 (163.840us/ch)
s2305 HEC E versus HEC-LEC dt DSSSD#6 (163.840us/ch)

s2400 LEC E versus HEC-LEC dt DSSSD#1 (163.840us/ch)
s2405 LEC E versus HEC-LEC dt DSSSD#6 (163.840us/ch)

s2500 HEC xy versus HEC-LEC dt DSSSD#1 (163.840us/ch)
s2505 HEC xy versus HEC-LEC dt DSSSD#6 (163.840us/ch)

s2600 LEC x_time - e_time DSSSD#1 (10ns/ch)
s2605 LEC x_time - e_time DSSSD#6 (10ns/ch)

s2700 e_diff versus HEC-LEC dt DSSSD#1 (163.840us/ch)
s2705 e_diff versus HEC-LEC dt DSSSD#6 (163.840us/ch)

s2800 m/z_lec/ctr versus HEC-LEC dt DSSSD#1 (163.840us/ch)
s2805 m/z_lec/ctr versus HEC-LEC dt DSSSD#6 (163.840us/ch)

s2900 LEC-LEC dt (20.48us/ch) versus HEC-LEC dt DSSSD#1 (163.840us/ch)
s2905 LEC-LEC dt (20.48us/ch) versus HEC-LEC dt DSSSD#6 (163.840us/ch)

s3000 FEE64 versus FEE64
Attachment 7: variables.dat
  $variables
  dsssd_a( 1 ) = 17
  dsssd_a( 2 ) = 18
  dsssd_a( 3 ) = 9
  dsssd_a( 4 ) = 10
  dsssd_a( 5 ) = 1
  dsssd_a( 6 ) = 2
  dsssd_b( 1 ) = 20
  dsssd_b( 2 ) = 19
  dsssd_b( 3 ) = 12
  dsssd_b( 4 ) = 11
  dsssd_b( 5 ) = 4
  dsssd_b( 6 ) = 3
  dsssd_c( 1 ) = 22
  dsssd_c( 2 ) = 21
  dsssd_c( 3 ) = 14
  dsssd_c( 4 ) = 13
  dsssd_c( 5 ) = 6
  dsssd_c( 6 ) = 5
  dsssd_d( 1 ) = 23
  dsssd_d( 2 ) = 24
  dsssd_d( 3 ) = 15
  dsssd_d( 4 ) = 16
  dsssd_d( 5 ) = 7
  dsssd_d( 6 ) = 8
  offset = 4096*0.0
  gain = 4096*1.0
  c = 8*1.0
  w = 32*0
  w(1) = 2000
  w(2) = -2000
  w(3) = 210
  w(4) = 210
  $[end]
  764   Sun Nov 4 04:31:49 2018 TDOffline analysis R4
Offline analysis of Oct18/R4_0 to _2

alpha background and pulser

Attachments 1-24 - per FEE64 ADC spectra (32k->4k)

Attachments 25-26 - expanded view pulser peaks FEE64 #2 & #17 & #24
                    anomalously small pulser peak widths observed -> possible LH Kapton PCB issue

Attachements 27-28 - sort program (pulser.f) and variables (variables.dat)
Attachment 1: 101.png
101.png
Attachment 2: 102.png
102.png
Attachment 3: 103.png
103.png
Attachment 4: 104.png
104.png
Attachment 5: 105.png
105.png
Attachment 6: 106.png
106.png
Attachment 7: 107.png
107.png
Attachment 8: 108.png
108.png
Attachment 9: 109.png
109.png
Attachment 10: 110.png
110.png
Attachment 11: 111.png
111.png
Attachment 12: 112.png
112.png
Attachment 13: 113.png
113.png
Attachment 14: 114.png
114.png
Attachment 15: 115.png
115.png
Attachment 16: 116.png
116.png
Attachment 17: 117.png
117.png
Attachment 18: 118.png
118.png
Attachment 19: 119.png
119.png
Attachment 20: 120.png
120.png
Attachment 21: 121.png
121.png
Attachment 22: 122.png
122.png
Attachment 23: 123.png
123.png
Attachment 24: 124.png
124.png
Attachment 25: 130.png
130.png
Attachment 26: 131.png
131.png
Attachment 27: pulser.f
*trigger
128  
*oned
0..1535     s 4096
*twod
*vars
*sort
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, channel, channel_ident
      INTEGER count, data(0:4095), disc
      INTEGER dsssd_a( 6 ), dsssd_b( 6 ), dsssd_c( 6 ), dsssd_d( 6 )
      INTEGER*8 dt
      INTEGER events, events_old
      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 idata(0:4095), ierr, information, information_index
      INTEGER invalid_id, items, items_old, j, j1, j1_old, j2, j2_old, k
      INTEGER l, mbs_data(0:3), module, m_disc, m_p, m_n, ptr
      INTEGER sc_channel
      INTEGER scaler( 32 ), scaler_old( 32 ), range, time_warp, total
      INTEGER*8 ts
      INTEGER ts28
      INTEGER ts48
      INTEGER*8 ts_old
      INTEGER type, w(0:31), zzz
C
      LOGICAL first, first_ts, hit(0:4095), ihit(0:4095)
      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

C----67---------------------------------------------------------------72------80
C
C     Namelists
C
      NAMELIST /variables/ gain, offset, c, w,
     +                     dsssd_a, dsssd_b, dsssd_c, dsssd_d

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, * ) ' *** AIDA @ RIKEN - pulser - May 2017'
      WRITE ( 6, * ) ' *** Entry init commences'

C     Initialise counters

      events = 0
      events_old = 0

      time_warp = 0
      invalid_id = 0

      first_ts = .TRUE.
      t3 = 0.0D+00
      t4 = 0.0D+00
      t4_old = 0.0D+00

C     Initialise event data

      DO i = 0, 4095
       data( i ) = 0
       hit( i ) = .FALSE.
      ENDDO

      first = .true.
      ts_old = 0
      e_time_old = 0
      count = 0
      total = 0

      DO i = 0, 127
       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, 32
       scaler( i ) = 0
       scaler_old( i ) = 0
       rates( i ) = 0.0
      ENDDO

C     Initialise time

      t1 = SECNDS(0.0)

C     Read program variables via NAMELIST I/O

      OPEN( 1, FILE = '/homes/npg/td/GREAT/variables.dat',
     +         IOSTAT = ierr )
      IF ( ierr.NE.0 ) THEN
       WRITE( 6, * ) ' *** OPEN I/O error:', ierr
       RETURN
      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
       RETURN
      ENDIF  

C     Display program variables

      DO i = 0, 4095, 8
       WRITE( 6, 9001 ) i, i+7, ( gain( j ), j = i, i+7 )
      ENDDO

      DO i = 0, 4095, 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, 6
       WRITE( 6, 9005 ) i, dsssd_a(i), dsssd_b(i),
     +                     dsssd_c(i), dsssd_d(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
       WRITE( 6, 9010 ) items, FLOAT( items - items_old ) / delta_t,
     +                  events, FLOAT( events - events_old ) / delta_t
       items_old = items
       events_old = events
       t2_old = t2

      ENDIF

      IF ( MOD( items, 50000 ).EQ.0 ) THEN

       t4 = DFLOAT( ts ) * 10.0D-9 - t3
       delta_t = t4 - t4_old

       IF ( delta_t.GE.1.0D+00 ) THEN

        DO i = 1, 32
         rates( i ) = DFLOAT( scaler( i ) - scaler_old( i ) ) / delta_t
        ENDDO
C        WRITE( 6, * ) ' '
        DO i = 1, 6
C         WRITE( 6, 9020 ) i, scaler( i ), scaler_old( i ),
C     +                    delta_t, rates( i )
         sc_channel = INT ( MOD( t4, 65536.0 ) + 0.5 )
C         CALL set1d( i+300, sc_channel, INT( rates( i ) + 0.5 ) )
        ENDDO
        DO i = 1, 32
         scaler_old( i ) = scaler( i )
        ENDDO
        t4_old = t4
       ENDIF

      ENDIF

C----67---------------------------------------------------------------72------80
C     Event data

C     For Linux/x64 & Windows/x86 least significant 16-bit word i2(1)
C     For Solaris/SPARC least significant 16-bit word i2(2)

C     GREAT 3.2.1 format
C
C     Other information
C      gid[0]=-1
C     gdata[0]=module number
C     gdata[1]=information type (2=Pause, 3=Resume, 4=SYNC100, 6=FEE64 disc etc)
C     gdata[2]=information field (bits 16-19)
C     gdata[3]=information field (bits 0-15)
C     gdata[4]=time stamp (bits 16-27)
C     gdata[5]=time stamp (bits 0-15) 
C
C     ADC data format
C
C       gid[0]=channel ident (FEE64 module bits 6-11, channel bits 0-5)
C     gdata[0]=fail/veto bits (AIDA: fail=0, veto=0 LEC/MEC, veto=1 HEC)
C     gdata[1]=?
C     gdata[2]=ADC data
C     gdata[3]=time stamp (bits 16-27)
C     gdata[4]=time stamp (bits 0-15) 

C----67---------------------------------------------------------------72------80
C     gid = - 1   Other information

      IF ( gid(0).EQ.-1 ) THEN

       i2(2) = 0
       i2(1) = gdata(0)
       module = AND( i4, Z'0000003f' )
       i2(1) = gdata(1)
       type = AND( i4, Z'0000000f' )
       i2(2) = gdata(2)
       i2(1) = gdata(3)
       information = AND( i4, Z'000fffff' )
       i2(2) = gdata(4)
       i2(1) = gdata(5)
       ts28 = AND( i4, Z'0fffffff' )

       IF ( type.EQ.2 ) THEN

C----67---------------------------------------------------------------72------80
C     PAUSE timestamp

        ts48 = information

        i8 = 0
        i4_2( 1 ) = ts48
        ts = LSHIFT( i8, 28 )
        i8 = 0
        i4_2(1) = ts28
        ts = ts + i8

       ELSEIF( type.EQ.3 ) THEN

C----67---------------------------------------------------------------72------80
C     RESUME timestamp

        ts48 = information

        i8 = 0
        i4_2( 1 ) = ts48
        ts = LSHIFT( i8, 28 )
        i8 = 0
        i4_2(1) = ts28
        ts = ts + i8

... 322 more lines ...
Attachment 28: variables.dat
  $variables
  dsssd_a( 1 ) = 17
  dsssd_a( 2 ) = 18
  dsssd_a( 3 ) = 9
  dsssd_a( 4 ) = 10
  dsssd_a( 5 ) = 1
  dsssd_a( 6 ) = 2
  dsssd_b( 1 ) = 20
  dsssd_b( 2 ) = 19
  dsssd_b( 3 ) = 12
  dsssd_b( 4 ) = 11
  dsssd_b( 5 ) = 4
  dsssd_b( 6 ) = 3
  dsssd_c( 1 ) = 22
  dsssd_c( 2 ) = 21
  dsssd_c( 3 ) = 14
  dsssd_c( 4 ) = 13
  dsssd_c( 5 ) = 6
  dsssd_c( 6 ) = 5
  dsssd_d( 1 ) = 23
  dsssd_d( 2 ) = 24
  dsssd_d( 3 ) = 15
  dsssd_d( 4 ) = 16
  dsssd_d( 5 ) = 7
  dsssd_d( 6 ) = 8
  offset = 4096*0.0
  gain = 4096*1.0
  c = 8*1.0
  w = 32*0
  w(1) = 2000
  w(2) = -2000
  w(3) = 210
  w(4) = 210
  $[end]
  501   Wed Nov 30 05:30:20 2016 TDOffline analysis R39_6
Degrader 8mm Al fixed + ( 0.5mm W/1mm Al + 1mm W/1mm Al) variable

S O R T    C O M P L E T E ..... Wed Nov 30 14:41:39 2016


 *** ENTRY finish

 *** data items:  261088000 ( 107267.57 Hz)
 *** ADC events:   36298818 (  14913.31 Hz)
 *** time warps:          0 (      0.00 Hz)
 *** DSSSD # 1 count:   5026250 old count:   5011592 dt:    301.12 s  LEC rate:  16691.77 Hz
 *** DSSSD # 2 count:   4868135 old count:   4854528 dt:    301.12 s  LEC rate:  16166.68 Hz
 *** DSSSD # 3 count:   3233592 old count:   3224663 dt:    301.12 s  LEC rate:  10738.50 Hz
 *** DSSSD # 4 count:   2914929 old count:   2906700 dt:    301.12 s  LEC rate:   9680.24 Hz
 *** DSSSD # 5 count:   2009399 old count:   2003732 dt:    301.12 s  LEC rate:   6673.05 Hz
 *** DSSSD # 6 count:   3695168 old count:   3684926 dt:    301.12 s  LEC rate:  12271.36 Hz
 *** DSSSD # 1 count:        24 old count:        24 dt:    301.12 s  HEC rate:      0.08 Hz
 *** DSSSD # 2 count:        12 old count:        12 dt:    301.12 s  HEC rate:      0.04 Hz
 *** DSSSD # 3 count:        39 old count:        39 dt:    301.12 s  HEC rate:      0.13 Hz
 *** DSSSD # 4 count:        94 old count:        94 dt:    301.12 s  HEC rate:      0.31 Hz
 *** DSSSD # 5 count:       152 old count:       152 dt:    301.12 s  HEC rate:      0.50 Hz
 *** DSSSD # 6 count:       131 old count:       131 dt:    301.12 s  HEC rate:      0.44 Hz

 *** ENTRY finish ends

S O R T    S T O P P E D ..... Wed Nov 30 14:42:53 2016
Attachment 1: 30.png
30.png
Attachment 2: 31.png
31.png
Attachment 3: 32.png
32.png
Attachment 4: 33.png
33.png
Attachment 5: 34.png
34.png
Attachment 6: 35.png
35.png
Attachment 7: 36.png
36.png
Attachment 8: 37.png
37.png
  500   Wed Nov 30 04:56:40 2016 TDOffline analysis R38_33
Degrader 8mm Al fixed + ( 0.3mm W/1mm Al + 1mm W/1mm Al) variable

ADC & disc data synchronised

S O R T    C O M P L E T E ..... Wed Nov 30 13:53:38 2016


 *** ENTRY finish

 *** data items:  261088000 (  95225.14 Hz)
 *** ADC events:   44184529 (  16115.17 Hz)
 *** time warps:          0 (      0.00 Hz)
 *** DSSSD # 1 count:   5050939 old count:   5048005 dt:    362.93 s  LEC rate:  13917.19 Hz
 *** DSSSD # 2 count:   5030612 old count:   5028047 dt:    362.93 s  LEC rate:  13861.19 Hz
 *** DSSSD # 3 count:   2817584 old count:   2815988 dt:    362.93 s  LEC rate:   7763.48 Hz
 *** DSSSD # 4 count:   2482552 old count:   2481283 dt:    362.93 s  LEC rate:   6840.34 Hz
 *** DSSSD # 5 count:   1549349 old count:   1548540 dt:    362.93 s  LEC rate:   4269.03 Hz
 *** DSSSD # 6 count:   3960450 old count:   3958332 dt:    362.93 s  LEC rate:  10912.50 Hz
 *** DSSSD # 1 count:        18 old count:        18 dt:    362.93 s  HEC rate:      0.05 Hz
 *** DSSSD # 2 count:         6 old count:         6 dt:    362.93 s  HEC rate:      0.02 Hz
 *** DSSSD # 3 count:        21 old count:        20 dt:    362.93 s  HEC rate:      0.06 Hz
 *** DSSSD # 4 count:        64 old count:        63 dt:    362.93 s  HEC rate:      0.18 Hz
 *** DSSSD # 5 count:       157 old count:       157 dt:    362.93 s  HEC rate:      0.43 Hz
 *** DSSSD # 6 count:       200 old count:       200 dt:    362.93 s  HEC rate:      0.55 Hz

 *** ENTRY finish ends

S O R T    S T O P P E D ..... Wed Nov 30 13:55:06 2016
Attachment 1: 20.png
20.png
Attachment 2: 21.png
21.png
Attachment 3: 23.png
23.png
Attachment 4: 24.png
24.png
Attachment 5: 25.png
25.png
Attachment 6: 26.png
26.png
Attachment 7: 27.png
27.png
Attachment 8: 22.png
22.png
Attachment 9: 28.png
28.png
  503   Wed Nov 30 06:33:52 2016 TDOffline analysis R36_30
Following target change from (CH2)n to C

Degrader (8mm Al) fixed + ( 0.3mm W/1mm Al + 1mm W/1m Al) variable

See https://elog.ph.ed.ac.uk/AIDA/499

S O R T    C O M P L E T E ..... Wed Nov 30 15:27:25 2016


 *** ENTRY finish

 *** data items:  261088000 ( 116005.27 Hz)
 *** ADC events:   36396496 (  16171.50 Hz)
 *** time warps:          0 (      0.00 Hz)
 *** DSSSD # 1 count:   5261837 old count:   5244729 dt:    302.14 s  LEC rate:  17415.17 Hz
 *** DSSSD # 2 count:   4766457 old count:   4750965 dt:    302.14 s  LEC rate:  15775.61 Hz
 *** DSSSD # 3 count:   2750220 old count:   2741230 dt:    302.14 s  LEC rate:   9102.44 Hz
 *** DSSSD # 4 count:   2371078 old count:   2363331 dt:    302.14 s  LEC rate:   7847.59 Hz
 *** DSSSD # 5 count:   1938346 old count:   1932204 dt:    302.14 s  LEC rate:   6415.37 Hz
 *** DSSSD # 6 count:   3579949 old count:   3568500 dt:    302.14 s  LEC rate:  11848.61 Hz
 *** DSSSD # 1 count:        18 old count:        18 dt:    302.14 s  HEC rate:      0.06 Hz
 *** DSSSD # 2 count:         1 old count:         1 dt:    302.14 s  HEC rate:      0.00 Hz
 *** DSSSD # 3 count:        13 old count:        13 dt:    302.14 s  HEC rate:      0.04 Hz
 *** DSSSD # 4 count:        37 old count:        37 dt:    302.14 s  HEC rate:      0.12 Hz
 *** DSSSD # 5 count:       110 old count:       110 dt:    302.14 s  HEC rate:      0.36 Hz
 *** DSSSD # 6 count:       144 old count:       144 dt:    302.14 s  HEC rate:      0.48 Hz

 *** ENTRY finish ends

S O R T    S T O P P E D ..... Wed Nov 30 15:28:35 2016
Attachment 1: 40.png
40.png
Attachment 2: 41.png
41.png
Attachment 3: 42.png
42.png
Attachment 4: 43.png
43.png
Attachment 5: 44.png
44.png
Attachment 6: 45.png
45.png
Attachment 7: 46.png
46.png
Attachment 8: 48.png
48.png
  497   Tue Nov 29 12:21:35 2016 TDOffline analysis R35_1

S O R T    C O M P L E T E ..... Tue Nov 29 20:45:36 2016


 *** ENTRY finish

 *** data items:  261088000 (  86387.95 Hz)
 *** ADC events:   35941414 (  11892.18 Hz)
 *** time warps:          0 (      0.00 Hz)
 *** DSSSD # 1 count:   5276204 old count:   5262371 dt:    300.27 s  LEC rate:  17571.70 Hz
 *** DSSSD # 2 count:   4680694 old count:   4668369 dt:    300.27 s  LEC rate:  15588.43 Hz
 *** DSSSD # 3 count:   2774752 old count:   2767324 dt:    300.27 s  LEC rate:   9240.95 Hz
 *** DSSSD # 4 count:   2405451 old count:   2399010 dt:    300.27 s  LEC rate:   8011.04 Hz
 *** DSSSD # 5 count:   1940411 old count:   1935190 dt:    300.27 s  LEC rate:   6462.28 Hz
 *** DSSSD # 6 count:   3578504 old count:   3569034 dt:    300.27 s  LEC rate:  11917.74 Hz
 *** DSSSD # 1 count:        54 old count:        54 dt:    300.27 s  HEC rate:      0.18 Hz
 *** DSSSD # 2 count:        76 old count:        76 dt:    300.27 s  HEC rate:      0.25 Hz
 *** DSSSD # 3 count:       124 old count:       124 dt:    300.27 s  HEC rate:      0.41 Hz
 *** DSSSD # 4 count:       129 old count:       129 dt:    300.27 s  HEC rate:      0.43 Hz
 *** DSSSD # 5 count:        84 old count:        84 dt:    300.27 s  HEC rate:      0.28 Hz
 *** DSSSD # 6 count:        33 old count:        33 dt:    300.27 s  HEC rate:      0.11 Hz

 *** ENTRY finish ends

S O R T    S T O P P E D ..... Tue Nov 29 20:47:11 2016
Attachment 1: 80.png
80.png
Attachment 2: 81.png
81.png
Attachment 3: 82.png
82.png
Attachment 4: 83.png
83.png
Attachment 5: 84.png
84.png
Attachment 6: 85.png
85.png
Attachment 7: 86.png
86.png
Attachment 8: 87.png
87.png
  491   Mon Nov 28 11:55:20 2016 TDOffline analysis R30_2
Degrader (8mm Al) fixed + ( 0.3mm W/1mm Al + 1mm W/1mm Al)

S O R T    C O M P L E T E ..... Mon Nov 28 20:53:02 2016


 *** ENTRY finish

 *** data items:  261088000 ( 120002.96 Hz)
 *** ADC events:   34802098 (  15995.96 Hz)
 *** time warps:          0 (      0.00 Hz)
 *** DSSSD # 1 count:   5905997 old count:   5887846 dt:    290.51 s  LEC rate:  20329.92 Hz
 *** DSSSD # 2 count:   5530448 old count:   5513790 dt:    290.51 s  LEC rate:  19037.19 Hz
 *** DSSSD # 3 count:   2899015 old count:   2890171 dt:    290.51 s  LEC rate:   9979.14 Hz
 *** DSSSD # 4 count:   2612137 old count:   2604306 dt:    290.51 s  LEC rate:   8991.63 Hz
 *** DSSSD # 5 count:   1944405 old count:   1938459 dt:    290.51 s  LEC rate:   6693.13 Hz
 *** DSSSD # 6 count:   3424115 old count:   3413635 dt:    290.51 s  LEC rate:  11786.66 Hz
 *** DSSSD # 1 count:        52 old count:        52 dt:    290.51 s  HEC rate:      0.18 Hz
 *** DSSSD # 2 count:        78 old count:        78 dt:    290.51 s  HEC rate:      0.27 Hz
 *** DSSSD # 3 count:       147 old count:       147 dt:    290.51 s  HEC rate:      0.51 Hz
 *** DSSSD # 4 count:       130 old count:       130 dt:    290.51 s  HEC rate:      0.45 Hz
 *** DSSSD # 5 count:        84 old count:        84 dt:    290.51 s  HEC rate:      0.29 Hz
 *** DSSSD # 6 count:        40 old count:        40 dt:    290.51 s  HEC rate:      0.14 Hz

 *** ENTRY finish ends

S O R T    S T O P P E D ..... Mon Nov 28 20:54:07 2016
Attachment 1: 20.png
20.png
Attachment 2: 21.png
21.png
Attachment 3: 22.png
22.png
  489   Mon Nov 28 05:38:39 2016 DKOffline analysis R28_14
Conditions: Two Pb walls with empty windows. 3mm Pb degrader between MUSIC 1 and MUSIC 2.

Variable degraders: 1.0 mm W (+1 mm Al) and 0.3 mm W (+1 mm Al)

 *** DSSSD # 1 count:   5901917 old count:   5892055 dt:    286.62 s  LEC rate:  20591.57 Hz
 *** DSSSD # 2 count:   5564387 old count:   5555277 dt:    286.62 s  LEC rate:  19413.94 Hz
 *** DSSSD # 3 count:   2995556 old count:   2990474 dt:    286.62 s  LEC rate:  10451.38 Hz
 *** DSSSD # 4 count:   2453137 old count:   2448680 dt:    286.62 s  LEC rate:   8558.90 Hz
 *** DSSSD # 5 count:   1926370 old count:   1923142 dt:    286.62 s  LEC rate:   6721.03 Hz
 *** DSSSD # 6 count:   3550119 old count:   3544196 dt:    286.62 s  LEC rate:  12386.23 Hz
 *** DSSSD # 1 count:        76 old count:        76 dt:    286.62 s  HEC rate:      0.27 Hz
 *** DSSSD # 2 count:        94 old count:        94 dt:    286.62 s  HEC rate:      0.33 Hz
 *** DSSSD # 3 count:       118 old count:       118 dt:    286.62 s  HEC rate:      0.41 Hz
 *** DSSSD # 4 count:        74 old count:        74 dt:    286.62 s  HEC rate:      0.26 Hz
 *** DSSSD # 5 count:        48 old count:        48 dt:    286.62 s  HEC rate:      0.17 Hz
 *** DSSSD # 6 count:        17 old count:        17 dt:    286.62 s  HEC rate:      0.06 Hz

 *** ENTRY finish ends

S O R T    S T O P P E D ..... Mon Nov 28 14:20:29 2016

See attached histograms
Attachment 1: r28_14-2120.png
r28_14-2120.png
Attachment 2: r28_14-2110.png
r28_14-2110.png
Attachment 3: r28_14-307.png
r28_14-307.png
  490   Mon Nov 28 06:49:04 2016 DKOffline analysis R27_1
Conditions: One Pb wall with empty windows. 3mm Pb degrader closer to AIDA, on the front of PE wall

Variable degraders: 1.0 mm W (+1 mm Al) and 0.3 mm W (+1 mm Al)

 *** data items:  261088000 ( 110875.14 Hz)
 *** ADC events:   39427988 (  16743.72 Hz)
 *** time warps:          0 (      0.00 Hz)
 *** DSSSD # 1 count:   5617273 old count:   5603383 dt:    327.21 s  LEC rate:  17167.27 Hz
 *** DSSSD # 2 count:   6007495 old count:   5992656 dt:    327.21 s  LEC rate:  18359.85 Hz
 *** DSSSD # 3 count:   3154280 old count:   3146595 dt:    327.21 s  LEC rate:   9639.98 Hz
 *** DSSSD # 4 count:   2713744 old count:   2706934 dt:    327.21 s  LEC rate:   8293.63 Hz
 *** DSSSD # 5 count:   1977642 old count:   1972985 dt:    327.21 s  LEC rate:   6043.99 Hz
 *** DSSSD # 6 count:   2943065 old count:   2935788 dt:    327.21 s  LEC rate:   8994.47 Hz
 *** DSSSD # 1 count:       128 old count:       128 dt:    327.21 s  HEC rate:      0.39 Hz
 *** DSSSD # 2 count:       174 old count:       174 dt:    327.21 s  HEC rate:      0.53 Hz
 *** DSSSD # 3 count:       236 old count:       236 dt:    327.21 s  HEC rate:      0.72 Hz
 *** DSSSD # 4 count:       158 old count:       158 dt:    327.21 s  HEC rate:      0.48 Hz
 *** DSSSD # 5 count:        80 old count:        80 dt:    327.21 s  HEC rate:      0.24 Hz
 *** DSSSD # 6 count:        29 old count:        29 dt:    327.21 s  HEC rate:      0.09 Hz

 *** ENTRY finish ends

S O R T    S T O P P E D ..... Mon Nov 28 15:32:15 2016

See attached histograms
Attachment 1: r27_1-307.png
r27_1-307.png
Attachment 2: r27_1-2110.png
r27_1-2110.png
Attachment 3: r27_1-2120.png
r27_1-2120.png
  583   Mon May 22 08:29:26 2017 TDOffline analysis R25_12 to R25_22
Offline analysis of AIDA data files R25_12 to R25_22 using MIDASsort

Note

- no calibration (raw ADC data only)
- no thresholds 
- no clustering
- events delimited by > 2us time difference between successive ADC data words 
- HEC = high energy channel = 20GeV FSR range for high energy implants
- LEC = low energy channel = 20MeV FSR range for decays (and pulser)

Attachment 1-4
--------------
s2000-s2005 per DSSSD LEC m_lec_p versus m_lec_n

            following spectra require 0 < m_lec_p < 8 and 0 < m_lec_n < 8
            all combinations of x & y plotted

s2010-s2015 per DSSSD LEC x versus y
s2020-s2025 per DSSSD LEC E_p versus E_n        (20keV/ch nominal)
s2030-s2035 per DSSSD LEC-LEC time              (2.56us/channel)

Attachment 5-8
---------------
s2000-s2005 per DSSSD LEC m_hec_p versus m_hec_n

            following spectra require 0 < m_hec_p < 8 and 0 < m_hec_n < 8
            all combinations of x & y plotted

s2010-s2015 per DSSSD HEC x versus y           
s2020-s2025 per DSSSD HEC E_p versus E_n       (20MeV/channel nominal)
s2030-s2035 per DSSSD HEC-HEC time             (2.56us/channel)

Attachment 9-11
---------------
s2200-s2205 per pixel, per DSSSD HEC-LEC time  (2.56us/channel)

Attachment 12-14
----------------
s2210-s2215 per pixel, per DSSSD HEC-LEC time  (163.840us/channel)

Attachments 15-17
-----------------
s2300-s2305 per DSSSD HEC Energy versus HEC-LEC time (20MeV/channel, 163.840us/channel)
s2400-s2405 per DSSSD LEC Energy versus HEC-LEC time (20keV/channel, 163.840us/channel)
s2500-s2505 per DSSSD x & y strips versus HEC-LEC time (163.840us/channel)

Attachment 18
-------------
Projection s2503 y-strips 0-60 -> HEC-LEC time 163.804us/channel
Attachment 1: 100.png
100.png
Attachment 2: 101.png
101.png
Attachment 3: 102.png
102.png
Attachment 4: 103.png
103.png
Attachment 5: 104.png
104.png
Attachment 6: 105.png
105.png
Attachment 7: 106.png
106.png
Attachment 8: 107.png
107.png
Attachment 9: 108.png
108.png
Attachment 10: 109.png
109.png
Attachment 11: 110.png
110.png
Attachment 12: 111.png
111.png
Attachment 13: 112.png
112.png
Attachment 14: 113.png
113.png
Attachment 15: 115.png
115.png
Attachment 16: 116.png
116.png
Attachment 17: 118.png
118.png
Attachment 18: 119.png
119.png
Attachment 19: implantdecay.f
*trigger
128  
*oned
1..32       s 65536
101..132    s 65536
200..202    s 4096
301..332    s 65536
400         s 64
2130..2135  s 65536
2030..2035  s 65536
2200..2205  s 65536
2210..2215  s 65536
*twod
1000..1001  s 4096 4096
1100        s 64 64
2000..2005  s 128 128
2010..2015  s 128 128
2020..2025  s 512 512
2100..2105  s 128 128
2110..2115  s 128 128
2120..2125  s 512 512
2300..2305  s 1024 1024
2400..2405  s 1024 1024
2500..2505  s 1024 1024
*vars
*sort
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, channel, channel_ident
      INTEGER count, data(0:4095), disc
      INTEGER dsssd_a( 6 ), dsssd_b( 6 ), dsssd_c( 6 ), dsssd_d( 6 )
      INTEGER*8 dt
      INTEGER events, events_old
      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 idata(0:4095), ierr, information, information_index
      INTEGER invalid_id, items, items_old, j, j1, j1_old, j2, j2_old, k
      INTEGER l, m, mbs_data(0:3), module, m_disc, ptr
      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 ts
      INTEGER ts28
      INTEGER ts48
      INTEGER*8 old_ts( 8, 2 ), ts_old 
      INTEGER implant_e( 0:127, 0:127, 1:8 )
      INTEGER*8 implant_ts( 0:127, 0:127, 1:8 )
      INTEGER*8 decay_ts( 0:127, 0:127, 1:8 )
      INTEGER type, w(0:31), zzz
      INTEGER z_lec, z_hec
C
      LOGICAL first, first_ts, hit(0:4095), ihit(0:4095)
      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

C----67---------------------------------------------------------------72------80
C
C     Namelists
C
      NAMELIST /variables/ gain, offset, c, w,
     +                     dsssd_a, dsssd_b, dsssd_c, dsssd_d

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, * ) ' *** AIDA @ RIKEN - implants - November 2016'
      WRITE ( 6, * ) ' *** Entry init commences'

C     Initialise counters

      events = 0
      events_old = 0

      time_warp = 0
      invalid_id = 0

      first_ts = .TRUE.
      t3 = 0.0D+00
      t4 = 0.0D+00
      t4_old = 0.0D+00

C     Initialise event data

      DO i = 0, 4095
       data( i ) = 0
       hit( i ) = .FALSE.
      ENDDO

      first = .true.
      ts_old = 0
      e_time_old = 0
      count = 0
      total = 0

      DO i = 0, 127
       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, 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

C     Initialise time

      t1 = SECNDS(0.0)

C     Read program variables via NAMELIST I/O

      OPEN( 1, FILE = '/homes/npg/td/GREAT/variables.dat',
     +         IOSTAT = ierr )
      IF ( ierr.NE.0 ) THEN
       WRITE( 6, * ) ' *** OPEN I/O error:', ierr
       RETURN
      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
       RETURN
      ENDIF  

C     Display program variables

      DO i = 0, 4095, 8
       WRITE( 6, 9001 ) i, i+7, ( gain( j ), j = i, i+7 )
      ENDDO

      DO i = 0, 4095, 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, 6
       WRITE( 6, 9005 ) i, dsssd_a(i), dsssd_b(i),
     +                     dsssd_c(i), dsssd_d(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
       WRITE( 6, 9010 ) items, FLOAT( items - items_old ) / delta_t,
     +                  events, FLOAT( events - events_old ) / delta_t
       items_old = items
       events_old = events
       t2_old = t2

      ENDIF

      IF ( MOD( items, 50000 ).EQ.0 ) THEN

       t4 = DFLOAT( ts ) * 10.0D-9 - t3
       delta_t = t4 - t4_old

       IF ( delta_t.GE.1.0D+00 ) THEN

        DO i = 1, 32
         rates( i ) = DFLOAT( scaler( i ) - scaler_old( i ) ) / delta_t
        ENDDO
C        WRITE( 6, * ) ' '
        DO i = 1, 6
C         WRITE( 6, 9020 ) i, scaler( i ), scaler_old( i ),
C     +                    delta_t, rates( i )
         sc_channel = INT ( MOD( t4, 65536.0 ) + 0.5 )
         CALL set1d( i+300, sc_channel, INT( rates( i ) + 0.5 ) )
        ENDDO
        DO i = 7, 12
C         WRITE( 6, 9030 ) i-6, scaler( i ), scaler_old( i ),
C     +                    delta_t, rates( i )
         sc_channel = INT ( MOD( t4, 65536.0 ) + 0.5 )
         CALL set1d( i+300, sc_channel, INT( rates( i ) + 0.5 ) )
        ENDDO
        DO i = 1, 32
         scaler_old( i ) = scaler( i )
        ENDDO
        t4_old = t4
       ENDIF

      ENDIF

C----67---------------------------------------------------------------72------80
C     Event data

C     For Linux/x64 & Windows/x86 least significant 16-bit word i2(1)
C     For Solaris/SPARC least significant 16-bit word i2(2)

C     GREAT 3.2.1 format
C
C     Other information
C      gid[0]=-1
C     gdata[0]=module number
C     gdata[1]=information type (2=Pause, 3=Resume, 4=SYNC100, 6=FEE64 disc etc)
C     gdata[2]=information field (bits 16-19)
C     gdata[3]=information field (bits 0-15)
C     gdata[4]=time stamp (bits 16-27)
C     gdata[5]=time stamp (bits 0-15) 
C
C     ADC data format
C
C       gid[0]=channel ident (FEE64 module bits 6-11, channel bits 0-5)
C     gdata[0]=fail/veto bits (AIDA: fail=0, veto=0 LEC/MEC, veto=1 HEC)
C     gdata[1]=?
C     gdata[2]=ADC data
C     gdata[3]=time stamp (bits 16-27)
C     gdata[4]=time stamp (bits 0-15) 

C----67---------------------------------------------------------------72------80
C     gid = - 1   Other information

      IF ( gid(0).EQ.-1 ) THEN

       i2(2) = 0
       i2(1) = gdata(0)
       module = AND( i4, Z'0000003f' )
       i2(1) = gdata(1)
... 1039 more lines ...
Attachment 20: titles.dat
s1 nnaida1: ADC data items (20us/ch)
s2 nnaida2: ADC data items (20us/ch)
s3 nnaida3: ADC data items (20us/ch)
s4 nnaida4: ADC data items (20us/ch)
s5 nnaida5: ADC data items (20us/ch)
s6 nnaida6: ADC data items (20us/ch)
s7 nnaida7: ADC data items (20us/ch)
s8 nnaida8: ADC data items (20us/ch)
s9 nnaida9: ADC data items (20us/ch)
s10 nnaida10: ADC data items (20us/ch)
s11 nnaida11: ADC data items (20us/ch)
s12 nnaida12: ADC data items (20us/ch)
s13 nnaida13: ADC data items (20us/ch)
s14 nnaida14: ADC data items (20us/ch)
s15 nnaida15: ADC data items (20us/ch)
s16 nnaida16: ADC data items (20us/ch)
s17 nnaida17: ADC data items (20us/ch)
s18 nnaida18: ADC data items (20us/ch)
s19 nnaida19: ADC data items (20us/ch)
s20 nnaida20: ADC data items (20us/ch)
s21 nnaida21: ADC data items (20us/ch)
s22 nnaida22: ADC data items (20us/ch)
s23 nnaida23: ADC data items (20us/ch)
s24 nnaida24: ADC data items (20us/ch)
s25 nnaida25: ADC data items (20us/ch)
s26 nnaida26: ADC data items (20us/ch)
s27 nnaida27: ADC data items (20us/ch)
s28 nnaida28: ADC data items (20us/ch)
s29 nnaida29: ADC data items (20us/ch)
s30 nnaida30: ADC data items (20us/ch)
s31 nnaida31: ADC data items (20us/ch)
s32 nnaida32: ADC data items (20us/ch)

s101 nnaida1: Disc data items (20us/ch)
s102 nnaida2: Disc data items (20us/ch)
s103 nnaida3: Disc data items (20us/ch)
s104 nnaida4: Disc data items (20us/ch)
s105 nnaida5: Disc data items (20us/ch)
s106 nnaida6: Disc data items (20us/ch)
s107 nnaida7: Disc data items (20us/ch)
s108 nnaida8: Disc data items (20us/ch)
s109 nnaida9: Disc data items (20us/ch)
s110 nnaida10: Disc data items (20us/ch)
s111 nnaida11: Disc data items (20us/ch)
s112 nnaida12: Disc data items (20us/ch)
s113 nnaida13: Disc data items (20us/ch)
s114 nnaida14: Disc data items (20us/ch)
s115 nnaida15: Disc data items (20us/ch)
s116 nnaida16: Disc data items (20us/ch)
s117 nnaida17: Disc data items (20us/ch)
s118 nnaida18: Disc data items (20us/ch)
s119 nnaida19: Disc data items (20us/ch)
s120 nnaida20: Disc data items (20us/ch)
s121 nnaida21: Disc data items (20us/ch)
s122 nnaida22: Disc data items (20us/ch)
s123 nnaida23: Disc data items (20us/ch)
s124 nnaida24: Disc data items (20us/ch)
s125 nnaida25: Disc data items (20us/ch)
s126 nnaida26: Disc data items (20us/ch)
s127 nnaida27: Disc data items (20us/ch)
s128 nnaida28: Disc data items (20us/ch)
s129 nnaida29: Disc data items (20us/ch)
s130 nnaida30: Disc data items (20us/ch)
s131 nnaida31: Disc data items (20us/ch)
s132 nnaida32: Disc data items (20us/ch)

s200 LEC & HEC channels (DSSSD order)
s201 LEC & HEC channels (ASIC order)

s301 DSSSD #1 LEC rate (Hz) versus time 
s306 DSSSD #6 LEC rate (Hz) versus time
s307 DSSSD #1 HEC rate (Hz) versus time 
s312 DSSSD #6 HEC rate (Hz) versus time

s1000 4*(module # - 1) + asic # (0-95 ADC; 100-195 Disc) v. timestamp (20us/ch)
s1001 HEC channel (2048-3584) versus energy

s2000 m_p versus m_n DSSSD #1 (LEC)
s2005 m_p versus m_n DSSSD #6 (LEC)
s2010 x versus y DSSSD #1 (LEC)
s2015 x versus y DSSSD #6 (LEC)
s2020 E_p versus E_n DSSSD #1 (LEC)
s2025 E_p versus E_n DSSSD #6 (LEC)
s2030 LEC-LEC time DSSSD #1 (2.56us/ch)
s2035 LEC-LEC time DSSSD #6 (2.56us/ch)

s2100 m_p versus m_n DSSSD #1 (HEC)
s2105 m_p versus m_n DSSSD #6 (HEC)
s2110 x versus y DSSSD #1 (HEC)
s2115 x versus y DSSSD #6 (HEC)
s2120 E_p versus E_n DSSSD #1 (HEC)
s2125 E_p versus E_n DSSSD #6 (HEC)

s2130 HEC-HEC time DSSSD #1 (2.56us/ch)
s2135 HEC-HEC time DSSSD #6 (2.56us/ch)

s2200 per pixel HEC-LEC time DSSSD #1 (2.56us/ch)
s2205 per pixel HEC-LEC time DSSSD #6 (2.56us/ch)

s2210 per pixel HEC-LEC time DSSSD #1 (163.840us/ch)
s2215 per pixel HEC-LEC time DSSSD #6 (163.840us/ch)

s2300 HEC E versus HEC-LEC dt DSSSD#1 (163.840us/ch)
s2305 HEC E versus HEC-LEC dt DSSSD#6 (163.840us/ch)

s2400 LEC E versus HEC-LEC dt DSSSD#1 (163.840us/ch)
s2405 LEC E versus HEC-LEC dt DSSSD#6 (163.840us/ch)

s2500 HEC xy versus HEC-LEC dt DSSSD#1 (163.840us/ch)
s2505 HEC xy versus HEC-LEC dt DSSSD#6 (163.840us/ch)


Attachment 21: variables.dat
  $variables
  dsssd_a( 1 ) = 17
  dsssd_a( 2 ) = 18
  dsssd_a( 3 ) = 9
  dsssd_a( 4 ) = 10
  dsssd_a( 5 ) = 1
  dsssd_a( 6 ) = 2
  dsssd_b( 1 ) = 20
  dsssd_b( 2 ) = 19
  dsssd_b( 3 ) = 12
  dsssd_b( 4 ) = 11
  dsssd_b( 5 ) = 4
  dsssd_b( 6 ) = 3
  dsssd_c( 1 ) = 22
  dsssd_c( 2 ) = 21
  dsssd_c( 3 ) = 14
  dsssd_c( 4 ) = 13
  dsssd_c( 5 ) = 6
  dsssd_c( 6 ) = 5
  dsssd_d( 1 ) = 23
  dsssd_d( 2 ) = 24
  dsssd_d( 3 ) = 15
  dsssd_d( 4 ) = 16
  dsssd_d( 5 ) = 7
  dsssd_d( 6 ) = 8
  offset = 4096*0.0
  gain = 4096*1.0
  c = 8*1.0
  w = 32*0
  $[end]
  485   Sun Nov 27 12:53:04 2016 TDOffline analysis R24_3

S O R T    C O M P L E T E ..... Sun Nov 27 20:32:48 2016


 *** ENTRY finish

 *** data items:  261088000 ( 102453.05 Hz)
 *** ADC events:   37054863 (  14540.63 Hz)
 *** time warps:          0 (      0.00 Hz)
 *** DSSSD # 1 count:   5327829 old count:   5309834 dt:    306.37 s  LEC rate:  17389.94 Hz
 *** DSSSD # 2 count:   4754694 old count:   4738448 dt:    306.37 s  LEC rate:  15519.24 Hz
 *** DSSSD # 3 count:   3419847 old count:   3408339 dt:    306.37 s  LEC rate:  11162.32 Hz
 *** DSSSD # 4 count:   3012429 old count:   3002271 dt:    306.37 s  LEC rate:   9832.52 Hz
 *** DSSSD # 5 count:   1976094 old count:   1969554 dt:    306.37 s  LEC rate:   6449.94 Hz
 *** DSSSD # 6 count:   3682932 old count:   3670389 dt:    306.37 s  LEC rate:  12021.03 Hz
 *** DSSSD # 1 count:        37 old count:        37 dt:    306.37 s  HEC rate:      0.12 Hz
 *** DSSSD # 2 count:        49 old count:        48 dt:    306.37 s  HEC rate:      0.16 Hz
 *** DSSSD # 3 count:       174 old count:       172 dt:    306.37 s  HEC rate:      0.57 Hz
 *** DSSSD # 4 count:       228 old count:       226 dt:    306.37 s  HEC rate:      0.74 Hz
 *** DSSSD # 5 count:       210 old count:       209 dt:    306.37 s  HEC rate:      0.69 Hz
 *** DSSSD # 6 count:       125 old count:       125 dt:    306.37 s  HEC rate:      0.41 Hz

 *** ENTRY finish ends

S O R T    S T O P P E D ..... Sun Nov 27 20:34:06 2016
Attachment 1: 20.png
20.png
Attachment 2: 21.png
21.png
Attachment 3: 22.png
22.png
Attachment 4: 24.png
24.png
  771   Fri Nov 9 02:26:58 2018 TDOffline analysis R18, R27, R29 - alpha background
Offline analysis of Oct18/R18_0 to 2, R27_0, 4 & 9, R29_0, 7 & 30

alpha background, BNC PB-4 pulser OFF, LEC/MEC slow comparator 0x64 (~1MeV)

*** ENTRY finish

 *** data items: 2043072278 ( 246610.11 Hz)
 *** ADC events:   10666744 (   1287.53 Hz)
 *** time warps:          2 (      0.00 Hz)
 *** DSSSD # 1 count:      1446 old count:      1446 dt: 231900.06 s  LEC rate:      0.01 Hz
 *** DSSSD # 2 count:       568 old count:       568 dt: 231900.06 s  LEC rate:      0.00 Hz
 *** DSSSD # 3 count:       644 old count:       644 dt: 231900.06 s  LEC rate:      0.00 Hz
 *** DSSSD # 4 count:       612 old count:       612 dt: 231900.06 s  LEC rate:      0.00 Hz
 *** DSSSD # 5 count:       591 old count:       591 dt: 231900.06 s  LEC rate:      0.00 Hz
 *** DSSSD # 6 count:       796 old count:       796 dt: 231900.06 s  LEC rate:      0.00 Hz
 *** DSSSD # 1 count:         6 old count:         6 dt: 231900.06 s  HEC rate:      0.00 Hz
 *** DSSSD # 2 count:         3 old count:         3 dt: 231900.06 s  HEC rate:      0.00 Hz
 *** DSSSD # 3 count:         0 old count:         0 dt: 231900.06 s  HEC rate:      0.00 Hz
 *** DSSSD # 4 count:         1 old count:         1 dt: 231900.06 s  HEC rate:      0.00 Hz
 *** DSSSD # 5 count:         0 old count:         0 dt: 231900.06 s  HEC rate:      0.00 Hz
 *** DSSSD # 6 count:         1 old count:         1 dt: 231900.06 s  HEC rate:      0.00 Hz

 *** ENTRY finish ends

attachment 1 - per DSSSD m_p versus m_n
           2 - per DSSSD x versus y
           3 - per DSSSD E_p versus E_n
           4 - DSSSD channel versus E_p
           5 - DSSSD channel versus E_p (expanded)
           6 - hit pattern (ASIC & DSSSD order)
           7 - titles
           8 - variables 
               N.B. ADC offsets not set, 0 by default
                    | E_p - E_n | < 2000


DSSSD configuration looks OK. Conclude problems identified and addressed November 5 are confirmed
as fixed - see https://elog.ph.ed.ac.uk/AIDA/766 and https://elog.ph.ed.ac.uk/AIDA/761
Attachment 1: 1.png
1.png
Attachment 2: 2.png
2.png
Attachment 3: 3.png
3.png
Attachment 4: 4.png
4.png
Attachment 5: 5.png
5.png
Attachment 6: 6.png
6.png
Attachment 7: titles.dat
s1 nnaida1: ADC data items (20us/ch)
s2 nnaida2: ADC data items (20us/ch)
s3 nnaida3: ADC data items (20us/ch)
s4 nnaida4: ADC data items (20us/ch)
s5 nnaida5: ADC data items (20us/ch)
s6 nnaida6: ADC data items (20us/ch)
s7 nnaida7: ADC data items (20us/ch)
s8 nnaida8: ADC data items (20us/ch)
s9 nnaida9: ADC data items (20us/ch)
s10 nnaida10: ADC data items (20us/ch)
s11 nnaida11: ADC data items (20us/ch)
s12 nnaida12: ADC data items (20us/ch)
s13 nnaida13: ADC data items (20us/ch)
s14 nnaida14: ADC data items (20us/ch)
s15 nnaida15: ADC data items (20us/ch)
s16 nnaida16: ADC data items (20us/ch)
s17 nnaida17: ADC data items (20us/ch)
s18 nnaida18: ADC data items (20us/ch)
s19 nnaida19: ADC data items (20us/ch)
s20 nnaida20: ADC data items (20us/ch)
s21 nnaida21: ADC data items (20us/ch)
s22 nnaida22: ADC data items (20us/ch)
s23 nnaida23: ADC data items (20us/ch)
s24 nnaida24: ADC data items (20us/ch)
s25 nnaida25: ADC data items (20us/ch)
s26 nnaida26: ADC data items (20us/ch)
s27 nnaida27: ADC data items (20us/ch)
s28 nnaida28: ADC data items (20us/ch)
s29 nnaida29: ADC data items (20us/ch)
s30 nnaida30: ADC data items (20us/ch)
s31 nnaida31: ADC data items (20us/ch)
s32 nnaida32: ADC data items (20us/ch)

s101 nnaida1: HEC ADC data items (20us/ch)
s124 nnaida1: HEC ADC data items (20us/ch)

s201 nnaida1: Disc data items (20us/ch)
s202 nnaida2: Disc data items (20us/ch)
s203 nnaida3: Disc data items (20us/ch)
s204 nnaida4: Disc data items (20us/ch)
s205 nnaida5: Disc data items (20us/ch)
s206 nnaida6: Disc data items (20us/ch)
s207 nnaida7: Disc data items (20us/ch)
s208 nnaida8: Disc data items (20us/ch)
s209 nnaida9: Disc data items (20us/ch)
s210 nnaida10: Disc data items (20us/ch)
s211 nnaida11: Disc data items (20us/ch)
s212 nnaida12: Disc data items (20us/ch)
s213 nnaida13: Disc data items (20us/ch)
s214 nnaida14: Disc data items (20us/ch)
s215 nnaida15: Disc data items (20us/ch)
s216 nnaida16: Disc data items (20us/ch)
s217 nnaida17: Disc data items (20us/ch)
s218 nnaida18: Disc data items (20us/ch)
s219 nnaida19: Disc data items (20us/ch)
s220 nnaida20: Disc data items (20us/ch)
s221 nnaida21: Disc data items (20us/ch)
s222 nnaida22: Disc data items (20us/ch)
s223 nnaida23: Disc data items (20us/ch)
s224 nnaida24: Disc data items (20us/ch)
s225 nnaida25: Disc data items (20us/ch)
s226 nnaida26: Disc data items (20us/ch)
s227 nnaida27: Disc data items (20us/ch)
s228 nnaida28: Disc data items (20us/ch)
s229 nnaida29: Disc data items (20us/ch)
s230 nnaida30: Disc data items (20us/ch)
s231 nnaida31: Disc data items (20us/ch)
s232 nnaida32: Disc data items (20us/ch)

s250 LEC & HEC channels (DSSSD order)
s251 LEC & HEC channels (ASIC order)

s301 DSSSD #1 LEC rate (Hz) versus time 
s306 DSSSD #6 LEC rate (Hz) versus time
s307 DSSSD #1 HEC rate (Hz) versus time 
s312 DSSSD #6 HEC rate (Hz) versus time

s1000 4*(module # - 1) + asic # (0-95 ADC; 100-195 Disc) v. timestamp (20us/ch)
s1001 HEC channel (2048-3584) versus energy

s2000 m_p versus m_n DSSSD #1 (LEC)
s2005 m_p versus m_n DSSSD #6 (LEC)
s2010 x versus y DSSSD #1 (LEC)
s2015 x versus y DSSSD #6 (LEC)
s2020 E_p versus E_n DSSSD #1 (LEC)
s2025 E_p versus E_n DSSSD #6 (LEC)
s2030 LEC-LEC time DSSSD #1 (2.56us/ch)
s2035 LEC-LEC time DSSSD #6 (2.56us/ch)

s2100 m_p versus m_n DSSSD #1 (HEC)
s2105 m_p versus m_n DSSSD #6 (HEC)
s2110 x versus y DSSSD #1 (HEC)
s2115 x versus y DSSSD #6 (HEC)
s2120 E_p versus E_n DSSSD #1 (HEC)
s2125 E_p versus E_n DSSSD #6 (HEC)

s2130 HEC-HEC time DSSSD #1 (2.56us/ch)
s2135 HEC-HEC time DSSSD #6 (2.56us/ch)

s2200 per pixel HEC-LEC time DSSSD #1 (2.56us/ch)
s2205 per pixel HEC-LEC time DSSSD #6 (2.56us/ch)

s2210 per pixel HEC-LEC time DSSSD #1 (163.840us/ch)
s2215 per pixel HEC-LEC time DSSSD #6 (163.840us/ch)

s2300 HEC E versus HEC-LEC dt DSSSD#1 (163.840us/ch)
s2305 HEC E versus HEC-LEC dt DSSSD#6 (163.840us/ch)

s2400 LEC E versus HEC-LEC dt DSSSD#1 (163.840us/ch)
s2405 LEC E versus HEC-LEC dt DSSSD#6 (163.840us/ch)

s2500 HEC xy versus HEC-LEC dt DSSSD#1 (163.840us/ch)
s2505 HEC xy versus HEC-LEC dt DSSSD#6 (163.840us/ch)

s2600 LEC x_time - e_time DSSSD#1 (10ns/ch)
s2605 LEC x_time - e_time DSSSD#6 (10ns/ch)

s2700 e_diff versus HEC-LEC dt DSSSD#1 (163.840us/ch)
s2705 e_diff versus HEC-LEC dt DSSSD#6 (163.840us/ch)

s2800 m/z_lec/ctr versus HEC-LEC dt DSSSD#1 (163.840us/ch)
s2805 m/z_lec/ctr versus HEC-LEC dt DSSSD#6 (163.840us/ch)

s2900 LEC-LEC dt (20.48us/ch) versus HEC-LEC dt DSSSD#1 (163.840us/ch)
s2905 LEC-LEC dt (20.48us/ch) versus HEC-LEC dt DSSSD#6 (163.840us/ch)

s3000 FEE64 versus FEE64
Attachment 8: variables.dat
  $variables
  dsssd_a( 1 ) = 17
  dsssd_a( 2 ) = 18
  dsssd_a( 3 ) = 9
  dsssd_a( 4 ) = 10
  dsssd_a( 5 ) = 1
  dsssd_a( 6 ) = 2
  dsssd_b( 1 ) = 20
  dsssd_b( 2 ) = 19
  dsssd_b( 3 ) = 12
  dsssd_b( 4 ) = 11
  dsssd_b( 5 ) = 4
  dsssd_b( 6 ) = 3
  dsssd_c( 1 ) = 22
  dsssd_c( 2 ) = 21
  dsssd_c( 3 ) = 14
  dsssd_c( 4 ) = 13
  dsssd_c( 5 ) = 6
  dsssd_c( 6 ) = 5
  dsssd_d( 1 ) = 23
  dsssd_d( 2 ) = 24
  dsssd_d( 3 ) = 15
  dsssd_d( 4 ) = 16
  dsssd_d( 5 ) = 7
  dsssd_d( 6 ) = 8
  offset = 4096*0.0
  gain = 4096*1.0
  c = 8*1.0
  w = 32*0
  w(1) = 2000
  w(2) = -2000
  w(3) = 210
  w(4) = 210
  $[end]
  766   Mon Nov 5 10:02:52 2018 TDOffline analysis R16 - beam background
Offline analysis of Oct18/R16_0 to 11


BNC PB-4 pulser OFF, LEC/MEC slow comparator 0x64 (~1MeV)
background (presumably high energy, low Z particles) from beam stopping in F11 E detector upstream 
of AIDA and BRIKEN - ~10-15cm Pb shielding installed in BRIKEN to stop heavy opns reaching AIDA

*** ENTRY finish

 *** data items: ********** ( -13850.89 Hz)
 *** ADC events: 2089104652 (  24903.76 Hz)
 *** time warps:          0 (      0.00 Hz)
 *** DSSSD # 1 count:     81173 old count:     81171 dt: 106837.14 s  LEC rate:      0.76 Hz
 *** DSSSD # 2 count:     73671 old count:     73668 dt: 106837.14 s  LEC rate:      0.69 Hz
 *** DSSSD # 3 count:     74114 old count:     74109 dt: 106837.14 s  LEC rate:      0.69 Hz
 *** DSSSD # 4 count:     73324 old count:     73320 dt: 106837.14 s  LEC rate:      0.69 Hz
 *** DSSSD # 5 count:     74312 old count:     74308 dt: 106837.14 s  LEC rate:      0.70 Hz
 *** DSSSD # 6 count:     74653 old count:     74649 dt: 106837.14 s  LEC rate:      0.70 Hz
 *** DSSSD # 1 count:         0 old count:         0 dt: 106837.14 s  HEC rate:      0.00 Hz
 *** DSSSD # 2 count:         0 old count:         0 dt: 106837.14 s  HEC rate:      0.00 Hz
 *** DSSSD # 3 count:         0 old count:         0 dt: 106837.14 s  HEC rate:      0.00 Hz
 *** DSSSD # 4 count:         0 old count:         0 dt: 106837.14 s  HEC rate:      0.00 Hz
 *** DSSSD # 5 count:         0 old count:         0 dt: 106837.14 s  HEC rate:      0.00 Hz
 *** DSSSD # 6 count:         0 old count:         0 dt: 106837.14 s  HEC rate:      0.00 Hz

 *** ENTRY finish ends


attachment 1 - per DSSSD m_p versus m_n
           2 - per DSSSD x versus y
           3 - per DSSSD E_p versus E_n
           4 - DSSSD hit pattern (ASIC & DSSSD order) expanded view of LEC events

DSSSD configuration looks OK. With improved statistics (cf alpha background) some missing strips identifiable
- see also Elog https://elog.ph.ed.ac.uk/AIDA/764
Attachment 1: 200.png
200.png
Attachment 2: 201.png
201.png
Attachment 3: 202.png
202.png
Attachment 4: 203.png
203.png
  761   Sat Nov 3 00:49:22 2018 TDOffline analysis R15 - beam background
Offline analysis of Oct18/R15_0 to 7


BNC PB-4 pulser OFF, LEC/MEC slow comparator 0x64 (~1MeV)
background (presumably high energy, low Z particles) from beam stopping in F11 E detector upstream 
of AIDA and BRIKEN - ~10-15cm Pb shielding installed in BRIKEN to stop heavy opns reaching AIDA

 *** ENTRY finish

 *** data items: 2088704000 (  42315.52 Hz)
 *** ADC events: 1236420713 (  25048.92 Hz)
 *** time warps:          0 (      0.00 Hz)
 *** DSSSD # 1 count:     12557 old count:     12557 dt:  63395.57 s  LEC rate:      0.20 Hz
 *** DSSSD # 2 count:     10332 old count:     10332 dt:  63395.57 s  LEC rate:      0.16 Hz
 *** DSSSD # 3 count:      9644 old count:      9644 dt:  63395.57 s  LEC rate:      0.15 Hz
 *** DSSSD # 4 count:     10024 old count:     10024 dt:  63395.57 s  LEC rate:      0.16 Hz
 *** DSSSD # 5 count:     11035 old count:     11035 dt:  63395.57 s  LEC rate:      0.17 Hz
 *** DSSSD # 6 count:     11684 old count:     11684 dt:  63395.57 s  LEC rate:      0.18 Hz
 *** DSSSD # 1 count:         0 old count:         0 dt:  63395.57 s  HEC rate:      0.00 Hz
 *** DSSSD # 2 count:         0 old count:         0 dt:  63395.57 s  HEC rate:      0.00 Hz
 *** DSSSD # 3 count:         0 old count:         0 dt:  63395.57 s  HEC rate:      0.00 Hz
 *** DSSSD # 4 count:         0 old count:         0 dt:  63395.57 s  HEC rate:      0.00 Hz
 *** DSSSD # 5 count:         0 old count:         0 dt:  63395.57 s  HEC rate:      0.00 Hz
 *** DSSSD # 6 count:         0 old count:         0 dt:  63395.57 s  HEC rate:      0.00 Hz

 *** ENTRY finish ends

attachment 1 - per DSSSD m_p versus m_n
           2 - per DSSSD x versus y
           3 - per DSSSD E_p versus E_n
           4 - DSSSD hit pattern (ASIC & DSSSD order) LEC & HEC events (note HEC events from ASIC check loads)
           5 - DSSSD hit pattern (ASIC & DSSSD order) expanded view of LEC events

DSSSD configuration looks OK. With improved statistics (cf alpha background) some missing
strips identifiable.
Attachment 1: 50.png
50.png
Attachment 2: 51.png
51.png
Attachment 3: 52.png
52.png
Attachment 4: 54.png
54.png
Attachment 5: 55.png
55.png
  717   Thu Jun 21 05:23:16 2018 TD, OHOffline analysis R1 - alpha background
Analysis of overnight alpha background run (June18/R1_0 ... _2) indicates that all FEE64s and DSSSDs are
operating OK.

Spectra for DSSSD #1-6 - see attachments 1-5

Sort program, variables and spectra titles - see attachments 6-7

Sort program statistics

 *** ENTRY finish

 *** data items:  671431805 ( 363493.31 Hz)
 *** ADC events:     743540 (    402.53 Hz)
 *** time warps:          0 (      0.00 Hz)
 *** DSSSD # 1 count:       266 old count:       266 dt:  73278.38 s  LEC rate:      0.00 Hz
 *** DSSSD # 2 count:       119 old count:       119 dt:  73278.38 s  LEC rate:      0.00 Hz
 *** DSSSD # 3 count:       131 old count:       131 dt:  73278.38 s  LEC rate:      0.00 Hz
 *** DSSSD # 4 count:        97 old count:        97 dt:  73278.38 s  LEC rate:      0.00 Hz
 *** DSSSD # 5 count:       103 old count:       103 dt:  73278.38 s  LEC rate:      0.00 Hz
 *** DSSSD # 6 count:       122 old count:       122 dt:  73278.38 s  LEC rate:      0.00 Hz
 *** DSSSD # 1 count:         0 old count:         0 dt:  73278.38 s  HEC rate:      0.00 Hz
 *** DSSSD # 2 count:         0 old count:         0 dt:  73278.38 s  HEC rate:      0.00 Hz
 *** DSSSD # 3 count:         0 old count:         0 dt:  73278.38 s  HEC rate:      0.00 Hz
 *** DSSSD # 4 count:         0 old count:         0 dt:  73278.38 s  HEC rate:      0.00 Hz
 *** DSSSD # 5 count:         0 old count:         0 dt:  73278.38 s  HEC rate:      0.00 Hz
 *** DSSSD # 6 count:         0 old count:         0 dt:  73278.38 s  HEC rate:      0.00 Hz

 *** ENTRY finish ends
Attachment 1: 40.png
40.png
Attachment 2: 41.png
41.png
Attachment 3: 42.png
42.png
Attachment 4: 43.png
43.png
Attachment 5: 45.png
45.png
Attachment 6: implantdecay4.f
*trigger
128  
*oned
1..32       s 65536
101..132    s 65536
201..232    s 65536
250..254    s 4096
301..332    s 65536
400         s 64
500..511    s 4096
2130..2135  s 65536
2030..2035  s 65536
2200..2205  s 65536
2210..2215  s 65536
*twod
1000..1001  s 4096 4096
1100..1105  s 512 512
2000..2005  s 128 128
2010..2015  s 128 128
2020..2025  s 512 512
2100..2105  s 128 128
2110..2115  s 128 128
2120..2125  s 512 512
2300..2305  s 1024 1024
2310..2315  s 1024 1024
2400..2405  s 1024 1024
2410..2415  s 1024 1024
2500..2505  s 1024 1024
2510..2515  s 1024 1024
2600..2605  s 256 256
2700..2705  s 1024 1024
2710..2715  s 1024 1024
2800..2805  s 1024 1024
2810..2815  s 1024 1024
2900..2905  s 1024 1024
2910..2915  s 1024 1024
3000..3005  s 1024 1024
3010..3015  s 1024 1024
3100..3105  s 1024 1024
4000        s 64 64
*vars
*sort
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, channel, channel_ident
      INTEGER count, ctr, data(0:4095), disc
      INTEGER dsssd_a( 6 ), dsssd_b( 6 ), dsssd_c( 6 ), dsssd_d( 6 )
      INTEGER*8 dt, dt1, dt2
      INTEGER ediff, events, events_old, fee(32)
      INTEGER hec_e_sum_x, hec_e_sum_y
      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 idata(0:4095), ierr, information, information_index
      INTEGER invalid_id, items, items_old, j, j1, j1_old, j2, j2_old, k
      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 next_ch, next_adc_data, next_module
      INTEGER sc_channel
      INTEGER scaler( 32 ), scaler_old( 32 ), range, time_warp, total
      INTEGER*8 ts
      INTEGER ts28
      INTEGER ts48
      INTEGER*8 old_ts( 8, 2 ), ts_old, next_ts
      INTEGER implant_e( 0:127, 0:127, 1:8 )
      INTEGER*8 implant_ts( 0:127, 0:127, 1:8 )
      INTEGER*8 decay_ts( 0:127, 0:127, 1:8 )
      INTEGER*8 old_decay_ts( 0:127, 0:127, 1:8 )
      INTEGER*8 older_decay_ts( 0:127, 0:127, 1:8 )
      INTEGER type, w(0:31), zzz
      INTEGER z_lec, z_lec_veto, z_hec
      INTEGER*8 hit_ts(0:4095), ihit_ts(0:4095), tdiff1, tdiff2
       
C
      LOGICAL first, first_ts, hit(0:4095), ihit(0:4095)
      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

C----67---------------------------------------------------------------72------80
C
C     Namelists
C
      NAMELIST /variables/ gain, offset, c, w,
     +                     dsssd_a, dsssd_b, dsssd_c, dsssd_d

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, * ) ' *** AIDA @ RIKEN - implant-decay - June 2017'
      WRITE ( 6, * ) ' *** Entry init commences'

C     Initialise counters

      events = 0
      events_old = 0

      time_warp = 0
      invalid_id = 0

      first_ts = .TRUE.
      t3 = 0.0D+00
      t4 = 0.0D+00
      t4_old = 0.0D+00

C     Initialise event data

      DO i = 0, 4095
       data( i ) = 0
       hit( i ) = .FALSE.
       hit_ts( i ) = 0
       ihit_ts( i ) = 0
      ENDDO
      DO i = 1, 32
       fee(i) = 0
      ENDDO

      first = .true.
      ts_old = 0
      e_time_old = 0
      count = 0
      total = 0
      ctr = 0

      DO i = 0, 127
       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, 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

C     Initialise time

      t1 = SECNDS(0.0)

C     Read program variables via NAMELIST I/O

      OPEN( 1, FILE = '/homes/npg/td/GREAT/variables.dat',
     +         IOSTAT = ierr )
      IF ( ierr.NE.0 ) THEN
       WRITE( 6, * ) ' *** OPEN I/O error:', ierr
       RETURN
      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
       RETURN
      ENDIF  

C     Display program variables

      DO i = 0, 4095, 8
       WRITE( 6, 9001 ) i, i+7, ( gain( j ), j = i, i+7 )
      ENDDO

      DO i = 0, 4095, 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, 6
       WRITE( 6, 9005 ) i, dsssd_a(i), dsssd_b(i),
     +                     dsssd_c(i), dsssd_d(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
       WRITE( 6, 9010 ) items, FLOAT( items - items_old ) / delta_t,
     +                  events, FLOAT( events - events_old ) / delta_t
       items_old = items
       events_old = events
       t2_old = t2

      ENDIF

      IF ( MOD( items, 50000 ).EQ.0 ) THEN

       t4 = DFLOAT( ts ) * 10.0D-9 - t3
       delta_t = t4 - t4_old

       IF ( delta_t.GE.1.0D+00 ) THEN

        DO i = 1, 32
         rates( i ) = DFLOAT( scaler( i ) - scaler_old( i ) ) / delta_t
        ENDDO
C        WRITE( 6, * ) ' '
        DO i = 1, 6
C         WRITE( 6, 9020 ) i, scaler( i ), scaler_old( i ),
C     +                    delta_t, rates( i )
         sc_channel = INT ( MOD( t4, 65536.0 ) + 0.5 )
         CALL set1d( i+300, sc_channel, INT( rates( i ) + 0.5 ) )
        ENDDO
        DO i = 7, 12
C         WRITE( 6, 9030 ) i-6, scaler( i ), scaler_old( i ),
C     +                    delta_t, rates( i )
         sc_channel = INT ( MOD( t4, 65536.0 ) + 0.5 )
         CALL set1d( i+300, sc_channel, INT( rates( i ) + 0.5 ) )
        ENDDO
        DO i = 1, 32
         scaler_old( i ) = scaler( i )
        ENDDO
        t4_old = t4
       ENDIF

      ENDIF

C----67---------------------------------------------------------------72------80
C     Event data

C     For Linux/x64 & Windows/x86 least significant 16-bit word i2(1)
C     For Solaris/SPARC least significant 16-bit word i2(2)

C     GREAT 3.2.1 format
... 1616 more lines ...
Attachment 7: variables.dat
  $variables
  dsssd_a( 1 ) = 1
  dsssd_a( 2 ) = 2
  dsssd_a( 3 ) = 9
  dsssd_a( 4 ) = 10
  dsssd_a( 5 ) = 17
  dsssd_a( 6 ) = 18
  dsssd_b( 1 ) = 20
  dsssd_b( 2 ) = 19
  dsssd_b( 3 ) = 12
  dsssd_b( 4 ) = 11
  dsssd_b( 5 ) = 4
  dsssd_b( 6 ) = 3
  dsssd_c( 1 ) = 22
  dsssd_c( 2 ) = 21
  dsssd_c( 3 ) = 14
  dsssd_c( 4 ) = 13
  dsssd_c( 5 ) = 6
  dsssd_c( 6 ) = 5
  dsssd_d( 1 ) = 23
  dsssd_d( 2 ) = 24
  dsssd_d( 3 ) = 15
  dsssd_d( 4 ) = 16
  dsssd_d( 5 ) = 7
  dsssd_d( 6 ) = 8
  offset = 4096*0.0
  gain = 4096*1.0
  c = 8*1.0
  w = 32*0
  offset(   0) =    16.10
  offset(   1) =   -14.93
  offset(   2) =    -4.03
  offset(   3) =    24.20
  offset(   4) =    -1.49
  offset(   5) =    18.78
  offset(   6) =     2.57
  offset(   7) =    -0.06
  offset(   8) =     8.52
  offset(   9) =    17.08
  offset(  10) =    -8.18
  offset(  11) =    11.06
  offset(  12) =     1.45
  offset(  13) =    -0.11
  offset(  14) =    -5.27
  offset(  15) =    12.08
  offset(  16) =   -13.45
  offset(  17) =   -23.82
  offset(  18) =   -18.33
  offset(  19) =     9.29
  offset(  20) =     1.15
  offset(  21) =   -20.83
  offset(  22) =   -14.00
  offset(  23) =     0.04
  offset(  24) =    -0.18
  offset(  25) =    -0.13
  offset(  26) =    -7.15
  offset(  27) =    -0.09
  offset(  28) =     8.33
  offset(  29) =   -16.03
  offset(  30) =     9.93
  offset(  31) =    -0.08
  offset(  32) =    -8.71
  offset(  33) =   -28.08
  offset(  34) =    -1.57
  offset(  35) =     4.32
  offset(  36) =    -0.98
  offset(  37) =     9.03
  offset(  38) =    10.18
  offset(  39) =    -8.64
  offset(  40) =    15.82
  offset(  41) = -9999.99
  offset(  42) =    25.30
  offset(  43) =    -1.59
  offset(  44) =     1.77
  offset(  45) =    10.58
  offset(  46) =    16.35
  offset(  47) =     3.63
  offset(  48) =    -3.42
  offset(  49) =   -14.59
  offset(  50) =    35.18
  offset(  51) =     4.42
  offset(  52) =    12.15
  offset(  53) = -9999.99
  offset(  54) = -9999.99
  offset(  55) = -9999.99
  offset(  56) =    15.79
  offset(  57) = -9999.99
  offset(  58) = -9999.99
  offset(  59) =    37.06
  offset(  60) =     0.63
  offset(  61) =    22.36
  offset(  62) = -9999.99
  offset(  63) =    11.76
  offset(  64) =    26.68
  offset(  65) =   321.40
  offset(  66) =   -16.48
  offset(  67) =    11.08
  offset(  68) =    13.34
  offset(  69) =    10.44
  offset(  70) =    -1.68
  offset(  71) =    52.27
  offset(  72) =    16.87
  offset(  73) = -9999.99
  offset(  74) =     0.25
  offset(  75) =    11.73
  offset(  76) =    -4.28
  offset(  77) =    26.30
  offset(  78) =   -14.03
  offset(  79) =   -19.91
  offset(  80) =     3.36
  offset(  81) =    -9.93
  offset(  82) =    -7.04
  offset(  83) =    14.20
  offset(  84) =     4.07
  offset(  85) =   -33.64
  offset(  86) =     1.75
  offset(  87) =    -8.25
  offset(  88) =     0.12
  offset(  89) = -9999.99
  offset(  90) =    -3.26
  offset(  91) =    14.34
  offset(  92) =   -24.36
  offset(  93) =    -8.27
  offset(  94) =    -0.30
  offset(  95) =     5.13
  offset(  96) =    47.52
  offset(  97) =    29.22
  offset(  98) =    10.99
  offset(  99) =    38.34
  offset( 100) =    25.33
  offset( 101) =     9.46
  offset( 102) =    29.87
  offset( 103) =    31.58
  offset( 104) =     8.87
  offset( 105) =    32.96
  offset( 106) =    -8.22
  offset( 107) =    12.51
  offset( 108) =    28.56
  offset( 109) =    34.60
  offset( 110) =     9.83
  offset( 111) =    50.37
  offset( 112) =     1.92
  offset( 113) =   102.85
  offset( 114) = -9999.99
  offset( 115) =   -39.65
  offset( 116) =   -40.48
  offset( 117) = -9999.99
  offset( 118) = -9999.99
  offset( 119) = -9999.99
  offset( 120) = -9999.99
  offset( 121) =   -31.23
  offset( 122) =   -34.30
  offset( 123) =  -165.01
  offset( 124) = -9999.99
  offset( 125) =   -25.43
  offset( 126) =   -12.98
  offset( 127) =   -24.50
  offset( 128) = -9999.99
  offset( 129) =   -30.13
  offset( 130) =   -14.02
  offset( 131) =     8.01
  offset( 132) =    -2.04
  offset( 133) =     5.37
  offset( 134) =   -13.81
  offset( 135) =   -32.82
  offset( 136) =    -2.94
  offset( 137) =   -19.71
  offset( 138) =   -15.45
  offset( 139) =    28.13
  offset( 140) =   -10.23
  offset( 141) =   -19.78
  offset( 142) =     6.51
  offset( 143) =   -22.02
  offset( 144) =   -35.77
  offset( 145) =   -20.93
  offset( 146) =    -8.63
  offset( 147) =   -17.27
  offset( 148) =    -4.88
  offset( 149) =     2.34
  offset( 150) =     9.69
  offset( 151) =   124.36
  offset( 152) =   -10.54
  offset( 153) =   -32.31
  offset( 154) =    -1.74
  offset( 155) =     6.40
  offset( 156) =   -41.53
  offset( 157) =   -12.64
  offset( 158) = -9999.99
  offset( 159) =    -2.22
  offset( 160) =   -17.31
  offset( 161) = -9999.99
  offset( 162) =     5.38
  offset( 163) =    -1.76
  offset( 164) =   -18.30
  offset( 165) =     9.68
  offset( 166) =     6.79
  offset( 167) =  -126.74
  offset( 168) =    -9.77
  offset( 169) =     8.47
  offset( 170) =     7.11
  offset( 171) =    28.72
  offset( 172) =   -29.96
  offset( 173) =   -27.75
  offset( 174) =    -2.43
  offset( 175) =    11.08
  offset( 176) =     7.05
  offset( 177) =   -27.55
  offset( 178) =    -1.09
  offset( 179) =     5.26
  offset( 180) =   -12.73
  offset( 181) =   -15.53
  offset( 182) =    -6.55
  offset( 183) =    34.68
  offset( 184) =    19.55
  offset( 185) =   -14.17
  offset( 186) =    -8.89
  offset( 187) =     8.12
  offset( 188) =    -9.06
  offset( 189) =     8.78
  offset( 190) =   -31.35
  offset( 191) =     5.97
  offset( 192) =   -13.90
  offset( 193) =   -10.36
  offset( 194) =    18.66
  offset( 195) =    -9.45
  offset( 196) =     8.94
  offset( 197) =     8.83
  offset( 198) =    -0.29
  offset( 199) =   -20.45
  offset( 200) =    -8.42
  offset( 201) =    -7.99
  offset( 202) =   -12.23
  offset( 203) =    -2.14
  offset( 204) =    -3.89
  offset( 205) =    -0.34
  offset( 206) =    25.49
  offset( 207) =    11.81
  offset( 208) =     7.79
  offset( 209) =    -2.57
  offset( 210) =    29.34
  offset( 211) =    -5.80
  offset( 212) =    14.64
  offset( 213) =    -2.74
  offset( 214) =    15.44
  offset( 215) =   -13.24
  offset( 216) =     3.93
  offset( 217) =     8.62
  offset( 218) =   -15.60
  offset( 219) =     6.02
  offset( 220) =    20.74
  offset( 221) =    17.48
  offset( 222) =   -23.06
  offset( 223) =    -3.13
  offset( 224) =    30.14
  offset( 225) =     0.82
  offset( 226) =     0.93
  offset( 227) =    11.01
  offset( 228) =    18.31
  offset( 229) =    26.80
  offset( 230) =    28.42
  offset( 231) =     8.85
  offset( 232) =    33.87
  offset( 233) =   -22.41
  offset( 234) =     9.61
  offset( 235) =   -21.32
  offset( 236) =    13.79
  offset( 237) =    16.58
  offset( 238) =    24.98
  offset( 239) =   -12.59
  offset( 240) =   -10.85
  offset( 241) =    -4.30
  offset( 242) =    15.44
  offset( 243) =     6.72
  offset( 244) =   -19.95
  offset( 245) =   -11.12
  offset( 246) =   -15.90
  offset( 247) =   -17.53
  offset( 248) =   -15.94
  offset( 249) =   -21.35
  offset( 250) =     5.89
  offset( 251) =     7.21
  offset( 252) =   -17.58
  offset( 253) =     4.70
  offset( 254) =     8.84
  offset( 255) =   -36.92
  offset( 256) =    18.26
  offset( 257) =   -13.97
  offset( 258) =     3.43
  offset( 259) =   -16.64
  offset( 260) =    23.27
  offset( 261) =    11.47
  offset( 262) =    -3.04
  offset( 263) =    11.25
  offset( 264) =    18.61
  offset( 265) =    -5.86
  offset( 266) =     1.15
  offset( 267) =    -1.39
  offset( 268) =     2.92
  offset( 269) =     7.21
  offset( 270) =    -0.52
... 1271 more lines ...
Attachment 8: titles.dat
s1 nnaida1: ADC data items (20us/ch)
s2 nnaida2: ADC data items (20us/ch)
s3 nnaida3: ADC data items (20us/ch)
s4 nnaida4: ADC data items (20us/ch)
s5 nnaida5: ADC data items (20us/ch)
s6 nnaida6: ADC data items (20us/ch)
s7 nnaida7: ADC data items (20us/ch)
s8 nnaida8: ADC data items (20us/ch)
s9 nnaida9: ADC data items (20us/ch)
s10 nnaida10: ADC data items (20us/ch)
s11 nnaida11: ADC data items (20us/ch)
s12 nnaida12: ADC data items (20us/ch)
s13 nnaida13: ADC data items (20us/ch)
s14 nnaida14: ADC data items (20us/ch)
s15 nnaida15: ADC data items (20us/ch)
s16 nnaida16: ADC data items (20us/ch)
s17 nnaida17: ADC data items (20us/ch)
s18 nnaida18: ADC data items (20us/ch)
s19 nnaida19: ADC data items (20us/ch)
s20 nnaida20: ADC data items (20us/ch)
s21 nnaida21: ADC data items (20us/ch)
s22 nnaida22: ADC data items (20us/ch)
s23 nnaida23: ADC data items (20us/ch)
s24 nnaida24: ADC data items (20us/ch)
s25 nnaida25: ADC data items (20us/ch)
s26 nnaida26: ADC data items (20us/ch)
s27 nnaida27: ADC data items (20us/ch)
s28 nnaida28: ADC data items (20us/ch)
s29 nnaida29: ADC data items (20us/ch)
s30 nnaida30: ADC data items (20us/ch)
s31 nnaida31: ADC data items (20us/ch)
s32 nnaida32: ADC data items (20us/ch)

s101 nnaida1: HEC ADC data items (20us/ch)
s124 nnaida1: HEC ADC data items (20us/ch)

s201 nnaida1: Disc data items (20us/ch)
s202 nnaida2: Disc data items (20us/ch)
s203 nnaida3: Disc data items (20us/ch)
s204 nnaida4: Disc data items (20us/ch)
s205 nnaida5: Disc data items (20us/ch)
s206 nnaida6: Disc data items (20us/ch)
s207 nnaida7: Disc data items (20us/ch)
s208 nnaida8: Disc data items (20us/ch)
s209 nnaida9: Disc data items (20us/ch)
s210 nnaida10: Disc data items (20us/ch)
s211 nnaida11: Disc data items (20us/ch)
s212 nnaida12: Disc data items (20us/ch)
s213 nnaida13: Disc data items (20us/ch)
s214 nnaida14: Disc data items (20us/ch)
s215 nnaida15: Disc data items (20us/ch)
s216 nnaida16: Disc data items (20us/ch)
s217 nnaida17: Disc data items (20us/ch)
s218 nnaida18: Disc data items (20us/ch)
s219 nnaida19: Disc data items (20us/ch)
s220 nnaida20: Disc data items (20us/ch)
s221 nnaida21: Disc data items (20us/ch)
s222 nnaida22: Disc data items (20us/ch)
s223 nnaida23: Disc data items (20us/ch)
s224 nnaida24: Disc data items (20us/ch)
s225 nnaida25: Disc data items (20us/ch)
s226 nnaida26: Disc data items (20us/ch)
s227 nnaida27: Disc data items (20us/ch)
s228 nnaida28: Disc data items (20us/ch)
s229 nnaida29: Disc data items (20us/ch)
s230 nnaida30: Disc data items (20us/ch)
s231 nnaida31: Disc data items (20us/ch)
s232 nnaida32: Disc data items (20us/ch)

s250 LEC & HEC channels (DSSSD order)
s251 LEC & HEC channels (ASIC order)

s301 DSSSD #1 LEC rate (Hz) versus time 
s306 DSSSD #6 LEC rate (Hz) versus time
s307 DSSSD #1 HEC rate (Hz) versus time 
s312 DSSSD #6 HEC rate (Hz) versus time

s1000 4*(module # - 1) + asic # (0-95 ADC; 100-195 Disc) v. timestamp (20us/ch)
s1001 HEC channel (2048-3584) versus energy

s2000 m_p versus m_n DSSSD #1 (LEC)
s2005 m_p versus m_n DSSSD #6 (LEC)
s2010 x versus y DSSSD #1 (LEC)
s2015 x versus y DSSSD #6 (LEC)
s2020 E_p versus E_n DSSSD #1 (LEC)
s2025 E_p versus E_n DSSSD #6 (LEC)
s2030 LEC-LEC time DSSSD #1 (2.56us/ch)
s2035 LEC-LEC time DSSSD #6 (2.56us/ch)

s2100 m_p versus m_n DSSSD #1 (HEC)
s2105 m_p versus m_n DSSSD #6 (HEC)
s2110 x versus y DSSSD #1 (HEC)
s2115 x versus y DSSSD #6 (HEC)
s2120 E_p versus E_n DSSSD #1 (HEC)
s2125 E_p versus E_n DSSSD #6 (HEC)

s2130 HEC-HEC time DSSSD #1 (2.56us/ch)
s2135 HEC-HEC time DSSSD #6 (2.56us/ch)

s2200 per pixel HEC-LEC time DSSSD #1 (2.56us/ch)
s2205 per pixel HEC-LEC time DSSSD #6 (2.56us/ch)

s2210 per pixel HEC-LEC time DSSSD #1 (163.840us/ch)
s2215 per pixel HEC-LEC time DSSSD #6 (163.840us/ch)

s2300 HEC E versus HEC-LEC dt DSSSD#1 (163.840us/ch)
s2305 HEC E versus HEC-LEC dt DSSSD#6 (163.840us/ch)

s2400 LEC E versus HEC-LEC dt DSSSD#1 (163.840us/ch)
s2405 LEC E versus HEC-LEC dt DSSSD#6 (163.840us/ch)

s2500 HEC xy versus HEC-LEC dt DSSSD#1 (163.840us/ch)
s2505 HEC xy versus HEC-LEC dt DSSSD#6 (163.840us/ch)

s2600 LEC x_time - e_time DSSSD#1 (10ns/ch)
s2605 LEC x_time - e_time DSSSD#6 (10ns/ch)

s2700 e_diff versus HEC-LEC dt DSSSD#1 (163.840us/ch)
s2705 e_diff versus HEC-LEC dt DSSSD#6 (163.840us/ch)

s2800 m/z_lec/ctr versus HEC-LEC dt DSSSD#1 (163.840us/ch)
s2805 m/z_lec/ctr versus HEC-LEC dt DSSSD#6 (163.840us/ch)

s2900 LEC-LEC dt (20.48us/ch) versus HEC-LEC dt DSSSD#1 (163.840us/ch)
s2905 LEC-LEC dt (20.48us/ch) versus HEC-LEC dt DSSSD#6 (163.840us/ch)

s3000 FEE64 versus FEE64
  686   Fri Oct 27 08:35:33 2017 TDOffline analysis Oct17/R17_31-R17_34
Offline analysis of High Energy Channel (HEC) 20GeV FSR data from AIDA data file Oct17/R17_31-R17_34

Degrader OUT

Attachment 1 - per DSSSD #1-6 m_p versus m_n
Attachment 2 - 2D spectra contour levels/colours
Attachment 3 - per DSSSD #1-6 x versus y 
Attachment 4 - per DSSSD #1-6 E_p versus E_n (20MeV/channel nominal)

Attachments 5-7 - as above but with all statistics from AIDA data files Oct17/R17_34

Attachment 8 - Sum of E_p energies for DSSSD #1-6 (5.6MeV/channel nominal)

               Peak at channel 1450 corresponds to an energy of 8.1GeV
               End of spectrum channel 1850 corresponds to 10.4GeV

               Note - HEC data not gain-matched
                      Using nominal HEC calibration accurate to ~ +/-20%

               
Attachment 1: 51.png
51.png
Attachment 2: 52.png
52.png
Attachment 3: 54.png
54.png
Attachment 4: 55.png
55.png
Attachment 5: 10.png
10.png
Attachment 6: 11.png
11.png
Attachment 7: 12.png
12.png
Attachment 8: 5.png
5.png
  688   Sat Oct 28 06:09:44 2017 TDOffline analysis Oct17/R17_31-R17_34
Attachment 1: 20.png
20.png
Attachment 2: 21.png
21.png
Attachment 3: 22.png
22.png
Attachment 4: 23.png
23.png
Attachment 5: 24.png
24.png
Attachment 6: 25.png
25.png
Attachment 7: 27.png
27.png
Attachment 8: 28.png
28.png
Attachment 9: 29.png
29.png
Attachment 10: 30.png
30.png
Attachment 11: 31.png
31.png
Attachment 12: 32.png
32.png
Attachment 13: 33.png
33.png
Attachment 14: 34.png
34.png
Attachment 15: 35.png
35.png
Attachment 16: 36.png
36.png
Attachment 17: 37.png
37.png
Attachment 18: 38.png
38.png
Attachment 19: 39.png
39.png
Attachment 20: 40.png
40.png
  702   Sat Nov 18 03:53:43 2017 TDOffline analysis Nov17/R2_8-R2_16 data
MIDASsort offline analysis Nov17/R2_8-R2_16 data

N.B. 
1) no software thresholds
2) no PID correlation

n+n FEE64 nnaida15 Ibias LF feedback 0xf

Implantation rate (DSSSD # 1-6) c. 15Hz
Elapsed time 3441s

 *** ENTRY finish

 *** data items: 2088746796 (  80001.28 Hz)
 *** ADC events:  363022502 (  13904.16 Hz)
 *** time warps:          0 (      0.00 Hz)
 *** DSSSD # 1 count:  15477751 old count:  15473939 dt:   3441.15 s  LEC rate:   4497.84 Hz
 *** DSSSD # 2 count:  12001205 old count:  11998220 dt:   3441.15 s  LEC rate:   3487.55 Hz
 *** DSSSD # 3 count:   5162182 old count:   5160784 dt:   3441.15 s  LEC rate:   1500.13 Hz
 *** DSSSD # 4 count:  11035612 old count:  11032823 dt:   3441.15 s  LEC rate:   3206.95 Hz
 *** DSSSD # 5 count:    395233 old count:    395139 dt:   3441.15 s  LEC rate:    114.85 Hz
 *** DSSSD # 6 count:    168148 old count:    168117 dt:   3441.15 s  LEC rate:     48.86 Hz
 *** DSSSD # 1 count:      1087 old count:      1086 dt:   3441.15 s  HEC rate:      0.32 Hz
 *** DSSSD # 2 count:      4233 old count:      4233 dt:   3441.15 s  HEC rate:      1.23 Hz
 *** DSSSD # 3 count:     12370 old count:     12364 dt:   3441.15 s  HEC rate:      3.59 Hz
 *** DSSSD # 4 count:     24528 old count:     24527 dt:   3441.15 s  HEC rate:      7.13 Hz
 *** DSSSD # 5 count:      7256 old count:      7252 dt:   3441.15 s  HEC rate:      2.11 Hz
 *** DSSSD # 6 count:       509 old count:       509 dt:   3441.15 s  HEC rate:      0.15 Hz

 *** ENTRY finish ends
Attachment 1: 110.png
110.png
Attachment 2: 111.png
111.png
Attachment 3: 112.png
112.png
Attachment 4: 113.png
113.png
Attachment 5: 114.png
114.png
Attachment 6: 115.png
115.png
Attachment 7: 116.png
116.png
Attachment 8: 117.png
117.png
Attachment 9: 118.png
118.png
Attachment 10: 119.png
119.png
Attachment 11: 120.png
120.png
Attachment 12: 121.png
121.png
Attachment 13: 122.png
122.png
Attachment 14: 123.png
123.png
Attachment 15: 124.png
124.png
Attachment 16: 125.png
125.png
Attachment 17: 126.png
126.png
Attachment 18: 127.png
127.png
Attachment 19: 128.png
128.png
Attachment 20: 129.png
129.png
Attachment 21: 130.png
130.png
Attachment 22: 131.png
131.png
  811   Fri May 24 07:21:31 2019 TDOffline analysis May19/R1, R2, R7, R8 & R9
Offline analysis of alpha background data files May19/R1, R2, R7, R8, R9 ( c. 310h )

1  per DSSSD m_p versus m_n
2  per DSSSD x versus y
10 per DSSSD E_p versus E_n
3  2D contours
4  hit pattern ( DSSSD order ) versus E
5  hit pattern ( DSSSD order)
6  hit pattern ( ASIC order )
7  sort program
8  sort variables 
9  sort spectra titles
11 x & y projections of 5 & 7MeV groups from hit pattern ( DSSSD order ) versus E
Attachment 1: 100.png
100.png
Attachment 2: 101.png
101.png
Attachment 3: 103.png
103.png
Attachment 4: 104.png
104.png
Attachment 5: 105.png
105.png
Attachment 6: 106.png
106.png
Attachment 7: implantdecay4.f
*trigger
128  
*oned
1..32       s 65536
101..132    s 65536
201..232    s 65536
250..254    s 4096
301..332    s 65536
400         s 64
500..511    s 4096
2130..2135  s 65536
2030..2035  s 65536
2200..2205  s 65536
2210..2215  s 65536
*twod
1000..1001  s 4096 4096
1100..1105  s 512 512
2000..2005  s 128 128
2010..2015  s 128 128
2020..2025  s 512 512
2100..2105  s 128 128
2110..2115  s 128 128
2120..2125  s 512 512
2300..2305  s 1024 1024
2310..2315  s 1024 1024
2400..2405  s 1024 1024
2410..2415  s 1024 1024
2500..2505  s 1024 1024
2510..2515  s 1024 1024
2600..2605  s 256 256
2700..2705  s 1024 1024
2710..2715  s 1024 1024
2800..2805  s 1024 1024
2810..2815  s 1024 1024
2900..2905  s 1024 1024
2910..2915  s 1024 1024
3000..3005  s 1024 1024
3010..3015  s 1024 1024
3100..3105  s 1024 1024
4000        s 64 64
*vars
*sort
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, channel, channel_ident
      INTEGER count, ctr, data(0:4095), disc
      INTEGER dsssd_a( 6 ), dsssd_b( 6 ), dsssd_c( 6 ), dsssd_d( 6 )
      INTEGER*8 dt, dt1, dt2
      INTEGER ediff, events, events_old, fee(32)
      INTEGER hec_e_sum_x, hec_e_sum_y
      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 idata(0:4095), ierr, information, information_index
      INTEGER invalid_id, items, items_old, j, j1, j1_old, j2, j2_old, k
      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 next_ch, next_adc_data, next_module
      INTEGER sc_channel
      INTEGER scaler( 32 ), scaler_old( 32 ), range, time_warp, total
      INTEGER*8 ts
      INTEGER ts28
      INTEGER ts48
      INTEGER*8 old_ts( 8, 2 ), ts_old, next_ts
      INTEGER implant_e( 0:127, 0:127, 1:8 )
      INTEGER*8 implant_ts( 0:127, 0:127, 1:8 )
      INTEGER*8 decay_ts( 0:127, 0:127, 1:8 )
      INTEGER*8 old_decay_ts( 0:127, 0:127, 1:8 )
      INTEGER*8 older_decay_ts( 0:127, 0:127, 1:8 )
      INTEGER type, w(0:31), zzz
      INTEGER z_lec, z_lec_veto, z_hec
      INTEGER*8 hit_ts(0:4095), ihit_ts(0:4095), tdiff1, tdiff2
       
C
      LOGICAL first, first_ts, hit(0:4095), ihit(0:4095)
      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

C----67---------------------------------------------------------------72------80
C
C     Namelists
C
      NAMELIST /variables/ gain, offset, c, w,
     +                     dsssd_a, dsssd_b, dsssd_c, dsssd_d

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, * ) ' *** AIDA @ RIKEN - implant-decay - June 2017'
      WRITE ( 6, * ) ' *** Entry init commences'

C     Initialise counters

      events = 0
      events_old = 0

      time_warp = 0
      invalid_id = 0

      first_ts = .TRUE.
      t3 = 0.0D+00
      t4 = 0.0D+00
      t4_old = 0.0D+00

C     Initialise event data

      DO i = 0, 4095
       data( i ) = 0
       hit( i ) = .FALSE.
       hit_ts( i ) = 0
       ihit_ts( i ) = 0
      ENDDO
      DO i = 1, 32
       fee(i) = 0
      ENDDO

      first = .true.
      ts_old = 0
      e_time_old = 0
      count = 0
      total = 0
      ctr = 0

      DO i = 0, 127
       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, 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

C     Initialise time

      t1 = SECNDS(0.0)

C     Read program variables via NAMELIST I/O

      OPEN( 1, FILE = '/homes/npg/td/GREAT/variables.dat',
     +         IOSTAT = ierr )
      IF ( ierr.NE.0 ) THEN
       WRITE( 6, * ) ' *** OPEN I/O error:', ierr
       RETURN
      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
       RETURN
      ENDIF  

C     Display program variables

      DO i = 0, 4095, 8
       WRITE( 6, 9001 ) i, i+7, ( gain( j ), j = i, i+7 )
      ENDDO

      DO i = 0, 4095, 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, 6
       WRITE( 6, 9005 ) i, dsssd_a(i), dsssd_b(i),
     +                     dsssd_c(i), dsssd_d(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
       WRITE( 6, 9010 ) items, FLOAT( items - items_old ) / delta_t,
     +                  events, FLOAT( events - events_old ) / delta_t
       items_old = items
       events_old = events
       t2_old = t2

      ENDIF

      IF ( MOD( items, 50000 ).EQ.0 ) THEN

       t4 = DFLOAT( ts ) * 10.0D-9 - t3
       delta_t = t4 - t4_old

       IF ( delta_t.GE.1.0D+00 ) THEN

        DO i = 1, 32
         rates( i ) = DFLOAT( scaler( i ) - scaler_old( i ) ) / delta_t
        ENDDO
C        WRITE( 6, * ) ' '
        DO i = 1, 6
C         WRITE( 6, 9020 ) i, scaler( i ), scaler_old( i ),
C     +                    delta_t, rates( i )
         sc_channel = INT ( MOD( t4, 65536.0 ) + 0.5 )
         CALL set1d( i+300, sc_channel, INT( rates( i ) + 0.5 ) )
        ENDDO
        DO i = 7, 12
C         WRITE( 6, 9030 ) i-6, scaler( i ), scaler_old( i ),
C     +                    delta_t, rates( i )
         sc_channel = INT ( MOD( t4, 65536.0 ) + 0.5 )
         CALL set1d( i+300, sc_channel, INT( rates( i ) + 0.5 ) )
        ENDDO
        DO i = 1, 32
         scaler_old( i ) = scaler( i )
        ENDDO
        t4_old = t4
       ENDIF

      ENDIF

C----67---------------------------------------------------------------72------80
C     Event data

C     For Linux/x64 & Windows/x86 least significant 16-bit word i2(1)
C     For Solaris/SPARC least significant 16-bit word i2(2)

C     GREAT 3.2.1 format
... 1616 more lines ...
Attachment 8: variables.dat
  $variables
  dsssd_a( 1 ) = 17
  dsssd_a( 2 ) = 18
  dsssd_a( 3 ) = 9
  dsssd_a( 4 ) = 10
  dsssd_a( 5 ) = 1
  dsssd_a( 6 ) = 2
  dsssd_b( 1 ) = 20
  dsssd_b( 2 ) = 19
  dsssd_b( 3 ) = 12
  dsssd_b( 4 ) = 11
  dsssd_b( 5 ) = 4
  dsssd_b( 6 ) = 3
  dsssd_c( 1 ) = 22
  dsssd_c( 2 ) = 21
  dsssd_c( 3 ) = 14
  dsssd_c( 4 ) = 13
  dsssd_c( 5 ) = 6
  dsssd_c( 6 ) = 5
  dsssd_d( 1 ) = 23
  dsssd_d( 2 ) = 24
  dsssd_d( 3 ) = 15
  dsssd_d( 4 ) = 16
  dsssd_d( 5 ) = 7
  dsssd_d( 6 ) = 8
  offset = 4096*0.0
  gain = 4096*1.0
  c = 8*1.0
  w = 32*0
  w(1) = 4000
  w(2) = -4000
  w(3) = 210
  w(4) = 210
  $[end]
Attachment 9: titles.dat
s1 nnaida1: ADC data items (20us/ch)
s2 nnaida2: ADC data items (20us/ch)
s3 nnaida3: ADC data items (20us/ch)
s4 nnaida4: ADC data items (20us/ch)
s5 nnaida5: ADC data items (20us/ch)
s6 nnaida6: ADC data items (20us/ch)
s7 nnaida7: ADC data items (20us/ch)
s8 nnaida8: ADC data items (20us/ch)
s9 nnaida9: ADC data items (20us/ch)
s10 nnaida10: ADC data items (20us/ch)
s11 nnaida11: ADC data items (20us/ch)
s12 nnaida12: ADC data items (20us/ch)
s13 nnaida13: ADC data items (20us/ch)
s14 nnaida14: ADC data items (20us/ch)
s15 nnaida15: ADC data items (20us/ch)
s16 nnaida16: ADC data items (20us/ch)
s17 nnaida17: ADC data items (20us/ch)
s18 nnaida18: ADC data items (20us/ch)
s19 nnaida19: ADC data items (20us/ch)
s20 nnaida20: ADC data items (20us/ch)
s21 nnaida21: ADC data items (20us/ch)
s22 nnaida22: ADC data items (20us/ch)
s23 nnaida23: ADC data items (20us/ch)
s24 nnaida24: ADC data items (20us/ch)
s25 nnaida25: ADC data items (20us/ch)
s26 nnaida26: ADC data items (20us/ch)
s27 nnaida27: ADC data items (20us/ch)
s28 nnaida28: ADC data items (20us/ch)
s29 nnaida29: ADC data items (20us/ch)
s30 nnaida30: ADC data items (20us/ch)
s31 nnaida31: ADC data items (20us/ch)
s32 nnaida32: ADC data items (20us/ch)

s101 nnaida1: HEC ADC data items (20us/ch)
s124 nnaida1: HEC ADC data items (20us/ch)

s201 nnaida1: Disc data items (20us/ch)
s202 nnaida2: Disc data items (20us/ch)
s203 nnaida3: Disc data items (20us/ch)
s204 nnaida4: Disc data items (20us/ch)
s205 nnaida5: Disc data items (20us/ch)
s206 nnaida6: Disc data items (20us/ch)
s207 nnaida7: Disc data items (20us/ch)
s208 nnaida8: Disc data items (20us/ch)
s209 nnaida9: Disc data items (20us/ch)
s210 nnaida10: Disc data items (20us/ch)
s211 nnaida11: Disc data items (20us/ch)
s212 nnaida12: Disc data items (20us/ch)
s213 nnaida13: Disc data items (20us/ch)
s214 nnaida14: Disc data items (20us/ch)
s215 nnaida15: Disc data items (20us/ch)
s216 nnaida16: Disc data items (20us/ch)
s217 nnaida17: Disc data items (20us/ch)
s218 nnaida18: Disc data items (20us/ch)
s219 nnaida19: Disc data items (20us/ch)
s220 nnaida20: Disc data items (20us/ch)
s221 nnaida21: Disc data items (20us/ch)
s222 nnaida22: Disc data items (20us/ch)
s223 nnaida23: Disc data items (20us/ch)
s224 nnaida24: Disc data items (20us/ch)
s225 nnaida25: Disc data items (20us/ch)
s226 nnaida26: Disc data items (20us/ch)
s227 nnaida27: Disc data items (20us/ch)
s228 nnaida28: Disc data items (20us/ch)
s229 nnaida29: Disc data items (20us/ch)
s230 nnaida30: Disc data items (20us/ch)
s231 nnaida31: Disc data items (20us/ch)
s232 nnaida32: Disc data items (20us/ch)

s250 LEC & HEC channels (DSSSD order)
s251 LEC & HEC channels (ASIC order)

s301 DSSSD #1 LEC rate (Hz) versus time 
s306 DSSSD #6 LEC rate (Hz) versus time
s307 DSSSD #1 HEC rate (Hz) versus time 
s312 DSSSD #6 HEC rate (Hz) versus time

s1000 4*(module # - 1) + asic # (0-95 ADC; 100-195 Disc) v. timestamp (20us/ch)
s1001 HEC channel (2048-3584) versus energy

s2000 m_p versus m_n DSSSD #1 (LEC)
s2005 m_p versus m_n DSSSD #6 (LEC)
s2010 x versus y DSSSD #1 (LEC)
s2015 x versus y DSSSD #6 (LEC)
s2020 E_p versus E_n DSSSD #1 (LEC)
s2025 E_p versus E_n DSSSD #6 (LEC)
s2030 LEC-LEC time DSSSD #1 (2.56us/ch)
s2035 LEC-LEC time DSSSD #6 (2.56us/ch)

s2100 m_p versus m_n DSSSD #1 (HEC)
s2105 m_p versus m_n DSSSD #6 (HEC)
s2110 x versus y DSSSD #1 (HEC)
s2115 x versus y DSSSD #6 (HEC)
s2120 E_p versus E_n DSSSD #1 (HEC)
s2125 E_p versus E_n DSSSD #6 (HEC)

s2130 HEC-HEC time DSSSD #1 (2.56us/ch)
s2135 HEC-HEC time DSSSD #6 (2.56us/ch)

s2200 per pixel HEC-LEC time DSSSD #1 (2.56us/ch)
s2205 per pixel HEC-LEC time DSSSD #6 (2.56us/ch)

s2210 per pixel HEC-LEC time DSSSD #1 (163.840us/ch)
s2215 per pixel HEC-LEC time DSSSD #6 (163.840us/ch)

s2300 HEC E versus HEC-LEC dt DSSSD#1 (163.840us/ch)
s2305 HEC E versus HEC-LEC dt DSSSD#6 (163.840us/ch)

s2400 LEC E versus HEC-LEC dt DSSSD#1 (163.840us/ch)
s2405 LEC E versus HEC-LEC dt DSSSD#6 (163.840us/ch)

s2500 HEC xy versus HEC-LEC dt DSSSD#1 (163.840us/ch)
s2505 HEC xy versus HEC-LEC dt DSSSD#6 (163.840us/ch)

s2600 LEC x_time - e_time DSSSD#1 (10ns/ch)
s2605 LEC x_time - e_time DSSSD#6 (10ns/ch)

s2700 e_diff versus HEC-LEC dt DSSSD#1 (163.840us/ch)
s2705 e_diff versus HEC-LEC dt DSSSD#6 (163.840us/ch)

s2800 m/z_lec/ctr versus HEC-LEC dt DSSSD#1 (163.840us/ch)
s2805 m/z_lec/ctr versus HEC-LEC dt DSSSD#6 (163.840us/ch)

s2900 LEC-LEC dt (20.48us/ch) versus HEC-LEC dt DSSSD#1 (163.840us/ch)
s2905 LEC-LEC dt (20.48us/ch) versus HEC-LEC dt DSSSD#6 (163.840us/ch)

s3000 FEE64 versus FEE64
Attachment 10: 107.png
107.png
Attachment 11: 115.png
115.png
  6   Mon Sep 15 13:51:28 2014 PatrickOdd spectra from running nnaida2 for 3 hours

 I'm running nnaida2 with the PB-5 and an old adaptor.(1pF)

I've seen spectra like this before but previously they have shown apparent digital relationship between the troughs. 

So maybe they show 3 peaks ? 

If so I don't quite understand the mechanism......

I just realised that the shaping time is 0.5uS so maybe that is why ??

 

Attachment 1: nnaida2_3hours
ELOG V3.1.3-7933898