1 |
* |
* |
2 |
* $Id: gpucal.F,v 3.2 2003/12/17 11:32:48 pamela Exp $ |
* $Id: gpucal.F,v 3.4 2005/12/14 03:13:53 cafagna Exp $ |
3 |
* |
* |
4 |
* $Log: gpucal.F,v $ |
* $Log: gpucal.F,v $ |
5 |
|
* Revision 3.4 2005/12/14 03:13:53 cafagna |
6 |
|
* Neutron detector added. Geometry and GPCALOR package |
7 |
|
* |
8 |
|
* Revision 3.3 2005/06/21 02:42:04 cafagna |
9 |
|
* Major modification to the geometry and to the random number chain |
10 |
|
* |
11 |
* Revision 3.2 2003/12/17 11:32:48 pamela |
* Revision 3.2 2003/12/17 11:32:48 pamela |
12 |
* CALO SIMULATION COMPLETED: geometry and special tracking parameters updated and simulation checked by a comparison with the Trieste's standalone Monte Carlo simulation |
* CALO SIMULATION COMPLETED: geometry and special tracking parameters updated and simulation checked by a comparison with the Trieste's standalone Monte Carlo simulation |
13 |
* |
* |
57 |
#include "gcflag.inc" |
#include "gcflag.inc" |
58 |
#include "gpencal.inc" |
#include "gpencal.inc" |
59 |
*end ml |
*end ml |
60 |
*EM: |
*EM: |
61 |
PARAMETER (CALIB=0.0001085200) |
PARAMETER (CALIB=0.0001085200) |
62 |
*END EM. |
*END EM. |
63 |
INTEGER IACT,IS,ID,IPA,IT,NUMVOL(20),IMEC,IDHIT |
INTEGER IACT,IS,ID,IPA,IT,NUMVOL(20),IMEC,IDHIT |
72 |
* |
* |
73 |
* It's a new event |
* It's a new event |
74 |
* |
* |
75 |
IF(IEVOLD.NE.IEVENT) THEN |
C # IF(IEVOLD.NE.IEVENT) THEN |
76 |
c print*,'IS,ID,NUMVOL=',IS,ID,NUMVOL |
C # c print*,'IS,ID,NUMVOL=',IS,ID,NUMVOL |
77 |
c print*,'GPUCAL:last event=',ievold |
C # c print*,'GPUCAL:last event=',ievold |
78 |
|
C # IEVOLD=IEVENT |
79 |
|
C # CALL VZERO(ENEMAT,44*9*32) |
80 |
|
C # ENDIF |
81 |
|
|
82 |
|
C # Please note that ENEMAT is zeroed just on the very first event here. |
83 |
|
C # Then is cleared into GUTREV on an event by event basis. |
84 |
|
IF(IEVOLD.EQ.-1) THEN |
85 |
IEVOLD=IEVENT |
IEVOLD=IEVENT |
86 |
CALL VZERO(ENEMAT,44*9*32) |
CALL VZERO(ENEMAT,44*9*32) |
87 |
ENDIF |
ENDIF |
88 |
|
|
89 |
*end ml. |
*end ml. |
90 |
* |
* |
91 |
* Track is inside a volume |
* Track is inside a volume |