/[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.12 by pam-de, Mon Aug 20 14:22:33 2007 UTC revision 1.14 by mocchiut, Thu Aug 30 07:59:25 2007 UTC
# Line 284  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           PRINT *,' TOF - WARNING F77: track with R = 0, discarded'  c         PRINT *,' TOF - WARNING F77: track with R = 0, discarded'
288           GOTO 969           GOTO 969
289        ENDIF        ENDIF
290  *     -------- *** tracking routine *** --------  *     -------- *** tracking routine *** --------
# Line 295  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           print *,' TOF - WARNING F77: tracking failed '  c         print *,' TOF - WARNING F77: tracking failed '
299           goto 969           goto 969
300        endif        endif
301  *     ------------------------------------------  *     ------------------------------------------
# Line 908  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            if (xhelp_a .eq.0) write (*,*) 'trk 11 ',i,xhelp_a  c          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 924  C--------------------------------------- Line 924  C---------------------------------------
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            if (xhelp_a .eq.0) write (*,*) 'trk 12 ',i,xhelp_a  c          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 939  C---- Line 939  C----
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            if (xhelp_a .eq.0) write (*,*) ' trk 21 ',i,xhelp_a  c          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 955  C---- Line 955  C----
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            if (xhelp_a .eq.0) write (*,*) ' trk 22 ',i,xhelp_a  c          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 971  C---- Line 971  C----
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            if (xhelp_a .eq.0) write (*,*) ' trk 31 ',i,xhelp_a  c          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 987  C---- Line 987  C----
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            if (xhelp_a .eq.0) write (*,*) ' trk 32 ',i,xhelp_a  c          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 1211  C      IF (tof11_i.GT.none_find.AND.tof3 Line 1211  C      IF (tof11_i.GT.none_find.AND.tof3
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         write(*,*) 'S11-S31 ',c1,c2,F  c       write(*,*) 'S11-S31 ',c1,c2,F
1215         write(*,*) 'S11-S31 ',xhelp1,xhelp2, beta_a(1)  c       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)) =
Legend:
Removed from v.1.12  
changed lines
  Added in v.1.14
  ViewVC Help
Powered by ViewVC 1.1.23