1 |
* |
* |
2 |
* $Id$ |
* $Id: gpdcas.F,v 3.1.1.1 2002/07/11 16:02:14 cafagna Exp $ |
3 |
* |
* |
4 |
* $Log$ |
* $Log: gpdcas.F,v $ |
5 |
|
* Revision 3.1.1.1 2002/07/11 16:02:14 cafagna |
6 |
|
* First GPAMELA release on CVS |
7 |
* |
* |
8 |
|
* |
9 |
|
* 19/09/2002 11.35.11 by Jens Lund |
10 |
*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 |
11 |
*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 |
12 |
*CMZ : 1.02/00 22/10/99 16.17.58 by Francesco Cafagna |
*CMZ : 1.02/00 22/10/99 16.17.58 by Francesco Cafagna |
107 |
IPARCAS(I) = VHITS(10,I) |
IPARCAS(I) = VHITS(10,I) |
108 |
P0CAS(I) = VHITS(11,I) |
P0CAS(I) = VHITS(11,I) |
109 |
ENDDO |
ENDDO |
110 |
|
*JeL: |
111 |
* |
* |
112 |
* Next detector |
* Next detector |
113 |
* |
* |
114 |
ID = IDCASY |
C ID = IDCASY |
115 |
CALL VZERO(VHITS,NH*NHM) |
C CALL VZERO(VHITS,NH*NHM) |
116 |
CALL VZERO(ITRA ,NHM) |
C CALL VZERO(ITRA ,NHM) |
117 |
CALL VZERO(NUMBV,NVOL*NHM) |
C CALL VZERO(NUMBV,NVOL*NHM) |
118 |
NHITS = -1 |
C NHITS = -1 |
119 |
* |
* |
120 |
* Update the max number of hits to be retrived |
* Update the max number of hits to be retrived |
121 |
* |
* |
122 |
NNHM = NHM - NTHCAS |
C NNHM = NHM - NTHCAS |
123 |
IF(NNHM.LE.0) THEN |
C IF(NNHM.LE.0) THEN |
124 |
WRITE(CHMAIL,10100) CSET,CNSED(IS,ID),NTHCAS,NHM, |
C WRITE(CHMAIL,10100) CSET,CNSED(IS,ID),NTHCAS,NHM, |
125 |
+ IDRUN,IDEVT |
C + IDRUN,IDEVT |
126 |
CALL GMAIL(1,0) |
C CALL GMAIL(1,0) |
127 |
GO TO 10 |
C GO TO 10 |
128 |
ENDIF |
C ENDIF |
129 |
CALL GPRHIT(CSET,CNSED(IS,ID),IS,ID,NVOL,NH,NNHM,1, |
C CALL GPRHIT(CSET,CNSED(IS,ID),IS,ID,NVOL,NH,NNHM,1, |
130 |
+ IT,NUMD,ITRA,NUMBV,VHITS,NHITS) |
C + IT,NUMD,ITRA,NUMBV,VHITS,NHITS) |
131 |
* |
* |
132 |
* If there was a problem in GPRHIT then exit |
* If there was a problem in GPRHIT then exit |
133 |
* |
* |
134 |
IF(NHITS.LT.0) THEN |
C IF(NHITS.LT.0) THEN |
135 |
WRITE(CHMAIL,10200) CSET,CNSED(IS,ID),NHITS,NHM, |
C WRITE(CHMAIL,10200) CSET,CNSED(IS,ID),NHITS,NHM, |
136 |
+ IDRUN,IDEVT |
C + IDRUN,IDEVT |
137 |
CALL GMAIL(1,0) |
C CALL GMAIL(1,0) |
138 |
GO TO 10 |
C GO TO 10 |
139 |
ENDIF |
C ENDIF |
140 |
* |
* |
141 |
* >>>>>>>>>>>>>>>>>>>>>> DIGI action |
* >>>>>>>>>>>>>>>>>>>>>> DIGI action |
142 |
* |
* |
143 |
* |
* |
144 |
* Fill the CWN common |
* Fill the CWN common |
145 |
* |
* |
146 |
DO I=1,NHITS |
C DO I=1,NHITS |
147 |
NTHCAS = NTHCAS + 1 |
C NTHCAS = NTHCAS + 1 |
148 |
ICAS(NTHCAS) = 2 + NUMBV(NVCAS,I) |
C ICAS(NTHCAS) = 2 + NUMBV(NVCAS,I) |
149 |
#if defined(GPAMELA_INTER) |
C#if defined(GPAMELA_INTER) |
150 |
XINCAS(NTHCAS) = 2*VHITS(1,I)-VHITS(4,I) |
C XINCAS(NTHCAS) = 2*VHITS(1,I)-VHITS(4,I) |
151 |
YINCAS(NTHCAS) = 2*VHITS(2,I)-VHITS(5,I) |
C YINCAS(NTHCAS) = 2*VHITS(2,I)-VHITS(5,I) |
152 |
ZINCAS(NTHCAS) = 2*VHITS(3,I)-VHITS(6,I) |
C ZINCAS(NTHCAS) = 2*VHITS(3,I)-VHITS(6,I) |
153 |
#endif |
C#endif |
154 |
#if !defined(GPAMELA_INTER) |
C#if !defined(GPAMELA_INTER) |
155 |
XINCAS(NTHCAS) = VHITS(1,I) |
C XINCAS(NTHCAS) = VHITS(1,I) |
156 |
YINCAS(NTHCAS) = VHITS(2,I) |
C YINCAS(NTHCAS) = VHITS(2,I) |
157 |
ZINCAS(NTHCAS) = VHITS(3,I) |
C ZINCAS(NTHCAS) = VHITS(3,I) |
158 |
#endif |
C#endif |
159 |
XOUTCAS(NTHCAS) = VHITS(4,I) |
C XOUTCAS(NTHCAS) = VHITS(4,I) |
160 |
YOUTCAS(NTHCAS) = VHITS(5,I) |
C YOUTCAS(NTHCAS) = VHITS(5,I) |
161 |
ZOUTCAS(NTHCAS) = VHITS(6,I) |
C ZOUTCAS(NTHCAS) = VHITS(6,I) |
162 |
ERELCAS(NTHCAS) = VHITS(7,I) |
C ERELCAS(NTHCAS) = VHITS(7,I) |
163 |
TIMECAS(NTHCAS)= VHITS(8,I) |
C TIMECAS(NTHCAS)= VHITS(8,I) |
164 |
PATHCAS(NTHCAS) = VHITS(9,I) |
C PATHCAS(NTHCAS) = VHITS(9,I) |
165 |
IPARCAS(NTHCAS) = VHITS(10,I) |
C IPARCAS(NTHCAS) = VHITS(10,I) |
166 |
P0CAS(NTHCAS) = VHITS(11,I) |
C P0CAS(NTHCAS) = VHITS(11,I) |
167 |
ENDDO |
C ENDDO |
168 |
|
*END: JeL. |
169 |
* |
* |
170 |
10000 FORMAT(' GPDCAS: VOLUME level too big: SET=',A4,', DET=',A4,', |
10000 FORMAT(' GPDCAS: VOLUME level too big: SET=',A4,', DET=',A4,', |
171 |
+ NVOL=',I3,', NVMAX=',I3,', RUN=',I5,', EVT=',I8) |
+ NVOL=',I3,', NVMAX=',I3,', RUN=',I5,', EVT=',I8) |