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,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 = 2
         IF (A.GE.NCHA) A = NCHA - 1
         DO J = A-1,A+1
            IF (DEXY(1,I,J).GE.EMIN) RQT1 = RQT1 + DEXY(1,I,J)
         ENDDO
C
         IF ((A-2).GE.1) THEN
            DO J=1,A-2
               PXY = DEXY(1,I,J)
               IF (PXY.GE.EMIN) RQT2 = RQT2 + PXY
            ENDDO
         ENDIF
C
         IF ((A+2).LE.NCHA) THEN
            DO J=A+2,NCHA
               PXY = DEXY(1,I,J)
               IF (PXY.GE.EMIN) RQT2 = RQT2 + PXY
            ENDDO
         ENDIF
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 = 2
         IF (B.GE.NCHA) B = NCHA - 1
         DO J = B-1,B+1
            IF (DEXY(2,I,J).GE.EMIN) RQT1 = RQT1 + DEXY(2,I,J)
         ENDDO
C     
         IF ((B-2).GE.1) THEN
            DO J=1,B-2
               PXY = DEXY(2,I,J)
               IF (PXY.GE.EMIN) RQT2 = RQT2 + PXY
            ENDDO
         ENDIF
C     
         IF ((B+2).LE.NCHA) THEN
            DO J=B+2,NCHA
               PXY = DEXY(2,I,J)
               IF (PXY.GE.EMIN) RQT2 = RQT2 + PXY
            ENDDO
         ENDIF
C
 700     CONTINUE
      ENDDO
C     
C
 800  RETURN
      END

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,NPLAV)
13	REAL BAR(2,NPLAV)
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 = 2
33	IF (A.GE.NCHA) A = NCHA - 1
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	IF ((A-2).GE.1) THEN
39	DO J=1,A-2
40	PXY = DEXY(1,I,J)
41	IF (PXY.GE.EMIN) RQT2 = RQT2 + PXY
42	ENDDO
43	ENDIF
44	C
45	IF ((A+2).LE.NCHA) THEN
46	DO J=A+2,NCHA
47	PXY = DEXY(1,I,J)
48	IF (PXY.GE.EMIN) RQT2 = RQT2 + PXY
49	ENDDO
50	ENDIF
51	C
52	600 CONTINUE
53	C
54	B = IBAR(2,I)
55	C
56	C EM
57	C
58	IF (B.EQ.-1) GOTO 700
59	C
60	C END EM
61	C
62	IF (B.LE.2) B = 2
63	IF (B.GE.NCHA) B = NCHA - 1
64	DO J = B-1,B+1
65	IF (DEXY(2,I,J).GE.EMIN) RQT1 = RQT1 + DEXY(2,I,J)
66	ENDDO
67	C
68	IF ((B-2).GE.1) THEN
69	DO J=1,B-2
70	PXY = DEXY(2,I,J)
71	IF (PXY.GE.EMIN) RQT2 = RQT2 + PXY
72	ENDDO
73	ENDIF
74	C
75	IF ((B+2).LE.NCHA) THEN
76	DO J=B+2,NCHA
77	PXY = DEXY(2,I,J)
78	IF (PXY.GE.EMIN) RQT2 = RQT2 + PXY
79	ENDDO
80	ENDIF
81	C
82	700 CONTINUE
83	ENDDO
84	C
85	C
86	800 RETURN
87	END
88