]>
Commit | Line | Data |
---|---|---|
da0e9ce3 | 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 |