/[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.1 by cafagna, Thu Jul 11 16:02:14 2002 UTC revision 3.2 by pamela, Wed Dec 17 11:32:48 2003 UTC
# Line 1  Line 1 
1  *  *
2  * $Id$  * $Id: gpdcal.F,v 3.1.1.1 2002/07/11 16:02:14 cafagna Exp $
3    *
4    * $Log: gpdcal.F,v $
5    * Revision 3.1.1.1  2002/07/11 16:02:14  cafagna
6    * First GPAMELA release on CVS
7  *  *
 * $Log$  
8  *  *
9  *CMZ :  2.03/00 31/10/2000  11.18.51  by  Francesco Cafagna  *CMZ :  2.03/00 31/10/2000  11.18.51  by  Francesco Cafagna
10  *CMZ :  2.02/00 10/10/2000  14.09.00  by  Francesco Cafagna  *CMZ :  2.02/00 10/10/2000  14.09.00  by  Francesco Cafagna
# Line 28  Line 31 
31  #include "gpsed.inc"  #include "gpsed.inc"
32  #include "gpkey.inc"  #include "gpkey.inc"
33  #include "gpgeo.inc"  #include "gpgeo.inc"
34    *ml:21/02/02
35    #include "gpencal.inc"
36    *end ml.
37  *  *
38        INTEGER IS,ID,NV,NVM,NNHM        INTEGER IS,ID,NV,NVM,NNHM
39        CHARACTER*4 CSET        CHARACTER*4 CSET
# Line 45  Line 51 
51        CALL VZERO(ERCASI,NHMCAL)        CALL VZERO(ERCASI,NHMCAL)
52        CALL VZERO(NUMD  ,NVOL)        CALL VZERO(NUMD  ,NVOL)
53        CALL VZERO(NUMD1 ,NVOL1)        CALL VZERO(NUMD1 ,NVOL1)
54    c ml:16/4/2003
55          CALL VZERO(ICASTRIP,NHMCAL)
56          CALL VZERO(ICAPLANE,NHMCAL)
57          CALL VZERO(ICAMOD,NHMCAL)
58          CALL VZERO(ENESTRIP,NHMCAL)
59    c end 16/4/2003
60  *  *
61        NV = NVOL        NV = NVOL
62        NVM= NVMAX        NVM= NVMAX
       NTHCAL = 0  
63  *  *
64  * The tungsten plane case  * The tungsten plane case
65  *  *
# Line 96  Line 107 
107  *  *
108  * Fill the CWN common  * Fill the CWN common
109  *  *
110          NTHCAL=0
111        DO I=1,NHITS        DO I=1,NHITS
112           NTHCAL = NTHCAL + 1           NTHCAL=NTHCAL+1
113           ICAPL(NTHCAL) = NUMBV((NVCAL-2),I)           ICASTRIP(NTHCAL)=VHITS(2,I)
114           ICASI(NTHCAL) = NUMBV((NVCAL-1),I)           ICAPLANE(NTHCAL)=VHITS(1,I)
115           ICAST(NTHCAL) = NUMBV((NVCAL),I)           ICAMOD(NTHCAL)=VHITS(3,I)
116           XINCAL(NTHCAL) = VHITS(1,I)           ENESTRIP(NTHCAL)=VHITS(4,I)
117           YINCAL(NTHCAL) = VHITS(2,I)        ENDDO
118           ZINCAL(NTHCAL) = VHITS(3,I)  cl      DO I=1,NHITS
119           ERELCAL(NTHCAL) = VHITS(4,I)  cl         NTHCAL = NTHCAL + 1
120    cl         ICAPL(NTHCAL) = NUMBV((NVCAL-2),I)
121    cl         ICASI(NTHCAL) = NUMBV((NVCAL-1),I)
122    cl         ICAST(NTHCAL) = NUMBV((NVCAL),I)
123    cl         XINCAL(NTHCAL) = VHITS(1,I)
124    cl         YINCAL(NTHCAL) = VHITS(2,I)
125    cl         ZINCAL(NTHCAL) = VHITS(3,I)
126    cl         ERELCAL(NTHCAL) = VHITS(4,I)
127  c         IPARCAL(NTHCAL) = VHITS(5,I)  c         IPARCAL(NTHCAL) = VHITS(5,I)
128  c         P0CAL(NTHCAL) = VHITS(6,I)  c         P0CAL(NTHCAL) = VHITS(6,I)
129        ENDDO  cl      ENDDO
130  10000  FORMAT(' GPDCAL: VOLUME level too big: SET=',A4,', DET=',A4,',  c    
131       +        NVOL=',I3,', NVMAX=',I3,', RUN=',I5,', EVT=',I8)  10000 FORMAT(' GPDCAL: VOLUME level too big: SET=',A4,', DET=',A4,',
132  10100  FORMAT(' GPDCAL: MAX number of hits retrived: SET=',A4,       +     NVOL=',I3,', NVMAX=',I3,', RUN=',I5,', EVT=',I8)
133       +        ', DET=',A4,',NHREAD=',I5,', NHMAX=',I5,  10100 FORMAT(' GPDCAL: MAX number of hits retrived: SET=',A4,
134       +        ', RUN=',I5,', EVT=',I8)       +     ', DET=',A4,',NHREAD=',I5,', NHMAX=',I5,
135  10200  FORMAT(' GPDCAL: Error in GPRHIT: SET=',A4,       +     ', RUN=',I5,', EVT=',I8)
136       +        ', DET=',A4,',NHITS=',I5,', NHMAX=',I5,  10200 FORMAT(' GPDCAL: Error in GPRHIT: SET=',A4,
137       +        ', RUN=',I5,', EVT=',I8)       +     ', DET=',A4,',NHITS=',I5,', NHMAX=',I5,
138     10 RETURN       +     ', RUN=',I5,', EVT=',I8)
139     10   RETURN
140        END        END
Legend:
Removed from v.3.1  
changed lines
  Added in v.3.2
  ViewVC Help
Powered by ViewVC 1.1.23