/[PAMELA software]/DarthVader/ToFLevel2/src/toftrk.for
ViewVC logotype

Diff of /DarthVader/ToFLevel2/src/toftrk.for

Parent Directory Parent Directory | Revision Log Revision Log | View Patch Patch

revision 1.1 by mocchiut, Sat Jun 17 12:14:56 2006 UTC revision 1.12 by pam-de, Mon Aug 20 14:22:33 2007 UTC
# Line 1  Line 1 
1  *****************************************************************************  *****************************************************************************
2        INTEGER FUNCTION TOFTRK()        INTEGER FUNCTION TOFTRK()
3  c      
4    C****************************************************************************
5    C     31-08-06  WM
6    C     Changed to use DOTRACK2
7    C     Beta calculation: now the flightpath (instead of cos(theta)) is used
8    C     Beta calculation: all 4 TDV measurements must be < 4095 (in the old
9    C     routine it was (t1+t2)<8000
10    C
11    C     08-12-06 WM:
12    C     adc_c-bug :  The raw ADC value was multiplied with cos(theta)
13    C     and AFTER that there was an if statement "if tof32(right,i,iadc) < 4095"
14    C
15    C  jan-07 GF: ADC/TDCflags(4,12) inserted to flag artificial ADC/TDC
16    C             values
17    C  jan-07 WM: artificial ADC values created using attenuation calibration
18    C  jan-07 WM: artificial TDC values created using xy_coor calibration
19    C  jan-07 WM: modified xtofpos flag "101". xtofpos must be inside physical
20    C             dimension of the paddle +/- 10 cm
21    C  jan-07 WM: if xtofpos=101 then this paddle is not used for beta
22    C             calculation
23    C  jan-07 WM: in the xtofpos calculation a check for TDC.ne.4095 was
24    C             inserted. In the old code one would still calculate a
25    C             xtofpos-value even if the TDC information was missing
26    C  jan-07 WM: flag for PMTs #10 and #35 added, TDC=819 due to bit-shift
27    C  jan-05 WM: bug fixed: calculation of zenith angles using DOTRACK2
28    C             was incorrect
29    C  jan-07 WM: bug fixed: in some cases tdc_tw was calculated due to a
30    C             leftover "xhelp" value
31    C  apr-07 WM: attenuation fit curve is now a double exponential fit
32    C             conversion from raw ADC to pC using calibration function
33    C             new variables xtr_tof(6) and ytr_tof(6) give track position
34    C             at ToF layers
35    C  aug-07 WM: artificial ADC creation revised: Now an ADC value is created
36    C             only if there is a TDC value (before ADC was created in ANY
37    C             case)
38    C
39    C****************************************************************************
40        IMPLICIT NONE        IMPLICIT NONE
41  C      C
42        include  'input_tof.txt'        include  'input_tof.txt'
43        include  'output_tof.txt'        include  'output_tof.txt'
44        include  'tofcomm.txt'        include  'tofcomm.txt'
45  C      C
46          
47  c     =======================================  c     =======================================
48  c     variables for tracking routine  c     variables for tracking routine
49  c     =======================================  c     =======================================
# Line 20  c     define TOF Z-coordinates Line 56  c     define TOF Z-coordinates
56        DOUBLE PRECISION ZTOF(NPTOF)        DOUBLE PRECISION ZTOF(NPTOF)
57        DATA ZTOF/53.74,53.04,23.94,23.44,-23.49,-24.34/ !Sergio 9.05.2006        DATA ZTOF/53.74,53.04,23.94,23.44,-23.49,-24.34/ !Sergio 9.05.2006
58    
59        integer itof        integer itof,pmt_id
60    
61        DOUBLE PRECISION al_p(5),        DOUBLE PRECISION al_p(5),
62       &     xout(NPOINT_MAX),yout(NPOINT_MAX),zin(NPTOF)       &     xout(NPOINT_MAX),yout(NPOINT_MAX),zin(NPTOF),
63         &     THXOUT(NPOINT_MAX),THYOUT(NPOINT_MAX),TLOUT(NPOINT_MAX)
64    
65    
66        INTEGER IFAIL        INTEGER IFAIL
67        REAL dx,dy,dr,ds  c      REAL dx,dy,dr
68        REAL yhelp,xdummy,xkorr0,xhelp,xhelp1,xhelp2        REAL ds
69          REAL t1,t2,t3,t4
70          REAL yhelp,xhelp,xhelp1,xhelp2
71        REAL c1,c2,sw,sxw,w_i        REAL c1,c2,sw,sxw,w_i
72        INTEGER icount        REAL dist,dl,F
73          INTEGER icount,ievent
74          REAL xhelp_a,xhelp_t
75    
76        REAL beta_mean        REAL beta_mean
77          REAL hepratio
78    
79          INTEGER j
80    
81          real atten,pc_adc
82    
83        REAL theta12,theta13,theta23  
84          REAL theta,phi
85  C--   DATA ZTOF/53.74,53.04,23.94,23.44,-23.49,-24.34/ !Sergio 9.05.2006  C--   DATA ZTOF/53.74,53.04,23.94,23.44,-23.49,-24.34/ !Sergio 9.05.2006
86        REAL tofarm12        REAL tofarm12
87        PARAMETER (tofarm12 = 29.70) ! from 53.39 to 23.69        PARAMETER (tofarm12 = 29.70) ! from 53.39 to 23.69
88        REAL tofarm23                REAL tofarm23
89        PARAMETER (tofarm23 = 47.61)  ! from 23.69 to -23.92        PARAMETER (tofarm23 = 47.61) ! from 23.69 to -23.92
90        REAL tofarm13        REAL tofarm13
91        PARAMETER (tofarm13 = 77.31)  ! from 53.39 to -23.92        PARAMETER (tofarm13 = 77.31) ! from 53.39 to -23.92
92    
93    
         
         
         
94        INTEGER ihelp        INTEGER ihelp
95        REAL xkorr        REAL xkorr,xpos
96    
97        REAL yl,yh,xl,xh        REAL yl,yh,xl,xh
98  C      C
99        REAL hmemor(9000000)        REAL hmemor(9000000)
100        INTEGER Iquest(100)        INTEGER Iquest(100)
101  C      C
102          DATA ievent / 0 /
103    
104        COMMON / pawcd / hmemor        COMMON / pawcd / hmemor
105        save / pawcd /        save / pawcd /
106  C      C
107        Common / QUESTd / Iquest        Common / QUESTd / Iquest
108        save / questd /        save / questd /
109  C      C
110  C     Begin !  C     Begin !
111  C      C
112        TOFTRK = 0        TOFTRK = 0
113    
114  *******************************************************************  *******************************************************************
115    
116          ievent = ievent +1
117    
118    C  ratio helium to proton ca. 4
119          hepratio = 4.
120    
121        offset = 1        offset = 1
122        slope = 2        slope = 2
123        left = 1        left = 1
# Line 81  C     Line 133  C    
133           beta_a(i) = 100.           beta_a(i) = 100.
134        enddo        enddo
135    
136  C------ read  tracking routine        do i=1,4
137             do j=1,12
138                adc_c(i,j) = 1000.
139             enddo
140          enddo
141    
142          do i=1,12
143             do j=1,4
144                tofmask(j,i) = 0
145             enddo
146          enddo
147    
148          do i=1,4
149             do j=1,12
150                adcflag(i,j) = 0
151             enddo
152          enddo
153    
154          do i=1,4
155             do j=1,12
156                tdcflag(i,j) = 0
157             enddo
158          enddo
159    
160          pmt_id=0
161    
162          do j=1,6
163          THXOUT(j) = 0.
164          THYOUT(j) = 0.
165          enddo
166    
167          do j=1,6
168          xtr_tof(j) = 100.
169          ytr_tof(j) = 100.
170          enddo
171    
172    C----------------------------------------------------------------------
173    C-------------------------get ToF data --------------------------------
174    C     we cannot use the tofxx(x,x,x)  data  from tofl2com since it is
175    C     manipulated (Time-walk, artificila ADc and TDC values using ToF
176    C     standalone information
177    C----------------------------------------------------------------------
178    
179    c     put the adc and tdc values from ntuple into tofxx(i,j,k) variables
180    
181          do j=1,8
182             tof11(1,j,2) = pc_adc(adc(ch11a(j),hb11a(j)))
183             tof11(2,j,2) = pc_adc(adc(ch11b(j),hb11b(j)))
184             tof11(1,j,1) = (tdc(ch11a(j),hb11a(j)))
185             tof11(2,j,1) = (tdc(ch11b(j),hb11b(j)))
186          enddo
187    
188    
189          do j=1,6
190             tof12(1,j,2) = pc_adc(adc(ch12a(j),hb12a(j)))
191             tof12(2,j,2) = pc_adc(adc(ch12b(j),hb12b(j)))
192             tof12(1,j,1) = (tdc(ch12a(j),hb12a(j)))
193             tof12(2,j,1) = (tdc(ch12b(j),hb12b(j)))
194          enddo
195    
196          do j=1,2
197             tof21(1,j,2) = pc_adc(adc(ch21a(j),hb21a(j)))
198             tof21(2,j,2) = pc_adc(adc(ch21b(j),hb21b(j)))
199             tof21(1,j,1) = (tdc(ch21a(j),hb21a(j)))
200             tof21(2,j,1) = (tdc(ch21b(j),hb21b(j)))
201          enddo
202    
203          do j=1,2
204             tof22(1,j,2) = pc_adc(adc(ch22a(j),hb22a(j)))
205             tof22(2,j,2) = pc_adc(adc(ch22b(j),hb22b(j)))
206             tof22(1,j,1) = (tdc(ch22a(j),hb22a(j)))
207             tof22(2,j,1) = (tdc(ch22b(j),hb22b(j)))
208          enddo
209    
210          do j=1,3
211             tof31(1,j,2) = pc_adc(adc(ch31a(j),hb31a(j)))
212             tof31(2,j,2) = pc_adc(adc(ch31b(j),hb31b(j)))
213             tof31(1,j,1) = (tdc(ch31a(j),hb31a(j)))
214             tof31(2,j,1) = (tdc(ch31b(j),hb31b(j)))
215          enddo
216    
217          do j=1,3
218             tof32(1,j,2) = pc_adc(adc(ch32a(j),hb32a(j)))
219             tof32(2,j,2) = pc_adc(adc(ch32b(j),hb32b(j)))
220             tof32(1,j,1) = (tdc(ch32a(j),hb32a(j)))
221             tof32(2,j,1) = (tdc(ch32b(j),hb32b(j)))
222          enddo
223    
224    C----------------------------------------------------------------------
225    
226          DO i = 1,8
227             if (abs(tof11(1,i,itdc)).gt.10000.) tof11(1,i,itdc)= 10000.
228             if (abs(tof11(2,i,itdc)).gt.10000.) tof11(2,i,itdc)= 10000.
229             if (abs(tof11(1,i,iadc)).gt.10000.) tof11(1,i,iadc)= 10000.
230             if (abs(tof11(2,i,iadc)).gt.10000.) tof11(2,i,iadc)= 10000.
231          ENDDO
232    
233          DO i = 1,6
234             if (abs(tof12(1,i,itdc)).gt.10000.) tof12(1,i,itdc)= 10000.
235             if (abs(tof12(2,i,itdc)).gt.10000.) tof12(2,i,itdc)= 10000.
236             if (abs(tof12(1,i,iadc)).gt.10000.) tof12(1,i,iadc)= 10000.
237             if (abs(tof12(2,i,iadc)).gt.10000.) tof12(2,i,iadc)= 10000.
238          ENDDO
239    
240    
241          DO i = 1,2
242             if (abs(tof21(1,i,itdc)).gt.10000.) tof21(1,i,itdc)= 10000.
243             if (abs(tof21(2,i,itdc)).gt.10000.) tof21(2,i,itdc)= 10000.
244             if (abs(tof21(1,i,iadc)).gt.10000.) tof21(1,i,iadc)= 10000.
245             if (abs(tof21(2,i,iadc)).gt.10000.) tof21(2,i,iadc)= 10000.
246          ENDDO
247    
248          DO i = 1,2
249             if (abs(tof22(1,i,itdc)).gt.10000.) tof22(1,i,itdc)= 10000.
250             if (abs(tof22(2,i,itdc)).gt.10000.) tof22(2,i,itdc)= 10000.
251             if (abs(tof22(1,i,iadc)).gt.10000.) tof22(1,i,iadc)= 10000.
252             if (abs(tof22(2,i,iadc)).gt.10000.) tof22(2,i,iadc)= 10000.
253          ENDDO
254    
255          DO i = 1,3
256             if (abs(tof31(1,i,itdc)).gt.10000.) tof31(1,i,itdc)= 10000.
257             if (abs(tof31(2,i,itdc)).gt.10000.) tof31(2,i,itdc)= 10000.
258             if (abs(tof31(1,i,iadc)).gt.10000.) tof31(1,i,iadc)= 10000.
259             if (abs(tof31(2,i,iadc)).gt.10000.) tof31(2,i,iadc)= 10000.
260          ENDDO
261    
262          DO i = 1,3
263             if (abs(tof32(1,i,itdc)).gt.10000.) tof32(1,i,itdc)= 10000.
264             if (abs(tof32(2,i,itdc)).gt.10000.) tof32(2,i,itdc)= 10000.
265             if (abs(tof32(1,i,iadc)).gt.10000.) tof32(1,i,iadc)= 10000.
266             if (abs(tof32(2,i,iadc)).gt.10000.) tof32(2,i,iadc)= 10000.
267          ENDDO
268    
269    C----------------------------------------------------------------------
270    
271    C------read  tracking routine
272  *     igoodevent = igoodevent+1  *     igoodevent = igoodevent+1
273  *     assigned input  parameters for track routine  *     assigned input  parameters for track routine
274  *     1) Z-coordinates where the trajectory is evaluated  *     1) Z-coordinates where the trajectory is evaluated
# Line 93  C     COPY THE ALFA VECTOR FROM AL_PP TO Line 280  C     COPY THE ALFA VECTOR FROM AL_PP TO
280        do i=1,5        do i=1,5
281           AL_P(i) = al_pp(i)           AL_P(i) = al_pp(i)
282        enddo        enddo
283    
284    c      write(*,*) AL_P
285    
286        if (al_p(5).eq.0.) THEN        if (al_p(5).eq.0.) THEN
287           PRINT *,' TOF - WARNING F77: track with R = 0, discarded'           PRINT *,' TOF - WARNING F77: track with R = 0, discarded'
288           GOTO 969           GOTO 969
289        ENDIF        ENDIF
290  *     -------- *** tracking routine *** --------  *     -------- *** tracking routine *** --------
291        IFAIL = 0        IFAIL = 0
292        call DOTRACK(NPTOF,ZIN,XOUT,YOUT,AL_P,IFAIL)  C     call DOTRACK(NPTOF,ZIN,XOUT,YOUT,AL_P,IFAIL)
293          call DOTRACK2(NPTOF,ZIN,XOUT,YOUT,THXOUT,THYOUT,TLOUT,AL_P,IFAIL)
294    
295    C     write(*,*) (TLOUT(i),i=1,6)
296    
297        if(IFAIL.ne.0)then        if(IFAIL.ne.0)then
298           print *,' TOF - WARNING F77: tracking failed '           print *,' TOF - WARNING F77: tracking failed '
299           goto 969           goto 969
300        endif        endif
301  *     ------------------------------------------  *     ------------------------------------------
302    
303  969   continue   969  continue
304    
305    C---  Fill xtr_tof  and ytr_tof: positions from tracker at ToF layers
306          do j=1,6
307          xtr_tof(j) = XOUT(j)
308          ytr_tof(j) = YOUT(j)
309          enddo
310    
311    
312    C---  convert  angles to radian
313          do j=1,6
314          THXOUT(j) = 3.1415927*THXOUT(j)/180.
315          THYOUT(j) = 3.1415927*THYOUT(j)/180.
316          enddo
317    
318          do j=1,6
319    c      write (*,*) j,THXOUT(j),THYOUT(j)
320          enddo
321    
322    
323    C----------------------------------------------------------------------
324    C------------------  set ADC & TDC flag = 0    ------------------------
325    C----------------------------------------------------------------------
326    
327          do j=1,8
328          if (adc(ch11a(j),hb11a(j)).LT.4096)adcflagtof(ch11a(j),hb11a(j))=0
329          if (adc(ch11b(j),hb11b(j)).LT.4096)adcflagtof(ch11b(j),hb11b(j))=0
330          if (tdc(ch11a(j),hb11a(j)).LT.4096)tdcflagtof(ch11a(j),hb11a(j))=0
331          if (tdc(ch11b(j),hb11b(j)).LT.4096)tdcflagtof(ch11b(j),hb11b(j))=0
332          enddo
333          do j=1,6
334          if (adc(ch12a(j),hb12a(j)).LT.4096)adcflagtof(ch12a(j),hb12a(j))=0
335          if (adc(ch12b(j),hb12b(j)).LT.4096)adcflagtof(ch12b(j),hb12b(j))=0
336          if (tdc(ch12a(j),hb12a(j)).LT.4096)tdcflagtof(ch12a(j),hb12a(j))=0
337          if (tdc(ch12b(j),hb12b(j)).LT.4096)tdcflagtof(ch12b(j),hb12b(j))=0
338          enddo
339          do j=1,2
340          if (adc(ch21a(j),hb21a(j)).LT.4096)adcflagtof(ch21a(j),hb21a(j))=0
341          if (adc(ch21b(j),hb21b(j)).LT.4096)adcflagtof(ch21b(j),hb21b(j))=0
342          if (tdc(ch21a(j),hb21a(j)).LT.4096)tdcflagtof(ch21a(j),hb21a(j))=0
343          if (tdc(ch21b(j),hb21b(j)).LT.4096)tdcflagtof(ch21b(j),hb21b(j))=0
344          enddo
345          do j=1,2
346          if (adc(ch22a(j),hb22a(j)).LT.4096)adcflagtof(ch22a(j),hb22a(j))=0
347          if (adc(ch22b(j),hb22b(j)).LT.4096)adcflagtof(ch22b(j),hb22b(j))=0
348          if (tdc(ch22a(j),hb22a(j)).LT.4096)tdcflagtof(ch22a(j),hb22a(j))=0
349          if (tdc(ch22b(j),hb22b(j)).LT.4096)tdcflagtof(ch22b(j),hb22b(j))=0
350          enddo
351          do j=1,3
352          if (adc(ch31a(j),hb31a(j)).LT.4096)adcflagtof(ch31a(j),hb31a(j))=0
353          if (adc(ch31b(j),hb31b(j)).LT.4096)adcflagtof(ch31b(j),hb31b(j))=0
354          if (tdc(ch31a(j),hb31a(j)).LT.4096)tdcflagtof(ch31a(j),hb31a(j))=0
355          if (tdc(ch31b(j),hb31b(j)).LT.4096)tdcflagtof(ch31b(j),hb31b(j))=0
356          enddo
357          do j=1,3
358          if (adc(ch32a(j),hb32a(j)).LT.4096)adcflagtof(ch32a(j),hb32a(j))=0
359          if (adc(ch32b(j),hb32b(j)).LT.4096)adcflagtof(ch32b(j),hb32b(j))=0
360          if (tdc(ch32a(j),hb32a(j)).LT.4096)tdcflagtof(ch32a(j),hb32a(j))=0
361          if (tdc(ch32b(j),hb32b(j)).LT.4096)tdcflagtof(ch32b(j),hb32b(j))=0
362          enddo
363    
364    
365    C----------------------------------------------------------------
366    C---------- Check PMTs 10 and 35 for strange TDC values----------
367    C----------------------------------------------------------------
368    
369    C---- S116A TDC=819
370           if (tof11(1,6,1).EQ.819) then
371                 tof11(1,6,1) = 4095
372                 tdcflagtof(ch11a(6),hb11a(6))=2
373           endif
374    
375    C---- S222B TDC=819
376           if (tof22(2,2,1).EQ.819) then
377                 tof22(2,2,1) = 4095
378                 tdcflagtof(ch22b(2),hb22b(2))=2
379           endif
380    
381  C-------------------------------------------------------------  C-------------------------------------------------------------
382  C-------  check which paddle penetrated the track  -----------  C-------check which paddle penetrated the track  -----------
383  C-------------------------------------------------------------  C-------------------------------------------------------------
384  c middle y (or x) position of the upper and middle ToF-Paddle  c     middle y (or x) position of the upper and middle ToF-Paddle
385  c       DATA tof11_x/ -17.85,-12.75,-7.65,-2.55,2.55,7.65,12.75,17.85/  c     DATA tof11_x/ -17.85,-12.75,-7.65,-2.55,2.55,7.65,12.75,17.85/
386  c       DATA tof12_y/ -13.75,-8.25,-2.75,2.75,8.25,13.75/  c     DATA tof12_y/ -13.75,-8.25,-2.75,2.75,8.25,13.75/
387  c       DATA tof21_y/ -3.75,3.75/  c     DATA tof21_y/  3.75,-3.75/  ! paddles in different order
388  c       DATA tof22_x/ -4.5,4.5/  c     DATA tof22_x/ -4.5,4.5/
389  c       DATA tof31_x/ -6.0,0.,6.0/  c     DATA tof31_x/ -6.0,0.,6.0/
390  c       DATA tof32_y/ -5.0,0.0,5.0/  c     DATA tof32_y/ -5.0,0.0,5.0/
391  c  c
392  c  S11 8 paddles  33.0 x 5.1 cm  c     S11 8 paddles  33.0 x 5.1 cm
393  c  S12 6 paddles  40.8 x 5.5 cm  c     S12 6 paddles  40.8 x 5.5 cm
394  c  S21 2 paddles  18.0 x 7.5 cm  c     S21 2 paddles  18.0 x 7.5 cm
395  c  S22 2 paddles  15.0 x 9.0 cm  c     S22 2 paddles  15.0 x 9.0 cm
396  c  S31 3 paddles  15.0 x 6.0 cm  c     S31 3 paddles  15.0 x 6.0 cm
397  c  S32 3 paddles  18.0 x 5.0 cm  c     S32 3 paddles  18.0 x 5.0 cm
398    
399    c     write(*,*) xout(1),xout(2),xout(3),xout(4),xout(5),xout(6)
400    c     write(*,*) yout(1),yout(2),yout(3),yout(4),yout(5),yout(6)
401    
402  C--------------     S11 --------------------------------------  C--------------S11 --------------------------------------
403    
404        tof11_i = none_find        tof11_i = none_find
405    
# Line 135  C--------------     S11 ---------------- Line 407  C--------------     S11 ----------------
407        yh =  33.0/2.        yh =  33.0/2.
408    
409        if ((yout(1).gt.yl).and.(yout(1).lt.yh)) then        if ((yout(1).gt.yl).and.(yout(1).lt.yh)) then
410        do i=1,8           do i=1,8
411        xl = tof11_x(i) - 5.1/2.              xl = tof11_x(i) - 5.1/2.
412        xh = tof11_x(i) + 5.1/2.              xh = tof11_x(i) + 5.1/2.
413        if ((xout(1).gt.xl).and.(xout(1).le.xh)) then              if ((xout(1).gt.xl).and.(xout(1).le.xh)) then
414        tof11_i=i                 tof11_i=i
415        endif              endif
416        enddo           enddo
417        endif        endif
418    
419  C--------------     S12 --------------------------------------  C--------------S12 --------------------------------------
420    
421        tof12_i = none_find        tof12_i = none_find
422    
423        xl = -40.8/2.        xl = -40.8/2.
424        xh =  40.8/2.        xh =  40.8/2.
425    
426        if ((xout(2).gt.xl).and.(xout(2).lt.xh)) then        if ((xout(2).gt.xl).and.(xout(2).lt.xh)) then
427        do i=1,6           do i=1,6
428        yl = tof12_y(i) - 5.5/2.              yl = tof12_y(i) - 5.5/2.
429        yh = tof12_y(i) + 5.5/2.              yh = tof12_y(i) + 5.5/2.
430        if ((yout(2).gt.yl).and.(yout(2).le.yh)) then              if ((yout(2).gt.yl).and.(yout(2).le.yh)) then
431        tof12_i=i                 tof12_i=i
432        endif              endif
433        enddo           enddo
434        endif        endif
435    
436  C--------------     S21 --------------------------------------  C--------------S21 --------------------------------------
437    
438        tof21_i = none_find        tof21_i = none_find
439    
440        xl = -18./2.        xl = -18./2.
441        xh =  18./2.        xh =  18./2.
442    
443        if ((xout(3).gt.xl).and.(xout(3).lt.xh)) then        if ((xout(3).gt.xl).and.(xout(3).lt.xh)) then
444        do i=1,2           do i=1,2
445        yl = tof21_y(i) - 7.5/2.              yl = tof21_y(i) - 7.5/2.
446        yh = tof21_y(i) + 7.5/2.              yh = tof21_y(i) + 7.5/2.
447        if ((yout(3).gt.yl).and.(yout(3).le.yh)) then              if ((yout(3).gt.yl).and.(yout(3).le.yh)) then
448        tof21_i=i                 tof21_i=i
449        endif              endif
450        enddo           enddo
451        endif        endif
452    
453  C--------------     S22 --------------------------------------  C--------------S22 --------------------------------------
454    
455        tof22_i = none_find        tof22_i = none_find
456    
457        yl = -15./2.        yl = -15./2.
458        yh =  15./2.        yh =  15./2.
459    
460        if ((yout(4).gt.yl).and.(yout(4).lt.yh)) then        if ((yout(4).gt.yl).and.(yout(4).lt.yh)) then
461        do i=1,2           do i=1,2
462        xl = tof22_x(i) - 9.0/2.              xl = tof22_x(i) - 9.0/2.
463        xh = tof22_x(i) + 9.0/2.              xh = tof22_x(i) + 9.0/2.
464        if ((xout(4).gt.xl).and.(xout(4).le.xh)) then              if ((xout(4).gt.xl).and.(xout(4).le.xh)) then
465        tof22_i=i                 tof22_i=i
466        endif              endif
467        enddo           enddo
468        endif        endif
469    
470  C--------------     S31 --------------------------------------  C--------------S31 --------------------------------------
471    
472        tof31_i = none_find        tof31_i = none_find
473    
474        yl = -15.0/2.        yl = -15.0/2.
475        yh =  15.0/2.        yh =  15.0/2.
476    
477        if ((yout(5).gt.yl).and.(yout(5).lt.yh)) then        if ((yout(5).gt.yl).and.(yout(5).lt.yh)) then
478        do i=1,3           do i=1,3
479        xl = tof31_x(i) - 6.0/2.              xl = tof31_x(i) - 6.0/2.
480        xh = tof31_x(i) + 6.0/2.              xh = tof31_x(i) + 6.0/2.
481        if ((xout(5).gt.xl).and.(xout(5).le.xh)) then              if ((xout(5).gt.xl).and.(xout(5).le.xh)) then
482        tof31_i=i                 tof31_i=i
483        endif              endif
484        enddo           enddo
485        endif        endif
486    
487  C--------------     S32 --------------------------------------  C--------------S32 --------------------------------------
488    
489        tof32_i = none_find        tof32_i = none_find
490    
491        xl = -18.0/2.        xl = -18.0/2.
492        xh =  18.0/2.        xh =  18.0/2.
493    
494        if ((xout(6).gt.xl).and.(xout(6).lt.xh)) then        if ((xout(6).gt.xl).and.(xout(6).lt.xh)) then
495        do i=1,3           do i=1,3
496        yl = tof32_y(i) - 5.0/2.              yl = tof32_y(i) - 5.0/2.
497        yh = tof32_y(i) + 5.0/2.              yh = tof32_y(i) + 5.0/2.
498        if ((yout(6).gt.yl).and.(yout(6).le.yh)) then              if ((yout(6).gt.yl).and.(yout(6).le.yh)) then
499        tof32_i=i                 tof32_i=i
500                endif
501             enddo
502        endif        endif
503    
504    
505    C     write(*,*) tof11_i,tof12_i,tof21_i,tof22_i,tof31_i,tof32_i
506    
507    C-----------------------------------------------------------------------
508    C--------------------Insert Artifical TDC Value  ---------------------
509    C     For each Paddle perform check:
510    C     if left paddle=4095  and right paddle OK => create TDC value left
511    C     if right paddle=4095  and left paddle OK => create TDC value right
512    C-----------------------------------------------------------------------
513    
514    C-----------------------S11 -----------------------------------------
515    
516          IF (tof11_i.GT.none_find) THEN
517             xpos = yout(1)
518             i = tof11_i
519             if ((tof11(1,tof11_i,itdc).EQ.4095).AND.
520         &        (tof11(2,tof11_i,itdc).LT.4095)) THEN
521    
522    c       write(*,*)'11lb4 ',i,tof11(1,tof11_i,itdc),tof11(2,tof11_i,itdc)
523    
524                tof11(1,tof11_i,itdc) = tof11(2,tof11_i,itdc)
525         &           + 2*(y_coor_lin11(tof11_i,offset)
526         &           + xpos*y_coor_lin11(tof11_i,slope))
527    
528    c       write(*,*)'11laf ',i,tof11(1,tof11_i,itdc),tof11(2,tof11_i,itdc)
529    
530                tdcflag(ch11a(i),hb11a(i)) = 1
531    
532            ENDIF
533             if ((tof11(2,tof11_i,itdc).EQ.4095).AND.
534         &        (tof11(1,tof11_i,itdc).LT.4095)) THEN
535    
536    c       write(*,*)'11rb4 ',i,tof11(1,tof11_i,itdc),tof11(2,tof11_i,itdc)
537    
538                tof11(2,tof11_i,itdc) = tof11(1,tof11_i,itdc)
539         &           - 2*(y_coor_lin11(tof11_i,offset)
540         &           + xpos*y_coor_lin11(tof11_i,slope))
541    c       write(*,*)'11raf ',i,tof11(1,tof11_i,itdc),tof11(2,tof11_i,itdc)
542    
543                tdcflag(ch11b(i),hb11b(i)) = 1
544             ENDIF
545          ENDIF
546    
547    C-----------------------S12 -----------------------------------------
548    
549          IF (tof12_i.GT.none_find) THEN
550             xpos = xout(2)
551             i = tof12_i
552             if ((tof12(1,tof12_i,itdc).EQ.4095).AND.
553         &        (tof12(2,tof12_i,itdc).LT.4095)) THEN
554                tof12(1,tof12_i,itdc) = tof12(2,tof12_i,itdc)
555         &           + 2*(x_coor_lin12(tof12_i,offset)
556         &           + xpos*x_coor_lin12(tof12_i,slope))
557                tdcflag(ch12a(i),hb12a(i)) = 1
558             ENDIF
559             if ((tof12(2,tof12_i,itdc).EQ.4095).AND.
560         &        (tof12(1,tof12_i,itdc).LT.4095)) THEN
561                tof12(2,tof12_i,itdc) = tof12(1,tof12_i,itdc)
562         &           - 2*(x_coor_lin12(tof12_i,offset)
563         &           + xpos*x_coor_lin12(tof12_i,slope))
564                tdcflag(ch12b(i),hb12b(i)) = 1
565             ENDIF
566          ENDIF
567    
568    C-----------------------S21 -----------------------------------------
569    
570          IF (tof21_i.GT.none_find) THEN
571             xpos = xout(3)
572             i = tof21_i
573             if ((tof21(1,tof21_i,itdc).EQ.4095).AND.
574         &        (tof21(2,tof21_i,itdc).LT.4095)) THEN
575                tof21(1,tof21_i,itdc) = tof21(2,tof21_i,itdc)
576         &           + 2*(x_coor_lin21(tof21_i,offset)
577         &           + xpos*x_coor_lin21(tof21_i,slope))
578                tdcflag(ch21a(i),hb21a(i)) = 1
579             ENDIF
580             if ((tof21(2,tof21_i,itdc).EQ.4095).AND.
581         &        (tof21(1,tof21_i,itdc).LT.4095)) THEN
582                tof21(2,tof21_i,itdc) = tof21(1,tof21_i,itdc)
583         &           - 2*(x_coor_lin21(tof21_i,offset)
584         &           + xpos*x_coor_lin21(tof21_i,slope))
585                tdcflag(ch21b(i),hb21b(i)) = 1
586             ENDIF
587          ENDIF
588    
589    C-----------------------S22 -----------------------------------------
590    
591          IF (tof22_i.GT.none_find) THEN
592             xpos = yout(4)
593             i = tof22_i
594             if ((tof22(1,tof22_i,itdc).EQ.4095).AND.
595         &        (tof22(2,tof22_i,itdc).LT.4095)) THEN
596                tof22(1,tof22_i,itdc) = tof22(2,tof22_i,itdc)
597         &           + 2*(y_coor_lin22(tof22_i,offset)
598         &           + xpos*y_coor_lin22(tof22_i,slope))
599                tdcflag(ch22a(i),hb22a(i)) = 1
600             ENDIF
601             if ((tof22(2,tof22_i,itdc).EQ.4095).AND.
602         &        (tof22(1,tof22_i,itdc).LT.4095)) THEN
603                tof22(2,tof22_i,itdc) = tof22(1,tof22_i,itdc)
604         &           - 2*(y_coor_lin22(tof22_i,offset)
605         &           + xpos*y_coor_lin22(tof22_i,slope))
606                tdcflag(ch22b(i),hb22b(i)) = 1
607             ENDIF
608          ENDIF
609    
610    C-----------------------S31 -----------------------------------------
611    
612          IF (tof31_i.GT.none_find) THEN
613             xpos = yout(5)
614             i = tof31_i
615             if ((tof31(1,tof31_i,itdc).EQ.4095).AND.
616         &        (tof31(2,tof31_i,itdc).LT.4095)) THEN
617                tof31(1,tof31_i,itdc) = tof31(2,tof31_i,itdc)
618         &           + 2*(y_coor_lin31(tof31_i,offset)
619         &           + xpos*y_coor_lin31(tof31_i,slope))
620                tdcflag(ch31a(i),hb31a(i)) = 1
621             ENDIF
622             if ((tof31(2,tof31_i,itdc).EQ.4095).AND.
623         &        (tof31(1,tof31_i,itdc).LT.4095)) THEN
624                tof31(2,tof31_i,itdc) = tof31(1,tof31_i,itdc)
625         &           - 2*(y_coor_lin31(tof31_i,offset)
626         &           + xpos*y_coor_lin31(tof31_i,slope))
627                tdcflag(ch31b(i),hb31b(i)) = 1
628             ENDIF
629          ENDIF
630    
631    C-----------------------S32 -----------------------------------------
632    
633          IF (tof32_i.GT.none_find) THEN
634             xpos = xout(6)
635             i = tof32_i
636             if ((tof32(1,tof32_i,itdc).EQ.4095).AND.
637         &        (tof32(2,tof32_i,itdc).LT.4095)) THEN
638                tof32(1,tof32_i,itdc) = tof32(2,tof32_i,itdc)
639         &           + 2*(x_coor_lin32(tof32_i,offset)
640         &           + xpos*x_coor_lin32(tof32_i,slope))
641                tdcflag(ch32a(i),hb32a(i)) = 1
642             ENDIF
643             if ((tof32(2,tof32_i,itdc).EQ.4095).AND.
644         &        (tof32(1,tof32_i,itdc).LT.4095)) THEN
645                tof32(2,tof32_i,itdc) = tof32(1,tof32_i,itdc)
646         &           - 2*(x_coor_lin32(tof32_i,offset)
647         &           + xpos*x_coor_lin32(tof32_i,slope))
648                tdcflag(ch32b(i),hb32b(i)) = 1
649             ENDIF
650          ENDIF
651    
652    C--------------------------------------------------------------------
653    C---- if paddle hit: if we have TDC value but no ADC, create ADC value
654    C--------------------------------------------------------------------
655    c     middle y (or x) position of the upper and middle ToF-Paddle
656    c     DATA tof11_x/ -17.85,-12.75,-7.65,-2.55,2.55,7.65,12.75,17.85/
657    c     DATA tof12_y/ -13.75,-8.25,-2.75,2.75,8.25,13.75/
658    c     DATA tof21_y/  3.75,-3.75/  ! paddles in different order
659    c     DATA tof22_x/ -4.5,4.5/
660    c     DATA tof31_x/ -6.0,0.,6.0/
661    c     DATA tof32_y/ -5.0,0.0,5.0/
662    
663    C----------------------------S1 -------------------------------------
664    
665          yhelp=yout(1)
666          IF (tof11_i.GT.none_find.AND.abs(yout(1)).lt.100) THEN
667             i = tof11_i
668             if ((tof11(left,tof11_i,itdc).LT.4095).AND.
669         &       (adc(ch11a(i),hb11a(i)).eq.4095)) then
670                phi   = atan(tan(THYOUT(1))/tan(THXOUT(1)))
671                theta = atan(tan(THXOUT(1))/cos(phi))
672    c            xkorr=adcx11(left,i,1)*exp(-yhelp/adcx11(left,i,2))
673                xkorr = atten(left,11,i,yhelp)
674                xkorr=xkorr/hepratio
675                tof11(left,i,iadc)=xkorr/cos(theta)
676                adcflag(ch11a(i),hb11a(i)) = 1
677             endif
678             if ((tof11(right,tof11_i,itdc).LT.4095).AND.
679         &       (adc(ch11b(i),hb11b(i)).eq.4095)) then
680                phi   = atan(tan(THYOUT(1))/tan(THXOUT(1)))
681                theta = atan(tan(THXOUT(1))/cos(phi))
682    c            xkorr=adcx11(right,i,1)*exp(yhelp/adcx11(right,i,2))
683                xkorr = atten(right,11,i,yhelp)
684                xkorr=xkorr/hepratio
685                tof11(right,i,iadc)=xkorr/cos(theta)
686                adcflag(ch11b(i),hb11b(i)) = 1
687             endif
688          ENDIF
689    
690          xhelp=xout(2)
691          IF (tof12_i.GT.none_find.AND.abs(xout(2)).lt.100) THEN
692             i = tof12_i
693             if ((tof12(left,tof12_i,itdc).LT.4095).AND.
694         &       (adc(ch12a(i),hb12a(i)).eq.4095)) then
695                phi   = atan(tan(THYOUT(2))/tan(THXOUT(2)))
696                theta = atan(tan(THXOUT(2))/cos(phi))
697    c            xkorr=adcx12(left,i,1)*exp(-xhelp/adcx12(left,i,2))
698                xkorr = atten(left,12,i,xhelp)
699                xkorr=xkorr/hepratio
700                tof12(left,i,iadc) = xkorr/cos(theta)
701                adcflag(ch12a(i),hb12a(i)) = 1
702             endif
703             if ((tof12(right,tof12_i,itdc).LT.4095).AND.
704         &       (adc(ch12b(i),hb12b(i)).eq.4095)) then
705                phi   = atan(tan(THYOUT(2))/tan(THXOUT(2)))
706                theta = atan(tan(THXOUT(2))/cos(phi))
707    c            xkorr=adcx12(right,i,1)*exp(xhelp/adcx12(right,i,2))
708                xkorr = atten(right,12,i,xhelp)
709                xkorr=xkorr/hepratio
710                tof12(right,i,iadc) = xkorr/cos(theta)
711                adcflag(ch12b(i),hb12b(i)) = 1
712             endif
713          ENDIF
714    
715    C-----------------------------S2 --------------------------------
716    
717          xhelp=xout(3)
718          IF (tof21_i.GT.none_find.AND.abs(xout(3)).lt.100) THEN
719             i = tof21_i
720             if ((tof21(left,tof21_i,itdc).LT.4095).AND.
721         &       (adc(ch21a(i),hb21a(i)).eq.4095)) then
722                phi   = atan(tan(THYOUT(3))/tan(THXOUT(3)))
723                theta = atan(tan(THXOUT(3))/cos(phi))
724    c            xkorr=adcx21(left,i,1)*exp(-xhelp/adcx21(left,i,2))
725                xkorr = atten(left,21,i,xhelp)
726                xkorr=xkorr/hepratio
727                tof21(left,i,iadc) = xkorr/cos(theta)
728                adcflag(ch21a(i),hb21a(i)) = 1
729             endif
730             if ((tof21(right,tof21_i,itdc).LT.4095).AND.
731         &       (adc(ch21b(i),hb21b(i)).eq.4095)) then
732                phi   = atan(tan(THYOUT(3))/tan(THXOUT(3)))
733                theta = atan(tan(THXOUT(3))/cos(phi))
734    c            xkorr=adcx21(right,i,1)*exp(xhelp/adcx21(right,i,2))
735                xkorr = atten(right,21,i,xhelp)
736                xkorr=xkorr/hepratio
737                tof21(right,i,iadc) = xkorr/cos(theta)
738                adcflag(ch21b(i),hb21b(i)) = 1
739             endif
740          ENDIF
741    
742    
743          yhelp=yout(4)
744          IF (tof22_i.GT.none_find.AND.abs(yout(4)).lt.100) THEN
745             i = tof22_i
746             if ((tof22(left,tof22_i,itdc).LT.4095).AND.
747         &       (adc(ch22a(i),hb22a(i)).eq.4095)) then
748                phi   = atan(tan(THYOUT(4))/tan(THXOUT(4)))
749                theta = atan(tan(THXOUT(4))/cos(phi))
750    c            xkorr=adcx22(left,i,1)*exp(-yhelp/adcx22(left,i,2))
751                xkorr = atten(left,22,i,yhelp)
752                xkorr=xkorr/hepratio
753                tof22(left,i,iadc) = xkorr/cos(theta)
754                adcflag(ch22a(i),hb22a(i)) = 1
755             endif
756             if ((tof22(right,tof22_i,itdc).LT.4095).AND.
757         &       (adc(ch22b(i),hb22b(i)).eq.4095)) then
758                phi   = atan(tan(THYOUT(4))/tan(THXOUT(4)))
759                theta = atan(tan(THXOUT(4))/cos(phi))
760    c            xkorr=adcx22(right,i,1)*exp(yhelp/adcx22(right,i,2))
761                xkorr = atten(right,22,i,yhelp)
762                xkorr=xkorr/hepratio
763                tof22(right,i,iadc) = xkorr/cos(theta)
764                adcflag(ch22b(i),hb22b(i)) = 1
765             endif
766          ENDIF
767    
768    C-----------------------------S3 --------------------------------
769    
770          yhelp=yout(5)
771          IF (tof31_i.GT.none_find.AND.abs(yout(5)).lt.100) THEN
772             i = tof31_i
773             if ((tof31(left,tof31_i,itdc).LT.4095).AND.
774         &       (adc(ch31a(i),hb31a(i)).eq.4095)) then
775                phi   = atan(tan(THYOUT(5))/tan(THXOUT(5)))
776                theta = atan(tan(THXOUT(5))/cos(phi))
777    c            xkorr=adcx31(left,i,1)*exp(-yhelp/adcx31(left,i,2))
778                xkorr = atten(left,31,i,yhelp)
779                xkorr=xkorr/hepratio
780                tof31(left,i,iadc) = xkorr/cos(theta)
781                adcflag(ch31a(i),hb31a(i)) = 1
782             endif
783             if ((tof31(right,tof31_i,itdc).LT.4095).AND.
784         &       (adc(ch31b(i),hb31b(i)).eq.4095)) then
785                phi   = atan(tan(THYOUT(5))/tan(THXOUT(5)))
786                theta = atan(tan(THXOUT(5))/cos(phi))
787    c            xkorr=adcx31(right,i,1)*exp(yhelp/adcx31(right,i,2))
788                xkorr = atten(right,31,i,yhelp)
789                xkorr=xkorr/hepratio
790                tof31(right,i,iadc) = xkorr/cos(theta)
791                adcflag(ch31b(i),hb31b(i)) = 1
792             endif
793          ENDIF
794    
795    
796          xhelp=xout(6)
797          IF (tof32_i.GT.none_find.AND.abs(xout(6)).lt.100) THEN
798             i = tof32_i
799             if ((tof32(left,tof32_i,itdc).LT.4095).AND.
800         &       (adc(ch32a(i),hb32a(i)).eq.4095)) then
801                phi   = atan(tan(THYOUT(6))/tan(THXOUT(6)))
802                theta = atan(tan(THXOUT(6))/cos(phi))
803    c            xkorr=adcx32(left,i,1)*exp(-xhelp/adcx32(left,i,2))
804                xkorr = atten(left,32,i,xhelp)
805                xkorr=xkorr/hepratio
806                tof32(left,i,iadc) = xkorr/cos(theta)
807                adcflag(ch32a(i),hb32a(i)) = 1
808             endif
809             if ((tof32(right,tof32_i,itdc).LT.4095).AND.
810         &       (adc(ch32b(i),hb32b(i)).eq.4095)) then
811                phi   = atan(tan(THYOUT(6))/tan(THXOUT(6)))
812                theta = atan(tan(THXOUT(6))/cos(phi))
813    c            xkorr=adcx32(right,i,1)*exp(xhelp/adcx32(right,i,2))
814                xkorr = atten(right,32,i,xhelp)
815                xkorr=xkorr/hepratio
816                tof32(right,i,iadc) = xkorr/cos(theta)
817                adcflag(ch32b(i),hb32b(i)) = 1
818             endif
819          ENDIF
820    
821    
822    C------------------------------------------------------------------
823    C---  calculate track position in paddle using timing difference
824    C------------------------------------------------------------------
825    
826          do i=1,3
827             xtofpos(i)=100.
828             ytofpos(i)=100.
829        enddo        enddo
830    C-----------------------------S1 --------------------------------
831    
832          IF (tof11_i.GT.none_find) THEN
833             IF ((tof11(1,tof11_i,itdc).NE.4095).AND.
834         &       (tof11(2,tof11_i,itdc).NE.4095)) THEN
835             ytofpos(1)  = ((tof11(1,tof11_i,itdc)-tof11(2,tof11_i,itdc))/2.
836         +        -y_coor_lin11(tof11_i,offset))/y_coor_lin11(tof11_i,slope)
837            if (abs(ytofpos(1)).gt.26.)  ytofpos(1)=101.
838          endif
839        endif        endif
840    
841  C----------------------------------------------------------------------        IF (tof12_i.GT.none_find) THEN
842  C--------------------Corrections on ADC-data -------------------------           IF ((tof12(1,tof12_i,itdc).NE.4095).AND.
843  C---------------------zenith angle theta  ---------------------------       &       (tof12(2,tof12_i,itdc).NE.4095)) THEN
844  C----------------------------------------------------------------------           xtofpos(1)  = ((tof12(1,tof12_i,itdc)-tof12(2,tof12_i,itdc))/2.
845         +        -x_coor_lin12(tof12_i,offset))/x_coor_lin12(tof12_i,slope)
846            if (abs(xtofpos(1)).gt.31.)  xtofpos(1)=101.
847          endif
848          endif
849    
850    C-----------------------------S2 --------------------------------
851    
852        dx=0.        IF (tof21_i.GT.none_find) THEN
853        dy=0.           IF ((tof21(1,tof21_i,itdc).NE.4095).AND.
854        dr=0.       &       (tof21(2,tof21_i,itdc).NE.4095)) THEN
855        theta13 = 0.           xtofpos(2) = ((tof21(1,tof21_i,itdc)-tof21(2,tof21_i,itdc))/2.
856         +        -x_coor_lin21(tof21_i,offset))/x_coor_lin21(tof21_i,slope)
857        if (xout(1).lt.100.) then          if (abs(xtofpos(2)).gt.19.)  xtofpos(2)=101.
858           dx = xout(1)-xout(6)        endif
          dy = yout(1)-yout(6)  
          dr = sqrt(dx*dx+dy*dy)  
          theta13 = atan(dr/tofarm13)  
