GEANT321/fluka/twopad.F

   1 *
   2 * $Id$
   3 *
   4 * $Log$
   5 * Revision 1.1.1.1  1995/10/24 10:20:02  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 TWOPAD.FOR
  13 *COPY TWOPAD
  14 *
  15 *=== twopad ===========================================================*
  16 *
  17       SUBROUTINE TWOPAD(UMO,ECM1,ECM2,PCM1,PCM2,COD1,COF1,SIF1,
  18      *COD2,COF2,SIF2,AM1,AM2)
  19
  20 #include "geant321/dblprc.inc"
  21 #include "geant321/dimpar.inc"
  22 #include "geant321/iounit.inc"
  23       REAL RNDM(1)
  24 C*****TWO PARTICLE DECAY IN THE CM - SYSTEM
  25       COMMON /FKGAMR/ REDU,AMO,AMM(15 )
  26       UMO2=UMO*UMO
  27       AM11=AM1*AM1
  28       AM22=AM2*AM2
  29       ECM1=(UMO2+AM11-AM22)/(2.D0*UMO)
  30       ECM2=UMO-ECM1
  31       WAU=ECM1**2-AM11
  32       IF (WAU.LT.0.D0) REDU=REDU-5.D0
  33       WAU=ABS(WAU)
  34       PCM1=SQRT(WAU)
  35       PCM2=PCM1
  36       CALL COSI(SIF1,COF1)
  37       CALL GRNDM(RNDM,1)
  38       COD1=RNDM(1)
  39       COD1=2.D0*COD1-1.D0
  40       COD2=-COD1
  41       COF2=-COF1
  42       SIF2=-SIF1
  43       RETURN
  44       END