Added SetUserDecay routine. When a particle decays the standard MC decay
[u/mrichter/AliRoot.git] / TGeant3 / galicef.F
CommitLineData
fe4da5cc 1*CMZ : 2.03/01 20/08/98 17.46.44 by Federico Carminati
2*CMZ : 2.00/01 20/04/98 10.51.01 by Unknown
3*CMZ : 2.00/00 16/04/98 12.54.34 by Unknown
4*CMZ : 1.05/00 08/06/94 09.19.45 by Nick van Eijndhoven (RUU/CERN)
5*-- Author : Nick van Eijndhoven (RUU/CERN) 08/06/94
6 SUBROUTINE SXPART
7C
8C *** DEFINITION OF ADDITIONAL PARTICLES AND DECAY MODES ***
9C *** NVE 08-JUL-1994 SAP UTRECHT ***
10C
11C CALLED BY : UGINIT
12C ORIGIN : NICK VAN EIJNDHOVEN
13C
14 DIMENSION BRATIO(6),MODE(6)
15C
16C --- Define additional particles
17 CALL GSPART(33,'OMEGA(783)',3,0.782, 0.,7.836E-23,0.,0)
18 CALL GSPART(34,'PHI(1020)' ,3,1.019, 0.,1.486E-22,0.,0)
19 CALL GSPART(35,'D +' ,4,1.870, 1.,1.066E-12,0.,0)
20 CALL GSPART(36,'D -' ,4,1.870,-1.,1.066E-12,0.,0)
21 CALL GSPART(37,'D 0' ,3,1.865, 0.,4.2E-13 ,0.,0)
22 CALL GSPART(38,'ANTI D 0' ,3,1.865, 0.,4.2E-13 ,0.,0)
23 CALL GSPART(42,'RHO +' ,4,0.768, 1.,4.353E-24,0.,0)
24 CALL GSPART(43,'RHO -' ,4,0.768,-1.,4.353E-24,0.,0)
25 CALL GSPART(44,'RHO 0' ,3,0.768, 0.,4.353E-24,0.,0)
26C
27C --- Define additional decay modes ---
28C --- omega(783) ---
29 DO KZ=1,6
30 BRATIO(KZ)=0
31 MODE(KZ)=0
32 ENDDO
33 IPA=33
34 BRATIO(1)=89.
35 BRATIO(2)=8.5
36 BRATIO(3)=2.5
37 MODE(1)=70809
38 MODE(2)=107
39 MODE(3)=908
40 CALL GSDK(IPA,BRATIO,MODE)
41C --- phi(1020) ---
42 DO KZ=1,6
43 BRATIO(KZ)=0
44 MODE(KZ)=0
45 ENDDO
46 IPA=34
47 BRATIO(1)=49.
48 BRATIO(2)=34.4
49 BRATIO(3)=12.9
50 BRATIO(4)=2.4
51 BRATIO(5)=1.3
52 MODE(1)=1112
53 MODE(2)=1610
54 MODE(3)=4407
55 MODE(4)=90807
56 MODE(5)=1701
57 CALL GSDK(IPA,BRATIO,MODE)
58C --- D+ ---
59 DO KZ=1,6
60 BRATIO(KZ)=0
61 MODE(KZ)=0
62 ENDDO
63 IPA=35
64 BRATIO(1)=25.
65 BRATIO(2)=25.
66 BRATIO(3)=25.
67 BRATIO(4)=25.
68 MODE(1)=80809
69 MODE(2)=120808
70 MODE(3)=111208
71 MODE(4)=110809
72 CALL GSDK(IPA,BRATIO,MODE)
73C --- D- ---
74 DO KZ=1,6
75 BRATIO(KZ)=0
76 MODE(KZ)=0
77 ENDDO
78 IPA=36
79 BRATIO(1)=25.
80 BRATIO(2)=25.
81 BRATIO(3)=25.
82 BRATIO(4)=25.
83 MODE(1)=90908
84 MODE(2)=110909
85 MODE(3)=121109
86 MODE(4)=120908
87 CALL GSDK(IPA,BRATIO,MODE)
88C --- D0 ---
89 DO KZ=1,6
90 BRATIO(KZ)=0
91 MODE(KZ)=0
92 ENDDO
93 IPA=37
94 BRATIO(1)=33.
95 BRATIO(2)=33.
96 BRATIO(3)=33.
97 MODE(1)=809
98 MODE(2)=1208
99 MODE(3)=1112
100 CALL GSDK(IPA,BRATIO,MODE)
101C --- Anti D0 ---
102 DO KZ=1,6
103 BRATIO(KZ)=0
104 MODE(KZ)=0
105 ENDDO
106 IPA=38
107 BRATIO(1)=33.
108 BRATIO(2)=33.
109 BRATIO(3)=33.
110 MODE(1)=809
111 MODE(2)=1109
112 MODE(3)=1112
113 CALL GSDK(IPA,BRATIO,MODE)
114C --- rho+ ---
115 DO KZ=1,6
116 BRATIO(KZ)=0
117 MODE(KZ)=0
118 ENDDO
119 IPA=42
120 BRATIO(1)=100.
121 MODE(1)=807
122 CALL GSDK(IPA,BRATIO,MODE)
123C --- rho- ---
124 DO KZ=1,6
125 BRATIO(KZ)=0
126 MODE(KZ)=0
127 ENDDO
128 IPA=43
129 BRATIO(1)=100.
130 MODE(1)=907
131 CALL GSDK(IPA,BRATIO,MODE)
132C --- rho0 ---
133 DO KZ=1,6
134 BRATIO(KZ)=0
135 MODE(KZ)=0
136 ENDDO
137 IPA=44
138 BRATIO(1)=100.
139 MODE(1)=707
140 CALL GSDK(IPA,BRATIO,MODE)
141C --- jpsi ---
142 DO KZ=1,6
143 BRATIO(KZ)=0
144 MODE(KZ)=0
145 ENDDO
146 IPA=113
147 BRATIO(1)=50.
148 BRATIO(2)=50.
149 MODE(1)=506
150 MODE(2)=605
151 CALL GSDK(IPA,BRATIO,MODE)
152C --- upsilon ---
153 IPA=114
154 CALL GSDK(IPA,BRATIO,MODE)
155C --- phi ---
156 IPA=115
157 CALL GSDK(IPA,BRATIO,MODE)
158C
159 RETURN
160 END
fe4da5cc 161*CMZ : 30/12/98 16.23.44 by Rene Brun
162*-- Author : Rene Brun 30/12/98
163 subroutine setclip(name,xmin,xmax,ymin,ymax,zmin,zmax)
164*
165* define the clip box parameters for volume name
166*
167#undef CERNLIB_GEANT321_GCUNIT_INC
168#include "geant321/gcunit.inc"
169#undef CERNLIB_GEANT321_GCMUTR_INC
170#include "geant321/gcmutr.inc"
171#undef CERNLIB_GEANT321_GCGOBJ_INC
172#include "geant321/gcgobj.inc"
173*KEND.
174 character *(*) name
175 real xmin,xmax,ymin,ymax,zmin,zmax
176*
177 IHOLE=0
178 NCVOLS=NCVOLS+1
179 IF(NCVOLS.EQ.MULTRA)THEN
180 WRITE(CHMAIL, 10000)
18110000 FORMAT(' *** GXDRAW ***:',
182 + ' No more space to store MCVOL information.')
183 CALL GMAIL(0,0)
184 WRITE(CHMAIL, 10100)
18510100 FORMAT(' *** GXDRAW ***: Please reset MCVOL')
186 CALL GMAIL(0,0)
187 return
188 ENDIF
189 IF(XMIN.GE.XMAX.OR.YMIN.GE.YMAX.OR.ZMIN.GE.ZMAX)THEN
190 WRITE(CHMAIL,10200)
19110200 FORMAT(' Wrong Box limits. Check values ')
192 CALL GMAIL(0,0)
193 return
194 ENDIF
195****SG
196 GNNVV(NCVOLS)=NAME
197 GNASH(NCVOLS)='BOX'
198 GXMIN(NCVOLS)=XMIN
199 GXMAX(NCVOLS)=XMAX
200 GYMIN(NCVOLS)=YMIN
201 GYMAX(NCVOLS)=YMAX
202 GZMIN(NCVOLS)=ZMIN
203 GZMAX(NCVOLS)=ZMAX
204 IF(GXMIN(NCVOLS).GT.-99999.)IHOLE=1
205* Resetting Mcvol mode
206 IF(GNNVV(NCVOLS).EQ.'.')THEN
207 IHOLE=0
208 DO 10 JJ=1,NCVOLS
209 GNNVV(JJ)=' '
210 GXMIN(JJ)=-100000
211 GXMAX(JJ)=-99999
212 GYMIN(JJ)=-100000
213 GYMAX(JJ)=-99999
214 GZMIN(JJ)=-100000
215 GZMAX(JJ)=-99999
216 10 CONTINUE
217 NCVOLS=0
218 ENDIF
219 end
220
221 subroutine setbomb(boom)
222*
223* set the bomb factor
224*
225#undef CERNLIB_GEANT321_GCMUTR_INC
226#include "geant321/gcmutr.inc"
227*KEND.
228*
229 real boom
230*
231 gboom = boom
232 end