*--------1---------2---------3---------4---------5---------6---------7-- * Trigger rate vs. SPACAL energy * to study the topological cuts * declare variables *---------------------------------- real function TRIG(L22_ON_OFF) include './Include/common.inc' include './Include/padgeom.inc' integer L22_ON_OFF ! 0 = OFF, 1 = ON real RAD, PHI, Pi, NLIET, HPAR parameter (Pi = 3.14159265358979) logical S7, S8, L2_22, PAD, ACC, SV, STT integer XY(20, 20) * book histograms *---------------------------------- if (IDNEVT.eq.1) then call HBOOK1(1, 'S8 trigger alone', 20, 0., 40.,0) call HBOOK1(2, 'S8 trigger & pad', 20, 0.,40.,0) call HBOOK1(11, '', 20, 0., 40.,0) call HBOOK1(12, 'S7 L2 validated', 20, 0., 40.,0) call HBOOK1(13, 'S7 L2 validated', 20, 0., 40.,0) call HBOOK2(21, 'XY acceptance', 20, -80.,80., 20,-80.,80., 0) call HBOOK2(22, 'XY spread', 20, -80.,80., 20,-80.,80., 0) call HBOOK2(23, 'XY spread valid', 20, -80.,80., 20,-80.,80., 0) call HBOOK1(31, 'N clusters', 20, 0., 80., 0) call HBOOK1(32, 'N clusters', 20, 0., 80., 0) end if if (L22_ON_OFF.eq.0) LIET(13) = 0 * L1_S8 condition *---------------------------------- if (NSEL.gt.0.and.ESEL(1).gt.2.and.S8ACT.eq.1) then S8 = .TRUE. else S8 = .FALSE. end if * L2_[22] condition *---------------------------------- RAD = sqrt(XSEL(1)**2+YSEL(1)**2) if (RAD.gt.20) then L2_22 = .TRUE. else L2_22 = .FALSE. end if * pad acceptance validation *---------------------------------- PHI = 180. * atan2(YSEL(1), XSEL(1)) / Pi if (PHI.lt.0.) PHI = 360. + PHI if (RAD.lt.32.and.(PHI.lt.135.or.PHI.gt.225)) then ACC = .TRUE. else ACC = .FALSE. end if * L1_PAD condition *---------------------------------- if (PADTRG.eq.1) then PAD = .TRUE. else PAD = .FALSE. end if * L1_S7 condition *---------------------------------- if (NSEL.gt.0.and.ESEL(1).gt.2.and.S7ACT.eq.1) then S7 = .TRUE. else S7 = .FALSE. end if * offline sector validation for S7 *---------------------------------- SV = .FALSE. do I = 1, 6 K = 2 * I - 1 L = 2 * I if (PHI.gt.EvenMinPhi(I).and. + PHI.lt.EvenMaxPhi(I).and. + Top_Trg(K).eq.1) SV = .TRUE. if (PHI.gt.OddMinPhi(I).and. + PHI.lt.OddMaxPhi(I).and. + Top_Trg(L).eq.1) SV = .TRUE. end do * quasi-online sector validation for S7 *---------------------------------- STT = .FALSE. * Sector #0: if (Top_Trg(1).eq.1) then if ((iand(LIET(12), 61152).ne.0)) STT = .TRUE. ! bits (05, 06, 07,) 09, 10, 11, 13, 14, 15 if ((iand(LIET(13), 4368).ne.0)) STT = .TRUE. ! bits 04, 08, 12 if ((iand(LIET(16), 14).ne.0)) STT = .TRUE. ! bits 01, 02, 03 end if * Sector #1: if (Top_Trg(2).eq.1) then if ((iand(LIET(12), 60416).ne.0)) STT = .TRUE. ! bits 10, 11, 13, 14, 15 if ((iand(LIET(13), 4352).ne.0)) STT = .TRUE. ! bits 08, 12 if ((iand(LIET(16), 238).ne.0)) STT = .TRUE. ! bits 01, 02, 03, 05, 06, 07 if ((iand(LIET(17), 1).ne.0)) STT = .TRUE. ! bits 00 end if * Sector #2: if (Top_Trg(3).eq.1) then if ((iand(LIET(12), 49152).ne.0)) STT = .TRUE. ! bits 14, 15 if ((iand(LIET(13), 12288).ne.0)) STT = .TRUE. ! bits 12, 13 if ((iand(LIET(16), 2252).ne.0)) STT = .TRUE. ! bits 02, 03, 06, 07, 11 if ((iand(LIET(17), 307).ne.0)) STT = .TRUE. ! bits 00, 01, 04, 05, 08 end if * Sector #3: if (Top_Trg(4).eq.1) then if ((iand(LIET(13), 12288).ne.0)) STT = .TRUE. ! bits 12, 13 if ((iand(LIET(16), 2184).ne.0)) STT = .TRUE. ! bits 03, 07, 11 if ((iand(LIET(17), 887).ne.0)) STT = .TRUE. ! bits 00, 01, (02), 04, 05, (06), 08, 09 end if * Sector #4: if (Top_Trg(5).eq.1) then if ((iand(LIET(13), 24576).ne.0)) STT = .TRUE. ! bits 13, 14 if ((iand(LIET(17), 3822).ne.0)) STT = .TRUE. ! bits (01), 02, 03, (05), 06, 07, (09), 10, 11 end if * Sector #5: if (Top_Trg(6).eq.1) then if ((iand(LIET(13), 49152).ne.0)) STT = .TRUE. ! bits 14, 15 if ((iand(LIET(17), 3278).ne.0)) STT = .TRUE. ! bits (01), 02, 03, 06, 07, 10, 11 if ((iand(LIET(18), 273).ne.0)) STT = .TRUE. ! bits 00, 04, 08 end if * Sector #10: if (Top_Trg(7).eq.1) then if ((iand(LIET(12), 52224).ne.0)) STT = .TRUE. ! bits 10, 11, 14, 15 if ((iand(LIET(13), 4352).ne.0)) STT = .TRUE. ! bits 08, 12 if ((iand(LIET(22), 204).ne.0)) STT = .TRUE. ! bits 02, 03, 06, 07 if ((iand(LIET(23), 1).ne.0)) STT = .TRUE. ! bits 00 end if * Sector #11: if (Top_Trg(8).eq.1) then if ((iand(LIET(8), 61024).ne.0)) STT = .TRUE. ! bits (05), 06, (09), 10, 11, (13), 14, 15 if ((iand(LIET(13), 14).ne.0)) STT = .TRUE. ! bits 01, 02, 03 end if * Sector #12: if (Top_Trg(9).eq.1) then if ((iand(LIET(7), 34816).ne.0)) STT = .TRUE. ! bits 11, 15 if ((iand(LIET(8), 30512).ne.0)) STT = .TRUE. ! bits 04, 05, 08, 09, (10), 12, 13, (14) if ((iand(LIET(13), 7).ne.0)) STT = .TRUE. ! bits 00, 01, 02 end if * Sector #13: if (Top_Trg(10).eq.1) then if ((iand(LIET(7), 52224).ne.0)) STT = .TRUE. ! bits 10, 11, 14, 15 if ((iand(LIET(8) , 12544).ne.0)) STT = .TRUE. ! bits 08, 12, 13 if ((iand(LIET(12), 140).ne.0)) STT = .TRUE. ! bits 02, 03, 07 if ((iand(LIET(13), 19).ne.0)) STT = .TRUE. ! bits 00, 01, 04 end if * Sector #14: if (Top_Trg(11).eq.1) then if ((iand(LIET(7), 57344).ne.0)) STT = .TRUE. ! bits 13, 14, 15 if ((iand(LIET(12), 238).ne.0)) STT = .TRUE. ! bits 01, 02, 03, 05, 06, 07 if ((iand(LIET(13), 17).ne.0)) STT = .TRUE. ! bits 00, 04 end if * Sector #15: if (Top_Trg(12).eq.1) then if ((iand(LIET(12), 3822).ne.0)) STT = .TRUE. ! bits 01, 02, 03, 05, 06, 07, (09, 10, 11) if ((iand(LIET(13), 272).ne.0)) STT = .TRUE. ! bits 04, 08 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) NLIET = 0. do I = 1, 20 do J = 1, 20 if (XY(I, J).eq.1) NLIET = NLIET + 1. end do end do HPAR = real(NSEL) * filling histograms *---------------------------------- if (S8.and.L2_22.and.ACC) call HF1(1, ESEL(1), 1.) if (S8.and.L2_22.and.ACC.and.PAD) call HF1(2, ESEL(1), 1.) if (S7.and.L2_22) call HF1(11, ESEL(1), 1.) if (S7.and.L2_22.and.STT) call HF1(12, ESEL(1), 1.) if (S7.and.L2_22.and.SV) call HF1(13, ESEL(1), 1.) if (S8.and.L2_22.and.ACC) call HF2(21, XSEL(1), YSEL(1), 1.) if (S7.and.L2_22) call HF2(22, XSEL(1), YSEL(1), 1.) if (S7.and.L2_22.and.STT) call HF2(23, XSEL(1), YSEL(1), 1.) if (S7) call HF1(31, NLIET, 1.) if (S7) call HF1(32, HPAR, 1.) return end