Better protection

diff --git a/GEANT321/peanut/bimsel.F b/GEANT321/peanut/bimsel.F
index ef77f5e9a760f1421ddc77ead68f654284567327..2e91084773af92c817e813809c4cfa4c84306f9b 100644 (file)
--- a/GEANT321/peanut/bimsel.F
+++ b/GEANT321/peanut/bimsel.F
@@ -2,6 +2,9 @@
 * $Id$
 *
 * $Log$
+* Revision 1.3  2002/06/20 15:58:37  hristov
+* Protection in case of invalid ACOS argument. The reason for such argument has to be investigated.
+*
 * Revision 1.2  2002/05/13 12:40:58  hristov
 * Dummy subroutines to avoid files with no code in
 *
@@ -290,26 +293,24 @@
             IF ( BIMPTR .GE. RADTOT ) THEN
                X1 = BIMPTR - RADTOT
                DUMMY = 2.D+00 * X1 / ( RADPRO + X1 )
-               IF (ABS(DUMMY).LE.1.0) THEN
-                  ANGRED = ACOS ( DUMMY ) / PI
-               ELSE
+               IF (ABS(DUMMY).GT.1.0) THEN
                   PRINT *,"Warning in GEANT321/peanut/bimsel.F "
                   PRINT *,"Illegal ACOS argument ",DUMMY
-                  ANGRED = SIGN(1.0, DUMMY )
+                  DUMMY = SIGN(1.0, DUMMY )
                ENDIF
+               ANGRED = ACOS ( DUMMY ) / PI
                X1 = X1 / ( R0PROT * RPRONU * RADCO2 )
                DSKRED = ( 0.5D+00 * X1 * X1 + X1 + 1.D+00 ) * EXP (-X1)
      &                * ANGRED
             ELSE
                X1 = RADPRO + BIMPTR - RADTOT
                DUMMY = 2.D+00 * X1 / ( RADPRO + X1 )
-               IF (ABS(DUMMY).LE.1.0) THEN
-                  ANGRED = ACOS ( DUMMY ) / PI
-               ELSE
+               IF (ABS(DUMMY).GT.1.0) THEN
                   PRINT *,"Warning in GEANT321/peanut/bimsel.F "
                   PRINT *,"Illegal ACOS argument ",DUMMY
-                  ANGRED = SIGN(1.0, DUMMY )
+                  DUMMY = SIGN(1.0, DUMMY )
                ENDIF
+               ANGRED = ACOS ( DUMMY ) / PI
                X1 = X1 / ( R0PROT * RPRONU * RADCO2 )
                DSKRED = 1.D+00 - ( 0.5D+00 * X1 * X1 + X1 + 1.D+00 )
      &                * EXP (-X1) * ANGRED