--- DarthVader/CalorimeterLevel2/src/laterale.for 2006/05/19 13:15:50 1.1.1.1 +++ DarthVader/CalorimeterLevel2/src/laterale.for 2009/08/04 14:01:20 1.5 @@ -9,8 +9,8 @@ INCLUDE 'INTEST.TXT' C INTEGER A,B, INPIA - INTEGER IBAR(2,NPLA) - REAL BAR(2,NPLA) + INTEGER IBAR(2,NPLAV) + REAL BAR(2,NPLAV) REAL RQT1, RQT2 COMMON / ANGOLO / BAR, IBAR SAVE / ANGOLO / @@ -29,22 +29,27 @@ C C END EM C - IF (A.LE.2) A = 3 - IF (A.GE.(NCHA-1)) A = NCHA - 2 + 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 - 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 + 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 +c print *,' I ',I,' A ',RQT1 600 CONTINUE C B = IBAR(2,I) @@ -55,26 +60,33 @@ C C END EM C - IF (B.LE.2) B = 3 - IF (B.GE.(NCHA-1)) B = NCHA - 2 + 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) +c print *,' ==> J ',j,' dexy ',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 + 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 +c print *,' I ',I,' B ',RQT1 700 CONTINUE ENDDO C C - 800 RETURN +c 800 RETURN + RETURN END