TEvtGen/EvtGenModels/begevtgenstore.F

   1 C--------------------------------------------------------------------------
   2 C
   3 C Environment:
   4 C      This software is part of the EvtGen package developed jointly
   5 C      for the BaBar and CLEO collaborations.  If you use all or part
   6 C      of it, please give an appropriate acknowledgement.
   7 C
   8 C Copyright Information: See EvtGen/COPYRIGHT
   9 C      Copyright (C) 1998      Caltech, UCSB
  10 C
  11 C Module: begevtgenstore.F
  12 C
  13 C Description:
  14 C
  15 C Modification history:
  16 C
  17 C    DJL/RYD     August 11, 1998         Module created
  18 C
  19 C------------------------------------------------------------------------
  20       subroutine begevtgenstore(nument,eventnum,numparticle,istat,
  21      +                          partnum,motherfirst,motherlast,
  22      +                          daugfirst,dauglast,
  23      +                          px,py,pz,e,m,x,y,z,t)
  24 *
  25 * routine to fill the stdhep common blocks from
  26 * evtgen (C++). This routine allows the C++ program not to
  27 * have to mess with common blocks.
  28 *
  29 * Anders Ryd,  Dec 96   Created.
  30 *
  31 *
  32       implicit none
  33 #include "EvtGenModels/stdhep.inc"
  34
  35         integer nument,eventnum,numparticle,istat,partnum
  36         integer motherfirst,motherlast,daugfirst,dauglast
  37
  38         integer stdhepnum
  39
  40         double precision px,py,pz,e,m,x,y,z,t
  41
  42         stdhepnum=partnum
  43
  44         if (istat.eq.1) then
  45
  46 c       make the quarks unstable so that geant will not try
  47 c       to put them through the detector simulation!
  48
  49          if (
  50      +        stdhepnum.eq.-1.or.
  51      +        stdhepnum.eq. 1.or.
  52      +        stdhepnum.eq.-2.or.
  53      +        stdhepnum.eq. 2.or.
  54      +        stdhepnum.eq.-3.or.
  55      +        stdhepnum.eq. 3.or.
  56      +        stdhepnum.eq.-4.or.
  57      +        stdhepnum.eq. 4.or.
  58      +        stdhepnum.eq.-5.or.
  59      +        stdhepnum.eq. 5.or.
  60      +        stdhepnum.eq.-6.or.
  61      +        stdhepnum.eq. 6) then
  62            istat=2
  63          endif
  64
  65
  66
  67
  68 c         if (
  69 c     +        stdhepnum.ne.-1.and.
  70 c     +        stdhepnum.ne. 1.and.
  71 c     +        stdhepnum.ne.-2.and.
  72 c     +        stdhepnum.ne. 2.and.
  73 c     +        stdhepnum.ne.-3.and.
  74 c     +        stdhepnum.ne. 3.and.
  75 c     +        stdhepnum.ne.-4.and.
  76 c     +        stdhepnum.ne. 4.and.
  77 c     +        stdhepnum.ne.-5.and.
  78 c     +        stdhepnum.ne. 5.and.
  79 c     +        stdhepnum.ne.-6.and.
  80 c     +        stdhepnum.ne. 6.and.
  81 c     +        stdhepnum.ne.-11.and.
  82 c     +        stdhepnum.ne.11.and.
  83 c     +        stdhepnum.ne.-13.and.
  84 c     +        stdhepnum.ne.13.and.
  85 c     +        stdhepnum.ne.12.and.
  86 c     +        stdhepnum.ne.-12.and.
  87 c     +        stdhepnum.ne.14.and.
  88 c     +        stdhepnum.ne.-14.and.
  89 c     +        stdhepnum.ne.16.and.
  90 c     +        stdhepnum.ne.-16.and.
  91 c     +        stdhepnum.ne.2212.and.
  92 c     +        stdhepnum.ne.-2212.and.
  93 c     +        stdhepnum.ne.2112.and.
  94 c     +        stdhepnum.ne.-2112.and.
  95 c     +        stdhepnum.ne.211.and.
  96 c     +        stdhepnum.ne.-211.and.
  97 c     +        stdhepnum.ne.321.and.
  98 c     +        stdhepnum.ne.-321.and.
  99 c     +        stdhepnum.ne.130.and.
 100 c     +        stdhepnum.ne.310.and.
 101 c     +        stdhepnum.ne.22) then
 102 c
 103 c            print *,'particle stdhep:',stdhepnum,' is stable'
 104 c
 105 c         endif
 106
 107         endif
 108
 109         nevhep=eventnum
 110         nhep=numparticle
 111         isthep(nument)=istat
 112         idhep(nument)=stdhepnum
 113         jmohep(1,nument)=motherfirst
 114         jmohep(2,nument)=motherlast
 115         jdahep(1,nument)=daugfirst
 116         jdahep(2,nument)=dauglast
 117         phep(1,nument)=px
 118         phep(2,nument)=py
 119         phep(3,nument)=pz
 120         phep(4,nument)=e
 121         phep(5,nument)=m
 122         vhep(1,nument)=x
 123         vhep(2,nument)=y
 124         vhep(3,nument)=z
 125         vhep(4,nument)=t
 126
 127          return
 128
 129       end
 130
 131
 132
 133
 134