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