--- DarthVader/CalorimeterLevel2/src/nuclei.for 2007/01/23 11:52:27 1.1 +++ DarthVader/CalorimeterLevel2/src/nuclei.for 2009/08/04 14:01:22 1.5 @@ -31,12 +31,12 @@ C REAL LUNGH REAL RPIANO(2) - REAL BAR(2,NPLA) - REAL NQQ(2*NPLA) + REAL BAR(2,NPLAV) + REAL NQQ(2*NPLAV) REAL RIL(NCHA/2) REAL REL(NCHA/2) REAL TG(2) - INTEGER IBAR(2,NPLA) + INTEGER IBAR(2,NPLAV) INTEGER NTOT(2),LH COMMON / ANGOLO / BAR, IBAR SAVE / ANGOLO / @@ -68,10 +68,19 @@ RPIANO(M) = 0. C DO J = 1,NPLA +C IF (MOD(J,2).EQ.0) THEN - SHIFT = -0.5 + IF (REVERSE.EQ.0) THEN + SHIFT = -0.5 + ELSE + SHIFT = +0.5 + ENDIF ELSE - SHIFT = 0.5 + IF (REVERSE.EQ.0) THEN + SHIFT = +0.5 + ELSE + SHIFT = -0.5 + ENDIF ENDIF C DO K = 1,NGIR @@ -157,7 +166,8 @@ QQQQ = QQQ / FLOAT(NPIA) ! Truncated mean . ENDIF C - 400 CONTINUE +c 400 CONTINUE + CONTINUE RETURN END