CalorimeterLevel2/src/laterale.for

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,NPLA)
      REAL BAR(2,NPLA)
      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

1	mocchiut	1.1	C
2			C---------------------------------------------------------------------
3			SUBROUTINE LATERALE(RQT1,RQT2)
4			C---------------------------------------------------------------------
5			C RQT1 (IT WILL BE CALLED QTRACK IN THE N-TUPLE) IS THE SUM OF THE DETECTED
6			C ENERGY IN THE STRIP ALONG THE TRACK AND THE TWO CLOSEST STRIPS . FOR ALL THE
7			C LAYERS . RQT2 (IS NOT USED IN THE N-TUPLA) IS THE TOTAL ENERGY MINUS RQT1 .
8			C
9			INCLUDE 'INTEST.TXT'
10			C
11			INTEGER A,B, INPIA
12			INTEGER IBAR(2,NPLA)
13			REAL BAR(2,NPLA)
14			REAL RQT1, RQT2
15			COMMON / ANGOLO / BAR, IBAR
16			SAVE / ANGOLO /
17			C
18			RQT1 = 0.
19			RQT2 = 0.
20			C
21			INPIA = 1
22			C
23			DO I = INPIA, NPLA
24			A = IBAR(1,I)
25			C
26			C EM
27			C
28			IF (A.EQ.-1) GOTO 600
29			C
30			C END EM
31			C
32			IF (A.LE.2) A = 3
33			IF (A.GE.(NCHA-1)) A = NCHA - 2
34			DO J = A-1,A+1
35			IF (DEXY(1,I,J).GE.EMIN) RQT1 = RQT1 + DEXY(1,I,J)
36			ENDDO
37			C
38			DO J=1,A-2
39			PXY = DEXY(1,I,J)
40			IF (PXY.GE.EMIN) RQT2 = RQT2 + PXY
41			ENDDO
42			C
43			DO J=A+2,NCHA
44			PXY = DEXY(1,I,J)
45			IF (PXY.GE.EMIN) RQT2 = RQT2 + PXY
46			ENDDO
47			C
48			600 CONTINUE
49			C
50			B = IBAR(2,I)
51			C
52			C EM
53			C
54			IF (B.EQ.-1) GOTO 700
55			C
56			C END EM
57			C
58			IF (B.LE.2) B = 3
59			IF (B.GE.(NCHA-1)) B = NCHA - 2
60			DO J = B-1,B+1
61			IF (DEXY(2,I,J).GE.EMIN) RQT1 = RQT1 + DEXY(2,I,J)
62			ENDDO
63			C
64			DO J=1,B-2
65			PXY = DEXY(2,I,J)
66			IF (PXY.GE.EMIN) RQT2 = RQT2 + PXY
67			ENDDO
68			C
69			DO J=B+2,NCHA
70			PXY = DEXY(2,I,J)
71			IF (PXY.GE.EMIN) RQT2 = RQT2 + PXY
72			ENDDO
73			C
74			700 CONTINUE
75			ENDDO
76			C
77			C
78			800 RETURN
79			END
80