This commit was generated by cvs2svn to compensate for changes in r2,
[u/mrichter/AliRoot.git] / GEANT321 / gcons / gpart.F
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