]>
Commit | Line | Data |
---|---|---|
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 | |
13 | C. | |
14 | C. ******************************************************************* | |
15 | C. * * | |
16 | C. * Define standard GEANT particles plus selected decay modes * | |
17 | C. * and branching ratios. * | |
18 | C. * * | |
19 | C. * All data taken from : M. AGUILAR-BENITEZ et al, * | |
20 | C. * Review of Particle Properties, * | |
21 | C. * Rev. Mod. Phys. 56(1984). * | |
22 | C. * * | |
23 | C. * Updated to: Phys.Rev. D50,1173(1994) * | |
24 | C. * S. Ravndal * | |
25 | C. * * | |
26 | C. * * | |
27 | C. * In the case of W and Z, the lifetimes are calculated from * | |
28 | C. * quoted upper limits on widths. * | |
29 | C. * * | |
30 | C. * NPAR Number of parent particles defined for decay. * | |
31 | C. * IPAR List of parent partilces allowed to decay. * | |
32 | C. * Currently set up for pi0,pi+,pi-,K0long,K+,K-, * | |
33 | C. * K0short,eta,lambda,sigma+,sigma0,sigma-,xi0, * | |
34 | C. * xi-,omega-,antilambda,antisigma -,antisigma 0, * | |
35 | C. * antisigma +,antixi 0,antixi +,antiomega + * | |
36 | C. * decays. * | |
37 | C. * MODE(I,J) I'th decay mode of J'th particle defined in * | |
38 | C. * IPAR. * | |
39 | C. * BRATIO(I,J) Branching ratio for I'th decay mode of J'th * | |
40 | C. * particle in IPAR. * | |
41 | C. * * | |
42 | C. * ==>Called by : <USER>, UGINIT * | |
43 | C. * Author G.Patrick ********* * | |
44 | C. * * | |
45 | C. ******************************************************************* | |
46 | C. | |
47 | #include "geant321/gcbank.inc" | |
48 | #include "geant321/gcphys.inc" | |
49 | #include "geant321/gconsp.inc" | |
50 | C | |
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) | |
55 | C | |
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)) | |
58 | C | |
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./ | |
111 | C. | |
112 | C. ------------------------------------------------------------------- | |
113 | C. | |
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) | |
167 | C | |
168 | C Define decay modes. | |
169 | C | |
170 | DO 10 I=1,NPAR | |
171 | CALL GSDK(IPAR(I),BRATIO(1,I),MODE(1,I)) | |
172 | 10 CONTINUE | |
173 | C | |
174 | 99 RETURN | |
175 | END |