C C--------------------------------------------------------------------- SUBROUTINE LATERALE(RQT1,RQT2) C--------------------------------------------------------------------- C RQT1 (IT WILL BE CALLED QTRACK IN THE N-TUPLE) IS THE SUM OF THE DETECTED C ENERGY IN THE STRIP ALONG THE TRACK AND THE TWO CLOSEST STRIPS . FOR ALL THE C LAYERS . RQT2 (IS NOT USED IN THE N-TUPLA) IS THE TOTAL ENERGY MINUS RQT1 . C INCLUDE 'INTEST.TXT' C INTEGER A,B, INPIA INTEGER IBAR(2,NPLAV) REAL BAR(2,NPLAV) REAL RQT1, RQT2 COMMON / ANGOLO / BAR, IBAR SAVE / ANGOLO / C RQT1 = 0. RQT2 = 0. C INPIA = 1 C DO I = INPIA, NPLA A = IBAR(1,I) C C EM C IF (A.EQ.-1) GOTO 600 C C END EM C IF (A.LE.2) A = 3 IF (A.GE.(NCHA-1)) A = NCHA - 2 DO J = A-1,A+1 IF (DEXY(1,I,J).GE.EMIN) RQT1 = RQT1 + DEXY(1,I,J) ENDDO C DO J=1,A-2 PXY = DEXY(1,I,J) IF (PXY.GE.EMIN) RQT2 = RQT2 + PXY ENDDO C DO J=A+2,NCHA PXY = DEXY(1,I,J) IF (PXY.GE.EMIN) RQT2 = RQT2 + PXY ENDDO C 600 CONTINUE C B = IBAR(2,I) C C EM C IF (B.EQ.-1) GOTO 700 C C END EM C IF (B.LE.2) B = 3 IF (B.GE.(NCHA-1)) B = NCHA - 2 DO J = B-1,B+1 IF (DEXY(2,I,J).GE.EMIN) RQT1 = RQT1 + DEXY(2,I,J) ENDDO C DO J=1,B-2 PXY = DEXY(2,I,J) IF (PXY.GE.EMIN) RQT2 = RQT2 + PXY ENDDO C DO J=B+2,NCHA PXY = DEXY(2,I,J) IF (PXY.GE.EMIN) RQT2 = RQT2 + PXY ENDDO C 700 CONTINUE ENDDO C C 800 RETURN END