/[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.2 by pamela, Thu Dec 5 10:17:41 2002 UTC
# Line 1  Line 1 
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
# Line 103  Line 107 
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)

Legend:
Removed from v.3.1.1.1  
changed lines
  Added in v.3.2

  ViewVC Help
Powered by ViewVC 1.1.23