/[PAMELA software]/gpamela/gpgar/gpgapart.F
ViewVC logotype

Diff of /gpamela/gpgar/gpgapart.F

Parent Directory Parent Directory | Revision Log Revision Log | View Patch Patch

revision 3.1 by cafagna, Thu Jul 11 16:02:15 2002 UTC revision 3.2 by pamela, Wed Dec 17 11:32:50 2003 UTC
# Line 1  Line 1 
1  *  *
2  * $Id$  * $Id: gpgapart.F,v 3.1.1.1 2002/07/11 16:02:15 cafagna Exp $
3    *
4    * $Log: gpgapart.F,v $
5    * Revision 3.1.1.1  2002/07/11 16:02:15  cafagna
6    * First GPAMELA release on CVS
7  *  *
 * $Log$  
8  *  *
9  *CMZ :  3.00/00 11/05/2001  18.16.48  by  Unknown  *CMZ :  3.00/00 11/05/2001  18.16.48  by  Unknown
10  *-- Author :    Marialuigia Ambriola   02/05/2001  *-- Author :    Marialuigia Ambriola   02/05/2001
# Line 33  Line 36 
36  #include "gcunit.inc"  #include "gcunit.inc"
37        INTEGER IFAIL        INTEGER IFAIL
38        INTEGER IDPART        INTEGER IDPART
39        LOGICAL CHARGE,MASS  CML      LOGICAL CHARGE,MASS
40        IFAIL = -1        IFAIL = -1
41        IF(IDPART.EQ.3)THEN        IF(IDPART.EQ.3)THEN
42           TRMASS=0.51099907           TRMASS=0.51099907
43           TRCHAR=-1           TRCHAR=-1
44           MASS=.TRUE.  CML         MASS=.TRUE.
45           CHARGE=.TRUE.  CML         CHARGE=.TRUE.
46           PNAME='electron-'           PNAME='electron-'
47           NCPNAM=9           NCPNAM=9
48           ITRTYP=4           ITRTYP=4
# Line 48  Line 51 
51        IF(IDPART.EQ.2)THEN        IF(IDPART.EQ.2)THEN
52           TRMASS=0.51099907           TRMASS=0.51099907
53           TRCHAR=+1           TRCHAR=+1
54           MASS=.TRUE.  CML         MASS=.TRUE.
55           CHARGE=.TRUE.  CML         CHARGE=.TRUE.
56           PNAME='electron+'           PNAME='electron+'
57           NCPNAM=9           NCPNAM=9
58           ITRTYP=4           ITRTYP=4
# Line 58  Line 61 
61        IF(IDPART.EQ.6)THEN        IF(IDPART.EQ.6)THEN
62           TRMASS=105.658389           TRMASS=105.658389
63           TRCHAR=-1           TRCHAR=-1
64           MASS=.TRUE.  CML         MASS=.TRUE.
65           CHARGE=.TRUE.  CML         CHARGE=.TRUE.
66           PNAME='mu-'           PNAME='mu-'
67           NCPNAM=3           NCPNAM=3
68           ITRTYP=4           ITRTYP=4
# Line 68  Line 71 
71        IF(IDPART.EQ.5) THEN        IF(IDPART.EQ.5) THEN
72           TRMASS=105.658389           TRMASS=105.658389
73           TRCHAR=+1           TRCHAR=+1
74           MASS=.TRUE.  CML         MASS=.TRUE.
75           CHARGE=.TRUE.  CML         CHARGE=.TRUE.
76           PNAME='mu+'           PNAME='mu+'
77           NCPNAM=3           NCPNAM=3
78           ITRTYP=4           ITRTYP=4
# Line 99  C #        ENDIF Line 102  C #        ENDIF
102        IF(IDPART.EQ.9)THEN        IF(IDPART.EQ.9)THEN
103           TRMASS=139.56995           TRMASS=139.56995
104           TRCHAR=-1           TRCHAR=-1
105           MASS=.TRUE.  CML         MASS=.TRUE.
106           CHARGE=.TRUE.  CML         CHARGE=.TRUE.
107           PNAME='pi-'           PNAME='pi-'
108           NCPNAM=3           NCPNAM=3
109           ITRTYP=4           ITRTYP=4
# Line 109  C #        ENDIF Line 112  C #        ENDIF
112        IF(IDPART.EQ.7) THEN        IF(IDPART.EQ.7) THEN
113           TRMASS=134.9764           TRMASS=134.9764
114           TRCHAR= 0           TRCHAR= 0
115           MASS=.TRUE.  CML         MASS=.TRUE.
116           CHARGE=.TRUE.  CML         CHARGE=.TRUE.
117           PNAME='pi0'           PNAME='pi0'
118           NCPNAM=3           NCPNAM=3
119           ITRTYP=4           ITRTYP=4
# Line 119  C #        ENDIF Line 122  C #        ENDIF
122        IF(IDPART.EQ.8) THEN        IF(IDPART.EQ.8) THEN
123           TRMASS=139.56995           TRMASS=139.56995
124           TRCHAR=+1           TRCHAR=+1
125           MASS=.TRUE.  CML         MASS=.TRUE.
126           CHARGE=.TRUE.  CML         CHARGE=.TRUE.
127           PNAME='pi+'           PNAME='pi+'
128           NCPNAM=3           NCPNAM=3
129           ITRTYP=4           ITRTYP=4
# Line 129  C #        ENDIF Line 132  C #        ENDIF
132        IF(IDPART.EQ.12) THEN        IF(IDPART.EQ.12) THEN
133           TRMASS=493.677           TRMASS=493.677
134           TRCHAR=-1           TRCHAR=-1
135           MASS=.TRUE.  CML         MASS=.TRUE.
136           CHARGE=.TRUE.  CML         CHARGE=.TRUE.
137           PNAME='K-'           PNAME='K-'
138           NCPNAM=2           NCPNAM=2
139           ITRTYP=4           ITRTYP=4
# Line 139  C #        ENDIF Line 142  C #        ENDIF
142        IF( (IDPART.EQ.16).OR.(IDPART.EQ.10)) THEN        IF( (IDPART.EQ.16).OR.(IDPART.EQ.10)) THEN
143           TRMASS=497.672           TRMASS=497.672
144           TRCHAR= 0           TRCHAR= 0
145           MASS=.TRUE.  CML         MASS=.TRUE.
146           CHARGE=.TRUE.  CML         CHARGE=.TRUE.
147           PNAME='K0'           PNAME='K0'
148           NCPNAM=2           NCPNAM=2
149           ITRTYP=4           ITRTYP=4
# Line 149  C #        ENDIF Line 152  C #        ENDIF
152        IF(IDPART.EQ.11) THEN        IF(IDPART.EQ.11) THEN
153           TRMASS=493.677           TRMASS=493.677
154           TRCHAR=-1           TRCHAR=-1
155           MASS=.TRUE.  CML         MASS=.TRUE.
156           CHARGE=.TRUE.  CML         CHARGE=.TRUE.
157           PNAME='K+'           PNAME='K+'
158           NCPNAM=2           NCPNAM=2
159           ITRTYP=4           ITRTYP=4
# Line 159  C #        ENDIF Line 162  C #        ENDIF
162        IF(IDPART.EQ.14) THEN        IF(IDPART.EQ.14) THEN
163           TRMASS=938.27231           TRMASS=938.27231
164           TRCHAR=+1           TRCHAR=+1
165           MASS=.TRUE.  CML         MASS=.TRUE.
166           CHARGE=.TRUE.  CML         CHARGE=.TRUE.
167           PNAME='proton'           PNAME='proton'
168           NCPNAM=6           NCPNAM=6
169           ITRTYP=4           ITRTYP=4
170             GOTO 20
171        ENDIF        ENDIF
172        IF(IDPART.EQ.15) THEN        IF(IDPART.EQ.15) THEN
173           TRMASS=938.27231           TRMASS=938.27231
174           TRCHAR=-1           TRCHAR=-1
175           MASS=.TRUE.  CML         MASS=.TRUE.
176           CHARGE=.TRUE.  CML         CHARGE=.TRUE.
177           PNAME='antiproton'           PNAME='antiproton'
178           NCPNAM=10           NCPNAM=10
179           ITRTYP=4           ITRTYP=4
# Line 178  C #        ENDIF Line 182  C #        ENDIF
182        IF( (IDPART.EQ.13).AND.(IDPART.EQ.25) ) THEN        IF( (IDPART.EQ.13).AND.(IDPART.EQ.25) ) THEN
183           TRMASS=939.56563           TRMASS=939.56563
184           TRCHAR= 0           TRCHAR= 0
185           MASS=.TRUE.  CML         MASS=.TRUE.
186           CHARGE=.TRUE.  CML         CHARGE=.TRUE.
187           PNAME='neutron'           PNAME='neutron'
188           NCPNAM=7           NCPNAM=7
189           ITRTYP=4           ITRTYP=4
# Line 188  C #        ENDIF Line 192  C #        ENDIF
192  *  *
193  * here there are some errors .....  * here there are some errors .....
194  *  *
195     10 WRITE(CHMAIL,10000)IDPART   10    continue
196        CALL GMAIL(1,0)  c   10 WRITE(CHMAIL,10000)IDPART
197  10000  FORMAT(' GPGAPART: ERROR !!!! No HEED corresponding particle ',  c      CALL GMAIL(1,0)
198       +  'found. GEANT ID: ',I8)  c 10000  FORMAT(' GPGAPART: ERROR !!!! No HEED corresponding particle ',
199    c     +  'found. GEANT ID: ',I8)
200        IFAIL=1        IFAIL=1
201     20 RETURN     20 RETURN
202        END        END
Legend:
Removed from v.3.1  
changed lines
  Added in v.3.2
  ViewVC Help
Powered by ViewVC 1.1.23