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 jetsetcdf(ip,m,ndaug,kf,km,px,py,pz,e)
23 C interface to JETSET 7.4 to have one particle decayed
24 C including possibly fragmentation, if the decay products include
32 COMMON/LYDAT1/MSTU(200),PARU(200),MSTJ(200),PARJ(200)
35 common/lyjets/n,k(4000,5),p(4000,5),v(4000,5)
46 integer ip,kf(ndmax),i,ndaug,km(ndmax)
47 integer kp,kid,ipart1,ipart
48 real*8 m,px(ndmax),py(ndmax),pz(ndmax),e(ndmax)
53 c used to use ly1ent both since it does not set the mass
54 c of the daughter we have to manipulate the common blocks our
55 c self. ryd April 25-1999
57 c call ly1ent(1,ip,0.0,0.0,0.0)
69 c now we can decay this particle
72 C code copied from LUEXEC to avoid error with shower
74 C...Decay products may develop a shower.
75 IF(MSTJ(92).GT.0) THEN
77 QMAX=SQRT(MAX(0.,(P(IP1,4)+P(IP1+1,4))**2-(P(IP1,1)+P(IP1+1,
78 & 1))**2-(P(IP1,2)+P(IP1+1,2))**2-(P(IP1,3)+P(IP1+1,3))**2))
79 CALL LYSHOW(IP1,IP1+1,QMAX)
82 ELSEIF(MSTJ(92).LT.0) THEN
84 CALL LYSHOW(IP1,-3,P(IP1,5))
97 c find partons, delete secondary partons, set mother pointers
106 if (abs(kid) .ge. 1 .and. abs(kid) .le. 8
108 2 .or. kid .ge. 91 .and. kid .le. 94) then
109 if (ipart1 .eq. 1) ipart1 = i
111 if (kp .ne. 1) goto 10
114 if (kp .gt. ipart) then
116 elseif (kp .ge. ipart1) then
131 c print '( 2I5,I12,4F12.4 )',ndaug,km(ndaug),kf(ndaug),
132 c 1 px(ndaug),py(ndaug),pz(ndaug),e(ndaug)