coverity

[u/mrichter/AliRoot.git] / HIJING / hijing1_36 / attrad.F

diff --git a/HIJING/hijing1_36/attrad.F b/HIJING/hijing1_36/attrad.F
index 6fdc61b41529e6746c844105aa929f1e4923e245..97af736d5ee81ff4f3a74af936a9878eebdeaba7 100644 (file)
--- a/HIJING/hijing1_36/attrad.F
+++ b/HIJING/hijing1_36/attrad.F
@@ -8,6 +8,7 @@ C       conduct soft radiation according to dipole approxiamtion
 C****************************************************************
        SUBROUTINE ATTRAD(IERROR)
 C
+#define BLANKET_SAVE
 #include "hiparnt.inc"
 #include "hijdat.inc"
 #include "lujets_hijing.inc"
@@ -18,6 +19,7 @@ C.....S INVARIANT MASS-SQUARED BETWEEN PARTONS I AND I+1......
 C.....SM IS THE LARGEST MASS-SQUARED....
 
 40     SM=0.
+       JL = 1
        DO 30 I=1,N-1
           S=2.*(P(I,4)*P(I+1,4)-P(I,1)*P(I+1,1)-P(I,2)*P(I+1,2)
      &         -P(I,3)*P(I+1,3))+P(I,5)**2+P(I+1,5)**2
@@ -62,7 +64,12 @@ C.....BOOST TO REST SYSTEM FOR PARTICLES JL AND JL+1.....
        CALL ATROBO(0.,0.,BEX,BEY,BEZ,IMIN,IMAX,IERROR)
        IF(IERROR.NE.0) RETURN
 C.....ROTATE TO Z-AXIS....
-       CTH=P(JL,3)/SQRT(P(JL,4)**2-P(JL,5)**2)
+        PABS=(P(JL,4)-P(JL,5))*(P(JL,4)+P(JL,5))
+        IF (PABS .GT. 0) THEN
+          CTH=P(JL,3)/SQRT(PABS)
+        ELSE
+          CTH=SIGN(1.,P(JL,3))
+        ENDIF
        IF(ABS(CTH).GT.1.0)  CTH=MAX(-1.,MIN(1.,CTH))
        THETA=ACOS(CTH)
        PHI=ULANGL_HIJING(P(JL,1),P(JL,2))