GEANT321/gphys/gannir.F

   1 *
   2 * $Id$
   3 *
   4 * $Log$
   5 * Revision 1.1.1.1  1995/10/24 10:21:21  cernlib
   6 * Geant
   7 *
   8 *
   9 #include "geant321/pilot.h"
  10 *CMZ :  3.21/02 29/03/94  15.41.21  by  S.Giani
  11 *-- Author :
  12       SUBROUTINE GANNIR
  13 C.
  14 C.    ******************************************************************
  15 C.    *                                                                *
  16 C.    *  Generates positron annihilation in rest                       *
  17 C.    *                                                                *
  18 C.    *    ==>Called by : GTELEC                                       *
  19 C.    *       Author    L.Urban, M.Maire *********                     *
  20 C.    *                                                                *
  21 C.    ******************************************************************
  22 C.
  23 #include "geant321/gcphys.inc"
  24 #include "geant321/gctrak.inc"
  25 #include "geant321/gccuts.inc"
  26 #include "geant321/gcking.inc"
  27 #include "geant321/gconsp.inc"
  28       DIMENSION RNDM(2)
  29 C.
  30 C.    ------------------------------------------------------------------
  31 C.
  32       KCASE = NAMEC(11)
  33       IF((IANNI.NE.1).OR.(EMASS.LE.CUTGAM)) THEN
  34           ISTOP = 2
  35           DESTEP = DESTEP + 2*EMASS
  36           RETURN
  37       ENDIF
  38          ISTOP = 1
  39          CALL GRNDM(RNDM,2)
  40          CT=-1.+2.*RNDM(1)
  41          ST=SQRT(ABS(1.-CT*CT))
  42          PHI=TWOPI*RNDM(2)
  43          NGKINE=NGKINE+1
  44          GKIN(1,NGKINE)=EMASS*ST*COS(PHI)
  45          GKIN(2,NGKINE)=EMASS*ST*SIN(PHI)
  46          GKIN(3,NGKINE)=EMASS*CT
  47          GKIN(4,NGKINE)=EMASS
  48          GKIN(5,NGKINE)=1.
  49          TOFD(NGKINE)=0.
  50          GPOS(1,NGKINE) = VECT(1)
  51          GPOS(2,NGKINE) = VECT(2)
  52          GPOS(3,NGKINE) = VECT(3)
  53          NGKINE=NGKINE+1
  54          DO 30 J=1,3
  55   30     GKIN(J,NGKINE)=-GKIN(J,NGKINE-1)
  56          GKIN(4,NGKINE)=EMASS
  57          TOFD(NGKINE)=0.
  58          GKIN(5,NGKINE)=1.
  59          GPOS(1,NGKINE) = VECT(1)
  60          GPOS(2,NGKINE) = VECT(2)
  61          GPOS(3,NGKINE) = VECT(3)
  62       END