[u/mrichter/AliRoot.git] / GEANT321 / fluka / xlamb.F

*
* $Id$
*
* $Log$
* Revision 1.1.1.1  1995/10/24 10:20:03  cernlib
* Geant
*
*
#include "geant321/pilot.h"
*CMZ :  3.21/02 29/03/94  15.41.44  by  S.Giani
*-- Author :
*$ CREATE XLAMB.FOR
*COPY XLAMB
*
*=== xlamb ============================================================*
*
      FUNCTION XLAMB(X,Y,Z)
 
#include "geant321/dblprc.inc"
#include "geant321/dimpar.inc"
#include "geant321/iounit.inc"
      COMMON /FKIDGB/ IDGB,IDG
      COMMON /FKGAMR/ REDU,AMO,AMM(15 )
      COMMON/FKDREI/TEST(12)
      COMMON/FKPRUN/ISYS
      YZ=Y-Z
      XLAMB=X*X-2.D0*X*(Y+Z)+YZ*YZ
      XLAM =XLAMB
      IF (IDGB.LE.0) GO TO 11
      IF(XLAMB.GT.1.D-12) GOTO 11
      WRITE(ISYS,12)
      WRITE(ISYS,10) XLAM,X,Y,Z,TEST
      WRITE(ISYS,13)
   12 FORMAT(/,10X,12H XLAMB PRINT)
   13 FORMAT(10X,60(1H*))
   10 FORMAT(4E20.8,5HXLAMB,/,12F10.5)
   11 CONTINUE
      IF(XLAMB.LE.0.D0)XLAMB=ABS(XLAMB)
      XLAMB=SQRT(XLAMB)
      RETURN
      END
Commit	Line	Data
fe4da5cc	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=XX-2.D0X(Y+Z)+YZYZ
	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