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

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

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

revision 1.1 by mocchiut, Sat Jun 17 12:14:56 2006 UTC revision 1.3 by mocchiut, Thu Aug 10 06:32:04 2006 UTC
# Line 92  C     Line 92  C    
92        enddo        enddo
93    
94    
95          do i=1,12
96             do j=1,4
97                tofmask(j,i) = 0
98             enddo
99          enddo
100    
101    
102  c the calibration files are read in the main program from xxx_tofcalib.rz  c the calibration files are read in the main program from xxx_tofcalib.rz
103    
104    
# Line 611  c       DATA tof32_y/ -5.0,0.0,5.0/ Line 618  c       DATA tof32_y/ -5.0,0.0,5.0/
618           i = tof11_i           i = tof11_i
619           xdummy=tof11(left,i,iadc)           xdummy=tof11(left,i,iadc)
620           tof11(left,i,iadc) = tof11(left,i,iadc)*cos(theta13)           tof11(left,i,iadc) = tof11(left,i,iadc)*cos(theta13)
621           if (tof11(left,i,iadc).lt.1000) then           if (tof11(left,i,iadc).lt.4095) then
622              xkorr=adcx11(left,i,1)*exp(-yhelp/adcx11(left,i,2))              xkorr=adcx11(left,i,1)*exp(-yhelp/adcx11(left,i,2))
623              xkorr0=adcx11(left,i,1)              xkorr0=adcx11(left,i,1)
624              adctof_c(ch11a(i),hb11a(i))=tof11(left,i,iadc)/xkorr              adctof_c(ch11a(i),hb11a(i))=tof11(left,i,iadc)/xkorr
625           endif           endif
626    
627           tof11(right,i,iadc) = tof11(right,i,iadc)*cos(theta13)           tof11(right,i,iadc) = tof11(right,i,iadc)*cos(theta13)
628           if (tof11(right,i,iadc).lt.1000) then           if (tof11(right,i,iadc).lt.4095) then
629              xkorr=adcx11(right,i,1)*exp(yhelp/adcx11(right,i,2))              xkorr=adcx11(right,i,1)*exp(yhelp/adcx11(right,i,2))
630              xkorr0=adcx11(right,i,1)              xkorr0=adcx11(right,i,1)
631              adctof_c(ch11b(i),hb11b(i))=tof11(right,i,iadc)/xkorr              adctof_c(ch11b(i),hb11b(i))=tof11(right,i,iadc)/xkorr
# Line 633  c       DATA tof32_y/ -5.0,0.0,5.0/ Line 640  c       DATA tof32_y/ -5.0,0.0,5.0/
640    
641           i = tof12_i           i = tof12_i
642           tof12(left,i,iadc) = tof12(left,i,iadc)*cos(theta13)           tof12(left,i,iadc) = tof12(left,i,iadc)*cos(theta13)
643           if (tof12(left,i,iadc).lt.1000) then           if (tof12(left,i,iadc).lt.4095) then
644              xkorr=adcx12(left,i,1)*exp(-xhelp/adcx12(left,i,2))              xkorr=adcx12(left,i,1)*exp(-xhelp/adcx12(left,i,2))
645              xkorr0=adcx12(left,i,1)              xkorr0=adcx12(left,i,1)
646              adctof_c(ch12a(i),hb12a(i))=tof12(left,i,iadc)/xkorr              adctof_c(ch12a(i),hb12a(i))=tof12(left,i,iadc)/xkorr
647           endif           endif
648    
649           tof12(right,i,iadc) = tof12(right,i,iadc)*cos(theta13)           tof12(right,i,iadc) = tof12(right,i,iadc)*cos(theta13)
650           if (tof12(right,i,iadc).lt.1000) then           if (tof12(right,i,iadc).lt.4095) then
651              xkorr=adcx12(right,i,1)*exp(xhelp/adcx12(right,i,2))              xkorr=adcx12(right,i,1)*exp(xhelp/adcx12(right,i,2))
652              xkorr0=adcx12(right,i,1)              xkorr0=adcx12(right,i,1)
653              adctof_c(ch12b(i),hb12b(i))=tof12(right,i,iadc)/xkorr              adctof_c(ch12b(i),hb12b(i))=tof12(right,i,iadc)/xkorr
# Line 657  C-----------------------------S2 ------- Line 664  C-----------------------------S2 -------
664    
665           i = tof21_i           i = tof21_i
666           tof21(left,i,iadc) = tof21(left,i,iadc)*cos(theta13)           tof21(left,i,iadc) = tof21(left,i,iadc)*cos(theta13)
667           if (tof21(left,i,iadc).lt.1000) then           if (tof21(left,i,iadc).lt.4095) then
668              xkorr=adcx21(left,i,1)*exp(-xhelp/adcx21(left,i,2))              xkorr=adcx21(left,i,1)*exp(-xhelp/adcx21(left,i,2))
669              xkorr0=adcx21(left,i,1)              xkorr0=adcx21(left,i,1)
670              adctof_c(ch21a(i),hb21a(i))=tof21(left,i,iadc)/xkorr              adctof_c(ch21a(i),hb21a(i))=tof21(left,i,iadc)/xkorr
671           endif           endif
672    
673           tof21(right,i,iadc) = tof21(right,i,iadc)*cos(theta13)           tof21(right,i,iadc) = tof21(right,i,iadc)*cos(theta13)
674           if (tof21(right,i,iadc).lt.1000) then           if (tof21(right,i,iadc).lt.4095) then
675              xkorr=adcx21(right,i,1)*exp(xhelp/adcx21(right,i,2))              xkorr=adcx21(right,i,1)*exp(xhelp/adcx21(right,i,2))
676              xkorr0=adcx21(right,i,1)              xkorr0=adcx21(right,i,1)
677              adctof_c(ch21b(i),hb21b(i))=tof21(right,i,iadc)/xkorr              adctof_c(ch21b(i),hb21b(i))=tof21(right,i,iadc)/xkorr
# Line 680  C-----------------------------S2 ------- Line 687  C-----------------------------S2 -------
687    
688           i = tof22_i           i = tof22_i
689           tof22(left,i,iadc) = tof22(left,i,iadc)*cos(theta13)           tof22(left,i,iadc) = tof22(left,i,iadc)*cos(theta13)
690           if (tof22(left,i,iadc).lt.1000) then           if (tof22(left,i,iadc).lt.4095) then
691              xkorr=adcx22(left,i,1)*exp(-yhelp/adcx22(left,i,2))              xkorr=adcx22(left,i,1)*exp(-yhelp/adcx22(left,i,2))
692              xkorr0=adcx22(left,i,1)              xkorr0=adcx22(left,i,1)
693              adctof_c(ch22a(i),hb22a(i))=tof22(left,i,iadc)/xkorr              adctof_c(ch22a(i),hb22a(i))=tof22(left,i,iadc)/xkorr
694           endif           endif
695    
696           tof22(right,i,iadc) = tof22(right,i,iadc)*cos(theta13)           tof22(right,i,iadc) = tof22(right,i,iadc)*cos(theta13)
697           if (tof22(right,i,iadc).lt.1000) then           if (tof22(right,i,iadc).lt.4095) then
698              xkorr=adcx22(right,i,1)*exp(yhelp/adcx22(right,i,2))              xkorr=adcx22(right,i,1)*exp(yhelp/adcx22(right,i,2))
699              xkorr0=adcx22(right,i,1)              xkorr0=adcx22(right,i,1)
700              adctof_c(ch22b(i),hb22b(i))=tof22(right,i,iadc)/xkorr              adctof_c(ch22b(i),hb22b(i))=tof22(right,i,iadc)/xkorr
# Line 704  C-----------------------------S3 ------- Line 711  C-----------------------------S3 -------
711    
712           i = tof31_i           i = tof31_i
713           tof31(left,i,iadc) = tof31(left,i,iadc)*cos(theta13)           tof31(left,i,iadc) = tof31(left,i,iadc)*cos(theta13)
714           if (tof31(left,i,iadc).lt.1000) then           if (tof31(left,i,iadc).lt.4095) then
715              xkorr=adcx31(left,i,1)*exp(-yhelp/adcx31(left,i,2))              xkorr=adcx31(left,i,1)*exp(-yhelp/adcx31(left,i,2))
716              xkorr0=adcx31(left,i,1)              xkorr0=adcx31(left,i,1)
717              adctof_c(ch31a(i),hb31a(i))=tof31(left,i,iadc)/xkorr              adctof_c(ch31a(i),hb31a(i))=tof31(left,i,iadc)/xkorr
718           endif           endif
719    
720           tof31(right,i,iadc) = tof31(right,i,iadc)*cos(theta13)           tof31(right,i,iadc) = tof31(right,i,iadc)*cos(theta13)
721           if (tof31(right,i,iadc).lt.1000) then           if (tof31(right,i,iadc).lt.4095) then
722              xkorr=adcx31(right,i,1)*exp(yhelp/adcx31(right,i,2))              xkorr=adcx31(right,i,1)*exp(yhelp/adcx31(right,i,2))
723              xkorr0=adcx31(right,i,1)              xkorr0=adcx31(right,i,1)
724              adctof_c(ch31b(i),hb31b(i))=tof31(right,i,iadc)/xkorr              adctof_c(ch31b(i),hb31b(i))=tof31(right,i,iadc)/xkorr
# Line 726  C-----------------------------S3 ------- Line 733  C-----------------------------S3 -------
733    
734           i = tof32_i           i = tof32_i
735           tof32(left,i,iadc) = tof32(left,i,iadc)*cos(theta13)           tof32(left,i,iadc) = tof32(left,i,iadc)*cos(theta13)
736           if (tof32(left,i,iadc).lt.1000) then           if (tof32(left,i,iadc).lt.4095) then
737              xkorr=adcx32(left,i,1)*exp(-xhelp/adcx32(left,i,2))              xkorr=adcx32(left,i,1)*exp(-xhelp/adcx32(left,i,2))
738              xkorr0=adcx32(left,i,1)              xkorr0=adcx32(left,i,1)
739              adctof_c(ch32a(i),hb32a(i))=tof32(left,i,iadc)/xkorr              adctof_c(ch32a(i),hb32a(i))=tof32(left,i,iadc)/xkorr
740           endif           endif
741    
742           tof32(right,i,iadc) = tof32(right,i,iadc)*cos(theta13)           tof32(right,i,iadc) = tof32(right,i,iadc)*cos(theta13)
743           if (tof32(right,i,iadc).lt.1000) then           if (tof32(right,i,iadc).lt.4095) then
744              xkorr=adcx32(right,i,1)*exp(xhelp/adcx32(right,i,2))              xkorr=adcx32(right,i,1)*exp(xhelp/adcx32(right,i,2))
745              xkorr0=adcx32(right,i,1)              xkorr0=adcx32(right,i,1)
746              adctof_c(ch32b(i),hb32b(i))=tof32(right,i,iadc)/xkorr              adctof_c(ch32b(i),hb32b(i))=tof32(right,i,iadc)/xkorr
# Line 760  C     S11 - S31 Line 767  C     S11 - S31
767           c1 = k_S11S31(1,ihelp)           c1 = k_S11S31(1,ihelp)
768           c2 = k_S11S31(2,ihelp)           c2 = k_S11S31(2,ihelp)
769           betatof_a(1) = c2/(cos(theta13)*(ds-c1))           betatof_a(1) = c2/(cos(theta13)*(ds-c1))
770    
771    C------- ToF Mask - S11 - S31
772    
773             tofmask(ch11a(tof11_i),hb11a(tof11_i)) =
774         $        tofmask(ch11a(tof11_i),hb11a(tof11_i)) + 1
775             tofmask(ch11b(tof11_i),hb11b(tof11_i)) =
776         $        tofmask(ch11b(tof11_i),hb11b(tof11_i)) + 1
777    
778             tofmask(ch31a(tof31_i),hb31a(tof31_i)) =
779         $        tofmask(ch31a(tof31_i),hb31a(tof31_i)) + 1
780             tofmask(ch31b(tof31_i),hb31b(tof31_i)) =
781         $        tofmask(ch31b(tof31_i),hb31b(tof31_i)) + 1
782    
783    C-------
784          
785        ENDIF        ENDIF
786                
787  C     S11 - S32  C     S11 - S32
# Line 771  C     S11 - S32 Line 793  C     S11 - S32
793           c1 = k_S11S32(1,ihelp)           c1 = k_S11S32(1,ihelp)
794           c2 = k_S11S32(2,ihelp)           c2 = k_S11S32(2,ihelp)
795           betatof_a(2) = c2/(cos(theta13)*(ds-c1))           betatof_a(2) = c2/(cos(theta13)*(ds-c1))
       ENDIF  
