/[PAMELA software]/gpamela/gpcas/gpdcas.F
ViewVC logotype

Diff of /gpamela/gpcas/gpdcas.F

Parent Directory Parent Directory | Revision Log Revision Log | View Patch Patch

revision 3.1.1.1 by cafagna, Thu Jul 11 16:02:14 2002 UTC revision 3.3 by lundq, Thu Aug 25 11:13:18 2005 UTC
# Line 1  Line 1 
1  *  *
2  * $Id$  * $Id: gpdcas.F,v 3.2 2002/12/05 10:17:41 pamela Exp $
3  *  *
4  * $Log$  * $Log: gpdcas.F,v $
5    * Revision 3.2  2002/12/05 10:17:41  pamela
6    * Update CAS and CALO geometries and positions. Makefile updated as well
7  *  *
8    * Revision 3.1.1.1  2002/07/11 16:02:14  cafagna
9    * First GPAMELA release on CVS
10    *
11    *
12    *                19/09/2002  11.35.11  by  Jens Lund
13  *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
14  *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
15  *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
# Line 83  Line 90 
90  *  *
91        NTHCAS = NHITS        NTHCAS = NHITS
92        DO I=1,NHITS        DO I=1,NHITS
93           ICAS(I)  = NUMBV(NVCAS,I)  *JoLu 05-08-25
94    *         ICAS(I)  = NUMBV(NVCAS,I)
95             ICAS(I)  = NUMBV(1,I)
96  #if defined(GPAMELA_INTER)  #if defined(GPAMELA_INTER)
97           XINCAS(I)  = 2*VHITS(1,I)-VHITS(4,I)           XINCAS(I)  = 2*VHITS(1,I)-VHITS(4,I)
98           YINCAS(I)  = 2*VHITS(2,I)-VHITS(5,I)           YINCAS(I)  = 2*VHITS(2,I)-VHITS(5,I)
# Line 103  Line 112 
112           IPARCAS(I) = VHITS(10,I)           IPARCAS(I) = VHITS(10,I)
113           P0CAS(I) = VHITS(11,I)           P0CAS(I) = VHITS(11,I)
114        ENDDO        ENDDO
115    *JeL:
116  *  *
117  * Next detector  * Next detector
118  *  *
119        ID   = IDCASY  C      ID   = IDCASY
120        CALL VZERO(VHITS,NH*NHM)  C      CALL VZERO(VHITS,NH*NHM)
121        CALL VZERO(ITRA ,NHM)  C      CALL VZERO(ITRA ,NHM)
122        CALL VZERO(NUMBV,NVOL*NHM)  C      CALL VZERO(NUMBV,NVOL*NHM)
123        NHITS = -1  C      NHITS = -1
124  *  *
125  * Update the max number of hits to be retrived  * Update the max number of hits to be retrived
126  *  *
127        NNHM = NHM - NTHCAS  C      NNHM = NHM - NTHCAS
128        IF(NNHM.LE.0) THEN  C      IF(NNHM.LE.0) THEN
129           WRITE(CHMAIL,10100) CSET,CNSED(IS,ID),NTHCAS,NHM,  C         WRITE(CHMAIL,10100) CSET,CNSED(IS,ID),NTHCAS,NHM,
130       +                        IDRUN,IDEVT  C     +                        IDRUN,IDEVT
131           CALL GMAIL(1,0)  C         CALL GMAIL(1,0)
132           GO TO 10  C         GO TO 10
133        ENDIF  C      ENDIF
134        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,
135       +            IT,NUMD,ITRA,NUMBV,VHITS,NHITS)  C     +            IT,NUMD,ITRA,NUMBV,VHITS,NHITS)
136  *  *
137  * If there was a problem in GPRHIT then exit  * If there was a problem in GPRHIT then exit
138  *  *
139        IF(NHITS.LT.0) THEN  C      IF(NHITS.LT.0) THEN
140           WRITE(CHMAIL,10200) CSET,CNSED(IS,ID),NHITS,NHM,  C         WRITE(CHMAIL,10200) CSET,CNSED(IS,ID),NHITS,NHM,
141       +                        IDRUN,IDEVT  C     +                        IDRUN,IDEVT
142           CALL GMAIL(1,0)  C         CALL GMAIL(1,0)
143           GO TO 10  C         GO TO 10
144        ENDIF  C      ENDIF
145  *  *
146  * >>>>>>>>>>>>>>>>>>>>>> DIGI action  * >>>>>>>>>>>>>>>>>>>>>> DIGI action
147  *  *
148  *  *
149  * Fill the CWN common  * Fill the CWN common
150  *  *
151        DO I=1,NHITS  C      DO I=1,NHITS
152           NTHCAS = NTHCAS + 1  C         NTHCAS = NTHCAS + 1
153           ICAS(NTHCAS)  = 2 + NUMBV(NVCAS,I)  C         ICAS(NTHCAS)  = 2 + NUMBV(NVCAS,I)
154  #if defined(GPAMELA_INTER)  C#if defined(GPAMELA_INTER)
155           XINCAS(NTHCAS)  = 2*VHITS(1,I)-VHITS(4,I)  C         XINCAS(NTHCAS)  = 2*VHITS(1,I)-VHITS(4,I)
156           YINCAS(NTHCAS)  = 2*VHITS(2,I)-VHITS(5,I)  C         YINCAS(NTHCAS)  = 2*VHITS(2,I)-VHITS(5,I)
157           ZINCAS(NTHCAS)  = 2*VHITS(3,I)-VHITS(6,I)  C         ZINCAS(NTHCAS)  = 2*VHITS(3,I)-VHITS(6,I)
158  #endif  C#endif
159  #if !defined(GPAMELA_INTER)  C#if !defined(GPAMELA_INTER)
160           XINCAS(NTHCAS)  = VHITS(1,I)  C         XINCAS(NTHCAS)  = VHITS(1,I)
161           YINCAS(NTHCAS)  = VHITS(2,I)  C         YINCAS(NTHCAS)  = VHITS(2,I)
162           ZINCAS(NTHCAS)  = VHITS(3,I)  C         ZINCAS(NTHCAS)  = VHITS(3,I)
163  #endif  C#endif
164           XOUTCAS(NTHCAS) = VHITS(4,I)  C         XOUTCAS(NTHCAS) = VHITS(4,I)
165           YOUTCAS(NTHCAS) = VHITS(5,I)  C         YOUTCAS(NTHCAS) = VHITS(5,I)
166           ZOUTCAS(NTHCAS) = VHITS(6,I)  C         ZOUTCAS(NTHCAS) = VHITS(6,I)
167           ERELCAS(NTHCAS) = VHITS(7,I)  C         ERELCAS(NTHCAS) = VHITS(7,I)
168           TIMECAS(NTHCAS)=  VHITS(8,I)  C         TIMECAS(NTHCAS)=  VHITS(8,I)
169           PATHCAS(NTHCAS) = VHITS(9,I)  C         PATHCAS(NTHCAS) = VHITS(9,I)
170           IPARCAS(NTHCAS) = VHITS(10,I)  C         IPARCAS(NTHCAS) = VHITS(10,I)
171           P0CAS(NTHCAS) = VHITS(11,I)  C         P0CAS(NTHCAS) = VHITS(11,I)
172        ENDDO  C      ENDDO
173    *END: JeL.
174  *  *
175  10000  FORMAT(' GPDCAS: VOLUME level too big: SET=',A4,', DET=',A4,',  10000  FORMAT(' GPDCAS: VOLUME level too big: SET=',A4,', DET=',A4,',
176       +        NVOL=',I3,', NVMAX=',I3,', RUN=',I5,', EVT=',I8)       +        NVOL=',I3,', NVMAX=',I3,', RUN=',I5,', EVT=',I8)
Legend:
Removed from v.3.1.1.1  
changed lines
  Added in v.3.3
  ViewVC Help
Powered by ViewVC 1.1.23