| 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) |