*--------1---------2---------3---------4---------5---------6---------7-- * Generating the STT LIET masks for * the PAD-SPACAL topological match * declare variables *---------------------------------- real function TOPO(void) integer void include './Include/common.inc' integer XY(20, 20) * book histograms *---------------------------------- if (IDNEVT.eq.1) then call HBOOK2(1, 'S0', 20, -10.,10., 20, -10.,10., 0) call HBOOK2(2, 'S1', 20, -10.,10., 20, -10.,10., 0) call HBOOK2(3, 'S2', 20, -10.,10., 20, -10.,10., 0) call HBOOK2(4, 'S3', 20, -10.,10., 20, -10.,10., 0) call HBOOK2(5, 'S4', 20, -10.,10., 20, -10.,10., 0) call HBOOK2(6, 'S5', 20, -10.,10., 20, -10.,10., 0) call HBOOK2(7, 'S10', 20, -10.,10., 20, -10.,10., 0) call HBOOK2(8, 'S11', 20, -10.,10., 20, -10.,10., 0) call HBOOK2(9, 'S12', 20, -10.,10., 20, -10.,10., 0) call HBOOK2(10, 'S13', 20, -10.,10., 20, -10.,10., 0) call HBOOK2(11, 'S14', 20, -10.,10., 20, -10.,10., 0) call HBOOK2(12, 'S15', 20, -10.,10., 20, -10.,10., 0) call HBOOK2(13, 'All', 20, -10.,10., 20, -10.,10., 0) end if * convert LIET bits to XY grid *---------------------------------- XY(1, 1) = IBITS(LIET(1), 0, 1) XY(1, 2) = IBITS(LIET(1), 1, 1) XY(1, 3) = IBITS(LIET(1), 2, 1) XY(1, 4) = IBITS(LIET(1), 3, 1) XY(1, 5) = IBITS(LIET(2), 0, 1) XY(1, 6) = IBITS(LIET(2), 1, 1) XY(1, 7) = IBITS(LIET(2), 2, 1) XY(1, 8) = IBITS(LIET(2), 3, 1) XY(1, 9) = IBITS(LIET(3), 0, 1) XY(1, 10) = IBITS(LIET(3), 1, 1) XY(1, 11) = IBITS(LIET(3), 2, 1) XY(1, 12) = IBITS(LIET(3), 3, 1) XY(1, 13) = IBITS(LIET(4), 0, 1) XY(1, 14) = IBITS(LIET(4), 1, 1) XY(1, 15) = IBITS(LIET(4), 2, 1) XY(1, 16) = IBITS(LIET(4), 3, 1) XY(1, 17) = IBITS(LIET(5), 0, 1) XY(1, 18) = IBITS(LIET(5), 1, 1) XY(1, 19) = IBITS(LIET(5), 2, 1) XY(1, 20) = IBITS(LIET(5), 3, 1) * XY(2, 1) = IBITS(LIET(1), 4, 1) XY(2, 2) = IBITS(LIET(1), 5, 1) XY(2, 3) = IBITS(LIET(1), 6, 1) XY(2, 4) = IBITS(LIET(1), 7, 1) XY(2, 5) = IBITS(LIET(2), 4, 1) XY(2, 6) = IBITS(LIET(2), 5, 1) XY(2, 7) = IBITS(LIET(2), 6, 1) XY(2, 8) = IBITS(LIET(2), 7, 1) XY(2, 9) = IBITS(LIET(3), 4, 1) XY(2, 10) = IBITS(LIET(3), 5, 1) XY(2, 11) = IBITS(LIET(3), 6, 1) XY(2, 12) = IBITS(LIET(3), 7, 1) XY(2, 13) = IBITS(LIET(4), 4, 1) XY(2, 14) = IBITS(LIET(4), 5, 1) XY(2, 15) = IBITS(LIET(4), 6, 1) XY(2, 16) = IBITS(LIET(4), 7, 1) XY(2, 17) = IBITS(LIET(5), 4, 1) XY(2, 18) = IBITS(LIET(5), 5, 1) XY(2, 19) = IBITS(LIET(5), 6, 1) XY(2, 20) = IBITS(LIET(5), 7, 1) * XY(3, 1) = IBITS(LIET(1), 8, 1) XY(3, 2) = IBITS(LIET(1), 9, 1) XY(3, 3) = IBITS(LIET(1), 10, 1) XY(3, 4) = IBITS(LIET(1), 11, 1) XY(3, 5) = IBITS(LIET(2), 8, 1) XY(3, 6) = IBITS(LIET(2), 9, 1) XY(3, 7) = IBITS(LIET(2), 10, 1) XY(3, 8) = IBITS(LIET(2), 11, 1) XY(3, 9) = IBITS(LIET(3), 8, 1) XY(3, 10) = IBITS(LIET(3), 9, 1) XY(3, 11) = IBITS(LIET(3), 10, 1) XY(3, 12) = IBITS(LIET(3), 11, 1) XY(3, 13) = IBITS(LIET(4), 8, 1) XY(3, 14) = IBITS(LIET(4), 9, 1) XY(3, 15) = IBITS(LIET(4), 10, 1) XY(3, 16) = IBITS(LIET(4), 11, 1) XY(3, 17) = IBITS(LIET(5), 8, 1) XY(3, 18) = IBITS(LIET(5), 9, 1) XY(3, 19) = IBITS(LIET(5), 10, 1) XY(3, 20) = IBITS(LIET(5), 11, 1) * XY(4, 1) = IBITS(LIET(1), 12, 1) XY(4, 2) = IBITS(LIET(1), 13, 1) XY(4, 3) = IBITS(LIET(1), 14, 1) XY(4, 4) = IBITS(LIET(1), 15, 1) XY(4, 5) = IBITS(LIET(2), 12, 1) XY(4, 6) = IBITS(LIET(2), 13, 1) XY(4, 7) = IBITS(LIET(2), 14, 1) XY(4, 8) = IBITS(LIET(2), 15, 1) XY(4, 9) = IBITS(LIET(3), 12, 1) XY(4, 10) = IBITS(LIET(3), 13, 1) XY(4, 11) = IBITS(LIET(3), 14, 1) XY(4, 12) = IBITS(LIET(3), 15, 1) XY(4, 13) = IBITS(LIET(4), 12, 1) XY(4, 14) = IBITS(LIET(4), 13, 1) XY(4, 15) = IBITS(LIET(4), 14, 1) XY(4, 16) = IBITS(LIET(4), 15, 1) XY(4, 17) = IBITS(LIET(5), 12, 1) XY(4, 18) = IBITS(LIET(5), 13, 1) XY(4, 19) = IBITS(LIET(5), 14, 1) XY(4, 20) = IBITS(LIET(5), 15, 1) * XY(5, 1) = IBITS(LIET(6), 0, 1) XY(5, 2) = IBITS(LIET(6), 1, 1) XY(5, 3) = IBITS(LIET(6), 2, 1) XY(5, 4) = IBITS(LIET(6), 3, 1) XY(5, 5) = IBITS(LIET(7), 0, 1) XY(5, 6) = IBITS(LIET(7), 1, 1) XY(5, 7) = IBITS(LIET(7), 2, 1) XY(5, 8) = IBITS(LIET(7), 3, 1) XY(5, 9) = IBITS(LIET(8), 0, 1) XY(5, 10) = IBITS(LIET(8), 1, 1) XY(5, 11) = IBITS(LIET(8), 2, 1) XY(5, 12) = IBITS(LIET(8), 3, 1) XY(5, 13) = IBITS(LIET(9), 0, 1) XY(5, 14) = IBITS(LIET(9), 1, 1) XY(5, 15) = IBITS(LIET(9), 2, 1) XY(5, 16) = IBITS(LIET(9), 3, 1) XY(5, 17) = IBITS(LIET(10), 0, 1) XY(5, 18) = IBITS(LIET(10), 1, 1) XY(5, 19) = IBITS(LIET(10), 2, 1) XY(5, 20) = IBITS(LIET(10), 3, 1) * XY(6, 1) = IBITS(LIET(6), 4, 1) XY(6, 2) = IBITS(LIET(6), 5, 1) XY(6, 3) = IBITS(LIET(6), 6, 1) XY(6, 4) = IBITS(LIET(6), 7, 1) XY(6, 5) = IBITS(LIET(7), 4, 1) XY(6, 6) = IBITS(LIET(7), 5, 1) XY(6, 7) = IBITS(LIET(7), 6, 1) XY(6, 8) = IBITS(LIET(7), 7, 1) XY(6, 9) = IBITS(LIET(8), 4, 1) XY(6, 10) = IBITS(LIET(8), 5, 1) XY(6, 11) = IBITS(LIET(8), 6, 1) XY(6, 12) = IBITS(LIET(8), 7, 1) XY(6, 13) = IBITS(LIET(9), 4, 1) XY(6, 14) = IBITS(LIET(9), 5, 1) XY(6, 15) = IBITS(LIET(9), 6, 1) XY(6, 16) = IBITS(LIET(9), 7, 1) XY(6, 17) = IBITS(LIET(10), 4, 1) XY(6, 18) = IBITS(LIET(10), 5, 1) XY(6, 19) = IBITS(LIET(10), 6, 1) XY(6, 20) = IBITS(LIET(10), 7, 1) * XY(7, 1) = IBITS(LIET(6), 8, 1) XY(7, 2) = IBITS(LIET(6), 9, 1) XY(7, 3) = IBITS(LIET(6), 10, 1) XY(7, 4) = IBITS(LIET(6), 11, 1) XY(7, 5) = IBITS(LIET(7), 8, 1) XY(7, 6) = IBITS(LIET(7), 9, 1) XY(7, 7) = IBITS(LIET(7), 10, 1) XY(7, 8) = IBITS(LIET(7), 11, 1) XY(7, 9) = IBITS(LIET(8), 8, 1) XY(7, 10) = IBITS(LIET(8), 9, 1) XY(7, 11) = IBITS(LIET(8), 10, 1) XY(7, 12) = IBITS(LIET(8), 11, 1) XY(7, 13) = IBITS(LIET(9), 8, 1) XY(7, 14) = IBITS(LIET(9), 9, 1) XY(7, 15) = IBITS(LIET(9), 10, 1) XY(7, 16) = IBITS(LIET(9), 11, 1) XY(7, 17) = IBITS(LIET(10), 8, 1) XY(7, 18) = IBITS(LIET(10), 9, 1) XY(7, 19) = IBITS(LIET(10), 10, 1) XY(7, 20) = IBITS(LIET(10), 11, 1) * XY(8, 1) = IBITS(LIET(6), 12, 1) XY(8, 2) = IBITS(LIET(6), 13, 1) XY(8, 3) = IBITS(LIET(6), 14, 1) XY(8, 4) = IBITS(LIET(6), 15, 1) XY(8, 5) = IBITS(LIET(7), 12, 1) XY(8, 6) = IBITS(LIET(7), 13, 1) XY(8, 7) = IBITS(LIET(7), 14, 1) XY(8, 8) = IBITS(LIET(7), 15, 1) XY(8, 9) = IBITS(LIET(8), 12, 1) XY(8, 10) = IBITS(LIET(8), 13, 1) XY(8, 11) = IBITS(LIET(8), 14, 1) XY(8, 12) = IBITS(LIET(8), 15, 1) XY(8, 13) = IBITS(LIET(9), 12, 1) XY(8, 14) = IBITS(LIET(9), 13, 1) XY(8, 15) = IBITS(LIET(9), 14, 1) XY(8, 16) = IBITS(LIET(9), 15, 1) XY(8, 17) = IBITS(LIET(10), 12, 1) XY(8, 18) = IBITS(LIET(10), 13, 1) XY(8, 19) = IBITS(LIET(10), 14, 1) XY(8, 20) = IBITS(LIET(10), 15, 1) * XY(9, 1) = IBITS(LIET(11), 0, 1) XY(9, 2) = IBITS(LIET(11), 1, 1) XY(9, 3) = IBITS(LIET(11), 2, 1) XY(9, 4) = IBITS(LIET(11), 3, 1) XY(9, 5) = IBITS(LIET(12), 0, 1) XY(9, 6) = IBITS(LIET(12), 1, 1) XY(9, 7) = IBITS(LIET(12), 2, 1) XY(9, 8) = IBITS(LIET(12), 3, 1) XY(9, 9) = IBITS(LIET(13), 0, 1) XY(9, 10) = IBITS(LIET(13), 1, 1) XY(9, 11) = IBITS(LIET(13), 2, 1) XY(9, 12) = IBITS(LIET(13), 3, 1) XY(9, 13) = IBITS(LIET(14), 0, 1) XY(9, 14) = IBITS(LIET(14), 1, 1) XY(9, 15) = IBITS(LIET(14), 2, 1) XY(9, 16) = IBITS(LIET(14), 3, 1) XY(9, 17) = IBITS(LIET(15), 0, 1) XY(9, 18) = IBITS(LIET(15), 1, 1) XY(9, 19) = IBITS(LIET(15), 2, 1) XY(9, 20) = IBITS(LIET(15), 3, 1) * XY(10, 1) = IBITS(LIET(11), 4, 1) XY(10, 2) = IBITS(LIET(11), 5, 1) XY(10, 3) = IBITS(LIET(11), 6, 1) XY(10, 4) = IBITS(LIET(11), 7, 1) XY(10, 5) = IBITS(LIET(12), 4, 1) XY(10, 6) = IBITS(LIET(12), 5, 1) XY(10, 7) = IBITS(LIET(12), 6, 1) XY(10, 8) = IBITS(LIET(12), 7, 1) XY(10, 9) = IBITS(LIET(13), 4, 1) XY(10, 10) = IBITS(LIET(13), 5, 1) XY(10, 11) = IBITS(LIET(13), 6, 1) XY(10, 12) = IBITS(LIET(13), 7, 1) XY(10, 13) = IBITS(LIET(14), 4, 1) XY(10, 14) = IBITS(LIET(14), 5, 1) XY(10, 15) = IBITS(LIET(14), 6, 1) XY(10, 16) = IBITS(LIET(14), 7, 1) XY(10, 17) = IBITS(LIET(15), 4, 1) XY(10, 18) = IBITS(LIET(15), 5, 1) XY(10, 19) = IBITS(LIET(15), 6, 1) XY(10, 20) = IBITS(LIET(15), 7, 1) * XY(11, 1) = IBITS(LIET(11), 8, 1) XY(11, 2) = IBITS(LIET(11), 9, 1) XY(11, 3) = IBITS(LIET(11), 10, 1) XY(11, 4) = IBITS(LIET(11), 11, 1) XY(11, 5) = IBITS(LIET(12), 8, 1) XY(11, 6) = IBITS(LIET(12), 9, 1) XY(11, 7) = IBITS(LIET(12), 10, 1) XY(11, 8) = IBITS(LIET(12), 11, 1) XY(11, 9) = IBITS(LIET(13), 8, 1) XY(11, 10) = IBITS(LIET(13), 9, 1) XY(11, 11) = IBITS(LIET(13), 10, 1) XY(11, 12) = IBITS(LIET(13), 11, 1) XY(11, 13) = IBITS(LIET(14), 8, 1) XY(11, 14) = IBITS(LIET(14), 9, 1) XY(11, 15) = IBITS(LIET(14), 10, 1) XY(11, 16) = IBITS(LIET(14), 11, 1) XY(11, 17) = IBITS(LIET(15), 8, 1) XY(11, 18) = IBITS(LIET(15), 9, 1) XY(11, 19) = IBITS(LIET(15), 10, 1) XY(11, 20) = IBITS(LIET(15), 11, 1) * XY(12, 1) = IBITS(LIET(11), 12, 1) XY(12, 2) = IBITS(LIET(11), 13, 1) XY(12, 3) = IBITS(LIET(11), 14, 1) XY(12, 4) = IBITS(LIET(11), 15, 1) XY(12, 5) = IBITS(LIET(12), 12, 1) XY(12, 6) = IBITS(LIET(12), 13, 1) XY(12, 7) = IBITS(LIET(12), 14, 1) XY(12, 8) = IBITS(LIET(12), 15, 1) XY(12, 9) = IBITS(LIET(13), 12, 1) XY(12, 10) = IBITS(LIET(13), 13, 1) XY(12, 11) = IBITS(LIET(13), 14, 1) XY(12, 12) = IBITS(LIET(13), 15, 1) XY(12, 13) = IBITS(LIET(14), 12, 1) XY(12, 14) = IBITS(LIET(14), 13, 1) XY(12, 15) = IBITS(LIET(14), 14, 1) XY(12, 16) = IBITS(LIET(14), 15, 1) XY(12, 17) = IBITS(LIET(15), 12, 1) XY(12, 18) = IBITS(LIET(15), 13, 1) XY(12, 19) = IBITS(LIET(15), 14, 1) XY(12, 20) = IBITS(LIET(15), 15, 1) * XY(13, 1) = IBITS(LIET(16), 0, 1) XY(13, 2) = IBITS(LIET(16), 1, 1) XY(13, 3) = IBITS(LIET(16), 2, 1) XY(13, 4) = IBITS(LIET(16), 3, 1) XY(13, 5) = IBITS(LIET(17), 0, 1) XY(13, 6) = IBITS(LIET(17), 1, 1) XY(13, 7) = IBITS(LIET(17), 2, 1) XY(13, 8) = IBITS(LIET(17), 3, 1) XY(13, 9) = IBITS(LIET(18), 0, 1) XY(13, 10) = IBITS(LIET(18), 1, 1) XY(13, 11) = IBITS(LIET(18), 2, 1) XY(13, 12) = IBITS(LIET(18), 3, 1) XY(13, 13) = IBITS(LIET(19), 0, 1) XY(13, 14) = IBITS(LIET(19), 1, 1) XY(13, 15) = IBITS(LIET(19), 2, 1) XY(13, 16) = IBITS(LIET(19), 3, 1) XY(13, 17) = IBITS(LIET(20), 0, 1) XY(13, 18) = IBITS(LIET(20), 1, 1) XY(13, 19) = IBITS(LIET(20), 2, 1) XY(13, 20) = IBITS(LIET(20), 3, 1) * XY(14, 1) = IBITS(LIET(16), 4, 1) XY(14, 2) = IBITS(LIET(16), 5, 1) XY(14, 3) = IBITS(LIET(16), 6, 1) XY(14, 4) = IBITS(LIET(16), 7, 1) XY(14, 5) = IBITS(LIET(17), 4, 1) XY(14, 6) = IBITS(LIET(17), 5, 1) XY(14, 7) = IBITS(LIET(17), 6, 1) XY(14, 8) = IBITS(LIET(17), 7, 1) XY(14, 9) = IBITS(LIET(18), 4, 1) XY(14, 10) = IBITS(LIET(18), 5, 1) XY(14, 11) = IBITS(LIET(18), 6, 1) XY(14, 12) = IBITS(LIET(18), 7, 1) XY(14, 13) = IBITS(LIET(19), 4, 1) XY(14, 14) = IBITS(LIET(19), 5, 1) XY(14, 15) = IBITS(LIET(19), 6, 1) XY(14, 16) = IBITS(LIET(19), 7, 1) XY(14, 17) = IBITS(LIET(20), 4, 1) XY(14, 18) = IBITS(LIET(20), 5, 1) XY(14, 19) = IBITS(LIET(20), 6, 1) XY(14, 20) = IBITS(LIET(20), 7, 1) * XY(15, 1) = IBITS(LIET(16), 8, 1) XY(15, 2) = IBITS(LIET(16), 9, 1) XY(15, 3) = IBITS(LIET(16), 10, 1) XY(15, 4) = IBITS(LIET(16), 11, 1) XY(15, 5) = IBITS(LIET(17), 8, 1) XY(15, 6) = IBITS(LIET(17), 9, 1) XY(15, 7) = IBITS(LIET(17), 10, 1) XY(15, 8) = IBITS(LIET(17), 11, 1) XY(15, 9) = IBITS(LIET(18), 8, 1) XY(15, 10) = IBITS(LIET(18), 9, 1) XY(15, 11) = IBITS(LIET(18), 10, 1) XY(15, 12) = IBITS(LIET(18), 11, 1) XY(15, 13) = IBITS(LIET(19), 8, 1) XY(15, 14) = IBITS(LIET(19), 9, 1) XY(15, 15) = IBITS(LIET(19), 10, 1) XY(15, 16) = IBITS(LIET(19), 11, 1) XY(15, 17) = IBITS(LIET(20), 8, 1) XY(15, 18) = IBITS(LIET(20), 9, 1) XY(15, 19) = IBITS(LIET(20), 10, 1) XY(15, 20) = IBITS(LIET(20), 11, 1) * XY(16, 1) = IBITS(LIET(16), 12, 1) XY(16, 2) = IBITS(LIET(16), 13, 1) XY(16, 3) = IBITS(LIET(16), 14, 1) XY(16, 4) = IBITS(LIET(16), 15, 1) XY(16, 5) = IBITS(LIET(17), 12, 1) XY(16, 6) = IBITS(LIET(17), 13, 1) XY(16, 7) = IBITS(LIET(17), 14, 1) XY(16, 8) = IBITS(LIET(17), 15, 1) XY(16, 9) = IBITS(LIET(18), 12, 1) XY(16, 10) = IBITS(LIET(18), 13, 1) XY(16, 11) = IBITS(LIET(18), 14, 1) XY(16, 12) = IBITS(LIET(18), 15, 1) XY(16, 13) = IBITS(LIET(19), 12, 1) XY(16, 14) = IBITS(LIET(19), 13, 1) XY(16, 15) = IBITS(LIET(19), 14, 1) XY(16, 16) = IBITS(LIET(19), 15, 1) XY(16, 17) = IBITS(LIET(20), 12, 1) XY(16, 18) = IBITS(LIET(20), 13, 1) XY(16, 19) = IBITS(LIET(20), 14, 1) XY(16, 20) = IBITS(LIET(20), 15, 1) * XY(17, 1) = IBITS(LIET(21), 0, 1) XY(17, 2) = IBITS(LIET(21), 1, 1) XY(17, 3) = IBITS(LIET(21), 2, 1) XY(17, 4) = IBITS(LIET(21), 3, 1) XY(17, 5) = IBITS(LIET(22), 0, 1) XY(17, 6) = IBITS(LIET(22), 1, 1) XY(17, 7) = IBITS(LIET(22), 2, 1) XY(17, 8) = IBITS(LIET(22), 3, 1) XY(17, 9) = IBITS(LIET(23), 0, 1) XY(17, 10) = IBITS(LIET(23), 1, 1) XY(17, 11) = IBITS(LIET(23), 2, 1) XY(17, 12) = IBITS(LIET(23), 3, 1) XY(17, 13) = IBITS(LIET(24), 0, 1) XY(17, 14) = IBITS(LIET(24), 1, 1) XY(17, 15) = IBITS(LIET(24), 2, 1) XY(17, 16) = IBITS(LIET(24), 3, 1) XY(17, 17) = IBITS(LIET(25), 0, 1) XY(17, 18) = IBITS(LIET(25), 1, 1) XY(17, 19) = IBITS(LIET(25), 2, 1) XY(17, 20) = IBITS(LIET(25), 3, 1) * XY(18, 1) = IBITS(LIET(21), 4, 1) XY(18, 2) = IBITS(LIET(21), 5, 1) XY(18, 3) = IBITS(LIET(21), 6, 1) XY(18, 4) = IBITS(LIET(21), 7, 1) XY(18, 5) = IBITS(LIET(22), 4, 1) XY(18, 6) = IBITS(LIET(22), 5, 1) XY(18, 7) = IBITS(LIET(22), 6, 1) XY(18, 8) = IBITS(LIET(22), 7, 1) XY(18, 9) = IBITS(LIET(23), 4, 1) XY(18, 10) = IBITS(LIET(23), 5, 1) XY(18, 11) = IBITS(LIET(23), 6, 1) XY(18, 12) = IBITS(LIET(23), 7, 1) XY(18, 13) = IBITS(LIET(24), 4, 1) XY(18, 14) = IBITS(LIET(24), 5, 1) XY(18, 15) = IBITS(LIET(24), 6, 1) XY(18, 16) = IBITS(LIET(24), 7, 1) XY(18, 17) = IBITS(LIET(25), 4, 1) XY(18, 18) = IBITS(LIET(25), 5, 1) XY(18, 19) = IBITS(LIET(25), 6, 1) XY(18, 20) = IBITS(LIET(25), 7, 1) * XY(19, 1) = IBITS(LIET(21), 8, 1) XY(19, 2) = IBITS(LIET(21), 9, 1) XY(19, 3) = IBITS(LIET(21), 10, 1) XY(19, 4) = IBITS(LIET(21), 11, 1) XY(19, 5) = IBITS(LIET(22), 8, 1) XY(19, 6) = IBITS(LIET(22), 9, 1) XY(19, 7) = IBITS(LIET(22), 10, 1) XY(19, 8) = IBITS(LIET(22), 11, 1) XY(19, 9) = IBITS(LIET(23), 8, 1) XY(19, 10) = IBITS(LIET(23), 9, 1) XY(19, 11) = IBITS(LIET(23), 10, 1) XY(19, 12) = IBITS(LIET(23), 11, 1) XY(19, 13) = IBITS(LIET(24), 8, 1) XY(19, 14) = IBITS(LIET(24), 9, 1) XY(19, 15) = IBITS(LIET(24), 10, 1) XY(19, 16) = IBITS(LIET(24), 11, 1) XY(19, 17) = IBITS(LIET(25), 8, 1) XY(19, 18) = IBITS(LIET(25), 9, 1) XY(19, 19) = IBITS(LIET(25), 10, 1) XY(19, 20) = IBITS(LIET(25), 11, 1) * XY(20, 1) = IBITS(LIET(21), 12, 1) XY(20, 2) = IBITS(LIET(21), 13, 1) XY(20, 3) = IBITS(LIET(21), 14, 1) XY(20, 4) = IBITS(LIET(21), 15, 1) XY(20, 5) = IBITS(LIET(22), 12, 1) XY(20, 6) = IBITS(LIET(22), 13, 1) XY(20, 7) = IBITS(LIET(22), 14, 1) XY(20, 8) = IBITS(LIET(22), 15, 1) XY(20, 9) = IBITS(LIET(23), 12, 1) XY(20, 10) = IBITS(LIET(23), 13, 1) XY(20, 11) = IBITS(LIET(23), 14, 1) XY(20, 12) = IBITS(LIET(23), 15, 1) XY(20, 13) = IBITS(LIET(24), 12, 1) XY(20, 14) = IBITS(LIET(24), 13, 1) XY(20, 15) = IBITS(LIET(24), 14, 1) XY(20, 16) = IBITS(LIET(24), 15, 1) XY(20, 17) = IBITS(LIET(25), 12, 1) XY(20, 18) = IBITS(LIET(25), 13, 1) XY(20, 19) = IBITS(LIET(25), 14, 1) XY(20, 20) = IBITS(LIET(25), 15, 1) * filling histograms * K = - X direction * J = Y direction *---------------------------------- do I = 1, 12 if (Top_Trg(I).eq.1) then do J = 1, 20 do K = 1, 20 YDIR = +(real(J) - 10.) XDIR = -(real(K) - 10.) if (XY(J, K).eq.1) call HF2(I, XDIR, YDIR, 1.) if (XY(J, K).eq.1) call HF2(13, XDIR, YDIR, 1.) end do end do end if end do return end