--- DarthVader/CalorimeterLevel2/src/nuclei.for 2007/01/23 11:52:27 1.1 +++ DarthVader/CalorimeterLevel2/src/nuclei.for 2014/01/16 15:29:15 1.6 @@ -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 @@ -92,7 +101,7 @@ NGIR = NCL(M,J) ! number of clusters in layer J of the M View IF (NGIR.EQ.0) GOTO 5 C IF (NGIR.NE.1) GOTO 5 - IPP = 0. + IPP = 0 DO L = 1,NGIR CALL LASTRISCIA(CLUS(M,J,L),IN) RIL(L) = FLOAT(ABS(NN - IN)) @@ -157,7 +166,8 @@ QQQQ = QQQ / FLOAT(NPIA) ! Truncated mean . ENDIF C - 400 CONTINUE +c 400 CONTINUE + CONTINUE RETURN END