GEANT321/fluka/xlamb.F

   1 *
   2 * $Id$
   3 *
   4 * $Log$
   5 * Revision 1.1.1.1  1995/10/24 10:20:03  cernlib
   6 * Geant
   7 *
   8 *
   9 #include "geant321/pilot.h"
  10 *CMZ :  3.21/02 29/03/94  15.41.44  by  S.Giani
  11 *-- Author :
  12 *$ CREATE XLAMB.FOR
  13 *COPY XLAMB
  14 *
  15 *=== xlamb ============================================================*
  16 *
  17       FUNCTION XLAMB(X,Y,Z)
  18
  19 #include "geant321/dblprc.inc"
  20 #include "geant321/dimpar.inc"
  21 #include "geant321/iounit.inc"
  22       COMMON /FKIDGB/ IDGB,IDG
  23       COMMON /FKGAMR/ REDU,AMO,AMM(15 )
  24       COMMON/FKDREI/TEST(12)
  25       COMMON/FKPRUN/ISYS
  26       YZ=Y-Z
  27       XLAMB=X*X-2.D0*X*(Y+Z)+YZ*YZ
  28       XLAM =XLAMB
  29       IF (IDGB.LE.0) GO TO 11
  30       IF(XLAMB.GT.1.D-12) GOTO 11
  31       WRITE(ISYS,12)
  32       WRITE(ISYS,10) XLAM,X,Y,Z,TEST
  33       WRITE(ISYS,13)
  34    12 FORMAT(/,10X,12H XLAMB PRINT)
  35    13 FORMAT(10X,60(1H*))
  36    10 FORMAT(4E20.8,5HXLAMB,/,12F10.5)
  37    11 CONTINUE
  38       IF(XLAMB.LE.0.D0)XLAMB=ABS(XLAMB)
  39       XLAMB=SQRT(XLAMB)
  40       RETURN
  41       END