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
11 C Module: begevtgenstore.F
15 C Modification history:
17 C DJL/RYD August 11, 1998 Module created
19 C------------------------------------------------------------------------
20 subroutine begevtgenstore(nument,eventnum,numparticle,istat,
21 + partnum,motherfirst,motherlast,
23 + px,py,pz,e,m,x,y,z,t)
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.
29 * Anders Ryd, Dec 96 Created.
33 #include "EvtGenModels/stdhep.inc"
35 integer nument,eventnum,numparticle,istat,partnum
36 integer motherfirst,motherlast,daugfirst,dauglast
40 double precision px,py,pz,e,m,x,y,z,t
46 c make the quarks unstable so that geant will not try
47 c to put them through the detector simulation!
61 + stdhepnum.eq. 6) then
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
103 c print *,'particle stdhep:',stdhepnum,' is stable'
112 idhep(nument)=stdhepnum
113 jmohep(1,nument)=motherfirst
114 jmohep(2,nument)=motherlast
115 jdahep(1,nument)=daugfirst
116 jdahep(2,nument)=dauglast