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 = 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	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 = 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