--- DarthVader/CalorimeterLevel2/src/nuclei.for 2007/01/23 11:52:27 1.1 +++ DarthVader/CalorimeterLevel2/src/nuclei.for 2007/07/20 08:24:55 1.2 @@ -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 / @@ -69,9 +69,17 @@ C DO J = 1,NPLA 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