TEvtGen/EvtGenModels/jetsetcont.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: jetsetcont.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 jetsetcont(energy,flavor,ndaug,kf,px,py,pz,e)
  21
  22
  23
  24       implicit none
  25
  26       common/lujets/n,k(4000,5),p(4000,5),v(4000,5)
  27
  28       integer n,k
  29       real p,v
  30
  31       real*8 energy
  32
  33       real senergy
  34
  35       integer lucomp
  36       external lucomp
  37
  38       integer flavor,kf(100),i,ndaug
  39       real*8  px(100),py(100),pz(100),e(100)
  40
  41       senergy=energy
  42
  43
  44       call lu2ent(0,flavor,-flavor,senergy)
  45
  46
  47       ndaug=0
  48
  49       do 10,i=4,n
  50
  51 c     only wants the particles that are created from the
  52 c     fragmentation, not the secondary daugthers.
  53
  54         if (k(i,3).eq.3) then
  55           kf(i-3)=k(i,2)
  56           px(i-3)=p(i,1)
  57           py(i-3)=p(i,2)
  58           pz(i-3)=p(i,3)
  59           e(i-3)=p(i,4)
  60           ndaug=ndaug+1
  61          endif
  62
  63 10    continue
  64
  65
  66       end
  67
  68
  69