1 |
* |
* |
2 |
* $Id$ |
* $Id: gpuspe.F,v 3.1.1.1 2002/07/11 16:02:14 cafagna Exp $ |
3 |
|
* |
4 |
|
* $Log: gpuspe.F,v $ |
5 |
|
* Revision 3.1.1.1 2002/07/11 16:02:14 cafagna |
6 |
|
* First GPAMELA release on CVS |
7 |
* |
* |
|
* $Log$ |
|
8 |
* |
* |
9 |
*CMZ : 2.01/00 06/03/2000 13.07.03 by Francesco Cafagna |
*CMZ : 2.01/00 06/03/2000 13.07.03 by Francesco Cafagna |
10 |
*CMZ : 2.00/00 03/03/2000 15.39.06 by Francesco Cafagna |
*CMZ : 2.00/00 03/03/2000 15.39.06 by Francesco Cafagna |
41 |
#include "gphit.inc" |
#include "gphit.inc" |
42 |
#include "gcunit.inc" |
#include "gcunit.inc" |
43 |
* |
* |
44 |
INTEGER IACT,IS,ID,IPA,IT,NUMVOL(20),IMEC,IDHIT |
INTEGER IACT,IS,ID,IPA,IT,NUMVOL(20),NVOLSPE(3),IMEC,IDHIT |
45 |
REAL DELOSS,PATH,TRAPAR(7),VHIT(NHSPE) |
REAL DELOSS,PATH,TRAPAR(7),VHIT(NHSPE) |
46 |
REAL DELTA2,SIGMA2,SIGMA,GRNDM0,GPGAUS |
REAL DELTA2,SIGMA2,SIGMA,GRNDM0,GPGAUS |
47 |
LOGICAL SAVE |
LOGICAL SAVE |
71 |
VHIT(9) = IPA |
VHIT(9) = IPA |
72 |
VHIT(10)= TRAPAR(7) |
VHIT(10)= TRAPAR(7) |
73 |
IDHIT = -1 |
IDHIT = -1 |
74 |
|
NVOLSPE(1) = NUMVOL(1) |
75 |
|
IF(NUMVOL(1).EQ.0) NVOLSPE(1)=6 |
76 |
|
NVOLSPE(2) = NUMVOL(2) |
77 |
|
NVOLSPE(3) = NUMVOL(3) |
78 |
ENDIF |
ENDIF |
79 |
* |
* |
80 |
* Track is exiting current volume |
* Track is exiting current volume |
114 |
* Store the hit |
* Store the hit |
115 |
* |
* |
116 |
IF(IDHIT.EQ.-1) THEN |
IF(IDHIT.EQ.-1) THEN |
117 |
IF(SAVE) CALL GSAHIT (IS,ID,IT,NUMVOL,VHIT,IDHIT) |
c$$$ IF(SAVE) CALL GSAHIT (IS,ID,IT,NUMVOL,VHIT,IDHIT) |
118 |
|
IF(SAVE) CALL GSAHIT (IS,ID,IT,NVOLSPE,VHIT,IDHIT) |
119 |
IF (IDHIT.EQ.0) THEN |
IF (IDHIT.EQ.0) THEN |
120 |
WRITE(CHMAIL,10000) IS,ID,IT |
WRITE(CHMAIL,10000) IS,ID,IT |
121 |
CALL GMAIL(1,0) |
CALL GMAIL(1,0) |
167 |
* Store the hit |
* Store the hit |
168 |
* |
* |
169 |
IF(IDHIT.EQ.-1) THEN |
IF(IDHIT.EQ.-1) THEN |
170 |
IF(SAVE) CALL GSAHIT (IS,ID,IT,NUMVOL,VHIT,IDHIT) |
c$$$ IF(SAVE) CALL GSAHIT (IS,ID,IT,NUMVOL,VHIT,IDHIT) |
171 |
|
IF(SAVE) CALL GSAHIT (IS,ID,IT,NVOLSPE,VHIT,IDHIT) |
172 |
IF (IDHIT.EQ.0) THEN |
IF (IDHIT.EQ.0) THEN |
173 |
WRITE(CHMAIL,10200) IS,ID,IT |
WRITE(CHMAIL,10200) IS,ID,IT |
174 |
CALL GMAIL(1,0) |
CALL GMAIL(1,0) |