1 C--------------------------------------------------------------------------
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.
8 C Copyright Information: See EvtGen/COPYRIGHT
9 C Copyright (C) 1998 Caltech, UCSB
15 C Modification history:
17 C DJL/RYD August 11, 1998 Module created
19 C------------------------------------------------------------------------
20 subroutine continuum(energy,flavor,ndaug,kf,px,py,pz,e)
24 common/lujets/n,k(4000,5),p(4000,5),v(4000,5)
29 real*8 pxsum,pysum,pzsum
38 integer flavor,kf(100),i,ndaug,istart
39 real*8 px(100),py(100),pz(100),e(100)
43 C if(flavor.ne.0) kf_user=flavor
47 C call lueevt(kf_user,senergy) ! e+e- --> qq shower parton model
48 4 call lueevt(flavor,senergy) ! e+e- --> qq shower parton model
54 C--- Save the primary quarks
57 if(abs(k(i,2)).gt.6) goto 1
58 if(k(i,2).eq.0) goto 1
59 if(k(i,3).ne.0) goto 1
68 C sum to check that we preserve momentum
76 if(k(i,3).le.0.and.k(i,2).ne.22) goto 2
77 c if particle is jet or cluster
79 if (k(i,2).eq.22.and.istart.eq.0) goto 22
80 if ( k(i,3).le.0) goto 2
82 if(k(k(i,3),2).eq.92.or.k(k(i,3),2).eq.91) goto 22
87 C Copy particle. If parent is cluster or string or photon if
88 C before first cluster or string. Uggly...
89 if (k(i,2).ne.22) istart=1
102 if (abs(pxsum).gt.0.001.or.
103 + abs(pysum).gt.0.001.or.
104 + abs(pzsum).gt.0.001) then
107 print *, 'Momentum not conserved in jetset fragmentation:'
108 print *,'dPx:',pxsum,' dPy:',pysum,' dPz:',pzsum