Moved sxsrot to geant library to avoid gcbank in shared libraries for the
[u/mrichter/AliRoot.git] / TGeant3 / galicef.F
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
7 C
8 C *** DEFINITION OF ADDITIONAL PARTICLES AND DECAY MODES ***
9 C *** NVE 08-JUL-1994 SAP UTRECHT ***
10 C
11 C CALLED BY : UGINIT
12 C ORIGIN    : NICK VAN EIJNDHOVEN
13 C
14       DIMENSION BRATIO(6),MODE(6)
15 C
16 C --- 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)
26 C
27 C --- Define additional decay modes ---
28 C --- 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)
41 C --- 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)
58 C --- 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)
73 C --- 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)
88 C --- 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)
101 C --- 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)
114 C --- 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)
123 C --- 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)
132 C --- 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)
141 C --- 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)
152 C --- upsilon ---
153       IPA=114
154       CALL GSDK(IPA,BRATIO,MODE)
155 C --- phi ---
156       IPA=115
157       CALL GSDK(IPA,BRATIO,MODE)
158 C
159       RETURN
160       END
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)
181 10000      FORMAT(' *** GXDRAW ***:',
182      +            ' No more space to store MCVOL information.')
183             CALL GMAIL(0,0)
184             WRITE(CHMAIL, 10100)
185 10100      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)
191 10200       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