859        endif        endif
860    
861          IF (tof22_i.GT.none_find) THEN
862             IF ((tof22(1,tof22_i,itdc).NE.4095).AND.
863         &       (tof22(2,tof22_i,itdc).NE.4095)) THEN
864             ytofpos(2) = ((tof22(1,tof22_i,itdc)-tof22(2,tof22_i,itdc))/2.
865         +        -y_coor_lin22(tof22_i,offset))/y_coor_lin22(tof22_i,slope)
866            if (abs(ytofpos(2)).gt.18.)  ytofpos(2)=101.
867          endif
868          endif
869    
870  C----------------------------------------------------------------------  C-----------------------------S3 --------------------------------
871  C------------------angle and ADC(x) correction  
872  C----------------------------------------------------------------------        IF (tof31_i.GT.none_find) THEN
873  C-----------------------------S1 --------------------------------           IF ((tof31(1,tof31_i,itdc).NE.4095).AND.
874         &       (tof31(2,tof31_i,itdc).NE.4095)) THEN
875             ytofpos(3)  = ((tof31(1,tof31_i,itdc)-tof31(2,tof31_i,itdc))/2.
876         +        -y_coor_lin31(tof31_i,offset))/y_coor_lin31(tof31_i,slope)
877            if (abs(ytofpos(3)).gt.18.)  ytofpos(3)=101.
878          endif
879          endif
880    
881          IF (tof32_i.GT.none_find) THEN
882             IF ((tof32(1,tof32_i,itdc).NE.4095).AND.
883         &       (tof32(2,tof32_i,itdc).NE.4095)) THEN
884             xtofpos(3)  = ((tof32(1,tof32_i,itdc)-tof32(2,tof32_i,itdc))/2.
885         +        -x_coor_lin32(tof32_i,offset))/x_coor_lin32(tof32_i,slope)
886            if (abs(xtofpos(3)).gt.19.)  xtofpos(3)=101.
887          endif
888          endif
889    
890    c      do i=1,3
891    c         if (abs(xtofpos(i)).gt.100.) then
892    c            xtofpos(i)=101.
893    c         endif
894    c         if (abs(ytofpos(i)).gt.100.) then
895    c            ytofpos(i)=101.
896    c         endif
897    c      enddo
898    
899    
900    
901    
902    C--------------------------------------------------------------------
903    C--------------------Time walk correction  -------------------------
904    C--------------------------------------------------------------------
905    
906    
907          DO i=1,8
908             xhelp= 0.
909             xhelp_a = tof11(left,i,iadc)
910             xhelp_t = tof11(left,i,itdc)
911              if (xhelp_a .eq.0) write (*,*) 'trk 11 ',i,xhelp_a
912    
913             if(xhelp_a<3786) xhelp = tw11(left,i)/sqrt(xhelp_a)
914             tof11(left,i,itdc) = xhelp_t  + xhelp
915             tdc_c(ch11a(i),hb11a(i))=tof11(left,i,itdc)
916             xhelp_a = tof11(right,i,iadc)
917             xhelp_t = tof11(right,i,itdc)
918             if(xhelp_a<3786) xhelp = tw11(right,i)/sqrt(xhelp_a)
919             tof11(right,i,itdc) = xhelp_t  + xhelp
920             tdc_c(ch11b(i),hb11b(i))=tof11(right,i,itdc)
921          ENDDO
922    
923          DO i=1,6
924             xhelp= 0.
925             xhelp_a = tof12(left,i,iadc)
926             xhelp_t = tof12(left,i,itdc)
927              if (xhelp_a .eq.0) write (*,*) 'trk 12 ',i,xhelp_a
928             if(xhelp_a<3786) xhelp = tw12(left,i)/sqrt(xhelp_a)
929             tof12(left,i,itdc) = xhelp_t  + xhelp
930             tdc_c(ch12a(i),hb12a(i))=tof12(left,i,itdc)
931             xhelp_a = tof12(right,i,iadc)
932             xhelp_t = tof12(right,i,itdc)
933             if(xhelp_a<3786) xhelp = tw12(right,i)/sqrt(xhelp_a)
934             tof12(right,i,itdc) = xhelp_t  + xhelp
935             tdc_c(ch12b(i),hb12b(i))=tof12(right,i,itdc)
936          ENDDO
937    C----
938          DO i=1,2
939             xhelp= 0.
940             xhelp_a = tof21(left,i,iadc)
941             xhelp_t = tof21(left,i,itdc)
942              if (xhelp_a .eq.0) write (*,*) ' trk 21 ',i,xhelp_a
943    
944             if(xhelp_a<3786) xhelp = tw21(left,i)/sqrt(xhelp_a)
945             tof21(left,i,itdc) = xhelp_t  + xhelp
946             tdc_c(ch21a(i),hb21a(i))=tof21(left,i,itdc)
947             xhelp_a = tof21(right,i,iadc)
948             xhelp_t = tof21(right,i,itdc)
949             if(xhelp_a<3786) xhelp = tw21(right,i)/sqrt(xhelp_a)
950             tof21(right,i,itdc) = xhelp_t  + xhelp
951             tdc_c(ch21b(i),hb21b(i))=tof21(right,i,itdc)
952          ENDDO
953    
954          DO i=1,2
955             xhelp= 0.
956             xhelp_a = tof22(left,i,iadc)
957             xhelp_t = tof22(left,i,itdc)
958              if (xhelp_a .eq.0) write (*,*) ' trk 22 ',i,xhelp_a
959             if(xhelp_a<3786) xhelp = tw22(left,i)/sqrt(xhelp_a)
960             tof22(left,i,itdc) = xhelp_t  + xhelp
961             tdc_c(ch22a(i),hb22a(i))=tof22(left,i,itdc)
962             xhelp_a = tof22(right,i,iadc)
963             xhelp_t = tof22(right,i,itdc)
964             if(xhelp_a<3786) xhelp = tw22(right,i)/sqrt(xhelp_a)
965             tof22(right,i,itdc) = xhelp_t  + xhelp
966             tdc_c(ch22b(i),hb22b(i))=tof22(right,i,itdc)
967          ENDDO
968    C----
969    
970          DO i=1,3
971             xhelp= 0.
972             xhelp_a = tof31(left,i,iadc)
973             xhelp_t = tof31(left,i,itdc)
974              if (xhelp_a .eq.0) write (*,*) ' trk 31 ',i,xhelp_a
975    
976             if(xhelp_a<3786) xhelp = tw31(left,i)/sqrt(xhelp_a)
977             tof31(left,i,itdc) = xhelp_t  + xhelp
978             tdc_c(ch31a(i),hb31a(i))=tof31(left,i,itdc)
979             xhelp_a = tof31(right,i,iadc)
980             xhelp_t = tof31(right,i,itdc)
981             if(xhelp_a<3786) xhelp = tw31(right,i)/sqrt(xhelp_a)
982             tof31(right,i,itdc) = xhelp_t  + xhelp
983             tdc_c(ch31b(i),hb31b(i))=tof31(right,i,itdc)
984          ENDDO
985    
986          DO i=1,3
987             xhelp= 0.
988             xhelp_a = tof32(left,i,iadc)
989             xhelp_t = tof32(left,i,itdc)
990              if (xhelp_a .eq.0) write (*,*) ' trk 32 ',i,xhelp_a
991    
992             if(xhelp_a<3786) xhelp = tw32(left,i)/sqrt(xhelp_a)
993             tof32(left,i,itdc) = xhelp_t  + xhelp
994             tdc_c(ch32a(i),hb32a(i))=tof32(left,i,itdc)
995             xhelp_a = tof32(right,i,iadc)
996             xhelp_t = tof32(right,i,itdc)
997             if(xhelp_a<3786) xhelp = tw32(right,i)/sqrt(xhelp_a)
998             tof32(right,i,itdc) = xhelp_t  + xhelp
999             tdc_c(ch32b(i),hb32b(i))=tof32(right,i,itdc)
1000          ENDDO
1001    
1002    
1003    C---------------------------------------------------------------------
1004    C--------------------Corrections on ADC-data -------------------------
1005    C-----------------angle and ADC(x) correction -----------------------
1006    
1007    C-----------------------------S1 -------------------------------------
1008    
1009        yhelp=yout(1)        yhelp=yout(1)
1010    
1011          phi   = atan(tan(THYOUT(1))/tan(THXOUT(1)))
1012          theta = atan(tan(THXOUT(1))/cos(phi))
1013    
1014        IF (tof11_i.GT.none_find.AND.yhelp.lt.100) THEN        IF (tof11_i.GT.none_find.AND.yhelp.lt.100) THEN
1015    
1016           i = tof11_i           i = tof11_i
1017           xdummy=tof11(left,i,iadc)  
1018           tof11(left,i,iadc) = tof11(left,i,iadc)*cos(theta13)           if (tof11(left,i,iadc).lt.3786) then
1019           if (tof11(left,i,iadc).lt.1000) then  c         if (adc(ch11a(i),hb11a(i)).lt.4095) then
1020              xkorr=adcx11(left,i,1)*exp(-yhelp/adcx11(left,i,2))              tof11(left,i,iadc) = tof11(left,i,iadc)*cos(theta)
1021              xkorr0=adcx11(left,i,1)  c            xkorr=adcx11(left,i,1)*exp(-yhelp/adcx11(left,i,2))
1022                xkorr = atten(left,11,i,yhelp)
1023                xkorr=xkorr/hepratio
1024              adc_c(ch11a(i),hb11a(i))=tof11(left,i,iadc)/xkorr              adc_c(ch11a(i),hb11a(i))=tof11(left,i,iadc)/xkorr
1025           endif           endif
1026    
1027           tof11(right,i,iadc) = tof11(right,i,iadc)*cos(theta13)  
1028           if (tof11(right,i,iadc).lt.1000) then           if (tof11(right,i,iadc).lt.3786) then
1029              xkorr=adcx11(right,i,1)*exp(yhelp/adcx11(right,i,2))  c         if (adc(ch11b(i),hb11b(i)).lt.4095) then
1030              xkorr0=adcx11(right,i,1)              tof11(right,i,iadc) = tof11(right,i,iadc)*cos(theta)
1031    c            xkorr=adcx11(right,i,1)*exp(yhelp/adcx11(right,i,2))
1032                xkorr = atten(right,11,i,yhelp)
1033                xkorr=xkorr/hepratio
1034              adc_c(ch11b(i),hb11b(i))=tof11(right,i,iadc)/xkorr              adc_c(ch11b(i),hb11b(i))=tof11(right,i,iadc)/xkorr
1035           endif           endif
1036        ENDIF        ENDIF
1037    
1038    
1039        xhelp=xout(2)        xhelp=xout(2)
1040          phi   = atan(tan(THYOUT(2))/tan(THXOUT(2)))
1041          theta = atan(tan(THXOUT(2))/cos(phi))
1042    c      write(*,*) 'theta12 ',theta
1043        IF (tof12_i.GT.none_find.AND.xhelp.lt.100) THEN        IF (tof12_i.GT.none_find.AND.xhelp.lt.100) THEN
1044    
1045           i = tof12_i           i = tof12_i
1046           tof12(left,i,iadc) = tof12(left,i,iadc)*cos(theta13)           if (tof12(left,i,iadc).lt.3786) then
1047           if (tof12(left,i,iadc).lt.1000) then  c         if (adc(ch12a(i),hb12a(i)).lt.4095) then
1048              xkorr=adcx12(left,i,1)*exp(-xhelp/adcx12(left,i,2))              tof12(left,i,iadc) = tof12(left,i,iadc)*cos(theta)
1049              xkorr0=adcx12(left,i,1)  c            xkorr=adcx12(left,i,1)*exp(-xhelp/adcx12(left,i,2))
1050                xkorr = atten(left,12,i,xhelp)
1051                xkorr=xkorr/hepratio
1052              adc_c(ch12a(i),hb12a(i))=tof12(left,i,iadc)/xkorr              adc_c(ch12a(i),hb12a(i))=tof12(left,i,iadc)/xkorr
1053           endif           endif
1054    
1055           tof12(right,i,iadc) = tof12(right,i,iadc)*cos(theta13)           if (tof12(right,i,iadc).lt.3786) then
1056           if (tof12(right,i,iadc).lt.1000) then  c         if (adc(ch12b(i),hb12b(i)).lt.4095) then
1057              xkorr=adcx12(right,i,1)*exp(xhelp/adcx12(right,i,2))              tof12(right,i,iadc) = tof12(right,i,iadc)*cos(theta)
1058              xkorr0=adcx12(right,i,1)  c            xkorr=adcx12(right,i,1)*exp(xhelp/adcx12(right,i,2))
1059                xkorr = atten(right,12,i,xhelp)
1060                xkorr=xkorr/hepratio
1061              adc_c(ch12b(i),hb12b(i))=tof12(right,i,iadc)/xkorr              adc_c(ch12b(i),hb12b(i))=tof12(right,i,iadc)/xkorr
1062           endif           endif
1063        ENDIF        ENDIF
# Line 297  C-----------------------------S1 ------- Line 1065  C-----------------------------S1 -------
1065  C-----------------------------S2 --------------------------------  C-----------------------------S2 --------------------------------
1066    
1067        xhelp=xout(3)        xhelp=xout(3)
1068          phi   = atan(tan(THYOUT(3))/tan(THXOUT(3)))
1069          theta = atan(tan(THXOUT(3))/cos(phi))
1070    c      write(*,*) 'theta21 ',theta
1071        IF (tof21_i.GT.none_find.AND.xhelp.lt.100) THEN        IF (tof21_i.GT.none_find.AND.xhelp.lt.100) THEN
1072    
1073           i = tof21_i           i = tof21_i
1074           tof21(left,i,iadc) = tof21(left,i,iadc)*cos(theta13)           if (tof21(left,i,iadc).lt.3786) then
1075           if (tof21(left,i,iadc).lt.1000) then  c         if (adc(ch21a(i),hb21a(i)).lt.4095) then
1076              xkorr=adcx21(left,i,1)*exp(-xhelp/adcx21(left,i,2))              tof21(left,i,iadc) = tof21(left,i,iadc)*cos(theta)
1077              xkorr0=adcx21(left,i,1)  c            xkorr=adcx21(left,i,1)*exp(-xhelp/adcx21(left,i,2))
1078                xkorr = atten(left,21,i,xhelp)
1079                xkorr=xkorr/hepratio
1080              adc_c(ch21a(i),hb21a(i))=tof21(left,i,iadc)/xkorr              adc_c(ch21a(i),hb21a(i))=tof21(left,i,iadc)/xkorr
1081           endif           endif
1082    
1083           tof21(right,i,iadc) = tof21(right,i,iadc)*cos(theta13)           if (tof21(right,i,iadc).lt.3786) then
1084           if (tof21(right,i,iadc).lt.1000) then  c         if (adc(ch21b(i),hb21b(i)).lt.4095) then
1085              xkorr=adcx21(right,i,1)*exp(xhelp/adcx21(right,i,2))              tof21(right,i,iadc) = tof21(right,i,iadc)*cos(theta)
1086              xkorr0=adcx21(right,i,1)  c            xkorr=adcx21(right,i,1)*exp(xhelp/adcx21(right,i,2))
1087                xkorr = atten(right,21,i,xhelp)
1088                xkorr=xkorr/hepratio
1089              adc_c(ch21b(i),hb21b(i))=tof21(right,i,iadc)/xkorr              adc_c(ch21b(i),hb21b(i))=tof21(right,i,iadc)/xkorr
1090           endif           endif
1091        ENDIF        ENDIF
1092    
1093        yhelp=yout(4)        yhelp=yout(4)
1094          phi   = atan(tan(THYOUT(4))/tan(THXOUT(4)))
1095          theta = atan(tan(THXOUT(4))/cos(phi))
1096    c      write(*,*) 'theta22 ',theta
1097    
1098        IF (tof22_i.GT.none_find.AND.yhelp.lt.100) THEN        IF (tof22_i.GT.none_find.AND.yhelp.lt.100) THEN
1099    
1100           i = tof22_i           i = tof22_i
1101           tof22(left,i,iadc) = tof22(left,i,iadc)*cos(theta13)           if (tof22(left,i,iadc).lt.3786) then
1102           if (tof22(left,i,iadc).lt.1000) then  c         if (adc(ch22a(i),hb22a(i)).lt.4095) then
1103              xkorr=adcx22(left,i,1)*exp(-yhelp/adcx22(left,i,2))              tof22(left,i,iadc) = tof22(left,i,iadc)*cos(theta)
1104              xkorr0=adcx22(left,i,1)  c            xkorr=adcx22(left,i,1)*exp(-yhelp/adcx22(left,i,2))
1105                xkorr = atten(left,22,i,yhelp)
1106                xkorr=xkorr/hepratio
1107              adc_c(ch22a(i),hb22a(i))=tof22(left,i,iadc)/xkorr              adc_c(ch22a(i),hb22a(i))=tof22(left,i,iadc)/xkorr
1108           endif           endif
1109    
1110           tof22(right,i,iadc) = tof22(right,i,iadc)*cos(theta13)           if (tof22(right,i,iadc).lt.3786) then
1111           if (tof22(right,i,iadc).lt.1000) then  c         if (adc(ch22b(i),hb22b(i)).lt.4095) then
1112              xkorr=adcx22(right,i,1)*exp(yhelp/adcx22(right,i,2))              tof22(right,i,iadc) = tof22(right,i,iadc)*cos(theta)
1113              xkorr0=adcx22(right,i,1)  c            xkorr=adcx22(right,i,1)*exp(yhelp/adcx22(right,i,2))
1114                xkorr = atten(right,22,i,yhelp)
1115                xkorr=xkorr/hepratio
1116              adc_c(ch22b(i),hb22b(i))=tof22(right,i,iadc)/xkorr              adc_c(ch22b(i),hb22b(i))=tof22(right,i,iadc)/xkorr
1117           endif           endif
1118        ENDIF        ENDIF
# Line 337  C-----------------------------S2 ------- Line 1120  C-----------------------------S2 -------
1120  C-----------------------------S3 --------------------------------  C-----------------------------S3 --------------------------------
1121    
1122        yhelp=yout(5)        yhelp=yout(5)
1123          phi   = atan(tan(THYOUT(5))/tan(THXOUT(5)))
1124          theta = atan(tan(THXOUT(5))/cos(phi))
1125    c      write(*,*) 'theta31 ',theta
1126    
1127        IF (tof31_i.GT.none_find.AND.yhelp.lt.100) THEN        IF (tof31_i.GT.none_find.AND.yhelp.lt.100) THEN
1128    
1129           i = tof31_i           i = tof31_i
1130           tof31(left,i,iadc) = tof31(left,i,iadc)*cos(theta13)           if (tof31(left,i,iadc).lt.3786) then
1131           if (tof31(left,i,iadc).lt.1000) then  c         if (adc(ch31a(i),hb31a(i)).lt.4095) then
1132              xkorr=adcx31(left,i,1)*exp(-yhelp/adcx31(left,i,2))              tof31(left,i,iadc) = tof31(left,i,iadc)*cos(theta)
1133              xkorr0=adcx31(left,i,1)  c            xkorr=adcx31(left,i,1)*exp(-yhelp/adcx31(left,i,2))
1134                xkorr = atten(left,31,i,yhelp)
1135                xkorr=xkorr/hepratio
1136              adc_c(ch31a(i),hb31a(i))=tof31(left,i,iadc)/xkorr              adc_c(ch31a(i),hb31a(i))=tof31(left,i,iadc)/xkorr
1137           endif           endif
1138    
1139           tof31(right,i,iadc) = tof31(right,i,iadc)*cos(theta13)           if (tof31(right,i,iadc).lt.3786) then
1140           if (tof31(right,i,iadc).lt.1000) then  c         if (adc(ch31b(i),hb31b(i)).lt.4095) then
1141              xkorr=adcx31(right,i,1)*exp(yhelp/adcx31(right,i,2))              tof31(right,i,iadc) = tof31(right,i,iadc)*cos(theta)
1142              xkorr0=adcx31(right,i,1)  c            xkorr=adcx31(right,i,1)*exp(yhelp/adcx31(right,i,2))
1143                xkorr = atten(right,31,i,yhelp)
1144                xkorr=xkorr/hepratio
1145              adc_c(ch31b(i),hb31b(i))=tof31(right,i,iadc)/xkorr              adc_c(ch31b(i),hb31b(i))=tof31(right,i,iadc)/xkorr
1146           endif           endif
1147        ENDIF        ENDIF
1148    
1149        xhelp=xout(6)        xhelp=xout(6)
1150          phi   = atan(tan(THYOUT(6))/tan(THXOUT(6)))
1151          theta = atan(tan(THXOUT(6))/cos(phi))
1152    c      write(*,*) 'theta32 ',theta
1153    
1154        IF (tof32_i.GT.none_find.AND.xhelp.lt.100) THEN        IF (tof32_i.GT.none_find.AND.xhelp.lt.100) THEN
1155    
1156           i = tof32_i           i = tof32_i
1157           tof32(left,i,iadc) = tof32(left,i,iadc)*cos(theta13)           if (tof32(left,i,iadc).lt.3786) then
1158           if (tof32(left,i,iadc).lt.1000) then  c         if (adc(ch32a(i),hb32a(i)).lt.4095) then
1159              xkorr=adcx32(left,i,1)*exp(-xhelp/adcx32(left,i,2))              tof32(left,i,iadc) = tof32(left,i,iadc)*cos(theta)
1160              xkorr0=adcx32(left,i,1)  c            xkorr=adcx32(left,i,1)*exp(-xhelp/adcx32(left,i,2))
1161                xkorr = atten(left,32,i,xhelp)
1162                xkorr=xkorr/hepratio
1163              adc_c(ch32a(i),hb32a(i))=tof32(left,i,iadc)/xkorr              adc_c(ch32a(i),hb32a(i))=tof32(left,i,iadc)/xkorr
1164           endif           endif
1165    
1166           tof32(right,i,iadc) = tof32(right,i,iadc)*cos(theta13)           if (tof32(right,i,iadc).lt.3786) then
1167           if (tof32(right,i,iadc).lt.1000) then  c         if (adc(ch32b(i),hb32b(i)).lt.4095) then
1168              xkorr=adcx32(right,i,1)*exp(xhelp/adcx32(right,i,2))              tof32(right,i,iadc) = tof32(right,i,iadc)*cos(theta)
1169              xkorr0=adcx32(right,i,1)  c            xkorr=adcx32(right,i,1)*exp(xhelp/adcx32(right,i,2))
1170                xkorr = atten(right,32,i,xhelp)
1171                xkorr=xkorr/hepratio
1172              adc_c(ch32b(i),hb32b(i))=tof32(right,i,iadc)/xkorr              adc_c(ch32b(i),hb32b(i))=tof32(right,i,iadc)/xkorr
1173           endif           endif
1174        ENDIF        ENDIF
# Line 384  C     DS = c1 + c2/beta*cos(theta) Line 1183  C     DS = c1 + c2/beta*cos(theta)
1183  C     c2 = 2d/c   gives c2 = 2d/(c*TDCresolution)  TDC=50ps/channel  C     c2 = 2d/c   gives c2 = 2d/(c*TDCresolution)  TDC=50ps/channel
1184  C     =>  c2 =  ca.60 for  0.45 m    c2 = ca.109  for 0.81 m  C     =>  c2 =  ca.60 for  0.45 m    c2 = ca.109  for 0.81 m
1185  C     since TDC resolution varies slightly c2 has to be calibrated  C     since TDC resolution varies slightly c2 has to be calibrated
1186    C     instead of cos(theta) use factor F:
1187    C     F =  pathlength/d
1188    C     => beta = c2*F/(DS-c1))
1189    
1190          dist = ZTOF(1) - ZTOF(5)
1191          dl = 0.
1192          DO I=1,5
1193             dl = dl + TLOUT(i)
1194          ENDDO
1195          F = dl/dist
1196    
1197  C     S11 - S31  C     S11 - S31
1198        IF (tof11_i.GT.none_find.AND.tof31_i.GT.none_find) THEN  C      IF (tof11_i.GT.none_find.AND.tof31_i.GT.none_find) THEN
1199           xhelp1 = tof11(1,tof11_i,itdc)+tof11(2,tof11_i,itdc)        IF ((tof11_i.GT.none_find).AND.(tof31_i.GT.none_find).AND.
1200           xhelp2 = tof31(1,tof31_i,itdc)+tof31(2,tof31_i,itdc)       &    (ytofpos(1).NE.101.).AND.(ytofpos(3).NE.101.)) THEN
1201           ds = xhelp1-xhelp2           t1 = tof11(1,tof11_i,itdc)
1202           ihelp=(tof11_i-1)*3+tof31_i           t2 = tof11(2,tof11_i,itdc)
1203           c1 = k_S11S31(1,ihelp)           t3 = tof31(1,tof31_i,itdc)
1204           c2 = k_S11S31(2,ihelp)                   t4 = tof31(2,tof31_i,itdc)
1205           if ((xhelp1.lt.8000.).and.(xhelp2.lt.8000))           IF ((t1.lt.4095).and.(t2.lt.4095).and.
1206       &   beta_a(1) = c2/(cos(theta13)*(ds-c1))       &        (t3.lt.4095).and.(t4.lt.4095)) THEN
1207                xhelp1 = tof11(1,tof11_i,itdc)+tof11(2,tof11_i,itdc)
1208                xhelp2 = tof31(1,tof31_i,itdc)+tof31(2,tof31_i,itdc)
1209                ds = xhelp1-xhelp2
1210                ihelp=(tof11_i-1)*3+tof31_i
1211                c1 = k_S11S31(1,ihelp)
1212                c2 = k_S11S31(2,ihelp)
1213                beta_a(1) = c2*F/(ds-c1)
1214           write(*,*) 'S11-S31 ',c1,c2,F
1215           write(*,*) 'S11-S31 ',xhelp1,xhelp2, beta_a(1)
1216    C-------ToF Mask - S11 - S31
1217    
1218                tofmask(ch11a(tof11_i),hb11a(tof11_i)) =
1219         $           tofmask(ch11a(tof11_i),hb11a(tof11_i)) + 1
1220                tofmask(ch11b(tof11_i),hb11b(tof11_i)) =
1221         $           tofmask(ch11b(tof11_i),hb11b(tof11_i)) + 1
1222    
1223                tofmask(ch31a(tof31_i),hb31a(tof31_i)) =
1224         $           tofmask(ch31a(tof31_i),hb31a(tof31_i)) + 1
1225                tofmask(ch31b(tof31_i),hb31b(tof31_i)) =
1226         $           tofmask(ch31b(tof31_i),hb31b(tof31_i)) + 1
1227    
1228             ENDIF
1229        ENDIF        ENDIF
1230          
1231          dist = ZTOF(1) - ZTOF(6)
1232          dl = 0.
1233          DO I=1,6
1234             dl = dl + TLOUT(i)
1235          ENDDO
1236          F = dl/dist
1237    
1238  C     S11 - S32  C     S11 - S32
1239        IF (tof11_i.GT.none_find.AND.tof32_i.GT.none_find) THEN  C      IF (tof11_i.GT.none_find.AND.tof32_i.GT.none_find) THEN
1240           xhelp1 = tof11(1,tof11_i,itdc)+tof11(2,tof11_i,itdc)         IF ((tof11_i.GT.none_find).AND.(tof32_i.GT.none_find).AND.
1241           xhelp2 = tof32(1,tof32_i,itdc)+tof32(2,tof32_i,itdc)       &    (ytofpos(1).NE.101.).AND.(xtofpos(3).NE.101.)) THEN
1242           ds = xhelp1-xhelp2           t1 = tof11(1,tof11_i,itdc)
1243           ihelp=(tof11_i-1)*3+tof32_i           t2 = tof11(2,tof11_i,itdc)
1244           c1 = k_S11S32(1,ihelp)           t3 = tof32(1,tof32_i,itdc)
1245           c2 = k_S11S32(2,ihelp)                   t4 = tof32(2,tof32_i,itdc)
1246           if ((xhelp1.lt.8000.).and.(xhelp2.lt.8000))           IF ((t1.lt.4095).and.(t2.lt.4095).and.
1247       &   beta_a(2) = c2/(cos(theta13)*(ds-c1))       &        (t3.lt.4095).and.(t4.lt.4095)) THEN
1248                xhelp1 = tof11(1,tof11_i,itdc)+tof11(2,tof11_i,itdc)
1249                xhelp2 = tof32(1,tof32_i,itdc)+tof32(2,tof32_i,itdc)
1250                ds = xhelp1-xhelp2
1251                ihelp=(tof11_i-1)*3+tof32_i
1252                c1 = k_S11S32(1,ihelp)
1253                c2 = k_S11S32(2,ihelp)
1254                beta_a(2) = c2*F/(ds-c1)
1255    C     write(*,*) 'S11-S32 ',xhelp1,xhelp2, beta_a(2)
1256    
1257    C-------ToF Mask - S11 - S32
1258    
1259                tofmask(ch11a(tof11_i),hb11a(tof11_i)) =
1260         $           tofmask(ch11a(tof11_i),hb11a(tof11_i)) + 1
1261                tofmask(ch11b(tof11_i),hb11b(tof11_i)) =
1262         $           tofmask(ch11b(tof11_i),hb11b(tof11_i)) + 1
1263    
1264                tofmask(ch32a(tof32_i),hb32a(tof32_i)) =
1265         $           tofmask(ch32a(tof32_i),hb32a(tof32_i)) + 1
1266                tofmask(ch32b(tof32_i),hb32b(tof32_i)) =
1267         $           tofmask(ch32b(tof32_i),hb32b(tof32_i)) + 1
1268    
1269    C-------
1270    
1271             ENDIF
1272        ENDIF        ENDIF
1273    
1274  C     S12 - S31  C     S12 - S31
1275        IF (tof12_i.GT.none_find.AND.tof31_i.GT.none_find) THEN        dist = ZTOF(2) - ZTOF(5)
1276           xhelp1 = tof12(1,tof12_i,itdc)+tof12(2,tof12_i,itdc)        dl = 0.
1277           xhelp2 = tof31(1,tof31_i,itdc)+tof31(2,tof31_i,itdc)        DO I=2,5
1278           ds = xhelp1-xhelp2           dl = dl + TLOUT(i)
1279           ihelp=(tof12_i-1)*3+tof31_i        ENDDO
1280           c1 = k_S12S31(1,ihelp)        F = dl/dist
1281           c2 = k_S12S31(2,ihelp)  
1282           if ((xhelp1.lt.8000.).and.(xhelp2.lt.8000))  C      IF (tof12_i.GT.none_find.AND.tof31_i.GT.none_find) THEN
1283       &   beta_a(3) = c2/(cos(theta13)*(ds-c1))         IF ((tof12_i.GT.none_find).AND.(tof31_i.GT.none_find).AND.
1284         &    (xtofpos(1).NE.101.).AND.(ytofpos(3).NE.101.)) THEN
1285             t1 = tof12(1,tof12_i,itdc)
1286             t2 = tof12(2,tof12_i,itdc)
1287             t3 = tof31(1,tof31_i,itdc)
1288             t4 = tof31(2,tof31_i,itdc)
1289             IF ((t1.lt.4095).and.(t2.lt.4095).and.
1290         &        (t3.lt.4095).and.(t4.lt.4095)) THEN
1291                xhelp1 = tof12(1,tof12_i,itdc)+tof12(2,tof12_i,itdc)
1292                xhelp2 = tof31(1,tof31_i,itdc)+tof31(2,tof31_i,itdc)
1293                ds = xhelp1-xhelp2
1294                ihelp=(tof12_i-1)*3+tof31_i
1295                c1 = k_S12S31(1,ihelp)
1296                c2 = k_S12S31(2,ihelp)
1297                beta_a(3) = c2*F/(ds-c1)
1298    C     write(*,*) 'S12-S31 ',xhelp1,xhelp2, beta_a(3)
1299    
1300    C-------ToF Mask - S12 - S31
1301    
1302                tofmask(ch12a(tof12_i),hb12a(tof12_i)) =
1303         $           tofmask(ch12a(tof12_i),hb12a(tof12_i)) + 1
1304                tofmask(ch12b(tof12_i),hb12b(tof12_i)) =
1305         $           tofmask(ch12b(tof12_i),hb12b(tof12_i)) + 1
1306    
1307                tofmask(ch31a(tof31_i),hb31a(tof31_i)) =
1308         $           tofmask(ch31a(tof31_i),hb31a(tof31_i)) + 1
1309                tofmask(ch31b(tof31_i),hb31b(tof31_i)) =
1310         $           tofmask(ch31b(tof31_i),hb31b(tof31_i)) + 1
1311    
1312    C-------
1313    
1314             ENDIF
1315        ENDIF        ENDIF
1316    
1317  C     S12 - S32  C     S12 - S32
1318        IF (tof12_i.GT.none_find.AND.tof32_i.GT.none_find) THEN  
1319           xhelp1 = tof12(1,tof12_i,itdc)+tof12(2,tof12_i,itdc)        dist = ZTOF(2) - ZTOF(6)
1320           xhelp2 = tof32(1,tof32_i,itdc)+tof32(2,tof32_i,itdc)        dl = 0.
1321           ds = xhelp1-xhelp2        DO I=2,6
1322           ihelp=(tof12_i-1)*3+tof32_i           dl = dl + TLOUT(i)
1323           c1 = k_S12S32(1,ihelp)        ENDDO
1324           c2 = k_S12S32(2,ihelp)        F = dl/dist
1325           if ((xhelp1.lt.8000.).and.(xhelp2.lt.8000))  
1326       &   beta_a(4) = c2/(cos(theta13)*(ds-c1))  C      IF (tof12_i.GT.none_find.AND.tof32_i.GT.none_find) THEN
1327           IF ((tof12_i.GT.none_find).AND.(tof32_i.GT.none_find).AND.
1328         &    (xtofpos(1).NE.101.).AND.(xtofpos(3).NE.101.)) THEN
1329             t1 = tof12(1,tof12_i,itdc)
1330             t2 = tof12(2,tof12_i,itdc)
1331             t3 = tof32(1,tof32_i,itdc)
1332             t4 = tof32(2,tof32_i,itdc)
1333             IF ((t1.lt.4095).and.(t2.lt.4095).and.
1334         &        (t3.lt.4095).and.(t4.lt.4095)) THEN
1335                xhelp1 = tof12(1,tof12_i,itdc)+tof12(2,tof12_i,itdc)
1336                xhelp2 = tof32(1,tof32_i,itdc)+tof32(2,tof32_i,itdc)
1337                ds = xhelp1-xhelp2
1338                ihelp=(tof12_i-1)*3+tof32_i
1339                c1 = k_S12S32(1,ihelp)
1340                c2 = k_S12S32(2,ihelp)
1341                beta_a(4) = c2*F/(ds-c1)
1342    C     write(*,*) 'S12-S32 ',xhelp1,xhelp2, beta_a(4)
1343    
1344    C-------ToF Mask - S12 - S32
1345    
1346                tofmask(ch12a(tof12_i),hb12a(tof12_i)) =
1347         $           tofmask(ch12a(tof12_i),hb12a(tof12_i)) + 1
1348                tofmask(ch12b(tof12_i),hb12b(tof12_i)) =
1349         $           tofmask(ch12b(tof12_i),hb12b(tof12_i)) + 1
1350    
1351                tofmask(ch32a(tof32_i),hb32a(tof32_i)) =
1352         $           tofmask(ch32a(tof32_i),hb32a(tof32_i)) + 1
1353                tofmask(ch32b(tof32_i),hb32b(tof32_i)) =
1354         $           tofmask(ch32b(tof32_i),hb32b(tof32_i)) + 1
1355    
1356    C-------
1357    
1358             ENDIF
1359        ENDIF        ENDIF
1360    
1361  C     S21 - S31  C     S21 - S31
1362        IF (tof21_i.GT.none_find.AND.tof31_i.GT.none_find) THEN  
1363           xhelp1 = tof21(1,tof21_i,itdc)+tof21(2,tof21_i,itdc)        dist = ZTOF(3) - ZTOF(5)
1364           xhelp2 = tof31(1,tof31_i,itdc)+tof31(2,tof31_i,itdc)        dl = 0.
1365           ds = xhelp1-xhelp2        DO I=3,5
1366           ihelp=(tof21_i-1)*3+tof31_i           dl = dl + TLOUT(i)
1367           c1 = k_S21S31(1,ihelp)        ENDDO
1368           c2 = k_S21S31(2,ihelp)        F = dl/dist
1369           if ((xhelp1.lt.8000.).and.(xhelp2.lt.8000))  
1370       &   beta_a(5) = c2/(cos(theta23)*(ds-c1))  C      IF (tof21_i.GT.none_find.AND.tof31_i.GT.none_find) THEN
1371           IF ((tof21_i.GT.none_find).AND.(tof31_i.GT.none_find).AND.
1372         &    (xtofpos(2).NE.101.).AND.(ytofpos(3).NE.101.)) THEN
1373             t1 = tof21(1,tof21_i,itdc)
1374             t2 = tof21(2,tof21_i,itdc)
1375             t3 = tof31(1,tof31_i,itdc)
1376             t4 = tof31(2,tof31_i,itdc)
1377             IF ((t1.lt.4095).and.(t2.lt.4095).and.
1378         &        (t3.lt.4095).and.(t4.lt.4095)) THEN
1379                xhelp1 = tof21(1,tof21_i,itdc)+tof21(2,tof21_i,itdc)
1380                xhelp2 = tof31(1,tof31_i,itdc)+tof31(2,tof31_i,itdc)
1381                ds = xhelp1-xhelp2
1382                ihelp=(tof21_i-1)*3+tof31_i
1383                c1 = k_S21S31(1,ihelp)
1384                c2 = k_S21S31(2,ihelp)
1385                beta_a(5) = c2*F/(ds-c1)
1386    
1387    C-------ToF Mask - S21 - S31
1388    
1389                tofmask(ch21a(tof21_i),hb21a(tof21_i)) =
1390         $           tofmask(ch21a(tof21_i),hb21a(tof21_i)) + 1
1391                tofmask(ch21b(tof21_i),hb21b(tof21_i)) =
1392         $           tofmask(ch21b(tof21_i),hb21b(tof21_i)) + 1
1393    
1394                tofmask(ch31a(tof31_i),hb31a(tof31_i)) =
1395         $           tofmask(ch31a(tof31_i),hb31a(tof31_i)) + 1
1396                tofmask(ch31b(tof31_i),hb31b(tof31_i)) =
1397         $           tofmask(ch31b(tof31_i),hb31b(tof31_i)) + 1
1398    
1399    C-------
1400    
1401             ENDIF
1402        ENDIF        ENDIF
1403    
1404  C     S21 - S32  C     S21 - S32
1405        IF (tof21_i.GT.none_find.AND.tof32_i.GT.none_find) THEN  
1406           xhelp1 = tof21(1,tof21_i,itdc)+tof21(2,tof21_i,itdc)        dist = ZTOF(3) - ZTOF(6)
1407           xhelp2 = tof32(1,tof32_i,itdc)+tof32(2,tof32_i,itdc)        dl = 0.
1408           ds = xhelp1-xhelp2        DO I=3,6
1409           ihelp=(tof21_i-1)*3+tof32_i           dl = dl + TLOUT(i)
1410           c1 = k_S21S32(1,ihelp)        ENDDO
1411           c2 = k_S21S32(2,ihelp)        F = dl/dist
1412           if ((xhelp1.lt.8000.).and.(xhelp2.lt.8000))  
1413       &   beta_a(6) = c2/(cos(theta23)*(ds-c1))  C      IF (tof21_i.GT.none_find.AND.tof32_i.GT.none_find) THEN
1414           IF ((tof21_i.GT.none_find).AND.(tof32_i.GT.none_find).AND.
1415         &    (xtofpos(2).NE.101.).AND.(xtofpos(3).NE.101.)) THEN
1416             t1 = tof21(1,tof21_i,itdc)
1417             t2 = tof21(2,tof21_i,itdc)
1418             t3 = tof32(1,tof32_i,itdc)
1419             t4 = tof32(2,tof32_i,itdc)
1420             IF ((t1.lt.4095).and.(t2.lt.4095).and.
1421         &        (t3.lt.4095).and.(t4.lt.4095)) THEN
1422                xhelp1 = tof21(1,tof21_i,itdc)+tof21(2,tof21_i,itdc)
1423                xhelp2 = tof32(1,tof32_i,itdc)+tof32(2,tof32_i,itdc)
1424                ds = xhelp1-xhelp2
1425                ihelp=(tof21_i-1)*3+tof32_i
1426                c1 = k_S21S32(1,ihelp)
1427                c2 = k_S21S32(2,ihelp)
1428                beta_a(6) = c2*F/(ds-c1)
1429    
1430    C-------ToF Mask - S21 - S32
1431    
1432                tofmask(ch21a(tof21_i),hb21a(tof21_i)) =
1433         $           tofmask(ch21a(tof21_i),hb21a(tof21_i)) + 1
1434                tofmask(ch21b(tof21_i),hb21b(tof21_i)) =
1435         $           tofmask(ch21b(tof21_i),hb21b(tof21_i)) + 1
1436    
1437                tofmask(ch32a(tof32_i),hb32a(tof32_i)) =
1438         $           tofmask(ch32a(tof32_i),hb32a(tof32_i)) + 1
1439                tofmask(ch32b(tof32_i),hb32b(tof32_i)) =
1440         $           tofmask(ch32b(tof32_i),hb32b(tof32_i)) + 1
1441    
1442    C-------
1443    
1444             ENDIF
1445        ENDIF        ENDIF
1446    
1447  C     S22 - S31  C     S22 - S31
1448        IF (tof22_i.GT.none_find.AND.tof31_i.GT.none_find) THEN  
1449           xhelp1 = tof22(1,tof22_i,itdc)+tof22(2,tof22_i,itdc)        dist = ZTOF(4) - ZTOF(5)
1450           xhelp2 = tof31(1,tof31_i,itdc)+tof31(2,tof31_i,itdc)        dl = 0.
1451           ds = xhelp1-xhelp2        DO I=4,5
1452           ihelp=(tof22_i-1)*3+tof31_i           dl = dl + TLOUT(i)
1453           c1 = k_S22S31(1,ihelp)        ENDDO
1454           c2 = k_S22S31(2,ihelp)        F = dl/dist
1455           if ((xhelp1.lt.8000.).and.(xhelp2.lt.8000))  
1456       &   beta_a(7) = c2/(cos(theta13)*(ds-c1))  C      IF (tof22_i.GT.none_find.AND.tof31_i.GT.none_find) THEN
1457           IF ((tof22_i.GT.none_find).AND.(tof31_i.GT.none_find).AND.
1458         &    (ytofpos(2).NE.101.).AND.(ytofpos(3).NE.101.)) THEN
1459             t1 = tof22(1,tof22_i,itdc)
1460             t2 = tof22(2,tof22_i,itdc)
1461             t3 = tof31(1,tof31_i,itdc)
1462             t4 = tof31(2,tof31_i,itdc)
1463             IF ((t1.lt.4095).and.(t2.lt.4095).and.
1464         &        (t3.lt.4095).and.(t4.lt.4095)) THEN
1465                xhelp1 = tof22(1,tof22_i,itdc)+tof22(2,tof22_i,itdc)
1466                xhelp2 = tof31(1,tof31_i,itdc)+tof31(2,tof31_i,itdc)
1467                ds = xhelp1-xhelp2
1468                ihelp=(tof22_i-1)*3+tof31_i
1469                c1 = k_S22S31(1,ihelp)
1470                c2 = k_S22S31(2,ihelp)
1471                beta_a(7) = c2*F/(ds-c1)
1472    
1473    C-------ToF Mask - S22 - S31
1474    
1475                tofmask(ch22a(tof22_i),hb22a(tof22_i)) =
1476         $           tofmask(ch22a(tof22_i),hb22a(tof22_i)) + 1
1477                tofmask(ch22b(tof22_i),hb22b(tof22_i)) =
1478         $           tofmask(ch22b(tof22_i),hb22b(tof22_i)) + 1
1479    
1480                tofmask(ch31a(tof31_i),hb31a(tof31_i)) =
1481         $           tofmask(ch31a(tof31_i),hb31a(tof31_i)) + 1
1482                tofmask(ch31b(tof31_i),hb31b(tof31_i)) =
1483         $           tofmask(ch31b(tof31_i),hb31b(tof31_i)) + 1
1484    
1485    C-------
1486    
1487             ENDIF
1488        ENDIF        ENDIF
1489          
1490  C     S22 - S32  C     S22 - S32
1491        IF (tof22_i.GT.none_find.AND.tof32_i.GT.none_find) THEN  
1492           xhelp1 = tof22(1,tof22_i,itdc)+tof22(2,tof22_i,itdc)        dist = ZTOF(4) - ZTOF(6)
1493           xhelp2 = tof32(1,tof32_i,itdc)+tof32(2,tof32_i,itdc)        dl = 0.
1494           ds = xhelp1-xhelp2        DO I=4,6
1495           ihelp=(tof22_i-1)*3+tof32_i           dl = dl + TLOUT(i)
1496           c1 = k_S22S32(1,ihelp)        ENDDO
1497           c2 = k_S22S32(2,ihelp)        F = dl/dist
1498           if ((xhelp1.lt.8000.).and.(xhelp2.lt.8000))  
1499       &   beta_a(8) = c2/(cos(theta13)*(ds-c1))  C      IF (tof22_i.GT.none_find.AND.tof32_i.GT.none_find) THEN
1500           IF ((tof22_i.GT.none_find).AND.(tof32_i.GT.none_find).AND.
1501         &    (ytofpos(2).NE.101.).AND.(xtofpos(3).NE.101.)) THEN
1502             t1 = tof22(1,tof22_i,itdc)
1503             t2 = tof22(2,tof22_i,itdc)
1504             t3 = tof32(1,tof32_i,itdc)
1505             t4 = tof32(2,tof32_i,itdc)
1506             IF ((t1.lt.4095).and.(t2.lt.4095).and.
1507         &        (t3.lt.4095).and.(t4.lt.4095)) THEN
1508                xhelp1 = tof22(1,tof22_i,itdc)+tof22(2,tof22_i,itdc)
1509                xhelp2 = tof32(1,tof32_i,itdc)+tof32(2,tof32_i,itdc)
1510                ds = xhelp1-xhelp2
1511                ihelp=(tof22_i-1)*3+tof32_i
1512                c1 = k_S22S32(1,ihelp)
1513                c2 = k_S22S32(2,ihelp)
1514                beta_a(8) = c2*F/(ds-c1)
1515    
1516    C-------ToF Mask - S22 - S32
1517    
1518                tofmask(ch22a(tof22_i),hb22a(tof22_i)) =
1519         $           tofmask(ch22a(tof22_i),hb22a(tof22_i)) + 1
1520                tofmask(ch22b(tof22_i),hb22b(tof22_i)) =
1521         $           tofmask(ch22b(tof22_i),hb22b(tof22_i)) + 1
1522    
1523                tofmask(ch32a(tof32_i),hb32a(tof32_i)) =
1524         $           tofmask(ch32a(tof32_i),hb32a(tof32_i)) + 1
1525                tofmask(ch32b(tof32_i),hb32b(tof32_i)) =
1526         $           tofmask(ch32b(tof32_i),hb32b(tof32_i)) + 1
1527    
1528    C-------
1529    
1530             ENDIF
1531        ENDIF        ENDIF
1532    
1533  C     S11 - S21  C     S11 - S21
1534        IF (tof11_i.GT.none_find.AND.tof21_i.GT.none_find) THEN  
1535           xhelp1 = tof11(1,tof11_i,itdc)+tof11(2,tof11_i,itdc)        dist = ZTOF(1) - ZTOF(3)
1536           xhelp2 = tof21(1,tof21_i,itdc)+tof21(2,tof21_i,itdc)        dl = 0.
1537           ds = xhelp1-xhelp2        DO I=1,3
1538           ihelp=(tof11_i-1)*2+tof21_i           dl = dl + TLOUT(i)
1539           c1 = k_S11S21(1,ihelp)        ENDDO
1540           c2 = k_S11S21(2,ihelp)        F = dl/dist
1541           if ((xhelp1.lt.8000.).and.(xhelp2.lt.8000))  
1542       &   beta_a(9) = c2/(cos(theta12)*(ds-c1))  C      IF (tof11_i.GT.none_find.AND.tof21_i.GT.none_find) THEN
1543           IF ((tof11_i.GT.none_find).AND.(tof21_i.GT.none_find).AND.
1544         &    (ytofpos(1).NE.101.).AND.(xtofpos(2).NE.101.)) THEN
1545             t1 = tof11(1,tof11_i,itdc)
1546             t2 = tof11(2,tof11_i,itdc)
1547             t3 = tof21(1,tof21_i,itdc)
1548             t4 = tof21(2,tof21_i,itdc)
1549             IF ((t1.lt.4095).and.(t2.lt.4095).and.
1550         &        (t3.lt.4095).and.(t4.lt.4095)) THEN
1551                xhelp1 = tof11(1,tof11_i,itdc)+tof11(2,tof11_i,itdc)
1552                xhelp2 = tof21(1,tof21_i,itdc)+tof21(2,tof21_i,itdc)
1553                ds = xhelp1-xhelp2
1554                ihelp=(tof11_i-1)*2+tof21_i
1555                c1 = k_S11S21(1,ihelp)
1556                c2 = k_S11S21(2,ihelp)
1557                beta_a(9) = c2*F/(ds-c1)
1558    
1559    C-------ToF Mask - S11 - S21
1560    
1561                tofmask(ch11a(tof11_i),hb11a(tof11_i)) =
1562         $           tofmask(ch11a(tof11_i),hb11a(tof11_i)) + 1
1563                tofmask(ch11b(tof11_i),hb11b(tof11_i)) =
1564         $           tofmask(ch11b(tof11_i),hb11b(tof11_i)) + 1
1565    
1566                tofmask(ch21a(tof21_i),hb21a(tof21_i)) =
1567         $           tofmask(ch21a(tof21_i),hb21a(tof21_i)) + 1
1568                tofmask(ch21b(tof21_i),hb21b(tof21_i)) =
1569         $           tofmask(ch21b(tof21_i),hb21b(tof21_i)) + 1
1570    
1571    C-------
1572    
1573             ENDIF
1574        ENDIF        ENDIF
1575          
1576  C     S11 - S22  C     S11 - S22
1577        IF (tof11_i.GT.none_find.AND.tof22_i.GT.none_find) THEN  
1578           xhelp1 = tof11(1,tof11_i,itdc)+tof11(2,tof11_i,itdc)        dist = ZTOF(1) - ZTOF(4)
1579           xhelp2 = tof22(1,tof22_i,itdc)+tof22(2,tof22_i,itdc)        dl = 0.
1580           ds = xhelp1-xhelp2        DO I=1,4
1581           ihelp=(tof11_i-1)*2+tof22_i           dl = dl + TLOUT(i)
1582           c1 = k_S11S22(1,ihelp)        ENDDO
1583           c2 = k_S11S22(2,ihelp)        F = dl/dist
1584           if ((xhelp1.lt.8000.).and.(xhelp2.lt.8000))  
1585       &   beta_a(10) = c2/(cos(theta12)*(ds-c1))  C      IF (tof11_i.GT.none_find.AND.tof22_i.GT.none_find) THEN
1586           IF ((tof11_i.GT.none_find).AND.(tof22_i.GT.none_find).AND.
1587         &    (ytofpos(1).NE.101.).AND.(ytofpos(2).NE.101.)) THEN
1588             t1 = tof11(1,tof11_i,itdc)
1589             t2 = tof11(2,tof11_i,itdc)
1590             t3 = tof22(1,tof22_i,itdc)
1591             t4 = tof22(2,tof22_i,itdc)
1592             IF ((t1.lt.4095).and.(t2.lt.4095).and.
1593         &        (t3.lt.4095).and.(t4.lt.4095)) THEN
1594                xhelp1 = tof11(1,tof11_i,itdc)+tof11(2,tof11_i,itdc)
1595                xhelp2 = tof22(1,tof22_i,itdc)+tof22(2,tof22_i,itdc)
1596                ds = xhelp1-xhelp2
1597                ihelp=(tof11_i-1)*2+tof22_i
1598                c1 = k_S11S22(1,ihelp)
1599                c2 = k_S11S22(2,ihelp)
1600                beta_a(10) = c2*F/(ds-c1)
1601    
1602    C-------ToF Mask - S11 - S22
1603    
1604                tofmask(ch11a(tof11_i),hb11a(tof11_i)) =
1605         $           tofmask(ch11a(tof11_i),hb11a(tof11_i)) + 1
1606                tofmask(ch11b(tof11_i),hb11b(tof11_i)) =
1607         $           tofmask(ch11b(tof11_i),hb11b(tof11_i)) + 1
1608    
1609                tofmask(ch22a(tof22_i),hb22a(tof22_i)) =
1610         $           tofmask(ch22a(tof22_i),hb22a(tof22_i)) + 1
1611                tofmask(ch22b(tof22_i),hb22b(tof22_i)) =
1612         $           tofmask(ch22b(tof22_i),hb22b(tof22_i)) + 1
1613    
1614    C-------
1615    
1616             ENDIF
1617        ENDIF        ENDIF
1618    
1619  C     S12 - S21  C     S12 - S21
1620        IF (tof12_i.GT.none_find.AND.tof21_i.GT.none_find) THEN  
1621           xhelp1 = tof12(1,tof12_i,itdc)+tof12(2,tof12_i,itdc)        dist = ZTOF(2) - ZTOF(3)
1622           xhelp2 = tof21(1,tof21_i,itdc)+tof21(2,tof21_i,itdc)        dl = 0.
1623           ds = xhelp1-xhelp2        DO I=2,3
1624           ihelp=(tof12_i-1)*2+tof21_i           dl = dl + TLOUT(i)
1625           c1 = k_S12S21(1,ihelp)        ENDDO
1626           c2 = k_S12S21(2,ihelp)        F = dl/dist
1627           if ((xhelp1.lt.8000.).and.(xhelp2.lt.8000))  
1628       &   beta_a(11) = c2/(cos(theta12)*(ds-c1))  C      IF (tof12_i.GT.none_find.AND.tof21_i.GT.none_find) THEN
1629           IF ((tof12_i.GT.none_find).AND.(tof21_i.GT.none_find).AND.
1630         &    (xtofpos(1).NE.101.).AND.(xtofpos(2).NE.101.)) THEN
1631             t1 = tof12(1,tof12_i,itdc)
1632             t2 = tof12(2,tof12_i,itdc)
1633             t3 = tof21(1,tof21_i,itdc)
1634             t4 = tof21(2,tof21_i,itdc)
1635             IF ((t1.lt.4095).and.(t2.lt.4095).and.
1636         &        (t3.lt.4095).and.(t4.lt.4095)) THEN
1637                xhelp1 = tof12(1,tof12_i,itdc)+tof12(2,tof12_i,itdc)
1638                xhelp2 = tof21(1,tof21_i,itdc)+tof21(2,tof21_i,itdc)
1639                ds = xhelp1-xhelp2
1640                ihelp=(tof12_i-1)*2+tof21_i
1641                c1 = k_S12S21(1,ihelp)
1642                c2 = k_S12S21(2,ihelp)
1643                beta_a(11) = c2*F/(ds-c1)
1644    
1645    C-------ToF Mask - S12 - S21
1646    
1647                tofmask(ch12a(tof12_i),hb12a(tof12_i)) =
1648         $           tofmask(ch12a(tof12_i),hb12a(tof12_i)) + 1
1649                tofmask(ch12b(tof12_i),hb12b(tof12_i)) =
1650         $           tofmask(ch12b(tof12_i),hb12b(tof12_i)) + 1
1651    
1652                tofmask(ch21a(tof21_i),hb21a(tof21_i)) =
1653         $           tofmask(ch21a(tof21_i),hb21a(tof21_i)) + 1
1654                tofmask(ch21b(tof21_i),hb21b(tof21_i)) =
1655         $           tofmask(ch21b(tof21_i),hb21b(tof21_i)) + 1
1656    
1657    C-------
1658    
1659             ENDIF
1660        ENDIF        ENDIF
1661    
1662  C     S12 - S22  C     S12 - S22
1663        IF (tof12_i.GT.none_find.AND.tof22_i.GT.none_find) THEN  
1664           xhelp1 = tof12(1,tof12_i,itdc)+tof12(2,tof12_i,itdc)        dist = ZTOF(2) - ZTOF(4)
1665           xhelp2 = tof22(1,tof22_i,itdc)+tof22(2,tof22_i,itdc)        dl = 0.
1666           ds = xhelp1-xhelp2        DO I=2,4
1667           ihelp=(tof12_i-1)*2+tof22_i           dl = dl + TLOUT(i)
1668           c1 = k_S12S22(1,ihelp)        ENDDO
1669           c2 = k_S12S22(2,ihelp)        F = dl/dist
1670           if ((xhelp1.lt.8000.).and.(xhelp2.lt.8000))  
1671       &   beta_a(12) = c2/(cos(theta12)*(ds-c1))  C      IF (tof12_i.GT.none_find.AND.tof22_i.GT.none_find) THEN
1672           IF ((tof12_i.GT.none_find).AND.(tof22_i.GT.none_find).AND.
1673         &    (xtofpos(1).NE.101.).AND.(ytofpos(2).NE.101.)) THEN
1674             t1 = tof12(1,tof12_i,itdc)
1675             t2 = tof12(2,tof12_i,itdc)
1676             t3 = tof22(1,tof22_i,itdc)
1677             t4 = tof22(2,tof22_i,itdc)
1678             IF ((t1.lt.4095).and.(t2.lt.4095).and.
1679         &        (t3.lt.4095).and.(t4.lt.4095)) THEN
1680                xhelp1 = tof12(1,tof12_i,itdc)+tof12(2,tof12_i,itdc)
1681                xhelp2 = tof22(1,tof22_i,itdc)+tof22(2,tof22_i,itdc)
1682                ds = xhelp1-xhelp2
1683                ihelp=(tof12_i-1)*2+tof22_i
1684                c1 = k_S12S22(1,ihelp)
1685                c2 = k_S12S22(2,ihelp)
1686                beta_a(12) = c2*F/(ds-c1)
1687    
1688    C-------ToF Mask - S12 - S22
1689    
1690                tofmask(ch12a(tof12_i),hb12a(tof12_i)) =
1691         $           tofmask(ch12a(tof12_i),hb12a(tof12_i)) + 1
1692                tofmask(ch12b(tof12_i),hb12b(tof12_i)) =
1693         $           tofmask(ch12b(tof12_i),hb12b(tof12_i)) + 1
1694    
1695                tofmask(ch22a(tof22_i),hb22a(tof22_i)) =
1696         $           tofmask(ch22a(tof22_i),hb22a(tof22_i)) + 1
1697                tofmask(ch22b(tof22_i),hb22b(tof22_i)) =
1698         $           tofmask(ch22b(tof22_i),hb22b(tof22_i)) + 1
1699    
1700    C-------
1701    
1702             ENDIF
1703        ENDIF        ENDIF
1704          
1705  C-------  C-------
1706    
1707        icount=0        icount=0
# Line 541  C------- Line 1714  C-------
1714              icount= icount+1              icount= icount+1
1715              if (i.le.4) w_i=1./(0.13**2.)              if (i.le.4) w_i=1./(0.13**2.)
1716              if ((i.ge.5).and.(i.le.8)) w_i=1./(0.16**2.)              if ((i.ge.5).and.(i.le.8)) w_i=1./(0.16**2.)
1717              if (i.ge.9) w_i=1./(0.25**2.)     ! to be checked              if (i.ge.9) w_i=1./(0.25**2.) ! to be checked
1718              sxw=sxw + beta_a(i)*w_i              sxw=sxw + beta_a(i)*w_i
1719              sw =sw + w_i              sw =sw + w_i
1720           endif           endif
1721        enddo        enddo
1722          
1723        if (icount.gt.0) beta_mean=sxw/sw        if (icount.gt.0) beta_mean=sxw/sw
1724        beta_a(13) = beta_mean        beta_a(13) = beta_mean
1725    
1726    
1727    c       IF (tof11_i.GT.none_find)
1728    c     &  write(*,*) '11 ',tof11(1,tof11_i,itdc),tof11(2,tof11_i,itdc)
1729    c       IF (tof12_i.GT.none_find)
1730    c     &  write(*,*) '12 ',tof12(1,tof12_i,itdc),tof12(2,tof12_i,itdc)
1731    
1732    c       IF (tof21_i.GT.none_find)
1733    c     &  write(*,*) '21 ',tof21(1,tof21_i,itdc),tof21(2,tof21_i,itdc)
1734    c       IF (tof22_i.GT.none_find)
1735    c     &  write(*,*) '22 ',tof22(1,tof22_i,itdc),tof22(2,tof22_i,itdc)
1736    
1737    c       IF (tof31_i.GT.none_find)
1738    c     &  write(*,*) '31 ',tof31(1,tof31_i,itdc),tof31(2,tof31_i,itdc)
1739    c       IF (tof32_i.GT.none_find)
1740    c     &  write(*,*) '32 ',tof32(1,tof32_i,itdc),tof32(2,tof32_i,itdc)
1741    
1742    c      write(*,*) xtofpos
1743    c      write(*,*) ytofpos
1744    C      write(*,*)'toftrk beta', beta_a
1745    C      write(*,*) adcflagtof
1746    C     write(*,*)'TOFTRK ',ievent,beta_a(1),beta_a(2),beta_a(3),beta_a(4)
1747    c      write(*,*) 'toftrk'
1748    c      write(*,*) xtofpos
1749    c      write(*,*) ytofpos
1750    c      write(*,*) xtr_tof
1751    c      write(*,*) ytr_tof
1752    
1753    
1754    
1755        RETURN        RETURN
1756        END        END
1757    
1758    
1759    
1760    
1761    C------------------------------------------------------------------
1762    
Legend:
Removed from v.1.1  
changed lines
  Added in v.1.12
  ViewVC Help
Powered by ViewVC 1.1.23