/[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.5 by cafagna, Thu May 11 00:15:34 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.4 2005/07/25 11:53:20 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
9    * A bug in hits retrival for calorimeter
10    *
11  * Revision 3.4  2005/07/25 11:53:20  cafagna  * Revision 3.4  2005/07/25 11:53:20  cafagna
12  * Several updates. See history for details  * Several updates. See history for details
13  *  *
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 107  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           GO TO 10           NHITS = NNHM
124        ENDIF        ENDIF
125  *  *
126  * >>>>>>>>>>>>>>>>>>>>>> DIGI action  * >>>>>>>>>>>>>>>>>>>>>> DIGI action
# Line 119  c$$$      ENDIF Line 128  c$$$      ENDIF
128  *  *
129  * Fill the CWN common  * Fill the CWN common
130  *  *
       NTHCAL=0  
131        DO I=1,NHITS        DO I=1,NHITS
132           NTHCAL = NTHCAL + 1           NTHCAL = NTHCAL + 1
133           ICAPL(NTHCAL) = NUMBV((NVCAL-2),I)           ICAPL(NTHCAL) = NUMBV((NVCAL-2),I)
# Line 129  c$$$      ENDIF Line 137  c$$$      ENDIF
137           YINCAL(NTHCAL) = VHITS(2,I)           YINCAL(NTHCAL) = VHITS(2,I)
138           ZINCAL(NTHCAL) = VHITS(3,I)           ZINCAL(NTHCAL) = VHITS(3,I)
139           ERELCAL(NTHCAL) = VHITS(4,I)           ERELCAL(NTHCAL) = VHITS(4,I)
 c         IPARCAL(NTHCAL) = VHITS(5,I)  
 c         P0CAL(NTHCAL) = VHITS(6,I)  
140        ENDDO        ENDDO
141  C  C
142  C The energy released per strip  C The energy released per strip
# Line 144  C Line 150  C
150        ID   = IDCASI        ID   = IDCASI
151        IT = -1        IT = -1
152        CSET = 'CAL '        CSET = 'CAL '
153  *  *    
154  * Check on the volume level  * Check on the volume level
155  *  *
156        IF (NV.GT.NVM) THEN        IF (NV.GT.NVM) THEN
# Line 172  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 185  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.5  
changed lines
  Added in v.3.7
  ViewVC Help
Powered by ViewVC 1.1.23