Bugfix in AliPoints2Memory
[u/mrichter/AliRoot.git] / GEANT321 / gcons / gpart.F
CommitLineData
fe4da5cc 1*
2* $Id$
3*
4* $Log$
5* Revision 1.1.1.1 1995/10/24 10:20:15 cernlib
6* Geant
7*
8*
9#include "geant321/pilot.h"
10*CMZ : 3.21/03 14/09/94 14.33.42 by S.Ravndal
11*-- Author :
12 SUBROUTINE GPART
13C.
14C. *******************************************************************
15C. * *
16C. * Define standard GEANT particles plus selected decay modes *
17C. * and branching ratios. *
18C. * *
19C. * All data taken from : M. AGUILAR-BENITEZ et al, *
20C. * Review of Particle Properties, *
21C. * Rev. Mod. Phys. 56(1984). *
22C. * *
23C. * Updated to: Phys.Rev. D50,1173(1994) *
24C. * S. Ravndal *
25C. * *
26C. * *
27C. * In the case of W and Z, the lifetimes are calculated from *
28C. * quoted upper limits on widths. *
29C. * *
30C. * NPAR Number of parent particles defined for decay. *
31C. * IPAR List of parent partilces allowed to decay. *
32C. * Currently set up for pi0,pi+,pi-,K0long,K+,K-, *
33C. * K0short,eta,lambda,sigma+,sigma0,sigma-,xi0, *
34C. * xi-,omega-,antilambda,antisigma -,antisigma 0, *
35C. * antisigma +,antixi 0,antixi +,antiomega + *
36C. * decays. *
37C. * MODE(I,J) I'th decay mode of J'th particle defined in *
38C. * IPAR. *
39C. * BRATIO(I,J) Branching ratio for I'th decay mode of J'th *
40C. * particle in IPAR. *
41C. * *
42C. * ==>Called by : <USER>, UGINIT *
43C. * Author G.Patrick ********* *
44C. * *
45C. *******************************************************************
46C.
47#include "geant321/gcbank.inc"
48#include "geant321/gcphys.inc"
49#include "geant321/gconsp.inc"
50C
51 DIMENSION BRATIO(6,24),BRTIO1(6,20),BRTIO2(6,4)
52 DIMENSION IPAR(24)
53 DIMENSION MODE(6,24),MODE1(6,20),MODE2(6,4)
54 DIMENSION UBUF(1)
55C
56 EQUIVALENCE(BRTIO1(1,1),BRATIO(1,1)),(BRTIO2(1,1),BRATIO(1,21))
57 EQUIVALENCE(MODE1(1,1),MODE(1,1)),(MODE2(1,1),MODE(1,21))
58C
59 DATA IPAR/ 7, 8, 9,10,11,12,16,17,18,19,20,21,22,23,24,26,27,
60 * 28,29,30,31,32,5,6/
61 DATA BRTIO1/ 98.802, 1.198, 0., 0., 0., 0.,
62 * 100.00, 0., 0., 0., 0., 0.,
63 * 100.00, 0., 0., 0., 0., 0.,
64 * 21.60, 19.35, 19.35, 13.50, 13.50, 12.38,
65 * 63.51, 21.17, 5.59, 4.82, 3.18, 1.73,
66 * 63.51, 21.17, 5.59, 4.82, 3.18, 1.73,
67 * 68.61, 31.39, 0., 0., 0., 0.,
68 * 38.80, 31.90, 23.6, 4.88, 0.5, 0.071,
69 * 63.90, 35.80, 0., 0., 0., 0.,
70 * 51.57, 48.30, 0., 0., 0., 0.,
71 * 100.00, 0., 0., 0., 0., 0.,
72 * 100.00, 0., 0., 0., 0., 0.,
73 * 100.00, 0., 0., 0., 0., 0.,
74 * 100.00, 0., 0., 0., 0., 0.,
75 * 67.80, 23.60, 8.60, 0., 0., 0.,
76 * 63.90, 35.80, 0., 0., 0., 0.,
77 * 51.57, 48.30, 0., 0., 0., 0.,
78 * 100.00, 0., 0., 0., 0., 0.,
79 * 100.00, 0., 0., 0., 0., 0.,
80 * 100.00, 0., 0., 0., 0., 0./
81 DATA BRTIO2/100.00, 0., 0., 0., 0., 0.,
82 * 67.80, 23.60, 8.60, 0., 0., 0.,
83 * 100.00, 0., 0., 0., 0., 0.,
84 * 100.00, 0., 0., 0., 0., 0./
85 DATA MODE1/ 101,30201, 0, 0, 0, 0,
86 * 405, 0, 0, 0, 0, 0,
87 * 406, 0, 0, 0, 0, 0,
88 * 70707,40308,40209,40608,40509,70908,
89 * 405, 708,90808,70402,70405,70708,
90 * 406, 709,80909,70403,70406,70709,
91 * 908, 707, 0, 0, 0, 0,
92 * 101,70707,70908,10908,10302,10107,
93 * 914, 713, 0, 0, 0, 0,
94 * 714, 813, 0, 0, 0, 0,
95 * 118, 0, 0, 0, 0, 0,
96 * 913, 0, 0, 0, 0, 0,
97 * 718, 0, 0, 0, 0, 0,
98 * 918, 0, 0, 0, 0, 0,
99 * 1218, 922, 723, 0, 0, 0,
100 * 815, 725, 0, 0, 0, 0,
101 * 715, 925, 0, 0, 0, 0,
102 * 126, 0, 0, 0, 0, 0,
103 * 825, 0, 0, 0, 0, 0,
104 * 726, 0, 0, 0, 0, 0/
105 DATA MODE2/ 826, 0, 0, 0, 0, 0,
106 * 1126, 830, 731, 0, 0, 0,
107 * 40402, 0, 0, 0, 0, 0,
108 * 40403, 0, 0, 0, 0, 0/
109 DATA NPAR/24/
110 DATA UBUF/0./
111C.
112C. -------------------------------------------------------------------
113C.
114 REMASS=EMASS
115 REMMU =EMMU
116 RPMASS=PMASS
117 CALL GSPART( 1,'GAMMA$ ',1,0. , 0.,1.E+15,UBUF,0)
118 CALL GSPART( 2,'POSITRON$ ',2,REMASS , 1.,1.E+15,UBUF,0)
119 CALL GSPART( 3,'ELECTRON$ ',2,REMASS ,-1.,1.E+15,UBUF,0)
120 CALL GSPART( 4,'NEUTRINO$ ',3,0. , 0.,1.E+15,UBUF,0)
121 CALL GSPART( 5,'MUON +$ ',5,REMMU , 1.,2.19703E-06,UBUF,0)
122 CALL GSPART( 6,'MUON -$ ',5,REMMU ,-1.,2.19703E-06,UBUF,0)
123 CALL GSPART( 7,'PION 0$ ',3,0.1349764,0.,0.84E-16,UBUF,0)
124 CALL GSPART( 8,'PION +$ ',4,0.1395700, 1.,2.603E-08,UBUF,0)
125 CALL GSPART( 9,'PION -$ ',4,0.1395700,-1.,2.603E-08,UBUF,0)
126 CALL GSPART(10,'KAON 0 LONG$',3,0.497672, 0.,5.17E-08,UBUF,0)
127 CALL GSPART(11,'KAON +$ ',4,0.493677, 1.,1.237E-08,UBUF,0)
128 CALL GSPART(12,'KAON -$ ',4,0.493677,-1.,1.237E-08,UBUF,0)
129 CALL GSPART(13,'NEUTRON$ ',3,0.93956563, 0.,8.870E+02,UBUF,0)
130 CALL GSPART(14,'PROTON$ ',4,RPMASS , 1.,1.000000E+15,UBUF,0)
131 CALL GSPART(15,'ANTIPROTON$ ',4,RPMASS ,-1.,1.000000E+15,UBUF,0)
132 CALL GSPART(16,'KAON 0 SHORT$',3,0.497672, 0.,8.926E-11,UBUF,0)
133 CALL GSPART(17,'ETA$ ',3,0.54745 , 0.,5.485E-19,UBUF,0)
134 CALL GSPART(18,'LAMBDA$ ',3,1.115684, 0.,2.632E-10,UBUF,0)
135 CALL GSPART(19,'SIGMA +$ ',4,1.18937 , 1.,0.799E-10,UBUF,0)
136 CALL GSPART(20,'SIGMA 0$ ',3,1.19255 , 0.,7.4E-20,UBUF,0)
137 CALL GSPART(21,'SIGMA -$ ',4,1.197436,-1.,1.479E-10,UBUF,0)
138 CALL GSPART(22,'XI 0$ ',3,1.3149 , 0.,2.9E-10,UBUF,0)
139 CALL GSPART(23,'XI -$ ',4,1.32132 ,-1.,1.639E-10,UBUF,0)
140 CALL GSPART(24,'OMEGA -$ ',4,1.67245 ,-1.,0.822E-10,UBUF,0)
141 CALL GSPART(25,'ANTINEUTRON$',3,0.93956563, 0.,8.870E+02,UBUF,0)
142 CALL GSPART(26,'ANTILAMBDA$ ',3,1.115684, 0.,2.632E-10,UBUF,0)
143 CALL GSPART(27,'ANTISIGMA -$',4,1.18937 ,-1.,0.799E-10,UBUF,0)
144 CALL GSPART(28,'ANTISIGMA 0$',3,1.19255 , 0.,7.4E-20,UBUF,0)
145 CALL GSPART(29,'ANTISIGMA +$',4,1.197436, 1.,1.479E-10,UBUF,0)
146 CALL GSPART(30,'ANTIXI 0$ ',3,1.3149 , 0.,2.9E-10,UBUF,0)
147 CALL GSPART(31,'ANTIXI +$ ',4,1.32132 , 1.,1.639E-10,UBUF,0)
148 CALL GSPART(32,'ANTIOMEGA +$',4,1.67245 , 1.,0.822E-10,UBUF,0)
149* CALL GSPART(33,'TAU +$ ',4,1.7771 , 1.,3.040000E-13,UBUF,0)
150* CALL GSPART(34,'TAU -$ ',4,1.7771 ,-1.,3.040000E-13,UBUF,0)
151* CALL GSPART(35,'D +$ ',4,1.8694 , 1.,1.062000E-12,UBUF,0)
152* CALL GSPART(36,'D -$ ',4,1.8694 ,-1.,1.062000E-12,UBUF,0)
153* CALL GSPART(37,'D 0$ ',3,1.8646 , 0.,4.280000E-13,UBUF,0)
154* CALL GSPART(38,'ANTI D 0$ ',3,1.8646 , 0.,4.280000E-13,UBUF,0)
155* CALL GSPART(39,'DS+$ ',4,1.9685 , 1.,4.360000E-13,UBUF,0)
156* CALL GSPART(40,'DS-$ ',4,1.9685 ,-1.,4.360000E-13,UBUF,0)
157* CALL GSPART(41,'LAMBDA C +$ ',4,2.2851 , 1.,1.790000E-13,UBUF,0)
158* CALL GSPART(42,'W +$ ',4,80.220 , 1.,3.160000E-25,UBUF,0)
159* CALL GSPART(43,'W -$ ',4,80.220 ,-1.,2.640000E-25,UBUF,0)
160* CALL GSPART(44,'Z 0$ ',3,91.187 , 0.,7.740000E-26,UBUF,0)
161 CALL GSPART(45,'DEUTERON$ ',8,1.875613,+1.,1.000000E+15,UBUF,0)
162 CALL GSPART(46,'TRITON$ ',8,2.80925 ,+1.,1.000000E+15,UBUF,0)
163 CALL GSPART(47,'ALPHA$ ',8,3.727417,+2.,1.000000E+15,UBUF,0)
164 CALL GSPART(48,'GEANTINO$ ',6,0. , 0.,1.000000E+15,UBUF,0)
165 CALL GSPART(49,'HE3$ ',8,2.80923 ,+2.,1.000000E+15,UBUF,0)
166 CALL GSPART(50,'Cerenkov$ ',7,0. , 0.,1.000000E+15,UBUF,0)
167C
168C Define decay modes.
169C
170 DO 10 I=1,NPAR
171 CALL GSDK(IPAR(I),BRATIO(1,I),MODE(1,I))
172 10 CONTINUE
173C
174 99 RETURN
175 END