796    
797    C------- ToF Mask - S11 - S32
798    
799             tofmask(ch11a(tof11_i),hb11a(tof11_i)) =
800         $        tofmask(ch11a(tof11_i),hb11a(tof11_i)) + 1
801             tofmask(ch11b(tof11_i),hb11b(tof11_i)) =
802         $        tofmask(ch11b(tof11_i),hb11b(tof11_i)) + 1
803    
804             tofmask(ch32a(tof32_i),hb32a(tof32_i)) =
805         $        tofmask(ch32a(tof32_i),hb32a(tof32_i)) + 1
806             tofmask(ch32b(tof32_i),hb32b(tof32_i)) =
807         $        tofmask(ch32b(tof32_i),hb32b(tof32_i)) + 1
808    
809    C-------
810    
811          ENDIF
812          
813  C     S12 - S31  C     S12 - S31
814        IF (tof12_i.GT.none_find.AND.tof31_i.GT.none_find) THEN        IF (tof12_i.GT.none_find.AND.tof31_i.GT.none_find) THEN
815           xhelp1 = tof12(1,tof12_i,itdc)+tof12(2,tof12_i,itdc)           xhelp1 = tof12(1,tof12_i,itdc)+tof12(2,tof12_i,itdc)
# Line 782  C     S12 - S31 Line 819  C     S12 - S31
819           c1 = k_S12S31(1,ihelp)           c1 = k_S12S31(1,ihelp)
820           c2 = k_S12S31(2,ihelp)           c2 = k_S12S31(2,ihelp)
821           betatof_a(3) = c2/(cos(theta13)*(ds-c1))           betatof_a(3) = c2/(cos(theta13)*(ds-c1))
       ENDIF  
