Bugfix in AliPoints2Memory

[u/mrichter/AliRoot.git] / GEANT321 / gtrak / gtckov.F

diff --git a/GEANT321/gtrak/gtckov.F b/GEANT321/gtrak/gtckov.F
index f7a7ecffa16e2abf6fc32516624b25638e44892b..3518a1b6f60162a7270ee827ab2b114214177684 100644 (file)
--- a/GEANT321/gtrak/gtckov.F
+++ b/GEANT321/gtrak/gtckov.F
@@ -104,7 +104,6 @@ C.    ------------------------------------------------------------------
 *            NMEC = 1
 *            LMEC(1)= 30
 *            GO TO 110
-
             GEKRAT=0.
             IEKBIN=1
          ELSE
@@ -131,7 +130,7 @@ C.    ------------------------------------------------------------------
 *
       IPROC  = 103
       STEP   = STEMAX
-*
+*     
 *  **   Step limitation due to in flight absorbtion ?
 *
       IF (ILABS.GT.0) THEN
@@ -329,12 +328,26 @@ C
                GEKRT1=1.-GEKRT2
                ABSCO2=Q(JABSC2+IEKBI2)*GEKRT1+Q(JABSC2+IEKBI2+1)*GEKRT2
                EFFIC2=Q(JEFFI2+IEKBI2)*GEKRT1+Q(JEFFI2+IEKBI2+1)*GEKRT2
+
                IF(JINDX2.GT.0) THEN
                  RIN2=Q(JINDX2+IEKBI2)*GEKRT1+Q(JINDX2+IEKBI2+1)*GEKRT2
                ELSE
                  RIN2=0.
                ENDIF
-               IPROC = 102
+               IF(ABS(RIN2-RIN1).LE.PREC) THEN
+*
+* *** The two media have the same refraction coefficient, i.e. nothing
+* *** has to be done
+*
+                  VECT(1)=VOUT(1)
+                  VECT(2)=VOUT(2)
+                  VECT(3)=VOUT(3)
+                  STEPLA = ABSCO2
+                  EFFIC  = EFFIC2
+                  IPROC=0
+               ELSE
+                  IPROC = 102
+               ENDIF
             ELSE
                ISTOP=2
                DESTEP=0
@@ -371,11 +384,13 @@ C
 *
       IF (IPROC.EQ.101) THEN
          ISTOP=2
+
          CALL GRNDM(RNDM,1)
          IF(RNDM.LT.EFFIC) THEN
 *
 * ***  Destep =/= 0 means that the photon has been detected
 *
+
             DESTEP=VECT(7)
          ELSE
             DESTEP=0.