/[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.9 by mocchiut, Mon Feb 5 15:36:44 2007 UTC revision 1.10 by mocchiut, Wed Feb 7 08:17:17 2007 UTC
# Line 24  C  jan-07 WM: in the xtofpos calculation Line 24  C  jan-07 WM: in the xtofpos calculation
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: claculation of zenith angles using DOTRACK2  C  jan-05 WM: bug fixed: calculation of zenith angles using DOTRACK2
28  C             was incorrect  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****************************************************************************  C****************************************************************************
32        IMPLICIT NONE        IMPLICIT NONE
33  C      C    
# Line 856  C--------------------------------------- Line 858  C---------------------------------------
858  C--------------------Time walk correction  -------------------------  C--------------------Time walk correction  -------------------------
859  C--------------------------------------------------------------------  C--------------------------------------------------------------------
860    
861    
862        DO i=1,8        DO i=1,8
863             xhelp= 0.
864           xhelp_a = tof11(left,i,iadc)           xhelp_a = tof11(left,i,iadc)
865           xhelp_t = tof11(left,i,itdc)           xhelp_t = tof11(left,i,itdc)
866           if(xhelp_a<4095) xhelp = tw11(left,i)/sqrt(xhelp_a)           if(xhelp_a<4095) xhelp = tw11(left,i)/sqrt(xhelp_a)
# Line 870  C--------------------------------------- Line 874  C---------------------------------------
874        ENDDO        ENDDO
875    
876        DO i=1,6        DO i=1,6
877             xhelp= 0.
878           xhelp_a = tof12(left,i,iadc)           xhelp_a = tof12(left,i,iadc)
879           xhelp_t = tof12(left,i,itdc)           xhelp_t = tof12(left,i,itdc)
880           if(xhelp_a<4095) xhelp = tw12(left,i)/sqrt(xhelp_a)           if(xhelp_a<4095) xhelp = tw12(left,i)/sqrt(xhelp_a)
# Line 883  C--------------------------------------- Line 888  C---------------------------------------
888        ENDDO        ENDDO
889  C----  C----
890        DO i=1,2        DO i=1,2
891             xhelp= 0.
892           xhelp_a = tof21(left,i,iadc)           xhelp_a = tof21(left,i,iadc)
893           xhelp_t = tof21(left,i,itdc)           xhelp_t = tof21(left,i,itdc)
894           if(xhelp_a<4095) xhelp = tw21(left,i)/sqrt(xhelp_a)           if(xhelp_a<4095) xhelp = tw21(left,i)/sqrt(xhelp_a)
# Line 896  C---- Line 902  C----
902        ENDDO        ENDDO
903    
904        DO i=1,2        DO i=1,2
905             xhelp= 0.
906           xhelp_a = tof22(left,i,iadc)           xhelp_a = tof22(left,i,iadc)
907           xhelp_t = tof22(left,i,itdc)           xhelp_t = tof22(left,i,itdc)
908           if(xhelp_a<4095) xhelp = tw22(left,i)/sqrt(xhelp_a)           if(xhelp_a<4095) xhelp = tw22(left,i)/sqrt(xhelp_a)
# Line 910  C---- Line 917  C----
917  C----  C----
918    
919        DO i=1,3        DO i=1,3
920             xhelp= 0.
921           xhelp_a = tof31(left,i,iadc)           xhelp_a = tof31(left,i,iadc)
922           xhelp_t = tof31(left,i,itdc)           xhelp_t = tof31(left,i,itdc)
923           if(xhelp_a<4095) xhelp = tw31(left,i)/sqrt(xhelp_a)           if(xhelp_a<4095) xhelp = tw31(left,i)/sqrt(xhelp_a)
# Line 923  C---- Line 931  C----
931        ENDDO        ENDDO
932    
933        DO i=1,3        DO i=1,3
934             xhelp= 0.
935           xhelp_a = tof32(left,i,iadc)           xhelp_a = tof32(left,i,iadc)
936           xhelp_t = tof32(left,i,itdc)           xhelp_t = tof32(left,i,itdc)
937           if(xhelp_a<4095) xhelp = tw32(left,i)/sqrt(xhelp_a)           if(xhelp_a<4095) xhelp = tw32(left,i)/sqrt(xhelp_a)
Legend:
Removed from v.1.9  
changed lines
  Added in v.1.10
  ViewVC Help
Powered by ViewVC 1.1.23