822    
823    C------- ToF Mask - S12 - S31
824    
825             tofmask(ch12a(tof12_i),hb12a(tof12_i)) =
826         $        tofmask(ch12a(tof12_i),hb12a(tof12_i)) + 1
827             tofmask(ch12b(tof12_i),hb12b(tof12_i)) =
828         $        tofmask(ch12b(tof12_i),hb12b(tof12_i)) + 1
829    
830             tofmask(ch31a(tof31_i),hb31a(tof31_i)) =
831         $        tofmask(ch31a(tof31_i),hb31a(tof31_i)) + 1
832             tofmask(ch31b(tof31_i),hb31b(tof31_i)) =
833         $        tofmask(ch31b(tof31_i),hb31b(tof31_i)) + 1
834    
835    C-------
836    
837          ENDIF
838                
839  C     S12 - S32  C     S12 - S32
840        IF (tof12_i.GT.none_find.AND.tof32_i.GT.none_find) THEN        IF (tof12_i.GT.none_find.AND.tof32_i.GT.none_find) THEN
841           xhelp1 = tof12(1,tof12_i,itdc)+tof12(2,tof12_i,itdc)           xhelp1 = tof12(1,tof12_i,itdc)+tof12(2,tof12_i,itdc)
# Line 793  C     S12 - S32 Line 845  C     S12 - S32
845           c1 = k_S12S32(1,ihelp)           c1 = k_S12S32(1,ihelp)
846           c2 = k_S12S32(2,ihelp)           c2 = k_S12S32(2,ihelp)
847           betatof_a(4) = c2/(cos(theta13)*(ds-c1))           betatof_a(4) = c2/(cos(theta13)*(ds-c1))
848    
849    C------- ToF Mask - S12 - S32
850    
851             tofmask(ch12a(tof12_i),hb12a(tof12_i)) =
852         $        tofmask(ch12a(tof12_i),hb12a(tof12_i)) + 1
853             tofmask(ch12b(tof12_i),hb12b(tof12_i)) =
854         $        tofmask(ch12b(tof12_i),hb12b(tof12_i)) + 1
855    
856             tofmask(ch32a(tof32_i),hb32a(tof32_i)) =
857         $        tofmask(ch32a(tof32_i),hb32a(tof32_i)) + 1
858             tofmask(ch32b(tof32_i),hb32b(tof32_i)) =
859         $        tofmask(ch32b(tof32_i),hb32b(tof32_i)) + 1
860    
861    C-------
862    
863        ENDIF        ENDIF
864    
865  C     S21 - S31  C     S21 - S31
# Line 804  C     S21 - S31 Line 871  C     S21 - S31
871           c1 = k_S21S31(1,ihelp)           c1 = k_S21S31(1,ihelp)
872           c2 = k_S21S31(2,ihelp)           c2 = k_S21S31(2,ihelp)
873           betatof_a(5) = c2/(cos(theta23)*(ds-c1))           betatof_a(5) = c2/(cos(theta23)*(ds-c1))
874    
875    C------- ToF Mask - S21 - S31
876    
877             tofmask(ch21a(tof21_i),hb21a(tof21_i)) =
878         $        tofmask(ch21a(tof21_i),hb21a(tof21_i)) + 1
879             tofmask(ch21b(tof21_i),hb21b(tof21_i)) =
880         $        tofmask(ch21b(tof21_i),hb21b(tof21_i)) + 1
881    
882             tofmask(ch31a(tof31_i),hb31a(tof31_i)) =
883         $        tofmask(ch31a(tof31_i),hb31a(tof31_i)) + 1
884             tofmask(ch31b(tof31_i),hb31b(tof31_i)) =
885         $        tofmask(ch31b(tof31_i),hb31b(tof31_i)) + 1
886    
887    C-------
888    
889        ENDIF        ENDIF
890    
891  C     S21 - S32  C     S21 - S32
# Line 815  C     S21 - S32 Line 897  C     S21 - S32
897           c1 = k_S21S32(1,ihelp)           c1 = k_S21S32(1,ihelp)
898           c2 = k_S21S32(2,ihelp)           c2 = k_S21S32(2,ihelp)
899           betatof_a(6) = c2/(cos(theta23)*(ds-c1))           betatof_a(6) = c2/(cos(theta23)*(ds-c1))
900                      
901    C------- ToF Mask - S21 - S32
902    
903             tofmask(ch21a(tof21_i),hb21a(tof21_i)) =
904         $        tofmask(ch21a(tof21_i),hb21a(tof21_i)) + 1
905             tofmask(ch21b(tof21_i),hb21b(tof21_i)) =
906         $        tofmask(ch21b(tof21_i),hb21b(tof21_i)) + 1
907    
908             tofmask(ch32a(tof32_i),hb32a(tof32_i)) =
909         $        tofmask(ch32a(tof32_i),hb32a(tof32_i)) + 1
910             tofmask(ch32b(tof32_i),hb32b(tof32_i)) =
911         $        tofmask(ch32b(tof32_i),hb32b(tof32_i)) + 1
912    
913    C-------
914    
915        ENDIF        ENDIF
916    
917  C     S22 - S31  C     S22 - S31
# Line 826  C     S22 - S31 Line 923  C     S22 - S31
923           c1 = k_S22S31(1,ihelp)           c1 = k_S22S31(1,ihelp)
924           c2 = k_S22S31(2,ihelp)           c2 = k_S22S31(2,ihelp)
925           betatof_a(7) = c2/(cos(theta13)*(ds-c1))           betatof_a(7) = c2/(cos(theta13)*(ds-c1))
926    
927    C------- ToF Mask - S22 - S31
928    
929             tofmask(ch22a(tof22_i),hb22a(tof22_i)) =
930         $        tofmask(ch22a(tof22_i),hb22a(tof22_i)) + 1
931             tofmask(ch22b(tof22_i),hb22b(tof22_i)) =
932         $        tofmask(ch22b(tof22_i),hb22b(tof22_i)) + 1
933    
934             tofmask(ch31a(tof31_i),hb31a(tof31_i)) =
935         $        tofmask(ch31a(tof31_i),hb31a(tof31_i)) + 1
936             tofmask(ch31b(tof31_i),hb31b(tof31_i)) =
937         $        tofmask(ch31b(tof31_i),hb31b(tof31_i)) + 1
938    
939    C-------          
940    
941        ENDIF        ENDIF
942                
943  C     S22 - S32  C     S22 - S32
# Line 837  C     S22 - S32 Line 949  C     S22 - S32
949           c1 = k_S22S32(1,ihelp)           c1 = k_S22S32(1,ihelp)
950           c2 = k_S22S32(2,ihelp)           c2 = k_S22S32(2,ihelp)
951           betatof_a(8) = c2/(cos(theta13)*(ds-c1))           betatof_a(8) = c2/(cos(theta13)*(ds-c1))
952    
953    C------- ToF Mask - S22 - S32
954    
955             tofmask(ch22a(tof22_i),hb22a(tof22_i)) =
956         $        tofmask(ch22a(tof22_i),hb22a(tof22_i)) + 1
957             tofmask(ch22b(tof22_i),hb22b(tof22_i)) =
958         $        tofmask(ch22b(tof22_i),hb22b(tof22_i)) + 1
959    
960             tofmask(ch32a(tof32_i),hb32a(tof32_i)) =
961         $        tofmask(ch32a(tof32_i),hb32a(tof32_i)) + 1
962             tofmask(ch32b(tof32_i),hb32b(tof32_i)) =
963         $        tofmask(ch32b(tof32_i),hb32b(tof32_i)) + 1
964    
965    C-------  
966    
967        ENDIF        ENDIF
968    
969  C     S11 - S21  C     S11 - S21
# Line 848  C     S11 - S21 Line 975  C     S11 - S21
975           c1 = k_S11S21(1,ihelp)           c1 = k_S11S21(1,ihelp)
976           c2 = k_S11S21(2,ihelp)           c2 = k_S11S21(2,ihelp)
977           betatof_a(9) = c2/(cos(theta13)*(ds-c1))           betatof_a(9) = c2/(cos(theta13)*(ds-c1))
978    
979    C------- ToF Mask - S11 - S21
980    
981             tofmask(ch11a(tof11_i),hb11a(tof11_i)) =
982         $        tofmask(ch11a(tof11_i),hb11a(tof11_i)) + 1
983             tofmask(ch11b(tof11_i),hb11b(tof11_i)) =
984         $        tofmask(ch11b(tof11_i),hb11b(tof11_i)) + 1
985    
986             tofmask(ch21a(tof21_i),hb21a(tof21_i)) =
987         $        tofmask(ch21a(tof21_i),hb21a(tof21_i)) + 1
988             tofmask(ch21b(tof21_i),hb21b(tof21_i)) =
989         $        tofmask(ch21b(tof21_i),hb21b(tof21_i)) + 1
990    
991    C-------  
992    
993        ENDIF        ENDIF
994                
995  C     S11 - S22  C     S11 - S22
# Line 859  C     S11 - S22 Line 1001  C     S11 - S22
1001           c1 = k_S11S22(1,ihelp)           c1 = k_S11S22(1,ihelp)
1002           c2 = k_S11S22(2,ihelp)           c2 = k_S11S22(2,ihelp)
1003           betatof_a(10) = c2/(cos(theta13)*(ds-c1))           betatof_a(10) = c2/(cos(theta13)*(ds-c1))
1004    
1005    C------- ToF Mask - S11 - S22
1006    
1007             tofmask(ch11a(tof11_i),hb11a(tof11_i)) =
1008         $        tofmask(ch11a(tof11_i),hb11a(tof11_i)) + 1
1009             tofmask(ch11b(tof11_i),hb11b(tof11_i)) =
1010         $        tofmask(ch11b(tof11_i),hb11b(tof11_i)) + 1
1011    
1012             tofmask(ch22a(tof22_i),hb22a(tof22_i)) =
1013         $        tofmask(ch22a(tof22_i),hb22a(tof22_i)) + 1
1014             tofmask(ch22b(tof22_i),hb22b(tof22_i)) =
1015         $        tofmask(ch22b(tof22_i),hb22b(tof22_i)) + 1
1016    
1017    C-------  
1018    
1019        ENDIF        ENDIF
1020    
1021  C     S12 - S21  C     S12 - S21
# Line 870  C     S12 - S21 Line 1027  C     S12 - S21
1027           c1 = k_S12S21(1,ihelp)           c1 = k_S12S21(1,ihelp)
1028           c2 = k_S12S21(2,ihelp)           c2 = k_S12S21(2,ihelp)
1029           betatof_a(11) = c2/(cos(theta13)*(ds-c1))           betatof_a(11) = c2/(cos(theta13)*(ds-c1))
1030    
1031    C------- ToF Mask - S12 - S21
1032    
1033             tofmask(ch12a(tof12_i),hb12a(tof12_i)) =
1034         $        tofmask(ch12a(tof12_i),hb12a(tof12_i)) + 1
1035             tofmask(ch12b(tof12_i),hb12b(tof12_i)) =
1036         $        tofmask(ch12b(tof12_i),hb12b(tof12_i)) + 1
1037    
1038             tofmask(ch21a(tof21_i),hb21a(tof21_i)) =
1039         $        tofmask(ch21a(tof21_i),hb21a(tof21_i)) + 1
1040             tofmask(ch21b(tof21_i),hb21b(tof21_i)) =
1041         $        tofmask(ch21b(tof21_i),hb21b(tof21_i)) + 1
1042    
1043    C-------  
1044    
1045        ENDIF        ENDIF
1046    
1047  C     S12 - S22  C     S12 - S22
# Line 881  C     S12 - S22 Line 1053  C     S12 - S22
1053           c1 = k_S12S22(1,ihelp)           c1 = k_S12S22(1,ihelp)
1054           c2 = k_S12S22(2,ihelp)           c2 = k_S12S22(2,ihelp)
1055           betatof_a(12) = c2/(cos(theta13)*(ds-c1))           betatof_a(12) = c2/(cos(theta13)*(ds-c1))
1056    
1057    C------- ToF Mask - S12 - S22
1058    
1059             tofmask(ch12a(tof12_i),hb12a(tof12_i)) =
1060         $        tofmask(ch12a(tof12_i),hb12a(tof12_i)) + 1
1061             tofmask(ch12b(tof12_i),hb12b(tof12_i)) =
1062         $        tofmask(ch12b(tof12_i),hb12b(tof12_i)) + 1
1063    
1064             tofmask(ch22a(tof22_i),hb22a(tof22_i)) =
1065         $        tofmask(ch22a(tof22_i),hb22a(tof22_i)) + 1
1066             tofmask(ch22b(tof22_i),hb22b(tof22_i)) =
1067         $        tofmask(ch22b(tof22_i),hb22b(tof22_i)) + 1
1068    
1069    C-------  
1070    
1071        ENDIF        ENDIF
1072          
1073  C-------  C---------------------------------------------------------
1074    
1075        icount=0        icount=0
1076        sw=0.        sw=0.
# Line 903  C------- Line 1090  C-------
1090                
1091        if (icount.gt.0) beta_mean=sxw/sw        if (icount.gt.0) beta_mean=sxw/sw
1092        betatof_a(13) = beta_mean        betatof_a(13) = beta_mean
 c      write(*,*) '------------- end  tofl2com ----------'  
1093    
1094   100  continue   100  continue
1095    
Legend:
Removed from v.1.1  
changed lines
  Added in v.1.3
  ViewVC Help
Powered by ViewVC 1.1.23