| 1 |
* |
* |
| 2 |
* $Id$ |
* $Id: gputof.F,v 3.1.1.1 2002/07/11 16:02:01 cafagna Exp $ |
| 3 |
|
* |
| 4 |
|
* $Log: gputof.F,v $ |
| 5 |
|
* Revision 3.1.1.1 2002/07/11 16:02:01 cafagna |
| 6 |
|
* First GPAMELA release on CVS |
| 7 |
* |
* |
|
* $Log$ |
|
| 8 |
* |
* |
| 9 |
*CMZ : 2.01/00 05/04/2000 14.35.18 by Marialuigia Ambriola |
*CMZ : 2.01/00 05/04/2000 14.35.18 by Marialuigia Ambriola |
| 10 |
*CMZ : 2.00/00 25/02/2000 19.34.00 by Francesco Cafagna |
*CMZ : 2.00/00 25/02/2000 19.34.00 by Francesco Cafagna |
| 50 |
IF(IACT.EQ.0) THEN |
IF(IACT.EQ.0) THEN |
| 51 |
IF(IMEC.GT.0) THEN |
IF(IMEC.GT.0) THEN |
| 52 |
SAVE = .TRUE. |
SAVE = .TRUE. |
| 53 |
*ML: |
VHIT(7) = DELOSS + VHIT(7) |
|
* VHIT(7) = DELOSS + VHIT(7) |
|
|
VHIT(11) = DELOSS + VHIT(11) |
|
|
*END ML. |
|
| 54 |
ENDIF |
ENDIF |
| 55 |
ENDIF |
ENDIF |
| 56 |
* |
* |
| 64 |
ENDDO |
ENDDO |
| 65 |
IF(IMEC.GT.0) THEN |
IF(IMEC.GT.0) THEN |
| 66 |
SAVE = .TRUE. |
SAVE = .TRUE. |
| 67 |
*ML: |
VHIT(7) = DELOSS |
|
VHIT(11) = DELOSS |
|
|
* VHIT(7) = DELOSS |
|
|
*END ML. |
|
| 68 |
ENDIF |
ENDIF |
| 69 |
VHIT(8) = TIME |
VHIT(8) = TIME |
| 70 |
VHIT(9) = PATH |
VHIT(9) = PATH |
| 71 |
VHIT(10) = IPA |
VHIT(10) = IPA |
| 72 |
*ML: |
VHIT(11) = TRAPAR(7) |
|
* VHIT(11) = TRAPAR(7) |
|
|
VHIT(7) = TRAPAR(7) |
|
|
*END ML. |
|
| 73 |
IDHIT = -1 |
IDHIT = -1 |
| 74 |
ENDIF |
ENDIF |
| 75 |
* |
* |
| 85 |
ENDDO |
ENDDO |
| 86 |
IF(IMEC.GT.0) THEN |
IF(IMEC.GT.0) THEN |
| 87 |
SAVE = .TRUE. |
SAVE = .TRUE. |
| 88 |
*ML: |
VHIT(7) = DELOSS +VHIT(7) |
|
* VHIT(7) = DELOSS + VHIT(7) |
|
|
VHIT(11) = DELOSS +VHIT(11) |
|
|
*END ML. |
|
| 89 |
ENDIF |
ENDIF |
| 90 |
VHIT(8) = (TIME + VHIT(8))/2. |
VHIT(8) = (TIME + VHIT(8))/2. |
| 91 |
VHIT(9) = PATH - VHIT(9) |
VHIT(9) = PATH - VHIT(9) |
| 92 |
*ML: |
c$$$ VHIT(11) = (TRAPAR(7) + VHIT(11))/2. |
|
* VHIT(11) = (TRAPAR(7) + VHIT(11))/2. |
|
|
VHIT(7) = (TRAPAR(7) + VHIT(7))/2. |
|
|
*END ML. |
|
| 93 |
* |
* |
| 94 |
* Store the hit |
* Store the hit |
| 95 |
* |
* |
| 96 |
IF(IDHIT.EQ.-1) THEN |
IF(IDHIT.EQ.-1) THEN |
|
*ML: |
|
|
* IF (SAVE) CALL GSAHIT (IS,ID,IT,NUMVOL,VHIT,IDHIT) |
|
|
* IF (IDHIT.EQ.0) THEN |
|
|
* WRITE(CHMAIL,10000) IS,ID,IT |
|
|
* CALL GMAIL(1,0) |
|
|
* ENDIF |
|
|
* ELSE |
|
|
* WRITE(CHMAIL,10100) IS,ID,IT |
|
|
* CALL GMAIL(1,0) |
|
|
* ENDIF |
|
| 97 |
IF (SAVE) THEN |
IF (SAVE) THEN |
| 98 |
IF(IS.EQ.ISS4) THEN |
IF(IS.EQ.ISS4) THEN |
| 99 |
CALL GSCHIT(IS,ID,1,NUMVOL,VHIT,1,IDHIT) |
CALL GSCHIT(IS,ID,1,NUMVOL,VHIT,1,IDHIT) |
| 109 |
WRITE(CHMAIL,10100) IS,ID,IT |
WRITE(CHMAIL,10100) IS,ID,IT |
| 110 |
CALL GMAIL(1,0) |
CALL GMAIL(1,0) |
| 111 |
ENDIF |
ENDIF |
|
* END ML. |
|
| 112 |
ENDIF |
ENDIF |
| 113 |
* |
* |
| 114 |
* Track is exiting the set up |
* Track is exiting the set up |
| 127 |
ENDDO |
ENDDO |
| 128 |
IF(IMEC.GT.0) THEN |
IF(IMEC.GT.0) THEN |
| 129 |
SAVE = .TRUE. |
SAVE = .TRUE. |
| 130 |
*ML: |
VHIT(7) = DELOSS + VHIT(7) |
|
* VHIT(7) = DELOSS + VHIT(7) |
|
|
VHIT(11) = DELOSS + VHIT(11) |
|
|
*END ML. |
|
| 131 |
ENDIF |
ENDIF |
| 132 |
VHIT(8) = (TIME + VHIT(8))/2. |
VHIT(8) = (TIME + VHIT(8))/2. |
| 133 |
VHIT(9) = PATH - VHIT(9) |
VHIT(9) = PATH - VHIT(9) |
| 134 |
*ML: |
c$$$ VHIT(11) = (TRAPAR(7) + VHIT(11))/2. |
|
* VHIT(11) = (TRAPAR(7) + VHIT(11))/2. |
|
|
VHIT(7) = (TRAPAR(7) + VHIT(7))/2. |
|
|
*END ML. |
|
| 135 |
* |
* |
| 136 |
* Store the hit |
* Store the hit |
| 137 |
* |
* |
| 138 |
IF(IDHIT.EQ.-1) THEN |
IF(IDHIT.EQ.-1) THEN |
|
*ML: |
|
|
* IF (SAVE) CALL GSAHIT (IS,ID,IT,NUMVOL,VHIT,IDHIT) |
|
|
* IF (IDHIT.EQ.0) THEN |
|
|
* WRITE(CHMAIL,10200) IS,ID,IT |
|
|
* CALL GMAIL(1,0) |
|
|
* ENDIF |
|
|
* ELSE |
|
|
* WRITE(CHMAIL,10300) IS,ID,IT |
|
|
* CALL GMAIL(1,0) |
|
|
* ENDIF |
|
| 139 |
IF (SAVE) THEN |
IF (SAVE) THEN |
| 140 |
IF(IS.EQ.ISS4) THEN |
IF(IS.EQ.ISS4) THEN |
| 141 |
CALL GSCHIT(IS,ID,1,NUMVOL,VHIT,1,IDHIT) |
CALL GSCHIT(IS,ID,1,NUMVOL,VHIT,1,IDHIT) |
| 151 |
WRITE(CHMAIL,10300) IS,ID,IT |
WRITE(CHMAIL,10300) IS,ID,IT |
| 152 |
CALL GMAIL(1,0) |
CALL GMAIL(1,0) |
| 153 |
ENDIF |
ENDIF |
|
*END ML. |
|
|
* IF (SAVE) CALL GSAHIT (IS,ID,IT,NUMVOL,VHIT,IDHIT) |
|
|
* IF (IDHIT.EQ.0) THEN |
|
|
* WRITE(CHMAIL,10200) IS,ID,IT |
|
|
* CALL GMAIL(1,0) |
|
|
* ENDIF |
|
|
* ELSE |
|
|
* WRITE(CHMAIL,10300) IS,ID,IT |
|
|
* CALL GMAIL(1,0) |
|
|
* ENDIF |
|
| 154 |
ENDIF |
ENDIF |
| 155 |
ENDIF |
ENDIF |
| 156 |
10000 FORMAT(' GPUTOF error: HIT have not be stored for SET=' |
10000 FORMAT(' GPUTOF error: HIT have not be stored for SET=' |