]>
Commit | Line | Data |
---|---|---|
fe4da5cc | 1 | * |
2 | * $Id$ | |
3 | * | |
4 | * $Log$ | |
57b2d457 | 5 | * Revision 1.1.1.1 1999/05/18 15:55:20 fca |
6 | * AliRoot sources | |
7 | * | |
fe4da5cc | 8 | * Revision 1.1.1.1 1995/10/24 10:21:18 cernlib |
9 | * Geant | |
10 | * | |
11 | * | |
12 | #include "geant321/pilot.h" | |
13 | *CMZ : 3.21/04 09/01/95 16.20.45 by S.Ravndal | |
14 | *-- Author : | |
15 | SUBROUTINE GLUDKY | |
16 | C. | |
17 | C. ****************************************************************** | |
18 | C. * * | |
19 | C. * GEANT INTERFACE TO THE LUND MONTE-CARLO * | |
20 | C. * FOR TAU DECAY GENERATION ONLY * | |
21 | C. * * | |
22 | C. * THIS INTERFACE TO LUND * | |
23 | C. * A) INHIBITS DECAYS OF K0S,SIG+,SIG-,XI0,XI-,LAM,OM- * | |
24 | C. * B) ACCEPTS STABLE PARTICLES ONLY * | |
25 | C. * C) CONVERTS THE PARTICLE CODES TO GEANT CONVENTIONS * | |
26 | C. * * | |
27 | C. * THUS THE PARTICLES FORWARDED TO GEANT ARE * | |
28 | C. * GAMMA,ELECTRON,NEUTRINO(E,MU,TAU),MU-,PI+,K+,K0S,K0L, * | |
29 | C. * PROTON,NEUTRON,SIG(+,-),XI(0,-),LAM,OM- * | |
30 | C. * .... AND THEIR ANTIPARTICLES * | |
31 | C. * * | |
32 | C. * That routine can be called from the user routine GUDCAY * | |
33 | C. * * | |
34 | C. * Author R.Hemingway 26 SEPT 1985 * | |
35 | C. * * | |
36 | C. ****************************************************************** | |
37 | C. | |
38 | C | |
39 | #include "geant321/gcflag.inc" | |
40 | #include "geant321/gckine.inc" | |
41 | #include "geant321/gcking.inc" | |
42 | #include "geant321/gctrak.inc" | |
43 | #include "geant321/gconsp.inc" | |
44 | #include "geant321/lujets.inc" | |
45 | #include "geant321/ludat3.inc" | |
46 | PARAMETER (NSEL=44) | |
47 | DIMENSION IPLUND(NSEL) | |
48 | DATA IPLUND/ | |
49 | + 22, -11, 11, 12, -13, 13, 111, 211, | |
50 | + -211, 130, 321, -321, 2112, 2212, -2212, 310, | |
51 | + 221, 3122, 3222, 3212, 3112, 3322, 3312, 3334, | |
52 | + -2112, -3122, -3112, -3212, -3222, -3322, -3312, -3334, | |
53 | + -15, 15, 411, -411, 421, -421, 431, -431, | |
54 | + 4122, 24, -24, 23/ | |
55 | *. | |
56 | *. ------------------------------------------------------------------ | |
57 | *. | |
58 | IF (IFINIT(6).EQ.0)THEN | |
59 | * | |
60 | *-- Set some particles to be stable in LUND so that GEANT can handle | |
61 | *-- their tracking and decays. (K0s, sigma+-, xi0-, lambda0, omega-) | |
62 | * | |
57b2d457 | 63 | MDCY(PYCOMP(310) ,1)=0 |
64 | MDCY(PYCOMP(3222),1)=0 | |
65 | MDCY(PYCOMP(3112),1)=0 | |
66 | MDCY(PYCOMP(3322),1)=0 | |
67 | MDCY(PYCOMP(3312),1)=0 | |
68 | MDCY(PYCOMP(3122),1)=0 | |
69 | MDCY(PYCOMP(3334),1)=0 | |
fe4da5cc | 70 | IFINIT(6) = 1 |
71 | ENDIF | |
72 | C | |
73 | IF(IPART.EQ.33.OR.IPART.EQ.34) THEN | |
74 | N=0 | |
75 | EDKY=GETOT | |
76 | THDKY=ACOS(VECT(6)) | |
77 | PHDKY=ATAN2(VECT(5),VECT(4)) | |
78 | * | |
79 | * for TAU+ TAU- | |
80 | * could be extended to all charm particles | |
81 | * | |
57b2d457 | 82 | CALL PY1ENT(1,IPLUND(IPART),EDKY,THDKY,PHDKY) |
83 | CALL PYEXEC | |
fe4da5cc | 84 | NGKINE=0 |
85 | DO 20 I = 1,N | |
86 | IF(K(I,1).GE.1.AND.K(I,1).LE.10) THEN | |
87 | DO 10 IGE=1, NSEL | |
88 | IF(IPLUND(IGE).EQ.K(I,2)) THEN | |
89 | NGKINE=NGKINE+1 | |
90 | GKIN(1,NGKINE)=P(I,1) | |
91 | GKIN(2,NGKINE)=P(I,2) | |
92 | GKIN(3,NGKINE)=P(I,3) | |
93 | GKIN(4,NGKINE)=P(I,4) | |
94 | GKIN(5,NGKINE)=IGE | |
95 | GPOS(1,NGKINE)=V(I,1) | |
96 | GPOS(2,NGKINE)=V(I,2) | |
97 | GPOS(3,NGKINE)=V(I,3) | |
98 | TOFD(NGKINE) =V(I,4)/(10.*CLIGHT) | |
99 | ENDIF | |
100 | 10 CONTINUE | |
101 | ENDIF | |
102 | 20 CONTINUE | |
103 | ENDIF | |
104 | * | |
105 | END |