--- DarthVader/ToFLevel2/src/tofl2com.for 2006/06/17 12:14:56 1.1 +++ DarthVader/ToFLevel2/src/tofl2com.for 2006/07/06 09:03:28 1.2 @@ -92,6 +92,13 @@ enddo + do i=1,12 + do j=1,4 + tofmask(j,i) = 0 + enddo + enddo + + c the calibration files are read in the main program from xxx_tofcalib.rz @@ -760,6 +767,21 @@ c1 = k_S11S31(1,ihelp) c2 = k_S11S31(2,ihelp) betatof_a(1) = c2/(cos(theta13)*(ds-c1)) + +C------- ToF Mask - S11 - S31 + + tofmask(ch11a(tof11_i),hb11a(tof11_i)) = + $ tofmask(ch11a(tof11_i),hb11a(tof11_i)) + 1 + tofmask(ch11b(tof11_i),hb11b(tof11_i)) = + $ tofmask(ch11b(tof11_i),hb11b(tof11_i)) + 1 + + tofmask(ch31a(tof31_i),hb31a(tof31_i)) = + $ tofmask(ch31a(tof31_i),hb31a(tof31_i)) + 1 + tofmask(ch31b(tof31_i),hb31b(tof31_i)) = + $ tofmask(ch31b(tof31_i),hb31b(tof31_i)) + 1 + +C------- + ENDIF C S11 - S32 @@ -771,8 +793,23 @@ c1 = k_S11S32(1,ihelp) c2 = k_S11S32(2,ihelp) betatof_a(2) = c2/(cos(theta13)*(ds-c1)) - ENDIF +C------- ToF Mask - S11 - S32 + + tofmask(ch11a(tof11_i),hb11a(tof11_i)) = + $ tofmask(ch11a(tof11_i),hb11a(tof11_i)) + 1 + tofmask(ch11b(tof11_i),hb11b(tof11_i)) = + $ tofmask(ch11b(tof11_i),hb11b(tof11_i)) + 1 + + tofmask(ch32a(tof32_i),hb32a(tof32_i)) = + $ tofmask(ch32a(tof32_i),hb32a(tof32_i)) + 1 + tofmask(ch32b(tof32_i),hb32b(tof32_i)) = + $ tofmask(ch32b(tof32_i),hb32b(tof32_i)) + 1 + +C------- + + ENDIF + C S12 - S31 IF (tof12_i.GT.none_find.AND.tof31_i.GT.none_find) THEN xhelp1 = tof12(1,tof12_i,itdc)+tof12(2,tof12_i,itdc) @@ -782,8 +819,23 @@ c1 = k_S12S31(1,ihelp) c2 = k_S12S31(2,ihelp) betatof_a(3) = c2/(cos(theta13)*(ds-c1)) - ENDIF +C------- ToF Mask - S12 - S31 + + tofmask(ch12a(tof12_i),hb12a(tof12_i)) = + $ tofmask(ch12a(tof12_i),hb12a(tof12_i)) + 1 + tofmask(ch12b(tof12_i),hb12b(tof12_i)) = + $ tofmask(ch12b(tof12_i),hb12b(tof12_i)) + 1 + + tofmask(ch31a(tof31_i),hb31a(tof31_i)) = + $ tofmask(ch31a(tof31_i),hb31a(tof31_i)) + 1 + tofmask(ch31b(tof31_i),hb31b(tof31_i)) = + $ tofmask(ch31b(tof31_i),hb31b(tof31_i)) + 1 + +C------- + + ENDIF + C S12 - S32 IF (tof12_i.GT.none_find.AND.tof32_i.GT.none_find) THEN xhelp1 = tof12(1,tof12_i,itdc)+tof12(2,tof12_i,itdc) @@ -793,6 +845,21 @@ c1 = k_S12S32(1,ihelp) c2 = k_S12S32(2,ihelp) betatof_a(4) = c2/(cos(theta13)*(ds-c1)) + +C------- ToF Mask - S12 - S32 + + tofmask(ch12a(tof12_i),hb12a(tof12_i)) = + $ tofmask(ch12a(tof12_i),hb12a(tof12_i)) + 1 + tofmask(ch12b(tof12_i),hb12b(tof12_i)) = + $ tofmask(ch12b(tof12_i),hb12b(tof12_i)) + 1 + + tofmask(ch32a(tof32_i),hb32a(tof32_i)) = + $ tofmask(ch32a(tof32_i),hb32a(tof32_i)) + 1 + tofmask(ch32b(tof32_i),hb32b(tof32_i)) = + $ tofmask(ch32b(tof32_i),hb32b(tof32_i)) + 1 + +C------- + ENDIF C S21 - S31 @@ -804,6 +871,21 @@ c1 = k_S21S31(1,ihelp) c2 = k_S21S31(2,ihelp) betatof_a(5) = c2/(cos(theta23)*(ds-c1)) + +C------- ToF Mask - S21 - S31 + + tofmask(ch21a(tof21_i),hb21a(tof21_i)) = + $ tofmask(ch21a(tof21_i),hb21a(tof21_i)) + 1 + tofmask(ch21b(tof21_i),hb21b(tof21_i)) = + $ tofmask(ch21b(tof21_i),hb21b(tof21_i)) + 1 + + tofmask(ch31a(tof31_i),hb31a(tof31_i)) = + $ tofmask(ch31a(tof31_i),hb31a(tof31_i)) + 1 + tofmask(ch31b(tof31_i),hb31b(tof31_i)) = + $ tofmask(ch31b(tof31_i),hb31b(tof31_i)) + 1 + +C------- + ENDIF C S21 - S32 @@ -815,6 +897,21 @@ c1 = k_S21S32(1,ihelp) c2 = k_S21S32(2,ihelp) betatof_a(6) = c2/(cos(theta23)*(ds-c1)) + +C------- ToF Mask - S21 - S32 + + tofmask(ch21a(tof21_i),hb21a(tof21_i)) = + $ tofmask(ch21a(tof21_i),hb21a(tof21_i)) + 1 + tofmask(ch21b(tof21_i),hb21b(tof21_i)) = + $ tofmask(ch21b(tof21_i),hb21b(tof21_i)) + 1 + + tofmask(ch32a(tof32_i),hb32a(tof32_i)) = + $ tofmask(ch32a(tof32_i),hb32a(tof32_i)) + 1 + tofmask(ch32b(tof32_i),hb32b(tof32_i)) = + $ tofmask(ch32b(tof32_i),hb32b(tof32_i)) + 1 + +C------- + ENDIF C S22 - S31 @@ -826,6 +923,21 @@ c1 = k_S22S31(1,ihelp) c2 = k_S22S31(2,ihelp) betatof_a(7) = c2/(cos(theta13)*(ds-c1)) + +C------- ToF Mask - S22 - S31 + + tofmask(ch22a(tof22_i),hb22a(tof22_i)) = + $ tofmask(ch22a(tof22_i),hb22a(tof22_i)) + 1 + tofmask(ch22b(tof22_i),hb22b(tof22_i)) = + $ tofmask(ch22b(tof22_i),hb22b(tof22_i)) + 1 + + tofmask(ch31a(tof31_i),hb31a(tof31_i)) = + $ tofmask(ch31a(tof31_i),hb31a(tof31_i)) + 1 + tofmask(ch31b(tof31_i),hb31b(tof31_i)) = + $ tofmask(ch31b(tof31_i),hb31b(tof31_i)) + 1 + +C------- + ENDIF C S22 - S32 @@ -837,6 +949,21 @@ c1 = k_S22S32(1,ihelp) c2 = k_S22S32(2,ihelp) betatof_a(8) = c2/(cos(theta13)*(ds-c1)) + +C------- ToF Mask - S22 - S32 + + tofmask(ch22a(tof22_i),hb22a(tof22_i)) = + $ tofmask(ch22a(tof22_i),hb22a(tof22_i)) + 1 + tofmask(ch22b(tof22_i),hb22b(tof22_i)) = + $ tofmask(ch22b(tof22_i),hb22b(tof22_i)) + 1 + + tofmask(ch32a(tof32_i),hb32a(tof32_i)) = + $ tofmask(ch32a(tof32_i),hb32a(tof32_i)) + 1 + tofmask(ch32b(tof32_i),hb32b(tof32_i)) = + $ tofmask(ch32b(tof32_i),hb32b(tof32_i)) + 1 + +C------- + ENDIF C S11 - S21 @@ -848,6 +975,21 @@ c1 = k_S11S21(1,ihelp) c2 = k_S11S21(2,ihelp) betatof_a(9) = c2/(cos(theta13)*(ds-c1)) + +C------- ToF Mask - S11 - S21 + + tofmask(ch11a(tof11_i),hb11a(tof11_i)) = + $ tofmask(ch11a(tof11_i),hb11a(tof11_i)) + 1 + tofmask(ch11b(tof11_i),hb11b(tof11_i)) = + $ tofmask(ch11b(tof11_i),hb11b(tof11_i)) + 1 + + tofmask(ch21a(tof21_i),hb21a(tof21_i)) = + $ tofmask(ch21a(tof21_i),hb21a(tof21_i)) + 1 + tofmask(ch21b(tof21_i),hb21b(tof21_i)) = + $ tofmask(ch21b(tof21_i),hb21b(tof21_i)) + 1 + +C------- + ENDIF C S11 - S22 @@ -859,6 +1001,21 @@ c1 = k_S11S22(1,ihelp) c2 = k_S11S22(2,ihelp) betatof_a(10) = c2/(cos(theta13)*(ds-c1)) + +C------- ToF Mask - S11 - S22 + + tofmask(ch11a(tof11_i),hb11a(tof11_i)) = + $ tofmask(ch11a(tof11_i),hb11a(tof11_i)) + 1 + tofmask(ch11b(tof11_i),hb11b(tof11_i)) = + $ tofmask(ch11b(tof11_i),hb11b(tof11_i)) + 1 + + tofmask(ch22a(tof22_i),hb22a(tof22_i)) = + $ tofmask(ch22a(tof22_i),hb22a(tof22_i)) + 1 + tofmask(ch22b(tof22_i),hb22b(tof22_i)) = + $ tofmask(ch22b(tof22_i),hb22b(tof22_i)) + 1 + +C------- + ENDIF C S12 - S21 @@ -870,6 +1027,21 @@ c1 = k_S12S21(1,ihelp) c2 = k_S12S21(2,ihelp) betatof_a(11) = c2/(cos(theta13)*(ds-c1)) + +C------- ToF Mask - S12 - S21 + + tofmask(ch12a(tof12_i),hb12a(tof12_i)) = + $ tofmask(ch12a(tof12_i),hb12a(tof12_i)) + 1 + tofmask(ch12b(tof12_i),hb12b(tof12_i)) = + $ tofmask(ch12b(tof12_i),hb12b(tof12_i)) + 1 + + tofmask(ch21a(tof21_i),hb21a(tof21_i)) = + $ tofmask(ch21a(tof21_i),hb21a(tof21_i)) + 1 + tofmask(ch21b(tof21_i),hb21b(tof21_i)) = + $ tofmask(ch21b(tof21_i),hb21b(tof21_i)) + 1 + +C------- + ENDIF C S12 - S22 @@ -881,9 +1053,24 @@ c1 = k_S12S22(1,ihelp) c2 = k_S12S22(2,ihelp) betatof_a(12) = c2/(cos(theta13)*(ds-c1)) + +C------- ToF Mask - S12 - S22 + + tofmask(ch12a(tof12_i),hb12a(tof12_i)) = + $ tofmask(ch12a(tof12_i),hb12a(tof12_i)) + 1 + tofmask(ch12b(tof12_i),hb12b(tof12_i)) = + $ tofmask(ch12b(tof12_i),hb12b(tof12_i)) + 1 + + tofmask(ch22a(tof22_i),hb22a(tof22_i)) = + $ tofmask(ch22a(tof22_i),hb22a(tof22_i)) + 1 + tofmask(ch22b(tof22_i),hb22b(tof22_i)) = + $ tofmask(ch22b(tof22_i),hb22b(tof22_i)) + 1 + +C------- + ENDIF - -C------- + +C--------------------------------------------------------- icount=0 sw=0. @@ -903,7 +1090,6 @@ if (icount.gt.0) beta_mean=sxw/sw betatof_a(13) = beta_mean -c write(*,*) '------------- end tofl2com ----------' 100 continue