PHOS/shaker/sh2bod.f

   1 *CMZ :          17/07/98  15.49.05  by  Federico Carminati
   2 *-- Author :
   3 c
   4 c       Federico Antinori Productions is proud to present...
   5 c
   6 c                 S H A K E R
   7 c
   8 c       Central Rapidity Phase Space Cocktail Event Generator
   9 c
  10 c       =====================================================
  11
  12 c       version 0/05
  13
  14 c       09.12.91, FA: 0/03      Pre-release
  15 c       28.02.92, FA: 0/04      SH2BOD and SHMTSC
  16
  17
  18       SUBROUTINE SH2BOD(AMP,AMDA,AMDB,PP,PDA,PDB)
  19 c       ===========================================
  20
  21 c       Two-body decay of a parent of mass AMP and lab 4-momentum PP to
  22 c       a daughter of mass AMDA and one of mass AMDB.
  23 c       PDA and PDB contain the lab 4-momenta of the daughters
  24
  25       DIMENSION PP(4),PDA(4),PDB(4)
  26       DIMENSION PACM(4),PBCM(4)
  27
  28       IF ((AMDA+AMDB).GT.AMP) RETURN
  29       PACM(4) = (AMP**2+AMDA**2-AMDB**2)/2./AMP
  30       PBCM(4) = (AMP**2+AMDB**2-AMDA**2)/2./AMP
  31       PCM = SQRT(PACM(4)**2-AMDA**2)
  32       CT = 2.*RLU(0.)-1.
  33       ST = SQRT(1.-CT**2)
  34       PHI = 2.*3.14159*RLU(0.)
  35       PACM(1) = PCM*ST*COS(PHI)
  36       PACM(2) = PCM*ST*SIN(PHI)
  37       PACM(4) = PCM*CT
  38       PBCM(1) = -PACM(1)
  39       PBCM(2) = -PACM(2)
  40       PBCM(3) = -PACM(3)
  41       CALL LORENB(AMP,PP,PACM,PDA)
  42       CALL LORENB(AMP,PP,PBCM,PDB)
  43
  44       RETURN
  45       END