/[PAMELA software]/gpamela/gpcal/gpdcal.F
ViewVC logotype

Diff of /gpamela/gpcal/gpdcal.F

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

revision 3.2 by pamela, Wed Dec 17 11:32:48 2003 UTC revision 3.3 by cafagna, Tue Jun 21 02:42:03 2005 UTC
# Line 1  Line 1 
1  *  *
2  * $Id: gpdcal.F,v 3.1.1.1 2002/07/11 16:02:14 cafagna Exp $  * $Id: gpdcal.F,v 3.2 2003/12/17 11:32:48 pamela Exp $
3  *  *
4  * $Log: gpdcal.F,v $  * $Log: gpdcal.F,v $
5    * Revision 3.2  2003/12/17 11:32:48  pamela
6    * CALO SIMULATION COMPLETED: geometry and special tracking parameters updated and simulation checked by a comparison with the Trieste's standalone Monte Carlo simulation
7    *
8  * Revision 3.1.1.1  2002/07/11 16:02:14  cafagna  * Revision 3.1.1.1  2002/07/11 16:02:14  cafagna
9  * First GPAMELA release on CVS  * First GPAMELA release on CVS
10  *  *
# Line 64  c end 16/4/2003 Line 67  c end 16/4/2003
67  * The tungsten plane case  * The tungsten plane case
68  *  *
69        IS   = ISCAL        IS   = ISCAL
70        ID   = IDCASI        ID   = IDCAST
71        IT = -1        IT = -1
72        CSET = 'CAL '        CSET = 'CAL '
73  *  *
# Line 109  c end 16/4/2003 Line 112  c end 16/4/2003
112  *  *
113        NTHCAL=0        NTHCAL=0
114        DO I=1,NHITS        DO I=1,NHITS
115           NTHCAL=NTHCAL+1           NTHCAL = NTHCAL + 1
116           ICASTRIP(NTHCAL)=VHITS(2,I)           ICAPL(NTHCAL) = NUMBV((NVCAL-2),I)
117           ICAPLANE(NTHCAL)=VHITS(1,I)           ICASI(NTHCAL) = NUMBV((NVCAL-1),I)
118           ICAMOD(NTHCAL)=VHITS(3,I)           ICAST(NTHCAL) = NUMBV((NVCAL),I)
119           ENESTRIP(NTHCAL)=VHITS(4,I)           XINCAL(NTHCAL) = VHITS(1,I)
120        ENDDO           YINCAL(NTHCAL) = VHITS(2,I)
121  cl      DO I=1,NHITS           ZINCAL(NTHCAL) = VHITS(3,I)
122  cl         NTHCAL = NTHCAL + 1           ERELCAL(NTHCAL) = VHITS(4,I)
 cl         ICAPL(NTHCAL) = NUMBV((NVCAL-2),I)  
 cl         ICASI(NTHCAL) = NUMBV((NVCAL-1),I)  
 cl         ICAST(NTHCAL) = NUMBV((NVCAL),I)  
 cl         XINCAL(NTHCAL) = VHITS(1,I)  
 cl         YINCAL(NTHCAL) = VHITS(2,I)  
 cl         ZINCAL(NTHCAL) = VHITS(3,I)  
 cl         ERELCAL(NTHCAL) = VHITS(4,I)  
123  c         IPARCAL(NTHCAL) = VHITS(5,I)  c         IPARCAL(NTHCAL) = VHITS(5,I)
124  c         P0CAL(NTHCAL) = VHITS(6,I)  c         P0CAL(NTHCAL) = VHITS(6,I)
125  cl      ENDDO        ENDDO
126    C
127    C The energy released per strip
128    C
129          NV = NVOL1
130          NVM= NVMAX
131    *
132    * The tungsten plane case
133    *
134          IS   = ISCAL
135          ID   = IDCASI
136          IT = -1
137          CSET = 'CAL '
138    *
139    * Check on the volume level
140    *
141          IF (NV.GT.NVM) THEN
142             WRITE(CHMAIL,10000) CSET,CNSED(IS,ID),NVOL,NVMAX,
143         +                        IDRUN,IDEVT
144             CALL GMAIL(1,0)
145             GO TO 10
146          ENDIF
147          CALL VZERO(VHITS,NH*NHM)
148          CALL VZERO(ITRA ,NHM)
149          NHITS = -1
150    *
151    * Update the max number of hits to be retrived
152    *
153          NNHM = NHM
154          IF(NNHM.LE.0) THEN
155             WRITE(CHMAIL,10100) CSET,CNSED(IS,ID),NTHCALI,NHM,
156         +                        IDRUN,IDEVT
157             CALL GMAIL(1,0)
158             GO TO 10
159          ENDIF
160          CALL VZERO(NUMBV,NVOL*NHM)
161          CALL GPRHIT(CSET,CNSED(IS,ID),IS,ID,NVOL,NH,NNHM,1,
162         +           IT,NUMD,ITRA,NUMBV,VHITS,NHITS)
163    *
164    * If there was a problem in GPRHIT then exit
165    *
166          IF(NHITS.LT.0) THEN
167             WRITE(CHMAIL,10200) CSET,CNSED(IS,ID),NHITS,NHM,
168         +                        IDRUN,IDEVT
169             CALL GMAIL(1,0)
170             GO TO 10
171          ENDIF
172    *
173    * >>>>>>>>>>>>>>>>>>>>>> DIGI action
174    *
175    *
176    * Fill the CWN common
177    *
178  c      c    
179          DO I=1,NHITS
180             NTHCALI=NTHCALI+1
181             ICASTRIP(NTHCALI)=NUMBV(NVCALI,I)
182             ICAPLANE(NTHCALI)=NUMBV((NVCALI-1),I)
183             ICAMOD(NTHCALI)=VHITS(1,I)
184             ENESTRIP(NTHCALI)=VHITS(2,I)
185          ENDDO
186    c
187  10000 FORMAT(' GPDCAL: VOLUME level too big: SET=',A4,', DET=',A4,',  10000 FORMAT(' GPDCAL: VOLUME level too big: SET=',A4,', DET=',A4,',
188       +     NVOL=',I3,', NVMAX=',I3,', RUN=',I5,', EVT=',I8)       +     NVOL=',I3,', NVMAX=',I3,', RUN=',I5,', EVT=',I8)
189  10100 FORMAT(' GPDCAL: MAX number of hits retrived: SET=',A4,  10100 FORMAT(' GPDCAL: MAX number of hits retrived: SET=',A4,
Legend:
Removed from v.3.2  
changed lines
  Added in v.3.3
  ViewVC Help
Powered by ViewVC 1.1.23