/[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.13 by mocchiut, Thu Aug 30 07:18:32 2007 UTC
# 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.13
  ViewVC Help
Powered by ViewVC 1.1.23