/[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.11 by mocchiut, Mon Apr 30 15:46:31 2007 UTC revision 1.12 by pam-de, Mon Aug 20 14:22:33 2007 UTC
# Line 7  C     Changed to use DOTRACK2 Line 7  C     Changed to use DOTRACK2
7  C     Beta calculation: now the flightpath (instead of cos(theta)) is used  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  C     Beta calculation: all 4 TDV measurements must be < 4095 (in the old
9  C     routine it was (t1+t2)<8000  C     routine it was (t1+t2)<8000
10  C      C
11  C     08-12-06 WM:  C     08-12-06 WM:
12  C     adc_c-bug :  The raw ADC value was multiplied with cos(theta)  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"  C     and AFTER that there was an if statement "if tof32(right,i,iadc) < 4095"
14  C      C
15  C  jan-07 GF: ADC/TDCflags(4,12) inserted to flag artificial ADC/TDC  C  jan-07 GF: ADC/TDCflags(4,12) inserted to flag artificial ADC/TDC
16  C             values  C             values
17  C  jan-07 WM: artificial ADC values created using attenuation calibration  C  jan-07 WM: artificial ADC values created using attenuation calibration
18  C  jan-07 WM: artificial TDC values created using xy_coor calibration  C  jan-07 WM: artificial TDC values created using xy_coor calibration
# Line 20  C  jan-07 WM: modified xtofpos flag "101 Line 20  C  jan-07 WM: modified xtofpos flag "101
20  C             dimension of the paddle +/- 10 cm  C             dimension of the paddle +/- 10 cm
21  C  jan-07 WM: if xtofpos=101 then this paddle is not used for beta  C  jan-07 WM: if xtofpos=101 then this paddle is not used for beta
22  C             calculation  C             calculation
23  C  jan-07 WM: in the xtofpos calculation a check for TDC.ne.4095 was  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  C             inserted. In the old code one would still calculate a
25  C             xtofpos-value even if the TDC information was missing  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  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  C  jan-05 WM: bug fixed: calculation of zenith angles using DOTRACK2
# Line 30  C  jan-07 WM: bug fixed: in some cases t Line 30  C  jan-07 WM: bug fixed: in some cases t
30  C             leftover "xhelp" value  C             leftover "xhelp" value
31  C  apr-07 WM: attenuation fit curve is now a double exponential fit  C  apr-07 WM: attenuation fit curve is now a double exponential fit
32  C             conversion from raw ADC to pC using calibration function  C             conversion from raw ADC to pC using calibration function
33  C             new variables xtr_tof(6) and ytr_tof(6) give track position  C             new variables xtr_tof(6) and ytr_tof(6) give track position
34  C             at ToF layers  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  C
39  C****************************************************************************  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 66  c      REAL dx,dy,dr Line 69  c      REAL dx,dy,dr
69        REAL t1,t2,t3,t4        REAL t1,t2,t3,t4
70        REAL yhelp,xhelp,xhelp1,xhelp2        REAL yhelp,xhelp,xhelp1,xhelp2
71        REAL c1,c2,sw,sxw,w_i        REAL c1,c2,sw,sxw,w_i
72        REAL dist,dl,F        REAL dist,dl,F
73        INTEGER icount,ievent        INTEGER icount,ievent
74        REAL xhelp_a,xhelp_t        REAL xhelp_a,xhelp_t
75    
76        REAL beta_mean        REAL beta_mean
77        REAL hepratio        REAL hepratio
78          
79        INTEGER j        INTEGER j
80    
81        real atten,pc_adc        real atten,pc_adc
# Line 82  c      REAL dx,dy,dr Line 85  c      REAL dx,dy,dr
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,xpos        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 /        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  *******************************************************************  *******************************************************************
# Line 169  C  ratio helium to proton ca. 4 Line 172  C  ratio helium to proton ca. 4
172  C----------------------------------------------------------------------  C----------------------------------------------------------------------
173  C-------------------------get ToF data --------------------------------  C-------------------------get ToF data --------------------------------
174  C     we cannot use the tofxx(x,x,x)  data  from tofl2com since it is  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  C     manipulated (Time-walk, artificila ADc and TDC values using ToF
176  C     standalone information  C     standalone information
177  C----------------------------------------------------------------------  C----------------------------------------------------------------------
178    
# Line 281  C     COPY THE ALFA VECTOR FROM AL_PP TO Line 284  C     COPY THE ALFA VECTOR FROM AL_PP TO
284  c      write(*,*) AL_P  c      write(*,*) AL_P
285    
286        if (al_p(5).eq.0.) THEN        if (al_p(5).eq.0.) THEN
287  c         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 *** --------
# Line 292  C     call DOTRACK(NPTOF,ZIN,XOUT,YOUT,A Line 295  C     call DOTRACK(NPTOF,ZIN,XOUT,YOUT,A
295  C     write(*,*) (TLOUT(i),i=1,6)  C     write(*,*) (TLOUT(i),i=1,6)
296    
297        if(IFAIL.ne.0)then        if(IFAIL.ne.0)then
298  c         print *,' TOF - WARNING F77: tracking failed '           print *,' TOF - WARNING F77: tracking failed '
299           goto 969           goto 969
300        endif        endif
301  *     ------------------------------------------  *     ------------------------------------------
# Line 302  c         print *,' TOF - WARNING F77: t Line 305  c         print *,' TOF - WARNING F77: t
305  C---  Fill xtr_tof  and ytr_tof: positions from tracker at ToF layers  C---  Fill xtr_tof  and ytr_tof: positions from tracker at ToF layers
306        do j=1,6        do j=1,6
307        xtr_tof(j) = XOUT(j)        xtr_tof(j) = XOUT(j)
308        ytr_tof(j) = YOUT(j)                ytr_tof(j) = YOUT(j)
309        enddo        enddo
310    
311    
# Line 368  C---- S116A TDC=819 Line 371  C---- S116A TDC=819
371               tof11(1,6,1) = 4095               tof11(1,6,1) = 4095
372               tdcflagtof(ch11a(6),hb11a(6))=2               tdcflagtof(ch11a(6),hb11a(6))=2
373         endif         endif
374      
375  C---- S222B TDC=819  C---- S222B TDC=819
376         if (tof22(2,2,1).EQ.819) then         if (tof22(2,2,1).EQ.819) then
377               tof22(2,2,1) = 4095               tof22(2,2,1) = 4095
378               tdcflagtof(ch22b(2),hb22b(2))=2               tdcflagtof(ch22b(2),hb22b(2))=2
379         endif         endif
380    
381  C-------------------------------------------------------------  C-------------------------------------------------------------
382  C-------check which paddle penetrated the track  -----------  C-------check which paddle penetrated the track  -----------
383  C-------------------------------------------------------------  C-------------------------------------------------------------
# Line 385  c     DATA tof21_y/  3.75,-3.75/  ! padd Line 388  c     DATA tof21_y/  3.75,-3.75/  ! padd
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
# Line 416  C--------------S11 --------------------- Line 419  C--------------S11 ---------------------
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    
# Line 433  C--------------S12 --------------------- Line 436  C--------------S12 ---------------------
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    
# Line 450  C--------------S21 --------------------- Line 453  C--------------S21 ---------------------
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    
# Line 467  C--------------S22 --------------------- Line 470  C--------------S22 ---------------------
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    
# Line 484  C--------------S31 --------------------- Line 487  C--------------S31 ---------------------
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    
# Line 511  C--------------------------------------- Line 514  C---------------------------------------
514  C-----------------------S11 -----------------------------------------  C-----------------------S11 -----------------------------------------
515    
516        IF (tof11_i.GT.none_find) THEN        IF (tof11_i.GT.none_find) THEN
517           xpos = yout(1)           xpos = yout(1)
518           i = tof11_i           i = tof11_i
519           if ((tof11(1,tof11_i,itdc).EQ.4095).AND.           if ((tof11(1,tof11_i,itdc).EQ.4095).AND.
520       &        (tof11(2,tof11_i,itdc).LT.4095)) THEN         &        (tof11(2,tof11_i,itdc).LT.4095)) THEN
521    
522  c       write(*,*)'11lb4 ',i,tof11(1,tof11_i,itdc),tof11(2,tof11_i,itdc)  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)              tof11(1,tof11_i,itdc) = tof11(2,tof11_i,itdc)
525       &           + 2*(y_coor_lin11(tof11_i,offset)       &           + 2*(y_coor_lin11(tof11_i,offset)
526       &           + xpos*y_coor_lin11(tof11_i,slope))       &           + xpos*y_coor_lin11(tof11_i,slope))
527    
528  c       write(*,*)'11laf ',i,tof11(1,tof11_i,itdc),tof11(2,tof11_i,itdc)  c       write(*,*)'11laf ',i,tof11(1,tof11_i,itdc),tof11(2,tof11_i,itdc)
529    
530              tdcflag(ch11a(i),hb11a(i)) = 1              tdcflag(ch11a(i),hb11a(i)) = 1
531      
532          ENDIF          ENDIF
533           if ((tof11(2,tof11_i,itdc).EQ.4095).AND.           if ((tof11(2,tof11_i,itdc).EQ.4095).AND.
534       &        (tof11(1,tof11_i,itdc).LT.4095)) THEN         &        (tof11(1,tof11_i,itdc).LT.4095)) THEN
535    
536  c       write(*,*)'11rb4 ',i,tof11(1,tof11_i,itdc),tof11(2,tof11_i,itdc)  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)              tof11(2,tof11_i,itdc) = tof11(1,tof11_i,itdc)
539       &           - 2*(y_coor_lin11(tof11_i,offset)       &           - 2*(y_coor_lin11(tof11_i,offset)
540       &           + xpos*y_coor_lin11(tof11_i,slope))       &           + xpos*y_coor_lin11(tof11_i,slope))
541  c       write(*,*)'11raf ',i,tof11(1,tof11_i,itdc),tof11(2,tof11_i,itdc)  c       write(*,*)'11raf ',i,tof11(1,tof11_i,itdc),tof11(2,tof11_i,itdc)
542    
543              tdcflag(ch11b(i),hb11b(i)) = 1              tdcflag(ch11b(i),hb11b(i)) = 1
544           ENDIF               ENDIF
545        ENDIF        ENDIF
546    
547  C-----------------------S12 -----------------------------------------  C-----------------------S12 -----------------------------------------
548    
549        IF (tof12_i.GT.none_find) THEN        IF (tof12_i.GT.none_find) THEN
550           xpos = xout(2)           xpos = xout(2)
551           i = tof12_i           i = tof12_i
552           if ((tof12(1,tof12_i,itdc).EQ.4095).AND.           if ((tof12(1,tof12_i,itdc).EQ.4095).AND.
553       &        (tof12(2,tof12_i,itdc).LT.4095)) THEN       &        (tof12(2,tof12_i,itdc).LT.4095)) THEN
# Line 561  C-----------------------S12 ------------ Line 564  C-----------------------S12 ------------
564              tdcflag(ch12b(i),hb12b(i)) = 1              tdcflag(ch12b(i),hb12b(i)) = 1
565           ENDIF           ENDIF
566        ENDIF        ENDIF
567          
568  C-----------------------S21 -----------------------------------------  C-----------------------S21 -----------------------------------------
569    
570        IF (tof21_i.GT.none_find) THEN        IF (tof21_i.GT.none_find) THEN
571           xpos = xout(3)           xpos = xout(3)
572           i = tof21_i           i = tof21_i
573           if ((tof21(1,tof21_i,itdc).EQ.4095).AND.           if ((tof21(1,tof21_i,itdc).EQ.4095).AND.
574       &        (tof21(2,tof21_i,itdc).LT.4095)) THEN       &        (tof21(2,tof21_i,itdc).LT.4095)) THEN
# Line 586  C-----------------------S21 ------------ Line 589  C-----------------------S21 ------------
589  C-----------------------S22 -----------------------------------------  C-----------------------S22 -----------------------------------------
590    
591        IF (tof22_i.GT.none_find) THEN        IF (tof22_i.GT.none_find) THEN
592           xpos = yout(4)           xpos = yout(4)
593           i = tof22_i           i = tof22_i
594           if ((tof22(1,tof22_i,itdc).EQ.4095).AND.           if ((tof22(1,tof22_i,itdc).EQ.4095).AND.
595       &        (tof22(2,tof22_i,itdc).LT.4095)) THEN       &        (tof22(2,tof22_i,itdc).LT.4095)) THEN
# Line 605  C-----------------------S22 ------------ Line 608  C-----------------------S22 ------------
608        ENDIF        ENDIF
609    
610  C-----------------------S31 -----------------------------------------  C-----------------------S31 -----------------------------------------
611          
612        IF (tof31_i.GT.none_find) THEN        IF (tof31_i.GT.none_find) THEN
613           xpos = yout(5)           xpos = yout(5)
614           i = tof31_i           i = tof31_i
615           if ((tof31(1,tof31_i,itdc).EQ.4095).AND.           if ((tof31(1,tof31_i,itdc).EQ.4095).AND.
616       &        (tof31(2,tof31_i,itdc).LT.4095)) THEN       &        (tof31(2,tof31_i,itdc).LT.4095)) THEN
# Line 628  C-----------------------S31 ------------ Line 631  C-----------------------S31 ------------
631  C-----------------------S32 -----------------------------------------  C-----------------------S32 -----------------------------------------
632    
633        IF (tof32_i.GT.none_find) THEN        IF (tof32_i.GT.none_find) THEN
634           xpos = xout(6)           xpos = xout(6)
635           i = tof32_i           i = tof32_i
636           if ((tof32(1,tof32_i,itdc).EQ.4095).AND.           if ((tof32(1,tof32_i,itdc).EQ.4095).AND.
637       &        (tof32(2,tof32_i,itdc).LT.4095)) THEN       &        (tof32(2,tof32_i,itdc).LT.4095)) THEN
# Line 647  C-----------------------S32 ------------ Line 650  C-----------------------S32 ------------
650        ENDIF        ENDIF
651    
652  C--------------------------------------------------------------------  C--------------------------------------------------------------------
653  C---- if TDCleft.and.TDCright and NO ADC insert artificial ADC  C---- if paddle hit: if we have TDC value but no ADC, create ADC value
 C---- values  
654  C--------------------------------------------------------------------  C--------------------------------------------------------------------
655  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
656  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/
# Line 663  C----------------------------S1 -------- Line 665  C----------------------------S1 --------
665        yhelp=yout(1)        yhelp=yout(1)
666        IF (tof11_i.GT.none_find.AND.abs(yout(1)).lt.100) THEN        IF (tof11_i.GT.none_find.AND.abs(yout(1)).lt.100) THEN
667           i = tof11_i           i = tof11_i
668  c         if (tof11(left,i,iadc).eq.4095) then           if ((tof11(left,tof11_i,itdc).LT.4095).AND.
669           if (adc(ch11a(i),hb11a(i)).eq.4095) then       &       (adc(ch11a(i),hb11a(i)).eq.4095)) then
670              phi   = atan(tan(THYOUT(1))/tan(THXOUT(1)))              phi   = atan(tan(THYOUT(1))/tan(THXOUT(1)))
671              theta = atan(tan(THXOUT(1))/cos(phi))              theta = atan(tan(THXOUT(1))/cos(phi))
672  c            xkorr=adcx11(left,i,1)*exp(-yhelp/adcx11(left,i,2))  c            xkorr=adcx11(left,i,1)*exp(-yhelp/adcx11(left,i,2))
673              xkorr = atten(left,11,i,yhelp)              xkorr = atten(left,11,i,yhelp)
674              xkorr=xkorr/hepratio                          xkorr=xkorr/hepratio
675              tof11(left,i,iadc)=xkorr/cos(theta)              tof11(left,i,iadc)=xkorr/cos(theta)
676              adcflag(ch11a(i),hb11a(i)) = 1              adcflag(ch11a(i),hb11a(i)) = 1
677           endif           endif
678  c         if (tof11(right,i,iadc).eq.4095) then           if ((tof11(right,tof11_i,itdc).LT.4095).AND.
679           if (adc(ch11b(i),hb11b(i)).eq.4095) then       &       (adc(ch11b(i),hb11b(i)).eq.4095)) then
680              phi   = atan(tan(THYOUT(1))/tan(THXOUT(1)))              phi   = atan(tan(THYOUT(1))/tan(THXOUT(1)))
681              theta = atan(tan(THXOUT(1))/cos(phi))              theta = atan(tan(THXOUT(1))/cos(phi))
682  c            xkorr=adcx11(right,i,1)*exp(yhelp/adcx11(right,i,2))  c            xkorr=adcx11(right,i,1)*exp(yhelp/adcx11(right,i,2))
683              xkorr = atten(right,11,i,yhelp)              xkorr = atten(right,11,i,yhelp)
684              xkorr=xkorr/hepratio                          xkorr=xkorr/hepratio
685              tof11(right,i,iadc)=xkorr/cos(theta)              tof11(right,i,iadc)=xkorr/cos(theta)
686              adcflag(ch11b(i),hb11b(i)) = 1              adcflag(ch11b(i),hb11b(i)) = 1
687           endif           endif
# Line 688  c            xkorr=adcx11(right,i,1)*exp Line 690  c            xkorr=adcx11(right,i,1)*exp
690        xhelp=xout(2)        xhelp=xout(2)
691        IF (tof12_i.GT.none_find.AND.abs(xout(2)).lt.100) THEN        IF (tof12_i.GT.none_find.AND.abs(xout(2)).lt.100) THEN
692           i = tof12_i           i = tof12_i
693  c         if (tof12(left,i,iadc).eq.4095) then           if ((tof12(left,tof12_i,itdc).LT.4095).AND.
694           if (adc(ch12a(i),hb12a(i)).eq.4095) then       &       (adc(ch12a(i),hb12a(i)).eq.4095)) then
695              phi   = atan(tan(THYOUT(2))/tan(THXOUT(2)))              phi   = atan(tan(THYOUT(2))/tan(THXOUT(2)))
696              theta = atan(tan(THXOUT(2))/cos(phi))              theta = atan(tan(THXOUT(2))/cos(phi))
697  c            xkorr=adcx12(left,i,1)*exp(-xhelp/adcx12(left,i,2))  c            xkorr=adcx12(left,i,1)*exp(-xhelp/adcx12(left,i,2))
698              xkorr = atten(left,12,i,xhelp)              xkorr = atten(left,12,i,xhelp)
699              xkorr=xkorr/hepratio                          xkorr=xkorr/hepratio
700              tof12(left,i,iadc) = xkorr/cos(theta)              tof12(left,i,iadc) = xkorr/cos(theta)
701              adcflag(ch12a(i),hb12a(i)) = 1              adcflag(ch12a(i),hb12a(i)) = 1
702           endif           endif
703  c         if (tof12(right,i,iadc).eq.4095) then           if ((tof12(right,tof12_i,itdc).LT.4095).AND.
704           if (adc(ch12b(i),hb12b(i)).eq.4095) then       &       (adc(ch12b(i),hb12b(i)).eq.4095)) then
705              phi   = atan(tan(THYOUT(2))/tan(THXOUT(2)))              phi   = atan(tan(THYOUT(2))/tan(THXOUT(2)))
706              theta = atan(tan(THXOUT(2))/cos(phi))              theta = atan(tan(THXOUT(2))/cos(phi))
707  c            xkorr=adcx12(right,i,1)*exp(xhelp/adcx12(right,i,2))  c            xkorr=adcx12(right,i,1)*exp(xhelp/adcx12(right,i,2))
708              xkorr = atten(right,12,i,xhelp)              xkorr = atten(right,12,i,xhelp)
709              xkorr=xkorr/hepratio                          xkorr=xkorr/hepratio
710              tof12(right,i,iadc) = xkorr/cos(theta)              tof12(right,i,iadc) = xkorr/cos(theta)
711              adcflag(ch12b(i),hb12b(i)) = 1              adcflag(ch12b(i),hb12b(i)) = 1
712           endif           endif
# Line 715  C-----------------------------S2 ------- Line 717  C-----------------------------S2 -------
717        xhelp=xout(3)        xhelp=xout(3)
718        IF (tof21_i.GT.none_find.AND.abs(xout(3)).lt.100) THEN        IF (tof21_i.GT.none_find.AND.abs(xout(3)).lt.100) THEN
719           i = tof21_i           i = tof21_i
720  c         if (tof21(left,i,iadc).eq.4095) then           if ((tof21(left,tof21_i,itdc).LT.4095).AND.
721           if (adc(ch21a(i),hb21a(i)).eq.4095) then       &       (adc(ch21a(i),hb21a(i)).eq.4095)) then
722              phi   = atan(tan(THYOUT(3))/tan(THXOUT(3)))              phi   = atan(tan(THYOUT(3))/tan(THXOUT(3)))
723              theta = atan(tan(THXOUT(3))/cos(phi))              theta = atan(tan(THXOUT(3))/cos(phi))
724  c            xkorr=adcx21(left,i,1)*exp(-xhelp/adcx21(left,i,2))  c            xkorr=adcx21(left,i,1)*exp(-xhelp/adcx21(left,i,2))
725              xkorr = atten(left,21,i,xhelp)              xkorr = atten(left,21,i,xhelp)
726              xkorr=xkorr/hepratio                          xkorr=xkorr/hepratio
727              tof21(left,i,iadc) = xkorr/cos(theta)              tof21(left,i,iadc) = xkorr/cos(theta)
728              adcflag(ch21a(i),hb21a(i)) = 1              adcflag(ch21a(i),hb21a(i)) = 1
729           endif           endif
730  c         if (tof21(right,i,iadc).eq.4095) then           if ((tof21(right,tof21_i,itdc).LT.4095).AND.
731           if (adc(ch21b(i),hb21b(i)).eq.4095) then       &       (adc(ch21b(i),hb21b(i)).eq.4095)) then
732              phi   = atan(tan(THYOUT(3))/tan(THXOUT(3)))              phi   = atan(tan(THYOUT(3))/tan(THXOUT(3)))
733              theta = atan(tan(THXOUT(3))/cos(phi))              theta = atan(tan(THXOUT(3))/cos(phi))
734  c            xkorr=adcx21(right,i,1)*exp(xhelp/adcx21(right,i,2))  c            xkorr=adcx21(right,i,1)*exp(xhelp/adcx21(right,i,2))
735              xkorr = atten(right,21,i,xhelp)              xkorr = atten(right,21,i,xhelp)
736              xkorr=xkorr/hepratio                          xkorr=xkorr/hepratio
737              tof21(right,i,iadc) = xkorr/cos(theta)              tof21(right,i,iadc) = xkorr/cos(theta)
738              adcflag(ch21b(i),hb21b(i)) = 1              adcflag(ch21b(i),hb21b(i)) = 1
739           endif           endif
# Line 741  c            xkorr=adcx21(right,i,1)*exp Line 743  c            xkorr=adcx21(right,i,1)*exp
743        yhelp=yout(4)        yhelp=yout(4)
744        IF (tof22_i.GT.none_find.AND.abs(yout(4)).lt.100) THEN        IF (tof22_i.GT.none_find.AND.abs(yout(4)).lt.100) THEN
745           i = tof22_i           i = tof22_i
746  c         if (tof22(left,i,iadc).eq.4095) then           if ((tof22(left,tof22_i,itdc).LT.4095).AND.
747           if (adc(ch22a(i),hb22a(i)).eq.4095) then       &       (adc(ch22a(i),hb22a(i)).eq.4095)) then
748              phi   = atan(tan(THYOUT(4))/tan(THXOUT(4)))              phi   = atan(tan(THYOUT(4))/tan(THXOUT(4)))
749              theta = atan(tan(THXOUT(4))/cos(phi))              theta = atan(tan(THXOUT(4))/cos(phi))
750  c            xkorr=adcx22(left,i,1)*exp(-yhelp/adcx22(left,i,2))  c            xkorr=adcx22(left,i,1)*exp(-yhelp/adcx22(left,i,2))
751              xkorr = atten(left,22,i,yhelp)              xkorr = atten(left,22,i,yhelp)
752              xkorr=xkorr/hepratio                          xkorr=xkorr/hepratio
753              tof22(left,i,iadc) = xkorr/cos(theta)              tof22(left,i,iadc) = xkorr/cos(theta)
754              adcflag(ch22a(i),hb22a(i)) = 1              adcflag(ch22a(i),hb22a(i)) = 1
755           endif           endif
756  c         if (tof22(right,i,iadc).eq.4095) then           if ((tof22(right,tof22_i,itdc).LT.4095).AND.
757           if (adc(ch22b(i),hb22b(i)).eq.4095) then       &       (adc(ch22b(i),hb22b(i)).eq.4095)) then
758              phi   = atan(tan(THYOUT(4))/tan(THXOUT(4)))              phi   = atan(tan(THYOUT(4))/tan(THXOUT(4)))
759              theta = atan(tan(THXOUT(4))/cos(phi))              theta = atan(tan(THXOUT(4))/cos(phi))
760  c            xkorr=adcx22(right,i,1)*exp(yhelp/adcx22(right,i,2))  c            xkorr=adcx22(right,i,1)*exp(yhelp/adcx22(right,i,2))
761              xkorr = atten(right,22,i,yhelp)              xkorr = atten(right,22,i,yhelp)
762              xkorr=xkorr/hepratio                          xkorr=xkorr/hepratio
763              tof22(right,i,iadc) = xkorr/cos(theta)              tof22(right,i,iadc) = xkorr/cos(theta)
764              adcflag(ch22b(i),hb22b(i)) = 1              adcflag(ch22b(i),hb22b(i)) = 1
765           endif           endif
# Line 768  C-----------------------------S3 ------- Line 770  C-----------------------------S3 -------
770        yhelp=yout(5)        yhelp=yout(5)
771        IF (tof31_i.GT.none_find.AND.abs(yout(5)).lt.100) THEN        IF (tof31_i.GT.none_find.AND.abs(yout(5)).lt.100) THEN
772           i = tof31_i           i = tof31_i
773  c         if (tof31(left,i,iadc).eq.4095) then           if ((tof31(left,tof31_i,itdc).LT.4095).AND.
774           if (adc(ch31a(i),hb31a(i)).eq.4095) then       &       (adc(ch31a(i),hb31a(i)).eq.4095)) then
775              phi   = atan(tan(THYOUT(5))/tan(THXOUT(5)))              phi   = atan(tan(THYOUT(5))/tan(THXOUT(5)))
776              theta = atan(tan(THXOUT(5))/cos(phi))              theta = atan(tan(THXOUT(5))/cos(phi))
777  c            xkorr=adcx31(left,i,1)*exp(-yhelp/adcx31(left,i,2))  c            xkorr=adcx31(left,i,1)*exp(-yhelp/adcx31(left,i,2))
778              xkorr = atten(left,31,i,yhelp)              xkorr = atten(left,31,i,yhelp)
779              xkorr=xkorr/hepratio                          xkorr=xkorr/hepratio
780              tof31(left,i,iadc) = xkorr/cos(theta)              tof31(left,i,iadc) = xkorr/cos(theta)
781              adcflag(ch31a(i),hb31a(i)) = 1              adcflag(ch31a(i),hb31a(i)) = 1
782           endif           endif
783  c         if (tof31(right,i,iadc).eq.4095) then           if ((tof31(right,tof31_i,itdc).LT.4095).AND.
784           if (adc(ch31b(i),hb31b(i)).eq.4095) then       &       (adc(ch31b(i),hb31b(i)).eq.4095)) then
785              phi   = atan(tan(THYOUT(5))/tan(THXOUT(5)))              phi   = atan(tan(THYOUT(5))/tan(THXOUT(5)))
786              theta = atan(tan(THXOUT(5))/cos(phi))              theta = atan(tan(THXOUT(5))/cos(phi))
787  c            xkorr=adcx31(right,i,1)*exp(yhelp/adcx31(right,i,2))  c            xkorr=adcx31(right,i,1)*exp(yhelp/adcx31(right,i,2))
788              xkorr = atten(right,31,i,yhelp)              xkorr = atten(right,31,i,yhelp)
789              xkorr=xkorr/hepratio                          xkorr=xkorr/hepratio
790              tof31(right,i,iadc) = xkorr/cos(theta)              tof31(right,i,iadc) = xkorr/cos(theta)
791              adcflag(ch31b(i),hb31b(i)) = 1              adcflag(ch31b(i),hb31b(i)) = 1
792           endif           endif
# Line 794  c            xkorr=adcx31(right,i,1)*exp Line 796  c            xkorr=adcx31(right,i,1)*exp
796        xhelp=xout(6)        xhelp=xout(6)
797        IF (tof32_i.GT.none_find.AND.abs(xout(6)).lt.100) THEN        IF (tof32_i.GT.none_find.AND.abs(xout(6)).lt.100) THEN
798           i = tof32_i           i = tof32_i
799  c         if (tof32(left,i,iadc).eq.4095) then           if ((tof32(left,tof32_i,itdc).LT.4095).AND.
800           if (adc(ch32a(i),hb32a(i)).eq.4095) then       &       (adc(ch32a(i),hb32a(i)).eq.4095)) then
801              phi   = atan(tan(THYOUT(6))/tan(THXOUT(6)))              phi   = atan(tan(THYOUT(6))/tan(THXOUT(6)))
802              theta = atan(tan(THXOUT(6))/cos(phi))              theta = atan(tan(THXOUT(6))/cos(phi))
803  c            xkorr=adcx32(left,i,1)*exp(-xhelp/adcx32(left,i,2))  c            xkorr=adcx32(left,i,1)*exp(-xhelp/adcx32(left,i,2))
804              xkorr = atten(left,32,i,xhelp)              xkorr = atten(left,32,i,xhelp)
805              xkorr=xkorr/hepratio                          xkorr=xkorr/hepratio
806              tof32(left,i,iadc) = xkorr/cos(theta)              tof32(left,i,iadc) = xkorr/cos(theta)
807              adcflag(ch32a(i),hb32a(i)) = 1              adcflag(ch32a(i),hb32a(i)) = 1
808           endif           endif
809  c         if (tof32(right,i,iadc).eq.4095) then           if ((tof32(right,tof32_i,itdc).LT.4095).AND.
810           if (adc(ch32b(i),hb32b(i)).eq.4095) then       &       (adc(ch32b(i),hb32b(i)).eq.4095)) then
811              phi   = atan(tan(THYOUT(6))/tan(THXOUT(6)))              phi   = atan(tan(THYOUT(6))/tan(THXOUT(6)))
812              theta = atan(tan(THXOUT(6))/cos(phi))              theta = atan(tan(THXOUT(6))/cos(phi))
813  c            xkorr=adcx32(right,i,1)*exp(xhelp/adcx32(right,i,2))  c            xkorr=adcx32(right,i,1)*exp(xhelp/adcx32(right,i,2))
814              xkorr = atten(right,32,i,xhelp)              xkorr = atten(right,32,i,xhelp)
815              xkorr=xkorr/hepratio                          xkorr=xkorr/hepratio
816              tof32(right,i,iadc) = xkorr/cos(theta)              tof32(right,i,iadc) = xkorr/cos(theta)
817              adcflag(ch32b(i),hb32b(i)) = 1              adcflag(ch32b(i),hb32b(i)) = 1
818           endif           endif
# Line 824  C--------------------------------------- Line 826  C---------------------------------------
826        do i=1,3        do i=1,3
827           xtofpos(i)=100.           xtofpos(i)=100.
828           ytofpos(i)=100.           ytofpos(i)=100.
829        enddo        enddo
830  C-----------------------------S1 --------------------------------  C-----------------------------S1 --------------------------------
831          
832        IF (tof11_i.GT.none_find) THEN        IF (tof11_i.GT.none_find) THEN
833           IF ((tof11(1,tof11_i,itdc).NE.4095).AND.           IF ((tof11(1,tof11_i,itdc).NE.4095).AND.
834       &       (tof11(2,tof11_i,itdc).NE.4095)) THEN       &       (tof11(2,tof11_i,itdc).NE.4095)) THEN
# Line 843  C-----------------------------S1 ------- Line 845  C-----------------------------S1 -------
845       +        -x_coor_lin12(tof12_i,offset))/x_coor_lin12(tof12_i,slope)       +        -x_coor_lin12(tof12_i,offset))/x_coor_lin12(tof12_i,slope)
846          if (abs(xtofpos(1)).gt.31.)  xtofpos(1)=101.          if (abs(xtofpos(1)).gt.31.)  xtofpos(1)=101.
847        endif        endif
848        endif              endif
849          
850  C-----------------------------S2 --------------------------------  C-----------------------------S2 --------------------------------
851          
852        IF (tof21_i.GT.none_find) THEN        IF (tof21_i.GT.none_find) THEN
853           IF ((tof21(1,tof21_i,itdc).NE.4095).AND.           IF ((tof21(1,tof21_i,itdc).NE.4095).AND.
854       &       (tof21(2,tof21_i,itdc).NE.4095)) THEN       &       (tof21(2,tof21_i,itdc).NE.4095)) THEN
855           xtofpos(2) = ((tof21(1,tof21_i,itdc)-tof21(2,tof21_i,itdc))/2.           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)       +        -x_coor_lin21(tof21_i,offset))/x_coor_lin21(tof21_i,slope)
857          if (abs(xtofpos(2)).gt.19.)  xtofpos(2)=101.          if (abs(xtofpos(2)).gt.19.)  xtofpos(2)=101.
858        endif        endif
859        endif        endif
860    
# Line 863  C-----------------------------S2 ------- Line 865  C-----------------------------S2 -------
865       +        -y_coor_lin22(tof22_i,offset))/y_coor_lin22(tof22_i,slope)       +        -y_coor_lin22(tof22_i,offset))/y_coor_lin22(tof22_i,slope)
866          if (abs(ytofpos(2)).gt.18.)  ytofpos(2)=101.          if (abs(ytofpos(2)).gt.18.)  ytofpos(2)=101.
867        endif        endif
868        endif              endif
869          
870  C-----------------------------S3 --------------------------------  C-----------------------------S3 --------------------------------
871          
872        IF (tof31_i.GT.none_find) THEN        IF (tof31_i.GT.none_find) THEN
873           IF ((tof31(1,tof31_i,itdc).NE.4095).AND.           IF ((tof31(1,tof31_i,itdc).NE.4095).AND.
874       &       (tof31(2,tof31_i,itdc).NE.4095)) THEN       &       (tof31(2,tof31_i,itdc).NE.4095)) THEN
# Line 883  C-----------------------------S3 ------- Line 885  C-----------------------------S3 -------
885       +        -x_coor_lin32(tof32_i,offset))/x_coor_lin32(tof32_i,slope)       +        -x_coor_lin32(tof32_i,offset))/x_coor_lin32(tof32_i,slope)
886          if (abs(xtofpos(3)).gt.19.)  xtofpos(3)=101.          if (abs(xtofpos(3)).gt.19.)  xtofpos(3)=101.
887        endif        endif
888        endif        endif
889          
890  c      do i=1,3  c      do i=1,3
891  c         if (abs(xtofpos(i)).gt.100.) then  c         if (abs(xtofpos(i)).gt.100.) then
892  c            xtofpos(i)=101.  c            xtofpos(i)=101.
# Line 892  c         endif Line 894  c         endif
894  c         if (abs(ytofpos(i)).gt.100.) then  c         if (abs(ytofpos(i)).gt.100.) then
895  c            ytofpos(i)=101.  c            ytofpos(i)=101.
896  c         endif  c         endif
897  c      enddo  c      enddo
898          
899    
900    
901    
# Line 906  C--------------------------------------- Line 908  C---------------------------------------
908           xhelp= 0.           xhelp= 0.
909           xhelp_a = tof11(left,i,iadc)           xhelp_a = tof11(left,i,iadc)
910           xhelp_t = tof11(left,i,itdc)           xhelp_t = tof11(left,i,itdc)
911  c          if (xhelp_a .eq.0) write (*,*) 'trk 11 ',i,xhelp_a            if (xhelp_a .eq.0) write (*,*) 'trk 11 ',i,xhelp_a
912    
913           if(xhelp_a<3786) xhelp = tw11(left,i)/sqrt(xhelp_a)           if(xhelp_a<3786) xhelp = tw11(left,i)/sqrt(xhelp_a)
914           tof11(left,i,itdc) = xhelp_t  + xhelp           tof11(left,i,itdc) = xhelp_t  + xhelp
# Line 922  c          if (xhelp_a .eq.0) write (*,* Line 924  c          if (xhelp_a .eq.0) write (*,*
924           xhelp= 0.           xhelp= 0.
925           xhelp_a = tof12(left,i,iadc)           xhelp_a = tof12(left,i,iadc)
926           xhelp_t = tof12(left,i,itdc)           xhelp_t = tof12(left,i,itdc)
927  c          if (xhelp_a .eq.0) write (*,*) 'trk 12 ',i,xhelp_a            if (xhelp_a .eq.0) write (*,*) 'trk 12 ',i,xhelp_a
928           if(xhelp_a<3786) xhelp = tw12(left,i)/sqrt(xhelp_a)           if(xhelp_a<3786) xhelp = tw12(left,i)/sqrt(xhelp_a)
929           tof12(left,i,itdc) = xhelp_t  + xhelp           tof12(left,i,itdc) = xhelp_t  + xhelp
930           tdc_c(ch12a(i),hb12a(i))=tof12(left,i,itdc)           tdc_c(ch12a(i),hb12a(i))=tof12(left,i,itdc)
# Line 932  c          if (xhelp_a .eq.0) write (*,* Line 934  c          if (xhelp_a .eq.0) write (*,*
934           tof12(right,i,itdc) = xhelp_t  + xhelp           tof12(right,i,itdc) = xhelp_t  + xhelp
935           tdc_c(ch12b(i),hb12b(i))=tof12(right,i,itdc)           tdc_c(ch12b(i),hb12b(i))=tof12(right,i,itdc)
936        ENDDO        ENDDO
937  C----  C----
938        DO i=1,2        DO i=1,2
939           xhelp= 0.           xhelp= 0.
940           xhelp_a = tof21(left,i,iadc)           xhelp_a = tof21(left,i,iadc)
941           xhelp_t = tof21(left,i,itdc)           xhelp_t = tof21(left,i,itdc)
942  c          if (xhelp_a .eq.0) write (*,*) ' trk 21 ',i,xhelp_a            if (xhelp_a .eq.0) write (*,*) ' trk 21 ',i,xhelp_a
943    
944           if(xhelp_a<3786) xhelp = tw21(left,i)/sqrt(xhelp_a)           if(xhelp_a<3786) xhelp = tw21(left,i)/sqrt(xhelp_a)
945           tof21(left,i,itdc) = xhelp_t  + xhelp           tof21(left,i,itdc) = xhelp_t  + xhelp
# Line 953  c          if (xhelp_a .eq.0) write (*,* Line 955  c          if (xhelp_a .eq.0) write (*,*
955           xhelp= 0.           xhelp= 0.
956           xhelp_a = tof22(left,i,iadc)           xhelp_a = tof22(left,i,iadc)
957           xhelp_t = tof22(left,i,itdc)           xhelp_t = tof22(left,i,itdc)
958  c          if (xhelp_a .eq.0) write (*,*) ' trk 22 ',i,xhelp_a            if (xhelp_a .eq.0) write (*,*) ' trk 22 ',i,xhelp_a
959           if(xhelp_a<3786) xhelp = tw22(left,i)/sqrt(xhelp_a)           if(xhelp_a<3786) xhelp = tw22(left,i)/sqrt(xhelp_a)
960           tof22(left,i,itdc) = xhelp_t  + xhelp           tof22(left,i,itdc) = xhelp_t  + xhelp
961           tdc_c(ch22a(i),hb22a(i))=tof22(left,i,itdc)           tdc_c(ch22a(i),hb22a(i))=tof22(left,i,itdc)
# Line 963  c          if (xhelp_a .eq.0) write (*,* Line 965  c          if (xhelp_a .eq.0) write (*,*
965           tof22(right,i,itdc) = xhelp_t  + xhelp           tof22(right,i,itdc) = xhelp_t  + xhelp
966           tdc_c(ch22b(i),hb22b(i))=tof22(right,i,itdc)           tdc_c(ch22b(i),hb22b(i))=tof22(right,i,itdc)
967        ENDDO        ENDDO
968  C----  C----
969    
970        DO i=1,3        DO i=1,3
971           xhelp= 0.           xhelp= 0.
972           xhelp_a = tof31(left,i,iadc)           xhelp_a = tof31(left,i,iadc)
973           xhelp_t = tof31(left,i,itdc)           xhelp_t = tof31(left,i,itdc)
974  c          if (xhelp_a .eq.0) write (*,*) ' trk 31 ',i,xhelp_a            if (xhelp_a .eq.0) write (*,*) ' trk 31 ',i,xhelp_a
975    
976           if(xhelp_a<3786) xhelp = tw31(left,i)/sqrt(xhelp_a)           if(xhelp_a<3786) xhelp = tw31(left,i)/sqrt(xhelp_a)
977           tof31(left,i,itdc) = xhelp_t  + xhelp           tof31(left,i,itdc) = xhelp_t  + xhelp
# Line 985  c          if (xhelp_a .eq.0) write (*,* Line 987  c          if (xhelp_a .eq.0) write (*,*
987           xhelp= 0.           xhelp= 0.
988           xhelp_a = tof32(left,i,iadc)           xhelp_a = tof32(left,i,iadc)
989           xhelp_t = tof32(left,i,itdc)           xhelp_t = tof32(left,i,itdc)
990  c          if (xhelp_a .eq.0) write (*,*) ' trk 32 ',i,xhelp_a            if (xhelp_a .eq.0) write (*,*) ' trk 32 ',i,xhelp_a
991    
992           if(xhelp_a<3786) xhelp = tw32(left,i)/sqrt(xhelp_a)           if(xhelp_a<3786) xhelp = tw32(left,i)/sqrt(xhelp_a)
993           tof32(left,i,itdc) = xhelp_t  + xhelp           tof32(left,i,itdc) = xhelp_t  + xhelp
# Line 1012  C-----------------------------S1 ------- Line 1014  C-----------------------------S1 -------
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            
1018           if (tof11(left,i,iadc).lt.3786) then           if (tof11(left,i,iadc).lt.3786) then
1019  c         if (adc(ch11a(i),hb11a(i)).lt.4095) then  c         if (adc(ch11a(i),hb11a(i)).lt.4095) then
1020              tof11(left,i,iadc) = tof11(left,i,iadc)*cos(theta)              tof11(left,i,iadc) = tof11(left,i,iadc)*cos(theta)
1021  c            xkorr=adcx11(left,i,1)*exp(-yhelp/adcx11(left,i,2))  c            xkorr=adcx11(left,i,1)*exp(-yhelp/adcx11(left,i,2))
1022              xkorr = atten(left,11,i,yhelp)              xkorr = atten(left,11,i,yhelp)
1023              xkorr=xkorr/hepratio              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            
1028           if (tof11(right,i,iadc).lt.3786) then           if (tof11(right,i,iadc).lt.3786) then
1029  c         if (adc(ch11b(i),hb11b(i)).lt.4095) then  c         if (adc(ch11b(i),hb11b(i)).lt.4095) then
1030              tof11(right,i,iadc) = tof11(right,i,iadc)*cos(theta)              tof11(right,i,iadc) = tof11(right,i,iadc)*cos(theta)
1031  c            xkorr=adcx11(right,i,1)*exp(yhelp/adcx11(right,i,2))  c            xkorr=adcx11(right,i,1)*exp(yhelp/adcx11(right,i,2))
1032              xkorr = atten(right,11,i,yhelp)              xkorr = atten(right,11,i,yhelp)
1033              xkorr=xkorr/hepratio                          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
# Line 1046  c         if (adc(ch12a(i),hb12a(i)).lt. Line 1048  c         if (adc(ch12a(i),hb12a(i)).lt.
1048              tof12(left,i,iadc) = tof12(left,i,iadc)*cos(theta)              tof12(left,i,iadc) = tof12(left,i,iadc)*cos(theta)
1049  c            xkorr=adcx12(left,i,1)*exp(-xhelp/adcx12(left,i,2))  c            xkorr=adcx12(left,i,1)*exp(-xhelp/adcx12(left,i,2))
1050              xkorr = atten(left,12,i,xhelp)              xkorr = atten(left,12,i,xhelp)
1051              xkorr=xkorr/hepratio                          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    
# Line 1055  c         if (adc(ch12b(i),hb12b(i)).lt. Line 1057  c         if (adc(ch12b(i),hb12b(i)).lt.
1057              tof12(right,i,iadc) = tof12(right,i,iadc)*cos(theta)              tof12(right,i,iadc) = tof12(right,i,iadc)*cos(theta)
1058  c            xkorr=adcx12(right,i,1)*exp(xhelp/adcx12(right,i,2))  c            xkorr=adcx12(right,i,1)*exp(xhelp/adcx12(right,i,2))
1059              xkorr = atten(right,12,i,xhelp)              xkorr = atten(right,12,i,xhelp)
1060              xkorr=xkorr/hepratio                          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 1074  c         if (adc(ch21a(i),hb21a(i)).lt. Line 1076  c         if (adc(ch21a(i),hb21a(i)).lt.
1076              tof21(left,i,iadc) = tof21(left,i,iadc)*cos(theta)              tof21(left,i,iadc) = tof21(left,i,iadc)*cos(theta)
1077  c            xkorr=adcx21(left,i,1)*exp(-xhelp/adcx21(left,i,2))  c            xkorr=adcx21(left,i,1)*exp(-xhelp/adcx21(left,i,2))
1078              xkorr = atten(left,21,i,xhelp)              xkorr = atten(left,21,i,xhelp)
1079              xkorr=xkorr/hepratio                          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    
# Line 1083  c         if (adc(ch21b(i),hb21b(i)).lt. Line 1085  c         if (adc(ch21b(i),hb21b(i)).lt.
1085              tof21(right,i,iadc) = tof21(right,i,iadc)*cos(theta)              tof21(right,i,iadc) = tof21(right,i,iadc)*cos(theta)
1086  c            xkorr=adcx21(right,i,1)*exp(xhelp/adcx21(right,i,2))  c            xkorr=adcx21(right,i,1)*exp(xhelp/adcx21(right,i,2))
1087              xkorr = atten(right,21,i,xhelp)              xkorr = atten(right,21,i,xhelp)
1088              xkorr=xkorr/hepratio                          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
# Line 1101  c         if (adc(ch22a(i),hb22a(i)).lt. Line 1103  c         if (adc(ch22a(i),hb22a(i)).lt.
1103              tof22(left,i,iadc) = tof22(left,i,iadc)*cos(theta)              tof22(left,i,iadc) = tof22(left,i,iadc)*cos(theta)
1104  c            xkorr=adcx22(left,i,1)*exp(-yhelp/adcx22(left,i,2))  c            xkorr=adcx22(left,i,1)*exp(-yhelp/adcx22(left,i,2))
1105              xkorr = atten(left,22,i,yhelp)              xkorr = atten(left,22,i,yhelp)
1106              xkorr=xkorr/hepratio                          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    
# Line 1110  c         if (adc(ch22b(i),hb22b(i)).lt. Line 1112  c         if (adc(ch22b(i),hb22b(i)).lt.
1112              tof22(right,i,iadc) = tof22(right,i,iadc)*cos(theta)              tof22(right,i,iadc) = tof22(right,i,iadc)*cos(theta)
1113  c            xkorr=adcx22(right,i,1)*exp(yhelp/adcx22(right,i,2))  c            xkorr=adcx22(right,i,1)*exp(yhelp/adcx22(right,i,2))
1114              xkorr = atten(right,22,i,yhelp)              xkorr = atten(right,22,i,yhelp)
1115              xkorr=xkorr/hepratio                          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 1121  C-----------------------------S3 ------- Line 1123  C-----------------------------S3 -------
1123        phi   = atan(tan(THYOUT(5))/tan(THXOUT(5)))        phi   = atan(tan(THYOUT(5))/tan(THXOUT(5)))
1124        theta = atan(tan(THXOUT(5))/cos(phi))        theta = atan(tan(THXOUT(5))/cos(phi))
1125  c      write(*,*) 'theta31 ',theta  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
# Line 1130  c         if (adc(ch31a(i),hb31a(i)).lt. Line 1132  c         if (adc(ch31a(i),hb31a(i)).lt.
1132              tof31(left,i,iadc) = tof31(left,i,iadc)*cos(theta)              tof31(left,i,iadc) = tof31(left,i,iadc)*cos(theta)
1133  c            xkorr=adcx31(left,i,1)*exp(-yhelp/adcx31(left,i,2))  c            xkorr=adcx31(left,i,1)*exp(-yhelp/adcx31(left,i,2))
1134              xkorr = atten(left,31,i,yhelp)              xkorr = atten(left,31,i,yhelp)
1135              xkorr=xkorr/hepratio                          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    
# Line 1139  c         if (adc(ch31b(i),hb31b(i)).lt. Line 1141  c         if (adc(ch31b(i),hb31b(i)).lt.
1141              tof31(right,i,iadc) = tof31(right,i,iadc)*cos(theta)              tof31(right,i,iadc) = tof31(right,i,iadc)*cos(theta)
1142  c            xkorr=adcx31(right,i,1)*exp(yhelp/adcx31(right,i,2))  c            xkorr=adcx31(right,i,1)*exp(yhelp/adcx31(right,i,2))
1143              xkorr = atten(right,31,i,yhelp)              xkorr = atten(right,31,i,yhelp)
1144              xkorr=xkorr/hepratio                          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
# Line 1157  c         if (adc(ch32a(i),hb32a(i)).lt. Line 1159  c         if (adc(ch32a(i),hb32a(i)).lt.
1159              tof32(left,i,iadc) = tof32(left,i,iadc)*cos(theta)              tof32(left,i,iadc) = tof32(left,i,iadc)*cos(theta)
1160  c            xkorr=adcx32(left,i,1)*exp(-xhelp/adcx32(left,i,2))  c            xkorr=adcx32(left,i,1)*exp(-xhelp/adcx32(left,i,2))
1161              xkorr = atten(left,32,i,xhelp)              xkorr = atten(left,32,i,xhelp)
1162              xkorr=xkorr/hepratio                          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    
# Line 1166  c         if (adc(ch32b(i),hb32b(i)).lt. Line 1168  c         if (adc(ch32b(i),hb32b(i)).lt.
1168              tof32(right,i,iadc) = tof32(right,i,iadc)*cos(theta)              tof32(right,i,iadc) = tof32(right,i,iadc)*cos(theta)
1169  c            xkorr=adcx32(right,i,1)*exp(xhelp/adcx32(right,i,2))  c            xkorr=adcx32(right,i,1)*exp(xhelp/adcx32(right,i,2))
1170              xkorr = atten(right,32,i,xhelp)              xkorr = atten(right,32,i,xhelp)
1171              xkorr=xkorr/hepratio                          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 1175  C--------------------------------------- Line 1177  C---------------------------------------
1177  C----------------------calculate Beta  ------------------------------  C----------------------calculate Beta  ------------------------------
1178  C-----------------------------------------------------------------------  C-----------------------------------------------------------------------
1179  C-------------------difference of sums  ---------------------------  C-------------------difference of sums  ---------------------------
1180  C      C
1181  C     DS = (t1+t2) - t3+t4)  C     DS = (t1+t2) - t3+t4)
1182  C     DS = c1 + c2/beta*cos(theta)  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:  C     instead of cos(theta) use factor F:
1187  C     F =  pathlength/d  C     F =  pathlength/d
1188  C     => beta = c2*F/(DS-c1))  C     => beta = c2*F/(DS-c1))
1189    
1190        dist = ZTOF(1) - ZTOF(5)        dist = ZTOF(1) - ZTOF(5)
# Line 1190  C     => beta = c2*F/(DS-c1)) Line 1192  C     => beta = c2*F/(DS-c1))
1192        DO I=1,5        DO I=1,5
1193           dl = dl + TLOUT(i)           dl = dl + TLOUT(i)
1194        ENDDO        ENDDO
1195        F = dl/dist        F = dl/dist
1196    
1197  C     S11 - S31  C     S11 - S31
1198  C      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
# Line 1200  C      IF (tof11_i.GT.none_find.AND.tof3 Line 1202  C      IF (tof11_i.GT.none_find.AND.tof3
1202           t2 = tof11(2,tof11_i,itdc)           t2 = tof11(2,tof11_i,itdc)
1203           t3 = tof31(1,tof31_i,itdc)           t3 = tof31(1,tof31_i,itdc)
1204           t4 = tof31(2,tof31_i,itdc)           t4 = tof31(2,tof31_i,itdc)
1205           IF ((t1.lt.4095).and.(t2.lt.4095).and.           IF ((t1.lt.4095).and.(t2.lt.4095).and.
1206       &        (t3.lt.4095).and.(t4.lt.4095)) THEN       &        (t3.lt.4095).and.(t4.lt.4095)) THEN
1207              xhelp1 = tof11(1,tof11_i,itdc)+tof11(2,tof11_i,itdc)              xhelp1 = tof11(1,tof11_i,itdc)+tof11(2,tof11_i,itdc)
1208              xhelp2 = tof31(1,tof31_i,itdc)+tof31(2,tof31_i,itdc)              xhelp2 = tof31(1,tof31_i,itdc)+tof31(2,tof31_i,itdc)
1209              ds = xhelp1-xhelp2              ds = xhelp1-xhelp2
1210              ihelp=(tof11_i-1)*3+tof31_i              ihelp=(tof11_i-1)*3+tof31_i
1211              c1 = k_S11S31(1,ihelp)              c1 = k_S11S31(1,ihelp)
1212              c2 = k_S11S31(2,ihelp)                      c2 = k_S11S31(2,ihelp)
1213              beta_a(1) = c2*F/(ds-c1)              beta_a(1) = c2*F/(ds-c1)
1214  C     write(*,*) 'S11-S31 ',xhelp1,xhelp2, beta_a(1)         write(*,*) 'S11-S31 ',c1,c2,F
1215           write(*,*) 'S11-S31 ',xhelp1,xhelp2, beta_a(1)
1216  C-------ToF Mask - S11 - S31  C-------ToF Mask - S11 - S31
1217    
1218              tofmask(ch11a(tof11_i),hb11a(tof11_i)) =              tofmask(ch11a(tof11_i),hb11a(tof11_i)) =
# Line 1219  C-------ToF Mask - S11 - S31 Line 1222  C-------ToF Mask - S11 - S31
1222    
1223              tofmask(ch31a(tof31_i),hb31a(tof31_i)) =              tofmask(ch31a(tof31_i),hb31a(tof31_i)) =
1224       $           tofmask(ch31a(tof31_i),hb31a(tof31_i)) + 1       $           tofmask(ch31a(tof31_i),hb31a(tof31_i)) + 1
1225              tofmask(ch31b(tof31_i),hb31b(tof31_i)) =              tofmask(ch31b(tof31_i),hb31b(tof31_i)) =
1226       $           tofmask(ch31b(tof31_i),hb31b(tof31_i)) + 1       $           tofmask(ch31b(tof31_i),hb31b(tof31_i)) + 1
1227    
1228           ENDIF           ENDIF
# Line 1230  C-------ToF Mask - S11 - S31 Line 1233  C-------ToF Mask - S11 - S31
1233        DO I=1,6        DO I=1,6
1234           dl = dl + TLOUT(i)           dl = dl + TLOUT(i)
1235        ENDDO        ENDDO
1236        F = dl/dist        F = dl/dist
1237          
1238  C     S11 - S32  C     S11 - S32
1239  C      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         IF ((tof11_i.GT.none_find).AND.(tof32_i.GT.none_find).AND.         IF ((tof11_i.GT.none_find).AND.(tof32_i.GT.none_find).AND.
# Line 1240  C      IF (tof11_i.GT.none_find.AND.tof3 Line 1243  C      IF (tof11_i.GT.none_find.AND.tof3
1243           t2 = tof11(2,tof11_i,itdc)           t2 = tof11(2,tof11_i,itdc)
1244           t3 = tof32(1,tof32_i,itdc)           t3 = tof32(1,tof32_i,itdc)
1245           t4 = tof32(2,tof32_i,itdc)           t4 = tof32(2,tof32_i,itdc)
1246           IF ((t1.lt.4095).and.(t2.lt.4095).and.           IF ((t1.lt.4095).and.(t2.lt.4095).and.
1247       &        (t3.lt.4095).and.(t4.lt.4095)) THEN       &        (t3.lt.4095).and.(t4.lt.4095)) THEN
1248              xhelp1 = tof11(1,tof11_i,itdc)+tof11(2,tof11_i,itdc)              xhelp1 = tof11(1,tof11_i,itdc)+tof11(2,tof11_i,itdc)
1249              xhelp2 = tof32(1,tof32_i,itdc)+tof32(2,tof32_i,itdc)              xhelp2 = tof32(1,tof32_i,itdc)+tof32(2,tof32_i,itdc)
1250              ds = xhelp1-xhelp2              ds = xhelp1-xhelp2
1251              ihelp=(tof11_i-1)*3+tof32_i              ihelp=(tof11_i-1)*3+tof32_i
1252              c1 = k_S11S32(1,ihelp)              c1 = k_S11S32(1,ihelp)
1253              c2 = k_S11S32(2,ihelp)                      c2 = k_S11S32(2,ihelp)
1254              beta_a(2) = c2*F/(ds-c1)              beta_a(2) = c2*F/(ds-c1)
1255  C     write(*,*) 'S11-S32 ',xhelp1,xhelp2, beta_a(2)  C     write(*,*) 'S11-S32 ',xhelp1,xhelp2, beta_a(2)
1256    
# Line 1260  C-------ToF Mask - S11 - S32 Line 1263  C-------ToF Mask - S11 - S32
1263    
1264              tofmask(ch32a(tof32_i),hb32a(tof32_i)) =              tofmask(ch32a(tof32_i),hb32a(tof32_i)) =
1265       $           tofmask(ch32a(tof32_i),hb32a(tof32_i)) + 1       $           tofmask(ch32a(tof32_i),hb32a(tof32_i)) + 1
1266              tofmask(ch32b(tof32_i),hb32b(tof32_i)) =              tofmask(ch32b(tof32_i),hb32b(tof32_i)) =
1267       $           tofmask(ch32b(tof32_i),hb32b(tof32_i)) + 1       $           tofmask(ch32b(tof32_i),hb32b(tof32_i)) + 1
1268    
1269  C-------  C-------
# Line 1274  C     S12 - S31 Line 1277  C     S12 - S31
1277        DO I=2,5        DO I=2,5
1278           dl = dl + TLOUT(i)           dl = dl + TLOUT(i)
1279        ENDDO        ENDDO
1280        F = dl/dist        F = dl/dist
1281    
1282  C      IF (tof12_i.GT.none_find.AND.tof31_i.GT.none_find) THEN  C      IF (tof12_i.GT.none_find.AND.tof31_i.GT.none_find) THEN
1283         IF ((tof12_i.GT.none_find).AND.(tof31_i.GT.none_find).AND.         IF ((tof12_i.GT.none_find).AND.(tof31_i.GT.none_find).AND.
# Line 1283  C      IF (tof12_i.GT.none_find.AND.tof3 Line 1286  C      IF (tof12_i.GT.none_find.AND.tof3
1286           t2 = tof12(2,tof12_i,itdc)           t2 = tof12(2,tof12_i,itdc)
1287           t3 = tof31(1,tof31_i,itdc)           t3 = tof31(1,tof31_i,itdc)
1288           t4 = tof31(2,tof31_i,itdc)           t4 = tof31(2,tof31_i,itdc)
1289           IF ((t1.lt.4095).and.(t2.lt.4095).and.           IF ((t1.lt.4095).and.(t2.lt.4095).and.
1290       &        (t3.lt.4095).and.(t4.lt.4095)) THEN       &        (t3.lt.4095).and.(t4.lt.4095)) THEN
1291              xhelp1 = tof12(1,tof12_i,itdc)+tof12(2,tof12_i,itdc)              xhelp1 = tof12(1,tof12_i,itdc)+tof12(2,tof12_i,itdc)
1292              xhelp2 = tof31(1,tof31_i,itdc)+tof31(2,tof31_i,itdc)              xhelp2 = tof31(1,tof31_i,itdc)+tof31(2,tof31_i,itdc)
1293              ds = xhelp1-xhelp2              ds = xhelp1-xhelp2
# Line 1303  C-------ToF Mask - S12 - S31 Line 1306  C-------ToF Mask - S12 - S31
1306    
1307              tofmask(ch31a(tof31_i),hb31a(tof31_i)) =              tofmask(ch31a(tof31_i),hb31a(tof31_i)) =
1308       $           tofmask(ch31a(tof31_i),hb31a(tof31_i)) + 1       $           tofmask(ch31a(tof31_i),hb31a(tof31_i)) + 1
1309              tofmask(ch31b(tof31_i),hb31b(tof31_i)) =              tofmask(ch31b(tof31_i),hb31b(tof31_i)) =
1310       $           tofmask(ch31b(tof31_i),hb31b(tof31_i)) + 1       $           tofmask(ch31b(tof31_i),hb31b(tof31_i)) + 1
1311    
1312  C-------  C-------
# Line 1318  C     S12 - S32 Line 1321  C     S12 - S32
1321        DO I=2,6        DO I=2,6
1322           dl = dl + TLOUT(i)           dl = dl + TLOUT(i)
1323        ENDDO        ENDDO
1324        F = dl/dist        F = dl/dist
1325    
1326  C      IF (tof12_i.GT.none_find.AND.tof32_i.GT.none_find) THEN  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.         IF ((tof12_i.GT.none_find).AND.(tof32_i.GT.none_find).AND.
# Line 1327  C      IF (tof12_i.GT.none_find.AND.tof3 Line 1330  C      IF (tof12_i.GT.none_find.AND.tof3
1330           t2 = tof12(2,tof12_i,itdc)           t2 = tof12(2,tof12_i,itdc)
1331           t3 = tof32(1,tof32_i,itdc)           t3 = tof32(1,tof32_i,itdc)
1332           t4 = tof32(2,tof32_i,itdc)           t4 = tof32(2,tof32_i,itdc)
1333           IF ((t1.lt.4095).and.(t2.lt.4095).and.           IF ((t1.lt.4095).and.(t2.lt.4095).and.
1334       &        (t3.lt.4095).and.(t4.lt.4095)) THEN       &        (t3.lt.4095).and.(t4.lt.4095)) THEN
1335              xhelp1 = tof12(1,tof12_i,itdc)+tof12(2,tof12_i,itdc)              xhelp1 = tof12(1,tof12_i,itdc)+tof12(2,tof12_i,itdc)
1336              xhelp2 = tof32(1,tof32_i,itdc)+tof32(2,tof32_i,itdc)              xhelp2 = tof32(1,tof32_i,itdc)+tof32(2,tof32_i,itdc)
1337              ds = xhelp1-xhelp2              ds = xhelp1-xhelp2
# Line 1347  C-------ToF Mask - S12 - S32 Line 1350  C-------ToF Mask - S12 - S32
1350    
1351              tofmask(ch32a(tof32_i),hb32a(tof32_i)) =              tofmask(ch32a(tof32_i),hb32a(tof32_i)) =
1352       $           tofmask(ch32a(tof32_i),hb32a(tof32_i)) + 1       $           tofmask(ch32a(tof32_i),hb32a(tof32_i)) + 1
1353              tofmask(ch32b(tof32_i),hb32b(tof32_i)) =              tofmask(ch32b(tof32_i),hb32b(tof32_i)) =
1354       $           tofmask(ch32b(tof32_i),hb32b(tof32_i)) + 1       $           tofmask(ch32b(tof32_i),hb32b(tof32_i)) + 1
1355    
1356  C-------  C-------
# Line 1362  C     S21 - S31 Line 1365  C     S21 - S31
1365        DO I=3,5        DO I=3,5
1366           dl = dl + TLOUT(i)           dl = dl + TLOUT(i)
1367        ENDDO        ENDDO
1368        F = dl/dist        F = dl/dist
1369    
1370  C      IF (tof21_i.GT.none_find.AND.tof31_i.GT.none_find) THEN  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.         IF ((tof21_i.GT.none_find).AND.(tof31_i.GT.none_find).AND.
# Line 1371  C      IF (tof21_i.GT.none_find.AND.tof3 Line 1374  C      IF (tof21_i.GT.none_find.AND.tof3
1374           t2 = tof21(2,tof21_i,itdc)           t2 = tof21(2,tof21_i,itdc)
1375           t3 = tof31(1,tof31_i,itdc)           t3 = tof31(1,tof31_i,itdc)
1376           t4 = tof31(2,tof31_i,itdc)           t4 = tof31(2,tof31_i,itdc)
1377           IF ((t1.lt.4095).and.(t2.lt.4095).and.           IF ((t1.lt.4095).and.(t2.lt.4095).and.
1378       &        (t3.lt.4095).and.(t4.lt.4095)) THEN       &        (t3.lt.4095).and.(t4.lt.4095)) THEN
1379              xhelp1 = tof21(1,tof21_i,itdc)+tof21(2,tof21_i,itdc)              xhelp1 = tof21(1,tof21_i,itdc)+tof21(2,tof21_i,itdc)
1380              xhelp2 = tof31(1,tof31_i,itdc)+tof31(2,tof31_i,itdc)              xhelp2 = tof31(1,tof31_i,itdc)+tof31(2,tof31_i,itdc)
1381              ds = xhelp1-xhelp2              ds = xhelp1-xhelp2
# Line 1390  C-------ToF Mask - S21 - S31 Line 1393  C-------ToF Mask - S21 - S31
1393    
1394              tofmask(ch31a(tof31_i),hb31a(tof31_i)) =              tofmask(ch31a(tof31_i),hb31a(tof31_i)) =
1395       $           tofmask(ch31a(tof31_i),hb31a(tof31_i)) + 1       $           tofmask(ch31a(tof31_i),hb31a(tof31_i)) + 1
1396              tofmask(ch31b(tof31_i),hb31b(tof31_i)) =              tofmask(ch31b(tof31_i),hb31b(tof31_i)) =
1397       $           tofmask(ch31b(tof31_i),hb31b(tof31_i)) + 1       $           tofmask(ch31b(tof31_i),hb31b(tof31_i)) + 1
1398    
1399  C-------  C-------
# Line 1405  C     S21 - S32 Line 1408  C     S21 - S32
1408        DO I=3,6        DO I=3,6
1409           dl = dl + TLOUT(i)           dl = dl + TLOUT(i)
1410        ENDDO        ENDDO
1411        F = dl/dist        F = dl/dist
1412    
1413  C      IF (tof21_i.GT.none_find.AND.tof32_i.GT.none_find) THEN  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.         IF ((tof21_i.GT.none_find).AND.(tof32_i.GT.none_find).AND.
# Line 1414  C      IF (tof21_i.GT.none_find.AND.tof3 Line 1417  C      IF (tof21_i.GT.none_find.AND.tof3
1417           t2 = tof21(2,tof21_i,itdc)           t2 = tof21(2,tof21_i,itdc)
1418           t3 = tof32(1,tof32_i,itdc)           t3 = tof32(1,tof32_i,itdc)
1419           t4 = tof32(2,tof32_i,itdc)           t4 = tof32(2,tof32_i,itdc)
1420           IF ((t1.lt.4095).and.(t2.lt.4095).and.           IF ((t1.lt.4095).and.(t2.lt.4095).and.
1421       &        (t3.lt.4095).and.(t4.lt.4095)) THEN       &        (t3.lt.4095).and.(t4.lt.4095)) THEN
1422              xhelp1 = tof21(1,tof21_i,itdc)+tof21(2,tof21_i,itdc)              xhelp1 = tof21(1,tof21_i,itdc)+tof21(2,tof21_i,itdc)
1423              xhelp2 = tof32(1,tof32_i,itdc)+tof32(2,tof32_i,itdc)              xhelp2 = tof32(1,tof32_i,itdc)+tof32(2,tof32_i,itdc)
1424              ds = xhelp1-xhelp2              ds = xhelp1-xhelp2
# Line 1433  C-------ToF Mask - S21 - S32 Line 1436  C-------ToF Mask - S21 - S32
1436    
1437              tofmask(ch32a(tof32_i),hb32a(tof32_i)) =              tofmask(ch32a(tof32_i),hb32a(tof32_i)) =
1438       $           tofmask(ch32a(tof32_i),hb32a(tof32_i)) + 1       $           tofmask(ch32a(tof32_i),hb32a(tof32_i)) + 1
1439              tofmask(ch32b(tof32_i),hb32b(tof32_i)) =              tofmask(ch32b(tof32_i),hb32b(tof32_i)) =
1440       $           tofmask(ch32b(tof32_i),hb32b(tof32_i)) + 1       $           tofmask(ch32b(tof32_i),hb32b(tof32_i)) + 1
1441    
1442  C-------  C-------
# Line 1448  C     S22 - S31 Line 1451  C     S22 - S31
1451        DO I=4,5        DO I=4,5
1452           dl = dl + TLOUT(i)           dl = dl + TLOUT(i)
1453        ENDDO        ENDDO
1454        F = dl/dist        F = dl/dist
1455    
1456  C      IF (tof22_i.GT.none_find.AND.tof31_i.GT.none_find) THEN  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.         IF ((tof22_i.GT.none_find).AND.(tof31_i.GT.none_find).AND.
# Line 1457  C      IF (tof22_i.GT.none_find.AND.tof3 Line 1460  C      IF (tof22_i.GT.none_find.AND.tof3
1460           t2 = tof22(2,tof22_i,itdc)           t2 = tof22(2,tof22_i,itdc)
1461           t3 = tof31(1,tof31_i,itdc)           t3 = tof31(1,tof31_i,itdc)
1462           t4 = tof31(2,tof31_i,itdc)           t4 = tof31(2,tof31_i,itdc)
1463           IF ((t1.lt.4095).and.(t2.lt.4095).and.           IF ((t1.lt.4095).and.(t2.lt.4095).and.
1464       &        (t3.lt.4095).and.(t4.lt.4095)) THEN       &        (t3.lt.4095).and.(t4.lt.4095)) THEN
1465              xhelp1 = tof22(1,tof22_i,itdc)+tof22(2,tof22_i,itdc)              xhelp1 = tof22(1,tof22_i,itdc)+tof22(2,tof22_i,itdc)
1466              xhelp2 = tof31(1,tof31_i,itdc)+tof31(2,tof31_i,itdc)              xhelp2 = tof31(1,tof31_i,itdc)+tof31(2,tof31_i,itdc)
1467              ds = xhelp1-xhelp2              ds = xhelp1-xhelp2
# Line 1476  C-------ToF Mask - S22 - S31 Line 1479  C-------ToF Mask - S22 - S31
1479    
1480              tofmask(ch31a(tof31_i),hb31a(tof31_i)) =              tofmask(ch31a(tof31_i),hb31a(tof31_i)) =
1481       $           tofmask(ch31a(tof31_i),hb31a(tof31_i)) + 1       $           tofmask(ch31a(tof31_i),hb31a(tof31_i)) + 1
1482              tofmask(ch31b(tof31_i),hb31b(tof31_i)) =              tofmask(ch31b(tof31_i),hb31b(tof31_i)) =
1483       $           tofmask(ch31b(tof31_i),hb31b(tof31_i)) + 1       $           tofmask(ch31b(tof31_i),hb31b(tof31_i)) + 1
1484    
1485  C-------  C-------
1486    
1487           ENDIF           ENDIF
1488        ENDIF        ENDIF
1489          
1490  C     S22 - S32  C     S22 - S32
1491    
1492        dist = ZTOF(4) - ZTOF(6)        dist = ZTOF(4) - ZTOF(6)
# Line 1491  C     S22 - S32 Line 1494  C     S22 - S32
1494        DO I=4,6        DO I=4,6
1495           dl = dl + TLOUT(i)           dl = dl + TLOUT(i)
1496        ENDDO        ENDDO
1497        F = dl/dist        F = dl/dist
1498    
1499  C      IF (tof22_i.GT.none_find.AND.tof32_i.GT.none_find) THEN  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.         IF ((tof22_i.GT.none_find).AND.(tof32_i.GT.none_find).AND.
# Line 1500  C      IF (tof22_i.GT.none_find.AND.tof3 Line 1503  C      IF (tof22_i.GT.none_find.AND.tof3
1503           t2 = tof22(2,tof22_i,itdc)           t2 = tof22(2,tof22_i,itdc)
1504           t3 = tof32(1,tof32_i,itdc)           t3 = tof32(1,tof32_i,itdc)
1505           t4 = tof32(2,tof32_i,itdc)           t4 = tof32(2,tof32_i,itdc)
1506           IF ((t1.lt.4095).and.(t2.lt.4095).and.           IF ((t1.lt.4095).and.(t2.lt.4095).and.
1507       &        (t3.lt.4095).and.(t4.lt.4095)) THEN       &        (t3.lt.4095).and.(t4.lt.4095)) THEN
1508              xhelp1 = tof22(1,tof22_i,itdc)+tof22(2,tof22_i,itdc)              xhelp1 = tof22(1,tof22_i,itdc)+tof22(2,tof22_i,itdc)
1509              xhelp2 = tof32(1,tof32_i,itdc)+tof32(2,tof32_i,itdc)              xhelp2 = tof32(1,tof32_i,itdc)+tof32(2,tof32_i,itdc)
1510              ds = xhelp1-xhelp2              ds = xhelp1-xhelp2
# Line 1519  C-------ToF Mask - S22 - S32 Line 1522  C-------ToF Mask - S22 - S32
1522    
1523              tofmask(ch32a(tof32_i),hb32a(tof32_i)) =              tofmask(ch32a(tof32_i),hb32a(tof32_i)) =
1524       $           tofmask(ch32a(tof32_i),hb32a(tof32_i)) + 1       $           tofmask(ch32a(tof32_i),hb32a(tof32_i)) + 1
1525              tofmask(ch32b(tof32_i),hb32b(tof32_i)) =              tofmask(ch32b(tof32_i),hb32b(tof32_i)) =
1526       $           tofmask(ch32b(tof32_i),hb32b(tof32_i)) + 1       $           tofmask(ch32b(tof32_i),hb32b(tof32_i)) + 1
1527    
1528  C-------  C-------
# Line 1534  C     S11 - S21 Line 1537  C     S11 - S21
1537        DO I=1,3        DO I=1,3
1538           dl = dl + TLOUT(i)           dl = dl + TLOUT(i)
1539        ENDDO        ENDDO
1540        F = dl/dist        F = dl/dist
1541    
1542  C      IF (tof11_i.GT.none_find.AND.tof21_i.GT.none_find) THEN  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.         IF ((tof11_i.GT.none_find).AND.(tof21_i.GT.none_find).AND.
# Line 1543  C      IF (tof11_i.GT.none_find.AND.tof2 Line 1546  C      IF (tof11_i.GT.none_find.AND.tof2
1546           t2 = tof11(2,tof11_i,itdc)           t2 = tof11(2,tof11_i,itdc)
1547           t3 = tof21(1,tof21_i,itdc)           t3 = tof21(1,tof21_i,itdc)
1548           t4 = tof21(2,tof21_i,itdc)           t4 = tof21(2,tof21_i,itdc)
1549           IF ((t1.lt.4095).and.(t2.lt.4095).and.           IF ((t1.lt.4095).and.(t2.lt.4095).and.
1550       &        (t3.lt.4095).and.(t4.lt.4095)) THEN       &        (t3.lt.4095).and.(t4.lt.4095)) THEN
1551              xhelp1 = tof11(1,tof11_i,itdc)+tof11(2,tof11_i,itdc)              xhelp1 = tof11(1,tof11_i,itdc)+tof11(2,tof11_i,itdc)
1552              xhelp2 = tof21(1,tof21_i,itdc)+tof21(2,tof21_i,itdc)              xhelp2 = tof21(1,tof21_i,itdc)+tof21(2,tof21_i,itdc)
1553              ds = xhelp1-xhelp2              ds = xhelp1-xhelp2
# Line 1562  C-------ToF Mask - S11 - S21 Line 1565  C-------ToF Mask - S11 - S21
1565    
1566              tofmask(ch21a(tof21_i),hb21a(tof21_i)) =              tofmask(ch21a(tof21_i),hb21a(tof21_i)) =
1567       $           tofmask(ch21a(tof21_i),hb21a(tof21_i)) + 1       $           tofmask(ch21a(tof21_i),hb21a(tof21_i)) + 1
1568              tofmask(ch21b(tof21_i),hb21b(tof21_i)) =              tofmask(ch21b(tof21_i),hb21b(tof21_i)) =
1569       $           tofmask(ch21b(tof21_i),hb21b(tof21_i)) + 1       $           tofmask(ch21b(tof21_i),hb21b(tof21_i)) + 1
1570    
1571  C-------  C-------
1572    
1573           ENDIF           ENDIF
1574        ENDIF        ENDIF
1575          
1576  C     S11 - S22  C     S11 - S22
1577    
1578        dist = ZTOF(1) - ZTOF(4)        dist = ZTOF(1) - ZTOF(4)
# Line 1577  C     S11 - S22 Line 1580  C     S11 - S22
1580        DO I=1,4        DO I=1,4
1581           dl = dl + TLOUT(i)           dl = dl + TLOUT(i)
1582        ENDDO        ENDDO
1583        F = dl/dist        F = dl/dist
1584    
1585  C      IF (tof11_i.GT.none_find.AND.tof22_i.GT.none_find) THEN  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.         IF ((tof11_i.GT.none_find).AND.(tof22_i.GT.none_find).AND.
# Line 1586  C      IF (tof11_i.GT.none_find.AND.tof2 Line 1589  C      IF (tof11_i.GT.none_find.AND.tof2
1589           t2 = tof11(2,tof11_i,itdc)           t2 = tof11(2,tof11_i,itdc)
1590           t3 = tof22(1,tof22_i,itdc)           t3 = tof22(1,tof22_i,itdc)
1591           t4 = tof22(2,tof22_i,itdc)           t4 = tof22(2,tof22_i,itdc)
1592           IF ((t1.lt.4095).and.(t2.lt.4095).and.           IF ((t1.lt.4095).and.(t2.lt.4095).and.
1593       &        (t3.lt.4095).and.(t4.lt.4095)) THEN       &        (t3.lt.4095).and.(t4.lt.4095)) THEN
1594              xhelp1 = tof11(1,tof11_i,itdc)+tof11(2,tof11_i,itdc)              xhelp1 = tof11(1,tof11_i,itdc)+tof11(2,tof11_i,itdc)
1595              xhelp2 = tof22(1,tof22_i,itdc)+tof22(2,tof22_i,itdc)              xhelp2 = tof22(1,tof22_i,itdc)+tof22(2,tof22_i,itdc)
1596              ds = xhelp1-xhelp2              ds = xhelp1-xhelp2
# Line 1605  C-------ToF Mask - S11 - S22 Line 1608  C-------ToF Mask - S11 - S22
1608    
1609              tofmask(ch22a(tof22_i),hb22a(tof22_i)) =              tofmask(ch22a(tof22_i),hb22a(tof22_i)) =
1610       $           tofmask(ch22a(tof22_i),hb22a(tof22_i)) + 1       $           tofmask(ch22a(tof22_i),hb22a(tof22_i)) + 1
1611              tofmask(ch22b(tof22_i),hb22b(tof22_i)) =              tofmask(ch22b(tof22_i),hb22b(tof22_i)) =
1612       $           tofmask(ch22b(tof22_i),hb22b(tof22_i)) + 1       $           tofmask(ch22b(tof22_i),hb22b(tof22_i)) + 1
1613    
1614  C-------  C-------
# Line 1620  C     S12 - S21 Line 1623  C     S12 - S21
1623        DO I=2,3        DO I=2,3
1624           dl = dl + TLOUT(i)           dl = dl + TLOUT(i)
1625        ENDDO        ENDDO
1626        F = dl/dist        F = dl/dist
1627    
1628  C      IF (tof12_i.GT.none_find.AND.tof21_i.GT.none_find) THEN  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.         IF ((tof12_i.GT.none_find).AND.(tof21_i.GT.none_find).AND.
# Line 1629  C      IF (tof12_i.GT.none_find.AND.tof2 Line 1632  C      IF (tof12_i.GT.none_find.AND.tof2
1632           t2 = tof12(2,tof12_i,itdc)           t2 = tof12(2,tof12_i,itdc)
1633           t3 = tof21(1,tof21_i,itdc)           t3 = tof21(1,tof21_i,itdc)
1634           t4 = tof21(2,tof21_i,itdc)           t4 = tof21(2,tof21_i,itdc)
1635           IF ((t1.lt.4095).and.(t2.lt.4095).and.           IF ((t1.lt.4095).and.(t2.lt.4095).and.
1636       &        (t3.lt.4095).and.(t4.lt.4095)) THEN       &        (t3.lt.4095).and.(t4.lt.4095)) THEN
1637              xhelp1 = tof12(1,tof12_i,itdc)+tof12(2,tof12_i,itdc)              xhelp1 = tof12(1,tof12_i,itdc)+tof12(2,tof12_i,itdc)
1638              xhelp2 = tof21(1,tof21_i,itdc)+tof21(2,tof21_i,itdc)              xhelp2 = tof21(1,tof21_i,itdc)+tof21(2,tof21_i,itdc)
1639              ds = xhelp1-xhelp2              ds = xhelp1-xhelp2
# Line 1648  C-------ToF Mask - S12 - S21 Line 1651  C-------ToF Mask - S12 - S21
1651    
1652              tofmask(ch21a(tof21_i),hb21a(tof21_i)) =              tofmask(ch21a(tof21_i),hb21a(tof21_i)) =
1653       $           tofmask(ch21a(tof21_i),hb21a(tof21_i)) + 1       $           tofmask(ch21a(tof21_i),hb21a(tof21_i)) + 1
1654              tofmask(ch21b(tof21_i),hb21b(tof21_i)) =              tofmask(ch21b(tof21_i),hb21b(tof21_i)) =
1655       $           tofmask(ch21b(tof21_i),hb21b(tof21_i)) + 1       $           tofmask(ch21b(tof21_i),hb21b(tof21_i)) + 1
1656    
1657  C-------  C-------
# Line 1663  C     S12 - S22 Line 1666  C     S12 - S22
1666        DO I=2,4        DO I=2,4
1667           dl = dl + TLOUT(i)           dl = dl + TLOUT(i)
1668        ENDDO        ENDDO
1669        F = dl/dist        F = dl/dist
1670    
1671  C      IF (tof12_i.GT.none_find.AND.tof22_i.GT.none_find) THEN  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.         IF ((tof12_i.GT.none_find).AND.(tof22_i.GT.none_find).AND.
# Line 1672  C      IF (tof12_i.GT.none_find.AND.tof2 Line 1675  C      IF (tof12_i.GT.none_find.AND.tof2
1675           t2 = tof12(2,tof12_i,itdc)           t2 = tof12(2,tof12_i,itdc)
1676           t3 = tof22(1,tof22_i,itdc)           t3 = tof22(1,tof22_i,itdc)
1677           t4 = tof22(2,tof22_i,itdc)           t4 = tof22(2,tof22_i,itdc)
1678           IF ((t1.lt.4095).and.(t2.lt.4095).and.           IF ((t1.lt.4095).and.(t2.lt.4095).and.
1679       &        (t3.lt.4095).and.(t4.lt.4095)) THEN       &        (t3.lt.4095).and.(t4.lt.4095)) THEN
1680              xhelp1 = tof12(1,tof12_i,itdc)+tof12(2,tof12_i,itdc)              xhelp1 = tof12(1,tof12_i,itdc)+tof12(2,tof12_i,itdc)
1681              xhelp2 = tof22(1,tof22_i,itdc)+tof22(2,tof22_i,itdc)              xhelp2 = tof22(1,tof22_i,itdc)+tof22(2,tof22_i,itdc)
1682              ds = xhelp1-xhelp2              ds = xhelp1-xhelp2
# Line 1691  C-------ToF Mask - S12 - S22 Line 1694  C-------ToF Mask - S12 - S22
1694    
1695              tofmask(ch22a(tof22_i),hb22a(tof22_i)) =              tofmask(ch22a(tof22_i),hb22a(tof22_i)) =
1696       $           tofmask(ch22a(tof22_i),hb22a(tof22_i)) + 1       $           tofmask(ch22a(tof22_i),hb22a(tof22_i)) + 1
1697              tofmask(ch22b(tof22_i),hb22b(tof22_i)) =              tofmask(ch22b(tof22_i),hb22b(tof22_i)) =
1698       $           tofmask(ch22b(tof22_i),hb22b(tof22_i)) + 1       $           tofmask(ch22b(tof22_i),hb22b(tof22_i)) + 1
1699    
1700  C-------  C-------
1701    
1702           ENDIF           ENDIF
1703        ENDIF        ENDIF
1704          
1705  C-------  C-------
1706    
1707        icount=0        icount=0
# Line 1711  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)  c       IF (tof11_i.GT.none_find)
1728  c     &  write(*,*) '11 ',tof11(1,tof11_i,itdc),tof11(2,tof11_i,itdc)  c     &  write(*,*) '11 ',tof11(1,tof11_i,itdc),tof11(2,tof11_i,itdc)
1729  c       IF (tof12_i.GT.none_find)  c       IF (tof12_i.GT.none_find)
1730  c     &  write(*,*) '12 ',tof12(1,tof12_i,itdc),tof12(2,tof12_i,itdc)  c     &  write(*,*) '12 ',tof12(1,tof12_i,itdc),tof12(2,tof12_i,itdc)
1731    
1732  c       IF (tof21_i.GT.none_find)  c       IF (tof21_i.GT.none_find)
1733  c     &  write(*,*) '21 ',tof21(1,tof21_i,itdc),tof21(2,tof21_i,itdc)  c     &  write(*,*) '21 ',tof21(1,tof21_i,itdc),tof21(2,tof21_i,itdc)
1734  c       IF (tof22_i.GT.none_find)  c       IF (tof22_i.GT.none_find)
1735  c     &  write(*,*) '22 ',tof22(1,tof22_i,itdc),tof22(2,tof22_i,itdc)  c     &  write(*,*) '22 ',tof22(1,tof22_i,itdc),tof22(2,tof22_i,itdc)
1736    
1737  c       IF (tof31_i.GT.none_find)  c       IF (tof31_i.GT.none_find)
1738  c     &  write(*,*) '31 ',tof31(1,tof31_i,itdc),tof31(2,tof31_i,itdc)  c     &  write(*,*) '31 ',tof31(1,tof31_i,itdc),tof31(2,tof31_i,itdc)
1739  c       IF (tof32_i.GT.none_find)  c       IF (tof32_i.GT.none_find)
1740  c     &  write(*,*) '32 ',tof32(1,tof32_i,itdc),tof32(2,tof32_i,itdc)  c     &  write(*,*) '32 ',tof32(1,tof32_i,itdc),tof32(2,tof32_i,itdc)
1741    
1742  c      write(*,*) xtofpos  c      write(*,*) xtofpos
# Line 1756  c      write(*,*) ytr_tof Line 1759  c      write(*,*) ytr_tof
1759    
1760    
1761  C------------------------------------------------------------------  C------------------------------------------------------------------
1762  C------------------------------------------------------------------  
Legend:
Removed from v.1.11  
changed lines
  Added in v.1.12
  ViewVC Help
Powered by ViewVC 1.1.23