/[PAMELA software]/gpamela/gpcat/gpdcat.F
ViewVC logotype

Diff of /gpamela/gpcat/gpdcat.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:42 2002 UTC
# Line 1  Line 1 
1  *  *
2  * $Id$  * $Id: gpdcat.F,v 3.1.1.1 2002/07/11 16:02:14 cafagna Exp $
3  *  *
4  * $Log$  * $Log: gpdcat.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.30.50  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 105  Line 109 
109           IPARCAT(I) = VHITS(10,I)           IPARCAT(I) = VHITS(10,I)
110           P0CAT(I) = VHITS(11,I)           P0CAT(I) = VHITS(11,I)
111        ENDDO        ENDDO
112    *JeL:
113  *  *
114  * Next detector  * Next detector
115  *  *
116        ID   = IDCATT  C      ID   = IDCATT
117        CALL VZERO(VHITS,NH*NHM)  C      CALL VZERO(VHITS,NH*NHM)
118        CALL VZERO(ITRA ,NHM)  C      CALL VZERO(ITRA ,NHM)
119        CALL VZERO(NUMBV,NVOL*NHM)  C      CALL VZERO(NUMBV,NVOL*NHM)
120        NHITS = -1  C      NHITS = -1
121  *  *
122  * Update the max number of hits to be retrived  * Update the max number of hits to be retrived
123  *  *
124        NNHM = NHM - NTHCAT  C      NNHM = NHM - NTHCAT
125        IF(NNHM.LE.0) THEN  C      IF(NNHM.LE.0) THEN
126           WRITE(CHMAIL,10100) CSET,CNSED(IS,ID),NTHCAT,NHM,  C         WRITE(CHMAIL,10100) CSET,CNSED(IS,ID),NTHCAT,NHM,
127       +                        IDRUN,IDEVT  C     +                        IDRUN,IDEVT
128           CALL GMAIL(1,0)  C         CALL GMAIL(1,0)
129           GO TO 10  C         GO TO 10
130        ENDIF  C      ENDIF
131        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,
132       +            IT,NUMD,ITRA,NUMBV,VHITS,NHITS)  C     +            IT,NUMD,ITRA,NUMBV,VHITS,NHITS)
133  *  *
134  * If there was a problem in GPRHIT then exit  * If there was a problem in GPRHIT then exit
135  *  *
136        IF(NHITS.LT.0) THEN  C      IF(NHITS.LT.0) THEN
137           WRITE(CHMAIL,10200) CSET,CNSED(IS,ID),NHITS,NHM,  C         WRITE(CHMAIL,10200) CSET,CNSED(IS,ID),NHITS,NHM,
138       +                        IDRUN,IDEVT  C     +                        IDRUN,IDEVT
139           CALL GMAIL(1,0)  C         CALL GMAIL(1,0)
140           GO TO 10  C         GO TO 10
141        ENDIF  C      ENDIF
142  *  *
143  * >>>>>>>>>>>>>>>>>>>>>> DIGI action  * >>>>>>>>>>>>>>>>>>>>>> DIGI action
144  *  *
145  *  *
146        DO I=1,NHITS  C      DO I=1,NHITS
147           NTHCAT = NTHCAT + 1  C         NTHCAT = NTHCAT + 1
148           ICAT(NTHCAT)  = 2 + NUMBV(NVCAT,I)  C         ICAT(NTHCAT)  = 2 + NUMBV(NVCAT,I)
149  #if defined(GPAMELA_INTER)  C#if defined(GPAMELA_INTER)
150           XINCAT(NTHCAT)  = 2*VHITS(1,I)-VHITS(4,I)  C         XINCAT(NTHCAT)  = 2*VHITS(1,I)-VHITS(4,I)
151           YINCAT(NTHCAT)  = 2*VHITS(2,I)-VHITS(5,I)  C         YINCAT(NTHCAT)  = 2*VHITS(2,I)-VHITS(5,I)
152           ZINCAT(NTHCAT)  = 2*VHITS(3,I)-VHITS(6,I)  C         ZINCAT(NTHCAT)  = 2*VHITS(3,I)-VHITS(6,I)
153  #endif  C#endif
154  #if !defined(GPAMELA_INTER)  C#if !defined(GPAMELA_INTER)
155           XINCAT(NTHCAT)  = VHITS(1,I)  C         XINCAT(NTHCAT)  = VHITS(1,I)
156           YINCAT(NTHCAT)  = VHITS(2,I)  C         YINCAT(NTHCAT)  = VHITS(2,I)
157           ZINCAT(NTHCAT)  = VHITS(3,I)  C         ZINCAT(NTHCAT)  = VHITS(3,I)
158  #endif  C#endif
159           XOUTCAT(NTHCAT) = VHITS(4,I)  C         XOUTCAT(NTHCAT) = VHITS(4,I)
160           YOUTCAT(NTHCAT) = VHITS(5,I)  C         YOUTCAT(NTHCAT) = VHITS(5,I)
161           ZOUTCAT(NTHCAT) = VHITS(6,I)  C         ZOUTCAT(NTHCAT) = VHITS(6,I)
162           ERELCAT(NTHCAT) = VHITS(7,I)  C         ERELCAT(NTHCAT) = VHITS(7,I)
163           TIMECAT(NTHCAT)=  VHITS(8,I)  C         TIMECAT(NTHCAT)=  VHITS(8,I)
164           PATHCAT(NTHCAT) = VHITS(9,I)  C         PATHCAT(NTHCAT) = VHITS(9,I)
165           IPARCAT(NTHCAT) = VHITS(10,I)  C         IPARCAT(NTHCAT) = VHITS(10,I)
166           P0CAT(NTHCAT) = VHITS(11,I)  C         P0CAT(NTHCAT) = VHITS(11,I)
167        ENDDO  C      ENDDO
168    *END: JeL.
169  10000  FORMAT(' GPDCAT: VOLUME level too big: SET=',A4,', DET=',A4,',  10000  FORMAT(' GPDCAT: VOLUME level too big: SET=',A4,', DET=',A4,',
170       +        NVOL=',I3,', NVMAX=',I3,', RUN=',I5,', EVT=',I8)       +        NVOL=',I3,', NVMAX=',I3,', RUN=',I5,', EVT=',I8)
171  10100  FORMAT(' GPDCAT: MAX number of hits retrived: SET=',A4,  10100  FORMAT(' GPDCAT: MAX number of hits retrived: SET=',A4,
Legend:
Removed from v.3.1.1.1  
changed lines
  Added in v.3.2
  ViewVC Help
Powered by ViewVC 1.1.23