/[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.6 by cafagna, Thu May 11 00:39:39 2006 UTC revision 3.7 by cafagna, Thu May 11 23:52:55 2006 UTC
# Line 1  Line 1 
1  *  *
2  * $Id: gpdcal.F,v 3.5 2006/05/11 00:15:34 cafagna Exp $  * $Id: gpdcal.F,v 3.6 2006/05/11 00:39:39 cafagna Exp $
3  *  *
4  * $Log: gpdcal.F,v $  * $Log: gpdcal.F,v $
5    * Revision 3.6  2006/05/11 00:39:39  cafagna
6    * A bug in hits retrival for calorimeter
7    *
8  * Revision 3.5  2006/05/11 00:15:34  cafagna  * Revision 3.5  2006/05/11 00:15:34  cafagna
9  * A bug in hits retrival for calorimeter  * A bug in hits retrival for calorimeter
10  *  *
# Line 11  Line 14 
14  * Revision 3.3  2005/06/21 02:42:03  cafagna  * Revision 3.3  2005/06/21 02:42:03  cafagna
15  * Major modification to the geometry and to the random number chain  * Major modification to the geometry and to the random number chain
16  *  *
17  * Revision 3.2  2003/12/17 11:32:48  pamela  *     Revision 3.2 2003/12/17 11:32:48 pamela CALO SIMULATION COMPLETED:
18  * CALO SIMULATION COMPLETED: geometry and special tracking parameters updated and simulation checked by a comparison with the Trieste's standalone Monte Carlo simulation  *     geometry and special tracking parameters updated and simulation
19    *     checked by a comparison with the Trieste's standalone Monte Carlo
20    *     simulation
21  *  *
22  * Revision 3.1.1.1  2002/07/11 16:02:14  cafagna  * Revision 3.1.1.1  2002/07/11 16:02:14  cafagna
23  * First GPAMELA release on CVS  * First GPAMELA release on CVS
# Line 110  c$$$      ENDIF Line 115  c$$$      ENDIF
115  *  *
116  * If there was a problem in GPRHIT then exit  * If there was a problem in GPRHIT then exit
117  *  *
118          NTHCAL=0
119        IF(NHITS.LT.0) THEN        IF(NHITS.LT.0) THEN
120           WRITE(CHMAIL,10200) CSET,CNSED(IS,ID),NHITS,NHM,           WRITE(CHMAIL,10200) CSET,CNSED(IS,ID),NHITS,NHM,
121       +                        IDRUN,IDEVT       +                        IDRUN,IDEVT
122           CALL GMAIL(1,0)           CALL GMAIL(1,0)
123        ELSE           NHITS = NNHM
124          ENDIF
125  *  *
126  * >>>>>>>>>>>>>>>>>>>>>> DIGI action  * >>>>>>>>>>>>>>>>>>>>>> DIGI action
127  *  *
128  *  *
129  * Fill the CWN common  * Fill the CWN common
130  *  *
131           NTHCAL=0        DO I=1,NHITS
132           DO I=1,NHITS           NTHCAL = NTHCAL + 1
133              NTHCAL = NTHCAL + 1           ICAPL(NTHCAL) = NUMBV((NVCAL-2),I)
134              ICAPL(NTHCAL) = NUMBV((NVCAL-2),I)           ICASI(NTHCAL) = NUMBV((NVCAL-1),I)
135              ICASI(NTHCAL) = NUMBV((NVCAL-1),I)           ICAST(NTHCAL) = NUMBV((NVCAL),I)
136              ICAST(NTHCAL) = NUMBV((NVCAL),I)           XINCAL(NTHCAL) = VHITS(1,I)
137              XINCAL(NTHCAL) = VHITS(1,I)           YINCAL(NTHCAL) = VHITS(2,I)
138              YINCAL(NTHCAL) = VHITS(2,I)           ZINCAL(NTHCAL) = VHITS(3,I)
139              ZINCAL(NTHCAL) = VHITS(3,I)           ERELCAL(NTHCAL) = VHITS(4,I)
140              ERELCAL(NTHCAL) = VHITS(4,I)        ENDDO
 c         IPARCAL(NTHCAL) = VHITS(5,I)  
 c     P0CAL(NTHCAL) = VHITS(6,I)  
          ENDDO  
       ENDIF  
141  C  C
142  C The energy released per strip  C The energy released per strip
143  C  C
# Line 175  C Line 178  C
178  *  *
179  * If there was a problem in GPRHIT then exit  * If there was a problem in GPRHIT then exit
180  *  *
181          NTHCALI=0
182        IF(NHITS.LT.0) THEN        IF(NHITS.LT.0) THEN
183           WRITE(CHMAIL,10200) CSET,CNSED(IS,ID),NHITS,NHM,           WRITE(CHMAIL,10200) CSET,CNSED(IS,ID),NHITS,NHM,
184       +                        IDRUN,IDEVT       +                        IDRUN,IDEVT
# Line 188  C Line 192  C
192  * Fill the CWN common  * Fill the CWN common
193  *  *
194  c      c    
       NTHCALI=0  
195        DO I=1,NHITS        DO I=1,NHITS
196           NTHCALI=NTHCALI+1           NTHCALI=NTHCALI+1
197           ICASTRIP(NTHCALI)=NUMBV1(NVCALI,I)           ICASTRIP(NTHCALI)=NUMBV1(NVCALI,I)
Legend:
Removed from v.3.6  
changed lines
  Added in v.3.7
  ViewVC Help
Powered by ViewVC 1.1.23